[논문 리뷰] Anomalous Stochastic Transport of Particles with Self-Reinforcement and Mittag–Leffler Distributed Rest Times

Daniel Han|arXiv (Cornell University)|2021. 11. 15.

Diffusion and Search Dynamics참고 문헌 65인용 수 11

한 줄 요약

이 논문은 자기강화와 미타그레플레르 분포를 가진 정지 시간을 갖는 지속적인 랜덤 워크 모델을 제안하며, 하이퍼볼릭 편미분방정식과 리만-리우빌 분수도수 도함수를 통해 모델링한다. 분석적·수치적으로 분석한 결과, 자기강화로 인해 중간 시간에 잠시 초산산운동이 발생하지만, 장기적으로는 무거운 尾 꼬리 분포를 가진 정지 시간으로 인해 초산산운동에서 초산산운동으로 회귀함을 확인하였으며, 이는 세포내 시스템에서 비정상 운동을 설명하는 생물학적으로 타당한 메커니즘을 제공한다.

ABSTRACT

We introduce a persistent random walk model for the stochastic transport of particles involving self-reinforcement and a rest state with Mittag–Leffler distributed residence times. The model involves a system of hyperbolic partial differential equations with a non-local switching term described by the Riemann–Liouville derivative. From Monte Carlo simulations, we found that this model generates superdiffusion at intermediate times but reverts to subdiffusion in the long time asymptotic limit. To confirm this result, we derived the equation for the second moment and find that it is subdiffusive in the long time limit. Analyses of two simpler models are also included, which demonstrate the dominance of the Mittag–Leffler rest state leading to subdiffusion. The observation that transient superdiffusion occurs in an eventually subdiffusive system is a useful feature for applications in stochastic biological transport.

연구 동기 및 목표

• 자기강화된 방향성과 비율분포를 가진 정지 시간을 갖는 생물학적 시스템에서의 비정상적인 확률적 운동을 모델링하기 위해.
• 장기적인 운동 행동을 결정짓는 자기강화와 비마르코프 정지 상태 간의 상호작용을 조사하기 위해.
• 초산산운동에서 초산산운동으로의 전이를 기술하는 분수도수 편미분방정식계를 유도하고 검증하기 위해.
• 세포내 소기관 이동에서 관찰되는 일시적인 초산산운동과 일치하는 메커니즘을 제공하기 위해.

제안 방법

• 활성 상태(±ν)와 정지 상태를 갖는 세 상태의 스토케스틱 모델을 제안하며, 리만-리우빌 분수도수 도함수를 통한 비국소적 전환을 갖는 하이퍼볼릭 편미분방정식으로 제어된다.
• 자기강화는 시간에 따라 변하는 전이 확률 r± = 1/3 ± α₀x/(2νt)를 통해 도입되며, 여기서 α₀는 자기강화 계수이다.
• 정지 시간은 리만-리우빌 분수도수 도함수 D₁₋ᵝₜ를 통해 미타그레플레르 분포를 가진 체류 시간으로 모델링된다.
• 라플라스 공간에서 두 번째 모멘트 방정식을 유도하고, 역라플라스 변환을 수행하여 장기적 스케일링을 분석한다.
• 분석적 예측을 검증하기 위해 JIT-컴파일된 파이썬을 사용한 몬테카를로 시뮬레이션을 병렬화하여 수행한다.
• 입자 위치의 확률밀도함수(PDF)를 분석하여 일시적인 이중모드 비대칭성과 장기적 라플라스 감쇠를 확인한다.

실험 결과

연구 질문

• RQ1 자기강화는 미타그레플레르 분포를 가진 정지 시간을 갖는 지속적인 랜덤 워크에서 초산산운동을 유도하는가?
• RQ2 무거운 꼬리 분포를 가진 정지 시간이 존재할 경우, 두 번째 모멘트의 장기적 점근적 행동은 어떠한가?
• RQ3 비마르코프 시스템에서 자기강화가 존재할 때, 일시적인 초산산운동과 최종적인 초산산운동이 공존할 수 있는가?
• RQ4 거듭제곱 법칙 정지 분포의 특성 시간 스케일 τ₀은 일시적인 초산산운동의 발생에 어떤 영향을 미치는가?
• RQ5 장기적 한계에서 시스템은 어느 정도 분수도수 확산 방정식으로 근사될 수 있는가?

주요 결과

• α₀ ≠ 0일 때 중간 시간에 두 번째 모멘트가 µ₂(t) ∼ t³|α₀|로 스케일링되며, 이는 일시적인 초산산운동을 나타낸다.
• 장기적 한계에서 두 번째 모멘트는 µ₂(t) ∼ tᵝ이며, 0 < β < 1로 스케일링되며, 이는 미타그레플레르 분포를 가진 정지 시간으로 인한 초산산운동임을 확인한다.
• 분수도수 확산 한계는 두 번째 예제에서만 존재하며, Dβ = ν²/(λ²Γ(β+1)τ₀ᵝ)를 갖는 분수도수 확산 방정식을 도출한다.
• 작은 τ₀일 경우, 중간 시간에 PDF는 비대칭적이며 이중모드이며, 자기강화로 인한 방향성 유지가 나타난다.
• 큰 τ₀일 경우, 장기적 한계에서 PDF는 라플라스 형태로 복귀하며, 이는 초산산운동 행동과 일치한다.
• 비마르코프 정지 상태는 자기강화 강도(α₀)에 관계없이 장기적 한계에서 초산산운동을 완전히 억제한다.