1
|
Curchod BFE, Glover WJ, Martínez TJ. SSAIMS-Stochastic-Selection Ab Initio Multiple Spawning for Efficient Nonadiabatic Molecular Dynamics. J Phys Chem A 2020; 124:6133-6143. [PMID: 32580552 DOI: 10.1021/acs.jpca.0c04113] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Ab initio multiple spawning provides a powerful and accurate way of describing the excited-state dynamics of molecular systems, whose strength resides in the proper description of coherence effects during nonadiabatic processes thanks to the coupling of trajectory basis functions. However, the simultaneous propagation of a large number of trajectory basis functions can be numerically inconvenient. We propose here an elegant and simple solution to this issue, which consists of (i) detecting uncoupled groups of coupled trajectory basis functions and (ii) selecting stochastically one of these groups to continue the ab initio multiple spawning dynamics. We show that this procedure can reproduce the results of full ab initio multiple spawning dynamics in cases where the uncoupled groups of trajectory basis functions stay uncoupled throughout the dynamics (which is often the case in high-dimensional problems). We present and discuss the aforementioned idea in detail and provide simple numerical applications on indole, ethylene, and protonated formaldimine, highlighting the potential of stochastic-selection ab initio multiple spawning.
Collapse
Affiliation(s)
- Basile F E Curchod
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - William J Glover
- NYU Shanghai, 1555 Century Ave., Shanghai 200122, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China.,Department of Chemistry, New York University, New York, New York 10003, United States
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| |
Collapse
|
2
|
Mališ M, Luber S. Trajectory Surface Hopping Nonadiabatic Molecular Dynamics with Kohn–Sham ΔSCF for Condensed-Phase Systems. J Chem Theory Comput 2020; 16:4071-4086. [DOI: 10.1021/acs.jctc.0c00372] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Momir Mališ
- University of Zurich, Department of Chemistry, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Sandra Luber
- University of Zurich, Department of Chemistry, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| |
Collapse
|
3
|
Nelson TR, White AJ, Bjorgaard JA, Sifain AE, Zhang Y, Nebgen B, Fernandez-Alberti S, Mozyrsky D, Roitberg AE, Tretiak S. Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials. Chem Rev 2020; 120:2215-2287. [PMID: 32040312 DOI: 10.1021/acs.chemrev.9b00447] [Citation(s) in RCA: 212] [Impact Index Per Article: 53.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.
Collapse
Affiliation(s)
- Tammie R Nelson
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Alexander J White
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Josiah A Bjorgaard
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Andrew E Sifain
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States.,U.S. Army Research Laboratory , Aberdeen Proving Ground , Maryland 21005 , United States
| | - Yu Zhang
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Benjamin Nebgen
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | | | - Dmitry Mozyrsky
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Adrian E Roitberg
- Department of Chemistry , University of Florida , Gainesville , Florida 32611 , United States
| | - Sergei Tretiak
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| |
Collapse
|
4
|
Wang L, Qiu J, Bai X, Xu J. Surface hopping methods for nonadiabatic dynamics in extended systems. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1435] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Linjun Wang
- Center for Chemistry of Novel & High‐Performance Materials, Department of Chemistry Zhejiang University Hangzhou China
| | - Jing Qiu
- Center for Chemistry of Novel & High‐Performance Materials, Department of Chemistry Zhejiang University Hangzhou China
| | - Xin Bai
- Center for Chemistry of Novel & High‐Performance Materials, Department of Chemistry Zhejiang University Hangzhou China
| | - Jiabo Xu
- Center for Chemistry of Novel & High‐Performance Materials, Department of Chemistry Zhejiang University Hangzhou China
| |
Collapse
|
5
|
Abstract
We present a subspace surface hopping strategy to deal with complex surface crossings in nonadiabatic dynamics. By focusing on only important adiabatic states, we make subspace crossing correction (SCC) in the framework of the standard fewest switches surface hopping (FSSH) and the global flux surface hopping (GFSH). The resulting SCC-FSSH and SCC-GFSH approaches show much better performance than the counterparts using all adiabatic states for surface hopping. As demonstrated in a series of Holstein models with up to over 1000 molecular sites, both SCC-FSSH and SCC-GFSH show excellent size independence with a large time step size of 1 fs. Especially, SCC-GFSH does not refer to nonadiabatic couplings at all and gives a more proper description of superexchange, and thus, it is promising for realistic applications with complex potential energy surfaces.
Collapse
Affiliation(s)
- Jing Qiu
- Center for Chemistry of Novel & High-Performance Materials , and Department of Chemistry , Zhejiang University , Hangzhou 310027 , China
| | - Xin Bai
- Center for Chemistry of Novel & High-Performance Materials , and Department of Chemistry , Zhejiang University , Hangzhou 310027 , China
| | - Linjun Wang
- Center for Chemistry of Novel & High-Performance Materials , and Department of Chemistry , Zhejiang University , Hangzhou 310027 , China
| |
Collapse
|
6
|
Qiu J, Bai X, Wang L. Crossing Classified and Corrected Fewest Switches Surface Hopping. J Phys Chem Lett 2018; 9:4319-4325. [PMID: 30011207 DOI: 10.1021/acs.jpclett.8b01902] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
In the traditional fewest switches surface hopping (FSSH), trivial crossings between uncoupled or weakly coupled states have highly peaked nonadiabatic couplings and thus are difficult to deal with in the preferred, adiabatic representation. Here, we classify surface crossings into four general types and propose a parameter-free crossing corrected FSSH (CC-FSSH) algorithm, which could treat multiple trivial crossings within a time interval. As examples, Holstein Hamiltonians with different parameters are adopted to mimic electron dynamics in tens to hundreds of molecules, which suffer from severe trivial crossing problems. Using existed surface hopping approaches as references, we show that CC-FSSH exhibits significantly fast time interval convergence and weak system size dependence. In all cases, a reliable description is achieved with a large time interval of 1 fs. With a simple formalism and the ability to describe complex surface crossings, CC-FSSH could potentially simulate general nonadiabatic dynamics in nanoscale materials with a high efficiency.
Collapse
Affiliation(s)
- Jing Qiu
- Department of Chemistry , Zhejiang University , Hangzhou 310027 , China
| | - Xin Bai
- Department of Chemistry , Zhejiang University , Hangzhou 310027 , China
| | - Linjun Wang
- Department of Chemistry , Zhejiang University , Hangzhou 310027 , China
| |
Collapse
|
7
|
Affiliation(s)
- Michael S. Schuurman
- National Research Council of Canada, Ottawa, Ontario K1A 06A, Canada
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1A 0R6, Canada
| | - Albert Stolow
- National Research Council of Canada, Ottawa, Ontario K1A 06A, Canada
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1A 0R6, Canada
- Department of Physics, University of Ottawa, Ottawa, Ontario K1A 0R6, Canada
| |
Collapse
|
8
|
Affiliation(s)
- Basile F. E. Curchod
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Todd J. Martínez
- Department of Chemistry and PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| |
Collapse
|
9
|
Ichikawa H, Takatsuka K. Chemical Modification of Conical Intersections in Photoisomerization Dynamics of Butadiene Derivatives. J Phys Chem A 2016; 121:315-325. [DOI: 10.1021/acs.jpca.6b10680] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Hiroki Ichikawa
- Department of Basic Science, The University of Tokyo, Komaba 3-8-1, Tokyo 153-8902, Japan
| | - Kazuo Takatsuka
- Fukui Institute for Fundamental
Chemistry, Kyoto University, Takano-Nishikiraki-cho 34-4, Sakyo-ku, Kyoto 606-8103, Japan
| |
Collapse
|
10
|
Meek GA, Levine BG. The best of both Reps—Diabatized Gaussians on adiabatic surfaces. J Chem Phys 2016; 145:184103. [DOI: 10.1063/1.4966967] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Affiliation(s)
- Garrett A. Meek
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Benjamin G. Levine
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
11
|
Gamiz-Hernandez AP, Kaila VRI. Conversion of light-energy into molecular strain in the photocycle of the photoactive yellow protein. Phys Chem Chem Phys 2016; 18:2802-9. [DOI: 10.1039/c5cp05244k] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
The Photoactive Yellow Protein (PYP) converts light energy into molecular strain, stored in the early pR0-photocycle intermediate.
Collapse
Affiliation(s)
| | - Ville R. I. Kaila
- Department Chemie
- Technische Universität München (TUM)
- D-85747 Garching
- Germany
| |
Collapse
|
12
|
Abstract
This Perspective addresses the use of coupled diabatic potential energy surfaces (PESs) together with rigorous quantum dynamics in full or reduced dimensional coordinate spaces to obtain accurate solutions to problems in nonadiabatic dynamics.
Collapse
Affiliation(s)
- Hua Guo
- Department of Chemistry and Chemical Biology
- University of New Mexico
- Albuquerque
- USA
| | | |
Collapse
|
13
|
Xia SH, Liu XY, Fang Q, Cui G. Photodissociation dynamics of CH3C(O)SH in argon matrix: A QM/MM nonadiabatic dynamics simulation. J Chem Phys 2015; 143:194303. [DOI: 10.1063/1.4935598] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
|
14
|
Meek GA, Levine BG. Accurate and efficient evaluation of transition probabilities at unavoided crossings in ab initio multiple spawning. Chem Phys 2015. [DOI: 10.1016/j.chemphys.2015.06.007] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
15
|
Kubař T, Welke K, Groenhof G. New QM/MM implementation of the DFTB3 method in the gromacs package. J Comput Chem 2015; 36:1978-89. [PMID: 26238364 DOI: 10.1002/jcc.24029] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2015] [Accepted: 07/09/2015] [Indexed: 01/07/2023]
Abstract
The approximate density-functional tight-binding theory method DFTB3 has been implemented in the quantum mechanics/molecular mechanics (QM/MM) framework of the Gromacs molecular simulation package. We show that the efficient smooth particle-mesh Ewald implementation of Gromacs extends to the calculation of QM/MM electrostatic interactions. Further, we make use of the various free-energy functionalities provided by Gromacs and the PLUMED plugin. We exploit the versatility and performance of the current framework in three typical applications of QM/MM methods to solve biophysical problems: (i) ultrafast proton transfer in malonaldehyde, (ii) conformation of the alanine dipeptide, and (iii) electron-induced repair of a DNA lesion. Also discussed is the further development of the framework, regarding mostly the options for parallelization.
Collapse
Affiliation(s)
- Tomáš Kubař
- Institute of Physical Chemistry and Center for Functional Nanostructures, Karlsruhe Institute of Technology, Karlsruhe, 76131, Germany
| | - Kai Welke
- Department of Chemistry, Nagoya University, Nagoya, 464-8602, Japan
| | - Gerrit Groenhof
- Nanoscience Center and Department of Chemistry, University of Jyväskylä, Jyväskylä, 40014, Finland
| |
Collapse
|
16
|
Glover WJ. Communication: Smoothing out excited-state dynamics: analytical gradients for dynamically weighted complete active space self-consistent field. J Chem Phys 2015; 141:171102. [PMID: 25381494 DOI: 10.1063/1.4901328] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
State averaged complete active space self-consistent field (SA-CASSCF) is a workhorse for determining the excited-state electronic structure of molecules, particularly for states with multireference character; however, the method suffers from known issues that have prevented its wider adoption. One issue is the presence of discontinuities in potential energy surfaces when a state that is not included in the state averaging crosses with one that is. In this communication I introduce a new dynamical weight with spline (DWS) scheme that mimics SA-CASSCF while removing energy discontinuities due to unweighted state crossings. In addition, analytical gradients for DWS-CASSCF (and other dynamically weighted schemes) are derived for the first time, enabling energy-conserving excited-state ab initio molecular dynamics in instances where SA-CASSCF fails.
Collapse
Affiliation(s)
- W J Glover
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| |
Collapse
|
17
|
Meek GA, Levine BG. Evaluation of the Time-Derivative Coupling for Accurate Electronic State Transition Probabilities from Numerical Simulations. J Phys Chem Lett 2014; 5:2351-2356. [PMID: 26279558 DOI: 10.1021/jz5009449] [Citation(s) in RCA: 79] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Spikes in the time-derivative coupling (TDC) near surface crossings make the accurate integration of the time-dependent Schrödinger equation in nonadiabatic molecular dynamics simulations a challenge. To address this issue, we present an approximation to the TDC based on a norm-preserving interpolation (NPI) of the adiabatic electronic wave functions within each time step. We apply NPI and two other schemes for computing the TDC in numerical simulations of the Landau-Zener model, comparing the simulated transfer probabilities to the exact solution. Though NPI does not require the analytical calculation of nonadiabatic coupling matrix elements, it consistently yields unsigned population transfer probability errors of ∼0.001, whereas analytical calculation of the TDC yields errors of 0.0-1.0 depending on the time step, the offset of the maximum in the TDC from the beginning of the time step, and the coupling strength. The approximation of Hammes-Schiffer and Tully yields errors intermediate between NPI and the analytical scheme.
Collapse
Affiliation(s)
- Garrett A Meek
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Benjamin G Levine
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|