Second-Order Rigidity: Theory & Applications

• Second-order rigidity is a property of discrete constraint systems where all nontrivial first-order flexes are blocked at the second order, ensuring local rigidity.
• It refines linear (first-order) rigidity by using dual criteria through self-stresses and energetic analyses, bridging the gap to prestress stability.
• Practical applications include bar-and-joint frameworks, rigid origami, and mechanical networks where energy-based methods confirm quartic growth and local stability.

Second-order rigidity is a geometric and energetic property of discrete constraint systems—such as bar-and-joint frameworks, origami crease patterns, and geometric incidence structures—which captures a class of local obstructions to finite deformations. A system is second-order rigid if every nontrivial first-order infinitesimal flex is obstructed—i.e., cannot be extended to an actual finite flex—at second order in a Taylor expansion of the constraint equations. This concept refines first-order (infinitesimal) rigidity by detecting hidden nonlinear rigidity not accessible to linear analyses, and is strictly weaker than prestress stability but still sufficient to guarantee local rigidity in many contexts. Second-order rigidity can be characterized via duality with self-stresses and is foundational in energy-based approaches to rigidity, with direct implications for mechanical stability in underconstrained or prestressed materials.

1. Definitions and Formalism

Second-order rigidity is defined for systems governed by a finite set of analytic constraints $f_\alpha(x) = 0$ in configuration space. The Taylor expansion of $f_\alpha$ along a path $x(t) = x_0 + v t + w t^2 + O(t^3)$ yields:

• First-order flex: A velocity $v$ is a first-order flex if $Df(x_0)\,v = 0$, i.e., it preserves the constraints to linear order.
• Second-order flex: A pair $(v, w)$ is a second-order flex if

$Df(x_0)\,w + \frac{1}{2} D^2 f(x_0)[v,v] = 0,$

i.e., the path can be extended to quadratic order while preserving the constraints.

• Second-order rigidity: The system is second-order rigid if the only solution to these equations with $v$ nontrivial is the trivial (rigid motion) flex (Gortler et al., 3 Jun 2025, Connelly et al., 2015).

In bar-and-joint frameworks, for configuration $p \in (\mathbb{R}^d)^n$, the edge-length constraints yield a rigidity matrix $R(p)$ and Hessian. In rigid origami, the folding angles $f_\alpha$0 are constrained by vertex and hole closure conditions $f_\alpha$1, with analogous first- and second-order analyses (He et al., 2020).

The dual criterion uses equilibrium stresses (cokernel elements of the rigidity matrix): a first-order flex $f_\alpha$2 can be extended to a second-order flex if and only if all stress-induced quadratic forms vanish: $f_\alpha$3 Second-order rigidity holds if for every nontrivial $f_\alpha$4 there exists some $f_\alpha$5 such that the above form is positive (He et al., 2020, Connelly et al., 2015, Damavandi et al., 2021).

2. Hierarchy and Duality: First-Order, Second-Order, and Prestress Stability

There is a strict implication chain among three central notions:

Concept	Condition	Implied by	Implies
First-order rigidity	All first-order flexes are trivial	—	Prestress stability
Prestress stability	There exists $f_\alpha$6 s.t. energy quadratic form is PD	First-order rigidity	Second-order rigidity
Second-order rigidity	All first-order flexes are blocked at second-order	Prestress stability	Local rigidity

Second-order rigidity sits strictly between prestress stability and mere local (continuous) rigidity. Prestress stability requires that some linear combination of self-stresses renders the quadratic form strictly positive-definite on the space of first-order flexes; second-order rigidity only demands, for each first-order flex, existence of some stress to block that flex (Connelly et al., 2015, He et al., 2020, Damavandi et al., 2021).

Second-order rigidity is also equivalent to "weak prestress stability": for each nontrivial first-order flex, there exists a stress which penalizes that flex, but not necessarily the same stress for all flexes (Power, 2022).

None of the above implications are reversible in general; examples are known where a system is second-order rigid but not prestress-stable, or locally rigid but not second-order rigid (He et al., 2020, Connelly et al., 2015).

3. Energetic and Analytical Criteria

Second-order rigidity can be rigorously detected via the energetic approach, by studying the growth order of a suitable energy function $f_\alpha$7 associated with the constraints: $f_\alpha$8 At an equilibrium, $f_\alpha$9, and the Hessian $x(t) = x_0 + v t + w t^2 + O(t^3)$0 may be degenerate along first-order flexes. The leading order of energy increase, in a direction of a first-order flex, depends on whether it can be extended to a second-order flex:

• If no nontrivial $x(t) = x_0 + v t + w t^2 + O(t^3)$1 solves the appropriate second-order system, energy grows at least quartically in $x(t) = x_0 + v t + w t^2 + O(t^3)$2 along all nontrivial infinitesimal directions (order $x(t) = x_0 + v t + w t^2 + O(t^3)$3), certifying a strict local minimum and hence rigidity (Gortler et al., 3 Jun 2025, Damavandi et al., 2021, Damavandi et al., 2021).
• If a higher-order flex exists, the corresponding energy growth is slower, and rigidity fails.

In the context of underconstrained or unstressed frameworks ($x(t) = x_0 + v t + w t^2 + O(t^3)$4), the Hessian of energy is positive semidefinite with a nontrivial nullspace: only second-order rigidity can guarantee a positive energy barrier against finite deformations. In cases with a unique positive self-stress, second-order rigidity is both necessary and sufficient for energetic rigidity (Damavandi et al., 2021).

A practical computational criterion is the positivity of a family of quadratic forms $x(t) = x_0 + v t + w t^2 + O(t^3)$5 associated with each self-stress basis vector $x(t) = x_0 + v t + w t^2 + O(t^3)$6. If, for every nontrivial flex $x(t) = x_0 + v t + w t^2 + O(t^3)$7, at least one $x(t) = x_0 + v t + w t^2 + O(t^3)$8, the system is second-order rigid (He et al., 2020).

In energy-based settings involving polyhedral origami, the order of vanishing of $x(t) = x_0 + v t + w t^2 + O(t^3)$9 along curves tracks the order of rigidity: second-order rigidity implies quartic energy growth (He, 27 Jun 2025).

4. Applications and Illustrative Examples

Second-order rigidity has deep consequences across discrete geometry and the mechanics of complex materials:

• Rigid origami: Planar degree-3 vertices and certain non-planar holes have linearized (first-order) flexes but are blocked at second order due to nontrivial self-stresses. Triangulated convex polyhedra are typically second-order rigid—and often prestress stable as well (He et al., 2020, Connelly et al., 2015).
• Bar-and-joint frameworks: Flat triangular prisms may admit infinitesimal flexes but are rigidified at second order; their energy landscape exhibits a strict local minimum only after quartic terms are included (Gortler et al., 3 Jun 2025).
• Mechanical/biological networks: Underconstrained spring networks, vertex models, and hypostatic packings can be second-order rigid (and hence energetically rigid at quartic order) even with a macroscopically vanishing linear shear modulus, due to emergent self-stress penalizing nonlinear flexes (Damavandi et al., 2021, Damavandi et al., 2021).
• Projective rigidity: Incidence structures in projective geometry can be flexible at the infinitesimal level yet second-order rigid, as in mechanisms where a hyperbola constraint allows nontrivial first-order motions that are non-extendable due to higher-order tangency (Berman et al., 10 Mar 2025).

In all cases, the detection and construction of suitable equilibrium stresses is central to certifying second-order rigidity.

5. Second-Order Rigidity in Infinite and Higher-Order Settings

For infinite frameworks (countably many joints and bars), second-order rigidity extends via asymptotic equilibrium stresses and suitably bounded forms of prestress stability. The equivalence between second-order rigidity and "weak prestress stability" holds for infinite frameworks under uniform boundedness hypotheses, providing sufficient conditions for the absence of directed continuous motions (Power, 2022).

Rigidity notions can also be extended to third or higher order, with the energetic growth order of the associated energy function diagnosing whether higher order flexes exist. Notably, there are no fractional rigidity orders between first and second: flexibility failures at lower order cannot bridge the gap between linear and quadratic energy vanishing (He, 27 Jun 2025). In frameworks where the nullity of the rigidity matrix is one, the existence of a $v$0-th order flex precisely controls rigidity order (Gortler et al., 3 Jun 2025).

6. Computation and Testing

Numerical and symbolic procedures for testing second-order rigidity typically involve:

1. Assembling the rigidity matrix and computing its nullspace (first-order flexes) and its left nullspace (self-stresses).
2. Forming the appropriate second-order tensors (Hessian or constraint second derivatives) and quadratic forms induced by self-stress directions.
3. Restricting these quadratic forms to the first-order flex space and verifying, for all nontrivial flexes, the existence of a self-stress making the form positive.
4. Formulating the positivity test as a semidefinite feasibility problem, sometimes solvable by convex optimization (He et al., 2020).

This approach is effective in bar-joint frameworks, rigid origami, and polyhedral constraint systems, with higher computational complexity than first-order tests but yielding finer resolution of marginal rigidity (He, 27 Jun 2025).

7. Broader Significance and Open Directions

Second-order rigidity provides a bridge between combinatorial/geometric rigidity and energetic/mechanical stability, especially in highly redundant or underconstrained networked systems. It underlies phenomena where apparent linearized floppiness paradoxically gives rise to nonlinear robustness and has been directly linked to emergent material rigidity in biological and synthetic systems with vanishing shear modulus at linear order. In the theory of geometric PDEs, second-order rigidity conjectures (e.g., for Einstein–Weyl conformal structures) underpin contact-equivalence classification of integrable classes (Berjawi et al., 2021).

Open problems include intrinsic algebraic characterization of second-order rigidity (for instance, via invariant polynomials), computational scalability for large or infinite systems, and generalized frameworks for higher-order rigidity and associated energy landscapes (Berman et al., 10 Mar 2025, Gortler et al., 3 Jun 2025, He, 27 Jun 2025).