{"id":"8b4b8ded-47b1-4170-b31c-f3b743fbdbd0","arxiv_id":"2501.16244","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"For the isothermal Euler equations, every bounded-variation entropy solution is unique and L2-stable in the class of weak entropy solutions with strong traces, without any smallness assumption on the variation.","lead":"This paper proves that entropy-compatible BV solutions of the one-dimensional isothermal gas equations are stable against much more irregular perturbations, even when the perturbations are not of bounded variation. It extends a known weak-strong uniqueness principle from small shocks to arbitrarily large shocks by building a carefully weighted front-tracking comparison.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Global a-contraction-with-shifts (Prop. 4.1) is imported but not proved; if its uniform constants fail for large shocks, the weight construction and Theorem 1.1 collapse.","rationale":"The paper's central theorem is a substantial extension of the weak-small-BV stability of [15] to large-BV data for the isothermal p-system. The modified front tracking construction (Sections 5-6) is carefully adapted to the TVD field w = -ln tau, and the weight a is explicitly designed so that the dissipation inequalities hold at each interaction. However, the argument is only as strong as Proposition 4.1, which supplies the shock-level a-contraction inequality. The paper imports this proposition from [30] and [27] and, crucially, the global version with uniform constants is not proved: Remark 4.2 defers the local-to-global compactness step to the author's own preprint [18] in the scalar setting. Since every shock in the front tracking, including large shocks produced by interactions, is assigned a weight ratio C1 or 1/C1 based on the big-shock case of Proposition 4.1, any failure of uniformity (e.g., a ratio that must degenerate with shock strength, or a shift whose existence horizon shrinks to zero) would break the weight monotonicity (32) and the boundary-flux estimates in Section 7, and with them the L2 contraction of Proposition 3.1. I could not find an internal inconsistency in the case-by-case weight construction, and the TVD-based bounds on L(t) and Q(t) (Lemmas 6.1-6.2) are plausible, but those bounds presuppose Proposition 4.1. Thus the appropriate verdict is CONDITIONAL: accept only if the global a-contraction statement is supplied, either by an independent proof or by a detailed extension argument from [30, 27]. This matches the reader's assessment; no verdict change is recommended.","tokens_in":24201,"tokens_out":15975,"duration_ms":143483,"concrete_test":"Verify the local-to-global step in Remark 4.2 for the p-system by computing, for a family of large 1-shocks with tau_L = 1 and tau_R in {2, 10, 100, 1000} (with v_R fixed by the Rankine-Hugoniot condition), the largest C1 < 1 such that the dissipation functional in Proposition 4.1 is non-positive for some shift with speed in [-lambda-hat, 0]. If the required C1 tends to 0, or the shift speed bound lambda-hat blows up, as tau_R grows, the uniform big-shock constants in Proposition 4.1 fail; if C1 stays bounded below by a positive constant independent of tau_R, the deferred compactness argument is credible and the concern does not invalidate the result.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The proof of Proposition 3.1, and hence Theorem 1.1, depends on Proposition 4.1: a global a-contraction-with-shifts estimate for every shock with states in a compact set D, with constants C0, C1, epsilon, and lambda-hat uniform over D. The cited inputs ([30, Prop. 5.1] and [27, Thm. 1]) are local statements around a fixed shock, with [27] specifically for small extremal shocks. Remark 4.2 explicitly defers the local-to-global compactness extension to the author's own unpublished preprint [18], which is in the scalar setting. If that extension fails—for example, if for large shocks the admissible weight ratio a2/a1 must approach 0 or infinity, or if the shift can be constructed only on a time horizon that shrinks to 0 as the shock varies over D—then the fixed ratios C1 and 1/C1 used in the weight xi_i (Section 6, (33)) cannot yield the dissipation inequality F_i^+ - F_i^- <= 0 in Section 7. The weight bounds in Proposition 6.1 and the final L2 contraction would collapse, leaving Theorem 1.1 unsupported. The paper provides no proof of the global statement, so this is the single most load-bearing unverified input.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proves Theorem 1.1: for the one-dimensional isothermal p-system, any BV entropy solution is stable and unique within the class S_weak of weak entropy solutions with strong traces, with no smallness assumption on the BV solution. The proof uses a modified front tracking algorithm together with a-weighted relative entropy with shifts. The weight is constructed so that it decays across every wave interaction and remains bounded away from zero and infinity uniformly in the approximation parameter. The paper also transfers the result to Eulerian coordinates in Theorem 1.2.","tokens_in":24504,"tokens_out":8874,"duration_ms":78257,"significance":"If correct, the result is a significant advance: it extends the weak-small-BV stability principle of Chen–Krupa–Vasseur to large BV data for isothermal gas dynamics, exploiting the TVD field w = -ln(tau). The paper contains a detailed interaction table (Section 5.2), a quantitative weight construction with explicit exponential bounds (Lemmas 6.1--6.3 and Proposition 6.1), and a largely self-contained front-tracking well-posedness proof (Appendix B). These are genuine strengths. The main theorem is falsifiable in the usual mathematical sense and does not rely on data fitting or numerical experiments.","major_comments":[{"comment":"Proposition 4.1 is the only input that supplies a-contraction-with-shifts inequalities for every shock with states in the compact set D, and it is used globally in the Riemann solver (Section 5.1, (19)) and for the shock contribution to the cone estimate in Section 7. The cited results [30, Prop. 5.1] and [27, Thm. 1] are local statements around a fixed shock, and [27] is specifically for small extremal shocks. Remark 4.2 explicitly defers the local-to-global compactness extension to the author's own unpublished preprint [18], which is in the scalar setting. The proof of Theorem 1.1 therefore rests on a global statement that is not proved or stated with hypotheses in the present manuscript. If, for large shocks in D, the admissible weight ratio a2/a1 cannot be kept within the fixed bounds C1 and 1/C1, or the shift can be constructed only on a time horizon that shrinks as the shock varies over D, then the inequalities F_i^+ - F_i^- <= 0 in Section 7 and the weight bounds in Proposition 6.1 would collapse. Please either prove the global version here or replace the reference to [18] with a published global result.","section":"§4, Proposition 4.1 and Remark 4.2"},{"comment":"The verification of the weight monotonicity (32) for interactions of type 4B and 9B is not carried out correctly as written. For type 4B the displayed chain reads 'a(t+,x) = C1 a(t+,x*−) = C1 a(t−,x*−) = C1(1−C0(wm−wl))(1−C0(wr−wm)) a(t+,x) < a(t+,x)', which has the same quantity a(t+,x) on both sides and is therefore circular; the type 9B line similarly mixes t+ and t− on both sides. Since (32) is used in Section 7 together with Lemma 7.1 to pass from t_j^- to t_j^+, the monotonicity at these interactions is a load-bearing step. Please rewrite the computation with the correct left/right states and factors of C1 and 1/C1.","section":"§6, Proposition 6.1"},{"comment":"The chain of inequalities after 'Consider any 0<t<T' contains several apparent notation errors that make the telescoping argument hard to verify: inside the ap-limits the integrand is written with u(t,x) and ψ(t,x) even though the limits are in s at times t_j, and the same symbol is used for the time argument in the integral and the time of the ap-limit (e.g. 'ap lim_{s→t^+_j} ∫ ... η(u(t,x)|ψ(t,x)) dx'). The argument needs to use η(u(s,x)|ψ(t_j,x)) with s→t_j^± and the correct intervals [−R+lt_j, R−lt_j]. Please correct these displayed formulas; the intended telescoping argument appears sound, but as written it is not checkable.","section":"§7, proof of Proposition 3.1"}],"minor_comments":[{"comment":"There are several OCR-style typos in the title and abstract ('ST ABILITY', 'Stro ng Trace', 'entropic'); please clean these up.","section":"Title and abstract"},{"comment":"The notation uses t both for the initial time and for the current time in 'for almost every t ∈ [t, ∞)'; please use e.g. t0 for the initial time.","section":"§4, Proposition 4.1"},{"comment":"The approximating initial data ψ0_δ is defined on R, but the construction by averages on a partition of [−R,R] should specify how the approximation is extended outside [−R,R] so that the BV and L∞ bounds (23) hold on all of R.","section":"§5.3, (23)"},{"comment":"The symbol R is used both for the radius of the cone and for the set of rarefaction indices ('Denote R to be the set of i corresponding to rarefactions'); use a different symbol, e.g. script R, for the set of indices.","section":"§7"},{"comment":"The line 'by taking ǫ ≤ V' should be justified more carefully: since ǫ can be chosen arbitrarily small by Remark 4.3, one can choose ǫ ≤ V before fixing the scheme; please state this explicitly.","section":"Lemma 5.3"},{"comment":"Lemma 7.1 is quoted from [31] without proof; since it is used at every ap-lim step, please state the precise hypotheses (e.g. that u has strong traces) and either give a proof or point to the exact lemma in [31].","section":"§7, Lemma 7.1"},{"comment":"The sentence 'Taking R′ sufficiently large in (16)' should be made quantitative: choosing R′ = R + lT suffices to obtain convergence in L∞([0,T];L2([−R,R])).","section":"Proof of Theorem 1.1"}],"recommendation":"major_revision","confidential_remarks":"The core strategy is coherent and the weight construction is detailed, but the reliance on a global version of Proposition 4.1 that is deferred to the author's own unpublished preprint [18] is the main obstacle. If the author can supply the local-to-global compactness proof (or cite a published global result), and fix the garbled verification of (32) for types 4B and 9B, the paper would be a strong candidate for acceptance. I do not see grounds for rejection, but the current version should not be accepted without those repairs."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a credible extension of Chen–Krupa–Vasseur's small-BV weak-strong stability result to large BV data for isothermal gas dynamics, with the weight construction as the real technical novelty. It deserves serious refereeing, though the referee should press on the deferred global a-contraction statement.\n\nWhat is new: the author uses the TVD field w = -ln(tau) to control how often shocks can change between 'small' and 'big' classifications, then builds a weight a(t,x) whose L(t) and Q(t) components stay uniformly bounded away from zero and infinity. Lemmas 5.2, 5.3, 6.1, and 6.3 form a coherent argument for those bounds, and the final relative-entropy cone estimate in Section 7 is a careful adaptation of the standard tool. Appendix A, converting back to Eulerian coordinates, is a useful addition. The interaction tables are detailed and mostly check out, and the paper is honest about what is imported.\n\nWhere the soft spots are: Proposition 4.1, the global a-contraction-with-shifts estimate for shocks with states in a compact set D, is the single most load-bearing input. The cited results are local around a fixed shock, and Remark 4.2 explicitly defers the local-to-global extension to the author's own unpublished preprint [18], which is in the scalar setting. The paper does not prove that the constants C0, C1, epsilon, and lambda-hat remain uniform over all large shocks in D, or that the shift can be constructed on a time horizon that does not shrink as the shock varies. If that extension fails, the weight ratio bounds in Proposition 6.1 and the dissipation inequality in Section 7 collapse, and Theorem 1.1 with them. On the evidence in front of me, this is a gap in what is proved here, not evidence of a false theorem, but it is exactly where the referee should concentrate.\n\nAlso worth noting: Lemma 7.1 is quoted with a proof omitted, and a few interaction estimates (e.g., type 5B and type 10) are asserted with 'one may show' rather than demonstrated. Those are minor by comparison.\n\nThe paper is for researchers working in conservation laws and weak-strong stability. It extends an established program rather than opening a new one, so the novelty is incremental but substantial for this specific system. The central argument hangs together if the deferred global contraction estimate is true. It deserves a serious referee rather than a desk rejection.","headline":"Extends the Chen-Krupa-Vasseur weak-BV stability framework to large BV data for isothermal Euler; the weight construction is genuine new work, but one load-bearing estimate is deferred to an unpublished preprint.","tokens_in":24994,"tokens_out":1563,"would_cite":true,"duration_ms":17106,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35L65","35B35","35L45","76N15"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proves that every entropic BV solution of the isothermal Euler system is stable and unique inside the much larger class of weak entropy solutions with strong traces, with no smallness assumption on the BV solution.","keywords":["isothermal gas dynamics","p-system","weak-BV stability","front tracking","relative entropy","a-contraction with shifts","total variation diminishing","uniqueness"],"falsifier":"Test Proposition 4.1 directly: fix a shock $(u_l,u_r)$ with both states in the compact set $D$, take $u\\in S_{\\mathrm{weak}}$ to be a nontrivial wild solution, and check whether, for every Lipschitz shift $h$ with $h(0)=x_0$, the weighted flux functional $a_1(\\dot h\\,\\eta(u_-|u_l)-q(u_-;u_l))-a_2(\\dot h\\,\\eta(u_+|u_r)-q(u_+;u_r))$ is positive on a set of positive measure; one shock and one $u$ for which it is always positive disproves the global $a$-contraction estimate and collapses the proof.","tokens_in":24009,"feed_emoji":"","tokens_out":9104,"duration_ms":81274,"temperature":0.7,"pith_summary":"This paper proves that a shock solution of the one-dimensional isothermal Euler equations, once its total variation is finite, is the only entropy-respecting weak solution that can carry the same initial data. The perturbations are allowed to be much wilder than the solution itself: they need not have bounded variation, only strong traces on every Lipschitz curve. The method builds a piecewise-constant 'shadow' solution by a modified front-tracking algorithm and assigns each shock a weight; as long as the weight stays bounded between two positive constants, the relative-entropy inequality forces the shadow and the wild solution to remain close in $L^2$ on every cone of information. The main contribution is showing that this weight can be constructed when the shocks are large, not merely small.","feed_headline":"Large-BV shocks: stable against wild entropy waves","feed_subtitle":"Front-tracking weights extend weak-strong stability to shocks of any finite size, no smallness required.","key_machinery":"The argument is carried by a weight function $a(t,x)=C^{L(t)}Q(t)\\prod_i \\xi_i(t,x)$ attached to a modified front-tracking approximant $\\psi$. The factor $L(t)$ counts the finitely many interactions in which big shocks merge or change classification; $Q(t)$ multiplies the 'small-shock numbers' accumulated when small shocks interact; each $\\xi_i$ encodes, for one shock, the a-contraction-with-shifts coefficient (a weighted relative-entropy contraction with a moving shift): $1\\pm C_0\\sigma_i$ for a small shock of strength $\\sigma_i$, and $C_1$ or $1/C_1$ for a big shock. Proposition 6.1 shows the weight is uniformly bounded above and below, is non-increasing at every wave interaction, and has exactly the jump ratios required by the shock contraction estimates. The companion input is the TVD field $w=-\\ln\\tau$: by Lemma 2.1 its variation $D(w_1,w_3)\\le D(w_1,w_2)+D(w_2,w_3)$ never increases across interactions, which limits how often shocks can switch between 'big' and 'small' labels and keeps $L$ and $Q$ under control.","core_discovery":"The central claim is Theorem 1.1. Let $u$ be an entropy solution of the $p$-system $\\partial_t\\tau-\\partial_x v=0$, $\\partial_t v+\\partial_x(1/\\tau)=0$ with bounded variation, initial data $u_0$, and strong traces; let $\\{u_n\\}$ be any sequence in $S_{\\mathrm{weak}}$, the class of weak entropy solutions with strong traces, whose initial data converge to $u_0$ in $L^2(\\mathbb{R})$. Then $u_n\\to u$ in $L^\\infty([0,T];L^2([-R,R]))$ for every $T,R>0$, and $u$ is the unique solution in $S_{\\mathrm{weak}}$ with data $u_0$. No smallness is assumed on the BV solution: the only regularity demanded of the competitors is the strong-trace property. A parallel statement for the Eulerian isothermal Euler system follows by the standard equivalence of the two coordinate frames.","pith_inferences":["If the deferred local-to-global compactness step in Remark 4.2 goes through, the same weight construction is likely portable to any 2x2 system with a TVD field and the same wave-curve structure; the author supplies isothermal Euler, but the role of the TVD field is general enough to suggest other candidates such as isothermal Euler-Poisson.","A sharper version of Theorem 1.1 might quantify the rate of convergence: the constants in Lemmas 6.1–6.3 and Proposition 6.1 depend explicitly on $\\|u_0\\|_{BV}$, $\\|u_0\\|_{L^\\infty}$, and $\\beta$, so the proof likely yields an explicit modulus of continuity in the initial $L^2$ error, though the paper does not track it.","The strong-trace assumption on the wild solutions is likely close to necessary: without it, the boundary terms in the relative-entropy computation are uncontrolled, so uniqueness could fail for merely $L^\\infty$ entropy solutions; an explicit counterexample of that kind would delineate the true boundary of $S_{\\mathrm{weak}}$.","The classification threshold $\\epsilon$ between 'big' and 'small' shocks is chosen so the two contraction regimes in Proposition 4.1 overlap; testing the overlap numerically for large shocks would provide a cheap check of the whole construction before attempting the deferred compactness proof."],"forward_implications":["Entropic BV solutions of the isothermal p-system are unique inside $S_{\\mathrm{weak}}$: any weak entropy solution with strong traces and the same initial data must coincide with the BV solution.","The stability estimate transfers to Eulerian coordinates (Theorem 1.2), so the same uniqueness holds for the isothermal Euler system in its original variables.","The result holds for arbitrarily large but finite total variation, extending the weak-small-BV stability principle of the prior 2x2 theory to a genuinely large-BV statement for this system.","Because $S_{\\mathrm{weak}}$ contains solutions produced by compensated compactness, BV solutions are stable against that entire family of wild solutions whenever strong traces are available.","The front-tracking approximants are not exact solutions; they serve only as a BV shadow that propagates $L^2$ closeness, so the theorem does not require the approximating scheme itself to converge to the BV solution."],"supporting_citations":[{"why":"Supplies the a-contraction-with-shifts estimate for large/extremal shocks used as Proposition 4.1's big-shock regime.","marker":"[30]"},{"why":"Supplies the sharp a-contraction estimate for small extremal shocks used as Proposition 4.1's small-shock regime.","marker":"[27]"},{"why":"Establishes the modified front-tracking plus weight strategy and the weak-small-BV stability framework that this paper extends to large BV.","marker":"[15]"},{"why":"Provides the TVD inequality for the field w = -ln tau recorded as Lemma 2.1, which controls shock classification changes.","marker":"[37]"},{"why":"Gives the equivalence of Eulerian and Lagrangian weak solutions used to transfer the main theorem to Eulerian coordinates.","marker":"[43]"},{"why":"Part of the L1-theory synthesis in Theorem 2.1 that identifies the limit of the front-tracking shadows with the unique BV solution.","marker":"[19]"},{"why":"Completes Theorem 2.1, the uniqueness of BV entropy solutions, which fixes the target of the stability argument.","marker":"[10]"},{"why":"Provides the relative-entropy approximate-limit lemma (Lemma 7.1) that allows summing dissipation across interaction times.","marker":"[31]"}],"fun_headline_variants":["Large-BV shocks: unique and stable against any weak wave","No smallness condition: large-BV solutions win in gas dynamics","Isothermal Euler: uniqueness holds for all finite BV data","Huge shocks, no smallness: uniqueness proven for Euler","Strong traces guarantee uniqueness for large-BV entropy solutions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole argument hinges on a global contraction estimate, imported from two prior papers, which says that every large shock (both states in a fixed compact set) admits a moving shift making a weighted relative-entropy flux non-positive; the paper uses the known local versions of this estimate and defers the local-to-global compactness step to a separate preprint by the same author, and if that deferred step fails for large shocks the weight construction and dissipation argument collapse.","fun_headline_variants_meta":{"raw":{"variants":["Large-BV shocks: unique and stable against any weak wave","No smallness condition: large-BV solutions win in gas dynamics","Isothermal Euler: uniqueness holds for all finite BV data","Huge shocks, no smallness: uniqueness proven for Euler","Strong traces guarantee uniqueness for large-BV entropy solutions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000284,"raw_usage":{"total_tokens":1623,"prompt_tokens":840,"completion_tokens":783,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":456,"completion_tokens_details":{"reasoning_tokens":699}},"tokens_in":456,"tokens_out":783,"duration_ms":7249,"temperature":1.0,"reasoning_tokens":699,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T13:36:37.143380+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Test Proposition 4.1 directly: fix a shock $(u_l,u_r)$ with both states in the compact set $D$, take $u\\in S_{\\mathrm{weak}}$ to be a nontrivial wild solution, and check whether, for every Lipschitz shift $h$ with $h(0)=x_0$, the weighted flux functional $a_1(\\dot h\\,\\eta(u_-|u_l)-q(u_-;u_l))-a_2(\\dot h\\,\\eta(u_+|u_r)-q(u_+;u_r))$ is positive on a set of positive measure; one shock and one $u$ for which it is always positive disproves the global $a$-contraction estimate and collapses the proof.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the a-contraction-with-shifts estimate for large/extremal shocks used as Proposition 4.1's big-shock regime."},{"cited_title":"Golding, Sam G","cited_arxiv_id":null,"evidence_quote":"Supplies the sharp a-contraction estimate for small extremal shocks used as Proposition 4.1's small-shock regime."},{"cited_title":"Krupa, and Alexis F","cited_arxiv_id":null,"evidence_quote":"Establishes the modified front-tracking plus weight strategy and the weak-small-BV stability framework that this paper extends to large BV."},{"cited_title":"G lobal solutions to the isothermal Euler-Poisson system with arbitrarily large data","cited_arxiv_id":null,"evidence_quote":"Provides the TVD inequality for the field w = -ln tau recorded as Lemma 2.1, which controls shock classification changes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the equivalence of Eulerian and Lagrangian weak solutions used to transfer the main theorem to Eulerian coordinates."},{"cited_title":"Colombo and Nils H","cited_arxiv_id":null,"evidence_quote":"Part of the L1-theory synthesis in Theorem 2.1 that identifies the limit of the front-tracking shadows with the unique BV solution."},{"cited_title":"Unique solutions to hyperbolic conservation laws with a strictly convex entropy","cited_arxiv_id":null,"evidence_quote":"Completes Theorem 2.1, the uniqueness of BV entropy solutions, which fixes the target of the stability argument."},{"cited_title":"Krupa and Alexis F","cited_arxiv_id":null,"evidence_quote":"Provides the relative-entropy approximate-limit lemma (Lemma 7.1) that allows summing dissipation across interaction times."}],"review_version":1}