tecktom.com

この章のまとめ

2026年10月15日8分で読めますAI生成教科書『虎の巻』第3部 第2章 比熱と熱伝導:調和近似の外へ

AI 生成・木下が未勉強

ハミルトニアンの選び方で固体を考える 第3部 第2章 §3.2.7 ★★★ 3次の項1つから、高温の比熱のずれ・熱膨張・有限の熱伝導率がまとめて出る。熱抵抗には結晶運動量を格子へ渡す U 過程が要る

この章では、原子間ポテンシャルを3次まで(比熱では4次まで)展開したハミルトニアン (3.2.1)・(3.2.2) を選んだ。3次の項は、伸ばす向きと縮める向きでポテンシャルを非対称にし、基準座標で書くと3つのフォノンの波数の和を逆格子ベクトルの分だけずらしてよい形 (3.2.5) になった。古典のカノニカル分布の累積展開を当てると、比熱には kBT/Ebk_{\mathrm B} T/E_{\mathrm b} に比例する補正 (3.2.9) が付き(3次の項は増やし、硬くなる向きの4次の項は減らす)、結合の平均の長さは温度に比例して伸びた (3.2.11)。3次の項を「ばね定数が間隔に依存すること」と読み替えた準調和近似からは、熱膨張率が比熱と Grüneisen 定数の積に比例するという Grüneisen の関係 (3.2.16) が出て、1次元の鎖では直接の計算と一致した。熱伝導は、フォノンの Boltzmann 方程式を緩和時間近似で解いて κ=13Cvℓ\kappa=\frac13Cv\ell の形になり、平均自由行程を有限にする散乱のうち、3フォノン過程の N 過程は全結晶運動量を保存するため熱抵抗を生じさせず、U 過程だけが流れを止めることを示した。これから、高温の 1/T1/T、中間の温度の指数関数的な増加、低温の境界散乱による T3T^3 という熱伝導率の3つの領域が出た。

捨てたものは3つで、展開の高次の項(B 型。展開は漸近的で、解離のため分配関数そのものは発散する)、一定体積での非調和性(準調和近似。熱膨張率は最低次で正しいが、比熱の補正は Morse の鎖で 12.5%12.5\% ずれる)、衝突の項の中身(フォノンの Boltzmann 方程式の導出は未解決で、緩和時間近似と Matthiessen の規則は C 型)である。次の章(第3部第3章)では、ここで使った分散関係 ωs(k)\omega_s(\bm{k}) とその圧力依存を、中性子と Raman 散乱でどう測るかを扱う。

総合問題 Morse の鎖の熱膨張を絶対零度から高温までたどる ★★★

オリジナル質量 M=1.0×10−25 kgM=1.0\times10^{-25}\,\mathrm{kg} の原子が間隔 a=2.5×10−10 ma=2.5\times10^{-10}\,\mathrm m で1列に並び、隣どうしが Morse ポテンシャル (3.2.3)(Eb=0.4 eVE_{\mathrm b}=0.4\,\mathrm{eV}、b=1.5×1010 m−1b=1.5\times10^{10}\,\mathrm{m^{-1}})で結合した鎖を考える。鎖は張力 0 に保たれているとする。

  1. (1) ばね定数 KK と、フォノンの最大の角振動数 ωm=2K/M\omega_{\mathrm m}=2\sqrt{K/M}、それに対応する温度 ℏωm/kB\hbar\omega_{\mathrm m}/ k_{\mathrm B} を求めよ。(使う道具:3.2.1 節、解析力学から様々な理論分野を組み立てる 第2部第4章の分散関係)
  2. (2) すべてのモードの Grüneisen 定数が等しく γ=3ba/2\gamma=3ba/2 になることを示せ。(使う道具:3.2.3 節)
  3. (3) 準調和近似で、線膨張率が αL(T)=α∞ C(T)/(NkB)\alpha_L(T)=\alpha_\infty\,C(T)/(N k_{\mathrm B})(α∞\alpha_\infty は古典の高温の値 (3.2.12))と書けることを示せ。低温(T≪ℏωm/kBT\ll\hbar\omega_{\mathrm m}/ k_{\mathrm B})で C=2π3NkB kBT/(ℏωm)C=\frac{2\pi}{3}N k_{\mathrm B}\, k_{\mathrm B} T/(\hbar\omega_{\mathrm m}) になることを使い、10 K10\,\mathrm K の αL\alpha_L を求めよ。3次元の結晶と温度依存の形がどう違うか述べよ。(使う道具:3.2.3 節、解析力学から様々な理論分野を組み立てる 第5部第3章の1次元の比熱)
  4. (4) この鎖で、3次の項による3フォノン過程は熱抵抗を生じさせるか。3次の項だけを摂動として残した Peierls の理論で、熱伝導率はどうなるか。(使う道具:3.2.4 節、3.2.5 節)
  5. (5) 高温の極限で、張力 0 の比熱 CfC_f と長さ一定の比熱 CLC_L の差を NkBN k_{\mathrm B} を単位に kBT/Ebk_{\mathrm B} T/E_{\mathrm b} で表し、T=300 KT=300\,\mathrm K を代入した値を求めよ(ここでは古典の式をそのまま使う)。(使う道具:3.2.2 節、3.2.5 節)
解答

(1) K=2Ebb2=2×0.4×1.602×10−19×2.25×1020 N m−1=28.8 N m−1K=2E_{\mathrm b}b^2=2\times0.4\times1.602\times10^{-19}\times2.25\times10^{20}\,\mathrm{N\,m^{-1}}=28.8\,\mathrm{N\,m^{-1}}。ωm=228.84/1.0×10−25 s−1=3.40×1013 s−1\omega_{\mathrm m}=2\sqrt{28.84/1.0\times10^{-25}}\,\mathrm{s^{-1}}=3.40\times10^{13}\,\mathrm{s^{-1}}、ℏωm/kB=259 K\hbar\omega_{\mathrm m}/ k_{\mathrm B}=259\,\mathrm K。室温はこの温度と同じ程度なので、古典の扱いは室温では目安にとどまる。

(2) 長さを L=NrL=Nr に変えると、分散関係は ωk=2ϕ′′(r)/M∣sin⁡(πν/N)∣\omega_k=2\sqrt{\phi''(r)/M}\lvert \sin(\pi\nu/N)\rvert(ka=2πν/Nka=2\pi\nu/N)で、モードの番号 ν\nu を固定すると sin⁡\sin の因子は rr によらない。したがって γk=−∂ln⁡ωk/∂ln⁡L=−(r/2)ϕ′′′(r)/ϕ′′(r)\gamma_k=-\partial\ln\omega_k/\partial\ln L=-(r/2)\phi'''(r)/\phi''(r) はすべてのモードで同じで、r=ar=a で −aK3/(2K)=6Ebb3a/(4Ebb2)=3ba/2=5.63-aK_3/(2K)=6E_{\mathrm b}b^3a/(4E_{\mathrm b}b^2)=3ba/2=5.63 である。

(3) 例題 3.2.3 (2) の αL=γC/(NKa2)\alpha_L=\gamma C/(NKa^2) で、γ\gamma が温度によらないので αL∝C(T)\alpha_L\propto C(T)。C=NkBC=N k_{\mathrm B} のとき α∞=4.31×10−5 K−1\alpha_\infty=4.31\times10^{-5}\,\mathrm{K^{-1}} なので、αL(T)=α∞C(T)/(NkB)\alpha_L(T)=\alpha_\infty C(T)/(N k_{\mathrm B}) である。10 K10\,\mathrm K では C/(NkB)=2π3×10/259.4=0.0807C/(N k_{\mathrm B})=\frac{2\pi}{3}\times10/259.4=0.0807(モードの和を数値でとると 0.08120.0812)で、αL=3.48×10−6 K−1\alpha_L=3.48\times10^{-6}\,\mathrm{K^{-1}}。3次元の結晶では低温の比熱が T3T^3 なので熱膨張率も T3T^3 で 0 に近づくが、1次元の鎖では低温の比熱が TT に比例するので、熱膨張率も TT に比例して 0 に近づく。どちらも絶対零度で熱膨張率が 0 になること(熱力学の第三法則と整合)は共通である。

(4) 単原子の鎖は音響の枝1本だけなので、例題 3.2.5 のとおり3フォノン過程はエネルギーと結晶運動量の保存を同時に満たせず、N 過程も U 過程も起きない。3次の項だけを摂動として残した Peierls の理論では衝突の項が 0 で、フォノンの緩和時間は無限大、熱伝導率は無限大のままである。フォノンが散乱されるようにするには、4次の項による4フォノン過程を入れる必要がある。ただし 3.2.5 節で見たとおり、1次元の非調和な鎖では4フォノン過程を入れても熱伝導率が長さとともに発散する例が知られていて、Fourier の法則の前提そのものが成り立たない。

(5) Cf/(NkB)=1+kBT/EbC_f/(N k_{\mathrm B})=1+ k_{\mathrm B} T/E_{\mathrm b}(式 (3.2.10))、CL/(NkB)=1−kBT/(8Eb)C_L/(N k_{\mathrm B})=1- k_{\mathrm B} T/(8E_{\mathrm b})(式 (3.2.25))なので、差は 98kBT/Eb\frac98 k_{\mathrm B} T/E_{\mathrm b}。これは TLα∞2B1/(NkB)TL\alpha_\infty^2B_1/(N k_{\mathrm B}) と一致する。300 K300\,\mathrm K では kBT/Eb=0.0646k_{\mathrm B} T/E_{\mathrm b}=0.0646 で、差は 0.07270.0727、つまり CfC_f と CLC_L は 7%7\% 程度違う。(1) で見たように 300 K300\,\mathrm K は ℏωm/kB\hbar\omega_{\mathrm m}/ k_{\mathrm B} と同じ程度なので、量子の効果まで入れると両方とも値は下がる。この数字は古典の式の中での大きさの目安である。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ