d .A0(*k, 从光线到物理光学直到现在，所有的效应和量都能够在几何光学的框架中进行表示。下一步中，我们探索在两种介质间界面的能量效应，例如，一个透镜的表面。众所周知，在界面处，一部分光被反射回去因此会造成透射部分能量的损失。4%是空气和玻璃介质间透射能量损失的典型值。似乎我们可以直接将这个值对每条光线的作用考虑进去，进而减少在探测平面的探测能量。然而，在我们简单的将此损失包含在光线追迹中之前，我们应该考虑其在光线光学，即费马原理中的正确性。在介质间界面4%的能量损失符合费马原理吗？答案是否定的，由于此原理仅处理光程，因此我们无法在光线光学的框架中找到这种表面效应的合理解释。在各种光学教科书中，你可以找到菲涅尔方程的推导，其给出了能量透射率T（透射比）和能量反射率R（反射比）的数学表达式[1]。此推导考虑的是理想电磁场平面波穿过两种不同折射率介质间的理想平面界面。这个结论使用了在平面界面处电场和磁场分量是连续函数的事实。由此直接推导出菲涅尔方程。而理想平面波以及连续横向场分量则来自于麦克斯韦方程组[2]。与完全基于费马原理的几何光学相比，我们是基于麦克斯韦方程组来考虑物理光学的。因此，应该依据物理光学来解释在两种介质界面处光能量的损失，并将其附加到光线追迹路径，已经引出了一种联合了光线和物理光学的算法。然而，当将传统的光线追迹强行的与一种基于物理光学的效应，如表面处的菲涅尔效应或者光线透过光栅的传输，联合起来的时候，我们会面临一个典型且严重的问题。即，除了入射光角度很小的时候，菲涅尔效应都是与局部偏振相关的。因此，为了精确地包含菲涅尔效应，简单的光线不够，我们还需要其偏振信息。读者可能会问道，那么光线追迹软件是以什么标准来处理那些问题。事实上是其不可能精确地处理光滑的或者光栅类型的表面效应，并且也没有包含偏振效应的模型。
$}fA ;BP 传统光线追迹有许多限制，上文提到的仅是我们想要去克服的其中一种。我们所需要的是使用物理光学来表征光线。接下来我们会讨论从以几何光学为基础的传统光线过渡到在物理光学框架中使用的智能光线。
$(U|JR@ 几何场追迹 3I+pe;
Fig.2光线与光滑表面（左）、散射表面（中）和局部周期表面（右）的相互作用。使用局部平面波来表征光线
Fig.2演示了一束光线与曲面的相互作用。如果我们以电磁场平面波来解释一束局部光线且将界面局部地作为平面，那么可以局部地应用从物理光学所获得的理想平面波与理想平面界面相互作用的结果[3]。因此，我们开发出了一种用于这种局部电磁场平面波概念的算法，即几何场追迹[4]。此理论来自于《光学原理》[5]第三章所提出的结论，并且我们已经探讨了在几何场近似条件下使用几何场追迹[4]来求解麦克斯韦方程组。对于局部平面波，使用这种近似的方法来求解麦克斯韦方程组，可以给出波前为主导的场的空间演化区域的精确解。相反，若一个场的尺寸没有远大于波长，则其开始受衍射主导且不能使用几何场追迹来传播该场，而需要使用一个衍射场追迹的方法来进行传播。换句话说，一个场的散度可以用其波前的曲率（几何场追迹）或者通过横向和k空间维度间的不确定性原理的扩展（衍射场追迹）来控制。空间中的任意场都有其衍射的和几何的主导区域，对于不同的区域需要使用合适的建模技术。如Fig.3，对于传播经过其焦点的球面波。我们想着重的强调的是，我们已经应用了能够自动测试几何场近似有效性的算法以确保选取最合适的建模技术以用于传播的各个区域。
J7xmf,76w
Fig.3 透镜系统出瞳处一个截锥球面波通过其焦点传播的严格一维计算。下图是上图焦点区域的放大。球面波的f/#为20。尽管其是一个旁轴案例，但衍射主要在焦点区域
在实际中通过追迹智能光线来实现几何场追迹，当在VirtualLab Fusion软件中使用时，可以获得以下特性：
智能光线知道其位置处的所有电磁场信息，三个电场分量和三个磁场分量的振幅，相位以及偏振信息。
智能光线知道并记住光源平面中附近光线的波前。这通过一种合适的光线索引概念（波前索引）来实现。这个方法联合了不同的横向插值技术以用于所用的场量并将这些场量分配到一条光线上去。插值技术包含了样条曲线插值和质心坐标网格插值。
伴随智能光线的还有另外一个索引概念（空间相干索引），其能够给出智能光线间的相干模式和非相干模式以及其组合模式。这使得我们可以完成对部分相干光的建模。特殊的情况也被包含其中，如完全相干光和非相干光。
为了包含颜色，时间相干和超短脉冲，我们将频率分配给智能光线，使用一个索引概念来分辨稳态光和脉冲光的频率分布（频率索引）。
通过追迹智能光线我们获得麦克斯韦方程组几何场近似的解[4]。因此，智能光线将传统光线追迹推广，并用一种科学地可靠的方式将其结合到物理光学建模中。实际上，至关重要的是几何场追迹获得电磁场结果和传统光线追迹获得的结果，如点列图的速度是一样快的。
_L":Wux 高精度的模拟透镜系统可以主要依据几何场追迹，由于透镜系统旨在传输球面波，因此几何场追迹近似是有效的。如Fig.4所示。智能光线追迹提供了位于最后一个透镜后的所有电磁场信息。然而，如前面所讨论的（见Fig.3），位于焦面区域的场不满足几何场近似。因此，使用智能光线在最后一个透镜后获得的场必须使用衍射技术传播到焦点。在标准的光线追迹软件中，这样的一个过程或多或少的依赖于软件包中的点扩散函数（PSF）和调制传递函数（MTF）计算工具来精确地完成。场追迹为此以及其他类似的建模任务提供了一种快速和可靠的物理光学建模策略。
oV`sCr5%
Fig.4 上面一行显示的是x方向偏振光经过Fig.1中透镜系统的最后一个透镜后的电场分量E = (Ex,Ey,Ez)。结果包含了一个从x到y方向较弱的串扰和一个中等的z分量。此结果是通过几何场追迹法在近1秒内获得的。下面一行是将上面所获得的几何场追迹结果衍射传播到焦点所得到的结果
~@l NBF 模拟散射，光栅以及衍射透镜追迹通过两种介质间的界面的智能光线为几何近似的有效性提供了一个简单的标准。由于智能光线知道其临近光线的波前，界面后光线的相交很容易地就被自动探测到。相交光线表明波前被破坏了，因此几何场近似不在有效（Fig.2，中图）。很明显地，当光学表面不是光滑且仅小范围的进行调制则光线会相交，例如，带划痕的透镜或者一个散射粗糙表面。此种情况下，光学界面后的光场需要一种衍射方法以用于继续传播。此时，我们完全得益于统一场追迹概念，其能够让我们结合不同的电磁场建模技术以用于系统不同的区域[6]。
让我们仔细的看一下一个散射表面后的衍射传播技术。一般情况下我们可以使用傅里叶分析来将理想平面波分解以生成不同的传播积分来探讨这种情况。当然，每个平面波可以通过几何场追迹进行进一步传播，通过一个物理光学理由和策略将一条输入光线转换成大量的输出光线。如果我们假设光学界面至少是局部周期性的（Fig.2，右图），例如，一个光栅或者衍射透镜，依据光栅方程，那么仅在离散方向上才有平面波的傅里叶分析结果，如光栅级次。这就允许我们选择那些感兴趣的级次进行使用，如一个衍射透镜的第一级次。通过几何场追迹，就构建了一个光栅和衍射透镜的物理光学建模的简单说明：通过智能光线传播场到光栅或者衍射透镜，局部光栅效应使用电磁场分析，例如通过傅里叶模态法（FMM）[7]，然后由此产生的偏转的局部电磁场平面波通过智能光线追迹进一步传播。需要强调的是在模拟光栅和衍射透镜时，与传统光线光学的准确结合是无法实现的，而是要求智能光线提供必要的场信息以处理局部光栅效应。在VirtualLab Fusion中这项技术已经得到了应用，因此能够使用它来设计衍射透镜。然后，衍射透镜的结构可以以各种格式导出，如用于光刻的GDSII文件。Fig.5演示了一个用光栅将光耦合进和耦合出一个波板以用于色彩混合。
0J~4
Fig.5 通过一个光栅将RGB光源的光耦合进入一个波板，通过波板的全内反射进行传播，并使用三个光栅将光耦合出波板以获得一个部分重叠的混合色彩。本次模拟是通过FMM模拟光栅并使用几何场追迹的方法来完成的。上图显示了光线传播和几何场追迹结果
3
6
;hg # 干涉和相干建模 ,H*3_c&Q
Fig.6 Mach-Zehnder干涉仪对于一个时间部分相干光源的干涉图案。此建模是使用几何场追迹完成的，在PC上对超过100个波长进行计算，用时在1min以内
s? Kn,6Y
Fig.7 对于一个不同厚度的薄膜，通过多重干涉所获得的反射光的强度结果。使用光源为RGB白光。建模使用的是几何场追迹，用时少于1min
干涉和相干现象需要使用物理光学建模。因此在传统光线追迹中自然无法包含这些现象。由于几何场追迹是一种物理光学技术，因此包含了干涉和相干效应，并可对其进行研究。Fig.6显示的使用绿色LED光源的Mach-Zehnder干涉仪的模拟结果。在干涉仪的臂上分别有一个玻璃平板和一个透镜。对于完全地相干光，我们期望出现一个完美的对比度环形干涉图案。然而，LED光源有大约30nm的带宽，其会导致一个大约10um的时间相干长度。因此，干涉图案必定会有一个不同的对比度。正如干涉仪的几何场追迹所揭示的一样。在Fig.7中，显示的是一个RGB光源被不同厚度的油膜反射后的结果。这个建模是通过在薄膜结构中进行非序列几何场追迹完成的。
wG
O)!u 4 综述从麦克斯韦方程组出发，将传统光线追迹推广并引出一种物理光学建模技术，我们将其称之为几何场追迹。几何场追迹可以与包含了麦克斯韦方程组衍射求解器的衍射场追迹技术无缝的结合。几何场追迹建模包含了如偏振，干涉，散斑和相干等物理光学效应，且与传统光线追迹速度一样快。用于光学软件VirtualLab Fusion的所有技术同样用于[8]中例子。
# eYVZ=E 作者 kn s]P<g 1oPT8 )[U Frank Wyrowski是Jena的Friedrich Schiller大学的技术物理教授以及应用计量物理组的负责人。在1999年他合作创办了LightTrans公司，并于2014年联合创办了Wyrowski Photonics公司。他喜欢并致力于各种光学建模和设计任务，特别是物理光学领域。其研究成果源源不断地为光学设计软件VirtualLab Fusion的进一步开发做出贡献。 )Yml'?V" eAMT7 2_
,"o \_{<z il~,y8WTU{ Christian Hellmann是Wyrowski Photonics UG的CEO并且负责VirtualLab Fusion的开发工作。他代表Bosch股份有限公司在Horb的联合教育大学完成他的计算机科学学习后加入到了LightTrans股份有限公司。在LightTrans公司其负责光学建模的编程和算法开发。从2013至2014年，其是LightTrans首席产品官。
K@m^QioMj 参考文献[1] E. Hecht and A. Zajac. Optics.AddisonWesley, 1974
[2] J. D. Jackson. Classical electrodynamics.John Wiley & Sons, New York, 1975
[3] A.v. Pfeil, F. Wyrowski, A. Drauschke, andH. Aagedal. Analysis of optical elements with the local plane-interfaceapproximation. Appl. Opt., 39(19):3304{3313, 2000.article
[4] F. Wyrowski, H. Zhong, S. Zhang, and C.Hellmann. Approximate solution of maxwells equations by geometrical optics. InD.G. Smith, F. Wyrowski, and A. Erdmann, editors, Optical Systems Design2015:Computational Optics, volume 9360 of Proceedings of SPIE, page 963009,SPIE, Bellingham, WA, 2015, 2015
[5] M. Born and E. Wolf. Principles of optics.Cambridge University Press, 1999
[6] F. Wyrowski and M. Kuhn. Introduction to eldtracing. J. Mod. Optics, 58(5-6):449{466, March 2011
[7] J. Turunen. Diraction theory ofmicroreliefgratings. In H. P. Herzig, editor, Microoptics elements, systems andapplications, chapter 2, pages 31{52. Taylor & Francis, London, 1997
[8] Unied optical design software “Wyrowski VirtualLab Fusion”, developed by Wyrowski Photonics UG,distributed and supported by LightTrans GmbH, Jena, Germany
GSp1,E2J
<T ).+
M/
