디바이 모형 (Debye Model)

1. 고전적 비열의 한계

고전 통계역학에 의하면, $N$ 개의 원자로 이루어진 고체의 격자 비열은 온도에 무관한 상수이다.

법칙2.1뒤롱-프티 법칙

고전적 등분배 정리에 의하면, 3차원 고체의 몰비열은:

C_V = 3Nk_B = 3R \approx 24.9\,\text{J/(mol·K)}

이를 뒤롱-프티 법칙(Dulong-Petit law)이라 한다. 이는 고온에서 잘 성립하지만, 저온에서 비열이 0으로 감소하는 실험적 사실을 설명하지 못한다.

2. 아인슈타인 모형

아인슈타인(1907)은 최초로 격자 비열에 양자론을 적용하였다. 모든 원자가 동일한 진동수 $\omega_E$ 로 독립적으로 진동한다고 가정한다.

아인슈타인 모형의 비열:

C_V^{\text{Ein}} = 3Nk_B \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$ 는 아인슈타인 온도이다. 이 모형은 저온에서 비열의 감소를 정성적으로 설명하지만, $T \to 0$ 에서 $C_V \propto e^{-\Theta_E/T}$ 로 지수적으로 감소하여 실험( $C_V \propto T^3$ )과 맞지 않는다.

3. 디바이 모형

법칙2.2디바이 모형

디바이(1912)는 격자 진동을 연속 탄성체의 음파로 근사하되, 총 모드 수를 $3N$ 으로 제한하여 비열을 계산하였다. 핵심 가정:

분산 관계를 선형 $\omega = v_s q$ 로 근사 (등방적 디바이 모형에서 종파·횡파 속도의 적절한 평균 사용)
디바이 차단 진동수 $\omega_D$ 를 도입하여 총 모드 수를 $3N$ 으로 제한:

\int_0^{\omega_D} g(\omega)\,d\omega = 3N

디바이 상태 밀도:

g(\omega) = \frac{9N}{\omega_D^3}\omega^2, \quad 0 \leq \omega \leq \omega_D

디바이 온도:

\Theta_D = \frac{\hbar\omega_D}{k_B}

4. 디바이 비열의 유도

유도디바이 비열 공식

격자 진동의 내부 에너지는:

U = \int_0^{\omega_D} \hbar\omega\,\langle n(\omega)\rangle\,g(\omega)\,d\omega = \int_0^{\omega_D} \frac{\hbar\omega}{e^{\hbar\omega/k_BT}-1}\cdot\frac{9N\omega^2}{\omega_D^3}\,d\omega

$x = \hbar\omega/k_BT$ 로 치환하면, $x_D = \Theta_D/T$ 이고:

U = 9Nk_BT\left(\frac{T}{\Theta_D}\right)^3\int_0^{x_D}\frac{x^3}{e^x-1}\,dx

비열 $C_V = \partial U/\partial T$ :

\boxed{C_V = 9Nk_B\left(\frac{T}{\Theta_D}\right)^3\int_0^{\Theta_D/T}\frac{x^4 e^x}{(e^x-1)^2}\,dx}

이를 디바이 비열 공식이라 한다. 디바이 함수 $D(x_D)$ 를 정의하면:

C_V = 3Nk_B \cdot D\!\left(\frac{\Theta_D}{T}\right)

■

5. 극한 거동

고온 극한 ( $T \gg \Theta_D$ , $x_D \to 0$ ):

적분 변수가 모두 작아 $x^4e^x/(e^x-1)^2 \approx x^2$ 을 사용하면:

C_V \to 9Nk_B\left(\frac{T}{\Theta_D}\right)^3 \cdot \frac{1}{3}\left(\frac{\Theta_D}{T}\right)^3 = 3Nk_B

뒤롱-프티 법칙을 자연스럽게 회복한다.

저온 극한 ( $T \ll \Theta_D$ , $x_D \to \infty$ ):

유도디바이 $T^3$ 법칙

$T \ll \Theta_D$ 이면 $x_D = \Theta_D/T \to \infty$ 이므로 적분의 상한을 무한대로 보낼 수 있다:

\int_0^{\infty}\frac{x^3}{e^x-1}\,dx = \frac{\pi^4}{15}

따라서 내부 에너지는:

U = 9Nk_BT\left(\frac{T}{\Theta_D}\right)^3 \cdot \frac{\pi^4}{15} = \frac{3\pi^4 Nk_B}{5}\frac{T^4}{\Theta_D^3}

비열은:

\boxed{C_V = \frac{12\pi^4}{5}Nk_B\left(\frac{T}{\Theta_D}\right)^3 \approx 234\,Nk_B\left(\frac{T}{\Theta_D}\right)^3}

이것이 유명한 디바이 $T^3$ 법칙이다.

■

이 $T^3$ 의존성은 저온에서의 실험 결과와 매우 잘 일치하며, 아인슈타인 모형의 지수적 감소와 근본적으로 다르다. 그 물리적 원인은 저온에서 장파장 음향 포논만이 열적으로 여기되기 때문이다.

6. 디바이 온도의 물리적 의미

정의2.1디바이 온도의 의미

디바이 온도 $\Theta_D$ 는 격자 진동의 양자 효과가 중요해지는 특성 온도이다:

$T \gg \Theta_D$ : 모든 포논 모드가 고전적으로 여기됨 (고전 극한)
$T \ll \Theta_D$ : 고에너지 포논 모드가 "동결" (양자 극한)
$T \sim \Theta_D$ : 양자-고전 전이 영역

주요 물질의 디바이 온도:

| 물질 | $\Theta_D$ (K) | 물질 | $\Theta_D$ (K) | |------|-------|------|-------| | Pb | 105 | Cu | 343 | | Au | 165 | Si | 645 | | Ag | 225 | C (다이아몬드) | 2230 | | Al | 428 | Fe | 470 |

디바이 온도는 원자 결합의 강도와 원자 질량에 의존한다. 가벼운 원자와 강한 결합은 높은 $\Theta_D$ 를 준다. 다이아몬드의 극히 높은 디바이 온도는 매우 강한 $sp^3$ 공유 결합과 가벼운 탄소 원자의 결합 때문이다.

예제구리의 저온 비열

$\Theta_D = 343\,\text{K}$ 인 구리에 대해, $T = 10\,\text{K}$ 에서의 격자 비열을 디바이 $T^3$ 법칙으로 계산하면:

C_V = \frac{12\pi^4}{5}Nk_B\left(\frac{T}{\Theta_D}\right)^3 = \frac{12\pi^4}{5} \times 8.314 \times \left(\frac{10}{343}\right)^3

= 233.8 \times 8.314 \times 2.48 \times 10^{-5} = 0.0482\,\text{J/(mol·K)}

실제 실험값은 약 $0.05\,\text{J/(mol·K)}$ 로 잘 일치한다. 그러나 금속에서는 전자 비열 기여 $C_{\text{el}} = \gamma T$ 도 있어, 전체 저온 비열은:

C = \gamma T + \beta T^3

$C/T$ 대 $T^2$ 그래프의 기울기로부터 $\beta$ (따라서 $\Theta_D$ )를, 절편으로부터 $\gamma$ 를 결정할 수 있다.

참고디바이 모형의 한계

디바이 모형은 분산 관계를 선형으로 근사하고 등방성을 가정하므로, 광학 포논이 존재하는 다원자 결정이나 분산 관계의 비선형성이 중요한 영역에서는 정량적으로 부정확하다. 보다 정확한 비열 계산을 위해서는 실험적 또는 제일원리(ab initio) 계산으로 구한 정확한 포논 상태 밀도 $g(\omega)$ 를 사용해야 한다.