Standard density functionals without van der Waals interactions yield an unsatisfactory description of ice phases, specifically, high density phases occurring under pressure are too unstable compared to the common low density phase I h observed at ambient conditions. Although the description is improved by using functionals that include van der Waals interactions, the errors in relative volumes remain sizable. Here we assess the random phase approximation (RPA) for the correlation energy and compare our results to experimental data as well as diffusion Monte Carlo data for ice. The RPA yields a very balanced description for all considered phases, approaching the accuracy of diffusion Monte Carlo in relative energies and volumes. This opens a route towards a concise description of molecular water phases on surfaces and in cavities.

Several proton-disordered crystalline ice structures are known to proton order at sufficiently low temperatures, provided that the right preparation procedure is used. For cubic ice, ice Ic, however, no proton ordering has been observed so far. Here, we subject ice Ic to an experimental protocol similar to that used to proton order hexagonal ice. In situ FT-IR spectroscopy carried out during this procedure reveals that the librational band of the spectrum narrows and acquires a structure that is observed neither in proton-disordered ice Ic nor in ice XI, the proton-ordered variant of hexagonal ice. On the basis of vibrational spectra computed for ice Ic and four of its proton-ordered variants using classical molecular dynamics and ab initio simulations, we conclude that the features of our experimental spectra are due to partial proton ordering, providing the first evidence of proton ordering in cubic ice. We further find that the proton-ordered structure with the lowest energy is ferroelectric, while the structure with the second lowest energy is weakly ferroelectric. Both structures fit the experimental spectral similarly well such that no unique assignment of proton order is possible based on our results.