{'fraction_non_alphanumeric': 0.057854878478337445, 'fraction_numerical': 0.02333568157802043, 'mean_word_length': 3.981140350877193, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 14, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Classification of particle numbers with unique Heitmann-Radin minimizer

We show that minimizers of the Heitmann-Radin energy (R. C. Heitmann, C. Radin, J. Stat. Phys. 22, 281-287, 1980) are unique if and only if the particle number N belongs to an infinite sequence whose first thirty-five elements are 1, 2, 3, 4, 5, 7, 8, 10, 12, 14, 16, 19, 21, 24, 27, 30, 33, 37, 40, 44, 48, 52, 56, 61, 65, 70, 75, 80, 85, 91, 96, 102, 108, 114, 120 (see the paper for a closed-form description of this sequence). The proof relies on the discrete differential geometry techniques introduced in (L. De Luca, G. Friesecke, preprint).


Introduction
A fundamental problem in statistical and solid mechanics is to explain theoretically why atoms at low temperature self-assemble into subsets of periodic lattices, and why these subsets exhibit specific polyhedral shapes.
In previous studies, the specific shapes have been beautifully explained under two simplifications. First, one assumes crystallization, i.e. one restricts the admissible atomic positions to lattice sites. Second, one passes to a coarse-grained description in which atomistic energy minimization is replaced by minimization of an effective surface energy of the region Ω occupied by the atoms. Such a surface energy governing the shape Ω was first written down by Gibbs, and in modern notation has the form (1.1) ∂Ω ϕ(ν(x)) dA(x).
Here Ω is an -up to regularity requirements arbitrary -subset of R d , ν(x) is the outward unit normal to Ω at the point x, ϕ is a surface energy density which captures the fact that interfaces with certain orientations with respect to the crystal lattice are favoured over others, and dA(x) is the usual area element (Hausdorff measure H d−1 ) on ∂Ω. A fundamental result going back to Taylor [9] and Fonseca and Müller [5] states that minimizers of (1.1) among sets of finite perimeter and fixed volume are unique up to translation, and given by a dilation of the Wulff shape {x ∈ R d : x · ν(x) ≤ ϕ(x) for all x ∈ ∂Ω}. For a related uniqueness result at finite temperature in the context of the 2D Ising model see [4].
Microscopically, the situation regarding uniqueness is much more subtle. Our goal in this note is to settle the uniqueness question completely in case of the two-dimensional Heitmann-Radin model, the perhaps simplest model describing the self-assembly of atoms into crystalline order and special shapes despite allowing arbitrary particle positions in R 2 . The Heitmann-Radin energy for a system of N identical particles with positions x 1 , .., x N ∈ R 2 , introduced in [6], is where here and below we abbreviate the position vector of all particles by X ∈ R 2N and the interaction potential V HR is given by The potential (1.3) arises naturally from the Lennard-Jones potential V (r) = r −2p − 2r −p by passing to the limit p → ∞ [3]. In [6], Heitmann and Radin proved the fundamental result that for any fixed N ∈ N, the configurations X minimizing E HR are, up to rotation and translation, subsets of the triangular lattice Moreover the authors give an explicit formula for the ground state energy from which it is elementary to deduce examples of non-uniqueness. Such examples by no means contradicts the macroscopic uniqueness result for (1.1); indeed in [1] the associated macroscopic surface energy density ϕ is derived for the model (1.2), and it is shown that as the number of N of particles gets large, the empirical measure µ N associated with any sequence of microscopic minimizers converges after re-scaling to the characteristic function of the unique Wulff shape. Further insight into the uniqueness question was recently achieved by Schmidt [7]. He showed that microscopic regular hexagons ( Figure 1, left picture) are the unique microscopic minimizer for their particle number, by re-writing the discrete energy in the form (1.1) via introducing suitable tiles around each particle and applying the macroscopic uniqueness theorem for (1.1). Moreover in [7] it is shown that for particle numbers N exceeding those of a microscopic regular hexagon by 1, the amount of nonuniqueness can be surprisingly large (for a refinement of this result see [2]).
Here we completely settle the question of when microscopic uniqueness occurs. The result is as follows.  An explicit list of the numbers between 1 and 120 satisfying (a) or (b) was given in the Abstract. Case (a) corresponds to regular hexagons (see the leftmost picture in Figure 1), recovering the uniqueness result of [7]. The other uniqueness cases were neither previously conjectured, nor can they be established with the same method, as the microscopic minimizers are no longer regular hexagons and hence the associated continuum sets introduced in [7] do not minimize the correspoding continuum energy (1.1). Our proof of Theorem 1.1 relies on the discrete differential geometry approach introduced recently by us in [3]. In [3] we used this approach to give a new proof of the Heitmann-Radin crystallization theorem; here we employ it to settle the uniqueness question. This approach starts by associating, to each particle configuration X = (x 1 , .., x N ), its bond graph, a planar graph (see Figure 2) whose definition we recall here: vertices correspond to the particle positions x j , edges to line segments [x j , x k ] connecting two particle positions of distance 1, and faces to open bounded subsets of R 2 which are nonempty, do not contain any point in X, and whose boundary is given by a cycle ∪ k i=1 [x i−1 , x i ] of edges for some points x 0 , .., x k with x k = x 0 . A key result of [3] needed in the proof of Theorem 1.1 is the following geometric decomposition of the Heitmann-Radin energy on N -particle configurations: where χ(X) is the Euler characteristic of the bond graph of X, P (X) is its combinatorial perimeter as introduced in [3], namely the number of boundary edges with "wire edges" counted twice (see Figure 2), and µ(X) is the defect measure (1.5) µ(X) = quadrilaterals + 2 pentagons + 3 hexagons + . . . , which can be viewed as a distance measure between the bond graph and vacancy-free subsets of the triangular lattice L. See Figure 2. Particles are said to be connected by an edge (or bond) if their distance is 1. For the bond graph above, the Euler characteristic χ is 1 (since the graph is connected), the defect measure µ (see (1.5)) equals 3 (since the graph contains a square and a pentagon), and the combinatorial perimeter equals 17. The latter is because there are 13 regular boundary edges, i.e. edges lying on the boundary of precisely one face, and 2 wire edges, i.e. edges not lying on the boundary of any face, which must be counted twice. It follows that the right hand side of (1.4) equals −3 · 17 + 17 + 3 + 3 · 1 = −28, which indeed agrees with the Heitmann-Radin energy of the configuration.

Proof of the theorem
We begin by showing non-uniqueness in the case when the particle number N is not of the form (a) or (b) in Theorem 1.1. First notice that for such N , the "canonical" minimizer from [6, 3] (see Figure 3, left picture) has a non-convex angle at the boundary. . Two minimizers of E HR for N = 106 = 3 s 2 + 3 s + 1 + (s + 1)k + j, with s = 5, k = 2, and j = 3. Left: The "canonical" minimizer from [6,3]: starting from the grey particles constituting a regular hexagon with sidelength s, the remaining (red) particles are added by following a counterclockwise path around the hexagon starting from the encircled point. Right: Another minimizer, obtained from the one on the left by moving the boxed segment of atoms.
As in the Figure, let k ∈ {0, 1, 2, 3, 4, 5} be the number of complete sides of the canonical minimizer that were added to the central regular hexagon. When k ≤ 2, we can move the first not yet covered side of the central hexagon (boxed particles, left picture in Figure 3) on top of the last not yet covered side of the central hexagon (see Figure 3, right picture). This preserves µ = 0 and χ = 1, and does not change the perimeter, and hence, by (1.4), yields another minimizer. When k ≥ 3, we can instead move the first new side of the canonical minimizer obtained by covering the central hexagon on top of the third covered side, again leaving defect measure, Euler characteristic and perimeter -and hence the Heitmann-Radin energy -unchanged. Either way, the new minimizer is not a rotated translate of the canonical one.
Before moving to the proof of uniqueness when N is of form (a) or (b), we introduce some notation (which was also used in [3]). For an N -point configuration, i.e. X ⊂ R 2 with X = N , endowed with the planar graph structure described in the Introduction, we say that x ∈ X is a boundary particle if it is adjacent to a boundary edge, i.e. an edge lying on the boundary of at most one face. We denote the set of boundary particles by ∂X, and say that X has simply closed polygonal boundary if the union of boundary edges forms a simply closed curve.
Let us also recall the Heitmann-Radin crystallization theorem proved in [6] (see also [3]). This result says that minimizers of the Heitmann-Radin energy (1.2) among arbitrary N -particle configurations, i.e. X ⊂ R 2 with X = N , belong -up to translation and rotation and for N ≥ 3 -to the set of crystallized configurations Let X ∈ X N L and let X := X \ ∂X. Then P (X) ≥ P (X ) + 6.
We are now in a position to prove the uniqueness statement in Theorem 1.1. We establish the claim only for N of the form (b), the proof in case (a) being analogous (and an alternative proof in case (a) being already known [7]). We use induction on s = s(N ). It is easy to see that the claim is satisfied for s = 0, 1, for any value of k in {0, 1, 2, 3, 4}. Fix such a k. Assuming that the minimizer is unique for N s−1 := 3(s − 1) 2 + 3(s − 1) + 1 + s k + s − 1, we prove that it is unique for N s := 3 s 2 + 3 s + 1 + (s + 1)k + s. Let X be a minimizer of E HR among N -particle configurations; then by the Heitmann-Radin theorem X belongs (up to rotation and translation) to X N L , and by (2.1) it satisfies P (X) = 6 s + k + 1. Set X := X \ ∂X. We have that X = N s − 6 s − k − 1 = N s−1 . By Lemma 2.1, we conclude that 6 s + k + 1 = P (X) ≥ P (X ) + 6 ≥ min Y ∈X N L P (Y ) + 6 = 6 (s − 1) + k + 1 + 6, and hence all the above inequalities are equalities. Therefore X is a minimizer. By the inductive assumption, the minimizer X is unique (up to rotation and translation). Since, by the construction in Lemma 2.1, the set X is fully determined by X , we obtain the claim.
Moreover, using (2.1) one can easily check that for any n ≥ 3, a n has a simple geometric meaning: it is the maximum particle number whose Heitmann-Radin minimizer has perimeter n, in formulae: (3.2) a n = max{N ∈ N : min X∈X N L P (X) = n}.
It appears that the sequence {a n } is well-known in number theory (see for instance [8,10]). For example, a n seems to co-incide with the number of nonnegative integer solutions to the Diophantine equation x + 2 y + 3 z = n.
Nevertheless, to the best of our knowledge, neither the geometric characterizations of a n in Theorem 1.1 and (3.2) nor the explicit one given in (3.1) have been observed before. It would be very interesting both from a physical and a mathematical point of view if a direct argument linking the Heitmann-Radin model with the above Diophantine equation could be found.