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

    基于第一性原理計(jì)算CuxNi19-x(x<19)混合團(tuán)簇的結(jié)構(gòu)研究①

    2011-12-26 07:13:38馮翠菊
    關(guān)鍵詞:原子間第一性偏析

    馮翠菊 劉 英

    (1.華北科技學(xué)院基礎(chǔ)部,北京 東燕郊 101601;2.河北師范大學(xué)物理與信息工程學(xué)院,河北 石家莊 050800)

    基于第一性原理計(jì)算CuxNi19-x(x<19)混合團(tuán)簇的結(jié)構(gòu)研究①

    馮翠菊1②劉 英2

    (1.華北科技學(xué)院基礎(chǔ)部,北京 東燕郊 101601;2.河北師范大學(xué)物理與信息工程學(xué)院,河北 石家莊 050800)

    采用基于第一性原理計(jì)算的Chen-M?bius三維晶格反演基礎(chǔ)上獲得的原子間相互作用勢參數(shù),運(yùn)用最陡下降法和共軛梯度法結(jié)合起來的最優(yōu)化計(jì)算Smart方法得到了CuxNi19-x(x<19)二元過渡金屬混合團(tuán)簇的穩(wěn)態(tài)構(gòu)型,發(fā)現(xiàn)Cu-Ni混合團(tuán)簇中有明顯Cu元素偏析現(xiàn)象,即Cu和Ni并沒有互相混合形成有序結(jié)構(gòu),而是分別聚集在一起,所形成的體系對稱性很低,呈割據(jù)狀態(tài)。

    晶格反演,Cu-Ni混合團(tuán)簇,穩(wěn)態(tài)結(jié)構(gòu)

    0 引言

    團(tuán)簇結(jié)構(gòu)和性質(zhì)的研究對于理解物質(zhì)從微觀到宏觀的過渡具有重要作用,團(tuán)簇結(jié)構(gòu)研究的基本問題是弄清團(tuán)簇如何由原子、分子一步步發(fā)展而成,以及隨著這種發(fā)展,團(tuán)簇的結(jié)構(gòu)和性質(zhì)如何變化。混合團(tuán)簇可以極大的擴(kuò)展單質(zhì)團(tuán)簇的性質(zhì),而且通過混合可以得到和加強(qiáng)團(tuán)簇的某些特殊性質(zhì)。由于組成、結(jié)構(gòu)及性質(zhì)的豐富多樣性,混合團(tuán)簇在電子、工程、催化等領(lǐng)域具有廣泛的應(yīng)用前景。

    團(tuán)簇在理論計(jì)算中存在著多種計(jì)算方法。最精確的從頭計(jì)算方法,由于其計(jì)算量與基函數(shù)數(shù)目的五次方或六次方成正比,限制了其在電子數(shù)目較多的體系(即較大尺寸團(tuán)簇中)的應(yīng)用。局域密度泛函方法在小團(tuán)簇計(jì)算中取得了巨大的成功,但當(dāng)團(tuán)簇尺寸增大到20個原子以上時(shí),該方法同樣受到了計(jì)算量的限制,對大中型團(tuán)簇顯得無能為力。目前對于金屬團(tuán)簇的研究主要是基于Gupta 、Sutton-Chen、EAM 等經(jīng)驗(yàn)勢,Gupta勢已成功研究了許多團(tuán)簇的結(jié)構(gòu),現(xiàn)在已有關(guān)于Cu、Ni、Ag、Au 等團(tuán)簇[1]及 Ni-Al、Cu-Co、Cu-Au 團(tuán)簇[2,3]的報(bào)道。

    基于第一性原理計(jì)算的材料的結(jié)合能,運(yùn)用Chen-M?bius反演[4-6],可以嚴(yán)格導(dǎo)出的原子間相互作用對勢,該方法廣泛地應(yīng)用于研究材料的靜力學(xué)和動力學(xué)特性[7-10]。本文將Chen-M?bius晶格反演應(yīng)用推廣到過渡金屬二元團(tuán)簇的研究,運(yùn)用最優(yōu)化方法計(jì)算得到了CuxNi19-x(x<19)二元過渡金屬團(tuán)簇的穩(wěn)態(tài)結(jié)構(gòu),由此開辟了一條基于第一性原理計(jì)算的、較高效率地研究二元金屬團(tuán)簇的結(jié)構(gòu)及其特性的途徑。

    1 勢函數(shù)的獲得與計(jì)算方法

    1.1 勢函數(shù)的獲得

    基于第一性原理計(jì)算的Chen-M?bius晶格反演獲得的原子間相互作用對勢,不依賴于經(jīng)驗(yàn)的形式和方法,直接由金屬和金屬間化合物的第一性原理自洽能帶計(jì)算和嚴(yán)格的晶格反演過程得出,其準(zhǔn)確性僅僅取決于第一性原理計(jì)算的準(zhǔn)確度和對勢近似的可靠性,因此它在材料的物理性質(zhì)的理論計(jì)算中得到了廣泛的應(yīng)用。Chen-M?bius反演的理論要點(diǎn)如下:

    對于同種原子構(gòu)成的晶體。設(shè)從第一性原理算出的晶體結(jié)合能函數(shù)可表示為:

    式中x為最近鄰原子間距,r0(n)為第n階近鄰的配位數(shù),b0(n)是以最近鄰距離為單位的、第n級近鄰到參考原子的距離,Φ(x)為對勢函數(shù)。為了理論上處理的方便,通過{b0(n)}的自乘得到{b(n)},{b(n)}可認(rèn)為構(gòu)成乘法半群,這里多出的虛格點(diǎn)相應(yīng)的虛配位數(shù)均為零。這時(shí)有

    這里

    由此可反演出的原子間相互作用勢普遍公式

    式中I(n)可通過下式求出

    這里求和號下b(d)|b(n)表示對所有這樣的b(d)求和:即b(d)∈{b(n)},且可以在{b(n)}中找到一個{b(k)},使得b(d)b(k)=b(n)。

    可以證明,(2)式是(1)式的解:

    關(guān)于異種原子間相互作用勢問題可以通過以下的方案解決[6],即認(rèn)為通過純元素晶體反演出的原子間相互作用勢在二元合金中仍然有效,如此,可以從合金的總結(jié)合能中扣除同種原子之間相互作用對結(jié)合能的貢獻(xiàn)而得出異種原子的“偏結(jié)合能”曲線,從偏結(jié)合能曲線反演出異種原子之間相互作用勢。

    表一中給出了Chen-M?bius晶格反演得到的原子間相互作用對勢,然后利用Morse勢形式擬合得到的參數(shù),其中Rcut表示計(jì)算原子間相互作用對勢的截止距離。

    1.2 計(jì)算方法

    基于原子間相互作用勢,尋找團(tuán)簇的穩(wěn)態(tài)結(jié)構(gòu)可以通過團(tuán)簇的能量最小化過程實(shí)現(xiàn)。團(tuán)簇的能量最小化過程實(shí)質(zhì)上是一個多元函數(shù)的極值問題,我們采用最陡下降法(Steepest descents)和共軛梯度方法(Conjugate gradient)有機(jī)結(jié)合起來的Smart方法[11]。

    由于在極值點(diǎn)附近能量的梯度趨于零,最陡下降法在系統(tǒng)趨近極值點(diǎn)時(shí)收斂變慢,盡管如此,這種方法由于運(yùn)算穩(wěn)定,仍被廣泛采用或作結(jié)構(gòu)能量最小化的預(yù)處理方法。共軛梯度方法以及各種Newton-Raphson方法對于偏離局域極值較大的系統(tǒng)構(gòu)型(即勢能曲面偏離二次曲面較大)是不穩(wěn)定的。另外,Newton-Raphson方法需要的存儲空間較大,不適于大體系計(jì)算。

    隨著團(tuán)簇的尺度增加,團(tuán)簇可能的構(gòu)型數(shù)目增加非??欤@然我們不能用窮舉的方法得到團(tuán)簇的基態(tài)。實(shí)際計(jì)算中,判斷系統(tǒng)最小化后的構(gòu)型是局域極小還是全域極小,是一個難以回答的問題。在我們的計(jì)算中,采用不同的系統(tǒng)構(gòu)型進(jìn)行多次能量最小化,由各極值點(diǎn)的能量值判斷能量最優(yōu)化結(jié)構(gòu)。對每一個固定原子數(shù)的團(tuán)簇,我們隨機(jī)地選取系統(tǒng)1000~2000個不同的初始構(gòu)型進(jìn)行能量最小化。本文首先采用最陡下降法對團(tuán)簇的結(jié)構(gòu)進(jìn)行預(yù)處理,然后采用共軛梯度方法計(jì)算團(tuán)簇的穩(wěn)態(tài)結(jié)構(gòu),從而找到二元金屬團(tuán)簇的最小能量結(jié)構(gòu)。當(dāng)然所找到的每一固定原子數(shù)的穩(wěn)態(tài)結(jié)構(gòu)不止一種,每種原子數(shù)的團(tuán)簇都分別對應(yīng)著多種穩(wěn)態(tài)結(jié)構(gòu),但其中有些是亞穩(wěn)態(tài),通過對它們的能量比較,我們可以找到(近似的)最小能量結(jié)構(gòu),也就是要找的最穩(wěn)態(tài)結(jié)構(gòu)

    表1 Chen-M?bius晶格反演得到的原子間相互作用對勢參數(shù)

    2 結(jié)果與討論

    圖1列出了CuxNi19-x(x<19)各種成分的穩(wěn)態(tài)構(gòu)型,從中可以看到隨著Cu原子數(shù)的逐漸增大,混合團(tuán)簇CuxNi19-x的穩(wěn)態(tài)構(gòu)型都是三個五邊形金字塔結(jié)構(gòu)互相交叉構(gòu)成的對稱結(jié)構(gòu),即雙二十面體套接結(jié)構(gòu),但由于Cu原子和Ni原子的半徑不同而發(fā)生了輕微的形變。在CuNi18中,唯一的Cu原子占據(jù)了雙二十面體的一個頂點(diǎn),而不是占據(jù)中心位置,這是由于兩種元素具有不同的表面能造成的,Cu的表面能比Ni的表面能低,在Cu原子和Ni原子形成混合團(tuán)簇時(shí),Cu原子會自動占據(jù)團(tuán)簇的表面位置,而Ni卻占據(jù)了內(nèi)部位置,形成了偏析現(xiàn)象。同樣的現(xiàn)象也發(fā)生在x=2、3、4,以及 x=12~18的情況下,即在Cu、Ni原子數(shù)混合比例相差較大時(shí),有明顯的Cu偏析現(xiàn)象。當(dāng)兩者的混合比例相差不大時(shí),我們同樣觀察到偏析現(xiàn)象,Cu和Ni并沒有互相混合形成有序結(jié)構(gòu),而是分別聚集在一起,所形成的體系對稱性很低,呈割據(jù)狀態(tài)。這種割據(jù)狀態(tài)的形成,可以看成Cu、Ni合金微觀不均勻性的具體體現(xiàn),這導(dǎo)致了宏觀上Cu-Ni固溶體的體現(xiàn)。

    圖1 CuXNi19-x的穩(wěn)態(tài)構(gòu)型 深色球代表Ni,淺色代表Cu

    同樣我們也計(jì)算了CuxNi55-x的穩(wěn)態(tài)構(gòu)型,圖2我們給出了Cu8Ni47、Cu27Ni28的穩(wěn)態(tài)構(gòu)型,這些穩(wěn)態(tài)構(gòu)型都是二十面體結(jié)構(gòu),一個原子占據(jù)了二十面體的中心,外殼層包含了兩層原子,最外面一層包含了42個原子,次外層包含了12個原子。由圖可以明顯看出,Cu、Ni原子混合比例相差較大時(shí),Cu原子明顯的偏析,混合比例相近時(shí),Cu也偏析呈明顯的割據(jù)狀態(tài),這一計(jì)算結(jié)果與J.M.Montejano-Carrizales等人[12]計(jì)算結(jié)果相一致,這說明了我們計(jì)算的可靠性。

    圖2 Cu8Ni47Cu27Ni28的穩(wěn)態(tài)構(gòu)型深色球代表Ni,淺色代表Cu

    3 結(jié)論

    本文采用基于第一性原理計(jì)算的Chen-M?bius三維晶格反演基礎(chǔ)上獲得的原子間相互作用勢參數(shù),運(yùn)用最陡下降法和共軛梯度法結(jié)合起來的最優(yōu)化計(jì)算Smart方法得到了Cu-Ni二元過渡金屬混合團(tuán)簇的穩(wěn)態(tài)構(gòu)型,發(fā)現(xiàn)CuNi混合團(tuán)簇的穩(wěn)態(tài)構(gòu)型與純金屬Ni團(tuán)簇的基態(tài)構(gòu)型完全一樣,只是由于Cu、Ni原子半徑的不同而略有畸變;混合團(tuán)簇中Cu原子趨于占據(jù)團(tuán)簇的表面位置,而Ni原子趨于占據(jù)中心位置,即Cu-Ni混合團(tuán)簇中有明顯Cu元素偏析現(xiàn)象,Cu和Ni并沒有互相混合形成有序結(jié)構(gòu),而是分別聚集在一起,所形成的體系對稱性很低,呈割據(jù)狀態(tài)。

    [1] J.Oviedo,R.E.Palmer,J.Chem.Phys.,117(2002),9548.

    [2] J.Jellinek,E.B.Krissinel,Chem Phys,lett.258(1996),283; ib.id, Chem, Phys.lett.272(1997),30

    [3] J.Wang,F(xiàn).Ding,el.al.,Solid State Commun,119(2001),13.

    [4] M.Li and N.X.Chen,Phys.Rev.B52(2)(1995)997-1003.

    [5] M.Li and N.X.Chen,Z.Phys.,B100(2)(1996)169 - 172; M.Li, S.J.Liu, N.X.Chen,Y.G.Teng, W.X.Zhang, and S.Z.Zhang,Phys.Lett.,A169(5)(1992)364-370.

    [6] M.Li,S.J.Liu,and N.X.Chen,Phys.Lett.,A177(2)(1993)134-138.

    [7] N.X.Chen,Z.D.Chen,Y.N.Shen,S.J.Liu,M.Li,Phys.Lett.,A184(4-5)(1994)347-351.

    [8] M.Li and N.X.Chen,Phys.Rev.B52(2)(1995)997-1003.

    [9] M.Li and N.X.Chen,Z.Phys.,B100(2)(1996) 169 - 172; M.Li, S.J.Liu,N.X.Chen,Y.G.Teng,W.X.Zhang,and S.Z.Zhang, Phys.Lett., A169 (5)(1992)364-370.

    [10] M.Li,S.J.Liu,and N.X.Chen,Phys.Lett.,A177(2)(1993)134-138.

    [11] J.P.Perdew,J.A.Chevary,S.H.Vosko,K.A.Jackson,M.R.Pederson,D.J.Singh,C.Fiolhais,Phys.Rev.,B46(11)(1992)6671-6687;J.P.Perdew,Physics Letters,A165(1)(1992)79-82;J.P.Perdew, Y.Wang,Phys.Rev.B45(23)(1992)13244-13249.

    [12] J.M.Montejano-CarrizalesM.P.Iniguer and.A.Alonso,Phys.Rev.B49(1995)16649-16658.

    First-principles study of the stable structures of CuxNi19-x(x<19)bimetallic clusters

    FENG CuiJu

    (Department of basic courses,North China institute of Science And Technology,Yanjiao Beijing-East 101601)

    Using first principles total-energy calculations and the Mobius inversion method in metal Cu,Ni and Cu-Ni,We obtained the pair-potential from cohesive energy.Structural of Cu-Ni bimetallic clusters derived from the smart method are presented and discussed.In Cu-Ni clusters,the Cu atoms preferentially occupy the surface sites,whereas the Ni atoms occupy the interior of the cluster.Cu segregates to the surface in the Cu-Ni clusters.

    Chen’s 3D lattice inverse method;bimetallic clusters;stable structure

    O469

    A

    1672-7169(2011)01-0071-05

    2011-01-04?;痦?xiàng)目:華北科技學(xué)院科研基金(B09017)。

    馮翠菊(1970-),女,河北藁城人,碩士,華北科技學(xué)院基礎(chǔ)部講師,研究方向:材料的模擬與設(shè)計(jì)。

    猜你喜歡
    原子間第一性偏析
    AuBe5型新相NdMgNi4-xCox的第一性原理研究
    SO2和NO2在γ-Al2O3(110)表面吸附的第一性原理計(jì)算
    原子間相互作用勢對中Al濃度Ni75AlxV25?x合金沉淀序列的影響?
    與熱庫耦合的光學(xué)腔內(nèi)三原子間的糾纏動力學(xué)?
    團(tuán)簇Mn3BP的電子自旋密度
    W、Bi摻雜及(W、Bi)共摻銳鈦礦TiO2的第一性原理計(jì)算
    缺陷和硫摻雜黑磷的第一性原理計(jì)算
    改善60Si2MnA彈簧鋼小方坯中心碳偏析的研究
    上海金屬(2016年4期)2016-11-23 05:39:12
    原鋁旋轉(zhuǎn)偏析法凈化工藝
    改善82B 小方坯碳偏析的連鑄工藝研究
    新疆鋼鐵(2015年1期)2015-11-07 03:11:03
    久久精品国产亚洲av香蕉五月| 波多野结衣高清作品| 欧美激情在线99| 国产伦理片在线播放av一区 | 国产真实伦视频高清在线观看| 国产欧美日韩精品一区二区| 成人性生交大片免费视频hd| 亚洲欧洲日产国产| 高清午夜精品一区二区三区 | av免费观看日本| 国产亚洲精品久久久久久毛片| 国产成人精品一,二区 | 国产精品.久久久| 校园春色视频在线观看| 春色校园在线视频观看| 97超视频在线观看视频| 免费看光身美女| 日韩欧美在线乱码| 亚洲国产色片| 五月玫瑰六月丁香| 国产精品一区二区性色av| 国产高清三级在线| 亚洲国产色片| 久久久精品94久久精品| 村上凉子中文字幕在线| 99久久成人亚洲精品观看| 一级毛片久久久久久久久女| 精品一区二区三区人妻视频| 国产在线男女| 久久久久网色| 少妇丰满av| av天堂中文字幕网| 久久久久久久久久久丰满| 日日干狠狠操夜夜爽| 久久精品人妻少妇| 亚洲av成人精品一区久久| 欧美成人一区二区免费高清观看| 草草在线视频免费看| 99久久九九国产精品国产免费| 99久久人妻综合| 亚洲七黄色美女视频| 在线免费观看的www视频| 欧美激情久久久久久爽电影| 国产日本99.免费观看| 日韩强制内射视频| 夜夜看夜夜爽夜夜摸| 亚洲av成人精品一区久久| 免费看日本二区| 国产精品日韩av在线免费观看| 亚洲一区二区三区色噜噜| 精品久久久久久久久久免费视频| 一进一出抽搐动态| 欧美激情在线99| 亚洲成人中文字幕在线播放| 伦精品一区二区三区| 少妇熟女欧美另类| 色综合站精品国产| 大型黄色视频在线免费观看| 天天一区二区日本电影三级| 日韩一区二区三区影片| 亚洲av成人av| 欧美潮喷喷水| 一卡2卡三卡四卡精品乱码亚洲| 在现免费观看毛片| 免费不卡的大黄色大毛片视频在线观看 | 精品人妻熟女av久视频| 国产成人福利小说| 国产精品,欧美在线| 91aial.com中文字幕在线观看| 国产精品99久久久久久久久| 亚洲欧美日韩高清在线视频| 免费看日本二区| 三级毛片av免费| а√天堂www在线а√下载| 中文精品一卡2卡3卡4更新| 欧美日韩在线观看h| 夜夜看夜夜爽夜夜摸| 欧美高清性xxxxhd video| 春色校园在线视频观看| 两个人视频免费观看高清| 久久久精品94久久精品| 国产精品免费一区二区三区在线| 精品人妻偷拍中文字幕| 国产高清不卡午夜福利| 久久亚洲国产成人精品v| 日韩大尺度精品在线看网址| 国产精品美女特级片免费视频播放器| 久久这里只有精品中国| 久久99精品国语久久久| 欧美不卡视频在线免费观看| 夜夜夜夜夜久久久久| 我要搜黄色片| 久久精品久久久久久噜噜老黄 | 乱系列少妇在线播放| www.色视频.com| 女同久久另类99精品国产91| 精品久久久久久久久久久久久| 嫩草影院入口| 一级毛片电影观看 | 婷婷色av中文字幕| 岛国在线免费视频观看| 亚洲国产欧洲综合997久久,| 男人舔奶头视频| 国产黄片美女视频| 成人综合一区亚洲| 性欧美人与动物交配| 亚洲精品影视一区二区三区av| 日日撸夜夜添| 人体艺术视频欧美日本| 国产色婷婷99| 热99re8久久精品国产| 一级二级三级毛片免费看| 日本撒尿小便嘘嘘汇集6| 亚洲四区av| 日韩强制内射视频| 国产精品精品国产色婷婷| 夜夜看夜夜爽夜夜摸| 亚洲,欧美,日韩| 大又大粗又爽又黄少妇毛片口| 一卡2卡三卡四卡精品乱码亚洲| 国产片特级美女逼逼视频| 听说在线观看完整版免费高清| 国产一区二区在线观看日韩| 亚洲中文字幕一区二区三区有码在线看| 亚洲av中文字字幕乱码综合| 久久人妻av系列| 亚洲成人久久爱视频| 嫩草影院新地址| 激情 狠狠 欧美| 高清日韩中文字幕在线| 自拍偷自拍亚洲精品老妇| 边亲边吃奶的免费视频| 亚洲第一区二区三区不卡| 日日摸夜夜添夜夜添av毛片| 最近最新中文字幕大全电影3| 嘟嘟电影网在线观看| 51国产日韩欧美| 简卡轻食公司| 禁无遮挡网站| 97热精品久久久久久| 最近中文字幕高清免费大全6| or卡值多少钱| a级毛片免费高清观看在线播放| 狂野欧美激情性xxxx在线观看| 亚洲内射少妇av| 国产成年人精品一区二区| 国产精品久久久久久精品电影| 成人三级黄色视频| 国产精品福利在线免费观看| 国产av在哪里看| 成人高潮视频无遮挡免费网站| 午夜激情福利司机影院| 国产高清视频在线观看网站| 亚洲av二区三区四区| 99热全是精品| 久久久久久久久久黄片| 少妇熟女aⅴ在线视频| 亚洲成人久久爱视频| 国产一区亚洲一区在线观看| 变态另类成人亚洲欧美熟女| 国产单亲对白刺激| 91久久精品国产一区二区三区| 一个人看的www免费观看视频| 能在线免费观看的黄片| 久久精品国产清高在天天线| 国产成人freesex在线| 一个人免费在线观看电影| 神马国产精品三级电影在线观看| 91av网一区二区| 哪个播放器可以免费观看大片| 最新中文字幕久久久久| 精品久久久久久久久av| 少妇猛男粗大的猛烈进出视频 | 高清毛片免费看| 日本成人三级电影网站| 午夜福利在线在线| 午夜福利在线观看吧| av黄色大香蕉| ponron亚洲| 观看美女的网站| 欧美高清性xxxxhd video| 三级男女做爰猛烈吃奶摸视频| 波多野结衣高清作品| 99热全是精品| 欧美不卡视频在线免费观看| 成人特级av手机在线观看| 性欧美人与动物交配| 夫妻性生交免费视频一级片| 亚洲无线在线观看| 别揉我奶头 嗯啊视频| 我的老师免费观看完整版| 长腿黑丝高跟| 如何舔出高潮| 精品人妻偷拍中文字幕| videossex国产| 啦啦啦观看免费观看视频高清| av天堂在线播放| 女人十人毛片免费观看3o分钟| av黄色大香蕉| 少妇人妻一区二区三区视频| 美女cb高潮喷水在线观看| 天堂av国产一区二区熟女人妻| 97人妻精品一区二区三区麻豆| 99riav亚洲国产免费| 小蜜桃在线观看免费完整版高清| 日韩欧美一区二区三区在线观看| 最近2019中文字幕mv第一页| av天堂中文字幕网| 99九九线精品视频在线观看视频| 免费黄网站久久成人精品| 国内精品美女久久久久久| 免费大片18禁| 欧美潮喷喷水| 国产成人91sexporn| 欧美xxxx性猛交bbbb| 国产探花极品一区二区| 亚洲国产欧美人成| 直男gayav资源| 大型黄色视频在线免费观看| 一本久久精品| 久久久久国产网址| a级毛片a级免费在线| 精品日产1卡2卡| 亚洲一区二区三区色噜噜| 青青草视频在线视频观看| 乱码一卡2卡4卡精品| 亚洲精品久久久久久婷婷小说 | 久久九九热精品免费| 国产成人午夜福利电影在线观看| 日韩三级伦理在线观看| 桃色一区二区三区在线观看| 悠悠久久av| 校园人妻丝袜中文字幕| 看片在线看免费视频| 日韩中字成人| 简卡轻食公司| 全区人妻精品视频| 69人妻影院| 国产成人a区在线观看| 免费看a级黄色片| 一边摸一边抽搐一进一小说| 精品少妇黑人巨大在线播放 | 性色avwww在线观看| 日韩成人伦理影院| 长腿黑丝高跟| 精品99又大又爽又粗少妇毛片| 边亲边吃奶的免费视频| 日韩高清综合在线| 国产成人freesex在线| 国产成人精品婷婷| 精品久久久久久成人av| 亚洲欧美精品专区久久| 男人和女人高潮做爰伦理| 亚洲欧美日韩东京热| 亚洲av成人av| 中文字幕人妻熟人妻熟丝袜美| 2021天堂中文幕一二区在线观| 日韩av在线大香蕉| 成人午夜精彩视频在线观看| 草草在线视频免费看| 一个人观看的视频www高清免费观看| 特大巨黑吊av在线直播| 成人亚洲欧美一区二区av| 九九久久精品国产亚洲av麻豆| 国产精品久久久久久av不卡| 中文字幕制服av| 麻豆成人av视频| 1000部很黄的大片| 色吧在线观看| 成人av在线播放网站| 男女那种视频在线观看| 日韩一本色道免费dvd| 精品久久久久久久久av| 国产69精品久久久久777片| 免费一级毛片在线播放高清视频| 国产极品精品免费视频能看的| 国产精品一区二区三区四区久久| 亚洲欧美日韩高清专用| 中国美白少妇内射xxxbb| 黄片无遮挡物在线观看| 亚洲av熟女| 人妻久久中文字幕网| 成人毛片60女人毛片免费| 免费人成视频x8x8入口观看| 免费看光身美女| 长腿黑丝高跟| 亚洲国产精品久久男人天堂| 哪个播放器可以免费观看大片| 国产黄色视频一区二区在线观看 | 草草在线视频免费看| 国产麻豆成人av免费视频| 国产老妇女一区| 一级毛片久久久久久久久女| 日韩在线高清观看一区二区三区| 久久人人爽人人爽人人片va| 国产 一区 欧美 日韩| 亚洲精品日韩av片在线观看| 啦啦啦啦在线视频资源| 日本五十路高清| 男人的好看免费观看在线视频| 日本熟妇午夜| 欧美一区二区亚洲| 日韩成人伦理影院| 一级毛片电影观看 | 国产精品人妻久久久久久| 亚洲精品色激情综合| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品成人久久久久久| 最近的中文字幕免费完整| 夜夜看夜夜爽夜夜摸| 久久韩国三级中文字幕| 日韩,欧美,国产一区二区三区 | 日韩欧美三级三区| 欧洲精品卡2卡3卡4卡5卡区| 午夜a级毛片| 久久精品久久久久久久性| 边亲边吃奶的免费视频| 欧美色视频一区免费| 美女脱内裤让男人舔精品视频 | 久久这里有精品视频免费| 久久精品国产清高在天天线| 尤物成人国产欧美一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 青春草国产在线视频 | 好男人视频免费观看在线| 伊人久久精品亚洲午夜| 亚洲av一区综合| 欧美高清成人免费视频www| 丰满人妻一区二区三区视频av| 欧美高清性xxxxhd video| 91av网一区二区| 一级黄片播放器| 久久久久久九九精品二区国产| 久久国内精品自在自线图片| 成年女人看的毛片在线观看| 国产成人影院久久av| 白带黄色成豆腐渣| 日本av手机在线免费观看| 日韩欧美在线乱码| 九色成人免费人妻av| 99热这里只有是精品50| 久久久久久久久大av| 嘟嘟电影网在线观看| 2021天堂中文幕一二区在线观| 最近的中文字幕免费完整| 一区二区三区高清视频在线| 最后的刺客免费高清国语| 亚洲国产精品合色在线| 丝袜美腿在线中文| 九草在线视频观看| 天美传媒精品一区二区| 欧美人与善性xxx| 日韩强制内射视频| 国产亚洲精品av在线| 亚洲精品国产成人久久av| 大又大粗又爽又黄少妇毛片口| 91麻豆精品激情在线观看国产| 99久久人妻综合| 亚洲欧洲国产日韩| 人妻久久中文字幕网| 哪个播放器可以免费观看大片| 国产视频内射| 久久人人精品亚洲av| 一级黄片播放器| 国产三级中文精品| 国产中年淑女户外野战色| 欧美一区二区精品小视频在线| 99在线视频只有这里精品首页| 九草在线视频观看| 卡戴珊不雅视频在线播放| 天天一区二区日本电影三级| 一个人免费在线观看电影| 亚洲,欧美,日韩| 欧美一区二区国产精品久久精品| 国产午夜精品一二区理论片| 一边亲一边摸免费视频| 免费人成在线观看视频色| 国产精品人妻久久久影院| 国产综合懂色| 亚州av有码| 99热这里只有是精品50| 亚洲内射少妇av| 中文字幕久久专区| 久久久久久国产a免费观看| 青春草视频在线免费观看| 长腿黑丝高跟| 国产av麻豆久久久久久久| 国产成人精品婷婷| 啦啦啦观看免费观看视频高清| 黑人高潮一二区| 国产成人a区在线观看| 97超视频在线观看视频| 天堂网av新在线| av又黄又爽大尺度在线免费看 | 精品国内亚洲2022精品成人| 国产一区二区在线av高清观看| 日韩视频在线欧美| 色哟哟·www| 精品国内亚洲2022精品成人| 成年女人永久免费观看视频| 国产精品不卡视频一区二区| 国产黄片视频在线免费观看| 老熟妇乱子伦视频在线观看| 亚洲婷婷狠狠爱综合网| 国产熟女欧美一区二区| 精品欧美国产一区二区三| 边亲边吃奶的免费视频| 国产黄色视频一区二区在线观看 | 丰满乱子伦码专区| 啦啦啦观看免费观看视频高清| 亚洲国产欧洲综合997久久,| av专区在线播放| 久久久久久国产a免费观看| 看黄色毛片网站| 黄色视频,在线免费观看| 午夜免费男女啪啪视频观看| 一个人免费在线观看电影| 亚洲国产精品合色在线| 亚洲五月天丁香| 精品不卡国产一区二区三区| 久久久久九九精品影院| 91久久精品电影网| 国产精品一二三区在线看| 91久久精品电影网| 久久人人爽人人片av| 日韩一本色道免费dvd| 国产日本99.免费观看| 国产三级在线视频| 精品日产1卡2卡| 村上凉子中文字幕在线| 久久草成人影院| 久久精品夜夜夜夜夜久久蜜豆| 免费黄网站久久成人精品| 亚洲综合色惰| 99精品在免费线老司机午夜| 淫秽高清视频在线观看| 国产精品乱码一区二三区的特点| 青青草视频在线视频观看| 亚洲国产精品合色在线| 插阴视频在线观看视频| 全区人妻精品视频| 可以在线观看的亚洲视频| 亚洲精品自拍成人| a级毛色黄片| 欧美日韩精品成人综合77777| 欧美极品一区二区三区四区| 免费av毛片视频| 亚洲在线观看片| 国产亚洲91精品色在线| 中文字幕av在线有码专区| 色哟哟哟哟哟哟| 国内精品宾馆在线| 国产三级中文精品| 久久久国产成人精品二区| 精华霜和精华液先用哪个| 中文字幕免费在线视频6| 久久久久久久久大av| 日韩大尺度精品在线看网址| 综合色av麻豆| 麻豆成人av视频| 69人妻影院| 一级黄色大片毛片| 久久久色成人| 国产三级中文精品| 国产精品伦人一区二区| 亚洲国产精品sss在线观看| 亚洲久久久久久中文字幕| 久久久久久久久大av| 日韩大尺度精品在线看网址| 嫩草影院精品99| 国内精品一区二区在线观看| 一级毛片我不卡| 人妻夜夜爽99麻豆av| 99热6这里只有精品| 久久99精品国语久久久| 日本三级黄在线观看| 国产精品免费一区二区三区在线| 国产探花在线观看一区二区| 偷拍熟女少妇极品色| 不卡视频在线观看欧美| 日日摸夜夜添夜夜爱| 精品不卡国产一区二区三区| 国产淫片久久久久久久久| 亚洲精品日韩在线中文字幕 | 99热网站在线观看| 亚洲三级黄色毛片| 亚洲av一区综合| 青春草国产在线视频 | 在线免费观看不下载黄p国产| 美女大奶头视频| 久久精品国产鲁丝片午夜精品| 亚洲精品色激情综合| 国产午夜精品一二区理论片| 亚洲国产精品成人综合色| 亚洲中文字幕日韩| 午夜激情福利司机影院| 国产一区二区在线观看日韩| 又爽又黄a免费视频| 亚洲欧洲日产国产| 一个人看视频在线观看www免费| 亚洲aⅴ乱码一区二区在线播放| 大香蕉久久网| 亚洲av免费高清在线观看| 欧美又色又爽又黄视频| 99热这里只有精品一区| 深夜精品福利| 国产成人午夜福利电影在线观看| 亚洲欧美日韩卡通动漫| 我的老师免费观看完整版| 久久精品国产亚洲网站| 男女做爰动态图高潮gif福利片| 大又大粗又爽又黄少妇毛片口| 我要看日韩黄色一级片| 成人亚洲欧美一区二区av| 亚洲欧美日韩高清专用| 免费看a级黄色片| 在线播放国产精品三级| 亚洲精品粉嫩美女一区| 尾随美女入室| 一个人观看的视频www高清免费观看| 亚洲一区二区三区色噜噜| a级毛色黄片| 国产一区二区三区av在线 | 亚洲精品粉嫩美女一区| 亚洲最大成人中文| 美女cb高潮喷水在线观看| 蜜臀久久99精品久久宅男| 日本在线视频免费播放| 床上黄色一级片| 亚洲国产精品合色在线| 麻豆国产av国片精品| 亚洲自拍偷在线| 亚洲国产精品国产精品| 少妇的逼好多水| 亚洲最大成人中文| 男人和女人高潮做爰伦理| 波多野结衣高清无吗| 欧美高清性xxxxhd video| 床上黄色一级片| 岛国毛片在线播放| 亚洲国产色片| 少妇丰满av| 久久久精品大字幕| 一级二级三级毛片免费看| 99视频精品全部免费 在线| 国产黄色视频一区二区在线观看 | 久久韩国三级中文字幕| ponron亚洲| 人妻久久中文字幕网| 变态另类丝袜制服| 大香蕉久久网| 人妻少妇偷人精品九色| 国产高清激情床上av| 国产精品野战在线观看| 女人十人毛片免费观看3o分钟| 99久久久亚洲精品蜜臀av| 女人被狂操c到高潮| 美女黄网站色视频| 亚洲,欧美,日韩| 亚洲精品国产av成人精品| 中文欧美无线码| 国产精品麻豆人妻色哟哟久久 | 色吧在线观看| 欧美人与善性xxx| 男人和女人高潮做爰伦理| 久久99蜜桃精品久久| 国产 一区 欧美 日韩| 一级毛片久久久久久久久女| 精品熟女少妇av免费看| 91精品国产九色| 中文精品一卡2卡3卡4更新| 欧美人与善性xxx| 免费观看的影片在线观看| 97人妻精品一区二区三区麻豆| 欧美一级a爱片免费观看看| 一级毛片久久久久久久久女| 日本黄大片高清| 麻豆av噜噜一区二区三区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美日本亚洲视频在线播放| 国产乱人偷精品视频| 国产视频首页在线观看| 中文亚洲av片在线观看爽| 成人二区视频| 日韩三级伦理在线观看| 亚洲国产精品合色在线| 午夜爱爱视频在线播放| 日韩欧美一区二区三区在线观看| 91精品一卡2卡3卡4卡| 亚洲va在线va天堂va国产| 精品人妻视频免费看| 国产日韩欧美在线精品| 亚洲七黄色美女视频| 大型黄色视频在线免费观看| 精品久久久久久久久久久久久| 免费大片18禁| 身体一侧抽搐| 蜜桃久久精品国产亚洲av| 嫩草影院精品99| 成人特级黄色片久久久久久久| 韩国av在线不卡| 国产成年人精品一区二区| 黄色一级大片看看| 久99久视频精品免费| 久久精品国产自在天天线| 看黄色毛片网站| 日本一二三区视频观看| 精品人妻视频免费看| 欧美潮喷喷水| 午夜老司机福利剧场| 久久99蜜桃精品久久| 欧美最黄视频在线播放免费| 一个人观看的视频www高清免费观看| 日韩av不卡免费在线播放| 国产色爽女视频免费观看|