黄永超，张廷蓉

(1.内江师范学院工程技术学院，内江641112;2.四川师范大学物理与电子工程学院，成都610066)

引 言

非线性光学不仅从理论上丰富了人们对光与物质相互作用的认识，而且已经得到广泛的应用。因此，研究非均匀介质在激光束作用下产生的非线性现象及其应用是一个非常热门的课题。梯度折射率介质是比较典型的非均匀介质，它在光纤耦合器、光通信、光学设计应用开发等方面具有广泛的应用前景。目前人们对梯度折射率介质的特性及其应用做了深入的研究［1－5］。梯度折射率介质中余弦高斯光束、平顶高斯光束、洛伦兹高斯光束、双曲正弦高斯光束和空心高斯光束的传输特性也有广泛研究［6－10］。双曲正弦高斯光束(sinh－Gaussian beams，ShGB)是厄米余弦高斯光束的特例［11］。对双曲正弦高斯光束中双曲正弦函数的变量取复数，称之为复变量双曲正弦高斯光束(elegant sinh－Gaussian beams，EShGB)［12］，双曲正弦高斯光束在直角坐标系下各轴上光强为0，但是各象限中双曲正弦高斯光束的光强分布具有明显的对称性，并且双曲正弦高斯光束可看作是具有相同束宽的两偏心高斯光束合成。基于此，研究复变量双曲正弦高斯光束具有重要意义。作者将运用广义惠更斯－菲涅耳衍射积分法推导出EShGB在折射率径向分布介质中的传输光场，运用空间二阶矩的定义解出光斑尺寸的的表达式，并数值模拟介质中EShGB传输特性，为其在实际应用中提供一定的理论参考。

1 EShGB通过梯度折射率介质的传输公式

在z=0mm 的平面上，EShGB 的场强分布［11－12］为:

式中，w0表示EShGB的束腰宽，α和γ为与双曲正弦函数项相关的参量，x0是入射面上的横向坐标，E0是入射面的场分布。EShGB通过梯度折射率介质中从源平面到观察平面z之间的光学系统用传输矩阵表示为［1］:

式中，A，B，C，D为1阶光学系统的变换矩阵元，β为梯度折射率系数。当EShGB通过这一光学系统的传输时，可由广义惠更斯－菲涅耳衍射积分描述［13］:

式中，k=2π/λ为波数，λ为波长。选取梯度径向分布的介质作为研究对象，其梯度折射率径向分布介质的折射率可表示为［7］:

式中，n0表示介质轴上的折射率。将(1)式带入(3)式，得EShGB通过梯度折射率介质的传输场为:

式中，x'=x/w0为相对坐标，a=αw0和b=γw0为光束参量，Z0=πw02/λ为瑞利长度。将(2)式带入(5)式，由I(x'，z)=E(x'，z)×E*(x'，z)，可得 EShGB 在梯度折射率径向分布的介质中任意一点的光强分布。其中E*(x'，z)表示E(x'，z)的共轭。

分析(11)式可知，当a=0，b=0时，光斑尺寸在z=0mm平面上取最小值(即束腰宽度)。对(11)式进行微分，可以得到EShGB在梯度折射率介质中传输时，在观察面z上的EShGB的光斑尺寸的变化率:

式中，M=(1+b2+F)/(βz)2－S，S=1+a2+F，H=ab/(Z0β)是为简化引入的变量。

2 数值计算及分析

运用(1)式对EShGB在z=0mm平面上的光强分布作了数值模拟(本文中计算都取λ=632nm，w0=0.96mm)，由图1可以得出:当x=0时，光强为0，这是由EShGB性质所决定的;当a=0.2，b=15(a较小，b较大)时，光强分布具有正弦高斯光束的形状;当a=2，b=0.2(a较大，b较小)时，光强分布表现出来了双曲正弦高斯光束特征，这是由于b较小时，EShGB趋于退化为双曲正弦高斯光束;当a=2，b=15(a，b较大)时，光强分布既不具有正弦高斯光束性质，也不具有双曲正弦高斯光束特征。

Fig.1 Intensity distributions at the plane z=0mm

运用(6)式做数值模拟，取折射率参量 β=0.1mm－1，图2中给出了光束在观察面(z=15mm)上的光强分布，由图2可知，当EShGB在梯度折射率介质中传输时，光强分布保持了z=0mm平面上的特征，即梯度折射率介质不改变光束的形状。对比图1和图2可以知道，折射率的空间分布对光强峰值及光斑有影响。

Fig.2 Intensity distribution at the plane z=15mm

为了进一步说明折射率的空间分布对光束光斑的影响，利用(11)式进行数值模拟，取折射率参量β=0.1mm－1，结果如图3所示。当EShGB在梯度折射率介质中传输时，随着传输距离的增加，光斑尺寸出现周期性变化;随光束参量a增大，光斑尺寸振荡幅度增大;随偏心参量b增大，光斑尺寸振荡幅度减小。由此可知，光斑尺寸振荡的幅度由光束参量a和b决定。

Fig.3 Spot size versus propagation distance z

为了进一步说明折射率系数β对光斑的影响，图4中给出了当折射率系数β取不同值时，光斑尺寸随传输距离z变化，由图4可以知道，光束参量一定(a=2，b=2)，EShGB在梯度折射率介质中传输时，随折射率系数β的增大，光斑尺寸的振荡周期性发生了变化;并发现，折射率系数β越大，振荡周期越小，即随折射率β的增大，在传输方向上周期被压缩，但是光斑尺寸没有变化。这表明梯度折射率系数决定了光斑尺寸变化的周期，但是对光斑尺寸振荡的幅度没有影响。

Fig.4 Spot size versus propagation distance z for different gradient－index parameters β

运行(12)式做数值模拟，取折射率参量 β=0.1mm－1，结果如图5所示。当EShGB在梯度折射率介质中传输时，随传输距离的增加，光斑尺寸变化率在前半个周期内为负值，由0逐渐减小到最小，后半个周期突然跃迁到正的最大值再逐步减小到0。这说明前半个周期光斑尺寸逐渐减小，即光束在聚焦;后半个周期光斑尺寸在逐渐增大，即光束在发散。由图5还可以看出，在同一观察面上，当光束参量a的增大，光斑尺寸变化率的范围在扩展;当光束参量b的增大，光斑尺寸变化率的范围被压缩。

Fig.5 Change rate of spot size versus propagation distance z

为了进一步说明折射率系数β对光斑尺寸变化率的影响，图6中给出了当折射率系数β取不同值时，光斑尺寸变化率随传输距离z的变化。由图6可知，光束参量一定(a=0.2，b=0.2)，EShGB 在梯度折射率介质中传输时，随折射率系数β增加，光斑尺寸变化率的范围在扩展;同时发现光斑尺寸变化率的周期被压缩，表明折射率系数增大，光束扩展或聚焦都将变快。

Fig.6 Change rate of spot size versus propagation distance for different gradient－index parameters β

3 结论

研究了梯度折射率介质对EShGB传输的影响，结果表明:当复变量双曲正弦高斯光束在梯度折射率介质中传输时，光束出现了周期性变化，周期由折射率系数决定，并且光斑尺寸随传输距离的增加，光斑尺寸出现了余弦平方的变化规律，光斑尺寸振荡幅度由光束参量确定，当随光束参量a的增大，光斑尺寸振荡幅度增大;随偏心参量b的增大，光斑尺寸振荡幅度减小。因此，通过传输控制参量可以得到余弦平方变化的光信号，从而得到数字脉冲光信号;同时发现，光斑尺寸的变化率随传输距离变化也具有周期性，通过对比可知，光斑尺寸变化率跃迁的位置就是光斑尺寸最小的位置，在同一观察面上，当光束参量a的增大，光斑尺寸变化率的范围在扩展;当光束参量b的增大，光斑尺寸变化率的范围在被压缩。在实验中，EShGB可由厄米－高斯模相干合成而产生，主要用于表征激光器的远场分布，因此，对于描述大功率半导体激光器的光场有着重要的意义。这些结论有助于大功率半导体激光器等方面的开发应用。

［1］ SICGMAN A E.Lasers［M］.Los Angeles，California，USA:University Science Books，1986:586－589.

［2］ BALLAV M，CHOWDHURY－CHAOS A R.A generalized nonlinear Schrödinger equation and optical soliton in a gradient index cylindrical media［J］.Solitons ＆ Fractals，2007，31(4):794－803.

［3］ JOHN H，LORENZ S D.Web inspection using gradient－indexed optics［J］.Proceedings of the IEEE，2005，41(6):1476－1482.

［4］ BALLAV M，CHOWDHURY A R.A generalized nonlinear schodinger equation and optical soliton in a gradient index cylindrical media［J］.Chaos Solitons and Fractals，2007，3l(4):794－803.

［5］ YABLON A D，BISE R T.Low－loss high－strength microstructured fiber fusion splices using GRIN fiber lenses［J］.Proceedings of the IEEE，2005，17(1):118－120.

［6］ SONG H Y，ZHANG T R，CHENG S H，et al.Propagation properties of flattened Gaussian beams in gradient－index medium［J］.High Power Laser and Particle Beams，2011，23(10):2630－2633(in Chinese).

［7］ SONG H Y，ZHANG T R，CHENG S H，et al.Propagation properties of cosine－Gaussian beams in gradient－index medium［J］.High Power Laser and Particle Beams，2011，23(4):890－894(in Chinese).

［8］ ZHOU G Q.Propagation of Lorentz－Gaussian beams in gradient－index media［J］.High Power Laser and Particle Beams，2013，25(1):42－46(in Chinese).

［9］ LIU L，HAO Z Q.Propagation of sinh－Gaussian beams in gradient－index medium ［J］.Laser Technology，2013，37(1):126－129(in Chinese).

［10］ HUANG Y C，LI C J，ZHANG X L.Propagation properties of hollow Gaussian beams in Gradient－index media［J］.Laser ＆ Optoelectronics Progress，2014，51(3):032601(in Chinese).

［11］ CASPERSON L W，TOVAR A A.Hermite sinusoidal Gaussian beams in complex optical systems［J］.Journal of the Optical Society of America，1998，A15(4):954－963.

［12］ WANG L，KONG R X.Propagation properties of the elegant sinh－Gaussian beams［J］.Opto－Electronic Engineering，2006，33(2):45－49(in Chinese).

［13］ COLLINS S A.Lens－system diffraction integral written in terms of matrix optics［J］.Journal of the Optical Society of America，1970，60(9):1168－1177.

［14］ ERDELYI A，MAGNUS W，OBERHETTINGER F，et al.Tables of integral transforms［M］.New York，USA:McGraw Hill，1954:856－859.