arXiv:1504.00250 · Physical Review Letters

Coupled-Trajectory Quantum-Classical Approach to Electronic Decoherence in Nonadiabatic Processes

Seung Kyu Min, Federica Agostini, Eberhard K. U. Gross

Published Aug 10, 2015• 180 citations

Abstract

We present a novel quantum-classical approach to nonadiabatic dynamics, deduced from the coupled electronic and nuclear equations in the framework of the exact factorization of the electron-nuclear wave function. The method is based on the quasiclassical interpretation of the nuclear wave function, whose phase is related to the classical momentum and whose density is represented in terms of classical trajectories. In this approximation, electronic decoherence is naturally induced as an effect of the coupling to the nuclei and correctly reproduces the expected quantum behavior. Moreover, the splitting of the nuclear wave packet is captured as a consequence of the correct approximation of the time-dependent potential of the theory. This new approach offers a clear improvement over Ehrenfest-like dynamics. The theoretical derivation presented in this Letter is supported by numerical results that are compared to quantum mechanical calculations.

DOI: 10.1103/physrevlett.115.073001

Related papers

Coupled-Trajectory Mixed Quantum-Classical Algorithm: A Deconstruction

Gossel Graeme H., Federica Agostini, Neepa T. Maitra

We analyze a mixed quantum-classical algorithm recently derived from the exact factorization equations [Min, Agostini, Gross, PRL 115, 073001 (2015)] to show the role of the different terms in the algorithm in bringing about decoherence and wavepacket branching. The algorithm has the structure of Ehrenfest equations plus a "coupled-trajectory" term for both the electronic and nuclear equations, and we analyze the relative roles played by the different nonadiabatic terms in these equations, including how they are computed in practice. In particular, we show that while the coupled-trajectory term in the electronic equation is essential in yielding accurate dynamics, that in the nuclear equation has a much smaller effect. A decoherence time is extracted from the electronic equations and compared with that of augmented fewest-switches surface-hopping. We revisit a series of nonadiabatic Tully model systems to illustrate our analysis.

Jan 1, 201845 citations

Coupled-Trajectory Mixed Quantum-Classical Algorithm: A Deconstruction

G. H. Gossel, Federica Agostini, Neepa T. Maitra

arXiv:1805.03534Journal of Chemical Theory and ComputationJul 31, 20181 citation

Electronic decoherence along a single nuclear trajectory

Matisse Wei-Yuan Tu, Eberhard K. U. Gross

We describe an approach to subsystem decoherence without the usual tracing out of the environment. The subsystem of focus is described entirely by a pure state evolving nonunitarily along a single classical trajectory of its environment. The approach is deduced from the exact factorization framework for arbitrary systems of electrons and nuclei. The nonunitarity of the electronic dynamics arises exclusively from nonadiabatic correlations between electrons and nuclei. We demonstrate that the approach correctly describes the coherence gain and the subsequent decoherence for the example of a nuclear trajectory passing through an avoided crossing, the prototypical case where single-trajectory Ehrenfest dynamics fails to produce decoherence.

arXiv:2502.08247Physical Review ResearchSep 17, 20252 citations

Quantum Trajectories for the Dynamics in the Exact Factorization Framework: A Proof-of-Principle Test

Francesco Talotta, Federica Agostini, Giovanni Ciccotti

In the framework of the exact factorization of the time-dependent electron-nuclear wave function, we investigate the possibility of solving the nuclear time-dependent Schrödinger equation based on trajectories. The nuclear equation is separated in a Hamilton-Jacobi equation for the phase of the wave function, and a continuity equation for its (squared) modulus. For illustrative adiabatic and nonadiabatic one-dimensional models, we implement a procedure to follow the evolution of the nuclear density along the characteristics of the Hamilton-Jacobi equation. Those characteristics are referred to as quantum trajectories, since they are generated via ordinary differential equations similar to Hamilton's equations, but including the so-called quantum potential, and they can be used to reconstruct exactly the quantum-mechanical nuclear wave function, provided infinite initial conditions are propagated in time.

The Journal of Physical Chemistry AJul 30, 202024 citations

Thawed Gaussian Ehrenfest dynamics

Jiří J. L. Vaníček

Ehrenfest dynamics is a widely used mixed quantum--classical approach for nonadiabatic molecular dynamics, whereas thawed Gaussian wavepacket dynamics provides an efficient semiclassical description of adiabatic nuclear quantum dynamics. Here we describe thawed Gaussian Ehrenfest dynamics (TGED), which unifies and generalizes these two methods to capture both electronic nonadiabaticity and nuclear quantum effects within a single framework. The fully variational formulation of TGED is derived by applying the time-dependent variational principle to a Hartree product of electronic and Gaussian nuclear wavepackets. Replacing the effective locally quadratic molecular potential obtained from this variational treatment by alternative effective locally quadratic potentials yields an infinite family of TGED methods, of which we present several members. We analyze the limiting cases of the general formalism and show, in particular, that it reduces to conventional Ehrenfest dynamics in the classical limit for the nuclei and to thawed Gaussian wavepacket dynamics in the absence of electronic coupling. Finally, we present explicit geometric integrators for the entire family of methods and identify the conditions under which the different approximations become exact.

physics.chem-phCross-listedarXiv:2607.10847Jul 14, 2026

Ehrenfest and classical path dynamics with decoherence and detailed balance

Parmeet Nijjar, Joanna Jankowska, Oleg V. Prezhdo

We present a semiclassical approach for nonadiabatic molecular dynamics based on the Ehrenfest method with corrections for decoherence and detailed balance. Decoherence is described via a coherence penalty functional that drives dynamics away from regions in Hilbert space characterized by large values of coherences. Detailed balance is incorporated by modification of the off-diagonal matrix elements with a quantum correction factor used in semiclassical approximations to quantum time-correlation functions. Both decoherence and detailed balance corrections introduce nonlinear terms to the Schrödinger equation. At the same time, the simplicity of fully deterministic dynamics and a single trajectory for each initial condition is preserved. In contrast, surface hopping is stochastic and requires averaging over multiple realization of the stochastic process for each initial condition. The Ehrenfest-decoherence-detailed-balance (Ehrenfest-DDB) method is adapted to the classical path approximation and ab initio time-dependent density functional theory and applied to an experimentally studied nanoscale system consisting of a fluorophore molecule and an scanning tunneling microscopy tip and undergoing current-induced charge injection, cooling, and recombination. Ehrenfest-DDB produces time scales that are similar to those obtained with decoherence induced surface hopping, which is a popular nonadiabatic molecular dynamics technique applied to condensed matter. At long times, Ehrenfest-DDB dynamics slows down considerably because the detailed balance correction makes off-diagonal elements go to zero on approach to Boltzmann equilibrium. The Ehrenfest-DDB technique provides efficient means to study quantum dynamics in large systems.

The Journal of Chemical PhysicsMay 28, 201957 citations