arXiv:1711.08197 · Physical review. B./Physical review. B

Coupled forward-backward trajectory approach for nonequilibrium electron-ion dynamics

Shunsuke A. Sato, Aaron Kelly, Ángel Rubio

Published Apr 25, 2018• 14 citations

Abstract

We introduce a simple ansatz for the wave function of a many-body system based on coupled forward and backward propagating semiclassical trajectories. This method is primarily aimed at, but not limited to, treating nonequilibrium dynamics in electron-phonon systems. The time evolution of the system is obtained from the Euler-Lagrange variational principle, and we show that this ansatz yields Ehrenfest mean-field theory in the limit that the forward and backward trajectories are orthogonal, and in the limit that they coalesce. We investigate accuracy and performance of this method by simulating electronic relaxation in the spin-boson model and the Holstein model. Although this method involves only pairs of semiclassical trajectories, it shows a substantial improvement over mean-field theory, capturing quantum coherence of nuclear dynamics as well as electron-nuclear correlations. This improvement is particularly evident in nonadiabatic systems, where the accuracy of this coupled trajectory method extends well beyond the perturbative electron-phonon coupling regime. This approach thus provides an attractive route forward to the ab initio description of relaxation processes, such as thermalization, in condensed phase systems.

DOI: 10.1103/physrevb.97.134308

Related papers

Stochastic semiclassical theory for nonequilibrium electron-phonon coupled systems

Antonio Picano, Francesco Grandi, Philipp Werner +1 more

We discuss a semiclassical approach to solve the quantum impurity model within nonequilibrium dynamical mean-field theory for electron-lattice models. The effect of electronic fluctuations on the phonon is kept beyond Ehrenfest dynamics, leading to a stochastic phonon evolution with damping and noise terms that are self-consistently determined by the electronic correlation functions in the fluctuating phonon field. Together with a solution of the electronic model based on a nonperturbative quantum Boltzmann equation, the approach can be used to address the coupled dynamics of the electrons and the lattice during photoinduced phase transitions. Results for the Anderson-Holstein model are benchmarked against numerically exact quantum Monte Carlo data. We find good agreement for the phonon distribution function at temperatures comparable to the charge ordering temperature. The general formulation can be extended to models with electron-electron interactions or multiorbital systems.

arXiv:2209.00428Physical review. B./Physical review. BJul 10, 202311 citations

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

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

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.

arXiv:1504.00250Physical Review LettersAug 10, 2015180 citations

Ab-initio variational wave functions for the time-dependent many-electron Schrödinger equation

Jannes Nys, Gabriel Pescia, Alessandro Sinibaldi +1 more

Understanding the real-time evolution of many-electron quantum systems is essential for studying dynamical properties in condensed matter, quantum chemistry, and complex materials, yet it poses a significant theoretical and computational challenge. Our work introduces a variational approach for fermionic time-dependent wave functions, surpassing mean-field approximations by accurately capturing many-body correlations. We employ time-dependent Jastrow factors and backflow transformations, enhanced through neural networks parameterizations. To compute the optimal time-dependent parameters, we employ the time-dependent variational Monte Carlo technique and introduce a new method based on Taylor-root expansions of the propagator, enhancing the accuracy of our simulations. The approach is demonstrated in three distinct systems. In all cases, we show clear signatures of many-body correlations in the dynamics. The results showcase the ability of our variational approach to accurately describe the time evolution, providing insight into quantum dynamical effects in interacting electronic systems, beyond the capabilities of mean-field.

arXiv:2403.07447Nature CommunicationsOct 30, 202423 citations

Energy Flow during Relaxation in an Electron–Phonon System with Multiple Modes: A Nonequilibrium Green’s Function Study

Ken Inayoshi, Akihisa Koga, Yuta Murakami

We investigate an energy flow in an extended Holstein model describing electron systems coupled to hot-phonons and heat-bath phonons. To analyze the relaxation process after the photo-excitation of electrons, we employ the nonequilibrium dynamical mean-field theory (DMFT). We find the backward energy flow during the relaxation, where the direction of energy transfer between electrons and hot-phonons is reversed. To clarify the microscopic mechanism of the backward energy flow, we introduce the approximated energy flows, which are calculated with the gradient and quasiparticle approximations and are related to the nonequilibrium distribution functions. We compare these approximated energy flows with the full energy flows calculated from the nonequilibrium DMFT. We find that, in the weak electron-hot-phonon coupling regime, the full and approximated energy flows are almost the same, meaning that the relaxation dynamics can be correctly understood in terms of the nonequilibrium distribution functions. As the strength of the electron-hot-phonon coupling increases, the approximated energy flow fails to qualitatively reproduce the full energy flow. This indicates that the microscopic origin of the energy flow cannot be solely explained by the nonequilibrium distribution functions. By comparing the energy flows with different levels of approximation, we reveal the role of the gradient and quasiparticle approximations.

arXiv:2307.12543Journal of the Physical Society of JapanNov 1, 20231 citation

In and Out-of-Equilibrium Ab Initio Theory of Electrons and Phonons

Gianluca Stefanucci, Robert van Leeuwen, Enrico Perfetto

In this work, we lay down the ab initio many-body quantum theory of electrons and phonons in equilibrium as well as in steady-state or time-varying settings.Our focus is on the harmonic approximation, but the developed tools can readily incorporate anharmonic effects.We begin by showing the necessity of determining the ab initio Hamiltonian in a self-consistent manner to ensure the existence of an equilibrium state.We then identify the correct partitioning into a "noninteracting" and an "interacting" part to carry out diagrammatic expansions in terms of dressed propagators and screened interactions.The final outcome is the finite-temperature nonequilibrium extension of the Hedin equations, showcasing the emergence of the coupling between electrons and coherent phonons through the time-local Ehrenfest diagram.The Hedin equations have limited practical utility for real-time simulations of systems driven out of equilibrium by external fields.To overcome this limitation, we leverage the versatility of the diagrammatic approach to generate a closed system of differential equations for the dressed propagators and nuclear displacements.These are the Kadanoff-Baym equations for electrons and phonons.The formalism naturally merges with the theory of conserving approximations, which guarantee the satisfaction of the continuity equation and the conservation of total energy during time evolution.As an example, we show that the popular Born-Oppenheimer approximation is not conserving whereas its dynamical extension is conserving, provided that the electrons are treated in the Fan-Migdal approximation with a dynamically screened electronphonon coupling.We also derive the formal solution of the Kadanoff-Baym equations for nonequilibrium steady states, which is useful for studies in photovoltaics and optoelectronics.Interestingly, the expansion of the phononic Green's function around the quasiphonon energies points to a correlation-induced splitting of the phonon dispersion in materials with no time-reversal invariance.

Physical Review XSep 11, 202332 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