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

    多端口石墨烯系統(tǒng)中的非局域電阻?

    2017-11-10 08:25:52王孜博江華謝心澄
    物理學(xué)報 2017年21期
    關(guān)鍵詞:局域霍爾石墨

    王孜博 江華 謝心澄

    1)(中國工程物理研究院微系統(tǒng)與太赫茲研究中心,成都 610200)

    2)(中國工程物理研究院電子工程研究所,綿陽 621999)

    3)(蘇州大學(xué)物理與光電能源學(xué)部,蘇州 215006)

    4)(蘇州大學(xué)高等研究院,蘇州 215006)

    5)(北京大學(xué)物理學(xué)院量子材料科學(xué)中心,北京 100871)

    6)(量子物質(zhì)協(xié)同創(chuàng)新中心,北京 100871)

    多端口石墨烯系統(tǒng)中的非局域電阻?

    王孜博1)2)江華3)4)?謝心澄5)6)

    1)(中國工程物理研究院微系統(tǒng)與太赫茲研究中心,成都 610200)

    2)(中國工程物理研究院電子工程研究所,綿陽 621999)

    3)(蘇州大學(xué)物理與光電能源學(xué)部,蘇州 215006)

    4)(蘇州大學(xué)高等研究院,蘇州 215006)

    5)(北京大學(xué)物理學(xué)院量子材料科學(xué)中心,北京 100871)

    6)(量子物質(zhì)協(xié)同創(chuàng)新中心,北京 100871)

    (2017年9月21日收到;2017年10月7日收到修改稿)

    非局域電阻,石墨烯,自旋霍爾效應(yīng),谷霍爾效應(yīng)

    1 引 言

    非局域測量方法特指在電流流經(jīng)區(qū)域以外對電壓信號的測量,并以非局域電阻的形式體現(xiàn).一種產(chǎn)生非局域電阻的方法是使得電流路徑不同于經(jīng)典的歐姆方式[1?3].例如,在量子霍爾效應(yīng)中,電流可以在樣品表面輸運而體區(qū)域保持絕緣[4].另一種獲得非局域電阻信號的方法是引入除電荷以外的其他自由度,例如自旋.這時的自旋流會由于自旋軌道耦合相互作用而偏離外場所引導(dǎo)的方向.由于非局域電阻通常起源于某些難以直接探測的非平庸物理相互作用,所以非局域測量現(xiàn)在已經(jīng)成為一種用來在新材料中觀測某些電磁學(xué)現(xiàn)象的強有力工具,并在多端口石墨烯體系中得到廣泛的應(yīng)用[5?10].

    目前,已經(jīng)在含如下三種效應(yīng)的多端石墨烯體系中實驗觀測到非局域電阻.

    1)自旋霍爾效應(yīng)(spin Hall effect,SHE).SHE是一種由自旋軌道耦合導(dǎo)致的、使得電流在流過樣品時在其垂直方向產(chǎn)生自旋流的輸運現(xiàn)象[11?17].由于只有自旋而不是電荷在自旋輸運中積累,且石墨烯的自旋軌道耦合強度非常小,以至難以直接在石墨烯中觀測到SHE.借助于非局域測量方法,人們可以在外加垂直電場或磁場的單層石墨烯、氫摻雜的石墨烯等多端H形樣品的狄拉克點或電中性點(charge neutral point,CNP)處觀測到一個巨大的非局域電阻峰值[5,6].在排除了經(jīng)典歐姆輸運等方式后,基本可以確定自旋軌道耦合機理的存在.最近的理論研究還發(fā)現(xiàn),對非局域電阻峰值的貢獻不僅取決于SHE,經(jīng)典的準彈道輸運、準擴散輸運、甚至無序吸附形式等都可能會對非局域電阻的性質(zhì)產(chǎn)生重要的影響,其貢獻可比甚至超過SHE本身[18,19].因此,如何能夠在實驗中進一步壓制這些干擾源而得到純的完全由SHE誘導(dǎo)的非局域電阻非常值得深入思考.

    2)塞曼自旋霍爾效應(yīng)(Zeeman spin Hall effect,ZSHE).ZSHE類似SHE,區(qū)別在于ZSHE通過外磁場打開自旋簡并而不是自旋軌道耦合來產(chǎn)生自旋流[5,20?23].實驗發(fā)現(xiàn),相比于零磁場,外加垂直磁場的石墨烯樣品可以測到一個遠大于前者的非局域電阻[5],這意味著此時的非局域電阻極有可能是由ZSHE誘導(dǎo)的.

    3)谷霍爾效應(yīng)(valley Hall effect,VHE).在打破空間反演對稱性的石墨烯中,谷指標可以起到類似自旋在SHE中的作用[24].最近人們在與襯底公度的石墨烯和外加偏壓的雙層石墨烯樣品的CNP附近同樣觀測到了巨大的非局域電阻[7?9].但是由于此時的石墨烯在CNP處打開能隙,這意味著這些巨大的非局域電阻出現(xiàn)在了能帶的帶隙中,這是極其令人費解的現(xiàn)象[25,26].在排除了SHE、邊緣態(tài)輸運及熱輸運等機理的可能后,人們傾向于認為這個非局域電阻信號可能是由VHE導(dǎo)致的.目前,對于打開能隙的石墨烯中的谷電流探測缺乏直接的實驗證據(jù),因而人們非常希望通過這個可能由VHE誘導(dǎo)的非局域電阻來判斷VHE的存在性.

    綜合來看,對于多端石墨烯樣品中的非局域電阻,目前的研究不僅缺乏對其物理性質(zhì)的深入理解,甚至對其產(chǎn)生機理本身仍存在較大的爭議.本綜述按照SHE,ZSHE和VHE的順序詳細地介紹相關(guān)的多端口石墨烯非局域電阻實驗以及與之對應(yīng)的理論和數(shù)值分析,希望能夠幫助科研工作者了解目前非局域電阻研究的前沿進展,并對未來的非局域電阻研究起到幫助作用.

    2 SHE在石墨烯體系中誘導(dǎo)的非局域電阻

    SHE是一種由自旋軌道耦合導(dǎo)致的、使得自旋非極化的電流在流過樣品時在其垂直方向產(chǎn)生自旋流或者垂直邊界產(chǎn)生自旋積累的輸運現(xiàn)象[11?17].這里的自旋軌道耦合是一種粒子自旋和“磁場”的相互作用,有趣的是,這個“磁場”并非真實存在而是由粒子相對周圍電場的運動而導(dǎo)致.因此,不同于傳統(tǒng)霍爾效應(yīng),實現(xiàn)SHE不需要外加破壞時間反演對稱性的磁場.

    為了定義電荷流通過SHE產(chǎn)生自旋流的強弱,一般使用自旋霍爾角θsH的概念,即產(chǎn)生的自旋流與原電荷流的比值.如圖1(a)的四端系統(tǒng)所示,由端口1入射電荷流大小為I1=I↑+I↓0且自旋流大小為Is1=I↑?I↓=0的電流,SHE導(dǎo)致端口2,3產(chǎn)生電荷流大小為0而自旋流大小為Is2=I↑+I↓0的純自旋流.則自旋霍爾角定義為

    進一步地,(1)式可以利用電導(dǎo)的定義改寫為

    圖1 SHE誘導(dǎo)非局域電阻示意圖(紅色箭頭代表自旋指向) (a)四端體系中,由于SHE,電荷流I1轉(zhuǎn)化為電荷為0的自旋流Is2;(b)六端體系中一個完整的由電荷流經(jīng)SHE和ISHE轉(zhuǎn)化為遠端非局域電荷流的示意圖Fig.1.The schematic diagram for the nonlocal resistance induced by the SHE:(a)In a four-terminal system,the charge current I1can be converted to charge mediated spin current Is2due to the SHE;(b)in a six-terminal system,the charge current is converted to a nonlocal one in a remote region due to the SHE and the ISHE in turn.The arrows colored by red indicate the spin directions.

    其中σxy是自旋霍爾電導(dǎo),σxx是縱向電導(dǎo).由于Landauer-Buttiker公式可以直接給出電導(dǎo)表達式,所以(2)式更利于理論計算.

    類似于SHE,入射的純自旋流還可以通過逆自旋霍爾效應(yīng)(inverse spin Hall effect,ISHE)轉(zhuǎn)化為垂直方向的電荷流[27?29].在一個如圖1(b)所示的六端系統(tǒng)中,端口1的入射電流I1通過SHE首先產(chǎn)生在端口2和5方向輸運的自旋流,這個自旋流隨后通過ISHE可以在端口3和4方向重新產(chǎn)生電荷流.這意味著可以在端口3和4測到一個非局域電壓VNL=V6?V4和非局域電阻RNL=VNL/I1.由于自旋流不帶電,其通常難以直接測量.而通過探測VNL和RNL等非局域信號,可以在實驗上方便地確認SHE的產(chǎn)生.由于純凈的石墨烯材料中自旋軌道耦合強度非常小,通常難以觀測到本征石墨烯中的SHE.目前,非局域電阻實驗通常試圖通過摻雜金、銀、銅等重金屬以增強石墨烯的自旋軌道耦合強度[30].同時,也有文獻表明元素質(zhì)量小的氫摻雜也可以增強石墨烯中的自旋軌道耦合[6].

    2.1 SHE對非局域電阻性質(zhì)的影響

    借助于非局域電阻的測量,實驗上可以在H形石墨烯樣品的狄拉克點處觀測到一個巨大的非局域電阻RNL.Abanin等[31]發(fā)展了一套理論,能夠非常好地解釋這個非局域電阻的起源.更重要的是,他們還給出了一個非局域電阻RNL和局域電阻RL的簡單關(guān)系式:RNL.然而,這個關(guān)系式在解釋某些實驗現(xiàn)象時遇到了一些問題.具體而言,除了RNL在狄拉克點的巨大峰值以外,實驗中還發(fā)現(xiàn)了另一個有意思的現(xiàn)象:對比于RL,RNL在費米能偏離狄拉克點時會更加急劇地衰減到0.例如,在本文的圖13(a)中,紅線所代表的非局域電阻RNL在Vg=1 V時已經(jīng)衰減到幾乎為0,而黑線所代表的局域電阻RL仍在衰減過程中.這個新奇的現(xiàn)象看起來似乎與經(jīng)典的理論公式相矛盾,因為已經(jīng)等于0的RNL不可能正比于不為0的,因而非常值得研究.

    最近,在外加垂直電場引入自旋軌道耦合的情況下,如圖2所示,北京大學(xué)陳劍豪實驗組在H形石墨烯樣品中同樣探測到了一個巨大的非局域電阻信號.類似地,非局域電阻RNL相比RL依然迅速衰減.為了描述這個非局域輸運現(xiàn)象,首先考慮一個含Rashba效應(yīng)的H形四端單層石墨烯系統(tǒng)[18].如圖3所示,電荷流從電極1向電極2縱向入射,在中心區(qū)通過SHE轉(zhuǎn)化為橫向的自旋流,最終通過ISHE轉(zhuǎn)化為縱向電荷流,并作為非局域信號被電極3和4探測到.系統(tǒng)的哈密頓量可以寫為

    圖2 外加垂直電場的H形石墨烯中局域電阻RL和非局域電阻RNL實驗數(shù)據(jù)(Vbg代表樣品背門電壓,Vtg代表樣品頂門電壓,二者的差值引入垂直電場.在實驗中,針對每一個固定的Vtg,Vbg從?60 V掃到60 V.由于費米面EF可以近似地由(Vtg+Vbg)/2決定,圖中Vbg的變化等價于EF的變化.摘自北京大學(xué)陳劍豪教授的未發(fā)表數(shù)據(jù))Fig.2.The local resistance RLand the nonlocal resistance RNLin an H-shaped graphene system with perpendicular electric field.In this fi gure,Vbgand Vtg stand for the bottom gate and the top gate voltage,whose di ff erence induces the electric field.In the experiment,Vbgsweeps from?60 V to 60 V with a fi xed Vtg.Since the fermi energy EFis approximately determined by(Vtg+Vbg)/2,the change of Vtgin this fi gure is equivalent to that of EF.Adapted from the unpublished data from J.H.Chen in Peking University.

    圖3 H形四端石墨烯體系示意圖(計算中電流從電極1流入,電極2流出,黑色虛線的矩形處為中心區(qū))[18]Fig.3.The schematic diagram of the proposed H-shaped four-terminal graphene system[18].The current is injected into lead 1 and flows out from lead 2.The rectangle,denoted by the black dashed line,is the center region.

    其中,和ci分別代表格點i的產(chǎn)生和消滅算符,εi是在位能,λR是外加Rashba耦合強度.為了模擬金屬電極,四個端口的在位能均定義為εUD,這個在位能可以通過門電壓調(diào)節(jié).無序只存在于中心區(qū)并使用均勻分布在[?w/2,w/2]的安德森無序描述,其中w是無序強度.中心區(qū)的示意圖用參數(shù)M,N,P和Q標記.例如,在圖中的系統(tǒng)為M=3,N=2,P=4及Q=3.通過Landauer-Buttiker公式計算,在圖4中得到局域電阻RL=(V1?V2)/I1及非局域電阻RNL=(V3?V4)/I1.

    圖4中的藍線和紅線分別展示了局域電阻RL和非局域電阻RNL隨費米能EF的變化.在圖中,Rashba自旋軌道耦合強度λR從0增加到0.3,而其他參數(shù)保持不變.可以發(fā)現(xiàn)這四張圖主要體現(xiàn)出三個特點.

    首先,非局域電阻在圖4(a)中表現(xiàn)出負值,這里的“負”意味著當(dāng)電流從電極1流向電極2時,電極3測到的電壓竟然比電極4低.類似地,盡管不是嚴格的負數(shù),但也可以在圖4(b)中觀測到相同的趨勢(即一對出現(xiàn)在EF=±0.1附近的谷底).這個現(xiàn)象可以由以往實驗中預(yù)言的準彈道輸運機理來解釋[32].具體而言,從電極1入射的載流子可以直接從電極4射出而不需要返回電極2.也就是說,可以在電極4測到一個正電壓而在電極3測到一個負電壓,這意味著非局域電阻RNL=(V3?V4)/I1此時是一個負值.因此,如圖5所示,可以得到結(jié)論:非局域電阻RNL由以下三部分構(gòu)成

    圖4 局域電阻RL和非局域電阻RNL分別用藍色和紅色在圖中標出(為了增強對比度,非局域電阻RNL放大了7.5倍)[18] (a)λR=0;(b)λR=0.1;(c)λR=0.2;(d)λR=0.3;安德森無序取為w=1Fig.4.The local resistance RLand the nonlocal resistance RNLare drawn in the blue and red lines,respectively[18]:(a)λR=0;(b)λR=0.1;(c)λR=0.2;(d)λR=0.3.The Anderson disorder strength is chosen as w=1.In order to make the comparison clear enough,the value of RNLis ampli fi ed by 7.5 times.

    圖5 三類輸運機理的示意圖(電流從電極1入射,并從電極2出射;紅線代表準彈道輸運機理下的Rballistic,這對于RNL貢獻負值;黃線代表經(jīng)典擴散Rclassic;綠線代表SHE導(dǎo)致的RHall;箭頭代表從電極1入射的電流傳播方向)[18]Fig.5.The schematic diagram for three kinds of transport mechanisms[18].The electron current is injected into lead 1 and then flows out from lead 2.The red line stands for the quasiballistic transport mechanism Rballistic,which makes a negative contribution to RNL.The yellow line represents the classic diffusion Rclassic.The green line denotes the spin Hall transport RHall.Arrows indicate the direction of the electron current which is injected into lead 1.

    圖5中紅色標記的第一項Rballistic為剛討論過的準彈道輸運;棕色標記的第二項Rclassic代表經(jīng)典擴散,也就是歐姆電阻的產(chǎn)生機理;綠色標記的第三項RHall起源于(3)式中的Rashba自旋軌道耦合,電流通過左側(cè)縱向區(qū)域時由于SHE會產(chǎn)生一個與電流方向相垂直的自旋流,并最終在右側(cè)的縱向區(qū)域由于ISHE而重新轉(zhuǎn)化為電流.

    第二,局域和非局域電阻均關(guān)于原點對稱且在EF=0處達到極大值.重要的是,類似于在實驗中觀察到的情況,非局域電阻RNL相比局域電阻RL迅速衰減到0.例如,在圖4(c)中可以發(fā)現(xiàn)當(dāng)|EF|>0.2時,非局域電阻幾乎已經(jīng)完全衰減到0并保持RNL=0.04h/e2,相反局域電阻仍在衰減過程中.這個奇異的現(xiàn)象與實驗結(jié)果一致卻與已知公式[31]

    相矛盾.實際上,這個矛盾似乎可以由負數(shù)部分的非局域電阻來解釋.具體而言,非局域電阻的基礎(chǔ)值是RNL=Rclassic+RHall,其很可能滿足(5)式.然而,根據(jù)(4)式,RNL中還存在一個疊加項Rballistic.由于Rballistic值為負,很自然地RNL會比理論預(yù)言(5)式更迅速地衰減到0.

    第三,如圖4所示,可以觀察到的非局域電阻性質(zhì)和Rashba效應(yīng)間的惟一關(guān)聯(lián)是非局域電阻的負數(shù)值趨勢隨Rashaba耦合強度λR的增加而減弱.相對應(yīng)地,非局域電阻相比局域電阻的迅速衰減現(xiàn)象也逐漸消失.這對于理解非局域電阻與SHE間的內(nèi)在機理仍遠遠不夠.事實上,僅憑圖4,甚至不知道非局域電阻的迅速衰減與SHE間是否存在關(guān)系.因此,除了圖中所計算的RNL,得到一個純的完全由SHE誘導(dǎo)的非局域電阻RHall就顯得尤為重要.

    圖6 新設(shè)計的六端系統(tǒng)示意圖(新增加的電極3′和4′與電極3和4呈鏡面對稱;如果自旋向上(下)的自旋流從電極1沿著紅線輸運,則根據(jù)對稱性,自旋向下(上)的自旋流必定沿黃線輸運)[18]Fig.6.The schematic diagram of the proposed sixterminal system[18].Newly added lead 3′and 4′locate at a mirror symmetry to lead 3 and 4.If the spin-up(down)current injected into lead 1 transports along the red line,then the spin-down(up)current must follow the yellow line.

    圖6中設(shè)計了一個全新的六端系統(tǒng)來代替H形四端系統(tǒng).圖6和圖3的惟一區(qū)別在于電極1和2的左端新添加了兩個電極3′和4′,并使得它們與電極3和4呈鏡面對稱關(guān)系.現(xiàn)在,如果向電極1注入自旋向上極化的自旋流,例如通過引入鐵磁電極的方法,一個非局域電壓信號V34↑可以在電極3和4被探測到.此時的非局域電阻應(yīng)寫為

    類似地,如果注入自旋向下極化的自旋流,可以探測到電壓V34↓,且此時的非局域電阻可以寫為

    由于準彈道輸運和經(jīng)典擴散與自旋極化方向無關(guān),這兩種效應(yīng)誘導(dǎo)的非局域電阻Rballistic和Rclassic在不同自旋方向下必相等,即Rballistic↑=Rballistic↓且Rclassic↑=Rclassic↓. 因此, 在得到RNL↑和RNL↓后,可以排除掉準彈道輸運和經(jīng)典擴散對結(jié)果的影響,并最終得到只由SHE誘導(dǎo)的非局域電阻:

    實際上,圖6中六端系統(tǒng)最大的優(yōu)點在于當(dāng)自旋向上極化的自旋流從電極1入射時,除了電極34,電極3′4′端也可以測得電壓根據(jù)對稱性分析,應(yīng)該等于在電極3,4測到的自旋向下電壓V34↓. 因此,RHall↑?RHall↓可以通過一步測量得到而不需要改變電極1的磁化方向,這對于具體的實驗操作而言至關(guān)重要.但對于理論計算,為了簡單起見,依然選擇使用兩個自旋方向來討論問題.

    圖7 對完全由SHE誘導(dǎo)的非局域電阻RHall的計算[18] (a)自旋向上和自旋向下對應(yīng)的局域電阻RL;(b)自旋向上和自旋向下對應(yīng)的非局域電阻RNL,內(nèi)插圖為自旋霍爾電導(dǎo)σxy;(c)黑線基于RNL↑?RNL↓計算得到,由于數(shù)值誤差在EF=0處引發(fā)了劇烈的擾動,增加了一條紅色虛線來準確地描述RNL↑?RNL↓,考慮紅線的修正后,藍線為黑線的絕對值RNL↑?RNL↓(d)基于(5)式計算的RHallFig.7.The nonlocal resistance RHallwhich is totally induced by the SHE[18].(a)The local resistance RLfor two spin direction.(b)The nonlocal resistance RNLfor two spin directions,and the inset is the spin Hall conductance σxy.(c)The black line is calculated based on RNL↑?RNL↓.Since the numerical error results in dramatic oscillation of RNL↑?RNL↓near EF=0,we add a red dad line to describe the accurate behavior of RNL↑?RNL↓near the Dirac point.At last,the blue line representsRNL↑?RNL↓which equals the absolute value of the black line.And the fi x near the Dirac pint shown by the red line is also considered.(d)RHallcalculated based on Eq.(5).

    圖7(a)首先展示了根據(jù)上述方法計算得到的局域電阻RL.由于局域電阻不依賴于霍爾效應(yīng)而對自旋方向不敏感,所以兩個極化方向完全相反的局域電阻表現(xiàn)出了幾乎完全相同的行為.圖7(b)給出了非局域電阻RNL↑和RNL↓.如圖所示,藍色的RNL↑和紅色的RNL↓反對稱.根據(jù)(5)式中對SHE的解釋,這個反對稱的現(xiàn)象可以由圖6(b)內(nèi)插圖中自旋霍爾電導(dǎo)σxy的行為來理解.具體而言,由于σxy關(guān)于原點反對稱[33]:σxy(EF)=?σxy(?EF),且σxy符號的改變實際上等價于SHE方向的改變.這意味著圖7(b)中藍線和紅線的較小者(即EF<0時的自旋向上和EF>0時的自旋向下)可以被視為沒有SHE的基準線,而較大者即為有SHE.因此,通過在EF<0時從RNL↑中扣除RNL↓,并在EF>0時從RNL↓中扣除RNL↑,可以得到完全由SHE誘導(dǎo)的非局域電阻:

    回到經(jīng)典(5)式,由于它的原始推導(dǎo)實際上只考慮了霍爾效應(yīng),因此(5)式中的RNL理論上應(yīng)該為RHall.在圖7(c)中,首先給出基于RNL↑?RNL↓計算得到的非局域電阻RHall,其中EF=0附近的抖動主要源自數(shù)值誤差.根據(jù)圖7(b),σxy在EF=0處等于0,這意味著系統(tǒng)在狄拉克點處并不存在SHE.因此,RNL在EF=0時對自旋不敏感,即RNL↑必等于RNL↓. 也就是說,RNL↑ ?RNL↓的準確值應(yīng)該由圖7(c)中的紅色虛線表示,它連接兩個峰值且通過原點.為了進一步驗證結(jié)果,根據(jù)(5)式,在圖7(d)中利用局域電阻RL和自旋霍爾電阻σxy重新計算了RHall,并將其與圖7(c)中的RHall進行對比.如圖所示,圖7(c)和圖7(d)的主要行為看起來完全相同,其中的細微區(qū)別可能是由于(5)式中的正比符號含有隨費米面變化的參數(shù).尤為重要的是,可以看到當(dāng)EF等于0時RHall也等于0,這實際上意味著圖4中非局域電阻RNL的尖銳峰與Rashba自旋軌道耦合完全無關(guān).同時,更令人意外的是,這個結(jié)果還意味著非局域電阻RNL相比于局域電阻RL迅速衰減的性質(zhì)同樣不是由Rashba效應(yīng)導(dǎo)致.不僅如此,如圖7(c)和圖7(d)所示,在EF=±0.1附近還存在一對由Rashba自旋軌道耦合導(dǎo)致的RHall峰值.顯然,這對峰會隨著Rashba效應(yīng)的增強而愈發(fā)明顯.現(xiàn)在,對于前文所述的關(guān)于圖4的第三個特征,可以給出一個較完美的解釋.如圖4所示,非局域電阻的負數(shù)值部分同樣出現(xiàn)在EF=±0.1附近,這恰好也是圖7(c)中RHall峰值所在的位置.因此,這對RHall的峰值實際上對非局域電阻RNL的負值趨勢起抵消作用,即表現(xiàn)為RNL相比RL的迅速衰減趨勢隨Rashba效應(yīng)的增強而逐漸消失,也即決定非局域電阻RNL迅速衰減趨勢的實際上是SHE與準彈道輸運的競爭.

    通過計算RHall和進一步的分析可知非局域電阻在狄拉克點表現(xiàn)出的巨大峰值的性質(zhì)與Rashba效應(yīng)無關(guān).實際上,這個極大峰很有可能是由單層石墨烯在狄拉克點處的極小態(tài)密度導(dǎo)致的,這也意味著不能單純地由非局域電阻峰值的出現(xiàn)來判斷SHE的存在性.其次,非局域電阻相比于局域電阻迅速衰減的性質(zhì)并不如之前所預(yù)料的那樣由Rashba效應(yīng)導(dǎo)致.現(xiàn)在,已清楚這個有趣的性質(zhì)實際上是由準彈道輸運誘導(dǎo)的非局域電阻Rballistic所導(dǎo)致.更重要的是,Rashba效應(yīng)本身不僅沒有誘導(dǎo)非局域電阻的迅速衰減,反而對這個迅速衰減的性質(zhì)起到了某種屏蔽作用.

    2.2 對由SHE誘導(dǎo)的非局域電阻組分的進一步分析

    最新的研究細化了由SHE誘導(dǎo)的非局域電阻的組分,并重新驗證了2.1節(jié)的結(jié)論[19].在(3)式的基礎(chǔ)上,考慮內(nèi)稟自旋軌道耦合和無序類型,系統(tǒng)的哈密頓量可以寫為

    其中,第四項代表內(nèi)稟自旋軌道耦合;最后一項是一個隨機在位能,它被用來描述吸附原子帶來的不同無序類型.

    如圖8(a)所示,首先考慮中心區(qū)完全純凈的石墨烯體系.經(jīng)過Landauer-Buttiker公式的計算,圖8(b)首先給出了固定寬度W時,非局域電阻RNL在不同長度L下隨費米面的變化.如圖8(b)所示,當(dāng)W>L時,非局域電阻RNL在狄拉克點表現(xiàn)出了一個巨大的峰值,這種行為類似于某種純凈石墨烯特有的歐姆電導(dǎo)輸運.進一步地,當(dāng)考慮金原子隨機吸附帶來(10)式中不為0的隨機在位能μ0時,如圖9(a)所示,即使此時所有的自旋軌道耦合項仍然被人為地取為0,依然可以觀測到非局域電阻RNL的極大峰值.與2.1節(jié)中的分析相同,這個結(jié)果首先意味著實驗中觀測到的極大峰值實際上與SHE并無關(guān)系,同時暗示構(gòu)成RNL的組合項除了(4)式中提及的三項,還應(yīng)包含在狄拉克點附近產(chǎn)生的“準經(jīng)典擴散”輸運行為.這種輸運機理源于此時的波函數(shù)可以穿透零能隙的狄拉克錐.同時,類似于圖4,這里依然可以觀測到一個負號的RNL.因此,作為對(4)式的細化,總的RNL應(yīng)該考慮以下四項的貢獻:

    其中,Rballistic和RHall意義與(4)式相同.需要強調(diào)的是,Rclassic此時代表由于經(jīng)典擴散造成的平庸歐姆貢獻,而Rqd代表一個石墨烯特有的準擴散輸運.

    此外,如圖9(a)中的兩個插圖所示,在不含自旋軌道耦合且L>W(wǎng)的情況下,RNL表現(xiàn)出負值.這意味著在低吸附原子濃度的情況下,Rclassic可以被近似忽略.具體而言,由于在金原子吸附濃度nad=15%的情況下,自由程?可估計為300—400 nm,所以在?<L的擴散輸運區(qū),歐姆貢獻Rclassic∝exp(?πL/W)由于L/W?1而可以完全忽略.也就是說,在L>W(wǎng)時,非局域電阻的性質(zhì)主要由負號Rballistic和正號RHall間的競爭所決定,這與2.1節(jié)中提出的推測符合.此外,從圖9(a)的插圖還可知,當(dāng)L>W(wǎng)時,非局域電阻RNL沒有表現(xiàn)出極大峰值,即此時的Rqd可以忽略,這也與圖8(b)的結(jié)論一致.圖9(b)還給出了分散分布和團簇分布兩種不同的吸附原子類型所對應(yīng)的自旋霍爾角.顯然,除了自旋軌道耦合項之外,吸附雜質(zhì)的不同類型也會影響SHE的產(chǎn)生.同時,更高的溫度也會降低自旋霍爾角.

    圖8 (a)計算所使用的六端石墨烯示意圖,非局域電阻RNL=VNL/I1,自旋霍爾角θsH=/I1,放大的插圖描述了可能存在的無序雜質(zhì)吸附;(b)當(dāng)不存在無序時,基于純凈石墨烯計算的非局域電阻RNL及在固定寬度W=50 nm的情況下,RNL隨樣品長度L的變化[19]Fig.8.(a)The schematic diagram of a six-terminal graphene employed to compute the nonlocal resistance RNL=RNL/I1and the spin Hall angel θsH=/I1.The enlargement shows the possible adatoms.(b)The nonlocal resistance RNLcalculated in a perfectly clean graphene system,and how RNLvaries with the sample size L with a fi xed width W=50 nm[19].

    圖9 (a)當(dāng)取摻雜濃度nad=15%的均勻摻雜金原子時,固定寬度W=50 nm,非局域電阻RNL隨不同樣品長度的變化(L=10 nm(主圖),L=100 nm(左插圖),L=300 nm(右插圖)),虛線代表(10)式中所有的自旋軌道耦合強度為0;(b)在相同的摻雜濃度下,不同摻雜種類對應(yīng)的不同形態(tài)的自旋霍爾角[19]Fig.9.(a)The nonlocal resistance RNLcalculated based on Fig.8(a)with nad=15%of scattered Au atoms,fi xed channel width W=50 nm,and several channel lengths:L=10 nm(main frame),L=100 nm(left inset),L=300 nm(right inset).Dotted lines plot RNLwhen all SOC terms in Eq.(10)are swithed o ff.(b)The spin Hall angel for the same concentration of Au adatoms which are scattered or clustered[19].

    通過上述分析,Tuan等[19]提出了一種新的排除Rballistic,Rclassic和Rqd,而只保留RHall影響的設(shè)計機理.相比于圖8(a)的示意圖,這個設(shè)計的變化在于:此時連接兩個霍爾條的中心區(qū)不含雜質(zhì)而由完全純凈的石墨烯組成,同時還要保證這個中心區(qū)的輸運通道L足夠長.此時,由于L>W(wǎng),Rqd和Rclassic等于0.同時,由于缺乏傳輸通道中的雜質(zhì)散射,Rballistic也幾乎趨于0.因此,這時圖8(a)中的電中性自旋流完全起源于左側(cè)霍爾條處的SHE,并在右側(cè)霍爾條由ISHE轉(zhuǎn)化為非局域電壓VNL.相比于2.1節(jié)中分離RHall的六端裝置,這個設(shè)計雖然物理正確,但由于其要求中心區(qū)完全純凈,在實際操作中并不現(xiàn)實.需要強調(diào)的是,基于這個設(shè)計得到的RHall同樣在靠近狄拉克點的兩側(cè)表現(xiàn)出一對峰值,而在狄拉克點為0,這與圖7(c)相符合.

    事實上,盡管大量的實驗和理論研究都觀察到了非局域電阻RNL奇怪的負數(shù)性質(zhì),并推測其為準彈道輸運行為的結(jié)果,但對于這個負數(shù)解,尚缺乏足夠的證據(jù)進行有效的分析.例如,這個負數(shù)RNL究竟會在什么情況下出現(xiàn).令人信服的工作需要依賴于如何能從RNL中分離出完全由準彈道輸運誘導(dǎo)的非局域電阻Rballisitc,或者得到非局域電流的微觀圖像.值得注意的是,不同于小尺寸下利用Landauer-Buttiker公式的數(shù)值模擬,當(dāng)尺寸達到擴散極限時,最新的研究通過解析求解自旋流所滿足的連續(xù)性方程,在得到負數(shù)RNL的同時還給出了對這個負數(shù)解的不同理解[34].解析解的優(yōu)勢在于其可以計算出自旋軌道耦合的磁場長度?R和自旋弛豫長度?s.通過對二者進行比較,研究者提出由于自旋進動,RNL的負數(shù)值出現(xiàn)在當(dāng)?R和?s大小可比時.盡管與之前的理論工作存在尺寸上的區(qū)別,但這也許會為人們理解非局域電阻RNL的負數(shù)行為提供一種新的思路.

    3 ZSHE在石墨烯體系中誘導(dǎo)的非局域電阻

    不同于SHE,石墨烯系統(tǒng)中的ZSHE不需要自旋軌道耦合,它是一種在外磁場作用下使得縱向電荷流轉(zhuǎn)化為橫向自旋流的現(xiàn)象[5,20?23].在垂直于面外的磁場作用下,如圖10(a)所示,石墨烯的狄拉克錐由于塞曼相互作用而分裂.這意味著電子型和空穴型的載流子可以在CNP附近獲得相反的自旋.

    圖10 ZSHE的非局域測量示意圖[5] (a)CNP處的塞曼自旋劈裂導(dǎo)致電子和空穴擁有方向相反的自旋;(b)在洛倫茲力的作用下,電子和空穴感受到相反的作用力,電荷流I轉(zhuǎn)化為自旋流I↑和I↓.由于這個大小不為0的自旋流IS=I↑?I↓可以傳播相當(dāng)遠的距離,從而可以在遠端測到一個非局域電壓VNL;(c)非局域電阻RNL的實驗值,在零磁場時,無非局域信號Fig.10.The schematic diagram for the nonlocal experiment of ZSHE[5].(a)Zeeman splitting at the CNP produces two pockets fi lled with electrons and holes having opposite spin.(b)In the presence of the Lorentz force,which has opposite signs for electrons and holes,I gives rise to transverse spin current I↑ and I↓.Since the nonzero spin current IS=I↑?I↓can reach remote regions,a nonlocal voltage VNLcan be detected.(c)The nonlocal conductivity RNLin experiment.There exists no nonlocal signal in zero magnetic field.

    首先在外加垂直磁場的石墨烯系統(tǒng)中實驗觀測到了ZSHE誘導(dǎo)的非局域電阻[5].由于這個非局域信號的樣品長度遠大于傳統(tǒng)歐姆信號的傳播范圍,因而其一定是由一個非平庸的機理所導(dǎo)致.起初,這一類由ZSHE所誘導(dǎo)的非局域信號被認為是由量子霍爾效應(yīng)中的邊緣態(tài)輸運所導(dǎo)致.但實驗發(fā)現(xiàn),即使在室溫(T=300 K)和極弱的外磁場(B?1 T)等明顯超出量子霍爾效應(yīng)適用范圍的情況下,非局域信號依然很明顯.這意味著這種非局域特征是由完全不同于量子霍爾效應(yīng)邊緣態(tài)輸運的機理所導(dǎo)致.此外,由于實驗中使用的六角氮化硼襯底引入了電雜質(zhì)和晶格畸變,這也排除了第2節(jié)中SHE所導(dǎo)致的RNL.因此,這個非局域電阻只能產(chǎn)生于ZSHE.

    類似于SHE中非局域電阻的定義,在圖10(b)中,自旋非極化的電流I從電極6向電極2入射.由于ZSHE,這個電荷流先轉(zhuǎn)化為電中性的自旋流IS=I↑?I↓,然后在逆塞曼自旋霍爾效應(yīng)的作用下被轉(zhuǎn)換為可以在電極3和5探測到的非局域電壓VNL=V5?V3.相對應(yīng)地,非局域電阻的定義為RNL=VNL/I.

    則速度

    且加速度a為

    顯然,能量E+和E?的準粒子在垂直磁場的作用下會感受到兩個方向相反的洛倫茲力,并向相反的方向傳播.不僅如此,盡管塞曼劈裂ΔZ通常很小,但在低溫(kBT<ΔZ)時,ΔZ可以在石墨烯中通過狄拉克錐的自旋退簡并引入兩種不同方向的自旋載流子.如圖10(a)所示,能量E+的準粒子為自旋向上極化,而能量E?的準粒子為自旋向下極化.正是上述兩種效應(yīng)的疊加使得圖10(b)中縱向傳播的電荷流橫向偏轉(zhuǎn)為電中性的自旋流.

    類似于SHE,含ZSHE的多端石墨烯體系同樣可以通過Landauer-Buttiker公式給出一個全數(shù)值模擬[20].在CNP附近,石墨烯中外加磁場的緊束縛哈密頓量可以寫為

    圖11 在含ZSHE的六端系統(tǒng)中計算的電霍爾電阻RH和非局域電阻RNL[20](a)退相干較小時的RH;(b)退相干較小時的RNL;(c)退相干較大時的RH;(d)退相干較大時的RNLFig.11.The charge Hall resistance RHand the nonlocal resistance RNLcalculated based on the six-terminal graphene system with the ZSHE[20]:(a)RHwith a smaller momentum-relaxing dephasing;(b)RNLwith a smaller momentum-relaxing dephasing;(c)RHwith a larger momentum-relaxing dephasing;(d)RNLwith a larger momentum-relaxing dephasing.

    其中,εi是在位能;σ=±1對應(yīng)于自旋向上下,即塞曼劈裂的大小為2gμBσB;t是最近鄰躍遷項;φij是一個由于外磁場導(dǎo)致的躍遷相位.需要注意的是,接下來的計算中,由于樣品放置在量子化外磁場中,磁長度?B和圖10(b)中的輸運通道寬度W的比值需要保持W/?B>1.同時,由于非平衡格林函數(shù)計算的復(fù)雜性要求樣品尺寸不能太大(計算中對于四端體系的寬度W≈2.7 nm,六端體系的寬度W≈2.0 nm),因此需要外加一個極大的磁場以保持W/?B>1.

    圖11給出了在圖10(b)所示的六端體系中使用Landauer-Buttiker公式的非平衡格林函數(shù)計算結(jié)果.在圖11(b)中,非局域電阻RNL=(V5?V3)/I表現(xiàn)出了與實驗圖10(c)非常相似的結(jié)果:狄拉克點處巨大的非局域電阻峰值以及另一對CNP處的側(cè)峰值.圖11(d)相比于圖11(b)增強了相位的退相干,此時的非局域電阻RNL在狄拉克點的峰值降低了兩個數(shù)量級,同時兩個側(cè)峰也完全消失.

    相比于SHE,ZSHE誘導(dǎo)的非局域電阻研究還缺乏足夠深入的分析.例如,由于外磁場的存在,適用于SHE的分析并不完全適用于ZSHE.但ZSHE中的非局域電阻卻有著與SHE類似的形態(tài),那么這兩類非局域電阻是否有著相同的組成,需要進一步的理論和實驗探討.

    4 VHE在石墨烯體系中誘導(dǎo)的非局域電阻

    在打破空間反演對稱性的石墨烯中,谷自由度可以起到與自旋類似的作用[24].這意味著谷自由度也能夠表現(xiàn)出對應(yīng)的霍爾效應(yīng)、磁化、光學(xué)選擇定則甚至手性邊緣態(tài)等現(xiàn)象.最近在MoS2等VI族過渡金屬硫化物中的研究使人們對谷自由度的研究進一步深入[35?37].這類直接半導(dǎo)體材料的帶隙正落在K點,同時其低能電子和空穴可以由帶強自旋谷耦合的狄拉克費米子描述.因此,實驗中往往同時觀測到SHE和VHE.近年來,大量的研究試圖通過電磁學(xué)方法來控制谷自由度,從而尋找可能的谷控制方法[38?40],并催生了谷電子學(xué)這門學(xué)科.

    由于不同谷的電子帶有相反的貝里相位,這使得通過谷產(chǎn)生谷依賴的輸運成為可能.從一個半經(jīng)典的角度考慮,上述輸運現(xiàn)象起源于能帶結(jié)構(gòu)導(dǎo)致的反常群速度[41]:

    其中,k是波矢,?n(k)是第n個能帶的貝里相位,E是外電場.(14)式的最后一項代表了一個同時垂直于外電場和貝里相位的反常群速度.接下來考慮外加一個打破空間反演對稱性的微擾,石墨烯打開能隙Δ,這時的狄拉克點附近的低能哈密頓量寫為

    則兩個谷附近的貝里相位形式為

    顯然這里不同的谷分別對應(yīng)兩個相反的貝里相位.結(jié)合(14)式,此時兩個谷感受到兩個完全相反的類洛倫茲力,因而會向兩個相反的方向傳播,如圖12所示.

    圖12 由于貝里相位的作用,屬于不同谷的電子會感受到兩個完全相反的類洛倫茲力Fig.12.Due to the effect of Berry curvatures,the electrons in two different valleys will experience the Lorentz force in opposite directions.

    從量子力學(xué)的角度而言,在K點附近,(15)式中的哈密頓量可以寫為

    則此時的速度和加速度算符為

    最近的實驗首先在與六角硼化氮襯底公度的石墨烯樣品中觀測到了非局域電阻信號[7].如圖13(a)和圖13(b)所示,當(dāng)石墨烯和硼化氮公度放置時,石墨烯的空間反演對稱性被打破并在CNP打開能隙,此時樣品表現(xiàn)出一個巨大的非局域電阻值RNL;而當(dāng)石墨烯和硼化氮交錯一個角度時,石墨烯保持空間反演對稱性且不打開能隙,此時的非局域電阻值非常小,以至于幾乎觀測不到,這排除了源于SHE的自旋流可能對圖13(a)中非局域電阻的貢獻.同時,實驗還認真地排除了邊緣態(tài)等可能對非局域電阻的貢獻.因此,一個合理的推測是這時的非局域信號可能源于由VHE誘導(dǎo)的谷電流.

    圖13 (a)在與六角硼化氮襯底公度的石墨烯樣品中探測非局域電阻的示意圖(插圖)及非局域電阻RNL實驗數(shù)據(jù)(主圖);(b)實驗中在公度與非公度情況下分別測量的非局域電阻RNL值[7]Fig.13.(a)The schematic diagram of the experimental device(inset)and the measured nonlocal resistance(main panel)for a multi-terminal graphene system aligned over an h-BN substrate;(b)experimental nonlocal resistance for aligned and non-aligned h-BN substrates[7].

    更多的實驗在外加垂直電場的雙層石墨烯樣品中同樣觀察到了巨大的非局域電阻[8,9].類似地,由于雙層石墨烯此時也在CNP打開能隙,這個非局域電阻信號同樣被傾向于理解為是VHE的結(jié)果.實際上,這類實驗在理論上最難以理解的地方在于:所有的非局域電阻信號均出現(xiàn)在當(dāng)費米能落在帶隙中時,而這里本不應(yīng)該存在任何物理的結(jié)果.Lensky等[25]認為這里的VHE可能起源于能隙以下的費米海體態(tài).對于在帶隙中的費米能EF,兩個相反的谷極化流在系統(tǒng)中傳播并產(chǎn)生電中性的谷電流.由于谷電流由能隙以下費米海中的電子傳播,谷電流是非耗散的.因此,在這個圖像中,即使系統(tǒng)表現(xiàn)為電絕緣且化學(xué)勢落在帶隙中,由費米海導(dǎo)致的谷電流仍然可以產(chǎn)生一個看起來違反直覺的非局域信號.針對這個觀點,進一步的研究直接使用了參照圖13中打開能隙且谷霍爾電導(dǎo)0的石墨烯哈密頓量,非平衡格林函數(shù)的計算結(jié)果表明非局域電阻RNL在帶隙內(nèi)永遠為零.類似地,基于文獻[8,9]中實驗使用的雙層哈密頓量的計算給出了同樣的結(jié)果.最近,借助于第一性原理方法對哈密頓量的修正,Marmolejo-Tejada等[26]提出,這個非局域電阻應(yīng)該是由某種特殊的出現(xiàn)在狄拉克點附近的邊緣態(tài)所誘導(dǎo),而并不是費米海體態(tài).

    圖14 計算所使用的含子格子錯位勢Δ的四端石墨烯示意圖Fig.14.The schematic diagram for the calculated four-terminal graphene system with a staggered sublattice potential Δ.

    根據(jù)Landauer-Buttiker公式,Cresti等[42]在一個如圖14所示的四端石墨烯中模擬了由VHE誘導(dǎo)的非局域電阻.此時,一個打破空間反演對稱性的哈密頓量可以寫為

    這里,ai和bi對應(yīng)于兩種子格子,并將寬度固定為W=50 nm.首先,當(dāng)錯位勢Δ取為0時,石墨烯不打開能隙.相對應(yīng)的計算結(jié)果顯示,無論長度L如何取值,系統(tǒng)都不表現(xiàn)出可以被觀測到的非局域電阻RNL.接下來,當(dāng)錯位勢取為Δ=60.2 meV時,類似于實驗現(xiàn)象,此時的石墨烯打開能隙Δ并在狄拉克點表現(xiàn)出一個巨大的非局域電阻峰值.顯然,這個對比結(jié)果意味著基于(19)式得到的非局域電阻實際上源于VHE的貢獻.更細致的計算表明,非局域電阻的峰值會隨著長度L的減小而降低.這是由于能隙中的輸運一般由能隙外的體態(tài)向能隙內(nèi)的隧穿來實現(xiàn),而隧穿概率會隨著樣品尺寸的增大而強烈地衰減.

    綜上所示,盡管目前并沒有對谷電流直接測量的實驗證據(jù),但理論和實驗研究都傾向于認為在打破空間反演對稱性的石墨烯中觀測到的巨大的非局域電阻值來源于VHE.更令人信服的結(jié)論還需要進一步深入的理論和數(shù)值計算.

    5 結(jié) 論

    本綜述系統(tǒng)地回顧了實驗上觀測到的可能分別由SHE,ZSHE和VHE誘導(dǎo)的非局域電阻,并介紹了通過Landauer-Buttiker公式給出的上述三種霍爾效應(yīng)在多端系統(tǒng)中誘導(dǎo)非局域電阻的全數(shù)值結(jié)果.對于SHE,理論計算較好地給出了符合實驗現(xiàn)象的非局域電阻模擬值.同時,更細致的理論研究預(yù)言實驗中觀測到的非局域電阻峰值的產(chǎn)生很可能與SHE無關(guān),而非局域電阻相比局域電阻的迅速衰減很可能源于SHE與準彈道輸運的競爭關(guān)系.通過構(gòu)造一個可以產(chǎn)生ZSHE和VHE的哈密頓量,也可以得到與實驗結(jié)果符合較好的非局域電阻計算結(jié)果.但相比于SHE,對于ZSHE和VHE如何誘導(dǎo)非局域電阻的物理機理目前缺乏足夠令人信服的結(jié)論.特別是對于VHE,非局域電阻出現(xiàn)在本不應(yīng)該測到任何物理信號的體能隙中,這一直是一個令人費解的現(xiàn)象,并引發(fā)了激烈的爭論.總體而言,盡管非局域電阻在近年來已經(jīng)吸引了人們足夠高的研究興趣,但由于其非平庸的物理性質(zhì),對于非局域電阻仍需要更深入的理論研究和更直接的實驗測量結(jié)果.

    感謝中國北京大學(xué)陳劍豪教授、韓偉教授,美國特拉華大學(xué)Branislav K.Nikolic教授和中國臺灣新竹清華大學(xué)Chunli Huang博士的討論.

    [1]Brune C,Buhmann H,Molenkamp L W,Maciejko J,Qi X L,Zhang S C 2009Science325 294

    [2]Chang C Z,Zhao W,Kim K Y,Wei P,Jain J K,Liu C,Chan M H W,Moodera J S 2015Phys.Rev.Lett.115 057206

    [3]Parameswaran S A,Grover T,Abanin D A,Pesin D A,Vishwanath A 2014Phys.Rev.X4 031035

    [4]McEuen P L,Szafer A,Richter C A,Alphenaar B W,Jain J K,Stone A D,Wheeler R G,Sacks R N 1990Phys.Rev.Lett.64 2062

    [5]Abanin D A,Morozov S V,Ponomarenko L A,Gorbachev R V,Mayorov A S,Katsnelson M I,Watanabe K,Taniguchi T,Novoselov K S,Levitov L S,Geim A K 2011Science332 328

    [6]Balakroshnan J,Koon G K W,Jaiswal M,Castro Neto A H,Ozyilmaz B 2013Nat.Phys.9 284

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

    [8]Shimazaki Y,Yamamoto M,Borzenets I V,Watanabe K,Taniguchi T,Tarucha S 2015Nat.Phys.11 1032

    [9]Sui M,Chen G,Ma L,Shan W Y,Tian D,Watanabe K,Taniguchi T,Jin X,Yao W,Xiao D,Zhang Y 2015Nat.Phys.11 1027

    [10]Yamamoto M,Shimazaki Y,Borzenets I V,Tarucha S 2015J.Phys.Soc.Jpn.84 121006

    [11]Hirsch J E 1999Phys.Rev.Lett.83 1834

    [12]Murakami S,Nagaosa N,Zhang S C 2003Science301 1348

    [13]Sinova J,Culcer D,Niu Q,Sinitsyn N A,Jungwirth T,MacDonald A H 2004Phys.Rev.Lett.92 126603

    [14]Kato Y K,Myers R C,Gossard A C,Awsschalom D D 2004Science306 1910

    [15]Kimura T,Otani Y 2007Phys.Rev.Lett.99 196604

    [16]Brune C,Roth A,Novik E G,Konig M,Buhmann H,Hankiewicz E M,Hanke W,Sinova J,Molenkamp L W 2010Nat.Phys.6 448

    [17]Sheng L,Sheng D N,Ting C S,Haldane F D M 2005Phys.Rev.Lett.95 136602

    [18]Wang Z,Liu H,Jiang H,Xie X C 2016Phys.Rev.B94 035409

    [19]Tuan D V,Marmolejo-Tejada J M,Waiintal X,Nikolic B K,Valenzuela S O,Roche S 2016Phys.Rev.Lett.117 176602

    [20]Chen C L,Chang C R,Nikolic B K 2012Phys.Rev.B85 155414

    [21]Renard J,Studer M,Folk J A 2014Phys.Rev.Lett.112 116601

    [22]Wei P,Lee S,Lemaitre F,Pinel L,Cutaia D,Cha W,Katmis F,Zhu Y,Heiman D,Hone J,Moodera J S,Chen C T 2016Nat.Mater.15 711

    [23]Abanin D A,Gorbachev R V,Novoselov K S,Geim A K,Levitov L S 2011Phys.Rev.Lett.107 096601

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

    [25]Lensky Y D,Song J C W,Samutpraphoot P,Levitov L S 2015Phys.Rev.Lett.114 256601

    [26]Marmolejo-Tejada J M,Garcia J H,Chang P H,Sheng X L,Cresti A,Roche S,Nikolic B K 2017 arXiv:1706.09361v1[cond-mat.mes-hall]

    [27]Saitoh E,Ueda M,Miyajima H,Tatara G 2006Appl.Phys.Lett.88 182509

    [28]Valenzuela S O,Tinkham M 2006Nature442 176

    [29]Kimura T,Otani Y,Sato T,Takahashi S,Maekawa S 2007Phys.Rev.Lett.98 156601

    [30]Balakrishnan J,Koon G K W,Avsar A,Ho Y,Lee J H,Jaiswal M,Baeck S,Ahn J,Ferreira A,Cazalilla M A,Castro Neto A H,Ozyilmaz B 2014Nat.Commun.5 4748

    [31]Abanin D A,Shytov A V,Levitov L S,Halperin B I 2009Phys.Rev.B79 035304

    [32]Mihajlovic G,Pearson J E,Garcia M A,Bader S D,Ho ff mann A 2009Phys.Rev.Lett.103 166601

    [33]Sheng L,Sheng D N,Ting C S 2005Phys.Rev.Lett.94 016602

    [34]Huang C,Chong Y D,Cazalilla M A 2017Phys.Rev.Lett.119 136804

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

    [36]Lee J,Mak K F,Shan J 2016Nat.Nanotechnol.11 421

    [37]Xiao D,Liu G B,Feng W X,Xu X D,Yao W 2012Phys.Rev.Lett.108 196802

    [38]Zhai F,Zhao X,Chang K,Xu H Q 2010Phys.Rev.B82 115442

    [39]Gunlycke D,White C T 2011Phys.Rev.Lett.106 136806

    [40]Zhai F,Ma Y L,Chang K 2011New J.Phys.13 083029

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

    [42]Cresti A,Nikolic B K,Garcia J H,Roche S 2016La Rivista Del Nuovo Cimento12 587

    PACS:72.80.Vp,71.70.Ej,85.35.–p DOI:10.7498/aps.66.217201

    *Project supported by the National Basic Research Program of China(Grant Nos.2015CB921102,2014CB920901),the National Natural Science Foundation of China(Grant Nos.11704348,11374219,11534001,11404300),and the Science Challenge Project,China(Grant No.TZ2016003-1).

    ?Corresponding author.E-mail:jianghuaphy@suda.edu.cnc

    Nonlocal resistance in multi-terminal graphene system?

    Wang Zi-Bo1)2)Jiang Hua3)4)?Xie Xin-Cheng5)6)

    1)(Microsystems and Terahertz Research Center,China Academy of Engineering Physics,Chengdu 610200,China)
    2)(Institute of Electronic Engineering,China Academy of Engineering Physics,Mianyang 621999,China)
    3)(College of Physics,Optoelectronics and Energy,Soochow University,Suzhou 215006,China)
    4)(Institute for Advanced Study of Soochow University,Suzhou 215006,China)
    5)(International Center for Quantum Materials,School of Physics,Peking University,Beijing 100871,China)
    6)(Collaborative Innovation Center of Quantum Matter,Beijing 100871,China)

    d 21 September 2017;revised manuscript

    7 October 2017)

    Since the nonlocal measurement is helpful in discovering nontrivial physics that is too difficult to detect directly,the nonlocal measurement has now become one of the research focuses in condensed matter physics.Recent experiments find the signal of the giant nonlocal resistance in an H-shaped multi-terminal graphene system.After excluding other possible transport mechanisms,such as the classic Ohmic diffusion and the edge states,researchers tend to believe that the nonlocal resistance signal originates from the spin/valley Hall effect existing in graphene sample.Based on the Landauer-Buttiker formula,the numerical results make a relatively perfect match with the experimental data in the same multi-terminal graphene system.However,though the theoretic research has made certain progress in explaining the existence of the nonlocal resistance,it is still difficult to understand some exotic behaviors of the nonlocal resistance,which exhibits properties even contradictory to the known classical theories.For instance,the nonlocal resistance decreases to zero much more rapidly than the local one,and the giant peak of the nonlocal resistance appears inside the energy gap of the graphene.In this review,the experiments focusing on the nonlocal resistance in multi-terminal graphene system are carefully reviewed.Besides,this review also shows the associated theoretic studies,and an overlook of the future study is also provided.

    nonlocal resistance,graphene,spin Hall effect,valley Hall effect

    非局域測量方法由于其能夠間接探測某些難以直接俘獲的非平庸物理機理,近年來已逐漸成為凝聚態(tài)物理的研究熱點之一.最近的實驗在H形多端口石墨烯樣品中發(fā)現(xiàn)了巨大的非局域電阻信號.在排除了經(jīng)典歐姆、邊緣態(tài)等可能的輸運形式后,人們傾向于認為這類非局域電阻是由多端石墨烯系統(tǒng)中存在的自旋霍爾效益或谷霍爾效應(yīng)所導(dǎo)致.借助于非平衡格林函數(shù)輸運計算,目前的理論可以在同樣的多端石墨烯體系中得到部分與實驗符合較好的數(shù)值模擬結(jié)果.針對實驗中發(fā)現(xiàn)的某些難以理解的、甚至與經(jīng)典理論相矛盾的非局域電阻性質(zhì),例如非局域電阻相比局域電阻在偏離電中性點時的迅速衰減、出現(xiàn)在能隙中的非局域電阻峰值等,目前的理論研究取得了一定的進展,但對這些奇異現(xiàn)象的理解仍存在較大的爭議.本綜述詳細回顧了多端口石墨烯體系中非局域電阻的相關(guān)實驗,并針對性地介紹與之配套的理論進展及對未來研究的展望.

    10.7498/aps.66.217201

    ?國家重點基礎(chǔ)研究發(fā)展計劃(批準號:2015CB921102,2014CB920901)、國家自然科學(xué)基金(批準號:11704348,11374219,11534001,11404300)和科學(xué)挑戰(zhàn)計劃(批準號:TZ2016003-1)資助的課題.

    ?通信作者.E-mail:jianghuaphy@suda.edu.cnc

    ?2017中國物理學(xué)會Chinese Physical Society

    猜你喜歡
    局域霍爾石墨
    石墨系升溫球的實踐與應(yīng)用
    昆鋼科技(2022年1期)2022-04-19 11:36:14
    局域積分散列最近鄰查找算法
    電子測試(2018年18期)2018-11-14 02:30:34
    石墨烯的健康路
    這個大童話講貓(五)名偵探西尼·霍爾
    PET成像的高分辨率快速局域重建算法的建立
    離子推力器和霍爾推力器的異同
    太空探索(2015年6期)2015-07-12 12:48:42
    基于局域波法和LSSVM的短期負荷預(yù)測
    電測與儀表(2015年7期)2015-04-09 11:39:50
    基于非正交變換的局域波束空時自適應(yīng)處理
    石墨礦中固定碳的分析與探討
    道格拉斯·斯高特·霍爾(1940-2013)
    91在线精品国自产拍蜜月| 一级黄片播放器| 精品人妻熟女av久视频| 久久久久久久亚洲中文字幕| 亚洲性久久影院| 伦精品一区二区三区| 免费看光身美女| 久久久午夜欧美精品| 欧美精品国产亚洲| 国产精品人妻久久久影院| 亚洲高清免费不卡视频| 性欧美人与动物交配| 亚洲欧美日韩高清在线视频| 国产精品野战在线观看| 看免费成人av毛片| 欧美bdsm另类| 免费观看人在逋| 国产成人a∨麻豆精品| 女人被狂操c到高潮| 国产精品电影一区二区三区| 看片在线看免费视频| 国产精品电影一区二区三区| 久久精品91蜜桃| 波野结衣二区三区在线| 一级a爱片免费观看的视频| 国国产精品蜜臀av免费| 亚洲国产色片| 国产av麻豆久久久久久久| 欧美三级亚洲精品| 日韩欧美三级三区| 欧美xxxx黑人xx丫x性爽| 欧美成人免费av一区二区三区| 久久久久久久久中文| 免费av观看视频| 国产精品综合久久久久久久免费| 亚洲国产欧洲综合997久久,| 在线免费十八禁| 日本一本二区三区精品| 成人特级黄色片久久久久久久| 午夜a级毛片| 黄色视频,在线免费观看| 久久久国产成人免费| 国产在视频线在精品| 亚洲色图av天堂| 婷婷精品国产亚洲av| av国产免费在线观看| 黄色视频,在线免费观看| 午夜爱爱视频在线播放| 国产精品电影一区二区三区| 天天躁夜夜躁狠狠久久av| 免费高清视频大片| 麻豆成人午夜福利视频| 欧美高清成人免费视频www| 日韩人妻高清精品专区| 久久久久久大精品| 韩国av在线不卡| 小蜜桃在线观看免费完整版高清| 三级国产精品欧美在线观看| 久久久精品94久久精品| 欧美色欧美亚洲另类二区| 亚洲中文日韩欧美视频| 午夜福利18| 国产精品精品国产色婷婷| 非洲黑人性xxxx精品又粗又长| 啦啦啦观看免费观看视频高清| 女同久久另类99精品国产91| 99久久中文字幕三级久久日本| 两个人视频免费观看高清| 美女高潮的动态| 欧美日本视频| eeuss影院久久| 一区二区三区四区激情视频 | 如何舔出高潮| 婷婷亚洲欧美| 亚州av有码| 亚洲四区av| 久久久久久久亚洲中文字幕| 97碰自拍视频| 国产精品精品国产色婷婷| 成年免费大片在线观看| 亚洲色图av天堂| 婷婷色综合大香蕉| 欧美成人一区二区免费高清观看| 亚洲图色成人| 1000部很黄的大片| 18禁在线播放成人免费| 精品久久久久久成人av| 伊人久久精品亚洲午夜| 国内揄拍国产精品人妻在线| 99久国产av精品国产电影| 国内精品久久久久精免费| 插逼视频在线观看| 搡女人真爽免费视频火全软件 | 看片在线看免费视频| 一进一出抽搐gif免费好疼| 女生性感内裤真人,穿戴方法视频| 一本一本综合久久| 亚洲欧美成人精品一区二区| 无遮挡黄片免费观看| 日韩成人av中文字幕在线观看 | 国产高清激情床上av| 欧美在线一区亚洲| 成人美女网站在线观看视频| 人人妻,人人澡人人爽秒播| 九九热线精品视视频播放| 免费在线观看成人毛片| 日韩精品中文字幕看吧| 日日撸夜夜添| 非洲黑人性xxxx精品又粗又长| 国产精品人妻久久久久久| 自拍偷自拍亚洲精品老妇| 少妇裸体淫交视频免费看高清| 夜夜夜夜夜久久久久| 日日摸夜夜添夜夜添av毛片| 热99在线观看视频| 在线免费十八禁| 看黄色毛片网站| 中文字幕av成人在线电影| 露出奶头的视频| 一级黄片播放器| 免费看光身美女| 九九久久精品国产亚洲av麻豆| 夜夜爽天天搞| 国产 一区 欧美 日韩| 身体一侧抽搐| 亚洲精品粉嫩美女一区| 夜夜夜夜夜久久久久| 欧美激情国产日韩精品一区| 免费电影在线观看免费观看| 在线免费观看不下载黄p国产| 赤兔流量卡办理| 免费看a级黄色片| 观看美女的网站| 午夜福利成人在线免费观看| 永久网站在线| 综合色丁香网| 十八禁国产超污无遮挡网站| 亚洲真实伦在线观看| 91在线精品国自产拍蜜月| 亚洲国产精品国产精品| 男女啪啪激烈高潮av片| 日韩av在线大香蕉| 天堂网av新在线| 日韩一本色道免费dvd| 99九九线精品视频在线观看视频| 特大巨黑吊av在线直播| 亚洲成人久久性| 国产一区二区激情短视频| 在线免费观看不下载黄p国产| 欧美日本视频| 美女黄网站色视频| 美女高潮的动态| 网址你懂的国产日韩在线| 99热只有精品国产| 久久午夜亚洲精品久久| 久久亚洲国产成人精品v| 欧美bdsm另类| 午夜福利视频1000在线观看| 69av精品久久久久久| 亚洲国产精品久久男人天堂| 国产乱人视频| 一区二区三区四区激情视频 | 欧美激情久久久久久爽电影| 九色成人免费人妻av| a级毛色黄片| 国内精品一区二区在线观看| 日韩三级伦理在线观看| 中文字幕精品亚洲无线码一区| 日日摸夜夜添夜夜添av毛片| 我的女老师完整版在线观看| 女人被狂操c到高潮| 精品福利观看| 久久久久久国产a免费观看| 国产单亲对白刺激| 三级经典国产精品| 丰满乱子伦码专区| 天堂影院成人在线观看| 12—13女人毛片做爰片一| 91久久精品国产一区二区三区| 成人鲁丝片一二三区免费| 波野结衣二区三区在线| 变态另类丝袜制服| 高清毛片免费看| 两个人的视频大全免费| 一个人观看的视频www高清免费观看| 在线播放国产精品三级| 又粗又爽又猛毛片免费看| 亚洲av第一区精品v没综合| 日本三级黄在线观看| 国产精品久久久久久久电影| 亚洲av一区综合| 国内久久婷婷六月综合欲色啪| 午夜福利在线观看免费完整高清在 | 美女黄网站色视频| 韩国av在线不卡| 波野结衣二区三区在线| 久久久精品94久久精品| 又黄又爽又刺激的免费视频.| ponron亚洲| 免费观看人在逋| 最近最新中文字幕大全电影3| 久久热精品热| 一边摸一边抽搐一进一小说| 亚洲成av人片在线播放无| 国内揄拍国产精品人妻在线| 麻豆精品久久久久久蜜桃| 国语自产精品视频在线第100页| 97超视频在线观看视频| 成人亚洲精品av一区二区| 夜夜看夜夜爽夜夜摸| 中文字幕熟女人妻在线| 成人美女网站在线观看视频| 久久亚洲精品不卡| 国产69精品久久久久777片| 成人亚洲欧美一区二区av| or卡值多少钱| 乱码一卡2卡4卡精品| 婷婷精品国产亚洲av在线| 久久6这里有精品| 午夜精品一区二区三区免费看| 男人舔奶头视频| 99久久无色码亚洲精品果冻| 可以在线观看毛片的网站| 深夜a级毛片| 亚洲精品国产av成人精品 | 少妇裸体淫交视频免费看高清| 又爽又黄a免费视频| 蜜臀久久99精品久久宅男| 一区二区三区免费毛片| 亚洲欧美成人精品一区二区| 亚洲av免费在线观看| 性欧美人与动物交配| 国产亚洲精品av在线| 自拍偷自拍亚洲精品老妇| 最近手机中文字幕大全| 看黄色毛片网站| 国产爱豆传媒在线观看| 国产精华一区二区三区| 国产精品不卡视频一区二区| 免费不卡的大黄色大毛片视频在线观看 | 日韩,欧美,国产一区二区三区 | 亚洲美女黄片视频| 国产精品嫩草影院av在线观看| 免费av观看视频| 国产精品一及| 内地一区二区视频在线| 日韩欧美精品v在线| 99久国产av精品国产电影| 亚洲国产精品合色在线| 国产av不卡久久| 亚洲av电影不卡..在线观看| 亚洲美女黄片视频| 成人av在线播放网站| 看十八女毛片水多多多| 免费看美女性在线毛片视频| 国产伦精品一区二区三区四那| 欧美日韩精品成人综合77777| 亚洲欧美日韩高清专用| 一级黄色大片毛片| 日日摸夜夜添夜夜爱| 小蜜桃在线观看免费完整版高清| 国产精品国产三级国产av玫瑰| 久久99热6这里只有精品| 免费观看在线日韩| 成人无遮挡网站| 亚洲无线观看免费| 一级毛片电影观看 | 女人十人毛片免费观看3o分钟| 成人美女网站在线观看视频| 香蕉av资源在线| 国产精品伦人一区二区| 亚洲精品在线观看二区| 好男人在线观看高清免费视频| 国产探花极品一区二区| 激情 狠狠 欧美| 久久久a久久爽久久v久久| 一级a爱片免费观看的视频| 最近在线观看免费完整版| 免费看光身美女| 九九热线精品视视频播放| 国产大屁股一区二区在线视频| 最好的美女福利视频网| 最近在线观看免费完整版| 白带黄色成豆腐渣| 狠狠狠狠99中文字幕| 国产探花在线观看一区二区| 99热只有精品国产| 日本三级黄在线观看| 午夜福利在线在线| 国产真实乱freesex| 三级男女做爰猛烈吃奶摸视频| 亚洲最大成人av| 精品久久久久久久久久免费视频| 村上凉子中文字幕在线| 99久国产av精品国产电影| 两性午夜刺激爽爽歪歪视频在线观看| 九色成人免费人妻av| 一本精品99久久精品77| 亚洲综合色惰| 尤物成人国产欧美一区二区三区| 美女大奶头视频| 99在线视频只有这里精品首页| 日韩精品中文字幕看吧| 免费大片18禁| 色5月婷婷丁香| 亚洲国产欧美人成| 国产色婷婷99| 国产精品,欧美在线| 国产精品久久电影中文字幕| 久久人人精品亚洲av| 久久久色成人| 伦理电影大哥的女人| 中国国产av一级| 国产精品野战在线观看| 亚洲国产欧美人成| 日日摸夜夜添夜夜添小说| 精品一区二区免费观看| 蜜桃久久精品国产亚洲av| 日韩欧美精品v在线| 秋霞在线观看毛片| 嫩草影视91久久| 久久久色成人| 1000部很黄的大片| 别揉我奶头~嗯~啊~动态视频| 看十八女毛片水多多多| 国产精品久久久久久av不卡| 又黄又爽又免费观看的视频| 我要搜黄色片| 久久精品国产亚洲av涩爱 | 在线免费十八禁| 中文亚洲av片在线观看爽| 亚洲av熟女| 亚洲精华国产精华液的使用体验 | 国内精品宾馆在线| 不卡一级毛片| 国产成人影院久久av| .国产精品久久| 久久6这里有精品| 不卡视频在线观看欧美| 一个人免费在线观看电影| 久久亚洲国产成人精品v| 丰满乱子伦码专区| 97超碰精品成人国产| 久久人人爽人人爽人人片va| 国内精品久久久久精免费| 97超视频在线观看视频| 精品久久久久久久久久免费视频| 久久中文看片网| 亚洲人成网站在线播| 日本色播在线视频| 老司机午夜福利在线观看视频| 欧美三级亚洲精品| 精品欧美国产一区二区三| 国产熟女欧美一区二区| 亚洲精品国产av成人精品 | 国产亚洲精品久久久com| 99视频精品全部免费 在线| 18禁在线无遮挡免费观看视频 | 给我免费播放毛片高清在线观看| 插逼视频在线观看| 久久九九热精品免费| 亚洲电影在线观看av| 春色校园在线视频观看| 寂寞人妻少妇视频99o| 美女cb高潮喷水在线观看| 九九在线视频观看精品| 亚洲欧美清纯卡通| 亚洲欧美精品自产自拍| 国内精品宾馆在线| 国产精品综合久久久久久久免费| 别揉我奶头~嗯~啊~动态视频| 亚洲中文字幕日韩| 亚洲四区av| 身体一侧抽搐| 精品一区二区三区视频在线| 国内少妇人妻偷人精品xxx网站| 春色校园在线视频观看| 国产免费一级a男人的天堂| 国产欧美日韩一区二区精品| 国产免费一级a男人的天堂| 国产精品嫩草影院av在线观看| 国产熟女欧美一区二区| 欧美色视频一区免费| 九九在线视频观看精品| 日韩亚洲欧美综合| 成人毛片a级毛片在线播放| 好男人在线观看高清免费视频| 全区人妻精品视频| 久久久色成人| 嫩草影视91久久| 久久人人爽人人片av| 99国产极品粉嫩在线观看| 久久精品夜色国产| 久久人人精品亚洲av| 韩国av在线不卡| 午夜影院日韩av| 久久人人爽人人片av| av中文乱码字幕在线| 久久精品夜色国产| 国产精品无大码| 欧美成人精品欧美一级黄| 国产不卡一卡二| 极品教师在线视频| 男女做爰动态图高潮gif福利片| 99热这里只有精品一区| 一进一出好大好爽视频| 国产色爽女视频免费观看| 日本精品一区二区三区蜜桃| 99热这里只有是精品50| 国产成人福利小说| 日日摸夜夜添夜夜添av毛片| 两个人视频免费观看高清| 伦理电影大哥的女人| 一区二区三区免费毛片| 国内精品美女久久久久久| av女优亚洲男人天堂| 国产精品综合久久久久久久免费| 亚洲人与动物交配视频| 国产午夜精品久久久久久一区二区三区 | 真实男女啪啪啪动态图| 国产男人的电影天堂91| 高清毛片免费观看视频网站| 日韩,欧美,国产一区二区三区 | 国产精品野战在线观看| 舔av片在线| 国产精品99久久久久久久久| 激情 狠狠 欧美| 内射极品少妇av片p| 久久6这里有精品| 国产黄片美女视频| 久久久色成人| 人人妻,人人澡人人爽秒播| 中文字幕精品亚洲无线码一区| 此物有八面人人有两片| 欧美中文日本在线观看视频| 搡老岳熟女国产| 亚洲不卡免费看| 国产精品一及| 麻豆一二三区av精品| 久久久久久久午夜电影| 国产精品99久久久久久久久| 午夜福利在线观看吧| 亚洲成人av在线免费| 亚洲四区av| 国产乱人视频| 如何舔出高潮| 日韩欧美三级三区| 亚洲丝袜综合中文字幕| 免费大片18禁| 国产高清不卡午夜福利| 国产精品一区www在线观看| 日产精品乱码卡一卡2卡三| 国产精品久久久久久久久免| 香蕉av资源在线| 免费人成在线观看视频色| 中文字幕熟女人妻在线| 精品福利观看| 国产淫片久久久久久久久| 波多野结衣高清无吗| 亚洲成a人片在线一区二区| 久久国产乱子免费精品| 一本精品99久久精品77| 一级毛片电影观看 | 男人狂女人下面高潮的视频| 菩萨蛮人人尽说江南好唐韦庄 | 男插女下体视频免费在线播放| 国产午夜福利久久久久久| 乱系列少妇在线播放| 最好的美女福利视频网| 校园春色视频在线观看| 久久久久久伊人网av| 黄色欧美视频在线观看| 久久久精品大字幕| 国产精品日韩av在线免费观看| 国产精品人妻久久久久久| 两个人视频免费观看高清| 国产亚洲精品久久久久久毛片| 美女cb高潮喷水在线观看| 欧美+亚洲+日韩+国产| 晚上一个人看的免费电影| 久久久久久久久大av| av专区在线播放| 超碰av人人做人人爽久久| 精品福利观看| 能在线免费观看的黄片| 日韩中字成人| 国产在线精品亚洲第一网站| 久久久精品94久久精品| 色综合色国产| 麻豆成人午夜福利视频| 国产精品无大码| 中文字幕人妻熟人妻熟丝袜美| 亚洲18禁久久av| 搡老熟女国产l中国老女人| 69人妻影院| 国产精品久久久久久久电影| 麻豆乱淫一区二区| 搡女人真爽免费视频火全软件 | 久久这里只有精品中国| 联通29元200g的流量卡| 观看美女的网站| 国产熟女欧美一区二区| 日本撒尿小便嘘嘘汇集6| 搡老岳熟女国产| 男女那种视频在线观看| 乱人视频在线观看| 国产老妇女一区| 中文字幕精品亚洲无线码一区| 搡老岳熟女国产| 国产精品一区二区免费欧美| 国产一区二区亚洲精品在线观看| 伊人久久精品亚洲午夜| 亚洲av免费高清在线观看| 国产av在哪里看| 99热这里只有精品一区| 国产一区二区在线av高清观看| 尾随美女入室| 18禁裸乳无遮挡免费网站照片| 一进一出抽搐gif免费好疼| 如何舔出高潮| 亚洲熟妇熟女久久| 国产伦在线观看视频一区| 国模一区二区三区四区视频| 国产精品久久久久久亚洲av鲁大| 人妻久久中文字幕网| 在线免费观看的www视频| 成人精品一区二区免费| 日本爱情动作片www.在线观看 | 少妇熟女欧美另类| 婷婷六月久久综合丁香| 国产精品电影一区二区三区| 校园人妻丝袜中文字幕| 亚洲一区二区三区色噜噜| 一级av片app| 日韩三级伦理在线观看| 国产探花极品一区二区| 欧美色视频一区免费| 村上凉子中文字幕在线| 午夜久久久久精精品| 中文字幕熟女人妻在线| 国产熟女欧美一区二区| 啦啦啦啦在线视频资源| 国产精品人妻久久久影院| 国产成人一区二区在线| 日本欧美国产在线视频| 日韩中字成人| 天堂av国产一区二区熟女人妻| 日本a在线网址| 国内精品美女久久久久久| 91精品国产九色| 国产国拍精品亚洲av在线观看| 搡老熟女国产l中国老女人| 亚洲高清免费不卡视频| 欧美+亚洲+日韩+国产| 午夜福利成人在线免费观看| 国产一区二区亚洲精品在线观看| 国产探花极品一区二区| 亚洲欧美中文字幕日韩二区| 日本 av在线| 99热精品在线国产| 久久韩国三级中文字幕| 日韩人妻高清精品专区| 尤物成人国产欧美一区二区三区| 91久久精品国产一区二区三区| 伦理电影大哥的女人| 51国产日韩欧美| 菩萨蛮人人尽说江南好唐韦庄 | 毛片女人毛片| 五月伊人婷婷丁香| 天堂av国产一区二区熟女人妻| 最近中文字幕高清免费大全6| 国语自产精品视频在线第100页| 亚洲国产精品成人久久小说 | 亚洲七黄色美女视频| 一区二区三区高清视频在线| 内地一区二区视频在线| 久久久精品欧美日韩精品| 人人妻人人澡人人爽人人夜夜 | 变态另类成人亚洲欧美熟女| 国产精品一二三区在线看| 五月玫瑰六月丁香| 伦理电影大哥的女人| 一级黄片播放器| 亚洲av美国av| 亚洲国产精品合色在线| 别揉我奶头~嗯~啊~动态视频| 看非洲黑人一级黄片| 亚洲人成网站在线观看播放| 在线观看午夜福利视频| 日韩人妻高清精品专区| 中文在线观看免费www的网站| 日韩强制内射视频| 欧美日韩一区二区视频在线观看视频在线 | 日日干狠狠操夜夜爽| 1000部很黄的大片| 国内久久婷婷六月综合欲色啪| 国产精品爽爽va在线观看网站| 波野结衣二区三区在线| 欧美3d第一页| 国产一区二区激情短视频| 一个人看的www免费观看视频| 国产成人aa在线观看| 国产成人91sexporn| 人妻少妇偷人精品九色| 在线观看一区二区三区| 日韩一区二区视频免费看| 久久久久国内视频| 欧美日本亚洲视频在线播放| 十八禁国产超污无遮挡网站| 69人妻影院| 国产精品久久久久久亚洲av鲁大| 久久精品国产亚洲av天美| 亚洲最大成人av| 人妻夜夜爽99麻豆av| 国产乱人偷精品视频|