1
|
Li G, Shi Z, Huang L, Wang L. Multiconfigurational Surface Hopping: a Time-Dependent Variational Approach with Momentum-Jump Trajectories. J Chem Theory Comput 2024. [PMID: 39215702 DOI: 10.1021/acs.jctc.4c00842] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
The Ehrenfest mean field dynamics and trajectory surface hopping have been widely used in nonadiabatic dynamics simulations. Based on the time-dependent variational principle (TDVP), the multiconfigurational Ehrenfest (MCE) method has also been developed and can be regarded as a multiconfigurational extension of the traditional Ehrenfest dynamics. However, it is not straightforward to apply the TDVP to surface hopping trajectories because there exists momentum jump during surface hops. To solve this problem, we here propose a multiconfigurational surface hopping (MCSH) method, where continuous momenta are obtained by linear interpolation and the interpolated trajectories are used to construct the basis functions for TDVP in a postprocessing manner. As demonstrated in a series of representative spin-boson models, MCSH achieves high accuracy with only several hundred trajectory bases and can uniformly improve the performance of surface hopping. In principle, MCSH can be combined with all kinds of mixed quantum-classical trajectories and thus has the potential to properly describe general nonadiabatic dynamics.
Collapse
Affiliation(s)
- Guijie Li
- Zhejiang Key Laboratory of Excited-State Energy Conversion and Energy Storage, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Zhecun Shi
- Zhejiang Key Laboratory of Excited-State Energy Conversion and Energy Storage, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Lei Huang
- Zhejiang Key Laboratory of Excited-State Energy Conversion and Energy Storage, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Linjun Wang
- Zhejiang Key Laboratory of Excited-State Energy Conversion and Energy Storage, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| |
Collapse
|
2
|
Guo X, Li G, Shi Z, Wang L. Surface Hopping with Reliable Wave Function by Introducing Auxiliary Wave Packets to Trajectory Branching. J Phys Chem Lett 2024:3345-3353. [PMID: 38498301 DOI: 10.1021/acs.jpclett.4c00437] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/20/2024]
Abstract
It is well-known that the widely utilized fewest switches surface hopping method suffers from the severe overcoherence problem, and thus adiabatic populations calculated by wave functions are generally inferior to those based on active states. More importantly, to achieve a complete description of nonadiabatic dynamics, the density matrix is essential. In this paper, we present an auxiliary branching corrected surface hopping (A-BCSH) method that introduces auxiliary wave packets (WPs) on the adiabatic potential energy surfaces for trajectory branching. Both rapid and gradual separation of WP components on different surfaces are characterized, and thus the correct decoherence time along each trajectory is captured. As demonstrated in the three standard Tully models, A-BCSH exhibits excellent internal consistency. Namely, close adiabatic populations are obtained based on both wave functions and active states. In particular, A-BCSH successfully obtains a reliable time-dependent spatial distribution of the density matrix, which relies only on electronic wave functions. Due to its high performance, our A-BCSH method provides a new and highly promising perspective on further development of more consistent surface hopping with reliable wave function.
Collapse
Affiliation(s)
- Xin Guo
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Guijie Li
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Zhecun Shi
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| |
Collapse
|
3
|
Akimov AV. Energy-Conserving and Thermally Corrected Neglect of Back-Reaction Approximation Method for Nonadiabatic Molecular Dynamics. J Phys Chem Lett 2023; 14:11673-11683. [PMID: 38109379 DOI: 10.1021/acs.jpclett.3c03029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2023]
Abstract
In this work, the energy-conserving and thermally corrected neglect of the back-reaction approximation approach for nonadiabatic molecular dynamics in extended atomistic systems is developed. The new approach introduces three key corrections to the original method: (1) it enforces the total energy conservation, (2) it introduces an explicit coupling of the system to its environment, and (3) it introduces a renormalization of nonadiabatic couplings to account for a difference between the instantaneous nuclear kinetic energy and the kinetic energy of guiding trajectories. In the new approach, an auxiliary kinetic energy variable is introduced as an independent dynamical variable. The new approach produces nonzero equilibrium populations, whereas the original neglect of the back-reaction approximation method does not. It yields population relaxation time scales that are favorably comparable to the reference values, and it introduces an explicit and controllable way of dissipating energy into a bath without an assumption of the bath being at equilibrium.
Collapse
Affiliation(s)
- Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260 United States
| |
Collapse
|
4
|
Li G, Shi Z, Guo X, Wang L. What is Missing in the Mean Field Description of Spatial Distribution of Population? Important Role of Auxiliary Wave Packets in Trajectory Branching. J Phys Chem Lett 2023; 14:9855-9863. [PMID: 37890155 DOI: 10.1021/acs.jpclett.3c02690] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/29/2023]
Abstract
When the traditional Ehrenfest mean field approach is employed to simulate nonadiabatic dynamics, an effective wave packet (WP) on the average potential energy surface (PES) is utilized to describe the nuclear motion. In the fully quantum picture, however, the WP components on different adiabatic PESs gradually separate in space because they evolve under different velocities and forces. Due to trajectory branching of the WP components, proper decoherence needs to be taken into account, and the spatial distribution of population cannot be described by a single effective WP. Here, we propose an auxiliary branching corrected mean field (A-BCMF) method, where trajectories of auxiliary WPs on adiabatic PESs are introduced. As benchmarked in the three standard Tully models, A-BCMF not only gives correct channel populations but also captures an accurate time-dependent spatial distribution of population. Thereby, we reveal the important role of auxiliary WPs in solving intrinsic problems of the widely used mean field description of nonadiabatic dynamics.
Collapse
Affiliation(s)
- Guijie Li
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Zhecun Shi
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Xin Guo
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| |
Collapse
|
5
|
Shao C, Shi Z, Xu J, Wang L. Learning Decoherence Time Formulas for Surface Hopping from Quantum Dynamics. J Phys Chem Lett 2023; 14:7680-7689. [PMID: 37606199 DOI: 10.1021/acs.jpclett.3c02019] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2023]
Abstract
Surface hopping simulations have achieved great success in many different fields, but their reliability has long been limited by the overcoherence problem. We here present a general machine learning assisted approach to identify optimal decoherence time formulas for surface hopping using exact quantum dynamics as references. In order to avoid computationally expensive force calculations, we use the nuclear kinetic energy and the adiabatic energy difference to iteratively generate the descriptor space. Through multilayer screening of the candidate descriptors and discrete optimization of the relevant parameters, we obtain new energy-based decoherence time formulas. As benchmarked in thousands of diverse multilevel systems and six standard scattering models, surface hopping with our new decoherence time formulas nearly reproduces the exact quantum dynamics while maintaining high efficiency. Thereby, our approach provides a promising avenue for systematically improving the accuracy of surface hopping simulations in complex systems from quantum dynamics data.
Collapse
Affiliation(s)
- Cancan Shao
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Zhecun Shi
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Jiabo Xu
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| |
Collapse
|
6
|
Gardner J, Habershon S, Maurer RJ. Assessing Mixed Quantum-Classical Molecular Dynamics Methods for Nonadiabatic Dynamics of Molecules on Metal Surfaces. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2023; 127:15257-15270. [PMID: 37583439 PMCID: PMC10424245 DOI: 10.1021/acs.jpcc.3c03591] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/27/2023] [Revised: 06/30/2023] [Indexed: 08/17/2023]
Abstract
Mixed quantum-classical (MQC) methods for simulating the dynamics of molecules at metal surfaces have the potential to accurately and efficiently provide mechanistic insight into reactive processes. Here, we introduce simple two-dimensional models for the scattering of diatomic molecules at metal surfaces based on recently published electronic structure data. We apply several MQC methods to investigate their ability to capture how nonadiabatic effects influence molecule-metal energy transfer during the scattering process. Specifically, we compare molecular dynamics with electronic friction, Ehrenfest dynamics, independent electron surface hopping, and the broadened classical master equation approach. In the case of independent electron surface hopping, we implement a simple decoherence correction approach and assess its impact on vibrationally inelastic scattering. Our results show that simple, low-dimensional models can be used to qualitatively capture experimentally observed vibrational energy transfer and provide insight into the relative performance of different MQC schemes. We observe that all approaches predict similar kinetic energy dependence but return different vibrational energy distributions. Finally, by varying the molecule-metal coupling, we can assess the coupling regime in which some MQC methods become unsuitable.
Collapse
Affiliation(s)
- James Gardner
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Scott Habershon
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Reinhard J. Maurer
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
- Department
of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
7
|
Li B, Xu J, Li G, Shi Z, Wang L. A Mixed Deterministic-Stochastic Algorithm of the Branching Corrected Mean Field Method for Nonadiabatic Dynamics. J Chem Phys 2022; 156:114116. [DOI: 10.1063/5.0084013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a new algorithm of the branching corrected mean field (BCMF) method for nonadiabatic dynamics [J. Xu and L. Wang, J. Phys. Chem. Lett. 11, 8283 (2020)], which combines the key advantages of the two existed algorithms, i.e., the deterministic BCMF algorithm based on weights of trajectory branches (BCMF-w) and the stochastic BCMF algorithm with random collapse of the electronic wavefunction (BCMF-s). The resulting mixed deterministic-stochastic BCMF algorithm (BCMF-ws) is benchmarked in a series of standard scattering problems with potential wells on the excited-state surfaces, which are common in realistic systems. In all investigated cases, BCMF-ws holds the same high accuracy while the computational time is reduced about two orders of magnitude compared to the original BCMF-w and BCMF-s algorithms, thus promising for nonadiabatic dynamics simulations of general systems.
Collapse
Affiliation(s)
| | | | | | | | - Linjun Wang
- Department of Chemistry, Zhejiang University, China
| |
Collapse
|
8
|
Interpretation of Adiabatic and Diabatic Populations from Trajectories of Branching Corrected Surface Hopping. CHINESE J CHEM PHYS 2022. [DOI: 10.1063/1674-0068/cjcp2201023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
|
9
|
Shao C, Xu J, Wang L. Branching and phase corrected surface hopping: A benchmark of nonadiabatic dynamics in multilevel systems. J Chem Phys 2021; 154:234109. [PMID: 34241240 DOI: 10.1063/5.0056224] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Since the seminal work of Tully [J. Chem. Phys. 93, 1061 (1990)], two-level scattering models have been extensively adopted as the standard benchmark systems to assess the performance of different trajectory surface hopping methods for nonadiabatic dynamics simulations. Here, we extend the branching and phase corrections to multilevel systems and combine them with both the traditional fewest switches surface hopping (FSSH) and its variant global flux surface hopping (GFSH) algorithms. To get a comprehensive evaluation of the proposed methods, we construct a series of more challenging and diverse three-level and four-level scattering models and use exact quantum solutions as references. Encouragingly, both FSSH and GFSH with the branching and phase corrections produce excellent and nearly identical results in all investigated systems, indicating that the new surface hopping methods are robust to describe multilevel problems and the reliability is insensitive to the definition of self-consistent hopping probabilities in the adiabatic representation. Furthermore, the branching correction is found to be especially important when dealing with strongly repulsive potential energy surfaces, which are common in realistic systems, thus promising for general applications.
Collapse
Affiliation(s)
- Cancan Shao
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Jiabo Xu
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|