확산 제한 응집

계산물리 시뮬레이션 마지막 수정: 2026-07-31 05:41:26

1. 개요[편집]

확산 제한 응집(diffusion-limited aggregation, DLA)은 무작위걷기하는 입자를 멀리서 하나씩 풀어놓고, 이미 자란 덩어리에 닿는 순간 그 자리에 붙여 버리는 성장 모형이다. 규칙은 이 두 줄이 전부인데, 결과물은 가지가 사방으로 뻗은 성긴 나뭇가지 구조 — 이른바 프랙탈 클러스터가 된다.¹

1981년 위튼(T. A. Witten)과 샌더(L. M. Sander)가 금속 연무 입자의 응집을 설명하려고 제안했고, 발표되자마자 전착·유전파괴·점성 손가락·광물 침착까지 전혀 다른 현상들이 같은 모양을 낸다는 것이 알려지면서 비평형 성장 연구의 표준 모형이 됐다. 성장이 표면 반응이 아니라 물질의 확산 공급에 의해 제한된다는 것이 이름의 뜻이다. 반대쪽 극단, 즉 확산이 충분히 빨라 표면 어디든 재료가 남아도는 반응 제한 성장은 조밀한 덩어리를 만든다. 같은 재료라도 어느 쪽이 병목이냐에 따라 형태가 통째로 갈린다는 것이 이 모형의 첫 번째 교훈이다.

2. 라플라스 성장[편집]

무작위걷기 입자의 도달 확률장 $\phi$ 는 정상상태에서 라플라스 방정식을 만족한다.

\nabla^2 \phi = 0, \qquad \phi|_{\text{클러스터}} = 0, \qquad v_n \propto -\left.\frac{\partial \phi}{\partial n}\right|_{\text{표면}}

경계가 자기 자신의 해에 의해 움직이는 이동경계 문제다. 그런데 이 틀은 DLA 전용이 아니다.

현상	퍼텐셜 $\phi$	경계 속도
DLA	걷기 입자 농도	확산 플럭스
점성 손가락(헬레-쇼)	압력 (다르시 법칙)	$v \propto -\nabla p$
전착	전해질 이온 농도	이온 플럭스
유전파괴	정전 퍼텐셜	$\propto \lvert\nabla\phi\rvert^{\eta}$

즉 새프먼-테일러 불안정에서 점성이 낮은 유체가 손가락을 내미는 것과, 절연체에서 리히텐베르크 무늬가 뻗는 것과, DLA 가지가 자라는 것은 같은 방정식의 서로 다른 실현이다. 유전파괴 모형(DBM)은 성장확률을 $\lvert\nabla\phi\rvert^{\eta}$ 로 일반화해 이 계열을 한 손잡이로 묶는데, $\eta = 1$ 이 정확히 DLA고 $\eta \to 0$ 은 표면 어디나 똑같이 자라는 조밀한 에덴 성장이 된다.

3. 조화측도와 차폐[편집]

표면의 한 점이 다음 입자를 받을 확률은 무한원점에서 출발한 무작위걷기가 그 점에 처음 닿을 확률, 즉 조화측도(harmonic measure)와 같다. 여기에 DLA의 모든 것이 들어 있다.

돌출한 팁은 등퍼텐셜선이 몰려 $\lvert\nabla\phi\rvert$ 가 크고, 가지 사이 협곡(fjord) 안쪽은 입자가 들어오다 벽에 먼저 붙어 버려 사실상 굶는다. 이것이 차폐(screening)다. 그래서 조금 앞선 팁은 더 빨리 자라고, 뒤처진 곳은 영영 뒤처진다 — 되먹임이 양(+)이라 구조가 계속 성기게 갈라진다. 표면장력 같은 안정화 항이 없으면 이 불안정성을 잡아 줄 길이 스케일 자체가 없고,² DLA에서는 그 역할을 격자 간격과 입자 크기(=잡음)가 대신한다.

그래서 성장확률 $p_i$ 는 표면 위에서 균등하기는커녕 자릿수 단위로 벌어진다. 이 불균일을 정량화하는 표준 도구가 모멘트 $\sum_i p_i^{q} \sim \varepsilon^{\tau(q)}$ 와 그 르장드르 변환인 $f(\alpha)$ 스펙트럼이다. $q$ 를 키우면 팁만, 줄이면 협곡만 골라 보는 셈인데, DLA에서는 한쪽 방향으로만 잘 작동한다.

차폐가 얼마나 극단적이냐면, 깊은 협곡의 성장확률은 클러스터 크기의 멱함수보다 빠르게 줄어든다. 그래서 조화측도의 다중프랙탈 스펙트럼은 $q < 0$ 쪽에서 모멘트가 발산해 $f(\alpha)$ 곡선이 끊긴다 — 다중프랙탈 해석이 정직하게 무너지는 희귀한 사례다. 반대로 정보차원은 2차원 단순연결 영역에서 마카로프 정리에 의해 정확히 1이다.

4. 프랙탈 차원[편집]

클러스터의 질량 $N$ 과 회전반경 $R_g$ 사이에

N \sim R_g^{D}

가 성립하고, $\ln N$ 대 $\ln R_g$ 를 회귀하면 $D$ 가 나온다. 상자 세기로도 재지만, DLA는 가지가 성기고 유한 크기 보정이 커서 회전반경 쪽이 훨씬 안정적이다. 값은 2차원에서 $D \approx 1.71$ , 3차원에서 $D \approx 2.5$ . 각각 매립 차원보다 확실히 작아서, 클러스터의 평균 밀도가 $R^{D-d} \to 0$ 으로 자라면서 계속 묽어진다는 뜻이다.

실무적으로는 클러스터 하나만 가지고 $D$ 를 재면 안 된다. DLA는 실현마다 팔 개수와 굵기가 제법 달라지므로 씨앗을 바꿔 수십~수백 개를 굴린 뒤 $\ln R_g$ 를 평균 내고, 초기 수백 입자 구간(아직 자기유사가 아님)과 최종 몇 %(성장이 끝나 팁이 굳은 구간)를 회귀에서 빼는 것이 관례다. 이 처리를 안 하면 $D$ 가 쉽게 $\pm 0.05$ 씩 흔들린다.

여기에 유명한 함정이 둘 있다.

격자 이방성. 정사각 격자 위에서 굴리면 축 방향이 미세하게 유리한데, 이 편향이 성장의 양의 되먹임을 타고 증폭된다. 입자 $10^6$ 개를 넘기면 클러스터가 눈에 띄게 십자 모양으로 변하고 $D$ 추정값도 흔들린다. 그래서 차원을 재려면 격자 없는 오프-래티스 DLA를 쓰는 것이 기본이다.
점착 확률. 닿자마자 붙이는 대신 확률 $p < 1$ 로만 붙이면 입자가 표면을 여러 번 두드리며 안쪽까지 파고들어 클러스터가 조밀해진다. 다만 이건 길이 $\xi \sim 1/p$ 이하에서만 그렇고, 그보다 크게 키우면 같은 $D \approx 1.71$ 로 되돌아온다. 즉 점착 확률은 보편성 부류를 바꾸지 않고 크로스오버 길이만 옮긴다.

격자 301²에서 위튼-샌더 DLA를 실제로 성장시킨다. 발사 원에 놓인 입자가 무작위 걷기를 하다 클러스터에 인접하면 점착 확률 s 로 붙는다. 색은 도착 순서(보라 → 노랑)라, 늦게 온 입자가 팁에만 쌓이고 협곡은 굶는 차폐가 그대로 남는다. 아래 패널은 ln N 대 ln R_g 산점도와 최소제곱 적합선이고, 표시되는 D 는 하드코딩된 값이 아니라 바로 그 회귀의 기울기다(실측 D = 1.688 ± 0.045, R² > 0.998). 정사각 격자 위에서 굴리므로 본문이 말한 축 방향 이방성이 들어 있고, N ≤ 4000 규모에서는 시드 간 산포가 ±0.05 수준이라 소수 둘째 자리까지 읽을 값은 못 된다.

5. 이웃한 성장 모형들[편집]

DLA의 위치를 알려면 규칙 하나씩만 바꾼 사촌들과 나란히 놓는 게 빠르다.

모형	규칙 차이	2D 차원	성격
에덴 성장	표면 아무 곳이나 균등하게	2 (조밀)	차폐 없음, 표면만 거칠다
탄도 응집	직선으로 날아와 붙음	2 (조밀)	확산 대신 관성, 그림자 효과만
DLA	무작위걷기 후 첫 접촉	$\approx 1.71$	조화측도 = 성장확률
클러스터-클러스터 응집	덩어리끼리도 걷다가 합침	$\approx 1.4\sim1.5$	연무·콜로이드 침전
DBM ( $\eta > 1$ )	성장확률 $\propto \lvert\nabla\phi\rvert^{\eta}$	1.71 미만	커질수록 바늘에 가까움

입자를 무작위로 걷게 하느냐가 조밀과 프랙탈을 가르는 유일한 차이라는 게 요점이다. 직선으로 날아오면 안쪽까지 잘 들어가 조밀해지고, 무작위로 걸으면 팁이 먼저 잡아먹어 성기게 갈라진다. 클러스터-클러스터 응집이 DLA보다 더 성긴 것도 같은 논리 — 붙는 쪽까지 이미 성긴 덩어리라 서로의 팁끼리 먼저 만난다.

6. 구현과 응용[편집]

정직하게 무작위걷기를 시키면 입자 하나가 클러스터에 닿기까지 걸음 수가 어마어마해서 $N = 10^5$ 도 버겁다. 표준 가속 기법은 세 가지다. (1) 클러스터 최대 반경 바깥의 발사원에서 출발시키고, (2) 훨씬 먼 소멸반경을 넘어가면 버리고 새로 뽑고, (3) 빈 공간에서는 한 걸음씩 걷는 대신 닿지 않는 최대 원의 반지름만큼 한 번에 점프시킨다. 세 번째가 결정적이라, 이걸 넣으면 $10^8$ 입자급 클러스터도 현실적인 시간에 나온다.³ 붙은 입자를 격자 해시나 사분트리에 넣어 최근접 이웃 질의를 $O(\log N)$ 으로 줄이는 것, 그리고 클러스터 반경이 커질 때 발사원·소멸반경을 함께 키우는 것까지가 한 세트다. 소멸반경을 인색하게 잡으면 멀리 도는 입자를 자꾸 버리게 되어 성장확률이 균등 쪽으로 편향되고, 클러스터가 실제보다 조밀해진다.

응용은 “라플라스장이 성장을 공급하는” 모든 곳이다. 전기도금 결함과 리튬 전지의 덴드라이트 성장, 절연 파괴 경로 예측, 다공성 매질 유동에서 물이 기름을 밀어낼 때의 손가락화, 광물 맥과 망간 침착 무늬, 세균 콜로니 형태. 다만 DLA는 잡음이 지배하는 극한이고 실제 결정 성장에는 표면장력과 결정 이방성이 있어서, 눈송이 같은 규칙적 수지상 구조를 재현하려면 잡음을 인위적으로 줄인 변형(noise-reduced DLA)이나 위상장 모형으로 넘어가야 한다.⁴

7. 관련 문서[편집]

학부생에게 “난수 몇 줄로 만들 수 있는 가장 예쁜 그림”을 물으면 만델브로 집합 다음으로 나오는 것이 대개 DLA다. 실제로 코드 30줄이면 돌아가는데, 그 30줄에서 왜 하필 1.71이 나오는지는 40년 넘게 아무도 해석적으로 증명하지 못했다. ↩
표면장력이 0인 순수 라플라스 성장은 유한 시간에 경계에 첨점(cusp)이 생겨 해가 깨지는 것으로 알려져 있다. DLA가 안 깨지는 이유는 입자 하나 크기라는 최소 길이가 그 특이점을 미리 잘라 주기 때문. 잡음이 정칙화 역할을 하는 흔치 않은 사례다. ↩
초창기 논문들의 클러스터가 수천 입자 규모였던 것도 이 때문이다. “DLA 차원은 1.71”이라는 값이 안정적으로 굳는 데 20년 가까이 걸렸고, 그 사이 1.6대와 1.75 사이를 오가는 보고가 계속 나왔다. 유한 크기 보정을 얕본 대가. ↩
그래서 “눈송이는 DLA다”라는 설명은 절반만 맞다. 눈송이는 육방정 이방성과 표면장력이 무늬를 결정하는 반면, DLA는 그 둘이 전부 없는 극단이다. 굳이 말하면 DLA는 눈송이가 안 되는 이유를 보여 주는 모형에 가깝다. ↩