分步传播边界层问题：高斯光束长距离传播结果不符

可能的问题原因 #

• 分步傅里叶方法的累积数值误差：短距离传播时误差可忽略，但长距离多次迭代后，一阶分步傅里叶（SSFM）的每步O(Δz²)误差会持续累积，最终导致结果偏离解析解。
• 吸收边界层的参数适配不足：即便滤波器本身没问题，固定参数的吸收层可能无法匹配长距离后光束的扩散程度——光束边缘能量占比变化，固定强度的吸收要么抑制不足导致边界反射，要么过度吸收损耗有效能量。
• 网格采样的混叠效应：长距离传播后光束空间频率成分变丰富，原网格采样率若低于Nyquist频率，傅里叶变换时会出现混叠，误差随传播步数不断放大。
• 解析解与数值模拟的前提偏差：要确认解析解是否基于无限空间假设，而数值模拟的有限网格在长距离下的截断效应未被完全消除；或者初始场与解析解的假设分布存在细微偏差，长距离后被放大。

对应的解决方法 #

• 升级分步傅里叶拆分阶数：把一阶SSFM换成二阶对称拆分方案（先做半步傅里叶域传播，再做空域相位调制，最后完成后半步傅里叶域传播），将每步误差降到O(Δz⁴)，大幅减少长距离累积误差。
• 动态调整边界层参数：根据当前传播距离下的光束尺寸，动态提升边界层的吸收斜率，或者改用高阶渐变窗（如布莱克曼窗）替代简单的线性吸收，避免边界处的突变反射。
• 优化网格采样参数：预先根据最大传播距离计算光束的预期扩散范围，调整网格点数和采样间隔，确保Nyquist频率覆盖光束的最高空间频率；也可以在传播中监测频谱，若出现混叠及时调整网格。
• 对齐解析解与模拟的前提条件：严格复刻解析解的初始场（比如高斯光束的振幅、相位分布），对比数值化后的解析解结果，逐点定位误差区域；若解析解忽略了某些效应（如色散），需在模拟中同步屏蔽对应模块。
• 添加能量监测：每步传播后计算光束总能量，若能量波动超过阈值，说明边界层或数值精度存在问题——能量损失过大可能是吸收过度，能量增加则大概率是边界反射未被抑制。

内容的提问来源于stack exchange，提问作者Henrique Guerra