Reduced Density Matrix Functional Theory Research Articles

Invertibility of the retarded response functions for initial mixed states: application to one-body reduced density matrix functional theory

ABSTRACTRecently, I have derived conditions to characterise the kernel of the retarded response function, under the assumption that the initial state is a ground state. In this paper, I demonstrate its generalisation to mixed states (ensembles). To make the proof work, the weights in the ensemble need to be decreasing for increasing energies of the pure states from which the mixed state is constructed. The resulting conditions are not easy to verify, but under the additional assumptions that the ensemble weights are directly related to the energies and that the full spectrum of the Hamiltonian participates in the ensemble, it is shown that potentials only belong to the kernel of the retarded response function if they commute with the initial Hamiltonian. These additional assumptions are valid for thermodynamic ensembles, which makes this result physically relevant. The conditions on the potentials for the thermodynamic ensembles are much stronger than in the pure state (zero temperature) case, leading to a much less involved kernel when the conditions are applied to the retarded one-body reduced density matrix response function.

Using the GVB Ansatz to develop ensemble DFT method for describing multiple strongly correlated electron pairs.

Ensemble density functional theory (DFT) furnishes a rigorous theoretical framework for describing the non-dynamic electron correlation arising from (near) degeneracy of several electronic configurations. Ensemble DFT naturally leads to fractional occupation numbers (FONs) for several Kohn-Sham (KS) orbitals, which thereby become variational parameters of the methodology. The currently available implementation of ensemble DFT in the form of the spin-restricted ensemble-referenced KS (REKS) method was originally designed for systems with only two fractionally occupied KS orbitals, which was sufficient to accurately describe dissociation of a single chemical bond or the singlet ground state of biradicaloid species. To extend applicability of the method to systems with several dissociating bonds or to polyradical species, more fractionally occupied orbitals must be included in the ensemble description. Here we investigate a possibility of developing the extended REKS methodology with the help of the generalized valence bond (GVB) wavefunction theory. The use of GVB enables one to derive a simple and physically transparent energy expression depending explicitly on the FONs of several KS orbitals. In this way, a version of the REKS method with four electrons in four fractionally occupied orbitals is derived and its accuracy in the calculation of various types of strongly correlated molecules is investigated. We propose a possible scheme to ameliorate the partial size-inconsistency that results from perfect spin-pairing. We conjecture that perfect pairing natural orbital (NO) functionals of reduced density matrix functional theory (RDMFT) should also display partial size-inconsistency.

Spectrum for nonmagnetic mott insulators from power functional within reduced density matrix functional theory.

We demonstrate that reduced density matrix functional theory (RDMFT), in conjunction with the power functional, can successfully treat the nonmagnetic insulating state of the transition metal oxides NiO and MnO, finding for both a gapped single particle spectrum. While long-range spin order is thus not necessary for qualitative agreement with experiment, we find that it is required for good agreement with the X-ray photoemission spectroscopy and Bremsstrahlung isochromat spectroscopy data. We further examine the nature of the natural orbitals in the materials, finding that they display significant Hubbard localization and are, as a consequence, very far from the corresponding Kohn-Sham orbitals. This contrasts with the case of the band insulator Si, in which the Kohn-Sham orbitals are found to be very close to the RDMFT natural orbitals.

Doping induced metal-insulator phase transition in NiO—a reduced density matrix functional theory perspective

The insulator to metal phase transition in NiO is studied within the framework of reduced density matrix functional theory (RDMFT) and density functional theory (DFT). We find that the spectral density obtained using RDMFT is in good agreement with experiments both for undoped as well as doped NiO. We find that the physical description of the hole-doping induced phase transition qualitatively differs depending on whether NiO is calculated within DFT or reduced density matrix functional. In the former case the underlying mechanism of the phase transition is identified to be a rigid shift of chemical potential, while in the latter case a redistribution of the spectral weight drives the transition. These latter results are found to be in good agreement with both experiments and previous many-body calculations.

Orbitals from local RDMFT: Are they Kohn-Sham or natural orbitals?

Recently, an approximate theoretical framework was introduced, called local reduced density matrix functional theory (local-RDMFT), where functionals of the one-body reduced density matrix (1-RDM) are minimized under the additional condition that the optimal orbitals satisfy a single electron Schrödinger equation with a local potential. In the present work, we focus on the character of these optimal orbitals. In particular, we compare orbitals obtained by local-RDMFT with those obtained with the full minimization (without the extra condition) by contrasting them against the exact NOs and orbitals from a density functional calculation using the local density approximation (LDA). We find that the orbitals from local-RMDFT are very close to LDA orbitals, contrary to those of the full minimization that resemble the exact NOs. Since local RDMFT preserves the good quality of the description of strong static correlation, this finding opens the way to a mixed density/density matrix scheme, where Kohn-Sham orbitals obtain fractional occupations from a minimization of the occupation numbers using 1-RDM functionals. This will allow for a description of strong correlation at a cost only minimally higher than a density functional calculation.

Invertibility of retarded response functions for Laplace transformable potentials: Application to one-body reduced density matrix functional theory.

A theorem for the invertibility of arbitrary response functions is presented under the following conditions: the time dependence of the potentials should be Laplace transformable and the initial state should be a ground state, though it might be degenerate. This theorem provides a rigorous foundation for all density-functional-like theories in the time-dependent linear response regime. Especially for time-dependent one-body reduced density matrix (1RDM) functional theory, this is an important step forward, since a solid foundation has currently been lacking. The theorem is equally valid for static response functions in the non-degenerate case, so can be used to characterize the uniqueness of the potential in the ground state version of the corresponding density-functional-like theory. Such a classification of the uniqueness of the non-local potential in ground state 1RDM functional theory has been lacking for decades. With the aid of presented invertibility theorem presented here, a complete classification of the non-uniqueness of the non-local potential in 1RDM functional theory can be given for the first time.

Quasipinning and selection rules for excitations in atoms and molecules

Postulated by Pauli to explain the electronic structure of atoms and molecules, the exclusion principle establishes an upper bound of 1 for fermionic natural occupation numbers ${{n}_{i}}$. A recent analysis of the pure $N$-representability problem provides a wide set of inequalities for the ${{n}_{i}}$, leading to constraints on these numbers. This has a strong potential impact on reduced density matrix functional theory as we know it. In this work we continue our study of the nature of these inequalities for some atomic and molecular systems. The results indicate that (quasi)saturation of some of them leads to selection rules for the dominant configurations in configuration interaction expansions, in favorable cases providing means for significantly reducing their computational requirements.

Reduced density-matrix functional theory: Correlation and spectroscopy.

In this work, we explore the performance of approximations to electron correlation in reduced density-matrix functional theory (RDMFT) and of approximations to the observables calculated within this theory. Our analysis focuses on the calculation of total energies, occupation numbers, removal/addition energies, and spectral functions. We use the exactly solvable Hubbard dimer at 1/4 and 1/2 fillings as test systems. This allows us to analyze the underlying physics and to elucidate the origin of the observed trends. For comparison, we also report the results of the GW approximation, where the self-energy functional is approximated, but no further hypothesis is made concerning the approximations of the observables. In particular, we focus on the atomic limit, where the two sites of the dimer are pulled apart and electrons localize on either site with equal probability, unless a small perturbation is present: this is the regime of strong electron correlation. In this limit, using the Hubbard dimer at 1/2 filling with or without a spin-symmetry-broken ground state allows us to explore how degeneracies and spin-symmetry breaking are treated in RDMFT. We find that, within the used approximations, neither in RDMFT nor in GW, the signature of strong correlation is present, when looking at the removal/addition energies and spectral function from the spin-singlet ground state, whereas both give the exact result for the spin-symmetry broken case. Moreover, we show how the spectroscopic properties change from one spin structure to the other.

Generalized Pauli constraints in reduced density matrix functional theory.

Functionals of the one-body reduced density matrix (1-RDM) are routinely minimized under Coleman's ensemble N-representability conditions. Recently, the topic of pure-state N-representability conditions, also known as generalized Pauli constraints, received increased attention following the discovery of a systematic way to derive them for any number of electrons and any finite dimensionality of the Hilbert space. The target of this work is to assess the potential impact of the enforcement of the pure-state conditions on the results of reduced density-matrix functional theory calculations. In particular, we examine whether the standard minimization of typical 1-RDM functionals under the ensemble N-representability conditions violates the pure-state conditions for prototype 3-electron systems. We also enforce the pure-state conditions, in addition to the ensemble ones, for the same systems and functionals and compare the correlation energies and optimal occupation numbers with those obtained by the enforcement of the ensemble conditions alone.

Reduced Density Matrix Functional Theory (RDMFT) and Linear Response Time-Dependent RDMFT (TD-RDMFT).

Recent advances in reduced density matrix functional theory (RDMFT) and linear response time-dependent reduced density matrix functional theory (TD-RDMFT) are reviewed. In particular, we present various approaches to develop approximate density matrix functionals which have been employed in RDMFT. We discuss the properties and performance of most available density matrix functionals. Progress in the development of functionals has been paralleled by formulation of novel RDMFT-based methods for predicting properties of molecular systems and solids. We give an overview of these methods. The time-dependent extension, TD-RDMFT, is a relatively new theory still awaiting practical and generally useful functionals which would work within the adiabatic approximation. In this chapter we concentrate on the formulation of TD-RDMFT response equations and various adiabatic approximations. None of the adiabatic approximations is fully satisfactory, so we also discuss a phase-dependent extension to TD-RDMFT employing the concept of phase-including-natural-spinorbitals (PINOs). We focus on applications of the linear response formulations to two-electron systems, for which the (almost) exact functional is known.

Quasi-particle energy spectra in local reduced density matrix functional theory.

Recently, we introduced [N. N. Lathiotakis, N. Helbig, A. Rubio, and N. I. Gidopoulos, Phys. Rev. A 90, 032511 (2014)] local reduced density matrix functional theory (local RDMFT), a theoretical scheme capable of incorporating static correlation effects in Kohn-Sham equations. Here, we apply local RDMFT to molecular systems of relatively large size, as a demonstration of its computational efficiency and its accuracy in predicting single-electron properties from the eigenvalue spectrum of the single-particle Hamiltonian with a local effective potential. We present encouraging results on the photoelectron spectrum of molecular systems and the relative stability of C20 isotopes. In addition, we propose a modelling of the fractional occupancies as functions of the orbital energies that further improves the efficiency of the method useful in applications to large systems and solids.

Local reduced-density-matrix-functional theory: Incorporating static correlation effects in Kohn-Sham equations

We propose a novel scheme to bring reduced density matrix functional theory (RDMFT) into the realm of density functional theory (DFT) that preserves the accurate density functional description at equilibrium, while incorporating accurately static and left-right correlation effects in molecules and keeping the good computational performance of DFT-based schemes. The key ingredient is to relax the requirement that the local potential is the functional derivative of the energy with respect to the density. Instead, we propose to restrict the search for the approximate natural orbitals within a domain where these orbitals are eigenfunctions of a single-particle hamiltonian with a local effective potential. In this way, fractional natural occupation numbers are accommodated into Kohn-Sham equations allowing for the description of molecular dissociation without breaking spin symmetry. Additionally, our scheme provides a natural way to connect an energy eigenvalue spectrum to the approximate natural orbitals and this spectrum is found to represent accurately the ionization potentials of atoms and small molecules.

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

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.

Natural occupation numbers: When do they vanish?

The non-vanishing of the natural orbital (NO) 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 wavefunction 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 wavefunction leads, in general, to a power law decay of the natural occupations, whereas infinitely differentiable wavefunctions typically have natural occupations that decay exponentially. We analyze for a number of explicit examples of two-particle systems that in case the wavefunction is non-analytic at its spatial diagonal (for instance, due to the presence of a Coulomb cusp) the natural orbital occupations are non-vanishing. We further derive a more general criterium for the non-vanishing of NO occupations for two-particle wavefunctions with a certain separability structure. On the basis of this criterium we show that for a two-particle system of harmonically confined electrons with a Coulombic interaction (the so-called Hookium) the natural orbital occupations never vanish.

Spectral Density and Metal-Insulator Phase Transition in Mott Insulators within Reduced Density Matrix Functional Theory

We present a method for calculating the spectrum of periodic solids within reduced density matrix functional theory. This method is validated by a detailed comparison of the angular momentum projected spectral density with that of well-established many-body techniques, finding very good agreement in all cases. The physics behind the pressure induced insulator-metal phase transition in MnO is investigated. The driving mechanism of this transition is identified as increased crystal field splitting with pressure, resulting in a charge redistribution between the Mn e(g) and t(2)g symmetry projected states.

Minimization procedure in reduced density matrix functional theory by means of an effective noninteracting system

In this work, we propose a self-consistent minimization procedure for functionals in reduced density matrix functional theory. We introduce an effective noninteracting system at finite temperature which is capable of reproducing the groundstate one-reduced density matrix of an interacting system at zero temperature. By introducing the concept of a temperature tensor the minimization with respect to the occupation numbers is shown to be greatly improved.

Correlation measures as benchmarks in reduced density matrix functional theory

AbstractIn reduced density matrix functional theory, approximate total energy expressions in terms of the one‐body reduced density matrix (1RDM) are minimized under the appropriate N‐representability subsidiary conditions. The optimal 1RDM is in general nonidempotent. Measures of nonidempotency are the total charge of weakly occupied states and correlation entropy given in terms of the eigenvalues of the 1RDM. In this work, we calculate these quantities for several systems using different 1RDM functionals and compare with accurate results from a multideterminant method. Our goal is to assess 1RDM functionals with respect to the quality of the obtained 1RDM. © 2012 Wiley Periodicals, Inc.

Open-shell reduced density matrix functional theory

Open-shell reduced density matrix functional theory is established by investigating the domain of the exact functional. For spin states that are the ground state, a particularly simple set is found to be the domain. It cannot be generalized to other spin states. A number of conditions satisfied by the exact density matrix functional is formulated and tested for approximate functionals. The exact functional does not suffer from fractional spin error, which is the source of the static correlation error in dissociated molecules. We prove that a simple approximation (called the Buijse-Baerends functional, Müller or square root functional) has a non-positive fractional spin error. In the case of the H atom the error is zero. Numerical results for a few atoms are given for approximate density and density matrix functionals as well as a recently developed range-separated combination of both.

Constraints of reduced density-matrix functional theory for the two-dimensional homogeneous electron gas

Reduced density-matrix functional theory (RDMFT) has become an appealing alternative to density-functional theory to describe electronic properties of highly-correlated systems. Here we derive exact conditions for the suitability of RDMFT to describe the two-dimensional homogeneous electron gas, which is the base system for, e.g., semiconductor quantum dots and quantum Hall devices. Following the method of Cioslowski and Pernal [J. Chem. Phys. {\bf 111}, 3396 (1999)] we focus on the properties of power functionals of the form $f(n,n')=(n n')^\alpha$ for the scaling function in the exchange-correlation energy. We show that in order to have stable and analytic solutions, and for $f$ to satisfy the homogeneous scaling constraint, the power is restricted to $1/4 \leq \alpha \leq 3/4$. Applying a reasonable ansatz for the momentum distribution and the lower bound for the exchange-correlation energy tightens the physical regime further to $0.64 \lesssim \alpha \leq 0.75$.

Discontinuities of the Chemical Potential in Reduced Density Matrix Functional Theory

Abstract Using the discontinuity of the chemical potential as a function of excess charge, the fundamental gaps for finite systems and the band gaps of extended solids are determined within reduced density matrix functional theory. We also present the necessary and sufficient conditions for the one-body reduced density matrix of a system with fractional charge to be ensemble N-representable. The performance of most modern day reduced density matrix functionals is assessed for the gaps and the correlation energy of finite systems. Our results show that for finite systems the PNOF, BBC3, and power functionals yield very accurate correlation energies while for a correct description of the fundamental gap the removal of self-interaction terms is essential. For extended solids we find that the power functional captures the correct band gap behavior for conventional semiconductors as well as strongly correlated Mott insulators, where a gap is obtained in absence of any magnetic ordering.