One-electron densities of freely rotating Wigner molecules

A formalism enabling computation of the one-particle density of a freely rotating assembly of identical particles that vibrate about their equilibrium positions with amplitudes much smaller than their average distances is presented. It produces densities as finite sums of products of angular and radial functions, the length of the expansion being determined by the interplay between the point-group and permutational symmetries of the system in question. Obtaining from a convolution of the rotational and bosonic components of the parent wavefunction, the angular functions are state-dependent. On the other hand, the radial functions are Gaussians with maxima located at the equilibrium lengths of the position vectors of individual particles and exponents depending on the scalar products of these vectors and the eigenvectors of the corresponding Hessian as well as the respective eigenvalues. Although the new formalism is particularly useful for studies of the Wigner molecules formed by electrons subject to weak confining potentials, it is readily adaptable to species (such as ʻballiums' and Coulomb crystals) composed of identical particles with arbitrary spin statistics and permutational symmetry. Several examples of applications of the present approach to the harmonium atoms within the strong-correlation regime are given.