Crystalline boson phases in harmonic traps: Beyond the Gross-Pitaevskii mean field

Strongly-interacting bosons in two-dimensional harmonic traps are described through breaking of rotational symmetry at the Hartree-Fock level and subsequent symmetry restoration via projection techniques, thus incorporating correlations beyond the Gross-Pitaevskii (GP) solution. The bosons localize and form polygonal-ringlike crystalline patterns, both for a repulsive contact potential and a Coulomb interaction, as revealed via conditional-probability-distribution analysis. For neutral bosons, the total energy of the crystalline phase saturates in contrast to the GP solution, and its spatial extent becomes smaller than that of the GP condensate. For charged bosons, the total energy and dimensions approach the values of classical pointlike charges in their equilibrium configuration.