조화 포텐셜 안에서 d차원 런‑텀블 입자의 정확한 정상상태 해

초록

이 논문은 등방성 조화 포텐셜 V(r)=μr²/2에 갇힌 d차원 런‑텀블 입자(RTP)의 비평형 정상분포를 정확히 구한다. 회전 대칭성을 이용해 한 좌표의 마진 p_X(x)만을 분석하면 반경 분포 p_R(r)와 전체 다변량 밀도를 적분 변환으로 얻을 수 있다. 1차원에서 임의의 속도 분포 W(v)를 갖는 일반화된 RTP에 대해 Kesten‑type 재귀와 스틱‑브레이킹(Dirichlet) 과정으로 정밀한 폐쇄형 해를 도출하고, 이를 d=1,2,3에 적용한다. d=1,2에서는 p_R(r)이 베타 분포가 되며, d=3에서는 새로운 폐쇄형 식이 얻어진다. 또한 템플링 비율 α=γ/μ에 따른 “형태 전이”를 밝혀내고, 열확산 D>0이 포함될 경우 Gaussian 컨볼루션으로 정규화된 분포를 제시한다. 결과는 수치 시뮬레이션과 완벽히 일치한다.

상세 분석

본 연구는 활발히 움직이는 미생물이나 인공 스위머를 모델링하는 런‑텀블 입자(RTP)의 정적 특성을, 외부 조화 포텐셜이라는 가장 단순한 구속 환경에서 완전 해석적으로 풀어낸 점이 혁신적이다. 핵심 아이디어는 시스템이 회전 대칭을 갖기 때문에 d차원 문제를 단일 좌표 마진 p_X(x)로 축소할 수 있다는 사실이다. 이를 위해 저자들은 1차원에서 속도 재설정 분포 W(v)를 임의로 지정할 수 있는 ‘일반화된 RTP’를 정의하고, Kesten 재귀 관계를 이용해 위치 X의 고정점 분포를 무한 곱 형태로 전개한다. 이 곱은 스틱‑브레이킹(Dirichlet) 과정과 동일시될 수 있는데, 즉 X는 무한히 작은 조각들의 가중합으로 표현된다. 이 접근법은 기존에 비선형 비국소 Fokker‑Planck 방정식을 직접 푸는 것이 거의 불가능했던 문제를, 확률론적 고정점 이론과 베르누이·베타 함수의 조합으로 해결한다는 점에서 방법론적 의의가 크다.

특히, 속도 분포 W(v)를 등방성 RTP의 투영 분포 W_proj(v)로 지정하면, d차원에서 각 좌표는 동일한 1차원 과정에 의해 기술된다. d=2에서는 W_proj(v) 가 아크사인 법칙이 되고, 결과적으로 반경 분포 p_R(r) 가 베타 분포 형태(p_R∝r^{2α−1}(1−(μr/v₀)²)^{α−1})를 갖는다. 이는 기존에 모멘트만으로 추정되던 결과를 완전한 확률밀도 함수로 확장한 것이다. d=3에서는 투영 분포가 균등이므로, 베타 형태가 깨지고 보다 복잡한 적분식(p_R∝∫₀^{A(r)}(A(r)−u)^{α−1}f′(u)du)으로 나타난다. 여기서 A(r)와 f(u)는 α와 μ, v₀에 의존하는 명시적 함수이며, α=1에서 특수한 단순식(아크탄하와 삼각함수 조합)으로 축소된다.

형태 전이는 파라미터 α=γ/μ, 즉 평균 런 시간과 포텐셜 완화 시간의 비율에 의해 지배된다. α<1(고 지속성)에서는 확률 질량이 turning surface r=v₀/μ 근처에 집중되어 경계에서 적분 가능하지만 발산하는 특이점을 보이며, α>1(저 지속성)에서는 중심부에 집중된 단일 피크를 가진다. 이러한 전이는 d=1,2에서 베타 분포의 지수 변화로, d=3에서는 내부·외부 피크가 공존하는 복합 전이로 나타난다.

열확산 D>0을 포함하면, 선형성 때문에 정상분포는 D=0 경우와 Gaussian(평균 0, 분산 D/μ) 커널의 컨볼루션으로 정확히 표현된다. 이 과정은 turning surface에서의 발산을 부드럽게 정규화하고, 포텐셜 경계 밖에 가우시안 꼬리를 생성한다. 저자들은 무차원 스케일링 변수 θ=2μD/v₀²를 도입해, θ→0(활동 지배)와 θ→∞(열 확산 지배) 사이의 교차를 정량화한다. 특히 d=3에서는 θ가 임계값을 넘을 때 전역 최대가 r=0에서 껍질 근처로 불연속적으로 이동하는 현상이 관측된다.

마지막으로, N-상태 RTP(속도값이 이산 집합 {v_i}에 확률 p_i로 선택)에도 동일한 프레임워크를 적용해, p_X(x)를 Dirichlet 가중합의 밀도로 정확히 구한다. 이는 실험적으로 제한된 속도 레벨을 갖는 인공 마이크로스위머 모델링에 바로 활용될 수 있다. 전체적으로, 확률적 고정점 이론, Dirichlet 과정, 그리고 물리적 파라미터 매핑을 결합해 고차원 비평형 시스템의 정확 해를 제공한 점이 가장 큰 공헌이다.