Exact factorization of the time-dependent electron-nuclear wave function

Challenges in Simulating Light-Induced Processes in DNA

In this contribution, we give a perspective on the main challenges in performing theoretical simulations of photoinduced phenomena within DNA and its molecular building blocks. We distinguish the different tasks that should be involved in the simulation of a complete DNA strand subject to UV irradiation: (i) stationary quantum chemical computations; (ii) the explicit description of the initial excitation of DNA with light; (iii) modeling the nonadiabatic excited state dynamics; (iv) simulation of the detected experimental observable; and (v) the subsequent analysis of the respective results. We succinctly describe the methods that are currently employed in each of these steps. While for each of them, there are different approaches with different degrees of accuracy, no feasible method exists to tackle all problems at once. Depending on the technique or combination of several ones, it can be problematic to describe the stacking of nucleobases, bond breaking and formation, quantum interferences and tunneling or even simply to characterize the involved wavefunctions. It is therefore argued that more method development and/or the combination of different techniques are urgently required. It is essential also to exercise these new developments in further studies on DNA and subsystems thereof, ideally comprising simulations of all of the different components that occur in the corresponding experiments.

Exact Factorization-Based Density Functional Theory of Electrons and Nuclei

The ground state energy of a system of electrons (r=r_{1},r_{2},…) and nuclei (R=R_{1},R_{2},…) is proven to be a variational functional of the electronic density n(r,R) and paramagnetic current density j_{p}(r,R) conditional on R, the nuclear wave function χ(R), an induced vector potential A_{μ}(R) and a quantum geometric tensor T_{μν}(R). n, j_{p}, A_{μ} and T_{μν} are defined in terms of the conditional electronic wave function Φ_{R}(r). The ground state (n,j_{p},χ,A_{μ},T_{μν}) can be calculated by solving self-consistently (i) conditional Kohn-Sham equations containing effective scalar and vector potentials v_{s}(r) and A_{xc}(r) that depend parametrically on R, (ii) the Schrödinger equation for χ(R), and (iii) Euler-Lagrange equations that determine T_{μν}. The theory is applied to the E⊗e Jahn-Teller model.

Coupled Electron-Nuclear Dynamics on H2+ within Time-Dependent Born-Oppenheimer Approximation

Quantum dynamical behavior of H₂⁺ in the presence of a linearly polarized, ultrashort, intense, infrared laser pulse has been studied by numerically solving the time-dependent Schrödinger equation with nuclear motion restricted in one-dimension along the direction of laser polarization and electronic motion in three-dimensions. On the basis of the time-dependent Born-Oppenheimer approximation, we have constructed time-dependent potentials for the ground electronic state (1sσ_g) of H₂⁺. Subsequent nuclear dynamics is then carried out on these field-dressed potential energy surfaces, and the dissociation dynamics is investigated. Our analyses reveal that although the electronic longitudinal degree of freedom plays the major role in governing the dissociation dynamics, contributions from the electronic transverse degree of freedom should also have to be taken into account to obtain accurate results. Also, modeling electron-nuclei Coulomb interactions in a one-dimensional calculation with an artificially chosen constant softening parameter leads to a discrepancy with the exact results. Comparing our results with other quantum and classical dynamical studies showed a good agreement with exact results.

Quantum-Classical Nonadiabatic Dynamics: Coupled- vs Independent-Trajectory Methods

Trajectory-based mixed quantum-classical approaches to coupled electron-nuclear dynamics suffer from well-studied problems such as the lack of (or incorrect account for) decoherence in the trajectory surface hopping method and the inability of reproducing the spatial splitting of a nuclear wave packet in Ehrenfest-like dynamics. In the context of electronic nonadiabatic processes, these problems can result in wrong predictions for quantum populations and in unphysical outcomes for the nuclear dynamics. In this paper, we propose a solution to these issues by approximating the coupled electronic and nuclear equations within the framework of the exact factorization of the electron-nuclear wave function. We present a simple quantum-classical scheme based on coupled classical trajectories and test it against the full quantum mechanical solution from wave packet dynamics for some model situations which represent particularly challenging problems for the above-mentioned traditional methods.

Confronting surface hopping molecular dynamics with Marcus theory for a molecular donor–acceptor system

We investigate the performance of fewest switches surface hopping (SH) in describing electron transfer (ET) for a molecular donor–acceptor system. Computer simulations are carried out for a wide range of reorganisation energy (λ), electronic coupling strength (H_ab) and driving force using our recently developed fragment orbital-based SH approach augmented with a simple decoherence correction. This methodology allows us to compute SH ET rates over more than four orders of magnitude, from the sub-picosecond to the nanosecond time regime. We find good agreement with semi-classical ET theory in the non-adiabatic ET regime. The correct scaling of the SH ET rate with electronic coupling strength is obtained and the Marcus inverted regime is reproduced, in line with previously reported results for a spin-boson model. Yet, we find that the SH ET rate falls below the semi-classical ET rate in the adiabatic regime, where the free energy barrier is in the order of k_BT in our simulations. We explain this by first signatures of non-exponential population decay of the initial charge state. For even larger electronic couplings (H_ab = λ/2), the free energy barrier vanishes and ET rates are no longer defined. At this point we observe a crossover from ET on the vibronic time scale to charge relaxation on the femtosecond time scale that is well described by thermally averaged Rabi oscillations. The extension of the analysis from the non-adiabatic limit to large electronic couplings and small or even vanishing activation barriers is relevant for our understanding of charge transport in organic semiconductors.

Exact Potential Driving the Electron Dynamics in Enhanced Ionization of H(2)(+)

It was recently shown that the exact factorization of the electron-nuclear wave function allows the construction of a Schrödinger equation for the electronic system, in which the potential contains exactly the effect of coupling to the nuclear degrees of freedom and any external fields. Here we study the exact potential acting on the electron in charge-resonance enhanced ionization in a model one-dimensional H(2)(+) molecule. We show there can be significant differences between the exact potential and that used in the traditional quasistatic analyses, arising from nonadiabatic coupling to the nuclear system, and that these are crucial to include for accurate simulations of time-resolved ionization dynamics and predictions of the ionization yield.

Laser-induced electron localization in H₂⁺: mixed quantum-classical dynamics based on the exact time-dependent potential energy surface

We study the exact nuclear time-dependent potential energy surface (TDPES) for laser-induced electron localization with a view to eventually developing a mixed quantum-classical dynamics method for strong-field processes. The TDPES is defined within the framework of the exact factorization [A. Abedi, N. T. Maitra, and E. K. U. Gross, Phys. Rev. Lett., 2010, 105, 123002] and contains the exact effect of the couplings to the electronic subsystem and to any external fields within a scalar potential. We compare its features with those of the quasistatic potential energy surfaces (QSPES) often used to analyse strong-field processes. We show that the gauge-independent component of the TDPES has a mean-field-like character very close to the density-weighted average of the QSPESs. Oscillations in this component are smoothened out by the gauge-dependent component, and both components are needed to yield the correct force on the nuclei. Once the localization begins to set in, the gradient of the exact TDPES tracks one QSPES and then switches to the other, similar to the description provided by surface-hopping between QSPESs. We show that evolving an ensemble of classical nuclear trajectories on the exact TDPES accurately reproduces the exact dynamics. This study suggests that the mixed quantum-classical dynamics scheme based on evolving multiple classical nuclear trajectories on the exact TDPES will be a novel and useful method to simulate strong field processes.

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

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.

The exact forces on classical nuclei in non-adiabatic charge transfer

The decomposition of electronic and nuclear motion presented in Abedi et al. [Phys. Rev. Lett. 105, 123002 (2010)] yields a time-dependent potential that drives the nuclear motion and fully accounts for the coupling to the electronic subsystem. Here, we show that propagation of an ensemble of independent classical nuclear trajectories on this exact potential yields dynamics that are essentially indistinguishable from the exact quantum dynamics for a model non-adiabatic charge transfer problem. We point out the importance of step and bump features in the exact potential that are critical in obtaining the correct splitting of the quasiclassical nuclear wave packet in space after it passes through an avoided crossing between two Born-Oppenheimer surfaces and analyze their structure. Finally, an analysis of the exact potentials in the context of trajectory surface hopping is presented, including preliminary investigations of velocity-adjustment and the force-induced decoherence effect.

Electronic Friction-Based Vibrational Lifetimes of Molecular Adsorbates: Beyond the Independent-Atom Approximation

We assess the accuracy of vibrational damping rates of diatomic adsorbates on metal surfaces as calculated within the local-density friction approximation (LDFA). An atoms-in-molecules (AIM) type charge partitioning scheme accounts for intramolecular contributions and overcomes the systematic underestimation of the nonadiabatic losses obtained within the prevalent independent-atom approximation. The quantitative agreement obtained with theoretical and experimental benchmark data suggests the LDFA-AIM scheme as an efficient and reliable approach to account for electronic dissipation in ab initio molecular dynamics simulations of surface chemical reactions.

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.

Semiclassical quantization of nonadiabatic systems with hopping periodic orbits

We present a semiclassical quantization condition, i.e., quantum-classical correspondence, for steady states of nonadiabatic systems consisting of fast and slow degrees of freedom (DOFs) by extending Gutzwiller's trace formula to a nonadiabatic form. The quantum-classical correspondence indicates that a set of primitive hopping periodic orbits, which are invariant under time evolution in the phase space of the slow DOF, should be quantized. The semiclassical quantization is then applied to a simple nonadiabatic model and accurately reproduces exact quantum energy levels. In addition to the semiclassical quantization condition, we also discuss chaotic dynamics involved in the classical limit of nonadiabatic dynamics.

Non-adiabatic effects in thermochemistry, spectroscopy and kinetics: the general importance of all three Born–Oppenheimer breakdown corrections

Analytical and numerical solutions describing Born–Oppenheimer breakdown in a simple, widely applicable, model depict shortcomings in modern computational methods.

Is the molecular Berry phase an artifact of the Born-Oppenheimer approximation?

We demonstrate that the molecular Berry phase and the corresponding nonanalyticity in the electronic Born-Oppenheimer wave function is, in general, not a true topological feature of the exact solution of the full electron-nuclear Schrödinger equation. For a numerically exactly solvable model we show that a nonanalyticity, and the associated geometric phase, only appear in the limit of infinite nuclear mass, while a perfectly smooth behavior is found for any finite nuclear mass.

Correlated electron-nuclear dynamics with conditional wave functions

The molecular Schrödinger equation is rewritten in terms of nonunitary equations of motion for the nuclei (or electrons) that depend parametrically on the configuration of an ensemble of generally defined electronic (or nuclear) trajectories. This scheme is exact and does not rely on the tracing out of degrees of freedom. Hence, the use of trajectory-based statistical techniques can be exploited to circumvent the calculation of the computationally demanding Born-Oppenheimer potential-energy surfaces and nonadiabatic coupling elements. The concept of the potential-energy surface is restored by establishing a formal connection with the exact factorization of the full wave function. This connection is used to gain insight from a simplified form of the exact propagation scheme.

Dynamical steps that bridge piecewise adiabatic shapes in the exact time-dependent potential energy surface

We study the exact time-dependent potential energy surface (TDPES) in the presence of strong nonadiabatic coupling between the electronic and nuclear motion. The concept of the TDPES emerges from the exact factorization of the full electron-nuclear wave function [A. Abedi, N. T. Maitra, and E. K. U. Gross, Phys. Rev. Lett. 105, 123002 (2010)]. Employing a one-dimensional model system, we show that the TDPES exhibits a dynamical step that bridges between piecewise adiabatic shapes. We analytically investigate the position of the steps and the nature of the switching between the adiabatic pieces of the TDPES.

Theoretical methods for ultrafast spectroscopy

Time-resolved spectroscopy in the femtosecond and attosecond time domain is a tool to unravel the dynamics of nuclear and electronic motion in molecular systems. Theoretical insight into the underlying physical processes is ideally gained by solving the time-dependent Schrödinger equation. In this work, methods currently used to solve this equation are reviewed in a compact presentation. These methods involve numerical representations of wavefunctions and operators, the calculation of time evolution operators, the setting up of the Hamiltonian operators and the types of coordinates to be used hereto. The advantages and disadvantages of some methods are discussed.

Correlated electron-nuclear dynamics: exact factorization of the molecular wavefunction

It was recently shown [A. Abedi, N. T. Maitra, and E. K. U. Gross, Phys. Rev. Lett. 105, 123002 (2010)] that the complete wavefunction for a system of electrons and nuclei evolving in a time-dependent external potential can be exactly factorized into an electronic wavefunction and a nuclear wavefunction. The concepts of an exact time-dependent potential energy surface (TDPES) and exact time-dependent vector potential emerge naturally from the formalism. Here, we present a detailed description of the formalism, including a full derivation of the equations that the electronic and nuclear wavefunctions satisfy. We demonstrate the relationship of this exact factorization to the traditional Born-Oppenheimer expansion. A one-dimensional model of the H(2)(+) molecule in a laser field shows the usefulness of the exact TDPES in interpreting coupled electron-nuclear dynamics: we show how features of its structure indicate the mechanism of dissociation. We compare the exact TDPES with potential energy surfaces from the time-dependent Hartree-approach, and also compare traditional Ehrenfest dynamics with Ehrenfest dynamics on the exact TDPES.

Regulating energy transfer of excited carriers and the case for excitation-induced hydrogen dissociation on hydrogenated graphene

Understanding and controlling of excited carrier dynamics is of fundamental and practical importance, particularly in photochemistry and solar energy applications. However, theory of energy relaxation of excited carriers is still in its early stage. Here, using ab initio molecular dynamics (MD) coupled with time-dependent density functional theory, we show a coverage-dependent energy transfer of photoexcited carriers in hydrogenated graphene, giving rise to distinctively different ion dynamics. Graphene with sparsely populated H is difficult to dissociate due to inefficient transfer of the excitation energy into kinetic energy of the H. In contrast, H can easily desorb from fully hydrogenated graphane. The key is to bring down the H antibonding state to the conduction band minimum as the band gap increases. These results can be contrasted to those of standard ground-state MD that predict H in the sparse case should be much less stable than that in fully hydrogenated graphane. Our findings thus signify the importance of carrying out explicit electronic dynamics in excited-state simulations.

Theory of high harmonic generation for probing time-resolved large-amplitude molecular vibrations with ultrashort intense lasers

We present a theory that incorporates the vibrational degrees of freedom in a high-order harmonic generation (HHG) process with ultrashort intense laser pulses. In this model, laser-induced time-dependent transition dipoles for each fixed molecular geometry are added coherently, weighted by the laser-driven time-dependent nuclear wave packet distribution. We show that the nuclear distribution can be strongly modified by the HHG driving laser. The validity of this model is first checked against results from the numerical solution of the time-dependent Schrödinger equation for a simple model system. We show that in combination with the established quantitative rescattering theory this model is able to reproduce the time-resolved pump-probe HHG spectra of N(2)O(4) reported by Li et al. [Science 322, 1207 (2008)].