First-principles evolutionary simulations are used to systematically predict stable compounds in the Zr-O system at pressures up to 120 GPa. Zr-O compounds and Hf-O compounds share many similarities, but four new phases Cmmm-Zr 3 O, R 3c-Zr 3 O 2 , Pmma-ZrO, and Fe 2 P-type ZrO 2 (P 62m) appear in the Zr-O system: the latter two phases appear in the Hf-O system at higher pressure. OII ZrO 2 (Pnma) transforms into Fe 2 P-type ZrO 2 at 102 GPa. Meanwhile, Fe 2 P-type ZrO 2 and P 62m-ZrO have similar structures based on x-Zr. However, the calculated Vickers hardness of Fe 2 P-type ZrO 2 (5.6 GPa) is inferior to that of P 62m-ZrO (14.1 GPa). The hardness of P 62m-ZrO (14.1 GPa) is lower than that of P 62m-HfO (16.1 GPa) and P 62m-TiO (16.6 GPa). On the whole, Zr-O compounds exhibit lower hardnesses and bulk moduli than Hf-O compounds. Published by AIP Publishing. [http://dx.