Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

The Journal of Chemical Physics

The adiabatic approximation in time-dependent density matrix functional theory: Response properties from dynamics of phase-including natural orbitals

Klaas J. H. Giesbertz, O. V. Gritsenko, Evert Jan Baerends

Published Nov 4, 2010• 32 citations
Publisher

Abstract

The adiabatic approximation is problematic in time-dependent density matrix functional theory. With pure density matrix functionals (invariant under phase change of the natural orbitals) it leads to lack of response in the occupation numbers, hence wrong frequency dependent responses, in particular α(ω→0)≠α(0) (the static polarizability). We propose to relinquish the requirement that the functional must be a pure one-body reduced density matrix (1RDM) functional, and to introduce additional variables which can be interpreted as phases of the one-particle states of the independent particle reference system formed with the natural orbitals, thus obtaining so-called phase-including natural orbital (PINO) functionals. We also stress the importance of the correct choice of the complex conjugation in the two-electron integrals in the commonly used functionals (they should not be of exchange type). We demonstrate with the Löwdin-Shull energy expression for two-electron systems, which is an example of a PINO functional, that for two-electron systems exact responses (polarizabilities, excitation energies) are obtained, while writing this energy expression in the usual way as a 1RDM functional yields erroneous responses.

DOI
10.1063/1.3499601

Related papers

Response calculations based on an independent particle system with the exact one-particle density matrix: Polarizabilities

Klaas J. H. Giesbertz, O. V. Gritsenko, Evert Jan Baerends

Recently, we have demonstrated that the problems finding a suitable adiabatic approximation in time-dependent one-body reduced density matrix functional theory can be remedied by introducing an additional degree of freedom to describe the system: the phase of the natural orbitals [K. J. H. Giesbertz, O. V. Gritsenko, and E. J. Baerends, Phys. Rev. Lett. 105, 013002 (2010); K. J. H. Giesbertz, O. V. Gritsenko, and E. J. Baerends, J. Chem. Phys. 133, 174119 (2010)]. In this article we will show in detail how the frequency-dependent response equations give the proper static limit (ω → 0), including the perturbation in the chemical potential, which is required in static response theory to ensure the correct number of particles. Additionally we show results for the polarizability for H2 and compare the performance of two different two-electron functionals: the phase-including Löwdin-Shull functional and the density matrix form of the Löwdin-Shull functional.

arXiv:1311.3544•The Journal of Chemical Physics•Mar 7, 2014•8 citations

Response Calculations with an Independent Particle System with an Exact One-Particle Density Matrix

Klaas J. H. Giesbertz, O. V. Gritsenko, Evert Jan Baerends

We use the natural orbitals to define an independent particle system, from which the exact one-particle density matrix can be obtained with an ensemble of degenerate determinantal ground states. Also defining explicit phases for the orbitals, and admitting functionals that are dependent on those phases, time-dependent equations for the orbitals and occupation numbers are obtained from an action principle. The wrong polarizability and lack of double excitations of straightforward adiabatic time-dependent density matrix functional theory are then corrected, and the important symmetry χ(ω)=χ{*}(-ω), lost in previous ad hoc improvements, is restored. The extension of the response calculations beyond the occupied-virtual pairs, which are the only ones admitted in time-dependent density functional theory, leads to greatly improved response properties.

Physical Review Letters•Jul 2, 2010•44 citations

Comment on "Response calculations with an independent particle system with an exact one-particle density matrix''

Ryan Requist, Oleg Pankratov

We point out an error in the argument [PRL 105, 013002 (2010)] that the time independence of the occupation numbers in the adiabatic approximation follows from the invariance of the ground-state interaction energy functional with respect to changes in the phases of the natural orbitals.

arXiv:1011.1482•Nov 5, 2010•0 citations

Adiabatic approximation in time-dependent reduced-density-matrix functional theory

Ryan Requist, Oleg Pankratov

With the aim of describing real-time electron dynamics, we introduce an adiabatic approximation for the equation of motion of the one-body reduced density matrix (one-matrix). The eigenvalues of the one-matrix, which represent the occupation numbers of single-particle orbitals, are obtained from the constrained minimization of the instantaneous ground-state energy functional rather than from their dynamical equations. The performance of the approximation vis-\`a-vis nonadiabatic effects is assessed in real-time simulations of a two-site Hubbard model. Due to Landau-Zener-type transitions, the system evolves into a nonstationary state with persistent oscillations in the observables. The amplitude of the oscillations displays a strongly nonmonotonic dependence on the strength of the electron-electron interaction and the rate of variation of the external potential. We interpret an associated resonance behavior in the phase of the oscillations in terms of ``scattering'' with spectator energy levels. To clarify the motivation for the minimization condition, we derive a sequence of energy functionals ${E}_{v}^{(n)}$, for which the corresponding sequence of minimizing one-matrices is asymptotic to the exact one-matrix in the adiabatic limit.

arXiv:0911.0945•Physical Review A•Apr 26, 2010•46 citations

Time-Dependent Quantum Many-Body Systems: Linear Response, Electronic Transport, and Reduced Density Matrices

Heiko Appel

Time-dependent density functional theory (TDDFT) provides a successful approach to calculate excitation energies of atomic and molecular systems. In part I of this work we present a double-pole approximation (DPA) to the response equations of TDDFT. The double-pole approximation provides an exact description of systems with two strongly coupled excitations which are isolated from the rest of the spectrum. In contrast to the traditional single- pole approximation of TDDFT the DPA also yields corrections to the Kohn-Sham oscillator strengths. Several critical pole separations can be identified, e.g. we find that the pole coupling can cause transitions to vanish entirely from the optical spectrum. We also demonstrate how to invert the double-pole solution which allows us to predict matrix elements of the exchange- correlation kernel fxc from experimental input. This can serve as benchmark for the construction of future approximations for the kernel fxc. Reduced density matrix functional theory (RDMFT) has emerged recently as promising candidate to treat strongly correlated electronic many-body systems beyond traditional density functional theory (DFT). The research within RDMFT was so far focussed on the static theory. In this work we attempt some first steps towards a time-dependent generalization of RDMFT. In part II we derive equations of motion for natural orbitals and occupation numbers. Using the equation of motion for the occupation numbers we show that an adiabatic extension of presently known ground-state functionals of static RDMFT always leads to occupation numbers which are constant in time. From the stationary conditions of the equations of motion for the N-body correlations (correlated parts of the N-body matrices) we derive a new class of ground-state functionals which can be used in static RDMFT. Applications are presented for a one-dimensional model system where the time-dependent many-body Schrödinger equation can be propagated numerically. We use optimal control theory to find optimized laser pulses for transitions in a model for atomic Helium. From the numerically exact correlated wavefunction we extract the exact time evolution of natural orbitals and occupation numbers for (i) laser-driven Helium and (ii) electron-ion scattering. Part III of this work considers time-dependent quantum transport within TDDFT. We present an algorithm for the calculation of extended eigenstates of single- particle Hamiltonians which is especially tailored to a finite-difference discretization of the Schrödinger equation. We consider the propagation of finite mesoscopic systems and demonstrate the limitations of such an approach. To overcome the shortcomings of a description of quantum transport in terms of finite systems we develop a time-propagation scheme for extended states which utilizes a mixed basis representation. Our discretization scheme allows to treat central device and lead regions on the same footing thus preventing artificial reflections at grid boundaries.

•Jan 1, 2007•10 citations

Orbital Functionals in Density-Matrix- and Current-Density-Functional Theory

N. Helbig

Density-Functional Theory (DFT), although widely used and very successful in the calculation of several observables, fails to correctly describe strongly correlated materials. In the first part of this work we, therefore, introduce reduced-density-matrix-functional theory (RDMFT) which is one possible way to treat electron correlation beyond DFT. Within this theory the one-body reduced density matrix (1-RDM) is used as the basic variable. Our main interest is the calculation of the fundamental gap which proves very problematic within DFT. In order to calculate the fundamental gap we generalize RDMFT to fractional particle numbers M by describing the system as an ensemble of an N and an N+1 particle system (with N≤ M≤N+1). For each fixed particle number, M, the total energy is minimized with respect to the natural orbitals and their occupation numbers. This leads to the total energy as a function of M. The derivative of this function with respect to the particle number has a discontinuity at integer particle number which is identical to the gap. In addition, we investigate the necessary and sufficient conditions for the 1-RDM of a system with fractional particle number to be N-representable. Numerical results are presented for alkali atoms, small molecules, and periodic systems. Another problem within DFT is the description of non-relativistic many- electron systems in the presence of magnetic fields. It requires the paramagnetic current density and the spin magnetization to be used as basic variables besides the electron density. However, electron-gas-based functionals of current-spin-density-functional Theory (CSDFT) exhibit derivative discontinuities as a function of the magnetic field whenever a new Landau level is occupied, which makes them difficult to use in practice. Since the appearance of Landau levels is, intrinsically, an orbital effect it is appealing to use orbital-dependent functionals. We have developed a CSDFT version of the optimized effective potential (OEP) method which allows for the use of explicitly orbital-dependent functionals in the context of CSDFT. We present the derivation of corresponding equations and show results for a quantum dot in external magnetic fields in the second part of this thesis.

•Jan 1, 2006•3 citations
The adiabatic approximation is problematic in time-dependent density matrix functional theory. With pure density matrix functionals (invariant under phase change of the natural orbitals) it leads to lack of response in the occupation numbers, hence wrong frequency dependent responses, in particular α(ω→0)≠α(0) (the static polarizability). We propose to relinquish the requirement that the functional must be a pure one-body reduced density matrix (1RDM) functional, and to introduce additional variables which can be interpreted as phases of the one-particle states of the independent particle reference system formed with the natural orbitals, thus obtaining so-called phase-including natural orbital (PINO) functionals. We also stress the importance of the correct choice of the complex conjugation in the two-electron integrals in the commonly used functionals (they should not be of exchange type). We demonstrate with the Löwdin-Shull energy expression for two-electron systems, which is an example of a PINO functional, that for two-electron systems exact responses (polarizabilities, excitation energies) are obtained, while writing this energy expression in the usual way as a 1RDM functional yields erroneous responses.
Recently, we have demonstrated that the problems finding a suitable adiabatic approximation in time-dependent one-body reduced density matrix functional theory can be remedied by introducing an additional degree of freedom to describe the system: the phase of the natural orbitals [K. J. H. Giesbertz, O. V. Gritsenko, and E. J. Baerends, Phys. Rev. Lett. 105, 013002 (2010); K. J. H. Giesbertz, O. V. Gritsenko, and E. J. Baerends, J. Chem. Phys. 133, 174119 (2010)]. In this article we will show in detail how the frequency-dependent response equations give the proper static limit (ω → 0), including the perturbation in the chemical potential, which is required in static response theory to ensure the correct number of particles. Additionally we show results for the polarizability for H2 and compare the performance of two different two-electron functionals: the phase-including Löwdin-Shull functional and the density matrix form of the Löwdin-Shull functional.
We use the natural orbitals to define an independent particle system, from which the exact one-particle density matrix can be obtained with an ensemble of degenerate determinantal ground states. Also defining explicit phases for the orbitals, and admitting functionals that are dependent on those phases, time-dependent equations for the orbitals and occupation numbers are obtained from an action principle. The wrong polarizability and lack of double excitations of straightforward adiabatic time-dependent density matrix functional theory are then corrected, and the important symmetry χ(ω)=χ{*}(-ω), lost in previous ad hoc improvements, is restored. The extension of the response calculations beyond the occupied-virtual pairs, which are the only ones admitted in time-dependent density functional theory, leads to greatly improved response properties.
With the aim of describing real-time electron dynamics, we introduce an adiabatic approximation for the equation of motion of the one-body reduced density matrix (one-matrix). The eigenvalues of the one-matrix, which represent the occupation numbers of single-particle orbitals, are obtained from the constrained minimization of the instantaneous ground-state energy functional rather than from their dynamical equations. The performance of the approximation vis-\`a-vis nonadiabatic effects is assessed in real-time simulations of a two-site Hubbard model. Due to Landau-Zener-type transitions, the system evolves into a nonstationary state with persistent oscillations in the observables. The amplitude of the oscillations displays a strongly nonmonotonic dependence on the strength of the electron-electron interaction and the rate of variation of the external potential. We interpret an associated resonance behavior in the phase of the oscillations in terms of ``scattering'' with spectator energy levels. To clarify the motivation for the minimization condition, we derive a sequence of energy functionals Ev(n){E}_{v}^{(n)}Ev(n)​, for which the corresponding sequence of minimizing one-matrices is asymptotic to the exact one-matrix in the adiabatic limit.
Time-dependent density functional theory (TDDFT) provides a successful approach to calculate excitation energies of atomic and molecular systems. In part I of this work we present a double-pole approximation (DPA) to the response equations of TDDFT. The double-pole approximation provides an exact description of systems with two strongly coupled excitations which are isolated from the rest of the spectrum. In contrast to the traditional single- pole approximation of TDDFT the DPA also yields corrections to the Kohn-Sham oscillator strengths. Several critical pole separations can be identified, e.g. we find that the pole coupling can cause transitions to vanish entirely from the optical spectrum. We also demonstrate how to invert the double-pole solution which allows us to predict matrix elements of the exchange- correlation kernel fxc from experimental input. This can serve as benchmark for the construction of future approximations for the kernel fxc. Reduced density matrix functional theory (RDMFT) has emerged recently as promising candidate to treat strongly correlated electronic many-body systems beyond traditional density functional theory (DFT). The research within RDMFT was so far focussed on the static theory. In this work we attempt some first steps towards a time-dependent generalization of RDMFT. In part II we derive equations of motion for natural orbitals and occupation numbers. Using the equation of motion for the occupation numbers we show that an adiabatic extension of presently known ground-state functionals of static RDMFT always leads to occupation numbers which are constant in time. From the stationary conditions of the equations of motion for the N-body correlations (correlated parts of the N-body matrices) we derive a new class of ground-state functionals which can be used in static RDMFT. Applications are presented for a one-dimensional model system where the time-dependent many-body Schrödinger equation can be propagated numerically. We use optimal control theory to find optimized laser pulses for transitions in a model for atomic Helium. From the numerically exact correlated wavefunction we extract the exact time evolution of natural orbitals and occupation numbers for (i) laser-driven Helium and (ii) electron-ion scattering. Part III of this work considers time-dependent quantum transport within TDDFT. We present an algorithm for the calculation of extended eigenstates of single- particle Hamiltonians which is especially tailored to a finite-difference discretization of the Schrödinger equation. We consider the propagation of finite mesoscopic systems and demonstrate the limitations of such an approach. To overcome the shortcomings of a description of quantum transport in terms of finite systems we develop a time-propagation scheme for extended states which utilizes a mixed basis representation. Our discretization scheme allows to treat central device and lead regions on the same footing thus preventing artificial reflections at grid boundaries.
Density-Functional Theory (DFT), although widely used and very successful in the calculation of several observables, fails to correctly describe strongly correlated materials. In the first part of this work we, therefore, introduce reduced-density-matrix-functional theory (RDMFT) which is one possible way to treat electron correlation beyond DFT. Within this theory the one-body reduced density matrix (1-RDM) is used as the basic variable. Our main interest is the calculation of the fundamental gap which proves very problematic within DFT. In order to calculate the fundamental gap we generalize RDMFT to fractional particle numbers M by describing the system as an ensemble of an N and an N+1 particle system (with N≤ M≤N+1). For each fixed particle number, M, the total energy is minimized with respect to the natural orbitals and their occupation numbers. This leads to the total energy as a function of M. The derivative of this function with respect to the particle number has a discontinuity at integer particle number which is identical to the gap. In addition, we investigate the necessary and sufficient conditions for the 1-RDM of a system with fractional particle number to be N-representable. Numerical results are presented for alkali atoms, small molecules, and periodic systems. Another problem within DFT is the description of non-relativistic many- electron systems in the presence of magnetic fields. It requires the paramagnetic current density and the spin magnetization to be used as basic variables besides the electron density. However, electron-gas-based functionals of current-spin-density-functional Theory (CSDFT) exhibit derivative discontinuities as a function of the magnetic field whenever a new Landau level is occupied, which makes them difficult to use in practice. Since the appearance of Landau levels is, intrinsically, an orbital effect it is appealing to use orbital-dependent functionals. We have developed a CSDFT version of the optimized effective potential (OEP) method which allows for the use of explicitly orbital-dependent functionals in the context of CSDFT. We present the derivation of corresponding equations and show results for a quantum dot in external magnetic fields in the second part of this thesis.