pith. sign in

arxiv: 2506.22683 · v3 · submitted 2025-06-27 · ❄️ cond-mat.quant-gas · cond-mat.stat-mech

"Depletion" of Superfluid Density: Universal Low-temperature Thermodynamics of Superfluids

Viktor Berger , Nikolay Prokof'ev , Boris Svistunov This is my paper

Pith reviewed 2026-05-19 07:56 UTC · model grok-4.3

classification ❄️ cond-mat.quant-gas cond-mat.stat-mech
keywords superfluid density depletionlow-temperature thermodynamicsPopov hydrodynamic actionfinite-size correctionsuniversal scalingquantum superfluidslattice bosonsJ-current model
0
0 comments X

The pith

The low-temperature depletion of superfluid density in d dimensions maps onto finite-size corrections in a (d+1)-dimensional classical field theory.

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

The paper shows that the decrease in superfluid density with rising temperature at low T follows from a direct equivalence to finite-size effects in a higher-dimensional classical system. This equivalence is obtained through Popov's hydrodynamic action and holds for superfluids both with and without Galilean symmetry. It generalizes the standard Landau formula while also supplying a grand-canonical version that exposes the same power-law scaling for a wide set of thermodynamic quantities. The mapping therefore supplies a single framework that replaces separate phonon-based arguments in different symmetry classes. Numeric checks on lattice bosons and the J-current model confirm the resulting predictions.

Core claim

Using Popov's hydrodynamic action, the theory of low-temperature depletion in a d-dimensional quantum superfluid maps onto the problem of finite-size (L) corrections in a (d+1)-dimensional anisotropic (pseudo-)classical-field system with U(1)-symmetric complex-valued action. In addition to generalizing Landau's formula, the work develops the grand canonical theory, which reveals a universal scaling, T^{d+1} and 1/L^{d+1}, for finite-T and finite-L effects of many thermodynamic quantities. The results are validated with numeric simulations of interacting lattice bosons and the J-current model.

What carries the argument

The mapping, via Popov's hydrodynamic action, of low-temperature quantum depletion in d dimensions to finite-size corrections in a (d+1)-dimensional anisotropic U(1)-symmetric classical field theory.

If this is right

  • Depletion of superfluid density occurs even when phonon wind is absent.
  • Many thermodynamic quantities acquire universal T^{d+1} scaling at low temperature.
  • The same quantities acquire 1/L^{d+1} scaling in the corresponding finite-size classical problem.
  • The grand-canonical formulation extends the results beyond the usual canonical ensemble.
  • The equivalence applies to both Galilean and non-Galilean quantum superfluids.

Where Pith is reading between the lines

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

  • The same hydrodynamic-to-classical mapping may be used to import efficient Monte Carlo techniques from classical finite-size studies into quantum superfluid calculations.
  • Analogous reductions could connect other low-temperature quantum corrections to classical finite-size problems in neighboring condensed-matter settings.
  • The universal power-law scaling supplies a concrete target for precision measurements in ultracold-atom experiments across different dimensions.

Load-bearing premise

Popov's hydrodynamic action remains accurate for the superfluid at low but finite temperatures even when Galilean symmetry is absent.

What would settle it

High-precision simulations of a concrete non-Galilean superfluid model that yield a depletion exponent other than d+1 would falsify the mapping.

Figures

Figures reproduced from arXiv: 2506.22683 by Boris Svistunov, Nikolay Prokof'ev, Viktor Berger.

Figure 2
Figure 2. Figure 2: FIG. 2. Results for [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 1
Figure 1. Figure 1: FIG. 1. Extracting parameter [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Superfluid stiffness Λ [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Superfluid stiffness Λ [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Superfluid stiffness as a function of density for soft [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Superfluid stiffness as a function of density for a [PITH_FULL_IMAGE:figures/full_fig_p014_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Finite-temperature depletion of superfluid stiffness in [PITH_FULL_IMAGE:figures/full_fig_p015_9.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Winding number estimator (112) for [PITH_FULL_IMAGE:figures/full_fig_p015_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Period of torsional oscillator as a function of [PITH_FULL_IMAGE:figures/full_fig_p016_12.png] view at source ↗
read the original abstract

In a Galilean superfluid, the depletion of superfluid density with rising temperature can be attributed to thermally excited non-interacting phonons. For systems without Galilean symmetry, it has been shown [1] that ``phonon wind" is no longer responsible for the depletion of superfluid density. In this work, we develop the theory of superfluid density at low temperature ($T$) and provide detailed derivations of all results announced in [1]. Using Popov's hydrodynamic action, we show that the theory of low-temperature depletion in a $d$-dimensional quantum superfluid maps onto the problem of finite-size ($L$) corrections in a $(d+1)$-dimensional anisotropic (pseudo-)classical-field system with U(1)-symmetric complex-valued action. In addition to generalizing Landau's (canonical) formula, we develop the grand canonical theory, which in a broader context reveals a universal scaling, $T^{d+1}$ and $1/L^{d+1}$, for finite-$T$ and finite-$L$ effects of many thermodynamic quantities. We validate our theory with numeric simulations of interacting lattice bosons and the J-current model.

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

1 major / 1 minor

Summary. The paper develops a theory of low-temperature superfluid density depletion in d-dimensional quantum superfluids lacking Galilean invariance. Using Popov's hydrodynamic action for phase and density fluctuations, it maps the problem onto finite-size corrections in a (d+1)-dimensional anisotropic classical complex-field theory with U(1) symmetry. This yields a generalization of Landau's formula, a grand-canonical formulation, and a universal T^{d+1} (or 1/L^{d+1}) scaling for several thermodynamic quantities, with numerical support from interacting lattice bosons and the J-current model.

Significance. If the central mapping holds, the work supplies a parameter-free derivation of universal low-temperature scaling for superfluid stiffness and related quantities that applies beyond Galilean systems. The explicit numerical validation on lattice models and the J-current model provides concrete, falsifiable support for the predicted T^{d+1} depletion and strengthens the claim of broad applicability to quantum many-body systems.

major comments (1)
  1. [hydrodynamic mapping and derivation of the T^{d+1} scaling] The derivation of the hydrodynamic mapping (substitution of Popov's action into the partition function and subsequent integration) does not explicitly demonstrate that higher-order derivative operators generated by the underlying lattice Hamiltonian remain irrelevant at the T^{d+1} order for the superfluid stiffness in non-Galilean models. Because the paper validates the result on lattice realizations, this omission is load-bearing for the claim that the mapping captures the leading depletion without additional lattice corrections.
minor comments (1)
  1. [Abstract] The abstract states agreement with lattice simulations but does not specify the system sizes, fitting ranges, or statistical uncertainties used to extract the T^{d+1} coefficient.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their positive assessment of our work and for the constructive major comment. We address it point by point below and indicate the revisions we will make to strengthen the manuscript.

read point-by-point responses

  1. Referee: The derivation of the hydrodynamic mapping (substitution of Popov's action into the partition function and subsequent integration) does not explicitly demonstrate that higher-order derivative operators generated by the underlying lattice Hamiltonian remain irrelevant at the T^{d+1} order for the superfluid stiffness in non-Galilean models. Because the paper validates the result on lattice realizations, this omission is load-bearing for the claim that the mapping captures the leading depletion without additional lattice corrections.

    Authors: We agree that an explicit discussion of the irrelevance of higher-order derivative operators would improve the clarity of the hydrodynamic mapping. In the revised manuscript we will add a dedicated paragraph (in the section deriving the mapping from Popov's action) that analyzes the scaling dimensions of possible lattice-generated operators with four or more gradients. These operators are irrelevant at the Gaussian fixed point governing the long-wavelength phase fluctuations and enter the superfluid stiffness only at O(T^{d+2}) and higher. This argument is independent of Galilean invariance and follows directly from the gradient expansion underlying the hydrodynamic action. The numerical agreement with lattice models then serves as a consistency check rather than the sole justification. We believe this addition removes the load-bearing character of the omission while preserving the original derivation. revision: yes

Circularity Check

0 steps flagged

Minor self-citation to prior announcement; central mapping from Popov action remains independent

full rationale

The derivation chain begins from Popov's established hydrodynamic action (phase and density fluctuations with quadratic gradients) and maps the d-dimensional quantum depletion problem onto (d+1)-dimensional classical finite-size corrections. This step is presented as a direct substitution into the partition function and does not reduce the target observable to a fitted parameter or to a self-defined quantity. The citation [1] announces results whose detailed derivations are supplied here; it is not invoked as a uniqueness theorem or load-bearing premise that forbids alternatives. Validation against lattice-boson and J-current simulations supplies an external benchmark. No self-definitional, fitted-input, or ansatz-smuggling reductions are exhibited by the paper's own equations.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the applicability of Popov's hydrodynamic action at low but finite temperature and on the validity of the U(1)-symmetric complex-field representation for the classical analog.

axioms (1)
  • domain assumption Popov's hydrodynamic action accurately describes the low-temperature superfluid even without Galilean symmetry.
    Invoked to establish the mapping between the quantum depletion problem and the classical finite-size problem.

pith-pipeline@v0.9.0 · 5744 in / 1391 out tokens · 34026 ms · 2026-05-19T07:56:49.874602+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.

Reference graph

Works this paper leans on

33 extracted references · 33 canonical work pages

  1. [1]

    terms do not con- tribute to the second-order derivative (9). Linear in k0 term does produce a non-zero contribution to the free en- ergy but this contribution is proportional to the second power of k0 and thus is subleading (as a 1/L2 correction) to the contribution from quadratic in k0 terms because it contains two extra derivatives. From this point the...

    work page

  2. [2]

    Calculations are readily performed in the Fourier representation (the mode with zero wavevector is gauged out): φ(r) = λdc−1 Ldc X k̸=0 φk eik·r , (18) k ≡ kn = 2πn1 L , 2πn2 L⊥ ,

    General case Being interested in the leading contributions, we cal- culate the averages based on the first (leading) term of the action (6). Calculations are readily performed in the Fourier representation (the mode with zero wavevector is gauged out): φ(r) = λdc−1 Ldc X k̸=0 φk eik·r , (18) k ≡ kn = 2πn1 L , 2πn2 L⊥ , . . . , 2πndc L⊥ . (19) Here n = (n1...

    work page

  3. [3]

    Here we have ∆ B∥ = ∆B⊥, and from (15) and (23) we immediately see that S∥(λ = 1) = − 1 dc (30) and ∆Λs(L) = −(1 + 2/dc)θ ΛsLdc

    Hypercubic sample (λ = 1) This case is particularly simple. Here we have ∆ B∥ = ∆B⊥, and from (15) and (23) we immediately see that S∥(λ = 1) = − 1 dc (30) and ∆Λs(L) = −(1 + 2/dc)θ ΛsLdc . (31)

    work page

  4. [4]

    We will not mention λ = 0 as an argument of functions throughout this section

    The λ = 0 limit This case is of prime interest. We will not mention λ = 0 as an argument of functions throughout this section. Now the sum over n2, n 3, . . . , n dc in Eq. (26) is replaced with the corresponding (dc − 1)-dimensional integral and the remaining sum over n1 can be performed analytically. In Sec. A, we consider this and similar sums that eme...

    work page

  5. [5]

    depletion

    Sign of “depletion” Since all corrections are proportional to θ their sign is controlled by the sign of θ, which, as we will see later, can flip within one and the same model as a func- tion of control parameter(s). The other two circum- stances controlling the sign are: (i) the aspect ratio λ and (ii) whether we are considering ∆Λ (∥) s or ∆Λ (⊥) s . The...

    work page

  6. [6]

    (36) Here b labels bonds of the dc-dimensional hypercu- bic lattice, and Jb is the bond current that takes on three integer values

    The model The simplest model to simulate is the J-current model HJ /T = 1 2K X b X Jb |Jb| , J b = −1, 0, +1 . (36) Here b labels bonds of the dc-dimensional hypercu- bic lattice, and Jb is the bond current that takes on three integer values. Allowed values of the bond currents are also required to obey the zero divergence constraint on each site. The mod...

    work page

  7. [7]

    To relate this expression to the statistics of winding num- bers (cf

    Winding number estimator for θ Using F = − 1 V ln Z (T = 1) , (37) we have θ = 2 V " 1 Z ∂Z ∂(k2 0) 2 − 1 Z ∂2Z ∂(k2 0)2 # k0=0 . To relate this expression to the statistics of winding num- bers (cf. [11]), we rewrite it in terms of the phase twist φ0 in the direction x [k2 0 = (φ0/Lx)2]: θ = 2L4 x V " 1 Z ∂Z ∂(φ2 0) 2 − 1 Z ∂2Z ∂(φ2 0)2 # φ0=0 . Expressi...

    work page

  8. [8]

    The K = 0.4 data are presented in Figs

    Numeric results Simulations of the model (36) were performed for two values of K corresponding (respectively) to negative and positive sign of θ: K = 0.4 and K = 1.0. The K = 0.4 data are presented in Figs. 1–2. Figure 1 shows results for θ as a function of system size when it is computed using Eq. (39). From this plot we deduce that θ = −0.145(1) with la...

    work page

  9. [9]

    read off the page

    is Lφ = L(0) φ + L(1) φ , (54) L(0) φ = − κ 2 ˙φ2 − n(0) s 2 (∇φ)2 , (55) L(1) φ = i (νk0 − α) κ ˙φφx − κνk0αφ2 x . (56) This brings us to the following expression ∆n(ν) s = ν2 κ2 Z dτ ddr K(r, τ) − κ ⟨φ2 x⟩ , (57) K(r, τ) = ⟨ ˙φ(r, τ) φx(r, τ) ˙φ(0, 0) φx(0, 0) ⟩0 , (58) where the average is taken with respect to the action S0 = RR dτ ddrL(0) φ . Evaluat...

    work page

  10. [10]

    pres- sure

    and higher have been omitted be- cause ns is given only as a second derivative with respect to k0. The grand potential corresponding to the La- grangian (101) can be straightforwardly evaluated in the Fourier representation. Performing the functional inte- gral over the quadratic field φ, for the generalized “pres- sure” ˜p(µ, T, k0) we have ˜p(µ, T, k0) ...

    work page

  11. [11]

    L. Landau. Theory of the superfluidity of helium ii. Phys. Rev., 60:356–358, Aug 1941

    work page 1941

  12. [12]

    P. W. Anderson. Considerations on the flow of superfluid helium. Rev. Mod. Phys., 38:298–310, Apr 1966

    work page 1966

  13. [13]

    L. Onsager. Statistical hydrodynamics. Il Nuovo Ci- mento (1943-1954), 6(2):279–287, Mar 1949

    work page 1943

  14. [14]

    R.P. Feynman. Chapter ii application of quantum me- chanics to liquid helium. volume 1 of Progress in Low Temperature Physics, pages 17–53. Elsevier, 1955

    work page 1955

  15. [15]

    E. L. Andronikashvili and Yu. G. Mamaladze. Quan- tization of macroscopic motions and hydrodynamics of rotating helium ii. Rev. Mod. Phys. , 38:567–625, Oct 1966

    work page 1966

  16. [16]

    Josephson

    B.D. Josephson. Possible new effects in superconductive tunnelling. Physics Letters, 1(7):251–253, 1962

    work page 1962

  17. [17]

    J. B. Mehl and W. Zimmermann. Flow of superfluid helium in a porous medium. Phys. Rev. , 167:214–229, Mar 1968

    work page 1968

  18. [18]

    Boninsegni

    M. Boninsegni. Thin helium film on a glass substrate. Journal of Low Temperature Physics , 159(3):441–451, May 2010

    work page 2010

  19. [19]

    Khalatnikov

    I.M. Khalatnikov. An Introduction To The Theory Of Superfluidity. CRC Press, 2018

    work page 2018

  20. [20]

    Prokof’ev and Boris V

    Nikolay V. Prokof’ev and Boris V. Svistunov. Two defini- tions of superfluid density.Phys. Rev. B, 61:11282–11284, May 2000

    work page 2000

  21. [21]

    E. L. Pollock and D. M. Ceperley. Path-integral computa- tion of superfluid densities. Phys. Rev. B , 36:8343–8352, Dec 1987

    work page 1987

  22. [22]

    Worm algo- rithms for classical statistical models

    Nikolay Prokof’ev and Boris Svistunov. Worm algo- rithms for classical statistical models. Phys. Rev. Lett. , 87:160601, Sep 2001

    work page 2001

  23. [23]

    Svistunov, E

    B. Svistunov, E. Babaev, and N. Prokof’ev. Superfluid States of Matter . Taylor & Francis, 2015

    work page 2015

  24. [24]

    V.N. Popov. Functional Integrals in Quantum Field The- ory and Statistical Physics . Springer Dordrecht, 1983

    work page 1983

  25. [25]

    N. V. Prokof’ev, B. V. Svistunov, and I. S. Tupitsyn. Ex- act, complete, and universal continuous-time worldline monte carlo approach to the statistics of discrete quan- tum systems. Journal of Experimental and Theoretical Physics, 87(2):310–321, Aug 1998

    work page 1998

  26. [26]

    Boninsegni, N

    M. Boninsegni, N. V. Prokof’ev, and B. V. Svistunov. Worm algorithm and diagrammatic monte carlo: A new approach to continuous-space path integral monte carlo simulations. Phys. Rev. E , 74:036701, Sep 2006

    work page 2006

  27. [27]

    F. D. M. Pobell, H. W. Chan, L. R. Corruccini, R. P. Henkel, S. W. Schwenterly, and J. D. Reppy. Low- temperature superfluid density in a restricted geometry. Phys. Rev. Lett., 28:542–544, Feb 1972

    work page 1972

  28. [28]

    C. W. Kiewiet, H. E. Hall, and J. D. Reppy. Superfluid density in porous vycor glass. Phys. Rev. Lett., 35:1286– 1289, Nov 1975

    work page 1975

  29. [29]

    Cao Lie-zhao, D. F. Brewer, C. Girit, E. N. Smith, and J. D. Reppy. Flow and torsional oscillator measurements on liquid helium in restricted geometries under pressure. Phys. Rev. B , 33:106–117, Jan 1986

    work page 1986

  30. [30]

    Halperin, Jeevak M

    William P. Halperin, Jeevak M. Parpia, and James A. Sauls. Superfluid helium-3 in confined quarters. Physics Today, 71(11):30–36, 11 2018

    work page 2018

  31. [31]

    smoking gun

    Anthony Hegg, Ruoshi Jiang, Jie Wang, Jinning Hou, Tao Zeng, Yucel Yildirim, and Wei Ku. Universal low- temperature fluctuation of unconventional superconduc- tors revealed: “smoking gun” leaves proper bosonic su- perfluidity the last theory standing. arXiv:2402.08730 [cond-mat.supr-con], 2024. 18 Appendix A: Useful mathematical relations Each expression ...

    work page arXiv 2024

  32. [32]

    (C3) The next coefficient ξ3 is then given by ξ3 = ∂ξ1 ∂µ = κ 2 ∂ns ∂n = κν 2

    = − ns 2 (C2) and ξ2 = 1 2 ∂2p0(µ, k2 0) ∂µ2 = κ 2 . (C3) The next coefficient ξ3 is then given by ξ3 = ∂ξ1 ∂µ = κ 2 ∂ns ∂n = κν 2 . (C4) Coefficient ξ4 is given by the derivative of ξ3. Using Eq. (97), we have ξ4 = 1 2 ∂ξ3 ∂µ = 1 4 ∂κ ∂µ ν + κ2 ∂ν ∂n = κ2 2 γ − κνλ 2 . (C5) Determining the final coefficient ξ5 is most conveniently done with the Jacobian ...

    work page

  33. [33]

    = D(ns, µ) D(k2 0, µ) = D(ns, µ) D(k2 0, n) D(k2 0, µ) D(k2 0, n) −1 . Observing/recalling that ∂ns(n, k2 0) ∂n k0=0 = ν , ∂ns(n, k2 0) ∂(k2 0) k0=0 = σ , ∂µ(n, k2 0) ∂(k2 0) k0=0 = ∂2E(n, k2 0) ∂n ∂(k2 0) k0=0 = 1 2 ∂ns(n, k2 0) ∂n k0=0 = ν 2 , ∂µ(n, k2 0) ∂n k0=0 = κ−1 , D(k2 0, µ) D(k2 0, n) −1 = ∂n(µ, k2 0) ∂µ k0=0 = κ , we conclude that ξ5 = −1 4(σ −...

    work page