---
title: Hard-Core Fermi-Hubbard Model
url: https://www.emergentmind.com/topics/hard-core-fermi-hubbard-model
type: topic
---

# Hard-Core Fermi-Hubbard Model

The hard-core Fermi-Hubbard model is the strong-coupling limit of the spinful Fermi-Hubbard Hamiltonian in which on-site double occupancy is energetically suppressed, so the low-energy Hilbert space is restricted from \(\{|0\rangle,|\uparrow\rangle,|\downarrow\rangle,|\uparrow\downarrow\rangle\}\) to a constrained sector dominated by \(|0\rangle,|\uparrow\rangle,|\downarrow\rangle\). For repulsive interactions at half filling, this is the Mott-insulating regime: charge fluctuations are frozen, local moments remain, and the low-energy theory reduces to a spin-\(\tfrac12\) Heisenberg antiferromagnet with superexchange \(J=4t^2/U\) [1007.0012]. In contemporary cold-atom realizations, this regime is accessed in homogeneous three-dimensional optical lattices with \(U/t\) up to \(\sim 18\), where doublon fractions become rare and extended antiferromagnetic correlations and the antiferromagnetic phase transition can be resolved directly [2502.08104; 2402.14605].

## 1. Microscopic formulation

The standard spinful Fermi-Hubbard Hamiltonian relevant to the hard-core limit is
\[
\hat H = -t \sum_{\langle i,j\rangle,\sigma} \left( \hat c^{\dagger}_{i\sigma} \hat c_{j\sigma} + \mathrm{h.c.} \right)
+ U \sum_i \hat n_{i\uparrow} \hat n_{i\downarrow}
+ \sum_{i,\sigma} (V_i-\mu)\hat n_{i\sigma},
\]
or, in homogeneous bulk realizations, without a dominant site-dependent confinement term,
\[
H = -t \sum_{\langle i,j \rangle,\sigma} \left( c^{\dagger}_{i\sigma} c_{j\sigma} + \mathrm{h.c.} \right)
+ U \sum_i n_{i\uparrow} n_{i\downarrow}
- \mu \sum_i (n_{i\uparrow}+n_{i\downarrow}).
\]
Here \(t\) is the nearest-neighbor tunneling matrix element, \(U\) is the on-site interaction between opposite spins, and \(\mu\) is the chemical potential. In optical lattices, \(t\) originates from Wannier-function overlap between neighboring wells, while \(U=g\int d^3r\,|w(\mathbf r)|^4\) with \(g=4\pi\hbar^2 a/m\); \(a\) is tuned through a Feshbach resonance \(a(B)=a_{\rm bg}\left(1-\frac{\Delta}{B-B_0}\right)\). Three orthogonal standing waves generate
\[
V_{\rm lat}(\mathbf r)=V_{0x}\cos^2(kx)+V_{0y}\cos^2(ky)+V_{0z}\cos^2(kz),
\]
with lattice spacing \(d=\lambda/2\) and recoil energy \(E_r=\hbar^2k^2/2m\); increasing lattice depth exponentially suppresses \(t\), and in the lowest band the next-nearest-neighbor tunneling is suppressed by about an order of magnitude [1007.0012].

In trapped realizations, the last term produces spatially varying local chemical potential \(\mu_{\rm local}(\mathbf r)=\mu-V_{\rm trap}(\mathbf r)\), typically with harmonic confinement. In homogeneous flat-top plus box-trap implementations, the central region instead realizes an approximately uniform \(\mu\), \(t\), and \(U\), which removes the dominant LDA convolution present in Gaussian lattices and makes direct comparison to homogeneous Hubbard theory possible over \(\sim 8\times 10^5\) sites [2502.08104].

## 2. Strong-coupling projection and effective models

In the repulsive model, “hard-core” means \(U/t\gg 1\): a doublon costs energy \(U\), so low-energy states exclude simultaneous occupation by opposite spins on the same site. This does not change fermionic statistics; it truncates the local Hilbert space by suppressing \(|\uparrow\downarrow\rangle\). At half filling, the system becomes a Mott insulator: it is incompressible, double occupancy is suppressed, and each site carries a localized spin degree of freedom. The low-energy effective Hamiltonian is
\[
\hat H_{\rm spin}=J\sum_{\langle i,j\rangle}\mathbf S_i\cdot \mathbf S_j,\qquad J=\frac{4t^2}{U},
\]
obtained from second-order virtual hopping processes that momentarily create a doublon and favor antiferromagnetic correlations [1007.0012].

A recurrent misconception is to identify the hard-core limit with a band insulator. The review explicitly distinguishes them: both are incompressible, but the Mott insulator is interaction-driven, suppresses doublons, and can carry large spin entropy, whereas a band insulator has maximal double occupancy and interactions play a minor role [1007.0012].

Away from half filling, a plausible low-energy description is a \(t\)–\(J\)-type constrained-hopping model, although the review emphasizes primarily the Heisenberg limit and Mott physics. In the strong-coupling half-filled case on planar graphs, the same projection yields the \(t\)–\(J\) model with no double occupancy, and two-particle scattering is reflection-only with no transmission, “as would be the case for classical hard spheres,” which underlies the universality construction based on hard-core scattering [1409.3585].

The attractive branch supplies a complementary hard-core limit. For \(U<0\) and large \(|U|\), opposite-spin fermions bind into on-site doublons; these pairs “can be regarded as hard-core bosons,” with pair hopping dominated by second-order tunneling \(t^2/|U|\) and repulsive nearest-neighbor interaction also proportional to \(t^2/|U|\) [1007.0012]. In the mass-imbalanced extended Fermi-Hubbard model with nearest-neighbor interaction \(V\), the doublon subspace is governed by
\[
\hat H_{\rm eff}
= \sum_{\langle i,j\rangle}\left[
J_{\rm eff}(\hat d_i^\dagger \hat d_j + \hat d_j^\dagger \hat d_i)
+ \frac{1}{2}V_{\rm eff}\hat n_i\hat n_j
\right],
\]
with
\[
J_{\rm eff}=\frac{2J_\uparrow J_\downarrow}{U},
\qquad
V_{\rm eff}=V-\frac{J_\uparrow^2+J_\downarrow^2}{U}.
\]
At the resonance point \(V_{\rm eff}=0\), doublons behave like free hard-core bosons; under a tilt they exhibit Bloch oscillations with period \(\pi/F\), while bound doublon pairs yield period \(\pi/(2F)\) [2408.08194].

## 3. Realization in optical lattices and homogeneous traps

The canonical realization uses ultracold fermions in optical lattices. Deepening the lattice reduces \(t\); tuning the scattering length through a Feshbach resonance increases \(U\); together these drive the system into the large-\(U/t\) regime. Three-dimensional simple-cubic lattices, 2D planes, 1D tubes, superlattices, and triangular lattices are all obtained by modifying lattice depths along selected axes, so dimensionality and frustration are directly engineered rather than emergent [1007.0012].

The homogeneous realization reported in 2025 combines a flat-top optical lattice with a cylindrical optical box. Each axis is formed from two counter-propagating 1064 nm flat-top beams, producing an almost uniform central disk of diameter \(\sim 60\,\mu\mathrm m\), while blue-detuned 532 nm ring and line beams isolate the central region. The resulting 3D simple-cubic sample contains \(\sim 8\times 10^5\) lattice sites, with a flat density profile over \(\sim 40\,\mu\mathrm m\) diameter. At lattice depth \(V_0=6.25E_r\), with \(E_r/h=29.30\,\mathrm{kHz}\), the hopping is \(t/h=1.40(1)\,\mathrm{kHz}\); the interaction is tuned over \(0\le U/h\le 24.67(8)\,\mathrm{kHz}\), corresponding to \(0\lesssim U/t\lesssim 18\) [2502.08104].

This homogeneous geometry is consequential for the hard-core regime. In harmonic traps, metallic, Mott-insulating, and band-insulating shells coexist. In the flat-top plus box-trap system, by contrast, the bulk Hamiltonian is essentially uniform, so the metal-to-Mott crossover, Pomeranchuk physics, and AFM correlations can be read out without dominant trap averaging [2502.08104].

## 4. Phases and experimental signatures

In trapped repulsive systems, the phase structure is shell-like. Low atom number gives a metallic cloud; at intermediate atom number a Mott core with density \(\sim 1\) appears at the center and metallic shells survive outside; when the central chemical potential exceeds \(U\), double occupancy rises again and eventually a band insulator with two atoms per site forms in the core. The defining hard-core signature is the suppression of doubly occupied sites in the Mott region, not merely incompressibility [1007.0012].

Double occupancy is the central diagnostic. In early optical-lattice experiments, a rapid lattice-depth ramp froze motion, a Feshbach resonance shifted the energy of doubly occupied sites, and an RF pulse transferred one partner to an empty spin state, making doublons countable by absorption imaging. In the Mott regime, double occupancy below \(2\%\) was measured, with estimated hole fraction below \(3\%\) and \(k_BT/U\approx 0.11\); cloud-size measurements simultaneously indicated an incompressible core [1007.0012].

In the homogeneous implementation, the doublon fraction is defined by
\[
D=\frac{2\langle n_\uparrow n_\downarrow\rangle}
{\langle n_\uparrow\rangle+\langle n_\downarrow\rangle}.
\]
The measurement freezes tunneling by ramping from \(6.25E_r\) to \(20E_r\) in \(20\,\mu\mathrm s\), then uses a Gaussian-shaped microwave pulse of total duration \(800\,\mu\mathrm s\) to drive \(|1\rangle\rightarrow|2\rangle\), producing distinct resonance peaks for atoms on singly occupied and doubly occupied sites. At half filling and entropy per particle \(s\simeq 0.3k_{\mathrm B}\), \(D(U)\) decreases smoothly with increasing \(U/t\); around \(U/t\simeq 8\), \(D\sim 0.5\%\), and for the coldest samples at \(U=9.97(18)t\), \(D\) reaches \(\sim 0.08\) [2502.08104].

Modulation spectroscopy probes the high-energy sector excluded by the hard-core constraint. Periodic modulation \(V_0(t)=V_0+\delta V\sin\omega t\) creates doublon-hole excitations resonantly at \(\hbar\omega\approx U\), directly measuring the Hubbard gap and the separation between lower and upper Hubbard bands. Doublon lifetimes scale exponentially with \(U/6t\), consistent with the need to redistribute energy \(U\) into many low-energy excitations [1007.0012].

## 5. Antiferromagnetism, Pomeranchuk physics, and criticality

The hard-core Fermi-Hubbard model is not exhausted by doublon suppression; its central low-energy content is quantum magnetism. In the 3D repulsive model at half filling, the Néel transition interpolates between a weak-coupling SDW instability and the strong-coupling Heisenberg regime. In the latter, the paper uses the known Heisenberg result
\[
T_N \approx 0.9575\,J,\qquad J=\frac{4t^2}{U},
\]
so the AFM transition is set by superexchange rather than by the Hubbard scale \(U\) [2402.14605].

The homogeneous 2025 experiment inferred this crossover already from doublon thermodynamics. At fixed \(U=5.98(7)t\), \(D(s)\) is non-monotonic: for high entropy \(s\gtrsim 0.9k_{\mathrm B}\), cooling reduces \(D\); around \(s\simeq 0.9k_{\mathrm B}\), \(D\) reaches a minimum; for \(0.4k_{\mathrm B}\lesssim s\lesssim 0.9k_{\mathrm B}\), \(D\) increases upon further cooling, identifying the Pomeranchuk effect; and for \(s\lesssim 0.4k_{\mathrm B}\), \(D\) turns downward again, signaling the development of extended antiferromagnetic correlations. The same qualitative structure appears at \(U=7.92(14)t\), while at \(U=9.97(18)t\) thermally induced doublons are so suppressed that the remaining doublons are primarily quantum fluctuations \(\sim t^2/U\) [2502.08104].

The AFM phase transition itself was reported in 2024 in a 3D uniform optical lattice of approximately \(8\times10^5\) sites. The spin structure factor at
\[
\mathbf Q_\pi=\frac{2\pi}{a}\left(\tfrac12,\tfrac12,\tfrac12\right)
\]
was measured via spin-sensitive Bragg scattering and defined as
\[
S_{\boldsymbol{\pi}}
=\frac{4}{N}\sum_{i,j} e^{i\boldsymbol{\pi}\cdot(\mathbf R_i-\mathbf R_j)}
\langle \hat S_i^z \hat S_j^z\rangle.
\]
Near criticality, the data were fitted with
\[
S_{\boldsymbol{\pi}}-S_0 \propto |p/p_c-1|^{-\gamma},\qquad \gamma=1.396,
\]
for \(p=U\), \(s\), or \(n\), consistent with the Heisenberg universality class. At half filling and near optimal coupling, the measured spin structure factor reached \(123(8)\), and the AFM dome in interaction space was bounded by \(U_{c,s}=8.84(33)t\) and \(U_{c,l}=14.86(42)t\). At \(U/t\simeq 11.75\), the critical entropy was \(s_c=0.27(1)k_B\); at \(U/t\simeq 8.21\), \(s_c=0.20(1)k_B\). Doping suppressed AFM order rapidly, with fitted critical fillings \(n_{c,s}=0.95(1)\) and \(n_{c,l}=1.07(1)\) [2402.14605].

## 6. Computational formulations and broader significance

The hard-core limit is also a central organizing principle in quantum simulation and Hamiltonian complexity. A qudit formulation based on a Qudit Fermionic Mapping encodes the full spinful Hubbard site in one ququart with local basis \(|0\rangle,|\uparrow\rangle,|\downarrow\rangle,|\uparrow\downarrow\rangle\); the hard-core limit is then enforced either by large \(\upsilon/J\) so that the fourth level is energetically inaccessible, or by restricting to the \(|0\rangle,|\uparrow\rangle,|\downarrow\rangle\) subspace. This preserves nearest-neighbor connectivity and makes the on-site interaction a single-qudit operation [2402.01243].

Ground-state preparation protocols have also been adapted to 2D hardware. QETU-based preparation of the \(2\times2\) Fermi-Hubbard model on a 9-qubit grid uses controlled time evolution, fermionic swap networks, and a 2D topology tailored to the Jordan-Wigner layout; the same circuit structure remains applicable in the hard-core limit, where on-site \(Z_{i,\uparrow}Z_{i,\downarrow}\) terms are removed or assigned a very large penalty [2411.18535]. A distinct “Fermi Machine” construction maps Hubbard models to multi-component noninteracting fermions and, in the large-\(U\) limit, can enforce perfect exclusion of double occupancy through the hidden-fermion structure [2407.20546].

At the complexity-theoretic end, planar-graph Hubbard dynamics at strong coupling and half filling reduce to the \(t\)–\(J\) model, where only reflection and no transmission occur in scattering events; this reflection-only hard-core dynamics suffices to encode universal quantum computation [1409.3585]. For parameter identification rather than state preparation, a Hamiltonian-learning protocol for bounded-degree Fermi-Hubbard graphs achieves \(\tilde{\mathcal O}(\epsilon^{-1})\) total evolution time and uses only simple one- or two-site fermionic manipulations, which is directly relevant to calibrating analog hard-core Hubbard simulators [2312.17390].

Alternative solid-state simulators reach the same qualitative regime. A semiconductor triple-dot array realizing an extended Fermi-Hubbard Hamiltonian achieved \(U/t\approx 7.1(4)\) and \(t/k_BT\approx 54(5)\), allowing detailed characterization of the collective Coulomb blockade transition as the finite-size analogue of the interaction-driven Mott metal-to-insulator transition [1702.07511]. This suggests that the hard-core Fermi-Hubbard model is best viewed not as a single limiting Hamiltonian but as a unifying low-energy description of constrained charge motion, Mott insulation, superexchange magnetism, and, in attractive or extended variants, hard-core doublon dynamics across cold-atom, solid-state, and digital-quantum platforms.

Source: https://www.emergentmind.com/topics/hard-core-fermi-hubbard-model