pith. sign in

arxiv: 2606.26225 · v1 · pith:OD7BDWJEnew · submitted 2026-06-24 · ❄️ cond-mat.quant-gas · hep-th· nucl-th

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

classification ❄️ cond-mat.quant-gas hep-thnucl-th
keywords odd-even staggeringunitary Fermi gasedge quasiparticleBdG equationsAiry systemThomas-Fermi surfacelarge-charge EFT
0
0 comments X

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.

The paper studies odd-even energy staggering for large particle number Q in a harmonically trapped unitary Fermi gas. It shows that an odd total particle number places the unpaired fermion at the edge rather than as a bulk excitation. In the edge region the Bogoliubov-de Gennes equations reduce to a universal coupled Airy system whose lowest positive eigenvalue sets the leading splitting energy proportional to Q to the power 1/9. A boundary effective field theory that couples an edge fermion to the superfluid Goldstone mode reproduces the same scaling. Direct numerical solution of the spectrum confirms both the predicted scaling and the spatial localization of the mode.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

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)
  1. [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.
  2. [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)
  1. [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.
  2. [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

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

2 free parameters · 2 axioms · 0 invented entities

Only the abstract is available, so the ledger is necessarily incomplete. The paper introduces at least one universal coefficient χ extracted from the Airy eigenvalue and two low-energy constants in the EFT; these function as free parameters that set the numerical prefactor. The reduction to the Airy system is treated as a derived result rather than an axiom, but its validity is a load-bearing assumption.

free parameters (2)
  • χ
    Universal edge coefficient fixed by the lowest eigenvalue of the coupled Airy system; sets the overall scale of the odd-even splitting.
  • two low-energy constants
    Introduced in the EFT description of the boundary fermionic mode coupled to the Goldstone field.
axioms (2)
  • domain assumption The large-N BdG description remains valid in the edge region of the trapped gas.
    Invoked to justify the reduction to the Airy system.
  • domain assumption The large-charge EFT correctly captures the boundary mode and its coupling to the superfluid phase.
    Required for the EFT to reproduce the same Q scaling.

pith-pipeline@v0.9.1-grok · 5713 in / 1590 out tokens · 24444 ms · 2026-06-26T01:04:46.087016+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

40 extracted references · 32 canonical work pages · 10 internal anchors

  1. [1]

    Giorgini, L

    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. [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

  3. [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. [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]

  5. [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]

  6. [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. [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. [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

  9. [9]

    Henkel and J

    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

  10. [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]

  11. [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

  12. [12]

    Zürn et al.Precise Characterization of Feshbach Resonances Using Trap-Sideband- Resolved RF Spectroscopy of Weakly Bound Molecules.Physical Review Letters110.13 (Mar

    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. [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. [14]

    Hellerman and I

    S. Hellerman and I. Swanson.Droplet-Edge Operators in Nonrelativistic Conformal Field Theories(Oct. 2020). arXiv:2010.07967 [hep-th]

  15. [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. [16]

    Blume, J

    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. [17]

    Hellerman, D

    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. [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. [19]

    Orlando, V

    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. [20]

    Hellerman, D

    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. [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. [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]

  23. [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]

  24. [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. [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

  26. [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

  27. [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]

  28. [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]

  29. [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. [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. [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. [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]

  33. [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. [34]

    Gandolfi

    S. Gandolfi. Private communication. 2026

  35. [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. [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]

  37. [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]

  38. [38]

    Werner and X

    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. [39]

    Franke, J

    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. [40]

    Mukherjee and Y

    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