← Atlas

Operators · Fermionic Hamiltonians

Fermi-Hubbard dimer operator (Jordan-Wigner encoded)

The two-site Fermi-Hubbard dimer Hamiltonian, mapped to qubits via the Jordan-Wigner transformation: nearest-neighbor hopping competing with on-site Coulomb repulsion.

Exact & formalUpdated

View frozen graph

1↑1↓2↑2↓Step 1 of 4 · Hop (XX+YY) · on 1↑, 2↑Hop (XX+YY)Step 2 of 4 · Hop (XX+YY) · on 1↓, 2↓Hop (XX+YY)Step 3 of 4 · U n₁↑n₁↓ · on 1↑, 1↓U n₁↑n₁↓Step 4 of 4 · U n₂↑n₂↓ · on 2↑, 2↓U n₂↑n₂↓

4 steps across 4 wires. Hover a step, or press play to go through them in order.

Expected outcomes

  • Half-filling ground sector (N=2)100%
  • Steps4
  • Wires4
  • Sites × spins2 × 2 = 4 modes / qubits
  • Terms2 hopping (XX+YY) + 2 on-site U
  • ConservedTotal fermion number N
  • Papers2

The single-band Fermi-Hubbard model on a lattice with hopping tt and on-site interaction UU is

H=−t∑⟨i,j⟩,σ(ciσ†cjσ+cjσ†ciσ)+U∑ini↑ni↓,niσ=ciσ†ciσ.H = -t\sum_{\langle i,j\rangle,\sigma} \big(c^\dagger_{i\sigma} c_{j\sigma} + c^\dagger_{j\sigma} c_{i\sigma}\big) + U\sum_i n_{i\uparrow} n_{i\downarrow}, \qquad n_{i\sigma} = c^\dagger_{i\sigma}c_{i\sigma}.

Two-site dimer. With 2 sites and 2 spins there are 4 fermionic modes (1↑,1↓,2↑,2↓)(1{\uparrow},1{\downarrow},2{\uparrow},2{\downarrow}), mapped to 4 qubits. Using the Jordan-Wigner transform with mode ordering 1↑,1↓,2↑,2↓1{\uparrow},1{\downarrow},2{\uparrow},2{\downarrow},

ck†=(∏l<kZl)Xk−iYk2,ck=(∏l<kZl)Xk+iYk2,c^\dagger_k = \Big(\prod_{l<k} Z_l\Big)\frac{X_k - iY_k}{2}, \qquad c_k = \Big(\prod_{l<k} Z_l\Big)\frac{X_k + iY_k}{2},

so a hopping term ck†ck+1+h.c.c^\dagger_k c_{k+1} + \text{h.c.} becomes 12(XkXk+1+YkYk+1)\tfrac12(X_kX_{k+1}+Y_kY_{k+1}) acting on adjacent-index modes (with an extra ZZ-string factor if the modes are not adjacent in the chosen ordering), and each density term niσ=12(I−Ziσ)n_{i\sigma} = \tfrac12(I - Z_{i\sigma}) is diagonal.

Particle-number conservation. Both the hopping and interaction terms conserve total fermion number N=∑iσniσN=\sum_{i\sigma} n_{i\sigma}: hopping moves a fermion between sites without creating or destroying one, and the interaction term is purely diagonal in occupation. Consequently [H,N]=0[H, N] = 0, which block-diagonalizes the Hamiltonian by total filling and is the symmetry sector structure any correct encoding must reproduce.

Physical content. For U/t→0U/t \to 0 the model reduces to free fermions hopping on the lattice (band physics); for U/t→∞U/t \to \infty double occupancy is forbidden and, at half filling, the low-energy physics maps onto the antiferromagnetic Heisenberg model via a Schrieffer-Wolff-type superexchange argument (Jeff=4t2/UJ_{\text{eff}} = 4t^2/U) — directly connecting this operator to the Heisenberg XXZ operator elsewhere in this catalog. The model is believed (though not proven in 2D) to host a dd-wave superconducting regime near intermediate coupling, which is the primary reason it remains a central target for near-term quantum simulation.

References

Exact & formal
  • Direct mathematics
  • Verified by construction
  • Peer-reviewed paper
Method
Total particle number N = Σ n_iσ (per-spin-orbital occupation from Jordan-Wigner) was checked to commute exactly with the assembled 4-qubit (2-site, 2-spin) SparsePauliOp Hamiltonian by explicit matrix commutator, confirming the hopping and interaction terms both conserve fermion number as required by the JW encoding.
Result
Pass · [H, N] is the zero matrix to numerical precision, confirming particle-number conservation of the encoded operator.
Strategies for solving the Fermi-Hubbard model on near-term quantum computers ↗
Kind
curated reference
Reviewed by
Leona Quantum curation pass
License
CC BY 4.0-compatible reference metadata
hubbard modelfermionic simulationjordan-wignerstrongly correlated electrons

Atlas stars stay in the public catalog. Saving this entry to your workspace starts an unstarred private copy.