蔡建文

（1.常州工學院 光電工程學院，常州213002；2.常州現(xiàn)代光電技術研究院，常州213002；3.常州市光電子材料與器件重點實驗室，常州213002）

引言

飛秒激光是20世紀后期發(fā)展起來的一種新型激光技術，它瞬間功率高［1-2］。與傳統(tǒng)激光系統(tǒng)相比，它是一種冷加工激光，加工材料在飛秒激光作用下熱變形小，可獲得亞微米級寬度的直線和點，已成為激光微加工領域的發(fā)展方向之一。飛秒激光微加工是采用逐點掃描加工方式，其焦點光斑光場的分布直接影響到加工的質量。由于受到激光束腰及透鏡衍射效應的影響，焦斑光場在空間分布呈橢球形，在三維光致聚合微成形加工中，用高倍物鏡（100×，NA＝1.4）聚焦，焦斑軸向尺寸也約是其橫向尺寸的3倍［3］。將這樣的焦斑應用于激光微加工，將大大降低激光加工的精度。為了提高激光微加工質量，必須對焦斑進行整形，實現(xiàn)軸向超分辨。1954年Toraldo提出超分辨的概念［4］后，人們研究了一系列的振幅型和位相型超分辨元件和方法［5-9］，通過模擬證明這些振幅型和相位型元件可提高焦斑軸向或橫向的分辨率。國內在超分辨技術方面，上海光機所設計出復振幅光瞳濾波器可以分別實現(xiàn)焦斑的軸向或橫向超分辨，與共焦系統(tǒng)結合后在一定程度上提高了三維分辨率［10］；上海理工大學基于嚴格的光學成像矢量衍射理論，通過詳細研究非對稱三區(qū)復振幅型光瞳濾波器的內外環(huán)歸一化半徑、各環(huán)相位分布和第一層透過率對y方向分辨率增益比、斯特爾比和第一旁瓣與主瓣相對強度的影響，設計了一種非對稱三區(qū)相位型光瞳濾波器［11］。由于激光光斑軸向超分辨可以大大提高飛秒激光的微加工精度，本文對超衍射理論進行了詳細地分析，采用優(yōu)化算法設計了對應于縱向超衍射的2種二元相位元件；仿真結果表明，在飛秒激光微加工系統(tǒng)中，這種位相元件的超衍射效果比較好。

1 微加工軸向超分辨機理

在飛秒激光微加工系統(tǒng)中，為了獲取較好的沿光軸方向的超衍射效果，在微加工物鏡前加上一個位相調制的Toraldo型光瞳濾波器（稱為位相板）來進行光斑形狀和大小調制，如圖1所示。

圖1 激光微加工系統(tǒng)中加入相位板Fig.1 Adding phase plate in laser micromachining system

按照衍射理論，聚焦光斑的大小由波長、數(shù)值孔徑、孔徑形狀決定。但是，即使這些量已經(jīng)固定，光斑大小仍然可以進一步縮小，以獲得超分辨率。如果使用相位板整形光束，可以使聚焦光斑在橫向和縱向小于衍射極限。

把位相板放置于物鏡前方，可在入瞳處改變入射光的位相大小，從而可減小焦平面處光斑尺寸的目的，利用菲涅爾衍射公式，可獲得焦平面附近歸一化光衍射場的點擴散函數(shù)［12］：

由于位于焦平面附近，因此可以將上述點擴散函數(shù)進行泰勒展開，從而可得到：

光瞳函數(shù)的m階矩以Im來表征：

由公式（2）和（3），可以得到整形后的軸向光強表達式為

為了表征強度分布，可采用3個基本量來描述相位板的整形結果：縱向增益GA、峰值能量比S和旁瓣能量MA。其中，GA定義為整形后和整形前軸向的光斑尺寸之比；S定義為整形后和整形前中心峰值能量之比；MA定義為整形過后除去中心光強外的最大旁瓣光強與中心峰值強度之比。

峰值能量比S、縱向增益GA可用下面表達式來定義（a表示加相位板以后的值，c表示未加相位板時候的值）：

對光強分布函數(shù)（4）求導可得：

對表達式（5）進行簡化，從而獲得光斑峰值位置為

式中uF為離焦量，即為偏移焦平面的位移量。

因此沿光軸方向的強度分布的峰值能量比S為

焦斑尺寸GA為

峰值能量比S、縱向增益GA都是相位板設計的重要參考指標。通過這些參數(shù)的優(yōu)化設計，可以獲得較理想的參考模型。

要獲得理想的結果，設計相位板還必須考慮焦斑的大小，峰值能量比S和旁瓣能量MA等，并需要確定在滿足約束條件下，對相應的參數(shù)進行優(yōu)化，因此優(yōu)化問題可以表述為

需滿足約束：

S＞a（a為可以接受的最小強度值）

2 設計及仿真

相位板的性能取決于環(huán)帶數(shù)量、每個環(huán)的半徑和相位。設計相位板過程中必須要兼顧光斑尺寸、峰值能量比和旁瓣能量等，并且需要在滿足上述約束條件的情況下找到一組最佳的優(yōu)化參數(shù)，因此采用遺傳算法來獲取相位板參數(shù)。遺傳算法在適應度函數(shù)選擇不當?shù)那闆r下有可能收斂于局部最優(yōu)，而不能達到全局最優(yōu)，因此這里采用全局優(yōu)化算法［13］與遺傳算法相結合，來獲得最優(yōu)解。

通過在優(yōu)化過程中調整遺傳算法的交叉率、變異率以及約束條件，采用Matlab工具得到一種0-π結構和另一種非0-π結構的位相板設計方案，如表1所示。兩種位相板調制后的軸向光強分布和橫向光強分布分別如圖2和圖3所示。對這兩種方案進行比較，從兩種相位板的衍射效果來看，非0-π結構的設計方案1比0-π結構的設計方案2具有更高的峰值能量比，更小的旁瓣能量，但光斑縱向尺寸壓縮效果基本一致，艾利斑壓縮GA≈0.75，即光斑縱向大小可壓縮至艾麗斑的75%，而橫向尺寸基本保持不變。從圖2和圖3可以看出，縱向旁瓣能量比橫向旁瓣能量要大得多，并且當相位板設計參數(shù)發(fā)生微量改變時，縱向旁瓣能量也更會產生較大的波動。雖然在縱向會產生比較大的旁瓣，然而在飛秒激光微細加工系統(tǒng)中，由于存在雙光子吸收效應，并且雙光子吸收是一種非線性效應，這可以很好地抑制縱向相位板的旁瓣副作用。

表1 優(yōu)化設計參數(shù)和表征參數(shù)Fig.1 Optimization of design parameters and characterization parameters

圖3 沿橫向超衍射光斑與艾利斑的光強分布Fig.3 Intensity distribution of super-resolution pattern and Airy pattern along lateral direction

圖2 沿光軸方向超衍射光斑與艾利斑的光強分布Fig.2 Intensity distribution of super-resolution pattern and Airy pattern along axial direction

3 結論

為了在激光維納加工中提高加工質量，采用菲涅爾衍射公式對激光光斑整形進行了理論分析，基于遺傳算法和設計約束條件通過Matlab設計了兩種二元相位元件，一種為0-π結構四環(huán)相位板，另一種為非0-π結構四環(huán)相位板，采用這兩種位相板調制后激光縱向光斑均可實現(xiàn)壓縮，縱向光斑大小分別可壓縮至艾麗斑的76%和75%，峰值能量比分別為0.39和0.42，旁瓣能量分別為0.64和0.41。從仿真結果可以看出，這兩種相位板均可應用于飛秒激光微加工。

［1］ Liu Qing，Chen Junjun，Guo Lili.Micromachining of wide-band optical materials with ultra-short laser pulses［J］.Journal of Applied Optics，2006，27（5）：428-432.劉青，陳鈞均，郭麗麗.超短激光脈沖對寬帶光學物質的微加工［J］.應用光學，2006，27（5）：428-432.

［2］ Wang Xianhua，Chen Feng，Yang Qing，et al.Experiment research on femtosecond-laser-induced optical micro-cavity in transparent materials［J］.Journal of Applied Optics，2007，28（2）：191-194.王先華，陳烽，楊青，等.飛秒激光在透明介質中誘導光學微腔的實驗研究［J］.應用光學，2007，28（2）：191-194.

［3］ Sun Hongbo，Tanaka Tomokazu，Kawata Satoshi.Three-dimensional focal spots related to two-photon excitation［J］.Appl.Phys.Lett.，2002，80（20）：3673-3675.

［4］ Toraldo di Francia G.Super-gain antennas and optical resolving power［J］.Nuovo Cimento Supplemento，1952，9（3）：426-435.

［5］ Cheng Ya，Sugioka Koji，Midorikawa Katsumi，et al.Control of the cross-sectional shape of a hollow microchannel embedded in photostructurable glass by use of a femtosecond laser［J］.Opt.Lett.，2003，28：55-57.

［6］ Ams M ，Marshall G D，Spence D J，et al.Slit beam shaping method for femtosecond laser direct-write fabrication of symmetric waveguides in bulk glasses［J］.Opt.Express，2005，13：5676-5681.

［7］ Cerullo G ，Osellame R ，Taccheo S，et al.Femtosecond micromachining of symmetric waveguides at 1.5mm by astigmatic beam focusing［J］.Opt.Lett.，2002，27，1938-1940.

［8］ Boyer G.New class of axially apdizing filters for confocal scanning microscopy［J］.J.Opt.Soc.Am.（A），2002，19（3）：584-589.

［9］ de Juana D M，Oti J E，Canales V F，et al.Transverse or axial superresolution in a 4Pi-confocal microscope by phase-only filters［J］.J.Opt.Soc.Am（A），2003，20（11）：2172-2178.

［10］ Wang Mei，Yun Maojin，Liu Liren，et al.Design of the super-resolution pupil filter in the confocal system［J］.Acta Optica Sinica，2011，31（6）：189-193.王美，云茂金，劉立人，等.共焦顯微系統(tǒng)中光學超分辨光瞳濾波器的設計［J］.光學學報，2011，31（6）：189-193.

［11］ Guo Shuwen，Guo Hanming，Zhuang Songlin.Onedimensional transverse superresolution achieved by three-zone dissymmetrical pupil filter［J］.Acta Photonica Sinica，2008，37（11）：2222-2226.郭舒文，郭漢明，莊松林.非對稱三區(qū)光瞳濾波器實現(xiàn)一維橫向超分辨［J］.光子學報，2008，37（11）：2222-2226.

［12］ Born M，Wolf E.Principles of optics［M］.New York：Pergamon，1975.

［13］ Gu Jun.Global optimization for satisfiability（SAT）problem［J］.IEEE Transactions on Knowledge and Data Engineering，1994，6（3）：361-381.