{"id":"b9214606-c9b5-4686-bec4-50313819ec6b","arxiv_id":"2509.00399","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":1,"one_line_summary":"For a polytropic quantum gas, the paper claims a shock-front expansion with cloud size growing as a power of time, plus damped relaxation modes for a trapped gas; the vacuum shock result is physically problematic.","lead":"This paper uses hydrodynamic equations to describe how quantum gases released from traps expand and relax, deriving self-similar shock-like solutions in one and two dimensions. The headline vacuum-expansion result is questionable, because a cloud expanding into empty space should form a smooth rarefaction, not a shock.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The vacuum-blast self-similar solution violates energy conservation: E(t) ~ t^{-2δ/(δ+2)} for Eq. (12), so it cannot be the long-time Euler solution of an isolated gas; the shock-into-vacuum boundary condition is the fragile premise.","rationale":"The reader's weakest_assumption correctly identifies the shock-into-vacuum boundary condition as the load-bearing premise. My analysis goes further and shows that the paper's own scaling solution makes the failure explicit: total energy decays as a power law rather than remaining constant, which is an internal inconsistency with the Euler equations, not merely a disagreement with an external consensus. The artificial viscosity used in the numerics explains why simulations appear to match a dissipative scaling, but it does not validate the inviscid claim. The 2D case suffers the same energy-conservation failure, and the derivation of its jump conditions is independently questionable because ∫ r v dr is not a conserved quantity. The trapped-gas relaxation analysis (Sec. III) and finite-background results (Sec. II.B) are separable and may contain valid contributions, but the abstract and conclusions rest on the vacuum-blast shock result. Since the central claim fails, the appropriate verdict is REJECT. No judgment is made about the authors' intent; the issue is technical.","tokens_in":15603,"tokens_out":12811,"duration_ms":155034,"concrete_test":"One analytical check: insert the claimed 1D solution (Eq. (12)) into the conserved energy E = ∫ [ρv²/2 + ρ^{δ+1}/(δ+1)] dx. With ρ=t^{-b}f(ξ), v=bξ t^{b-1} (because g=bξ and b+c=1), E(t) = t^{2b-2} A + t^{-bδ} B, where A=(1/2)∫ξ²f dξ and B=(1/(δ+1))∫f^{δ+1}dξ. For b=2/(δ+2), 2b-2 = -bδ = -2δ/(δ+2), so E(t) = (A+B) t^{-2δ/(δ+2)}, not constant and vanishing at large t. This contradicts energy conservation of the Euler equations. A conservative shock-capturing simulation without artificial viscosity would show R(t) ~ t and E(t) = const, instead of the paper's exponents.","verdict_should_be":"REJECT","load_bearing_attack":"The central claim is that a finite-mass gas released into vacuum develops a shock front with scaling R(t) ~ t^b, b=2/(δ+2) in 1D and b=1/(δ+1) in 2D. This rests entirely on treating the vacuum boundary as a shock with Rankine–Hugoniot conditions (10)–(11) and (60)–(66). For a polytropic gas expanding into vacuum, the physical free-boundary condition is P=0 and ρ→0 continuously; a shock with ρ_- > 0 and v_- = U ahead of vacuum is not a valid Rankine–Hugoniot shock because the mass and momentum flux conditions force P_- = 0. More decisively, the authors' own solution does not conserve energy. In 1D, Eq. (12) gives ρ = t^{-b} f(ξ), v = b x/t (since g(ξ)=bξ and b+c=1). The total energy E = ∫[ρv²/2 + ρ^{δ+1}/(δ+1)] dx becomes E(t) = t^{2b-2} A + t^{-bδ} B, with A=(1/2)∫ξ²f dξ and B=(1/(δ+1))∫f^{δ+1}dξ. Since b=2/(δ+2), both exponents equal -2δ/(δ+2) < 0, so E(t) → 0 and is not conserved. The Euler equations with no external force conserve energy; hence the claimed scaling cannot describe the long-time evolution of an isolated cloud. The numerics in the paper use artificial viscosity (Eq. A1), which dissipates energy and can produce apparent collapse; the data do not test the inviscid claim. The same energy decay occurs for the 2D solution (Eq. (71)). Moreover, the 2D derivation of the second jump condition uses d/dt∫ r v dr=0, which is not a conserved quantity of the radial Euler equations, so the jump conditions themselves are internally suspect. The finite-background and trapped-gas sections may survive, but the headline vacuum-blast result collapses.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies long-time solutions of the zero-temperature Euler equations for a polytropic quantum gas with specific enthalpy w = ρ^δ. The central claim is that a finite-mass gas cloud released into vacuum develops a shock front and attains self-similar form with R(t) ~ t^b, where b = 2/(δ+2) in 1D and b = 1/(δ+1) in 2D, and the paper provides explicit scaling functions. The authors also analyze expansion into a finite background, where they find ballistic sound-wave propagation, and the relaxation of a trapped gas to a steady state through damped oscillations. Numerical simulations with artificial viscosity are presented as validation, and code is provided.","tokens_in":16147,"tokens_out":10704,"duration_ms":128859,"significance":"The finite-background analysis and the trapped-gas relaxation calculation are potentially useful and give concrete, testable predictions; the paper also has the virtue of providing reproducible code and detailed numerical comparisons. However, the central vacuum-blast claim—which forms the paper's main advertised result—is not a solution of the Euler equations: the claimed scaling violates energy conservation and the shock boundary condition is not a valid Rankine-Hugoniot condition. Because this flaw affects the central contribution, the significance of the manuscript as it stands is limited.","major_comments":[{"comment":"The claimed self-similar solution does not conserve energy. With ρ = t^{-b} f(x/t^b) and v = b x/t, the total energy E = ∫[ρv²/2 + ρ^{δ+1}/(δ+1)] dx equals t^{2b-2} A + t^{-bδ} B. For b = 2/(δ+2), both exponents equal -2δ/(δ+2) < 0, so E(t) → 0. The Euler equations with no external force conserve energy for an isolated finite-mass cloud, so this solution cannot be the long-time limit of the inviscid dynamics. After pressure becomes negligible, the correct long-time behavior is ballistic free streaming, R(t) ~ t.","section":"§II.A, Eqs. (5)-(14)"},{"comment":"The Rankine-Hugoniot conditions used to select the solution are not valid for expansion into vacuum. Across a shock with vacuum on the plus side, ρ_+ = v_+ = P_+ = 0, mass and momentum conservation impose v_- = U and P_- = 0, meaning a finite-pressure shock into vacuum cannot exist. The paper instead imposes ρ_-^δ = U²/2 (Eq. (11); see Eq. (10) and Eq. (65)). A free boundary of a polytropic gas in vacuum is a rarefaction/contact boundary with ρ → 0 continuously, not a shock. Thus the shock-front premise underlying the entire vacuum-blast construction is not a consequence of the Euler equations.","section":"§II.A, Eq. (10); §IV, Eqs. (60)-(66)"},{"comment":"The derivation of the second jump condition in 2D uses d/dt∫ r v dr = 0, but ∫ r v dr is not conserved by the radial Euler equations (54). Only the mass-weighted moment, d/dt∫ r ρ dr = 0, is a conservation law. Therefore Eq. (64) and the resulting relation f(ξ_f) = (b² ξ_f²/2)^{1/δ} are not derivable. This invalidates the 2D scaling solution independently of the 1D objections.","section":"§IV, Eqs. (58)-(64)"},{"comment":"The numerical validation solves dissipative equations with artificial viscosity η∂²ρ and ν∂²v. Since the claimed Euler scaling dissipates energy at a rate t^{-2δ/(δ+2)}, agreement between DNS and the analytic formulas demonstrates only that the viscous system possesses the corresponding self-similar attractor; it does not test the inviscid claim. The assertion that 'the dissipation terms become irrelevant in the long time scaling regime' is not supported: finite ν provides exactly the energy sink needed to reach the decaying-energy solution.","section":"Appendix A, Eq. (A1)"}],"minor_comments":[{"comment":"The expression for the second exponent is garbled. From b + c = 1, one infers c = 1 - b = δ/(δ+2). Please correct the equation and the surrounding text.","section":"Eq. (6)"},{"comment":"The parentheses in the formula for f(ξ) are ambiguous. It should read f(ξ) = [ (ξ_f²(2−δ) + δ ξ²) / (δ+2)² ]^{1/δ}, assuming that is the intended expression.","section":"Eq. (12a)"},{"comment":"The second Rankine-Hugoniot condition is written as v_-²/2 + ρ_-^δ/v_- = U, which is dimensionally unclear and is not a standard jump condition. Please rewrite it in terms of P_- and U (or remove it, since it is not a correct RH relation).","section":"Eq. (10)"},{"comment":"The limit ρ_B → 0 is singular: the linearized sound speed c_s → 0, so the finite-background analysis cannot be expected to recover the vacuum exponents. The curve in Fig. 2(c) should be presented with this caveat.","section":"§II.B, Fig. 2(c)"},{"comment":"The replacement of the density-dependent viscous term by ν∂²_x v is acknowledged, but the claim that this is a good approximation for shallow traps requires more support. It would be helpful to quantify the error or to compare with the density-dependent form in a simple test.","section":"§III, Eq. (32b)"}],"recommendation":"reject","confidential_remarks":"The central vacuum-blast claim is not a valid solution of the Euler equations; this is a load-bearing error that cannot be fixed by local edits. The finite-background and trapped-gas sections are internally more sound and could potentially form separate manuscripts, but the present title and abstract make the vacuum result the main contribution. I recommend rejection rather than major revision, because the advertised scaling exponents and scaling functions for free expansion into vacuum are wrong and the required replacement (a rarefaction/free-streaming solution with R~t) would change the paper's central conclusions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nWhat you should know: the paper's headline claim—that a gas released into vacuum develops a shock with R(t) ~ t^b, b=2/(δ+2) in 1D—does not hold up against the inviscid Euler equations. The finite-background and trapped-gas sections are considerably better and may stand on their own.\n\nWhat is new and good. The finite-background analysis (Sec II B) gives a clean ballistic scaling and a concrete 1/√t correction with the right coefficient; that looks correct and is a nice application to quantum gases. The trapped-gas relaxation (Sec III) maps the linearized eigenvalue problem onto associated Legendre functions and computes complex eigenfrequencies that match the nonlinear numerics; that is solid, checkable work. The paper also ships code and is honest about dissipative regularization, including the caveat that viscosity matters for the trapped gas.\n\nThe soft spot is the vacuum blast (Secs II A and IV). The free boundary of a polytropic gas expanding into vacuum has P=0 and ρ=0 at the front; a shock with ρ_->0 moving into vacuum violates the momentum Rankine–Hugoniot condition unless P_-=0, and Eq. (10) is not the standard jump condition. More decisively, the explicit scaling solution (Eq. 12) does not conserve energy. With ρ=t^{-b}f(x/t^b) and v=b x/t, the total energy scales as t^{-2δ/(δ+2)} and decays to zero; the inviscid Euler equations conserve energy, so this cannot be the late-time solution of an isolated release. The same defect appears in 2D, and there the derivation of the second jump condition uses ∫ r v dr as a conserved quantity, which it is not. The numerical evidence comes from simulations with artificial viscosity; those can produce shock-like fronts, but they do not test the inviscid claim.\n\nThe novelty is also overstated. For a polytropic equation of state, the TvNS scaling and exponents are classical; mapping δ onto quantum-gas equations of state is an application, not a new solution. The paper itself cites earlier quantum shock work (Refs. 14, 37, 38), so the claim that “all work has focused on classical systems” is too strong.\n\nBottom line: the central claim is wrong, and I would not rely on it in my own work. But the rest of the paper contains useful, correct pieces, and the authors are clearheaded about the numerics. I would send it to a referee rather than desk-reject it, with the expectation that the vacuum-blast sections be revised or removed. The finite-background and trapped-gas results deserve review and are likely citable.","headline":"The vacuum-blast headline does not conserve energy; the finite-background and trapped-gas sections are the real value.","tokens_in":16568,"tokens_out":4257,"would_cite":false,"duration_ms":52006,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A quantum gas released from a trap is claimed to expand as a self-similar shock front, with cloud radius growing as t^(2/(2+δ)) (1D) or t^(1/(1+δ)) (2D), set by the gas's equation of state.","keywords":["quantum gases","Euler hydrodynamics","self-similar solutions","shock front","blast wave","time-of-flight expansion","relaxation dynamics","power-law equation of state"],"falsifier":"Release a finite-energy gas with power-law enthalpy (e.g., δ = 1) into vacuum and integrate the Euler equations with vanishing artificial viscosity. The claim requires a genuine front shock with nonzero density jump (ρ_-^δ = U²/2) and sub-ballistic radius R ~ t^(2/(2+δ)); the classical alternative is a rarefaction edge with density → 0 and ballistic R ~ t fixed by energy conservation. Experimentally, time-of-flight expansion of a quasi-1D unitary Fermi gas (δ = 2/5) distinguishes b = 5/6 from ballistic b = 1 by measuring R(t).","tokens_in":15543,"feed_emoji":"💥","tokens_out":24706,"duration_ms":235204,"temperature":0.7,"pith_summary":"This paper studies what happens when a zero-temperature quantum gas — bosons or fermions — is suddenly released from a trap, as in standard time-of-flight experiments. Using Euler hydrodynamics with a power-law equation of state (specific enthalpy w = Dρ^δ), the authors claim that a localized cloud expanding into vacuum develops a shock front at its leading edge and settles into self-similar motion at long times: density and velocity take the forms ρ = t^(-b) f(x/t^b), v = t^(-c) g(x/t^b), with b = 2/(2+δ) in one dimension and b = 1/(1+δ) in two dimensions. The full scaling functions are derived in closed form and verified by direct numerical simulation. Expansion into a finite-density background is instead ballistic, with the front moving at the sound speed; and a trapped gas relaxes to its steady state through damped oscillations whose frequency and decay rate follow from a linearized spectral problem. If correct, the results make the expansion of a quantum gas a direct, quantitative probe of its equation of state and connect cold-atom dynamics to the classic blast-wave problem.","feed_headline":"Expanding quantum gases sharpen into shock fronts","feed_subtitle":"Cloud size should grow as t^(2/(2+δ)) in 1D and t^(1/(1+δ)) in 2D, a direct readout of the gas's equation of state.","key_machinery":"The carrying mechanism is the self-similar scaling Ansatz ρ = t^(-b) f(x/t^b), v = t^(-c) g(x/t^b) applied to the Euler equations with power-law enthalpy (in 2D the density prefactor is t^(-2b) and coordinates are radial). It converts the partial differential equations into ordinary differential equations for f and g, with time-independence fixing the exponents b and c. The Rankine-Hugoniot relations at the expanding front collapse to ρ_-^δ = U²/2, and mass conservation fixes the scaled front position ξ_f, closing the problem exactly — the same structure as classical blast waves.","core_discovery":"At zero temperature the hydrodynamics of a quantum gas reduces to Euler equations with power-law enthalpy w = Dρ^δ. Claim: a localized mass released into vacuum develops a long-time self-similar shock-fronted solution ρ = t^(-b) f(x/t^b), v = t^(-c) g(x/t^b), with b = 2/(2+δ) in 1D, 1/(1+δ) in 2D. Rankine-Hugoniot conditions at the front reduce to ρ_-^δ = U²/2; the scaling functions close in algebraic form, and mass conservation fixes the front position. Numerical solution of the Euler equations for δ = 2/5 in one and two dimensions collapses onto these curves. The same framework gives ballistic, sound-speed-limited expansion into a finite background, and damped-oscillatory relaxation to a s","pith_inferences":["The 1D and 2D exponents are both consistent with the general-dimensional form b = 2/(2 + dδ); if so, a 3D unitary Fermi gas (δ = 2/3) should expand with R ~ t^(1/2), a quantitative prediction the paper does not state.","For an isolated cloud the total energy is conserved, yet under the sub-ballistic scaling both kinetic and internal energy decay to zero at long times; how this closes — via heating, a breakdown of the zero-temperature equation of state, or a crossover to ballistic motion — is the question the analysis leaves implicit.","The relaxation spectrum depends on viscosity and trap frequency only through combinations such as ω²/λ, so precise measurements of the damped oscillations of a trapped gas could extract an effective hydrodynamic viscosity, which the paper notes is not yet fixed from theory."],"forward_implications":["In time-of-flight expansion into vacuum, the cloud radius should grow as R(t) ~ t^(2/(2+δ)) in 1D and t^(1/(1+δ)) in 2D — measurably sub-ballistic — with the exponent fixed purely by the equation-of-state exponent δ.","Density and velocity profiles at long times should collapse onto the paper's closed-form scaling functions after rescaling by t^b, independent of the shape of the initial localized cloud.","For release into a finite background density ρ_B, the excess density propagates at the sound speed c_s = sqrt(δρ_B^δ) with an explicit x/t profile, and the front perturbation height decays as 1/√t.","A trapped gas excited out of equilibrium relaxes to the steady state ρ_ss ∝ (1 − x²/a²)^(1/δ) with damped oscillations whose decay rate and frequency are set by the lowest eigenvalue of a Legendre-type spectral problem.","Because the enthalpy exponent δ takes different values for Lieb-Liniger bosons, unitary Fermi gases, and finite-range Riesz gases, the same scaling laws should appear across these distinct quantum fluids, with R(t) ~ t^b a direct signature of which gas one has."],"supporting_citations":[{"why":"The classic self-similar blast-wave solutions whose Ansatz and shock-matching logic the paper adapts to quantum gases.","marker":"[41–44]"},{"why":"The similarity-and-dimensional-methods framework that supplies the scaling construction for the asymptotic solutions.","marker":"[40]"},{"why":"Rankine-Hugoniot blast-front conditions in a one-dimensional cold gas, used to close the vacuum-expansion boundary problem.","marker":"[50, 51]"},{"why":"Experimental observation of shock waves in a strongly interacting Fermi gas and the hydrodynamic theory giving the δ = 2/5 equation of state for the quasi-1D unitary gas.","marker":"[14, 17]"},{"why":"Shows a finite-range Riesz gas maps to the same hydrodynamic form with δ = s, extending the power-law enthalpy class to arbitrary exponents.","marker":"[57]"},{"why":"Supplies the radial boundary conditions and Euler-versus-Navier-Stokes comparison used for the two-dimensional numerical solution.","marker":"[72]"}],"fun_headline_variants":["Quantum gases punch shock fronts on expansion","Expanding quantum gas clouds form shock fronts","Shock fronts emerge in quantum gas hydrodynamics","Quantum gas expansion obeys power-law scaling","Hydrodynamics predicts quantum gas shock fronts"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The construction hinges on the leading edge of a gas expanding into vacuum being a genuine shock with a nonzero density jump obeying ρ_-^δ = U²/2, rather than a free rarefaction edge where density and pressure vanish; if the edge is a rarefaction, the sub-ballistic exponents R ~ t^(2/(2+δ)) give way to ballistic growth set by energy conservation.","fun_headline_variants_meta":{"raw":{"variants":["Quantum gases punch shock fronts on expansion","Expanding quantum gas clouds form shock fronts","Shock fronts emerge in quantum gas hydrodynamics","Quantum gas expansion obeys power-law scaling","Hydrodynamics predicts quantum gas shock fronts"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000175,"raw_usage":{"total_tokens":1112,"prompt_tokens":725,"completion_tokens":387,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":469,"completion_tokens_details":{"reasoning_tokens":322}},"tokens_in":469,"tokens_out":387,"duration_ms":5083,"temperature":1.0,"reasoning_tokens":322,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T13:41:31.013252+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Release a finite-energy gas with power-law enthalpy (e.g., δ = 1) into vacuum and integrate the Euler equations with vanishing artificial viscosity. The claim requires a genuine front shock with nonzero density jump (ρ_-^δ = U²/2) and sub-ballistic radius R ~ t^(2/(2+δ)); the classical alternative is a rarefaction edge with density → 0 and ballistic R ~ t fixed by energy conservation. Experimentally, time-of-flight expansion of a quasi-1D unitary Fermi gas (δ = 2/5) distinguishes b = 5/6 from ballistic b = 1 by measuring R(t).","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The similarity-and-dimensional-methods framework that supplies the scaling construction for the asymptotic solutions."},{"cited_title":"Kumar, M","cited_arxiv_id":null,"evidence_quote":"Shows a finite-range Riesz gas maps to the same hydrodynamic form with δ = s, extending the power-law enthalpy class to arbitrary exponents."},{"cited_title":"Kumar and R","cited_arxiv_id":null,"evidence_quote":"Supplies the radial boundary conditions and Euler-versus-Navier-Stokes comparison used for the two-dimensional numerical solution."}],"review_version":1}