On minima of difference of theta functions and application to hexagonal crystallization
Let $z=x+iy \in \mathbb{H}:=\{z= x+ i y\in\mathbb{C}: y>0\}$ and $ θ(α;z)=\sum_{(m,n)\in\mathbb{Z}^2 } e^{-α\frac{π}{y }|mz+n|^2}$ be the theta function associated with the lattice $L ={\mathbb Z}\oplus z{\mathbb Z}$. In this paper we consider the following minimization problem of difference of two theta functions \begin{equation}\aligned\nonumber \min_{ \mathbb{H} } \Big(θ(α; z)-βθ(2α; z)\Big) \endaligned\end{equation} where $α\geq 1$ and $ β\in (-\infty, +\infty)$. We prove that there is a critical value $β_c=\sqrt2$ (independent of $α$) such that if $β\leqβ_c$, the minimizer is $\frac{1}{2}+i\frac{\sqrt3}{2}$ (up to translation and rotation) which corresponds to the hexagonal lattice, and if $β>β_c$, the minimizer does not exist. Our result partially answers some questions raised in \cite{Bet2016, Bet2018, Bet2020, Bet2019AMP} and gives a new proof in the crystallization of hexagonal lattice under Yukawa potential.