Imagine you're watching a crowded stadium empty through a single gate. You can predict the overall flow rate with simple fluid dynamics, but tracking which specific person is where — that's a different, harder problem. This paper does the equivalent for quantum fermions: it shows that a semiclassical phase-space picture doesn't just get the bulk entropy right, it actually reproduces the individual eigenvalues of the microscopic correlation matrix. That's the full fingerprint, not just the summary statistic. The committed claim: for a one-dimensional Fermi gas at nonzero temperature expanding freely from a finite box into a semi-infinite vacuum, a hydrodynamic (Euler-scale) description recovers the full set of correlation matrix eigenvalues — and through them, the von Neumann entanglement entropy via the Yang-Yang formula. Prior work (the same group's 2023 result) showed this only at zero temperature and only for the entropy, not the eigenvalue spectrum itself. The extension to finite temperature and to the full eigenvalue distribution is the new ground here. The architecture is generalized hydrodynamics (GHD) applied to free fermions — not interacting, so the full GHD machinery simplifies to semiclassical Euler equations for the Wigner function in phase space. The key structural move is different for the system versus the reservoir. For the system (the original box), the phase-space density directly gives the eigenvalues — a clean, satisfying result. For the reservoir (the vacuum the gas expands into), you need a modification: the naive phase-space density misses quantum coherences carried by particles that entered the reservoir at different times. The authors introduce a 'modified filling function' that captures these coherences within the Yang-Yang formula, making the reservoir entropy match the microscopic calculation. This asymmetry between system and reservoir is the paper's most interesting structural insight. The integrity picture is solid for a theoretical/analytical paper. The authors benchmark their hydrodynamic predictions against exact numerical diagonalization of the correlation matrix for free fermions on a lattice, showing quantitative agreement across system sizes of up to ~200 lattice sites and various temperatures. The comparison is apples-to-apples: same model, exact numerics versus the semiclassical prediction. No free parameters are being tuned. The supplementary material includes seven figures covering different regimes. The limitation is that this is free fermions only — the integrable, non-interacting case — so the validation is internal and exact, but the scope is narrow. A bonus result falls out naturally: the eigenvalue description leads to a Levitov-Lesovik-type formula for the full counting statistics (FCS) of particle number in the system during the transient expansion regime. This connects entanglement structure to a measurable transport quantity, which is a meaningful bridge between the entanglement theory community and mesoscopic physics experimentalists. The obvious next experiment the authors did not run: interacting fermions. The entire framework leans on the free-fermion simplification where the Wigner function obeys a classical Euler equation exactly. For interacting integrable systems (XXZ spin chain, Lieb-Liniger gas), GHD is well-developed but the connection between the GHD filling function and correlation matrix eigenvalues is not established. The authors gesture at this in their conclusion. My honest read: this is (c) — they're saving it for the next paper, because the interacting case is a substantially harder problem that likely requires new conceptual ingredients beyond the wave-packet picture used here. For the broader field, this paper participates in the ongoing effort to understand when and why semiclassical / hydrodynamic descriptions capture genuinely quantum information — entanglement — rather than just thermodynamic quantities. The fact that phase-space density encodes correlation matrix eigenvalues (not just entropy) is a stronger statement than what existed before, and the system-vs-reservoir asymmetry teaches something real about where coherences live in expanding quantum gases.
