張繼魁， 時家明， 苗 雷， 王啟超， 趙大鵬， 曾 杰

(電子工程學院 脈沖功率激光技術(shù)國家重點實驗室， 安徽 合肥 230037)

?

近中紅外與1.06 μm和1.54 μm激光兼容隱身光子晶體研究

張繼魁， 時家明*， 苗 雷， 王啟超， 趙大鵬， 曾 杰

(電子工程學院 脈沖功率激光技術(shù)國家重點實驗室， 安徽 合肥 230037)

為實現(xiàn)飛行器高溫部分的紅外與激光兼容偽裝，設(shè)計了一種基于一維光子晶體的近中紅外與1.06 μm和1.54 μm激光兼容隱身材料?；诒∧さ膫鬏斁仃嚪ê彤愘|(zhì)結(jié)構(gòu)理論，拓展了光子晶體的禁帶寬度，使之覆蓋近中紅外波段；隨后，利用摻雜原理，在光子晶體周期結(jié)構(gòu)中引入了兩種缺陷。結(jié)果顯示，在1～5 μm的帶隙中出現(xiàn)了波長分別為1.06 μm和1.54 μm的缺陷模，反射率分別為1.21%和1.79%，這種具有“光譜挖空”特性的光子晶體可以實現(xiàn)近中紅外與1.06 μm和1.54 μm激光兼容隱身。

光子晶體； 近中紅外； 激光； 兼容偽裝

1 引 言

火箭或噴氣飛機等飛行器在高速飛行時，由于其發(fā)動機工作和氣動加熱等原因，某些部位的溫度很高，例如火箭的燃燒室溫度可達2 000～3 000 ℃。由維恩位移定律[1]可知，其輻射峰值波長處于0.966～1.449 μm的近紅外范圍內(nèi)；而飛機尾噴管以及其表面反射的太陽光輻射通常處于中紅外波段[1]。紅外探測設(shè)備通常利用飛行器高溫部位的強烈紅外輻射進行目標識別與跟蹤。因此，有效抑制飛行器的紅外輻射，實現(xiàn)近中紅外波段的兼容是非常必要的。

此外，常見的用于測距和制導的激光器包括波長為1.06 μm的Nd∶YAG激光器和1.54 μm的摻鉺激光器，這兩類激光器通常以主動探測方式為主，基于目標表面的高反射特性，通過探測回波信號來確定目標的位置。而飛行器表面通常由金屬或合金構(gòu)成，其在上述波長上具有較高的反射率。因此，對于此類目標，一方面需要抑制其在近中紅外波段的發(fā)射率ε，而另一方面還需降低其表面在上述激光波長處的反射率ρ。通常，在不考慮吸收的情況下，ε+ρ=1，可以看出，發(fā)射率和反射率的同時降低不可能實現(xiàn)。

然而，光子晶體作為一種由不同介電常數(shù)的介質(zhì)周期性排列組成的人工晶體，其特有的“帶隙[2]”和“光子局域[3]”性質(zhì)使實現(xiàn)多波段兼容偽裝成為可能。帶隙又稱禁帶，是特定頻率范圍內(nèi)對電磁波的高反射帶。光子局域是光子晶體禁帶中出現(xiàn)的頻率極窄的缺陷態(tài)，它是由光子晶體的周期性或者對稱性結(jié)構(gòu)遭到破壞造成的。與缺陷態(tài)頻率一致的光子會被禁錮在缺陷附近，一旦離開缺陷位置，光就會大幅度衰減。

基于此，本文設(shè)計了一種在近中紅外波段具有低發(fā)射率，且在1.06 μm和1.54 μm波長處具有較低反射率的一維異質(zhì)結(jié)構(gòu)光子晶體，通過理論分析發(fā)現(xiàn)，所設(shè)計的光子晶體能夠較好地滿足在上述波段和波長處的偽裝要求。該研究結(jié)果對于下一步的實際應(yīng)用具有借鑒和參考價值。

2 傳輸矩陣法

圖1所示為理想光子晶體的反射率曲線。

圖1 理想光子晶體反射率

光子晶體特有的周期性結(jié)構(gòu)對入射電磁波有調(diào)制作用。電磁波在其中的傳播可以等效為在若干個基本單元中傳播的疊加，而每個基本單元都可以看作一個普通的雙端口器件。它的傳播特性采用一個2×2復矩陣描述，整個結(jié)構(gòu)的傳播特性也可以采用一個2×2的復矩陣進行表示，該矩陣為各單元矩陣的乘積。該方法稱為傳輸矩陣法(Transfer matrix method,TMM)[4]。

當橫電(或磁)波以角度θ0入射時，光子晶體中第k層材料的折射角θk滿足Snell折射定律，即

nksinθk=nk-1sinθk-1=…n0sinθ0，

(1)

其特征矩陣Mk可表示為[5]

(2)

其中

δk=2πλ-1nkdkcosθk，

式中，ηk和dk分別為第k層材料的折射率和幾何厚度。此時，光子晶體總的特征矩陣可以表示為

(3)

由此可以得到反射率R的表達式為

(4)

透射率T的表達式為

(5)

3 結(jié)構(gòu)設(shè)計

光子晶體的帶隙主要有3個影響因素：光子晶體空間結(jié)構(gòu)分布、介電常數(shù)比值以及光子晶體的幾何結(jié)構(gòu)[6]。一維光子晶體結(jié)構(gòu)在空間一個方向上延伸，結(jié)構(gòu)單一，因此空間結(jié)構(gòu)分布對一維光子晶體的帶隙影響不大。光子晶體的幾何機構(gòu)對帶隙的影響與具體結(jié)構(gòu)有關(guān)，例如在一維梳狀光子晶體[7]中，禁帶寬度隨接枝長度的變化而變化，但其實質(zhì)仍然是介質(zhì)折射率之比。由于一維光子晶體是層狀周期性結(jié)構(gòu)，因此，人們主要基于此類結(jié)構(gòu)中不同介質(zhì)介電常數(shù)之比對帶隙的影響進行光子晶體的帶隙設(shè)計。一般來說，兩種介質(zhì)的折射率比值越大，其對應(yīng)的光子晶體帶隙就越寬[8]。但由于常見材料的高低折射率比值的限制，通常制備的光子晶體帶隙很難覆蓋近中紅外波段。因此，必須對光子晶體的帶隙進行拓展。文獻[9-10]曾提出通過隨機或者漸變地改變一維光子晶體每個單元內(nèi)介質(zhì)層的厚度的方法來拓寬帶隙。雖然該方法理論上可行，但其在實際設(shè)計過程中的可行性欠佳。另一種方法是利用幾個一維光子晶體的組合[11-12]來獲得較大的相對帶寬。通過調(diào)整各個光子晶體的周期數(shù)，使每個晶體的禁帶互相銜接，從而擴展了相對帶寬。該方法獲得的光子晶體通常被稱為異質(zhì)結(jié)構(gòu)光子晶體。

評價禁帶寬度一般有兩種方式：一種是計算禁帶的絕對帶寬，即禁帶上下邊界的差值；另一種則是計算絕對帶寬和禁帶中心的比值，即相對帶寬D(Δω/ω0)。兩個一維光子晶體合成后的相對帶寬D0可表示為[13]

(6)

兩個光子晶體相對帶寬的簡單代數(shù)相加為

(7)

其中，a1、b1、a2、b2分別為歸一化的兩個晶體禁帶的下邊界和上邊界，以(2πc/Λ)為單位，Λ為周期。這種通過若干個晶體疊加的方式能夠擴展一維光子晶體的相對帶寬，且與疊加順序無關(guān)[14]。

圖2為本文設(shè)計的光子晶體的結(jié)構(gòu)示意圖。我們選用碲(Te)和氟化鎂(MgF2)分別作為高低折射率材料[15]，充分考慮這兩種介質(zhì)材料在近中紅外波段實部和虛部的影響，設(shè)計了SUB|(AB)5|(AB)2(ABC)AB|AIR異質(zhì)結(jié)構(gòu)的一維光子晶體。其中SUB是基底，AIR是空氣，A為高折射率材料Te，B為低折射率材料MgF2。前5個周期為PC1，A的厚度為161.73 nm，B的厚度為558.6 nm；后4個周期為PC2，A的厚度為80.56 nm，B的厚度為281.96 nm(第8周期A的厚度為10.93 nm，B的厚度為687.78 nm)。C為摻雜材料鍺(Ge)，厚度為21.89 nm。其中第8周期厚度不同的A和B以及摻雜材料C是對PC2引入的兩種缺陷。由于PC1和PC2是相互獨立的兩個對稱結(jié)構(gòu)，因此在PC2中引入的缺陷不會對PC1的帶隙產(chǎn)生影響[16]，而且可以在PC2的帶隙中產(chǎn)生兩個缺陷態(tài)。由缺陷模的波函數(shù)的迭加[17]可知，兩缺陷層間的光學厚度越小，其間的相互作用就越強[18]。由于文中設(shè)計的這兩種缺陷的光學距離很小，因此，微小的光學厚度以及相對位置的變化均會引起缺陷位置和透過率的急劇變化。所以，在制備過程中對膜層厚度和相對位置的控制提出了較高要求，應(yīng)盡量減小膜層位置及厚度對光子晶體特性的影響。通過理論計算，獲得該結(jié)構(gòu)光子晶體的光譜反射率和吸收率曲線分別如圖3和圖4所示。

圖2 光子晶體結(jié)構(gòu)圖

圖3 光子晶體的光譜反射率

圖4 光子晶體的光譜吸收率

4 結(jié)果分析

從圖3可以看出，該結(jié)構(gòu)光子晶體在1～5 μm波段出現(xiàn)了高反帶，在1.81～2.47 μm波段內(nèi)的反射率達到95.07%以上，在2.47～5 μm波段內(nèi)的反射率為1。覆蓋在飛行器表面的這種具有高反射帶光子晶體，可以使飛行器產(chǎn)生的處于高反射帶內(nèi)的大部分紅外輻射反射回去，從而有效抑制飛行器高溫表面的紅外輻射；同時，在1.06 μm和1.54 μm處的反射率分別為1.21%和1.79%，形成了兩個明顯的反射谷，這兩個常見軍用激光器波長處的低反射可大大降低對方探測激光的回波功率，達到激光隱身的目的，從而可使飛行器避免來自激光器的威脅。結(jié)合圖4的光譜吸收率曲線可知，在1.06 μm和1.54 μm處的吸收率分別為98.79%和96.48%。從1.54 μm之后，吸收率急劇降低。吸收率在1.81 μm處為4.72%；在1.81～2.47 μm范圍內(nèi)，吸收率保持在4.72%以下；在2.47～5 μm范圍內(nèi)，吸收率為0。說明這兩個激光波長處的低反射率是由光子晶體的高吸收造成的，而在其他波段，光子晶體的吸收率很低。在圖3的反射率曲線中，反射率從5 μm之后開始下降，5.16 μm處的反射率為95.04%，5.19 μm處下降到90.30%，之后開始低于90%，在5.28 μm處達到最低，為37.72%，之后雖然稍有上升，但是最高僅為90.03%，這之后又有更大幅度的下降?？傮w來說，反射率一直維持在一個較低的水平上。而由圖4可以看出，吸收率在這個波段也一直維持在極低的水平上。這些波段的低反射和低吸收必然對應(yīng)著高透射，這種特性使光子晶體在高反射帶和1.06 μm、1.54 μm高吸收所積聚的能量透射出去，起到散熱的作用，避免溫度的升高對光子晶體的壽命以及隱身效能產(chǎn)生影響。這部分能量大部分處在紅外大氣窗口之外，由于大氣中的水蒸氣、二氧化碳和臭氧的強烈吸收而迅速衰減，從而不能夠被探測方所接收?？梢钥闯觯疚脑O(shè)計的光子晶體既能夠抑制近中紅外波段的紅外輻射，又能吸收1.06 μm和1.54 μm的入射激光，同時還能保持其他波段的高透射，這種對反射率、吸收率和透射率的人為干預(yù)起到了光譜轉(zhuǎn)換的作用，在多光譜兼容偽裝方面具有一定的意義。

在設(shè)計中，缺陷膜層厚度變化對缺陷膜的影響十分明顯，且這種影響不是獨立的。保持其他層厚度不變，當?shù)?周期(缺陷周期)的Te層厚度增加10 nm時，1.06 μm激光的反射率變?yōu)?.83%，而1.54 μm激光的反射率增加到63.48%；當減少10 nm時，1.06 μm激光的反射率增加到4.68%，而1.54 μm激光的反射率激增至86.81%。同樣保持其他層厚度不變，改變摻雜層Ge的厚度，當其增加10 nm時，1.06 μm激光的反射率增至14.23%，1.54 μm激光的反射率增加到9.69%；當減少10 nm時，1.06 μm激光的反射率變?yōu)?7.22%，而1.54 μm激光的反射率增至8.45%。繼續(xù)增大改變量，增加20 nm時，1.06 μm激光反射率增加到37.30%，1.54 μm激光反射率增加到24.44%；當減少20 nm時，1.06 μm激光反射率變?yōu)?1.48%，而1.54 μm激光反射率增加到25.80%，這時該光子晶體已經(jīng)失去對激光的偽裝作用。由此可見，膜層厚度的控制是這種光子晶體設(shè)計的關(guān)鍵所在。

5 結(jié) 論

本文根據(jù)光子晶體的帶隙和光子局域特性，設(shè)計了一種可以實現(xiàn)近中紅外和1.06 μm與1.54 μm激光兼容隱身的光子晶體材料。首先利用薄膜的傳輸矩陣法，對電磁波在光子晶體中的傳播進行了理論分析，隨后利用異質(zhì)結(jié)構(gòu)拓展了光子晶體的帶隙，使之在1.81～2.47 μm波段的反射率達到95.07%以上，在2.47～5 μm波段的反射率達到1。針對PC2引入兩種缺陷，使得1.06 μm和1.54 μm處的反射率降低到1.21%和1.79%。 最后對設(shè)計過程當中應(yīng)該注意的問題進行了定量分析，說明了精確控制膜層物理厚度的重要性。文中設(shè)計光子晶體時運用的異質(zhì)結(jié)構(gòu)理論、摻雜理論等對實現(xiàn)其他波段的紅外兼容偽裝具有一定的借鑒意義。

[1] 張建奇，方小平. 紅外物理 [M]. 西安：西安電子科技大學出版社， 2007. ZHANG J Q, FANG X P.InfraredPhysics[M]. Xi’an: Xi’an Electronic Science and Technology University Press, 2007. (in Chinese)

[2] FINK Y， WINN J N， FAN S,etal.. A dielectric omnidirectional reflector [J].Science, 1998, 28(2):1679-1682.

[3] WIERSMA D S, BARTOLINI P, LAGENDIJK A,etal.. Localization of light in a disordered medium [J].Lett.Nat., 1997, 390:18-25

[4] 王丹，王勇，韓琳，等. 近紅外光區(qū)一維光子晶體高反射鏡的設(shè)計 [J]. 揚州大學學報(自然科學版), 2012, 15(1):29. WANG D， WANG Y, HAN L,etal..Design on a one-dimensional crystal high reflector in near-infrared region [J].J.YangzhouUniv.(Nat.Sci.Ed.), 2012, 15(1):29. (in Chinese)

[5] ZHAO X K, ZHAO Q W, ZHANG Q H,etal.. Numerical study on laser and infrared compatible stealth by forming “hole-digging spectrum” of doped photonic crystal [J].SPIE, 2015, 7847:784731-7.

[6] 趙芬芬. 基于缺陷態(tài)光子晶體的多通道濾波器設(shè)計與分析 [D]. 秦皇島：燕山大學， 2013. ZHAO F F.TheDesignandAnalysisofMulti-channelFilterBasedonDefectPhotonicCrystal[D]. Qinhuangdao: Yanshan University, 2013. (in Chinese)

[7] 翦知漸,彭景翠,張高明,等. 一維類梳狀波導光子晶體的頻率帶隙寬度 [J]. 湖南大學學報(自然科學版)， 2004，31(2):103-106. JIAN Z J, PENG J C, ZHNAG G M,etal.. The width of frequency band gap in one-dimensional comb-like waveguide photonic crystals [J].J.HunanUniv. (Nat.Sci.Ed.), 2004, 31(2):103-106 . (in Chinese)

[8] 丁明耀. 一維三元光子晶體的帶隙擴展方法研究 [D]. 長沙：湖南大學， 2012 DING M Y.DesignofBroadbandOmnidirectionalBandGapinTheOne-dimensionalTernaryPhotonicCrystals[D]. Changsha: Hunan University, 2012. (in Chinese)

[9] YABLONOVITCH E. Engineered omnidirectional external-reflectivity spectra from one-dimensional layered interference filters [J].Opt.Lett.,1998 , 23(21):1648-1649.

[10] LI H Q , CHEN H , QIU X J. Band -gap extension of disordered 1D binary photonic crystals [J].Physical, 2000, 279(1-3):164-167.

[11] MIR A, AKJOUJ A, BOUDOUTI E H,etal.. Large omnidirectional band gaps and selective transmission in one-dimensional multiplayer photonic structures [J].Vacuum， 2001, 63(1-2):197-203.

[12] WANG X, HU X H, LI Y Z,etal.. Enlargement of omnidirectional total reflection frequency range in one-dimensional photonic crystals by using photonic heterostructures [J].Appl.Phys.Lett., 2002, 80(23):4291-4293.

[13] 黃弼勤，顧培夫. 一維光子晶體禁帶的展寬 [J]. 光學學報， 2003, 23(12):1498-1499. HUANG B Q, GU P F. Extension of one-dimensional photonic crystal’s band gap [J].ActaOpticaSinica, 2003, 23(12):1498-1499. (in Chinese)

[14] 王充，彭同江，段濤,等. 異質(zhì)結(jié)構(gòu)光子晶體的制備與帶隙特性研究 [J]. 人工晶體學報， 2010， 39(2):475-480. WANG C, PENG T J, DUAN T,etal.. Preparation and band-gap properties of heterostructure photonic crystal [J].J.Synth.Cryst., 2010， 39(2):475-480. (in Chinese)

[15] PALIK E D.HandbookofOpticalConstantsofSolids[M]. Maryland: Academic Press, 1985.

[16] 田國勛，吳永剛，王占山,等. 一維光子晶體的雙通道位置設(shè)計及調(diào)整 [J]. 光學學報， 2005, 25(5):662-664. TIAN G X, WU Y G, WANG Z S,etal.. Design and modification of double channel of one-dimensional photonic crystal [J].ActaOpticaSinica, 2005, 25(5):662-664. (in Chinese)

[17] 徐海燕,陳鶴鳴. 大入射角一維磁光子晶體寬帶光隔離器 [J]. 光子學報, 2015, 44(10):1023004. XU H Y, CHEN H M. Broad-band optical isolator with wide-angle based on one-dimensional magneto-optical crystals [J].ActaPhoton.Sinica, 2015, 44(10):1023004. (in Chinese)

[18] 李文勝,付艷華,張琴,等. 含雙缺陷的一維準周期光子晶體的缺陷模特性 [J]. 半導體光電， 2010, 31(2):241-243. LI W S, FU Y H, ZHANG Q,etal.. Characteristics of defect mode of quasi-periodic 1-D photonic crystal containing double defects [J].Semicond.Optoelectron., 2010, 31(2):241-243. (in Chinese).

張繼魁(1992-)，男，河南周口人，碩士研究生，2015年于合肥電子工程學院獲得學士學位，主要從事光子晶體方面的研究。

E-mail: m17718130921@163.com

E時家明(1966-)，男，安徽巢湖人，教授，博士生導師，1996年于中國科學院等離子體物理研究所獲得博士學位，主要從事隱身與反隱身技術(shù)的研究。

E-mail: shijiaming66@163.com

Research on Compatible Stealth Photonic Crystal Against Near/Middle Infrared and 1.06 μm and 1.54 μm Lasers

ZHANG Ji-kui, SHI Jia-ming*, MIAO Lei, WANG Qi-chao, ZHAO Da-peng, ZENG Jie

(StateKeyLaboratoryofPulsedPowerLaserTechnology,ElectronicEngineeringInstitute,Hefei230037,China)*CorrespondingAuthor,E-mail:shujiaming66@163.com

A compatible stealth material against near/middle infrared as well as 1.06 μm and 1.54 μm lasers, based on the structure of one-dimensional photonic crystal (PC) was designed to achieve the compatible stealth of high temperature part of aircrafts against infrared and lasers. Based on the transfer matrix method (TMM) of thin-film optical theory and heterogeneous structure theory, the forbidden band can be broadened to cover the near/middle infrared wave band. Then, two defects were added to the periodical structure of the PC according to the doping theory. The results show that there are two defect modes that located in the wavelengths of 1.06 μm and 1.54 μm in the band gap of 1-5 μm, respectively, with the corresponding spectral reflectance of 1.21% and 1.79%. This PC with “hole-digging spectrum” can realize the compatible stealth against near/middle infrared as well as 1.06 μm and 1.54 μm lasers.

photonic crystal; near/middle infrared; laser; compatible stealth

1000-7032(2016)09-1130-05

2016-05-09；

2016-06-27

國家高技術(shù)研究發(fā)展計劃(863計劃)(2015AA0392)資助項目

TN29

A

10.3788/fgxb20163709.1130