Assessment of a new approach for the two-electron cumulant in natural-orbital-functional theory

Coupled cluster-inspired geminal wavefunctions

Electron pairs have an illustrious history in chemistry, from powerful concepts to understanding structural stability and reactive changes to the promise of serving as building blocks of quantitative descriptions of the electronic structure of complex molecules and materials. However, traditionally, two-electron wavefunctions (geminals) have not enjoyed the popularity and widespread use of the more standard single-particle methods. This has changed recently, with a renewed interest in the development of geminal wavefunctions as an alternative to describing strongly correlated phenomena. Hence, there is a need to find geminal methods that are accurate, computationally tractable, and do not demand significant input from the user (particularly via cumbersome and often ill-behaved orbital optimization steps). Here, we propose new families of geminal wavefunctions inspired by the pair coupled cluster doubles ansatz. We present a new hierarchy of two-electron wavefunctions that extends the one-reference orbital idea to other geminals. Moreover, we show how to incorporate single-like excitations in this framework without leaving the quasiparticle picture. We explore the role of imposing seniority restrictions on these wavefunctions and benchmark these new methods on model strongly correlated systems.

del Campo J. Charge delocalization error in Piris Natural Orbital Functionals

Piris Natural Orbital Functionals (PNOF) have been recognized as a low-scaling alternative to study strong correlated systems. In this work, we address the performance of the fifth functional (PNOF5) and the seventh functional (PNOF7) to deal with another common problem, the charge delocalization error. The effects of this problem can be observed in charged systems of repeated well-separated fragments, where the energy should be the sum of the charged and neutral fragments, regardless of how the charge is distributed. In practice, an energetic overstabilization of fractional charged fragments leads to a preference for having the charge delocalized throughout the system. To establish the performance of PNOF functionals regarding charge delocalization error, charged chains of helium atoms and the W4-17-MR set molecules were used as base fragments and their energy, charge distribution and correlation regime were studied. It was found that PNOF5 prefers localized charge distributions, while PNOF7 improves the treatment of interpair static correlation and tends to the correct energetic limit for several cases, although a preference for delocalized charge distributions may arise in highly strong correlation regimes. Overall, it is concluded that PNOF functionals can simultaneously deal with static correlation and charge delocalization errors, resulting in a promising choice to study charge-related problems.

Photoemission Spectra from the Extended Koopman's Theorem, Revisited

The Extended Koopman's Theorem (EKT) provides a straightforward way to compute charged excitations from any level of theory. In this work we make the link with the many-body effective energy theory (MEET) that we derived to calculate the spectral function, which is directly related to photoemission spectra. In particular, we show that at its lowest level of approximation the MEET removal and addition energies correspond to the so-called diagonal approximation of the EKT. Thanks to this link, the EKT and the MEET can benefit from mutual insight. In particular, one can readily extend the EKT to calculate the full spectral function, and choose a more optimal basis set for the MEET by solving the EKT secular equation. We illustrate these findings with the examples of the Hubbard dimer and bulk silicon.

Resolution of the identity approximation applied to PNOF correlation calculations

In this work, the required algebra to employ the resolution of the identity approximation within the Piris Natural Orbital Functional (PNOF) is developed, leading to an implementation named DoNOF-RI. The arithmetic scaling is reduced from fifth-order to fourth-order, and the memory scaling is reduced from fourth-order to third-order, allowing significant computational time savings. After the DoNOF-RI calculation has fully converged, a restart with four-center electron repulsion integrals can be performed to remove the effect of the auxiliary basis set incompleteness, quickly converging to the exact result. The proposed approach has been tested on cycloalkanes and other molecules of general interest to study the numerical results, as well as the speed-ups achieved by PNOF7-RI when compared with PNOF7.

Capturing static and dynamic correlation with ΔNO-MP2 and ΔNO-CCSD

The ΔNO method for static correlation is combined with second-order Møller-Plesset perturbation theory (MP2) and coupled-cluster singles and doubles (CCSD) to account for dynamic correlation. The MP2 and CCSD expressions are adapted from finite-temperature CCSD, which includes orbital occupancies and vacancies, and expanded orbital summations. Correlation is partitioned with the aid of damping factors incorporated into the MP2 and CCSD residual equations. Potential energy curves for a selection of diatomics are in good agreement with extrapolated full configuration interaction results and on par with conventional multireference approaches.

Reliability of interacting quantum atoms (IQA) data computed from post-HF densities: impact of the approximation used

The BBC1 approximation is recommended for IQA calculations; MP2/BBC1 and CCSD/BBC1 produced highly comparable FAMSEC-based interpretations of intramolecular interactions.

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.

Discontinuities of the Chemical Potential in Reduced Density Matrix Functional Theory

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.

Charge transfer, double and bond-breaking excitations with time-dependent density matrix functional theory

Time-dependent density functional theory (TDDFT) in its current adiabatic implementations exhibits three remarkable failures: (a) completely wrong behavior of the excited state surface along a bond-breaking coordinate; (b) lack of doubly excited configurations; (c) much too low charge transfer excitation energies. These TDDFT failure cases are all strikingly exhibited by prototype two-electron systems such as dissociating H2 and HeH+. We find for these systems with time-dependent density matrix functional theory that: (a) Within previously formulated simple adiabatic approximations, the bonding-to-antibonding excited state surface as well as charge transfer excitations are described without problems, but not the double excitations; (b) An adiabatic approximation is formulated in which also the double excitations are fully accounted for.

Dispersion interactions within the Piris natural orbital functional theory: The helium dimer

The authors have investigated the description of the dispersion interaction within the Piris natural orbital functional (PNOF) theory. The PNOF arises from an explicit antisymmetric approach for the two-particle cumulant in terms of two symmetric matrices, Delta and Lambda. The functional forms of these matrices are obtained from the generalization of the two-particle system expressions, except for the off-diagonal elements of Delta. The mean value theorem and the partial sum rule obtained for the off-diagonal elements of Delta provide a prescription for deriving practical functionals. In particular, the previous employed approximation {Jpp/2} for the mean values {Jp*} affords several molecular properties but it is incapable to account for dispersion effects. In this work, the authors analyze a new approach for Jp* obtained by factorization of the matrix Delta within the bounds on its off-diagonal elements imposed by the positivity conditions of the two-particle reduced density matrix. Additional terms for the matrix elements of Lambda proportional to the square root of the holes are again introduced to describe properly the occupation numbers of the lowest occupied levels. The authors have found that the cross products between weakly occupied orbitals must be removed from the functional form of Lambda to obtain a correct long-range asymptotic behavior. The PNOF is used to predict the binding energy as well as the equilibrium distance of the helium dimer. The results are compared with the full configuration-interaction calculations and the corresponding experimental data.