힘장

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

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

1. 개요[편집]

양자역학은 정확하지만 느리고, 힘장은 근사지만 빠르다. 단백질 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 하이브리드나 반응성 힘장을 꺼내야 한다. ↩

힘장