本文介绍了模拟光在均匀介质中传播的四种快速而严格的方法。结果表明，在自由空间传播中，对光滑强相位项的解析处理在减少计算量方面是非常有效的。因此，在不限制快速傅里叶变换算法应用的情况下，我们重新设计了平面波角谱（SPW）算子来处理线性、球形和一般光滑相位项。特别是对于非傍轴场传播，所提出的技术可以显著地减少所需的采样点数量。数值结果表明了新方法的有效性和准确性。
一.文章介绍
光学建模与设计是研究与开发中极其重要的一部分。由于人们对高质量光学系统（包括衍射光学和微光学、散射物体和部分相干源）的需求日益增加，基于几何光学和物理光学相结合的模拟方法，即场追迹变得非常重要。这种模拟技术的一个重要部分是谐波场在均匀介质中的传播。然而，能够快速、准确地模拟一般光场在自由空间中的传播仍然是一项具有挑战性的任务。常用的算法只能做到快速或者只是准确。
在本文中，我们没有进一步的物理近似，介绍了四种新的算法，基于平面波（SPW）算子的角谱，有效地计算包含平滑但强相位项的非傍轴矢量光场的传播。根据光滑相位项的形状，可以使用不同的传播算子。它们的共同点是避免了光滑相位项指数函数的采样。相反，平滑相位项是解析处理的，只需对残差进行采样即可执行传播操作；因此，称为半解析传播技术。
首先，在第二节中我们给出一个问题的描述并引入数学符号。然后，在第3节中，我们考虑了一个球面相位项，Mansuripur[6]为此引入了一种严格的技术，称为使用快速傅里叶变换（FFT）的扩展菲涅耳衍射积分。在本节中，通过应用Van der Avoort等人最初使用的数值合适的抛物线拟合技术改进了该概念。在另一种情况下[7]，详细讨论了扩展菲涅耳算子在数值上可行的参数空间。此外，我们还介绍了扩展的菲涅耳算符的快速反演方法，用于快速计算非傍轴场到焦点区域的传播。
在第四节中，我们描述了一个用于光场快速传播的半解析SPW算子，它包含一个光滑的线性相位项。该方法基于线性相位项和横向偏移量的解析处理。之后，我们将这两种技术结合起来，得到了一个数值有效的半解析SPW算子，它能够同时解析地处理线性和球形相位项。
最后，在第6节中，我们通过将光场分解成具有平滑线性相位项的子光场，将半解析SPW算子概念推广到平滑相位的通用形状。在目标平面上，所有传播子光场被相干地相加，其中解析已知的平滑线性相位项以数值有效的方式使用第7节中介绍的逆抛物面分解技术（PDT）进行处理。数值结果证明了新的传播方法的有效性和准确性。所有的模拟都是用光学软件VirtualLab完成的。
二.均匀介质中的场追迹
在光场追迹法中，光在线性、均匀和各向同性介质中快速而精确的传播是由谐波场的概念处理的。结果表明，任何电磁场都可以分解为一组谐波场[8,9]。在空间频率域中，以特定角频率ω0振荡的单次谐波场定义为
 (1)
用位置向量  和角频率ω分别表示。请注意，下列理论是完全矢量的，因为在式（1）中，谐波场分量代表三个电场分量和三个磁场分量，由于计算效率高，常用的谐波传播技术基于FFT算法[10]。一种严格的传播技术是SPW算子[5]，其中各谐波场分量的复振幅在与传播方向正交的平面边界上，通过傅里叶变换（FT）分解成一组平面波
 (2)
 是初始平面边界上的横向位置向量，是  对应的空间频率矢量。用  表示的平面波通过与传播因子相乘，在距离z上传播
 (3)
 表示折射率为n的均匀介质中的波数，c为光的真空速度。最后，利用逆傅里叶变换将所有平面波叠加，从而得到SPW传播算子，
 (4)
从物理角度来看，SPW算子对任何传播距离z和任何空间频率矢量k[5]都是有效的。然而，对于长的传播距离，采样公式（4）的数值工作量太大。对于非傍轴场，它包含高频分量，数值工作量将变得更高。图1示意性地示出了由于快速增长的数值工作量而导致的SPW算子的有限范围。
一篇文献中报道了两个近似公式（4），以克服这一局限性。在这两种情况下，即菲涅耳和远场积分，使用近似来分离球面相位项与积分的数值计算，并且球面相位项是被解析地处理的。这就大大减少了数值计算的工作量。然而，由于这些近似，两种解决方案的适用性受到限制。菲涅耳积分[11]使用空间频率分量的泰勒展开。用这种方法，将式（3）的球面相位函数替换为抛物线相位函数，从而得到式（4）中逆傅里叶变换积分的半解析解。如图1所示，该概念仅适用于具有低空间频率的傍轴场。对于长传播距离z，可以应用远场积分[5]，其中公式（4）的逆傅里叶变换积分可以用固定相位的想法来半解析地求解。
图1：场追迹中常用的基于FFT的传播算子的有效范围的示意图。灰色区域表示在数值上可行的参数空间。严格的SPW算子在长传播距离或大空间频率下很难突破数值工作量。由于物理近似的性质，远场和菲涅耳传播算子具有有限的有效范围。
由于普通的自由空间传播算子使用范围的限制，先进的传播技术的发展对于非傍轴场追迹是必不可少的。这些技术应覆盖图1中剩余的白参数空间，并在精度和数值计算之间取得最佳折衷。文献[2,3,12-14]中有几种方法。然而，这个问题的一般解决办法仍不清楚。在实践中，非傍轴谐波场振幅出现了采样问题，其中部分包含强而光滑的相位项  数学表达式为
 (5)
|
