Classification of particle numbers with unique Heitmann-Radin minimizer

TL;DR

This paper characterizes when minimizers of the Heitmann-Radin energy are unique, identifying a specific sequence of particle numbers for which this holds, using discrete differential geometry techniques.

Contribution

It provides a complete classification of particle numbers with unique Heitmann-Radin minimizers, extending understanding of energy minimization in discrete systems.

Findings

Unique minimizers occur for a specific infinite sequence of particle numbers.

The sequence of particle numbers with unique minimizers is explicitly characterized.

The proof employs advanced discrete differential geometry methods.

Abstract

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).

Equations

$${\int }_{\partial \Omega }\phi (\nu (x))dA(x).$$

$${\mathcal{E}}_{HR}(X):=\frac{1}{2}\sum _{1\le i<j\le N}{V}_{HR}(|x_i-x_j|),$$

$$V_{HR}(r)=\left\{\begin{array}[]{ll}+\infty&\textrm{if }r<1\\ -1&\textrm{if }r=1\\ 0&\textrm{if }r>1.\end{array}\right.$$

$$\mathcal{L}=\{i\mathbf{e}+j\mathbf{f}\,:\,i,j\in\mathbb{Z}\},\quad\mathbf{e}=\left(\begin{array}[]{l}1\\ 0\end{array}\right),\quad\mathbf{f}=\left(\begin{array}[]{l}1/2\\ {\sqrt{3}}/2\end{array}\right).$$

$$(a)N=3s^2+3s+1$$

$$(b)N=3s^2+3s+1+(s+1)k+s$$

$${\mathcal{E}}_{HR}(X)=-3N+P(X)+\mu (X)+3\chi (X),$$

$$\mu (X)=♯\text{quadrilaterals}+2♯\text{pentagons}+3♯\text{hexagons}+\dots ,$$

$$\begin{gathered} {\mathcal{X}}_{\mathcal{L}}^N:=\{X\subset \mathcal{L}:♯X=N,\text{mbox}allfacesofX\text{mbox}aretriangles, \\ X\text{ has simply closed polygonal boundary}\}. \end{gathered}$$

$$\underset{X\subset {\mathbb{R}}^2,♯X=N}{\min }{\mathcal{E}}_{HR}(X)=\underset{X\in {\mathcal{X}}_{\mathcal{L}}^N}{\min }{\mathcal{E}}_{HR}(X)=-3N+3+\underset{X\in {\mathcal{X}}_{\mathcal{L}}^N}{\min }P(X).$$

$$N=3s^2(N)+3s(N)+1+(s(N)+1)k(N)+j(N).$$

$$\min_{X\in\mathcal{X}_{\mathcal{L}}^{N}}P(X)=\left\{\begin{array}[]{ll}6\,s(N)&\textrm{if }k(N)=j(N)=0,\\ 6\,s(N)+k(N)+1&\textrm{otherwise}.\end{array}\right.$$

$$P(X)\ge P(X^{\prime })+6.$$

$$X:=X^{\prime }\cup \{x\in \mathcal{L}:\text{there exists }{x}^{\prime }\in X^{\prime }\text{ with }|x-x^{\prime }|=1\};$$

$$P(X)=P(X^{\prime })+6.$$

$$6s+k+1=P(X)\ge P(X^{\prime })+6\ge \underset{Y\in {\mathcal{X}}_{\mathcal{L}}^{N^{\prime }}}{\min }P(Y)+6=6(s-1)+k+1+6,$$

$$\{N\ge 3:\text{mbox}theminimizerofP\text{mbox}in{\mathcal{X}}_{\mathcal{L}}^N\text{mbox}isunique\}=\{a_n{\}}_{n\ge 3},$$

$$a_{n}:=\left\{\begin{array}[]{ll}3\,s^{2}+3\,s+1&\textrm{if }n=6\,s,\\ 3\,s^{2}+3\,s+1+(s+1)k+s&\textrm{if }n=6\,s+k+1,\,k\in\{0,1,2,3,4\}.\end{array}\right.$$

$$a_n=\max \{N\in \mathbb{N}:\underset{X\in {\mathcal{X}}_{\mathcal{L}}^N}{\min }P(X)=n\}.$$

$$x+2y+3z=n.$$

Full text

1. 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 $\Omega$ occupied by the atoms. Such a surface energy governing the shape $\Omega$ was first written down by Gibbs, and in modern notation has the form

[TABLE]

Here $\Omega$ is an – up to regularity requirements arbitrary – subset of $\mathbb{R}^{d}$, $\nu(x)$ is the outward unit normal to $\Omega$ at the point $x$, $\varphi$ 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 ${\mathcal{H}}^{d-1}$) on $\partial\Omega$. 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\in\mathbb{R}^{d}\,:\,x\cdot\nu(x)\leq\varphi(x)\mbox{ for all }x\in\partial\Omega\}$. 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 $\mathbb{R}^{2}$.

The Heitmann-Radin energy for a system of $N$ identical particles with positions $x_{1},..,x_{N}\in\mathbb{R}^{2}$, introduced in [6], is

[TABLE]

where here and below we abbreviate the position vector of all particles by $X\in\mathbb{R}^{2N}$ and the interaction potential $V_{HR}$ is given by

[TABLE]

The potential (1.3) arises naturally from the Lennard-Jones potential $V(r)=r^{-2p}-2r^{-p}$ by passing to the limit $p\to\infty$ [3]. In [6], Heitmann and Radin proved the fundamental result that for any fixed $N\in\mathbb{N}$, the configurations $X$ minimizing $\mathcal{E}_{HR}$ are, up to rotation and translation, subsets of the triangular lattice

[TABLE]

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 $\varphi$ is derived for the model (1.2), and it is shown that as the number of $N$ of particles gets large, the empirical measure $\mu_{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.

Theorem 1.1**.**

Let $N\in\mathbb{N}$. The minimizers of $\mathcal{E}_{HR}$ among $N$-particle configurations are unique up to translation and rotation if and only if either

[TABLE]

for some $s\in\mathbb{N}\cup\{0\}$, or

[TABLE]

for some $s\in\mathbb{N}\cup\{0\}$ and some $k\in\{0,1,2,3,4\}$.

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 $\mathbb{R}^{2}$ which are nonempty, do not contain any point in $X$, and whose boundary is given by a cycle $\cup_{i=1}^{k}[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:

[TABLE]

where $\chi(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 $\mu(X)$ is the defect measure

[TABLE]

which can be viewed as a distance measure between the bond graph and vacancy-free subsets of the triangular lattice $\mathcal{L}$. See Figure 2.

2. 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.

As in the Figure, let $k\in\{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\leq 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 $\mu=0$ and $\chi=1$, and does not change the perimeter, and hence, by (1.4), yields another minimizer. When $k\geq 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\subset\mathbb{R}^{2}$ with $\sharp X=N$, endowed with the planar graph structure described in the Introduction, we say that $x\in 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 $\partial 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\subset\mathbb{R}^{2}$ with $\sharp X=N$, belong – up to translation and rotation and for $N\geq 3$ – to the set of crystallized configurations

[TABLE]

This reduces the uniqueness question to uniqueness among configurations in $\mathcal{X}_{\mathcal{L}}^{N}$, and together with (1.4) shows that

[TABLE]

Moreover, to describe the numerical value of the ground state energy, we recall that for any $N\in\mathbb{N}$ there exists a uniquely determined triple

$(s(N),k(N),j(N))\in(\mathbb{N}\cup\{0\})^{3}$ with $k(N)\leq 5$ and $j(N)\leq s(N)$ such that

[TABLE]

In terms of the numbers $s$, $k$, and $j$, it follows from [6, 3] that

[TABLE]

Finally, we will need the following result which is an immediate consequence of [3, Lemma 4.1, Lemma 4.2].

Lemma 2.1**.**

Let $N,N^{\prime}\in\mathbb{N}$.

Let $X\in\mathcal{X}_{\mathcal{L}}^{N}$ and let $X^{\prime}:=X\setminus\partial X$. Then

[TABLE]

Conversely, let $X^{\prime}\in\mathcal{X}_{\mathcal{L}}^{N^{\prime}}$ and set

[TABLE]

then

[TABLE]

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 $\mathcal{E}_{HR}$ among $N$-particle configurations; then by the Heitmann-Radin theorem $X$ belongs (up to rotation and translation) to $\mathcal{X}_{\mathcal{L}}^{N}$, and by (2.1) it satisfies $P(X)=6\,s+k+1$. Set $X^{\prime}:=X\setminus\partial X$. We have that $\sharp X^{\prime}=N_{s}-6\,s-k-1=N_{s-1}$. By Lemma 2.1, we conclude that

[TABLE]

and hence all the above inequalities are equalities. Therefore $X^{\prime}$ is a minimizer. By the inductive assumption, the minimizer $X^{\prime}$ is unique (up to rotation and translation). Since, by the construction in Lemma 2.1, the set $X$ is fully determined by $X^{\prime}$, we obtain the claim.

3. Concluding remarks

The above arguments show that

[TABLE]

where

[TABLE]

Moreover, using (2.1) one can easily check that for any $n\geq 3$, $a_{n}$ has a simple geometric meaning: it is the maximum particle number whose Heitmann-Radin minimizer has perimeter $n$, in formulae:

[TABLE]

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

[TABLE]

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.

Acknowledgements: We are indebted to Michael Baake for pointing out to us the website [10] after attending a lecture on the work reported here. The reesearch of LDL was funded, and that of GF partially supported, by the DFG Collaborative Research Center CRC 109 “Discretization in Geometry and Dynamics”.