Harmonically confined particles with long-range repulsive interactions

We study an interacting system of $N$ classical particles on a line at thermal equilibrium. The particles are confined by a harmonic trap and repelling each other via pairwise interaction potential that behaves as a power law $\propto \sum_{\substack{i\neq j}}^N|x_i-x_j|^{-k}$ (with $k>-2$) of their mutual distance. This is a generalization of the well known cases of the one component plasma ($k=-1$), Dyson's log-gas ($k\to 0^+$), and the Calogero-Moser model ($k=2$). Due to the competition between harmonic confinement and pairwise repulsion, the particles spread over a finite region of space for all $k>-2$. We compute exactly the average density profile for large $N$ for all $k>-2$ and show that while it is independent of temperature for sufficiently low temperature, it has a rich and nontrivial dependence on $k$ with distinct behavior for $-2<k<1$, $k>1$ and $k=1$.