Minimal Soft Lattice Theta Functions

We study the minimality properties of a new type of"soft"theta functions. For a lattice $L\subset \mathbb{R}^d$, a $L$-periodic distribution of mass $\mu_L$ and an other mass $\nu_z$ centred at $z\in \mathbb{R}^d$, we define, for all scaling parameter $\alpha>0$, the translated lattice theta function $\theta_{\mu_L+\nu_z}(\alpha)$ as the Gaussian interaction energy between $\nu_z$ and $\mu_L$. We show that any strict local or global minimality result that is true in the point case $\mu=\nu=\delta_0$ also holds for $L\mapsto \theta_{\mu_L+\nu_0}(\alpha)$ and $z\mapsto \theta_{\mu_L+\nu_z}(\alpha)$ when the measures are radially symmetric with respect to the points of $L\cup \{z\}$ and sufficiently rescaled around them (i.e. at a low scale). The minimality at all scales is also proved when the radially symmetric measures are generated by a completely monotone kernel. The method is based on a generalized Jacobi transformation formula, some standard integral representations for lattice energies and an approximation argument. Furthermore, for the honeycomb lattice $\mathsf{H}$, the center of any primitive honeycomb is shown to minimize $z\mapsto \theta_{\mu_{\mathsf{H}}+\nu_z}(\alpha)$ and many applications are stated for other particular physically relevant lattices including the triangular, square, cubic, orthorhombic, body-centred-cubic and face-centred-cubic lattices.


Introduction
The mathematical justification of crystal's shape and formation is a very difficult problem which has been actively studied (see [19] and references therein). Indeed, the analytic or numerical investigation of static many-particle Hamiltonian's ground states plays a central role in the design of materials (see e.g. [59]), but the large number of critical points as well as the nonlinearities emerging from the corresponding systems make this mathematical investigation very challenging. Thus, a first natural step is to study systems that are already in a periodic order and where the interaction between points is given by a radially symmetric potential. These potentials arise in physics models of matter (see e.g. [47,49]) in the case of the Born-Oppenheimer adiabatic approximation: the electrons effects are neglected, the energy is reduced to the nuclei interactions (see [56, p. 33]), and two-body potentials are the simplest way to express the total potential energy of the system (see [56, p. 945]).
This problem of minimizing a potential energy per point of the form where L is a Bravais lattice (see also Definition 2.4), which is the most simple possible periodic configuration of points, has received a lot of attention, especially in the following cases: Lennard-Jones type potentials [7,8,9,18], Morse potential [10], two-dimensional Thomas-Fermi model for solid [16], Coulombian renormalized energy [15,61,62], completely monotone interaction potentials [7,26,28], Bose-Einstein Condensates [2,53], diblocks and 3-blocks copolymer interactions [24,50], vortices in quantum ferrofluids [51], inverse power laws (Epstein zeta function) [23,31,32,34,35,36,58], and also in more general settings [14,33]. An important mathematical object, which appears to be central in this theory (see e.g. [7,26]), is the lattice theta function. Given a d-dimensional Bravais lattice L, a scaling parameter α > 0 and a point z ∈ R d , we define the lattice theta function and the translated lattice theta function by (1.1) Physically, θ L (α) can be viewed as the Gaussian self-interaction of L and θ L+z (α) as the Gaussian interaction between point z and lattice L. They actually are the energies per point of the so-called Gaussian Core Model (GCM) restricted to lattices. This model was initially introduced by Stillinger [65] and motivated by the Flory-Krigbaum potential between the centers-of-mass of two polymer chains in an athermal solvent [42]. The phase diagram of the three-dimensional GCM has been numerically investigated for example in [69]. Furthermore, two different problems concerning (1.1) appear to be quite natural once α > 0 is fixed: • the minimization of L → θ L (α) among d-dimensional Bravais lattices with the same density; • the minimization of z → θ L+z (α) among vectors z ∈ R d , where L is fixed.
These minimization problems have been studied by many authors, see e.g. [3,9,11,13,25,26,28,33,34,37,38,39,40,41,52,54,60,64], and two of the most significant results are due to Montgomery [52] -who proved the minimality of the triangular lattice for L → θ L (α) among two-dimensional Bravais lattices of any fixed density -and Cohn, Kumar, Miller, Radchenko and Viazovska [28] -who recently proved the minimality among all periodic configurations with the same density of E 8 and the Leech lattice in dimensions d ∈ {8, 24}. Furthermore, in [11], the minimizer of the translated theta function z → θ L+z (α), for a fixed lattice L, has been also proved to be connected to the optimal electrostatic interaction (and more general long-range weighted interaction energies) between periodic distributions of charges located on L, solving a conjecture stated by Born in [21] about the optimality of the rock-salt structure.
Since polymer chains can be seen as soft interpenetrable spheres with an extent of the order of their radius of gyration (see [49,Sect. 3]), we propose a generalization of the periodic GCM to mass interactions, in the same spirit as it has been done in dimension d = 2 in [12] with Knüpfer. We want to study the minimality properties of lattices L and points z for a kind of "soft GCM" (SGCM), where the objects are smeared out. Furthermore, these mass interaction energies can also be viewed as the expectation values of the lattice theta function and translated lattice theta function defined by (1.1) where the position of the lattice points (resp. the position of z) follow a radially symmetric probability distribution µ L (resp. ν z ). We therefore define (see also Definition 2.8) the translated soft lattice theta function by The main goal of this paper is to derive some minimality properties of lattices L and points z for (L, z) → θ µ L +νz (α), where α > 0 is fixed, generalizing our previous work [12] to the Gaussian interaction potential in all dimensions. In particular, the idea is to link the critical points, strict local minima or global minima of (L, z) → θ L+z (α) with those of (L, z) → θ µ L +νz (α).
Minimizing L → θ µ L +ν 0 (α) and z → θ µ L +νz (α) can be interpreted in terms of "defects" in the periodic SGCM. Thus, we want to understand at which scales the type (i.e. the profile of µ, ν) or the size (i.e. the size of the ball containing almost all the mass of the measures) of the defects do not play any role in these minimization problems. Indeed, many defects appear in perfect crystals and they give generally to the material its properties (corrosion resistance, softness, thermal expansion, etc.). In Solid-State Physics, two kinds of point defects are important (see e.g. Kaxiras [48,Chap. 9]): 1. the extrinsic defects, such as a substitutional impurity, corresponding to z = 0 in our model.
An atom in a perfect crystal is substituted by another one of the same kind (µ = ν) or of a different kind (µ = ν). This impurity is usually chemically similar to the crystal's atom, with a similar size. We are looking for the minimizer of L → θ µ L +ν 0 (α), i.e. the Gaussian interaction between the mass ν 0 centred at the origin and all the masses µ L centred at lattice sites, among a class of lattices with the same density. This also includes the energy per point of the perfect crystal itself when µ = ν (see the top line of Figure 1); 2. the intrinsic defects, such as an interstitial defect, corresponding to z = 0 in our model. An additional atom is located somewhere in the unit cell of the crystal (but not at a lattice site) and is generally smaller and chemically different than the crystal's atoms. We are looking for the minimizer of z → θ µ L +νz (α), i.e. the Gaussian interaction between the mass ν z centred at z and all the masses µ L centred at lattice sites (see the bottom line of Figure 1). Figure 1: Example of different kinds of defects. The lattice is L = Z 2 and the masses are all Gaussian, with different variances. We have chosen to represent only the primitive cell Q Z 2 . On the top: the case without defect µ = ν and z = 0 (left), and the case of an extrinsic defect µ = ν and z = 0 (right). On the bottom: the case of a intrinsic defect µ = ν and z = 0.
Using the methods developed in [12] extended to the d-dimensional case and to the translated lattice theta function, we show that for any measures µ, ν and any α > 0, we can rescale the measures around the points by sufficiently small factors ε and δ, getting two new measures µ ε and ν δ , such that any strict local or global minimality result which is true for (L, z) → θ L+z (α) is also true for (L, z) → θ µ ε L +ν δ z (α). We will say that the minimality is true at a low scale. In other words, there is a difference of scale between the lattice spacing and the radius of the balls where most of the masses are concentrated such that the measures (or the "defects") do not play any role in terms of strict local or global minima. Furthermore, if µ and ν have densities of the form x → ρ(|x| 2 ) where ρ is the Laplace transform of a nonnegative Borel measure (ρ is a completely monotone function), then the optimality occurs at all scales, i.e. for all such measures. As in [12], our results are purely qualitative: we do not give any value of ε and δ such that the properties hold. However, it is interesting to see that, once (ε, δ) are below some threshold values (ε 0 , δ 0 ), then the desired minimality occurs whatever the quotient ε/δ is. Furthermore, a numerical investigation have been made in [12,Remark 14] for the triangular lattice case with uniform measures on disks, showing that the strict local minimality of the triangular lattice proved for small values of the parameters (ε, δ) certainly does not hold when these parameters are too large.
We notice that, as recalled in [12], this type of minimization problem involving smeared out particles appears in many physical and biological systems as condensed matter theory [46], quantum physics models [17], diblock copolymer systems in the low volume fraction limit [55], magnetized disks interactions [44] and swarming or flocking related models [4,5,22,67].
We have also devoted a complete section of the paper to the applications of our results. They mainly are corollaries of minimality results obtained in previous works in the point case µ = ν = δ 0 , i.e. for the lattice theta function and the translated lattice theta functions defined by (1.1), for some particular physically relevant lattices. We thought it was a good opportunity to review all those results in this paper in order to know what are the main open problems associated to these soft theta functions. Furthermore, we show that the minimum of z → θ H+z (α), where H is a honeycomb lattice (see (2.6) and Figure 4), is the centre of an honeycomb of H, and the same is therefore true for the smeared out cases previously stated.
In terms of generalization, it appears to be straightforward that all the results in this paper can be proved for more general energies of the form where f is a L 1 completely monotone summable function, as the one we have studied in [12]. Since our original goal was to study a new kind of theta functions that could have other applications in Number Theory and Mathematical Physics, we did not extend our results to these types of energies. The reader can refer to [12] for details.
Plan of the paper. After giving the definition of lattices, energies, measures and minimality at a low scale and at all scales in Section 2, we show some preliminary results in Section 3, including the generalized Jacobi transformation formula. Our main results are stated and proved in Section 4 and many applications are given in Section 5, where the minimality of the primitive honeycomb's center in the honeycomb lattice case is also proved.

Definitions
We start by defining the space of Bravais lattices with a fixed density as well as their unit cells and the notion of dual lattice. We call (e i ) 1≤i≤d the orthonormal basis of R d , |.| the euclidean norm on R d , u·v the associated scalar product of u, v ∈ R d and B r the closed ball of radius r > 0 and centred at the origin. Furthermore, we write M d (R) the space of n × n matrices with real coefficients. The unit cell (of volume 1) of such Bravais lattice L is defined by We also recall the definitions of the following important lattices belonging to L • d (see also Figure  2): The triangular lattice Λ 1 : These lattices are physically relevant because they correspond to the main crystal structures that exist in nature (see [48, p. 8-9]) which are Bravais lattices. Notice for instance that among the 118 known elements, 21 (resp. 26) have a BCC (resp. FCC) structure. We will also consider the honeycomb lattice H defined by and we call H one of its primitive hexagons (see Figure 4). This is a typical example of periodic configuration, i.e. an union of translated Bravais lattices, that arises in Physics (e.g. as the structure of a graphene sheet). Furthermore, we define D 4 , E 8 , D + d and the Leech lattice as in [30]. It turns out that these lattices have many interesting properties related to energy minimization or packing (see for instance [27,28,33,34,64,68]).
As explained in [66,Sect. 1.4] (see also [14, p. 14]), any lattice L ∈ L • d can be parametrized by a pointL = ( Thus, as in [12], the metric on D d is chosen as the euclidean metric on R n d , d(L, Λ) denotes the euclidean distance between two lattices L, Λ ∈ D d and B r (L) denotes the open ball of radius r centred at L. We write is defined as the n d × n d real matrix of second derivatives with respect to x 1 , ..., x n d . D 3 E is correspondingly the tensor of all third derivatives. The notions of strict local minimizer and critical point in D d is defined as follows: Let us now focus on the interaction potentials we want to consider throughout this paper. Even though the main interaction potential will be the Gaussian one, it turns out that we will use many properties of more general lattice energies. We define two classes of functions that can be written as the Laplace transform of a measure.
Furthermore we say that f ∈ CM d if f ∈ F d and µ f is nonnegative.
Furthermore, the energy per point of any L ∈ D d interacting through a radial potential f ∈ F d is defined by an absolutely convergent sum as follows.
Because our goal is to study masses interactions, we need to specify what kind of measures we are working with. We note P(R d ) the space of probability measures on R d and P r (R d ) the space of probability measures on R d that are rotationally symmetric with respect to the origin. Furthermore, we define the following subspace of P r (R d ): Remark 2.5 (Completely monotone functions). The notations CM d and "cm" mean that ρ (or f in Definition 2.3) is a completely monotone function, i.e. (−1) k ρ (k) (t) ≥ 0 for any t > 0 and any k ∈ N, which is indeed equivalent, by Hausdorff-Bernstein-Widder Theorem [6], for ρ to be the Laplace transform of a nonnegative Radon measure µ ρ .
We now define the periodic measure µ L that corresponds to the union of measures µ centred at all the points of a Bravais lattice L. Definition 2.6 (Periodized measure). For any L ∈ D d and any µ ∈ P(R d ), the periodized measure µ L is defined by In Figure 3, we have represented two different kinds of measures µ L for µ being in P r (R 2 ) and P cm r (R 2 ).
Furthermore, because we want to show the optimality of the previously defined lattices where the masses are sufficiently concentrated around the lattice sites, we define the rescaled measure as follows.
Definition 2.7 (Rescaled measure). Let d ≥ 1. For any µ ∈ P(R d ) and any ε > 0, the rescaled measure µ ε is defined, for any measurable set F ⊂ R d , by We write µ ε L and µ ε z the corresponding rescaled measures of µ L and µ z defined by (2.8).
We finally define the main energies that we are studying in this paper: the translated soft lattice theta function. It corresponds to the Gaussian interaction energy of the measure µ L with itself or with an other measure ν z located at z ∈ R d , and are therefore a generalization of the lattice and translated lattice theta functions θ L (α) and θ L+z (α) defined by (1.1).
Definition 2.8 (Translated soft lattice theta functions). Let d ≥ 1. For any L ∈ D d , any µ, ν ∈ P(R d ), any z ∈ R d and any α > 0, we define the translated soft lattice theta function of µ L by the measure ν z by Remark 2.9 (The honeycomb lattice case). Since the honeycomb lattice H = Λ 1 ∪ (Λ 1 + u) is a union of two translated triangular lattice, it is easy to compute its Gaussian energy per point. More precisely, we have Remark 2.10 (Special cases). We notice that: • For any µ ∈ P(R d ), any L ∈ D d (or any periodic configuration) and any α > 0, it turns out that the self-interaction of µ L is θ µ L +µ 0 (α) (i.e. when ν = µ and z = 0).
• The translated soft lattice theta function for the rescaled measures µ ε and ν δ is given by The goal of this work is to study minimization problems for L → θ µ L +ν 0 (α) in D d -which includes the self-interaction of µ L -and z → θ µ L +νz (α) in Q L for fixed L, at different scales. We therefore precise what we mean by being critical or minimal at a low scale and at all scales.
Definition 2.11 (Criticality and minimality at all scales and at a low scale). Let d ≥ 1 and µ, ν ∈ P(R d ), then: Furthermore, let L be fixed, then we define the same notions of critical point and (strict local) minimality for z → θ µ L +νz (α) of z 0 in Q L at all scales or at a low scale.

Preliminaries
It is well-known that the lattice theta function satisfies the following identity, called Jacobi Transformation Formula (see e.g. [45,Thm. A] or [20] for a general formula involving harmonic polynomials): for any L ∈ D d , any z ∈ R d and any α > 0, We generalize this formula to mass interaction in the following result.
Proposition 3.1 (Generalized Jacobi Transformation Formula). For any d ≥ 1, any L ∈ D d , any µ, ν ∈ P r (R d ), any z ∈ R d and any α > 0, we have

2)
where, for any m ∈ P r (R d ), and ψ m is the Lebesgue-Stieltjes measure of t → m(B t ), i.e. ψ m ([r 1 , r 2 )) = m(B r 2 ) − m(B r 1 ) and J β is the Bessel function of the first kind.
Proof. By the classical Jacobi Transformation Formula (3.1) and Fubini's Theorem, we get wherem is the notation for the Fourier transform of a measure m ∈ P r (R d ). We now recall that m is given by the Hankel-Stieltjes transform (see e.g. [29, Section 2]), i.e. for any x ∈ R d , where ψ m is the Lebesgue-Stieltjes measure of t → m(B t ), and the proof is completed. ) and the fact that ψ m ∈ P(R + ), we obtain We furthermore notice that for d = 2, z = 0 and µ = ν, we recover the formula proved in [12,Prop. 7], up to a factor 2 in the argument of J d 2 −1 that we have corrected and which does not influence the final result.
We now give some basic facts that will be used in the proofs of our main results.
Proof. If f and g have both this representation as the Laplace transform of a measure, the formula [57,Thm. 5.3.11]) shows that f g has also this representation. Moreover, we naturally have |f (r)g(r)| = O(r −d/2−σ ) for some σ > 0 as r → +∞, because it is the case for both functions.
The next lemma is a simple consequence of the Jacobi Transformation Formula (3.1). The following results come from the simple fact that, for any f ∈ F d , any L ∈ D d and any z ∈ R d , as explained e.g. in [7,Sect. 3.1].

Main results
In this part, we generalize the results of [12] about the (strict local) minimality of a lattice for L → θ µ L +ν 0 (α), for a given α > 0, to any dimension. Furthermore, using the same strategy based on an approximation of the quantity we sum by a completely monotone potential, we show the same kind of results for z → θ µ L +νz (α) where L is fixed.
It is important to notice that, by the generalized Jacobi Transformation Formula (3.2), if z = 0, then θ µ L +ν 0 (α) is the sum of a radial potential over L * , i.e. an energy of type E f [L * ] for some f , as defined in Definition 2.4. Therefore, as shown for instance in [8,9,33], any lattice L such that L * has enough symmetries is a critical point of the soft lattice theta function, e.g. Λ 1 , D 3 , D *  Let us prove the second part of the proposition. Since z 0 is a critical point of z → θ L+z (β) for all β > 0, we obtain, using Jacobi transformation formula (3.1) and computing the derivative with respect to z i , We recall that, by Proposition 3.1, the soft lattice theta function can be written as We now use the fact that there exists a Borel measure µ h such that h = L[µ h ] ∈ F d as explained above. Therefore, integrating (4.1) against µ h , where t = π/β is the variable of integration, gives ∀i ∈ {1, ..., d}, q∈L * h(|q| 2 )q i sin(2πq · z 0 ) = 0, for any α > 0. It follows that ∂ z i θ µ L +νz 0 (α) = 0 for all i ∈ {1, ..., d} and all α > 0, i.e. z 0 is a critical point of z → θ µ L +νz (α) in Q L for all α > 0. The next result generalizes [12,Prop. 11] to L → θ µ L +ν 0 (α) and z → θ µ L +νz (α) by giving a sufficient condition for the strict local minimality of a lattice or a point, which are already minimal for the theta functions (1.1), in the case of masses interactions. in D d for some α 0 > 0, then L 0 is a strict local minimizer of L → θ µ L +ν 0 (α 0 ) in D d at a low scale.
2. Let L ∈ D d . If z 0 is a critical point of z → θ L+z (α) in Q L for any α > 0 and a strict local minimizer of z → θ L+z 0 (α 0 ) for some α 0 > 0, then z 0 is a strict local minimizer of z → θ µ L +νz (α 0 ) on Q L at a low scale.
The strict local minimality of L 0 and z 0 also holds at all scales if µ, ν ∈ P cm r (R d ).
Proof. Let us prove the first part of the theorem. The proof is actually a straightforward generalization of [12,Prop. 11]. According to (3.2), it is equivalent to show the strict local minimality of where, for any measure m ∈ P r (R d ), and it is easy to check that where j 1 , j 2 are independent of ε, δ and j ε,δ is at least of order min{ε, δ} 2 . Since ψ µ and ψ ν are both probability measures, the expansion of the second derivative is straightforward. Therefore, the result follows by continuity of (ε, δ) → A ε,δ [L 0 ] and by the fact -following from the boundedness of the Bessel functions J m -that L → D 3 E h ε,δ [L] is bounded on any ball centred at L 0 , independently of (ε, δ) ∈ [0, 1] 2 .
Let us prove the second part of the theorem. Since z 0 is a critical point of z → θ L+z (α) in Q L for all α > 0, we have, by Proposition 4.1, that z 0 is a critical point of z → θ µ L +νz (α 0 ). By Proposition 3.1, we can write, as above, As explained in the proof of the first part of the theorem, by the Taylor expansion of J d/2−1 , it is straightforward to prove that where A ε,δ,z 0 [L] ∈ M d (R) and sup z∈Q L A ε,δ,z [L] → 0 as (ε, δ) → (0, 0), for any chosen norm . on M d (R). Since z 0 is a strict local minimizer of z → θ L+z (α 0 ), it follows that D 2 z θ µ ε L +ν δ z 0 (α 0 ) is positive definite for ε and δ sufficiently small, i.e. there exists ε 0 , δ 0 such that for any 0 ≤ ε < ε 0 and any 0 ≤ δ < δ 0 , z 0 is a strict local minimizer of z → θ µ ε L +ν δ z (α 0 ).
The next result is a generalization of [12,Thm. 2 and 3] in arbitrary dimension to our translated lattice theta function and gives a sufficient condition for the global minimality of a lattice in D d or a point in Q L . Theorem 4.5 (Global minimality). Let d ≥ 1 and µ, ν ∈ P r (R d ), then we have: 1. If L 0 is the unique global minimizer and a strict local minimizer of L → θ L (α) in D d for all α > 0, then, for any α 0 > 0, L 0 is the unique global minimizer of L → θ µ L +ν 0 (α 0 ) in D d at a low scale.
2. Let L ∈ D d . If z 0 is a global minimizer and a strict local minimizer of z → θ L+z (α) in Q L for all α > 0, then, for any α 0 > 0, z 0 is a global minimizer of z → θ µ L +νz (α 0 ) in Q L at a low scale.
The global minimality of L 0 and z 0 also hold at all scales if µ, ν ∈ P cm r (R d ).
Proof. Let us prove the first part of the theorem and let us start with the case µ, ν ∈ P cm r (R d ). If L 0 is the unique global minimizer of L → θ L (α) for all α > 0, then, by Lemmas 3.4 and 3.6, L 0 = L * 0 is the unique global minimizer of any lattice energy of the form E f where f ∈ CM d . Furthermore, it has been proved in [12,Lem. 10] that µ ∈ P cm r (R d ) ⇐⇒μ ∈ P cm r (R d ). Therefore, as the set of completely monotone functions is stable by product, we haveμν ∈ P cm r (R d ) and it follows from the complete monotonicity of t → e −βt for any β > 0 that h : is completely monotone. Therefore, L 0 is the unique global minimizer of L → θ µ L +ν 0 (α 0 ) for any fixed α 0 > 0.
For µ, ν ∈ P r (R d ), the proof is again a generalization of [12,Thm. 2]. We first remark that any minimizer of L → θ µ ε L +ν δ 0 (α 0 ) belongs to a ball of center L 0 with a finite radius. Indeed, this fact is proved in [12,Lemma 12] in dimension d = 2 and is directly generalizable to any dimension by bounding below |p| 2 , p ∈ L, by |p | 2 where p ∈ L λ = λ −1 Z ⊕ Λ for some λ ≥ 1 which is also a parameter of L, once the lattice is parametrized by n d parameters as explained in Section 2, and Λ ∈ L • d−1 . We therefore get e −πα 0 |λ −1 (m,0,...,0)+x−y| 2 dµ ε (x)dν δ (y) → +∞ as λ → +∞. The same can be done in all the unbounded directions of the fundamental domain D d , proving that the global minimizer necessarily belongs to a compact subset of D d , for example a closed ball containing L 0 that we will note K.
Definingh ε,δ (r) :=k ε,δ (r)e − πr α 0 , we then have by the exponential decay of the lattice theta function and where C does not depend on the lattice L ∈ K. Therefore, (4.3) is proved. For the second inequality, we have for a constant C independent of ε and δ. Indeed, it is a consequence of the fact thath ε,δ = L[m ε,δ ] and m ε,δ has a compact support and then follows by the strict local minimality of L 0 for L → θ L (α) for all α > 0, which implies the same strict local minimality for any E f where f ∈ CM d (see Lemma 3.6), and in particular for E h ε,δ (see [12,Section 2.2] for details).
Let us prove the second part of the theorem and let us begin again with µ, ν ∈ P cm r (R d ). The proof uses the same ingredient as the previous one. Indeed, if z 0 is a global minimizer of z → θ L+z (α) in Q L for all α > 0, then by Lemma 3.6 the same holds for where f ∈ CM d . Let α 0 > 0, then by the generalized Jacobi Transformation Formula (3.2), we have, for some constant C α 0 ,d > 0, for a completely monotone function h ∈ CM d , as explained previously. We therefore obtain, by Poisson Summation Formula (see e.g. [45, Appendix A]), Since h is completely monotone, thenĥ is completely monotone by [12,Lem. 10] and therefore z 0 is a global minimizer of z → p∈Lĥ (|p + z| 2 ) by Lemma 3.6, which concludes the proof.
Let us now consider µ, ν ∈ P r (R d ). For convenience, we define, for fixed α 0 > 0 and L ∈ D d , where h ε,δ is defined by (4.2). We again recall that minimizing z → θ µ L +νz (α 0 ) is equivalent with minimizing F h ε,δ . We remark that, for h ε,δ defined above, we have, for any z ∈ Q L , for some constant C > 0. The second inequality follows from the strict local minimality of z 0 for z → θ L+z (α 0 ). Therefore, for any z ∈ Q L , Thus, for any z ∈ Q L such that F h ε,δ (z) ≤ F h ε,δ (z 0 ), this implies |z − z 0 | ≤ C max{ε, δ} which is not true for sufficiently small ε and δ by the strict local optimality of z 0 in Q L previously shown in Theorem 4.3.
Remark 4.6 (Difference between ε and δ). An interesting property is the fact that ε and δ can be chosen with different scales -for instance ε = 1/n and δ = 1/ √ n -as long as they are below the critical values ε 0 and δ 0 .
Remark 4.7 (Global minimizer of the lattice theta function). Notice that the minimizer of L → θ L (α) in D d is known to be the same for all α > 0, so far, only in dimension d ∈ {2, 8, 24} as proved in [28,52]. In these dimensions, the minimizers are, respectively, the triangular lattice Λ 1 , E 8 and the Leech lattice. 2. For any d ≥ 1, we proved in [13,Prop. 1.3] (see also [11,Prop. 3.7]) that the center c a = 1 2 (a 1 , ..., a d ) of the primitive cell Q Z d a of the orthorhombic lattice Z d a defined by (2.3) is the unique global minimizer of z → θ Z d a +z (α) for all α > 0. Applications of these results will be shown in the next section, as well as the optimality of the center of the primitive hexagon in the honeycomb lattice case, using Baernstein's theorem.
Furthermore, if z 0 is a global minimizer of z → θ L+z (α) for all α > 0, then it is the case as α → +∞ and it turns out that z 0 is necessarily a deep hole of L, i.e. a solution of the following optimization problem: max

Applications to particular lattices
We finally apply the previous results to some particular lattices defined by (2.1)-(2.4) as well as D 4 , E 8 , the Leech lattice and D + d . This is the perfect opportunity to recall the main results that are currently known about the local and global minima of the (translated) lattice theta functions, which are now generalized to the masses interactions case.
Montgomery [52] proved the minimality of Λ 1 in D 2 for L → θ L (α) for all α > 0. Furthermore, as recalled in the previous section, Baernstein [3] proved the minimality of the two barycenters of the primitive triangles composing Q Λ 1 for z → θ Λ 1 +z (α) also for all α > 0. Therefore, applying Montgomery's theorem, Baernstein's theorem, and our main results, we extend the minimality of Λ 1 and its barycenters to Gaussian masses interactions. Notice that the µ = ν, z = 0 case has been already proved in [12].
Remark 5.2 (Importance of the triangular lattice). The triangular lattice arises in many physical models such as Bose-Einstein Condensates [2], Superconductivity [62], Coulomb Gases [63] or diblock copolymer interaction [24]. We also recall that Λ 1 is conjectured (see Cohn and Kumar [26,Conjecture 9.4]) to be the unique minimizer of the lattice theta function among periodic configurations (not only among Bravais lattices) of fixed density. Furthermore, Λ 1 can also be viewed as a layer of a FCC or BCC lattice potentially shifted by a vector parallel to z 1 or z 2 , as explained in [13,Sect. II.2]. Then, as already discussed in [9], Corollary 5.1 is of great interest for the understanding of BCC and FCC stability in the smeared out particle case, using dimension reduction techniques as in [13].
We now show the following result, using Baernstein's theorem [3, Thm. 1], for the honeycomb lattice H. 3 . Then, for any µ, ν ∈ P r (R 2 ), z 0 is a critical point of z → θ µ H +νz (α) in H at all scales. Furthermore, z 0 is the unique global minimizer of z → θ µ H +νz (α) in H at a low scale. Moreover, if µ, ν ∈ P cm r (R 2 ), then the minimality holds at all scales.
Proof. It is sufficient to show the µ = ν = δ 0 case and to apply the second part of Theorem 4.5 which remains true in this non-Bravais case by simply following the lines of its proof and using the fact that θ H+z (α) = 1 2 (θ Λ 1 +z (α) + θ Λ 1 +u+z (α)) .
We therefore have proved that θ H+z (α) ≥ θ H+z 1 (α), and z 1 is the unique global minimizer of z → θ H+z (α) in H for all α > 0, by symmetry. Figure 4: Minimizer z 0 of z → θ H+z (α) in a primitive hexagon H. The black dots represent Λ 1 and the grey one represent Λ 1 + u.
Remark 5.4 (Gaussian interaction between Lithium atom and Graphene sheet). The Gaussian interaction between a mass and a honeycomb structure is physically relevant. For instance, in [43], the authors have been designed a Gaussian approximation potential modelling the interaction energy between a Lithium atom and a Graphene sheet structure. This has been done by applying Machine Learning to Density Functional Theory reference data. Our result gives the exact location of the Lithium atom minimizing the Gaussian interaction with a Graphene sheet.
Because of the particular role of the cubic lattices Z d and their orthorhombic deformations Remark 5.8 (Conjecture for D 3 and D * 3 in the soft case). We also notice that, according to Sarnak and Strömbergsson [64,Eq. (43)], D 3 (resp. D * 3 ) is expected to be a global minimizer of L → θ L (α) in D 3 for any α ≥ 1 (resp. α ≤ 1). Corollary 5.7 and the fact that the only three-dimensional "volume-stationary" lattices are Z 3 , D 3 and D 3 * (see [10,Section 3]) supports this conjecture for the soft lattice theta function.
Remark 5.10 (Deep holes and BCC/FCC lattices). For L ∈ {D 3 , D * 3 }, the global minimizers of z → θ L+z (α) in Q L are expected to be the deep holes of L, solution of (4.5), i.e. z 0 = 2 − 1 3 (1, 1, 1) for D 3 and z 1 = 2 − 2 3 (1, 1, 0) for D * 3 as well as all their images in Q L by symmetry. These locations are the usual one where a different atom can be added to the structure, in order to create, for example, an ion. They are also called octahedral sites, as for the cubic lattice (see Remark 5.6).
Finally, the same kind of local results for L → θ µ L +ν 0 (α) can be stated for some other special lattices, based on [33,34,64]. 2. For all odd integer d ≥ 9, there exists α d such that for any α 0 > α d , D + d is a strict local minimizer of L → θ µ L +ν 0 (α 0 ) on D d at a low scale if µ, ν ∈ P r (R d ).
3. E 8 and the Leech lattice are global minimizers of L → θ µ L +ν 0 (α 0 ) on D d , for any α 0 > 0, at a low scale if µ, ν ∈ P r (R d ) and at all scales if µ, ν ∈ P cm r (R d ), for the corresponding dimensions d ∈ {8, 24}