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

Softmax parameterization of the occupation numbers for natural orbital functionals based on electron pairing approaches

Within the framework of natural orbital functional theory, having a convenient representation of the occupation numbers and orbitals becomes critical for the computational performance of the calculations. Recognizing this, we propose an innovative parametrization of the occupation numbers that takes advantage of the electron-pairing approach used in Piris natural orbital functionals through the adoption of the softmax function, a pivotal component in modern deep-learning models. Our approach not only ensures adherence to the N-representability of the first-order reduced density matrix (1RDM) but also significantly enhances the computational efficiency of 1RDM functional theory calculations. The effectiveness of this alternative parameterization approach was assessed using the W4-17-MR molecular set, which demonstrated faster and more robust convergence compared to previous implementations.

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.

Time evolution of natural orbitals in ab initio molecular dynamics

This work combines for the first time ab initio molecular dynamics (AIMD) within the Born-Oppenheimer approximation with a global natural orbital functional (GNOF), an approximate functional of the one-particle reduced density matrix. The most prominent feature of GNOF-AIMD is its ability to display the real-time evolution of natural orbitals, providing detailed information on the time-dependent electronic structure of complex systems and processes, including reactive collisions. The quartet ground-state reaction N(4S) + H2(1Σ) → NH(3Σ) + H(2S) is taken as a validation test. Collision energy influences on integral cross sections for different initial rovibrational states of H2 and rotational-state distributions of the NH product are discussed, showing a good agreement with previous high-quality theoretical results.

A unified density matrix functional construction of quantum baths in density matrix embedding theory beyond the mean-field approximation

The equivalence in one-electron quantum baths between the practical implementation of density matrix embedding theory (DMET) and the more recent Householder-transformed density matrix functional embedding theory has been shown previously in the standard but special case where the reference full-size (one-electron reduced) density matrix, from which the bath is constructed, is idempotent [S. Yalouz et al., J. Chem. Phys. 157, 214112 (2022)]. We prove mathematically that the equivalence remains valid when the density matrix is not idempotent anymore, thus allowing for the construction of correlated (one-electron) quantum baths. A density-matrix functional exactification of DMET is derived within the present unified quantum embedding formalism. Numerical examples reveal that the embedding cluster can be quite sensitive to the level of density-matrix functional approximation used for computing the reference density matrix.