[论文解读] Simulation of particle systems interacting through hitting times

Vadim Kaushansky, Christoph Reisinger|arXiv (Cornell University)|May 29, 2018

Markov Chains and Monte Carlo Methods参考文献 10被引用 12

一句话总结

本文提出了一种新颖的欧拉型粒子方法，用于通过首次 hitting 时间模拟具有平均场相互作用的 McKean–Vlasov 方程，其中粒子在边界处吸收会引发反馈。该方法在时间步长上实现了 1/2 阶的收敛率，通过使用布朗桥和局部网格加密进一步提升至 1 阶，数值验证表明其有效性，并成功应用于具有爆破行为的系统性风险和神经网络模型。

ABSTRACT

We develop an Euler-type particle method for the simulation of a McKean--Vlasov equation arising from a mean-field model with positive feedback from hitting a boundary. Under assumptions on the parameters which ensure differentiable solutions, we establish convergence of order $1/2$ in the time step. Moreover, we give a modification of the scheme using Brownian bridges and local mesh refinement, which improves the order to $1$. We confirm our theoretical results with numerical tests and empirically investigate cases with blow-up.

研究动机与目标

• 开发一种针对具有边界 hitting 反馈的 McKean–Vlasov 方程的构造性数值格式，该类方程此前缺乏严格的模拟方法。
• 在确保解具有可微性的假设下，分析基于粒子的时间步进格式的收敛速率。
• 通过布朗桥插值和局部网格加密，将收敛阶从 1/2 提升至 1。
• 通过实验研究系统性风险和神经网络模型中爆破现象的行为。
• 通过大规模粒子系统进行数值实验，验证理论收敛速率。

提出的方法

• 提出一种时间半隐式欧拉型粒子格式，用于模拟漂移项依赖于粒子首次 hitting 到零的 hitting 时间密度的 McKean–Vlasov SDE。
• 该格式利用粒子路径对 hitting 时间分布 $ L_t $ 的经验估计来计算反馈漂移。
• 改进的格式引入布朗桥插值，以提升在吸收边界附近的路径近似精度。
• 在 hitting 时间附近自适应地应用局部网格加密，以增强关键区域的精度。
• 在初始分布和参数具有正则性的假设下分析收敛性，确保解在爆破时间前可微。
• 通过方差缩减和交换性论证，建立经验 hitting 时间测度在概率和 $ L^2 $-范数下的收敛性。

实验结果

研究问题

• RQ1 标准欧拉型粒子格式在具有 hitting 时间反馈的 McKean–Vlasov 方程中，其收敛速率如何？
• RQ2 布朗桥插值和局部网格加密是否能将收敛速率提升至超过 1/2？
• RQ3 当 hitting 时间分布出现不连续性或爆破时，该格式的行为如何？
• RQ4 粒子系统规模对经验 hitting 时间估计精度有何影响？
• RQ5 该方法能否成功应用于现实世界模型，如金融网络中的系统性风险或神经元群体动力学？

主要发现

• 在正则性假设下，标准欧拉型粒子格式在时间步长上实现了 $ 1/2 $ 阶的收敛率。
• 采用布朗桥和局部网格加密的改进格式将收敛阶提升至 1 阶。
• 数值测试验证了理论收敛速率，并展示了在捕捉爆破行为方面的鲁棒性。
• 当 $ N \to \infty $ 时，经验 hitting 时间测度 $ ilde{L}_t $ 以概率收敛至真实 hitting 时间分布。
• 随着 $ N \to \infty $，经验 hitting 时间估计量 $ \hat{L}_t^N $ 的方差趋于零，确保了一致性。
• 该方法成功捕捉到了 $ L_t $ 中不连续性的形成，对应于金融和神经网络应用中的系统性冲击。