Melting of Ice under Pressure

The melting of ice under pressure is investigated with a series of first principles molecular dynamics simulations. In particular, a two-phase approach is used to determine the melting temperature of the ice-VII phase in the range of 10 to 50 GPa. Our computed melting temperatures are consistent with existing diamond anvil cell experiments. We find that for pressures between 10 to 40 GPa, ice melts as a molecular solid. For pressures above ∼45 GPa there is a sharp increase in the slope of the melting curve due to the presence of molecular dissociation and proton diffusion in the solid, prior to melting. The onset of significant proton diffusion in ice-VII as a function of increasing temperature is found to be gradual and bears many similarities to that of a type-II superionic solid