Nemato-Elasticity: Coupling Order & Mechanics

• Nemato-elasticity is the coupling between nematic orientational order and mechanical strain, leading to programmed deformations in materials such as elastomers and thin sheets.
• Theoretical frameworks, including metric formulations, bulk constitutive theories, and viscoelastic models, capture diverse phenomena from defect formation to actuation.
• This topic underpins practical advances in responsive materials design, enabling applications in stress mitigation, locomotion, and adaptive electronic systems.

Nemato-elasticity denotes the coupling between nematic orientational order and mechanical response. In the literature, that coupling appears in several distinct but mathematically related forms: as a bilinear interaction between a nematic tensor and elastic strain, as a spontaneous internal metric that evolves with domain structure, as a target metric for thin sheets undergoing nematic–isotropic transition, and as a coupling between electronic nematicity and lattice distortions. Across nematic elastomers and glasses, thin shells and ribbons, hybrid molecular–colloidal and biopolymer nematics, viscoelastic liquid crystals, and electronically driven nematic systems, the common theme is that orientational symmetry breaking modifies admissible deformations, while elastic constraints reshape the order parameter, its fluctuations, and the resulting morphology or dynamics (Pismen, 2014, Feng et al., 2014, Meese et al., 31 Jul 2025).

1. Foundational formulations

A standard continuum representation of in-plane nematic order is the traceless symmetric tensor

$$Q^\alpha{}_\beta = S\left(n^\alpha n_\beta -\frac12 \delta^\alpha_\beta\right),$$

which in thin-sheet theories directly determines the post-transition metric. For a nematic-to-isotropic transition in a thin sheet, the induced metric can be written as

$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$

and in local coordinates parallel and perpendicular to the director the spontaneous stretches reduce to

$$\xi^1=\frac{1}{\ell}x^1,\qquad \xi^2=\ell x^2,\qquad \ell=\sqrt{\frac{1+a}{1-a}}.$$

This is the geometric core of metric theories of thin nemato-elastic bodies: orientational order prescribes a non-Euclidean target metric, and a sufficiently thin body bends rather than storing large in-plane strain (Zakharov et al., 2015).

In bulk polydomain nematic elastomers, the same coupling is often written through an internal spontaneous metric rather than a target surface metric. A representative form is

$$\mathbf G = \mathbf P\, \mathrm{diag}\!\left(\Lambda^2,\frac{\Delta^2}{\Lambda^2},\frac{1}{\Delta^2}\right) \mathbf P^T,$$

with coarse-grained free energy

$$W = W_e + W_r,\qquad W_e = \frac{1}{2}\mu\left[\mathrm{tr}\!\left(\mathbf F^T\mathbf G^{-1}\mathbf F\right)-3\right].$$

Here the elastic response is measured relative to an evolving internal metric $\mathbf G$, so deformation can be accommodated not only by $\mathbf F$ but also by internal domain evolution parametrized by $\Lambda$ and $\Delta$ (Maghsoodi et al., 2023).

A distinct but closely related formulation appears in biopolymer gels, where nematic order is not embedded in the reference state but induced by deformation. The Landau-type free energy for shear deformations is

$$F_s = \int \left[ \bar{\mu}\,\operatorname{Tr}\tilde{\mathbf v}^2 - t\,\operatorname{Tr}(\tilde{\mathbf v}\cdot \mathbf Q) + V(\mathbf Q) \right]\,d^dr,$$

with $$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$0 the traceless strain and $$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$1 the induced alignment tensor. In the linear regime,

$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$2

This suggests a general distinction between embedded nematicity and strain-induced nematicity: in one case elasticity is measured relative to a nematic internal state, and in the other the elastic field itself generates nematic order (Feng et al., 2014).

2. Metric programming in thin sheets, shells, and ribbons

Thin-sheet nemato-elasticity is most explicit in metric theories of nematic elastomer and glass films undergoing nematic–isotropic transition. In the Lagrangian description of a very thin, incompressible sheet with frozen director field, an infinitesimal interval transforms as

$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$3

and the pre- and post-transition metrics satisfy

$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$4

For a charge-$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$5 vortex texture, the induced metric is realized by a cone, with

$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$6

For the opposite $$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$7 aster texture, the circumference-to-radius ratio is too large for a smooth axisymmetric cone, and the sheet forms an anticone assembled from monotone patches separated by radial creases. The creases occur at equally spaced azimuths

$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$8

so the number of creases is necessarily even. The minimal crease count obeys

$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$9

and increases with spontaneous strain. Radially varying $$\xi^1=\frac{1}{\ell}x^1,\qquad \xi^2=\ell x^2,\qquad \ell=\sqrt{\frac{1+a}{1-a}}.$$0 or $$\xi^1=\frac{1}{\ell}x^1,\qquad \xi^2=\ell x^2,\qquad \ell=\sqrt{\frac{1+a}{1-a}}.$$1, interpreted as doping, produces curved generatrices and distributed Gaussian curvature rather than curvature concentrated only at the apex; near the origin, if $$\xi^1=\frac{1}{\ell}x^1,\qquad \xi^2=\ell x^2,\qquad \ell=\sqrt{\frac{1+a}{1-a}}.$$2, then

$$\xi^1=\frac{1}{\ell}x^1,\qquad \xi^2=\ell x^2,\qquad \ell=\sqrt{\frac{1+a}{1-a}}.$$3

The apex singularity is regularized by defect-core melting of nematic order, modeled through a Landau–de Gennes amplitude $$\xi^1=\frac{1}{\ell}x^1,\qquad \xi^2=\ell x^2,\qquad \ell=\sqrt{\frac{1+a}{1-a}}.$$4 with $$\xi^1=\frac{1}{\ell}x^1,\qquad \xi^2=\ell x^2,\qquad \ell=\sqrt{\frac{1+a}{1-a}}.$$5, so that the anisotropy vanishes quadratically near the core (Pismen, 2014).

Half-charged defects generate a different geometry. For isolated $$\xi^1=\frac{1}{\ell}x^1,\qquad \xi^2=\ell x^2,\qquad \ell=\sqrt{\frac{1+a}{1-a}}.$$6 defects in a thin sheet, the transformed metric yields Gaussian curvature

$$\xi^1=\frac{1}{\ell}x^1,\qquad \xi^2=\ell x^2,\qquad \ell=\sqrt{\frac{1+a}{1-a}}.$$7

The $$\xi^1=\frac{1}{\ell}x^1,\qquad \xi^2=\ell x^2,\qquad \ell=\sqrt{\frac{1+a}{1-a}}.$$8 defect therefore carries a dipolar curvature singularity, while the $$\xi^1=\frac{1}{\ell}x^1,\qquad \xi^2=\ell x^2,\qquad \ell=\sqrt{\frac{1+a}{1-a}}.$$9 defect carries a hexapolar one. In elliptic domains, the two foci act as $$\mathbf G = \mathbf P\, \mathrm{diag}\!\left(\Lambda^2,\frac{\Delta^2}{\Lambda^2},\frac{1}{\Delta^2}\right) \mathbf P^T,$$0 defects, and the exact curvature field is

$$\mathbf G = \mathbf P\, \mathrm{diag}\!\left(\Lambda^2,\frac{\Delta^2}{\Lambda^2},\frac{1}{\Delta^2}\right) \mathbf P^T,$$1

Boundary anchoring selects the texture and flips the curvature map, while multiply connected domains naturally generate $$\mathbf G = \mathbf P\, \mathrm{diag}\!\left(\Lambda^2,\frac{\Delta^2}{\Lambda^2},\frac{1}{\Delta^2}\right) \mathbf P^T,$$2 defects. A triangulated discrete relaxation scheme then produces three-dimensional embeddings by minimizing an energy with stretching and bending penalties,

$$\mathbf G = \mathbf P\, \mathrm{diag}\!\left(\Lambda^2,\frac{\Delta^2}{\Lambda^2},\frac{1}{\Delta^2}\right) \mathbf P^T,$$3

with $$\mathbf G = \mathbf P\, \mathrm{diag}\!\left(\Lambda^2,\frac{\Delta^2}{\Lambda^2},\frac{1}{\Delta^2}\right) \mathbf P^T,$$4 taken small (Zakharov et al., 2015).

Ribbons provide a complementary thin-body problem. For twisted nematic elastomer strips, a Koiter-type plate model with incompatible reference forms

$$\mathbf G = \mathbf P\, \mathrm{diag}\!\left(\Lambda^2,\frac{\Delta^2}{\Lambda^2},\frac{1}{\Delta^2}\right) \mathbf P^T,$$5

captures the helicoid-to-spiral transition. The helicoid carries stretching but little bending mismatch,

$$\mathbf G = \mathbf P\, \mathrm{diag}\!\left(\Lambda^2,\frac{\Delta^2}{\Lambda^2},\frac{1}{\Delta^2}\right) \mathbf P^T,$$6

whereas the spiral is an isometry that pays bending,

$$\mathbf G = \mathbf P\, \mathrm{diag}\!\left(\Lambda^2,\frac{\Delta^2}{\Lambda^2},\frac{1}{\Delta^2}\right) \mathbf P^T,$$7

Comparison yields the threshold

$$\mathbf G = \mathbf P\, \mathrm{diag}\!\left(\Lambda^2,\frac{\Delta^2}{\Lambda^2},\frac{1}{\Delta^2}\right) \mathbf P^T,$$8

for the helicoid to be favored. The transition is therefore governed by the competition between stretching and bending in a thin incompatible plate (Tomassetti et al., 2016).

3. Bulk constitutive and variational theories

For isotropic-genesis polydomain nematic elastomers, the constitutive state space is explicitly enlarged by internal variables. The admissible states satisfy

$$\mathbf G = \mathbf P\, \mathrm{diag}\!\left(\Lambda^2,\frac{\Delta^2}{\Lambda^2},\frac{1}{\Delta^2}\right) \mathbf P^T,$$9

with limiting cases

$$W = W_e + W_r,\qquad W_e = \frac{1}{2}\mu\left[\mathrm{tr}\!\left(\mathbf F^T\mathbf G^{-1}\mathbf F\right)-3\right].$$0

$$W = W_e + W_r,\qquad W_e = \frac{1}{2}\mu\left[\mathrm{tr}\!\left(\mathbf F^T\mathbf G^{-1}\mathbf F\right)-3\right].$$1

$$W = W_e + W_r,\qquad W_e = \frac{1}{2}\mu\left[\mathrm{tr}\!\left(\mathbf F^T\mathbf G^{-1}\mathbf F\right)-3\right].$$2

The free energy

$$W = W_e + W_r,\qquad W_e = \frac{1}{2}\mu\left[\mathrm{tr}\!\left(\mathbf F^T\mathbf G^{-1}\mathbf F\right)-3\right].$$3

together with overdamped internal-variable kinetics

$$W = W_e + W_r,\qquad W_e = \frac{1}{2}\mu\left[\mathrm{tr}\!\left(\mathbf F^T\mathbf G^{-1}\mathbf F\right)-3\right].$$4

produces the semi-soft sequence of an initial elastic regime, a softened domain-evolution window, and hardening near saturation. In this framework, nemato-elasticity is a genuine internal-variable theory rather than a fixed anisotropic constitutive law (Maghsoodi et al., 2023).

A solvable mean-field counterpart is provided by the elastic Maier–Saupe–Zwanzig model, which combines a lattice Maier–Saupe interaction with Warner–Terentjev elasticity and quadrupolar quenched random fields. In the uniaxial specialization,

$$W = W_e + W_r,\qquad W_e = \frac{1}{2}\mu\left[\mathrm{tr}\!\left(\mathbf F^T\mathbf G^{-1}\mathbf F\right)-3\right].$$5

The model predicts a stress–strain coexistence curve below a freezing temperature, analogous to the $$W = W_e + W_r,\qquad W_e = \frac{1}{2}\mu\left[\mathrm{tr}\!\left(\mathbf F^T\mathbf G^{-1}\mathbf F\right)-3\right].$$6–$$W = W_e + W_r,\qquad W_e = \frac{1}{2}\mu\left[\mathrm{tr}\!\left(\mathbf F^T\mathbf G^{-1}\mathbf F\right)-3\right].$$7 diagram of a simple fluid, and below a critical disorder strength the tie lines resemble the experimental stress–strain plateau associated with the polydomain–monodomain transition. In the monodomain case, random fields may soften the first-order isotropic–nematic transition, provided samples are formed in the nematic state (Liarte et al., 2011).

A rigorous variational formulation treats the director as a field on the deformed configuration rather than the reference body. The energies are

$$W = W_e + W_r,\qquad W_e = \frac{1}{2}\mu\left[\mathrm{tr}\!\left(\mathbf F^T\mathbf G^{-1}\mathbf F\right)-3\right].$$8

Under the determinant constraint $$W = W_e + W_r,\qquad W_e = \frac{1}{2}\mu\left[\mathrm{tr}\!\left(\mathbf F^T\mathbf G^{-1}\mathbf F\right)-3\right].$$9, with $\mathbf G$0 for $\mathbf G$1 and no cavitation, the lower semicontinuous envelope is

$\mathbf G$2

Here the elastic term is replaced by its quasiconvexification and the nematic term by its tangential quasiconvexification. This provides a mathematically precise effective theory for nemato-elastic systems that develop fine-scale elastic or orientational microstructure (Mora-Corral et al., 2017).

4. Nonuniform loading, contact mechanics, and actuation

The consequences of nemato-elasticity are especially sharp under nonuniform loading. In the classical Hertz problem for a rigid spherical indenter of radius $\mathbf G$3 on an isotropic elastic half-space,

$\mathbf G$4

Polydomain nematic elastomers instead display a universal three-regime response: $\mathbf G$5 For the fitted material $\mathbf G$6, simulation gives $\mathbf G$7, while experiment gives $\mathbf G$8 to $\mathbf G$9 for the typical 10% crosslinked sample; a weaker-crosslinked sample gives $\mathbf F$0. The softened intermediate exponent is not a simple consequence of a softer scalar modulus, but of a transformed zone beneath the indenter in which $\mathbf F$1 evolves strongly toward a biaxial pattern and the contact pressure becomes lower and flatter. Above the nematic–isotropic transition temperature the same materials revert to Hertzian behavior, identifying domain evolution as the controlling mechanism (Maghsoodi et al., 2023).

A related redistribution of stress occurs in adhesion. For a conventional adhered linear elastic cylinder, the interfacial normal stress obeys

$\mathbf F$2

so the edge singularity controls failure. In an isotropic-genesis polydomain nematic elastomer cylinder under tensile load, the same singularity is recovered only at low load. As the soft regime activates, the edge singularity disappears, the edge stress drops dramatically, and additional load is transferred to the interior. For $\mathbf F$3, the shift of the maximum stress from edge to center occurs at

$\mathbf F$4

At high load near the edge,

$\mathbf F$5

and

$\mathbf F$6

indicating nearly uniaxial local deformation approaching a monodomain state. Soft elasticity thus acts as a stress-concentration mitigation mechanism, not merely as a reduction of effective modulus (Maghsoodi et al., 2024).

Nemato-elastic actuation can also generate locomotion. In overdamped Lagrangian finite-element models of slender rods and stripes, the scaled elastic energy density is

$\mathbf F$7

A traveling beam that switches a local region between states of different nematic order advects a window of altered intrinsic strain $\mathbf F$8 and, for cross-sectional director gradients, altered spontaneous curvature $\mathbf F$9. For an inflexible rod on a high-friction substrate, the net displacement per beam passage is

$\Lambda$0

against the beam propagation direction, giving speed

$\Lambda$1

With splayed director through thickness, the curvature radius obeys

$\Lambda$2

and the motion direction can reverse. The same framework predicts buckling-mediated slowdown, steering under oblique actuation, and migration on substrate-friction gradients (Zakharov et al., 2015).

5. Disorder, defects, and hybrid nematic materials

In isotropic-genesis nematic elastomers, the network acts not as a rigid quenched field but as a thermally compliant medium that both remembers the orientational pattern present at cross-linking and generates its own random anisotropy. The composite random field is

$\Lambda$3

and the quadratic free energy becomes

$\Lambda$4

Because $\Lambda$5 is strongest at small $\Lambda$6, the correlators can become oscillatory when

$\Lambda$7

The thermal correlator then has oscillation wavelength

$\Lambda$8

This places memory, localization length, and random anisotropy at the center of nemato-elastic order in disordered networks (Lu et al., 2012).

In quasi-two-dimensional actin nematics at an oil–water interface, defect morphology directly reports elastic anisotropy. The Frank free energy reduces to splay and bend terms,

$\Lambda$9

with

$\Delta$0

For $\Delta$1 defects, $\Delta$2 gives a U-shaped texture, whereas $\Delta$3 gives a V-shaped one. Experimentally, $\Delta$4 rises from about $\Delta$5 to $\Delta$6 as actin length doubles from $\Delta$7 to $\Delta$8, while longer filaments eventually reverse that trend and for $\Delta$9 actin the paper reports

$$F_s = \int \left[ \bar{\mu}\,\operatorname{Tr}\tilde{\mathbf v}^2 - t\,\operatorname{Tr}(\tilde{\mathbf v}\cdot \mathbf Q) + V(\mathbf Q) \right]\,d^dr,$$0

Sparse aligned microtubules act as elastic dopants, producing

$$F_s = \int \left[ \bar{\mu}\,\operatorname{Tr}\tilde{\mathbf v}^2 - t\,\operatorname{Tr}(\tilde{\mathbf v}\cdot \mathbf Q) + V(\mathbf Q) \right]\,d^dr,$$1

and allowing optical extraction of absolute moduli, including

$$F_s = \int \left[ \bar{\mu}\,\operatorname{Tr}\tilde{\mathbf v}^2 - t\,\operatorname{Tr}(\tilde{\mathbf v}\cdot \mathbf Q) + V(\mathbf Q) \right]\,d^dr,$$2

and

$$F_s = \int \left[ \bar{\mu}\,\operatorname{Tr}\tilde{\mathbf v}^2 - t\,\operatorname{Tr}(\tilde{\mathbf v}\cdot \mathbf Q) + V(\mathbf Q) \right]\,d^dr,$$3

The same continuum framework captures defect annihilation and long-time arrested structures (Zhang et al., 2017).

Hybrid molecular–colloidal nematics and composite colloidal nematics generalize this idea. For anti-nematic rods in a thermotropic nematic solvent, homeotropic anchoring produces a marked increase of the splay elasticity, while for co-aligned rods the dominant enhancement shifts to bend. In the hybrid 5CB–rod system, the measured $$F_s = \int \left[ \bar{\mu}\,\operatorname{Tr}\tilde{\mathbf v}^2 - t\,\operatorname{Tr}(\tilde{\mathbf v}\cdot \mathbf Q) + V(\mathbf Q) \right]\,d^dr,$$4 increases linearly with rod fraction, and at $$F_s = \int \left[ \bar{\mu}\,\operatorname{Tr}\tilde{\mathbf v}^2 - t\,\operatorname{Tr}(\tilde{\mathbf v}\cdot \mathbf Q) + V(\mathbf Q) \right]\,d^dr,$$5 the theory gives

$$F_s = \int \left[ \bar{\mu}\,\operatorname{Tr}\tilde{\mathbf v}^2 - t\,\operatorname{Tr}(\tilde{\mathbf v}\cdot \mathbf Q) + V(\mathbf Q) \right]\,d^dr,$$6

so the increase is overwhelmingly surface-anchoring-driven (Senyuk et al., 2021). In Onsager–Straley theory for composite colloidal nematics, anti-nematic rods obey

$$F_s = \int \left[ \bar{\mu}\,\operatorname{Tr}\tilde{\mathbf v}^2 - t\,\operatorname{Tr}(\tilde{\mathbf v}\cdot \mathbf Q) + V(\mathbf Q) \right]\,d^dr,$$7

whereas anti-nematic disks satisfy

$$F_s = \int \left[ \bar{\mu}\,\operatorname{Tr}\tilde{\mathbf v}^2 - t\,\operatorname{Tr}(\tilde{\mathbf v}\cdot \mathbf Q) + V(\mathbf Q) \right]\,d^dr,$$8

The logarithmic dependence on anti-nematic order and the large bend–splay ratio for anti-nematic discotic nematics show that nemato-elasticity in mixed or hybrid systems is controlled not only by concentration and aspect ratio but also by whether the embedded component is nematic or anti-nematic (Wensink, 2018).

6. Compatibility, viscoelasticity, and critical dynamics

In nematic liquids, a viscoelastic extension of hydrodynamics introduces an evolving natural configuration through the multiplicative split

$$F_s = \int \left[ \bar{\mu}\,\operatorname{Tr}\tilde{\mathbf v}^2 - t\,\operatorname{Tr}(\tilde{\mathbf v}\cdot \mathbf Q) + V(\mathbf Q) \right]\,d^dr,$$9

with strain energy

$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$00

and shape tensor

$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$01

The co-deformational derivative

$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$02

then evolves through

$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$03

In the low-frequency limit this reduces to compressible Ericksen–Leslie theory and yields Parodi’s relation

$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$04

together with

$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$05

A symmetry analysis of the relaxation tensor shows that uniaxial symmetry allows exactly four relaxation times (Turzi, 2016).

In electronically driven nematicity, the strain field is itself constrained. The defining bilinear coupling is

$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$06

but the elastic strain must satisfy the Saint Venant compatibility condition

$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$07

Using a co-rotating helical basis, the compatible strain occupies only three helical directions, while the five-component nematic field decomposes as

$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$08

After integrating out elastic degrees of freedom, the effective quadratic theory suppresses $$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$09, $$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$10, and $$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$11, leaving the helical doublet $$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$12 as the soft sector. In that sense elasticity bestows tensor compatibility upon the nematic order parameter by suppressing incompatible nematic fluctuations. In defect-free media this produces direction-selective criticality even without crystalline anisotropy; in the presence of defects, elastic pinning fields generate random longitudinal and transverse conjugate fields for the local nematic order parameter (Meese et al., 31 Jul 2025).

A microscopic metallic theory reaches a complementary conclusion from the opposite direction. Starting from a $$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$13 Pomeranchuk field $$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$14 and an impurity-assisted electron–transverse-acoustic-phonon coupling, integrating out electrons generates the mixed bubble

$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$15

The nematic quantum critical point is shifted to

$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$16

and the collective excitations become two hybridized nematic–phonon modes. Near criticality the nematic-like branch remains linear in $$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$17 and damped, whereas along the diagonal at $$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$18 the phonon-like branch obeys

$$g_{\alpha\beta}=(1-a^2)^{-1}\left[(1+a^2)g^0_{\alpha\beta} - 2a Q_{\alpha\beta}\right],$$19

The critical soft mode is therefore transferred to the acoustic sector rather than remaining a pure nematic fluctuation (Christensen et al., 4 Jul 2025).

Across these formulations, nemato-elasticity is not a single constitutive law but a family of coupled theories in which orientational order changes the metric, elastic moduli, or low-energy mode content of matter. Thin-sheet theories emphasize programmed intrinsic geometry; bulk elastomer theories emphasize evolving spontaneous metrics and semi-soft response; hybrid and disordered systems emphasize defect morphology, random anisotropy, and elastic hierarchy; and electronic and viscoelastic theories emphasize compatibility, relaxation spectra, and hybrid collective modes. Taken together, these works place nemato-elasticity at the intersection of symmetry, geometry, and mechanics.