정준 앙상블

편집 역사 토론
계산물리 물리 마지막 수정: 2026-08-01 04:27:15

1. 개요[편집]

정준 앙상블
Canonical Ensemble
다른 이름NVT 앙상블
고정량입자수 N, 부피 V, 온도 T
확률 가중볼츠만 인자 exp(−E/kT)
특성 함수헬름홀츠 자유에너지 F = −kT ln Z
주 구현메트로폴리스 MC, MD 서모스탯

정준 앙상블(canonical ensemble)은 입자수 NN·부피 VV·온도 TT가 고정된 계가 가질 수 있는 미시상태들의 확률분포로, 에너지 EiE_i인 미시상태 ii가 나타날 확률이 볼츠만 인자에 비례하는(pi∝e−βEip_i \propto e^{-\beta E_i}, β=1/kBT\beta = 1/k_B T) 앙상블이다. 통계역학과 시뮬레이션에서 압도적으로 자주 쓰이는 기본 설정이라, 논문에서 그냥 “NVT로 돌렸다”고 하면 이 얘기다.

범위 정리: 이 문서는 앙상블 자체의 정의·분배함수·요동 관계·앙상블 등가성을 다룬다. 온도를 실제로 붙잡는 MD 알고리즘의 비교는 서모스탯, 속도 성분 하나의 분포만 따로 보는 것은 맥스웰-볼츠만 분포, 임계점 근방의 거동은 상전이 문서에 있다.

2. 열저장조에서 나오는 볼츠만 인자[편집]

정준 앙상블은 공리로 던져지는 게 아니라 미시정준 앙상블에서 유도된다. 관심 있는 작은 계 SS와 거대한 열저장조 RR를 붙여 놓고, 둘을 합친 전체는 에너지 EtotE_\text{tot}의 고립계라고 하자. 고립계에서는 모든 미시상태가 등확률이므로, SS가 상태 ii에 있을 확률은 저장조가 나머지 에너지를 나눠 가질 수 있는 경우의 수에 비례한다.

pi∝ΩR(Etot−Ei)=exp⁡ ⁣[SR(Etot−Ei)kB]p_i \propto \Omega_R(E_\text{tot} - E_i) = \exp\!\left[\frac{S_R(E_\text{tot} - E_i)}{k_B}\right]

저장조가 충분히 커서 Ei≪EtotE_i \ll E_\text{tot}이면 엔트로피를 1차까지 전개할 수 있고, ∂SR/∂E=1/T\partial S_R/\partial E = 1/T이므로

SR(Etot−Ei)≈SR(Etot)−EiTS_R(E_\text{tot} - E_i) \approx S_R(E_\text{tot}) - \frac{E_i}{T}

가 되어 pi∝e−βEip_i \propto e^{-\beta E_i}가 떨어진다. 전개를 1차에서 끊었다는 게 핵심 가정이고, 2차 항이 살아남는 유한 저장조에서는 정준분포가 아니라 그 보정이 붙은 분포가 나온다.1 흔한 오해와 달리 “볼츠만 인자는 낮은 에너지를 선호한다”는 서술은 반쪽짜리다 — 실제 확률은 p(E)∝g(E)e−βEp(E) \propto g(E) e^{-\beta E}로 상태밀도 g(E)g(E)와의 곱이며, 거시계에서 에너지 분포가 좁은 봉우리를 이루는 이유가 바로 급증하는 g(E)g(E)와 급감하는 e−βEe^{-\beta E}의 경합이다.

3. 분배함수와 자유에너지[편집]

규격화 상수가 곧 분배함수다.

Z(N,V,T)=∑ie−βEi,Z=1N! h3N∫e−βH(q,p) d3Nq d3NpZ(N,V,T) = \sum_i e^{-\beta E_i}, \qquad Z = \frac{1}{N!\,h^{3N}} \int e^{-\beta H(\mathbf{q},\mathbf{p})}\, d^{3N}q\, d^{3N}p

여기서 모든 열역학이 나온다. 헬름홀츠 자유에너지가 특성 함수이고, 나머지는 미분 한 번씩이다.

F=−kBTln⁡Z,⟨E⟩=−∂ln⁡Z∂β,S=−(∂F∂T)N,V,p=−(∂F∂V)N,TF = -k_B T \ln Z, \qquad \langle E \rangle = -\frac{\partial \ln Z}{\partial \beta}, \qquad S = -\left(\frac{\partial F}{\partial T}\right)_{N,V}, \qquad p = -\left(\frac{\partial F}{\partial V}\right)_{N,T}

고전계에서 운동량 적분은 가우스 적분이라 해석적으로 떨어지고 열적 드브로이 파장만 남긴다. 진짜 어려운 건 구성 적분 Zconf=∫e−βU(q)d3NqZ_\text{conf} = \int e^{-\beta U(\mathbf{q})} d^{3N}q이며, 이게 3N3N차원이라 손으로 안 풀리기 때문에 몬테카를로 방법과 분자동역학이 밥벌이를 한다.

실무에서 뼈아픈 사실 하나. 시뮬레이션은 ZZ 자체를 직접 측정하지 못한다. 평균값 ⟨A⟩=∑iAie−βEi/Z\langle A \rangle = \sum_i A_i e^{-\beta E_i} / Z에서는 ZZ가 분모·분자로 약분되지만, FF에는 ln⁡Z\ln Z가 통째로 들어가기 때문이다. 그래서 자유에너지는 항상 두 상태의 차이로 구하며, 자유에너지 섭동이나 열역학적 적분 같은 우회로가 따로 존재한다.

4. 에너지 요동과 열용량[편집]

ln⁡Z\ln Z를 β\beta로 두 번 미분하면 요동-소산 관계의 교과서적 사례가 나온다.

⟨ΔE2⟩=⟨E2⟩−⟨E⟩2=∂2ln⁡Z∂β2=kBT2CV\langle \Delta E^2 \rangle = \langle E^2 \rangle - \langle E \rangle^2 = \frac{\partial^2 \ln Z}{\partial \beta^2} = k_B T^2 C_V

정준 앙상블에서 에너지는 고정값이 아니라 요동친다는 점, 그리고 그 요동의 크기가 응답계수인 열용량과 같다는 점이 요지다. EE와 CVC_V 모두 NN에 비례하므로 상대 요동은 ΔE/⟨E⟩∼N−1/2\Delta E/\langle E \rangle \sim N^{-1/2}이고, 아보가드로 수 규모에서는 10−1210^{-12} 수준이라 눈에 안 보인다. 반대로 입자 수백 개짜리 MD 상자에서는 총에너지가 몇 퍼센트씩 흔들리는 게 정상이며, 이걸 버그로 오해하고 서모스탯을 조여 버리면 통계가 망가진다.

이 식은 실전에서 매우 유용하다. 열용량을 구하는 방법이 두 가지가 되기 때문이다. (a) 여러 온도에서 돌려 ⟨E⟩(T)\langle E \rangle (T)를 수치미분하거나, (b) 한 온도에서 에너지 요동 분산만 재거나. (b)는 한 번의 실행으로 끝나지만 분산 추정이라 수렴이 느리고, 상전이 근방에서는 요동 자체가 발산해 유효표본크기가 급락한다.

5. 다른 앙상블과의 비교[편집]

앙상블고정량요동하는 양특성 함수
미시정준 (NVE)N,V,EN, V, E없음 (온도가 유도량)엔트로피 S=kBln⁡ΩS = k_B \ln \Omega
정준 (NVT)N,V,TN, V, T에너지헬름홀츠 F=−kBTln⁡ZF = -k_B T \ln Z
대정준 (μVT)μ,V,T\mu, V, T에너지, 입자수대퍼텐셜 Ω=−kBTln⁡Ξ\Omega = -k_B T \ln \Xi
등온등압 (NPT)N,p,TN, p, T에너지, 부피깁스 GG

열역학적 극한(N→∞N \to \infty)에서는 이들이 등가다 — 상대 요동이 0으로 가므로 어느 앙상블로 계산하든 같은 상태방정식이 나온다. 다만 유한계에서는 다르며, 1차 상전이의 공존 영역에서 미시정준 앙상블은 음의 비열을 가질 수 있지만 정준 앙상블은 ⟨ΔE2⟩≥0\langle \Delta E^2 \rangle \ge 0 때문에 원리적으로 불가능하다.2 화학 퍼텐셜을 고정하는 대정준 앙상블은 흡착·다공성 매질 충전 같은 문제에서 필수다.

6. 시뮬레이션에서의 구현[편집]

몬테카를로 쪽은 깔끔하다. 메트로폴리스 수용확률 min⁡(1,e−βΔE)\min(1, e^{-\beta \Delta E})는 정준분포를 정확히 만족하는 상세균형 조건에서 유도되므로, 마르코프 연쇄 몬테카를로로 이징 모형을 돌리면 별도 장치 없이 NVT다. 표본이 상관되는 문제는 병렬 템퍼링이나 해밀토니안 몬테카를로로 완화한다.

MD 쪽은 사정이 다르다. 뉴턴 방정식은 그냥 두면 NVE라서 온도를 붙잡는 장치를 덧대야 하고, 그 장치가 정말 정준분포를 만드는지는 알고리즘마다 다르다.

  • 노제-후버: 확장 라그랑지안에 열욕 자유도를 추가한 결정론적 방식. 정준분포를 생성한다는 증명이 있지만 에르고딕 가설에 의존해, 조화진동자처럼 단순한 계에서는 실제로 깨진다. 그래서 열욕을 사슬로 잇는 노제-후버 체인이 표준이 됐다.
  • 랑주뱅: 마찰항과 잡음항을 넣어 요동-소산 정리를 직접 강제한다. 에르고딕성이 강건해 요즘 가장 무난한 선택. 상세는 랑주뱅 동역학 참고.
  • 베렌트센: 속도를 지수적으로 목표 온도로 끌어당긴다. 평균 온도는 맞지만 운동에너지 분포의 폭이 정준값보다 좁아 정준 앙상블이 아니다. 평형화 단계에서만 쓰고 생산 구간에서는 갈아타는 것이 국룰.3

7. 관련 문서[편집]

8. Footnotes[편집]

  1. 저장조의 열용량이 유한하면 2차 항 ∝Ei2/(kBT2CR)\propto E_i^2 / (k_B T^2 C_R)가 살아남는다. 이게 이른바 “초정준(supercanonical)” 보정이고, 클러스터 물리처럼 저장조가 작을 때 실제로 문제가 된다. 대부분의 교과서는 “저장조는 무한히 크다”는 한 줄로 이 항을 조용히 묻는다. ↩

  2. 이 때문에 유한 클러스터의 녹는 전이에서 미시정준 열용량 곡선이 음수 구간을 갖는 실험 결과가 나왔을 때 한바탕 소동이 있었다. 앙상블 등가성이 성립하지 않는 영역이 있다는 것 자체는 이미 알려져 있었지만, 실물로 보는 건 또 다른 얘기였다. ↩

  3. 베렌트센 서모스탯이 “정준이 아니다”라는 사실은 알려진 지 오래인데도 논문 메서드 절에서 여전히 자주 목격된다. 심지어 결합상수 τ\tau를 작게 줄일수록 분포가 더 나빠지는데, 사람들은 “온도가 빨리 맞으니 좋은 것”으로 읽는다. 온도가 맞는 것과 분포가 맞는 것은 다른 얘기다. ↩