Editing Ewald sum
Jump to navigation
Jump to search
The edit can be undone. Please check the comparison below to verify that this is what you want to do, and then publish the changes below to finish undoing the edit.
Latest revision | Your text | ||
Line 1: | Line 1: | ||
{{stub-general}} | |||
The '''Ewald sum''' technique <ref>[http://dx.doi.org/10.1002/andp.19213690304 Paul Ewald "Die Berechnung Optischer und Electrostatischer Gitterpotentiale", Annalen der Physik '''64''' pp. 253-287 (1921)]</ref> was originally developed by Paul Ewald to evaluate the Madelung constant <ref>[http://dx.doi.org/10.1063/1.1727895 S. G. Brush, H. L. Sahlin and E. Teller "Monte Carlo Study of a One-Component Plasma. I", Journal of Chemical Physics '''45''' pp. 2102-2118 (1966)]</ref>. It is now widely used in order to simulate systems with | The '''Ewald sum''' technique <ref>[http://dx.doi.org/10.1002/andp.19213690304 Paul Ewald "Die Berechnung Optischer und Electrostatischer Gitterpotentiale", Annalen der Physik '''64''' pp. 253-287 (1921)]</ref> was originally developed by Paul Ewald to evaluate the Madelung constant <ref>[http://dx.doi.org/10.1063/1.1727895 S. G. Brush, H. L. Sahlin and E. Teller "Monte Carlo Study of a One-Component Plasma. I", Journal of Chemical Physics '''45''' pp. 2102-2118 (1966)]</ref>. It is now widely used in order to simulate systems with | ||
[[long range interactions]] (typically, [[Electrostatics |electrostatic interactions]]). Its aim is the computation of the interaction of a system with [[periodic boundary conditions]] with all its replicas. This is accomplished by the introduction of fictitious "charge clouds" that shield the charges. The interaction is then divided into a shielded part, which may be evaluated by the usual means, and a part that cancels the introduction of the clouds, which is evaluated in [[Fourier_analysis | Fourier space]]. | [[long range interactions]] (typically, [[Electrostatics |electrostatic interactions]]). Its aim is the computation of the interaction of a system with [[periodic boundary conditions]] with all its replicas. This is accomplished by the introduction of fictitious "charge clouds" that shield the charges. The interaction is then divided into a shielded part, which may be evaluated by the usual means, and a part that cancels the introduction of the clouds, which is evaluated in [[Fourier_analysis | Fourier space]]. | ||
Line 4: | Line 5: | ||
In a periodic system one wishes to evaluate the [[internal energy]] <math>U</math> (Eq. 1.1 <ref>[http://dx.doi.org/10.1098/rspa.1980.0135 S. W. de Leeuw, J. W. Perram and E. R. Smith "Simulation of Electrostatic Systems in Periodic Boundary Conditions. I. Lattice Sums and Dielectric Constants", Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences '''373''' pp. 27-56 (1980)]</ref>): | In a periodic system one wishes to evaluate the [[internal energy]] <math>U</math> (Eq. 1.1 <ref>[http://dx.doi.org/10.1098/rspa.1980.0135 S. W. de Leeuw, J. W. Perram and E. R. Smith "Simulation of Electrostatic Systems in Periodic Boundary Conditions. I. Lattice Sums and Dielectric Constants", Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences '''373''' pp. 27-56 (1980)]</ref>): | ||
:<math>U = \frac{1}{2} {\sum_{\mathbf n}}^ | :<math>U = \frac{1}{2} {\sum_{\mathbf n}}^' \left[ \sum_{i=1}^N \sum_{j=1}^N \phi \left({\mathbf r}_{ij} + L{\mathbf n}, {\mathbf \Omega_i}, {\mathbf \Omega_j} \right) \right] </math> | ||
where one sums over all the [[Building up a simple cubic lattice | simple cubic lattice]] points <math>{\mathbf n} = (l,m,n)</math>. The prime on the first summation indicates that if <math>i=j</math> then the <math>{\mathbf n} = 0</math> term is omitted. <math>L</math> is the length of the side of the cubic simulation box, <math>N</math> is the number of particles, and <math>{\mathbf \Omega}</math> represent the [[Euler angles]]. | where one sums over all the [[Building up a simple cubic lattice | simple cubic lattice]] points <math>{\mathbf n} = (l,m,n)</math>. The prime on the first summation indicates that if <math>i=j</math> then the <math>{\mathbf n} = 0</math> term is omitted. <math>L</math> is the length of the side of the cubic simulation box, <math>N</math> is the number of particles, and <math>{\mathbf \Omega}</math> represent the [[Euler angles]]. |