On the time evolution of fermionic occupation numbers

del Campo J. Excited States by Coupling Piris Natural Orbital Functionals with the Extended Random-Phase Approximation

In this work, we explore the use of Piris natural orbital functionals (PNOFs) to calculate excited-state energies by coupling their reconstructed second-order reduced density matrix with the extended random-phase approximation (ERPA). We have named the general method PNOF-ERPA, and specific approaches are referred to as PNOF-ERPA0, PNOF-ERPA1, and PNOF-ERPA2, according to the way the excitation operator is built. The implementation has been tested in the first excited states of H2, HeH+, LiH, Li2, and N2 showing good results compared to the configuration interaction (CI) method. As expected, an increase in accuracy is observed on going from ERPA0 to ERPA1 and ERPA2. We also studied the effect of electron correlation included by PNOF5, PNOF7, and the recently proposed global NOF (GNOF) on the predicted excited states. PNOF5 appears to be good and may even provide better results in very small systems, but including more electron correlation becomes important as the system size increases, where GNOF achieves better results. Overall, the extension of PNOF to excited states has been successful, making it a promising method for further applications.

Outstanding improvement in removing the delocalization error by global natural orbital functional

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 correlations 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 improve the results obtained previously with PNOF7.

Electron Correlation in the Iron(II) Porphyrin by Natural Orbital Functional Approximations

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 second-order 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.

Time dependent reduced density matrix functional theory at strong correlation: insights from a two-site Anderson impurity model

The one-body density matrix has recently attracted considerable attention as a promising key quantity for the description of systems out of equilibrium. Its time evolution is given in terms of the two-body density matrix, and thus the central challenge is to find approximations to the latter. An extra layer of difficulty is added when dealing with strong electron correlations. In this work, we explore precisely this regime by looking at the two-site Anderson impurity model as a case study. To address the system's dynamics, we use an adiabatic approximation based on the exact ground-state two-body density matrix. We find that this adiabatic extension does not reproduce the exact results even for a slow switch-on of the external perturbation, and we trace back this behavior to the lack of an accurate imaginary part of the adiabatic approximation to the two-body density matrix. The attempt to restore an approximate imaginary part through a Hilbert transform of the real part works well only for very short times, but quickly deteriorates for longer times, with the one-body density matrix being pushed out of its N-representability domain. Our results thus pose an important constraint on practical prescriptions to perform the time evolution of the one-body density matrix.