Fermionized Photonic TG Gas

• Tonks–Girardeau gas of fermionized photons is a 1D system where engineered strong photon-photon interactions transform bosons into hard-core particles exhibiting free fermion properties.
• Methodologies such as Bose–Fermi mapping, DMRG, and determinant techniques precisely predict experimental signatures like perfect antibunching and Friedel oscillations in correlation functions.
• Experimental setups in cavity QED, waveguides, and superconducting circuits enable controlled photon-photon repulsion, offering a platform for observing strongly correlated many-body quantum effects.

A Tonks-Girardeau (TG) gas of fermionized photons refers to a one-dimensional (1D) system where photons, through engineered strong interactions, behave as impenetrable (hard-core) bosons, and thereby exhibit many-body correlations characteristic of free fermions. The term derives from the original TG model for bosonic atoms with infinite delta-function repulsion, but is realized here in photonic systems through effective photon-photon interactions—typically mediated by nonlinearity in cavity QED architectures, waveguides, or superconducting transmission lines. The physical and analytic framework is underpinned by the Bose–Fermi mapping, enabling exact solutions and the prediction of experimentally robust “fermionization” signatures such as perfect antibunching and Friedel oscillations in correlation functions.

1. Physical Foundations and Hamiltonian Modeling

The TG regime for photons can be reached when effective photon-photon repulsion is large compared to kinetic hopping and the density is low. In a system of coupled optical cavities described by the Jaynes-Cummings-Hubbard (JCH) model, the effective Hamiltonian (at zero detuning, ℏ ≈ 1) is

$$H_{\mathrm{JCH}}/g = -(\kappa/g) \sum_{\langle i,j\rangle}(a_i^\dagger a_j + \text{h.c.}) + \sum_i [ a_i^\dagger \sigma_i^- + \text{h.c.} - ((\mu - \omega_c)/g) N_i ],$$

where $a_i$ is the photonic annihilation operator, $\sigma_i^-$ the atomic lowering operator, $\kappa$ the cavity photon hopping, $g$ the Jaynes-Cummings coupling, and $N_i$ the total excitation number on site $i$. The effective onsite photon-photon repulsion is $U_{\mathrm{eff}} \simeq (2-\sqrt{2})g$ due to the JC anharmonicity.

A TG-like phase emerges for $\kappa / U_{\mathrm{eff}} \ll 1$ and low densities $\rho \ll 1$. In practice, for $\kappa/g \lesssim 0.1$ and $\rho \lesssim 0.5$, observables from density-matrix renormalization group (DMRG) calculations quantitatively match TG predictions (D'Souza et al., 2013).

In nonlinear transmission lines with embedded Josephson junctions, the microscopic circuit Lagrangian leads, via coarse-graining and SVEA, to an effective Lieb–Liniger Hamiltonian,

$$H = \int d\zeta \left[ \Psi^\dagger(\zeta) \left(-\frac{\hbar^2}{2m}\partial_\zeta^2 \right) \Psi(\zeta) + \frac{g}{2} (\Psi^\dagger(\zeta))^2 (\Psi(\zeta))^2 \right]$$

where $\Psi(\zeta)$ is the bosonic field operator, and the dimensionless parameter $\gamma = mg/(\hbar^2 \rho)$ controls the crossover from weakly to strongly interacting regimes. $\gamma \gg 1$ realizes the TG (hard-core) limit (Fatis et al., 12 Jan 2026).

2. Bose–Fermi Mapping and Wavefunction Structure

The distinctive property of the TG regime is the exact mapping between the bosonic many-body ground state and the ground state of free fermions. The Girardeau mapping prescribes

$$\Psi_B(x_1,\dots,x_N) = A(x_1,\dots,x_N) \Psi_F(x_1,\dots,x_N),$$

where $$A(x_1,\dots,x_N) = \prod_{i<j} \mathrm{sgn}(x_i - x_j)$$ enforces bosonic symmetry, and $\Psi_F$ is a Slater determinant of single-particle fermionic orbitals.

On discrete lattices, Jordan-Wigner transformation provides the operator mapping: $$a_j = \exp\left(i\pi \sum_{\ell < j} c^\dagger_\ell c_\ell\right) c_j,$$ where $c_j$ are fermionic operators, yielding $(a_j)^2 = 0$ in the hard-core (TG) limit (D'Souza et al., 2013). This structure underpins the exact calculation of correlations and emergence of fermionic statistics in bosonic observables.

3. Correlation Functions and Spectral Properties

One- and Two-Body Correlations

The first-order correlation,

$g^{(1)}(r) = \langle a_i^\dagger a_{i+r} \rangle,$

exhibits a characteristic TG long-distance decay, $g^{(1)}(r) \sim r^{-1/2}$ for $1 \ll r \ll L$ on a ring of length $L$, more precisely: $g^{(1)}(r)\simeq [2 n L |\sin(\pi r/L)|]^{-1/2}$. DMRG simulations confirm this scaling with logarithmic fits deriving an exponent $\alpha \approx 1/2$ over a range of $\kappa/g$ in the JCH system (D'Souza et al., 2013).

The normalized second-order correlation,

$$g^{(2)}(r) = \frac{\langle a_i^\dagger a_{i+r}^\dagger a_{i+r} a_i \rangle}{\langle n \rangle^2},$$

takes the analytic form

$$g^{(2)}(r) = 1 - \left[\frac{\sin(\pi n r)}{\pi n r}\right]^2,$$

displaying a “correlation hole” $g^{(2)}(0) = 0$ (perfect antibunching) and Friedel oscillations of period $\lambda_F = 1/n$ due to effective Pauli exclusion (D'Souza et al., 2013, Fatis et al., 12 Jan 2026, Hao et al., 2022).

Spectral Function and Finite Temperature

In the TG limit, all finite-temperature and trap-dependent single-particle Green’s functions, $g^{(\pm)}$, can be expressed as Fredholm determinants constructed from single-particle wavefunctions and Fermi–Dirac weights (Patu, 2022). The spectral function $A(k,\omega)$ is obtained via Fourier transforms of the retarded Green’s function and reflects free-fermion branches and lattice effects. In the presence of harmonic confinement, the number of observable singular lines in $A(k,\omega)$ reduces, offering a spectroscopically sharp signature of Mott versus superfluid regions (Patu, 2022).

Higher-Order Correlations

The full $n$-body reduced density matrix for the TG state can be expressed as a Toeplitz determinant of order $N-n$, with explicit analytic expressions for two- and three-body reduced density matrices (Hao et al., 2022). These encode the many-body structure, with vanishing of the correlations when any two coordinates coincide (hard-core constraint). Natural orbital decomposition provides the eigenfunctions and occupation numbers relevant for BCS-type descriptions (Sabater et al., 2024).

4. Experimental Realizations and Observable Signatures

Cavity-QED, Waveguides, and Superconducting Platforms

TG physics has been theoretically and numerically established in coupled Jaynes–Cummings cavities via DMRG (D'Souza et al., 2013). In propagating geometries, superconducting transmission lines embedded with Josephson junctions can be spatially tapered to adiabatically transform an incident coherent field into a TG gas. Critical to this process is the controlled tuning of the effective interaction parameter $\gamma$ through circuit parameters (Josephson energy, capacitance, impedance, photon density) and system length (Fatis et al., 12 Jan 2026).

Measurement Protocols

Distinctive experimental signatures are anticipated in photonic correlation measurements:

• $g^{(2)}(0) \approx 0$, indicating perfect antibunching in the output photon statistics;
• Friedel oscillations in $g^{(2)}(\Delta t)$, with period $1/\rho$ in the temporal domain;
• Algebraic decay $$g^{(1)}(\Delta t) \sim (\pi \rho |\Delta t|)^{-1/2}$$ of first-order coherence;
• Power-law divergent momentum distribution $n(k) \sim |k|^{-1/2}$ near $k \to 0$.

These can be probed by standard photon correlation setups (e.g., coincidence counting, quadrature detection) in both microwave and optical domains (Fatis et al., 12 Jan 2026, D'Souza et al., 2013).

5. Quantum Entanglement, Pairing, and BCS-Analogues

The TG regime supports nontrivial pairing and entanglement in both atomic and photonic arenas. For the “fermionic TG gas,” the ground state is exactly recast as a number-conserving BCS (Bardeen-Cooper-Schrieffer) wavefunction, constructed from pair-creation operators on natural orbitals of the one-body density matrix (Sabater et al., 2024). The expectation values of pairing observables, such as

$$\langle a_k^\dagger a_{-k}^\dagger a_{-\ell} a_\ell + \text{h.c.} \rangle < 0,$$

are strictly negative, providing an unambiguous experimental witness for quantum pairing even in the absence of conventional superconductivity.

For TG photons, the same BCS analytic machinery applies: natural photonic orbitals and occupation numbers extracted from the one-body density matrix allow construction of pair-creation operators, and measurement of $g^{(2)}$ and related witnesses yields direct evidence for paired fermionized photonic states (Sabater et al., 2024).

6. Numerical and Analytic Methodologies

Computational studies of TG gases—atomic or photonic—benefit from the exact solvability of the model via Bose–Fermi mapping. Key tools include:

• DMRG for ground-state properties of lattice and cavity systems (D'Souza et al., 2013);
• Determinant formulae (Toeplitz and Fredholm) for arbitrary $n$-body correlators, spectral functions, and thermal averages (Patu, 2022, Hao et al., 2022);
• Tensor-network methods (TEBD) for real-time evolution under spatially varying interactions in transmission lines (Fatis et al., 12 Jan 2026).

The mapping to free fermions ensures that computational complexity is controlled, and enables direct evaluation of experimentally relevant observables, contingent only on explicit diagonalization of the effective single-particle Hamiltonian and evaluation of well-characterized determinant structures.

7. Phase Diagram and Regimes of Validity

In the JCH model and equivalent platforms, the TG phase occupies a broad region in the phase diagram between the vacuum ($n = 0$) and the one-photon Mott lobe. In this window, for $\kappa/g \lesssim 0.13$ and density $\rho \lesssim 0.5$, the following TG features are simultaneously satisfied: $g^{(1)}(r)\sim r^{-1/2}$, $g^{(2)}(0)\approx 0$, $n(k)\sim |k|^{-1/2}$, vanishing superfluid fraction, small condensate fraction, and ground-state energy matching free fermions (D'Souza et al., 2013).

As parameters are tuned outside this regime (e.g., increasing $\kappa/g$ or density), the system crosses over to quasi-condensed or fully superfluid photonic behavior, with correlations and momentum distributions deviating from TG predictions.

---

In summary, the Tonks–Girardeau gas of fermionized photons represents an archetypal realization of strongly correlated bosonic light in one dimension, with analytic tractability supported by the Bose–Fermi mapping, concrete proposals and observations in several experimental platforms, and rigorous predictions for correlation and spectral properties that distinguish it sharply from weakly interacting or dilute photon gases (D'Souza et al., 2013, Fatis et al., 12 Jan 2026, Sabater et al., 2024, Hao et al., 2022, Patu, 2022).