scholarly journals Kinetically-constrained ring-polymer molecular dynamics for non-adiabatic chemistries involving solvent and donor–acceptor dynamical effects

2016 ◽  
Vol 195 ◽  
pp. 191-214 ◽  
Author(s):  
Joshua S. Kretchmer ◽  
Thomas F. Miller III

We investigate the performance of the recently developed kinetically-constrained ring polymer molecular dynamics (KC-RPMD) method for the description of model condensed-phase electron transfer (ET) reactions in which solvent and donor–acceptor dynamics play an important role. Comparison of KC-RPMD with results from Golden-Rule rate theories and numerically exact quantum dynamics calculations demonstrates that KC-RPMD accurately captures the combination of electronic- and nuclear-dynamical effects throughout the Marcus (intermediate solvent friction) and Zusman (large solvent friction) regimes of ET. It is also demonstrated that KC-RPMD accurately describes systems in which the magnitude of the diabatic coupling depends on the position of a dynamical donor–acceptor mode. In addition to these successes, however, we present an unsurprising failure of KC-RPMD to capture the enhancement of the ET rate in the low solvent friction regime associated with nuclear coherence effects. In this analysis, we re-visit several aspects of the original KC-RPMD formulation, including the form of the kinetic constraint and the choice of the mass of the auxiliary electronic variable. In particular, we introduce a Langevin bath for the auxiliary electronic variable to correct for its unphysically low coupling with the nuclear degrees of freedom. This work demonstrates that the KC-RPMD method is well suited for the direct simulation of non-adiabatic donor–acceptor chemistries associated with many complex systems, including those for which solvent dynamics plays an important role in the reaction mechanism.

2020 ◽  
Author(s):  
Sutirtha N. Chowdhury ◽  
Pengfei Huo

We present the non-adiabatic Matsubara dynamics, a general framework for computing the time-correlation function (TCF) of electronically non-adiabatic systems. This new formalism is derived based on the generalized Kubo-transformed time-correlation function, using the Wigner representation for both the nuclear degrees of freedom (DOF) and the electronic mapping variables. By dropping the non-Matsubara nuclear normal modes in the quantum Liouvillian and explicitly integrate these modes out of the TCF, we derived the non-adiabatic Matsubara dynamics approach. Further making the approximation to drop the imaginary part of the Matsubara Liouvillian and enforce the nuclear momentum integral to be real, we arrived at the non-adiabatic ring-polymer molecular dynamics (NRPMD) approach. We have further justified the capability of NRPMD for simulating the non-equilibrium time-correlation function. This work provides the rigorous theoretical foundation for several recently proposed state-dependent RPMD approaches and offers a general framework for developing new non-adiabatic quantum dynamics approaches in the future.


2020 ◽  
Author(s):  
Sutirtha N. Chowdhury ◽  
Pengfei Huo

We present the non-adiabatic Matsubara dynamics, a general framework for computing the time-correlation function (TCF) of electronically non-adiabatic systems. This new formalism is derived based on the generalized Kubo-transformed time-correlation function, using the Wigner representation for both the nuclear degrees of freedom (DOF) and the electronic mapping variables. By dropping the non-Matsubara nuclear normal modes in the quantum Liouvillian and explicitly integrate these modes out of the TCF, we derived the non-adiabatic Matsubara dynamics approach. Further making the approximation to drop the imaginary part of the Matsubara Liouvillian and enforce the nuclear momentum integral to be real, we arrived at the non-adiabatic ring-polymer molecular dynamics (NRPMD) approach. We have further justified the capability of NRPMD for simulating the non-equilibrium time-correlation function. This work provides the rigorous theoretical foundation for several recently proposed state-dependent RPMD approaches and offers a general framework for developing new non-adiabatic quantum dynamics approaches in the future.


2021 ◽  
Vol 23 (10) ◽  
pp. 6141-6153
Author(s):  
Jianwei Cao ◽  
Yanan Wu ◽  
Haitao Ma ◽  
Zhitao Shen ◽  
Wensheng Bian

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


2021 ◽  
Author(s):  
Xinyang Li ◽  
Pengfei Huo

<div>We use the ab-initio ring polymer molecular dynamics (RPMD) approach to investigate tunneling controlled reactions in methylhydroxycarbene. Nuclear tunneling effects enable molecules to overcome the barriers which can not be overcome classically. Under low-temperature conditions, intrinsic quantum tunneling effects canfacilitate the chemical reaction in a pathway that is neither favored thermodynamically nor kinetically. This</div><div>behavior is referred to as the tunneling controlled chemical reaction and regarded as the third paradigm of chemical</div><div>reaction controls. In this work, we use the ab-initio RPMD approach to incorporate the tunneling effects in our quantum dynamics simulations. The reaction kinetics of two competitive reaction pathways at various temperatures are investigated with the Kohn-Sham density functional theory (KS-DFT) on-the-fly molecular dynamics simulations and the ring polymer quantization of the nuclei. The reaction rate constants obtained here agree extremely well with the experimentally measured rates. We demonstrate the feasibility of using ab-initio RPMD rate calculations in a realistic molecular system, and provide an interesting and important example for future investigations on reaction mechanisms dominated by quantum tunneling effects.</div>


2015 ◽  
Vol 142 (19) ◽  
pp. 191101 ◽  
Author(s):  
Timothy J. H. Hele ◽  
Michael J. Willatt ◽  
Andrea Muolo ◽  
Stuart C. Althorpe

2021 ◽  
Author(s):  
Xinyang Li ◽  
Pengfei Huo

<div>We use the ab-initio ring polymer molecular dynamics (RPMD) approach to investigate tunneling controlled reactions in methylhydroxycarbene. Nuclear tunneling effects enable molecules to overcome the barriers which can not be overcome classically. Under low-temperature conditions, intrinsic quantum tunneling effects canfacilitate the chemical reaction in a pathway that is neither favored thermodynamically nor kinetically. This</div><div>behavior is referred to as the tunneling controlled chemical reaction and regarded as the third paradigm of chemical</div><div>reaction controls. In this work, we use the ab-initio RPMD approach to incorporate the tunneling effects in our quantum dynamics simulations. The reaction kinetics of two competitive reaction pathways at various temperatures are investigated with the Kohn-Sham density functional theory (KS-DFT) on-the-fly molecular dynamics simulations and the ring polymer quantization of the nuclei. The reaction rate constants obtained here agree extremely well with the experimentally measured rates. We demonstrate the feasibility of using ab-initio RPMD rate calculations in a realistic molecular system, and provide an interesting and important example for future investigations on reaction mechanisms dominated by quantum tunneling effects.</div>


2020 ◽  
Vol 22 (41) ◽  
pp. 23657-23664
Author(s):  
Yang Liu ◽  
Hongwei Song ◽  
Jun Li

The kinetics of the title reaction is studied by running the ring polymer molecular dynamics and quantum dynamics on an accurate potential energy surface.


2016 ◽  
Vol 195 ◽  
pp. 253-268 ◽  
Author(s):  
Jessica Ryan Duke ◽  
Nandini Ananth

We present a mean field ring polymer molecular dynamics method to calculate the rate of electron transfer (ET) in multi-state, multi-electron condensed-phase processes. Our approach involves calculating a transition state theory (TST) estimate to the rate using an exact path integral in discrete electronic states and continuous Cartesian nuclear coordinates. A dynamic recrossing correction to the TST rate is then obtained from real-time dynamics simulations using mean field ring polymer molecular dynamics. We employ two different reaction coordinates in our simulations and show that, despite the use of mean field dynamics, the use of an accurate dividing surface to compute TST rates allows us to achieve remarkable agreement with Fermi's golden rule rates for nonadiabatic ET in the normal regime of Marcus theory. Further, we show that using a reaction coordinate based on electronic state populations allows us to capture the turnover in rates for ET in the Marcus inverted regime.


Sign in / Sign up

Export Citation Format

Share Document