Search

There are different ways to obtain an exact one-electron theory for a many-electron system, and the exact electron factorization (EEF) is one of them. In the EEF, the Schr\"odinger equation for one electron in the environment of other electrons is constructed. The environment provides the potentials that appear in this equation: A scalar potential $v^{\rm H}$ representing the energy of the environment and another scalar potential $v^{\rm G}$ as well as a vector potential that have geometric meaning. By replacing the interacting many-electron system with the non-interacting Kohn-Sham (KS) system, we show how the EEF is related to density functional theory (DFT) and we interpret the Hartree-exchange-correlation potential as well as the Pauli potential in terms of the EEF. In particular, we show that from the EEF viewpoint, the Pauli potential does not represent the difference between a fermionic and a bosonic non-interacting system, but that it corresponds to $v^{\rm G}$ and partly to $v^{\rm H}$ for the (fermionic) KS system. We then study the meaning of $v^{\rm G}$ in detail: Its geometric origin as a metric measuring the change of the environment is presented. Additionally, its behavior for a simple model of a homo- and heteronucler diatomic is investigated and interpreted with the help of a two-state model. In this way, we provide a physical interpretation for the one-electron potentials that appear in the EEF and in DFT., Comment: arXiv admin note: substantial text overlap with arXiv:2010.14885

Quantum trajectory calculations for electrons are a useful tool in the field of molecular dynamics, e.g. to understand processes in ultrafast spectroscopy. They have, however, two limitation: On the one hand, such calculations are typically based on the Born-Oppenheimer approximation (BOA) and the electron dynamics for stationary nuclei is considered, thus neglecting quantum effects of the nuclei. On the other hand, even if the quantum nuclear motion would be taken into account, a BOA dynamics on a single potential energy surface would not provide any electron trajectories because the electronic part is treated as a stationary problem. By using the exact factorization method, we overcome these limitations and generalize the theory of electronic quantum trajectories to a fully quantum-mechanical treatment of the nuclei. After reviewing the time-dependent theory of quantum hydrodynamics and quantum trajectories, we show that the nuclei can be viewed as a quantum clock for the electronic motion and we develop a fully quantum-mechanical clock-dependent version of quantum hydrodynamics. This theory is used to obtain electronic trajectories for quantum nuclei, as is exemplified for a model system of a proton-coupled electron transfer dynamics. Our work generalizes the concept of quantum trajectories and lays the foundations for the development of trajectory-based simulation methods of electron dynamics beyond the BOA., Comment: arXiv admin note: substantial text overlap with arXiv:1911.13033

When a molecule dissociates, the exact Kohn-Sham (KS) and Pauli potentials may form step structures. Reproducing these steps correctly is central for the description of dissociation and charge-transfer processes in density functional theory (DFT): The steps align the KS eigenvalues of the dissociating subsystems relative to each other and determine where electrons localize. While the step height can be calculated from the asymptotic behavior of the KS orbitals, this provides limited insight into what causes the steps. We give an explanation of the steps with an exact mapping of the many-electron problem to a one-electron problem, the exact electron factorizaton (EEF). The potentials appearing in the EEF have a clear physical meaning that translates to the DFT potentials by replacing the interacting many-electron system with the KS system. With a simple model of a diatomic, we illustrate that the steps are a consequence of spatial electron entanglement and are the result of a charge transfer. From this mechanism, the step height can immediately be deduced. Moreover, two methods to approximately reproduce the potentials during dissociation suggest themselves. One is based on the states of the dissociated system, while the other one is based on an analogy to the Born-Oppenheimer treatment of a molecule. The latter method also shows that the steps connect adiabatic potential energy surfaces. The view of DFT from the EEF thus provides a better understanding of how many-electron effects are encoded in a one-electron theory and how they can be modeled.

The theoretical and computational description of materials properties is a task of utmost scientific and technological importance. A first-principles description of electron-electron interactions poses an immense challenge that is usually approached by converting the many-electron problem to an effective one-electron problem. There are different ways to obtain an exact one-electron theory for a many-electron system. An emergent method is the exact electron factorization (EEF) -- one of the branches of the Exact Factorization approach to many-body systems. In the EEF, the Schr\"odinger equation for one electron, in the environment of all other electrons, is formulated. The influence of the environment is reflected in the potential $v^{\rm H}$, which represents the energy of the environment, and in a potential $v^{\rm G}$, which has a geometrical meaning. In this paper, we focus on $v^{\rm G}$ and study its properties in detail. We investigate the geometric origin of $v^{\rm G}$ as a metric measuring the change of the environment, exemplify how translation and scaling of the state of the environment are reflected in $v^{\rm G}$, and explain its shape for homo- and heteronuclear diatomic model systems. Based on the close connection between the EEF and density functional theory, we also use $v^{\rm G}$ to provide an alternative interpretation to the Pauli potential in orbital-free density functional theory.

Attosecond interferometry (AI) is an experimental technique based on ionizing a system with an attosecond pulse train in the presence of an assisting laser. This assisting laser provides multiple pathways for the photoelectron wave packet to reach the same final state, and interference of these pathways can be used to probe properties of matter. The mechanism of AI is well-understood for isolated atoms and molecules in the gas phase, but not so much in the condensed phase, especially if the substrate under study is transparent. Then, additional pathways open up for the electron due to scattering from neighbouring atoms. We investigate to what extent these additional pathways influence the measured photoionization delay with the help of one- and three-dimensional model systems. In both cases, we find that the total delay can be expressed as the sum of a local (photoionization) delay and a non-local delay which contains the effect of electron scattering during transport. The 1D system shows that the non-local delay is an oscillatory function of the distance between the sites where ionization and scattering take place. A similar result is obtained in 3D, but the modulation depth of the non-local delay is found to strongly depend on the effective scattering cross section. We conclude that attosecond interferometry of disordered systems like liquids at low photon energies (20-30 eV) is mainly sensitive to the local delay, i.e., to changes of the photoionization dynamics induced by the immediate environment of the ionized entity, and less to electron scattering during transport through the medium.

If one-electron observables of a many-electron system are of interest, a many-electron dynamics can be represented exactly by a one-electron dynamics with effective potentials. The formalism for this reduction is provided by the Exact Electron Factorization (EEF). We study the time-dependent features of the EEF effective potentials for a model of an atom ionized by an ultrastrong and ultrashort laser pulse, with the aim of understanding what is needed to develop computationally feasible approximations. It is found that the simplest approximation, the so-called time-independent conditional amplitude (TICA) approximation, is complementary to single-active electron (SAE) approaches as it reproduced the exact dynamics well for high photon frequencies of the laser field or large Keldysh parameter. For relatively low frequencies of the laser field or for smaller Keldysh parameters, we find that excited state dynamics in the core region of the atom leads to a time-dependent ionization barrier in the EEF potential. The time-dependence of the barrier needs to be described accurately to correctly model many-electron effects, and we conclude that a multi-state extension of the TICA approximation is a possible route how this can be achieved. In general, our study sheds a different light on one-electron pictures of strong-field ionization and shows that many-electron effects for such processes may be included by solving a one-electron Schr\"odinger equation, provided the core dynamics can be modeled successfully.

In a recent comment, Ruth Signorell raises a number of issues that she considers to question the validity of our approach to determine mean free paths for electron scattering in liquid water and our comparison with the results on amorphous ice by Michaud, Wen, and Sanche. Here, we show that these critiques are unjustified, being either unfounded or based on misconceptions by the author of the comment. We nevertheless welcome the opportunity to further clarify certain aspects of our work that we did not discuss in detail in our letter., Comment: The comment which is reply is referring to was rejected for publication, hence this reply will also not be published

The mean free paths of low-energy electrons in liquid water are of fundamental importance for modelling radiation damage and many related physico-chemical processes. Neither theoretical predictions nor experimental estimations have so far converged to yield reliable values for these parameters. We therefore introduce a new approach to determine the elastic and inelastic mean free paths (EMFP, IMFP) based on experimental data. We report extensive ab-initio calculations of electron quantum scattering with water clusters, which are brought to convergence with respect to the cluster size. This provides a first-principles approach to condensed-phase scattering that includes both multiple-scattering and condensation effects. The obtained differential cross sections are used in a detailed Monte-Carlo simulation to extract EMFP and IMFP from two recent liquid-microjet experiments that determined the effective attenuation length (EAL) and the photoelectron angular distribution (PAD) following oxygen 1s-ionization of liquid water. For electron kinetic energies from 10 eV to 300 eV, we find that the IMFP is noticeably larger than the EAL. The EMFP is longer than that of gas-phase water and the IMFP is longer compared to the latest theoretical estimations, but both the EMFP and IMFP are much shorter than suggested by experimental results for amorphous ice. The Python module developed for the analysis is available at https://gitlab.com/axelschild/CLstunfti and can be used to further refine our results when new experimental data become available.

Quantum hydrodynamics is a formulation of quantum mechanics based on the probability density and flux (current) density of a quantum system. It can be used to define trajectories which allow for a particle-based interpretation of quantum mechanics, commonly known as Bohmian mechanics. However, quantum hydrodynamics rests on the usual time-dependent formulation of quantum mechanics where time appears as a parameter. This parameter describes the correlation of the state of the quantum system with an external system -- a clock -- which behaves according to classical mechanics. With the Exact Factorization of a quantum system into a marginal and a conditional system, quantum mechanics and hence quantum hydrodynamics can be generalized for quantum clocks. In this article, the theory is developed and it is shown that trajectories for the quantum system can still be defined, and that these trajectories depend conditionally on the trajectory of the clock. Such trajectories are not only interesting from a fundamental point of view, but they can also find practical applications whenever a dynamics relative to an external time parameter is composed of "fast" and "slow" degrees of freedom and the interest is in the fast ones, while quantum effects of the slow ones (like a branching of the wavepacket) cannot be neglected. As an illustration, time- and clock-dependent trajectories are calculated for a model system of a non-adiabatic dynamics, where an electron is the quantum system, a nucleus is the quantum clock, and an external time parameter is provided, e.g. via an interaction with a laser field that is not treated explicitly.

For localized and oriented vibrationally excited molecules, the one-body probability density of the nuclei (one-nucleus density) is studied. 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. In contrast to the full many-dimensional nuclear probability density, the one-nucleus density contains less information and may thus be better accessible by experiment, especially for large molecules. It also provides a quantum-mechanical view of molecular vibrations that can easily be visualized. We study how the nodal structure of the wavefunctions of vibrationally excited states translates to the one-nucleus density. It is found that nodes are not necessarily visible: 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 some simple rules for predicting the shape of the one-nucleus density from the normal mode coordinates, and thus for predicting if a vibrational excitation is visible in a corresponding experiment.

The local conservation of a physical quantity whose distribution changes with time is mathematically described by the continuity equation. The corresponding time parameter, however, is defined with respect to an idealized classical clock. We consider what happens when this classical time is replaced by a non-relativistic quantum-mechanical description of the clock. From the clock-dependent Schr\"odinger equation (as analogue of the time-dependent Schr\"odinger equation) we derive a continuity equation, where, instead of a time-derivative, an operator occurs that depends on the flux (probability current) density of the clock. This clock-dependent continuity equation can be used to analyze the dynamics of a quantum system and to study degrees of freedom that may be used as internal clocks for an approximate description of the dynamics of the remaining degrees of freedom. As an illustration, we study a simple model for coupled electron-nuclear dynamics and interpret the nuclei as quantum clock for the electronic motion. We find that whenever the Born-Oppenheimer approximation is valid, the continuity equation shows that the nuclei are the only relevant clock for the electrons.

We present an exact single-electron picture that describes the correlated electron dynamics in strong laser fields. Our approach is based on the factorization of the electronic wavefunction as a product of a marginal and a conditional amplitude. The marginal amplitude, which depends only on one electronic coordinate and yields the exact one-electron density and current density, obeys a time-dependent Schr\"odinger equation with an effective time-dependent potential. The exact equations are used to derive an approximation that is a step towards a general and feasible ab-initio single-electron approximation for molecules. The derivation also challenges the usual interpretation of the single-active electron approximation. From the study of model systems, we find that the exact and approximate single-electron potentials for processes with negligible two-electron ionization lead to a qualitatively similar dynamics, but that the ionization barrier may be explicitly time-dependent.

We present a novel ab-initio single-electron approach to correlated electron dynamics in strong laser fields. By writing the electronic wavefunction as a product of a marginal one-electron wavefunction and a conditional wavefunction, we show that the exact harmonic spectrum can be obtained from a single-electron Schr\"odinger equation. To obtain the one-electron potential in practice, we propose an adiabatic approximation, i.e. a potential is generated that depends only on the position of one electron. This potential, together with the laser interaction, is then used to obtain the dynamics of the system. For a model Helium atom in a laser field, we show that by using our approach, the high-order harmonic generation spectrum can be obtained to a good approximation.

ORBKIT is a toolbox for post-processing electronic structure calculations based on a highly modular and portable Python architecture. The program allows computing a multitude of electronic properties of molecular systems on arbitrary spatial grids from the basis set representation of its electronic wavefunction, as well as several grid-independent properties. The required data can be extracted directly from the standard output of a large number of quantum chemical programs. ORBKIT can be used as a standalone program to determine standard quantities, for example, the electron density, molecular orbitals, and derivatives thereof. The cornerstone of ORBKIT is its modular structure. The existing basic functions can be arranged in an individual way and can be easily extended by user-written modules to determine any other desired quantities. ORBKIT offers multiple output formats that can be processed by common visualization tools (VMD, Molden, etc.). Additionally, ORBKIT possesses routines to order molecular orbitals computed at different nuclear configurations according to their electronic character and to interpolate the wavefunction between these configurations. The program is open-source under GNU-LGPLv3 license and freely available at http://sourceforge.net/projects/orbkit/. This article provides an overview of ORBKIT with particular focus on its capabilities and applicability, and includes several example calculations.

We consider coherent tunneling of one-dimensional model systems in non-cyclic or cyclic symmetric double well potentials. Generic potentials are constructed which allow for analytical estimates of the quantum dynamics in the non-relativistic deep tunneling regime, in terms of the tunneling distance, barrier height and mass (or moment of inertia). For cyclic systems, the results may be scaled to agree well with periodic potentials for which semi-analytical results in terms of Mathieu functions exist. Starting from a wavepacket which is initially localized in one of the potential wells, the subsequent periodic tunneling is associated with tunneling velocities. These velocities (or angular velocities) are evaluated as the ratio of the flux densities versus the probability densities. The maximum velocities are found under the top of the barrier where they scale as the square root of the ratio of barrier height and mass (or moment of inertia), independent of the tunneling distance. They are applied exemplarily to several prototypical molecular models of non-cyclic and cyclic tunneling, including ammonia inversion, Cope rearrangment of semibullvalene, torsions of molecular fragments, and rotational tunneling in strong laser fields. Typical maximum velocities and angular velocities are in the order of a few km/s and from 10 to 100 THz for our non-cyclic and cyclic systems, respectively, much faster than time-averaged velocities. Even for the more extreme case of an electron tunneling through a barrier of height of one Hartree, the velocity is only about one percent of the speed of light. Estimates of the corresponding time scales for passing through {the narrow domain just} below the potential barrier are in the domain from 2 to 40 fs, much shorter than the tunneling times.