General up to next-nearest neighbour elasticity of triangular lattices in three dimensions

Cyril Dubus (MSC; Jean-Baptiste Fournier (MSC; Ken Sekimoto (MSC; LPCT)

REVIEW

General up to next-nearest neighbour elasticity of triangular lattices in three dimensions

Cyril Dubus (MSC, LPCT), Ken Sekimoto (MSC, Jean-Baptiste Fournier (MSC

Soft Condensed Matter · Submitted 21 Feb 2006 · 2 references

Pith's verdict: Real result and a useful resolution of the Keating/Feng-Sen question, but the five-constant count rests on an unproven linear-independence claim that needs verification.

T0 · AI review

How this review was made: An AI referee (deepseek-v4-flash) read the paper on August 28, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

T0 AI review An AI read the paper and answered a fixed list of questions.
T1 Audit Automatic checks that the references and files the paper relies on exist and match.
T2 Compute The paper's code or calculations re-run with the same results.
T3 Formal statement The claims written out in a formal math language a computer can read.
T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 5 major objections · 6 minor · Ask about this paper · The compiled record · Read the full review →

Challenge this review Comment (0) Read on arXiv ↗

Cite ▾

arXiv IDcond-mat/0602494 (version 1)

CitationCyril Dubus (MSC, Jean-Baptiste Fournier (MSC, Ken Sekimoto (MSC, LPCT). "General up to next-nearest neighbour elasticity of triangular lattices in three dimensions." arXiv:cond-mat/0602494 (2006). https://pith.science/paper/cond-mat/0602494

Cite this reviewBibTeX for Pith's review of this paper

Claim paper

Why Pith says this

Highlight any text to comment on it.

Editor's note

read the letter

Short version: the paper is worth reading and probably right, but I would not fully trust the five-constant count until the linear-independence assertion in §3(b) is proved or checked numerically. The central equation (3.2) gives a clean, compact form for the most general harmonic energy of a triangular lattice in 3D with up to next-nearest neighbors, and the key conceptual claim—that boundary terms are mandatory for rotational invariance and that the Kt term is admissible—is well argued and, as far as I can tell, correct.

The paper does several things right. The symmetry-based enumeration is systematic and the reduction from the seven rigidities in figure 4 to five bulk constants is genuinely new. The continuum limit is a nice external check: for central-force springs it recovers λ = μ, which matches the classical result. The discussion of why edge terms contribute at order ω² under a rotation is clear and convincing. If I worked on cytoskeleton or gel elasticity, I would want to know this result.

The soft spots are concentrated in the proof apparatus. The five quadratic forms in (3.2) are said to be linearly independent “in their functional space,” but no proof or reference is given. That is the exact step that makes the five-constant statement. If two in-plane forms are dependent, the count shrinks; if the π/6 angle terms are not in the span, it grows. Either way the headline claim changes. So this is load-bearing, not cosmetic. The constraints D_i = 0, E_i = 0, G*_i = 0 are also waved at with “easy to show,” though the appendix does include a fair amount of algebra. The vertex classification into five types and eight doublets is done for one specific regular-hexagon shape; the “most general” claim assumes that classification is exhaustive, which is plausible but not proven.

None of this makes me think the result is wrong. The derivation is parameter-free, the logic is coherent, and the admissibility of Kt is supported by the explicit boundary analysis. But the missing proof is a genuine reproducibility gap. A rank test on a periodic cluster would settle it in an afternoon, and I would want that or a detailed proof before using the five-constant count as a foundation.

The paper deserves a serious referee and probably publication after revision. Send it to peer review, ask the authors to make the independence argument explicit and to expand the thin parts of the appendix. People in soft matter theory and lattice elasticity will benefit from having this result on the record.

Referee report

Significance

If the central claim holds, the paper resolves a longstanding controversy about admissible non-central terms in triangular-lattice elasticity and provides a systematic, parameter-free procedure for deriving lattice elasticity. The continuum check against the known central-force result λ = μ (Kantor et al.) is a useful external benchmark, and the mapping to measurable elastic constants is practically relevant for modeling cytoskeletons, percolating gels, and membranes. The main result, however, rests on unproven algebraic and classification claims; the five-constant count and the very form of Eq. (3.2) are not yet established with full rigor. The potential significance is high, but the current version is conditional on supplying the missing demonstrations.

Major objections (5)

[§3(b), Eq. (3.2)] The statement after Eq. (3.2) that the five quadratic forms 'can be shown to be linearly independent in their functional space' is given without proof. This assertion is the exact step that converts the seven (or more) rigidities of Fig. 4 into the headline claim of only five independent bulk elastic constants. If two of the forms are linearly dependent, the true count is smaller; if the π/6 angle terms are not in the span of the five forms, the count is larger. Because no derivation or machine-checked rank test is provided, the central quantitative claim is not yet verified. Please provide an explicit proof or an algorithmic verification (e.g., a rank computation on a finite cluster with periodic boundary conditions).
[§3(a), Eq. (3.1)] The identity (3.1), which expresses sums of π/6-angle deviations in type-IV clusters as a linear combination of π/3-type and 2π/3-type terms, is asserted without derivation. This identity is load-bearing because it is used to argue that no additional independent bulk term associated with the π/6 angles appears. If (3.1) is incorrect, the five-constant form (3.2) is incomplete. Please provide the full derivation of (3.1).
[Appendix §b(iii), §c(ii), §d(ii)] The translational and rotational invariance constraints Di = 0, \bar Di = 0, Ei = 0, G = 0, and Gi* = 0 are each justified only by 'it is easy to show' or 'properly chosen distorted states', without specifying the distorted states or carrying out the algebra. These constraints are essential to the final energy form (3.2). Since the entire derivation depends on these eliminations, please supply the complete calculations, including the explicit displacement patterns used and the resulting linear equations for the coefficients.
[§2(c)(i)] The classification of vertices into five types and vertex doublets into eight types is performed only for the geometry of a regular hexagonal lattice. The paper claims the 'most general' elasticity of triangular lattices, but it does not demonstrate that this classification is exhaustive for arbitrary boundary shapes. If a different boundary shape produces additional vertex or doublet environments, the rotational-invariance conditions could impose further constraints on the bulk coefficients, changing the five-constant result. Please clarify whether the bulk result is independent of the boundary shape, or restrict the claim to regular hexagonal lattices.
[§2(a), Abstract] The derivation explicitly assumes in §2(a) that the elastic properties are symmetric with respect to the lattice plane P, which excludes any coupling between in-plane and normal displacements at harmonic order. The abstract's unqualified claim of the 'most general form' is therefore too strong. The mirror-symmetry assumption should be stated in the abstract, and the result presented as 'most general under the assumption of reflection symmetry with respect to P'.

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

Triangular-lattice elasticity with next-nearest-neighbour couplings reduces to five independent bulk constants, and a non-central spring term once thought to break rotational symmetry is allowed.

The reading

This paper asks what the most general quadratic elastic energy of a two-dimensional triangular lattice embedded in three dimensions can be, when interactions reach only to next-nearest neighbours. Its answer is that the bulk part of that energy has exactly five independent constants, provided the energy stored at the lattice edges is included and used to enforce global rotational invariance. The compact form is equation (3.2): nearest-neighbour pairs contribute normal, transverse, and longitudinal stiffnesses $K_n$, $K_t$, $K_\ell$, while next-nearest-neighbour pairs contribute $-K_n/3$, $K'_t$, and $K'_\ell$, plus edge terms. A central consequence is that the non-central nearest-neighbour term with coefficient $K_t$, which earlier work in the percolation-network literature took to violate rotational invariance, is in fact admissible. This matters for any model of cytoskeletal networks, polymerized membranes, or disordered gels built on a triangular lattice, because it fixes which elastic parameters are independent and measurable.

What carries the argument

The machinery is a classification of local environments. For the interaction range fixed at the next-nearest-neighbour distance, vertices can be grouped into seven neighbour clusters; by their environment relative to those clusters the edge vertices fall into five types and oriented nearest- and next-nearest-neighbour vertex pairs into eight doublet types. Each type has its own symmetry-allowed quadratic terms. The paper then rewrites edge sums as bulk sums using a discrete Stokes-type identity, and imposes translational and rotational invariance on the full energy rather than on the bulk alone. This apparatus is what makes the enumeration exhaustive and forces the boundary terms that legitimise the $K_t$ term.

Extended reading notes

Core claim

The central result is equation (3.2). For a regular hexagonal triangular lattice lying in a plane, with elastic properties mirror-symmetric about that plane, the most general harmonic energy with nearest- and next-nearest-neighbour interactions that is invariant under global translations and rotations takes the form $$F = \frac{1}{2}\sum_{\langle\$\alpha$\$\beta$\rangle_I}\left\{K_n(u_n^\$\alpha$-u_n^\$\beta$)^2 + K_t[(u_\perp^\$\alpha$-u_\perp^\$\beta$)\cdot t_{\perp}^{\$\alpha$\$\beta$}]^2 + K_\ell[(u_\perp^\$\alpha$-u_\perp^\$\beta$)\cdot $t^{{\alpha\beta}}$]^2\right\} + \frac{1}{2}\sum_{\langle\$\alpha$\$\beta$\rangle_{II}}\left\{-\frac{K_n}{3}(u_n^\$\alpha$-u_n^\$\beta$)^2 + K'_t[\cdots]^2 + K'_\ell[\cdots]^2\right\} + \text{edge terms},$$ with only the five bulk constants $K_n$, $K_t$, $K_\ell$, $K'_t$, and $K'_\ell$ independent. The derivation shows that the boundary elasticity cannot be ignored: terms living on the edges cancel the non-invariant parts of bulk terms under rigid rotations, so the manifestly invariant bulk-only expression with $K_t$ is legitimate. This rescues the non-central nearest-neighbour term that invariance arguments in the literature had been taken to forbid and that percolation-network models had consequently excluded.

Load-bearing premise

The derivation assumes the elastic properties are symmetric with respect to the lattice plane and that the five vertex types and eight doublet types found on a regular hexagonal boundary exhaust every possible local environment; if either fails, the five-constant bulk form need not be the most general.

Editorial extensions

If this is right

Any harmonic triangular-lattice model with interactions up to next-nearest neighbours can be re-expressed in the five-constant form, so any apparent additional bulk rigidity is a combination of those five or an edge effect.
The non-central nearest-neighbour term $K_t$ can be kept in models of percolating gels or cytoskeletons without violating rotational invariance; it is an independent parameter rather than an artifact.
In the continuum limit only three combinations survive, the two Lamé coefficients and the bending modulus, so the discrete description has two extra coefficients visible only at finite lattice scale.
Stability of the flat lattice does not require all five bulk constants to be positive; a negative $K'_t$ is possible while the lattice remains stable.
Omitting edge terms makes otherwise allowed bulk terms fail rotational invariance, so simulations and fits of finite triangular networks should keep the boundary elasticity explicit.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

The neighbour-cluster enumeration should transfer to other lattices and interaction ranges, giving an automatic way to list the independent harmonic constants of any discrete network with a well-defined local environment.
Since the two extra discrete constants vanish in the long-wavelength limit, fitting phonon spectra or fluctuation modes of a triangulated cytoskeleton at wavelengths of a few lattice spacings could separately determine $K_t$ and $K'_t$.
The boundary-term cancellation is size-dependent in origin, so small patches may show apparent bulk moduli that drift with patch size if edge terms are omitted; measuring that drift would test the construction directly.

Original abstract

We establish the most general form of the discrete elasticity of a 2D triangular lattice embedded in three dimensions, taking into account up to next-nearest neighbour interactions. Besides crystalline system, this is relevant to biological physics (e.g., red blood cell cytoskeleton) and soft matter (e.g., percolating gels, etc.). In order to correctly impose the rotational invariance of the bulk terms, it turns out to be necessary to take into account explicitly the elasticity associated with the vertices located at the edges of the lattice. We find that some terms that were suspected in the litterature to violate rotational symmetry are in fact admissible

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 3,636 input tokens, $0.0002 reported.

What the abstract itself states, as probabilities: claims first 0.13 · states limitation 0.40 · human participants 0.03 · controlled comparison 0.02 · randomized assignment 0.01 · states quantitative result 0.06 · compares with prior results 0.43 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; show every role.

We establish the most general form of the discrete elasticity of a 2D triangular lattice embedded in three dimensions, taking into account up to next-nearest neighbour interactions. Besides crystalline system, this is relevant to biological physics (e.g., red blood cell cytoskeleton) and soft matter (e.g., percolating gels, etc.). In order to correctly impose the rotational invariance of the bulk terms, it turns out to be necessary to take into account explicitly the elasticity associated with the vertices located at the edges of the lattice. We find that some terms that were suspected in the litterature to violate rotational symmetry are in fact admissible

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#	sentence	role	the authors' own work	states the main conclusion	strength
s0	We establish the most general form of the discrete elasticity of a 2D triangular lattice embedded in three dimensions, taking into account up to next-nearest neighbour interactions.	contribution 0.97	0.96	0.73	asserted 1.00
s1	Besides crystalline system, this is relevant to biological physics (e.g., red blood cell cytoskeleton) and soft matter (e.g., percolating gels, etc.).	background 0.65	0.09	0.04	asserted 0.96
s2	In order to correctly impose the rotational invariance of the bulk terms, it turns out to be necessary to take into account explicitly the elasticity associated with the vertices located at the edges of the lattice.	result 0.48	0.82	0.16	asserted 1.00
s3	We find that some terms that were suspected in the litterature to violate rotational symmetry are in fact admissible	result 0.99	0.94	0.27	asserted 1.00

What kind of paper the classifier reports this is: theoretical result 1.00 · evidence from proof or derivation 1.00.

What was examined, and what was not

sources read

arXiv listing abstract v1

not examined

full text
figures and tables
references and citing sentences
relations to other papers

levels run

P abstract

not run

S sections
T body sentences
C citations
R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures

**Figure 1.** Triangular elastic lattice with the vertex

**Figure 2.** The seven types of neighbour clusters as indicated

**Figure 3.** The 5 types of vertices (left) and the 8 types of vert

**Figure 4.** Interpretation of the elastic terms in (2.22) in te

**Figure 5.** Number of vertices, Vi, and of pair of vertices, Pi, as a function of their type (see figure 3) for a regular hexagon of edge length Nℓ. (r.-h. s.), which is a boundary term, can be shown to yield an extensive contribution −2ω 2 A, where A is the area of the lattice, just as the first term of the r.-h. s, which is a bulk term. It follows that: i) this first term is not rotationally invariant by itself, ii) that boundary terms are therefore necessary to correctly ensure rotational invariance, and that, as announced, iii) the term with coefficient Kt is indeed allowed, even if we…

Figures reproduced from arXiv:cond-mat/0602494 by the authors.

References

Works cited (2)

[1]

Arbabi, S., Sahimi, M., 1993 Mechanics of disordered solids . I. Percolation on elastic networks with central forces. Phys.Rev. B 47, 695-702. Born, M. & von K´ arm´ an, T. 1912 ¨Uber Schwingungen in Raumgittern. Physik. Z. 13, 297-309. Born, K. 1914 Zur Raumgittertheorie des Diamanten. Ann. Physik 44,

View on Pith 1993
[605]

Born, M & Huang, K. 1954 Dynamical theory of crystal lattices , ch. V, §23. Oxford Uni- versity Press, New York. Coughlin, M. F., Stamenovic, D. 2003 A Prestressed Cable Net work Model of the Adherent Cell Cytoskeleton. Biophys. J 84, 1328 Discher, D. E., Boal, D. H. & Boey, S. K. 1997 Phase transition s and anisotropic responses of planar triangular nets...

View on Pith 1954

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.

[1] [1]

Arbabi, S., Sahimi, M., 1993 Mechanics of disordered solids . I. Percolation on elastic networks with central forces. Phys.Rev. B 47, 695-702. Born, M. & von K´ arm´ an, T. 1912 ¨Uber Schwingungen in Raumgittern. Physik. Z. 13, 297-309. Born, K. 1914 Zur Raumgittertheorie des Diamanten. Ann. Physik 44,

View on Pith 1993

[2] [605]

Born, M & Huang, K. 1954 Dynamical theory of crystal lattices , ch. V, §23. Oxford Uni- versity Press, New York. Coughlin, M. F., Stamenovic, D. 2003 A Prestressed Cable Net work Model of the Adherent Cell Cytoskeleton. Biophys. J 84, 1328 Discher, D. E., Boal, D. H. & Boey, S. K. 1997 Phase transition s and anisotropic responses of planar triangular nets...

View on Pith 1954