{"id":"c13e9633-ce5a-4797-b6a5-12ac96e9720d","arxiv_id":"2506.17148","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Small perturbations of a large shocking simple wave in a strictly hyperbolic 1D conservation law still form a shock, with a universal inverse-cubic leading profile and a Lipschitz maximal development boundary.","lead":"This math paper proves that shock formation is stable near large 'simple waves' in one-dimensional hyperbolic conservation laws, and describes the exact shape of the singularity. It shows that nearby perturbed solutions still develop shocks, and gives a detailed expansion of the solution near the shock, which will be used to study vanishing viscosity.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The higher-regularity bootstrap in Section 7.5 is omitted, but Theorem 4.2(iii) and the homogeneous expansion Theorem 8.10 depend on it: the C^k-in-(τ,u) bound is not derived from C^1-small data.","rationale":"The reader's weakest-assumption pick, Condition 4.1, identifies a genuine hypothesis on the background wave, and the reader also flags Sections 7.4 and 7.5. I agree that Condition 4.1 is a real restriction, but I do not see it as the most load-bearing point: the paper states it clearly and the construction in Section 4.1 shows admissible examples. The more pressing issue is the omitted higher-regularity induction in Section 7.5. Theorem 4.2(iii) is used directly in Section 8 to pass from smoothness in (τ,u) to the homogeneous expansion of ψ, and it is exactly the step that needs to convert C^1-small data plus C^{k+1} data into a C^k bound on P. The preceding bootstrap arguments assume full high-order smallness, so they cannot by themselves prove the stated estimate. Corollary 7.9 is the bridge, and its proof is explicitly deferred. The constant's dependence on ||P||_{C^{k-1}} is conventional in nonlinear hyperbolic bootstrap arguments, so the gap is likely repairable; hence I do not move the verdict to reject. But until the induction is supplied, Theorem 8.10, and with it the companion vanishing-viscosity application, rests on an unproved regularity propagation statement. A concrete re-derivation of the induction for k = 3 would settle whether the current hypotheses are sufficient or whether an extra smallness condition on high derivatives is needed.","tokens_in":42536,"tokens_out":16854,"duration_ms":169208,"concrete_test":"Write out the missing induction in Corollary 7.9 for k = 3: commute Proposition 5.2 up to Γ^3, apply Propositions 6.1 and 6.2 to the equations for Γ^α φ_{I'}, then integrate the equations for Γ^α Φ, Γ^α ψ, and Γ^α μ. Verify at each step that every Q_|α| term can be bounded using ||P||_{C^2_{τ,u}} and ε, with no occurrence of ||Γ^3 P|| on the right-hand side before the estimate closes. If the induction succeeds, the gap is filled; if closing requires an a priori bound on ||P||_{C^3_{τ,u}}, then Theorem 8.10 needs an additional smallness assumption on third derivatives of the initial perturbation v.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The single most load-bearing gap is Section 7.5. Theorem 4.2(iii) asserts ||P||_{C^k_{τ,u}(M_{t*+δ})} ≤ M ||ψ(0,·)||_{C^{k+1}_x} from an initial perturbation that is only C^1-small, and Section 8 uses exactly this statement to obtain Taylor expansions for ψ and μ in (τ,u) up to order k. But Sections 7.2–7.3 close a bootstrap in which the full C^k norm of P is already assumed small (for example ||P||_{C^k_{τ,u}} ≤ ε^{3/4} or ≤ ζ^{3/4}); they do not show that large higher derivatives of the data remain controlled. Corollary 7.9 is the step that would bridge this, but its proof is explicitly omitted: 'The propagation of higher regularity estimates after closing a nonlinear problem has been carried through in several contexts, so we omit the details.' The corollary's constant is allowed to depend on ||P||_{C^{k-1}}, which is a plausible structure, but the commuted equations in Proposition 5.2 contain quadratic terms Q_|α| involving derivatives of the same order, so it is not established that the C^k norm does not feed back into itself. Without this induction, Theorem 8.10's expansion to homogeneity k−2 is unsupported for data that are only C^1-small. This is an admitted omission rather than a demonstrated contradiction, but it is load-bearing for the paper's strongest advertised conclusion.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the 1D hyperbolic system (3.1) under strict hyperbolicity and genuine nonlinearity of one characteristic family. It proves that any sufficiently C^1-small perturbation of a shocking simple wave (with a mild nondegeneracy condition for intermediate speeds, Condition 4.1) forms a shock at time t* = 1 + O(ε), that the solution is C^k regular in the eikonal coordinates (τ,u) up to the boundary of a 'boxed' maximal globally hyperbolic development, and that that boundary is Lipschitz and decomposes into preshock, Cauchy horizon, singular, and extensible sets (Theorem 4.2). Under a stronger nondegeneracy condition (Condition 8.1) it gives a homogeneous expansion of ψ about the preshock with a universal leading inverse-cubic cusp p1 = −a0 ∂_{I0}λ(0)^{-1}𝔲 r_{I0}(0) (Theorem 8.10). The proof is based on renormalized unknowns P, eikonal coordinates, the structural cancellation in (5.3), and hyperbolic estimates in Sections 6–7.","tokens_in":42922,"tokens_out":5311,"duration_ms":52997,"significance":"If the results hold, this is a substantial contribution: it extends perturbative shock formation from small data to a large class of simple waves, handles intermediate characteristics under a mild hypothesis, and gives one of the few detailed descriptions of the MGHD boundary in a 1D system. The paper is transparent about its hypotheses, and the renormalization/cancellation in (5.3) is clean. The homogeneous expansion and its universal leading term are well motivated and ready for use in the companion vanishing-viscosity paper. Condition 4.1 is a genuine structural assumption on the background rather than an ad hoc fit, and the parameter dependencies are, for the most part, explicitly tracked.","major_comments":[{"comment":"The advertised conclusion Theorem 4.2(iii) asserts ||P||_{C^k_{τ,u}(M_{t*+δ})} ≤ M ||ψ(0,·)||_{C^{k+1}_x} from an initial perturbation that is only C^1-small. The proof given in Sections 7.2–7.3 is a bootstrap that assumes the C^k norm is already small (e.g. by ζ^{3/4} in (7.6)) and improves that assumption; it does not show that C^k smallness is propagated from C^{k+1}-bounded data. Corollary 7.9 is exactly the missing step, but its proof is omitted with the sentence 'The propagation of higher regularity estimates after closing a nonlinear problem has been carried through in several contexts, so we omit the details.' This is load-bearing: Theorem 8.10 and Proposition 8.2 use Theorem 4.2(iii) to Taylor-expand μ and ψ in (τ,u), so without a complete induction the homogeneous expansion is not justified at the stated level of regularity. Please supply the full induction or state Theorem 4.2(iii) under an explicit C^k-smallness assumption.","section":"Section 7.5, Corollary 7.9"},{"comment":"The proof of the actual zero-crossing of μ is only sketched. The text says 'This argument has appeared several times in the literature (see, for example, [55]), so we only sketch it.' The sketch asserts ∂τ μ = −1 + O(ε) near the Θ-shocking characteristic and integrates to μ ↘ 0 at t* = 1 + O(ε). To justify Theorem 4.2(i) and the construction of the MGHD one must prove that the relevant characteristic remains in the domain D*_0 where the C^k bounds have been established, that the O(ε) term is uniform as μ → 0, and that inf |(∂_u ψ)^{I0}| is bounded below on the zero set; the current sketch does not track these constants or the size of the neighborhood. Since the formation of a shock is the central assertion, this zero-crossing should be proved in full rather than referenced.","section":"Section 7.4"},{"comment":"The characterization of the boundary B_{1+δ} as eΣ*_0 and the decomposition into B^{pre}, B^{Cau}, B^{sing}, B^{ext} is asserted after a relatively short argument. In particular, the proof that B^{sing} is empty for intermediate characteristics uses (7.11) and states that following an integral curve of μL^{(N)} from q lowers μ to zero 'within time ≍ r'; this quantitative claim is not derived, and the choice of r and the uniform positivity of μL^{(N)}μ need to be established. Since the Lipschitz regularity of B_{t*+δ} and the component classification are part of Theorem 4.2(ii), this section should be expanded into a complete proof.","section":"Section 7.6"}],"minor_comments":[{"comment":"The heading 'Renormalized eqations' should read 'Renormalized equations'.","section":"Section 5 heading"},{"comment":"The sentence 'We therefore satisfy the required L1 smallness (α < 1)' should refer to the quantity κ from Proposition 6.2, not α.","section":"Section 6, after Proposition 6.2"},{"comment":"The expression '|τ − 1_e| ≤ δ' appears to contain a typo and should probably read '|τ − 1| ≤ δ'.","section":"Section 7.4"},{"comment":"The error exponent in Proposition 8.9 is stated as O(𝔡^{h−2k−3n+1}), whereas Theorem 8.10 has O(𝔡^{h−2m−3n+1}); if this is not a typo, the dependence on m and n should be clarified, since the two statements otherwise look parallel.","section":"Propositions 8.9 and 8.10"}],"recommendation":"major_revision","confidential_remarks":"The main obstacle to acceptance is the omitted proof of Corollary 7.9. The stress-test concern is accurate: Theorem 4.2(iii) currently over-claims relative to the proof, and Theorem 8.10 inherits the gap. The zero-crossing proof in Section 7.4 and the boundary characterization in Section 7.6 also need to be written out in full. These issues are fixable within the manuscript's scope—they are omissions rather than demonstrated contradictions—and the renormalized structure in Section 5 is sound, so I would not reject on current evidence."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the thing to know: this is a real step forward, not a routine extension. The authors show shock formation is stable near large-amplitude simple waves in general strictly hyperbolic 1D systems, including intermediate speeds under a mild nondegeneracy condition, and they give a homogeneous expansion with a universal cubic cusp. The nice surprise is the intermediate-wave instability mechanism in Section 5.2: they identify why the naive Grönwall argument fails and resolve it by confining degeneracy to a small interval. The renormalized unknowns and the cancellation in (5.3) are clean and convincing. The MGHD boundary decomposition into preshock/Cauchy horizon/singular parts is also new for this setting and carefully stated.\n\nThe soft spot is load-bearing. Theorem 4.2(iii) asserts C^k regularity in (τ,u) from C^1-small data, but the proof is delegated to Corollary 7.9, whose proof is explicitly omitted. The bootstrap in Sections 7.2–7.3 assumes the full C^k norm is already small, so it does not control large higher derivatives. And the commuted equations in Proposition 5.2 contain quadratic terms at top order, so the step is not innocent. Theorem 8.10 uses exactly this regularity to get Taylor expansions up to order k−2. As written, the advertised homogeneous expansion is unsupported for data that are only C^1-small. This is an admitted omission, and it may be fillable by standard regularization arguments, but it is not a detail.\n\nThe zero-crossing of μ is only sketched in 7.4, but that is minor; the mechanism is standard once C^1 control is available. The MGHD uniqueness caveat is open and stated honestly.\n\nWho the paper is for: people working on shock formation in conservation laws and hyperbolic systems. They will find the geometric framework and the intermediate-wave analysis valuable. My own verdict is conditional: the first-order stability theory is likely correct, and the expansion is plausible, but the proof as written does not close the advertised regularity.\n\nRecommendation: yes, send to a serious referee — the ideas are substantial and the authors are clearly honest about what they have and haven't shown. Insist that the higher-regularity bootstrap be supplied, or that Theorem 4.2(iii) and Theorem 8.10 be restated with the regularity assumed small.","headline":"A genuinely new and careful paper on shock formation for large data, but the advertised higher-regularity bootstrap is omitted, so the strongest theorems are conditional pending that proof.","tokens_in":43433,"tokens_out":2755,"would_cite":true,"duration_ms":28912,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35L65","35L67","35B44"],"pacs":[],"model":"deepseek-v4-flash","headline":"Shock formation is stable near shocking simple waves in 1D genuinely nonlinear systems, and the first singularity has a universal inverse-cubic cusp.","keywords":["shock formation","conservation laws","simple waves","eikonal function","inverse foliation density","maximal globally hyperbolic development","homogeneous expansion","nondegenerate singularity"],"falsifier":"Compute the spectrum of the linearized transverse operator in the idealized case $\\mu=1-\\tau$: if any eigenvalue $\\gamma$ is negative, perturbation amplitudes scale like $(1-\\tau)^\\gamma$ and diverge at the shock time. A system satisfying the structural hypotheses with an intermediate-speed simple wave that violates Condition 4.1 and exhibits such a negative eigenvalue would refute unconditional stability.","tokens_in":42361,"feed_emoji":"💥","tokens_out":13170,"duration_ms":128380,"temperature":0.7,"pith_summary":"One-dimensional conservation laws with a genuinely nonlinear wave family—meaning the relevant wave speed genuinely changes as the wave compresses—often turn smooth data into shocks. This paper proves that near any background simple wave that shocks, every sufficiently small $C^1$ perturbation shocks too, at a time differing by $O(\\varepsilon)$, and remains smooth in specially adapted coordinates up to that time. The background may have arbitrarily large amplitude, so the result goes beyond the classical small-data regime; the only extra assumption is a mild confinement condition when the shocking wave has intermediate speed. The paper also describes the boundary of the maximal classical development near the first singularity and, under a generic nondegeneracy, expands the solution into homogeneous functions whose leading term is a universal inverse-cubic cusp in the shocking component. This description is the input used in the companion analysis of the vanishing viscosity limit.","feed_headline":"Shock formation survives C^1 perturbations of simple waves","feed_subtitle":"A geometric eikonal proof shows the first singularity is a universal inverse-cubic cusp, opening large-data shock theory.","key_machinery":"The machinery is the eikonal description of the characteristic family. One defines $u$ to be constant along the genuinely nonlinear $I_0$-characteristics and sets $\\mu=(\\partial_x u)^{-1}$, the inverse foliation density; shock formation is exactly the vanishing $\\mu\\to 0$. Changing to $(\\tau,u)$ coordinates sends $\\partial_\\tau$ to $L=\\partial_t+\\lambda\\partial_x$ and $\\partial_u$ to $\\mu\\partial_x$, so the dangerous spatial derivative is tamed by $\\mu$. The fundamental unknowns $\\mathcal{P}$ satisfy equations in which the quadratic self-interaction term that drives ordinary blowup cancels, leaving $\\Phi=\\mu(\\partial_x\\psi)^{I_0}$ bounded in $(\\tau,u)$. For a simple wave, $\\mu$ decreases linearly to zero along the shocking characteristic, and the perturbative analysis shows the second $\\tau$-derivative of $\\mu$ stays small, forcing $\\mu$ to hit zero at time $1+O(\\varepsilon)$. For intermediate-speed shocks, Condition 4.1 confines the background's near-degeneracy to a short interval because otherwise the linearized transverse system has no preferred integration direction.","core_discovery":"The paper's central claim is that shock formation is stable near shocking simple waves: fixing a background simple wave $\\Theta$ that first shocks at time $1$ and satisfies the structural hypotheses, any initial datum $\\psi(0,\\cdot)=\\Theta(0,\\cdot)+v$ with $v$ small in $C^1$ produces a solution whose first blowup is a shock at time $t_*=1+O(\\varepsilon)$. In the eikonal coordinates $(\\tau,u)$ adapted to the shocking characteristic, the fundamental unknowns $\\mathcal{P}$—including $\\psi$, the renormalized derivative $\\Phi=\\mu(\\partial_x\\psi)^{I_0}$, the transverse derivatives, and the inverse foliation density $\\mu=(\\partial_x u)^{-1}$—are as smooth as the data, so the physical gradient blowup is entirely a degeneracy of the coordinate map. With a stronger generic nondegeneracy, the solution admits a homogeneous expansion whose leading singular term is $p_1=-a_0\\,\\partial_{I_0}\\lambda(0)^{-1}\\,\\mathfrak{u}\\,r_{I_0}(0)$, where $\\mathfrak{u}$ solves $x=-a_0 t\\,\\mathfrak{u}+b_0\\mathfrak{u}^3$. The same estimates characterize the maximal globally hyperbolic development in a spacetime box: its future boundary is uniformly Lipschitz and splits into preshock, Cauchy horizon, singular, and extensible sets, with the singular set empty for intermediate-speed shocks.","pith_inferences":["If, as the paper suggests, typical first singularities are isolated nondegenerate shocks, then shortly before blowup any such solution approximates a simple wave, and this perturbative theorem would cover generic large-data shock formation rather than only an open neighborhood of simple waves.","The universal inverse-cubic cusp indicates that the strong-norm vanishing viscosity profile should be the same for all such systems, not just for scalar conservation laws; the companion paper is the natural test of this transfer.","The linearized divergence described in Section 5.2 suggests the mild nondegeneracy condition is not cosmetic: constructing an explicit intermediate-speed system that violates Condition 4.1 with a negative linearized eigenvalue would show unconditional stability really fails outside the theorem's hypotheses."],"forward_implications":["Every $C^1$-small perturbation of a shocking simple wave forms a shock at time $1+O(\\varepsilon)$, so the set of shocking initial data is open around large-amplitude simple-wave data.","In eikonal coordinates the solution is as regular as its data up to the shock, so the singularity is carried entirely by the coordinate change rather than by the renormalized solution variables.","The boundary of the maximal development in a spacetime box is uniformly Lipschitz, and generically it is a crease of two curves with a point mass of curvature; the Lipschitz bound is sharp, since smooth data can produce boundary curvature that is not a Radon measure.","Under the generic nondegeneracy condition, the leading singularity is universal: after a shear-and-shift change of frame, the shocking component behaves like the inverse-cubic profile $\\mathfrak{u}$ with coefficients set by the eigenvalue geometry.","For intermediate-speed shocks, the singular set on the boundary of the maximal development is empty; the boundary consists of a preshock together with a Cauchy horizon."],"supporting_citations":[{"why":"Establishes the small-data singularity framework for one-dimensional genuinely nonlinear systems that this paper extends to large simple-wave data.","marker":"[45]"},{"why":"Introduces genuine nonlinearity as the mechanism that drives shock formation.","marker":"[50]"},{"why":"Supplies the density-of-characteristics idea underlying the inverse foliation density and eikonal coordinates.","marker":"[48]"},{"why":"Provides the geometric eikonal method in which blowup is treated as degeneracy of the coordinate map.","marker":"[29]"},{"why":"Develops the homogeneous expansion in the universal cusp profile that Section 8 adapts from scalar equations to systems.","marker":"[25]"},{"why":"Supplies the maximal-development-in-a-box construction and the boundary decomposition used in Theorem 4.2(ii).","marker":"[70]"},{"why":"Describes the one-dimensional MGHD boundary geometry, including the curvature crease that Proposition 4.3 shows is sharp.","marker":"[2]"},{"why":"Provides a nearby one-dimensional shock-formation setting for simple waves that the present paper generalizes to strictly hyperbolic systems.","marker":"[33]"}],"fun_headline_variants":["Stable shock formation: inverse-cubic cusp is universal","Perturbations of simple waves still shock: universal cusp","Shock formation stable; singularity follows inverse-cubic law","Universal inverse-cubic cusp in stable shock formation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"For intermediate-speed shocks, the proof assumes the background wave's flattening region is confined to a short interval; without this confinement, small perturbations can diverge at the shock time and stability is not known.","fun_headline_variants_meta":{"raw":{"variants":["Stable shock formation: inverse-cubic cusp is universal","Perturbations of simple waves still shock: universal cusp","Shock formation stable; singularity follows inverse-cubic law","Universal inverse-cubic cusp in stable shock formation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000442,"raw_usage":{"total_tokens":2225,"prompt_tokens":918,"completion_tokens":1307,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":534,"completion_tokens_details":{"reasoning_tokens":1238}},"tokens_in":534,"tokens_out":1307,"duration_ms":12032,"temperature":1.0,"reasoning_tokens":1238,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:10:41.239024+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the spectrum of the linearized transverse operator in the idealized case $\\mu=1-\\tau$: if any eigenvalue $\\gamma$ is negative, perturbation amplitudes scale like $(1-\\tau)^\\gamma$ and diverge at the shock time. A system satisfying the structural hypotheses with an intermediate-speed simple wave that violates Condition 4.1 and exhibits such a negative eigenvalue would refute unconditional stability.","supporting_citations":[{"cited_title":"Formation of singularities in one-dimensional nonlinear wave propagation.Comm","cited_arxiv_id":null,"evidence_quote":"Establishes the small-data singularity framework for one-dimensional genuinely nonlinear systems that this paper extends to large simple-wave data."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces genuine nonlinearity as the mechanism that drives shock formation."},{"cited_title":"Keller, Lu Ting","cited_arxiv_id":null,"evidence_quote":"Supplies the density-of-characteristics idea underlying the inverse foliation density and eikonal coordinates."},{"cited_title":"The formation of shocks in 3-dimensional fluids","cited_arxiv_id":null,"evidence_quote":"Provides the geometric eikonal method in which blowup is treated as degeneracy of the coordinate map."},{"cited_title":"The inviscid limit of viscous Burgers at nondegenerate shock formation","cited_arxiv_id":null,"evidence_quote":"Develops the homogeneous expansion in the universal cusp profile that Section 8 adapts from scalar equations to systems."},{"cited_title":"The geometry of maximal development and shock formation for the Euler equations in multiple space dimensions","cited_arxiv_id":null,"evidence_quote":"Supplies the maximal-development-in-a-box construction and the boundary decomposition used in Theorem 4.2(ii)."},{"cited_title":"The relativistic Euler equations: ESI notes on their geo- analytic structures and implications for shocks in1𝐷 and multi-dimensions","cited_arxiv_id":null,"evidence_quote":"Describes the one-dimensional MGHD boundary geometry, including the curvature crease that Proposition 4.3 shows is sharp."},{"cited_title":"On the formation of shocks of electromagnetic plane waves in non-linear crystals","cited_arxiv_id":null,"evidence_quote":"Provides a nearby one-dimensional shock-formation setting for simple waves that the present paper generalizes to strictly hyperbolic systems."}],"review_version":2}