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.