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

    疇壁系統(tǒng)中的非厄米趨膚效應(yīng)

    2022-09-14 10:08:22鄧天舒
    物理學(xué)報(bào) 2022年17期
    關(guān)鍵詞:布洛赫哈密頓量布里淵

    鄧天舒

    (清華大學(xué)高等研究院,北京 100084)

    非厄米趨膚效應(yīng)是近幾年非厄米物理研究領(lǐng)域中的熱點(diǎn)問題,它揭示了非厄米系統(tǒng)中體態(tài)波函數(shù)和能譜計(jì)算會敏感依賴于邊界條件的新奇現(xiàn)象.人們提出廣義布里淵區(qū)的概念用以刻畫非厄米系統(tǒng)中的體態(tài)波函數(shù)和能帶性質(zhì).基于廣義布里淵區(qū)計(jì)算的非布洛赫拓?fù)鋽?shù)可以重新構(gòu)建非厄米拓?fù)潴w邊對應(yīng)關(guān)系.然而,過去關(guān)于非厄米趨膚效應(yīng)的討論主要針對開放邊界條件,如果采用疇壁邊界條件,廣義布里淵區(qū)和非布洛赫拓?fù)鋽?shù)的計(jì)算都需要重新考慮.本文綜述了近幾年關(guān)于疇壁邊界條件下非厄米趨膚效應(yīng)的若干研究工作,首先從一般的一維非厄米單帶模型出發(fā),推導(dǎo)廣義布里淵區(qū)方程的一般形式;然后回顧了非厄米SSH (Su-Schieffer-Heeger)模型中廣義布里淵區(qū)和非布洛赫拓?fù)鋽?shù)的計(jì)算;最后在一維光量子行走的系統(tǒng)中,介紹了實(shí)驗(yàn)上非厄米趨膚效應(yīng)的實(shí)現(xiàn)和非厄米拓?fù)溥吘墤B(tài)的探測.

    1 引言

    在孤立的量子系統(tǒng)中利用厄米型的哈密頓量來研究物相的基本性質(zhì)是量子物理研究的基本范式.然而實(shí)際上,絕對孤立的系統(tǒng)是不存在的,因此人們發(fā)展了一系列理論用于描述與環(huán)境有相互作用的系統(tǒng)的演化,其中一個最為經(jīng)典的描述就是在馬爾科夫近似下得到的Lindblad 量子主方程[1,2].量子主方程的研究常常需要涉及到非厄米的哈密頓量: 一方面,在忽略量子躍遷效應(yīng)時(shí),Lindblad型主方程可以直接近似為非厄米哈密頓量下的薛定諤方程[3,4];另一方面,如果把密度矩陣映射為雙希爾伯特空間的波函數(shù),那么其在主方程下的演化可以嚴(yán)格映射為雙空間波函數(shù)在非厄米哈密頓量下的演化[5—7].相較于量子主方程,非厄米哈密頓量的描述在形式上顯得更為簡潔,將之與傳統(tǒng)厄米哈密頓量中的各種經(jīng)典理論結(jié)合,就有可能建立新的理論框架,做出新的理論預(yù)言.如非厄米拓?fù)淠軒Ю碚摰陌l(fā)展就是將非厄米哈密頓量與傳統(tǒng)拓?fù)淠軒ЫY(jié)合的典型范例[8—12].此外,近年來隨著實(shí)驗(yàn)技術(shù)的發(fā)展,人們有能力在越來越多的物理系統(tǒng)中模擬非厄米哈密頓量[13—19],這進(jìn)一步促使人們探索非厄米系統(tǒng)中的種種新奇現(xiàn)象.其中廣受關(guān)注的一個成果就是關(guān)于非厄米趨膚效應(yīng)的研究[20—25].

    非厄米趨膚效應(yīng)指的是在非厄米系統(tǒng)中,體態(tài)波函數(shù)局域在邊界附近的現(xiàn)象.為了刻畫這種效應(yīng),需要將傳統(tǒng)的布洛赫能帶理論修正為非布洛赫能帶理論.以一維系統(tǒng)為例,在傳統(tǒng)的布洛赫能帶理論中,無論是開邊界條件還是周期邊界條件,體態(tài)波函數(shù)都滿足布洛赫波的形式ψ(x)=eikx.然而在非厄米系統(tǒng)中,周期邊界條件下的體態(tài)波函數(shù)依然是布洛赫波;但是對于存在邊界的系統(tǒng),波函數(shù)有可能會局域在邊界附近,其體態(tài)波函數(shù)會滿足ψ(x)=βx的形式,其中β=reik,r和k都是實(shí)數(shù).這里的β會在復(fù)平面內(nèi)形成一條封閉曲線,被稱為廣義布里淵區(qū).在不同的系統(tǒng)中,不同邊界條件下求解這條封閉曲線的方程被稱為廣義布里淵區(qū)方程.基于廣義布里淵區(qū)可以正確地計(jì)算不同邊界下的實(shí)空間連續(xù)譜,即非布洛赫能帶.非布洛赫能帶理論最初的一個成功應(yīng)用就是在非厄米拓?fù)鋯栴}中.在廣義布里淵區(qū)上做積分得到的非布洛赫拓?fù)鋽?shù),可以正確地預(yù)測系統(tǒng)中拓?fù)溥吘墤B(tài)的產(chǎn)生與消失,而布洛赫拓?fù)鋽?shù)做不到這一點(diǎn)[21].這一案例表明,廣義布里淵區(qū)和非布洛赫能帶理論對理解和描述非厄米系統(tǒng)有著重要的意義.

    非厄米趨膚效應(yīng)被提出以來就吸引了凝聚態(tài)理論各界的廣泛關(guān)注,非布洛赫能帶理論取得了蓬勃發(fā)展,出現(xiàn)了很多重要成果[26—43],如開邊界條件下一維系統(tǒng)中廣義布里淵區(qū)方程的一般形式[26]量子主方程中的非厄米趨膚效應(yīng)[32,34,37,38]、非厄米趨膚效應(yīng)的拓?fù)淦鹪碵29,30]等.以上絕大多數(shù)的討論都是針對開邊界條件,而廣義布里淵區(qū)方程對邊界條件的依賴非常敏感,所以此時(shí)研究多種不同邊界條件下的廣義布里淵區(qū)是十分重要且必要的.本文重點(diǎn)關(guān)注疇壁邊界條件下的非厄米趨膚效應(yīng)[17,44].疇壁邊界條件指的是把兩塊不同參數(shù)的材料拼接在一起形成的一種空間結(jié)構(gòu),在拓?fù)淠軒Ю碚撗芯恐杏兄鴱V泛的應(yīng)用[45—47].一方面,在該系統(tǒng)中研究廣義布里淵區(qū)可以幫助構(gòu)建疇壁邊界條件下的非厄米體邊對應(yīng)關(guān)系;另一方面,在很多量子模擬的實(shí)驗(yàn)系統(tǒng)中,疇壁邊界條件比開邊界條件更容易實(shí)現(xiàn)(如光量子行走系統(tǒng)),因此對疇壁邊界條件下非厄米系統(tǒng)的深入研究也為實(shí)驗(yàn)實(shí)現(xiàn)非厄米趨膚效應(yīng)提供了更多的可能.

    本文分為以下幾個部分對疇壁邊界條件下非厄米趨膚效應(yīng)的相關(guān)研究展開綜述: 第2 節(jié)介紹在一般的一維單帶緊束縛模型中,廣義布里淵區(qū)方程的推導(dǎo);第3 節(jié)以非厄米SSH (Su-Schieffer-Heeger)模型為例,介紹疇壁邊界條件下非布洛赫能帶理論的應(yīng)用,以及非厄米拓?fù)潴w邊對應(yīng)關(guān)系的構(gòu)建;第4 節(jié)探討非厄米趨膚效應(yīng)在一維量子行走體系中的實(shí)現(xiàn)方案,并再次應(yīng)用非布洛赫能帶理構(gòu)建非厄米拓?fù)潴w邊對應(yīng)關(guān)系,最后討論實(shí)驗(yàn)上對拓?fù)溥吘墤B(tài)的探測.

    2 疇壁邊界條件下廣義布里淵區(qū)的計(jì)算

    本節(jié)從一維單帶模型出發(fā),推導(dǎo)疇壁邊界條件下的廣義布里淵區(qū)方程.考慮由兩個區(qū)域組成的一維單帶模型,兩個區(qū)域在邊界拼接到一起,形成疇壁邊界條件.如圖1 所示,其哈密頓量可寫作代表左右兩個區(qū)域內(nèi)部相差n個格點(diǎn)的躍遷,m代表最大躍遷幅度.而tM,n代表所有跨越邊界的躍遷.圖2 給出了該模型中哈密頓量的所有本征波函數(shù)分布圖,可以看到,當(dāng),哈密頓量是非厄米的,波函數(shù)局域在了兩個疇壁附近,即非厄米趨膚效應(yīng).在遠(yuǎn)離邊界的左邊(右邊)區(qū)域內(nèi)部,系統(tǒng)具有平移對稱性,實(shí)空間波函數(shù)可分別表示為這里的βx可視作厄米哈密頓量中布洛赫波形式 eikx的推廣.對于給定的能量E,通過分析實(shí)空間的薛定諤方程可以得到特征方程hα(βα)=E,其中

    圖1 疇壁邊界條件下的一維單帶模型示意圖,格點(diǎn)首尾相連形成環(huán)狀結(jié)構(gòu).躍遷距離最遠(yuǎn)為mFig.1.Illustration of single-band lattice model under domain-wall configuration.The hopping distance is m at most.

    圖2 波函數(shù)在實(shí)空間的分布圖.其中 m =2,tM,1=tM,2=tM,-1=tM,-2=1,t L,-2=4,t L,-1=-1,t L,1=3,tL,2=2,t R,-2=2,t R,-1=1,t R,1=3,t R,2=-2,N=30Fig.2.Wave-function distribution in real space.Here,m=2,t M,1=tM,2=tM,-1=tM,-2=1,t L,-2=4,tL,-1=-1,t L,1=3,t L,2=2,t R,-2=2,t R,-1=1,tR,1=3,t R,2=-2,N =30 .

    其中,Cn(Dn)是通解中的待定系數(shù).需要注意的是此時(shí)得到的通解形式只考慮了哈密頓量左右區(qū)域內(nèi)部的躍遷tL,n和tR,n,而沒有考慮跨越邊界的躍遷tM,n.這是因?yàn)閠M,n作用在了邊界附近的 4m個格點(diǎn),只有當(dāng)討論邊界條件時(shí)才會引入tM,n的影響.

    接下來把通解波函數(shù)代入到疇壁邊界條件下的哈密頓量中,得到 4m個邊界方程:

    這里n=1,2,3,···,N,因此(12)式在熱力學(xué)極限下可求得無窮多個解[21].綜上,|η2m|=|η2m+1|即為廣義布里淵區(qū)方程,據(jù)此可以求出疇壁邊界條件下的廣義布里淵區(qū)和非布洛赫能帶的能量.值得注意的一點(diǎn)是,廣義布里淵區(qū)的計(jì)算與tM,n參數(shù)的具體取值無關(guān).這體現(xiàn)出廣義布里淵區(qū)方程具有抗干擾性,只要兩塊材料在邊界上相連,相連處躍遷參數(shù)的擾動并不會影響兩塊體材內(nèi)部廣義布里淵區(qū)的計(jì)算結(jié)果.

    3 非厄米SSH 模型

    本節(jié)將回顧疇壁邊界條件下的非厄米SSH 模型[44].以此模型為例,具體闡釋廣義布里淵區(qū)的計(jì)算,并將之應(yīng)用到非厄米能帶理論中,計(jì)算非布洛赫能量和非布洛赫拓?fù)鋽?shù),驗(yàn)證疇壁邊界條件下非厄米系統(tǒng)中的體邊對應(yīng)關(guān)系.

    3.1 模 型

    如圖3 所示,考慮疇壁邊界條件下的非厄米SSH 模型.開邊界條件下的非厄米SSH 模型過去有過許多討論[48,49],而在文獻(xiàn)[21]中首次利用廣義布里淵區(qū)概念解釋了該模型中反常的體邊對應(yīng)現(xiàn)象.而疇壁邊界條件下非厄米SSH 模型的哈密頓量可寫作

    圖3 疇壁邊 界條件 下的非厄米SSH 模 型.在左邊(右邊)區(qū)域元胞內(nèi)部a 子格到b 子格的躍遷是-γ/2,而元胞內(nèi)部b 子格到a 子格的是 t L1(R)+γ/2 .此外,疇壁兩邊不同區(qū)域具有相同的 t2 和 γ [44]Fig.3.A non-Hermitian SSH model with two bulks.For left(right) bulk,the intra-cell hopping from a-site to b-site is-γ/2while the intra-cell hopping from b-site to asite is +γ/2 .Also,different bulks holds same t2 and γ[44].

    該模型具有手征對稱性σzh(k)σz=-h(k) 和非厄米趨膚效應(yīng),是一個探討非厄米趨膚效應(yīng)和非厄米拓?fù)淠軒Ю碚摰慕?jīng)典模型.開邊界條件下的非厄米SSH 模型已經(jīng)有過詳細(xì)的討論[21],在本節(jié)后面的部分,將集中討論疇壁邊界條件下非厄米SSH 模型中廣義布里淵區(qū)的計(jì)算方法,并基于此構(gòu)建非厄米拓?fù)潴w邊對應(yīng)關(guān)系.

    3.2 廣義布里淵區(qū)

    本節(jié)采用第2 節(jié)中的方法,具體探討廣義布里淵區(qū)的計(jì)算.首先可以將實(shí)空間哈密頓量中特定能量的本征態(tài)寫作

    (18)式可化作一個關(guān)于β的一元二次方程,方程有兩個根,所以(16)式中的nmax取 2 .

    將(15)式,(19)式和(20)式代入邊界條件:

    由第2 節(jié)中相似的邏輯便可判斷,體態(tài)的波函數(shù)解的數(shù)目一定正比于N,為了滿足這一點(diǎn),集合{ηi|i=1,2,3,4}中的元素一定滿足

    這就是疇壁邊界條件下非厄米SSH 模型的廣義布里淵區(qū)方程.(29)式與文獻(xiàn)[44]中所給的廣義布里淵區(qū)方程完全一致,是一個更為簡潔的形式.

    接下來詳細(xì)闡述如何由廣義布里淵區(qū)方程計(jì)算得到廣義布里淵區(qū)和非布洛赫能帶.形式上,(29)式是關(guān)于4 個變量βL,1,βL,2,βR,1,βR,2的方程.但是從(18)式可以看到,所有的ηi都由能量E決定,所以方程(29)中只有一個獨(dú)立的復(fù)變量(在非厄米系統(tǒng)中E和β一般都是復(fù)數(shù)).因此在方程(29) 的約束下,便可確定任意一個βα,n在復(fù)平面內(nèi)形成的封閉曲線.具體來講,可以從βL,1=reip出發(fā),先將完整的復(fù)變量βL,1=reip形式代入(18)式中,計(jì)算能量E,再通過能量計(jì)算βL,2,βR,1,βR,2.然后把求得的所有{1/βL,1,1/βL,2,βR,1,βR,2}進(jìn)行排序,代入方程(29)中.于是方程(29)將變成一個關(guān)于r和p方程,遍歷p∈[0,2π),便可得到r關(guān)于p的函數(shù),由此便可繪制左邊區(qū)域的廣義布里淵區(qū)的其中一支解.依照此法可求得所有的βα,n.

    圖4(c)和圖4(d)分別給出了左邊區(qū)域和右邊區(qū)域的廣義布里淵區(qū).在開邊界條件下的非厄米SSH 模型中,廣義布里淵區(qū)的方程寫作|β1|=|β2|,因此β的兩支解會在復(fù)平面內(nèi)重合.僅將邊界條件改為疇壁邊界條件后,就會得到完全不同的廣義布里淵區(qū),這也印證了非厄米趨膚效應(yīng)對邊界條件變化的敏感性.

    圖4 (a),(b)體態(tài)波函數(shù)的理論(黑線)和數(shù)值(綠點(diǎn))結(jié)果,(a)圖的鏈長為 N =20,(b) 圖的鏈長為 N =40 .(c) 復(fù)平面內(nèi)βL,1(紅色點(diǎn)劃線)和 β L,2 (藍(lán)色虛線)表示的左鏈的非布洛赫布里淵區(qū).(d) 復(fù)平面內(nèi) β R,1 (紅色點(diǎn)劃線)和 β R,2 (藍(lán)色虛線)表示的右鏈的非布洛赫布里淵區(qū).在圖(c)和圖(d)中,取熱力學(xué)極限 N →∞ 并且用黑色實(shí)線畫出了布洛赫布里淵區(qū)作為對比.在所有的子圖中,選取的參數(shù)為=-t2,=1.5t2,γ =1.33t2 [44]Fig.4.(a),(b) Theoretical (black lines) and numerical (green dots) results of bulk-state energy spectrum.The length of chain is N=20 for panel (a) and N =40 for panel (b).(c) Non-Bloch Brillouin zones of the left bulk,represented by β L,1 (red dash-dotted line) and β L,2 (blue dashed line) on the complex plane.(d) Non-Bloch Brillouin zones of the right bulk,represented byβR,1(red dash-dotted line) and β R,2 (blue dashed line) on the complex plane.In panels (c) and (d),we take the thermodynamic limit N →∞,and we also plot the Bloch Brillouin zones with black solid lines for comparison.For all subplots,we take =-t2,=1.5t2,γ =1.33t2 [44].

    我們又將廣義布里淵區(qū)與布里淵區(qū) eik在復(fù)平面上形成的單位圓做對比.需要注意的是,在兩種情況下,廣義布里淵區(qū)會退化為單位圓.一種情況是當(dāng)哈密頓量為厄米的γ=0,另一種情況是系統(tǒng)具有周期邊界條件=.這表明非厄米趨膚效應(yīng)是系統(tǒng)非厄米性與邊界條件共同作用的結(jié)果.

    此外,本文還基于廣義布里淵區(qū)和特征方程,繪制了非布洛赫能帶(見圖4(a)和圖4(b)),并將之與實(shí)空間嚴(yán)格對角化哈密頓量的能譜作對比.可以看到在N=20 和N=40 的情況下,數(shù)值計(jì)算的能譜和廣義布里淵區(qū)的理論得到的非布洛赫能帶符合得很好,并且隨著N的增大,兩者的差距變得更小.這充分說明本文計(jì)算得到的廣義布里淵區(qū)在疇壁邊界條件下可以精確刻畫熱力學(xué)極限下實(shí)空間的本征能譜.

    3.3 非布洛赫拓?fù)鋽?shù)與拓?fù)潴w邊對應(yīng)關(guān)系

    本節(jié)基于3.2 節(jié)得到的廣義布里淵區(qū),計(jì)算疇壁邊界條件下的非布洛赫拓?fù)洳蛔兞?構(gòu)建非厄米系統(tǒng)中的拓?fù)潴w邊對應(yīng)關(guān)系.由系統(tǒng)的布洛赫哈密頓量(14)可知,系統(tǒng)具有手征對稱性,其拓?fù)洳蛔兞靠梢杂蒢ak 相定義:

    這里νL和νR分別代表左邊區(qū)域和右邊區(qū)域的拓?fù)鋽?shù),兩者之差如果非0,疇壁上就會出現(xiàn)邊緣態(tài),這就是傳統(tǒng)意義上疇壁邊界條件下的體邊對應(yīng)關(guān)系.但是對于非厄米系統(tǒng)來說,需采用非布洛赫能帶理論,將積分區(qū)間從布里淵區(qū)修正為廣義布里淵區(qū).因此將β=r(p)eip和h(β) 代入(30)式便可得到非布洛赫拓?fù)鋽?shù)的計(jì)算公式:

    需要注意的是,從圖4(c)和圖4(d)可以看出,廣義布里淵區(qū)方程一般會求得r(p) 的兩支解,也就是說(32)式中有兩個可以選擇的積分區(qū)間分別對應(yīng)βα,1和βα,2在復(fù)平面上形成的封閉曲線.在非厄米SSH 模型中可以證明,同一個參數(shù)下求得的兩個廣義布里淵區(qū)上積分得到的非布洛赫拓?fù)鋽?shù)是一致的[44].

    圖5(d)—(f)分別把布洛赫和非布洛赫拓?fù)鋽?shù)和數(shù)值計(jì)算的能譜做對比,確認(rèn)了非布洛赫體邊對應(yīng)性.在厄米的情況下γ=0,非布洛赫和布洛赫拓?fù)鋽?shù)的結(jié)果一致,即當(dāng)時(shí)系統(tǒng)存在零模.隨著γ增加,布洛赫和非布洛赫拓?fù)鋽?shù)之間的區(qū)別便顯現(xiàn)出來,左右區(qū)域非布洛赫拓?fù)鋽?shù)的差值可以正確地預(yù)測邊緣零模的出現(xiàn)與消失,而布洛赫拓?fù)鋽?shù)無法做到這一點(diǎn).這進(jìn)一步說明,相比于布洛赫拓?fù)鋽?shù),基于廣義布里淵區(qū)的非布洛赫拓?fù)鋽?shù)更能正確地判斷系統(tǒng)拓?fù)湫再|(zhì)的變化.

    圖5 (a)—(c)顏色欄在 - 平面能譜絕對值的最小值.這里選取 N =40 .圖(a)中 γ =0,圖(b)中 γ =0.67 t2,圖(c)中γ=1.33 t2 .黑色線對應(yīng)參數(shù)=,此時(shí)疇壁構(gòu)型會退化回周期邊界條件的均勻格點(diǎn)模型.玫紅色的虛線對應(yīng)的參數(shù)為=-.圖(a)—(c)里的紅色點(diǎn)劃線給出了用在圖(d)—(f)里的參數(shù),即=1.5 t2 .(d)—(f)能譜絕對值(上部)和拓?fù)鋽?shù)(下部).圖中下部給出了左鏈的布洛赫拓?fù)鋽?shù) νL (玫紅虛線) 和右鏈的非布洛赫拓?fù)鋽?shù) ν R (黑色虛線),以及左鏈的非布洛赫拓?fù)鋽?shù) (紅色實(shí)線)和右鏈的非布洛赫拓?fù)鋽?shù) (藍(lán)色實(shí)線)[44]Fig.5.(a)—(c) Contour plots of absolute values of the energy-spectrum minimum on the -plane.Here,we take N =40 .We also take γ =0 for (a),γ =0.67 t2 for (b) and γ =1.33 t2 for (c).The black solid lines are given by =,where the domainwall configuration is reduced to single homogeneous bulk with a periodic boundary condition.The magenta dashed lines are given by =-.The red dashed-dotted line in panels (a)—(c) correspond to parameters we use in panels (d)—(f) with =1.5 t2 .(d)—(f) The absolute values of the energy spectrum (upper panels) and various winding numbers (lower panels).In the lower panel,we show the Bloch winding numbers for the left bulk νL (magenta dashed lines) and the right bulk ν R (black dashed lines),as well as non-Bloch winding numbers for the left bulk (red solid lines) and the right bulk (blue solid lines)[44].

    此外,我們還發(fā)現(xiàn),即使疇壁某一側(cè)的參數(shù)不變化,其非布洛赫拓?fù)鋽?shù)的計(jì)算卻依賴于廣義布里淵區(qū),進(jìn)而會受到另一側(cè)參數(shù)的影響.例如,圖5(f)中右側(cè)區(qū)域的非布洛赫拓?fù)鋽?shù)會在0.16 附近從 0 變化為 0.5,而此時(shí)右鏈的參數(shù)始終保持不變.這是疇壁邊界條件下的非厄米系統(tǒng)又一個獨(dú)特的性質(zhì).

    4 一維量子行走體系中的非厄米趨膚效應(yīng)

    在第3 節(jié)討論的非厄米SSH 模型中,已經(jīng)從理論上推導(dǎo)了廣義布里淵區(qū)和非布洛赫拓?fù)鋽?shù),并基于此重新構(gòu)建了非厄米系統(tǒng)中的拓?fù)潴w邊對應(yīng)關(guān)系.但是如何在實(shí)驗(yàn)系統(tǒng)中模擬非厄米趨膚效應(yīng)并驗(yàn)證非厄米體邊對應(yīng)關(guān)系,依然是一個重要的問題.本節(jié)將介紹如何在一維光量子行走的實(shí)驗(yàn)平臺中模擬具有非厄米趨膚效應(yīng)的有效哈密頓量,并在該模型中討論廣義布里淵區(qū)和非布洛赫拓?fù)鋽?shù)的計(jì)算,再次驗(yàn)證非厄米體邊對應(yīng)關(guān)系[17].

    4.1 模 型

    量子行走是一類高度可控的周期驅(qū)動系統(tǒng),被廣泛應(yīng)用于量子算法[50,51]和量子模擬[52—55]的研究.近些年來,人們開始致力于在離散量子行走平臺中模擬和研究非厄米拓?fù)湎郲56—58],這為研究具有非厄米趨膚效應(yīng)的拓?fù)浞瞧接鼓P吞峁┝艘粋€良好的基礎(chǔ).

    在量子行走系統(tǒng)中,選定如下的非幺正時(shí)間演化算符:

    通過對角化實(shí)空間的時(shí)間演化算符,得到了所有的本征態(tài)波函數(shù),如圖6(b)所示.可以看到,隨著γ的增加,體態(tài)波函數(shù)會漸漸局域在疇壁邊界附近.由此可以判斷系統(tǒng)具有非厄米趨膚效應(yīng).該模型還具有手征對稱性,σxUσx=U-1,所以也可以在其中探討非厄米拓?fù)湎嘧兒腕w邊對應(yīng)關(guān)系.接下來將在第4.2 節(jié)討論該模型中的廣義布里淵區(qū)的計(jì)算,在第4.3 節(jié)討論非布洛赫拓?fù)鋽?shù)的計(jì)算.

    圖6 具有非厄米趨膚效應(yīng)的量子行走 (a) 疇壁邊界條件示意圖;(b) 體態(tài)波函數(shù)(黑色)和邊緣態(tài)波函數(shù)(紅色)的空間分布,其中是手征對稱算符的本征態(tài);(c) 左邊(右邊)區(qū)域的廣義布里淵區(qū) β L(R),i,i=1,2 代表 β 有兩支解.標(biāo)準(zhǔn)的廣義布里淵區(qū)是單位圓(黑色實(shí) 線).參數(shù)選取為 N =30 ,=0.1875π,=-0.3333π,=0.2π,=-0.6667π,以及 γ =0.5 [17]Fig.6.Quantum walks with non-Hermitian skin effect: (a) Schematic illustration of the domain-wall configuration.(b) Spatial distribution of the projected norms of bulk (black) and edge (red) states,withare eigenstates of the chiral-symmetry operator.(c) Generalized Brillouin zones for left (right) bulk β L(R),i,i =1,2 indicates two solutions of β .Whereas the standard Brillouin zones are indicated by unit circles (solid black line).Parameters: N =30,=0.1875π,=-0.3333π,=0.2π,=-0.6667π,γ=0.5 [17].

    4.2 廣義布里淵區(qū)的計(jì)算

    為了計(jì)算廣義布里淵區(qū),先給出波函數(shù)本征態(tài)的通解形式:

    于是便可得到形式上與(29)式完全一樣的廣義布里淵區(qū)方程,并可按照3.2 節(jié)中的方法求解廣義布里淵區(qū).圖6(c)給出了左右區(qū)域廣義布里淵區(qū)的例子.

    4.3 非布洛赫拓?fù)鋽?shù)與體邊對應(yīng)關(guān)系

    接下來,基于廣義布里淵區(qū),計(jì)算量子行走模型中的非布洛赫拓?fù)鋽?shù).需要注意的是,量子行走是周期驅(qū)動系統(tǒng),即使是對于傳統(tǒng)的幺正時(shí)間演化算符,其拓?fù)湫再|(zhì)與靜態(tài)哈密頓量的情況也有所區(qū)別.對于靜態(tài)哈密頓量,如系統(tǒng)滿足手征對稱性,那么其能譜一定是關(guān)于E=0 對稱分布,能隙也一定出現(xiàn)在零能處.但是對于周期驅(qū)動的系統(tǒng)來說,其準(zhǔn)能量分布在 [-π,π) 區(qū)間,如果系統(tǒng)有手征對稱性,那么準(zhǔn)能量不僅關(guān)于E=0 對稱,也關(guān)于E=π對稱,所以拓?fù)浣^緣體的能隙會出現(xiàn)在準(zhǔn)能量ε=0 和ε=π 兩個位置.因此當(dāng)系統(tǒng)處于拓?fù)浞瞧接箙^(qū)時(shí),其邊緣態(tài)能量會有 0 和 π 兩個模式,相應(yīng)地,也需計(jì)算 0 模和 π 模兩種拓?fù)鋽?shù)ν0和νπ.關(guān)于Floquet 系統(tǒng)中的拓?fù)淅碚?在過往的文章中有過很多討論[59—63].這里先對布洛赫拓?fù)鋽?shù)ν0和νπ的計(jì)算方法做簡單介紹.首先可以定義另一個時(shí)間框架下的時(shí)間演化算符:

    U′和U相當(dāng)于選取了不同的時(shí)間起始點(diǎn),他們具有相同的準(zhǔn)能量能譜,都滿足手征對稱性,對稱性算符均為σx.然后參照(30)式中的方法,分別計(jì)算不同時(shí)間框架下的布洛赫拓?fù)鋽?shù)ν和ν′.注意(30)式中的χk和ψk分別是時(shí)間演化算符的左右本征矢,定義為

    對于傳統(tǒng)的幺正量子行走系統(tǒng),ν0(π)即可預(yù)測系統(tǒng)0(π)模的拓?fù)溥吘墤B(tài),但是在本文討論的非幺正量子行走系統(tǒng)中,布洛赫拓?fù)鋽?shù)并不能正確地刻畫非厄米系統(tǒng)中的拓?fù)湫再|(zhì),因此,還需根據(jù)廣義布里淵區(qū)的計(jì)算結(jié)果和拓?fù)鋽?shù)的計(jì)算公式(32),計(jì)算不同框架下的非布洛赫拓?fù)鋽?shù)和,進(jìn)而得到0(π)模的非布洛赫拓?fù)鋽?shù):

    圖7 黑色實(shí)線為準(zhǔn)能量能譜,紅色實(shí)線和藍(lán)色實(shí)線分別代表零模拓?fù)洳蛔兞康牟钪岛?π 模拓?fù)洳蛔兞康牟钪?灰色虛線和橙色虛線分別代表左右兩邊零模拓?fù)洳蛔兞康牟钪岛?π 模拓?fù)洳蛔兞康牟钪?藍(lán)色點(diǎn)代表圖8 選取的參數(shù).參數(shù)選 擇: =9π/16,=π/4,=3π/4,以 及γ=0.2746[17]Fig.7.Black solid line represents quasienergy spectrum.Red solid line and blue solid line represents zero-mode and π-mode topological invariant difference,respectively.Grey dashed line and orange dashed line represents zero-mode and π -mode topological invariant difference,respectively.The blue dot indicates the parameter for Fig.8.Parameters:=9π/16,=π/4,=3π/4,and γ =0.2746 [17].

    4.4 實(shí)驗(yàn)探測邊緣態(tài)

    前面幾節(jié)已經(jīng)在理論上驗(yàn)證了非布洛赫拓?fù)鋽?shù)與非厄米系統(tǒng)中邊緣態(tài)的對應(yīng)關(guān)系,但是如何在實(shí)驗(yàn)上探測非厄米系統(tǒng)中的邊緣態(tài)依然是一個具有挑戰(zhàn)性的問題.由于非厄米趨膚效應(yīng)的存在,邊緣態(tài)的波函數(shù)分布和體態(tài)波函數(shù)一樣,都是局域在邊界上,因此很難通過直接測量波函數(shù)在實(shí)空間分布的方法識別出拓?fù)溥吘墤B(tài).為了克服這一困難,我們通過對每一步時(shí)間演化的波函數(shù)加權(quán)求和的方法直接將邊緣態(tài)模式提取出來.具體方法如下.首先定義

    然后,把含時(shí)演化的波函數(shù)在本征態(tài)上做展開,得到

    如果準(zhǔn)能量能譜為純實(shí)數(shù)時(shí),那么對于Enε的本征態(tài),其加權(quán)系數(shù)也就是說只要時(shí)間演化足夠長,加權(quán)求和的波函數(shù)|Φ0,π(t)〉中只會留下能量為En=0,π 的本征態(tài),也就是0 模和 π 模的拓?fù)溥吘墤B(tài).需要注意的是,在非幺正量子行走系統(tǒng)中,準(zhǔn)能量能譜不一定總是純實(shí)數(shù),所以需要選取準(zhǔn)能量為純實(shí)數(shù)的參數(shù)區(qū)間進(jìn)行探測.

    圖8 給出了將該方法應(yīng)用到實(shí)驗(yàn)上得到的觀測結(jié)果.在拓?fù)洳蛔兞繚M足=(0,-1) 的參數(shù)點(diǎn),對于不同的初態(tài)選取,實(shí)驗(yàn)上都探測到了Φπ,+(x)的模式,并且其他模式占據(jù)很小.圖8(d)—(f)中也將觀測所得到的Φπ,+(x) 與實(shí)空間數(shù)值對角化得到的邊緣態(tài)模式波函數(shù)做了對比.由兩者相符可以看出Φπ,+(x) 的確提取到了ε=π 的本征邊緣態(tài)模式.在實(shí)驗(yàn)上觀測得到的拓?fù)溥吘墤B(tài)與非布洛赫拓?fù)洳蛔兞恳恢?這也再一次印證了基于廣義布里淵區(qū)構(gòu)建的非厄米拓?fù)潴w邊對應(yīng)關(guān)系.

    圖8 (a)—(c) 不同初態(tài)選擇下,7 步演化之后實(shí)驗(yàn)測量得到的 Φ ε,μ(x) ; (d)—(f) 數(shù)值模擬得到 Φ π,+(x) 和實(shí)驗(yàn)結(jié)果的對比.圖中還給出了數(shù)值對角化 N =15 的疇壁系統(tǒng)對峰值進(jìn)行縮放之后的結(jié)果.對角化得到的邊緣態(tài)進(jìn)行縮放的比例通過擬合Φπ,+(x)的中心峰值得到.其他參數(shù)與圖7 相同[17]Fig.8.(a)—(c) Experimentally measured Φ ε,μ(x) after the seventh step with different initial states;(d)—(f) simulated Φ π,+(x) compared with experimental result.We also show the scaled norms of the edge state with ε =π by diagonalizing a domain-wall system with N =15 .Norms of the edge states from diagonalization are scaled to fit the central peak of the numerically-simulated Φ π,+(x) .Other parameters are as same as Fig.7[17].

    5 總結(jié)與展望

    本文主要介紹了疇壁邊界條件下具有非厄米趨膚效應(yīng)的系統(tǒng)中廣義布里淵區(qū)方程的推導(dǎo)與求解,并由此出發(fā)討論了非布洛赫能帶理論在非厄米SSH 模型和非幺正一維量子行走實(shí)驗(yàn)?zāi)P椭械膽?yīng)用.基于廣義布里淵區(qū)計(jì)算的非布洛赫拓?fù)鋽?shù)可以正確預(yù)測疇壁系統(tǒng)中的拓?fù)溥吘墤B(tài),這印證了非布洛赫能帶理論在不同邊界條件下非厄米系統(tǒng)中的有效性.在不同模型中,疇壁邊界條件下的廣義布里淵區(qū)方程有一個統(tǒng)一的方程形式,對該方程背后物理內(nèi)涵的充分理解,或許可以幫助我們找到一個對多種邊界條件都普適的廣義布里淵區(qū)方程.此外,不同邊界條件下的高維系統(tǒng)中,廣義布里淵區(qū)的計(jì)算分析也尚未被充分討論,有待于日后進(jìn)一步探索.

    猜你喜歡
    布洛赫哈密頓量布里淵
    哈密頓量宇稱-時(shí)間對稱性的刻畫*
    幾種哈密頓量的寫法與變換
    基于布里淵散射的光纖溫度和應(yīng)變快速感知
    光通信研究(2022年1期)2022-02-18 11:58:42
    基于金剛石中不同軸向NV色心的磁力計(jì)的探討
    科技視界(2019年19期)2019-08-29 02:58:06
    意林·作文素材(2019年8期)2019-05-12 15:12:34
    能量均分定理的一種證明
    平行四邊形晶格結(jié)構(gòu)中布里淵區(qū)的研究
    布洛赫“黑暗瞬間”的深意
    人民論壇(2016年22期)2016-12-13 11:04:09
    面心立方、體心立方晶格第一布里淵區(qū)的畫法
    論馬克?布洛赫的比較史情懷
    教師·中(2016年6期)2016-08-23 21:21:51
    国产乱人偷精品视频| 大片免费播放器 马上看| 制服人妻中文乱码| 欧美另类一区| 国产精品三级大全| 亚洲色图 男人天堂 中文字幕| 最新的欧美精品一区二区| 亚洲精品aⅴ在线观看| 色精品久久人妻99蜜桃| 两性夫妻黄色片| 精品久久久久久电影网| 一级毛片 在线播放| xxx大片免费视频| 制服人妻中文乱码| 少妇猛男粗大的猛烈进出视频| 丝袜美足系列| 亚洲欧美精品综合一区二区三区| 丝袜喷水一区| 亚洲av日韩在线播放| 国产97色在线日韩免费| 人妻人人澡人人爽人人| 国产在线一区二区三区精| 看免费成人av毛片| 欧美日韩福利视频一区二区| 美女福利国产在线| 韩国精品一区二区三区| 色94色欧美一区二区| 国产极品粉嫩免费观看在线| 在线观看免费高清a一片| 国产一级毛片在线| 尾随美女入室| 观看av在线不卡| 精品一区二区三区av网在线观看 | 国产xxxxx性猛交| 亚洲av电影在线进入| 日本一区二区免费在线视频| 咕卡用的链子| 久久久久久久久免费视频了| 免费日韩欧美在线观看| 一区二区日韩欧美中文字幕| 这个男人来自地球电影免费观看 | 啦啦啦 在线观看视频| 两性夫妻黄色片| 国产一区二区 视频在线| 亚洲欧美日韩另类电影网站| 精品人妻熟女毛片av久久网站| 欧美精品一区二区免费开放| 亚洲 欧美一区二区三区| 久久 成人 亚洲| 18禁动态无遮挡网站| 国产欧美亚洲国产| 精品一区二区免费观看| 老司机影院成人| videos熟女内射| 亚洲人成77777在线视频| 亚洲av日韩精品久久久久久密 | 亚洲国产成人一精品久久久| 纵有疾风起免费观看全集完整版| 日韩av不卡免费在线播放| 中文精品一卡2卡3卡4更新| 国产精品熟女久久久久浪| 波多野结衣av一区二区av| 国产精品av久久久久免费| e午夜精品久久久久久久| 制服人妻中文乱码| av卡一久久| 国产精品国产三级国产专区5o| 999久久久国产精品视频| 国产成人系列免费观看| 涩涩av久久男人的天堂| a级片在线免费高清观看视频| av有码第一页| 久久免费观看电影| 久久久亚洲精品成人影院| 久久这里只有精品19| 日韩一卡2卡3卡4卡2021年| 看十八女毛片水多多多| 亚洲天堂av无毛| 国产高清国产精品国产三级| 免费日韩欧美在线观看| 亚洲精品一区蜜桃| 深夜精品福利| 国产成人精品无人区| 久久久久精品久久久久真实原创| 欧美中文综合在线视频| 91精品国产国语对白视频| 黄网站色视频无遮挡免费观看| www.自偷自拍.com| 久久毛片免费看一区二区三区| 久久亚洲国产成人精品v| 亚洲,欧美,日韩| 亚洲欧美精品自产自拍| 国产精品国产三级专区第一集| 精品国产一区二区久久| 中国国产av一级| av不卡在线播放| 如何舔出高潮| 免费少妇av软件| 老汉色av国产亚洲站长工具| 男女免费视频国产| 日韩精品有码人妻一区| 国产精品 欧美亚洲| 亚洲国产最新在线播放| 亚洲欧美一区二区三区国产| 亚洲精品国产av蜜桃| 国产成人一区二区在线| 日本色播在线视频| 无限看片的www在线观看| 国产精品麻豆人妻色哟哟久久| 久久人妻熟女aⅴ| 自拍欧美九色日韩亚洲蝌蚪91| 色视频在线一区二区三区| 老司机亚洲免费影院| 人人妻,人人澡人人爽秒播 | 欧美精品高潮呻吟av久久| 人成视频在线观看免费观看| 日本av手机在线免费观看| 国产爽快片一区二区三区| 亚洲av日韩精品久久久久久密 | 熟女av电影| 国产精品蜜桃在线观看| 久久天堂一区二区三区四区| 日韩,欧美,国产一区二区三区| 亚洲精品,欧美精品| 欧美精品av麻豆av| 欧美乱码精品一区二区三区| 99热全是精品| 国产在线一区二区三区精| 青青草视频在线视频观看| 2021少妇久久久久久久久久久| 妹子高潮喷水视频| 午夜免费观看性视频| av视频免费观看在线观看| 看十八女毛片水多多多| 亚洲激情五月婷婷啪啪| 久久久久精品久久久久真实原创| 精品一区在线观看国产| 五月天丁香电影| 国产野战对白在线观看| 不卡视频在线观看欧美| 日韩av不卡免费在线播放| 十分钟在线观看高清视频www| 国精品久久久久久国模美| 建设人人有责人人尽责人人享有的| 久久精品熟女亚洲av麻豆精品| 99精国产麻豆久久婷婷| 欧美日韩一级在线毛片| 无限看片的www在线观看| 两个人免费观看高清视频| 91aial.com中文字幕在线观看| 日日摸夜夜添夜夜爱| 女人高潮潮喷娇喘18禁视频| 十八禁人妻一区二区| 国产成人免费无遮挡视频| 黄色怎么调成土黄色| 亚洲国产日韩一区二区| 亚洲欧洲日产国产| 自线自在国产av| 国产福利在线免费观看视频| 少妇猛男粗大的猛烈进出视频| 午夜免费男女啪啪视频观看| 性少妇av在线| 久久久精品区二区三区| 国产成人精品久久二区二区91 | 精品国产乱码久久久久久男人| 色94色欧美一区二区| 精品国产超薄肉色丝袜足j| 精品亚洲成a人片在线观看| 一区二区日韩欧美中文字幕| 亚洲欧美一区二区三区黑人| 精品国产一区二区久久| 男男h啪啪无遮挡| 欧美日韩亚洲高清精品| 亚洲久久久国产精品| 又大又爽又粗| 男女国产视频网站| 午夜激情av网站| 国产xxxxx性猛交| 满18在线观看网站| 爱豆传媒免费全集在线观看| 亚洲精品国产av成人精品| 十八禁高潮呻吟视频| 美女脱内裤让男人舔精品视频| 天天操日日干夜夜撸| 啦啦啦在线观看免费高清www| 麻豆乱淫一区二区| 亚洲欧美精品自产自拍| 久久韩国三级中文字幕| 国产精品国产av在线观看| 99热国产这里只有精品6| 亚洲精品日本国产第一区| 在线精品无人区一区二区三| 免费久久久久久久精品成人欧美视频| 亚洲精品一区蜜桃| 亚洲熟女精品中文字幕| 老司机影院成人| 亚洲精品美女久久久久99蜜臀 | 多毛熟女@视频| 悠悠久久av| 男女免费视频国产| 精品一区在线观看国产| 日韩一区二区三区影片| 蜜桃在线观看..| 国产精品 国内视频| 久久久精品94久久精品| 久久久久国产精品人妻一区二区| 国产爽快片一区二区三区| 国产一区二区三区综合在线观看| 国产精品一国产av| 国产精品女同一区二区软件| 久久精品久久久久久久性| 久久亚洲国产成人精品v| 欧美亚洲 丝袜 人妻 在线| 在线观看国产h片| 精品国产超薄肉色丝袜足j| 免费在线观看完整版高清| 亚洲精品视频女| 免费av中文字幕在线| 亚洲精品第二区| 1024香蕉在线观看| 青草久久国产| 高清不卡的av网站| 午夜精品国产一区二区电影| 午夜福利,免费看| 母亲3免费完整高清在线观看| 日本vs欧美在线观看视频| 99久久99久久久精品蜜桃| 自线自在国产av| 国产高清不卡午夜福利| 午夜老司机福利片| 午夜福利视频精品| 亚洲美女搞黄在线观看| 久久国产精品大桥未久av| 色婷婷av一区二区三区视频| 人妻人人澡人人爽人人| 亚洲欧美色中文字幕在线| 美女大奶头黄色视频| 欧美另类一区| 亚洲三区欧美一区| 亚洲成av片中文字幕在线观看| 亚洲一区二区三区欧美精品| 成人影院久久| 国产无遮挡羞羞视频在线观看| 亚洲久久久国产精品| 亚洲欧美成人综合另类久久久| 一本久久精品| 成年动漫av网址| 老鸭窝网址在线观看| 一级片'在线观看视频| 日韩免费高清中文字幕av| 成年av动漫网址| 亚洲av电影在线进入| 少妇人妻 视频| 精品国产国语对白av| 日本91视频免费播放| 高清黄色对白视频在线免费看| 一本—道久久a久久精品蜜桃钙片| 亚洲精品日本国产第一区| 赤兔流量卡办理| 日韩大码丰满熟妇| 久久97久久精品| 亚洲精品在线美女| 亚洲精品美女久久久久99蜜臀 | 亚洲av男天堂| 亚洲四区av| 美女中出高潮动态图| 麻豆av在线久日| 久久午夜综合久久蜜桃| 美女视频免费永久观看网站| 国产在线视频一区二区| 最新在线观看一区二区三区 | 精品第一国产精品| 老司机靠b影院| 亚洲男人天堂网一区| 人体艺术视频欧美日本| 最近最新中文字幕免费大全7| 色播在线永久视频| 精品人妻一区二区三区麻豆| 中文字幕人妻丝袜制服| 国产激情久久老熟女| 性色av一级| 国产成人精品福利久久| 多毛熟女@视频| av在线app专区| av免费观看日本| 秋霞在线观看毛片| a级毛片在线看网站| av国产久精品久网站免费入址| 亚洲成人免费av在线播放| 国产免费福利视频在线观看| 9色porny在线观看| 国产成人午夜福利电影在线观看| bbb黄色大片| 国产精品一区二区精品视频观看| 久久精品国产亚洲av涩爱| 亚洲av欧美aⅴ国产| 国产欧美亚洲国产| 美女主播在线视频| 又粗又硬又长又爽又黄的视频| 啦啦啦视频在线资源免费观看| 99久久人妻综合| 国产精品久久久久久人妻精品电影 | 久久国产精品大桥未久av| 成人漫画全彩无遮挡| 久久精品熟女亚洲av麻豆精品| 国产精品 国内视频| av国产精品久久久久影院| a级片在线免费高清观看视频| 国产精品免费大片| h视频一区二区三区| 日韩大片免费观看网站| 久久久国产欧美日韩av| 一本—道久久a久久精品蜜桃钙片| 九九爱精品视频在线观看| 国产一区二区三区av在线| 99热全是精品| 超色免费av| 最近2019中文字幕mv第一页| 美女大奶头黄色视频| 国产精品久久久久成人av| 亚洲av电影在线进入| 一二三四中文在线观看免费高清| 在线天堂中文资源库| 黑人欧美特级aaaaaa片| 午夜福利,免费看| 国产人伦9x9x在线观看| 欧美老熟妇乱子伦牲交| 久久热在线av| 亚洲国产精品999| 国产精品欧美亚洲77777| 天天影视国产精品| 亚洲av福利一区| 国产亚洲精品第一综合不卡| 最黄视频免费看| 国产精品99久久99久久久不卡 | 国产一区二区 视频在线| 老熟女久久久| 美女脱内裤让男人舔精品视频| 777米奇影视久久| 女的被弄到高潮叫床怎么办| 中文字幕人妻丝袜制服| 成人18禁高潮啪啪吃奶动态图| 日韩一区二区视频免费看| 高清av免费在线| 美女午夜性视频免费| 日韩中文字幕视频在线看片| 日韩 欧美 亚洲 中文字幕| 成年女人毛片免费观看观看9 | 你懂的网址亚洲精品在线观看| 国产男人的电影天堂91| 999精品在线视频| 亚洲第一区二区三区不卡| 最近中文字幕高清免费大全6| 国产成人系列免费观看| 人人妻人人澡人人看| 亚洲熟女毛片儿| 女性生殖器流出的白浆| 伦理电影免费视频| 国产精品香港三级国产av潘金莲 | av片东京热男人的天堂| 中文字幕亚洲精品专区| 人人妻人人澡人人看| 成年av动漫网址| 视频区图区小说| 18禁观看日本| 99热网站在线观看| 在线观看一区二区三区激情| av网站在线播放免费| 久久精品熟女亚洲av麻豆精品| 日韩av免费高清视频| 精品视频人人做人人爽| 国产视频首页在线观看| 亚洲成av片中文字幕在线观看| 精品一区二区三区av网在线观看 | 精品国产一区二区三区四区第35| 制服诱惑二区| 一边摸一边抽搐一进一出视频| 你懂的网址亚洲精品在线观看| 亚洲精品国产av成人精品| 中文字幕人妻丝袜制服| 亚洲欧美中文字幕日韩二区| 亚洲精品视频女| 日韩一区二区三区影片| av免费观看日本| 亚洲精品国产av蜜桃| 亚洲欧美成人综合另类久久久| av线在线观看网站| 国产精品亚洲av一区麻豆 | 日韩伦理黄色片| 久久久久网色| 欧美日韩一区二区视频在线观看视频在线| 亚洲精品国产av成人精品| 成年美女黄网站色视频大全免费| 欧美日韩精品网址| 国产精品偷伦视频观看了| 婷婷成人精品国产| 午夜老司机福利片| 欧美精品高潮呻吟av久久| 在线观看三级黄色| 午夜免费男女啪啪视频观看| 天堂8中文在线网| 日韩,欧美,国产一区二区三区| 亚洲,一卡二卡三卡| 伦理电影大哥的女人| 99久久99久久久精品蜜桃| 一级毛片黄色毛片免费观看视频| 下体分泌物呈黄色| 久久精品亚洲av国产电影网| 天天躁夜夜躁狠狠躁躁| 亚洲成国产人片在线观看| 精品一区二区三卡| 啦啦啦在线观看免费高清www| 午夜精品国产一区二区电影| 交换朋友夫妻互换小说| av卡一久久| 国产精品国产av在线观看| 日韩精品有码人妻一区| 秋霞在线观看毛片| 在线看a的网站| 9191精品国产免费久久| 在线精品无人区一区二区三| 国产精品三级大全| 91成人精品电影| 亚洲av成人精品一二三区| 精品亚洲成a人片在线观看| 中文精品一卡2卡3卡4更新| 日韩欧美一区视频在线观看| 久久人人爽av亚洲精品天堂| 午夜福利免费观看在线| 国产精品熟女久久久久浪| 日韩人妻精品一区2区三区| 建设人人有责人人尽责人人享有的| 国产精品偷伦视频观看了| 欧美日韩亚洲高清精品| 成人国产av品久久久| 欧美激情极品国产一区二区三区| 亚洲av电影在线观看一区二区三区| 精品一区二区三卡| 黑人欧美特级aaaaaa片| 人人妻人人澡人人爽人人夜夜| 日韩免费高清中文字幕av| 亚洲av欧美aⅴ国产| 国产精品一二三区在线看| 观看av在线不卡| 99久久人妻综合| 777米奇影视久久| 国产午夜精品一二区理论片| 国精品久久久久久国模美| 人体艺术视频欧美日本| 亚洲欧美精品综合一区二区三区| 男女之事视频高清在线观看 | 超碰97精品在线观看| 欧美精品人与动牲交sv欧美| 欧美国产精品va在线观看不卡| 成人免费观看视频高清| 成人三级做爰电影| 久久精品久久久久久久性| 日韩人妻精品一区2区三区| 成人影院久久| 精品亚洲乱码少妇综合久久| 欧美日韩视频精品一区| 亚洲精品第二区| 18在线观看网站| 日韩 亚洲 欧美在线| 欧美xxⅹ黑人| 咕卡用的链子| 久久国产精品大桥未久av| 五月天丁香电影| 久久久国产精品麻豆| 亚洲 欧美一区二区三区| 国产熟女欧美一区二区| 一级a爱视频在线免费观看| 亚洲熟女精品中文字幕| 免费黄频网站在线观看国产| 激情五月婷婷亚洲| 九草在线视频观看| 色婷婷久久久亚洲欧美| 如日韩欧美国产精品一区二区三区| 亚洲美女搞黄在线观看| av又黄又爽大尺度在线免费看| 亚洲视频免费观看视频| 午夜日韩欧美国产| 日韩人妻精品一区2区三区| 欧美成人精品欧美一级黄| 国产精品亚洲av一区麻豆 | 国产97色在线日韩免费| 欧美97在线视频| 91老司机精品| 免费观看人在逋| 亚洲精品第二区| 免费在线观看黄色视频的| tube8黄色片| 中国国产av一级| 90打野战视频偷拍视频| 国产精品一区二区精品视频观看| 操出白浆在线播放| 波多野结衣av一区二区av| svipshipincom国产片| 九草在线视频观看| 男人爽女人下面视频在线观看| 91aial.com中文字幕在线观看| 老熟女久久久| 2018国产大陆天天弄谢| 啦啦啦 在线观看视频| 久久精品亚洲熟妇少妇任你| 免费高清在线观看视频在线观看| 十八禁高潮呻吟视频| 国产片内射在线| 日韩熟女老妇一区二区性免费视频| 免费日韩欧美在线观看| 免费观看a级毛片全部| 女人高潮潮喷娇喘18禁视频| av一本久久久久| 成人免费观看视频高清| 最近2019中文字幕mv第一页| 亚洲av成人精品一二三区| 欧美黑人精品巨大| 少妇猛男粗大的猛烈进出视频| 大话2 男鬼变身卡| 国产一级毛片在线| 啦啦啦在线观看免费高清www| 91精品三级在线观看| 欧美 亚洲 国产 日韩一| 两性夫妻黄色片| 国产精品麻豆人妻色哟哟久久| 9色porny在线观看| 成年美女黄网站色视频大全免费| 蜜桃国产av成人99| 国产男人的电影天堂91| 国产欧美亚洲国产| a 毛片基地| 日韩精品有码人妻一区| 考比视频在线观看| 国产免费一区二区三区四区乱码| 丝袜人妻中文字幕| 97在线人人人人妻| 9色porny在线观看| 精品免费久久久久久久清纯 | 亚洲欧美精品自产自拍| 丝袜在线中文字幕| 日韩av免费高清视频| 国产成人精品福利久久| 九色亚洲精品在线播放| 日韩制服骚丝袜av| 免费看不卡的av| 欧美日韩国产mv在线观看视频| 汤姆久久久久久久影院中文字幕| 日日啪夜夜爽| 久久久久久久久免费视频了| 亚洲欧美精品自产自拍| 久久久久久人人人人人| 黑人猛操日本美女一级片| 欧美最新免费一区二区三区| 免费在线观看完整版高清| 国产在视频线精品| 午夜av观看不卡| 男女无遮挡免费网站观看| 国产免费视频播放在线视频| 色精品久久人妻99蜜桃| 高清在线视频一区二区三区| 久久综合国产亚洲精品| 制服丝袜香蕉在线| 悠悠久久av| √禁漫天堂资源中文www| 免费高清在线观看日韩| 建设人人有责人人尽责人人享有的| 亚洲精品乱久久久久久| 99国产精品免费福利视频| 亚洲一区中文字幕在线| 国产国语露脸激情在线看| 免费看不卡的av| 夜夜骑夜夜射夜夜干| 欧美人与性动交α欧美精品济南到| 久久久久久久国产电影| 一本色道久久久久久精品综合| 在线观看免费视频网站a站| 在线观看国产h片| 十八禁网站网址无遮挡| 中文字幕高清在线视频| 黑丝袜美女国产一区| 国产片特级美女逼逼视频| 这个男人来自地球电影免费观看 | 香蕉国产在线看| 韩国精品一区二区三区| 一本色道久久久久久精品综合| 精品一品国产午夜福利视频| 久久综合国产亚洲精品| 亚洲中文av在线| 不卡av一区二区三区| a级毛片在线看网站| 丝袜在线中文字幕| 色婷婷av一区二区三区视频| 成年动漫av网址| 青青草视频在线视频观看| 国产成人精品在线电影| 男人爽女人下面视频在线观看| 大香蕉久久成人网| 伊人久久大香线蕉亚洲五| 狂野欧美激情性xxxx| 尾随美女入室| 国产成人a∨麻豆精品| 国产精品av久久久久免费| 亚洲国产精品成人久久小说| 91国产中文字幕| 日日啪夜夜爽| 国产高清国产精品国产三级| 国产成人精品久久久久久| 九九爱精品视频在线观看| 9色porny在线观看| 欧美精品一区二区大全| 涩涩av久久男人的天堂| 免费观看性生交大片5| 午夜日本视频在线| 久久久国产欧美日韩av| 国产1区2区3区精品|