Study on the Mechanism of Mach-Zehnder Interferometer for Suppressing the Delay Time Signature
1 引言
混沌激光具有高带宽、大幅度、类噪声以及对初值条件敏感等特性,被广泛应用于混沌保密通信[1-4]、高精度雷达测距、混沌光时域反射仪、光学相干层析、高速随机数发生器等领域[5-10]。作为B类激光器的半导体激光器,通过附加自由度(如:光反馈、外部光注入)才会具有丰富的非线性动力学特性,而外腔光反馈半导体激光器因其结构简单、动力学特性丰富、稳定性好,易于调节和控制,成为广泛使用的混沌光源。而这种光反馈半导体激光器具有明显的时延特性,即半导体激光器输出的混沌激光会呈现与反馈腔长或注入延时有关的弱周期性。弱周期性的存在会降低混沌通信系统的安全性,削弱产生的随机数的无序性。因此,混沌半导体激光器时延信息的隐藏或抑制已成为混沌应用中亟待解决的问题。
Rontani等[11]首先对单光反馈系统的混沌延迟时间特征问题进行了理论分析,在低反馈强度下,调节注入偏置电流使激光器的弛豫振荡周期接近于外腔反馈延迟时间时,激光器弛豫振荡频率将掩盖反馈时延特征而被隐藏,但是此时混沌信号的维度和带宽很难得到提高。针对上述问题,提出了多种方法来优化光反馈半导体激光器,从而达到抑制或隐藏其时延特性的效果。例如:采用布拉格光纤光栅的分布反馈[12]或啁啾光纤光栅的色散反馈抑制时延信息[13],滤波光反馈隐藏半导体激光器的时延信息[14],利用双激光器互耦合[15]或三个激光器组成的系统(两个单光反馈的激光器产生的混沌激光注入到第三个激光器[16]或级联模式[17])来抑制时延信息。但是,研究较多的是双光反馈或双光注入的模式,2009年西南大学Wu等[18]研究双光反馈半导体激光器对时延特征的抑制,2013年Wang等[19]采用延迟自相干来抑制延时特征。2013年Xiang等[20]提出单光反馈激光器产生的混沌激光双光注入到一个从激光器中来抑制时延特征信息,次年提出相位调制双路反馈抑制时延信息[21],2016年本课题组研究双光反馈激光器产生的混沌激光注入到一个从激光器来抑制时延信息[22-23]。然而,以上研究都只是说明通过调节相关参数可以实现延时特征峰的抑制,而对这些抑制方法的内在机理并没有深入分析。因此在前期工作的基础上经过深入分析发现,双光反馈或双光注入抑制延时特征峰方案本质上都存在马赫-曾德尔干涉仪(MZI)的结构,本文通过两种方案的对比,总结了MZI在抑制时延特性中的普遍规律,分析了MZI在抑制时延特性中的作用与机理,以期揭示这些方案抑制延时特性的内在机制。
2 理论模型
为了分析MZI抑制延迟特性的作用机理,对比两种抑制延时特征峰的常见方案。方案一:如图1(a)所示,由一个主激光器和从激光器构成,主激光器是一个单光反馈激光器,产生的混沌激光经分束器②分为两路光Ⅰ、Ⅱ(Ⅰ、Ⅱ延时不同),光束Ⅰ、Ⅱ再经耦合器③耦合后注入到从激光器,从而实现双光注入,相当于通过由分束器②、两路光纤和耦合器③构成的MZI把光注入到从激光器中,激光从端口1输出。方案二:如图1(b)所示,为双光反馈方案,激光器输出的光经分束器②分为两路光Ⅰ、Ⅱ(Ⅰ、Ⅱ延时不同),光束Ⅰ、Ⅱ再经耦合器③耦合后反馈到激光器,相当于通过MZI把光反馈到激光器中。
图 1. 模拟原理图。(a)双光注入方案;(b)双光反馈方案
Fig. 1. Schematic of simulation. (a) Double light injection scheme; (b) dual optical feedback scheme
下载图片 查看所有图片
在上述两种方案中,直接分析MZI对延时特性的抑制机理比较复杂,为此分析了单光反馈主激光器产生的混沌激光只经MZI作用的方案,相当于方案一去掉了从激光器,激光直接从端口③输出,如图1(a)所示。假设从主激光器输出光的电场强度表示为E(t)=2A(t)exp[i(ωt+ϕ)],经MZI的作用输出的光强表示为
式中τd为在MZI中由于两臂Ⅰ、Ⅱ存在光程差产生的时间延迟。为了便于表述,把此分析方案称为方案三。
2.1 速率方程模型
对于带有光反馈的半导体激光器,采用著名的Lang-Kobayashi速率方程模拟动力学特性[24]。本文中,为使模拟的激光器输出与真实的激光器输出相类似,在模拟时需考虑增益饱和效应的影响[25]。
其中
式中GN为微分增益,τp为光子寿命,τn为载流子寿命,ε为增益饱和系数,E和Ei为电场强度,ϕi和ϕ为电场相位,Ni和N为载流子浓度。(2)~(4)式为双光反馈的激光器,(2)式和(3)式的后两项代表反馈项,kf1和kf2为反馈强度,τ1和τ2为反馈延迟时间。(5)~(7)式为双光注入的从激光器,kin1和kin2为注入强度,τin1和τin2为注入延迟时间,Δf=Δω/2π=(ωi-ω)/2π为失谐频率。利用(2)~(7)式对方案一进行模拟,利用(2)~(4)式对方案二和方案三进行模拟。
2.2 时延信息的提取
通常采用自相关函数、互信息、填充因子分析以及局部线性模型等方法来分析光反馈混沌系统的延迟时间特征。采用自相关函数进行分析,通过混沌时间序列的自相关函数获得外腔反馈延迟时间的信息:
式中I(t)=E2(t)为混沌激光光强,τ为延迟时间。
为了量化对时延特征的隐藏,引入ρ表征自相关函数在区间W(τ)上的峰值,ρ值越小,对时延特征的抑制越好。由于实际的延迟时间与理论的延迟时间有一定的偏差,设定W(τ)=[τ-1 ns,τ+1 ns],τ=10 ns。
图 2. 单光反馈激光器输出混沌信号的延时特征峰ρ随反馈强度kf1的变化
Fig. 2. Peak value ρ of the delay signature for chaos laser with single path optical feedback as function of feedback strength kf1
下载图片 查看所有图片
3 结果与讨论
利用4阶龙格-库塔算法对(2)~(7)式进行数值求解,模拟中采用的参数为[26]:阈值电流Jth=[N0+1/GNτP]/τn,微分增益GN=8.4×10-13 m3·s-1,透明载流子浓度N0=1.4×1024 m-3,载流子寿命τn=2.04×10-9 s,光子寿命τp=1.927×10-12 s,增益饱和系数ε=2.5×10-23,线性展宽因子α=5.0。本文中不考虑失谐的影响,令Δf=0。
3.1 MZI抑制延时特性的规律
在方案一中,主激光器为单光反馈激光器,模拟时令kf2=0,设定τ1=10 ns。为达到对延时特征峰的最佳抑制,主激光器输出的混沌信号应具有最小的延时特征峰,此时反馈强度kf1=11.4 ns-1,如图2所示。注入强度kin1=kin2=16 ns-1,注入延迟时间τin1=20 ns,τin2=τin1+τd,通过调节MZI中光路的臂长来调节τd的变化。在方案二中,设定τ1=10 ns,τ2=τ1+τd,两路反馈光的反馈强度取相等[18],反馈强度[22]kf1=kf2=6.9 ns-1,此时能很好地抑制延时特征峰。在方案三中,参数的选择与方案一相同,输出的光强利用(1)式模拟计算获得。利用自相关函数对三种方案输出的光强进行模拟分析,延时特征峰ρ随τd的变化如图3所示。图3(a)~(c)所示为延迟特征峰随τd变化的整体情况,图3(d)~(f)所示为τd较小(两臂差别很小)时的情况。
图 3. 延时特征峰ρ随τd的变化。(a)(b)方案一;(c)(d)方案二;(e)(f)方案三
Fig. 3. Peak value ρ of the delay signature as function of τd. (a)(b) Scheme 1; (c)(d) scheme 2; (e)(f) scheme 3
下载图片 查看所有图片
由于反馈延迟时间τ1=10 ns,为了全面展示MZI中τd对延时特征的影响,设定τd变化的范围为[0,12 ns],τd变化的间隔为Δτd=0.1 ns,如图3(a)~(c)所示。当τd=0时,MZI中的两臂相等,对于方案一,即一个单光反馈的主激光器产生的混沌激光单光注入到一个从激光器,方案二相当于一个普通的单光反馈的激光器,而方案三为一个单光反馈的主激光器,此时延时特征都很强。当τd=10 ns时,对于方案一和方案三,由于τ1=10 ns,输出激光的电场强度E(t)与E(t-τd)具有很强的相关性,经过MZI作用,两束光叠加后结果并不能抑制延时特征。对于方案二,讨论当τ1=10 ns、τ2=20 ns的情况,相当于一个外腔长等于另一个外腔长的二倍,此时延时特征不能得到抑制,与文献[
16]的实验结果一致。在0<τd<10 ns的很大区域内,对于上面三种方案,延时特征得到较大的抑制,基本上处于一个稳定值,但是方案一明显好于方案二和方案三。
比较特殊的情形是:τd在区域(0,0.2 ns)变化时,延时特征峰ρ随τd变化显著。当τd在10 ns附近时,有着相似的变化规律。为了清晰地刻画这些区域的变化规律,对τd在(0,0.2 ns)区域进行了更精细的模拟,Δτd=0.002 ns,如图3(d)~(f)所示。对于方案一和方案三,ρ值在τd接近于0的区域剧烈振荡,ρ值随τd呈现一定的周期变化,虽然振荡的幅度有所差别,但是振荡的周期基本一致。对于方案二,双光反馈系统,ρ值随τd的变化规律与方案一和方案三有一定的差别,特别是在τd=0.01 ns或τd=0.02 ns时,ρ=1;τd=0.106 ns时,ρ=0.89,如图3(e)黑色箭头所示的位置。在模拟时,为了消除初值设定对混沌振荡的影响,首先尝试一路反馈光反馈,另一路反馈光断开,当系统进入到混沌振荡后再接通另一路反馈光,当τd=0.02 ns时,在双光反馈的状态下,激光器很快进入稳态,不能进行混沌激光输出,如图4(a)、(b)所示;当τd=0.106 ns时,在双光反馈的状态下,激光器处于准周期振荡状态,如图4(c)、(d)所示。除去这些特殊的位置,ρ值随τd的变化规律与方案一和方案三基本相似。综上所述,通过三种方案的对比分析,对于延时特征峰的抑制,MZI在三种方案中起主要作用。对于方案一,输出的混沌激光不仅受到MZI的作用,还受到从激光器的非线性放大;对于方案二,输出的混沌激光来自MZI与激光器本身多次循环作用的结果。因此,在这两种方案中,直接分析MZI的作用机理相当困难,而对于方案三,输出的混沌激光仅受到MZI的作用,有利于分析MZI对抑制延时特征的作用机理。
图 4. 时间序列和自相关函数。(a)(b) τd=0.02 ns;(c)(d) τd=0.106 ns
Fig. 4. Intensity time series and autocorrelation function. (a)(b) τd=0.02 ns; (c)(d) τd=0.106 ns
下载图片 查看所有图片
单光反馈激光器直接输出的混沌激光的光强时间序列,如图5(a)所示,经MZI作用后输出的混沌激光的时间序列,τd=0.08 ns,如图5(b)所示,从图中可以看出,在MZI的作用下,输出的混沌激光的光强在时域得到重新的分布。方案一输出的混沌激光的光强时间序列,τd=0.08 ns,如图5(g)所示,对比图5(b),光强变得连续光滑,因为方案一与方案三的区别是增加了一个从激光器,从激光器相当于一个非线性放大器,对第一种方案输出混沌激光进行非线性放大处理,使光强变得光滑连续。方案二输出的混沌激光的光强时间序列,τd=0.098 ns,如图5(h)所示。相应的自相关函数如图5(d)、(e)、(i)、(j)所示,通过对比发现,当τd<0.1 ns时,选择合适的延迟时间,延时特征峰得到很大的抑制。
对于方案三,由图5(c)可以看出,当τd=4 ns(τd>0.1 ns)时,混沌激光的时间序列发生了明显变化。由图5(f)可以看出,原有的延时特征峰(τ=10 ns)不仅没有得到有效的抑制,而且在τ=4,6,14 ns等处出现了新的延时特征峰。为进一步总结新出现的延时特征峰的规律,任选了τd=1,3,5 ns时三种方案的自相关函数的情况,如图6所示。同样地,在τ=10 ns的位置仍然存在明显的延时特征峰,同时出现了新的延时特征峰。
以方案三中τd=3 ns的情形进行分析,如图7所示。经MZI作用输出的混沌激光是由这两路光叠加而成,光路2相当于由光路1延迟3 ns产生,因此A与B有很强的相关性,对应于图6(a)中τd=3 ns的延时特征峰。而混沌激光中本身A与C具有相关性,对应于图6(a)中τ=10 ns的延时特征峰,所以B与C也有一定的相关性,对应于图6(a)中τ=7 ns的延时特征峰,但是峰值较低。设单光反馈激光器的反馈延时时间τ1=10 ns,则在
(m为正整数)的位置出现新的延时特征峰。比较特殊的是:在τd=τ1/2=5 ns时,此时A与B和B与C相差的时间相等,相当于两个延时特征峰的位置重合,因此在此处出现很强的延时特征峰,对应于图6(a)中τ=5 ns的延时特征峰,并且大于τ=10 ns的延时特征峰。另一种特殊的情形是:τd很小(τd<0.1 ns)时,产生新的延时特征峰的位置与原有的延时特征峰的位置基本重合,因此观察不到新的延时特征峰,如图5(e)所示。方案一的规律与方案三相同,只是整个延时特征峰都偏小,如图6(b)所示。方案二是双光反馈,反馈延时时间为τ1=10 ns,τ2=τ1+τd,在
(m1和m2为非负整数)出现延时特征峰,如图6(c)所示。τ1和τ2位置的延时特征峰随τd的变化如图8所示,在τd变化的过程中,两个位置的延时特征峰的大小基本一致。在三种方案中,除了一些特殊的位置(MZI中一个臂长是另一个臂长的1/2),τ=10 ns位置的延时特征峰值最大或基本相当,所以τ=10 ns位置的ρ(主延时特征峰)最具有代表性。
图 5. 时间序列与自相关函数。方案三(a)(d) τd=0 ns,(b)(e) τd=0.08 ns,(c)(f) τd=4 ns;(g)(i)方案一,τd=0.08 ns;(h)(j)方案二,τd=0.098 ns
Fig. 5. Time series and autocorrelation function. Curves of scheme 3 when (a)(d) τd=0 ns, (b)(e) τd=0.08 ns, (c)(f) τd=4 ns; curves of scheme 1 when (g)(i) τd=0.08 ns; curves of scheme 2 when (h)(j) τd=0.098 ns
下载图片 查看所有图片
图 6. 不同τd下的自相关函数。(a)方案三 ;(b)方案一 ;(c)方案二
Fig. 6. Autocorrelation for different τd. (a) Scheme 3; (b) scheme 1; (c) scheme 2
下载图片 查看所有图片
图 7. 分析方案三的延时特征峰
Fig. 7. Analysis of the delay signature for scheme 3
下载图片 查看所有图片
3.2 MZI抑制延时特性的作用机理分析
以方案三为研究对象,为了不失一般性,单光反馈激光器的反馈延迟时间任选两种情形,反馈延迟时间τ1=4 ns和τ1=7 ns,其他参数与前面的一致。反馈延迟时间τ1=4 ns时,如图9(a)、(b)所示,在τd=0 ns和τd=4 ns时,延时特征峰ρ很大,当0<τd<4 ns时,延时特征得到较大的抑制,特别是τd在区域(0,0.2 ns)变化时,延时特征峰ρ随τd周期性变化,整个变化的规律与τ1=10 ns的一致。反馈延迟时间τ1=7 ns时,如图9(c)、(d)所示,变化规律与τ1=4 ns相似,MZI对延时特征峰ρ的抑制具有一定的普遍规律。
图 8. 方案二中τ1=10 ns和τ2=τ1+τd位置的延时特征峰ρ随τd的变化
Fig. 8. Peak value ρ of the delay signature corresponding to τ1=10 ns and τ2=τ1+τd as function of τd in the scheme 2
下载图片 查看所有图片
图 9. 延时特征峰ρ随τd的变化。(a)(b) τ1=4 ns;(c)(d) τ1=7 ns
Fig. 9. Peak value ρ of the delay signature as function of τd. (a)(b) τ1=4 ns; (c)(d) τ1=7 ns
下载图片 查看所有图片
下面分析当τ1=10 ns、反馈强度kf1=11.4 ns-1时MZI对延时特征峰的影响因素。图10(a)所示为混沌激光经MZI作用后的延时特征峰ρ随τd的变化情况。当τd=0时,ρ对应于单光反馈激光器的延时特征峰,设ρ0=0.267。随着τd的增大,ρ逐渐趋于一个稳定值ρ=0.125,如图10(a)中的黑色点划线所示。对于(1)式,前两项是MZI两臂混沌信号的光强,最后一项是干涉项。当忽略干涉项时,I(t)=A2(t)+A2(t-τd)。图10(b)所示为混沌信号的延时特征峰ρ随τd的变化情况,当τd在区间(0,0.4 ns)变化时,ρ有一定的起伏,因为此时A(t)与A(t-τd)具有较强的相关性,如图11(a)所示;当τd>0.4 ns时,A(t)与A(t-τd)的相关性逐渐消失,在随τd增大的过程中,ρ逐渐趋于ρ0,如图10(b)中的黑色点划线所示。说明非相干项对抑制延时特征峰的贡献较小,仅在区间(0,0.4 ns)对ρ值有较小的影响。当忽略非相干项时,I(t)=2A(t)A(t-τd)cos[ωτd+ϕ(t)-ϕ(t-τd)],图10(c)所示为混沌信号的延时特征峰ρ随τd的变化情况,随着τd的增大,ρ同样逐渐趋于一个稳定值ρ=0.125,如图10(c)中的黑色点画线所示。对比图10(a)与图10(c)可以发现,当τd大于一定值时,对延时特征峰的抑制主要由相干项决定。当τd在(0,0.06 ns)变化时,ρ随τd周期性的剧烈变化,但是对比图10(a),其变化的周期明显不同。综上所述,(1)式中的所有项对ρ值都有影响。在τd<0.1 ns时,ρ表现出周期性的起伏变化。针对这个问题,对cos[ωτd+ϕ(t)-ϕ(t-τd)]统计求平均,令该值为α,α随τd的变化如图10(d)所示。对比图10(a),两者变化的周期基本一致。在τd<0.1 ns时,ϕ(t)与ϕ(t-τd)具有较强的相关性,如图11(b)所示,τd的变化对α起主要作用。τd越小,ϕ(t)与ϕ(t-τd)的相关性越强,α的幅值越大。当τd>0.1 ns时,ϕ(t)与ϕ(t-τd)的相关性逐渐消失,α趋于0,对应ρ趋于稳定,不随τd变化。方案三中经MZI作用输出的混沌激光的平均功率Pave随τd 的变化如图12所示,对比图12与图10(d),具有相同的振荡周期。当τd <0.1 ns时,两路光具有较强的相关性,Pave随τd周期性振荡,此时两路光的叠加表现为相干叠加;当τd >0.1 ns时,两路光相关性逐渐消失,Pave不再随τd变化,而是趋于一个稳定值,MZI中一个臂上的平均光强Pave1如图12中的红色虚线所示,由于两个臂上的光相同,Pave=2Pave1,所以此时光的叠加为非相干叠加。综上所述,当MZI的两臂接近相等时(τd <0.1 ns),两路光的叠加为相干叠加,输出的混沌激光的延时特征峰随τd周期性变化;当MZI的两臂臂长相差较大且不成比例时,两路光的叠加为非相干叠加,延时特征峰得到一定的抑制同时不随τd变化,此时对延时特征峰的抑制也主要由干涉项决定。
图 10. 影响延时特征峰ρ的因素
Fig. 10. Factors influencing the delay time signature peak ρ
下载图片 查看所有图片
图 11. 自相关函数图。(a)混沌光强的自相关函数;(b)混沌光相位的自相关函数
Fig. 11. Autocorrelation function diagram. (a) Autocorrelation function for chaotic light intensity; (b) autocorrelation function for chaotic light phase
下载图片 查看所有图片
图 12. 混沌光的平均功率Pave随τd的变化
Fig. 12. Average power of the output of chaos versus τd
下载图片 查看所有图片
最后对方案二的特殊位置给出解释,如图3(d)所示,在τd 较小的区域(τd <0.1 ns),经MZI作用反馈到激光器的的反馈光功率也要随着τd变化,而反馈光功率的大小直接影响着激光器的状态。由于在方案二中,混沌激光是MZI与激光器本身多次循环作用的结果,情形相对复杂,经MZI作用的反馈光功率随τd的变化规律与方案三的情形不能完全对应。当τd处在某些特殊位置时(τd=0.01,0.02,0.106 ns),反馈强度很小,激光器处于稳态或者准周期振荡。
4 结论
针对双光反馈和双光注入两种方案,分析了MZI对抑制延时特征峰的规律。在MZI的两臂接近相等的区域,延时特征峰随τd进行一定的周期振荡,τd设定合适时,延时特征峰可以得到最佳的抑制效果。MZI的两臂明显不相等且不成比例时,延时特征峰得到一定的抑制,并且趋于一个稳定值,但是会出现新的延时特征峰。利用方案三分析了MZI的作用机理,当MZI的两臂接近相等时,两路混沌激光的电场强度和相位都具有较强的相关性,光的叠加为相干叠加,叠加得到的混沌激光的延时特征峰与τd的变化密切相关;当MZI的两臂明显不相等且不成比例时,两路混沌激光的电场强度和相位的相关性消失,光的叠加为非相干叠加,但是对延时特征峰的抑制还是主要由干涉项决定。通过探寻MZI在抑制时延特性中的规律和作用机理,为后续的研究提供一定的参考。
吴天安, 张晓旭, 张胜海, 兰淑静. 马赫-曾德尔干涉仪对抑制延时特征的作用机理研究[J]. 中国激光, 2018, 45(5): 0501003. Wu Tianan, Zhang Xiaoxu, Zhang Shenghai, Lan Shujing. Study on the Mechanism of Mach-Zehnder Interferometer for Suppressing the Delay Time Signature[J]. Chinese Journal of Lasers, 2018, 45(5): 0501003.