劉　華

(中國科學(xué)院 長春光精密機(jī)械與物理研究所 光電技術(shù)研發(fā)中心，吉林 長春 130033)

?

利用漢克爾變換設(shè)計高斯光束整形衍射元件的應(yīng)用研究

劉華

(中國科學(xué)院 長春光精密機(jī)械與物理研究所 光電技術(shù)研發(fā)中心，吉林 長春 130033)

衍射光學(xué)元件由于可以實現(xiàn)對高斯光束的整形而被重視，其通常的設(shè)計方法為G-S算法，由于使用傅里葉變換運算量大、費時長，將快速漢克爾變換應(yīng)用到這些算法中可以極大地提高運算速度，節(jié)省運算時間，為設(shè)計復(fù)雜的光束整形元件提供了高效、可行的方法。本文利用該種方法設(shè)計針對中心波長為775 nm、光束束腰口徑為6 cm的激光器，成功設(shè)計了一個具有二階相位的折衍混合光學(xué)元件。僅單獨這一片元件，既可在距離其35 m處得到一半徑為200 μm的圓形平頂光斑，均方根誤差D<0.021。當(dāng)抽樣值取215時，在普通PC機(jī)上運行時間僅為20.05 s，大大節(jié)省了優(yōu)化設(shè)計時間(整個優(yōu)化設(shè)計過程往往需要幾十次甚是上百次這種運算)。同時利用離子刻蝕技術(shù)加工了該折衍混合元件，并進(jìn)行了實際測試，結(jié)果與設(shè)計值基本相符，整形效果較好。這種單片的整形元件不僅整形效果好，還有利于與激光器的集成，簡化系統(tǒng)的調(diào)節(jié)。

衍射元件；光束整形；高斯光束；漢克爾變換

1　引　言

伴隨著激光技術(shù)的發(fā)展，激光器的功率不斷提升，應(yīng)用也日益廣泛，如激光焊接、刻蝕、掃描、打標(biāo)[1-4]等。這些應(yīng)用都需要滿足激光光束會聚以后，光斑能量呈均勻分布，即形成平頂形狀的光斑。然而實際上由激光器直接出射的激光光束，其能量分布往往呈高斯?fàn)?，中心?qiáng)、邊緣弱，不能滿足使用需求，需要利用光學(xué)方法對其進(jìn)行光束整形。近年來，利用衍射光學(xué)元件進(jìn)行光束勻化整形取得了極大的進(jìn)展，人們已經(jīng)利用位相片、全息片、二元光學(xué)器件等衍射元件獲得了光強(qiáng)分布均勻的光斑[5-9]。用于光束整形的衍射元件其傳統(tǒng)的設(shè)計方法有幾何光學(xué)法、模擬退火法、遺傳算法、G-S算法等[10-11]。其中幾何光學(xué)方法計算不夠精確，求得的結(jié)果往往作為G-S算法的初始值。模擬退火法、遺傳算法是基于搜索極值的優(yōu)化算法，魯棒性強(qiáng)，要想找到最優(yōu)解, 需要消耗大量的搜索時間，收斂速度很慢。G-S算法雖然收斂速度快，但是由于使用傅里葉變換，計算量大，計算周期長，很難實現(xiàn)快速、準(zhǔn)確的設(shè)計和結(jié)果分析。另外在光束整形的過程中，很多情況除了需要對光束進(jìn)行能量調(diào)節(jié)外，還需要在遠(yuǎn)場將能量聚焦。為了實現(xiàn)這兩者的作用，通常做法是在衍射元件的后方再加入一片傅里葉變換透鏡[12]。這樣會導(dǎo)致光路中元件增多，調(diào)整不便，集成度不高。即使是做成單片，也基本是在平面基底上制作，導(dǎo)致衍射原件的位相周期比較小，為加工制作帶來不便[13]。

旋轉(zhuǎn)對稱性質(zhì)的傅里葉變換可以用快速漢克爾變換來代替[14-15]，但大部分都是用于激光光束傳播過程的計算，目前還沒有系統(tǒng)的將該算法應(yīng)用于G-S整形算法中的文獻(xiàn)報道。本文將其應(yīng)用于G-S優(yōu)化算法中，大大提高了計算速度，節(jié)省了計算時間。同時將衍射元件基底的另一表面制作成帶有曲率的球面，形成一個折衍混合光學(xué)元件，能夠同時完成光斑的勻化和聚焦功能。針對中心波長λ=775 nm的激光器，利用數(shù)學(xué)計算軟件Matlab成功設(shè)計了一個具有二階相位的折衍混合光學(xué)元件，在距離光束束腰35 cm處得到一個半徑為200 μm的圓形平頂光斑，均方根誤差D小于0.021。當(dāng)抽樣值取215時，在普通PC機(jī)上運行時間僅為20.05 s，如果使用Matlab軟件中自帶的傅里葉變換和逆變換函數(shù)，運算時間會超過10 min。而優(yōu)化設(shè)計過程往往會需要進(jìn)行幾十次甚至上百次的運算，這會極大地節(jié)省優(yōu)化設(shè)計時間。最后利用離子刻蝕技術(shù)將衍射位相刻制在了以平凸透鏡為基底的平面上，形成了理想的高斯光束整形折衍混合元件。在實驗室搭建了光斑測試系統(tǒng)，利用激光斑探測器接收光斑，表明設(shè)計結(jié)果理想，元件整形效果很好。同時系統(tǒng)只有一個光學(xué)元件，簡單易調(diào)，便于集成。

2　設(shè)計方法理論分析

圖1　光束整形的原理圖 Fig.1　Schematic diagram of beam shaping

(1)

式中，r為入射光線經(jīng)過P1元件表面的徑向坐標(biāo)，w0為高斯光束束腰半徑。|f0(r)|為輸入光波的實振幅，φ0(r)為輸入光波的位相。衍射光學(xué)元件P1的透過率函數(shù)為

(2)

則入射光波在經(jīng)過衍射光學(xué)元件P1以及透鏡后的復(fù)振幅分布為

(3)

其中

對f1(r)進(jìn)行傅里葉變換變換，得輸出函數(shù)為

(4)

式中，p為入射光線經(jīng)過P2面的徑向坐標(biāo)，J0為零階貝塞爾(Besse1)函數(shù)，z為輸入面和輸出面的間距，C為一常數(shù)。|g1(p)|為輸出光波的實振幅，φg1(p)為輸出光波的位相。設(shè)計要求輸出的實振幅為|g0(p)|，利用G-S迭代算法進(jìn)行優(yōu)化計算，便可以求得衍射光學(xué)元件的透過率函數(shù)。G-S迭代算法是一種常見的位相恢復(fù)優(yōu)化計算，它的每次循環(huán)均由下列4個基本步驟組成

(5)

(6)

(7)

(8)

式中，β為實常數(shù)，一般取0<β<1。

以上計算過程是一個迭代過程，其中每次迭代都需要進(jìn)行傅里葉變化和傅里葉逆變化，如果直接利用數(shù)學(xué)計算軟件中的對對應(yīng)函數(shù)來計算，伴隨著抽樣數(shù)值的增大，其計算時間以指數(shù)倍增大，甚至?xí)斐伤罊C(jī)。在實際計算中，如果使用快速漢克爾變換，可以大大提高計算速度，節(jié)約計算時間。根據(jù)式(5)可得

(9)

當(dāng)對輸入輸出函數(shù)進(jìn)行N次抽樣時，

(10)

(11)

(12)

式中，fk(rN)≡0，k0=[2exp(a)+exp(2a)]/{[1+exp(a)]2[1-exp(-2a)]}，r0=[1+exp(a)]exp(-aN)/2。a可以通過解方程a=-ln[1-exp(-a)]/(N-1)求得。

循環(huán)過程中，利用均方根誤差D作為評價函數(shù)，一旦D滿足要求，循環(huán)結(jié)束，求得衍射光學(xué)元件的透過率函數(shù)。

(13)

式中，|gk(R)|i為計算中產(chǎn)生的輸出實振幅分布的第i個抽樣點值，|g0(R)|i為設(shè)計要求輸出的實振幅分布的第i個抽樣點值，N為抽樣點數(shù)。

3　實驗與結(jié)果

圖2　衍射元件2階位相分布圖 Fig.2　Binary phase of the diffractive optical element

圖3　接收面上能量分布的設(shè)計結(jié)果 Fig.3　Designed energy distribution on the detector plane

圖4　掩膜版結(jié)構(gòu)圖與折衍混合元件實物 Fig.4　Configure of diffractive templet, binary refractive and diffractive optical element

我們在實驗室的光學(xué)平臺上搭建了測試系統(tǒng)，其光路如圖5所示，實物圖如圖6所示，能量呈高斯形狀分布的775 nm激光器經(jīng)過擴(kuò)束鏡后形成束腰口徑為60 mm的準(zhǔn)直光束。將折衍混合光學(xué)元件置于束腰處，在距離其350 mm光斑會聚的位置，利用CCD類型的光電探測器測量了整形后的光能量分布情況。由于只用一片折衍混合元件即可實現(xiàn)光束的整形和聚光，該系統(tǒng)搭建過程簡單，調(diào)節(jié)方便。同時這個元件非常容易與激光器集成起來應(yīng)用。

圖5　測試光路圖 Fig.5　Test optical system

通過使用CCD類型的光電探測器，我們首先直接測量了775 nm的激光光斑能量分布，其分布結(jié)果見圖7，該圖為三維能量分布圖，平面坐標(biāo)為光斑的徑向尺寸，單位為100 μm，縱向坐標(biāo)為能量值。可以明顯看出中心能量最高，邊緣能量低，呈高斯形狀分布。然后按照實物光路結(jié)構(gòu)測量了經(jīng)過折衍混合光學(xué)元件整形后的光斑，其能量具體分布狀況見圖8(該圖為三維能量分布圖，平面坐標(biāo)為光斑的徑向尺寸，單位為100 μm，縱向坐標(biāo)為能量值)，其頂部基本平坦，只是在光斑邊緣處有一個小的尖峰，這主要是由于衍射元件的環(huán)帶位置加工存在誤差以及實驗裝調(diào)時各個器件相對位置存在誤差等原因造成。經(jīng)過計算，該光斑的尺寸為200 μm，均方根誤差D=0.025，結(jié)果與設(shè)計值基本相符，整形效果很好，能夠滿足應(yīng)用需求。說明本方法設(shè)計的整形元件不僅速度快，而且效果理想。就是說，在利用數(shù)學(xué)計算軟件Matlab優(yōu)化設(shè)計激光光束整形元件的過程中，將快速漢克爾變換替代其原有的傅里葉變換，可以大大提高G-S優(yōu)化算法的速度，節(jié)省計算時間，同時能夠獲得理想的整形效果，把高斯分布的激光光束變換會聚后，得到平頂分布的光斑。

圖7　高斯光束的能量分布圖 Fig.7　Energy distribution of Gaussian beam

圖8　整形后能量分布圖 (3維立體) Fig.8　Energy distribution after beam shaping(3D)

4　結(jié)　論

本文針對中心波長λ=775 nm的激光器，成功設(shè)計了一個具有二階相位的衍射光學(xué)元件，可以在距離光束束腰350 mm處得到一半徑為200 μm的圓形平頂光斑。均方根誤差D<0.021，當(dāng)抽樣值取215時，在普通PC機(jī)上運行時間僅為20.05 s。而使用傅里葉變換的計算時間

超過10 min。這種單次的計算時間的縮短，對于優(yōu)化過程非常重要，因為整個優(yōu)化設(shè)計過程往往需要反反復(fù)復(fù)進(jìn)行幾十次甚至上百次的計算，尤其是經(jīng)驗匱乏的新手，需要的計算次數(shù)就更多。同時將衍射元件與折射元件相結(jié)合，形成折衍射混合光學(xué)元件，單片就能完成對激光整形和會聚的作用，有利于及系統(tǒng)的集成和調(diào)節(jié)。最后通過對該衍射元件的制作和測試，表明設(shè)計結(jié)果理想，元件整形效果較好。

[1]DOPARRE.M,COLUB M A,LUDGE B,etal.. Investigation of computer-generated diffractive beam shapers for flattening of single-modal CO2laser beams[J].Appl.Opt.,2005,34(14):2487-2497.

[2]吳環(huán)寶,王肇圻,傅汝廉,等.CO2激光雕刻機(jī)的折/衍混雜f-θ鏡設(shè)計[J]. 紅外與激光工程,2006,35(3):351-358.

WU H B,WANG Z Q,FU R L,etal.. Design of hybrid refractive-diffractivef-θlens working on CO2laser drilling system[J].InfraredandLaserEngineering,2006,35(3):351-358.(in Chinese)

[3]HERMSDORF J,KAIERLE S. Benefits of optogalvanic effect and beam shaping for laser stabilized GMA welding[J].ChineseOptics,2014,7(1):112-117.

[4]肖龍勝,唐霞輝,秦應(yīng)雄,等.2 kW射頻板條CO2激光器輸出光束整形特性研究[J].中國激光,2014,41(4)：1-7.

XIAO L SH,TANG X H,QIN Y X,etal.. Shaping characteristic of output beam of 2 kW radio frequency slab CO2laser[J].ChineseJ.Laser,2014,41(4)：1-7.(in Chinese)

[5]許超，張靜娟.用于圓對稱光束波前變換的位相型光學(xué)系統(tǒng)[J].物理學(xué)報,1993,42(8)：1245-1251.

XU CH,ZHANG J J. A phase only optical system for wavefront transmission of rotational symmetric beam[J].ActaPhysicSinica,1993,42(8):1245-1251.(in Chinese)

[6]包紅春,鄔敏賢,金國藩，等.二元光學(xué)波面變形器件的研究[J].光子學(xué)報,1994,14(9)：988-991.

BAO H CH,WU M X,JIN G F,etal.. Binary optical element for laser beam shaping[J].ActaOpticaSinica,1994,14(9):988-991.(in Chinese)

[7]劉強(qiáng)，李珂，楊佳，等.衍射光學(xué)元件設(shè)計中的插值迭代方法[J].計算物理,2009,26(5):731-736.

LIU Q,LI K,YANG J,etal.. Interpolation iterative arithmetic in design of diffractive optical elements[J].ChineseJ.ComputationalPhysics,2009,26(5):731-736.(in Chinese)

[8]于曉晨，胡家升，王連寶.提高激光束整形質(zhì)量的新方法[J].中國激光,2012,39(1):1-6.

YU X CH,HU J SH,WANG L B. New methods for improving the quality of laser beam shaping [J].ChineseJ.Laser,2012,39(1):1-6.(in Chinese)

[9]龐輝,應(yīng)朝福 ,范長江,等.用于光束整形的衍射光學(xué)元件設(shè)計的混合算法[J].光子學(xué)報,2010,39(6):977-981.

PANG H,YING CH F,FAN CH J,etal.. Design diffractive optical elements for beam shaping with hybrid algorithm[J].ActaPhoTonIcaSInica,2004,31(12):8-11.(in Chinese)

[10]陳燾,羅崇泰,王多書,等.遺傳算法用于衍射光學(xué)元件的優(yōu)化設(shè)計[J].光電工程,2004,31(12):8-11.

CHEN S,LUO C T,WANG D SH,etal.. Genetic algorithms used for optimization design of diffractive optical elements[J].Opto-ElectronicEngineering,2004,31(12):8-11.(in Chinese)

[11]ZHOU G Y,YUAN X C,PHILIP D,etal.. Design of diffractive phase elements for beam shaping:hybrid approach[J].J.Opt.Soc.Am.A,2001,18(4):791-799.

[12]李鳳友,李紅軍,盧振武,等.激光直寫系統(tǒng)焦斑整形的研究[J].光學(xué) 精密工程,2001,9(1):14-18.

LI F Y,LI H J,LU ZH W,etal.. Focused spot shaping in laser writing system[J].Opt.PrecisionEng.,2001,9(1):14-18.(in Chinese)

[13]徐兵,王雪輝,陳林森,等.激光束整形為正方框形光束的DOE的設(shè)計及實驗[J].激光技術(shù),2005,29(2):169-182.

XU B,WANG X H,CHEN L S,etal.. The designation and experiment of the DOE applied to shape the laser beam to square frame shape[J].LaserTechnology,2005,29(2)：169-182.(in Chinese)

[14]ZHANG D W,YUAN X C,NGO N Q,etal.. Fast Hankel transform and its application for studying the propagation of cylindrical electromagnetic fields[J].OpticsExpress,2002,10(12) :521-525.

[15]MAGNI V,CERULLO G,SILVESTRI S D. High-accuracy fast Hankel translation for optical beam propagation[J].J.Opt.Soc.Am.A,1992,9(11)：2031-2033.

Gaussian beam shaping diffractive optical element designed by Hankel transformation

LIU Hua

(Opto-electronicsTechnologyCenter,ChangchunInstituteofOptics,FineMechanicsandPhysics,ChineseAcademyofSciences,Changchun130033,China)

，E-mail：girlliuhua@sohu.com

Because diffractive optical element can shape Gaussian beam, it is attached great importance. The conventional design methods such as G-S algorithm tend to need a large amount of calculation and a long time, while Hankel transformation can increase calculation speed and save time. We successfully designed a binary refractive and diffractive optical element applying Hankel transformation into G-S algorithm. When the diffractive optical element was placed at the waist of the laser beam with the wavelength of 775 nm and the waist diameter of 6 cm, a uniform spot could be gotten at the detector plane 35 cm behind it. Its radius was 200 μm and root-mean-square error was smaller than 0.021. The calculation time was only 20.05 s with PC, saving much time during the process of optimization, in which hundreds of calculation were needed. We fabricated the refractive and diffractive optical element by ion beam etching, and tested it in lab. The experiment results agreed well with the caculation. This kind of element not only can get good spot, but also is suitable of integration with laser.

diffractive optical element;beam shaping;Gaussian beam;Hankel transformation

2015-12-15；

2016-01-11

國家自然科學(xué)基金重點資助項目(No.61137001)

2095-1531(2016)02-0277-07

O436

A

10.3788/CO.20160902.0277

劉華(1976—)，女，遼寧撫順人，博士，副研究員，2001年于東北師范大學(xué)獲得學(xué)士學(xué)位，2006年于中國科學(xué)院長春光學(xué)精密機(jī)械與物理研究所獲得博士學(xué)位，主要從事光學(xué)檢測、光學(xué)設(shè)計、衍射光學(xué)、非成像光學(xué)方面的研究。E-mail:girlliuhua@sohu.com

Supported by Major Project of National Natural Science Foundation of China(No.61137001)