{"id":"e6971257-c580-479d-a60f-b40e118df904","arxiv_id":"2501.09415","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A slowly decaying strain soliton solution is derived for a damped Boussinesq equation with retarded dissipation, and its decay agrees with axisymmetric simulations of the full viscoelastic wave equations.","lead":"This paper models how strain pulses (solitons) fade as they travel along rods made of soft, memory-rich polymers. It derives a wave equation with memory effects and shows that the frequency dependence of the material's nonlinear elasticity must be included to predict how fast the pulses shrink.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (57) gives the initial soliton width as R*|nu0|sqrt(A0*beta00/6) ~ 0.010 mm, but Eq. (38) requires sqrt(6/(A0*beta00)) ~ 5.45 mm; the printed validation input violates the long-wave ordering and is not reproducible as written.","rationale":"The reader's weakest assumption, the small-dissipation ordering against real polymer parameters, is a legitimate caveat and is explicitly disclosed by the authors in the stress tests of Fig. 6. However, the more immediate load-bearing problem is that the manuscript cannot be reproduced as written at the level of the central numerical validation. Substituting the Table 1 values into Eq. (57) gives L0 ~ 0.010 mm, whereas Eq. (38) gives L0 ~ 5.45 mm for the same amplitude. The printed width is two orders of magnitude below the long-wave requirement L >> R* and is smaller than the reported axial grid spacing, so a simulation started with that value would not be testing the reduced Boussinesq model. This is an internal inconsistency in the paper's own setup, not a disagreement with external consensus. If the simulation code used the correct L0, the defect is confined to a typographical error and the central theoretical result likely survives; if not, the Fig. 5 agreement is unexplained. A conditional acceptance is therefore the appropriate verdict: the authors should confirm which L0 was used and correct Eq. (57) before the validation claim is taken as established.","tokens_in":15263,"tokens_out":20002,"duration_ms":196952,"concrete_test":"Re-run the 3D simulation of Sec. 5 exactly as written, with A0 = -2e-3 and L0 = 0.010 mm from Eq. (57), and compare the emerging strain profile to Eq. (38). If no slowly decaying soliton with width ~5 mm forms and no match to Fig. 5 is obtained, the published comparison rests on an unstated different L0. Alternatively, recompute L0 from Eq. (38) and check the stored initial condition in the simulation code; if the stored value is 5.45 mm, the typo is confined to the text but must still be corrected.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The most load-bearing defect is in the validation protocol, not in the asymptotic derivation itself. Section 5 states that the 3D simulation is initialized with the solitary-wave profile of Eq. (38) at amplitude A0 = -2e-3 and width L0 = R*|nu0|sqrt(A0*beta00/6). With the reported values R* = 0.7 mm, nu0 = 0.33, beta00 = -5.5, this printed formula gives L0 ~ 0.010 mm. Eq. (38), however, requires A0 = 6 nu0^2 R*^2 / (beta00 L0^2), so the correct relation is L0 = R*|nu0|sqrt(6/(A0*beta00)) ~ 5.45 mm. The printed expression is the reciprocal of the correct one inside the square root. With L0 = 0.010 mm, the initial condition is two orders of magnitude narrower than the rod gyration radius, violating the long-wave assumption L >> R* that underpins the damped Boussinesq reduction; it is also smaller than the stated axial grid spacing of 0.1 mm, so such a spike could not be resolved by the reported mesh. If the simulations actually used 0.010 mm, Fig. 5 could not exhibit the slowly decaying KdV-type soliton claimed. If they used the correct 5.45 mm, then Eq. (57) misstates the run and the central comparison cannot be independently reproduced from the manuscript as written. Either way, the published validation input is internally inconsistent with the model being validated, and the headline agreement rests on an unstated correction unless the formula is confirmed.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript derives a damped Boussinesq-type equation for longitudinal strain waves in a thin nonlinear viscoelastic rod, starting from a 3D multiple-integral constitutive model with several relaxation processes. Using asymptotic multiple-scales reduction under the assumptions of long waves, small strains, and small dissipation, the authors obtain a damped KdV equation with retarded dissipation and then apply adiabatic soliton perturbation theory to derive ordinary differential equations for the slow decay of a strain soliton's width and velocity. The predicted decay is compared with direct 3D numerical simulations, both with and without nonlinear viscoelastic dissipation. The paper also tests the model outside its formal validity range by increasing the viscoelastic moduli, and it claims that the nonlinear dissipation theory gives better agreement with the 3D simulations than the linear theory.","tokens_in":15536,"tokens_out":10031,"duration_ms":88153,"significance":"Provided the validation input is corrected, the paper offers a significant methodological advance: it gives a systematic route from a nonlinear viscoelastic constitutive law to a reduced 1D model with explicit soliton decay laws, including nontrivial frequency-dependent third-order elastic constants. The inclusion of nonlinear dissipation in the decay ODEs is new and is motivated by experimental measurements on polymers. The stress tests in Fig. 6 are a strong feature, as they probe the limits of the asymptotic model honestly, and the comparison is not fitted: the soliton decay law is derived analytically and the 3D run serves as an independent benchmark. However, the central quantitative claim of agreement with 3D simulations is currently compromised by an inconsistent initial-condition formula in Eq. (57), which makes the reported comparison non-reproducible as written.","major_comments":[{"comment":"The initial-width formula is inconsistent with the soliton relation in Eq. (38). From Eq. (38), A0 = 6 ν0^2 R*^2 / (β00 L0^2), so the correct expression is L0 = R* |ν0| sqrt(6 / (A0 β00)). With the reported values R* = 0.7 mm, ν0 = 0.33, β00 = -5.5, and A0 = -2e-3, this gives L0 ≈ 5.4 mm. Equation (57) instead gives L0 = R* |ν0| sqrt(A0 β00 / 6) ≈ 0.010 mm, which is the reciprocal of the correct ratio inside the square root. A 0.010 mm-wide initial pulse violates the long-wave assumption L >> R* and is smaller than the stated axial grid spacing of 0.1 mm; such an initial condition could not produce the slowly decaying KdV-type soliton shown in Fig. 5. If the simulations actually used the correct 5.4 mm width, then Eq. (57) misreports the initial condition and the comparison is not independently reproducible from the manuscript as written. The authors must correct Eq. (57) and state explicitly which width was used in the 3D runs.","section":"Section 5, Eq. (57)"}],"minor_comments":[{"comment":"There are typos: 'preset study' should be 'present study', and 'ecxept' should be 'except'.","section":"Section 4.1"},{"comment":"The phrase 'the system (39) and (40) can be solved exactly' overstates the status of the solution; the solution relies on the asymptotic approximations I1(θ) ≈ 8θ/15 and I2(θ) ≈ 0. Please say 'solved asymptotically' or 'solved approximately'.","section":"Section 4.1, Eqs. (45)-(48)"},{"comment":"The listed γ_s values (0.33·10^-3) do not exactly equal E_s/E0 = 1.3·10^-3 / 4 = 3.25·10^-4 unless all ν_s equal ν0. Since Eq. (30) allows for frequency-dependent Poisson ratios, the table should state the values of ν_s used, or clarify that the γ_s are the directly prescribed parameters.","section":"Section 5, Table 1"},{"comment":"The notation R*|ν0| sqrt(A0 β00 / 6) is ambiguous because A0 is negative; the product A0 β00 is positive, but the expression should be parenthesized or rewritten with an explicitly positive quantity such as |A0 β00| / 6. The text also refers to 'the amplitude 2·10^-3' while A0 = -2·10^-3; please clarify that this is a compressive amplitude.","section":"Section 5, Eq. (57)"},{"comment":"The label 'nonphysical (unstable)' appears as a fragment in the figure; the caption should explain that the negative-dissipation regime corresponds to the steady-state solution of Eq. (52) and is outside the model's range of validity.","section":"Section 5, Fig. 7"}],"recommendation":"major_revision","confidential_remarks":"The paper is well within the scope of the journal and the asymptotic derivation appears carefully done. The single most important issue is the reproducibility of the validation protocol: Eq. (57) is internally inconsistent with Eq. (38) and the reported agreement with 3D simulations cannot be verified without a correction. If the authors confirm the correct initial width and amend the equation, I would be willing to accept a revised version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this paper deserves a serious referee, but the validation section as printed has an internal inconsistency that must be fixed before the headline comparison can be trusted as reported.\n\nWhat is genuinely new: the authors extend strain-soliton theory from elastic rods to rods of nonlinear viscoelastic material with retarded linear and nonlinear dissipation, covering arbitrary relaxation times instead of only the high- and low-frequency limits treated in Refs. [19-21]. Including frequency-dependent third-order moduli in the dissipation is physically motivated by their earlier polystyrene measurements and is the piece I had not seen before. The reduction from the 3D equations to the damped Boussinesq-type equation (35) is systematic, with the asymptotics clearly stated, and the Karpman-Maslov perturbation derivation of the decay laws (39)-(40) and (52)-(53) is standard and internally consistent. The paper is also honest about its own limits: the small-dissipation ordering epsilon^2 << gamma_s << epsilon is violated when realistic polymer parameters are used, the three-moduli constitutive restriction is stated, and Fig. 7 shows the nonphysical negative dissipation at large amplitudes.\n\nThe load-bearing soft spot is in Section 5, and the stress-test note is right: Eq. (57) as printed contradicts Eq. (38). With the table values, Eq. (57) gives L0 ~ 0.010 mm, while Eq. (38) requires L0 ~ 5.4 mm; the ratio under the square root is inverted. A 0.010 mm spike is two orders of magnitude narrower than the gyration radius, violating the long-wave scaling the whole reduction depends on, and it is smaller than the stated 0.1 mm axial spacing, so the simulation as described could not resolve or evolve it as a soliton. Either the runs actually used ~5.4 mm and the formula is a typo, or the central theory/simulation comparison is not reproducible from the manuscript. This looks fixable, but the authors need to state the actual L0 and correct Eq. (57) before the paper is accepted; the agreement in Fig. 5 cannot be judged otherwise.\n\nMinor points: the '3D' simulation is axisymmetric with one angular point, and the material parameters are synthetic rather than measured. Both reduce the strength of the applied claims without touching the theoretical core.\n\nFor readers working on strain solitons, viscoelastic waveguides, or nonlinear NDE, this is a useful and citable contribution. Send it to review, and make the Eq. (57) correction a condition of acceptance.","headline":"A solid, citable extension of strain-soliton theory to nonlinear viscoelastic rods, but the printed initial-width formula in Eq. (57) inverts the correct relation and makes the central simulation comparison unreproducible as written.","tokens_in":16171,"tokens_out":4999,"would_cite":true,"duration_ms":40340,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["74J30","35Q53","35C08"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that longitudinal strain solitons in nonlinear viscoelastic rods decay according to explicit widening and deceleration laws, and that the reduced one-dimensional model reproduces full 3D simulations of the decay.","keywords":["strain soliton","viscoelastic waveguide","nonlinear viscoelasticity","damped Boussinesq equation","Korteweg-de Vries equation","retarded operators","third-order elastic constants","solitary wave decay"],"falsifier":"The theory predicts a specific critical amplitude at which the nonlinear decay rate in Eq. (52) vanishes, so that larger solitons would stop decaying and become unstable. A decisive test is to simulate or measure the phase portrait $dA/dt$ versus $A$ while increasing the nonlinear viscoelastic moduli $\\beta_{03}$ and $\\beta_{33}$ in small steps; if the zero-crossing does not occur at the predicted amplitude, or if the decay law deviates systematically from the integrated ODEs, the retarded-nonlinear-dissipation model is wrong.","tokens_in":14953,"feed_emoji":"🌊","tokens_out":8963,"duration_ms":87517,"temperature":0.7,"pith_summary":"This paper is trying to establish that long longitudinal strain solitons in a thin rod made of a nonlinear viscoelastic material can be described by a damped Boussinesq-type equation with memory, and that the slow decay of such solitons---widening, amplitude decrease, and deceleration---follows explicit formulas. The authors derive the equation from the full three-dimensional equations of motion under long-wave, small-strain, small-dissipation scaling, then use perturbation theory for nearly integrable systems to obtain the decaying-soliton solution. They verify the predicted decay against direct numerical simulations of the original 3D equations. If the claim is right, a one-dimensional model with retarded linear and nonlinear dissipation is a valid tool for predicting soliton attenuation in viscoelastic waveguides, and frequency-dependent third-order elastic constants cannot be ignored.","feed_headline":"Strain solitons in viscoelastic rods follow a predicted decay law","feed_subtitle":"A one-dimensional model with memory effects matches full 3D simulations of solitary wave attenuation.","key_machinery":"The load-bearing object is the retarded operator $\\hat R_s[f](t)=\\int_{-\\infty}^{t} e^{-(t-t_1)/\\tau_s}\\dot f(t_1)\\,dt_1$ for $s\\ge 1$, with $\\hat R_0$ the identity, which encodes each relaxation process as an exponential memory kernel. Inserted into the stress-strain relation, these operators produce both linear dissipation ($\\gamma_s\\hat R_s[u]$) and bilinear nonlinear dissipation ($\\beta_{su}\\hat R_s[u]\\hat R_u[u]$) in the wave equation. The argument then reduces the Boussinesq equation to a damped Korteweg-de Vries equation by multiple scales and applies adiabatic soliton perturbation theory, which converts the small dissipative terms into quadratures over the soliton shape and yields the ODEs for $L(t)$ and $x_0(t)$.","core_discovery":"The paper's central claim is that the damped Boussinesq-type equation (35) possesses a slowly decaying solitary wave solution of the form $u(x,t)=A(t)\\cosh^{-2}((x-x_0(t))/L(t))$, with $A(t)=6\\nu_0^2R_*^2/(\\beta_{00}L(t)^2)$, where the width $L(t)$ and position $x_0(t)$ obey the ordinary differential equations (39)--(40) for linear dissipation and (52)--(53) when nonlinear dissipative terms are included. The rates are set by universal functions $I_{1,\\mathrm{lin}}$, $I_{2,\\mathrm{lin}}$, $I_{1,\\mathrm{nl}}$, and $I_{2,\\mathrm{nl}}$ of the ratios of relaxation times to the instantaneous soliton width; the linear ones involve the digamma function, and the nonlinear ones involve hypergeometric functions. The paper reports that this solution reproduces the decay observed in full 3D simulations, and that including the frequency dependence of the third-order elastic moduli is essential for the agreement.","pith_inferences":["The paper does not develop this, but the same retarded-operator machinery could be matched to broadband viscoelastic data by increasing the number of relaxation processes $q$, turning soliton decay measurements into a probe of the relaxation spectrum.","The nonzero steady-state amplitude where nonlinear dissipation vanishes, shown in Fig. 7, implies a stability boundary for strain solitons: above a critical amplitude the soliton would grow rather than decay, and measuring that boundary would test the model independently of the small-dissipation ordering.","The paper's stress tests suggest that even when real polymer parameters violate the assumed ordering, the nonlinear theory may still capture the qualitative shape of $dA/dt$ versus $A$; a practical consequence is that soliton-based nondestructive evaluation in such materials should treat nonlinear dissipation as a separately calibrated effect rather than compute it directly from measured moduli."],"forward_implications":["The damped Boussinesq-type equation (35) can serve as a one-dimensional predictive tool for longitudinal strain soliton attenuation in viscoelastic waveguides, replacing full 3D simulations for design and analysis.","Soliton decay is not a simple exponential: with retarded dissipation the width grows as $\\sqrt{1+D_1t}$ in the long-soliton limit and exponentially in the short-soliton limit, with a smooth crossover encoded in $I_{1,\\mathrm{lin}}$.","Nonlinear, frequency-dependent dissipation changes the decay rate substantially; for parameters in the range estimated for polystyrene, the decay is several times slower than the linear-dissipation prediction, so linear models overestimate attenuation.","The same formulas give the soliton speed correction during decay, relating the deceleration to the same memory functions $I_{2,\\mathrm{lin}}$ and $I_{2,\\mathrm{nl}}$, which can be checked by tracking the soliton position in experiments.","For small enough dissipation, the truncated theory matches full 3D simulations, including the relation between decay rate and amplitude, so it extends existing soliton theory for elastic rods to viscoelastic materials with arbitrary relaxation times."],"supporting_citations":[{"why":"Supplies the asymptotic derivation of Boussinesq-type models for longitudinal waves in elastic rods that this paper generalizes to viscoelastic materials.","marker":"[5]"},{"why":"Provides the nonlinear generalization of the standard linear solid constitutive model with retarded operators used in Eq. (4).","marker":"[26]"},{"why":"Supplies the adiabatic perturbation theory for solitons that yields the decay ODEs for width and position.","marker":"[39]"},{"why":"Provides the multiple-scales method used to reduce the wave equation to damped KdV form.","marker":"[37]"},{"why":"Reports the measured frequency dependence of third-order elastic moduli of polystyrene that motivates the nonlinear dissipative terms.","marker":"[25]"},{"why":"Describes the multidomain pseudospectral method used for the full 3D simulations that validate the theory.","marker":"[41]"}],"fun_headline_variants":["Solitons in viscoelastic rods follow a slow decay law","Slowly decaying strain solitons match 3D simulations","Damped Boussinesq equation predicts soliton decay","Memory of elasticity sets soliton decay rate","Frequency-dependent constants drive soliton decay"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole derivation rests on the ordering $\\varepsilon^2\\ll\\gamma_s\\ll\\varepsilon$ together with the long-wave limit $L\\gg R_*$: dissipation must be much weaker than nonlinearity, and the soliton much wider than the rod's cross-section. Real polymers such as polystyrene violate this ordering by one to two orders of magnitude, so for them the predicted decay rates are approximate at best.","fun_headline_variants_meta":{"raw":{"variants":["Solitons in viscoelastic rods follow a slow decay law","Slowly decaying strain solitons match 3D simulations","Damped Boussinesq equation predicts soliton decay","Memory of elasticity sets soliton decay rate","Frequency-dependent constants drive soliton decay"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000521,"raw_usage":{"total_tokens":2495,"prompt_tokens":895,"completion_tokens":1600,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":511,"completion_tokens_details":{"reasoning_tokens":1523}},"tokens_in":511,"tokens_out":1600,"duration_ms":15242,"temperature":1.0,"reasoning_tokens":1523,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:03:51.300066+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"The theory predicts a specific critical amplitude at which the nonlinear decay rate in Eq. (52) vanishes, so that larger solitons would stop decaying and become unstable. A decisive test is to simulate or measure the phase portrait $dA/dt$ versus $A$ while increasing the nonlinear viscoelastic moduli $\\beta_{03}$ and $\\beta_{33}$ in small steps; if the zero-crossing does not occur at the predicted amplitude, or if the decay law deviates systematically from the integrated ODEs, the retarded-nonlinear-dissipation model is wrong.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the asymptotic derivation of Boussinesq-type models for longitudinal waves in elastic rods that this paper generalizes to viscoelastic materials."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the adiabatic perturbation theory for solitons that yields the decay ODEs for width and position."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the multiple-scales method used to reduce the wave equation to damped KdV form."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Describes the multidomain pseudospectral method used for the full 3D simulations that validate the theory."}],"review_version":1}