EwaldSummation
Ewald summation program for computing the long range Coulomb interactions in 3D Periodic systems.
Because the equation in Ewald summation program is complicated, it must be checked and verified. Here we do this by calculating the Madelung constant for the crystal structure of SiO2(beta-cristobalite) , TiO2(rutile), Al2O3 (corrundum) and ZnS(Sphalerite). If the Madeulng constant is correct, then the Ewald sum works correctly. The calculated Madelung constants for the crystal structures are shown in the below figures:
1. Silica (Beta-cristobalite)
Madelung constant for silica (beta-cristobalite) = 2.22
2. TiO2 (Rutile)
Madelung constant for TiO2 (Rutile) = 2.40
3. ZnS (Sphalerite)
Madelung constant for ZnS (Sphalerite) = 1.63
4. Al2O3 (corrundum)
Madelung constant for Al2O3 (corrundum) = 4.17
