High pressure ices

<jats:p> H <jats:sub>2</jats:sub> O will be more resistant to metallization than previously thought. From computational evolutionary structure searches, we find a sequence of new stable and meta-stable structures for the ground state of ice in the 1–5 TPa (10 to 50 Mbar) regime, in the static approximation. The previously proposed <jats:italic>Pbcm</jats:italic> structure is superseded by a <jats:italic>Pmc</jats:italic> 2 <jats:sub>1</jats:sub> phase at <jats:italic>p</jats:italic>  = 930 GPa, followed by a predicted transition to a <jats:italic>P</jats:italic> 2 <jats:sub>1</jats:sub> crystal structure at <jats:italic>p</jats:italic>  = 1.3 TPa. This phase, featuring higher coordination at O and H, is stable over a wide pressure range, reaching 4.8 TPa. We analyze carefully the geometrical changes in the calculated structures, especially the buckling at the H in O-H-O motifs. All structures are insulating—chemistry burns a deep and (with pressure increase) lasting hole in the density of states near the highest occupied electronic levels of what might be component metallic lattices. Metallization of ice in our calculations occurs only near 4.8 TPa, where the metallic <jats:italic>C</jats:italic> 2/ <jats:italic>m</jats:italic> phase becomes most stable. In this regime, zero-point energies much larger than typical enthalpy differences suggest possible melting of the H sublattice, or even the entire crystal. </jats:p>