“…We carried out MD simulations using LAMMPS [25] with the extended Finnis-Sinclair potential [26] which gives the correct elastic constants at zero pressure and reproduces the pressure dependence of the specific volume up to values close to 500 GPa. The extended Finnis-Sinclair potential was tested by verifying the gamma surface [27] and comparing it with ab initio results.…”