Expected Riesz energy of some determinantal processes on flat tori

We compute the expected Riesz energy of random points on flat tori drawn from certain translation invariant determinantal processes and determine the process in the family providing the optimal asymptotic expected Riesz energy.


Introduction
Our objective is to study the asymptotics of the expected Riesz energy of certain point processes (random finite point configurations) in a flat torus Ω ⊂ R d . If Λ is a lattice in R d (i.e. Λ = AZ d for some nonsingular square matrix A) we identify the fundamental domain where the column vectors of A = [v 1 , . . . , v d ] is a Z basis of the lattice Λ, with the flat torus R d /Λ.
For the sphere S 2 , the authors in [1] estimate asymptotically the expected energy of points of the, so-called, spherical ensemble. In [3], the authors study the harmonic ensemble in S d and prove, in some cases, the optimality of the expected asymptotic energy of this process among rotation invariant determinantal processes. In both cases, the expected asymptotic energy was used to get upper bounds for the minimal Riesz energy. Here, we study also the optimality of the expected asymptotic energy among a collection of determinantal processes invariant under translations and it turns out that the best process can be found as an easy consequence of Riesz's rearrangement inequality.
This provides explicit examples with the lowest energy bounds known on the torus in high dimensions. zeta function for Λ is defined, for s > d, as Observe that ζ Λ (s; x) is the Λ-periodic potential generated by the Riesz s-energy |x| −s .
When s ≤ d the sum above is infinite for all Then, F s,Λ (x) is an entire function of s and therefore by the relation we obtain a analytic continuation of ζ Λ (s; x) to s ∈ C \ {d}. Observe that the function 1/Γ(s) is entire and that all the sums in (1) converge uniformly. We are interested in the range 0 < s < d.
We denote as X a (simple) random point process in a compact set Ω ⊂ R d . And let µ be the normalized Lebesgue measure. A way to describe the process is to specify the random variable counting the number of points of the process in D, for all Borel sets D ⊂ Ω. We denote this random variable as X (D).
These point processes are characterized by their joint intensity functions ρ k in Ω k satisfying that for any family of mutually disjoint subsets D 1 , . . . , D k ⊂ Ω. We assume that ρ k (x 1 , . . . , x k ) = 0 when x i = x j for i = j.
A random point process is called determinantal with kernel K : Ω × Ω → C, if it is simple and the joint intensities with respect to a background measure µ are given by for every k ≥ 1 and x 1 , . . . , x k ∈ Ω.
To define the processes we will consider only projection kernels.
Definition 1. We say that K is a projection kernel if it is a Hermitian projection kernel, i.e. the integral operator in L 2 (µ) with kernel K is self-adjoint and has eigenvalues 1 and 0.
By Macchi-Soshnikov's theorem [8,Theorem 4.5.5], a projection kernel K(x, y) defines a determinantal process and it has N points almost surely if the trace for the corresponding integral operator equals N , i.e. if Observe that the random vector in Ω N generated with density is a determinantal process with the right marginals i.e. the joint intensities are given by determinants of the kernel [2, Remark 4.2.6].
Given now a function f : Ω × Ω → [0, ∞) it is easy to compute the expected pair potential energy, [8, Formula (1.2.2)]: Proposition 1. Let K(x, y) be a projection kernel with trace N in Ω and let ω = (x 1 , . . . , x N ) ∈ Ω N be N random points generated by the corresponding determinantal point process. Then, for any measurable f :

Flat torus.
In our setting we take as Ω ⊂ R d the flat torus R d /Λ, for some lattice Λ with dual Λ * .
To construct the kernel we consider for w ∈ Λ * , the Laplace-Beltrami eigenfunctions f w (u) = e 2πi u,w of eigenvalue −4π 2 w, w . Then has compact support and we define the kernels and the corresponding determinantal point processes on the flat torus Ω. For these processes (we are not going to distinguish between the integral operator defined by the kernel K N and the kernel itself) we get points almost surely.
is the Dirichlet kernel.

1.3.
Some known results about minimal periodic Riesz s-energy. It was shown in [7] that for 0 < s < d there exists a constant C s,d independent of Λ such that for N → ∞ The constant C s,d above is not known (unless d = 1). In [7] the authors found an upper bound in terms of the Epstein zeta function. Recall that for a lattice Λ ⊂ R d , the Epstein zeta function ζ Λ (s) defined by can be extended analytically (as in (1)) to C \ {d}. One can see easily that ζ Λ (0) = −1 and the residue of where Λ runs on the lattices with |Λ| = 1. It has been conjectured (see [4]) that if d = 2, 4, 8, or 24, then C s,d = ζ Λ d (s) where Λ d denotes (respectively) the hexagonal lattice, the D 4 lattice, the E 8 lattice and the Leech lattice (scaled to have |Λ d | = 1). When d = 1 indeed C s,1 = ζ Z (s) = 2ζ(s). For d = 2 it is known, due to the work of several authors, that inf Λ ζ Λ (s) is attained for the triangular lattice, see [11] where the result is deduced from the corresponding result for theta functions. It is observed in [12] that from Siegel's integration formula it follows that where dλ d is the volume measure in the space of lattices, [14, p. 172]. One deduces then that C s,d < 0, although for large dimensions there are no examples providing negative bounds. Indeed, from [12], see also [13,Theorem 1], all explicitly known lattices in large dimensions are such that the corresponding Epstein zeta functions have a zero in 0 < s < d, i.e. the analogue of the Riemann hypothesis fails for Epstein zeta functions, see Remark 2.

Expected energies
By Proposition 1 the expected periodic Riesz s-energy of t N = tr(K N ) random points ω = (x 1 , . . . , x t N ) drawn from the determinantal process defined by the kernel It is easy to see, [6], that for 0 < s < d , and therefore by translation invariance Our first result is a nice closed expression for the integral above.
Theorem 1. Let ω = (x 1 , . . . , x t N ) be drawn from the determinantal process on the flat torus R d /Λ given by the kernel Then, for 0 < s < d, Remark 1. Observe that for N random points chosen independently and uniformly in Ω (i.e. for the Poisson point process) the expected energy is given by so the improvement (lowering) in the determinantal case comes from the last summand above. Therefore, to get a good upper bound for the minimal energy we want to maximize the sum This is not an easy task in general. For example, when t N = 2 this would lead to find the shortest non-zero vector in the lattice Λ * i.e. m(Λ * ) = min{|w| : w ∈ Λ * \ {0}}, or equivalently, the density of the densest lattice sphere packing.

Proof. To compute the integral in (3) we write
and using the expression for F s,Λ (u), where the sums converge uniformly, we get v∈Λ Ω Observe that Ω e 2πi u,w−w +η dµ(u) = δ w −w,η , Putting all together, and using that for s < d This last integral converges for all s < d, and using that (for α < −1) we get the result.
Now we define a way to get different invariant kernels (by choosing different sequences of functions κ) and we estimate the corresponding expected energies.  Then, for 0 < s < d, if ω = (x 1 , . . . , x t N ) ∈ Ω t N are t N points drawn from the determinantal process defined by κ D,Λ The measure ν is a multiple of the Lebesgue measure such that if Ω * is a fundamental domain for Λ * then ν(Ω * ) = 1.
Proof. From Theorem 1 we get that For the other term, A natural question is now, given a fixed lattice Λ, to find the optimal D ⊂ R d in the definition of κ D,Λ . By the result above, all we have to check is what is the domain giving the larger potential I D ν . The result follows from Riesz rearrangement inequality [10, p. 87].
Theorem 2 (Riesz rearrangement inequality). Given f, g, H nonnegative functions in R d with h(x) = H(|x|) symmetrically decreasing. Then where f * , g * are the symmetric decreasing rearrangements of f and g.
We apply the result for f = g = χ D getting f * = g * = χ D * where D * is the open ball centered at the origin with |D * | = |D|.

2.
1. An upper bound for the minimal periodic Riesz s-energy. It is clear that E s,Λ (N ) ≤ E ω∈(R d ) N (E s,Λ (ω)), and therefore from the results above we get that C s,d in (2) satisfies, for 0 < s < d, where D is any bounded domain such that |D||Λ| = 1. By the discussion above the best choice is to take the set D to be a ball.  Proof. The proof is an easy computation. We take the normalization of the Fourier transformf