李曉楠 周璐 趙國忠?

1) (首都師范大學物理系,北京 100048)

2) (北京市成像理論與技術高精尖創(chuàng)新中心,北京 100048)

3) (太赫茲光電子學教育部重點實驗室,北京 100048)

具有螺旋波前的電磁波是攜帶軌道角動量的渦旋波束,渦旋波束存在的相位奇點,使其在微粒操控和通訊等領域有特殊的應用.本文提出了一種基于反射型超表面的太赫茲寬帶渦旋波束產生器,該器件由超表面-電介質-金屬三層結構構成,頂層為兩個正交I形金屬結構單元組成的超表面,中間層是聚酰亞胺介質,最底層為金屬.通過對超表面單元結構參數(shù)的優(yōu)化設計,可以實現(xiàn)在不同旋轉角度下反射波振幅高達90%以上,同時反射波相位隨旋轉角線性變化的目的.進一步利用這些單元結構按照Pancharatnam-Berry相位原理進行超表面布陣,在0.8—1.4 THz頻率范圍內,可以將圓偏振太赫茲波束轉換為具有軌道角動量的渦旋波束,這一器件的工作帶寬相對較寬,結構簡單,轉換效率高,在太赫茲渦旋波束產生方面具有潛在的應用價值.

1 引 言

太赫茲波具有穿透性、吸水性、高分辨率等諸多優(yōu)越特性,在光譜[1-3]、成像[4-6]、無損檢測[7,8]等領域已經(jīng)展示了重要的應用前景.攜帶軌道角動量的渦旋光束,具有由exp(ilΦ)描述的螺旋相位結構,其中Φ是方位角,l 是拓撲荷數(shù),它在不增加帶寬的情況下通過對光通訊信息進行編碼,可以極大地擴展通信容量,在光通信[9]領域具有廣闊的應用前景.傳統(tǒng)上人們利用螺旋相位板[10]、全息衍射光柵[11]以及天線陣列[12,13]等方法來產生渦旋波束,然而這些方法存在體積大,不利于器件集成化、成本高等問題.作為一種人工合成材料,超表面具有自然材料所不具備的性質,因而受到廣泛關注.人們通過在金屬表面設計微結構來實現(xiàn)對電磁波的相位和振幅的調制.近年來應用超表面實現(xiàn)的功能器件有偏振轉換器[14-16]、高方向輻射天線[17]、高增益透鏡陣列[18,19]以及渦旋光束產生器[20,21]等.

利用超表面產生渦旋波束不僅解決了體積大、不利于集成的問題,還大大降低了制造成本.近年來,在可見光和微波波段基于超表面產生渦旋波束的研究較多.2011年,哈佛大學Yu等[22]在光學薄界面上設計V字形天線單元,通過選擇八個相位間隔45°的V字形單元結構排布超表面,實現(xiàn)了在紅外波段產生拓撲荷數(shù)為1的渦旋光.2014年,Zhang等[23]在微波頻段內基于開口諧振環(huán)提出了單層透射超表面,圓偏振光入射超表面產生了拓撲荷數(shù)為1的渦旋光,但是透過率僅有24.7%.2017年,Lou等[24]基于三層超表面結構在10.1—10.9 GHz產生了渦旋光,但是存在制備復雜等問題.

目前為止,關于渦旋波束的產生,普遍存在效率不高、制備困難等問題,且大部分工作是在微波、紅外頻段,在太赫茲頻段產生渦旋波束的研究甚少.本文提出了一種新型反射式超表面結構,實現(xiàn)了在太赫茲頻段產生高效率渦旋光的目的.

2 理論分析

為了實現(xiàn)利用超表面產生渦旋光的目的,需要設計能夠獨立地控制入射和反射波的單元結構.根據(jù)Pancharatnam-Berry (P-B)相位原理,超表面上可以進行單元結構布陣,各個單元應該能夠獲得振幅相近、相位可以獨立調節(jié)的反射波.我們提出了一個由正交I形金屬結構單元組成的結構,可以滿足上述控制反射波振幅和相位的設計要求.

當一束平面波垂直照射到反射型超表面單元結構,單元結構以其中心為旋轉中心,波束傳播方向為旋轉軸,逆時針旋轉α角時,如圖1所示,反射矩陣可以寫作[25]：

其中 rxx表示沿x 軸正方向的線偏振光入射時,沿x 軸負方向出射的反射系數(shù),同理 rxy,ryx,ryy有類似的意義.

對于圓偏振光入射到超表面,反射矩陣表示為

其中參數(shù)α為單元結構相對于其中心的轉角,LR代表入射光的偏振態(tài),L表示左旋偏振,R表示右旋偏振,RαLR表示圓偏振光入射時的反射矩陣.

從(3)式可以看出,只有 rLL,rRR攜帶軌道角動量,當單元結構滿足 |rxx|≈|ryy|≈ 1 或者|rxy|≈|ryx|≈ 0且相位差為π時,兩個交叉極化之間不存在相位差,因此當單元結構旋轉角α時,相位φ=±2α,“+”和“—”表示右旋和左旋入射光.

基于上述理論可知,在設計P-B超表面時,通過旋轉單元結構找出具有不同相位的單元,使得單元結構具有相同且近似為1的反射振幅,同時具有180°的相位差,可實現(xiàn)在圓偏振光入射時,反射超表面可以使入射光成為攜帶軌道角動量的渦旋光.但是,目前渦旋光產生器件存在帶寬窄的問題,為了增加超表面的工作帶寬,需要確保在不同頻率處具有相同的斜率：

其中 ?φ(R,R) 表示右旋入射光的反射相位變化,fi和 fj表示任意兩處不同的頻率,為了保證超表面有合適的工作帶寬,需要保證 (fj-fi) 足夠大.

3 單元結構設計

在超表面結構設計過程中,首先考慮的主要問題為如何設計單元結構、單元結構參數(shù)如何確定、襯底如何選擇等.我們設計了一個由超表面-電介質-金屬三層結構組成的反射型超表面器件,如圖1所示.超表面由三層組成.其中頂層為正交I型金屬結構,厚度為d1；中間層為介質層,材料為聚酰亞胺,厚度為d2；底層為金膜,厚度為d1.考慮實驗室材料、各個仿真結果以及樣品制備,最終得出單元結構具體參數(shù)如下：周期p = 90 μm,金屬線寬w = 3 μm,l1= 60 μm,l2= 28 μm,s = 27 μm,金屬層厚度d1= 0.2 μm,襯底厚度d2= 35 μm.入射光為偏振方向沿x 或y 軸方向的線偏光,波矢沿z 軸方向.

圖1 超表面單元結構示意圖 (a)頂視圖；(b)側視圖Fig.1.Schematic diagram of the unit cell of metasurface：(a) Top view；(b) side view.

4 單元結構仿真結果與分析

利用CST microwave studio (CST)軟件,以周期性邊界條件對單元結構反射振幅和反射相位進行仿真,仿真頻率范圍設為0.5—2.0 THz,結果如圖2.在線性x 偏振和線性y 偏振入射下,在0.8—1.4 THz頻率范圍內共極化反射幅值均大于0.9且相位差接近180°,在0.93 THz時共極化相位差剛好為180°且反射幅值近似為1.因此可以通過改變轉角實現(xiàn)對反射波振幅和相位的調制.

圖2 線偏振光入射單元結構產生的反射相位和振幅譜 (a)相位譜；(b)振幅譜Fig.2.Reflected phase and amplitude spectrum produced by linearly polarized light incident unit call：(a) Phase；(b) amplitude.

以0.85 THz為例,對單元結構反射振幅和相位進行數(shù)值仿真,可得出不同轉角與相位及振幅的一一對應關系,如圖3.不同轉角下單元結構的反射振幅近似相等,轉角與相位呈現(xiàn)線性關系,與上述P-B相位原理符合得很好.因此可得到不同相位值對應的不同轉角的單元結構.

入射光為圓偏振光時,對參數(shù)α進行掃描,掃描范圍為0°—157.5°,掃描步長設置為22.5°,仿真結果如圖4所示.

圖3 仿真得到的不同轉角的相位和反射振幅Fig.3.The phase and reflection amplitude of different corners simulated.

圖4 不同轉角單元結構的反射相位和振幅譜 (a) LCP入射的相位譜；(b) RCP入射的相位譜；(c) LCP入射的振幅譜；(d) RCP入射的振幅譜Fig.4.Reflective phase and amplitude spectra of unit cell structure under different rotation angle：(a) Phase spectra at LCP incident；(b) phase spectra at RCP incident；(c) amplitude spectra at LCP incident；(d) amplitude spectra at RCP incident.

圖4(a)和圖4(b)分別為左旋圓偏光(LCP)和右旋圓偏光(RCP)入射到超表面單元結構時,轉角α在0°—157.5°范圍內以22.5°變化時得到的相位譜.從圖4(a)和圖4(b)可以看出,相位以二倍轉角(45°)間隔增大(減小),不同轉角對應的相位覆蓋了2π.圖4(c)和圖4(d)分別為左旋圓偏光和右旋圓偏光入射到超表面單元結構時,轉角α在0°—157.5°范圍內以22.5°的變化改變時得到的振幅譜,可以看出,反射振幅幾乎相同且高達0.9,滿足控制反射波振幅和相位的設計要求.

5 P-B超表面設計及仿真結果分析

基于文中提出的反射型超表面單元,可以通過相位分布設計產生不同拓撲荷數(shù)的渦旋波束.以超表面中心為原點,超表面每個位置(x ,y)的相位需要滿足

其中 φ(x,y) 表示超表面 (x,y) 位置需要滿足的相位,l 表示拓撲荷數(shù).通過改變l 的值可以產生不同拓撲荷數(shù)的渦旋波束.

本文設計了拓撲荷數(shù)l = 1和l = 2的反射型超表面.根據(jù)(4)式確定的相位-位置關系以及單元結構相位是轉角α的2倍關系,排布了拓撲荷數(shù)為1和2的超表面結構,如圖5所示.

對l = 1和l = 2的兩種超表面反射波束特性進行仿真.不失一般性,假設入射波束為左旋圓偏振波,半徑r1,r2,r3分別設為120,200,300 μm,其中,半徑的選擇考慮到各個單元結構之間不能夠重疊,以及實驗測量中的尺寸、加工難度的因素.以0.93 THz為例,圖6給出了反射圓偏振波的振幅和相位分布.

圖5 兩種用于產生拓撲荷數(shù)分別為 (a) l = 1和(b) l = 2的渦旋波束超表面結構Fig.5.Two kinds of metasurface structures for generating the vortex beam with topological charge numbers (a) l = 1 and (b) l = 2.

圖6 超表面產生l = 1和l = 2的渦旋波束反射振幅和相位分布 LCP入射l = 1超表面的(a)振幅分布和(b)相位分布；RCP入射l = 2超表面的(c)振幅分布和(d)相位分布Fig.6.Reflective amplitude and phase distributions of vortex beams with l = 1 and l = 2 generated by metasurface.LCP incident l = 1 metasurface：(a) amplitude distribution and (b) phase distribution；RCP incident l = 2 metasurface：(c) amplitude distribution and (d) phase distribution.

圖6(a)和圖6(b)是0.93 THz處左旋圓偏振波垂直入射到拓撲荷數(shù)為1的反射型超表面上,反射光場距離超表面500 μm處的反射振幅和相位分布.從圖中可以看出振幅呈現(xiàn)中心為暗環(huán)且中心強度始終為零的分布,相位圍繞一周改變2π.圖6(c)和圖6(d)是0.93 THz的右旋圓偏振波垂直入射到拓撲荷數(shù)為2的反射型超表面上,反射光場距離超表面500 μm處的反射振幅和相位分布.可以看出振幅分布呈現(xiàn)中心場強為零,外環(huán)場強相對較大,相位覆蓋4π.表明了在圓偏振波束入射到超表面可以產生拓撲荷數(shù)為1和2的渦旋波束.

圖7 不同頻率下l = 1和l = 3超表面產生的反射渦旋波束振幅、相位分布 l = 1超表面(a) 0.8 THz頻率下振幅分布,(b) 0.8 THz頻率下相位分布,(c) 1.4 THz頻率下振幅分布,(d) 1.4 THz頻率下相位分布；l = 3超表面(e) 1.1 THz頻率下振幅分布,(f) 1.1 THz頻率下相位分布Fig.7.The amplitude and phase distribution of reflective vortex beam generated by the LCP incident l = 1 and l = 3 metasurface at different frequencies.l = 1：(a) amplitude distribution at 0.8 THz,(b) phase distribution at 0.8 THz,(c) amplitude distribution at 1.4 THz,(d) phase distribution at 1.4 THz.l = 3：(e) amplitude distribution at 1.1 THz,(f) phase distribution at 1.1 THz.

從圖6(a)和圖6(c)可以看出拓撲荷數(shù)為1的渦旋光束的中心暗環(huán)半徑相對于拓撲荷數(shù)為2的渦旋光中心暗環(huán)半徑較小,這是因為隨著拓撲荷數(shù)的增大,超表面相鄰的相位梯度也隨之增大導致的.

在0.8—1.4 THz頻段范圍內,圓偏光入射到超表面時均可產生渦旋光.圖7給出了在0.8 THz和1.4 THz兩個頻率下,左旋圓偏振光入射在l = 1的超表面時,反射波束在傳播一定距離處的相位分布.

圖7(a)和圖7(c)分別為左旋圓偏光在0.8,1.4 THz下入射到l = 1的超表面上產生的反射光場距離超表面700 μm處的振幅分布,可以看出,振幅分布相對均勻,呈現(xiàn)中間暗環(huán),外環(huán)場強相對較大的甜甜圈分布,圖7(b)和圖7(d)分別為左旋圓偏光在0.8,1.4 THz下入射到l = 1的超表面上產生的反射光場距離超表面700 μm處的相位分布,可以看出相位變化了2π,綜上符合渦旋光的相位性質,并且反射率達到86%.

此外,為證明設計的超表面可以產生更高階的渦旋光束,設計并仿真了拓撲荷數(shù)為3的超表面振幅和相位圖像.以1.1 THz頻率下為例,圖7(e)和圖7(f)給出了反射場距離超表面300 μm處的反射振幅和相位分布.從圖中可以看出振幅呈現(xiàn)中間暗環(huán)的甜甜圈形狀,相位覆蓋了6π,仍然符合渦旋光束的特性.

結合圖6和圖7的仿真結果,在0.8,0.93,1.4 THz頻率下,圓偏振光垂直入射到l = 1和l = 2的超表面時,反射波束具有螺旋相位,可產生中心強度為零的渦旋光.證明提出的反射型超表面可以在0.8—1.4 THz頻率范圍內產生渦旋光,并且可以產生更高階的渦旋光束.

6 結 論

基于相位突變的結構單元,利用P-B相位原理設計了一個簡單而高效的反射型超表面,該超表面由正交I型結構-介質-金屬三層結構組成.仿真結果表明,我們設計的超表面結構可以在0.8—1.4 THz頻段圓偏振光正入射條件下,產生攜帶軌道角動量的渦旋波束,反射率達到86%,最高可達92%.我們設計的超表面具有效率高、寬帶并且結構簡單等優(yōu)點,對基于超表面的太赫茲渦旋波束產生具有參考價值.