The evaluation of the mechanical properties was carried out by calculating the elastic properties by performing finite distortions of the lattice. The elastic constants C ij were then derived from the strainstress relationship, assuming statistically isotropic solid on a large scale. The C ij were then used to obtain the bulk (B) and shear (G) moduli using the Voigt averaging method .