# Spectral Gaps for Reversible Markov Processes with Chaotic Invariant
Measures: The Kac Process with Hard Sphere Collisions in Three Dimensions

Eric<EMAIL_ADDRESS>[ Maria
<EMAIL_ADDRESS>[ Michael Loss
<EMAIL_ADDRESS>[ Department of Mathematics, Rutgers
University,m1
Department of Mathematics and CMAF-CIO, University of Lisbon,m2
School of Mathematics, Georgia Institute of Technologym3  Address of Eric
Carlen
Department of Mathematics, Rutgers University
110 Frelinghuysen Road, Piscataway NJ 08854 USA
Address of Maria Carvalho
DEpartment of Mathematics and CMAF-CIO
University of Lisbon 1649-003 Lisbon, Portugal
Adress of Michael Loss
School of Mathematics, Georgia Institute ofTechnology
Atlanta, GA 30332 USA

###### Abstract

We develop a method for producing estimates on the spectral gaps of reversible
Markov jump processes with chaotic invariant measures, that is effective in
the case of degenerate jump rates, and we apply it to prove the Kac conjecture
for hard sphere collision in three dimensions.

60G07,

39B62,

35Q20,

spectral gap,

degenerate rates,

kinetic theory,

###### keywords:

[class=MSC]

###### keywords:

arXiv:1812.03874

, and

t1Partially supported by U.S. NSF grant DMS 1501007 and DMS 1764254
t2Partially supported by F.C.T. grant UID/MAT/04561/2013 t3Partially supported
by F.C.T. grant SFRH/BSAB/113685/2015 t4Partially supported by U.S. NSF grant
DMS-1600560

##  1 Introduction

In a seminal paper of 1956, Mark Kac [12] introduced a family of continuous
time reversible Markov jump processes on the sphere $S^{N-1}(\sqrt{N})$ of
radius $\sqrt{N}$ in ${\mathord{\mathbb{R}}}^{N}$. This family of processes,
and its generalizations, have drawn the attention of many researchers. Kac was
motivated by a connection, in the large $N$ limit, to the non-linear Boltzmann
equation. The connection arises through a particular “asymptotic independence”
property of sequences $\\{{\rm d}\mu_{N}\\}$, where ${\rm d}\mu_{N}$ is a
probability measure on $S^{N-1}$. This property is possessed, in particular,
by the sequence $\\{{\rm d}\sigma_{N}\\}$ of uniform probability measures on
$S^{N-1}(\sqrt{N})$. Let $\vec{v}=(v_{1},\dots,v_{N})$ denote a generic point
on $S^{N-1}(\sqrt{N})$ of radius $\sqrt{N}$. Let $\phi$ be any bounded
continuous function on ${\mathord{\mathbb{R}}}^{k}$ and ${\rm
d}\gamma=(2\pi)^{-1/2}e^{-v^{2}/2}{\rm d}v$ be the unit Gaussian probability
measure on ${\mathord{\mathbb{R}}}$. As is well known, going back at least to
Mehler [14],

$\lim_{N\to\infty}\int_{S^{N-1}(\sqrt{N})}\phi(v_{1},\dots,v_{k}){\rm
d}\sigma_{N}=\int_{{\mathord{\mathbb{R}}}^{k}}\phi(v_{1},\dots,v_{k}){\rm
d}\gamma^{\otimes k}\ .$

As long as one only looks at coordinates belonging to a fixed, finite set, in
the large $N$ limit, the coordinates in this set are asymptotically
independent. The main result of [12] concerned sequences of probability
measures $\\{{\rm d}\mu_{N}\\}$ on $S^{N-1}(\sqrt{N})$ with the property that,
for some probability density $f$ on ${\mathord{\mathbb{R}}}$ with zero mean
and unit variance,

$\lim_{N\to\infty}\int_{S^{N-1}(\sqrt{N})}\phi(v_{1},\dots,v_{k}){\rm
d}\mu_{N}=\int_{{\mathord{\mathbb{R}}}^{k}}\phi(v_{1},\dots,v_{k})\prod_{j=1}^{k}f(v_{j}){\rm
d}v_{j},$

in which case the sequence $\\{{\rm d}\mu_{N}\\}$ was said by Kac to be
$f(v){\rm d}v$ chaotic. He proved that chaoticity was propagated in time by
solutions of the forward Kolmogrov equations associated to the Kac processes.
Moreover, if $\\{{\rm d}\mu_{N}(t)\\}$ is the sequence of laws at time $t$
starting from an $f(v){\rm d}v$ chaotic sequence, $\\{{\rm d}\mu_{N}(t)\\}$ is
$f(t,v){\rm d}v$ chaotic where $f(t,v)$ is the solution of the Kac- Boltzmann
equation with initial data $f(v)$. (The Kac Boltzmann equation is a simple
model of the Boltzmann equation for a gas in one dimension.) He also made a
conjecture, that went unsolved for a long time, concerning the spectral gap of
the generator of this family of processes. Since the processes are reversible,
their generators are self adjoint, and it is not hard to see that the null
space is spanned by the constants. Kac conjectured a gap $\Delta_{N}$
separating $0$ from the rest of the spectrum that is bounded below uniformly
in $N$. That is, $\lim_{N\to\infty}\Delta_{N}>0$. This was finally proved by
Janvresse in 2000 [11], and shortly afterwards the exact value of $\Delta_{N}$
for all $N$ was determined in [4].

A few years after his original work, Kac returned to these problems [13], but
this time for a physically realistic model of a gas in three dimensions
undergoing “hard sphere” collisions that conserve energy and momentum. As he
showed, this physical model would have, through propagation of chaos, a direct
connection to the actual Boltzmann equation for hard sphere collisions, and
not only a toy model of it. However, in the physical model, the rates at which
different pairs of molecules collide depend on their velocities: The rates are
not bounded away from $0$, and there is no bound from above that is uniform in
$N$. It is much harder to estimate spectral gaps for the generators of jump
processes with rates that are not bounded from below, and the lack of an upper
bound that is uniform in $N$ makes it much harder to prove propagation of
chaos.

In this paper, we prove the Kac conjecture for the Kac model with hard sphere
collisions in ${\mathord{\mathbb{R}}}^{3}$. We do so using a method that has
three essential components. These are:

(1) The introduction of a conjugate process, in which at each step all but one
of the velocities are updated. The rates in this process are still not bounded
below, but they depend only on the one velocity that is left fixed during the
jump. There is also a simple connection between the spectral gaps of the
original process and the conjugate process, and the central problem becomes
the determination of the spectral gap for the conjugate process.

(2) Quantitative estimates on the chaoticity of the sequence of invariant
measures: We prove and apply estimates quantitatively expressing the near
independence of any finite set of coordinates for large $N$.

(3) A trial function decomposition: We decompose any trial function $f$ for
the spectral gap problem into 3 pieces, $f=s+g+h$ that are mutually orthogonal
in the $L^{2}$ space for the invariant measure, and, due to quantitative chaos
estimates, are nearly orthogonal with respect to the inner product given by
the Dirichlet form of the conjugate process. Each of these pieces has a
particular special structure that facilitates the proof of estimates of the
type we seek.

The first two components have been present in our work on Kac type models
since our early papers [4, 5] on the models (as in [12]) with uniform jump
rates, though in the early papers, the conjugate process is not considered
explicitly as a process. However, the connection between its spectral gap and
the spectral gap for the Kac process has been central to the approach from the
beginning. Work by two of us and Jeff Geronimo [7] dealt with the quantitative
chaos estimates needed for the three dimensional energy and momentum
conserving collision considered here, but applied them to “Maxwellian
molecules” models which, unlike to hard sphere model, has rates that are
bounded below. There too, the approach yielded the exact value of the the
spectral gap for a wide class of “Maxwellian molecules” models.

Finally in [6] we proved the Kac conjecture for a “hard sphere” model with one
dimensional velocities, and introduced a somewhat simpler version of component
(3), the trial function decomposition. In application to kinetic theory, as
explained in [6], the spectral gap in the symmetric sector, i.e., for
functions that are invariant under permutations of coordinates, is especially
important. It is this quantity that can be related to the spectral gap for the
linearized Boltzmann equation, and one would like to have explicit estimates
on this gap. Therefore, in [6] we worked hard to render all estimates as sharp
and explicit as possible, and to treat only the symmetric sector for which
fewer estimates were required.

It was clear to us at the time we wrote [6] that we had a general method that
would prove the existence of a spectral gap, uniformly in $N$, for the
physical three dimensional hard sphere Kac model, and we announced this in
several lectures. The result is quoted in reference 9 of [15], as a personal
communication, and used in the development of the quantitative treatment of
propagation of chaos that is provided there. After our paper [6] appeared with
the details provided only for the symmetric sector and the one dimensional
model, Stéphane Mischler and Clément Mouhot asked us several times to provide
the details. This paper answers their request, and moreover, in the course of
preparing this answer, it has provided a clearer picture of the how the method
explained [6] in can be extended and applied to more complicated models, such
as the main example treated here.

The method to be explained here may be applied to a wide class of sequences of
reversible Markov jump processes whose sequence of invariant measures
satisfies certain “quantitative chaos” estimates that are specified here. The
method is not at all restricted to the treatment of the symmetric sector, and
perhaps had we explained the method in [6] without obscuring it behind the
detail of so many explicit computations, necessary for the precise
quantitative estimates obtained there, this would have been clear some years
ago.

Therefore, in the present paper, we prove the Kac conjecture for hard sphere
collisions in three dimensions without any symmetry condition in as simple a
manner as possible to provide a clear view of the method. To do this, we make
use of constants $C$ that change from line to line but are independent of $N$
that are not explicitly evaluated here, but easily could be – at the expense
of more pages and less clarity.

In addition to the applications to quantitative propagation of chaos developed
in [15], uniform bounds on the spectral gap are important in certain problems
concerning the hydrodynamic limits of certain kinetic models, as explained in
[10]. These authors considered a one dimensional model essentially equivalent
to the one considered in [6], and asked for the spectral gap. Sasada [17]
provided the answer to the question they raised, noting that she could not
simply apply the result of [6] as it applied to the symmetric sector only.
This is true, but as shown here, the method used in [6] may readily extended
to answer a much broader range of questions. Much beautiful work has been done
on the question of estimating spectral gaps for Kac type processes, and we
refer to the papers of Sasada [17] and Caputo [1, 2], in addition to our own
papers cited here, for significant contributions. However, it is not clear to
us that any of the other methods that have been developed for this class of
models applies to the main example at hand which is considerably more complex
than the models considered in most other work.

###  1.1 The Kac collision process

For $N\in{\mathord{\mathbb{N}}}$, $p\in{\mathord{\mathbb{R}}}^{3}$ and
$E>|p|^{2}$, let ${\mathcal{S}}_{N,E,p}$ be the set consisting of $N$–tuples
$\vec{v}=(v_{1},\dots,v_{N})$ of vectors $v_{j}$ in
${\mathord{\mathbb{R}}}^{3}$ with
${\displaystyle\frac{1}{N}\sum_{j=1}^{N}|v_{j}|^{2}=E}$ and
${\displaystyle\frac{1}{N}\sum_{j=1}^{N}v_{j}=p}$. In what follows, a point
$\vec{v}\in{\mathcal{S}}_{N,E,p}$ specifies the velocities of a collection of
$N$ particles with mass $2$, so that $E$ is the kinetic energy per particle,
and $p$ is one-half the momentum per particle. The Markov jump process
introduced by Mark Kac [13] describes a random binary collision process for
the $N$ particles, in which the collisions conserve both energy and momentum,
and thus if the process starts on ${\mathcal{S}}_{N,E,p}$, it will remain on
${\mathcal{S}}_{N,E,p}$ for all time.

Recall that a random variable $T$ with values in $(0,\infty)$ is exponential
with parameter $\lambda$ in case ${\rm Pr}(T\geq t)=e^{-\lambda t}$. When the
collision process begins, associated to each pair $(v_{i},v_{j})$, $i<j$, is
an exponential random variable $T_{i,j}$ with parameter

$\lambda_{i,j}=N\left(\begin{matrix}N\cr
2\cr\end{matrix}\right)^{-1}|v_{i}-v_{j}|^{\alpha}\ ,$ (1)

where $0\leq\alpha\leq 2$, and $\alpha=1$ is the case of main interest: As
explained in [13], (1) is motivated by a connection between the Kac process
and the Boltzmann equation, and $\alpha=1$ corresponds to “hard-sphere
collisions”.

$T_{i,j}$ represents the waiting time for particles $i$ and $j$ to collide,
and the set of these random times is taken to be independent. The first
collision occurs at time

$T=\min_{i<j}\\{T_{i,j}\\}\ .$ (2)

As is well known [8], the minimum of an independent set of exponential random
variables is itself exponential, and the parameter of the minimum is the sum
of the parameters of the random variables in the set. In particular, if
$\alpha=0$, $T$ is exponential with parameter $N$, and the expected waiting
time for the first collision of some pair to occur is $1/N$.

At the time $T$, the pair $(i,j)$ furnishing the minimum collide: The state of
the process “jumps” from $(v_{1},\dots,v_{N})$ to
$(v_{1},v_{2},\dots,v_{i}^{*},\dots,v_{j}^{*},\dots,v_{N})$, where only
$v_{i}$ and $v_{j}$ have changed. Since the process is conceived to model
momentum and energy conserving collisions we require that

$v_{i}^{*}+v_{j}^{*}=v_{i}+v_{j}\qquad{\rm
and}\qquad|v_{i}^{*}|^{2}+|v_{j}^{*}|^{2}=|v_{i}|^{2}+|v_{j}|^{2}\ .$ (3)

Then, by the parallelogram law, it follows that

$|v_{i}^{*}-v_{j}^{*}|=|v_{i}-v_{j}|\ .$ (4)

Given $v_{i}$ and $v_{j}$, the kinematically possible collisions of particles
$i$ and $j$; i.e., those satisfying (3), may be parameterized in term of a
unit vector $\sigma\in S^{2}$, the unit sphere in ${\mathord{\mathbb{R}}}^{3}$
as follows:

$\displaystyle v_{i}^{*}(\sigma)$ $\displaystyle=$
$\displaystyle\frac{v_{i}+v_{j}}{2}+\frac{|v_{i}-v_{j}|}{2}\sigma$
$\displaystyle v_{j}^{*}(\sigma)$ $\displaystyle=$
$\displaystyle\frac{v_{i}+v_{j}}{2}-\frac{|v_{i}-v_{j}|}{2}\sigma$ (5)

The particular kinematically possible collision that occurs at time $T$ is
selected according to the following rule: There is given, in the specification
of the process, a non-negative, even function $b$ on $[-1,1]$ such that for
any fixed $\sigma^{\prime}\in S^{2}$, with ${\rm d}\sigma$ denoting the
uniform probability measure on $S^{2}$

$\int_{S^{2}}b(\sigma\cdot\sigma^{\prime}){\rm d}\sigma=1\quad{\rm
or,equivalently},\qquad\frac{1}{2}\int_{-1}^{1}b(t){\rm d}t=1\ .$ (6)

The example of main interest turns out to be

$b(x)=1\ .$ (7)

When $\alpha=1$ and $b$ is given by (7), the Kac process models “hard sphere”
or “billiard ball” collisions [13]. (There are two standard parameterizations
of the set of energy and momentum conserving collisions, the “$\sigma$
parameterization” given by (1.1), and the “$n$ parameterization”. While the
latter is often used in physics texts and is used in [13], the former, used
here, has advantages. One is that in this parameterization, $b$ is constant,
while in the other it is not due to a non-constant Jacobian relating the two
parameterizations. See Appendix A.1 of [3] for more information; equation
(A.18) of [3] is the formula relating the $b$ functions for the two
representations.)

In any case, as long as $v_{i}\neq v_{j}$,
$b(\sigma\cdot(v_{i}-v_{j})/|v_{i}-v_{j}|)$ is a probability density on
$S^{2}$. At time $T$, $\sigma$ is selected from the law
$b(\sigma\cdot(v_{i}-v_{j})/|v_{i}-v_{j}|){\rm d}\sigma$, and then the process
executes the collision step in which $v_{i}^{*}$ and $v_{j}^{*}$ are given by
(1.1). (If $v_{i}=v_{j}$, no jump is made.) Then, all of the waiting times are
“reset” and the process begins afresh. This completes the probabilistic
description of the one parameter family of Kac collision process.

This one parameter family of Kac collision process is a little more general
than the one considered by Kac: There is an extra parameter $\alpha$ that
ranges between $0$ and $2$. The case $\alpha=0$ corresponds to Maxwellian
molecules as in [12] or [7]. The case $\alpha=1$ is the hard sphere case that
is our main focus. The case $\alpha=2$ is the case of “super hard spheres” and
estimates for this case will be useful in our study of $\alpha=1$. Villani
[20] discovered in the context of entropy production estimates that analysis
of the non-physical case $\alpha=2$ could provide very helpful information on
the physical cases $\alpha\leq 1$, and we make essential use of this insight
in our analysis of spectral gaps.

###  1.2 The generator of the Kac process

The object of our investigation is the spectral gap for the generator of the
Markov semigroup associated to this process. For any continuous function $f$
on ${\mathcal{S}}_{N,E,p}$, in particular without any symmetry assumption,
define

$L_{N,\alpha}f(\vec{v})=\frac{1}{h}\lim_{h\to 0}{\rm E}\\{f(\vec{v}(h))\ |\
\vec{v}(0)=\vec{v}\ \\}\ .$

We can write this more explicitly as

$L_{N,\alpha}f(\vec{v})=-{N}{\left(\begin{matrix}N\cr
2\cr\end{matrix}\right)}^{-1}\sum_{i<j}|v_{i}-v_{j}|^{\alpha}\left[f(\vec{v})-[f]^{(i,j)}(\vec{v})\right]$
(8)

where

$[f]^{(i,j)}(\vec{v})=\int_{S^{2}}b\left(\sigma\cdot\frac{v_{i}-v_{j}}{|v_{i}-v_{j}|}\right)f(R_{i,j,\sigma}\vec{v}){\rm
d}\sigma$ (9)

and
${\displaystyle(R_{i,j,\sigma}\vec{v})_{k}=\begin{cases}v_{i}^{*}(\sigma)&k=i\\\
v_{j}^{*}(\sigma)&k=j\\\ v_{k}&k\neq i,j\end{cases}}$.

By (4) and (1.1),

$\cos\theta:=\sigma\cdot\frac{v_{i}-v_{j}}{|v_{i}-v_{j}|}=\frac{v_{i}^{*}-v_{j}^{*}}{|v_{i}^{*}-v_{j}^{*}|}\cdot\frac{v_{i}-v_{j}}{|v_{i}-v_{j}|}\
.$

By this and (4) once again, rates for the jump from $\vec{v}$ to
$R_{i,j,\sigma}\vec{v}$ and from $R_{i,j,\sigma}\vec{v}$ to $\vec{v}$ are
equal. This is the property of “detailed balance” or “microscopic
reversibility”. The analytic expression of this is self-adjointness of the
generator $L_{N,\alpha}$:

Let ${\rm d}\sigma_{N}$ denote the uniform probability measure on
${\mathcal{S}}_{N,E,p}$. (Note that ${\mathcal{S}}_{N,E,p}$ is isometric to a
sphere of radius $\sqrt{N(E-|p|^{2})}$ in ${\mathord{\mathbb{R}}}^{3N-4}$, and
by uniform, we mean uniform with respect to the symmetries of this sphere.)

For any two unit vectors $\sigma$ and $\omega$, one sees from (1.1) that

$R_{i,j,\sigma}(R_{i,j,\omega}\vec{v})=R_{i,j,\sigma}\vec{v}\ .$ (10)

From this and the fact that the measure ${\rm d}\sigma_{N}\otimes{\rm
d}\sigma$ is invariant under

$(\vec{v},\sigma)\mapsto(R_{i,j,\sigma}\vec{v},(v_{i}-v_{j})/|v_{i}-v_{j}|)\
,$

it follows that for any two continuous functions $f$ and $g$ on
${\mathcal{S}}_{N,E,p}$,

$\langle g,L_{N,\alpha}f\rangle_{L^{2}(\sigma_{N})}=\langle L_{N,\alpha}g,\
f\rangle_{L^{2}(\sigma_{N})}\ ,$

where $\langle\cdot,\cdot\rangle_{L^{2}(\sigma_{N})}$ denotes the inner
product on $L^{2}({\mathcal{S}}_{N,E,p},\sigma_{N})$. Thus, $L_{N,\alpha}$ is
a self adjoint operator on $L^{2}({\mathcal{S}}_{N,E,p},\sigma_{N})$. Notice
that the formulas (8) and (9) do not involve the parameters $E$ and $p$, and
hence our notation references only $N$ and $\alpha$.

Define the quadratic form ${\cal E}_{N,\alpha}$ by ${\cal
E}_{N,\alpha}(f,f)=-\langle f,L_{N,\alpha}f\rangle_{L^{2}(\sigma_{N})}$. A
simple computation using (10) shows that

${\cal E}_{N,\alpha}(f,f)=\\\ \frac{N}{2}{\left(\begin{matrix}N\cr
2\cr\end{matrix}\right)}^{-1}\sum_{i<j}\int_{{\mathcal{S}}_{N,E,p}}\int_{S^{2}}|v_{i}-v_{j}|^{\alpha}b\left(\sigma\cdot\frac{v_{i}-v_{j}}{|v_{i}-v_{j}|}\right)\left[f(\vec{v})-f(R_{i,j,\sigma}\vec{v})\right]^{2}{\rm
d}\sigma{\rm d}\sigma_{N}\ .$ (11)

One sees from this expression that $L_{N,\alpha}$ is a negative semi-definite
operator, and that provided $b$ is continuous at $1$, $L_{N,\alpha}f=0$ if and
only if $f$ is constant. We are interested in the spectral gap of the operator
$L_{N,\alpha}$ on $L^{2}({\mathcal{S}}_{N,E,p},\sigma_{N})$:

$\Delta_{N,\alpha}(E,p)=\inf\left\\{{\cal E}_{N,\alpha}(f,f)\ :\ \langle
f,1\rangle_{L^{2}(\sigma_{N})}=0\quad{\rm
and}\quad{\|}f{\|}_{L^{2}(\sigma_{N})}^{2}=1\ \right\\}\ .$ (12)

For fixed $N$, the dependence of $\Delta_{N,E,p}$ on $E$ and $p$ is quite
simple: Consider the point transformation

$\phi_{E,p}(v_{1},\dots,v_{N}):=\frac{1}{\sqrt{E-|p|^{2}}}(v_{1}-p,\dots,v_{N}-p)\
$

that identifies ${\mathcal{S}}_{N,E,p}$ with ${\mathcal{S}}_{N,1,0}$. The
induced transformation $U_{E,p}$ from
$L^{2}({\mathcal{S}}_{N,1,0},\sigma_{N})$ to
$L^{2}({\mathcal{S}}_{N,E,p},\sigma_{N})$ given by $U_{E,p}f=f\circ\phi_{E,p}$
is evidently unitary. A simple computation then shows that

${\cal E}_{N,\alpha}(U_{E,p}f,U_{E,p}f)=(E-|p|^{2})^{\alpha/2}{\cal
E}_{N,\alpha}(f,f,)\ .$ (13)

As an immediate consequence,

$\Delta_{N,\alpha}(E,p)=(E-|p|^{2})^{\alpha/2}\Delta_{N,\alpha}(1,0)\ .$ (14)

The dependence of $\Delta_{N,\alpha}(E,p)$ on $N$ is not so simple.
Nonetheless, we have seen that the problem of estimating the quantity
$\Delta_{N,\alpha}(E,p)$ is essentially the same as the problem of estimating
$\Delta_{N,\alpha}(1,0)$. We therefore simplify our notation:

######  DEFINITION 1.1 (Spectral gap).

The spectral gap for the $N$ particle Kac model is the quantity

$\Delta_{N,\alpha}:=\Delta_{N,\alpha}(1,0)\ .$ (15)

In what follows, we write ${\mathcal{S}}_{N}$ to denote
${\mathcal{S}}_{N,1,0}$, and consider the Kac process on ${\mathcal{S}}_{N}$
unless other values of $E$ and $p$ are explicitly specified. The Kac
conjecture for hard sphere collisions [13] is that
$\liminf_{N\to\infty}\Delta_{N,1}>0\ .$ Our main result shows somewhat more:

######  THEOREM 1.2 (Spectral gap for the Kac Model with $0\leq\alpha\leq 2$).

For each continuous non-negative even function $b$ on $[-1,1]$ statisfying
(6), and for each $\alpha\in[0,2]$, there is a strictly positive constant $K$
depending only on $b$ and $\alpha$, and explicitly computable, such that

$\Delta_{N,\alpha}\geq K>0$

for all $N$. In particular, this is true with $b$ given by (7) and $\alpha=1$,
the $3$ dimensional hard sphere Kac model.

###  1.3 The conjugate Kac process and its generator

Our method involves the introduction of another family of reversible Markov
jump processes on ${\mathcal{S}}_{N}$ that are conjugate to the Kac process.
For fixed $N$ and $\alpha$, this process is described as follows: Given
$\vec{v}\in{\mathcal{S}}_{N}$, Let
$\\{\widehat{T}_{1},\dots,\widehat{T}_{N}\\}$ be $N$ independent exponential
variables such that the parameter $\lambda_{k}(\vec{v})$ of $\widehat{T}_{k}$
is

$\lambda_{k}(\vec{v})=\frac{1}{N}\left[\frac{N^{2}-(1+|v_{k}|^{2})N}{(N-1)^{2}}\right]^{\alpha/2}\
.$

Since the total energy is $N$ and the total momentum is zero, the maximum
possible value of $|v_{k}|^{2}$ on ${\mathcal{S}}_{N}$ is $N-1$, and thus
$\lambda_{k}\geq 0$, with equality only when $|v_{k}|^{2}$ takes on its
maximal value.

The first jump time is
$\widehat{T}=\min\\{\widehat{T}_{1},\dots,\widehat{T}_{N}\\}$. At the jump
time, if $k$ is the index furnishing the minimum, $\vec{v}$ jumps to a new
point on ${\mathcal{S}}_{N}$ such that $v_{k}$ is unchanged, but conditional
on $v_{k}$, the other coordinates are redistributed uniformly. That is, the
process makes a conditional jump to uniform, conditional on $v_{k}$ which is
held fixed. After the jump, the process starts afresh. This completes the
description of the conjugate Kac process.

Note that the conjugate process is trivial for $N=2$, since then
$v_{2}=-v_{1}$, so that given one velocity, the other is known exactly, and
the “conditional jump to uniform” is no jump at all in this case. However,
alsready for $N=3$, the process is far from trivial.

######  Remark 1.3.

If $|v_{k}|^{2}$ is close to its expected value of $1$, then
$\lambda_{k}(\vec{v})\approx\frac{1}{N}$, which is exact for $\alpha=0$. In
this case, we have $N$ independent Poisson clocks with rate $\tfrac{1}{N}$
each, so that the mean waiting time for some jump is $1$.

For $\alpha>0$, the rates $\lambda_{k}(\vec{v})$ are not bounded away from
$0$. However, at most one of them can be very close to zero for any given
state $\vec{v}$. This is because, $\lambda_{k}(\vec{v})=0$ if and only if
$|v_{k}|^{2}$ takes on its maximum value, $N-1$. For at most one value of $k$
is it possible that $|v_{k}|^{2}>\tfrac{1}{2}N$, and for
$|v_{j}|^{2}\leq\tfrac{1}{2}N$,
$\lambda_{j}(\vec{v})=\frac{1}{2N}+{\mathcal{O}}(\frac{1}{N^{2}})$. Thus, for
all $\alpha\in[0,2]$, for large $N$, the expected waiting time for a jump is
very close to $1/N$, and this one jump will bring $N-1$ of the particles very
close to equilibrium. If the expected waiting time were exactly $1/N$ and the
jump took all $N$ particles to equilibrium, the spectral gap would be exactly
$1-1/N$. This is not misleading; we shall show that for the conjugate Kac
process, the spectral gap is indeed $1-1/N$ plus lower order corrections.

To write down the generator, introduce the conditional expectation operators
$P_{k}$, $k=1,\dots,N$, defined as follows:

For any function $\phi$ in $L^{2}({\mathcal{S}}_{N})$, and any $k$ with $1\leq
k\leq N$, define $P_{k}(\phi)$ to be the orthogonal projection of $\phi$ onto
the subspace of $L^{2}({\mathcal{S}}_{N})$ consisting of square integrable
functions that depend on $\vec{v}$ through $v_{k}$ alone. That is,
$P_{k}(\phi)$ is the unique element of $L^{2}({\mathcal{S}}_{N})$ of the form
$f(v_{k})$ such that

$\int_{{\mathcal{S}}_{N}}\phi(\vec{v})g(v_{k}){\rm
d}\sigma_{N}=\int_{{\mathcal{S}}_{N}}f(v_{k})g(v_{k}){\rm d}\sigma_{N}$ (16)

for all continuous functions $g$ on ${\mathord{\mathbb{R}}}^{3}$. In
probabilistic language, $P_{k}\phi$ is the conditional expectation of $\phi$
given $v_{k}$:

$P_{k}\phi=E\\{\phi\ :\ v_{k}\\}\ .$ (17)

The generator of the conjugate Kac process is then given by

$\widehat{L}_{N,\alpha}f=-\frac{1}{N}\sum_{k=1}^{N}\left[\frac{N^{2}-(1+|v_{k}|^{2})N}{(N-1)^{2}}\right]^{\alpha/2}[f-P_{k}f]\
,$ (18)

which is the analog of (8). Define the quadratic form
${\mathcal{D}}_{N,\alpha}$ by

${\mathcal{D}}_{N,\alpha}(f,f)=-\langle
f,\widehat{L}_{N,\alpha}f\rangle_{L^{2}(\sigma_{N})}\ .$

A simple computation using (10) shows that

${\cal
D}_{N,\alpha}(f,f)=\frac{1}{N}\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}\left[\frac{N^{2}-(1+|v_{k}|^{2})N}{(N-1)^{2}}\right]^{\alpha/2}\left[f^{2}-fP_{k}f\right]{\rm
d}\sigma_{N}\ .$ (19)

The spectral gap for the conjugate Kac process is the quantity defined by

$\widehat{\Delta}_{N,\alpha}=\inf\left\\{{\mathcal{D}}_{N,\alpha}(f,f)\ :\ ,\
\langle f,1\rangle_{L^{2}({\mathcal{S}}_{N})}=0\quad{\rm
and}\quad{\|}f{\|}_{L^{2}({\mathcal{S}}_{N})}^{2}=1\ \right\\}\ .$ (20)

The following theorem bears out the heuristic discussion in Remark 1.3

######  THEOREM 1.4.

For all $N\geq 3$, and all $\alpha\in[0,2]$, $\widehat{\Delta}_{N,\alpha}>0$.
Moreover, there is a constant $C$ independent of $N$ such that

$\widehat{\Delta}_{N,\alpha}\geq 1-\frac{1}{N}-\frac{C}{N^{3/2}}\ .$ (21)

######  Remark 1.5.

The constant $C$ is large enough that the first statement does not follow from
(21) which is only a meaningful bound when $N$ is large enough that the right
side is positive.

###  1.4 The link between the Kac process and its conjugate

The following theorem provides the link between the Kac process and its
conjugate:

######  THEOREM 1.6.

For all $N\geq 3$,

$\Delta_{N,\alpha}\geq\frac{N}{N-1}\Delta_{N-1,\alpha}\widehat{\Delta}_{N,\alpha}\
.$ (22)

Before proving Theorem 1.6, we recall some explicit formulas that will be
useful here and elsewhere. The proof of Theorem 1.6 uses the methods
introduced in [4, 5, 6]. The estimation of $\Delta_{N,\alpha}$ in terms of
$\Delta_{N-1,\alpha}$ is based on a parameterization of ${\mathcal{S}}_{N}$,
for $N\geq 3$, in terms of ${\mathcal{S}}_{N-1}\times B$ where $B$ is the unit
ball. For each $k=1,\dots,N$, define $\pi_{k}:{\mathcal{S}}_{N}\to B$ by

$\pi_{k}(\vec{v})=\frac{1}{\sqrt{N-1}}v_{k}\ .$ (23)

(Note that because of the constraints $\sum_{j=1}^{N}v_{j}=0$ and
$\sum_{j=1}^{N}|v_{j}|^{2}=N$, the largest value of $|v_{k}|$ on
${\mathcal{S}}_{N}$ is $\sqrt{N-1}$.)

Define a map $T_{1}:{\mathcal{S}}_{N-1}\times B\to{\mathcal{S}}_{N}$ as
follows:

$T_{1}(\vec{y},v)=\left(\sqrt{N-1}v\ ,\
\beta(v)y_{1}-{1\over\sqrt{N-1}}v,\dots,\beta(v)y_{N-1}-{1\over\sqrt{N-1}}v\right)\
,$ (24)

where

$\beta^{2}(v)=\frac{N}{N-1}(1-|v|^{2})\ .$ (25)

The subscript $1$ in $T_{1}$ indicates that the vector $v$ from $B$ went into
the first place. We likewise define $T_{2},\dots,T_{N}$ by placing this
coordinate in the corresponding position.

In the coordinates $(\vec{y},v)$ on ${\mathcal{S}}_{N}$ induced by any of the
maps $T_{k}$, one has the integral factorization formula

$\int_{{\mathcal{S}}_{N}}\phi(\vec{v}){\rm
d}\sigma_{N}=\int_{B}\left[\int_{{\mathcal{S}}_{N-1}}\phi(T_{k}(\vec{y},v)){\rm
d}\sigma_{N-1}\right]{\rm d}\nu_{N}(v)\ .$ (26)

where for all $N\geq 3$,

${\rm d}\nu_{N}(v)={|S^{3N-7}|\over|S^{3N-4}|}(1-|v|^{2})^{(3N-8)/2}{\rm d}v\
.$ (27)

Also, note that for $i\neq k,j\neq k$,

$R_{i,j,\sigma}(T_{k}(\vec{y},v))=T_{k}(R_{i,j,\sigma}(\vec{y}),v)\ .$ (28)

We now have the means to relate ${\cal E}_{N,\alpha}$ to ${\cal
E}_{N-1,\alpha}$.

For each $k=1,\dots,N$, define the conditional Dirichlet form ${\cal
E}_{N,\alpha}(f,f|v_{k})$ on $L^{2}({\mathcal{S}}_{N},\sigma_{N})$ by

${\cal E}_{N,\alpha}(f,f|v_{k})=(N-1){\left(\begin{matrix}N-1\cr
2\cr\end{matrix}\right)}^{-1}\times\\\ \sum_{i<j;i,j\neq
k}\int_{{\mathcal{S}}_{N-1}}\int_{S^{2}}|y_{i}-y_{j}|^{\alpha}b\left(\sigma\cdot\frac{y_{i}-y_{j}}{|y_{i}-y_{j}|}\right)F^{2}(\vec{v},y){\rm
d}\sigma{\rm d}\sigma_{N-1}(y)\ .$ (29)

where
$F(\vec{v},y):=\left[f(T_{k}(\pi_{k}(\vec{v}),y)-f(R_{i,j,\sigma}T_{k}(\pi_{k}(\vec{v}),y))\right]$.

As the integration on the right is only over the “slices” of
${\mathcal{S}}_{N}$ at constant values of $v_{k}$, the result is still a non-
trivial function of $v_{k}$. For each fixed $v_{k}$, the conditional Dirichlet
form is simply the $N-1$ particle Dirichlet form acting in the $\vec{y}$
variables.

Note that by (24) and (25), when $\vec{v}=T_{k}(\vec{y},v)$ and $i,j\neq k$,

$|v_{i}-v_{j}|^{2}=\beta^{2}(\pi_{k}(\vec{v}))|y_{i}-y_{j}|^{2}=\frac{N^{2}-(1+|v_{k}|^{2})N}{(N-1)^{2}}|y_{i}-y_{j}|^{2}\
.$

We define, for $v\in{\mathord{\mathbb{R}}}^{3}$, $|v|^{2}\leq N-1$,

######  DEFINITION 1.7.

$w_{N}(v):=\frac{N^{2}-(1+|v|^{2})N}{(N-1)^{2}}\ .$ (30)

We therefore have that

$|v_{i}-v_{j}|^{\alpha}b\left(\sigma\cdot\frac{v_{i}-v_{j}}{|v_{i}-v_{j}|}\right)=w^{\alpha/2}_{N}(v_{k})|y_{i}-y_{j}|^{\alpha}b\left(\sigma\cdot\frac{y_{i}-y_{j}}{|y_{i}-y_{j}|}\right)\
.$

Then, using (28), one easily checks that

${\cal
E}_{N,\alpha}(f,f)=\frac{N}{N-1}\left(\frac{1}{N}\sum_{k=1}^{N}\int_{B}w^{\alpha/2}_{N}(v_{k}){\cal
E}_{N,\alpha}(f,f|v_{k}){\rm d}\nu_{N}(v_{k}/\sqrt{N-1})\right)\ .$ (31)

###  1.5 Proof of Theorem 1.6

###### Proof of Theorem 1.6.

To estimate the right hand side of (31) in terms of $\Delta_{N-1,\alpha}$, we
must take into account that for fixed $v_{k}$, $f$ need not be orthogonal to
the constants as a function of the remaining variables $\vec{y}$. To take this
into account, we use the projection operators already introduced in (16) and
(17). Using the factorization formula (26), we have an explicit formula:

$P_{k}\phi(\vec{v})=\int_{{\mathcal{S}}_{N-1}}\phi(T_{k}(\vec{y},v_{k}/\sqrt{N-1})){\rm
d}\sigma_{N-1}\ ,$

Now note that ${\cal E}_{N,\alpha}(f,f|v_{k})={\cal
E}_{N,\alpha}(f-P_{k}f,f-P_{k}f|v_{k})$, and then using the spectral gap for
$N-1$ particles and (31), one has

$\displaystyle{\cal E}_{N,\alpha}(f,f)$ $\displaystyle\geq$
$\displaystyle\frac{N}{N-1}\Delta_{N-1,\alpha}\left(\frac{1}{N}\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}w^{\alpha/2}_{N}(v_{k})\left[f-P_{k}f\right]^{2}{\rm
d}\sigma_{N}\right)$ (32) $\displaystyle=$
$\displaystyle\frac{N}{N-1}\Delta_{N-1,\alpha}{\cal D}_{N,\alpha}(f,f)$

since

$\frac{1}{N}\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}w^{\alpha/2}_{N}(v_{k})\left[f-P_{k}f\right]^{2}{\rm
d}\sigma_{N}=\\\
\frac{1}{N}\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}w^{\alpha/2}_{N}(v_{k})\left[f^{2}-fP_{k}f\right]{\rm
d}\sigma_{N}={\cal D}_{N,\alpha}(f,f)\ .$

The theorem follows directly from (32) and the variational characterizations
of $\Delta_{N,\alpha}$ and $\widehat{\Delta}_{N,\alpha}$ ∎

###  1.6 Proof of the main theorem

Combining Theorem 1.6 and Theorem 1.4 yields, for a constant $C$, independent
of $N>3$,

$\Delta_{N,\alpha}\geq\left(1-\frac{C}{N^{3/2}}\right)\Delta_{N-1,\alpha}$
(33)

The main result will follow easily from this, and a bound on
$\Delta_{2,\alpha}$, and our next task is to prove such a bound.

Because $|v_{i}-v_{j}|$ can be arbitrarily small on ${\mathcal{S}}_{N}$ for
any $N>2$, for any given $C>0$, there will be functions $f\in
L^{2}(\sigma_{N})$ that satisfy $\langle f,1\rangle_{L^{2}(\sigma_{N})}=0$ and
${\|}f{\|}^{2}_{L^{2}(\sigma_{N})}=1$ such that
$f(\vec{v})L_{N,\alpha}f(\vec{v})<Cf(\vec{v})L_{N,0}f(\vec{v})$ for some
$\vec{v}\in{\mathcal{S}}_{N}$. This precludes a simple and direct comparison
of the Dirichlet forms ${\cal E}_{N,\alpha}$ and ${\mathcal{E}}_{N,0}$.

For $N=2$, things are much better: Then by definition of
${\mathcal{S}}_{2}:={\mathcal{S}}_{2,1,0}$, for all
$(v_{1},v_{2})\in{\mathcal{S}}_{2}$, $v_{2}=-v_{1}$, and $|v_{1}|=|v_{2}|=1$,
so that $|v_{1}-v_{2}|=2$ everywhere on ${\mathcal{S}}_{2}$. That is, for
$N=2$, there is no significant difference between $\alpha=0$ and $\alpha>0$.
For $\alpha=0$ and a number of choices of $b$, $\Delta_{2}$ has been computed
in [7]. The following is proved in Lemma 2.1 of [7]

######  LEMMA 1.8 (Spectral gap for $N=2$ and hard sphere collisions ).

With $b(x)=1$,

$\Delta_{2,1}=2\ .$ (34)

The proof given in [7] is fairly simple, and it is easy to apply the formulas
there to other choices for the probability density $b$, and to show that as
long as $b$ is even and continuous on $[-1,1]$, $\Delta_{2,\alpha}>0$ for all
$\alpha\in[0,2]$.

We are now ready to prove the main theorem:

###### Proof of Theorem 1.2.

Since $\Delta_{2,\alpha}>0$ by Lemma 1.8, Theorem 1.6 and the first part of
Theorem 1.4 yield

$\Delta_{3,\alpha}\geq\frac{3}{2}\Delta_{2,\alpha}\widehat{\Delta}_{3,\alpha}>0\
,$

and then the obvious iteration yields $\Delta_{N,\alpha}>0$ for all $N\geq 2$.
To go further and prove that $\inf_{N\geq 2}\Delta_{N,\alpha}>0$, we use the
second part of Theorem 1.4:

Let $N_{0}$ be such that ${\displaystyle 1-CN_{0}^{-3/2}>0}$. Then
${\displaystyle
K_{0}:=\prod_{j=N_{0}}^{\infty}\left(1-\frac{C}{j^{3/2}}\right)>0}$ and for
all $N\geq N_{0}$, $\Delta_{N,\alpha}\geq K_{0}\Delta_{N_{0},\alpha}$. ∎

######  Remark 1.9.

As we shall see, it is possible to explicitly compute the constant $C$ in
Theorem 1.4. To keep the presentation free of clutter, we have not carried
this through here, but it would be a simple, if tedious, exercise to track the
constants step by step. As for the first part of Theorem 1.4, it is easy to
give an explicit lower bound on $\widehat{\Delta}_{N,\alpha}$ for all $N\geq
4$, and we do so below. The case $N=3$ is more difficult, and we use a simple
compactness argument to prove $\widehat{\Delta}_{3,\alpha}>0$. However we do
sketch a method for explicitly estimating $\widehat{\Delta}_{3,\alpha}$. Thus,
the method we employ to prove Theorem 1.2 can be used to prove explicit
bounds.

It remains to prove Theorem 1.4, and we prepare the way for this in the next
section. Throughout the rest of the paper, we are concerned soley with the
conjugate Kac process. All of the analysis that directly involves the Kac
process itself is complete at this point.

##  2 Estimates for the conjugate process

It is in the proof of Theorem 1.4 that new ideas are required to deal with the
non-uniform jump rates of the conjugate process, and we begin with a heuristic
discussion of these ideas.

As in the case $\alpha=0$, we rely in part on the fact that the invariant
measure $\sigma_{N}$ (of both processes) is chaotic in the sense of Kac. More
specifically, it is $\gamma$ chaotic where

${\rm d}\gamma=(2\pi/3)^{-3/2}e^{-3|v|^{2}/2}{\rm d}v$

is the isotropic Gaussian distribution on ${\mathord{\mathbb{R}}}^{3}$ with
unit variance. This means that for any $k\in{\mathord{\mathbb{N}}}$ and any
bounded continuous function $\psi(v_{1},\dots,v_{k})$ on
${\mathord{\mathbb{R}}}^{3k}$,

$\lim_{N\to\infty}\int_{{\mathcal{S}}_{N}}\psi(v_{1},\dots,v_{k}){\rm
d}\sigma_{N}=\int_{{\mathord{\mathbb{R}}}^{3k}}\psi(v_{1},\dots,v_{k}){\rm
d}\gamma^{\otimes k}\ .$

That is, as long as $k$ is much less than $N$, the random variables
$v_{1},\dots,v_{k}$ are nearly independent, and by symmetry this is true of
any set of $k$ distinct coordinate functions on ${\mathcal{S}}_{N}$. The
notion of chaos was also introduced by Kac in [12], and the main result of
that paper was that for the model with one dimensional velocities and
$\alpha=0$, chaos is propagated by the dynamics. Propagation of chaos for
$\alpha>0$ is much harder, and this was only proved later by Sznitman [19],
also in $3$ dimensions.

In case $\alpha=0$, the range of $I-\widehat{L}_{N,0}$ has a special structure
that facilitates the study of the spectral gap for $\widehat{L}_{N,0}$. The
subspace of the range that is orthogonal to the constants consists of
functions $f$ of the form: $f(\vec{v})=\sum_{j=1}^{N}\varphi_{j}(v_{j})$ such
that each $\varphi_{j}(v_{j})$ is square integrable and such that $f$ is
orthognal to the constants. One choice for the $\varphi_{j}$’s is
$\varphi_{j}=P_{j}f$, but there are other choices: Since for any fixed
$a\in{\mathord{\mathbb{R}}}^{3}$ and $b\in{\mathord{\mathbb{R}}}$,

$\sum_{j=1}^{N}(a\cdot v_{j}+b(|v_{j}|^{2}-1))=0\,,$ (35)

we may make the replacement
$\varphi_{j}(v_{j})\longrightarrow\varphi_{j}(v_{j})+a\cdot
v_{j}+b(|v_{j}|^{2}-1)$ without changing $f(\vec{v})$. There is however, a
prefered choice of the functions $\varphi_{j}$ that plays an important role in
what follows. As we shall show, there is a unique choice that minimizes
$\sum_{j=1}^{N}{\|}\varphi_{j}(v_{j}){\|}_{2}^{2}$ which has a number of
useful properties.

######  DEFINITION 2.1.

Let ${\mathcal{A}}_{N}$ denote the subspace of $L^{2}({\mathcal{S}}_{N})$ that
is the closure of the span of functions of the form

$f(\vec{v})=\sum_{j=1}^{N}\varphi_{j}(v_{j})$ (36)

for bounded continuous functions $\varphi_{1},\dots,\varphi_{N}$ in
${\mathord{\mathbb{R}}}^{3}$ such that $\int_{{\mathcal{S}}_{N}}f{\rm
d}\sigma_{N}=0$.

When $\alpha\neq 0$, ${\mathcal{A}}_{N}$ is not an invariant subspace of
$\widehat{L}_{N,\alpha}$. Nonetheless, as we explain, the gap may be bounded
using a trial function decomposition based on ${\mathcal{A}}_{N}$, and for
this the approximate independence that comes along with the chaoticity of
$\sigma_{N}$ is essential.

To see how this works, suppose that one replaces the state space
${\mathcal{S}}_{N}$ with ${\mathord{\mathbb{R}}}^{3N}$, and replaces the
conjugate Kac process with the “conditional jump to uniform” process with
respect to ${\rm d}\gamma^{\otimes N}$. In this case, with the invariant
measure being a product measure, the corresponding conditional expectation
operators $P_{k}$ will all commute. One might therefore expect that the
operators $P_{k}$ figuring in the definition (18) of $\widehat{L}_{N,\alpha}$
almost commute for large $N$. Suppose that they exactly commute, or, what is
the same thing, that the coordinate functions $v_{1},\dots,v_{N}$ are exactly
independent.

Since $0=\int_{{\mathcal{S}}_{N}}f{\rm
d}\sigma_{N}=\sum_{j=1}^{N}\int_{{\mathcal{S}}_{N}}\varphi_{j}(v_{j}){\rm
d}\sigma_{N}=0$, replacing $\varphi_{j}(v_{j})$ by
$\varphi_{j}(v_{j})-\int_{{\mathcal{S}}_{N}}\varphi_{j}(v_{j}){\rm
d}\sigma_{N}$, we may assume without loss of generality in (36) that
$\int_{{\mathcal{S}}_{N}}\varphi_{j}(v_{j}){\rm d}\sigma_{N}=0$ for each $j$.
Granted the exact independence, we would then have that for $k\neq j$,
$P_{k}\varphi_{j}(v_{k})=0$, while
$P_{k}\varphi_{k}(v_{k})=\varphi_{k}(v_{k})$. Thus, for
$f\in{\mathcal{A}}_{N}$, $f-P_{k}f=\sum_{j\neq k}\varphi_{j}(v_{j})$, and
then, again using the independence,

$\displaystyle{\cal D}_{N,\alpha}(f,f)$ $\displaystyle=$
$\displaystyle\frac{1}{N}\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}w_{N}^{\alpha/2}(v_{k})\sum_{j\neq
k}\varphi_{j}^{2}(v_{j}){\rm d}\sigma_{N}$ (37) $\displaystyle=$
$\displaystyle\frac{1}{N}\sum_{k=1}^{N}\left(\int_{{\mathcal{S}}_{N}}w_{N}^{\alpha/2}(v_{k}){\rm
d}\sigma_{N}\right)\left(\int_{{\mathcal{S}}_{N}}\sum_{j\neq
k}\varphi_{j}^{2}(v_{j}){\rm d}\sigma_{N}\right)$

As is shown below, the integral over the rate, which is evidently independent
of $k$, is bounded below by $1-C/N^{2}$ for some constant $C$ that is
independent of $N$. Thus, we would have

$\displaystyle{\cal D}_{N,\alpha}(f,f)$ $\displaystyle\geq$
$\displaystyle\left(1-\frac{C}{N^{2}}\right)\frac{N-1}{N}\sum_{j=1}^{N}{\|}\varphi_{j}(v_{j}){\|}_{2}^{2}$
(38) $\displaystyle=$
$\displaystyle\left(1-\frac{C}{N^{2}}\right)\frac{N-1}{N}{\|}f{\|}_{2}^{2}\ $

which is even better than (21).

The equality in (38) comes from the identity
$\sum_{j=1}^{N}{\|}\varphi_{j}(v_{j}){\|}_{2}^{2}={\|}f{\|}_{2}^{2}$ which is
true when there is exact independence of the coordinate functions. In our
setting, we do not have exact independence, and must prove and use appropriate
quantitative chaos estimates. For instance, in (50) of Theorem 2.5, it is
shown that in our setting, for a particular decomposition
$f(\vec{v})=\sum_{j=1}^{N}\varphi(v_{j})$ – such decompositions are not
unique, even if one requires each $\varphi_{j}$ to be orthogonal to the
constants – one has

$\sum_{j=1}^{N}{\|}\varphi_{j}(v_{j}){\|}_{2}^{2}\leq\left(1+\frac{C}{N^{2}}\right){\|}f{\|}_{2}^{2}$
(39)

for all $N\geq 3$ and with $C$ independent of $N$. Using this after the first
inequality in (38) still yeilds someting even better than (21). .

Of course, one must consider trial functions that are not in
${\mathcal{A}}_{N}$, and for trial functions $f$ that are in
${\mathcal{A}}_{N}^{\perp}$, things are better still. Such functions are shown
to belong to the null space of $P_{k}$ for each $k$. Therefore, for
$f\in{\mathcal{A}}_{N}^{\perp}$, we would have from (19)

${\cal
D}_{N,\alpha}(f,f)=\int_{{\mathcal{S}}_{N}}\left(\frac{1}{N}\sum_{k=1}^{N}w_{N}^{\alpha/2}(v_{k})\right)f^{2}{\rm
d}\sigma_{N}\ ,$ (40)

which is a significant simplification of (19). It is shown below (see Lemma
2.12 and Remark 2.13) that for some constant $C$ independent of $N$,

$\frac{1}{N}\sum_{k=1}^{N}w_{N}^{\alpha/2}(v_{k})\geq
1-\left(1-\frac{\alpha}{2}\right)\frac{1}{N}-\frac{C}{N^{2}}\ .$

Combining this with (40) would then yield

${\cal
D}_{N,\alpha}(f,f)\geq\left(1-\left(1-\frac{\alpha}{2}\right)\frac{1}{N}-\frac{C}{N^{2}}\right){\|}f{\|}_{2}^{2}\
.$

For $\alpha>0$, this is much stronger than (21), and for this bound we do not
even use the approximate independence.

Since ${\mathcal{A}}_{N}$ is not an invariant subspace for
$\widehat{L}_{N,\alpha}$, one has to show that for $g\in{\mathcal{A}}_{N}$ and
$h\in{\mathcal{A}}_{N}^{\perp}$, $D(g,h)$ is small. We shall show, again using
the approximate independence, that

$|D(g,h)|\leq\frac{C}{N^{3/2}}{\|}g{\|}_{2}{\|}h{\|}_{2}\ .$

It is the estimate in this step that is responsible for the $N^{3/2}$ term in
(33). A more refined argument, like the one provided for this step in [6] for
the model with one dimensional velocities, would presumably improve $N^{3/2}$
to $N^{2}$, but since we have elected not to keep track of constants, there is
no point in pursuing this here.

Our proof will closely follow these heuristics, but of course we must
carefully control the departures from exact independence wherever it was used
above. There is one significant twist. Though we have the estimate (39), it is
much easier to prove the weaker analog of it with $N^{2}$ replaced by $N$, and
perhaps there are other models that are “less chaotic” in which the weaker
bound is all that one has. The weaker bound cannot be used directly in (38) to
obtain anything useful, but the simple device of defining
${\mathcal{F}}_{N,\alpha}:={\|}f{\|}_{2}^{2}-{\cal D}_{N,\alpha}(f,f)$ reduces
the problem of estimating the gap for ${\cal D}_{N,\alpha}$ to that of
obtaining an appropriate of an upper bound for ${\mathcal{F}}_{N,\alpha}$, and
with the $1$ out of the way, the weaker version of (39) becomes useful. This
is what we do in Section 3.3 to complete the proof.

The next subsection presents the “quantitative chaos” estimates that are used
to control the weak dependence of the coordinate function for large $N$. It is
important that some of the results turn out to be meaningful even for small
$N$, such as $N=3$.

###  2.1 Quantitative chaos

A number of the quantitative chaos bounds that we need may be expressed in
terms of the correlation operator $K$-operator that we now define:

Let $B$ denote the unit ball in ${\mathord{\mathbb{R}}}^{3}$. Let $N\geq 3$
and let $\nu_{N}$ be given by (27), so that for any function $\psi$ on $B$,
and any $k$,

$\int_{B}\psi(v){\rm
d}\nu_{N}=\int_{{\mathcal{S}}_{N}}\psi(\pi_{k}(\vec{v})){\rm d}\sigma_{N}\ ,$

where $\pi_{k}$ is given by (23), and $\nu_{N}$ is given by (27). We define
the operator $K$ on $L^{2}(B,\nu_{N})$ by

$\langle\psi_{1},K\psi_{2}\rangle_{L^{2}(B,\nu_{N})}=\int_{{\mathcal{S}}_{N}}\psi_{1}^{*}(\pi_{1}(\vec{v}))\psi_{2}(\pi_{2}(\vec{v})){\rm
d}\sigma_{N}\ .$ (41)

$K$ is evidently self adjoint.

######  DEFINITION 2.2.

For $j=0,\dots,4$ define fnctions $\xi_{j}(v)$ on $B$ by

$\xi_{0}(v)=1\quad\xi_{j}(v)=v_{j}\ ,j=1,2,3,\quad{\rm
and}\quad\xi_{4}(v)=(|v|^{2}-1)/(N-1)\ .$ (42)

The spectrum of $K$ is determined in [7], where the following facts are
proven:

######  LEMMA 2.3.

Let $N\geq 3$. The operator $K$ is compact. The function $\xi_{0}$ is an
eigenfunction of $K$ with eigenvalue $1$, and it spans the corresponding
eigenspace. The functions $\xi_{j}$, $j=1,2,3,4$ are eigenfunctions of $K$
with eigenvalue $-1/(N-1)$, and they are an orthogonal basis for this
eigenspace. No other eigenvalues of $K$ are larger in absolute value than
$\tfrac{5N-3}{3(N-1)^{3}}$. Therefore, for all $\psi_{1},\psi_{2}\in
L^{2}(B,\nu_{N})$ that are orthogonal to the constants, the three components
of $v$ and $v^{2}$,

$\left|\int_{{\mathcal{S}}_{N}}\psi_{1}^{*}(\pi_{1}(\vec{v}))\psi_{2}(\pi_{2}(\vec{v})){\rm
d}\sigma_{N}\right|\leq\frac{5N-3}{3(N-1)^{3}}{\|}\psi_{2}\circ\pi_{1}{\|}_{2}{\|}\psi_{2}\circ\pi_{2}{\|}_{2}\
.$ (43)

Equivalently, for all functions $\psi\in L^{2}(B,\nu_{N})$, that are
orthogonal to $1$, the three components of $v$ and $v^{2}$,

${\|}K\psi{\|}_{2}\leq\frac{5N-3}{3(N-1)^{3}}{\|}\psi{\|}_{2}\ .$ (44)

Finally, every eigenvalues $\kappa$ of $K$, other than $1$,
$\frac{5N-3}{3(N-1)^{3}}$ and $\frac{1}{N-1}$ staisfies

$-\frac{7N-3}{3(N-1)^{4}}\leq\kappa<\frac{5N-3}{3(N-1)^{3}}\ .$ (45)

######  Remark 2.4.

The number on the left in (45) is the eigenvalue denoted by $\kappa_{1,2}$ in
Section 8 of [7].

Fix some $k$, and let $\mathcal{H}$ the subspace of $L^{2}({\mathcal{S}}_{N})$
spanned by functions of the form $\varphi(v_{k})$ for some $k$. Since $v_{k}$
ranges over the ball of radius $\sqrt{N-1}$ in ${\mathord{\mathbb{R}}}^{3}$,
one may think of $\mathcal{H}$ as a Hilbert space consisting of square
integrable functions on this ball, with respect to a scaled version of the
measure $\nu_{N}$.

It will be convenient in what follows to think of $K$ as an operator on
$\mathcal{H}$. Note that $\varphi(v_{k})=\tilde{\varphi}(\pi_{k}(\vec{v}))$
where $\tilde{\varphi}(v)=\varphi(\sqrt{N-1}v)$. Define

$K\varphi(v_{k})=(K\tilde{\varphi})(\pi_{k}(\vec{v}))\ .$ (46)

The spectrum of $K$, including multiplictiy, thought of this way is naturally
the same, but the eigenfunctions change by scaling. For example, now
$|v|^{2}-1$ is an eigenfunction with eigenvalue $-1/(N-1)$. In this notation,
we have that for any function $\xi$ on ${\mathord{\mathbb{R}}}^{3}$ so that
$\xi(v_{1})$ is in $L^{2}(\sigma_{N})$,

$E\\{\xi(v_{1})\ |\ v_{2}=v\\}=K\xi(v_{N})\ .$ (47)

The $K$ operator defined by (47) is simply a “scaled” version of the $K$
operator defined in (41), scaled so it operates on functions on $\mathcal{H}$.
For some computations, particularly in the computation of eigenvalues of $K$,
the definition (41) is more convenient. For other computations, more directly
connected the the Kac process, (47) has advantages. This slight abuse of
notation will simplify many formulas that follow without introducing any
ambiguity.

Since $K$ is compact, there is a orthonormal basis of $\mathcal{H}$ consisting
of eigenvectors of $K$. This orthonormal basis is determined explicitly in
[7], but all we need to know is that is can be written as
$\\{\eta_{\iota}\\}_{\iota\geq 0}$ where

$\eta_{0}(v)=1\ ,\ \eta_{j}(v)=\sqrt{3}{\bf e_{j}}\cdot v\ ,\ 1\leq j\leq 3\
{\rm and}\ \eta_{4}(v)=C_{N}(|v|^{2}-1)\ ,$ (48)

with $C_{N}$ being a normalization constant. This follows directly from Lemma
2.3 and (46).

Let $\kappa_{\iota}$ denote the eigenvalue corresponding to $\eta_{\iota}$, so
that $K\eta_{\iota}=\kappa_{\iota}\eta_{\iota}$. Our first application of
Lemma 2.3 concerns the norm of functions in ${\mathcal{A}}_{N}$:

######  THEOREM 2.5.

Let $N\geq 3$, and let $f\in{\mathcal{A}}_{N}$ be orthogonal to $1$. Then
there is a unique choice of $\varphi_{1},\dots,\varphi_{N}$ with
$f=\sum_{j=1}^{N}\varphi_{j}(v_{j})$ and each $\varphi_{j}(v_{j})$ orthogonal
to the constants that minimizes $\sum_{k=1}^{N}{\|}\varphi_{k}{\|}_{2}^{2}$
where ${\|}\varphi_{k}{\|}_{2}^{2}$ denotes
$\int_{{\mathcal{S}}_{N}}|\varphi_{k}(v_{k})|^{2}{\rm d}\sigma_{N}$. Let

$\varphi_{j}(v_{j})=\sum_{i=1}^{\infty}a_{j,i}\eta_{i}(v_{j})$ (49)

be the expansion of $\varphi_{j}$ in the orthonormal basis consisting of
eigenfunctions of $K$ that is specified above. Then this minimizer is
characterized by

$\sum_{j=1}^{N}a_{j,i}=0\qquad{\rm for}\qquad 1\leq i\leq 4\ .$ (50)

For this choice,

$\left(1-\frac{7N-3}{3(N-1)^{3}}\right)\sum_{k=1}^{N}{\|}\varphi_{k}{\|}_{2}^{2}\leq{\|}f{\|}_{2}^{2}\leq\left(1+\frac{5N-3}{3(N-1)^{2}}\right)\sum_{k=1}^{N}{\|}\varphi_{k}{\|}_{2}^{2}\
,$ (51)

In particular, let ${\mathcal{H}}_{N,k}$ denote the subspace of
$L^{2}({\mathcal{S}}_{N})$ consisting of functions of the form
$\varphi(v_{k})$. Define ${\mathcal{B}}_{N}$ to be the subspace of
$\bigoplus_{k=1}^{N}{\mathcal{H}}_{N,k}$ consisiting of
$(\varphi_{1},\dots,\varphi_{N})$ such that (50) is satisfied. Then the
operator $T:{\mathcal{B}}_{N}\to{\mathcal{A}}_{N}$ defined by by

$T(\varphi_{1}(v_{1}),\dots,\varphi_{N}(v_{N}))=\sum_{k=1}^{N}\varphi_{k}(v_{k})$

is bounded with a bounded inverse.

######  Remark 2.6.

Define $c_{N}:=\left(1-\frac{7N-3}{3(N-1)^{3}}\right)$ and
$C_{N}=\left(1+\frac{5N-3}{3(N-1)^{2}}\right)$. Note the different exponents
in the denominator, and that $c_{N}>0$ for all $N\geq 3$. Also note that
$c_{N}=1-{\mathcal{O}}(1/N)$, and $C_{N}=1+{\mathcal{O}}(1/N)$, and then we
can rewrite (51) as

$c_{N}\sum_{k=1}^{N}{\|}\varphi_{k}{\|}_{2}^{2}\leq{\|}f{\|}_{2}^{2}\leq
C_{N}\sum_{k=1}^{N}{\|}\varphi_{k}{\|}_{2}^{2}\ ,$

and of course we would have equality here with $c_{N}=C_{N}=1$ if the
coordinate functions were exactly independent. Theorem 2.5 gives a
quantitative expression of the fact that for large $N$, the coordinate
functions are approximately pairwise independent.

###### Proof of Theorem 2.5.

As noted above, we may assume that each $\varphi_{j}$ is orthogonal to the
constants. We expand each $\varphi_{j}$ in the eigenbasis of $K$ as follows:

$\varphi_{j}(v_{j})=\sum_{\iota=1}^{\infty}a_{j,\iota}\eta_{\iota}(v_{j})\ .$
(52)

Then evidently
${\displaystyle{\|}\varphi_{j}{\|}_{2}^{2}=\sum_{\iota=1}^{\infty}|a_{j,\iota}|^{2}}$.
On account of (35), for $j=1,2,3,4$, we may replace $a_{j,i}$ by
$a_{j,i}-t_{i}$ without changing
$f(\vec{v})=\sum_{j=1}^{N}\varphi_{j}(v_{j})$. With this modification,

$\sum_{j=1}^{N}{\|}\varphi_{j}{\|}_{2}^{2}=\sum_{j=1}^{N}\sum_{i=1}^{4}(a_{j,i}-t_{i})^{2}+\sum_{j=1}^{N}\sum_{i=5}^{\infty}|a_{i,j}|^{2}\
,$

which is evidently minimized by taking
$t_{i}=-\frac{1}{N}\sum_{j=1}^{N}a_{j,i}$. and then making this replacement,
(50) is satisfied.

Next, for $j\neq k$,

$\int_{{\mathcal{S}}_{N}}\varphi^{*}_{j}(v_{j})\varphi_{k}(v_{k}){\rm
d}\sigma_{N}=\sum_{\iota,\iota^{\prime}=1}^{\infty}a^{*}_{j,\iota}a_{k,\iota^{\prime}}\langle\eta_{j,\iota},K\eta_{k,\iota^{\prime}}\rangle_{\mathcal{H}}=\sum_{\iota=1}^{\infty}\kappa_{\iota}a^{*}_{j,\iota}a_{k,\iota}\
.$

Therefore, when $f$ is given by (36) and (52) with (50) satisfied,

${\|}f{\|}_{2}^{2}=\sum_{\iota=1}^{\infty}\left(\sum_{j=1}^{N}|a_{j,\iota}|^{2}+\kappa_{\iota}\sum_{j\neq
k,j,k=1}^{N}\Re a^{*}_{j,\iota}a_{k,\iota}\right)$ (53)

For $\iota=1,\dots,4$, we have, using (50), the identity

$\sum_{j\neq k,j,k=1}^{N}\Re
a^{*}_{j,\iota}a_{k,\iota}=\left|\sum_{j=1}^{N}a_{j,\iota}\right|^{2}-\sum_{j=1}^{N}a_{i,j}^{2}=-\sum_{j=1}^{N}a_{i,j}^{2}\
,$

and then since $\kappa_{i}=-\frac{1}{N-1}$ for $i=1,\dots,4$,

$\sum_{\iota=1}^{4}\left(\sum_{j=1}^{N}|a_{j,\iota}|^{2}+\kappa_{\iota}\sum_{j\neq
k,j,k=1}^{N}\Re
a^{*}_{j,\iota}a_{k,\iota}\right)=\frac{N}{N-1}\sum_{\iota=1}^{4}\sum_{j=1}^{N}|a_{j,\iota}|^{2}\
.$ (54)

For $\iota>4$, we simply use the fact for such $\iota$, $\kappa_{\iota}$ is
${\cal O}(1/N^{2})$ or smaller, and this takes the place of (50), which is not
satisfied for such $\iota$, in eliminating a factor of $N$. Then since the
$N\times N$ matrix that has $0$ in every diagonal entry, and $1$ elsewhere has
eigenvalues $N-1$ and $-1$,

$-\sum_{j=1}^{N}a_{j,i}^{2}\leq\sum_{j\neq k,j,k=1}^{N}\Re
a^{*}_{j,\iota}a_{k,\iota}\leq(N-1)\sum_{j=1}^{N}a_{j,i}^{2}\ .$ (55)

Hence, for $\iota>4$, an upper bound on
${\displaystyle\left(\sum_{j=1}^{N}|a_{j,\iota}|^{2}+\kappa_{\iota}\sum_{j\neq
k,j,k=1}^{N}\Re a^{*}_{j,\iota}a_{k,\iota}\right)}$ is

$\left(1+\max\left\\{\frac{7N-3}{3(N-1)^{4}}\ ,\
(N-1)\frac{5N-3}{3(N-1)^{3}}\right\\}\right)\sum_{j=1}^{N}a_{j,i}^{2}$

where we have used (45). Evidently the maximum is furnished by the second
quantitiy in the braces. Summing on $\iota>4$ and combining this with (54)
yields the upper bound in (51).

For $\iota>4$, a lower bound on
${\displaystyle\left(\sum_{j=1}^{N}|a_{j,\iota}|^{2}+\kappa_{\iota}\sum_{j\neq
k,j,k=1}^{N}\Re a^{*}_{j,\iota}a_{k,\iota}\right)}$ is

$\left(1-\max\left\\{(N-1)\frac{7N-3}{3(N-1)^{4}}\ ,\
\frac{5N-3}{3(N-1)^{3}}\right\\}\right)\sum_{j=1}^{N}a_{j,i}^{2}$

where we have again used (45). Evidently the maximum is furnished by the first
quantitiy in the braces. Summing on $\iota>4$ and combining this with (54)
yields the lower bound in (51).

By Lemma 2.3, $\min\\{1-\kappa_{\iota}\\}=\tfrac{5N-3}{3(N-1)^{3}}$, and thus
we have the lower bound. For the upper bound, we use $\kappa_{\iota}=-1/(N-1)$
for $\iota=1,2,3,4$, and note that for all $N\geq 3$,
$\frac{N}{N-1}\leq\left(1+\frac{5N-3}{3(N-1)^{2}}\right)$. This gives us the
upper bound. The rest is now clear, including the fact that (35) is the only
source of non-uniqueness in the representation of $f\in{\mathcal{A}}_{N}$ in
the form $f(\vec{v})=\sum_{j=1}^{N}\varphi_{j}(v_{j})$. ∎

There is another type of quantitative chaos estimate that we need. For any
functions $\xi$ on ${\mathord{\mathbb{R}}}^{3}$ such that $\xi(v_{k})\in
L^{2}(\sigma_{N})$ for some (and hence all) $k$, consider the conditional
expectation

$E\\{\xi(v_{k})\ |\ v_{j}=v\\}=K\xi(v)$ (56)

for $j\neq k$. If the coordinate functions were exactly independent, this
would simply be the expectation of $\xi(v_{k})$, which is a finite constant.
It turns out that when $\xi(v_{k})$ is a polynomial in $|v_{k}|^{2}$, the
conditional expectation is at least bounded – not only on ${\mathcal{S}}_{N}$,
which is trivial, but the bound is independent of $N$. Here is one such
estimate:

######  LEMMA 2.7.

For $\psi(v)=|v|^{8}$, there is a constant $C<\infty$ such that
${\|}K\psi{\|}_{\infty}\leq C$ for all $N$.

###### Proof of Lemma 2.7.

The formula (24) gives us

$\displaystyle K\psi(v)$ $\displaystyle=$
$\displaystyle\int_{{\mathcal{S}}_{N-1}}\left|\sqrt{\frac{N-v^{2}}{N-1}}\vec{y}-\frac{1}{\sqrt{N(N-1)}}v\right|^{8}{\rm
d}\sigma_{N-1}$ $\displaystyle\leq$ $\displaystyle
2^{7}\left(\left(\frac{N-v^{2}}{N-1}\right)^{4}\int_{{\mathcal{S}}_{N-1}}|\vec{y}|^{8}{\rm
d}\sigma_{N-1}-\frac{|v|^{8}}{N^{4}(N-1)^{4}}\right)$

It is evident that $\int_{{\mathcal{S}}_{N-1}}|\vec{y}|^{8}{\rm
d}\sigma_{N-1}$ is bounded uniformly in $N$, and in fact,

$\lim_{N\to\infty}\int_{{\mathcal{S}}_{N-1}}|\vec{y}|^{8}{\rm
d}\sigma_{N-1}=(2\pi/3)^{-3/2}\int_{{\mathord{\mathbb{R}}}^{3}}|y|^{8}e^{-3|y|^{2}/2}\
.$

∎

In the remainder of this section we collect the other estimates of this type
that we need. Their proofs, which are more intricate but still largely
computational, are presented in Appendix A.

######  LEMMA 2.8.

There is a finite constant $C$ such that for all $N>3$ and all $v$ such that
$v=v_{N}$ for some $\vec{v}\in{\mathcal{S}}_{N}$,

$|E\\{|v_{1}|^{4}\ |\ v_{N}=v\\}-S(v)|\leq\frac{C}{N}$ (57)

where

$S(v)=\frac{N^{2}+|v|^{4}-2N|v|^{2}}{(N-1)^{2}}\ .$ (58)

######  LEMMA 2.9.

There is a finite constant $C$ such that for all $N>3$ and all $(v,w)$ such
that $(v,w)=(v_{N-1},v_{N})$ for some $\vec{v}\in{\mathcal{S}}_{N}$,

$|E\\{|v_{1}|^{4}\ |\ (v_{N-1},v_{N})=(v,w)\\}-S(v,w)|\leq\frac{C}{N}$ (59)

where

$S(v,w)=\frac{N^{2}+|v|^{4}+|w|^{4}+2N|v|^{2}+2N|w|^{2}+2|v|^{2}|w|^{2}}{(N-2)^{2}}\
.$ (60)

###  2.2 The operators $W^{(\alpha)}$ and $P^{(\alpha)}$

Let $\alpha\in[0,2]$, and define the self adjoint operator $P^{({\alpha})}$ by

$P^{({\alpha})}=\frac{1}{N}\sum_{k=1}^{N}w_{N}^{\alpha/2}(v_{k})P_{k}\ ,$ (61)

recalling that $w_{N}(v)$ is defined in (30), and $P_{k}$ is defined in (16),
or equivalently (17). For each $k$, both $P_{k}$ and the multiplication
operator ${\displaystyle w_{N}^{\alpha/2}(v_{k})}$ are commuting and self
adjoint, and hence $P^{({\alpha})}$ is indeed self adjoint, and non-negative.
Since each $P_{k}$ is a projection,

$\frac{1}{N}\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}w_{N}^{\alpha/2}(v_{k})|P_{k}f|^{2}{\rm
d}\sigma_{N}=\langle
f,P^{({\alpha})}f\rangle_{L^{2}({\mathcal{S}}_{N},\sigma_{N})}\ .$ (62)

Define the function $W^{({\alpha})}$ by

$W^{({\alpha})}=\frac{1}{N}\sum_{k=1}^{N}w_{N}^{\alpha/2}(v_{k})\ .$ (63)

Then

$\frac{1}{N}\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}w_{N}^{\alpha/2}(v_{k})f^{2}{\rm
d}\sigma=\int_{{\mathcal{S}}_{N}}W^{({\alpha})}f^{2}{\rm d}\sigma_{N}\ ,$ (64)

and we can write:

${\cal D}_{N,\alpha}(f,f):=\int_{{\mathcal{S}}_{N}}W^{({\alpha})}f^{2}{\rm
d}\sigma_{N}-\langle
f,P^{({\alpha})}f\rangle_{L^{2}({\mathcal{S}}_{N},\sigma_{N})}\ .$ (65)

Equivalently, by the computations just below (32),

${\cal
D}_{N,\alpha}(f,f)=\frac{1}{N}\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}w_{N}^{\alpha/2}(v_{k})\left[f-P_{k}f\right]^{2}{\rm
d}\sigma_{N}\ ,$ (66)

and hence ${\cal D}_{N,\alpha}(f,f)\geq 0$ for all $f$ since for each $k$,
$|v_{k}|^{2}\leq N-1$. (Recall that because of the momentum constrain, not all
of the energy can reside in a single particle.) It follows that ${\cal
D}_{N,\alpha}(f,f)=0$ if and only if $f-P_{k}f=0$ almost everywhere for each
$k$, and then in this case

${\|}f{\|}_{2}^{2}-\langle
f,P^{(0)}f\rangle_{L^{2}({\mathcal{S}}_{N})}={\mathcal{D}}_{N,0}(f,f)=\frac{1}{N}\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}|f-P_{k}f|^{2}{\rm
d}\sigma_{N}=0\ .$

Evidently, $P^{(0)}$ is a contraction, and $1$ is an eigenvalue of
multiplicity one, and the eigenspace is spanned by the constant function $1$
[7]. This proves:

######  LEMMA 2.10.

For all $N\geq 2$ and all $\alpha\in[0,2]$, and all non-zero $f\in
L^{2}({\mathcal{S}}_{N})$ that are orthogonal to the constants, ${\cal
D}_{N,\alpha}(f,f)>0$.

We use the following lemma proved in [6, Lemma 3.5]:

######  LEMMA 2.11.

For all $0<\alpha\leq 2$ and all $x>-1$,

$(1+x)^{\alpha/2}\geq 1+\tfrac{\alpha}{2}x-(1-\tfrac{\alpha}{2})x^{2}\ .$ (67)

######  LEMMA 2.12.

For all $N$, all $0<\alpha\leq 2$, and for all $\vec{v}\in{\mathcal{S}}_{N}$,

$1-(1-\tfrac{\alpha}{2})\frac{N((N-1)^{2}+1)}{(N-1)^{4}}-\tfrac{\alpha}{2}\frac{1}{(N-1)^{2}}+(1-\tfrac{\alpha}{2})\frac{N+1}{(N-1)^{3}}\
\leq\ W^{(\alpha)}(\vec{v})\ \leq\\\
\left(1-\frac{1}{(N-1)^{2}}\right)^{\alpha/2}\ .$ (68)

Furthermore, for all $\vec{v}\in{\mathcal{S}}_{N}$, with $W^{(\alpha)}$ given
by (63),

$W^{(0)}(\vec{v})\ =\ 1\qquad{\rm and}\qquad\ W^{(2)}(\vec{v})\ =\
1-\frac{1}{(N-1)^{2}}\ .$ (69)

###### Proof.

Repeated use will be made of

$\frac{1}{N}\sum_{k=1}^{N}|v_{k}|^{2}=1\;$ (70)

that identity is part of the definition of ${\mathcal{S}}_{N}$.

Because of (70),
${\displaystyle\frac{1}{N}\sum_{k=1}^{N}\left(\frac{N^{2}-(1+|v_{k}|^{2})N}{(N-1)^{2}}\right)=1-\frac{1}{(N-1)^{2}}}$.
Since $x\mapsto x^{\tfrac{\alpha}{2}}$ is concave on
${\mathord{\mathbb{R}}}_{+}$ for $0\leq\alpha\leq 2$, Jensen’s inequality
yields the upper bound.

To prove the lower bound, use the inequality (67): Writing

$\frac{N^{2}-(1+|v_{k}|^{2})N}{(N-1)^{2}}=1+\frac{N(1-|v_{k}|^{2})-1}{(N-1)^{2}}\
,$ (71)

and applying (67) and (70) yields

$W^{({\alpha})}(\vec{v})\geq
1-\tfrac{\alpha}{2}\frac{1}{(N-1)^{2}}-(1-\tfrac{\alpha}{2})\frac{1}{N}\sum_{k=1}^{N}\left(\frac{N(1-|v_{k}|^{2})-1}{(N-1)^{2}}\right)^{2}\
.$

Expanding the square on the right and applying (70) twice more, we find

$W^{({\alpha})}(\vec{v})\geq
1-\tfrac{\alpha}{2}\frac{1}{(N-1)^{2}}-\frac{1-\tfrac{\alpha}{2}}{(N-1)^{4}}\left[1-N^{2}+N\sum_{k=1}^{N}|v_{k}|^{4}\right]\
.$ (72)

The maximum of ${\displaystyle\sum_{k=1}^{N}|v_{k}|^{4}}$ on
${\mathcal{S}}_{N}$ is no greater than the maximum of the convex function
$\sum_{k=1}^{N}x_{k}^{2}$ on the convex set of $(x_{1},\dots,x_{N})$
satisfying

$0\leq x_{j}\leq N-1\quad{\rm for\ all}\quad j=1,\dots,N\qquad{\rm
and}\qquad\sum_{j=1}^{N}x_{j}=N\ .$ (73)

The extreme points are obtained by permuting the coordinates of
$(N-1,1,0,\dots,0)$. Evaluating the sum at such a point yields the stated
bound, The final statement is obvious. ∎

######  Remark 2.13.

Lemma 2.12 shows that for large $N$,

$W^{(\alpha)}(\vec{v})\geq
1-\left(1-\frac{\alpha}{2}\right)\frac{1}{N}+\mathcal{O}\left(\frac{1}{N^{2}}\right)\
.$ (74)

The fact that the coefficient of $1/N$ is no less than $-1$ is essential for
the result that we shall prove.

We are particularly concerned with the case $\alpha=1$, and shall provide all
the details in this case only. For $\alpha=1$, the lower bound simplifies
further to

$W^{(1)}(\vec{v})\geq
C_{N}:=1-\tfrac{1}{2}\frac{1}{N-1}-\tfrac{1}{2}\frac{1}{(N-1)^{2}}+\tfrac{1}{2}\frac{1}{(N-1)^{3}}-\tfrac{1}{2}\frac{1}{(N-1)^{4}}$
(75)

It is easily seen that for all $N\geq 2$, $C_{N}$ increases as $N$ increases.
For small $N$, we have the explicit values

$C_{3}=\frac{21}{32}\quad{\rm and}\quad C_{4}=\frac{64}{81}\ .$

We now turn to $P^{(\alpha)}$. By (71), for each $k$ for all $v_{k}$,

$w_{N}^{\alpha/2}(v_{k})=\left[\frac{N^{2}-(1+|v_{k}|^{2})N}{(N-1)^{2}}\right]^{\alpha/2}=(1+x_{N}(v_{k}))^{\alpha/2}\
,$ (76)

where

$x_{N}(v_{k})=\frac{1}{N-1}-\frac{N}{(N-1)^{2}}|v_{k}|^{2}\ .$ (77)

Note that ${\displaystyle-1\leq x_{N}(v_{k})\leq\frac{1}{N-1}}$. Then

$w_{N}^{\alpha/2}(v_{k})\leq\left(\frac{N}{N-1}\right)^{\alpha/2}$ (78)

and by (76) and the bounds from Lemma 2.11,
$1+\tfrac{\alpha}{2}x+(\frac{\alpha}{2}-1)x^{2}\leq(1+x)^{\alpha/2}\leq
1+\tfrac{\alpha}{2}x$,

$|w_{N}^{\alpha/2}(v_{k})-1|\leq\frac{1}{N-1}+\frac{3N}{(N-1)^{2}}|v_{k}|^{2}+\frac{N^{2}}{(N-1)^{4}}|v_{k}|^{4}$
(79)

for all $\alpha\in[0,2]$, where we have made estimates to simplify the right
hand side. Thus, while $W^{(\alpha)}$ is only constant for $\alpha=0,2$, it is
nearly constant for all $\alpha\in(0,2)$ when $N$ is large. However, its
range, and hence the spectrum of the multiplication operator specified by
$W^{(\alpha)}$, is a closed interval of positive length. At this point we
record a simple lemma that will be useful later.

######  LEMMA 2.14.

For all $p\geq 1$, there a constant $C$ depending only on $p$, so that for an
$N\geq 3$ and all $\alpha\in[0,2]$,

$\left(\int_{{\mathcal{S}}_{N}}|w_{N}^{\alpha/2}(v_{k})-1|^{p}{\rm
d}\sigma_{N}\right)^{1/p}\leq\frac{C}{N}\ .$ (80)

###### Proof.

This is an immediate consequnce of (79), the triangle inequality, and the fact
that for all $m\in{\mathord{\mathbb{N}}}$,

$\lim_{N\to\infty}\int_{{\mathcal{S}}_{N}}|v_{k}|^{mp}{\rm
d}\sigma_{N}=\int_{{\mathord{\mathbb{R}}}^{3}}|v|^{mp}{\rm d}\gamma\ .$

∎

######  LEMMA 2.15.

For all $\alpha\in[0,2]$, the null space of $P^{(\alpha)}$ is independent of
$\alpha$. If $h$ belongs to the null space of $P^{(0)}$, then $P_{k}h=0$ for
each $k=1,\dots N$. For all $\alpha\in[0,2]$, the closure of the range of
$P^{(\alpha)}$ is the subspace ${\mathcal{A}}_{N}$ of
$L^{2}({\mathcal{S}}_{N})$ defined in Definition 2.1.

###### Proof.

Since $P^{(\alpha)}\geq 0$, $h$ belongs to the null space of $P^{(\alpha)}$ if
and only if $\langle h,P^{(\alpha)}h\rangle=0$. But ${\displaystyle 0=\langle
h,P^{({\alpha})}h\rangle=\frac{1}{N}\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}w_{N}^{\alpha/2}(v_{k})|P_{k}h|^{2}{\rm
d}\sigma_{N}}$. Since $w_{N}^{\alpha/2}(v_{k})\geq 0$ almost everywhere, it
must be the case that $|P_{k}h|^{2}$ vanishes identically. Thus is $h$ is in
the null space of $P^{(\alpha)}$, $P_{k}h=0$ for each $k$, and $h$ is in the
null space of $P^{(0)}$. Conversely if $h$ is in the null space of $P^{(0)}$,
then $P_{k}h=0$ for each $k$, and then clearly $P^{(\alpha)}h=0$.

Since each $P^{(\alpha)}\geq 0$, the closure of its range is the orthogonal
complement of its null space. Since the null space does not depend on
$\alpha$, neither does the range. Evidently, ${\mathcal{A}}_{N}$ is the
closure of the range of $P^{(0)}$. ∎

###  2.3 The spectrum of $\widehat{L}_{N,0}$

Already in our paper [5] we have proved results that specify the exact
spectral gap of $\widehat{L}^{N,\alpha}$ for $\alpha=0$. This case is
especially amenable for several reasons. First, since $W^{(0)}=1$,

$\widehat{L}_{N,0}f=f-P^{({0})}f\ ,$

and hence the problem is to determine the spectrum of $P^{({0})}$. Second,
there is an orthonormal basis of $L^{2}({\mathcal{S}}_{N})$ consisting of
eigenfunctions of $P^{(0}$. This is the case because each $P_{k}$ is an
average of rotations, so the finite dimensional spaces spanned by spherical
harmonics of given maximal degree are invariant under $P^{(0)}$, and therefore
one can study the spectrum of $P^{(0)}$ by studying the eigenvalue equation
$P^{(0)}f=\lambda f$. This is the approach we took in our previous work.
However, this approach cannot work even for $\alpha=2$, the next simplest
case: In this case, $P^{(2)}$ has an interval of continuous spectrum, as we
shall see below. Therefore, we now give another argument that determines the
spectral gap of $\widehat{L}_{N,0}$ that does extend to $\alpha=2$ at least.

######  LEMMA 2.16.

For all $N\geq 3$,

$\widehat{\Delta}_{N,0}=1-\frac{3N-1}{3(N-1)^{2}}=1-\frac{1}{N}+\mathcal{O}\left(\frac{1}{N^{2}}\right)\
.$ (81)

and the second largest eigenvalue of $P^{(0)}$, $\mu^{(0)}$, is given by

$\mu^{(0)}=\frac{3N-1}{3(N-1)^{2}}\ .$ (82)

###### Proof.

The range of $P^{(0)}$ is ${\mathcal{A}}_{N}$, and it suffices to determine
the spectrum of $P^{(0)}$ as an operator on ${\mathcal{A}}_{N}$. For
${\displaystyle
f(\vec{v})=\sum_{j=1}^{N}\varphi_{j}(v_{j})\in{\mathcal{A}}_{N}}$, we compute

$P^{(0)}\left(\sum_{j=1}^{N}\varphi_{j}(v_{j})\right)=\frac{1}{N}\sum_{k=1}^{N}\left(\varphi_{k}(v_{k})+\sum_{j\neq
k,j=1}^{N}K\varphi_{j}(v_{k})\right)\ .$

By this computation, with
$T:\bigoplus_{j=1}^{N}{\mathcal{H}}_{N,j}\to{\mathcal{A}}_{N}$ defined as in
Theorem 2.5,

$T^{-1}P^{(0)}T={\bf M}^{(0)}$

where ${\bf M}^{(0)}=[M_{i,j}^{(0)}]$ is the $N\times N$ block matrix operator
on $\bigoplus_{j=1}^{N}{\mathcal{H}}_{N,j}$ given by

$M_{i,j}^{(0)}=\frac{1}{N}I\quad{\rm if}\ i=j\quad{\rm and}\quad
M_{i,j}^{(0)}=\frac{1}{N}K\quad{\rm if}\ i\neq j\ .$

Note that ${\bf M}^{(0)}$ is unitarily equivalent to the block matrix operator
in $\bigoplus_{j=1}^{N}{\mathcal{H}}_{N,j}$ given by

$\frac{1}{N}\left[\begin{array}[]{ccccc}I+(N-1)K&0&0&\cdots&0\\\
0&I-K&0&\cdots&0\\\ 0&0&I-K&\cdots&0\\\ \vdots&\vdots&\vdots&\ddots&\vdots\\\
0&0&0&0&I-K\end{array}\right]\ .$

It follows that if $\lambda\not=0$ is an eigenvalue of $P^{(0)}$ then either
$\lambda$ is an eigenvalue of $\frac{1}{N}(I+(N-1)K)$ or else $\lambda$ is an
eigenvalue of $\frac{1}{N}(I-K)$. Thus, the second largest eigenvalue of ${\bf
M}^{(0)}$, and hence $P^{(0)}$, is either $1+(N-1)\kappa$, where $\kappa$ is
the second largest eigenvalue of $K$, or else $1-\kappa$ where $\kappa$ is the
least eigenvalue of $K$. From the information on the spectrum of $K$ provided
in Lemma 2.3, one immediately deduces (82), and then (81) follows directly. ∎

###  2.4 The spectrum of $\widehat{L}_{N,\alpha}$, $\alpha\in(0,2]$

After $\alpha=0$, the next simplest case is $\alpha=2$ since then at least
$W^{(2)}$ is constant; as we have seen $W^{(2)}=1-(N-1)^{-2}$. It follows that
$1$ is an eigenfunction for $P^{(2)}$ with eigenvalue $1-(N-1)^{-2}$, and it
spans the eigenspace. That is, $1$ spans the null space of
$\widehat{L}_{N,2}$.

######  LEMMA 2.17.

For all $N\geq 3$, $\widehat{\Delta}_{N,2}>0$.

###### Proof.

The range of $P^{(2)}$ is ${\mathcal{A}}_{N}$, and as with $\alpha=0$, ${\bf
M}^{(2)}:=T^{-1}P^{(2)}T$ has a simple block matrix structure:

$\displaystyle P^{(2)}\left(\sum_{j=1}^{N}\varphi_{j}(v_{j})\right)$
$\displaystyle=$
$\displaystyle\frac{1}{N}\sum_{k=1}^{N}w_{N,2}(v_{k})P_{k}\left(\sum_{j=1}^{N}\varphi_{j}(v_{j})\right)$
$\displaystyle=$
$\displaystyle\sum_{k=1}^{N}\frac{1}{N}w_{N,2}(v_{k})\left(\varphi_{k}(v_{k})+\sum_{j\neq
k,j=1}^{N}K\varphi_{j}(v_{k})\right)$

By Theorem 2.5, ${\bf M}^{(2)}=T^{-1}P^{(2)}T={\bf W}^{(2)}({\bf I}+{\bf C})$,
where

${\bf
W}^{(2)}=\frac{1}{N}\left[\begin{array}[]{cccc}w_{N,2}(v_{1})&0&\cdots&0\\\
0&w_{N,2}(v_{2})&\cdots&0\\\ \vdots&\cdots&\ddots&\vdots\\\
0&\cdots&0&w_{N,N}(v_{N})\end{array}\right]\ ,$

${\bf I}$ is the identity on $\bigoplus_{j=1}^{N}{\mathcal{H}}_{N,j}$, and
${\bf C}$ is given by

${\bf C}=\left[\begin{array}[]{ccccc}0&K&K&\cdots&K\\\ K&0&K&\cdots&K\\\
\vdots&\cdots&\cdots&\ddots&\vdots\\\ K&\cdots&K&0&K\\\
K&\cdots&K&K&0\end{array}\right]\ ,$

Since ${\bf M}^{(2)}$ and $P^{(2)}$ are similar, they have the same spectrum,
and in particular, the spectrum of ${\bf M}^{(2)}$ is real. (This is also
evident from the identity ${\bf M}^{(2)}={\bf W}^{(2)}({\bf I}+{\bf C})$, and
the fact that for bounded operators $A$ and $B$ on any Hilbert space, $AB$ and
$BA$ have the same spectrum.)

Since the range of $\frac{1}{N}w_{N,2}$ is $[0,(N-1)^{-1}]$, this interval is
the spectrum of ${\bf W}^{(2)}$. Note that ${\bf C}$, and hence ${\bf
W}^{(2)}{\bf C}$ is compact. By Weyl’s lemma, the essential spectrum of
$TP^{(2)}T^{-1}$, and hence of $P^{(2)}$, is the essential spectrum of ${\bf
W}^{(2)}$, which is the interval $[0,(N-1)^{-1}]$. Hence any spectrum of
$P^{(2)}$ in $(N-1)^{-1},1-(N-1)^{-2})$ consists of isolated eigenvalues, and
the isolated eigenvalues can only accumulate at a point in $[0,(N-1)^{-1}]$.
In particular, $1-(N-1)^{-2}$ cannot be an accumulation point, and hence
$P^{(2)}$ has a spectral gap below its top eigenvalue $1-(N-1)^{-2}$. This
proves that $\widehat{\Delta}_{N,2}>0$ for all $N\geq 3$. ∎

For $\alpha\in(0,2)$, $W^{(\alpha)}$ is not constant – although for large $N$
it is nearly constant. This means that for such $\alpha$, one cannot determine
the spectrum of $\widehat{L}_{N,\alpha}$ simply by determining the spectrum of
$P^{(\alpha)}$, and moreover, ${\mathcal{A}}_{N}$ is not invariant under
$\widehat{L}_{N,\alpha}$. However, there is a simple comparison that one can
make between ${\cal D}_{N,\alpha}$ and ${\mathcal{D}}_{N,2}$ that provides the
bound on $\widehat{\Delta}_{N,\alpha}$ that we seek.

######  LEMMA 2.18.

For all $N\geq 3$, and all $\alpha\in[0,2]$,

$\widehat{\Delta}_{N,\alpha}\geq\left(\frac{N-1}{N}\right)^{1-\alpha/2}\widehat{\Delta}_{N,2}>0\
.$

###### Proof.

By (78), for all $f$ and $k$, and all $\alpha\in(0,2)$,

$\left[\frac{N^{2}-(1+|v_{k}|^{2})N}{(N-1)^{2}}\right]^{\alpha/2}\left[f-P_{k}f\right]^{2}\geq\left(\frac{N-1}{N}\right)^{1-\alpha/2}\left[\frac{N^{2}-(1+|v_{k}|^{2})N}{(N-1)^{2}}\right]\left[f-P_{k}f\right]^{2}$

It follows immediately that ${\displaystyle{\cal
D}_{N,\alpha}(f,f)\geq\left(\frac{N-1}{N}\right)^{1-\alpha/2}{\mathcal{D}}_{N,2}(f,f)}$,
and then that
${\displaystyle\widehat{\Delta}_{N,\alpha}\geq\left(\frac{N-1}{N}\right)^{1-\alpha/2}\widehat{\Delta}_{N,2}>0}.$
∎

At this point, we have proved the first part of Theorem 1.4, and all that
remains is to prove the second part.

##  3 A sharper lower bound on $\widehat{\Delta}_{N,1}$ for large $N$

In this section we obtain lower bounds on ${\mathcal{D}}_{N,1}(f,f)$ for $f$
orthogonal to the constants that become sharper and sharper as $N$ increases.
To keep the computations simple, we do this explicitly for $\alpha=1$, though
the method applies to all $\alpha\in(0,2)$. We shall prove the following,
which is simply a specialization of Theorem 1.4:

######  THEOREM 3.1.

There is a constant $C$ independent of $N$ such that whenever $f$ is
orthogonal to the constants,

${\mathcal{D}}_{N,1}(f,f)\geq\left(1-\frac{1}{N}-\frac{C}{N^{3/2}}\right){\|}f{\|}_{2}^{2}\
.$ (83)

The bound (83) is meaningless for $N$ such that the right side is negative.
However, no matter what $C$ is, there is an $N_{0}\in{\mathord{\mathbb{N}}}$
such that for all $N\geq N_{0}$, the right side is positive. From that point
on, we have what we need for our induction. The rest of this section is
devoted to the proof of Theorem 3.1

###  3.1 The trial function decomposition

We begin by specifying a trial function decomposition that we shall use. Let
${\mathcal{A}}_{N}$ be the subspace of $L^{2}(\sigma_{N})$ defined in
Definition 2.1. For any $f\in L^{2}({\mathcal{S}}_{N})$ orthogonal to the
constants, define $p$ and $h$ to be the orthogonal projections of $f$ onto
${\mathcal{A}}_{N}$ and ${\mathcal{A}}_{N}^{\perp}$ respectively. Then since
$1\in{\mathcal{A}}_{N}$, $h$ is orthogonal to the constant, and then $p=f-h$
is orthogonal to the constants.

By Lemma 2.15, $h$ is the component of $f$ in the null space of
$P^{({\alpha})}$ for each $\alpha\in[0,2]$, and hence

$\langle f,P^{({\alpha})}f\rangle=\langle p,P^{({\alpha})}p\rangle$ (84)

which yields

$\int_{{\mathcal{S}}_{N}}W^{({\alpha})}f^{2}{\rm d}\sigma-\langle
f,P^{({\alpha})}f\rangle_{L^{2}({\mathcal{S}}_{N})}=\int_{{\mathcal{S}}_{N}}W^{({\alpha})}f^{2}{\rm
d}\sigma-\langle p,P^{({\alpha})}p\rangle_{L^{2}({\mathcal{S}}_{N})}\ .$ (85)

Since $p\in{\mathcal{A}}_{N}$, there are $N$ functions
$\phi_{1},\dots,\phi_{N}$ of a single variable such that $\phi_{j}(v_{j})\in
L^{2}({\mathcal{S}}_{N})$ for each $j$, and

$p(\vec{v})=\sum_{j=1}^{N}\phi_{j}(v_{j})\ ,$ (86)

and we shall always choose the particular representation of this form that is
specificed in Theorem 2.5. That is, the eigenfunctions expansion

$\phi_{j}=\sum_{i=1}^{\infty}a_{j,i}\eta_{i}(v_{j})$ (87)

given in (49) is such that (50) is satsified; i.e., $\sum_{j=1}^{N}a_{j,i}=0$
for $j=1,\dots,4$. We make a further decomposition of $\phi_{j}(v_{j})$ as
follows:

######  DEFINITION 3.2.

Let $p$ be a function given by a sum of the form (86) where for each $j$,
$\phi_{j}(v_{j})$ is orthogonal to the constants, and moreover, (87) is
satisfied with $\sum_{j=1}^{N}a_{j,i}=0$ for $j=1,\dots,4$. Define

$\psi_{j}=\sum_{i=1}^{4}a_{j,i}\eta_{i}(v_{j})\quad{\rm
and}\quad\varphi_{j}(v_{j})=\sum_{i=5}^{\infty}a_{j,i}\eta_{i}(v_{j})$ (88)

so that $\phi_{j}=\psi_{j}+\varphi_{j}$. Next, define

$g(\vec{v})=\sum_{j=1}^{N}\varphi_{j}(v_{j})\quad{\rm and}\quad
s(\vec{v})=\sum_{j=1}^{N}\psi_{j}(v_{j})$ (89)

Finally the trial function decomposition of any $f\in L^{2}(\sigma_{N})$ that
is orthogonal to the constants is given by

$f=g+s+h$ (90)

where $h$ is the component of $f$ in the null space of $P^{(\alpha)}$, $p$ is
the component of $f$ in the closure of the range of $P^{(\alpha)}$, and
$p=g+s$ is the decomposition of $p$ defined in (90).

######  Remark 3.3.

It is easy to see that when $p$ is symmetric under coordinate permutations,
one can take the functions $\phi_{j}$ in (86) to be all the same. In
particular, each $\psi_{j}$ has the form $\psi_{j}(v_{j})=a\cdot
v_{j}+b(|v_{j}|^{2}-1)$ for some fixed $a\in{\mathord{\mathbb{R}}}^{3}$ and
$b\in{\mathord{\mathbb{R}}}$. Then

$s(\vec{v})=a\cdot\left(\sum_{j=1}^{N}v_{j}\right)+b\left(\sum_{j=1}^{N}(|v_{j}|^{2}-1)\right)=0$

on account of the constraints on the momentum and energy. Hence when $p$ is
symmetric $s=0$, and in this case the trial function decomposition simplifies
to $f=g+h$, as in [6].

We have seen in Lemma 2.15 that for each $k$, $P_{k}h=0$. The next lemma shows
that each $P_{k}$ also has a simple action on $s$:

######  LEMMA 3.4.

For each $k$, the function $s$ in the trial function decomposition satisfies

$P_{k}s(\vec{v})=\frac{N-2}{N-1}\psi_{k}(v_{k})\ .$ (91)

###### Proof.

Note that ${\displaystyle
P_{k}s(\vec{v})=\psi_{k}(v_{k})-\frac{1}{N-1}\sum_{j\neq k}\psi_{j}(v_{k})}$.
Writing $\psi_{j}=\sum_{i=1}^{4}a_{j,i}\eta_{i}$ and recalling that
$\sum_{j=1}^{N}a_{j,i}=0$ for $i=1,\dots,4$, for any fixed $v$

$\sum_{j=1}^{N}\psi_{j}(v)=\sum_{j=1}^{N}\sum_{i=0}^{4}a_{j,i}\eta_{i}(v)=\sum_{i=0}^{4}\left(\sum_{j=1}^{N}a_{j,i}\right)\eta_{i}(v)=0\
,$

from which (91) follows. ∎

Each of the components $g$, $s$ and $h$ have their own special properties that
we shall repeatedly use.

(1) A very useful feature of $g(\vec{v})=\sum_{j=1}^{N}\varphi_{j}(v_{j})$ is
that, by Lemma 2.3 for each $j$

${\|}K\varphi_{j}{\|}_{2}\leq\frac{5N-3}{3(N-1)^{3}}{\|}\varphi_{j}{\|}_{2}\
.$

This gives us something almost like Lemma 3.4 for $g$:

$P_{k}g(\vec{v})=\varphi_{k}(v_{k})+\sum_{j\neq k}K\varphi_{j}(v_{k})\ ,$

and hence

${\|}P_{k}g-\varphi_{k}(v_{k}){\|}_{2}\leq\left\|\sum_{j\neq
k}K\varphi_{j}(v_{k})\right\|_{2}\leq\frac{C}{N^{2}}\sum_{k}\|\varphi_{k}\|^{2}\
,$

where in the last inequality we have used Theorem 2.5.

(2) The key feature of $s(\vec{v})=\sum_{j=1}^{N}\psi_{j}(v_{j})$ is that
$P_{k}$ has a very simple action on $s$, given in Lemma 3.4.

Another is that each $\psi_{j}(v_{j})$ belongs to $L^{4}(\sigma_{N})$, and for
a constant $C$ independent of $N$, ${\|}\psi_{j}{\|}_{4}\leq
C{\|}\psi_{j}{\|}_{2}$. This is essentially because the integrals
$\int_{{\mathcal{S}}_{N}}|v|^{2m}{\rm d}\sigma_{N}$ are bounded uniformly in
$N$ for each $m$. In particular, if we wish to estimate the
$L^{2}(\sigma_{N})$ norm of $|v_{k}|^{2}\psi_{k}(v_{k})$, we can apply
Schwarz’s inequality to bound this by $C{\|}\psi_{k}{\|}_{4}$, and then,
changing $C$, to $C{\|}\psi_{k}{\|}_{2}$. This will be used in estimating the
quantity in (100) below.

(3) A very useful feature of $h(\vec{v})$ is that, by Lemma 2.15, $P_{k}h=0$
for each $k$, and in particular, $P^{(1)}h=0$.

###  3.2 Lower bound on ${\mathcal{D}}_{N,1}(f,f)$

For $\alpha=1$, the lower bound (72) simplifies to

$W^{({1})}(\vec{v})\geq\widetilde{W}^{({1})}(\vec{v}):=1+\frac{1}{(N-1)^{3}}-\frac{1}{2}\frac{N}{(N-1)^{4}}\sum_{k=1}^{N}|v_{k}|^{4}\
.$ (92)

Define

$\widetilde{\cal
D}_{N,1}(f,f)=\int_{{\mathcal{S}}_{N}}\widetilde{W}^{({1})}f^{2}{\rm
d}\sigma_{N}-\langle f,P^{(1)}f\rangle\ .$ (93)

By (92), ${\cal D}_{N,1}(f,f)\geq\widetilde{\cal D}_{N,1}(f,f)$.

Now let $f$ be orthogonal to the constants, and let $f=g+s+h$ be the trial
function decomposition of $f$ as specified above. This notation will be used
throughout this subsection. Note that

$\displaystyle\widetilde{\cal D}_{N,1}(f,f)$ $\displaystyle=$
$\displaystyle\widetilde{\cal D}_{N,1}(g,g)+\widetilde{\cal
D}_{N,1}(s,s)+\widetilde{\cal D}_{N,1}(h,h)$ $\displaystyle+$ $\displaystyle
2\widetilde{\cal D}_{N,1}(g,h)+2\widetilde{\cal D}_{N,1}(s,h)+2\widetilde{\cal
D}_{N,1}(g,s)$

The next lemma says that $g$, $s$ and $h$ are almost mutually orthogonal with
respect to the inner product given by $\widetilde{\cal D}_{N,1}$, and hence
the last three terms above make a negligible contribution. This decouples the
contributions of $g$, $s$ and $h$, which may then be analyzed separately,
taking advantage of their different helpful properties.

######  LEMMA 3.5.

There is a constant $C$ independent of $N$ such that for any $f\in
L^{2}(\sigma_{N})$ that is orthogonal to the constants, if $f=g+s+h$ is the
trial function decomposition as specified above, then

$2|{\widetilde{\cal D}_{N,1}}(g,h)|+2|{\widetilde{\cal
D}_{N,1}}(s,h)|+2|{\widetilde{\cal
D}_{N,1}}(g,s)|\leq\frac{C}{N^{3/2}}{\|}f{\|}_{2}^{2}\ .$

###### Proof.

Since $P^{(1)}h=0$, and since $g$ and $h$ are orthogonal, recalling that we
may write $g(\vec{v})=\sum_{j=1}^{N}\varphi_{j}(v_{j})$,

$\displaystyle\widetilde{\cal D}_{N,1}(g,h)$ $\displaystyle=$
$\displaystyle\int_{{\mathcal{S}}_{N}}\widetilde{W}^{({1})}gh{\rm
d}\sigma_{N}=-\frac{1}{2}\frac{N}{(N-1)^{4}}\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}|v_{k}|^{4}gh{\rm
d}\sigma_{N}$ (94) $\displaystyle=$
$\displaystyle-\frac{1}{2}\frac{N}{(N-1)^{4}}\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}|v_{k}|^{4}\varphi_{k}(v_{k})h{\rm
d}\sigma_{N}$ $\displaystyle-$
$\displaystyle\frac{1}{2}\frac{N}{(N-1)^{4}}\sum_{j\neq
k}^{N}\int_{{\mathcal{S}}_{N}}|v_{k}|^{4}\varphi_{j}(v_{j})h{\rm
d}\sigma_{N}{\rm d}\sigma_{N}$ (95)

The integral in (94) vanishes since $P_{k}h=0$. Next consider the integral in
(95). It will be convenient to introduce the notation $\xi(x)=x^{8}$ for the
eighth power. Then, with this definition, the Schwarz inequality, and then
application of the $K$ operator,

$\displaystyle\left|\int_{{\mathcal{S}}_{N}}|v_{k}|^{4}\varphi_{j}(v_{j})h{\rm
d}\sigma_{N}\right|$ $\displaystyle\leq$
$\displaystyle{\|}h{\|}_{2}\left(\int_{{\mathcal{S}}_{N}}|v_{k}|^{8}\varphi_{j}^{2}(v_{j}){\rm
d}\sigma_{N}\right)^{1/2}$ (96) $\displaystyle=$
$\displaystyle{\|}h{\|}_{2}\left(\int_{{\mathcal{S}}_{N}}K\xi(v_{j})\varphi_{j}^{2}(v_{j}){\rm
d}\sigma_{N}\right)^{1/2}\ .$

By Lemma 2.7, there is a constant $C$ so that, independent of $N$,
${\|}K\xi{\|}_{\infty}\leq C$. Therefore,
${\displaystyle\left|\int_{{\mathcal{S}}_{N}}h|v_{k}|^{4}\varphi_{j}(v_{j}){\rm
d}\sigma_{N}\right|\leq C{\|}h{\|}_{2}{\|}\varphi_{j}{\|}_{2}}$. Using this in
(95) gives us

$\left|\frac{N}{(N-1)^{4}}\int_{{\mathcal{S}}_{N}}\left(\sum_{k=1}^{N}|v_{k}|^{4}\right)gh{\rm
d}\sigma\right|\leq\frac{N}{(N-1)^{3}}C{\|}h{\|}_{2}\left(\sum_{j=1}^{N}{\|}\varphi_{j}{\|}_{2}\right)\
,$ (97)

and then since Theorem 2.5 gives us

$\sum_{j=1}^{N}{\|}\varphi{\|}_{2}\leq\left(1-\frac{5N-3}{3(N-1)^{2}}\right)^{-1/2}\sqrt{N}{\|}g{\|}_{2}\
,$

we have that the left hand side of (97) is bounded by
$\frac{C}{N^{3/2}}{\|}g{\|}_{2}{\|}h{\|}_{2}$ for a constant $C$ independent
of $N$. We conclude that $|\widetilde{\cal D}_{N,1}(s,h)|\leq CN^{-3/2}$.

Finally, we consider $\widetilde{\cal D}_{N,1}(s,g)$. This time we must also
estimate $\langle s,P^{(1)}g\rangle$. Because the span of $\\{\eta_{j}(v_{j})\
:\ 1\leq j\leq N\\}$ is invariant under $P^{(0)}$, and every function in it is
orthogonal to $g$,

$\langle s,P^{(1)}g\rangle=\langle s,P^{(1)}g\rangle-\langle
P^{(0)}s,g\rangle=\langle s,(P^{(1)}-P^{(0)})g\rangle\ .$

Introducing the short notation $\tilde{w}(v):=w_{N,1}(v)-1$ to be used in this
proof only, and writing

$s(\vec{v})=\sum_{\j=1}^{N}\psi_{j}(v_{j})\quad{\rm and}\quad
g(\vec{v})=\sum_{j=1}^{N}\varphi_{\ell}(v_{\ell})\ ,$

we have

$\langle
s,P^{(1)}g\rangle=\sum_{j,k,\ell=1}^{N}\int_{{\mathcal{S}}_{N}}\psi_{j}(v_{j})\tilde{w}(v_{k})P_{k}\varphi_{\ell}(v_{\ell})\
.$ (98)

We now split the sum over $j$, $k$ and $\ell$, into five parts

${\it(i)}\ j=\ell=k\,\quad{\it(ii)}\ j\neq k,\ell=k\,\quad{\it(iii)}\
j=k,\ell\neq k\quad{\rm and}\quad{\it(iv)}\ j=\ell,\ell\neq k\ ,$

and finally, ${\it(v)}\ j\neq\ell,\ell\neq k,k\neq j$

$\displaystyle\langle s,P^{(1)}g\rangle$ $\displaystyle=$
$\displaystyle\langle s,(P^{(1)}-P^{(0)})g\rangle$ (99) $\displaystyle=$
$\displaystyle\frac{1}{N}\sum_{k=1}^{N}\langle\tilde{w}(v_{k})\psi_{k}(v_{k}),\varphi_{k}(v_{k})\rangle$
(100) $\displaystyle+$ $\displaystyle\frac{1}{N}\sum_{k=1}^{N}\sum_{j\neq
k}^{N}\langle\tilde{w}(v_{k})\psi_{j}(v_{j}),\varphi_{k}(v_{k})\rangle$ (101)
$\displaystyle+$ $\displaystyle\frac{1}{N}\sum_{k=1}^{N}\sum_{\ell\neq
k}^{N}\langle\tilde{w}(v_{k})\psi_{k}(v_{k}),P_{k}\varphi_{\ell}(v_{k})\rangle$
(102) $\displaystyle+$ $\displaystyle\frac{1}{N}\sum_{k=1}^{N}\sum_{\ell\neq
k}^{N}\langle\tilde{w}(v_{k})\psi_{\ell}(v_{k}),P_{k}\varphi_{\ell}(v_{k})\rangle$
(103) $\displaystyle+$ $\displaystyle\frac{1}{N}\sum_{j\neq
k,k\neq\ell,\ell\neq
j}^{N}\langle\tilde{w}(v_{k})\psi_{j}(v_{j}),P_{k}\varphi_{\ell}(v_{k})\rangle$
(104)

We estimate (100) as follows, using Lemma 2.14 to bound
${\|}\tilde{w}(v_{k}){\|}_{4}$:

$\displaystyle\frac{1}{N}\left|\sum_{k=1}^{N}\langle\tilde{w}(v_{k})\psi_{k}(v_{k}),\varphi_{k}(v_{k})\rangle\right|$
$\displaystyle\leq$
$\displaystyle\frac{1}{N}\sum_{k=1}^{N}{\|}\tilde{w}(v_{k}){\|}_{4}{\|}\psi_{k}{\|}_{4}{\|}\varphi_{k}{\|}_{2}\leq\frac{C}{N^{2}}\sum_{k=1}^{N}{\|}\psi_{k}{\|}_{2}{\|}\varphi_{k}{\|}_{2}$
$\displaystyle\leq$
$\displaystyle\leq\frac{C}{N^{2}}({\|}s{\|}_{2}^{2}+{\|}g{\|}_{2}^{2})\ .$

Since $P_{k}\psi_{j}(v_{k})=-\frac{1}{N-1}\psi_{j}(v_{k})$, the argument used
to estimate (100) shows that the absolute value of the sum in (101) is bounded
above by

$\frac{C}{N^{3}}\sum_{k=1}^{N}\sum_{j\neq
k}^{N}{\|}\psi_{j}{\|}_{2}{\|}\varphi_{k}{\|}_{2}\leq\frac{C}{N^{2}}({\|}s{\|}_{2}^{2}+{\|}g{\|}_{2}^{2})\
,$

as we found for (100). Since for $k\neq j$, ${\|}P_{k}\varphi_{j}{\|}_{2}\leq
CN^{-2}{\|}\varphi_{j}{\|}_{2}$ (by (44)), the argument used to estimate (100)
shows that the absolute value of the sum in (102) is bounded above by

$\frac{C}{N^{4}}\sum_{k=1}^{N}\sum_{j\neq
k}^{N}{\|}\psi_{k}{\|}_{2}{\|}\varphi_{j}{\|}_{2}\leq\frac{C}{N^{3}}({\|}s{\|}_{2}^{2}+{\|}g{\|}_{2}^{2})\
,$

even better than the previous bounds. Finally, for the terms in (104),

$\displaystyle|\langle\tilde{w}(v_{k})\psi_{j}(v_{j}),P_{k}\varphi_{\ell}(v_{\ell})\rangle|$
$\displaystyle=$ $\displaystyle|\langle
P_{k}\varphi_{\ell}(v_{\ell})\tilde{w}(v_{k}),P_{k}\psi_{j}(v_{j})\rangle|$
$\displaystyle=$ $\displaystyle\frac{1}{N-1}|\langle
P_{k}\varphi_{\ell}(v_{\ell})\tilde{w}(v_{k}),\psi_{j}(v_{j})\rangle|$
$\displaystyle\leq$
$\displaystyle\frac{1}{N-1}{\|}K\varphi_{\ell}{\|}_{2}{\|}\psi_{j}{\|}_{2}\leq\frac{C}{N^{3}}{\|}\varphi_{\ell}{\|}_{2}{\|}\psi_{j}{\|}_{2}\
,$

where in the last inequality we have used Lemma 2.3, and the fact that for
each $j$, $\psi_{j}$ is an eigenfunction of $K$ considered as an operator on
$L^{2}({\mathcal{S}}_{N})$, with eigenvalue $-\frac{1}{N-1}$. Thus,

$\frac{C}{N^{4}}\sum_{j\neq k,k\neq\ell,\ell\neq
j}^{N}{\|}\psi_{j}{\|}_{2}{\|}\varphi_{\ell}{\|}_{2}\leq\frac{C}{N^{2}}({\|}s{\|}_{2}^{2}+{\|}g{\|}_{2}^{2})\
.$

This proves $|\langle
s,P^{(1)}g\rangle|\leq\frac{C}{N^{2}}({\|}s{\|}_{2}^{2}+{\|}g{\|}_{2}^{2})$. ∎

We now turn to the estimation of $\widetilde{\cal D}_{N,1}(g,g)$ and
$\widetilde{\cal D}_{N,1}(s,s)$.

######  LEMMA 3.6.

There is a constant $C$ independent of $N\geq 3$ such that for all $g$ and $s$
as above,

$\langle
g,P^{({1})}g\rangle\leq\frac{1}{N}\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}\left[\frac{N^{2}-(1+|v_{k}|^{2})N}{(N-1)^{2}}\right]^{1/2}\varphi^{2}_{k}(v_{k}){\rm
d}\sigma+\frac{C}{N^{2}}{\|}g{\|}_{2}^{2}\ .$ (105)

and

$\langle
s,P^{({1})}s\rangle\leq\frac{1}{N}\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}\left[\frac{N^{2}-(1+|v_{k}|^{2})N}{(N-1)^{2}}\right]^{1/2}\psi_{k}^{2}(v_{k}){\rm
d}\sigma\ .$ (106)

###### Proof.

Note first of all that $P_{k}g=\varphi_{k}(v_{k})+\sum_{j\neq
k}K\varphi_{j}(v_{k})$, and thus

$\displaystyle\langle g,P^{({1})}g\rangle$ $\displaystyle=$
$\displaystyle\frac{1}{N}\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}\left[\frac{N^{2}-(1+|v_{k}|^{2})N}{(N-1)^{2}}\right]^{1/2}|P_{k}g|^{2}{\rm
d}\sigma_{N}$ (107) $\displaystyle=$
$\displaystyle\frac{1}{N}\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}\left[\frac{N^{2}-(1+|v_{k}|^{2})N}{(N-1)^{2}}\right]^{1/2}\varphi^{2}_{k}(v_{k}){\rm
d}\sigma_{N}$ $\displaystyle+$
$\displaystyle\frac{2}{N}\sum_{k=1}^{N}\sum_{j\neq
k}\int_{{\mathcal{S}}_{N}}\left[\frac{N^{2}-(1+|v_{k}|^{2})N}{(N-1)^{2}}\right]^{1/2}\varphi_{k}(v_{k})K\varphi_{j}(v_{k}){\rm
d}\sigma_{N}$ $\displaystyle+$
$\displaystyle\frac{1}{N}\sum_{k=1}^{N}\sum_{j\neq k}\sum_{\ell\neq
k}\int_{{\mathcal{S}}_{N}}\left[\frac{N^{2}-(1+|v_{k}|^{2})N}{(N-1)^{2}}\right]^{1/2}K\varphi_{j}(v_{k})K\varphi_{\ell}(v_{k}){\rm
d}\sigma_{N}$ (108)

By the Schwarz inequality and (44), the sum of integrals in (107) is bounded
above by

$\frac{C}{N^{3}}\sum_{k=1}^{N}\sum_{j\neq
k}{\|}\varphi_{j}{\|}_{2}{\|}\varphi_{k}{\|}_{2}\leq\frac{C}{N^{2}}{\|}g{\|}_{2}^{2}\
.$

Similarly, by (44) and Lemma 2.5, the sum of integrals in (108) is bounded
above by

$\frac{C}{N^{4}}\sum_{j,k=1}^{N}{\|}\varphi_{j}{\|}_{2}{\|}\varphi_{k}{\|}_{2}\leq\frac{C}{N^{3}}{\|}g{\|}_{2}^{2}\
.$

Using the two bounds we have just derived on (107) and (108) respectively,
yields (105).

$\langle
s,P^{({1})}s\rangle=\frac{1}{N}\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}\left[\frac{N^{2}-(1+|v_{k}|^{2})N}{(N-1)^{2}}\right]^{1/2}|P_{k}s|^{2}{\rm
d}\sigma_{N}\ ,$

(106) follows directly from Lemma 3.4. ∎

######  LEMMA 3.7.

There is a constant $C$ such that for all $N$ and all $g$ and $s$ as above,

$\int_{{\mathcal{S}}_{N}}\sum_{k=1}^{N}|v_{k}|^{4}g^{2}{\rm
d}\sigma_{N}\leq\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}\varphi_{k}(v_{k})^{2}|v_{k}|^{4}{\rm
d}\sigma+{C}{N}{\|}g{\|}_{2}^{2}\ ,$ (109)

and

$\int_{{\mathcal{S}}_{N}}\sum_{k=1}^{N}|v_{k}|^{4}s^{2}{\rm
d}\sigma_{N}\leq\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}\psi_{k}(v_{k})^{2}|v_{k}|^{4}{\rm
d}\sigma+{C}{N}{\|}s{\|}_{2}^{2}\ ,$ (110)

###### Proof.

$\displaystyle\int_{{\mathcal{S}}_{N}}\sum_{k=1}^{N}|v_{k}|^{4}g^{2}{\rm
d}\sigma$ $\displaystyle=$
$\displaystyle\sum_{i,j,k=1}^{N}\int_{ST}\varphi_{i}(v_{i})\varphi_{j}(v_{j})|v_{k}|^{4}{\rm
d}\sigma_{N}$ (111) $\displaystyle=$
$\displaystyle\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}\varphi_{k}(v_{k})^{2}|v_{k}|^{4}{\rm
d}\sigma$ $\displaystyle+$ $\displaystyle 2\sum_{k=1}^{N}\sum_{j\neq
k}\int_{{\mathcal{S}}_{N}}\varphi_{j}(v_{j})\varphi_{k}(v_{k})|v_{k}|^{4}{\rm
d}\sigma$ (112) $\displaystyle+$ $\displaystyle\sum_{k=1}^{N}\sum_{j\neq
k}\int_{{\mathcal{S}}_{N}}\varphi_{j}(v_{j})^{2}|v_{k}|^{4}{\rm d}\sigma$
(113) $\displaystyle+$ $\displaystyle\sum_{i\neq j,j\neq k,k\neq
i}\int_{{\mathcal{S}}_{N}}\varphi_{i}(v_{i})\varphi_{j}(v_{j})|v_{k}|^{4}{\rm
d}\sigma$ (114)

By Lemma 3.8, Lemma 3.9 and Lemma 3.10 below the terms in (112), (113) and
(114) add up to no more than $CN{\|}g{\|}_{2}^{2}$, which proves (109). The
same argument using the same lemmas proves (110). ∎

######  LEMMA 3.8.

There is a constant $C$ such that for all $N$ and all $g$ and $s$ as above,

$2\sum_{k=1}^{N}\sum_{j\neq
k}\int_{{\mathcal{S}}_{N}}\varphi_{j}(v_{j})\varphi_{k}(v_{k})|v_{k}|^{4}{\rm
d}\sigma\leq{C}{N}{\|}g{\|}_{2}^{2}\ .$ (115)

and

$2\sum_{k=1}^{N}\sum_{j\neq
k}\int_{{\mathcal{S}}_{N}}\psi_{j}(v_{j})\psi_{k}(v_{k})|v_{k}|^{4}{\rm
d}\sigma\leq{C}{\|}s{\|}_{2}^{2}\ .$ (116)

###### Proof.

For $j\neq k$, using the pointwise bound $|v_{k}|^{4}\leq(N-1)^{2}$ and then
(44),

$\int_{{\mathcal{S}}_{N}}\varphi_{j}(v_{j})\varphi_{k}(v_{k})|v_{k}|^{4}{\rm
d}\sigma\leq(N-1)^{2}{\|}K\varphi_{j}{\|}_{2}{\|}\varphi_{k}{\|}_{2}\leq\frac{(5N-3)(N-1)^{2}}{3(N-1)^{3}}{\|}\varphi_{j}{\|}_{2}{\|}\varphi_{k}{\|}_{2}\
.$ (117)

Then by Theorem 2.5, (115) follows. Next,

$\int_{{\mathcal{S}}_{N}}\psi_{j}(v_{j})\psi_{k}(v_{k})|v_{k}|^{4}{\rm
d}\sigma\leq{\|}K\psi_{j}{\|}_{2}{\|}|v_{k}|^{4}\psi_{k}{\|}_{2}\leq\frac{1}{N-1}{\|}\psi_{j}{\|}_{2}C{\|}\psi_{k}{\|}_{4}\leq\frac{C}{N}{\|}\psi_{j}{\|}_{2}C{\|}\psi_{k}{\|}_{2}\
.$ (118)

Then by Theorem 2.5 again, (116) follows. ∎

######  LEMMA 3.9.

There is a constant $C$ such that for all $N$ and all $g$ and $s$ as above,

$\sum_{k=1}^{N}\sum_{j\neq
k}\int_{{\mathcal{S}}_{N}}\varphi_{j}(v_{j})^{2}|v_{k}|^{4}{\rm
d}\sigma_{N}\leq CN{\|}g{\|}_{2}^{2}\ .$ (119)

and

$\sum_{k=1}^{N}\sum_{j\neq
k}\int_{{\mathcal{S}}_{N}}\psi_{j}(v_{j})^{2}|v_{k}|^{4}{\rm d}\sigma_{N}\leq
CN{\|}s{\|}_{2}^{2}\ .$ (120)

###### Proof.

By Lemma 2.8, there is a finite constant $C$ independent of $N$ such that

$\sum_{k=1}^{N}\sum_{j\neq
k}\int_{{\mathcal{S}}_{N}}\varphi_{j}(v_{j})^{2}|v_{k}|^{4}{\rm
d}\sigma_{N}\leq
N\sum_{j=1}^{N}\int_{{\mathcal{S}}_{N}}\frac{N^{2}+|v_{j}|^{4}-2N|v_{j}|^{2}}{(N-1)^{2}}\varphi^{2}_{j}(v_{j}){\rm
d}\sigma_{N}+CN\sum_{j=1}^{N}{\|}\varphi_{j}{\|}_{2}^{2}\ .$

Then by Lemma 2.5, (119) follows. The same analysis yields (120) . ∎

######  LEMMA 3.10.

There is a constant $C$ such that for all $N$ and all $g$ and $s$ as above,

$\sum_{i\neq j,j\neq k,k\neq
i}\int_{{\mathcal{S}}_{N}}\varphi_{i}(v_{i})\varphi_{j}(v_{j})|v_{k}|^{4}{\rm
d}\sigma\leq C{\|}g{\|}_{2}^{2}\ .$ (121)

and

$\sum_{i\neq j,j\neq k,k\neq
i}\int_{{\mathcal{S}}_{N}}\psi_{i}(v_{i})\psi_{j}(v_{j})|v_{k}|^{4}{\rm
d}\sigma\leq C{\|}s{\|}_{2}^{2}\ .$ (122)

###### Proof.

By Lemma 2.9, there is a finite constant $C$ independent of $N$ (but changing
from line to line) such that

$\displaystyle\sum_{i\neq j,j\neq k,k\neq
i}\int_{{\mathcal{S}}_{N}}\varphi_{i}(v_{i})\varphi_{j}(v_{j})|v_{k}|^{4}{\rm
d}\sigma\leq\frac{C}{N}\sum_{i\not=j}{\|}\varphi_{i}{\|}_{2}{\|}\varphi_{j}{\|}_{2}+$
$\displaystyle(N-2)\sum_{i\neq
j}\int_{{\mathcal{S}}_{N}}\frac{N^{2}+|v_{i}|^{4}+|v_{j}|^{4}+2N|v_{i}|^{2}+2N|v_{j}|^{2}+2|v_{i}|^{2}|v_{j}|^{2}}{(N-2)^{2}}\varphi_{i}(v_{i})\varphi_{j}(v_{j}){\rm
d}\sigma_{N}\ .$

Note that for $i\neq j$,

$\int_{{\mathcal{S}}_{N}}\frac{N^{2}+|v_{i}|^{4}+|v_{j}|^{4}+2N|v_{i}|^{2}+2N|v_{j}|^{2}}{N-2}\varphi_{i}(v_{i})\varphi_{j}(v_{j}){\rm
d}\sigma_{N}\leq\frac{C}{N}{\|}\varphi_{i}{\|}_{2}{\|}\varphi_{j}{\|}_{2}\ $

since in each term we may either replace $\varphi_{i}$ by $K\varphi_{i}$ or
$\varphi_{j}$ by $K\varphi_{j}$, and this gives a factor of $CN^{-2}$. Then by
Lemma 2.5, $\sum_{i\neq j}{\|}\varphi_{i}{\|}_{2}{\|}\varphi_{j}{\|}_{2}\leq
CN{\|}g{\|}_{2}^{2}$.

The remaining terms must be handled differently. For $j=1,\dots,N$, let
$\xi_{j}$ denote the function $\xi_{j}(v)=|v_{j}|^{2}\varphi_{j}(v_{j})$, and
note that $\xi_{j}$ is orthogonal to the constants. Therefore,

$\displaystyle\int_{{\mathcal{S}}_{N}}\frac{|v_{i}|^{2}|v_{j}|^{2}}{N-2}\varphi_{i}(v_{i})\varphi_{j}(v_{j}){\rm
d}\sigma_{N}$ $\displaystyle=$
$\displaystyle\frac{1}{N-2}\langle\xi_{i},K\xi_{j}\rangle$ $\displaystyle\leq$
$\displaystyle\frac{1}{N-2}\frac{1}{N-1}{\|}\xi_{i}{\|}_{2}{\|}\xi_{j}{\|}_{2}\leq
C{\|}\varphi_{i}{\|}_{2}{\|}\varphi_{j}{\|}_{2}$

Then by Lemma 2.5, $\sum_{i\neq
j}{\|}\varphi_{i}{\|}_{2}{\|}\varphi_{j}{\|}_{2}\leq CN{\|}g{\|}_{2}^{2}$, and
(121) follows.

Next,

$\displaystyle\sum_{i\neq j,j\neq k,k\neq
i}\int_{{\mathcal{S}}_{N}}\psi_{i}(v_{i})\psi_{j}(v_{j})|v_{k}|^{4}{\rm
d}\sigma\leq CN^{2}{\|}\psi_{i}{\|}_{2}{\|}\psi_{j}{\|}_{2}+$
$\displaystyle(N-2)\sum_{i\neq
j}\int_{{\mathcal{S}}_{N}}\frac{N^{2}+|v_{i}|^{4}+|v_{j}|^{4}+2N|v_{i}|^{2}+2N|v_{j}|^{2}+2|v_{i}|^{2}|v_{j}|^{2}}{(N-2)^{2}}\psi_{i}(v_{i})\psi_{j}(v_{j}){\rm
d}\sigma_{N}\ .$

The main term is

$\frac{N^{2}}{N-2}\sum_{i\neq
j}\int_{{\mathcal{S}}_{N}}\psi_{i}(v_{i})\psi_{j}(v_{j}){\rm
d}\sigma_{N}=\frac{N^{2}}{N-2}\sum_{i\neq
j}\langle\psi_{1},K\psi_{j}\rangle=\frac{N^{2}}{(N-2)(N-1)}\sum_{i\neq
j}{\|}\psi_{i}{\|}_{2}{\|}\psi_{j}{\|}_{2}\ ,$

and simple estimates show that all remaining terms are smaller. ∎

###  3.3 Lower bounds on $\widetilde{\cal D}_{N,1}(g,g)$ and $\widetilde{\cal
D}_{N,1}(s,s)$

We are now ready to estimate $\widetilde{\cal D}_{N,1}(g,g)$ and
$\widetilde{\cal D}_{N,1}(s,s)$. We first define a quadratic form
${\mathcal{F}}$ on $L^{2}(\sigma_{N})$ as follows: For all functions $r$ in
$L^{2}(\sigma_{N})$, define

${\mathcal{F}}(r,r):=\frac{1}{2}\frac{N}{(N-1)^{4}}\int_{{\mathcal{S}}_{N}}\sum_{k=1}^{N}|v_{k}|^{4}r^{2}{\rm
d}\sigma+\langle r,P^{({1})}r\rangle$ (123)

######  LEMMA 3.11.

For all $g$ and $s$ as above,

$\widetilde{\cal
D}_{N,1}(g,g)\geq{\|}g{\|}_{2}^{2}-{\mathcal{F}}(g,g)\quad{\rm
and}\quad\widetilde{\cal
D}_{N,1}(s,s)\geq{\|}s{\|}_{2}^{2}-{\mathcal{F}}(s,s)\ .$ (124)

###### Proof.

This is immediate from (72), (93) and the definition of ${\mathcal{F}}$. ∎

######  LEMMA 3.12.

There is a finite constant $C$ independent of $N$ such that for all $g$ and
$s$ as above, with ${\mathcal{F}}$ defined by (123)

${\mathcal{F}}(g,g)\leq\left(\frac{1}{N}+\frac{C}{N^{2}}\right){\|}g{\|}_{2}^{2}\qquad{\rm
and}\qquad{\mathcal{F}}(s,s)\leq\left(\frac{1}{N}+\frac{C}{N^{2}}\right){\|}s{\|}_{2}^{2}\
.$ (125)

###### Proof.

By Lemma 3.6 and Lemma 3.7,

$\displaystyle{\mathcal{F}}(g,g)$ $\displaystyle\leq$
$\displaystyle\frac{1}{N}\sum_{k=1}^{N}\left[\int_{{\mathcal{S}}_{N}}\left[\frac{N^{2}-(1+|v_{k}|^{2})N}{(N-1)^{2}}\right]^{1/2}+\frac{1}{2}\frac{N^{2}}{(N-1)^{4}}\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}|v_{k}|^{4}\right]\varphi(v_{k})^{2}{\rm
d}\sigma$ $\displaystyle+$ $\displaystyle\frac{C}{N^{2}}{\|}g{\|}_{2}^{2}$

Define ${\displaystyle y_{k}:=\frac{N}{(N-1)^{2}}|v_{k}|^{2}}.$ Then $0\leq
y_{k}\leq N/(N-1)$, and

${\mathcal{F}}(g,g)\leq\frac{1}{N}\sum_{k=1}^{N}\int_{{\mathcal{S}}_{N}}w(y_{k})\varphi^{2}(v_{k}){\rm
d}\sigma_{N}+\frac{C}{N^{2}}{\|}g{\|}_{2}^{2}\ $ (126)

where ${\displaystyle
w(y)=\left(\frac{N}{N-1}-y\right)^{1/2}+\frac{1}{2}y^{2}}$. Simple
calculations show that $w(y)\leq\sqrt{N/(N-1)}$ for all $0\leq y\leq N/(N-1)$,
and in fact, for $N\geq 7$, $w(y)$ is monotone decreasing on this interval.
Then (126) becomes

${\mathcal{F}}(g,g)\leq\sqrt{N/(N-1)}\frac{1}{N}\sum_{k=1}^{N}{\|}\varphi_{k}{\|}_{2}^{2}+\frac{C}{N^{2}}{\|}g{\|}_{2}^{2}$

Now (125) follows directly from Theorem 2.5. The proof of (125) is the same. ∎

###  3.4 Proof of Theorem 3.1

###### Proof of Theorem 3.1.

By Lemma 3.5,

${\cal D}_{N,1}(f,f)\geq\widetilde{\cal D}_{N,1}(f,f)\geq\widetilde{\cal
D}_{N,1}(g,g)+\widetilde{\cal D}_{N,1}(s,s)+\widetilde{\cal
D}_{N,1}(h,h)-CN^{-3/2}\|f\|^{2}s\ .$

By Lemma 3.11 and Lemma 3.12,

$\widetilde{\cal D}_{N,1}(g,g)+\widetilde{\cal
D}_{N,1}(s,s)\geq\left(1-\frac{1}{N}-\frac{C}{N^{2}}\right)({\|}g{\|}_{2}^{2}+{\|}s{\|}_{2}^{2})\
.$

Since $P^{(1)}h=0$, (74) yields ${\displaystyle\widetilde{\cal
D}_{N,1}(h,h)\geq\left(1-\frac{1}{2N}-\frac{C}{N^{2}}\right){\|}h{\|}_{2}^{2}}$,
adding the estimates completes the proof since
${\|}f{\|}_{2}^{2}={\|}g{\|}_{2}^{2}+{\|}s{\|}_{2}^{2}+{\|}h{\|}_{2}^{2}$. ∎

##  Appendix A Some computational proofs

###### Proof of Lemma 2.8.

By (24),

$E\\{|v_{1}|^{4}\ |\ v_{N}=v\\}=\\\
\int_{{\mathcal{S}}_{N-1}}\left(\eta^{4}(v)|\vec{y}|^{4}+\frac{\eta^{2}(v)}{(N-1)^{2}}|\vec{y}\cdot
v|^{2}+\frac{|v|^{4}}{(N-1)^{4}}+2\frac{\eta^{2}(v)}{(N-1)^{2}}|\vec{y}|^{2}\right){\rm
d}\sigma_{N-1}\ ,$ (127)

where

$\eta^{2}(v,w)=\frac{N-|v|^{2}-|v|^{2}/(N-1)}{N-1}\ ,$

Define $M_{N}:=\int_{{\mathcal{S}}_{N-2}}|\vec{y}|^{4}{\rm d}\sigma_{N-2}$
which is bounded uniformly in $N$:

$\lim_{N\to\infty}\int_{{\mathcal{S}}_{N-1}}|\vec{y}|^{4}{\rm
d}\sigma_{N-1}=(2\pi/3)^{-3/2}\int_{{\mathord{\mathbb{R}}}^{3}}|y|^{4}e^{-3|y|^{2}/2}\
.$

Then the right hand side of (127) becomes

$M_{N}\eta^{4}(v)+\frac{1}{3(N-1)^{2}}\eta^{2}(v)|v|^{2}+\frac{|v|^{4}}{(N-1)^{4}}+2\frac{\eta^{2}(v)}{(N-1)^{2}}$
(128)

Note that for some constant $C$ independent of $N$,

$\frac{1}{3(N-1)^{2}}\eta^{2}(v)|v|^{2}+\frac{|v|^{4}}{(N-1)^{4}}+2\frac{\eta^{2}(v)}{(N-1)^{2}}\leq\frac{C}{N}\
.$ (129)

Next,

$\displaystyle\eta^{4}(v)=\frac{N^{2}+|v|^{4}-2N|v|^{2}}{(N-1)^{2}}+\frac{|v|^{4}}{(N-1)^{4}}+2\frac{(N-|v|^{2})|v|^{2}}{(N-1)^{3}}\
.$

Again, for some constant $C$ independent of $N$,

$\frac{|v|^{4}}{(N-1)^{4}}+2\frac{(N-|v|^{2})|v|^{2}}{(N-1)^{3}}\leq\frac{C}{N}\
.$

∎

###### Proof of Lemma 2.9.

By a simple adaptation of (24),

$E\\{|v_{1}|^{4}\ |\ (v_{N-1},v_{N})=(v,w)\\}=\\\
\int_{{\mathcal{S}}_{N-2}}\left(\beta^{4}(v,w)|\vec{y}|^{4}+\frac{\beta^{2}(v,w)}{(N-2)^{2}}|\vec{y}\cdot(v+w)|^{2}+\frac{|v+w|^{4}}{(N-2)^{4}}+2\frac{\beta^{2}(v,w)}{(N-2)^{2}}|\vec{y}|^{2}\right){\rm
d}\sigma_{N-2}\ ,$ (130)

where

$\beta^{2}(v,w)=\frac{N-|v|^{2}-|w|^{2}-|v+w|^{2}/(N-2)}{N-2}\ ,$

which is non-negative on the allowed values for $(v,w)$. Note that
$\beta^{2}(v,w)\leq N/(N-2)$. Define
$M_{N}:=\int_{{\mathcal{S}}_{N-2}}|\vec{y}|^{4}{\rm d}\sigma_{N-2}$ which is
bounded uniformly in $N$:

$\lim_{N\to\infty}\int_{{\mathcal{S}}_{N-2}}|\vec{y}|^{4}{\rm
d}\sigma_{N-2}=(2\pi/3)^{-3/2}\int_{{\mathord{\mathbb{R}}}^{3}}|y|^{4}e^{-3|y|^{2}/2}\
.$

Then the right hand side of (130) becomes

$M_{N}\beta^{4}(v,w)+\frac{1}{3(N-2)^{2}}\beta^{2}(v,w)|v+w|^{2}+\frac{|v+w|^{4}}{(N-2)^{4}}+2\frac{\beta^{2}(v,w)}{(N-2)^{2}}$
(131)

Note that for some constant $C$ independent of $N$,

$\frac{1}{3(N-2)^{2}}\beta^{2}(v,w)|v+w|^{2}+\frac{|v+w|^{4}}{(N-2)^{4}}+2\frac{\beta^{2}(v,w)}{(N-2)^{2}}\leq\frac{C}{N}\
.$ (132)

Next,

$\displaystyle\beta^{4}(v,w)$ $\displaystyle=$
$\displaystyle\frac{N^{2}+|v|^{4}+|w|^{4}+2N|v|^{2}+2N|w|^{2}+2|v|^{2}|w|^{2}}{(N-2)^{2}}$
$\displaystyle+$
$\displaystyle\frac{|v+w|^{4}}{(N-2)^{4}}+2\frac{(N-|v|^{2}-|w|^{2})|v+w|^{2}}{(N-2)^{3}}\
.$

Again, for some constant $C$ independent of $N$,

$\frac{|v+w|^{4}}{(N-2)^{4}}+2\frac{(N-|v|^{2}-|w|^{2})|v+w|^{2}}{(N-2)^{3}}\leq\frac{C}{N}\
.$

∎

##  Appendix B Quantitiative estimates on $\widehat{\Delta}_{N,2}$

###  B.1 An explicit bound for $N\geq 4$

By (78), $P^{(\alpha)}$ defined in (61) satisfies

$0\leq P^{(\alpha)}\leq\left(\frac{N}{N-1}\right)^{\alpha/2}P^{(0)}\ $ (134)

for all $\alpha\in[0,2]$. As we have seen, the second largest eigenvalue of
$P^{(0)}$, denoted $\mu^{(0)}_{N}$, is given by

$\mu^{(0)}_{N}=\frac{3N-1}{3(N-1)^{2}}\ .$ (135)

It follows from (134) and (135) that for all $f$ orthogonal to the constants,

$\langle
f,P^{(\alpha)}f\rangle\leq\left(\frac{N}{N-1}\right)^{\alpha/2}\frac{3N-1}{3(N-1)^{2}}{\|}f{\|}_{2}^{2}\
,$ (136)

for all $\alpha\in[0,2]$. For $\alpha=2$, we have ${\displaystyle\langle
f,P^{(2)}f\rangle\leq\frac{N(3N-1)}{3(N-1)^{3}}{\|}f{\|}_{2}^{2}}$. Note that

$\frac{N(3N-1)}{3(N-1)^{3}}=\frac{1}{N-1}+\frac{5}{3}\frac{1}{(N-1)^{2}}+\frac{2}{3}\frac{1}{(N-1)^{3}}\
,$

which evidently decreases monotonically as $N$ increases. Next, since
$W^{(2)}(\vec{v})=1-\frac{1}{(N-1)^{2}}$, we have that for all $f$ orthogonal
to the constants

$-\langle f,\widehat{L}_{N,2}f\rangle=\langle
f,(W^{(2)}-P^{(2)})f\rangle\geq\left(1-\frac{1}{N-1}-\frac{8}{3}\frac{1}{(N-1)^{2}}-\frac{2}{3}\frac{1}{(N-1)^{3}}\
\right){\|}f{\|}_{2}^{2}\ .$ (137)

For $N=3$, this yields only the useless bound $-\langle
f,\widehat{L}_{3,2}f\rangle\geq-\frac{1}{4}{\|}f{\|}_{2}^{2}$. But already for
$N=4$, it yields

$-\langle f,\widehat{L}_{4,2}f\rangle\geq\frac{28}{81}{\|}f{\|}_{2}^{2}\ .$

Since the right hand side of (137) increases as $N$ increases, this, together
with the comparison from Lemma 2.18, proves:

######  THEOREM B.1.

For all $N\geq 4$,

$\widehat{\Delta}_{N,2}\geq
1-\frac{1}{N-1}-\frac{8}{3}\frac{1}{(N-1)^{2}}-\frac{2}{3}\frac{1}{(N-1)^{3}}>0\
,$

and for all $\alpha\in(0,2)$,

$\widehat{\Delta}_{N,\alpha}\geq\left(\frac{N-1}{N}\right)^{1-\alpha/2}\left(1-\frac{1}{N-1}-\frac{8}{3}\frac{1}{(N-1)^{2}}-\frac{2}{3}\frac{1}{(N-1)^{3}}\
\right)>0\ .$

At this point, the only estimate we lack for a fully quantitative result is a
quantitative estimate on $\widehat{\Delta}_{3,2}$.

###  B.2 An explicit bound for $N=3$

By what has been explained earlier,
$\widehat{\Delta}_{3,2}=\frac{3}{4}-\nu_{3}$ where

$\nu_{3}=\sup\left\\{\langle f,P^{(2)}f\rangle_{L^{2}({\mathcal{S}}_{3})}\ :\
{\|}f{\|}_{2}=1\ ,\quad\langle f,1\rangle_{L^{2}({\mathcal{S}}_{N})}=0\
\right\\}\ ,$ (138)

and by Lemma 2.17, $\widehat{\Delta}_{3,2}>0$, or, what is the same
$\nu_{3}<\frac{3}{4}$.

If $\nu_{3}\leq\frac{1}{2}$, then evidently
$\widehat{\Delta}_{3,2}\geq\frac{1}{4}$. Therefore, we need only consider the
possibility that $\nu_{3}>\frac{1}{2}$, and as we have seen, in in this case
$\nu_{3}$ is an eigenvalue of $P^{(2)}$, and necessarily
$\nu_{3}<\frac{3}{4}$.

In seeking the second largest eigenvalue of $P^{(2)}$, we need only consider
functions $f$ of the form

$f(\vec{v})=\sum_{j=1}^{N}\varphi(v_{j})$ (139)

or

$f(\vec{v})=\varphi(v_{1})-\varphi(v_{2})\ ,$ (140)

where in the second case we have taken advantage of the the symmetry of
$P^{(2)}$ to assume without loss of generality that $f$ is antisymmetric under
interchange of $v_{1}$ and $v_{2}$.

######  LEMMA B.2.

For $N=3$, the largest eigenvalue of $P^{(2)}$ on the orthogonal complement of
the symmetric sector is no greater than $0.735$. Thus, either
$\widehat{\Delta}_{3,2}\geq 0.015$, or else the gap eigenfucntion is
symmetric.

###### Proof of Lemma B.2.

For later use, we begin the proof for $N\geq 3$, and specialize to $N=3$
later. Let $f$ be given by (140), where we may assume that $\varphi$ is
orthogonal to the constants. Then

$\frac{1}{N}w_{N,2}(v_{1})(1-K)\varphi(v_{1})-\frac{1}{N}w_{N,2}(v_{2})(1-K)\varphi(v_{2})=\lambda(\varphi(v_{1})-\varphi(v_{2})\
.$

Multiplying by $\varphi(v_{1})$ and integrating,

$\frac{1}{N}\int_{{\mathcal{S}}_{N}}w_{N,2}(v_{1})|(1-K)\varphi(v_{1})|^{2}=\lambda\langle\varphi,(1-K)\varphi\rangle\
.$ (141)

By (76),

$w_{N,2}(v)=\frac{N}{N-1}-\frac{N}{(N-1)^{2}}|v_{k}|^{2}\ .$ (142)

Using (142) in (141) yields

$\frac{1}{N-1}\langle\varphi,(1-K)^{2}\varphi\rangle-\frac{1}{(N-1)^{2}}\langle(1-K)\varphi,|v|^{2}(1-K)\varphi\rangle=\lambda\langle\varphi,(1-K)\varphi\rangle\
.$ (143)

Now write $\sqrt{1-K}\varphi=\psi+\zeta$ where $\psi$ is orthogonal to the
constants, the three components of $v$ and $|v|^{2}$. Then $\zeta$ is an
eigenvector of $K$ with eigenvalue $-1/(N-1)$, and hence

$\frac{1}{N-1}\langle\varphi,(1-K)^{2}\varphi\rangle=\frac{1}{N-1}\langle\psi,(1-K)\psi\rangle\
\ ,$ (144)

and

$\displaystyle\langle(1-K)\varphi,|v|^{2}(1-K)\varphi\rangle$ $\displaystyle=$
$\displaystyle\langle\sqrt{1-K}\psi,|v|^{2}\sqrt{1-K}\psi\rangle+\frac{N-2}{(N-1)^{2}}{\|}|v|\zeta{\|}_{2}^{2}$
$\displaystyle-$ $\displaystyle
2{\|}|v|\sqrt{1-K}\psi{\|}_{2}\frac{\sqrt{N-1}}{N-2}{\|}|v|\zeta{\|}_{2}^{2}$
$\displaystyle\geq$
$\displaystyle\left(1-\frac{1}{t}\right)\langle\sqrt{1-K}\psi,|v|^{2}\sqrt{1-K}\psi\rangle+(1-t)\frac{N-2}{(N-1)^{2}}{\|}|v|\zeta{\|}_{2}^{2}\
,$

for all $t>0$, where we have used the arithmetic-geometric mean inequality.

At this point we specialize to $N=3$, and carry out some explicit computations
that could be done for all $N\geq 3$, but are then more cumbersome.

Write $\zeta=\sum_{j=1}^{4}a_{j}\eta_{j}(v)$ where, as before, $\eta(j)v)={\bf
e}\cdot v$ for $j=1,2,3$, and where $\eta_{4}(v)=|v|^{2}-1$. One readily
computes that

${\|}|v|\zeta{\|}_{2}^{2}=\sum_{j=1}^{4}|a_{j}|^{2}{\|}|v|\eta_{j}{\|}_{2}^{2}$
(146)

and that ${\displaystyle\int_{{\mathcal{S}}_{3}}|v_{1}|^{4}{\rm
d}\sigma_{3}=\frac{5}{4}}$ and
${\displaystyle\int_{{\mathcal{S}}_{3}}|v_{1}|^{6}{\rm
d}\sigma_{3}=\frac{7}{4}}$. From here it follows that

${\|}|v|\eta_{j}{\|}_{2}^{2}=\frac{5}{4}\quad{\rm for}\quad j=1,2,3\quad{\rm
and}\quad{\|}|v|\eta_{4}{\|}_{2}^{2}=1\ .$

Using this in (146) finally yields
${\|}|v|\zeta{\|}_{2}^{2}\geq{\|}\zeta{\|}_{2}^{2}$, and evidently
$\langle\sqrt{1-K}\psi,|v|^{2}\sqrt{1-K}\psi\rangle\leq
3{\|}\sqrt{1-K}\psi{\|}_{2}^{2}$. Therefore, for $0<t<1$, we have from (B.2)
that

$\langle(1-K)\varphi,|v|^{2}(1-K)\varphi\rangle\geq
3\left(1-\frac{1}{\lambda}\right){\|}\sqrt{1-K}\psi{\|}_{2}^{2}+(1-\lambda){\|}\zeta{\|}_{2}^{2}\
.$

Using this estimate in (143) yields

$\frac{1}{2}{\|}\sqrt{1-K}\psi{\|}_{2}^{2}+\frac{3}{4}{\|}\zeta{\|}_{2}^{2}-3\left(1-\frac{1}{t}\right){\|}\sqrt{1-K}\psi{\|}_{2}^{2}-(1-t){\|}\zeta{\|}_{2}^{2}\geq\lambda({\|}\psi{\|}_{2}^{2}+{\|}\zeta{\|}_{2}^{2})\
.$

The second most negative eigenvalue of $K$ for $N=3$ is $-\frac{3}{8}$; see
[7, Section 8], where this eigenvalue is denoted $\kappa_{1,2}$. It follows
that ${\|}\sqrt{1-K}\psi{\|}_{2}^{2}\leq\frac{11}{8}{\|}\psi{\|}_{2}^{2}$.
Therefore,

$\left(\frac{1}{16}-3+\frac{3}{t}\right){\|}\psi{\|}_{2}^{2}+\left(\frac{3}{4}-1+t\right){\|}\zeta{\|}_{2}^{2}\geq\lambda({\|}\psi{\|}_{2}^{2}+{\|}\zeta{\|}_{2}^{2})\
.$

Choosing $t=0.985$, we have that
$0.735({\|}\psi{\|}_{2}^{2}+{\|}\zeta{\|}_{2}^{2})\geq\lambda({\|}\psi{\|}_{2}^{2}+{\|}\zeta{\|}_{2}^{2})$.
∎

The remaining task is to bound the second largest eigenvalue of $P^{(2)}$ in
the symmetric sector. We begin considering general $N\geq 3$ and shall
specialize to $N=3$ later.

Let $f$ be given as in (139). Then $P^{(2)}f=\lambda f$ becomes

$\frac{1}{N}\sum_{k=1}^{N}w_{N,2}(v_{k})(\varphi(v_{k})+(N-1)K\varphi(v_{k}))=\lambda\sum_{j=1}^{N}\varphi(v_{j})\
.$ (147)

By Theorem 2.5,

$\frac{1}{N}w_{N,2}(v)(\varphi(v)+(N-1)K\varphi(v))=\lambda\varphi(v)\ .$
(148)

we have that Therefore, multiplying both sides of (148) by $\varphi(v)$ and
integrating, we obtain

$\frac{1}{N}\int_{{\mathcal{S}}_{N}}\varphi(v_{1})w_{N,2}(v_{1})(\varphi(v_{1})+(N-1)K\varphi(v_{1})){\rm
d}\sigma_{N}=\lambda{\|}\varphi{\|}_{2}^{2}\ .$ (149)

By (76), (149) becomes

$\langle\varphi,K\varphi\rangle-\frac{1}{N-1}\int_{{\mathcal{S}}_{N}}\varphi(v_{1})|v_{1}|^{2}K\varphi(v_{1})){\rm
d}\sigma_{N}=\\\
\left(\lambda-\frac{1}{N-1}\right){\|}\varphi{\|}_{2}^{2}+\frac{1}{(N-1)^{2}}\int_{{\mathcal{S}}_{N}}|v_{1}|^{2}\varphi^{2}(v_{1})){\rm
d}\sigma_{N}\ .$ (150)

Define an operator $M$ by $M\phi(v)=|v|^{2}(1+(N-1)K)\phi(v)$ .

Then (150) becomes

$\left(\lambda-\frac{1}{N-1}\right)=\frac{\langle\varphi,K\varphi\rangle-(N-1)^{-2}\langle\varphi,M\varphi\rangle}{{\|}\varphi{\|}_{2}^{2}}\
,$

Thus $\lambda-(N-1)^{-1}$ can be computed by computing the supremum of the
right hand side as $\varphi$ ranges over functions that are orthogonal to $1$,
the three components of $v$ and $|v|^{2}$.

Note that $M$ commutes with rotations so the different angular momentum
sectors are mutually orthogonal, and can be considered separately. In each
sector, by the usual recursion relations for orthogonal polynomials, the
matrix representing $M$ in the eigenbasis of $K$ is tri-diagonal and
explicitly computable, and the bounds proved in [7, Section 8] can be used to
limit the number of angular momentum sectors that need to be considered. Hence
one could obtain explicit bounds this way.

## References

  * [1] P. Caputo, Spectral gap inequalities in product spaces with conservation laws, In Stochastic Analysis on Large Interacting Systems T. Funaki and H. Osada, eds., Math. Soc. Japan, Tokyo. (2004), 53 88.
  * [2] P. Caputo,  On the spectral gap of the Kac walk and other binary collision processes, ALEA Lat. Am. J. Probab. Math. Stat., 4, (2008), 205-222.
  * [3] E. A. Carlen, J. A. Carrillo and M. C. Carvalho:  Strong convergence towards homogeneous cooling states for dissipative Maxwell models, Ann. l’Institut H. Poincare (C) Non Linear Analysis 26, 5, 2009, 1675-1700 in Journes Equations aux derivees partielles, Nantes, 5-9 Juin 2000.
  * [4] E. A. Carlen, M. C. Carvalho and M. Loss:  Many body aspects of approach to equilibrium, in Journes Equations aux derivees partielles, Nantes, 5-9 Juin 2000.
  * [5] E. A. Carlen, M. C. Carvalho and M. Loss:  Determination of the spectral gap for Kac’s master equation and related stochastic evolution, Acta Mathematica 191, (2003) 1-54
  * [6] E. A. Carlen, M. C. Carvalho and M. Loss:  Spectral gap for the Kac model with hard sphere collisions. J. Funct. Anal. 266 (2014), no. 3, 1787-1832
  * [7] E. A. Carlen, J. Geronimo and M. Loss: Determination of the spectral gap in the Kac model for physical momentum and energy conserving collisions, SIAM Jour. Mathematical Analysis 40, no. 1, (2008) 327-364
  * [8] W. Feller, An introduction to Probability Theory and its Applications, Volume I, Wiley and sons, New York, 1950
  * [9] G. Giroux and R. Ferland, Global spectral gap for Dirichlet Kac random motions. J. Stat. Phys. 132, (2008) 561-567.
  * [10] A. Grigo, A. Khanin and D. Szász, Mixing rates of particle systems with energy exchange, Nonlinearity 25 (2012), 2349-2376.
  * [11] E. Janvresse,  Spectral Gap for Kac’s model of Boltzmann Equation, Annals. of Prob., 29 (2001) 288-304.
  * [12] M. Kac, Foundations of kinetic theory, Proc. 3rd Berkeley symp. Math. Stat. Prob., J. Neyman, ed. Univ. of California, vol 3, (1956) 171–197.
  * [13] M. Kac  Probability and Related Topics in Physical Sciences, Interscience Publ. LTD., London, New York (1959)
  * [14] F.G. Mehler  Über die Entwicklung einer Function von beliebig vielen Variablen nach Laplaschen Functionen höherer Ordnungn. Crelle’s Journal 66 (1866), 161–176.
  * [15] S. Mischler and C. Mouhot, Kac’s Program in Kinetic Theory, Invent. Math., 193 (2013) 1-147.
  * [16] M. Sasada, On the spectral gap of the Kac walk and other binary collision processes on $d$-dimensional lattice, In Symmetries, Integrable Systems and Representations, K. Iohara et al., eds. Springer Proc. Math. Stat. 40, 543 560. Springer, Heidelberg, 2013.
  * [17] M. Sasada, Spectral gap for particle systems with degenerate energy exchange rates, Ann. Probab. 43 (2015),1663-1711.
  * [18] G. Szegö, Orthogonal Polynomials, Vol. 23 of A.M.S. Colloquium Series Publications, A.M.S., Rovidence, 1967.
  * [19] A. S. Sznitman, Topics in propagation of chaos. In: École d Été de Probabilités de Saint-Flour XIX 1989. Lecture Notes in Math., 1464, pp. 165 251. Springer, Berlin, 1991.
  * [20] C. Villani, Cercignani’s Conjecture is Sometimes True and Always Almost True, Comm. Math. Phys., 234, (2003) 455–490.