Madelungkonstante

Die Madelung-Konstante α ist ein einheitenloser Faktor, der als das Verhältnis der durchschnittlichen Bindungsenergie pro Ion im Kristallgitter zur durchschnittlichen Bindungsenergie pro Ion bei einem einzelnen Ionenpaar definiert ist. Er hängt dabei nur vom Strukturtyp ab und ist unabhängig von der Ionenladung und den Zellparametern. Typische Vertreter dieser Kristallgruppe sind die Alkalimetallhalogenide (Beispiele sind NaCl, KBr und CsCl), bei denen die Bindung durch Coulombkräfte entsteht. Dabei gibt das Metallatom ein Elektron an das Halogenatom ab und es entsteht eine kugelsymmetrische Ladungsverteilung.

Die Konstante wurde nach Erwin Madelung, der sie zuerst ableitete, benannt^[1].

Weil die Madelung-Konstante vom Coulomb-Gesetz für Punktladungen abgeleitet ist, verliert sie ihre Gültigkeit bei nicht-punktförmigen Ionen (Ionen mit kovalenten Bindungen wie z. B. im Pyritkristall) und bei Ionen mit unterschiedlicher Polarität (z. B. in der Reihe ZnS, TiO₂, CdCl₂, CdI₂).

Berechnung der Bindungsenergie im Gitter

Die Bindungsenergie für ein Ionenpaar lässt sich mittels Coulomb-Gesetz wie folgt berechnen:

$E_{IP}=\int_\infty^d F\,\mathrm{d}r =\int_\infty^d\frac{1}{4\pi\varepsilon_0}\cdot \frac{z^2e^2}{r^2}\,\mathrm{d}r =-\frac{1}{4\pi\varepsilon_0}\cdot \frac{z^2e^2}{d}$

z – Ladung der Ionen

e – Elementarladung

d – (kleinster) Abstand der Ionen

ε₀ – Dielektrizitätskonstante des Vakuums

Da in einem Kristallgitter nicht nur ein Ionenpaar vorhanden ist, sondern im Raum weitere Kat- und Anionen, wird bei der Kristallbildung weitere Energie frei, allerdings auch wieder benötigt, um gleich geladene Ionen anzunähern. Die folgende Gleichung soll dies erläutern:

$E_{IG}=n_1\cdot \frac{1}{c_1}\cdot E_{IP}-n_2\cdot \frac{1}{c_2}\cdot E_{IP}+n_3\cdot \frac{1}{c_3}\cdot E_{IP}-\ldots$

Die im Gitter gespeicherte Energie E_IG ergibt sich dabei als die Summe der bei der Gitterbildung frei gewordenen und benötigten Energien zu jedem Ion. Dabei ist n die Anzahl, wie oft ein bestimmtes Ion vorkommt, und c ein Faktor, der den Abstand des Ions angibt. Diese Faktoren können zu einem, vom Kristall abhängigen Faktor α – der Madelung-Konstante – zusammengefasst werden, so dass sich für die Bindungsenergie eines Ions im Gitter folgende Gleichung ergibt:

$E_{IG}=\alpha\cdot E_{IP}=-\alpha\frac{1}{4\pi\varepsilon_0}\cdot \frac{z^2e^2}{d}$

Diese Gleichung beschreibt die Bindungsenergie nur eines Ions im Gitter. Um die Energie zu erhalten, die bei der Bildung einer bestimmten Stoffmenge n frei wird, muss diese Gleichung noch mit der Avogadro-Konstanten N_A sowie der Stoffmenge n multipliziert werden:

$E_{G}=-\alpha\cdot \frac{N_A}{4\pi\varepsilon_0}\cdot \frac{z^2\cdot e^2\cdot n}{d}$

Wie man sehen kann, ist dieser Wert negativ, da die Gitterbildung exotherm ist. Zur genaueren Berechnung von Gitterenergien reicht die alleinige Betrachtung von regelmäßig angeordneten Coulomb-Punktladungen nicht aus. Eine Erweiterung des Modells führt zur Born-Landé-Gleichung.

Berechnung der Madelung-Konstante für NaCl

Beim Ionengitter von NaCl handelt es sich um eine kubisch flächenzentrierte Kristallstruktur, wie sie rechts abgebildet ist (rot sind die Anionen und grün die Kationen). Der Abstand der beiden Ionen beträgt bei NaCl etwa d = 0,3 nm.

Aus der oben genannten Gleichung müssen nun die Anzahl der jeweiligen benachbarten Ionen n, deren relativer Abstand als Vielfaches vom Abstand d c sowie deren Ladung, also, ob sie sich anziehen oder abstoßen, bestimmt werden.

Wenn wir von einem Ion, z. B. dem rot dargestellten mit der Nummer 0 ausgehen, haben wir so als erstes n=6 Ionen (grün mit der Nummer 1 dargestellt) im Abstand von 1d, die angezogen werden. Danach folgen 12 rote Ionen (2), deren Abstand mittels Satz des Pythagoras den Wert$\sqrt{2}d$ ergibt und die vom gleichartigen Ion abgestoßen werden, sowie 8 grüne Ionen im Abstand von$\sqrt{3}d$, die angezogen werden. Die folgende Tabelle setzt diese Zahlen fort:

Nr.	n	c	Ladung
1	6	$\sqrt{1}=1$	+
2	12	$\sqrt{2}$	-
3	8	$\sqrt{3}$	+
4	6	$\sqrt{4}=2$	-
5	24	$\sqrt{5}$	+
…	…	…	…

Setzt man diese Prozedur fort, so gelangt man zu folgender, in der Literatur häufig anzutreffender Reihendarstellung der Madelung-Konstante^[2]:

$\alpha=\sum_{k=1}^\infty (-1)^{k+1}\frac{n_k}{c_k}=\frac{6}{1}-\frac{12}{\sqrt{2}}+\frac{8}{\sqrt{3}}-\frac{6}{2}+\frac{24}{\sqrt{5}}\pm\ldots\approx 1{,}7476$

Das ist jedoch falsch, weil diese Reihe divergiert, wie erst 1951 bewiesen wurde^[3][4]. Die Summe über die Punkte des Kristallgitters ist bedingt konvergent, hängt also von der Reihenfolge der Summanden ab. Die obige Reihe würde einer Summation über konzentrische Kugelschalen entsprechen, was auch physikalisch nicht sinnvoll ist. Der „richtige“ Wert ergibt sich, indem über die Gitterpunkte innerhalb eines Würfels mit Kantenlänge 2n + 1 summiert und der Grenzwert für $n \to \infty$ gebildet wird. Eine mathematische Begründung, dass dieser Wert der „richtige“ ist, erfolgte 1985 durch David Borwein, Jonathan Borwein und Keith F. Taylor^[4]: Sie definieren eine Funktion einer komplexen Variablen s:

$b(2s) = \sum_{i,j,k=-\infty}^{\infty}\prime {{(-1)^{i+j+k}} \over { (i^2 + j^2 + k^2)^s}}$,

wobei der Strich bedeutet, dass der Term i = j = k = 0 wegzulassen ist. Für hinreichend große s ist die Reihe absolut konvergent, die Madelung-Konstante (s = 1 / 2) ergibt sich durch analytische Fortsetzung.

Die Summation über Würfel konvergiert so langsam, dass sie für praktische Berechnungen unbrauchbar ist. Mit einem kleinen Trick können jedoch wesentlich bessere Näherungen α_n gefunden werden: Es wird über alle Punkte innerhalb eines Würfels mit Kantenlänge 2n + 1 summiert, die Punkte auf den Seiten werden nur halb gezählt, die an den Kanten zu einem Viertel und die an den Ecken zu einem Achtel. Dieses Verfahren ist als Evjen-Methode bekannt^[5]. Bereits α₃ liefert mit 1,7470 einen sehr guten Näherungswert; α₅ = 1,74750 und α₁₀ = 1,7475686. Die physikalische Motivation dieses Verfahrens ergibt sich aus der Forderung, dass die Partialsummen α_n über einen elektrisch neutralen (endlichen) Kristall gebildet werden sollen.

Bereits vor Evjen veröffentlichte Paul Peter Ewald das heute als Ewald-Methode bekannte Verfahren zur Berechnung der Madelung-Konstante^[6]. Seine Methode ist ein Spezialfall der Poissonschen Summenformel, was Ewald jedoch nicht bekannt war.

Heute stehen zahlreiche numerische Methoden und leistungsfähige Computer zur Verfügung, sodass die Berechnung der Madelung-Konstanten für beliebige Gitter mit hoher Genauigkeit kein Problem mehr darstellt.

Weitere Kristallstrukturen

• Für die NaCl-Struktur (kubisch-flächenzentriert) hat die Madelung-Konstante den Wert α = 1,747564594633…
• Für die CsCl-Struktur (kubisch-primitiv) gilt α = 1,7627…
• Für die Zinkblende-Struktur (kubisch-flächenzentriert mit zwei Atomen in der Einheitszelle) gilt α = 1,6381…
• Titan(IV)-oxid α = 2,40…
• Calciumfluorid α = 2,52…
• Kupfer(II)-oxid α = 4,12…
• Saphir α = 4,17…

Einzelnachweise

1. ↑ E. Madelung, Phys. Zs. 19, (1918), 542
2. ↑ Charles Kittel: Einführung in die Festkörperphysik, Oldenbourg, München und Wien 1999, ISBN 3486238434
3. ↑ O. Emersleben: Math. Nachr 4 (1951), 468
4. ↑ ^{a b} D. Borwein, J. M. Borwein, K. F. Taylor: Convergence of Lattice Sums and Madelung’s Constant, J. Math. Phys. 26 (1985), 2999–3009, doi:10.1063/1.526675
5. ↑ H. M. Evjen: On the Stability of Certain Heteropolar Crystals, Phys. Rev. 39 (1932), 675–687, http://link.aps.org/abstract/PR/v39/p675
6. ↑ P. P. Ewald: Die Berechnung optischer und elektrostatischer Gitterpotentiale, Ann. Phys. 64 (1921), 253–287, doi:10.1002/andp.19213690304