크라머스 탈출률

편집 역사 토론
계산물리 통계 분자동역학 마지막 수정: 2026-09-15 04:32:15
상위 문서: 포커-플랑크 방정식

1. 개요[편집]

크라머스 탈출률(Kramers escape rate)은 열잡음에 흔들리는 입자가 퍼텐셜 우물에 갇혀 있다가 장벽을 넘어 탈출하는 단위시간당 확률이다. 1940년 한드릭 크라머스가 포커-플랑크 방정식을 장벽 위에서 정상 흐름으로 풀어 얻었고, 화학반응 속도론을 동역학과 마찰의 문제로 다시 정의한 결과가 됐다.

핵심 결과는 아레니우스 형태에 전인자(prefactor)가 마찰에 의존한다는 점이다.

k=ω02π⋅γ2/4+ωb2−γ/2ωb  e−ΔU/kBTk = \frac{\omega_0}{2\pi} \cdot \frac{\sqrt{\gamma^2/4 + \omega_b^2} - \gamma/2}{\omega_b} \; e^{-\Delta U / k_B T}

ω0\omega_0는 우물 바닥의 각진동수, ωb\omega_b는 장벽 정점의 곡률에서 오는 허진동수, γ\gamma는 마찰계수, ΔU\Delta U는 장벽 높이다. 지수 인자는 아레니우스 방정식과 같고, 앞의 계수가 크라머스의 기여다.1

이중우물 U(x)=ΔU(x²−1)² 에서 과감쇠 랑주뱅 워커 720개를 x=−1 에서 출발시켜 x=+1 에 닿는 첫 통과 시간을 실측한다(Euler–Maruyama, kT=1, dt=0.0015). 화면의 k 측정은 사건 수를 총 워커-시간으로 나눈 갱신률, k 구적은 같은 MFPT 이중적분을 격자 2400으로 정확히 푼 값, k 크라머스는 안장점 점근식이다. 측정은 구적값을 1% 안에서 재현하며, 점근식은 ΔU=2.5 에서 22% 과대평가하다 ΔU=7 에서 7%로 좁혀진다. 과감쇠 서술이므로 ζ 는 시간축을 되파는 역할만 하고(k ∝ 1/ζ), 저마찰 가지와 크라머스 터노버는 여기서 나오지 않는다.

2. 세 가지 마찰 영역[편집]

같은 장벽인데도 마찰 크기에 따라 탈출 속도의 물리가 완전히 달라진다.

영역조건거동
저마찰 (에너지 확산)γ≪ωb\gamma \ll \omega_bk∝γk \propto \gamma
중간 (전이상태 이론 극한)γ∼ωb\gamma \sim \omega_bk≈(ω0/2π)e−βΔUk \approx (\omega_0/2\pi)e^{-\beta\Delta U}
고마찰 (공간 확산)γ≫ωb\gamma \gg \omega_bk∝1/γk \propto 1/\gamma

마찰이 커지면 속도가 느려지는 건 직관적이다 — 점성이 클수록 장벽을 넘기 힘들다. 반직관적인 쪽은 저마찰에서도 느려진다는 것이다. 마찰이 거의 없으면 열 저장소와 에너지를 주고받는 경로 자체가 없어서, 입자가 장벽을 넘을 만큼의 에너지를 얻는 데 시간이 걸린다. 두 영역 사이에 최대값이 존재하고, 이 비단조 형태를 크라머스 터노버(turnover)라고 부른다. 멜니코프-메시코프가 두 극한을 잇는 통일식을 1986년에 완성했다.

고마찰 극한의 식은 실무에서 가장 많이 쓰인다.

k≈ω0ωb2πγ e−ΔU/kBTk \approx \frac{\omega_0 \omega_b}{2\pi\gamma}\, e^{-\Delta U/k_B T}

이게 스몰루호프스키 극한이며, 랑주뱅 동역학의 관성을 버린 과감쇠(overdamped) 서술과 짝을 이룬다. 용액 속 생체분자나 콜로이드는 거의 항상 이쪽이다.

3. 전이상태 이론과의 관계[편집]

전이상태 이론(TST)은 장벽 정점을 지나는 플럭스를 세면서 한 번 넘은 것은 돌아오지 않는다고 가정한다. 크라머스 식의 중간 마찰 극한이 정확히 TST와 일치하고, 나머지 영역에서는 TST가 속도를 과대평가한다. 그 비율을 투과계수 κ=k/kTST≤1\kappa = k/k_{\text{TST}} \le 1 로 정의한다.

κ<1\kappa < 1의 물리적 원인은 재교차(recrossing)다. 마찰이 크면 장벽을 넘자마자 용매에 얻어맞아 되돌아오는 궤적이 생기고, 그만큼 유효 반응 속도가 깎인다. 그래서 용매 점성이 높은 조건에서 반응 속도가 지수 인자 예측보다 느려지는 현상 — 용매 마찰 효과 — 이 관측되며, 이 편차 자체가 크라머스 이론의 실험적 증거로 쓰였다.

거꾸로, 반응좌표를 잘못 고르면 인위적인 재교차가 늘어 κ\kappa가 급격히 작아진다. 자유에너지 장벽은 반응좌표를 어떻게 잡아도 비슷하게 나오는 반면 전인자는 민감하므로, κ\kappa는 사실상 반응좌표 품질의 지표로도 읽힌다.2

4. 수치적으로 구하는 방법[편집]

  1. 직접 랑주뱅 적분. 우물에서 시작해 장벽을 넘을 때까지의 시간을 수천 번 모아 히스토그램을 만든다. 탈출 시간은 지수분포를 따르고, 그 평균의 역수가 속도다. ΔU/kBT=5\Delta U/k_BT = 5 정도면 현실적이지만, 10을 넘으면 기다리는 게 불가능해진다 — 이게 희귀 사건 문제의 원형이다.

  2. 평균 첫 통과 시간의 해석적 적분. 과감쇠 1차원에서는 평균 첫 통과 시간이 이중적분으로 정확히 풀린다.

    ⟨t⟩=1D∫ab ⁣eβU(y)[∫−∞y ⁣e−βU(z) dz]dy\langle t \rangle = \frac{1}{D}\int_{a}^{b}\! e^{\beta U(y)} \left[\int_{-\infty}^{y}\! e^{-\beta U(z)}\,dz\right] dy

    이 적분을 안장점 근사로 처리하면 크라머스 식이 그대로 나온다. 수치적으로는 구적법으로 바로 계산할 수 있어, 시뮬레이션 코드를 검증하는 기준 답으로 쓰기 좋다.

  3. 고급 샘플링. 자유에너지 장벽이 높으면 우산 샘플링·메타동역학으로 ΔU\Delta U를, 전이 경로 샘플링·포워드 플럭스 샘플링으로 전인자와 κ\kappa를 따로 구한다. 분자동역학 실무에서 반응 속도를 보고할 때 거의 표준 조합이다.

검증 팁: 랑주뱅 적분기의 시간 간격이 크면 장벽 정점 근처에서 궤적이 튀어 탈출률이 과대평가된다. Δt\Delta t를 절반으로 줄여 속도가 변하지 않는지 보는 것이 첫 점검이고, 그다음이 열욕 온도 확인(⟨v2⟩=kBT/m\langle v^2\rangle = k_BT/m)이다.

5. 응용[편집]

  • 화학반응 속도론. 용매 내 이성질화, 리간드 결합·해리. 아이링 식의 전인자에 용매 의존성을 넣는 물리적 근거가 된다.
  • 단백질 접힘. 접힘 속도의 상한이 확산 한계로 정해지며, 이걸 크라머스 꼴로 다룬다. “접힘 깔때기”에서 유효 확산계수를 재는 실험이 이 틀에 기반한다.
  • 자성 입자의 자화 반전. 네엘-브라운 완화 시간 τ=τ0eKV/kBT\tau = \tau_0 e^{KV/k_BT} 이 그대로 크라머스 문제이고, 하드디스크의 초상자성 한계와 기록 밀도 상한을 결정한다.
  • 조지프슨 접합. 위상 입자가 씻겨진 세탁판 퍼텐셜에서 탈출하는 것이 스위칭 전류의 통계적 분포를 만든다. 저온에서는 열 탈출이 양자 터널링으로 넘어간다.
  • 확률 공명. 주기적 구동에 탈출률이 동기화되면 잡음이 오히려 신호를 키운다. 크라머스 속도가 시간에 따라 변하는 경우로 다뤄진다.

6. 한계[편집]

크라머스 이론은 (i) 장벽이 kBTk_BT보다 충분히 높고, (ii) 마찰이 시간 상관 없는 백색잡음이며, (iii) 반응좌표가 1차원일 때 성립한다. 셋 다 현실에서 자주 깨진다. 마찰에 메모리가 있으면 그로테-헤인스-폴락 이론으로 확장되고(유효 마찰이 장벽 통과 시간 스케일에서 평가된다), 다차원에서는 랭거의 안장점 이론이 필요하다. 장벽이 kBTk_BT 수준으로 낮아지면 “탈출”과 “진동”의 구분 자체가 무의미해져서, 속도상수라는 개념을 포기하고 전체 확률분포를 다루는 편이 낫다.

7. 관련 문서[편집]

8. Footnotes[편집]

  1. 크라머스 논문은 1940년 Physica에 실렸는데 30년 가까이 거의 인용되지 않았다. 용매 효과를 재는 실험 기술이 따라온 1970~80년대에 와서 폭발적으로 인용된다. 좋은 이론이 시대를 잘못 만난 대표 사례. ↩

  2. κ\kappa가 10−310^{-3}처럼 나오면 마찰이 큰 게 아니라 반응좌표를 잘못 잡았을 확률이 훨씬 높다. 이 경우 답은 “샘플을 더 모으자”가 아니라 “좌표를 다시 고르자”다. ↩