Adiabatic approximation of time-dependent density matrix functional response theory

Time-dependent density matrix functional theory can be formulated in terms of coupled-perturbed response equations, in which a coupling matrix K(omega) features, analogous to the well-known time-dependent density functional theory (TDDFT) case. An adiabatic approximation is needed to solve these equations, but the adiabatic approximation is much more critical since there is not a good "zero order" as in TDDFT, in which the virtual-occupied Kohn-Sham orbital energy differences serve this purpose. We discuss a simple approximation proposed earlier which uses only results from static calculations, called the static approximation (SA), and show that it is deficient, since it leads to zero response of the natural orbital occupation numbers. This leads to wrong behavior in the omega-->0 limit. An improved adiabatic approximation (AA) is formulated. The two-electron system affords a derivation of exact coupled-perturbed equations for the density matrix response, permitting analytical comparison of the adiabatic approximation with the exact equations. For the two-electron system also, the exact density matrix functional (2-matrix in terms of 1-matrix) is known, enabling testing of the static and adiabatic approximations unobscured by approximations in the functional. The two-electron HeH(+) molecule shows that at the equilibrium distance, SA consistently underestimates the frequency-dependent polarizability alpha(omega), the adiabatic TDDFT overestimates alpha(omega), while AA improves upon SA and, indeed, AA produces the correct alpha(0). For stretched HeH(+), adiabatic density matrix functional theory corrects the too low first excitation energy and overpolarization of adiabatic TDDFT methods and exhibits excellent agreement with high-quality CCSD ("exact") results over a large omega range.     

The structure of the second-order reduced density matrix in density matrix functional theory and its construction from formal criteria

Some formal requirements for the second-order reduced density matrix are discussed in the context of density matrix functional theory. They serve as a basis for the ad hoc construction of the second-order reduced density matrix in terms of the first-order reduced density matrix and lead to implicit functionals where the occupation numbers of the natural orbitals are obtained as diagonal elements of an idempotent matrix the elements of which represent the variational parameters to be optimized. The numerical results obtained from a first realization of such an implicit density matrix functional give excellent agreement with the results of full configuration interaction calculations for four-electron systems like LiH and Be. Results for H2O taken as an example for a somewhat larger molecule are numerically less satisfactory but still give reasonable occupation numbers of the natural orbitals and indicate the capability of density matrix functional theory to cope with static electron correlation.     

Analytic second-order energy derivatives in natural orbital functional theory

The analytic energy gradients in the atomic orbital representation have recently been published (Mitxelena and Piris in J Chem Phys 146:014102, 2017) within the framework of the natural orbital functional theory (NOFT). We provide here an alternative expression for them in terms of natural orbitals, and use it to derive the analytic second-order energy derivatives with respect to nuclear displacements in the NOFT. The computational burden is shifted to the calculation of perturbed natural orbitals and occupancies, since a set of linear coupled-perturbed equations obtained from the variational Euler equations must be solved to attain the analytic Hessian at the perturbed geometry. The linear response of both natural orbitals and occupation numbers to nuclear geometry displacements need only specify the reconstruction of the second-order reduced density matrix in terms of occupation numbers.     

Time-dependent density-functional and reduced density-matrix methods for few electrons: Exact versus adiabatic approximations

To address the impact of electron correlations in the linear and non-linear response regimes of interacting many-electron systems exposed to time-dependent external fields, we study one-dimensional (1D) systems where the interacting problem is solved exactly by exploiting the mapping of the 1D N-electron problem onto an N-dimensional single electron problem. We analyze the performance of the recently derived 1D local density approximation as well as the exact-exchange orbital functional for those systems. We show that the interaction with an external resonant laser field shows Rabi oscillations which are detuned due to the lack of memory in adiabatic approximations. To investigate situations where static correlations play a role, we consider the time-evolution of the natural occupation numbers associated to the reduced one-body density matrix. Those studies shed light on the non-locality and time-dependence of the exchange and correlation functionals in time-dependent density and density-matrix functional theories.     

Higher-order response in O(N) by perturbed projection

Perturbed projection for linear scaling solution of the coupled-perturbed self-consistent-field equations [V. Weber, A.M.N. Niklasson, and M. Challacombe, Phys. Rev. Lett. 92, 193002 (2004)] is extended to the computation of higher-order static response properties. Although generally applicable, perturbed projection is further developed here in the context of the self-consistent first and second electric hyperpolarizabilities at the Hartree-Fock level of theory. Nonorthogonal, density-matrix analogs of Wigner's 2n+1 rule valid for linear one-electron perturbations are given up to fourth order. Linear scaling and locality of the higher-order response densities under perturbation by a global electric field are demonstrated for three-dimensional water clusters.     

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

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 Lo?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.     

New constraints upon the electron-electron repulsion energy functional of the one-electron reduced density matrix

Three strict constraints upon the electron-electron repulsion energy functional of the one-electron reduced density matrix (the 1-matrix) are obtained by combining its invariance and stationary properties with the extended Koopmans' theorem. The constraints relate the quantities derived from the functional pertaining to an N-electron system with those of its (N-1)-electron congener. Together with the N-representability requirement for the 1-matrix of the congener, identities involving the electron-electron repulsion energies of the two systems and their derivatives with respect to the 1-matrices seriously narrow down the choices for potential approximate density-matrix functionals. This fact is well illustrated in the case of two-electron systems, where the validity of the new constraints is confirmed and found to originate from a nontrivial cancellation among different terms. Thus, the constraints provide a new tool for the construction and testing of new functionals that complements the previously known conditions such as the reproduction of the homogeneous gas energies and momentum distributions, convexity, and the N-representability of the associated 2-matrices.     

Excitation energies from range-separated time-dependent density and density matrix functional theory

Time-dependent density functional theory (TD-DFT) in the adiabatic formulation exhibits known failures when applied to predicting excitation energies. One of them is the lack of the doubly excited configurations. On the other hand, the time-dependent theory based on a one-electron reduced density matrix functional (time-dependent density matrix functional theory, TD-DMFT) has proven accurate in determining single and double excitations of H(2) molecule if the exact functional is employed in the adiabatic approximation. We propose a new approach for computing excited state energies that relies on functionals of electron density and one-electron reduced density matrix, where the latter is applied in the long-range region of electron-electron interactions. A similar approach has been recently successfully employed in predicting ground state potential energy curves of diatomic molecules even in the dissociation limit, where static correlation effects are dominating. In the paper, a time-dependent functional theory based on the range-separation of electronic interaction operator is rigorously formulated. To turn the approach into a practical scheme the adiabatic approximation is proposed for the short- and long-range components of the coupling matrix present in the linear response equations. In the end, the problem of finding excitation energies is turned into an eigenproblem for a symmetric matrix. Assignment of obtained excitations is discussed and it is shown how to identify double excitations from the analysis of approximate transition density matrix elements. The proposed method used with the short-range local density approximation (srLDA) and the long-range Buijse-Baerends density matrix functional (lrBB) is applied to H(2) molecule (at equilibrium geometry and in the dissociation limit) and to Be atom. The method accounts for double excitations in the investigated systems but, unfortunately, the accuracy of some of them is poor. The quality of the other excitations is in general much better than that offered by TD-DFT-LDA or TD-DMFT-BB approximations if the range-separation parameter is properly chosen. The latter remains an open problem.     

Projected gradient algorithms for Hartree-Fock and density matrix functional theory calculations

We present projected gradient algorithms designed for optimizing various functionals defined on the set of N-representable one-electron reduced density matrices. We show that projected gradient algorithms are efficient in minimizing the Hartree-Fock or the Muller-Buijse-Baerends functional. On the other hand, they converge very slowly when applied to the recently proposed BBk (k=1,2,3) functionals [O. Gritsenko et al., J. Chem. Phys. 122, 204102 (2005)]. This is due to the fact that the BBk functionals are not proper functionals of the density matrix.     

Direct minimization of the energy in density functional theory

The minimization of the energy functional of the first-order density matrix γ( r , r ') is achieved using unitary transformations applied to γ. Equivalently, such transformations can be carried out also on one-electron orbitals (natural orbitals) and their occupation (integer or non-integer) numbers. The conventional local density approximation based on the electron density p( r ) is then considered as a special case. The direct minimization of the energy functional of p with respect to the parameters of the unitary transformation leads to stationary conditions that are all equivalent to the Kohn–Sham equations. Preliminary numerical tests show that the proposed algorithms for the direct minimization of the energy work in a satisfactory manner. © John Wiley & Sons, Inc.     

Performance of the density matrix functional theory in the quantum theory of atoms in molecules

The generalization to arbitrary molecular geometries of the energetic partitioning provided by the atomic virial theorem of the quantum theory of atoms in molecules (QTAIM) leads to an exact and chemically intuitive energy partitioning scheme, the interacting quantum atoms (IQA) approach, that depends on the availability of second-order reduced density matrices (2-RDMs). This work explores the performance of this approach in particular and of the QTAIM in general with approximate 2-RDMs obtained from the density matrix functional theory (DMFT), which rests on the natural expansion (natural orbitals and their corresponding occupation numbers) of the first-order reduced density matrix (1-RDM). A number of these functionals have been implemented in the promolden code and used to perform QTAIM and IQA analyses on several representative molecules and model chemical reactions. Total energies, covalent intra- and interbasin exchange-correlation interactions, as well as localization and delocalization indices have been determined with these functionals from 1-RDMs obtained at different levels of theory. Results are compared to the values computed from the exact 2-RDMs, whenever possible.     

Inequalities for fermion density matrices

A partial trace over the occupation numbers of all but k states in the density matrix of an ensemble with an arbitrary number of single-particle states is defined as the (reduced) k-state density matrix. This matrix is used to obtain a complete, practical solution to the problem of determining the representability of the diagonal elements of the one- and two-particle (reduced) density matrices. This solution is expressed as a series of linear inequalities involving the density-matrix elements; the inequalities are identical with those derived previously by Davidson and McCrae by a different method. In addition, our method is used to obtain nonlinear, matrix inequalities on the off-diagonal elements of the density matrices.     

A one-electron approximation to domain-averaged Fermi hole analysis

In general, full domain-averaged Fermi hole (DAFH) analysis for correlated wavefunctions requires explicit use of the correlated pair density, but such a quantity is not always readily available. We propose instead a simple one-electron approximation, which we call pseudo-DAFH or pDAFH, and which requires instead only the natural orbitals (and their occupation numbers). From comparisons of the DAFH and pDAFH modes of analysis for the bond dissociation processes in H2, N2 and LiH, as well as for the electronic structure of more complex bonding patterns, such as in CH2Li2 and Li4, we conclude that pDAFH analysis could indeed prove to be very useful when the correlated pair density is not available. Detailed comparisons are also presented of values of the shared-electron distribution index (SEDI), a proposed one-electron approximation to it (pSEDI) and a generalized Wiberg index.     

Generalized density-functional theory: Conquering the N -representability problem with exact functionals for the electron pair density and the second-order reduced density matrix

Using the constrained search and Legendre-transform formalisms, one can derive “generalized” density-functional theories, in which the fundamental variable is either the electron pair density or the second-order reduced density matrix. In both approaches, theN-representability problem is solved by the functional, and the variational principle is with respect to all pair densities (density matrices) that are nonnegative and appropriately normalized. The Legendre-transform formulation provides a lower bound on the constrained-search functional. Noting that experience in density-functional and density-matrix theories suggests that it is easier to approximate functionals than it is to approximate the set ofN-representable densities sheds some light on the significance of this work.     

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, Mu?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.     

Wigner’s (2<Emphasis Type="Italic">n</Emphasis> + 1) rule for nonlinear Schrödinger equations

Wigner’s (2n + 1) rule is proven for general total energy functionals, constructed from the expectation value of a linear Hamiltonian and a nonlinear functional of the electron density or of the first order density matrix. Such functionals are common in independent particle models, like the Kohn–Sham density functional or Hartree–Fock theories, but they occur as well in reaction-field type solvent effect models and range-separated hybrid density functional approaches. The fulfillment of the (2n + 1) rule is crucial for the development of efficient perturbation approaches for the treatment of one-electron and two-electron perturbations. A general transformation formula is derived, that removes some of the (2n + 1)-rule violating matrix elements of the perturbation operator from the general expression of the arbitrary order perturbational energy correction of the nonlinear problem.     

Benchmark calculations for reduced density-matrix functional theory

Reduced density-matrix functional theory (RDMFT) is a promising alternative approach to the problem of electron correlation. Like standard density functional theory, it contains an unknown exchange-correlation functional, for which several approximations have been proposed in the last years. In this article, we benchmark some of these functionals in an extended set of molecules with respect to total and atomization energies. Our results show that the most recent RDMFT functionals give very satisfactory results compared to standard quantum chemistry and density functional approaches.     

Projected quasiparticle theory for molecular electronic structure

We derive and implement symmetry-projected Hartree-Fock-Bogoliubov (HFB) equations and apply them to the molecular electronic structure problem. All symmetries (particle number, spin, spatial, and complex conjugation) are deliberately broken and restored in a self-consistent variation-after-projection approach. We show that the resulting method yields a comprehensive black-box treatment of static correlations with effective one-electron (mean-field) computational cost. The ensuing wave function is of multireference character and permeates the entire Hilbert space of the problem. The energy expression is different from regular HFB theory but remains a functional of an independent quasiparticle density matrix. All reduced density matrices are expressible as an integration of transition density matrices over a gauge grid. We present several proof-of-principle examples demonstrating the compelling power of projected quasiparticle theory for quantum chemistry.     

Correlation strength and information entropy

The correlation present in the nondegenerate ground state of an interacting Fermi system is discussed in terms of reduced density matrices and their cumulant expansion. By generalizing a result obtained for the interacting uniform electron gas (correlation induced exchange-hole narrowing), possible measures of the correlation strength in terms of natural occupation numbers (the eigenvalues of the true one-particle density matrix) are introduced. These quantities-the v-order nonidempotency and the information entropy of the natural occupation numbers-result from the correlated many-body wave function and characterize the ground-state correlation in addition to the usual correlation energy. The uniform electron gas serves as a first illustrative example. © 1995 John Wiley & Sons, Inc.