Discovering chemistry with an ab initio nanoreactor

Chemical understanding is driven by the experimental discovery of new compounds and reactivity, and is supported by theory and computation that provides detailed physical insight. While theoretical and computational studies have generally focused on specific processes or mechanistic hypotheses, recent methodological and computational advances harken the advent of their principal role in discovery. Here we report the development and application of the ab initio nanoreactor – a highly accelerated, first-principles molecular dynamics simulation of chemical reactions that discovers new molecules and mechanisms without preordained reaction coordinates or elementary steps. Using the nanoreactor we show new pathways for glycine synthesis from primitive compounds proposed to exist on the early Earth, providing new insight into the classic Urey-Miller experiment. These results highlight the emergence of theoretical and computational chemistry as a tool for discovery in addition to its traditional role of interpreting experimental findings.

[1]  H. Urey,et al.  Organic compound synthesis on the primitive earth. , 1959, Science.

[2]  Stanley L. Miller,et al.  Organic Compound Synthes on the Primitive Eart: Several questions about the origin of life have been answered, but much remains to be studied , 1959 .

[3]  V. R. Saunders,et al.  A “Level–Shifting” method for converging closed shell Hartree–Fock wave functions , 1973 .

[4]  P. Pulay Convergence acceleration of iterative sequences. the case of scf iteration , 1980 .

[5]  J. Kasting,et al.  Earth's early atmosphere , 1987, Science.

[6]  K. Aoki,et al.  FT-IR Study of the Solid State Polymerization of Acetylene under Pressure , 1996 .

[7]  A. Becke Density-functional thermochemistry. , 1996 .

[8]  Michele Parrinello,et al.  Solid-State Polymerization of Acetylene under Pressure: Ab Initio Simulation , 1997 .

[9]  K. Peterson,et al.  An examination of intrinsic errors in electronic structure methods using the Environmental Molecular Sciences Laboratory computational results database and the Gaussian-2 set , 1998 .

[10]  A. Wahner,et al.  Gas‐phase reaction of N2O5 with water vapor: Importance of heterogeneous hydrolysis of N2O5 and surface desorption of HNO3 in a large Teflon chamber , 1998 .

[11]  Mark A. Miller,et al.  Archetypal energy landscapes , 1998, Nature.

[12]  Martin Head-Gordon,et al.  Benchmark variational coupled cluster doubles results , 2000 .

[13]  J. Badding,et al.  Solid State Polymerization of Acetylene at High Pressure and Low Temperature , 2000 .

[14]  Alessandro Laio,et al.  Efficient exploration of reactive potential energy surfaces using Car-Parrinello molecular dynamics. , 2003, Physical review letters.

[15]  Andrew G. Leach,et al.  A Standard Set of Pericyclic Reactions of Hydrocarbons for the Benchmarking of Computational Methods: The Performance of ab Initio, Density Functional, CASSCF, CASPT2, and CBS-QB3 Methods for the Prediction of Activation Barriers, Reaction Energetics, and Transition State Geometries , 2003 .

[16]  A. Chakraborty,et al.  A growing string method for determining transition states: comparison to the nudged elastic band and string methods. , 2004, The Journal of chemical physics.

[17]  Gabor A. Somorjai,et al.  Formation of Hollow Nanocrystals Through the Nanoscale Kirkendall Effect , 2004, Science.

[18]  Bernd Ensing,et al.  Metadynamics as a Tool for Exploring Free Energy Landscapes of Chemical Reactions , 2006 .

[19]  Martin Head-Gordon,et al.  Advances in Methods and Algorithms in a Modern Quantum Chemistry Program Package , 2006 .

[20]  Pedro Alexandrino Fernandes,et al.  General performance of density functionals. , 2007, The journal of physical chemistry. A.

[21]  H. James Cleaves,et al.  A Reassessment of Prebiotic Organic Synthesis in Neutral Planetary Atmospheres , 2008, Origins of Life and Evolution of Biospheres.

[22]  F. Matthias Bickelhaupt,et al.  Energy landscapes of nucleophilic substitution reactions: A comparison of density functional theory and coupled cluster methods , 2007, J. Comput. Chem..

[23]  Aric Hagberg,et al.  Exploring Network Structure, Dynamics, and Function using NetworkX , 2008, Proceedings of the Python in Science Conference.

[24]  Ivan S Ufimtsev,et al.  Quantum Chemistry on Graphical Processing Units. 1. Strategies for Two-Electron Integral Evaluation. , 2008, Journal of chemical theory and computation.

[25]  Branko Ruscic,et al.  High-accuracy extrapolated ab initio thermochemistry. III. Additional improvements and overview. , 2008, The Journal of chemical physics.

[26]  Gaël Varoquaux,et al.  Proceedings of the 20th Python in Science Conference 2021 (SciPy 2021), Virtual Conference, July 12 - July 18, 2021 , 2008, SciPy.

[27]  Todd J. Martinez,et al.  Graphical Processing Units for Quantum Chemistry , 2008, Computing in Science & Engineering.

[28]  Ivan S Ufimtsev,et al.  Quantum Chemistry on Graphical Processing Units. 3. Analytical Energy Gradients, Geometry Optimization, and First Principles Molecular Dynamics. , 2009, Journal of chemical theory and computation.

[29]  E. Reed,et al.  Ab initio simulation of the equation of state and kinetics of shocked water. , 2009, The Journal of chemical physics.

[30]  Ivan S Ufimtsev,et al.  Quantum Chemistry on Graphical Processing Units. 2. Direct Self-Consistent-Field Implementation. , 2009, Journal of chemical theory and computation.

[31]  B. A. Lindquist,et al.  Photodynamics in complex environments: ab initio multiple spawning quantum mechanical/molecular mechanical dynamics. , 2009, The journal of physical chemistry. B.

[32]  Weitao Yang,et al.  Accelerating self-consistent field convergence with the augmented Roothaan-Hall energy function. , 2010, The Journal of chemical physics.

[33]  Amitesh Maiti,et al.  Synthesis of glycine-containing complexes in impacts of comets on early Earth. , 2010, Nature chemistry.

[34]  Douglas Thain,et al.  Work Queue + Python: A Framework For Scalable Scientific Ensemble Applications , 2011 .

[35]  Todd J. Martínez,et al.  Charge Transfer and Polarization in Solvated Proteins from Ab Initio Molecular Dynamics , 2011 .

[36]  Gaël Varoquaux,et al.  Scikit-learn: Machine Learning in Python , 2011, J. Mach. Learn. Res..

[37]  Koichi M. T. Yamada,et al.  Interstellar molecules : their laboratory and interstellar habitat , 2011 .

[38]  K. Menten,et al.  Molecules Detected in Interstellar Space , 2011 .

[39]  Christine M. Isborn,et al.  Excited-State Electronic Structure with Configuration Interaction Singles and Tamm–Dancoff Time-Dependent Density Functional Theory on Graphical Processing Units , 2011, Journal of chemical theory and computation.

[40]  B. Viskolcz,et al.  Chemical evolution of biomolecule building blocks. Can thermodynamics explain the accumulation of glycine in the prebiotic ocean? , 2011, Physical chemistry chemical physics : PCCP.

[41]  Fabio Pietrucci,et al.  Graph theory meets ab initio molecular dynamics: atomic structures and transformations at the nanoscale. , 2011, Physical review letters.

[42]  Ivan S. Ufimtsev,et al.  Dynamic Precision for Electron Repulsion Integral Evaluation on Graphical Processing Units (GPUs). , 2011, Journal of chemical theory and computation.

[43]  Edward F. Valeev,et al.  Prediction of Reaction Barriers and Thermochemical Properties with Explicitly Correlated Coupled-Cluster Methods: A Basis Set Assessment. , 2012, Journal of chemical theory and computation.

[44]  A comparison of density functional theory and coupled cluster methods for the calculation of electric dipole polarizability gradients of methane , 2012 .

[45]  Satoshi Maeda,et al.  Toward Predicting Full Catalytic Cycle Using Automatic Reaction Path Search Method: A Case Study on HCo(CO)3-Catalyzed Hydroformylation. , 2012, Journal of chemical theory and computation.

[46]  Pathways for the Formation and Evolution of Peptides in Prebiotic Environments , 2012 .

[47]  Heather J Kulik,et al.  Ab initio quantum chemistry for protein structures. , 2012, The journal of physical chemistry. B.

[48]  P. Wipf,et al.  Stochastic voyages into uncharted chemical space produce a representative library of all possible drug-like compounds. , 2013, Journal of the American Chemical Society.

[49]  Todd J. Martínez,et al.  Generating Efficient Quantum Chemistry Codes for Novel Architectures. , 2013, Journal of chemical theory and computation.

[50]  Paul M. Zimmerman,et al.  Automated discovery of chemically reasonable elementary reaction steps , 2013, J. Comput. Chem..

[51]  Cooper J. Galvin,et al.  Complex Chemical Reaction Networks from Heuristics-Aided Quantum Chemistry. , 2014, Journal of chemical theory and computation.