KPZ 방정식

1. 개요[편집]

KPZ 방정식 Kardar–Parisi–Zhang equation
제안	카르다르 · 파리시 · 장 (1986)
형태	∂th = ν∇²h + (λ/2)(∇h)² + η
1+1차원 지수	α = 1/2, β = 1/3, z = 3/2 (정확)
스케일링 관계	β = α/z, α + z = 2 (갈릴레이 불변성)
λ → 0 극한	에드워즈-윌킨슨 보편성류 (1+1차원 β = 1/4, z = 2)
수학적 지위	고전적으로 ill-posed — 재규격화 필요 (헤어러, 2014)

계면은 왜 거칠어지는가. 그리고 왜 모두 똑같은 방식으로 거칠어지는가.

KPZ 방정식(Kardar–Parisi–Zhang equation)은 무작위로 성장하는 계면의 높이 $h(\mathbf{x},t)$ 를 기술하는

$$\frac{\partial h}{\partial t} = \nu\nabla^{2}h + \frac{\lambda}{2}(\nabla h)^{2} + \eta(\mathbf{x},t)$$

꼴의 확률 편미분방정식이다. 우변은 순서대로 매끄럽게 펴는 표면 완화, 성장이 계면의 법선 방향으로 일어난다는 기하가 만드는 비선형항, 그리고 성장 잡음이며, $\eta$ 는 보통 $\langle\eta(\mathbf{x},t)\eta(\mathbf{x}',t')\rangle = 2D\,\delta^{d}(\mathbf{x}-\mathbf{x}')\delta(t-t')$ 인 백색잡음으로 둔다. 1986년 메흐란 카르다르, 조르조 파리시, 장이청이 단 4쪽짜리 논문에서 제시했다.¹

이 방정식이 특별한 이유는 방정식 자체가 아니라 그 뒤에 딸려온 보편성류다. 에덴 모형의 곰팡이 껍질, 탄도 침착으로 쌓은 박막, 종이가 젖어 올라오는 자국, 대장균 콜로니 가장자리, 비대칭 단순 배제 과정의 밀도장, 무작위 매질 속 방향성 중합체의 자유에너지 — 미시 규칙이 전혀 다른 이 계들이 전부 같은 지수와 같은 요동 분포를 낸다. 이 문서는 그 연속체 방정식과 지수·보편성의 구조를 다룬다. 격자 모형 쪽 이야기와 실제 계면 폭 측정은 에덴 모형 문서에 있고, 보편성류·고정점 개념 자체는 재규격화군이 담당한다.

2. 비선형항은 어디서 오는가[편집]

$(\nabla h)^{2}$ 는 손으로 끼워 넣은 항이 아니다. 성장은 계면에 수직인 방향으로 일어난다는 기하학적 사실 하나에서 곧바로 나온다.

계면이 법선 방향으로 일정 속도 $v$ 로 전진한다고 하자. 높이 함수 $h$ 로 보면 법선 방향 변위 $v\,\delta t$ 의 수직 성분은 $v\,\delta t\sqrt{1+(\nabla h)^{2}}$ 이므로

$$\frac{\partial h}{\partial t} = v\sqrt{1+(\nabla h)^{2}} \;\approx\; v + \frac{v}{2}(\nabla h)^{2} + O\!\big((\nabla h)^{4}\big)$$

가 된다. 상수 $v$ 는 좌표계를 함께 움직이면 사라지고, 남는 것이 정확히 $\lambda = v$ 인 KPZ 항이다. 즉 λ는 계면의 평균 성장 속도 그 자체이고, 성장이 있는 한 이 항은 꺼지지 않는다. 대칭성 논리로도 같은 결론이 나온다 — $h \to h + c$ 병진 불변성 때문에 $h$ 자체는 못 나오고, 계면이 위아래로 대칭이 아니므로 ($h \to -h$ 가 금지) 최저차 허용 항이 $(\nabla h)^{2}$ 다.

부호는 물리적으로 의미가 있다. $\lambda > 0$ 이면 튀어나온 부분이 옆으로도 자라 돌출부가 더 빨리 퍼지고, $\lambda < 0$ 이면 식각(etching)처럼 반대가 된다. 지수는 $|\lambda|$ 의 부호에 무관하지만, 뒤에 나올 트레이시-위덤 분포의 왜도 방향은 부호를 따라 뒤집힌다.

$u = \nabla h$ 로 바꾸면 방정식이 잡음 있는 버거스 방정식

$$\partial_t \mathbf{u} + \lambda(\mathbf{u}\cdot\nabla)\mathbf{u} = \nu\nabla^{2}\mathbf{u} + \nabla\eta$$

이 된다. KPZ의 $\lambda$ 항이 버거스의 이류항과 같은 물건이라는 뜻이고, 충격파 형성·점성해의 언어가 그대로 옮겨온다. 한편 콜-호프 변환 $W = \exp(\lambda h/2\nu)$ 를 쓰면 곱셈 잡음이 붙은 확률 열방정식 $\partial_t W = \nu\nabla^{2}W + (\lambda/2\nu)W\eta$ 로 선형화된다. 이 세 방정식이 같은 대상의 세 얼굴이라는 것이 KPZ 이론의 출발점이다.

3. 거칠기 지수와 패밀리-비식 스케일링[편집]

계면의 거칠기는 폭

$$W(L,t) = \Big\langle \big[h(\mathbf{x},t) - \bar h(t)\big]^{2}\Big\rangle^{1/2}$$

로 잰다. 유한한 계 크기 $L$ 에서 $W$ 는 처음에 멱으로 자라다가 상관길이가 $L$ 을 넘는 순간 포화하며, 두 국면이 하나의 스케일링 함수로 묶인다.

$$W(L,t) \sim L^{\alpha} f\!\left(\frac{t}{L^{z}}\right), \qquad f(u) \sim \begin{cases} u^{\beta}, & u \ll 1 \\ \text{const}, & u \gg 1\end{cases}$$

이것이 패밀리-비식 스케일링(Family–Vicsek, 1985)이고, 세 지수는 $\beta = \alpha/z$ 로 묶여 둘만 독립이다. $\alpha$ 는 거칠기 지수(포화 폭이 계 크기에 어떻게 붙는가), $\beta$ 는 성장 지수(초기에 얼마나 빨리 거칠어지는가), $z$ 는 동역학 지수(상관이 얼마나 빨리 퍼지는가)다.

1+1차원에서는 세 값이 정확히 알려져 있다.

$$\alpha = \tfrac12, \qquad \beta = \tfrac13, \qquad z = \tfrac32$$

두 개의 독립된 논증이 이 값을 못 박는다. 첫째, 1차원 KPZ의 정상 상태 분포는 $\lambda$ 와 무관하게 브라운 운동(즉 선형 이론과 같은 것)이라는 요동-소산 관계가 성립해 $\alpha = 1/2$ 가 정확히 고정된다. 둘째, 갈릴레이 불변성이 남은 하나를 준다.

4. 갈릴레이 불변성과 α + z = 2[편집]

방정식은 기울임 변환

$$h(\mathbf{x},t) \;\to\; h(\mathbf{x}+\lambda\boldsymbol{\epsilon}t,\,t) + \boldsymbol{\epsilon}\cdot\mathbf{x}, \qquad |\boldsymbol{\epsilon}| \ll 1$$

아래에서 형태가 변하지 않는다. 계면을 살짝 기울이면 그 기울기가 만드는 여분의 성장 속도가 정확히 $\lambda\boldsymbol{\epsilon}$ 만큼의 이동으로 흡수된다는 뜻이고, 버거스 방정식 쪽에서 보면 이름 그대로 갈릴레이 변환이다. 이 대칭성이 $\lambda$ 를 재규격화 흐름에서 불변으로 고정하므로, 조대화 아래에서 $\lambda$ 의 스케일링 차원 $\alpha + z - 2$ 는 0이어야 한다.

$$\boxed{\;\alpha + z = 2\;}$$

이 관계는 차원과 무관하게 성립한다. 1+1차원에서 $\alpha = 1/2$ 와 결합하면 $z = 3/2$, 따라서 $\beta = 1/3$ 이다.² 2+1차원에서는 정확해가 없어 수치에 의존하며, 대규모 시뮬레이션이 주는 값은 $\alpha \approx 0.39$, $\beta \approx 0.24$, $z \approx 1.61$ 근처다 — $\alpha + z = 2$ 는 여전히 잘 맞는다.

5. EW와의 대비 — λ를 끄면 무슨 일이 일어나는가[편집]

$\lambda = 0$ 으로 두면 방정식은 선형 확률 확산방정식, 즉 에드워즈-윌킨슨(EW) 방정식이 된다. 푸리에 변환 한 번으로 정확히 풀리고 지수는

$$\alpha_{\rm EW} = \frac{2-d}{2}, \qquad \beta_{\rm EW} = \frac{2-d}{4}, \qquad z_{\rm EW} = 2$$

다. 1+1차원($d=1$)에서 $\alpha_{\rm EW} = 1/2$ 로 KPZ와 값이 같다. 실무적으로 아주 고약한 함정인데, 포화 폭만 재서 $\alpha \approx 0.5$ 를 얻고 “KPZ 확인”이라고 쓰면 아무것도 확인하지 못한 것이 된다. 두 류를 가르는 것은 $\beta$($1/3$ 대 $1/4$)와 $z$($3/2$ 대 $2$)이지 $\alpha$ 가 아니다.

재규격화군으로 보면 사정이 이렇다. 결합상수 $g = \lambda^{2}D/\nu^{3}$ 의 흐름을 1루프로 계산하면 $d < 2$ 에서 비선형항이 관련(relevant) 연산자라 아무리 작은 $\lambda$ 도 결국 강결합 KPZ 고정점으로 흘러간다. $d = 2$ 는 주변적이고, $d > 2$ 에서는 약결합 EW 고정점이 국소적으로 안정해져 거칠어짐 전이(roughening transition)가 생긴다 — 잡음이 약하면 계면이 EW처럼 매끈하게 남고, 임계값을 넘겨야 KPZ 강결합 상으로 넘어간다. 상부 임계 차원이 4인지 아예 없는지는 30년 넘게 논쟁 중이다.

6. 트레이시-위덤 — 지수 너머의 보편성[편집]

2000년대에 밝혀진 것은 지수뿐 아니라 높이 요동의 확률분포 전체가 보편적이라는 사실이다. 1+1차원에서

$$h(x,t) \simeq v_{\infty}t + (\Gamma t)^{1/3}\,\chi$$

로 쓰면 재척도된 변수 $\chi$ 의 분포가 계에 무관하게 정해진다. 다만 초기조건의 기하가 어느 분포인지를 가른다.

초기조건	계면 기하	$\chi$ 의 분포
좁은 쐐기(점 씨앗)	곡면, 방사상 성장	트레이시-위덤 GUE
평평한 계면	평면	트레이시-위덤 GOE
브라운 초기조건	정상 상태	바이크-레인스 $F_{0}$

트레이시-위덤 분포는 원래 무작위 행렬 앙상블의 최대 고유값 분포로 발견된 것이다. 자라는 곰팡이의 가장자리 요동이 거대 무작위 행렬의 최대 고유값과 같은 분포를 따른다는 이 대응은, 처음 보고됐을 때 그럴듯한 우연으로 취급됐다가 지금은 KPZ 이론의 중심이 되었다. 실험적 확인은 2010년 다케우치와 사노가 네마틱 액정의 전기대류 난류 영역 경계에서 두 기하(원형·평면)를 모두 만들어 GUE와 GOE를 갈라 본 것이 결정적이었다. 수학 쪽에서는 2010년 아미르-코윈-쿼스텔, 사사모토-슈폰, 칼라브레제-르두살-로소가 좁은 쐐기 초기조건의 한 점 분포를 정확히 계산해 냈다.

7. Ill-posed 방정식과 정칙 구조[편집]

여기서 불편한 사실을 짚어야 한다. KPZ 방정식은 쓰인 그대로는 수학적으로 말이 되지 않는다. 백색잡음이 구동하는 1차원 해는 시간마다 브라운 운동만큼 거칠어서 $h(\cdot,t)$ 가 공간에서 미분 불가능하고, 따라서 $(\nabla h)^{2}$ 는 분포의 제곱이라 정의되지 않는다. 매끄러운 잡음으로 정칙화한 뒤 극한을 취하면 발산하는 상수가 나오고, 무한대를 빼 주는 재규격화를 해야 유한한 극한이 남는다.

• 콜-호프 해(베르티니-지아코민, 1997). $W = e^{\lambda h/2\nu}$ 가 만족하는 확률 열방정식은 곱셈 잡음이 선형이라 이토 적분으로 잘 정의된다. 이 $W$ 로부터 $h = (2\nu/\lambda)\log W$ 를 해로 정의하는 우회로. 잘 작동하지만 비선형항을 조금만 바꿔도(예: $(\nabla h)^{2}$ 대신 다른 함수) 변환이 통하지 않는다.
• 정칙 구조 이론(마르틴 헤어러, 2013–14). 거친 신호를 국소적으로 “테일러 전개하되, 기저를 다항식이 아니라 잡음이 만든 대상들로 잡는” 대수적 틀을 세워, 발산하는 항을 체계적으로 상쇄하는 재규격화 절차를 정의한다. KPZ를 직접, 그리고 우회 없이 푸는 데 성공했고 헤어러는 2014년 필즈상을 받았다.³ 비슷한 시기에 구비넬리 계열의 패러컨트롤드 미적분이 독립적으로 같은 문제를 다뤘다.

이 이야기는 순수하게 추상적이지 않다. “무한대를 빼야 한다”는 사실은 이산화가 유효 파라미터를 바꾼다는 수치적 현상의 다른 얼굴이기 때문이다.

8. 수치적으로 다루기 — 함정 모음[편집]

KPZ는 “그냥 전진 오일러로 돌리면 되겠지” 하고 접근했다가 여러 번 물리는 방정식이다.

• 이산화가 $\lambda$ 를 바꾼다. 잡음의 진폭 $D$, 격자 간격 $\Delta x$, 시간 스텝 $\Delta t$ 를 정하는 순간 유효 $\nu,\lambda,D$ 가 연속 방정식의 값과 달라진다(위의 재규격화 상수가 격자에서는 유한한 이동으로 나타난다). 그래서 서로 다른 이산화가 서로 다른 유효 결합 $g$ 로 출발하고, 점근 영역에 도달하는 시간이 크게 달라진다.
• 비선형항 이산화가 정상 분포를 결정한다. 순진한 중심차분 $[(h_{i+1}-h_{i-1})/2\Delta x]^{2}$ 는 1차원 KPZ가 가져야 할 브라운 정상 분포를 정확히 재현하지 못한다. 람-신(1998)이 제시한 $$\frac{1}{3\Delta x^{2}}\Big[(h_{i+1}-h_{i})^{2} + (h_{i+1}-h_{i})(h_{i}-h_{i-1}) + (h_{i}-h_{i-1})^{2}\Big]$$ 형태가 정상 상태를 보존해 사실상의 표준이 되었다.
• 강한 비선형에서 발산한다. $\lambda$ 를 키우면 명시적 스킴이 유한 시간에 폭발하는 것이 반복 보고됐다(다스굽타·다스 사르마·김, 1996–97). 원인은 국소적으로 큰 기울기가 생기면 $(\nabla h)^{2}$ 가 그 기울기를 더 키우는 양의 되먹임이고, 확산항의 $\Delta t \le \Delta x^{2}/2\nu$ 조건을 지켜도 막히지 않는다. 대책은 $\Delta t$ 를 훨씬 줄이거나, 잡음 감축(noise reduction) 기법을 쓰거나, 애초에 방정식 대신 같은 보편성류의 격자 모형을 돌리는 것이다. 실무에서 KPZ 지수를 재고 싶을 때 에덴 모형이나 단일 스텝 모형을 쓰는 이유가 여기 있다.
• 점근 영역은 늦게 온다. 유효 지수 $\beta_{\rm eff}(t) = d\log W/d\log t$ 는 보정항 때문에 오랫동안 $1/3$ 아래에서 기어간다. 상수 성분을 빼지 않고 기울기를 뽑으면 $\beta \approx 0.25$ 가 나와 EW로 오독하기 딱 좋다. 여러 $L$ 에서 재고 $t/L^{z}$ 축으로 곡선 붕괴(collapse)를 확인하는 것이 최소한의 검증이다.
• 잡음 스케일링. 백색잡음의 이산 대응은 $\eta_i^{n} = \sqrt{2D/(\Delta x^{d}\Delta t)}\,\xi$ ($\xi$ 는 표준정규)다. $\sqrt{\Delta t}$ 를 빠뜨리는 것이 이 바닥의 전형적인 첫 버그.

9. 관련 문서[편집]

• 에덴 모형 · 확산 제한 응집 · 프랙탈
• 버거스 방정식 · 쿠라모토-시바신스키 방정식 · 점성해
• 재규격화군 · 상전이 · 멱법칙 · 허스트 지수
• 확률미분방정식 · 이토 적분 · 브라운 운동
• 극값 통계 · 무작위 행렬 · 방향성 중합체 · 비대칭 단순 배제 과정
• 몬테카를로 방법 · 유한차분법

10. Footnotes[편집]

1. Kardar, Parisi, Zhang (1986), Phys. Rev. Lett. 56, 889 — “Dynamic scaling of growing interfaces”. 본문 4쪽에 그림 한 장 없고, 지금 이 문서보다 짧다. 인용 수는 만 회를 넘겼다. 파리시는 나중에 전혀 다른 업적(무질서계의 복잡성)으로 2021년 노벨 물리학상을 받는다. ↩

2. “그래서 $\beta$ 가 왜 하필 $1/3$ 이냐”에 대한 짧은 답은 없다. $\alpha=1/2$ 는 정상 분포에서, $\alpha+z=2$ 는 대칭성에서 나오고, 나머지는 나눗셈이다. 이 두 개의 독립된 사실이 우연처럼 맞물려 정확한 유리수 지수를 뱉는다는 점이 1차원 KPZ의 매력이자, 2차원에서 아무도 같은 재주를 못 부리는 이유다. ↩

3. 필즈상 인용문에 “확률 편미분방정식”이 등장한 것이 이때가 처음이다. 물리학자들이 1986년에 써 놓고 30년간 잘 쓰던 방정식이 실은 정의되지도 않은 식이었다는 사실은, 이 바닥에서 “일단 돌려”가 얼마나 자주 옳은지를 보여 주는 사례로도, 반대로 얼마나 위태로운지를 보여 주는 사례로도 인용된다. ↩