{"id":"ff8e7bf5-d81a-44b2-a2d2-c1b0bf00d71e","arxiv_id":"1908.06522","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"At late times, a weak non-integrable interaction makes a quantum shock decay exponentially and adds a universal low-energy correction that changes the density profile in a window behind the shock front.","lead":"What happens to a quantum shock wave when weak interactions are switched on? The authors show that at late times the shock decays exponentially with a universal correction that preserves its shape away from the front.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified.","rationale":"The reader's verdict ACCEPT is reasonable. The paper derives the central claim transparently from the free-fermion shock solution plus a weak-interaction kinetic closure. The main assumptions are the Boltzmann linearization, the single-particle decay rate, and the high-branch-only source approximation; the paper checks the source consistency in Appendix C and the other corrections are suppressed by the stated small parameters. The algebraic path from Eq. (29) to Eq. (36) is internally consistent under the stated λ≪1, tSΓ(km)≪1, and λτ≫1 conditions. The residual approximations affect only subleading regions of the γ integral, and the numerical consistency estimate, while not an asymptotic control, is sufficiently small for the claimed model. I found no specific technical step on which the central claim fails.","tokens_in":11671,"tokens_out":51039,"duration_ms":546017,"concrete_test":"Directly solve the linearized Boltzmann equation (16) with the kernel (24) and the free-shock initial condition, without the midpoint replacement of Eq. (25), and compare δρ(λ,t) at λtΓ(km)∼1 with Eq. (38). Agreement to within a few percent would close the last residual concern; disagreement would pinpoint the source-evaluation step as the failure point.","verdict_should_be":"UNCHANGED","load_bearing_attack":"No load-bearing flaw found in the derivation of Eq. (38). The least secure step is the linearized Boltzmann closure in Eqs. (16)–(21): the high-energy branch is treated as decaying independently with the single-particle rate Γ(k), and the source J is restricted to the high-energy branch. The paper itself flags, at the end of Sec. III and in App. C, that the residual low-energy feedback into J is not asymptotically small but is numerically small, of order 2–3% for the kernel of Eq. (24). That check, together with the parametric controls tSΓ(km)≪1 and λτ≫1 used to justify the approximations in Eqs. (25) and (33), is adequate within the stated validity window. The midpoint/averaging replacement in Eq. (25) is inaccurate only in the boundary layer γ≲√(λ/τ), whose contribution to δρ_low is suppressed by powers of λ; the dominant γ∼1 region is smooth and accurately captured. The stated validity limits, including the eventual encroachment of quantum ripples, are explicitly acknowledged, and the central claim is confined to that regime.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the late-time fate of a quantum shock wave in a one-dimensional spinless fermion system after weak generic (non-integrable) interactions are turned on. For free fermions, the authors recap the known result that the shock front carries quantum Airy-function ripples only in a parametrically small window, while the bulk of the shock is classical. They then add weak interactions through a linearized Boltzmann equation with a single-particle decay rate Γ(k) and compute the resulting density profile. The central result, Eq. (38), states that the density deviation is approximately (k_m−k_0)e^{−tΓ(k_m)}(t_S√λ)/(tπ)[1+F(λtΓ(k_m))], with F(z)=(8/(5√z))∫_0^z dy√y e^{−y}. Thus the shock height decays exponentially with the single-particle decay rate, while a universal low-energy correction reshapes the profile only in a window λtΓ(k_m)∼1; away from that window the free-fermion shape is preserved on a time-independent spatial scale. The result is explicitly restricted to t≫t_S, 1≫λ≫λ_cr, t_SΓ(k_m)≪1, and tΓ(k_m)≪(ΔN)^{2/3}. The paper includes a self-consistency check (Appendix C) showing that neglecting the low-energy feedback into the source term is numerically small, at most about 2%.","tokens_in":11867,"tokens_out":15955,"duration_ms":159576,"significance":"If the result holds, this is a valuable analytic prediction for a regime that is not accessible by the free-fermion Airy-function theory: it describes how a shock wave decays and deforms under generic weak interactions. The scaling function F(z) in Eq. (38) is parameter-free, and the only external input, the single-particle decay rate Γ(k_m), is taken from independent perturbation theory (Ref. 15). The derivation is explicit and the paper carefully states and respects its validity window, including the eventual breakdown when quantum ripples catch up with the kinetic correction. The main residual risk is that the linearized Boltzmann closure, especially the source approximation in Eq. (21), is controlled only by an a posteriori numerical check rather than a fully rigorous asymptotic estimate. However, the authors themselves flag this limitation at the end of Section III and in Appendix C, and their estimate that the neglected low-energy source contributes at most about 2% is adequate within the stated regime. Overall, the paper is a solid, self-contained contribution that identifies a falsifiable late-time scaling prediction.","major_comments":[],"minor_comments":[{"comment":"The symbol τ is used for two different quantities: in Eq. (30) it denotes (t′−t_1)/t_1, while in Eq. (33) it is reused as a dummy integration variable. Please rename one of them to avoid confusion.","section":"Equations (30) and (33)"},{"comment":"The integration limits in Eq. (33) are printed ambiguously as “∫λ/γ 0”; they should be typeset as ∫_0^{λ/γ}.","section":"Equation (33)"},{"comment":"Eq. (A12) contains an extra closing parenthesis after (k_m−k_0) and unbalanced parentheses in the oscillatory term; the typesetting should be corrected.","section":"Appendix A, Eq. (A12)"},{"comment":"The kernel W_{p→k} in Eq. (24) appears singular at p=k, while Γ(p) in Eq. (17) is defined as an integral of W up to p and is finite. A sentence clarifying the region of validity of the expression in Eq. (24) or its regularization would prevent an apparent inconsistency.","section":"Equation (24)"},{"comment":"The midpoint replacement in Eq. (25) is justified only by the condition t≫t_S√λ; it would be helpful to state explicitly that the error is confined to a boundary layer and is suppressed by powers of λ, as can be inferred from the estimate I_1=λ^5G in Appendix B.","section":"Equation (25)"}],"recommendation":"minor_revision","confidential_remarks":"The manuscript is a competent and self-contained theoretical derivation, with clear attribution of prior results. The self-citations (Refs. 6 and 15) are appropriate because the free-fermion shock formula and the single-particle decay rate are originally derived in those works. The main risk, the linearized kinetic closure, is explicitly acknowledged by the authors and checked to be numerically small, which I consider sufficient for publication. I recommend acceptance after minor revisions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this paper genuinely extends the free-fermion shock results to weakly interacting systems and produces a concrete, testable prediction. The new piece is Eq. (38): the late-time density profile is the exponentially decaying ballistic shock times a universal factor [1 + F(lambda t Gamma(k_m))]. The scaling function F(z) is new and, as the authors argue, its form is tied to the phase-space structure of the scattering kernel rather than details of the interaction. The basic picture—the shock front decays with the single-particle rate while low-energy back-reaction reshapes a window of width ~ w/(Gamma t_S) behind the front—is plausible and clearly stated.\n\nThe derivation is honest. They start from Bettelheim-Glazman, show that for t >> t_S the ripples occupy only a fraction (Delta N)^-2/3 of the shock, and then linearize a Boltzmann equation with a single-particle decay rate. They explicitly list the necessities: t_S Gamma << 1, lambda >> lambda_cr, and the eventual breakdown when ripples catch up at t Gamma ~ (Delta N)^2/3. Appendix C checks the consistency of dropping the low-energy source, finding corrections of at most 2%. The stress-test note agrees that the midpoint/averaging step is controlled outside a boundary layer suppressed by powers of lambda.\n\nSoft spots: the whole edifice assumes the linearized Boltzmann closure, which is a physical assumption and not derived from the microscopic Hamiltonian. No numerics or exact limits confirm the final F(z). The decay rate Gamma(k) comes from perturbation theory for a particular density-density model (Ref. [15]), so the universality of F is a conjecture, albeit a reasonable one. Also, F(z) being O(1) at z ~ 1 means the 'correction' is not a small perturbation in that window; the linearized equation is doing serious work there. But the paper flags this and checks numerically the smallness of the neglected source. These caveats are proportionate and do not undermine the main claim within the stated validity window.\n\nWho would benefit: anyone working on quantum hydrodynamics, generalized hydrodynamics, or thermalization in 1D. It is a solid theory paper that deserves a serious referee. I would accept it after a moderate revision that sharpens the validity domain of Eq. (38) and, if possible, adds a numerical check of the Boltzmann closure. The central result is likely correct within the approximations, and the paper is more transparent than most.","headline":"A careful analytic derivation that gives the first concrete picture of how weak interactions dissolve a free-fermion shock; the linearized Boltzmann closure is the softest step, but it is handled honestly.","tokens_in":12358,"tokens_out":3370,"would_cite":true,"duration_ms":33391,"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":"Weak generic interactions make a fermionic shock wave's amplitude decay exponentially while preserving its shape away from the front.","keywords":["shock waves","quantum hydrodynamics","fermions","Boltzmann equation","integrability breaking","single-particle decay","Wigner function","late-time dynamics"],"falsifier":"A numerical simulation of weakly interacting one-dimensional fermions (or of the full three-particle Boltzmann equation) initialized as a smooth density bump would settle the claim: plot $\\delta\\rho(\\lambda,t)/[(k_m-k_0)e^{-t\\Gamma(k_m)}t_S\\sqrt{\\lambda}/(t\\pi)]$ against $\\lambda t\\Gamma(k_m)$ for $t\\gg t_S$. Collapse onto $1+F(\\lambda t\\Gamma(k_m))$ supports Eq. (38); a systematic deviation, or a low-energy source contribution growing beyond the stated $\\sim 2\\%$, falsifies the linearization and decay-rate closure.","tokens_in":11484,"feed_emoji":"🌊","tokens_out":12010,"duration_ms":105255,"temperature":0.7,"pith_summary":"This paper studies the late-time fate of a shock wave formed in a one-dimensional Fermi gas by a smooth density bump, when generic non-integrable interactions are added to the free-fermion problem. By linearizing the three-particle Boltzmann equation around the free shock and using a perturbative single-particle decay rate $\\Gamma(k_m)$, the authors derive a closed expression for the density deviation: $\\delta\\rho \\approx (k_m-k_0)e^{-t\\Gamma(k_m)}\\frac{t_S\\sqrt{\\lambda}}{t\\pi}[1+F(\\lambda t\\Gamma(k_m))]$. The result says that the shock height shrinks exponentially at the single-particle decay rate, while the profile, in the dimensionless coordinate $\\lambda$, keeps the free-fermion shape except in a window $\\lambda t\\Gamma(k_m)\\sim 1$ where the universal correction $F$ matters. This matters because it gives a concrete analytic picture of how integrability-breaking interactions destroy a hydrodynamic singularity, and it shows that the naive exponentially decaying kinematic shock is valid over most of the structure.","feed_headline":"Quantum shocks decay exponentially but keep their shape","feed_subtitle":"A linearized Boltzmann equation fixes how weak interactions erase a fermionic shock: amplitude decays, shape survives.","key_machinery":"The carrying object is the linearized Boltzmann equation $(\\partial_t+\\frac{k}{m}\\partial_x)f=J-\\Gamma(k)f$, obtained by splitting the three-particle collision integral into a source $J$ from high-energy particles decaying into low energies and a loss term with the perturbative single-particle decay rate. The high-energy branch is the free-fermion step function times $e^{-t\\Gamma(k)}$, and $J$ is approximated by integrating that branch against $W_{p\\to k}\\propto (k-k_0)^2/(p-k)^5$. This turns the collisional problem into an integral along characteristics; the scaling variable $\\lambda t\\Gamma(k_m)$ emerges from the competition between ballistic motion and decay, and the function $F$ is the integral that survives averaging over the shock's momentum interval.","core_discovery":"The central claim is Eq. (38). Starting from the Wigner function and a local Fermi surface, the authors argue that for times well after shock formation the high-energy branch of the distribution is $\\theta(k-k_-(x,t))\\theta(k_+(x,t)-k)e^{-t\\Gamma(k)}$, and that this branch alone sources the low-energy particles through the golden-rule rate $W_{p\\to k}$. Solving the linearized kinetic equation along ballistic characteristics gives $\\delta\\rho \\approx (k_m-k_0)e^{-t\\Gamma(k_m)}\\frac{t_S\\sqrt{\\lambda}}{t\\pi}[1+F(\\lambda t\\Gamma(k_m))]$, with $F(z)=\\frac{8}{5\\sqrt{z}}\\int_0^z dy\\,\\sqrt{y}e^{-y}$. In plain terms: away from the front the shock keeps its free-fermion shape and decays exponentially, and only in the window $\\lambda t\\Gamma(k_m)\\sim 1$ does the low-energy correction reshape it, on a time-independent physical scale $x_+(t)-x\\sim w/(\\Gamma(k_m)t_S)$.","pith_inferences":["Extending the authors' stated conjecture, the same scaling function should control the late-time density of weakly interacting spinful fermions, because the $(p-k)^{-5}$ suppression comes from the small low-energy density of states and survives spin; a numerical simulation with spin would test whether $F$ appears with the same scaling argument.","Equation (38) implies that the shape of the shock survives for a time set by the decay rate rather than by diffusion, so in a weakly interacting one-dimensional gas the structure is erased on a scale $\\sim\\Gamma(k_m)^{-1}\\ln(\\Delta N)$ once the exponentially small amplitude is accounted for.","A cold-atom quench experiment that images the density after creating a localized bump could measure $F$ directly: plotting the rescaled deviation against $\\lambda t\\Gamma(k_m)$ should collapse data taken at different times onto one curve.","If a non-perturbative treatment in the spirit of the nonlinear Luttinger liquid framework mentioned in the paper produced a different scaling function, that would mark the point where the perturbative decay-rate closure breaks down at stronger coupling."],"forward_implications":["At $\\Gamma(k_m)t\\gg 1$ the shock amplitude is exponentially small, but away from the front the profile follows the same $t^{-1}$ ballistic decay as free fermions, so interactions erase the shock by decay rather than by diffusive spreading.","Inside the window $\\lambda t\\Gamma(k_m)\\sim 1$, the low-energy contribution is comparable to the high-energy one, so the profile there is set by the universal function $F$ rather than by the naive ballistic formula.","The region of quantum ripples occupies a fixed fraction $\\sim(\\Delta N)^{-2/3}$ of the shock, so Eq. (38) is the legitimate classical description until $t\\Gamma(k_m)\\sim(\\Delta N)^{2/3}$, by which time the amplitude is exponentially small.","The total density remains a monotonic function of $\\lambda$, so the shock front retains its ordering even while its amplitude decays."],"supporting_citations":[{"why":"Provides the free-fermion shock solution, the Airy-function description of quantum ripples, and the initial-state construction that the interacting calculation starts from.","marker":"6"},{"why":"Supplies the single-particle decay rate $\\Gamma(k)$ and the golden-rule rates $W_{p\\to k}$ used in the linearized collision integral.","marker":"15"},{"why":"Establishes three-particle collisions as the relaxation mechanism for generic integrability-breaking interactions in one dimension.","marker":"14"},{"why":"Provides the kinetic (Boltzmann) equation whose linearization yields the central equation of motion.","marker":"12,13"}],"fun_headline_variants":["Quantum shocks decay, shape persists at late times","Fermionic shocks: amplitude fades, form survives","Shock waves lose punch, not their profile","Late-time shocks: exponential decay, stable shape"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The late-time profile rests on the assumption that the interacting shock is accurately described by the linearized Boltzmann equation with a perturbative single-particle decay rate $\\Gamma(k)$ and with the source $J$ fed only by the overhanging high-energy branch; the authors' check of this closure in Appendix C is perturbative and estimates the neglected contribution as at most about 2%.","fun_headline_variants_meta":{"raw":{"variants":["Quantum shocks decay, shape persists at late times","Fermionic shocks: amplitude fades, form survives","Shock waves lose punch, not their profile","Late-time shocks: exponential decay, stable shape"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000145,"raw_usage":{"total_tokens":1121,"prompt_tokens":832,"completion_tokens":289,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":448,"completion_tokens_details":{"reasoning_tokens":228}},"tokens_in":448,"tokens_out":289,"duration_ms":3915,"temperature":1.0,"reasoning_tokens":228,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:41:38.806721+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A numerical simulation of weakly interacting one-dimensional fermions (or of the full three-particle Boltzmann equation) initialized as a smooth density bump would settle the claim: plot $\\delta\\rho(\\lambda,t)/[(k_m-k_0)e^{-t\\Gamma(k_m)}t_S\\sqrt{\\lambda}/(t\\pi)]$ against $\\lambda t\\Gamma(k_m)$ for $t\\gg t_S$. Collapse onto $1+F(\\lambda t\\Gamma(k_m))$ supports Eq. (38); a systematic deviation, or a low-energy source contribution growing beyond the stated $\\sim 2\\%$, falsifies the linearization and decay-rate closure.","supporting_citations":[{"cited_title":"Bettelheim, L","cited_arxiv_id":null,"evidence_quote":"Provides the free-fermion shock solution, the Airy-function description of quantum ripples, and the initial-state construction that the interacting calculation starts from."},{"cited_title":"Khodas, M","cited_arxiv_id":null,"evidence_quote":"Supplies the single-particle decay rate $\\Gamma(k)$ and the golden-rule rates $W_{p\\to k}$ used in the linearized collision integral."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes three-particle collisions as the relaxation mechanism for generic integrability-breaking interactions in one dimension."}],"review_version":1}