1
|
Cook L, Runeson JE, Richardson JO, Hele TJH. Which Algorithm Best Propagates the Meyer-Miller-Stock-Thoss Mapping Hamiltonian for Non-Adiabatic Dynamics? J Chem Theory Comput 2023; 19:6109-6125. [PMID: 37704193 PMCID: PMC10536990 DOI: 10.1021/acs.jctc.3c00709] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2023] [Indexed: 09/15/2023]
Abstract
A common strategy to simulate mixed quantum-classical dynamics is by propagating classical trajectories with mapping variables, often using the Meyer-Miller-Stock-Thoss (MMST) Hamiltonian or the related spin-mapping approach. When mapping the quantum subsystem, the coupled dynamics reduce to a set of equations of motion to integrate. Several numerical algorithms have been proposed, but a thorough performance comparison appears to be lacking. Here, we compare three time-propagation algorithms for the MMST Hamiltonian: the Momentum Integral (MInt) (J. Chem. Phys., 2018, 148, 102326), the Split-Liouvillian (SL) (Chem. Phys., 2017, 482, 124-134), and the algorithm in J. Chem. Phys., 2012, 136, 084101 that we refer to as the Degenerate Eigenvalue (DE) algorithm due to the approximation required during derivation. We analyze the accuracy of individual trajectories, correlation functions, energy conservation, symplecticity, Liouville's theorem, and the computational cost. We find that the MInt algorithm is the only rigorously symplectic algorithm. However, comparable accuracy at a lower computational cost can be obtained with the SL algorithm. The approximation implicitly made within the DE algorithm conserves energy poorly, even for small timesteps, and thus leads to slightly different results. These results should guide future mapping-variable simulations.
Collapse
Affiliation(s)
- Lauren
E. Cook
- Department
of Chemistry, University College London, Christopher Ingold Building, London WC1H 0AJ, U.K.
| | - Johan E. Runeson
- Department
of Chemistry and Applied Biosciences, ETH
Zürich, Zürich 8093, Switzerland
| | - Jeremy O. Richardson
- Department
of Chemistry and Applied Biosciences, ETH
Zürich, Zürich 8093, Switzerland
| | - Timothy J. H. Hele
- Department
of Chemistry, University College London, Christopher Ingold Building, London WC1H 0AJ, U.K.
| |
Collapse
|
2
|
Malpathak S, Church MS, Ananth N. A Semiclassical Framework for Mixed Quantum Classical Dynamics. J Phys Chem A 2022; 126:6359-6375. [PMID: 36070472 DOI: 10.1021/acs.jpca.2c03467] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Semiclassical (SC) approximations for quantum dynamic simulations in complex chemical systems range from rigorously accurate methods that are computationally expensive to methods that exhibit near-classical scaling with system size but are limited in their ability to describe quantum effects. In practical studies of high-dimensional reactions, neither extreme is the best choice: frequently a high-level quantum mechanical description is only required for a handful of modes, while the majority of environment modes that do not play a key role in the reactive event of interest are well served with a lower level of theory. In this feature, we introduce modified Filinov filtration as a powerful tool to construct mixed quantum-classical SC theories where different subsystems can be quantized to different extents without introducing ad hoc intersubsystem interaction terms. We demonstrate that these Filinov-based SC methods can systematically tune between quantum and classical limit SC behavior, offering a practical way forward to accurate and computationally efficient simulations of high-dimensional quantum processes.
Collapse
Affiliation(s)
- Shreyas Malpathak
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University Ithaca, New York 14853, United States
| | - Matthew S Church
- Department of Chemistry, Brown University, Providence, Rhode Island 02906, United States
| | - Nandini Ananth
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University Ithaca, New York 14853, United States
| |
Collapse
|
3
|
Church MS, Ananth N. Semiclassical dynamics in the mixed quantum-classical limit. J Chem Phys 2019; 151:134109. [PMID: 31594341 DOI: 10.1063/1.5117160] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The semiclassical double Herman-Kluk initial value representation is an accurate approach to computing quantum real time correlation functions, but its applications are limited by the need to evaluate an oscillatory integral. In previous work, we have shown that this "sign problem" can be mitigated using the modified Filinov filtration technique to control the extent to which individual modes of the system contribute to the overall phase of the integrand. Here, we follow this idea to a logical conclusion: we analytically derive a general expression for the mixed quantum-classical limit of the semiclassical correlation function-analytical mixed quantum-classical-initial value representation (AMQC-IVR), where the phase contributions from the "classical" modes of the system are filtered while the "quantum" modes are treated in the full semiclassical limit. We numerically demonstrate the accuracy and efficiency of the AMQC-IVR formulation in calculations of quantum correlation functions and reaction rates using three model systems with varied coupling strengths between the classical and quantum subsystems. We also introduce a separable prefactor approximation that further reduces computational cost but is only accurate in the limit of weak coupling between the quantum and classical subsystems.
Collapse
Affiliation(s)
- Matthew S Church
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Nandini Ananth
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
4
|
Cendagorta JR, Bačić Z, Tuckerman ME. An open-chain imaginary-time path-integral sampling approach to the calculation of approximate symmetrized quantum time correlation functions. J Chem Phys 2018; 148:102340. [PMID: 29544313 DOI: 10.1063/1.5005543] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023] Open
Abstract
We introduce a scheme for approximating quantum time correlation functions numerically within the Feynman path integral formulation. Starting with the symmetrized version of the correlation function expressed as a discretized path integral, we introduce a change of integration variables often used in the derivation of trajectory-based semiclassical methods. In particular, we transform to sum and difference variables between forward and backward complex-time propagation paths. Once the transformation is performed, the potential energy is expanded in powers of the difference variables, which allows us to perform the integrals over these variables analytically. The manner in which this procedure is carried out results in an open-chain path integral (in the remaining sum variables) with a modified potential that is evaluated using imaginary-time path-integral sampling rather than requiring the generation of a large ensemble of trajectories. Consequently, any number of path integral sampling schemes can be employed to compute the remaining path integral, including Monte Carlo, path-integral molecular dynamics, or enhanced path-integral molecular dynamics. We believe that this approach constitutes a different perspective in semiclassical-type approximations to quantum time correlation functions. Importantly, we argue that our approximation can be systematically improved within a cumulant expansion formalism. We test this approximation on a set of one-dimensional problems that are commonly used to benchmark approximate quantum dynamical schemes. We show that the method is at least as accurate as the popular ring-polymer molecular dynamics technique and linearized semiclassical initial value representation for correlation functions of linear operators in most of these examples and improves the accuracy of correlation functions of nonlinear operators.
Collapse
Affiliation(s)
| | - Zlatko Bačić
- Department of Chemistry, New York University, New York, New York 10003, USA
| | - Mark E Tuckerman
- Department of Chemistry, New York University, New York, New York 10003, USA
| |
Collapse
|
5
|
Aieta C, Ceotto M. A quantum method for thermal rate constant calculations from stationary phase approximation of the thermal flux-flux correlation function integral. J Chem Phys 2017; 146:214115. [DOI: 10.1063/1.4984099] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Chiara Aieta
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133 Milano, Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133 Milano, Italy
| |
Collapse
|
6
|
Petersen J, Kay KG. Wave packet propagation across barriers by semiclassical initial value methods. J Chem Phys 2015; 143:014107. [PMID: 26156465 DOI: 10.1063/1.4923221] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Semiclassical initial value representation (IVR) formulas for the propagator have difficulty describing tunneling through barriers. A key reason is that these formulas do not automatically reduce, in the classical limit, to the version of the Van Vleck-Gutzwiller (VVG) propagator required to treat barrier tunneling, which involves trajectories that have complex initial conditions and that follow paths in complex time. In this work, a simple IVR expression, that has the correct tunneling form in the classical limit, is derived for the propagator in the case of one-dimensional barrier transmission. Similarly, an IVR formula, that reduces to the Generalized Gaussian Wave Packet Dynamics (GGWPD) expression [D. Huber, E. J. Heller, and R. Littlejohn, J. Chem. Phys. 89, 2003 (1988)] in the classical limit, is derived for the transmitted wave packet. Uniform semiclassical versions of the IVR formulas are presented and simplified expressions in terms of real trajectories and WKB penetration factors are described. Numerical tests show that the uniform IVR treatment gives good results for wave packet transmission through the Eckart and Gaussian barriers in all cases examined. In contrast, even when applied with the proper complex trajectories, the VVG and GGWPD treatments are inaccurate when the mean energy of the wave packet is near the classical transmission threshold. The IVR expressions for the propagator and wave packet are cast as contour integrals in the complex space of initial conditions and these are generalized to potentially allow treatment of a larger variety of systems. A steepest descent analysis of the contour integral formula for the wave packet in the present cases confirms its relationship to the GGWPD method, verifies its semiclassical validity, and explains results of numerical calculations.
Collapse
Affiliation(s)
- Jakob Petersen
- Department of Chemical Physics, Weizmann Institute of Science, Rehovot 76100, Israel
| | - Kenneth G Kay
- Department of Chemistry, Bar-Ilan University, Ramat-Gan 52900, Israel
| |
Collapse
|
7
|
Kryvohuz M. Calculation of chemical reaction rate constants using on-the-fly high level electronic structure computations with account of multidimensional tunneling. J Chem Phys 2012; 137:234304. [DOI: 10.1063/1.4769195] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
8
|
Tao G, Miller WH. Time-dependent importance sampling in semiclassical initial value representation calculations for time correlation functions. II. A simplified implementation. J Chem Phys 2012; 137:124105. [DOI: 10.1063/1.4752206] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
9
|
Mazzuca J, Garashchuk S, Jakowski J. Description of proton transfer in soybean lipoxygenase-1 employing approximate quantum trajectory dynamics. Chem Phys Lett 2012. [DOI: 10.1016/j.cplett.2012.06.019] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
10
|
Head-Gordon M. Publications. Mol Phys 2012. [DOI: 10.1080/00268976.2012.693252] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
11
|
Tao G, Miller WH. Time-dependent importance sampling in semiclassical initial value representation calculations for time correlation functions. J Chem Phys 2011; 135:024104. [DOI: 10.1063/1.3600656] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
12
|
Goussev A, Schubert R, Waalkens H, Wiggins S. The flux-flux correlation function for anharmonic barriers. J Chem Phys 2010; 133:244113. [DOI: 10.1063/1.3518425] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
|
13
|
|
14
|
Garashchuk S, Vazhappilly T. Wavepacket approach to the cumulative reaction probability within the flux operator formalism. J Chem Phys 2009; 131:164108. [PMID: 19894928 DOI: 10.1063/1.3251333] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Expressions for the singular flux operator eigenfunctions and eigenvalues are given in terms of the Dirac delta-function representable as a localized Gaussian wavepacket. This functional form enables computation of the cumulative reaction probability N(E) from the wavepacket time-correlation functions. The Gaussian based form of the flux eigenfunctions, which is not tied to a finite basis of a quantum-mechanical calculation, is particularly useful for approximate calculation of N(E) with the trajectory based wavepacket propagation techniques. Numerical illustration is given for the Eckart barrier using the conventional quantum-mechanical propagation and the quantum trajectory dynamics with the approximate quantum potential. N(E) converges with respect to the Gaussian width parameter, and the convergence is faster at low energy. The approximate trajectory calculation overestimates tunneling in the low energy regime, but gives a significant improvement over the parabolic estimate of the tunneling probability.
Collapse
Affiliation(s)
- Sophya Garashchuk
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, South Carolina 29208, USA.
| | | |
Collapse
|
15
|
Dawes R, Passalacqua A, Wagner AF, Sewell TD, Minkoff M, Thompson DL. Interpolating moving least-squares methods for fitting potential energy surfaces: Using classical trajectories to explore configuration space. J Chem Phys 2009; 130:144107. [DOI: 10.1063/1.3111261] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
16
|
Sibert EL, Ramesh SG, Gulmen TS. Vibrational Relaxation of OH and CH Fundamentals of Polar and Nonpolar Molecules in the Condensed Phase. J Phys Chem A 2008; 112:11291-305. [DOI: 10.1021/jp8068442] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Edwin L. Sibert
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, Madison, Wisconsin 53706
| | - Sai G. Ramesh
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, Madison, Wisconsin 53706
| | - Tolga S. Gulmen
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, Madison, Wisconsin 53706
| |
Collapse
|
17
|
Issack BB, Roy PN. Semiclassical initial value representation treatment of a hydrogen bonded complex of rigid water molecules from a single trajectory in Cartesian coordinates. J Chem Phys 2007; 127:144306. [DOI: 10.1063/1.2786456] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
|