분극성 힘장

분자동역학 계산화학 마지막 수정: 2026-09-02 04:51:18

상위 문서: 힘장

1. 개요[편집]

분극성 힘장 Polarizable Force Field
해결하려는 것	고정전하 모형이 못 담는 환경 의존 분극
세 갈래	유도쌍극자 · 드루드 진동자 · 변동전하
감쇠	Thole 감쇠(분극 파탄 방지)
대표 힘장	AMOEBA · CHARMM Drude · AMBER ff02
비용	고정전하 대비 대략 2~10배
효과가 큰 곳	이온 용매화 · 막 계면 · 양이온-π

물 분자의 쌍극자 모멘트는 기체에서 1.85 D, 액체에서 대략 2.6~3.0 D다. 고정전하 힘장은 이 둘 중 하나만 고른다.

분극성 힘장(polarizable force field)은 각 원자의 전하 분포가 주변 전기장에 반응해 변하도록 허용한 분자 힘장이다. 일반적인 힘장은 원자마다 부분전하 $q_i$ 를 못박아 두는데, 이 값은 애초에 특정 환경(대개 액체 물)의 평균적 분극 상태를 미리 반영해 튜닝된 값이다. 그래서 계산 대상이 그 환경 안에 있는 동안은 잘 맞지만, 분자가 환경을 갈아타는 순간 틀린다 — 수용액에서 지질 이중층 내부로 들어가는 이온, 벌크 물에서 단백질 소수성 주머니로 들어가는 리간드, 기체상에서 결정으로 가는 분자가 전부 그렇다.

이 문서는 분극을 어떻게 모형에 넣느냐만 다룬다. 결합·각·이면각·레너드-존스 퍼텐셜 같은 나머지 항, 파라미터화 철학, 힘장 계열의 텃밭 이야기는 힘장 문서에 있으니 반복하지 않는다.

2. 고정전하가 놓치는 것[편집]

고정전하 모형은 분극을 평균장으로 흡수한다. TIP3P 물의 쌍극자가 2.35 D인 것은 실측이 아니라 “액체 물의 밀도와 기화열을 맞추려면 이 정도로 과분극시켜야 한다”는 역산 결과다. 이 전략의 대가는 세 가지로 나타난다.

전이 가능성 부족: 물 안에서 맞춘 전하가 저유전 환경(막 내부, 단백질 코어, 기체상)에서 틀린다. 유전율이 80에서 2로 바뀌는데 전하는 그대로다.
다체 효과 부재: $\text{Na}^+$ 를 물이 여섯 개 둘러싸면 첫 번째 물과 여섯 번째 물의 결합 에너지가 다르다(협동적 감쇠, cooperativity). 쌍 상호작용의 합으로는 이 비가법성을 만들 수 없다.
이방성 응답 부재: 벤젠 고리는 면 수직 방향과 면 내 방향의 분극률이 다르다. 양이온-π 상호작용의 세기가 여기서 나온다.

특히 이온이 문제다. 고정전하 힘장에서 할로겐 음이온의 물-공기 계면 흡착이나 이온쌍 형성 자유에너지가 계통적으로 틀리는 것은 오래된 병이다. 실무적 미봉책으로 이온의 전하를 통째로 줄이는 전하 스케일링(electronic continuum correction, 대개 $1/\sqrt{\varepsilon_\infty}\approx 0.75$ 배)이 널리 쓰이는데, 이것 자체가 “분극을 평균적으로 흉내 내겠다”는 자백이다.

3. 유도쌍극자 모형[편집]

가장 직접적인 처방. 원자 $i$ 에 등방 분극률 $\alpha_i$ 를 주고, 그 자리의 전기장에 비례하는 유도쌍극자를 만든다.

\boldsymbol{\mu}_i = \alpha_i \mathbf{E}_i = \alpha_i\left(\mathbf{E}_i^{(0)} + \sum_{j \ne i} \mathbf{T}_{ij}\,\boldsymbol{\mu}_j\right)

$\mathbf{E}_i^{(0)}$ 은 영구 전하(또는 다중극 전개를 사각극까지 쓴 영구 다중극)가 만드는 장이고, $\mathbf{T}_{ij}$ 는 쌍극자장 텐서다. 유도쌍극자가 서로의 장에 다시 기여하므로 이 식은 자기무모순(self-consistent) 연립방정식이고, 이것이 분극성 힘장 비용의 정체다.

에플퀴스트(Applequist, 1972)가 원자 단위 분극률을 상호작용시키는 이 모형을 정식화했는데, 곧바로 유명한 병이 드러났다. 두 분극 중심이 가까워지면 상호 유도가 폭주해서 유도쌍극자가 발산하는 분극 파탄(polarization catastrophe)이다. 거리가 $r < (4\alpha_i\alpha_j)^{1/6}$ 정도로 줄면 실제로 터진다.

톨레 감쇠(Thole, 1981)가 표준 해법이다. 점 분극률 대신 유한한 폭으로 퍼진 전하 분포를 가정해 근거리에서 $\mathbf{T}_{ij}$ 를 부드럽게 죽인다. 감쇠 폭은 무차원 거리 $u = r_{ij}/(\alpha_i\alpha_j)^{1/6}$ 로 정의되고, 감쇠 세기 $a$ 는 분자 분극률 실험값에 맞춰 하나로 정한다(AMOEBA는 $a=0.39$ ). 파탄을 막을 뿐 아니라 분자 분극률 텐서의 이방성을 등방 원자 분극률의 조합만으로 재현하는 효과까지 있다는 것이 이 처방의 아름다운 점이다.

4. 드루드 진동자[편집]

두 번째 계보는 전자 구름을 용수철에 매달린 가상 입자로 치환한다. 원자 $i$ 에 전하 $q_D$ 를 갖는 질량 작은 드루드 입자를 붙이고 강성 $k_D$ 의 조화 용수철로 연결하면, 정적 극한에서 이 계의 분극률은

\alpha_i = \frac{q_D^2}{k_D}

이 되어 유도쌍극자 모형과 정확히 같은 물리를 준다. 다른 점은 구현이다. 드루드 입자는 그냥 하나의 추가 입자이므로, 기존 분자동역학 코드의 정전기 루틴(컷오프, 에발트 합산/PME, 셀 리스트)을 한 줄도 안 고치고 재사용할 수 있다. 그래서 쉘 모델(shell model), 전하-용수철 모델(charge-on-spring)이라는 다른 이름으로도 불리며 고체물리에서는 훨씬 오래된 물건이다.

문제는 시간 적분이다. 드루드 입자를 매 스텝 에너지 최소화하면(SCF) 결국 유도쌍극자와 같은 비용이 든다. 대신 쓰는 것이 확장 라그랑지안(extended Lagrangian) 전략이다. 드루드 입자에 아주 작은 질량(보통 0.4 amu)을 주고 실제 자유도처럼 적분하되, 드루드 자유도만 별도의 저온 열욕(약 1 K)에 묶어 실제 원자의 열운동과 에너지 교환이 일어나지 않게 한다. 그러면 드루드 입자는 SCF 해 근방에서 미세하게 진동하며 따라다니고, 계는 진짜 SCF 궤적과 거의 구별되지 않는다. 카-파리넬로 분자동역학이 전자 궤도함수에 쓰는 것과 정확히 같은 아이디어다.¹

라무뢰와 루(Lamoureux & Roux, 2003)가 이 이중 열욕 방식을 확립했고, SWM4-NDP 물 모형과 함께 CHARMM Drude 힘장으로 이어졌다. 실무적 대가는 두 가지다. 시간 스텝을 1 fs로 줄여야 하고, 드루드 온도가 슬금슬금 오르는 것을 감시해야 한다 — 드루드 열욕이 새면 분극이 물리적 의미를 잃고, 심하면 드루드 입자가 이웃 원자 쪽으로 떨어져 나가는 극화 붕괴가 난다.

5. 변동전하[편집]

세 번째 계보는 쌍극자를 도입하는 대신 원자 전하 자체를 변수로 만든다. 전자음성도 균등화(electronegativity equalization) 원리에 기대어, 계 전체의 화학 퍼텐셜이 같아지도록 전하를 분자 내에서 흘려보낸다. 에너지를 전하에 대해 2차까지 전개하면

E(\{q\}) = \sum_i \left(\chi_i^0 q_i + \tfrac{1}{2}\eta_i^0 q_i^2\right) + \sum_{i<j} J_{ij}(r_{ij})\,q_i q_j

$\chi^0$ 는 전자음성도, $\eta^0$ 는 화학적 경도, $J_{ij}$ 는 쿨롱 상호작용이다. 전하 보존 제약 아래 최소화하면 선형계 하나로 전하가 결정된다. 라페와 고다드의 QEq(1991)가 대표 정식화이고, CHARMM 계열의 변동전하 힘장(Patel & Brooks, 2004)도 같은 뿌리다.

장점은 명확하다 — 추가 입자도, 쌍극자 텐서도 없고 선형계 하나만 풀면 된다. 단점도 명확하다.

면 밖 분극이 안 된다. 전하가 원자 위치에만 있으므로 평면 분자의 면 수직 방향 분극률이 구조적으로 0이다. 벤젠의 양이온-π를 못 다룬다.
장거리 전하 이동이 비물리적이다. 긴 사슬 분자에서 한쪽 끝을 건드리면 반대쪽 끝까지 전하가 흘러, 분자가 금속처럼 행동한다. 분극률이 길이에 초선형으로 증가하는 병으로 나타난다.

그래서 순수 변동전하 힘장은 생체분자 주류에서 밀렸고, 유도쌍극자와 결합한 하이브리드나 반응성 힘장(ReaxFF의 전하 평형 부분)에서 살아남았다. 장거리 이동 문제는 전하 이동을 결합 단위로 제한하는 분할전하(split-charge) 계열이나 ACKS2 같은 확장으로 보완한다.

6. 대표 힘장[편집]

힘장	분극 방식	정전기	비고
AMOEBA	Thole 감쇠 유도쌍극자	원자 다중극(사각극까지)	Tinker/OpenMM, 가장 비싸고 가장 정교
CHARMM Drude	드루드 진동자	점전하 + 이방성 보정	확장 라그랑지안, 생체막 강점
AMBER ff02	유도쌍극자	RESP 점전하	초기 세대, 후속 흐름은 다른 계열로 이동
CHARMM FQ	변동전하	동적 점전하	개념 검증 성격이 강함

AMOEBA(Ponder & Ren)가 가장 급진적이다. 영구 정전기를 점전하가 아니라 원자 중심 다중극으로 표현하고(단극+쌍극자+사각극), 여기에 톨레 감쇠 유도쌍극자를 얹는다. 그래서 수소결합의 방향성이나 할로겐 결합의 시그마 홀 같은 이방성 현상이 임의의 가상 자리(virtual site) 없이 자연스럽게 나온다. 대가는 비용이다 — 다중극 상호작용은 점전하보다 항이 훨씬 많고, PME도 다중극 버전을 따로 구현해야 한다. AMBER ff02(Cieplak, Caldwell, Kollman, 2001)는 기존 점전하 위에 유도쌍극자를 얹은 보수적 설계로 일찍 나왔지만, 파라미터 생태계가 뒤따르지 못해 널리 쓰이지는 못했다.

7. SCF의 비용과 함정[편집]

유도쌍극자 계열의 실행 시간을 결정하는 것은 자기무모순 반복이다. 보통 야코비/가우스-자이델 계열이나 전처리 켤레기울기법으로 푸는데, 수렴 판정을 RMS 쌍극자 변화 $10^{-5}$ D 정도로 두면 스텝당 4~8회 반복이 든다. 반복 한 번이 사실상 정전기 계산 한 번이므로, 분극성 계산이 고정전하보다 비싼 이유의 거의 전부가 여기 있다.

여기 아주 고약한 함정이 하나 있다. SCF를 느슨하게 수렴시키면 힘이 보존력이 아니게 된다. 유도쌍극자에 대한 에너지의 변분 조건이 정확히 만족되어야 헬만-파인만식 힘 공식이 성립하는데, 반복을 일찍 끊으면 그 조건이 깨져서 힘이 에너지의 정확한 기울기가 아니게 된다. 결과는 꾸준한 에너지 드리프트로 나타나고, NVE 검증을 안 하고 랑주뱅 동역학이나 온도조절기로 돌리면 열욕이 드리프트를 흡수해 버려서 문제를 알아채지도 못한다. 분극성 MD에서 “일단 돌려”가 특히 위험한 지점.²

우회로가 몇 가지 있다.

확장 라그랑지안: 드루드 방식처럼 쌍극자에 가상 질량을 주고 저온 열욕에 묶는다. 스텝당 SCF가 사라진다.
예측자-교정자(ASPC 등): 이전 스텝들의 쌍극자를 외삽해 초기값을 잘 잡아 반복 수를 절반 이하로 줄인다. 시간 가역성을 유지하도록 설계된 버전을 써야 드리프트가 안 난다.
절단 섭동 전개(OPT $n$ , iAMOEBA 등): SCF를 아예 포기하고 유도 방정식을 고정된 차수까지만 전개한다. 반복 수가 고정이므로 힘이 해석적으로 정의되고, 결과적으로 완전히 보존적이다. 정확도를 조금 내주고 재현성을 산다.

비용을 정리하면 대략 이렇다. 드루드는 고정전하 대비 2~~4배, AMOEBA는 5~~10배 수준이며, GPU 구현과 다중 시간 스텝으로 상당 부분을 되찾을 수 있다. 이 배수는 자유에너지 섭동 같은 샘플링 집약적 계산에서 그대로 실험 기간으로 환산된다.

8. 어디서 이득이 큰가[편집]

분극을 넣는다고 모든 것이 좋아지지는 않는다. 벌크 물에서 잘 튜닝된 고정전하 힘장은 여전히 강력하고, 비용 대비로 보면 지는 싸움인 경우가 많다. 이득이 명확한 곳은 환경의 유전 특성이 급격히 변하는 문제들이다.

이온 용매화와 이온쌍: 수화 자유에너지, 이온 선택성, 계면 흡착. 특히 $\text{Mg}^{2+}$ , $\text{Ca}^{2+}$ 같은 다가 양이온은 고정전하로는 배위 구조부터 틀린다.
이온 채널과 막 단백질: 유전율 80의 물에서 2의 지질 내부로 이온이 들어가는 과정 자체가 분극 문제다. 채널의 선택 필터에서 이온-카보닐 배위 에너지는 다체 효과가 지배한다.
양이온-π, 할로겐 결합: 방향성 있는 이방성 상호작용. AMOEBA 같은 다중극 힘장의 홈그라운드.
핵산과 다가 이온: DNA 인산 골격 주변의 이온 응축은 고정전하 모형에서 계통 오차가 크다.
결정 구조 예측·용매화 자유에너지: 기체상에서 응집상으로 넘어가는 에너지를 다루므로 전이 가능성이 곧 정확도다.

반대로 큰 단백질의 접힘 경로처럼 마이크로초~밀리초 샘플링이 병목인 문제에서는, 같은 계산 자원을 분극에 쓰는 것보다 샘플링에 쓰는 게 나은 경우가 흔하다. “정확한 힘장으로 짧게” 대 “덜 정확한 힘장으로 길게” 는 이 분야의 영원한 논쟁이고, 정답은 관심 있는 물리량이 무엇이냐에 달렸다.

9. 현업에서의 현실[편집]

파라미터 생태계가 압도적으로 얇다. 새 리간드 하나에 AMOEBA 파라미터를 만드는 일은 고정전하보다 몇 배 더 노동집약적이고, 자동화 도구가 있어도 검증은 여전히 손이다.
드루드는 시간 스텝 1 fs가 사실상 강제다. 수소 결합 구속(구속 동역학)으로 2 fs를 벌던 습관이 안 통한다.
분극성 힘장 결과를 발표할 때 SCF 수렴 기준과 드루드 열욕 설정을 안 쓰면 재현이 불가능하다. 검증 및 확인의 기본이며, 고정전하 시절의 “컷오프와 물 모형만 쓰면 된다”보다 항목이 늘었다.
GROMACS, OpenMM, LAMMPS, NAMD, Tinker 모두 어떤 형태로든 분극성을 지원하지만 지원 범위와 성능이 제각각이라, 힘장을 고르는 순간 코드가 반쯤 정해진다.

10. 관련 문서[편집]

그래서 드루드 MD를 “가난한 자의 카-파리넬로”라고 부르기도 한다. 전자 파동함수 대신 용수철에 매달린 가짜 입자를 쓰는 것뿐인데, 전자 자유도를 이온 자유도보다 훨씬 차갑게 유지해 단열 분리를 흉내 낸다는 발상은 완전히 같다. ↩
이 함정이 고약한 이유는 결과가 그럴듯해 보이기 때문이다. 구조 분포도 멀쩡하고 RDF도 그럴싸한데 에너지만 조용히 흐른다. NVE로 수십 ps 돌려 총에너지 기울기를 보는 30분짜리 검사가 논문 하나를 구할 수 있다. ↩