{"id":"afef6840-6b04-4068-ad2c-36bb1c71a0fb","arxiv_id":"1908.08716","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A linear PDE with a piecewise-constant moving coefficient reproduces several short-time features of KdV dispersive shocks, while its long-time behavior differs.","lead":"The authors present a linear wave equation with a moving step that mimics the oscillating front, or dispersive shock, normally associated with the nonlinear Korteweg-de Vries equation. A generalist might read it because it shows some shock-like features do not require nonlinearity, and a simple linear model can approximate them at short times.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"All numerical evidence for the claimed shock features depends on the unproved UTM representation in Proposition 3.1, with unspecified quadrature; this gap leaves the central claim conditional.","rationale":"The reader's weakest assumption matches my own reading: Proposition 3.1 is the load-bearing tool, and it is unproved, with the accompanying numerical methodology underspecified. I read the full text and found no code, no data, and no error estimates for the UTM evaluation; the paper itself flags that the method is not asymptotically accurate and that the proof is deferred. The central claim is plausible and the model is interesting, so I do not see a ground for rejection, but the evidence as presented is not independently verifiable. The c=4a choice is a genuine interpretive weakness—the model is given the nonlinear shock speed rather than deriving it—but even if that were addressed, the plotted comparisons still depend on the correctness of the UTM representation. Therefore the appropriate verdict remains conditional: the qualitative claims should be accepted only after Proposition 3.1 is proved or independently verified numerically. I am not raising an objection based on disagreement with prevailing views; the concern is internal to the paper's support for its own claim, and a direct numerical comparison would settle it.","tokens_in":5804,"tokens_out":7915,"duration_ms":78333,"concrete_test":"Run an independent direct numerical simulation of (4) with the same step data and interface conditions (4e), using a Chebyshev collocation discretization in x with an interface treatment at 0 and a high-order time integrator, at the (x,t) values of Figures 1–4 and for a=1, c=4; also recompute E_model and E_LKdV at t=0.1 for a=1, 1/2, 1/4. Compare pointwise against the UTM-based values. If the maximum discrepancy is below plotting tolerance (about 1e-3) and converges under mesh refinement, Proposition 3.1 is validated for the plotted regime; if not, the numerical claims need revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that model (3) captures KdV dispersive-shock features and is a better short-time linearization than the naive one—is supported entirely by Figures 1–5, all computed from Proposition 3.1. That proposition is stated without proof; the only justification is 'The proof of this proposition may be obtained by following the method of [6], suitably modified in light of the results of [4]' (§3). The numerical evaluation is also not reproducible: the text says the methodology 'is not asymptotically accurate' and provides no contours, quadrature rules, or code. If the contour deformation is wrong for the plotted times, or if the τ-dependence in h and g_j is mishandled, the displayed oscillations and the error comparison in Figure 5 could be artifacts. The claim in §5 that 'our results are not contingent on the unified transform method' does not remove this dependency, since every displayed q is produced by that representation. This is a verification gap in the evidence for the main claim, not a disagreement with consensus; supplying a proof and an independent numerical check would close it. The separate issue that c=4a is chosen to match the known KdV shock speed is real but secondary: even if c were independently selected, the plotted evidence would still require validating the UTM formula.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a linear dispersive PDE with a moving-interface coefficient (equation (3)) that is intended to model early-time dispersive shock behavior of the Korteweg-de Vries equation with step-like initial data. The authors claim that, by linearizing about the step initial condition rather than about zero, the model reproduces several qualitative features of KdV dispersive shocks: growing oscillation amplitude, frequency increase with amplitude, transient linear peak arrangement, and a low-amplitude error smaller than that of the naive linearization. The model is solved formally using the unified transform method (UTM), leading to an explicit integral representation (Proposition 3.1). The paper compares this solution against numerics for KdV and against the stationary long-time limit, concluding that the linear model is a better short-time linearization than linearizing about zero while deviating at long times.","tokens_in":6080,"tokens_out":1711,"duration_ms":18057,"significance":"If the central claim is correct, the paper offers a conceptually interesting observation: a genuinely linear model with a simple moving interface can mimic several hallmarks of nonlinear dispersive shocks, which are usually regarded as intrinsically nonlinear phenomena. The idea of linearizing about the initial condition to obtain a piecewise-constant linear operator is suggestive and could be useful for cheap short-time estimates in other dispersive equations. The paper also includes explicit long-time asymptotics for the model and a comparison with a numerical KdV solver. However, the significance is currently conditional because the main supporting evidence is numerical and depends on an unproved solution representation, and the one free parameter (front speed c) is set to the known KdV shock speed in the comparisons.","major_comments":[{"comment":"Proposition 3.1, which provides the explicit solution representation used for every figure in the paper, is stated without proof. The only justification is the sentence: 'The proof of this proposition may be obtained by following the method of [6], suitably modified in light of the results of [4].' Since all qualitative claims—oscillation growth, peak arrangement, and the error comparison in Figure 5—are computed from this representation, the central claim of the paper is unsupported unless this proposition is proved or at least verified independently. The authors should include a proof or a precise citation to a fully worked argument, and ideally an independent numerical check (e.g., direct time-stepping of the linear PDE (3)) for the plotted parameter values.","section":"§3, Proposition 3.1"},{"comment":"The numerical evaluation of the UTM integrals is not reproducible from the text. The paper states that the methodology 'is not asymptotically accurate' but does not specify the contour deformation, the quadrature rules, or the discretization parameters used to produce Figures 1–5. Without this information, the reader cannot assess whether the displayed oscillations, the transient peak alignment, or the error comparison are genuine properties of the model or numerical artifacts. The authors should provide full numerical details (contours, tolerance, number of quadrature points) or release code, and ideally compare with a direct numerical solution of (3) at a few (x,t) values.","section":"§4 and Figures 1–5"},{"comment":"The short-time analysis is heuristic. The paper argues that analyzing U(x,t) = Q(x - gamma t, t) at fixed t for large c is equivalent to small time because of the scaling t -> c^{-3/2}t, but no rigorous asymptotic statement is made, and no error bounds are provided. Since the claim that the linear model is a 'better short-time linearization' is one of the central conclusions, the authors should either state a rigorous asymptotic result (with error estimates) or clearly label this as a numerical observation and support it with quantitative convergence data.","section":"§4.1"},{"comment":"The front speed c is introduced as a free parameter but is set to c = 4a in the comparisons with KdV, which is precisely the known speed of the KdV dispersive shock. This means the model is given the velocity of the phenomenon it is claimed to reproduce. The other qualitative features (oscillations, amplitude growth, peak arrangement) are not forced by this choice, so the issue is secondary, but it should be addressed explicitly: the authors should either justify c = 4a on independent grounds or demonstrate that the qualitative agreement is robust to variations in c over a range of values.","section":"§2 and §4.4"}],"minor_comments":[{"comment":"The paper contains several typos and minor notation inconsistencies, including 'In Figure 1a we plot' while the figure caption refers to 'A plot of U(x, 0.1)' with no subfigure label, and the use of 'ν' versus 'nu' in the displayed system for the functions g_j in Section 3. A careful proofreading is recommended.","section":"General"},{"comment":"The stationary solution formulas in Section 4.3 are not fully derived; in particular, the constants c1, c2, b1, b2, b3 are stated without explanation. Adding a brief derivation or a reference would improve clarity.","section":"§4.3"},{"comment":"In the small-amplitude discussion, the error comparison in Figure 5 is presented only at a single time t = 0.1 and without error bars or convergence data. A short discussion of the dependence on t and on the numerical resolution of u_KdV would make the comparison more convincing.","section":"§4.4"},{"comment":"The concluding sentence 'Our results are not contingent on the unified transform method' is misleading: while the model itself is independent of UTM, every displayed solution and all numerical evidence in the paper rely on the UTM representation of Proposition 3.1. The authors should clarify that the qualitative conclusions are not contingent on UTM as a method, but the numerical verification is.","section":"§5"}],"recommendation":"major_revision","confidential_remarks":"The central idea is original and the numerical evidence, if validated, would be of interest to the dispersive-shock community. However, the unproved Proposition 3.1 is the sole basis for all figures, and the paper currently does not provide enough numerical detail for reproducibility. This is fixable within the manuscript's scope by adding a proof (or a precise reference to a fully worked derivation), an independent numerical check, and full numerical parameters. I recommend major revision rather than rejection because the core qualitative claim appears plausible and the missing pieces are within reach. I would also suggest the editors ensure that the comparison with c=4a is explicitly framed as using the known KdV shock speed, since this is a potential circularity that the authors should address head-on."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a paper with a genuinely nice idea and an honest statement of its limits, but the main technical engine is asserted rather than proved. The model—linear KdV with a piecewise constant coefficient that jumps at x = ct—does mimic several short-time features of the KdV dispersive shock: oscillation amplitude grows, frequency increases with a, and the peak arrangement is transiently linear. The authors also correctly say the long-time behavior is wrong, so they are not overclaiming.\n\nWhat's new here is the move to linearize about the initial step rather than about zero. That turns the problem into an interface problem for a third-order linear PDE, and the UTM representation in Prop 3.1 is apparently new for this equation. The rescaling that reduces a and c to the single ratio a/c is neat, and the cautionary remark about comparing nonlinear theory to experiment is fair and useful.\n\nThe soft spots are real and load-bearing. Prop 3.1 is stated without proof; the sentence 'may be obtained by following the method of [6]' is not a proof, and every figure in the paper is computed from that representation. The numerical evaluation is described as 'not asymptotically accurate,' but no contours, quadrature rules, or code are given. So the central qualitative claim is conditional on a formula we cannot check. The choice c = 4a is an understandable way to match the known KdV shock speed, but it means the front speed is not emergent; that is secondary, since the oscillation structure is not forced by that choice. The 'better than naive linearization' conclusion rests on Figure 5, which compares errors at a single time t = 0.1 for three amplitudes, with no error bars. It is suggestive, not conclusive.\n\nI think this deserves a serious referee, not a desk reject. A referee should ask for a derivation of Prop 3.1 (or a pointer to a paper with it), for reproducible numerics, and for a more systematic error study. If those are supplied, this becomes a useful short-time surrogate for KdV shocks. I wouldn't cite it yet for my own work, but I'd follow what comes out of revision.","headline":"A clever linear surrogate for short-time KdV dispersive shocks, but the central formula is unproved and the numerics unreproducible; needs a serious referee, not a desk reject.","tokens_in":6578,"tokens_out":2640,"would_cite":false,"duration_ms":26519,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q53","35C15","35B40"],"pacs":[],"model":"deepseek-v4-flash","headline":"A moving-jump linear model reproduces KdV dispersive-shock features at short times.","keywords":["dispersive shocks","Korteweg-de Vries equation","linearization about initial data","moving interface problem","unified transform method","step initial data","short-time asymptotics","long-time behavior"],"falsifier":"At $a=1$, $c=4$, solve the moving-interface problem (3)-(4) with an independent numerical method, such as a spectral or finite-difference scheme that enforces continuity of $q,q_x,q_{xx}$ at $x=0$, and compare with the plotted profiles at $t=0.1,0.5,1.5,2.75,4,8,20$: the reported growth of the maximum, the near-linear alignment of the first three peaks at $t=2.75$, and the approach to the stationary solution at $t=20$ are directly checkable. The representation in Proposition 3.1 can also be tested by substituting the integral expressions into the PDE and interface conditions.","tokens_in":5583,"feed_emoji":"🌊","tokens_out":12429,"duration_ms":116493,"temperature":0.7,"pith_summary":"The paper claims that a dispersive shock, normally viewed as a nonlinear phenomenon, can be reproduced at early times by a deliberately chosen linear PDE. The construction replaces the KdV nonlinear term $6u u_x$ by a moving jump in the coefficient of $u_x$, freezing the initial profile's left and right amplitudes $6a$ and $0$ on either side of the front $x=ct$. Using the unified transform method, the authors obtain an explicit integral representation for this moving-interface problem and show numerically that its maximum grows in time, its oscillation frequency increases with jump amplitude, and its early peaks line up almost linearly, all features of the KdV dispersive shock. They also show that this linear model is a better short-time approximation to KdV than linearizing about zero, while its long-time stationary profile is explicitly different from KdV's. If the claim is right, it provides a cheap linear route to short-time shock estimates and a caution that early agreement between a linear model and nonlinear observations does not by itself confirm nonlinear mechanisms.","feed_headline":"A linear PDE mimics a KdV dispersive shock at short times","feed_subtitle":"Freezing the nonlinearity into a moving step recreates shock oscillations and beats linearizing about zero.","key_machinery":"The load-bearing object is the piecewise-constant linear differential operator obtained by freezing the nonlinear coefficient in the KdV equation at its initial left/right values across a moving front. In the traveling frame $q(x,t)=u(x+ct,t)$, the model is $q_t+q_{xxx}=q_x(c-6a\\,\\mathbf{1}_{x<0})$ with continuity of $q,q_x,q_{xx}$ at $x=0$. The solution is represented by the unified transform method, a contour-integral representation for linear PDEs with piecewise-constant coefficients: Proposition 3.1 expresses $q$ as integrals over the boundary of a sector $\\mathbb{D}$ in the spectral $\\lambda$-plane, involving the three roots $\\nu_0,\\nu_1,\\nu_2$ of the dispersion relations $\\nu^3+c\\nu=\\lambda^3$ and $\\nu^3+(c-6a)\\nu=\\lambda^3$. Deforming these contours into regions of exponential decay and applying standard quadrature gives the high-accuracy profiles plotted in the paper. This representation is also what allows the authors to extract the long-time stationary behavior by solving the ordinary differential equation $q_{xxx}=c q_x$ on each side.","core_discovery":"On its own terms, the paper's discovery is that the moving-interface linear problem\n$$u_t+u_{xxx}=-6a\\,\\mathbf{1}_{x<ct}\\,u_x,\\qquad u(x,0)=a\\,\\mathbf{1}_{x<0},$$\nwith continuity of $u,u_x,u_{xx}$ at $x=ct$, shares the short-time qualitative signatures of the KdV dispersive shock: growing oscillation amplitude, higher frequencies for larger $a$, and a transient near-linear arrangement of the leading peaks. The authors obtain the needed solution formula (Proposition 3.1) by adapting the unified transform method for interface problems, use it to compute high-accuracy profiles, and compare errors against the KdV solution computed by numerical inverse scattering. Their error measure shows the model dominates the naive linearization $F=0$ in the small-amplitude limit, with errors that decay in $|x|$. At long times the model tends to a stationary solution with oscillations on the left side of the front, which is not the KdV long-time behavior; the paper presents this difference explicitly as the boundary of the model's validity.","pith_inferences":["The same freezing recipe should transfer to other dispersive equations whose nonlinear term acts as a coefficient of $u_x$; replacing that coefficient by its left/right limiting values across a moving front should produce the same short-time oscillation growth.","A practical extension is iterative re-linearization: solve the linear interface problem up to a positive time, linearize about that newer profile, and solve again; this could yield numerical schemes that track dispersion better, though it would require a fast multi-interface solver.","Because the model reaches a stationary oscillatory state on the left at long times, an interesting test is whether a similar stationary pattern appears in KdV at intermediate times before solitons dominate; the paper does not pursue this."],"forward_implications":["For step initial data, the linear model's maximum profile amplitude grows for small times, a signature the authors trace to the nonlinearity in KdV and show is preserved in the linear moving-interface model.","Increasing the jump amplitude $a$ increases the frequency of the oscillations in the model, matching the amplitude-frequency relation of KdV dispersive shocks.","After rescaling, the model is the universal small-amplitude limit for general initial data of the form $a(v(x)+h(x))$, so integrable corrections to the step contribute only lower-order effects.","In the small-amplitude limit the error between the model and KdV is $O(a^{-1})$, with a smaller constant than for the naive linearization, and the error decays for large $|x|$.","Because the model is faithful at short times but not at long times, comparisons of nonlinear PDE predictions with experiments should use sufficiently long times for the nonlinearity to manifest; the paper states this as a cautionary consequence."],"supporting_citations":[{"why":"Supplies the unified transform method for interface problems that Proposition 3.1 generalizes to the moving-front KdV-type equation.","marker":"[6]"},{"why":"Extends the interface method to the linear KdV equation; the proof of Proposition 3.1 is said to follow from [6] suitably modified in light of this work.","marker":"[4]"},{"why":"Provides the numerical inverse-scattering computation of the KdV solution with step-like data used for the error comparison in Figure 5.","marker":"[2]"},{"why":"Establishes the KdV dispersive-shock phenomena, including front speed $4a$ and faster solitons, that the linear model is designed to reproduce.","marker":"[1]"},{"why":"Gives the long-time asymptotics of KdV with step-like initial data, the reference point for the paper's statement that its linear model deviates at long times.","marker":"[5]"},{"why":"Provides the similarity solution of the naive linearization about zero, the baseline whose time-independent maximum amplitude the paper contrasts with its model's growing amplitude.","marker":"[3]"}],"fun_headline_variants":["Linear step PDE mimics KdV shock at short times","Moving-interface linear model reproduces KdV shock","Short-time shock mimicry without nonlinearity","Linear trick: KdV-style shock from a step","Linear PDE briefly echoes KdV dispersive shock"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"All the plotted conclusions rest on the stated solution formula for the moving-interface problem, which the paper gives without proof; if that formula is not correct or the numerical evaluation of its integrals is not accurate at the times shown, the claimed reproduction of shock features is unsupported.","fun_headline_variants_meta":{"raw":{"variants":["Linear step PDE mimics KdV shock at short times","Moving-interface linear model reproduces KdV shock","Short-time shock mimicry without nonlinearity","Linear trick: KdV-style shock from a step","Linear PDE briefly echoes KdV dispersive shock"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000153,"raw_usage":{"total_tokens":1149,"prompt_tokens":827,"completion_tokens":322,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":443,"completion_tokens_details":{"reasoning_tokens":249}},"tokens_in":443,"tokens_out":322,"duration_ms":3994,"temperature":1.0,"reasoning_tokens":249,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:31:40.657431+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"At $a=1$, $c=4$, solve the moving-interface problem (3)-(4) with an independent numerical method, such as a spectral or finite-difference scheme that enforces continuity of $q,q_x,q_{xx}$ at $x=0$, and compare with the plotted profiles at $t=0.1,0.5,1.5,2.75,4,8,20$: the reported growth of the maximum, the near-linear alignment of the first three peaks at $t=2.75$, and the approach to the stationary solution at $t=20$ are directly checkable. The representation in Proposition 3.1 can also be tested by substituting the integral expressions into the PDE and interface conditions.","supporting_citations":[{"cited_title":"3, 253–275","cited_arxiv_id":null,"evidence_quote":"Supplies the unified transform method for interface problems that Proposition 3.1 generalizes to the moving-front KdV-type equation."},{"cited_title":"2, 489–509","cited_arxiv_id":null,"evidence_quote":"Extends the interface method to the linear KdV equation; the proof of Proposition 3.1 is said to follow from [6] suitably modified in light of this work."},{"cited_title":"On numerical inverse scattering for the Korteweg-de Vries equation with discontinuous step-like data","cited_arxiv_id":"1809.09263","evidence_quote":"Provides the numerical inverse-scattering computation of the KdV solution with step-like data used for the error comparison in Figure 5."},{"cited_title":"2, 022906","cited_arxiv_id":null,"evidence_quote":"Establishes the KdV dispersive-shock phenomena, including front speed $4a$ and faster solitons, that the linear model is designed to reproduce."},{"cited_title":"7, 1839–1864","cited_arxiv_id":null,"evidence_quote":"Gives the long-time asymptotics of KdV with step-like initial data, the reference point for the paper's statement that its linear model deviates at long times."},{"cited_title":"3, 813–837","cited_arxiv_id":null,"evidence_quote":"Provides the similarity solution of the naive linearization about zero, the baseline whose time-independent maximum amplitude the paper contrasts with its model's growing amplitude."}],"review_version":1}