Lowest-known Energy Configuration of $N=100\,000$ Coulomb Charges in a Disk: Breaking the $10^{5}$ Barrier

Abstract

We report the first calculation of the lowest-known energy configuration for $N=100\,000$ classical point charges confined to a disk and interacting via the $1/r$ Coulomb potential - a system size never before achieved for the Thomson problem in a disk, more than doubling the previous record of $N=40\,886$ reported by Amore and Zarate. An adaptive defect-targeting subdomain optimization strategy yields an approximately linear growth of the wall-clock time per cycle with $N$; for $N=100\,000$ a single $24$-core CPU workstation without GPU acceleration required $\approx31$ hours. The energy $E_{\min}=7.80466624157\times10^{9}$ deviates from the standard empirical asymptotic expansion by only $2.525\times10^{-7}$. We demonstrate that this remarkable agreement provides a stringent test of the rigorous next-order asymptotic theory for Riesz interactions, while the residual deviation is consistent with an $O(N^{4/3})$ boundary-layer correction whose absence from the standard fitting basis biases the empirical $N^{3/2}$ coefficient away from its exact crystalline value, thereby tying the deviation to the two-dimensional crystallization conjecture. Bond-orientational order analysis reveals a polycrystalline bulk threaded by radial grain boundaries, establishing a new benchmark for two-dimensional Coulomb systems.