垂直腔面发射激光器的温度特性研究

戚向涛，顾亚平*，张　曼，方斯喆

(1.中国科学院大学，北京 100049；2.中国科学院 声学研究所，北京 100190)

引　言

激光信号具有远高于电信号传输速率的特点，更适合于未来高速率的传输网络[1-3]。激光器作为光纤通信系统的核心器件是系统仿真中考虑的一个重要因素。垂直腔面发射激光器(vertical cavity surface emitting laser,VCSEL)具有使用简单、功耗较低等特点，目前广泛应用于光纤通信和光信号处理领域[4-5]。VCSEL激光器输出的光功率与器件的温度相关，当器件温度改变后，激光器输出的光功率也会相应发生变化[6-8]。准确地分析激光器的温度特性可以保证激光器在合理的外界环境温度范围内使用[9]。

国内外很多学者研究仿真温度对VCSEL性能的影响，并建立了大量的模型。NAKWASKI建立了2维热效应模型[10-11]。YU等人建立用于单模模拟的VCSEL的速率方程模型[12]。MOROZOV等人考虑多模行为的速率方程模型[13]。SU等人建立的静态热特性模拟的温度相关模型[14]。

本文中通过建立光功率与工作电流关系(P-I)模型，并简化模型参量及引入电压电流关系(U-I)特性曲线来优化分析模型。使用Levenberg-Marquardt(LM)算法来实现模型的求解。对比分析20℃下的测量数据与拟合数据的相似度，进一步来预测不同温度下的P-I特性曲线，并分析用户可正常使用网络的温度范围。

1　VCSEL的P-I模型及优化

1.1　P-I模型建立

VCSEL的输出光功率P0和外部驱动电流I满足可以用下式来表示,这里P0与I的建模关系称为P-I建模[15]:

P0=η(T)[I-Ith(N,T)]

(1)

式中,P0表示激光器输出的光功率；I表示注入到激光器的外部驱动电流，包含外部加载的偏置电流Ib和信号电流，在没有信号输入时驱动电流为偏置电流Ib；η(T)是P-I曲线的斜率，其大小与温度T相关，从能量转换角度看，斜率相当于转换效率：斜率越高，相同电流I对应的输出光功率越高，那么相同电能可转换的光能越多;Ith(N,T)表示阈值电流，其大小主要和载流子数N及温度T相关。

为简化(1)式的表达，可以假设转换效率η(T)受温度影响较小[16]，即η(T)近似于常数η。同时假设忽略空间烧孔效应的影响。阈值电流可以表示为只与温度相关的函数[14]，如下式所示：

Ith(N,T)=Ith,0+Ioff(T)

(2)

式中，Ith,0为常数，Ioff(T)是与温度相关的经验热偏置电流，该电流随激光器的变化而变化。Ioff(T)可以由温度的多项式表示：

(3)

式中，an是多项式的参量，n表示阶数。

此时(1)式可以表示为：

(4)

式中，η(T)表示转换效率，是P-I曲线的斜率。本文中假设η(T)受温度影响较小，η(T)近似于常数η，模型求解中，η同样为待求参量。温度T受环境温度T0以及自身的温度影响，自身的温度与器件产生的瞬时功率Pi相关，瞬时功率Pi=UI，即受U-I(电压-电流)特性影响：

(5)

式中，Rth表示VCSEL热阻抗；τth表示热时间常数；T0代表环境温度；此处I为偏置电流Ib，U表示输入电压。

假设温度不发生变化时，(5)式后面微分项可忽略不计。结合(2)式～(5)式可以得到P0和I的隐式关系：

P0=η{I-Ith，0-

(6)

1.2　P-I模型参量的优化

模型的参量多少决定该模型整体的运行速度。参量太多会使得模型计算量增大，并且容易产生过拟。所以适当地优化参量个数在优化网络的同时会提高模型的计算精度。

将(4)式中的高次多项式展开可得：

P0=η(I-Ith，0-a0-

a1T-a2T2-a3T3-a4T4-…)

(7)

根据(7)式可以看出，参量Ith,0(Ith,0为常数，模型求解中Ith,0为待求参量)和参量a0皆为1阶常量，这里可以将这两个参量合并为一个参量,即可以将(3)式优化为：

(8)

故而可以在模型的整体计算上减少一个参量。随着n的增大，的值越来越小，并趋于0。这里假设n>4时，趋近为0,所以取n=1,2,3,4。

此时，P0和I的隐式关系为：

P0=η(I-Ith，0-a1T-a2T2-

a3T3-a4T4)

(9)

式中，T=T0+(UI-P0)Rth。

1.3　VCSEL的U-I特性曲线

温度T受外界环境温度T0和自身的温度影响，自身的温度与器件产生的瞬时功率Pi相关，瞬时功率Pi=UI相关，即受电压-电流(U-I)特性影响：

T=T0+(UI-P0)Rth

(10)

电压的变化关系和VCSEL的特性相关。为了简化模型，可以认为电压主要受温度和电流的影响，即：

U=f(I,T0)

(11)

在实际应用中，电压会随着温度和电流的变化而变化，所以可以用多项式来拟合电压U的输出[6]，如下式所示:

(12)

式中，bi和cj为多项式的拟合参量；n，m为多项式的拟合的阶数，0≤i≤n，0≤j≤m。利用(12)式，根据在20℃的环境温度下得到的实验数据U及I，拟合其U-I特性曲线。在20℃下，当n=2，m=4时,U-I特性曲线与实测数据的对比图如图1所示。

Fig.1a—U-Icurve of measurement data and fit data at 20℃b—U-Icurve of different temperature

图1中的横坐标代表驱动电流，纵坐标代表电压。图1a中的实线是温度为20℃时的实测U-I数据，虚线是经过(12)式且n=2，m=4时拟合后的曲线。从图1a可以看出,拟合的数据可以较好地反映实际数据的分布。图1b为根据拟合参量求得的在10℃，20℃，30℃，…，90℃温度下的U-I特性曲线。通过图1b可以看出：(1)在固定温度下随着驱动电流的增加，电压也会随之增加；(2)在固定的驱动电流下，随着温度的增加，电压也会随之增加。

此时，模型整体修正为式：

P0=η(I-Ith,0-a1T-a2T2-

a3T3-a4T4-…)

(13)

式中，

(14)

2　基于LM算法的模型求解

LM算法结合了Gaussian-Newton法和梯度下降法，利用标准值优化初值，是一种快速算法。具有稳定性强和快速收敛等优点，是快速有效求解目标函数的常用方法[17]。所以选择LM算法来求解(13)式。

在LM算法中，存在如下关系式[18]：

uk+1=uk+Δu

(15)

式中，uk为第k次迭代时的输入向量，uk+1为第k+1次迭代时的输入向量，Δu为两次迭代时向量的变化量。此处u=[η,Ith,0,Rth,a1,a2,a3,a4]是7维向量，是所需要求解的参量。

对于牛顿法有：

Δu=-[▽2E(u)]-1▽E(u)

(16)

(17)

式中，E(u)是误差指数函数，▽E(u)为梯度，▽2E(u)为Hessian矩阵，e(u)为期望功率和实际功率误差误差向量，P0为期望输出功率值，P为实际输出功率值。在此处，P是已知的，由(14)式得到关于未知数η,Ith,0,Rth,a1,a2,a3,a4的隐函数关系式。

对于Hessian矩阵：

▽2E(u)=[A(u)]Te(u)+S(u)

(18)

(19)

式中，A(u)为Jacobian矩阵。当目标接近极值时，S(u)很小，可忽略不计，此时：

▽2E(u)≈[A(u)]Te(u)

(20)

对Gaussian-Newton法：

▽2E(u)≈[A(u)]TA(u)

(21)

迭代过程中，对E(u)的Hessian矩阵无法求逆时，需对E(u)的Hessian矩阵进行变形：

▽2E(u)≈[A(u)]TA(u)+μI(u)

(22)

式中，I为单位矩阵；μ为阻尼系数，且μ>0为常数。故有：

▽u=-(AkTAk+μI)-1AkTek

(23)

uk+1=uk-(AkTAk+μI)-1AkTek

(24)

上面(15)式～(24)式即为LM算法的表达式，ek表示e(uk)。每次迭代时重新计算误差函数E(uk+1)，若其小于E(uk)，则保留此次迭代结果，减小阻尼系数μ来提高计算精度，并继续进行迭代，令μ=μ/β(β>1，β为大于1的常数)；反之则增大阻尼系数μ后重新计算，令μ=μ·β，直到满足E(uk+1)

3　实验结果及分析

在20℃的环境温度下得到的输出光功率P0和外部驱动电流I的实验数据，根据改进的P-I模型，使用LM算法来求解模型，拟合P-I模型曲线。并给出不同温度情况下(10℃，20℃，30℃，…，90℃)P-I曲线的预测结果。

模型的初始值和最终模型求解出的参量值如表1所示。

Table 1　The initial value and the solved value of the improved P-I model

3.1　20℃时P-I曲线

图2是温度为20℃时的实测数据与拟合数据的P-I曲线结果对比图。图2中横坐标代表驱动电流，纵坐标代表光功率，其中实线为拟合数据、虚线为实测数据。从图2中可以看出,模型几乎可以拟合出实测数据的趋势，并且拟合曲线相对测试曲线更加平滑。

Fig.2　P-I curve of measurement data and fit data at 20℃

3.2　不同温度下P-I曲线

将模型求解的参量值η,Ith,0,Rth,a1,a2,a3,a4代入(14)式,并求解温度在10℃，20℃，30℃，…，90℃下的P-I曲线,如图3所示。图3a表示完全按照模型拟合出的各P-I曲线，不考虑物理情况下光功率不能小于0。图3b显示为光功率大于0 时的P-I曲线。

由图3a可以看出,光功率随着温度和驱动电流变化的整体走势。同样在固定温度下，输出光功率随着驱动电流的增加先增后减；在固定的驱动电流下，输出光功率随着温度增加而减小。因为光功率不可能小于0，所以光功率小于0时都强制置为0。从图3b可以看出,当温度不大于60℃时才会有光功率输出。所以想要有光功率输出，则温度不能高于60℃。

目前假定电信机房里VCSEL激光器在直流输入且输出的平均光功率低于2mW时，用户的光猫无法检测到信号。分析图3b，当温度为30℃时,光功率在2mW左右浮动，所以细分位于30℃附近(24℃～32℃)的P-I曲线图，近一步精确地推测激光器正常工作的温度。

Fig.3　P-I curve of fit dataa—10℃～90℃(power unlimited)　b—10℃～60℃(power>0mW)

图4显示了24℃～32℃温度时的P-I曲线。从图上可以看出，只有在光功率大于2mW时才能保证激光器被检测到，保证用户正常使用网络;当温度大于31℃时，任何驱动电流下都不会有光功率大于2mW。所以要保证用户正常上网，电机房里VCSEL激光器工作的环境温度最多不能高于31℃。

Fig.4　P-I curve of fit data at 24℃～32℃

4　结　论

主要研究了VCSEL的P-I模型的温度特性。首先通过对P-I模型进行建模，并通过简化参量及引入U-I特性曲线来优化分析模型。模型使用LM算法计算模型参量。根据在20℃下的实测数据，计算出模型的拟合参量,最终在20℃下采用P-I模型拟合P-I曲线效果非常好，与实测的P-I曲线非常吻合。同时根据建立的模型对不同温度下的P-I曲线进行了分析，指出了用户可正常使用网络的温度范围。