H 정리

1. 개요[편집]

H 정리 H-theorem
발표	루트비히 볼츠만, 1872
대상	볼츠만 방정식의 해 $f(\mathbf x,\mathbf v,t)$
H 함수	$H(t)=\int f\ln f\, d^3v$
주장	$dH/dt\le 0$
등호 조건	$f$ 가 국소 맥스웰-볼츠만 분포
숨은 가정	분자 혼돈 가정 (Stosszahlansatz)
주요 반론	로슈미트 역설 (1876) · 체르멜로 재귀 반론 (1896)

볼츠만은 제2법칙을 뉴턴 역학에서 유도했다고 믿었다. 사실은 유도한 것이 아니라, 가정 한 줄에 미리 넣어 둔 것을 되찾은 것이었다.

H 정리(H-theorem)는 볼츠만 방정식을 따르는 기체의 속도분포함수에 대해 정의된 범함수 $H=\int f\ln f\,d^3v$ 가 시간에 대해 결코 증가하지 않으며, 등호는 분포가 맥스웰-볼츠만 분포일 때만 성립한다는 정리다. $S=-k_B H$ 로 두면 이것이 곧 엔트로피 증가이므로, 1872년 볼츠만이 발표했을 때의 의미는 명확했다 — 열역학 제2법칙을 분자 운동에서 유도했다는 선언이었다.

그 선언은 절반만 옳았다. 정리 자체는 참이지만, 증명의 어딘가에서 시간 방향이 이미 입력되어 있었기 때문이다. 이 문서는 (1) $H$ 의 정의와 $dH/dt\le0$ 증명의 구조, (2) 등호 조건이 왜 맥스웰 분포인지, (3) 시간 비대칭이 정확히 어느 줄에서 들어오는지, 그리고 (4) DSMC와 격자 볼츠만이 이 정리를 이산 수준에서 어떻게 물려받는지를 다룬다. 반론의 세부(속도를 뒤집으면 어떻게 되는가, 푸앵카레 재귀 시간은 얼마인가)는 로슈미트 역설이 전담하므로 여기서는 구조만 짚는다.

2. H 함수의 정의[편집]

공간 균일한 기체를 생각하자. 속도분포함수 $f(\mathbf v,t)$ 에 대해

$$H(t) \;=\; \int_{\mathbb R^3} f(\mathbf v,t)\,\ln f(\mathbf v,t)\; d^3v$$

공간 의존성이 있으면 위치에 대해서도 적분해 $H(t)=\int\!\!\int f\ln f\,d^3v\,d^3x$ 로 정의한다. 이 경우 볼츠만 방정식의 수송항 $\mathbf v\cdot\nabla_x f$ 가 주기·경면반사 경계에서 전 적분에 0을 기여하고, 속도에 의존하지 않는 외력항도 마찬가지로 0이 되므로 충돌항만 $H$ 를 움직인다. 이것이 정리 전체의 무대 설정이다.

부호 규약을 먼저 정리해 두는 게 좋다.

양	정의	평형에서	시간에 따라
$H$	$\int f\ln f$	최소	감소
$S=-k_B H$	단위 부피당 엔트로피	최대	증가

즉 $H$ 는 부호를 뒤집은 엔트로피 밀도다. 실제로 $f$ 에 맥스웰 분포를 넣고 적분하면 $-k_BH$ 가 이상기체의 엔트로피 밀도와 $h$ 와 $N!$ 에서 오는 상수 차이만 남기고 일치한다. 자세한 대응은 엔트로피 문서 참고.¹

3. dH/dt ≤ 0 증명의 구조[편집]

증명은 세 단계이고, 각 단계가 딱 한 가지 사실만 쓴다.

1단계 — 미분한다. $\int f\,d^3v=n$ 이 보존되므로 $\int \partial_t f\,d^3v=0$ 이고, 따라서

$$\frac{dH}{dt} = \int (1+\ln f)\,\frac{\partial f}{\partial t}\,d^3v = \int \ln f\;Q(f,f)\,d^3v$$

2단계 — 대칭화한다. 충돌적분 $Q(f,f)$ 는 네 개의 속도 $(\mathbf v,\mathbf v_*)\leftrightarrow(\mathbf v',\mathbf v'_*)$ 를 짝짓는 구조이고, 미시 가역성(단면적이 정·역 충돌에 대해 같음)과 적분변수의 교환 $\mathbf v\leftrightarrow\mathbf v_*$, $(\mathbf v,\mathbf v_*)\leftrightarrow(\mathbf v',\mathbf v'_*)$ 을 이용해 같은 적분을 네 가지 방식으로 쓴 뒤 평균하면

$$\frac{dH}{dt} = -\frac14\int\!\!\int\!\!\int \Big[\ln(f'f'_*)-\ln(f f_*)\Big]\Big[f'f'_* - f f_*\Big]\, |\mathbf v-\mathbf v_*|\,\sigma\;d\Omega\,d^3v_*\,d^3v$$

3단계 — 부호를 읽는다. 로그는 단조증가 함수이므로 임의의 양수 $a,b$ 에 대해 $(\ln a-\ln b)(a-b)\ge 0$ 이다. 대괄호 두 개의 곱이 언제나 음이 아니고, 나머지 인자($|\mathbf v-\mathbf v_*|$, $\sigma$)도 음이 아니며, 앞에 마이너스가 붙어 있으므로

$$\frac{dH}{dt}\ \le\ 0$$

증명의 전체 무게가 “로그가 단조증가한다”는 한 줄에 실려 있다는 점이 이 정리의 인상적인 부분이다. 어려운 곳은 3단계가 아니라 2단계의 대칭화이고, 그 대칭화가 성립하려면 충돌적분이 $f$ 의 곱 꼴로 쓰여 있어야 한다. 바로 그 곱이 문제의 씨앗이다.

4. 등호 조건 — 왜 하필 맥스웰 분포인가[편집]

$dH/dt=0$ 이 되려면 적분의 피적분함수가 거의 모든 곳에서 0이어야 하고, 로그의 단조성 때문에 그것은

$$f(\mathbf v')f(\mathbf v'_*) = f(\mathbf v)f(\mathbf v_*)$$

가 모든 충돌 가능한 조합에 대해 성립한다는 뜻이다. 양변에 로그를 씌우면

$$\ln f(\mathbf v') + \ln f(\mathbf v'_*) = \ln f(\mathbf v)+\ln f(\mathbf v_*)$$

즉 $\ln f$ 가 충돌 불변량이라는 조건이 된다. 탄성충돌에서 보존되는 가법적 양은 질량(상수 1), 운동량 $\mathbf v$, 운동에너지 $|\mathbf v|^2$ 다섯 개뿐이므로

$$\ln f = a + \mathbf b\cdot\mathbf v + c\,|\mathbf v|^2$$

이고, 정규화 가능하려면 $c<0$ 이어야 한다. 제곱을 완성해 정리하면 정확히

$$f^{\rm eq}(\mathbf v) = n\left(\frac{m}{2\pi k_B T}\right)^{3/2}\exp\!\left(-\frac{m|\mathbf v-\mathbf u|^2}{2k_BT}\right)$$

평형통계역학이 가정으로 깔고 들어가는 맥스웰-볼츠만 분포를, 운동론에서는 결과로 얻는다. 이것이 H 정리의 가장 논쟁의 여지가 없는 성과다. 게다가 “다섯 개의 충돌 불변량 $\Rightarrow$ 다섯 개의 보존형 거시 방정식”이라는 대응까지 여기서 함께 나온다.

BGK 모형에서도 같은 결론이 훨씬 쉽게 나온다. $Q=-(f-f^{\rm eq})/\tau$ 를 넣으면 $f$ 와 $f^{\rm eq}$ 가 다섯 모멘트를 공유하므로 $\int(f-f^{\rm eq})\ln f^{\rm eq}\,d^3v=0$ 이고

$$\frac{dH}{dt} = -\frac1\tau\int (f-f^{\rm eq})\left(\ln f-\ln f^{\rm eq}\right)d^3v \ \le\ 0$$

이 즉시 따른다. 피적분함수가 점별로 음이 아니기 때문이다. 충돌적분을 뺄셈 하나로 바꿔도 H 정리가 살아남는다는 사실이 BGK가 장난감 모형 이상으로 대접받는 이유 중 하나다. 이 형태는 이완의 시간 규모까지 노골적으로 보여 준다 — $H-H^{\rm eq}$ 가 대략 $e^{-t/\tau}$ 로 잦아들고, $\tau$ 는 평균 충돌 시간 규모다.²

5. 상대 엔트로피로 다시 읽기[편집]

$H$ 를 절대적인 양으로 보면 로그 안의 차원 때문에 찜찜하지만, 평형과의 차이로 보면 깔끔해진다. $f$ 와 $f^{\rm eq}$ 가 같은 다섯 모멘트를 가질 때 $\int f\ln f^{\rm eq}=\int f^{\rm eq}\ln f^{\rm eq}$ 이므로

$$H(f)-H(f^{\rm eq}) = \int f\ln\frac{f}{f^{\rm eq}}\,d^3v \;=\; n\,D_{\rm KL}\!\left(\tfrac{f}{n}\,\Big\Vert\,\tfrac{f^{\rm eq}}{n}\right) \ \ge\ 0$$

즉 H 정리는 “분포가 평형 분포로부터의 쿨백-라이블러 발산을 단조 감소시킨다”는 진술이다. 이 형태의 좋은 점이 셋이다. 첫째, 좌표 변환에 대해 불변이라 차원 문제가 사라진다. 둘째, 리아푸노프 함수로서의 역할이 명시적이 된다 — 아래로 유계이고 단조 감소하므로 평형이 유일한 도달점임이 따라온다. 셋째, 최대 엔트로피 원리와의 관계가 정확히 드러난다 — 주어진 모멘트 구속 아래 $H$ 를 최소화하는 분포가 곧 맥스웰 분포다.

정량적 이완 속도까지 묻는 것이 현대 연구 주제다. 체르파티·빌라니 계열의 작업은 $H(f)-H(f^{\rm eq})$ 와 그 소산률 사이의 함수 부등식(엔트로피-엔트로피 소산 부등식)을 세워 대수적 혹은 지수적 수렴률을 증명하는데, 하드 스피어 같은 실제 단면적에서 최적 상수를 얻는 문제는 아직 완전히 닫혀 있지 않다.

6. 시간 비대칭은 어디서 들어왔나[편집]

이제 급소다. 미시 역학은 시간 반전에 대칭인데 어떻게 $dH/dt\le0$ 이라는 방향성 있는 결론이 나왔을까?

범인은 충돌적분을 세울 때 쓴 분자 혼돈 가정(Stosszahlansatz)이다. 이체 분포를 단일입자 분포의 곱으로 인수분해하는

$$f_2(\mathbf v,\mathbf v_*,t)\ \approx\ f(\mathbf v,t)\,f(\mathbf v_*,t)$$

이 가정은 “충돌 직전에” 두 분자의 속도가 무상관이라고 요구한다. 여기서 “직전”이 전부다. 충돌 후의 두 분자는 산란각을 통해 서로의 속도를 알고 있으므로 명백히 상관되어 있는데, 가정은 그 상관을 배제하지 않는다. 전은 무상관, 후는 상관 — 시간에 대해 비대칭인 조건을 넣었으니 비대칭인 결론이 나오는 것이 당연하다.

증명 구조를 되짚으면 이 사실이 더 선명하다. 2단계의 대칭화가 성립한 이유는 충돌항이 $f'f'_*-ff_*$ 라는 곱의 차이 꼴이었기 때문이고, 그 곱 꼴 자체가 분자 혼돈 가정의 산물이다. 즉 비가역성은 3단계의 로그 부등식에서 나온 것처럼 보이지만, 실제로는 그 부등식을 쓸 수 있게 판을 깔아 준 1단계 이전의 가정에서 왔다.

따라서 H 정리는 역학 정리가 아니라 통계적 진술이다. “모든 미시상태에서 $H$ 가 감소한다”가 아니라 “압도적 다수의 미시상태에서 감소한다”이며, 로슈미트가 만든 반전 상태나 체르멜로가 지적한 재귀 시점은 그 다수에 속하지 않는 예외다. 볼츠만 본인이 1877년에 이 입장으로 정리해 답했고, 오늘날 이 답이 표준으로 받아들여진다. 반전 상태에서 실제로 무슨 일이 일어나는지, 재귀 시간이 왜 $e^{cN}$ 규모라 물리적으로 무의미한지, 조대화와 거시상태 셈이 어떻게 그림을 완성하는지는 로슈미트 역설 문서가 전담한다.

수학 쪽 결론도 같은 방향을 가리킨다. 랜퍼드(1975)는 볼츠만-그라드 극한에서 거의 모든 초기조건에 대해 짧은 시간 동안 볼츠만 방정식이 엄밀히 유도됨을 증명했다. 두 개의 단서가 각각 예외 집합을 측도 0 으로 밀어내는 장치와, 상관이 아직 축적되지 않은 구간이라는 제약에 대응한다. 그 시간 제약은 50년이 지난 지금도 본질적으로 풀리지 않았다.

7. 이산 수준의 H 정리[편집]

정리가 연속 방정식의 성질이라면, 그것을 이산화한 코드에서도 살아남을까? 이 질문이 실무에서 꽤 중요하다.

7.1. DSMC[편집]

직접 시뮬레이션 몬테카를로는 분포함수를 대표입자 표본으로 대신하고, 충돌적분을 확률적 충돌 표본추출로 대체한다. 여기서 분자 혼돈 가정은 알고리즘에 구조적으로 내장된다 — 충돌 상대를 같은 셀 안에서 위치 상관 없이 무작위로 고르기 때문이다. 그래서 DSMC는 볼츠만 방정식의 비가역성을 그대로 물려받는다. 충돌 표본추출이 상세 균형을 만족하도록 설계되어 있으므로 평형 분포가 맥스웰 분포로 고정되고, 속도 히스토그램에서 계산한 $H$ 는 평균적으로 감소한다.

다만 통계 잡음 때문에 $H(t)$ 가 매끄럽게 내려오지 않는다. 유한한 대표입자 수 $N_p$ 때문에 $H$ 추정치에 $O(1/N_p)$ 급의 편향과 요동이 얹히고, 평형 근처에서는 그 잡음 바닥 위에서 $H$ 가 진동한다. 실무 감각으로는 이렇다 — $H(t)$ 곡선이 단조 감소 뒤 잡음 수준에서 평평해지면 정상, 평형 도달 후 계속 표류하면 셀 크기·시간간격·충돌 짝짓기 중 하나가 잘못된 신호다. 셀 크기가 평균자유행로보다 크면 실제보다 강한 충돌을 세게 되어 수치적 과잉 소산이 생기는데, 이것이 $H$ 감소율의 이상으로 먼저 드러난다.

7.2. 격자 볼츠만[편집]

격자 볼츠만 방법은 사정이 더 미묘하다. 이산 속도 집합 $\{\mathbf c_i\}$ 와 가중치 $\{w_i\}$ 에 대한 자연스러운 이산 H 함수는

$$H = \sum_i f_i \ln\frac{f_i}{w_i}$$

이고, 이것이 모멘트 구속 아래 최소화될 때의 분포가 이산 평형이다. 문제는 표준 LBGK가 쓰는 평형분포가 이 최소화 문제의 해가 아니라, 맥스웰 분포를 마하수에 대해 2차까지 자른 다항식이라는 점이다. 다항식 평형은 음수 값을 가질 수도 있고, 결과적으로 표준 LBM에는 이산 H 정리가 없다. 저완화시간(고 레이놀즈수) 영역에서 LBM이 유명하게 불안정한 이유가 여기에 있다.

이를 정면으로 고친 것이 엔트로피 격자 볼츠만(ELBM)이다.

• 평형분포를 다항식 대신 위 $H$ 의 실제 최소화 해로 잡는다. D2Q9·D3Q19에서는 닫힌 형태가 알려져 있다.
• 충돌을 $f_i \to f_i + \alpha\beta\,(f_i^{\rm eq}-f_i)$ 로 쓰고, $H(f+\alpha(f^{\rm eq}-f))=H(f)$ 를 만족하는 $\alpha$ 를 셀마다 비선형 방정식으로 풀어 잡는다. 이러면 이산 엔트로피가 절대 잘못된 방향으로 가지 않는다.
• 대가는 매 셀·매 스텝의 스칼라 뉴턴 반복이다. 안정성을 계산량으로 산 셈이고, 그래서 실무에서는 필요한 영역에만 걸거나 다중완화시간(MRT)·정칙화 같은 더 싼 안정화로 대체하는 경우가 많다.³

여기서 얻는 교훈이 이 문서에서 가장 실용적인 대목이다. 연속 방정식이 가진 부등식이 이산화 후에도 성립하는지는 별도로 증명해야 하는 사항이며, 공짜로 따라오지 않는다. 같은 구조가 압축성 CFD 쪽에서는 이산 엔트로피 부등식과 엔트로피 조건으로, 분자 시뮬레이션 쪽에서는 심플렉틱 적분기의 구조 보존으로 반복해서 나타난다.

7.3. 그 밖의 이산 판본[편집]

• 채프먼-엔스코그 관점. 채프먼-엔스코그 전개의 1차 항을 H 정리에 넣으면 $-k_B\,dH/dt$ 가 점성 소산 항과 열전도 소산 항의 합으로 정리되고, 각 계수가 음이 아니라는 조건이 곧 $\mu\ge0$, $\kappa\ge0$ 이다. 나비에-스토크스의 소산 부등식이 H 정리의 크누센수 1차 근사인 셈이다.
• 격자 기체·이산 속도법. 이산 속도법(DVM)은 속도공간만 자르므로 잘 설계하면 이산 H 정리를 보존할 수 있고, 이것이 저속 미세유동에서 DVM이 신뢰받는 근거 중 하나다.
• 확률 모형. 랑주뱅 동역학이나 포커-플랑크 방정식 계열에서도 같은 구조의 정리가 있다. 상대 엔트로피가 감소하고 그 소산률이 피셔 정보로 표현되는 형태이며, H 정리의 확산 판본으로 봐도 무방하다.

8. 여담[편집]

• 볼츠만은 원래 이 양을 $E$ 로 적었다. 지금의 $H$ 표기는 1890년대 영국 쪽 논문에서 굳어진 것이고, 대문자 그리스 문자 에타(Η)였다는 설과 그냥 라틴 문자 H 였다는 설이 있는데 확실한 결론은 없다. 정리 이름이 “H”라는 알파벳 하나인 탓에 검색이 유난히 어려운 정리이기도 하다.
• “볼츠만이 제2법칙을 증명했다”와 “볼츠만의 증명은 순환논법이다”는 둘 다 흔한 요약인데 둘 다 부정확하다. 정확한 문장은 **“볼츠만은 하나의 통계적 가정으로부터 제2법칙을 유도했고, 그 가정이 어디서 오는지는 그 뒤 150년의 연구 주제가 되었다”**이다. 순환논법은 아니다 — 가정과 결론이 같은 문장이 아니기 때문이다.
• 이 정리를 처음 배우는 사람이 가장 많이 하는 오해가 ”$H$ 는 계의 엔트로피”라는 것이다. $-k_BH$ 는 단일입자 분포로만 계산한 엔트로피라서, 분자 사이 상관에 저장된 몫이 통째로 빠져 있다. 평형 근처의 희박기체에서는 그 몫이 작아 문제가 되지 않지만, 조밀한 액체나 강한 상관이 있는 계에서는 $-k_BH$ 를 엔트로피라고 부르면 안 된다. 로슈미트가 뒤집은 상태에서 정보가 “사라진” 것처럼 보이는 이유도 정확히 이것이다 — 정보는 상관으로 옮겨갔고 $H$ 는 그것을 보지 못한다.
• 실무 CFD·운동론 코드에서 $H(t)$ 를 찍어 보는 것은 값싼 검증 도구다. 초기 비평형 분포에서 출발해 $H$ 가 단조 감소하며 맥스웰 값으로 수렴하는지만 봐도 충돌 연산자·경계조건·시간적분의 버그가 상당수 걸러진다. 정리가 곧 회귀 테스트가 되는 드문 사례.

9. 관련 문서[편집]

• 볼츠만 방정식 · 기체운동론 · 맥스웰-볼츠만 분포 · 채프먼-엔스코그 전개
• 엔트로피 · 열역학 제2법칙 · 로슈미트 역설 · 시간의 화살
• 직접 시뮬레이션 몬테카를로 · 격자 볼츠만 방법 · 크누센수 · 평균자유행로
• 쿨백-라이블러 발산 · 최대 엔트로피 원리 · 젠센 부등식
• 리우빌 정리 · 에르고딕 가설 · 미시정준 앙상블
• 엔트로피 조건 · 랑주뱅 동역학 · 포커-플랑크 방정식 · 요동정리

10. Footnotes[편집]

1. 상수 차이가 나는 이유는 고전 운동론이 위상공간의 “칸 크기”를 정할 수 없기 때문이다. $\ln f$ 의 $f$ 는 차원을 가진 양이라 로그 안에 들어가면 안 되는데, 어차피 $H$ 의 차이만 물리적 의미를 갖기 때문에 관행적으로 넘어간다. 찜찜하면 상대 엔트로피 형태로 쓰면 되고, 실제로 현대 문헌은 대부분 그렇게 쓴다. ↩

2. 이 이완이 얼마나 빠른가 하면, 평균 충돌 몇 번이면 초기 분포의 기억이 거의 지워진다. 그래서 연속체 CFD가 성립하는 것이고, 동시에 재진입체 앞의 충격층처럼 충돌이 몇 번밖에 안 일어나는 영역에서는 나비에-스토크스가 무너지는 것이다. “몇 번”이라는 아주 작은 수가 두 세계를 가른다. ↩

3. ELBM 논문이 처음 나왔을 때 LBM 커뮤니티의 반응이 미지근했던 이유가 재미있다. “안정성 문제는 격자를 늘리면 해결된다”는 반론이었는데, 격자를 늘릴 여유가 있으면 애초에 LBM을 쓸 이유의 절반이 사라진다는 점에서 자기모순에 가까웠다. 지금은 고 레이놀즈수 LBM에서 어떤 형태로든 엔트로피성 안정화를 쓰는 것이 표준에 가깝다. ↩