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

    耦合Frenkel-Kontorova 雙鏈的格波解及其色散關(guān)系*

    2022-08-12 14:27:56蘇瑞霞黃霞鄭志剛
    物理學(xué)報(bào) 2022年15期
    關(guān)鍵詞:雙鏈聲子色散

    蘇瑞霞 黃霞 鄭志剛

    1) (中國礦業(yè)大學(xué)(北京) 理學(xué)院物理系,北京 100083)

    2) (華北電力大學(xué)數(shù)理學(xué)院,北京 102206)

    3) (華僑大學(xué)信息科學(xué)與工程學(xué)院,廈門 361021)

    求解了兩條相同的耦合Frenkel-Kontorova (FK) 鏈在低溫、有限溫和高溫情況下的格波解及色散關(guān)系,進(jìn)而研究了耦合FK 雙鏈的晶格振動(dòng)特點(diǎn).結(jié)果表明,耦合FK 雙鏈的色散關(guān)系包含一個(gè)聲學(xué)支和一個(gè)光學(xué)支,兩者的頻譜范圍和頻率禁帶與FK 鏈的恢復(fù)力系數(shù)、鏈間耦合強(qiáng)度系數(shù)均有關(guān)聯(lián),低溫和有限溫的情況還與外勢深度有關(guān)系.并且研究發(fā)現(xiàn)當(dāng)鏈間耦合強(qiáng)度較小時(shí),不存在頻率禁帶;當(dāng)鏈間耦合強(qiáng)度逐漸增加到某一臨界值后,頻率禁帶出現(xiàn),且隨著鏈間耦合強(qiáng)度增加,頻隙不斷變大,這是因?yàn)楣鈱W(xué)支隨著鏈間耦合強(qiáng)度增加不斷向高頻方向移動(dòng).此外,還發(fā)現(xiàn)帶隙結(jié)構(gòu)出現(xiàn)的臨界鏈間耦合強(qiáng)度始終為FK 鏈恢復(fù)力系數(shù)的2 倍,并不受溫度的影響.本文還研究了給定鏈間耦合強(qiáng)度下溫度對耦合FK 雙鏈色散關(guān)系的影響規(guī)律.本研究內(nèi)容可為分析鏈間界面耦合和溫度對晶格的振動(dòng)特點(diǎn)和物理性質(zhì)的影響提供理論依據(jù),從而對于能量輸運(yùn)、熱調(diào)控等實(shí)際應(yīng)用發(fā)揮重要的指導(dǎo)作用.

    1 引言

    目前中國芯片雖已步入世界領(lǐng)先,但仍存在芯片散熱等卡脖子問題[1?4].而解決這些問題的關(guān)鍵步驟一定會(huì)涉及到固體中的熱傳導(dǎo)問題和高性能導(dǎo)熱材料.對于固體中的熱傳導(dǎo)可以分為電子導(dǎo)熱和晶格導(dǎo)熱.對于非金屬和一般半導(dǎo)體結(jié)構(gòu)中的熱傳導(dǎo)主要是晶格的導(dǎo)熱,其來源于晶格振動(dòng)的格波.格波是固體物理中典型的集體激發(fā),它代表系統(tǒng)內(nèi)所有原子在各自的平衡位置附近振動(dòng)的一種集體模式.在簡諧近似下格波是相互獨(dú)立的,聲子是指格波的量子.一個(gè)格波,即一種振動(dòng)模,稱為一個(gè)聲子模式.它是晶格原子集體運(yùn)動(dòng)狀態(tài)的激發(fā)單元,因此系統(tǒng)一般狀態(tài)可看作是這些獨(dú)立的基本激發(fā)單元(聲子)的集合.引入聲子概念后,N個(gè)原子的耦合振蕩問題在簡諧近似下就被約化為無相互作用的理想聲子氣體來進(jìn)行理論分析和處理,由此可以很好地解釋晶體熱力學(xué)性質(zhì)的微觀機(jī)制.

    最早物理學(xué)家們以一維單原子鏈為固體材料的簡化模型,在簡諧近似下研究了一維單原子鏈的格波理論及色散關(guān)系[5,6].考慮到實(shí)際固體往往是非線性的,需加入非線性作用,而非線性作用的出現(xiàn)使得聲子與聲子相互作用,導(dǎo)致色散關(guān)系在一般情況下很難解析求解.有些情況溫度比較低時(shí),可將系統(tǒng)的非線性項(xiàng)看作是一個(gè)小的微擾,運(yùn)用標(biāo)準(zhǔn)力矩陣擾動(dòng)法可以解析得到色散關(guān)系[7];而溫度比較高時(shí),小的非線性微擾近似已不再成立.近些年來研究學(xué)者們提出了聲子譜能量密度(spectral energy density,SED)方法[8?10]和有效聲子理論(effective phonon theory,EPT)[11?13]來預(yù)測非線性系統(tǒng)的色散關(guān)系.這兩種方法在高溫下也可以準(zhǔn)確地預(yù)測非線性系統(tǒng)的色散關(guān)系.但值得注意的是,上述預(yù)測色散關(guān)系的理論大多適用于獨(dú)立的固體結(jié)構(gòu)和獨(dú)立的晶格鏈模型.但在實(shí)際應(yīng)用中,固體材料往往不是孤立存在的,最常見的為耦合式結(jié)構(gòu),例如廣泛應(yīng)用于芯片散熱、航空航天等領(lǐng)域的碳納米管陣列[14?17]和多層石墨烯[18?22]等納米材料.耦合式的結(jié)構(gòu)會(huì)導(dǎo)致界面的出現(xiàn),這時(shí)除了非線性導(dǎo)致的聲子與聲子相互作用,還會(huì)出現(xiàn)聲子與界面的作用,從而影響到系統(tǒng)的色散關(guān)系.目前關(guān)于這方面的研究還十分匱乏,因此非常有必要研究耦合系統(tǒng)的格波理論及色散關(guān)系.

    本文選取耦合Frenkel-Kontorova (FK) 雙鏈模型來研究耦合系統(tǒng)的格波理論.FK 模型是在1938 年由Frenkel 和Kontorovaza[23]兩位物理學(xué)家在研究晶體位錯(cuò)運(yùn)動(dòng)時(shí)首次提出.此模型形式簡單,是一個(gè)置于外部周期勢場中和最近鄰粒子間存在簡諧相互作用的原子鏈模型.它可以很好地描述廣泛的非線性物理現(xiàn)象,例如電荷密度波的傳播、在晶格表面原子吸附層的動(dòng)力學(xué),以及磁有序結(jié)構(gòu)中的磁疇壁等;在凝聚態(tài)物理和統(tǒng)計(jì)物理研究中也被廣泛地用來模擬系統(tǒng)的熱傳導(dǎo)及聲子輸運(yùn)過程.因此我們采用耦合FK 雙鏈模型.本文主要計(jì)算了在低溫小振動(dòng)近似、有限溫、高溫近似情況下耦合FK 雙鏈的格波解及色散關(guān)系.具體討論了色散關(guān)系曲線的分布情況、頻譜范圍和帶隙結(jié)構(gòu)的影響因素等.并且給出了溫度對耦合FK 雙鏈色散關(guān)系的影響規(guī)律.研究結(jié)果可為分析鏈間界面耦合對晶格的振動(dòng)特點(diǎn)和物理性質(zhì)的影響提供理論依據(jù).

    2 耦合FK 雙鏈的色散關(guān)系

    2.1 耦合FK 雙鏈模型

    考慮兩條相同的耦合FK 鏈模型(如圖1 所示),選取最簡單的鏈間耦合作用,即鏈1 與鏈2 通過簡諧耦合作用相連接,耦合強(qiáng)度為kc.為計(jì)算方便,取兩條FK 鏈的原子質(zhì)量相同均為m.un,j表示鏈n(n=1,2)中第j個(gè)粒子偏離平衡位置的位移.假設(shè)每條鏈中相鄰原子之間的平衡距離為a,兩條鏈當(dāng)中的第j個(gè)粒子對為一個(gè)原胞,所以原胞的大小為a.

    圖1 耦合FK 雙鏈模型示意圖Fig.1.Schematics of the coupled FK chains model.

    2.2 低溫近似下耦合FK 雙鏈的色散關(guān)系

    2.2.1 色散關(guān)系計(jì)算方法

    FK 鏈中相鄰原子間的相互作用力是線性恢復(fù)力,設(shè)恢復(fù)力常數(shù)為kn(n=1,2),外勢Un(un,j)?Vsub,n(2π)?2cos(2πun,j),Vsub,n為鏈n的外勢深度.則耦合FK 雙鏈中第j個(gè)原胞內(nèi)兩個(gè)原子的運(yùn)動(dòng)方程為

    (1)式很難解析求解,在低溫情況下,原子的動(dòng)能較小,其活動(dòng)范圍主要限制在 FK 鏈外勢的勢阱中.因此,在低溫小振動(dòng)近似下,耦合系統(tǒng)相應(yīng)的運(yùn)動(dòng)方程可以寫為

    設(shè)該方程組有行波解:

    其中,a為每條鏈中相鄰原子之間的平衡距離,ω為頻率,q為波數(shù).將(3)式代入(2)式中可以得到關(guān)于振幅A1和A2的方程組:

    方程組(4)式有非零解的條件是其系數(shù)矩陣的行列式等于零,即

    這里令

    由此可得

    這是關(guān)于ω2的一元二次方程,取相同的外勢深度并令Vsub,1Vsub,2V求解(7)式可得

    因?yàn)閮蓷lFK 鏈相同,所以k1k2,可以得到

    2.2.2 色散關(guān)系結(jié)果分析

    (9)式給出了兩條相同的耦合FK 鏈系統(tǒng)在低溫小振動(dòng)近似下的色散關(guān)系.可以發(fā)現(xiàn)聲學(xué)波的色散關(guān)系(ω?)始終與單條FK 鏈結(jié)果[24]完全相同,并不隨鏈間耦合作用而發(fā)生變化;而光學(xué)波的色散關(guān)系(ω+)隨鏈間耦合作用增強(qiáng)向高頻方向移動(dòng).這是因?yàn)閷τ诼晫W(xué)支ω?,第j個(gè)原胞中的粒子對在做同向運(yùn)動(dòng),運(yùn)動(dòng)完全一致,像一個(gè)剛體一樣整體地振動(dòng),因此對應(yīng)的色散關(guān)系與單條FK 鏈結(jié)果完全相同.但值得注意的是這里只是色散關(guān)系表現(xiàn)形式與單鏈相同,仍然是兩條鏈所有原子集體運(yùn)動(dòng)的結(jié)果.對于光學(xué)支ω+,第j個(gè)原胞中的粒子對在做相對振動(dòng),運(yùn)動(dòng)完全相反,運(yùn)動(dòng)中保持它們的質(zhì)心不變,對應(yīng)的光學(xué)支會(huì)受到鏈間耦合強(qiáng)度的調(diào)制.由(9)式的色散關(guān)系可以發(fā)現(xiàn),聲學(xué)支和光學(xué)支的頻譜范圍 Δω?、Δω+和頻率禁帶 Δω′與FK 鏈的恢復(fù)力系數(shù)、外勢深度以及鏈間耦合強(qiáng)度系數(shù)均有關(guān)聯(lián).

    現(xiàn)對(9)式得到的ω+和ω?進(jìn)行數(shù)值擬合,得色散關(guān)系曲線如圖2(k1k21,V1,m1)所示.圖2 呈現(xiàn)了不同鏈間耦合強(qiáng)度kc下的色散關(guān)系,可以看到頻率禁帶的存在與否取決于鏈間耦合強(qiáng)度的大小.當(dāng)鏈間耦合強(qiáng)度較小時(shí),不存在頻率禁帶(見圖2(a)),但光學(xué)支位置隨著鏈間耦合強(qiáng)度的增加不斷明顯升高;當(dāng)鏈間耦合強(qiáng)度增加到某一臨界值后,頻率禁帶將出現(xiàn)(見圖2(b)),并隨著鏈間耦合強(qiáng)度增加,頻隙不斷變大(見圖2(c)).通過以下計(jì)算,可以得到頻隙為0 的臨界鏈間耦合強(qiáng)度.令

    圖2 耦合FK 雙鏈在低溫小振動(dòng)近似情況下的色散關(guān)系(k1 k2 1,V1,m1) (a) kc 1;(b) kc 2 ;(c) kc 5Fig.2.Dispersion relationship of the coupled FK chains under low temperature approximation with small vibration(k1 k2 1,V1,m1):(a) kc 1;(b) kc 2 ;(c) kc 5 .

    由(9)式得

    將(11)式和(12)式代入(10)式計(jì)算可得

    由以上計(jì)算結(jié)果可知,兩條相同耦合FK 鏈系統(tǒng)在低溫下的臨界鏈間耦合強(qiáng)度為FK 鏈恢復(fù)力系數(shù)的2 倍,與外勢深度無關(guān).取k1k21,可得臨界鏈間耦合強(qiáng)度kc2,與圖2(b)數(shù)值計(jì)算結(jié)果一致.

    2.3 有限溫度下耦合FK 雙鏈的色散關(guān)系

    2.3.1 色散關(guān)系計(jì)算方法

    第2.2 節(jié)考慮的是耦合FK 雙鏈模型在低溫小振動(dòng)近似下的情況,外勢泰勒展開的高階項(xiàng)可以忽略,而在有限溫度,即溫度并不是很低的情況下,這樣的近似是不合理的.這時(shí)需要借助耦合系統(tǒng)的EPT[25]來具體討論耦合FK 雙鏈的色散關(guān)系及鏈間耦合對此耦合系統(tǒng)色散關(guān)系的影響.

    首先給出兩條相同的耦合非線性鏈模型的哈密頓量:

    這里H1,H2和Hc分別表示鏈1、鏈2 和耦合項(xiàng)的哈密頓量.un,j仍然表示鏈n(n=1,2)中第j個(gè)粒子偏離平衡位置的位移,pn,j是此粒子的動(dòng)量.N表示鏈長,mn是每條鏈粒子的質(zhì)量(為方便,設(shè)所有的粒子質(zhì)量mn1).為了不失一般性,鏈內(nèi)相鄰粒子相互作用勢Vn(un,j ?un,j?1)、外勢Un(un,j)和鏈間相互作用勢V12(u1,j ?u2,j) 的表達(dá)式為

    其中,gs為鏈內(nèi)相鄰粒子相互作用強(qiáng)度,σs為外勢作用強(qiáng)度,ls為鏈間相互作用強(qiáng)度,s2,3,4,···.根據(jù)耦合非線性雙鏈的EPT[25],其重整化聲子譜的解析表達(dá)式為

    其中,ωn,k,+為耦合系統(tǒng)聲子譜的光學(xué)支,ωn,k,?為聲學(xué)支;ωk2 sin(k/2) 是單條簡諧鏈的聲子譜(這里kaq).重整化的參數(shù)αn,γn,+和γn,?為

    其中,對于鏈1 (n=1) 取“+”,對于鏈2 (n=2)取“–”;

    上述理論給出的是一般的兩條相同耦合非線性鏈模型的有效聲子譜.對于兩條相同的耦合FK鏈模型(鏈間耦合作用仍然通過簡諧耦合作用相連接),表達(dá)式為

    在有限溫度下,粒子在外勢中振動(dòng)還不太劇烈,振動(dòng)較小,可以對外勢Un(un,j) 進(jìn)行泰勒展開,從而得到耦合FK 雙鏈模型的重整化的參數(shù)αn,γn,?和γn,+為

    其中,對于鏈1 (n=1) 取“+”,對于鏈2 (n=2)取“–”.

    2.3.2 色散關(guān)系結(jié)果分析

    對2.3.1 節(jié)的EPT 計(jì)算結(jié)果進(jìn)行數(shù)值擬合,溫度T=0.05 時(shí)的耦合FK 雙鏈(k1k21,V1,m1)的色散關(guān)系如圖3 所示,可以看出色散關(guān)系和圖2 相比有所下移.這是由于有限溫度下,外勢的高階項(xiàng)作用開始明顯,第3 節(jié)會(huì)進(jìn)一步說明溫度對色散關(guān)系的影響.同樣圖3 也展示了色散關(guān)系隨鏈間耦合作用的變化情況:當(dāng)鏈間耦合強(qiáng)度較小時(shí),不存在頻率禁帶(見圖3(a));當(dāng)鏈間耦合強(qiáng)度增加到某一臨界值后,頻率禁帶將出現(xiàn)(見圖3(b)),并隨著鏈間耦合強(qiáng)度增加,頻隙不斷變大(見圖3(c)).通過以下計(jì)算,可得到頻隙為0的臨界鏈間耦合強(qiáng)度.令

    圖3 耦合FK 雙鏈(k1 k2 1,V1,m1)在有限溫度(T=0.05)下的色散關(guān)系 (a) kc 1;(b) kc 2 ;(c) kc 5Fig.3.Dispersion relationship of the coupled FK chains(k1 k2 1,V1,m1) under finite temperature(T=0.05):(a) kc 1;(b) kc 2;(c) kc 5 .

    由(18)式結(jié)合(25)式—(27)式得

    將(29)式和(30)式代入(28)式計(jì)算可得

    上述計(jì)算結(jié)果說明對于兩條相同的耦合FK鏈系統(tǒng),在有限溫度下,臨界鏈間耦合強(qiáng)度仍為FK 鏈恢復(fù)力系數(shù)的2 倍(圖3(b)數(shù)值結(jié)果也驗(yàn)證了這一結(jié)論),與外勢深度無關(guān),和低溫小振動(dòng)近似下的結(jié)論相同.

    2.4 高溫近似下耦合FK 雙鏈的色散關(guān)系

    2.4.1 色散關(guān)系計(jì)算方法

    在高溫區(qū)域,由于粒子動(dòng)能很大,會(huì)脫離外勢的束縛,這時(shí)外勢可以忽略,所以FK 鏈又回歸到了簡諧鏈,兩條相同耦合FK 鏈的色散關(guān)系和兩條相耦合簡諧鏈的結(jié)果相同(只需令(9)式的外勢深度V=0),即

    2.4.2 色散關(guān)系結(jié)果分析

    由(32)式可知,高溫情況下頻譜范圍和頻率禁帶僅與FK 鏈的恢復(fù)力系數(shù)、鏈間耦合強(qiáng)度系數(shù)有關(guān)聯(lián),與外勢深度無關(guān),這點(diǎn)與低溫及有限溫的情況有所不同.圖4 給出了高溫下耦合FK 雙鏈(k1k21,V1,m1)的色散關(guān)系.發(fā)現(xiàn)在長波極限下,聲學(xué)支頻率為0,光學(xué)支也比低溫和有限溫情況有明顯下移,主要由于高溫下粒子運(yùn)動(dòng)劇烈,脫離了外勢的約束.此外,從圖4 可以看出,當(dāng)鏈間耦合強(qiáng)度逐漸增加到某一臨界值后,頻率禁帶出現(xiàn),并隨著鏈間耦合強(qiáng)度增加,頻隙不斷變大,這在低溫、有限溫的討論中已提到,結(jié)論均類似.按照前面類似的計(jì)算方法,也可得到頻隙為0 的臨界鏈間耦合強(qiáng)度為kc2k1.結(jié)合低溫和有限溫的討論可知:兩條相同的耦合FK 鏈系統(tǒng),在鏈間耦合作用為簡諧耦合作用時(shí),臨界鏈間耦合強(qiáng)度始終為FK 鏈的恢復(fù)力系數(shù)的2 倍,即kc2k1,并不受溫度的影響.

    圖4 耦合FK 雙鏈在高溫近似下的色散關(guān)系(k1 k2 1,V 1,m1) (a) kc 1;(b) kc 2;(c) kc 5Fig.4.Dispersion relation of coupled FK chains under high temperature approximation (k1 k2 1,V1,m1):(a) kc 1;(b) kc 2;(c) kc 5 .

    上述結(jié)論表明FK 鏈的恢復(fù)力系數(shù)越大,所對應(yīng)臨界鏈間耦合強(qiáng)度越大.換言之,要想使恢復(fù)力系數(shù)很大的耦合雙鏈出現(xiàn)頻隙結(jié)構(gòu),鏈間耦合強(qiáng)度必須足夠強(qiáng).鏈間耦合強(qiáng)度的強(qiáng)弱和很多因素有關(guān),如果耦合固體材料缺乏良好的接觸或表面不是很平整,界面耦合強(qiáng)度系數(shù)就會(huì)比較弱;如果耦合固體材料能夠彼此很好地接觸,在界面處就可能形成化學(xué)鍵,導(dǎo)致界面耦合強(qiáng)度系數(shù)會(huì)很強(qiáng).從上述結(jié)論發(fā)現(xiàn),通過改變耦合雙鏈的恢復(fù)力系數(shù),進(jìn)而能改變帶隙結(jié)構(gòu).對帶隙結(jié)構(gòu)的改變和控制可以指導(dǎo)人們設(shè)計(jì)一定的聲子禁帶和允帶(處于禁帶頻率范圍內(nèi)的振動(dòng)將被禁止在晶體中傳播[26,27];處于通帶頻率范圍內(nèi)的振動(dòng)會(huì)傳遍整個(gè)結(jié)構(gòu)),可用于模式選擇、能量傳輸、傳播的指向性等諸多領(lǐng)域,對新型功能材料在工程上的應(yīng)用很有幫助[28,29].

    3 溫度對耦合FK 雙鏈色散關(guān)系的影響

    第2 節(jié)給出了在低溫、有限溫和高溫情況下耦合FK 雙鏈色散關(guān)系的解析表達(dá)式,并分析了色散關(guān)系的分布情況、頻譜范圍和帶隙結(jié)構(gòu)的影響因素等.本節(jié)對于相同的耦合FK 雙鏈,固定鏈間耦合強(qiáng)度kc,具體討論溫度對耦合FK 雙鏈色散關(guān)系的影響規(guī)律.但什么溫度為低溫? 什么溫度為高溫?第2.2 節(jié)和第2.4 節(jié)的討論只是表明在溫度足夠低和足夠高時(shí),可分別用(9)式和(32)式計(jì)算色散關(guān)系,但并不明確溫度具體要低于或高于什么值.考慮到EPT 可適用于低溫情況,我們可以用EPT 來計(jì)算具體低溫?cái)?shù)值下的色散關(guān)系,并與低溫小振動(dòng)近似下的結(jié)果對比.而在高溫時(shí),采用SED 方法[8?10]計(jì)算具體高溫?cái)?shù)值下的色散關(guān)系,并與2.4 節(jié)高溫解析的結(jié)果對比.下面簡單介紹SED 方法[8?10],SED 的表達(dá)式為

    圖5 展示了耦合FK 雙鏈在kc1.0 時(shí)的色散關(guān)系隨溫度的變化情況(考慮到SED 的計(jì)算時(shí)間成本,本文只計(jì)算了波數(shù)0— π/a(鏈常數(shù)a=1)范圍的色散關(guān)系,但數(shù)據(jù)足夠說明問題).圖5 中紅實(shí)線和黑實(shí)線是EPT 的計(jì)算結(jié)果;紅色和黑色圓點(diǎn)是在低溫小振動(dòng)近似下計(jì)算的結(jié)果;紅色和黑色方框是高溫下計(jì)算的解析結(jié)果;等高線圖是SED 方法的計(jì)算結(jié)果.從圖5(a)和圖5(b)可以發(fā)現(xiàn),溫度T低于0.01 時(shí),EPT 計(jì)算的色散關(guān)系與低溫小振動(dòng)下的解析結(jié)果幾乎無差別,完全吻合.這說明溫度低于0.01 時(shí),就可以用前面低溫小振動(dòng)近似得到的解析結(jié)果去計(jì)算色散關(guān)系,并且色散關(guān)系分布在此低溫區(qū)域基本不受溫度影響.當(dāng)溫度升高,如溫度高于0.01 時(shí)(見圖5(c)和圖5(d)),色散關(guān)系與圖5(a)和圖5(b)相比有明顯下移,且溫度越高,下移幅度越大.這是由于溫度升高,外勢高階項(xiàng)開始起作用,并且發(fā)現(xiàn)聲學(xué)支下移幅度較光學(xué)支下移幅度大.因此可以推論在鏈間耦合強(qiáng)度足夠大,出現(xiàn)頻隙結(jié)構(gòu)的耦合系統(tǒng)中,溫度升高頻隙會(huì)變大.當(dāng)溫度進(jìn)一步升高,例如溫度高于0.1 時(shí)(見圖5(e)和圖5(f)),色散關(guān)系較圖5(c)和圖5(d)有大幅度下移,聲學(xué)支在長波極限下頻率幾乎降為0,光學(xué)支也有明顯下移.說明在溫度高于0.1 時(shí),外勢對粒子運(yùn)動(dòng)影響已經(jīng)可以忽略,并且SED 的計(jì)算結(jié)果和前面高溫解析計(jì)算的色散關(guān)系完全吻合.因此溫度高于0.1 時(shí),就可以用2.4 節(jié)高溫解析結(jié)果來計(jì)算色散關(guān)系,并且色散關(guān)系分布在高溫區(qū)域基本不受溫度影響.

    圖5 耦 合FK 雙鏈在不同溫 度下的色散關(guān) 系(k1 k2 1,V1,m1 ,kc 1.0) (a) T0.005;(b) T0.01 ;(c) T0.02;(d) T0.06;(e) T0.10;(f) T 0.20Fig.5.Dispersion relation of coupled FK chains at different temperature (k1 k2 1,V1,m1 ,kc 1.0):(a) T0.005 ;(b) T0.01;(c) T0.02;(d) T0.06;(e) T0.10;(f) T0.20 .

    4 總結(jié)

    本文計(jì)算了兩條相同的耦合FK 鏈系統(tǒng)在低溫、有限溫和高溫情況下色散關(guān)系的解析表達(dá)式.研究發(fā)現(xiàn)頻譜范圍、頻率禁帶與FK 鏈的恢復(fù)力系數(shù)、鏈間耦合強(qiáng)度系數(shù)均有關(guān)聯(lián),低溫和有限溫的情況還與外勢深度有關(guān)系.帶隙結(jié)構(gòu)出現(xiàn)的臨界鏈間耦合強(qiáng)度始終為FK 鏈恢復(fù)力系數(shù)的2 倍,并且不受溫度的影響.進(jìn)一步研究了給定鏈間耦合強(qiáng)度下溫度對耦合FK 雙鏈色散關(guān)系的影響規(guī)律.發(fā)現(xiàn)溫度在低于0.01 時(shí),色散關(guān)系分布基本不受溫度影響;溫度在高于0.01 時(shí),色散關(guān)系有明顯下移,且溫度越高,下移幅度越大,直至溫度上升到0.1,聲學(xué)支在長波極限下頻率幾乎降為0.至此耦合系統(tǒng)脫離外勢的影響,溫度再升高色散關(guān)系分布將不受溫度影響.此外還給出了低溫小振動(dòng)近似及高溫下解析得到的色散關(guān)系的溫度適用范圍.上述研究結(jié)果讓我們對鏈間耦合作用和溫度使晶格振動(dòng)色散關(guān)系產(chǎn)生的變化有了具體的認(rèn)識(shí),有助于拓寬和加深對耦合材料物理性質(zhì)的了解,可以指導(dǎo)實(shí)現(xiàn)設(shè)計(jì)一定的聲子禁帶和允帶,可用于模式選擇、能量傳輸、傳播的指向性等諸多領(lǐng)域,這對新型功能材料在工程上的應(yīng)用很有幫助[28,29].此外,鏈間耦合作用的存在及溫度的變化導(dǎo)致系統(tǒng)色散關(guān)系的改變可以幫助我們理解很多傳熱現(xiàn)象.例如碳納米管陣列和多層石墨烯等納米材料的熱導(dǎo)率明顯低于單臂碳納米管[30]和單層石墨烯[31]的熱導(dǎo)率;溫度升高使得納米材料熱導(dǎo)率降低等.因此本文研究的結(jié)果對于能量輸運(yùn)、熱調(diào)控等實(shí)際應(yīng)用有著重要的指導(dǎo)意義.

    猜你喜歡
    雙鏈聲子色散
    “光的折射”“光的色散”知識(shí)鞏固
    “光的折射”“光的色散”知識(shí)鞏固
    半無限板類聲子晶體帶隙仿真的PWE/NS-FEM方法
    “光的折射”“光的色散”知識(shí)鞏固
    納米表面聲子 首次實(shí)現(xiàn)三維成像
    聲子晶體覆蓋層吸聲機(jī)理研究
    昆蟲共生細(xì)菌活體制造雙鏈RNA
    海外星云 (2021年21期)2021-01-19 14:17:31
    『光的折射』『光的色散』隨堂練
    基于聲子晶體理論的導(dǎo)線防舞方法及數(shù)值驗(yàn)證
    高新區(qū)科技企業(yè)孵化網(wǎng)絡(luò)“雙層雙鏈”結(jié)構(gòu)研究
    中文字幕另类日韩欧美亚洲嫩草| 黄片播放在线免费| 99久久99久久久精品蜜桃| 国产精品久久久久久精品古装| 色在线成人网| 亚洲熟妇中文字幕五十中出 | 无人区码免费观看不卡| 精品第一国产精品| 亚洲欧美精品综合一区二区三区| 丝瓜视频免费看黄片| 少妇被粗大的猛进出69影院| 啪啪无遮挡十八禁网站| 大片电影免费在线观看免费| 曰老女人黄片| 国产在线观看jvid| 国产成人精品无人区| 男人的好看免费观看在线视频 | 久久久国产精品麻豆| 欧美av亚洲av综合av国产av| 国产99白浆流出| 欧美人与性动交α欧美软件| 高潮久久久久久久久久久不卡| 免费日韩欧美在线观看| 一级毛片高清免费大全| 真人做人爱边吃奶动态| 好看av亚洲va欧美ⅴa在| 9热在线视频观看99| 久久久精品区二区三区| 精品电影一区二区在线| 一级a爱片免费观看的视频| 热99久久久久精品小说推荐| 日韩三级视频一区二区三区| 狠狠婷婷综合久久久久久88av| 亚洲七黄色美女视频| a级毛片在线看网站| 别揉我奶头~嗯~啊~动态视频| 高清毛片免费观看视频网站 | 一级片'在线观看视频| 在线观看舔阴道视频| 天天躁夜夜躁狠狠躁躁| 天天添夜夜摸| 亚洲成a人片在线一区二区| 脱女人内裤的视频| 亚洲国产毛片av蜜桃av| 亚洲精品中文字幕一二三四区| av片东京热男人的天堂| 欧美人与性动交α欧美软件| 80岁老熟妇乱子伦牲交| 大香蕉久久成人网| 成人18禁高潮啪啪吃奶动态图| 水蜜桃什么品种好| 99re在线观看精品视频| 国产成人影院久久av| 国产精品 欧美亚洲| 手机成人av网站| 午夜福利欧美成人| 亚洲av日韩在线播放| 久久精品国产99精品国产亚洲性色 | 国产xxxxx性猛交| 日本a在线网址| 首页视频小说图片口味搜索| 亚洲精品一卡2卡三卡4卡5卡| 久久狼人影院| 亚洲五月婷婷丁香| 国产99久久九九免费精品| 亚洲精品中文字幕在线视频| 亚洲精品粉嫩美女一区| 两人在一起打扑克的视频| 伦理电影免费视频| av片东京热男人的天堂| 搡老岳熟女国产| 无遮挡黄片免费观看| a级片在线免费高清观看视频| 少妇被粗大的猛进出69影院| 村上凉子中文字幕在线| 国产成人系列免费观看| 欧美日韩瑟瑟在线播放| 人人澡人人妻人| 色综合婷婷激情| 一级黄色大片毛片| 精品亚洲成国产av| a级片在线免费高清观看视频| 免费观看人在逋| 久久性视频一级片| 成人国语在线视频| 久久人妻av系列| 亚洲av成人一区二区三| 午夜两性在线视频| 男女午夜视频在线观看| 黄片小视频在线播放| 极品教师在线免费播放| 国产精品成人在线| 777米奇影视久久| 超碰成人久久| 精品国产美女av久久久久小说| 久久99一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 免费日韩欧美在线观看| 日韩欧美一区视频在线观看| 后天国语完整版免费观看| 国产成人av激情在线播放| 老汉色∧v一级毛片| 老汉色av国产亚洲站长工具| 久久精品国产99精品国产亚洲性色 | 亚洲视频免费观看视频| 中文欧美无线码| 亚洲精品久久午夜乱码| 亚洲精品在线观看二区| 国产成人精品无人区| 国产片内射在线| 亚洲三区欧美一区| 国产精品二区激情视频| 国产极品粉嫩免费观看在线| 久久精品亚洲熟妇少妇任你| 宅男免费午夜| 精品国产一区二区久久| 国产在视频线精品| 一区福利在线观看| 青草久久国产| 女人爽到高潮嗷嗷叫在线视频| 亚洲成人免费av在线播放| 男女之事视频高清在线观看| 国产aⅴ精品一区二区三区波| 精品国产亚洲在线| 在线天堂中文资源库| 国产国语露脸激情在线看| 正在播放国产对白刺激| 亚洲五月色婷婷综合| 99热只有精品国产| 日韩一卡2卡3卡4卡2021年| 一本大道久久a久久精品| 亚洲国产看品久久| 久久天躁狠狠躁夜夜2o2o| 夜夜躁狠狠躁天天躁| 亚洲全国av大片| 99精国产麻豆久久婷婷| 久久久国产一区二区| 国产精品1区2区在线观看. | 欧美精品av麻豆av| 国产深夜福利视频在线观看| 免费女性裸体啪啪无遮挡网站| 婷婷丁香在线五月| 三上悠亚av全集在线观看| 精品国产乱码久久久久久男人| 国产成人啪精品午夜网站| 在线观看免费视频日本深夜| 国产免费现黄频在线看| 欧美国产精品一级二级三级| 好男人电影高清在线观看| 丝袜美腿诱惑在线| 午夜免费鲁丝| 国产不卡一卡二| 日韩大码丰满熟妇| 老司机深夜福利视频在线观看| 日韩欧美一区视频在线观看| 久久久精品国产亚洲av高清涩受| 老司机亚洲免费影院| 黄色女人牲交| 两人在一起打扑克的视频| 午夜激情av网站| 亚洲av成人av| 免费不卡黄色视频| 一区二区三区激情视频| 亚洲九九香蕉| 9色porny在线观看| 久久狼人影院| 久久这里只有精品19| 亚洲午夜精品一区,二区,三区| 亚洲成人国产一区在线观看| 国产成人免费观看mmmm| 又紧又爽又黄一区二区| 极品教师在线免费播放| 亚洲在线自拍视频| 日本欧美视频一区| 久久香蕉激情| 国产免费男女视频| 人人妻人人爽人人添夜夜欢视频| 国产区一区二久久| 免费在线观看影片大全网站| 少妇裸体淫交视频免费看高清 | 一级毛片高清免费大全| 亚洲精品一二三| 免费在线观看黄色视频的| 欧美乱色亚洲激情| 久久香蕉国产精品| av网站在线播放免费| 国产一区二区三区视频了| 欧美日韩亚洲综合一区二区三区_| 国产激情久久老熟女| 国产一区二区三区视频了| 久久中文字幕一级| 成人影院久久| 日韩成人在线观看一区二区三区| 成人黄色视频免费在线看| 精品久久久精品久久久| 国产精品二区激情视频| 亚洲国产中文字幕在线视频| 韩国精品一区二区三区| 欧美日韩成人在线一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 一二三四社区在线视频社区8| 日本vs欧美在线观看视频| 国产乱人伦免费视频| 久久香蕉国产精品| 成人18禁高潮啪啪吃奶动态图| av中文乱码字幕在线| 老鸭窝网址在线观看| 亚洲色图综合在线观看| av线在线观看网站| 视频在线观看一区二区三区| 国产成人免费无遮挡视频| 国产xxxxx性猛交| 欧美日韩视频精品一区| 国产精品乱码一区二三区的特点 | 在线永久观看黄色视频| 欧美 日韩 精品 国产| 欧美日韩亚洲综合一区二区三区_| 欧美大码av| av天堂在线播放| 午夜福利在线观看吧| 少妇猛男粗大的猛烈进出视频| 精品久久久久久,| 久久九九热精品免费| 高清欧美精品videossex| 黄片小视频在线播放| 18禁裸乳无遮挡动漫免费视频| 久久香蕉精品热| 欧美激情高清一区二区三区| 人妻丰满熟妇av一区二区三区 | 美女 人体艺术 gogo| 亚洲欧美精品综合一区二区三区| 51午夜福利影视在线观看| 精品久久久久久久久久免费视频 | 久久精品国产亚洲av香蕉五月 | 97人妻天天添夜夜摸| 国产激情久久老熟女| 欧美国产精品va在线观看不卡| 日韩 欧美 亚洲 中文字幕| 免费观看人在逋| 正在播放国产对白刺激| 午夜视频精品福利| 男女之事视频高清在线观看| 成人18禁高潮啪啪吃奶动态图| 每晚都被弄得嗷嗷叫到高潮| 岛国毛片在线播放| av网站在线播放免费| 麻豆乱淫一区二区| 亚洲精品美女久久av网站| 在线观看免费午夜福利视频| 搡老熟女国产l中国老女人| 午夜福利视频在线观看免费| 男女之事视频高清在线观看| 精品第一国产精品| 国产在线一区二区三区精| 大型av网站在线播放| 精品一区二区三区av网在线观看| 纯流量卡能插随身wifi吗| 国产一区二区激情短视频| 99久久国产精品久久久| 亚洲国产精品合色在线| 国产真人三级小视频在线观看| 这个男人来自地球电影免费观看| 免费观看精品视频网站| 天天躁夜夜躁狠狠躁躁| 狠狠婷婷综合久久久久久88av| 高清在线国产一区| 亚洲 国产 在线| 久久中文字幕人妻熟女| 国产黄色免费在线视频| 多毛熟女@视频| 欧美精品啪啪一区二区三区| 国产精品香港三级国产av潘金莲| 久久精品亚洲av国产电影网| 窝窝影院91人妻| 99国产极品粉嫩在线观看| 国产精品av久久久久免费| 亚洲成国产人片在线观看| 老司机在亚洲福利影院| 日韩欧美一区视频在线观看| 午夜精品在线福利| 窝窝影院91人妻| 久久久国产欧美日韩av| 亚洲一码二码三码区别大吗| 亚洲欧洲精品一区二区精品久久久| 日本黄色视频三级网站网址 | 色精品久久人妻99蜜桃| 日韩欧美在线二视频 | 一区二区三区精品91| 免费在线观看视频国产中文字幕亚洲| 国产精品.久久久| av天堂久久9| 天天影视国产精品| 性少妇av在线| 日本五十路高清| 一级,二级,三级黄色视频| 不卡av一区二区三区| 一个人免费在线观看的高清视频| 91精品国产国语对白视频| 男女床上黄色一级片免费看| 中文字幕最新亚洲高清| 1024香蕉在线观看| 黑人猛操日本美女一级片| 天天影视国产精品| 久久久久国产一级毛片高清牌| 99久久人妻综合| 国产熟女午夜一区二区三区| 国产亚洲精品久久久久久毛片 | 女人被躁到高潮嗷嗷叫费观| 免费观看a级毛片全部| 成人特级黄色片久久久久久久| 夜夜爽天天搞| 成年人午夜在线观看视频| 国产精品av久久久久免费| 在线观看免费高清a一片| 欧美黄色片欧美黄色片| 99久久综合精品五月天人人| 久久国产精品男人的天堂亚洲| 精品一区二区三区四区五区乱码| 18禁观看日本| 一级作爱视频免费观看| 久久久久久人人人人人| 后天国语完整版免费观看| 国产亚洲精品久久久久久毛片 | 精品电影一区二区在线| 国产精品成人在线| 免费在线观看视频国产中文字幕亚洲| 亚洲av美国av| 很黄的视频免费| 久久国产精品男人的天堂亚洲| 18禁美女被吸乳视频| 色尼玛亚洲综合影院| 老司机亚洲免费影院| 黄片小视频在线播放| 在线观看免费午夜福利视频| 9色porny在线观看| 国产精品免费视频内射| 亚洲av日韩精品久久久久久密| 18禁裸乳无遮挡免费网站照片 | 新久久久久国产一级毛片| 国产午夜精品久久久久久| 日韩欧美在线二视频 | 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品在线美女| 交换朋友夫妻互换小说| 国产精品欧美亚洲77777| 伊人久久大香线蕉亚洲五| 亚洲国产精品一区二区三区在线| 国产精品电影一区二区三区 | 欧美成狂野欧美在线观看| av中文乱码字幕在线| 51午夜福利影视在线观看| 中文字幕最新亚洲高清| 欧美亚洲日本最大视频资源| 精品乱码久久久久久99久播| 波多野结衣一区麻豆| 热99re8久久精品国产| 欧美最黄视频在线播放免费 | cao死你这个sao货| 露出奶头的视频| 伦理电影免费视频| 久久婷婷成人综合色麻豆| 男女之事视频高清在线观看| 久久影院123| 男人舔女人的私密视频| 在线看a的网站| 久久久国产一区二区| 国产精品 国内视频| 午夜老司机福利片| 久久精品国产综合久久久| 一级a爱片免费观看的视频| 久久狼人影院| 国产蜜桃级精品一区二区三区 | 免费日韩欧美在线观看| 两个人免费观看高清视频| 在线av久久热| 变态另类成人亚洲欧美熟女 | 国产区一区二久久| 亚洲av熟女| 无人区码免费观看不卡| 国产亚洲精品一区二区www | 国产又色又爽无遮挡免费看| 十八禁网站免费在线| 正在播放国产对白刺激| 精品一品国产午夜福利视频| 日本vs欧美在线观看视频| 国产精品乱码一区二三区的特点 | 9191精品国产免费久久| 国产精品影院久久| 男人舔女人的私密视频| 啪啪无遮挡十八禁网站| 美女国产高潮福利片在线看| 国产精品久久视频播放| 亚洲熟女毛片儿| 久久久久国内视频| 精品国产亚洲在线| 欧美日韩福利视频一区二区| 极品少妇高潮喷水抽搐| 欧美一级毛片孕妇| 亚洲精品久久成人aⅴ小说| 熟女少妇亚洲综合色aaa.| 亚洲国产看品久久| 国产高清国产精品国产三级| 久久久久精品人妻al黑| 亚洲av成人一区二区三| 免费观看人在逋| 老司机亚洲免费影院| 叶爱在线成人免费视频播放| 色播在线永久视频| 欧美日韩乱码在线| 不卡av一区二区三区| 交换朋友夫妻互换小说| 日韩中文字幕欧美一区二区| 成人黄色视频免费在线看| 精品电影一区二区在线| 国产区一区二久久| 日韩精品免费视频一区二区三区| 黑人猛操日本美女一级片| 亚洲成人国产一区在线观看| 啦啦啦在线免费观看视频4| 91成年电影在线观看| 色综合婷婷激情| 香蕉丝袜av| 99国产精品一区二区三区| 狠狠婷婷综合久久久久久88av| 久久久久国产精品人妻aⅴ院 | 最新美女视频免费是黄的| 日本欧美视频一区| tocl精华| 在线播放国产精品三级| 国产成人欧美| 亚洲av电影在线进入| 男女高潮啪啪啪动态图| 91字幕亚洲| 午夜精品久久久久久毛片777| 精品少妇一区二区三区视频日本电影| 中文字幕另类日韩欧美亚洲嫩草| 天天影视国产精品| 亚洲精品美女久久久久99蜜臀| 色综合欧美亚洲国产小说| 亚洲 国产 在线| 男人操女人黄网站| 激情视频va一区二区三区| 久久久久久久精品吃奶| 一级a爱视频在线免费观看| 国产精品成人在线| 久久精品国产a三级三级三级| 久久久久视频综合| avwww免费| 亚洲精品国产色婷婷电影| 亚洲av成人av| 国产在视频线精品| 视频在线观看一区二区三区| 中文字幕制服av| 欧美丝袜亚洲另类 | 91国产中文字幕| 丝袜美足系列| 久久中文看片网| 美女扒开内裤让男人捅视频| tube8黄色片| 免费不卡黄色视频| 十八禁高潮呻吟视频| 别揉我奶头~嗯~啊~动态视频| 日本黄色视频三级网站网址 | 欧美日韩精品网址| 18禁美女被吸乳视频| 欧美乱码精品一区二区三区| 亚洲中文av在线| 99精品在免费线老司机午夜| 99re6热这里在线精品视频| av天堂久久9| 国产免费av片在线观看野外av| 色婷婷av一区二区三区视频| 自线自在国产av| 日韩有码中文字幕| 国产精品综合久久久久久久免费 | 午夜福利免费观看在线| 久久ye,这里只有精品| 色在线成人网| 18禁黄网站禁片午夜丰满| 少妇粗大呻吟视频| 国产精品偷伦视频观看了| 午夜免费鲁丝| 久久中文看片网| 色尼玛亚洲综合影院| 亚洲精品久久成人aⅴ小说| 91老司机精品| 大陆偷拍与自拍| 欧美精品亚洲一区二区| 亚洲熟女精品中文字幕| videos熟女内射| 国内久久婷婷六月综合欲色啪| 日本黄色日本黄色录像| 成人永久免费在线观看视频| 超色免费av| 国产精品影院久久| 黑人猛操日本美女一级片| 黄色视频不卡| 又黄又粗又硬又大视频| 国产精品成人在线| 色在线成人网| 9色porny在线观看| 欧美激情高清一区二区三区| 亚洲三区欧美一区| 亚洲成a人片在线一区二区| 黄色女人牲交| 看黄色毛片网站| 夜夜爽天天搞| 自线自在国产av| 国产精品久久电影中文字幕 | 色尼玛亚洲综合影院| 亚洲精品粉嫩美女一区| 成年人午夜在线观看视频| 后天国语完整版免费观看| 午夜福利视频在线观看免费| 最近最新中文字幕大全电影3 | 丁香六月欧美| 亚洲av欧美aⅴ国产| 国产精品综合久久久久久久免费 | av超薄肉色丝袜交足视频| 欧美亚洲日本最大视频资源| 久久久久国产一级毛片高清牌| 欧美黑人欧美精品刺激| 99国产精品一区二区三区| 在线观看免费日韩欧美大片| 久久久国产欧美日韩av| 人人澡人人妻人| 欧美日韩乱码在线| 亚洲精品中文字幕在线视频| 777米奇影视久久| 亚洲av成人av| 又黄又爽又免费观看的视频| 波多野结衣一区麻豆| 啦啦啦免费观看视频1| xxxhd国产人妻xxx| 亚洲av日韩在线播放| 丝袜在线中文字幕| 欧美另类亚洲清纯唯美| 国产精品.久久久| 午夜激情av网站| 热99re8久久精品国产| 久久精品国产亚洲av香蕉五月 | 一级片免费观看大全| 国产高清视频在线播放一区| 亚洲精品国产一区二区精华液| 亚洲精品国产精品久久久不卡| 9热在线视频观看99| 一级作爱视频免费观看| 亚洲片人在线观看| 91精品国产国语对白视频| 国产精品 国内视频| 99久久人妻综合| 欧美大码av| 日韩欧美免费精品| 美女视频免费永久观看网站| 亚洲avbb在线观看| 90打野战视频偷拍视频| 国产精品免费视频内射| av片东京热男人的天堂| 99精品欧美一区二区三区四区| 热99国产精品久久久久久7| 一级作爱视频免费观看| 精品人妻在线不人妻| 亚洲午夜理论影院| 久久性视频一级片| av中文乱码字幕在线| 亚洲情色 制服丝袜| 亚洲欧美一区二区三区黑人| 国产精品综合久久久久久久免费 | 亚洲人成电影观看| av线在线观看网站| 啦啦啦视频在线资源免费观看| 69av精品久久久久久| 精品国产超薄肉色丝袜足j| 国产97色在线日韩免费| 老汉色av国产亚洲站长工具| 亚洲色图av天堂| 性少妇av在线| 亚洲成av片中文字幕在线观看| 岛国在线观看网站| 美女 人体艺术 gogo| 国产精品国产av在线观看| 久久精品亚洲精品国产色婷小说| 免费观看a级毛片全部| 精品人妻在线不人妻| 天天躁日日躁夜夜躁夜夜| 如日韩欧美国产精品一区二区三区| 午夜福利欧美成人| 午夜日韩欧美国产| a在线观看视频网站| 国产精品久久久av美女十八| 啦啦啦视频在线资源免费观看| 这个男人来自地球电影免费观看| 波多野结衣av一区二区av| 窝窝影院91人妻| 这个男人来自地球电影免费观看| 亚洲国产毛片av蜜桃av| 少妇猛男粗大的猛烈进出视频| 麻豆av在线久日| 亚洲国产毛片av蜜桃av| 久久国产亚洲av麻豆专区| 亚洲人成电影观看| 日韩 欧美 亚洲 中文字幕| 亚洲国产欧美网| 在线观看免费日韩欧美大片| 一a级毛片在线观看| 精品国产国语对白av| 男人操女人黄网站| 黄色女人牲交| 欧美人与性动交α欧美软件| 国产精品美女特级片免费视频播放器 | 久久热在线av| 精品免费久久久久久久清纯 | 国产精品久久电影中文字幕 | 在线观看免费高清a一片| 一区二区三区精品91|