{"id":"b5b135e4-848e-46a6-98c4-0dfab15e4b83","arxiv_id":"1908.01291","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"This paper derives exact finite-temperature dynamics and stability boundaries for a Tonks-Girardeau gas in a periodically modulated harmonic trap, with both quantum and hydrodynamic descriptions reducing to the same Mathieu equation.","lead":"Physicists worked out exactly how a one-dimensional gas of hard-core particles, the Tonks-Girardeau gas, responds when the trap holding it is modulated at a fixed frequency. The results predict when the gas stays stable or grows without bound, and show that a single classical equation governs the quantum dynamics at any temperature.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified","rationale":"I re-derived the Mathieu mapping and checked Eq. (26) against the Pinney initial conditions and constraint; no algebraic inconsistency appears. I also checked the stability bijection: in unstable Mathieu bands, no choice of coefficients A, B, C can cancel the growing Floquet component, because the constraint AB-C^2 = omega0^2/W^2 would require a negative squared term. The hydrodynamic comparison is structurally identical because both approaches use the same lambda equation. The only substantive limitation is the physical breakdown of the TG model in unstable regimes, exactly as the authors state in Sec. VII; this is a scope condition on the central claim, not a flaw in the claimed exact TG dynamics. The numerical examples are illustrative and do not affect the stability diagram's validity. No critical red flags were found, and the reader's verdict of ACCEPT remains appropriate.","tokens_in":19181,"tokens_out":19843,"duration_ms":220000,"concrete_test":"Verify Fig. 3 numerically: integrate Eq. (11) directly with a high-accuracy adaptive integrator over a grid in (Omega/omega0, alpha) for long times, and compare the bounded/unbounded classification of lambda(t) with the Mathieu Floquet exponent computed from Eqs. (30)-(31). A mismatch in any region would falsify the claimed bijection.","verdict_should_be":"UNCHANGED","load_bearing_attack":"After checking the central construction, I do not find a load-bearing defect. The Ermakov-Pinney to Mathieu mapping (Eqs. 22-26) is internally consistent: Eq. (26) satisfies the initial conditions and the Pinney constraint AB-C^2 = omega0^2/W^2, and in unstable Mathieu bands the growing Floquet component cannot be cancelled by the coefficients A, B, C because the constraint forbids the required zero amplitude. Thus the lambda-stability bijection holds. The acknowledged breakdown of the TG model in the unstable regime (Sec. VII) is a physical applicability limit, not an internal inconsistency, and the paper states it explicitly.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a one-dimensional gas of hard-core bosons in a harmonic trap whose frequency is modulated as ω²(t)=ω0²(1−α sin Ωt). The authors use the Fermi–Bose mapping and the scaling solution for harmonic trapping to show that the finite-temperature reduced one-body density matrix, and hence the density and momentum distributions, are governed by a single scaling parameter λ(t) satisfying the Ermakov–Pinney equation. They map the linearized equation for λ to Mathieu's equation with parameters a=(2ω0/Ω)² and q=2ω0²α/Ω², construct λ(t) from even and odd Mathieu functions via Pinney's formula, and compute Floquet exponents to obtain the stability diagram in the (a,α) plane. They show that the finite-temperature hydrodynamic equations admit the same scaling parameter and therefore the same stability diagram. Numerical examples in stable and unstable regimes compare the exact many-body dynamics with the hydrodynamic approximation at two temperatures, and the paper discusses parametric resonances and the many-body bounce.","tokens_in":19239,"tokens_out":20991,"duration_ms":212863,"significance":"This is a clean and useful result. If accepted, it reduces a nontrivial finite-temperature interacting many-body problem to a single linear ODE with a well-known stability chart, with no fitted parameters. The central mapping is internally consistent: Eqs. (22)–(26) follow from the Schrödinger equation and Pinney's construction, and the positive-definite constraint AB−C²=ω0²/W² in Eq. (14) rules out cancellation of the growing Floquet component in unstable Mathieu bands, so the λ-stability bijection is sound. The paper also makes a valuable methodological point that the hydrodynamic approach reproduces the exact stability diagram exactly. The numerical work is honest about the Lorentzian approximation in the hydrodynamic momentum distribution (Appendix A) and about the breakdown of the TG description at exponentially growing densities in the unstable regime (Sec. VII). The latter is a physical applicability limit explicitly acknowledged by the authors, not an internal inconsistency. The stability diagram in Fig. 3 and the resonance condition Ω_j=2ω0/j are crisp falsifiable predictions for experiments.","major_comments":[],"minor_comments":[{"comment":"The phase factor in the scaling solution for the single-particle orbitals should read −iE_j(t)t/ℏ; as typeset, −iE_j(t)t is dimensionally incorrect.","section":"II.B, Eq. (10)"},{"comment":"The relation between Mathieu parameters is written q=α/2 a, which is ambiguous; it should be q=α a/2.","section":"V, Eq. (36)"},{"comment":"The symbol C is used both for the Pinney coefficient in Eq. (14) and for the even Mathieu function C(a,q,τ); please disambiguate for clarity.","section":"II.B, Eq. (14), and III, Eqs. (24)–(25)"},{"comment":"The sentence that stable solutions correspond to ν being real is imprecise at integer ν=0,1, where the second solution is nonperiodic and unbounded; the surrounding discussion handles this, but the statement should be qualified.","section":"IV.B"},{"comment":"The reference to Fig. 3(a) should be to Fig. 3, since that figure has no panel (a).","section":"VI.A"}],"recommendation":"minor_revision","confidential_remarks":"The manuscript is a solid contribution; the derivation is sound and the novelty relative to Ref. [20] is clear. My requests are cosmetic; after a minor revision it should be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hey,\n\nQuick take on Atas, Simmons, Kheruntsyan (arXiv:1908.01291): this is a genuinely useful paper, and the central result holds up. The authors give an exact finite-temperature treatment of a Tonks-Girardeau gas in a sinusoidally modulated harmonic trap, reducing the many-body dynamics to a single Ermakov–Pinney equation and mapping its stability to Mathieu's equation. The resulting stability diagram in (a, α) space, with the primary and higher-order parametric resonances, is clean and physically well-posed. I checked the construction: the scaling solution, the Ermakov–Pinney to Mathieu mapping, and the Floquet-based stability analysis are internally consistent. The stress-test note is right that the AB−C² constraint forbids cancellation of the growing Floquet component in unstable bands, so the λ-stability bijection is solid.\n\nWhat is actually new: for this driven protocol, they provide explicit exact dynamics of the real-space density and momentum distribution at finite temperature, including the many-body bounce phenomenon in both stable and unstable regimes. Quinn and Haque had zero-T energy absorption; Atas et al. had quenches. The periodic-modulation finite-T case with explicit stability boundaries is a real addition. The hydrodynamic comparison is also a plus: same Ermakov–Pinney equation, same Mathieu mapping, same stability diagram. No fitted parameters, no circularity.\n\nSoft spots, in proportion: the hydrodynamic momentum distribution uses a low-temperature Lorentzian approximation for the local momentum distribution (Appendix A). It captures the qualitative physics and the bulk of n(k,t), but quantitative accuracy degrades at higher temperatures and at large momentum; the authors are candid about that. The other caveat is in Section VII: in the unstable regime, the peak density grows exponentially, so the local Lieb–Liniger parameter γ(x,t) eventually drops to order unity and the TG model stops being a good description. The authors state this explicitly. It's a physical applicability limit, not an internal contradiction. Minor: the Mathieu-function sections are dense and somewhat encyclopedic, but the essential derivation is standard and the appendices are complete.\n\nThe citation pattern looks honest: the scaling formalism comes from Ref. [23] and the quench paper [24], and self-citation there is justified because those are the sources of the technique. No evidence of overclaiming.\n\nBottom line: solid paper, clear value for ultracold-atom and out-of-equilibrium many-body theory. It deserves a serious referee, and I would expect acceptance after minor revisions. I'd cite it if I were working in TG dynamics.\n\nRecommend sending to peer review.","headline":"Clean exact finite-T solution for driven TG gas with a solid Mathieu stability mapping; worth refereeing, with the unstable-regime validity issue properly flagged.","tokens_in":19741,"tokens_out":2368,"would_cite":true,"duration_ms":24894,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves that the finite-temperature dynamics of a periodically modulated Tonks-Girardeau gas collapses into a single scaling parameter whose stability boundaries are those of Mathieu's equation.","keywords":["Tonks-Girardeau gas","finite-temperature dynamics","parametric resonance","Mathieu equation","Ermakov-Pinney equation","scaling solution","momentum distribution","hydrodynamic approach"],"falsifier":"Drive a finite-$\\gamma$ Lieb-Liniger gas (or a time-dependent matrix-product-state simulation of the same protocol) with the same sinusoidal trap modulation, fix a parameter point $(a,\\alpha)$ predicted to be stable and one predicted unstable, and compare the measured $\\lambda(t)$ and momentum width against the Mathieu-based prediction. A clear failure is if the $\\Omega=2\\omega_0$ resonance requires a nonzero threshold amplitude or if the $j=3$ resonance appears at $a=9$ rather than shifted upward; either observation would show that the Tonks-Girardeau scaling mapping stops capturing the boundaries before the gas leaves the infinite-repulsion regime.","tokens_in":18998,"feed_emoji":"⚛️","tokens_out":7179,"duration_ms":68556,"temperature":0.7,"pith_summary":"This paper seeks to show that the out-of-equilibrium dynamics of a finite-temperature, harmonically trapped Tonks-Girardeau gas under sinusoidal modulation of the trap frequency is governed entirely by a single scaling parameter $\\lambda(t)$, and that the stability of the many-body system is exactly the stability of Mathieu's equation. Because both the exact quantum many-body solution and the finite-temperature hydrodynamic equations collapse onto the same $\\lambda(t)$, both approaches produce the identical stability diagram and identical parametric resonances. That allows the paper to give explicit exact formulas for the density and momentum distributions, and to interpret the driven gas as a parametric oscillator whose resonant frequencies align with the breathing modes of the trap. A sympathetic reader would value this as a rare example where a genuinely interacting many-body problem at finite temperature is solved exactly and reduced to a single classical stability problem.","feed_headline":"Driven quantum gas maps exactly onto Mathieu's equation","feed_subtitle":"For a Tonks-Girardeau gas, quantum and hydrodynamic theories share the exact same resonance boundaries.","key_machinery":"The load-bearing device is the scaling solution for a time-dependent harmonic trap: $\\varphi_j(x,t)=\\lambda^{-1/2}\\varphi_j(x/\\lambda,0)\\exp\\left[\\frac{imx^2}{2\\hbar}\\frac{\\dot\\lambda}{\\lambda}-iE_j(t)t\\right]$, together with the Ermakov-Pinney equation and Pinney's construction of $\\lambda$ as $\\sqrt{A\\lambda_1^2+B\\lambda_2^2+2C\\lambda_1\\lambda_2}$ from two independent solutions of the homogeneous oscillator equation. Under the change of variable $\\Omega t=\\pi/2-2\\tau$, that homogeneous equation becomes Mathieu's equation, an ordinary linear oscillator with periodic coefficients of the form $\\ddot z+(a-2q\\cos 2\\tau)z=0$. The Floquet stability chart of this equation then determines whether the gas's density and momentum distributions stay bounded or grow exponentially.","core_discovery":"The central claim is that for $\\omega^2(t)=\\omega_0^2(1-\\alpha\\sin\\Omega t)$, the exact finite-temperature many-body dynamics of the trapped Tonks-Girardeau gas is fully determined by $\\lambda(t)$, the positive solution of the Ermakov-Pinney equation $\\ddot\\lambda+\\omega^2(t)\\lambda=\\omega_0^2/\\lambda^3$. The paper constructs $\\lambda(t)$ from two independent Mathieu functions $C$ and $S$ via Pinney's formula and proves that the long-time behavior is in direct bijection with the stability of Mathieu's equation with $a=(2\\omega_0/\\Omega)^2$ and $q=2\\omega_0^2\\alpha/\\Omega^2$. The same $\\lambda(t)$ appears in the finite-temperature hydrodynamic theory, so the stability diagram and the structure of parametric resonances are exactly the same in both descriptions. The paper also provides closed-form expressions for the Floquet exponent $\\nu(a,q)$, which sets the rate of exponential growth in the unstable regions.","pith_inferences":["If the same Ermakov-Pinney structure survives in a weakly anharmonic trap, the Mathieu stability boundaries themselves could serve as an observable signature of integrability; testing them with finite-$g$ Lieb-Liniger gases would quantify how quickly beyond-Tonks-Girardeau corrections destroy the exact resonance geometry.","The closed-form Floquet exponent could be used to define a heating rate per driving cycle in the unstable regime, giving a quantitative figure of merit for future quantum heat engine proposals on this system; the paper gestures at that direction but does not develop it.","Cosine modulation, obtainable from these solutions by a time shift and $q\\to -q$, makes the same stability diagram apply to experiments that switch the drive on smoothly, so the resonance structure is likely insensitive to the exact turn-on waveform."],"forward_implications":["A gas driven at $\\Omega=2\\omega_0$, the primary breathing-mode resonance, is unstable for every modulation amplitude $\\alpha>0$; the higher resonances $\\Omega_j=2\\omega_0/j$ require a finite $\\alpha$ and their centers shift downward in frequency.","Because the hydrodynamic equations share the same scaling parameter and the same Mathieu map, the stability boundaries and resonance widths are identical between the classical-fluid and exact-quantum descriptions, not merely approximately equal.","The collective many-body bounce, an extra narrowing of the momentum distribution within each breathing cycle, persists under periodic driving in both stable and unstable regimes, and increasing temperature blurs it.","The Floquet exponent $\\nu(a,q)$ is the single number controlling the long-time fate of every observable: a real $\\nu$ means bounded oscillations, while a complex $\\nu$ means exponential growth with a known rate."],"supporting_citations":[{"why":"Supplies the exact finite-temperature Fredholm-determinant formalism used to compute the one-body density matrix and observables.","marker":"[23]"},{"why":"Introduced the quench dynamics and the many-body bounce that this work extends to periodic driving.","marker":"[24]"},{"why":"Earlier zero-temperature study of periodically modulated trapped bosons whose resonance structure is here extended to the finite-temperature Tonks-Girardeau limit.","marker":"[20]"},{"why":"Established the scaling solution for a trapped Tonks-Girardeau gas in a time-dependent harmonic trap, giving Eq. (10).","marker":"[27]"},{"why":"Original source of the Ermakov equation that becomes the central scaling equation.","marker":"[28]"},{"why":"Pinney's general solution formula that constructs $\\lambda(t)$ from two independent homogeneous solutions.","marker":"[29]"},{"why":"Mathieu's equation in its canonical form, the target of the mapping.","marker":"[21]"},{"why":"Standard theory of Mathieu functions and Floquet exponents used throughout the stability analysis.","marker":"[22]"},{"why":"Source of the finite-temperature hydrodynamic equations that share the same scaling parameter and Ermakov-Pinney equation.","marker":"[26]"}],"fun_headline_variants":["Exact Mathieu map for driven Tonks-Girardeau gas","Quantum and fluid theories agree on resonance boundaries","Finite-temp Tonks-Girardeau matches Mathieu stability","Parametric resonances from Mathieu equation exact","Driven gas: Mathieu equation governs dynamics exactly"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central assumption is that the gas remains in the Tonks-Girardeau (infinite-repulsion) regime at all times; in unstable parameter regions the density peaks grow exponentially, so the local dimensionless interaction strength can fall to order one or below and the model stops being valid there.","fun_headline_variants_meta":{"raw":{"variants":["Exact Mathieu map for driven Tonks-Girardeau gas","Quantum and fluid theories agree on resonance boundaries","Finite-temp Tonks-Girardeau matches Mathieu stability","Parametric resonances from Mathieu equation exact","Driven gas: Mathieu equation governs dynamics exactly"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00015,"raw_usage":{"total_tokens":1177,"prompt_tokens":903,"completion_tokens":274,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":519,"completion_tokens_details":{"reasoning_tokens":193}},"tokens_in":519,"tokens_out":274,"duration_ms":3505,"temperature":1.0,"reasoning_tokens":193,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:17:02.540154+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Drive a finite-$\\gamma$ Lieb-Liniger gas (or a time-dependent matrix-product-state simulation of the same protocol) with the same sinusoidal trap modulation, fix a parameter point $(a,\\alpha)$ predicted to be stable and one predicted unstable, and compare the measured $\\lambda(t)$ and momentum width against the Mathieu-based prediction. A clear failure is if the $\\Omega=2\\omega_0$ resonance requires a nonzero threshold amplitude or if the $j=3$ resonance appears at $a=9$ rather than shifted upward; either observation would show that the Tonks-Girardeau scaling mapping stops capturing the boundaries before the gas leaves the infinite-repulsion regime.","supporting_citations":[{"cited_title":"Thorwart, P","cited_arxiv_id":null,"evidence_quote":"Introduced the quench dynamics and the many-body bounce that this work extends to periodic driving."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier zero-temperature study of periodically modulated trapped bosons whose resonance structure is here extended to the finite-temperature Tonks-Girardeau limit."},{"cited_title":"McLachlan, Theory and Application of Mathieu Func- tions, Dover Publications (Dover, 1964)","cited_arxiv_id":null,"evidence_quote":"Established the scaling solution for a trapped Tonks-Girardeau gas in a time-dependent harmonic trap, giving Eq. (10)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Original source of the Ermakov equation that becomes the central scaling equation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Pinney's general solution formula that constructs $\\lambda(t)$ from two independent homogeneous solutions."},{"cited_title":"Hnggi and C","cited_arxiv_id":null,"evidence_quote":"Mathieu's equation in its canonical form, the target of the mapping."},{"cited_title":"Mathieu, Journal de Math´ ematiques Pures et Ap- pliqu´ ees13, 137 (1868)","cited_arxiv_id":null,"evidence_quote":"Source of the finite-temperature hydrodynamic equations that share the same scaling parameter and Ermakov-Pinney equation."}],"review_version":1}