• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    超穎表面原理與研究進展

    2017-10-23 10:27:20李天佑黃玲玲王涌天
    中國光學 2017年5期
    關鍵詞:全息圖偏振光渦旋

    李天佑,黃玲玲,王涌天

    (北京理工大學 光電學院,北京 100081)

    超穎表面原理與研究進展

    李天佑,黃玲玲*,王涌天

    (北京理工大學 光電學院,北京 100081)

    超穎表面是一種特殊的二維亞波長陣列結(jié)構,具有很強的光場調(diào)控能力,兼有超薄、低損耗、易加工等優(yōu)勢,表現(xiàn)出廣闊應用前景,近些年吸引了廣泛的研究興趣。本文將綜述超穎表面的原理與研究進展,簡要分析已報道的超穎表面類型,并重點介紹其在廣義折反射、偏振轉(zhuǎn)換、光束軌道角動量操控以及計算全息中的應用。超穎表面設計靈活、功能強大,有望在諸多應用中取代傳統(tǒng)光電器件。未來基于超穎表面的新型、寬帶、高效率、多功能以及主動可調(diào)功能器件等將是其重要的發(fā)展方向。

    超穎表面;偏振轉(zhuǎn)換;光束軌道角動量;計算全息

    1 引 言

    光學元器件和設備通過以指定方式調(diào)控光的振幅、相位以及偏振,從而得以實現(xiàn)對光束傳播方向與波前的操縱,并在顯示、傳感與通信等應用中賦予所需的信息。傳統(tǒng)的光學器件往往通過光在特定介質(zhì)中累積傳播的光程來達到對光的調(diào)控[1-2],具有較大的體積,難以適用于未來集成化、微型化器件發(fā)展趨勢。

    超穎材料(Metamaterial)是一種由亞波長周期微納結(jié)構構成,具有自然界材料所不具備的超常物理特性的人造電磁材料[3]。因其獨特的特性,超穎材料在過去的二十年間引起了廣泛的研究[4-12],并在超分辨成像[13-14],光學隱形[15-16]等領域有著潛在的應用前景。傳統(tǒng)的超穎材料可視為等效介質(zhì)(亦稱為體超穎材料,Bulk Metamaterial),對光場的調(diào)控仍依賴于光在材料中的傳播以及相互作用的積累。超穎材料通常使用金屬材料并利用了結(jié)構的諧振響應,具有難以避免的高損耗與色散特性;同時亞波長三維結(jié)構不易加工制造,這些因素限制了超穎材料的進一步研究與應用。

    為了克服超穎材料的諸多缺點,人們想到了“降維”的辦法,將研究重點轉(zhuǎn)向平面二維結(jié)構。超穎表面(Metasurface)是一種二維平面超穎材料,通常由亞波長尺寸的周期,準周期或者隨機的單元構成,其厚度遠小于波長,可視為二維表面[17-18]。超穎表面利用每個單元結(jié)構對入射光場的強烈響應來改變局部光場的振幅和相位,從而以亞波長尺度對光場振幅與相位進行調(diào)制,進而實現(xiàn)對近場與遠場的調(diào)控。這種遠小于入射波長的厚度使得超穎表面不再依賴光在介質(zhì)中傳播的光程累積,因此介電常數(shù)、磁導率、折射率等參數(shù)不再被用來描述超穎表面,而往往是用表面阻抗、振幅、相位、偏振等參量來描述其物理性質(zhì)。超穎表面在具有很強的光場調(diào)控能力的同時,兼有超薄、低損耗、易加工等優(yōu)勢,在奇異光學[19]、光束整形、偏振調(diào)控[20]、光學全息技術[21]、光束軌道角動量操控[19]等方面具有廣闊的應用前景。

    超穎表面本質(zhì)上是一種有著不同幾何結(jié)構以及亞波長空間間隔的光學天線陣列,通過設計每個納米天線輻射場的相位延遲,利用惠更斯原理,能夠?qū)崿F(xiàn)對出射波前的任意調(diào)控。這與微波波段的頻率選擇表面(Frequency-Selective Surface,F(xiàn)SS)有著本質(zhì)的不同,頻率選擇表面沒有空間可變的電磁響應。相比傳統(tǒng)光學元件以及體超穎表面,超穎表面具有一些顯著優(yōu)點[22]:(1)超薄化。超穎表面能夠在遠小于入射波長的尺度上實現(xiàn)對入射光場振幅、偏振與相位的調(diào)控,顯著減小了光在介質(zhì)中傳播的損耗。(2)亞波長像素。超穎表面能夠以亞波長尺度的像素對光場進行調(diào)控,能有效避免在傳統(tǒng)衍射光學元件中多級衍射的串擾,在微型集成光電系統(tǒng)中有著廣闊的應用前景。(3)寬帶特性。通過特定的設計,超穎表面能夠在寬波段對光場進行同樣的調(diào)制,具有弱色散或者無色散的特性。上述優(yōu)點再加上超穎表面在設計靈活性以及低成本制造方面的優(yōu)勢使得其有望成為新一代超薄、小型光場調(diào)控器件。

    本文將從以下兩方面對超穎表面的研究進展進行了概述,并結(jié)合作者對超穎表面的認識對超穎表面的發(fā)展提出了展望。

    (1)超穎表面的基本分類。按照光場透反射空間、陣列結(jié)構類型以及相位調(diào)控機理,超穎表面可分為不同大類。

    (2)超穎表面的功能應用研究現(xiàn)狀,包括利用超穎表面的異常折射/反射現(xiàn)象、基于超穎表面的任意波前調(diào)控以及基于超穎表面的偏振轉(zhuǎn)換器件研究。

    2 超穎表面調(diào)控機理研究現(xiàn)狀

    2.1基于光學天線的諧振響應

    光學天線[22]通常具有多種多樣的結(jié)構,包括金屬/介質(zhì)顆粒,金屬薄膜上的小孔以及上述二者的多層疊加。光學天線能夠利用其局域表面等離激元諧振響應調(diào)控天線輻射場的相位,從而在亞波長尺度實現(xiàn)對光場的操控。設計特定結(jié)構的光學天線是此類超穎表面的核心思想。

    2011年,哈佛大學的Yu等人[19]首次使用V型天線實現(xiàn)了對表面相位的不連續(xù)調(diào)控。V型天線由兩個長度均為h的臂構成,他們之間的夾角為Δ,如圖1(a)所示。V型天線具有所謂的雙頻諧振特性,即該天線可以支持“對稱”和“反對稱”模式。對稱模式近似成長度為h的單臂天線的情形,因此其諧振波長發(fā)生在h≈λeff/2處,其中λeff為等效波長;反對稱模式近似為長度為2h的直臂天線,其諧振波長發(fā)生在2h≈λeff/2處。在兩個諧振波長附近,天線對相位的改變非常顯著。兩種諧振模式能夠被入射光獨立激發(fā),并且一般會導致散射場與入射場具有不同的偏振態(tài)。通過設計天線的結(jié)構參數(shù),可以使得在線偏振光入射時,其正交偏振態(tài)的出射光擁有對相位和振幅的調(diào)控能力,相位的調(diào)制范圍可以涵蓋0到2π,散射場的振幅則由天線的結(jié)構確定。

    圖1 (a)V型天線結(jié)構示意圖;(b)Y型天線結(jié)構示意圖[22] Fig.1 (a)Schematic of the V-shaped antenna. (b)Schematic of the Y-shaped antenna[22]

    V型天線的一種變型是Y型天線[25],如圖1(b)所示。V型天線臂長和夾角的改變會使得“對稱模式”和“反對稱模式”對應的諧振頻率產(chǎn)生漂移,為了改變這一特性,可以加入另一臂L作為一個新的調(diào)控自由度。當改變L的長度時,可以保證沿Y型天線對稱軸方向的偏振光諧振頻率不變,而垂直于對稱軸方向的諧振頻率可以獨立調(diào)節(jié)。當Y型天線具有C3對稱性時,即天線各臂長相等且各臂之間夾角為120°時,不存在偏振轉(zhuǎn)換特性。在線偏振光的激發(fā)下,V型天線與Y型天線的散射場一般來說是橢圓偏振光,當選擇合適的入射偏振態(tài)與出射偏振態(tài)的正交組合,可以實現(xiàn)出射光相位在0到2π范圍內(nèi)的精確控制。這類基于光學天線的相位調(diào)制機制依賴天線的諧振響應,因此本質(zhì)上具有窄帶的特性,并由于振幅與相位并不能獨立調(diào)制,需要仔細設計結(jié)構的參數(shù)。

    2.2惠更斯表面

    為了提高透射光的效率,Pfeiffer等人[26]提出了一種有別于天線諧振方法的表面結(jié)構。這種結(jié)構可以通過控制表面電、磁極化率(αe與αm)來達到完全消除反射的目的。表面等效電流與磁流正比于表面電極化率與磁極化率。通過調(diào)控表面極化率,并結(jié)合邊界條件,散射波前能夠獲得任意形式。當某一表面滿足

    式中,η0為自由空間阻抗,則其復振幅透射系數(shù)可以表示成[26]:

    圖2 (a)惠更斯表面樣品實物圖;(b)上半部分:實驗測得的輻射圖案;下半部分:測得的x-y平面磁場強度分布[26]。(c)光波段各向同性惠更斯超穎表面示意圖與單元結(jié)構示意圖;(d)在1.5 μm波長光照射下的異常折射光的仿真電場分布[27] Fig.2 (a)Photograph of the fabricated metasurfaces with anomalous refraction. (b)Upper panel: the measured radiative pattern; bottom panel: the distribution of measured magnetic field in x-y plane[26]. (c)Schematic of an optically thin isotropic Huygens′ metasurface and its unit cell. (d)Simulated electric field distribution of the anomalous beam at 1.5 μm[27]

    式中,ω為入射光頻率。當αe主要為實數(shù)時,則透射率接近1,并且透射相位可以涵蓋0到2π的范圍。Pfeiffer等人在微波領域進行了了原理驗證[26],如圖2(a)、2(b)所示。圖2(a)為樣品實物圖,圖中底部為樣品單元結(jié)構,基底頂側(cè)的微型銅線路可提供所需的電極化電流,底側(cè)的開口諧振環(huán)可以提供必要的磁極化電流。圖2(b)為所得實驗結(jié)果,在達到光線偏折效果的同時獲得了高達86%的透過率?;莞贡砻娴脑O計原理在其他頻段同樣適用,光學領域的表面結(jié)構如圖2(c)所示[27]。Pfeiffer等人設計了一種各向同性結(jié)構,在1.5 μm波長處實現(xiàn)了對入射光的調(diào)制。這種各向同性結(jié)構能夠有效降低樣品加工方面的困難。可獨立對振幅與相位進行調(diào)控的,利用全介質(zhì)結(jié)構以及應用于全息領域的惠更斯超穎表面也于近期被報道[28-31]。

    2.3基于貝里相位(Pancharatnam-Berry Phase)的超穎表面

    上述兩種機理本質(zhì)上都是色散的,其相位、振幅調(diào)制特性在入射光頻率發(fā)生變化時往往不再適用。而基于貝里相位的超穎表面具有寬頻帶的特性,通過相同單元結(jié)構的不同空間朝向即可實現(xiàn)連續(xù)的相位調(diào)制,本質(zhì)上是一種無色散的調(diào)制機理。

    圖3 貝里相位與偏振演化路徑。(a)初態(tài)沿不同路徑演化到同一終態(tài)所對應立體角;(b)當圓偏振入射,經(jīng)過納米棒天線陣列,并取相反圓偏振時,光的偏振態(tài)演化路徑。偏振態(tài)演化路徑分別經(jīng)歷σ→φ1→σ與σ→φ2→σ[22],其中φ1和φ2分別為兩不同天線的方位角 Fig.3 Pancharatnam-Berry phase and the evolution of polarization. (a) Different polarization evolution path in Poincaré sphere. (b)The polarization evolution path for metasurfaces composed of nanorods array when illuminated with circularly polarized light and selected orthogonal handedness circularly polarized light. The polarization state undergo σ→φ1→σ and σ→φ2→σ respectively, where φ1 and φ2 are the azimuthal angle for two different nanorods, respectively

    貝里相位也稱為幾何相位,1956年,Pancharatnam在研究中指出,當一束偏振光在其偏振態(tài)經(jīng)過多次改變后,其初末偏振態(tài)之間的相位改變與其偏振態(tài)改變過程在邦加球表面劃過的軌跡有關[32],1984年,Berry[33]重新研究了這一現(xiàn)象,并指出其具有幾何性質(zhì)的一面。如圖3(a)所示,考慮邦加球上某一初始偏振態(tài)A在沿兩個不同路徑變化到達同一最終態(tài)D和D′。D和D′雖然偏振態(tài)相同,但二者相差一個相位差,即貝里相位,可由下式表示:

    式中,ΩABCDA為邦加球面上兩條演化路徑所圍立體角。當某一手性為σ的圓偏振光入射,經(jīng)天線散射作用后終態(tài)變?yōu)橄喾词中缘膱A偏振光σ。若中途經(jīng)過兩個空間電極化方向分為φ1和φ2的納米天線,比如納米棒,其空間朝向就是其空間電極化方向,它們產(chǎn)生與本身電極化方向一致的線偏振態(tài),偏振方向分別為φ1和φ2,分別對應邦加球赤道上兩個點。經(jīng)過兩段不同的演化路徑之后,根據(jù)貝里相位原理,最終相同的終態(tài)之間存在 2(φ2-φ1)的相位差,如圖3(b)所示。偏振變化與相位差的關系亦可以由瓊斯矩陣清楚得到[34-35]。一般來說,一個各向異性散射體的瓊斯矩陣可以表示為:

    其中

    為旋轉(zhuǎn)角為θ的旋轉(zhuǎn)矩陣,t1和t2分別為入射光沿著該各向異性散射體主軸方向的前向散射系數(shù)。若給定右旋圓偏振光入射,通過這一各向異性散射體后電矢量可以表示為[36]

    等號右邊第一項表示與入射圓偏振光旋向相同的散射光,第二項表示擁有與入射旋向相反的散射光,但是附帶了2θ的貝里相位。Kang等人[36]利用這一原理,設計了一種U型小孔相位梯度超穎表面。通過仔細設計U型孔的結(jié)構以及旋轉(zhuǎn)U型孔的朝向,可以使得t1與t2振幅相同但是相位差正好為π,于是入射圓偏振光可完全轉(zhuǎn)化為反旋向的圓偏振光并攜帶所需相位,并實現(xiàn)光束偏折與匯聚等作用,如圖4中(a)、(b)所示。改變?nèi)肷涔獾男蚩梢愿淖児饩€偏折方向以及將光線由匯聚變?yōu)閿U散效果。Huang等人利用貝里相位原理,采用棒型納米天線實現(xiàn)廣義折反射現(xiàn)象[37]、寬帶渦旋光生成(圖4(c))、以及三維納米全息[21]等應用。利用該原理的彩色全息方法也得到了證實[38]。

    圖4 基于貝里相位的幾類超穎表面及其應用。(a)手性可調(diào)異常透射;(b)手性選擇性雙極性透鏡[36];(c)手性選擇性寬帶渦旋光束發(fā)生器[37] Fig.4 Several metasurfaces and their applications based on Berry phase. (a)Helicity-dependent tunable anomalous transmission. (b)Helicity-dependent dual-polarity metalens[36]. (c)Helicity-dependent broadband vortex beam generator[37]

    值得注意的是,貝里相位其本質(zhì)機理是無色散的,不會因為入射光頻率變化引入相位失真,但是對于特定結(jié)構總是存在一定的操作帶寬。一方面,入射波長須數(shù)倍于超穎表面結(jié)構尺寸以保證亞波長的要求;另一方面,入射波長不能太大以保證天線的散射效率。其結(jié)構的散射效率仍取決于局域表面等離激元諧振(Localized Surface Plasmon Resonance)。

    2.4其它設計方法

    其中,E和B,Ec和Bc分別為兩互補屏的后向散射場。靈活運用互補結(jié)構可以增大超穎表面的設計靈活性[40-41]。

    圖5 兩類結(jié)構互補的超穎表面:(a)顆粒性天線陣列;(b)小孔型天線陣列[22] Fig.5 Two types of complementary metasurfaces.(a)An particle antenna array. (b)A complementary aperture array

    超穎表面設計亦可以使用反射陣列,反射陣列由金屬天線,薄介質(zhì)層以及金屬基底構成。這類反射陣列受到早先有關微波與厘米波研究的啟發(fā)[42],該方法利用了天線與其鏡像天線的耦合,從而達到更廣的相位響應范圍。

    相當多的研究關注亞波長金屬結(jié)構,但是金屬在光頻段的歐姆損耗限制了超穎表面性能的提高。而低損耗,高折射率介質(zhì)構成的超穎表面受到了越來越多的關注。根據(jù)米氏散射理論,介質(zhì)可在光波段發(fā)生強烈的電磁散射,散射后的電磁場可以寫為電磁多級子散射場的級數(shù)和的形式[43]。與金屬結(jié)構主要由電子諧振影響結(jié)構性質(zhì)不同,在幾何結(jié)構規(guī)則的介質(zhì)諧振子中,最低階的諧振模式是磁偶極諧振,二階模式是電偶極諧振[44-46]。更高階的模式對總散射場貢獻很小,一般可以被忽略[47]。

    利用介質(zhì)作為超穎材料的結(jié)構單元也曾吸引了相當多的研究興趣,并在不同波段實現(xiàn)了負的電磁響應[48-50]。一般來說,電、磁諧振分別處于不同的頻率時,介質(zhì)諧振子只能提供最多 π的相位變化[2];為了使相位變化涵蓋2π的范圍,可以改變介質(zhì)單元的結(jié)構[51]或者是結(jié)構之間的間距[52],使得電、磁偶極子諧振重合,實現(xiàn)完全的相位調(diào)控。一系列結(jié)構諸如塊狀[53]、圓柱形[52,54]、橢圓柱形[55]等納米結(jié)構均可以實現(xiàn)對相位的完全調(diào)控甚至同時對偏振與相位完全調(diào)控[55]。而基于貝里相位的原理,也可利用介質(zhì)材料實現(xiàn)對光場相位的完全調(diào)控,例如采用具有連續(xù)方位角變化的介質(zhì)條結(jié)構[56]。

    3 超穎表面功能應用研究現(xiàn)狀

    3.1基于超穎表面的廣義折射和反射效應

    光在經(jīng)過兩種均勻材料截面時,根據(jù)斯涅耳定律以及菲涅爾公式可以確定光的反射角、透射角、反射系數(shù)以及透射系數(shù)。若光經(jīng)過界面時,在亞波長尺度內(nèi)產(chǎn)生了顯著的相位突變,則光線在經(jīng)過該界面時便不再遵守傳統(tǒng)的斯涅耳定律,而新的傳播規(guī)律仍然可以由費馬原理(Fermat′s Principle)推導而得[19]。

    式中,λ0是真空中的波長,dΦ/dx是平面上的相位梯度。如果dΦ/dx=0,則廣義折射與反射定律退化為傳統(tǒng)折射與反射定律,當dΦ/dx≠0時,偏離傳統(tǒng)定律的折射光與反射光稱為異常折射光與異常反射光。

    對于更一般的情況,即入射波矢與反射、折射波矢不處于同一平面時,費馬原理同樣適用,由此可以得到更一般的廣義折射與反射定律[57]:

    公式中的角度詳見圖6,式(8)與(9)分別代表平行或垂直于入射面的相位梯度。

    圖6 一般情況的廣義折射、反射示意圖,沿界面表面的相位梯度為dφ/dr[1,57] Fig.6 Schematic of generalized anomalous refraction and reflection; the interfacial phase gradient is dφ/dr[1,57]

    圖7 (a)V型天線陣列SEM樣品圖像;(b)實驗測量的異常折射數(shù)據(jù)。(c)納米棒陣列SEM樣品圖像;(d)入射手性σ=1時正常反射現(xiàn)象與異常反射現(xiàn)象 Fig.7 (a)SEM image of array of V-shaped antennas. (b)Measured data of anomalous refraction. (c)SEM image of array of nanorods. (d)Ordinary and anomalous reflection with incident helicity σ=1

    Yu等人[19]最早驗證了這種廣義折反射現(xiàn)象,如圖7(a)、7(b)所示,通過仔細選擇V型天線的結(jié)構,保證線偏振光垂直入射下各個V型天線沿其正交偏振方向的散射光振幅近似相等且相鄰單元散射光之間擁有恒定的相位差,從而保證恒定的相位梯度,則可以實現(xiàn)線偏振垂直入射下相應正交偏振折射光線任意偏折。Huang等人采用納米棒陣列驗證了依賴手性的廣義折射與反射定律,如圖7(c)、7(d)所示,即利用貝里相位原理控制相鄰單位之間的相位差,并提出了廣義光柵方程[37],即當相鄰納米棒的間隔與波長相近時,該周期結(jié)構會產(chǎn)生高衍射級次,各衍射級次傳播方向滿足:

    ntsin(θm)-nisin(θi)=

    式中,m代表衍射級次,σ代表入射圓偏振光的手性,可取±1。當m≠0時,存在異常衍射光。

    這種廣義反射與折射現(xiàn)象同樣可以應用于其它波段。Zhang等人[58]利用一種C型天線,在太赫茲波段實現(xiàn)了寬頻段的波束偏折。多種利用反射陣列結(jié)構實現(xiàn)廣義反射現(xiàn)象也得到證實[59-60]:Sun等人實現(xiàn)了H型反射天線[59]與塊狀反射天線[60]的廣義反射,理想情況下,所有入射能量都會轉(zhuǎn)化為異常反射能量,利用塊狀反射天線可測得反射陣列超穎表面效率達80%,遠大于效率僅為10%~20%的透射性V型天線。

    3.2基于超穎表面的任意相位梯度調(diào)制

    3.2.1 基于超穎表面的平板透鏡

    透鏡是光學系統(tǒng)中的重要元件,具有成像,照明,聚焦等功能[43]。傳統(tǒng)透鏡需要利用光波在透鏡材料中傳播的相位積累來達到控制波前的目的,因此需要設計透鏡的表面曲率以及采用合適的厚度,故傳統(tǒng)透鏡一般比較大。使用衍射光學元件[61]和梯度折射率結(jié)構[62]可以為透鏡設計帶來新的思路,可實現(xiàn)平板透鏡,但由于自然界可用折射率材料有限,以及面臨加工方面的困難,這種類型的平板透鏡仍受到集成化、超薄化的挑戰(zhàn)。

    為實現(xiàn)平板透鏡,超穎表面需要施加一個相位輪廓

    用于將入射平面波轉(zhuǎn)化為匯聚在焦距f處的球面波。顯然若超穎表面相位滿足上述分布,可以在實現(xiàn)高數(shù)值孔徑的同時不會帶來球面像差。Aieta等人[63]利用V型天線實現(xiàn)了1.55 μm波長處的平板透鏡,如圖8(a)所示。但這種平板透鏡的效率相當?shù)?,僅有1%左右的入射光被匯聚。這是因為實驗中僅采用了單層結(jié)構,且樣品填充率比較低,并只采集了正交偏振光。為提高操作波長帶寬以及轉(zhuǎn)換效率,可以考慮使用反射型超穎平板透鏡[64],如圖8(b)。同樣,出于抑制背景光的目的,可以通過巴比涅原理[65]設計V型小孔天線陣列匯聚入射光。Sergei等人[66]通過計算全息圖的原理,設計了一種反射式匯聚超穎表面,該超穎表面相位分布與計算全息圖所得相位一致,從而可以獲得單點以及多點聚焦的功能。

    圖8 (a)V型天線陣列構成的平板透鏡[63];(b)反射型超穎表面平板透鏡結(jié)構單元及其仿真所得焦平面處的電場強度分布[64] Fig.8 (a)Flat lens consisting of an array of V-shaped antennas[63]. (b) The unit cell of flat lens made of reflect-array and the calculated electric field distribution on the focus plane[64]

    基于貝里相位,可以利用納米棒[67]或者U型小孔[36]設計平板透鏡,這種平板透鏡一般可以根據(jù)入射光的手性變更匯聚或者分散的功能。同樣可以設計反射陣列[68]以達到增大轉(zhuǎn)換效率的目的。

    通過介質(zhì)超穎表面也能夠?qū)崿F(xiàn)高效率平板透鏡。通過采用高對比透射陣列(High-contrast Transmit Arrays,HCTAs),利用單晶硅圓柱的半徑變化實現(xiàn)對透射相位的調(diào)制,僅僅采用單層結(jié)構便實現(xiàn)了85%的透射效率以及72%的聚焦效率[54]。通過將圓柱換為橢圓結(jié)構可以獲得更多的調(diào)控自由度,可以實現(xiàn)x、y入射偏振光分別匯聚到不同側(cè)的效果[55]。最近,Mohammadreza等人[69]報道了一種工作在可見光波段,成像質(zhì)量可匹敵目前最先進商業(yè)物鏡的平板透鏡。該平板透鏡基本原理基于貝里相位,通過設計鰭狀TiO2柱的參數(shù)使其達到半波片的作用從而最大化轉(zhuǎn)換效率。精妙的加工工藝進一步完善了該平板透鏡的性能。該方法制得的透鏡最高獲得了86%的聚焦效率并可分辨亞波長間距的納米尺度特征,放大率可達170倍[69]。

    一般來說上述平板透鏡雖然沒有球面像差,但是單色像差(Monochromatic Aberration)總是存在。當入射光不再垂直于平板透鏡表面時,由于非垂直入射會引入一個附加的線性相位分布,透射光或者是反射光波前便不再是球面。此外平板透鏡通常也具有色差(Chromatic Aberration),通過設計耦合方型諧振子陣列(Rectangular Dielectric Resonators,RDRs),可以使超穎表面擁有一個依賴于波長的相位分布來補償這種色散累積的相位差,從而實現(xiàn)在1.3、1.5以及1.8 μm三個波長均能匯聚至同一焦點[70],但是在其他波長處仍然存在著色差。寬波帶消色差平板透鏡的實現(xiàn)仍然具有一定的挑戰(zhàn)性。

    3.2.2 基于超穎表面的計算全息

    全息圖記錄的不是一幅圖像,而是看似雜亂的一系列振幅與相位信息,利用參考光照射全息圖即可重構出虛擬的三維圖像[71]。計算全息[72-74](Computer Generated Hologram, CGH)技術可以方便的獲得重構目標物體所需的相位分布,超穎表面則可提供一系列對局部振幅、相位以及偏振的調(diào)控方法,因此是實現(xiàn)亞波長像素計算全息的一種有效工具。

    Walther等人[75]提出一種金屬孔徑陣列,通過調(diào)節(jié)孔徑的大小能夠調(diào)控局部單元的振幅,從而實現(xiàn)遠場全息圖像。利用孔徑的色散,同一樣品在905 nm以及1 385 nm處可以獲得兩個完全不同的全息圖,如圖9(a)所示。Ni等人[76]設計一系列V型小孔,如圖9(b),使其可獲得8級相位分布以及2級振幅分布。通過使用這小孔陣列在超穎表面近似地引入所需的振幅和相位,能夠在遠場獲得所需的圖像。

    基于貝里相位的全息也被廣泛研究。由棒型納米天線的相位調(diào)控特性可知[37],通過簡單旋轉(zhuǎn)納米天線的朝向,即可記錄三維物體計算全息圖所需的相位信息,其方位角編碼了連續(xù)可調(diào)的局域相位信息。對某一手性入射光,取相反手性的出射圓偏振光即可獲得所需光學再現(xiàn)像[21],如圖9(c)所示。實驗上可觀察到實心三維飛機模型和一系列空心螺旋點陣。這種方法可以有效增大三維全息的視場角,實現(xiàn)零級光束同軸再現(xiàn),并可以避免多級衍射像。單層透射型納米棒陣列的成像效率往往不夠高,適當優(yōu)化單元結(jié)構并且采用反射陣列,可以獲得更高的轉(zhuǎn)換效率。Zheng等人[77]設計了一種反射陣列(圖9(d)),使得沿納米棒長軸與沿納米棒短軸方向的線偏振光反射系數(shù)相位差在600~1 000 nm的范圍內(nèi)為π,從而該結(jié)構能起到半波片的作用,對于反射陣列可使得反射光完全轉(zhuǎn)化為同旋向圓偏振光(反射陣列同旋向圓偏振光攜帶貝里相位信息),極大地提高效率。同時Zheng等人采用了16級相位分布。實驗證實了這一結(jié)構在825 nm波長處效率可達85%,并且在630~1 050 nm波長范圍內(nèi)擁有相當高的成像質(zhì)量。

    圖9 基于超穎表面的全息。(a)金屬小孔陣列,以及入射波長分別為λ1=905 nm以及λ2=1 385 nm的遠場再現(xiàn)圖[75];(b)利用V型小孔實現(xiàn)的全息的樣品SEM圖[78];(c)基于貝里相位納米棒陣列三維全息示意圖[21];(d)基于貝里相位的高效率反射陣列全息單元結(jié)構及其實驗測得的遠場全息圖[77] Fig.9 Hologram based on metasurfaces. (a) an array of metal aperture and the far-field image of metasurface at λ1=905 nm and λ2=1 385 nm, respectively[75]. (b) SEM image of metasurface consisting of V-shaped apertures to realize holography[78]. (c)Schematic of metasurface hologram composed of nanorods array based on Pancharatnam-Berry phase principle[21]. (d)Unit cell of reflect-array metasurfacehologram with high efficiency base on Pancharatnam-Berry phase and the measured image in the far field[77]

    在此基礎上,為進一步提高全息圖信息容量,Huang等人提出了一種基于幾何相位的全息圖復用方法[79]。通過具有不同相移全息圖的線性疊加,可以得到總的合成全息圖復振幅為:

    Htot(u,v;λ,z,θ)=

    (12)

    圖10 全息圖的不同復用方式[79]:(a)可在不同聚焦位置以及旋轉(zhuǎn)角觀察全息圖的實驗裝置以及固定觀察距離分別在不同圓偏振光入射時所成再現(xiàn)像。(b)同一樣品不同z平面再現(xiàn)像;(c)同一樣品不同距離以及離軸角觀測的再現(xiàn)像 Fig.10 Different forms of holographic multiplexing[79]. (a)Experimental setup for observing the holographic images at separate focus positions and rotational angles and the reconstructions of image for different incident helicity. (b)The measured images for different incident helicity at different z planes. (c)Schematic of the observation for the four encoded objects with separate off-axis angles and reconstruction distances and the experimental reconstruction images

    式中,Htot是總合成全息圖的復振幅,Hl是單個全息圖復振幅或者其復共軛,p為全息圖總數(shù),u和v是空間頻率,xl和yl分別為每幅圖的空間位移。理論上p可以是任何整數(shù),但是實際上由于超穎表面有限的分辨率以及空間帶寬積,p不能夠過大。根據(jù)上述概念,Huang等人[79]實驗證明了該復用方法。例如:在同一位置,不同入射圓偏振態(tài)下可分別再現(xiàn)“BIT”以及?;請D像(圖10(a));右旋光入射,左旋光出射,在不同成像距離可得到不同再現(xiàn)像(圖10(b));以及在不同離軸角度與成像距離(圖10(c))均可以實現(xiàn)全息圖復用。這種方法能有效減少加工制造的壓力并且大幅提高單一樣品的信息容量。

    利用惠更斯超穎表面也可以實現(xiàn)納米全息應用。Zhao等人[31]利用惠更斯超穎表面的原理設計了一種生成全息圖的方式,通過在一定范圍掃描結(jié)構高度與半徑,找出介質(zhì)圓柱適宜的高度與半徑以滿足相應的表面極化率要求以引入再現(xiàn)全息圖所需的相位,并保持振幅一致。為了減少樣品制備難度,設計時采用固定介質(zhì)圓柱半徑差來對應相位突變的思路,而不是根據(jù)固定相位差選擇結(jié)構參數(shù)。實驗最終獲得了高達86%的透射效率以及23.6%的光學效率。

    近期有關彩色全息的研究逐漸引起了人們的興趣。Huang等人[80]利用鋁納米棒實現(xiàn)了多色反射式超穎全息。鋁作為有著更高等離子體諧振頻率的材料,可以將表面等離子諧振頻率從可見光區(qū)擴展到紫外區(qū)[81-82],這對于生成全彩全息圖有著重要意義。通過單一調(diào)節(jié)鋁納米棒的長度,可實現(xiàn)在不同波長處相位的二元調(diào)節(jié),即產(chǎn)生0、π的相位差[80]。這樣能夠顯著降低建造難度,但也造成了多級衍射的問題。通過在一個單元中劃分4個子單元,每個子單元對某一波長獨立響應,這樣在白光入射時便形成了多色全息圖。

    Wang等人[38]在可見光頻段分別實現(xiàn)了消色差以及高色散的介質(zhì)超穎表面全息應用,其相位操控基本原理是利用貝里相位,通過旋轉(zhuǎn)納米塊的朝向來獲得所需相位突變;但Wang等人通過結(jié)構設計使得各個納米塊起到半波片的作用以提高偏振轉(zhuǎn)化效率。該超穎表面的一個單元結(jié)構同樣包含有4個子單元,分別對紅、黃和藍光響應,其中有兩個子單元對應藍光響應以提高效率。在消色差應用中,各子單元擁有相同的相位編碼,不同入射光可生成相同的再現(xiàn)像,通過調(diào)節(jié)紅黃藍3種不同顏色入射光的能量即可獲得任意顏色的全息圖。而在高色散應用中,各子單元相位分布不一致,針對3種入射波長可分別生成不同圖案。Wang等人[38]利用該思路生成了一幅擁有多種色彩的花盆圖案。

    3.2.3 基于超穎表面的光束軌道角動量操控

    渦旋光束是指渦旋中心光強為零且具有螺旋等相位面分布的一種奇異光束[83]。在近軸傳播條件下,渦旋光束的波前可用exp(ilφ)描述[84],l被稱為拓撲電荷數(shù),可取分數(shù)或者整數(shù),其攜帶的軌道角動量L=l?,?為普朗克常數(shù)。以軌道角動量作為信息載體可顯著提高信息容量,在與自旋或者角動量相關的量子信息或粒子操縱等領域有廣泛的應用前景[85-86]。通常生成渦旋光束的方法包括空間光調(diào)制器、計算全息片、螺旋相位板,達曼光柵等[83]。這些方法所用器件通常較大或者僅適用于窄帶范圍。利用超穎表面生成渦旋光的方法可以實現(xiàn)超薄化,并能具有寬帶的特性。

    圖11 利用超穎表面生成渦旋光:(a)用V型天線超穎表面構成的相位輪廓以及生成的干涉圖案[19];(b)可見光頻率空間復用超穎表面示意圖[90];(c)開口環(huán)寬帶渦旋光生成超穎表面SEM圖像以及其左旋光與線偏振光入射下測得的透射圖案[87];(d)從左至右依次為分割式(Segmented)、插入式(Interleaved)與諧波響應(Harmonic Response)超穎表面及其形成的遠場強度分布示意圖[91] Fig.11 Vortex beam generation based on metasurfaces: (a)A phase profile created using V-shaped antenna metasurfaces and the interference patterns[19]; (b)Schematic of visible-frequency metasurfaces for spatially multiplexing optical vortices[90]; (c)SEM image of the split-ring metasurface which was designed for generating an optical vortex beam and the measured transmitted patterns under the incident beam with left-handed circular polarization and linear polarization, respectively[87]; (d)Schematic of far-field intensity distributions from segmented, interleaved, and harmonic response metasurfaces(from left to right)[91]

    Yu等人[19]最早關于V型天線的報道中即驗證了利用超穎表面實現(xiàn)渦旋光束的生成,如圖11(a)所示。通過將超穎表面四象限按照方位角平分為8份,按方位角的增加每部分增加π/4的相位差,可獲得拓撲電荷數(shù)l=1的渦旋光,由相位的奇異性可獲得中心暗區(qū)域[83]。通過與同方向傳播的高斯光束共軸傳播,可以獲得螺旋狀干涉條紋;若二者傾斜一定角度,則可獲得叉型干涉條紋?;谪惱锵辔辉淼募{米棒[37](圖4(c))或者開口環(huán)結(jié)構[87](圖11(c)),均可以實現(xiàn)寬光譜的渦旋光束生成,且在翻轉(zhuǎn)入射圓偏振光時,所得的渦旋光拓撲電荷數(shù)正負異號。Yang等人[53]采用介質(zhì)結(jié)構,獲得了更高的渦旋光轉(zhuǎn)換效率。對渦旋光的并行處理,包括不同渦旋光本征態(tài)的空間復用與解復用[88, 89]等應用,受到了一系列的關注。通過適當結(jié)構設計,單個樣品可實現(xiàn)同時生成與調(diào)控多種渦旋光,實現(xiàn)渦旋光的多路復用。Mehmood等人[90]利用同心圓環(huán)分布的納米棒狀孔實現(xiàn)了可見光波段渦旋光復用,該結(jié)構可同時生成3種拓撲電荷數(shù)的渦旋光并分別在距樣品表面不同距離聚焦,如圖11(b)。Maguid討論了超穎表面信息容量的問題[91]。實現(xiàn)復用時,單元結(jié)構的表面分布方式,如圖11(d),會影響復用時超穎表面的信息容納能力。

    近年來利用達曼光柵的原理制作相位板對渦旋光進行調(diào)控受到了人們的關注。在傳統(tǒng)的二元相位光柵中,達曼相位編碼方法可生成能量均勻的各衍射級次[92]。達曼渦旋光柵(Dammann Optical Vortex Grating)可同時實現(xiàn)不同拓撲電荷數(shù)渦旋光的復用與解復用[93]。該結(jié)構的設計思路是首先得到各級次不同拓撲電荷數(shù)渦旋光表面純相位疊加,再利用達曼光柵編碼方法進行二元相位分布設計。進一步利用達曼空間波帶片[94](Spiral Dammann Zone Plate)則可以實現(xiàn)三維渦旋陣列的生成與調(diào)控[95-96]。上述方法中,若采用超穎表面,則有望擺脫二元相位帶來的限制,降低設計難度。

    3.3偏振轉(zhuǎn)換

    對于許多現(xiàn)代光電應用,控制電磁波的偏振轉(zhuǎn)換有著重要的意義。由于超穎表面超薄、緊湊以及在設計靈活性上的特點,利用超穎表面實現(xiàn)偏振轉(zhuǎn)換引起了越來越多研究者的興趣。一般來說,高度的對稱性有利于超穎結(jié)構單元維持偏振態(tài),而對稱性的破缺則可以實現(xiàn)對偏振態(tài)的調(diào)控。

    圖12 (a)由V型天線構成的寬波段四分之一波片結(jié)構示意圖;(b)所得圓偏振光的偏振分析[20];(c)THz波段實現(xiàn)線偏振到圓偏振轉(zhuǎn)換樣品圖;(d)實驗所得(c)中樣品txx和tyy振幅,相位差以及橢偏度[98];(e)用以實現(xiàn)圓-圓偏振轉(zhuǎn)換的樣品示意圖;(f)仿真與實驗所得透射率[102] Fig.12 (a)Schematic of a broadband quarter-wave plate consisting of V-shaped antennas. (b)State-of-polarization analysis for the obtained circular polarization light[20]. (c)An image of the fabricated sample which produces the linear to circular polarization conversion at THz band. (d) Experimentally measured transmission amplitude, phase retardation, and ellipticity for the sample exhibited in (c)[98]. (e)Schematic of sample to realize asymmetric transmission and polarization conversion. (f)Simulated and measured transmission[102]

    Wang等人[97]利用單層單元結(jié)構為十字形孔的超穎表面制得了厚度僅為0.29λ的線偏振到圓偏振的轉(zhuǎn)換器件,但這種結(jié)構利用了結(jié)構的窄帶諧振特性,工作帶寬很小。利用V型天線,Yu等人[20]設計了一種寬波段四分之一波片,并能產(chǎn)生高質(zhì)量的圓偏振光,其工作波段可從5 μm覆蓋至12 μm,如圖12(a)所示。其單元結(jié)構包含兩個子單元,各個子單元當線偏振光入射時可以產(chǎn)生兩束同方向傳播的,偏振方向相互正交的線偏振光,并通過子單元之間的錯位提供π/2的相位差。除了寬波帶的特性,這種結(jié)構還利用廣義折射現(xiàn)象,使得所獲得圓偏振光與原垂直入射光之間存在一定偏折,便于他們之間的分離。其生成圓偏振光的偏振分析如圖12(b)所示,可通過旋轉(zhuǎn)置于探測器前的偏振片并測量透射能量獲得。

    采用雙層[98-99]或者多層[100]結(jié)構,利用近場或者類法布里-珀羅(Fabry-Pérot)耦合,能有效提高偏振轉(zhuǎn)換效率。Cong等人[98]提出了一種THz波段雙層堆疊金屬線網(wǎng)格構成的偏振轉(zhuǎn)換器件,該器件可實現(xiàn)線偏振到圓偏振的轉(zhuǎn)換。其第一層的金屬線網(wǎng)格與x軸成45°角以使得其垂直方向的透射系數(shù)|txx|和|tyy|振幅一致,但是其相位差依賴頻率變化;第二層的網(wǎng)格則能補償這種相位差,使得x、y分量近似保持90°相位差,如圖12(c)、12(d)所示。Pfeiffer等人[101-102]設計另一種用三層結(jié)構超穎表面實現(xiàn)非對稱圓偏振的轉(zhuǎn)化方法,通過構建各層各向異性結(jié)構參數(shù),可以實現(xiàn)對整個透射瓊斯矩陣的設計,其結(jié)構與實驗如圖12(e)、12(f)所示。這種結(jié)構可以在寬波長范圍內(nèi)對于右旋到左旋的透射率可達到50%,而其他透射分量不到2.5%。

    圖13 基于超穎表面旋光應用。(a)實現(xiàn)線偏振旋轉(zhuǎn)的反射型超穎表面示意圖;(b)實現(xiàn)旋光的透射型超穎表面示意圖;(c)實驗測得反射型結(jié)構同偏振與正交偏振反射率;(d)實驗測得透射型旋光結(jié)構正交偏振透射率和同偏振反射率[103] Fig.13 Linear polarization rotation based on metasurfaces. (a)Schematic of the metasurface to rotate the incident linearly polarized light in reflection; (b)Schematic of the metasurface to rotate the incident linearly polarized light in transmission. (c)Experimentally measured co- and cross-polarized reflection for the reflection type. (d)Experimentally measured cross-polarized transmission and co-polarized reflection for the transmission type[103]

    除了線偏振到圓偏振的轉(zhuǎn)換,線偏振態(tài)的旋轉(zhuǎn),即旋光性,也受到一系列關注。Grady等人[103]提出了一種高效的反射與透射陣列如圖13所示。反射型線偏振旋轉(zhuǎn)結(jié)構由一層金屬短線(Cut-Wire)陣列、一層介質(zhì)中間層以及金屬基底層構成,這3層結(jié)構形成一種類法布里-珀羅腔的結(jié)構,通過設計各結(jié)構參數(shù)使得同偏振方向達到干涉相消條件,可以實現(xiàn)正交偏振方向達到80%的反射率而同偏振方向反射率極低,如圖13(a)、13(c)所示。在金屬短線陣列兩側(cè)放置方向垂直的光柵結(jié)構則可以實現(xiàn)透射型的90°偏振旋轉(zhuǎn),如圖13(b)、13(d)。值得注意的是,在透射型結(jié)構中將金屬短線陣列換成相位操控結(jié)構,則可以在偏振轉(zhuǎn)換的同時發(fā)生廣義折射等現(xiàn)象,極大提高設計靈活性。利用該思路,超寬頻帶[104],工作在可見光波段[105]以及利用介質(zhì)超穎表面[53]的線偏振旋轉(zhuǎn)得到驗證。特別是,若采用介質(zhì)結(jié)構,即一層長方體硅結(jié)構與金屬基平面中間間隔一層PMMA,由于介質(zhì)結(jié)構的低損耗,可以在近紅外波長處實現(xiàn)高達98%的偏振轉(zhuǎn)換效率以及200 nm的帶寬[53]。

    4 結(jié)束語

    超穎表面自提出之日起就備受關注。超穎表面將體超穎材料的維度從三維降到二維,有效降低了加工難度和歐姆損耗,不再依賴光在介質(zhì)中傳播帶來的相位積累,而是通過平面亞波長天線進行相位調(diào)制,使其可稱為超薄器件。通過合適的單元結(jié)構,超穎表面可以對光場振幅、相位、偏振進行調(diào)控,展示了對光場強大的操縱能力,其在光儲存、集成光學、量子信息處理等方面具有廣闊的應用潛力。本文簡要回顧了最近幾年發(fā)展起來的超穎表面不同調(diào)制機理,以及在不同波段有關任意波前調(diào)控、偏振轉(zhuǎn)換等方面功能應用。有關超穎表面的研究領域非常廣泛,除了本文所述的內(nèi)容,還包含非線性[106]、電磁吸波材料[107]、表面波全息[108]以及主動可調(diào)超穎表面等研究領域。例如,通過物理性質(zhì)可調(diào)材料的引入,研究人員得以實現(xiàn)對表面結(jié)構諧振性質(zhì)的調(diào)控,進而實現(xiàn)可調(diào)超穎表面。報道較多的途徑有利用半導體[109-110]、石墨烯[111-113]或者相變材料[114-115]以實現(xiàn)諧振的調(diào)控。強烈的光與物質(zhì)相互作用可以使得超穎表面僅需要少量可調(diào)材料,這有助于實現(xiàn)超穎表面器件的快速調(diào)制。綜上,超穎表面與生俱來的超薄、緊湊的性質(zhì),以及空間逐點光場調(diào)控能力,使得其在未來集成化、微型化空間光場調(diào)制方面具有重要應用前景。

    [1] YU N,CAPASSO F. Flat optics with designer metasurfaces[J].NatureMaterials,2014,13(2):139-150.

    [2] CHEN H,TAYLOR A J,YU N. A review of metasurfaces:physics and applications[J].ArXiv:1605.07672,2016.

    [3] SMITH D R,PENDRY J B,WILTSHIRE M. Metamaterials and negative refractive index[J].Science,2004,305(5685):788-792.

    [4] PENDRY J B,HOLDEN A J,STEWART W J,etal.. Extremely low frequency plasmons in metallic mesostructures[J].PhysicalReviewLetters,1996,76(25):4773-4776.

    [5] PENDRY J B,HOLDEN A J,ROBBINS D J,etal.. Magnetism from conductors and enhanced nonlinear phenomena[J].IEEETransactionsonMicrowaveTheoryandTechniques,1999,47(11):2075-2084.

    [6] SMITH D R,PADILLA W J,VIER D C,etal.. Composite medium with simultaneously negative permeability and permittivity[J].PhysicalReviewLetters,2000,84(18):4184-4187.

    [7] SHELBY R A,SMITH D R,SCHULTZ S. Experimental verification of a negative index of refraction[J].Science,2001,292(5514):77-79.

    [8] ZHANG S,FAN W J,PANOIU N C,etal.. Experimental demonstration of near-infrared negative-index metamaterials[J].PhysicalReviewLetters,2005,95(13):137404.

    [9] VALENTINE J,ZHANG S,ZENTGRAF T,etal.. Three-dimensional optical metamaterial with a negative refractive index[J].Nature,2008,455(7211):376-379.

    [10] XIAO S,CHETTIAR U K,KILDISHEV A V,etal.. Yellow-light negative-index metamaterials[J].OpticsLetters,2009,34(22):3478-3480.

    [11] SMITH D R,SCHULTZ S,MARKOS P,etal.. Determination of effective permittivity and permeability of metamaterials from reflection and transmission coefficients[J].PhysicalReviewB,2002,65(19):195104.

    [12] SMITH D R,VIER D C,KOSCHNY T,etal.. Electromagnetic parameter retrieval from inhomogeneous metamaterials[J].PhysicalReviewEStatisticalNonlinear&SoftMatterPhysics,2005,71(3):142-154.

    [13] PENDRY J B. Negative refraction makes a perfect lens[J].PhysicalReviewLetters,2000,85(18):3966-3969.

    [14] FANG N,LEE H,SUN C,etal.. Sub-diffraction-limited optical imaging with a silver superlens[J].Science,2005,308(5721):534-537.

    [15] PENDRY J B,SCHURIG D,SMITH D R. Controlling electromagnetic fields[J].Science,2006,312(5781):1780-1782.

    [16] LEONHARDT U. Optical conformal mapping[J].Science,2006,312(5781):1777-1780.

    [17] HOLLOWAY C L,KUESTER E F,GORDON J A,etal.. An overview of the theory and applications of metasurfaces:the two-dimensional equivalents of metamaterials[J].IEEEAntennasandPropagationMagazine,2012,54(2):10-35.

    [18] KILDISHEV A V,BOLTASSEVA A,SHALAEV V M. Planar photonics with metasurfaces[J].Science,2013,339(6125):1232009.

    [19] YU N,GENEVET P,KATS M A,etal.. Light propagation with phase discontinuities:generalized laws of reflection and refraction[J].Science,2011,334(6054):333-337.

    [20] YU N,AIETA F,GENEVET P,etal.. A broadband, background-free quarter-wave plate based on plasmonic metasurfaces[J].NanoLetters,2012,12(12):6328-6333.

    [21] HUANG L,CHEN X,MüHLENBERND H,etal.. Three-dimensional optical holography using a plasmonic metasurface[J].NatureCommunications,2013,4:2808.

    [22] 黃玲玲.基于手性光場作用的超穎表面的相位調(diào)控特性及其應用[D].北京:清華大學,2014.

    HUANG L L. The phase modulation property of metasurfaces based on chiral field interaction and its applications[D]. Beijing:Tsinghua University,2014.(inChinese)

    [23] BHARADWAJ P,DEUTSCH B,NOVOTNY L. Optical antennas[J].AdvancesinOpticsandPhotonics,2009,1(3):438-483.

    [24] NOVOTNY L,VAN HULST N. Antennas for light[J].NaturePhotonics,2011,5(2):83-90.

    [25] KATS M A,GENEVET P,AOUST G,etal.. Giant birefringence in optical antenna arrays with widely tailorable optical anisotropy[J].ProceedingsoftheNationalAcademyofSciences,2012,109(31):12364-12368.

    [26] PFEIFFER C,GRBIC A. Metamaterial Huygens' surfaces:tailoring wave fronts with reflectionless sheets[J].PhysicalReviewLetters,2013,110(19):197401.

    [27] PFEIFFER C,EMANI N K,SHALTOUT A M,etal.. Efficient light bending with isotropic metamaterial Huygens' surfaces[J].NanoLetters,2014,14(5):2491-2497.

    [28] KIM M,WONG A M H,ELEFTHERIADES G V. Optical Huygens metasurfaces with independent control of the magnitude and phase of the local reflection coefficients[J].PhysicalReviewX,2014,4(4):041042.

    [29] WANG Z,SHI J,CHEN J. High-efficiency electromagnetic wave controlling with all-dielectric Huygens'metasurfaces[J].InternationalJ.AntennasandPropagation,2015:1-7.

    [30] DECKER M,STAUDE I,FALKNER M,etal.. High-efficiency dielectric Huygens' surfaces[J].AdvancedOpticalMaterials,2015,3(6):813-820.

    [31] ZHAO W,JIANG H,LIU B,etal.. Dielectric Huygens'metasurface for high-efficiency hologram operating in transmission mode[J].ScientificReports,2016,6:30613.

    [32] Generalized Theory of Interference and its Applications.Part 2:Partially Coherent Pencils[J].ProceedingsofIndianAcademyofSciences,1956,section A,4(6):398-417.

    [33] BERRY M V. Quantal phase-factors accompanying adiabatic changes[J].ProceedingsoftheRoyalSocietyofLondonSeriesA-MathematicalandPhysicalSciences,1984,392(1802):45-57.

    [34] MENZEL C,ROCKSTUHL C,LEDERER F. An advanced jones calculus for the classification of periodic metamaterials[J].Phys.Rev.A,2010,82(5):53811.

    [35] ARMITAGE N P. Constraints on jones transmission matrices from time-reversal invariance and discrete spatial symmetries[J].PhysicalReviewB,2014,90(3):35135.

    [36] KANG M,FENG T,WANG H T,etal.. Wave front engineering from an array of thin aperture antennas[J].OpticsExpress,2012,20(14):15882-15890.

    [37] HUANG L,CHEN X,M HLENBERND H,etal.. Dispersionless phase discontinuities for controlling light propagation[J].NanoLetters,2012,12(11):5750-5755.

    [38] WANG B,DONG F,LI Q,etal.. Visible-frequency dielectric metasurfaces for multiwavelength achromatic and highly dispersive holograms[J].NanoLetters,2016,16(8):5235-5240.

    [39] FALCONE F,LOPETEGI T,LASO M,etal.. Babinet principle applied to the design of metasurfaces and metamaterials[J].PhysicalReviewLetters,2004,93(19):197401.

    [40] ZENTGRAF T,MEYRATH T P,SEIDEL A,etal.. Babinet′S principle for optical frequency metamaterials and nanoantennas[J].PhysicalReviewB,2007,76(3):033407.

    [41] CHEN H T,O'HARA J F,TAYLOR A J,etal.. Complementary planar terahertz metamaterials[J].OptExpress,2007,15(3):1084-1095.

    [42] POZAR D M,TARGONSKI S D,SYRIGOS H D. Design of millimeter wave microstrip reflectarrays[J].IEEETransactionsonAntennasandPropagation,1997,45(2):287-296.

    [43] 馬科斯 玻恩,埃米爾沃爾夫.光學原理——光的傳播、干涉和衍射的電磁理論(第七版)[M].北京:電子工業(yè)出版社,2009.

    BORN M,WOLF E.PrinciplesofOptics:ElectromagneticTheoryofPropagation,InterferenceandDiffractionofLight[M]. Beijing:Electronic Industry Press,2009.(in Chinese)

    [44] EVLYUKHIN A B,REINHARDT C,CHICHKOV B N. Multipole light scattering by nonspherical nanoparticles in the discrete dipole approximation[J].PhysicalReviewB,2011,84(23):235429.

    [45] EVLYUKHIN A B,REINHARDT C,SEIDEL A,etal.. Optical response features of Si-nanoparticle arrays[J].PhysicalReviewB,2010,82(4):045404.

    [46] ZHAO Q,KANG L,DU B,etal.. Experimental demonstration of isotropic negative permeability in a three-dimensional dielectric composite[J].PhysicalReviewLetters,2008,101(2):027402.

    [47] VYNCK K,FELBACQ D,CENTENO E,etal.. All-dielectric rod-type metamaterials at optical frequencies[J].PhysicalReviewLetters,2009,102(13):133901.

    [48] ZHAO Q,ZHOU J,ZHANG F,etal.. Mie resonance-based dielectric metamaterials[J].MaterialsToday,2009,12(12):60-69.

    [49] PENG L,RAN L,CHEN H,etal.. Experimental observation of left-handed behavior in an array of standard dielectric resonators[J].PhysicalReviewLetters,2007,98(15):157403.

    [50] GINN J C,BRENER I,PETERS D W,etal.. Realizing optical magnetism from dielectric metamaterials[J].PhysicalReviewLetters,2012,108(9):097402.

    [51] CHENG J,ANSARI-OGHOL-BEIG D,MOSALLAEI H. Wave manipulation with designer dielectric metasurfaces[J].OpticsLetters,2014,39(21):6285-6288.

    [52] CHONG K E,STAUDE I,JAMES A,etal.. Polarization-independent silicon metadevices for efficient optical wavefront control[J].NanoLetters,2015,15(8):5369-5374.

    [53] YANG Y,WANG W,MOITRA P,etal.. Dielectric meta-reflectarray for broadband linear polarization conversion and optical vortex generation[J].NanoLetters,2014,14(3):1394-1399.

    [54] ARBABI A,HORIE Y,BALL A J,etal.. Subwavelength-thick lenses with high numerical apertures and large efficiency based on high-contrast transmitarrays[J].NatureCommunications,2015,6:7069.

    [55] ARBABI A,HORIE Y,BAGHERI M,etal.. Dielectric metasurfaces for complete control of phase and polarization with subwavelength spatial resolution and high transmission[J].NatureNanotechnology,2015,10(11):937-943.

    [56] LIN D,FAN P,HASMAN E,etal.. Dielectric gradient metasurface optical elements[J].Science,2014,345(6194):298-302.

    [57] AIETA F,GENEVET P,YU N,etal.. Out-of-plane reflection and refraction of light by anisotropic optical antenna metasurfaces with phase discontinuities[J].NanoLetters,2012,12(3):1702-1706.

    [58] ZHANG X,TIAN Z,YUE W,etal.. Broadband terahertz wave deflection based on c-shape complex metamaterials with phase discontinuities[J].AdvancedMaterials,2013,25(33):4567-4572.

    [59] SUN S,HE Q,XIAO S,etal.. Gradient-index meta-surfaces as a bridge linking propagating waves and surface waves[J].NatureMaterials,2012,11(5):426-431.

    [60] SUN S,YANG K,WANG C,etal.. High-efficiency broadband anomalous reflection by gradient meta-surfaces[J].NanoLetters,2012,12(12):6223-6229.

    [61] 金國藩,嚴瑛白,鄔敏賢.二元光學[M].北京:國防工業(yè)出版社,1998.

    JIN G F,YAN Y B,WU M X.BinaryOptics[M]. Beijing:National Defense Industry Press,1998.(inChinese)

    [62] MARCHAND E W.GradientIndexOptics[M]. New York:New York Academic Press,1978.

    [63] AIETA F,GENEVET P,KATS M A,etal.. Aberration-free ultrathin flat lenses and axicons at telecom wavelengths based on plasmonic metasurfaces[J].NanoLetters,2012,12(9):4932-4936.

    [64] LI X,XIAO S,CAI B,etal.. Flat metasurfaces to focus electromagnetic waves in reflection geometry[J].OpticsLetters,2012,37(23):4940-4942.

    [65] NI X,ISHII S,KILDISHEV A V,etal.. Ultra-thin, planar, babinet-inverted plasmonic metalenses[J].Light:Science&Applications,2013,2(4):e72.

    [66] KUZNETSOV S A,ASTAFEV M A,BERUETE M,etal.. Planar holographic metasurfaces for terahertz focusing[J].ScientificReports,2015,5:7738.

    [67] CHEN X,HUANG L,M HLENBERND H,etal.. Dual-polarity plasmonicmetalens for visible light[J].NatureCommunications,2012,3:1198.

    [68] PORS A,NIELSEN M G,ERIKSEN R L,etal.. Broadband focusing flat mirrors based on plasmonic gradient metasurfaces[J].NanoLetters,2013,13(2):829-834.

    [69] KHORASANINEJAD M,CHEN W T,DEVLIN R C,etal.. Metalenses at visible wavelengths: diffraction-limited focusing and subwavelength resolution imaging[J].Science,2016,352(6290):1190-1194.

    [70] KHORASANINEJAD M,AIETA F,KANHAIYA P,etal.. Achromatic metasurface lens at telecommunication wavelengths[J].NanoLetters,2015,15(8):5358-5362.

    [71] 謝敬輝,廖寧放,曹良才.傅里葉光學與現(xiàn)代光學基礎[M].北京:北京理工大學出版社,2007.

    XIE J H,LIAO N F,CAO L C.FundamentalsofFourierOpticsandContemporaryOptics[M]. Beijing:Beijing Institute of Technology Press,2007.(inChinese)

    [72] SLINGER C,CAMERON C,STANLEY M. Computer-generated holography as a generic display technology[J].Computer,2005,38(8):46-53.

    [73] KELLY D P,MONAGHAN D S,PANDEY N,etal.. Digital holographic capture and optoelectronic reconstruction for 3D displays[J].InternationalJ.DigitalMultimediaBroadcasting,2010,2010:1-14.

    [74] GENG J. Three-dimensional display technologies[J].AdvancesinOpticsandPhotonics,2013,5(4):456-535.

    [75] WALTHER B,HELGERT C,ROCKSTUHL C,etal.. Spatial and spectral light shaping with metamaterials[J].AdvancedMaterials,2012,24(47):6300-6304.

    [76] NI X,KILDISHEV A V,SHALAEV V M. Metasurface holograms for visible light[J].Nat.Comm ,2013,4:2807.

    [77] ZHENG G,M HLENBERND H,KENNEY M,etal.. Metasurface holograms reaching 80% efficiency[J].NatureNanotechnology,2015,10(4):308-312.

    [78] NI X,KILDISHEV A V,SHALAEV V M. Metasurface holograms for visible light[J].NatureCommunications,2013,4:2807.

    [79] HUANG L,M HLENBERND H,LI X,etal.. Broadband hybrid holographic multiplexing with geometric metasurfaces[J].AdvancedMaterials,2015,27(41):6444-6449.

    [80] HUANG Y,CHEN W T,TSAI W,etal.. Aluminum plasmonic multicolor meta-hologram[J].NanoLetters,2015,15(5):3122-3127.

    [81] KNIGHT M W,KING N S,LIU L,etal.. Aluminum for plasmonics[J].ACSNano,2014,8(1):834-840.

    [82] KNIGHT M W,LIU L,WANG Y,etal.. Aluminum plasmonic nanoantennas[J].NanoLetters,2012,12(11):6000-6004.

    [83] PADGETT M,COURTIAL J,ALLEN L. Light's orbital angular momentum[J].PhysicsToday,2004,57(5):35-40.

    [84] ALLEN L,BEIJERSBERGEN M W,SPREEUWR,etal.. Orbital angular-momentum of light and the transformation of Laguerre-Gaussian laser modes[J].PhysicalReviewA,1992,45(11):8185-8189.

    [85] MOLINA-TERRIZA G,TORRES J P,TORNER L. Twisted photons[J].NaturePhysics,2007,3(5):305-310.

    [86] VERBEECK J,TIAN H,SCHATTSCHNEIDER P. Production and application of electron vortex beams[J].Nature,2010,467(7313):301-304.

    [87] ZENG J,LI L,YANG X,etal.. Generating and separating twisted light by gradient rotation split-ring antenna metasurfaces[J].NanoLetters,2016,16(5):3101-3108.

    [88] LI S,WANG J. Simultaneous demultiplexing and steering of multiple orbital angular momentum modes[J].ScientificReports,2015,5:15406.

    [89] REN H,LI X,ZHANG Q,etal.. On-chip noninterference angular momentum multiplexing of broadband light[J].Science,2016,352(6287):805-809.

    [90] MEHMOOD M Q,MEI S,HUSSAIN S,etal.. Visible-frequency metasurface for structuring and spatially multiplexing optical vortices[J].AdvancedMaterials,2016,28(13):2533-2539.

    [91] MAGUID E,YULEVICH I,VEKSLER D,etal.. Photonic spin-controlled multifunctional shared-aperture antenna array[J].Science,2016,352(6290):1202-1206.

    [92] ZHANG N,YUAN X C,BURGE R E. Extending the detection range of optical vortices by dammann vortex gratings[J].Opt.Lett.,2010,35(20):3495-3497.

    [93] LEI T,ZHANG M,LI Y,etal.. Massive individual orbital angular momentum channels for multiplexing enabled by dammann gratings[J].Light:Science&Applications,2015,4(3):e257.

    [94] YU J,ZHOU C,JIA W,etal.. Generation of dipole vortex array using spiral dammann zone plates[J].Appl.Opt.,2012,51(28):6799-6804.

    [95] YU J,ZHOU C,JIA W,etal.. Three-dimensional dammann array[J].Appl.Opt.,2012,51(10):1619-1630.

    [96] YU J,ZHOU C,JIA W,etal.. Three-dimensional dammann vortex array with tunable topological charge[J].Appl.Opt.,2012,51(13):2485-2490.

    [97] WANG D,GU Y,GONG Y,etal.. An ultrathin terahertz quarter-wave plate using planar babinet-inverted metasurface[J].OpticsExpress,2015,23(9):11114-11122.

    [98] CONG L,XU N,GU J,etal.. Highly flexible broadband terahertz metamaterial quarter-wave plate[J].Laser&PhotonicsReviews,2014,8(4):626-632.

    [99] LI Y,ZHANG J,QU S,etal.. Achieving wide-band linear-to-circular polarization conversion using ultra-thin bi-layered metasurfaces[J].J.AppliedPhysic,2015,117(4):44501.

    [100] SUN W,HE Q,HAO J,etal.. A transparent metamaterial to manipulate electromagnetic wave polarizations[J].OpticsLetters,2011,36(6):927-929.

    [101] PFEIFFER C,GRBIC A. Bianisotropic metasurfaces for optimal polarization control:analysis and synthesis[J].PhysicalReviewApplied,2014,2(4):044011.

    [102] PFEIFFER C,ZHANG C,RAY V,etal.. High performance bianisotropicmetasurfaces:asymmetric transmission of light[J].PhysicalReviewLetters,2014,113(2):023902.

    [103] GRADY N K,HEYES J E,CHOWDHURY D R,etal.. Terahertz metamaterials for linear polarization conversion and anomalous refraction[J].Science,2013,340(6138):1304-1307.

    [104] CHEN H,WANG J,MA H,etal.. Ultra-wideband polarization conversion metasurfaces based on multiple plasmon resonances[J].J.AppliedPhysics,2014,115(15):154504.

    [105] DAI Y,REN W,CAI H,etal.. Realizing full visible spectrum metamaterial half-wave plates with patterned metal nanoarray/insulator/metal film structure[J].OpticsExpress,2014,22(7):7465-7472.

    [106] MINOVICH A E,MIROSHNICHENKO A E,BYKOV A Y,etal.. Functional and nonlinear optical metasurfaces[J].Laser&PhotonicsReviews,2015,9(2):195-213.

    [107] CUI Y,HE Y,JIN Y,etal.. Plasmonic and metamaterial structures as electromagnetic absorbers[J].Laser&PhotonicsReviews,2014,8(4):495-520.

    [108] GENEVET P,CAPASSO F. Holographic optical metasurfaces:a review of current progress[J].ReportsonProgressinPhysics,2015,78(2):24401.

    [109] GU J,SINGH R,LIU X,etal.. Active control of electromagnetically induced transparency analogue in terahertz metamaterials[J].NatureCommunications,2012,3:1151.

    [110] WATTS C M,SHREKENHAMER D,MONTOYA J,etal.. Terahertz compressive imaging with metamaterial spatial light modulators[J].NaturePhotonics,2014,8(8):605-609.

    [111] HUIDOBRO P A,KRAFT M,MAIER S A,etal.. Graphene as a tunable anisotropic or isotropic plasmonic metasurface[J].ACSNano,2016,10(5):5499-5506.

    [112] DABIDIAN N,DUTTA-GUPTA S,KHOLMANOV I,etal.. Experimental demonstration of phase modulation and motion sensing using graphene-integrated metasurfaces[J].NanoLetters,2016,16(6):3607-3615.

    [113] LI Z,YU N. Modulation of mid-infrared light using graphene-metal plasmonic antennas[J].AppliedPhysicsLetters,2013,102(13):131108.

    [114] KATS M A,BLANCHARD R,GENEVET P,etal.. Thermal tuning of mid-infrared plasmonic antenna arrays using a phase change material[J].OpticsLetters,2013,38(3):368-370.

    [115] YIN X,SCH FERLING M,MICHEL A U,etal.. Active chiral plasmonics[J].NanoLetters,2015,15(7):4255-4260.

    Theprincipleandresearchprogressofmetasurfaces

    LI Tian-you, HUANG Ling-ling*, WANG Yong-tian

    (SchoolofOptoelectronics,BeijingInstituteofTechnology,Beijing100081,China)

    Metasurfaces, a specific type of subwavelength-scale two-dimensional(2D) array, have shown great flexibilities to control the light field. Metasurfaces possess the unique advantages of ultrathin, low-loss and easy fabrication, therefore attract immense attentions and show broad prospects for applications. This paper reviews recent progress in the physics of metasurfaces and the applications for generalized refraction and reflection, polarization conversion, controlling the orbit angular moment of light and computer generated holography(CGH). Due to the extra design freedom and powerful functionality, metasurfaces allow the possibility to substitute the conventional optical components in a variety of applications. In the future, novel types of metasurfaces with broadband, high efficiency, multi-task and active tunability are still highly desired and of great value.

    metasurface;polarization conversion;orbit angular momentum;computer generated holography(CGH)

    TP394.1; TH691.9

    A

    10.3788/CO.20171005.0523

    李天佑(1994—),男,湖南岳陽人,碩士研究生,2015年于北京理工大學獲得學士學位,主要從事超穎表面功能器件方面的研究。E-mail:19940416@bit.edu.cn

    黃玲玲(1986—),女,福建莆田人,特別研究員,博士生導師,北京市科技新星,2009年于天津大學、南開大學獲得雙學士學位,2014年于清華大學獲得博士學位,主要從事超穎表面物理機制與功能器件方面的研究。E-mail:huanglingling@bit.edu.cn

    王涌天(1957—),男,教授、博導,“長江學者獎勵計劃”特聘教授,主要從事頭盔立體顯示、增強現(xiàn)實虛實融合顯示、體三維和全息三維顯示及其相關光學系統(tǒng)設計方面的研究。E-mail:wyt@bit.edu.cn

    2017-05-11;

    2017-08-13

    國家自然科學基金資助項目(No.61505007);北京市科技新星計劃項目資助

    Supported by National Natural Science Foundation of China(No.61505007);Beijing Nova Program(No.Z171100001117047)

    2095-1531(2017)05-0523-18

    *Correspondingauthor,E-mail:huanglingling@bit.edu.cn

    猜你喜歡
    全息圖偏振光渦旋
    基于PM算法的渦旋電磁波引信超分辨測向方法
    部分偏振光分解的研究①
    光渦旋方程解的存在性研究
    線偏振光振動合成矢量端點軌跡方程的推導
    能觸摸的全息圖
    變截面復雜渦旋型線的加工幾何與力學仿真
    能觸摸的全息圖
    針刺結(jié)合偏振光疼痛治療儀治療三叉神經(jīng)痛30例
    抖動誤差對光學掃描全息圖像重建的影響
    河南科技(2014年8期)2014-02-27 14:07:49
    應該重視感生(渦旋)電場的方向性教學
    物理與工程(2013年6期)2013-03-11 16:06:19
    久久精品亚洲精品国产色婷小说| 中文资源天堂在线| 亚洲第一电影网av| 97超级碰碰碰精品色视频在线观看| 欧美黑人巨大hd| 亚洲av第一区精品v没综合| 亚洲av免费在线观看| 伊人久久大香线蕉亚洲五| 动漫黄色视频在线观看| 欧美乱码精品一区二区三区| 91在线观看av| 老汉色∧v一级毛片| 国产三级在线视频| 日韩欧美一区二区三区在线观看| or卡值多少钱| 国产精品,欧美在线| 成人三级黄色视频| 欧美+亚洲+日韩+国产| 一个人看视频在线观看www免费 | 男女下面进入的视频免费午夜| 午夜福利高清视频| 高潮久久久久久久久久久不卡| 搞女人的毛片| 亚洲中文字幕日韩| 国产精品爽爽va在线观看网站| 久久精品国产清高在天天线| 老司机福利观看| 十八禁网站免费在线| 免费看a级黄色片| 熟妇人妻久久中文字幕3abv| 高清在线国产一区| 国产欧美日韩精品一区二区| netflix在线观看网站| 亚洲在线自拍视频| 天堂av国产一区二区熟女人妻| av片东京热男人的天堂| 18美女黄网站色大片免费观看| 精品久久久久久久久久免费视频| 久9热在线精品视频| 国产成人欧美在线观看| 国产成人a区在线观看| 国产 一区 欧美 日韩| 精品福利观看| 男女午夜视频在线观看| 十八禁网站免费在线| or卡值多少钱| 精品一区二区三区视频在线 | 少妇熟女aⅴ在线视频| 国产在线精品亚洲第一网站| 无限看片的www在线观看| 脱女人内裤的视频| 免费av不卡在线播放| 亚洲专区中文字幕在线| 91麻豆精品激情在线观看国产| 19禁男女啪啪无遮挡网站| 免费观看人在逋| 日韩高清综合在线| 在线看三级毛片| 99热这里只有精品一区| 久久久久久久久中文| 国产在线精品亚洲第一网站| 少妇熟女aⅴ在线视频| 国产精品1区2区在线观看.| 99久久久亚洲精品蜜臀av| 999久久久精品免费观看国产| 亚洲一区二区三区色噜噜| 亚洲欧美日韩无卡精品| 18禁国产床啪视频网站| 性色av乱码一区二区三区2| 午夜免费激情av| 国产成人av教育| 成人18禁在线播放| 又黄又粗又硬又大视频| 91久久精品电影网| 久久中文看片网| 亚洲18禁久久av| 欧美午夜高清在线| 91久久精品国产一区二区成人 | 好男人在线观看高清免费视频| 无人区码免费观看不卡| 最新美女视频免费是黄的| 成人精品一区二区免费| 每晚都被弄得嗷嗷叫到高潮| 九九久久精品国产亚洲av麻豆| 国产亚洲欧美在线一区二区| 99热只有精品国产| 男人舔女人下体高潮全视频| 一个人免费在线观看的高清视频| 免费看美女性在线毛片视频| 欧美性猛交╳xxx乱大交人| av视频在线观看入口| 一进一出抽搐动态| 可以在线观看毛片的网站| 国语自产精品视频在线第100页| 两个人的视频大全免费| 国内揄拍国产精品人妻在线| 成年女人永久免费观看视频| 91久久精品国产一区二区成人 | 夜夜躁狠狠躁天天躁| 国产免费男女视频| 日韩精品青青久久久久久| 久久精品国产亚洲av香蕉五月| 中文字幕熟女人妻在线| 偷拍熟女少妇极品色| 首页视频小说图片口味搜索| 国产亚洲精品综合一区在线观看| 国产爱豆传媒在线观看| 国产爱豆传媒在线观看| 国产精品久久久人人做人人爽| 2021天堂中文幕一二区在线观| 成年女人看的毛片在线观看| 亚洲精品亚洲一区二区| 久久人人精品亚洲av| 一级a爱片免费观看的视频| 国产成年人精品一区二区| 国产免费av片在线观看野外av| 免费在线观看影片大全网站| 热99在线观看视频| 欧美一区二区亚洲| 女警被强在线播放| 午夜福利成人在线免费观看| 91在线观看av| 欧美又色又爽又黄视频| 国产精品久久久人人做人人爽| 天美传媒精品一区二区| 成年女人毛片免费观看观看9| 国内精品久久久久精免费| 91在线观看av| 国产三级中文精品| 麻豆国产97在线/欧美| 亚洲av日韩精品久久久久久密| 精品久久久久久久末码| 人妻久久中文字幕网| 亚洲国产中文字幕在线视频| 老司机在亚洲福利影院| 又粗又爽又猛毛片免费看| 日韩欧美 国产精品| 久久久成人免费电影| 国内毛片毛片毛片毛片毛片| 久久久久九九精品影院| 国产成人啪精品午夜网站| 免费一级毛片在线播放高清视频| 精品一区二区三区视频在线观看免费| 亚洲精品日韩av片在线观看 | 深夜精品福利| 国产淫片久久久久久久久 | www国产在线视频色| 一级黄片播放器| 成人午夜高清在线视频| 老司机午夜福利在线观看视频| 国产欧美日韩精品一区二区| 国产精品综合久久久久久久免费| 天天躁日日操中文字幕| 桃红色精品国产亚洲av| a在线观看视频网站| 婷婷丁香在线五月| 美女被艹到高潮喷水动态| 欧美日韩精品网址| 老熟妇仑乱视频hdxx| 精品国产美女av久久久久小说| 超碰av人人做人人爽久久 | 丁香六月欧美| 免费在线观看日本一区| 99精品在免费线老司机午夜| 亚洲专区国产一区二区| av天堂中文字幕网| 亚洲中文字幕日韩| 国产野战对白在线观看| 淫秽高清视频在线观看| 在线观看免费午夜福利视频| 国产探花在线观看一区二区| 老汉色∧v一级毛片| 国产伦精品一区二区三区视频9 | 高清毛片免费观看视频网站| 两个人视频免费观看高清| 成人av在线播放网站| 成熟少妇高潮喷水视频| 内射极品少妇av片p| 日韩免费av在线播放| 波多野结衣高清作品| 露出奶头的视频| 亚洲av第一区精品v没综合| 国产97色在线日韩免费| 久久精品国产综合久久久| 好男人电影高清在线观看| 天堂√8在线中文| 欧美黑人巨大hd| 给我免费播放毛片高清在线观看| 男女那种视频在线观看| 亚洲中文字幕日韩| 男女床上黄色一级片免费看| 久久久久九九精品影院| av黄色大香蕉| 久久国产精品影院| 99久久精品热视频| 69av精品久久久久久| 久久久久久久精品吃奶| www.熟女人妻精品国产| e午夜精品久久久久久久| 国产探花在线观看一区二区| 国产精品电影一区二区三区| 久久精品人妻少妇| 精品久久久久久久毛片微露脸| 欧美成人a在线观看| 一本久久中文字幕| 午夜激情福利司机影院| 啦啦啦观看免费观看视频高清| 成人特级黄色片久久久久久久| 热99在线观看视频| 国产一区二区三区在线臀色熟女| 老司机午夜十八禁免费视频| 国产欧美日韩一区二区精品| 亚洲av不卡在线观看| 好看av亚洲va欧美ⅴa在| 神马国产精品三级电影在线观看| 90打野战视频偷拍视频| 久久久久性生活片| 18+在线观看网站| 国产高清视频在线播放一区| 国产精品三级大全| 欧美性猛交黑人性爽| 成人性生交大片免费视频hd| 欧美最黄视频在线播放免费| 亚洲精品在线观看二区| 黄色视频,在线免费观看| 热99re8久久精品国产| 免费观看精品视频网站| 欧美又色又爽又黄视频| 制服丝袜大香蕉在线| 夜夜夜夜夜久久久久| 18+在线观看网站| 人人妻,人人澡人人爽秒播| 露出奶头的视频| 一进一出抽搐动态| 天堂av国产一区二区熟女人妻| 国产综合懂色| 波野结衣二区三区在线 | 亚洲电影在线观看av| 日本 av在线| 久久久久久大精品| 毛片女人毛片| av专区在线播放| 在线观看av片永久免费下载| a级毛片a级免费在线| 亚洲成a人片在线一区二区| 国产精品亚洲av一区麻豆| 99精品在免费线老司机午夜| 人人妻人人澡欧美一区二区| 天堂√8在线中文| 无人区码免费观看不卡| 老司机深夜福利视频在线观看| 国产真实伦视频高清在线观看 | 成人性生交大片免费视频hd| 欧美成人a在线观看| 亚洲第一电影网av| 悠悠久久av| 国产蜜桃级精品一区二区三区| 精品国产超薄肉色丝袜足j| 精品一区二区三区av网在线观看| 国产精品一区二区三区四区免费观看 | 欧美一区二区国产精品久久精品| 很黄的视频免费| 亚洲国产色片| 免费av观看视频| 国产男靠女视频免费网站| 精品免费久久久久久久清纯| xxx96com| 老汉色∧v一级毛片| 日韩免费av在线播放| 黄色视频,在线免费观看| 精品无人区乱码1区二区| 人人妻人人看人人澡| 久久精品国产综合久久久| 99热这里只有是精品50| 午夜福利在线观看免费完整高清在 | 久久九九热精品免费| 久久九九热精品免费| 色尼玛亚洲综合影院| 日本成人三级电影网站| 有码 亚洲区| 中文字幕精品亚洲无线码一区| 中文资源天堂在线| 757午夜福利合集在线观看| 亚洲专区国产一区二区| 国产一区二区在线av高清观看| 少妇人妻精品综合一区二区 | 别揉我奶头~嗯~啊~动态视频| 老司机午夜十八禁免费视频| 一个人看视频在线观看www免费 | 久久天躁狠狠躁夜夜2o2o| 精品久久久久久久久久久久久| 中文亚洲av片在线观看爽| 国产高清视频在线观看网站| 中文字幕精品亚洲无线码一区| 内射极品少妇av片p| 男人舔女人下体高潮全视频| 久久久色成人| 国产av一区在线观看免费| 日韩欧美在线乱码| 亚洲国产精品sss在线观看| 最近视频中文字幕2019在线8| 夜夜夜夜夜久久久久| 757午夜福利合集在线观看| 国产私拍福利视频在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲午夜理论影院| 色在线成人网| 成人无遮挡网站| 欧美黄色淫秽网站| 亚洲中文字幕日韩| 91麻豆精品激情在线观看国产| 三级毛片av免费| 一本综合久久免费| 男女床上黄色一级片免费看| 中文字幕精品亚洲无线码一区| 老司机深夜福利视频在线观看| 午夜福利18| 欧美成人a在线观看| 亚洲av成人精品一区久久| 国产精品久久久久久久电影 | 亚洲成人久久爱视频| 国产亚洲精品久久久久久毛片| 免费高清视频大片| 国产成人啪精品午夜网站| 精品一区二区三区视频在线观看免费| 岛国在线免费视频观看| 亚洲天堂国产精品一区在线| 国产高潮美女av| 精品久久久久久久毛片微露脸| 制服丝袜大香蕉在线| www国产在线视频色| 一个人看的www免费观看视频| 亚洲av免费在线观看| 精品无人区乱码1区二区| 免费无遮挡裸体视频| 免费高清视频大片| 国产精品嫩草影院av在线观看 | 18禁在线播放成人免费| 亚洲久久久久久中文字幕| 婷婷亚洲欧美| 国产精品久久久久久精品电影| 一区二区三区免费毛片| 国产在线精品亚洲第一网站| 99国产极品粉嫩在线观看| 欧美黑人巨大hd| 国产欧美日韩精品一区二区| 午夜福利18| 久久这里只有精品中国| 久久6这里有精品| 日本一本二区三区精品| 午夜影院日韩av| 最近最新免费中文字幕在线| av专区在线播放| 一级a爱片免费观看的视频| 午夜免费激情av| 一本综合久久免费| 97人妻精品一区二区三区麻豆| 老鸭窝网址在线观看| 国产精品一区二区免费欧美| 搡老岳熟女国产| 最近在线观看免费完整版| 一个人看视频在线观看www免费 | 精品无人区乱码1区二区| 特级一级黄色大片| 国产精品一区二区三区四区久久| 久久精品国产99精品国产亚洲性色| 熟女电影av网| 99久久九九国产精品国产免费| 国内毛片毛片毛片毛片毛片| 国产av麻豆久久久久久久| 91久久精品电影网| 91字幕亚洲| 久久久久久国产a免费观看| 国产69精品久久久久777片| 亚洲国产高清在线一区二区三| 国产精品精品国产色婷婷| 尤物成人国产欧美一区二区三区| 亚洲av电影在线进入| 精品国内亚洲2022精品成人| 97人妻精品一区二区三区麻豆| 精品欧美国产一区二区三| 久久精品亚洲精品国产色婷小说| 久久久久久久午夜电影| 天堂影院成人在线观看| 欧美日本亚洲视频在线播放| 亚洲欧美日韩高清专用| 99国产精品一区二区三区| 波多野结衣高清无吗| www.www免费av| 国产伦在线观看视频一区| 丰满人妻一区二区三区视频av | 免费观看精品视频网站| 成人三级黄色视频| 99国产综合亚洲精品| 最新美女视频免费是黄的| 国产成人系列免费观看| 操出白浆在线播放| 国产精品1区2区在线观看.| 亚洲无线在线观看| 亚洲欧美精品综合久久99| 国产午夜精品论理片| ponron亚洲| 国产精品综合久久久久久久免费| 亚洲五月天丁香| 国产伦精品一区二区三区四那| 久久亚洲精品不卡| xxxwww97欧美| 在线天堂最新版资源| 久久久精品欧美日韩精品| 亚洲不卡免费看| 亚洲 国产 在线| 一个人免费在线观看电影| 成年女人毛片免费观看观看9| 国产97色在线日韩免费| 少妇高潮的动态图| 哪里可以看免费的av片| 99热这里只有精品一区| 亚洲自拍偷在线| 1024手机看黄色片| 欧美不卡视频在线免费观看| 黄色视频,在线免费观看| 天天躁日日操中文字幕| 观看免费一级毛片| 五月伊人婷婷丁香| 国产日本99.免费观看| 三级国产精品欧美在线观看| 亚洲精品色激情综合| 国产精品久久电影中文字幕| 久久国产乱子伦精品免费另类| 男女视频在线观看网站免费| 亚洲avbb在线观看| 又紧又爽又黄一区二区| 久久久精品欧美日韩精品| 夜夜躁狠狠躁天天躁| avwww免费| 日本撒尿小便嘘嘘汇集6| 午夜福利在线观看免费完整高清在 | 色综合婷婷激情| 深夜精品福利| 欧美乱妇无乱码| 悠悠久久av| 亚洲aⅴ乱码一区二区在线播放| 亚洲av五月六月丁香网| 毛片女人毛片| 亚洲国产高清在线一区二区三| 97超视频在线观看视频| 亚洲色图av天堂| 国产精品av视频在线免费观看| 免费搜索国产男女视频| 成年免费大片在线观看| 国产熟女xx| 国产成人系列免费观看| 搡老妇女老女人老熟妇| 色综合欧美亚洲国产小说| 国产精品国产高清国产av| 美女大奶头视频| 母亲3免费完整高清在线观看| 一卡2卡三卡四卡精品乱码亚洲| 9191精品国产免费久久| 亚洲欧美日韩东京热| 亚洲美女视频黄频| 欧美成人一区二区免费高清观看| 欧美大码av| 国产午夜福利久久久久久| 亚洲乱码一区二区免费版| 久久精品综合一区二区三区| 亚洲国产色片| 免费人成在线观看视频色| 少妇丰满av| 中文字幕精品亚洲无线码一区| 成人特级黄色片久久久久久久| 超碰av人人做人人爽久久 | 很黄的视频免费| 色老头精品视频在线观看| 十八禁网站免费在线| 日韩欧美在线二视频| 日本与韩国留学比较| 亚洲av免费高清在线观看| 老司机深夜福利视频在线观看| 久久久久久国产a免费观看| 国产 一区 欧美 日韩| www日本黄色视频网| 很黄的视频免费| 国产高潮美女av| 国模一区二区三区四区视频| 国产黄片美女视频| 亚洲国产精品成人综合色| 欧美日本视频| 国产精品久久久人人做人人爽| 99久国产av精品| 亚洲 欧美 日韩 在线 免费| 亚洲精华国产精华精| 黑人欧美特级aaaaaa片| 日本黄大片高清| 日韩欧美在线二视频| 丁香六月欧美| 国产高清视频在线观看网站| 午夜久久久久精精品| 久久人人精品亚洲av| 99热6这里只有精品| 久久精品人妻少妇| 观看美女的网站| 久久精品国产综合久久久| 亚洲精品美女久久久久99蜜臀| 亚洲精品乱码久久久v下载方式 | 99久久久亚洲精品蜜臀av| 国模一区二区三区四区视频| 精品福利观看| 天天躁日日操中文字幕| 国产精品国产高清国产av| 一夜夜www| 韩国av一区二区三区四区| 国产欧美日韩一区二区三| 欧美极品一区二区三区四区| 成人av一区二区三区在线看| 精品熟女少妇八av免费久了| 国产一区在线观看成人免费| 两个人看的免费小视频| АⅤ资源中文在线天堂| 欧美日韩国产亚洲二区| 母亲3免费完整高清在线观看| 一进一出抽搐gif免费好疼| 最新在线观看一区二区三区| 国产精品免费一区二区三区在线| 91九色精品人成在线观看| 日韩欧美精品免费久久 | 啦啦啦免费观看视频1| 麻豆国产av国片精品| 男人舔奶头视频| 亚洲第一电影网av| 噜噜噜噜噜久久久久久91| 欧美不卡视频在线免费观看| 熟女人妻精品中文字幕| 欧美日本亚洲视频在线播放| 国产 一区 欧美 日韩| 国产精品自产拍在线观看55亚洲| 在线观看免费视频日本深夜| 97超级碰碰碰精品色视频在线观看| 免费av毛片视频| netflix在线观看网站| 久久久久国内视频| 人人妻人人澡欧美一区二区| 国产aⅴ精品一区二区三区波| av在线天堂中文字幕| 日韩中文字幕欧美一区二区| 久久久国产成人精品二区| 国产精品一及| 无限看片的www在线观看| 亚洲最大成人手机在线| 18禁黄网站禁片免费观看直播| 国内精品久久久久精免费| 日日干狠狠操夜夜爽| 999久久久精品免费观看国产| 免费看日本二区| 国产欧美日韩精品亚洲av| 色老头精品视频在线观看| 久99久视频精品免费| 久久这里只有精品中国| 免费av观看视频| 国产成人影院久久av| 午夜免费成人在线视频| 天堂影院成人在线观看| 在线观看美女被高潮喷水网站 | 欧美日本亚洲视频在线播放| 99热这里只有是精品50| 狠狠狠狠99中文字幕| ponron亚洲| 国产精品永久免费网站| 国产一区二区激情短视频| 黄色成人免费大全| 精品无人区乱码1区二区| 免费电影在线观看免费观看| 内射极品少妇av片p| 99在线视频只有这里精品首页| 在线观看一区二区三区| av专区在线播放| 一区二区三区免费毛片| 亚洲乱码一区二区免费版| 淫妇啪啪啪对白视频| 欧美激情在线99| 日韩精品青青久久久久久| 国产综合懂色| 狠狠狠狠99中文字幕| 亚洲专区中文字幕在线| 日韩欧美国产在线观看| 亚洲av第一区精品v没综合| 欧美一区二区精品小视频在线| 欧美另类亚洲清纯唯美| 国产精品爽爽va在线观看网站| 日韩中文字幕欧美一区二区| 精品午夜福利视频在线观看一区| 国产精品日韩av在线免费观看| 国产老妇女一区| 国产成人欧美在线观看| 亚洲片人在线观看| 国产激情偷乱视频一区二区| 中亚洲国语对白在线视频| 亚洲美女黄片视频| 免费人成在线观看视频色| 午夜福利在线观看吧| 少妇人妻一区二区三区视频| 免费一级毛片在线播放高清视频| 性色avwww在线观看| 一级a爱片免费观看的视频| 一级毛片高清免费大全| 久久久久国内视频| 精华霜和精华液先用哪个| 一区二区三区国产精品乱码| 美女cb高潮喷水在线观看| 在线观看一区二区三区| 91在线观看av| 精品一区二区三区视频在线 | 天堂av国产一区二区熟女人妻| 国产成人影院久久av| 久久精品国产自在天天线|