---
title: 'Q-SENSE: Seniority-Based Quantum Subspace Expansion'
url: https://www.emergentmind.com/topics/quantum-seniority-based-subspace-expansion-q-sense
type: topic
---

# Q-SENSE: Seniority-Based Quantum Subspace Expansion

Searching arXiv for recent papers on Q-SENSE and related seniority-based quantum subspace methods.
Quantum SENiority-based Subspace Expansion (Q-SENSE) denotes a seniority-structured family of hybrid quantum-classical methods for electronic-structure simulation in which the quantum device is used on a compressed paired-electron representation and the resulting information is expanded classically into a larger determinant space. In the 2026 realization "Doubling the size of quantum selected configuration interaction based on seniority-zero space and its application to QC-QSCI-AFQMC," Q-SENSE is implemented as three coupled ideas: sampling in the seniority-zero sector through DOCI-QSCI, expanding the sampled configurations to seniority-breaking determinants through a Cartesian product of spin strings, and then using the expanded multi-determinant state as the trial wave function in phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) to recover dynamical correlation in the full orbital space [2602.07912]. Related preprints use the same label for closely allied seniority-constrained subspace-expansion and state-preparation strategies, but the common principle is the use of seniority symmetry to reduce circuit depth or qubit count while retaining a systematically improvable classical post-processing stage [2509.01061, 2504.19760].

## 1. Seniority as the organizing principle

Seniority is the number of unpaired particles. In molecular electronic structure, for spatial orbital \(i\) with spin occupations \(n_{i\uparrow}, n_{i\downarrow} \in \{0,1\}\), the seniority operator used in the 2026 Q-SENSE realization is
\[
\Omega = \sum_i \big(n_{i\uparrow} + n_{i\downarrow} - 2 n_{i\uparrow} n_{i\downarrow}\big).
\]
A Slater determinant has \(\Omega=0\) when each spatial orbital is either empty or doubly occupied, so that no singly occupied orbitals appear [2602.07912].

This modern electronic-structure usage is consistent with the older many-body notion of seniority as the number of particles not paired to angular momentum \(J=0\). In the single-shell formalism, seniority labels blocks of the many-body Hilbert space and can be expressed through quasi-spin algebra; that earlier literature also established the importance of exact and partial conservation of seniority in pairing-dominated Hamiltonians [1406.4094]. This historical connection is significant because Q-SENSE exploits the same structural fact: low-seniority sectors often capture the dominant static or pairing correlation in a much smaller subspace than the full configuration space.

Doubly occupied configuration interaction (DOCI) is the CI problem restricted to the \(\Omega=0\) sector,
\[
|\Psi_{\mathrm{DOCI}}\rangle = \sum_I c_I |\Phi_I\rangle,
\]
with determinants \( |\Phi_I\rangle \) spanning the seniority-zero subspace. The central attraction of DOCI is that it captures pair correlation efficiently and uses far fewer configurations than full CI [2602.07912]. Q-SENSE inherits this compactness, but it is not identical to DOCI. Its distinctive step is to treat seniority-zero as a compressed sampling manifold rather than as the final variational space.

A separate but related formalism appears in the orthogonal subspace-expansion variant of Q-SENSE, where orbital seniority is represented by commuting operators \(\hat{\Omega}_i\) or, under Jordan–Wigner, by stabilizers \(\hat{S}_i = \hat{Z}_{2i}\hat{Z}_{2i+1}\). In that setting, basis states from different seniority sectors are orthogonal by construction [2509.01061]. This suggests that “Q-SENSE” is best understood not as one immutable algorithm but as a seniority-centered design pattern spanning several hybrid methods.

## 2. Seniority-zero compression and the DOCI-QSCI construction

The 2026 implementation begins from quantum selected configuration interaction (QSCI), which constructs a classically diagonalizable effective Hamiltonian on a sampled subspace \(S\),
\[
\hat{H}_{\mathrm{eff}} = P_S \hat{H} P_S,\qquad
P_S = \sum_{i=1}^{R} |\Phi_i\rangle \langle \Phi_i|.
\]
The computational-basis bitstrings \( \{|\Phi_i\rangle\} \) are obtained by sampling a prepared quantum state in the computational basis, retaining \(R\) distinct bitstrings [2602.07912].

The seniority-zero restriction changes the qubit accounting. Conventional spin-orbital mappings use one qubit per spin orbital, so a device with \(Q\) qubits represents \(Q\) spin orbitals, or \(Q/2\) spatial orbitals. In seniority-zero, each spatial orbital is binary—empty or doubly occupied—so one qubit per spatial orbital suffices. With \(N\) spatial orbitals, the seniority-zero mapping therefore uses \(N\) qubits instead of \(2N\), and the accessible spatial-orbital space doubles at fixed device size. In the 2026 work, sampling is performed directly on \(N\) qubits representing spatial-orbital pair occupancies [2602.07912].

State preparation for this compressed sampling stage uses a spinless local unitary cluster Jastrow (LUCJ) ansatz with one repetition and Jastrow couplings restricted to adjacent orbitals or qubits. Circuit parameters are initialized from frozen-core CCSD amplitudes in the chosen active space. Sampling is carried out either on the IBM Quantum device `ibm_kobe` or on a noiseless simulator, and samples that violate the target electron number are discarded [2602.07912].

The restriction to \(\Omega=0\) is advantageous but incomplete. The 2026 paper states directly that this sector restriction can compromise quantitative accuracy. The reason is that the \(\Omega=0\) sector efficiently represents pair or static correlation, whereas much of the dynamical correlation resides in seniority-breaking configurations. Q-SENSE addresses this deficiency by expanding the sampled seniority-zero configurations into a larger spinful determinant set rather than by increasing the quantum-state complexity at the sampling stage [2602.07912].

## 3. Cartesian-product subspace expansion

The defining expansion step in the 2026 realization starts from \(R\) distinct \(N\)-bit seniority-zero strings in the pair-occupation representation,
\[
\{\Phi_i^{(0)}\}, \qquad i=1,\dots,R.
\]
Each bitstring is interpreted as both an \(\alpha\)-spin string and a \(\beta\)-spin string. Because DOCI imposes identical \(\alpha\) and \(\beta\) occupations, the pools are
\[
\{\Phi_i^{(\alpha)}\} = \{\Phi_i^{(\beta)}\} = \{\Phi_i^{(0)}\}.
\]
Q-SENSE then constructs a larger determinant set in the spinful space through the Cartesian product
\[
\{\tilde{\Phi}_k\} =
\{\Phi_i^{(\alpha)} \Phi_j^{(\beta)} \mid i,j \in \{1,\dots,R\}\},
\]
which yields up to \(R^2\) spin-resolved Slater determinants [2602.07912].

The crucial observation is that when \(i \neq j\), the \(\alpha\) and \(\beta\) strings differ, so singly occupied orbitals appear and the determinant has \(\Omega>0\). The effective Hamiltonian is therefore built on a space that contains both seniority-zero and seniority-breaking determinants, even though the quantum sampling cost was incurred only in the compact seniority-zero representation [2602.07912].

An optional heat-bath-like enlargement can then be applied classically. Starting from \(\{\tilde{\Phi}_k\}\), an external determinant \(l\) is added if
\[
\max_k |H_{lk} c_k| > \epsilon,
\]
where \(H_{lk}=\langle \Phi_l|\hat{H}|\Phi_k\rangle\), \(c_k\) are the current CI coefficients, and \(\epsilon\) is a user-set threshold implemented as `select_cutoff` in PySCF [2602.07912]. This step follows a conventional selected-CI growth and is conceptually separate from the seniority-based Cartesian expansion itself.

The computational consequence is a deliberate shift of cost from the quantum to the classical side. If seniority-zero sampling returns \(R\) bitstrings, the determinant count grows from \(O(R)\) to \(O(R^2)\) after the Cartesian product. The qubit count remains \(N\) at the sampling stage, but the expanded Hamiltonian must be built and diagonalized classically. This suggests a characteristic Q-SENSE trade-off: qubit compression and shallow sampling circuits are obtained at the price of a larger post-sampling linear-algebra problem [2602.07912].

## 4. End-to-end workflow and ph-AFQMC recovery of dynamical correlation

In the 2026 workflow, the electronic Hamiltonian in the active space is
\[
\hat{H} = \sum_{pq} h_{pq} a_p^\dagger a_q
+ \frac{1}{2}\sum_{pqrs} g_{pqrs} a_p^\dagger a_q^\dagger a_s a_r,
\]
with standard one-electron integrals \(h_{pq}\) and antisymmetrized two-electron integrals \(g_{pqrs}\) [2602.07912].

The algorithm proceeds in a fixed sequence. One first builds an \(N\)-qubit spinless LUCJ circuit, initializes parameters from frozen-core CCSD, executes \(S\) shots, discards electron-number-violating samples, and collects \(R\) unique seniority-zero bitstrings. One then forms the \(\alpha\)- and \(\beta\)-pools, generates all Cartesian-product determinants, assembles the effective Hamiltonian \(H_{\mathrm{eff}} = P_S \hat{H} P_S\) in that basis, and solves
\[
H_{\mathrm{eff}} c = E c
\]
classically for the ground-state eigenpair. Optional selected-CI enlargement repeats the importance-based growth and rediagonalization. The resulting multi-determinant QSCI state,
\[
|\Psi_{\mathrm{QSCI}}\rangle = \sum_k c_k |\tilde{\Phi}_k\rangle,
\]
is then used as the trial state \( \Psi_T \) in ph-AFQMC [2602.07912].

The AFQMC stage is not an auxiliary refinement but a central component of the Q-SENSE realization. The local energy for a walker Slater determinant \(|\Phi\rangle\) is
\[
E_L(\Phi) = \frac{\langle \Psi_T|\hat{H}|\Phi\rangle}
{\langle \Psi_T|\Phi\rangle}.
\]
In the phaseless formulation, the walker-weight update over time step \(\Delta\tau\) is
\[
w' = w \exp[-\Delta\tau \,\mathrm{Re}(E_L(\Phi)-E_T)] \times \max(0,\cos \Delta\theta),
\]
where \(E_T\) is a reference energy and \(\Delta\theta\) is the phase of the overlap ratio [2602.07912]. The expanded QSCI trial is intended to mitigate phaseless bias by supplying a trial wave function that already contains a broad, explicitly multi-determinant description of seniority-breaking configurations.

The implementation details are concrete. Qiskit v2.2.3 is used for the sampling circuits; `ipie` v0.6.2 is used for AFQMC; the default AFQMC time step is \(\Delta\tau=0.005\,E_h^{-1}\); and the standard settings are 50 steps and 3000 blocks unless otherwise noted. For hardware runs on `ibm_kobe`, no additional error mitigation is applied beyond postselection on valid electron number [2602.07912]. This is important because it shows that the reported behavior is not contingent on a large error-mitigation stack.

## 5. Benchmark systems and reported performance

The method is evaluated on three classes of systems: the H\(_6\) linear chain dissociation, N\(_2\) dissociation, and the addition of singlet O\(_2\) to a BODIPY dye. The H\(_6\) calculations use a cc-pVTZ active-space basis, `ibm_kobe` hardware and a noiseless simulator, \(10^5\) shots per geometry point, and 400 AFQMC walkers. The N\(_2\) calculations use cc-pVQZ, \(10^5\) shots per point, and 640 walkers, with comparison to reduced multireference CCSD and CCSD(T). The BODIPY–O\(_2\) calculations use an RB3LYP/6-31G(d)-optimized reaction path and IRC, energies at 6-31G(d,p), QSCI trials built in natural orbitals from CISD, active spaces up to \((20e,20o)\), \(10^6\) shots, and AFQMC with 512 walkers and 1000 blocks [2602.07912].

For the H\(_6\) chain, DOCI-QSCI-AFQMC on `ibm_kobe` reproduces CASCI-AFQMC and matches HCI references across the dissociation coordinate, within chemical accuracy except for slight deviations at three pre-minimum points. Simulator runs without enlargement may undersample because the distribution is skewed toward HF, but adding selected-CI expansion recovers high accuracy. DOCI-AFQMC, which uses a seniority-zero trial without Cartesian expansion, is reported to be significantly less accurate [2602.07912].

For N\(_2\), a fixed subspace without enlargement causes both DOCI-QSCI and DOCI-QSCI-AFQMC to fail to reproduce dissociation quantitatively because of insufficient sampling. After enlarged-subspace construction, DOCI-QSCI qualitatively captures dissociation, while DOCI-QSCI-AFQMC closely matches RMR-CCSD and RMR-CCSD(T) over the full curve. By contrast, single-reference RCCSD and RCCSD(T) are qualitatively incorrect in the stretched regime [2602.07912]. The comparison is important because it isolates the regime in which seniority-structured multireference information matters.

For BODIPY–O\(_2\), the reported activation and reaction energies in kcal/mol are: DOCI-QSCI\((20e,20o)\)-AFQMC, \(E_a=22.1\pm 2.4\), \(E_r=-17.5\pm 2.4\); DOCI-QSCI\((10e,10o)\)-AFQMC, \(E_a=15.9\pm 2.0\), \(E_r=-20.6\pm 2.1\); DOCI-QSCI\((20e,20o)\), \(E_a=9.1\), \(E_r=-50.3\); DOCI-QSCI\((10e,10o)\), \(E_a=18.1\), \(E_r=-35.6\); RCCSD, \(E_a=29.4\), \(E_r=-22.1\); RCCSD(T), \(E_a=3.2\), \(E_r=-18.4\); SCF, \(E_a=30.9\), \(E_r=-26.1\); RB3LYP, \(E_a=22.5\), \(E_r=-12.6\) [2602.07912]. The same study reports \(T_1\) diagnostics of 0.012 for the reactant, 0.068 for the transition state, and 0.012 for the product, with transition-state \(T_1>0.02\) indicating single-reference CC breakdown. This provides the stated rationale for why CCSD and CCSD(T) disagree qualitatively on the barrier while DOCI-QSCI-AFQMC yields reasonable reaction energetics [2602.07912].

Taken together, the benchmarks support three specific claims made in the 2026 paper: DOCI-QSCI doubles the orbital space accessible to conventional QSCI at fixed device size, the Cartesian expansion is necessary to inject seniority-breaking determinants, and subsequent ph-AFQMC post-processing delivers reasonably high accuracy [2602.07912].

## 6. Broader formulations, misconceptions, and limitations

The broader literature shows that Q-SENSE is not yet a single standardized protocol. One formulation presents it as an orthogonal, seniority-symmetry-based subspace method that interpolates between VQE and CI by constructing basis states in distinct seniority sectors, measuring \(H_{ij}=\langle \psi_i|\hat{H}|\psi_j\rangle\), and, because orthogonality is enforced by construction, solving a standard Hermitian eigenvalue problem with \(S=I\) rather than a generalized one [2509.01061]. Another formulation describes a seniority-driven operator-selection and state-preparation framework in which compact ansätze are built from ordered seniority-zero pair excitations and sparsely admitted rank-one exchanges, with subspace-expansion and unitary-ansatz realizations treated as equivalent viewpoints [2504.19760]. A review of quantum subspace methods supplies the general Rayleigh–Ritz and generalized-eigenvalue machinery into which seniority-adapted constructions naturally fit, even though that review does not explicitly define Q-SENSE [2312.00178].

A common misconception is that Q-SENSE requires exact seniority conservation by the Hamiltonian. The literature does not support that claim. The 2026 DOCI-QSCI-AFQMC realization explicitly uses a seniority-zero sampling stage precisely because the full problem contains important seniority-breaking contributions, and the orthogonal-subspace formulation states that the full Coulombic electronic Hamiltonian generally does not commute with total seniority [2602.07912, 2509.01061]. Q-SENSE is therefore better understood as exploiting seniority as a compression and selection principle, not as assuming an exact symmetry of the full electronic Hamiltonian.

The limitations are also explicit. Restricting the initial sampling to \(\Omega=0\) biases determinant selection toward pair-occupied patterns and makes open-shell and high-spin targets inaccessible within that sampling scheme. Finite-shot sampling can skew the determinant distribution, for example by overweighting HF, which degrades trial quality unless the subspace is enlarged. The Cartesian product produces \(R^2\) growth in the determinant count, and the remaining dynamical correlation may still include subtle high-seniority contributions not fully recovered by the initial expansion [2602.07912]. In the broader variants, the gains from seniority structure diminish when seniority breaking is strong enough that many inter-sector couplings must be retained [2509.01061].

The natural directions for further development are stated in the sources. The 2026 paper identifies adaptive expansion strategies, improved ansätze or ADAPTively optimized sampling preparations, nonorthogonal multi-reference trials, more sophisticated selection criteria beyond \(\max_k |H_{lk} c_k|\), tighter integration with tailored CC or perturbative corrections, and extensions beyond \(\Omega=0\) for open-shell targets [2602.07912]. This suggests that Q-SENSE is evolving toward a modular framework in which seniority organizes the quantum stage, while increasingly flexible classical or projector-based corrections restore the parts of correlation that seniority-zero sampling cannot represent on its own.

Source: https://www.emergentmind.com/topics/quantum-seniority-based-subspace-expansion-q-sense