0000000000076313
AUTHOR
Klaas J. H. Giesbertz
Long-range interactions and the sign of natural amplitudes in two-electron systems
In singlet two-electron systems the natural occupation numbers of the one-particle reduced density matrix are given as squares of the natural amplitudes which are defined as the expansion coefficients of the two-electron wave function in a natural orbital basis. In this work we relate the sign of the natural amplitudes to the nature of the two-body interaction. We show that long-range Coulomb-type interactions are responsible for the appearance of positive amplitudes and give both analytical and numerical examples that illustrate how the long-distance structure of the wave function affects these amplitudes. We further demonstrate that the amplitudes show an avoided crossing behavior as func…
Density-potential mappings in quantum dynamics
In a recent letter [Europhys. Lett. 95, 13001 (2011)] the question of whether the density of a time-dependent quantum system determines its external potential was reformulated as a fixed point problem. This idea was used to generalize the existence and uniqueness theorems underlying time-dependent density functional theory. In this work we extend this proof to allow for more general norms and provide a numerical implementation of the fixed-point iteration scheme. We focus on the one-dimensional case as it allows for a more in-depth analysis using singular Sturm-Liouville theory and at the same time provides an easy visualization of the numerical applications in space and time. We give an ex…
Oscillator Strengths of Electronic Excitations with Response Theory using Phase Including Natural Orbital Functionals
The key characteristics of electronic excitations of many-electron systems, the excitation energies ωα and the oscillator strengths fα, can be obtained from linear response theory. In one-electron models and within the adiabatic approximation, the zeros of the inverse response matrix, which occur at the excitation energies, can be obtained from a simple diagonalization. Particular cases are the eigenvalue equations of time-dependent density functional theory (TDDFT), time-dependent density matrix functional theory, and the recently developed phase-including natural orbital (PINO) functional theory. In this paper, an expression for the oscillator strengths fα of the electronic excitations is…
Approximate energy functionals for one-body reduced density matrix functional theory from many-body perturbation theory
We develop a systematic approach to construct energy functionals of the one-particle reduced density matrix (1RDM) for equilibrium systems at finite temperature. The starting point of our formulation is the grand potential $\Omega [\mathbf{G}]$ regarded as variational functional of the Green's function $G$ based on diagrammatic many-body perturbation theory and for which we consider either the Klein or Luttinger-Ward form. By restricting the input Green's function to be one-to-one related to a set on one-particle reduced density matrices (1RDM) this functional becomes a functional of the 1RDM. To establish the one-to-one mapping we use that, at any finite temperature and for a given 1RDM $\…
Time-Dependent Reduced Density Matrix Functional Theory
In this chapter we will give an introduction into one-body reduced density matrix functional theory (RDMFT). This is a rather new method to deal with the quantum many-body problem. Especially the development of a time-dependent version, TDRDMFT , is very recent. Therefore, there are many open questions and the formalism has not crystalized yet into a standard form such as in (TD)DFT. Although RDMFT has similarities with DFT, there are many more differences. This chapter is too short for a full introduction into the wondrous world of RDMFT, but we hope to give an idea what (TD)RDMFT might bring.
Natural occupation numbers: When do they vanish?
The non-vanishing of the natural orbital occupation numbers of the one-particle density matrix of many-body systems has important consequences for the existence of a density matrix-potential mapping for nonlocal potentials in reduced density matrix functional theory and for the validity of the extended Koopmans' Theorem. On the basis of Weyl's theorem we give a connection between the differentiability properties of the ground state wave function and the rate at which the natural occupations approach zero when ordered as a descending series. We show, in particular, that the presence of a Coulomb cusp in the wave function leads, in general, to a power law decay of the natural occupations, whe…
Towards nonlocal density functionals by explicit modelling of the exchange-correlation hole in inhomogeneous systems
We put forward new approach for the development of a non-local density functional by a direct modeling of the shape of exchange-correlation (xc) hole in inhomogeneous systems. The functional is aimed at giving an accurate xc-energy and an accurate corresponding xc-potential even in difficult near-degeneracy situations such as molecular bond breaking. In particular we demand that: (1) the xc hole properly contains -1 electron, (2) the xc-potential has the asymptotic -1/r behavior outside finite systems and (3) the xc-potential has the correct step structure related to the derivative discontinuities of the xc-energy functional. None of the currently existing functionals satisfies all these re…
Compact two-electron wave function for bond dissociation and Van der Waals interactions: A natural amplitude assessment
Electron correlations in molecules can be divided in short range dynamical correlations, long range Van der Waals type interactions and near degeneracy static correlations. In this work we analyze for a one-dimensional model of a two-electron system how these three types of correlations can be incorporated in a simple wave function of restricted functional form consisting of an orbital product multiplied by a single correlation function $f(r_{12})$ depending on the interelectronic distance $r_{12}$. Since the three types of correlations mentioned lead to different signatures in terms of the natural orbital (NO) amplitudes in two-electron systems we make an analysis of the wave function in t…
Response calculations based on an independent particle system with the exact one-particle density matrix: Excitation energies
Adiabatic response time-dependent density functional theory (TDDFT) suffers from the restriction to basically an occupied → virtual single excitation formulation. Adiabatic time-dependent density matrix functional theory allows to break away from this restriction. Problematic excitations for TDDFT, viz. bonding-antibonding, double, charge transfer, and higher excitations, are calculated along the bond-dissociation coordinate of the prototype molecules H2 and HeH+ using the recently developed adiabatic linear response phase-including (PI) natural orbital theory (PINO). The possibility to systematically increase the scope of the calculation from excitations out of (strongly) occupied into wea…