{"id":"138e0962-3de2-4291-9618-974f53e3b993","arxiv_id":"2509.02965","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For scalar viscous conservation laws with polynomial flux u^p, p in [2,4], arbitrarily large viscous shock profiles are L2-contractive and time-asymptotically stable under small H1 perturbations, with t^{-1/4} decay for L1 perturbations.","lead":"This math paper proves that for a scalar viscous conservation law with flux u^p (2≤p≤4), even very strong shock waves are stable: any small smooth perturbation shrinks to zero over time in a weighted L2 sense, and vanishes with a t^{-1/4} decay rate if the perturbation is integrable. It answers an open question about whether the a-contraction method can work for large shocks.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Lemma 3.2 as written fails at p=2: the asserted l1(u_-,u_+)<0 is 0, so the proof of the uniform dissipation constant β is incomplete; the lemma may be repairable but needs verification.","rationale":"The paper's central claim is an L^2-contraction and asymptotic stability result for arbitrarily large viscous shocks with f(u)=u^p, 2≤p≤4. The proof is a weighted relative-entropy argument, and the whole contraction estimate rests on Lemma 3.2: the weight-dependent coefficient g(U) must be uniformly negative on [u_+,u_-]. The reader identified exactly this as the weakest assumption, and my independent reading agrees. The specific failure at p=2 is real: the inequality l1<0 asserted in (3.18) is false at p=2, since both terms vanish. Therefore, as written, the proof of Lemma 3.2 does not establish the strict negativity at u_- for p=2. This is a genuine proof gap, not merely a stylistic issue, because the uniform β is what makes (3.28) dissipative. However, the gap is localized and appears fixable: for p=2 the full expression for g(u_-) is -a(u_-)<0, and for p∈(2,4) the omitted l1 inequality can be verified by a short calculus argument. I did not find a more fundamental obstruction: the weighted Poincaré inequality, the shift ODE, the smallness absorption of the cubic terms, and the time-asymptotic decay argument all appear structurally sound once Lemma 3.2 is repaired. For this reason, the concern does not move the verdict from the reader's CONDITIONAL assessment; it reinforces it. A single concrete algebraic re-check of Lemma 3.2 with the full expression would settle whether the concern is merely a fixable omission or an actual failure of the method.","tokens_in":20969,"tokens_out":19774,"duration_ms":200773,"concrete_test":"Compute g(u_-) from the full expression (3.14)-(3.16) without dropping the quadratic term, for p=2 and arbitrary 0<u_+<u_-; verify that g(u_-) = -a(u_-) < 0. More generally, symbolically verify that the full g(U) is strictly negative on [u_+,u_-] for all p∈[2,4] (or find a counterexample), and check that the cubic term in (3.26) is bounded by C∥ϕ∥_{H^1}∫ϕ²|Uξ| with a constant independent of the shock. If these checks pass, Lemma 3.2 and the contraction estimate are repairable with the same final theorem.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing step is Lemma 3.2's claim that the effective diffusion coefficient satisfies g(U) ≤ -β on the whole interval [u_+,u_-]. The proof at the boundary u_- in (3.18) discards the nonpositive quadratic term (p^2-5p+4)/4 and then asserts the bound l1(u_-,u_+) < 0 for all p ∈ [2,4]. This assertion is false at p=2: l1 = 0 exactly. Hence the displayed chain gives only g(u_-) ≤ 0, not the strict negativity needed to define β. The uniform β is what converts the weighted Poincaré estimate (3.12) into the dissipative bound (3.28); without it the term β∫(ϕ^X)^2|Uξ| cannot be used to control the right-hand sides in Proposition 3.1, and the proof of Theorem 1.1 fails as written. The gap appears repairable: retaining the quadratic term for p=2 gives g(u_-) = -a(u_-) < 0, and for p ∈ (2,4) the inequality l1<0 can be proved directly, but that repair is not present in the manuscript. Relatedly, the passage from (3.26) to (3.28) silently drops the cubic remainder ∫O(1)(ϕ^X)^3Uξ dξ; this can also be absorbed by the H^1 smallness, but the absorption should be displayed. Both issues sit in the same key estimate, so a careful re-derivation of Lemma 3.2 and of (3.26)-(3.28) is the right focal point.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the scalar viscous conservation law u_t+(u^p)_x=u_{xx} with 2≤p≤4 and end states 0<u_+<u_-. It constructs an explicit time-dependent shift X(t) and a shock-strength-dependent weight a(U) so that, for H^1-small perturbations of any (possibly arbitrarily large) viscous shock, the weighted L^2 norm of the perturbation is non-increasing in time. From this contraction property the authors derive L∞ convergence of the solution to the shifted shock profile, and, when the initial perturbation is in L1, a t^{-1/4} L2 decay rate. The proof follows the a-contraction framework: a weighted Poincaré inequality on the shock layer, a time-dependent shift chosen to remove a nonlocal term, and a priori H1 estimates. The central algebraic step is Lemma 3.2, which asserts that a certain effective diffusion coefficient g(U) is uniformly negative on [u_+,u_-], yielding the dissipative term in the weighted energy identity.","tokens_in":21370,"tokens_out":18668,"duration_ms":197182,"significance":"If the identified gaps are repaired, this is a meaningful advance: it extends the a-contraction approach for large viscous shocks from the Burgers case to the family of polynomial fluxes f(u)=u^p for 2≤p≤4, thereby addressing the open question raised by Blochas-Cheng about possible viscous destabilization for large shocks. The construction of the shift and weight is explicit and parameter-free, and the stability/decay statements are concrete and falsifiable. The paper uses prior results (L1-contraction, the weighted Poincaré inequality) as black boxes rather than fitting anything to the target conclusion, which is a genuine strength. The main concern is not the overall strategy but the correctness of several displayed estimates inside Lemma 3.2 and the passage from the weighted energy inequality to the final contraction estimate.","major_comments":[{"comment":"The assertion l1(u_-,u_+)<0 for all p∈[2,4] is false at p=2. Direct substitution gives l1=(2-p)u_+(u_- -u_+)u_-^{p-2}+u_+^2(u_-^{p-2}-u_+^{p-2})=0 when p=2. Hence the displayed chain in (3.18) yields only g(u_-)≤0 at p=2, not the strict negativity needed to define β. The uniform negativity of g is exactly what converts the weighted Poincaré estimate into the dissipative bound (3.28); without it the contraction proof fails as written. The gap is repairable—for p=2 one can retain the discarded quadratic term, which equals -1/2 in the bracket and gives g(u_-)≤-a(u_-)<0, or treat p=2 separately using the known Burgers result—but this repair is not present in the manuscript.","section":"Lemma 3.2, Eq. (3.18)"},{"comment":"In passing from (3.26) to (3.28), the cubic remainder ∫ O(1)(φ^X)^3 Uξ dξ is dropped without justification. The estimate (3.27) controls only the shift term 1/2 Ẋ∫(φ^X)^2 aξ, not the cubic term. The cubic term can be absorbed by the H^1 smallness: it is bounded by C∥φ^X∥_{L∞}∫(φ^X)^2|Uξ| ≤ C N(T)∫(φ^X)^2|Uξ|, and this contribution should be included in the coefficient β−Cϵ1. Since (3.28) is the core contraction inequality, the absorption of both the shift term and the cubic remainder should be displayed explicitly.","section":"Section 3, Eqs. (3.26)-(3.28)"}],"minor_comments":[{"comment":"The statement h'''(Ubar)>0 for all p∈[2,4] is false at p=2, where h'''=0. The inequality in (3.21) remains valid with h'''≥0; please adjust the wording to distinguish the equality case.","section":"Eq. (3.21)"},{"comment":"The inequality d/dt ∫(Ca+1)(φ^X)^2 ≤ -2∫(φξ^X)^2 is asserted without proof. It can be obtained by taking C large and combining (3.28) with (3.34), but the argument should be shown, especially the choice of the large constant that dominates the positive right-hand side of (3.34).","section":"Section 2.5, Eq. (2.28)"},{"comment":"The citation 'Kang [20]' in the introduction does not match the bibliography, where [20] is Kružkov; the relevant reference appears to be [10]. Please correct the citation/numbering.","section":"Introduction and references"},{"comment":"The second term on the right-hand side should be C0∥φ0∥²_{H1} to be consistent with (2.23). If the weaker form C0∥φ0∥_{H1} is intended, this should be stated explicitly.","section":"Eq. (2.24)"},{"comment":"For non-integer p, f(u)=u^p is not a polynomial on R and is not smooth at u=0. The local existence statement for arbitrary M>0 needs a smallness/positivity hypothesis, or p should be specified as an integer. Since the continuity argument only uses small M, this does not affect the main result, but the proposition as stated is overbroad.","section":"Proposition 2.1 and Appendix"},{"comment":"There are several minor typos: 'intepolation' for 'interpolation', 'satisfing' for 'satisfying', an incomplete sentence at the end of Section 3 ('complete the proof of Proposition 2.2 by .'), and the symbol β is used for two different constants in Lemma 3.2 and Section 2.5.","section":"General presentation"}],"recommendation":"major_revision","confidential_remarks":"The two major issues identified are local and repairable: Lemma 3.2 can be fixed by handling p=2 separately or retaining the discarded term, and the cubic term in (3.26)-(3.28) can be absorbed by the existing smallness assumption. Once those are repaired and the derivation of (2.28) is expanded, the paper should be publishable. The false inequality at p=2 and the dropped cubic term appear to be genuine but fixable gaps rather than signs of a flawed strategy."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a genuine advance. It is the first L2-contraction result for arbitrarily large viscous shocks in scalar conservation laws with polynomial flux u^p for 2≤p≤4, and it directly addresses the open question raised by Blochas-Cheng. The construction of the weight and the time-dependent shift is explicit, the estimates are derived from the equation, and the local well-posedness argument in the appendix is substantial. No fitting, no circularity; the result stands on the a-contraction framework plus external tools like the weighted Poincaré inequality and L1-contraction. The t^{-1/4} decay under extra L1 initial data is a nice bonus.\n\nThe soft spots are real but not fatal. The stress-test note is essentially correct: in Lemma 3.2, the assertion l1(u_-,u_+)<0 in (3.18) is false at p=2—the expression is exactly zero there. So the displayed proof does not establish the strict negative bound on g(u_-) at p=2, and hence does not yet give the uniform constant β that powers the contraction estimate. The gap is repairable: if the discarded quadratic term is kept, g(u_-)<0 follows immediately for p=2, and for p∈(2,4) the inequality as written can be checked. But the repair is not in the manuscript.\n\nThere are smaller blemishes in the same region. In (3.21), the claim h'''(U)>0 is false at p=2 because h''' is identically zero; the subsequent inequality still holds, so this is minor. The H1 norm in (2.24) appears to be missing a square. The cubic remainder in (3.26) is dropped without explicit absorption before reaching (3.28); it can be absorbed by the H1 smallness, but the step should be displayed. These are all fixable with careful rewriting, not signs of a wrong central argument.\n\nThe citation pattern looks honest: the paper builds on Kang-Vasseur, Kang, and the recent Blochas-Cheng counterexample, and does not overclaim. The reference numbering has small inconsistencies, but nothing misleading.\n\nWho should read it: specialists in scalar conservation laws and the a-contraction method. It settles a concrete open case and shows the method survives large shock strength for a natural class of fluxes. I would send it to peer review rather than desk reject; the referee report should focus on Lemma 3.2 and the cubic-term absorption. With those repaired, the paper is publishable.","headline":"This paper proves L2 contraction and asymptotic stability for arbitrarily large viscous shocks for polynomial fluxes u^p, 2≤p≤4, and the main result is new and credible, but the proof of the key dissipation lemma has a repairable gap at p=2.","tokens_in":21868,"tokens_out":2923,"would_cite":true,"duration_ms":35220,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35L65","35B35","35B40","35L67"],"pacs":[],"model":"deepseek-v4-flash","headline":"For the scalar viscous conservation law u_t+(u^p)_x=u_xx with 2≤p≤4, every viscous shock of arbitrary amplitude is L2-contracting and asymptotically stable under small H1 perturbations, and L1 perturbations decay at the rate t^{-1/4}.","keywords":["L2 contraction","viscous shock wave","large-amplitude shock","a-contraction method","time-dependent shift","scalar viscous conservation law","polynomial flux","asymptotic stability"],"falsifier":"For fixed p∈[2,4] and end states 0<u_+<u_-, evaluate m(U) in (3.16) on a dense grid of U∈[u_+,u_-]. If m(U)≤0 at any interior U, then g(U)≥0 and the dissipation term in (3.28) loses its sign, so the weighted contraction estimate cannot hold as stated; a concrete numerical search would settle the lemma's validity, especially at the endpoint p=2.","tokens_in":20843,"feed_emoji":"📉","tokens_out":11430,"duration_ms":108352,"temperature":0.7,"pith_summary":"The paper settles a question left open by a recent counterexample: whether the a-contraction property, known for arbitrary shock amplitudes at the inviscid level and for small viscous shocks, also holds for large viscous shocks. For the strictly convex polynomial flux f(u)=u^p with 2≤p≤4, the authors prove that for any end states 0<u_+<u_- and any sufficiently small H1 perturbation of the shock profile, the weighted L2 distance to a suitably shifted profile is non-increasing in time. Thus the shock is stable no matter how large its amplitude, provided the initial perturbation is small. The same argument gives L∞ convergence to the shifted profile and, when the initial perturbation is integrable, an explicit L2 decay rate of order t^{-1/4}.","feed_headline":"Arbitrarily large viscous shocks are L2-contracting","feed_subtitle":"The a-contraction method with a shock-strength weight gives stability for u_t + (u^p)_x = u_xx when 2 ≤ p ≤ 4.","key_machinery":"The proof is carried by a weighted relative-entropy functional rather than by the usual anti-derivative or spectral method. The weight is a(U)= (f(U)-f(u_-))/(U-u_-) - (f(U)-f(u_+))/(U-u_+), positive and bounded on the shock layer, and the shift X(t) is chosen by the ODE (2.5) so that the transport term created by shifting becomes a square with favorable sign. The decisive algebraic step is Lemma 3.2, which asserts that the effective diffusion coefficient g(U) in (3.14) is strictly negative on [u_+,u_-]; combined with the weighted Poincaré inequality of Lemma 3.1, this negativity converts the weighted energy identity into pure dissipation plus small cubic terms controlled by the H1 smallness","core_discovery":"The central claim of Theorem 1.1 is that the viscous shock profile U of u_t+(u^p)_x=u_xx is nonlinearly stable in H1 for arbitrary shock strength when p∈[2,4]. More precisely, there is a critical ϵ* so that if the initial H1 distance to the profile is below ϵ*, then a time-dependent shift X(t) and a positive weight a(U) exist with d/dt ∫ a(U(x-st-X(t))) |u(t,x)-U(x-st-X(t))|^2 dx ≤ 0 for all t>0. The shift obeys an ODE that compensates the longitudinal drift of the perturbation; the weight, which depends on the shock strength, is trivial exactly for Burgers flux p=2. As a corollary of the monotonicity and the H1 a priori estimates, the solution converges to the shifted shock in L∞, the shift","pith_inferences":["One endpoint issue in the paper as written: Eq. (3.18) uses the assertion l1(u_-,u_+)<0 for all p∈[2,4], but this quantity is 0 at p=2; the strict negativity needed for dissipation appears recoverable by retaining a quadratic term the proof drops, so the p=2 case likely survives with a small correction.","Because ϵ* and the constants in the estimates depend on u± and p, “arbitrarily large shock” means arbitrary amplitude with small perturbation; extending the contraction to genuinely large perturbations remains open in general.","The same shift-and-weight mechanism may transfer to planar multidimensional viscous shocks or to degenerate Oleinik shocks, with the algebra of the negativity lemma as the main obstacle to overcome."],"forward_implications":["Global existence and uniform H1 bounds hold for small H1 perturbations of arbitrarily large shocks in (1.1) for p∈[2,4].","The weighted L2 distance to the shifted profile is a Lyapunov function, so contraction holds for all time, not only asymptotically.","L∞ convergence of the solution to the shifted shock holds, and the shift velocity tends to zero, so the shift grows at most sublinearly.","If the initial perturbation is in L1, the L2 error decays as C||φ0||/(1+C t^{1/4}||φ0||).","The theorem answers the open question raised by a recent counterexample to a-contraction for large viscous shocks, at least for polynomial fluxes of degree 2 through 4."],"supporting_citations":[{"why":"Establishes L2-contraction for large shocks of the Burgers equation; the time-dependent shift and weight construction used here is modeled on it.","marker":"[13]"},{"why":"Supplies the weighted Poincaré inequality (Lemma 3.1) that converts the weighted energy identity into dissipation.","marker":"[14]"},{"why":"Raises the open question addressed here by showing a-contraction can fail for large viscous shocks of other equations.","marker":"[1]"},{"why":"Provides the weighted-energy strategy for large Oleinik shocks that motivates the proof for polynomial fluxes.","marker":"[7]"},{"why":"Gives the L1-stability of scalar viscous shock waves used to control the shift and to reach the L2 decay rate.","marker":"[27]"},{"why":"Proves a-contraction criteria at the inviscid level for arbitrary shock amplitudes, the property this paper extends to the viscous level for p∈[2,4].","marker":"[12]"}],"fun_headline_variants":["L2 contraction for all shock strengths, p∈[2,4]","Arbitrarily large shocks are L2-contracting","Large shocks stable: L2 contraction proven","Shock stability for arbitrary strength, 2≤p≤4","Viscous shocks: L2 contraction at all amplitudes"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The argument's load-bearing premise is that a certain algebraic expression built from the flux and the shock profile—the effective dissipation coefficient g(U)—is strictly negative on the whole interval [u_+,u_-]; the proof of this lemma as written uses a strict inequality that is actually zero at p=2, though retaining a discarded quadratic term appears to restore the needed negativity.","fun_headline_variants_meta":{"raw":{"variants":["L2 contraction for all shock strengths, p∈[2,4]","Arbitrarily large shocks are L2-contracting","Large shocks stable: L2 contraction proven","Shock stability for arbitrary strength, 2≤p≤4","Viscous shocks: L2 contraction at all amplitudes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000778,"raw_usage":{"total_tokens":3255,"prompt_tokens":705,"completion_tokens":2550,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":449,"completion_tokens_details":{"reasoning_tokens":2466}},"tokens_in":449,"tokens_out":2550,"duration_ms":23416,"temperature":1.0,"reasoning_tokens":2466,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T11:16:13.059073+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For fixed p∈[2,4] and end states 0<u_+<u_-, evaluate m(U) in (3.16) on a dense grid of U∈[u_+,u_-]. If m(U)≤0 at any interior U, then g(U)≥0 and the dissipation term in (3.28) loses its sign, so the weighted contraction estimate cannot hold as stated; a concrete numerical search would settle the lemma's validity, especially at the endpoint p=2.","supporting_citations":[{"cited_title":"Kang and A","cited_arxiv_id":null,"evidence_quote":"Establishes L2-contraction for large shocks of the Burgers equation; the time-dependent shift and weight construction used here is modeled on it."},{"cited_title":"Kang and A","cited_arxiv_id":null,"evidence_quote":"Supplies the weighted Poincaré inequality (Lemma 3.1) that converts the weighted energy identity into dissipation."},{"cited_title":"Viscous Destabilization for Large Shocks of Conservation Laws","cited_arxiv_id":"2501.01537","evidence_quote":"Raises the open question addressed here by showing a-contraction can fail for large viscous shocks of other equations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the weighted-energy strategy for large Oleinik shocks that motivates the proof for polynomial fluxes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the L1-stability of scalar viscous shock waves used to control the shift and to reach the L2 decay rate."},{"cited_title":"Kang and A","cited_arxiv_id":null,"evidence_quote":"Proves a-contraction criteria at the inviscid level for arbitrary shock amplitudes, the property this paper extends to the viscous level for p∈[2,4]."}],"review_version":1}