perturbation theory is, like in the name, an approximation of an unknown system by the assumption that it is a little different from a known system.
In our case, is multipole expansion the known system to approximate the charge distribution.
iirc, we dont have a “correct” model of the nucleus which means we cannot use that as reference.
Multipole expansion is an unbiased way to “model” it.