Horizontal and Vertical Regularity of Elastic Wave Geometry

Basic Information

• Paper ID: 2511.16466
• Title: Horizontal and Vertical Regularity of Elastic Wave Geometry
• Authors: Joonas Ilmavirta, Pieti Kirkkopelto, Antti Kykkänen
• Classification: math.DG (Differential Geometry), math.AG (Algebraic Geometry), math.AP (Analysis of Partial Differential Equations)
• Publication Date: November 21, 2025 (arXiv preprint)
• Paper Link: https://arxiv.org/abs/2511.16466

Abstract

The elastic properties of materials are encoded by a stiffness tensor field, and the propagation of elastic waves is modeled by the elastic wave equation. This paper characterizes the analytical and algebraic properties that general anisotropic stiffness tensor fields must satisfy in order for Finsler geometric methods to be applicable to inverse problems related to elastic wave imaging.

Research Background and Motivation

Problem Background

1. Core Problem: Determine which regularity conditions on the stiffness tensor field ensure that elastic wave geometry possesses sufficient smoothness to allow Finsler geometric methods to be applied to inverse problems such as seismic imaging.
2. Significance:
   • Seismic wave propagation can be precisely described by geodesics in Finsler geometry
   • In practical applications, materials are often anisotropic and non-smooth
   • Reliable geometric models must be established under low regularity conditions
3. Limitations of Existing Methods:
   • Riemannian geometric methods can be used in the isotropic case
   • Classical Finsler geometry theory requires Finsler functions to be at least C^{1,1}
   • Microlocal analysis fails under finite regularity
   • Completely general Finsler geometry is too broad
4. Research Motivation:
   • Distinguish "horizontal regularity" (along manifold directions) from "vertical regularity" (along fiber directions)
   • Horizontal regularity depends on the regularity of the stiffness tensor along the manifold
   • Vertical regularity depends on the algebraic properties of the stiffness tensor
   • Provide precise criteria for the applicability of geometric methods to physical applications

Core Contributions

1. Horizontal Regularity Theorem (Theorem 1): Proves that if the stiffness tensor field c ∈ C^k(M) and the slowness surface of the qP branch is globally separated, then the associated Finsler function F^c_ is C^k along the manifold and smooth along the fiber.
2. Complete Characterization of Two-Dimensional Vertical Regularity (Theorem 2): Provides necessary and sufficient conditions for the slowness surface of a two-dimensional stiffness tensor to be real and separated: R(c) ≠ 0 and D(c) ≥ 0, where R and D are explicit polynomials in the stiffness tensor components.
3. High-Dimensional Vertical Regularity Results (Theorem 3): For n ∉ {2,4,8} dimensions, proves:
   • There exists an open neighborhood around each isotropic stiffness tensor where the qP branch is separated
   • The slowness surface of every complex stiffness tensor is singular as a complex scheme
4. Anisotropic Function Space Theory: Introduces new anisotropic regularity classes C^k_l and proves anisotropic implicit and inverse function theorems.
5. Inverse Problem Applications:
   • Corollary 5: Proves injectivity of the geodesic X-ray transform in low regularity elastic geometry
   • Corollary 6: Determines low regularity elastic wave geometry from travel time data up to Finsler isometry

Detailed Methodology

Task Definition

Input: Stiffness tensor field c: M → E_K(n), where M ⊂ ℝ^n is a smooth domain and E_K(n) is the space of stiffness tensors

Output: Determine the regularity class of the associated elastic wave geometry (Finsler function F^c_)

Constraints:

• Stiffness tensor satisfies elastic symmetry: c_ = c_ = c_
• The slowness surface of the qP branch is globally separated

Core Mathematical Framework

1. Anisotropic Regularity Classes

Definition: A function f: V × U → ℝ^d (V ⊂ ℝ^m, U ⊂ ℝ^n) belongs to the C^k_l class (1 ≤ k ≤ l) if and only if:

• f ∈ C^k(V × U) ∩ C^l(V × U)
• ∂^α_u f ∈ C^{min{k,l-|α|}}(V × U) for all |α| ≤ l
• ∂^α_v f ∈ C^{l-|α|}(V × U) for all |α| ≤ k

Physical Meaning: Captures different regularity along the manifold direction (v variable) and fiber direction (u variable).

2. Christoffel Matrix and Slowness Polynomial

The Christoffel matrix is defined as: $(\Gamma^c)_{il}(x,p) = c_{ijkl}(x)p_j p_k$

Its characteristic polynomial P_c(p) = det(Γ^c(p) - λId) is called the slowness polynomial.

The qP wave corresponds to the largest eigenvalue λ^c_, and the Finsler function is defined as: $F^c_{qP} = (\ell^{-1}_{1/2λ^c_{qP}})^*(λ^c_{qP})^{1/2}$

where ℓ_H is the map induced by the Legendre transform.

3. Horizontal Regularity Proof Strategy

Step 1: Prove the anisotropic inverse function theorem (Proposition 19)

• If f: A → B is a C^1 bundle diffeomorphism and f ∈ C^k_l(A;B)
• Then f^{-1} ∈ C^k_l(B;A)

Key Lemma 20: For functions on Euclidean domains, if:

• f ∈ C^k_l(V × U; ℝ^n)
• f(v,·) is a C^l-diffeomorphism
• f is fiber-independent in the v direction

Then f^{-1} ∈ C^k_l(V × f(U); ℝ^n)

Step 2: Prove the anisotropic implicit function theorem (Proposition 22)

• Apply to the slowness polynomial equation χ(Γ^c)(x,p,λ) = 0
• Since the qP branch is separated, ∂λχ(Γ^c)|{λ=λ^c_} ≠ 0
• Obtain λ^c_ ∈ C^k_∞(T^*M \ {0})

Step 3: Regularity of the Legendre transform (Proposition 24)

• If F̃ ∈ C^k_l(T^*M \ {0}) and l ≥ 3
• Then F = (ℓ^{-1}{1/2F̃^2})^*F̃ ∈ C^k(TM \ {0})

Complete Proof: Combining the above results yields F^c_ ∈ C^k_∞(TM \ {0}).

4. Vertical Regularity Analysis

Two-Dimensional Case: Uses resultant theory to analyze the system: