브라운 운동

편집 역사 토론
계산물리 통계 물리 마지막 수정: 2026-08-02 05:04:18

1. 개요[편집]

브라운 운동
Brownian Motion
최초 관찰로버트 브라운 (1827)
이론아인슈타인 (1905), 스몰루호프스키 (1906), 랑주뱅 (1908)
검증장 페랭 (1908~1913), 노벨물리학상 1926
핵심 관계식<r²> = 2dDt
수학적 모형위너 과정 (허스트 지수 H = 1/2)
궤적의 차원하우스도르프 차원 2 (d ≥ 2)

원자를 본 사람은 아무도 없었다. 그런데 꽃가루가 떨리는 걸 보고 원자의 개수를 셌다.

브라운 운동(Brownian motion)은 유체에 떠 있는 미세 입자가 주변 분자들의 무질서한 충돌 때문에 끊임없이 불규칙하게 움직이는 현상이며, 동시에 그 현상을 이상화한 확률과정(위너 과정)을 가리키는 이름이기도 하다. 물리 현상으로서는 열요동의 가장 눈에 잘 보이는 얼굴이고, 수학 대상으로서는 연속이지만 어디서도 미분되지 않는 괴물이며, 시뮬레이션에서는 확률미분방정식 적분기의 표준 시험대다.

범위 정리를 먼저 해 두자. 관성이 살아 있는 유한 마찰 방정식과 그 적분기(BAOAB)는 랑주뱅 동역학이, 기체 분자 자체의 충돌 간격은 평균자유행로가, 자기유사성 지수의 일반론은 허스트 지수가 맡는다. 이 문서는 왜 변위가 t\sqrt{t} 로 자라는가, 그 확률과정은 정확히 무엇인가, 그걸 어떻게 수치적으로 굴리는가에 집중한다.

2. 꽃가루에서 아보가드로 수까지[편집]

1827년 식물학자 로버트 브라운은 물에 띄운 꽃가루에서 튀어나온 미세 입자들이 쉬지 않고 떨리는 것을 현미경으로 관찰했다. 브라운이 위대한 지점은 발견 자체가 아니라 대조 실험이다. 그는 “살아 있으니 움직인다”는 설명을 배제하려고 죽은 표본, 화석 조각, 심지어 스핑크스에서 긁어낸 돌가루까지 같은 방식으로 관찰해 똑같은 떨림을 확인했다. 생명 현상이 아니라는 것까지는 확정했지만, 원인은 이후 80년 가까이 미궁이었다.1

1905년 아인슈타인은 이 떨림을 원자의 실재성을 증명할 도구로 재해석했다. 그의 논문은 “지금 다루는 운동이 이른바 브라운 운동과 같은 것일 수도 있다”는 조심스러운 문장으로 시작하는데, 실제 목표는 처음부터 하나였다 — 분자 가설이 옳다면 관측 가능한 양이 얼마여야 하는지 숫자로 못 박아, 실험으로 판정 가능하게 만드는 것. 여기서 나온 것이 평균제곱변위 공식과 확산계수 표현이다. 1906년 스몰루호프스키가 독립적으로 같은 결론에 도달했고, 1908년 랑주뱅은 확률적 힘을 뉴턴 방정식에 직접 집어넣는 “무한히 더 단순한” 유도를 제시했다.

판정은 장 페랭의 몫이었다. 1908년부터 그는 입자 크기가 균일한 수지 에멀션을 만들어 궤적을 30초 간격으로 찍고, 침강 평형과 평균제곱변위 두 경로에서 각각 아보가드로 수를 뽑아 약 6.7×10236.7 \times 10^{23} 이라는 값을 얻었다. 서로 완전히 다른 원리의 측정이 같은 숫자로 수렴하자 원자론에 남아 있던 저항은 사실상 끝났고, 페랭은 1926년 노벨물리학상을 받았다.2

3. 왜 tt 가 아니라 t\sqrt{t} 인가[편집]

dd 차원에서 브라운 입자의 평균제곱변위는 시간에 선형으로 자란다.

⟨∣r(t)−r(0)∣2⟩=2dDt\langle |\mathbf{r}(t) - \mathbf{r}(0)|^2 \rangle = 2dDt

1차원이면 2Dt2Dt, 2차원이면 4Dt4Dt, 3차원이면 6Dt6Dt 다. 따라서 전형적 변위 자체는 2dDt∝t\sqrt{2dDt} \propto \sqrt{t} 이고, 이 “제곱근” 하나가 확산과 이동(drift)을 가르는 경계선이다.

이유는 통계다. 시간 tt 동안 방향이 독립적으로 갱신되는 걸음을 N∝tN \propto t 번 밟았다고 하면, 총 변위는 독립 벡터들의 합이다. 교차항 ⟨Δri⋅Δrj⟩\langle \Delta\mathbf{r}_i \cdot \Delta\mathbf{r}_j \rangle 가 i≠ji \ne j 에서 0이 되므로 분산은 더해지지만 진폭은 더해지지 않는다. 즉 ⟨r2⟩∝N∝t\langle r^2 \rangle \propto N \propto t, ∣r∣∝N|r| \propto \sqrt{N}. 일정한 속도로 흘러가는 이류(advection)가 tt 에 비례하는 것과 대비된다 — 이 두 항이 한 방정식 안에서 경쟁하는 것이 대류-확산 방정식이고, 그 비율이 페클레 수다. 확산은 짧은 시간엔 이기고 긴 시간엔 반드시 진다.3

확산계수를 물성으로 이어 주는 것이 아인슈타인 관계다. 이동도 μ\mu(단위 힘당 종단속도)를 쓰면

D=μkBT,μ=16πηa  ⇒  D=kBT6πηaD = \mu k_B T, \qquad \mu = \frac{1}{6\pi\eta a} \;\Rightarrow\; D = \frac{k_B T}{6\pi \eta a}

가 되고, 후자를 스토크스-아인슈타인 식이라 부른다. 좌변은 요동(무작위 걸음의 크기), 우변의 μ\mu 는 소산(끌림 저항)이며 둘을 kBTk_BT 가 묶는다. 랑주뱅 동역학 문서가 다루는 요동-소산 정리의 가장 오래된 구체적 사례가 바로 이 식이다. 숫자를 넣어 보면 물속(η=10−3 Pa⋅s\eta = 10^{-3}\,\mathrm{Pa\cdot s})의 지름 1 μm 콜로이드가 300 K에서 D≈0.44 μm2/sD \approx 0.44\ \mu\mathrm{m}^2/\mathrm{s}, 1초 동안 2차원 변위가 4Dt≈1.3 μm\sqrt{4Dt} \approx 1.3\ \mu\mathrm{m} 정도다. 현미경으로 딱 보기 좋은 스케일이라 페랭이 손으로 셀 수 있었던 것.

4. 위너 과정 — 연속인데 미분이 안 된다[편집]

물리에서 이상화를 끝까지 밀면 수학의 위너 과정 W(t)W(t) 가 나온다. 정의는 네 줄이다.

  • W(0)=0W(0) = 0.
  • 독립 증분: 겹치지 않는 구간의 증분끼리 독립.
  • 정상 가우시안 증분: W(t)−W(s)∼N(0, t−s)W(t) - W(s) \sim \mathcal{N}(0,\, t-s).
  • 경로가 거의 확실히 연속.

증분이 가우시안인 것은 가정이라기보다 결론에 가깝다. 수많은 분자 충돌의 합이므로 중심극한정리가 개별 충돌의 분포 모양을 지워 버린다. 이 때문에 브라운 운동은 미시 세부에 놀랄 만큼 둔감하고, 거꾸로 말하면 궤적만 봐서는 유체의 분자 구조를 알아낼 수 없다.

가장 유명한 성질은 거의 확실히 어디서도 미분 불가능이라는 것이다(페일리-위너-지그문트, 1933). 증분의 크기가 Δt\sqrt{\Delta t} 스케일이므로 차분몫은 Δt/Δt\Delta t/\sqrt{\Delta t} 의 역수, 즉 1/Δt1/\sqrt{\Delta t} 로 발산한다. 물리적으로도 정직한 결론이다 — “브라운 입자의 순간 속도”를 시간 해상도를 높여 가며 재면 값이 계속 커지고, 진짜 열속도 kBT/m\sqrt{k_BT/m} 는 관성 완화 시간(마이크로초 이하) 밑으로 내려가야 비로소 보인다.4

또 하나는 자기유사성이다.

W(at)=da W(t)W(at) \overset{d}{=} \sqrt{a}\,W(t)

시간을 100배 늘이고 공간을 10배 줄이면 통계적으로 같은 그림이 나온다. 지수 1/21/2 가 허스트 지수 HH 의 특수값이고, 이를 H∈(0,1)H \in (0,1) 로 일반화해 증분에 장기 상관을 넣은 것이 분수 브라운 운동이다(H>1/2H > 1/2 이면 지속성, H<1/2H < 1/2 이면 반지속성). 궤적의 기하도 자기유사적이라, 2차원 이상에서 브라운 경로의 하우스도르프 차원은 정확히 2다 — 선이지만 면적을 채울 기세로 헤매고 다닌다는 뜻이며, 이 감각은 프랙탈 문서와 통한다. 확산하는 입자가 서로 달라붙으면 확산 제한 응집의 성긴 가지 구조가 나오는데, 그 프랙탈 차원이 1.71 언저리로 낮은 이유도 결국 브라운 궤적이 표면부터 훑고 지나가기 때문이다.

5. 수치 시뮬레이션 — Δt\sqrt{\Delta t} 함정[편집]

가장 단순한 브라운 동역학(과감쇠 극한)은 위치만 갱신한다.

rn+1=rn+DkBTF(rn) Δt+2D Δt  ξn,ξn∼N(0,I)\mathbf{r}_{n+1} = \mathbf{r}_n + \frac{D}{k_B T}\mathbf{F}(\mathbf{r}_n)\,\Delta t + \sqrt{2D\,\Delta t}\;\boldsymbol{\xi}_n, \qquad \boldsymbol{\xi}_n \sim \mathcal{N}(0, I)

이것이 오일러-마루야마 도식이다. 결정론적 오일러법과 겉모습이 같지만 결정적 차이가 하나 있다. 잡음항이 Δt\Delta t 가 아니라 Δt\sqrt{\Delta t} 로 스케일된다. 확률미분방정식 입문자가 가장 자주 밟는 지뢰가 여기로, 잡음에도 Δt\Delta t 를 곱해 버리면 스텝을 줄일수록 확산이 사라져 계가 조용히 0 K로 얼어붙는다. 반대로 Δt\sqrt{\Delta t} 를 쓰면 스텝을 줄여도 ⟨r2⟩=2dDt\langle r^2\rangle = 2dDt 가 유지된다. 검증도 쉽다 — 힘을 끄고 돌린 뒤 로그-로그 평면에서 기울기 1이 나오는지 보면 된다.

가우시안 증분은 균등난수에서 만들어야 하고, 그 표준 도구가 박스-뮐러 변환이다(성능 때문에 지구랫 방법으로 갈아타기도 한다). 균등난수를 그냥 쓰면 한 스텝 통계는 틀리지만 여러 스텝을 합치면 다시 중심극한정리가 구해 주는데, 짧은 시간 통계를 볼 거면 얌전히 가우시안을 쓰는 게 국룰이다.

2차원 무작위걷기 입자 800개를 동시에 굴리며 앙상블 평균제곱변위를 매 스텝 누적한다. 아래 로그-로그 패널의 기울기 α와 확산계수 D는 회귀로 실측한 값이고, 축당 스텝 분산이 σ²이므로 이론값 D = σ²/2Δt 옆에 상대오차와 함께 표시된다. 2차원에서 MSD 추정량의 상대 표준오차는 1/√N이라 N=800이면 3.5% 수준이고, 화면의 어긋남 대부분은 그 통계 산포다.

정확도 차수도 결정론적 감각과 다르다. 오일러-마루야마는 강수렴 차수 0.5, 약수렴 차수 1이다. 궤적 하나하나를 정확히 따라가려면(강수렴) 형편없지만, 평균값·분산 같은 분포 통계만 필요하면(약수렴) 1차라 쓸 만하다. 분자 시뮬레이션은 거의 항상 후자만 필요하므로 오일러-마루야마가 여전히 현역이다. 궤적 자체가 필요하면 밀스타인 도식 이상으로 올라가야 한다.

6. 이토냐 스트라토노비치냐[편집]

잡음 세기가 위치에 의존하면(D=D(r)D = D(\mathbf{r}), 예컨대 벽 근처에서 유체역학적 저항이 커지는 경우) ∫b(Xt) dWt\int b(X_t)\,dW_t 를 어느 시점의 bb 로 평가할지가 답을 바꾼다. 리만 적분에서는 구간 안 어디를 찍든 같은 값에 수렴하지만, WW 의 2차 변동이 유한하기 때문에 확률적분에서는 그렇지 않다.

  • 이토 해석: 구간 왼쪽 끝에서 평가. 미래를 안 보는 비예측적(non-anticipating) 구성이라 마팅게일 성질이 살아 있고, 금융공학과 필터링 이론의 기본값이다. 대가로 연쇄법칙이 이토 보조정리로 바뀌어 12b2f′′\frac12 b^2 f'' 항이 추가로 붙는다.
  • 스트라토노비치 해석: 구간 중점에서 평가. 보통의 연쇄법칙이 그대로 성립하고, 상관시간이 유한한 물리적 잡음을 백색 극한으로 보낼 때 나오는 자연스러운 극한(웡-자카이 정리)이라 물리 쪽이 선호한다.

둘은 표류항 12b ∂xb\frac12 b\,\partial_x b 만큼 차이 나므로 서로 정확히 변환된다. 즉 틀린 쪽은 없고, 안 밝히는 쪽이 틀렸다 — 잡음이 상수(가법적)면 두 해석이 일치하므로 대부분의 MD 코드는 이 논쟁을 만날 일이 없지만, 위치 의존 확산계수를 넣는 순간 논문 메서드 절에 해석을 명시해야 한다. 상세는 이토 적분과 확률미분방정식 문서로 미룬다.

7. 어디에 쓰이나[편집]

  • 연속체 ↔ 입자 대응. 위너 과정의 확률밀도가 따르는 방정식이 정확히 확산 방정식이라, 같은 물리를 대류-확산 방정식으로 풀 수도 있고 입자를 뿌려 풀 수도 있다. 고차원·복잡 형상에서는 후자(무작위걷기 몬테카를로)가 유리한 경우가 많다.
  • 표본추출 엔진. 과감쇠 랑주뱅 방정식에 힘 −∇U-\nabla U 를 넣으면 정상분포가 볼츠만 분포가 되고, 이걸 이산화한 뒤 메트로폴리스 보정을 얹은 것이 MALA다. 마르코프 연쇄 몬테카를로와 몬테카를로 방법의 확산 기반 계열이 전부 여기서 갈라져 나온다.
  • 스케일 감각. 같은 무작위걸음이라도 기체 분자는 평균자유행로가 수십 nm에 충돌 간격이 나노초 이하라 걸음 하나가 개별 사건으로 식별되지만, 콜로이드는 초당 101910^{19} 회 규모의 충돌을 뭉뚱그린 결과만 본다. 둘 다 확산이되 “한 걸음”의 의미가 완전히 다르다.
  • 실험 도구. 입자 추적으로 DD 를 재고 스토크스-아인슈타인을 뒤집으면 유체의 점도나 입자 크기가 나온다(미세유변학, 동적광산란). 세포질처럼 ⟨r2⟩∝tα\langle r^2\rangle \propto t^\alpha 에서 α<1\alpha < 1 인 이상확산이 보이면, 그때부터는 분수 브라운 운동이나 연속시간 무작위걷기의 영역이다.

8. 관련 문서[편집]

9. Footnotes[편집]

  1. 그래서 영어권에서는 “Brownian motion”이 브라운의 이름을 딴 것치고 브라운이 한 일이 적다는 농담이 있다. 사실 떨림 자체는 현미경을 들여다본 사람이면 다 봤고, 1785년 얀 잉엔하우스가 알코올 위 목탄 가루에서 같은 것을 기록했다. 브라운의 공은 “살아 있어서 그렇다”는 가장 그럴듯한 오답을 실험으로 사살했다는 데 있다. ↩

  2. 아인슈타인의 1905년 논문보다 5년 앞서 루이 바슐리에가 박사논문에서 사실상 같은 확률과정을 옵션 가격에 적용했다. 지도교수 푸앵카레의 평가는 “주제가 다소 특이함”이었고, 논문은 반세기 가까이 묻혔다가 금융공학이 뜨면서 발굴됐다. 브라운 운동은 물리학보다 파리 증권거래소에서 먼저 수식이 됐던 셈이다. ↩

  3. 확산에만 의존하는 수송이 얼마나 답답한지는 방 크기로 실감할 수 있다. 공기 중 향 분자의 DD 가 대략 10−5 m2/s10^{-5}\ \mathrm{m^2/s} 이므로 5 m를 순수 확산으로 건너오려면 t∼L2/2Dt \sim L^2/2D, 즉 열흘이 넘는다. 실제로 몇 초 만에 냄새가 오는 건 전적으로 공기의 대류 덕이다. 확산 하나만 믿고 설계한 반응기가 왜 안 돌아가는지도 같은 계산으로 설명된다. ↩

  4. 열속도 스케일의 “탄도적” 구간 ⟨r2⟩≃(kBT/m)t2\langle r^2 \rangle \simeq (k_BT/m)t^2 은 아인슈타인 본인이 관측 불가능하다고 못 박았던 영역인데, 2010년대에 광집게로 붙잡은 마이크로 구슬을 나노초·옹스트롬 해상도로 추적해 결국 측정됐다. 100년 묵은 “이건 못 잰다”가 장비의 발전으로 깨진 흔치 않은 사례. ↩