Pith. sign in

REVIEW 2 major objections 1 minor 44 references

Riemann Rarefaction Waves in a Strongly Interacting Fermi Gas

T0 review · 2 major / 1 minor · reviewed 2026-06-27 · grok-4.3

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read At unitarity a Fermi gas expands into a rarefaction wave whose shape and speed agree with the ideal Euler Riemann solution at all temperatures.

desk verdict The experiment realizes self-similar rarefaction waves matching the 1D Euler Riemann solution in a unitary Fermi gas, with deviations appearing as viscosity increases. read the letter →

arxiv 2606.06659 v1 pith:KQODGI5A submitted 2026-06-04 cond-mat.quant-gas cond-mat.stat-mechcond-mat.str-el

Eric A. Wolf , Martin Zwierlein This is my paper · ORCID
classification cond-mat.quant-gascond-mat.stat-mechcond-mat.str-el
keywords stronglyinteractingFermigasunitarityrarefactionwaveRiemannsolutionEulerequationBEC-BCScrossoverhydrodynamicsviscosity
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper examines the expansion of a homogeneous strongly interacting Fermi gas released into vacuum in a shock tube geometry. At unitarity the gas is scale invariant and nearly inviscid, so the resulting rarefaction wave dynamics become self-similar and match Riemann's solution of the Euler equation. This agreement holds for every temperature probed. Away from unitarity in the BEC-BCS crossover, deviations from the ideal solution grow as viscosity increases, yet approximate self-similarity persists even when sound diffusivity rises twentyfold. The work shows these gases can serve as a controllable setting for nonlinear hydrodynamics.

What carries the argument

Riemann's solution of the Euler equation applied to the rarefaction wave dynamics in the shock tube geometry

What would settle it

A measurement of the density profile or velocity field at unitarity that deviates systematically from the predicted self-similar form after accounting for imaging resolution and trap effects would falsify the agreement.

Watch

Extended reading notes

Core claim

In a shock tube geometry a unitary Fermi gas released into vacuum forms a rarefaction wave that is self-similar and agrees with the solution of the one-dimensional Euler equations for an ideal fluid, for every temperature examined. Deviations from this ideal behavior increase as the interaction strength moves away from unitarity and viscosity rises, although self-similarity is still roughly preserved on the BCS side.

Load-bearing premise

The observed expansion can be directly compared to the ideal 1D Euler Riemann solution without significant corrections from the trapping potential, finite imaging resolution, or three-dimensional effects.

Editorial extensions

If this is right

  • The rarefaction dynamics remain self-similar at unitarity independent of temperature.
  • Deviations from the Riemann solution grow with increasing viscosity away from unitarity.
  • Approximate self-similarity persists on the BCS side even when sound diffusivity increases twentyfold.
  • Strongly interacting Fermi gases provide a controllable platform for studying nonlinear hydrodynamics.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The setup could be extended to other initial conditions to test whether the same self-similar match appears for shock waves.
  • If three-dimensional effects remain negligible, the same gas could be used to study hydrodynamic instabilities in controlled geometries.
  • The persistence of self-similarity at long times may indicate a general property of low-viscosity flows in one dimension.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 1 minor

Summary. The paper investigates rarefaction wave dynamics in a homogeneous strongly interacting Fermi gas released from a shock-tube geometry into vacuum. At unitarity the flow is reported to be self-similar and in excellent agreement with the exact Riemann solution of the 1D inviscid Euler equations (using the unitary Fermi equation of state) for all temperatures studied. Away from unitarity in the BEC-BCS crossover, deviations from the Riemann solution grow with increasing viscosity, yet approximate self-similarity is still observed even when sound diffusivity increases by a factor of twenty; this is interpreted as the long-time approach of 1D Navier-Stokes rarefaction flows to the Euler self-similar solution.

Significance. If the direct comparison to the parameter-free Riemann solution holds after systematics are quantified, the work supplies a clean, tunable experimental platform for nonlinear hydrodynamics in a scale-invariant quantum fluid. The observation that self-similarity persists even with substantially elevated viscosity supplies a concrete test of how viscous corrections decay at long times, a point of broader interest in fluid dynamics.

major comments (2)
  1. [Abstract / main results section] Abstract and main text (central claim): the assertion of 'excellent agreement' with the ideal 1D Euler Riemann solution is load-bearing for the paper's conclusion, yet no quantitative bounds are given on residual harmonic confinement, finite imaging resolution, or transverse (3D) expansion effects. These corrections could shift the effective initial discontinuity or smooth the rarefaction fan at a level comparable to the claimed agreement; explicit upper limits on each contribution (e.g., via auxiliary measurements or simulations) are required before the match can be regarded as conclusive rather than suggestive.
  2. [Results / comparison to Riemann solution] Main text (comparison procedure): the manuscript states that the observed density and velocity profiles are compared directly to the self-similar Riemann solution without fitted parameters. It is not shown how the experimental initial discontinuity is mapped onto the theoretical step, nor how the integration time is chosen relative to the sound-crossing time; any ambiguity here would undermine the parameter-free character of the test.
minor comments (1)
  1. [Figures and methods] Figure captions and text should explicitly state the imaging resolution, trap frequencies after release, and the time window over which self-similarity is assessed.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments. We agree that quantitative bounds on systematics are needed to make the agreement conclusive and will revise the manuscript accordingly. Our responses to the major comments follow.

read point-by-point responses

  1. Referee: [Abstract / main results section] Abstract and main text (central claim): the assertion of 'excellent agreement' with the ideal 1D Euler Riemann solution is load-bearing for the paper's conclusion, yet no quantitative bounds are given on residual harmonic confinement, finite imaging resolution, or transverse (3D) expansion effects. These corrections could shift the effective initial discontinuity or smooth the rarefaction fan at a level comparable to the claimed agreement; explicit upper limits on each contribution (e.g., via auxiliary measurements or simulations) are required before the match can be regarded as conclusive rather than suggestive.

    Authors: We agree that the current manuscript lacks explicit quantitative upper limits on these effects. In the revision we will add estimates derived from auxiliary trap-frequency measurements (for residual confinement), point-spread-function characterization (for imaging resolution), and 3D hydrodynamic simulations (for transverse expansion). These bounds will be shown to lie below the level of the observed agreement with the Riemann solution. revision: yes

  2. Referee: [Results / comparison to Riemann solution] Main text (comparison procedure): the manuscript states that the observed density and velocity profiles are compared directly to the self-similar Riemann solution without fitted parameters. It is not shown how the experimental initial discontinuity is mapped onto the theoretical step, nor how the integration time is chosen relative to the sound-crossing time; any ambiguity here would undermine the parameter-free character of the test.

    Authors: The initial density and velocity profiles measured at the earliest post-release time serve as the step-function initial condition, with the discontinuity location fixed by the known optical-trap geometry. Time is the laboratory expansion time normalized by the sound-crossing time t_s = L/c_s (L = initial cloud length, c_s from the unitary EOS). We will expand the methods section with an explicit description and a schematic of this mapping to remove any ambiguity. revision: yes

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity; experimental comparison to independent classical Riemann solution.

full rationale

The paper reports direct experimental observations of rarefaction waves in a unitary Fermi gas released in a shock-tube geometry and compares the measured self-similar profiles to Riemann's established analytical solution of the 1D inviscid Euler equations (using the known unitary equation of state). No derivation chain, parameter fitting, or self-citation is invoked to generate the predicted profiles; the claim is an empirical match to an external classical result. The provided abstract and context contain no equations or steps that reduce the reported agreement to inputs by construction.

Assumptions & free parameters 0 free parameters · 1 assumptions · 0 invented entities

The central claim rests on standard hydrodynamic assumptions (Euler equations apply in the inviscid limit) and the domain assumption of scale invariance at unitarity; no free parameters or new entities are introduced in the abstract.

assumptions (1)
  • domain assumption The unitary Fermi gas obeys the scale-invariant Euler equations with negligible viscosity
    Invoked when stating agreement with Riemann's solution for all probed temperatures.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Riemann Rarefaction Waves in a Strongly Interacting Fermi Gas." pith.science (2026). https://pith.science/paper/KQODGI5A

@misc{pith2026260606659,
  author       = {Pith},
  title        = {Pith review of: Riemann Rarefaction Waves in a Strongly Interacting Fermi Gas},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/KQODGI5A}},
  note         = {Machine review of arXiv:2606.06659}
}
read the original abstract

We investigate the expansion of a homogeneous, strongly interacting Fermi gas released into vacuum in a ``shock tube'' geometry. At unitarity, where the gas is scale invariant and nearly inviscid, we find that the resulting rarefaction wave dynamics are self-similar and in excellent agreement with Riemann's solution of the Euler equation for all temperatures probed. Probing interactions away from unitarity within the BEC-BCS crossover, we observe increasing deviations from the Riemann solution as viscosity increases. However, even on the BCS side, where the sound diffusivity is increased twenty-fold, self-similarity is still approximately preserved. This may reflect how 1D Navier-Stokes rarefaction flows approach Euler self-similar solutions at long times. Our work demonstrates the utility of strongly interacting Fermi gases for the study of nonlinear hydrodynamics in a highly controllable setting.

Figures

Figures reproduced from arXiv: 2606.06659 by the authors.

Figure 1
Figure 1. Experimental protocol. (a) A schematic 3D render [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Self-similar expansion. (a) A collection of normal [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. Expansion at different c0. (a) A collection of normalized densities ρ˜ vs. x, for expansion times from 1.0 to 4.0 ms. Plotted in blue is the same data as [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Expansion at different kF a values. (a) & (b): Ridgeline plots of normalized density ρ˜ vs. (a) position x and (b) ˜ξ for different initial values of the interaction parameter (kF a) −1 . Displayed data is for expansion times from t = 1.0 ms (light) to t = 4.0 ms (dark…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

44 extracted references · 2 canonical work pages

  1. [1]

    Deviations are calculated using a piecewise linear best-fit curve chosen for each dataset [37]

    (c) RMS deviations from self-similarity of the data plot- ted in (b). Deviations are calculated using a piecewise linear best-fit curve chosen for each dataset [37]. (d) A plot of the sound diffusivities measured via box shaking resonances for each of the(kF a)−1 values plotted in (a) and (b). ∝a 2 decreases away from unitarity, the mean free path l= (nσ)...

  2. [2]

    Keimer, S

    B. Keimer, S. A. Kivelson, M. R. Norman, S. Uchida, and J. Zaanen, From quantum matter to high-temperature superconductivity in copper oxides, Nature518, 179 (2015)

    2015

  3. [3]

    Schmitt and P

    A. Schmitt and P. Shternin, Reaction rates and trans- port in neutron stars, inThe Physics and Astrophysics of Neutron Stars, edited by L. Rezzolla, P. Pizzochero, D.I.Jones, N.Rea,andI.Vidaña(SpringerInternational Publishing, Cham, 2018) pp. 455–574

    2018

  4. [4]

    Thomas and T

    S. Thomas and T. Derek, Nearly perfect fluidity: from cold atomic gases to hot quark gluon plasmas, Reports on Progress in Physics72, 126001 (2009)

    2009

  5. [5]

    Shuryak, Strongly coupled quark-gluon plasma in heavy ion collisions, Rev

    E. Shuryak, Strongly coupled quark-gluon plasma in heavy ion collisions, Rev. Mod. Phys.89, 035001 (2017)

    2017

  6. [6]

    Giorgini, L

    S. Giorgini, L. P. Pitaevskii, and S. Stringari, Theory of ultracold atomic Fermi gases, Reviews of Modern Physics 80, 1215 (2008)

    2008

  7. [7]

    Ketterle and M

    W. Ketterle and M. Zwierlein, Making, probing and un- derstanding ultracold Fermi gases, Rivista del Nuovo Ci- mento31(2008)

    2008

  8. [8]

    Zwerger, ed.,The BCS-BEC crossover and the uni- tary Fermi gas, Vol

    W. Zwerger, ed.,The BCS-BEC crossover and the uni- tary Fermi gas, Vol. 836 (Springer, 2012)

    2012

Show all 44 references
  1. [9]

    Enrico Fermi

    M. W. Zwierlein, Thermodynamics of strongly inter- acting Fermi gases, inProceedings of the International School of Physics "Enrico Fermi", Vol. 191 (IOS Press,

  2. [10]

    L. P. Pitaevskii and S. Stringari,Bose-Einstein Conden- sation and Superfluidity(Oxford University Press, Ox- ford, 2016)

    2016

  3. [11]

    Krinner, T

    S. Krinner, T. Esslinger, and J.-P. Brantut, Two-terminal transport measurements with cold atoms, Journal of Physics: Condensed Matter29, 343003 (2017)

    2017

  4. [12]

    Sommer, M

    A. Sommer, M. Ku, G. Roati, and M. W. Zwierlein, Uni- versal spin transport in a strongly interacting Fermi gas, Nature472, 201 (2011)

    2011

  5. [13]

    Luciuk, S

    C. Luciuk, S. Smale, F. Böttcher, H. Sharum, B. A. Olsen, S. Trotzky, T. Enss, and J. H. Thywissen, Obser- vation of quantum-limited spin transport in strongly in- teracting two-dimensional Fermi gases, Phys. Rev. Lett. 118, 130405 (2017)

    2017

  6. [14]

    C. Cao, E. Elliott, J. Joseph, H. Wu, J. Petricka, T. Schäfer, and J. E. Thomas, Universal quantum vis- cosity in a unitary Fermi gas, Science331, 58 (2011)

    2011

  7. [15]

    P. B. Patel, Z. Yan, B. Mukherjee, R. J. Fletcher, J. Struck, and M. W. Zwierlein, Universal sound diffu- sion in a strongly interacting Fermi gas, Science370, 1222 (2020)

    2020

  8. [16]

    Huang, Y

    S. Huang, Y. Ji, T. Repplinger, G. G. T. Assumpção, J. Chen, G. L. Schumacher, F. J. Vivanco, H. Kurkjian, and N. Navon, Emergence of sound in a tunable Fermi fluid, Phys. Rev. X15, 011074 (2025)

    2025

  9. [17]

    X. Li, X. Luo, S. Wang, K. Xie, X.-P. Liu, H. Hu, Y.-A. Chen, X.-C. Yao, and J.-W. Pan, Second sound attenua- tion near quantum criticality, Science375, 528 (2022)

    2022

  10. [18]

    Z. Yan, P. B. Patel, B. Mukherjee, C. J. Vale, R. J. Fletcher, and M. W. Zwierlein, Thermography of the su- perfluid transition in a strongly interacting Fermi gas, Science383, 629 (2024)

    2024

  11. [19]

    X. Li, J. Huang, and J. E. Thomas, Universal density shift coefficients for the thermal conductivity and shear viscosity of a unitary Fermi gas, Phys. Rev. Res.6, L042021 (2024)

    2024

  12. [20]

    M. J. H. Ku, A. T. Sommer, L. W. Cheuk, and M. W. Zwierlein, Revealing the superfluid lambda transition in the universal thermodynamics of a unitary Fermi gas, Science335, 563 (2012)

    2012

  13. [21]

    C.Menotti, P.Pedri,andS.Stringari,Expansionofanin- teracting Fermi gas, Phys. Rev. Lett.89, 250402 (2002)

    2002

  14. [22]

    K. M. O’Hara, S. L. Hemmer, M. E. Gehm, S. R. Granade, and J. E. Thomas, Observation of a strongly interacting degenerate Fermi gas of atoms, Science298, 2179 (2002)

    2002

  15. [23]

    Schäfer and C

    T. Schäfer and C. Chafin, Scaling flows and dissipation in the dilute Fermi gas at unitarity, inThe BCS-BEC Crossover and the Unitary Fermi Gas, edited by W. Zw- erger (Springer Berlin Heidelberg, Berlin, Heidelberg,

  16. [24]

    M. K. Tey, L. A. Sidorenkov, E. R. S. Guajardo, R. Grimm, M. J. H. Ku, M. W. Zwierlein, Y.-H. Hou, L. Pitaevskii, and S. Stringari, Collective modes in a unitary Fermi gas across the superfluid phase transition, Phys. Rev. Lett.110, 055303 (2013)

    2013

  17. [25]

    Y.-H. Hou, L. P. Pitaevskii, and S. Stringari, Scaling solutions of the two-fluid hydrodynamic equations in a harmonically trapped gas at unitarity, Phys. Rev. A87, 033620 (2013)

    2013

  18. [26]

    E. A. Kuznetsov, M. Y. Kagan, and A. V. Turlapov, Ex- pansion of the strongly interacting superfluid Fermi gas: Symmetries and self-similar regimes, Phys. Rev. A101, 043612 (2020)

    2020

  19. [27]

    Hou and T

    J. Hou and T. Schäfer, Dissipative superfluid hydrody- namics for the unitary Fermi gas, Phys. Rev. A104, 023313 (2021)

    2021

  20. [28]

    J. A. Joseph, J. E. Thomas, M. Kulkarni, and A. G. Abanov, Observation of shock waves in a strongly inter- acting Fermi gas, Phys. Rev. Lett.106, 150401 (2011)

    2011

  21. [29]

    J. A. Joseph, E. Elliott, and J. E. Thomas, Shear vis- cosity of a unitary Fermi gas near the superfluid phase transition, Phys. Rev. Lett.115, 020401 (2015)

    2015

  22. [30]

    Elliott, J

    E. Elliott, J. A. Joseph, and J. E. Thomas, Observa- tionofconformalsymmetrybreakingandscaleinvariance in expanding Fermi gases, Phys. Rev. Lett.112, 040405 (2014)

    2014

  23. [31]

    L. Wang, X. Yan, J. Min, D. Sun, X. Xie, S.-G. Peng, M. Zhan, and K. Jiang, Scale invariance of a spherical unitary Fermi gas, Phys. Rev. Lett.132, 243403 (2024)

    2024

  24. [32]

    Mukherjee, Z

    B. Mukherjee, Z. Yan, P. B. Patel, Z. Hadzibabic, T. Yef- sah, J. Struck, and M. W. Zwierlein, Homogeneous atomic Fermi gases, Phys. Rev. Lett.118, 123401 (2017)

    2017

  25. [33]

    Riemann, Über die Fortpflanzung ebener Luftwellen von endlicher Schwingungsweite, Abhandlungen der Königlichen Gesellschaft der Wissenschaften zu Göttin- gen8(1860)

    B. Riemann, Über die Fortpflanzung ebener Luftwellen von endlicher Schwingungsweite, Abhandlungen der Königlichen Gesellschaft der Wissenschaften zu Göttin- gen8(1860)

  26. [34]

    Bartenstein, A

    M. Bartenstein, A. Altmeyer, S. Riedl, R. Geursen, S. Jochim, C. Chin, J. H. Denschlag, R. Grimm, A. Si- moni, E. Tiesinga, C. J. Williams, and P. S. Juli- enne, Precise determination of6Licold collision param- eters by radio-frequency spectroscopy on weakly bound molecules, P...

    2005

  27. [35]

    C. H. Schunck, Y. Shin, A. Schirotzek, and W. Ketterle, 6 Determination of the fermion pair size in a resonantly interacting superfluid, Nature454, 739 (2008)

    2008

  28. [36]

    G. Zürn, T. Lompe, A. N. Wenz, S. Jochim, P. S. Julienne, and J. M. Hutson, Precise characterization of 6LiFeshbach resonances using trap-sideband-resolved rf spectroscopyofweaklyboundmolecules,Phys.Rev.Lett. 110, 135301 (2013)

    2013

  29. [37]

    L. D. Landau and E. M. Lifshitz,Fluid Mechanics, 2nd ed. (Pergamon Press, 1987)

    1987

  30. [38]

    See Supplemental Material for details on experimental methods and Navier–Stokes simulations

  31. [39]

    Y. Liao, A. Grankin, A. Panigrahi, V. Galitski, and L. Levitov, Tunable viscosity across the BCS-BEC crossover (2026), arXiv:2604.10759

    2026 arXiv

  32. [40]

    H. Hu, A. Minguzzi, X.-J. Liu, and M. P. Tosi, Collective modes and ballistic expansion of a Fermi gas in the BCS- BEC crossover, Phys. Rev. Lett.93, 190403 (2004)

    2004

  33. [41]

    M. F. Fernández, D. H. Rajkov, G. D. Pace, N. Grani, M. Inguscio, F. Scazza, S. Stringari, and G. Roati, An- gular momentum of rotating fermionic superfluids by Sagnac phonon interferometry (2025), arXiv:2511.02664

    2025

  34. [42]

    Patel,Quantum transport in strongly interacting, ul- tracold Fermi gases in box potentials, Ph.D

    P. Patel,Quantum transport in strongly interacting, ul- tracold Fermi gases in box potentials, Ph.D. thesis, Mas- sachusetts Institute of Technology (2022)

    2022

  35. [43]

    Ketterle and M

    W. Ketterle and M. W. Zwierlein, Making, probing and understanding ultracold Fermi gases, La Rivista del Nuovo Cimento31, 247–422 (2008)

    2008

  36. [44]

    X. Li, J. Huang, and J. E. Thomas, Universal density shift coefficients for the thermal conductivity and shear viscosity of a unitary Fermi gas, Phys. Rev. Res.6, L042021 (2024). 7 SUPPLEMENT AR Y INFORMA TION A. Polarization rotation imaging Due to the high optical density of...

    2024

Pith tools

Reviewed June 27, 2026 · model on record in the stance chip above.