吳 瓊， 任志君， 李曉東， 杜林岳， 毛和法， 胡海華

(1.浙江師范大學(xué) 浙江省光信息檢測(cè)與顯示技術(shù)研究重點(diǎn)實(shí)驗(yàn)室,浙江 金華 321004;2.浙江師范大學(xué) 信息光學(xué)研究所 浙江 金華 321004)

0 引 言

1987年，美國(guó)羅徹斯特大學(xué)的Durnin等[1-2]首次提出了具有實(shí)用價(jià)值的“無衍射光束”的概念，產(chǎn)生的第一種無衍射貝塞爾光束是一種中心光斑直徑很小且不隨傳播距離而改變的光束.隨后，Gutiérrez-Vega的研究小組[3]指出，波動(dòng)方程的無衍射解有4種形式，分別對(duì)應(yīng)直角坐標(biāo)系下的余弦光束(平行光模型)、極坐標(biāo)系下的貝塞爾光束、拋物坐標(biāo)下的拋物線光束及橢圓坐標(biāo)系下的馬丟光束.

相對(duì)于直角坐標(biāo)系下的余弦光束模型和極坐標(biāo)系下的貝塞爾光束模型，橢圓坐標(biāo)系下的馬丟光束模型在數(shù)學(xué)上更為復(fù)雜，光學(xué)形態(tài)也豐富得多.因此,馬丟光束在實(shí)際中有新的用途，比如,馬丟光束可以用于光學(xué)格子的產(chǎn)生等[4-5].馬丟光束包括第一、第二類奇型和第一、第二類偶型分布.與余弦函數(shù)和貝塞爾函數(shù)相同，4種類型馬丟函數(shù)也不能平方可積，因此,理想馬丟光束只能被看作是一個(gè)理論模型.實(shí)驗(yàn)上能產(chǎn)生的馬丟光束常需要用高斯函數(shù)剪切得到，這就是馬丟-高斯光束[6-7]，它具有準(zhǔn)無衍射傳輸特性.因此，實(shí)際中人們更多關(guān)注的是馬丟-高斯光束的產(chǎn)生和傳輸特性.

衍射問題是光學(xué)中極難處理的問題之一，極少的衍射問題能求得嚴(yán)格解.為研究光束的傳輸，在非傍軸近似理論中,虛源法是由Deschamps[8]提出的，被Felsen的研究小組[9]進(jìn)行系統(tǒng)的發(fā)展，用來精確求解亥姆霍茲方程，進(jìn)而精確研究光束在自由空間非傍軸傳輸?shù)囊环N方法.截止目前，虛源法[9]已被用于研究幾種具有實(shí)用價(jià)值光束的傳輸特性，比如余弦-高斯光束[10]、貝塞爾-高斯光束[11]、加速艾里光束[12]、厄米-拉蓋爾-高斯光束[13]等，都得到了很有價(jià)值的研究結(jié)果.然而，對(duì)于橢圓坐標(biāo)系下的馬丟函數(shù)，由于其數(shù)學(xué)上的復(fù)雜性，至今還沒有研究者能夠運(yùn)用虛源法來研究馬丟-高斯光束的傳輸特性.基于光波傳播的獨(dú)立性原理和疊加性原理，利用數(shù)學(xué)上的級(jí)數(shù)展開，復(fù)雜光束可展開為一系列簡(jiǎn)單光束的疊加.據(jù)此，筆者以第一類奇型馬丟-高斯光束為例，將級(jí)數(shù)展開與虛源法相結(jié)合，推導(dǎo)得到了這種光束傳輸?shù)姆e分表達(dá)式.利用類似的處理方法，還分別得到了第二類奇型馬丟-高斯光束和第一、第二類偶型馬丟-高斯光束的積分表達(dá)式.這一組積分表達(dá)式為精確認(rèn)識(shí)馬丟-高斯光束的傳輸特性提供了方便.

1 用虛源法產(chǎn)生奇型馬丟高斯光束

馬丟光束有4種類型.在橢圓坐標(biāo)系下，根據(jù)馬丟函數(shù)的定義，馬丟光束的4種類型分別為：第一類2m+2階奇型馬丟光束se2m+2(η,q)Jo2m+2(ξ,q)，第二類2m+1階奇型馬丟光束se2m+1(η,q)Jo2m+1(ξ,q)；第一類2m階偶型馬丟光束ce2m(η,q)Je2m(ξ,q)，第二類2m+1階偶型馬丟光束ce2m+1(η,q)Je2m+1(ξ,q).其中，m=1，2，3,….表達(dá)式中，ce和Je分別代表偶型角向和徑向的馬丟函數(shù)，se和Jo分別代表奇型角向和徑向的馬丟函數(shù)[6-7,14-18].

橢圓坐標(biāo)系在數(shù)學(xué)上較難理解.既然光波在自由空間中可獨(dú)立傳播和疊加，那么筆者就采用與文獻(xiàn)[6-7]相同的數(shù)學(xué)手段，利用數(shù)學(xué)級(jí)數(shù)展開的方式，將橢圓坐標(biāo)系下的馬丟函數(shù)展開為極坐標(biāo)系下貝塞爾函數(shù)的形式來處理馬丟-高斯光束.兩類奇型馬丟函數(shù)可展開為:

(1)

(2)

2類偶型馬丟函數(shù)可展開為:

(3)

(4)

本文以第一類2m+2階奇型馬丟-高斯光束為例進(jìn)行研究.為產(chǎn)生馬丟-高斯光束，設(shè)標(biāo)量波函數(shù)為E(r,θ,z).光束的傳輸需滿足波動(dòng)方程沿z方向傳播的馬丟-高斯光束，在z=0處，有

(5)

式(5)可以看作第一類奇型馬丟-高斯光束傳輸?shù)倪吔鐥l件.

在圓柱坐標(biāo)系下，假定有無窮多個(gè)位于z=zex且強(qiáng)度為Sex(2l+2)的電環(huán)，環(huán)半徑為l=lex.在z>0空間產(chǎn)生第一類奇型馬丟-高斯光束，該光源的標(biāo)量波函數(shù)可寫為

E2m+2,2l+2(r,θ,z)=M2m+2,2l+2(r,z)sin[(2l+2)θ].

(6)

在式(5)所示的邊界條件下，式(6)滿足非齊次亥姆霍茲方程

(7)

根據(jù)極坐標(biāo)系下的傅里葉-貝塞爾變換

(8)

(9)

(10)

(11)

求解式(11)可得

(12)

在z>0處，為了獲得馬丟高斯光束，必須滿足式(5)給出的z=0處的邊界條件.即在傍軸近似下，式(12)需滿足如下的邊界條件:

(13)

由式(13)給定的邊界條件可得到式(12)中的參量zex,rex和Sex(2l+2)分別為:

(14)

(15)

(16)

將式(14)～式(16)中的zex,rex,Sex代入式(10)中，得到M2m+2,2l+2(r,z)的積分表達(dá)式

(17)

將zex,rex,Sex代入式(12)，也可得到M2m+2,2l+2(r,z)的傍軸近似解

(18)

式(18)中，M2m+2,2l+2,p中的下標(biāo)p表示傍軸近似.相對(duì)于式(17)的非傍軸傳輸?shù)姆e分解，式(18)所給出的傍軸近似解析解形式更為簡(jiǎn)單.在要求不是特別嚴(yán)格時(shí)，在滿足傍軸近似的條件下，式(18)能更簡(jiǎn)單地處理馬丟-高斯光束的傳輸問題.

(19)

通過式(19)即可精確計(jì)算第一類奇型馬丟-高斯光束在自由空間傳輸時(shí)的光場(chǎng)分布.如當(dāng)r=0時(shí)，通過式(19)可求出第一類奇階馬丟-高斯光束的軸上的光場(chǎng)分布

(20)

圖1 第一類奇型馬丟-高斯光束的軸上光強(qiáng)分布

很多時(shí)候，光束傳輸時(shí)的軸上光場(chǎng)分布最值得人們關(guān)注.故用于計(jì)算軸上光場(chǎng)分布的式(20)有一定的理論價(jià)值.

根據(jù)虛源法推導(dǎo)得到了馬丟-高斯光束在自由空間傳輸?shù)姆e分表達(dá)式，筆者數(shù)值模擬了第一類奇型馬丟-高斯光束的軸上光強(qiáng)分布,具體如圖1所示.從圖1中可以看出，當(dāng)傳輸距離較小時(shí)，非傍軸光強(qiáng)分布的結(jié)果與傍軸理論計(jì)算的結(jié)果有明顯差異;隨著傳輸距離的增加，傍軸的近似計(jì)算結(jié)果與考慮校正項(xiàng)的非傍軸計(jì)算結(jié)果越來越接近.說明根據(jù)衍射理論來研究光束的近場(chǎng)傳輸時(shí)，利用非傍軸理論才能得出正確的結(jié)果.在研究光束的遠(yuǎn)場(chǎng)傳輸時(shí)，利用傍軸理論的近似法才能得到相對(duì)正確的研究結(jié)果.

筆者用相同的方法得到第二類奇型馬丟-高斯光束在自由空間傳輸?shù)木_積分表達(dá)式為

(21)

第一類偶型馬丟-高斯光束在自由空間傳輸?shù)木_積分表達(dá)式為

(22)

第二類偶型馬丟-高斯光束在自由空間傳輸?shù)木_積分表達(dá)式為

(23)

綜上所述，式(21)～式(23)即為另外3種馬丟-高斯光束傳輸?shù)姆e分表達(dá)式，相應(yīng)的傍軸近似解和軸上光場(chǎng)分布解也可相應(yīng)得到，限于篇幅，本文就不再贅述.

2 結(jié) 論

相對(duì)于貝塞爾光束、艾里光束等(準(zhǔn))無衍射光束，馬丟光束最大的特點(diǎn)是具有多種光學(xué)形態(tài).這種多樣性并不僅僅體現(xiàn)在馬丟光束本身就有4種形態(tài)，而且它的4種形態(tài)的每一種光學(xué)形態(tài)都隨著橢圓離焦量的改變而改變.因此，馬丟光束的構(gòu)造和產(chǎn)生在實(shí)際的科學(xué)研究中具有更大的實(shí)用性和靈活性.

但在數(shù)學(xué)上描述馬丟光束的馬丟函數(shù)也更加復(fù)雜，這大大阻礙了人們研究馬丟光束簇的傳輸特性.為解決這一光學(xué)難題，筆者根據(jù)光波傳播的獨(dú)立性原理和疊加性原理，在數(shù)學(xué)上將橢圓坐標(biāo)系下的馬丟函數(shù)展開為極坐標(biāo)系下的貝塞爾函數(shù).在貝塞爾函數(shù)關(guān)系的基礎(chǔ)上，利用虛源法精確推導(dǎo)了4種馬丟-高斯光束的積分表達(dá)式.該組積分表達(dá)式為精確、定量地研究馬丟-高斯光束的傳輸特性奠定了理論基礎(chǔ).另外，本文采用級(jí)數(shù)展開(對(duì)應(yīng)光波傳播的獨(dú)立性原理和疊加性原理)與虛源技術(shù)相結(jié)合的手段，為研究其他復(fù)雜光束的傳輸問題提供了很好的方法與思路.

[1]Durnin J,Miceli J J,Eberly J H.Diffraction-free beams[J].Physical Review Letters,1987,58(15):1499-1501.

[2]Durnin J.Exact solutions for nondiffracting beams[J].Journal of the Optical Society of America A,1987,4(4):651-654.

[3]Gutiérrez-Vega J C,Iturbe-Castillo M D,Chávez-Cerda S.Alternative formulation for invariant optical fields:Mathieu beams[J].Optics Letters,2000,25(20):1493-1495.

[4]Ye Fangwei,Mihalache D,Hu B.Elliptic vortices in composite Mathieu lattices[J].Physical Review A,2009,79(5):1744-1747.

[5]Kartashov Y V,Egorov A A,Vysloukh V A,et al.Shaping soliton properties in Mathieu lattices[J].Optics Letters,2006,31(2):238-240.

[6]Chafiq A,Hricha Z,Belafhal A.A detailed study of Mathieu-Gauss beams propagation through an apertured ABCD optical system[J].Optics Communications,2006,265(2):594-602.

[7]Gutiérrez-Vega J C,Bandres M A.Normalization of the Mathieu-Gauss optical beams[J].Journal of the Optical Society of America A,2007,24(1):215-220.

[8]Deschamps G A.Gaussian beam as a bundle of complex rays[J].Electronics Letters,1971,7(23):684-685.

[9]Felsen L B.Evanescent waves[J].Journal of the Optical Society of America,1976,66(8):751-760.

[10]Zhang Y,Song Y,Chen Z,et al.Virtual sources for a cosh-Gaussian beam[J].Optics Letters,2007,32(3):292-294.

[11]Li Y,Lee H,Wolf E.New generalized Bessel-Gaussian beams[J].Journal of the Optical Society of America,2004,21(4):640-646.

[12]Yan S,Yao B,Lei M,et al.Virtual source for an Airy beam[J].Optics Letters,2012,37(22):4774-4776.

[13]Deng D,Guo Q.Elegant Hermite-Laguerre-Gaussian beams[J].Optics Letters,2008,33(11):1225-1227.

[14]Chafiq A,Hricha Z,Belafhal A.Parametric characterization of Mathieu-Gauss beams[J].Optics Communications,2008,282(13):2590-2594.

[15]Couture M,Belanger P.From Gaussian beam to complex-source-point spherical wave[J].Physical Review A,1981,24(1):355-359.

[16]Seshadri S R.Virtual source for the Bessel-Gauss beam[J].Optics Letters,2002,27(12):988-1000.

[17]王碩琛,吳逢鐵,馮聰.一種產(chǎn)生Mathieu 光束的光學(xué)系統(tǒng)[J].光學(xué)學(xué)報(bào),2017,37(2):0205002.

[18]謝曉霞,李冬,吳逢鐵.橢圓孔徑與軸棱錐系統(tǒng)產(chǎn)生帶狀近似無衍射光束[J].光學(xué)學(xué)報(bào),2015,35(1):0126001.