Introduction

The Leech lattice $\Lambda_{24}\subset \mathbb{R}^{24}$ represents an extraordinary confluence of sphere packing optimality, exceptional group theory, and modular forms . Normalized such that its minimal non-zero squared norm is $\min_{\mathbf{R}\in\Lambda_{24}\setminus\{\mathbf{0}\}}\left\lVert \mathbf{R} \right\rVert^2 = 4$, $\Lambda_{24}$ is the unique even unimodular lattice in dimension 24 containing no roots. Its universal optimality, established by Cohn, Kumar, Miller, Radchenko, and Viazovska , proves that $\Lambda_{24}$ minimizes potential energy among all unit-density 24-dimensional point configurations for every completely monotonic function of squared Euclidean distance.

In geometric analysis, every non-empty metric shell of $\Lambda_{24}$ forms a spherical $11$-design . In lattice dynamics, this design property forces discrete lattice sums in the Born–von Kárman expansion to match their continuous $\mathrm{O}(24)$-spherical averages for all polynomial terms of degree at most $11$, producing an acoustic tensor that is strictly isotropic through order $\|\mathbf{k}\|^8$.

The purpose of this paper is to establish the complete geometric, harmonic, and dynamical architecture of $\Lambda_{24}$ from its contact vectors to infinity. We address five foundational questions:

1. Polytopal Voronoi Facets: While literature frequently cites the kissing number $196{,}560$, we determine the complete facet census of $\operatorname{Vor}(\Lambda_{24})$ and prove why vectors of norm $\|R\|^2 \ge 8$ are strictly excluded.

2. Modular Channel Classification: We establish the double-cusp isomorphism $S_{12+d}^0 \cong M_{d-12}$ on $\mathrm{SL}_2(\mathbb{Z})$ governing all harmonic shell moments, and we prove why degree-$14$ moments vanish on every shell.

3. Completed $L$-Function Positivity: We evaluate the Fourier coefficients of $\Delta^2$ to order $25$, identify where sign alternation breaks down, and prove via Mellin transformation that $L(\Delta^2, s) > 0$ for all $s > 0$.

4. Acoustic Elasticity and Cauchy Reduction: We formalize the macroscopic elasticity under virial prestress and prove that zero-stress equilibrium algebraically forces $\lambda = \mu = S_2$, fixing the sound velocity ratio to $v_L / v_T = \sqrt{3}$.

5. Higher-Order Tensor Anisotropy: We derive the complete irreducible tensor decomposition at orders $\|\mathbf{k}\|^{10}$ and $\|\mathbf{k}\|^{12}$. We prove that while a two-shell model cancels order-$\|\mathbf{k}\|^{10}$ via $f''(6) = \frac{1}{48}f''(4)$, it leaves a residual transverse splitting at order $\|\mathbf{k}\|^{12}$. We then construct exact three-shell ($1080 : 45 : 1$) and infinite-crystal multi-Gaussian potentials that achieve simultaneous ultraisotropy through order $\|\mathbf{k}\|^{12}$.

Leech Lattice Geometry and Normalization

Let $\Lambda_{24}\subset \mathbb{R}^{24}$ denote the standard even unimodular Leech lattice normalized such that: $$\begin{equation} \min_{\mathbf{R} \in \Lambda_{24}\setminus \{\mathbf{0}\}} \left\lVert \mathbf{R} \right\rVert^2 = 4. \end{equation}$$

Because $\det(\Lambda_{24}) = 1$, the volume of the fundamental unit cell is $v_c = 1$. The equilibrium mass density with unit particle mass $M = 1$ is $\rho = M/v_c = 1$. The reciprocal lattice is $\Lambda_{\mathrm{rec}} = 2\pi \Lambda_{24}^* = 2\pi \Lambda_{24}$, and the first Brillouin zone is the Voronoi cell $\mathcal{B}_{24} \coloneqq \mathbb{R}^{24}/(2\pi \Lambda_{24})$.

The theta series of $\Lambda_{24}$ is: $$\begin{equation} \Theta_{\Lambda_{24}}(q) = \sum_{m=0}^\infty |S_m| q^m = E_{12}(q) - \frac{65520}{691}\Delta(q) = 1 + 196560 q^2 + 16773120 q^3 + 398034000 q^4 + \dots, \end{equation}$$ where $S_m \coloneqq \{ \mathbf{R} \in \Lambda_{24}: \left\lVert \mathbf{R} \right\rVert^2 = 2m \}$. Because $\Lambda_{24}$ contains no roots, $S_1 = \emptyset$ ($|S_1| = 0$). For every $m \ge 2$, the exact shell population is: $$\begin{equation} \label{eq:shell_pop} |S_m| = \frac{65520}{691} \left( \sigma_{11}(m) - \tau(m) \right), \end{equation}$$ where $\sigma_{11}(m) = \sum_{d|m} d^{11}$ and $\tau(m)$ is Ramanujan’s tau function.

Exact Voronoi Facet Geometry of $\operatorname{Vor}(\Lambda_{24})$

The Voronoi cell of $\Lambda_{24}$ is the convex polytope: $$\begin{equation} \operatorname{Vor}(\Lambda_{24}) \coloneqq \left\{ \mathbf{x} \in \mathbb{R}^{24} : \mathbf{x} \cdot \mathbf{R} \le \frac{1}{2}\left\lVert \mathbf{R} \right\rVert^2 \quad \forall \mathbf{R} \in \Lambda_{24}\setminus \{\mathbf{0}\} \right\}. \end{equation}$$

A vector $\mathbf{v} \in \Lambda_{24}\setminus \{\mathbf{0}\}$ is Voronoi-relevant if its bounding hyperplane $$H_{\mathbf{v}} \coloneqq \left\{ \mathbf{x} \in \mathbb{R}^{24} : \mathbf{x} \cdot \mathbf{v} = \frac{1}{2}\left\lVert \mathbf{v} \right\rVert^2 \right\}$$ intersects $\operatorname{Vor}(\Lambda_{24})$ in a 23-dimensional facet $$F_{\mathbf{v}} = \operatorname{Vor}(\Lambda_{24}) \cap H_{\mathbf{v}}.$$

Covering-Radius Upper Bound

Voronoi Coset Uniqueness and Facet Census

Harmonic Multipoles to Infinity: $S_{12+d}^0 \cong M_{d-12}$

Let $\operatorname{Harm}_d(\mathbb{R}^{24})$ denote the space of homogeneous harmonic polynomials of degree $d$. For $P_d\in\operatorname{Harm}_d(\mathbb{R}^{24})$, the weighted theta series is $$\begin{equation} \Theta_{\Lambda_{24},P_d}(\tau) \coloneqq \sum_{\mathbf{R}\in\Lambda_{24}} P_d(\mathbf{R})q^{\left\lVert \mathbf{R} \right\rVert^2/2}, \qquad q=e^{2\pi i\tau}, \quad \tau\in\mathbb{H}. \end{equation}$$

By the Hecke–Schoeneberg theorem , for $d\ge1$, $\Theta_{\Lambda_{24},P_d}$ is a cusp form of weight $w=12+d$ on $\mathrm{SL}_2(\mathbb{Z})$.

Because $\min_{\mathbf{R}\neq0}\left\lVert \mathbf{R} \right\rVert^2=4$, the $q^1$ Fourier coefficient vanishes: $$\begin{equation} \Theta_{\Lambda_{24},P_d}(\tau) = \sum_{m=2}^{\infty} \left( \sum_{\mathbf{R}\in S_m}P_d(\mathbf{R}) \right)q^m = a_2q^2+a_3q^3+\dots \implies a_1=0. \end{equation}$$

We define the double-cusp subspace: $$\begin{equation} S_w^0(\mathrm{SL}_2(\mathbb{Z})) \coloneqq \left\{ F\in S_w(\mathrm{SL}_2(\mathbb{Z})):\ a_1(F)=0 \right\}. \end{equation}$$

Degree-12 Fourier Coefficients and Completed $L$-Function

For $d=12$, $M_0(\mathrm{SL}_2(\mathbb{Z}))=\mathbb{C}$, so $\dim S_{24}^0(\mathrm{SL}_2(\mathbb{Z}))=1$.

The unique generator is $$\Delta^2(\tau) = \sum_{m=2}^{\infty}c_m(\Delta^2)q^m,$$ where $$c_m(\Delta^2) = \sum_{j=1}^{m-1}\tau(j)\tau(m-j).$$ Consequently, on every shell: $$\begin{equation} M_m(12)=c_m(\Delta^2)M_2(12). \end{equation}$$

Continuum Phonon Dynamics and Cauchy–Lamé Reduction

Consider particles interacting via an admissible pair potential $V(r)=f(r^2)$.

The Born–von Kárman dynamical matrix $D(\mathbf{k})\in\operatorname{Mat}_{24}(\mathbb{C})$ is: $$\begin{equation} \label{eq:dynamical_matrix_def} \boxed{ D_{ab}(\mathbf{k}) \coloneqq \sum_{\mathbf{R}\in\Lambda_{24}\setminus\{\mathbf{0}\}} \bigl(1-\cos(\mathbf{k}\cdot\mathbf{R})\bigr) \left[ 2\delta_{ab}f'(\left\lVert \mathbf{R} \right\rVert^2) + 4R_aR_bf''(\left\lVert \mathbf{R} \right\rVert^2) \right]. } \end{equation}$$

Expanding $$1-\cos(\mathbf{k}\cdot\mathbf{R}) = \frac12(\mathbf{k}\cdot\mathbf{R})^2 + O(\left\lVert \mathbf{k} \right\rVert^4),$$ the quadratic acoustic tensor is: $$\begin{equation} D_{ab}^{(2)}(\mathbf{k}) = \sum_{c,d=1}^{24} C_{acbd}k_ck_d, \end{equation}$$ where $$\begin{equation} C_{acbd} = \sum_{\mathbf{R}\neq\mathbf{0}} R_cR_d \left[ \delta_{ab}f'(\left\lVert \mathbf{R} \right\rVert^2) + 2R_aR_bf''(\left\lVert \mathbf{R} \right\rVert^2) \right]. \end{equation}$$

Define the radial shell sums: $$\begin{equation} S_1 \coloneqq \frac{1}{24} \sum_{\mathbf{R}\neq\mathbf{0}} \left\lVert \mathbf{R} \right\rVert^2f'(\left\lVert \mathbf{R} \right\rVert^2), \end{equation}$$ and $$\begin{equation} S_2 \coloneqq \frac{1}{312} \sum_{\mathbf{R}\neq\mathbf{0}} \left\lVert \mathbf{R} \right\rVert^4f''(\left\lVert \mathbf{R} \right\rVert^2). \end{equation}$$

Design-Controlled Ultraisotropy Through $\|\mathbf{k}\|^8$

Expanding $$D(\mathbf{k}) = \sum_{j=1}^{\infty}D^{(2j)}(\mathbf{k}),$$ the order-$\|\mathbf{k}\|^{2j}$ coefficient is: $$\begin{equation} D_{ab}^{(2j)}(\mathbf{k}) = \frac{(-1)^{j-1}}{(2j)!} \sum_{\mathbf{R}\neq\mathbf{0}} (\mathbf{k}\cdot\mathbf{R})^{2j} \left[ 2\delta_{ab}f'(\left\lVert \mathbf{R} \right\rVert^2) + 4R_aR_bf''(\left\lVert \mathbf{R} \right\rVert^2) \right]. \end{equation}$$

The summands contain homogeneous polynomials in $\mathbf{R}$ of degrees $2j$ and $2j+2$.

Order $\|\mathbf{k}\|^{10}$ and the Exact Two-Shell Cancellation

At order $\|\mathbf{k}\|^{10}$ ($j=5$), the $f'$ term has degree $10\le11$ and is isotropic. The $f''$ term has degree $12$.

Using $$R_aR_b(\mathbf{k}\cdot\mathbf{R})^{10} = \frac{1}{132} \partial_a\partial_b (\mathbf{k}\cdot\mathbf{R})^{12}$$ and the series prefactor $\frac{4}{10!}$: $$\begin{equation} \label{eq:D10_exact} \boxed{ D_{\mathrm{aniso},ab}^{(10)}(\mathbf{k}) = \frac{1}{33\cdot10!} \mathcal{A}_{12}[f''] \, \partial_a\partial_bH_{12}^{(2)}(\mathbf{k}), } \end{equation}$$ where $$\begin{equation} \mathcal{A}_{12}[f''] \coloneqq \sum_{m=2}^{\infty} c_m(\Delta^2)f''(2m). \end{equation}$$

Minimal-Shell Non-12-Design Certificate: Exact Combinatorial Derivation

In the integer realization $$\Lambda_{24}^{\mathrm{int}}=\sqrt{8}\,\Lambda_{24},$$ vectors in the minimal shell $X_{\mathrm{int}}$ have squared norm $32$.

The $196{,}560$ vectors partition into three classical shapes:

• Type 1: $(\pm4^2,0^{22})$, multiplicity $$2^2\binom{24}{2}=1{,}104.$$

• Type 2: $(\pm2^8,0^{16})$, supported on the $759$ octads of the extended binary Golay code $\mathcal{G}_{24}$, multiplicity $$2^7\times759=97{,}152.$$

• Type 3: $(\mp3^1,\pm1^{23})$, odd vectors with sign parity, multiplicity $$24\times2^{12}=98{,}304.$$

We evaluate the moment $$M_{12}^{\mathrm{int}}(\mathbf{u}) \coloneqq \sum_{\mathbf{R}\in X_{\mathrm{int}}} (\mathbf{u}\cdot\mathbf{R})^{12}$$ along $$u_A=e_1$$ and $$u_B=\frac{e_1+e_2}{\sqrt2}.$$

1. Evaluation along $u_A=e_1$:

   • Type 1: Coordinate $1$ is non-zero in $2\times23\times2=92$ vectors, contributing $$92\times4^{12} = 1{,}543{,}503{,}872.$$

   • Type 2: In the Steiner system $S(5,8,24)$, the number of octads containing a fixed point is $$759\times\frac{8}{24}=253.$$ With $128$ sign configurations per octad, $$253\times128=32{,}384$$ vectors have $R_1=\pm2$, contributing $$32{,}384\times2^{12} = 132{,}644{,}864.$$

   • Type 3: In $4{,}096$ vectors, coordinate $1$ is $\pm3$: $$4{,}096\times3^{12} = 2{,}176{,}782{,}336.$$ In the remaining $$23\times4{,}096=94{,}208$$ vectors, coordinate $1$ is $\pm1$, contributing $$94{,}208.$$ Thus Type 3 contributes $$2{,}176{,}876{,}544.$$

   Summing all three types: $$\begin{equation} M_{12}^{\mathrm{int}}(u_A) = 1{,}543{,}503{,}872 + 132{,}644{,}864 + 2{,}176{,}876{,}544 = 3{,}853{,}025{,}280. \end{equation}$$

2. Evaluation along $u_B=(e_1+e_2)/\sqrt2$:

   Here $$(\mathbf{u}\cdot\mathbf{R})^{12} = \frac{1}{64}(R_1+R_2)^{12}.$$

   • Type 1: Vectors with support on $\{1,2\}$ contribute $$2\times8^{12}.$$ Vectors with support on $\{1,k\}$ or $\{2,k\}$, $k\ge3$, contribute $$176\times4^{12}.$$ Dividing by $64$ gives $$2{,}193{,}620{,}992.$$

   • Type 2: In $S(5,8,24)$, exactly $$759\times\frac{28}{276}=77$$ octads contain both $1$ and $2$, and $$2(253-77)=352$$ octads contain exactly one of them. Summing over sign parities gives $$1{,}294{,}729{,}216.$$

   • Type 3: Summing over coordinate parities gives $$1{,}076{,}887{,}552.$$

   Hence $$\begin{equation} M_{12}^{\mathrm{int}}(u_B) = 2{,}193{,}620{,}992 + 1{,}294{,}729{,}216 + 1{,}076{,}887{,}552 = 4{,}565{,}237{,}760. \end{equation}$$

Taking the difference yields the exact harmonic moment: $$\begin{equation} \boxed{ H_{12}^{\mathrm{int}}(u_B) - H_{12}^{\mathrm{int}}(u_A) = M_{12}^{\mathrm{int}}(u_B) - M_{12}^{\mathrm{int}}(u_A) = 712{,}212{,}480 \neq0. } \end{equation}$$

For $$\mathbf{u} = \frac{1}{\sqrt5}(1,2,0,\dots),$$ the transverse Hessian splits into eigenspaces of multiplicities $1$ and $22$ with exact eigenvalues $$\begin{equation} \lambda_{\mathrm{in}} = \frac{254700453888}{625}, \qquad \lambda_{\mathrm{out}} = \frac{209217945600}{625}, \end{equation}$$ giving the exact non-zero gap: $$\begin{equation} \Delta\lambda_{\mathrm{int}} = \lambda_{\mathrm{in}} - \lambda_{\mathrm{out}} = \frac{45482508288}{625} \neq0. \end{equation}$$

Exact Two-Shell Cancellation

Restricting Eq. [eq:D10_exact] to the first two shells $S_2\cup S_3$ ($c_2=1$, $c_3=-48$): $$\begin{equation} \mathcal{A}_{12}^{(2\text{-shell})} = f''(4)-48f''(6) = 0. \end{equation}$$

Therefore $$\begin{equation} \boxed{ f''(6)=\frac{1}{48}f''(4). } \end{equation}$$

Under this condition, the order-$\|\mathbf{k}\|^{10}$ anisotropic tensor vanishes identically.

Order $\|\mathbf{k}\|^{12}$ Irreducible Tensor Decomposition

At $j=6$, $$D_{ab}^{(12)}(\mathbf{k}) = T_{ab}^{(12,1)}(\mathbf{k}) + T_{ab}^{(12,2)}(\mathbf{k}).$$

Generic Transverse Splitting

For any unit transverse polarization $\mathbf{v}\in\mathbf{k}^{\perp}$, $$\mathbf{v}\cdot\mathbf{k}=0, \qquad \left\lVert \mathbf{v} \right\rVert=1,$$ we have $$\begin{equation} \mathbf{v}^{\top} D_{\mathrm{aniso}}^{(12)}(\mathbf{k}) \mathbf{v} = -\frac{1}{12!} \left[ 2\mathcal{F}_1 + \frac16\mathcal{F}_2 \right] H_{12}^{(2)}(\mathbf{k}) - \frac{\left\lVert \mathbf{k} \right\rVert^2}{12\cdot12!} \mathcal{F}_2[f''] \, \mathbf{v}^{\top} \bigl(\nabla\nabla H_{12}^{(2)}(\mathbf{k})\bigr) \mathbf{v}. \end{equation}$$

Multiscale Annihilation: Two-Shell, Three-Shell, and Multi-Gaussian

Residual Splitting in the Two-Shell Model

In the two-shell model ($S_2\cup S_3$): $$\begin{equation} \mathcal{A}_{12} = f''(4)-48f''(6) = 0 \implies f''(6)=\frac{1}{48}f''(4). \end{equation}$$

Evaluating $\mathcal{F}_2$ under this condition: $$\begin{align} \mathcal{F}_2 &= 2f''(4)-144f''(6) \\ &= 2f''(4) - 144\left(\frac{1}{48}f''(4)\right) \\ &= 2f''(4)-3f''(4) \\ &= \boxed{-f''(4)\neq0}. \end{align}$$