Wigner bosonic molecules with repulsive interactions and harmonic confinement

In earlier work, the momentum distribution of harmonically confined bosons has been calculated for interparticle interactions which are also harmonic for an arbitrary number N of repelling Bose particles. After a brief discussion of the pair function in this model, attention is also given, but now for N=2, to the single-particle density for the repulsive interaction potential λ/r122, where r12 is the interboson distance and λ the repulsive coupling strength. The single-particle density is displayed for specific values of λ, with attention focusing on the limit of strong repulsive coupling. Throughout, the focus is on the fingerprints of Wigner bosonic molecules.