{"id":"b5dc5cd6-7d6b-4868-b427-157873949d1d","arxiv_id":"2508.06324","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"An effective nonlinear dispersive wave equation for shear waves in magneto-tunable soft laminates is derived by homogenization, showing that solitary wave speed limits can be adjusted with magnetic fields and microstructure.","lead":"This paper derives an approximate equation describing how shear waves move through a soft layered material that can be tuned with a magnetic field, including nonlinear and dispersive effects that create solitary waves. The result gives designers a recipe for controlling wave speeds and band gaps in magneto-active composites.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Solitary-wave predictions violate the homogenization's own scale-separation assumption: at s=1.026c, Eq. (72) gives L/ℓ≈0.3–0.44 and δ≈1.8, so ε=ℓ/L is O(2–3) and δ is O(1). The speed bound (83) is reached exactly as L/ℓ→0, opposite to the long-wave limit used to derive Eq. (2).","rationale":"The reader's weakest_assumption focused on the truncated layer-wise model and fixed remnant magnetization, and noted the lack of direct validation of the full homogenized equation. I agree that the absence of a layer-resolved test of Eq. (64) is an important gap. But the more load-bearing and more specific problem is that the paper's own solitary-wave examples are computed in a regime where the homogenization expansion has already broken down: at s=1.026c, the predicted soliton width is L/ℓ≈0.3–0.44 and the strain amplitude is δ≈1.8, whereas the derivation assumes ε=ℓ/L≪1 and δ≪1. The claimed magnetically tunable speed bound is reached as L/ℓ→0, the exact opposite of the long-wave limit. This is an internal-consistency concern: the asymptotic equation is used outside its formal validity, and the paper provides no direct numerical evidence that the layer-resolved medium actually supports these solitary waves. I do not see a reason to reject the derivation itself; the formulas for ζ and η are internally coherent, and the mKdV validation at wavelengths 8ℓ–16ℓ is a useful check. However, acceptance of the central claim should be conditional on either demonstrating, via layer-resolved simulations, that the predicted solitary waves and the speed bound actually occur in the laminate for the quoted parameters, or on substantially restricting the claims to the regime where ε≪1 and δ≪1 hold. For those reasons I maintain the conditional verdict, with the condition sharpened to require scale-consistent validation of the solitary-wave predictions.","tokens_in":33942,"tokens_out":28720,"duration_ms":291674,"concrete_test":"Use the finite-volume layer-resolved code of Section 5.4 (with Table 3 parameters) to initialize the exact travelling-wave strain profile (72) from the full homogenized model at s=1.026c (L/ℓ≈0.30, δ≈1.85) and, for comparison, at s=1.001c (L/ℓ≈2.2, δ≈0.36). Propagate for at least 20 periods. If the small-amplitude long-wavelength profile (s=1.001c) propagates with constant shape and speed matching (72), the asymptotic model is trustworthy in its valid regime; if the s=1.026c profile also persists with negligible radiation and no change of speed, the solitary-wave/bound claims may survive outside the formal regime. If the narrow large-amplitude profile distorts, sheds radiation, or fails to persist, then the predicted speed-amplitude relation and the upper bound (83) are artifacts of extrapolating Eq. (64) beyond its derivation.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The homogenized equation (51) is derived under two explicit smallness assumptions: ε=ℓ/L≪1 (long waves) and δ=a/L≪1 (moderate amplitude), see Section 4.1. The subsequent travelling-wave analysis of Section 5.3 defines the soliton width via Eq. (72) as L=ℓ/√c1, with c1 given by Eq. (69) for the full model (64) or Eq. (70) for the mKdV model (65). For the paper's own parameters (Table 3) and the illustrative speed s=1.026c, the full model gives c1≈10.96, hence L/ℓ≈0.30 and δ≈1.85; the mKdV model gives L/ℓ≈0.44 and δ≈1.77. These numbers are not small: the soliton is narrower than the microstructure period and its strain amplitude is order one. The magnetically tunable upper bound reported in Eq. (83) is realized in the limit L/ℓ→0, which is the extreme violation of the long-wave assumption. Thus the solitary-wave branch (72) and the speed bound (83) are asymptotic predictions applied outside their derivation's domain of validity. Direct numerical evidence for these solitary waves in the layered medium is not provided; the only layer-resolved validation in Section 5.4 concerns the mKdV approximation of the impact problem, at wavelengths 8ℓ–16ℓ, and in the nonlinear dispersive case even that comparison is only qualitative. This is an internal consistency issue, not merely a missing experiment: the effective equation is used in a regime where the two-scale expansion that produced it is not justified.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies nonlinear shear waves in periodic hard-magnetic soft laminates. It first revisits the hard-magnetic constitutive framework to use a symmetric total stress and, for generalized neo-Hookean layers under a static magneto-deformation, reduces the layer dynamics to the nonlinear equation (1) with cubic stiffness. Using two-scale homogenization along the lines of Andrianov et al., it derives the effective Boussinesq-type equation (2) with closed-form effective coefficients (Table 1), and an approximate effective strain-energy function (56). It then improves the linear dispersive representation by means of three parameters, reduces the equation to an mKdV model, studies sech solitary waves, derives speed-amplitude relations and an upper bound on the solitary-wave speed, and demonstrates magnetic tunability of band gaps and of the solitary-wave speed bound. The mKdV model is compared with finite-volume simulations of an impact problem in a layered medium.","tokens_in":34455,"tokens_out":14439,"duration_ms":151053,"significance":"If taken at face value, the paper provides a compact, closed-form description of nonlinear dispersive shear waves in magneto-active laminates, with no fitting of the central homogenized coefficients, and it gives concrete, falsifiable predictions for magnetically tunable solitary-wave speeds. The comparison with the exact Floquet-Bloch dispersion relation is a strong point, as is the explicit derivation of an effective strain energy. The main reservation is that the solitary-wave analysis is presented in a regime where the homogenization's scale-separation parameters are not small; the predictions are interesting but currently not supported by layer-resolved simulations in that regime. With appropriate restrictions or additional validation, the paper would be a useful contribution.","major_comments":[{"comment":"The solitary-wave branch violates the scale-separation hypotheses ε=ℓ/L≪1 and δ=a/L≪1 under which Eq. (2) was derived. For Table 3 with s=1.026c, Eq. (69) gives c1≈10.9, so L/ℓ≈0.30 and δ≈1.85; the mKdV branch (70) gives L/ℓ≈0.44 and δ≈1.77. Thus ε and δ are O(1), not small. The speed bound (83) is approached exactly as L/ℓ→0, since at s²/c²=√(1+η/ηt) the coefficient c3 in Eq. (69) diverges and L→0, which is the extreme violation of the long-wave limit. The paper's own Conclusion states that Eq. (64) is restricted to the low-frequency range and to waves of moderate amplitude. The layer-resolved validation in §5.4 covers only the mKdV model at wavelengths 8ℓ–16ℓ, and in the nonlinear dispersive case the comparison is qualitative. Please either restrict the soliton claims to the valid asymptotic regime, provide direct layer-resolved simulations of the solitary waves in that regime, or clea","section":"§4.1, Eq. (47); §5.3, Eqs. (69)-(72), (83)"},{"comment":"The numerical validation in §5.4 does not exercise the full homogenized model used for the solitary-wave branch. The mKdV reduction (65) is independent of the splitting (ηy,ηm,ηt) as long as Eq. (58)2 holds, so the agreement in Figs. 5–6 says little about the modified dispersion model (58)/(62) or about the solitary waves (69)-(72) that underlie the speed bound (83). In the nonlinear dispersive case (iii), the authors themselves state that the comparison is only qualitative. Moreover, the coefficients (62) are chosen to match the first band gap and are not derived from the two-scale expansion. To support the central speed-amplitude claim, the full model (64) should be compared with layer-resolved simulations, or the claim should be presented as a prediction of the homogenized model only.","section":"§5.2-§5.4"}],"minor_comments":[{"comment":"In the change of variables {ŷ=ε²y, t̂=t−y/c, u=εû}, the parameter ε is said to be 'of the same order as the microstructure's characteristic length ℓ'. Since y and u are dimensional, ε must be dimensionless. Please clarify by writing ε=ℓ/L or by introducing normalized variables, so that the ordering is unambiguous.","section":"§5.2, scaling paragraph"},{"comment":"The nested square roots in these equations are typeset ambiguously. Add explicit parentheses, e.g. maxδ = sqrt( (sqrt(1+η/ηt)−1)/(ζ/6) ), so the reader can verify the numerical values.","section":"Eqs. (74), (84)"},{"comment":"It would help to state explicitly that h(α) has the same units as g(α), so that ζ is dimensionless, and that the expression for η is dimensionless only after using c²=⟨g⟩/⟨ρ⟩.","section":"Table 1"},{"comment":"The artificial mass density choices are clearly flagged as a thought experiment. A short remark on whether the shock-formation distance (76) changes when c is modified by those density choices would avoid possible confusion.","section":"§5.4, footnotes 2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the journal's scope and the derivation is generally careful. The main obstacle is the scale-separation issue in the solitary-wave analysis: the central speed bound is presented in a regime where the homogenization hypotheses are violated, and the existing layer-resolved simulations do not validate the full nonlinear dispersive model. If the authors can provide layer-resolved simulations of the solitary-wave branch in a valid regime, or explicitly reframe the results as properties of the truncated model outside its homogenization domain, I would support publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper delivers a genuinely new effective model: a cubic nonlinear, fourth-order dispersive equation for shear waves in hard-magnetic soft laminates, with coefficients in closed form (Table 1). I checked the homogenization logic and it follows Andrianov et al. cleanly; the constitutive revision (symmetric stress) is a real correction to earlier work. The linear dispersion comparison with the exact Floquet–Bloch relation is convincing, and the authors are honest that their improved dispersion model is a calibration, not a derivation. The impact simulations showing mKdV's limitations are useful and appropriately skeptical.\n\nBut the solitary-wave section has a load-bearing consistency problem. The homogenized equation was derived under ε=ℓ/L≪1 and δ=a/L≪1. The traveling-wave analysis then produces solitons whose width is L/ℓ≈0.3–0.44 and strain amplitude δ≈1.8 for the paper's own parameters (s=1.026c). Those are not small; the soliton is narrower than the microstructure period. The magnetic tuning bound in Eq (83) is reached in the limit L/ℓ→0, which is the extreme violation of the long-wave assumption. So the paper's headline claim—an adjustable upper bound on solitary wave speed—is a prediction of the PDE, not a justified prediction about the laminate. No direct numerical simulation of these solitons in the layered medium is provided. The only layer-resolved validation is the mKdV impact problem at wavelengths 8ℓ–16ℓ, and in the nonlinear dispersive case even that is qualitative.\n\nThis is not a fatal flaw in the homogenization itself, but it is more than a minor caveat. The authors need to either (a) simulate solitary waves in the layered equations to demonstrate that they exist in the regime they advertise, or (b) explicitly present the solitary-wave analysis as a formal property of the effective equation, with a clear discussion of why it may extrapolate beyond the derivation. They already state the general restrictions in the conclusion; they should apply the same care to the solitary-wave claims.\n\nBottom line: the effective equation is worth having, and the paper deserves peer review. A serious referee should push on the solitary-wave regime and the unvalidated full equation (64), but the core derivation is sound.","headline":"New effective equation is a solid homogenization result; the solitary-wave section applies it where the two-scale expansion does not hold, so the speed-bound claims need reframing or direct simulation.","tokens_in":34887,"tokens_out":2583,"would_cite":true,"duration_ms":28686,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["74Q05","74J30","74F15","35Q51"],"pacs":[],"model":"deepseek-v4-flash","headline":"Nonlinear shear waves across a soft magneto-active laminate are governed by a single effective dispersive wave equation, with explicit coefficients and magnetically tunable solitary-wave solutions.","keywords":["homogenization","nonlinear shear waves","magneto-active laminates","hard-magnetic soft solids","solitary waves","mKdV equation","band gaps","dispersion"],"falsifier":"Measure the maximum speed of shear solitary waves in a hard-magnetic Gent laminate under a range of magnetic inductions and compare with the predicted bound; exceeding the bound, or finding the bound insensitive to the field, would falsify the homogenised model. Alternatively, compute the next-order amplitude correction in the homogenisation expansion and check whether its coefficient is negligible compared with the nonlinearity coefficient at moderate amplitudes.","tokens_in":33872,"feed_emoji":"🧲","tokens_out":6004,"duration_ms":67280,"temperature":0.7,"pith_summary":"The paper derives a single effective wave equation that governs finite-amplitude shear waves travelling across a periodic stack of soft magneto-active layers. Starting from a layer-wise nonlinear wave equation with cubic stiffness, asymptotic homogenisation in the long-wave, moderate-amplitude regime yields a dispersive nonlinear equation whose coefficients are explicit functions of the two phases' volume fractions, stiffnesses, and densities. The equation supports solitary waves with a computed speed–amplitude relation, and the paper shows that the maximum admissible solitary-wave speed can be moved by changing an applied magnetic field or the microstructure. If correct, this gives a design tool for tunable acoustic filters and wave control in soft magneto-active composites.","feed_headline":"Magnetic field tunes the speed limit of shear solitons in soft laminates","feed_subtitle":"One effective dispersive wave equation with explicit coefficients governs nonlinear shear waves in layered magneto-active materials.","key_machinery":"The load-bearing object is the homogenized nonlinear dispersive wave equation $$$c^{2}$(1+\\zeta $u_y^{2}$)u_{yy}+\\eta\\$ell^{2}$ $c^{2}$ u_{yyyy}=u_{tt},$$ with effective coefficients $\\zeta$ and $\\eta$ given by explicit formulas in terms of the layer properties. This equation carries the argument by replacing the alternating-layer microstructure with a homogeneous dispersive medium, so that nonlinearity, dispersion, and magneto-elastic tunability are all encoded in a handful of coefficient formulas.","core_discovery":"The paper's central claim is that a periodically layered soft magneto-active solid can be replaced, in the long-wave moderate-amplitude regime, by a homogeneous nonlinear dispersive medium governed by $$$c^{2}$(1+\\zeta $u_y^{2}$)u_{yy}+\\eta\\$ell^{2}$ $c^{2}$ u_{yyyy}=u_{tt},$$ where $c^2=\\langle g\\rangle/\\langle\\rho\\rangle$ and the dimensionless coefficients $\\zeta$ and $\\eta$ are explicit functions of the two phases' volume fractions, shear stiffnesses, and densities. The paper further claims that this equation admits solitary waves with an explicit sech profile, a speed–amplitude relation, and a magnetically tunable upper bound on speed, and that a unidirectional modified Korteweg–de Vries reduction captu","pith_inferences":["A natural extension is to relax the fixed-remnant-magnetisation assumption; allowing magnetisation fluctuations would likely add magnetic coupling terms to the effective equation and could shift the solitary-wave speed bound.","Because the dispersion coefficient is an explicit function of the impedance contrast between phases, measuring the first band-gap edge could serve as a non-destructive route to infer layer properties.","Applying the same homogenisation scheme to the two-polarisation shear system would likely yield a vector dispersive equation supporting polarised solitary-wave families.","The optimised homogenised model, being much cheaper than finite-volume simulations of the layered medium, could act as a surrogate in design optimisation of magnetically tunable acoustic filters."],"forward_implications":["The first shear band gap of a soft magneto-active laminate is approximated by the homogenised dispersion relation with optimised coefficients, so band-gap edges can be estimated directly from phase properties.","In Gent-type laminates the effective dispersion coefficient and the solitary-wave speed bound change with the applied magnetic field through the static stretch, giving a post-fabrication tuning knob.","Solitary waves in the homogenised model exist only for relative speeds satisfying a polynomial inequality; approaching the bound shrinks the wavelength to zero while saturating the strain amplitude.","The mKdV reduction describes unidirectional propagation accurately only for speeds close to the linear wave speed, and should not be used for quantitative high-amplitude predictions.","The derived effective strain energy links dynamic homogenisation to a static homogenised material theory, so static and dynamic laminate responses can be described within one model in the small-nonlinearity limit."],"supporting_citations":[{"why":"Supplies the two-scale asymptotic homogenisation procedure used to derive the effective nonlinear dispersive wave equation.","marker":"[7]"},{"why":"Provides the hard-magnetic laminate model and the linear layer-wise wave equation that this work extends to nonlinear amplitudes.","marker":"[75]"},{"why":"Gives the revised hard-magnetic energy formulation with symmetric total Cauchy stress used for the constitutive model.","marker":"[29]"},{"why":"Supplies higher-order cell solutions and the exact Floquet–Bloch dispersion relation used for validation of the homogenised model.","marker":"[6]"},{"why":"Provides the optimised dispersion coefficients and the band-gap matching strategy used to improve the homogenised dispersion relation.","marker":"[69]"},{"why":"Supplies the finite-volume layered-medium simulations and the solitary-wave phenomenology used for the impact-problem comparison.","marker":"[73]"},{"why":"Gives the neo-Hookean effective strain energy for laminates that the newly derived effective energy generalises.","marker":"[35]"},{"why":"Provides the remnant-magnetisation constitutive law for hard-magnetic soft materials that the paper revises and builds upon.","marker":"[76]"}],"fun_headline_variants":["Magnetic field tunes soliton speed limit in soft laminates","Layered soft magnets cap shear soliton speeds","Magnet-tunable shear solitons in soft laminates","Speed ceiling for shear solitons set by magnetic field","Magneto-elastic laminates: magnetic tuning of soliton speed"],"cache_read_input_tokens":2816,"weakest_assumption_plain":"The load-bearing premise is that each layer's shear stiffness is exactly the quadratic truncation $g^{(\\alpha)}+\\tfrac{1}{3}h^{(\\alpha)}u_y^2$ and that the remnant magnetisation stays constant during wave propagation; if the neglected magnetisation fluctuations or fourth-order strain terms matter at the amplitudes considered, the effective equation, the mKdV reduction, and the solitary-wave speed bounds all change.","fun_headline_variants_meta":{"raw":{"variants":["Magnetic field tunes soliton speed limit in soft laminates","Layered soft magnets cap shear soliton speeds","Magnet-tunable shear solitons in soft laminates","Speed ceiling for shear solitons set by magnetic field","Magneto-elastic laminates: magnetic tuning of soliton speed"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000215,"raw_usage":{"total_tokens":1309,"prompt_tokens":828,"completion_tokens":481,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":572,"completion_tokens_details":{"reasoning_tokens":397}},"tokens_in":572,"tokens_out":481,"duration_ms":5167,"temperature":1.0,"reasoning_tokens":397,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T22:47:43.823002+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the maximum speed of shear solitary waves in a hard-magnetic Gent laminate under a range of magnetic inductions and compare with the predicted bound; exceeding the bound, or finding the bound insensitive to the field, would falsify the homogenised model. Alternatively, compute the next-order amplitude correction in the homogenisation expansion and check whether its coefficient is negligible compared with the nonlinearity coefficient at moderate amplitudes.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the two-scale asymptotic homogenisation procedure used to derive the effective nonlinear dispersive wave equation."},{"cited_title":"Zhang and S","cited_arxiv_id":null,"evidence_quote":"Provides the hard-magnetic laminate model and the linear layer-wise wave equation that this work extends to nonlinear amplitudes."},{"cited_title":"Dorfmann and R","cited_arxiv_id":null,"evidence_quote":"Gives the revised hard-magnetic energy formulation with symmetric total Cauchy stress used for the constitutive model."},{"cited_title":"Wautier and B","cited_arxiv_id":null,"evidence_quote":"Provides the optimised dispersion coefficients and the band-gap matching strategy used to improve the homogenised dispersion relation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the finite-volume layered-medium simulations and the solitary-wave phenomenology used for the impact-problem comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the neo-Hookean effective strain energy for laminates that the newly derived effective energy generalises."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the remnant-magnetisation constitutive law for hard-magnetic soft materials that the paper revises and builds upon."}],"review_version":1}