{"id":"a0d33cf8-2d45-4e4a-8ca8-f4b9fb1f25ce","arxiv_id":"2507.02224","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Small-amplitude viscous shock profiles exist and are unique for the one-dimensional Brenner-Navier-Stokes-Fourier system with arbitrary positive C2 transport coefficients, with explicit quantitative decay and shape estimates.","lead":"This paper proves that weak shock waves in a one-dimensional gas model with temperature-dependent viscosity and heat conduction exist as smooth, unique traveling wave profiles. The accompanying quantitative estimates on the profiles are the exact tools needed to prove stability of these shocks in the vanishing-viscosity limit.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Reduced vector field f is not shown to be smooth at ε=0 under C^2 coefficients; Fenichel's theorem may be inapplicable.","rationale":"The reader's weakest-assumption analysis focused on the ideal polytropic gas law, which is an explicit modeling restriction and not an internal flaw. My stress-test found a more immediate technical gap: the proof requires the reduced vector field f to be smooth across ε=0, but the stated C^2 regularity of τ, μ, κ only guarantees that g2 is C^2. Division by ε^2 of a C^2 function with a second-order zero typically yields only a continuous function, not a differentiable one, so the application of Fenichel's first theorem is not justified. The critical manifold computation and the estimates in Section 4 depend on this smoothness, either through the existence of the invariant manifold or through higher derivatives used to bound v''', u'', and θ''. The central construction is plausible and the theorem may be true, but the proof as written needs either a stronger regularity hypothesis (e.g., C^3 coefficients) or a separate argument establishing sufficient regularity of the reduced vector field. Hence I recommend CONDITIONAL rather than unconditional ACCEPT.","tokens_in":15,"tokens_out":47358,"duration_ms":1065913,"concrete_test":"Let τ(θ)=1+(θ−θ−)^{5/2} (C^2 but not C^3), μ=κ=1, γ=7/5, R=1, v−=θ−=1. Differentiate (2.6) three times with respect to ε at (v−,0,0,0) and solve for ∂^3g2/∂ε^3. If this third derivative is unbounded or undefined because it involves τ'''(θ−), then the Taylor expansion (3.8) cannot be extended to a C^1 quotient, so 1/ε^2 g2 and hence f in (3.4) fail to be C^1 at ε=0, invalidating the invocation of Proposition 2.1 for C^2 coefficients.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The proof's reduction to geometric singular perturbation theory hinges on the claim that the function f in (3.4) is smooth on a full neighborhood of ε=0, justified only by the Taylor expansion (3.8) showing that g2(v−+εw0, ε^2w1, ε^2w2, ε) has zero constant and first-order terms. But the transport coefficients are only C^2, so the implicit function theorem gives g2 merely C^2. A C^2 function with a second-order zero does not become smooth after division by ε^2: e.g., F(ε)=ε^{5/2} is C^2, while F/ε^2=ε^{1/2} is not differentiable at 0. The coefficients may be chosen to create such terms (e.g., τ(θ)=1+(θ−θ−)^{5/2}), so there is no reason for 1/ε^2 g2 to be C^1. Consequently f may be only continuous at ε=0, and Fenichel's first theorem (Proposition 2.1), which requires smoothness, does not apply as stated. Moreover, Section 4 differentiates the invariant manifold up to order three to prove (1.8)-(1.9); this requires more regularity of the manifold than the C^2 data provides. The theorem might still be true, but the proof as written requires C^3 (or a finer argument controlling the remainder) for the stated class.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies traveling wave solutions (viscous shocks) for the one-dimensional Brenner-Navier-Stokes-Fourier system in Lagrangian coordinates, with temperature-dependent transport coefficients τ(θ), μ(θ), κ(θ). The main results, Theorem 1.1 and Theorem 1.2, assert that for any positive C^2 coefficients and any left state, sufficiently small-amplitude right states satisfying the Rankine-Hugoniot and Lax entropy conditions admit a unique monotone traveling wave profile, with quantitative exponential decay and derivative estimates, including the sharp profile law ṽ' = A(v+−ṽ)(ṽ−v−)+O(ε)(v+−ṽ)(ṽ−v−). The proof uses integration of the ODE system, elimination of the velocity variable, the implicit function theorem to express temperature and its derivative in terms of v and its derivatives, and then geometric singular perturbation theory via Fenichel's first theorem. This generalizes the authors' earlier constant-coefficient work [13] and is motivated by the stability analysis in [12].","tokens_in":19313,"tokens_out":9420,"duration_ms":106392,"significance":"If the proof is correct, the result is a useful and nontrivial generalization of the constant-coefficient existence theory to physically realistic temperature-dependent coefficients, and the explicit constant A in (1.15) and the quantitative estimates in Theorem 1.2 are precisely the type of information needed for the contraction-stability approach of [12]. The geometric reduction using the implicit function theorem is a sound strategy and could apply to other systems with nonlinear coefficients. However, the manuscript's central regularity claim—that the reduced vector field f in (3.4) extends smoothly to ε=0 under only C^2 assumptions—is not justified and is load-bearing for the application of Fenichel's theorem and for the derivative estimates in Section 4.","major_comments":[{"comment":"The assertion that f can be smoothly extended to ε=0 is not supported by the arguments given. The implicit function theorem applied to the C^2 functions h1 and h2 yields g1 and g2 only of class C^2. The Taylor expansion (3.8) shows only that the limit of ε^{-2} g2(v_- + εw0, ε^2 w1, ε^2 w2, ε) exists as ε→0; it does not establish that the quotient is differentiable, let alone smooth. For example, φ(ε)=ε^{5/2} is C^2 with φ(0)=φ'(0)=0, while φ(ε)/ε^2=ε^{1/2} is not differentiable at 0. Since Proposition 2.1 requires the vector field to be smooth, the application of Fenichel's first theorem to (3.1) is not justified as written. This is a load-bearing gap: the invariant manifold parametrization (3.11) and all subsequent estimates depend on the smoothness of f near ε=0.","section":"Section 3, Eq. (3.4) and (3.8)"},{"comment":"The estimates |w''_1(z)| ≤ C|w'_0(z)| and |w'''_1(z)| ≤ C|w'_0(z)| are obtained by differentiating (4.3) up to three times, which requires the function s1(w0,ε) to be of class C^3 in both variables. Even if the reduced vector field f were smooth, this is a strong regularity requirement; under the stated C^2 hypotheses on τ, μ, κ, the manuscript does not establish the needed regularity of the invariant manifold. Please either prove the required regularity of s1 or rework the argument so that (1.8)-(1.9) follow from the lower-regularity information that is actually available. As written, the proof of the quantitative estimates requires more regularity than the hypotheses provide.","section":"Section 4, after Eq. (4.3)"},{"comment":"The theorem is stated for arbitrary positive C^2 coefficients, but the proof appears to require considerably more regularity: at minimum, the reduced vector field must be C^1 for a finite-regularity Fenichel-type result, and the derivative estimates in Section 4 seem to require at least C^3 regularity of the invariant manifold. It would be helpful if the authors specified precisely which regularity of τ, μ, κ is needed for each step; if the intended scope is indeed C^2, a substantially different argument or an approximation argument is required.","section":"General scope"}],"minor_comments":[{"comment":"The phrase 'admits tor viscous shocks' appears to be a typo; it should likely read 'admits viscous shocks' or 'admits two viscous shocks'.","section":"Section 1.1, paragraph 1"},{"comment":"The arguments of h1 and h2 are written as (ṽ, θ̃, ṽ′, ṽ′′, θ̃, ε), but from the context and from (2.6) the fourth argument should be θ̃′; please correct the notation.","section":"Eq. (2.4)"},{"comment":"The term 'l.o.t.' in the Taylor expansion is not precise. Given the C^2 regularity of g2, the remainder is o(ε^2), not necessarily O(ε^3); specifying the exact remainder order would clarify why the limit, but not the smoothness, of the quotient follows.","section":"Eq. (3.8)"},{"comment":"'It also requires to obtain the bounds' is grammatically awkward; consider rephrasing to 'It is also necessary to obtain the bounds for |ũ''_ε| and |θ̃''_ε|'.","section":"Section 4, text after (4.1)"},{"comment":"'Sincey(ξ)' is missing a space; it should read 'Since y(ξ)'.","section":"Proof of Corollary 1.1"}],"recommendation":"major_revision","confidential_remarks":"The central geometric idea is sound and the paper fills a real gap in the literature, but the regularity issue in Section 3 is the key obstacle. If the authors can supply a rigorous finite-regularity version of the reduction, or alternatively restate the theorem for C^∞ (or sufficiently high finite regularity) coefficients, the paper would be a solid contribution. The current manuscript, however, cannot be accepted with the proof as written."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, the paper genuinely does something new: it proves existence, uniqueness, and quantitative profile estimates for small-amplitude viscous shocks in the one-dimensional BNSF system with arbitrary positive C^2 temperature-dependent coefficients, generalizing the constant-coefficient result [13] and supplying the exact estimates the companion stability paper [12] needs. The use of the implicit function theorem to eliminate θ and reduce to a singular perturbation problem is substantive and sensible, and the critical manifold calculation giving the constant A is careful. The paper reads as honest, workmanlike applied PDE, not a cosmetic generalization.\n\nThe soft spot is a regularity gap in the geometric singular perturbation reduction. The proof claims that f in (3.4) extends smoothly to ε=0, based on the Taylor expansion (3.8) of g2. But the coefficients are only C^2, and because h1 contains τ'(θ), the implicit function theorem gives g2 only C^1 (or at best C^2 in some directions). A C^2 function with a second-order zero does not become smooth after division by ε^2—ε^2|ε| is a counterexample. So the conclusion that (1/ε^2)g2 is smooth on V is not justified, and the application of Fenichel's theorem as stated (Proposition 2.1, which requires smooth f) is not legitimate. Section 4 then differentiates the invariant manifold up to third order, which demands even more regularity. This is not a manufactured objection; the gap is in the written proof. It may be fixable by assuming C^∞ coefficients, by a finite-regularity GSP argument, or by approximation, but none of that is in the paper.\n\nThe rest of the proof—uniqueness via the unstable manifold, exponential decay, and the ratio estimates—is standard and appears sound once the invariant manifold exists. The ideal polytropic gas law is a stated restriction, not a hidden flaw.\n\nWho is this for? Mathematical fluid dynamicists working on viscous shock profiles and their stability, especially those using contraction-with-shifts methods. I would send it to a serious referee, but the referee should insist on a fix for the regularity issue before acceptance. I would not cite it as a finished result until that is sorted.","headline":"Genuinely new result and a clever reduction, but the proof has a real regularity gap at the heart of the GSP argument.","tokens_in":19861,"tokens_out":5530,"would_cite":false,"duration_ms":61044,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76N15","35Q30","35C07","35K65"],"pacs":[],"model":"deepseek-v4-flash","headline":"For any positive temperature-dependent viscosity and heat-conduction laws, the Brenner–Navier–Stokes–Fourier system has a unique small-amplitude viscous shock connecting admissible end states.","keywords":["traveling wave solutions","viscous shocks","Brenner-Navier-Stokes-Fourier system","temperature-dependent transport coefficients","geometric singular perturbation theory","Rankine-Hugoniot condition","Lax entropy condition","polytropic gas"],"falsifier":"Fix an admissible case, for example a 3-shock with the polytropic law, $\\gamma=1.4$, $\\tau(\\theta)=1+\\theta^2$, $\\mu(\\theta)=\\theta^2$, $\\kappa(\\theta)=\\theta^2$, and amplitude $\\varepsilon=10^{-3}$. Numerically integrate the traveling-wave ODE system (3.15) from near the left state toward the right state and compute the ratio $\\tilde v'/((v_+-\\tilde v)(\\tilde v-v_-))$ across the layer. If the ratio departs from the constant $A$ in (1.15) by more than $O(\\varepsilon)$, or if two distinct monotone connecting orbits emerge for the same end states, the sharp-profile law or the uniqueness claim is false.","tokens_in":18836,"feed_emoji":"🌊","tokens_out":19515,"duration_ms":181289,"temperature":0.7,"pith_summary":"The paper proves that the one-dimensional Brenner–Navier–Stokes–Fourier (BNSF) system — a two-velocity reformulation of gas dynamics — admits a unique monotone traveling wave (viscous shock) connecting any two nearby states that obey the Rankine–Hugoniot and Lax entropy conditions, for arbitrary positive temperature-dependent viscosity, heat conductivity, and the Brenner coefficient. This goes beyond the previous constant-coefficient result by allowing the transport coefficients to be any positive $C^2$ functions of temperature, as kinetic theory predicts for real gases. The proof reduces the traveling-wave ODEs to a single higher-order equation using the implicit function theorem, then applies geometric singular perturbation theory to a one-dimensional critical manifold. Along with existence and uniqueness, the paper establishes quantitative estimates — exponential decay at a rate proportional to the shock amplitude, comparable derivative magnitudes, and a sharp logistic-type profile law — that a companion stability study requires. If the theorem is correct, these profiles are the building blocks for studying inviscid limits from this class of viscous systems.","feed_headline":"Viscous shocks exist and are unique for temperature-dependent BNSF gas","feed_subtitle":"Any positive viscosity and heat-conduction laws now admit exact shock profiles, at small amplitude.","key_machinery":"The central mechanism is a geometric singular perturbation reduction. The traveling-wave ODE system is integrated to a first-order system; using the polytropic gas law, the implicit function theorem eliminates the temperature variable and its derivative in favor of $v$, $v'$, $v''$, producing a single third-order ODE. A scaling $v=v_-+\\varepsilon w_0(\\varepsilon\\xi)$, $v'=\\varepsilon^2 w_1(\\varepsilon\\xi)$, $v''=\\varepsilon^2 w_2(\\varepsilon\\xi)$ turns this into the singularly perturbed system (3.1), whose critical manifold — the set where the fast variables are equilibrated — is the one-dimensional curve $M_0=\\{w_2=0,\\ w_1=A(w_0-w_0^2)\\}$. A standard geometric singular perturbation result, Fenichel's first theorem, asserts that this normally hyperbolic manifold persists as a nearby locally invariant manifold for small $\\varepsilon$; the two critical points $(0,0,0)$ and $(1,0,0)$ lie on it, yielding the connecting orbit. The quantitative estimates come from the same manifold: the profile law $\\tilde v'/((v_+-\\tilde v)(\\tilde v-v_-))=A+O(\\varepsilon)$ is literally the parametrization $w_1=A w_0(1-w_0)+\\varepsilon s_1(w_0,\\varepsilon)$.","core_discovery":"The paper's central claim is that for any positive $C^2$ functions $\\tau(\\theta)$, $\\mu(\\theta)$, $\\kappa(\\theta)$ and any left state $(v_-,u_-,\\theta_-)$, every right state satisfying the Rankine–Hugoniot and Lax entropy conditions with sufficiently small amplitude $\\varepsilon=|v_--v_+|$ is connected to the left state by a unique monotone traveling wave $(\\tilde v,\\tilde u,\\tilde\\theta)(x-\\sigma t)$. Theorem 1.2 quantifies the profile: the wave approaches its end states exponentially at rate $C\\varepsilon$; the three derivatives are comparable, $|\\tilde v'|\\sim|\\tilde u'|\\sim|\\tilde\\theta'|$; and the shock layer obeys the sharp law $\\tilde v'=A(v_+-\\tilde v)(\\tilde v-v_-)+O(\\varepsilon)(v_+-\\tilde v)(\\tilde v-v_-)$, with an explicit constant $A$ built from the transport coefficients at the left state and the sound speed. The wave is monotone in each component, with signs depending on whether it is a 1-shock or a 3-shock. The quantity $A$ is not an afterthought: it is exactly the slope of the critical manifold in the reduced system, so the theorem identifies the precise shape of the viscous shock at small amplitude.","pith_inferences":["Beyond the paper: the logistic profile law implies the small-amplitude shock thickness is of order $1/(A\\varepsilon)$, so measuring the layer width at two different amplitudes would give a direct numerical estimate of the constant $A$.","Beyond the paper: because the construction needs only an invertible implicit-function Jacobian and a normally hyperbolic critical manifold, the same reduction is likely to work for other temperature-dependent viscous systems, such as Navier–Stokes–Korteweg, provided an analogous reduced ODE can be formed.","Beyond the paper: the one-dimensional unstable-manifold argument used for uniqueness does not by itself classify non-monotone heteroclinic orbits; ruling those out would require controlling trajectories that leave the invariant manifold and return, which the present theorem does not address."],"forward_implications":["Every admissible right state with small amplitude is connected to the left state by exactly one monotone viscous shock, up to translation, so the small-amplitude shock family is fully classified near the end state.","The profile law makes the shock layer approximately logistic: with $y=(\\tilde v-v_-)/\\varepsilon$, one has $dy/d\\xi=A\\varepsilon y(1-y)+O(\\varepsilon^2)$, so the wave shape is a rescaled logistic curve whose thickness is set by $A\\varepsilon$.","The exponential decay and derivative-ratio estimates supply exactly the uniform bounds required by the companion contraction-stability argument, making these waves usable as reference solutions for stability of Riemann shocks.","For a 1-shock the constructed profile satisfies $\\tilde v'<0$, $\\tilde u'<0$, $\\tilde\\theta'>0$, and for a 3-shock the signs are opposite; the fixed monotonicity pattern is part of the uniqueness statement."],"supporting_citations":[{"why":"Supplies the constant-coefficient predecessor whose geometric singular perturbation framework this paper extends to temperature-dependent transport coefficients.","marker":"[13]"},{"why":"Companion stability paper whose required quantitative estimates motivate Theorem 1.2 and give the Lagrangian-coordinate formulation of the BNSF system.","marker":"[12]"},{"why":"Provides the geometric singular perturbation persistence theorem (Fenichel's first theorem) used to construct the invariant manifold near the critical curve.","marker":"[20]"},{"why":"Supplies the unstable manifold theorem for ODEs on which the uniqueness proof rests.","marker":"[33]"},{"why":"Provides the Gronwall-type lemma that converts the profile-law estimate into exponential decay toward the end states.","marker":"[21]"},{"why":"Introduces the volume-velocity constitutive relation that defines the BNSF system studied here.","marker":"[4]"}],"fun_headline_variants":["Unique viscous shocks for BNSF with any positive coefficients","Small-amplitude shocks proven for temperature-dependent BNSF","Exact shock profiles for BNSF with temperature-dependent coefficients","BNSF viscous shocks: existence and uniqueness at small amplitude","Temperature-dependent BNSF admits unique viscous shock waves"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The proof assumes the gas obeys the ideal polytropic law $p=R\\theta/v$ with internal energy $e=R\\theta/(\\gamma-1)$; the shock-speed formula, the implicit-function reduction, and the explicit constant $A$ are all derived from this law, so a different equation of state would require redoing the whole geometric construction.","fun_headline_variants_meta":{"raw":{"variants":["Unique viscous shocks for BNSF with any positive coefficients","Small-amplitude shocks proven for temperature-dependent BNSF","Exact shock profiles for BNSF with temperature-dependent coefficients","BNSF viscous shocks: existence and uniqueness at small amplitude","Temperature-dependent BNSF admits unique viscous shock waves"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00048,"raw_usage":{"total_tokens":2409,"prompt_tokens":1015,"completion_tokens":1394,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":631,"completion_tokens_details":{"reasoning_tokens":1312}},"tokens_in":631,"tokens_out":1394,"duration_ms":98794,"temperature":1.0,"reasoning_tokens":1312,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:36:44.259886+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fix an admissible case, for example a 3-shock with the polytropic law, $\\gamma=1.4$, $\\tau(\\theta)=1+\\theta^2$, $\\mu(\\theta)=\\theta^2$, $\\kappa(\\theta)=\\theta^2$, and amplitude $\\varepsilon=10^{-3}$. Numerically integrate the traveling-wave ODE system (3.15) from near the left state toward the right state and compute the ratio $\\tilde v'/((v_+-\\tilde v)(\\tilde v-v_-))$ across the layer. If the ratio departs from the constant $A$ in (1.15) by more than $O(\\varepsilon)$, or if two distinct monotone connecting orbits emerge for the same end states, the sharp-profile law or the uniqueness claim is false.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the constant-coefficient predecessor whose geometric singular perturbation framework this paper extends to temperature-dependent transport coefficients."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the geometric singular perturbation persistence theorem (Fenichel's first theorem) used to construct the invariant manifold near the critical curve."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the unstable manifold theorem for ODEs on which the uniqueness proof rests."},{"cited_title":"Kang and A","cited_arxiv_id":null,"evidence_quote":"Provides the Gronwall-type lemma that converts the profile-law estimate into exponential decay toward the end states."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the volume-velocity constitutive relation that defines the BNSF system studied here."}],"review_version":1}