Centroidal Voronoi Tessellations as Electrostatic Equilibria: A Generalized Thomson Problem in Convex Domains

Abstract: We present a variational framework in which Centroidal Voronoi Tessellations (CVTs) arise as local minimizers of a generalized electrostatic energy functional. By modeling interior point distributions in a convex domain as repelling charges balanced against a continuous boundary charge, we show that the resulting equilibrium configurations converge to CVT structures. We prove this by showing that CVTs minimize both the classical centroidal energy and the electrostatic potential, establishing a connection between geometric quantization and potential theory. Finally, we introduce a thermodynamic annealing scheme for global CVT optimization, rooted in Boltzmann statistics and random walk dynamics. By introducing a scheme for varying time steps (faster or slower cooling) we show that the set of minima of the centroid energy functional (and therefore the electrostatic potential) can be recovered. By recovering a set of generator locations corresponding to each minimum we can create a lattice continuation that allows for a customizable framework for individual minimum seeking.