{"id":"79f27090-ae4a-458d-9a5e-5beb95b8df0a","arxiv_id":"1908.01640","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Closed-form constitutive tensors are obtained for a second-gradient elastic continuum equivalent to a three-stiffness hexagonal truss, and validated against lattice solutions for shear and uniaxial strain.","lead":"This paper derives a continuous 'second-gradient elastic' description for a hexagonal lattice of hinged springs, including internal length and direction-dependent effects. It compares the continuum model with the actual lattice under shear and compression, showing the model captures the lattice response better than the classical elastic model.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The standard SGE extension is underdetermined: Eq. (51) fixes four free continuation coefficients by an unproved relaxation rule, and the four aligned validation cases (errors up to 19.3%) do not establish general energetic equivalence; an off-axis BVP would test the continuation.","rationale":"The reader's weakest assumption is precisely the point where I find the argument least secure. Equations (14)-(16) leave an underdetermined passage from the 4D condensed SGE to the 6D standard SGE; Eq. (51) collapses it by a rule stated but not derived. The four validation BVPs are genuine out-of-sample checks and give real support, including the non-centrosymmetric transverse displacement of Fig. 11 that the Cauchy model cannot produce. However, the errors in Tables 2 and 3 are not uniformly small, and all loadings are aligned with x1/x2, so the off-axis anisotropic and null-space components of the continuation are not probed. The undefined a45 in Eq. (71) and the absence of the coupled closed-form solutions and Mathematica code weaken reproducibility but do not by themselves disprove the central claim. I therefore keep the reader's CONDITIONAL verdict: the paper should be accepted only after the continuation rule is justified (or validated on an off-axis case) and the missing artifacts are supplied.","tokens_in":23996,"tokens_out":13359,"duration_ms":148167,"concrete_test":"Compute the energy mismatch (U_lat - U_SGE)/U_lat for a uniaxial strain imposed at 30 degrees to the x1 axis on the same strip geometry used in Section 3.2, solving the lattice exactly and the SGE with the C, M, A matrices from Eq. (1). If the SGE mismatch is not substantially below the Cauchy mismatch, as it is in the aligned cases of Tables 2 and 3, then the relaxation rule in Eq. (51) does not define a generally equivalent SGE material.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that C, M, and A in Eq. (1) define a second-gradient continuum energetically equivalent to the hexagonal lattice for equilibrium deformation fields beyond the quadratic identification fields. Part I fixes only the restriction of the energy to the 4D 'condensed' curvature subspace (Eqs. (6)-(7)); the extension to the full 6D curvature space is underdetermined. After imposing Eq. (15) and the symmetry classes, four coefficients remain free: Q35, DeltaA11, DeltaM15, DeltaM16. These are set in Eq. (51) by a 'relaxation' rule described in words rather than derived, so it is not shown that the resulting M and A reproduce the lattice energy on the two non-condensed curvature modes. The Section 3 validation uses four BVPs, all aligned with the x1/x2 lattice axes, and the improvement over Cauchy is real but not uniform: Ex2 simple shear parallel to x2 retains 19.3% energy error (Table 2), and three uniaxial cases have negative errors, meaning the SGE is stiffer than the lattice (Table 3). The validation is evidence but not proof that the chosen continuation is the physically correct SGE for general loads; if it is not, the closed-form model is not the claimed equivalent material.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper completes the identification of a second-gradient elastic (SGE) continuum equivalent to a planar hexagonal lattice of axially deformable bars. Starting from the 'condensed' constitutive matrices C, M*, A* derived in Part I under equilibrium-satisfying quadratic displacement fields, the authors construct a 'standard' SGE material with constitutive matrices C, M, A given in closed form in Eq. (1). They analyze the symmetry classes of these tensors (isotropic C, Z3-symmetric M, D6-symmetric A), study positive definiteness of the condensed energy and its dependence on the bar stiffness ratios, and validate the model by comparing the lattice response with the SGE and Cauchy continuum responses under simple shear and uniaxial strain in four aligned loading configurations. The validation is reported as a reduction of strain-energy mismatch when the SGE model replaces the classical Cauchy model.","tokens_in":24233,"tokens_out":4158,"duration_ms":40950,"significance":"If the identification is accepted, the paper provides a closed-form, parameter-free second-gradient continuum model of a hexagonal lattice, with an internal length and non-centrosymmetric coupling terms, which is of genuine interest for micromechanics and metamaterial design. The authors give exact analytical lattice solutions (obtained with Mathematica) and explicit formulas for all constitutive coefficients, so the paper is checkable and reproducible in principle. The validation is not circular: the constitutive parameters are obtained from energy matching, not from fitting the validation responses, and the paper demonstrates a clear non-centrosymmetric effect (transverse displacement under uniaxial strain) that a Cauchy model cannot capture. The main weaknesses are the underdetermined extension from the condensed to the standard form and the mixed quality of the validation results, which fall short of the claimed 'excellent agreement'.","major_comments":[{"comment":"The extension from the condensed to the standard SGE form is underdetermined. Equations (14)-(16) leave four coefficients (Q35, ΔA11, ΔM15, ΔM16) unspecified by energy matching, and Eq. (51) fixes them by a rule described in words ('imposing an equivalence ... in which the equilibrium constraint on the coefficients βijk is neglected') rather than derived from an explicit energy-matching or minimization criterion. Since the matrices M and A in Eq. (1) depend on this choice, the claim that the resulting SGE material is energetically equivalent to the lattice for arbitrary displacement fields is not established. A concrete test would be an off-axis boundary value problem (e.g., simple shear or bending with a loading direction not aligned with the x1/x2 lattice axes) that activates the non-condensed curvature modes; the current validation only uses four aligned cases (Tables 2 and 3).","section":"§2.4, Eq. (51)"},{"comment":"The reported errors undercut the 'excellent agreement' language in the Abstract and Section 1. For Ex2, simple shear parallel to the x2-axis, the SGE energy error is 19.3%, compared with 26.74% for Cauchy, so the improvement is modest and the absolute error remains large. For the same loading, Ex3 shows an SGE error of only 1.14%, indicating strong configuration dependence. The paper should either soften the global validation claim or explain why the x2-aligned shear case is an outlier.","section":"§3.1.3, Table 2"},{"comment":"Negative discrepancies (e.g., -6.83%, -7.47%, -8.23% for the uniaxial parallel to x1 cases) show that the SGE model is stiffer than the lattice in these configurations. This contradicts the impression of a consistently improved match and means the energy comparison is not one-sided. The statement that the advantage 'becomes clear' for uniaxial strain parallel to x2 is not supported by the last column of Table 3, where the SGE errors are -1.24%, 0.57%, and -0.72%. The authors should address why the equivalent SGE solid can be stiffer than the lattice and what this implies for the claimed energetic equivalence.","section":"§3.2.3, Table 3"},{"comment":"The expression for τ211(x2) contains the coefficient a45, which is not defined in Eq. (1) or elsewhere in the paper; the A matrix in Eq. (1) uses only a11, a12, a16, a26, and a34. As written, the closed-form solution for uniaxial strain parallel to the x2-axis is not reproducible. Please define a45 or correct the typo (possibly it should be a26 or a16).","section":"§3.2.1, Eq. (71)"},{"comment":"The assertion that the chosen values of Q35, ΔA11, ΔM15, and ΔM16 provide a positive-definiteness domain identical for the condensed and standard forms is not proved and is not obvious, because the standard form operates on a larger curvature space (six curvature components versus four in the condensed form). A symbolic or numerical verification of this equality should be provided, or the statement should be qualified.","section":"§2.4, after Eq. (51)"}],"minor_comments":[{"comment":"The word 'verifed' should be 'verified'.","section":"§2.1, Eq. (33)"},{"comment":"The word 'isoropic' in the sentence about the C matrix should be 'isotropic'.","section":"§2.2"},{"comment":"The phrase 'first indentification step' contains a typo: 'indentification' should be 'identification'.","section":"§2.4"},{"comment":"The header lists the ratios as 'ˆk/k and ˆk/k'; the second ratio should be '~k/k'.","section":"Table 1 header"},{"comment":"The statement that the dimensionless energy difference coincides with the dimensionless difference of the resultant tractions is non-obvious and should be derived or referenced, since it is used to interpret Tables 2 and 3.","section":"§3.1.3, footnote 2"},{"comment":"In the last row of R[q](θ), there is a misplaced period after '( c(θ) +c(3θ))/2.'; the notation should be cleaned up.","section":"Eq. (20)"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the authors' reliance on their own Part I is legitimate, as Part I is published separately. The main risk is the underdetermined relaxation rule in Eq. (51); the validation tests are all aligned with the lattice axes and therefore do not exercise the non-condensed curvature modes. If the authors can provide an off-axis boundary value problem or another test that discriminates among the possible continuations, the central claim would be much stronger. The issues with Tables 2 and 3 and the undefined a45 are fixable but need to be addressed honestly."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First: this is a solid continuation of Part I, not a standalone discovery. The new material is the closed-form passage from the condensed SGE to the standard C–M–A form, the symmetry classification (C isotropic, M Z3, A D6), and the validation against lattice solutions under simple shear and uniaxial strain. The non-centrosymmetric term m is genuinely activated in the lattice response, and the transverse displacement it produces is a nice check.\n\nWhat the paper does well: the underdetermination is handled honestly. After imposing symmetry and equilibrium, four coefficients remain free, and the authors fix them with an explicit rule: they drop the equilibrium constraint on the identification fields. That rule is a choice, not a derivation, but it is clearly stated and the consequences are shown. The algebra is heavy but internally coherent. The validation shows a real reduction in energy mismatch relative to Cauchy in most cases, and the sign change of the transverse displacement with the permutation of the two non-perimeter springs is persuasive evidence that the non-centrosymmetric effect is physical.\n\nWhere it is soft: the chosen relaxation is not derived from lattice mechanics, and the validation exercises only loadings aligned with the lattice axes. That does not test the two non-condensed curvature modes, so the claim of general energetic equivalence is not established. An off-axis test would help. The tables support an improvement, not 'excellent agreement': errors go up to 19.3% in shear (Ex2, x2-direction), and the uniaxial cases include negative mismatches, meaning the SGE is stiffer than the lattice. Eq. (71) uses a coefficient a45 that never appears in the definitions. The equality of positive-definiteness domains between condensed and standard forms is stated without proof. The coupled BVP solutions are omitted, which limits reproducibility.\n\nBottom line: the paper is a useful contribution to lattice-to-continuum identification, and the methodological gap is real but not disqualifying. The central claim should be read as holding under the chosen relaxation and for the tested loadings. It deserves a serious referee; I would not desk reject it. The authors should define a45, temper the validation wording, and ideally test an off-axis case or provide the code.","headline":"Useful closed-form SGE identification from a hexagonal lattice with a transparent but underdetermined continuation rule; the validation is decent but stops short of 'excellent'.","tokens_in":24830,"tokens_out":5551,"would_cite":true,"duration_ms":50657,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper establishes that a periodic hexagonal lattice of axially deformable bars with three stiffnesses is, at scales large relative to its cell, energetically equivalent to a homogeneous second-gradient elastic continuum given by…","keywords":["strain gradient elasticity","non-centrosymmetric material","internal length","homogenization","hexagonal lattice","second-gradient elasticity","lattice energy equivalence","simple shear validation"],"falsifier":"Subject a finite hexagonal lattice strip and its equivalent second-gradient continuum to three-point bending or a concentrated edge load, where the displacement field is not the quadratic equilibrium family used in the identification, and compare cell-averaged displacements and stored energies; if the SGE mismatch is no better than the Cauchy mismatch, the sufficiency assumption that the relaxed model works for arbitrary fields is refuted.","tokens_in":23743,"feed_emoji":"🕸️","tokens_out":16116,"duration_ms":141471,"temperature":0.7,"pith_summary":"The paper completes a two-part identification programme by turning a discrete hexagonal lattice, built from hinged bars with three different axial stiffnesses, into a homogeneous second-gradient elastic continuum written in closed form through three constitutive matrices. At first order the lattice is just an isotropic Cauchy solid, but at second order the equivalent material acquires an internal length, a strain-curvature coupling that makes it non-centrosymmetric, and higher-order anisotropy, even though the local response remains isotropic. The authors validate the model by solving simple shear and uniaxial strain problems exactly for the lattice and for the continuum, showing that the second-gradient material tracks the exact lattice energy far more closely than the classical Cauchy model. They also map the bar-stiffness ratios for which the equivalent energy is positive definite, and show that even an indefinite equivalent can have positive stored energy once the body is large enough compared with the cell size.","feed_headline":"Honeycomb lattice acts like a solid with an internal length","feed_subtitle":"Closed-form constitutive matrices reproduce shear and uniaxial strain far better than the classical Cauchy model.","key_machinery":"The central object is the trio of constitutive matrices $\\mathbf{C}(k,\\hat{k},\\tilde{k})$, $\\mathbf{M}(k,\\hat{k},\\tilde{k})$, and $\\mathbf{A}(k,\\hat{k},\\tilde{k})$ of the standard second-gradient material, obtained from the condensed tensors $\\mathbf{M}^*$ and $\\mathbf{A}^*$ of Part I through the transformation $\\mathbf{q}^* = \\mathbf{Q}\\,\\mathbf{q}_{\\mathrm{SGE}}$, $\\mathbf{M} = \\mathbf{M}^*\\mathbf{Q} + \\Delta\\mathbf{M}$, $\\mathbf{A} = \\mathbf{Q}^T \\mathbf{A}^*\\mathbf{Q} + \\Delta\\mathbf{A}$, where the equilibrium constraints fix only part of $\\mathbf{Q}$, $\\Delta\\mathbf{M}$, and $\\Delta\\mathbf{A}$. The identification leaves four coefficients undetermined ($Q_{35}$, $\\Delta A_{11}$, $\\Delta M_{15}$, $\\Delta M_{16}$); the paper fixes them with the relaxation rule of Eq. (51), which corresponds to dropping the equilibrium constraint at the first identification step and yields closed-form rational expressions in the bar stiffnesses while preserving the positive-definiteness domain. The cell side $\\ell$ enters explicitly, linearly in $\\mathbf{M}$ and quadratically in $\\mathbf{A}$, which is what gives the continuum its internal length. Symmetry analysis assigns $\\mathbf{C}$ to $O(2)$, $\\mathbf{M}$ to $Z_3$, and $\\mathbf{A}$ to $D_6$; these symmetry restrictions together with the relaxation rule determine the remaining entries.","core_discovery":"The paper claims that the 'condensed' second-gradient material identified in Part I—by matching the lattice strain energy under quadratic displacement fields whose stress fields satisfy equilibrium—can be relaxed into a standard second-gradient elastic material with constitutive matrices $\\mathbf{C}$, $\\mathbf{M}$, and $\\mathbf{A}$ written in closed form in Eq. (1). The Cauchy matrix $\\mathbf{C}$ is isotropic with Lamé constants in Eq. (2); the coupling matrix $\\mathbf{M}$ is non-centrosymmetric with $Z_3$ symmetry and vanishes exactly when the two non-perimeter bar stiffnesses are equal or when $k = \\hat{k}\\tilde{k}/(2(\\hat{k}+\\tilde{k}))$; the curvature matrix $\\mathbf{A}$ has $D_6$ symmetry. The resulting material carries an internal length $\\ell$, making it non-local and anisotropic in its higher-order response while locally isotropic. The paper validates this identification by solving simple shear and uniaxial strain boundary value problems exactly for both the lattice and the continuum, reporting that the second-gradient model tracks the exact lattice energy far more closely than the Cauchy model and captures a transverse displacement in uniaxial strain that a Cauchy solid cannot produce.","pith_inferences":["Inference: the same matching procedure could be applied to other periodic truss topologies, such as triangular, kagome, square, or chiral lattices, with the symmetry classes of $\\mathbf{M}$ and $\\mathbf{A}$ encoding the lattice point group, so lattice geometry would translate directly into second-gradient constitutive data.","Inference: the relaxation rule of Eq. (51) is one admissible choice among several, since $Q_{35}$ and $\\Delta A_{11}$ remain undetermined by symmetry; alternative standard materials with identical condensed behaviour exist, and comparing them against the same benchmarks would isolate how much of the reported accuracy comes from the identification rather than from the specific rule.","Inference: if the internal-length picture is right, elastic wave dispersion in a finite hexagonal lattice should deviate from the Cauchy prediction at wavenumbers near $1/\\ell$, with the $\\mathbf{M}$ term coupling shear and dilatational branches; measuring dispersion curves would validate the model in dynamics.","Inference: because swapping $\\hat{k}$ and $\\tilde{k}$ flips the sign of $m$ while leaving the other constitutive matrices unchanged, the same lattice geometry could act as a mechanically switchable non-centrosymmetric medium whose bending direction under uniaxial strain is set by the stiffness ordering."],"forward_implications":["A finite hexagonal truss can be replaced by the closed-form continuum of Eq. (1) at scales a few times the cell size, cutting computational cost while retaining size-dependent behaviour through the internal length $\\ell$.","The model predicts a measurable non-centrosymmetric effect: under uniaxial strain along one axis, the lattice develops a transverse displacement whenever $\\hat{k}\\ne\\tilde{k}$ and $k\\ne\\hat{k}\\tilde{k}/(2(\\hat{k}+\\tilde{k}))$, an effect invisible to a Cauchy or centrosymmetric continuum.","Strain-energy comparisons in Tables 2 and 3 show the equivalent SGE material reduces the error from roughly 8-27% for the Cauchy model to about 1-9% in shear and below about 2% for uniaxial strain along the cell direction.","Because the higher-order energy scales as $(\\ell/\\rho)^2$ relative to the classical energy, larger specimens should behave like the Cauchy solid while small specimens should show non-classical effects; the crossover ratio $\\rho/\\ell$ depends on the bar stiffness ratios.","Positive definiteness of the equivalent material is not guaranteed: some stiffness ratios, including the regular equal-stiffness honeycomb, give an indefinite energy, but the energy stored in a finite disk becomes positive once the disk radius is large enough, so finite-size stability is the physically relevant check."],"supporting_citations":[{"why":"Supplies the condensed constitutive matrices M* and A* identified in Part I by energy matching on quadratic equilibrated displacement fields; the standard SGE material in this paper is derived from them.","marker":"[28]"},{"why":"Provide the transformation rules for strain and curvature tensors under rotations and reflections used to assign the O(2), Z3, and D6 symmetry classes.","marker":"[1, 2]"},{"why":"Define the isotropic Cauchy moduli to which the equivalent material reduces when the cell size ell vanishes, giving the first-order baseline for the higher-order comparisons.","marker":"[12, 31]"},{"why":"Earlier results establishing that higher-order solids equivalent to heterogeneous Cauchy materials are not always positive definite, motivating the positive-definiteness analysis and the finite-domain energy check.","marker":"[3, 4, 9, 34]"}],"fun_headline_variants":["Hexagonal lattice maps to a second-gradient solid with a hidden length","Lattice bars yield a continuum with an internal length scale","Second-gradient model captures shears and strains classical misses","From honeycomb bars to non-local elasticity: a validated route","Honeycomb lattice: an internal length emerges from bar interactions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that matching the lattice energy on the restricted family of equilibrium quadratic displacement fields, together with the specific rule chosen for the undetermined coefficients in Eq. (51), fixes a standard second-gradient material that also matches the lattice under arbitrary displacement fields; if that sufficiency step fails, the four validation problems would not be a genuine test.","fun_headline_variants_meta":{"raw":{"variants":["Hexagonal lattice maps to a second-gradient solid with a hidden length","Lattice bars yield a continuum with an internal length scale","Second-gradient model captures shears and strains classical misses","From honeycomb bars to non-local elasticity: a validated route","Honeycomb lattice: an internal length emerges from bar interactions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000201,"raw_usage":{"total_tokens":1397,"prompt_tokens":982,"completion_tokens":415,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":598,"completion_tokens_details":{"reasoning_tokens":331}},"tokens_in":598,"tokens_out":415,"duration_ms":4755,"temperature":1.0,"reasoning_tokens":331,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:06:41.882574+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Subject a finite hexagonal lattice strip and its equivalent second-gradient continuum to three-point bending or a concentrated edge load, where the displacement field is not the quadratic equilibrium family used in the identification, and compare cell-averaged displacements and stored energies; if the SGE mismatch is no better than the Cauchy mismatch, the sufficiency assumption that the relaxed model works for arbitrary fields is refuted.","supporting_citations":[{"cited_title":"Rizzi, D","cited_arxiv_id":null,"evidence_quote":"Supplies the condensed constitutive matrices M* and A* identified in Part I by energy matching on quadratic equilibrated displacement fields; the standard SGE material in this paper is derived from them."}],"review_version":1}