Federica Agostini

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.

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.

Mixed quantum-classical dynamics from the exact decomposition of electron-nuclear motion

We present a novel mixed quantum-classical approach to the coupled electron-nuclear dynamics based on the exact factorisation of the electron-nuclear wave function, recently proposed in Abedi A., Maitra N. T. and Gross E. K. U., Phys. Rev. Lett., 105 (2010) 123002. In this framework, the correct classical limit of the nuclear dynamics is worked out by taking the classical limit of the exact time-dependent Schr?dinger equation satisfied by the nuclear wave function. The effect of the time-dependent scalar and vector potentials, representing the exact electronic back-reaction on the nuclear subsystem, is consistently derived within the classical approximation. We examine with an example the performance of the proposed mixed quantum-classical scheme in comparison with exact calculations.

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.

Classical nuclear motion coupled to electronic non-adiabatic transitions

Based on the exact factorization of the electron-nuclear wave function, we have recently proposed a mixed quantum-classical scheme [A. Abedi, F. Agostini, and E. K. U. Gross, Europhys. Lett. 106, 33001 (2014)] to deal with non-adiabatic processes. Here we present a comprehensive description of the formalism, including the full derivation of the equations of motion. Numerical results are presented for a model system for non-adiabatic charge transfer in order to test the performance of the method and to validate the underlying approximations.

Mixed quantum-classical dynamics on the exact time-dependent potential energy surface: a fresh look at non-adiabatic processes

The exact nuclear time-dependent potential energy surface arises from the exact decomposition of the electronic and nuclear motions, recently presented in [A. Abedi, N.T. Maitra, and E.K.U. Gross, Phys. Rev. Lett. 105, 123002 (2010)]. Such time-dependent potential drives nuclear motion and fully accounts for the coupling to the electronic subsystem. We investigate the features of the potential in the context of electronic non-adiabatic processes and employ it to study the performance of the classical approximation on nuclear dynamics. We observe that the potential, after the nuclear wave packet splits at an avoided crossing, develops dynamical steps connecting different regions, along the nuclear coordinate, in which it has the same slope as one or the other adiabatic surface. A detailed analysis of these steps is presented for systems with different non-adiabatic coupling strength. The exact factorisation of the electron-nuclear wave function is at the basis of the decomposition. In particular, the nuclear part is the true nuclear wave function, solution of a time-dependent Schrödinger equation and leading to the exact many-body density and current density. As a consequence, the Ehrenfest theorem can be extended to the nuclear subsystem and Hamiltonian, as discussed here with an analytical derivation and numerical results.

Nuclear velocity perturbation theory for vibrational circular dichroism: An approach based on the exact factorization of the electron-nuclear wave function.

The nuclear velocity perturbation theory (NVPT) for vibrational circular dichroism (VCD) is derived from the exact factorization of the electron-nuclear wave function. This new formalism offers an exact starting point to include correction terms to the Born-Oppenheimer (BO) form of the molecular wave function, similar to the complete-adiabatic approximation. The corrections depend on a small parameter that, in a classical treatment of the nuclei, is identified as the nuclear velocity. Apart from proposing a rigorous basis for the NVPT, we show that the rotational strengths, related to the intensity of the VCD signal, contain a new contribution beyond-BO that can be evaluated with the NVPT and that only arises when the exact factorization approach is employed. Numerical results are presented for chiral and non-chiral systems to test the validity of the approach.

Ab Initio Nonadiabatic Dynamics With Coupled Trajectories: A Rigorous Approach to Quantum (De)Coherence

We report the first nonadiabatic molecular dynamics study based on the exact factorization of the electron-nuclear wave function. Our approach (a coupled-trajectory mixed quantum-classical, CT-MQC, scheme) is based on the quantum-classical limit derived from systematic and controlled approximations to the full quantum-mechanical problem formulated in the exact-factorization framework. Its strength is the ability to correctly capture quantum (de)coherence effects in a trajectory-based approach to excited-state dynamics. We show this by benchmarking CT-MQC dynamics against a revised version of the popular fewest-switches surface-hopping scheme that is able to fix its well-documented overcoherence issue. The CT-MQC approach is successfully applied to investigation of the photochemistry (ring-opening) of oxirane in the gas phase, analyzing in detail the role of decoherence. This work represents a significant step forward in the establishment of the exact factorization as a powerful tool to study excited-state dynamics, not only for interpretation purposes but mainly for nonadiabatic ab initio molecular dynamics simulations.

Electronic Flux Density beyond the Born–Oppenheimer Approximation

In the Born-Oppenheimer approximation, the electronic wave function is typically real-valued and hence the electronic flux density (current density) seems to vanish. This is unfortunate for chemistry, because it precludes the possibility to monitor the electronic motion associated with the nuclear motion during chemical rearrangements from a Born-Oppenheimer simulation of the process. We study an electronic flux density obtained from a correction to the electronic wave function. This correction is derived via nuclear velocity perturbation theory applied in the framework of the exact factorization of electrons and nuclei. To compute the correction, only the ground state potential energy surface and the electronic wave function are needed. For a model system, we demonstrate that this electronic flux density approximates the true one very well, for coherent tunneling dynamics as well as for over-the-barrier scattering, and already for mass ratios between electrons and nuclei that are much larger than the true mass ratios.

The adiabatic limit of the exact factorization of the electron-nuclear wave function

We propose a procedure to analyze the relation between the exact factorization of the electron-nuclear wave function and the Born-Oppenheimer approximation. We define the adiabatic limit as the limit of infinite nuclear mass. To this end, we introduce a unit system that singles out the dependence on the electron-nuclear mass ratio of each term appearing in the equations of the exact factorization. We observe how non-adiabatic effects induced by the coupling to the nuclear motion affect electronic properties and we analyze the leading term, connecting it to the classical nuclear momentum. Its dependence on the mass ratio is tested numerically on a model of proton-coupled electron transfer in different non-adiabatic regimes.