高斯型空间光孤子相互作用的数值研究

强激光与粒子束

HIGHPOWERLASERANDPARTICLEBEAMS

Vol.18,No.3

,Mar.2006

文章编号:()001-4322200603-0381-04 1

高斯型空间光孤子相互作用的数值研究

林晓东1, 吴正茂1, 夏光琼1, 陈建国2

()1.西南师范大学物理学院,重庆400715;10064 2.四川大学电子信息学院,成都6

*

对高斯孤子在对数型饱和非线性介质中的相互作用进行了研究,考查了 摘 要: 通过数值模拟的方法,两光束间的相对振幅和相对相位对其相互作用的影响。结果表明:高斯孤子之间的相互作用敏感地依赖于两光束间的主要作用交替地表现为相互排斥和相互光束间的相对振幅和相对相位。在不同的振幅差异范围内,

吸引,并由于高斯孤子的不稳定性,导致了光束在碰撞后以一种尺寸周期性变化的呼吸模式传输。随着相对相两光束间始终持续地表现出强烈的排斥作用,直到相对相位增加到一个2位的增大,π周期之后。而且碰撞之后,光束也都以呼吸模式进行传输,其分离的角度越大,呼吸就越明显。 关键词: 高斯孤子; 对数型饱和非线性介质; 分步傅里叶变换法437 文献标识码: 中图分类号: O A

在全光通信、数据处理以及光传输等方面有着巨大的应用前 空间光孤子由于其独特而优越的传输性能,

]1-4

景,自从人们首次发现光束的自陷现象以来,一直是研究的热点课题之一[。孤子之间的相互作用在很多方

5-7]面类似于粒子,因此在孤子的所有特性中,最吸引人的就是这种相互作用,或称为孤子间的相互碰撞[。孤子

间的相互作用就象粒子一样,可以发生弹性碰撞、非弹性碰撞以及完全非弹性碰撞等行为,但这些研究成果多集中于克尔介质、液晶等非线性介质中的孤子碰撞效应。

在克尔非线性介质中2维孤子是不稳定的,最终形成灾变自聚焦甚至光束分裂,而在具 前人的研究表明,

]8-10

有饱和非线性的介质里,空间孤子能够稳定地存在[。本文在分析非线性薛定谔方程(的基础上,利NLSE)

11]12-14]用分步傅里叶变换方法[,数值模拟了在一种对数型饱和非线性介质中的高斯型空间光孤子[的碰撞行

为,并研究了不同的入射参数对两个孤子碰撞行为的影响。研究表明,在这样的饱和非线性介质中,孤子碰撞行为比克尔介质中的弹性碰撞情况要复杂得多,会产生孤子融合、分裂以及湮灭等新现象。

1 数值模型

我们考虑比较简单的光束在对数饱和介质中(1+1)D的传输情况。假设空间光束以连续 为了便于分析,

波方式入射y方向导引的平板波导,这时非线性薛定谔方程就与y无关,因此可简化为

2

nkA1∂AΔ0∂iA=0+2+z2knx∂0∂0

()1

式中:波数k/式中第一项代表沿z轴传A为电场时空慢变包络;nc;n为光强引起的折射率变化。(1)ωΔ0=00第二项代表衍射效应,第三项代表非线性效应。输,

其非线性折射率变化部分为 对于对数型饱和非线性介质来说,

(/n=nlnIIΔ2th)

2

这里I=|AIn|,th为阈值光强,2为非线性折射系数。于是可以得到光束在对数型饱和介质中的传输方程

2

IA∂A∂ilnA=0+2+IZ2thX∂∂

()2

()()3

2

式中://Z=zNkX=xヘN,N=nkn0;200。

)具有高斯型分布的孤子解3 方程(

[12]

X[()Z)Z)X2]+i+β(α(A(X,Z)=B(Z)ex4p-2()

2σZ式中:/反映了位相弯曲的情况。当高斯B是振幅;e的宽度;σ为光束振幅1α为位相;β为光束波阵面的曲率,

{2

}2005-05-26;2005-11-02*收稿日期: 修订日期:

基金项目:西南师范大学校青年科技基金资助课题()220-413044

作者简介:林晓东(,男,西南师范大学讲师,主要从事非线性光学、半导体激光器与光纤通信研究;1975-)linxdwu.edu.cn@s

382

强激光与粒子束

第18卷

光束以束腰入射到对数型饱和介质中且腰斑半径σ=/高斯光束获得自陷,形成孤子。22时,σs=ヘ)式孤子解的相互作用。我们应用分步傅里叶变换方法对()式进行数值43 本文采用数值模拟的方法讨论(

A求解得到∂/对方程∂/求解。第一步只考虑前Δ2的衍射过程,ZAZ=i∂

2X2∂

1~2-/)=F{)[(/]}20ex24A1(ZA(Z-iΔπpfX)Δ

2

输入光,解方程∂/]计算整个Δ得到[(A1Z=-ilnA1IA1,Z的非线性效应,∂|/|th)2

)][(A2=A1exlnI-i|A1|/pth+1

/得到第三步再计算一次后Δ2的衍射过程,Z21~-{/]}[(4A(X,Z)=FA2ex2Z-iΔπpfX)Δ

()5

~()表示初始位置处光场的傅里叶变换,这里,0F-1表示逆傅里叶变换,AfX是光场的频谱。第二步以A1作为

()6

()7

只需把Z分解为若干个这就是在介质中传输了ΔZ距离之后的光场分布。要得到传输距离为Z的光场分布,

将()式计算的结果反复迭代就可以得到。只要Δ7Z,Z足够小,分步傅里叶变换法所得到的结果就是准确可Δ

]11

靠的[。

入射处光束中心与Z轴 我们假设两个高斯孤子分别从介质中心Z轴两侧以倾角α1和α2入射到介质中,

两光束的相对相位为θ,振幅分别为B1和B2,因此我们研究的初始入射高斯孤子对为分别相距为Xd1和Xd2,

22

((iiX+XdX0X-XdX01-1)2-2)(,)()[(,)()[)],](A1X0=Bxxexi A2X0=Bθp-ip-ip10e20e((ii2ZZ2ZZin1+01)in2+02)

()8

式中:瑞利长度Z/(其中j=1,参量,2,2)Zσσ000inj=jj为高斯光束的腰斑半径;j为光束入射点与束腰的距离;

/。(X0arctanX0Z0j与入射角αj的关系为αj=jj)

2 计算结果及分析

计算中采用的参数为I0.0001,Xd2, 图1表示不同振幅的高斯孤子相互碰撞的结果,σσth=1=Xd2=0=s/即两光束中心相距为4,以孤子条件平行入射到对数非线性介质中,且同相位。=ヘ22,θ=0,α=0,

当两孤子的振幅相同时,表现出周期性吸引和分离的现象,两光束就这样周而复始地相互 计算结果表明,

缠绕在一起(见图1());一旦两光束的振幅有了一点微小的差别,就会产生不同的碰撞效果。光束相互周期a性缠绕几次后最终还是分离开。

Fi.1 InfluenceofrelativeamlitudeontheGaussiansolitoninteractionsgp

图1 相对振幅对高斯孤子间相互作用的影响

随着非等幅比例的不断增大,两孤子间融合的次数逐渐减少(但发生融合 在振幅差别小于5%的范围内,

的位置不变),很快就完全分离开,两光束之间主要表现为相互排斥。当振幅差别达到5%时,两光束就不再发生融合,而是发生类似于弹性碰撞的现象(见图1())。在这之后,随着振幅差别的逐渐增大,两光束之间再次b表现出相互吸引的作用(见图1())。当振幅差别达到2两光束的吸引力和排斥力达到平衡,c5%左右的时候,

第3期

林晓东等:高斯型空间光孤子相互作用的数值研究

383

互不干扰地平行传输(见图1())。随着振幅差别进一步增大,两光束间的吸引力也增大,导致两光束在平行d传输一段距离之后,再次融合。当振幅差别达到8两光束就表现为相互穿越的现象了,随着振0%左右的时候,幅差别的增大,其穿越角度也越大,而且穿越之后的两孤子分别以呼吸模式进行传输(见图1())。这种现象e一直会持续到振幅差别为3这时由于两孤子之间的振幅差异太大,小幅孤子对大幅孤子几乎就没有00%左右,什么影响了(见图1())。我们用不同的振幅基数进行了研究,结果表明:对于同样的振幅差别百分比,两光束f的相互作用同样遵循以上规律。这正好说明了高斯孤子不同于克尔孤子的最大特点:高斯孤子的传输行为与其振幅无关,而只与初始位相和光束尺寸有关。

0)0.0001,B1( 图2表示两光束之间的相对相位对高斯孤子相互作用的影响。计算中采用的参数为Ith=)/即两光束中心相距为4,以孤子条件平行入射,且振幅相等。0=1,=B2(Xd2,22,σσα=0,1=Xd2=0=s=ヘFi.2 InfluenceofrelativephaseontheGaussiansolitoninteractionsg

图2 相对相位对高斯孤子间相互作用的影响

当两孤子相位相同(相对相位为零)时,两光束间同样表现出周期性缠绕现象(见图1())。a 计算结果表明,

在θ<π/两光束最初都有一点吸引,然后马但两光束的相互作用对相对相位十分敏感,一旦θ≠0,4的范围内,虽然两光束间会存在持续上就一直表现为相互排斥,永远分开(图2())。而在π/ad)4<θ<2π的范围内,～(

的吸引力,但它们不再发生融合,而永不相遇地独立传播下去(图2()和())。这表明,随着相对相位的逐渐增ef大,两光束间的排斥作用也越来越强,直到θ=2又回到了相对相位为零的状态。我们看到,在以孤子条件π时,平行入射且振幅相等相对相位不为零的情况下,两光束发生碰撞之后,各自也都以呼吸模式进行传输,而且其分离的角度越大,光束的呼吸就越明显。

可以定性解释为:我们把两个孤子间的相互作用理解为一个孤子对另一个 对于以上高斯孤子的碰撞行为,孤子的干扰。孤子本身是一个导引模,当它独自在介质中传输时,背景折射率是平坦的。当它附近出现了另一个孤子时,则该孤子的背景折射率将发生改变。一方面,当两孤子振幅相等时,它们对相互的背景折射率的影响是完全一样的,即这种作用是对称的,因此也就出现了孤子间的周期性缠绕现象(如图1())。当两孤子振a幅出现差别时,这种对称作用失去平衡,就出现了排斥以及穿越的现象(如图1())。当两孤子振幅相差b)e～(非常大时,小幅孤子对大幅孤子基本上没有什么影响。另一方面,共相时,两孤子交叠处折射率因干涉加强而)处高,两边低,因此两孤子相互吸引,孤子间增加,导致孤子的背景折射率发生倾斜,两孤子之间的中心(X=0的相互作用表现为吸引碰撞;失相时,两孤子交叠处折射率因干涉相消而减小,导致孤子的背景折射率在两孤子之间的中心处低,而两边高,因此两孤子相互排斥,孤子相互作用表现为排斥碰撞。

我们还可以看到一个值得注意的现象:和克尔介质有所不同,当两个高斯孤子相互融合、 从前面的分析中,

12]

,即光束尺寸以准周期变化的一种空间孤子形态。这穿越之后,其各自的传输行为就表现为一种呼吸模式[是由于高斯孤子不象克尔孤子那样稳定,一旦受到碰撞,高斯孤子将具有一个偏离孤子条件的微小变化,孤子

15]

就象处于势阱底部的粒子受到扰动而在势阱中来回周期性运动一样,变成了所谓的呼吸模式[。

384

强激光与粒子束

第18卷

3 结 论

高斯孤子之间的相 我们数值模拟了高斯型空间孤子在对数型饱和非线性介质中的相互作用。结果表明,

微小的振幅差异会导致高斯孤子互作用敏感地依赖于两光束间的相对振幅和相对相位。在同相位的情况下,

间不同的相互作用,在不同的振幅差异范围内,光束间的主要作用交替地表现为相互排斥和相互吸引,并由于高斯孤子的不稳定性,光速在碰撞后以一种呼吸模式传输。在同振幅的情况下,其相对相位的变化同样会影响高斯孤子间的碰撞行为,但与相对振幅的影响不同的是,随着相对相位的增大,两光束间持续地表现出强烈的排斥作用,直到相对相位增加到一个2再重复上一个周期的行为。总之,高斯孤子间的碰撞行为比π周期之后,克尔孤子等的碰撞行为要丰富得多,因此,如何有效地利用高斯孤子间这些独特的碰撞行为,是值得我们深入研究的课题。

致 谢:在此表示衷心感谢！ 本论文的工作得到唐永林博士的热心指导,

参考文献:

[][]():1hiaoRY,GarmireE,TownesCH.Self-trainfoticalbeamsJ.PhsRevLett,1964,1315479. Cppgopy[][]2arthelemaneufS,FrochlProaationsolitonetauto-confinementdefaisceauxlaserparnonlineariteotiuedekerrJ.Ot ByA,MyC.pgpqp():1985,553201.Commun,

[][]3ureeGC,ShultzJL,Shar.Observationofself-trainfanoticalbeamduetophotorefractiveeffectJ.PhsRevLett,1993,71 DpGJppgopy():4533.

[]]():4heWL,LeeKK,LeeWK.Observationoftwo-dimensionalbrihthotovoltaicsatialsolitons[J.PhsRevLett,1999,83163182. Sgppy[],:[]5teemanGISeevM.OticalsatialsolitonsandtheirinteractionsuniversalitnddiversitJ.Science,1999,286:1518. Sggppyay

[][]:6ecciantiM,BrzdakiewiczKA,AssantoG.NonlocalsatialsolitoninteractionsinnematicliuidcrstalsJ.OtLett,2002,27(16) Ppqyp1460.

[],[]():7ltanirEA,SteemanGILaneCH.InteractionsofdissiativesatialsolitonsJ.OtLett,2004,293283. Uggppp[]8atzS,HerrmannJ.Proaationofoticalbeamsandtheproertiesoftwo-dimensionalsatialsolitonsinmediawithalocalsaturablenon- Gpgppp

[]():J.JOSAB,1997,1471795.linearrefractiveindex

[][]():9idalF,JohnstonTW.Asmtoticbehaviorofradiallmmetricself-focusedbeamsJ.PhsRevE,1997,5533571. Vypysyy[][]():10hristodoulidesDN,CoskunTH,JosehRI.IncoherentsatialsolitonsinsaturablenonlinearmediaJ.OtLett,1997,22141080. Cppp[]:11ahaTR,AblowitzMJ.AnalticalandnumericalasectsofcertainnonlinearevolutioneuationsII.NumericalnonlinearSchrödiner Typqg

[]euationJ.Joomhs,1984,55:203.qfCpPy[][]():12nderAW,MitchellJD.MihtorhinatialsolitonsandbulletsJ.OtLett,1997,22116. Sygympgspp[][]13anChenJG,LiDY,etal.SecificGaussianbeamsproaatininsideloarithmicallaturablenonlinearmediaJ.JModOt, TgYL,ppgggysp():1999,4681177.

[][]():14ikhonenkoV.ObservationofmihtorhinatialsolitonsJ.OtLett,1998,238594. Tgympgspp[]]():(15J.物理学报,1999,4871248.TanLi 唐永林,李大义,陈建国,等.对数型折射率饱和非线性介质中光孤子高斯型呼吸模式[gYL,

DY,ChenJG,etal.Gaussian-tebreathmodesofsatialsolitonformedinloarithmicallaturablenonlinearmedia.AcataPhsicayppgysy():)Sinica,1999,4871248

NumericalresearchofinteractionsbetweenGaussian-shaedsatialsolitonspp

1112

,LINXiao-donhen-maoIAGuan-ionNJian-uo WUZ X CHEg,ggqg,g

(1.Schoolohsics,SouthwestNormalUniversitChonin00715,China;fPyy,gqg42.Deartmentotoelectronics,SichuanUniversitChendu610064,China)pfOpy,g:bstractheinteractionsbetweenGaussian-shaedsolitonsintheloarithmicsaturablenonlinearmediawerenumericall A Tpgyinvestiated.TheresultsshowthattheinteractionsbetweenGaussian-shaedsolitonsaresensitivetorelativeamlitudeandrela-gpp,,tivephase.Indifferentreionsofrelativeamlitudethemaineffectisreulsionorattractionalternatelandduetotheinstabili-gppy,,tfGaussiansolitonsthebeamsbreatheperiodicallfterinteractions.Aloniththerisinftherelativephasethereareyoyagwgo,,stronerreulsionsbetweenthebeamsuntilattheendofeachccle.Thebeamsbreathealsointhissituationandthelareran-gpyg,letheearatethemorestronlhereathe.gyspgytyb

:;;eordsaussiansolitonoarithmicsaturablenonlinearmedialit-steouriertransform K G L SgppFyw