미시정준 앙상블

계산물리 분자동역학 통계 마지막 수정: 2026-08-02 05:11:42

1. 개요[편집]

미시정준 앙상블 Microcanonical Ensemble
다른 이름	NVE 앙상블
고정량	입자수 N, 부피 V, 전체 에너지 E
기본 공준	등확률 원리 (equal a priori probability)
특성 함수	엔트로피 S = k_B ln Ω(E)
온도	유도량 — 1/T = ∂S/∂E
주 구현	서모스탯 없는 MD + 심플렉틱 적분기

계에 아무것도 해 주지 않는 것. 그게 NVE다. 그래서 가장 정직하고, 그래서 가장 무섭다.

미시정준 앙상블(microcanonical ensemble)은 입자수 $N$ ·부피 $V$ ·전체 에너지 $E$ 가 모두 고정된 고립계의 미시상태 분포로, 주어진 에너지에서 접근 가능한 모든 미시상태가 똑같은 확률을 갖는다고 놓는 앙상블이다. 시뮬레이션 쪽에서는 흔히 NVE로 부르며, “서모스탯을 끄고 뉴턴 방정식만 그대로 적분한 MD”가 곧 이것이다.

범위 정리. 열저장조가 붙은 $NVT$ 계와 볼츠만 인자·분배함수 이야기는 정준 앙상블이, 온도를 실제로 붙잡는 알고리즘 비교는 서모스탯이, 입자 교환까지 허용하는 확장은 대정준 앙상블이 담당한다. 이 문서는 NVE 그 자체의 통계역학과, MD 실무에서 NVE가 왜 품질 검사 도구로 쓰이는가를 다룬다.

2. 등확률 원리와 볼츠만 엔트로피[편집]

미시정준 앙상블의 출발점은 유도가 아니라 공준이다. 에너지 $E$ 인 고립계에서 접근 가능한 미시상태가 $\Omega(E)$ 개일 때, 각 상태의 확률은 $1/\Omega(E)$ 로 전부 같다. 이 등확률 가정을 정당화하는 동역학적 근거가 리우빌 정리다 — 위상공간 흐름이 부피를 보존하므로 에너지 초곡면 위의 균등 측도가 해밀토니안 역학의 불변 측도가 된다. 즉 “균등 분포는 시간이 지나도 균등하게 남는” 유일하게 자연스러운 선택이다.

여기서 통계역학 전체의 주춧돌이 나온다.

S(N,V,E) = k_B \ln \Omega(N,V,E)

볼츠만의 묘비에 새겨진 그 식이며, 나머지 열역학은 전부 미분이다.

\frac{1}{T} = \left(\frac{\partial S}{\partial E}\right)_{N,V}, \qquad \frac{p}{T} = \left(\frac{\partial S}{\partial V}\right)_{N,E}, \qquad \frac{\mu}{T} = -\left(\frac{\partial S}{\partial N}\right)_{V,E}

온도가 입력이 아니라 출력이라는 점이 정준 앙상블과의 결정적 차이다. NVE MD를 돌려서 “온도 300 K로 설정했다”는 말은 엄밀히는 성립하지 않는다. 초기 속도로 에너지를 넣어 두면 계가 알아서 어떤 온도로 정착하는 것이며, 얼음이 녹거나 단백질이 풀리면서 퍼텐셜 에너지를 빨아들이면 온도는 그만큼 내려간다.

연속계에서는 미시상태를 셀 수 없으니 위상공간 부피로 대신하고, 폭 $\Delta E$ 짜리 에너지 껍질을 잡는다.

\Omega(E) = \frac{1}{N!\,h^{3N}} \int \delta\big(H(\mathbf{q},\mathbf{p}) - E\big)\, d^{3N}q\, d^{3N}p

껍질 폭 $\Delta E$ 를 얼마로 잡느냐는 임의적인데, 이 자의성이 열역학 극한에서 자동으로 사라진다. $\Omega$ 는 $N$ 에 대해 지수적으로 크므로 $\ln \Omega \sim O(N)$ 인 반면 $\ln(\Delta E/E)$ 는 기껏해야 $O(\ln N)$ 이고, 입자당 엔트로피로 나누면 $O(\ln N / N) \to 0$ 이다. 다르게 말하면 $\Omega(E)$ 를 껍질로 세든 그 이하 전체 부피로 세든 같은 답이 나온다.¹

3. 앙상블 등가성과 그 붕괴[편집]

열역학 극한에서 NVE와 NVT가 같은 상태방정식을 준다는 것이 앙상블 등가성이다. 근거는 단순하다 — 정준 앙상블의 에너지 상대 요동이 $\Delta E/\langle E\rangle \sim N^{-1/2}$ 이므로 $N$ 이 커지면 에너지가 사실상 고정되고, 결국 미시정준과 구분이 안 된다. 그래서 거시계에서는 “계산하기 편한 앙상블을 고르면 된다”가 통한다.

문제는 등가성이 보장이 아니라 조건부라는 것이다. 세 곳에서 깨진다.

유한계. 입자 수백 개짜리 클러스터에서는 $N^{-1/2}$ 가 몇 퍼센트라 앙상블 차이가 결과에 그대로 남는다.
장거리 상호작용. 중력처럼 $1/r$ 로 감쇠하는 계는 에너지가 시성적(extensive)이지 않아 열역학 극한 자체가 존재하지 않는다.
1차 상전이 근방. 미시정준 엔트로피 $S(E)$ 가 볼록한 구간(convex intruder)을 가질 수 있는데, 정준 앙상블은 이 구간을 르장드르 변환 과정에서 통째로 지워 버린다(맥스웰 구성).

그리고 등가성 붕괴의 가장 화려한 결과가 음의 열용량이다. 미시정준에서 $C = (\partial T/\partial E)^{-1}$ 이므로 $S(E)$ 가 볼록하면 에너지를 넣었는데 온도가 내려가는 구간이 생긴다. 자기중력계가 교과서적 예로, 비리얼 정리에 의해 $\langle E \rangle = -\langle K \rangle$ 이므로 별 무리가 에너지를 잃으면 수축하면서 오히려 뜨거워진다 — 별이 스스로 진화하는 이유이자 중력 열역학이 골치 아픈 이유다. 나트륨 클러스터의 녹는 전이에서도 음의 열용량 구간이 실측됐다.

반면 정준 앙상블에서는 $\langle \Delta E^2\rangle = k_BT^2 C_V \ge 0$ 이므로 음의 열용량이 원리적으로 불가능하다. 분산이 음수일 수 없기 때문이다. 같은 물리계인데 앙상블에 따라 나오는 답이 다르다는 사실을 처음 만나면 대개 버그를 의심하는데, 이건 진짜다.²

4. 요동으로 열용량 재기[편집]

정준 앙상블에서 $C_V$ 를 에너지 요동으로 재는 것은 유명하지만, NVE에서는 전체 에너지가 상수라 요동이 0이다. 대신 운동에너지와 퍼텐셜 에너지가 서로 주고받는 요동이 남고, 여기서 $C_V$ 를 뽑는 것이 르보비츠-퍼커스-베를레(1967) 관계다. 3차원 $N$ 입자계에서

\frac{\langle \delta K^2 \rangle}{\langle K \rangle^2} = \frac{2}{3N}\left(1 - \frac{3Nk_B}{2C_V}\right)

이 성립한다. 우변 괄호의 구조를 보면 이 식이 무슨 말을 하는지 바로 보인다. 이상기체는 $C_V = \frac{3}{2}Nk_B$ 라 괄호가 정확히 0 — 상호작용이 없으면 $K = E$ 라 운동에너지가 아예 요동칠 수 없다. 상호작용이 강해져 $C_V$ 가 커질수록 괄호가 1에 가까워지고, 그때 운동에너지 요동은 정준 앙상블 값 $\langle \delta K^2\rangle_{NVT} = \frac{3}{2}Nk_B^2T^2$ 에 접근한다. NVE의 운동에너지 요동은 항상 NVT보다 작다는 것이 결론이며, 이걸 모르고 NVE 궤적에 정준 공식을 그대로 들이대면 $C_V$ 가 체계적으로 틀린다.

실무적으로는 위 식을 뒤집어 $C_V$ 를 얻는다. 온도를 여러 개 돌려 $\langle E\rangle(T)$ 를 수치미분하는 것보다 한 번의 긴 NVE 궤적으로 끝나므로 싸지만, 분산 추정이라 수렴이 느리고 상전이 근방에서는 요동 자체가 커져 유효표본크기가 폭락한다.

5. MD에서의 NVE — 에너지 보존이 곧 성적표[편집]

서모스탯을 끄면 그게 NVE다. 이 단순한 사실 덕분에 NVE는 앙상블이면서 동시에 적분 품질의 유일한 실전 지표로 쓰인다. NVT로 돌리면 서모스탯이 에너지 오차를 조용히 흡수해 버려서, 시간 스텝이 두 배 크든 컷오프 처리가 엉망이든 온도 그래프는 태연하게 300 K를 가리킨다. 서모스탯을 떼고 총에너지를 그려 봐야 비로소 코드가 진짜로 뉴턴 방정식을 풀고 있는지가 드러난다.

여기서 베를레 적분·립프로그·속도 베를레가 왜 60년째 표준인지가 나온다. 이들은 심플렉틱 적분기이자 시간 가역적이라, 실제 해밀토니안 $H$ 대신 그것과 $O(\Delta t^2)$ 만큼 떨어진 섀도 해밀토니안을 (지수적으로 작은 오차를 빼면) 정확히 보존한다. 결과적으로 에너지 오차가 단조 증가하지 않고 유계로 진동한다. 4차 룽게-쿠타는 한 스텝 정확도는 훨씬 좋지만 심플렉틱하지 않아 에너지가 서서히 표류하고, 나노초 단위로 돌리는 MD에서는 이 표류가 정확도를 압도한다. “정확한 적분기”가 아니라 “구조를 보존하는 적분기”를 고르는 대표적인 사례.

실전 판정 기준은 상대 에너지 표류다. 총에너지의 표류를 선형 회귀로 뽑아 $|\Delta E|/E$ 가 나노초당 $10^{-5}$ 수준 이하면 대체로 통과로 본다(요구 정밀도에 따라 더 조이기도 한다). 이 값이 나쁘면 원인은 거의 정해져 있다.

시간 스텝. 가장 빠른 진동(X–H 신축, 주기 10 fs 안팎)의 1/10 이하가 안전선이라 무구속 계는 1 fs, 결합을 구속하면 2 fs가 국룰이다.
컷오프 불연속. 퍼텐셜을 컷오프에서 뚝 자르면 입자가 경계를 넘을 때마다 힘이 튀면서 에너지가 계단식으로 주입된다. 시프트·스위치 함수나 부드러운 감쇠가 필수이며, 이건 레너드-존스 퍼텐셜 같은 단거리 항에서 특히 티가 난다.
정전기 파라미터. 에발트 합산/PME의 격자 간격, 실공간 컷오프, 스플라인 차수, 상대 허용오차가 헐거우면 힘 오차가 잡음처럼 들어와 표류로 나타난다.
구속 알고리즘 허용오차. SHAKE/LINCS의 수렴 기준을 느슨하게 두면 결합 길이가 미세하게 어긋나며 에너지가 샌다.
이웃 목록 갱신 주기와 정밀도. 버퍼 없이 목록을 오래 재사용하면 상호작용 쌍이 누락되고, 단정밀도 누적에서는 부동소수점 연산 오차가 무작위걷기처럼 쌓인다.

즉 에너지 표류는 하나의 숫자에 온갖 설정 실수가 모여 드는 통합 진단 지표다. 새 힘장이나 새 코드를 세팅하면 짧은 NVE를 먼저 돌려 이 숫자부터 확인하는 것이 순서다.

6. 순간 온도, 에르고딕성, 그리고 함정[편집]

NVE에도 온도는 있다. 등분배 정리로 정의되는 순간 온도

T(t) = \frac{2K(t)}{N_f k_B}

이며, $N_f$ 는 자유도 수다. 여기서 자유도를 세는 게 은근한 함정으로, 무게중심 운동량이 보존량이라 병진 3자유도를 빼야 하고( $N_f = 3N-3$ ), 구속 결합이 있으면 그 수만큼 더 뺀다. 이걸 틀리면 온도가 체계적으로 몇 퍼센트 어긋난다. 그리고 $T(t)$ 는 요동치는 양이며, 그 폭이 앞 절의 요동 관계로 정해진다 — 입자 1000개짜리 상자에서 순간 온도가 ±10 K 흔들리는 것은 지극히 정상이다.

NVE MD가 미시정준 평균을 준다는 보장은 에르고딕 가설에 기대고 있다. 하나의 궤적이 에너지 초곡면을 충분히 골고루 훑어야 시간평균 = 앙상블평균이 성립하는데, 이건 정리가 아니라 희망이다. KAM 정리가 말해 주듯 약한 비선형 계에서는 불변 원환면이 살아남아 궤적이 위상공간의 일부에만 갇히고, 페르미-파스타-울람-칭구 문제가 바로 그 반례다. 실제 MD에서는 유리 상태·단백질 접힘처럼 장벽이 높은 계가 같은 병을 앓으며, 이때는 앙상블을 바꾸든(병렬 템퍼링) 표본추출 편향을 넣든 해야 한다.

마지막으로 자주 헷갈리는 것 하나. “플라잉 아이스 큐브”는 NVE의 병이 아니다. 그건 균일 속도 스케일링 서모스탯이 무게중심 운동을 열운동으로 착각해 내부 자유도를 얼려 버리는 현상이라 NVT 쪽 사고이고(등분배 정리 참고), 순수 NVE에서는 총운동량 자체가 보존량이라 그런 누적이 일어나지 않는다. 다만 초기 속도를 뽑을 때 무게중심 운동량을 0으로 맞추지 않으면 그 운동에너지가 영원히 계에 남아 온도를 부풀리므로, 초기화 단계에서 한 번 제거하고 $N_f$ 에서 3을 빼는 처리는 NVE에서도 그대로 필요하다.

앙상블 선택의 실무 요약은 간단하다. 적분기 검증·에너지 보존 확인·충격/충돌처럼 열욕이 물리적으로 없는 문제는 NVE, 실험 조건을 흉내 내는 생산 계산은 NVT나 등온등압 앙상블, 입자 수가 변해야 하면 대정준 앙상블이다.³⁴

7. 관련 문서[편집]

“껍질을 어떻게 잡아도 상관없다”는 이 무신경함이 처음엔 사기처럼 보이지만, 지수함수의 위력을 생각하면 당연하다. $\Omega \sim e^{10^{23}}$ 짜리 수를 다루는 판에 앞에 곱해진 상수 몇 개는 소수점 아래 스물몇째 자리에서 소멸한다. 참고로 볼츠만 엔트로피(껍질)와 깁스 엔트로피(누적 부피) 중 어느 쪽이 “옳은” 정의인가는 음의 절대온도 논쟁에서 다시 불붙었고, 아직도 완전히 정리되지 않았다. ↩
음의 열용량을 처음 들으면 “그럼 에너지를 계속 빼면 무한히 뜨거워지나” 싶은데, 실제로 자기중력계는 그 방향으로 폭주한다(중력 열재앙, gravothermal catastrophe). 구상성단 중심부가 붕괴하는 시나리오가 이것이며, 통계역학이 “안정 평형이 없는 계”를 진지하게 다뤄야 하는 몇 안 되는 무대다. ↩
NVE를 “설정할 게 없으니 제일 쉬운 앙상블”로 여기는 오해가 흔한데, 실상은 정반대다. 서모스탯이 있으면 어지간한 에너지 오차는 열욕이 먹어 치우고 결과가 그럴싸하게 나오지만, NVE는 코드의 모든 부실을 총에너지 곡선 하나에 정직하게 적어 놓는다. NVE가 잘 돌지 않는 세팅은 NVT에서도 이미 틀려 있고, 다만 안 보일 뿐이다. ↩
논문 메서드 절에서 “NVE 앙상블에서 100 ps 평형화 후”라는 문장을 보면 한 번쯤 갸웃해도 좋다. NVE는 원리적으로 온도를 조절할 수단이 없으므로 평형화 단계와는 궁합이 나쁘다. 대개는 NVT로 온도를 잡은 뒤 그 마지막 상태에서 NVE 생산 구간을 시작한다는 뜻인데, 그 과정을 뭉뚱그려 쓴 것이다. ↩