{"id":"3afce45d-8cc7-4ed9-9156-945f6f59f7a3","arxiv_id":"2607.24994","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":3,"one_line_summary":"For symmetric data, arbitrarily small diffusion selects a unique linearly stable time-periodic solution from the infinite family of relaxation oscillations present at zero viscosity.","lead":"Small viscosity collapses an infinite family of relaxation oscillations in a conservation law coupled to a fast ODE down to one linearly stable periodic orbit. The result gives a rigorous selection mechanism for time-periodic dust-void-like profiles in a simplified plasma-inspired PDE–ODE model.","discovery_kind":"extension","skeptic_critique":{"model":"moonshotai/kimi-k3","headline":"The a-posteriori shock-free slope bound (eqs. 4.19/4.23) is the real soft spot, but it limits the basin rather than undermining the stated theorem — same concern as the reader, and it does not move the verdict.","rationale":"The reader identified the same load-bearing point — the a-posteriori closure of the slope lower bound — and correctly classified it as a scoping limitation rather than an internal inconsistency. My independent pass confirms this and finds no stronger objection. The proof's logical structure is sound: the bound (4.19)/(4.23) is assumed for the GSPT and averaging steps, the periodic orbit and its spectral gap are constructed under it (Props. 5.6–5.10, Cor. 5.11), and §5.3 then proves the bound persists for data in the stated H² neighbourhood, which is exactly what Theorem 2.7 claims. The condition η < c/log(ε^{-1}) is consistently derived from requiring e^{-κ/η} = O(ε²) in Prop. 4.10, and the η^{1/3}ε remainder is genuinely lower order than the O(ε) spectral gap, so the stability conclusion is not threatened by the truncation. The numerics are finite-domain and code-free, which keeps reproducibility modest, but they are corroborating rather than load-bearing. The authors themselves flag the open territory (large-slope, non-monotone, non-symmetric data; possible longer-period orbits and MMOs), so nothing is hidden. Because the concern limits the basin but does not undermine the stated claim, and the reader's ACCEPT already reflects this, I recommend the verdict stand unchanged. The proposed numerical basin-mapping test would nonetheless be worth running, as it would determine whether the restriction is essential or an artifact of the proof technique.","tokens_in":45026,"tokens_out":2059,"duration_ms":65542,"concrete_test":"Numerically map the actual basin against the proved one. Using the paper's finite-domain scheme (§3, parameters as in Fig. 2), take symmetric monotone tanh-type initial profiles n_0(x) = ½(1+tanh(β(x−L/2))) and sweep the slope β from small values up past the Proposition 2.5 / Remark 2.6 threshold C ≈ 1/(8∆M). For each run, track min over one period of the discrete analogue of ∂_m X(θ,0) (equivalently max_x ∂_x n at the symmetry point, since ∂_m X(θ,0) = 1/(2∂_x n(t,0))). If this minimum stays bounded below by an ε-independent constant and the trajectory converges to the same selected orbit, the H²-closeness restriction is merely technical and the theorem's basin is larger than proved. If instead min ∂_m X(θ,0) → O(ε) or below and the (E,M) trajectory departs from the relaxation cycle (folded-node signatures, MMOs as in Fig. 9), the restriction is essential and the theorem's scoping is sh","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (Theorem 2.7) is that small viscosity selects a unique, linearly stable relaxation-periodic orbit. Every step of the proof architecture — the reciprocal-function formulation (Prop. 4.1), the GSPT control of the radial coordinate (Prop. 4.5, via [40]), the θ-time change (Prop. 4.8), and the averaging (Props. 4.4/4.11) — is conditioned on ∂_m X(θ,0) ≥ q_0 > 0 uniformly along the cycle. This bound is never proved a priori from the dynamics; it is only closed a posteriori in §5.3 (Corollary 5.14) for initial data already close to the selected orbit in the strong norm (5.27), which is then translated back to H²-closeness. This is exactly the reader's weakest_assumption, and on re-reading I confirm it is the genuinely load-bearing point: the normal-switching condition (4.21) fails precisely when ε/∂_m X(θ,0) ceases to be small, i.e., near shock formation, and Remark 4.7 explicitly notes that folded-node/MMO behaviour becomes possible there. However, I stress-test whether this invalidates the theorem and conclude it does not: Theorem 2.7 states the H²-neighbourhood hypothesis (2.9) as an assumption, so the conditional closure is logically sound — the bound is assumed, the orbit is constructed, and then the bound is verified to persist inside the basin. There is no circularity: §5.3 proves persistence forward in time for data in the basin, which is all the claim requires. The honest cost is scoping, not correctness: nothing is said about data with slope near or above the Proposition 2.5 threshold C ≈ 1/(8∆M), and the \"selection from an infinite-dimensional family\" narrative in the abstract is strictly stronger than what is proved (selection holds only among orbits reachable from the H² ball; longer-period orbits and non-monotone data are openly excluded). One secondary, lesser point: Prop. 4.5 treats p(y) = 1/∂_m X(θ,0) as a given bounded perturbation when invoking [40], which is legitimate only because (4.19) is assumed — another manifestation of the same","agreement_with_reader":"agree"},"referee_report":{"model":"moonshotai/kimi-k3","summary":"The paper studies a viscous scalar conservation law on R, ∂_t n = −a(E,M)∂_x[n(1−n)] + ε∂_{xx}n, globally coupled through M(t) = −∫_{−∞}^0 n dx to a fast ODE ηĖ = g(E,M) of van der Pol type. For ε = 0 the (E,M) dynamics closes and exhibits a relaxation limit cycle, and each monotone initial profile n_0 produces its own periodic modulation: an uncountable family of cycles with a common period (Prop. 2.5). The main result (Thm. 2.7) is that for 0 < ε ≪ 1 and 0 ≤ η < c/log(ε⁻¹), diffusion selects an isolated periodic orbit, linearly asymptotically stable with rate O(ε) for symmetric, monotone, H²-nearby data, with period converging to the inviscid one. The proof uses a reciprocal-function formulation (Prop. 4.1), polar-like coordinates and GSPT control of the radial variable (Prop. 4.5, via [40]), a θ-time change (Prop. 4.8), averaging to an autonomous equation for Ȳ (Props. 4.4/4.11), a unique stationary state selected by the conserved quantity K (Prop. 5.3), a spectral-gap analysis of a rank-one Sturm–Liouville perturbation (Props. 5.6–5.10), and an implicit-function/Poincaré argument for the untruncated system (§6). Global existence is proved in Appendix A; numerics in §3 support the O(ε) decay rate.","tokens_in":45407,"tokens_out":8519,"duration_ms":137141,"significance":"If correct, this is a strong and, to my knowledge, first rigorous result showing that vanishing viscosity collapses an infinite-dimensional continuum of relaxation-oscillation cycles to a single stable orbit in a PDE–ODE coupled conservation law, a selection mechanism relevant to void oscillations in dusty plasmas. The argument is a forward derivation: the averaged nonlinearity G (eqs. (4.9)/(4.38)) and the conserved quantity K are computed from the model, not fitted. The spectral analysis of §5.2 — reducing stability of a non-self-adjoint rank-one perturbation to the Laplace transform of a positive measure — is elegant and self-contained. The paper ships a quantitative, falsifiable prediction (decay rate O(ε) in (2.10)) that is tested numerically (Fig. 7: fitted slopes s±(ε) ≈ −2ε, −2.3ε). Limitations (basin restriction, possible MMOs near shock formation) are disclosed honestly in §2.3 and Remark 4.7.","major_comments":[{"comment":"The application of the Szmolyan–Wechselberger results [40] is not directly justified: [40] concerns autonomous slow–fast ODEs, whereas here the slow vector field g_1 = a/4 − εp(y) contains p(y) = 1/∂_m X(θ,0), which is not a function of (x,y,z) alone but is slaved to the PDE component X. The lower bound (4.19) controls |p|, but the GSPT estimates (slow-manifold attraction, fold extension, Poincaré map) are C^1-robust statements and require control of ∂_θ p as well, i.e. of ∂_θ∂_m X(θ,0). Since the radial asymptotics (4.18) underpins both the time change (Prop. 4.8) and the averaging (Prop. 4.11), the uniformity of [40]'s estimates in p (and its θ-derivative, controllable via (4.26) and the H²-closeness) should be stated and verified explicitly.","section":"§4.2.2, Proposition 4.5"},{"comment":"This remark performs load-bearing work: Corollary 5.14 establishes persistence of the bound ∂_m X(θ,0) ≥ q_0 > 0 (equivalently (5.23)) only for the truncated averaged equation (5.1), yet assumptions (4.19)/(4.23) are imposed on the full O(ε²)- (resp. O(ε²+η^{1/3}ε)-) perturbed dynamics. The transfer of the H²-ball attraction and of the pointwise lower bound to the untruncated system is exactly what closes the a-posteriori shock-free hypothesis inside the basin of Thm. 2.7. A one-sentence 'perturbation argument' is not sufficient here; please supply a proof or a detailed sketch (e.g. a Duhamel estimate for the perturbed semigroup combined with Lemma 5.12).","section":"§5.3–§6, Remark 6.1"}],"minor_comments":[{"comment":"Notation collision: n∗ denotes the fixed tanh-type reference profile in §2.1 (Thm. 2.1) and the periodic orbit n∗(t,x) in Thm. 2.7. Please distinguish them.","section":"§2.1 vs §2.3"},{"comment":"The translation of K(0) > z_0/2 appears to have the inequality reversed: K = −2M− + O(ε) > z_0/2 = 2(M+ − M−) + O(η^{1/3}) yields 0 > 2M+, i.e. M+ < O(ε) + O(η^{1/3}). The displayed '−2M− + O(ε) ≤ 2(M+−M−) + O(η^{1/3})' should be '≥'. The final conclusion (M+ negative and small) is correct.","section":"§5.1, after Proposition 5.3"},{"comment":"'Since q(0)=1' should read 'q(0)>0' (Prop. 5.9 gives q(0)=H(0)>0). The argument is unaffected.","section":"§5.2, proof of Proposition 5.10"},{"comment":"'The integral (5.3) is well-defined' — (5.3) is the list of limit properties; the reference should be to (4.9)/(4.38).","section":"§5, proof of Proposition 5.1"},{"comment":"The exclusion of eigenvalue sequences accumulating on the imaginary axis works via the analytic function F, but a sentence confirming that eigenvalues of A outside spec(A_0) are exactly the zeros of F−1/2 (and hence isolated) would close the argument, given the degeneracy H(1)=0 of the Sturm–Liouville weight.","section":"§5.2, Corollary 5.11"},{"comment":"R = 2(‖ũ_0‖_{H²} ∨ ẽ_0 ∨ m̃_0) should use |ẽ_0|, |m̃_0|, since these can be negative.","section":"Appendix A, Proposition A.2"},{"comment":"The normal switching condition (4.21) is said to hold 'if ε|p(y)| is small enough'; since a(E±,M∓) ≠ 0 by Assumption 2.4, please state the explicit smallness condition on ε/q_0.","section":"§4.2.2, point 2"},{"comment":"Suggest adding one sentence after Thm. 2.7 clarifying that the constants C, M_0, Δ_0 in (2.9) depend on the uniform lower bound q_0 in (4.19), and that persistence of the shock-free bound is proved only within the basin (forward reference to §5.3/Remark 6.1). The conditional logic is correct as stated, but easy to miss on first reading.","section":"§2.3"},{"comment":"Typos: 'Lipshitz' → 'Lipschitz' throughout (Thm. 2.1, §A); the abstract says 'a unique periodic solution' while Thm. 2.7 asserts an isolated one — align wording; Fig. 4 axis labels render as placeholder glyphs in the arXiv PDF; Δx, Δt are specified only for Figs. 6–7, please add them for Figs. 2–5.","section":"General"}],"recommendation":"minor_revision","confidential_remarks":"The manuscript is carefully written and self-contained, with an appropriate citation pattern and a clearly stated gap to the prior literature. The stress-test concern about the a-posteriori slope bound (4.19)/(4.23) is real but, on my reading, does not undermine Theorem 2.7, whose hypotheses are explicitly conditional; the cost is scoping of the basin, which the authors acknowledge. The two major comments ask for justification of steps where robustness is asserted rather than shown (Props. 4.5 and Remark 6.1); both look repairable with standard estimates. Numerics are illustrative and do not carry the claims. Well suited to the journal."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The new result is Theorem 2.7: small viscosity collapses an infinite-dimensional family of inviscid relaxation cycles (Prop 2.5) to one isolated, linearly stable periodic orbit for symmetric monotone data, with rate O(ε) and period converging to the inviscid one. That combination—non-propagating, globally mass-coupled viscous conservation law plus fast relaxation ODE—is not covered by classical viscous-front theory, periodic travelling-wave stability, or standard FHN slow–fast ODEs. The architecture is standard pieces newly assembled: reciprocal profile X, polar-like (r,θ) coordinates outside the central rectangle, GSPT radial control via Szmolyan–Wechselberger, cycle averaging to a stationary problem for G, conserved K fixing the unique averaged state, rank-one Sturm–Liouville spectral gap, then Poincaré + IFT through the O(ε², η^{1/3}ε) remainders.\n\nWhat works: global existence and symmetry are clean; the averaging and spectral analysis are careful; the hypotheses (S-shaped critical manifold, slow nullcline geometry) are structural rather than fitted; numerics (finite-domain Lax scheme) illustrate the selection, the O(ε) decay, and the canard-like transition when the slow nullcline nears a fold. The authors flag the open questions (non-symmetric data, large slope, possible longer-period orbits, MMOs) instead of hiding them.\n\nThe real soft spot is exactly the one the stress-test flags: every averaging/GSPT step needs ∂_m X(θ,0) ⩾ q_0 > 0 so shocks never form and normal switching stays valid. That bound is closed only a posteriori inside an H² ball around the selected orbit (Cor 5.14). Logically the theorem is fine—it assumes the neighbourhood (2.9) and verifies persistence—but the abstract’s “selects a unique” language is stronger than the proved basin. Large-slope or non-monotone data remain open, and Remark 4.7 correctly notes folded-node/MMO possibilities once the bound fails. Minor: numerics are code-free and finite-domain, so they support rather than prove.\n\nThis is for people working on slow–fast PDEs, viscous conservation laws with global coupling, or mathematical plasma models. The math is solid enough for a serious referee. I would send it to peer review; the scoped limitations are honest, not fatal.","headline":"Clean viscosity-selection theorem for a globally coupled conservation-law–relaxation oscillator; the a-posteriori slope bound scopes the basin but does not break the stated claim.","tokens_in":44715,"tokens_out":610,"would_cite":true,"duration_ms":17293,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q49","34C26","34K33","35P15"],"pacs":[],"model":"grok-4.5","headline":"Small diffusion collapses an infinite family of relaxation oscillations into one linearly stable periodic orbit.","keywords":["slow-fast system","relaxation oscillation","conservation law","averaging","spectral stability","viscous regularisation","transport-ODE coupling","plasma void"],"falsifier":"Simulate or analyse the system from a symmetric monotone initial profile whose midpoint slope is large enough that the viscous term cannot prevent steepening: if a shock forms or the orbit fails to approach the unique averaged steady state predicted by the theory, the selection-and-stability claim fails for that regime.","tokens_in":44134,"feed_emoji":"🌊","tokens_out":925,"duration_ms":17186,"temperature":0.7,"pith_summary":"A viscous conservation law for a particle density is coupled to a fast ODE that can sustain relaxation oscillations. When viscosity vanishes, every admissible initial density profile evolves into its own time-periodic solution, all sharing the same period set by the slow–fast ODE. The paper proves that a small positive viscosity, acting on symmetric monotone initial data, selects a single isolated periodic orbit and makes it linearly asymptotically stable, with attraction rate proportional to the viscosity. The argument rewrites the profile via its reciprocal function, averages over the relaxation cycle, and obtains a stationary averaged equation whose unique admissible steady state is shown stable by spectral methods. The result matters because it shows how weak diffusion can lift a massive degeneracy and pick a preferred oscillatory pattern in a transport–ODE model motivated by dusty-plasma voids.","feed_headline":"Tiny diffusion picks one stable oscillation from infinitely many","feed_subtitle":"A viscous conservation law coupled to a relaxation oscillator collapses a continuum of cycles to a single linearly stable orbit","key_machinery":"Averaging of the reciprocal-profile equation over one relaxation cycle. After a slow-time change of variables the leading-order averaged dynamics is autonomous; its unique stationary solution compatible with a conserved mass integral is linearly stable by a Sturm–Liouville spectral-gap argument on an invariant hyperplane, and a Poincaré-map/implicit-function argument lifts the conclusion back to the original system.","core_discovery":"Under the paper’s structural assumptions on the nullclines and for symmetric monotone initial data, sufficiently small positive viscosity ε (with time-scale ratio η at most order 1/log(1/ε)) produces an isolated time-periodic solution of the coupled system that is linearly asymptotically stable in H²×ℝ, attracts nearby solutions at rate O(ε), and whose period converges to the common inviscid period as ε→0. Without viscosity the same system admits an uncountable family of such periodic solutions, one for each initial profile.","pith_inferences":["The same averaging-plus-spectral strategy may apply to other scalar conservation laws globally coupled to a relaxation oscillator, not only the dusty-plasma idealisation.","Non-symmetric initial data appear numerically to select different limit cycles labelled by total mass; a full uniqueness theory without symmetry remains open.","If the midpoint slope is allowed to approach zero, folded-node or canard mechanisms could generate mixed-mode oscillations, consistent with the paper’s numerical canard-like transients."],"forward_implications":["Weak diffusion is enough to select a unique relaxation-oscillation density profile among the continuum present at zero viscosity.","The attraction rate to that orbit scales linearly with the diffusion coefficient ε.","The selected period converges to the period of the inviscid slow–fast limit cycle as ε→0.","Symmetric monotone data near the orbit remain strictly increasing for all time, so the reciprocal-profile representation never breaks.","Changing the relative placement of the slow and fast nullclines can replace large-amplitude relaxation cycles by small-amplitude or mixed-mode patterns."],"fun_headline_variants":["Small viscosity selects one stable cycle from infinitely many","Diffusion isolates a unique linearly stable periodic orbit","Vanishing viscosity picks single stable oscillation via averaging","Positive viscosity collapses continuum of cycles to one stable solution","Tiny diffusion stabilises unique time-periodic solution"],"cache_read_input_tokens":32896,"weakest_assumption_plain":"The inverse density profile must keep a strictly positive slope at the midpoint for all time, so that no shock forms and the averaging hypotheses stay valid; this is only guaranteed when the initial data are already close to the selected orbit.","fun_headline_variants_meta":{"raw":{"variants":["Small viscosity selects one stable cycle from infinitely many","Diffusion isolates a unique linearly stable periodic orbit","Vanishing viscosity picks single stable oscillation via averaging","Positive viscosity collapses continuum of cycles to one stable solution","Tiny diffusion stabilises unique time-periodic solution"]},"model":"grok-4.5","effort":"low","cost_usd":0.003206,"raw_usage":{"total_tokens":1008,"prompt_tokens":656,"num_sources_used":0,"completion_tokens":58,"cost_in_usd_ticks":32064000,"prompt_tokens_details":{"text_tokens":656,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":294,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":656,"tokens_out":58,"duration_ms":5190,"temperature":1.0,"reasoning_tokens":294,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T04:05:49.148008+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Simulate or analyse the system from a symmetric monotone initial profile whose midpoint slope is large enough that the viscous term cannot prevent steepening: if a shock forms or the orbit fails to approach the unique averaged steady state predicted by the theory, the selection-and-stability claim fails for that regime.","supporting_citations":[],"review_version":1}