フォノンとは - 理数の散策路

フォノンとは、結晶内の格子振動による弾性波を量子化したものです。

古典論の弾性波
量子論の弾性波

古典論の弾性波

古典的な弾性波を導きます。簡単のため、間隔 $a$ で並ぶ１次元の $N$ 個の原子配列（格子）を考えます。ある原子に働く力は両隣りの原子との距離に比例すると仮定すると以下で表されます。

$F_{s} = C (q_{s + 1} - q_{s}) + C (q_{s - 1} - q_{s}) - ①$

ここで、 $q_{s}$ と $m$ は原子 $s$ の変位と質量、 $C$ は比例定数（弾性定数）を表します。 $q_{s}$ は以下で定義します。

$q_{s} \propto e^{i (s k a - ω t)} - ②$

また、 $k$ は波数で以下で定義します。 $n$ は波の数で、 $n$ が大きくなるほど波長は短くなります。原子は間隔 $a$ で並んでいるため、意味のある最も短い波長は $2 a$ （ $k = 2 π / λ = π / a$ ）となります。

$k = \frac{2 π n}{N a}, n = 0, \pm 1, \pm 2, \dots, \pm \frac{N}{2}$

このとき、分散関係は、

$ω_{k}^{2} = \frac{2 C}{m} (1 - \cos k a) - ③$ $= \frac{4 C}{m} \sin^{2} \frac{1}{2} k a - ④$

群速度は④から求められます。

$v_{g} \equiv \frac{d ω_{k}}{d k} = \sqrt{\frac{C a^{2}}{m}} \cos \frac{1}{2} k a - ⑤$

古典的ハミルトニアンは以下で表すことができます。

$H = \sum_{s = 1}^{N} (\frac{1}{2 m} p_{s}^{2} + \frac{C}{2} (q_{s + 1} - q_{s})^{2}) - ⑥$ $= \sum_{k} (\frac{1}{2 m} P_{k}^{2} + \frac{m ω_{k}^{2}}{2} Q_{k}^{2}) - ⑦$

ここで、座標と運動量を離散フーリエ変換した基準座標 $Q_{s}$ と共役運動量 $P_{k}$ を以下で定義します。

$Q_{k} = \frac{1}{\sqrt{N}} \sum_{s} q_{s} e^{- i k s a} - ⑧$ $P_{k} = \frac{1}{\sqrt{N}} \sum_{s} p_{s} e^{- i k s a} － ⑧$

③を導く

①より運動方程式は以下になります。

$m \frac{d^{2} q_{s}}{d t^{2}} = F_{s} = C (q_{s + 1} + q_{s - 1} - 2 q_{s}) - (1)$

これに②を代入すると、

$- ω_{k}^{2} m q_{s} = C (q_{s + 1} + q_{s - 1} - 2 q_{s})$

さらに②より $q_{s \pm 1} = q_{s} e^{\pm i k a}$ であるから、

$- ω_{k}^{2} m = C (e^{i k a} + e^{- i k a} - 2) = 2 C (\cos k a - 1)$

これより③が導かれます。

⑥を導く

⑥をハミルトン方程式に代入すると、 $q_{s}$ を含まない項は０になるため、

$\frac{d p_{s}}{d t} = - \frac{\partial H}{\partial q_{s}}$ $= C (q_{s + 1} + q_{s - 1} - 2 q_{s})$

これより、運動方程式 (1) が導かれることが分かります。

⑦を導く

⑥の第１項に⑧の逆変換

$p_{s} = \frac{1}{\sqrt{N}} \sum_{k} P_{k} e^{i k s a}$

を代入すると、

$\sum_{s = 1}^{N} \frac{1}{2 m} p_{s}^{2} = \sum_{s = 1}^{N} \frac{1}{2 m} \frac{1}{N} (\sum_{k} P_{k} e^{i k s a})^{2}$ $= \frac{1}{2 m} \sum_{k} P_{k}^{2}$

一方、⑥の第２項について、 $q_{s}$ が含まれる項は以下になりますが、

$\frac{C}{2} (q_{s + 1} - q_{s})^{2} + \frac{C}{2} (q_{s} - q_{s - 1})^{2}$ $= \frac{C}{2} (q_{s + 1}^{2} - 2 q_{s + 1} q_{s} + 2 q_{s}^{2} - 2 q_{s} q_{s - 1} + q_{s - 1}^{2})$

$q_{s + 1} q_{s}$ と $q_{s} q_{s - 1}$ の項はそれぞれ $q_{s + 1}$ と $q_{s - 1}$ へも寄与しているので、それらは $\times 1 / 2$ とすると、 $q_{s}$ に寄与する項は以下になります。

$\frac{C}{2} (- q_{s + 1} q_{s} + 2 q_{s}^{2} - q_{s} q_{s - 1})$

これに⑧の逆変換

$q_{s} = \frac{1}{\sqrt{N}} \sum_{k} Q_{k} e^{i k s a}$

を代入すると、

$\sum_{s = 1}^{N} \frac{C}{2} (- q_{s + 1} q_{s} + 2 q_{s}^{2} - q_{s} q_{s - 1})$

$= \sum_{s = 1}^{N} \frac{C}{2} \frac{1}{N} (\sum_{k} Q_{k} e^{i k s a})^{2} (- e^{i k a} - e^{- i k a} + 2)$ $= \sum_{k} C Q_{k}^{2} (1 - \cos k a) = \sum_{k} \frac{m ω_{k}^{2}}{2} Q_{k}^{2}$