一种计算扩散存在下缪子弛豫函数的数值方法

我们在强碰撞近似框架内，提出了一种准确且高效的方法，用于计算缪子处局部场随机波动的影响，例如在缪子扩散的情况下。该方法基于在离散时间基上对马尔可夫过程的重构，导出了缪子极化函数的求和方程，并通过离散傅里叶变换求解。后者在形式上类似于（尽管不完全相同于）原始连续时间模型的积分方程，后者通过拉普拉斯变换求解。使用实际案例参数值，发现离散时间强碰撞模型的解即使在粗粒度时间采样下，也能以极高的精度逼近连续时间解。通过快速傅里叶变换算法进行的计算非常高效，即使在慢速计算机上也适合实验数据的实时拟合。