Incorporating Berry Force Effects into the Fewest Switches Surface-Hopping Algorithm: Intersystem Crossing and the Case of Electronic Degeneracy

A Phase-Space Electronic Hamiltonian for Molecules in a Static Magnetic Field II: Quantum Chemistry Calculations with Gauge Invariant Atomic Orbitals

In a companion paper, we have developed a phase-space electronic structure theory of molecules in magnetic fields, whereby the electronic energy levels arise from diagonalizing a phase-space Hamiltonian ĤPS(X, P, G, B) that depends parametrically on nuclear position and momentum. The resulting eigenvalues are translationally invariant; moreover, if the magnetic field is in the z-direction, then the eigenvalues are also invariant to rotations around the z-direction. However, like all Hamiltonians in a magnetic field, the theory has a gauge degree of freedom (corresponding to the position of the magnetic origin in the vector potential), and requires either (i) formally, a complete set of electronic states or (ii) in practice, gauge-invariant atomic orbitals (GIAOs) in order to realize such translational and rotational invariance. Here we describe how to implement a phase-space electronic Hamiltonian using GIAOs within a practical electronic structure package (in our case, Q-Chem). We further show that novel phenomena can be observed with finite B-fields, including minimum energy structures with Πmin ≠ 0, indicating nonzero electronic motion in the ground-state.

Generalization of Quantum-Trajectory Surface Hopping to Multiple Quantum States

In this work, we present a generalization of the quantum trajectory surface hopping (QTSH) to multiple states and its implementation in the Libra package for nonadiabatic dynamics. In lieu of the ad hoc velocity rescaling used in many trajectory-based surface hopping approaches, QTSH utilizes quantum forces to evolve nuclear degrees of freedom continuously. It also lifts the unphysical constraint of enforcing the total energy conservation at the individual trajectory level and rather conserves the total energy at the trajectory ensemble level. Leveraging our new implementation of the multistate QTSH, we perform a comparative analysis of this method with the conventional fewest switches surface hopping approach. We combine the QTSH and decoherence corrections based on the simplified decay of mixing (SDM) and exact factorization (XF), leading to the QTSH-SDM and QTSH-XF schemes. Using the Holstein, superexchange, and phenol model Hamiltonians, we assess the relative accuracy of the resulting combined schemes in reproducing branching ratios, population, and coherence dynamics for a broad range of initial conditions. We observe that the decoherence correction in QTSH is crucial to improve energy conservation as well as the internal consistency between the population from the quantum probability and active state.

Mixed Quantum-Classical Dynamics under Arbitrary Unitary Basis Transformations

A common approach to minimizing the cost of quantum computations is by unitarily transforming a quantum system into a basis that can be optimally truncated. Here, we derive classical equations of motion subjected to similar unitary transformations and propose their integration into mixed quantum-classical dynamics, allowing this class of methods to be applied within arbitrary bases for both the quantum and classical coordinates. To this end, canonical positions and momenta of the classical degrees of freedom are combined into a set of complex-valued coordinates amenable to unitary transformations. We demonstrate the potential of the resulting approach by means of surface hopping calculations of an electronic carrier scattering onto a single impurity in the presence of phonons. Appropriate basis transformations, capturing both the localization of the impurity and the delocalization of higher-energy excitations, are shown to faithfully capture the dynamics within a fraction of the classical and quantum basis sets.

Spin-Dependent Stereochemistry: A Nonadiabatic Quantum Dynamics Case Study of S + H2 → SH + H Reaction

We study the spin-dependent stereodynamics of the S + H2 → SH + H reaction by using full-dimensional quantum dynamics calculations with zero total nuclear angular momentum along the triplet 3A″ states and singlet 1A' states. We find that the interplay between the electronic spin direction and the molecular geometry has a measurable influence on the singlet-triplet intersystem crossing reaction probabilities. Our results show that for some incident scattering angles in the body-fixed frame, the relative difference in intersystem crossing reaction probabilities (as determined between spin up and spin down initial states) can be as large as 15%. Our findings are an ab initio demonstration of spin-dependent nonadiabatic dynamics, which we hope will shine light as far as understanding the chiral-induced spin selectivity effect.

On the circularly polarized luminescence of individual triplet sublevels

We discuss the possibility of using circularly polarized luminescence (CPL) as a tool to probe individual triplet spin sublevels that are populated nonadiabatically following photoexcitation. This study is motivated by a mechanism proposed for chirality-induced spin selectivity in which coupled electronic-nuclear dynamics may lead to a non-statistical population of the three triplet sublevels in chiral systems. We find that low-temperature CPL should aid in quantifying the exact spin state/s populated through coupled electronic-nuclear motion in chiral molecules.

Nonadiabatic Dynamics in a Continuous Circularly Polarized Laser Field with Floquet Phase-Space Surface Hopping

Nonadiabatic chemical reactions involving continuous circularly polarized light (cw CPL) have not attracted as much attention as dynamics in unpolarized/linearly polarized light. However, including circularly (in contrast to linearly) polarized light allows one to effectively introduce a complex-valued time-dependent Hamiltonian, which offers a new path for control or exploration through the introduction of Berry forces. Here, we investigate several inexpensive semiclassical approaches for modeling such nonadiabatic dynamics in the presence of a time-dependent complex-valued Hamiltonian, beginning with a straightforward instantaneous adiabatic fewest-switches surface hopping (IA-FSSH) approach (where the electronic states depend on position and time), continuing to a standard Floquet fewest switches surface hopping (F-FSSH) approach (where the electronic states depend on position and frequency), and ending with an exotic Floquet phase-space surface hopping (F-PSSH) approach (where the electronic states depend on position, frequency, and momentum). Using a set of model systems with time-dependent complex-valued Hamiltonians, we show that the Floquet phase-space adiabats are the optimal choice of basis as far as accounting for Berry phase effects and delivering accuracy. Thus, the F-PSSH algorithm sets the stage for future modeling of nonadiabatic dynamics under strong externally pumped circular polarization.

A quantum-classical Liouville formalism in a preconditioned basis and its connection with phase-space surface hopping

We revisit a recent proposal to model nonadiabatic problems with a complex-valued Hamiltonian through a phase-space surface hopping (PSSH) algorithm employing a pseudo-diabatic basis. Here, we show that such a pseudo-diabatic PSSH (PD-PSSH) ansatz is consistent with a quantum-classical Liouville equation (QCLE) that can be derived following a preconditioning process, and we demonstrate that a proper PD-PSSH algorithm is able to capture some geometric magnetic effects (whereas the standard fewest switches surface hopping approach cannot capture such effects). We also find that a preconditioned QCLE can outperform the standard QCLE in certain cases, highlighting the fact that there is no unique QCLE. Finally, we also point out that one can construct a mean-field Ehrenfest algorithm using a phase-space representation similar to what is done for PSSH. These findings would appear extremely helpful as far as understanding and simulating nonadiabatic dynamics with complex-valued Hamiltonians and/or spin degeneracy.

Modeling Spin-Dependent Nonadiabatic Dynamics with Electronic Degeneracy: A Phase-Space Surface-Hopping Method

Nuclear Berry curvature effects emerge from electronic spin degeneracy and can lead to nontrivial spin-dependent (nonadiabatic) nuclear dynamics. However, such effects are not captured fully by any current mixed quantum-classical method such as fewest-switches surface hopping. In this work, we present a phase-space surface-hopping (PSSH) approach to simulate singlet-triplet intersystem crossing dynamics. We show that with a simple pseudodiabatic ansatz, a PSSH algorithm can capture the relevant Berry curvature effects and make predictions in agreement with exact quantum dynamics for a simple singlet-triplet model Hamiltonian. Thus, this approach represents an important step toward simulating photochemical and spin processes concomitantly, as relevant to intersystem crossing and spin-lattice relaxation dynamics.

A phase-space semiclassical approach for modeling nonadiabatic nuclear dynamics with electronic spin

Chemical relaxation phenomena, including photochemistry and electron transfer processes, form a vigorous area of research in which nonadiabatic dynamics plays a fundamental role. However, for electronic systems with spin degrees of freedom, there are few if any applicable and practical quasiclassical methods. Here, we show that for nonadiabatic dynamics with two electronic states and a complex-valued Hamiltonian that does not obey time-reversal symmetry (as relevant to many coupled nuclear-electronic-spin systems), the optimal semiclassical approach is to generalize Tully’s surface hopping dynamics from coordinate space to phase space. In order to generate the relevant phase-space adiabatic surfaces, one isolates a proper set of diabats, applies a phase gauge transformation, and then diagonalizes the total Hamiltonian (which is now parameterized by both R and P). The resulting algorithm is simple and valid in both the adiabatic and nonadiabatic limits, incorporating all Berry curvature effects. Most importantly, the resulting algorithm allows for the study of semiclassical nonadiabatic dynamics in the presence of spin–orbit coupling and/or external magnetic fields. One expects many simulations to follow as far as modeling cutting-edge experiments with entangled nuclear, electronic, and spin degrees of freedom, e.g., experiments displaying chiral-induced spin selectivity.

On the meaning of Berry force for unrestricted systems treated with mean-field electronic structure

We show that the Berry force as computed by an approximate, mean-field electronic structure can be meaningful if properly interpreted. In particular, for a model Hamiltonian representing a molecular system with an even number of electrons interacting via a two-body (Hubbard) interaction and a spin-orbit coupling, we show that a meaningful nonzero Berry force emerges whenever there is spin unrestriction-even though the Hamiltonian is real-valued and formally the on-diagonal single-surface Berry force must be zero. Moreover, if properly applied, this mean-field Berry force yields roughly the correct asymptotic motion for scattering through an avoided crossing. That being said, within the context of a ground-state calculation, several nuances do arise as far interpreting the Berry force correctly, and as a practical matter, the Berry force diverges near the Coulson-Fischer point (which can lead to numerical instabilities). We do not address magnetic fields here.