[论文解读] Mean square displacement and instantaneous diffusion coefficient of charged particles in stochastic motion

Gabriela Raluca Mocanu|arXiv (Cornell University)|Jun 25, 2019

High-Energy Particle Collisions Research被引用 4

一句话总结

本文提出一种数值方法，用于计算带电粒子在各种天体物理构型中经历随机运动时的均方位移（MSD）和瞬时扩散系数（D(t)）。通过数值求解随机微分方程（SDE）——采用欧拉法、龙格-库塔法及混合方法——发现D(t)在中间时间区间表现出强烈的非单调行为，其峰值可比长时间极限高出一个数量级，从而能够对等离子体环境中的物理条件进行差异诊断。

ABSTRACT

The mean square displacement and instantaneous diffusion coefficient for different configurations of charged particles in stochastic motion are calculated by numerically solving the associated equations of motion. The method is suitable for obtaining accurate descriptions of diffusion in both intermediate and long time regimes. It is also appropriate for studying a variety of astrophysical configurations since it may incorporate microscopic physics that analytical methods cannot cope with. The results show that, in the intermediary time regime, the diffusion coefficient has an irregular behavior, which can be described in terms of the complex interplay appearing between the physical parameters describing the configuration. The main conclusion is that such an approach may serve at differential diagnosis of different astrophysical configurations.

研究动机与目标

• 开发一种数值框架，用于计算在复杂非平衡天体物理等离子体系统中瞬时扩散系数D(t)的值，这些系统中解析解难以获得。
• 研究物理参数（如电场强度、阻尼力、谐振势和磁场）如何影响D(t)和MSD的时间演化。
• 探索D(t)的中间时间行为，该行为在解析处理中常被忽略，但可能对瞬态天体物理事件携带关键诊断信息。
• 评估D(t)作为诊断工具的潜力，依据其独特的扩散特征区分不同天体物理构型。
• 证明SDE的数值积分能够捕捉复杂、非平凡的动力学行为——包括瞬时扩散增强——这些行为超出了标准解析近似的范围。

提出的方法

• 数值求解带电粒子运动的随机微分方程（SDE），表达式为 $ \frac{d^2x}{dt^2} = a_{\text{det}} + a_{\text{stoch}} $，其中 $ a_{\text{det}} $ 为确定性加速度，$ a_{\text{stoch}} $ 为零均值白噪声加速度，满足 $ \langle \xi_A(t_1)\xi_A(t_2) \rangle = A\delta(t_1 - t_2) $。
• 将粒子运动分解为平均分量 $ \bar{x}(t) $ 和涨落分量 $ \delta x(t) $，初始条件为 $ \delta x(0) = \delta v(0) = 0 $，以相对于平均轨迹计算扩散。
• 根据不同构型采用不同的数值积分器：常数电场（A）和具有常数阻尼的谐振势（B）使用欧拉法，具有记忆依赖阻尼的构型（C）使用二阶龙格-库塔法，磁场构型（D）采用混合欧拉-莱芒斯与考夫曼方法。
• 将所有方程无量纲化，以提高数值稳定性，并实现跨物理尺度的比较。
• 从系综平均轨迹计算均方位移（MSD） $ \langle x^2(t) \rangle $，然后通过 $ D(t) = \frac{1}{2} \partial_t \langle x^2(t) \rangle $ 推导瞬时扩散系数。
• 分析D(t)的时间演化，特别是其中间时间区间的特征及其长时间极限 $ D = \lim_{t \to \infty} D(t) $，涵盖四种不同的物理构型。

实验结果

研究问题

• RQ1 在不同物理条件下，带电粒子系统在随机运动中，瞬时扩散系数D(t)如何随时间演化？
• RQ2 D(t)的中间时间行为与长时间极限的偏离程度如何？哪些物理参数控制这些偏离？
• RQ3 D(t)在中间时间区间的形状和峰值幅度能否作为区分不同天体物理等离子体构型的诊断特征？
• RQ4 随机加速度与确定性力（如电场、谐振势、阻尼力）之比如何影响瞬时扩散动力学？
• RQ5 与马尔可夫过程或无磁场情况相比，阻尼力的记忆效应或磁场的存在如何改变D(t)的时间演化？

主要发现

• 在中间时间区间，瞬时扩散系数D(t)表现出非平凡且不规则的行为，其峰值可比长时间极限D高出一个数量级。
• 对于具有记忆依赖阻尼的谐振势中的粒子（情况C），最大D(t)及其达到时间均随 $ W^2 $ 增大而减小，表明势阱束缚作用增强。
• 在固定 $ W^2 $ 条件下，最大D(t)及其峰值时间也随 $ \bar{\alpha} $ 增大而减小，证实更高的阻尼力会抑制瞬时扩散。
• 在恒定磁场中（情况D），三维均方位移满足线性时间依赖关系（$ \langle q^2(n) \rangle_0 \sim n $），这是由于沿磁场方向运动不受限，而垂直方向运动被限制。
• 长时间极限D始终为常数且有限，但D(t)的中间行为对构型的物理参数高度敏感，因而具备差异诊断潜力。
• 不同参数集下D(t)的曲线具有显著差异且互不重叠，证实时间分辨的扩散谱可用于推断天体物理等离子体中潜在的物理条件。