Debye-Waller-Faktor

Der Debye-Waller-Faktor (DWF, nach Peter Debye und Ivar Waller) beschreibt, wie die Intensität der an einem Kristallgitter kohärent elastisch gestreuten Strahlung von der Temperatur abhängt. Nur diese elastische Streuung unterliegt den Laue-Bedingungen; die komplementäre, inelastische Streuung wird als thermisch-diffus bezeichnet.

Durch thermische Bewegung der Atome werden die Reflexe der elastischen Streuung nicht verbreitert, sondern ihre Intensität herabgesetzt. Es erscheint allerdings ein diffuser Untergrund zwischen den Reflexen als Folge der Energieerhaltung.

In der Neutronenstreuung wird der Begriff Debye-Waller-Faktor teilweise unterschiedslos auf kohärente und inkohärente Streuung angewandt; teilweise wird für letztere aber auch der genauere Begriff Lamb-Mößbauer-Faktor benutzt.

Definition

Die Intensität $I_0$ der einfallenden Welle wird durch Multiplikation mit dem Debye-Waller-Faktor auf die Intensität $I$ der gestreuten Welle reduziert, und zwar um den Faktor $(1-DWF)$:

$I=I_0\cdot \underset{DWF<1}{\underbrace{\exp (-\frac{1}{3}{\left|\overrightarrow{G}\right|}^2\overline{u^2})}}$

mit

• der natürlichen Exponentialfunktion $\exp ()$
• einem Gittervektor $G$ des reziproken Gitters
• der temperaturabhängigen Oszillations­amplitude $u=u(T)$ der Atome.

Die Bragg-Beugungsreflexe werden also aufgrund der Gitterschwingungen umso mehr gedämpft, je höher die Temperatur und je höher ihre Ordnung ist:

• Der DWF ist maximal, wenn die Atome in der Nähe des absoluten Nullpunkts nicht schwingen (entspricht dem statischen Fall):

$T\approx 0K\Rightarrow u\approx 0\Rightarrow DWF\approx 1$.

• Bei größerer Temperatur wird $\overline{|\overrightarrow{u}{|}^2}$ größer und somit der Exponentialfaktor kleiner.
• Der DWF und somit die Reflex-Intensität ist außerdem umso kleiner, je größer $|\overrightarrow{G}{|}^2$ ist, also je höher die Millerschen Indizes der Netzebenenschar sind, an der die Bragg-Reflexion stattfindet.

Bei Betrachtung eines harmonischen Oszillators mit der Energie:

$\begin{array}{rlrl}\bar{E}& =\frac{1}{2}M{\omega }^2& & \overline{u^2}=\frac{3}{2}{k}_{b}T\\ & \iff & & \overline{u^2}=\frac{3k_{b}T}{M{\omega }^2}\end{array}$

mit

• der Boltzmann-Konstante $k_b$

lässt sich der temperaturabhängige Debye-Waller-Faktor auch schreiben als:

$DWF=\exp (-\frac{k_{b}T{\left|\overrightarrow{G}\right|}^2}{M{\omega }^2})$

Herleitung

Der Strukturfaktor $F_{hkl}$ ist ein Maß für die relative Intensität eines durch die Millerschen Indizes $h$, $k$, $l$ bestimmten Beugungsreflexes:

$F_{hkl}=\sum _i{f}_i\exp [i\overrightarrow{G}\cdot {\overrightarrow{r}}_i]$

Die Summe läuft über alle Atome der Basis. Dabei ist

• $i$ die imaginäre Einheit
• $\overrightarrow{G}=h{\overrightarrow{b}}_1+k{\overrightarrow{b}}_2+l{\overrightarrow{b}}_3$ ein reziproker Gittervektor
• ${\overrightarrow{r}}_i$ ein Ortsvektor, der von einem festen Bezugspunkt innerhalb der Elementarzelle zum Kern des $i$-ten Atom zeigt
• $f_i$ der atomare Streufaktor des $i$-ten Atoms:

$f_i=\underset{V_{A_i}}{\int }{n}_i(\overrightarrow{\tilde{r}})\exp [i\overrightarrow{G}\cdot \overrightarrow{\tilde{r}}]{\mathrm{d}}^3\tilde{r}$

• $V_{A_i}$ das Volumen
• $n_i(\overrightarrow{\tilde{r}})$ das Streuvermögen (z. B. Elektronendichte bei Röntgenbeugung, Ladungsdichte bei Elektronenbeugung) des $i$-ten Atoms.

Betrachtet man die thermische Bewegung der Atome, so ist ${\overrightarrow{r}}_i$ zeitabhängig. Nun zerlegt man ${\overrightarrow{r}}_i$ in einen mittleren Aufenthaltsort ${\overrightarrow{r}}_{i,0}$ (Gleichgewichtslage, ruhend) und die Auslenkung ${\overrightarrow{u}}_i(t)$ (zeitabhängig):

${\overrightarrow{r}}_i(t)={\overrightarrow{r}}_{i,0}+{\overrightarrow{u}}_i(t)$

Die Schwingungsperioden sind sehr kurz ($<10^{-10}$s) gegenüber der Beobachtungsdauer, sodass immer ein zeitlicher Mittelwert gemessen wird:

$\begin{array}{rl}\overline{F_{hkl}}& =\sum _i{f}_i\overline{\exp [i\overrightarrow{G}\cdot ({\overrightarrow{r}}_{i.0}+{\overrightarrow{u}}_i(t))]}\\ & =\sum _i{f}_i\exp [i\overrightarrow{G}\cdot {\overrightarrow{r}}_{i.0}]\overline{\exp [i\overrightarrow{G}\cdot {\overrightarrow{u}}_i(t)]}\end{array}$

Für kleine Auslenkungen entwickelt man die Exponentialfunktion bis zur zweiten Ordnung

$\overline{\exp [i\overrightarrow{G}\cdot {\overrightarrow{u}}_i(t)]}\approx 1+i\overline{\overrightarrow{G}\cdot {\overrightarrow{u}}_i(t)}-\frac{1}{2}\stackrel{‾}{{(\overrightarrow{G}\cdot {\overrightarrow{u}}_i(t))}^2}$

Die erste Ordnung verschwindet $(\overline{\overrightarrow{G}\cdot {\overrightarrow{u}}_i(t)}=0)$, da die Auslenkungen ${\overrightarrow{u}}_i(t)$ statistisch in alle Raumrichtungen erfolgen $(\overline{{\overrightarrow{u}}_i(t)}=0)$ und nicht mit der Richtung von $\overrightarrow{G}$ korreliert sind.

Die zweite Ordnung ist

$\stackrel{‾}{{(\overrightarrow{G}\cdot {\overrightarrow{u}}_i(t))}^2}=|\overrightarrow{G}{|}^2\overline{|{\overrightarrow{u}}_i(t){|}^2}\overline{{\cos }^2\theta }$

Dabei ist $\theta$ der Winkel zwischen $\overrightarrow{G}$ und ${\overrightarrow{u}}_i(t)$. Man mittelt ${\cos }^2\theta$ über alle Richtungen im dreidimensionalen Raum, also Integration über die Einheitskugel:

$\overline{{\cos }^2\theta }=\frac{\underset{0}{\overset{2\pi }{\int }}\mathrm{d}\varphi \underset{0}{\overset{\pi }{\int }}\mathrm{d}\theta \sin \theta {\cos }^2\theta }{\underset{0}{\overset{2\pi }{\int }}\mathrm{d}\varphi \underset{0}{\overset{\pi }{\int }}\mathrm{d}\theta \sin \theta }=\frac{1}{3}$

In die Exponentialfunktion eingesetzt ergibt dies:

$\begin{array}{rl}\overline{\exp [i\overrightarrow{G}\cdot {\overrightarrow{u}}_i(t)]}& \approx 1-\frac{1}{6}|\overrightarrow{G}{|}^2\overline{|{\overrightarrow{u}}_i(t){|}^2}\\ & \approx \exp [-\frac{1}{6}|\overrightarrow{G}{|}^2\overline{|{\overrightarrow{u}}_i(t){|}^2}]\end{array}$

Der Strukturfaktor schreibt sich nun:

$\overline{F_{hkl}}=\sum _i{f}_i\exp [i\overrightarrow{G}\cdot {\overrightarrow{r}}_{i.0}]\exp [-\frac{1}{6}|\overrightarrow{G}{|}^2\overline{|{\overrightarrow{u}}_i(t){|}^2}]$

Für gleichartige Atome ist $\overline{|{\overrightarrow{u}}_i(t){|}^2}$ für alle $i$ annähernd gleich $\overline{|\overrightarrow{u}{|}^2}$. Somit kann man den zweiten Exponentialfaktor vor die Summe ziehen:

$\begin{array}{rl}\overline{F_{hkl}}& =\exp [-\frac{1}{6}|\overrightarrow{G}{|}^2\overline{|\overrightarrow{u}(t){|}^2}]\sum _i{f}_i\exp [i\overrightarrow{G}\cdot {\overrightarrow{r}}_{i.0}]\\ & =\exp [-\frac{1}{6}|\overrightarrow{G}{|}^2\overline{|\overrightarrow{u}(t){|}^2}]F_{hkl}^0\end{array}$

Darin ist $F_{hkl}^0$ der Strukturfaktor des statischen Falls (starres Gitter, keine Bewegung der Atome).

Die Intensität ist proportional zum Betragsquadrat des Strukturfaktors: $I=c|F_{hkl}{|}^2$. Die zeitlich gemittelte Intensität ist somit

$\begin{array}{rl}\bar{I}& =c|\overline{F_{hkl}}{|}^2\\ & =c|F_{hkl}^0{|}^2\exp [-\frac{1}{3}|\overrightarrow{G}{|}^2\overline{|\overrightarrow{u}{|}^2}]\\ & =I_0\exp [-\frac{1}{3}|\overrightarrow{G}{|}^2\overline{|\overrightarrow{u}{|}^2}]\end{array}$

Die gemittelte Intensität ist also gegenüber dem statischen Fall $I_0=c|F_{hkl}^0{|}^2$ um den Debye-Waller-Faktor $\exp [-\frac{1}{3}|\overrightarrow{G}{|}^2\overline{|\overrightarrow{u}{|}^2}]$ erniedrigt.