Distorted Rayleigh Operator in Elasticity

• Distorted Rayleigh Operator is a microlocal tool that reduces elastic surface wave analysis to a scalar eigenvalue problem, capturing material anisotropy and boundary effects.
• Its construction involves microlocal diagonalization and explicit computation of principal and subprincipal symbols using boundary normal coordinates and the impedance tensor.
• The operator’s discrete spectrum governs surface wave frequencies, enabling precise quantization conditions with applications in seismology and materials science.

The distorted Rayleigh operator arises in the precise microlocal analysis of Rayleigh-type surface waves for linear elasticity in bodies with traction-free boundaries. This operator reduces the problem of constructing surface wave quasimodes in elastic media—possibly anisotropic—to an eigenvalue problem for a classical, selfadjoint, scalar first-order pseudo-differential operator defined intrinsically on the boundary. Its construction involves a reduction of the elastic wave equation via boundary normal coordinates, a microlocal diagonalization of the displacement-to-traction map, and the explicit computation of both its principal and subprincipal symbols. The operator’s discrete spectrum determines the frequencies of exponentially localized elastic surface waves, and its detailed structure reflects the anisotropy and curvature of the boundary and elasticity tensor (Hansen, 2010).

1. Elastic Wave Equation and Displacement–Traction Map

Let $M$ be a compact Riemannian manifold with smooth boundary $X = \partial M$, equipped with an elasticity tensor $C^{ijk\ell}$ (satisfying symmetry and strong convexity) and mass density $\rho > 0$. The vibrational analysis centers on the elastic wave equation,

$$L\,u = \lambda^2\,\rho\,u \text{ in } M, \quad T\,u = 0 \text{ on } X,$$

where $L$ is the elasticity operator and $T$ is the boundary traction, defined via Green’s formula: $$\int_M(C\,\operatorname{Def}\,u \mid \operatorname{Def}\,v)\,dV = \int_M(Lu \mid v)\,dV + \int_X(Tu \mid v)\,dA.$$ Introducing the semiclassical parameter $h = \lambda^{-1}$, one examines the shifted operator $h^2L - \rho$ and seeks a boundary-parametrix $B_h$ such that

$$(h^2L-\rho)\,B_h \equiv 0, \qquad B_h\big|_X = \operatorname{Id}$$

modulo rapidly decaying errors, for covectors in the elliptic (subsonic) region. The displacement–traction, or Neumann, map is then

$Z_h = h\,T\,B_h: C^\infty(X) \to C^\infty(X).$

2. Microlocal Reduction and the Scalar First-Order Operator

Microlocally, in the elliptic region, $Z_h$ acts as a first-order classical $h$-pseudodifferential operator on vector-valued boundary functions, with principal symbol given by the Hermitian impedance matrix $z(x,\xi)$. Under global assumptions (detailed in Section 5), $z(x,\xi)$ has a one-dimensional kernel along its characteristic set $\Sigma = \{ \det z(x,\xi) = 0 \}$. By microlocal conjugation, $Z_h$ may be block-diagonalized on the characteristic set so that, on the Rayleigh mode, it corresponds to a scalar operator: $Z_h \simeq (P_h - h^{-1})$ where $P_h$ is a selfadjoint operator on half-densities, and its leading symbol is independent of $h$ up to negligible corrections. In local coordinates, $P_h$ has the semiclassical pseudodifferential structure

$$P_h f(x) = (2\pi h)^{-(n-1)} \iint e^{\frac{i}{h}\langle x-y, \eta \rangle} \big(p(x,\eta) + h p_{\mathrm{sub}}(x,\eta)\big)\, f(y)\, dy\, d\eta + O(h^2).$$

$P_h$ is a classical, selfadjoint, elliptic first-order pseudodifferential operator on $X$.

3. Surface-Impedance Tensor and Principal Symbol

The principal symbol $p(x,\xi)$ of $P_h$ is constructed from the surface impedance tensor,

$$z(x, \xi) = i\,\big(c(x, \nu) q(x, \xi) + c(x, \nu, \xi)\big),$$

where $$c(x, \alpha, \beta)^i{}_k = C^{ij\ell m}(x)\, \alpha_j\, \beta_m\, \delta^i{}_k$$, $c(x, \alpha) = c(x,\alpha,\alpha)$, and $q(x,\xi)$ is the unique spectral factor solving

$$c(x, \xi + s\nu) - \rho\,\mathrm{Id} = (s-q^*)\,c(x,\nu)\,(s-q),\quad \operatorname{spec} q \subset \{ \Im s < 0 \}.$$

The characteristic set $\Sigma = \{\xi \in T^* X : \det z(x, \xi) = 0 \}$ is assumed to be intersected exactly once by each radial line $t \xi$ ($t>0$), permitting the definition of a function $p(x, \xi)$, positively homogeneous of degree one, with $\Sigma = p^{-1}(1)$. Thus,

$$\sigma_{\mathrm{prin}}(P_h)(x,\xi) = p(x,\xi) > 0, \quad p(x, t\xi) = t p(x,\xi).$$

This symbol encodes the local (geometric and material) propagation properties of Rayleigh-type surface waves.

4. Subprincipal Symbol and Eigenvector Bundle

The subprincipal symbol $p_{\mathrm{sub}}(x,\xi)$ encapsulates finer geometric information and involves the subprincipal datum $z_-$ of $Z_h$, the radial derivative $\dot{z}$ of $z$ with respect to homogeneity, and a globally defined unit eigenvector $v(x,\xi)$ of $z(x,\xi)$ on $\Sigma$: $$p_{\mathrm{sub}}(x,\xi) = (\dot{z} v \mid v)^{-1} \left\{ \Re(z_- v \mid v) + \Im \operatorname{tr}(v^*\, {}^h\nabla z \, \circ \, {}^v\nabla v) \right\}, \quad (x, \xi) \in \Sigma.$$ Here, ${}^h\nabla$ and ${}^v\nabla$ are horizontal and vertical derivatives on $T^*X$ respectively. The existence of a global nowhere-vanishing section $v(x,\xi)$ depends on the topological triviality of the kernel bundle $\ker\,z|_\Sigma$.

5. Global Hypotheses and Barnett–Lothe Conditions

The diagonalization and global construction of the distorted Rayleigh operator require:

• (U): For each $\xi \in E$, $z(x,\xi)$ has at most one non-positive eigenvalue (automatic in three dimensions).
• (E1): Every radial line $\{ t\xi : t > 0 \}$ meets the characteristic set $\Sigma$.
• (E2): The bundle $\ker\,z|_\Sigma$ is topologically trivial (admits a global unit section $v$).

These conditions are global variants of the Barnett–Lothe conditions. Under (U) and (E1), the zero of $\det z$ is simple and unique on each ray; (E2) enables continuous selection of $v$.

6. Eigenvalue Problem and Surface Quasimodes

The operator $P$ thus obtained is a classical, selfadjoint, elliptic first-order operator on $X$, whose discrete spectrum $\{\mu_j\}$ tends to infinity. For a complete orthonormal system $\{f_j\}$ of eigenfunctions of $P$,

$P f_j = \mu_j f_j$

one constructs corresponding quasimodes of the original elasticity problem,

$$u_j(r,x) = h_j^{-1/2} B_{h_j} J_{h_j} f_j, \quad h_j = \mu_j^{-1},$$

with $J_{h_j}$ a microlocal isometry, such that

$$(L - \mu_j^2 \rho) u_j = O(h_j^\infty) \text{ in } M, \quad T u_j = 0 \text{ on } X.$$

These $u_j$ are exponentially localized near the boundary and the exact quantization condition for Rayleigh surface quasimodes is $P f = \mu f$, with frequencies $\lambda_j = \mu_j$.

7. Implications and Applications

The formulation and explicit construction of the distorted Rayleigh operator provide an exact semi-classical reduction of elastic surface wave propagation to a spectral problem on the boundary, valid for general anisotropic elasticity tensors and Riemannian geometry, subject to the stated hypotheses. The computation of principal and subprincipal symbols enables direct analysis of surface wave dispersion and attenuation, with direct implications for the study of surface phenomena in seismology, materials science, and mechanical engineering (Hansen, 2010). The operator framework is particularly suited for study of spectral asymptotics, multiplicity, and localization properties of surface waves in anisotropic and inhomogeneous media.