Pith. sign in

REVIEW 3 major objections 6 minor 36 references

Observation of Dynamical Fermionization

T0 review · 3 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash

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

Pith's one-line read Expanding a one-dimensional gas of strongly interacting bosons flips its momentum distribution from bosonic to fermionic, directly exposing the system's rapidities.

desk verdict First clean experimental observation of dynamical fermionization and the first measurement of a rapidity distribution, with strong no-free-parameter simulation agreement; the 'direct' claim is slightly model-mediated and the main curves lack error bars. read the letter →

arxiv 1908.05364 v1 pith:TJEU3TAO submitted 2019-08-14 cond-mat.quant-gas physics.atom-ph

Joshua M. Wilson , Neel Malvania , Yuan Le , Yicheng Zhang , Marcos Rigol , David S. Weiss This is my paper · ORCID
classification cond-mat.quant-gasphysics.atom-ph
keywords dynamicalfermionizationTonks-Girardeaugasrapidities1DBosemomentumdistributiontime-of-flightimagingintegrablequantumsystemsBose-Fermioscillations
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 reports the first direct measurement of rapidities—the conserved quantities that govern integrable many-body quantum systems—by watching the momentum distribution of a one-dimensional Tonks-Girardeau gas of strongly interacting bosons transform from a bosonic peak into a fermionic, rounded shape after the axial trap is switched off. In the Tonks-Girardeau limit, the asymptotic momentum distribution after free 1D expansion is exactly the distribution of rapidities, so a time-of-flight image taken late in the expansion amounts to a readout of those conserved quantities. The measured profiles agree almost perfectly with hard-core-boson simulations that contain no free parameters, and the same apparatus also shows the predicted oscillation between bosonic and fermionic momentum shapes after a sudden change of trap depth. The paper's contribution is to bring integrability, previously a mostly theoretical structure, into direct experimental view.

What carries the argument

The load-bearing object is the Tonks-Girardeau (T-G) gas: a 1D Bose gas whose contact interaction is so strong that no two atoms can occupy the same point. Its many-body wavefunction is the absolute value of a noninteracting spinless-fermion wavefunction, so local observables match those of fermions while the momentum distribution remains bosonic; the distribution of rapidities, however, is identical to the noninteracting Fermi momentum distribution. The experiment uses a 2D optical lattice to create many independent 1D tubes, suddenly flattens the axial potential to start 1D expansion, then shuts off the transverse lattice in 32 microseconds at a chosen time to freeze interactions and let the cloud expand for 70 milliseconds before absorption imaging. The numerical twin is a continuum-limit lattice of hard-core bosons, evolved exactly through a Jordan-Wigner mapping onto free fermions; it includes the initial tube populations, the Gaussian axial trap with anti-trap, the time-of-flight propagation, instrumental resolution, and the sum over tubes, with no adjustable parameters.

What would settle it

Take the same initial T-G gas and repeat the fermionization measurement with several lattice shutoff durations (0, 16, 32, 108, 270 microseconds as in the paper's control). If the asymptotic 1D profiles differ by more than the small broadening already reported, the claim that the late-time profile equals the momentum distribution—and hence the rapidity distribution—is falsified. A second decisive check is to compare the late-time profile directly with the computed momentum distribution of the noninteracting Fermi gas that shares the initial density; any systematic mismatch beyond experimental noise would falsify the identification.

Watch

Extended reading notes

Core claim

On its own terms, the paper establishes that a strongly interacting 1D Bose gas, when released from its axial confinement while remaining transversely confined, dynamically fermionizes: the momentum distribution, initially peaked like that of a bosonic condensate, deforms within about 12 ms into a rounded distribution that matches the momentum distribution of a noninteracting Fermi gas. Because the Tonks-Girardeau wavefunction equals the absolute value of the noninteracting Fermi wavefunction, the asymptotic momentum distribution of the expanding gas is precisely the distribution of rapidities, the conserved quantities of the integrable Lieb-Liniger model. The experiment therefore claims to have measured rapidities in an interacting many-body system for the first time, with no-free-parameter simulations reproducing the data at long expansion times. It further claims that after a sudden change in axial trap depth, the momentum distribution oscillates between fermionic and bosonic shapes with the predicted period, with small deviations from T-G theory traced to the finite interaction parameter $\gamma$.

Load-bearing premise

The measurement assumes that the 70 ms time-of-flight spatial distribution is the momentum distribution at the moment the transverse lattice is switched off, which requires the 32 microsecond lattice shutoff to remove all interaction energy before the axial wavefunction changes while keeping the flight long enough for the asymptotic distribution to form.

Editorial extensions

If this is right

  • Rapidities become experimentally accessible observables, turning integrable models from theoretical constructs into testable systems.
  • The no-free-parameter agreement validates the hard-core-boson description of 1D expansion dynamics in strongly interacting gases.
  • The observed bosonic-fermionic oscillations after trap quenches confirm the predicted coherent breathing of momentum space and locate where finite-$\gamma$ corrections become visible.
  • The same time-of-flight protocol can be applied after more general quenches and to other integrable models, such as the 1D Fermi-Hubbard model, where rapidity distributions feed into generalized hydrodynamics.
  • With rapidity distributions known, integrable-system theory can predict correlation functions and subsequent dynamics, not just the momentum profile.

Reading between the lines

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

  • If the asymptotic profile is truly the rapidity distribution, the same protocol can map how rapidity distributions evolve after arbitrary quenches, offering a direct test of generalized-hydrodynamics predictions for non-thermalizing dynamics.
  • An unstated corollary is that varying the transverse-lattice shutoff speed and extrapolating to zero shutoff time would quantify the systematic error of the sudden-freeze assumption; this is a simple extension of the control already shown in the supplement.
  • Because rapidities encode the full many-body state in integrable systems, this readout could be combined with local density or noise-correlation measurements to distinguish integrable from thermalizing dynamics in one and the same experimental run.
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, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The manuscript reports time-of-flight (TOF) measurements of a 1D Tonks-Girardeau (T-G) gas of strongly interacting ultracold bosons after the axial confining potential is suddenly removed, observing the predicted evolution of the momentum distribution from a peaked bosonic shape to a rounded fermionic shape (dynamical fermionization). The authors also study the momentum dynamics after a quench to a new finite axial trap depth, observing bosonic-fermionic oscillations. They compare their measured TOF profiles with no-free-parameter exact numerical simulations of hard-core bosons in the T-G limit, reporting good agreement, and claim that the asymptotic TOF distributions constitute the first direct measurement of the distribution of rapidities in an interacting many-body quantum system.

Significance. If the central claim holds, this is a significant result: it provides a direct experimental window into the conserved quantities (rapidities) of an integrable many-body system, a long-sought observable in quantum integrable physics. The paper's strengths are the no-free-parameter comparison between experiment and exact hard-core-boson simulations, the use of predictions made in prior independent theoretical work (Refs. [7,8]), and the clean 1D expansion protocol that isolates the momentum dynamics. The experimental observation of dynamical fermionization itself appears robust and well supported, as is the observation of the bosonic-fermionic oscillations. However, the stronger claim of 'directly measuring rapidities' is not fully supported by the data as presented, because the measured TOF profiles are taken at finite evolution times and the asymptotic limit is inferred through the simulation rather than demonstrated experimentally.

major comments (3)
  1. [Page 6, paragraph after Fig. 2D ('We have thus measured the distribution of rapidities')] The claim that the asymptotic rapidity distribution has been directly measured is not supported by the data shown in Fig. S1B, where the FWHM of the TOF profiles is still increasing at tev = 15 ms, the latest usable time before the axial potential becomes non-flat. Because the total detection time tdet is fixed and the remaining TOF is tdet - tev, a fully converged momentum distribution would produce a decreasing FWHM as tev increases; the observed increase indicates that the momentum distribution at tev is still evolving. The connection between the measured finite-time TOF profiles and the asymptotic rapidity distribution is therefore made by the no-free-parameter hard-core-boson simulation, not by a directly observed asymptotic limit. The authors should either provide an experimental demonstration that the normalized profiles have converged (for example, by showing that the profile shape is stationary after rescaling by tTOF) or explicitly state that the rapidity distribution is inferred from agreement with the T-G simulation, thereby tempering the word 'directly'.
  2. [Fig. 2D and the associated discussion of agreement] The 'almost perfect agreement' and 'essentially indistinguishable' claims are based on overlaid curves without error bars on the experimental profiles and without any quantitative goodness-of-fit metric. Since the central claim rests on this agreement, the authors should provide error bars or confidence intervals on the measured TOF profiles and a quantitative measure of the deviation between experiment and theory at each tev, especially at the longest evolution times. This would allow the reader to assess whether the agreement is indeed statistically consistent.
  3. [SI Section E (Quench from low to high ωz)] There is an internal inconsistency between the main text and the Supplemental Material regarding the breathing period of the T-G gas after the quench. The main text states 'The theoretical period is ∼9% shorter than in the experiment' and explains that finite γ makes the experimental period longer, while SI Section E states 'We observe that the experimental breathing period is 9±0.6% smaller than the T-G theory.' These two statements have opposite signs. Because this period comparison is used to support the finite-γ interpretation of the quench data, the sign error must be corrected and the claim reconciled with the data shown in Fig. 3A.
minor comments (6)
  1. [Page 5, text near Fig. 3] There is a typographical error: 'experimentaly' should be 'experimentally.'
  2. [Page 4, text near Fig. 2] The phrase 'just starts to effect the TOF distribution' uses 'effect' as a verb; it should be 'affect.'
  3. [Fig. 2D caption] The caption lists 'the first six times shown in A' but does not identify which color corresponds to which tev value; please add a legend or describe the color order.
  4. [Abstract and introduction] The statement 'Rapidities have not previously been measured in any interacting many-body quantum system' is a strong negative claim; it should be qualified with 'to the best of our knowledge' unless an exhaustive literature search is provided.
  5. [SI Section C] There is a typographical error: 'shutoff' appears to be a misspelling of 'shutoff.'
  6. [Page 4, discussion of finite-γ agreement] Because the initial dimensionless coupling is only γ ≈ 8.5 (with values as low as 4.2), the statement 'The agreement at long times suggests that the T-G gas model is sufficient for our finite γ system' would be strengthened by a brief quantitative argument or a finite-γ calculation showing that corrections are expected to be small in the asymptotic regime.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: parameter-free theory and simulations provide external benchmarks; the self-citations are not load-bearing.

full rationale

The paper's central derivation chain is an experimental measurement compared against parameter-free hard-core boson simulations and prior analytical predictions. The key theoretical input, that the asymptotic momentum distribution after 1D expansion equals the distribution of rapidities, is imported from independent published work (Refs. [6-11]) rather than derived or fitted within this paper. The numerical simulations explicitly incorporate the experimental details -- initial size, axial potential, evolution up to tev, time of flight, imaging resolution, and sum over tubes -- with no free parameters, so the agreement shown in Fig. 2D is an external falsifiable benchmark, not a fit renamed as a prediction. The presence of M. Rigol as a coauthor and the citation of his earlier papers (Refs. [7, 8, 18]) constitutes a minor self-citation, but those cited results are prior, parameter-free, and not adjusted to this dataset; Ref. [8] also has no overlapping authorship. The finite-tev concern raised by Fig. S1B -- that the FWHM is still broadening at the latest usable time -- is a legitimate experimental limitation on how directly the asymptotic distribution is observed, and the paper's reliance on a simulation of the full protocol to bridge that gap is a correctness/robustness issue rather than a circularity: the simulation is not defined in terms of the measured quantity and was not fitted to it. Therefore no specific circular reduction can be exhibited, and the appropriate finding is no significant circularity, with score 2 only to reflect the minor non-load-bearing self-citation.

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

The central claim rests on standard integrable-system mappings and on experimental approximations (flat axial potential, sudden shutoff, zero temperature, T-G limit). No parameter was fitted to the data; all simulation inputs are measured experimental values. The paper explicitly flags finite-gamma and anharmonic-trap effects as sources of deviations in the quench experiments.

assumptions (6)
  • domain assumption The 2D lattice tubes are independent 1D systems with negligible tunneling, described by the Lieb-Liniger model.
    Used throughout; if inter-tube coupling or 2D effects were significant, the T-G mapping would fail. The paper states negligible tunneling among tubes.
  • domain assumption The hard-core boson (gamma goes to infinity) model accurately represents the experimental gas with finite gamma approximately 4.2 to 8.5.
    The numerical theory assumes the T-G limit; experiment uses finite gamma. The paper argues finite-gamma corrections are small for the asymptotic distribution but sees deviations in quench dynamics.
  • domain assumption The axial potential is approximately flat over the 40 micrometer expansion region, so expansion is effectively free.
    Residual anti-trap is partially canceled by a shallow dipole trap; deviations appear past 15 ms and are excluded from analysis.
  • domain assumption The initial state in each tube is the zero-temperature ground state.
    Temperature effects are invoked to explain early-time discrepancies, so the T=0 assumption is a premise of the exact numerical match.
  • domain assumption For a T-G gas expanding in 1D, the asymptotic momentum distribution equals the rapidity distribution, which is the momentum distribution of noninteracting fermions.
    Cited Sutherland and Rigol-Muramatsu; it is the theoretical basis for identifying the measured asymptotic distribution with rapidities.
  • standard math The Jordan-Wigner/Bose-Fermi mapping exactly relates hard-core boson correlations to noninteracting fermions in 1D.
    Used for all numerical calculations in the SI.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Observation of Dynamical Fermionization." pith.science (2026). https://pith.science/paper/TJEU3TAO

@misc{pith2026190805364,
  author       = {Pith},
  title        = {Pith review of: Observation of Dynamical Fermionization},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/TJEU3TAO}},
  note         = {Machine review of arXiv:1908.05364}
}
read the original abstract

We observe dynamical fermionization, where the momentum distribution of a Tonks-Girardeau (T-G) gas of strongly interacting bosons in 1D evolves from bosonic to fermionic after its axial confinement is removed. The asymptotic momentum distribution after expansion in 1D is the distribution of rapidities, which are the conserved quantities associated with many-body integrable systems. Rapidities have not previously been measured in any interacting many-body quantum system. Our measurements agree well with T-G gas theory. We also study momentum evolution after the trap depth is suddenly changed to a new non-zero value. We observe the predicted bosonic-fermionic oscillations and see deviations from the theory outside of the T-G gas limit.

Figures

Figures reproduced from arXiv: 1908.05364 by the authors.

Figure 1
Figure 1. FIG. 1. Timing and measurement. ( [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Dynamical fermionization. ( [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Bose-Fermi oscillations (quench from low to high [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Bose-Fermi oscillations (quench from high to low [PITH_FULL_IMAGE:figures/full_fig_p010_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

36 extracted references · 24 canonical work pages

  1. [1]

    Bethe, Zeitschrift f¨ ur Physik71, 205 (1931)

    H. Bethe, Zeitschrift f¨ ur Physik71, 205 (1931)

    1931

  2. [2]

    Sutherland, Beautiful Models (World Scientific, Singapore, 2004)

    B. Sutherland, Beautiful Models (World Scientific, Singapore, 2004)

    work page 2004

  3. [3]

    M. A. Cazalilla, R. Citro, T. Giamarchi, E. Orignac, and M. Rigol, Rev. Mod. Phys. 83, 1405 (2011)

    work page 2011

  4. [4]

    X.-W. Guan, M. T. Batchelor, and C. Lee, Rev. Mod. Phys. 85, 1633 (2013)

    2013

  5. [5]

    Calabrese, F

    P. Calabrese, F. H. L. Essler, and G. Mussardo, J. Stat. Mech. 2016, 064001 (2016)

    work page 2016

  6. [6]

    Sutherland, Phys

    B. Sutherland, Phys. Rev. Lett. 80, 3678 (1998)

    work page 1998

  7. [7]

    Rigol and A

    M. Rigol and A. Muramatsu, Phys. Rev. Lett. 94, 240403 (2005)

    work page 2005

  8. [8]

    Minguzzi and D

    A. Minguzzi and D. M. Gangardt, Phys. Rev. Lett. 94, 240404 (2005)

    work page 2005

Show all 36 references
  1. [9]

    C. J. Bolech, F. Heidrich-Meisner, S. Langer, I. P. McCulloch, G. Orso, and M. Rigol, Phys. Rev. Lett. 109, 110602 (2012)

    2012

  2. [10]

    A. S. Campbell, D. M. Gangardt, and K. V. Kheruntsyan, Phys. Rev. Lett. 114, 125302 (2015)

    2015

  3. [11]

    Z. Mei, L. Vidmar, F. Heidrich-Meisner, and C. J. Bolech, Phys. Rev. A 93, 021607 (2016)

    2016

  4. [12]

    E. H. Lieb and W. Liniger, Phys. Rev. 130, 1605 (1963)

    1963

  5. [13]

    Kinoshita, T

    T. Kinoshita, T. Wenger, and D. S. Weiss, Science 305, 1125 (2004)

    2004

  6. [14]

    Paredes, A

    B. Paredes, A. Widera, V. Murg, O. Mandel, S. F¨ olling, I. Cirac, G. V. Shlyapnikov, T. W. H¨ ansch, and I. Bloch, Nature429, 277 (2004)

    2004

  7. [15]

    Kinoshita, T

    T. Kinoshita, T. Wenger, and D. S. Weiss, Phys. Rev. Lett. 95, 190406 (2005)

    2005

  8. [16]

    Bloch, J

    I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008)

    2008

  9. [17]

    Olshanii, Phys

    M. Olshanii, Phys. Rev. Lett. 81, 938 (1998)

    1998

  10. [18]

    Xu and M

    W. Xu and M. Rigol, Phys. Rev. A 92, 063623 (2015)

    2015

  11. [19]

    Rigol, Phys

    M. Rigol, Phys. Rev. A 72, 063607 (2005)

    2005

  12. [20]

    B. Fang, G. Carleo, A. Johnson, and I. Bouchoule, Phys. Rev. Lett. 113, 035301 (2014)

    2014

  13. [21]

    Menotti and S

    C. Menotti and S. Stringari, Phys. Rev. A 66, 043610 (2002)

    2002

  14. [22]

    Moritz, T

    H. Moritz, T. St¨ oferle, M. K¨ ohl, and T. Esslinger, Phys. Rev. Lett.91, 250402 (2003)

    2003

  15. [23]

    Vidmar, J

    L. Vidmar, J. P. Ronzheimer, M. Schreiber, S. Braun, S. S. Hodgman, S. Langer, F. Heidrich- Meisner, I. Bloch, and U. Schneider, Phys. Rev. Lett. 115, 175301 (2015)

    2015

  16. [24]

    Schemmer, I

    M. Schemmer, I. Bouchoule, B. Doyon, and J. Dubail, Phys. Rev. Lett. 122, 090601 (2019). 11

    2019

  17. [25]

    Bertini, M

    B. Bertini, M. Collura, J. De Nardis, and M. Fagotti, Phys. Rev. Lett. 117, 207201 (2016)

    2016

  18. [26]

    O. A. Castro-Alvaredo, B. Doyon, and T. Yoshimura, Phys. Rev. X 6, 041065 (2016)

    2016

  19. [27]

    Kinoshita, T

    T. Kinoshita, T. Wenger, and D. S. Weiss, Nature 440, 900 (2006)

    2006

  20. [28]

    Y. Tang, W. Kao, K.-Y. Li, S. Seo, K. Mallayya, M. Rigol, S. Gopalakrishnan, and B. L. Lev, Phys. Rev. X 8, 021030 (2018)

    2018

  21. [29]

    J.-S. Caux, B. Doyon, J. Dubail, R. Konik, and T. Yoshimura, SciPost Phys. 6, 70 (2019)

    2019

  22. [30]

    Kinoshita, T

    T. Kinoshita, T. Wenger, and D. S. Weiss, Physical Review A 71, 011602 (2005)

    2005

  23. [31]

    S. R. Segal, Q. Diot, E. A. Cornell, A. A. Zozulya, and D. Z. Anderson, Phys. Rev. A 81, 053601 (2010)

    2010

  24. [32]

    B. L. Tolra, K. M. O’Hara, J. H. Huckans, W. D. Phillips, S. L. Rolston, and J. V. Porto, Phys. Rev. Lett. 92, 190401 (2004)

    2004

  25. [33]

    Rigol and A

    M. Rigol and A. Muramatsu, Phys. Rev. A 70, 031603(R) (2004)

    2004

  26. [34]

    Rigol and A

    M. Rigol and A. Muramatsu, Phys. Rev. A 72, 013604 (2005)

    2005

  27. [35]

    Rigol and A

    M. Rigol and A. Muramatsu, Phys. Rev. Lett. 93, 230404 (2004)

    2004

  28. [36]

    Rigol and A

    M. Rigol and A. Muramatsu, Mod. Phys. Lett. 19, 861 (2005). ACKNOWLEDGMENTS This work was supported by the National Science Foundation under Grant No. PHY- 1707482 (Y.Z. and M.R.) and No. PHY-1707576 (D.S.W., J.M.W., N.M., and Y.L.) and by the U.S. Army Research Office No. W911N...

    2005

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.