Q-SMAKAS Tooltip Failure Mode Search Project Results

Table of Contents

1. Experimental process
2. Results format explanation
3. Results
    a. Loaded Tooltips
         DC10c C100GeATD C100GeATS C100GeATSr5 C100GeATSr6 C100GeCTS
         C110GeS C111Ge3 AdamGe22 AdamGe33 Diad3Ge22 DCB55AGe
         DCB55BGe DCB55CGe DCB57Ge DCB5Ge DCB65Ge DCB6Ge
         DCB75Ge DCBIceane7Ge MCB55Ge MCB57Ge MCB5Ge MCB75Ge
         TwistaneGe
    b. Unloaded Tooltips
         DC10c C100GeATD C100GeATS C100GeATSr5 C100GeATSr6 C100GeCTS
         C110GeS C111Ge3 AdamGe22 AdamGe33 Diad3Ge22 DCB55AGe
         DCB55BGe DCB55CGe DCB57Ge DCB5Ge DCB65Ge DCB6Ge
         DCB75Ge DCBIceane7Ge MCB55Ge MCB57Ge MCB5Ge MCB75Ge
         TwistaneGe
    c. Q-SMAKAS Method Evaluation Molecules
         Biotin
         Heptane
         Octane

Experimental process

We performed an initial geometry optimization (aka energy minimization) on each tooltip using Sandia's Massively Parallel Quantum Chemistry program (MPQC) with the RHF/3-21G level of theory. Some optimizations were performed in-house, but most were computed with the NHAH network.

Some of the larger tooltips took too long (several days) to converge with RHF/3-21G, so RHF/STO-3G was used instead. These lower levels of theory are both sufficiently accurate to detect failure modes (see Note on High Temperature/Low Level Logic.) The tooltips processed with RHF/STO-3G were: MCB5Ge, MCB57Ge, DCB75Ge, C100GeATD, C100GeATS, C100GeCTS, and MCB75Ge.

Using the optimized tooltip geometries, we then performed GROMACS molecular dynamics (MD) simulations on each of the tooltips at a high temperature for 50,000 iterations, taking snapshots of the tooltip every 40 or 50 iterations. The thermal energy of the high temperature simulations perturbed each tooltip wildly and the snapshots generated (kinetically accessible) tooltip deformations.

Initially, simulation at 1500 Kelvin was performed with each tooltip as a baseline level of tooltip deformation. The root mean square deviation (RMSD) of the tooltips at this temperature was around 0.02. An RHF/3-21G|STO-3G geometry optimization was then performed for each deformation on the NHAH network. No failure modes were discovered with this level of deformation.

Next, we performed simulations at a temperature yeilding an RMSD of at least 0.04 in 500 Kelvin steps. The range of temperatures used were: 2000 K, 2500 K, 3000 K, 3500 K, 4000 K, and 4500 K. An RHF/3-21G|STO-3G geometry optimization was then performed for each deformation on the NHAH network. Several potential failure modes were discovered at this level of deformation; they are under analysis.

Results format explanation

The first column shows each of the tooltips in their normal, non-failure modes. There is also a link to the tooltip's PDB (Protein Data Bank) file which describes its geometry. (There are several graphical PDB viewers. The Swiss-PDB Viewer runs on Mac, Windows, and Linux.) An "L" in the PDB filename indicates that the tooltip is loaded, a "U" in the PDB filename indicates that the tooltip is un-loaded. The tooltip's energy (see Note on Energy Units) is also shown in this column.

The second column shows potential failure modes for the tooltip if any have been discovered so far. There is also a link to the mode's PDB file. The number and letter after the L/U indicates the deformed structure index and input dataset the mode came from. The mode's energy is also shown. If no potential failure modes have been discovered yet, the second column contains links to the input and output files (described below.)

The third column shows the number of initial deforms that minimized to the defect structure, a summary of the MD simulation performed to generate the initial deforms, and notes/observations about the defect. The MD simulation run summary line, eg, "MD run: 3500 K (0.04 RMSD)", shows the simulation environment temperature needed to produce tooltip deformations with the indicated RMSD.

The final row for each tooltip contains links to the input and results data files. The letter after the L/U indicates which simulation environment temperature was used to generate the tooltip deformations. The following table maps the letter to the temperature:

z1500 K    y2000 K    x2500 K    w3000 K    v3500 K    u4000 K    t4500 K

Notes

1. High Temperature/Low Level Logic
If these tooltips were to be treated as organic molecules, the identification of defect structures would correspond to the identification of reaction pathways between starting and final structures. Based on these structures, one would identify an organic "mechanism" that explains the rearrangement that must occur to a structure to change from starting to final forms. In order to do this correctly, a series of geometry optimizations would need to be performed on the reactant (starting structure), product (final structure) and one or more transition structure(s) along the reaction coordinate defined by the mechanism. In contrast, the high temperatures and (relatively) low levels of theory employed in this defect structure analysis are employed in order to "brute-force" the identification of defect structures, simply providing final geometries that may correspond to realizable minimum energy forms that can be used as the base structures for developing plausible mechanisms.

The low levels of theory (RHF/3-21G/STO-3G) are appropriate for these studies because this study is looking for strongly bound (covalent) tooltip structures. The strengths of covalent bonds are such that even low theory levels predict their connectivity correctly. The relative energies of these structures, however, will differ based on the level of theory. The generation of reasonable minimum energy candidate structures for the defect calculations requires more computational resources than the higher-level geometry optimizations of the identified defect structures because, for instance, 5000 defect starting structures taken from the GROMACS calculations may optimize to only 15 unique tooltip geometries, a consequence of not actively directing the generation of the defect structures.

2. Energy Units
The standard unit used in quantum chemistry for reporting the TOTAL ENERGY of molecules is the Hartree. The standard view in quantum chemistry is that a system of completely dissociated ("infinite separation") nuclei and electrons corresponds to zero energy (0 Hartree). This means the separation of a molecule into its component atoms AND the separation of each atom into its individual electrons and nucleus (the individual protons and neutrons in the nuclei are NOT included in this separation). As electrons and nuclei are brought closer to one another to form the molecule, the energy, as defined, decreases (more negative).

The only appropriate way to compare energies in quantum chemical calculations is to compare the SAME numbers of atoms at the SAME level of theory. Therefore, each tooltip and its defect structures are considered as distinct groups with no analysis overlap between different base tooltips (named structures).

There are two types of minimum energy structures considered in these (and all) calculations. The first is the GLOBAL minimum, the structure that yields to lowest possible energy for a given set of atoms at a given level of theory. The other type is the LOCAL minimum, a structure that is stable (as determined by performing quantum chemical vibrational analyses on the molecule) and will not spontaneously rearrange to a more stable form, but is higher in energy than the GLOBAL minimum (basically, all other stable forms). In theory (although part of the Q-SMAKAS project is to determine the validity of this statement), the "loaded" tooltip is the GLOBAL minimum energy form for these structures and each defect structure is a LOCAL minimum.


Results

1. The "Catalytic Handwaving" of Inverted Hydrogen Atoms
The vast majority of identified defects that have made it through the quantum chemical process are "defective" in the positions of their hydrogen atoms. A typical example is shown for the DC10c tooltip below (click to enlarge.) There are two good reasons why these defects are of little importance to the overall scope of the project and are ignored as part of the "reasonable defect study" despite being important consequences of the methodology.

The "Look where it's going" reason - the most common hydrogen defect finds a C-H group on the outer surface of a tooltip inverting such that the C-H now points INTO the tooltip frame. One could imagine the surprise of reaching for your mouse with your right arm and instead firmly grasping your pancreas. While the tooltips are hollow structures, there is little-to-no room for stray hydrogen atoms.

Which brings up an important question. How do you get that crazy result in the first place?! In Q-SMAKAS, there are two reasons. The first reason is the easier of the two to understand. The temperature of the defect structure MD runs is hovering between 1500 and 4500 K. By comparison, room temperature is 300 K. These tooltips are being subjected to enormous energies, the result of which is the attainment of otherwise chemically inaccessible geometries.

The second reason is methodological. Molecular dynamics simulations require vectors and time steps. The vector is, of course, the momentum of the atom in question. The atoms will bounce around based on the energy of the system and the spring constants holding the atoms together. Simulation balances the energy of the system (wanting to move the atoms around) and the restoring forces of the bonds between atoms. In real life, the forces and motion "just are" and everything "just occurs." We can't calculate the instantaneous forces in the system and monitor it in "real time" because computers aren't yet fast enough. Instead, we perform the recalculations in the simulations at "time steps," intervals between which forces are allowed to do their thing bounded by instances when the system is recalculated to see where everything goes and where everything SHOULD go based on Netwon's Laws of motion.

Imagine driving down the highway - but you can only open your eyes for one second every five seconds. You open your eyes, correct your course, close your eyes and wait/hope, open again, correct, close, etc. If there's a curve in the road, or a barrier, you may over-shoot or under-shoot your course depending on how good your correction is. The situation in the molecular dynamics calculation is a further point complicated from this driving analogue. Because the simulator, and not the Earth, is responsible for all of the physics in the system, you and your car have no concept of what happens during those five second intervals of closed eyes, meaning you're just as likely to drive on a straight path as you would be through a brick wall. In the case of the simulator, since the wall isn't visible in those 5 seconds, you simply pass right through it. The same even CAN happen in the molecular dynamics simulation, where a sufficiently energetic hydrogen will simply "tunnel" through the nucleus of its attached carbon atom and come out the other side to form the new inverted geometry.

The dynamics simulator does not know that a hydrogen atom isn't allowed to pass through a carbon atom, it only knows that the hydrogen atom has a vector that says "go this far a second and await further instructions." Accordingly, a hydrogen can end up in some pretty unconventional positions if (a) the temperature is too high (huge vectors) or (b) the time step is too long (plenty of time move in a direction before the simulator stops it and recalculates the forces). This is a particular problem for the hydrogen because it is the lightest of the atoms, meaning it's also the one that moves around the most.

Our inverted hydrogen atoms, I suspect, are suffering more from the former, which means we could just back off on the time step. Alternatively, we can just increase the mass of the hydrogen atoms so that they don't move as freely as they seem to in the current runs. As our only goal is to find defects, there's nothing wrong with using a little magic in the simulator.

The "Hydrogen? What Hydrogen?" reason - the purpose of hydrogen in the universe is to satisfy unfilled valences in molecules. The goal of quantum chemical studies of dimer tooltips for diamondoid mechanosynthesis is to remove as much of the structure as possible so that the calculations (a) finish as quickly as possible and (b) can be run at as high a level of theory as possible to get the best results possible. The tooltips currently being studied are "bare minimum" structures that contain dimer, binding atoms (Ge or C), the nearest carbon atoms to the binding atoms to recreate the geometry of the binding atom/dimer complex, and hydrogen atoms everywhere an atom has a valence that needs filling. These tooltips, then, are bare minimum structures that run reasonable quickly on reasonable computers. In silico, the hydrogen atoms complete stable electron octets around the peripheral carbon atoms in the tooltip to keep the computational cost down. In reality (it is believed), these tooltips will be the leading edge of a complete tooltip assembly. That is, the hydrogen atoms causing the inverted defect forms simply won't be in the actual product and inverting a C-C group is, well, I don't want to be in the room when a factory full of it happens.

To demonstrate, I've put together some examples of extended tooltip-diamond frameworks in the images below (click to enlarge.) It is quite difficult for the hydrogen to invert when it isn't there. The process of embedding a tooltip into the diamond framework does remove many of the hydrogens in each of the tooltips (although there isn't a single tooltip that can be easily embedded that removes all of the hydrogen atoms in the immediate vicinity of the dimer binding position). The remaining hydrogen atoms in the images of the extended tooltips are shown in green around the blue core tooltip. It should not be assumed that these hydrogen atoms may be the ones to undergo inversion as found in the tooltip calculations. Inversion, as observed in the Q-SMAKAS runs, is an extremely high energy event and is most likely just a methodological artifact. It is, however, conceivable that these nearest hydrogen atoms may migrate on the tooltips under more mundane "defect" circumstances, a possibility that is part of the defect analysis. If this is found to be particularly easy for one tooltip or another, that clearly indicates that a particular tooltip may be too easily broken in a real manufacturing setting and should therefore be considered a possibly unusable structure.

It is worth noting/reminding that ALL of these tooltips are designed to do the same thing - bind and transfer carbon dimers. Accordingly, with a few caveats, the best one wins.

2.

Loaded Tooltips

DC10c

DC10c.L.pdb
Mode: Normal
Energy: -723.28856 Hartrees

DC10c.L.228v.pdb
Mode: Hydrogen Inversion
Energy: -723.13909 Hartrees
Deforms that minimized to the alternate mode:
1 in MD run: 3000 K (0.0407 RMSD) using Deuterium
22 in MD run: 3500 K (0.0348 RMSD)
57 in MD run: 5500 K (0.04 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms


DC10c.L.563v.pdb
Mode: Failure
Energy: -723.16632 Hartrees
Deforms that minimized to the alternate mode:
1 in MD run: 3000 K (0.0407 RMSD) using Deuterium
2 in MD run: 3500 K (0.0348 RMSD)
10 in MD run: 5500 K (0.04 RMSD)

Notes: Tooltip base opened. The second defect (563) is a crack in the tooltip due to a hydrogen migration. This type of defect wouldn't otherwise occur in a Ge-based structure, as the double bond of the dimer binding section stabilizes the otherwise bare carbon. If more of the tool were present this may not happen?


DC10c.L.719v.pdb
Mode: Failure
Energy: -723.28976 Hartrees
Deforms that minimized to the alternate mode:
1 in MD run: 3500 K (0.0348 RMSD)

Notes: Failure mode (10) from the DC10c paper. RHF/3-21G puts this defect 3.15 kJ/mol more stable than the "normal" tooltip. This is one of those potentially problematic cases, as the nature of the defect is not one that's easily deflected by building out the tooltip structure (but there's a lot to be said for the barrier and the higher-level calculations.)


DC10c.L.928r.pdb
Mode: Un-classified
Energy: -723.17136 Hartrees
Deforms that minimized to the alternate mode:
2 in MD run: 3000 K (0.0407 RMSD) using Deuterium
3 in MD run: 5500 K (0.04 RMSD)

Notes: The four carbons opposite the dimer form two extra bonds.

Input files:
DC10c.L.w.tar.gz [1.4 Mb]3000 K (0.0407 RMSD)1000 samples, using Deuterium
DC10c.L.v.tar.gz [1.4 Mb]3500 K (0.0348 RMSD)1000 samples
DC10c.L.r.tar.gz [1.4 Mb]5500 K (0.04 RMSD)1000 samples
Result files:
DC10c.L.results.w.tar.gz [35 Mb]
DC10c.L.results.v.tar.gz [35 Mb]
DC10c.L.results.r.tar.gz [26 Mb]

C100GeATD

1_C100GeATD.L.pdb
Mode: Normal
Energy: -5280.93599 Hartrees

1_C100GeATD.L.929t.pdb
Mode: Hydrogen Inversion
Energy: -5280.68948 Hartrees
Deforms that minimized to the alternate mode:
86 in MD run: 4500 K (~0.04 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms

Here are some "sleek" renderings of this tooltip mode courtesy of Andy Haveland-Robinson

solid

schematic

Andy also rendered some ultra-cool videos of the 4500 Kelvin MD run:
  1_C100GeATD.L.t.spheres.xvid.avi [19 Mb] (2.0 ps at 720x576)
  1_C100GeATD.L.t.solid.xvid.avi [28 Mb] (10.0 ps at 360x288) Xvid.org


1_C100GeATD.L.928t.pdb
Mode:
  Thermally Inaccessible
Energy: -5280.87534 Hartrees
Deforms that minimized to the alternate mode:
1 in MD run: 4500 K (~0.04 RMSD)

Notes: One hydrogen migrated to the dimer. The RHF/STO-3G energy difference from the "normal" tip is 159.51 kJ/mol and represents, by far, the lowest energy defect. Room temperature corresponds to about 2.5 kJ/mol, so we're talking about a defect high above thermal accessibility.

Input files:
1_C100GeATD.L.z.tar.gz [2.9 Mb]1500 K (~0.02 RMSD)1200 samples
1_C100GeATD.L.t.tar.gz [2.3 Mb]4500 K (~0.04 RMSD)1000 samples
Result files:
1_C100GeATD.L.results.z.tar.gz [64 Mb]
1_C100GeATD.L.results.t.tar.gz [27 Mb]

C100GeATS

1_C100GeATS.L.pdb
Mode: Normal
Energy: -5280.98163 Hartrees

1_C100GeATS.L.647t.pdb
Mode: Hydrogen Inversion
Energy: -5280.50773 Hartrees
Deforms that minimized to the alternate mode:
98 in MD run: 4500 K (~0.04 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms

Input files:
1_C100GeATS.L.z.tar.gz [2.8 Mb]1500 K (~0.02 RMSD)1200 samples
1_C100GeATS.L.t.tar.gz [2.3 Mb]4500 K (~0.04 RMSD)1000 samples
Result files:
1_C100GeATS.L.results.z.tar.gz [70 Mb]
1_C100GeATS.L.results.t.tar.gz [34 Mb]

C100GeATSr5

1_C100GeATSr5.L.pdb
Mode: Normal
Energy: -4705.93019 Hartrees

1_C100GeATSr5.L.609v.pdb
Mode: Hydrogen Inversion
Energy: -4705.72718 Hartrees
Deforms that minimized to the alternate mode:
97 in MD run: 3500 K (~0.04 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms


1_C100GeATSr5.L.650v.pdb
Mode: Un-classified
Energy: -4705.82137 Hartrees
Deforms that minimized to the alternate mode:
3 in MD run: 3500 K (~0.04 RMSD)

Notes: Dimer cracked.

Input files:
1_C100GeATSr5.L.z.tar.gz [1.8 Mb]1500 K (~0.02 RMSD)1200 samples
1_C100GeATSr5.L.v.tar.gz [1.5 Mb]3500 K (~0.04 RMSD)1000 samples
Result files:
1_C100GeATSr5.L.results.z.tar.gz [34 Mb]
1_C100GeATSr5.L.results.v.tar.gz [31 Mb]

C100GeATSr6

1_C100GeATSr6.L.pdb
Mode: Normal
Energy: -4861.28173 Hartrees

1_C100GeATSr6.L.424x.pdb
Mode: Hydrogen Inversion
Energy: -4861.18756 Hartrees
Deforms that minimized to the alternate mode:
107 in MD run: 2500 K (~0.04 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms

Input files:
1_C100GeATSr6.L.z.tar.gz [2.3 Mb]1500 K (~0.02 RMSD)1200 samples
1_C100GeATSr6.L.x.tar.gz [1.8 Mb]2500 K (~0.04 RMSD)1000 samples
Result files:
1_C100GeATSr6.L.results.z.tar.gz [50 Mb]
1_C100GeATSr6.L.results.x.tar.gz [43 Mb]

C100GeCTS

1_C100GeCTS.L.pdb
Mode: Normal
Energy: -5283.35665 Hartrees

1_C100GeCTS.L.123v.pdb
Mode: Hydrogen Inversion
Energy: -5283.11688 Hartrees
Deforms that minimized to the alternate mode:
26 in MD run: 3500 K (~0.04 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms


1_C100GeCTS.L.511v.pdb
Mode: Un-classified
Energy: -5283.17814 Hartrees
Deforms that minimized to the alternate mode:
4 in MD run: 3500 K (~0.04 RMSD)

Notes: Dimer cracked.

Input files:
1_C100GeCTS.L.v.tar.gz [2.4 Mb]3500 K (~0.04 RMSD)1000 samples
Result files:
1_C100GeCTS.L.results.v.tar.gz [63.5 Mb]

C110GeS

2_C110GeS.L.pdb
Mode: Normal
Energy: -4515.37009 Hartrees

2_C110GeS.L.424v.pdb
Mode: Hydrogen Inversion
Energy: -4515.23658 Hartrees
Deforms that minimized to the alternate mode:
102 in MD run: 3500 K (~0.04 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms


2_C110GeS.L.558v.pdb
Mode: Un-classified
Energy: -4515.36179 Hartrees
Deforms that minimized to the alternate mode:
2 in MD run: 3500 K (~0.04 RMSD)

Notes: Dimer cracked.


2_C110GeS.L.783v.pdb
Mode: Un-classified
Energy: -4515.32278 Hartrees
Deforms that minimized to the alternate mode:
1 in MD run: 3500 K (~0.04 RMSD)

Notes: Dimer is held by just one carbon with the other carbon bound back onto the tooltip; and a hydrogen migrated to the dimer.


2_C110GeS.L.613v.pdb
Mode: Un-classified
Energy: -4515.37451 Hartrees
Deforms that minimized to the alternate mode:
1 in MD run: 3500 K (~0.04 RMSD)

Notes: Dimer is held by just one carbon with the other carbon protruding from the tooltip; lower energy mode.

Input files:
2_C110GeS.L.z.tar.gz [1.6 Mb]1500 K (~0.02 RMSD)1200 samples
2_C110GeS.L.v.tar.gz [1.2 Mb]3500 K (~0.04 RMSD)1000 samples
Result files:
2_C110GeS.L.results.z.tar.gz [30 Mb]
2_C110GeS.L.results.v.tar.gz [24 Mb]

C111Ge3

3_C111Ge3.L.pdb
Mode: Normal
Energy: -4821.32074 Hartrees

3_C111Ge3.L.814v.pdb
Mode: Hydrogen Inversion
Energy: -4821.18351 Hartrees
Deforms that minimized to the alternate mode:
84 in MD run: 3500 K (~0.04 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms

Input files:
3_C111Ge3.L.z.tar.gz [2.1 Mb]1500 K (~0.02 RMSD)1200 samples
3_C111Ge3.L.v.tar.gz [1.6 Mb]2500 K (~0.04 RMSD)1000 samples
Result files:
3_C111Ge3.L.results.z.tar.gz [37 Mb]
3_C111Ge3.L.results.v.tar.gz [31 Mb]

AdamGe22

5_AdamGe22.L.pdb
Mode: Normal
Energy: -4476.53808 Hartrees

5_AdamGe22.L.23v.pdb
Mode: Hydrogen Inversion
Energy: -4476.38819 Hartrees
Deforms that minimized to the alternate mode:
1 in MD run: 3000 K (0.0402 RMSD) using Deuterium
6 in MD run: 3500 K (0.0413 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms


5_AdamGe22.L.208v.pdb
Mode: Failure
Energy: -4476.53628 Hartrees
Deforms that minimized to the alternate mode:
1 in MD run: 3000 K (0.0402 RMSD) using Deuterium
1 in MD run: 3500 K (0.0413 RMSD)

Notes: Strange dimer orientation. A subtle structural change, with the arrangement of the dimer on the Ge...Ge displaced from its preferred arrangement, which is only 4.7 kJ/mol lower (so we're looking at wobble in the dimer.)


5_AdamGe22.L.787v.pdb
Mode: Failure
Energy: -4476.52428 Hartrees
Deforms that minimized to the alternate mode:
1 in MD run: 3000 K (0.0402 RMSD) using Deuterium
1 in MD run: 3500 K (0.0413 RMSD)

Notes: Tooltip opened up. A hydrogen migration from a Ge atom onto the nearest C atom that the Ge was attached to. The Ge-C bond is broken with this.


5_AdamGe22.L.701w.pdb
Mode: Failure
Energy: -4476.54660 Hartrees
Deforms that minimized to the alternate mode:
1 in MD run: 3000 K (0.0402 RMSD) using Deuterium

Notes: Dimer is held by just one carbon with the other carbon protruding from the tooltip. Lower energy configuration.

Input files:
5_AdamGe22.L.z.tar.gz [1.4 Mb]1500 K (~0.02 RMSD)1200 samples
5_AdamGe22.L.v.tar.gz [1.1 Mb]3500 K (0.0413 RMSD)1000 samples
5_AdamGe22.L.w.tar.gz [1.2 Mb]3000 K (0.0402 RMSD) with Deuterium1000 samples
Result files:
5_AdamGe22.L.results.z.tar.gz [28 Mb]
5_AdamGe22.L.results.v.tar.gz [24 Mb]
5_AdamGe22.L.results.w.tar.gz [25 Mb]

AdamGe33

5_AdamGe33.L.pdb
Mode: "Normal"
Energy: -4477.73819 Hartrees

5_AdamGe33.L.275w.pdb
Mode: Hydrogen Inversion
Energy: -4477.66301 Hartrees
Deforms that minimized to the alternate mode:
93 in MD run: 3000 K (~0.04 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms


5_AdamGe33.L.568w.pdb
Mode: Un-classified
Energy: -4477.70575 Hartrees
Deforms that minimized to the alternate mode:
2 in MD run: 3000 K (~0.04 RMSD)

Notes: Dimer cracked.

Input files:
5_AdamGe33.L.z.tar.gz [1.6 Mb]1500 K (~0.02 RMSD)1200 samples
5_AdamGe33.L.w.tar.gz [1.3 Mb]3000 K (~0.04 RMSD)1000 samples
Result files:
5_AdamGe33.L.results.z.tar.gz [29 Mb]
5_AdamGe33.L.results.w.tar.gz [24 Mb]

Diad3Ge22

6_Diad3Ge22.L.pdb
Mode: Normal
Energy: -4783.59967 Hartrees

6_Diad3Ge22.L.911x.pdb
Mode: Hydrogen Inversion
Energy: -4783.42440 Hartrees
Deforms that minimized to the alternate mode:
3 in MD run: 2500 K (~0.04 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms


6_Diad3Ge22.L.152x.pdb
Mode: Failure
Energy: -4783.59967 Hartrees
Deforms that minimized to the alternate mode:
?? in MD run: 2500 K (0.04 RMSD)

Notes: Dimer changed sides. Detected via visual inspection of the MD run - otherwise it has the same energy as the normal mode tooltip. Counting them would be a pain because we'd have to do it manually/visually since energy gives no indication. Come up with an algorithm to automatically detect symmetric deformations by labeling a select triple/triangle of atoms, then do a fuzzy comparison to atoms in the immediate vicinity of the binding pocket.


6_Diad3Ge22.L.156x.pdb
Mode: Failure
Energy: -4783.55711 Hartrees
Deforms that minimized to the alternate mode:
13 in MD run: 2500 K (0.04 RMSD)

Notes: Dimer twisted. One class of interesting defect here, resulting from the flip of one of the H atoms on the Ge...Ge region (so the H's aren't both on the same side anymore). The flipped case is 111 kJ/mol higher than the "normal" tooltip.

Input files:
6_Diad3Ge22.L.z.tar.gz [2.1 Mb]1500 K (~0.02 RMSD)1200 samples
6_Diad3Ge22.L.x.tar.gz [1.7 Mb]3000 K (~0.04 RMSD)1000 samples
Result files:
6_Diad3Ge22.L.results.z.tar.gz [39 Mb]
6_Diad3Ge22.L.results.x.tar.gz [31 Mb]

DCB55AGe

7_DCB55AGe.L.pdb
Mode: Normal
Energy: -4743.61394 Hartrees

7_DCB55AGe.L.953t.pdb
Mode: Hydrogen Inversion
Energy: -4743.54031 Hartrees
Deforms that minimized to the alternate mode:
121 in MD run: 4500 K (~0.04 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms

Input files:
7_DCB55AGe.L.z.tar.gz [1.9 Mb]1500 K (~0.02 RMSD)1200 samples
7_DCB55AGe.L.t.tar.gz [1.5 Mb]4500 K (~0.04 RMSD)1000 samples
Result files:
7_DCB55AGe.L.results.z.tar.gz [34 Mb]
7_DCB55AGe.L.results.t.tar.gz [19 Mb]

DCB55BGe

7_DCB55BGe.L.pdb
Mode: Normal
Energy: -4821.21100 Hartrees

7_DCB55BGe.L.432w.pdb
Mode: Hydrogen Inversion
Energy: -4743.54031 Hartrees
Deforms that minimized to the alternate mode:
20 in MD run: 3000 K (0.0610 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms

Input files:
7_DCB55BGe.L.z.tar.gz [2.1 Mb]1500 K (0.0432 RMSD)1200 samples
7_DCB55BGe.L.w.tar.gz [1.7 Mb]3000 K (0.0610 RMSD)1000 samples
Result files:
7_DCB55BGe.L.results.z.tar.gz [39 Mb]
7_DCB55BGe.L.results.w.tar.gz [10 Mb]

DCB55CGe

7_DCB55CGe.L.pdb
Mode: Normal
Energy: -4858.87444 Hartrees

7_DCB55CGe.L.398v.pdb
Mode: Hydrogen Inversion
Energy: -4858.84148 Hartrees
Deforms that minimized to the alternate mode:
125 in MD run: 3500 K (~0.04 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms

Input files:
7_DCB55CGe.L.z.tar.gz [2.1 Mb]1500 K (~0.02 RMSD)1200 samples
7_DCB55CGe.L.v.tar.gz [1.7 Mb]3500 K (~0.04 RMSD)1000 samples
Result files:
7_DCB55CGe.L.results.z.tar.gz [38 Mb]
7_DCB55CGe.L.results.v.tar.gz [23 Mb]

DCB57Ge

7_DCB57Ge.L.pdb
Mode: Normal
Energy: -4858.94184 Hartrees

7_DCB57Ge.L.573u.pdb
Mode: Hydrogen Inversion
Energy: -4858.85935 Hartrees
Deforms that minimized to the alternate mode:
132 in MD run: 4000 K (~0.04 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms


7_DCB57Ge.L.302u.pdb
Mode: Un-classified
Energy: -4858.76556 Hartrees
Deforms that minimized to the alternate mode:
2 in MD run: 4000 K (~0.04 RMSD)

Notes: Dimer cracked.

Input files:
7_DCB57Ge.L.z.tar.gz [2.1 Mb]1500 K (~0.02 RMSD)1200 samples
7_DCB57Ge.L.u.tar.gz [1.7 Mb]4000 K (~0.04 RMSD)1000 samples
Result files:
7_DCB57Ge.L.results.z.tar.gz [36 Mb]
7_DCB57Ge.L.results.u.tar.gz [26 Mb]

DCB5Ge

7_DCB5Ge.L.pdb
Mode: Normal
Energy: -4821.22754 Hartrees

7_DCB5Ge.L.445w.pdb
Mode: Hydrogen Inversion
Energy: -4821.17042 Hartrees
Deforms that minimized to the alternate mode:
46 in MD run: 3000 K (~0.04 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms


7_DCB5Ge.L.111w.pdb
Mode: Un-classified
Energy: -4821.17948 Hartrees
Deforms that minimized to the alternate mode:
1 in MD run: 3000 K (~0.04 RMSD)

Notes: Dimer cracked.

Input files:
7_DCB5Ge.L.z.tar.gz [2.1 Mb]1500 K (~0.02 RMSD)1200 samples
7_DCB5Ge.L.w.tar.gz [1.7 Mb]4000 K (~0.04 RMSD)1000 samples
Result files:
7_DCB5Ge.L.results.z.tar.gz [43 Mb]
7_DCB5Ge.L.results.w.tar.gz [31 Mb]

DCB65Ge

7_DCB65Ge.L.pdb
Mode: Normal
Energy: -4936.57467 Hartrees

7_DCB65Ge.L.320w.pdb
Mode: Hydrogen Inversion
Energy: -4936.54164 Hartrees
Deforms that minimized to the alternate mode:
192 in MD run: 3000 K (~0.04 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms

Input files:
7_DCB65Ge.L.z.tar.gz [2.3 Mb]1500 K (~0.02 RMSD)1200 samples
7_DCB65Ge.L.w.tar.gz [1.9 Mb]3000 K (~0.04 RMSD)1000 samples
Result files:
7_DCB65Ge.L.results.z.tar.gz [43 Mb]
7_DCB65Ge.L.results.w.tar.gz [30 Mb]

DCB6Ge

7_DCB6Ge.L.pdb
Mode: Normal
Energy: -4897.82372 Hartrees

7_DCB6Ge.L.711v.pdb
Mode: Hydrogen Inversion
Energy: -4897.70150 Hartrees
Deforms that minimized to the alternate mode:
90 in MD run: 3500 K (~0.04 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms

Input files:
7_DCB6Ge.L.z.tar.gz [2.2 Mb]1500 K (~0.02 RMSD)1200 samples
7_DCB6Ge.L.v.tar.gz [1.8 Mb]3500 K (~0.04 RMSD)1000 samples
Result files:
7_DCB6Ge.L.results.z.tar.gz [39 Mb]
7_DCB6Ge.L.results.v.tar.gz [30 Mb]

DCB75Ge

7_DCB75Ge.L.pdb
Mode: Normal
Energy: -5096.09096 Hartrees

7_DCB75Ge.L.274x.pdb
Mode: Hydrogen Inversion
Energy: -5096.03682 Hartrees
Deforms that minimized to the alternate mode:
478 in MD run: 2500 K (~0.04 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms

Input files:
7_DCB75Ge.L.z.tar.gz [2.8 Mb]1500 K (~0.02 RMSD)1200 samples
7_DCB75Ge.L.x.tar.gz [2.3 Mb]2500 K (~0.04 RMSD)1000 samples
Result files:
7_DCB75Ge.L.results.z.tar.gz [81 Mb]
7_DCB75Ge.L.results.x.tar.gz [56 Mb]

DCBIceane7Ge

7_DCBIceane7Ge.L.pdb
Mode: Normal
Energy: -4630.64599 Hartrees

7_DCBIceane7Ge.L.721v.pdb
Mode: Hydrogen Inversion
Energy: -4630.55822 Hartrees
Deforms that minimized to the alternate mode:
31 in MD run: 3500 K (~0.04 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms

Input files:
7_DCBIceane7Ge.L.z.tar.gz [1.8 Mb]1500 K (~0.02 RMSD)1200 samples
7_DCBIceane7Ge.L.v.tar.gz [1.4 Mb]3500 K (~0.04 RMSD)1000 samples
Result files:
7_DCBIceane7Ge.L.results.z.tar.gz [34 Mb]
7_DCBIceane7Ge.L.results.v.tar.gz [30 Mb]

MCB55Ge

7_MCB55Ge.L.pdb
Mode: Normal
Energy: -4860.03772 Hartrees

7_MCB55Ge.L.373q.pdb
Mode: Hydrogen Inversion
Energy: -4860.08270 Hartrees
Deforms that minimized to the alternate mode:
10 in MD run: 6000 K (~0.08 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms


7_MCB55Ge.L.800q.pdb
Mode: Un-classified
Energy: -4860.13980 Hartrees
Deforms that minimized to the alternate mode:
3 in MD run: 6000 K (~0.08 RMSD)

Notes: Hydrogen re-arrangement to significantly lower energy - suggests a new, more stable tooltip?

Input files:
7_MCB55Ge.L.z.tar.gz [2.2 Mb]1500 K (~0.04 RMSD)1200 samples
7_MCB55Ge.L.q.tar.gz [1.8 Mb]6000 K (~0.08 RMSD)1000 samples
Result files:
7_MCB55Ge.L.results.z.tar.gz [47 Mb]
7_MCB55Ge.L.results.q.tar.gz [1.1 Mb] (Very few deforms actually converged during energy minimization.)

MCB57Ge

7_MCB57Ge.L.pdb
Mode: Normal
Energy: -5020.12499 Hartrees
No alternate modes detected.
Input files:
7_MCB57Ge.L.z.tar.gz [2.2 Mb]1500 K (~0.02 RMSD)1200 samples
7_MCB57Ge.L.y.tar.gz [2.1 Mb]2000 K (~0.04 RMSD)1000 samples
Result files:
7_MCB57Ge.L.results.z.tar.gz [77 Mb]
7_MCB57Ge.L.results.y.tar.gz [62 Mb]

MCB5Ge

7_MCB5Ge.L.pdb
Mode: Normal
Energy: -4942.96584 Hartrees
No alternate modes detected.
Input files:
7_MCB5Ge.L.z.tar.gz [2.5 Mb]1500 K (~0.02 RMSD)1200 samples
7_MCB5Ge.L.y.tar.gz [2.0 Mb]2000 K (~0.04 RMSD)1000 samples
Result files:
7_MCB5Ge.L.results.z.tar.gz [70 Mb]
7_MCB5Ge.L.results.y.tar.gz [61 Mb]

MCB75Ge

7_MCB75Ge.L.pdb
Mode: Normal
Energy: -5361.42659 Hartrees
No alternate modes detected.
Input files:
7_MCB75Ge.L.x.tar.gz [2.7 Mb]2500 K (~0.04 RMSD)1000 samples
Result files:
7_MCB75Ge.L.results.x.tar.gz [17 Mb] (Very few deforms actually converged during energy minimization.)

TwistaneGe

9_TwistaneGe.L.pdb
Mode: Normal
Energy: -4438.88518 Hartrees

9_TwistaneGe.L.701u.pdb
Mode: Hydrogen Inversion
Energy: -4438.73775 Hartrees
Deforms that minimized to the alternate mode:
24 in MD run: 4000 K (~0.04 RMSD)

Notes: One or more hydrogens inverted to inside the tooltip. Ignorable - see The "Catalytic Handwaving" of Inverted Hydrogen Atoms


9_TwistaneGe.L.238u.pdb
Mode: Un-classified
Energy: -4438.90750 Hartrees
Deforms that minimized to the alternate mode:
1 in MD run: 4000 K (~0.04 RMSD)

Notes: Dimer is held by just one carbon with the other carbon protruding from the tooltip and with a hydrogen attached. Lower energy configuration.

Input files:
9_TwistaneGe.L.z.tar.gz [1.4 Mb]1500 K (~0.02 RMSD)1200 samples
9_TwistaneGe.L.u.tar.gz [1.1 Mb]4000 K (~0.04 RMSD)1000 samples
Result files:
9_TwistaneGe.L.results.z.tar.gz [26 Mb]
9_TwistaneGe.L.results.u.tar.gz [18 Mb]


Copyright © 2014 Helcorp Interactive, LLC. All rights reserved.
NanoHive-1, NanoHive@Home, and their logos are trademarks of Helcorp Interactive, LLC.