Electron flux in molecules induced by nuclear motion

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.

Spin current in the early stage of radical reactions and its mechanisms

We study the electronic spin flux (atomic-scale flow of the spin density in molecules) by a perturbation analysis and ab initio nonadiabatic calculations. We derive a general perturbative expression of the charge and spin fluxes and identify the driving perturbation of the fluxes to be the time derivative of the electron-nucleus interaction term in the Hamiltonian. We then expand the expression in molecular orbitals so as to identify relevant components of the fluxes. Our perturbation theory describes the electronic fluxes in the early stage of reactions in an intuitively clear manner. The perturbation theory is then applied to an analysis of the spin flux obtained in ab initio calculations of the radical reaction of O2 and CH3· starting from three distinct spin configurations; (a) CH3· and triplet O2 with total spin of the system set Stot=1/2 (b) CH3· and singlet O2, Stot=1/2, and (c) CH3· and triplet O2, Stot=3/2. Further analysis of the time-dependent behaviors of the spin flux in these numerical simulations reveals (i) the spin flux induces rearrangement of the local spin structure, such as reduction of the spin polarization arising from the triplet O2 and (ii) the spin flux flows from O2 to CH3· in the reaction starting from spin configuration (a) and from CH3· to O2 in that starting from configuration (b), whereas no major intermolecular spin flux was observed in that starting from configuration (c). Our study thus establishes the mechanism of the spin flux that rearranges the local spin structures associated with chemical bonds.

Electronic-state chaos, intramolecular electronic energy redistribution, and chemical bonding in persisting multidimensional nonadiabatic systems

We study the chaotic, huge fluctuation of electronic state, resultant intramolecular energy redistribution, and strong chemical bonding surviving the fluctuation with exceedingly long lifetimes of highly excited boron clusters. Those excited states constitute densely quasi-degenerate state manifolds. The huge fluctuation is induced by persisting multidimensional nonadiabatic transitions among the states in the manifold. We clarify the mechanism of their coexistence and its physical significance. In doing so, we concentrate on two theoretical aspects. One is quantum chaos and energy randomization, which are to be directly extracted from the properties of the total electronic wavefunctions. The present dynamical chaos takes place through frequent transitions from adiabatic states to others, thereby making it very rare for the system to find dissociation channels. This phenomenon leads to the concept of what we call intramolecular nonadiabatic electronic-energy redistribution, which is an electronic-state generaliztion of the notion of intramolecular vibrational energy redistribution. The other aspect is about the peculiar chemical bonding. We investigate it with the energy natural orbitals (ENOs) to see what kind of theoretical structures lie behind the huge fluctuation. The ENO energy levels representing the highly excited states under study appear to have four robust layers. We show that the energy layers responsible for chaotic dynamics and those for chemical bonding are widely separated from each other, and only when an event of what we call "inter-layer crossing" happens to burst can the destruction of these robust energy layers occur, resulting in molecular dissociation. This crossing event happens only rarely because of the large energy gaps between the ENO layers. It is shown that the layers of high energy composed of complex-valued ENOs induce the turbulent flow of electrons and electronic-energy in the cluster. In addition, the random and fast time-oscillations of those high energy ENOs serve as a random force on the nuclear dynamics, which can work to prevent a concentration of high nuclear kinetic energy in the dissociation channels.

Energy natural orbital characterization of nonadiabatic electron wavepackets in the densely quasi-degenerate electronic state manifold

Dynamics and energetic structure of largely fluctuating nonadiabatic electron wavepackets are studied in terms of Energy Natural Orbitals (ENOs) [K. Takatsuka and Y. Arasaki, J. Chem. Phys. 154, 094103 (2021)]. Such huge fluctuating states are sampled from the highly excited states of clusters of 12 boron atoms (B₁₂), which have densely quasi-degenerate electronic excited-state manifold, each adiabatic state of which gets promptly mixed with other states through the frequent and enduring nonadiabatic interactions within the manifold. Yet, the wavepacket states are expected to be of very long lifetimes. This excited-state electronic wavepacket dynamics is extremely interesting but very hard to analyze since they are usually represented in large time-dependent configuration interaction wavefunctions and/or in some other complicated forms. We have found that ENO gives an invariant energy orbital picture to characterize not only the static highly correlated electronic wavefunctions but also those time-dependent electronic wavefunctions. Hence, we first demonstrate how the ENO representation works for some general cases, choosing proton transfer in water dimer and electron-deficient multicenter chemical bonding in diborane in the ground state. We then penetrate with ENO deep into the analysis of the essential nature of nonadiabatic electron wavepacket dynamics in the excited states and show the mechanism of the coexistence of huge electronic fluctuation and rather strong chemical bonds under very random electron flows within the molecule. To quantify the intra-molecular energy flow associated with the huge electronic-state fluctuation, we define and numerically demonstrate what we call the electronic energy flux.

Quantum Chaos in the Dynamics of Molecules

Quantum chaos is reviewed from the viewpoint of "what is molecule?", particularly placing emphasis on their dynamics. Molecules are composed of heavy nuclei and light electrons, and thereby the very basic molecular theory due to Born and Oppenheimer gives a view that quantum electronic states provide potential functions working on nuclei, which in turn are often treated classically or semiclassically. Therefore, the classic study of chaos in molecular science began with those nuclear dynamics particularly about the vibrational energy randomization within a molecule. Statistical laws in probabilities and rates of chemical reactions even for small molecules of several atoms are among the chemical phenomena requiring the notion of chaos. Particularly the dynamics behind unimolecular decomposition are referred to as Intra-molecular Vibrational energy Redistribution (IVR). Semiclassical mechanics is also one of the main research fields of quantum chaos. We herein demonstrate chaos that appears only in semiclassical and full quantum dynamics. A fundamental phenomenon possibly giving birth to quantum chaos is "bifurcation and merging" of quantum wavepackets, rather than "stretching and folding" of the baker's transformation and the horseshoe map as a geometrical foundation of classical chaos. Such wavepacket bifurcation and merging are indeed experimentally measurable as we showed before in the series of studies on real-time probing of nonadiabatic chemical reactions. After tracking these aspects of molecular chaos, we will explore quantum chaos found in nonadiabatic electron wavepacket dynamics, which emerges in the realm far beyond the Born-Oppenheimer paradigm. In this class of chaos, we propose a notion of Intra-molecular Nonadiabatic Electronic Energy Redistribution (INEER), which is a consequence of the chaotic fluxes of electrons and energy within a molecule.

Real-time electronic energy current and quantum energy flux in molecules

Intra- and inter-molecular electronic energy current is formulated by defining the probability current of electronic energy, called the energy flux. Among vast possible applications to electronic energy transfer phenomena, including chemical reaction dynamics, here we present a first numerical example from highly excited nonadiabatic electron wavepacket dynamics of a boron cluster B₁₂.

Nature of chemical bond and potential barrier in an invariant energy-orbital picture

Physical nature of the chemical bond and potential barrier is studied in terms of energy natural orbitals (ENOs), which are extracted from highly correlated electronic wavefunctions. ENO provides an objective one-electron picture about energy distribution in a molecule, just as the natural orbitals (NOs) represent one electron view about electronic charge distribution. ENO is invariant in the same sense as NO is, that is, ENOs converge to the exact ones as a series of approximate wavefunctions approach the exact one, no matter how the methods of approximation are adopted. Energy distribution analysis based on ENO can give novel insights about the nature of chemical bonding and formation of potential barriers, besides information based on the charge distribution alone. With ENOs extracted from full configuration interaction wavefunctions in a finite yet large enough basis set, we analyze the nature of chemical bonding of three low-lying electronic states of a hydrogen molecule, all being in different classes of the so-called covalent bond. The mechanism of energy lowering in bond formation, which gives a binding energy, is important, yet not the only concern for this small molecule. Another key notion in chemical bonding is whether a potential basin is well generated stiff enough to support a vibrational state(s) on it. Based on the virial theorem in the adiabatic approximation and in terms of the energy and force analyses with ENOs, we study the roles of the electronic kinetic energy and its nuclear derivative(s) on how they determine the curvature (or the force constant) of the potential basins. The same idea is applied to the potential barrier and, thereby, the transition states. The rate constant within the transition-state theory is formally shown to be described in terms of the electronic kinetic energy and the nuclear derivatives only. Thus, the chemical bonding and rate process are interconnected behind the scenes. Besides this aspect, we pay attention to (1) the effects of electron correlation that manifests itself not only in the orbital energy but also in the population of ENOs and (2) the role of nonadiabaticity (diabatic state mixing), resulting in double basins and a barrier on a single potential curve in bond formation. These factors differentiate a covalent bond into subclasses.

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.

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.

Charge separation and successive reconfigurations of electronic and protonic states in a water-splitting catalytic cycle with the Mn4CaO5 cluster. On the mechanism of water splitting in PSII

Charge separation, reloading of electrons and protons, and O₂ generation in a catalytic cycle for water splitting with Mn₄CaO₅ in PSII.

On the Elementary Chemical Mechanisms of Unidirectional Proton Transfers: A Nonadiabatic Electron-Wavepacket Dynamics Study

We propose a set of chemical reaction mechanisms of unidirectional proton transfers, which may possibly work as an elementary process in chemical and biological systems. Being theoretically derived based on our series of studies on charge separation dynamics in water splitting by Mn oxides, the present mechanisms have been constructed after careful exploration over the accumulated biological studies on cytochrome c oxidase (CcO) and bacteriorhodopsin. In particular, we have focused on the biochemical findings in the literature that unidirectional transfers of approximately two protons are driven by one electron passage through the reaction center (binuclear center) in CcO, whereas no such dissipative electron transfer is believed to be demanded in the proton transport in bacteriorhodopsin. The proposed basic mechanisms of unidirectional proton transfers are further reduced to two elementary dynamical processes, namely, what we call the coupled proton and electron-wavepacket transfer (CPEWT) and the inverse CPEWT. To show that the proposed mechanisms can indeed be materialized in a molecular level, we construct model systems with possible molecules that are rather familiar in biological chemistry, for which we perform the ab initio calculations of full-dimensional nonadiabatic electron-wavepacket dynamics coupled with all nuclear motions including proton transfers.

Na2 Vibrating in the Double-Well Potential of State 2 1Σu+ (JM = 00): A Pulsating "Quantum Bubble" with Antagonistic Electronic Flux

The theory of concerted electronic and nuclear flux densities associated with the vibration and dissociation of a multielectron nonrotating homonuclear diatomic molecule (or ion) in an electronic state ^2S+1Σ_g,u⁺ (JM = 00) is presented. The electronic population density, nuclear probability density, and nuclear flux density are isotropic. A theorem of Barth , presented in this issue, shows that the electronic flux density (EFD) is also isotropic. Hence, the evolving system appears as a pulsating, or exploding, "quantum bubble". Application of the theory to Na₂ vibrating in the double-minimum potential of the 2  ¹Σ_u⁺ (JM = 00) excited state reveals that the EFD consists of two antagonistic components. One arises from electrons that flow essentially coherently with the nuclei. The other, which is oppositely directed (i.e., antagonistic) and more intense, is due to the transition in electronic structure from "Rydberg" to "ionic" type as the nuclei traverse the potential barrier between inner and outer potential wells. This "transition" component of the EFD rises and falls sharply as the nuclei cross the barrier.

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.

Concerted Electronic and Nuclear Fluxes During Coherent Tunnelling in Asymmetric Double-Well Potentials

The quantum theory of concerted electronic and nuclear fluxes (CENFs) during coherent periodic tunnelling from reactants (R) to products (P) and back to R in molecules with asymmetric double-well potentials is developed. The results are deduced from the solution of the time-dependent Schrödinger equation as a coherent superposition of two eigenstates; here, these are the two states of the lowest tunnelling doublet. This allows the periodic time evolutions of the resulting electronic and nuclear probability densities (EPDs and NPDs) as well as the CENFs to be expressed in terms of simple sinusodial functions. These analytical results reveal various phenomena during coherent tunnelling in asymmetric double-well potentials, e.g., all EPDs and NPDs as well as all CENFs are synchronous. Distortion of the symmetric reference to a system with an asymmetric double-well potential breaks the spatial symmetry of the EPDs and NPDs, but, surprisingly, the symmetry of the CENFs is conserved. Exemplary application to the Cope rearrangement of semibullvalene shows that tunnelling of the ideal symmetric system can be suppressed by asymmetries induced by rather small external electric fields. The amplitude for the half tunnelling, half nontunnelling border is as low as 0.218 × 10(-8) V/cm. At the same time, the delocalized eigenstates of the symmetric reference, which can be regarded as Schrödinger's cat-type states representing R and P with equal probabilities, get localized at one or the other minima of the asymmetric double-well potential, representing either R or P.

Staircase patterns of nuclear fluxes during coherent tunneling in excited doublets of symmetric double well potentials

Step-by-step flux for one-by-one transfers of the lobes of the density, from the reactant (left) to the product (right) in the excited tunneling doublet.

Imaging the Ultrafast Photoelectron Transfer Process in Alizarin-TiO2

In this work, we adopt a quantum mechanical approach based on time-dependent density functional theory (TDDFT) to study the optical and electronic properties of alizarin supported on TiO2 nano-crystallites, as a prototypical dye-sensitized solar cell. To ensure proper alignment of the donor (alizarin) and acceptor (TiO2 nano-crystallite) levels, static optical excitation spectra are simulated using time-dependent density functional theory in response. The ultrafast photoelectron transfer from the dye to the cluster is simulated using an explicitly time-dependent, one-electron TDDFT ansatz. The model considers the δ-pulse excitation of a single active electron localized in the dye to the complete set of energetically accessible, delocalized molecular orbitals of the dye/nano-crystallite complex. A set of quantum mechanical tools derived from the transition electronic flux density is introduced to visualize and analyze the process in real time. The evolution of the created wave packet subject to absorbing boundary conditions at the borders of the cluster reveal that, while the electrons of the aromatic rings of alizarin are heavily involved in an ultrafast charge redistribution between the carbonyl groups of the dye molecule, they do not contribute positively to the electron injection and, overall, they delay the process.

Quantum theory of concerted electronic and nuclear fluxes associated with adiabatic intramolecular processes

Example of concerted electronic (right) and nuclear (left) fluxes: isomerization of B₄.