Electronic Flux Density beyond the Born–Oppenheimer Approximation

Practical phase-space electronic Hamiltonians for ab initio dynamics

Modern electronic structure theory is built around the Born-Oppenheimer approximation and the construction of an electronic Hamiltonian Ĥel(X) that depends on the nuclear position X (and not the nuclear momentum P). In this article, using the well-known theory of electron translation (Γ') and rotational (Γ″) factors to couple electronic transitions to nuclear motion, we construct a practical phase-space electronic Hamiltonian that depends on both nuclear position and momentum, ĤPS(X,P). While classical Born-Oppenheimer dynamics that run along the eigensurfaces of the operator Ĥel(X) can recover many nuclear properties correctly, we present some evidence that motion along the eigensurfaces of ĤPS(X,P) can better capture both nuclear and electronic properties (including the elusive electronic momentum studied by Nafie). Moreover, only the latter (as opposed to the former) conserves the total linear and angular momentum in general.

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.

Representation and conservation of angular momentum in the Born-Oppenheimer theory of polyatomic molecules

This paper concerns the representation of angular momentum operators in the Born-Oppenheimer theory of polyatomic molecules and the various forms of the associated conservation laws. Topics addressed include the question of whether these conservation laws are exactly equivalent or only to some order of the Born-Oppenheimer parameter κ = (m/M)^1/4 and what the correlation is between angular momentum quantum numbers in the various representations. These questions are addressed in both problems involving a single potential energy surface and those with multiple, strongly coupled surfaces and in both the electrostatic model and those for which fine structure and electron spin are important. The analysis leads to an examination of the transformation laws under rotations of the electronic Hamiltonian; of the basis states, both adiabatic and diabatic, along with their phase conventions; of the potential energy matrix; and of the derivative couplings. These transformation laws are placed in the geometrical context of the structures in the nuclear configuration space that are induced by rotations, which include the rotational orbits or fibers, the surfaces upon which the orientation of the molecule changes but not its shape, and the section, an initial value surface that cuts transversally through the fibers. Finally, it is suggested that the usual Born-Oppenheimer approximation can be replaced by a dressing transformation, that is, a sequence of unitary transformations that block-diagonalize the Hamiltonian. When the dressing transformation is carried out, we find that the angular momentum operator does not change. This is a part of a system of exact equivalences among various representations of angular momentum operators in Born-Oppenheimer theory. Our analysis accommodates large-amplitude motions and is not dependent on small-amplitude expansions about an equilibrium position. Our analysis applies to noncollinear configurations of a polyatomic molecule; this covers all but a subset of measure zero (the collinear configurations) in the nuclear configuration space.

Exact Factorization Adventures: A Promising Approach for Non-Bound States

Modeling the dynamics of non-bound states in molecules requires an accurate description of how electronic motion affects nuclear motion and vice-versa. The exact factorization (XF) approach offers a unique perspective, in that it provides potentials that act on the nuclear subsystem or electronic subsystem, which contain the effects of the coupling to the other subsystem in an exact way. We briefly review the various applications of the XF idea in different realms, and how features of these potentials aid in the interpretation of two different laser-driven dissociation mechanisms. We present a detailed study of the different ways the coupling terms in recently-developed XF-based mixed quantum-classical approximations are evaluated, where either truly coupled trajectories, or auxiliary trajectories that mimic the coupling are used, and discuss their effect in both a surface-hopping framework as well as the rigorously-derived coupled-trajectory mixed quantum-classical approach.

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'.

Adiabatic electronic flux in molecules and in condensed matter

The theory of adiabatic electron transport in a correlated condensed-matter system is rooted in a seminal paper by Niu and Thouless [J. Phys A {\bf 17}, 2453 (1984)]; I adopt here an analogous logic in order to { retrieve the known expression for the adiabatic electronic flux in a molecular system [L. A. Nafie, J. Chem. Phys. {\bf 79}, 4950 (1983)]. Its derivation here is considerably simpler than those available in the current quantum-chemistry literature; it also explicitly identifies the adiabaticity parameter, in terms of which the adiabatic flux and the electron density are both exact to first order. It is shown that continuity equation is conserved to the same order. For the sake of completeness,} I also briefly outline the relevance of the macroscopic electronic flux to the physics of solids and liquids.

Quantum-classical nonadiabatic dynamics of Floquet driven systems

We develop a trajectory-based approach for excited-state molecular dynamics simulations of systems subject to an external periodic drive. We combine the exact-factorization formalism, allowing us to treat electron-nuclear systems in nonadiabatic regimes, with the Floquet formalism for time-periodic processes. The theory is developed starting with the molecular time-dependent Schrödinger equation with the inclusion of an external periodic drive that couples to the system dipole moment. With the support of the Floquet formalism, quantum dynamics is approximated by combining classical-like, trajectory-based, nuclear evolution with electronic dynamics represented in the Floquet basis. The resulting algorithm, which is an extension of the coupled-trajectory mixed quantum-classical scheme for periodically driven systems, is applied to a model study, exactly solvable, with different field intensities.

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.

Relaxation dynamics through a conical intersection: Quantum and quantum-classical studies

We study the relaxation process through a conical intersection of a photo-excited retinal chromophore model. The analysis is based on a two-electronic-state two-dimensional Hamiltonian developed by Hahn and Stock [J. Phys. Chem. B 104 1146 (2000)] to reproduce, with a minimal model, the main features of the 11-cis to all-trans isomerization of the retinal of rhodopsin. In particular, we focus on the performance of various trajectory-based schemes to nonadiabatic dynamics, and we compare quantum-classical results to the numerically exact quantum vibronic wavepacket dynamics. The purpose of this work is to investigate, by analyzing electronic and nuclear observables, how the sampling of initial conditions for the trajectories affects the subsequent dynamics.

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

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.

Probing Electronic Fluxes via Time-Resolved X-Ray Scattering

The current flux density is a vector field that can be used to describe theoretically how electrons flow in a system out of equilibrium. In this work, we unequivocally demonstrate that the signal obtained from time-resolved x-ray scattering does not only map the time evolution of the electronic charge distribution, but also encodes information about the associated electronic current flux density. We show how the electronic current flux density qualitatively maps the distribution of electronic momenta and reveals the underlying mechanism of ultrafast charge migration processes, while also providing quantitative information about the timescales of electronic coherences.

On the Probability Density of the Nuclei in a Vibrationally Excited Molecule

For localized and oriented vibrationally excited molecules, the qualitative features of the one-body probability density of the nuclei (one-nucleus density) are investigated. Like the familiar and widely used one-electron density that represents the probability of finding an electron at a given location in space, the one-nucleus density represents the probability of finding a nucleus at a given position in space independent of the location of the other nuclei and independent of their type. In contrast to the electrons, however, the nuclei are comparably localized. Due to this localization of the individual nuclei, the one-nucleus density provides a quantum-mechanical representation of the "chemical picture" of the molecule as an object that can largely be understood in a three-dimensional space, even though its full nuclear probability density is defined on the high-dimensional configuration space of all the nuclei. We study how the nodal structure of the wavefunctions of vibrationally excited states translates to the one-nucleus density. It is found that nodes do not necessarily lead to visible changes in the one-nucleus density: Already for relatively small molecules, only certain vibrational excitations change the one-nucleus density qualitatively compared to the ground state. It turns out that there are simple rules for predicting the shape of the one-nucleus density from the normal mode coordinates. A Python module for the computation of the one-nucleus density is provided at https://gitlab.com/axelschild/mQNMc.

On the numerical solution of the exact factorization equations

The exact factorization (EF) approach to coupled electron-ion dynamics recasts the time-dependent molecular Schrödinger equation as two coupled equations, one for the nuclear wavefunction and one for the conditional electronic wavefunction. The potentials appearing in these equations have provided insight into non-adiabatic processes, and new practical non-adiabatic dynamics methods have been formulated starting from these equations. Here, we provide a first demonstration of a self-consistent solution of the exact equations, with a preliminary analysis of their stability and convergence properties. The equations have an unprecedented mathematical form, involving a Hamiltonian outside the class of Hermitian Hamiltonians usually encountered in time-propagation, and so the usual numerical methods for time-dependent Schrödinger fail when applied in a straightforward way to the EF equations. We find an approach that enables stable propagation long enough to witness non-adiabatic behavior in a model system before non-trivial instabilities take over. Implications for the development and analysis of EF-based methods are discussed.

Curly arrows, electron flow, and reaction mechanisms from the perspective of the bonding evolution theory

Despite the usefulness of curly arrows in chemistry, their relationship with real electron density flows is still imprecise, and even their direct connection to quantum chemistry is still controversial. The paradigmatic description - from first principles - of the mechanistic aspects of a given chemical process is based mainly on the relative energies and geometrical changes at the stationary points of the potential energy surface along the reaction pathway; however, it is not sufficient to describe chemical systems in terms of bonding aspects. Probing the electron density distribution during a chemical reaction can provide important insights, enabling us to understand and control chemical reactions. This aim has required an extension of the relationships between the concepts of traditional chemistry and those of quantum mechanics. Bonding evolution theory (BET), which combines the topological analysis of the electron localization function (ELF) and Thom's catastrophe theory (CT), provides a powerful method that offers insight into the molecular mechanism of chemical rearrangements. In agreement with the laws of physical and aspects of quantum theory, BET can be considered an appropriate tool to tackle chemical reactivity with a wide range of possible applications. In this work, BET is applied to address a long-standing problem: the ability to monitor the flow of electron density. BET analysis shows a connection between quantum mechanics and bond making/forming processes. Likewise, the present approach retrieves the classical curly arrows used to describe the rearrangements of chemical bonds and provides detailed physical grounds for this type of representation. We demonstrate this procedure using the test set of prototypical examples of thermal ring apertures, and the degenerated Cope rearrangement of semibullvalene.

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.