如何在Zemax下模拟单模光纤的光束耦合

1、如何在Zemax下模拟单模光纤的光束耦合设计前的准备Zemax公司感谢SussMicroOpticsSA公司的ReinhardVoelkel博士提供本文使用到的实验数据。我们同时提供本文的的日文版本本文描述了一种商用的光纤耦合器，系统使用SUSSMicroOpticsFC-Q-250微透镜阵列来耦合两根康宁(Corning)SMF-28e光纤。如下图所示：0.9mm2min0.9mm供应商提供的上述元件的参数如下：单模光纤，康宁SMF-28e数值孔径0.14纤芯直径8.3pm模场直径1.31pm9.20.4pm微透镜阵列，SUSSMicroOpticsSMO39920基片材料熔融石英基片厚度0

2、.9mm内部透过率0.99透镜直径240pm透镜节距250pm曲率半径330pm圆锥常数(Conicconstant)0数值孔径0.17附件中的文件singlemodecoupler.zmx是整个系统的Zemax文件。请注意一下几点：物面到透镜的距离和透镜到像面的距离设定为0.1mm,是因为这比较接近实际情况。后面经过优化过程时候，这个尺寸还会发生变化;透镜到像面的距离使用了Pick-upsolve，以确保和前面的物面到透镜的距离之间相等。既然两组透镜和光纤之间是完全一致的(在制造公差之内)，因而整个系统也就应该是空间反演对称和轴对称的(eitherwayround)；两个透镜之间的距离设定为

3、2mm，因为这个是实验中使用的数据。同样地，这个距离后面也将会被严格的优化；系统孔径光阑设定为根据光阑尺寸浮动(floatbystopsize)，而光阑设定在第一个透镜的后表面。这就意味着系统的孔径光阑由透镜的实际孔径决定。因而光纤的模式在这个系统中传输的过程中，就有可能受限于透镜的实际孔径。在这个例子中，光纤的模式要比透镜的实际孔径小很多。当心“数值孔径”的多种不同定义。它有可能指的是边缘光束倾角的正弦值，有可能是光强降低到1/e2时的光束倾角的正弦值(我们将会看到Zemax会在不同的场合使用这两种定义)，也有可能定义为光强降到1%峰值强度时光束倾角的正弦值，康宁便使用这种定义。这些非常重要

4、！孔径上定义了高斯切趾(Gaussianapodization)，用来产生光束的高斯分布。当前这只是一种近似，后面将会做进一步的精确的计算。透镜孔径的大部分区域是衍射受限的光学质量的，并且被光纤模式照射到的区域是衍射受限的。Mj-ILawutUpdateSettingsPrintWindowZoom使用旁轴高斯光束计算旁轴高斯光束算法是最简单可以用来分析光纤耦合的分析方法。不过，这种方法只能获得对系统性能初步的了解。根据康宁的产品参数表，光纤在1.31pm波长下的模场直径为9.20.4pm。因此，我们按照下图所示的情形设置旁轴高斯光束计算(AnalysisPhysicalOpticsParax

5、ialGaussianBeam):图中光束的束腰直径Waist总是相对于表面1来计算的，在本例中它和物面出于同一个位置。因此，高斯光束的束腰直径4.6pm就位于源光纤的位置。光束然后传输经过光学系统。从上图我们可以看出表面3上的1/e2光束直径是65pm,而表面4是70pm。这些表面的实际的物理半口径为120pm。也就是说大约两个光束直径以外的光将会被阻隔掉。另外需要注意的是像面并非位于光束最佳聚焦聚焦点处：像面处光斑的大小为5.3pm，而其实根据系统的对称性的假定，高斯束腰直径应该是4.6pm。我们将会优化表面1的厚度（同时也会通过Pick-upsolve来控制表面5的厚度）来改进这些。请注

6、意表面5的厚度是通过Pick-upsolve来控制的，因为我们希望系统倒过来使用时也能得到同样的耦合效果：我们使用了完全相同的两组光纤和透镜（在制造公差范围内），因而我们期望最好的系统是输入输出互易的。Zemax中有一项优化操作数（operand）GBPS，指的是旁轴高斯光束尺寸，可以用来优化光纤和耦合镜头之间的距离。根据系统的上述对称性，我们知道，高斯光束束腰的最佳尺寸是4.6pm，因此，优化函数就只有简单的一行，如下图：优化后给出的光纤到透镜之间的距离为0.117mm，下面是相应的旁轴高斯光束的数据:上述便是旁轴高斯光束计算所能给我们提供的信息。相关的Zemax文件为optimizatio

7、n.zmx。返回目录使用单模光纤耦合计算单模光纤耦合计算方法(位于AnalysisCalculationsFiberCouplingEfficiency)提供了更加有力的用来分析具有高斯分布的光纤模式的能力。它会执行两种计算：能量传输计算(energytransportcalculation)和模式匹配计算(modematchingcalculation)。系统效率(systemefficiency：S)是用通过入瞳(entrancepupil)的所有光能量，并且考虑了渐晕(vignetting)和偏振光传输情况下，经过系统之后的能量的总和除以从光纤辐射的所有能量得到的：这里Fs(x,y)指的

8、是源光纤的振幅函数，分子是仅在光学系统的入瞳处的积分，而t(x,y)是光学系统的振幅传输函数。传输过程受到体吸收和光学镀膜(打开偏振传输的情况下)的影响。光学系统中的像差所引起的位相差也会影响光纤的耦合效率。当向接收光纤传输的汇聚波前的各个点上的模式完全和光纤的模式(包括振幅和位相)想匹配的时候，耦合效率达到最大。它的数学描述是通过光纤和波前振幅之间的重叠积分(overlapintegral)来定义的：这里Fr(x,y)用来描述接收光纤复振幅函数，W(x,y)是光学系统出瞳处的波前的复振幅函数，而符号代表了复数共轭操作。注意这些函数都是复数形式的，因而这个表达式是相干重叠积分。T的最大值为1.

9、0，并且随着光纤的振幅和位相和波前振幅位相之间的失配的增加而降低。Zemax会计算上述的S和T的值。总的功率耦合效率也是从这些数据得来的。Zemax也会计算理论上的最大耦合效率，这个计算是在胡略了像差但是考虑模式的渐晕、传输和其他振幅失配因素之后得来的。在计算中，源光纤模式和接收光纤模式是通过高斯光束的数值孔径NA(定义为物方或者像方的介质折射率和光束上功率降到1/e2处的半张角的正弦的乘积)这个角度可以通过下面两种方法计算获得：通过高斯光束计算得来的发散角，使用模场直径来定义光束束腰直径(参见本文前面的计算)；根据康宁产品参数表中提供的1%功率处的NA计算得到1/e2处的数值孔径。通常，NA

10、大约为0.09,因此耦合效率的计算设置如下图所示：得到的计算结果如下:我们也可以一行优化函数来优化系统的耦合效率，该操作数为FICL：OptimizationOperagd1经过几个优化循环之后，光纤到透镜的距离变成了0.11mm(而在旁轴高斯光束计算中为0.117)，详细的结果如下:注意一下几点：系统耦合效率(systemefficiency)并未显著改变，因为这主要是由表面的孔径和模式尺寸决定的，而轻微的离焦对其的影响并不大；接收端耦合效率(Receiverefficiency)得到了提升，因为重新聚焦使得源光纤模式在经过系统传输之后更好地和接收光纤的模式匹配；最大耦合效率(maximum

11、efficiency)可以通过下述方法得到提升：增加非球面，增加额外的表面等等。本例中，这个耦合效率基本达到了极限。本节使用的zemax源文件为afterFICLoptimization.zmx返回目录使用物理光学计算(Usingthephsicalopticscalculation)在单模光纤耦合计算的基础上运用物理光学传输算法可以极大的扩展这个方法的应用。类似地，我们需要计算基于物理光学传输的重叠积分，这样做有一下几点好处：可以定义任何复杂模式，而不仅限于高斯模式；当接收光纤位置的时候，光纤耦合重叠积分可以在任何一个表面上计算，这包括当不仅限于代表光纤的表面；外部程序，比如光束传输(Bea

12、mPropagation和时域有限差分方法(Finite-Difference-TimeDomain)代码，可以使用成.zbf文件格式或者生成DLL接口，用来计算光纤(或者任何集成光学元件)的模式结构，并且可以将其表达成和室的复杂振幅分布函数的形式。关于这一点，可以参看这篇文章；由于传输受光阑限制或者远距离传输产生的衍射效应也可以得到准确的模拟。POP计算可以通过以下窗口设置，点击AnalysisPhysicalOpticsPhusicalOpticsPropagation打开设置窗口:|POPSetting?区POPGeneralBeamD&finitioriDi护layFiberDataT

13、iltAboutY(deg):o(/ComputeFiberCouplingIntegralTiltAboutX(degJ:cTOK.CancelFih已Typ已：|GaussianWaist二|、File:1Waistk0.0046DecenterX|0.0046DecenterY|aJApertureXApertureYOrderXOrder丫FiberPosition:|ChiefRay3以上设置了表面1上的束腰半口径为4.6pm的高斯模式，并且经过系统传输至像面，在像面上我们使用完全相同的模式来计算重叠积分。在光束的参数设置表中，我们词啊用了256x256的采样率，然后需要点击Auto

14、按钮来设定合适的起始矩阵；物理光学窗口中关于光纤耦合的报告，可以参看下图中红色框选区域。优化函数编辑器中的POPD操作数可以给出所有的物理光学数据，这通常是更有用一些。具体的关于POPD操作数的叙述可以参考用户手册中的优化那一章。POPD操作数直接调用先前保存的POPAnalysis窗口中的设置，因此，在优化之前先要去设置一下这些参数。412PhysicalOptjfsPopagatjpQUpdateSettingsPrintWindowTextZoom13.0QQ03,-0.0167目.QDOSYCOQRDTNFlTEVl=lLIJEQ.0167PHFI3EY-CR口三爲SECTI叮卜SUR

15、FACE七SMF2E：Effl.FRDDND即煦33.P2-Q-2EH.LM昭.ROD33M_MaFZHS12QLFLIM3SfflMfl._.B_ENGTH.3SEII炬INIMZEJtL.EIZIZI93HFEI.HQQODEGE闌陋讪心虫唧临册臥山用乩卩皿此I山HPLHf亠尸正民EFFICIENCi-iSY37EHEh冇門7乩FBZEI/ER0r?557；3aCDUFIh7J站召筈韶冷id牡Hg铝亍申64詔POPD数据（通过优化函数得到）描述Description数值Value26M21.0820TotalFiberCoupling0.9941SystemEfficiency0.9982

16、ReceiverEfficiency0.99610PilotBeamWaist4.57pm23EffectiveBeamWidth4.84pm位相是你需要关注的最有用的信息，因为IrradianceProfile几乎是理想的高斯形（M2为1.082）。接收端的模式实际上在每个地方都几乎为零，因此位相也直接告诉我们模式失配的程度。要显示位相信息，打开POP设置窗口的显示（Display）那一栏，按照下图所示的方式设置：注意位相图的抛物线和四次曲线的形状，这等价于聚焦和球差。另外也要注意透镜边缘对位相曲线产生的影响。根据系统的效率（Systemefficiency），我们知道由于透镜外相尺寸的限制

17、，系统约有小于1%的能量损耗。如果我们进一步对光纤耦合效率进行优化（记住光纤和透镜的距离是唯一的变量），我们可以得到少量的提升:412:PhysicalOptjfsPopagatjpQIWmsMOE6出至dUpdateSettinaePrintWindowTextZoom0.016b-O.066PHFI5EY-CR口三巧SECTI叮卜5URFHCE65MF3E：Effl.PFDDND阳煦33.K-D-2E0.LK昭.ROC33H.I1.FZBBinUFLIMJSfflMft-.aENCTM3SEIEI炬INIMZEIL.EIO3EI3HFELIZQQODEG=0-JTtRa：（13BSEiEE

18、tQQQ:B-ITEBPHFtE=-D.ElSIRRDIRJS4FH=tEREFTliflPLRtEFIES!EFFICIENCYiSTO7EHEi.ViSLL.KBZEIVER0C吒兀1LCDUFLZJj0.?WEEVPILai：IZE=H.NE-QBS,URZST=耳弓&?9爾1PDE=；h韭E圧-3日盒FhVLE1G：H=5QEE=2E-Q32光纤耦合效率并未得到显著提升，这是因为位相差产生的位置对应的能量非常低:上图可以通过下述方法得到：从POP窗口中输出cross-section形式的irradiance图，然后克隆它(windowsclone)，然后将克隆后的窗口设置成显示位相，而

19、非光强，介质在windows菜单中选择overlay.将这两个曲线重叠成上图的样子。上述内容涉及到的Zemax源文件为afterPOP.zmx到现在位置，我们讨论的旁轴高丝光束，单模光纤耦合和完整的POP算法之间只有一些比较微妙的差异，他们的计算过都比较相似,要得到更多的关于他们的本质属性的结果，我们就需要更加复杂的分析。不过，随着耦合器长度的增加，衍射效应也变得越来越重要,这几种分析方法的差异也会更加明显，而POP方法的优势也会逐渐显现出来。我们如果将透镜到透镜的距离设置为20mm，基于光线的单模光纤耦合算法得到的结果几乎不会发生变化：FICL计算出来的耦合效率为0.99。这是因为光束在两组

20、透镜之间是近似准直传输的。然而，POP算法得到了0.57的耦合效率：几乎减半。这是由高斯模式在两个透镜之间的空间中的衍射和尺寸的变化引起的。传输了20mm之后，高斯模式的1/e2宽度增加到了0.15mm，这已经达到了透镜的0.12mm的尺度。结果，位数不少的一部分光能量便在第二个透镜的物理光阑出产生衍射，下图显示了经过第二个透镜之前和之后的光强分布和位相分布的重叠图。明显聚焦到接收光纤上的光束已经明显不是高斯形式，这时它的M2为2.45.412:PhysicalOptjfsPopagatjgQUpdateSettingsPrintWindowTextZaom3.10E-tO0H231UJ百Tl

21、lz=w-tSdA匚玉L.55E+O0H.EI0E+CI0DSC020QQ0000YC0DRDIMA7EURLUE工良FI匸|工RMCE-CRUSSW；E匚TIONSLIRFH匚E7S-F2：3曲0,ffiUJMDFC-Q-S3.DOIELRDC23UM.FUfBHUPLnNL昂HklFfuaENGTHLi刖I3EWgINIMDEC1.CBEiEiDHTQ.EiDEQOEGDB-ITER.X=CL胛日Ei曰CMFEFKIRHRDIR-GE-3.LECTEtEiEiU廻门旨用11_口朮吒甘N1UTBLFHUEB-已丽MEEC1UflITSFIES!EFFIQHJEt!EYTEnD.7S395dE

22、ELECVERQ.7025*?匚LFLHNGD.DSTr尸JLRISIZE-CiitlEKE-S,trasi-ras-3.ldE-BPS.RrllLEIH-|-1i日治相-阳2这和用光线计算出来的结果是不一样的。将2D外形图的Y方向放大，便可以得到如下图所示的光线分布WLIjLawut实际上，由于球差的存在，光线计算的结果是在第二个透镜上的光束口径要小于在第一个透镜上的光束口径。对于这种近似准直光束远距离传输的问题，POP的算法要比光线光学的方法准确。并且，POP可以队耦合器进行严格的优化。将我们已经优化过的光纤到透镜的距离设置为固定，然后将20mm的透镜间隔设置为变量，经过几个优化循环之后，

23、便会得到2.072mm的最优化透镜间隔。在实验中采用的2mm数值，其结果和它之间的差别是微乎其微的。使用UniversalPlot，我们可以看到透镜间隔对光纤耦合效率的影响。:UniversalPlot1UpdateSettingsPrintWindowTextZoom!-JKEKn-D也LHMUJU10耳曰0.9t.0.S2Q.?0.SE!Q.st.0.Hi0.82THICKNESSINMILLIMETERSON5URFHCE3FIBERCULIPLIIG-.-ERSUSLEI-SSEFHRHTiaNSKI.PGLONO320.PJ-0-2EB.LXHB.0CC训儿FIBEBDJUPLIUJMlRRB?ENFffiLE阳LPT3NHRTJEHZNJ(CONFICURRTION1DF1类似地，当光束在两个光纤之间传输的过程中