---
title: 'k-UpCCGSD: Unitary Pair Coupled-Cluster Ansatz'
url: https://www.emergentmind.com/topics/k-upccgsd
type: topic
---

# k-UpCCGSD: Unitary Pair Coupled-Cluster Ansatz

Searching arXiv for recent papers on k-UpCCGSD and foundational references.
k-UpCCGSD, short for \(k\)-unitary pair coupled-cluster generalized singles and doubles, is a layered unitary coupled-cluster ansatz used in variational quantum eigensolver (VQE) workflows for electronic-structure problems. It combines generalized single excitations with a restricted class of paired double excitations and repeats the corresponding unitary block \(k\) times, yielding a variational family that is systematically improvable in \(k\) while remaining substantially more compact than full generalized doubles constructions [1810.02327]. The ansatz has been studied as a near-term quantum computing wavefunction for ground and excited states, spin-state energetics, transition-metal oxides, reaction pathways, molecular property optimization, and larger active-space simulations [2504.08494].

## 1. Formal definition and operator content

In its standard form, k-UpCCGSD starts from a reference Slater determinant \( |\Phi_0\rangle \), usually Hartree–Fock or a restricted/open-shell Hartree–Fock determinant in an active space, and applies \(k\) successive anti-Hermitian cluster layers. A commonly used expression is
\[
|\Psi_{k\text{-UpCCGSD}}\rangle
=\prod_{\ell=1}^k
\exp\!\Bigl[
\bigl(\hat T_\text{pCCD}^{(\ell)}+\hat T_S^{(\ell)}\bigr)-\mathrm{H.c.}
\Bigr]\,
|\Phi_0\rangle,
\]
where \(\hat T_S\) is a generalized single-excitation operator and \(\hat T_\text{pCCD}\) is a pair-double operator [1810.02327]. In later active-space formulations the same structure is written as
\[
|\Psi_k(\boldsymbol{\theta})\rangle
=
\prod_{i=1}^k
\exp\!\bigl(\hat T^{(i)}(\boldsymbol{\theta})-\hat T^{(i)\dagger}(\boldsymbol{\theta})\bigr)
|\Phi_0\rangle,
\]
with \(\hat T^{(i)}=\hat T_1^{(i)}+\hat T_{2,\mathrm{pair}}^{(i)}\) [2504.08494].

The defining restriction is on the doubles sector. Standard UCC doubles allow all two-electron excitations, whereas UpCCGSD retains only paired moves between spatial orbitals. In the language used across the cited studies, these are described as seniority-zero double excitations, pair coupled-cluster doubles, or simultaneous excitation of an \(\alpha/\beta\) electron pair from one spatial orbital to another [2208.07977]. Generalized singles remain unrestricted over the active space.

The parameter \(k\) counts how many pair-GSD layers are stacked. As \(k\) increases, the ansatz becomes more expressive; one study states that in the limit \(k\to\infty\), and in the absence of Trotterization error, it spans the full unitary manifold generated by pair-GSD [2504.08494]. This makes \(k\) the principal systematic-improvement parameter.

## 2. Relation to UCCSD and UCCGSD

k-UpCCGSD is typically positioned between UCCSD and UCCGSD. UCCSD includes singles and all occupied-to-virtual doubles; UCCGSD generalizes both singles and doubles over all orbitals; k-UpCCGSD retains generalized singles but restricts doubles to pair excitations [1810.02327]. The resulting compression is the central reason for its use in VQE.

| Ansatz | Doubles sector | Reported asymptotic cost |
|---|---|---|
| UCCSD | all doubles | depth \(\mathcal O((N-\eta)^2\eta)\) |
| UCCGSD | generalized all-orbital doubles | depth \(\mathcal O(N^3)\), amplitudes \(\mathcal O(N^4)\) |
| k-UpCCGSD | paired doubles + generalized singles | depth \(\mathcal O(kN)\), amplitudes \(\mathcal O(kN^2)\) |

The foundational benchmark paper reports that k-UpCCGSD requires circuit depth \(\mathcal O(kN)\), compared with \(\mathcal O(N^3)\) for UCCGSD and \(\mathcal O((N-\eta)^2\eta)\) for UCCSD, where \(N\) is the number of spin orbitals and \(\eta\) is the number of electrons [1810.02327]. In the heme-related study, the same scaling is expressed in spatial-orbital notation: per layer, generalized singles contribute \(\mathcal O(M^2)\) parameters, paired doubles contribute \(\mathcal O(M^2)\), total parameters are \(\mathcal O(kM^2)\), and depth scales as \(O(kM)\) for \(M\) spatial orbitals [2504.08494].

System-specific resource counts illustrate the same trend. For the \(M=20\) spin-orbital active space used for Li\(_2\)Co\(_2\)O\(_4\)/Co\(_2\)O\(_4\), singles plus paired doubles give 135 parameters per layer, so k-UpCCGSD has 405 parameters for \(k=3\) and 675 for \(k=5\), with 5,850 and 9,750 CNOTs respectively; the same study reports 20,950 CNOTs for UCCSD and 78,150 for UCCGSD [2208.07977]. This suggests that the ansatz is designed as an accuracy–resource compromise rather than as a replacement for the full expressibility of UCCGSD.

## 3. Circuit realization and VQE workflows

In practical implementations, fermionic excitation generators are mapped to qubit operators, commonly with the Jordan–Wigner transformation. The heme-related study states that each \(\hat a_p^\dagger \hat a_q\) and each paired double becomes a Pauli string under Jordan–Wigner, and that each layer \(\exp(\hat T^{(i)}-\hat T^{(i)\dagger})\) is Trotter-approximated with a single Trotter step,
\[
\exp\!\bigl(\sum_j\theta_j(\hat g_j-\hat g_j^\dagger)\bigr)
\approx
\prod_j \exp\!\bigl[\theta_j(\hat g_j-\hat g_j^\dagger)\bigr],
\]
with each small exponential expanded analytically as a \(2\times2\) or \(4\times4\) rotation [2504.08494].

A complementary circuit description appears in the dibenzothiophene study. There, PennyLane’s k-UpCCGSD template maps each single-excitation generator to an \(e^{-i\theta P/2}\) rotation on a weight-two Pauli string and each paired-double generator to rotations on weight-four Pauli strings. Gate synthesis is described as basis rotation, a single-qubit \(R_Z(\theta)\), and reversal of the basis change; a weight-\(m\) Pauli string thus costs roughly \(2(m-1)\) CNOTs plus one \(R_Z\) [2512.04322].

The ansatz has also been embedded in more elaborate optimization loops. For spin-state energetics in a heme-related model, it was used in a state-averaged, orbital-optimized VQE targeting singlet, triplet, and quintet states simultaneously through
\[
\min_{\boldsymbol{\theta},\kappa}\;\sum_i w_i \big\langle\Phi_i\big|
\hat U^\dagger(\boldsymbol{\theta})\,\tilde H_q(\kappa)\,\hat U(\boldsymbol{\theta})
\big|\Phi_i\big\rangle,
\]
with equal weights \(w_i\), ADAM for the \(\theta\) parameters, and a final orbital update by PySCF’s classical orbital optimizer using 1- and 2-RDMs from the last VQE step [2504.08494]. Initial states were either single-reference ROHF determinants for each spin state or, for the triplet only, a small two-determinant multi-reference superposition denoted T1.

A different extension couples k-UpCCGSD to a two-phase optimization of the one-particle reduced density matrix. In that setting, the first phase minimizes the energy, while the second adds an RMSD-based 1-RDM penalty to the loss. For an active space \((4,4)\) and \(k=1\), the study reports 22 parameters and a CNOT depth per layer on the order of twice the number of parameters, approximately 50–60 gates; energy changes are modest, but electron density, dipole moments, and atomic charges improve substantially [2507.07667]. This indicates that k-UpCCGSD can serve not only as an energy ansatz but also as a source of reduced-density-matrix observables.

## 4. Canonical molecular benchmarks

The original comparative benchmarks on H\(_4\), H\(_2\)O, and N\(_2\) established the main empirical profile of k-UpCCGSD. For H\(_4\) in STO-3G with \(N=8\), 1-UpCCGSD uses 72 amplitudes and has a ground-state non-parallelity error \( \mathrm{NPE}_{gs}\approx 2.1 \) m\(E_h\), while 2-UpCCGSD uses 144 amplitudes and reaches \( \mathrm{NPE}_{gs}=0.0 \); UCCGSD also reaches \(0.0\) but with 3,136 amplitudes. For H\(_4\) in 6-31G with \(N=16\), 2-UpCCGSD gives \( \mathrm{NPE}_{gs}\approx 0.7 \), 3-UpCCGSD gives \( \approx 0.1 \), and UCCGSD gives \(0.0\). For H\(_2\)O double dissociation in STO-3G, 2-UpCCGSD gives \( \mathrm{NPE}_{gs}\approx 0.4 \) and 3-UpCCGSD gives \( \approx 0.02 \). For N\(_2\) dissociation in STO-3G, 4-UpCCGSD gives \( \mathrm{NPE}_{gs}\approx 3.8 \), 5-UpCCGSD gives \( \approx 1.6 \), and UCCGSD gives \( \approx 1.3 \), with chemical accuracy stated as \(\simeq 1.6\) m\(E_h\) [1810.02327].

The same work also examined excited states with orthogonally constrained VQE. Typical first-excited-state NPEs for k-UpCCGSD were 0.0 for H\(_4\)/STO-3G at 2-UpCCGSD and higher, approximately 0.9 m\(E_h\) for H\(_2\)O/STO-3G at 2-UpCCGSD and approximately 0.0 at 3-UpCCGSD, and approximately 2.6 m\(E_h\) and 1.1 m\(E_h\) for N\(_2\)/STO-3G at 4-UpCCGSD and 5-UpCCGSD, respectively [1810.02327].

A notable result in the excited-state setting is that a specialized multi-determinantal reference obtained from classical linear-response calculations improved excited-state energetics. In the N\(_2\) \(\pi\to\pi^\*\) example, improving the OC-VQE reference to a four-determinant state reduced the UCCGSD error at \(1.8\) Å from approximately \(7.9\) m\(E_h\) to approximately \(0.45\) m\(E_h\) [1810.02327]. Although that example concerns UCCGSD, it is directly relevant to k-UpCCGSD because the same paper frames multi-determinantal reference construction as a mechanism for improving excited-state calculations within constrained VQE.

## 5. Strong correlation, transition metals, and application domains

The heme-related spin-state study is one of the most detailed examinations of k-UpCCGSD in a strongly correlated transition-metal context. Using an in-house statevector simulator and single- and multi-reference trial wavefunctions, it computed singlet, triplet, and quintet energetics for active spaces from 5 to 10 spatial orbitals, equivalent to 10–20 qubits. For \(k=4\), the VQE spin-state energetics were found to agree with CASSCF to within 1–4 kcal/mol; for 5- and 6-orbital active spaces, all three spin energies and their gaps agreed with CASSCF to within 1 kcal/mol. For 7–9 orbitals, the single-reference T0 ansatz drifted up to approximately 5 kcal/mol in the worst case, whereas the multi-reference triplet T1 kept most errors within 2 kcal/mol and the singlet–quintet and triplet–quintet gaps within 4 kcal/mol. In the 10-orbital run, limited to \(k=3\) and 32-bit precision, the singlet error reached up to 15 kcal/mol for T1, although the quintet–triplet gap remained within chemical accuracy [2504.08494].

That study also reported multi-reference diagnostics \(Z_{s(1)}\). The singlet states had \(Z_{s(1)}\approx 0.11\)–0.29, triplets lay near or above \(0.1\) for most \((e,o)\), and quintets had \(Z_{s(1)}\lesssim 0.03\). The interpretation given there is that singlets show strong static correlation, triplets moderate multi-reference character, and quintets are essentially single-reference [2504.08494]. These results delimit a regime in which k-UpCCGSD is capable of reproducing spin-state energetics of strongly correlated systems but remains sensitive to reference quality and resource limits.

In transition-metal oxides, VQE simulations of Li\(_2\)Co\(_2\)O\(_4\) and Co\(_2\)O\(_4\) found that k-UpCCGSD with \(k=5\) produces results similar to UCCSD but at a lower cost. For the Li\(_2\)Co\(_2\)O\(_4\rightarrow\)Co\(_2\)O\(_4\) energy difference relative to CASCI, the reported values were \(+0.49\) kcal/mol for UCCSD, \(-3.81\) for UCCGSD, \(+1.69\) for k-UpCCGSD(\(k=3\)), and \(-0.78\) for k-UpCCGSD(\(k=5\)). The same work reports absolute ground-state errors versus CASCI of approximately \(+7.17\) kcal/mol at \(k=3\) and \(+3.46\) kcal/mol at \(k=5\), with near-equilibrium Co\(_2\)O\(_4\) potential-energy-curve agreement within approximately 0.5 kcal/mol relative to CCSD for both UCCSD and \(k=5\), but errors larger than 1 kcal/mol at stretched geometries beyond \(4.5\) Å [2208.07977].

For a reaction-pathway application, a study of chloride attack on chloromethane used a 4-qubit HOMO–LUMO active-space model and found that in noiseless simulations UCCSD and k-UpCCGSD both reproduced the full configuration interaction potential-energy surface within chemical accuracy. The reported RMSEs relative to FCI were \(7.8\times 10^{-4}\) kcal/mol for UCCSD and \(2.93\times 10^{-2}\), \(2.55\times 10^{-2}\), \(3.20\times 10^{-2}\), \(2.04\times 10^{-2}\), and \(2.10\times 10^{-2}\) kcal/mol for \(k=1,\dots,5\) k-UpCCGSD, respectively. Under a Qulacs arbitrary noise model, UCCSD had energy-error bounds of approximately 1.26 to 1.54 mHa, whereas k-UpCCGSD gave smaller error ranges, with maxima from 0.889 to 0.696 mHa over \(k=1\) to \(4\) and 0.721 mHa at \(k=5\) [2112.15314]. This is one of the clearest demonstrations in the literature that the reduced structure of k-UpCCGSD can translate into lower noise sensitivity.

At larger scale, the dibenzothiophene study used a 14-qubit \((8e,8o)\) active space and reported 252 total variational parameters, circuit depth 9,398 layers, and 15,375 total parameterized gates. The final ground-state energy was \(-864.69062\) Ha after 114 VQE iterations and 112.4 s wall time on the state-vector simulator, corresponding to a recovered correlation energy of \(-9.08\) Ha relative to the Hartree–Fock reference \(-855.6110\) Ha. The same paper contrasts this with ADAPT-VQE, which reached \(-855.5711\) Ha with a 41-layer circuit and concludes that the k-UpCCGSD circuit, although chemically accurate in simulation, is infeasible for hardware execution [2512.04322].

## 6. Limitations, misconceptions, and research directions

A recurring limitation is that the favorable asymptotic scaling of k-UpCCGSD does not imply universally shallow hardware circuits. The dibenzothiophene case, with a 9,398-layer circuit and 15,375 parameterized gates, is an explicit counterexample: the ansatz remained viable in simulation but was judged infeasible for hardware execution [2512.04322]. A common misconception is therefore to read the \(\mathcal O(kN)\) or \(\mathcal O(kM)\) depth statements as guarantees of near-term executability; the application studies show that prefactors and active-space size remain decisive.

Another limitation is expressibility in strongly correlated regimes. The heme study notes that a single-layer Trotter step introduces approximation error, although it is variationally suppressed, and reports some overstabilization of high-spin states by up to approximately 2.5 kcal/mol in a few cases [2504.08494]. The transition-metal-oxide study states that the ansatz misses higher-order three- and higher-body excitations and that a single-reference HF state may have insufficient overlap in strongly correlated regimes such as stretched bonds or open-shell Co\(^{4+}\) states [2208.07977]. The SN2 study likewise emphasizes that larger \(k\) may be required for very strongly correlated or larger active spaces, but that this can exceed NISQ coherence times [2112.15314].

The literature also shows that reference preparation matters. In the heme calculations, a small two-determinant multi-reference triplet initial state improved convergence relative to the single-reference T0 setup [2504.08494]. In the foundational excited-state benchmarks, a specialized multi-determinantal reference improved constrained-VQE energetics [1810.02327]. This suggests that the restriction to paired doubles does not eliminate the need for careful state preparation when the target state has substantial multi-reference character.

Several research directions are identified explicitly in the cited works. The heme-related study recommends choosing \(k\) as large as resources permit, enforcing spin symmetry by using pair-GSD excitation operators that commute with \(\hat S^2\) and \(\hat S_z\), and considers adaptive methods such as ADAPT-VQE, memory-streaming of cluster operators, perturbative MRPT(2) corrections on top of k-UpCCGSD, and near-term hardware experiments on strongly correlated transition-metal active spaces as future improvements [2504.08494]. The 1-RDM optimization study adds a different direction: when energies are already near classical references, augmenting the VQE objective with a density-matrix penalty can substantially improve molecular properties without changing the ansatz class itself [2507.07667].

Taken together, these results place k-UpCCGSD in a specific methodological niche. It is not the most expressive coupled-cluster ansatz and not always the most hardware-efficient in absolute terms, but it is systematically improvable, compact relative to UCCGSD, often more noise-robust than UCCSD, and already capable of chemically accurate results in several benchmark and application settings when \(k\), the active space, and the reference state are chosen appropriately [1810.02327].

Source: https://www.emergentmind.com/topics/k-upccgsd