{"id":"5612e440-7da4-4e87-af94-3e966f309d62","arxiv_id":"2507.07571","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Two generalized KdV/BBM-type continuum models reproduce the dispersive shock waves of the granular chain, including in the zero-precompression sonic vacuum regime.","lead":"This paper derives two new first-order wave equations that approximate the motion of a chain of touching beads, and shows they capture shock-like waves of the discrete chain even when the chain is not pre-compressed. The models give researchers a simpler, more analytically tractable way to study these nonlinear waves.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Zero-precompression validity rests on an extrapolation: at r_+=0 the linearized dispersion used to select the unidirectional branch vanishes identically, and the only supporting numerics modify the initial data and use r_+=1e-5 for the non-regularized model.","rationale":"The reader identified the singular zero-precompression limit of the unidirectional reduction as the weakest assumption; the full text supports this. Section 2.1 builds (2.10) from a positive branch whose coefficient A^{(p-1)/2} vanishes at A=0, so the first-order character of the model is not derived in the sonic vacuum. Section 6.2 also shows the practical consequences: the non-regularized model was not run at r_+=0 but at r_+=10^{-5}, and the initial condition was smoothed with delta=1 rather than delta=50; the regularized model at exact zero shows trailing-edge deviations. These facts weaken the abstract's unconditional zero-precompression statement. I nevertheless agree with the conditional verdict: the precompressed case is carefully derived and validated, the derivation of the regularized model is a reasonable BBM-type analogue, and the zero-precompression failure would weaken but not invalidate the paper's core modeling contribution. No independent verification or code is supplied, so the concern cannot be dismissed without the proposed computational check.","tokens_in":22933,"tokens_out":5313,"duration_ms":63553,"concrete_test":"Fix the Riemann data as in Fig. 10 (r_-=0.05, r_+=0, p=3/2, delta=50) and integrate the non-regularized model (2.12) at exact zero background with a positivity-preserving conservative scheme; measure leading/trailing edge speeds and trailing wavenumber using the same Appendix procedure, and compare with the discrete chain. If the exact zero-background run cannot be performed with sharp data, or if the measured edges move outside the scatter of Fig. 10, the zero-precompression clause of the central claim is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires models (2.12) and (2.14) to describe DSWs of the discrete chain (2.2) at r_+=0. The derivation selects the positive branch of the linearized dispersion relation (2.10)/(2.11), whose prefactor is sqrt(p) A^{(p-1)/2}. At A=r_+=0 the linearized dispersion (2.3) is identically zero, phase/group velocities vanish, and there is no distinguished right-going wave direction; a first-order unidirectional reduction is therefore not a justified limit but a formal extrapolation. Eq. (2.12) itself becomes degenerate at r=0: the nonlinear flux derivative at zero vanishes, so the model's dynamics near the zero background rest entirely on finite-amplitude regions where the linearization used to choose the branch does not apply. The numerical support in Section 6.2 is not decisive: the non-regularized model was integrated with r_+=10^{-5} and with a smoother initial transition (delta=1 instead of delta=50), and the regularized model at r_+=0 shows visible trailing-edge deviations from the lattice DSW in Fig. 10, with good profile agreement only for the small jump r_-=0.05. Thus the abstract's 'good agreement ... even in cases where no precompression is present' is stronger than the evidence presented, and the central claim is load-bearing on an unverified singular limit.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives two first-order continuum models for the granular crystal lattice (2.2): a generalized KdV-type model (2.12) and a regularized BBM-type model (2.14). It analyzes their solitary and periodic traveling waves, derives conservation laws and Whitham modulation systems, applies the DSW fitting method to obtain leading- and trailing-edge predictions, and compares those predictions with direct numerical simulations of the discrete lattice. The principal claims are that the new models outperform the classical KdV approximation for small precompression and remain accurate even in the zero-precompression (sonic vacuum) limit.","tokens_in":23224,"tokens_out":4538,"duration_ms":49177,"significance":"If the claims hold, the paper provides a useful intermediate-level unidirectional description of a canonical nonlinear lattice, with parameter-free DSW edge predictions obtained from an established fitting method. The derivation from the discrete equations is transparent, the comparison against the KdV benchmark is appropriate, and the precompression results are convincing. The main weakness is that the zero-precompression claim, which is central to the abstract, rests on a singular limit and on numerical evidence that the authors themselves describe as showing deviations; this part of the claim needs to be either strengthened or carefully qualified.","major_comments":[{"comment":"The unidirectional reduction selects the positive branch of the linearized dispersion relation (2.10), whose prefactor sqrt(p) A^{(p-1)/2} vanishes at A = r_+ = 0. At r_+ = 0 the linearized dispersion relation (2.3) is identically zero, so there is no distinguished right-going wave direction and the KdV-type scaling used for the reduction becomes singular. Applying (2.12) and (2.14) at r_+ = 0 is therefore a formal extrapolation rather than a justified limit. This is load-bearing because the abstract's central claim concerns exactly this case; please either supply an explicit asymptotic or numerical justification for the limit, or restrict the claim to finite precompression plus a clearly labeled empirical extrapolation.","section":"§2.1, Eqs. (2.10)–(2.12)"},{"comment":"The numerical support for the zero-precompression claim is partial. The text states that the trailing-edge DSW features of both continuum models deviate from those of the granular lattice, that solitonic amplitudes deviate as r_- increases, and that agreement is good only for the small jump r_- = 0.05. In addition, the non-regularized model required r_+ = 10^{-5} and δ = 1 instead of the standard δ = 50. The abstract's unconditional statement of 'good agreement ... even in cases where no precompression is present' is therefore stronger than the evidence presented. Please add quantitative error metrics for the edge quantities and revise the summary claims to match the demonstrated regime.","section":"§6.2, Figs. 10–11"},{"comment":"The DSW fitting formulas for the discrete lattice are derived from the linearized dispersion relation around r_+ and contain factors such as (r_-/r_+)^{(p-1)/4}, which diverge as r_+ → 0; the paper acknowledges in §6.2 that the fitting formulas are invalid at r_+ = 0. Consequently, the zero-precompression comparisons rely entirely on direct numerical simulation of the continuum PDEs, without an independent theoretical prediction for the discrete chain. This limitation should be stated in the abstract and conclusions, or a separate zero-precompression prediction should be supplied.","section":"§5.3, Eqs. (5.15)–(5.16)"}],"minor_comments":[{"comment":"The formula for the p = 3/2 solitary wave of the non-regularized model is garbled: the argument of the sine contains malformed square-root factors. Please correct the typesetting so that the expression can be verified against Eq. (3.5c).","section":"Eq. (3.10)"},{"comment":"References [18] and [33] cite the same paper (Chong, Geisler, Kevrekidis, Biondini, Wave Motion 130, 2024) and should be consolidated to avoid duplicate citation.","section":"References"},{"comment":"The use of K_m for the complete elliptic integral of the first kind alongside K for the wavenumber is potentially confusing even though it is explained in the text; consider using a different symbol such as mathcal{K}(m).","section":"§4.3, Eq. (4.15)"},{"comment":"The sentence reporting numerical instability in the non-regularized model would benefit from stating whether the instability occurs before or after the DSW forms, and whether the choice δ = 1 and r_+ = 10^{-5} changes the measured edge features in a quantifiable way.","section":"§6.2, first paragraph"}],"recommendation":"major_revision","confidential_remarks":"No concerns about novelty or scope. The main issue is the overstatement of the zero-precompression validity in the abstract and conclusions relative to the evidence in Section 6.2; my assessment is that this can be addressed by a careful revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nRead this one. The genuinely new thing is the pair of first-order continuum models, (2.12) and (2.14): a generalized KdV and a BBM-type regularization, derived from the granular strain equations rather than fitted. That is a real step beyond the standard KdV reduction, which requires precompression, and beyond the authors' earlier bidirectional second-order model. The linear dispersion comparison, solitary and periodic wave analysis, conservation laws, and Whitham modulation systems are worked out in a mostly clean, standard way. In the precompressed regime, the DSW edge predictions agree well with the discrete lattice, and they improve on KdV at small precompression. That part holds up.\n\nThe soft spot is exactly where the abstract is strongest. The derivation selects the positive branch of a linearized dispersion relation whose prefactor is sqrt(p) A^{(p-1)/2}; at A = r_+ = 0 that linearization vanishes identically, so the unidirectional reduction is an extrapolation, not a justified limit. The numerics in Section 6.2 reflect this: the non-regularized model is run with r_+ = 1e-5 and a smoother initial transition (delta = 1), and the regularized model at r_+ = 0 shows visible trailing-edge deviations except for the smallest jump. The paper itself admits the trailing-edge features deviate, and the conclusions say \"reasonable\" rather than \"good.\" So the abstract overstates the zero-precompression evidence. This is a fixable overclaim, not a load-bearing crack in the precompression story. I would not call the central idea invalid; I would call the sonic-vacuum claim a formal extrapolation supported by partial numerics.\n\nAlso worth noting: well-posedness of (2.12) is not established, and the Whitham modulation system is justified only formally. The authors say so themselves, so those are acknowledged limitations, not hidden flaws. I did not find invented entities or fitted simulations passing as predictions, and the citation pattern looks appropriate.\n\nWho is this for? Anyone working on granular crystals, discrete dispersive hydrodynamics, or DSW fitting in lattice models. The models give an analytically tractable intermediate description that KdV cannot provide near the sonic vacuum. I would send this to a serious referee, with the request that the authors either soften the abstract and conclusions or provide additional zero-precompression evidence; code and data would help.","headline":"New first-order generalized KdV and BBM-type models for the granular chain are a real step beyond the KdV reduction; the zero-precompression claim is an extrapolation that should be softened, but the paper deserves serious refereeing.","tokens_in":23780,"tokens_out":2369,"would_cite":true,"duration_ms":27215,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q53","35L67","37K40"],"pacs":[],"model":"deepseek-v4-flash","headline":"Two first-order continuum PDEs can stand in for the discrete granular crystal lattice and reproduce its solitary waves, periodic waves, and dispersive shock waves even without precompression.","keywords":["granular crystals","dispersive shock waves","continuum models","generalized KdV equation","BBM regularization","Whitham modulation theory","sonic vacuum","nonlinear lattices"],"falsifier":"Run the discrete lattice and both continuum PDEs with exactly zero precompression and a fixed jump, then compare the measured trailing-edge wavenumber and edge speeds; if the PDE values do not approach the lattice values as the numerical smoothing and the artificial $r_+=10^{-5}$ background are removed, the claimed sonic-vacuum validity fails.","tokens_in":22728,"feed_emoji":"🌊","tokens_out":11452,"duration_ms":118075,"temperature":0.7,"pith_summary":"The paper aims to show that two first-order-in-time continuum PDEs can stand in for the discrete granular crystal chain and reproduce its nonlinear wave behavior, including dispersive shock waves. The first is a generalized Korteweg–de Vries equation, and the second is its BBM-style regularization; both are derived from the positive branch of the chain's linearized dispersion relation. Using conservation laws, Whitham modulation theory, and the DSW fitting method, the authors obtain explicit predictions for the edges of dispersive shock waves and check them against direct simulations of the lattice. The agreement holds even in the no-precompression limit where the lattice is purely nonlinear and the classical KdV reduction has no linear dispersion to build on. A sympathetic reader would care because these PDEs are far more tractable analytically than the discrete lattice and could support rigorous existence and stability results.","feed_headline":"Two continuum models match granular crystal shock waves","feed_subtitle":"The agreement survives the purely nonlinear, zero-precompression limit where the KdV approximation fails.","key_machinery":"The load-bearing object is the positive branch of the linearized dispersion relation, $\\Omega = \\sqrt{p}\\,A^{(p-1)/2}K\\sqrt{1-\\epsilon^2K^2/12}$, expanded in its long-wave form (2.11), which converts the ill-posed second-order-in-time PDE into the first-order generalized KdV equation (2.12). Inverting the operator $1+\\epsilon^2\\partial_X^2/24$ regularizes the model into the BBM-type equation (2.14), whose dispersion relation stays bounded for large wavenumbers. This pair of models carries the argument because their solitary and periodic traveling waves, two conservation laws, and averaged Whitham modulation equations can be written explicitly enough to run the DSW fitting procedure and to predict how the leading and trailing edges of a dispersive shock move.","core_discovery":"On the paper's own terms, the central discovery is that the strain dynamics of the granular chain, $\\ddot r_n = (r_{n+1})^p - 2(r_n)^p + (r_{n-1})^p$, admit two unidirectional continuum limits that keep dispersive effects: the generalized KdV model (2.12) and the regularized BBM analogue (2.14). Both models come from taking the positive branch of the linearized dispersion relation, so they describe right-going waves, and both have enough conservation laws to close a Whitham modulation system by averaging. The DSW fitting procedure then yields explicit formulas for the trailing-edge wavenumber and the leading- and trailing-edge speeds, namely (5.10), (5.12), and (5.16). Numerical comparison shows that the PDE predictions agree with the discrete lattice for precompressed chains and continue to give a reasonable description of the DSW spatial profile at zero precompression, where the KdV model is not available.","pith_inferences":["A testable extension left implicit in the paper is to use these PDEs as design tools for laboratory granular-chain experiments at zero precompression, where the predicted edge speeds and wavenumbers could be measured directly.","The same two-step route—selecting the positive branch of the linearized dispersion relation and then regularizing—should transfer to other power-law lattices, including dimers, decorated chains, or two-dimensional packings, giving first-order continuum models for settings where no linear dispersion exists.","The artificial $r_+=10^{-5}$ required by the non-regularized model in the sonic-vacuum limit suggests a genuine singular limit; a matched-asymptotics analysis as $r_+\\to 0$ could show whether the DSW edge quantities obey universal power laws shared by the lattice and the PDEs.","If the models hold up, equation (5.2) gives a concrete way to prepare lattice initial data so that a laboratory shock experiment follows the continuum prediction from the earliest times."],"forward_implications":["For granular chains with zero precompression, where the KdV reduction is unavailable, the two PDEs provide quantitative predictions for DSW edge speeds, amplitudes, and wavenumbers.","For finite but small precompression, the new models approximate the discrete chain noticeably better than KdV; as precompression grows, all models become comparable.","The explicit edge formulas from DSW fitting, such as $k_- = 4\\,\\mathrm{arcsec}[(r_-/r_+)^{(p-1)/4}]$ for the lattice, can be used to initialize or interpret numerical and experimental Riemann problems.","The regularized model has a bounded dispersion relation and runs stably at $r_+=0$, while the non-regularized model needs a small positive background and smoother initial data; this distinction matters for choosing a model in practice.","The Whitham modulation systems derived here open a route to studying hyperbolicity, genuine nonlinearity, and rarefaction-wave structure in the continuum descriptions of lattices."],"supporting_citations":[{"why":"supplies the granular-chain setting and the sonic-vacuum concept that the new models are designed to cover","marker":"[1]"},{"why":"provides the Whitham modulation equations for the discrete FPUT-type chain used for lattice DSW fitting","marker":"[16]"},{"why":"supplies the DSW fitting method from which the edge-speed and wavenumber predictions are derived","marker":"[20]"},{"why":"is the earlier second-order regularized continuum model that the present first-order models extend and compare with","marker":"[22]"},{"why":"provides the BBM regularization idea used to obtain the regularized model (2.14)","marker":"[32]"},{"why":"gives the KdV-based DSW approximation used as the benchmark that the new models must beat","marker":"[33]"},{"why":"provides the KdV dispersive-shock edge formulas that the paper rescales and compares in the precompressed case","marker":"[34]"}],"fun_headline_variants":["Two continuum models match granular shocks","Continuum pair captures granular shock waves","Granular shock waves fit two continuum models","Beyond KdV: models match granular DSWs","Two models describe granular shock waves"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a one-directional, right-going wave description remains valid when the chain has zero precompression, even though the derivation expands around a nonzero background and the linear wave speed vanishes in that limit; the paper relies on numerical agreement there rather than on a controlled asymptotic argument.","fun_headline_variants_meta":{"raw":{"variants":["Two continuum models match granular shocks","Continuum pair captures granular shock waves","Granular shock waves fit two continuum models","Beyond KdV: models match granular DSWs","Two models describe granular shock waves"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000464,"raw_usage":{"total_tokens":2267,"prompt_tokens":845,"completion_tokens":1422,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":461,"completion_tokens_details":{"reasoning_tokens":1358}},"tokens_in":461,"tokens_out":1422,"duration_ms":11794,"temperature":1.0,"reasoning_tokens":1358,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T18:37:23.960067+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the discrete lattice and both continuum PDEs with exactly zero precompression and a fixed jump, then compare the measured trailing-edge wavenumber and edge speeds; if the PDE values do not approach the lattice values as the numerical smoothing and the artificial $r_+=10^{-5}$ background are removed, the claimed sonic-vacuum validity fails.","supporting_citations":[{"cited_title":"Nesterenko","cited_arxiv_id":null,"evidence_quote":"supplies the granular-chain setting and the sonic-vacuum concept that the new models are designed to cover"},{"cited_title":"Dreyer, M","cited_arxiv_id":null,"evidence_quote":"provides the Whitham modulation equations for the discrete FPUT-type chain used for lattice DSW fitting"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the DSW fitting method from which the edge-speed and wavenumber predictions are derived"},{"cited_title":"Kevrekidis","cited_arxiv_id":null,"evidence_quote":"is the earlier second-order regularized continuum model that the present first-order models extend and compare with"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the BBM regularization idea used to obtain the regularized model (2.14)"},{"cited_title":"Kevrekidis, and Gino Biondini","cited_arxiv_id":null,"evidence_quote":"gives the KdV-based DSW approximation used as the benchmark that the new models must beat"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the KdV dispersive-shock edge formulas that the paper rescales and compares in the precompressed case"}],"review_version":1}