級聯(lián)單模光纖中初始啁啾對高階孤子脈沖壓縮的影響

徐永釗， 張 耿， 葉 海， 劉敏霞

(東莞理工學(xué)院 電子工程學(xué)院， 廣東 東莞 523808)

?

級聯(lián)單模光纖中初始啁啾對高階孤子脈沖壓縮的影響

徐永釗， 張 耿， 葉 海， 劉敏霞*

(東莞理工學(xué)院 電子工程學(xué)院， 廣東 東莞 523808)

采用兩段級聯(lián)單模光纖對高階孤子脈沖進(jìn)行壓縮。兩段光纖具有不同的反常色散值，當(dāng)高階孤子脈沖在第一段光纖中獲得最大程度壓縮時(shí)，通過轉(zhuǎn)換色散值不同的光纖，使壓縮脈沖繼續(xù)以高階孤子的形式在第二段光纖中再次被壓縮。每段光纖的長度都進(jìn)行了優(yōu)化，使得脈沖在每段光纖中都獲得最大程度的壓縮。基于非線性薛定諤方程，數(shù)值研究了初始啁啾對高階孤子脈沖壓縮的影響。研究結(jié)果表明，初始啁啾對高階孤子脈沖的壓縮有重要影響。與無初始啁啾時(shí)的情形相比，正的初始啁啾能增強(qiáng)每段光纖中脈沖的壓縮效果，降低壓縮脈沖的基座能量，而負(fù)初始啁啾的影響則相反。隨著初始啁啾參量Cp的增大，脈沖在每段光纖中的壓縮因子均增加，而基座能量、最優(yōu)光纖長度均減小。

非線性光學(xué)； 光脈沖壓縮； 光孤子； 啁啾

1 引 言

超短光脈沖在高速光通信、超快光子學(xué)、非線性物理以及生物醫(yī)學(xué)等眾多領(lǐng)域有重要的應(yīng)用，脈沖壓縮技術(shù)的研究一直受到人們廣泛的關(guān)注[1-7]。利用光纖壓縮光脈沖是獲得超短光脈沖的主要手段之一，常用的方案有兩種：絕熱孤子壓縮和高階孤子壓縮[8]。絕熱孤子壓縮技術(shù)通常使用色散漸減光纖作為壓縮光纖，可以獲得基座能量很低的壓縮脈沖[9]，但是這種壓縮方法需要較長的光纖長度。另外，由于色散漸減光纖色散漸減率的限制，脈沖最大壓縮因子的典型值限制在20左右[10]。而利用高階孤子壓縮技術(shù)可以獲得很高的脈沖壓縮因子，但壓縮過程中有很多能量轉(zhuǎn)移到基座中，導(dǎo)致壓縮脈沖具有寬且高的基座[11]。為了避免這兩種壓縮技術(shù)的顯著缺點(diǎn)，Li等[12-13]提出了一種基于級聯(lián)單模光纖的高階孤子脈沖壓縮技術(shù)。他們的方案采用兩段或三段具有不同色散的單模光纖級聯(lián)在一起，對高階孤子脈沖進(jìn)行壓縮。每段光纖的長度都進(jìn)行了優(yōu)化，以使脈沖在每段光纖內(nèi)都獲得最大程度的壓縮。研究結(jié)果顯示，每段光纖的最優(yōu)長度都小于半個(gè)孤子周期，脈沖經(jīng)過兩個(gè)或三個(gè)階段的級聯(lián)壓縮，可以獲得很高程度的壓縮，且基座相對較小。例如，兩個(gè)階段的三階孤子壓縮可以獲得70.8的壓縮因子，基座能量僅為44.8%；三個(gè)階段的三階孤子壓縮則可以獲得高達(dá)599.7的壓縮因子，而基座能量也僅為58.8%?；诩壜?lián)單模光纖的高階孤子脈沖壓縮技術(shù)優(yōu)于常規(guī)的絕熱孤子壓縮和高階孤子壓縮技術(shù)，但是，Li等的研究僅考慮了對無啁啾高階孤子脈沖壓縮的情形。據(jù)我們所知，脈沖的初始啁啾對級聯(lián)高階孤子脈沖壓縮的影響目前還沒有這方面的研究報(bào)道。眾所周知，實(shí)際的激光器輸出的脈沖經(jīng)常具有啁啾，作為脈沖的一個(gè)重要參量，研究啁啾對脈沖壓縮效應(yīng)的影響具有現(xiàn)實(shí)意義。

本文通過數(shù)值求解非線性薛定諤方程，研究了初始啁啾對級聯(lián)單模光纖中高階孤子脈沖壓縮的影響。結(jié)果發(fā)現(xiàn)，初始啁啾對級聯(lián)單模光纖中高階孤子脈沖的壓縮有重要影響。

2 理論模型

光脈沖在單模光纖中的傳輸采用非線性薛定諤方程描述，表示如下[14]：

(1)

其中A為光場振幅，z為傳輸距離，α為光纖損耗系數(shù)，τ為以群速度為移動(dòng)參考系的時(shí)間參量，βm為在中心頻率ω0處m階色散系數(shù)，γ為光纖的非線性系數(shù)。方程左邊第二項(xiàng)為損耗，第三項(xiàng)為色散效應(yīng)；方程右邊第一項(xiàng)為自相位調(diào)制(SPM)，后面兩項(xiàng)為高階非線性效應(yīng)項(xiàng)，分別為自變陡和脈沖內(nèi)拉曼散射。

假設(shè)初始脈沖為雙曲正割脈沖，脈沖的光場表達(dá)為：

(2)

其中，P0是脈沖的峰值功率；T0為脈沖1/e強(qiáng)度處的半寬度，對于雙曲正割脈沖，T0與脈沖的半高全寬(TFWHM)的關(guān)系為TFWHM≈1.763T0；Cp為脈沖的初始線性啁啾參量。

本文研究高階孤子脈沖在單模光纖中的壓縮，脈沖的孤子階數(shù)N定義為：

(3)

其中LD和LNL分別為色散長度和非線性長度。對于高階孤子脈沖，N>1。高階孤子脈沖在光纖中傳輸時(shí)，其波形會發(fā)生周期性的演化，孤子周期為：

(4)

基于級聯(lián)光纖的高階孤子脈沖壓縮技術(shù)的原理是當(dāng)高階孤子在一段光纖中獲得最大程度的壓縮時(shí)，立刻輸入到下一段光纖，通過轉(zhuǎn)換光纖的色散，使輸入脈沖在接下來的光纖中重新作為高階孤子繼續(xù)被壓縮。在級聯(lián)壓縮的過程中，每段光纖都選取最優(yōu)的光纖長度和適當(dāng)?shù)纳⒅?，使脈沖作為高階孤子在每個(gè)階段都達(dá)到最大程度的壓縮。經(jīng)過兩段或三段級聯(lián)光纖，即通過兩個(gè)或三個(gè)階段的級聯(lián)壓縮，脈沖最終可以獲得很高的壓縮因子。為了衡量脈沖壓縮后的質(zhì)量，引入壓縮因子FC和基座能量(%)兩個(gè)參量進(jìn)行描述，其中壓縮因子FC定義為初始脈沖的寬度與壓縮后脈沖的寬度之比：

(5)

壓縮脈沖的基座能量定義為：

(6)

其中，Etotal是壓縮脈沖的總能量，Esech是理想雙曲正割脈沖的能量，該雙曲正割脈沖具有與壓縮脈沖相同的寬度(TFWHM)和峰值強(qiáng)度。認(rèn)為壓縮脈沖有基座，雙曲正割脈沖無基座，則兩者能量之差相對于壓縮脈沖的總能量即為壓縮脈沖的基座能量，代表基座能量所占壓縮脈沖總能量的百分比。

3 結(jié)果與討論

本文考慮基于兩段級聯(lián)光纖的高階孤子壓縮的情形。設(shè)第一段光纖初始輸入脈沖的孤子階數(shù)為N1，當(dāng)脈沖在第一段光纖中壓縮到最大程度時(shí)從光纖中輸出，并輸入到第二段光纖中。設(shè)第二段光纖輸入脈沖的孤子階數(shù)為N2。假設(shè)N1=N2=N，且兩段光纖的非線性系數(shù)相同(γ=2.5 W/km)，其他參量的取值分別為：初始脈沖的寬度TFWHM=30 ps，第一段光纖的二階色散系數(shù)β21=-20 ps2/km(反常色散)。在數(shù)值計(jì)算中，忽略光纖衰減的影響。另外，為便于與文獻(xiàn)[12]的結(jié)果進(jìn)行比較，我們參考其做法，忽略了高階色散和高階非線性效應(yīng)的影響。實(shí)際上，對于T0>5 ps的脈沖，高階色散和高階非線性效應(yīng)對脈沖傳輸?shù)挠绊懞苄　?/p>

圖1 無初始啁啾的三階孤子(N=3)脈沖的峰值功率在一個(gè)孤子周期內(nèi)的演化圖。(a)第一段光纖中一個(gè)孤子周期z01內(nèi)的演化圖；(b)第二段光纖中一個(gè)孤子周期z02內(nèi)的演化圖。

Fig.1 Peak power evolution of the third order soliton (N=3) pulse within one soliton period. (a) Peak power evolution in the first fiber within one soliton periodz01. (b) Peak power evolution in the second fiber within one soliton periodz02.

首先考慮無初始啁啾的三階孤子(N=3)脈沖在光纖中壓縮的情形。圖1(a)所示為第一段光纖內(nèi)，脈沖的峰值功率在一個(gè)孤子周期z01內(nèi)的演化圖。圖中虛線所標(biāo)記的位置為脈沖峰值功率達(dá)到最大值時(shí)的相對位置L1/z01，在該位置脈沖獲得最大程度的壓縮，相應(yīng)的光纖長度L1為第一段光纖的最優(yōu)長度。脈沖在L1處輸出后，隨即作為第二段光纖的輸入脈沖。若第二段光纖的二階色散系數(shù)β22取值為-1.769 ps2/km時(shí)，可使得輸入脈沖的孤子階數(shù)仍然為N=3。圖1(b)所示為第二段光纖內(nèi)，脈沖的峰值功率在一個(gè)孤子周期z02內(nèi)的演化圖?？梢钥吹?，在一個(gè)孤子周期內(nèi)，脈沖的峰值功率在兩段光纖中演化的情形極為相似，而且L1/z01≈L2/z02。與圖1(a)類似，脈沖在相對位置L2/z02處獲得最大程度的壓縮，第二段

光纖的最優(yōu)長度為L2。表1給出了無初始啁啾的三階孤子壓縮時(shí)的光纖設(shè)計(jì)。脈沖經(jīng)過第一段光纖的壓縮，獲得的壓縮因子為8.4，基座能量為26.0%(基座能量所占脈沖總能量的百分比)；經(jīng)過第二段光纖的壓縮，最終獲得的壓縮因子為70.8，基座能量為44.8%。

表1 無初始啁啾(Cp=0)的三階孤子(N=3)脈沖壓縮時(shí)的光纖設(shè)計(jì)

Tab.1 Fiber design in pulse compression forN=3 andCp=0

β2i/(ps2·km-1)Li/kmLi/Z0i第一段光纖-205.3910.237第二段光纖-1.7690.8650.235

接下來，我們改變脈沖的初始啁啾參量Cp，研究初始啁啾對三階孤子脈沖壓縮的影響。圖2(a)和(b)分別給出了脈沖經(jīng)過第一段光纖和第二段光纖壓縮后，壓縮因子和基座能量隨初始啁啾參量Cp的變化曲線。由圖可見，隨著啁啾參量Cp的增大，壓縮因子單調(diào)增加，而基座能量則單調(diào)減小。在兩段光纖中，壓縮因子和基座能量隨初始啁啾參量Cp的變化規(guī)律相同。與無初始啁啾(Cp=0)時(shí)的情形相比，可以發(fā)現(xiàn)，正初始啁啾(Cp>0)可以增強(qiáng)脈沖的壓縮效果，減小壓縮脈沖的基座能量。當(dāng)正啁啾較大時(shí)，這種效果非常顯著。例如，與Cp=0時(shí)的結(jié)果相比，當(dāng)Cp=4時(shí)，脈沖經(jīng)過第一段光纖的壓縮，壓縮因子由8.4增加到12.0，而基座能量由26.0%減少到僅9.0%；經(jīng)過第二段光纖的壓縮，壓縮因子由70.8增加到102.0，而基座能量則由44.8%減少到31.8%。相反，負(fù)初始啁啾(Cp<0)會使脈沖的壓縮效果減弱，基座能量增加。負(fù)初始啁啾值越大(負(fù)初始啁啾參量的絕對值Cp越大)，脈沖壓縮的質(zhì)量越差。

為比較壓縮的效果，圖3(a)和(b)分別給出了三階孤子在第一段和第二段光纖中，當(dāng)Cp=0和Cp=4時(shí)，脈沖壓縮程度最大時(shí)的波形。通過圖3中脈沖波形的對比，可以清晰地看出正初始啁啾對脈沖壓縮效果的改善。

圖2 初始啁啾對壓縮因子和基座能量的影響。(a)脈沖在第一段光纖的壓縮因子和基座能量隨初始啁啾參量Cp的變化曲線；(b)脈沖在第二段光纖的壓縮因子和基座能量隨初始啁啾參量Cp的變化曲線。

Fig.2 Effect of initial chirp on compression factor and pedestal energy. (a) Calculated compression factor and pedestal energy as a function ofCpin the first fiber. (b) Calculated compression factor and pedestal energy as a function ofCpin the second fiber.

圖3 當(dāng)Cp=0和Cp=4時(shí)，三階孤子脈沖在第一段和第二段光纖中壓縮程度最大時(shí)的波形。(a)脈沖在第一段光纖中壓縮程度最大時(shí)的波形；(b)脈沖在第二段光纖中壓縮程度最大時(shí)的波形。

Fig.3 Pulse shapes where compression is maximized in both the first and second fibers of the third-order soliton whenCp=0 andCp=4. (a) Pulse shapes where compression is maximized in the first fiber. (b) Pulse shapes where compression is maximized in the second fiber.

圖4 第二段光纖的色散值與第一段光纖的色散值的比值β22/β21隨初始啁啾參量Cp的變化曲線

Fig.4β22/β21as a function ofCp

圖5 第一段光纖(a)和第二段光纖(b)的最優(yōu)光纖長度和最優(yōu)光纖長度與孤子周期的比值隨初始啁啾參量Cp的變化曲線

Fig.5 Optimal length and the ratio of optimal length to soliton period of the first fiber (a) and the second fiber (b) as a function ofCp

為便于比較，圖5中同時(shí)也給出了兩段光纖的最優(yōu)長度與孤子周期的比值Li/z0i(i=1,2)隨初始啁啾參量Cp的變化曲線。在第一段光纖中，L1/z01和L1隨啁啾參量Cp的變化趨勢相同；但在第二段光纖中，L2/z02基本保持不變，與初始啁啾參量Cp無關(guān)。L2/z02僅與輸入第二段光纖的孤子階數(shù)N有關(guān)。

保持初始脈沖的孤子階數(shù)不變，我們計(jì)算了不同寬度脈沖的壓縮效果。結(jié)果表明，在不考慮損耗、高階色散和高階非線性效應(yīng)的情況下，當(dāng)初始啁啾相同時(shí)，對于不同寬度的初始脈沖，其壓縮因子、基座能量的占比以及兩段光纖最優(yōu)長度與孤子周期的比值Li/z0i均保持不變。

表2 初始啁啾參量Cp=0和Cp=4時(shí)，不同孤子階數(shù)的脈沖在兩段級聯(lián)光纖中的壓縮因子和基座能量

我們同時(shí)也研究了其他階數(shù)的高階孤子脈沖壓縮的情形。結(jié)果發(fā)現(xiàn)，初始啁啾對不同階數(shù)的高階孤子的影響有相同的規(guī)律。表2給出了初始啁啾參量Cp=0和Cp=4時(shí)，不同孤子階數(shù)的脈沖在兩段級聯(lián)光纖中的壓縮因子和基座能量的對比，其中初始啁啾參量Cp=0時(shí)的結(jié)果即文獻(xiàn)[12]的計(jì)算結(jié)果。表2的結(jié)果顯示，與無初始啁啾時(shí)的情形相比，當(dāng)Cp=4時(shí)，不同N值的壓縮因子都有顯著的提高，而基座能量都明顯降低。另外，N值越小，初始啁啾的影響越顯著。例如，當(dāng)Cp=4時(shí)，對于二階孤子(N=2)，與無初始啁啾時(shí)的情況相比，經(jīng)過第一級壓縮后，壓縮因子增大了一倍，而基座能量減少67%；經(jīng)過第二級壓縮后，壓縮因子也接近增大一倍，而基座能量減少48%。隨著孤子階數(shù)的增加，壓縮因子的增幅減小，基座能量下降的幅度也減小。

表3給出了Cp=4時(shí)，不同N值的高階孤子的兩段級聯(lián)光纖的設(shè)計(jì)參數(shù)。結(jié)果顯示，在相同初始啁啾的情況下，隨著N的增加，兩段光纖色散值的比值β22/β21減小。另外，隨N值的增加，L1/z01和L2/z02同時(shí)減小。這表明脈沖的孤子階數(shù)越高，在一個(gè)孤子周期內(nèi)發(fā)生最大程度壓縮的相對距離越短。例如，當(dāng)N=2時(shí)，第一段光纖和第二段光纖分別在L1=0.146z01處和L2=0.475z02處發(fā)生最大程度的壓縮；當(dāng)N=5時(shí)，第一段光纖和第二段光纖分別在L1=0.0735z01處和L2=0.109z02處發(fā)生最大程度的壓縮。當(dāng)Cp=4時(shí)，對于不同的N值，兩段光纖的最優(yōu)長度L1和L2均小于半個(gè)孤子周期。

表3 Cp=4時(shí)，兩段級聯(lián)光纖的設(shè)計(jì)

4 結(jié) 論

本文數(shù)值研究了初始啁啾對兩段級聯(lián)單模光纖中高階孤子脈沖壓縮的影響。研究結(jié)果表明，與無初始啁啾時(shí)的情形相比，正的初始啁啾能增強(qiáng)每段光纖中脈沖的壓縮效果，降低壓縮脈沖的基座能量，而負(fù)初始啁啾的影響則相反。隨著初始啁啾參量Cp的增大，脈沖在每段光纖中的壓縮因子均增加，而基座能量均減小。當(dāng)正初始啁啾較大時(shí)，與無初始啁啾的情形相比，脈沖壓縮的質(zhì)量可以獲得十分顯著的改善。另外，隨著初始啁啾參量Cp的增大，兩段光纖的最優(yōu)長度以及第一段光纖的最優(yōu)長度與孤子周期的比值L1/z01均減小，而第二段光纖的最優(yōu)長度與孤子周期的比值L2/z02則保持不變，與初始啁啾無關(guān)。對于不同階數(shù)的高階孤子，初始啁啾對脈沖壓縮的影響也不同。當(dāng)脈沖具有相同正初始啁啾時(shí)，隨著孤子階數(shù)的增加，與無初始啁啾的情形相比，壓縮因子的增幅減小，基座能量下降的幅度也減小。

[1] 曹文華,徐平,劉頌豪. 馬赫-曾德爾型色散漸減光纖干涉儀的孤子效應(yīng)脈沖壓縮研究 [J]. 光學(xué)學(xué)報(bào), 2011, 31(4):419001-1-8. CAO W H, XU P, LIU S H. Soliton-effect pulse compression in a dispersion-decreasing fiber-based Mach-Zehnder interferometer [J].ActaOpt.Sinica, 2011, 31(4):419001-1-8. (in Chinese)

[2] 李博,婁淑琴,譚中偉,等. 兩種基于交叉相位調(diào)制的光脈沖壓縮方案 [J]. 物理學(xué)報(bào), 2012, 61(19):194203-1-7. LI B, LOU S Q, TAN Z W,etal.. Two kinds of optical pulse compression approaches based on cross phase modulation [J].ActaPhys.Sinica, 2012, 61(19):194203-1-7. (in Chinese)

[3] WANG D, LENG Y X, XU Z Z. Optical pulse compression of ultrashort laser pulses in a multi-hollow-core fiber [J].Opt.Commun., 2012, 285(9):2418-2421.

[4] VORONIN A A, MIKHAILOVA J M, GORJAN M,etal.. Pulse compression to sub-cycle field waveforms with split-dispersion cascaded hollow fibers [J].Opt.Lett., 2013, 38(21):4354-4357.

[5] LEFORT C, KALASHYAN M, DUCOURTHIAL G,etal.. Sub-30-fs pulse compression and pulse shaping at the output of a 2-m-long optical fiber in the near-infrared range [J].J.Opt.Soc.Am. B, 2014, 31(10):2317-2324.

[6] BANANEJ A, AGHAEE A K, MADJABADI A,etal.. Adjustable enhancement of higher order soliton effect pulse compression by using optimized active dispersion decreased nonlinear fiber loop mirror [J].Optik-Int.J.LightElectronOpt., 2014, 125(19):5476-5479.

[7] RUBENCHIK A M, CHEKHOVSKOY I S, FEDORUK M P,etal.. Nonlinear pulse combining and pulse compression in multi-core fibers [J].Opt.Lett., 2015, 40(5):721-724.

[8] AGRAWAL G P.ApplicationofNonlinearFiberOptics[M]. 2nd ed. New York: Academic Press, 2008.

[9] MOSTOFI A, HATAMI-HANZA H, CHU P L. Optimum dispersion profile for compression of fundamental solitons in dispersion decreasing fibers [J].IEEEJ.Quant.Electron., 1997, 33(4):620-628.

[10] PELUSI M D, LIU H F. Higher order soliton pulse compression in dispersion-decreasing optical fibers [J].IEEEJ.Quant.Electron., 1997, 33(8):1430-1439.

[11] CHAN K, LIU H F. Short pulse generation by higher order soliton-effect compression: effects of optical fiber characteristics [J].IEEEJ.Quant.Electron., 1996, 31(12):2226-2235.

[12] LI Q, KUTZ J N, WAI P K A. Cascaded higher-order soliton for non-adiabatic pulse compression [J].J.Opt.Soc.Am. B, 2010, 27(11):2180-2189.

[13] LI Q, KUTZ J N, WAI P K A. High-degree pulse compression and high-coherence supercontinuum generation in a convex dispersion profile [J].Opt.Commun., 2013, 301-302:29-33.

[14] AGRAWAL G P.NonlinearFiberOptics[M]. San Diego: Academic Press, 2007.

徐永釗(1972-)，男，廣東廣州人，博士，副教授，2007年于北京郵電大學(xué)獲得博士學(xué)位，主要從事光纖通信技術(shù)、非線性光纖光學(xué)的研究。

E-mail： yongzhaoxu@126.com劉敏霞(1979-)，女，河北石家莊人，博士，副教授，2007年于北京大學(xué)獲得博士學(xué)位，主要從事光電材料、非線性光學(xué)的研究。

E-mail： liumxdgut@qq.com

Effect of Initial Frequency Chirp on Pulse Compression of Higher-order Solitons in Cascaded Single-mode Fibers

XU Yong-zhao， ZHANG Geng， YE Hai， LIU Min-xia*

(SchoolofElectronicEngineering,DongguanUniversityofTechnology,Dongguan523808,China)

We consider the pulse compression of higher-order solitons in two cascaded single-mode fibers. The two single-mode fibers have different anomalous dispersion value. When the higher-order soliton pulse achieves maximal compression in the first fiber, by switching the fiber with different dispersion, the pulse is compressed again as a new higher-order soliton in the second fiber. Each fiber length is optimized in order to achieve maximal compression inside each fiber segment. Based on the generalized nonlinear Schr?dinger equation, we numerically study the effect of initial frequency chirp on pulse compression of higher-order solitons. Our numerical results show that initial frequency chirp has significant influence on higher-order solitons pulse compression. In comparison with the case of without initial chirp, a positive chirp can enhance the pulse compression and can reduce the pedestal energy in each fiber segment, while a negative chirp does the opposite. In each fiber segment, as the chirp parameterCpincreases, the compression factor increases, while pedestal energy and optimal fiber length decrease.

nonlinear optics; optical pulse compression； optical solitons； chirp

1000-7032(2016)11-1360-07

2016-05-07；

2016-06-20

國家自然科學(xué)基金(61501118)； 廣東省自然科學(xué)基金(2014A030310262)資助項(xiàng)目

TN929．11

A

10.3788/fgxb20163711.1360

*CorrespondingAuthor,E-mail:liumxdgut@qq.com