1
|
Zhu Y, Peng J, Xu C, Lan Z. Unsupervised Machine Learning in the Analysis of Nonadiabatic Molecular Dynamics Simulation. J Phys Chem Lett 2024; 15:9601-9619. [PMID: 39270134 DOI: 10.1021/acs.jpclett.4c01751] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/15/2024]
Abstract
The all-atomic full-dimensional-level simulations of nonadiabatic molecular dynamics (NAMD) in large realistic systems has received high research interest in recent years. However, such NAMD simulations normally generate an enormous amount of time-dependent high-dimensional data, leading to a significant challenge in result analyses. Based on unsupervised machine learning (ML) methods, considerable efforts were devoted to developing novel and easy-to-use analysis tools for the identification of photoinduced reaction channels and the comprehensive understanding of complicated molecular motions in NAMD simulations. Here, we tried to survey recent advances in this field, particularly to focus on how to use unsupervised ML methods to analyze the trajectory-based NAMD simulation results. Our purpose is to offer a comprehensive discussion on several essential components of this analysis protocol, including the selection of ML methods, the construction of molecular descriptors, the establishment of analytical frameworks, their advantages and limitations, and persistent challenges.
Collapse
Affiliation(s)
- Yifei Zhu
- MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Jiawei Peng
- MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Chao Xu
- MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Zhenggang Lan
- MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| |
Collapse
|
2
|
Duston T, Tao Z, Bian X, Bhati M, Rawlinson J, Littlejohn RG, Pei Z, Shao Y, Subotnik JE. A Phase-Space Electronic Hamiltonian For Vibrational Circular Dichroism. J Chem Theory Comput 2024. [PMID: 39226223 DOI: 10.1021/acs.jctc.4c00662] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/05/2024]
Abstract
We show empirically that a phase-space non-Born-Oppenheimer electronic Hamiltonian approach to quantum chemistry (where the electronic Hamiltonian is parametrized by both nuclear position and momentum, ĤPS(R,P)) is both a practical and accurate means to recover vibrational circular dichroism spectra. We further hypothesize that such a phase-space approach may lead to very new dynamical physics beyond spectroscopic circular dichroism, with potential implications for understanding chiral induced spin selectivity (CISS), noting that classical phase-space approaches conserve the total nuclear plus electronic momentum, whereas classical Born-Oppenheimer approaches do not (they conserve only the nuclear momentum).
Collapse
Affiliation(s)
- Titouan Duston
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Zhen Tao
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Xuezhi Bian
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Mansi Bhati
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Jonathan Rawlinson
- Department of Mathematics, Nottingham Trent University, Nottingham NG1 4FQ, United Kingdom
| | - Robert G Littlejohn
- Department of Physics, University of California, Berkeley, California 94720, United States
| | - Zheng Pei
- Department of Chemistry, The University of Oklahoma, Norman, Oklahoma 73104, United States
| | - Yihan Shao
- Department of Chemistry, The University of Oklahoma, Norman, Oklahoma 73104, United States
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
3
|
London N, Bu S, Johnson B, Ananth N. Mean-Field Ring Polymer Rates Using a Population Dividing Surface. J Phys Chem A 2024; 128:5730-5739. [PMID: 38976564 DOI: 10.1021/acs.jpca.4c00005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/10/2024]
Abstract
Mean-field ring polymer molecular dynamics offers a computationally efficient method for the simulation of reaction rates in multilevel systems. Previous work has established that, to model a nonadiabatic state-to-state reaction accurately, it is necessary to ensure reactive trajectories form kinked ring polymer configurations at the dividing surface. Building on this idea, we introduce a population difference coordinate and a reactive flux expression modified to only include contributions from kinked configurations. We test the accuracy of the resulting mean-field rate theory on a series of linear vibronic coupling model systems. We demonstrate that this new kMF-RP rate approach is efficient to implement and quantitatively accurate for models over a wide range of driving forces, coupling strengths, and temperatures.
Collapse
Affiliation(s)
- Nathan London
- Department of Chemistry, Cornell University, Ithaca, New York 14853, United States
| | - Siyu Bu
- Department of Chemistry, Cornell University, Ithaca, New York 14853, United States
| | - Britta Johnson
- Department of Chemistry, Cornell University, Ithaca, New York 14853, United States
| | - Nandini Ananth
- Department of Chemistry, Cornell University, Ithaca, New York 14853, United States
| |
Collapse
|
4
|
Lawrence JE, Mannouch JR, Richardson JO. A size-consistent multi-state mapping approach to surface hopping. J Chem Phys 2024; 160:244112. [PMID: 38940540 DOI: 10.1063/5.0208575] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2024] [Accepted: 06/05/2024] [Indexed: 06/29/2024] Open
Abstract
We develop a multi-state generalization of the recently proposed mapping approach to surface hopping (MASH) for the simulation of electronically nonadiabatic dynamics. This new approach extends the original MASH method to be able to treat systems with more than two electronic states. It differs from previous approaches in that it is size consistent and rigorously recovers the original two-state MASH in the appropriate limits. We demonstrate the accuracy of the method by applying it to a series of model systems for which exact benchmark results are available, and we find that the method is well suited to the simulation of photochemical relaxation processes.
Collapse
Affiliation(s)
- Joseph E Lawrence
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8093 Zurich, Switzerland
- Simons Center for Computational Physical Chemistry, New York University, New York, New York 10003, USA
- Department of Chemistry, New York University, New York, New York 10003, USA
| | - Jonathan R Mannouch
- Hamburg Center for Ultrafast Imaging, Universität Hamburg and the Max Planck Institute for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Jeremy O Richardson
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8093 Zurich, Switzerland
| |
Collapse
|
5
|
Huang H, Peng J, Zhang Y, Gu FL, Lan Z, Xu C. The development of the QM/MM interface and its application for the on-the-fly QM/MM nonadiabatic dynamics in JADE package: Theory, implementation, and applications. J Chem Phys 2024; 160:234101. [PMID: 38884395 DOI: 10.1063/5.0215036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2024] [Accepted: 05/15/2024] [Indexed: 06/18/2024] Open
Abstract
Understanding the nonadiabatic dynamics of complex systems is a challenging task in computational photochemistry. Herein, we present an efficient and user-friendly quantum mechanics/molecular mechanics (QM/MM) interface to run on-the-fly nonadiabatic dynamics. Currently, this interface consists of an independent set of codes designed for general-purpose use. Herein, we demonstrate the ability and feasibility of the QM/MM interface by integrating it with our long-term developed JADE package. Tailored to handle nonadiabatic processes in various complex systems, especially condensed phases and protein environments, we delve into the theories, implementations, and applications of on-the-fly QM/MM nonadiabatic dynamics. The QM/MM approach is established within the framework of the additive QM/MM scheme, employing electrostatic embedding, link-atom inclusion, and charge-redistribution schemes to treat the QM/MM boundary. Trajectory surface-hopping dynamics are facilitated using the fewest switches algorithm, encompassing classical and quantum treatments for nuclear and electronic motions, respectively. Finally, we report simulations of nonadiabatic dynamics for two typical systems: azomethane in water and the retinal chromophore PSB3 in a protein environment. Our results not only illustrate the power of the QM/MM program but also reveal the important roles of environmental factors in nonadiabatic processes.
Collapse
Affiliation(s)
- Haiyi Huang
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
- Center for Advanced Materials Research, Beijing Normal University, Zhuhai 519087, China
- MOE Key Laboratory of Theoretical and Computational Photochemistry, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Jiawei Peng
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Yulin Zhang
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Feng Long Gu
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Chao Xu
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
6
|
Mannouch JR, Kelly A. Quantum Quality with Classical Cost: Ab Initio Nonadiabatic Dynamics Simulations Using the Mapping Approach to Surface Hopping. J Phys Chem Lett 2024; 15:5814-5823. [PMID: 38781480 PMCID: PMC11163471 DOI: 10.1021/acs.jpclett.4c00535] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2024] [Revised: 05/07/2024] [Accepted: 05/09/2024] [Indexed: 05/25/2024]
Abstract
Nonadiabatic dynamics methods are an essential tool for investigating photochemical processes. In the context of employing first-principles electronic structure techniques, such simulations can be carried out in a practical manner using semiclassical trajectory-based methods or wave packet approaches. While all approaches applicable to first-principles simulations are necessarily approximate, it is commonly thought that wave packet approaches offer inherent advantages over their semiclassical counterparts in terms of accuracy and that this trait simply comes at a higher computational cost. Here we demonstrate that the mapping approach to surface hopping (MASH), a recently introduced trajectory-based nonadiabatic dynamics method, can be efficiently applied in tandem with ab initio electronic structure. Our results even suggest that MASH may provide more accurate results than on-the-fly wave packet techniques, all at a much lower computational cost.
Collapse
Affiliation(s)
- Jonathan R. Mannouch
- Hamburg Center for Ultrafast
Imaging, Universität Hamburg and
the Max Planck Institute
for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Aaron Kelly
- Hamburg Center for Ultrafast
Imaging, Universität Hamburg and
the Max Planck Institute
for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| |
Collapse
|
7
|
Hutton L, Moreno Carrascosa A, Prentice AW, Simmermacher M, Runeson JE, Paterson MJ, Kirrander A. Using a multistate mapping approach to surface hopping to predict the ultrafast electron diffraction signal of gas-phase cyclobutanone. J Chem Phys 2024; 160:204307. [PMID: 38814011 DOI: 10.1063/5.0203667] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2024] [Accepted: 05/05/2024] [Indexed: 05/31/2024] Open
Abstract
Using the recently developed multistate mapping approach to surface hopping (multistate MASH) method combined with SA(3)-CASSCF(12,12)/aug-cc-pVDZ electronic structure calculations, the gas-phase isotropic ultrafast electron diffraction (UED) of cyclobutanone is predicted and analyzed. After excitation into the n-3s Rydberg state (S2), cyclobutanone can relax through two S2/S1 conical intersections, one characterized by compression of the CO bond and the other by dissociation of the α-CC bond. Subsequent transfer into the ground state (S0) is then achieved via two additional S1/S0 conical intersections that lead to three reaction pathways: α ring-opening, ethene/ketene production, and CO liberation. The isotropic gas-phase UED signal is predicted from the multistate MASH simulations, allowing for a direct comparison to the experimental data. This work, which is a contribution to the cyclobutanone prediction challenge, facilitates the identification of the main photoproducts in the UED signal and thereby emphasizes the importance of dynamics simulations for the interpretation of ultrafast experiments.
Collapse
Affiliation(s)
- Lewis Hutton
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - Andrés Moreno Carrascosa
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - Andrew W Prentice
- Institute of Chemical Sciences, School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh EH14 4AS, United Kingdom
| | - Mats Simmermacher
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - Johan E Runeson
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - Martin J Paterson
- Institute of Chemical Sciences, School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh EH14 4AS, United Kingdom
| | - Adam Kirrander
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
8
|
Liu XY, Wang SR, Fang WH, Cui G. Nuclear Quantum Effects on Nonadiabatic Dynamics of a Green Fluorescent Protein Chromophore Analogue: Ring-Polymer Surface-Hopping Simulation. J Chem Theory Comput 2024; 20:3426-3439. [PMID: 38656202 DOI: 10.1021/acs.jctc.4c00068] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2024]
Abstract
Herein, we have used the "on-the-fly" ring-polymer surface-hopping simulation method with the centroid approximation (RPSH-CA), in combination with the multireference OM2/MRCI electronic structure calculations to study the photoinduced dynamics of a green fluorescent protein (GFP) chromophore analogue in the gas phase, i.e., o-HBI, at 50, 100, and 300 K with 1, 5, 10, and 15 beads (3600 1 ps trajectories). The electronic structure calculations identified five new minimum-energy conical intersection (MECI) structures, which, together with the previous one, play crucial roles in the excited-state decay dynamics of o-HBI. It is also found that the excited-state intramolecular proton transfer (ESIPT) occurs in an ultrafast manner and is completed within 20 fs in all the simulation conditions because there is no barrier associated with this ESIPT process in the S1 state. However, the other excited-state dynamical results are strongly related to the number of beads. At 50 and 100 K, the nuclear quantum effects (NQEs) are very important; therefore, the excited-state dynamical results change significantly with the bead number. For example, the S1 decay time deduced from time-dependent state populations becomes longer as the bead number increases. Nevertheless, an essentially convergent trend is observed when the bead number is close to 10. In contrast, at 300 K, the NQEs become weaker and the above dynamical results converge very quickly even with 1 bead. Most importantly, the NQEs seriously affect the excited-state decay mechanism of o-HBI. At 50 and 100 K, most trajectories decay to the S0 state via perpendicular keto MECIs, whereas, at 300 K, only twisted keto MECIs are responsible for the excited-state decay. The present work not only comprehensively explores the temperature-dependent photoinduced dynamics of o-HBI, but also demonstrates the importance and necessity of NQEs in nonadiabatic dynamics simulations, especially at relatively low temperatures.
Collapse
Affiliation(s)
- Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Sheng-Rui Wang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
9
|
Makri N. Kink Sum for Long-Memory Small Matrix Path Integral Dynamics. J Phys Chem B 2024. [PMID: 38437738 DOI: 10.1021/acs.jpcb.3c08282] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/06/2024]
Abstract
The small matrix decomposition of the real-time path integral (SMatPI) allows for numerically exact and efficient propagation of the reduced density matrix (RDM) for system-bath Hamiltonians. Its high efficiency lies in the small size of the SMatPI matrices employed in the iterative algorithm, whose size is equal to that of the full RDM. By avoiding the storage and multiplication of large tensors, the SMatPI algorithm is applicable in multistate systems under a variety of conditions. The main computational effort is the evaluation of path sums within the entangled memory length to construct the SMatPI matrices. A number of methods are available for this task, each with its own favorable parameter regime, but calculations with strong system-bath coupling and long memory at low temperatures remain out of reach. The present paper evaluates the path sums by binning the paths (in forward time only) based on their amplitudes, which depend on the number and type of kinks they contain. The algorithm is very efficient, leading to a dramatic acceleration of path sums and significantly extending the accessible memory length in the most challenging regimes.
Collapse
Affiliation(s)
- Nancy Makri
- Departments of Chemistry and Physics, University of Illinois, 505 S. Mathews Avenue, Urbana, Illinois 61801, United States
| |
Collapse
|
10
|
Amati G, Mannouch JR, Richardson JO. Detailed balance in mixed quantum-classical mapping approaches. J Chem Phys 2023; 159:214114. [PMID: 38054513 DOI: 10.1063/5.0176291] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2023] [Accepted: 11/07/2023] [Indexed: 12/07/2023] Open
Abstract
The violation of detailed balance poses a serious problem for the majority of current quasiclassical methods for simulating nonadiabatic dynamics. In order to analyze the severity of the problem, we predict the long-time limits of the electronic populations according to various quasiclassical mapping approaches by applying arguments from classical ergodic theory. Our analysis confirms that regions of the mapping space that correspond to negative populations, which most mapping approaches introduce in order to go beyond the Ehrenfest approximation, pose the most serious issue for reproducing the correct thermalization behavior. This is because inverted potentials, which arise from negative electronic populations entering the nuclear force, can result in trajectories unphysically accelerating off to infinity. The recently developed mapping approach to surface hopping (MASH) provides a simple way of avoiding inverted potentials while retaining an accurate description of the dynamics. We prove that MASH, unlike any other quasiclassical approach, is guaranteed to describe the exact thermalization behavior of all quantum-classical systems, confirming it as one of the most promising methods for simulating nonadiabatic dynamics in real condensed-phase systems.
Collapse
Affiliation(s)
- Graziano Amati
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| | - Jonathan R Mannouch
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| | - Jeremy O Richardson
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| |
Collapse
|
11
|
Miyazaki K, Ananth N. Nonadiabatic simulations of photoisomerization and dissociation in ethylene using ab initio classical trajectories. J Chem Phys 2023; 159:124110. [PMID: 38127384 DOI: 10.1063/5.0163371] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2023] [Accepted: 09/05/2023] [Indexed: 12/23/2023] Open
Abstract
We simulate the nonadiabatic dynamics of photo-induced isomerization and dissociation in ethylene using ab initio classical trajectories in an extended phase space of nuclear and electronic variables. This is achieved by employing the linearized semiclassical initial value representation method for nonadiabatic dynamics, where discrete electronic states are mapped to continuous classical variables using either the Meyer-Miller-Stock-Thoss representation or a more recently introduced spin mapping approach. Trajectory initial conditions are sampled by constraining electronic state variables to a single initial excited state and by drawing nuclear phase space configurations from a Wigner distribution at a finite temperature. An ensemble of classical ab initio trajectories is then generated to compute thermal population correlation functions and analyze the mechanisms of isomerization and dissociation. Our results serve as a demonstration that this parameter-free semiclassical approach is computationally efficient and accurate, identifying mechanistic pathways in agreement with previous theoretical studies and also uncovering dissociation pathways observed experimentally.
Collapse
Affiliation(s)
- K Miyazaki
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - N Ananth
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
12
|
Abstract
We present a nonadiabatic classical-trajectory approach that offers the best of both worlds between fewest-switches surface hopping (FSSH) and quasiclassical mapping dynamics. This mapping approach to surface hopping (MASH) propagates the nuclei on the active adiabatic potential-energy surface, such as in FSSH. However, unlike in FSSH, transitions between active surfaces are deterministic and occur when the electronic mapping variables evolve between specified regions of the electronic phase space. This guarantees internal consistency between the active surface and the electronic degrees of freedom throughout the dynamics. MASH is rigorously derivable from exact quantum mechanics as a limit of the quantum-classical Liouville equation (QCLE), leading to a unique prescription for momentum rescaling and frustrated hops. Hence, a quantum-jump procedure can, in principle, be used to systematically converge the accuracy of the results to that of the QCLE. This jump procedure also provides a rigorous framework for deriving approximate decoherence corrections similar to those proposed for FSSH. We apply MASH to simulate the nonadiabatic dynamics in various model systems and show that it consistently produces more accurate results than FSSH at a comparable computational cost.
Collapse
|
13
|
Mousavi SZ, Shadman HR, Habibi M, Didandeh M, Nikzad A, Golmohammadi M, Maleki R, Suwaileh WA, Khataee A, Zargar M, Razmjou A. Elucidating the Sorption Mechanisms of Environmental Pollutants Using Molecular Simulation. Ind Eng Chem Res 2023. [DOI: 10.1021/acs.iecr.2c02333] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/22/2023]
Affiliation(s)
- Seyedeh Zahra Mousavi
- Department of Chemical Engineering, Tarbiat Modares University, Tehran, 1411944961, Iran
| | - Hamid Reza Shadman
- Department of Polymer Engineering & Color Technology, Amirkabir University of Technology, Tehran, 6351713178, Iran
| | - Meysam Habibi
- Department of Chemical Engineering, University of Tehran, Tehran, 6718773654, Iran
| | - Mohsen Didandeh
- Department of Chemical Engineering, Tarbiat Modares University, Tehran, 1411944961, Iran
| | - Arash Nikzad
- Mechanical Engineering Department, University of British Columbia, Vancouver, BC V6T 1Z4, Canada
| | - Mahsa Golmohammadi
- Department of Polymer Engineering & Color Technology, Amirkabir University of Technology, Tehran, 6351713178, Iran
| | - Reza Maleki
- Department of Chemical Technologies, Iranian Research Organization for Science and Technology (IROST), P.O. Box 33535111, Tehran, 3313193685, Iran
| | - Wafa Ali Suwaileh
- Chemical Engineering Program, Texas A&M University at Qatar, Education City, Doha 23874, Qatar
| | - Alireza Khataee
- Research Laboratory of Advanced Water and Wastewater Treatment Processes, Department of Applied Chemistry, Faculty of Chemistry, University of Tabriz, 51666-16471 Tabriz, Iran
- Department of Materials Science and Nanotechnology Engineering, Faculty of Engineering, Near East University, 99138 Nicosia, Mersin 10 Turkey
| | - Masoumeh Zargar
- Mineral Recovery Research Center (MRRC), School of Engineering, Edith Cowan University, Joondalup, Perth WA 6027, Australia
| | - Amir Razmjou
- Mineral Recovery Research Center (MRRC), School of Engineering, Edith Cowan University, Joondalup, Perth WA 6027, Australia
- UNESCO Centre for Membrane Science and Technology, School of Chemical Engineering, University of New South Wales, Sydney, NSW 2052, Australia
| |
Collapse
|
14
|
Amati G, Runeson JE, Richardson JO. On detailed balance in nonadiabatic dynamics: From spin spheres to equilibrium ellipsoids. J Chem Phys 2023; 158:064113. [PMID: 36792511 DOI: 10.1063/5.0137828] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Abstract
Trajectory-based methods that propagate classical nuclei on multiple quantum electronic states are often used to simulate nonadiabatic processes in the condensed phase. A long-standing problem of these methods is their lack of detailed balance, meaning that they do not conserve the equilibrium distribution. In this article, we investigate ideas for restoring detailed balance in mixed quantum-classical systems by tailoring the previously proposed spin-mapping approach to thermal equilibrium. We find that adapting the spin magnitude can recover the correct long-time populations but is insufficient to conserve the full equilibrium distribution. The latter can however be achieved by a more flexible mapping of the spin onto an ellipsoid, which is constructed to fulfill detailed balance for arbitrary potentials. This ellipsoid approach solves the problem of negative populations that has plagued previous mapping approaches and can therefore be applied also to strongly asymmetric and anharmonic systems. Because it conserves the thermal distribution, the method can also exploit efficient sampling schemes used in standard molecular dynamics, which drastically reduces the number of trajectories needed for convergence. The dynamics does however still have mean-field character, as is observed most clearly by evaluating reaction rates in the golden-rule limit. This implies that although the ellipsoid mapping provides a rigorous framework, further work is required to find an accurate classical-trajectory approximation that captures more properties of the true quantum dynamics.
Collapse
Affiliation(s)
- Graziano Amati
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | - Johan E Runeson
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | |
Collapse
|
15
|
Bossion D, Chowdhury SN, Huo P. Non-adiabatic ring polymer molecular dynamics in the phase space of the SU(N) Lie group. J Chem Phys 2023; 158:044123. [PMID: 36725494 DOI: 10.1063/5.0133970] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
We derive the non-adiabatic ring polymer molecular dynamics (RPMD) approach in the phase space of the SU(N) Lie Group. This method, which we refer to as the spin mapping non-adiabatic RPMD (SM-NRPMD), is based on the spin-mapping formalism for the electronic degrees of freedom (DOFs) and ring polymer path-integral description for the nuclear DOFs. Using the Stratonovich-Weyl transform for the electronic DOFs and the Wigner transform for the nuclear DOFs, we derived an exact expression of the Kubo-transformed time-correlation function (TCF). We further derive the spin mapping non-adiabatic Matsubara dynamics using the Matsubara approximation that removes the high frequency nuclear normal modes in the TCF and derive the SM-NRPMD approach from the non-adiabatic Matsubara dynamics by discarding the imaginary part of the Liouvillian. The SM-NRPMD method has numerical advantages compared to the original NRPMD method based on the Meyer-Miller-Stock-Thoss (MMST) mapping formalism due to a more natural mapping using the SU(N) Lie Group that preserves the symmetry of the original system. We numerically compute the Kubo-transformed position auto-correlation function and electronic population correlation function for three-state model systems. The numerical results demonstrate the accuracy of the SM-NRPMD method, which outperforms the original MMST-based NRPMD. We envision that the SM-NRPMD method will be a powerful approach to simulate electronic non-adiabatic dynamics and nuclear quantum effects accurately.
Collapse
Affiliation(s)
- Duncan Bossion
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| | - Sutirtha N Chowdhury
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| |
Collapse
|
16
|
Lin K, Peng J, Xu C, Gu FL, Lan Z. Trajectory Propagation of Symmetrical Quasi-classical Dynamics with Meyer-Miller Mapping Hamiltonian Using Machine Learning. J Phys Chem Lett 2022; 13:11678-11688. [PMID: 36511563 DOI: 10.1021/acs.jpclett.2c02159] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
The long short-term memory recurrent neural network (LSTM-RNN) approach is applied to realize the trajectory-based nonadiabatic dynamics within the framework of the symmetrical quasi-classical dynamics method based on the Meyer-Miller mapping Hamiltonian (MM-SQC). After construction, the LSTM-RNN model allows us to propagate the entire trajectory evolutions of all involved degrees of freedoms (DOFs) from initial conditions. The proposed idea is proven to be reliable and accurate in the simulations of the dynamics of several site-exciton electron-phonon coupling models and three Tully's scattering models. It indicates that the LSTM-RNN model perfectly captures the dynamical information on the trajectory evolution in the MM-SQC dynamics. Our work proposes a novel machine learning approach in the simulation of trajectory-based nonadiabatic dynamic of complex systems with a large number of DOFs.
Collapse
Affiliation(s)
- Kunni Lin
- School of Chemistry, South China Normal University, Guangzhou 510006, P. R. China
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| | - Jiawei Peng
- School of Chemistry, South China Normal University, Guangzhou 510006, P. R. China
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| | - Chao Xu
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Feng Long Gu
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Zhenggang Lan
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| |
Collapse
|
17
|
Amati G, Saller MAC, Kelly A, Richardson JO. Quasiclassical approaches to the generalized quantum master equation. J Chem Phys 2022; 157:234103. [PMID: 36550031 DOI: 10.1063/5.0124028] [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/29/2022] Open
Abstract
The formalism of the generalized quantum master equation (GQME) is an effective tool to simultaneously increase the accuracy and the efficiency of quasiclassical trajectory methods in the simulation of nonadiabatic quantum dynamics. The GQME expresses correlation functions in terms of a non-Markovian equation of motion, involving memory kernels that are typically fast-decaying and can therefore be computed by short-time quasiclassical trajectories. In this paper, we study the approximate solution of the GQME, obtained by calculating the kernels with two methods: Ehrenfest mean-field theory and spin-mapping. We test the approaches on a range of spin-boson models with increasing energy bias between the two electronic levels and place a particular focus on the long-time limits of the populations. We find that the accuracy of the predictions of the GQME depends strongly on the specific technique used to calculate the kernels. In particular, spin-mapping outperforms Ehrenfest for all the systems studied. The problem of unphysical negative electronic populations affecting spin-mapping is resolved by coupling the method with the master equation. Conversely, Ehrenfest in conjunction with the GQME can predict negative populations, despite the fact that the populations calculated from direct dynamics are positive definite.
Collapse
Affiliation(s)
- Graziano Amati
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | | - Aaron Kelly
- The Hamburg Center for Ultrafast Imaging, Universität Hamburg, Max Planck Institute for the Structure and Dynamics of Matter and Center for Free-Electron Laser Science, 22761 Hamburg, Germany
| | | |
Collapse
|
18
|
Rivlin T, Pollak E. Nonadiabatic Couplings Can Speed Up Quantum Tunneling Transition Path Times. J Phys Chem Lett 2022; 13:10558-10566. [PMID: 36342976 PMCID: PMC9677498 DOI: 10.1021/acs.jpclett.2c03008] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2022] [Accepted: 11/02/2022] [Indexed: 06/16/2023]
Abstract
Quantum tunneling is known to play an important role in the dynamics of systems with nonadiabatic couplings. However, until recently, the time-domain properties of nonadiabatic scattering have been severely under-explored. Using numerically exact quantum methods, we study the impact that nonadiabatic couplings have on the time it takes to tunnel through a barrier. We find that the Wigner phase time is the appropriate measure to use when determining the tunneling flight time also when considering nonadiabatic systems. The central result of the present study is that in an avoided crossing system in one dimension, the nonadiabatic couplings speed up the tunneling event, relative to the adiabatic case in which all nonadiabatic coupling is ignored. This has implications for both the study of quantum tunneling times and for the field of nonadiabatic scattering and chemistry.
Collapse
|
19
|
Lin K, Peng J, Xu C, Gu FL, Lan Z. Automatic Evolution of Machine-Learning-Based Quantum Dynamics with Uncertainty Analysis. J Chem Theory Comput 2022; 18:5837-5855. [DOI: 10.1021/acs.jctc.2c00702] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Kunni Lin
- School of Chemistry, South China Normal University, Guangzhou510006, P. R. China
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
| | - Jiawei Peng
- School of Chemistry, South China Normal University, Guangzhou510006, P. R. China
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
| | - Chao Xu
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou510006, P. R. China
| | - Feng Long Gu
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou510006, P. R. China
| | - Zhenggang Lan
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou510006, P. R. China
| |
Collapse
|
20
|
He X, Wu B, Rivlin T, Liu J, Pollak E. Transition Path Flight Times and Nonadiabatic Electronic Transitions. J Phys Chem Lett 2022; 13:6966-6974. [PMID: 35877977 PMCID: PMC9358656 DOI: 10.1021/acs.jpclett.2c01425] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2022] [Accepted: 07/18/2022] [Indexed: 06/15/2023]
Abstract
Transition path flight times are studied for scattering on two electronic surfaces with a single crossing. These flight times reveal nontrivial quantum effects such as resonance lifetimes and nonclassical passage times and reveal that nonadiabatic effects often increase flight times. The flight times are computed using numerically exact time propagation and compared with results obtained from the Fewest Switches Surface Hopping (FSSH) method. Comparison of the two methods shows that the FSSH method is reliable for transition path times only when the scattering is classically allowed on the relevant adiabatic surfaces. However, where quantum effects such as tunneling and resonances dominate, the FSSH method is not adequate to accurately predict the correct times and transition probabilities. These results highlight limitations in methods which do not account for quantum interference effects, and suggest that measuring flight times is important for obtaining insights from the time-domain into quantum effects in nonadiabatic scattering.
Collapse
Affiliation(s)
- Xin He
- Beijing
National Laboratory for Molecular Sciences, Institute of Theoretical
and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Baihua Wu
- Beijing
National Laboratory for Molecular Sciences, Institute of Theoretical
and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Tom Rivlin
- Chemical
and Biological Physics Department, Weizmann
Institute of Science, 76100 Rehovot, Israel
| | - Jian Liu
- Beijing
National Laboratory for Molecular Sciences, Institute of Theoretical
and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Eli Pollak
- Chemical
and Biological Physics Department, Weizmann
Institute of Science, 76100 Rehovot, Israel
| |
Collapse
|
21
|
Bossion DL, Ying W, Chowdhury S, Huo P. Non-adiabatic Mapping Dynamics in the Phase Space of the ${SU}(N)$ Lie Group. J Chem Phys 2022; 157:084105. [DOI: 10.1063/5.0094893] [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 the rigorous theoretical framework of the generalized spin mapping representation for non-adiabatic dynamics. Our work is based up a new mapping formalism recently introduced by Runeson and Richardson in [J. Chem. Phys. 152, 084110 (2020)], which uses the generators of the su(N) Lie algebra to represent N discrete electronic states, thus preserving the size of the original Hilbert space. Following this interesting idea, the Stratonovich-Weyl transform is used to map an operator in the Hilbert space to a continuous function on the SU(N) Lie group, i.e., a smooth manifold which is a phase space of continuous variables. We further use the Wigner representation to describe the nuclear degrees of freedom, and derived an exact expression of the time-correlation function as well as the exact quantum Liouvillian for the non-adiabatic system. Making the linearization approximation, this exact Liouvillian is reduced to the Liouvillian of several recently proposed methods, and the performance of this Linearized method is tested using non-adiabatic models. We envision that the theoretical work presented here provides a rigorous and unified framework to formally derive non-adiabatic quantum dynamics approaches with continuous variables and connect the previous methods in a clear and concise manner.
Collapse
Affiliation(s)
| | | | - Sutirtha Chowdhury
- Department of Chemistry, University of Rochester, United States of America
| | - Pengfei Huo
- Department of Chemsitry, University of Rochester Department of Chemistry, United States of America
| |
Collapse
|
22
|
He X, Wu B, Shang Y, Li B, Cheng X, Liu J. New phase space formulations and quantum dynamics approaches. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1619] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Xin He
- Beijing National Laboratory for Molecular Sciences Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University Beijing China
| | - Baihua Wu
- Beijing National Laboratory for Molecular Sciences Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University Beijing China
| | - Youhao Shang
- Beijing National Laboratory for Molecular Sciences Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University Beijing China
| | - Bingqi Li
- Beijing National Laboratory for Molecular Sciences Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University Beijing China
| | - Xiangsong Cheng
- Beijing National Laboratory for Molecular Sciences Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University Beijing China
| | - Jian Liu
- Beijing National Laboratory for Molecular Sciences Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University Beijing China
| |
Collapse
|
23
|
Runeson JE, Lawrence JE, Mannouch JR, Richardson JO. Explaining the Efficiency of Photosynthesis: Quantum Uncertainty or Classical Vibrations? J Phys Chem Lett 2022; 13:3392-3399. [PMID: 35404611 PMCID: PMC9036581 DOI: 10.1021/acs.jpclett.2c00538] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2022] [Accepted: 04/01/2022] [Indexed: 06/14/2023]
Abstract
Photosynthetic organisms are known to use a mechanism of vibrationally assisted exciton energy transfer to efficiently harvest energy from light. The importance of quantum effects in this mechanism is a long-standing topic of debate, which has traditionally focused on the role of excitonic coherences. Here, we address another recent claim: that the efficient energy transfer in the Fenna-Matthews-Olson complex relies on nuclear quantum uncertainty and would not function if the vibrations were classical. We present a counter-example to this claim, showing by trajectory-based simulations that a description in terms of quantum electrons and classical nuclei is indeed sufficient to describe the funneling of energy to the reaction center. We analyze and compare these findings to previous classical-nuclear approximations that predicted the absence of an energy funnel and conclude that the key difference and the reason for the discrepancy is the ability of the trajectories to properly account for Newton's third law.
Collapse
Affiliation(s)
- Johan E. Runeson
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | | | | | |
Collapse
|
24
|
Abstract
The quantum-classical path integral (QCPI) is a rigorous formulation of nonadiabatic dynamics, where the dynamical interaction between a quantum system and its environment is captured consistently through classical trajectories driven by forces along quantum paths of the system. In this Letter, we develop a small matrix decomposition (SMatQCPI) that eliminates the tensor storage requirements of the iterative QCPI algorithm. In the case of a system coupled to a harmonic bath, SMatQCPI provides fully quantum mechanical propagation, which also reduces the computational cost to that of a single QCPI step. Further, the SMatQCPI matrices only need to account for quantum contributions to decoherence, allowing high efficiency in challenging regimes of incoherent dynamics. Overall, this new composite algorithm combines the best features of two powerful path integral formulations and offers a versatile tool for simulating condensed phase quantum dynamics.
Collapse
Affiliation(s)
- Sohang Kundu
- Department of Chemistry, University of Illinois, Urbana, Illinois 61801, United States
| | - Nancy Makri
- Department of Chemistry, University of Illinois, Urbana, Illinois 61801, United States
- Department of Physics, University of Illinois, Urbana, Illinois 61801, United States
- Illinois Quantum Information Science and Technology Center, University of Illinois, Urbana, Illinois 61801, United States
| |
Collapse
|
25
|
Ionic Liquid@Metal-Organic Framework as a Solid Electrolyte in a Lithium-Ion Battery: Current Performance and Perspective at Molecular Level. NANOMATERIALS 2022; 12:nano12071076. [PMID: 35407194 PMCID: PMC9000457 DOI: 10.3390/nano12071076] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/08/2022] [Revised: 02/27/2022] [Accepted: 03/07/2022] [Indexed: 01/27/2023]
Abstract
Searching for a suitable electrolyte in a lithium-ion battery is a challenging task. The electrolyte must not only be chemically and mechanically stable, but also be able to transport lithium ions efficiently. Ionic liquid incorporated into a metal-organic framework (IL@MOF) has currently emerged as an interesting class of hybrid material that could offer excellent electrochemical properties. However, the understanding of the mechanism and factors that govern its fast ionic conduction is crucial as well. In this review, the characteristics and potential use of IL@MOF as an electrolyte in a lithium-ion battery are highlighted. The importance of computational methods is emphasized as a comprehensive tool to investigate the atomistic behavior of IL@MOF and its interaction in electrochemical environments.
Collapse
|
26
|
Talbot JJ, Head-Gordon M, Miller WH, Cotton SJ. Dynamic signatures of electronically nonadiabatic coupling in sodium hydride: a rigorous test for the symmetric quasi-classical model applied to realistic, ab initio electronic states in the adiabatic representation. Phys Chem Chem Phys 2022; 24:4820-4831. [PMID: 35156112 DOI: 10.1039/d1cp04090a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Sodium hydride (NaH) in the gas phase presents a seemingly simple electronic structure making it a potentially tractable system for the detailed investigation of nonadiabatic molecular dynamics from both computational and experimental standpoints. The single vibrational degree of freedom, as well as the strong nonadiabatic coupling that arises from the excited electronic states taking on considerable ionic character, provides a realistic chemical system to test the accuracy of quasi-classical methods to model population dynamics where the results are directly comparable against quantum mechanical benchmarks. Using a simulated pump-probe type experiment, this work presents computational predictions of population transfer through the avoided crossings of NaH via symmetric quasi-classical Meyer-Miller (SQC/MM), Ehrenfest, and exact quantum dynamics on realistic, ab initio potential energy surfaces. The main driving force for population transfer arises from the ground vibrational level of the D1Σ+ adiabatic state that is embedded in the manifold of near-dissociation C1Σ+ vibrational states. When coupled through a sharply localized first-order derivative coupling most of the population transfers between t = 15 and t = 30 fs depending on the initially excited vibronic wavepacket. While quantum mechanical effects are expected due to the reduced mass of NaH, predictions of the population dynamics from both the SQC/MM and Ehrenfest models perform remarkably well against the quantum dynamics benchmark. Additionally, an analysis of the vibronic structure in the nonadiabatically coupled regime is presented using a variational eigensolver methodology.
Collapse
Affiliation(s)
- Justin J Talbot
- Department of Chemistry, University of California, Berkeley, California 94720, USA.
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA. .,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California, 94720, USA
| | - William H Miller
- Department of Chemistry, University of California, Berkeley, California 94720, USA.
| | - Stephen J Cotton
- Department of Chemistry, University of California, Berkeley, California 94720, USA.
| |
Collapse
|
27
|
Mannouch JR, Richardson JO. A partially linearized spin-mapping approach for simulating nonlinear optical spectra. J Chem Phys 2022; 156:024108. [PMID: 35032975 DOI: 10.1063/5.0077744] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a partially linearized method based on spin-mapping for computing both linear and nonlinear optical spectra. As observables are obtained from ensembles of classical trajectories, the approach can be applied to the large condensed-phase systems that undergo photosynthetic light-harvesting processes. In particular, the recently derived spin partially linearized density matrix method has been shown to exhibit superior accuracy in computing population dynamics compared to other related classical-trajectory methods. Such a method should also be ideally suited to describing the quantum coherences generated by interaction with light. We demonstrate that this is, indeed, the case by calculating the nonlinear optical response functions relevant for the pump-probe and 2D photon-echo spectra for a Frenkel biexciton model and the Fenna-Matthews-Olsen light-harvesting complex. One especially desirable feature of our approach is that the full spectrum can be decomposed into its constituent components associated with the various Liouville-space pathways, offering a greater insight beyond what can be directly obtained from experiments.
Collapse
|
28
|
Cho KH, Rhee YM. Computational elucidations on the role of vibrations in energy transfer processes of photosynthetic complexes. Phys Chem Chem Phys 2021; 23:26623-26639. [PMID: 34842245 DOI: 10.1039/d1cp04615b] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Coupling between pigment excitations and nuclear movements in photosynthetic complexes is known to modulate the excitation energy transfer (EET) efficiencies. Toward providing microscopic information, researchers often apply simulation techniques and investigate how vibrations are involved in EET processes. Here, reports on such roles of nuclear movements are discussed from a theory perspective. While vibrations naturally present random thermal fluctuations that can affect energy transferring characteristics, they can also be intertwined with exciton structures and create more specific non-adiabatic energy transfer pathways. For reliable simulations, a bath model that accurately mimics a given molecular system is required. Methods for obtaining such a model in combination with quantum chemical electronic structure calculations and molecular dynamics trajectory simulations are discussed. Various quantum dynamics simulation tools that can handle pigment-to-pigment energy transfers together with their vibrational characters are also touched on. Behaviors of molecular vibrations often deviate from ideality, especially when all-atom details are included, which practically forces us to treat them classically. We conclude this perspective by considering some recent reports that suggest that classical descriptions of bath effects with all-atom details may still produce valuable information for analyzing sophisticated contributions by vibrations to EET processes.
Collapse
Affiliation(s)
- Kwang Hyun Cho
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| | - Young Min Rhee
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| |
Collapse
|
29
|
Liu J, He X, Wu B. Unified Formulation of Phase Space Mapping Approaches for Nonadiabatic Quantum Dynamics. Acc Chem Res 2021; 54:4215-4228. [PMID: 34756027 DOI: 10.1021/acs.accounts.1c00511] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Nonadiabatic dynamical processes are one of the most important quantum mechanical phenomena in chemical, materials, biological, and environmental molecular systems, where the coupling between different electronic states is either inherent in the molecular structure or induced by the (intense) external field. The curse of dimensionality indicates the intractable exponential scaling of calculation effort with system size and restricts the implementation of "numerically exact" approaches for realistic large systems. The phase space formulation of quantum mechanics offers an important theoretical framework for constructing practical approximate trajectory-based methods for quantum dynamics. This Account reviews our recent progress in phase space mapping theory: a unified framework for constructing the mapping Hamiltonian on phase space for coupled F-state systems where the renowned Meyer-Miller Hamiltonian model is a special case, a general phase space formulation of quantum mechanics for nonadiabatic systems where the electronic degrees of freedom are mapped onto constraint space and the nuclear degrees of freedom are mapped onto infinite space, and an isomorphism between the mapping phase space approach for nonadiabatic systems and that for nonequilibrium electron transport processes. While the zero-point-energy parameter is conventionally assumed to be positive, we show that the constraint implied in the conventional Meyer-Miller mapping Hamiltonian requires that such a parameter can be negative as well and lies in (-1/F, +∞) for each electronic degree of freedom. More importantly, the zero-point-energy parameter should be interpreted as a special case of a commutator matrix in the comprehensive phase space mapping Hamiltonian for nonadiabatic systems. From the rigorous formulation of mapping phase space, we propose approximate but practical trajectory-based nonadiabatic dynamics methods. The applications to both gas phase and condensed phase problems include the spin-boson model for condensed phase dissipative two-state systems, the three-state photodissociation models, the seven-site model of the Fenna-Matthews-Olson monomer in photosynthesis of green sulfur bacteria, the strongly coupled molecular/atomic matter-optical cavity systems designed for controlling and manipulating chemical dynamical processes, and the Landauer model for a quantum dot state coupled with two electrodes. In these applications the overall performance of our phase space mapping dynamics approach is superior to two prevailing trajectory-based methods, Ehrenfest dynamics and fewest switches surface hopping.
Collapse
Affiliation(s)
- Jian Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Xin He
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Baihua Wu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
30
|
Lin K, Peng J, Gu FL, Lan Z. Simulation of Open Quantum Dynamics with Bootstrap-Based Long Short-Term Memory Recurrent Neural Network. J Phys Chem Lett 2021; 12:10225-10234. [PMID: 34647736 DOI: 10.1021/acs.jpclett.1c02672] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
The recurrent neural network with the long short-term memory cell (LSTM-NN) is employed to simulate the long-time dynamics of open quantum systems. The bootstrap method is applied in the LSTM-NN construction and prediction, which provides a Monte Carlo estimation of a forecasting confidence interval. Within this approach, a large number of LSTM-NNs are constructed by resampling time-series sequences that were obtained from the early stage quantum evolution given by numerically exact multilayer multiconfigurational time-dependent Hartree method. The built LSTM-NN ensemble is used for the reliable propagation of the long-time quantum dynamics, and the simulated result is highly consistent with the exact evolution. The forecasting uncertainty that partially reflects the reliability of the LSTM-NN prediction is also given. This demonstrates the bootstrap-based LSTM-NN approach is a practical and powerful tool to propagate the long-time quantum dynamics of open systems with high accuracy and low computational cost.
Collapse
Affiliation(s)
- Kunni Lin
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| | - Jiawei Peng
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Feng Long Gu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| | - Zhenggang Lan
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| |
Collapse
|
31
|
Weight BM, Mandal A, Huo P. Ab initio symmetric quasi-classical approach to investigate molecular Tully models. J Chem Phys 2021; 155:084106. [PMID: 34470343 DOI: 10.1063/5.0061934] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We perform on-the-fly non-adiabatic molecular dynamics simulations using the symmetrical quasi-classical (SQC) approach with the recently suggested molecular Tully models: ethylene and fulvene. We attempt to provide benchmarks of the SQC methods using both the square and triangle windowing schemes as well as the recently proposed electronic zero-point-energy correction scheme (the so-called γ correction). We use the quasi-diabatic propagation scheme to directly interface the diabatic SQC methods with adiabatic electronic structure calculations. Our results showcase the drastic improvement of the accuracy by using the trajectory-adjusted γ-corrections, which outperform the widely used trajectory surface hopping method with decoherence corrections. These calculations provide useful and non-trivial tests to systematically investigate the numerical performance of various diabatic quantum dynamics approaches, going beyond simple diabatic model systems that have been used as the major workhorse in the quantum dynamics field. At the same time, these available benchmark studies will also likely foster the development of new quantum dynamics approaches based on these techniques.
Collapse
Affiliation(s)
- Braden M Weight
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA
| | - Arkajit Mandal
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
32
|
Lang H, Vendrell O, Hauke P. Generalized discrete truncated Wigner approximation for nonadiabatic quantum-classical dynamics. J Chem Phys 2021; 155:024111. [PMID: 34266254 DOI: 10.1063/5.0054696] [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/15/2022] Open
Abstract
Nonadiabatic molecular dynamics occur in a wide range of chemical reactions and femtochemistry experiments involving electronically excited states. These dynamics are hard to treat numerically as the system's complexity increases, and it is thus desirable to have accurate yet affordable methods for their simulation. Here, we introduce a linearized semiclassical method, the generalized discrete truncated Wigner approximation (GDTWA), which is well-established in the context of quantum spin lattice systems, into the arena of chemical nonadiabatic systems. In contrast to traditional continuous mapping approaches, e.g., the Meyer-Miller-Stock-Thoss and the spin mappings, GDTWA samples the electron degrees of freedom in a discrete phase space and thus forbids an unphysical unbounded growth of electronic state populations. The discrete sampling also accounts for an effective reduced but non-vanishing zero-point energy without an explicit parameter, which makes it possible to treat the identity operator and other operators on an equal footing. As numerical benchmarks on two linear vibronic coupling models and Tully's models show, GDTWA has a satisfactory accuracy in a wide parameter regime, independent of whether the dynamics is dominated by relaxation or by coherent interactions. Our results suggest that the method can be very adequate to treat challenging nonadiabatic dynamics problems in chemistry and related fields.
Collapse
Affiliation(s)
- Haifeng Lang
- Theoretical Chemistry, Institute of Physical Chemistry, Heidelberg University, Im Neuenheimer Feld 229, 69120 Heidelberg, Germany
| | - Oriol Vendrell
- Theoretical Chemistry, Institute of Physical Chemistry, Heidelberg University, Im Neuenheimer Feld 229, 69120 Heidelberg, Germany
| | - Philipp Hauke
- INO-CNR BEC Center and Department of Physics, University of Trento, Via Sommarive 14, I-38123 Trento, Italy
| |
Collapse
|
33
|
Koopman J, Grimme S. From QCEIMS to QCxMS: A Tool to Routinely Calculate CID Mass Spectra Using Molecular Dynamics. JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY 2021; 32:1735-1751. [PMID: 34080847 DOI: 10.1021/jasms.1c00098] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
Mass spectrometry (MS) is a powerful tool in chemical research and substance identification. For the computational modeling of electron ionization MS, we have developed the quantum-chemical electron ionization mass spectra (QCEIMS) program. Here, we present an extension of QCEIMS to calculate collision-induced dissociation (CID) spectra. The more general applicability is accounted for by the new name QCxMS, where "x" refers to EI or CID. To this end, fragmentation and rearrangement reactions are computed "on-the-fly" in Born-Oppenheimer molecular dynamics (MD) simulations with the semiempirical GFN2-xTB Hamiltonian, which provides an efficient quantum mechanical description of all elements up to Z = 86 (Rn). Through the explicit modeling of multicollision processes between precursor ions and neutral gas atoms as well as temperature-induced decomposition reactions, QCxMS provides detailed insight into the collision kinetics and fragmentation pathways. In combination with the CREST program to determine the preferential protonation sites, QCxMS becomes the first standalone MD-based program that can predict mass spectra based solely on molecular structures as input. We demonstrate this for six organic molecules with masses ranging from 159 to 296 Da, for which QCxMS yields CID spectra in reasonable agreement with experiments.
Collapse
Affiliation(s)
- Jeroen Koopman
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| |
Collapse
|
34
|
Hantz ER, Lindert S. Adaptative Steered Molecular Dynamics Study of Mutagenesis Effects on Calcium Affinity in the Regulatory Domain of Cardiac Troponin C. J Chem Inf Model 2021; 61:3052-3057. [PMID: 34080877 DOI: 10.1021/acs.jcim.1c00419] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Calcium-dependent cardiac muscle contraction is regulated by the protein complex troponin (cTn) and specifically by the regulatory N-terminal domain (N-cTnC) which contains one active Ca2+ binding site (site II). It has been previously shown that cardiac muscle contractility and functionality is affected by mutations in N-cTnC which alter calcium binding affinity. Here, we describe the application of adaptive steered molecular dynamics to characterize the influence of N-cTnC mutations on site II calcium binding affinity. We observed the correct trends for all of the studied calcium sensitizing and desensitizing mutants, in conjunction with loop II perturbations. Additionally, the potential of mean force accuracy was shown to increase substantially with increasingly slower speeds and using fewer trajectories. This study presents a novel approach to computationally estimate the Ca2+ binding affinity of N-cTnC structures and is a valuable potential tool to support the design and characterization of novel mutations with potential therapeutic benefits.
Collapse
Affiliation(s)
- Eric R Hantz
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Steffen Lindert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
35
|
Hu D, Xie Y, Peng J, Lan Z. On-the-Fly Symmetrical Quasi-Classical Dynamics with Meyer-Miller Mapping Hamiltonian for the Treatment of Nonadiabatic Dynamics at Conical Intersections. J Chem Theory Comput 2021; 17:3267-3279. [PMID: 34028268 DOI: 10.1021/acs.jctc.0c01249] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The on-the-fly version of the symmetrical quasi-classical dynamics method based on the Meyer-Miller mapping Hamiltonian (SQC/MM) is implemented to study the nonadiabatic dynamics at conical intersections of polyatomic systems. The current on-the-fly implementation of the SQC/MM method is based on the adiabatic representation and the dressed momentum. To include the zero-point energy (ZPE) correction of the electronic mapping variables, we employ both the γ-adjusted and γ-fixed approaches. Nonadiabatic dynamics of the methaniminium cation (CH2NH2+) and azomethane are simulated using the on-the-fly SQC/MM method. For CH2NH2+, both ZPE correction approaches give reasonable and consistent results. However, for azomethane, the γ-adjusted version of the SQC/MM dynamics behaves much better than the γ-fixed version. Further analysis indicates that it is always recommended to use the γ-adjusted SQC/MM dynamics in the on-the-fly simulation of photoinduced dynamics of polyatomic systems, particularly when the excited state is well separated from the ground state in the Franck-Condon region. This work indicates that the on-the-fly SQC/MM method is a powerful simulation protocol to deal with the nonadiabatic dynamics of realistic polyatomic systems.
Collapse
Affiliation(s)
- Deping Hu
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China.,School of Environment, South China Normal University, Guangzhou 510006, China
| | - Yu Xie
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China.,School of Environment, South China Normal University, Guangzhou 510006, China
| | - Jiawei Peng
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China.,School of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China.,School of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
36
|
Chowdhury SN, Huo P. Non-adiabatic Matsubara dynamics and non-adiabatic ring-polymer molecular dynamics. J Chem Phys 2021; 154:124124. [PMID: 33810665 DOI: 10.1063/5.0042136] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
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 TCF using the Wigner representation for both the nuclear degrees of freedom and the electronic mapping variables. By dropping the non-Matsubara nuclear normal modes in the quantum Liouvillian and explicitly integrating these modes out from the expression 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 TCF. 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 methods in the future.
Collapse
Affiliation(s)
- Sutirtha N Chowdhury
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| |
Collapse
|
37
|
He X, Gong Z, Wu B, Liu J. Negative Zero-Point-Energy Parameter in the Meyer-Miller Mapping Model for Nonadiabatic Dynamics. J Phys Chem Lett 2021; 12:2496-2501. [PMID: 33667108 DOI: 10.1021/acs.jpclett.1c00232] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
The celebrated Meyer-Miller mapping model has been a useful approach for generating practical trajectory-based nonadiabatic dynamics methods. It is generally assumed that the zero-point-energy (ZPE) parameter is positive. The constraint implied in the conventional Meyer-Miller mapping Hamiltonian for an F-electronic-state system actually requires γ∈(-1/F, ∞) for the ZPE parameter for each electronic degree of freedom. Both negative and positive values are possible for such a parameter. We first establish a rigorous formulation to construct exact mapping models in the Cartesian phase space when the constraint is applied. When nuclear dynamics is approximated by the linearized semiclassical initial value representation, a negative ZPE parameter could lead to reasonably good performance in describing dynamic behaviors in typical spin-boson models for condensed-phase two-state systems, even at challenging zero temperature.
Collapse
Affiliation(s)
- Xin He
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Zhihao Gong
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Baihua Wu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Jian Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
38
|
Peng J, Xie Y, Hu D, Lan Z. Analysis of bath motion in MM-SQC dynamics via dimensionality reduction approach: Principal component analysis. J Chem Phys 2021; 154:094122. [PMID: 33685149 DOI: 10.1063/5.0039743] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The system-plus-bath model is an important tool to understand the nonadiabatic dynamics of large molecular systems. Understanding the collective motion of a large number of bath modes is essential for revealing their key roles in the overall dynamics. Here, we applied principal component analysis (PCA) to investigate the bath motion in the basis of a large dataset generated from the symmetrical quasi-classical dynamics method based on the Meyer-Miller mapping Hamiltonian nonadiabatic dynamics for the excited-state energy transfer in the Frenkel-exciton model. The PCA method clearly elucidated that two types of bath modes, which either display strong vibronic coupling or have frequencies close to that of the electronic transition, are important to the nonadiabatic dynamics. These observations were fully consistent with the physical insights. The conclusions were based on the PCA of the trajectory data and did not involve significant pre-defined physical knowledge. The results show that the PCA approach, which is one of the simplest unsupervised machine learning dimensionality reduction methods, is a powerful one for analyzing complicated nonadiabatic dynamics in the condensed phase with many degrees of freedom.
Collapse
Affiliation(s)
- Jiawei Peng
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| | - Yu Xie
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| | - Deping Hu
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
39
|
Brown SE, Shakib FA. Recent progress in approximate quantum dynamics methods for the study of proton-coupled electron transfer reactions. Phys Chem Chem Phys 2021; 23:2535-2556. [PMID: 33367437 DOI: 10.1039/d0cp05166g] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
Proton-coupled electron transfer (PCET) reactions are ubiquitous natural processes at the heart of energy conversion reactions in photosynthesis and respiration, DNA repair, and diverse enzymatic reactions. Theoretical formulation and computational method developments have eyed modeling of thermal and photoinduced PCET for the last three decades. The accumulation of these studies, collected in dozens of reviews, accounts, and perspectives, has firmly established the influence of quantum effects, including non-adiabatic electronic transitions, vibrational relaxation, zero-point energy, and proton tunneling, on the rate and mechanism of PCET reactions. Here, we focus on some recently-developed methods, spanning the last eight years, that can quantitatively capture these effects in the PCET context and provide efficient means for their qualitative description in complex systems. The theoretical background of each method and their accuracy with respect to exact results are discussed and the results of relevant PCET simulations based on each method are presented.
Collapse
Affiliation(s)
- Sandra E Brown
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, CA 90095, USA
| | - Farnaz A Shakib
- Department of Chemistry and Environmental Science, New Jersey Institute of Technology, Newark, NJ 07102, USA.
| |
Collapse
|
40
|
Tang D, Shen L, Fang WH. Evaluation of mixed quantum-classical molecular dynamics on cis-azobenzene photoisomerization. Phys Chem Chem Phys 2021; 23:13951-13964. [PMID: 34142685 DOI: 10.1039/d1cp01374b] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The quantitative prediction of nonadiabatic transitions between different electronic states is important to understand ultrafast processes in photochemistry. A variety of mixed quantum-classical molecular dynamics methods such as surface hopping and Ehrenfest mean-field have been developed. However, how to choose an appropriate one from a wide diversity of dynamics algorithms to study a realistic photochemical process is still unclear. In this work, we implemented 30 combinations of different mixed quantum-classical dynamics methods, including 24 surface hopping models with 8 decoherence corrections and 3 momentum rescaling strategies as well as 6 mean-field models. Then we performed numerical investigations by simulating the photoisomerization of cis-azobenzene combined with on-the-fly electronic structure calculations. Predictions of the S1 lifetime and the quantum yield of the photoproduct using different models are distinct. Surface hopping is more robust than mean-field in our test system. Moreover, the choice of momentum rescaling methods in surface hopping brings more significant changes than decoherence corrections, while a large discrepancy between simulation results with different mean-field algorithms has been observed.
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.
| | - Lin Shen
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China.
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China.
| |
Collapse
|
41
|
Gao X, Geva E. Improving the Accuracy of Quasiclassical Mapping Hamiltonian Methods by Treating the Window Function Width as an Adjustable Parameter. J Phys Chem A 2020; 124:11006-11016. [DOI: 10.1021/acs.jpca.0c09750] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Xing Gao
- School of Materials, Sun Yat-sen University, Shenzhen, Guangdong 518100, China
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
42
|
Mannouch JR, Richardson JO. A partially linearized spin-mapping approach for nonadiabatic dynamics. I. Derivation of the theory. J Chem Phys 2020; 153:194109. [DOI: 10.1063/5.0031168] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
43
|
Mannouch JR, Richardson JO. A partially linearized spin-mapping approach for nonadiabatic dynamics. II. Analysis and comparison with related approaches. J Chem Phys 2020; 153:194110. [DOI: 10.1063/5.0031173] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
44
|
Gao X, Lai Y, Geva E. Simulating Absorption Spectra of Multiexcitonic Systems via Quasiclassical Mapping Hamiltonian Methods. J Chem Theory Comput 2020; 16:6465-6480. [DOI: 10.1021/acs.jctc.0c00709] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Xing Gao
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Yifan Lai
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
45
|
Gao X, Geva E. A Nonperturbative Methodology for Simulating Multidimensional Spectra of Multiexcitonic Molecular Systems via Quasiclassical Mapping Hamiltonian Methods. J Chem Theory Comput 2020; 16:6491-6502. [DOI: 10.1021/acs.jctc.0c00843] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xing Gao
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
46
|
Braver Y, Valkunas L, Gelzinis A. Benchmarking the forward–backward trajectory solution of the quantum-classical Liouville equation. J Chem Phys 2020; 152:214116. [DOI: 10.1063/5.0006538] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yakov Braver
- Institute of Chemical Physics, Faculty of Physics, Vilnius University, Saulėtekio Ave. 9-III, LT-10222 Vilnius, Lithuania
- Department of Molecular Compound Physics, Center for Physical Sciences and Technology, Saulėtekio Ave. 3, LT-10257 Vilnius, Lithuania
| | - Leonas Valkunas
- Institute of Chemical Physics, Faculty of Physics, Vilnius University, Saulėtekio Ave. 9-III, LT-10222 Vilnius, Lithuania
- Department of Molecular Compound Physics, Center for Physical Sciences and Technology, Saulėtekio Ave. 3, LT-10257 Vilnius, Lithuania
| | - Andrius Gelzinis
- Institute of Chemical Physics, Faculty of Physics, Vilnius University, Saulėtekio Ave. 9-III, LT-10222 Vilnius, Lithuania
- Department of Molecular Compound Physics, Center for Physical Sciences and Technology, Saulėtekio Ave. 3, LT-10257 Vilnius, Lithuania
| |
Collapse
|
47
|
Liu Y, Gao X, Lai Y, Mulvihill E, Geva E. Electronic Dynamics through Conical Intersections via Quasiclassical Mapping Hamiltonian Methods. J Chem Theory Comput 2020; 16:4479-4488. [DOI: 10.1021/acs.jctc.0c00177] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Yudan Liu
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Xing Gao
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Yifan Lai
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Ellen Mulvihill
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
48
|
Gao X, Saller MAC, Liu Y, Kelly A, Richardson JO, Geva E. Benchmarking Quasiclassical Mapping Hamiltonian Methods for Simulating Electronically Nonadiabatic Molecular Dynamics. J Chem Theory Comput 2020; 16:2883-2895. [DOI: 10.1021/acs.jctc.9b01267] [Citation(s) in RCA: 37] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Xing Gao
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | | | - Yudan Liu
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Aaron Kelly
- Department of Chemistry, Dalhousie University, 15000 Halifax, Nova Scotia, Canada
| | | | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
49
|
Ranya S, Ananth N. Multistate ring polymer instantons and nonadiabatic reaction rates. J Chem Phys 2020; 152:114112. [DOI: 10.1063/1.5132807] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Srinath Ranya
- 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
|
50
|
Mandal A, Hunt KLC. Variance of the energy of a quantum system in a time-dependent perturbation: Determination by nonadiabatic transition probabilities. J Chem Phys 2020; 152:104110. [DOI: 10.1063/1.5140009] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Affiliation(s)
- Anirban Mandal
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Katharine L. C. Hunt
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|