Efficient sampling of ice structures by electrostatic switching

An electrostatic switching procedure is introduced that enables the systematic generation of high-quality ice configurations at various temperatures. Proton disordered ice Ih configurations were generated for the TIP4P water model at temperatures from 50 to 240 K, for the SPC/Fw water model from 100 to 240 K, and for the DC97 water model at 240 K. The resulting configurations were found to properly sample the canonical ensemble. The dielectric constant of ice Ih was determined from the net dipole fluctuation of the ice configurations. The calculated dielectric constant compares favorably with the study by Rick and Haymet [J. Chem. Phys. 2003, 118, 9291]. However, our method gives smaller error bars, especially at lower temperatures. At temperatures above 200 K, a type of hydrogen-bond defect is identified that cannot be categorized as a D or L type defect.