[论文解读] 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

一句话总结

本文提出了两种蒙特卡洛模拟方法，用于模拟在扩散系数界面处的粒子扩散，通过灵活选择随机微积分解释（伊藤、斯特拉托诺维奇、等温），解决了‘伊藤-斯特拉托诺维奇困境’。结果表明，不同解释下粒子分布与首次通过动力学存在显著差异，关键发现包括：非均匀的粒子驻留时间、非伊藤情况下向低黏度壁面的偏置，以及与微积分形式相关的幂律首次通过时间分布。

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.

研究动机与目标

• 解决在扩散系数不连续的非均匀介质中过阻尼朗之万动力学的‘伊藤-斯特拉托诺维奇困境’。
• 开发两种蒙特卡洛模拟方法，使此类系统中可灵活选择随机微积分解释（伊藤、斯特拉托诺维奇、等温）。
• 研究不同微积分解释对粒子分布、首次通过时间及界面两侧平均驻留时间的影响。
• 将蒙特卡洛方法与理论预测以及现有的布朗运动/朗之万动力学模拟结果进行验证。

提出的方法

• 开发两种不同的蒙特卡洛算法，模拟在一维区间内具有两个不同扩散系数区域的粒子运动，采用离散时间步长。
• 通过在界面处及各区域内修改转移概率，实现不同随机微积分解释（伊藤、斯特拉托诺维奇、等温）。
• 采用反射边界条件研究稳态粒子分布，采用吸收边界条件分析首次通过时间。
• 将粒子轨迹分解为多次访问界面之间的片段，以计算在每一侧的平均驻留时间，同时考虑片段起始时的偏置。
• 将蒙特卡洛结果与首次通过时间的解析解以及基于格点的理论推导（如 $ n = d^2/4 $）进行对比，以验证精度。
• 利用轨迹分解方法计算在每一侧的平均驻留时间，区分条件与非条件下的平均首次通过时间（MFPT）。

实验结果

研究问题

• RQ1 随机微积分的选择（伊藤、斯特拉托诺维奇、等温）如何影响扩散系数界面处的稳态粒子分布？
• RQ2 在具有中心界面的一维区间中，不同微积分解释对粒子到达壁面的平均首次通过时间有何影响？
• RQ3 首次通过时间分布在不同壁面之间有何差异？是否表现出幂律行为，且幂指数依赖于所用微积分？
• RQ4 在特定一侧的平均驻留时间是否与到达该壁面的概率相关？这种关系如何依赖于轨迹平均方法？
• RQ5 粒子驻留时间分布如何随时间演化？不同微积分形式导致的瞬态偏置现象由何引起？

主要发现

• 对于非等温微积分，粒子分布的稳态在界面处出现不连续，唯独在等温情况下分布为均匀。
• 在伊藤解释下，粒子到达任一壁面无偏好；但在其他所有情况下，均存在向低黏度侧的偏置。
• 粒子到达低黏度壁面的平均首次通过时间随对侧黏度升高而增加，但在等温情况下趋于恒定。
• 当黏度比值较高时，粒子到达低黏度壁面的首次通过时间分布显著变宽，并表现出幂律依赖关系，其指数取决于所用微积分形式。
• 在特定一侧的平均驻留时间可能与到达该壁面的概率显著不同，这种差异在很大程度上取决于轨迹平均方式。
• 对于某些微积分形式，粒子群体在早期时间短暂偏好低黏度侧，但随后因动力学更慢而转向偏好高黏度侧。