晶体衍射与Ш函数 [学习笔记]

本文先证明一个重要的引理：泊松求和公式，然后给出一个实例：一维晶体衍射，并在计算过程中引入Ш函数，并介绍Ш函数的一些性质。

首先，我们先证明泊松求和公式：

若为速降函数，则

$\begin{matrix} (1) & \sum_{k = - \infty}^{\infty} φ (k) = \sum_{k = - \infty}^{\infty} F φ (k) \end{matrix}$ $\sum_{k=-\infty}^{\infty}\varphi(k) = \sum_{k=-\infty}^{\infty}{\scr F}\varphi(k) \tag1$

证明如下：

我们先将函数 $\varphi(x)$ 周期化，即在一维直线上每隔 $k$ 放一个 $\varphi$ ，这样就可以得到

\begin{matrix} (2) & Φ (x) = \sum_{k = - \infty}^{\infty} φ (x - k) \end{matrix}

$\Phi(x) = \sum_{k=-\infty}^{\infty}\varphi(x-k)\tag 2$

接着，我们将 $\Phi(x)$ 用傅里叶级数展开，有：

\begin{matrix} (3) & Φ (x) = \sum_{k = - \infty}^{\infty} \hat{Φ} (k) e^{2 π i k x} \end{matrix}

$\Phi(x) = \sum_{k=-\infty}^{\infty} \hat{\Phi}(k) e^{2\pi ikx}\tag 3$

其中，

\hat{Φ} (s) = \int_{0}^{1} Φ (x) e^{- 2 π i s x} d x

$\hat{\Phi}(s) = \int_{0}^{1}\Phi(x) e^{-2\pi isx}dx$

我们将 $(2)$ 式带入其中，并令 $u=x-k$ ，则有

\begin{aligned} \hat{Φ} (s) & = \int_{0}^{1} \sum_{k = - \infty}^{\infty} φ (x - k) e^{- 2 π i s x} d x \\ = \sum_{k = - \infty}^{\infty} \int_{- k}^{1 - k} φ (u) e^{- 2 π i s (u + k)} d u \\ = \sum_{k = - \infty}^{\infty} \int_{- k}^{1 - k} φ (u) e^{- 2 π i s u} e^{- 2 π i s k} d u \end{aligned}

$\begin{align*} \hat{\Phi}(s) &= \int_{0}^{1} \sum_{k=-\infty}^{\infty}\varphi(x-k) e^{-2\pi isx}dx\\ &= \sum_{k=-\infty}^{\infty}\int_{-k}^{1-k}\varphi(u) e^{-2\pi is(u+k)}du\\ &=\sum_{k=-\infty}^{\infty}\int_{-k}^{1-k}\varphi(u) e^{-2\pi isu} e^{-2\pi isk}du \end{align*}$

注意到傅里叶级数的原函数周期为 $1$ ，因此积分区间 ${k,1+k}$ 对于任何实数 $k$ 积分的值都是相同的，同时，因为是级数展开式，因此式中的 $s$ 和 $k$ 都是整数，于是 $e^{-2\pi isk}=1$ 。最后，我们注意到， $k$ 的取值从 $-\infty$ 到 $\infty$ ，相当于积分上下限的增加，最后将增加到 $(-\infty, \infty)$ 。于是，我们得到：

\hat{Φ} = F φ

$\hat{\Phi} = {\scr F}\varphi$

将其代回 $(3)$ 式，我们有

\begin{matrix} (4) & Φ (x) = \sum_{k = - \infty}^{\infty} F φ (k) e^{2 π i k x} \end{matrix}

$\Phi(x) = \sum_{k=-\infty}^{\infty} {\scr F}\varphi(k) e^{2\pi ikx}\tag 4$

使 $(2)$ 、 $(4)$ 两式中的 $k=0$ ，可得

Φ (0) = \sum_{k = - \infty}^{\infty} φ (- k) = \sum_{k = - \infty}^{\infty} φ (k)

$\Phi(0) = \sum_{k=-\infty}^{\infty} \varphi(-k) = \sum_{k=-\infty}^{\infty} \varphi(k)$

Φ (0) = \sum_{k = - \infty}^{\infty} F φ (k)

$\Phi(0) = \sum_{k=-\infty}^{\infty} {\scr F}\varphi(k)$

对比求和项，得证

接下来，我们给出一个分布傅里叶变换应用的实例：一维晶体衍射

方便起见，我们研究一维情况，即原子分布在一条无限上的直线上

这里写图片描述

我们设单个原子的密度为 $\rho(x)$ ，考虑到量子力学，其实这就是概率密度。接着，我们选定某个原子作为零点，设这些原子两两间隔为 $a$ ，则整个一维直线上的原子的密度可以表示为所有位置的原子的密度之和，即

ρ_{a} (x) = \sum_{k = - \infty}^{\infty} ρ (x - k a)

$\rho_a(x) = \sum_{k=-\infty}^{\infty} \rho(x-ka)$

注意到右边可以表示成， $\delta$ 函数的卷积形式，有

\begin{matrix} (5) & ρ_{a} (x) = \sum_{k = - \infty}^{\infty} ρ (x) * δ (x - k a) = ρ (x) * \sum_{k = - \infty}^{\infty} δ (x - k a) \end{matrix}

$\rho_a(x) = \sum_{k=-\infty}^{\infty} \rho(x) * \delta(x-ka) = \rho(x) * \sum_{k=-\infty}^{\infty} \delta(x-ka) \tag5$

此时，我们引入Ш函数，我们定义：

$Ш_{a} (x) = \sum_{k = - \infty}^{\infty} δ (x - k a)$ $Ш_a(x) = \sum_{k=-\infty}^{\infty} \delta(x-ka)$

$(5)$ 式就可以表示为

ρ_{a} (x) = ρ (x) * Ш_{a} (x)

$\rho_a(x) = \rho(x) * Ш_a(x)$

两边取傅里叶变换，有：

\begin{matrix} (6) & F ρ_{a} (x) = (F ρ (x)) (F Ш_{a} (x)) \end{matrix}

${\scr F} \rho_a(x) = \left({\scr F} \rho(x) \right) \left({\scr F} Ш_a(x)\right) \tag 6$

想解出 $(6)$ 式，就必须知道Ш函数的傅里叶变换，我们先看Ш函数的傅里叶变换

首先，当 $a=0$ 时，

Ш = \sum_{k = - \infty}^{\infty} δ (x - k)

$Ш = \sum_{k=-\infty}^{\infty} \delta(x-k)$

于是，

\begin{aligned} F Ш = ⟨ F Ш, φ ⟩ & = ⟨ Ш, F φ ⟩ \\ = ⟨ \sum_{k = - \infty}^{\infty} δ (x - k), F φ ⟩ \\ = \sum_{k = - \infty}^{\infty} F φ (k) \end{aligned}

$\begin{align*} {\scr F}Ш = \left<{\scr F}Ш, \varphi \right> &= \left<Ш, {\scr F}\varphi \right>\\ &= \left< \sum_{k=-\infty}^{\infty} \delta(x-k), {\scr F}\varphi \right>\\ &=\sum_{k=-\infty}^{\infty}{\scr F}\varphi(k) \end{align*}$

由泊松求和公式 $(1)$ 式，可得

\begin{aligned} ⟨ Ш, F φ ⟩ & = \sum_{k = - \infty}^{\infty} φ (k) \\ = ⟨ Ш, φ ⟩ \\ = Ш \end{aligned}

$\begin{align*} \left<Ш, {\scr F}\varphi \right> &=\sum_{k=-\infty}^{\infty} \varphi(k)\\ &=\left<Ш, \varphi \right>\\ &=Ш \end{align*}$

于是我们就得到了

\begin{matrix} (7) & F Ш = Ш \end{matrix}

${\scr F}Ш = Ш\tag 7$

接下来看更普通的情况，一下过程用到了 $\delta(ax)$ 的性质

\begin{aligned} Ш_{a} (x) & = \sum_{k = - \infty}^{\infty} δ (x - k a) \\ = \sum_{k = - \infty}^{\infty} δ (a (\frac{x}{a} - k)) \\ = \frac{1}{a} \sum_{k = - \infty}^{\infty} δ (\frac{x}{a} - k) \\ = \frac{1}{a} Ш (\frac{x}{a}) \end{aligned}

$\begin{align*} Ш_a(x)&=\sum_{k=-\infty}^{\infty} \delta(x-ka)\\ &=\sum_{k=-\infty}^{\infty} \delta \left( a \left(\frac{x}{a}-k \right) \right)\\ &=\frac{1}{a} \sum_{k=-\infty}^{\infty} \delta \left( \frac{x}{a}-k \right)\\ &=\frac{1}{a} Ш( \frac{x}{a} ) \end{align*}$

那么

\begin{aligned} F Ш_{a} (x) = F (\frac{1}{a} Ш (\frac{x}{a})) & = \frac{1}{a} F (Ш (\frac{x}{a})) \\ = \frac{1}{a} (a (F Ш) (a x)) \\ = Ш (a x) \\ = \sum_{k = - \infty}^{\infty} δ (a x - k) \\ = \sum_{k = - \infty}^{\infty} δ (a (x - \frac{k}{a})) \\ = \frac{1}{a} \sum_{k = - \infty}^{\infty} δ (x - \frac{1}{a} k) \\ = \frac{1}{a} Ш_{\frac{1}{a}} \end{aligned}

$\begin{align*} {\scr F}Ш_a(x) = {\scr F} \left( \frac{1}{a} Ш( \frac{x}{a} ) \right) &= \frac{1}{a} {\scr F} \left( Ш( \frac{x}{a} ) \right)\\ &=\frac{1}{a} \left(a \left( {\scr F} Ш \right) ( ax ) \right)\\ &=Ш(ax)\\ &=\sum_{k=-\infty}^{\infty} \delta (ax-k)\\ &=\sum_{k=-\infty}^{\infty} \delta \left(a \left(x-\frac{k}{a} \right) \right)\\ &=\frac{1}{a}\sum_{k=-\infty}^{\infty} \delta \left(x-\frac{1}{a}k \right)\\ &=\frac{1}{a} Ш _ {\frac{1}{a}} \end{align*}$

即

\begin{matrix} (8) & F Ш_{a} = \frac{1}{a} Ш_{\frac{1}{a}} \end{matrix}

${\scr F}Ш_a = \frac{1}{a} Ш _ {\frac{1}{a}}\tag8$

$Ш_a$ 与 ${\scr F}Ш_a$ 的部分图像如上图所示，可见他们的间隔互为倒数， ${\scr F}Ш_a$ 的高度是 $Ш_a$ 的 $\frac{1}{a}$ 倍（虽然其实箭头指代的都是正无穷）

现在我们回到待解决的 $(6)$ 式，我们将 $(8)$ 式代入 $(6)$ 式中，即可得到：

\begin{aligned} F ρ_{a} & = (F ρ) (\frac{1}{a} Ш_{\frac{1}{a}}) \\ = \frac{1}{a} (F ρ) \sum_{k = - \infty}^{\infty} δ (x - \frac{1}{a} k) \\ = \frac{1}{a} \sum_{k = - \infty}^{\infty} F ρ (\frac{k}{a}) δ (x - \frac{k}{a}) \end{aligned}

$\begin{align*} {\scr F} \rho_a &= \left({\scr F} \rho \right) \left( \frac{1}{a} Ш _ {\frac{1}{a}} \right)\\ &=\frac{1}{a} \left({\scr F} \rho \right) \sum_{k=-\infty}^{\infty} \delta \left(x-\frac{1}{a}k \right)\\ &=\frac{1}{a} \sum_{k=-\infty}^{\infty} {\scr F} \rho \left( \frac{k}{a} \right) \delta \left(x-\frac{k}{a} \right) \end{align*}$

从衍射的公式可以看出，一维晶格的衍射图案为一系列间距为 $\frac{1}{p}$ 的点，各点的强度如上式所示

最后，我们给出Ш函数的一些常用性质