1
|
Li S, Xie BB, Yin BW, Liu L, Shen L, Fang WH. Construction of Highly Accurate Machine Learning Potential Energy Surfaces for Excited-State Dynamics Simulations Based on Low-Level Data Sets. J Phys Chem A 2024; 128:5516-5524. [PMID: 38954640 DOI: 10.1021/acs.jpca.4c02028] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/04/2024]
Abstract
Machine learning is capable of effectively predicting the potential energies of molecules in the presence of high-quality data sets. Its application in the construction of ground- and excited-state potential energy surfaces is attractive to accelerate nonadiabatic molecular dynamics simulations of photochemical reactions. Because of the huge computational cost of excited-state electronic structure calculations, the construction of a high-quality data set becomes a bottleneck. In the present work, we first built two data sets. One was obtained from surface hopping dynamics simulations at the semiempirical OM2/MRCI level. Another was extracted from the dynamics trajectories at the CASSCF level, which was reported previously. The ground- and excited-state potential energy surfaces of ethylene-bridged azobenzene at the CASSCF computational level were constructed based on the former low-level data set. Although non-neural network machine learning methods can achieve good or modest performance during the training process, only neural network models provide reliable predictions on the latter external test data set. The BPNN and SchNet combined with the Δ-ML scheme and the force term in the loss functions are recommended for dynamics simulations. Then, we performed excited-state dynamics simulations of the photoisomerization of ethylene-bridged azobenzene on machine learning potential energy surfaces. Compared with the lifetimes of the first excited state (S1) estimated at different computational levels, our results on the E isomer are in good agreement with the high-level estimation. However, the overestimation of the Z isomer is unimproved. It suggests that smaller errors during the training process do not necessarily translate to more accurate predictions on high-level potential energies or better performance on nonadiabatic dynamics simulations, at least in the present case.
Collapse
Affiliation(s)
- Shuai Li
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Bin-Bin Xie
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, Hangzhou 311231, Zhejiang, P. R. China
| | - Bo-Wen Yin
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, Hangzhou 311231, Zhejiang, P. R. China
| | - Lihong Liu
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Lin Shen
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
- Yantai-Jingshi Institute of Material Genome Engineering, Yantai 265505, Shandong, P. R. China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
- Shandong Laboratory of Yantai Advanced Materials and Green Manufacturing, Yantai 264006, Shandong, P. R. China
| |
Collapse
|
2
|
Dupuy L, Rikus A, Maitra NT. Exact-Factorization-Based Surface Hopping without Velocity Adjustment. J Phys Chem Lett 2024:2643-2649. [PMID: 38422391 DOI: 10.1021/acs.jpclett.4c00115] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/02/2024]
Abstract
While surface hopping has emerged as a powerful method for simulating non-adiabatic dynamics in large molecules, the ad hoc nature of the necessary velocity adjustments and decoherence corrections in the algorithm somewhat reduces its reliability. Here we propose a new scheme that eliminates these aspects by combining the nuclear equation from the quantum-trajectory surface-hopping approach with the electronic equation derived from the exact-factorization approach. The resulting method, denoted QTSH-XF, yields a surface-hopping method on firmer ground than previous and is shown to successfully capture dynamics in Tully models and in a linear vibronic coupling model of the photoexcited uracil cation.
Collapse
Affiliation(s)
- Lucien Dupuy
- Department of Physics, Rutgers University, Newark, New Jersey 07102, United States
| | - Anton Rikus
- Department of Physics, Rutgers University, Newark, New Jersey 07102, United States
- University of Münster, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, 48149 Münster, Germany
| | - Neepa T Maitra
- Department of Physics, Rutgers University, Newark, New Jersey 07102, United States
| |
Collapse
|
3
|
Arribas EV, Ibele LM, Lauvergnat D, Maitra NT, Agostini F. Significance of Energy Conservation in Coupled-Trajectory Approaches to Nonadiabatic Dynamics. J Chem Theory Comput 2023; 19:7787-7800. [PMID: 37853509 DOI: 10.1021/acs.jctc.3c00845] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2023]
Abstract
Through approximating electron-nuclear correlation terms in the exact factorization approach, trajectory-based methods have been derived and successfully applied to the dynamics of a variety of light-induced molecular processes, capturing quantum (de)coherence effects rigorously. These terms account for the coupling among the trajectories, recovering the nonlocal nature of quantum nuclear dynamics that is completely overlooked in traditional independent-trajectory algorithms. Nevertheless, some of the approximations introduced in the derivation of some of these methods do not conserve the total energy. We analyze energy conservation in the coupled-trajectory mixed quantum-classical (CTMQC) algorithm and explore the performance of a modified algorithm, CTMQC-E, where some of the terms are redefined to restore energy conservation. A set of molecular models is used as a test, namely, 2-cis-penta-2,4-dienimium cation, bis(methylene) adamantyl radical cation, butatriene cation, uracil radical cation, and neutral pyrazine.
Collapse
Affiliation(s)
| | - Lea M Ibele
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| | - Neepa T Maitra
- Department of Physics, Rutgers University, Newark, New Jersey 07102, United States
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| |
Collapse
|
4
|
Villaseco Arribas E, Vindel-Zandbergen P, Roy S, Maitra NT. Different flavors of exact-factorization-based mixed quantum-classical methods for multistate dynamics. Phys Chem Chem Phys 2023; 25:26380-26395. [PMID: 37750820 DOI: 10.1039/d3cp03464j] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/27/2023]
Abstract
The exact factorization approach has led to the development of new mixed quantum-classical methods for simulating coupled electron-ion dynamics. We compare their performance for dynamics when more than two electronic states are occupied at a given time, and analyze: (1) the use of coupled versus auxiliary trajectories in evaluating the electron-nuclear correlation terms, (2) the approximation of using these terms within surface-hopping and Ehrenfest frameworks, and (3) the relevance of the exact conditions of zero population transfer away from nonadiabatic coupling regions and total energy conservation. Dynamics through the three-state conical intersection in the uracil radical cation as well as polaritonic models in one dimension are studied.
Collapse
Affiliation(s)
| | - Patricia Vindel-Zandbergen
- Department of Physics, Rutgers University, Newark 07102, New Jersey, USA.
- Department of Chemistry, New York University, New York, New York 10003, USA
| | - Saswata Roy
- Department of Physics, Rutgers University, Newark 07102, New Jersey, USA.
| | - Neepa T Maitra
- Department of Physics, Rutgers University, Newark 07102, New Jersey, USA.
| |
Collapse
|
5
|
Tang D, Jia L, Shen L, Fang WH. Fewest-Switches Surface Hopping with Long Short-Term Memory Networks. J Phys Chem Lett 2022; 13:10377-10387. [PMID: 36317657 DOI: 10.1021/acs.jpclett.2c02299] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
The mixed quantum-classical dynamical simulation is essential for studying nonadiabatic phenomena in photophysics and photochemistry. In recent years, many machine learning models have been developed to accelerate the time evolution of the nuclear subsystem. Herein, we implement long short-term memory (LSTM) networks as a propagator to accelerate the time evolution of the electronic subsystem during the fewest-switches surface hopping (FSSH) simulations. A small number of reference trajectories are generated using the original FSSH method, and then the LSTM networks can be built, accompanied by careful examination of typical LSTM-FSSH trajectories that employ the same initial condition and random numbers as the corresponding reference. The constructed network is applied to FSSH to further produce a trajectory ensemble to reveal the mechanism of nonadiabatic processes. Taking Tully's three models as test systems, we qualitatively reproduced the collective results. This work demonstrates that LSTM can be applied to the most popular surface hopping simulations.
Collapse
Affiliation(s)
- Diandong Tang
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Luyang Jia
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Lin Shen
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
- Yantai-Jingshi Institute of Material Genome Engineering, Yantai 265505, Shandong, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
- Yantai-Jingshi Institute of Material Genome Engineering, Yantai 265505, Shandong, China
| |
Collapse
|
6
|
Villaseco Arribas E, Agostini F, Maitra NT. Exact Factorization Adventures: A Promising Approach for Non-Bound States. Molecules 2022; 27:molecules27134002. [PMID: 35807246 PMCID: PMC9267945 DOI: 10.3390/molecules27134002] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2022] [Revised: 06/07/2022] [Accepted: 06/09/2022] [Indexed: 11/29/2022] Open
Abstract
Modeling the dynamics of non-bound states in molecules requires an accurate description of how electronic motion affects nuclear motion and vice-versa. The exact factorization (XF) approach offers a unique perspective, in that it provides potentials that act on the nuclear subsystem or electronic subsystem, which contain the effects of the coupling to the other subsystem in an exact way. We briefly review the various applications of the XF idea in different realms, and how features of these potentials aid in the interpretation of two different laser-driven dissociation mechanisms. We present a detailed study of the different ways the coupling terms in recently-developed XF-based mixed quantum-classical approximations are evaluated, where either truly coupled trajectories, or auxiliary trajectories that mimic the coupling are used, and discuss their effect in both a surface-hopping framework as well as the rigorously-derived coupled-trajectory mixed quantum-classical approach.
Collapse
Affiliation(s)
| | - Federica Agostini
- Institut de Chimie Physique UMR8000, Université Paris-Saclay, CNRS, 91405 Orsay, France;
| | - Neepa T. Maitra
- Department of Physics, Rutgers University, Newark, NJ 07102, USA;
- Correspondence:
| |
Collapse
|
7
|
Vindel-Zandbergen P, Ibele LM, Ha JK, Min SK, Curchod BFE, Maitra NT. Study of the Decoherence Correction Derived from the Exact Factorization Approach for Nonadiabatic Dynamics. J Chem Theory Comput 2021; 17:3852-3862. [PMID: 34138553 PMCID: PMC8280698 DOI: 10.1021/acs.jctc.1c00346] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
![]()
We present a detailed
study of the decoherence correction to surface
hopping that was recently derived from the exact factorization approach.
Ab initio multiple spawning calculations that use the same initial
conditions and the same electronic structure method are used as a
reference for three molecules: ethylene, the methaniminium cation,
and fulvene, for which nonadiabatic dynamics follows a photoexcitation.
A comparison with the Granucci–Persico energy-based decoherence
correction and the augmented fewest-switches surface-hopping scheme
shows that the three decoherence-corrected methods operate on individual
trajectories in a qualitatively different way, but the results averaged
over trajectories are similar for these systems.
Collapse
Affiliation(s)
| | - Lea M Ibele
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, U.K
| | - Jong-Kwon Ha
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulsan 44919, Republic of Korea
| | - Seung Kyu Min
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulsan 44919, Republic of Korea
| | - Basile F E Curchod
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, U.K
| | - Neepa T Maitra
- Department of Physics, Rutgers University, Newark, New Jersey 07102, United States
| |
Collapse
|