The calculation of the equation of state of each phase was done as follows. First, we determined the equilibrium structure at zero pressure and at 50 GPa. Convergence criteria of 10 –4 Ry/bohr in the ...