Dipole field issue in particle-mesh Ewald method with periodic boundary conditions I am working on a thesis that makes a great use of molecular-dynamics simulations, and I am trying to understand how the particle-mesh Ewald method  works. The problem is, I have difficulties understanding its very premise; now I will explain what I think I have learnt:
Long-range electrostatic forces do not converge if periodic boundary conditions are enforced; thus we cannot obtain them by summing pairwise interactions between every charge in the system if we take the periodic images of charges into account. The problem is, I see this as an issue with serious practical consequences, and I cannot imagine how a mathematical rearranging could solve it. 
I’ll make an example:
Let’s suppose we have an electrically neutral unit cell that gains an electric dipole moment during the simulation. This unit cell, and its dipole moment, would be instantly reproduced in all image cells.  Let’s consider the electric potential in an arbitrary point: it would have an infinite number of shells of dipoles around, with area increasing as $r^2$ ($r$ is the shell radius) and the electric field dipole component would decrease as $r^{-2}$; the sum does not converge, so we would have an infinite electric potential in every point! Am I missing something?
I cannot see how PME or Ewald Summation, or any other algorithm, can solve a physical issue, unless that in some way those methods put additional boundary conditions. But I don’t see how. Can you help me understand? Thank you in advance.
EDIT: I was wrong about the infinite potential, because there is a cosine term in the dipole component that zeroes the potential in my proposed shell-by-shell calculation. Anyway, if we consider the electric field instead, we have it falling at $r^{-3}$; by adding the field produced by each shell we obtain a series of $1/r$ terms, that is still diverging at infinity, so my problem is still unsolved.
 A: Well, the starting point is the fact that $1/r$ Coulomb interaction which diverges at infinity in k-space is only $1/k^2$ which converges at infinity, but diverges at 0. If one takes only parts of both which do converge, then it is possible to avoid divergences at all. 
Do you agree with this and the problem appears at the next step or you do not agree with this starting point? 
In fact, Ewald summation does not make magically divergent integrals converge. If you try to take the periodical system with "+" and "-" charges, the divergency at infinity stemming from Coulomb divergency is unphysical: at large distance, the average charge seen by distant probe charge is zero. However, if you break it in parts, the energy from interaction with "+" and "-" charges separately does diverge. Ewald summation is the rigorous mathematical trick to avoid calculation of interaction with "+"s and "-"s separately. One could as well compute the interaction with dipoles which would not diverge. 
In your question, you are stressing the role of periodical boundary conditions. Are they really important? Have you tried to consider the infinite lattice instead? 
A: I solved the issue. I was actually missing something. I didn't notice that the electric field at the center of a spherical shell of ideal dipoles is zero. I think it is not obvious if you look at the formula, but it can be demonstrated by using the continuous limit and integrating on the entire shell. Since we work with neutral simulation boxes (of course), and the other multipole components of the electric field are short-range, I think we can say that the total field converges. So the divergence is a mathematical issue indeed, not a physical one.
