미시정준·정준·대정준 앙상블과 분배함수에서 열역학 상태량과 요동을 체계적으로 유도합니다.

CHEMISTRY · PART 9 · 대학원 브리지: 양자·통계·계면03 / 07

실험은 온도와 압력을 제어하지만 시뮬레이션은 에너지·부피·입자수를 고정하거나 thermostat·barostat로 분포를 샘플링합니다.

어떤 변수를 고정한 확률분포를 선택해야 실제 실험 조건과 올바른 평균을 비교할 수 있을까요?

이번 장의 답고립계는 NVE, 열저장고와 접촉한 계는 NVT, 입자 저장고와 교환하는 흡착계는 μVT가 자연스럽습니다. 분배함수의 로그 미분으로 평균에너지·압력·입자수와 요동을 얻습니다.

01기초부터 계산까지

CHEMISTRY, FROM ZERO

앙상블·분배함수: 처음 배우는 사람을 위한 출발점

분배함수는 상태를 단순히 세는 합이 아니라 각 에너지 상태의 Boltzmann 가중치를 모아 열역학 전체를 생성하는 함수입니다. 어떤 양이 고정되는지에 따라 미시정준·정준·대정준 앙상블을 고르고, 평균뿐 아니라 요동과 유한크기 효과를 함께 해석합니다.

START HERE

앙상블

같은 거시 제약을 만족하는 미시상태들의 확률적 모음입니다.

START HERE

분배함수

상태의 Boltzmann 가중치를 합해 열역학량을 생성하는 함수입니다.

START HERE

대정준 앙상블

온도·부피·화학퍼텐셜을 고정하고 입자수 변동을 허용하는 앙상블입니다.

CONCEPT ROUTE

이 글에서 생각이 이어지는 순서

용어를 외우는 대신, 관찰에서 모형과 계산을 거쳐 실제 해석으로 가는 흐름을 먼저 잡습니다.

01실험 제약

출발점

02앙상블

핵심 변수

03분배함수

모델 적용

04평균·요동

검증과 해석

$$ constraints → distribution → thermodynamics $$

READ THE EQUATION

공식을 계산기에 넣기 전에 읽는 법

$$ Z=\sum_i e^{-\beta E_i}, \qquad A=-k_BT\ln Z $$
β=1/kBT이며 상태합에는 에너지 준위의 축퇴와 고전 연속상태의 위상공간 측도가 포함됩니다.
01

정준앙상블에서 β=1/kBT이고 Z=Σe⁻βEi로 각 미시상태의 확률 pi=e⁻βEi/Z를 정규화합니다.

02

Helmholtz 자유에너지 A=−kBTlnZ, 평균에너지 U=−∂lnZ/∂β, 열용량은 에너지 요동과 연결됩니다.

03

입자 교환이 있으면 대정준 분배함수 Ξ와 화학퍼텐셜 μ를 사용하며 ensemble equivalence는 열역학한계 조건을 확인합니다.

WORKED FOUNDATIONS

기본 계산과 판단을 세 단계로 완전히 전개합니다

숫자 대입, 단위, 결과의 물리적 뜻을 한 번에 확인하세요.

01

에너지 0과 ε=kBT인 두 준위의 Z

$$ Z=1+e^{-\beta\varepsilon}=1+e^{-1}=1.368 $$

바닥상태 가중치는 1, 들뜬상태는 e⁻¹입니다. 상태 축퇴가 있다면 각 항에 degeneracy를 곱합니다.

02

두 준위의 들뜬상태 확률과 평균에너지

$$ p_1=\frac{e^{-1}}{1+e^{-1}}=0.269,\qquad U=p_1\varepsilon=0.269\,\varepsilon $$

온도가 낮아 ε≫kBT가 되면 들뜬상태 확률이 지수적으로 작아지고, 고온에서는 두 비축퇴 상태가 각각 1/2에 접근합니다.

03

두 준위계의 Helmholtz 자유에너지

$$ A=-k_BT\ln Z=-k_BT\ln(1.368)=-0.313\,k_BT $$

A가 바닥에너지 0보다 낮아지는 것은 상태 선택의 엔트로피 기여 때문입니다. 에너지 평균 U와 자유에너지 A를 혼동하지 않습니다.

스스로 확인하기 · 어떤 변수를 고정한 확률분포를 선택해야 실제 실험 조건과 올바른 평균을 비교할 수 있을까요?

고립계는 NVE, 열저장고와 접촉한 계는 NVT, 입자 저장고와 교환하는 흡착계는 μVT가 자연스럽습니다. 분배함수의 로그 미분으로 평균에너지·압력·입자수와 요동을 얻습니다.

02왜 이 개념이 필요한가

실험은 온도와 압력을 제어하지만 시뮬레이션은 에너지·부피·입자수를 고정하거나 thermostat·barostat로 분포를 샘플링합니다. 이 장면을 풀려면 용어를 외우는 것보다 어떤 물리량을 고정하고 무엇을 변화시키는지 먼저 정해야 합니다. 미시정준·정준·대정준 앙상블과 분배함수에서 열역학 상태량과 요동을 체계적으로 유도합니다.

읽기 전에 구분할 것
  • 앙상블은 시뮬레이션 기법 이름보다 제약조건과 확률분포입니다.
  • 시간평균과 앙상블평균의 동일성은 ergodicity·충분한 샘플링을 요구합니다.
  • 평균값만 맞아도 분포와 요동이 틀릴 수 있습니다.

03핵심 모델과 식

정준분포에서 상태 i의 확률은 e−βEi/Z입니다. Z는 단순 정규화 상수가 아니라 Helmholtz 에너지 A=−kBT lnZ를 생성합니다. 대정준분배함수 Ξ는 입자수 변동을 허용해 흡착·전자계와 연결됩니다.

$$ Z=\sum_i e^{-\beta E_i}, \qquad A=-k_BT\ln Z $$
β=1/kBT이며 상태합에는 에너지 준위의 축퇴와 고전 연속상태의 위상공간 측도가 포함됩니다.

04공학 예제로 계산하기

MOF 기체흡착 GCMC는 μVT에서 분자 삽입·삭제·이동을 샘플링합니다. 실험 압력을 μ로 바꾸려면 실제 기체 fugacity와 힘장 일관성이 필요합니다.

풀이와 해석의 순서
  • 실험에서 고정된 세기변수와 교환 가능한 보존량을 확인합니다.
  • 그 조건에 맞는 앙상블과 열역학 퍼텐셜을 선택합니다.
  • 평형화·상관시간·독립표본 수를 평가하고 평균과 신뢰구간을 보고합니다.

05오해·한계·다음 연결

유한계·장거리상호작용·1차 상전이에서는 앙상블 등가성과 샘플링이 까다로울 수 있습니다. thermostat이 목표분포는 맞춰도 실제 동역학을 왜곡할 수 있습니다.

재료·DFT로 이어지는 연결
  • MD의 NVE/NVT/NPT 선택과 열역학 출력 해석에 필수입니다.
  • GCMC 흡착과 grand-canonical DFT의 공통 언어입니다.
  • 자유에너지 계산·향상샘플링은 희귀상태 확률을 복원합니다.

06이 글의 용어 사전

본문을 읽다가 “그래서 이 말이 정확히 뭐지?”라는 질문이 생기면 여기에서 다시 확인할 수 있습니다.

앙상블ensemble
같은 거시 제약을 만족하는 미시상태들의 확률적 모음입니다.
분배함수partition function
상태의 Boltzmann 가중치를 합해 열역학량을 생성하는 함수입니다.
대정준 앙상블grand canonical ensemble
온도·부피·화학퍼텐셜을 고정하고 입자수 변동을 허용하는 앙상블입니다.
근사approximation
정확해를 직접 얻기 어려울 때 통제된 가정으로 문제를 계산 가능하게 바꾸는 선택입니다.

07참고자료와 더 읽을거리

  1. KAIST 화학과 대학원 교과목 개요 — 양자화학·통계열역학·반응동역학·표면화학·전기화학의 공식 범위를 따릅니다.
  2. MIT OpenCourseWare, Statistical Mechanics I — 앙상블·분배함수·양자통계를 대학원 수준의 문제와 함께 다룹니다.
  3. NIST CODATA Fundamental Physical Constants — 양자·통계 계산에 쓰는 기본 상수의 최신 권고값을 제공합니다.
NEXT CHEMISTRY NOTE양자통계와 phonon →