We report fundamental progress in the description of an archetypical hydrogen-bonded system: H2O-ice. We calculated by first principles the ice VII-VIII phase boundary, an order-disorder phase boundary, by performing a complete sampling of configurations in a 16-molecule supercell. Vibrational effects are included within quasiharmonic theory. Several important features were successfully described: negative Clapeyron slope, isotope effect, and nearly-constant Tc at low pressures. The equation of state was very well described below 40 GPa. Inclusion of tunneling effects and increase of supercell size should decrease the remaining discrepancies, i.e., the overestimation of volume above ∼40 GPa and the underestimation of Tc.
Bibliographical noteFunding Information:
KU and RMW thank SISSA for support during the early stage of this research. Calculations have been carried out using the Quantum-ESPRESSO distribution  . This research was supported by NSF grants EAR-0757903 , EAR-0635990 , EAR-1047629 , and ATM-0428774 (VLab). Computations were performed at the Minnesota Supercomputing Institute.