Active Learning Accelerates Ab Initio Molecular Dynamics on Pericyclic Reactive Energy Surfaces

2020 ◽  
Author(s):  
Shi Jun Ang ◽  
Wujie Wang ◽  
Daniel Schwalbe-Koda ◽  
Simon Axelrod ◽  
Rafael Gomez-Bombarelli
Keyword(s):  
Molecular Dynamics ◽  
Potential Energy ◽  
Transition State ◽  
Potential Energy Surface ◽  
Ab Initio ◽  
Energy Surface ◽  
Computational Cost ◽  
Reactive Trajectories ◽  
Dynamics Simulations ◽  
Energy Surfaces

<div>Modeling dynamical effects in chemical reactions, such as post-transition state bifurcation, requires <i>ab initio</i> molecular dynamics simulations due to the breakdown of simpler static models like transition state theory. However, these simulations tend to be restricted to lower-accuracy electronic structure methods and scarce sampling because of their high computational cost. Here, we report the use of statistical learning to accelerate reactive molecular dynamics simulations by combining high-throughput ab initio calculations, graph-convolution interatomic potentials and active learning. This pipeline was demonstrated on an ambimodal trispericyclic reaction involving 8,8-dicyanoheptafulvene and 6,6-dimethylfulvene. With a dataset size of approximately</div><div>31,000 M062X/def2-SVP quantum mechanical calculations, the computational cost of exploring the reactive potential energy surface was reduced by an order of magnitude. Thousands of virtually costless picosecond-long reactive trajectories suggest that post-transition state bifurcation plays a minor role for the reaction in vacuum. Furthermore, a transfer-learning strategy effectively upgraded the potential energy surface to higher</div><div>levels of theory ((SMD-)M06-2X/def2-TZVPD in vacuum and three other solvents, as well as the more accurate DLPNO-DSD-PBEP86 D3BJ/def2-TZVPD) using about 10% additional calculations for each surface. Since the larger basis set and the dynamic correlation capture intramolecular non-covalent interactions more accurately, they uncover longer lifetimes for the charge-separated intermediate on the more accurate potential energy surfaces. The character of the intermediate switches from entropic to thermodynamic upon including implicit solvation effects, with lifetimes increasing with solvent polarity. Analysis of 2,000 reactive trajectories on the chloroform PES shows a qualitative agreement with the experimentally-reported periselectivity for this reaction. This overall approach is broadly applicable and opens a door to the study of dynamical effects in larger, previously-intractable reactive systems.</div>

Active Learning Accelerates Ab Initio Molecular Dynamics on Pericyclic Reactive Energy Surfaces

2020 ◽  
Author(s):  
Shi Jun Ang ◽  
Wujie Wang ◽  
Daniel Schwalbe-Koda ◽  
Simon Axelrod ◽  
Rafael Gomez-Bombarelli
Keyword(s):  
Molecular Dynamics ◽  
Potential Energy ◽  
Transition State ◽  
Potential Energy Surface ◽  
Ab Initio ◽  
Energy Surface ◽  
Computational Cost ◽  
Reactive Trajectories ◽  
Dynamics Simulations ◽  
Energy Surfaces

<div>Modeling dynamical effects in chemical reactions, such as post-transition state bifurcation, requires <i>ab initio</i> molecular dynamics simulations due to the breakdown of simpler static models like transition state theory. However, these simulations tend to be restricted to lower-accuracy electronic structure methods and scarce sampling because of their high computational cost. Here, we report the use of statistical learning to accelerate reactive molecular dynamics simulations by combining high-throughput ab initio calculations, graph-convolution interatomic potentials and active learning. This pipeline was demonstrated on an ambimodal trispericyclic reaction involving 8,8-dicyanoheptafulvene and 6,6-dimethylfulvene. With a dataset size of approximately</div><div>31,000 M062X/def2-SVP quantum mechanical calculations, the computational cost of exploring the reactive potential energy surface was reduced by an order of magnitude. Thousands of virtually costless picosecond-long reactive trajectories suggest that post-transition state bifurcation plays a minor role for the reaction in vacuum. Furthermore, a transfer-learning strategy effectively upgraded the potential energy surface to higher</div><div>levels of theory ((SMD-)M06-2X/def2-TZVPD in vacuum and three other solvents, as well as the more accurate DLPNO-DSD-PBEP86 D3BJ/def2-TZVPD) using about 10% additional calculations for each surface. Since the larger basis set and the dynamic correlation capture intramolecular non-covalent interactions more accurately, they uncover longer lifetimes for the charge-separated intermediate on the more accurate potential energy surfaces. The character of the intermediate switches from entropic to thermodynamic upon including implicit solvation effects, with lifetimes increasing with solvent polarity. Analysis of 2,000 reactive trajectories on the chloroform PES shows a qualitative agreement with the experimentally-reported periselectivity for this reaction. This overall approach is broadly applicable and opens a door to the study of dynamical effects in larger, previously-intractable reactive systems.</div>

scholarly journals Active Learning Accelerates Ab Initio Molecular Dynamics on Pericyclic Reactive Energy Surfaces

2020 ◽  
Author(s):  
Shi Jun Ang ◽  
Wujie Wang ◽  
Daniel Schwalbe-Koda ◽  
Simon Axelrod ◽  
Rafael Gomez-Bombarelli
Keyword(s):  
Molecular Dynamics ◽  
Active Learning ◽  
Potential Energy ◽  
Transition State ◽  
Potential Energy Surface ◽  
Ab Initio ◽  
Energy Surface ◽  
Computational Cost ◽  
Ab Initio Molecular Dynamics ◽  
Dynamics Simulations

Modeling dynamical effects in chemical reactions, such as post-transition state bifurcation, requires ab initio molecular dynamics simulations due to the breakdown of simpler static models like transition state theory. However, these simulations tend to be restricted to lower-accuracy electronic structure methods and scarce sampling because of their high computational cost. Here, we report the use of statistical learning to accelerate reactive molecular dynamics simulations by combining high-throughput ab initio calculations, graph-convolution interatomic potentials and active learning. This pipeline was demonstrated on an ambimodal trispericyclic reaction involving 8,8-dicyanoheptafulvene and 6,6-dimethylfulvene. With a training dataset size of approximately 31,000 M062X/def2-SVP quantum mechanical calculations, the computational cost of exploring the reactive potential energy surface was reduced by an order of magnitude. These tools are then used to automatically predict a reaction mechanism that is in agreement with the experimentally-reported product distribution. In addition, thousands of virtually costless picosecond-long reactive trajectories suggest that post-transition state bifurcation plays a very minor role in the reaction. Furthermore, a transfer-learning strategy effectively allowed to upgrade the potential energy surface to higher levels of theory (def2-TZVPD basis set and double hybrid functional) using less than 10% additional calculations. Since these methods capture intramolecular non-covalent interactions more accurately, they uncover longer lifetimes for entropic intermediates. This overall approach is broadly applicable and opens the door to the study of dynamical effects in larger, previously-intractable reactive systems.

Dynamics and kinetics of the Si(1D) + H2/D2 reactions on a new global ab initio potential energy surface

2021 ◽  
Vol 23 (10) ◽  
pp. 6141-6153
Author(s):  
Jianwei Cao ◽  
Yanan Wu ◽  
Haitao Ma ◽  
Zhitao Shen ◽  
Wensheng Bian
Keyword(s):  
Molecular Dynamics ◽  
Potential Energy ◽  
Potential Energy Surface ◽  
Ab Initio ◽  
Quantum Dynamics ◽  
Energy Surface ◽  
Molecular Dynamics Calculations ◽  
Ring Polymer ◽  
Ring Polymer Molecular Dynamics ◽  
Kinetics Of

Quantum dynamics and ring polymer molecular dynamics calculations reveal interesting dynamical and kinetic behaviors of an endothermic complex-forming reaction.

scholarly journals A [1 + 2] cycloaddition instead of usual [2 + 3] cycloaddition between the B12N12 cluster and methyl azide: Potential energy surface calculations and Born–Oppenheimer molecular dynamics simulations

2020 ◽  
Vol 45 ◽  
pp. 146867831990058
Author(s):  
Parvaneh Pakravan ◽  
Seyyed Amir Siadati
Keyword(s):  
Molecular Dynamics ◽  
Potential Energy ◽  
Potential Energy Surface ◽  
Molecular Dynamics Simulations ◽  
Energy Surface ◽  
Cycloaddition Reaction ◽  
Reliable Data ◽  
Reaction Channel ◽  
Dynamics Simulations

We have examined here the possibility of functionalization of the B12N12 cluster by methyl azide by means of a [2 + 3] cycloaddition reaction in analogy with the spontaneous functionalization of C20 fullerene using the same reaction. To achieve more reliable data, all possible interactions at different positions and orientations were considered by reaction channel study and potential energy surface calculations. Also, Born–Oppenheimer molecular dynamics simulations were used to find probable species which could emerge during the reactions.

Collisional excitation of CH(X2Π) by He: new ab initio potential energy surfaces and scattering calculations

2015 ◽  
Vol 17 (33) ◽  
pp. 21583-21593 ◽  
Author(s):  
Sarantos Marinakis ◽  
Indigo Lily Dean ◽  
Jacek Kłos ◽  
François Lique
Keyword(s):  
Potential Energy ◽  
Potential Energy Surface ◽  
Ab Initio ◽  
Potential Energy Surfaces ◽  
Energy Surface ◽  
Experimental Results ◽  
Collisional Excitation ◽  
Energy Surfaces

We present a new CH(X)–He potential energy surface which is able to reproduce all the available experimental results.

Photorelaxation of imidazole and adenine via electron-driven proton transfer along H2O wires

2016 ◽  
Vol 195 ◽  
pp. 237-251 ◽  
Author(s):  
Rafał Szabla ◽  
Robert W. Góra ◽  
Mikołaj Janicki ◽  
Jiří Šponer
Keyword(s):  
Molecular Dynamics ◽  
Excited State ◽  
Potential Energy ◽  
Proton Transfer ◽  
Potential Energy Surface ◽  
Molecular Dynamics Simulations ◽  
Gas Phase ◽  
Energy Surface ◽  
Water Molecules ◽  
Dynamics Simulations

Photochemically created πσ* states were classified among the most prominent factors determining the ultrafast radiationless deactivation and photostability of many biomolecular building blocks. In the past two decades, the gas phase photochemistry of πσ* excitations was extensively investigated and was attributed to N–H and O–H bond fission processes. However, complete understanding of the complex photorelaxation pathways of πσ* states in the aqueous environment was very challenging, owing to the direct participation of solvent molecules in the excited-state deactivation. Here, we present non-adiabatic molecular dynamics simulations and potential energy surface calculations of the photoexcited imidazole–(H2O)5 cluster using the algebraic diagrammatic construction method to the second-order [ADC(2)]. We show that electron driven proton transfer (EDPT) along a wire of at least two water molecules may lead to the formation of a πσ*/S0 state crossing, similarly to what we suggested for 2-aminooxazole. We expand on our previous findings by direct comparison of the imidazole–(H2O)5 cluster to non-adiabatic molecular dynamics simulations of imidazole in the gas phase, which reveal that the presence of water molecules extends the overall excited-state lifetime of the chromophore. To embed the results in a biological context, we provide calculations of potential energy surface cuts for the analogous photorelaxation mechanism present in adenine, which contains an imidazole ring in its structure.

Sign in / Sign up

Export Citation Format

Share Document