{"id":"7df01ec2-acde-4732-8cbf-747235a39a67","arxiv_id":"2501.07831","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"low","formal_verification":"none","parameter_count":2,"one_line_summary":"The authors construct self-similar solutions of 1D compressible Euler in which a stationary vacuum boundary waits a finite time, then moves with the physical vacuum condition.","lead":"This paper constructs explicit solutions of the 1D compressible Euler equations where the vacuum boundary stays fixed for a finite time, then suddenly starts moving with physical vacuum behavior. It is the first rigorous construction of such waiting-time solutions in gas dynamics.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Lemma 4.7's B-to-D orbit is the crux, and its barrier trap is asserted without the required global comparison; if the ordering fails for some (γ,μ), the waiting-time construction collapses.","rationale":"The reader and I identify the same structural weak point: the global barrier comparison in Lemma 4.7 is stated rather than proved. I partially agree because the precise gap is sharper than just 'the ordering': the proof also needs the sign of F on the G=0 curve to justify the asserted divergence H'→+∞ at a first intersection. The concern is not an internal inconsistency; the inequalities are plausible and are likely supplied by a short algebraic lemma. But because Lemma 4.7 is the only mechanism producing the B-D orbit that terminates at the moving vacuum boundary, the theorem's correctness is conditional on this comparison. I therefore recommend conditional acceptance: add the missing algebraic verification before full acceptance. I do not see a reason to reject: the Poincaré-Dulac analysis at D is standard, Lemma 3.7's weak-solution selection is coherent, and the local expansions at B and D are consistent. The sign typo in (4.5) flagged by the reader is harmless. The precise gap is the unstated verification that the corridor H_G<H<H_sp is genuinely trapping on the whole interval (1,U_D).","tokens_in":20327,"tokens_out":16813,"duration_ms":174522,"concrete_test":"Verify by exact symbolic computation (Sturm sequences or cylindrical algebraic decomposition) the two inequalities for all γ∈(1,3), μ∈(0,1), U∈(1,2/(3-γ)): H_sp-H_G>0, equivalently (γ-1)^2 U(U+k2)/4 - (U-1)(U-μ)>0; and F>0 on G=0, equivalently U(U-1)(U-μ) - (U^2-k1U+μ)(U+k2)>0, with k1=((γ+1)+μ(3-γ))/2 and k2=2(1-μ)/(γ-1). If both hold, Lemma 4.7's trap is valid and the reader's concern is a presentation gap; if either fails, produce the explicit (γ,μ,U) counterexample and the B-D connection is not established.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The existence part of Theorem 2.3 rests on Lemma 4.7: a unique integral curve from B=(1,0) to D=(2/(3-γ), ((γ-1)/(3-γ))^2). The proof encloses this curve between H_G(U)=U(U-1)(U-μ)/(U+k2) and H_sp(U)=(γ-1)^2U^2/4 on 1<U<U_D. As written it checks only the slope comparison at B and then asserts the ordering H_G<H<H_sp. Two global facts are needed but not established: (i) H_sp>H_G throughout (1,U_D) so the upper barrier cannot be crossed before D; (ii) on the G=0 curve H_G, F>0 on (1,U_D), so a hypothetical first intersection with H_G would force H'=F/G→+∞ and contradict H-H_G→0 from above. Without (ii), a first intersection could have H'→-∞ and the crossing would not be ruled out. These are algebraic inequalities in three variables; they are not written out, and they are exactly what keeps the unique orbit from leaving the corridor. The rest of the phase-portrait construction is coherent, so this is the load-bearing soft spot, not a physical assumption.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs a two-parameter family of self-similar waiting-time solutions for the one-dimensional isentropic compressible Euler equations with vacuum, for adiabatic exponents γ∈(1,3) and similarity exponents μ∈(0,1). In the constructed solutions the vacuum boundary is stationary for t<0, the velocity and sound speed have at least C^1 regularity up to the boundary, and at t=0 the boundary begins to move according to b(t)=y_B t^{1/μ} while satisfying the physical vacuum condition for every t>0. A weak discontinuity emerges from the singular point along the sonic curve. The proof reduces the PDE to the self-similar ODE system (2.7), follows the explicit parabolic special solution H_sp from C to A to E to D, and then uses a phase-portrait analysis and a barrier argument to construct the unique orbit from B to D. Theorem 2.3 summarizes the existence and regularity statements.","tokens_in":20578,"tokens_out":23568,"duration_ms":220159,"significance":"If the proof gaps identified below are repaired, this is the first rigorous construction of waiting-time solutions in gas dynamics in which a stationary vacuum interface changes after finite time into a moving interface with the physical vacuum condition. The paper has substantial technical content: the explicit Burgers special solution is used in a clean way, Lemma 3.7 gives a genuine weak-solution uniqueness statement for the continuation through the point A, and the Poincaré-Dulac analysis at D is detailed and coherent. The construction also yields a genuinely two-parameter family of solutions, with μ and an integration constant as free parameters. The main weakness is localized in Lemma 4.7: the global barrier ordering that is the crux of the B-to-D connection is asserted rather than proved. This is a proof obligation, not a physical assumption, and it appears fixable by direct algebra.","major_comments":[{"comment":"The proof of the B–D connection is incomplete at the crucial step where a first intersection with the lower barrier H_G is ruled out. The text states that at a first intersection U_r one has lim_{U→U_r} H'(U)=+∞, but this limit is +∞ only if F(H_G(U),U)>0 on (1,U_D). Since G is positive just above G=0, the sign of H'=F/G is the sign of F; if F were negative at the intersection, one would have H'→−∞ and a transversal crossing of G=0 would not be excluded. The needed positivity of F on the curve G=0 is a finite algebraic inequality in (γ,μ) and is not supplied. The upper-barrier claim H<H_sp is likewise asserted without proof; it can be justified from H(1)<H_sp(1), uniqueness of solutions of (2.9), and the fact that the two curves meet only at D, but none of this appears. Finally, the conclusion H(U_D)=H_D presupposes that the solution exists up to U_D; this can be obtained from the two barriers together with the observation that H_G(U)>(U−1)^2 on (1,U_D), but that observation is also absent. These omissions are load-bearing because Lemma 4.7 is the only argument that the unique curve from B reaches D.","section":"Lemma 4.7"},{"comment":"As printed, the displayed formulas for the special solution have the exponent of |U_sp| with the opposite sign from the one obtained by integrating (3.21). Integration of dz/dU=−1/(μU)+(1/μ−1)/(U−U_C) gives |y|=K|U_sp|^{-1/μ}|U_sp−U_C|^{1/μ−1}, not K|U_sp|^{1/μ}|U_sp−U_C|^{1/μ−1}. With the printed exponent, the limits as y→0 and y→∞ in Lemmas 3.6 and 3.9 are interchanged, so those lemmas do not follow from the formulas as stated. The asymptotic formulas (3.26)–(3.27) and the calculation in Section 5 are consistent with the corrected exponent, so this appears to be a typographical error, but it should be fixed throughout.","section":"Equations (3.15), (3.23), (3.25), (3.35), (3.36)"}],"minor_comments":[{"comment":"The denominator in the formula for ds/dy is written as y((U−1)^2+H); to match the expansion of Δ=(U−1)^2−H and the subsequent singularity it should presumably be y((U−1)^2−H).","section":"Section 4, Eq. (4.5)"},{"comment":"The displayed derivative of H_G(U) has a typo: the second term in the numerator should be k_2(3U^2−2(1+μ)U+μ), not k_2(3U−2(1+μ)U+μ k_2). With the corrected formula the claimed positivity for U≥1 is straightforward to verify.","section":"Lemma 4.6"},{"comment":"The sentence 'Property 3 is a consequence of Property 2 and Lemma 4.5' is imprecise; the Hölder continuity up to O rests on the behavior at A and E as well, and the reader would benefit from an explicit reference to the expansions in Section 3 and Lemma 4.8.","section":"Section 5, Proof of Theorem 2.3"}],"recommendation":"major_revision","confidential_remarks":"The paper is significant and the overall strategy is credible, but the barrier lemma in Section 4 is genuinely load-bearing and is not proved as written. The missing algebraic inequalities are likely true and can be supplied in a revision, so I see no grounds for rejection. My recommendation differs from the reader's accept because the gap is in the central connection B-to-D rather than a purely local presentation issue."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is a genuinely new result. It gives the first rigorous self-similar waiting-time solutions for compressible Euler with physical vacuum after the boundary starts moving, and it does so for all γ∈(1,3) and μ∈(0,1). The construction is concrete: a Burgers/simple-wave leg C→A→E→D, then a new D→B trajectory that meets the vacuum boundary. The weak discontinuity across the sonic curve is a nice qualitative feature.\n\nWhat's good. The paper does real work. The special solution reduces to Burgers, Lemma 3.7 uses the weak formulation to argue the continuation through A is unique among self-similar ODE continuations, and the Poincaré–Dulac analysis at D correctly identifies the two branches and shows only the non-special branch can reach B. The expansions near B and D are carefully derived. I checked the main line: it's coherent. I also agree the proof is self-contained modulo standard ODE theory; the free parameters are explicit.\n\nThe soft spot is Lemma 4.7, the connection from B to D. The proof puts H(U) between H_G (the G=0 curve) and H_sp and says both are barriers. The upper bound H_G<H<H_sp is the load-bearing fact, but the paper only checks slopes at B and then asserts the ordering. For the lower barrier, the contradiction at a first intersection U_r relies on H'→+∞, which requires F>0 along H_G on (1,U_D). That positivity is not proven anywhere; if F were negative there, H'→−∞ and the curve could cross. Likewise H_sp>H_G on the whole interval is not shown. These are algebraic inequalities in γ and μ, likely true, but they are not written out. It's a proof obligation, not a physical assumption, and it is exactly what keeps the unique orbit from leaving the corridor. A referee should ask for this as a required revision.\n\nMinor: there's a sign slip in the display of ds/dy in Lemma 4.2 (plus H instead of minus), but the subsequent M computation uses the correct Δ derivative.\n\nBottom line: the central claim is probably right, and the structure is publishable. The paper deserves a serious referee; the gap in 4.7 is repairable with an explicit algebraic lemma. I'd send it to review but ask for the barrier verification before acceptance.","headline":"A solid, novel construction of waiting-time vacuum boundaries for 1D Euler, but Lemma 4.7's barrier argument has a real gap that needs an algebraic verification.","tokens_in":21114,"tokens_out":5634,"would_cite":true,"duration_ms":46716,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q31","35L65","76N10","35B40"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper constructs self-similar solutions in which a vacuum interface stays stationary for finite time, then starts moving with physical-vacuum behavior.","keywords":["waiting time solutions","vacuum free boundary","compressible Euler equations","self-similar solutions","physical vacuum condition","sonic curve","phase portrait analysis","adiabatic exponent"],"falsifier":"Numerically integrate the effective ODE (2.9) for a specific pair such as $\\gamma=1.1$, $\\mu=0.1$, starting from $B=(1,0)$ with slope $C_2=\\gamma(1-\\mu)/(1+k_2)$, and test whether the resulting curve crosses either the lower boundary $G=0$ or the upper special solution $H_{sp}$ before reaching $U_D=2/(3-\\gamma)$; a crossing before $U_D$ would falsify Lemma 4.7 and Theorem 2.3.","tokens_in":20112,"feed_emoji":"⏳","tokens_out":12108,"duration_ms":102225,"temperature":0.7,"pith_summary":"This paper proves that one-dimensional compressible Euler flows with a vacuum boundary can exhibit a waiting time: the interface remains fixed at x=0 for t<0, then at t=0 begins moving along a power-law curve and immediately satisfies the physical vacuum condition. The authors construct, for every adiabatic exponent $\\gamma\\in(1,3)$ and every similarity exponent $\\mu\\in(0,1)$, a continuum family of self-similar solutions realizing this transition. These are the first rigorous waiting-time solutions in gas dynamics, confirming a behavior long expected from formal asymptotic studies. The solutions are smooth away from the vacuum boundary and a sonic curve, Hölder continuous at the singular point, and Lipschitz across the sonic curve.","feed_headline":"Gas waits, then moves: first waiting-time Euler solutions","feed_subtitle":"For all γ∈(1,3), a stationary vacuum interface turns into a moving physical-vacuum boundary after finite time.","key_machinery":"The central object is the self-similar ODE system (2.7) in the similarity variables $(U,H)$, obtained from the Euler equations via the scaling (2.5)–(2.6). The proof follows a single trajectory in the $(U,H)$ phase plane connecting five critical points: $C$ (stationary boundary for $t<0$), $A$ (the origin in profile variables), $E$ (the point $x=0$ for $t>0$), $D$ (the first sonic point, a weak discontinuity), and $B$ (the moving vacuum boundary, the second sonic point). The trajectory starts on the special solution $H_{sp}(U)=\\frac{(\\gamma-1)^2}{4}U^2$ — the self-similar simple wave — which carries it from $C$ through $A$ to $D$. At $D$ the Poincaré–Dulac theorem classifies all integral curves through the node: the special solution is the unique curve with slope $C_1$, and every other curve enters with slope $C_2$, producing the Lipschitz (not $C^1$) crossing. The final leg from $B$ to $D$ is obtained by a barrier argument: the curve with the positive unstable slope at $B$ is trapped between the curve $G=0$ and the special solution, forcing it to reach $D$ and thereby close the trajectory.","core_discovery":"Theorem 2.3 asserts that for each $\\gamma\\in(1,3)$ and each $\\mu\\in(0,1)$ there exist self-similar $(\\rho,u)$ solving the vacuum free-boundary Euler equations with $b(t)=0$ for $t<0$ and $b(t)=y_B t^{1/\\mu}$ for $t\\ge 0$, with $y_B<0$. Before $t=0$ the solution is a simple wave, the special solution $H_{sp}(U)=\\frac{(\\gamma-1)^2}{4}U^2$, and the stationary boundary has at least $C^1$ velocity and sound speed up to the boundary. After $t=0$ the same special solution extends through the origin, then crosses the first sonic point $D$ with only Lipschitz regularity (a weak discontinuity), and terminates at a second sonic point $B$ that is the moving vacuum boundary, where the physical vacuum condition holds for every $t>0$. The passage through $D$ is forced by the Poincaré–Dulac classification: only one integral curve has the smooth slope at $D$, all others enter with a different slope, and the construction selects the latter to reach the vacuum boundary.","pith_inferences":["The same phase-portrait mechanism may extend to other similarity exponents $\\mu$ outside $(0,1)$ or to radial symmetry, where the analogous sonic-point classification would determine whether waiting-time behavior persists.","A direct numerical check of the barrier inequality in Lemma 4.7 over the full parameter range would either confirm the global ordering or uncover a counterexample near $\\gamma$ close to 1; the paper only performs local slope comparisons.","The construction suggests that the waiting time is not a special accident but a robust phenomenon selected by the requirement that the continuation be a weak solution across $t=0$; Lemma 3.7 shows any other extension through $A$ creates a jump that violates the weak formulation.","One could test whether the constructed solutions are stable under small perturbations of the initial profile, a question not addressed by the existence proof."],"forward_implications":["If the theorem is correct, the well-known breakdown of simple waves in Burgers' equation can be resolved inside the full Euler equations: the self-similar Burgers solution continues after its singularity as a Hölder continuous Euler solution with a moving physical-vacuum boundary.","The construction provides a rigorous example where the physical vacuum condition emerges dynamically from initial data that do not satisfy it, supporting the conjecture that such transitions occur generically for vacuum states.","Because finite speed of propagation and far-field cutoffs apply, the local self-similar solutions can be glued to produce finite-energy solutions exhibiting the same waiting-time transition.","The sonic curve emanating from the singular point is a new type of weak discontinuity in these flows, giving an explicit mechanism for loss of regularity."],"supporting_citations":[{"why":"Supplies the dichotomy of instantaneous versus stationary self-similar vacuum solutions that frames the waiting-time question.","marker":"[23]"},{"why":"Provides the formal asymptotic analysis of waiting-time solutions for the shallow water equations that this paper makes rigorous for the full Euler system.","marker":"[7]"},{"why":"Establishes well- and ill-posedness of various vacuum states and conjectures the transition to physical vacuum, which motivates the construction.","marker":"[19]"},{"why":"Contains the Poincaré–Dulac theorem used to classify integral curves through the sonic node D.","marker":"[1]"},{"why":"Gives vacuum dam-break solutions with non-smooth waiting-time behavior, providing a comparative example in the literature.","marker":"[3]"},{"why":"Supplies the well-posedness theory for the physical vacuum condition that the constructed moving boundary is required to satisfy.","marker":"[18]"}],"fun_headline_variants":["Gas waits, then moves: new Euler waiting-time solutions","Stationary to moving vacuum: Euler waiting-time solutions","Waiting gas, then motion: explicit Euler solutions for γ∈(1,3)","Hold on, then go: Euler equations admit waiting-time solutions","From pause to flow: continuum of Euler waiting-time solutions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The existence proof depends on the unproven-in-detail claim that the trajectory from the boundary point remains strictly between two comparison curves for the entire interval of $U$ values up to the sonic point; if that global ordering failed for some $\\gamma$ near 1, the connecting curve would not reach its target.","fun_headline_variants_meta":{"raw":{"variants":["Gas waits, then moves: new Euler waiting-time solutions","Stationary to moving vacuum: Euler waiting-time solutions","Waiting gas, then motion: explicit Euler solutions for γ∈(1,3)","Hold on, then go: Euler equations admit waiting-time solutions","From pause to flow: continuum of Euler waiting-time solutions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000744,"raw_usage":{"total_tokens":3311,"prompt_tokens":931,"completion_tokens":2380,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":547,"completion_tokens_details":{"reasoning_tokens":2294}},"tokens_in":547,"tokens_out":2380,"duration_ms":19566,"temperature":1.0,"reasoning_tokens":2294,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:36:35.462306+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically integrate the effective ODE (2.9) for a specific pair such as $\\gamma=1.1$, $\\mu=0.1$, starting from $B=(1,0)$ with slope $C_2=\\gamma(1-\\mu)/(1+k_2)$, and test whether the resulting curve crosses either the lower boundary $G=0$ or the upper special solution $H_{sp}$ before reaching $U_D=2/(3-\\gamma)$; a crossing before $U_D$ would falsify Lemma 4.7 and Theorem 2.3.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the dichotomy of instantaneous versus stationary self-similar vacuum solutions that frames the waiting-time question."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Contains the Poincaré–Dulac theorem used to classify integral curves through the sonic node D."},{"cited_title":"Pure Appl","cited_arxiv_id":null,"evidence_quote":"Supplies the well-posedness theory for the physical vacuum condition that the constructed moving boundary is required to satisfy."}],"review_version":1}