---
title: 'Excess Solid Stress: Theory and Applications'
url: https://www.emergentmind.com/topics/excess-solid-stress
type: topic
---

# Excess Solid Stress: Theory and Applications

Excess solid stress denotes a stress-like excess quantity localized at a solid interface or defect after subtraction of bulk contributions, and it governs the in-plane forces transmitted by surfaces, solid–liquid interfaces, coherent solid–solid interfaces, and atomic steps. In soft matter, it directly controls capillarity-induced elastic deformation and contact-line force balance; in interfacial thermodynamics it appears through Shuttleworth-type relations linking interfacial free energy to strain; and in a distinct recent biomedical usage it denotes an imaging-derived stress proxy in glioma defined from strain and stiffness contrast rather than an interfacial excess in the Gibbs sense [1310.3941] [1702.00684] [1304.0144] [2510.00009].

## 1. Definitions and thermodynamic basis

For soft solid–liquid interfaces, the Gibbs excess framework defines surface stress as the integral of the interfacial excess traction relative to the bulk value. For a planar interface with normal along \(z\),
\[
\Upsilon_{ij}=\int_{-\infty}^{\infty}\left[\sigma_{ij}(z)-\sigma_{ij}^{\text{bulk}}\right]\,dz,
\]
where \(\sigma_{ij}(z)\) is the local stress and \(\sigma_{ij}^{\text{bulk}}\) is the stress far from the interface. In this form, \(\Upsilon_{ij}\) is the integrated lateral force density that the interfacial zone exerts on the adjoining bulk, and it therefore directly enters the mechanics of soft wetting and elastocapillary deformation [1310.3941].

A central distinction between solids and liquids is that solid surface stress need not equal surface free energy. For solids, the Shuttleworth relation states
\[
\tau_{ij} = \gamma \delta_{ij} + \frac{\partial \gamma}{\partial \varepsilon^s_{ij}},
\]
and for isotropic surfaces,
\[
\Upsilon = \gamma + \frac{d\gamma}{d\varepsilon}.
\]
The excess contribution is the strain-dependent term \(\partial \gamma/\partial \varepsilon^s_{ij}\), or \(d\gamma/d\varepsilon\) in the isotropic scalar form. This term vanishes for simple liquids, for which \(\Upsilon_{ij}=\gamma\delta_{ij}\), but not for elastic solids because the surface can store and release energy under in-plane strain [1702.00684].

The same thermodynamic structure reappears for coherent solid–solid interfaces. Frolov and Mishin define the interface stress tensor as an excess over bulk stresses,
\[
\tau_{ij}^{XY} := \frac{1}{A}[V' \bar Q_{ij}]_{XY},
\]
with the specific excess depending on the chosen reference extensive properties \(X\) and \(Y\). Their generalized Shuttleworth-type equation is
\[
\left(\frac{\partial \gamma}{\partial e_{ij}}\right)^{XY} = \tau_{ij}^{XY} - \delta_{ij}\gamma.
\]
In this setting, excess solid stress is explicitly a thermodynamic excess conjugate to in-plane elastic strain at the interface rather than a simple local stress concentration [1304.0144].

## 2. Soft solid–liquid interfaces and capillarity

The microscopic model developed for amorphous soft solids combines enthalpic interactions and entropic elasticity in a density functional theory with a sharp-interface approximation. Molecular interactions are split into short-range repulsion and long-range attraction: the repulsion contributes an isotropic contact pressure in liquids or an elastic stress tensor in solids, while attraction is represented by effective potentials \(\phi_{\alpha\beta}\) and their planar integrals \(\Pi_{\alpha\beta}(h)\). These quantities recover the classical adhesion and surface energies through
\[
\mathcal A_{AB}=-\int_0^\infty \Pi_{AB}(z)\,dz=\gamma_A+\gamma_B-\gamma_{AB}.
\]
Mechanical equilibrium is written microscopically as
\[
\nabla\cdot\bar\sigma-\nabla\phi=0,
\]
and for amorphous soft solids the constitutive law is small-strain isotropic Hooke’s law [1310.3941].

Within that framework, liquid–vapor and liquid–liquid interfaces satisfy \(\Upsilon=\gamma\), but a solid–liquid interface does not in general. Under plane strain parallel to the interface, the in-plane elastic stress is tied to the normal elastic stress through
\[
\sigma_{xx}^{\text{el}}=\frac{\nu}{1-\nu}\,\sigma_{zz}^{\text{el}},
\]
so the local Poisson ratio \(\nu\) of the interfacial solid layer controls the difference between stress and energy. The resulting closed-form expression for a solid–liquid interface is
\[
\Upsilon_{SL}=\frac{\nu}{1-\nu}\,\gamma_{SL}+\frac{1-2\nu}{1-\nu}\left(\gamma_{SV}+\gamma_{LV}\right),
\]
with excess
\[
\Delta\Upsilon\equiv \Upsilon_{SL}-\gamma_{SL}=\frac{1-2\nu}{1-\nu}\left(\gamma_{SV}+\gamma_{LV}-\gamma_{SL}\right).
\]
The bracketed term is the work of adhesion between the solid and the liquid. Hence \(\Upsilon_{SL}-\gamma_{SL}\) is controlled by interfacial compressibility and vanishes in the incompressible limit \(\nu=1/2\) [1310.3941].

This has immediate consequences for wetting. The elastocapillary length,
\[
\ell_{ec}=\frac{\gamma}{E},
\]
sets the scale of capillarity-induced deformation, but when \(\Upsilon\neq\gamma\) it is the surface stress \(\Upsilon\) that enters the traction boundary condition. At a three-phase contact line on a soft substrate, the transmitted tangential elastic force per unit length becomes
\[
f_t^{\text{el}}=\left(\Upsilon_{SL}-\Upsilon_{SV}\right)-\left(\gamma_{SL}-\gamma_{SV}\right)=\frac{1-2\nu}{1-\nu}\,\gamma\,\left(1+\cos\theta\right).
\]
For \(\nu=1/2\), \(f_t^{\text{el}}=0\), recovering a liquid-like Neumann construction; for \(\nu<1/2\), a finite tangential traction appears and modifies wetting ridge geometry [1310.3941].

The quantitative comparison reported in that work links theory to both simulation and experiment. For a lattice solid with springs yielding \(\nu=0.2\) and \(\gamma_{SL}=\gamma_{SV}\), the predicted \(f_t^{\text{el}}/\gamma=0.75\) agrees with measured values \(0.81\pm0.17\). Experiments on soft elastomeric wires likewise showed significant tangential strain near contact lines, consistent with a compressible interfacial zone even when the bulk material behaves nearly incompressibly [1310.3941].

## 3. Direct measurement and strain-dependent surface stress in soft gels

Direct measurement of strain-dependent solid surface stress in silicone gels established a linear constitutive form over accessible surface strains,
\[
\Upsilon(\varepsilon)=\Upsilon_0+\Lambda\varepsilon,
\]
with \(\Upsilon_0 = 19 \pm 3\) mN/m and \(\Lambda = 126 \pm 6\) mN/m, valid up to \(\varepsilon \approx 0.25\). Here \(\Lambda\) is the dilatational surface elastic modulus, and the measured ratio \(\Lambda/\Upsilon_0 \approx 6.6\) shows that the excess term can greatly exceed the zero-strain surface stress. The surface stress increases by about \(2.5\times\) at \(25\%\) strain and doubles with only \(\approx 17\%\) strain [1702.00684].

The measurement protocol relied on wetting ridges under large glycerol droplets on biaxially stretched silicone gels. The macroscopic contact angle remained \(\theta=90.8^\circ\) and was unchanged by stretching, indicating \(\gamma_{SV}\approx\gamma_{SL}\) and negligible hysteresis, while the microscopic ridge opening angle \(\alpha\) changed strongly with strain. In the symmetric case the local force balance is
\[
2\Upsilon \cos(\alpha/2)=\Upsilon_{LV},
\]
with \(\Upsilon_{LV}=41\pm1\) mN/m. This relation yields \(\Upsilon\) directly from the ridge geometry, independent of bulk elasticity. At zero applied strain, \(\alpha=91.2^\circ\) and \(\Upsilon\approx29\) mN/m because the local strain at the ridge is already \(\approx 8\%\); at \(\varepsilon_\infty=0.18\), \(\alpha=126.3^\circ\) and \(\Upsilon\approx44\) mN/m [1702.00684].

These experiments also established the mechanical boundary condition for a solid free surface with surface stress,
\[
\sigma \cdot n = \nabla_s \cdot \tau_s,
\]
which reduces to
\[
\sigma \cdot n = \Upsilon \kappa n
\]
for isotropic homogeneous \(\Upsilon\), and to
\[
\sigma \cdot n = \nabla_s \Upsilon + \Upsilon \kappa n
\]
when \(\Upsilon\) varies spatially. The measured increase in \(\Upsilon\) with strain directly enlarges the elastocapillary length \(\ell_{ec}=\Upsilon/E\), extending the capillary-dominated near field around contact lines [1702.00684].

A complementary adhesion study connected this excess stress to contact stiffness. In quasi-static pull-off experiments between untreated rigid silica spheres of radius \(R\in[7.9,32.0]\) \(\mu\)m and compliant PDMS gels with \(E=5.6\) kPa and \(\nu=0.48\), the zero-strain surface stress was \(\Upsilon_0\approx0.02\) N/m and the directly measured surface modulus on similar gels satisfied \(\Lambda\approx6\Upsilon_0\). The zero-strain elastocapillary length was
\[
\ell_{ec0}=\Upsilon_0/E \approx 3.6\ \mu\text{m}.
\]
Near contact, the geometry exhibited capillary-like features consistent with this scale [1707.03089].

In that contact problem, the measured adhesive force was represented as
\[
F_{\text{total}} \approx F_{EL}+F_{CL},
\]
where the circumferential surface-stress contribution is
\[
F_{CL}=2\pi a \sin(\Theta)\Upsilon,
\]
and under total wetting,
\[
F_{CL}=2\pi(a^2/R)\Upsilon.
\]
With fixed \(\Upsilon=\Upsilon_0\), the capillary term increases the initial force but adds no stiffness. The missing stiffness arises from the excess, strain-dependent part of the surface stress:
\[
K_{CL}\approx 2\pi(a_0^2/R)\,\frac{d\Upsilon}{dD}.
\]
Using the experimentally observed exponential growth of the constant-curvature domain \(\Delta S\), the study approximated \(\Upsilon(D)\approx \Upsilon_0(\Delta S/\Delta S_0)\), leading to a nonzero \(K_{CL}\) that was required to match the measured initial stiffness \(k_0\) [1707.03089].

The two gel studies together established a common point. Macroscopic observables such as Young–Dupré contact angles report surface energies, but microscopic ridge angles, near-contact curvature, and contact stiffness are sensitive to the strain-dependent excess component of solid surface stress. Neglecting \(d\gamma/d\varepsilon\) therefore removes a mechanically important restoring force [1702.00684] [1707.03089].

## 4. Coherent solid–solid interfaces under nonhydrostatic stress

For coherent solid–solid interfaces, excess solid stress is formulated within a fully thermodynamic theory for multicomponent solids under nonhydrostatic mechanical stresses. The system comprises two homogeneous phases, \(\alpha\) and \(\beta\), separated by a plane coherent interface with unit normal along \(x_3\). Coherency means that the same lattice penetrates both phases, lattice sites are conserved during interface migration, and sliding is prohibited. Because coherent interfaces can sustain shear stresses parallel to the interface, the thermodynamics contains not only interface free energy and interface stress but also interface excess shear [1304.0144].

All interface excess quantities are defined through Cahn’s generalized excess method rather than a geometric dividing surface. For any extensive property \(Z\),
\[
[Z]_{XY}:=
\frac{
\det
\begin{pmatrix}
Z & X & Y\\
Z^\alpha & X^\alpha & Y^\alpha\\
Z^\beta & X^\beta & Y^\beta
\end{pmatrix}
}{
\det
\begin{pmatrix}
X^\alpha & Y^\alpha\\
X^\beta & Y^\beta
\end{pmatrix}
},
\]
which ensures \([X]_{XY}=[Y]_{XY}=0\). The interface free energy is then written as an excess quantity, and the full adsorption equation contains thermal, chemical, normal-stress, shear-stress, and in-plane strain terms [1304.0144].

In small-strain form, the adsorption equation is
\[
d\gamma = -\frac{[S]_{XY}}{A}\,dT - \sum_{k\ge2}\frac{[N_k]_{XY}}{A}\,dM_{k1} - \frac{[N]_{XY}}{A}\,d\phi_1 - \sum_l\frac{[n_l]_{XY}}{A}\,d\mu_l
-\sum_{i=1,2,3}\frac{[V\bar F_{i3}/\bar F_{33}]_{XY}}{A}\,d\sigma_{3i}
+\sum_{i,j=1,2}\left(\tau_{ij}^{XY}-\delta_{ij}\gamma\right)\,de_{ji}.
\]
This identifies the interface stress tensor \(\tau_{ij}^{XY}\) as the quantity thermodynamically conjugate to in-plane elastic strain, while the excess shear measures \([V\bar F_{13}/\bar F_{33}]_{XY}\) and \([V\bar F_{23}/\bar F_{33}]_{XY}\) are conjugate to \(\sigma_{31}\) and \(\sigma_{32}\) [1304.0144].

A distinctive feature of this formalism is nonuniqueness. Because the generalized excess \([\,\cdot\,]_{XY}\) depends on the chosen reference pair \(X,Y\), the numerical value of \(\tau_{ij}^{XY}\) also depends on that choice. Frolov and Mishin interpret this as path dependence on the coherent phase coexistence hypersurface rather than as a failure of definition. Under hydrostatic stress, the interface stress reduces to the unique form
\[
\tau_{ij}=\frac{V}{A}\left(\bar\sigma_{ij}+p\delta_{ij}\right), \qquad i,j=1,2.
\]
The theory also yields Maxwell relations that couple interface stress to segregation, excess volume, excess shear, temperature, and nonhydrostatic loading [1304.0144].

This coherent-interface formulation clarifies a common misconception inherited from liquid interfaces. For coherent solids, interface stress is not simply a scalar tension attached to area creation; it is a tensorial excess whose value depends on strain state, composition variables, and the constrained thermodynamic path used to define the excess [1304.0144].

## 5. Step excess stress at faceted crystalline surfaces

At atomic steps on faceted crystalline surfaces, excess solid stress appears as the line-defect analog of surface stress. The thermodynamic theory of steps on faceted surfaces defines the step free energy per unit length \(\gamma^{\mathrm{st}}\) through an Euler relation for a step-containing region and expresses it in terms of step excess quantities using Cahn’s determinant formalism. The resulting step adsorption equation is
\[
d(\gamma^{\mathrm{st}}L)
=
-[S]_{XY}\,dT
-[N]_{XY}\,d\mu
-[A]_{XY}\,d\gamma
+\sum_{i,j}^{x,y}[\sigma_{ij}V-\delta_{ij}\gamma A]_{XY}\,d\varepsilon_{ij},
\]
where \(x\) is parallel to the step and \(y\) is perpendicular to it [1701.07413].

The corresponding step excess stress tensor is
\[
[\tau_{ij}]_{XY}
\equiv
\frac{1}{L}\frac{\partial(\gamma^{\mathrm{st}}L)}{\partial\varepsilon_{ij}}
=
\frac{1}{L}[\sigma_{ij}V-\delta_{ij}\gamma A]_{XY},
\qquad i,j=x,y,
\]
with Shuttleworth-like relation
\[
[\tau_{ij}]_{XY}=
\delta_{ix}\delta_{jx}\gamma^{\mathrm{st}}
+\frac{\partial\gamma^{\mathrm{st}}}{\partial\varepsilon_{ij}}.
\]
For the practical dividing-line choice \(X=A\), \(Y=N\), the stress excess simplifies to
\[
[\tau_{ij}]_{AN}=\frac{1}{L}[\sigma_{ij}V]_{AN}.
\]
Thus, as for surfaces and coherent interfaces, the excess stress is distinct from the free energy itself and is intrinsically strain dependent [1701.07413].

Atomistic calculations for \(\langle 110\rangle\) steps on Cu(111) using the Mishin EAM potential show that both excess energy and excess step stress have weak temperature dependence up to a homologous temperature of approximately \(0.6\), above which they increase strongly and the step stress becomes more isotropic. At \(0\) K,
\[
[\tau_\perp]_{AN} = - 38.3 \,\mathrm{meV/}, \qquad
[\tau_\parallel]_{AN} = + 34.3 \,\mathrm{meV/},
\]
so the low-temperature step stress is anisotropic, being compressive perpendicular to the step and tensile parallel to it. The step free energy decreases from \((103.61 \pm 0.02)\) meV/ at \(0\) K to \((45.8 \pm 0.4)\) meV/ at the melting point \(T_m=1327\) K, but remains finite up to \(T_m\), indicating the absence of a roughening temperature for this \(\{111\}\) facet in the model [1701.07413].

This line-defect formalism extends the meaning of excess solid stress from two-dimensional interfaces to one-dimensional defects. It also shows that the excess-stress concept is not restricted to soft matter or continuum interfaces: it remains well defined in atomistic crystalline thermodynamics through appropriate excess constructions [1701.07413].

## 6. Surface-layer mechanics and generalized Young–Laplace equations for solids

A different but related reconstruction of excess solid stress arises when the surface of a solid is modeled as a finite-thickness layer with membrane and transverse shear responses. In Zaixing Huang’s formulation, the natural boundary condition derived from a bulk-plus-surface Lagrangian reduces in quasi-static conditions to
\[
\sigma_{kj}n_j = D_A \sigma^s_{Ak},
\]
or in tensor notation,
\[
\boldsymbol{\sigma}\,\boldsymbol{n} = \nabla_s \cdot \boldsymbol{\sigma}^s.
\]
Here \(\sigma^s_{Ak}=(1-\chi H)\tilde\sigma_{Ak}\) is a curvature-renormalized surface stress, with \(\chi\) the Tolman length and \(H\) the mean curvature [1806.02582].

The surface stress is decomposed into membrane and transverse shear parts,
\[
\sigma^s_{Ak}=\tau^s_{AB}g^B{}_k + q_A n_k,
\]
where \(\tau^s_{AB}\) is the in-plane membrane stress tensor and \(q_A\) is the surface transverse shear stress. Projecting the traction balance gives tangential and normal equilibrium:
\[
P\,\boldsymbol{\sigma}\,\boldsymbol{n}
=
\left\{
D_B[(1-\chi H)\tau^s_{AB}]
-(1-\chi H)b_A{}^B q_B
\right\}g^A,
\]
\[
\boldsymbol{n}\cdot\boldsymbol{\sigma}\,\boldsymbol{n}
=
D_A[(1-\chi H)q_A]
+(1-\chi H)b^{AB}\tau^s_{AB}.
\]
These equations generalize the liquid Young–Laplace equation by including in-plane shearing and transverse shearing within the surface layer [1806.02582].

For the intrinsic constitutive content of the surface layer, the paper adopts a Shuttleworth–Herring-type membrane law,
\[
\tau^s_{AB}=\gamma g_{AB}+\frac{\partial\gamma}{\partial\varepsilon^s_{AB}},
\]
and, for the purely intrinsic part without explicit in-plane strain,
\[
\tau^s_{AB}=\gamma g_{AB}.
\]
The intrinsic transverse shear stress is proposed as
\[
q_A=t\,D_A(K^{-1}), \qquad t:=\mu\varepsilon,
\]
so it depends on the gradient of the inverse Gaussian curvature and a strain-dependent coefficient. This term vanishes on a sphere because \(K\) is constant, and then the normal balance reduces to a Tolman-corrected Laplace form; if also \(\chi=0\), the classical liquid result \(\Delta p=2\gamma H\) is recovered [1806.02582].

In this framework, excess solid stress refers to the residual bulk stress induced by intrinsic surface membrane stress and transverse shear. Even without external loading, nonzero surface tractions generated by \(\nabla_s\cdot\sigma^s\) must be balanced by internal stresses satisfying \(\nabla\cdot\sigma=0\). The formulation therefore shifts emphasis from excess as a purely thermodynamic subtraction to excess as a surface-originated loading that generates residual stress fields in the interior of a solid [1806.02582].

## 7. Imaging-derived excess solid stress in glioma

In the glioma study, “excess solid stress” is defined operationally rather than as an interfacial excess. The deformation field \(u(r)\) maps undeformed coordinates \(r\) to \(r' = r + u(r)\), with deformation gradient
\[
F=I+\nabla u,
\]
and Green–Lagrange strain
\[
E=\frac12(F^T F-I).
\]
Two regional strain measures are used: volumetric strain \(\mathrm{tr}(E)\) and octahedral shear strain
\[
\mathrm{OSS}=\sqrt{\frac{2}{3}\left[(E_{11}-E_{22})^2 + (E_{22}-E_{33})^2 + (E_{33}-E_{11})^2 + 6(E_{12}^2 + E_{23}^2 + E_{13}^2)\right]}.
\]
The shear modulus is derived from magnetic resonance elastography as
\[
\mu=\rho\,\mathrm{SWS}^2,
\]
with tissue density \(\rho=1000\) kg/m\(^3\), and the stiffness differential is
\[
\Delta\mu=\mu_o-\mu_1,
\]
where \(\mu_o\) is the modulus of apparently unaffected brain and \(\mu_1\) that of peritumoral tissue [2510.00009].

Excess solid stress is then defined as the product of strain and stiffness differential:
\[
E_{\mathrm{vol}}=\Delta\mu\cdot\mathrm{tr}(E), \qquad
E_{\mathrm{shear}}=\Delta\mu\cdot\mathrm{OSS}.
\]
These quantities have units of Pa. In the same study, volumetric stress and shear stress are estimated from an isotropic constitutive law,
\[
\sigma=\lambda\,\mathrm{tr}(E)\,I + 2\mu E,
\]
with \(\nu=0.4\) used for volumetric stress computation [2510.00009].

The measured deformation patterns in patients were spatially heterogeneous and extended well beyond visible tumor margins. Non-tumor brain displacement was \(3.7\pm0.5\) mm in glioma patients versus \(3.0\pm0.3\) mm in healthy volunteers. Peritumoral volumetric strain was \(-17\pm9\%\), unaffected-brain volumetric strain in patients was \(-9\pm5\%\), and healthy-volunteer volumetric strain was \(0\pm2\%\). Peritumoral OSS was \(23\pm9\%\), compared with \(13\pm6\%\) in unaffected brain and \(0\pm2\%\) in healthy volunteers [2510.00009].

The principal clinical result concerned survival. Excess volumetric stress had mean \(4.3\pm11.6\) Pa, range \(-13\pm6\) Pa to \(30\pm7\) Pa, and showed a linear correlation with survival days with slope \(-0.017\pm0.006\) Pa/day, \(R=-0.70\), and \(p=0.02\). Excess shear stress had mean \(-7.2\pm14.2\) Pa and weaker prognostic value (\(p=0.13\)). Neither deformation magnitude alone, tumor size, nor any single-region stiffness parameter correlated with survival. In this biomedical usage, excess solid stress is therefore an imaging-derived biomarker that integrates cause (strain) and consequence (stiffness degradation), rather than a surface or interface thermodynamic excess [2510.00009].

Across these literatures, excess solid stress consistently signifies a contribution that cannot be reduced to bulk elasticity alone. In soft interfaces it is the strain-dependent increment beyond surface energy; in coherent interfaces it is a generalized excess conjugate to in-plane strain; at steps it is a line excess linked to step free energy; in surface-layer mechanics it is the residual stress generated by intrinsic surface tractions; and in glioma imaging it is a regional stress proxy built from strain and stiffness contrast. The shared theme is that solids support mechanically consequential excess stresses that liquids, in their simplest description, do not.

Source: https://www.emergentmind.com/topics/excess-solid-stress