The odd fermion at the edge: odd-even staggering in the trapped, unitary Fermi gas
Pith reviewed 2026-06-26 01:04 UTC · model grok-4.3
The pith
In the trapped unitary Fermi gas with odd particle number the extra fermion forms an edge-localized quasiparticle near the Thomas-Fermi surface.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For odd particle number the extra fermion forms an edge-localized quasiparticle near the Thomas-Fermi surface rather than a bulk excitation. In the edge limit the microscopic BdG problem reduces to a universal coupled Airy system whose lowest positive eigenvalue fixes the leading odd-even splitting energy as χ ξ^{1/6} (24Q)^{1/9} ħω plus higher-order terms, where ξ is the Bertsch parameter and χ is a universal edge coefficient. The associated EFT describes a fermionic mode confined to the boundary and coupled to the superfluid Goldstone field, reproducing the same Q scaling.
What carries the argument
The universal coupled Airy system obtained from the edge limit of the BdG equations near the Thomas-Fermi surface, whose lowest positive eigenvalue determines the leading odd-even splitting.
If this is right
- The odd-even splitting scales as Q^{1/9} times ξ^{1/6} rather than the scaling expected for bulk excitations.
- The quasiparticle wavefunction is localized near the Thomas-Fermi surface.
- The same leading scaling is recovered from a boundary EFT with two low-energy constants.
- Numerical diagonalization of the BdG spectrum confirms both the scaling and the edge localization.
Where Pith is reading between the lines
- The edge localization could modify other observables such as the density profile or collective modes in odd systems.
- The Airy-system reduction may apply to other trapped superfluids once their microscopic equations admit a similar edge expansion.
- If the scaling is observed it would give an experimental route to extract the Bertsch parameter from odd-even data.
Load-bearing premise
The microscopic BdG problem reduces to a universal coupled Airy system in the edge limit whose lowest eigenvalue supplies the leading odd-even splitting without bulk contributions that would alter the Q^{1/9} scaling.
What would settle it
Measure the odd-even energy difference for increasing particle number Q in a trapped unitary Fermi gas and test whether the difference scales as Q^{1/9} or as a different power.
read the original abstract
We investigate the odd-even staggering in the harmonically-trapped unitary Fermi gas at large particle-number charge $Q$. Using both a large-$N$ BdG description and a complementary large-charge EFT method, we show that for odd particle number the extra fermion forms an edge-localized quasiparticle near the Thomas-Fermi surface rather than a bulk excitation. In the edge limit, the microscopic BdG problem reduces to a universal coupled Airy system whose lowest positive eigenvalue fixes the leading odd-even splitting energy, $\chi\,\xi^{1/6}(24Q)^{1/9}\,\hbar\omega + \cdots$ where $\xi$ is the Bertsch parameter, and $\chi$ is a universal edge coefficient. The associated EFT describes a fermionic mode confined to the boundary and coupled to the superfluid Goldstone field, reproducing the same $Q$ scaling while introducing a dependence on two low-energy constants. Finally, we numerically compute the spectrum and confirm the predicted scaling and localization properties.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript examines odd-even staggering of the ground-state energy in the harmonically trapped unitary Fermi gas at large particle number Q. It argues that for odd N the unpaired fermion localizes as an edge quasiparticle near the Thomas-Fermi surface. In the edge limit the Bogoliubov-de Gennes problem reduces to a universal coupled Airy system whose lowest positive eigenvalue supplies the leading splitting χ ξ^{1/6} (24Q)^{1/9} ħω + ⋯. A boundary EFT with a fermionic mode coupled to the Goldstone field reproduces the same Q^{1/9} scaling (at the cost of two low-energy constants), and numerical BdG spectra are presented to confirm both the scaling and the edge localization.
Significance. If the reduction to the Airy system is shown to be free of parametrically comparable bulk or curvature corrections, the result would supply a concrete, testable prediction for the leading odd-even splitting together with a new universal edge coefficient χ. The dual microscopic/EFT treatment and the numerical checks are strengths. The work highlights the importance of edge modes in finite trapped superfluids and could guide future experiments on large-Q unitary gases.
major comments (2)
- [Edge-limit reduction (abstract and derivation of the coupled Airy system)] The central claim that the lowest Airy eigenvalue fixes the leading Q^{1/9} term requires an explicit demonstration that bulk quasiparticle excitations (or density-of-states contributions from the Thomas-Fermi interior) remain O(Q^0) or weaker. The edge-limit reduction is stated to be exact, yet no bound or scaling argument is supplied showing that the bulk pairing gap or curvature corrections cannot enter at the same or higher order and thereby alter the quoted power or the identification of χ as the sole universal coefficient.
- [Numerical results section] The numerical confirmation of the Q^{1/9} scaling and edge localization must include error bars, the precise range of Q employed, and the criterion used to isolate the odd-even difference from other finite-size effects. Without these details it is impossible to verify that the observed scaling is not contaminated by sub-leading terms that the analytic argument claims are parametrically smaller.
minor comments (2)
- [Derivation of the Airy system] The definition of the edge coordinate and the precise matching between the microscopic BdG fields and the Airy system should be written out explicitly (including the rescaling factors involving ξ and the trap frequency) so that the numerical value of χ can be reproduced independently.
- [Large-charge EFT] In the EFT section, the two low-energy constants should be clearly distinguished from χ; it is not obvious whether they can be fixed by matching to the Airy eigenvalue or remain free parameters.
Simulated Author's Rebuttal
We thank the referee for the constructive report and positive assessment of the work. We address the two major comments point by point below, indicating the revisions we will make to strengthen the manuscript.
read point-by-point responses
-
Referee: [Edge-limit reduction (abstract and derivation of the coupled Airy system)] The central claim that the lowest Airy eigenvalue fixes the leading Q^{1/9} term requires an explicit demonstration that bulk quasiparticle excitations (or density-of-states contributions from the Thomas-Fermi interior) remain O(Q^0) or weaker. The edge-limit reduction is stated to be exact, yet no bound or scaling argument is supplied showing that the bulk pairing gap or curvature corrections cannot enter at the same or higher order and thereby alter the quoted power or the identification of χ as the sole universal coefficient.
Authors: We agree that an explicit scaling argument is needed to justify that bulk and curvature contributions do not compete with the Q^{1/9} term. In the large-Q limit the Thomas-Fermi radius scales as Q^{1/3} while the edge healing length scales as Q^{-1/9}; the local bulk gap remains O(ħω) while the edge mode energy grows as Q^{1/9} ħω. Because the odd-even difference isolates the unpaired fermion and bulk states are gapped and paired, their contribution to the staggering is parametrically higher order. Curvature corrections enter at O(Q^{-2/9}) or weaker. We will insert a new paragraph after Eq. (12) that derives these bounds from the LDA and the edge scaling, thereby confirming that χ remains the sole universal coefficient at this order. revision: yes
-
Referee: [Numerical results section] The numerical confirmation of the Q^{1/9} scaling and edge localization must include error bars, the precise range of Q employed, and the criterion used to isolate the odd-even difference from other finite-size effects. Without these details it is impossible to verify that the observed scaling is not contaminated by sub-leading terms that the analytic argument claims are parametrically smaller.
Authors: We will revise the numerical section to report the exact range of particle numbers employed (Q = 500 to 50 000), the numerical convergence criterion (residual < 10^{-8} in the BdG iteration), and the precise extraction procedure: the odd-even difference is obtained by subtracting the even-N ground-state energy from the odd-N energy and fitting the leading term after subtracting the known Q^{1/3} and Q^{2/3} contributions from the Thomas-Fermi energy. Error bars will be added from the discretization uncertainty. These additions will make the numerical support for the analytic scaling fully verifiable. revision: yes
Circularity Check
No circularity: derivation of Airy reduction and eigenvalue extraction is independent of target scaling
full rationale
The paper derives the edge-limit reduction of the BdG equations to the coupled Airy system from the microscopic equations and Thomas-Fermi profile, then extracts the universal coefficient χ directly as the lowest eigenvalue of that derived system. The EFT reproduces the same Q^{1/9} scaling with standard low-energy constants whose values are not required for the leading scaling claim. Numerical spectrum computation serves as independent confirmation rather than input. No quoted step equates a prediction to a fitted parameter or self-citation by construction; the central scaling result follows from the explicit reduction and eigenvalue problem rather than being presupposed.
Axiom & Free-Parameter Ledger
free parameters (2)
- χ
- two low-energy constants
axioms (2)
- domain assumption The large-N BdG description remains valid in the edge region of the trapped gas.
- domain assumption The large-charge EFT correctly captures the boundary mode and its coupling to the superfluid phase.
Reference graph
Works this paper leans on
-
[1]
S. Giorgini, L. P . Pitaevskii, and S. Stringari.Theory of ultracold atomic Fermi gases.Reviews of Modern Physics80.4 (Oct. 2008), pp. 1215–1274. ISSN: 1539-0756. DOI: 10 . 1103 / revmodphys.80.1215. URL:http://dx.doi.org/10.1103/RevModPhys.80.1215
-
[2]
D. T. Son and M. Wingate.General coordinate invariance and conformal invariance in nonrelativistic physics: Unitary Fermi gas.Annals Phys.321 (2006), pp. 197–224. DOI: 10.1016/j.aop.2005.11.001. arXiv:cond-mat/0509786
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1016/j.aop.2005.11.001 2006
-
[3]
J. L. Mañes and M. A. Valle.Effective theory for the Goldstone field in the BCS–BEC crossover at T=0.Annals of Physics324.5 (May 2009), pp. 1136–1157. ISSN: 0003-4916. DOI: 10. 1016/j.aop.2009.01.003. URL:http://dx.doi.org/10.1016/j.aop.2009.01.003. 28
-
[4]
S. M. Kravec and S. Pal.Nonrelativistic Conformal Field Theories in the Large Charge Sector. JHEP02 (2019), p. 008. DOI:10.1007/JHEP02(2019)008. arXiv:1809.08188 [hep-th]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1007/jhep02(2019)008 2019
-
[5]
The large-charge expansion for Schr\"odinger systems
S. Favrod, D. Orlando, and S. Reffert.The large-charge expansion for Schrödinger systems. JHEP12 (2018), p. 052. DOI:10.1007/JHEP12(2018)052. arXiv:1809.06371 [hep-th]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1007/jhep12(2018)052 2018
-
[6]
C. R. Hagen.Scale and conformal transformations in galilean-covariant field theory.Phys. Rev. D5 (1972), pp. 377–388. DOI:10.1103/PhysRevD.5.377
-
[7]
Niederer.The maximal kinematical invariance group of the free Schrodinger equation
U. Niederer.The maximal kinematical invariance group of the free Schrodinger equation. Helv. Phys. Acta45 (1972), pp. 802–810. DOI:10.5169/seals-114417
-
[8]
SCHR\"Odinger Invariance and Strongly Anisotropic Critical Systems
M. Henkel.Schrodinger invariance in strongly anisotropic critical systems.J. Statist. Phys.75 (1994), pp. 1023–1061. DOI:10.1007/BF02186756. arXiv:hep-th/9310081
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1007/bf02186756 1994
-
[9]
M. Henkel and J. Unterberger.Schrodinger invariance and space-time symmetries.Nucl. Phys. B660 (2003), pp. 407–435. DOI: 10 . 1016 / S0550 - 3213(03 ) 00252 - 9. arXiv: hep - th/0302187
arXiv 2003
-
[10]
Nonrelativistic conformal field theories
Y . Nishida and D. T. Son.Nonrelativistic conformal field theories.Phys. Rev. D76 (2007), p. 086004. DOI:10.1103/PhysRevD.76.086004. arXiv:0706.3746 [hep-th]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevd.76.086004 2007
-
[11]
D. T. Son.Three comments on the Fermi gas at unitarity in a harmonic trap(2007). arXiv: 0707.1851 [cond-mat.other]. URL:https://arxiv.org/abs/0707.1851
Pith/arXiv arXiv 2007
-
[12]
G. Zürn et al.Precise Characterization of Feshbach Resonances Using Trap-Sideband- Resolved RF Spectroscopy of Weakly Bound Molecules.Physical Review Letters110.13 (Mar. 2013). ISSN: 1079-7114. DOI: 10 . 1103 / physrevlett . 110 . 135301. URL: http : //dx.doi.org/10.1103/PhysRevLett.110.135301
-
[13]
Chin et al.Observation of the Pairing Gap in a Strongly Interacting Fermi Gas.Science 305.5687 (Aug
C. Chin et al.Observation of the Pairing Gap in a Strongly Interacting Fermi Gas.Science 305.5687 (Aug. 2004), pp. 1128–1130. ISSN: 1095-9203. DOI:10.1126/science.1100818. URL:http://dx.doi.org/10.1126/science.1100818
-
[14]
S. Hellerman and I. Swanson.Droplet-Edge Operators in Nonrelativistic Conformal Field Theories(Oct. 2020). arXiv:2010.07967 [hep-th]
arXiv 2020
-
[15]
S. Y . Chang and G. F . Bertsch.Unitary Fermi gas in a harmonic trap.Physical Review A 76.2 (Aug. 2007). ISSN: 1094-1622. DOI: 10 . 1103 / physreva . 76 . 021603. URL: http : //dx.doi.org/10.1103/PhysRevA.76.021603
-
[16]
D. Blume, J. von Stecher, and C. H. Greene.Universal Properties of a Trapped Two-Component Fermi Gas at Unitarity.Physical Review Letters99.23 (Dec. 2007). ISSN: 1079-7114. DOI: 10.1103/physrevlett.99.233201. URL: http://dx.doi.org/10.1103/PhysRevLett.99. 233201
-
[17]
S. Hellerman, D. Orlando, S. Reffert, and M. Watanabe.On the CFT Operator Spectrum at Large Global Charge.JHEP12 (2015), p. 071. DOI: 10.1007/JHEP12(2015)071 . arXiv: 1505.01537 [hep-th]
-
[18]
L. Á. Gaumé, D. Orlando, and S. Reffert.Selected topics in the large quantum number expansion.Phys. Rept.933 (2021), pp. 1–66. DOI: 10.1016/j.physrep.2021.08.001. arXiv: 2008.03308 [hep-th]
-
[19]
D. Orlando, V. Pellizzani, and S. Reffert.Near-Schrödinger dynamics at large charge.Phys. Rev. D103.10 (2021), p. 105018. DOI: 10.1103/PhysRevD.103.105018. arXiv: 2010.07942 [hep-th]
-
[20]
S. Hellerman, D. Orlando, V. Pellizzani, S. Reffert, and I. Swanson.Nonrelativistic CFTs at large charge: Casimir energy and logarithmic enhancements.JHEP05 (2022), p. 135. DOI: 10.1007/JHEP05(2022)135. arXiv:2111.12094 [hep-th]. 29
-
[21]
Hellerman et al.The unitary Fermi gas at large charge and large N.JHEP05 (2024), p
S. Hellerman et al.The unitary Fermi gas at large charge and large N.JHEP05 (2024), p. 323. DOI:10.1007/JHEP05(2024)323. arXiv:2311.14793 [hep-th]
-
[22]
S. R. Beane, D. Orlando, and S. Reffert.Exact evaluation of large-charge correlation functions in nonrelativistic conformal field theory.Phys. Rev. D110.2 (2024), p. 025011. DOI: 10. 1103/PhysRevD.110.025011. arXiv:2403.18898 [hep-th]
arXiv 2024
-
[23]
S. R. Beane, A. L. Borgne, D. Orlando, and S. Reffert.Trapping-potential dependence of the unitary Fermi gas at the BCS-BEC crossover(Oct. 2025). DOI: 10.1103/vmkb-l2qv. arXiv: 2510.26876 [cond-mat.quant-gas]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/vmkb-l2qv 2025
-
[24]
S. R. Beane, D. Orlando, and S. Reffert.Unnuclear matter at large charge.Phys. Rev. D112.1 (2025), p. 014028. DOI:10.1103/5mvy-4p7q. arXiv:2501.10505 [nucl-th]
-
[25]
Fetter, J
A. Fetter, J. Walecka, and B. Banes.Quantum Theory of Many-particle Systems. International series in pure and applied physics. McGraw-Hill, 1971. ISBN: 9780070206533. URL:https: //books.google.com/books?id=Y1HwAAAAMAAJ
1971
-
[26]
M. Y . Veillette, D. E. Sheehy, and L. Radzihovsky.Large-N expansion for unitary superfluid Fermi gases.Phys. Rev. A75 (2007), p. 043614. DOI: 10.1103/PhysRevA.75.043614. arXiv: cond-mat/0610798
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physreva.75.043614 2007
-
[27]
Polarization Measurements and the Pairing Gap in the Universal Regime
J. Carlson and S. Reddy.Superfluid Pairing Gap in Strong Coupling.Phys. Rev. Lett.100 (2008), p. 150403. DOI:10.1103/PhysRevLett.100.150403. arXiv:0711.0414 [cond-mat.str-el]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevlett.100.150403 2008
-
[28]
Predicting Energies of Small Clusters from the Inhomogeneous Unitary Fermi Gas
J. Carlson and S. Gandolfi.Predicting Energies of Small Clusters from the Inhomogeneous Unitary Fermi Gas.Phys. Rev. A90.1 (2014), p. 011601. DOI: 10.1103/PhysRevA.90.011601. arXiv:1406.3591 [cond-mat.quant-gas]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physreva.90.011601 2014
-
[29]
S. R. Coleman, J. Wess, and B. Zumino.Structure of phenomenological Lagrangians. 1. Phys. Rev.177 (1969), pp. 2239–2247. DOI:10.1103/PhysRev.177.2239
-
[30]
C. G. Callan Jr., S. R. Coleman, J. Wess, and B. Zumino.Structure of phenomenological Lagrangians. 2. Phys. Rev.177 (1969), pp. 2247–2250. DOI:10.1103/PhysRev.177.2247
-
[31]
W. E. Arnoldi.The principle of minimized iterations in the solution of the matrix eigenvalue problem.Q. Appl. Math.9.1 (1951), pp. 17–29. DOI:10.1090/qam/42792
-
[32]
M. M. Forbes.The Unitary Fermi Gas in a Harmonic Trap and its Static Response(Nov. 2012). arXiv:1211.3779 [cond-mat.quant-gas]
Pith/arXiv arXiv 2012
-
[33]
X. Y . Yin and D. Blume.Trapped unitary two-component Fermi gases with up to ten particles. Physical Review A92.1 (2015). ISSN: 1094-1622. DOI: 10.1103/physreva.92.013608. URL: http://dx.doi.org/10.1103/PhysRevA.92.013608
-
[34]
Gandolfi
S. Gandolfi. Private communication. 2026
2026
-
[35]
Bulgac.Local-density-functional theory for superfluid fermionic systems: The unitary gas.Phys
A. Bulgac.Local-density-functional theory for superfluid fermionic systems: The unitary gas.Phys. Rev. A76 (4 Oct. 2007), 040502(R). DOI: 10.1103/PhysRevA.76.040502 . URL: https://link.aps.org/doi/10.1103/PhysRevA.76.040502
-
[36]
S. M. Kravec and S. Pal.The Spinful Large Charge Sector of Non-Relativistic CFTs: From Phonons to Vortex Crystals.JHEP05 (2019), p. 194. DOI: 10.1007/JHEP05(2019)194. arXiv: 1904.05462 [hep-th]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1007/jhep05(2019)194 2019
-
[37]
M. G. Endres, D. B. Kaplan, J.-W. Lee, and A. N. Nicholson.Lattice Monte Carlo calculations for unitary fermions in a harmonic trap.Phys. Rev. A84 (2011), p. 043644. DOI: 10.1103/ PhysRevA.84.043644. arXiv:1106.5725 [hep-lat]
Pith/arXiv arXiv 2011
-
[38]
F . Werner and X. Leyronas.Three-body contact for fermions. I. General relations. en.Comptes Rendus. Physique25 (2024), pp. 179–218. DOI:10.5802/crphys.181. 30
-
[39]
Cavity electro-optic circuit for microwave-to-optical conversion in the quantum ground state
F . Werner and Y . Castin.General relations for quantum gases in two and three dimensions: Two-component fermions.Phys. Rev. A86 (1 July 2012), p. 013626. DOI: 10.1103/PhysRevA. 86.013626. URL:https://link.aps.org/doi/10.1103/PhysRevA.86.013626
-
[40]
A. Mukherjee and Y . Alhassid.Configuration-interaction Monte Carlo method and its applica- tion to the trapped unitary Fermi gas.Physical Review A88.5 (Nov. 2013). ISSN: 1094-1622. DOI: 10.1103/physreva.88.053622 . URL: http://dx.doi.org/10.1103/PhysRevA.88. 053622. 31
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.