**A. Python code used to calculated**

```
from ase.build import bulk
from gpaw import GPAW, PW, FermiDirac