第1章：フォノンとは何か？

格子振動の量子化と準粒子概念の導入

📖 読了時間: 約60分 📊 難易度: 入門レベル 💻 コード例: 0個 📝 演習問題: 0問

🌐 EN | 🇯🇵 JP | Last sync: 2025-12-20

材料科学道場 > フォノン入門 > 第1章

第1章：フォノンとは何か？

格子振動の量子化と準粒子概念の導入

学習目標

1.1 はじめに：なぜフォノンを学ぶのか？

固体材料の多くの性質は、原子がどのように振動しているかによって決まります。熱伝導率、比熱、熱膨張、さらには超伝導や相転移といった現象まで、すべて格子振動と深く関わっています。**フォノン（phonon）**は、この格子振動を量子力学的に記述するための基本概念です。

本章では、古典的な格子振動の描像から始めて、それを量子化することでフォノンという準粒子の概念がどのように生まれるのかを学びます。この理解は、材料科学や固体物理学における多くの現象を理解するための基礎となります。

1.2 格子振動の古典的描像

1.2.1 バネで結ばれた原子モデル

結晶固体では、原子は規則正しく配列しています。しかし、原子は静止しているわけではなく、平衡位置の周りで常に振動しています。この振動を理解するために、最も単純なモデルとして一次元原子鎖を考えましょう。

一次元原子鎖モデル

質量 \(m\) の同一な原子が、格子定数 \(a\) の間隔で一列に並んでおり、隣接する原子同士がバネ定数 \(K\) のバネで結ばれていると考えます。各原子は平衡位置からのずれ \(u_n\) だけ変位できます。

この単純なモデルは、実際の結晶における原子間相互作用を表現しています。原子間の化学結合が「バネ」に相当し、平衡位置からずれると復元力が働きます。

1.2.2 調和近似

原子が平衡位置から変位したとき、その原子が感じるポテンシャルエネルギー \(V(u)\) は、一般には複雑な関数です。しかし、変位が小さい場合、ポテンシャルを平衡位置の周りでTaylor展開できます：

調和近似（Harmonic Approximation）

平衡位置 \(u = 0\) の周りでポテンシャルエネルギーを展開すると、

\[ V(u) = V(0) + \left(\frac{dV}{du}\right){u=0} u + \frac{1}{2}\left(\frac{d^2V}{du^2}\right){u=0} u^2 + \cdots \]

平衡位置では \(\frac{dV}{du} = 0\) であり、\(V(0)\) は定数として無視できるため、最低次の項は二次の項になります。この近似を調和近似といい、ポテンシャルは

\[ V(u) \approx \frac{1}{2}Ku^2 \]

となります。ここで \(K = \frac{d^2V}{du^2}\) はバネ定数に対応します。

調和近似の妥当性

調和近似は以下の条件で良い近似になります：

室温以下の多くの結晶材料では、調和近似は驚くほどよく成立します。ただし、高温や相転移近傍、非線形光学現象などでは、非調和項（三次以上の項）を考慮する必要があります。

1.2.3 運動方程式と基準振動

\(n\) 番目の原子の運動方程式は、Newton の第二法則から次のように書けます：

\[ m\frac{d^2u_n}{dt^2} = K(u_{n+1} - u_n) + K(u_{n-1} - u_n) = K(u_{n+1} + u_{n-1} - 2u_n) \]

この連立微分方程式を解くために、基準振動（normal mode）の解を仮定します。すべての原子が同じ角振動数 \(\omega\) で振動する波動解を考えると：

\[ u_n(t) = A e^{i(kna - \omega t)} \]

ここで、\(k\) は波数、\(A\) は振幅、\(na\) は \(n\) 番目の原子の平衡位置です。この解を運動方程式に代入すると、

\[ -m\omega^2 = K(e^{ika} + e^{-ika} - 2) = 2K(\cos(ka) - 1) \]

整理すると、分散関係（dispersion relation）が得られます：

一次元単原子格子の分散関係

\[ \omega(k) = 2\sqrt{\frac{K}{m}} \left|\sin\left(\frac{ka}{2}\right)\right| \]

この式は、波数 \(k\) と角振動数 \(\omega\) の関係を与えます。最大振動数は \(\omega_{\text{max}} = 2\sqrt{K/m}\) です。

例：銅の格子振動

銅（Cu）の場合、実験的に測定された最大フォノン振動数は約 \(\omega_{\text{max}} \approx 5 \times 10^{13}\) rad/s です。これは、振動周期が約 \(10^{-13}\) 秒、つまり0.1ピコ秒程度に対応します。このような高速な振動が、材料の熱的性質を決定しています。

1.3 格子振動の量子化：フォノンの導入

1.3.1 量子力学的な調和振動子

古典力学では、格子振動のエネルギーは連続的な値をとることができます。しかし、量子力学では状況が異なります。ミクロな系では、エネルギーは離散的な値（量子化された値）しかとることができません。

一次元調和振動子のHamiltonian（エネルギー演算子）は：

\[ \hat{H} = \frac{\hat{p}^2}{2m} + \frac{1}{2}m\omega^2\hat{x}^2 \]

ここで、\(\hat{p}\) は運動量演算子、\(\hat{x}\) は位置演算子です。このHamiltonianの固有値（許されるエネルギー準位）は、量子力学の標準的な計算により次のようになります：

調和振動子のエネルギー準位

\[ E_n = \left(n + \frac{1}{2}\right)\hbar\omega, \quad n = 0, 1, 2, 3, \ldots \]

ここで、\(n\) は量子数、\(\hbar\) はプランク定数を \(2\pi\) で割った換算プランク定数（\(\hbar = 1.055 \times 10^{-34}\) J·s）です。

この式から、重要な点が2つ読み取れます：

1.3.2 フォノンの概念

格子振動を量子化すると、振動のエネルギーが \(\hbar\omega\) の単位でやり取りされることがわかりました。この状況を、あたかもエネルギー \(\hbar\omega\) を持つ粒子が \(n\) 個存在すると解釈することができます。

フォノン（Phonon）

フォノンとは、格子振動の量子化によって生じる準粒子（quasiparticle）です。波数 \(k\) のフォノンは、以下の性質を持ちます：

基準振動モード \((k, \omega)\) に \(n\) 個のフォノンが励起されていることは、そのモードのエネルギーが \(E_n = (n + 1/2)\hbar\omega\) であることと等価です。

「準粒子」とは？

フォノンは電子や光子のような「真の粒子」ではなく、準粒子（quasiparticle）と呼ばれます。これは、多数の原子の集団的な運動を、あたかも一つの粒子であるかのように扱うための便利な概念です。準粒子の考え方は、複雑な多体系を理解するための強力な道具となっています。

準粒子の例：

1.3.3 ボソン的性質

フォノンはボソン（boson）です。これは、次のような性質を意味します：

\[ \langle n \rangle = \frac{1}{e^{\hbar\omega/k_B T} - 1} \]

ここで \(k_B\) はBoltzmann定数です。

例：室温でのフォノン占有数

室温（\(T = 300\) K）において、\(k_B T \approx 26\) meV です。振動数 \(\omega = 2 \times 10^{13}\) rad/s のフォノンは \(\hbar\omega \approx 13\) meV なので、

\[ \langle n \rangle = \frac{1}{e^{13/26} - 1} \approx \frac{1}{1.65 - 1} = 1.54 \]

つまり、室温では平均して約1.5個のフォノンがそのモードに励起されています。これが材料の熱容量や熱伝導に寄与します。

1.4 ゼロ点エネルギーとその物理的意義

1.4.1 ゼロ点エネルギー

調和振動子のエネルギー \(E_n = (n + 1/2)\hbar\omega\) において、\(n = 0\) の基底状態でも \(E_0 = \frac{1}{2}\hbar\omega\) というエネルギーが残ります。これをゼロ点エネルギー（zero-point energy）といいます。

ゼロ点エネルギー

絶対零度（\(T = 0\) K）においても、量子力学的な不確定性原理により、原子は完全に静止することができません。基底状態でのエネルギー \(E_0 = \frac{1}{2}\hbar\omega\) がゼロ点エネルギーです。

結晶全体のゼロ点エネルギーは、すべてのモード \(\mathbf{k}\) について和をとって：

\[ E_{\text{ZP}} = \sum_{\mathbf{k}} \frac{1}{2}\hbar\omega(\mathbf{k}) \]

1.4.2 ゼロ点エネルギーの重要性

ゼロ点エネルギーは、単なる理論的な興味を超えて、実際の物理現象に重要な影響を与えます：

例：水素の同位体効果

水素（H₂）と重水素（D₂）の結合解離エネルギーは、質量の違いによるゼロ点エネルギーの差により異なります。H₂の方が軽いため、ゼロ点エネルギーが大きく、結合エネルギーは約5%低くなります。この効果は化学反応速度にも影響を与え、kinetic isotope effectとして知られています。

1.5 光子との類比

フォノンの概念は、電磁波の量子化である光子（photon）との類比で理解するとより明確になります。

性質光子（Photon）フォノン（Phonon）
元の波電磁波格子振動（弾性波）
媒質真空中でも伝播結晶格子が必要
エネルギー\(E = \hbar\omega\)\(E = \hbar\omega\)
運動量\(\mathbf{p} = \hbar\mathbf{k}\)\(\mathbf{p} = \hbar\mathbf{k}\)（準運動量）
スピン1（ベクトルボソン）0（スカラーボソン）
統計性Bose-Einstein統計Bose-Einstein統計
分散関係\(\omega = c\mathbf{k}
静止質量00（準粒子として）

重要な違い：準運動量

光子の運動量 \(\mathbf{p} = \hbar\mathbf{k}\) は真の運動量で、並進対称性に対応する保存量です。一方、フォノンの \(\hbar\mathbf{k}\) は準運動量または結晶運動量と呼ばれ、離散的な並進対称性に対応します。そのため、フォノンの運動量保存則は逆格子ベクトル \(\mathbf{G}\) を除いて保存されるという修正を受けます：

\[ \mathbf{k}_1 + \mathbf{k}_2 = \mathbf{k}_3 + \mathbf{G} \]

これは、フォノンの相互作用やフォノン散乱過程を理解する上で重要です。

1.6 Python による格子振動の可視化

理論を具体的に理解するために、Pythonを使って一次元原子鎖の振動を可視化してみましょう。

一次元原子鎖の波動解の可視化

import numpy as np
import matplotlib.pyplot as plt
from matplotlib.animation import FuncAnimation

# パラメータ設定
N = 20           # 原子数
a = 1.0          # 格子定数
k = 2*np.pi/10   # 波数
omega = 1.0      # 角振動数
A = 0.3          # 振幅

# 原子の平衡位置
n = np.arange(N)
x_eq = n * a

# 時間発展
fig, ax = plt.subplots(figsize=(12, 4))

def update(frame):
    ax.clear()
    t = frame * 0.1

    # 各原子の変位
    u = A * np.sin(k * x_eq - omega * t)
    x_actual = x_eq + u

    # プロット
    ax.plot(x_eq, np.zeros_like(x_eq), 'o',
            color='lightgray', markersize=8, label='平衡位置')
    ax.plot(x_actual, np.zeros_like(x_actual), 'o',
            color='blue', markersize=12, label='実際の位置')

    # 変位を矢印で表示
    for i in range(N):
        if abs(u[i]) > 0.01:
            ax.arrow(x_eq[i], 0, u[i], 0, head_width=0.1,
                    head_length=0.05, fc='red', ec='red', alpha=0.6)

    ax.set_xlim(-1, N*a)
    ax.set_ylim(-0.5, 0.5)
    ax.set_xlabel('位置', fontsize=12)
    ax.set_title(f'一次元原子鎖の振動（k={k:.2f}, t={t:.2f}）', fontsize=14)
    ax.legend()
    ax.grid(True, alpha=0.3)

anim = FuncAnimation(fig, update, frames=100, interval=50)
plt.tight_layout()
plt.show()

分散関係のプロット

import numpy as np
import matplotlib.pyplot as plt

# パラメータ
K = 1.0  # バネ定数
m = 1.0  # 質量
a = 1.0  # 格子定数

# 波数の範囲（第一Brillouin帯）
k = np.linspace(-np.pi/a, np.pi/a, 500)

# 分散関係
omega = 2 * np.sqrt(K/m) * np.abs(np.sin(k*a/2))

# プロット
plt.figure(figsize=(10, 6))
plt.plot(k, omega, linewidth=2, color='blue')
plt.axhline(y=2*np.sqrt(K/m), color='red', linestyle='--',
            label=r'$\omega_{\max} = 2\sqrt{K/m}$')
plt.axvline(x=0, color='gray', linestyle='-', alpha=0.3)
plt.axvline(x=-np.pi/a, color='gray', linestyle='--', alpha=0.3)
plt.axvline(x=np.pi/a, color='gray', linestyle='--', alpha=0.3)

plt.xlabel(r'波数 $k$ [1/a]', fontsize=14)
plt.ylabel(r'角振動数 $\omega$ [rad/s]', fontsize=14)
plt.title('一次元単原子格子の分散関係', fontsize=16)
plt.legend(fontsize=12)
plt.grid(True, alpha=0.3)
plt.xlim(-np.pi/a, np.pi/a)
plt.ylim(0, 2.2*np.sqrt(K/m))
plt.tight_layout()
plt.show()

# 群速度の計算とプロット
dk = k[1] - k[0]
v_g = np.gradient(omega, dk)

plt.figure(figsize=(10, 6))
plt.plot(k, v_g, linewidth=2, color='green')
plt.axhline(y=0, color='black', linestyle='-', linewidth=0.5)
plt.xlabel(r'波数 $k$ [1/a]', fontsize=14)
plt.ylabel(r'群速度 $v_g = d\omega/dk$ [m/s]', fontsize=14)
plt.title('群速度（フォノンの伝播速度）', fontsize=16)
plt.grid(True, alpha=0.3)
plt.tight_layout()
plt.show()

コードの解説

1.7 歴史的背景

1.7.1 Einstein の貢献

Albert Einstein（1907年）は、固体の比熱を説明するために、すべての原子が同じ振動数 \(\omega_E\) で独立に振動するEinsteinモデルを提唱しました。これは量子論を固体物理学に応用した最初の成功例の一つです。

Einstein モデルの比熱

Einstein モデルでは、モル比熱は次式で与えられます：

\[ C_V = 3R \left(\frac{\Theta_E}{T}\right)^2 \frac{e^{\Theta_E/T}}{(e^{\Theta_E/T} - 1)^2} \]

ここで \(\Theta_E = \hbar\omega_E/k_B\) はEinstein温度、\(R\) は気体定数です。この式は高温で古典的なDulong-Petit則（\(C_V = 3R\)）に一致し、低温で指数関数的にゼロに近づきます。

1.7.2 Debye の改良

Peter Debye（1912年）は、原子が独立に振動するのではなく、連続的な弾性波として振る舞うというより現実的なモデルを提案しました。Debyeモデルでは、振動数に上限（Debye振動数 \(\omega_D\)）を設け、それ以下の振動数を持つモードが連続的に分布すると仮定します。

Debye モデルの低温比熱

低温（\(T \ll \Theta_D\)、\(\Theta_D = \hbar\omega_D/k_B\) はDebye温度）では、比熱は温度の3乗に比例します：

\[ C_V = \frac{12\pi^4}{5}R\left(\frac{T}{\Theta_D}\right)^3 \]

このDebye \(T^3\) 則は、多くの結晶材料で実験的に確認されており、現代的なフォノン理論の基礎となっています。

1.7.3 「フォノン」という名前の由来

「phonon」という用語は、ロシアの物理学者Igor Tamm（1932年）によって初めて使われました。ギリシャ語の「phone」（音）に由来し、音波・弾性波の量子化を意味します。「光子（photon）」との類推で命名されました。

1.8 まとめ

本章の要点

演習問題

問題1：基礎概念の確認

以下の空欄を埋めなさい：

  1. 調和近似では、ポテンシャルエネルギーを \(V(u) = \frac{1}{2}K\)_______ と近似する。
  2. 調和振動子のエネルギー準位は \(E_n = (n + \)_______\()\hbar\omega\) である。
  3. フォノンは_______統計に従うボソンである。

問題2：数値計算

銅（Cu）のDebye温度は \(\Theta_D = 343\) K です。

  1. Debye振動数 \(\omega_D\) を計算しなさい。（\(\hbar = 1.055 \times 10^{-34}\) J·s、\(k_B = 1.381 \times 10^{-23}\) J/K）
  2. 室温（300 K）における平均フォノン数 \(\langle n \rangle\) を、\(\omega = \omega_D\) として計算しなさい。

問題3：応用問題

一次元原子鎖の分散関係 \(\omega(k) = 2\sqrt{K/m}|\sin(ka/2)|\) において、

  1. 群速度 \(v_g = d\omega/dk\) を計算しなさい。
  2. \(k = 0\) 付近での群速度（音速）を求めなさい。
  3. Brillouin帯端（\(k = \pm\pi/a\)）での群速度を求め、その物理的意味を説明しなさい。

問題4：概念理解

以下の問いに答えなさい：

  1. フォノンが「準粒子」と呼ばれる理由を、光子との違いを踏まえて説明しなさい。
  2. ゼロ点エネルギーが物質の性質にどのような影響を与えるか、具体例を挙げて説明しなさい。
  3. 調和近似が破綻するのはどのような状況か、3つ例を挙げなさい。

← シリーズトップへ | 次章：フォノン分散関係 →

免責事項

この教育コンテンツは、橋本研究室のナレッジベース用にAIの支援を受けて作成されました。内容の正確性には十分注意を払っていますが、学術的な引用や重要な意思決定には、査読付き論文や標準的な教科書を参照することをお勧めします。

シリーズ目次に戻る 次の章 →