Wavevector-Dependent Equilibrium Transverse Modulus

• The topic defines M_T(k) as the coefficient linking Fourier shear stress and strain, revealing scale-dependent and nonlocal elastic behavior.
• Effective medium theories and dynamic homogenization models extract M_T(k) to analyze shear responses in disordered, composite, and network materials.
• Microscopic and statistical analyses capture the transition from macroscopic shear modulus to k-dependent regimes, informing design in soft matter and amorphous systems.

The wavevector-dependent equilibrium transverse modulus, variably denoted as $M_T(k)$, $G_T(k)$, or $\mu(k)$, quantifies the shear elasticity of a material as a function of deformation wavelength or reciprocal length $k$. In contrast to the conventional (macroscopic) shear modulus, which applies in the limit $k \to 0$, the wavevector dependence encodes the nonlocal or scale-dependent character of shear response, including dispersion, heterogeneity, and structural inhomogeneities. The $k$-dependent transverse modulus is central to understanding both dynamic (propagating or diffusive) and static (equilibrium) responses in ordered and disordered solids, fluids, and network materials. It arises naturally in effective medium theories, dynamic homogenization frameworks, Green–Kubo analyses, and field-theoretic elasticity of mesoscale and amorphous systems.

1. Fundamental Definitions and Mathematical Formulations

The equilibrium wavevector-dependent transverse modulus $M_T(k)$ (or $G_T(k)$, $\mu(k)$ depending on context) is defined as the coefficient in the stress–strain response to a transverse (shear) deformation at spatial frequency $k$. In the most general terms, if a plane-wave shear displacement field $u_y(x) = u_0 \cos(kx)$ is imposed, the resulting stress satisfies

$\sigma_{xy}(k) = M_T(k)\, \varepsilon_{xy}(k),$

where $\varepsilon_{xy}(k)$ is the corresponding Fourier component of the strain field. For isotropic media, only two independent elastic moduli (shear and bulk) are needed to describe linear response at each $k$.

The transverse modulus can equivalently be extracted from correlation functions:

• In molecular systems, via the static transverse projection of the stress correlation tensor $$C_{xy,xy}(k)=\langle\hat{\sigma}_{xy}(k)\hat{\sigma}_{xy}(-k)\rangle$$ (Uneyama, 2022, Wittmer et al., 2023).
• In network and periodic materials, via effective-medium or homogenization theories yielding $M_T(k)$ from the energy density or cell-averaged stress under imposed deformations (Katz et al., 2012, Norris et al., 2012).

2. Effective Medium Theory and Network Models

For random rod networks, the effective-medium theory yields a closed-form expression for the wavevector-dependent equilibrium transverse modulus. In a system of randomly distributed nodes connected by straight, slender, linearly elastic rods of Young's modulus $E$ and cross-section $A$, with isotropic rod orientation and length distribution $P(\ell)$, the modulus is (Katz et al., 2012):

$$M_T(k) = \frac{n C A E}{k^2} \int_0^\infty \frac{P(\ell)}{\ell}\, d\ell \int_{-1}^{1} (1-\mu^2)\, \sin^2\biggl(\frac{k\ell\mu}{2}\biggr)\, d\mu, \quad \mu \equiv \cos\theta$$

Here, $n$ is node density, $C$ is mean coordination number, and the brackets denote averages over orientation and phase. The modulus is nonlocal and exhibits characteristic crossover:

• Long-wavelength limit ($k \to 0$):

$M_T(0) = \mathcal{F} E / 15$

where $\mathcal{F}$ is the volume fraction. This recovers the classical continuum shear modulus.

• Short-wavelength regime ($k\ell \gg 1$):

$M_T(k) \propto k^{-2}$

as the rods' contribution decouples at shorter length scales.

This framework captures the reduction of shear stiffness at finite $k$ due to the phase mismatch among rods of different orientations and lengths. The wavevector dependence is strictly a manifestation of the nonlocal nature of stress transfer in the elastic network.

3. Dynamic Homogenization in Periodic and Heterogeneous Solids

Advanced homogenization techniques such as plane-wave expansion (PWE) yield $k$-dependent moduli for periodic elastic media. In the formalism of Willis-constitutive equations (Norris et al., 2012), the effective moduli are written as

$$C^{\mathrm{eff}}_{ijkl}(k) = \langle C_{ijkl} \rangle - q^{\dagger}_{ijp}(k) G_{\neq0}(0, k) q_{klp}(k)$$

with $q_{ijp}(g;k) = (k+g)_m\, C_{ijpm}(g)$ and $G_{\neq0}$ the Green's operator over nonzero reciprocal lattice vectors $g$. The $k$-dependent transverse modulus, $M_T(k) = C^{\mathrm{eff}}_{1212}(k)$, encodes the nonlocal shear response of the periodic composite. In 1D multilayer systems, $M(k)$ shows band-structure-type oscillations and, in some exceptional cases, diverges at band edges but generally recovers the harmonic mean at $k \to 0$ (Norris et al., 2012).

4. Statistical and Microscopic Theories for Fluids, Amorphous Solids, and Soft Materials

In molecular fluids and amorphous systems, the nonlocal modulus is accessed via scale-dependent stress–stress correlations—Green–Kubo formalism generalized to wavevector space (Levashov, 2014, Uneyama, 2022, Zhou et al., 2023, Wittmer et al., 2023). For atomic or Brownian systems, the transverse modulus is obtained as:

$$G_T(k) = \frac{1}{k_B T V} \langle \hat{\sigma}_{yx}(k)\, \hat{\sigma}_{yx}(-k) \rangle,$$

where $\hat\sigma_{\alpha\beta}(k)$ is the Fourier transform of the Irving–Kirkwood stress tensor (Uneyama, 2022, Wittmer et al., 2023).

Key behaviors include:

• Liquids: $G_T(k=0) = 0$ as expected, but $G_T(k)$ is nonzero at finite $k$, reflecting finite-wavelength shear rigidity due to transient stress propagation, with $G_T(k)$ rising toward a plateau value at large $k$ (Levashov, 2014).
• Amorphous solids: Replica mean-field analysis relates $G_T(q)$ to the distribution of particle localization lengths $P(\xi)$; $G_T(q)$ exhibits $q^2$ softening for small $q$ and super-exponential decay at large $q$, directly reflecting the universal heterogeneity in the underlying structure (Zhou et al., 2023).

For overdamped Brownian systems, scale-dependent elastic moduli are defined via the linear response of the stress tensor to applied strain fields at finite $k$ (Uneyama, 2022). The simple harmonic dumbbell model exhibits a detailed, analytic expression for $G(k,0)$ with nontrivial $k$-dependence.

5. Extraction from Correlation Functions and Simulation

Direct extraction of $\mu(k)$ or $G_T(k)$ from simulations proceeds by analysis of strain or stress correlation functions in $k$-space. The transverse strain–strain correlation $C_T(q)$ in two-dimensional isotropic solids approaches a $q$-independent plateau for $q\to0$ proportional to $1/(4\mu)$, allowing the extraction of $\mu$ or, extended to finite $q$, a $k$-dependent $\mu(k)$ (Wittmer et al., 2023):

$\mu(k) = \frac{1}{4\beta V C_T(k)}$

Finite-$k$ effects are manifest in deviations that track the underlying structure factor or disorder, and angle dependence in tensor components is a coordinate artifact consistent with the 4th-order isotropic tensor analysis, not an indication of intrinsic anisotropy.

6. Wavevector Dependence in Viscoelastic and Non-Equilibrium Materials

In viscoelastic fluids, the Maxwell model can be extended to account for $k$-dependence by introducing two key $k$-dependent parameters: an effective Maxwell time $\tilde{\tau}_M(k)$ and a viscous response correction factor $w(k)$. The wavevector-dependent (complex) shear modulus is then (Hansen, 2024):

$$G^*(k, \omega) = \frac{i\omega \eta_0 w(k)}{1 + i\omega \tilde{\tau}_M(k)}$$

In the equilibrium (zero-frequency) limit, $M(k) = \lim_{\omega \to 0} G^*(k, \omega) = 0$ (for liquids), but for high frequency or at finite $k$, the plateau modulus $M_\infty(k)=\eta_0 w(k)/\tilde{\tau}_M(k)$ reflects the ability to sustain high-frequency shear waves. The $w(k)$ factor typically shows a pronounced minimum at a characteristic $k^*$, which defines a length scale $\ell^*=2\pi/k^*$ for maximal reduction in viscous dissipation.

7. Physical Interpretation and Relevance

The $k$-dependence of the equilibrium transverse modulus encodes:

• The crossover from macroscopic (hydrodynamic) elasticity to localized or nonlocal responses at mesoscopic and microscopic scales.
• The microscopic mechanisms—such as network connectivity, heterogeneity, relaxation time distributions, or viscoelasticity—that limit or enable shear rigidity at finite $k$.
• The scale over which shear waves, stress correlations, or particle localization are coherent and elastic, providing a probe of heterogeneity and the spatial range of elasticity in disordered, composite, or soft matter systems.

In summary, the wavevector-dependent equilibrium transverse modulus serves as a unifying descriptor for the nonlocal shear response of materials across network, periodic, fluid, amorphous, and soft-matter systems. It bridges continuum elasticity and microscopic statistical or molecular theories, with key analytical forms traced to effective medium theories (Katz et al., 2012), dynamic homogenization (Norris et al., 2012), atomic stress correlations (Levashov, 2014, Uneyama, 2022), replica mean-field elasticity (Zhou et al., 2023), and viscoelastic model refinements (Hansen, 2024).