Communication: On the calculation of time-dependent electron flux within the Born-Oppenheimer approximation: A flux-flux reflection principle

Diagonalizing the Born-Oppenheimer Hamiltonian via Moyal perturbation theory, nonadiabatic corrections, and translational degrees of freedom

This article describes a method for calculating higher order or nonadiabatic corrections in Born-Oppenheimer theory and its interaction with the translational degrees of freedom. The method uses the Wigner-Weyl correspondence to map nuclear operators into functions on the classical phase space and the Moyal star product to represent operator multiplication on those functions. These are explained in the body of the paper. The result is a power series in κ2, where κ = (m/M)1/4 is the usual Born-Oppenheimer parameter. The lowest order term is the usual Born-Oppenheimer approximation, while higher order terms are nonadiabatic corrections. These are needed in calculations of electronic currents, momenta, and densities. The separation of nuclear and electronic degrees of freedom takes place in the context of the exact symmetries (for an isolated molecule) of translations and rotations, and these, especially translations, are explicitly incorporated into our discussion. This article presents an independent derivation of the Moyal expansion in molecular Born-Oppenheimer theory. We show how electronic currents and momenta can be calculated within the framework of Moyal perturbation theory; we derive the transformation laws of the electronic Hamiltonian, the electronic eigenstates, and the derivative couplings under translations; we discuss in detail the rectilinear motion of the molecular center of mass in the Born-Oppenheimer representation; and we show how the elimination of the translational components of the derivative couplings leads to a unitary transformation that has the effect of exactly separating the translational degrees of freedom.

Quantum flux densities for electronic-nuclear motion: exact versus Born-Oppenheimer dynamics

We study the coupled electronic-nuclear dynamics in a model system to compare numerically exact calculations of electronic and nuclear flux densities with those obtained from the Born-Oppenheimer (BO) approximation. Within the adiabatic expansion of the total wave function, we identify the terms which contribute to the flux densities. It is found that only off-diagonal elements that involve the interaction between different electronic states contribute to the electronic flux whereas in the nuclear case the major contribution belongs to the BO electronic state. New flux densities are introduced where in both, the electronic and the nuclear case, the main contribution is contained in the component corresponding to the BO state. As a consequence, they can be determined within the BO approximation, and an excellent agreement with the exact results is found. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.

On the molecular electronic flux: Role of nonadiabaticity and violation of conservation

Analysis of electron flux within and in between molecules is crucial in the study of real-time dynamics of molecular electron wavepacket evolution such as those in attosecond laser chemistry and ultrafast chemical reaction dynamics. We here address two mutually correlated issues on the conservation law of molecular electronic flux, which serves as a key consistency condition for electron dynamics. The first one is about a close relation between "weak" nonadiabaticity and the electron dynamics in low-energy chemical reactions. We show that the electronic flux in adiabatic reactions can be consistently reproduced by taking account of nonadiabaticity. Such nonadiabaticity is usually weak in the sense that it does not have a major effect on nuclear dynamics, whereas it plays an important role in electronic dynamics. Our discussion is based on a nonadiabatic extension of the electronic wavefunction similar in idea to the complete adiabatic formalism developed by Nafie [J. Chem. Phys. 79, 4950 (1983)], which has also recently been reformulated by Patchkovskii [J. Chem. Phys. 137, 084109 (2012)]. We give straightforward proof of the theoretical assertion presented by Nafie using a time-dependent mixed quantum-classical framework and a standard perturbation expansion. Explicitly taking account of the flux conservation, we show that the nonadiabatically induced flux realizes the adiabatic time evolution of the electronic density. In other words, the divergence of the nonadiabatic flux equals the time derivative of the electronic density along an adiabatic time evolution of the target molecule. The second issue is about the accurate computationability of the flux. The calculation of flux needs an accurate representation of the (relative) quantum phase, in addition to the amplitude factor, of a total wavefunction and demands special attention for practical calculations. This paper is the first one to approach this issue directly and show how the difficulties arise explicitly. In doing so, we reveal that a number of widely accepted truncation techniques for static property calculations are potential sources of numerical flux non-conservation. We also theoretically propose alternative strategies to realize better flux conservation.

Time-dependent momentum expectation values from different quantum probability and flux densities

Based on the Ehrenfest theorem, the time-dependent expectation value of a momentum operator can be evaluated equivalently in two ways. The integrals appearing in the expressions are taken over two different functions. In one case, the integrand is the quantum mechanical flux density j̲, and in the other, a different quantity j̲̃ appears, which also has the units of a flux density. The quantum flux density j̲ is related to the probability density ρ via the continuity equation, and j̲̃ may as well be used to define a density ρ̃ that fulfills a continuity equation. Employing a model for the coupled dynamics of an electron and a proton, we document the properties of the densities and flux densities. It is shown that although the mean momentum derived from the two quantities is identical, the various functions exhibit a very different coordinate and time-dependence. In particular, it is found that the flux density j̲̃ directly monitors temporal changes in the probability density, and the density ρ̃ carries information about wave packet dispersion occurring in different spatial directions.