亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        拓?fù)渎曌优c聲子霍爾效應(yīng)?

        2017-12-05 02:34:56邢玉恒徐錫方張力發(fā)
        物理學(xué)報 2017年22期
        關(guān)鍵詞:聲子能帶對稱性

        邢玉恒 徐錫方 張力發(fā)

        (南京師范大學(xué)物理科學(xué)與技術(shù)學(xué)院,南京 210023)

        拓?fù)渎曌优c聲子霍爾效應(yīng)?

        邢玉恒 徐錫方 張力發(fā)?

        (南京師范大學(xué)物理科學(xué)與技術(shù)學(xué)院,南京 210023)

        (2017年9月28日收到;2017年11月6日收到修改稿)

        拓?fù)鋵W(xué)與物理的結(jié)合是近幾十年物理學(xué)蓬勃發(fā)展的一個新領(lǐng)域,它不僅活躍在量子場理論以及高能物理中,更廣泛地存在于凝聚態(tài)物理體系中,包括量子(反常、自旋)霍爾效應(yīng)和拓?fù)浣^緣體(超導(dǎo)體)等.聲子是凝聚態(tài)體系中熱輸運的主要載體;最近由于各種聲子器件的發(fā)現(xiàn),聲子學(xué)得到了廣泛的關(guān)注.本文介紹了聲子的拓?fù)湫再|(zhì)以及聲子的霍爾效應(yīng)現(xiàn)象,分別評述了在破壞時間反演對稱、破壞空間反演對稱、以及同時破壞時間和空間反演對稱三種情況下所產(chǎn)生的聲子霍爾效應(yīng)、聲子谷霍爾效應(yīng)等相關(guān)物理研究進(jìn)展.最后對拓?fù)鋵W(xué)在其他聲學(xué)體系中的應(yīng)用做了簡單介紹,并進(jìn)一步討論了其未來的發(fā)展方向.

        拓?fù)鋵W(xué),Berry相位,Berry曲率,聲子霍爾效應(yīng)

        1 引 言

        拓?fù)鋵W(xué)是近代發(fā)展起來的一個數(shù)學(xué)分支,主要研究幾何圖形或空間在連續(xù)改變形狀后還能保持不變的一些性質(zhì)的學(xué)科.早在十七世紀(jì)科學(xué)家就提出了“位置幾何學(xué)”和“相位分析”學(xué)說.1984年,英國布里斯托大學(xué)的科學(xué)家Michael Berry就提出了關(guān)于量子力學(xué)中的“幾何相”概念[1],主要描述了當(dāng)一個量子體系在參數(shù)空間沿著閉合路徑緩慢變化時,系統(tǒng)本征態(tài)所發(fā)生的變化,即系統(tǒng)的哈密頓量在參數(shù)空間發(fā)生絕熱演化時體系的本征態(tài)的變化.Berry相在拓?fù)湮锢韺W(xué)領(lǐng)域是一個重要的概念[2],Berry曲率由于揭示了量子力學(xué)中絕熱過程的奧秘且最近因其在量子計算機(jī)的實現(xiàn)中的重要用途而引起了科學(xué)家們的廣泛興趣[3].Berry曲率最早出現(xiàn)在帶有布洛赫電子的單帶的動力學(xué)中,在應(yīng)用于量子霍爾體系時,其半經(jīng)典理論可以解釋霍爾電流和量子霍爾電導(dǎo)率[4].近幾年來,Berry相效應(yīng)在解決鐵磁材料中的反?;魻栃?yīng)時起到了重要的作用[5,6],這也使Berry曲率在拓?fù)洳牧蠎?yīng)用中的研究得到科學(xué)家們的廣泛關(guān)注.由于在拓?fù)湮锢矸矫嬖瓌?chuàng)性的貢獻(xiàn),Thouless,Haldane和Kosterlitz被授予2016年諾貝爾物理學(xué)獎.拓?fù)浜臀锢淼那擅罱Y(jié)合,帶動了凝聚態(tài)物理的快速發(fā)展,比較有代表性的就是其近十年來逐漸興起的一個分支——拓?fù)浣^緣體[7,8],即在絕緣體中帶隙中存在拓?fù)浔Wo(hù)的邊界態(tài).

        由于熱二極管、熱晶體管、熱邏輯門等[9?11]的相繼提出,最近幾十年我們見證了聲子學(xué)的飛速發(fā)展,從而使得利用聲子進(jìn)行熱流控制和信息處理成為可能.由于聲子作為一種中性的準(zhǔn)粒子不能像電子一樣可以通過洛倫茲力直接和磁場耦合,因此當(dāng)Rikken等[12]科學(xué)家在順磁介質(zhì)樣品中發(fā)現(xiàn)聲子霍爾效應(yīng)現(xiàn)象時,物理學(xué)家們普遍感到很驚奇,即在有熱流通過的樣品薄膜的垂直方向上施加磁場時可以觀察到橫向熱流,由于電子對熱流貢獻(xiàn)可以忽略,所以把這種現(xiàn)象稱之為聲子霍爾效應(yīng).此實驗被Inyushkin和Taldenkov[13]所重復(fù).從那以后,科學(xué)家們提出了幾種理論來解釋這一現(xiàn)象[14].對于各種量子、自旋或反常霍爾效應(yīng)的電子輸運性質(zhì),拓?fù)淅碚撘呀?jīng)成功地解釋了一些潛在的機(jī)理,那么對于聲子的輸運性質(zhì)是否也可以運用拓?fù)淅碚搧磉M(jìn)行研究呢?

        2010年,Zhang等[15]將聲子的霍爾熱導(dǎo)率和聲子譜的Berry曲率聯(lián)系起來,揭示了聲子的拓?fù)湫再|(zhì),并發(fā)現(xiàn)了其中的拓?fù)湎嘧?從而開創(chuàng)了拓?fù)渎曌訉W(xué)的研究.其后,在2012年,Qin等[16]從能量磁化的角度對聲子霍爾熱導(dǎo)提出了修正.當(dāng)磁場作用在具有自旋聲子相互作用的體系上時,時間反演對稱性會發(fā)生破缺,從而產(chǎn)生聲子霍爾效應(yīng).在非磁性體系中,空間反演對稱破缺而時間反演對稱性得到保護(hù),可以觀察到谷聲子霍爾效應(yīng)[17].在2016年,Zhang研究組又繼續(xù)通過半經(jīng)典波包動力學(xué)方法統(tǒng)一研究了電子、聲子、磁子等熱霍爾效應(yīng),并得出了霍爾熱導(dǎo)率在拓?fù)淅碚摽蚣芟碌耐ㄓ霉絒18].最近Liu等[19]研究了時間反演對稱性和空間反演對稱性同時破缺的體系中的拓?fù)渎曌虞斶\現(xiàn)象.

        本文以Berry曲率為切入點,著重討論Berry相以及Berry曲率的基本理論與性質(zhì).通過研究在單獨破壞時間反演對稱、破壞空間反演對稱以及同時破壞時間反演對稱和空間反演對稱的聲子體系中Berry曲率的變化對聲子輸運的影響.在本文的第二部分介紹了Berry相的基本理論知識,了解Berry曲率與拓?fù)浣Y(jié)構(gòu)之間的關(guān)系;第三部分介紹了破壞時間反演對稱時所產(chǎn)生的聲子霍爾效應(yīng),研究其中的Berry曲率對拓?fù)浣Y(jié)構(gòu)性質(zhì)的一些影響;第四部分介紹了破壞空間反演對稱后形成的聲子谷霍爾效應(yīng);第五部分主要介紹同時破壞時間和空間反演對稱時對拓?fù)浣Y(jié)構(gòu)的影響;第六部分簡單介紹了Berry相效應(yīng)在其他聲學(xué)體系中的應(yīng)用;第七部分對拓?fù)渎曌訉W(xué)的應(yīng)用做了簡單展望.

        2 Berry相及拓?fù)浣Y(jié)構(gòu)

        Berry相位由三個至關(guān)重要的部分構(gòu)成,這三個關(guān)鍵的特性決定了其在拓?fù)鋵W(xué)中的重要地位.第一,Berry相位是規(guī)范不變的,其本征波函數(shù)由線性齊次方程定義,這個性質(zhì)使得Berry相位更加的物理化,早期的實驗研究也是通過這個性質(zhì)進(jìn)行的.第二,Berry相位具有獨特的幾何性質(zhì),Berry相可以表示為在閉合參數(shù)空間的線性積分,且與實際路徑的變化率無關(guān),Berry曲率的這種獨特的性質(zhì)使得其得到了廣泛的應(yīng)用.第三,Berry相位與規(guī)范場論和微分幾何有著密切的聯(lián)系,這使得Berry相成為一種直觀且意義十分重大的物理概念,尤其在今天的霍爾效應(yīng)以及拓?fù)湮锢韺W(xué)的研究中具有重要的意義[7].

        2.1 絕熱循環(huán)演化

        波函數(shù)的整體相位不會對物理觀測造成影響,然而量子態(tài)之間的相對相位與很多的物理現(xiàn)象有關(guān)系.1984年,英國布里斯托大學(xué)的科學(xué)家Michael Berry提出了關(guān)于量子力學(xué)中的“幾何相”概念,為了更好地詮釋這個相位的物理意義,以下對這個相位進(jìn)行簡單的推導(dǎo).

        當(dāng)物理系統(tǒng)沿著參數(shù)空間中任意閉合回路做絕熱演化,其非簡并本征態(tài)的變化除了一個額定的相位之外,必然回到其本身,該相位就稱為Berry相位.設(shè)體系的哈密頓量算符是一組參量的R(R1,R2,···,RD)函數(shù):

        而R隨著時間做周期性變化

        R(t)的周期變化在參量空間定義了一條閉合曲線C,假設(shè)此時周期T足夠大,以致哈密頓算符隨時間的變化非常緩慢,從而使得這一過程被稱為絕熱演化過程,致使系統(tǒng)在每一瞬間都是靜止的.于是,對于某一瞬時t,瞬時定態(tài)薛定諤方程成立.在絕熱條件下,瞬時本征波函數(shù)的含時薛定諤方程為

        乘以〈m(R(t))|化簡后得

        考慮一般含時薛定諤方程

        在絕熱近似條件下,將本征函數(shù)歸一化條件求導(dǎo)可得

        所以絕熱近似下式含時薛定諤方程的解為

        在引入R(t)的空間“矢勢”后,Berry相位可以寫成

        其中Bm(R)=?R×Am(R)稱為參數(shù)空間的“磁場強(qiáng)度”.

        從上式可以看到:第一,Berry相位γm是實量;第二,Berry相位是一個規(guī)范不變的量,即在物理上是一個可以觀測的量;第三,Berry相位不依賴于系統(tǒng)參數(shù)演化的快慢,只與演化的路徑有關(guān);第四,對于一維的參數(shù)空間,Berry相位為零,但是對于高維的參數(shù)空間,Berry相位不為零;第五,Berry相位是不可積的,也就是說它不能寫成一個關(guān)于參數(shù)R的連續(xù)函數(shù)[20].故Berry相位僅與閉合回路或者閉合曲面的幾何性質(zhì)有關(guān),與參數(shù)隨時間的變化率無關(guān),因此,Berry相位又稱為幾何相位.

        2.2 Berry曲率

        通過(8)式以及下面恒等式:

        所以Berry曲率可以寫成本征態(tài)的求和[21]:

        在這里,通常將A稱為Berry聯(lián)絡(luò)(Berry connection),將B稱為Berry曲率(Berry curvature).

        在上述的公式推導(dǎo)過程中,最重要的就是絕熱假設(shè),即參數(shù)變化的足夠緩慢從而使體系一直約束在某個能級上.也就是說此時的取值空間已經(jīng)從整個完備的希爾伯特空間投影到某一個特定的能級上,這就是產(chǎn)生Berry相位的本質(zhì).

        2.3 Berry曲率與對稱性的關(guān)系

        上述已經(jīng)寫出的Berry曲率的表達(dá)式

        滿足時間反演對稱時

        滿足空間反演對稱時

        若同時滿足空間和時間反演對稱

        2.4 布洛赫能帶中的Berry相位

        以上介紹了用哈密頓算子描述通用體系Berry相位的基本概念,現(xiàn)在考慮它在晶體中的實現(xiàn).我們知道,晶體的能帶結(jié)構(gòu)為研究Berry相效應(yīng)提供了一個天然的平臺.在孤立電子近似下,首先可以寫出單電子情況下的哈密頓量:

        其中V(r+a)=V(r),a是布拉伐晶格的格矢,通過布洛赫定理,周期性的哈密頓量的本征態(tài)滿足以下條件:

        這里n是能帶的數(shù)目,?q是晶體的動量,通過進(jìn)一步變換可得

        在這里,unq(r)=e?iq·rψnq(r) 是布洛赫函數(shù)的周期性部分,滿足嚴(yán)格的周期性邊界條件unq(r+a)=unq(r).Berry相位以及Berry曲率可以寫成:

        3 聲子霍爾效應(yīng)

        3.1 霍爾效應(yīng)的發(fā)現(xiàn)

        霍爾效應(yīng)是電磁效應(yīng)的一種,這一現(xiàn)象是美國物理學(xué)家霍爾于1879年在研究金屬的導(dǎo)電機(jī)理時發(fā)現(xiàn)的:即當(dāng)電流垂直于外磁場通過導(dǎo)體時,載流子發(fā)生偏轉(zhuǎn),在垂直于電流和磁場的方向會產(chǎn)生一附加電場,從而在導(dǎo)體的兩端產(chǎn)生電勢差,這個電勢差也被稱為霍爾電勢差.

        圖1 (a)聲子霍爾效應(yīng)示意圖,熱流方向在磁場作用下發(fā)生偏轉(zhuǎn)[24];(b)在Tb3Ga5O12薄膜樣品測得的橫向溫差隨磁場大小和方向的變化關(guān)系[12]Fig.1.(a)Phonon Hall e ff ect:the direction of the heat flow is de flected by the magnetic field[24];(b)magnetotransverse temperature di ff erence in a nonoriented sample of Tb3Ga5O15for heat currents perpendicular(circles)and parallel to the field(squares)[12].

        近些年來,半導(dǎo)體材料和低維物理學(xué)的發(fā)展推動了許多霍爾效應(yīng)的發(fā)展.在電絕緣的晶體中,聲子是熱流的載體,由于聲子不能直接和磁場相互作用,所以通常認(rèn)為不存在聲子霍爾效應(yīng).最近,Strohm等[12]在順磁介質(zhì)Tb3Ga5O12薄膜樣品中發(fā)現(xiàn)了聲子霍爾效應(yīng)現(xiàn)象,如圖1所示,即在有縱向溫度梯度的順磁絕緣體薄膜上,加上垂直于薄膜平面的磁場后,測得橫向溫度差.當(dāng)磁場方向平行于縱向熱流時,橫向溫度差可以忽略,然而當(dāng)磁場方向垂直于縱向熱流(縱向溫差為1 K)時,可以觀察到200μK的溫度差.由于電子不參與熱輸運,所以Strohm等將之命名為聲子霍爾效應(yīng).此后該實驗并被Inyushkin和Taldenkov[13]所重復(fù).科學(xué)家們試圖提出幾種理論來解釋這一現(xiàn)象[14,22],但是相對復(fù)雜并且有相互矛盾之處.另外,文獻(xiàn)[23]利用非平衡態(tài)格林函數(shù)方法研究了具有自旋聲子相互作用的四終端的二維薄膜的聲子霍爾溫差,從另一角度解釋了聲子霍爾效應(yīng).Zhang等[15]的工作從拓?fù)涮匦缘慕嵌仁状螌β曌踊魻栃?yīng)進(jìn)行了揭示,下面首先從聲子霍爾體系的哈密頓量講起.

        3.2 哈密頓量以及二次量子化

        在有磁場存在的情況下[14],離子晶體中晶格格點的動能可以寫成

        所以,上式可以寫成

        如果此時沿z方向施加強(qiáng)度為B的磁場,前提是只考慮二維情況下的運動(即只沿x和y方向),此時離子α的動能可以寫成:

        之前的理論中,科學(xué)家們考慮到拉曼相互作用的影響,基于量子理論和基本對稱性的自旋-聲子相互作用[26?28],得到了HI=h.(u×p),這里h=gcM,M與磁場B成正比.如果磁場沿z方向的作用,此時自旋-聲子相互作用可以寫成所以整個晶格的哈密頓量可以寫成[29]

        這和(25)式類似.然而此哈密頓量并不總是正定的,我們應(yīng)該選用(25)式,其始終是正定的.

        哈密頓量是關(guān)于u和p二次型方程,可以寫出線性運動方程:

        上式并不是標(biāo)準(zhǔn)的本征值問題,重新定義一個新的向量x=(μ,ε)T.然后,通過布洛赫定理,(28)式可以寫成[15]:

        I是單位矩陣.因此,上式的本征值問題可以寫成

        在這里,由于該有效哈密頓量是非厄米的,需要計(算右本征矢量為xk和左本征矢量對于本征矢εk的歸一化條件即為

        對于上式中的本征值問題,可知完備基包括負(fù)的分支等.從聲子霍爾效應(yīng)的拓?fù)湫再|(zhì)可知[15],聲子霍爾熱導(dǎo)率的公式可以寫成包括正負(fù)頻的所有分支的貢獻(xiàn)形式,為了簡化符號,對于所有的分支,定義含時算子變換為變換關(guān)系為故由以上可以得到[15]:

        這里的f(ωk)為玻色-愛因斯坦分布方程.

        在上述變化后,位移和動量運算符可以用二次量子化算符寫成[30]:

        這里|ωk|=ωksign(σ).通過上述,系統(tǒng)的哈密頓量可以寫成

        通過一系列變換后,哈密頓量最終可以化為[30]

        系統(tǒng)的熱流可以寫為[30]

        這里V是N個基本原胞的總體積.通過對上述式的求導(dǎo)并引入Berry聯(lián)絡(luò)得出最終的熱流公式:

        通過格林-庫珀公式,可以計算得到霍爾熱導(dǎo)線性響應(yīng)貢獻(xiàn)κKubo:

        其中:

        是σ?guī)У溅摇鋷У膶erry曲率的貢獻(xiàn).另外我們還必須考慮能量磁化貢獻(xiàn)ME[16],總的霍爾熱導(dǎo)率可以寫為

        最終可以得到聲子熱導(dǎo)率公式[16]

        這里[16]

        從上述公式中可以看出聲子的霍爾熱導(dǎo)率由其色散關(guān)系和Berry曲率所決定,這與量子霍爾電導(dǎo)有一定的類比性[4].

        此外,通過半經(jīng)典的波包動力學(xué)方法,可以研究熱霍爾效應(yīng),包括電子熱霍爾效應(yīng)、磁振子霍爾效應(yīng)以及聲子霍爾效應(yīng).在熱霍爾效應(yīng)中,可以定義粒子磁矩為

        故體系的能量和熱磁矩分別為

        這里εn是布洛赫能帶的能量,μ是化學(xué)勢.類比于電子軌道磁化,熱軌道磁化MQ有兩個部分,一個來自于波包的自轉(zhuǎn)有關(guān),記為;另一個是波包中心的運動的貢獻(xiàn),將其表示為[18].

        可以得到總的熱流為

        此公式是電子霍爾熱導(dǎo)、磁振子霍爾熱導(dǎo)以及聲子霍爾熱導(dǎo)的通用公式.其中對于聲子可以進(jìn)一步寫成:

        其中n,d分別是原胞內(nèi)的原子數(shù)和每一個原子的自由度.

        3.3 聲子霍爾效應(yīng)的拓?fù)湫再|(zhì)

        由于本證矢量的波矢依賴性,可以計算其Berry相效應(yīng).通過 Berry近似,首先設(shè)x(t)=, 并代入上式中. 由于 Berry相,所以Berry曲率可寫成

        然后通過將第一布里淵區(qū)域上的Berry曲率整合得到相關(guān)的拓?fù)潢悢?shù):

        如果沒有拉曼自旋聲子的影響,也就是說當(dāng)h=0時,此時 Berry曲率處處為零,聲子的霍爾熱導(dǎo)率消失.當(dāng)外加磁場作用時,Berry曲率不為零,因此,聲子霍爾效應(yīng)也隨之出現(xiàn).另外研究發(fā)現(xiàn),如果體系出現(xiàn)對稱性滿足SDS?1=D,SAS?1=?A,此時聲子的霍爾熱導(dǎo)為零[31],此時如果磁場改變了,Berry曲率也會改變,但是陳數(shù)在很大的范圍內(nèi)并不改變.盡管陳數(shù)被量化了,但是由于(40)式中f(ωσ)(ωσ+ωσ′)2項的存在,聲子的霍爾熱導(dǎo)率并沒有被量化.

        圖2 (a)—(d)當(dāng)臨界磁場hc? =hc?10?2rad/ps時能帶1到4的Berry曲率輪廓圖;(e)—(h)當(dāng)臨界磁場hc+=hc+10?2rad/ps時能帶1到4的Berry曲率輪廓圖;從(a)—(h)水平軸和垂直軸分別對應(yīng)于波矢kx和ky;(i)在不同磁場中的Berry曲率,實線和虛線分別對應(yīng)于在臨界磁場為hc?時?2和?3,相應(yīng)的實點和虛點分別對應(yīng)于臨界磁場為hc+時?2和?3;(j)四種能帶的陳數(shù)C1(實線),C2(虛線),C3(點線),C4(虛點線);(k)在臨界磁場為hc附近的不同磁場下能帶2和能帶3的色散關(guān)系圖,其中虛線、實線以及點線分別對應(yīng)于在臨界磁場為hc?,hc,hc+時的帶[15]Fig.2.(a)–(d)The contour map of Berry curvatures for bands 1–4 at hc? =hc? 10?2rad/ps;(e)–(h)the contour map of Berry curvatures for bands 1–4 at hc+=hc+10?2rad/ps;for(a)–(h),the horizontal and vertical axes correspond to wave vector kxand ky;(i) ? at di ff erent magnetic fields,the solid and dashed lines correspond to ?2and ?3at hc?,respectively,while dotted and dash-dotted lines correspond to those at hc+;(j)Chern numbers of four bands,C1(solid line),C2(dashed line),C3(dotted line),and C4(dashdotted line);(k)the dispersion relation of bands 2 and 3 at di ff erent magnetic fields in the vicinity of hc,the dashed,solid and dotted lines correspond to the bands at hc?,hcand hc+,respectively[15].

        隨著磁場的增加,在臨界磁場hc的附近,發(fā)現(xiàn)存在拓?fù)湎嘧?不同頻帶在靠近臨界磁場的Berry曲率由圖2(a)—(f)給出,當(dāng)hc附近磁場變化非常小時,在Γ點附近的能帶2和能帶3的Berry曲率幾乎是完全不同的,而能帶1和能帶4卻保持不變.如圖2(i)所示,在臨界磁場上方和下方Berry曲率發(fā)生了急劇的變化,從而導(dǎo)致能帶2和3的陳數(shù)發(fā)生跳躍,如圖2(j)所示,這個跳躍表明兩條帶的拓?fù)浣Y(jié)構(gòu)在臨界磁場下突然改變,這與能帶結(jié)構(gòu)的相變有著重要的聯(lián)系[15].

        為了進(jìn)一步研究聲子拓?fù)鋷ЫY(jié)構(gòu)變化的機(jī)理,同樣給出了在臨界磁場附近的聲子的色散關(guān)系圖,如圖2(k)所示,可以看到當(dāng)磁場增加到hc時,在Γ點附近能帶2和能帶3彼此的相互靠近.同樣,在臨界磁場區(qū)域,兩條能帶發(fā)生了簡并且兩條帶的形狀出現(xiàn)了錐形.而在臨界點hc上方,兩條能帶發(fā)生了分裂.也就是說在臨界磁場附近的聲子的色散關(guān)系直接影響了Berry曲率,從而導(dǎo)致拓?fù)湎嘧?

        這里提出的聲子霍爾熱導(dǎo)率的拓?fù)鋵W(xué)方法是一般的,可以應(yīng)用于低溫下的真實材料.最近的研究發(fā)現(xiàn)它也可以用于磁霍爾效應(yīng)[32].從拓?fù)湫再|(zhì)和色散關(guān)系來解釋的聲子霍爾效應(yīng)相變也可以推廣到研究其他霍爾效應(yīng)或非平衡運輸中的相變.

        4 聲子谷霍爾效應(yīng)

        4.1 谷電子霍爾效應(yīng)

        由于空間反演對稱性的破壞,動量空間里分開的不等價的能谷K,K′,通過Berry曲率或軌道磁矩可以來描述這種不等價性.因此這兩處不等價的能谷便構(gòu)成了除電荷和自旋以外電子的另一個自由度,這也導(dǎo)致了能谷電子學(xué)的出現(xiàn).

        自然界大量存在類似于石墨烯的六角蜂窩狀晶格層狀材料,例如氮化硼、硒化鎵、過渡金屬二硫化物等.不同于石墨烯,這些材料中的A和B子格是由不同的原子構(gòu)成,從而體系破壞了空間反演對稱性.這些材料也為研究能谷電子學(xué)提供了條件.

        在A/B堆積的雙層石墨烯材料中,通過施加垂直與二維平面的電場來打破石墨烯的空間反演對稱性,從而在K和K′能谷處產(chǎn)生大小相等但反向的并且與能谷相關(guān)的Berry曲率及軌道磁矩[33].因此也帶來了一系列新奇的效應(yīng),如能谷霍爾效應(yīng)[34].如圖3,在非零的Berry曲率和外加面內(nèi)電場的作用下,谷電子獲得反常的大小與Berry曲率成正比[35,36]的橫向速度,不同的谷電子偏轉(zhuǎn)方向相反.因此在樣品兩端可探測到谷極化電流,圖中Berry曲率的反差使得K和K′能谷電子在垂直外加電場方向的運動發(fā)生橫向相反方向的偏轉(zhuǎn),這就是谷電子霍爾效應(yīng).最近,實驗上已經(jīng)在單層二硫化鉬三極管[37]和石墨烯超晶格[38]中觀測到了谷電子霍爾效應(yīng).

        通過對稱性分析可知,破壞時間反演對稱性或破壞空間反演對稱性的系統(tǒng)中聲子攜帶非零的角動量.在這一部分,我們引入非磁性的六角AB晶格體系,即破壞空間反演對稱性,但保留時間反演對稱性.在兩個不等價的能谷處手性的聲子大量集中,即谷上的聲子是左旋或者右旋圓極化的,且非簡并的聲子模可以用量子的贗角動量標(biāo)記[17].所以,類比谷電子霍爾效應(yīng),對于谷聲子,自然要提出兩個問題:谷聲子的Berry曲率是否為零?如果Berry曲率是非零的,是否存在谷聲子霍爾效應(yīng)?

        圖3 谷電子霍爾效應(yīng)示意圖[34]Fig.3.Valley electronic Hall e ff ect diagram[34].

        4.2 聲子Berry曲率

        之前計算聲子的Berry曲率,采用傳統(tǒng)的晶格動力學(xué)周期解

        這里的u和ε是維度為2(二維運動)或3(三維運動)的列向量.接著可以得到運動方程

        而實際上,為了準(zhǔn)確地得到聲子的Berry曲率,必須考慮到原胞內(nèi)格點的相位差,因此這里采取布洛赫波形式的聲子波函數(shù)

        其中dα是第α原子相較于第l個原胞的平衡位置,因此可以得出

        所以動力學(xué)矩陣變?yōu)?/p>

        此時的運動方程為

        這與之前的運動方程是等價的,它們具有相同的色散關(guān)系,僅僅是波函數(shù)相差了一個相位.

        在第二節(jié)已經(jīng)計算過聲子的Berry曲率,這里由于具有時間反演對稱性,只需考慮所有正的聲子模.

        由于ε(k,σ)=ε?(?k,?σ), 在時間反演對稱性的保護(hù)下, 有ε(k,σ)=ε?(?k,σ), 所以因此

        類比于電子的Berry曲率計算Hψn=Enψn,聲子的Berry曲率為

        這一形式也與第二節(jié)所給的聲子Berry曲率相符合.

        4.3 谷聲子霍爾效應(yīng)

        非磁性的六角A/B晶格體系破壞了空間反演對稱性,如圖4(a),在能谷處觀察到了非零的聲子Berry曲率.能帶1和能帶4在能谷處擁有較大的Berry曲率,而能帶2和能帶3在能谷處的Berry曲率較小.正因為能谷處存在非零的Berry曲率,施加縱向的應(yīng)變梯度Estrain,類比于谷電子霍爾效應(yīng),會出現(xiàn)與應(yīng)變梯度和Berry曲率的叉矢成正比的反常速度,νanom∝?Estrain×B.因此,被左旋或者右旋極化光子激發(fā)的位于不同能谷處的聲子會沿橫向移動.如圖4(b)和圖4(c),如果光子的極性反轉(zhuǎn),那么橫向聲子流的方向也會翻轉(zhuǎn).隨著聲子在某條邊的積累,可以測得橫向的溫度差.很顯然,如果圓極化的受激光子極性反轉(zhuǎn),則橫向溫差也會相反.

        圖4 六角A/B晶格的聲子貝里曲率和谷聲子霍爾效應(yīng) (a)能帶1(底部二維圖)和能帶2(頂部三維圖)的Berry曲率;(b),(c)谷聲子被右旋或左旋極化光激發(fā)(霍爾熱流由紅色箭頭標(biāo)出),在應(yīng)力梯度下產(chǎn)生相應(yīng)的聲子[17]Fig.4.Phonon Berry curvature and valleyphonon Hall e ff ect in a honeycomb of A/B lattice:(a)Berry curvatureof band 1(bottom contour plot)and band 2(top 3D plot);(b),(c)schematic of the valley phonon hall e ff ect(the Hall current denoted by the olive curve arrows)undera strain gradient(the orange arrows),where valley phononsare excited by a ray of right-handed or left-handed polarizedlight(the red wave lines)[17].

        聲子霍爾效應(yīng)已在實驗上證實[39],在順磁絕緣體中,磁場可以改變聲子的輸運,因此可以測得橫向的溫度差,但是其效應(yīng)受限于及其微弱的自旋聲子相互作用.在破壞空間反演對稱性而保持時間反演對稱性的體系中,非零的Berry曲率誘發(fā)的谷聲子霍爾效應(yīng)確實廣泛存在于六角晶格中,將會有更大范圍的應(yīng)用.

        5 同時破壞時間和空間反演對稱的聲子模型

        由于時間反演對稱破缺會產(chǎn)生具有非零拓?fù)潢悢?shù)的聲子態(tài),即聲子霍爾效應(yīng),其具有單向?qū)ǖ牟皇苌⑸涞穆曌舆吔缒J?空間反演對稱破缺會產(chǎn)生具有非零的Berry相的谷,可用作調(diào)控聲子輸運的新型量子自由度,即谷聲子霍爾效應(yīng);如果兩種對稱性破缺同時存在,則能產(chǎn)生新的拓?fù)湎嘧儾⒀萆鲐S富的聲子拓?fù)淞孔討B(tài).

        圖5 (a)由A/B晶格組成的蜂窩狀晶格示意圖,以及在K和K′聲子谷附近的狄拉克錐;(b)當(dāng)同時破壞掉時間反演和空間反演對稱兩個谷的頻率不變,其中A/B晶格聲子態(tài)A由紅色標(biāo)出,B由藍(lán)色標(biāo)出[19]Fig.5.(a)A honeycomb lattice composed of A/B sublattices,and the Dirac cones of phonons near the K and K′valleys;(b)the two valleys that arenondegenerate in frequency caused by breaking P and T simultaneously[19].

        Liu等[19]的最新研究中,利用類似于Haldane蜂窩晶格模型,研究并探索Berry相與聲子的拓?fù)湫?yīng).在這種晶格中,縱向光學(xué)和聲學(xué)模式在K處形成線性交叉條帶和 Dirac點,如圖5所示.當(dāng)破壞空間反演對稱性以及時間反演對稱時,哈密頓量可以寫成:

        在上述的介紹中,我們已經(jīng)了解了單獨破壞時間反演對稱或者空間反演對稱對Berry曲率的影響,此時Berry通量在每個谷中的值為±π,且陳數(shù)CK(K′)=±1/2. 由于C=CK+CK′,故在單獨破壞時間反演時,Berry曲率B(k)=B(?k),陳數(shù)為C=±1,此系統(tǒng)中可以觀測到聲子霍爾效應(yīng);在單獨破壞空間反演對稱時,B(k)=?B(?k),陳數(shù)C=0,可以觀測到谷聲子霍爾效應(yīng).如圖5所示,當(dāng)同時破壞時間反演和空間對稱性時,兩聲子谷之間的頻率簡并消失,K和K′之間的能帶帶隙打開且相互獨立,即在極端情況下,只能在一個谷中實現(xiàn)零帶隙能帶,這種發(fā)現(xiàn)可以應(yīng)用在聲子谷濾波器中.

        6 Berry相效應(yīng)在其他聲學(xué)體系中的應(yīng)用

        近幾年,關(guān)于Berry曲率在聲學(xué)體系中的應(yīng)用的研究很多,例如在低溫的情況下,磁絕緣體中的自旋激發(fā)的熱導(dǎo)率可超過聲子的自發(fā)激發(fā)的熱導(dǎo)率,但由于其電荷為中性的,故自旋波不會產(chǎn)生相應(yīng)的熱霍爾效應(yīng).然而,在戈薇晶格,Berry曲率與熱導(dǎo)率之間有著緊密的聯(lián)系,在實驗中觀察到κ隨著溫度或磁場的變化而發(fā)生符號反轉(zhuǎn),其與金屬帶之間的陳數(shù)通量的符號交替相關(guān)[40].其次,在電子和聲子相互作用的系統(tǒng)中,量子霍爾效應(yīng)與無耗散電流情況下霍爾電導(dǎo)率之間也有著緊密的聯(lián)系,最近已經(jīng)觀察到,當(dāng)電流增加到一個臨界值時,對應(yīng)于聲速的載波漂移速度出現(xiàn)了一個戲劇性的耗散[41].在反轉(zhuǎn)對稱性破壞的非磁性系統(tǒng)中,Berry曲率誘導(dǎo)橫向谷聲子霍爾效應(yīng)將產(chǎn)生一些新的應(yīng)用.

        對于谷贗自旋,由于其作為新型信息載體的潛力巨大且在動量空間中能標(biāo)注能量極值的量子狀態(tài)而引起了人們的關(guān)注[37,38,42,43].在最近的一項研究中,武漢大學(xué)劉正猷研究組從實驗上在聲波晶體中觀察到了聲音的拓?fù)涔容斶\現(xiàn)象[44],他們除了通過聲場的空間掃描直接觀察谷選擇性邊緣模式之外,還在急劇彎曲的界面中觀察到抗反射現(xiàn)象.與傳統(tǒng)聲波導(dǎo)相比,拓?fù)浔Wo(hù)的聲音傳輸與其明顯不同,這一發(fā)現(xiàn)可以作為設(shè)計具有非常規(guī)功能設(shè)備的基礎(chǔ).

        另外,近期劉正猷教授研究組在聲子晶體中首次引入谷自由度的概念,預(yù)測了聲學(xué)谷態(tài)所具有的渦旋聲場分布并提出了產(chǎn)生特定手性聲渦旋場的方法,在此基礎(chǔ)上,他們進(jìn)一步實現(xiàn)了谷束分離和基于谷態(tài)的粒子操控.可以預(yù)期其潛在的應(yīng)用與聲學(xué)渦旋及物質(zhì)的異乎尋常的相互作用,例如觸發(fā)俘獲的微粒的旋轉(zhuǎn)而不接觸.這一發(fā)現(xiàn)具有重要的科學(xué)意義[45].

        此外,通過剛性鍵或者彈簧連接起來的質(zhì)點組成的框架或格子模型構(gòu)造在結(jié)構(gòu)工程、建筑和材料科學(xué)等不同領(lǐng)域有著重要的應(yīng)用.在這類晶格中,剛性鍵的數(shù)量與格子自由度之間的差異決定了它們的“零頻率”軟盤模式的數(shù)量.當(dāng)這些平衡時,系統(tǒng)處于機(jī)械不穩(wěn)定的邊緣,被稱為等靜壓.由于某些擴(kuò)展的等靜壓晶格在其邊界的軟盤模式呈現(xiàn)局部化,這些邊界模式對局部擾動不敏感,似乎具有拓?fù)淦鹪?可聯(lián)想到在量子霍爾效應(yīng)和拓?fù)浣^緣子中發(fā)生的受保護(hù)的電子邊界模式.由此科學(xué)家們建立了拓?fù)錂C(jī)械模式和電子系統(tǒng)的拓?fù)鋷Ю碚撝g的聯(lián)系,并預(yù)測了具有不同邊界模式的新的拓?fù)浯篌w積機(jī)械相的存在[46].

        最近,科學(xué)家在機(jī)械拓?fù)浣^緣子中觀察到聲子的螺旋邊緣態(tài),這一發(fā)現(xiàn)使得利用表面聲子的穩(wěn)定性設(shè)計可靠的波導(dǎo)的拓?fù)渎晫W(xué)材料成為可能[47].在具有單向彈性邊緣波的拓?fù)渎曌泳w中,科學(xué)家們通過使用陀螺慣性效應(yīng)來打破時間反演對稱,實現(xiàn)電子量子(反常)霍爾效應(yīng)的聲子模擬,研究結(jié)果激發(fā)了新型表面波器件的設(shè)計熱潮,將廣泛應(yīng)用于電子、電信和聲學(xué)成像[48].

        7 展 望

        拓?fù)湮锢韺W(xué)是一個方興未艾的領(lǐng)域,有著很廣泛的前景等待著科學(xué)家們?nèi)ラ_發(fā).自從量子霍爾效應(yīng)發(fā)現(xiàn)以來,凝聚態(tài)物理中的拓?fù)洮F(xiàn)象層出不窮,如拓?fù)浣^緣體、拓?fù)涑瑢?dǎo)體等.拓?fù)渑c聲子學(xué)的結(jié)合才剛剛開始,已有的相關(guān)研究工作中已可看到很大的發(fā)展空間和廣闊的應(yīng)用前景.

        在磁性體系中,聲子霍爾效應(yīng)的發(fā)現(xiàn)給聲子輸運與控制提供了新的方法,但至今沒有更多的實驗研究聲子霍爾效應(yīng).由于絕大部分順磁材料自旋聲子相互作用很小,所以測得霍爾效應(yīng)很弱,如何尋找更大效應(yīng)的順磁材料是一個挑戰(zhàn).另外鐵磁絕緣體,由于更大的磁化率和更大的自旋聲子相互作用,可以存在更強(qiáng)的的聲子霍爾效應(yīng),但是由于其中的磁振子霍爾效應(yīng)也很強(qiáng),所以區(qū)分鐵磁絕緣材料中熱霍爾效應(yīng)中磁振子和聲子的貢獻(xiàn)也是一個很有意義的課題,由于磁振子和聲子的拓?fù)湫再|(zhì)不一樣,拓?fù)渎曌訉W(xué)在鐵磁材料中聲子霍爾效應(yīng)可以發(fā)揮更大的作用.

        塊體能譜的拓?fù)涮匦詫嶋H上也對應(yīng)著非平凡的邊界態(tài),那么對于這里的聲子霍爾效應(yīng),理論上應(yīng)該存在相應(yīng)的拓?fù)浔Wo(hù)的邊界態(tài)或者邊界模式,和磁振子一樣聲子等玻色子體系中都沒有費米面的概念,所以真正意義上的玻色子拓?fù)浣^緣體不存在,但是可以通過波導(dǎo)選擇帶隙中的模式,這樣磁振子和聲子拓?fù)浣^緣體也可以存在.拓?fù)浯耪褡咏^緣體已經(jīng)發(fā)現(xiàn),但是拓?fù)渎曌咏^緣體還沒有真正發(fā)現(xiàn),尋找拓?fù)浔Wo(hù)的邊界聲子模對于無耗散的聲子輸運控制及相關(guān)能量科學(xué)將具有潛在的應(yīng)用價值.

        谷聲子霍爾效應(yīng)為聲子的控制提供了一個新的自由度.由于電子的谷間散射與谷聲子緊密聯(lián)系,如最近發(fā)現(xiàn)的谷聲子具有確定的手性并且決定電子谷間散射的選擇定則[16],因而谷聲子對谷電子學(xué)的發(fā)展意義重大.對于谷聲子相關(guān)理論,以及不同模型和材料中的手性聲子計算還需廣泛深入的研究.另外實驗直接觀測谷聲子以及谷聲子霍爾效應(yīng)也是一項具有重大意義的研究.

        對于同時破壞時空對稱性的體系,谷聲子的能量簡并打開,為實現(xiàn)單谷聲子的選擇與應(yīng)用提供了思路.而如何找到這樣的材料還有待進(jìn)一步的探索研究.

        拓?fù)渎曌訉W(xué)才剛剛開始.而拓?fù)渑c自旋電子學(xué)的結(jié)合已經(jīng)引領(lǐng)了凝聚態(tài)物理幾十年.我們相信,由于聲子在凝聚態(tài)物理中廣泛存在,如超導(dǎo)、電聲相互作用、布里淵區(qū)及拉曼散射、熱電效應(yīng)、熱效應(yīng)等,拓?fù)渎曌釉谶@些領(lǐng)域中不可回避并起到關(guān)鍵作用,也將引領(lǐng)凝聚態(tài)物理學(xué)的蓬勃發(fā)展.

        [1]Berry M V 1984Math.Phys.Sci.392 45

        [2]Wilczek F,Shapere A 1989Geometric Phases Phys.5 05857

        [3]Prabhakar S,Melnik R,Bonilla L L 2014Phys.Rev.B89 245310

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

        [5]Karplus R,Luttinger J M 1954Phys.Rev.95 1154

        [6]Zeng C,Yao Y,Niu Q,Weitering H H 2006Phys.Rev.Lett.96 037204

        [7]Hasan M Z,Kane C L 2010Rev.Mod.Phys.82 3045

        [8]Qi X L,Zhang S C 2011Rev.Mod.Phys.83 1057

        [9]Li B,Wang L,Casati G 2004Phys.Rev.Lett.93 184301

        [10]Chang C W,Okawa D,Majumdar A,Zettl A 2006Science314 1121

        [11]Li B,Wang L,Casati G 2006Appl.Phys.Lett.88 143501

        [12]Rikken G L J A,Strohm C,Wyder P 2002Phys.Rev.Lett.89 133005

        [13]Inyushkin A V,Taldenkov A N 2007JETP Lett.86 379

        [14]Sheng L,Sheng D N,Ting C S 2006Phys.Rev.Lett.96 155901

        [15]Zhang L,Ren J,Wang J S,Li B W 2010Phys.Rev.Lett.105 225901

        [16]Qin T,Zhou J,Shi J 2012Phys.Rev.B86 104305

        [17]Zhang L,Niu Q 2015Phys.Rev.Lett.115 115502

        [18]Zhang L 2016New J.Phys.18 103039

        [19]Liu Y,Xu Y,Duan W 2017 arXiv preprint arXiv:1707.07142

        [20]Zhou J H 2012Ph.D.Dissertation(Beijing:Institute of Theoretical Physics,Chinese Academy of Sciences)(in Chinese)[周建輝 2012博士學(xué)位論文 (北京:中國科學(xué)院理論物理研究所)]

        [21]Nagaosa N,Sinova J,Onoda S,MacDonald A H,Ong N P 2010Rev.Modern Phys.82 1539

        [22]Kagan Y,Maksimov L A 2008Phys.Rev.Lett.100 145902

        [23]Wang J S,Zhang L 2009Phys.Rev.B80 012301

        [24]Zhang L 2011Ph.D.Dissertation(Singapore:National University of Singapore)

        [25]Holz A 1972Nuovo Cimento B9 83

        [26]Strohm C,Rikken G,Wyder P 2005Phys.Rev.Lett.95 155901

        [27]Kronig R L 1939Physica6 33

        [28]van Vleck J H 1940Phys.Rev.57 426

        [29]Wang L and Li B 2007Phys.Rev.Lett.99 177208

        [30]Zhang L,Ren J,Wang J S,Li B W 2011J.Phys.Condens.Matter23 305402

        [31]Zhang L,Wang J S,Li B 2009New J.Phys.11 113038

        [32]Onose Y,Ideue T,Katsura H,Shiomi Y,Nagaosa N,Tokura Y 2010Science329 297

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

        [34]Zeng H L,Cui X D 2016Acta Phys.Sin.45 505

        [35]Chang M C,Niu Q 1996Phys.Rev.B53 7010

        [36]Xiao D,Chang M C,Niu Q 2010Rev.Mod.Phys.82 1959

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

        [38]Gorbachev R V,Song S J C,Yu G L,Kretinin A V,Withers F,Cao Y,Mishchenko A,Grigorieva I V,Novoselov K S,Levitov L S,Geim A K 2014Science346 448

        [39]Cao T,Wang G,Han W P,Ye H Q,Zhu C R,Shi J R,Niu Q,Tan P H,Wang E,Liu B L,Feng J 2012Nat.Commun.3 887

        [40]Hirschberger M,Chisnell R,Young S,Lee N P 2015Phys.Rev.Lett.115 106603

        [41]Heinonen O,Taylor P L,Girvin S M 1984Phys.Rev.B30 3016

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

        [43]Rycerz A,Tworzydlo J,Beenakker C W J 2007Europhys.Lett.79 57003

        [44]Lu J Y,Qiu C Y,Ye L P,Fan X Y,Ke M Z,Zhang F,Liu Z Y 2016Nature Phys.13 369

        [45]Lu J Y,Qiu C Y,Ke M Z,Liu Z Y 2016Phys.Rev.Lett.116 093901

        [46]KaneC L,LubenskyT C 2013arXivpreprint arXiv:1308.0554

        [47]Süsstrunk R,Huber S D 2015Science349 47

        [48]Wang P,Lu L,Bertoldi K 2015Phys.Rev.Lett.115 104302

        PACS:66.70.–f,03.65.Vf,72.10.Bg,72.15.GdDOI:10.7498/aps.66.226601

        *Project supported by the National Natural Science Foundation of China(Grant No.11574154).

        ?Corresponding author.E-mail:phyzlf@njnu.edu.cn

        Topological phonons and phonon Hall e ff ects?

        Xing Yu-Heng Xu Xi-Fang Zhang Li-Fa?

        (School of Physics and Technology,Nanjing Normal University,Nanjing 210023,China)

        28 September 2017;revised manuscript

        6 November 2017)

        The combination of topology and physics is a new field of physics development in recent decades.It is not only active in quantum field theory and high energy physics,but also widely exists in condensed matter physics,including quantum(anomalous,spin)Hall e ff ect and topological insulators(superconductors)etc.Phonon,as the main carrier of heat transport in the crystal,recently,due to the discovery of various phonon devices,phonons has been widely concerned by scientist.In this paper,we introduce the topological properties of phonons and the phonon hall e ff ect.We have reviewed the related physical research progress of phonon hall e ff ect,phonon valley hall e ff ect and so on,which are generated by breaking the time reversal symmetry,spatial inversion symmetry,both breaking the time and spatial inversion symmetry.Finally,the application of topology in other acoustic systems is brie fly introduced,and the future development direction is discussed too.

        topology,Berry phase,Berry curvature,phonon Hall e ff ect

        10.7498/aps.66.226601

        ?國家自然科學(xué)基金(批準(zhǔn)號:11574154)資助的課題.

        ?通信作者.E-mail:phyzlf@njnu.edu.cn

        猜你喜歡
        聲子能帶對稱性
        一類截斷Hankel算子的復(fù)對稱性
        半無限板類聲子晶體帶隙仿真的PWE/NS-FEM方法
        巧用對稱性解題
        吃東西時注意多
        橫向不調(diào)伴TMD患者髁突位置及對稱性
        納米表面聲子 首次實現(xiàn)三維成像
        聲子晶體覆蓋層吸聲機(jī)理研究
        汽車轉(zhuǎn)向管柱吸能帶變形研究和仿真優(yōu)化
        Life OR Death Decision
        基于聲子晶體理論的導(dǎo)線防舞方法及數(shù)值驗證
        青青青伊人色综合久久亚洲综合| 久久久久久久波多野结衣高潮 | 亚洲欧洲日产国码久在线观看| 中文字幕一区二区三区喷水| 真实夫妻露脸爱视频九色网| 国产莉萝无码av在线播放| 五月天久久国产你懂的| 女女同性av一区二区三区免费看| 国产精品国产自产自拍高清av| 乱子伦一区二区三区| 粗大的内捧猛烈进出在线视频| 太大太粗太爽免费视频| 人妻少妇精品视中文字幕免费| 亚洲日韩av一区二区三区中文| 国产免费破外女真实出血视频 | 亚洲人成网7777777国产| 欧美日本视频一区| 香港三级日本三韩级人妇久久| 日本少妇高潮喷水xxxxxxx| 国产精品久久久久久久久鸭| 日韩女优一区二区视频| 中国老熟女露脸老女人| 久久精品黄aa片一区二区三区| 在线精品国内视频秒播| 蜜臀久久久精品国产亚洲av| 亚洲精品乱码久久久久久| 中国老妇女毛茸茸bbwbabes| 最新欧美一级视频| 国产内射一级一片高清内射视频| 久久久久久九九99精品| 午夜片无码区在线| 亚洲av日韩一区二三四五六七| 亚洲精品1区2区在线观看| 天天爽夜夜爽夜夜爽| 2017天天爽夜夜爽精品视频| 蜜桃网站免费在线观看视频| 无码国产69精品久久久久孕妇| 99热这里只有精品4| 久久久精品久久久国产| 国内精品久久久久久99| 国产日韩欧美亚洲精品中字|