For a rank-nn Euclidean lattice LL, we write min⁡L=min⁡0≠x∈L⟨x,x⟩\min L=\min_{0\ne x\in L}\langle x,x\rangle for the minimum squared length of nonzero vectors in LL; the kissing number

τ(L)=∣{x∈L:⟨x,x⟩=min⁡L}∣\tau(L)=\left|\left\{x\in L:\langle x,x\rangle=\min L\right\}\right|

counts the number of such vectors, including both signs. The maximal lattice kissing number in dimension nn is

Klat(n)=max⁡rank⁡L=nτ(L).K_{\mathrm{lat}}(n)=\max_{\operatorname{rank}L=n}\tau(L).

In this note, we obtain a new asymptotic upper bound for Klat(n)K_{\mathrm{lat}}(n), using explicit auxiliary functions constructed in OpenAI’s preprint [15].

Theorem 1. As n→∞n\to\infty,

Klat(n)≤(2eπ+o(1))n=2(12log⁡2(2e/π)+o(1))n.K_{\mathrm{lat}}(n)\leq\left(\sqrt{\frac{2e}{\pi}}+o(1)\right)^n =2^{\left(\frac{1}{2}\log_{2}(2e/\pi)+o(1)\right)n}.

Here 12log⁡2(2e/π)=0.395599455708322…\frac{1}{2}\log_{2}(2e/\pi)=0.395599455708322\ldots.

For a packing of congruent balls, the contact degree of a ball is the number of other balls tangent to it; the average contact degree of a nonempty finite packing is the arithmetic mean of these degrees. After scaling such balls to radius 1/21/2, their centers form a finite set C⊂RnC\subset\mathbb{R}^{n} with minimum separation at least 1, and the average contact degree is

τ‾(C)=#{(x,y)∈C2:∥x−y∥=1}∣C∣.\overline{\tau}(C)=\frac{\#\left\{(x,y)\in C^{2}:\lVert x-y\rVert=1\right\}}{|C|}.

Corollary 4 gives the same upper bound as Theorem 1 for τ‾(C)\overline{\tau}(C), uniformly over all such CC.

To place the bound in context, let K(n)K(n) be the ordinary kissing number—i.e., the largest number of nonoverlapping unit balls that can touch one unit ball in Rn\mathbb{R}^{n}. Equivalently, K(n)K(n) is the largest size of a set of unit vectors with pairwise inner products at most 1/21/2. Normalizing the shortest vectors of a lattice gives such a set, so Klat(n)≤K(n)K_{\mathrm{lat}}(n) \le K(n). The classical spherical-code bound of Kabatiansky and Levenshtein [12], from 1978, gives the exponent 0.40141…0.40141\ldots for K(n)K(n). Combining it with Sidel’nikov’s spherical-cap reduction [18] gives the optimized classical bound [12], with exponent 0.400944…0.400944\ldots; subsequent work of Sardari and Zargar [17], and of Zargar [22], improved multiplicative factors.

Fourier methods offer a second route, through the contacts in a packing. In 2020, Afkhami-Jeddi, Cohn, Hartman, de Laat, and Tajdini [1] extrapolated an implied kissing exponent of 0.39560.3956 (the d=∞d=\infty row of [1]) from their numerically optimized density bound. Its interpretation as an average-contact bound relied on interpolation and sign conjectures. Their Theorem 5.3 gives an unconditional inequality for admissible auxiliary functions; we apply this result to the explicit functions of [15], obtaining an exponent that agrees with their extrapolation to the displayed precision, without relying on the intermediate conjectures.

The cap-optimized spherical-code hierarchy of [15] gives

K(n)≤2(0.39661+o(1))nK(n) \le2^{(0.39661+o(1))n}

at level 2; see [15]. Our Theorem 1 and Corollary 4 give a smaller exponent for lattice kissing numbers and average contact degree, respectively. In fact, Lemma A.1 in Appendix A shows that the infimum of the cap-optimized fixed-level asymptotic kissing exponents of the hierarchy in [15] remains strictly above ours. To our knowledge, Theorem 1 gives the first proven asymptotic upper bound for lattice kissing numbers whose exponential rate is strictly smaller than the best known upper-bound rate for the unrestricted kissing number.2

Meanwhile, there remains a large gap between the upper and lower bounds for lattice kissing numbers. The Barnes–Wall lattices [3] in dimensions n=2mn=2^{m} have kissing numbers

∏i=1m(2i+2)=exp⁡(Θ((log⁡n)2));\prod_{i=1}^{m}(2^{i}+2)=\exp(\Theta((\log n)^{2}));

see Nebe [14]. Recent work of Laarhoven and the present author [13] improves this to

Klat(n)≥(12πe2+o(1))ne2nK_{\mathrm{lat}}(n) \ge\left(\frac{1}{2\pi e^{2}}+o(1)\right)\sqrt{n}e^{2\sqrt{n}}

as n→∞n\to\infty through prime powers.3

Our analytic input comes from Chapter 1 of [15]; all references to that preprint use the version updated August 6, 2026. That chapter determines the exponential strength of the Fourier linear program for sphere-packing density, developed by Gorbachev [11] and Cohn–Elkies [6]. Writing LPn\mathrm{LP}_{n} for the infimum of the density upper bounds supplied by this program, [15] states

lim⁡n→∞LPn1/n=e2π.\lim_{n\to\infty}\mathrm{LP}_{n}^{1/n}=\sqrt{\frac{e}{2\pi}}.

Consequently, the maximal packing density Δn\Delta_{n} satisfies Δn≤2−(0.604400544…+o(1))n\Delta_{n}\le2^{-(0.604400544\ldots+o(1))n}; the equality above determines the rate of the linear program itself. A comparison theorem of Cohn and Zhao [7] then implies that the ordinary spherical linear program for K(n)K(n) cannot attain an exponent smaller than 12log⁡2(2e/π)\frac{1}{2}\log_{2}(2e/\pi).

For our bounds, we use the explicit auxiliary functions ff constructed in [15] to establish the asymptotic density upper bound, rather than the density-rate theorem alone. The density bound depends on f(0)/f^(0)f(0)/\widehat{f}(0) together with a volume factor, whereas the kissing and contact bounds depend on f(0)/[−f(R)]f(0)/[-f(R)], where RR is the center-to-center distance at contact. We evaluate the latter ratio using the saddle-point estimate at the target parameter u0u_{0} in [15], Chapter 1, §4. Lemma 3 below collects the required properties of these functions and the resulting contact-value estimate, with precise references to the global sign and saddle-point estimates on which our argument relies.

All logarithms without a subscript are natural.

The Fourier contact bound

The following bound for average contact degree appears in [1], Theorem 5.3. Its finite-packing assertion also follows from the positive-type bound of Dostert, Kolpakov, and de Oliveira Filho [9], Theorem 4.1, specialized to congruent balls. As noted in [1], it is the Euclidean analogue of Fortier Bourque and Petri’s trace-formula bound [10], Proposition 4.1.

We use the normalization f^(ξ)=∫Rnf(x)e−2πi⟨x,ξ⟩ dx\widehat{f}(\xi)=\int_{\mathbb{R}^{n}}f(x)e^{-2\pi i\langle x,\xi\rangle}\,\mathrm{d}x. For radial ff, write f(r)f(r) for its value at a vector of norm rr.

Lemma 2 ([1], Theorem 5.3; [9], Theorem 4.1). Let ff be a real radial Schwartz function on Rn\mathbb{R}^{n}, and suppose that, for some R>0R>0,

f^≥0,f(r)≤0(r≥R),f(R)<0.\widehat{f}\ge0,\qquad f(r)\le0\quad(r\ge R),\qquad f(R)<0.

Then every rank-nn lattice satisfies τ(L)≤f(0)/[−f(R)]\tau(L)\le f(0)/[-f(R)]. The same ratio bounds the average contact degree of every nonempty finite packing of balls of radius R/2R/2, with contacts at center distance RR.

Proof. We scale LL so that its shortest nonzero vectors have length RR, equivalently min⁡L=R2\min L=R^{2}. Writing L∗L^{*} for the dual lattice, Poisson summation gives

0≤1covol⁡(L)∑ξ∈L∗f^(ξ)=∑x∈Lf(x)≤f(0)+τ(L)f(R);0\le\frac{1}{\operatorname{covol}(L)}\sum_{\xi\in L^{*}}\widehat{f}(\xi) =\sum_{x\in L}f(x) \le f(0)+\tau(L)f(R);

since f(R)<0f(R)<0, this proves the lattice bound.

For a nonempty finite set CC of minimum separation at least RR, Fourier inversion gives

0≤∫Rnf^(ξ)∣∑x∈Ce2πi⟨x,ξ⟩∣2 dξ=∑x,y∈Cf(x−y).0\le\int_{\mathbb{R}^{n}}\widehat{f}(\xi)\left|\sum_{x\in C}e^{2\pi i\langle x,\xi\rangle}\right|^{2}\,\mathrm{d}\xi =\sum_{x,y\in C}f(x-y).

If MM counts the ordered pairs at distance RR, then the last sum is at most ∣C∣f(0)+Mf(R)|C|f(0)+Mf(R). Thus M/∣C∣≤f(0)/[−f(R)]M/|C|\le f(0)/[-f(R)]. □\square

The auxiliary functions

Fix sufficiently small ε>0\varepsilon>0, and put

λ=n2,u=1+ε4,m=λ(1+u)2,ψ=Γ′Γ.\lambda=\frac{n}{2},\qquad u=1+\frac{\varepsilon}{4},\qquad m=\frac{\lambda(1+u)}{2},\qquad \psi=\frac{\Gamma'}{\Gamma}.

Let wεw_{\varepsilon} be the compactly supported signed density in [15], Chapter 1, §4,4 and define the real function

gε(t)=∫0∞wε(a)(cosh⁡(at)−1) da.g_{\varepsilon}(t)=\int_{0}^{\infty}w_{\varepsilon}(a)(\cosh(at)-1)\,\mathrm{d}a.

Thus gε(t)=hε(it)g_{\varepsilon}(t)=h_{\varepsilon}(it) in the notation of [15], Chapter 1. We collect the analytic input, including the contact value extracted from the [15], Chapter 1 saddle estimate, in the following lemma.

Lemma 3. For every fixed sufficiently small ε>0\varepsilon>0, there is n0(ε)n_{0}(\varepsilon) such that, for every n≥n0(ε)n\ge n_{0}(\varepsilon), there exist a real radial Schwartz function f=fn,εf=f_{n,\varepsilon} on Rn\mathbb{R}^{n} and a radius R=Rn,ε>0R=R_{n,\varepsilon}>0 satisfying the hypotheses of Lemma 2, with

log⁡R=−12log⁡π+12ψ(m)+gε′(u),(1)\log R=-\frac{1}{2}\log\pi+\frac{1}{2}\psi(m)+g_{\varepsilon}'(u), \tag*{(1)}
f(0)=2(u−1)πλ/2eλgε(1),(2)f(0)=2(u-1)\pi^{\lambda/2}e^{\lambda g_{\varepsilon}(1)}, \tag*{(2)}
−f(R)=An,επ−λu/2Γ(m)eλgε(u)R−2m,An,ε≍εn,(3)-f(R)=A_{n,\varepsilon}\pi^{-\lambda u/2}\Gamma(m)e^{\lambda g_{\varepsilon}(u)}R^{-2m}, \qquad A_{n,\varepsilon}\asymp_{\varepsilon}\sqrt{n}, \tag*{(3)}
gε′(u)=−12log⁡(π/2)+O(ε).(4)g_{\varepsilon}'(u)=-\frac{1}{2}\log(\pi/2)+O(\varepsilon). \tag*{(4)}

In particular, An,ε>0A_{n,\varepsilon}>0. The constants implicit in ≍ε\asymp_{\varepsilon} may depend on ε\varepsilon but not on nn.5 More precisely, put

pε=(u−1)[1−(1+u)2]<0,Vn,ε=λ4ψ′(m)+gε′′(u).p_{\varepsilon}=(u-1)[1-(1+u)^{2}]<0, \qquad V_{n,\varepsilon}=\frac{\lambda}{4}\psi'(m)+g_{\varepsilon}''(u).

Then Vn,ε≍ε1V_{n,\varepsilon}\asymp_{\varepsilon}1 and

An,ε=(−pε)λ2πVn,ε(1+oε(1)).(5)A_{n,\varepsilon}=(-p_{\varepsilon})\sqrt{\frac{\lambda}{2\pi V_{n,\varepsilon}}}\left(1+o_{\varepsilon}(1)\right). \tag*{(5)}

Proof. We take ff to be the function denoted f−f_{-} in [15], Chapter 1, §4. The origin value (2) and radius (1) are given in [15], Chapter 1, (42) and (44), respectively. The derivative estimate (4) is [15], Chapter 1, Lemma 4.2, and the global sign conditions are established in [15], Chapter 1, Theorem 4.1. For the contact value, the target parameter is u=u0u=u_{0} in that preprint’s notation, and this endpoint is included in [15], Chapter 1, Lemma 4.8. The saddle coefficient pεp_{\varepsilon} is given in [15], Chapter 1, (38)–(39). In [15], Chapter 1, (70)–(71), the inverse Mellin factor λ/(2π)\lambda/(2\pi) and Gaussian factor combine as

λ2π2πλVn,ε=λ2πVn,ε,\frac{\lambda}{2\pi}\sqrt{\frac{2\pi}{\lambda V_{n,\varepsilon}}} = \sqrt{\frac{\lambda}{2\pi V_{n,\varepsilon}}},

giving (3) and (5). The curvature bounds in [15], Chapter 1, Lemma 4.5, (56) give Vn,ε≍ε1V_{n,\varepsilon}\asymp_{\varepsilon}1. Since pε<0p_{\varepsilon}<0, the contact value is negative and An,ε≍εnA_{n,\varepsilon}\asymp_{\varepsilon}\sqrt{n}. □

We also need the estimate

gε(u)−gε(1)=O(ε).(6)g_{\varepsilon}(u)-g_{\varepsilon}(1)=O(\varepsilon). \tag*{(6)}

To verify this, let wε,+w_{\varepsilon,+} and wε,−w_{\varepsilon,-} be the positive and negative parts of wεw_{\varepsilon}, so that wε=wε,+−wε,−w_{\varepsilon}=w_{\varepsilon,+}-w_{\varepsilon,-}. The density construction and [15], Chapter 1, Lemma 4.2 give

wε,−(a)≤e−2a2a2cosh⁡a,∫0∞wε,+(a)asinh⁡(ua) da=O(e−c/ε2)w_{\varepsilon,-}(a)\le\frac{e^{-2a}}{2a^{2}\cosh a}, \qquad \int_{0}^{\infty}w_{\varepsilon,+}(a)a\sinh(ua)\,\mathrm{d}a = O(e^{-c/\varepsilon^{2}})

for an absolute c>0c>0. For 1≤t≤u≤5/41\le t\le u\le5/4, we use the inequalities sinh⁡(ta)≤tacosh⁡(ta)\sinh(ta)\le ta\cosh(ta) and cosh⁡(ta)≤e(t−1)acosh⁡a\cosh(ta)\le e^{(t-1)a}\cosh a to obtain

∣gε′(t)∣≤t2∫0∞e−(3−t)a da+O(e−c/ε2)≤514+O(e−c/ε2).|g_{\varepsilon}'(t)| \le \frac{t}{2}\int_{0}^{\infty}e^{-(3-t)a}\,\mathrm{d}a + O(e^{-c/\varepsilon^{2}}) \le \frac{5}{14} + O(e^{-c/\varepsilon^{2}}).

Integration over [1,u][1,u], of length ε/4\varepsilon/4, proves (6).

The kissing number bound

We now estimate the ratio f(0)/[−f(R)]f(0)/[-f(R)] and apply Lemma 2 to obtain the lattice kissing number bound.

Proof of Theorem 1. Fix sufficiently small ε>0\varepsilon>0. Taking logarithms in (2) and (3) and substituting (1) gives

log⁡f(0)−f(R)=λ[gε(1)−gε(u)+(1+u)gε′(u)]+mψ(m)−log⁡Γ(m)+log⁡2(u−1)An,ε.(7)\begin{aligned} \log\frac{f(0)}{-f(R)} &= \lambda\left[g_{\varepsilon}(1)-g_{\varepsilon}(u)+(1+u)g_{\varepsilon}'(u)\right] \\ &\quad+m\psi(m)-\log\Gamma(m)+\log\frac{2(u-1)}{A_{n,\varepsilon}}. \tag*{(7)} \end{aligned}

(The powers of π\pi cancel because 2m=λ(1+u)2m=\lambda(1+u).) Stirling’s formula yields

mψ(m)−log⁡Γ(m)=m+12log⁡m+O(1).m\psi(m)-\log\Gamma(m)=m+\frac{1}{2}\log m+O(1).

Since log⁡An,ε=12log⁡n+Oε(1)\log A_{n,\varepsilon}=\frac{1}{2}\log n+O_{\varepsilon}(1) and m/n=(1+u)/4m/n=(1+u)/4, we have

12log⁡m+log⁡2(u−1)An,ε=Oε(1).\frac{1}{2}\log m+\log\frac{2(u-1)}{A_{n,\varepsilon}}=O_{\varepsilon}(1).

Dividing (7) by n=2λn=2\lambda therefore gives

1nlog⁡f(0)−f(R)=12[gε(1)−gε(u)+(1+u)gε′(u)+1+u2]+Oε(n−1).\frac{1}{n}\log\frac{f(0)}{-f(R)} =\frac{1}{2}\left[g_{\varepsilon}(1)-g_{\varepsilon}(u)+(1+u)g_{\varepsilon}'(u)+\frac{1+u}{2}\right]+O_{\varepsilon}(n^{-1}).

Since u=1+O(ε)u=1+O(\varepsilon), equations (4) and (6) then give

1nlog⁡f(0)−f(R)=12log⁡2eπ+O(ε)+Oε(n−1).\frac{1}{n}\log\frac{f(0)}{-f(R)} =\frac{1}{2}\log\frac{2e}{\pi}+O(\varepsilon)+O_{\varepsilon}(n^{-1}).

Applying Lemma 2 and taking the upper limit as n→∞n\to\infty with ε\varepsilon fixed, we obtain

lim sup⁡n→∞1nlog⁡Klat(n)≤12log⁡2eπ+O(ε).\limsup_{n\to\infty}\frac{1}{n}\log K_{\mathrm{lat}}(n) \leq\frac{1}{2}\log\frac{2e}{\pi}+O(\varepsilon).

Letting ε↓0\varepsilon\downarrow0 proves the theorem. □\square

Average contact degree

The same ratio estimate also bounds average contact degree.

Corollary 4. As n→∞n\to\infty, uniformly over nonempty finite C⊂RnC\subset\mathbb{R}^{n} with minimum separation at least 11,

τˉ(C)≤(2eπ+o(1))n.\bar{\tau}(C)\leq\left(\sqrt{\frac{2e}{\pi}}+o(1)\right)^{n}.

Proof. Fix sufficiently small ε>0\varepsilon>0, and take ff and RR from Lemma 3. The scaled set RCRC has minimum separation at least RR, and its ordered pairs at distance RR correspond exactly to the ordered pairs of CC at distance 11. The finite-packing part of Lemma 2 therefore gives

τˉ(C)≤f(0)−f(R);\bar{\tau}(C)\leq\frac{f(0)}{-f(R)};

the estimate for this ratio in the proof of Theorem 1 is independent of CC. Taking n→∞n\to\infty with ε\varepsilon fixed and then letting ε↓0\varepsilon\downarrow0 proves the stated uniform bound. □\square

Concluding remarks

Uniformity and contact counts

Explicitly, for every η>0\eta> 0 there is n0(η)n_{0}(\eta) such that, for every n≥n0(η)n \ge n_{0}(\eta), both Klat(n)K_{\mathrm{lat}}(n) and τ‾(C)\overline{\tau}(C) for every eligible C⊂RnC \subset\mathbb{R}^{n} are at most 2(12log⁡2(2e/π)+η)n2^{(\frac{1}{2}\log_{2}(2e/\pi)+\eta)n}. For a packing of NN congruent balls, the contact number—the number of unordered tangent pairs—is therefore at most

N2(2eπ+o(1))n,\frac{N}{2}\left(\sqrt{\frac{2e}{\pi}}+o(1)\right)^{n},

uniformly in NN.

Our bounds apply to lattice kissing numbers and to the average contact degree of arbitrary finite packings of congruent balls. They do not establish the same exponential bound for the maximum contact degree in such a packing, and therefore do not give the same bound for the unrestricted kissing number K(n)K(n); they also do not address packings with unequal radii. For background on the contact-number problem, see the Bezdek and Khan [5] survey.

Comparison with the spherical linear program

The exponent in Theorem 1 also gives a lower bound on the exponential rate of the optimal ordinary kissing-number linear-programming bound. Let KnLPK_{n}^{\mathrm{LP}} denote the infimum of the upper bounds for K(n)K(n) supplied by the two-point spherical linear program of Delsarte, Goethals, and Seidel [8]. Cohn and Zhao [7] specialized to the kissing angle π/3\pi/3, proved

LPn≤2−nKnLP;(8)\mathrm{LP}_{n} \le2^{-n}K_{n}^{\mathrm{LP}}; \tag*{(8)}

their comparison applies in the same dimension nn.6 Combining (8) with lim⁡n→∞LPn1/n=e/(2π)\lim_{n\to\infty}\mathrm{LP}_{n}^{1/n}=\sqrt{e/(2\pi)} gives

lim inf⁡n→∞1nlog⁡2KnLP≥12log⁡22eπ.(9)\liminf_{n\to\infty}\frac{1}{n}\log_{2}K_{n}^{\mathrm{LP}} \ge\frac{1}{2}\log_{2}\frac{2e}{\pi}. \tag*{(9)}

Thus our bounds for lattice kissing numbers and average contact degree attain an exponent below which the ordinary spherical linear program cannot bound the unrestricted kissing number.7 The inequality (9) is a lower bound on the optimal spherical linear programming upper bound KnLPK_{n}^{\mathrm{LP}}, not a lower bound on the kissing number K(n)K(n) itself. Earlier asymptotic lower bounds for the optimal spherical linear-programming bound were obtained by Samorodnitsky [16].

The same constant appears within the spherical-code hierarchy of [15], whose certificates belong to this two-point linear program; see [15]. Specifically, write σ=2Γr\sigma=2\Gamma_{r} for a certificate’s threshold and Φ=Φr\Phi=\Phi_{r} for its code-size exponent. Then [15] gives

Φ=12log⁡221−σ−12log⁡22πe+DKL(ν∥ρK)2log⁡2,\Phi=\frac{1}{2}\log_{2}\frac{2}{1-\sigma}-\frac{1}{2}\log_{2}\frac{2\pi}{e}+\frac{D_{\mathrm{KL}}(\nu\Vert\rho_{K})}{2\log2},

where ν\nu and ρK\rho_{K} are the probability measures defined there and DKLD_{\mathrm{KL}} denotes relative entropy. For parameters feasible at the kissing angle, σ≥1/2\sigma\ge1/2; nonnegativity of relative entropy therefore gives Φ≥12log⁡2(2e/π)\Phi\ge\frac{1}{2}\log_{2}(2e/\pi). This identity makes the same lower bound explicit for those certificates, while the Cohn–Zhao comparison applies to the entire two-point spherical linear program.

Appendix A gives a stronger limitation for the spherical-code hierarchy of [15]: Lemma A.1 shows that, even after spherical-cap optimization and taking the infimum over fixed levels, the hierarchy’s asymptotic kissing exponent remains strictly above ours. This statement concerns that certificate family; it does not establish a strict gap in (9) for the full two-point spherical linear program.

The sharpness argument in [15] Chapter 2, Theorem 8.3 sends σ→1\sigma\to1, establishing the optimized sphere-packing density exponent. This is compatible with Lemma A.1, which fixes the kissing angle at π/3\pi/3.

References

Appendix A. Separation from the OpenAI spherical-code hierarchy

In this appendix, we show that the asymptotic kissing-number exponents supplied by the hierarchy of [15] Chapter 2 remain uniformly separated from the exponent in Theorem 1, even after spherical-cap optimization and taking the infimum over fixed levels. Starting from the relative-entropy identity of [15], Chapter 2, Proposition 8.1, we use Jensen’s inequality to bound the mass of a fixed interval and thereby obtain an explicit positive lower bound on the relative entropy.

Lemma A.1. Put

α∗=12log⁡22eπ,β=23−2πarctan⁡123.\alpha_{*}=\frac{1}{2}\log_{2}\frac{2e}{\pi},\qquad \beta=\frac{2}{3}-\frac{2}{\pi}\arctan\frac{1}{2\sqrt{3}}.

The infimum of the fixed-level asymptotic kissing-number exponents supplied by the hierarchy of [15], Chapter 2, Theorem 1.2, over all finite levels and the spherical-cap parameters in [15], Chapter 2, (89), is at least

α∗−14log⁡2(4β(1−β))=0.395815863377…>α∗.\alpha_{*}-\frac{1}{4}\log_{2}\left(4\beta(1-\beta)\right)=0.395815863377\ldots>\alpha_{*}.

Here the level is fixed before the dimension tends to infinity, and the infimum is then taken over the resulting asymptotic rates. The argument does not address levels that vary with the dimension.

Proof of Lemma A.1. For a finite-level certificate, let xi≥0x_i\geq0 and Ri>0R_i>0 be the quadratic coordinates and residues in [15], Chapter 2, Proposition 8.1, (102), so that ∑iRi=1\sum_i R_i=1. Define

wi=xixi+1/4,μ=∑iRiδwi,w_i=\sqrt{\frac{x_i}{x_i+1/4}},\qquad \mu=\sum_i R_i\delta_{w_i},

where δw\delta_w denotes unit point mass at ww, and put

kw(v)=(1−v2)(1−w2)1−v2+v2w2,F(v)=∫01kw(v) dμ(w)(0≤v<1).k_w(v)=\frac{(1-v^2)(1-w^2)}{1-v^2+v^2w^2},\qquad F(v)=\int_0^1 k_w(v)\,\mathrm{d}\mu(w)\qquad(0\leq v<1).

In particular, k0(v)=1k_0(v)=1 and k1(v)=0k_1(v)=0. Under t=(1−v2)/4t=(1-v^2)/4, the reference measures in [15], Chapter 2, Proposition 8.1 become

dν(v)=dv,dρ(v)=2 dvπ1−v2.\mathrm{d}\nu(v)=\mathrm{d}v,\qquad \mathrm{d}\rho(v)=\frac{2\,\mathrm{d}v}{\pi\sqrt{1-v^2}}.

The kernel satisfies

∫01kw(v) dρ(v)=1−w(0≤w≤1).\int_0^1 k_w(v)\,\mathrm{d}\rho(v)=1-w\qquad(0\leq w\leq1).

We define the normalizing constant

Z≔∫01F(v) dρ(v);Z\coloneqq\int_0^1 F(v)\,\mathrm{d}\rho(v);

the kernel identity gives

Z=1−∫01w dμ(w).Z=1-\int_0^1 w\,\mathrm{d}\mu(w).

By [15], Chapter 2, (102)–(104), the certificate’s threshold σ\sigma and exponent Φ\Phi satisfy

σ=∫01w dμ(w),Z=1−σ,Φ=−12log⁡2∫01log⁡F(v) dv.(10)\begin{aligned} \sigma&=\int_0^1 w\,\mathrm{d}\mu(w),& Z&=1-\sigma,& \Phi&=-\frac{1}{2\log2}\int_0^1\log F(v)\,\mathrm{d}v. \tag*{(10)} \end{aligned}

We now enlarge the class under consideration by allowing arbitrary probability measures μ\mu on [0,1][0,1], including endpoint atoms, rather than only those arising from finite-level certificates. For these measures, define FF and ZZ as above and σ\sigma and Φ\Phi by the same formulas in (10). The kernel identity ensures that Z=1−σZ=1-\sigma continues to hold. We set Φ=+∞\Phi=+\infty if F≡0F\equiv0.

We first include cap optimization and reduce to σ=1/2\sigma=1/2. A certificate used at an inner-product threshold 0≤s≤1/20\leq s\leq1/2 is feasible when σ≥s\sigma\geq s; cap reduction adds 12log⁡2(2(1−s))\frac{1}{2}\log_2(2(1-s)) to its exponent.

Set q=1/(2(1−s))q=1/(2(1-s)) and replace μ\mu by qμ+(1−q)δ1q\mu+(1-q)\delta_{1}. This replaces FF by qFqF, adds exactly that cap cost, and replaces the normalizing constant by

Z′=qZ=1−σ2(1−s)≤12;Z' = qZ = \frac{1-\sigma}{2(1-s)} \le\frac{1}{2};

hence, the new measure has mean at least 1/21/2. Direct use at the kissing angle is the case s=1/2s=1/2, q=1q=1; the endpoint s=0s=0 is also covered by taking μ=δ0\mu=\delta_{0}.

We may further restrict to measures of mean 1/21/2. Indeed, for any probability measure μ\mu with mean σ≥1/2\sigma\ge1/2, put

a=12σ,μ~=aμ+(1−a)δ0.a=\frac{1}{2\sigma}, \qquad\widetilde{\mu}=a\mu+(1-a)\delta_{0}.

Then μ~\widetilde{\mu} has mean 1/21/2. Since k0(v)=1k_{0}(v)=1 and 0≤F(v)≤10\le F(v)\le1, its associated function F~\widetilde{F} satisfies

F~(v)=aF(v)+(1−a)≥F(v).\widetilde{F}(v)=aF(v)+(1-a)\ge F(v).

Thus its associated exponent satisfies Φ~≤Φ\widetilde{\Phi}\le\Phi. A lower bound for the exponent of every probability measure of mean 1/21/2 therefore also bounds the exponent of every probability measure of mean at least 1/21/2, including the measures obtained by the cap reduction above.

For such a measure, P=2FρP=2F\rho is a probability measure, and the relative-entropy identity becomes

Φ=α∗+DKL(ν∥P)2log⁡2.(11)\Phi=\alpha_{*}+\frac{D_{\mathrm{KL}}(\nu\lVert P)}{2\log2}. \tag*{(11)}

Indeed, ∫01log⁡(dν/dρ) dν=log⁡π−1\int_{0}^{1}\log(\mathrm{d}\nu/\mathrm{d}\rho)\,\mathrm{d}\nu=\log\pi-1, which gives (11) directly from (10). The logarithms are integrable: kw(v)≥(1−w)(1−v2)k_{w}(v)\ge(1-w)(1-v^{2}) gives F(v)≥(1−v2)/2F(v)\ge(1-v^{2})/2.

For 0≤v≤1/20\le v\le1/2, the function w↦kw(v)w\mapsto k_{w}(v) is concave, since

∂2kw(v)∂w2=−2(1−v2)(1−v2−3v2w2)(1−v2+v2w2)3≤0.\frac{\partial^{2}k_{w}(v)}{\partial w^{2}} = -\frac{2(1-v^{2})(1-v^{2}-3v^{2}w^{2})}{(1-v^{2}+v^{2}w^{2})^{3}} \le0.

Jensen’s inequality gives

F(v)≤k1/2(v)=3(1−v2)4−3v2(0≤v≤1/2).F(v)\le k_{1/2}(v)=\frac{3(1-v^{2})}{4-3v^{2}} \qquad(0\le v\le1/2).

For E=[0,1/2]E=[0,1/2], we have ν(E)=1/2\nu(E)=1/2 and consequently

P(E)≤3π∫01/21−v21−3v2/4 dv=23−2πarctan⁡123=β<12.P(E)\le\frac{3}{\pi}\int_{0}^{1/2}\frac{\sqrt{1-v^{2}}}{1-3v^{2}/4}\,\mathrm{d}v = \frac{2}{3}-\frac{2}{\pi}\arctan\frac{1}{2\sqrt{3}} = \beta<\frac{1}{2}.

The integral follows by setting v=sin⁡θv=\sin\theta; the strict inequality follows from the fact that tan⁡(π/12)=2−3<1/(23)\tan(\pi/12)=2-\sqrt{3}<1/(2\sqrt{3}). Applying the log-sum inequality to the partition {E,Ec}\{E,E^{c}\} then gives

DKL(ν∥P)≥−12log⁡(4P(E)(1−P(E)))≥−12log⁡(4β(1−β)).D_{\mathrm{KL}}(\nu\lVert P) \ge -\frac{1}{2}\log\left(4P(E)(1-P(E))\right) \ge -\frac{1}{2}\log\left(4\beta(1-\beta)\right).

Substitution in (11) proves the claimed bound. (It is independent of the level and cap parameter, so taking their infimum preserves it.) □\square

References

  1. [1]Nima Afkhami-Jeddi, Henry Cohn, Thomas Hartman, David de Laat, and Amirhossein Tajdini, High-dimensional sphere packing and the modular bootstrap, J. High Energy Phys. 2020 (2020), no. 12, 066, doi:10.1007/JHEP12(2020)066.
  2. [2]Kurt M. Anstreicher, Improved linear programming bounds for antipodal spherical codes, Discrete Comput. Geom. 28 (2002), no. 1, 107–114, doi:10.1007/s00454-001-0080-5.DOI
  3. [3]E. S. Barnes and G. E. Wall, Some extreme forms defined in terms of Abelian groups, J. Austral. Math. Soc. 1 (1959), no. 1, 47–63, doi:10.1017/S1446788700025064.
  4. [4]Huck Bennett, Alexander Golovnev, and Noah Stephens-Davidowitz, Difficulties constructing lattices with exponential kissing number from codes, IEEE Trans. Inform. Theory 71 (2025), no. 10, 7644–7648, doi:10.1109/TIT.2025.3593195.DOI
  5. [5]Károly Bezdek and Muhammad A. Khan, Contact numbers for sphere packings, New Trends in Intuitive Geometry, Bolyai Soc. Math. Stud., vol. 27, Springer, Berlin, Heidelberg, 2018, doi:10.1007/978-3-662-57413-3_2, pp. 25–47.DOI
  6. [6]Henry Cohn and Noam Elkies, New upper bounds on sphere packings I, Ann. of Math. (2) 157 (2003), no. 2, 689–714, doi:10.4007/annals.2003.157.689.DOI
  7. [7]Henry Cohn and Yufei Zhao, Sphere packing bounds via spherical codes, Duke Math. J. 163 (2014), no. 10, 1965–2002, doi:10.1215/00127094-2738857.DOI
  8. [8]P. Delsarte, J. M. Goethals, and J. J. Seidel, Spherical codes and designs, Geom. Dedicata 6 (1977), no. 3, 363–388, doi:10.1007/BF03187604.DOI
  9. [9]Maria Dostert, Alexander Kolpakov, and Fernando Mário de Oliveira Filho, Semidefinite programming bounds for the average kissing number, Israel J. Math. 247 (2022), no. 2, 635–659, doi:10.1007/s11856-022-2288-4.DOI
  10. [10]Maxime Fortier Bourque and Bram Petri, Kissing numbers of closed hyperbolic manifolds, Amer. J. Math. 144 (2022), no. 4, 1067–1085, doi:10.1353/ajm.2022.0023.DOI
  11. [11]D. V. Gorbachev, Extremal problem for entire functions of exponential spherical type, connected with the Levenshtein bound on the sphere packing density in Rⁿ, Izv. Tul. Gos. Univ. Ser. Mat. Mekh. Inform. 6 (2000), no. 1, 71–78 (Russian).
  12. [12]G. A. Kabatiansky and V. I. Levenshtein, On bounds for packings on a sphere and in space, Problems Inform. Transmission 14 (1978), no. 1, 1–17.
  13. [13]Thijs Laarhoven and Scott Duke Kominers, A lattice family with kissing numbers τ(Lₙ) ≥ e²√ⁿ, arXiv:2609.13608v2, 2026.
  14. [14]Gabriele Nebe, The second minimum of Barnes–Wall lattices, arXiv:2603.23133v2, 2026.arxiv.org/abs/2603.23133
  15. [15]OpenAI, Ten advances in mathematics and theoretical computer science, 2026, technical manuscript, version of August 6, 2026. https://cdn.openai.com/pdf/ten-proofs-oai.pdf (accessed September 26, 2026).
  16. [16]Alex Samorodnitsky, On linear programming bounds for spherical codes and designs, Discrete Comput. Geom. 31 (2004), no. 3, 385–394, doi:10.1007/s00454-003-2858-0.DOI
  17. [17]Naser Talebizadeh Sardari and Masoud Zargar, New upper bounds for spherical codes and packings, Math. Ann. 389 (2024), 3653–3703, doi:10.1007/s00208-023-02738-z.DOI
  18. [18]V. M. Sidel’nikov, New bounds for densest packing of spheres in n-dimensional Euclidean space, Math. USSR-Sb. 24 (1974), no. 1, 147–157, doi:10.1070/SM1974v024n01ABEH001911.DOI
  19. [19]Nils-Peter Skoruppa, Quick asymptotic upper bounds for lattice kissing numbers, Mathematika 49 (2002), no. 1–2, 51–57, doi:10.1112/S0025579300016041.DOI
  20. [20]Serge Vlăduţ, Lattices with exponentially large kissing numbers, Mosc. J. Comb. Number Theory 8 (2019), no. 2, 163–177, retracted in 2025.
  21. [21]G. L. Watson, The number of minimum points of a positive quadratic form, Dissertationes Math. (Rozprawy Mat.), no. 84, Instytut Matematyczny Polskiej Akademii Nauk, Warszawa, 1971.
  22. [22]Masoud Zargar, Stiefel manifolds and upper bounds for spherical codes and packings, arXiv:2407.10697, 2024.arxiv.org/abs/2407.10697

Paper details

Contents