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

    一維準周期鑲嵌模型的擴展態(tài)

    2024-12-31 00:00:00張彪張艷陽
    科技創(chuàng)新與應用 2024年26期
    關鍵詞:電導

    摘" 要:最近的理論研究發(fā)現(xiàn),一維準周期鑲嵌模型存在遷移率邊。在該文中,研究該模型的量子輸運,尤其是擴展態(tài)的兩端口電導,模型的輸運相圖、平均值和統(tǒng)計分布。擴展態(tài)的電導并不恒為1,而是有干涉共振條紋,用波函數(shù)的行為解釋這些干涉條紋。該模型的擴展態(tài)能在極強的準周期勢能下存活,通過電導的統(tǒng)計分布,以及與波函數(shù)分形維度的比較,解釋其中的物理圖像。

    關鍵詞:準周期鑲嵌模型;擴展態(tài);電導;介觀輸運;量子輸運

    中圖分類號:O413" " " 文獻標志碼:A" " " " " 文章編號:2095-2945(2024)26-0032-06

    Abstract: Recent theoretical studies have found that there is a mobility edge in the one-dimensional quasi-periodic mosaic model. This paper studies the quantum transport of the model, especially the two-port conductance of the extended state, its transport phase diagram, average value and statistical distribution. The conductance of the extended state is not always 1, but there are interference resonance fringes. We explain these interference fringes by the behavior of wave function. The extended state of the model can survive under extremely strong quasi-periodic potential energy. Explain the physical images through the statistical distribution of conductance and the comparison with the fractal dimension of wave function.

    Keywords: quasiperiodic mosaic model; extended state; conductance; mesoscopic transport; quantum transport

    量子輸運是量子相干效應明顯時發(fā)生的輸運現(xiàn)象,其可以發(fā)生于低溫的普通晶體材料中,也可以發(fā)生于光子晶體等現(xiàn)代量子材料中。量子輸運是微觀粒子(如電子)的波函數(shù)經(jīng)歷晶格和雜質(zhì)相干散射的結(jié)果,因此晶體的細節(jié),例如周期性或非周期性都會對量子輸運產(chǎn)生很大的影響。根據(jù)標準的固體能帶理論,周期晶格中的波函數(shù)都是擴展的布洛赫波,有良好的導電性。然而,無序會破壞晶格的周期性,同時影響波函數(shù)的擴展性。這個影響與系統(tǒng)的空間維度密切相關[1-2]。在三維晶格中,弱無序不會改變波函數(shù)的擴展性,但強無序會導致金屬-絕緣體相變(或叫“擴展態(tài)-局域態(tài)轉(zhuǎn)變”),能量軸上的相變點稱為遷移率邊。不過,在一維或二維晶格中,任何有限的無序都會使波函數(shù)局域化,從而導致全部狀態(tài)都處于絕緣相[1-2]。這些都是凝聚態(tài)物理中熟知的結(jié)論。

    準周期結(jié)構微妙地介于無序和周期之間,是傳統(tǒng)固體電子理論較少研究的對象,其對量子態(tài)的影響是現(xiàn)代凝聚態(tài)物理感興趣的話題。近期最著名的準周期結(jié)構出現(xiàn)在轉(zhuǎn)角石墨烯中,兩層石墨烯之間存在一個非公度角度的時候,晶格形成準周期結(jié)構[3-4]。這個二維準周期結(jié)構會在狄拉克點附近形成平帶,從而有助于超導、強關聯(lián)等豐富物理現(xiàn)象的發(fā)生。

    在一維系統(tǒng)中,準周期對波函數(shù)的作用也是近期的一個研究熱點[5]。在著名的Aubry-André-Harper模型(簡稱AAH模型)中,人們在一維緊束縛模型上施加了一個周期勢,但其周期與晶格的周期不公度,從而構成準周期勢。理論研究發(fā)現(xiàn),如果準周期勢足夠弱的話,即便在一維也能出現(xiàn)擴展態(tài)[6]。當準周期勢能大于臨界值后,本征態(tài)會被局域化。但這個模型中的擴展態(tài)-局域態(tài)轉(zhuǎn)變對能量軸上的所有本征態(tài)是同時發(fā)生的,因此沒有遷移率邊。

    最近,Wang等[7]對AAH模型進行了改進,設計了一個被稱為鑲嵌模型的一維結(jié)構,即周期性(周期大于1)地讓格點的勢能為零。這個看似簡單的改變就能讓系統(tǒng)出現(xiàn)遷移率邊:能帶中心為擴展態(tài),能帶邊緣為局域態(tài)。

    在本論文中,對這個鑲嵌模型的量子輸運進行了更細致的理論研究,尤其是其中擴展態(tài)(即金屬態(tài))的行為。計算了該模型不同參數(shù)下的分形維度和兩端口電導,并詳細討論了其統(tǒng)計漲落和尺寸標度效應。研究發(fā)現(xiàn)擴展態(tài)確實對應著非零的電導,但電導的行為有很多有趣的細節(jié),例如干涉共振起伏與波函數(shù)的關系。在能帶中心,該模型的擴展態(tài)可以在非常大的勢能強度下存活,但這個過程中,其能帶結(jié)構和電導的統(tǒng)計性質(zhì)會有對應的改變。我們分析和解釋了擴展態(tài)的特征。

    1" 模型和方法

    許多原來只出現(xiàn)在理論研究中的一維晶格模型,目前已經(jīng)能夠通過冷原子、光晶格等現(xiàn)代實驗方法實現(xiàn)。其中的周期、無序或準周期結(jié)構也能夠通過越來越成熟的技術進行設計和調(diào)控。在本文中,考慮一維晶格中的準周期鑲嵌模型,其緊束縛哈密頓量可以描述為[7]

    H=t(cc+H.c.)+2λjnj , (1)

    式中:c(cj)是費米子在第j個格點處的產(chǎn)生(湮滅)算符;nj是該處的粒子數(shù)算符(或勢能算符)。t是最近相鄰格點的跳躍系數(shù),把其設為1作為能量單位。同時把晶格常數(shù)(最近鄰格點間的距離)設為1作為長度單位。λj是第j個格點上的準周期鑲嵌勢能,由下面的式子給出

    λj=λcos[2π(ωj+θ)],j=mκ

    0,其他, (2)

    式中:λ是勢能的強度,θ是勢能的空間初相位。參數(shù)ω決定了勢能的周期,如果其是無理數(shù)就能保證勢能周期與晶格周期非公度,形成研究的準周期勢能。如無明確指出,設置θ和ω的值分別為: θ=0,ω=(-1)/2。κ是決定鑲嵌調(diào)制周期的整數(shù),當κ=1時,即為最初的AAH模型,不具備遷移率邊。當κ≠1時,此時模型有鑲嵌,即可出現(xiàn)遷移率邊[7]。本文中κ=2。取 m=1,2,...,N,此時系統(tǒng)長度為L=κN。對于這個模型,文獻[7]中已經(jīng)給出了不同勢能強度時候的遷移率邊為Ec=±1/λ,也就是說,當能量|E|小于1/λ時,本征態(tài)處于擴展態(tài),而|E|大于1/λ時的本征態(tài)處于局域態(tài)。有了哈密頓量,就可以通過標準的對角化方法求得本征能量和相應的本征態(tài)。為了表征本征態(tài)的局域化性質(zhì),引入逆參與率(IPR),其定義式為[7-8]

    IPR(m)=Ψm,j, (3)

    式中:Ψm,j表示第m個本征態(tài)。典型的擴展態(tài)和局域態(tài)分別對應著IPR接近0和1。通過IPR還可以定義另一個更直觀的量,即分形維度,表達式如下

    Γ=- 。 (4)

    在一維晶格中,趨于0和1分別對應局域態(tài)和擴展態(tài)。

    以上物理量研究的都是孤立晶格中波函數(shù)的性質(zhì)。波函數(shù)的局域化程度會直接影響量子輸運。設想把一定長度的該模型晶格(即研究對象)接上電極,就可以測量其電導。在零溫時,可以利用著名的Landauer 公式來模擬計算此時的電導,其正比于費米子的量子透射率,用格林函數(shù)表達是[9-10]

    (5)

    式中:e是電子電荷;h是普朗克常數(shù);τ是透射率;ES是測量時研究對象中的費米能。ΓL和ΓR分別是左右電極的譜函數(shù)矩陣,其計算公式為

    ΓL/R(Elead)=i(ΣL/R (Elead)-Σ(Elead),(6)

    式中:ΣL/R(Elead)是描述左右電極自能的矩陣;Elead是電極中的費米能。Gr(E)和Ga(E)分別是推遲和超前格林函數(shù),計算公式為

    Gr/a(E)=[E-H-Σr/a(Elead)]-1,(7)

    式中:H為研究對象的哈密頓量;E為本征能量。

    實驗中測到的電導是經(jīng)歷了多個相位相干長度后的平均結(jié)果,這在無序系統(tǒng)中相當于取了不同的無序構型進行平均[11-12]。在這里的準周期結(jié)構中,為了類比無序系統(tǒng)中的無序構型樣本,在其他關鍵模型參數(shù)取定的情況下,采用2種抽樣和平均方式。第一種是在公式(2)中取一組隨機的θ,簡稱為“相位平均”。第二種是在公式(5)中取一組隨機的電極費米能Elead,簡稱為“能量平均”。這2個參數(shù)的改變都能導致電導的漲落。然后就可以得到相應電導的平均值、方差和統(tǒng)計分布,從而獲得量子輸運的更豐富信息。

    2" 結(jié)果和討論

    首先來研究孤立體系的性質(zhì),即本征波函數(shù)的分形維數(shù)。由之前的介紹得知,當分形維度Γ趨于1時,模型處于擴展態(tài),這引發(fā)了筆者對于分形維度與樣本長度之間關系的深入探討。

    在圖1中,取勢能強度λ=1,畫出了本征函數(shù)的分形維度Γ作為本征能量E的函數(shù),不同形狀的點代表了不同的樣品長度。由于計算結(jié)果對于能量原點是對稱的,因此只畫了正能量部分的結(jié)果。圖中第一個明顯的特征是:本征能量的分布并不完全連續(xù),而是被分割成了幾個子能帶。這是鑲嵌勢能的一個特點[7]。

    現(xiàn)在仔細觀察數(shù)值點的細節(jié)??梢钥吹剑谶w移率邊Ec=1的兩側(cè),分形維度的尺寸標度確實具有不同的特征。在遷移率邊左側(cè)的本征態(tài),其分形維度數(shù)值都隨著尺寸N的增加而增大。例如,菱形數(shù)據(jù)點(N=25 600)位于圓形數(shù)據(jù)點(N=400)的上方。因此可以預見在熱力學極限N→∞時,分形維度將趨于1,顯示出明顯的擴展態(tài)行為。另一方面,在遷移率邊右側(cè)的本征態(tài),其標度行為則是相反的:菱形數(shù)據(jù)點(N=25 600)位于圓形數(shù)據(jù)點(N=400)的下方,也就是分形維度隨著尺寸增大而減小,從而在熱力學極限下趨于0,顯示出標準的局域態(tài)行為。

    以上是一個固定λ值時孤立體系的波函數(shù)局域化性質(zhì),應證了解析計算的結(jié)果。現(xiàn)在轉(zhuǎn)到輸運性質(zhì)。在圖2中,畫出了平均電導隨著樣品能量與勢能的變化的相圖。圖2(a)、圖2(c)、圖2(e)表示電極能量EL的平均的結(jié)果,圖2(b)、圖2(d)、圖2 (f)表示相位偏移θ的平均的結(jié)果。圖2(a)、圖2(b)中的Nx=1 000, 圖2(c)、圖2(d)中的Nx=2 000,圖2(e)、圖2(f)中的Nx=4 000。圖中的淺色和深色分別代表電導接近1和0。首先可以看到,2種平均下亮色的區(qū)域都是一致的,而且這些區(qū)域的范圍并不隨著長度增長而縮小,因此代表了穩(wěn)定的擴展態(tài)。這些擴展態(tài)的分布區(qū)域與分形維度[5]的計算結(jié)果一致。因此,量子輸運確實忠實地體現(xiàn)了波函數(shù)的性質(zhì)。

    Nx=1 000,2 000,4 000,朗道電導的平均值L作為樣本中心能量ES和準周期勢強度λ的函數(shù),圖2(a)(c)(e)表示電極能量EL的改變,圖2(b)(d)(f)表示相位偏移θ的改變。

    進一步觀察發(fā)現(xiàn),淺色區(qū)域的電導盡管都遠大于0,但并不都等于1,而是有起伏,形成精致的干涉共振條紋。其中最顯眼的是ES=0的中心線上的淺色的點,代表了最明顯的共振峰。在相同長度的情況下,共振峰的位置對2種平均都是相同的,但共振峰的個數(shù)隨著長度的增長而增加。

    為了更好地理解電導和波函數(shù)之間的對應關系,例如圖2所示的共振峰的來源,比較實空間波函數(shù)的起伏。與電導隨長度變化的關聯(lián)性。為了反映波函數(shù)對于樣品空間的依賴細節(jié),選取一個位于擴展態(tài)區(qū)域的本征能量,計算波函數(shù)模平方的傅里葉變換,其包含了波函數(shù)實空間漲落的特征性信息。與此同時計算該樣品處于這個費米能時,電導作為樣品長度的依賴關系,并且對其也作傅里葉變換。在圖3中,選取了2個典型的能量,分別用黑線和灰線畫出了這2個量的傅里葉變換??梢钥吹?,2條曲線在多數(shù)情況下具有相同的峰值位置,這就說明了電導的起伏確實來源于波函數(shù)的起伏,兩者是高度相關的。剩下有少數(shù)峰值不相關,可以歸結(jié)為樣品邊界的散射影響,因為計算波函數(shù)時采用的是固定邊界條件,但計算電導時樣品邊界則連著半無限電極。這2種不同邊界對波函數(shù)的散射,會導致一些不同的相干空間漲落,于是形成不同的傅里葉峰。

    上面討論的是電導的平均值,現(xiàn)在來研究其統(tǒng)計漲落。眾所周知,電導的統(tǒng)計性質(zhì)與局域化性質(zhì)密切相關[2,11-12]。在圖4中展現(xiàn)了電極能量改變時電導的統(tǒng)計直方圖,4個子圖對應著依次增大的勢函數(shù)強度λ。在圖4(a)中,極弱勢能時(λ=0.125 63),電導集中分布高度集中在接近1的區(qū)域,也就是接近完美透射的擴展態(tài)。隨著勢能增大,如圖4(b)和圖4(c)所示,電導分布逐漸展寬,直到占滿區(qū)間(0,1)。與此同時,統(tǒng)計分布的重心往低電導方向移動。最后,強勢能時,如圖4(d)所示,電導出現(xiàn)一個峰值接近0的長尾分布。此時盡管電導的平均值比較小,但即便在很長樣品的情況下(例如N=100 000)這個平均值仍然不是0,而且時不時還會出現(xiàn)一些電導遠大于零的統(tǒng)計樣本。這對于只有一個通道、背散射風險很大的一維時間反演模型已經(jīng)非常不容易了[2,11]??傊藭r系統(tǒng)仍然處于良好定義的擴展態(tài)。與此相反,在|E|gt;1/λ時的局域態(tài),即便在弱勢能、樣品長度不大的情況下,所有統(tǒng)計樣品的兩端口電導都非常接近零了。

    對于周期系統(tǒng),良好定義的能帶結(jié)構與輸運性質(zhì)之間有良好的對應關系[2-3]。雖然嚴格來說,準周期結(jié)構不具備周期性,但仍然可以用周期結(jié)構來逼近準周期結(jié)構,如斐波那契結(jié)構[5,8],從而獲得能帶結(jié)構。為了更加深入地理解不同勢能強度下輸運特性與量子態(tài)的內(nèi)在聯(lián)系,利用斐波那契晶格來計算準周期勢能系統(tǒng)的能帶結(jié)構,具體計算方法如下。取一個長度為第n個斐波那契數(shù)Fn的樣本,即L=Fn,此時取公式(2)中勢能的周期參數(shù)ω=ωn=Fn-1/Fn,其是一個有理數(shù)。于是就可以以這個長度為Fn的樣品作為超元胞,按照標準的傅里葉變換來計算能帶結(jié)構,即色散關系E(k)。當n趨于無窮時,ωn也趨于無理數(shù)(-1)/2,從而使得樣品趨于本文研究的準周期勢模型[5,8]。

    圖5就是按照這種方法計算的能帶結(jié)構,超元胞長度L=F17=1 597。只畫出了關心的能帶中心附近的幾個子能帶。4個子圖分別對應著圖4中4個子圖的勢能強度λ。在圖5(a)展示的弱勢能情形,子能帶顯示出接近線性的色散關系,并且子帶之間的子帶隙非常小。根據(jù)Thouless公式,電導正比于群速度(即色散關系的斜率)與態(tài)密度的乘積[13]。因此,圖5(a)中的這種能帶結(jié)構,即接近于常數(shù)的群速度與極小的帶隙,在2個因素上都有利于形成較大的電導,并且攜帶較大電導的本征態(tài)稠密地分別在能量軸上。這就是本文在前面的圖4(a)中看到的現(xiàn)象。

    隨著勢能的增強,在圖5(b)—圖5(d)中,子能帶變得越來越平,子帶隙也變得越來越大。這種變化趨勢可能會導致量子態(tài)的輸運能力下降。但是,系統(tǒng)通過2個方面的努力來讓自己避免被完全局域化。首先,子能帶在強勢能的時候雖然被壓制得很窄,但其寬度并不是零,這與局域態(tài)有數(shù)量上的明顯差別。其次,子能帶之間變得更密集(阻止態(tài)密度過度降低),以此來彌補群速度降低而給Thouless電導帶來的損失。因此在強勢能情況下,最后的總效果就是,能帶中心仍然存在著擴展態(tài),但擴展態(tài)之間的子帶隙仍然比弱勢能情況下的更大。這就導致一個現(xiàn)象:當費米面穿過頑強的子能帶時,電導遠大于0;當費米面穿過子帶隙時,電導接近于0。這就造成了在圖4(b)—圖4(d)中顯示的電導漲落增強。

    上面這個物理圖像,還可以在更強的勢能強度下得到進一步的印證??梢钥吹?,在準周期勢λ遠大于帶寬的情形下,這個一維模型在ES=0附近仍然存在擴展態(tài),這是比較特別的。為了更清楚地看到大λ情況下發(fā)生了什么,在圖6中畫出了E=0附近的分形維度(圖6(a))和電導圖(圖6(b)),但擴大了λ的掃描范圍。從圖6(a)的分形維度可以看到,數(shù)值接近1的擴展態(tài)(淺灰色)在λ接近100的情況下仍然非常頑強地存在著。但與此同時,圖6(b)的電導在λ不到10的時候就接近0(黑色),似乎變成局域態(tài)了。但放大后就能看到,此時是在黑色(電導接近0)的海洋里點綴著淺灰色(電導接近1)的孤立點。

    這些現(xiàn)象再一次反映了在圖4和圖5中看到的物理圖像。擴展態(tài)能帶分裂為多個擴展態(tài)子能帶。隨著λ的增加,子能帶逐漸變窄,它們之間的子帶隙則變寬。于是,當掃描能量的時候,如果掃到了子能帶里,就能得到電導接近1;如果掃到了帶隙,就得到電導接近0。這就是圖3中電導分布越來越寬、平均值越來越接近0的原因。但無論如何,只要能量在子能帶里,電導就不會是0,因此也仍然是良好定義的擴展態(tài)。有理由相信,在λ趨于無窮時,仍然會在ES=0附近存在一個無限細的擴展態(tài)子能帶。

    3" 結(jié)論

    本文研究了一維準周期嵌套模型的量子輸運性質(zhì),得到以下結(jié)論。

    1)擴展態(tài)的電導不為零,且具有隨著能量和勢能強度改變的干涉共振起伏。這些起伏的大部分都與波函數(shù)的起伏對應,少部分來源于界面處的反射與干涉。

    2)在主能帶中心,隨著勢能增長,擴展態(tài)一直存在,但能量分布范圍縮小。這種分布可以緩解由于子能隙擴大而導致的態(tài)密度下降,并彌補子能帶變平而導致電導的過度下降。

    3)在此過程中,電導的漲落逐漸增大,平均值也在降低,但不為零。原因是擴展態(tài)子帶的帶寬逐漸變?。ǖ粫?yōu)榱悖訋吨饾u變大。于是當費米能碰到子帶(帶隙)的時候就表現(xiàn)為接近1(0)的數(shù)值。

    參考文獻:

    [1] E. ABRAHAMS, P.W. ANDERSON, D. C. LICCIARDELLO, et al.. Scaling theory of localization: Absence of quantum diffusion in two dimensions[J].Phys.Rev.Lett.42,673 (1979).

    [2] P. A. LEE, T. V. RAMAKRISHNAN. Disordered electronic systems[J]. Rev. Mod. Phys.57,287(1985).

    [3] B. BISTRITZER, A. H. MACDONALD. Moiré butterflies in twisted bilayer graphene[J].Phys. Rev. B 84,035440 (2011).

    [4] Y. CAO, V. FATEMI, A. DEMIR, et al. Correlated insulator behaviour at half-filling in magic-angle graphene superlattices[J]. Nature 556,80(2018).

    [5] A. JAGANNATHAN. The Fibonacci quasicrystal: Case study of hidden dimensions and multifractality[J]. Rev. Mod. Phys. 93,045001(2021).

    [6] S. AUBRY, G. ANDRé. Analyticity breaking and Anderson localization in incommensurate lattices[J].Ann. Isr. Phys. Soc. 3,133(1980).

    [7] Wang Y, Xia X, Zhang L, et al. One-dimensional quasiperiodic mosaic lattice with exact mobility edges[J]. Phys. Rev. Lett.125,196604 (2020).

    [8] WANG Y, ZHANG L, SUN W, et al. Quantum phase with coexisting localized, extended, and critical zones[J]. Phys. Rev. B 106,L140203 (2022).

    [9] S. DATTA. Electronic transport in mesoscopic systems[M]. Cambridge: Cambridge University Press, 1995.

    [10] S. DATTA. Quantum transport: Atom to transistor [M].Cambridge: Cambridge University Press, 2005.

    [11] F. EVERS,A. D. MIRLIN. Anderson transitions[J]. Rev. Mod. Phys. 80, 1355 (2008).

    [12] LUO X,T. OHTSUKI. Universality classes of the Anderson transitions driven by quasiperiodic potential in the three-dimensional Wigner-Dyson symmetry classes[J]. Phys. Rev. B 106,104205 (2022).

    [13] J. T. Edwards,D. G. Thoules." Numerical studies of localization in disordered systems[J]. J. Phys.C 5,807 (1972).

    基金項目:廣東省自然科學基金(2023A1515010698)

    第一作者簡介:張彪(1997-),男,碩士研究生。研究方向為量子輸運。

    *通信作者:張艷陽(1978-),男,博士,副教授。研究方向為量子輸運。

    猜你喜歡
    電導
    超聲脈沖電導結(jié)合關節(jié)松動術治療肩周炎的臨床研究
    健康護理(2022年3期)2022-05-26 02:27:49
    金寶AK96血液透析機的電導度故障分析與處理
    基于IEC標準的電阻表(阻抗表)和電導表的技術要求研究
    電子制作(2018年14期)2018-08-21 01:38:38
    Cd組分x對非晶態(tài)Hg1-xCdxTe薄膜暗電導的影響
    紅外技術(2017年1期)2017-03-27 09:02:20
    基于電導增量法的模型預測控制光伏MPPT算法
    鹽酸羥考酮緩釋片聯(lián)合超聲電導儀治療中重度癌痛的臨床研究
    RNA干擾HeLa細胞IKCa1基因?qū)χ须妼р}激活鉀通道電流的影響
    金寶AK200S透析機透析液濃度故障維修一例
    多磺酸黏多糖乳膏聯(lián)合超聲電導儀治療靜脈炎30例
    電導測試在蓄電池測試中的應用
    科技資訊(2012年3期)2012-04-29 22:39:01
    亚洲精品色激情综合| 757午夜福利合集在线观看| 成人国语在线视频| 国产精品自产拍在线观看55亚洲| 久久精品人妻少妇| 中文字幕av在线有码专区| av福利片在线| 高清毛片免费观看视频网站| 中文字幕人妻丝袜一区二区| 日韩中文字幕欧美一区二区| 啦啦啦观看免费观看视频高清| 亚洲成人久久爱视频| 精品国产乱码久久久久久男人| 午夜福利欧美成人| 777久久人妻少妇嫩草av网站| 久久人妻av系列| 青草久久国产| 久久久久久国产a免费观看| 午夜福利高清视频| 久久久久免费精品人妻一区二区| 精品午夜福利视频在线观看一区| 欧美黑人巨大hd| 男女床上黄色一级片免费看| 久久精品人妻少妇| 欧美性长视频在线观看| 伦理电影免费视频| 欧美日韩福利视频一区二区| 精品国产乱码久久久久久男人| 久久婷婷人人爽人人干人人爱| 在线观看www视频免费| 一区二区三区国产精品乱码| 久久久国产成人精品二区| a级毛片在线看网站| 老汉色av国产亚洲站长工具| 午夜福利在线在线| 激情在线观看视频在线高清| 久久热在线av| 18禁黄网站禁片午夜丰满| 91av网站免费观看| 国产精品久久久人人做人人爽| 久久婷婷成人综合色麻豆| 一本大道久久a久久精品| 成熟少妇高潮喷水视频| 精品久久久久久成人av| 欧美色视频一区免费| 国产成人精品无人区| videosex国产| 女生性感内裤真人,穿戴方法视频| 男女下面进入的视频免费午夜| www.自偷自拍.com| 欧美成人午夜精品| 久久精品亚洲精品国产色婷小说| 国产精品一区二区免费欧美| 成人国语在线视频| 午夜精品久久久久久毛片777| 欧美zozozo另类| 19禁男女啪啪无遮挡网站| 一区福利在线观看| 国产三级中文精品| 久久精品91蜜桃| 99国产精品一区二区三区| 国产精品野战在线观看| 国产精品 欧美亚洲| 欧美3d第一页| 在线观看日韩欧美| 一边摸一边抽搐一进一小说| 香蕉丝袜av| 国产一区二区在线av高清观看| 国产免费av片在线观看野外av| 村上凉子中文字幕在线| 久久婷婷成人综合色麻豆| 亚洲人成网站在线播放欧美日韩| 天天一区二区日本电影三级| netflix在线观看网站| 怎么达到女性高潮| 国产激情久久老熟女| 国产成人啪精品午夜网站| 又黄又粗又硬又大视频| 18禁观看日本| 午夜福利18| 久久精品综合一区二区三区| 99riav亚洲国产免费| 日本撒尿小便嘘嘘汇集6| 18禁美女被吸乳视频| 久久香蕉国产精品| 最近最新中文字幕大全免费视频| 日本一二三区视频观看| 97超级碰碰碰精品色视频在线观看| 精品久久久久久久久久久久久| 成人一区二区视频在线观看| 日韩国内少妇激情av| 日韩欧美在线乱码| 欧美成狂野欧美在线观看| 99久久综合精品五月天人人| 丰满人妻熟妇乱又伦精品不卡| 亚洲av成人一区二区三| 亚洲欧美精品综合一区二区三区| 91九色精品人成在线观看| 久久香蕉精品热| 999精品在线视频| 丁香六月欧美| 国产精品电影一区二区三区| 日本三级黄在线观看| 免费观看精品视频网站| 最近最新中文字幕大全电影3| 可以在线观看毛片的网站| 亚洲免费av在线视频| 亚洲欧美精品综合久久99| 成熟少妇高潮喷水视频| 91大片在线观看| 人人妻人人看人人澡| 久久精品亚洲精品国产色婷小说| 亚洲 欧美 日韩 在线 免费| 18禁观看日本| 999久久久国产精品视频| 国产午夜福利久久久久久| 精品不卡国产一区二区三区| 精品一区二区三区视频在线观看免费| 不卡av一区二区三区| 少妇粗大呻吟视频| 亚洲精品在线美女| 夜夜看夜夜爽夜夜摸| 精品少妇一区二区三区视频日本电影| 成人高潮视频无遮挡免费网站| 国产精品免费视频内射| 婷婷精品国产亚洲av在线| 日韩精品免费视频一区二区三区| 少妇被粗大的猛进出69影院| 波多野结衣高清作品| 亚洲av成人av| 日本免费a在线| 日韩av在线大香蕉| 日韩高清综合在线| 可以在线观看的亚洲视频| 国产欧美日韩一区二区精品| 欧美丝袜亚洲另类 | 成人永久免费在线观看视频| 国产免费男女视频| 香蕉久久夜色| 亚洲一区中文字幕在线| 不卡一级毛片| 亚洲自拍偷在线| 97人妻精品一区二区三区麻豆| 国产av不卡久久| 90打野战视频偷拍视频| 国产亚洲精品久久久久5区| 麻豆av在线久日| 最近在线观看免费完整版| 亚洲欧美激情综合另类| 99久久综合精品五月天人人| 最新美女视频免费是黄的| 国产激情偷乱视频一区二区| 亚洲欧美日韩东京热| 激情在线观看视频在线高清| 亚洲精品中文字幕在线视频| 国产精品久久久久久亚洲av鲁大| 午夜福利18| 他把我摸到了高潮在线观看| 一二三四社区在线视频社区8| 淫秽高清视频在线观看| 狠狠狠狠99中文字幕| 最近最新免费中文字幕在线| 欧美人与性动交α欧美精品济南到| 首页视频小说图片口味搜索| 成人亚洲精品av一区二区| www.精华液| 久久国产精品人妻蜜桃| 免费在线观看日本一区| 国产1区2区3区精品| 亚洲最大成人中文| 国产亚洲精品第一综合不卡| av在线播放免费不卡| 一进一出抽搐gif免费好疼| 亚洲欧美精品综合一区二区三区| 国产精品电影一区二区三区| 99久久99久久久精品蜜桃| 成人国产综合亚洲| 亚洲色图av天堂| av福利片在线观看| 国产一区二区三区视频了| 国产三级中文精品| 国产一区二区在线观看日韩 | 日韩高清综合在线| 极品教师在线免费播放| 一个人免费在线观看的高清视频| 日韩欧美免费精品| 99久久无色码亚洲精品果冻| 亚洲精品久久国产高清桃花| 欧美色视频一区免费| x7x7x7水蜜桃| 中出人妻视频一区二区| 国产精品 欧美亚洲| 国产在线精品亚洲第一网站| 看免费av毛片| 99久久99久久久精品蜜桃| 国产精品乱码一区二三区的特点| 国产精品98久久久久久宅男小说| 国产午夜福利久久久久久| 91字幕亚洲| 国产不卡一卡二| 久久久久久久久久黄片| 久久精品综合一区二区三区| 日日干狠狠操夜夜爽| 日韩高清综合在线| 国产麻豆成人av免费视频| 女人被狂操c到高潮| 级片在线观看| 精品高清国产在线一区| 亚洲国产精品合色在线| 国产一级毛片七仙女欲春2| 午夜影院日韩av| 香蕉av资源在线| 中文字幕熟女人妻在线| 啪啪无遮挡十八禁网站| 别揉我奶头~嗯~啊~动态视频| 欧美色欧美亚洲另类二区| 欧美另类亚洲清纯唯美| 国内精品久久久久精免费| 欧美 亚洲 国产 日韩一| 国产精品久久久久久久电影 | 首页视频小说图片口味搜索| 中文在线观看免费www的网站 | АⅤ资源中文在线天堂| 51午夜福利影视在线观看| 国产成人影院久久av| 热99re8久久精品国产| 一卡2卡三卡四卡精品乱码亚洲| 成人18禁高潮啪啪吃奶动态图| 一区二区三区国产精品乱码| 少妇熟女aⅴ在线视频| 久久精品91蜜桃| 午夜影院日韩av| 国产私拍福利视频在线观看| 精品久久久久久久毛片微露脸| 国产成人精品无人区| 别揉我奶头~嗯~啊~动态视频| 精品久久蜜臀av无| 黄色a级毛片大全视频| 亚洲自偷自拍图片 自拍| 精品一区二区三区视频在线观看免费| 女人高潮潮喷娇喘18禁视频| 久久天堂一区二区三区四区| 国产精品爽爽va在线观看网站| 一个人观看的视频www高清免费观看 | 搡老熟女国产l中国老女人| 黑人巨大精品欧美一区二区mp4| 人人妻人人澡欧美一区二区| 岛国视频午夜一区免费看| 黄频高清免费视频| 真人一进一出gif抽搐免费| 成人三级做爰电影| 国产激情偷乱视频一区二区| 国产精品一区二区精品视频观看| 好男人电影高清在线观看| 两个人免费观看高清视频| 亚洲av日韩精品久久久久久密| 久久精品成人免费网站| 久久久久久亚洲精品国产蜜桃av| 中文字幕高清在线视频| 小说图片视频综合网站| 村上凉子中文字幕在线| 亚洲精品美女久久av网站| 校园春色视频在线观看| 午夜精品一区二区三区免费看| www日本在线高清视频| 我的老师免费观看完整版| 亚洲国产精品sss在线观看| 免费高清视频大片| 夜夜爽天天搞| 午夜老司机福利片| 天堂影院成人在线观看| 国产免费男女视频| 午夜福利高清视频| 欧美午夜高清在线| 日本免费a在线| 欧美日韩福利视频一区二区| 日韩欧美国产在线观看| 日本熟妇午夜| 国产三级中文精品| 中文字幕人妻丝袜一区二区| av超薄肉色丝袜交足视频| 国产区一区二久久| 丁香六月欧美| 精品国产超薄肉色丝袜足j| 国产精品香港三级国产av潘金莲| 亚洲无线在线观看| 久久午夜亚洲精品久久| 黄色丝袜av网址大全| 亚洲一区中文字幕在线| 午夜影院日韩av| 色精品久久人妻99蜜桃| 国内精品久久久久久久电影| 午夜精品一区二区三区免费看| 好男人电影高清在线观看| 国产精品久久久久久精品电影| 90打野战视频偷拍视频| 男插女下体视频免费在线播放| 久久草成人影院| 日日夜夜操网爽| 亚洲专区国产一区二区| 午夜福利在线在线| 深夜精品福利| 亚洲aⅴ乱码一区二区在线播放 | 亚洲熟女毛片儿| 天天躁狠狠躁夜夜躁狠狠躁| 99精品欧美一区二区三区四区| 欧美成人性av电影在线观看| 操出白浆在线播放| 一本精品99久久精品77| 久久久久国产精品人妻aⅴ院| 好男人电影高清在线观看| 宅男免费午夜| 久久精品国产亚洲av高清一级| 一进一出抽搐动态| 观看免费一级毛片| 午夜激情av网站| 村上凉子中文字幕在线| 国产一区二区激情短视频| 日本一二三区视频观看| 亚洲av电影不卡..在线观看| 人人妻人人看人人澡| 国产高清视频在线播放一区| 久久亚洲精品不卡| 丝袜美腿诱惑在线| 午夜精品久久久久久毛片777| 国产片内射在线| 日韩国内少妇激情av| 国产精品久久久久久亚洲av鲁大| 免费在线观看完整版高清| 免费搜索国产男女视频| 国产单亲对白刺激| 亚洲午夜理论影院| 老司机靠b影院| 免费在线观看完整版高清| 男女视频在线观看网站免费 | 亚洲 国产 在线| 亚洲国产欧洲综合997久久,| 日本 欧美在线| 成人亚洲精品av一区二区| 国产成+人综合+亚洲专区| 日本免费一区二区三区高清不卡| 美女高潮喷水抽搐中文字幕| a级毛片在线看网站| 色av中文字幕| 亚洲一区二区三区不卡视频| 嫩草影视91久久| 久久精品影院6| 91国产中文字幕| 日韩欧美国产在线观看| 精品福利观看| 中文字幕人成人乱码亚洲影| 国产精品精品国产色婷婷| 可以在线观看的亚洲视频| 91老司机精品| 国产亚洲精品久久久久久毛片| 日本 av在线| 欧美在线一区亚洲| 欧美一级毛片孕妇| 久久久久免费精品人妻一区二区| 欧美+亚洲+日韩+国产| 国产一区二区激情短视频| 国产熟女午夜一区二区三区| 91老司机精品| 亚洲电影在线观看av| 国产精品,欧美在线| 母亲3免费完整高清在线观看| 精品午夜福利视频在线观看一区| 无限看片的www在线观看| 欧美日本亚洲视频在线播放| a级毛片a级免费在线| 18禁黄网站禁片午夜丰满| 一级毛片精品| 久久久久久人人人人人| 变态另类丝袜制服| 午夜精品久久久久久毛片777| 久久久精品大字幕| 妹子高潮喷水视频| 亚洲,欧美精品.| 男插女下体视频免费在线播放| 午夜久久久久精精品| 国产亚洲精品久久久久久毛片| 他把我摸到了高潮在线观看| 午夜福利在线观看吧| 久久久久精品国产欧美久久久| 国产精品1区2区在线观看.| 久久这里只有精品中国| 国产精品 国内视频| 成人18禁在线播放| 亚洲精品一卡2卡三卡4卡5卡| 久久这里只有精品中国| 母亲3免费完整高清在线观看| svipshipincom国产片| 国语自产精品视频在线第100页| 男女下面进入的视频免费午夜| 亚洲国产精品成人综合色| 亚洲男人的天堂狠狠| av免费在线观看网站| 黄色女人牲交| 免费看十八禁软件| 正在播放国产对白刺激| 在线观看www视频免费| 亚洲av成人av| 一进一出抽搐动态| 在线观看舔阴道视频| 欧美日本亚洲视频在线播放| 国产伦人伦偷精品视频| 欧美成人免费av一区二区三区| 欧美激情久久久久久爽电影| 国产不卡一卡二| 1024视频免费在线观看| 男女那种视频在线观看| 人妻夜夜爽99麻豆av| 亚洲av日韩精品久久久久久密| 精品无人区乱码1区二区| 久久久久免费精品人妻一区二区| 亚洲国产精品999在线| 免费搜索国产男女视频| 最新美女视频免费是黄的| 日本成人三级电影网站| 亚洲成人国产一区在线观看| 国产欧美日韩一区二区精品| 人人妻人人看人人澡| 亚洲国产精品sss在线观看| 啪啪无遮挡十八禁网站| 久久热在线av| 亚洲 国产 在线| 黄频高清免费视频| 久久久久精品国产欧美久久久| 精品国产亚洲在线| 国产三级黄色录像| 两人在一起打扑克的视频| 少妇熟女aⅴ在线视频| 丰满的人妻完整版| 久久午夜亚洲精品久久| 欧美高清成人免费视频www| 他把我摸到了高潮在线观看| or卡值多少钱| 熟女少妇亚洲综合色aaa.| 欧美成人午夜精品| 女生性感内裤真人,穿戴方法视频| 国内精品久久久久久久电影| 亚洲va日本ⅴa欧美va伊人久久| 两个人免费观看高清视频| 欧美色视频一区免费| 香蕉久久夜色| 亚洲片人在线观看| 亚洲欧美日韩高清专用| 亚洲av成人一区二区三| 99re在线观看精品视频| 91国产中文字幕| 欧美高清成人免费视频www| 国产一区二区在线观看日韩 | 国产av一区在线观看免费| 12—13女人毛片做爰片一| 欧美日韩亚洲综合一区二区三区_| 美女扒开内裤让男人捅视频| 在线观看www视频免费| 久久久精品欧美日韩精品| 亚洲国产中文字幕在线视频| 老汉色∧v一级毛片| 久久人妻福利社区极品人妻图片| 国产精品 欧美亚洲| 黄色 视频免费看| 国产伦人伦偷精品视频| 91大片在线观看| 午夜免费成人在线视频| 亚洲国产欧美一区二区综合| 国产精品电影一区二区三区| 丁香欧美五月| 两性夫妻黄色片| 国产成+人综合+亚洲专区| 欧美午夜高清在线| 久久久精品国产亚洲av高清涩受| 曰老女人黄片| 亚洲精品国产精品久久久不卡| 人妻久久中文字幕网| 少妇粗大呻吟视频| 色播亚洲综合网| 亚洲av片天天在线观看| 丰满人妻一区二区三区视频av | 欧美不卡视频在线免费观看 | 久久久精品大字幕| av在线播放免费不卡| 久久久精品国产亚洲av高清涩受| 国产精品亚洲一级av第二区| 欧美日韩乱码在线| 99精品在免费线老司机午夜| 国产高清激情床上av| 在线看三级毛片| 妹子高潮喷水视频| 国产成人精品无人区| 男男h啪啪无遮挡| 宅男免费午夜| 欧美一级毛片孕妇| 日本免费一区二区三区高清不卡| av超薄肉色丝袜交足视频| 中文资源天堂在线| 一二三四社区在线视频社区8| 一卡2卡三卡四卡精品乱码亚洲| 伊人久久大香线蕉亚洲五| 黄片大片在线免费观看| 在线十欧美十亚洲十日本专区| 国产一区二区激情短视频| 精品久久久久久久久久久久久| 午夜福利18| 日韩有码中文字幕| 一二三四在线观看免费中文在| 丝袜美腿诱惑在线| 亚洲av中文字字幕乱码综合| 18禁国产床啪视频网站| 成人三级做爰电影| 亚洲真实伦在线观看| 麻豆成人午夜福利视频| 国产精品亚洲一级av第二区| 两个人视频免费观看高清| 黑人操中国人逼视频| 国产av不卡久久| 中文在线观看免费www的网站 | 国产亚洲欧美98| 亚洲专区中文字幕在线| 嫁个100分男人电影在线观看| 欧美激情久久久久久爽电影| 最近视频中文字幕2019在线8| 看免费av毛片| 亚洲成人国产一区在线观看| 亚洲国产中文字幕在线视频| 亚洲欧美日韩无卡精品| 亚洲人成伊人成综合网2020| 12—13女人毛片做爰片一| 桃红色精品国产亚洲av| 国产精品98久久久久久宅男小说| 一夜夜www| 日韩国内少妇激情av| 我要搜黄色片| 成人精品一区二区免费| 天堂√8在线中文| 男女那种视频在线观看| 99久久精品热视频| 午夜福利成人在线免费观看| 日韩欧美精品v在线| 桃色一区二区三区在线观看| 欧美+亚洲+日韩+国产| 三级毛片av免费| 日本免费一区二区三区高清不卡| 制服诱惑二区| 精品国产亚洲在线| 桃红色精品国产亚洲av| 在线a可以看的网站| 国产黄色小视频在线观看| 欧美日韩国产亚洲二区| 国产精品精品国产色婷婷| 亚洲av片天天在线观看| 久久精品综合一区二区三区| 啦啦啦韩国在线观看视频| 亚洲国产欧洲综合997久久,| 九色国产91popny在线| 丰满的人妻完整版| 色尼玛亚洲综合影院| 亚洲片人在线观看| 精品久久久久久成人av| 美女午夜性视频免费| 久久草成人影院| 少妇人妻一区二区三区视频| 亚洲免费av在线视频| av天堂在线播放| 窝窝影院91人妻| 国产黄a三级三级三级人| 50天的宝宝边吃奶边哭怎么回事| 国产精品免费一区二区三区在线| 精品久久久久久久久久久久久| 两个人免费观看高清视频| 日本三级黄在线观看| 国产一区二区在线观看日韩 | 白带黄色成豆腐渣| 亚洲国产精品成人综合色| 免费无遮挡裸体视频| 亚洲国产欧美一区二区综合| 婷婷精品国产亚洲av| 久久天躁狠狠躁夜夜2o2o| 国产精品98久久久久久宅男小说| 老鸭窝网址在线观看| 国产精品98久久久久久宅男小说| 日本 欧美在线| 精品少妇一区二区三区视频日本电影| 一边摸一边抽搐一进一小说| av福利片在线观看| a在线观看视频网站| 在线观看www视频免费| 成年人黄色毛片网站| 亚洲成人久久性| 99久久国产精品久久久| 99国产极品粉嫩在线观看| 夜夜夜夜夜久久久久| 在线视频色国产色| 欧美另类亚洲清纯唯美| 国产蜜桃级精品一区二区三区| 天堂影院成人在线观看| 亚洲精品色激情综合| 精品国内亚洲2022精品成人| 国产精品日韩av在线免费观看| 欧美高清成人免费视频www| 香蕉丝袜av| 色老头精品视频在线观看| 最近最新免费中文字幕在线| 亚洲成人免费电影在线观看| 国产真实乱freesex| 18禁观看日本| 亚洲成a人片在线一区二区| 色在线成人网| 亚洲av电影不卡..在线观看| 99久久精品热视频| 91老司机精品| 亚洲无线在线观看| 国产免费男女视频| av超薄肉色丝袜交足视频|