에르고딕 가설

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

1. 개요[편집]

통계역학 교과서 1장에서 “그렇다고 치자”로 넘어간 그 가정. 그런데 시뮬레이션을 돌리는 사람은 그 가정을 매일 실제로 검사해야 한다.

에르고딕 가설(ergodic hypothesis)은 한 궤도를 따라 충분히 오래 잰 시간 평균이 위상공간 전체에 대한 앙상블 평균과 같다는 가정이다.

\lim_{T\to\infty}\frac{1}{T}\int_0^T A(\mathbf{x}(t))\,dt \;=\; \int A(\mathbf{x})\,d\mu(\mathbf{x})

이 등식이 통계역학 전체의 접합부다. 왼쪽은 실험과 분자동역학이 실제로 재는 것이고, 오른쪽은 정준 앙상블이나 미시정준 앙상블 같은 이론이 계산하는 것이다. 둘이 같지 않으면 우리가 재는 값과 이론이 예측하는 값이 서로 다른 물건이 된다.

시뮬레이션 입장에서 이 가정이 중요한 이유는 철학이 아니라 비용 때문이다. 100만 개 배위를 독립적으로 뽑아 앙상블 평균을 내는 것은 불가능하지만, 궤적 하나를 100만 스텝 굴리는 것은 매일 하는 일이다. 후자로 전자를 대신할 수 있다는 면허가 곧 에르고딕성이며, 그 면허가 실제로 유효한지 확인하지 않고 평균을 내는 순간 결과는 그냥 틀린 값이다.

2. 볼츠만의 원안이 왜 틀렸나[편집]

볼츠만이 처음 쓴 형태는 훨씬 강했다. “궤도가 시간이 지나면 에너지 초곡면 위의 모든 점을 통과한다”는 것이었고, 그리스어 ergon(일) + hodos(길)에서 이름을 땄다. 문제는 이게 위상수학적으로 불가능하다는 것이다. 매끄러운 궤도는 1차원 곡선인데 $N$ 입자계의 에너지 초곡면은 $6N-1$ 차원이고, 1차원 곡선이 자기교차 없이 고차원 다양체를 전부 덮을 수는 없다. 1913년 로젠탈과 플랑슈렐이 이를 엄밀히 못 박았다.

그래서 에렌페스트 부부가 1911년에 내놓은 수정판이 준에르고딕 가설이다 — 모든 점을 지나는 게 아니라 모든 점에 임의로 가까이 간다(궤도가 조밀하다). 이건 가능하다. 다만 조밀성만으로는 평균이 수렴한다는 보장이 안 나와서, 진짜 해결은 20년 뒤 측도론에서 왔다.

3. 버코프 정리와 정확한 정의[편집]

1931년 버코프의 개별 에르고딕 정리가 결정판이다. 측도 보존 변환 $T$ 와 유한 측도 $\mu$ , $f \in L^1$ 에 대해 시간 평균

\bar f(x) = \lim_{n\to\infty}\frac{1}{n}\sum_{k=0}^{n-1} f(T^k x)

가 거의 모든 초기조건에서 존재한다는 것이 전반부다. “거의 모든”이 핵심이다 — 수렴하지 않는 초기조건은 반드시 있지만 측도가 0이라 무시해도 좋다. 후반부가 조건부다. 계가 에르고딕하면, 그리고 오직 그때만 $\bar f$ 가 초기조건에 의존하지 않는 상수 $\int f\,d\mu$ 가 된다. (폰 노이만의 평균 에르고딕 정리는 같은 결론을 $L^2$ 노름 수렴으로 준다.)

그리고 에르고딕성의 정확한 정의는 이것이다.

임의의 불변집합 $A$ ( $T^{-1}A = A$ )에 대해 $\mu(A) = 0$ 또는 $\mu(A) = 1$ .

즉 “궤도가 온 공간을 훑는다”가 아니라 **“위상공간을 둘 다 양의 측도인 두 조각으로 쪼갤 수 없다”**가 정의다. 이 형태가 왜 중요하냐면, 반례를 만드는 방법을 곧바로 알려 주기 때문이다. 양의 측도를 갖는 불변 영역을 하나만 찾아내면 그 계는 에르고딕하지 않다. 뒤에 나오는 반례들이 전부 이 방식이다.

4. 혼합과 그 위의 계층[편집]

에르고딕성은 이 동네에서 가장 약한 조건이다. 위로 갈수록 강해지는 계층이 있다.

\text{베르누이} \subset \text{콜모고로프}(K) \subset \text{혼합} \subset \text{에르고딕}

에르고딕: 시간 평균이 앙상블 평균으로 간다. 그게 전부다. 얼마나 빨리 가는지는 아무 말도 안 해 준다.
혼합(mixing): $\mu(A \cap T^{-n}B) \to \mu(A)\mu(B)$ . 초기 영역이 시간이 지나며 위상공간 전체에 균일하게 번진다는 뜻이고, 이게 있어야 상관함수가 감쇠하고 평형으로의 완화를 말할 수 있다. 에르고딕하지만 혼합이 아닌 표준 예가 원 위의 무리수 회전이다 — 궤도가 원 전체에 조밀하게 채워지지만 뭉친 덩어리는 영원히 뭉친 채 돌기만 한다.
콜모고로프·베르누이: 상관이 지수적으로 죽고, 극단적으로는 계가 동전 던지기와 동형이 된다. 시나이 당구나 하드스피어 기체가 이 부류에 속한다는 결과들이 있다.

실무에서 필요한 것은 사실 에르고딕성이 아니라 혼합, 그것도 유한 시간 안에 일어나는 혼합이다. “무한 시간이면 수렴한다”는 보장은 마감이 있는 사람에게 아무 쓸모가 없다.

5. 증명된 쪽[편집]

에르고딕성이 실제로 증명된 물리계는 놀랄 만큼 적고, 그 목록이 거의 다 당구공 계다. 시나이가 1970년에 원형 산란체가 박힌 사각 당구대(시나이 당구)의 에르고딕성과 혼합성을 증명한 것이 출발점이다. 볼록한 산란체의 곡면이 평행하게 들어온 궤도 다발을 매번 발산시켜서, 국소적으로 지수적 불안정성이 강제되는 구조다. 부니모비치는 곡면이 오목해도 초점 거리가 충분히 짧으면 같은 결론이 난다는 것을 보였다(스타디움 당구).

물리적으로 중요한 것은 이것이 강체구 기체와 같은 부류라는 점이다. $N$ 개의 강체구가 상자 안에서 탄성 충돌하는 계가 에르고딕하다는 “볼츠만-시나이 가설”은 오랫동안 부분적으로만 증명되다가 시마니의 작업으로 상당히 일반적인 경우까지 확장됐다. 분자동역학에서 쓰는 레너드-존스 유체는 강체구가 아니라 매끄러운 퍼텐셜이라 이 증명이 그대로 적용되지는 않지만, “짧은 거리에서 강하게 밀어내는 다체 계는 잘 섞인다”는 직관의 유일한 엄밀한 근거가 여기다. 우리가 MD 평균을 믿는 근거의 상당 부분이 사실 당구공 정리에 기대고 있는 셈이다.

6. 반례들 — 이게 예외가 아니라는 게 문제[편집]

에르고딕성이 성립하지 않는 예는 병적인 구성물이 아니라 우리가 매일 시뮬레이션하는 계들이다.

KAM 정리의 불변 토러스. 근적분가능 해밀토니안 계에서 살아남는 토러스들의 합집합은 불변집합이면서 측도가 $1 - O(\sqrt\epsilon)$ 이다. 0도 1도 아니므로 정의상 에르고딕하지 않다. 그것도 측도 0짜리 예외가 아니라 거의 전부가 반례라는 게 결정타였다. 에농-하일레스 계에서 이 그림을 눈으로 볼 수 있다.
적분가능계. 보존량이 자유도만큼 있으면 궤도는 애초에 저차원 토러스에 갇힌다. 조화진동자 사슬은 정규모드 에너지가 각각 보존돼서 에너지가 절대 재분배되지 않는다.
페르미-파스타-울람-칭구 문제. 1953년, 비선형 사슬을 굴리면 에너지가 모든 모드에 골고루 퍼지리라 기대하고 시작한 최초의 컴퓨터 수치실험이 정반대 결과를 냈다. 에너지가 몇 개 모드 사이에서 왔다 갔다 하다가 초기 상태로 거의 되돌아온 것이다. 등분배가 안 일어났다. 낮은 에너지에서는 KAM 토러스가, 어떤 문턱을 넘으면 공명 겹침이 일어나 비로소 등분배가 시작된다는 것이 이후의 정리된 그림이다.
앤더슨 국소화. 양자 쪽 반례. 무질서 격자에서 고유상태가 지수적으로 갇혀 버리면 파속이 초기 위치 근처를 영원히 벗어나지 못한다. 시간 평균한 밀도 분포가 균일 분포와 완전히 다르고, 다체 국소화 문제에서는 이것이 “고립계가 스스로 열평형에 도달하지 못한다”는 형태로 나타난다.

7. 분자동역학: 이 문서의 실전 부분[편집]

MD가 정준 앙상블 평균을 궤적의 시간 평균으로 대신 계산하는 근거가 정확히 에르고딕 가설이다. 그리고 실전에서 이 근거가 무너지는 방식은 대개 정리의 반례가 아니라 시간이 모자란 것이다.

계가 원리적으로 에르고딕해도, 자유에너지 지형에 높은 장벽이 있으면 궤적이 하나의 골짜기 안에서만 논다. 그 안에서는 관측량이 예쁘게 수렴하는 것처럼 보이고 오차 막대도 작게 나오지만, 값 자체가 틀렸다. 이걸 유사에르고딕성 결여(effective non-ergodicity)라 부르며, 단백질 접힘·유리 전이·결정 다형처럼 시간 척도가 갈라지는 계에서는 기본값에 가깝다. 무서운 점은 이 실패가 조용하다는 것 — 발산도 경고도 없이 그냥 잘 수렴한 것처럼 보인다.

서모스탯 선택도 여기에 직결된다. 노제-후버 서모스탯은 정준분포를 정확히 생성한다는 증명이 있지만, 그 증명이 에르고딕성을 가정한다. 그리고 그 가정이 하필 가장 단순한 계인 조화진동자에서 대놓고 깨진다 — 단일 노제-후버 진동자의 위상공간에는 규칙적 토러스와 카오스 영역이 공존해서, 초기조건에 따라 온도 분포가 정준분포와 다르게 나온다. 이 사실이 알려지면서 마르티나·클라인·터커먼이 1992년 열욕을 사슬로 잇는 노제-후버 체인을 내놓았고, 그게 지금 표준이 됐다. 서모스탯을 데우는 서모스탯을 다는 눈물겨운 구조인데, 목적은 하나 — 남는 불변 토러스를 깨서 에르고딕성을 억지로 확보하는 것이다. 랑주뱅 서모스탯이 요즘 인기 있는 이유도 무작위력이 결정론적 불변구조를 원천적으로 부숴 주기 때문이다.

관련해서 등분배 정리를 진단으로 쓰는 관행을 기억할 만하다. 자유도별 평균 운동에너지가 $k_BT/2$ 에서 고르게 나오는지 확인하는 것인데, 특정 모드만 온도가 다르면 그건 에너지가 재분배되지 않고 있다는 직접 증거다. “flying ice cube” 현상이 이 진단에 걸리는 대표 사례다.

8. 이산 판본: MCMC[편집]

마르코프 연쇄 몬테카를로에서 에르고딕성은 두 조건으로 번역된다. 기약성(어느 상태에서 어느 상태로든 유한 단계에 갈 수 있다)과 비주기성(주기적으로 갇혀 도는 부분집합이 없다)이다. 여기에 상세균형이 주는 불변성을 더하면 연쇄가 초기값과 무관하게 목표분포로 수렴하고 시간 평균이 기댓값으로 간다. 즉 MCMC의 수렴 정리는 버코프 정리의 마르코프 연쇄 판본이다.

그리고 실패 양상도 MD와 똑같다. 이론상 기약이어도(제안 분포에 양의 확률이 있으니까) 봉우리 사이 골짜기가 깊으면 현생에서는 못 넘어간다. 처방도 같은 계열이다.

병렬 템퍼링: 고온 복제본이 장벽을 넘어 다니고 교환으로 저온 사슬에 배위를 넘겨준다. 유효 에르고딕성을 사는 가장 직접적인 방법.
자유에너지 섭동: 여기서는 문제가 더 고약하다. $\Delta F = -k_BT\ln\langle e^{-\beta\Delta U}\rangle$ 의 지수 평균은 드물게 나타나는 낮은 $\Delta U$ 배위가 값을 지배한다. 즉 평균에 기여하는 배위가 표본에서 거의 안 나오는 꼬리에 몰려 있어서, 겉보기 수렴이 진짜 수렴과 특히 무관하다. 상태 겹침을 확인하지 않은 FEP 결과는 그냥 못 믿는다.

실전 점검은 결국 두 가지로 수렴한다. 겔만-루빈 진단은 서로 다른 초기점에서 출발한 여러 사슬(MD라면 여러 독립 궤적)의 사슬 간 분산과 사슬 내 분산을 비교한다. 초기조건 정보가 남아 있으면 $\hat R$ 이 1보다 크게 나오는데, 이건 정확히 “불변집합이 쪼개져 있는가”를 실험적으로 묻는 것이다. 단일 궤적으로는 원리적으로 잡을 수 없는 실패라 독립 복제본을 여러 개 돌리는 것이 유일한 방어다. 유효표본크기는 상관 때문에 실제로 몇 개의 독립 표본에 해당하는지를 세어, 100만 스텝짜리 궤적이 사실은 독립 표본 12개짜리였다는 사실을 알려 준다. 요즘 권장선은 순위 정규화 $\hat R < 1.01$ 과 파라미터당 ESS $> 400$ 정도다.¹²

9. 정리[편집]

에르고딕 가설은 세 층으로 나눠 이해하는 게 실용적이다. 수학적으로는 불변집합의 측도가 0 또는 1이라는 명확한 정의가 있고, 버코프 정리가 그 아래에서 시간 평균의 수렴을 보장한다. 물리적으로는 KAM 토러스라는 양의 측도 반례 때문에 일반 해밀토니안 계에서 성립하지 않으며, 자유도가 커질수록 반례의 측도가 줄어들어 실질적으로 성립하는 쪽으로 간다. 계산적으로는 성립 여부보다 유한한 계산 시간 안에 성립하는가가 유일하게 중요한 질문이고, 그 답은 정리가 아니라 $\hat R$ ·ESS·복제본 간 산포로 매번 새로 검사해야 한다.³⁴

10. 관련 문서[편집]

MD 논문에서 “100 ns 시뮬레이션 3회 반복”이라는 표기를 보면 그게 그냥 관행이 아니라 겔만-루빈을 손으로 하는 것이다. 세 궤적의 관측량이 서로 다른 값에서 평평해지면, 오차 막대를 아무리 예쁘게 그려도 그 결과는 폐기 대상이다. ↩
반대로 단일 궤적의 표준편차만 보고하는 습관이 위험한 이유도 같다. 한 골짜기 안에 갇힌 궤적의 요동은 작고 예쁘다. 정밀해 보이는 틀린 값만큼 리뷰를 잘 통과하는 것도 없다. ↩
볼츠만이 에르고딕 가설을 세운 게 1870년대인데, 그게 일반적으로 거짓임이 밝혀진 건 1950년대 KAM이다. 80년 동안 통계역학은 틀린 가정 위에서 대단히 정확한 예측을 내놓고 있었다. 요즘은 에르고딕성보다 “전형성”(typicality) — 거의 모든 미시상태가 이미 열평형처럼 보인다 — 쪽 설명이 더 근본적이라고 보는 시각이 우세하다. ↩
이 동네 최고의 농담: “에르고딕 가설은 시간 평균이 앙상블 평균과 같다는 가정이고, 실무에서는 내 시뮬레이션이 충분히 길다는 가정이다.” 둘 다 검증 없이 쓰인다는 점까지 정확히 같다. ↩