Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

The Journal of Chemical Physics

An adiabatic connection for doubly-occupied configuration interaction wave functions

Nam Vu, Ion Mitxelena, A. Eugene DePrince

Published Dec 28, 2019• 11 citations
Publisher

Abstract

An adiabatic connection (AC) is developed as an electron correlation correction for doubly occupied configuration interaction (DOCI) wave functions. Following the work of Pernal [Phys. Rev. Lett. 120, 013001 (2018)], the working equations of the approach, termed AC-DOCI, are rooted in the extended random phase approximation (ERPA) and require knowledge of only the ground-state two-electron reduced density matrix (2RDM) from the DOCI. As such, the AC is naturally suited to modeling electron correlation in variational 2RDM (v2RDM)-based approximations to the DOCI. The v2RDM-driven AC-DOCI is applied to the dissociation of molecular nitrogen and the double dissociation of water; the approach yields energies that are similar in quality to those from second-order multireference perturbation theory near equilibrium, but the quality of the AC-DOCI energy degrades at stretched geometries. The exact adiabatic connection path suggests the assumption that the one-electron reduced-density matrix is constant along the AC path is invalid at stretched geometries, but this deficiency alone cannot explain the observed behavior. Rather, it appears that the ERPA's single-particle-transition ansatz cannot, in general, provide good approximations to the 2RDM along the AC path. The AC-DOCI is also applied to a set of 45 reaction energies; for these systems, the approach has an average accuracy that is comparable to that of single-reference second-order many-body perturbation theory.

DOI
10.1063/1.5130660

Related papers

Exact and approximate adiabatic connection formulae for the correlation energy in multireference ground and excited states

Katarzyna Pernal

Recently it has been shown how to employ the adiabatic connection (AC) formalism to obtain correlation energy for multireference wavefunctions [K. Pernal, Phys. Rev. Lett. 120, 013001 (2018)]. Approximations to the exact AC formulation have been based on assuming that a one-electron reduced density matrix is constant along the AC path and by employing the extended random phase approximation. In this paper, the importance of these approximations is examined by comparing approximate AC integrands with their exact counterparts obtained for the hydrogen molecule in its ground and excited states. Encouraging results obtained for H2 indicate that AC is a viable and promising approach to a correlation energy problem not only for ground but also for excited states of electronic systems.

The Journal of Chemical Physics•Nov 26, 2018•33 citations

Toward more accurate adiabatic connection approach for multireference wavefunctions

Mikuláš Matoušek, Michał Hapka, Libor Veis +1 more

A multiconfigurational adiabatic connection (AC) formalism is an attractive approach to compute the dynamic correlation within the complete active space self-consistent field and density matrix renormalization group (DMRG) models. Practical realizations of AC have been based on two approximations: (i) fixing one- and two-electron reduced density matrices (1- and 2-RDMs) at the zero-coupling constant limit and (ii) extended random phase approximation (ERPA). This work investigates the effect of removing the "fixed-RDM" approximation in AC. The analysis is carried out for two electronic Hamiltonian partitionings: the group product function- and the Dyall Hamiltonians. Exact reference AC integrands are generated from the DMRG full configuration interaction solver. Two AC models are investigated, employing either exact 1- and 2-RDMs or their second-order expansions in the coupling constant in the ERPA equations. Calculations for model molecules indicate that lifting the fixed-RDM approximation is a viable way toward improving the accuracy of existing AC approximations.

arXiv:2210.09087•The Journal of Chemical Physics•Jan 10, 2023•13 citations

Size-extensive seniority-zero energy functionals derived from configuration interaction with double excitations

Nam Vu, A. Eugene DePrince

The doubly occupied configuration interaction (DOCI) approach can provide an accurate black-box description of nondynamic electron correlation at a computational cost that increases combinatorially with the system size. Remarkably, a pair coupled-cluster doubles (pCCD) approach (also known as the antisymmetrized product of one-reference orbital geminals) can reproduce DOCI energies with only a quadratic number of wave function parameters, and, when neglecting the cost associated with any two-electron integral transformations, these parameters can be determined at a cubic computational cost. Other simpler seniority-zero approaches derived from size-extensive modified configuration interaction doubles functionals can also provide approximations to DOCI energies at similar computational costs. We develop seniority-zero formulations of the coupled-electron pair approximation, the averaged coupled-pair functional, averaged quadratic coupled-cluster, and the parametric two-electron reduced density matrix (p2RDM) approach. These methods are Hermitian and thus offer several potential advantages over pCCD theory, including a reduction in the number of variable parameters and simplified definitions of reduced density matrices. Of the methods investigated, only the pair p2RDM (pp2RDM) approach yields energies that are comparable in quality to pCCD and DOCI. For the molecular systems investigated, pp2RDM-derived RDMs are found to be better approximations to DOCI ones than those obtained from pCCD.

The Journal of Chemical Physics•Jun 22, 2020•14 citations

Electron Correlation from the Adiabatic Connection for Multireference Wave Functions

Katarzyna Pernal

An adiabatic connection (AC) formula for the electron correlation energy is derived for a broad class of multireference wave functions. The AC expression recovers dynamic correlation energy and assures a balanced treatment of the correlation energy. Coupling the AC formalism with the extended random phase approximation allows one to find the correlation energy only from reference one- and two-electron reduced density matrices. If the generalized valence bond perfect pairing model is employed a simple closed-form expression for the approximate AC formula is obtained. This results in the overall M^{5} scaling of the computation cost making the method one of the most efficient multireference approaches accounting for dynamic electron correlation also for the strongly correlated systems.

Physical Review Letters•Jan 5, 2018•85 citations

Multi-reference perturbation theories based on the DOCI wavefunction

Yutaro Otani, Kaho Nakatani, Naoki Nakatani

We developed an alternative approach to access the full configuration interaction (FCI) energy based on perturbation expansions using the seniority-number basis. In this study, we started from the ground state of the seniority-zero Hamiltonian, also known as the doubly occupied configuration interaction (DOCI) wavefunction, to derive the Rayleigh-Schrödinger perturbation theories up to the fourth-order. We also derived the second-order perturbation theory by taking the DOCI-Fock operator as the zeroth-order Hamiltonian. From demonstrative calculations for small molecules, the DOCI-SCF-RSPT4 and the DOCI-SCF-MP2 quantitatively reproduced the FCI energy. Although solving the DOCI wavefunction is an NP-hard problem in classical computing, connections to the density matrix renormalization group or quantum computing will provide a promising way to accurately solve many-electron problems.

The Journal of Chemical Physics•Dec 4, 2025•0 citations

Polynomial scaling approximations and dynamic correlation corrections to doubly occupied configuration interaction wave functions

Mario Van Raemdonck, Diego R. Alcoba, Ward Poelmans +6 more

A class of polynomial scaling methods that approximate Doubly Occupied Configuration Interaction (DOCI) wave functions and improve the description of dynamic correlation is introduced. The accuracy of the resulting wave functions is analysed by comparing energies and studying the overlap between the newly developed methods and full configuration interaction wave functions, showing that a low energy does not necessarily entail a good approximation of the exact wave function. Due to the dependence of DOCI wave functions on the single-particle basis chosen, several orbital optimisation algorithms are introduced. An energy-based algorithm using the simulated annealing method is used as a benchmark. As a computationally more affordable alternative, a seniority number minimising algorithm is developed and compared to the energy based one revealing that the seniority minimising orbital set performs well. Given a well-chosen orbital basis, it is shown that the newly developed DOCI based wave functions are especially suitable for the computationally efficient description of static correlation and to lesser extent dynamic correlation.

The Journal of Chemical Physics•Sep 11, 2015•43 citations
An adiabatic connection (AC) is developed as an electron correlation correction for doubly occupied configuration interaction (DOCI) wave functions. Following the work of Pernal [Phys. Rev. Lett. 120, 013001 (2018)], the working equations of the approach, termed AC-DOCI, are rooted in the extended random phase approximation (ERPA) and require knowledge of only the ground-state two-electron reduced density matrix (2RDM) from the DOCI. As such, the AC is naturally suited to modeling electron correlation in variational 2RDM (v2RDM)-based approximations to the DOCI. The v2RDM-driven AC-DOCI is applied to the dissociation of molecular nitrogen and the double dissociation of water; the approach yields energies that are similar in quality to those from second-order multireference perturbation theory near equilibrium, but the quality of the AC-DOCI energy degrades at stretched geometries. The exact adiabatic connection path suggests the assumption that the one-electron reduced-density matrix is constant along the AC path is invalid at stretched geometries, but this deficiency alone cannot explain the observed behavior. Rather, it appears that the ERPA's single-particle-transition ansatz cannot, in general, provide good approximations to the 2RDM along the AC path. The AC-DOCI is also applied to a set of 45 reaction energies; for these systems, the approach has an average accuracy that is comparable to that of single-reference second-order many-body perturbation theory.
Recently it has been shown how to employ the adiabatic connection (AC) formalism to obtain correlation energy for multireference wavefunctions [K. Pernal, Phys. Rev. Lett. 120, 013001 (2018)]. Approximations to the exact AC formulation have been based on assuming that a one-electron reduced density matrix is constant along the AC path and by employing the extended random phase approximation. In this paper, the importance of these approximations is examined by comparing approximate AC integrands with their exact counterparts obtained for the hydrogen molecule in its ground and excited states. Encouraging results obtained for H2 indicate that AC is a viable and promising approach to a correlation energy problem not only for ground but also for excited states of electronic systems.
A multiconfigurational adiabatic connection (AC) formalism is an attractive approach to compute the dynamic correlation within the complete active space self-consistent field and density matrix renormalization group (DMRG) models. Practical realizations of AC have been based on two approximations: (i) fixing one- and two-electron reduced density matrices (1- and 2-RDMs) at the zero-coupling constant limit and (ii) extended random phase approximation (ERPA). This work investigates the effect of removing the "fixed-RDM" approximation in AC. The analysis is carried out for two electronic Hamiltonian partitionings: the group product function- and the Dyall Hamiltonians. Exact reference AC integrands are generated from the DMRG full configuration interaction solver. Two AC models are investigated, employing either exact 1- and 2-RDMs or their second-order expansions in the coupling constant in the ERPA equations. Calculations for model molecules indicate that lifting the fixed-RDM approximation is a viable way toward improving the accuracy of existing AC approximations.
The doubly occupied configuration interaction (DOCI) approach can provide an accurate black-box description of nondynamic electron correlation at a computational cost that increases combinatorially with the system size. Remarkably, a pair coupled-cluster doubles (pCCD) approach (also known as the antisymmetrized product of one-reference orbital geminals) can reproduce DOCI energies with only a quadratic number of wave function parameters, and, when neglecting the cost associated with any two-electron integral transformations, these parameters can be determined at a cubic computational cost. Other simpler seniority-zero approaches derived from size-extensive modified configuration interaction doubles functionals can also provide approximations to DOCI energies at similar computational costs. We develop seniority-zero formulations of the coupled-electron pair approximation, the averaged coupled-pair functional, averaged quadratic coupled-cluster, and the parametric two-electron reduced density matrix (p2RDM) approach. These methods are Hermitian and thus offer several potential advantages over pCCD theory, including a reduction in the number of variable parameters and simplified definitions of reduced density matrices. Of the methods investigated, only the pair p2RDM (pp2RDM) approach yields energies that are comparable in quality to pCCD and DOCI. For the molecular systems investigated, pp2RDM-derived RDMs are found to be better approximations to DOCI ones than those obtained from pCCD.
An adiabatic connection (AC) formula for the electron correlation energy is derived for a broad class of multireference wave functions. The AC expression recovers dynamic correlation energy and assures a balanced treatment of the correlation energy. Coupling the AC formalism with the extended random phase approximation allows one to find the correlation energy only from reference one- and two-electron reduced density matrices. If the generalized valence bond perfect pairing model is employed a simple closed-form expression for the approximate AC formula is obtained. This results in the overall M^{5} scaling of the computation cost making the method one of the most efficient multireference approaches accounting for dynamic electron correlation also for the strongly correlated systems.
We developed an alternative approach to access the full configuration interaction (FCI) energy based on perturbation expansions using the seniority-number basis. In this study, we started from the ground state of the seniority-zero Hamiltonian, also known as the doubly occupied configuration interaction (DOCI) wavefunction, to derive the Rayleigh-Schrödinger perturbation theories up to the fourth-order. We also derived the second-order perturbation theory by taking the DOCI-Fock operator as the zeroth-order Hamiltonian. From demonstrative calculations for small molecules, the DOCI-SCF-RSPT4 and the DOCI-SCF-MP2 quantitatively reproduced the FCI energy. Although solving the DOCI wavefunction is an NP-hard problem in classical computing, connections to the density matrix renormalization group or quantum computing will provide a promising way to accurately solve many-electron problems.
A class of polynomial scaling methods that approximate Doubly Occupied Configuration Interaction (DOCI) wave functions and improve the description of dynamic correlation is introduced. The accuracy of the resulting wave functions is analysed by comparing energies and studying the overlap between the newly developed methods and full configuration interaction wave functions, showing that a low energy does not necessarily entail a good approximation of the exact wave function. Due to the dependence of DOCI wave functions on the single-particle basis chosen, several orbital optimisation algorithms are introduced. An energy-based algorithm using the simulated annealing method is used as a benchmark. As a computationally more affordable alternative, a seniority number minimising algorithm is developed and compared to the energy based one revealing that the seniority minimising orbital set performs well. Given a well-chosen orbital basis, it is shown that the newly developed DOCI based wave functions are especially suitable for the computationally efficient description of static correlation and to lesser extent dynamic correlation.