The Emergence of the Hexagonal Lattice in Two-Dimensional Wigner Fragments

At very low density, the electrons in a uniform electron gas spontaneously break symmetry and form a crystalline lattice called a Wigner crystal. But which type of crystal will the electrons form? We report a numerical study of the density profiles of fragments of Wigner crystals from first principles. To simulate Wigner fragments, we use Clifford periodic boundary conditions and a renormalized distance in the Coulomb potential. Moreover, we show that high-spin restricted open-shell Hartree-Fock theory becomes exact in the low-density limit. We are thus able to accurately capture the localization in two-dimensional Wigner fragments with many electrons. No assumptions about the positions where the electrons will localize are made. The density profiles we obtain emerge naturally when we minimize the total energy of the system. We clearly observe the emergence of the hexagonal crystal structure, which has been predicted to be the ground-state structure of the two-dimensional Wigner crystal.