[논문 리뷰] Monte Carlo Approaches for Simulating a Particle at a Diffusivity Interface and the "Ito--Stratonovich Dilemma"

Hendrick W. de Haan, Mykyta V. Chubynsky외 1명|arXiv (Cornell University)|2012. 08. 24.

Stochastic processes and statistical mechanics참고 문헌 8인용 수 8

한 줄 요약

이 논문은 확산계수 인터페이스에서 입자 확산을 모의하기 위해 몬테카를로 시뮬레이션 방법 두 가지를 제안하며, 스토케스틱 미적분학의 해석 유연성을 제공함으로써 '이토–스트라토노비치 딜레마'를 해결한다 (이토, 스트라토노비치, 등온). 입자 분포 및 제1도달 시간 역학은 다양한 해석에 따라 크게 달라지며, 주요 발견으로는 비균일한 입자 체류 시간, 비이토 해석에서 낮은 점성 벽 측면으로의 편향, 그리고 해석에 따라 결정되는 지수를 가진 등급법적 제1도달 시간 분포가 있다.

ABSTRACT

The possibility of different interpretations of the stochastic term (or calculi) in the overdamped Langevin equation for the motion of a particle in an inhomogeneous medium is often referred to as the "Ito--Stratonovich dilemma," although there is, in fact, a continuum of choices. We introduce two Monte Carlo (MC) simulation approaches for studying such systems, with both approaches giving the choice between different calculi (in particular, Ito, Stratonovich, and "isothermal"). To demonstrate these approaches, we study the diffusion on a 1D interval of a particle released at an interface (in the middle of the system) between two media where this particle has different diffusivities (e.g., two fluids with different viscosities). For reflecting boundary conditions at the ends of the interval, a discontinuity at the interface in the stationary-state particle distribution is found, except for the isothermal case, as expected. We also study the first-passage problem using absorbing boundary conditions. Good agreement is found when comparing the MC approaches to theory as well as Brownian and Langevin dynamics simulations. Also, the results themselves turn out to be interesting. For instance: 1) for some calculi, there are more particles on the low-viscosity side (LVS) at earlier times and then fewer particles at later times; 2) there is no preference to end up on a particular wall for the Ito variant, but a bias towards the wall on the LVS in all other cases; 3) the mean first-passage time to the wall on the LVS grows as the viscosity on the other side is increased, except for the isothermal case where it approaches a constant; 4) when the viscosity ratio is high, the first-passage-time distribution for the wall on the LVS is much broader than for the other wall, with a power-law dependence of the former in a certain time interval whose exponent depends on the calculus.

연구 동기 및 목표

• 비연속적인 확산계수를 가진 비균질 매질에서 과다된 랑주방 정역학의 '이토–스트라토노비치 딜레마'를 해결하기 위해.
• 이토, 스트라토노비치, 등온 해석을 자유롭게 선택할 수 있도록 하는 두 가지 몬테카를로 시뮬레이션 접근법을 개발하기 위해.
• 다양한 미적분학 해석이 입자 분포, 제1도달 시간, 인터페이스 양측에서의 시간 평균 체류 시간에 미치는 영향을 조사하기 위해.
• MC 방법의 정확성을 검증하기 위해 이론적 예측 및 기존의 브라운운동/랑주방 역학 시뮬레이션과 비교하기 위해.

제안 방법

• 이산 시간 단계를 사용하여 서로 다른 확산계수를 가진 두 영역이 있는 1차원 간격에서 입자 운동을 시뮬레이션하는 두 가지 별개의 몬테카를로 알고리즘을 개발한다.
• 인터페이스 및 각 영역 내에서 전이 확률를 수정함으로써 이토, 스트라토노비치, 등온 해석과 같은 다양한 스토케스틱 미적분학 해석을 구현한다.
• 정적 상태에서 입자 분포를 연구하기 위해 반사 경계 조건을 사용하고, 제1도달 시간을 분석하기 위해 흡수 경계 조건을 사용한다.
• 입자 궤적을 인터페이스에 대한 방문 간의 세그먼트로 분해하여 각 측면에서 평균 체류 시간을 계산하며, 세그먼트 시작 시의 편향을 고려한다.
• 평균 제1도달 시간에 대한 해석적 해와 격자 기반 이론 유도 결과(예: $ n = d^2/4 $)와의 비교를 통해 정확도를 검증한다.
• 궤적 분해를 통해 각 측면에서의 평균 체류 시간을 계산하며, 조건부 및 무조건적 평균 제1도달 시간(MFPT)을 구분한다.

실험 결과

연구 질문

• RQ1 스토케스틱 미적분학 해석(이토, 스트라토노비치, 등온)의 선택이 확산계수 인터페이스에서 정적 상태 입자 분포에 어떤 영향을 미치는가?
• RQ2 중앙 인터페이스가 있는 1차원 간격에서 벽으로의 평균 제1도달 시간에 다양한 미적분학 해석 해석이 어떤 영향을 미치는가?
• RQ3 두 벽으로의 제1도달 시간 분포는 어떻게 다를지, 그리고 계산에 따라 결정되는 지수를 가진 등급법적 행동을 보이는가?
• RQ4 특정 측면에서의 평균 체류 시간이 그 벽에 도달할 확률과 상관이 있는가, 그리고 이 관계는 평균화 방법에 따라 어떻게 달라지는가?
• RQ5 입자 체류 시간 분포는 시간이 지남에 따라 어떻게 변화하는가, 그리고 다양한 계산법에서 발생하는 일시적 편향의 원인은 무엇인가?

주요 결과

• 비등온 계산법에서는 정적 상태에서 인터페이스에서 입자 분포에 불연속성이 나타나지만, 등온 해석에서는 균일한 분포를 보인다.
• 이토 해석에서는 어느 벽으로 도달할지에 대한 선호도가 없지만, 그 외 모든 경우에서 낮은 점성 측면으로의 편향이 나타난다.
• 저점성 벽으로의 평균 제1도달 시간은 반대 측면의 점성이 높아질수록 증가하지만, 등온 해석에서는 일정한 값으로 수렴한다.
• 점성 비율이 높을 경우, 저점성 벽으로의 제1도달 시간 분포는 상당히 넓어지며, 사용된 계산법에 따라 결정되는 지수를 가진 등급법적 의존성을 보인다.
• 특정 측면에서의 평균 체류 시간은 그 측면으로 도달할 확률과 다를 수 있으며, 이 격차는 궤적에 대한 평균화 방법에 따라 결정된다.
• 일부 계산법에서는 초기 시점에 낮은 점성 측면을 일시적으로 선호하는 입자 집단이 나타나지만, 후속적으로 느린 동역학으로 인해 고점성 측면을 선호하게 된다.