Time-Derivative Couplings for Self-Consistent Electronically Nonadiabatic Dynamics

Generalized Semiclassical Ehrenfest Method: A Route to Wave Function-Free Photochemistry and Nonadiabatic Dynamics with Only Potential Energies and Gradients

We reconsider recent methods by which direct dynamics calculations of electronically nonadiabatic processes can be carried out while requiring only adiabatic potential energies and their gradients. We show that these methods can be understood in terms of a new generalization of the well-known semiclassical Ehrenfest method. This is convenient because it eliminates the need to evaluate electronic wave functions and their matrix elements along the mixed quantum-classical trajectories. The new approximations and procedures enabling this advance are the curvature-driven approximation to the time-derivative coupling, the generalized semiclassical Ehrenfest method, and a new gradient correction scheme called the time-derivative matrix (TDM) scheme. When spin-orbit coupling is present, one can carry out dynamics calculations in the fully adiabatic basis using potential energies and gradients calculated without spin-orbit coupling plus the spin-orbit coupling matrix elements. Even when spin-orbit coupling is neglected, the method is useful because it allows calculations by electronic structure methods for which nonadiabatic coupling vectors are unavailable. In order to place the new considerations in context, the article starts out with a review of background material on trajectory surface hopping, the semiclassical Ehrenfest scheme, and methods for incorporating decoherence. We consider both internal conversion and intersystem crossing. We also review several examples from our group of successful applications of the curvature-driven approximation.

Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

Nonadiabatic Coupling in Trajectory Surface Hopping: Accurate Time Derivative Couplings by the Curvature-Driven Approximation

Trajectory surface hopping (TSH) is a widely used mixed quantum-classical dynamics method that is used to simulate molecular dynamics with multiple electronic states. In TSH, time-derivative coupling is employed to propagate the electronic coefficients and in that way to determine when the electronic state on which the nuclear trajectory is propagated switches. In this work, we discuss nonadiabatic TSH dynamics algorithms employing the curvature-driven approximation and overlap-based time derivative couplings, and we report test calculations on six photochemical reactions where we compare the results to one another and to calculations employing analytic nonadiabatic coupling vectors. We correct previous published results thanks to a bug found in the software. We also provide additional, more detailed studies of the time-derivative couplings. Our results show good agreement between curvature-driven algorithms and overlap-based algorithms.

Semiclassical Multistate Dynamics for Six Coupled 5A' States of O + O2

Dynamics simulations of high-energy O₂-O collisions play an important role in simulating thermal energy content and heat flux in flows around hypersonic vehicles. To carry out such dynamics simulations efficiently requires accurate global potential energy surfaces and (in most algorithms) state couplings for many energetically accessible electronic states. The ability to treat collisions involving many coupled electronic states has been a challenge for decades. Very recently, a new diabatization method, the parametrically managed diabatization by deep neural network (PM-DDNN), has been developed. The PM-DDNN method uses a deep neural network architecture with an activation function parametrically dependent on input data to discover and fit the diabatic potential energy matrix (DPEM) as a function of geometry, and the adiabatic potential energy surfaces are obtained by diagonalization of a small matrix with analytic matrix elements. Here, we applied the PM-DDNN method to the six lowest-energy potential energy surfaces in the ⁵A' manifold of O₃ to perform simultaneous diabatization and fitting; the data are obtained by extended multistate complete-active-space second-order perturbation theory. We then used the adiabatic surfaces for dynamics calculations with three methods: coherent switching with decay of mixing (CSDM), curvature-driven CSDM (κCSDM), and electronically curvature-driven CSDM (eκCSDM). The κCSDM calculations require only adiabatic potential energies and gradients. The three dynamical methods are in good agreement. We then calculated electronically nonadiabatic, electronically inelastic, and dissociative cross sections for seven initial collision energies, five initial vibrational levels, and four initial rotational levels. Trends in the electronically inelastic cross sections as functions of the initial collision energy and vibrational level were rationalized in terms of the coordinate ranges where the gaps between the second and third potential energy surfaces are small.

Direct Nonadiabatic Dynamics of Ammonia with Curvature-Driven Coherent Switching with Decay of Mixing and with Fewest Switches with Time Uncertainty: An Illustration of Population Leaking in Trajectory Surface Hopping Due to Frustrated Hops

Mixed quantum-classical nonadiabatic dynamics is a widely used approach to simulate molecular dynamics involving multiple electronic states. There are two main categories of mixed quantum-classical nonadiabatic dynamics algorithms, namely, trajectory surface hopping (TSH) in which the trajectory propagates on a single potential energy surface, interrupted by hops, and self-consistent-potential (SCP) methods, such as semiclassical Ehrenfest, in which propagation occurs on a mean-field surface without hops. In this work, we will illustrate an example of severe population leaking in TSH. We emphasize that such leaking is a combined effect of frustrated hops and long-time simulations that drive the final excited-state population toward zero as a function of time. We further show that such leaking can be alleviated-but not eliminated-by the fewest switches with time uncertainty TSH algorithm (here implemented in the SHARC program); the time uncertainty algorithm slows down the leaking process by a factor of 4.1. The population leaking is not present in coherent switching with decay of mixing (CSDM), which is an SCP method with non-Markovian decoherence included. Another result in this paper is that we find very similar results with the original CSDM algorithm, with time-derivative CSDM (tCSDM), and with curvature-driven CSDM (κCSDM). Not only do we find good agreement for electronically nonadiabatic transition probabilities but also we find good agreement of the norms of the effective nonadiabatic couplings (NACs) that are derived from the curvature-driven time-derivative couplings as implemented in κCSDM with the time-dependent norms of the nonadiabatic coupling vectors computed by state-averaged complete-active-space self-consistent field theory.

Decoherence and Its Role in Electronically Nonadiabatic Dynamics

Decoherence is the tendency of a time-evolved reduced density matrix for a subsystem to assume a form corresponding to a statistical ensemble of states rather than a coherent combination of pure-state wave functions. When a molecular process involves changes in the electronic state and the coordinates of the nuclei, as in ultraviolet or visible light photochemistry or electronically inelastic collisions, the reduced density matrix of the electronic subsystem suffers decoherence, due to its interaction with the nuclear subsystem. We present the background necessary to conceptualize this decoherence; in particular, we discuss the density matrix description of pure states and mixed states, and we discuss pointer states and decoherence time. We then discuss how decoherence is treated in the coherent switching with decay of mixing algorithm and the trajectory surface hopping method for semiclassical calculations of electronically nonadiabatic processes.

Nonadiabatic Dynamics of 1,3-Cyclohexadiene by Curvature-Driven Coherent Switching with Decay of Mixing

The photoinduced ring-opening reaction of 1,3-cyclohexadiene to produce 1,3,5-hexatriene is a classic electrocyclic reaction and is also a prototype for many reactions of biological and synthetic importance. Here, we simulate the ultrafast nonadiabatic dynamics of the reaction in the manifold of the three lowest valence electronic states by using extended multistate complete-active-space second-order perturbation theory (XMS-CASPT2) combined with the curvature-driven coherent switching with decay of mixing (κCSDM) dynamical method. We obtain an excited-state lifetime of 79 fs, and a product quantum yield of 40% from the 500 trajectories initiated in the S₁ excited state. The obtained lifetime and quantum yield values are very close to previously reported experimental and computed values, showing the capability of performing a reasonable nonadiabatic ring-opening dynamics with the κCSDM method that does not require nonadiabatic coupling vectors, time derivatives, or diabatization. In addition, we study the ring-opening reaction by initiating the trajectories in the dark state S₂. We also optimize the S₀/S₁ and S₁/S₂ minimum-energy conical intersections (MECIs) by XMS-CASPT2; for S₁/S₂, we optimized both an inner and an outer local-minimum-energy conical intersections (LMECIs). We provide the potential energy profile along the ring-opening coordinate by joining selected critical points via linear synchronous transit paths. We find the inner S₁/S₂ LMECI to be more crucial than the outer one.

Diabatic States of Molecules

Quantitative simulations of electronically nonadiabatic molecular processes require both accurate dynamics algorithms and accurate electronic structure information. Direct semiclassical nonadiabatic dynamics is expensive due to the high cost of electronic structure calculations, and hence it is limited to small systems, limited ensemble averaging, ultrafast processes, and/or electronic structure methods that are only semiquantitatively accurate. The cost of dynamics calculations can be made manageable if analytic fits are made to the electronic structure data, and such fits are most conveniently carried out in a diabatic representation because the surfaces are smooth and the couplings between states are smooth scalar functions. Diabatic representations, unlike the adiabatic ones produced by most electronic structure methods, are not unique, and finding suitable diabatic representations often involves time-consuming nonsystematic diabatization steps. The biggest drawback of using diabatic bases is that it can require large amounts of effort to perform a globally consistent diabatization, and one of our goals has been to develop methods to do this efficiently and automatically. In this Feature Article, we introduce the mathematical framework of diabatic representations, and we discuss diabatization methods, including adiabatic-to-diabatic transformations and recent progress toward the goal of automatization.

Nonadiabatic Dynamics Algorithms with Only Potential Energies and Gradients: Curvature-Driven Coherent Switching with Decay of Mixing and Curvature-Driven Trajectory Surface Hopping

Direct dynamics by mixed quantum-classical nonadiabatic methods is an important tool for understanding processes involving multiple electronic states. Very often, the computational bottleneck of such direct simulation comes from electronic structure theory. For example, at every time step of a trajectory, nonadiabatic dynamics requires potential energy surfaces, their gradients, and the matrix elements coupling the surfaces. The need for the couplings can be alleviated by employing the time derivatives of the wave functions, which can be evaluated from overlaps of electronic wave functions at successive time steps. However, evaluation of overlap integrals is still expensive for large systems. In addition, for electronic structure methods for which the wave functions or the coupling matrix elements are not available, nonadiabatic dynamics algorithms become inapplicable. In this work, building on recent work by Baeck and An, we propose new nonadiabatic dynamics algorithms that only require adiabatic potential energies and their gradients. The new methods are named curvature-driven coherent switching with decay of mixing (κCSDM) and curvature-driven trajectory surface hopping (κTSH). We show how powerful these new methods are in terms of computation time and accuracy as compared to previous mixed quantum-classical nonadiabatic dynamics algorithms. The lowering of the computational cost will allow longer nonadiabatic trajectories and greater ensemble averaging to be affordable, and the ability to calculate the dynamics without electronic structure coupling matrix elements extends the dynamics capability to new classes of electronic structure methods.

Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections

Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.

Direct coherent switching with decay of mixing for intersystem crossing dynamics of thioformaldehyde: The effect of decoherence

We evaluate the effect of electronic decoherence on intersystem crossing in the photodynamics of thioformaldehyde. First, we show that the state-averaged complete-active-space self-consistent field electronic structure calculations with a properly chosen active space of 12 active electrons in 10 active orbitals can predict the potential energy surfaces and the singlet-triplet spin-orbit couplings quite well for CH₂S, and we use this method for direct dynamics by coherent switching with decay of mixing (CSDM). We obtain similar dynamical results with CSDM or by adding energy-based decoherence to trajectory surface hopping, with the population of triplet states tending to a small steady-state value over 500 fs. Without decoherence, the state populations calculated by the conventional trajectory surface hopping method or the semiclassical Ehrenfest method gradually increase. This difference shows that decoherence changes the nature of the results not just quantitatively but qualitatively.