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.

The relative stability of the singlet, triplet, and quintet spin states of iron(II) porphyrin (FeP) represents a challenging problem for electronic structure methods. While it is currently accepted that the ground state is a triplet, multiconfigurational wave function-based methods predict a quintet, and density functional approximations vary between triplet and quintet states, leading to a prediction that highly depends on the features of the method employed. The recently proposed Global Natural Orbital Functional (GNOF) aims to provide a balanced treatment between static and dynamic correlation, and together with the previous Piris Natural Orbital Functionals (PNOFs), allowed us to explore the importance of each type of correlation in the stability order of the states of FeP with a method that conserves the spin of the system. It is noteworthy that GNOF correlates all electrons in all available orbitals for a given basis set; in the case of the FeP with a double-ζ basis set as used in this work, this means that GNOF can properly correlate 186 electrons in 465 orbitals, significantly increasing the sizes of systems amenable to multiconfigurational treatment. Results show that PNOF5, PNOF7s, and PNOF7 predict the quintet to have a lower energy than the triplet state; however, the addition of dynamic correlation via secondorder Møller−Plesset corrections (NOF-MP2) turns the triplet state to be lower than the quintet state, a prediction also reproduced by GNOF that incorporates much more dynamic correlation than its predecessors.

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.

Computation of molecular orbital electron repulsion integrals (MO-ERIs) as a transformation from atomic orbital ERIs (AO-ERIs) is the bottleneck of second-order electron propagator calculations when a single orbital is studied. In this contribution, asymmetric density fitting is combined with modified Cholesky decomposition to generate efficiently the required MO-ERIs. The key point of the presented algorithms is to keep track of integrals through partial contractions performed on three-center AO-ERIs; these contractions are stored in RAM instead of the AO-ERIs. Two implementations are provided, an in-core, which reduces the arithmetic and memory scaling factors as compared to the four-center AO-ERIs contraction method, and a semidirect, which overcomes memory limitations by evaluating antisymmetrized MO-ERIs in batches. On the numerical side, the proposed approach is fast and stable. The bad effects due to ill conditioning, namely, several negative and close to zero eigenvalues due to machine round off errors of the matrix associated with the density fitting process, are effectively controlled by means of a modified Cholesky factorization that avoids the matrix inversion needed to perform the asymmetrical density fitting implementation. The numerical experience presented shows that the in-core implementation is highly competitive to perform calculations on medium and large basis sets, while the semidirect implementation has small variations in time by changes in the available memory. The general applicability is illustrated on a set of selected relatively large-size molecules.

This work assesses the performance of the recently proposed global natural orbital functional (GNOF) against the charge delocalization error. GNOF provides a good balance between static and dynamic electronic correlation leading to accurate total energies while preserving spin, even for systems with a highly multi-configurational character. Several analyses were applied to the functional, namely i) how the charge is distributed in super-systems of two fragments, ii) the stability of ionization potentials while increasing the system size, and iii) potential energy curves of a neutral and charged diatomic system. GNOF was found to practically eliminate the charge delocalization error in many of the studied systems or greatly improves the results obtained previously with PNOF7.