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

    類石墨烯復(fù)雜晶胞光子晶體中的確定性界面態(tài)?

    2017-12-05 02:35:00賈子源楊玉婷季立宇杭志宏
    物理學(xué)報 2017年22期
    關(guān)鍵詞:性質(zhì)界面結(jié)構(gòu)

    賈子源 楊玉婷 季立宇 杭志宏

    (蘇州大學(xué)物理與光電.能源學(xué)部,蘇州納米科技協(xié)同創(chuàng)新中心,蘇州 215006)

    類石墨烯復(fù)雜晶胞光子晶體中的確定性界面態(tài)?

    賈子源 楊玉婷 季立宇 杭志宏?

    (蘇州大學(xué)物理與光電.能源學(xué)部,蘇州納米科技協(xié)同創(chuàng)新中心,蘇州 215006)

    (2017年7月24日收到;2017年8月13日收到修改稿)

    拓?fù)浣^緣體是當(dāng)前凝聚態(tài)物理領(lǐng)域研究的熱點(diǎn)問題.利用石墨烯材料的特殊能帶特性來實(shí)現(xiàn)拓?fù)漭斶\(yùn)特性在設(shè)計下一代電子和能谷電子器件方面具有較廣泛的應(yīng)用前景.基于光子與電子的類比,利用光子拓?fù)洳牧蠈?shí)現(xiàn)了確定性界面態(tài);構(gòu)建了具有C6v對稱性的類似石墨烯結(jié)構(gòu)的的光子晶體復(fù)雜晶格;通過多種方式降低晶格對稱性來獲得具有C3v,C3,C2v和C2對稱的晶體,從而打破能谷簡并實(shí)現(xiàn)全光子帶隙結(jié)構(gòu);將體拓?fù)湫再|(zhì)不同的兩種光子晶體擺放在一起,在此具有反轉(zhuǎn)體能帶性質(zhì)的界面上,實(shí)現(xiàn)了具有單向傳輸特性的拓?fù)浯_定性界面態(tài)的傳輸.利用光子晶體結(jié)構(gòu)的容易加工性,可以簡便地調(diào)控拓?fù)浣缑鎽B(tài)控制光的傳播,可為未來光拓?fù)浣^緣體的研究提供良好的平臺.

    光子晶體,拓?fù)?谷自由度,界面態(tài)

    1 引 言

    拓?fù)鋵W(xué)是近代發(fā)展起來的數(shù)學(xué)概念,而隨著凝聚態(tài)物理的發(fā)展,人們發(fā)現(xiàn)了一類具有相變、非平庸特性的新型材料——拓?fù)浣^緣體.早在1980年,Klitzing發(fā)現(xiàn)了整數(shù)量子霍爾效應(yīng)[1],此系統(tǒng)的特性完全符合拓?fù)浣^緣體的概念.隨后Thouless-Kohmoto-Nightingale-den Nijs關(guān)系[2]指出其中量子化的根源是由于拓?fù)洳蛔兞康拇嬖?更是證明了量子化整數(shù)對應(yīng)的是非零陳數(shù).在此系統(tǒng)中出現(xiàn)了一種新的物態(tài)拓?fù)鋺B(tài),在其中電子輸運(yùn)的表現(xiàn)為完全的背散射抑制和單向傳輸.不同于需要外加磁場破壞時間反演對稱性來實(shí)現(xiàn)整數(shù)量子霍爾效應(yīng),Kane和Mele[3]最早提出在低能量級的單層石墨烯中存在自旋軌道耦合導(dǎo)致的量子自旋霍爾效應(yīng),即在石墨烯材料邊界存在一對共軛且自旋方向相反的無能隙拓?fù)溥吔鐟B(tài).我們已經(jīng)知道以上觀測到的量子化霍爾效應(yīng)的根本原因是材料的拓?fù)涮匦?與具體波動方程無關(guān).從布洛赫定理和能帶性質(zhì)來看,光子可與電子的性質(zhì)相比擬,那么光學(xué)領(lǐng)域也可能實(shí)現(xiàn)類似的光學(xué)拓?fù)浣^緣體.

    石墨烯是從石墨中剝離出來的、單層碳原子構(gòu)成的二維晶體結(jié)構(gòu)[4].因其優(yōu)越的導(dǎo)電導(dǎo)熱高透明度等性能,石墨烯廣泛地應(yīng)用在很多領(lǐng)域.而且石墨烯具有很多同素異形體結(jié)構(gòu),如具有更加復(fù)雜晶胞結(jié)構(gòu)的石墨炔材料.受石墨烯材料具有的許多新奇物理性質(zhì)的啟發(fā),人們開始研究與石墨烯具有相同對稱性、相似能帶結(jié)構(gòu)的光子晶體(photonic crystal,PC).構(gòu)成光子晶體的介電材料具有不同介電常數(shù),它在空間中按一定周期排列[5,6].已經(jīng)證明可以利用光子晶體來操縱光的傳播行為[7,8],如光子晶體波導(dǎo)[9,10].隨著光子晶體能帶理論的研究,尤其是設(shè)計和實(shí)驗(yàn)方法的發(fā)展,將拓?fù)涞母拍钜牍鈱W(xué)系統(tǒng)特別是用光子晶體實(shí)現(xiàn)拓?fù)浣^緣體,將會不斷拓寬人們構(gòu)造新型元器件的方法和思路.

    迄今已經(jīng)開展了許多波導(dǎo)傳輸性質(zhì)的研究工作.其中一方面是利用光子晶體界面態(tài)去控制光的傳輸[11].已經(jīng)研究的構(gòu)建光子晶體界面態(tài)的方法有:構(gòu)建光子晶體異質(zhì)結(jié)構(gòu)[12,13]、修飾光子晶體表面[14,15]、堆積擁有共同帶隙的光子晶體等[16].但是這些方法并不能保證界面態(tài)一定會產(chǎn)生,需要仔細(xì)地調(diào)節(jié)結(jié)構(gòu)才會產(chǎn)生界面態(tài).不同于之前的研究成果,香港科技大學(xué)Chan課題組[17,18]利用一種新的物理機(jī)理在一維和二維光子晶體中實(shí)現(xiàn)了確定性界面態(tài),即通過操控光子晶體結(jié)構(gòu)參數(shù),簡并的狄拉克錐色散會破缺,從而出現(xiàn)帶隙,能帶性質(zhì)發(fā)生拓?fù)湎嘧?構(gòu)成界面的這兩種光子晶體的體能帶拓?fù)湫再|(zhì)反轉(zhuǎn),則在兩者共同帶隙處就一定會存在界面態(tài).而且微波段實(shí)驗(yàn)已經(jīng)在二維正方晶格光子晶體中證明了確定性界面態(tài)的存在[19,20].并且還可以同時在多個頻段實(shí)現(xiàn)界面態(tài),基于此性質(zhì)可以用作多帶波導(dǎo)的信息傳輸,甚至是將應(yīng)用拓展于光頻段實(shí)現(xiàn)波分復(fù)用器[21].

    然而以往的研究工作中,主要研究確定性界面態(tài)存在于光子晶體簡單晶格結(jié)構(gòu)中如正方晶格和三角晶格,卻沒有關(guān)于復(fù)雜晶胞的探討.在本文中,類似于石墨烯材料的復(fù)雜晶胞結(jié)構(gòu),我們構(gòu)建了具有復(fù)雜晶格的光子晶體,嘗試研究復(fù)雜晶胞所具備的多種對稱性破缺的可能性.比之類石墨烯的的碳同素異形體在實(shí)驗(yàn)上極難獲取,我們嘗試構(gòu)造的復(fù)雜晶胞光子晶體在實(shí)驗(yàn)上便于實(shí)現(xiàn),而且可以通過調(diào)節(jié)結(jié)構(gòu)參數(shù)和對稱性來操縱色散關(guān)系.我們通過四種方式打破光子晶體結(jié)構(gòu)的C6v對稱性,破缺能谷簡并的狄拉克點(diǎn)構(gòu)建具有反轉(zhuǎn)體能帶性質(zhì)的兩個光子晶體;基于界面態(tài)存在的拓?fù)錂C(jī)理,在光子晶體復(fù)雜晶胞中實(shí)現(xiàn)確定性界面態(tài).本文的研究工作拓寬了操控光子晶體能帶實(shí)現(xiàn)拓?fù)浣缑鎽B(tài)的思路和方法,這些方法簡單易行,而且可以為光波導(dǎo)應(yīng)用方面提供更廣闊的平臺.

    2 類石墨烯復(fù)雜晶胞

    人們對電子的兩個內(nèi)稟自由度電荷和自旋的研究促進(jìn)了自旋電子學(xué)的成熟發(fā)展和廣泛應(yīng)用.隨著石墨烯材料的發(fā)現(xiàn)和研究,科研工作者操控電子新的自由度能谷,則促進(jìn)了能谷電子學(xué)(valleytronics)的理論和應(yīng)用探索[22,23].能谷電子學(xué)旨在利用能谷作為載體在現(xiàn)代器件中加工處理信息.許多有趣的物理性質(zhì)已經(jīng)在理論上得到證明并在實(shí)驗(yàn)中實(shí)現(xiàn),如谷分離(valley filter)[24]、谷霍爾效應(yīng)(valley Hall e ff ect)[25,26].石墨烯材料具有蜂窩狀晶格,一個原胞結(jié)構(gòu)由兩個碳原子構(gòu)成,具有C6v對稱性.同時石墨烯具有許多同素異形體,即具有更加復(fù)雜多樣的晶胞結(jié)構(gòu)[27].類石墨烯復(fù)雜晶胞就是用若干碳原子簇代替原胞上的一個原子.石墨烯材料優(yōu)異的應(yīng)用特性源于其特殊的能帶性質(zhì).在動量空間中存在的一對不等價的能量局部極值叫做能谷,表現(xiàn)為在第一布里淵區(qū)六個角處(K和K′點(diǎn))石墨烯的能帶簡并在一起,簡并點(diǎn)稱為狄拉克點(diǎn)(Dirac point).K和K′點(diǎn)的能谷是不等價的,其波函數(shù)滿足時間反演對稱性.兩者作為能谷自由度的兩個態(tài)可以像電子自旋向上、向下態(tài)一樣作為信息輸運(yùn)的載體,被稱作谷贗自旋.我們同樣考慮二維光子晶體.蜂巢晶格中的每個碳原子被一根無限高的介質(zhì)柱取代,此種光子晶體也具有前述的類似能帶性質(zhì);再用三個介質(zhì)圓柱替代原本的一個圓柱構(gòu)建類石墨烯的復(fù)雜晶胞結(jié)構(gòu).對比于石墨烯結(jié)構(gòu)在可調(diào)控性方面的局限性,此復(fù)雜晶胞光子晶體在對稱性方面具有更多可利用的自由度,并且更加容易操縱產(chǎn)生界面態(tài)色散關(guān)系.

    二維光子晶體中確定性界面態(tài)的產(chǎn)生與界面通道兩邊結(jié)構(gòu)的不同拓?fù)湫再|(zhì)有關(guān).具體體現(xiàn)為如果構(gòu)建此界面的兩側(cè)光子晶體體能帶性質(zhì)反轉(zhuǎn),那么就一定存在確定性界面態(tài).以最簡單的二維正方晶格光子晶體為例(如圖1(a)),通過調(diào)控有效介電常數(shù)、晶格常數(shù)和介質(zhì)柱半徑,可以讓能帶在布里淵區(qū)中心Γ點(diǎn)發(fā)生簡并,這就是通常所說的偶然簡并[28].簡并的能帶色散滿足線性關(guān)系,它和電子體系中的狄拉克方程具有對應(yīng)關(guān)系.如此便在光子晶體系統(tǒng)中實(shí)現(xiàn)了類似石墨烯結(jié)構(gòu)中的狄拉克錐.通過繼續(xù)調(diào)節(jié)結(jié)構(gòu)參數(shù)則可以打開簡并的能帶獲得帶隙.如圖1(b)中,在光子晶體1(PC1)的能帶示意圖中,三條能帶在在Γ點(diǎn)的本征模式分布依次為類似電子系統(tǒng)的s,py,px軌道.然而當(dāng)我們將結(jié)構(gòu)調(diào)整為光子晶體2(PC2),它的體能帶結(jié)構(gòu)則會與PC1發(fā)生反轉(zhuǎn),如圖1(c)所示.三條能帶中下面簡并的兩條在Γ點(diǎn)為本征模式px,py,而最上面一條為s模式.因此PC1與PC2能帶性質(zhì)反轉(zhuǎn),在兩者反轉(zhuǎn)能帶共同帶隙中則一定存在界面態(tài).此界面態(tài)的存在與否與光子晶體的體性質(zhì)有關(guān),只需將不同體性質(zhì)光子晶體擺放一起就一定存在,因此稱為確定性界面態(tài).圖1中的能帶結(jié)構(gòu)均為示意圖,構(gòu)成此系統(tǒng)的PC1與PC2的參數(shù)以及界面態(tài)具體傳播性質(zhì)可以參考文獻(xiàn)[19],這里不再詳述.

    圖1 (a)光子晶體在布里淵區(qū)中心的狄拉克錐色散;(b),(c)兩個體能帶反轉(zhuǎn)的光子晶體Fig.1.(a)Dirac cone dispersion of PC around the center of Brillouin zone center;(b),(c)two inverted band structure of PC1 and PC2.

    在本文工作中,基于確定性界面態(tài)存在的物理機(jī)理,構(gòu)建擁有更多自由度的復(fù)雜晶胞光子晶體,并實(shí)現(xiàn)拓?fù)浣缑鎽B(tài).圖2(a)為復(fù)雜結(jié)構(gòu)的蜂窩狀光子晶體示意圖,深藍(lán)色實(shí)心圓表示介質(zhì)柱,淺藍(lán)色的區(qū)域?yàn)楸尘安牧?在紅色實(shí)線標(biāo)示的六角晶胞中可以清晰地看到三個介質(zhì)圓柱(用綠色凸顯)簇?fù)碓谝黄饦?gòu)成正三角形(綠色實(shí)線),依然具有C6v對稱性,此即為三個介質(zhì)柱代替一個介質(zhì)柱構(gòu)成的復(fù)雜晶胞光子晶體.其中相對介電常數(shù)ε=7.5的介質(zhì)圓柱嵌在背景為空氣的材料中.圖中紅色菱形虛線框內(nèi)的原胞結(jié)構(gòu),晶格常數(shù)為a,介質(zhì)柱的半徑r=0.1a,到三角形中心的距離為h=0.16a.用COMSOL Multiphysics仿真軟件計算此結(jié)構(gòu)的TM(transverse magnetic)模式能帶,即電場Ez沿著介質(zhì)柱z軸方向.計算的能帶圖見圖2(b),第一和第二條能帶在布里淵區(qū)(插圖中六角形示意圖)中K和K′能谷簡并,簡并點(diǎn)歸一化頻率為0.362c/a(c為真空中光速).由此可見雖然原胞結(jié)構(gòu)變復(fù)雜,但是保持對稱性不變光子晶體的能帶性質(zhì)依然不變.

    圖2 (a)復(fù)雜光子晶體結(jié)構(gòu)示意圖;(b)光子晶體能帶結(jié)構(gòu)Fig.2.(a)The schematic of the complex lattice PC;(b)PC’s band diagram.

    3 C3v對稱性光子晶體中的確定性界面態(tài)

    通過打破光子晶體的C6v對稱性,可以很容易破缺K和K′點(diǎn)的能谷簡并從而形成帶隙.以圖2原胞中三個介質(zhì)柱構(gòu)成的正三角形(灰色區(qū)域)中心為基點(diǎn),在分別對每個三角形拉伸和壓縮柱子之間的距離,調(diào)控后得到結(jié)構(gòu)A,如圖3(a)所示.光子晶體A原胞中左下角三個介質(zhì)圓柱拉伸后到中心的距離為h1=0.18a,而右上角三個柱子壓縮后到中心的距離為h2=0.12a.通過這種方式實(shí)現(xiàn)了具有C3v對稱性的光子晶體.從圖3(b)中可以看到在第一條和第二條能帶中間出現(xiàn)完全帶隙,原本在布里淵邊界K點(diǎn)的簡并變成帶隙的下限K1和上限K2點(diǎn).同時K1和K2點(diǎn)的性質(zhì)是不同的,如圖3(c)的相位分布.兩者位相分布的手征性是相反的,前者是沿逆時針方向(紫色箭頭所示)相位遞減而后者是順時針方向(綠箭頭所示)相位遞減,旋轉(zhuǎn)一周回復(fù)初始相位,對應(yīng)角動量m=1和m=?1.而且相位奇點(diǎn)的位置也不同,分別位于拉伸的柱子中心(K1)和壓縮的柱子中心(K2).將拉伸和壓縮柱子到中心距離的方式互換,得到光子晶體結(jié)構(gòu)B,能帶如圖4(a),它的體能帶性質(zhì)與結(jié)構(gòu)A的能帶是反轉(zhuǎn)的.可以從圖4(b)K點(diǎn)的相位手征性上很清晰地看到,K1(K2)點(diǎn)相位以順(逆)時針方向遞減.

    圖3 (a)兩個相反結(jié)構(gòu)C3v對稱的光子晶體結(jié)構(gòu)示意圖;(b)光子晶體結(jié)構(gòu)A全帶隙能帶;(c)在K1和K2點(diǎn)相反手征性的相位分布Fig.3.(a)The schematic of two inverted PC structures with C3vsymmetry;(b)band diagram with complete band gap for photonic crystal structure A;(d),(e)the inverted chirality of phase distributions at K1 and K2points.

    圖4 (a)光子晶體結(jié)構(gòu)B的全帶隙能帶;(b)在K1和K2點(diǎn)相反手征性的相位分布Fig.4.(a)Band diagram with complete band gap for photonic crystal structure B;(b)the inverted chirality of phase distributions at K1and K2points.

    可以看到光子晶體A和B的體能帶性質(zhì)類似,在布里淵區(qū)邊界K點(diǎn)附近有完全帶隙,但此帶隙的上下帶邊的能帶手性正好反轉(zhuǎn),即在此帶隙光子晶體A和B的體拓?fù)湫再|(zhì)相反.可以預(yù)期,如果將光子晶體A和B擺放在一起,在兩者的界面就會有拓?fù)浯_定性界面態(tài).我們用超原胞的方法計算了光子晶體A和B沿著kx方向的投影能帶來找尋此拓?fù)浯_定性界面態(tài).計算時在界面兩邊各用10個光子晶體原胞構(gòu)成一個超原胞,周期性邊界條件應(yīng)用于光子晶體x和y方向.能帶結(jié)構(gòu)如圖5(a)所示,灰色部分是光子晶體的體能帶,綠色和紫色線部分表示界面態(tài)的色散關(guān)系.能帶一部分是在兩個光子晶體共同帶隙處而另有一部分延伸到體帶里.從色散曲線的斜率可以看出在界面處存在兩個相反方向傳輸?shù)耐負(fù)鋺B(tài).而且兩個態(tài)的自旋極化方向相反,分別為自旋向上(紫色線)和自旋向下(綠色線).為了更形象地表征界面態(tài)的本征模態(tài),圖5(b)展示了頻率為0.384c/a同時波矢kx= ?0.38(2π/a)的Ez本征場分布(圖5(a)中黑色五角星處).從圖中可以很明顯地看到場局域在兩個光子晶體界面處,而向兩邊衰減很快.

    利用光子晶體拓?fù)鋺B(tài)的傳輸性質(zhì),可以設(shè)計并實(shí)現(xiàn)界面態(tài)傳輸波導(dǎo).我們用有限元模擬的方法計算光子晶體界面態(tài)的Ez電場分布,如圖6.模擬的上下區(qū)域分別由光子晶體結(jié)構(gòu)A和B構(gòu)成,并包含24×16個復(fù)合晶胞結(jié)構(gòu).區(qū)域四周的邊界條件用完美匹配層,這樣可以避免來自邊界的散射.需要用具有手征性的光學(xué)激發(fā)源才能激發(fā)單向的界面態(tài).我們用四個點(diǎn)源陣列(2×2)構(gòu)建手征性激發(fā)源[29].其中,四個點(diǎn)源發(fā)射的TM波波函數(shù)的初始相位按照不同旋轉(zhuǎn)方向依次遞減.初始相位按照逆(順)時針方向遞減的激發(fā)源和光學(xué)渦旋指數(shù)即拓?fù)潆姾蓴?shù)m=1(m=?1)的電場分布一致.將拓?fù)潆姾蓴?shù)為m=1的手征性源放在光子晶體A和B界面中間,位置如圖6中白色圓所示,并且光子晶體A置于結(jié)構(gòu)B的下方.在源發(fā)射頻率為0.384c/a時激發(fā)的單向傳輸界面態(tài)場分布如圖6所示,可以觀察到激發(fā)的拓?fù)浣缑鎽B(tài)沿著光子晶體結(jié)構(gòu)A和B的界面通道向左(圖中白色箭頭方向)單向傳輸.從圖5的投影能帶中可以看到,在頻率為0.384c/a時,界面態(tài)色散能帶的頻率和體能帶會有部分重疊.所以在此時,界面在傳播的過程中會激發(fā)一些體的電場分布.但是由于這些體電場分布相較于界面局域場非常弱,因而并不能對拓?fù)鋺B(tài)的傳輸與應(yīng)用造成很大的影響.如果用拓?fù)潆姾蔀閙=?1的源則可以激發(fā)向右傳播的拓?fù)鋺B(tài).在兩個光子晶體界面通道傳輸?shù)慕缑鎽B(tài)是受拓?fù)浔Wo(hù)的.界面態(tài)具有傳輸單向特性且不受傳播路徑上可能存在的雜質(zhì)的影響.即使光子晶體材料有雜質(zhì)不純凈或者是加工制備不完美的情況,拓?fù)鋺B(tài)依然能夠很好地在界面?zhèn)鬏敹粫a(chǎn)生背散射,因此能夠?qū)崿F(xiàn)拓?fù)湫畔⑤斶\(yùn)從而制備光學(xué)器件.這在光學(xué)波導(dǎo)方面將具有很大的應(yīng)用前景.

    圖5 (a)光子晶體投影能帶;(b)界面態(tài)Ez本征場分布Fig.5.(a)The projected band structure of two PCs;(b)the Ezeigen field distribution of interface state.

    圖6 模擬的光子晶體單向傳輸界面態(tài)Ez電場分布Fig.6.The simulated Ez field distribution of topological interface state with unidirectional propagation for two PCs.

    4 C3對稱性光子晶體中的確定性界面態(tài)

    在原子尺度上操控石墨烯實(shí)現(xiàn)拓?fù)浣缑鎽B(tài)對技術(shù)和條件有一定要求.而人工介質(zhì)材料微結(jié)構(gòu)單元構(gòu)建的光子晶體則更容易操作和可調(diào).并且光子晶體可調(diào)控的參數(shù)有更多的自由度,比如材料的介電常數(shù)、晶格常數(shù)、占空比等.通過調(diào)控結(jié)構(gòu)還可獲得一種C3對稱性的光子晶體,即分別將圖3中的結(jié)構(gòu)A和B圍繞三角形中心逆時針旋轉(zhuǎn)角度θ=15?,如圖7(a)所示.通過這種打破晶格對稱性的方式依然可以破缺能帶結(jié)構(gòu)中K和K′狄拉克點(diǎn)簡并.圖7(b)和圖7(c)中兩個能帶在K點(diǎn)其性質(zhì)反轉(zhuǎn).調(diào)節(jié)的兩個結(jié)構(gòu)之間有共同帶隙,經(jīng)過投影能帶的計算(計算方式同圖5),在兩個光子晶體界面處依然存在確定性界面態(tài),如圖7(d)所示.圖7(e)中五角星標(biāo)示的是在歸一化頻率為0.384c/a、波矢kx= ?0.4(2π/a)時的界面態(tài)Ez本征場分布.利用手征性源激發(fā),基于此結(jié)構(gòu)的界面態(tài)色散性質(zhì),依然可以像圖6一樣設(shè)計單向拓?fù)鋺B(tài)傳輸波導(dǎo),從而實(shí)現(xiàn)信息的輸運(yùn)與傳播.但通過調(diào)節(jié)介質(zhì)柱簇的旋轉(zhuǎn)角度,可以調(diào)控此拓?fù)浣缑鎽B(tài)的色散關(guān)系.

    圖7 (a)具有C3對稱的兩個光子晶體結(jié)構(gòu)示意圖;(b),(c)光子晶體能帶結(jié)構(gòu);(d)投影能帶;(e)界面態(tài)Ez本征場分布Fig.7.(a)The schematics of two PCs with C3symmetry;(b),(c)band diagrams of PCs;(d)the projected band structure;(e)the Ezeigen field distribution of interface state.

    5 C2v和C2對稱性光子晶體中分別實(shí)現(xiàn)確定性界面態(tài)

    圖8 (a)具有C2v對稱性的兩個光子晶體C和D示意圖;(b)光子晶體C能帶結(jié)構(gòu);(c)K1和K2點(diǎn)本征場分布Fig.8.(a)The schematics of two PCs with C2vsymmetry;(b)band diagram of PC with structure C;(c)the eigen field distribution at K1and K2points.

    進(jìn)一步打破光子晶體對稱性,還可以用具有C2v對稱的光子晶體中實(shí)現(xiàn)確定性界面態(tài).結(jié)構(gòu)如圖8(a)所示,光子晶體介質(zhì)柱半徑r=0.1a,相對介電常數(shù)ε=7.5.將原胞中原本構(gòu)成正三角形底邊處的兩個介質(zhì)圓柱,沿y方向(箭頭所示方向)增加距離d1=0.104a,而三角形頂點(diǎn)處圓柱依然保持到中心的距離h=0.16a不變,三個介質(zhì)柱重新形成等腰三角形,從而得到結(jié)構(gòu)C.結(jié)構(gòu)C的能帶結(jié)構(gòu)如圖8(b),原本K點(diǎn)簡并的狄拉克點(diǎn)被打破因而出現(xiàn)完全帶隙.與之相反的調(diào)控得到結(jié)構(gòu)D,即只沿y方向減小底邊處的介質(zhì)柱的距離d2=0.08a,仍然可以破缺狄拉克點(diǎn),如圖9(a)所示能帶結(jié)構(gòu).光子晶體C和D的拓?fù)湫再|(zhì)依然是不同的,兩者能帶反轉(zhuǎn).從圖8和圖9(c)中所示的K1和K2點(diǎn)電場振幅本征場分布可以很清楚地看到,在這兩點(diǎn)處能帶性質(zhì)是反轉(zhuǎn)的.則兩個光子晶體共同帶隙處存在確定性界面態(tài),如圖10(a).圖10(b)為波矢kx=0.15(2π/a)、歸一化頻率為0.360c/a的Ez本征場分布.

    接著將圖8光子晶體C和D結(jié)構(gòu)的原胞中構(gòu)成兩個等腰三角形的介質(zhì)柱群以各自三角形中心為原點(diǎn),均順時針旋轉(zhuǎn)角度θ=15?,得到光子晶體E和F,示意圖如圖11(a).計算的能帶結(jié)構(gòu)如圖11(b)和圖11(c)所示.雖然進(jìn)一步降低了晶格對稱性,但是光子晶體的拓?fù)湫再|(zhì)保持不變,依然有拓?fù)鋺B(tài)在兩個能帶反轉(zhuǎn)的光子晶體界面通道處傳輸,如圖11(d).其中五角星標(biāo)記處的波矢kx=0.2(2π/a),頻率為0.361c/a,此點(diǎn)的本征場分布如圖11(e).

    圖9 (a)光子晶體D能帶結(jié)構(gòu);(b)K1和K2點(diǎn)本征場分布Fig.9.(a)Band diagram of PC with structure D;(b)the eigen field distribution at K1and K2points.

    圖10 (a)光子晶體C和D的投影能帶結(jié)構(gòu);(b)界面態(tài)本征場分布Fig.10.(a)The projected band structure;(b)the Ezeigenif eld distribution of interface state.

    圖11 (a)具有C2對稱性的兩個光子晶體示意圖;(b),(c)光子晶體能帶結(jié)構(gòu);(d)投影能帶;(e)界面態(tài)Ez本征場分布Fig.11.(a)The schematics of two PCs with C2symmetry;(b),(c)band diagrams of PCs;(d)the projected band structure;(e)the Ezeigen field distribution of interface state.

    6 C3v和C2對稱性光子晶體存在確定性界面態(tài)

    值得強(qiáng)調(diào)的是,我們利用復(fù)雜晶胞光子晶體產(chǎn)生的拓?fù)浣缑鎽B(tài),不依賴于破缺的對稱性方式,只要構(gòu)建兩個拓?fù)湫再|(zhì)不同即體能帶反轉(zhuǎn)的光子晶體,那么在界面通道就一定會存在確定性界面態(tài).在具有C3v,C3,C2v和C2對稱性的這四組光子晶體中,任意將體能帶反轉(zhuǎn)的兩個光子晶體組建在一起,就會有界面態(tài)存在.例如C3v對稱的光子晶體A分別與C2v對稱的結(jié)構(gòu)D以及C2對稱的結(jié)構(gòu)F擺放在一起,在共同帶隙處都可以實(shí)現(xiàn)拓?fù)浣缑鎽B(tài).這里以一組光子晶體為例來證明界面態(tài)的存在.如圖12所示,將光子晶體A和F擺放在一起,這兩者的體能帶是反轉(zhuǎn)的,通過計算得到的投影能帶可知存在確定性界面態(tài).這表明光子晶體界面態(tài)不依靠具體結(jié)構(gòu)而取決于體能帶的性質(zhì),這在信息輸運(yùn)方面具有重要的研究意義和價值,因而具有更加廣闊的應(yīng)用空間.

    圖12 (a)光子晶體A和F的原胞結(jié)構(gòu)示意圖;(b)投影能帶Fig.12.(a)The schematics of unit cell for PC A and F;(b)the projected band structure.

    7 討 論

    用光子晶體構(gòu)建了類石墨烯復(fù)雜晶胞結(jié)構(gòu).通過多種方式打破蜂窩狀晶格光子晶體的C6v對稱性,獲得了C3v,C3,C2v和C2對稱性的光子晶體,從而可以破缺原本簡并的狄拉克點(diǎn)得到全帶隙的能谷.基于存在確定性界面態(tài)的物理機(jī)理,就是不依賴于破缺光子晶體對稱性的方式,只取決于界面兩邊兩個光子晶體拓?fù)湫再|(zhì)不同即體能帶性質(zhì)反轉(zhuǎn).我們用兩個復(fù)雜晶胞光子晶體構(gòu)建界面通道并實(shí)現(xiàn)了拓?fù)浣缑鎽B(tài).利用拓?fù)鋺B(tài)的缺陷免疫和單向性質(zhì)可以靈活調(diào)控光在晶體內(nèi)部和界面的傳播.相較于石墨烯實(shí)現(xiàn)拓?fù)浣^緣體,利用光子晶體的宏觀特性可以更加直接地在實(shí)驗(yàn)中激發(fā)、探測光拓?fù)鋺B(tài),從而提供嶄新的光調(diào)控方式.基于界面態(tài)在界面的拓?fù)漭斶\(yùn),還可以為設(shè)計新型信息傳播器件提供良好的平臺.另外,光子晶體能谷還有更多可調(diào)控的自由度需要去研究與探討.若能夠直接而有效地選擇性激發(fā)不等價的光子能谷,進(jìn)而方便地探測和觀察到能谷偏振現(xiàn)象,這將對拓?fù)漭斶\(yùn)信息具有重大的意義.

    感謝蘇州大學(xué)物理與光電.能源學(xué)部蔣建華教授及武漢大學(xué)物理學(xué)院邱春印教授對本文的有益討論.

    [1]Klitzing K V,Dorda G,Pepper M 1980Phys.Rev.Lett.45 494

    [2]Thouless D J,Kohmoto M,Nightingale M P,den Nijs M 1982Phys.Rev.Lett.49 405

    [3]Kane C L,Mele E J 2005Phys.Rev.Lett.95 146802

    [4]Novoselov K S,Geim A K,Morozov S V,Jiang D,Zhang Y,Dubonos S V,Grigorieva I V,Firsov A A 2004Science306 666

    [5]Yablonovitch E 1987Phys.Rev.Lett.58 2059

    [6]John S 1987Phys.Rev.Lett.58 2486

    [7]Sakoda K 2004Optical Properties of Photonic Crystals(2nd Ed.)(Berlin:Springer)

    [8]Joannopoulos J D,Johnson S G,Winn J N,Meade R D 2008Photonic Crystals:Molding the Flow of Light(2nd Ed.)(New Jersey:Princeton University Press)

    [9]Mekis A,Chen J C,Kurland I,Fan S,Villeneuve P R,Joannopoulos J D 1996Phys.Rev.Lett.77 3787

    [10]Lin S Y,Chow E,Hietala V,Villeneuve P R,Joannopoulos J D 1998Science282 274

    [11]Robertson W M,Arjavalingam G,Meade R D,Brommer K D,Rappe A M,Joannopoulos J D 1993Opt.Lett.18 528

    [12]Istrate E,Sargent E H 2006Rev.Mod.Phys.78 455

    [13]Guo J,Sun Y,Zhang Y,Li H,Jiang H,Chen H 2008Phys.Rev.E78 026607

    [14]Meade R D,Brommer K D,Rappe A M,Joannopoulos J D 1991Phys.Rev.B44 10961

    [15]Ramos-Mendieta F,Halevi P 1999Phys.Rev.B59 15112

    [16]Choi H G,Oh S S,Lee S G,Kim M W,Kim J E,Park H Y,Kee C S 2006J.Appl.Phys.100 123105

    [17]Xiao M,Zhang Z Q,Chan C T 2014Phys.Rev.X4 021017

    [18]Huang X Q,Xiao M,Zhang Z Q,Chan C T 2014Phys.Rev.B90 075423

    [19]Yang Y T,Huang X Q,Hang Z H 2016Phys.Rev.Appl.5 034009

    [20]Huang X Q,Yang Y T,Hang Z H,Zhang Z Q,Chan C T 2016Phys.Rev.B93 085415

    [21]Yang Y T,Xu T,Xu X F,Hang Z H 2017Opt.Lett.42 3085

    [22]Rycerz A,Jakub T J,Beenakker C W J 2007Nature Phys.3 172

    [23]Xu X D,Yao W,Xiao D,Heinz T F 2014Nature Phys.10 343

    [24]Garcia-Pomar J L,Cortijo A,Nieto-Vesperinas M 2008Phys.Rev.Lett.100 236801

    [25]Xiao D,Yao W,Niu Q 2007Phys.Rev.Lett.99 236809

    [26]Mak K F,McGill K L,Park J,McEuen P L 2014Science344 1489

    [27]Enyashin A N,Ivanovskii A L 2011Phys.Status Solidi248 1879

    [28]Huang X Q,Lai Y,Hang Z H,Zheng H H,Chan C T 2011Nature Mater.10 582

    [29]Yang Y T,Xu Y F,Xu T,Wang H X,Jiang J H,Hu X,Hang Z H 2016 arXiv:1610.07780v1

    PACS:78.67.Pt,73.20.At,42.70.Qs,61.48.GhDOI:10.7498/aps.66.227802

    *Project supported by the National Natural Science Foundation of China(Grant No.11574226),the Natural Science Foundation of Jiangsu Province,China(Grant No.BK20170058),and the Priority Academic Program Development(PAPD)of Jiangsu Higher Education Institutions,China.

    ?Corresponding author.E-mail:zhhang@suda.edu.cn

    Deterministic interface states in photonic crystal with graphene-allotrope-like complex unit cells?

    Jia Zi-Yuan Yang Yu-Ting Ji Li-Yu Hang Zhi-Hong?

    (College of Physics,Optoelectronics and Energy and Collaborative Innovation Center of Suzhou Nano Science and Technology,Soochow University,Suzhou 215006,China)

    24 July 2017;revised manuscript

    13 August 2017)

    Topological insulators have aroused much research interest in condensed matter physics in recent years.Topological protected edge states can propagate unidirectionally and backscattering free along the boundaries of the topological insulators’which will be important for future electronic devices for its immunity to defects.Topology is dependent only on the symmetry of lattice of the system rather than its speci fic wave form.Thus,based on the analogy between electronics and photons,photonic topological insulator has also been demonstrated both theoretically and experimentally.Graphene,composed of a monolayer of carbon atoms in honeycomb lattice,exhibits unusual properties due to its intriguing band diagram.Many types of graphene allotropes have been proposed theoretically.However,due to fabrication difficulties,most of graphene allotropes are unavailable.Here,we propose to study two dimensional(2D)photonic crystal(PC)with complex lattices,similar to that of graphene allotrope.The complex PC structure provides more degrees of freedom in manipulating its symmetry.

    Interface states can also exist in the interface region between two PCs,if they have di ff erent topological properties.Without any surface decoration,deterministic interface states can be created when bulk photonic band inversion can be induced and are demonstrated theoretically and experimentally in 2D PCs with square lattice.By controlling the parameters of PCs,their bulk photonic band properties are engineered and topological phase transition occurs.By inverting the bulk photonic band properties,interface states exist in the common band gaps for two PC systems in the gapped region.Similarly,we proceed to complex honeycomb lattice of PCs.By lowering its originalC6vsymmetry toC3v,C3,C2vand evenC2symmetry,the degeneracies of valley Dirac dispersion at the corners of Brillouin zone are lifted.Photonic band inversion occurs in all four symmetries and the deterministic interface states are numerically realized in the interface region between two PCs.Unidirectional propagation of interface state immune to backscattering along the interface channels is demonstrated if a source with proper optical vortex index is utilized.Due to its easy fabrication,PC is a perfect platform to explore the topological properties of complex lattice and these acquired topological optical states can be of bene fit to the control the propagation of light in the photonic waveguide.

    photonic crystal,topology,valley,interface state

    10.7498/aps.66.227802

    ?國家自然科學(xué)基金(批準(zhǔn)號:11574226)、江蘇省自然科學(xué)基金(批準(zhǔn)號:BK20170058)和江蘇省高校優(yōu)勢學(xué)科建設(shè)工程資助的課題.

    ?通信作者.E-mail:zhhang@suda.edu.cn

    猜你喜歡
    性質(zhì)界面結(jié)構(gòu)
    《形而上學(xué)》△卷的結(jié)構(gòu)和位置
    隨機(jī)變量的分布列性質(zhì)的應(yīng)用
    完全平方數(shù)的性質(zhì)及其應(yīng)用
    國企黨委前置研究的“四個界面”
    論結(jié)構(gòu)
    中華詩詞(2019年7期)2019-11-25 01:43:04
    九點(diǎn)圓的性質(zhì)和應(yīng)用
    厲害了,我的性質(zhì)
    基于FANUC PICTURE的虛擬軸坐標(biāo)顯示界面開發(fā)方法研究
    人機(jī)交互界面發(fā)展趨勢研究
    論《日出》的結(jié)構(gòu)
    国产伦精品一区二区三区四那| 免费高清视频大片| 欧美成人性av电影在线观看| 成人av在线播放网站| 亚洲人成电影免费在线| www.熟女人妻精品国产| 免费无遮挡裸体视频| 久久国产乱子伦精品免费另类| 久久精品影院6| 桃红色精品国产亚洲av| 欧美高清性xxxxhd video| 啦啦啦韩国在线观看视频| 亚洲无线观看免费| 国产一级毛片七仙女欲春2| 91字幕亚洲| eeuss影院久久| 长腿黑丝高跟| 久久久久久大精品| 51午夜福利影视在线观看| 欧美日韩乱码在线| 亚洲黑人精品在线| 一区二区三区免费毛片| 高清毛片免费观看视频网站| 国产亚洲av嫩草精品影院| 亚洲av免费高清在线观看| 18禁黄网站禁片午夜丰满| aaaaa片日本免费| 亚洲欧美激情综合另类| 国产v大片淫在线免费观看| 直男gayav资源| 又黄又爽又刺激的免费视频.| 波多野结衣巨乳人妻| 免费在线观看影片大全网站| 丰满乱子伦码专区| 狂野欧美白嫩少妇大欣赏| 啦啦啦观看免费观看视频高清| 99视频精品全部免费 在线| 亚洲成人精品中文字幕电影| 亚州av有码| 嫩草影院新地址| 禁无遮挡网站| 久久久国产成人精品二区| 国产精品一区二区免费欧美| 9191精品国产免费久久| 黄色日韩在线| 午夜免费成人在线视频| 午夜影院日韩av| 国产高清激情床上av| 天堂动漫精品| 免费大片18禁| 俺也久久电影网| 国产成人aa在线观看| aaaaa片日本免费| 亚洲精品乱码久久久v下载方式| 国产精品一区二区免费欧美| 日韩欧美在线乱码| 伦理电影大哥的女人| 国内精品美女久久久久久| 人人妻人人澡欧美一区二区| 国产私拍福利视频在线观看| 性插视频无遮挡在线免费观看| 精品久久久久久久末码| 久99久视频精品免费| 99热这里只有是精品在线观看 | 国产在线男女| 亚洲 国产 在线| 亚洲内射少妇av| 真人一进一出gif抽搐免费| 亚洲人成网站在线播| 日韩大尺度精品在线看网址| 国产成+人综合+亚洲专区| 成年女人看的毛片在线观看| 国产午夜福利久久久久久| 神马国产精品三级电影在线观看| 国产主播在线观看一区二区| 色哟哟·www| 免费人成视频x8x8入口观看| 1024手机看黄色片| 成年女人毛片免费观看观看9| 国产探花在线观看一区二区| 12—13女人毛片做爰片一| avwww免费| 88av欧美| 日本a在线网址| 亚洲人与动物交配视频| 欧美高清成人免费视频www| 婷婷精品国产亚洲av在线| 露出奶头的视频| xxxwww97欧美| 一个人看的www免费观看视频| 黄色女人牲交| 国产真实乱freesex| 精品久久久久久久人妻蜜臀av| 亚洲精品亚洲一区二区| 免费看美女性在线毛片视频| 久久精品国产亚洲av天美| 精品不卡国产一区二区三区| 男女床上黄色一级片免费看| x7x7x7水蜜桃| 看十八女毛片水多多多| 我的女老师完整版在线观看| 国产成人av教育| 在线播放无遮挡| 精品人妻视频免费看| 日本a在线网址| 免费黄网站久久成人精品 | 最近最新免费中文字幕在线| 九色国产91popny在线| 免费在线观看影片大全网站| 精品免费久久久久久久清纯| 天堂影院成人在线观看| 三级国产精品欧美在线观看| 免费大片18禁| 最近最新中文字幕大全电影3| 永久网站在线| 性色av乱码一区二区三区2| 一卡2卡三卡四卡精品乱码亚洲| 国产精品98久久久久久宅男小说| 亚洲国产欧洲综合997久久,| 午夜激情欧美在线| 特级一级黄色大片| 精品福利观看| 在线a可以看的网站| 国产真实伦视频高清在线观看 | 男女之事视频高清在线观看| 欧洲精品卡2卡3卡4卡5卡区| 两个人的视频大全免费| 亚洲性夜色夜夜综合| 国产精品伦人一区二区| 久久国产精品影院| 99在线人妻在线中文字幕| 美女 人体艺术 gogo| 国产成+人综合+亚洲专区| 国产主播在线观看一区二区| 国产人妻一区二区三区在| 欧美日本视频| 少妇丰满av| 十八禁人妻一区二区| 窝窝影院91人妻| 99国产极品粉嫩在线观看| 国产一区二区在线av高清观看| 亚洲av第一区精品v没综合| 成人国产综合亚洲| 国产亚洲精品久久久久久毛片| 亚洲一区二区三区色噜噜| 婷婷色综合大香蕉| a级毛片免费高清观看在线播放| 精品欧美国产一区二区三| 色尼玛亚洲综合影院| 天天一区二区日本电影三级| 国产精品影院久久| 久久人人爽人人爽人人片va | 麻豆av噜噜一区二区三区| 国产伦精品一区二区三区四那| 午夜两性在线视频| 在线观看一区二区三区| 欧美日韩乱码在线| 在线看三级毛片| 成人欧美大片| 色哟哟哟哟哟哟| 美女高潮喷水抽搐中文字幕| 中出人妻视频一区二区| 精品国产亚洲在线| 免费大片18禁| 欧美最黄视频在线播放免费| 欧美日本亚洲视频在线播放| 久久午夜亚洲精品久久| 国产成+人综合+亚洲专区| 他把我摸到了高潮在线观看| 亚洲成a人片在线一区二区| 网址你懂的国产日韩在线| 一二三四社区在线视频社区8| 久久久久性生活片| 99视频精品全部免费 在线| a级毛片免费高清观看在线播放| 特大巨黑吊av在线直播| 国产极品精品免费视频能看的| 国产精品自产拍在线观看55亚洲| 直男gayav资源| 国产一级毛片七仙女欲春2| 97超视频在线观看视频| 欧美成人a在线观看| 日本黄色片子视频| 亚洲中文日韩欧美视频| 一区二区三区高清视频在线| 在线免费观看的www视频| 亚洲国产精品成人综合色| 久久99热6这里只有精品| 久久久久免费精品人妻一区二区| 夜夜看夜夜爽夜夜摸| 五月伊人婷婷丁香| 国产精品精品国产色婷婷| 日日摸夜夜添夜夜添av毛片 | 亚洲专区国产一区二区| 村上凉子中文字幕在线| 亚洲av电影在线进入| 一级作爱视频免费观看| 国产在线精品亚洲第一网站| 亚洲专区国产一区二区| 在线观看一区二区三区| а√天堂www在线а√下载| 女人被狂操c到高潮| 一本综合久久免费| 亚洲成人精品中文字幕电影| 美女xxoo啪啪120秒动态图 | 欧美又色又爽又黄视频| or卡值多少钱| 亚洲欧美激情综合另类| 9191精品国产免费久久| 男人舔奶头视频| 舔av片在线| 欧美色欧美亚洲另类二区| 亚洲午夜理论影院| 国产av在哪里看| 亚洲av一区综合| 欧美日韩福利视频一区二区| 老熟妇乱子伦视频在线观看| 国产精品自产拍在线观看55亚洲| 怎么达到女性高潮| 久久亚洲精品不卡| 国产日本99.免费观看| 国产视频一区二区在线看| 我要看日韩黄色一级片| 中文字幕人妻熟人妻熟丝袜美| 神马国产精品三级电影在线观看| 日韩亚洲欧美综合| 成人国产综合亚洲| 亚洲最大成人av| 精品久久久久久久久久久久久| 国产精品久久久久久亚洲av鲁大| 在现免费观看毛片| 激情在线观看视频在线高清| 偷拍熟女少妇极品色| 亚洲人成网站在线播| 在线a可以看的网站| 草草在线视频免费看| 给我免费播放毛片高清在线观看| 少妇高潮的动态图| 乱人视频在线观看| 亚洲av成人精品一区久久| 成人三级黄色视频| 国产成人av教育| 一级黄色大片毛片| 一夜夜www| 可以在线观看毛片的网站| 国产高清激情床上av| 国内精品久久久久久久电影| 亚洲国产高清在线一区二区三| 少妇裸体淫交视频免费看高清| 亚洲人成网站高清观看| 久久这里只有精品中国| 最新在线观看一区二区三区| 国产乱人视频| 一本一本综合久久| 在线观看66精品国产| 嫩草影院精品99| 日本一本二区三区精品| 最新中文字幕久久久久| 97超级碰碰碰精品色视频在线观看| 夜夜看夜夜爽夜夜摸| 亚洲av熟女| 中文字幕高清在线视频| 久久中文看片网| 婷婷丁香在线五月| 午夜免费成人在线视频| 日韩亚洲欧美综合| 国产视频一区二区在线看| 一夜夜www| 午夜视频国产福利| 我的女老师完整版在线观看| 国产精品永久免费网站| 村上凉子中文字幕在线| 午夜激情欧美在线| 精品一区二区三区人妻视频| 国产精华一区二区三区| 亚洲专区国产一区二区| 成人午夜高清在线视频| 成人永久免费在线观看视频| 亚洲最大成人av| 国产色婷婷99| 五月玫瑰六月丁香| 精品人妻一区二区三区麻豆 | 久久精品国产亚洲av涩爱 | 激情在线观看视频在线高清| а√天堂www在线а√下载| 看黄色毛片网站| 97超视频在线观看视频| 精品久久久久久久久久久久久| 18禁黄网站禁片免费观看直播| bbb黄色大片| 一区二区三区激情视频| 欧美成人一区二区免费高清观看| 国产白丝娇喘喷水9色精品| 俄罗斯特黄特色一大片| 赤兔流量卡办理| 亚洲18禁久久av| www.熟女人妻精品国产| 国产主播在线观看一区二区| 亚洲久久久久久中文字幕| 色在线成人网| 欧美色欧美亚洲另类二区| 一进一出抽搐gif免费好疼| 亚洲精品久久国产高清桃花| .国产精品久久| 日本a在线网址| netflix在线观看网站| 哪里可以看免费的av片| 国产精品一区二区免费欧美| 久久久久久久久久黄片| 国产在视频线在精品| 午夜福利欧美成人| 蜜桃亚洲精品一区二区三区| 精品人妻视频免费看| 很黄的视频免费| 国产高清视频在线播放一区| 国产人妻一区二区三区在| 亚洲第一区二区三区不卡| 中国美女看黄片| 老熟妇乱子伦视频在线观看| 毛片一级片免费看久久久久 | 亚洲第一区二区三区不卡| 欧美黑人巨大hd| 亚洲国产精品合色在线| 免费观看人在逋| 午夜日韩欧美国产| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av中文字字幕乱码综合| 女生性感内裤真人,穿戴方法视频| 成人美女网站在线观看视频| 免费高清视频大片| 看免费av毛片| 亚洲国产欧美人成| 丁香六月欧美| 综合色av麻豆| 亚洲成av人片在线播放无| 成年女人永久免费观看视频| 日韩大尺度精品在线看网址| 国产黄a三级三级三级人| 久久国产精品人妻蜜桃| 日韩欧美三级三区| 淫妇啪啪啪对白视频| 亚洲av.av天堂| www.熟女人妻精品国产| 亚洲av.av天堂| 直男gayav资源| 国产主播在线观看一区二区| 尤物成人国产欧美一区二区三区| 美女黄网站色视频| 久久久久久久亚洲中文字幕 | 午夜免费男女啪啪视频观看 | 国产精品免费一区二区三区在线| 十八禁国产超污无遮挡网站| 久9热在线精品视频| 亚洲人成网站在线播放欧美日韩| 日韩高清综合在线| 欧美三级亚洲精品| 亚洲av免费高清在线观看| 亚洲第一区二区三区不卡| 亚洲欧美清纯卡通| 身体一侧抽搐| 亚洲人成电影免费在线| 男人舔奶头视频| 国产中年淑女户外野战色| 日韩免费av在线播放| 国产人妻一区二区三区在| 99久久久亚洲精品蜜臀av| 亚洲欧美日韩高清专用| 悠悠久久av| 一级作爱视频免费观看| 嫩草影视91久久| 少妇人妻精品综合一区二区 | 好男人在线观看高清免费视频| 成人特级黄色片久久久久久久| 久久久久九九精品影院| 成年女人看的毛片在线观看| 一级作爱视频免费观看| 此物有八面人人有两片| 夜夜夜夜夜久久久久| 91午夜精品亚洲一区二区三区 | 日韩欧美精品v在线| 精品乱码久久久久久99久播| 欧美色欧美亚洲另类二区| 亚洲专区国产一区二区| 一进一出好大好爽视频| 婷婷色综合大香蕉| 亚洲成a人片在线一区二区| 中文字幕av成人在线电影| 中国美女看黄片| 日韩精品中文字幕看吧| 乱人视频在线观看| 亚洲中文日韩欧美视频| 亚洲精品色激情综合| 老司机福利观看| 日本与韩国留学比较| 久久久久九九精品影院| 好看av亚洲va欧美ⅴa在| 少妇的逼好多水| 日本三级黄在线观看| 午夜影院日韩av| 亚洲精品久久国产高清桃花| 天天一区二区日本电影三级| 噜噜噜噜噜久久久久久91| 看免费av毛片| 成人三级黄色视频| 国产亚洲av嫩草精品影院| 伦理电影大哥的女人| 丰满乱子伦码专区| 乱码一卡2卡4卡精品| 亚洲中文字幕一区二区三区有码在线看| 搡老妇女老女人老熟妇| 少妇的逼水好多| 日本精品一区二区三区蜜桃| 老司机福利观看| 国产三级中文精品| 日本一二三区视频观看| 成人国产综合亚洲| 国产精品爽爽va在线观看网站| 日韩欧美国产在线观看| 亚洲成人中文字幕在线播放| 男人和女人高潮做爰伦理| 又黄又爽又免费观看的视频| 欧美日韩中文字幕国产精品一区二区三区| av黄色大香蕉| 深夜a级毛片| 亚洲人成电影免费在线| 女同久久另类99精品国产91| 偷拍熟女少妇极品色| 最近最新中文字幕大全电影3| 精品国产亚洲在线| 69av精品久久久久久| 人妻制服诱惑在线中文字幕| 国产精品自产拍在线观看55亚洲| 小蜜桃在线观看免费完整版高清| 午夜福利免费观看在线| 日韩国内少妇激情av| 欧美成人a在线观看| 亚洲国产高清在线一区二区三| 一级作爱视频免费观看| 久久国产乱子伦精品免费另类| 欧美另类亚洲清纯唯美| 亚洲自偷自拍三级| 9191精品国产免费久久| 丝袜美腿在线中文| 午夜福利在线观看免费完整高清在 | 午夜福利视频1000在线观看| 久久热精品热| 91麻豆av在线| 日韩欧美三级三区| 日韩欧美在线乱码| 免费人成视频x8x8入口观看| 不卡一级毛片| 夜夜夜夜夜久久久久| 久久精品久久久久久噜噜老黄 | 全区人妻精品视频| 少妇人妻精品综合一区二区 | 免费在线观看成人毛片| 特级一级黄色大片| 午夜亚洲福利在线播放| 在线a可以看的网站| 精品一区二区三区视频在线| 久久国产乱子伦精品免费另类| 亚洲美女搞黄在线观看 | 久久精品久久久久久噜噜老黄 | 国内精品美女久久久久久| 一级黄片播放器| 婷婷色综合大香蕉| 99热这里只有是精品50| or卡值多少钱| 亚洲av二区三区四区| 在现免费观看毛片| 亚洲美女黄片视频| 日韩欧美国产在线观看| 成人特级av手机在线观看| 内射极品少妇av片p| 精品久久久久久久末码| 久久热精品热| 校园春色视频在线观看| 一二三四社区在线视频社区8| 在线观看舔阴道视频| 免费搜索国产男女视频| 1000部很黄的大片| 亚洲精品456在线播放app | 十八禁人妻一区二区| 中文字幕av成人在线电影| 91麻豆av在线| 天堂动漫精品| 国产精品99久久久久久久久| 一本久久中文字幕| 亚洲真实伦在线观看| 亚洲男人的天堂狠狠| 深夜精品福利| 久久久久久久久久成人| 美女cb高潮喷水在线观看| 亚洲av成人不卡在线观看播放网| 欧美高清性xxxxhd video| 久久精品影院6| 亚洲,欧美,日韩| 性色avwww在线观看| 国产三级中文精品| 色尼玛亚洲综合影院| 国产精品三级大全| 伦理电影大哥的女人| 1024手机看黄色片| 美女免费视频网站| 国内精品久久久久精免费| 欧美性感艳星| 九九在线视频观看精品| aaaaa片日本免费| 久久婷婷人人爽人人干人人爱| 99热精品在线国产| 成人午夜高清在线视频| 成人av一区二区三区在线看| 中出人妻视频一区二区| 亚洲国产欧洲综合997久久,| 免费看光身美女| 欧美成狂野欧美在线观看| 午夜福利18| 最近最新中文字幕大全电影3| or卡值多少钱| 在线观看免费视频日本深夜| 很黄的视频免费| 舔av片在线| 91在线观看av| 久久久久久久午夜电影| 国产伦精品一区二区三区四那| 亚洲无线在线观看| 人人妻人人看人人澡| 亚洲狠狠婷婷综合久久图片| a级一级毛片免费在线观看| 综合色av麻豆| 国产伦一二天堂av在线观看| 最近中文字幕高清免费大全6 | 看片在线看免费视频| а√天堂www在线а√下载| 身体一侧抽搐| 人妻久久中文字幕网| 长腿黑丝高跟| 伦理电影大哥的女人| 国产精品美女特级片免费视频播放器| 精品国产三级普通话版| 中文字幕高清在线视频| 18美女黄网站色大片免费观看| 嫩草影视91久久| www.熟女人妻精品国产| 最近视频中文字幕2019在线8| 变态另类成人亚洲欧美熟女| 亚洲人成伊人成综合网2020| 国产私拍福利视频在线观看| 国产精品亚洲美女久久久| 亚洲最大成人手机在线| 欧美在线黄色| 看片在线看免费视频| 中亚洲国语对白在线视频| 国产精品久久久久久亚洲av鲁大| 亚洲天堂国产精品一区在线| 偷拍熟女少妇极品色| 国产极品精品免费视频能看的| 国产精品久久久久久久久免 | 乱码一卡2卡4卡精品| 欧美一区二区亚洲| 国产一级毛片七仙女欲春2| 亚洲欧美精品综合久久99| 国产黄a三级三级三级人| 精品福利观看| 久久久久性生活片| 88av欧美| 国产精品国产高清国产av| bbb黄色大片| 久久久久久久久大av| 国产伦在线观看视频一区| 好男人在线观看高清免费视频| 久久精品国产自在天天线| 亚洲av中文字字幕乱码综合| 在线观看午夜福利视频| 精品久久久久久,| 免费观看精品视频网站| 性插视频无遮挡在线免费观看| 99国产极品粉嫩在线观看| av在线老鸭窝| 一区二区三区四区激情视频 | 欧美另类亚洲清纯唯美| 九九热线精品视视频播放| 熟女电影av网| 精品午夜福利视频在线观看一区| 床上黄色一级片| 国产黄片美女视频| 老司机午夜十八禁免费视频| 美女高潮喷水抽搐中文字幕| 国产精品国产高清国产av| 一级黄片播放器| 国产成年人精品一区二区| 99久久九九国产精品国产免费| 精品无人区乱码1区二区| 亚洲av美国av| 白带黄色成豆腐渣| aaaaa片日本免费| 国产精品一区二区三区四区久久| 亚洲av五月六月丁香网| 99久久精品国产亚洲精品| 亚洲精品色激情综合| 免费无遮挡裸体视频| 亚洲成a人片在线一区二区| 免费高清视频大片| 变态另类成人亚洲欧美熟女| 男人狂女人下面高潮的视频| 久久久久久国产a免费观看| 亚洲中文字幕日韩| 欧美极品一区二区三区四区| 九九热线精品视视频播放| 午夜福利在线观看免费完整高清在 | 91久久精品电影网| av黄色大香蕉| av视频在线观看入口| 亚洲真实伦在线观看| 欧美最黄视频在线播放免费|