물 모형

분자동역학 계산화학 마지막 수정: 2026-08-18 04:37:26

1. 개요[편집]

물 모형 Water Model
분야	분자동역학 × 계산화학 × 통계역학
대표 계열	SPC/E · TIP3P · TIP4P/2005 · TIP5P
표준 구성	강체 기하 + 점전하 + O 자리 LJ 하나
필수 동반자	SETTLE + PME
업계 상수	모든 물성을 동시에 맞추는 모형은 없다

시뮬레이션에서 원자 수의 8할이 물인 경우가 흔하다. 즉 계산 시간의 8할과 오차의 8할도 물이다.

물 모형(water model)은 분자 시뮬레이션에서 물 분자 하나를 몇 개의 상호작용 자리(site)와 어떤 파라미터로 표현할지를 규정한 유효 퍼텐셜이며, 사실상 힘장의 부분집합이면서도 독립된 연구 분야를 이루고 있다. 표준형은 놀랄 만큼 단순하다 — 기하를 강체로 고정하고, 몇 개의 점전하를 얹고, 산소 자리에 레너드-존스 퍼텐셜 하나를 둔다. 이 소박한 물체가 사면체 그물을 짜고, 4 °C에서 밀도 최대를 내고, 단백질을 접히게 만든다.

물이 이렇게 따로 취급받는 이유는 물이 유난히 이상하기 때문이다. 밀도 최대, 비정상적으로 높은 끓는점과 비열, 음의 열팽창 구간, 압력에 따라 점도가 내려가는 구간까지 — 이 이상성 목록의 뿌리는 전부 수소결합의 방향성이고, 그것을 점전하 서너 개로 흉내 내는 것이 물 모형의 기예다.

2. 강체 점전하 모형의 지형[편집]

2.1. 3점 모형[편집]

산소 하나와 수소 둘, 자리가 곧 원자다. 음전하는 산소 위에 앉는다.

SPC(Simple Point Charge, 1981) — H–O–H 각도를 실험값이 아니라 정사면체 각 109.47°로 놓았다. 액체 구조를 맞추려는 의도적 선택이다.
SPC/E(Extended, 1987) — SPC의 전하를 키워 쌍극자 모멘트를 올리고, 대신 자체 분극 에너지를 상수로 빼주는 보정을 도입했다. 응집상 물 분자는 기체상보다 훨씬 큰 유효 쌍극자를 갖는데, 고정전하 모형은 그 분극 비용을 지불하지 않으므로 손으로 빼준 것이다. 확산계수와 밀도가 SPC보다 확연히 개선됐다.
TIP3P(1983) — 기하는 기체상 실험값(O–H 0.9572 Å, 각도 104.52°)을 쓴다. 오늘날 가장 널리 쓰이는 물 모형이며, 뒤에서 보듯 그 이유는 성능이 아니라 역사다.

2.2. 4점 모형과 M 자리[편집]

TIP4P 계열은 음전하를 산소에서 떼어내 H–O–H 이등분선을 따라 수소 쪽으로 약 0.15 Å 옮긴 가상 자리(M-site, 질량도 LJ도 없는 유령 자리)에 놓는다. 왜 이런 짓을 하는가?

세 점전하로 물의 정전기장을 표현할 때, 음전하를 산소에 두면 쌍극자는 맞출 수 있어도 사중극자 성분이 부족하다. 그런데 물의 사면체 배위를 만드는 것은 쌍극자가 아니라 그보다 높은 다중극이다(다중극 전개 참고). 음전하를 앞으로 당기면 같은 쌍극자를 유지하면서 사중극자를 키울 수 있고, 그 결과 산소-산소 동경분포함수의 첫 봉우리와 수소결합의 각도 분포가 눈에 띄게 개선된다. 자리 하나를 추가해서 산 것은 정전기장의 모양이지 전하량이 아니다.

부수적 이득도 있다. LJ 중심(산소)과 음전하 중심(M)이 분리되면서 반발과 인력의 기하학을 따로 조절할 수 있다.

TIP4P(1983) — 원형.
TIP4P-Ew(2004) — 원본 TIP4P는 구면 컷오프 조건에서 맞춰졌는데, 사람들이 에발트 합산/PME로 돌리기 시작하면서 물성이 어긋났다. 그래서 에발트를 전제로 다시 맞춘 판본이다.
TIP4P/2005 — 밀도, 밀도 최대점, 얼음 상들의 상대 안정성 등 넓은 물성 집합을 목표로 재파라미터화했다. 강체 비분극 모형 중 종합 성능이 가장 좋다는 평가가 굳어져 있다.
TIP4P/Ice — 얼음의 녹는점과 상평형을 맞추는 데 특화. 대신 상온 액체 물성은 양보한다.
OPC(2014) — 물의 실험 다중극 모멘트를 우선 목표로 잡아 4점 틀에서 재최적화한 계열.

2.3. 5점 모형[편집]

TIP5P는 여기에 산소의 비공유 전자쌍 방향으로 음전하 자리 둘을 더 둔다. 받개 방향성이 훨씬 선명해져서, 밀도 최대점을 상온 근처에서 제대로 재현하는 데 성공한 첫 사례로 유명하다. 대가는 비용이다. 전하 자리가 4개가 되면 물-물 쌍마다 정전기 계산이 $4\times4=16$ 쌍으로, 3점·4점의 9쌍 대비 약 1.8배가 된다. 물이 계의 대부분인 시뮬레이션에서 이 배율은 그대로 벽시계 시간이다.¹

3. 무엇을 맞추고 무엇을 틀리는가[편집]

여기가 이 문서의 핵심이다. 강체 비분극 모형은 모든 물성을 동시에 맞출 수 없다. 이것은 파라미터 튜닝을 게을리해서가 아니라 함수 형태의 한계에서 오는 구조적 제약이다.

모형	녹는점	밀도 최대점	확산계수	유전상수
실험	273 K	277 K (4 °C)	2.3	78.4
TIP3P	약 146 K	상온보다 한참 아래	약 5.5 (2배 이상 빠름)	약 90대 (과대)
SPC/E	약 215 K	약 240 K	약 2.5	약 70
TIP4P/2005	약 252 K	약 278 K	약 2.1	약 60 (과소)
TIP5P	약 274 K	약 277 K	약 2.6	약 80

(확산계수 단위는 $10^{-5}\,\mathrm{cm^2/s}$ , 298 K 기준. 문헌값은 컷오프·정전기 처리·계 크기에 따라 흔들리므로 자릿수와 경향만 읽는 것이 안전하다.)

표에서 읽어야 할 것은 개별 숫자가 아니라 엇갈림이다.

TIP5P는 밀도 최대점과 녹는점을 훌륭하게 잡지만 비싸고, 압력·온도 범위를 넓히면 무너진다.
TIP4P/2005는 밀도·상평형·표면장력·점도를 폭넓게 잘 맞추는데, 유전상수는 실험의 3/4 수준이다. 이온 용매화나 전기장 응답을 다룬다면 이 결손이 치명적일 수 있다.
TIP3P는 유전상수를 과대평가하고 확산이 2배 이상 빠르다. 즉 동역학이 통째로 빠르다. TIP3P 물속 단백질의 완화 시간을 실험과 비교할 때 이 사실을 잊으면 안 된다.

특히 녹는점이 하나같이 실험보다 낮다는 점을 눈여겨보자. 상온에서 맞춘 모형들이 상온 근처의 액체이기만 하면 된다는 조건으로 최적화됐기 때문이며, 그래서 “0 °C에서 언다”를 기대하고 냉각 시뮬레이션을 돌리면 낭패를 본다. 얼음을 다뤄야 하면 TIP4P/Ice처럼 그 목적으로 만든 모형을 골라야 한다.

3.1. 왜 구조적으로 불가능한가[편집]

세 가지 한계가 겹친다.

고정전하 — 액체 물 분자의 유효 쌍극자는 기체상 1.85 D보다 크지만 그 값이 환경에 따라 변한다. 고정전하 모형은 평균값 하나를 박아넣으므로, 기체상 이량체와 액체를 동시에 맞출 수 없다. SPC/E의 자체 분극 보정은 액체를 택하고 기체를 포기하겠다는 선언이었다.
비가법성 부재 — 수소결합의 협동성(사슬로 이어질수록 강해지는 성질)은 원리적으로 다체 효과인데, 이체 퍼텐셜은 그것을 담을 그릇이 없다. 응집상 평균값으로 이체 파라미터에 미리 녹여 넣는 것이 최선이다.
강체 + 고전역학 — 진동을 얼려버렸으니 적외선 스펙트럼을 낼 수 없고, 수소의 영점 운동과 터널링도 없다. 동위원소 효과( $\mathrm{D_2O}$ 의 녹는점이 3.8 °C로 올라가는 것 등)는 원리적으로 재현 불가다.

4. 그 위의 세대[편집]

4.1. 분극 모형[편집]

전하가 환경에 반응하게 만들면 위 한계 1과 2가 완화된다. 두 가지 구현 방식이 표준이다.

드루드 진동자(Drude oscillator) — 무거운 원자에 용수철로 매달린 가상의 음전하 입자를 붙인다. 외부 전기장이 걸리면 입자가 밀려나 유도 쌍극자가 생긴다. 물 쪽 대표가 SWM4-NDP(4점 + 음의 드루드 입자)이며, CHARMM 계열의 분극 힘장과 짝을 이룬다. 확장 라그랑주 방식으로 드루드 입자에 작은 질량을 주고 저온 열욕에 묶어 함께 적분하는 것이 실무 관행이다.
유도 쌍극자 + 다중극 — AMOEBA가 대표. 점전하 대신 원자별 영구 다중극(전하·쌍극자·사중극자)을 두고 상호 유도 쌍극자를 자기무모순으로 푼다. 정전기장의 표현력이 근본적으로 높지만, 매 스텝 유도 쌍극자 SCF 반복이 들어가 비용이 크게 오른다.

성능은 확실히 좋아진다. 특히 이온 주변의 물, 계면, 강한 전기장 아래에서 고정전하 모형이 실패하는 지점을 메운다. 그럼에도 널리 퍼지지 않은 이유는 비용과 생태계다. 단백질 힘장 전체가 함께 분극형이어야 의미가 있는데, 그 파라미터화가 아직 고정전하 힘장만큼 성숙하지 않았다. 분극성 힘장 참고.

4.2. 다체·기계학습 퍼텐셜[편집]

최근의 방향은 아예 참조 데이터를 결합 클러스터 방법 수준으로 만들고 그걸 학습시키는 것이다. MB-pol은 이체·삼체 항을 CCSD(T) 계산에 맞춰 적합하고 그 이상은 고전 분극 모형으로 처리하는 다체 전개 방식으로, 이량체·삼량체 에너지부터 액체 구조, 진동 스펙트럼, 밀도 이상까지 한꺼번에 재현한다는 점에서 판을 바꿨다는 평가를 받는다. 신경망 기반 기계학습 퍼텐셜도 같은 노선이며, 비용은 고전 모형보다 한두 자릿수 비싸지만 DFT보다는 훨씬 싸다.

여기까지 오면 남는 문제는 핵 양자 효과다. 참조 전자구조가 아무리 정확해도 수소를 고전 입자로 굴리면 오차가 남는다. 그래서 정밀 계산에는 경로적분 분자동역학이 동원되고, q-TIP4P/F처럼 애초에 경로적분과 함께 쓰도록 설계된 유연 모형도 따로 있다.

5. 실행 규약 — 짝을 맞춰라[편집]

물 모형을 쓴다는 것은 파라미터 파일 하나를 고르는 일이 아니라 시뮬레이션 설정 전체를 함께 고정하는 일이다.

강체 구속은 필수. 물의 O–H 신축은 파수 3000 cm⁻¹대라 이걸 살려두면 시간 스텝이 0.5 fs 수준으로 묶인다. 구속 동역학의 SHAKE/RATTLE로 굳힐 수도 있지만, 물처럼 구속 3개짜리 강체 삼각형은 해석적으로 정확한 해가 존재한다. 그것이 SETTLE이며, 반복 없이 닫힌 형태로 풀리므로 SHAKE보다 빠르고 수렴 실패도 없다. GROMACS·NAMD 등이 물에는 기본으로 SETTLE을 쓰는 이유다.
장거리 정전기는 에발트/PME로. 컷오프로 자르면 물의 유전 응답과 계면 구조가 눈에 띄게 왜곡된다. 더 중요한 것은 모형이 어떤 정전기 처리 아래에서 맞춰졌는지 일치시키는 것이다. TIP4P-Ew가 따로 만들어진 이유가 정확히 이것이고, 반대로 컷오프 시절에 맞춰진 모형을 PME로 돌리면 원 논문의 물성이 그대로 나오지 않는다. 주기경계조건 아래 상자 크기도 함께 보고해야 한다.
컷오프·꼬리 보정도 규약의 일부. LJ 컷오프 거리와 분산 꼬리 보정 여부에 따라 밀도가 소수점 셋째 자리에서 흔들린다. 반데르발스 힘 문서에서 지적했듯 계면이 있는 계에서는 균일 밀도 가정이 깨져 더 예민해진다.
힘장과의 궁합. AMBER의 단백질 파라미터는 TIP3P와, CHARMM은 수정된 TIP3P와, OPLS는 TIP4P/SPC 계열과 함께 검증됐다. 물만 좋은 걸로 바꾸면 단백질 거동이 함께 바뀌며, 그 조합은 아무도 검증한 적이 없다.

6. TIP3P는 왜 아직도 살아 있는가[편집]

객관적 성능만 보면 TIP3P는 밀리는 모형이다. 그런데 여전히 가장 많이 쓰인다. 이유는 셋이다.

역사적 종속성. AMBER와 CHARMM의 생체분자 파라미터가 1980~90년대에 TIP3P 물속에서 맞춰졌다. 단백질 힘장의 부분전하와 이면각 항이 그 용매 환경을 전제하고 있으므로, 물만 갈아 끼우면 균형이 깨진다.
비교 가능성. 30년치 문헌이 TIP3P로 쌓여 있다. 새 결과를 그 축적과 비교하려면 같은 물을 써야 한다.
싸다. 3점, 자리 간 상호작용 9쌍. 빠른 확산은 심지어 의도치 않은 이득이 되기도 한다 — 구조 표본을 더 빨리 훑기 때문에 접힘이나 결합 자유에너지 표본화가 빨라진다. 물론 시간 척도를 실험과 직접 비교하면 안 된다는 조건이 붙는다.

즉 TIP3P의 생존은 물리적 우수성이 아니라 **생태계 잠김(lock-in)**이다. 비디오 규격 전쟁이나 자판 배열과 같은 종류의 이야기이며, 실제로 근래에는 AMBER 계열도 OPC나 TIP4P-Ew로 옮겨가는 흐름이 뚜렷하다.²

7. 모형 고르기[편집]

단백질 구조·동역학이 목적이면 기존 힘장과 검증된 조합(흔히 TIP3P 또는 TIP4P-Ew)을 따르고, 물 자체의 물성이 목적이면 TIP4P/2005가 기본값이며, 얼음·저온이면 TIP4P/Ice로 간다. 이온 용매화나 전기장 응답이 걸려 있으면 고정전하 모형의 유전상수 결손부터 확인하고 필요하면 분극 모형으로 넘어가야 하고, 진동 스펙트럼이나 양성자 전달은 강체 모형이 애초에 자격이 없다.

무엇을 고르든 논문에 모형 이름, 정전기 처리, 컷오프, 구속 알고리즘, 온도·압력 조절기를 함께 적는 것이 검증 및 확인의 최소 요건이다. “물은 SPC 썼습니다” 한 줄로는 아무도 재현할 수 없다.³

8. 관련 문서[편집]

가상 자리(M, L)는 질량이 없어서 운동방정식에 직접 들어가지 않는다. 그래서 매 스텝 가상 자리에 걸린 힘을 실제 원자들에 기하학적으로 재분배해야 하는데, 이 재분배를 틀리게 구현하면 에너지가 조용히 새어나간다. 물 모형 구현 버그의 고전적 산지다. ↩
“왜 아직도 TIP3P냐”는 질문은 이 바닥의 영원한 떡밥이고, 답은 대개 “그럼 우리 힘장 전체를 다시 파라미터화해줄래?”로 끝난다. 물 하나 바꾸는 데 수년치 검증이 딸려오는 구조라, 개선을 알면서도 못 옮기는 전형적인 기술 부채다. ↩
실제로 같은 계·같은 모형인데 다른 그룹의 밀도가 0.5% 어긋나는 일이 벌어지고, 추적해보면 대개 LJ 꼬리 보정을 한쪽만 켰거나 PME 격자 간격이 달랐다. 물성 재현 논문에서 설정 표가 본문보다 긴 이유가 있다. ↩