주기경계조건

분자동역학 계산물리 마지막 수정: 2026-07-11 04:21:11

1. 개요[편집]

주기경계조건(Periodic Boundary Conditions, PBC)은 유한한 개수의 입자만 담을 수 있는 시뮬레이션 상자를 사방으로 무한히 복제하여, 마치 무한히 넓은 계(bulk)를 다루는 것처럼 근사하는 경계 처리 기법이다. 분자동역학과 밀도범함수이론 계산의 사실상 표준 경계조건으로, 상자의 한쪽 면으로 빠져나간 입자가 반대쪽 면으로 다시 들어오게 만들어 벽(wall)이라는 인위적 경계를 없앤다.

핵심 동기는 단순하다. 컴퓨터가 다룰 수 있는 원자는 아무리 많아야 수백만~수억 개인데, 이는 각설탕 하나에 든 원자 수( $\sim 10^{22}$ )에 비하면 먼지 한 톨이다. 이 먼지를 그냥 진공에 띄워두면 표면에 노출된 원자 비율이 압도적으로 커져서¹, 우리가 알고 싶은 “덩어리 내부”의 물성이 표면 효과에 완전히 묻혀버린다. PBC는 이 표면을 통째로 없애버리는 마법이다.

2. 무한 복제라는 발상[편집]

PBC의 그림은 이렇다. 한 변의 길이가 $L$ 인 정육면체 상자(primary cell)를 생각하고, 이 상자를 $x$ , $y$ , $z$ 세 방향으로 빈틈없이 타일링한다. 각 복제 상자(image cell)에는 원본과 완벽히 똑같은 입자들이, 똑같은 상대 위치에, 똑같은 속도로 존재한다. 원본 상자 안의 입자 $i$ 가 위치 $\mathbf{r}_i$ 에 있다면, 그 이미지들은 다음 위치에 무한히 늘어서 있다.

\mathbf{r}_i + \mathbf{n} L, \quad \mathbf{n} = (n_x, n_y, n_z) \in \mathbb{Z}^3

어떤 입자가 상자 오른쪽 벽을 넘어가면, 그 순간 왼쪽 벽으로 동일한 입자가 들어온다. 그래서 상자 안의 입자 수 $N$ 은 항상 보존된다. 위상수학적으로 이 공간은 3차원 토러스(torus)와 같다 — 팩맨이 화면 오른쪽으로 사라지면 왼쪽에서 나타나는 바로 그 세계다.² 좌표를 상자 안으로 되돌리는 연산(wrapping)은 보통 다음처럼 구현한다.

x_i \leftarrow x_i - L \left\lfloor \frac{x_i}{L} \right\rfloor

3. 최소 이미지 규약[편집]

무한 복제를 곧이곧대로 받아들이면 곤란한 질문이 하나 생긴다. 입자 $i$ 와 $j$ 사이의 거리를 잴 때, $j$ 의 원본을 봐야 하나 아니면 그 수많은 이미지 중 하나를 봐야 하나?

여기서 등장하는 것이 최소 이미지 규약(minimum image convention)이다. 규칙은 간단하다. 입자 $i$ 는 상대 입자 $j$ 의 무한한 이미지들 중 가장 가까운 하나하고만 상호작용한다. 성분별로 쓰면 이렇게 된다.

\Delta x = \Delta x - L \, \mathrm{round}\!\left( \frac{\Delta x}{L} \right)

이 한 줄이 PBC의 심장이다. 덕분에 각 입자쌍의 거리는 항상 $[-L/2, L/2)$ 범위로 접힌다. 다만 이 규약이 성립하려면 컷오프 반경( $r_c$ )이 상자 절반보다 작아야 한다는 절대 조건이 붙는다.

r_c < \frac{L}{2}

이걸 어기면 한 입자가 다른 입자의 원본과 이미지 두 개를 동시에 “보는” 참사가 벌어진다 — 자기 자신과 상호작용하는 원자가 튀어나오는 셈이다.³ 레너드-존스 퍼텐셜처럼 짧은 거리에서 급격히 감쇠하는 힘장은 컷오프가 자연스럽지만, 쿨롱 같은 장거리 힘은 컷오프로 끊으면 안 되기 때문에 에발트 합산 같은 별도 처방이 필요하다.

4. 표면 효과 제거와 유한 크기 오차[편집]

PBC의 최대 미덕은 표면을 지운다는 것이다. 상자 벽이 없으니 벽에 붙은 원자도 없고, 모든 원자는 동등하게 “내부” 원자처럼 취급된다. 덕분에 밀도, 압력, 확산계수 같은 벌크 물성을 표면 오염 없이 측정할 수 있다.

하지만 공짜 점심은 없다. PBC는 표면 효과를 유한 크기 효과(finite-size effect)로 바꿔치기할 뿐이다. 상자가 너무 작으면 다음과 같은 인위적 부작용이 생긴다.

자기상관 유령: 입자가 자기 자신의 이미지와 상호작용하면서, 실제로는 없는 주기적 질서가 강제된다. 파장이 $L$ 보다 긴 요동(long-wavelength fluctuation)은 아예 담기지 못한다.
상전이 왜곡: 임계점 근처의 발산하는 상관 길이(correlation length)가 상자 크기에 잘려, 상전이 온도나 임계 지수가 틀어진다.
유체 동역학 결합: 확산계수는 상자 크기 $L$ 에 대해 $1/L$ 로 수렴하는 유명한 보정항(Yeh–Hummer 보정)을 가진다.

그래서 실무에서는 여러 상자 크기로 계산을 돌려 $L \to \infty$ 로 외삽하거나, 물성이 상자 크기에 둔감해질 때까지 상자를 키운다. “격자가 다 했다”의 분자 버전으로, 여기서는 “상자 크기가 다 했다”가 된다.

5. 변형과 특수한 상자[편집]

정육면체만 고집할 이유는 없다. 계의 대칭성과 목적에 따라 다양한 상자가 쓰인다.

직육면체 / 사방정계 상자: 세 변의 길이가 다른( $L_x \neq L_y \neq L_z$ ) 가장 흔한 일반화. 계면이나 슬래브(slab) 계산에 유용하다.
삼사정계(triclinic) 상자: 결정의 단위격자를 그대로 반영하는 평행육면체. 밴드 구조 계산이나 응력 하 결정 시뮬레이션에서 필수.
절단 팔면체(truncated octahedron) / 마름모 십이면체: 구형 분자(예: 용매 속 단백질)를 감쌀 때 모서리 공간 낭비를 줄여 물 분자 수를 아낄 수 있는 절약형 상자.
슬래브 + 진공층: 표면·박막 계산에서는 일부러 한 방향에만 두꺼운 진공을 넣어 “2.5차원” 주기성을 만든다. 표면을 살리고 싶을 때 PBC를 반쯤 끄는 트릭이다.

한편 비평형 계산에서는 상자 자체를 흐르게 만드는 리스-에드워즈(Lees–Edwards) 경계조건으로 전단(shear) 유동을 구현하기도 한다. LAMMPS나 GROMACS 같은 주류 코드는 이 모든 상자 형태를 기본으로 지원한다.

6. 관련 문서[편집]

원자 1000개를 정육면체로 쌓으면 한 변이 10개, 표면에 노출된 원자가 전체의 절반에 가깝다. “내부 물성”을 보려는 사람에게 이건 재앙이다. ↩
그래서 PBC를 처음 배우면 다들 “아, 팩맨 우주였구나” 하고 무릎을 친다. 위상수학적으로 완전히 정확한 비유라는 게 더 킹받는 지점. ↩
$r_c \geq L/2$ 인 상태로 돌린 시뮬레이션은 조용히 틀린 답을 뱉는다. 에너지 보존은 멀쩡해 보이는데 물성만 이상하게 나오는, 디버깅 지옥의 단골 원인이다. ↩