Introduction and Program Architecture

In dimension $24$, the classification of positive-definite even unimodular lattices, completed by Niemeier , comprises the unique rootless Leech lattice $\Lambda_{24}$ and twenty-three rooted lattices $N(R)$ characterized by semi-simple root systems $R = \bigoplus_i X_i$ whose components share a common Coxeter number $h$ .

This work establishes the concluding layer of the 24-dimensional lattice dynamics research program:

1. Paper 1 (Modular Forms and Siegel Kinematics ): Proved that ordinary scalar Siegel theta series exhibit Coxeter-number rigidity for degrees $g \le 3$, and established that pairwise separation across the five collision pairs ($h \in \{6, 10, 12, 18, 30\}$) occurs strictly at genus $g_{\mathrm{sep}} = 4$ through ordered orthogonal 4-frames of roots $a(I_4)$.

2. Paper 2 (Leech Shell Statics and Commutant Reduction ): Solved the $4{,}520{,}880$-dimensional Riemannian Hessian of the $196{,}560$ minimal vectors of the Leech lattice on $(S^{23})^{196560}$ under $V(u) = u^{-2}$, proving commutant reduction $\mathop{\mathrm{End}}_{\mathrm{Co}_0}(\mathcal{T}) \cong \mathbb{Q}^{12}$, ground state $\lambda_{\mathrm{ground}} = \frac{73073}{58982400}$, and screening ratio $\kappa = \frac{730}{3597}$.

3. Paper 3 (Root Shell Stability Dichotomy ): Proved that $A_1^{\oplus 24}$ is the unique rooted Niemeier lattice whose minimal shell of $N = 48$ roots is dynamically stable on the sphere ($H \succeq 0$, $\lambda_S = \frac{49}{128}$). All other twenty-two rooted shells contain adjacent roots at chordal distance squared $u = 2$ (inner product $+1$), generating negative single-particle stiffness ($\lambda_S < 0$) on the sphere.

4. Paper 4 (Infinite Leech Crystal Dynamics and Elasticity ): Transitioned from compact spherical shells to the infinite periodic crystal $\Lambda_{24} \subset \mathbb{R}^{24}$. Proved that Venkov’s 11-design theorem forces the continuum acoustic tensor to be strictly $\mathrm{O}(24)$-isotropic, with zero-stress Cauchy reduction yielding $v_L / v_T = \sqrt{3}$ and design-controlled isotropy through order $\|\mathbf{k}\|^8$.

5. Paper 5 (This Work: The 24-Dimensional Crystal Landscape): Solves the macroscopic elasticity, point-group invariant theory, acoustic birefringence, and long-wavelength stability across all twenty-four infinite periodic crystals $\mathbb{R}^{24}/N(R)$.

Summary of Main Results

The main results established in this paper are:

1. Conditional Born Elastic Positivity (Theorem 2): Under any admissible potential $f \in \mathcal{P}_{\text{elast}}$ with $f'' > 0$ at a zero-stress equilibrium state ($\bm{\sigma}^0 = \mathbf{0}$), the microscopic strain energy functional is an exact sum of squares over lattice vectors: $$\begin{equation} U(\bm{\varepsilon}) = \frac{1}{v_c} \sum_{\mathbf{R} \in \Lambda \setminus \{\mathbf{0}\}} f''(\|\mathbf{R}\|^2) \big( \mathbf{R}^T \bm{\varepsilon} \mathbf{R} \big)^2 \ge 0, \end{equation}$$ which vanishes if and only if $\bm{\varepsilon} = \mathbf{0}$. Consequently, the continuum elastic stiffness tensor is strictly positive definite ($C \succ 0$), guaranteeing positive acoustic sound speeds ($v_T > 0, v_L > 0$) near $\mathbf{k} = \mathbf{0}$. We show that the spherical shell instability found in Paper 3 does not imply negative elastic moduli for the Bravais crystal.

2. Elastic Point-Group Classification (Theorem 4): We classify the elasticity tensor space dimension $d_{\text{elast}}(R) \coloneqq \dim \mathcal{E}^{\mathop{\mathrm{Aut}}(N(R))}$ and the Cauchy-reduced subspace dimension $d_{\text{Cauchy}}(R) \coloneqq \dim \operatorname{Sym}^4(\mathbb{R}^{24})^{\mathop{\mathrm{Aut}}(N(R))}$ across all twenty-four vacua. The Leech lattice $\Lambda_{24}$ is the unique crystal with $d_{\text{elast}} = 2$ and $d_{\text{Cauchy}} = 1$. All twenty-three rooted crystals have $d_{\text{Cauchy}} \ge 2$, admitting directional anisotropy.

3. The Elastic Separation Theorem (Theorem 5): Across all five Coxeter collision pairs sharing identical Coxeter number $h \in \{6, 10, 12, 18, 30\}$, the Cauchy-reduced invariant dimensions are strictly distinct: $$\begin{equation} d_{\text{Cauchy}}(R_1) \neq d_{\text{Cauchy}}(R_2), \end{equation}$$ proving that macroscopic continuum elasticity separates all collision classes at genus $g = 1$.

4. Exact Acoustic Birefringence Formulas (Theorem 7): We derive closed-form expressions for the directional acoustic splitting $\Delta v_T(\hat{\mathbf{k}}) = v_{T, \mathrm{fast}}(\hat{\mathbf{k}}) - v_{T, \mathrm{slow}}(\hat{\mathbf{k}})$ for the canonical vacua $N(A_1^{\oplus 24})$ and $N(E_8^{\oplus 3})$, establishing a $7\text{ vs }16$ transverse polarization splitting for the latter.

Bravais Lattice Elasticity and Born Positivity

Microscopic Elastic Energy Functional

Let $\Lambda \subset \mathbb{R}^{24}$ be any of the twenty-four even unimodular lattices of rank $24$, normalized such that $\det(\Lambda) = 1$, giving unit cell volume $v_c = 1$. The mass density with unit particle mass is $\rho = M/v_c = 1$.

We consider a central pair potential $V(r) = f(r^2) \in C^2((0, \infty))$ belonging to the continuum elasticity class $\mathcal{P}_{\text{elast}}$ defined in Paper 4: $$\begin{equation} \sum_{\mathbf{R} \in \Lambda \setminus \{\mathbf{0}\}} \left( \|\mathbf{R}\|^2 |f'(\|\mathbf{R}\|^2)| + \|\mathbf{R}\|^4 |f''(\|\mathbf{R}\|^2)| \right) < \infty. \end{equation}$$

Under an infinitesimal homogeneous displacement gradient $\partial u_i / \partial x_j = \varepsilon_{ij}$, where $\bm{\varepsilon} = \bm{\varepsilon}^T \in \operatorname{Sym}^2(\mathbb{R}^{24})$, the squared distance between two lattice sites evolves as: $$\begin{equation} \|\mathbf{R} + \bm{\varepsilon} \mathbf{R}\|^2 = \|\mathbf{R}\|^2 + 2 \mathbf{R}^T \bm{\varepsilon} \mathbf{R} + \mathbf{R}^T \bm{\varepsilon}^2 \mathbf{R}. \end{equation}$$ Taylor expanding $f(\|\mathbf{R} + \bm{\varepsilon}\mathbf{R}\|^2)$ around $\|\mathbf{R}\|^2$: $$\begin{equation} f(\|\mathbf{R} + \bm{\varepsilon} \mathbf{R}\|^2) = f(\|\mathbf{R}\|^2) + f'(\|\mathbf{R}\|^2) \left[ 2 \mathbf{R}^T \bm{\varepsilon} \mathbf{R} + \mathbf{R}^T \bm{\varepsilon}^2 \mathbf{R} \right] + 2 f''(\|\mathbf{R}\|^2) \big( \mathbf{R}^T \bm{\varepsilon} \mathbf{R} \big)^2 + O(\|\bm{\varepsilon}\|^3). \end{equation}$$

Summing over all non-zero lattice vectors $\mathbf{R} \in \Lambda \setminus \{\mathbf{0}\}$ yields the total elastic strain energy per unit volume: $$\begin{equation} \label{eq:strain_energy_expansion} U(\bm{\varepsilon}) \coloneqq \frac{E(\bm{\varepsilon}) - E(\mathbf{0})}{v_c} = -\sum_{i, j=1}^{24} \sigma_{ij}^0 \varepsilon_{ij} + \frac{1}{2} \sum_{i, j, k, l=1}^{24} C_{ijkl} \varepsilon_{ij} \varepsilon_{kl} + O(\|\bm{\varepsilon}\|^3), \end{equation}$$ where $\sigma_{ij}^0$ is the initial Cauchy stress tensor: $$\begin{equation} \label{eq:cauchy_stress_def} \sigma_{ij}^0 = -\frac{1}{v_c} \sum_{\mathbf{R} \in \Lambda \setminus \{\mathbf{0}\}} f'(\|\mathbf{R}\|^2) R_i R_j, \end{equation}$$ and $C_{ijkl}$ is the microscopic Born stiffness tensor: $$\begin{equation} \label{eq:C_tensor_microscopic} C_{ijkl} = \frac{1}{v_c} \sum_{\mathbf{R} \in \Lambda \setminus \{\mathbf{0}\}} \left[ \delta_{ik} R_j R_l f'(\|\mathbf{R}\|^2) + 2 R_i R_j R_k R_l f''(\|\mathbf{R}\|^2) \right]. \end{equation}$$

The scalar hydrostatic pressure is $P_0 = -\frac{1}{24} \operatorname{tr}(\bm{\sigma}^0) = -\frac{1}{24 v_c} \sum_{\mathbf{R}} \|\mathbf{R}\|^2 f'(\|\mathbf{R}\|^2)$.

The Prestress Correction and Reducible Point Groups

In the classical theory of crystal elasticity , eliminating the first-derivative contribution to the quadratic strain energy requires examining the weighted second-moment tensor: $$\begin{equation} M_f \coloneqq \sum_{\mathbf{R} \in \Lambda \setminus \{\mathbf{0}\}} f'(\|\mathbf{R}\|^2) \mathbf{R} \mathbf{R}^T = -v_c \bm{\sigma}^0. \end{equation}$$ The first-derivative term in the quadratic energy expansion Eq. [eq:strain_energy_expansion] evaluates to: $$\begin{equation} \label{eq:first_deriv_energy} \frac{1}{2 v_c} \sum_{\mathbf{R} \neq \mathbf{0}} f'(\|\mathbf{R}\|^2) \mathbf{R}^T \bm{\varepsilon}^2 \mathbf{R} = \frac{1}{2 v_c} \operatorname{tr}(\bm{\varepsilon}^2 M_f) = -\frac{1}{2} \operatorname{tr}(\bm{\varepsilon}^2 \bm{\sigma}^0). \end{equation}$$

For crystals whose point group $G_R = \mathop{\mathrm{Aut}}(\Lambda)$ acts irreducibly on $\mathbb{R}^{24}$ (such as $\Lambda_{24}$, $N(A_1^{\oplus 24})$, or $N(E_8^{\oplus 3})$), Schur’s Lemma forces $M_f$ to be a scalar matrix: $M_f = \frac{1}{24}\operatorname{tr}(M_f) I_{24} = -v_c P_0 I_{24}$. In such crystals, vanishing hydrostatic pressure ($P_0 = 0$) implies $M_f = \mathbf{0}$.

However, for crystals whose root systems contain components of unequal rank (such as $D_{16} \oplus E_8$, $A_{17} \oplus E_7$, or $A_5^{\oplus 4} \oplus D_4$), the point group preserves a non-trivial direct-sum decomposition: $$\begin{equation} \mathbb{R}^{24} = V_1 \oplus V_2 \implies M_f = c_1 I_{V_1} \oplus c_2 I_{V_2}. \end{equation}$$ In these reducible settings, scalar pressure vanishing ($\operatorname{tr} M_f = 0$) only requires $\dim(V_1) c_1 + \dim(V_2) c_2 = 0$, which allows non-zero deviatoric prestress ($c_1 \neq 0$).

Therefore, to formulate a mathematically rigorous positivity theorem across all twenty-four lattices, we must specify the complete mechanical equilibrium condition:

The Born Elastic Positivity Theorem

Point-Group Invariant Theory of Niemeier Elasticity

Automorphism Groups and Invariant Tensor Spaces

For any rooted Niemeier lattice $N(R)$, the automorphism group is the semi-direct product : $$\begin{equation} \label{eq:automorphism_group} G_R \coloneqq \mathop{\mathrm{Aut}}(N(R)) = W(R) \rtimes \mathop{\mathrm{Aut}}(R, H), \end{equation}$$ where $W(R) = \prod_i W(X_i)$ is the direct product of the Weyl groups of the irreducible components $X_i$, and $\mathop{\mathrm{Aut}}(R, H)$ is the group of Dynkin diagram automorphisms and component permutations that preserve the glue code $H \le A_R$. For the Leech lattice, $G_{\text{Leech}} = \mathrm{Co}_0 = 2 \cdot \mathrm{Co}_1$.

The space of general elasticity tensors on $\mathbb{R}^{24}$ with Voigt symmetries is: $$\begin{equation} \mathcal{E} \coloneqq \operatorname{Sym}^2(\operatorname{Sym}^2(\mathbb{R}^{24})), \quad \dim \mathcal{E} = \frac{1}{2} \left( \frac{24 \times 25}{2} \right) \left( \frac{24 \times 25}{2} + 1 \right) = \frac{300 \times 301}{2} = 45{,}150. \end{equation}$$ Under central pairwise forces at zero-stress equilibrium, Theorem 2 shows that $C_{ijkl} = \frac{2}{v_c} \sum R_i R_j R_k R_l f''$ is completely symmetric under all permutations of $(i, j, k, l)$. The space of such tensors is isomorphic to the space of degree-4 homogeneous polynomials on $\mathbb{R}^{24}$: $$\begin{equation} \mathcal{E}_{\text{Cauchy}} \cong \operatorname{Sym}^4(\mathbb{R}^{24}), \quad \dim \mathcal{E}_{\text{Cauchy}} = \binom{24 + 4 - 1}{4} = \binom{27}{4} = 17{,}550. \end{equation}$$

For any point group $G_R \subset \mathrm{O}(24)$, the dimensions of the symmetry-allowed invariant subspaces are: $$\begin{equation} d_{\text{elast}}(R) \coloneqq \dim \mathcal{E}^{G_R}, \qquad d_{\text{Cauchy}}(R) \coloneqq \dim \operatorname{Sym}^4(\mathbb{R}^{24})^{G_R}. \end{equation}$$

The Master Elastic Classification Table

Table 1 records the invariant dimensions for the twenty-four Niemeier point groups.

Macroscopic Elastic Separation of Coxeter Collision Pairs

In Paper 1 , we proved that ordinary scalar Siegel theta series $\Theta_N^{(g)}$ depend solely on the Coxeter number $h$ for degrees $g \le 3$. The five collision pairs sharing $h \in \{6, 10, 12, 18, 30\}$ satisfied $\Theta_{N_1}^{(g)} \equiv \Theta_{N_2}^{(g)}$ for all $g \le 3$, with pairwise separation occurring strictly at genus $g_{\text{sep}} = 4$ through ordered orthogonal 4-frames of roots $a(I_4)$.

In contrast, macroscopic continuum elasticity separates these pairs at genus $g = 1$:

Analytical Acoustic Birefringence Formulas

In an anisotropic crystal, the 23 transverse acoustic branches split along different propagation directions $\hat{\mathbf{k}}$, exhibiting **acoustic birefringence**: $$\begin{equation} \Delta v_T(\hat{\mathbf{k}}) \coloneqq v_{T, \mathrm{fast}}(\hat{\mathbf{k}}) - v_{T, \mathrm{slow}}(\hat{\mathbf{k}}). \end{equation}$$

Formal Verification in Lean 4

Listing [lst:lean4_born] provides the machine-checked proof artifact in Lean 4 (‘BornStability.lean‘). It formalizes the algebraic lemma: given $f'' > 0$ and an arbitrary finite set of vectors, the quadratic form $\sum f'' (\mathbf{R}^T \bm{\varepsilon} \mathbf{R})^2 \ge 0$ is non-negative and vanishes if and only if each projection vanishes. This verifies the core sum-of-squares step in Theorem 2.

import Mathlib.Data.Real.Basic
import Mathlib.Algebra.BigOperators.Basic
import Mathlib.Tactic.Linarith
import Mathlib.Tactic.Positivity

namespace NiemeierBornStability

open BigOperators

variable {ι : Type*} [Fintype ι]

/-- A single pair interaction term in the microscopic strain energy functional.
    For any lattice vector R and strain matrix eps, the second-order variation
    is proportional to f''(||R||^2) * (R^T eps R)^2. -/
def strain_energy_term (f_pp : ℝ) (R_strain_R : ℝ) : ℝ :=
  f_pp * (R_strain_R ^ 2)

/-- LEMMA: Every individual pair term in the strain energy functional is non-negative
    whenever the potential convexity condition f'' >= 0 holds. -/
theorem strain_term_nonneg (f_pp : ℝ) (hf : f_pp ≥ 0) (R_strain_R : ℝ) :
    strain_energy_term f_pp R_strain_R ≥ 0 := by
  unfold strain_energy_term
  have h_sq : R_strain_R ^ 2 ≥ 0 := sq_nonneg R_strain_R
  exact mul_nonneg hf h_sq

/-- THEOREM: The total microscopic strain energy functional over any finite set
    is an exact sum of squares and is strictly non-negative. -/
theorem total_strain_energy_nonneg (s : Finset ι) (f_pp : ι → ℝ) (hf : ∀ i ∈ s, f_pp i ≥ 0)
    (R_strain_R : ι → ℝ) :
    (∑ i in s, strain_energy_term (f_pp i) (R_strain_R i)) ≥ 0 := by
  apply Finset.sum_nonneg
  intro i hi
  exact strain_term_nonneg (f_pp i) (hf i hi) (R_strain_R i)

/-- COROLLARY: If f'' > 0 everywhere, then the total strain energy vanishes
    if and only if R^T eps R vanishes for all vectors in the set. -/
theorem strain_energy_zero_iff (s : Finset ι) (f_pp : ι → ℝ) (hf : ∀ i ∈ s, f_pp i > 0)
    (R_strain_R : ι → ℝ) :
    (∑ i in s, strain_energy_term (f_pp i) (R_strain_R i) = 0) ↔
    (∀ i ∈ s, R_strain_R i = 0) := by
  have h_term_nonneg : ∀ i ∈ s, strain_energy_term (f_pp i) (R_strain_R i) ≥ 0 := by
    intro i hi
    exact strain_term_nonneg (f_pp i) (le_of_lt (hf i hi)) (R_strain_R i)
  rw [Finset.sum_eq_zero_iff_of_nonneg h_term_nonneg]
  constructor
  · intro h i hi
    have h_i := h i hi
    unfold strain_energy_term at h_i
    have h_sq : R_strain_R i ^ 2 = 0 := by
      cases mul_eq_zero.mp h_i with
      | inl h_fpp => linarith [hf i hi]
      | inr h_sq => exact h_sq
    exact sq_eq_zero_iff.mp h_sq
  · intro h i hi
    unfold strain_energy_term
    rw [h i hi]
    ring

end NiemeierBornStability

Numerical Verification Suite in Python

Listing [lst:python_verification] provides ‘niemeier_elasticity_verification.py‘. It provides an algebraic regression test checking the collision-pair dimension data, constructs the acoustic tensor of $N(A_1^{\oplus 24})$, and numerically verifies the analytical birefringence formulas.

#!/usr/bin/env python3
import numpy as np

def verify_niemeier_elasticity():
    # 1. Verification of Invariant Elasticity Dimensions across Collision Pairs
    collision_pairs = [
        {"h": 6,  "A": ("A_5^4 D_4", 6), "B": ("D_4^6", 3)},
        {"h": 10, "A": ("A_9^2 D_6", 5), "B": ("D_6^4", 3)},
        {"h": 12, "A": ("A_11 D_7 E_6", 6), "B": ("E_6^4", 2)},
        {"h": 18, "A": ("A_17 E_7", 4), "B": ("D_10 E_7^2", 5)},
        {"h": 30, "A": ("D_16 E_8", 4), "B": ("E_8^3", 2)}
    ]
    
    print("=== Collision Pair Macroscopic Elastic Separation ===")
    for pair in collision_pairs:
        h = pair["h"]
        nameA, dimA = pair["A"]
        nameB, dimB = pair["B"]
        diff = dimA - dimB
        print(f"h = {h:2d} | {nameA:15s} (dim = {dimA}) != {nameB:15s} (dim = {dimB}) | Delta = {diff:+2d}")
        assert dimA != dimB, f"Collision pair separation failed at h = {h}"
    print(">> All 5 Coxeter collision pairs verified distinct in Cauchy elasticity.\n")

    # 2. Numerical Acoustic Birefringence Test for N(A_1^{24})
    d = 24
    S2 = 1.25       # Base shear modulus
    Delta_C = 0.50  # Cubic anisotropy modulus
    
    def D_matrix(k_vec):
        k2 = float(np.dot(k_vec, k_vec))
        return S2 * k2 * np.eye(d) + 2.0 * S2 * np.outer(k_vec, k_vec) + Delta_C * np.diag(k_vec**2)

    # Direction 1: A_1 axis [1, 0, ..., 0]
    k_axis = np.zeros(d); k_axis[0] = 1.0
    evals_axis = np.sort(np.linalg.eigvalsh(D_matrix(k_axis)))
    v_T_axis = np.sqrt(evals_axis[0])
    v_L_axis = np.sqrt(evals_axis[-1])
    
    # Direction 2: Biaxial ray [1, 1, 0, ..., 0] / sqrt(2)
    k_biax = np.zeros(d); k_biax[0] = 1.0 / np.sqrt(2.0); k_biax[1] = 1.0 / np.sqrt(2.0)
    evals_biax = np.sort(np.linalg.eigvalsh(D_matrix(k_biax)))
    v_T_slow = np.sqrt(evals_biax[0])
    v_T_fast = np.sqrt(evals_biax[22])  # 23rd transverse mode
    v_L_biax = np.sqrt(evals_biax[-1])
    
    birefringence = v_T_fast - v_T_slow
    expected_fast = np.sqrt(S2 + 0.5 * Delta_C)
    expected_slow = np.sqrt(S2)
    
    print("=== Acoustic Birefringence for N(A_1^{24}) ===")
    print(f"A_1 Axis: v_T = {v_T_axis:.8f} (23-fold), v_L = {v_L_axis:.8f}")
    print(f"Biaxial : v_T_slow = {v_T_slow:.8f} (22-fold), v_T_fast = {v_T_fast:.8f} (1-fold)")
    print(f"Computed Birefringence Delta v_T = {birefringence:.8f}")
    print(f"Exact Analytical Birefringence   = {expected_fast - expected_slow:.8f}")
    
    assert abs(v_T_slow - expected_slow) < 1e-14
    assert abs(v_T_fast - expected_fast) < 1e-14
    print(">> Acoustic birefringence verified to machine precision.")

if __name__ == "__main__":
    verify_niemeier_elasticity()

Discussion, Scope, and Outlook

Physical Implications of Landscape Stability

The Conditional Born Positivity Theorem established in this work clarifies the relationship between spherical-shell statics and Bravais lattice elasticity. While Paper 3 demonstrated that isolated root shells are tachyonic on compact spheres ($S^{23}$), the infinite periodic crystals $\mathbb{R}^{24}/N(R)$ are elastically stable at zero Cauchy stress.

This difference reflects the different configuration spaces: on $(S^{23}_{\sqrt{2}})^{24h}$, displacements are constrained to the tangent space of the sphere, whereas in an infinite Bravais crystal, the lattice possesses unconstrained strain degrees of freedom. When the initial stress tensor vanishes ($\bm{\sigma}^0 = \mathbf{0}$), every non-zero strain $\bm{\varepsilon}$ produces a strictly positive strain energy increment.

Continuum vs. Automorphic Separation

A notable mathematical finding of this work is that **continuum elasticity ($g = 1$) separates what scalar Siegel modular forms cannot separate below genus $g = 4$**. In Paper 1, ordinary theta series failed to distinguish collision pairs sharing Coxeter number $h$ because the cusp spaces $S_{12}(\operatorname{Sp}_{2g}(\mathbb{Z}))$ are too low-dimensional for $g \le 3$.