[論文レビュー] Monte Carlo Approaches for Simulating a Particle at a Diffusivity Interface and the "Ito--Stratonovich Dilemma"

Hendrick W. de Haan, Mykyta V. Chubynsky+1 more|arXiv (Cornell University)|Aug 24, 2012

Stochastic processes and statistical mechanics参考文献 8被引用数 8

ひとこと要約

本稿では、不連続な拡散係数を持つ界面における粒子拡散をモデル化する2つのモンテカルロシミュレーション手法を導入し、確率積分の解釈（伊藤、ストラトニヴィチ、等温）を柔軟に可能にすることで「伊藤＝ストラトニヴィチのジレンマ」を解決する。粒子の分布および初到達ダイナミクスは、解釈の種別によって顕著に異なることが示され、主な発見として、非一様な粒子滞在時間、非伊藤解釈における低粘度壁へのバイアス、および解釈依存の指数を有するパワーロー法則に従う初到達時間分布が得られた。

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.

研究の動機と目的

• 不連続な拡散係数を有する不均一媒質における過減衰ランジュバン力学の「伊藤＝ストラトニヴィチのジレンマ」を解決すること。
• このような系において、確率積分の解釈（伊藤、ストラトニヴィチ、等温）を柔軟に選択可能な2つのモンテカルロシミュレーション手法を開発すること。
• 異なる微積分解釈が、粒子分布、初到達時間、および界面の各側での時間平均滞在時間に与える影響を調査すること。
• MC手法の正確性を理論的予測および既存のブラウン運動／ランジュバン力学シミュレーションと照合すること。

提案手法

• 1次元区間内で異なる拡散係数を有する2領域を持つ系において、離散時間ステップを用いて粒子運動をシミュレートする2つの異なるモンテカルロアルゴリズムを開発する。
• 界面および各領域内での遷移確率を変更することで、異なる確率微積分解釈（伊藤、ストラトニヴィチ、等温）を実装する。
• 定常状態における粒子分布を調査するため反射境界条件を、初到達時間を分析するため吸収境界条件を用いる。
• 粒子の軌道を界面への訪問間のセグメントに分解し、各側での平均滞在時間を計算する。この際、セグメントの開始に起因するバイアスを考慮する。
• 平均初到達時間の解析的解と、格子ベースの理論的導出（例：$ n = d^2/4 $）と照合することで、精度を検証する。
• 軌道分解を用いて、各側での平均滞在時間を計算し、条件付きおよび無条件のMFPTを区別する。

実験結果

リサーチクエスチョン

• RQ1 確率微積分の解釈（伊藤、ストラトニヴィチ、等温）の選択が、拡散係数界面における定常状態の粒子分布に与える影響は何か？
• RQ2 中心に界面を持つ1次元区間において、異なる微積分解釈が壁への平均初到達時間に与える影響は何か？
• RQ3 両壁間の初到達時間分布はどのように異なり、パワーロー法則的挙動を示すか？その指数は微積分の種別に依存するか？
• RQ4 特定の側での平均滞在時間は、その壁に到達する確率と相関するか？この相関関係は、平均化の方法にどのように依存するか？
• RQ5 粒子の滞在時間分布は時間経過とともにどのように変化するか？また、異なる微積分法による一時的バイアスの原因は何か？

主な発見

• 非等温微積分法では、定常状態において界面で粒子分布に不連続性が生じるが、等温法では分布が一様である。
• 伊藤解釈では、どちらの壁にも到達する傾向に偏りがないが、他のすべての解釈では低粘度側へのバイアスが生じる。
• 低粘度側への平均初到達時間は、反対側の粘度が高くなるに従って増加するが、等温法では一定値に近づく。
• 粘度比が高くなると、低粘度側への初到達時間分布は著しく広がり、微積分法に依存する指数を有するパワーロー法則的依存性を示す。
• 特定の側での平均滞在時間は、その壁に到達する確率とは著しく異なる場合があり、この乖離は、軌道に沿った平均化の方法に強く依存する。
• 一部の微積分法では、初期段階において粒子が一時的に低粘度側を優遇するが、その後、運動が遅いために高粘度側を優越するよう変化する。