Meine Merkliste

Home
Lexikon
Madelung-Konstante

Madelung-Konstante

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₂).

Weiteres empfehlenswertes Fachwissen

Richtiges Wägen mit Laborwaagen: Die Wägefibel

Empfindlichkeitsprüfung von Laborwaagen

Tägliche Sichtkontrolle der Laborwaagen

Inhaltsverzeichnis

1 Berechnung der Bindungsenergie im Gitter
2 Berechnung der Madelung-Konstante für NaCl
3 Weitere Kristallstrukturen
4 Einzelnachweise

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 $2 n + 1$ summiert und der Grenzwert für $n \to \infty$ gebildet wird. Eine mathematische Begründung, dass dieser Wert der „richtige“ ist, erfolgte durch Borwein, Borwein und Taylor: 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 $2 n + 1$ summiert, die Punkte an der Oberfläche 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 $α 3$ liefert mit 1,7470 einen sehr guten Näherungswert; $α 5$ = 1,74750 und $α 10$ = 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...

Einzelnachweise

↑ E. Madelung, Phys. Zs. 19, (1918), 542
↑ Charles Kittel: Einführung in die Festkörperphysik, Oldenbourg, München und Wien 1999, ISBN 3486238434
↑ O. Emersleben: Math. Nachr 4 (1951), 468
↑ 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
↑ H. M. Evjen: On the Stability of Certain Heteropolar Crystals, Phys. Rev. 39 (1932), 675–687, http://link.aps.org/abstract/PR/v39/p675
↑ P. P. Ewald: Die Berechnung optischer und elektrostatischer Gitterpotentiale, Ann. Phys. 64 (1921), 253-287, doi:10.1002/andp.19213690304

Kategorien: Kristallographie | Festkörperphysik

Dieser Artikel basiert auf dem Artikel Madelung-Konstante aus der freien Enzyklopädie Wikipedia und steht unter der GNU-Lizenz für freie Dokumentation. In der Wikipedia ist eine Liste der Autoren verfügbar.