구속 동역학

편집 역사 토론
분자동역학 수치해석 계산화학 마지막 수정: 2026-08-04 04:42:15

1. 개요[편집]

구속 동역학
Constrained Dynamics
분야분자동역학 × 수치적분
목적빠른 결합진동 제거 → 시간 스텝 확대
수학적 정체홀로노믹 구속이 붙은 지표 3 미분대수방정식
대표 알고리즘SHAKE, RATTLE, LINCS, SETTLE
대안수소 질량 재분배(HMR)

구속 동역학(constrained dynamics)은 원자 사이의 결합 길이나 각도처럼 계에 부과된 홀로노믹 구속 σk(r)=0\sigma_k(\mathbf{r}) = 0 을 매 적분 스텝마다 정확히 만족시키도록 구속력을 함께 풀어 운동방정식을 적분하는 기법이다. 분자동역학에서는 사실상 “C–H 결합을 강체 막대로 굳혀 시간 스텝을 두 배로 늘리는 기술”의 다른 이름이다.

동기는 시간 스텝의 저주다. C–H 신축 진동은 파수 약 3000 cm⁻¹, 주기로 환산하면 11 fs 남짓이다. 베를레 적분이 이 진동을 안정적으로 따라가려면 주기의 1/10 이하, 즉 1 fs를 써야 한다. 그런데 우리가 보고 싶은 단백질 접힘이나 리간드 결합은 μs~ms 스케일이다. 아무도 관심 없는 수소 떨림 때문에 전체 계산 속도가 결정되는 이 상황이 견딜 수 없어서, 사람들은 그 자유도를 아예 없애기로 했다. 결합을 전부 고정하면 다음으로 빠른 운동은 수소가 낀 각도 굽힘(~1500 cm⁻¹, 주기 22 fs)이 되고, 스텝을 2 fs로 올릴 수 있다.1

2. 홀로노믹 구속과 라그랑주 승수[편집]

결합 kk 가 원자 i,ji,j 를 길이 dkd_k 로 묶는다면 구속식은 이렇게 쓴다.

σk(r)=∣ri−rj∣2−dk2=0\sigma_k(\mathbf{r}) = |\mathbf{r}_i - \mathbf{r}_j|^2 - d_k^2 = 0

제곱으로 쓰는 건 제곱근의 미분을 피하려는 실용적 선택이다. 이제 라그랑주 승수법으로 각 구속에 승수 λk\lambda_k 를 달면 운동방정식이 다음이 된다.

mir¨i=Fi−∑kλk∇iσk(r)m_i \ddot{\mathbf{r}}_i = \mathbf{F}_i - \sum_k \lambda_k \nabla_i \sigma_k(\mathbf{r})

승수 λk\lambda_k 가 곧 그 결합이 감당하는 장력이다. 미분방정식과 대수 구속식이 섞였으므로 이 계는 미분대수방정식이고, 위치 구속만 놓고 보면 지표 3이다. 구속식을 시간으로 한 번 미분한 숨은 구속

∇σk⋅r˙=0\nabla \sigma_k \cdot \dot{\mathbf{r}} = 0

까지 함께 강제하면 지표가 2로 내려간다. 아래 알고리즘들의 차이는 결국 “이 두 구속 중 무엇을 어떻게 강제하는가”로 요약된다.

3. SHAKE — 반복으로 밀어 넣기[편집]

SHAKE(Ryckaert–Ciccotti–Berendsen, 1977)는 가장 오래되고 가장 널리 이식된 방식이다. 발상은 무식할 정도로 단순하다. 일단 구속을 무시하고 베를레 한 스텝을 밟은 뒤, 위치를 조금씩 되밀어 모든 구속을 만족시킨다.

구속 kk 하나만 보면 보정량은 σk\sigma_k 에 대한 2차 방정식이 되는데, 2차항을 버리고 선형화하면 승수의 근사값이 닫힌 형태로 나온다. 문제는 한 원자가 여러 구속에 동시에 걸려 있어 하나를 고치면 다른 게 깨진다는 것. 그래서 SHAKE는 구속들을 가우스-자이델식으로 순회하며 모든 구속의 상대오차가 허용값(보통 10−4∼10−610^{-4} \sim 10^{-6}) 아래로 떨어질 때까지 반복한다.

장점은 일반성이다. 구속의 위상이 어떻게 얽혀 있든 그냥 돈다. 단점도 명확하다. 반복 횟수가 계에 따라 들쭉날쭉해 병렬화가 어렵고, 구속이 강하게 결합된 계에서는 수렴이 느리다. 결정적인 약점은 따로 있다 — SHAKE는 위치만 구속면에 올려놓고 속도는 내버려 둔다. 속도가 구속 다양체의 접평면을 조금씩 벗어나면 시간 가역성이 깨지고, 그토록 애써 얻은 심플렉틱 적분기의 장기 에너지 보존이 새기 시작한다.

4. RATTLE과 LINCS[편집]

RATTLE(Andersen, 1983)이 그 구멍을 막는다. 위치 갱신 후 위치 구속을 사영하고, 속도 갱신 후 숨은 구속 ∇σ⊤M−1p=0\nabla\sigma^\top M^{-1}\mathbf{p} = 0 을 두 번째 사영으로 한 번 더 강제한다. 그 결과 RATTLE은 velocity Verlet의 구속판에 정확히 대응하며, 구속 다양체 위에서 심플렉틱임이 증명되어 있다. 장시간 NVE 궤적을 돌릴 거면 SHAKE가 아니라 RATTLE이 정답이다. 이 계열의 기하학적 배경은 기하 수치적분에 정리돼 있다.

LINCS(Hess 외, 1997)는 접근이 다르다. 반복 대신 직접 해법을 쓴다. 구속 방향 성분을 사영으로 걷어내는 연산은 결합 행렬 BM−1B⊤\mathbf{B}M^{-1}\mathbf{B}^\top 의 역을 필요로 하는데, LINCS는 이 역행렬을 노이만 급수로 전개해 유한 항(기본 4차)에서 자른다. 결합이 약하게 얽혀 있으면 급수가 빨리 죽어서 몇 번의 희소 행렬-벡터 곱으로 끝난다. 사영 후에는 결합이 미세하게 짧아지는 기하학적 효과가 남으므로 회전 보정을 한 번 더 넣는다.

LINCS는 반복 횟수가 고정이라 결정론적이고 병렬화가 쉽다 — GROMACS가 기본값으로 쓰는 이유다. 대신 각도까지 구속해 구속들이 강하게 얽히면 급수 수렴이 나빠지므로 전개 차수를 올리거나 SHAKE로 돌아가야 한다.

5. SETTLE — 물은 그냥 손으로 푼다[편집]

생체 시뮬레이션 상자의 원자 대부분은 물이다. 그런데 TIP3P·SPC 같은 강체 3점 물 모형은 원자 3개에 구속 3개가 걸린 최악의 조합이라, 반복 해법도 급수 전개도 잘 안 먹힌다.

SETTLE(Miyamoto–Kollman, 1992)은 여기에 반복 대신 해석해를 준다. 구속을 무시한 갱신 위치의 무게중심을 그대로 두고, 원래 삼각형을 강체로 회전시켜 그 위치에 최적으로 맞추는 문제로 재해석하면 회전각들이 닫힌 형태로 떨어진다. 반복이 없으니 수렴 판정도, 허용오차도 없다. 정확하고 빠르다. 오늘날 물 분자는 SETTLE, 나머지 결합은 LINCS나 SHAKE로 나눠 처리하는 것이 표준 조합이다.2

6. 구속이 통계에 남기는 흔적[편집]

구속은 공짜가 아니다. 최소한 세 군데에 영향을 준다.

  • 자유도. 구속 NcN_c 개는 자유도를 그만큼 깎는다. Nf=3N−Nc−3N_f = 3N - N_c - 3 (병진 3 제외). 온도는 ⟨K⟩=12NfkBT\langle K\rangle = \frac{1}{2}N_f k_BT 로 읽으므로 이 수를 틀리면 온도가 체계적으로 어긋난다. 강체 물 100개짜리 테스트 계에서 구속을 안 빼면 30%씩 틀어진다는 얘기가 등분배 정리에 있다.
  • 압력. 구속력도 원자에 작용하는 힘이므로 비리얼 정리의 비리얼 합에 기여해야 한다. 이걸 빠뜨리면 온도는 멀쩡한데 압력만 크게 틀리며, 등온등압 앙상블 계산의 밀도가 조용히 어긋난다.
  • 배치 분포 자체. 이게 가장 미묘하다. “무한히 뻣뻣한 스프링”의 극한과 “처음부터 강체 구속”은 같은 분포를 주지 않는다. 구속계의 구성 적분에는 질량가중 구속 야코비안 Zkl=∑imi−1∇iσk⋅∇iσlZ_{kl} = \sum_i m_i^{-1}\nabla_i\sigma_k\cdot\nabla_i\sigma_l 의 행렬식이 붙고, 그 차이를 메우는 보정이 픽스만 퍼텐셜 −12kBTln⁡det⁡Z-\frac{1}{2}k_BT\ln\det \mathbf{Z} 다. 결합 길이만 고정하는 통상적 사용에서는 무시할 만하지만, 각도까지 구속하면 무시할 수 없다.3

7. 수소 질량 재분배라는 다른 길[편집]

같은 목표를 정반대 방향에서 공략하는 것이 수소 질량 재분배(HMR)다. 구속을 거는 대신, 무거운 원자의 질량 일부를 붙어 있는 수소로 옮겨 수소를 1 amu에서 3~4 amu로 살찌운다. 분자 전체 질량은 보존된다. 진동수는 k/μ\sqrt{k/\mu} 이므로 무거워진 수소의 운동이 느려지고, 결합 구속과 함께 쓰면 4 fs 스텝까지 열린다.

핵심은 질량이 구성 적분에 들어가지 않는다는 사실이다. 운동량 적분은 가우스 적분으로 분리되므로, 밀도·자유에너지·구조 같은 평형 물성은 질량을 바꿔도 원리적으로 불변이다. 대신 동역학은 왜곡된다 — 확산계수와 속도상관, 진동 스펙트럼은 믿으면 안 된다.4 자유에너지 계산에는 부담 없이 쓰고, 수송계수를 재는 계산에는 쓰지 않는 것이 상식선이다.

8. 관련 문서[편집]

9. Footnotes[편집]

  1. “결합을 굳혀도 되나?”라는 질문에 대한 변명은 이렇다. C–H 신축의 영점 에너지는 kBTk_BT 보다 훨씬 커서 상온에서 사실상 바닥 상태에 얼어 있고, 그러면 그 자유도는 애초에 고전 MD가 제대로 다루지 못하는 양자 자유도다. 잘못 다루느니 없애는 게 낫다는, 다소 뻔뻔하지만 실제로 잘 작동하는 논리. ↩

  2. SETTLE이 나오기 전에는 물 상자 하나 돌리는 데 SHAKE 반복이 전체 시간의 20%를 넘게 잡아먹는 일이 흔했다. 상자의 90%가 물인 계에서 물만 해석적으로 처리한 것은, 병목을 정확히 겨냥한 모범적 최적화다. ↩

  3. 그래서 “결합 길이는 구속해도 되지만 각도는 웬만하면 건드리지 말라”는 경험칙이 생겼다. 각도까지 굳히면 스텝은 좀 더 벌지만, 자유에너지 지형이 미묘하게 뒤틀린 채 며칠을 돌게 된다. 벌어들인 시간보다 비싼 대가다. ↩

  4. 수소를 네 배 무겁게 만들어 놓고 물의 확산계수를 보고했다가 리뷰에서 잡히는 사례가 주기적으로 반복된다. 평형 물성이 불변이라는 정리는 정확하지만, 그 정리에 “평형”이라는 단어가 들어 있다는 점을 다들 마지막에 읽는다. ↩