高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于迭代算法的半导体光放大器动态模型构建

徐贵勇 胡立发 邓灿冉 张士勋 楚广勇

引用本文:
Citation:

基于迭代算法的半导体光放大器动态模型构建

    作者简介: 徐贵勇(1996-), 男, 硕士研究生, 现主要从事光纤通信方面的研究.
    通讯作者: 胡立发, hulifa@jiangnan.edu.cn
  • 基金项目:

    国家重点研发计划资助项目 2017YFC1502203

    中央高校基础研究专项经费资助项目 JUSRP11829

    国家自然科学基金面上资助项目 61475152

    江苏省自然科学基金资助项目 BK20170199

    江苏省"六大人才高峰"资助项目 2016-GDZB-011

    中央高校基础研究专项经费资助项目 JUSRP11835

  • 中图分类号: TN929.11

Construction of dynamic model of semiconductor optical amplifier based on iterative algorithm

    Corresponding author: HU Lifa, hulifa@jiangnan.edu.cn ;
  • CLC number: TN929.11

  • 摘要: 为了分析并精确预测半导体光放大器的性能, 对InP-InGaAsP均匀掩埋的半导体光放大器建立了一种有效的数学模型, 考虑了自发辐射与受激辐射之间的关系, 实时模拟分析了偏置电流、输入功率对增益和噪声指数的影响。结果表明, 在偏置电流为120mA、输入功率为-10dBm时, 半导体光放大器的性能最佳。该模型能够对半导体光放大器的设计提供一定的借鉴。
  • Figure 1.  Structure diagram of SOA

    Figure 2.  Internal structure diagram of SOA

    Figure 3.  Model algorithm of SOA

    Figure 4.  Optical spectrum and noise spectrum of SOA output with bandwidth of 0.1nm

    Figure 5.  Spatial distribution of forward- and backward- propagating signal and noise

    a—carrier density  b—photon rate

    Figure 6.  The effect of bias current with different SOA input powers

    Figure 7.  The effect of input power on fiber-to-fiber gain and noise figure for different SOA bias currents

  • [1]

    AHSAN M S. Past, present and future of passive optical network[J]. Trends in Opto-Electro and Optical Communications, 2019, 1(3): 16-31. 
    [2]

    CLIKAN N N, AKSOY M. A review of self-seeded RSOA based on WDM PON [J]. Canadian Journal of Electrical and Computer Engineering, 2019, 42 (1): 2-9.
    [3]

    CUI L L, WANG H L, LI W, et al. Study on gain recovery time of wavelength conversion based on single-port-coupled QD-SOA [J]. Laser Technology, 2016, 40(5):742-745(in Chinese). 
    [4]

    KIM K S. A research framework for the clean-slate design of next-ge-neration optical access [J]. Fiber and Integrated Optics, 2012, 31(2): 90-110. doi: 10.1080/01468030.2012.654597
    [5]

    WEN H, CHEN X, LI G F, et al. Four-mode semiconductor optical amplifier[J].Proceedings of the SPIE, 2016, 9774:977406. 
    [6]

    KHARRAZ O M, MAHDI M A, ATIEH A, et al. Enhanced gain saturation model of non-linear semiconductor optical amplifiers[J]. IET Optoelectronics, 2018, 12(6): 7263-7268. 
    [7]

    MEEHAN A, CONNELLY M. Experimental characterization and modeling of the improved low frequency response of a current modula- ted bulk RSOA slow light based microwave phase shifter[J]. Optics Communications, 2015, 341: 241-244. doi: 10.1016/j.optcom.2014.11.087
    [8]

    ZORAN V, ROGERIO P D, SHAHPARI A. et al. Efficient dynamic modeling of the reflective semiconductor optical amplifier [J]. IEEE Journal of Selected Topics in Quantum Electronics, 2013, 19(5): 1-10. 
    [9]

    CONNELLY M. Dynamic model of slow light in a tensile-strained semiconductor optical amplifier[J]. Optical and Quantum Electronics, 2014, 46(10): 1247-1252. doi: 10.1007/s11082-013-9831-9
    [10]

    ENGEL T, RIZOU Z V. Semiconductor optical amplifier direct mo-dulation with double-stage birefringent fiber loop [J]. Applied Phy-sics, 2016, B122(6): 158-161. 
    [11]

    CLEVELAND E R, RUPPALT L B, PROKES S M. Effect of an in situ hydrogen plasma pre-treatment on the reduction of GaSb native oxides prior to atomic layer deposition[J]. Applied Surface Science, 2013, 277: 167-175. doi: 10.1016/j.apsusc.2013.04.018
    [12]

    DANIEL M, CHIN S. Near-field scanning optical microscopy probes with heat sinks for higher power operations[J]. Materials Research Express, 2015, 2(10): 105202. doi: 10.1088/2053-1591/2/10/105202
    [13]

    DAVENPORT M L, JARED C H, NICOLAS V, et al. Heterogeneous silicon/Ⅲ-Ⅴ semiconductor optical amplifiers[J]. IEEE Journal of Selected Topics in Quantum Electronics, 2016, 22(6): 78-88. doi: 10.1109/JSTQE.2016.2593103
    [14]

    BERNHARD S. Characterization and design of multifunction photonic devices for next-generation fiber-to-the-home optical network units [D].Vieuna, Austria: Austrian Institute of Technology, 2011: 36.
    [15]

    AKIF A C, ABDULLAH A. Noise figure degradation in balanced amplifiers[J]. IEEE Microwave and Wireless Components Letters, 2017, 27(9): 848-850. doi: 10.1109/LMWC.2017.2734745
    [16]

    DONG J J. High-speed all-optical signal processing based on semiconductor optical amplifiers and optical filters [D]. Wuhan: Huazhong University of Science and Technology, 2008: 18-29(in Ch-inese).
  • [1] 詹仪郑义赵玉辉徐云峰 . 不同抽运方式下掺Yb3+光纤放大器的增益和噪声特性. 激光技术, 2007, 31(1): 50-53.
    [2] 周维军王荣波李泽仁 . 分布式光纤喇曼放大器实验研究. 激光技术, 2010, 34(3): 373-376. doi: 10.3969/j.issn.1001-3806.2010.03.025
    [3] 张铁英王红星何伍福马杰 . 光通信中基于脉冲调制的联合编码调制研究. 激光技术, 2010, 34(6): 843-846. doi: 10.3969/j.issn.1001-3806.2010.06.033
    [4] 汪锋饶炯辉 . 水下无线光通信系统中可变增益放大器的实现. 激光技术, 2012, 36(1): 99-102. doi: 10.3969/j.issn.1001-3806.2012.01.025
    [5] 许江衡罗斌潘炜贾习坤汪帆 . 基于耦合腔的VCSOAs增益带宽优化. 激光技术, 2006, 30(1): 60-63.
    [6] 贾习坤罗斌潘炜姚海峰曹昌胜 . 传输矩阵法研究垂直腔半导体光放大器增益特性. 激光技术, 2005, 29(4): 377-379,406.
    [7] 秦张淼罗斌潘炜 . 垂直腔半导体光放大器增益的理论分析. 激光技术, 2006, 30(5): 452-454.
    [8] 黄光焰杨亚培戴基智李晓惠刘震成娟娟 . 增益玻璃波导光场的数值模拟. 激光技术, 2005, 29(1): 35-37.
    [9] 车继波杨亚培刘爽官周国薛辉 . Er3+/Yb3+共掺磷酸盐玻璃光纤放大器的增益综述. 激光技术, 2006, 30(1): 82-85.
    [10] 陈海燕刘威李莉 . 相对相位噪声对相干光通信系统性能影响. 激光技术, 2016, 40(1): 94-98. doi: 10.7510/jgjs.issn.1001-3806.2016.01.021
    [11] 王刚罗斌潘炜 . 垂直腔半导体光放大器带宽特性研究. 激光技术, 2007, 31(2): 137-140.
    [12] 王欣罗斌潘炜李建平 . 传输矩阵法研究MEMS可调谐垂直腔半导体光放大器. 激光技术, 2007, 31(6): 630-633.
    [13] 金伟张海涛巩马理闫平杨欣张凯姜丰 . 无线光通信中的非成像集中器. 激光技术, 2003, 27(4): 311-316.
    [14] 白菊蓉郭宇成王彦本 . 一种改进的OFDM水下可见光无线通信系统. 激光技术, 2021, 45(5): 647-653. doi: 10.7510/jgjs.issn.1001-3806.2021.05.019
    [15] 张新陆程元丽王骐 . X射线激光在柱状等离子体中的传播与放大. 激光技术, 2002, 26(3): 221-222.
    [16] 刘剑峰于思源韩琦琦高宠马晶谭立英 . 空间光通信的时间平滑实验研究. 激光技术, 2008, 32(1): 11-14.
    [17] 王博吴琼刘立奇王涛朱仁江张鹏汪丽杰 . 水下无线光通信系统研究进展. 激光技术, 2022, 46(1): 99-109. doi: 10.7510/jgjs.issn.1001-3806.2022.01.010
    [18] 任广军赵杰林姚建铨 . 光通信波段液晶双折射效应的研究. 激光技术, 2011, 35(2): 242-244. doi: 10.3969/j.issn.1001-3806.2011.02.027
    [19] 姚文明饶炯辉张晓晖熊天林于洋 . 水下无线光通信中的FDPIM性能研究. 激光技术, 2013, 37(5): 605-609. doi: 10.7510/jgjs.issn.1001-3806.2013.05.010
    [20] 方水平 . 可调RF信号光传输系统相位噪声优化设计. 激光技术, 2016, 40(4): 582-585. doi: 10.7510/jgjs.issn.1001-3806.2016.04.026
  • 加载中
图(7)
计量
  • 文章访问数:  5537
  • HTML全文浏览量:  4436
  • PDF下载量:  18
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-04-24
  • 录用日期:  2019-08-16
  • 刊出日期:  2020-03-25

基于迭代算法的半导体光放大器动态模型构建

    通讯作者: 胡立发, hulifa@jiangnan.edu.cn
    作者简介: 徐贵勇(1996-), 男, 硕士研究生, 现主要从事光纤通信方面的研究
  • 1. 江南大学 理学院, 无锡 214122
  • 2. 国家超级计算无锡中心, 无锡 214100
  • 3. 江苏省轻工光电工程技术研究中心, 无锡 214122
  • 4. 江苏省模式识别与计算智能工程实验室, 无锡 214122
基金项目:  国家重点研发计划资助项目 2017YFC1502203中央高校基础研究专项经费资助项目 JUSRP11829国家自然科学基金面上资助项目 61475152江苏省自然科学基金资助项目 BK20170199江苏省"六大人才高峰"资助项目 2016-GDZB-011中央高校基础研究专项经费资助项目 JUSRP11835

摘要: 为了分析并精确预测半导体光放大器的性能, 对InP-InGaAsP均匀掩埋的半导体光放大器建立了一种有效的数学模型, 考虑了自发辐射与受激辐射之间的关系, 实时模拟分析了偏置电流、输入功率对增益和噪声指数的影响。结果表明, 在偏置电流为120mA、输入功率为-10dBm时, 半导体光放大器的性能最佳。该模型能够对半导体光放大器的设计提供一定的借鉴。

English Abstract

    • 基于光纤接入技术的宽带接入网技术发展迅速,用户端数据业务和视频业务已经成为用户宽带消费的主体,未来的高宽带技术急需一种大容量、高速率、低成本的接入方案[1],5G网络需求的超密集接入网技术成为目前研究的热点。集成化的光网络单元利用分布反馈激光器和半导体光放大器(semiconductor optical amplifier,SOA)集成作为一个上行信号源,而SOA性能的优劣将直接影响到终端用户的网络体验,因此对SOA进行精确预测是非常重要的[2]

    • SOA是双端口行波放大器,结构如图 1所示,其输入、输出的两个端面均镀有增透膜。SOA是在有外部电流注入情况下,使外来光子获得增益和放大的器件[3]。SOA的工作带宽一般比较宽,工作波段为1550nm时,其增益带宽一般在60nm左右。SOA有源区长度一般为1mm,厚度一般为0.4μm,中心有源区宽度约为0.4μm[4]

      Figure 1.  Structure diagram of SOA

    • 为分析并精确预测SOA的性能,需要对SOA进行数学建模。最近,WEN等人提出了一种四模式SOA(four mode SOA, FM-SOA)的理论模型,该模型的不足在于高阶模态模型对空间非均匀性更加敏感[5]。KHARRAZ等人提出了一种非线性半导体光放大器增益饱和增强模型,然而该模型没有考虑增益与波长的关系,同时忽略了放大自发辐射(amplified spontaneous emission, ASE)噪声[6],它们均限制了SOA建模的操作条件范围。作者以SOA宽带理论模型为基础,提出一种宽带InP-InGaAsP的SOA模型,通过建立合适的增益模型,将其扩展为具有量子阱有源区的SOA模型,放大器内的自发辐射由行波功率方程模拟,并忽略自发信号的相位。在该模型设计中,阐明了自发辐射与受激辐射之间的关系,分析了增益和噪声与偏置电流、输入功率之间的关系。

    • SOA是一种双端口行波放大器,其有源区是由InGaAsP直接禁带体材料制作的,呈脊行波导结构,周围材料是InP质半导体,有源区掩埋在InP质基底内,如图 2所示[7]

      Figure 2.  Internal structure diagram of SOA

      图 2a为SOA内部结构侧视图,图 2b为内部结构俯视图。SOA中心有源区的长度为Lc,有源区中心区域横向锥度的长度为Lt,因此SOA建模的平均长度可由L=Lc+Lt近似给出[8]。在放大器端面,横向锥度将光学限制因子Γ从其中心有源区域逐渐减小到0,这使得模型横向轮廓得以扩展,从而改善了单模光纤的输入和输出耦合效率。

      InGaAsP直接带隙体材料活性区域具有的材料增益系数gm(ν, n)[9]为:

      $ \begin{array}{*{20}{c}} {{g_{\rm{m}}}\left( {\nu ,n} \right) = \frac{{{c^2}}}{{4\sqrt 2 {{\rm{ \mathsf{ π} }}^3}/2n_1^2\tau {\nu ^2}}}{{\left[ {\frac{{2{m_{\rm{e}}}{m_{\rm{h}}}}}{{\hbar \left( {{m_{\rm{e}}} + {m_{\rm{h}}}} \right)}}} \right]}^{3/2}} \times }\\ {\int_0^\infty {\left\{ {\sqrt {{\nu ^\prime } - \frac{{{E_{\rm{g}}}\left( n \right)}}{h}} \left[ {{f_{\rm{c}}}\left( {{\nu ^\prime }} \right) - {f_{\rm{v}}}\left( {{\nu ^\prime }} \right)} \right] \times } \right.} }\\ {\left. {\left[ {\frac{{2{T_0}}}{{1 + {{\left( {2{\rm{ \mathsf{ π} }}{T_0}} \right)}^2}{{\left( {{\nu ^\prime } - \nu } \right)}^2}}}} \right]} \right\}{\rm{d}}{\nu ^\prime }} \end{array} $

      (1)

      式中, c为真空中的光速,ν为光频率,ν′表示激发后的光频率,n1表示有源区的折射率,h为普朗克常数,ħ为普朗克常数的1/2π,τ表示辐射载流子的复合寿命,me, mh分别为导带和价带的空穴有效质量,n表示导带的载流子密度, fcfv分别为能带和价带的费米能级分布, Eg表示带隙能量, T0表示电子与单色场的相干相互作用的平均寿命。

    • 在该模型中,Ns信号在发生耦合损耗之前被注入了光频νk和能量Pin, k,在嵌入式波导的辅助下,信号通过放大器传播,并在相对的小平面上射出。SOA模型的建立基于一组耦合微分方程,该方程描述了放大器内部变量之间的相互作用,即载流子密度和光子速率。在放大器中,每个输入信号引起的场的空间变化分量可以分解成两个复杂的行波Es, k+Es, k-来表示在z轴正方向和负方向上的传播,Es, k+Es, k-遵循行波方程[10]

      $ \begin{array}{*{20}{c}} {\frac{{{\rm{d}}E_{{\rm{s}},k}^ + \left( z \right)}}{{{\rm{d}}z}} = \left\{ {{\rm{i}}{\beta _k} + \frac{1}{2}\left[ {\mathit{\Gamma }{g_{\rm{m}}}\left( {{\nu _k},n} \right) - } \right.} \right.}\\ {\left. {\left. {\alpha \left( n \right)} \right]} \right\}E_{{\rm{s}},k}^ + \left( z \right)} \end{array} $

      (2)

      $ \begin{array}{*{20}{c}} {\frac{{{\rm{d}}E_{{\rm{s}},k}^ - \left( z \right)}}{{{\rm{d}}z}} = \left\{ {{\rm{i}}{\beta _k} + \frac{1}{2}\left[ {\mathit{\Gamma }{g_{\rm{m}}}\left( {{\nu _k},n} \right) - } \right.} \right.}\\ {\left. {\left. {\alpha \left( n \right)} \right]} \right\}E_{{\rm{s}},k}^ - \left( z \right)} \end{array} $

      (3)

      式中,α(n)表示材料损耗系数, $\mathrm{i}=\sqrt{-1}$, k=1, 2, …, Ns, 行波方程受边界条件影响:

      $ E_{{\rm{s}},k}^ + \left( 0 \right) = \left( {1 - {r_1}} \right){E_{{\rm{in}},k}} + {r_1}{E_{{\rm{s}},k}}\left( 0 \right) $

      (4)

      $ E_{{\rm{s}},k}^ - \left( L \right) = {r_2}E_{{\rm{s}},k}^ + \left( L \right) $

      (5)

      输入面左侧的输入信号字段是:

      $ {E_{{\rm{in}},k}} = \sqrt {\frac{{{\eta _{{\rm{in}}}}{P_{{\rm{in}},k}}}}{{h{\nu _k}}}} $

      (6)

      输出面右侧的输出信号字段是:

      $ {E_{{\rm{out}},k}} = \left( {1 - {r_2}} \right)E_{{\rm{s}},k}^ + \left( L \right) $

      (7)

      耦合损耗后的输出信号功率为:

      $ {P_{{\rm{out}},k}} = h{v_k}{\eta _{{\rm{out}}}}{\left| {{E_{{\rm{out}},k}}} \right|^2} $

      (8)

      式中, ηin, ηout分别为输入和输出的耦合效率,振幅反射率系数$r_{1}=\sqrt{R_{1}}, r_{2}=\sqrt{R_{2}}$,R1R2分别为输入面和输出面的反射率。信号传播系数$\beta_{k}=\frac{2 \pi n_{\mathrm{eq}} \nu_{k}}{c}$, neq是给定放大器波导的等效指数[11]

      $ {n_{{\rm{eq}}}} = {\left[ {\left({n_1^2 - n_2^2} \right)\frac{\mathit{\Gamma }}{{2 - \mathit{\Gamma }}} + n_2^2} \right]^{1/2}} $

      (9)

      式中, n2是有源区周围InP材料的折射率。

    • 信号放大还取决于放大器自身产生的自发辐射噪声的数量,这是因为噪声功率不仅参与了可用载流子的排放,而且还有助于使增益饱和。然而无需将自发辐射视为相干信号,因为它在相对宽的波段上连续分布,相邻波长分量之间具有随机相位。当存在反射小平面时,自发辐射出的噪声表明纵向腔模式的存在,由于这个原因,可以假设噪声光子仅存在于谐振腔整数倍的离散频率处,这些频率由下式给出:

      $ {\nu _j} = {\nu _{\rm{c}}} + \Delta {\nu _{\rm{c}}} + j{K_{\rm{m}}}\Delta {\nu _{\rm{m}}} $

      (10)

      式中, j=0, …, Nm-1;$\nu_\text{c}=\frac{E_\text{g, 0}}{h}$; Δνc为与共振ν0相匹配的频率偏移量;Eg, 0表示无注入载流子的带隙能量;Km, Nm是两个正整数,其值取决于SOA增益带宽和模型方程数值解所要求的精度。

      纵模频率间隔为:

      $ \Delta {\nu _{\rm{m}}} = \frac{c}{{2\int_0^L {{n_{{\rm{eq}}}}} {\rm{d}}z}} $

      (11)

      可以证明,在两个相邻谐振腔上对相干信号求平均与对在行波功率(或光子速率)方程方面相干处理信号是相同的。自发辐射可用功率来描述,信号必须用具有明确幅度和相位的波来处理。Nj+Nj-被定义为特定偏振的自发辐射光子率,服从以下行波方程:

      $ \begin{array}{*{20}{c}} {\frac{{{\rm{d}}N_j^ + \left( z \right)}}{{{\rm{d}}z}} = \left[ {\mathit{\Gamma }{g_{\rm{m}}}\left( {{\nu _j},n} \right) - } \right.}\\ {\left. {\alpha \left( n \right)} \right]N_j^ + + {R_{{\rm{sp}}}}\left( {{\nu _j},n} \right)} \end{array} $

      (12)

      $ \begin{array}{*{20}{c}} {\frac{{{\rm{d}}N_j^ - \left( z \right)}}{{{\rm{d}}z}} = \left[ {\mathit{\Gamma }{g_{\rm{m}}}\left( {{\nu _j},n} \right) - } \right.}\\ {\left. {\alpha \left( n \right)} \right]N_j^ - + {R_{{\rm{sp}}}}\left( {{\nu _j},n} \right)} \end{array} $

      (13)

      式中,Rsp表示耦合到噪声信号中的自发辐射噪声。

      在符合边界条件的情况下:

      $ \left\{ {\begin{array}{*{20}{l}} {N_j^ + (0) = {R_1}N_j^ - (0)}\\ {N_j^ - (L) = {R_2}N_j^ + (L)} \end{array}} \right. $

      (14)

      理想的放大器没有增益饱和(这意味着整个放大器的载流子密度恒定)、零损耗系数、小平面反射率和耦合损耗。其材料增益系数为gm=gm′(>0),在这种情况下Nj+可从下式解得:

      $ \frac{{{\rm{d}}N_j^ + }}{{{\rm{d}}z}} = \mathit{\Gamma }g_{\rm{m}}^\prime \left( {{\nu _j},n} \right)N_j^ + + {R_{{\rm{sp}}}}\left( {{\nu _j},.n} \right) $

      (15)

      同时在KmΔνm的中心频率带νj中,在单极化处给出输出噪声功率:

      $ {N_{{\rm{out}}}} = \frac{{G - 1}}{{\mathit{\Gamma }g_{\rm{m}}^\prime \left( {{\nu _j},n} \right)}}h{\nu _j}{R_{{\rm{sp}}}}\left( {{\nu _j},n} \right) $

      (16)

      式中, G是放大器增益,用行波功率方程来描述Nj+Nj-,并假设间距中所有的自发光子都处于共振频率中。在实际装置中,源自Rsp注入的自发光子将均匀分布在KmΔνm上。如果在νj的单程增益为Gs,则间隔Δνm的频率信号增益为[9]

      $ G(\nu ) = \frac{{\left( {1 - {R_1}} \right)\left( {1 - {R_2}} \right){G_{\rm{s}}}}}{{{{\left( {1 - \sqrt {{R_1}{R_2}} {G_{\rm{s}}}} \right)}^2} + 4{G_{\rm{s}}}\sqrt {{R_1}{R_2}} {{\sin }^2}\phi }} $

      (17)

      式中, 单通相移为:

      $ \phi = \frac{{2{\rm{ \mathsf{ π} }}\nu }}{c}\int_0^L {{n_{{\rm{eq}}}}} \left( z \right){\rm{d}}z $

      (18)

      共振时信号增益为:

      $ G\left( {{v_j}} \right) = \frac{{\left( {1 - {R_1}} \right)\left( {1 - {R_2}} \right){G_{\rm{s}}}}}{{{{\left( {1 - \sqrt {{R_1}{R_2}} {G_{\rm{s}}}} \right)}^2}}} $

      (19)

      假设放大器的噪声输入光谱密度σin均匀分布在以νj为中心的Δνm处,那么Δνm中的总输出噪声为:

      $ \begin{array}{*{20}{c}} {{N_{{\rm{out}}}} = {\sigma _{{\rm{in}}}}\int_{\nu = {\nu _j} - \Delta {\nu _{\rm{m}}}/2}^{\nu = {\nu _j} + \Delta {\nu _{\rm{m}}}/2} G (\nu ){\rm{d}}\nu = }\\ {{\sigma _{{\rm{in}}}}\Delta {\nu _{\rm{m}}}\frac{1}{{\rm{ \mathsf{ π} }}}\int_0^{\rm{ \mathsf{ π} }} G \left( \phi \right){\rm{d}}\phi } \end{array} $

      (20)

      如果输入噪声功率集中在νj处,输出噪声光子率将为:

      $ {N_{{\rm{res}},{\rm{out}}}} = {\sigma _{{\rm{in}}}}\Delta {\nu _{\rm{m}}}G\left( {{\nu _j}} \right){K_j} $

      (21)

      式中, Kj使得Nres, out=Nout,则其值为:

      $ {K_j} = \frac{{\int_0^{\rm{ \mathsf{ π} }} G (\phi ){\rm{d}}\phi }}{{{\rm{ \mathsf{ π} }}G\left( {{\nu _j}} \right)}} = \frac{1}{{\rm{ \mathsf{ π} }}}\int_0^{\rm{ \mathsf{ π} }} {\frac{{{\rm{d}}\phi }}{{1 + {\gamma ^2}{{\sin }^2}\phi }}} $

      (22)

      式中,$ \gamma=\frac{4 G\left(\nu_{j}\right) \sqrt{R_{1} R_{2}}}{\left[1-\sqrt{R_{1} R_{2}} G_{\mathrm{s}}\left(\nu_{j}\right)\right]^{2}} $, 且Kj零面反射率的单位一致。

    • 载流子密度服从速率方程:

      $ \begin{array}{*{20}{c}} {\frac{{{\rm{d}}n(z)}}{{{\rm{d}}t}} = \frac{I}{{edLW}} - R(n(z)) - \frac{\mathit{\Gamma }}{{dW}}\left\{ {\sum\limits_{k = 1}^{{N_{\rm{s}}}} {{g_{\rm{m}}}} \left[ {{\nu _k},n(z)} \right] \times } \right.}\\ {\left. {\left[ {N_{{\rm{s}},k}^ + (z) + N_{{\rm{s}},k}^ - (z)} \right]} \right\} - }\\ {\frac{{2\mathit{\Gamma }}}{{dW}}\left\{ {\sum\limits_{j = 0}^{{N_{\rm{m}}} - 1} {{g_{\rm{m}}}} \left[ {{\nu _j},n(z)} \right]{K_j}\left[ {N_j^ + (z) + N_j^ - (z)} \right]} \right\}} \end{array} $

      (23)

      式中, I是偏置电流, W是中心有源区宽度。在(23)式中,假设所有偏置电流仅通过有源区而不通过周围的InP区,偏置电流在有源区宽度上将均匀分布。

      重组速率项R(n)为[12]

      $ R\left( n \right) = {R_{{\rm{rad}}}}\left( n \right) + {R_{{\rm{n}} - {\rm{rad}}}}\left( n \right) $

      (24)

      式中, RradRn-rad(n)分别是辐射和非辐射载流子复合率,两者都可以表示为载流子宽度n的多项式函数:

      $ {R_{{\rm{rad}}}}(n) = {A_{{\rm{rad}}}}(n) + {B_{{\rm{rad}}}}{n^2} $

      (25)

      $ {R_{{\rm{n}} - {\rm{rad}}}}(n) = {A_{{\rm{n}} - {\rm{rad}}}}n + {B_{{\rm{n}} - {\rm{rad}}}}{n^2} + {C_{{\rm{aug}}}}{n^3} + {D_{{\rm{leak}}}}{n^{5.5}} $

      (26)

      式中,AradBrad为线性双分子辐射复合系数[13], An-rad表示线性非辐射复合系数[14]Bn-rad表示非辐射双分子重组,Caug表示俄歇复合系数,由于泄漏效应,Dleak表示重组系数,在具有高载流子密度的SOA中,俄歇复合和泄漏效应可能是十分明显的,对于该模型所考虑的特定器件,假设从有源区到周围InP区的载流子泄漏忽略不计。

    • 由于SOA模型方程不能解析求解,因此需要进行数值求解,该算法的流程图如图 3所示。

      Figure 3.  Model algorithm of SOA

      在该模型中, 为了预测SOA特性,使用一种调整载流子密度的算法。放大器模型被分割为i=1~Nz部分。Q(i)表示第i部分,W(i)表示一个小于整体的分量, 使得整个放大器中的Q(i)值接近零,然后在每节运算中计算Q(i)来进行数值求解。算法的第1步是将信号场和自发辐射光子率初始化为零,其次利用(2)式、(3)式和(12)式、(13)式的有限差分解估计信号场和噪声光子密度。迭代持续进行,直到整个SOA的信号场、噪声光子率和载流子密度在连续迭代之间的百分比变化小于期望的容忍度。当迭代停止时,计算出输出自发辐射功率谱密度,并计算诸如信号增益,噪声系数和输出自发噪声功率之类的参量,该算法在多种操作条件下具有良好的收敛性和稳定性。

      SOA可以用于同时放大不同波长的信号,图 4为两个波长为1550nm和1555nm的输入信号功率系数谱以及信号波长1550nm时,偏置电流为120mA下SOA的噪声频谱图。当在SOA中放大调制信号时,必须注意放大器以线性方式操作,否则将发生交叉增益调制。

      Figure 4.  Optical spectrum and noise spectrum of SOA output with bandwidth of 0.1nm

      噪声系数Fn的计算公式为[15]

      $ {F_{\rm{n}}} = 10\lg \left[ {\frac{{{\sigma _{{\rm{ASE}}}}(\nu )}}{{h\nu G(\nu )}} + \frac{{{\eta _{{\rm{out}}}}}}{{G(\nu )}}} \right] $

      (27)

      式中, ηout表示输出耦合损耗,假定放大器输入信号噪声是可忽略的,则G(ν)是光频ν处光纤间的增益。

      SOA载流子密度、前向和后向传播的总ASE光子速率以及低输入信号功率和高输入信号功率下信号光子速率的空间分布情况如图 5所示。由于低的小平面反射率,后向传播的信号可以忽略不计。在低输入功率时,载流子密度空间分布具有对称性,在SOA的中心处达到峰值,并朝向输入和输出小平面拖尾,如图 5a所示,这是因为ASE在这些区域达到峰值,对于图 5b,在高输入功率下,载流子密度空间分布变得更不对称, 对于空间分布的分析可用于辅助SOA的设计。

      Figure 5.  Spatial distribution of forward- and backward- propagating signal and noise

      图 6中描述了SOA在输入信号波长为1550nm下增益和噪声与偏置电流、输入功率之间的关系。对于图 6a,在偏置电流较小的情况下,SOA工作在非饱和状态,即线性区,增益随着偏置电流近似呈线性变化,随着偏置电流的增大,增益变化缓慢,并且逐渐达到上限值,在120mA左右达到饱和状态,即非线性区。在偏置电流一定的情况下,增益随着输入功率的增大而减小。噪声变化情况如图 6b所示,偏置电流较小时,噪声急剧下降,在偏置电流约为30mA时,SOA开始工作,并在120mA附近变化缓慢,趋于稳定。在偏置电流一定的情况下,输入功率增大时噪声呈减小状态。从图中可以看出,偏置电流在120mA左右时,增益饱和,噪声较低,SOA性能达到最优。

      Figure 6.  The effect of bias current with different SOA input powers

      图 7可以看出, 增益随着输入功率的增大呈下降趋势:输入功率较小时,增益下降缓慢; 当功率增加到-20dBm时,下降趋势非常明显,反映了SOA的增益饱和性。同样的,噪声在较小输入功率下变化缓慢:当输入功率大于-20dBm时,下降趋势明显;当输入功率大于-5dBm时,可以看出,噪声随着输入功率的增大而增大,不再呈下降趋势,造成这一反常现象是由于光纤非线性效应[16]。在输入功率一定时,增益随着偏置电流的增大而增大,噪声则随着偏置电流的增大而减小。从图中可以看出,输入功率在-10dBm左右,SOA的性能最优,此时的增益相对稳定且较高,噪声稳定且较小。

      Figure 7.  The effect of input power on fiber-to-fiber gain and noise figure for different SOA bias currents

    • 通过分析SOA自发与受激辐射之间的关系,建立了一种有效的SOA理论模型。在对SOA的建模过程中,通过求解自发辐射系数的数值解,采用多阶迭代算法来分析SOA的实时增益以及噪声指数。从结果可知,SOA的增益和噪声指数在输入光功率不变时,随偏置电流的增大而逐渐趋于饱和, 在偏置电流一定的情况下,SOA的增益以及噪声指数随输入功率的增大而增大, 在偏置电流为120mA、输入功率为-10dBm时,SOA的性能最佳。该模型可用于研究不同材料和几何参量对SOA特性的影响,为改善半导体光放大器模型提供了一定的借鉴意义。

参考文献 (16)

目录

    /

    返回文章
    返回