[論文レビュー] Lattice vibrations in the harmonic approximation

G. V. Paradezhenko, Nikolai B. Melnikov+1 more|arXiv (Cornell University)|Nov 12, 2017

Electrostatics and Colloid Interactions参考文献 3被引用数 3

ひとこと要約

本稿は、正規モードの量子化と原子の振動のガウス統計を用いて、調和近似におけるデバイ・ワラー因子の簡潔な理論的導出を提供する。デバイモデルにおけるデバイ・ワラー因子の計算に適した効率的な公式を導出し、高温においては $ 2W(\boldsymbol{k}) \approx \frac{3\hbar^{2}\kappa^{2}T}{M\Theta_{\mathrm{D}}^{2}} $ とスケーリングすることを示した。ここで $ \Theta_{\mathrm{D}} $ はデバイ温度である。

ABSTRACT

We present some theoretical results on the lattice vibrations that are necessary for a concise derivation of the Debye-Waller factor in the harmonic approximation. First we obtain an expression for displacement of an atom in a crystal lattice from its equilibrium position. Then we show that an atomic displacement has the Gaussian distribution. Finally, we obtain the computational formula for the Debye-Waller factor in the Debye model.

研究の動機と目的

• デバイ・ワラー因子を導出するために必要な、調和近似における格子振動の簡潔な理論的基盤を提供すること。
• 正準平均としての二乗原子変位の計算がガウス確率密度を用いて行えることを確立し、量子統計的取り扱いを簡略化すること。
• 金属におけるフォノン補正の計算に用いることを見越して、デバイモデルにおけるデバイ・ワラー因子の計算に適した効率的な公式を導出すること。
• メルミンの線形演算子の指数平均に対する公式が、複雑な量子統計的計算を回避して特性関数を一言で導出可能であることを示すこと。

提案手法

• ハミルトニアン形式を用いて、古典的正規モードから出発し、それらの第二量子化を行う。
• 変位に関する線形演算子の指数関数の正準平均を扱うためにメルミンの公式を適用し、原子変位の特性関数を一言で導出する。
• 得られた原子変位の分布がガウス分布であることを示し、統計的平均の簡略化を可能にする。
• 線形分散関係 $ \omega = vq $ を仮定することでデバイモデルを導入し、$ 0 \leq \omega \leq \omega_{\mathrm{D}} $ の範囲でフォノン状態密度 $ n(\omega) = \frac{9N\omega^{2}}{\omega_{\mathrm{D}}^{3}} $ を導出する。
• 状態密度を用いて正規モードの和を周波数に関する積分に変換し、デバイ・ワラー因子を周波数 $ \omega $ に関する積分で表す式を得る。
• $ \coth(x) \approx 1/x $ を小刻みな $ x $ に対して用いて、高温極限における積分を評価し、デバイ・ワラー因子の閉形式の式を得る。

実験結果

リサーチクエスチョン

• RQ1 長時間の平均二乗変位の計算を避け、調和近似においてデバイ・ワラー因子を効率的に導出する方法は何か？
• RQ2 ガウス分布は、結晶格子内での原子変位の統計的記述において果たす役割は何か？
• RQ3 メルミンの公式は、原子変位の特性関数の導出をどのように簡略化するか？
• RQ4 デバイモデルにおけるデバイ・ワラー因子の解析的形は何か？また、高温における挙動はいかなるものか？
• RQ5 デバイモデルにおけるフォノン状態密度とデバイ周波数の関係は何か？

主な発見

• 調和近似における正準平均としての二乗原子変位は、ガウス確率密度を用いた平均と等価であり、統計的取り扱いが簡略化可能である。
• メルミンの公式を用いることで、原子変位の特性関数を一言で導出でき、複雑な量子統計的変形を回避できる。
• デバイモデルにおけるデバイ・ワラー因子は $ 2W(\kappa) = \frac{3\hbar\kappa^{2}}{2M\omega_{\mathrm{D}}^{3}} \int_{0}^{\omega_{\mathrm{D}}} \coth\left(\frac{\hbar\omega}{2T}\right) \omega \, d\omega $ と表され、計算に適している。
• 高温（$ T \gg \Theta_{\mathrm{D}} $）においては、デバイ・ワラー因子は $ 2W(\kappa) \approx \frac{3\hbar^{2}\kappa^{2}T}{M\Theta_{\mathrm{D}}^{2}} $ に簡略化され、ここで $ \Theta_{\mathrm{D}} = \hbar\omega_{\mathrm{D}} $ はデバイ温度である。
• デバイモデルにおけるフォノン状態密度は、$ 0 \leq \omega \leq \omega_{\mathrm{D}} $ の範囲で $ n(\omega) = \frac{9N\omega^{2}}{\omega_{\mathrm{D}}^{3}} $ であり、球対称のブリユアンゾーンと線形分散関係の仮定から導出される。
• 導出結果から、特に局所的な磁気的性質に対するフォノン補正を推定するにあたり、調和近似が十分であることが確認された。