Next Order Asymptotics and Renormalized Energy for Riesz Interactions

We study systems of $n$ points in the Euclidean space of dimension $d \ge 1$ interacting via a Riesz kernel $|x|^{-s}$ and confined by an external potential, in the regime where $d-2\le s