{"id":"dc287476-e784-402c-858e-1ab83c8fcf7d","arxiv_id":"2505.15078","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Small isothermal Euler shocks are stable and unique in the class of vanishing viscosity limits, even under large perturbations of finite relative entropy.","lead":"This paper proves that small-amplitude shock waves of the one-dimensional isothermal Euler equations are stable against large perturbations, provided the solutions are obtained as vanishing viscosity limits of Navier-Stokes flows. It is the first such stability result for the isothermal case, and it also gives global existence of strong solutions for degenerate viscosity isothermal Navier-Stokes equations.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 1.3(i) well-prepared initial data may be impossible for α≥1/2: the deferred construction must satisfy the one-sided derivative condition (1.9), which conflicts with approximating arbitrary finite-entropy u0 near large-v0 regions.","rationale":"The reader identified the missing well-prepared initial-data construction as the weakest assumption; I agree with that location. However, the paper's own Theorem 1.2 adds a specific technical condition, (1.9), and the isothermal BD functional is not a direct special case of the barotropic construction in [34]. The tension between (1.9) and approximating arbitrary finite-relative-entropy data is concrete and potentially fatal: a jump in u0 inside a region where v0 is very large cannot be approximated in L2 while keeping the one-sided derivative bound required by (1.9) and maintaining convergence of the relative entropy. This is not a disagreement with consensus; it is an internal gap in the argument, since the paper asserts Theorem 1.3(i) without proof. The contraction estimates in Sections 3–4 are detailed and appear coherent, but they only apply to data that the theorem must first produce; without the well-prepared sequence the uniform-in-ν bound (5.1) has no starting point. The proposed test isolates the obstruction: for a fixed large M, the derivative bound M^{-1/2} makes the jump approximation error nonvanishing, so (1.15) cannot hold. If the test reveals that a construction exists (for example by exploiting the freedom to modify vν0 on a set of vanishing measure while keeping relative entropy close), then the concern is resolved and the current conditional verdict can stand. As it stands, the paper should not be fully accepted unless the well-prepared data lemma for the isothermal functional, including (1.9), is supplied.","tokens_in":41042,"tokens_out":25268,"duration_ms":253080,"concrete_test":"Fix α=1/2, v̄=1, ū=0, and take U0=(v0,u0) with v0(x)=1+(M−1)χ_{(0,1)}(x), u0(x)=χ_{(0,1)}(x) for M=10^4 (so E0=Φ(M)+1/2<∞). Check whether there exists any sequence (vν0,uν0) satisfying (1.15) and (1.9). A decisive analytical check: prove that every absolutely continuous uν0 with ∂x uν0 ≤ M^{-1/2} has ‖uν0−u0‖_2 ≥ c M^{-1/4} for some c>0, while (1.15) forces vν0≈M on most of (0,1) and hence forces exactly that derivative bound. If the lower bound holds, no well-prepared sequence exists and Theorem 1.3(i) is false; if a construction with error→0 is exhibited, the concern is resolved.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim rests on Theorem 1.3(i), which asserts that every finite-relative-entropy datum U0 can be approximated by smooth well-prepared data satisfying (1.15). The proof is deferred to [34], whose barotropic construction does not cover the isothermal BD functional. There is an internal consistency problem with the existence theorem used for α∈[1/2,1]: Theorem 1.2 requires (1.9), i.e., ∂x uν0 ≤ (ρν0)^{1−α}. A finite-relative-entropy U0 is allowed to have an upward jump in u0 inside a region where v0 is very large, so ρ0 is very small. For any approximating sequence satisfying (1.15), relative-entropy convergence forces vν0→v0 in L1 on bounded sets, hence ρν0≈1/M in that region and (1.9) gives ∂x uν0 ≲ M^{-(1−α)}. Approximating a unit step of u0 in L2 by functions with derivative bounded above by K requires transition width at least c/K; the L2 error is bounded below by a positive constant depending on K. As M grows, K→0 and the minimal error does not vanish. Thus the asserted well-prepared construction is not merely omitted; it appears to be impossible for a concrete class of admissible initial data. If no such construction exists, the class of inviscid limits in Theorem 1.3 may be empty or exclude the stated \"any large perturbation\" case, and the stability/uniqueness conclusion collapses.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the one-dimensional isothermal Navier-Stokes system (1.1) with viscosity coefficient μ(v)=b v^{-α}, α∈[0,1]. It claims (i) global existence of large strong solutions for smooth initial data with positive density lower bound (Theorems 1.1 and 1.2), (ii) a contraction property for large perturbations of viscous shocks (Theorem 1.4), and (iii) stability and uniqueness of small-amplitude Riemann shocks of the isothermal Euler system in the class of vanishing viscosity limits (Theorem 1.3). The contraction proof in Sections 3–4 uses a relative entropy with a BD-type modified velocity, a carefully designed weight function, and a shift ODE. The final inviscid-limit argument in Section 5 is mostly deferred to the isentropic paper [34], with only the shift bound (1.18) proved in detail.","tokens_in":41392,"tokens_out":10534,"duration_ms":90003,"significance":"If the main theorem holds, it would provide the first stability-and-uniqueness result for Riemann shocks of the isothermal Euler system in a class of physical vanishing-viscosity limits, covering large perturbations with finite relative entropy. The a-contraction estimate of Theorem 1.4 is a substantial technical contribution, adapting the method of [33,34] to the isothermal BD functional. However, the paper contains several load-bearing gaps: the well-prepared initial data construction, the inviscid limit passage, and the global existence theorem for α∈[1/2,1] are not proved in the text and are only deferred to references. A concrete obstruction to the well-prepared construction for α≥1/2 is described below. The significance is therefore conditional on these gaps being closed.","major_comments":[{"comment":"The well-prepared initial data asserted in Theorem 1.3(i) are not constructed in the paper; the proof is deferred to [34], which treats the isentropic system. This is not a routine translation, because for α∈[1/2,1] the only existence theorem available (Theorem 1.2) requires the one-sided derivative bound (1.9). Consider an admissible initial datum U0 with v0=M on a set of positive measure and u0 containing an upward jump of order one in that region. Finite relative entropy allows such U0 for arbitrarily large M. The convergence in (1.15) forces vν0→M, hence ρν0≈1/M, so (1.9) gives ∂x uν0 ≤ M^{-(1-α)}. Approximating a unit jump with derivative bounded by this value in L2 requires a transition layer of length at least M^{1-α}, and the minimal L2 error is at least c M^{1-α}, which grows without bound as M increases. Thus no sequence satisfying (1.15) and (1.9) can exist for this datum. The paper neither proves nor cites an isothermal version of the well-prepared construction, and the cited barotropic construction in [34] does not address (1.9). This is load-bearing, since (1.15) is the bridge from arbitrary finite-relative-entropy data to the contraction estimate of Theorem 1.4.","section":"Theorem 1.3(i), Eq. (1.15), with the condition (1.9)"},{"comment":"The proof of Theorem 1.3 is not self-contained. After deriving the shift bound (1.18), the paper states that the existence of inviscid limits (1.16) and the stability estimate (1.17) follow 'essentially the same' as in [34], and gives no details. Because the present system is isothermal with a different BD-modified velocity (1.20), the weak convergence of (vν,uν) to a measure-valued limit and the lower-semicontinuity arguments for the extended relative entropy dΦ(v∞/v̄) require verification. The compactness results in [34] are adapted to the isentropic functional and may not carry over verbatim, especially in the presence of the singular part dvs. Since (1.16)–(1.17) form the main conclusion of the paper, this is a major gap.","section":"Section 5, proof of Theorem 1.3"},{"comment":"The global existence theorem for α∈[1/2,1] is not proved; the final paragraph of Section 2 states that the argument follows the barotropic case [32] via the active potential and then says 'Thus, we omit the proof.' Similarly, parts of Theorem 1.1 (the BD entropy estimate and the standard parabolic compactness) are deferred to [40]. While these techniques are standard, the isothermal case has a distinct density-bound argument (Proposition 2.1 uses logr instead of the barotropic ρ^{(γ-1)/2}), and the derivative condition (1.9) is a new structural restriction. Since the class X_T is the ambient space in which the contraction Theorem 1.4 is applied, an incomplete existence proof weakens the foundation of the paper.","section":"Section 2, Theorem 1.2"}],"minor_comments":[{"comment":"The notation ∫_R dΦ(v∞/v̄(x−X∞(·))) is ambiguous: the measure dΦ depends on the decomposition of v∞ into absolutely continuous and singular parts with respect to the shifted reference, and this dependence should be made explicit.","section":"Theorem 1.3, Eq. (1.17)"},{"comment":"The passage to the limit ε→0 in the test-function argument is terse; the text asserts continuity of ∫ ψ vν dx but does not verify all hypotheses of the dominated convergence theorem for the term involving hν and the viscous flux.","section":"Section 5, Eq. (5.3)"},{"comment":"The phrase 'unique in the class of inviscid limits' is not formalized; the paper should define the class of inviscid limits (e.g., subsequential limits for which (1.16) holds) and state the uniqueness assertion as a precise theorem.","section":"Theorem 1.3, statement"},{"comment":"Several key lemmas (Lemma 3.3, Lemma 3.4, Proposition 4.2) are quoted from [33] and [20] without proof; the paper should explicitly indicate which statements are new and which are directly taken from these references, and it should be noted that [20] is an arXiv preprint by overlapping authors.","section":"Sections 3–4, citations to [33] and [20]"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern about the well-prepared initial data appears genuine and is not addressed anywhere in the manuscript. For α∈[1/2,1], the condition (1.9) seems to obstruct approximating admissible data with upward velocity jumps in regions of large specific volume. This is not merely a missing proof of a routine construction; it may force a restriction of the admissible class in Theorem 1.3. The editor may wish to ask the authors to provide a complete construction or to state the theorem under a one-sided Lipschitz condition on u0. The paper also leans very heavily on [20,33,34], to the point that the main inviscid-limit proof is not present in the text."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know about this paper. First, it proves a genuine new contraction property for large perturbations of small-amplitude viscous shocks for the isothermal Navier-Stokes system, covering α∈[0,1]. The proof in Sections 3–4 is detailed and looks coherent; the expansion in shock size plus truncation is an honest extension of the Kang–Vasseur machinery. Second, the advertised stability theorem for isothermal Euler shocks, Theorem 1.3, leans on two deferred pieces that are more than routine: the well-prepared initial-data construction (i) and part of the inviscid-limit argument. Take Theorem 1.3 on faith at your peril.\n\nWhat is genuinely new: global existence for isothermal NS with degenerate viscosity (Theorem 1.1 fully proved; Theorem 1.2 sketched, with details deferred to [32]) and the first a-contraction estimate for the isothermal BD functional with shift. If you extract only the contraction theorem, the paper is a solid contribution.\n\nThe soft spot is the well-prepared data. For α∈[1/2,1], Theorem 1.2 requires ∂_x u0 ≤ ρ0^{1−α}. If the Euler datum has an upward jump in u0 inside a region where v0 is very large (ρ0 tiny), any smooth approximating sequence satisfying that one-sided bound has slope at most M^{−(1−α)}, which tends to 0. Approximating such a jump in L2 is then impossible; the minimal L2 error stays bounded away from zero. So the asserted construction (1.15) looks not merely omitted but impossible for a concrete class of finite-relative-entropy data. The authors defer to [34], but that is the barotropic case, and the isothermal BD functional is different. Unless there is an extra argument, Theorem 1.3(i) is false as stated for α≥1/2. This does not necessarily sink the α<1/2 case, where Theorem 1.1 needs no such derivative condition, but the theorem as written makes no such restriction.\n\nThe rest of Theorem 1.3 — the measure-valued limits, the shift bound — mostly follows [34] and [20] with the isothermal relative entropy; the shift estimate (1.18) is actually proved in the paper and checks out.\n\nVerdict: send it to a serious referee, with instructions to focus ruthlessly on Theorem 1.3(i). If the authors can prove the well-prepared construction in the isothermal setting, or restrict the theorem to α<1/2 and the existence theorem they actually prove, the contraction part deserves publication.","headline":"Strong contraction machinery for isothermal viscous shocks, but the advertised Euler stability theorem rests on a well-prepared-data construction that looks impossible for the degenerate-viscosity range.","tokens_in":41904,"tokens_out":2816,"would_cite":true,"duration_ms":26239,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q30","76N10","35B35"],"pacs":[],"model":"deepseek-v4-flash","headline":"Small-amplitude Riemann shocks of the one-dimensional isothermal Euler equations are stable and unique within the class of vanishing viscosity limits from large Navier–Stokes perturbations.","keywords":["Isothermal Euler","Isothermal Navier-Stokes","vanishing viscosity limit","relative entropy","shock stability","uniqueness","degenerate viscosity","large perturbations"],"falsifier":"The central claim would be falsified by finding a finite-relative-entropy initial datum for which no sequence of smooth data satisfying (1.15) exists, since Theorem 1.3(i) would then have an empty domain. It would also be falsified by producing $W$ in the setting of Proposition 4.2 with $\\int_0^1 W^2\\,dy\\le C_1$ and $\\sqrt{y(1-y)}\\,\\partial_y W\\in L^2(0,1)$ but $R_\\delta(W)>0$ for arbitrarily small $\\delta$, because that nonlinear Poincaré inequality is the step that closes the contraction estimate.","tokens_in":40865,"feed_emoji":"⚡","tokens_out":14518,"duration_ms":117904,"temperature":0.7,"pith_summary":"This paper establishes that small-amplitude shock waves of the one-dimensional isothermal Euler equations are stable and unique, provided solutions are understood as vanishing-viscosity limits of the associated isothermal Navier–Stokes system. The perturbations around the shock may be arbitrarily large in amplitude, as long as their initial relative entropy is finite. To reach this conclusion, the paper first proves global-in-time existence of strong solutions to the Navier–Stokes system with degenerate, density-dependent viscosity coefficients and possibly different far-field states, then proves a contraction property for large perturbations of viscous shocks that is uniform in the viscosity parameter. If the main theorem is correct, it supplies a physically natural class of Euler shock solutions in which stability and uniqueness hold without the usual small-BV or strong-trace assumptions.","feed_headline":"Isothermal gas shocks stay stable under arbitrary-size disturbances","feed_subtitle":"Global Navier-Stokes flows with large finite-entropy perturbations converge to the inviscid shock, up to a small shift.","key_machinery":"The load-bearing object is the BD-type relative entropy $$E((v_1,u_1)|(v_2,u_2))=\\Phi(v_1/v_2)+\\frac12\\left(u_1-\\frac{(v_1)_x}{$v_1^{{\\alpha+1}}$}-u_2+\\frac{(v_2)_x}{$v_2^{{\\alpha+1}}$}\\right)^2,$$ with $\\Phi(z)=z-1-\\log z$; it measures distance to a viscous shock using the effective velocity $h=u-v_x/v^{\\alpha+1}$. The contraction theorem shows that a weighted version of this functional, with monotone weight $a(\\xi)=1-\\frac{\\lambda}{\\varepsilon}(p(\\tilde v(\\xi))-p(v_-))$ and a shift $X(t)$, decays up to dissipative bulk terms. The shift is chosen by an ODE driven by the functional $Y$, and the decisive estimate is a nonlinear Poincaré inequality applied after rescaling the shock profile by $y=(p(v_-)-p(\\tilde v))/\\varepsilon$, which converts the dangerous quadratic terms into controllable ones.","core_discovery":"The central result, Theorem 1.3, states that for any small shock amplitude $\\varepsilon=|p(v_+)-p(v_-)|$ and any initial datum with finite relative entropy $E_0=\\int_{\\mathbb R}\\eta((v_0,u_0)|(\\bar v,\\bar u))dx$, the Navier–Stokes solutions converge, up to subsequence and a time-dependent shift $X_\\infty$, to an inviscid limit $(v_\\infty,u_\\infty)$ satisfying $$\\int_{\\mathbb R} d\\Phi(v_\\infty/\\bar v(\\cdot-X_\\infty(t)))+\\int_{\\mathbb R}\\frac{|u_\\infty-\\bar u(\\cdot-X_\\infty(t))|^2}{2}\\,dx \\le C E_0$$ and $|X_\\infty(t)-\\sigma t|\\le C(T)|v_--v_+|^{-1}(\\sqrt{E_0}+E_0)$. When $E_0=0$, the bound forces $v_\\infty=\\bar v$ and $u_\\infty=\\bar u$ almost everywhere, which gives uniqueness of the Riemann shock in this class. The smallness of the shock amplitude is needed for the contraction theorem, not for the inviscid-limit passage itself.","pith_inferences":["Editorial inference: the a-contraction machinery could plausibly be adapted to multidimensional planar shocks or to other pressure laws with logarithmic entropy, but the one-dimensional nonlinear Poincaré inequality would need to be re-derived for each new geometry.","Editorial inference: the $|v_--v_+|^{-1}$ factor suggests that very weak shocks may have slowly converging shift; the paper does not address whether this scaling is sharp.","Editorial inference: uniqueness in the paper is relative to the constructed vanishing-viscosity class, so other entropy solutions not obtained as such limits could in principle coexist; this is a limitation of the solution class, not a claim about all weak solutions."],"forward_implications":["A well-posedness class for small isothermal Riemann shocks emerges: perturbations may be arbitrarily large, with stability controlled solely by $E_0$, the initial relative entropy.","The shift bound $|X_\\infty(t)-\\sigma t|\\le C(T)|v_--v_+|^{-1}(\\sqrt{E_0}+E_0)$ gives quantitative control of the shock location; in particular $E_0=0$ forces $X_\\infty(t)=\\sigma t$ and recovers the shock exactly.","The global existence theorems supply uniform-in-$\\nu$ large strong solutions for $\\mu(\\rho)=\\rho^\\alpha$, $\\alpha\\in[0,1]$, with degenerate viscosity near vacuum and different far-field states, giving the solutions whose vanishing-viscosity limits Theorem 1.3 controls.","Stability and uniqueness hold without imposing BV or strong-trace conditions on the perturbation, because the relative-entropy framework is nonlinear and the shift absorbs the shock location uncertainty."],"supporting_citations":[{"why":"Supplies the inviscid-limit stability framework and the deferred construction of well-prepared initial data assumed in Theorem 1.3(i).","marker":"[34]"},{"why":"Develops the a-contraction method for large perturbations of barotropic Navier-Stokes shocks, adapted here to the isothermal case.","marker":"[33]"},{"why":"Provides the relative-entropy lemmas, weight-function estimates, shift ODE existence, and the proof strategy for the shift bound.","marker":"[20]"},{"why":"Gives the global large strong-solution proof structure for one-dimensional degenerate-viscosity Navier-Stokes used in Theorem 1.1.","marker":"[40]"},{"why":"Extends global strong solutions to degenerate viscosities alpha at least 1/2, supplying the argument behind Theorem 1.2.","marker":"[32]"},{"why":"Introduces the active potential used to remove the viscosity lower-bound dependence in the existence proof.","marker":"[14]"},{"why":"Foundational relative-entropy method on which the contraction and stability arguments are built.","marker":"[15]"},{"why":"Establishes uniqueness and L2-stability of Lipschitz solutions, the classical result this shock-stability theorem extends.","marker":"[17]"},{"why":"Introduces the BD entropy that underlies the effective-velocity functional E.","marker":"[3]"}],"fun_headline_variants":["Shock stability under arbitrary-size disturbances","Large perturbations cannot break isothermal shock stability","Inviscid limit yields Riemann shock stability for large noise","Even huge finite entropy preserves shock convergence","Small shock amplitude: any disturbance still converges"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The proof of Theorem 1.3(i) assumes that any initial datum with finite relative entropy can be approximated by smooth well-prepared data satisfying (1.15), including convergence of the modified relative entropy; the construction is not carried out here and is deferred to a reference.","fun_headline_variants_meta":{"raw":{"variants":["Shock stability under arbitrary-size disturbances","Large perturbations cannot break isothermal shock stability","Inviscid limit yields Riemann shock stability for large noise","Even huge finite entropy preserves shock convergence","Small shock amplitude: any disturbance still converges"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00087,"raw_usage":{"total_tokens":3765,"prompt_tokens":939,"completion_tokens":2826,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":555,"completion_tokens_details":{"reasoning_tokens":2759}},"tokens_in":555,"tokens_out":2826,"duration_ms":24596,"temperature":1.0,"reasoning_tokens":2759,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T15:24:08.060889+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"The central claim would be falsified by finding a finite-relative-entropy initial datum for which no sequence of smooth data satisfying (1.15) exists, since Theorem 1.3(i) would then have an empty domain. It would also be falsified by producing $W$ in the setting of Proposition 4.2 with $\\int_0^1 W^2\\,dy\\le C_1$ and $\\sqrt{y(1-y)}\\,\\partial_y W\\in L^2(0,1)$ but $R_\\delta(W)>0$ for arbitrarily small $\\delta$, because that nonlinear Poincaré inequality is the step that closes the contraction estimate.","supporting_citations":[{"cited_title":"Kang and A","cited_arxiv_id":null,"evidence_quote":"Supplies the inviscid-limit stability framework and the deferred construction of well-prepared initial data assumed in Theorem 1.3(i)."},{"cited_title":"Kang and A","cited_arxiv_id":null,"evidence_quote":"Develops the a-contraction method for large perturbations of barotropic Navier-Stokes shocks, adapted here to the isothermal case."},{"cited_title":"Mellet and A","cited_arxiv_id":null,"evidence_quote":"Gives the global large strong-solution proof structure for one-dimensional degenerate-viscosity Navier-Stokes used in Theorem 1.1."},{"cited_title":"Kang and A","cited_arxiv_id":null,"evidence_quote":"Extends global strong solutions to degenerate viscosities alpha at least 1/2, supplying the argument behind Theorem 1.2."},{"cited_title":"Constantin, T","cited_arxiv_id":null,"evidence_quote":"Introduces the active potential used to remove the viscosity lower-bound dependence in the existence proof."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Foundational relative-entropy method on which the contraction and stability arguments are built."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes uniqueness and L2-stability of Lipschitz solutions, the classical result this shock-stability theorem extends."},{"cited_title":"Bresch and B","cited_arxiv_id":null,"evidence_quote":"Introduces the BD entropy that underlies the effective-velocity functional E."}],"review_version":1}