To allow direct comparison to our ab initio results, we first computed K with the finite strain difference method. The reduced computational cost allowed for smoother volumetric deformations (1.5% change over 250 ps) and a longer equilibration (250 ps, shown on Figure ) than was possible with AIMD. This approach led to well-behaved P -V plots and to the values presented in Table .