Practical and reliable approximation of nonadiabatic coupling terms between triplet electronic states using only adiabatic potential energies

Dynamics of FHCl Including Two Bidirectional Dissociation Channels: Comparative Study Using Quantum Nuclear Wavepackets and Semiclassical Trajectory Surface Hopping

The dynamics starting from the charge transfer excited state of neutral FHCl includes two bidirectional dissociation channels, producing "FH" and "Cl" fragments or "F" and "HCl" fragments by electron-transfer (ET) or proton-transfer (PT) processes, respectively. The quantum features of its dynamics were studied using the quantum dynamics of nuclear wavepacket propagation and the semiclassical dynamics of trajectory surface hopping propagation. The MS-CASPT2(17,11)/aug-cc-pVTZ method was used for calculating the energies of electronic states. Two critical quantum features identified in previous studies were the dominance of PT over ET and the ∼80 fs time gap between the onset of the earlier PT and the later ET processes. These features, in contrast to classical anticipation, were confirmed again, and their classical interpretations were developed based on the results of semiclassical dynamics. The relative location of nonadiabatic regions with respect to the starting point of dynamics and the activation of angular motion altering the HFCl angle play crucial roles in governing the two quantum features.

Nonadiabatic Derivative Couplings Calculated Using Information of Potential Energy Surfaces without Wavefunctions: Ab Initio and Machine Learning Implementations

In this work, we implemented an approximate algorithm for calculating nonadiabatic coupling matrix elements (NACMEs) of a polyatomic system with ab initio methods and machine learning (ML) models. Utilizing this algorithm, one can calculate NACMEs using only the information of potential energy surfaces (PESs), i.e., energies, and gradients as well as Hessian matrix elements. We used a realistic system, namely CH₂NH, to compare NACMEs calculated by this approximate PES-based algorithm and the accurate wavefunction-based algorithm. Our results show that this approximate PES-based algorithm can give very accurate results comparable to the wavefunction-based algorithm except at energetically degenerate points, i.e., conical intersections. We also tested a machine learning (ML)-trained model with this approximate PES-based algorithm, which also supplied similarly accurate NACMEs but more efficiently. The advantage of this PES-based algorithm is its significant potential to combine with electronic structure methods that do not implement wavefunction-based algorithms, low-scaling energy-based fragment methods, etc., and in particular efficient ML models, to compute NACMEs. The present work could encourage further research on nonadiabatic processes of large systems simulated by ab initio nonadiabatic dynamics simulation methods in which NACMEs are always required.

Nuclear wave-packet-propagation-based study of the electron-coupled, proton-transfer process in the charge-transfer state of FHCl exhibiting three electronic states in full-dimensional space

The charge-transfer (CT) excited state of FHCl (F⁺H-Cl^-), generated by the photodetachment of an electron from its precursor anion (FHCl^-) by a photon energy of ∼9.5 eV, is a realistic prototype of two bidirectional-coupled reaction pathways, namely the proton-transfer (PT) and electron-transfer (ET) channels, that produce F + HCl and FH + Cl combinations, respectively. The early-time dynamics of the CT was studied via the time-dependent propagations of nuclear wave packets comprising three nonadiabatically coupled electronic states defined within a three-dimensional space. The detailed analyses of the early-time dynamics revealed an interesting phenomenon in which the onset of PT was ∼80 fs earlier than that of ET, indicating that PT dominated ET in this case. A more significant finding was that the proper adjustment of the electronic-charge distribution for the onset of ET was obtained ∼80 fs after the onset of PT; this adjustment was mediated by the initial movement of the H atom, i.e., the F-H vibration mode. To avail experimental observables, the branching ratio, χ = PT/(PT + ET), and absorption spectrum generating the neutral FHCl molecule from its precursor anion were also simulated. The results further demonstrated the dependences of the χs and spectrum on the change in the initial vibration level of the precursor anion, as well as the isotopic substitution of the connecting H atom with deuterium, tritium, and muonium.

Fewest switches surface hopping with Baeck-An couplings

In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.

Fewest switches surface hopping with Baeck-An couplings

In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.

Wavepacket propagations for the early time dynamics of proton-coupled electron transfer in the charge-transfer state of NH3Cl complex

A charge-transfer (CT) excited state of NH₃Cl, generated by photo-detachment of an electron from the anionic NH₃Cl^- precursor, can be represented as H₂N⁺-H-Cl^- and proceeds to two chemical reactions: one reaction generating NH₂ and HCl resulting from a proton transfer (PT) and the other reaction producing NH₃ and a Cl atom resulting from an electron transfer (ET); both are coupled to form a typical proton-coupled electron transfer (PCET) process. The early time dynamics of this CT were studied using time-dependent wavepacket propagation on three nonadiabatically coupled electronic states in a reduced three-dimensional space. The electronic states were treated using the XMS-CASPT2/aug-cc-pVTZ ab initio methodology. The population dynamics of the three coupled electronic states were analyzed in detail to reveal the initial stage of the PCET process up to ∼100 fs, while the branching ratio, χ = PT/(ET+PT), was determined after wavepacket propagations of up to 2000 fs. Another main result is the dependence of χ on the vibration levels of the initial precursor anion and the isotope substitution of the connecting H atom with deuterium and tritium. Our study reveals the detailed microscopic features of the PCET process embedded in the CT state of the NH₃Cl complex and certain systematic dependences of the branching ratio χ on the above factors.

Combining SchNet and SHARC: The SchNarc Machine Learning Approach for Excited-State Dynamics

In recent years, deep learning has become a part of our everyday life and is revolutionizing quantum chemistry as well. In this work, we show how deep learning can be used to advance the research field of photochemistry by learning all important properties-multiple energies, forces, and different couplings-for photodynamics simulations. We simplify such simulations substantially by (i) a phase-free training skipping costly preprocessing of raw quantum chemistry data; (ii) rotationally covariant nonadiabatic couplings, which can either be trained or (iii) alternatively be approximated from only ML potentials, their gradients, and Hessians; and (iv) incorporating spin-orbit couplings. As the deep-learning method, we employ SchNet with its automatically determined representation of molecular structures and extend it for multiple electronic states. In combination with the molecular dynamics program SHARC, our approach termed SchNarc is tested on two polyatomic molecules and paves the way toward efficient photodynamics simulations of complex systems.