Potential energy surface stationary points and dynamics of the F−+ CH3I double inversion mechanism

2017 ◽  
Vol 19 (30) ◽  
pp. 20127-20136 ◽  
Author(s):  
Yong-Tao Ma ◽  
Xinyou Ma ◽  
Anyang Li ◽  
Hua Guo ◽  
Li Yang ◽  
...  
Keyword(s):  
Potential Energy ◽  
Potential Energy Surface ◽  
Energy Surface ◽  
Stationary Points ◽  
Direct Dynamics ◽  
Dynamics Simulations ◽  
Inversion Mechanism

Direct dynamics simulations were performed to study the SN2 double inversion mechanism SN2-DI, with retention of configuration, for the F−+ CH3I reaction.

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 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.

The search for stationary points on a quantum mechanical/molecular mechanical potential-energy surface

2002 ◽  
Vol 107 (3) ◽  
pp. 147-153 ◽  
Author(s):  
Xavier Prat-Resina ◽  
Mireia Garcia-Viloca ◽  
Gerald Monard ◽  
Angels González-Lafont ◽  
José M. Lluch ◽  
...  
Keyword(s):  
Potential Energy ◽  
Potential Energy Surface ◽  
Energy Surface ◽  
Stationary Points ◽  
Quantum Mechanical

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.

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>

DUAL-LEVEL DIRECT DYNAMICS STUDIES ON THE HYDROGEN ABSTRACTION REACTION OF CF3CHFCF3 (HFC-227ea) WITH Cl ATOM

2013 ◽  
Vol 12 (02) ◽  
pp. 1250119
Author(s):  
MASOUD ARABIEH ◽  
MANSOUR ZAHEDI
Keyword(s):  
Potential Energy ◽  
Potential Energy Surface ◽  
Energy Surface ◽  
Hydrogen Abstraction ◽  
Rate Coefficients ◽  
Temperature Range ◽  
Direct Dynamics ◽  
Abstraction Reaction ◽  
Hydrogen Abstraction Reaction ◽  
Cl Atom

The hydrogen abstraction reaction of CF3CHFCF3 (HFC-227ea) by Cl atom is studied theoretically over a temperature range of 50–2000 K. The reaction path information such as geometries, gradients, etc. are initially obtained at the BB1K/6-31+G (d,p) level. In order to improve the energies, high level single point energy (SPE) calculations are carried out at BB1K/MG3S on BB1K/6-31+G(d,p) geometries. The rate coefficients are obtained by dual-level direct dynamics approach with the interpolated SPE method at BB1K/MG3S//BB1K/6-31+G(d,p) level. Aforementioned rate coefficients are calculated based on canonical variational transition state theory (CVT) with small-curvature tunneling (SCT) correction. The results show that for the title reaction, variational effect is small while the tunneling effects are significant at low temperature. In order to investigate the reliability of improved BB1K potential energy surface in predicting barrier height of the reaction, different quantum models built such as QCISD(T)/6-311++G(d,p)//MP2/6-31G(d,p) , QCISD(T)/aug-cc-pv-DZ//MP2/6-31G(d,p) , and MP4SDQ/aug-cc-pv-DZ//MP2/6-31G(d,p) are applied to calculate the rate coefficients of the reaction. The results show the rate coefficients calculated based on MP2 potential energy surface are underestimated as compared to the experimental value. Finally three models are applied to fit the calculated rate coefficients over the temperature range of 50–2000 K. It is shown that the new four-parameter model provided a better temperature dependence than the other expressions. It is expected that obtained theoretical results can be useful for modeling the kinetics nature of the reaction of HFC-227ea and Cl atom over temperature range where no experimental data is accessible.

Sign in / Sign up

Export Citation Format

Share Document