Quantum-classical description of environmental effects on electronic dynamics at conical intersections

A first principles derivation of energy-conserving momentum jumps in surface hopping simulations

The fewest switches surface hopping (FSSH) method proposed by Tully in 1990 [Tully, J. Chem. Phys. 93, 1061 (1990)]-along with its many later variations-forms the basis for most practical simulations of molecular dynamics with electronic transitions in realistic systems. Despite its popularity, a rigorous formal derivation of the algorithm has yet to be achieved. In this paper, we derive the energy-conserving momentum jumps employed by FSSH from the perspective of quantum trajectory surface hopping (QTSH) [Martens, J. Phys. Chem. A 123, 1110 (2019)]. In the limit of localized nonadiabatic transitions, simple mathematical and physical arguments allow the FSSH algorithm to be derived from first principles. For general processes, the quantum forces characterizing the QTSH method provide accurate results for nonadiabatic dynamics with rigorous energy conservation, at the ensemble level, within the consistency of the underlying stochastic surface hopping without resorting to the artificial momentum rescaling of FSSH.

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.

Geometric Phase Effects in Nonadiabatic Dynamics near Conical Intersections

Dynamical consideration that goes beyond the common Born-Oppenheimer approximation (BOA) becomes necessary when energy differences between electronic potential energy surfaces become small or vanish. One of the typical scenarios of the BOA breakdown in molecules beyond diatomics is a conical intersection (CI) of electronic potential energy surfaces. CIs provide an efficient mechanism for radiationless electronic transitions: acting as "funnels" for the nuclear wave function, they enable rapid conversion of the excessive electronic energy into the nuclear motion. In addition, CIs introduce nontrivial geometric phases (GPs) for both electronic and nuclear wave functions. These phases manifest themselves in change of the wave function signs if one considers an evolution of the system around the CI. This sign change is independent of the shape of the encircling contour and thus has a topological character. How these extra phases affect nonadiabatic dynamics is the main question that is addressed in this Account. We start by considering the simplest model providing the CI topology: two-dimensional two-state linear vibronic coupling model. Selecting this model instead of a real molecule has the advantage that various dynamical regimes can be easily modeled in the model by varying parameters, whereas any fixed molecule provides the system specific behavior that may not be very illustrative. After demonstrating when GP effects are important and how they modify the dynamics for two sets of initial conditions (starting from the ground and excited electronic states), we give examples of molecular systems where the described GP effects are crucial for adequate description of nonadiabatic dynamics. Interestingly, although the GP has a topological character, the extent to which accounting for GPs affect nuclear dynamics profoundly depends on topography of potential energy surfaces. Understanding an extent of changes introduced by the GP in chemical dynamics poses a problem of capturing GP effects by approximate methods of simulating nonadiabatic dynamics that can go beyond simple models. We assess the performance of both fully quantum (wave packet dynamics) and quantum-classical (surface-hopping, Ehrenfest, and quantum-classical Liouville equation) approaches in various cases where GP effects are important. It has been identified that the key to success in approximate methods is a method organization that prevents the quantum nuclear kinetic energy operator to act directly on adiabatic electronic wave functions.

Accurate Long-Time Mixed Quantum-Classical Liouville Dynamics via the Transfer Tensor Method

In this Letter, we combine the recently introduced transfer tensor method with the mixed quantum-classical Liouville method. The resulting protocol provides an accurate, general, flexible and robust new route for simulating the reduced dynamics of the quantum subsystem for arbitrarily long times, starting with computationally feasible short-time mixed quantum-classical Liouville dynamical maps. The accuracy and feasibility of the methodology are demonstrated on a spin-boson benchmark model.

Simulation of femtosecond two-dimensional electronic spectra of conical intersections

We have simulated femtosecond two-dimensional (2D) electronic spectra for an excited-state conical intersection using the wave-function version of the equation-of-motion phase-matching approach. We show that 2D spectra at fixed values of the waiting time provide information on the structure of the vibronic eigenstates of the conical intersection, while the evolution of the spectra with the waiting time reveals predominantly ground-state wave-packet dynamics. The results show that 2D spectra of conical intersection systems differ significantly from those obtained for chromophores with well separated excited-state potential-energy surfaces. The spectral signatures which can be attributed to conical intersections are discussed.

Quantum dynamics in open quantum-classical systems

Often quantum systems are not isolated and interactions with their environments must be taken into account. In such open quantum systems these environmental interactions can lead to decoherence and dissipation, which have a marked influence on the properties of the quantum system. In many instances the environment is well-approximated by classical mechanics, so that one is led to consider the dynamics of open quantum-classical systems. Since a full quantum dynamical description of large many-body systems is not currently feasible, mixed quantum-classical methods can provide accurate and computationally tractable ways to follow the dynamics of both the system and its environment. This review focuses on quantum-classical Liouville dynamics, one of several quantum-classical descriptions, and discusses the problems that arise when one attempts to combine quantum and classical mechanics, coherence and decoherence in quantum-classical systems, nonadiabatic dynamics, surface-hopping and mean-field theories and their relation to quantum-classical Liouville dynamics, as well as methods for simulating the dynamics.

Analysis of geometric phase effects in the quantum-classical Liouville formalism

We analyze two approaches to the quantum-classical Liouville (QCL) formalism that differ in the order of two operations: Wigner transformation and projection onto adiabatic electronic states. The analysis is carried out on a two-dimensional linear vibronic model where geometric phase (GP) effects arising from a conical intersection profoundly affect nuclear dynamics. We find that the Wigner-then-Adiabatic (WA) QCL approach captures GP effects, whereas the Adiabatic-then-Wigner (AW) QCL approach does not. Moreover, the Wigner transform in AW-QCL leads to an ill-defined Fourier transform of double-valued functions. The double-valued character of these functions stems from the nontrivial GP of adiabatic electronic states in the presence of a conical intersection. In contrast, WA-QCL avoids this issue by starting with the Wigner transform of single-valued quantities of the full problem. As a consequence, GP effects in WA-QCL can be associated with a dynamical term in the corresponding equation of motion. Since the WA-QCL approach uses solely the adiabatic potentials and non-adiabatic derivative couplings as an input, our results indicate that WA-QCL can capture GP effects in two-state crossing problems using first-principles electronic structure calculations without prior diabatization or introduction of explicit phase factors.

Can we derive Tully's surface-hopping algorithm from the semiclassical quantum Liouville equation? Almost, but only with decoherence

In this article, we demonstrate that Tully's fewest-switches surface hopping (FSSH) algorithm approximately obeys the mixed quantum-classical Liouville equation (QCLE), provided that several conditions are satisfied--some major conditions, and some minor. The major conditions are: (1) nuclei must be moving quickly with large momenta; (2) there cannot be explicit recoherences or interference effects between nuclear wave packets; (3) force-based decoherence must be added to the FSSH algorithm, and the trajectories can no longer rigorously be independent (though approximations for independent trajectories are possible). We furthermore expect that FSSH (with decoherence) will be most robust when nonadiabatic transitions in an adiabatic basis are dictated primarily by derivative couplings that are presumably localized to crossing regions, rather than by small but pervasive off-diagonal force matrix elements. In the end, our results emphasize the strengths of and possibilities for the FSSH algorithm when decoherence is included, while also demonstrating the limitations of the FSSH algorithm and its inherent inability to follow the QCLE exactly.

Adiabatic and nonadiabatic contributions to the energy of a system subject to a time-dependent perturbation: complete separation and physical interpretation

When a time-dependent perturbation acts on a quantum system that is initially in the nondegenerate ground state ∣0> of an unperturbed Hamiltonian H(0), the wave function acquires excited-state components ∣k> with coefficients c(k)(t) exp(-iE(k)t/ℏ), where E(k) denotes the energy of the unperturbed state ∣k>. It is well known that each coefficient c(k)(t) separates into an adiabatic term a(k)(t) that reflects the adjustment of the ground state to the perturbation--without actual transitions--and a nonadiabatic term b(k)(t) that yields the probability amplitude for a transition to the excited state. In this work, we prove that the energy at any time t also separates completely into adiabatic and nonadiabatic components, after accounting for the secular and normalization terms that appear in the solution of the time-dependent Schrödinger equation via Dirac's method of variation of constants. This result is derived explicitly through third order in the perturbation. We prove that the cross-terms between the adiabatic and nonadiabatic parts of c(k)(t) vanish, when the energy at time t is determined as an expectation value. The adiabatic term in the energy is identical to the total energy obtained from static perturbation theory, for a system exposed to the instantaneous perturbation λH'(t). The nonadiabatic term is a sum over excited states ∣k> of the transition probability multiplied by the transition energy. By evaluating the probabilities of transition to the excited eigenstates ∣k'(t)> of the instantaneous Hamiltonian H(t), we provide a physically transparent explanation of the result for E(t). To lowest order in the perturbation parameter λ, the probability of finding the system in state ∣k'(t)> is given by λ(2) ∣b(k)(t)∣(2). At third order, the transition probability depends on a second-order transition coefficient, derived in this work. We indicate expected differences between the results for transition probabilities obtained from this work and from Fermi's golden rule.