Madelung constant

Last updated December 18, 2023

The Madelung constant is used in determining the electrostatic potential of a single ion in a crystal by approximating the ions by point charges. It is named after Erwin Madelung, a German physicist.^[1]

Because the anions and cations in an ionic solid attract each other by virtue of their opposing charges, separating the ions requires a certain amount of energy. This energy must be given to the system in order to break the anion–cation bonds. The energy required to break these bonds for one mole of an ionic solid under standard conditions is the lattice energy.

Formal expression

The Madelung constant allows for the calculation of the electric potential $V i$ of the ion at position $r i$ due to all other ions of the lattice

V_{i}={\frac {e}{4\pi \varepsilon _{0}}}\sum _{j\neq i}{\frac {z_{j}}{r_{ij}}}\,\!

where $r_{ij}=|r_{i}-r_{j}|$ is the distance between the $i$ th and the $j$ th ion. In addition,

$z j =$ number of charges of the $j$ th ion
$e =$ the elementary charge, 1.6022×10⁻¹⁹ C
$4 πε 0 =$ 1.112×10⁻¹⁰ C²/(J⋅m); $ε 0$ is the permittivity of free space.

If the distances $r ij$ are normalized to the nearest neighbor distance $r 0$ , the potential may be written

V_{i}={\frac {e}{4\pi \varepsilon _{0}r_{0}}}\sum _{j}{\frac {z_{j}r_{0}}{r_{ij}}}={\frac {e}{4\pi \varepsilon _{0}r_{0}}}M_{i}

with $M i$ being the (dimensionless) Madelung constant of the $i$ th ion

M_{i}=\sum _{j}{\frac {z_{j}}{r_{ij}/r_{0}}}.

Another convention is to base the reference length on the cubic root $w$ of the unit cell volume, which for cubic systems is equal to the lattice constant. Thus, the Madelung constant then reads

{\overline {M}}_{i}=\sum _{j}{\frac {z_{j}}{r_{ij}/w}}=M_{i}{\frac {w}{r_{0}}}.

The electrostatic energy of the ion at site $r i$ then is the product of its charge with the potential acting at its site

E_{el,i}=z_{i}eV_{i}={\frac {e^{2}}{4\pi \varepsilon _{0}r_{0}}}z_{i}M_{i}.

There occur as many Madelung constants $M i$ in a crystal structure as ions occupy different lattice sites. For example, for the ionic crystal NaCl, there arise two Madelung constants – one for Na and another for Cl. Since both ions, however, occupy lattice sites of the same symmetry they both are of the same magnitude and differ only by sign. The electrical charge of the Na⁺ and Cl⁻ ion are assumed to be onefold positive and negative, respectively, $z Na = 1$ and $z Cl = -1$ . The nearest neighbour distance amounts to half the lattice constant of the cubic unit cell $r_{0}={\tfrac {a}{2}}$ and the Madelung constants become

M_{\text{Na}}=-M_{\text{Cl}}={\sum _{j,k,\ell =-\infty }^{\infty }}\!\!\!\!^{\prime }\ \ {{(-1)^{j+k+\ell }} \over {\sqrt {j^{2}+k^{2}+\ell ^{2}}}}.

This graph demonstrates the non-convergence of the expanding spheres method for calculating the Madelung constant for NaCl as compared to the expanding cubes method, which is convergent. NaClMadelungConstant.png — This graph demonstrates the non-convergence of the expanding spheres method for calculating the Madelung constant for NaCl as compared to the expanding cubes method, which is convergent.

The prime indicates that the term $j=k=\ell =0$ is to be left out. Since this sum is conditionally convergent it is not suitable as definition of Madelung's constant unless the order of summation is also specified. There are two "obvious" methods of summing this series, by expanding cubes or expanding spheres. Although the latter is often found in the literature,^[2]

M=-6+{\frac {12}{\sqrt {2}}}-{\frac {8}{\sqrt {3}}}+{\frac {6}{2}}-{\frac {24}{\sqrt {5}}}+\dots

it fails to converge, as was shown by Emersleben in 1951.^[3] The summation over expanding cubes converges to the correct value, although very slowly. An alternative summation procedure, presented by Borwein, Borwein and Taylor, uses analytic continuation of an absolutely convergent series.^[4]

There are many practical methods for calculating Madelung's constant using either direct summation (for example, the Evjen method^[5]) or integral transforms, which are used in the Ewald method.^[6] A fast converging formula for the Madelung constant of NaCl is

12\,\pi \sum _{m,n\geq 1,\,\mathrm {odd} }\operatorname {sech} ^{2}\left({\frac {\pi }{2}}(m^{2}+n^{2})^{1/2}\right)

^[7]

**Examples of Madelung constants**
Ion in crystalline compound	$M$ (based on $r 0$ )	${\overline {M}}$ (based on $w$ )
Cl⁻ and Cs⁺ in CsCl	±1.762675	±2.035362
Cl⁻ and Na⁺ in rocksalt NaCl	±1.747565	±3.495129
S²⁻ and Zn²⁺ in sphalerite ZnS	±3.276110	±7.56585
F⁻ in fluorite CaF₂	1.762675	4.070723
Ca²⁺ in fluorite CaF₂	-3.276110	−7.56585

The continuous reduction of $M$ with decreasing coordination number $Z$ for the three cubic AB compounds (when accounting for the doubled charges in ZnS) explains the observed propensity of alkali halides to crystallize in the structure with highest $Z$ compatible with their ionic radii. Note also how the fluorite structure being intermediate between the caesium chloride and sphalerite structures is reflected in the Madelung constants.

Generalization

It is assumed for the calculation of Madelung constants that an ion's charge density may be approximated by a point charge. This is allowed, if the electron distribution of the ion is spherically symmetric. In particular cases, however, when the ions reside on lattice site of certain crystallographic point groups, the inclusion of higher order moments, i.e. multipole moments of the charge density might be required. It is shown by electrostatics that the interaction between two point charges only accounts for the first term of a general Taylor series describing the interaction between two charge distributions of arbitrary shape. Accordingly, the Madelung constant only represents the monopole-monopole term.

The electrostatic interaction model of ions in solids has thus been extended to a point multipole concept that also includes higher multipole moments like dipoles, quadrupoles etc.^[8]^[9]^[10] These concepts require the determination of higher order Madelung constants or so-called electrostatic lattice constants. The proper calculation of electrostatic lattice constants has to consider the crystallographic point groups of ionic lattice sites; for instance, dipole moments may only arise on polar lattice sites, i. e. exhibiting a C₁, C_1h, C_n or C_nv site symmetry (n = 2, 3, 4 or 6).^[11] These second order Madelung constants turned out to have significant effects on the lattice energy and other physical properties of heteropolar crystals.^[12]

Application to organic salts

The Madelung constant is also a useful quantity in describing the lattice energy of organic salts. Izgorodina and coworkers have described a generalised method (called the EUGEN method) of calculating the Madelung constant for any crystal structure.^[13]

Related Research Articles

Electrostatics is a branch of physics that studies slow-moving or stationary electric charges.

<span class="mw-page-title-main">Polaron</span> Quasiparticle in condensed matter physics

A polaron is a quasiparticle used in condensed matter physics to understand the interactions between electrons and atoms in a solid material. The polaron concept was proposed by Lev Landau in 1933 and Solomon Pekar in 1946 to describe an electron moving in a dielectric crystal where the atoms displace from their equilibrium positions to effectively screen the charge of an electron, known as a phonon cloud. This lowers the electron mobility and increases the electron's effective mass.

The DLVO theory explains the aggregation and kinetic stability of aqueous dispersions quantitatively and describes the force between charged surfaces interacting through a liquid medium. It combines the effects of the van der Waals attraction and the electrostatic repulsion due to the so-called double layer of counterions. The electrostatic part of the DLVO interaction is computed in the mean field approximation in the limit of low surface potentials - that is when the potential energy of an elementary charge on the surface is much smaller than the thermal energy scale, $. For two spheres of radius each having a charge separated by a center-to-center distance in a fluid of dielectric constant containing a concentration of monovalent ions, the electrostatic potential takes the form of a screened-Coulomb or Yukawa potential,$

A multipole expansion is a mathematical series representing a function that depends on angles—usually the two angles used in the spherical coordinate system for three-dimensional Euclidean space, $. Similarly to Taylor series, multipole expansions are useful because oftentimes only the first few terms are needed to provide a good approximation of the original function. The function being expanded may be real- or complex-valued and is defined either on, or less often on for some other .$

Jellium, also known as the uniform electron gas (UEG) or homogeneous electron gas (HEG), is a quantum mechanical model of interacting electrons in a solid where the positive charges are assumed to be uniformly distributed in space; the electron density is a uniform quantity as well in space. This model allows one to focus on the effects in solids that occur due to the quantum nature of electrons and their mutual repulsive interactions without explicit introduction of the atomic lattice and structure making up a real material. Jellium is often used in solid-state physics as a simple model of delocalized electrons in a metal, where it can qualitatively reproduce features of real metals such as screening, plasmons, Wigner crystallization and Friedel oscillations.

The chemists Peter Debye and Erich Hückel noticed that solutions that contain ionic solutes do not behave ideally even at very low concentrations. So, while the concentration of the solutes is fundamental to the calculation of the dynamics of a solution, they theorized that an extra factor that they termed gamma is necessary to the calculation of the activities of the solution. Hence they developed the Debye–Hückel equation and Debye–Hückel limiting law. The activity is only proportional to the concentration and is altered by a factor known as the activity coefficient $. This factor takes into account the interaction energy of ions in solution.$

Surface states are electronic states found at the surface of materials. They are formed due to the sharp transition from solid material that ends with a surface and are found only at the atom layers closest to the surface. The termination of a material with a surface leads to a change of the electronic band structure from the bulk material to the vacuum. In the weakened potential at the surface, new electronic states can be formed, so called surface states.

In chemistry, the lattice energy is the energy change upon formation of one mole of a crystalline ionic compound from its constituent ions, which are assumed to initially be in the gaseous state. It is a measure of the cohesive forces that bind ionic solids. The size of the lattice energy is connected to many other physical properties including solubility, hardness, and volatility. Since it generally cannot be measured directly, the lattice energy is usually deduced from experimental data via the Born–Haber cycle.

Ewald summation, named after Paul Peter Ewald, is a method for computing long-range interactions in periodic systems. It was first developed as the method for calculating the electrostatic energies of ionic crystals, and is now commonly used for calculating long-range interactions in computational chemistry. Ewald summation is a special case of the Poisson summation formula, replacing the summation of interaction energies in real space with an equivalent summation in Fourier space. In this method, the long-range interaction is divided into two parts: a short-range contribution, and a long-range contribution which does not have a singularity. The short-range contribution is calculated in real space, whereas the long-range contribution is calculated using a Fourier transform. The advantage of this method is the rapid convergence of the energy compared with that of a direct summation. This means that the method has high accuracy and reasonable speed when computing long-range interactions, and it is thus the de facto standard method for calculating long-range interactions in periodic systems. The method requires charge neutrality of the molecular system to accurately calculate the total Coulombic interaction. A study of the truncation errors introduced in the energy and force calculations of disordered point-charge systems is provided by Kolafa and Perram.

In condensed matter physics and crystallography, the static structure factor is a mathematical description of how a material scatters incident radiation. The structure factor is a critical tool in the interpretation of scattering patterns obtained in X-ray, electron and neutron diffraction experiments.

In condensed-matter physics, channelling (or channeling) is the process that constrains the path of a charged particle in a crystalline solid.

In statistical mechanics, the two-dimensional square lattice Ising model is a simple lattice model of interacting magnetic spins. The model is notable for having nontrivial interactions, yet having an analytical solution. The model was solved by Lars Onsager for the special case that the external magnetic field H = 0. An analytical solution for the general case for $has yet to be found.$

The Poisson–Boltzmann equation is a useful equation in many settings, whether it be to understand physiological interfaces, polymer science, electron interactions in a semiconductor, or more. It aims to describe the distribution of the electric potential in solution in the direction normal to a charged surface. This distribution is important to determine how the electrostatic interactions will affect the molecules in solution. The Poisson–Boltzmann equation is derived via mean-field assumptions. From the Poisson–Boltzmann equation many other equations have been derived with a number of different assumptions.

The Debye–Hückel theory was proposed by Peter Debye and Erich Hückel as a theoretical explanation for departures from ideality in solutions of electrolytes and plasmas. It is a linearized Poisson–Boltzmann model, which assumes an extremely simplified model of electrolyte solution but nevertheless gave accurate predictions of mean activity coefficients for ions in dilute solution. The Debye–Hückel equation provides a starting point for modern treatments of non-ideality of electrolyte solutions.

The Kapustinskii equation calculates the lattice energy U_L for an ionic crystal, which is experimentally difficult to determine. It is named after Anatoli Fedorovich Kapustinskii who published the formula in 1956.

The Born–Landé equation is a means of calculating the lattice energy of a crystalline ionic compound. In 1918 Max Born and Alfred Landé proposed that the lattice energy could be derived from the electrostatic potential of the ionic lattice and a repulsive potential energy term.

The isotopic shift is the shift in various forms of spectroscopy that occurs when one nuclear isotope is replaced by another.

The electric dipole moment is a measure of the separation of positive and negative electrical charges within a system, that is, a measure of the system's overall polarity. The SI unit for electric dipole moment is the coulomb-meter (C⋅m). The debye (D) is another unit of measurement used in atomic physics and chemistry.

Double layer forces occur between charged objects across liquids, typically water. This force acts over distances that are comparable to the Debye length, which is on the order of one to a few tenths of nanometers. The strength of these forces increases with the magnitude of the surface charge density. For two similarly charged objects, this force is repulsive and decays exponentially at larger distances, see figure. For unequally charged objects and eventually at shorted distances, these forces may also be attractive. The theory due to Derjaguin, Landau, Verwey, and Overbeek (DLVO) combines such double layer forces together with Van der Waals forces in order to estimate the actual interaction potential between colloidal particles.

The Born–Mayer equation is an equation that is used to calculate the lattice energy of a crystalline ionic compound. It is a refinement of the Born–Landé equation by using an improved repulsion term.

References

↑ Madelung E (1918). "Das elektrische Feld in Systemen von regelmäßig angeordneten Punktladungen". Phys. Z. XIX: 524–533.
↑ Charles Kittel: Introduction to Solid State Physics , Wiley 1995, ISBN 0-471-11181-3
↑ Emersleben, O. (1951). "Das Selbstpotential einer endlichen Reihe neutraler äquidistanter Punktepaare". Mathematische Nachrichten . 4 (3–4): 468. doi:10.1002/mana.3210040140.
↑ Borwein, D.; Borwein, J. M.; Taylor, K. F. (1985). "Convergence of Lattice Sums and Madelung's Constant". J. Math. Phys. 26 (11): 2999–3009. Bibcode:1985JMP....26.2999B. doi:10.1063/1.526675. hdl: 1959.13/1043576 .
↑ Evjen, H. M. (1932). "On the Stability of Certain Heteropolar Crystals" (PDF). Phys. Rev. 39 (4): 675–687. Bibcode:1932PhRv...39..675E. doi:10.1103/physrev.39.675.
↑ Ewald, P. P. (1921). "Die Berechnung optischer und elektrostatischer Gitterpotentiale". Ann. Phys. 64 (3): 253–287. Bibcode:1921AnP...369..253E. doi:10.1002/andp.19213690304.
↑ Bailey, David; Borwein, Jonathan; Kapoor, Vishaal; Weisstein, Eric (March 9, 2006). "Ten Problems in Experimental Mathematics" (PDF). The American Mathematical Monthly. 113 (6): 481. doi:10.2307/27641975. JSTOR 27641975.
↑ J. Kanamori; T. Moriya; K. Motizuki & T. Nagamiya (1955). "Methods of Calculating the Crystalline Electric Field". J. Phys. Soc. Jpn. 10 (2): 93–102. Bibcode:1955JPSJ...10...93K. doi:10.1143/JPSJ.10.93.
↑ B. R. A. Nijboer & F. W. de Wette (1957). "On the calculation of lattice sums". Physica. 23 (1–5): 309–321. Bibcode:1957Phy....23..309N. doi:10.1016/S0031-8914(57)92124-9. hdl: 1874/15643 . S2CID 122383484.
↑ E. F. Bertaut (1978). "The equivalent charge concept and its application to the electrostatic energy of charges and multipoles". J. Phys. (Paris). 39 (2): 1331–48. Bibcode:1978JPCS...39...97B. doi:10.1016/0022-3697(78)90206-8.
↑ M. Birkholz (1995). "Crystal-field induced dipoles in heteropolar crystals – I. concept". Z. Phys. B. 96 (3): 325–332. Bibcode:1995ZPhyB..96..325B. CiteSeerX 10.1.1.424.5632 . doi:10.1007/BF01313054. S2CID 122527743.
↑ M. Birkholz (1995). "Crystal-field induced dipoles in heteropolar crystals – II. physical significance". Z. Phys. B. 96 (3): 333–340. Bibcode:1995ZPhyB..96..333B. doi:10.1007/BF01313055. S2CID 122393358.
↑ E. Izgorodina; et al. (2009). "The Madelung Constant of Organic Salts". Crystal Growth & Design. 9 (11): 4834–4839. doi:10.1021/cg900656z.

External links

Glasser, Leslie (2012). "Solid-state energetics and electrostatics: Madelung constants and Madelung energies". Inorg. Chem. 51 (4): 2420–2424. doi:10.1021/ic2023852. PMID 22242970.
Sakamoto, Y. (1958). "Madelung constants of simple crystals expressed in terms of Born's basic potentials of 15 figures". J. Chem. Phys. 28 (1): 164–165. Bibcode:1958JChPh..28..164S. doi:10.1063/1.1744060.
Sakamoto, Y. (1958). "Errata 2: Madelung constants of simple crystals expressed in terms of Born's basic potentials of 15 figures". J. Chem. Phys. 28 (6): 1253. Bibcode:1958JChPh..28.1253S. doi: 10.1063/1.1744387 .
Zucker, I. J. (1975). "Madelung constants and lattice sums for invariant cubic lattice complexes and certain tetragonal structures". J. Phys. A: Math. Gen. 8 (11): 1734–1745. Bibcode:1975JPhA....8.1734Z. doi:10.1088/0305-4470/8/11/008.
Zucker, I. J. (1976). "Functional equations for poly-dimensional zeta functions and the evaluation of Madelung constants". J. Phys. A: Math. Gen. 9 (4): 499–505. Bibcode:1976JPhA....9..499Z. doi:10.1088/0305-4470/9/4/006.
Weisstein, Eric W. "Madelung Constants". MathWorld .
OEIS sequenceA085469(Decimal expansion of Madelung constant (negated) for NaCl structure)

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[1] Madelung E (1918). "Das elektrische Feld in Systemen von regelmäßig angeordneten Punktladungen". Phys. Z. XIX: 524–533.

[2] Charles Kittel: Introduction to Solid State Physics , Wiley 1995, ISBN 0-471-11181-3

[3] Emersleben, O. (1951). "Das Selbstpotential einer endlichen Reihe neutraler äquidistanter Punktepaare". Mathematische Nachrichten . 4 (3–4): 468. doi:10.1002/mana.3210040140.

[4] Borwein, D.; Borwein, J. M.; Taylor, K. F. (1985). "Convergence of Lattice Sums and Madelung's Constant". J. Math. Phys. 26 (11): 2999–3009. Bibcode:1985JMP....26.2999B. doi:10.1063/1.526675. hdl: 1959.13/1043576 .

[5] Evjen, H. M. (1932). "On the Stability of Certain Heteropolar Crystals" (PDF). Phys. Rev. 39 (4): 675–687. Bibcode:1932PhRv...39..675E. doi:10.1103/physrev.39.675.

[6] Ewald, P. P. (1921). "Die Berechnung optischer und elektrostatischer Gitterpotentiale". Ann. Phys. 64 (3): 253–287. Bibcode:1921AnP...369..253E. doi:10.1002/andp.19213690304.

[7] Bailey, David; Borwein, Jonathan; Kapoor, Vishaal; Weisstein, Eric (March 9, 2006). "Ten Problems in Experimental Mathematics" (PDF). The American Mathematical Monthly. 113 (6): 481. doi:10.2307/27641975. JSTOR 27641975.

[Kana1955-8] J. Kanamori; T. Moriya; K. Motizuki & T. Nagamiya (1955). "Methods of Calculating the Crystalline Electric Field". J. Phys. Soc. Jpn. 10 (2): 93–102. Bibcode:1955JPSJ...10...93K. doi:10.1143/JPSJ.10.93.

[Nijb1957-9] B. R. A. Nijboer & F. W. de Wette (1957). "On the calculation of lattice sums". Physica. 23 (1–5): 309–321. Bibcode:1957Phy....23..309N. doi:10.1016/S0031-8914(57)92124-9. hdl: 1874/15643 . S2CID 122383484.

[Bert1978-10] E. F. Bertaut (1978). "The equivalent charge concept and its application to the electrostatic energy of charges and multipoles". J. Phys. (Paris). 39 (2): 1331–48. Bibcode:1978JPCS...39...97B. doi:10.1016/0022-3697(78)90206-8.

[ZPB1995a-11] M. Birkholz (1995). "Crystal-field induced dipoles in heteropolar crystals – I. concept". Z. Phys. B. 96 (3): 325–332. Bibcode:1995ZPhyB..96..325B. CiteSeerX 10.1.1.424.5632 . doi:10.1007/BF01313054. S2CID 122527743.

[ZPB1995b-12] M. Birkholz (1995). "Crystal-field induced dipoles in heteropolar crystals – II. physical significance". Z. Phys. B. 96 (3): 333–340. Bibcode:1995ZPhyB..96..333B. doi:10.1007/BF01313055. S2CID 122393358.

[Izgorodina2009-13] E. Izgorodina; et al. (2009). "The Madelung Constant of Organic Salts". Crystal Growth & Design. 9 (11): 4834–4839. doi:10.1021/cg900656z.

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]

[11]

[12]

[13]