pith. sign in

arxiv: 2504.19310 · v1 · pith:AOWR6UA3new · submitted 2025-04-27 · ❄️ cond-mat.quant-gas · cond-mat.stat-mech

Fluctuation-induced first-order superfluid transition in unitary SU(N) Fermi gases

Georgii Kalagov This is my paper

Pith reviewed 2026-05-22 18:33 UTC · model grok-4.3

classification ❄️ cond-mat.quant-gas cond-mat.stat-mech
keywords SU(N) Fermi gassuperfluid transitionfunctional renormalization groupfirst-order phase transitionunitary regimefluctuation effectscritical temperature
0
0 comments X

The pith

Fluctuations induce a first-order superfluid transition in unitary SU(N) Fermi gases for N at least 4.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper applies the functional renormalization group to an SU(N)-symmetric Fermi gas with an auxiliary bosonic field to go beyond mean-field descriptions of pairing. It finds that fluctuations change the superfluid transition from second-order to first-order when N reaches 4 or higher. In the unitary limit this yields concrete values showing that the critical temperature falls as N grows while the jumps in superfluid gap and entropy density become larger. The calculation discretizes the effective potential on a grid and integrates the flow equations numerically to extract these thermodynamic discontinuities.

Core claim

Using the leading-order derivative expansion of the partially bosonized effective average action and grid discretization of the effective potential, the superfluid transition in unitary SU(N) Fermi gases is a fluctuation-induced first-order transition for N ≥ 4, absent at mean-field level, with the critical temperature decreasing and the discontinuities in the superfluid gap and entropy density increasing as functions of N.

What carries the argument

Leading-order derivative expansion of the partially bosonized effective average action with grid discretization of the effective potential to integrate the renormalization-group flow equations.

If this is right

  • The critical temperature in the unitary regime decreases with increasing N.
  • Discontinuities in the superfluid gap grow larger for higher N.
  • Jumps in entropy density become more pronounced, strengthening the first-order character.
  • The transition remains second-order at the mean-field level for all N.

Where Pith is reading between the lines

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

  • Similar renormalization-group treatments could be used to check whether the first-order character survives at other interaction strengths away from unitarity.
  • Higher-order terms in the derivative expansion would provide a test of how sensitive the discontinuity sizes are to the truncation.
  • The N-dependence offers a controllable knob for experiments to tune the strength of the first-order jump.

Load-bearing premise

The leading-order derivative expansion together with grid discretization of the effective potential is sufficient to determine the order of the transition and the size of the discontinuities.

What would settle it

Direct measurement of a continuous superfluid transition with no jump in the gap or entropy density for an N=4 unitary Fermi gas would contradict the first-order prediction.

Figures

Figures reproduced from arXiv: 2504.19310 by Georgii Kalagov.

Figure 1
Figure 1. Figure 1: FIG. 1. (Color online). Scale dependence of the running po [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (Color online). (a) Temperature dependence of the [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (Color online). Thermodynamic parameters of the [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (Color online). Typical temperature dependence [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
read the original abstract

We investigate the superfluid phase transition in an $\mathrm{SU}(N)$-symmetric Fermi gas with $N$ distinct spin states using the functional renormalization group. To capture pairing phenomena beyond mean-field theory, we introduce an auxiliary bosonic field and employ the leading order of the derivative expansion of the partially bosonized effective average action. By discretizing the effective potential on a grid and numerically integrating the flow equations, we resolve the thermodynamic behavior near the transition. Our results reveal a fluctuation-induced first-order phase transition for $N \geq 4$, which is absent at the mean-field level. In the unitary regime, we provide quantitative predictions for the critical temperature, as well as for the discontinuities in the superfluid gap and entropy density as functions of $N$. With increasing $N$, the critical temperature decreases, while the discontinuities become more pronounced, indicating a stronger first-order transition.

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 applies the functional renormalization group to unitary SU(N) Fermi gases in the partially bosonized formulation. Using the leading-order derivative expansion and a grid discretization of the effective potential, the authors numerically integrate the Wetterich equation and report that fluctuations drive a first-order superfluid transition for N ≥ 4, in contrast to the continuous transition found at mean-field level. Quantitative results are given for the critical temperature, the jump in the superfluid gap, and the discontinuity in entropy density as functions of N.

Significance. If the truncation is controlled, the work supplies concrete, falsifiable predictions for the order of the transition and the size of thermodynamic jumps in multi-component unitary Fermi gases, which are directly relevant to ongoing cold-atom experiments with SU(N) symmetry. The numerical FRG implementation with grid evolution of the potential is a clear technical strength that allows resolution of the barrier structure beyond analytic mean-field approximations.

major comments (2)
  1. [§3.2, Eq. (12)] §3.2 and Eq. (12): the leading-order derivative expansion (local-potential approximation) is used throughout; the barrier height that determines the first-order character is generated solely by the flow of the four-point vertex. Higher-derivative terms or momentum-dependent vertices, which are omitted, can flatten or eliminate this barrier, potentially converting the reported weak first-order jump at N=4 into a continuous transition. A quantitative estimate of the truncation error on the order parameter discontinuity is required.
  2. [§4.1, Figure 2] §4.1, Figure 2: the grid discretization of U_k(ρ) and the choice of regulator are not accompanied by convergence tests with respect to grid spacing or regulator parameter. Because the location and height of any potential barrier are sensitive to these numerical details, the reported discontinuities for N=4–6 could shift under refinement; explicit error bars or stability checks against grid size are needed to support the central claim.
minor comments (2)
  1. [Table I] The notation for the renormalized chemical potential and the definition of the unitary limit should be stated explicitly in the caption of Table I to avoid ambiguity when comparing to other FRG studies.
  2. [Figure 1] A brief discussion of how the mean-field limit is recovered by freezing the flow at the UV scale would clarify the comparison shown in Figure 1.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the careful reading of our manuscript and for highlighting both the potential impact of our results and the technical aspects of the FRG implementation. We address the two major comments point by point below. While we maintain that the leading-order calculation provides valuable new insights into fluctuation effects in SU(N) Fermi gases, we acknowledge the limitations of the truncation and will strengthen the presentation of numerical controls in the revision.

read point-by-point responses

  1. Referee: [§3.2, Eq. (12)] §3.2 and Eq. (12): the leading-order derivative expansion (local-potential approximation) is used throughout; the barrier height that determines the first-order character is generated solely by the flow of the four-point vertex. Higher-derivative terms or momentum-dependent vertices, which are omitted, can flatten or eliminate this barrier, potentially converting the reported weak first-order jump at N=4 into a continuous transition. A quantitative estimate of the truncation error on the order parameter discontinuity is required.

    Authors: We agree that the local-potential approximation constitutes a truncation and that momentum-dependent vertices or higher-derivative terms could in principle reduce the barrier height. Prior FRG work on the N=2 unitary Fermi gas indicates that the qualitative order of the transition is robust under such extensions, although quantitative shifts in the discontinuity can occur. A controlled quantitative error estimate on the jump would indeed require a next-to-leading-order derivative expansion, which leads to a substantially more involved set of flow equations and is outside the scope of the present study. In the revised manuscript we will add an explicit paragraph in Sec. 3.2 discussing this truncation uncertainty and citing relevant benchmarks from the FRG literature on fermionic systems. revision: partial

  2. Referee: [§4.1, Figure 2] §4.1, Figure 2: the grid discretization of U_k(ρ) and the choice of regulator are not accompanied by convergence tests with respect to grid spacing or regulator parameter. Because the location and height of any potential barrier are sensitive to these numerical details, the reported discontinuities for N=4–6 could shift under refinement; explicit error bars or stability checks against grid size are needed to support the central claim.

    Authors: We have performed additional numerical checks that were not reported in the original submission. Specifically, we repeated the flow for N=4 with a doubled number of grid points and with two different values of the regulator parameter; the resulting change in the order-parameter discontinuity is below 5 %. We will include these convergence tests, together with explicit error estimates, in a new subsection of Sec. 4.1 and will add error bars to the data points in Figure 2 and the tabulated values. revision: yes

standing simulated objections not resolved
  • Quantitative estimate of the truncation error on the order-parameter discontinuity without extending the derivative expansion to higher orders

Circularity Check

0 steps flagged

No significant circularity; results from numerical FRG flow integration

full rationale

The derivation proceeds by introducing an auxiliary boson, writing the leading-order derivative expansion of the partially bosonized Wetterich equation, discretizing the effective potential U_k(ρ) on a grid, and numerically integrating the flow down to k=0. The mean-field limit is recovered simply by halting the flow at the ultraviolet scale; the infrared discontinuities for N≥4 emerge from the integrated flow equations themselves rather than from any fitted parameter or self-referential definition. No load-bearing self-citation, uniqueness theorem, or ansatz smuggling is invoked to force the first-order character. The procedure is therefore self-contained against external benchmarks and does not reduce the reported critical temperature or discontinuities to the input assumptions by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The calculation rests on the validity of the leading-order derivative expansion and on the numerical stability of the grid discretization for the effective potential; no new particles or forces are introduced.

axioms (1)
  • domain assumption Leading order of the derivative expansion of the partially bosonized effective average action is sufficient near the transition.
    Stated in the abstract as the truncation employed to capture pairing beyond mean-field.

pith-pipeline@v0.9.0 · 5685 in / 1329 out tokens · 42817 ms · 2026-05-22T18:33:46.089688+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Functional renormalization group equations for antisymmetric tensor field models at finite temperature

    hep-th 2025-05 unverdicted novelty 4.0

    Derives FRG flow equations for antisymmetric tensor field models at finite temperature under SU(n) to USp(n) and SO(n) to SU(n/2) symmetry breaking.

Reference graph

Works this paper leans on

49 extracted references · 49 canonical work pages · cited by 1 Pith paper · 3 internal anchors

  1. [1]

    Scalar” modes α and σ ) and V -terms (“Vector

    Bosonic contribution The inhomogeneous part of the Hessian, obtained as the second functional derivative of Γ k with respect to the bosonic fluctuations, Eq. (19), and expanded to linear order in a, s, (sufficient for deriving the desired flows), is given by [eΓ(2) k,p(Q, Q′)](B) =β(2π)dδq0+q′ 0,0 a δ(q + q′ − p) + a∗ δ(q + q′ + p) AS 0 0 AV , (D4) 13 for...

    work page

  2. [2]

    Fermionic contribution The inhomogeneous additives to the fermionic part of the Hessian are given by [eΓ(2) k,p(Q, Q′)](F) = igk√ 2N β(2π)dδq0,q′ 0 a δ(q − q′ − p) + a∗ δ(q − q′ + p) I2 ⊗ J, (D10) for the background Eq. (34), and [eΓ(2) k,p0 (Q, Q′)](F) = gk√ 2N β(2π)dδ(q − q′) (s∗ − ia∗) δq0−q′ 0−p0,0 + (s − ia) δq0−q′ 0+p0,0 − J 0 0 0 + (s∗ + ia∗) δq0−q...

    work page

  3. [3]

    Exotic many-body physics with large-spin Fermi gases

    C. Wu, Exotic many-body physics with large-spin fermi gases, Physics 3, 92 (2010), 1103.1933

    work page internal anchor Pith review Pith/arXiv arXiv 2010

  4. [4]

    M. A. Cazalilla and A. M. Rey, Ultracold Fermi gases with emergent SU(N) symmetry, Reports on Progress in Physics 77, 124401 (2014)

    work page 2014

  5. [5]

    A. V. Gorshkov, The high-symmetry switch, Nature Physics 2014 10:10 10, 708 (2014)

    work page 2014

  6. [6]

    Mukherjee and J

    B. Mukherjee and J. M. Hutson, SU( N) symmetry with ultracold alkali dimers: Weak dependence of scattering properties on hyperfine state, Phys. Rev. Res. 7, 013099 (2025)

    work page 2025

  7. [7]

    B. J. DeSalvo, M. Yan, P. G. Mickelson, Y. N. Mar- tinez de Escobar, and T. C. Killian, Degenerate Fermi Gas of 87Sr, Phys. Rev. Lett. 105, 030402 (2010)

    work page 2010

  8. [8]

    S. Taie, R. Yamazaki, S. Sugawa, and Y. Takahashi, An SU(6) Mott insulator of an atomic Fermi gas realized by large-spin Pomeranchuk cooling, Nature Phys. 8, 825 (2012), 1208.4883

    work page internal anchor Pith review Pith/arXiv arXiv 2012

  9. [9]

    S. Taie, Y. Takasu, S. Sugawa, R. Yamazaki, T. Tsuji- moto, R. Murakami, and Y. Takahashi, Realization of a SU(2) × SU(6) system of fermions in a cold atomic gas, Phys. Rev. Lett. 105, 190401 (2010)

    work page 2010

  10. [10]

    J. S. Krauser, J. Heinze, N. Fl¨ aschner, S. G¨ otze, O. J¨ urgensen, D. S. L¨ uhmann, C. Becker, and K. Sen- gstock, Coherent multi-flavour spin dynamics in a fermionic quantum gas, Nature Physics 2012 8:11 8, 813 (2012), 1203.0948

    work page internal anchor Pith review Pith/arXiv arXiv 2012

  11. [11]

    Pagano, M

    G. Pagano, M. Mancini, G. Cappellini, P. Lombardi, F. Sch¨ afer, H. Hu, X. J. Liu, J. Catani, C. Sias, M. Ingus- cio, and L. Fallani, A one-dimensional liquid of fermions with tunable spin, Nature Physics 2014 10:3 10, 198 (2014)

    work page 2014

  12. [12]

    Zhang, M

    X. Zhang, M. Bishof, S. L. Bromley, C. V. Kraus, M. S. Safronova, P. Zoller, A. M. Rey, and J. Ye, Spectroscopic observation of SU( N)-symmetric interactions in Sr or- bital magnetism, Science 345, 1467 (2014). 15

    work page 2014

  13. [13]

    Ho and S

    T.-L. Ho and S. Yip, Pairing of Fermions with Arbitrary Spin, Phys. Rev. Lett. 82, 247 (1999)

    work page 1999

  14. [14]

    Honerkamp and W

    C. Honerkamp and W. Hofstetter, BCS pairing in Fermi systems with N different hyperfine states, Physical Re- view B - Condensed Matter and Materials Physics 70, 094521 (2004)

    work page 2004

  15. [15]

    R. W. Cherng, G. Refael, and E. Demler, Superfluidity and Magnetism in Multicomponent Ultracold Fermions, Phys. Rev. Lett. 99, 130406 (2007)

    work page 2007

  16. [16]

    Capponi, G

    S. Capponi, G. Roux, P. Lecheminant, P. Azaria, E. Boulat, and S. R. White, Molecular superfluid phase in systems of one-dimensional multicomponent fermionic cold atoms, Phys. Rev. A 77, 013624 (2008)

    work page 2008

  17. [17]

    L. He, M. Jin, and P. Zhuang, Superfluidity in a three- flavor Fermi gas with SU(3) symmetry, Phys. Rev. A 74, 033604 (2006)

    work page 2006

  18. [18]

    A. Rapp, G. Zar´ and, C. Honerkamp, and W. Hofstetter, Color Superfluidity and “Baryon” Formation in Ultracold Fermions, Phys. Rev. Lett. 98, 160405 (2007)

    work page 2007

  19. [19]

    Ozawa and G

    T. Ozawa and G. Baym, Population imbalance and pair- ing in the BCS-BEC crossover of three-component ultra- cold fermions, Phys. Rev. A 82, 063615 (2010)

    work page 2010

  20. [20]

    Schlottmann and A

    P. Schlottmann and A. A. Zvyagin, Superfluid instabil- ity in ultra-cold gas of fermionic atoms with attractive potential in a one-dimensional trap, Journal of Physics: Conference Series 529, 012013 (2014)

    work page 2014

  21. [21]

    Y. C. Zhang, Superfluid density, Josephson relation and pairing fluctuations in a multi-component fermion super- fluid, Scientific Reports 2021 11:1 11, 1 (2021)

    work page 2021

  22. [22]

    M. A. Perlin, D. Barberena, M. Mamaev, B. Sun- dar, R. J. Lewis-Swan, and A. M. Rey, Engineering infinite-range SU(N) interactions with spin-orbit-coupled fermions in an optical lattice, Phys. Rev. A 105, 023326 (2022)

    work page 2022

  23. [23]

    M. V. Komarova, M. Y. Nalimov, and J. Honkonen, Tem- perature Green’s functions in Fermi systems: The super- conducting phase transition, Theoretical and Mathemat- ical Physics 2013 176:1 176, 906 (2013)

    work page 2013

  24. [24]

    Kalagov, M

    G. Kalagov, M. Kompaniets, and M. Nalimov, Renormalization-group investigation of a superconduct- ing U( N)-phase transition using five loops calculations, Nuclear Physics B 905, 16 (2016)

    work page 2016

  25. [25]

    M. Y. Nalimov and A. V. Ovsyannikov, Convergent per- turbation theory for studying phase transitions, Theoret- ical and Mathematical Physics 204, 1033 (2020)

    work page 2020

  26. [26]

    Bednyakov and A

    A. Bednyakov and A. Pikelner, Six-loop beta functions in general scalar theory, Journal of High Energy Physics 2021 2021:4 2021, 1 (2021), 2102.12832

    work page arXiv 2021

  27. [27]

    M. J. Ku, A. T. Sommer, L. W. Cheuk, and M. W. Zwier- lein, Revealing the superfluid lambda transition in the universal thermodynamics of a unitary Fermi gas, Sci- ence 335, 563 (2012)

    work page 2012

  28. [28]

    Burovski, N

    E. Burovski, N. Prokof’ev, B. Svistunov, and M. Troyer, Critical temperature and thermodynamics of attractive fermions at unitarity, Phys. Rev. Lett.96, 160402 (2006)

    work page 2006

  29. [29]

    Bulgac, J

    A. Bulgac, J. E. Drut, and P. Magierski, Quantum Monte Carlo simulations of the BCS-BEC crossover at finite temperature, Phys. Rev. A 78, 023625 (2008)

    work page 2008

  30. [30]

    Haussmann, W

    R. Haussmann, W. Rantner, S. Cerrito, and W. Zwerger, Thermodynamics of the BCS-BEC crossover, Phys. Rev. A 75, 023610 (2007)

    work page 2007

  31. [31]

    Diehl, S

    S. Diehl, S. Floerchinger, H. Gies, J. Pawlowkski, and C. Wetterich, Functional renormalization group approach to the BCS-BEC crossover, Annalen der Physik 522, 615 (2010)

    work page 2010

  32. [32]

    Boettcher, J

    I. Boettcher, J. M. Pawlowski, and S. Diehl, Ultra- cold atoms and the Functional Renormalization Group, Nuclear Physics B - Proceedings Supplements 228, 63 (2012)

    work page 2012

  33. [33]

    Boettcher, J

    I. Boettcher, J. M. Pawlowski, and C. Wetterich, Criti- cal temperature and superfluid gap of the unitary Fermi gas from functional renormalization, Phys. Rev. A 89, 053630 (2014)

    work page 2014

  34. [34]

    Tanizaki, G

    Y. Tanizaki, G. Fej¨ os, and T. Hatsuda, Functional renor- malization group approach to conventional theory of su- perfluidity and beyond, Journal of Physics: Conference Series 529, 012027 (2014)

    work page 2014

  35. [35]

    A. N. Vasiliev, Functional Methods in Quantum Field Theory and Statistical Physics (Gordon and Breach, Am- sterdam, 1998) p. 326

    work page 1998

  36. [36]

    Zumino, Normal forms of complex matrices, Journal of Mathematical Physics 3, 1055 (1962)

    B. Zumino, Normal forms of complex matrices, Journal of Mathematical Physics 3, 1055 (1962)

    work page 1962

  37. [37]

    A. A. Abrikosov, L. P. Gorkov, and I. E. Dzialoshin- skii, Methods of quantum field theory in statistical physics (Dover Publications, 1975) p. 352

    work page 1975

  38. [38]

    Dupuis, L

    N. Dupuis, L. Canet, A. Eichhorn, W. Metzner, J. Pawlowski, M. Tissier, and N. Wschebor, The nonper- turbative functional renormalization group and its appli- cations, Physics Reports 910, 1 (2021)

    work page 2021

  39. [39]

    Berges, N

    J. Berges, N. Tetradis, and C. Wetterich, Non- perturbative renormalization flow in quantum field the- ory and statistical physics, Physics Reports 363, 223 (2002)

    work page 2002

  40. [40]

    J. M. Pawlowski, Aspects of the functional renormalisa- tion group, Annals of Physics 322, 2831 (2007)

    work page 2007

  41. [41]

    D. F. Litim, Optimisation of the exact renormalisation group, Physics Letters B 486, 92 (2000)

    work page 2000

  42. [42]

    D. F. Litim, Optimized renormalization group flows, Phys. Rev. D 64, 105007 (2001)

    work page 2001

  43. [43]

    Braaten, M

    E. Braaten, M. Kusunoki, and D. Zhang, Scattering mod- els for ultracold atoms, Annals of Physics 323, 1770 (2008)

    work page 2008

  44. [44]

    Hnatiˇ c and G

    M. Hnatiˇ c and G. Kalagov, Superfluidity in multicompo- nent fermions via the functional renormalization group, Nuclear Physics B 991, 116192 (2023)

    work page 2023

  45. [45]

    Kalagov, Supplemental material (2025), includes a Mathematica notebook (.nb) with explicit matrix expres- sions used in Appendix B

    G. Kalagov, Supplemental material (2025), includes a Mathematica notebook (.nb) with explicit matrix expres- sions used in Appendix B

    work page 2025

  46. [46]

    Grossi and N

    E. Grossi and N. Wink, Resolving phase transitions with discontinuous Galerkin methods, SciPost Phys. Core 6, 071 (2023)

    work page 2023

  47. [47]

    Koenigstein, M

    A. Koenigstein, M. J. Steil, N. Wink, E. Grossi, J. Braun, M. Buballa, and D. H. Rischke, Numerical fluid dynam- ics for FRG flow equations: Zero-dimensional QFTs as numerical test cases. I. The O( N) model, Phys. Rev. D 106, 065012 (2022)

    work page 2022

  48. [48]

    Koenigstein, M

    A. Koenigstein, M. J. Steil, N. Wink, E. Grossi, and J. Braun, Numerical fluid dynamics for FRG flow equa- tions: Zero-dimensional QFTs as numerical test cases. II. Entropy production and irreversibility of RG flows, Phys. Rev. D 106, 065013 (2022)

    work page 2022

  49. [49]

    Knoll and D

    D. Knoll and D. Keyes, Jacobian-free Newton–Krylov methods: a survey of approaches and applications, Jour- nal of Computational Physics 193, 357 (2004)

    work page 2004