Dear Juerg:<div><br></div><div>That is very helpful.  Thank you.  I will now go back to my old calculations (where I got the Berry Phase dipole, used periodic boundary conditions, and calculated values for the dielectric constant that were off by many orders of magnitude) and double-check my units.  However, I checked them previously so think I probably have a different problem.  <br><br>This is the fluctuation formula I have been using:<br>https://www.cp2k.org/exercises:2015_ethz_mmm:monte_carlo_ice<br></div><div><br></div><div>If anyone knows of any reason it would be inappropriate/inconsistent with the input file I posted, it would be great to know.  Thank you again for your help; this has been troubling me for a long time.</div>