An exact analytical solution in strongly non-local media with a harmonic potential has been studied. Two-dimensional Bessel solitary wave clusters have been obtained by a self-similar method. The intensity distributions of optical beam with different parameters have been discussed in detail. It is found that the solitary waves have a symmetric necklace distribution and the number of facular points is double the value of the quantum number $m$. The modulation of the external potential field to the width of light beam is also discussed.
