힘장

분자동역학 계산화학 마지막 수정: 2026-07-07 14:07:55

1. 개요[편집]

힘장 Force Field
분야	분자동역학 × 계산화학 × 통계역학
정체	원자 좌표 → 퍼텐셜 에너지 함수
대표 계열	AMBER, CHARMM, OPLS, GROMOS
쓰이는 곳	분자동역학, 도킹, 자유에너지 계산
대표 소프트웨어	GROMACS, LAMMPS, NAMD

양자역학은 정확하지만 느리고, 힘장은 근사지만 빠르다. 단백질 10만 원자를 마이크로초 동안 돌리려면 근사와 손을 잡아야 한다.

힘장(Force Field)은 원자들의 좌표를 입력받아 계의 퍼텐셜 에너지 $U(\mathbf{r}_1, \dots, \mathbf{r}_N)$ 를 돌려주는, 물리적으로 동기화된 해석적 함수의 집합과 그에 딸린 파라미터 전체를 가리킨다. 분자동역학에서 각 원자에 작용하는 힘은 이 퍼텐셜의 음의 기울기 $\mathbf{F}_i = -\nabla_i U$ 로 얻어지므로, 힘장이 곧 시뮬레이션의 물리 법칙 그 자체가 된다.

전자를 명시적으로 다루는 양자역학은 원리적으로 옳지만, 수만~수십만 원자를 나노초 이상 돌리기엔 턱없이 느리다. 힘장은 전자 자유도를 통째로 지워버리고 원자를 고전적 입자로 취급하는 대담한 근사를 통해, 뉴턴 운동방정식을 시간 적분할 수 있게 만든다. 정확도를 팔아 속도를 산 것이며, 이 거래가 계산생물학과 재료과학 전체를 떠받친다.¹

2. 퍼텐셜 함수의 구조[편집]

전형적인 생체분자 힘장의 퍼텐셜은 결합(bonded) 항과 비결합(non-bonded) 항으로 나뉜다.

U = \underbrace{\sum_{\text{bonds}} k_b (r - r_0)^2 + \sum_{\text{angles}} k_\theta (\theta - \theta_0)^2 + \sum_{\text{dihedrals}} \frac{V_n}{2}[1 + \cos(n\phi - \gamma)]}_{\text{결합 항}}

+ \underbrace{\sum_{i<j} \left[ 4\varepsilon_{ij}\left( \left(\frac{\sigma_{ij}}{r_{ij}}\right)^{12} - \left(\frac{\sigma_{ij}}{r_{ij}}\right)^{6} \right) + \frac{q_i q_j}{4\pi\varepsilon_0 r_{ij}} \right]}_{\text{비결합 항}}

각 항의 의미는 명료하다. 결합 신축과 각도 변형은 조화 진동자(용수철)로, 이면각(dihedral)은 회전 장벽을 주기 함수로 표현한다. 비결합 항의 앞부분이 그 유명한 레너드-존스 퍼텐셜로 반데르발스 인력과 반발을 담고, 뒷부분이 쿨롱 정전기 상호작용이다.

주목할 점은 결합을 용수철로 근사했다는 것이다. 이 때문에 표준 힘장에서는 화학 결합이 끊어지거나 생기지 못한다. 반응을 다루려면 ReaxFF 같은 반응성 힘장이 따로 필요하다.²

3. 파라미터화: 힘장의 진짜 알맹이[편집]

함수 형태보다 중요한 것이 파라미터다. $k_b$ , $r_0$ , $\varepsilon$ , $\sigma$ , 부분전하 $q_i$ 같은 수천 개의 숫자를 어떻게 정하느냐가 힘장의 품질을 결정한다. 파라미터는 두 가지 소스에서 온다.

양자화학 계산 — 밀도범함수이론이나 하트리-폭 방법으로 작은 분자 조각의 에너지 표면과 전하 분포를 계산해 맞춘다.
실험 데이터 — 밀도, 기화열, 용해도, 결정 구조 같은 거시 물성에 재현되도록 튜닝한다.

이 파라미터화 과정은 대단히 노동집약적이며, 힘장 개발 논문 한 편이 수년의 산물인 경우가 흔하다. 그래서 새 분자(특히 신약 후보나 새 리간드)에 기존 힘장을 확장하는 일 자체가 하나의 전문 분야다.

4. 주요 힘장 계열[편집]

생체분자 세계는 몇몇 힘장 “가문”으로 나뉘며, 각자 철학과 텃밭이 다르다.

힘장	태생	강점 영역	특징
AMBER	핵산·단백질	DNA/RNA, 단백질	부분전하를 QM 정전 퍼텐셜(RESP)로
CHARMM	생체막·단백질	지질, 막단백질	방대한 지질 파라미터, 개선 항 풍부
OPLS	유기 액체	소분자, 용매	액체 물성 재현에 특화
GROMOS	생체분자	자유에너지	일부 수소 통합(united-atom)

이들은 서로 호환되지 않는다. AMBER 파라미터를 CHARMM 물 모델과 섞으면 물리가 어긋난다. 물 모델(TIP3P, SPC/E 등)조차 힘장과 짝이 정해져 있어서, 조합을 함부로 바꾸면 결과가 미묘하게 틀어진다.

5. 해상도의 선택: 전원자에서 조립까지[편집]

같은 힘장 철학 안에서도 원자를 얼마나 세밀하게 표현하느냐는 또 하나의 결정이다. 이 선택이 다룰 수 있는 계의 크기와 시간 규모를 직접 좌우한다.

전원자(all-atom) — 수소를 포함한 모든 원자를 명시한다. 가장 정확하지만 그만큼 비싸다. 현대 단백질 시뮬레이션의 기본값.
통합원자(united-atom) — 메틸기( $\text{CH}_3$ ) 같은 곳의 수소를 무거운 원자에 흡수시켜 입자 수를 줄인다. 지질 이중층처럼 큰 계에서 여전히 유용하다.
거친 알갱이(coarse-grained) — 여러 원자를 하나의 비드로 뭉친다. MARTINI 힘장이 대표적이며, 세포막 규모의 마이크로초 시뮬레이션을 가능하게 한다.

정확도와 규모는 늘 상충한다. 세밀할수록 옳지만 느리고, 뭉칠수록 크게 볼 수 있지만 화학적 디테일을 잃는다. 무엇을 보려는지가 해상도를 결정한다.

6. 장거리 상호작용과 계산 비용[편집]

비결합 항이 계산 비용의 대부분을 차지한다. 원자 쌍이 $N$ 개면 소박하게 계산할 때 $O(N^2)$ 이 되어 대형 계에서는 감당이 안 된다. 그래서 두 가지 가속 기법이 표준이다.

레너드-존스 퍼텐셜처럼 빨리 감쇠하는 상호작용은 컷오프 반경 바깥을 무시하는 근거리 목록(cell list / neighbor list)으로 $O(N)$ 에 가깝게 만든다. 반면 $1/r$ 로 느리게 감쇠하는 정전기 항은 컷오프로 자르면 심각한 오차가 생기므로, 주기 경계에서 고속 푸리에 변환을 활용하는 Ewald 합산, 실무적으로는 PME(Particle Mesh Ewald)로 처리한다. 이 PME 덕분에 오늘날 GROMACS나 LAMMPS가 수백만 원자를 GPU에서 돌릴 수 있다.

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

“힘장 하나 잘못 골라서 논문 결과가 뒤집힌다”는 게 과장이 아니다. 같은 단백질도 AMBER냐 CHARMM이냐에 따라 접힘(folding) 경로가 다르게 나올 수 있다.
신약 개발에서 새 리간드의 파라미터를 못 구해 좌초하는 프로젝트가 흔하다. GAFF 같은 범용 힘장이 구원투수지만 정확도는 각오해야 한다.
물 모델은 힘장의 숨은 주인공이다. 물을 어떻게 모델링하느냐가 용매화 자유에너지 전체를 좌우한다.
분자동역학 결과를 발표할 때 힘장과 물 모델, 컷오프를 명시하지 않으면 재현이 불가능하다. 검증 및 확인의 기본 중 기본이다.

8. 관련 문서[편집]

대략적인 감을 잡자면, 밀도범함수이론으로 수백 원자를 한 스텝 미는 비용으로 힘장은 수십만 원자를 수천 스텝 밀 수 있다. 규모의 차이가 곧 다룰 수 있는 문제의 차이다. ↩
그래서 표준 힘장으로 “효소가 기질을 절단하는 순간”을 볼 수는 없다. 결합이 안 끊어지니까. 이럴 땐 QM/MM 하이브리드나 반응성 힘장을 꺼내야 한다. ↩

힘장