Time-dependent density functional theory for open quantum systems with unitary propagation

Synthesis of Hidden Subgroup Quantum Algorithms and Quantum Chemical Dynamics

We describe a general formalism for quantum dynamics and show how this formalism subsumes several quantum algorithms, including the Deutsch, Deutsch-Jozsa, Bernstein-Vazirani, Simon, and Shor algorithms as well as the conventional approach to quantum dynamics based on tensor networks. The common framework exposes similarities among quantum algorithms and natural quantum phenomena: we illustrate this connection by showing how the correlated behavior of protons in water wire systems that are common in many biological and materials systems parallels the structure of Shor's algorithm.

Shortcut to Self-Consistent Light-Matter Interaction and Realistic Spectra from First Principles

We introduce a simple approach to how an electromagnetic environment can be efficiently embedded into state-of-the-art electronic structure methods, taking the form of radiation-reaction forces. We demonstrate that this self-consistently provides access to radiative emission, natural linewidth, Lamb shifts, strong coupling, electromagnetically induced transparency, Purcell-enhanced and superradiant emission. As an example, we illustrate its seamless integration into time-dependent density-functional theory with virtually no additional cost, presenting a convenient shortcut to light-matter interactions.

An Application of the Madelung Formalism for Dissipating and Decaying Systems

This paper is concerned with the modeling and analysis of quantum dissipation and diffusion phenomena in the Schrödinger picture. We derive and investigate in detail the Schrödinger-type equations accounting for dissipation and diffusion effects. From a mathematical viewpoint, this equation allows one to achieve and analyze all aspects of the quantum dissipative systems, regarding the wave equation, Hamilton–Jacobi and continuity equations. This simplification requires the performance of “the Madelung decomposition” of “the wave function”, which is rigorously attained under the general Lagrangian justification for this modification of quantum mechanics. It is proved that most of the important equations of dissipative quantum physics, such as convection-diffusion, Fokker–Planck and quantum Boltzmann, have a common origin and can be unified in one equation.

Nonlinear light absorption in many-electron systems excited by an instantaneous electric field: a non-perturbative approach

Applications of low-cost non-perturbative approaches in real time, such as time-dependent density functional theory, for the study of nonlinear optical properties of large and complex systems are gaining increasing popularity. However, their assessment still requires the analysis and understanding of elementary dynamical processes in simple model systems. Motivated by the aim of simulating optical nonlinearities in molecules, here exemplified by the case of the quaterthiophene oligomer, we investigate light absorption in many-electron interacting systems beyond the linear regime by using a single broadband impulse of an electric field; i.e. an electrical impulse in the instantaneous limit. We determine non-pertubatively the absorption cross section from the Fourier transform of the time-dependent induced dipole moment, which can be obtained from the time evolution of the wavefunction. We discuss the dependence of the resulting cross section on the magnitude of the impulse and we highlight the advantages of this method in comparison with perturbation theory by working on a one-dimensional model system for which numerically exact solutions are accessible. Thus, we demonstrate that the considered non-pertubative approach provides us with an effective tool for investigating fluence-dependent nonlinear optical excitations.

Beyond the State of the Art: Novel Approaches for Thermal and Electrical Transport in Nanoscale Devices

Almost any interaction between two physical entities can be described through the transfer of either charge, spin, momentum, or energy. Therefore, any theory able to describe these transport phenomena can shed light on a variety of physical, chemical, and biological effects, enriching our understanding of complex, yet fundamental, natural processes, e.g., catalysis or photosynthesis. In this review, we will discuss the standard workhorses for transport in nanoscale devices, namely Boltzmann's equation and Landauer's approach. We will emphasize their strengths, but also analyze their limits, proposing theories and models useful to go beyond the state of the art in the investigation of transport in nanoscale devices.

Time-Dependent Density Functional Theory for Open Systems and Its Applications

Photovoltaic devices, electrochemical cells, catalysis processes, light emitting diodes, scanning tunneling microscopes, molecular electronics, and related devices have one thing in common: open quantum systems where energy and matter are not conserved. Traditionally quantum chemistry is confined to isolated and closed systems, while quantum dissipation theory studies open quantum systems. The key quantity in quantum dissipation theory is the reduced system density matrix. As the reduced system density matrix is an O(M! × M!) matrix, where M is the number of the particles of the system of interest, quantum dissipation theory can only be employed to simulate systems of a few particles or degrees of freedom. It is thus important to combine quantum chemistry and quantum dissipation theory so that realistic open quantum systems can be simulated from first-principles. We have developed a first-principles method to simulate the dynamics of open electronic systems, the time-dependent density functional theory for open systems (TDDFT-OS). Instead of the reduced system density matrix, the key quantity is the reduced single-electron density matrix, which is an N × N matrix where N is the number of the atomic bases of the system of interest. As the dimension of the key quantity is drastically reduced, the TDDFT-OS can thus be used to simulate the dynamics of realistic open electronic systems and efficient numerical algorithms have been developed. As an application, we apply the method to study how quantum interference develops in a molecular transistor in time domain. We include electron-phonon interaction in our simulation and show that quantum interference in the given system is robust against nuclear vibration not only in the steady state but also in the transient dynamics. As another application, by combining TDDFT-OS with Ehrenfest dynamics, we study current-induced dissociation of water molecules under scanning tunneling microscopy and follow its time dependent dynamics. Given the rapid development in ultrafast experiments with atomic resolution in recent years, time dependent simulation of open electronic systems will be useful to gain insight and understanding of such experiments. This Account will mainly focus on the practical aspects of the TDDFT-OS method, describing the numerical implementation and demonstrating the method with applications.

Revealing electronic open quantum systems with subsystem TDDFT

Open quantum systems (OQSs) are perhaps the most realistic systems one can approach through simulations. In recent years, describing OQSs with Density Functional Theory (DFT) has been a prominent avenue of research with most approaches based on a density matrix partitioning in conjunction with an ad-hoc description of system-bath interactions. We propose a different theoretical approach to OQSs based on partitioning of the electron density. Employing the machinery of subsystem DFT (and its time-dependent extension), we provide a novel way of isolating and analyzing the various terms contributing to the coupling between the system and the surrounding bath. To illustrate the theory, we provide numerical simulations on a toy system (a molecular dimer) and on a condensed phase system (solvated excimer). The simulations show that non-Markovian dynamics in the electronic system-bath interactions are important in chemical applications. For instance, we show that the superexchange mechanism of transport in donor-bridge-acceptor systems is a non-Markovian interaction between the donor-acceptor (OQS) with the bridge (bath) which is fully characterized by real-time subsystem time-dependent DFT.

The role of interparticle interaction and environmental coupling in a two-particle open quantum system

The effects of bath coupling on an interacting two-particle quantum system are studied using tools from information theory. Shannon entropies of the one (reduced) and two-particle distribution functions in position, momentum and separable phase-space are examined. Results show that the presence of the bath leads to a delocalization of the distribution functions in position space, and a localization in momentum space. This can be interpreted as a loss of information in position space and a gain of information in momentum space. The entropy sum of the system, in the presence of a bath, is shown to be dependent on the strength of the interparticle potential and also on the strength of the coupling to the bath. The statistical correlation between the particles, and its dependence on the bath and interparticle potential, is examined using mutual information. A stronger repulsive potential between particles, in the presence of the bath, yields a smaller correlation between the particles positions, and a larger one between their momenta.

Dynamical Processes in Open Quantum Systems from a TDDFT Perspective: Resonances and Electron Photoemission

We present a review of different computational methods to describe time-dependent phenomena in open quantum systems and their extension to a density-functional framework. We focus the discussion on electron emission processes in atoms and molecules addressing excited-state lifetimes and dissipative processes. Initially we analyze the concept of an electronic resonance, a central concept in spectroscopy associated with a metastable state from which an electron eventually escapes (electronic lifetime). Resonances play a fundamental role in many time-dependent molecular phenomena but can be rationalized from a time-independent context in terms of scattering states. We introduce the method of complex scaling, which is used to capture resonant states as localized states in the spirit of usual bound-state methods, and work on its extension to static and time-dependent density-functional theory. In a time-dependent setting, complex scaling can be used to describe excitations in the continuum as well as wave packet dynamics leading to electron emission. This process can also be treated by using open boundary conditions which allow time-dependent simulations of emission processes without artificial reflections at the boundaries (i.e., borders of the simulation box). We compare in detail different schemes to implement open boundaries, namely transparent boundaries using Green functions, and absorbing boundaries in the form of complex absorbing potentials and mask functions. The last two are regularly used together with time-dependent density-functional theory to describe the electron emission dynamics of atoms and molecules. Finally, we discuss approaches to the calculation of energy and angle-resolved time-dependent pump-probe photoelectron spectroscopy of molecular systems.

Representing the thermal state in time-dependent density functional theory

Classical molecular dynamics (MD) provides a powerful and widely used approach to determining thermodynamic properties by integrating the classical equations of motion of a system of atoms. Time-Dependent Density Functional Theory (TDDFT) provides a powerful and increasingly useful approach to integrating the quantum equations of motion for a system of electrons. TDDFT efficiently captures the unitary evolution of a many-electron state by mapping the system into a fictitious non-interacting system. In analogy to MD, one could imagine obtaining the thermodynamic properties of an electronic system from a TDDFT simulation in which the electrons are excited from their ground state by a time-dependent potential and then allowed to evolve freely in time while statistical data are captured from periodic snapshots of the system. For a variety of systems (e.g., many metals), the electrons reach an effective state of internal equilibrium due to electron-electron interactions on a time scale that is short compared to electron-phonon equilibration. During the initial time-evolution of such systems following electronic excitation, electron-phonon interactions should be negligible, and therefore, TDDFT should successfully capture the internal thermalization of the electrons. However, it is unclear how TDDFT represents the resulting thermal state. In particular, the thermal state is usually represented in quantum statistical mechanics as a mixed state, while the occupations of the TDDFT wavefunctions are fixed by the initial state in TDDFT. We work to address this puzzle by (A) reformulating quantum statistical mechanics so that thermodynamic expectations can be obtained as an unweighted average over a set of many-body pure states and (B) constructing a family of non-interacting (single determinant) TDDFT states that approximate the required many-body states for the canonical ensemble.

Existence, uniqueness, and construction of the density-potential mapping in time-dependent density-functional theory

In this work we review the mapping from densities to potentials in quantum mechanics, which is the basic building block of time-dependent density-functional theory and the Kohn-Sham construction. We first present detailed conditions such that a mapping from potentials to densities is defined by solving the time-dependent Schrödinger equation. We specifically discuss intricacies connected with the unboundedness of the Hamiltonian and derive the local-force equation. This equation is then used to set up an iterative sequence that determines a potential that generates a specified density via time propagation of an initial state. This fixed-point procedure needs the invertibility of a certain Sturm-Liouville problem, which we discuss for different situations. Based on these considerations we then present a discussion of the famous Runge-Gross theorem which provides a density-potential mapping for time-analytic potentials. Further we give conditions such that the general fixed-point approach is well-defined and converges under certain assumptions. Then the application of such a fixed-point procedure to lattice Hamiltonians is discussed and the numerical realization of the density-potential mapping is shown. We conclude by presenting an extension of the density-potential mapping to include vector-potentials and photons.

Time-dependent density functional theory for open systems with a positivity-preserving decomposition scheme for environment spectral functions

Understanding electronic dynamics on material surfaces is fundamentally important for applications including nanoelectronics, inhomogeneous catalysis, and photovoltaics. Practical approaches based on time-dependent density functional theory for open systems have been developed to characterize the dissipative dynamics of electrons in bulk materials. The accuracy and reliability of such approaches depend critically on how the electronic structure and memory effects of surrounding material environment are accounted for. In this work, we develop a novel squared-Lorentzian decomposition scheme, which preserves the positive semi-definiteness of the environment spectral matrix. The resulting electronic dynamics is guaranteed to be both accurate and convergent even in the long-time limit. The long-time stability of electronic dynamics simulation is thus greatly improved within the current decomposition scheme. The validity and usefulness of our new approach are exemplified via two prototypical model systems: quasi-one-dimensional atomic chains and two-dimensional bilayer graphene.

Determining excitation-energy transfer times and mechanisms from stochastic time-dependent density functional theory

We developed an approach for calculating excitation-energy transfer times in supermolecular arrangements based on stochastic time-dependent density functional theory (STDDFT). The combination of real-time propagation and the stochastic Schrödinger equation with a Kohn-Sham Hamiltonian allows for simulating how an excitation spreads through an assembly of molecular systems. The influence that approximations, such as the dipole-dipole coupling approximation of Förster theory, have on energy-transfer times can be checked explicitly. As a first application of our approach we investigate a light-harvesting-inspired model ring system, calculating the time it takes for an excitation to travel from one side of the ring to the opposite side under ideal and perturbed conditions. Among other things we find that completely removing a molecule from the ring may inhibit energy transfer less than having an energetically detuned molecule in the ring. In addition, Förster's dipole coupling approximation may noticeably overestimate excitation-energy transfer efficiency.

Force-field functor theory: classical force-fields which reproduce equilibrium quantum distributions

Feynman and Hibbs were the first to variationally determine an effective potential whose associated classical canonical ensemble approximates the exact quantum partition function. We examine the existence of a map between the local potential and an effective classical potential which matches the exact quantum equilibrium density and partition function. The usefulness of such a mapping rests in its ability to readily improve Born-Oppenheimer potentials for use with classical sampling. We show that such a map is unique and must exist. To explore the feasibility of using this result to improve classical molecular mechanics, we numerically produce a map from a library of randomly generated one-dimensional potential/effective potential pairs then evaluate its performance on independent test problems. We also apply the map to simulate liquid para-hydrogen, finding that the resulting radial pair distribution functions agree well with path integral Monte Carlo simulations. The surprising accessibility and transferability of the technique suggest a quantitative route to adapting Born-Oppenheimer potentials, with a motivation similar in spirit to the powerful ideas and approximations of density functional theory.

Dividing line between quantum and classical trajectories in a measurement problem: Bohmian time constant

This Letter proposes an answer to a challenge posed by Bell on the lack of clarity in regards to the dividing line between the quantum and classical regimes in a measurement problem. To this end, a generalized logarithmic nonlinear Schrödinger equation is proposed to describe the time evolution of a quantum dissipative system under continuous measurement. Within the Bohmian mechanics framework, a solution to this equation reveals a novel result: it displays a time constant that should represent the dividing line between the quantum and classical trajectories. It is shown that continuous measurements and damping not only disturb the particle but compel the system to converge in time to a Newtonian regime. While the width of the wave packet may reach a stationary regime, its quantum trajectories converge exponentially in time to classical trajectories. In particular, it is shown that damping tends to suppress further quantum effects on a time scale shorter than the relaxation time of the system. If the initial wave packet width is taken to be equal to 2.8×10(-15) m (the approximate size of an electron), the Bohmian time constant is found to have an upper limit, i.e., τ(Bmax)=10(-26) s.

Reduced purities as measures of decoherence in many-electron systems

A hierarchy of measures of decoherence for many-electron systems that is based on the purity and the hierarchy of reduced electronic density matrices is presented. These reduced purities can be used to characterize electronic decoherence in the common case when the many-body electronic density matrix is not known and only reduced information about the electronic subsystem is available. Being defined from reduced electronic quantities, the interpretation of the reduced purities is more intricate than the usual (many-body) purity. This is because the nonidempotency of the r-body reduced electronic density matrix that is the basis of the reduced purity measures can arise due to decoherence or due to electronic correlations. To guide the interpretation, explicit expressions are provided for the one-body and two-body reduced purities for a general electronic state. Using them, the information content and structure of the one-body and two-body reduced purities is established, and limits on the changes that decoherence can induce are elucidated. The practical use of the reduced purities to understand decoherence dynamics in many-electron systems is exemplified through an analysis of the electronic decoherence dynamics in a model molecular system.

Fragment-based time-dependent density functional theory

Using the Runge-Gross theorem that establishes the foundation of time-dependent density functional theory, we prove that for a given electronic Hamiltonian, choice of initial state, and choice of fragmentation, there is a unique single-particle potential (dubbed time-dependent partition potential) which, when added to each of the preselected fragment potentials, forces the fragment densities to evolve in such a way that their sum equals the exact molecular density at all times. This uniqueness theorem suggests new ways of computing the time-dependent properties of electronic systems via fragment-time-dependent density functional theory calculations. We derive a formally exact relationship between the partition potential and the total density, and illustrate our approach on a simple model system for binary fragmentation in a laser field.

Time-dependent density functional theory for many-electron systems interacting with cavity photons

Time-dependent (current) density functional theory for many-electron systems strongly coupled to quantized electromagnetic modes of a microcavity is proposed. It is shown that the electron-photon wave function is a unique functional of the electronic (current) density and the expectation values of photonic coordinates. The Kohn-Sham system is constructed, which allows us to calculate the above basic variables by solving self-consistent equations for noninteracting particles. We suggest possible approximations for the exchange-correlation potentials and discuss implications of this approach for the theory of open quantum systems. In particular we show that it naturally leads to time-dependent density functional theory for systems coupled to the Caldeira-Leggett bath.

Quantum transport through an array of quantum dots

The transient current through an array of as many as 1000 quantum dots is simulated with two newly developed quantum mechanical methods. To our surprise, upon switching on the bias voltage, the current increases linearly with time before reaching its steady state value. And the time required for the current to reach its steady state value is proportional to the length of the array, and more interestingly, is exactly the time for a conducting electron to travel through the array at the Fermi velocity. These quantum phenomena can be understood by a simple analysis on the energetics of an equivalent classical circuit. An experimental design is proposed to confirm the numerical findings.

A stochastic approach to open quantum systems

Stochastic methods are ubiquitous to a variety of fields, ranging from physics to economics and mathematics. In many cases, in the investigation of natural processes, stochasticity arises every time one considers the dynamics of a system in contact with a somewhat bigger system, an environment with which it is considered in thermal equilibrium. Any small fluctuation of the environment has some random effect on the system. In physics, stochastic methods have been applied to the investigation of phase transitions, thermal and electrical noise, thermal relaxation, quantum information, Brownian motion and so on. In this review, we will focus on the so-called stochastic Schrödinger equation. This is useful as a starting point to investigate the dynamics of open quantum systems capable of exchanging energy and momentum with an external environment. We discuss in some detail the general derivation of a stochastic Schrödinger equation and some of its recent applications to spin thermal transport, thermal relaxation, and Bose-Einstein condensation. We thoroughly discuss the advantages of this formalism with respect to the more common approach in terms of the reduced density matrix. The applications discussed here constitute only a few examples of a much wider range of applicability.

Effect of cusps in time-dependent quantum mechanics

Within nonrelativistic quantum mechanics, spatial cusps in initial wave functions can lead to nonanalytic behavior in time. We suggest a method for calculating the short-time behavior in such situations. For these cases, the density does not match its Taylor expansion in time, but the Runge-Gross proof of time-dependent density functional theory still holds, as it requires only the potential to be time analytic.

Simplified Sum-Over-States Approach for Predicting Resonance Raman Spectra. Application to Nucleic Acid Bases

Resonance Raman spectra provide a valuable probe into molecular excited-state structures and properties. Moreover, resonance enhancement is of importance for the chemical contribution to surface-enhanced Raman scattering. In this work, we introduce a simplified sum-over-states scheme for computing Raman spectra and Raman excitation profiles. The proposed sum-over-states approach uses derivatives of electronic excitation energies and transition dipole moments, which can be efficiently computed from time-dependent density functional theory. We analyze and interpret the resonance Raman spectra and Raman excitation profiles of nucleic acid bases using the present approach. Contributions of individual excited states under strictly resonant and nonresonant conditions are investigated, and smooth interpolation between both limiting cases is obtained.