/EwaldSummation

Ewald summation program for computing the long range Coulomb interactions in 3D Periodic systems

Primary LanguageC++

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 alt tag

2. TiO2 (Rutile)

Madelung constant for TiO2 (Rutile) = 2.40 alt tag

3. ZnS (Sphalerite)

Madelung constant for ZnS (Sphalerite) = 1.63 alt tag

4. Al2O3 (corrundum)

Madelung constant for Al2O3 (corrundum) = 4.17 alt tag