Minimal Soft Lattice Theta Functions

We study the minimality properties of a new type of “soft” theta functions. For a lattice L⊂Rd\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L\subset {\mathbb {R}}^d$$\end{document}, an L-periodic distribution of mass μL\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _L$$\end{document}, and another mass νz\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu _z$$\end{document} centered at z∈Rd\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$z\in {\mathbb {R}}^d$$\end{document}, we define, for all scaling parameters α>0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha >0$$\end{document}, the translated lattice theta function θμL+νz(α)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta _{\mu _L+\nu _z}(\alpha )$$\end{document} as the Gaussian interaction energy between νz\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu _z$$\end{document} and μL\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _L$$\end{document}. We show that any strict local or global minimality result that is true in the point case μ=ν=δ0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu =\nu =\delta _0$$\end{document} also holds for L↦θμL+ν0(α)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L\mapsto \theta _{\mu _L+\nu _0}(\alpha )$$\end{document} and z↦θμL+νz(α)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$z\mapsto \theta _{\mu _L+\nu _z}(\alpha )$$\end{document} when the measures are radially symmetric with respect to the points of L∪{z}\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L\cup \{z\}$$\end{document} 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 H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathsf {H}}$$\end{document}, the center of any primitive honeycomb is shown to minimize z↦θμH+νz(α)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$z\mapsto \theta _{\mu _{{\mathsf {H}}}+\nu _z}(\alpha )$$\end{document}, and many applications are stated for other particular physically relevant lattices including the triangular, square, cubic, orthorhombic, body-centered-cubic, and face-centered-cubic lattices.


Introduction
The mathematical justification of a crystal's shape and formation is a very difficult problem that 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], twodimensional Thomas-Fermi model for a 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 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-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 et al. [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 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 centered at the origin and all the masses μ L centered 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 Fig. 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 centered at z and all the masses μ L centered at lattice sites (see the bottom line of Fig. 1).
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 that 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 has 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 with these soft theta functions. Furthermore, we show that the minimum of z → θ H+z (α), where H is a honeycomb lattice (see (2.6) and Fig. 4), is the center of a 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 Sect. 2, we show some preliminary results in Sect. 3, including the generalized Jacobi transformation formula. Our main results are stated and proved in Sect. 4, and many applications are given in Sect. 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 centered at the origin. Furthermore, we write M d (R) for the space of n × n matrices with real coefficients.
Furthermore, the dual lattice of L ∈ L • d is defined by We also recall the definitions of the following important lattices belonging to L • d (see also Fig. 2

):
The triangular lattice 1 (2.5) 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 Fig. 4). This is a typical example of periodic configuration, i.e., a 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] 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 centered at L.
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 are 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.
as r → +∞ for some σ > 0 and if f can be represented as the Laplace transform of a Radon measure μ f , i.e., 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, 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 notation 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 the 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 μ centered 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 (2.8) In Fig. 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 for 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 another measure ν z located at z ∈ R d , and is therefore a generalization of the lattice and translated lattice theta functions θ L (α) and θ L+z (α) defined by (1.1).
, 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
is a union of two translated triangular lattices, it is easy to compute its Gaussian energy per point. More precisely, we have
• 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 state precisely 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: 1. We say that L 0 ∈ D d is a critical point or a (strict local) minimum of L → θ μ L +ν 0 (α) at all scales on D d if, for any α > 0, L 0 is a critical point or a (strict local) minimum on D d of L → θ μ L +ν 0 (α). 2. Let α 0 > 0 be fixed. We say that L 0 ∈ D d is a critical point or a (strict local) minimum of L → θ μ L +ν 0 (α 0 ) at a low scale on D d if there exist ε 0 > 0 and δ 0 > 0 such that, for any ε ∈ [0, ε 0 ) and any δ ∈ [0, δ 0 ), L 0 is a critical point or a (strict local) minimum of L → θ μ ε Furthermore, let L be fixed, and 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 the 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 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 thatm 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.
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 Jd 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.

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   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 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. 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 L 0 for where, for any measure m ∈ P r (R d ), First, because L 0 is a critical point of L → θ L (α) in D d for all α > 0, it implies, by Proposition 4.1, that L 0 is a critical point of E h ε,δ in D d . Second, we also know that   ε j 1 (s, t, |q|) + δ j 2 (s, t, |q|) + j ε,δ (s, t, |q|) , 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 centered 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 . 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 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 that is also a parameter of L, once the lattice is parametrized by n d parameters as explained in Sect. 2, and ∈ L • d−1 . We therefore get 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 denote by K . In the following, we write, as in (4.2), and we recall that minimizing L → θ μ ε because, by Lemma 3.4, L 0 = L * 0 . As in [12], we claim there exists h ε,δ such that, for any L ∈ K , for some constant C > 0 independent of ε, δ. If (4.3)-(4.4) hold, then, for any L ∈ K , Thus, for any L ∈ K such that E h ε,δ [L] ≤ E h ε,δ [L 0 ], this implies that d(L, L 0 ) ≤ C max{ε, δ}, which contradicts the strict local minimality of L 0 for sufficiently small ε and δ. To constructh ε,δ , the idea is to approximate by a bounded completely monotone function k ε,δ = L[m ε,δ ] such that m ε,δ is a positive measure with a compact support, k ε,δ (0) = 1 and k ε,δ ≤ C. It is indeed sufficient to apply the method described in [12, proof of Thm. 2] to k 1 ε and k 2 δ , which are then approximated respectively by k 1 ε = L[m ε ] and k 2 ε = L[m δ ] and where m ε,δ = m ε * m δ . We therefore get, for any r > 0, 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 that h ε,δ = 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 ingredients 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 the Poisson summation formula (see, e.g., [45,Appendix A]), Since h is completely monotone,ĥ 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 to 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. 1. In dimension d = 2, if L = 1 , then Baernstein [3] proved that the barycenters z 1 =  [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:

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.
Moreover, for all α 0 > 0, 1 is a global minimizer of L → θ μ L +ν 0 (α 0 ) in D  be the barycenters of the primitive triangles of Q 1 . Then, for any μ, ν ∈ P r (R 2 ), z 1 and z 2 are critical points of z → θ μ 1 +ν z (α) in Q 1 at all scales. Furthermore, for any α 0 > 0, z 1 and z 2 are the unique global minimizers of z → θ μ 1 +ν z (α 0 ) in Q 1 at a low scale. Moreover, if μ, ν ∈ P cm r (R 2 ), then the minimality of z 1 and z 2 holds at all scales.
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 an 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. as depicted in Fig. 4. 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 (α) .

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 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  1. For any μ, ν ∈ P r (R d ), Z d is a critical point of L → θ μ L +ν 0 (α) in D d at all scales. 2. For any μ, ν ∈ P r (R d ) and any α 0 > 0, a = (1, 1, ..., 1) is the unique minimizer of a → θ μ Z d a +ν 0 (α 0 ) in {(a 1 , ..., a d ) ∈ (0, +∞) : d i=1 a i = 1} at a low scale; i.e., Z d is the unique minimizer of L → θ μ L +ν 0 (α 0 ) among orthorhombic lattices at a low scale. Moreover, if μ, ν ∈ P cm r (R d ), the minimality holds at all scales. 3. Let (a 1 , ..., a d ) ∈ (0, +∞) be such that d i=1 a i = 1, and let c a := 1 2 (a 1 , ..., a d ) be the center of Q Z d a . Then, for any μ, ν ∈ P r (R d ), c a is a critical point of z → θ μ Z d a +ν z (α) in Q Z d a at all scales. Furthermore, for any α 0 > 0, c a is the unique global minimizer of z → θ μ Z d a +ν z (α 0 ) in Q Z d a at a low scale. Moreover, if μ, ν ∈ P cm r (R d ), the minimality holds at all scales.
Remark 5.6 (Octahedral site of the cubic lattice) In the cubic case (a 1 , ..., a d ) = (1, ..., 1), the center c = (1/2, ..., 1/2), also called the octahedral site of Z d , is indeed the preferred location to add an extra atom in order to create, for example, an ion (like the CsCl which has a BCC structure).
In [13], it has been proved, using a computer assisted method, that D * 3 (resp. D 3 ) is a strict local minimizer of L → θ L (α) for any α ∈ A (resp. α ∈ B) defined by