We consider a system of N bosons in the mean-field scaling regime for a class of interactions including the repulsive Coulomb potential. We derive an asymptotic expansion of the low-energy eigenstates and the corresponding energies, which provides corrections to Bogoliubov theory to any order in
$1/N$.