準調和近似

固体物理学における準調和近似（じゅんちょうわきんじ、テンプレート:Lang-en-short）とは、熱膨張のような体積依存性のある熱現象を説明する際にもちいられる、フォノンに基くモデルである。 格子定数のそれぞれを調節可能なパラメータとみなし、かつ調和振動子近似できることを仮定する。

概要

準調和近似は結晶格子の動力学を記述するフォノンモデルを拡張したものである。フォノンモデルは全ての原子間に働く力が調和的であることを仮定するが、すると平衡距離が温度によらず一定になってしまうため熱膨張のような現象を記述することができない。

ここで、準調和近似モデルは、フォノン振動数に体積依存性をくわえ、体積が一定の条件下では調和近似が成り立つものとする。

熱力学

準調和近似のもとでは、ある結晶格子のヘルムホルツエネルギーテンプレート:Mvarは以下のように書ける。

$F(T,V)=E_{\mathrm{l}\mathrm{a}\mathrm{t}}(V)+U_{\mathrm{v}\mathrm{i}\mathrm{b}}(T,V)-TS(T,V)$

ここでテンプレート:Mathは静的な格子内部エネルギー、テンプレート:Mathは格子の振動内部エネルギー、つまりフォノン系のもつエネルギー、テンプレート:Mvarは体積、テンプレート:Mvarは系の振動自由度に由来するエントロピーである。振動エネルギーは以下の式で与えられる。

$U_{\mathrm{v}\mathrm{i}\mathrm{b}}(T,V)=\frac{1}{N}\sum _{𝒌,i}\frac{1}{2}\hslash {\omega }_{𝒌,i}(V)+\frac{1}{N}\sum _{𝒌,i}\frac{\hslash {\omega }_{𝒌,i}(V)}{\exp ({\mathrm{\Theta }}_{𝒌,i}(V)/T)-1}=\frac{1}{N}\sum _{𝒌,i}[\frac{1}{2}+n_{𝒌,i}(T,V)]\hslash {\omega }_{𝒌,i}(V)$

ここで、テンプレート:Mvarは項の総数、${\mathrm{\Theta }}_{𝒌,i}(V)=\hslash {\omega }_{𝒌,i}(V)/k_B$は波数ベクトルテンプレート:Mvarにおけるテンプレート:Mvar番目のフォノンの体積テンプレート:Mvarの条件下での特性温度、$n_{𝒌,i}(T,V)$は温度テンプレート:Mvarおよび体積テンプレート:Mvarにおけるテンプレート:Math-フォノンの数である。慣例のとおり、$\hslash$は換算プランク定数およびテンプレート:Mathはボルツマン定数である。テンプレート:Mathの初項はフォノン系の零点エネルギーであり、熱膨張に零点熱圧力として貢献する。

ヘルムホルツ自由エネルギーテンプレート:Mvarは以下のとおり与えられる。

$F=E_{\mathrm{l}\mathrm{a}\mathrm{t}}(V)+\frac{1}{N}\sum _{𝒌,i}\frac{1}{2}\hslash {\omega }_{𝒌,i}(V)+\frac{1}{N}\sum _{𝒌,i}{k}_{B}T\ln [1-\exp (-{\mathrm{\Theta }}_{𝒌,i}(V)/T)]$

また、エントロピーテンプレート:Mvarは次のように与えられる。

$S=-{\left(\frac{\partial F}{\partial T}\right)}_V=-\frac{1}{N}\sum _{𝒌,i}{k}_B\ln [1-\exp (-{\mathrm{\Theta }}_{𝒌,i}(V)/T)]+\frac{1}{NT}\sum _{𝒌,i}\frac{\hslash {\omega }_{𝒌,i}(V)}{\exp ({\mathrm{\Theta }}_{𝒌,i}(V)/T)-1}$

簡単に確かめられるとおり、これらはテンプレート:Mathを満たす。

テンプレート:Mvarの関数としての周波数テンプレート:Mvarは分散関係と呼ばれる。体積テンプレート:Mvarを一定値に保つ条件の下では、これらの方程式は調和近似に相当することに注意されたい。

ルジャンドル変換を適用することにより系のギブズ自由エネルギーテンプレート:Mvarを温度と圧力の関数として書くこともできる。

$G(T,P)=\underset{V}{\min }[E_{\mathrm{l}\mathrm{a}\mathrm{t}}(V)+U_{\mathrm{v}\mathrm{i}\mathrm{b}}(V,T)-TS(T,V)+PV]$

ここでテンプレート:Mvarは圧力である。テンプレート:Mvarは所与のテンプレート:Mvarおよびテンプレート:Mvarにおける平衡体積において最小値をとる。

派生物理量

ギブズ自由エネルギーが得られたならば、そこから派生して多くの熱力学量が得られる。以下に、調和近似だけでは決定することのできない物理量をいくつか示す。

平衡体積

圧力および温度の関数としての体積テンプレート:Mathはギブズ自由エネルギーを最低化することにより得られる。

熱膨張係数

熱膨張係数テンプレート:Mathはテンプレート:Mathから以下のように得られる。

${\alpha }_V=\frac{1}{V}{\left(\frac{\partial V}{\partial T}\right)}_P$

グリュナイゼン定数

グリュナイゼン定数テンプレート:Mvarは各フォノンモードごとに以下のように定義される。

${\gamma }_i=-\frac{\partial \ln {\omega }_i}{\partial \ln V}$

ここでテンプレート:Mvarはフォノンモードを表わす添字である。総グリュナイゼン定数はすべてのテンプレート:Mvarの総和として得られる。グリュナイゼン定数は系の非調和性の尺度であり、熱膨張に密接に関連する。

参照文献

• Dove, Martin T. (1993). Introduction to lattice dynamics, Cambridge university press. テンプレート:ISBN2.