{"id":"d66e1e5f-8600-4145-81b0-be4fe4b4364f","arxiv_id":"2607.25602","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A Boussinesq-corrected second-order homogenization procedure produces well-posed strain-gradient elasticity models for elastodynamics in arbitrary periodic media, with original reciprocity reductions for the effective inertial tensors.","lead":"Periodic materials like architected lattices can be replaced by simpler effective materials, but the usual improvement step often creates equations that misbehave. This paper adds a corrective weight and proves the resulting strain-gradient wave model is well-posed, while also reducing the amount of microscale computation needed.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Well-posedness proof is internally solid; the load-bearing gap is the unproven O(ϵ^3) asymptotic validity of the modified SGE equation for arbitrary (esp. rough, high-contrast, perforated) periodic media.","rationale":"Good-faith reading: the paper aims to replace ill-posed higher-order homogenized elastodynamic operators with a well-posed, microstructure-dependent SGE model. The well-posedness part is proven: I checked Prop. 4's positivity proofs, the finiteness of ϑe and ϑk (D2 and I2 are positive definite), and the Hille-Yosida application (monotonicity and surjectivity). The reciprocity identities are correct applications of Lemma 2; the reduced cell-problem count is legitimate. Independent support: explicit threshold formulas, no data fitting, and numerical dispersion errors that show the expected O(ϵ^4) accuracy for the three tested cells. The single load-bearing gap is the missing rigorous justification that the modified PDE is the effective equation. This is not an internal inconsistency—it is an unproven formal step. The manuscript openly acknowledges the lack of field-convergence studies and the open boundary-layer problem. My concern matches the reader's weakest_assumption; I therefore recommend no change to the CONDITIONAL verdict, but I emphasize that the burden of proof for the asymptotic claim is on deriving a posteriori estimates or at least a systematic residual computation. Without that, a reader cannot distinguish a genuinely effective well-posed model from a well-posed PDE that merely shares the first terms of a formal expansion.","tokens_in":56756,"tokens_out":25356,"duration_ms":207325,"concrete_test":"Run a transient wave-propagation convergence study on a high-contrast, non-centrosymmetric, perforated 2D cell (e.g., the Z3 cell with a 10× larger modulus contrast and a void). Fix the source and observation time, and compute the relative L2 error between the SG(ϑ) solution and the fully-resolved microstructured solution for ϵ = 1/8, 1/16, 1/32 (with ϑ just above ϑ0 recomputed for each geometry). If the error does not scale as O(ϵ^2) or better, the O(ϵ^3) residual claim is false for such media; if it does, the formal expansion is supported in the regime that matters. An analytic companion: derive the leading-order remainder in (4.29) by substituting the two-scale ansatz (2.16) into the original wave equation and verify it is bounded by C ϵ^3 with C uniform over the cell contrast.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim has two components: (i) the corrected SGE PDE (5.22) is well-posed, and (ii) it is the effective equation of the periodic medium, i.e. U_ϵ satisfies it up to an O(ϵ^3) residual. Component (i) is a theorem that I find internally sound: Proposition 4 follows from the explicit eigenvalue problems (8.6) and (8.8), whose matrices are finite and whose leading tensors D2 and I2 are positive definite, so ϑ0 is finite; the Hille-Yosida argument in §8.5 correctly verifies maximal monotonicity and the energy identity. No circular fitting is involved. The fragile step is component (ii). Equations (3.39)–(3.40) and (4.29)–(4.30) assert the O(ϵ^3) residual on the basis of the formal two-scale expansion (2.7)–(2.8), truncated at O(ϵ^3). The paper provides no a posteriori error estimate for full 2D/3D elastodynamics; the cited justification [21] is leading-order only and [5,3] are scalar/1D. For rough, high-contrast, or perforated cells, the residual may fail to be genuinely O(ϵ^3), in which case the 'well-posed SGE model' is well-posed but not the asymptotic model. The manuscript itself defers field-convergence studies and boundary conditions (§9), and the numerical evidence covers only three 2D cells. Thus the well-posedness claim stands, but the homogenization claim remains conditional.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives second-order two-scale homogenized models for linear elastostatics and elastodynamics in periodic media in arbitrary dimension (d=2,3). It develops reciprocity identities that give alternative expressions for the effective stiffness and inertia tensors, reducing the number of cell problems, and introduces a one-parameter ('Boussinesq trick') family of fourth-order-in-space effective PDEs. The main results are: (i) for a weight ϑ above an explicitly computed threshold, the homogenized strain-gradient elasticity (SGE) energy densities are positive definite; (ii) the corresponding transient wave equation (5.22) is well-posed in free space, proved by verifying maximal monotonicity and applying Hille-Yosida; and (iii) the dispersion relation of the SGE model reproduces the expected second/fourth-order asymptotic behavior in numerical tests on three 2D cells (square, hexagonal, and a non-centrosymmetric Z3 lattice). The well-posedness theorem is internally sound; the fragile part is the claim that the corrected SGE equation is the actual effective equation, which rests on a formal two-scale expansion truncated at O(ε^3) with no a posteriori error estimate supplied for full elasticity.","tokens_in":57168,"tokens_out":4902,"duration_ms":52661,"significance":"If the asymptotic validity were established, this would be a substantial contribution: it provides a systematic, microstructure-based construction of well-posed strain-gradient effective models for arbitrary periodic media, with a transparent, computable threshold for the tunable weight. The reciprocity reductions are practically valuable (e.g., 27 cell problems instead of 90 in 3D dynamics). The well-posedness proof for the corrected model is genuine and does not rely on fitting dispersion data; the effective tensors are computed from cell problems, and the Boussinesq weight is selected from explicit eigenvalue problems. The numerical experiments give honest evidence of the expected asymptotic rates and illustrate the ill-posedness of the uncorrected SG(0) model. However, the claim that the SGE equation is the homogenized model with an O(ε^3) residual is not proven for full 2D/3D elasticity; the manuscript itself defers field-convergence studies to future work. This is the main gap between the paper's title and what is demonstrated.","major_comments":[{"comment":"The central quantitative claim is that U_ϵ satisfies the effective PDEs up to an O(ϵ^3) residual. This is asserted on the basis of the formal two-scale expansion (2.7)-(2.8) truncated at O(ϵ^3). No a posteriori error estimate is given for full 2D/3D elastodynamics; the cited [21] is leading-order only and [5,3] treat scalar or 1D problems. For rough, high-contrast, or perforated cells, the residual may fail to be genuinely O(ϵ^3). The manuscript itself acknowledges in §9 that a quantitative field-convergence study is still needed. This missing support is load-bearing because, without it, Proposition 5 proves well-posedness of an SGE-like PDE but not that this PDE is the effective equation of the microstructure.","section":"§3.3, §4.3, Eqs. (3.39)–(3.40), (4.29)–(4.30)"},{"comment":"The paper claims the entire homogenization procedure carries over to periodically perforated media by replacing Y with Y^s and citing [22]. For the leading-order problem this is standard, but for second-order strain-gradient terms the effect of traction-free holes on the higher-order cell problems and on the O(ϵ^3) residual is not analyzed. Given that all three numerical examples are perforated, the absence of a rigorous justification for this extension weakens the 'arbitrary periodic media' claim exactly in the regime where the central residual estimate is most questionable.","section":"§6.2 (perforated media) and §2.1"},{"comment":"The Hille-Yosida proof is internally consistent and a genuine contribution. However, the statement is restricted to the unbounded-domain problem in R^d with the particular variational definition of R_ϑ^{-1}. The paper correctly notes in §9 that boundary and interface conditions for the SGE model in finite domains remain open. This is a scope limitation rather than an error, but it should be clearly separated from the well-posedness claim in free space, which is fully proved.","section":"§8.5, Proposition 5"},{"comment":"The numerical validation covers only three two-dimensional, relatively smooth periodic cells. Proposition 6 establishes the O(δ^2) expansion of phase velocities under a formal small-δ assumption, but the agreement with Floquet-Bloch results is only empirical and is not used to validate the O(ϵ^3) residual. In particular, the paper does not test the model on rough, high-contrast, or genuinely three-dimensional cells, which are the cases where the missing error estimate is most relevant.","section":"§7 and §6.7, Proposition 6"}],"minor_comments":[{"comment":"The running title contains a typo ('ELASTOST A TICS'); the abstract also renders 'R^d' as 'R d' in several places. Please proofread.","section":"Title/Abstract"},{"comment":"In the proof of Proposition 5, the scalar parameter in (8.15) is written as 'some µ∈R^d'; it should be µ∈R. This is a typographical error and does not affect the argument.","section":"§8.5, surjectivity argument"},{"comment":"The remark refers to 'Prop. 4(b)' for the positivity of the Christoffel eigenvalues, but the relevant property is the positive definiteness of the Fourier symbols Q_ϑ(iξ), R_ϑ(iξ), which is Prop. 4(c). Please correct the cross-reference.","section":"Remark 18"},{"comment":"The discussion of the fourth-order-in-time equation is candid about the lack of well-posedness guarantees and the Ostrogradsky instability. This is useful context, but the paragraph could be shortened, as it is not used in the main line of the paper.","section":"§6.6"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid contribution to higher-order homogenization and fits the journal's scope. The well-posedness theorem (Prop. 5) is sound and interesting, and the reciprocity identities provide a practical computational reduction. The main issue is not an internal inconsistency but an unproven load-bearing claim: the O(ϵ^3) asymptotic validity of the corrected SGE equation for full elasticity. I would encourage the authors to add at least a rigorous statement (or a clearly delimited conjecture with supporting numerical evidence) of the residual estimate, or to reframe the central claim as 'well-posed SGE model formally derived by homogenization plus numerical evidence.' A major revision is appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the good news: Proposition 5 is a real theorem. The Hille-Yosida argument in §8.5 checks out — maximal monotonicity and surjectivity are handled cleanly, the energy identity is correct, and the thresholds ϑ₀ are computable. The reciprocity-based formulas (Props. 1 and 2) are genuinely new for the inertial tensors and give a practical reduction in cell problems. The Boussinesq trick is a useful, honest way to restore sign-definiteness without fitting to dispersion data, and the numerical dispersion curves match the expected 4th-order accuracy. The contrast with the ill-posed SG(0) model is a nice illustration. The authors also know their limitations: they explicitly defer 3D, boundary conditions, and field-convergence studies.\n\nThe soft spot is exactly where the stress-test note puts it: the homogenization claim — that U_ϵ satisfies (5.22) up to O(ε³) — rests on a formal two-scale expansion. No a posteriori error estimate is given for full 2D/3D elastodynamics; the cited works [21] are leading-order, and [5,3] are scalar or 1D. That doesn't sink the paper, but it means the well-posedness theorem is about the corrected SGE operator, not about the actual fine-scale solution. For rough, high-contrast, or perforated cells, the residual might not be genuinely O(ε³). I'd like to see either a rigorous estimate (even in a simplified setting) or a much more careful statement that the effective-model claim is conditional.\n\nMinor issues: no code or data released; only three 2D cells tested; the choice of ϑ affects accuracy (SG(0) is more accurate but ill-posed, SG(ϑ) is less accurate), so practical guidance for picking ϑ beyond the threshold is thin. Also, the paper itself notes that the fourth-order-in-time model has unresolved well-posedness — that's fine, but it's another loose end.\n\nOverall, this deserves a serious referee. I'd accept with major revisions: the core theorem is solid and the reciprocity results are useful. The authors should be pushed to either prove or explicitly bound the residual, or at least to frame the asymptotic claim as conjectural under less regular data. I'd also ask for code and data to back the numerical section. This is a paper I'd cite for the well-posedness construction.","headline":"Solid well-posedness theorem for a Boussinesq-corrected strain-gradient elastodynamic model; the formal asymptotic validity of the effective equation is the main caveat.","tokens_in":57636,"tokens_out":2968,"would_cite":true,"duration_ms":31073,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35B27","74Q05","74J05","74B05"],"pacs":[],"model":"deepseek-v4-flash","headline":"A single tunable weight turns higher-order homogenized wave equations into well-posed strain-gradient models for any periodic elastic medium.","keywords":["homogenization","strain-gradient elasticity","elastodynamics","two-scale asymptotic expansion","Boussinesq trick","well-posedness","reciprocity identities","dispersion relations"],"falsifier":"For a specific non-centrosymmetric high-contrast 2D cell, compute the residual term Oϵ in (3.39)–(3.40) from the computed cell functions and measure how it scales with ϵ; if the L² norm of the residual does not decay as ϵ³, the central asymptotic claim fails. Alternatively, run the well-posed SG(ϑ) model and a Floquet-Bloch reference for a 3D perforated cell and check whether the relative phase-velocity error follows the predicted O(ϵ⁴) trend; deviation would indicate that the formal model, though well-posed, is not the actual asymptotic model.","tokens_in":56651,"feed_emoji":"🌊","tokens_out":11218,"duration_ms":95824,"temperature":0.7,"pith_summary":"This paper establishes that the higher-order homogenized equations obtained by pushing two-scale asymptotic expansions to order ϵ² can be made mathematically sound for any periodic elastic medium in two or three dimensions. The key move is a sign-restoring 'Boussinesq trick': adding a tunable scalar weight ϑ to a Laplacian of the leading-order balance equation replaces the potentially sign-indefinite higher-order effective tensors with ϑ-dependent tensors, yielding a genuine strain-gradient elasticity (SGE) model whose strain and kinetic energy densities are positive definite. With ϑ above an explicitly computable threshold, the resulting transient wave equation has a unique solution for given initial data and forcing, and conserves energy when unforced. Reciprocity identities for the cell problems reduce the required computations to two static cell functions and one inertial cell function. The upshot is a practical, well-posed effective model that captures anisotropic and dispersive wave propagation in architectured materials at a fraction of the cost of resolving the microstructure, with dispersion errors of fourth order in the scale parameter.","feed_headline":"One tunable weight tames ill-posed homogenized wave equations","feed_subtitle":"Makes strain-gradient elastic wave models mathematically sound and computable from a few cell solutions.","key_machinery":"The central mechanism is the Boussinesq trick: adding to the O(ϵ²) macroscopic balance equation a tunable scalar weight ϑ times the Laplacian of the leading-order balance equation (and, in dynamics, an additional identity obtained by applying a filtered differential operator to the leading-order equation). This converts the raw effective tensors C2 and B2 — which can be sign-indefinite and are responsible for ill-posedness — into ϑ-dependent tensors A2(ϑ)=ϑD2−C2+c2 and J2(ϑ)=ϑI2−B2+b2 with the symmetry and positivity needed for a valid SGE model. The threshold ϑ0 is obtained from the symmetric eigenvalue problems (8.2)/(8.6)/(8.8), and the reciprocity identities of Propositions 1 and 2 suppl","core_discovery":"Central claim: the fourth-order PDE from second-order homogenization of elastodynamics, generically ill-posed due to wrong-sign higher-order tensors, can be recast as a well-posed strain-gradient (SGE) model. The Boussinesq trick adds a tunable weight ϑ times the Laplacian of the leading-order balance equation, producing tensors A2(ϑ)=ϑD2−C2+c2 and J2(ϑ)=ϑI2−B2+b2. For ϑ above a threshold from small eigenvalue problems, strain and kinetic energies are positive definite and Hille-Yosida gives a unique solution with conserved energy; statics follow the same route. Reciprocity identities compute all effective tensors from χ1, χ2, ζ2, skipping the third-order cell problems.","pith_inferences":["A testable extension is to choose ϑ just above the threshold, or to use a direction- or frequency-dependent weight, to recover some of the dispersion accuracy that the numerical examples show is lost when ϑ is taken larger than necessary.","Because the paper proves well-posedness but not a rigorous a posteriori error estimate for full elasticity, the practical validity for high-contrast, perforated, or 3D cells remains open; a direct residual computation for such cells would settle whether the O(ϵ³) claim holds.","The paper's outlook identifies boundary and interface conditions preserving O(ϵ³) accuracy as an open problem in higher dimensions; without such conditions, finite-domain simulations with the SGE model will require ad-hoc closures, limiting quantitative use.","The reciprocity-based reduction suggests a cheap computational pipeline — solve χ1, χ2, ζ2, then two eigenvalue problems — that could be embedded in topological optimization of microstructures for target dispersion, matching the paper's stated de-homogenization motivation."],"forward_implications":["For any periodic medium satisfying the minimal ellipticity and boundedness assumptions, the transient SGE effective wave equation has a unique solution for ϑ above the threshold, with energy conservation in the unforced case.","The effective dispersion relations come from a generalized Christoffel equation whose matrices are Hermitian positive definite for ϑ>ϑ0, and the phase-velocity expansion has a real O(δ²) correction; the O(δ) term vanishes except for non-centrosymmetric cells with a double leading-order eigenvalue.","Only the first and second static cell functions (χ1, χ2) and the second inertial cell function (ζ2) are needed to build the model; the third-order cell problems are not required for evaluating the effective tensors, cutting the 3D cell-solution count from 90 to 27.","The dynamic strain-gradient elasticity tensor differs from its static counterpart whenever the mass density is heterogeneous, so using a static SG tensor in dynamics requires the fourth-order-in-time term or the modified tensor A2(ϑ).","The well-posedness statement applies beyond this homogenization context to any SGE material whose operators satisfy the positivity and coercivity properties of Proposition 4."],"fun_headline_variants":["Tunable weight restores well-posedness in homogenized elastodynamics","Boussinesq trick yields well-posed strain-gradient wave models","One scalar weight fixes ill-posed homogenized wave PDEs","Well-posed strain-gradient elasticity from a tunable weight","Arbitrary periodic media now admit well-posed gradient-wave models"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the truncated two-scale expansion, together with the Boussinesq-trick modifications, yields an effective PDE with a true O(ϵ³) residual for full elasticity; the paper supplies no rigorous a posteriori error estimate for general 2D/3D elastic media, so if the residual is worse than O(ϵ³) for rough, high-contrast, or perforated cells, the well-posed model may be approximating the wrong asymptotic limit.","fun_headline_variants_meta":{"raw":{"variants":["Tunable weight restores well-posedness in homogenized elastodynamics","Boussinesq trick yields well-posed strain-gradient wave models","One scalar weight fixes ill-posed homogenized wave PDEs","Well-posed strain-gradient elasticity from a tunable weight","Arbitrary periodic media now admit well-posed gradient-wave models"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000183,"raw_usage":{"total_tokens":1206,"prompt_tokens":855,"completion_tokens":351,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":599,"completion_tokens_details":{"reasoning_tokens":275}},"tokens_in":599,"tokens_out":351,"duration_ms":4159,"temperature":1.0,"reasoning_tokens":275,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T01:57:45.310824+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For a specific non-centrosymmetric high-contrast 2D cell, compute the residual term Oϵ in (3.39)–(3.40) from the computed cell functions and measure how it scales with ϵ; if the L² norm of the residual does not decay as ϵ³, the central asymptotic claim fails. Alternatively, run the well-posed SG(ϑ) model and a Floquet-Bloch reference for a 3D perforated cell and check whether the relative phase-velocity error follows the predicted O(ϵ⁴) trend; deviation would indicate that the formal model, though well-posed, is not the actual asymptotic model.","supporting_citations":[],"review_version":1}