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

arXiv:1503.04878 · The Journal of Chemical Physics

Pair extended coupled cluster doubles

Thomas M. Henderson, Ireneusz W. Bulik, Gustavo E. Scuseria

Published Jun 5, 2015• 70 citations
arXivHTMLPDF

Abstract

The accurate and efficient description of strongly correlated systems remains an important challenge for computational methods. Doubly occupied configuration interaction (DOCI), in which all electrons are paired and no correlations which break these pairs are permitted, can in many cases provide an accurate account of strong correlations, albeit at combinatorial computational cost. Recently, there has been significant interest in a method we refer to as pair coupled cluster doubles (pCCD), a variant of coupled cluster doubles in which the electrons are paired. This is simply because pCCD provides energies nearly identical to those of DOCI, but at mean-field computational cost (disregarding the cost of the two-electron integral transformation). Here, we introduce the more complete pair extended coupled cluster doubles (pECCD) approach which, like pCCD, has mean-field cost and reproduces DOCI energetically. We show that unlike pCCD, pECCD also reproduces the DOCI wave function with high accuracy. Moreover, pECCD yields sensible albeit inexact results even for attractive interactions where pCCD breaks down.

DOI
10.1063/1.4921986

Related papers

Seniority-based coupled cluster theory

Thomas M. Henderson, Ireneusz W. Bulik, Tamar Stein +1 more

Doubly occupied configuration interaction (DOCI) with optimized orbitals often accurately describes strong correlations while working in a Hilbert space much smaller than that needed for full configuration interaction. However, the scaling of such calculations remains combinatorial with system size. Pair coupled cluster doubles (pCCD) is very successful in reproducing DOCI energetically, but can do so with low polynomial scaling (N(3), disregarding the two-electron integral transformation from atomic to molecular orbitals). We show here several examples illustrating the success of pCCD in reproducing both the DOCI energy and wave function and show how this success frequently comes about. What DOCI and pCCD lack are an effective treatment of dynamic correlations, which we here add by including higher-seniority cluster amplitudes which are excluded from pCCD. This frozen pair coupled cluster approach is comparable in cost to traditional closed-shell coupled cluster methods with results that are competitive for weakly correlated systems and often superior for the description of strongly correlated systems.

arXiv:1410.6529•The Journal of Chemical Physics•Dec 22, 2014•166 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

Influence of broken-pair excitations on the exact pair wavefunction

Jacob M. Wahlen-Strothman, Thomas M. Henderson, Gustavo E. Scuseria

Doubly occupied configuration interaction (DOCI), the exact diagonalisation of the Hamiltonian in the paired (seniority zero) sector of the Hilbert space, is a combinatorial cost wave function that can be very efficiently approximated by pair coupled cluster doubles (pCCD) at mean-field computational cost. As such, it is a very interesting candidate as a starting point for building the full configuration interaction (FCI) ground state eigenfunction belonging to all (not just paired) seniority sectors. The true seniority zero sector of FCI (referred to here as FCI0) includes the effect of coupling between all seniority sectors rather than just seniority zero, and is, in principle, different from DOCI. We here study the accuracy with which DOCI approximates FCI0. Using a set of small Hubbard lattices, where FCI is possible, we show that DOCI ∼ FCI0 under weak correlation. However, in the strong correlation regime, the nature of the FCI0 wavefunction can change significantly, rendering DOCI and pCCD a less than ideal starting point for approximating FCI.

arXiv:1708.05696•Molecular Physics•Sep 8, 2017•13 citations

Using full configuration interaction quantum Monte Carlo in a seniority zero space to investigate the correlation energy equivalence of pair coupled cluster doubles and doubly occupied configuration interaction

James J. Shepherd, Thomas M. Henderson, Gustavo E. Scuseria

Over the past few years, pair coupled cluster doubles (pCCD) has shown promise for the description of strong correlation. This promise is related to its apparent ability to match results from doubly occupied configuration interaction (DOCI), even though the latter method has exponential computational cost. Here, by modifying the full configuration interaction quantum Monte Carlo algorithm to sample only the seniority zero sector of Hilbert space, we show that the DOCI and pCCD energies are in agreement for a variety of 2D Hubbard models, including for systems well out of reach for conventional configuration interaction algorithms. Our calculations are aided by the sign problem being much reduced in the seniority zero space compared with the full space. We present evidence for this and then discuss the sign problem in terms of the wave function of the system which appears to have a simplified sign structure.

arXiv:1602.03543•The Journal of Chemical Physics•Mar 7, 2016•34 citations

Advancement of the pair coupled cluster doubles-based molecular property calculations with quantum embedding

Rahul Chakraborty

The works included in this thesis aim to extend the applications of pair-coupled cluster doubles (pCCD) theory. pCCD produces reliable results for strongly correlated chemical systems and performs better than other traditional quantum chemical methods. Still, achieving greater accuracy requires post-pCCD corrections, which steeply increase computational cost, limiting its widespread application. Quantum embedding approaches, which reserve accurate but high-cost methods for the active site of chemical change and approximated methods for other regions of a chemical system, can be used to circumvent this problem. In this context, two new pCCD-based embedding formulations have been implemented and tested. Both the pCCD-in-DFT and pCCD-in-pCCD methods produce reliable results across a wide range of intermolecular interactions. In particular, the latter method shows significantly better performance, as it does not depend on approximated density functionals. Additionally, estimation of molecular properties is also pertinent for practical applications with pCCD. As the traditional response-theory-based formulation for pCCD analyzed here did not meet the required level of accuracy, an alternate approach based on the expectation-value method was implemented for pCCD and post-pCCD methods. This approach produces excellent agreement with the references at a much lower computational cost, as we bypass the need to solve the response equations. Thus, the pCCD-in-pCCD embedding method and the expectation-value-pCCD approach provide routes for large-scale simulations of challenging, complex systems.

•Jun 4, 2026•0 citations

Variational coupled cluster for ground and excited states

Antoine Marie, Fábris Kossoski, Pierre‐François Loos

In single-reference coupled-cluster (CC) methods, one has to solve a set of non-linear polynomial equations in order to determine the so-called amplitudes that are then used to compute the energy and other properties. Although it is of common practice to converge to the (lowest-energy) ground-state solution, it is also possible, thanks to tailored algorithms, to access higher-energy roots of these equations that may or may not correspond to genuine excited states. Here, we explore the structure of the energy landscape of variational CC and we compare it with its (projected) traditional version in the case where the excitation operator is restricted to paired double excitations (pCCD). By investigating two model systems (the symmetric stretching of the linear H4 molecule and the continuous deformation of the square H4 molecule into a rectangular arrangement) in the presence of weak and strong correlations, the performance of variational pCCD (VpCCD) and traditional pCCD is gauged against their configuration interaction (CI) equivalent, known as doubly occupied CI, for reference Slater determinants made of ground- or excited-state Hartree–Fock orbitals or state-specific orbitals optimized directly at the VpCCD level. The influence of spatial symmetry breaking is also investigated.

arXiv:2106.11305•The Journal of Chemical Physics•Sep 9, 2021•53 citations
The accurate and efficient description of strongly correlated systems remains an important challenge for computational methods. Doubly occupied configuration interaction (DOCI), in which all electrons are paired and no correlations which break these pairs are permitted, can in many cases provide an accurate account of strong correlations, albeit at combinatorial computational cost. Recently, there has been significant interest in a method we refer to as pair coupled cluster doubles (pCCD), a variant of coupled cluster doubles in which the electrons are paired. This is simply because pCCD provides energies nearly identical to those of DOCI, but at mean-field computational cost (disregarding the cost of the two-electron integral transformation). Here, we introduce the more complete pair extended coupled cluster doubles (pECCD) approach which, like pCCD, has mean-field cost and reproduces DOCI energetically. We show that unlike pCCD, pECCD also reproduces the DOCI wave function with high accuracy. Moreover, pECCD yields sensible albeit inexact results even for attractive interactions where pCCD breaks down.
Doubly occupied configuration interaction (DOCI) with optimized orbitals often accurately describes strong correlations while working in a Hilbert space much smaller than that needed for full configuration interaction. However, the scaling of such calculations remains combinatorial with system size. Pair coupled cluster doubles (pCCD) is very successful in reproducing DOCI energetically, but can do so with low polynomial scaling (N(3), disregarding the two-electron integral transformation from atomic to molecular orbitals). We show here several examples illustrating the success of pCCD in reproducing both the DOCI energy and wave function and show how this success frequently comes about. What DOCI and pCCD lack are an effective treatment of dynamic correlations, which we here add by including higher-seniority cluster amplitudes which are excluded from pCCD. This frozen pair coupled cluster approach is comparable in cost to traditional closed-shell coupled cluster methods with results that are competitive for weakly correlated systems and often superior for the description of strongly correlated systems.
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.
Doubly occupied configuration interaction (DOCI), the exact diagonalisation of the Hamiltonian in the paired (seniority zero) sector of the Hilbert space, is a combinatorial cost wave function that can be very efficiently approximated by pair coupled cluster doubles (pCCD) at mean-field computational cost. As such, it is a very interesting candidate as a starting point for building the full configuration interaction (FCI) ground state eigenfunction belonging to all (not just paired) seniority sectors. The true seniority zero sector of FCI (referred to here as FCI0) includes the effect of coupling between all seniority sectors rather than just seniority zero, and is, in principle, different from DOCI. We here study the accuracy with which DOCI approximates FCI0. Using a set of small Hubbard lattices, where FCI is possible, we show that DOCI ∼ FCI0 under weak correlation. However, in the strong correlation regime, the nature of the FCI0 wavefunction can change significantly, rendering DOCI and pCCD a less than ideal starting point for approximating FCI.
Over the past few years, pair coupled cluster doubles (pCCD) has shown promise for the description of strong correlation. This promise is related to its apparent ability to match results from doubly occupied configuration interaction (DOCI), even though the latter method has exponential computational cost. Here, by modifying the full configuration interaction quantum Monte Carlo algorithm to sample only the seniority zero sector of Hilbert space, we show that the DOCI and pCCD energies are in agreement for a variety of 2D Hubbard models, including for systems well out of reach for conventional configuration interaction algorithms. Our calculations are aided by the sign problem being much reduced in the seniority zero space compared with the full space. We present evidence for this and then discuss the sign problem in terms of the wave function of the system which appears to have a simplified sign structure.
The works included in this thesis aim to extend the applications of pair-coupled cluster doubles (pCCD) theory. pCCD produces reliable results for strongly correlated chemical systems and performs better than other traditional quantum chemical methods. Still, achieving greater accuracy requires post-pCCD corrections, which steeply increase computational cost, limiting its widespread application. Quantum embedding approaches, which reserve accurate but high-cost methods for the active site of chemical change and approximated methods for other regions of a chemical system, can be used to circumvent this problem. In this context, two new pCCD-based embedding formulations have been implemented and tested. Both the pCCD-in-DFT and pCCD-in-pCCD methods produce reliable results across a wide range of intermolecular interactions. In particular, the latter method shows significantly better performance, as it does not depend on approximated density functionals. Additionally, estimation of molecular properties is also pertinent for practical applications with pCCD. As the traditional response-theory-based formulation for pCCD analyzed here did not meet the required level of accuracy, an alternate approach based on the expectation-value method was implemented for pCCD and post-pCCD methods. This approach produces excellent agreement with the references at a much lower computational cost, as we bypass the need to solve the response equations. Thus, the pCCD-in-pCCD embedding method and the expectation-value-pCCD approach provide routes for large-scale simulations of challenging, complex systems.
In single-reference coupled-cluster (CC) methods, one has to solve a set of non-linear polynomial equations in order to determine the so-called amplitudes that are then used to compute the energy and other properties. Although it is of common practice to converge to the (lowest-energy) ground-state solution, it is also possible, thanks to tailored algorithms, to access higher-energy roots of these equations that may or may not correspond to genuine excited states. Here, we explore the structure of the energy landscape of variational CC and we compare it with its (projected) traditional version in the case where the excitation operator is restricted to paired double excitations (pCCD). By investigating two model systems (the symmetric stretching of the linear H4 molecule and the continuous deformation of the square H4 molecule into a rectangular arrangement) in the presence of weak and strong correlations, the performance of variational pCCD (VpCCD) and traditional pCCD is gauged against their configuration interaction (CI) equivalent, known as doubly occupied CI, for reference Slater determinants made of ground- or excited-state Hartree–Fock orbitals or state-specific orbitals optimized directly at the VpCCD level. The influence of spatial symmetry breaking is also investigated.