랑주뱅 동역학

분자동역학 계산물리 마지막 수정: 2026-07-20 04:18:07

1. 개요[편집]

랑주뱅 동역학 Langevin Dynamics
분야	분자동역학 × 통계역학 × 확률미분방정식
지배 식	랑주뱅 방정식 (마찰 + 무작위 힘)
핵심 정리	요동-소산 정리
역할	온도조절기 · 열욕 결합
대표 적분기	BAOAB

뉴턴 방정식에 마찰이라는 브레이크와 무작위 발길질을 더하면, 계는 저절로 목표 온도로 수렴한다.

랑주뱅 동역학(Langevin dynamics)은 뉴턴의 운동방정식에 마찰항과 무작위 힘을 더해, 명시적으로 모델링하지 않은 주변 환경(열욕, heat bath)과의 상호작용을 확률적으로 모사하는 분자동역학 기법이다. 용매 분자 수만 개를 일일이 계산하는 대신, 그 효과를 “느리게 끌어당기는 마찰 + 마구잡이로 때리는 힘”이라는 두 항으로 요약해 버린다.¹

핵심은 이 두 항이 무관하지 않다는 것이다. 마찰은 계에서 에너지를 빼앗고, 무작위 힘은 에너지를 넣어준다. 둘의 균형이 정확히 맞도록 강도를 묶어두면, 계는 어떤 초기 상태에서 출발하든 목표 온도의 정준 앙상블(NVT)로 저절로 이완된다. 그래서 랑주뱅 동역학은 그 자체로 훌륭한 서모스탯이 된다.

2. 랑주뱅 방정식[편집]

입자 $i$ 에 대한 지배 방정식은 뉴턴 식에 두 항이 붙은 확률미분방정식이다:

m_i \frac{d^2 \mathbf{r}_i}{dt^2} = -\nabla_i U - \gamma m_i \frac{d\mathbf{r}_i}{dt} + \sqrt{2 \gamma m_i k_B T}\,\boldsymbol{\eta}_i(t)

우변 첫 항은 평소처럼 힘장에서 오는 보존력 $-\nabla_i U$ 다. 둘째 항 $-\gamma m_i \dot{\mathbf{r}}_i$ 는 속도에 비례해 저항하는 마찰력으로, $\gamma$ 는 마찰 계수(단위: 1/시간)다. 셋째 항이 무작위 힘으로, $\boldsymbol{\eta}_i(t)$ 는 평균 0, 델타 상관을 갖는 백색 잡음이다:

\langle \boldsymbol{\eta}_i(t) \rangle = 0, \qquad \langle \eta_{i,\alpha}(t)\,\eta_{j,\beta}(t') \rangle = \delta_{ij}\,\delta_{\alpha\beta}\,\delta(t - t')

$\gamma$ 가 0이면 순수 뉴턴 역학(베를레 적분으로 푸는 NVE 앙상블)으로 되돌아가고, $\gamma$ 가 크면 관성이 무의미해지는 과감쇠 극한, 즉 브라운 동역학으로 넘어간다.

3. 요동-소산 정리[편집]

무작위 힘의 세기 앞에 붙은 계수 $\sqrt{2\gamma m_i k_B T}$ 가 우연이 아니라는 점이 이 이론의 정수다. 마찰(소산, dissipation)의 세기 $\gamma$ 와 무작위 힘(요동, fluctuation)의 세기가 온도 $T$ 를 매개로 정확히 묶여 있어야만 계가 올바른 볼츠만 분포로 수렴한다. 이 관계가 요동-소산 정리(fluctuation-dissipation theorem)다.

직관은 이렇다. 마찰만 있으면 계는 에너지를 계속 잃고 0 K로 얼어붙는다. 무작위 힘만 있으면 계는 에너지를 계속 받아 무한히 뜨거워진다. 둘의 강도를 온도로 저울질해 딱 맞추면, 평균적으로 들어오는 에너지 = 나가는 에너지가 되어 계의 온도가 $T$ 에 붙박인다. 자연에서 열평형이 유지되는 원리 그 자체를 알고리즘으로 옮긴 셈이다.²

4. 브라운 동역학과의 대비[편집]

마찰 계수 $\gamma$ 를 아주 크게 잡으면 가속도 항 $m\ddot{\mathbf{r}}$ 을 무시할 수 있는 과감쇠(overdamped) 극한이 되고, 이때 남는 1차 방정식이 브라운 동역학(Brownian dynamics)이다:

\gamma m_i \frac{d\mathbf{r}_i}{dt} = -\nabla_i U + \sqrt{2 \gamma m_i k_B T}\,\boldsymbol{\eta}_i(t)

관성이 사라졌으니 속도를 추적할 필요 없이 위치만 갱신하면 된다. 콜로이드 입자, 고분자 사슬, 세포 내 확산처럼 물속에서 느리게 움직여 관성이 무의미한 계에 잘 맞고, 관성 시간 스케일을 건너뛰므로 시간 스텝을 크게 쓸 수 있다. 정리하면 뉴턴(관성 O, 마찰 X) → 랑주뱅(관성 O, 마찰 O) → 브라운(관성 X, 마찰 O) 순으로 감쇠가 세지는 스펙트럼이다.

5. 적분기: BAOAB[편집]

랑주뱅 방정식은 무작위 항 때문에 일반 베를레 적분을 그냥 못 쓴다. 그래서 방정식을 세 조각으로 쪼개 번갈아 적분하는 연산자 분할(operator splitting) 기법이 쓰인다. 세 조각은 관례적으로 이렇게 부른다.

A (drift) — 위치를 속도로 밀기: $\mathbf{r} \mathrel{+}= \tfrac{1}{2}\mathbf{v}\,\Delta t$
B (kick) — 속도를 힘으로 차기: $\mathbf{v} \mathrel{+}= \tfrac{1}{2}(\mathbf{F}/m)\,\Delta t$
O (Ornstein-Uhlenbeck) — 마찰·무작위 힘으로 속도를 열화: $\mathbf{v} \leftarrow c_1 \mathbf{v} + c_2\,\boldsymbol{\xi}$

이들을 대칭으로 배열한 BAOAB 순서가 특히 유명하다. 여러 분할 순서(BAOAB, OBABO, ABOBA 등) 중 BAOAB가 큰 시간 스텝에서도 위치 관련 통계량(구조·자유에너지)의 오차가 가장 작다는 것이 이론·실험으로 밝혀져, 오늘날 랑주뱅 서모스탯의 사실상 표준이 됐다.³ 마찰이 속도를 확률적으로 감쇠시키는 O 단계 덕에 랑주뱅 열욕은 노제-후버 같은 결정론적 서모스탯보다 국소적이고 잘 발산하지 않는다.

6. 활용[편집]

랑주뱅 동역학은 쓰임새가 넓다.

온도조절기 — NVT 앙상블을 뽑는 가장 견고한 방법 중 하나. 특히 초기 이완이나 담금질 과정에서 계를 빠르게 목표 온도로 붙잡는다. 담금질 모사와도 궁합이 좋다.
암시적 용매(implicit solvent) — 용매 분자를 명시적으로 두는 대신 마찰·무작위 힘으로 대체하면, 물 수만 개를 계산할 자원으로 용질 자체를 더 오래 돌릴 수 있다.
거친입자(coarse-grained) 모델 — 원자 여럿을 하나의 알갱이로 묶으면 잃어버린 자유도의 요동이 사라지는데, 랑주뱅 항이 그 빠진 자유도의 열적 효과를 되살려준다. 소산 입자 동역학(DPD)이 대표 사례.
샘플링 가속 — 위치 정보만 필요한 경우 큰 시간 스텝으로 상태 공간을 넓게 훑는 데 유리하다.

GROMACS, LAMMPS, OpenMM 등 주요 코드가 모두 랑주뱅(및 브라운) 적분기를 기본으로 제공한다. “일단 NVT로 안정화부터”라는 실무 국룰에서 가장 먼저 손이 가는 도구다.

7. 관련 문서[편집]

이름은 1908년 폴 랑주뱅(Paul Langevin)이 브라운 운동을 기술하려고 세운 방정식에서 왔다. 아인슈타인이 1905년 확산 계수로 브라운 운동을 설명했다면, 랑주뱅은 같은 결과를 “마찰 + 무작위 힘”이라는 훨씬 직관적인 운동방정식으로 다시 유도했다. 확률미분방정식의 역사적 원조 격. ↩
그래서 마찰 계수 $\gamma$ 는 물리적 정확도가 아니라 조절 손잡이에 가깝다. 값을 크게 하면 온도 제어는 빨라지지만 실제 동역학(확산계수·상관시간)이 왜곡된다. 확산계수처럼 진짜 동역학 물성을 뽑고 싶다면 $\gamma$ 를 작게 두거나 아예 NVE로 돌려야 한다 — 서모스탯은 구조를 얻는 도구지 동역학을 얻는 도구가 아니라는 격언이 여기서 나온다. ↩
BAOAB의 우수성은 Leimkuhler와 Matthews가 2013년 전후로 정리했다. 핵심은 O 단계를 A 단계들 사이 가운데에 끼워 넣은 대칭 구조가, 큰 시간 스텝에서 위치 분포의 유효 온도 오차를 상쇄한다는 것. 같은 재료(A·B·O)를 어떤 순서로 쌓느냐가 정확도를 가른다는 게 은근히 심오하다. ↩