Barkoutsos et al. put the transformation before the circuit: they "propose a transformation of the electronic structure Hamiltonian in the second quantization framework into the particle-hole (p/h) picture, which offers a better starting point for the expansion of the trial wavefunction", so that the state "is parametrized in a way to efficiently explore the sector of the molecular Fock space that contains the desired solution". The circuit family follows from that choice — "a new family of quantum circuits based on exchange-type gates that enable accurate calculations while keeping the gate count (i.e., the circuit depth) low" — and the resulting method is named in the paper: "the particle-hole implementation of the Unitary Coupled Cluster (UCC) method within the Variational Quantum Eigensolver approach ... named q-UCC". The claim to read carefully is the last one, because it is about Trotterisation rather than about accuracy in general: they "show how a single Trotter step can accurately and efficiently reproduce the ground state energies of simple molecular systems", and *simple* is the paper's own word. No hardware run is reported.

The Hamiltonian whose ground state is wanted, together with whatever structure is to be respected — particle number, spin, point-group symmetry, a reference determinant — and the connectivity and native gate set of the device the family has to run on.

1. Hamiltonian whose eigenvalues are wanted → Parameterised circuit familyexponentiate excitations, Trotter-split

   The p/h picture redefines creation/annihilation operators as $\hat b_i^\dagger=\hat a_i$ for holes and $\hat b_m^\dagger=\hat a_m^\dagger$ for particles (Eqs. 5-8), so normal-ordering $\hat N_b[\cdot]$ against $|\Phi_0\rangle$ gives $\hat H^{p/h}=E_{HF}+\sum_{rs}\langle r|\hat F|s\rangle \hat N_b[\hat a_r^\dagger \hat a_s]+\tfrac12\sum_{rstu}\langle rs|\hat g|tu\rangle \hat N_b[\hat a_r^\dagger \hat a_s^\dagger \hat a_u \hat a_t]$ (Eq. 9). The end state is $|\Psi(\vec\theta)\rangle=e^{\hat T(\vec\theta)-\hat T^\dagger(\vec\theta)}|\Phi_0\rangle$, with $\hat T=\hat T_1+\hat T_2$ parameterised by amplitudes $\theta_i^m,\theta_{i,j}^{m,n}$ in $\hat T_1=\sum_{i;m}\theta_i^m\hat a_m^\dagger \hat a_i$ and $\hat T_2=\tfrac12\sum_{i,j;m,n}\theta_{i,j}^{m,n}\hat a_n^\dagger \hat a_m^\dagger \hat a_j \hat a_i$. approximation: the circuit realises this exponential by a first-order Trotter split $e^{\hat A+\hat B}=\lim_{n\to\infty}(e^{\hat A/n}e^{\hat B/n})^n$ (Eq. 16) truncated at a finite number of steps n. assumption: the redefined operators still obey fermionic anti-commutation, so $\hat N_b$ defines a valid quasiparticle vacuum at $|\Phi_0\rangle$.

   approximationassumption

A circuit family with a fixed structure and free real parameters, together with the number of those parameters — which is the size of the classical search problem handed to the next layer.

given  electronic Hamiltonian in second quantization, HF orbital basis
           H_el = sum_ij h_ij a_i^dag a_j + sum_ijkl g_ijkl a_i^dag a_j^dag a_l a_k    (Eq. 1)
       HF reference determinant  |Phi_0> = prod_{i=1}^N a_i^dag |vac>                    (Eq. 4)
           equivalently the bit-string ket |1...1 0...0>, N occupied orbitals set to 1
                                                        (unlabeled form, Sec. III, step iii)
       index convention:  i,j,k,l occupied;  m,n,p,q virtual;  r,s,t,u either
       fermion-to-qubit map: Jordan-Wigner                                    (Sec. III, step ii)

requires  cluster operator truncated at singles+doubles,  T = T1 + T2
    # higher T_n need only longer circuits, not new theory -- the paper restricts
    #   its own implementation to UCCSD as "unpractical" beyond that (Sec. II B 1)

# --- rewrite the Hamiltonian in the particle-hole picture ----------------------
redefine   b_i^dag = a_i,       b_i = a_i^dag      (hole creation / annihilation)    (Eqs. 5, 7)
           b_m^dag = a_m^dag,   b_m = a_m          (particle creation / annihilation) (Eqs. 6, 8)
    # b's still obey fermionic anti-commutation -- N_b[.] then defines a valid
    #   quasiparticle vacuum at |Phi_0>
normal-order  N_b[.]  against |Phi_0>  to get
    H_p/h = E_HF + sum_rs <r|F|s> N_b[a_r^dag a_s]
                  + (1/2) sum_rstu <rs|g|tu> N_b[a_r^dag a_s^dag a_u a_t]         (Eq. 9)
    # this is the operator the device measures -- the qubit register itself is
    #   not re-encoded, only the Hamiltonian is rewritten

# --- build the cluster operator --------------------------------------------------
T1 = sum_{i;m}           theta_i^m         a_m^dag a_i                              (Eq. 13)
T2 = (1/2) sum_{i,j;m,n} theta_{i,j}^{m,n} a_n^dag a_m^dag a_j a_i                  (Eq. 14)
    # sums run over ALL occupied/virtual pairs -- an active-space restriction to
    #   primed subsets i',m' shrinks them directly, same functional form (Eqs. 23-24)
|Psi(theta)> = exp( (T1 - T1^dag) + (T2 - T2^dag) ) |Phi_0>                         (Eq. 12)

# --- exponentiate by first-order Trotter splitting, n steps ----------------------
exp( (T1-T1^dag) + (T2-T2^dag) )
    ~=  [ exp((T1-T1^dag)/n)  exp((T2-T2^dag)/n) ]^n                                (Eq. 25)
    # exact only as n -> infinity, general two-term case  (Eq. 16)
    # each factor realised by the excitation circuits of Fig. 1;
    #   the underlying Pauli-string form is given in App. A (Eq. AA.6)

# --- optimise theta by VQE, at the chosen n --------------------------------------
loop
    prepare |Psi(theta)>  on the device via the n-step circuit above
    measure  <Psi(theta)| H_p/h |Psi(theta)>
    theta <- BFGS step on the measured energy                                       (Sec. III)
until convergence

return  |Psi(theta_min)>,   E_GS = E_HF + E_corr_UCCSD(theta_min)                    (Eq. 15)
# n=1, WITH theta fully re-optimised by VQE AT that n, measures error < 1e-10 Ha
#   at equilibrium H2 (Eq. 27; Fig. 5 red triangles) -- a DIFFERENT, larger
#   'Trotter error' (Eq. 26, amplitudes fixed at their exact-exponentiation
#   values) is plotted in the same figure for comparison and is NOT the
#   quantity claimed small here. The paper calls the n=1 result "surprising",
#   explained by VQE's flexibility 'absorbing' the error, and explicitly not
#   a proven bound                                          (Sec. IV B 2; Sec. V)

Barkoutsos et al. give one scaling: q-UCCSD circuit depth as $\mathcal{O}\left( \binom{N_{\mathrm{occ}}}{2} \times \binom{N_{\mathrm{virt}}}{2} \times N_{\mathrm{qubits}} \right)$ — binomial coefficients — for the cluster operator truncated at $\hat{T}=\hat{T}_1+\hat{T}_2$, where $N_{\mathrm{occ}}$ and $N_{\mathrm{virt}}$ count only orbitals taking part in the excitations, so an active space shrinks it directly. Everything else is measured on $\mathrm{H_2}$ in 6-31G (8 qubits, no Hamiltonian reduction), in simulation, with nearest-neighbour connectivity and exchange gates counted as one native two-qubit operation. q-UCCSD at a single Trotter step: 15 parameters, 708 one-qubit and 608 two-qubit gates, identical under both Hamiltonians — the particle-hole rewrite buys optimiser work, BFGS iterations falling 53 to 27 and Pauli-string evaluations $3383 \times N_s$ to $1471 \times N_s$ with $N_s$ never fixed, an overall speed-up of "about a factor 2 to 3". For the exchange-gate heuristic ansatz, reaching $10^{-7}$ Ha at equilibrium takes $D=6$ blocks (84 parameters, 42 two-qubit gates, no one-qubit gates) against $D=8$ (112, 56); demanding chemical accuracy ($0.5 \times 10^{-2}$ Ha) along the whole dissociation path erases that saving — $D=10$, 140, 70 under both. Nothing bounds the VQE iteration count; $D$ is fixed empirically. No hardware run.

• Particle-conserving VQE circuits · Qiskit

  From the repository — run, not written up from a paper · unsupported

  About

  Givens-style or excitation-preserving blocks keep evolution inside a fixed-particle-number sector.

  Methods

  None found yet.

  Data

  None found yet.

  Code

  Qiskit

  Results

  Literature-backed method record; algorithmic scope and evidence boundary reviewed, with no benchmark run claimed.

• Particle-conserving VQE circuits

  Givens-style or excitation-preserving blocks keep evolution inside a fixed-particle-number sector.