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

    乙腈在石墨烯極板表面的分子動力學(xué)模擬*

    2022-04-12 03:15:36姚眉捷王永朋梁斐劉書樂
    關(guān)鍵詞:極板擴散系數(shù)電荷

    姚眉捷,王永朋,梁斐,劉書樂

    中山大學(xué)材料科學(xué)與工程學(xué)院,廣東 廣州 510275

    超級電容器,又稱為電化學(xué)電容器,是一種主要依靠雙電層結(jié)構(gòu)或氧化還原贗電容電荷來儲存電能的新型儲能裝置,它的性能介于物理電容器和二次電池之間,兼有物理電容器的功率密度大和二次電池的能量密度高的特點[1-10]。按照反應(yīng)機理來劃分主要分為雙電層電容器和法拉第贗電容電容器,其結(jié)構(gòu)主要由電極材料、隔膜和電解質(zhì)組成[9,11]。對于雙電層電容器,它一般通過活性電極表面對電解液中粒子的可逆性物理吸附來儲存電荷,進而實現(xiàn)電能的儲存。其中,離子液體是常用的電解質(zhì)[4,12-15]。

    離子液體在服役工況溫度下具有獨特的物理化學(xué)性質(zhì),如熱穩(wěn)定性好、揮發(fā)性較低和工作溫度寬等優(yōu)點,但離子液體高黏度的特點是雙電層超級電容器功率性能較差的主要原因[16-17]。通??梢栽陔x子液體中加入有機溶劑改善離子傳輸,進而提高超級電容器的功率密度[18]。乙腈ACN(Acetonitrile)作為最簡單的腈類分子,其具有高介電常數(shù)而且可以和水以任意比例互溶,是常用電容器的溶劑,同時也是液相色譜中的流動相溶劑[19-20]。Chaban 等[21]在先前的模擬工作中報道了乙腈的加入能夠顯著提高咪唑類離子液體的電導(dǎo)率,其原因在于乙腈溶劑穩(wěn)定了離子對,而更大尺寸的離子簇被破壞。Frolov 等[22]針對碳納米管與幾種不同室溫離子液體在乙腈混合物中的相互作用做了大量分子動力學(xué)模擬工作,研究中發(fā)現(xiàn)離子液體的所有成分都在碳納米管表面富集,甚至在未充電的碳納米管表面也形成了幾個不同的層。Wang 等[23-24]在乙腈分子和二氧化硅表面所構(gòu)成的界面體系的模擬中,通過修整二氧化硅表面極性,發(fā)現(xiàn)乙腈分子的結(jié)構(gòu)和動力學(xué)會隨著表面的極性的改變而發(fā)生顯著變化。

    綜上,基于乙腈和離子液體混合溶劑作為電解質(zhì),使用碳納米管等作為電極的雙電層電容器的模擬研究已經(jīng)有一部分展開。但是到目前為止,尚未有和乙腈在雙電層電容器碳極板表面的結(jié)構(gòu)動力學(xué)相關(guān)的詳盡研究。在服役工況下的雙電層電容器中,極板通常帶有一定量的電荷,且通常極板上會有微小的雜質(zhì)[22,25]或多孔結(jié)構(gòu)形成的團簇[6,26],這些因素對極板表面乙腈結(jié)構(gòu)和動力學(xué)的影響值得深入探究。從乙腈在平板石墨烯電極表面的體系出發(fā),本文采用分子動力學(xué)模擬手段,依次對極板帶電和極板修飾兩種條件下乙腈在極板表面的結(jié)構(gòu)和動力學(xué)性質(zhì)進行了模擬。

    1 分子動力學(xué)模擬

    模擬體系由兩塊石墨烯極板和乙腈分子構(gòu)成,如圖1(a)所示。每塊石墨烯極板由4 層碳原子構(gòu)成,石墨烯極板的尺寸為5.15 nm × 4.25 nm。在最初的體系設(shè)定中,兩塊石墨烯極板均為電中性,它們之間的距離為12 nm??偣灿? 787 個乙腈分子限制在兩組石墨烯極板之間,這樣體系中乙腈的平均密度等于343 K 和一個大氣壓下乙腈的體相密度。模擬中,在x、y方向施加周期性邊界條件;而由于界面的存在,非周期性邊界條件施加在z方向。所有的分子動力學(xué)模擬都是在溫度為343 K 的正則系綜(NVT) 下使用Nose-Hoover 熱?。?7-28]進行。

    圖1 平板石墨烯電極-乙腈和團簇修飾石墨烯電極-乙腈體系示意圖Fig.1 Schematic diagram of planar graphene-acetonitrile system and cluster modified graphene-acetonitrile system (white,cyan and blue spheres represent hydrogen,carbon and nitrogen atoms)

    在模擬的力場參數(shù)選擇方面,乙腈選用的是Nikitin 和Lyubartsev 開發(fā)的全原子模型[29],該模型能夠很好的給出乙腈的密度、擴散系數(shù)等物性。石墨烯極板上的碳原子力場采用Amber96 力場中sp2雜化碳原子的參數(shù)[30]。計算模擬過程中,體系中非鍵相互作用的截斷半徑是1.4 nm。長程的靜電相互作用采用PPPM(Particle-particle Particlemesh)[31]算法處理。同時考慮到系統(tǒng)在z方向是非周期性的,在PPPM 計算中我們使用了Berkowitz等[32-33]提出的修正項。全部的分子動力學(xué)模擬均使用LAMMPS[34]軟件運行。

    在電中性平板石墨烯電極的基礎(chǔ)上,我們對極板做了兩種不同設(shè)置以模擬不同條件下乙腈在極板附近的結(jié)構(gòu)以及動力學(xué)性質(zhì)。第一種設(shè)置是在每塊極板上修飾兩個寬度為1 nm 的納米團簇,團簇由碳原子構(gòu)成,兩個團簇之間的距離為1.22 nm,如圖1(b)所示。我們將團簇對應(yīng)的碳原子層數(shù)從一層逐漸增加到五層,這對應(yīng)于團簇的高度從0.334 nm 到1.67 nm。第二種設(shè)置是讓石墨烯極板帶一定量電荷,這樣可以模擬在不同充放電狀態(tài)下乙腈的結(jié)構(gòu)和動力學(xué)性質(zhì)。在帶電體系中,左邊的石墨烯極板帶正電,右邊的石墨烯極板帶等量電荷的負電,這樣體系整體保持電中性。極板的電荷密度從0 開始,以2.4μC/cm2為間隔逐漸增加到12.0μC/cm2,總共模擬了5個帶電體系。

    對于所有的模擬體系,我們首先預(yù)平衡5 ns(時間步長為1 fs),之后再運行10 ns,并收集軌跡做數(shù)據(jù)分析,其時間步長為2 fs。

    2 分析與討論

    2.1 平板與團簇修飾體系

    2.1.1 密度分析我們計算了平板體系和團簇修飾體系的乙腈氮原子的密度(數(shù)密度)分布。不同體系氮原子的密度分布,如圖2 所示。其中,z= 0 nm對應(yīng)于石墨烯極板最靠近乙腈的頂層碳原子的位置。對于平板體系和修飾體系,密度分布的第一個峰出現(xiàn)在z= 0.295 nm 左右(圖2)。對于平板體系的最靠近極板的第一個峰,其峰值約為體相密度的5倍,這個密度值比較接近之前親水二氧化硅表面乙腈的密度[23,25]。這說明盡管本研究中使用的石墨烯電極表面并未被任何極性基團修飾,但由于極板與乙腈之間存在范德華相互作用的吸引力,因此乙腈在平板電極表面仍然能夠有比較強的吸附。隨著乙腈分子逐漸遠離極板,極板對于乙腈的吸引逐漸減弱,當(dāng)z>3 nm 后,乙腈密度值逐漸趨向于體相密度。當(dāng)石墨烯極板表面被團簇修飾后,第一個峰的位置雖然沒有發(fā)生變化,但是峰值較平板體系已經(jīng)減小了近2/3。很明顯,這種現(xiàn)象出現(xiàn)的原因是修飾團簇對乙腈存在空間位阻效應(yīng),并且從圖2中第一個峰也可以看出,無論修飾團簇高度如何,其修飾體系的第一個峰的數(shù)密度都幾乎不變,這表明修飾團簇對石墨烯表面附近的空間位阻效應(yīng)不受其團簇高度的影響。隨著乙腈向體相延伸,團簇高度與乙腈的數(shù)密度呈現(xiàn)出強烈的同步效應(yīng)。當(dāng)修飾高度為一層時,這個體系在z= 0.685 nm 左右(第二個峰)出現(xiàn)最高峰,當(dāng)修飾高度為二層時,這個體系在z= 1.015 nm 左右(第三個峰)出現(xiàn)最高峰,依次類推。簡而言之,修飾團簇層數(shù)為n時,它在n+1 個峰會出現(xiàn)最大值(與其他修飾體系和平板體系對比)。

    圖2 平板電極(黑色曲線)與不同層數(shù)團簇修飾電極體系(彩色曲線)的乙腈氮原子數(shù)密度Fig.2 ACN nitrogen atom number densities of planar electrode(black curve)and modified electrode with different number of layers of clusters(colored curve)

    2.1.2 取向分析為了表征乙腈分子在界面上的分子取向,我們計算了乙腈分子從氮原子指向氰基碳原子的N—C 向量在距離石墨烯極板表面不同距離內(nèi)的區(qū)域的取向分布。因為乙腈分子的N—C向量與其偶極向量相平行,所以可用其來表征乙腈偶極在界面上的排布。圖3給出了不同修飾層數(shù)下乙腈N—C 向量的取向在不同區(qū)域的分布,其中θ是N—C 向量和z軸,也就是垂直于碳極板的法向的夾角。為了更好地分析界面不同區(qū)域的取向分布,我們將原來界面區(qū)域分為三個層,其z軸坐標(biāo)范圍依次為:第一層,z= 0.296~1.015 nm;第二層,z=1.015~1.775nm;第三層,z= 1.775~2.62 nm。如圖3(a)所示,在最靠近石墨烯極板的第一層,幾乎所有的取向分布都在三個角度值(cosθ=-1.0,0.1,1.0)附近有明顯的分布,這三個角度分別對應(yīng)于N—C 向量垂直于表面指向碳極板、平行于極板表面(略微傾斜一定角度)、垂直于表面指向乙腈體相。cosθ= ±1.0 處的峰值說明存在排列取向相反的乙腈分子,乙腈可能在碳極板表面形成了類似磷脂雙分子層的雙層結(jié)構(gòu)[23,25]。

    對于平板體系,cosθ= 1.0 處的峰值并不算太高,大多數(shù)第一層的乙腈分子的N—C 向量仍然傾向于平行于表面。隨著修飾團簇的層數(shù)的增加,cosθ= ±1.0 處的峰值出現(xiàn)了明顯的增加,說明由于團簇的空間位阻效應(yīng),有更多的乙腈分子傾向于垂直于極板排列。

    如圖3(b)和3(c)所示,隨著乙腈分子逐漸遠離極板,在第二層和第三層的乙腈分子的取向分布相對變得更均勻。對于一層修飾和二層修飾的體系,乙腈分子在第二層和第三層的取向分布逐漸趨近于一條水平線,也就是各個角度的分布都存在,這非常類似于體相中乙腈分子的取向分布。從圖3(c)可以看出,對于五層原子團簇修飾的體系,即便在距離表面較遠的第三層,乙腈分子的取向分布仍然顯示出強烈的不均勻性,說明修飾團簇擾動了其周圍的乙腈分子的排布。

    圖3 平板(黑色)與團簇修飾體系(彩色)的乙腈N—C向量取向分布Fig.3 ACN N—C vector orientation distributions for planar(black)electrode system and cluster modified electrode system(colored)

    2.1.3 擴散系數(shù)在體相液體中,分子的擴散一般使用一個擴散系數(shù)來表征,即使用經(jīng)典的Einstein 公式計算均方位移MSD(Mean Squared Displacement)對時間的斜率或者通過Green-Kubo 公式對速度自關(guān)聯(lián)函數(shù)進行積分。由于在界面上分子的擴散是各向異性的,因此無法如體相那樣用一個單一的擴散系數(shù)進行表征,而需要考慮擴散張量的對角元[35-37]。在界面體系中擴散張量僅包含三個對角元Dxx,Dyy,Dzz,其中Dxx=Dyy通常表示為D||,用于表征平行于表面的擴散,而Dzz用于表征垂直于表面的擴散。在本項研究中,我們主要計算了在固液界面中乙腈分子平行于表面的二維擴散系數(shù),可通過計算乙腈分子的二維均方位移得到。二維均方位移的計算公式為

    在獲得固液界面中乙腈分子的二維均方位移后,我們通過式(2)計算固液界面中乙腈分子的二維擴散系數(shù)D||[37-38],其公式為

    圖4給出了不同修飾表面體系中平行于表面的擴散系數(shù)隨碳極板距離變化的趨勢。為了更好地得到極板附近的擴散性質(zhì),我們分別計算了六個亞層中的擴散系數(shù)。亞層的定義如下:第一亞層z= 1.28~1.56 nm;第二亞層z= 1.56~1.95 nm,第三亞層z= 1.95~2.35 nm,第四亞層z= 2.35~2.72 nm,第五亞層z=2.72~3.21nm,第六亞層z=3.21~3.60 nm。這些亞層的定義大體上遵循圖2 中各密度峰的范圍。如圖所示,擴散系數(shù)隨著修飾層數(shù)以及距離碳極板表面的距離都呈單調(diào)變化的趨勢。首先,對于所有的表面層,隨著修飾團簇的尺寸逐步增加,擴散系數(shù)呈逐漸下降的趨勢,這說明團簇的修飾限制了表面上乙腈分子的活動空間;其次,隨著乙腈分子逐漸遠離碳極板,擴散系數(shù)逐漸增加,說明修飾團簇對于乙腈分子的擴散有阻礙作用。即便對于無修飾的平板表面,其最靠近極板的乙腈分子受到極板的吸引力,其擴散系數(shù)也低于遠離表面。

    圖4 平板體系(黑色)與修飾體系(彩色)擴散系數(shù)Fig.4 Diffusion coefficient of planar electrode system(black)and cluster modified electrode system(colored)

    2.2 平板帶電體系

    2.2.1 密度分析我們分別對左右極板計算了帶電石墨烯極板體系的腈基氮原子密度分布,如圖5所示。乙腈的第一個峰出現(xiàn)在z= 0.28 nm,這說明由于極板的范德華力和靜電相互作用的吸引,使得乙腈能夠在極板上有一定的吸附。值得注意的是,第一個峰在z=0.47 nm 處出現(xiàn)了一個肩峰,這與之前Hu 等[18]研究乙腈在二氧化硅表面的密度的特征類似,說明在石墨烯表面很可能存在兩層取向相反的乙腈分子。進一步地,我們觀察到第二個峰出現(xiàn)在z=0.73 nm處,隨著極板表面正電荷密度的增加,峰值呈現(xiàn)不斷擴大的趨勢,如圖5(a)所示。與此對應(yīng)的,在帶負電荷極板,我們發(fā)現(xiàn)了一些更讓人出乎意料的結(jié)果。如圖5(b)所示,對于靠近右側(cè)極板的第一個峰,均出現(xiàn)了雙重峰的現(xiàn)象,且當(dāng)電荷密度0、2.4、4.8μC/cm2時,肩峰均在主峰的左側(cè);而隨著極板電荷量的增加,主峰和肩峰的位置發(fā)生了對換。這說明隨著極板負電荷的增加,對帶負電的乙腈氮原子的排斥增強,使更多乙腈氮原子遠離極板表面。一方面,肩峰的出現(xiàn)如前面所述,石墨烯的表面很可能存在兩層取向相反的乙腈分子;另一方面肩峰的出現(xiàn)和極板表面的電荷量存在極強的耦合性,且當(dāng)電荷密度為-9.6 μC/cm2和-12.0 μC/cm2時,發(fā)現(xiàn)其左邊的肩峰幾乎可以忽略不計,我們猜想由于電荷密度的增大所帶來的排斥作用,導(dǎo)致乙腈分子無法在z= 0.28 nm 左右處駐留。密度曲線繼續(xù)向體相延伸,在第二重峰時肩峰特征已經(jīng)消失,此時隨著電荷的增加,峰的位置向體相中間偏移。整體上來看,峰偏移的同時也有峰值增大的趨勢。

    圖5 不同電荷密度下氮原子數(shù)密度圖Fig.5 ACN nitrogen atom number densities under different charge densities

    通過對正負電極板不同電荷密度下表面乙腈密度的分析,發(fā)現(xiàn)密度峰的位置和峰高隨電荷密度都有一定程度的改變,這說明乙腈的取向隨電場強度增加發(fā)生了變化。在接下來的分析中,我們將通過取向分析揭示這一變化。

    2.2.2 取向分析為了更進一步地揭示帶電體系界面處的結(jié)構(gòu),我們計算了帶電體系在極板界面第一層乙腈分子的N—C 向量與垂直于極板表面的z軸夾角的分布,如圖6 所示。當(dāng)極板所帶正電荷密度逐漸增加,在cosθ= 1.0 處的概率密度越來越大,這說明乙腈分子傾向于垂直于極板表面排列,同時乙腈的氮原子更加靠近表面。值得我們注意的是,當(dāng)電荷密度為9.6 μC/cm2和12.0 μC/cm2時,取向分布在cosθ=1.0 近乎是單調(diào)上升的趨勢。這正好和我們在2.2.1節(jié)的推測呼應(yīng):電荷密度的增大不僅讓乙腈分子取向變化并且也使分子遠離極板駐留。其中最直接的原因是乙腈的分子取向發(fā)生改變,例如在正極板附近,遠離腈基的帶正電的一端被極板強烈排斥,導(dǎo)致分子發(fā)生旋轉(zhuǎn),進而取向發(fā)生改變。相應(yīng)地,當(dāng)極板所帶的負電荷密度逐漸增加時,在cosθ= -1.0 處的概率密度越來越大,這說明乙腈分子傾向于垂直于極板表面排列,同時乙腈的帶正電的甲基更加靠近表面,能夠被帶正電的表面吸引。

    圖6 不同電荷密度下第一層乙腈分子的取向分布Fig.6 Orientation distribution of the first layer of ACN under different charge densities

    2.2.3 擴散系數(shù)圖7 給出了不同帶電極板體系中平行于表面的擴散系數(shù)在距離碳極板表面不同距離的亞層中的分布圖,類似2.1.3節(jié)擴散系數(shù)的分析,我們在這里計算了最靠近極板表面的前三個亞層的擴散系數(shù)。對于所有的體系,隨著乙腈逐漸遠離表面,擴散系數(shù)都出現(xiàn)了一個逐漸增加的趨勢。這說明在帶電極板上,由于受到極板靜電力的吸引,乙腈分子的平移受到了一定程度的限制。另外一個很明顯的趨勢就是,隨著極板電荷密度的增強,乙腈的擴散系數(shù)也出現(xiàn)了下降,這同樣可以歸因為極板吸引力的增強。比較有趣的是,對于帶正電荷的極板,當(dāng)極板從不帶電荷增加到電荷密度為2.4μC/cm2時,擴散系數(shù)反而出現(xiàn)了一個增加。這與帶負電的極板非常不同,乙腈在所有帶負電極板表面的擴散系數(shù)均低于不帶電荷極板的擴散系數(shù)。這種現(xiàn)象的可能原因是在該電荷密度的情況下,乙腈的分子取向從原來的腈基氮靠近極板轉(zhuǎn)變成遠離極板,此時靜電作用強烈排斥乙腈分子,間接導(dǎo)致其擴散速度變快,而對于其他帶電體系,乙腈分子取向雖然也發(fā)生改變,但是由于距離極板較遠,排斥力減小,故擴散速度也相應(yīng)變慢。

    圖7 不同電荷密度下乙腈的擴散系數(shù)Fig.7 Diffusion coefficient of ACN under different charge densities

    3 結(jié) 論

    在本文中,我們采用分子動力學(xué)模擬研究了超級電容器中團簇修飾以及極板電荷密度對極板表面乙腈溶劑的結(jié)構(gòu)和動力學(xué)的影響。對于團簇修飾,模擬發(fā)現(xiàn):團簇會對乙腈的吸附起一定的阻礙作用。且隨著修飾團簇高度的增加,這種阻礙作用就越強。這種阻礙效應(yīng)會降低乙腈分子在極板表面的吸附,并使分子取向從平行于極板表面轉(zhuǎn)變?yōu)榇怪庇跇O板表面。同時,團簇的存在會使界面上分子的擴散動力學(xué)減慢。當(dāng)極板處于帶電狀態(tài)下時,極板電荷密度的增強不僅增強了乙腈分子的吸附,同時使得乙腈分子的取向發(fā)生了明顯改變,進而影響了乙腈在極板表面的擴散。模擬工作揭示了乙腈在復(fù)雜極板表面上的結(jié)構(gòu)和動力學(xué)性質(zhì)與平板表面非常不同,這將有助于未來深入理解乙腈離子液體混合溶液在復(fù)雜極板表面上的電化學(xué)性能。

    猜你喜歡
    極板擴散系數(shù)電荷
    連續(xù)分布電荷體系電荷元的自能問題*
    物理通報(2024年4期)2024-04-09 12:41:28
    電荷知識知多少
    一種新型微球聚焦測量極板設(shè)計及應(yīng)用*
    電荷守恒在化學(xué)解題中的應(yīng)用
    電容器與帶電粒子在電場中的運動測試題
    燃料電池雙極板制備工藝及其熱壓機設(shè)計研究
    基于Sauer-Freise 方法的Co- Mn 體系fcc 相互擴散系數(shù)的研究
    上海金屬(2015年5期)2015-11-29 01:13:59
    FCC Ni-Cu 及Ni-Mn 合金互擴散系數(shù)測定
    上海金屬(2015年6期)2015-11-29 01:09:09
    靜電現(xiàn)象有什么用?
    關(guān)于鉛酸蓄電池中正極板柵溶解問題的研究
    国产黄色免费在线视频| 欧美3d第一页| 日韩在线高清观看一区二区三区| 制服丝袜香蕉在线| 亚洲国产av新网站| 午夜福利视频精品| 亚洲精品中文字幕在线视频| 亚洲精品一二三| 亚洲欧美精品自产自拍| 成人国语在线视频| 一本久久精品| 亚洲av在线观看美女高潮| 国产一区二区三区综合在线观看 | 日韩大片免费观看网站| 一级毛片我不卡| 纯流量卡能插随身wifi吗| 欧美97在线视频| 九九在线视频观看精品| 日韩三级伦理在线观看| 欧美日韩国产mv在线观看视频| 特大巨黑吊av在线直播| 亚洲性久久影院| 亚洲精品久久久久久婷婷小说| 亚洲综合色惰| 免费不卡的大黄色大毛片视频在线观看| 国产欧美另类精品又又久久亚洲欧美| 男女边摸边吃奶| 人妻夜夜爽99麻豆av| 日本午夜av视频| 久久久久久久久大av| 精品久久久久久久久亚洲| 亚洲欧洲国产日韩| 免费不卡的大黄色大毛片视频在线观看| 亚洲欧洲精品一区二区精品久久久 | 黑丝袜美女国产一区| 91午夜精品亚洲一区二区三区| 亚洲色图 男人天堂 中文字幕 | 街头女战士在线观看网站| 桃花免费在线播放| 亚洲国产最新在线播放| 精品少妇内射三级| 亚洲经典国产精华液单| 日韩免费高清中文字幕av| 久久韩国三级中文字幕| 久久午夜综合久久蜜桃| 国产有黄有色有爽视频| 久久免费观看电影| 日韩中字成人| 大片电影免费在线观看免费| av国产精品久久久久影院| 一级毛片我不卡| 久久人人爽人人片av| 亚洲成人一二三区av| 在线观看人妻少妇| 亚洲国产精品999| 最近中文字幕2019免费版| 成人毛片a级毛片在线播放| 久久免费观看电影| 黑人猛操日本美女一级片| 曰老女人黄片| 熟女av电影| 汤姆久久久久久久影院中文字幕| 日本爱情动作片www.在线观看| 亚洲精品第二区| 国产成人午夜福利电影在线观看| 如日韩欧美国产精品一区二区三区 | 久久影院123| 日韩一区二区三区影片| 国产成人精品福利久久| av免费在线看不卡| 下体分泌物呈黄色| 精品视频人人做人人爽| 国产一区二区三区av在线| 久久99热这里只频精品6学生| 国产无遮挡羞羞视频在线观看| 伦精品一区二区三区| 中文欧美无线码| 日韩成人av中文字幕在线观看| 妹子高潮喷水视频| .国产精品久久| 亚洲国产av影院在线观看| 亚洲国产最新在线播放| 久久婷婷青草| 亚洲美女黄色视频免费看| 亚洲精品久久午夜乱码| 亚洲美女搞黄在线观看| 成年美女黄网站色视频大全免费 | 精品少妇黑人巨大在线播放| 亚洲国产精品专区欧美| 黄色视频在线播放观看不卡| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久这里有精品视频免费| 日韩人妻高清精品专区| 久久久久久久精品精品| 午夜精品国产一区二区电影| 久久鲁丝午夜福利片| 国产一区二区在线观看日韩| 精品少妇内射三级| 日本欧美视频一区| 22中文网久久字幕| 国精品久久久久久国模美| 国产深夜福利视频在线观看| 伊人亚洲综合成人网| 中国国产av一级| 91精品国产九色| 国产欧美亚洲国产| 日本vs欧美在线观看视频| 最近手机中文字幕大全| 亚洲欧美色中文字幕在线| 免费观看在线日韩| 国产视频首页在线观看| 18禁观看日本| 亚洲精品aⅴ在线观看| 一级毛片电影观看| 亚洲成人一二三区av| 国产精品一区二区三区四区免费观看| 亚洲精品色激情综合| 中文乱码字字幕精品一区二区三区| 亚洲精品久久成人aⅴ小说 | 亚洲国产av新网站| 欧美日韩精品成人综合77777| 69精品国产乱码久久久| 国国产精品蜜臀av免费| 精品久久久精品久久久| 国产精品人妻久久久久久| 日韩在线高清观看一区二区三区| 大话2 男鬼变身卡| 又大又黄又爽视频免费| 久久精品夜色国产| 极品少妇高潮喷水抽搐| 亚洲精品久久成人aⅴ小说 | 女性生殖器流出的白浆| 国产成人精品久久久久久| av在线观看视频网站免费| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品久久久久久av不卡| 性色av一级| av专区在线播放| 亚洲国产精品一区二区三区在线| 精品国产一区二区久久| 视频中文字幕在线观看| 97精品久久久久久久久久精品| 国产亚洲av片在线观看秒播厂| 国产伦理片在线播放av一区| 91久久精品电影网| 久久青草综合色| 18禁动态无遮挡网站| 高清午夜精品一区二区三区| 九色成人免费人妻av| 99热6这里只有精品| 国产在线免费精品| 69精品国产乱码久久久| 久久久久久久久大av| 飞空精品影院首页| 婷婷色av中文字幕| 极品少妇高潮喷水抽搐| 人妻制服诱惑在线中文字幕| 国产精品偷伦视频观看了| 久久久久久伊人网av| 18禁在线播放成人免费| 久久国产亚洲av麻豆专区| 午夜激情福利司机影院| 欧美97在线视频| 99国产综合亚洲精品| 日韩中字成人| 校园人妻丝袜中文字幕| 国产69精品久久久久777片| 赤兔流量卡办理| 亚洲欧美精品自产自拍| 99久久中文字幕三级久久日本| 91久久精品国产一区二区成人| 蜜桃国产av成人99| 中国美白少妇内射xxxbb| 三级国产精品欧美在线观看| 日韩三级伦理在线观看| 久久久久网色| 亚洲三级黄色毛片| 免费人妻精品一区二区三区视频| 国产免费一区二区三区四区乱码| 欧美老熟妇乱子伦牲交| 欧美精品国产亚洲| 18禁观看日本| 免费人妻精品一区二区三区视频| 亚洲精品日本国产第一区| 亚洲无线观看免费| 国产精品人妻久久久影院| 一级黄片播放器| 黄色一级大片看看| a级毛片免费高清观看在线播放| 中文字幕免费在线视频6| 免费观看av网站的网址| 久久97久久精品| 精品一区二区三卡| 69精品国产乱码久久久| 亚洲五月色婷婷综合| 伦理电影大哥的女人| 国产男女内射视频| 中文字幕亚洲精品专区| 两个人的视频大全免费| 精品国产国语对白av| 久久婷婷青草| 视频区图区小说| 欧美精品亚洲一区二区| 亚洲成人av在线免费| 亚洲欧洲精品一区二区精品久久久 | 少妇精品久久久久久久| 日韩伦理黄色片| 亚洲欧美中文字幕日韩二区| av有码第一页| 青春草亚洲视频在线观看| 十八禁高潮呻吟视频| 亚洲激情五月婷婷啪啪| 免费日韩欧美在线观看| 国产白丝娇喘喷水9色精品| 97在线视频观看| 久久久久久久国产电影| 蜜臀久久99精品久久宅男| 亚洲国产最新在线播放| 国产一区二区三区综合在线观看 | 国产精品一区www在线观看| 大香蕉久久网| 中国国产av一级| 国产亚洲av片在线观看秒播厂| 免费人妻精品一区二区三区视频| a 毛片基地| www.av在线官网国产| 爱豆传媒免费全集在线观看| 午夜精品国产一区二区电影| 国产午夜精品久久久久久一区二区三区| 99久久精品国产国产毛片| 欧美人与善性xxx| 亚洲av福利一区| 丝袜在线中文字幕| 亚洲av中文av极速乱| 国产免费一区二区三区四区乱码| 十八禁网站网址无遮挡| 精品一区在线观看国产| 国产亚洲精品第一综合不卡 | 免费观看性生交大片5| 亚洲精品久久成人aⅴ小说 | 国产在线一区二区三区精| 汤姆久久久久久久影院中文字幕| 简卡轻食公司| 熟女电影av网| 嫩草影院入口| 国产在线视频一区二区| 18禁在线播放成人免费| 久久这里有精品视频免费| 国产伦精品一区二区三区视频9| 亚洲欧美日韩另类电影网站| 亚洲精品aⅴ在线观看| 亚洲成人手机| 黄色一级大片看看| 欧美日韩视频高清一区二区三区二| 免费少妇av软件| 视频区图区小说| a级毛色黄片| 欧美三级亚洲精品| 国精品久久久久久国模美| 考比视频在线观看| 有码 亚洲区| 老熟女久久久| 又大又黄又爽视频免费| 亚洲国产欧美日韩在线播放| 91国产中文字幕| 亚洲国产色片| 亚洲第一区二区三区不卡| 免费观看无遮挡的男女| 在线精品无人区一区二区三| 九九在线视频观看精品| 精品卡一卡二卡四卡免费| a级片在线免费高清观看视频| 97超视频在线观看视频| 久久精品久久久久久久性| 91久久精品电影网| 国产老妇伦熟女老妇高清| 国产成人精品在线电影| 最近中文字幕高清免费大全6| 成年女人在线观看亚洲视频| 欧美精品国产亚洲| 国产日韩欧美在线精品| 爱豆传媒免费全集在线观看| 老女人水多毛片| 精品国产乱码久久久久久小说| 少妇熟女欧美另类| 日韩人妻高清精品专区| 黄色视频在线播放观看不卡| 91久久精品电影网| 夜夜骑夜夜射夜夜干| 欧美日本中文国产一区发布| 日韩成人av中文字幕在线观看| 考比视频在线观看| 欧美另类一区| 一级毛片aaaaaa免费看小| 亚洲av.av天堂| 人妻少妇偷人精品九色| 国产免费福利视频在线观看| 2021少妇久久久久久久久久久| 国产精品.久久久| 青青草视频在线视频观看| 精品一品国产午夜福利视频| 精品午夜福利在线看| 乱人伦中国视频| 91aial.com中文字幕在线观看| 最近2019中文字幕mv第一页| 精品人妻偷拍中文字幕| 九色亚洲精品在线播放| 日本-黄色视频高清免费观看| 99久久综合免费| 国产精品99久久久久久久久| 黄片播放在线免费| 国产成人午夜福利电影在线观看| 伦理电影免费视频| kizo精华| 十分钟在线观看高清视频www| 免费黄色在线免费观看| 赤兔流量卡办理| 高清黄色对白视频在线免费看| 777米奇影视久久| 99热全是精品| 亚洲成色77777| 亚洲一级一片aⅴ在线观看| 日韩成人av中文字幕在线观看| 亚洲欧美日韩另类电影网站| 日韩在线高清观看一区二区三区| 精品熟女少妇av免费看| √禁漫天堂资源中文www| 婷婷成人精品国产| 99精国产麻豆久久婷婷| 多毛熟女@视频| 中文字幕亚洲精品专区| 制服丝袜香蕉在线| 国产精品女同一区二区软件| 精品亚洲乱码少妇综合久久| 精品人妻偷拍中文字幕| 啦啦啦中文免费视频观看日本| 国产色爽女视频免费观看| 日产精品乱码卡一卡2卡三| 免费日韩欧美在线观看| 两个人的视频大全免费| 国产精品99久久99久久久不卡 | 97精品久久久久久久久久精品| 成人毛片a级毛片在线播放| 热99国产精品久久久久久7| 大陆偷拍与自拍| 嫩草影院入口| 插逼视频在线观看| 春色校园在线视频观看| 三级国产精品欧美在线观看| a 毛片基地| 中文天堂在线官网| 狂野欧美激情性bbbbbb| 日韩一本色道免费dvd| 精品一区二区三区视频在线| 久久免费观看电影| 国产乱来视频区| 亚洲精品日韩在线中文字幕| 国产日韩欧美视频二区| 91aial.com中文字幕在线观看| 丰满饥渴人妻一区二区三| av在线观看视频网站免费| 国产精品国产三级国产专区5o| freevideosex欧美| 亚洲成人一二三区av| 欧美日韩在线观看h| 天堂俺去俺来也www色官网| 亚洲av中文av极速乱| 亚洲欧美日韩另类电影网站| 日本与韩国留学比较| 热re99久久精品国产66热6| 91aial.com中文字幕在线观看| 国产69精品久久久久777片| 中文欧美无线码| 婷婷色av中文字幕| 日韩精品有码人妻一区| av在线app专区| 中国三级夫妇交换| 久久国产精品男人的天堂亚洲 | 精品一区二区三区视频在线| 久久久亚洲精品成人影院| 午夜激情av网站| 日韩av在线免费看完整版不卡| 午夜激情福利司机影院| 国产一区二区三区av在线| 美女视频免费永久观看网站| 亚洲av成人精品一二三区| 亚洲精品一区蜜桃| 一本—道久久a久久精品蜜桃钙片| 中文字幕人妻熟人妻熟丝袜美| 日日撸夜夜添| 大又大粗又爽又黄少妇毛片口| 2018国产大陆天天弄谢| 热99久久久久精品小说推荐| 9色porny在线观看| 国产精品一区二区三区四区免费观看| 久久午夜福利片| 一级毛片黄色毛片免费观看视频| 亚洲欧美中文字幕日韩二区| 91午夜精品亚洲一区二区三区| 一边摸一边做爽爽视频免费| 亚洲色图 男人天堂 中文字幕 | 男人添女人高潮全过程视频| 国产片内射在线| 久久av网站| 男女高潮啪啪啪动态图| 国产成人精品一,二区| 亚洲精华国产精华液的使用体验| 人妻制服诱惑在线中文字幕| 精品午夜福利在线看| 亚洲精品日本国产第一区| 国产伦理片在线播放av一区| 国内精品宾馆在线| 哪个播放器可以免费观看大片| 视频中文字幕在线观看| 美女中出高潮动态图| 欧美日韩精品成人综合77777| 少妇人妻精品综合一区二区| 久久人妻熟女aⅴ| 久久久久人妻精品一区果冻| 老女人水多毛片| 亚洲av成人精品一二三区| 狂野欧美白嫩少妇大欣赏| 大香蕉久久网| 黑人猛操日本美女一级片| 国产一级毛片在线| 黄片无遮挡物在线观看| 亚洲成人手机| 99国产精品免费福利视频| 欧美国产精品一级二级三级| 亚洲精品乱码久久久v下载方式| 日韩人妻高清精品专区| 国产精品女同一区二区软件| 亚洲无线观看免费| 亚洲欧洲日产国产| 国产日韩欧美在线精品| 国产成人免费观看mmmm| 日日啪夜夜爽| 日韩精品有码人妻一区| 婷婷成人精品国产| 这个男人来自地球电影免费观看 | 亚洲精品久久久久久婷婷小说| 亚洲欧洲日产国产| 新久久久久国产一级毛片| 亚洲美女视频黄频| 高清毛片免费看| 久久久久精品久久久久真实原创| 亚洲国产欧美日韩在线播放| 亚洲av欧美aⅴ国产| 搡老乐熟女国产| 在线观看免费视频网站a站| 十分钟在线观看高清视频www| 久久久久久久国产电影| 男人操女人黄网站| 免费黄色在线免费观看| 91午夜精品亚洲一区二区三区| 国产老妇伦熟女老妇高清| 肉色欧美久久久久久久蜜桃| 午夜老司机福利剧场| 五月开心婷婷网| 日本av手机在线免费观看| 秋霞伦理黄片| 桃花免费在线播放| 99久久人妻综合| 大香蕉久久成人网| 51国产日韩欧美| 成年人免费黄色播放视频| 国产毛片在线视频| www.色视频.com| 日韩 亚洲 欧美在线| 国产免费福利视频在线观看| 日本欧美国产在线视频| 久久国产精品大桥未久av| 精品久久久精品久久久| 日本免费在线观看一区| xxxhd国产人妻xxx| 91精品三级在线观看| 国产亚洲欧美精品永久| 欧美成人午夜免费资源| 国产国语露脸激情在线看| 三上悠亚av全集在线观看| 国产一区二区在线观看日韩| 亚洲美女搞黄在线观看| 免费播放大片免费观看视频在线观看| 狠狠精品人妻久久久久久综合| 99热国产这里只有精品6| 国产精品久久久久久久久免| 最黄视频免费看| 久久久久久人妻| 国产成人91sexporn| 午夜福利视频在线观看免费| 成人国产av品久久久| 搡老乐熟女国产| 免费人成在线观看视频色| xxx大片免费视频| 99久久人妻综合| 美女主播在线视频| 秋霞在线观看毛片| 欧美精品一区二区大全| 一本一本综合久久| 国产精品久久久久久久电影| 少妇猛男粗大的猛烈进出视频| 欧美人与善性xxx| 中国国产av一级| 国产不卡av网站在线观看| 多毛熟女@视频| 国产av一区二区精品久久| 国产男女内射视频| 一级毛片aaaaaa免费看小| 2018国产大陆天天弄谢| 日本欧美视频一区| 少妇丰满av| 亚洲精品成人av观看孕妇| 男女免费视频国产| 99九九在线精品视频| 久久久久久久久久久免费av| 午夜免费鲁丝| 亚洲久久久国产精品| 成人18禁高潮啪啪吃奶动态图 | 男女边摸边吃奶| 亚洲精品成人av观看孕妇| 日韩一本色道免费dvd| 国产亚洲一区二区精品| 日本vs欧美在线观看视频| 日韩电影二区| 啦啦啦中文免费视频观看日本| 中文字幕人妻丝袜制服| 青春草国产在线视频| 一区二区av电影网| 久久毛片免费看一区二区三区| 国产成人a∨麻豆精品| 中文天堂在线官网| 亚洲精品国产av蜜桃| 2018国产大陆天天弄谢| 大香蕉久久网| 狂野欧美激情性bbbbbb| 国产伦精品一区二区三区视频9| 两个人免费观看高清视频| 国产成人a∨麻豆精品| av网站免费在线观看视频| 久久99热6这里只有精品| 99热国产这里只有精品6| 亚洲欧美色中文字幕在线| 国产片内射在线| 亚洲成人av在线免费| 亚洲精品成人av观看孕妇| 精品久久久久久电影网| a级片在线免费高清观看视频| 午夜91福利影院| 大香蕉久久成人网| 国产日韩欧美在线精品| 久久精品久久精品一区二区三区| 日本黄大片高清| 婷婷色av中文字幕| 十分钟在线观看高清视频www| a级片在线免费高清观看视频| 午夜福利影视在线免费观看| 亚洲精品久久久久久婷婷小说| 涩涩av久久男人的天堂| 亚洲情色 制服丝袜| 夜夜看夜夜爽夜夜摸| 国产日韩欧美视频二区| 亚洲伊人久久精品综合| 国产成人一区二区在线| 五月天丁香电影| 久久精品久久久久久噜噜老黄| 午夜视频国产福利| 一区二区三区乱码不卡18| 91aial.com中文字幕在线观看| 夜夜爽夜夜爽视频| 一区二区三区四区激情视频| 久久久亚洲精品成人影院| 97在线人人人人妻| 国产成人av激情在线播放 | av卡一久久| 日韩一本色道免费dvd| av黄色大香蕉| 午夜免费鲁丝| 久久久久久久久久久久大奶| 亚洲精品一区蜜桃| 黄片播放在线免费| 大香蕉久久网| 久久精品国产亚洲网站| 国产精品人妻久久久影院| 岛国毛片在线播放| 久久人人爽人人片av| 十分钟在线观看高清视频www| 91精品三级在线观看| 下体分泌物呈黄色| 亚洲美女视频黄频| 亚洲国产精品999| 女性被躁到高潮视频| 精品国产乱码久久久久久小说| 另类亚洲欧美激情| 美女cb高潮喷水在线观看| 观看av在线不卡| 国产午夜精品一二区理论片| 美女cb高潮喷水在线观看| 日韩伦理黄色片| 久久亚洲国产成人精品v| 五月天丁香电影| 国产一区二区在线观看av| 国产精品人妻久久久久久| 热99国产精品久久久久久7| 视频区图区小说| 午夜福利,免费看| 热99国产精品久久久久久7| 视频区图区小说| 人妻一区二区av| 女性被躁到高潮视频| 日韩精品有码人妻一区| 国产伦理片在线播放av一区| a级毛片黄视频| 国产成人免费观看mmmm| 九草在线视频观看| 青春草国产在线视频|