| 1 |
Madelung Energy Sample Calculation |
| 2 |
|
| 3 |
The sample in this directory provides a way of checking the value of |
| 4 |
the Madelung Energy for a perfect crystal of NaCl. The relevant |
| 5 |
quantities are: |
| 6 |
|
| 7 |
M (Madelung constant) = 1.74756 |
| 8 |
a (lattice constant) = 5.65 Angstroms |
| 9 |
q^2 / (4 pi e0 a) = 58.77233 kcal / mol |
| 10 |
M q^2 / (4 pi e0 a) = 102.708173 kcal / mol |
| 11 |
|
| 12 |
The file NaCl.md contains 8000 ions, so the total electrostatic energy |
| 13 |
of the perfect crystal in this file should be: |
| 14 |
|
| 15 |
V_electrostatic = -821665.38 kcal / mol |
| 16 |
|
| 17 |
Using different electrostatic calculation methods, we can get quite |
| 18 |
close to this value. |
| 19 |
|
| 20 |
For example, with : |
| 21 |
cutoffMethod = "shifted_force"; |
| 22 |
electrostaticScreeningMethod = "damped"; |
| 23 |
cutoffRadius = 28; |
| 24 |
dampingAlpha = 0.14159292; |
| 25 |
|
| 26 |
The resultant electrostatic potential is: -821667.68 kcal / mol |
| 27 |
|
| 28 |
To obtain values for the electrostatic potential in OpenMD, we add the |
| 29 |
ELECTROSTATIC_POTENTIAL keyword to the end of the statFileFormat: |
| 30 |
|
| 31 |
statFileFormat = "TIME|TOTAL_ENERGY|POTENTIAL_ENERGY|KINETIC_ENERGY|TEMPERATURE|PRESSURE|VOLUME|CONSERVED_QUANTITY|ELECTROSTATIC_POTENTIAL"; |
