Elastic: Multidisciplinary Concepts

• Elastic is a phenomenon defined by recoverable deformation in materials, characterized by load-strain relations in isotropic and critically elastic systems.
• Recent studies extend elasticity to encompass odd, spin, and conformal variants, revealing nonconservative behaviors, tunable responses, and applications in metamaterials.
• Elastic concepts also underpin energy minimization in network science, optimization algorithms, and even astrophysical models, highlighting its versatility across disciplines.

“Elastic” most commonly denotes recoverable deformation and the constitutive relations that connect load, strain, and stored energy, but contemporary research uses the term across a much wider technical range. In condensed matter and mechanics it spans conventional linear elasticity, critically elastic disordered media, active and nonreciprocal solids, hydroelastic waves, and elastic turbulence; in geometry it labels variational theories of strips and curves; in computation it appears in first-principles tensor extraction, automated toolkits, and optimization methods inspired by the elastic net; and in network science it denotes robustness metrics and energy-like graph formalisms (Hagh et al., 2022, Scheibner et al., 2019, Chubelaschwili et al., 2010, Bates et al., 2015, Lin et al., 2024).

1. Recoverable deformation and critically elastic matter

In its conventional mechanical sense, elasticity concerns the response of a medium that deforms under load and recovers when the load is removed. In isotropic media, bulk and shear moduli provide the standard coarse-grained description of resistance to volumetric and shape deformations. A particularly constrained regime arises in critically elastic materials, defined as mechanically rigid systems possessing exactly one state of self-stress (SSS). These systems can be generated from a parent with two SSS by removing a single constraint or bond (Hagh et al., 2022).

For such daughter systems, the bulk modulus $K'$ and shear modulus $G'$ are not independent. Instead, the normalized pair $(K'/K,\, G'/G)$ lies on a universal ellipse parameterized entirely by properties of the parent system. The parent-dependent parameter is

$$\eta=\frac{(a_1 b_2-a_2 b_1)^2}{(a_1^2+a_2^2)(b_1^2+b_2^2)},$$

where $a_i$ and $b_i$ are projections of parent SSS vectors onto bulk and shear affine deformations. The ellipse can be oriented so that the two daughter moduli are either positively correlated or negatively correlated, dividing parent systems into two behavioral classes (Hagh et al., 2022).

The framework was verified in simulations of spring networks and soft-sphere packings near marginal rigidity. Varying the parent geometry or contact arrangement changes $\eta$, and therefore changes the accessible elliptic locus of daughter moduli. This enables deterministic modulus selection by bond pruning: for a desired $K'$ or $G'$, there can be up to two distinct removable bonds that produce the target daughter modulus. In two dimensions, the corresponding Poisson ratio is

$\nu'=\frac{K'-G'}{K'+G'}.$

The same deterministic restriction can therefore drive sign changes in Poisson ratio, including transitions between jammed and auxetic responses (Hagh et al., 2022).

This critically elastic setting is unusual because it replaces the generic independence of elastic constants with a parent-controlled functional interrelation. A plausible implication is that near-marginal disordered media provide a design regime in which microscopic constraint editing is sufficient to prescribe macroscopic isotropic response with unusually fine control.

2. Nonconservative, spin, and conformal generalizations

A major recent extension of elasticity relaxes the assumption that stress derives from a stored-energy functional. In odd elasticity, the usual major symmetry of the elastic modulus tensor is broken, so the constitutive law contains an antisymmetric component under exchange of index pairs. In a two-dimensional isotropic solid, the conventional bulk modulus $G'$0 and shear modulus $G'$1 are supplemented by two additional moduli, $G'$2 and $G'$3, where $G'$4 couples compression to internal torque density and $G'$5 couples one shear mode to an orthogonal shear stress. Because these odd terms are nonconservative, a closed quasi-static strain cycle can perform nonzero work; the solid behaves as a distributed engine rather than a purely passive elastic body (Scheibner et al., 2019).

The odd-elastic constitutive matrix for a 2D isotropic medium can be written as

$G'$6

Within this framework, activity can drive auxetic behavior, alter Poisson response, and support propagating elastic waves even in overdamped media. The review literature further situates odd elasticity alongside odd viscosity, emphasizing transverse responses, modified dislocation dynamics, topological waves, and microscopic realizations through transverse noncentral forces, active hinges, and nonsymmetric mobility matrices (Fruchart et al., 2022).

A second generalization transfers elasticity from positional degrees of freedom to internal order-parameter space. Spin elasticity treats collective spin textures as objects that deform under spin torques and recover when those torques are removed. For a one-dimensional domain-wall pair assembly, the prototype “spin spring” obeys the Hooke-like law

$G'$7

The theory introduces a spin strain tensor $G'$8, a spin stress tensor $G'$9, and a local constitutive relation $(K'/K,\, G'/G)$0, while emphasizing that the spin modulus tensor depends on geometry, topology, and position within the texture. Reported consequences include a Poisson effect in spin space, elastic oscillations, spin stress waves, and current-driven tuning of the elastic state of spin textures (Gao et al., 23 Mar 2026).

A third extension emerges in mechanism-based metamaterials. Conformal elasticity describes low-energy deformations of designer dilational metamaterials whose soft pathway penalizes shear strongly while permitting local dilation and rotation. In two dimensions, the compatible low-energy deformations are conformal maps,

$(K'/K,\, G'/G)$1

and the compatibility condition reduces to the Cauchy–Riemann structure. Experiments and finite-element simulations on rotating-square architectures show that generic loads can indeed produce nonuniform, nonlinear deformations that are well fit by conformal maps, and the resulting bulk-boundary principle implies that boundary dilation determines the interior deformation field (Czajkowski et al., 2021).

Together, these developments show that “elastic” no longer implies purely conservative Hookean response in Euclidean space. It can also signify nonreciprocal constitutive structure, recoverable deformation in spin morphology, or mechanism-constrained conformal kinematics.

3. Elastic strips and elastica in differential geometry

In differential geometry, “elastic” often refers to variational problems for curves and developable surfaces. For narrow developable strips, the governing functional is the Sadowsky functional

$(K'/K,\, G'/G)$2

where $(K'/K,\, G'/G)$3 is curvature and $(K'/K,\, G'/G)$4 is modified torsion. Elastic strips are critical points of this functional under fixed-length constraints. A central result is the derivation of two conservation laws: a constant force vector

$(K'/K,\, G'/G)$5

and a constant torque vector

$(K'/K,\, G'/G)$6

These first integrals characterize equilibrium and lead to two integrable classes, elastic momentum strips and force-free strips, each connected to spherical elastic curves traced by the binormal or tangent field. The theory also relates force-free strips to Hopf tori and introduces a $(K'/K,\, G'/G)$7-functional that reduces the variational problem in $(K'/K,\, G'/G)$8 to one on the unit sphere (Chubelaschwili et al., 2010).

The related elastica problem studies curves in $(K'/K,\, G'/G)$9 that are stationary for the integral of squared curvature,

$$\eta=\frac{(a_1 b_2-a_2 b_1)^2}{(a_1^2+a_2^2)(b_1^2+b_2^2)},$$0

Written as a second-order variational system, the corresponding Lagrangian is

$$\eta=\frac{(a_1 b_2-a_2 b_1)^2}{(a_1^2+a_2^2)(b_1^2+b_2^2)},$$1

Reparametrization invariance produces conserved quantities through Noether-type arguments, including explicit linear and angular momenta, and leads naturally to a constrained Hamiltonian formulation using Ostrogradski momenta and the Cartan form. The same work develops a Bleuler–Gupta quantization procedure in which admissible quantum states are annihilated by operator versions of the classical constraints (Bates et al., 2015).

These geometric theories use “elastic” in a strict variational sense: equilibrium shapes are selected by curvature-based functionals rather than by continuum stress-strain tensors. The shared theme is recoverable bending encoded through conserved geometric quantities.

4. Waves, turbulence, and nanoscale elastic emission

Elastic phenomena also appear in wave propagation and flow-mediated instabilities. For linearized flow past submerged obstacles under an elastic sheet, with gravitational effects included in two dimensions, exponential asymptotics reveals that the waves are exponentially small in the limit of small bending length relative to obstacle depth. In the two-dimensional problem, the relative strength of elastic and gravitational restoring forces divides behavior into two classes: in one regime, constant-amplitude elastic waves and gravity waves extend indefinitely upstream and downstream; in the other, all waves decay exponentially away from the obstacle. The nonlinear two-dimensional problem predicts a third intermediate regime with one-sided persistent waves, while the three-dimensional geometry yields elastic waves that extend ahead of the submerged source and decay algebraically in space (Lustri, 2022).

At the nanoscale, the elastic Purcell effect provides an elastic analogue of cavity-enhanced electromagnetic emission. A localized elastic source couples to the local density of elastic states through the imaginary part of the dyadic elastic Green function, and resonant elastic nanoparticles act as antennas that modify emission rates. The formalism introduces an elastic Purcell factor and an effective mode volume,

$$\eta=\frac{(a_1 b_2-a_2 b_1)^2}{(a_1^2+a_2^2)(b_1^2+b_2^2)},$$2

For a submicron gold sphere in silicon, low-order shear and mixed modes were found to provide considerable elastic Purcell factors, indicating that phonon emission can be selectively enhanced by resonant elastic modes (Schmidt et al., 2018).

In viscoelastic fluids, elastic turbulence denotes a disordered flow state at vanishing inertia and large elasticity. Extensive numerical simulations of two-dimensional Taylor–Couette flow found that the purely elastic instability is supercritical, with order parameter scaling

$$\eta=\frac{(a_1 b_2-a_2 b_1)^2}{(a_1^2+a_2^2)(b_1^2+b_2^2)},$$3

The fully nonlinear state is weakly anisotropic and strongly nonhomogeneous: fluctuations are confined to a dynamically active region adjacent to the inner wall, with distinct elastic and kinetic boundary layers. This result clarifies earlier claims of subcriticality by showing that previously reported hysteresis can arise from insufficient resolution (Hou et al., 15 Oct 2025).

Across these examples, elasticity governs wave selection, localization, and dissipation control in settings ranging from hydroelastic free surfaces to resonant nanoantennas and inertialess turbulence.

5. Measurement, first-principles calculation, and automated elasticity workflows

High-throughput materials informatics has made elastic tensors and derived observables routine screening targets. A large first-principles study using the vdW-DF-optB88 functional computed elastic properties for 11,067 bulk materials and 257 monolayer materials, and reported that elastic constants generally decrease with reduced dimensionality. The same dataset identified a correlation between low minimum elastic constants and low exfoliation energies in layered materials, offering a prescreen for van der Waals bonding, and predicted several auxetic candidates with negative Poisson ratio (Choudhary et al., 2018).

Automated pipelines now systematize such computations. MyElas is an integrated framework that automates preprocessing, first-principles calculation, extraction of second- and third-order elastic constants by the energy-strain method, derivation of related quantities, and visualization of anisotropy in elastic modulus and sound velocity (Wang et al., 2022). ElasTool automates second-order elastic-constant calculations for any 2D and 3D crystal system using three strain-matrix sets—OHESS, ULICS, and ASESS—and supports both zero- and high-temperature workflows (Liu, 2020). ElATools provides terminal-based analysis of anisotropic elastic properties for 2D and 3D materials, computes bulk, Young’s, shear, and other derived moduli under Voigt, Reuss, and Hill schemes, and integrates an online and offline Materials Project database with more than 13,000 elastic stiffness constants for 3D materials (Yalameha et al., 2021).

Long-wavelength perturbation theory has also produced a more microscopic route to elastic tensors. A recent formulation extracts bulk elastic tensors and 2D bending rigidities from interatomic force constants, while treating macroscopic electric fields and long-range electrostatics separately to enforce the correct electrical boundary conditions. The reported convergence behavior is especially noteworthy: multipolar interactions up to at least octupoles are required for accurate short-circuit elastic tensors in bulk materials, while orders beyond octupole are needed to converge bending-rigidity tensors of 2D crystals (Lin et al., 2024).

Experimental inference methods complement computation. Pendant capsule elastometry fits axisymmetric shell-theory solutions to images of deflated capsules and thereby extracts elastic parameters from shape alone. For nonlinear Hookean elasticity it determines the Young’s surface modulus and Poisson’s ratio; for wrinkled capsules, an additional wavelength measurement yields the bending modulus through

$$\eta=\frac{(a_1 b_2-a_2 b_1)^2}{(a_1^2+a_2^2)(b_1^2+b_2^2)},$$4

from which layer thickness can be derived (Hegemann et al., 2017).

The resulting landscape is one in which elastic characterization is increasingly automated, high-throughput, and multimodal, spanning direct imaging, density-functional workflows, anisotropy analysis, and long-wavelength lattice perturbation theory.

6. Network, graph, optimization, and relativistic usages

Outside continuum mechanics, “elastic” has acquired specialized meanings grounded in robustness, energy minimization, and regularization. In complex-network analysis, elasticity is a topological robustness metric defined as the area under the curve of normalized throughput $$\eta=\frac{(a_1 b_2-a_2 b_1)^2}{(a_1^2+a_2^2)(b_1^2+b_2^2)},$$5 versus the fraction of nodes remaining during failure or attack: $$\eta=\frac{(a_1 b_2-a_2 b_1)^2}{(a_1^2+a_2^2)(b_1^2+b_2^2)},$$6 The metric is normalized to $$\eta=\frac{(a_1 b_2-a_2 b_1)^2}{(a_1^2+a_2^2)(b_1^2+b_2^2)},$$7, applies even after disconnection, and was used to compare Internet and social-network topologies under random and targeted removal. Reported examples include targeted/random elasticity values of $$\eta=\frac{(a_1 b_2-a_2 b_1)^2}{(a_1^2+a_2^2)(b_1^2+b_2^2)},$$8 for Abilene and $$\eta=\frac{(a_1 b_2-a_2 b_1)^2}{(a_1^2+a_2^2)(b_1^2+b_2^2)},$$9 for MySpace, while scale-free networks remained robust to random removal but vulnerable to targeted attacks (0811.4040).

In graph theory, an elastic graph is a finite graph whose edges carry positive elastic constants $a_i$0. For a map $a_i$1 into a length graph, the Dirichlet or elastic energy is

$a_i$2

A homotopy class is loosening if it decreases elastic energy for every target graph. The fundamental theorem identifies three equivalent quantities: $a_i$3 so that loosening is characterized by $a_i$4 (Thurston, 2016).

In statistics and optimization, Elastic Gradient Descent generalizes both gradient descent and forward stagewise regression in analogy with the elastic net. With gradient $a_i$5 and selection matrix $a_i$6, the update direction is

$a_i$7

The method generates sparse solution paths that resemble those of the elastic net, while reported experiments show similar solutions with speedups of up to three orders of magnitude on the investigated data (Allerbo et al., 2022).

Relativistic astrophysics extends the term back to material response, but in strongly gravitating regimes. A general-relativistic framework for self-gravitating elastic matter introduces a stored-energy function $a_i$8 and anisotropic radial and tangential pressures. Numerical mass-radius diagrams indicate that elasticity can increase maximum mass and compactness by up to approximately 22\%, and stable elastic stars can reach compactness $a_i$9 while satisfying energy conditions and subluminal propagation constraints (Alho et al., 2021).

Taken together, these usages indicate that “elastic” functions both as a literal descriptor of recoverable deformation and as a transferable formal concept: it names energy-minimizing geometries, tunable regularization paths, throughput-preserving robustness measures, and even viable ultracompact stars. This suggests that elasticity has become less a single theory than a family of structurally related frameworks centered on constrained response, stability, and recoverable change.