This is a question triggered by this post
Madelung's constant is defined to the coefficient of electrostatic potential energy in a ionic crystal. In the example of $NaCl$, \begin{equation} M = \sum_{ijk}{}^{'}\frac{(-1)^{i+j+k}}{\sqrt{i^2+j^2+k^2}} \end{equation} is conditionally convergent.
Since this sum is conditionally convergent it is not suitable as definition of Madelung's constant unless the order of summation is also specified. There are two "obvious" methods of summing this series, by expanding cubes or expanding spheres. The latter, though devoid of a meaningful physical interpretation (there are no spherical crystals) is rather popular because of its simplicity. Thus, the following expansion is often found in the literature:[2] $$ M = -6 +12/ \sqrt{2} -8/ \sqrt{3} +6/2 - 24/ \sqrt{5} + \dotsb = -1.74756\dots. $$ However, this is wrong as this series diverges as was shown by Emersleben in 1951.[3][4] The summation over expanding cubes converges to the correct value. An unambiguous mathematical definition is given by Borwein, Borwein and Taylor by means of analytic continuation of an absolutely convergent series.
I have the following questions.
1) Expanding sphere leads to a divergent series. OK, what if I make a perfect spherical sample of $NaCl$, will the experimentally measured Madelung's constant to be infinity?
(EDIT: I maybe have some misunderstanding of divergence: it could be the case that the series is bounded but don't have a definite limit. So is the series for expanding sphere bounded? and if there is no definite limit, what's the experimentally measured value for a perfect spherical crystal?)
2) What physical principle dictates the order of summation? or why does finite number obtained by analytic continuation should be consistent with the observed value?
3) What is the role of charge neutrality here? I ever programed to compute Madelung's constant using the fractional charge idea((assign $\frac{1}{8}$ charge to the corner, $\frac{1}{4}$ to the edge, $\frac{1}{2}$ to the face ), which makes the expanding cube charge neutral. Does that mean charge neutrality is one of conditions that must be enforced? or it is just for the sake of computational efficiency?