简介：本文讨论了如何使用FRED对球透镜封装的半导体激光二极管耦合到单模光纤进行准确的建模，这是在光纤通信领域很常见的一个光学系统。该模型演示了FRED传播相干光场的能力、它的精确激光二极管束(Laser Diode Beam)光源模型以及准确的计算光纤耦合效率。 s%@HchZ 1 +5ue) ` 模型 ;s w3MRJ 4@"n7/< 在FRED模型中使用的半导体激光二极管是Mitsubishi(三菱) ML725C8F，这是一个InGaAsP / InP多量子阱(MQW)激光器，工作波长是1310nm。Mitsubishi光源说明书定义了输出光束的在x和y方向的发散角分别是25和30度（远场功率分布的全1/e宽度)。没有提及在x和y焦点位置的任何偏移，所以我们假定它们和光源处的分布是一致的。 B 0 ZLGB ~uw eBp~O 我们在FRED中使用激光二极管束光源类型对激光二极管光源建模，以及设置光源产生相干输出。 M<3m/l%`Y 图1. 激光二极管光源编辑
$m0-IyXcv Y@ 'ahxF 注意到在激光二极管光束光源的设置里面，发散角由功率的1/e2标准定义。这就要求制造商提供的发散角要乘以一个开方因子。 >5bd !b, 图2. 球透镜封装的激光二极管耦合到光纤系统原理图（侧视图）
k^Uk= )9 b-)3MR:4 直径为1.5mm的球透镜是Mitsubishi激光二极管集成的一部分，它的位置在距离激光二极管发射表面1.88mm处。 +KHk`2{y~ i&zJwUr(< 在FRED中使用球形元件基元，就可以创建该透镜。为方便起见，全局坐标原点选在球透镜的输出表面与光轴的交点处。 J85Kgd1
\a 图3. 全局坐标原点的定义
g?e-D.p SF RRJN@|" 值得注意的是，我们使用了FRED的N-BK7模型来定义球透镜的材料，在1310nm波长处折射率大小是1.5036。 ka [%p, H j WerX -$ 模型中使用的单模光纤(SMF)位于距离全局坐标原点1.9mm处，它的结构(由下图定义)基于单模光纤的典型值。光纤纤心的半径是5μm，且由直径为125μm包层包裹着。纤心和包层的折射率大小分别是1.465和1.47，它们之间的折射率差为0.36%。 \++#adN:K ul ag$ ge 图4. 单模光纤示意图
'sN
(=CQ zK ir 模型中还包含了一个吸收涂敷层，或者是夹层，覆盖在光纤表面。 (# VF>;;L : t6.J 在FRED中定义的光纤是一个组件，它包含了多个元件基元：一个圆柱体用于纤芯、光管用于包层和涂敷层。 ]?j[P=\ Avo"jN*<d 
~e R6 [; bPo*L~xdk 注意到“Fiber Cladding”管道的内壁恰好与“Fiber Core”圆柱体的外壁是重合的。为了正确的建模，用户需要手动的设置包层管道的内壁为不可追迹(Never Traceable)。不这样做的话将会导致光线追迹错误，因为两个表面放置在空间里完全一样的位置，而且它们具有两个不同的材料设置。对于“Fiber Coating”的内壁需要同样的设置。 tC'@yX w\Bx =a>vc 
,?> :Cdz4 *Q:EICDE7 在这一模型中光纤涂层认为是吸收的，且拥有停止所有(Halt All)光线追迹控制。所有其它的表面是不加涂层的。 GeCyq%dN x\!Uk! fM 仿真 .5YIf~!59 b}5hqI y FRED使用如下的方程来计算光纤耦合效率(CE)： H2D j`0 ^f bw0 
4XsKOv ZHW|P 其中Einc是入射场分布，Efiber是光纤基模的场分布(由FRED根据光纤规格参数自动计算)。 *<T,F yc| F/zbb 一般来说，CE是一个复数，所以耦合功率实际上是：CEpower = Re[CE]2 + Im[CE]2 <aEY=IF4 } l4d/I 因此，我们要想精确的计算光纤耦合，需要在光纤入口的后面放置一个分析面来保证该表面的反射系数能够准确的纳入考虑之中。 4.0JgX > aV
Q 非常重要的是，分析面是大于我们所期望的基模的模场直径(MFD)，以便进行精确的重叠积分。同样重要的是，我们应该意识到数值积分的精确性依赖于分析面中划分网格的数目。在本例中，50μm宽的分析面上251×251的网格，可认为是足够的。 K#oF=4_/| 图5. 分析面放置在光纤界面的后面
L?p,Sy<RI 图6. 光源用128*128采样点光线追迹与渲染
,T3_*:0hk! K h;jiK ! 由FRED光纤耦合效率计算得出的返回值是两个场分布之间的重叠部分，且没有考虑入射场的功率。因此要想知道多少功率耦合到该模式中一定要做到以下两步： 6=xbi{m$ }Qo:;&"3 1.通过辐射照度的计算确定分析面处的功率值(P) czHbdEh 2.通过光纤耦合效率分析确定CE的值 JYU0&nZl4 ^GN |}W 耦合到光纤模式中的功率大小可以简单的表示为P * CEpower。 YQk<1./}I ^9PB+mz 追迹完从具有2048×2048个样本点的光源发出的光线后，当我们计算辐射照度时，输出窗口里就会显示出到达光纤接口后面的分析面处的光源功率值。 2V 9vS 图7. 分析面处的积分功率值
H3Z "u _&
qM^ 可以看出，26.55%的光功率到达了分析面。为了确定到光纤模式中的耦合，这里使用了FRED光纤耦合效率分析。注意到0.005mm的光纤纤芯半径在这里需要准确的输入。 <xWBS/K 图.8 光纤耦合效率分析对话框
} I;5yk,o }v?_.MtS 点击完OK后，结果会显示在输出窗口中。 Sx%vJYH0 图9. 光纤耦合效率显示在输出窗口
auP6\kp Me zIjUfgO/M 可以看出，耦合效率为71.44%。因此，在这个系统总的耦合功率百分比为71.44%*26.55% = 19.0%。 235wl V`MV_zA2 ML725C8F激光二极管工作光源是在5mW，因此在该配置中，光纤传输的信号差点不到1mW。 #aV2+ `d 'WW[ ' 对齐灵敏度 @hE$x-TP0 y $K#M 对于测定设计公差以及激光二极管/光纤包的可行性，理解光纤对齐灵敏度是非常有必要的。使用FRED脚本功能可以很容易的完成这件事。 b#<@&0KE < BNCo5* 与该FRED文件相关联的共有三个内置脚本： | p1 pa4%} 纵向距离扫描 KoPhPH 横向偏移扫描 M)oJ06`K 倾斜扫描 !@<>S>uGG ^%q h E8 这三个脚本之间是相似的：通过用户控制的步长，每个脚本调整了光纤的位置、计算了耦合系数并打印到输出窗口或者到Microsoft Excel电子表格中(如果有需要)。 u LX V, e~, /Z\i 纵向对齐灵敏度 N7`<t&T@ j 4B|ktf 在距离扫描脚本文本的顶端，用户输入光纤的开始和结束位置，以及希望运行的扫描分辨率（步长）。 xe3t_y 0$9I.%4jAJ 如果用户希望FRED将数据打印到Microsoft Excel电子表格中并绘图，就要设置exportToExcel标签值为True。 /:j9 #kj 1li1& 
a ges-Z_X &E>zvRBQ 就在这定义了光纤的参数，这只是用于光纤耦合效率的计算。 xgeKz^, \dL# PI3 
+|w-1&- jJmg9&^R 头部打印出来后，脚本的主循环就开始了。这是一个“for”循环，它会一步一步的改变光纤的位置-[1]，追迹光线-[2]，计算照度并确定总功率-[3]，计算光纤耦合效率-[4]，最后计算模式功率-[5]。 1 JU1XQi 图10. 位置扫描脚本的主循环
){.J `X5r Y C
uuj$ 注意到函数FiberCoupleStepIndex返回了两个值-“coupleReal” 和“coupleImag”，这些变量是耦合系数的实部和虚部。 ?*~Pgh >uL m j{ /' 下图表示的是，对于球透镜到光纤的距离从1.5mm到2.5mm变化的结果。 n?QpVROo\ 图11. 光纤耦合vs距离
: Jh !:[kS1s>M 激光二极管的制造商Mitsubishi指定了在距球透镜1.9mm位置处，光纤耦合功率的最大值为0.8mW(16%的效率)，FRED在耦合中计算出了稍微偏大的值。这种差异可以解释为：耦合对光纤模式尺寸和折射率分布极为敏感。很遗憾的是，Mitsubishi没有给出使用光纤的具体细节。 I &e ,R +r+H`cT@ 横向准直灵敏度 Q\DD^Pbq Mp V6Vbp “横向偏移扫描”脚本与之前十分相似，除了用户为扫描定义了如下的参数： F+`DfI]/m 0-^wY8n-= 
图12. 在Z=1.86mm位置处：光纤耦合vs横向偏移
l|[8'*]r! OudD1( )W 方向灵敏度 cN> z`x l 7x%0 ^~/n 该脚本同样与先前的脚本十分相似，这里用户定义了取向的角度范围。注意到该脚本只是在水平方向倾斜了光纤，并不是一个任意的角度。 9 xvE?8;M# 图13. 在Z=1.86mm位置处：光纤耦合vs水平方向旋转
'M*+HY\.0 =~&Fq$$ 结束语 |xTf:@hgHf ` NC{+A 在本文中，FRED展现出了从激光二极管到光纤耦合准确计算的能力。其计算结果与激光二极管生产商提供的耦合信息一致。FRED的相干传输能力以及高散射相干的精确定义对于这种类型问题的仿真是很关键的。 'CDRb3w}B }@pe `AF^ 本例系统数据(单位是mm) G B+U>nf
