• <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
    午夜免费男女啪啪视频观看 | 乱码一卡2卡4卡精品| av在线天堂中文字幕| 欧美不卡视频在线免费观看| 久久6这里有精品| 精品福利观看| av女优亚洲男人天堂| 深夜a级毛片| 少妇高潮的动态图| 国产激情偷乱视频一区二区| 一级黄色大片毛片| 久久精品国产亚洲av天美| 观看美女的网站| 男女下面进入的视频免费午夜| 精品一区二区三区视频在线| 日韩欧美一区二区三区在线观看| 午夜免费男女啪啪视频观看 | 久久人人精品亚洲av| 看片在线看免费视频| 亚洲人与动物交配视频| 国模一区二区三区四区视频| 亚洲美女视频黄频| 九九在线视频观看精品| 美女高潮的动态| 欧美激情在线99| 久久精品夜夜夜夜夜久久蜜豆| 亚洲精品乱码久久久v下载方式| 极品教师在线免费播放| 99热网站在线观看| av.在线天堂| 级片在线观看| 男女之事视频高清在线观看| 99久久久亚洲精品蜜臀av| 中文字幕免费在线视频6| 真人一进一出gif抽搐免费| av视频在线观看入口| 婷婷六月久久综合丁香| 精品人妻一区二区三区麻豆 | 亚洲,欧美,日韩| 国产高清三级在线| 久久亚洲真实| av天堂在线播放| 亚洲三级黄色毛片| 欧美一区二区亚洲| 国产高潮美女av| 男人狂女人下面高潮的视频| 国内揄拍国产精品人妻在线| 中文字幕av成人在线电影| 亚洲精品一卡2卡三卡4卡5卡| 亚洲久久久久久中文字幕| 在线免费十八禁| 乱人视频在线观看| 国产免费一级a男人的天堂| 黄色配什么色好看| 成人二区视频| 国产精品精品国产色婷婷| 成人国产一区最新在线观看| 国产精品98久久久久久宅男小说| 人妻丰满熟妇av一区二区三区| 午夜精品久久久久久毛片777| 欧美一区二区亚洲| 长腿黑丝高跟| 91久久精品电影网| 黄色欧美视频在线观看| 国内精品宾馆在线| 麻豆国产97在线/欧美| 亚洲av中文av极速乱 | 国产熟女欧美一区二区| 免费在线观看日本一区| 久久99热6这里只有精品| 毛片一级片免费看久久久久 | 精品免费久久久久久久清纯| 91久久精品国产一区二区三区| 一进一出抽搐gif免费好疼| 成人国产综合亚洲| 又爽又黄a免费视频| 国产精品嫩草影院av在线观看 | 国产成人福利小说| 男女啪啪激烈高潮av片| 精品不卡国产一区二区三区| 伊人久久精品亚洲午夜| 伦精品一区二区三区| 又黄又爽又刺激的免费视频.| 日本一二三区视频观看| 美女xxoo啪啪120秒动态图| 久久精品国产自在天天线| 中文字幕免费在线视频6| a级毛片a级免费在线| 国产熟女欧美一区二区| 人妻丰满熟妇av一区二区三区| 免费搜索国产男女视频| 男女做爰动态图高潮gif福利片| 不卡视频在线观看欧美| 日本一本二区三区精品| 免费看日本二区| 免费搜索国产男女视频| 精品国产三级普通话版| 国产男靠女视频免费网站| 两个人的视频大全免费| 看十八女毛片水多多多| 国产精品野战在线观看| 一边摸一边抽搐一进一小说| 久久久久久大精品| 精品久久久久久久久亚洲 | 亚洲成人中文字幕在线播放| 日本五十路高清| 一级a爱片免费观看的视频| 日韩欧美一区二区三区在线观看| 热99在线观看视频| 99热只有精品国产| 成年女人永久免费观看视频| 俄罗斯特黄特色一大片| 亚洲精品粉嫩美女一区| 99热网站在线观看| 亚洲熟妇中文字幕五十中出| 我的女老师完整版在线观看| 美女高潮喷水抽搐中文字幕| 色哟哟·www| 成人特级av手机在线观看| 听说在线观看完整版免费高清| 久久人人爽人人爽人人片va| 在线观看美女被高潮喷水网站| 男女之事视频高清在线观看| 欧美一区二区亚洲| 成人午夜高清在线视频| 国产伦人伦偷精品视频| 日本精品一区二区三区蜜桃| 草草在线视频免费看| 3wmmmm亚洲av在线观看| 深夜a级毛片| 黄片wwwwww| 久久这里只有精品中国| a级毛片a级免费在线| 欧美在线一区亚洲| 岛国在线免费视频观看| 精品人妻熟女av久视频| 一卡2卡三卡四卡精品乱码亚洲| 狂野欧美激情性xxxx在线观看| 亚洲精品在线观看二区| 99久久精品热视频| 日本 欧美在线| 如何舔出高潮| 日本熟妇午夜| 全区人妻精品视频| 久久精品国产清高在天天线| 性色avwww在线观看| 国产69精品久久久久777片| 最近最新中文字幕大全电影3| 国产成人一区二区在线| 日本-黄色视频高清免费观看| 亚洲精品一区av在线观看| 天堂√8在线中文| 一边摸一边抽搐一进一小说| 啦啦啦啦在线视频资源| 搞女人的毛片| 欧美+亚洲+日韩+国产| 欧美精品国产亚洲| 春色校园在线视频观看| 国产精品伦人一区二区| 日本熟妇午夜| 久久国产乱子免费精品| 男女之事视频高清在线观看| 最新在线观看一区二区三区| 欧美绝顶高潮抽搐喷水| 国产精品98久久久久久宅男小说| 日韩精品中文字幕看吧| 日本与韩国留学比较| 久久久成人免费电影| 69av精品久久久久久| 亚洲精品456在线播放app | 国产男靠女视频免费网站| 亚洲国产欧美人成| 午夜福利在线观看吧| 色尼玛亚洲综合影院| 成人永久免费在线观看视频| 又粗又爽又猛毛片免费看| 久久久久久大精品| 中文字幕精品亚洲无线码一区| 亚洲不卡免费看| 免费看光身美女| 欧美bdsm另类| 天天躁日日操中文字幕| 精品人妻一区二区三区麻豆 | 国产又黄又爽又无遮挡在线| 亚洲欧美日韩东京热| 日本与韩国留学比较| 亚洲熟妇熟女久久| 亚洲av免费高清在线观看| 亚洲午夜理论影院| 亚洲va日本ⅴa欧美va伊人久久| 中文亚洲av片在线观看爽| 美女大奶头视频| 午夜亚洲福利在线播放| 51国产日韩欧美| 丰满人妻一区二区三区视频av| 日本五十路高清| 色吧在线观看| 女生性感内裤真人,穿戴方法视频| 天美传媒精品一区二区| 亚洲四区av| 欧美xxxx性猛交bbbb| 国产精品伦人一区二区| 中文字幕免费在线视频6| 99热这里只有是精品50| 内射极品少妇av片p| av福利片在线观看| 小说图片视频综合网站| 不卡一级毛片| 免费看a级黄色片| 国产欧美日韩精品亚洲av| 深夜精品福利| 亚洲内射少妇av| 亚洲综合色惰| 色吧在线观看| 床上黄色一级片| 精品国产三级普通话版| 国产一区二区亚洲精品在线观看| 久久精品国产清高在天天线| 亚洲av一区综合| 18禁裸乳无遮挡免费网站照片| a在线观看视频网站| 国产精品久久久久久久久免| 亚洲在线观看片| 国产探花极品一区二区| av在线蜜桃| 日韩精品有码人妻一区| 国产一区二区亚洲精品在线观看| 精品免费久久久久久久清纯| 精品久久久久久久末码| 一边摸一边抽搐一进一小说| 亚洲国产色片| 国产麻豆成人av免费视频| 可以在线观看毛片的网站| 熟妇人妻久久中文字幕3abv| 日韩精品中文字幕看吧| 一区二区三区激情视频| 波多野结衣高清无吗| www日本黄色视频网| or卡值多少钱| 简卡轻食公司| 美女 人体艺术 gogo| 欧美3d第一页| 99热这里只有精品一区| 99热这里只有是精品50| 国产av麻豆久久久久久久| 99热精品在线国产| 日韩大尺度精品在线看网址| 天堂av国产一区二区熟女人妻| 成人综合一区亚洲| 国产亚洲精品综合一区在线观看| 亚洲中文日韩欧美视频| 中国美白少妇内射xxxbb| 久久久久久久亚洲中文字幕| 国产在线男女| 内射极品少妇av片p| 亚洲最大成人中文| 网址你懂的国产日韩在线| 日韩中文字幕欧美一区二区| 国产女主播在线喷水免费视频网站 | 国产又黄又爽又无遮挡在线| 人妻丰满熟妇av一区二区三区| 亚洲精品乱码久久久v下载方式| 免费看日本二区| 一本精品99久久精品77| 国产精品人妻久久久久久| 国产精品综合久久久久久久免费| 色噜噜av男人的天堂激情| 久久久久久久精品吃奶| 老司机午夜福利在线观看视频| 18+在线观看网站| 婷婷精品国产亚洲av在线| 久久久精品欧美日韩精品| 非洲黑人性xxxx精品又粗又长| 97超级碰碰碰精品色视频在线观看| 在线看三级毛片| 又爽又黄a免费视频| 久久久久免费精品人妻一区二区| 成人永久免费在线观看视频| 免费无遮挡裸体视频| 国产精品一区二区三区四区免费观看 | 性色avwww在线观看| 国产一区二区在线av高清观看| 精品人妻1区二区| 少妇熟女aⅴ在线视频| 国产伦精品一区二区三区四那| 老司机午夜福利在线观看视频| 最近中文字幕高清免费大全6 | 久久精品夜夜夜夜夜久久蜜豆| av天堂在线播放| 伊人久久精品亚洲午夜| 成人一区二区视频在线观看| 亚洲欧美精品综合久久99| 男插女下体视频免费在线播放| 国产一区二区亚洲精品在线观看| 午夜福利在线在线| 波野结衣二区三区在线| 亚洲av中文字字幕乱码综合| 亚洲图色成人| 俺也久久电影网| 日本a在线网址| 在线观看午夜福利视频| 中文字幕精品亚洲无线码一区| 如何舔出高潮| 精品一区二区三区av网在线观看| 久久久精品大字幕| 亚洲专区国产一区二区| 2021天堂中文幕一二区在线观| 免费电影在线观看免费观看| 91久久精品国产一区二区三区| 日本免费a在线| 成年人黄色毛片网站| 色av中文字幕| 色5月婷婷丁香| 色哟哟哟哟哟哟| 热99在线观看视频| 亚洲av日韩精品久久久久久密| 男女之事视频高清在线观看| 18禁黄网站禁片午夜丰满| 亚洲男人的天堂狠狠| 国产午夜精品论理片| 三级国产精品欧美在线观看| 久久精品国产自在天天线| 欧美黑人欧美精品刺激| 欧美一区二区亚洲| 美女cb高潮喷水在线观看| 午夜福利成人在线免费观看| 免费人成在线观看视频色| 搞女人的毛片| 成人av一区二区三区在线看| 国产精品av视频在线免费观看| 精品午夜福利在线看| 一夜夜www| 在线免费十八禁| 亚洲不卡免费看| 欧美潮喷喷水| 欧美最新免费一区二区三区| 亚洲乱码一区二区免费版| 99热只有精品国产| 亚洲va在线va天堂va国产| 两人在一起打扑克的视频| 精品久久久久久久久久久久久| 淫妇啪啪啪对白视频| 欧美激情国产日韩精品一区| 日本-黄色视频高清免费观看| 午夜免费成人在线视频| 国产成年人精品一区二区| 性插视频无遮挡在线免费观看| 国产精品嫩草影院av在线观看 | 麻豆av噜噜一区二区三区| 免费在线观看日本一区| 国产精品久久久久久久电影| 色综合亚洲欧美另类图片| 免费人成视频x8x8入口观看| 亚洲av免费高清在线观看| 国内久久婷婷六月综合欲色啪| 小蜜桃在线观看免费完整版高清| 亚洲va日本ⅴa欧美va伊人久久| 99视频精品全部免费 在线| 精品久久久久久久久亚洲 | 九九爱精品视频在线观看| 91狼人影院| 成人毛片a级毛片在线播放| 97超视频在线观看视频| 日本爱情动作片www.在线观看 | 免费看美女性在线毛片视频| 毛片女人毛片| 亚洲avbb在线观看| 日韩一区二区视频免费看| 女同久久另类99精品国产91| 丰满的人妻完整版| 久久久久九九精品影院| 国产色爽女视频免费观看| 22中文网久久字幕| 亚洲精品色激情综合| 一区二区三区免费毛片| 久99久视频精品免费| 性插视频无遮挡在线免费观看| 69人妻影院| 一本久久中文字幕| 高清毛片免费观看视频网站| 狂野欧美白嫩少妇大欣赏| 欧美中文日本在线观看视频| 欧美日韩亚洲国产一区二区在线观看| 免费看av在线观看网站| 观看免费一级毛片| 国产精品99久久久久久久久| av.在线天堂| 国产三级中文精品| 丰满的人妻完整版| 亚洲精品色激情综合| 男女边吃奶边做爰视频| 色播亚洲综合网| 国产精品一区www在线观看 | 啦啦啦观看免费观看视频高清| 成年女人看的毛片在线观看| 成人性生交大片免费视频hd| 国产色婷婷99| 成人av一区二区三区在线看| a级毛片免费高清观看在线播放| 久久精品影院6| 国产欧美日韩一区二区精品| 99久久精品一区二区三区| 亚洲专区国产一区二区| 久久99热这里只有精品18| 男插女下体视频免费在线播放| 久久久久久久久久久丰满 | 久久精品人妻少妇| 国产精品自产拍在线观看55亚洲| 人妻少妇偷人精品九色| 韩国av一区二区三区四区| 亚洲美女黄片视频| 欧美激情在线99| 欧美极品一区二区三区四区| 伦理电影大哥的女人| 美女黄网站色视频| 人妻少妇偷人精品九色| 欧美日本视频| 精华霜和精华液先用哪个| 欧美成人一区二区免费高清观看| 亚洲精品乱码久久久v下载方式| 国内精品一区二区在线观看| 色噜噜av男人的天堂激情| 午夜福利在线观看吧| 久久久久久久亚洲中文字幕| 久久久久久九九精品二区国产| 九色成人免费人妻av| 麻豆国产av国片精品| 国产探花在线观看一区二区| 永久网站在线| 男女边吃奶边做爰视频| 日本黄色片子视频| 亚洲精品一卡2卡三卡4卡5卡| 搡女人真爽免费视频火全软件 | 婷婷精品国产亚洲av在线| 欧美成人性av电影在线观看| 日本在线视频免费播放| 三级国产精品欧美在线观看| 亚洲狠狠婷婷综合久久图片| 国产高清三级在线| 一区福利在线观看| 嫩草影院新地址| 欧美高清性xxxxhd video| 少妇高潮的动态图| 免费人成在线观看视频色| 黄色一级大片看看| 色吧在线观看| 51国产日韩欧美| 国产老妇女一区| 日本爱情动作片www.在线观看 | 国产69精品久久久久777片| 免费av不卡在线播放| 亚洲狠狠婷婷综合久久图片| 成年女人毛片免费观看观看9| 九九在线视频观看精品| 18禁在线播放成人免费| av黄色大香蕉| 国产男靠女视频免费网站| 春色校园在线视频观看| 欧美绝顶高潮抽搐喷水| 亚洲内射少妇av| 亚洲第一电影网av| 久久欧美精品欧美久久欧美| 欧美黑人巨大hd| 日日摸夜夜添夜夜添av毛片 | x7x7x7水蜜桃| 亚洲专区国产一区二区| 极品教师在线免费播放| 国产精品一区二区免费欧美| 免费看a级黄色片| 国产淫片久久久久久久久| 国产毛片a区久久久久| 国产av不卡久久| 日韩欧美在线二视频| 又粗又爽又猛毛片免费看| 精品久久久久久久末码| 在线国产一区二区在线| 国产一区二区在线观看日韩| 亚洲国产色片| 精品福利观看| 制服丝袜大香蕉在线| 国产探花在线观看一区二区| 一区福利在线观看| 无人区码免费观看不卡| 亚洲人成网站在线播放欧美日韩| 国产精品一区二区三区四区久久| 乱系列少妇在线播放| 变态另类成人亚洲欧美熟女| 国产精品女同一区二区软件 | 一进一出抽搐gif免费好疼| 久久亚洲真实| 少妇的逼水好多| 97超视频在线观看视频| 中文字幕免费在线视频6| 最好的美女福利视频网| www.www免费av| 不卡视频在线观看欧美| 琪琪午夜伦伦电影理论片6080| 亚洲美女搞黄在线观看 | 亚洲专区中文字幕在线| 99久久无色码亚洲精品果冻| 国产亚洲精品综合一区在线观看| 99热网站在线观看| 欧美成人免费av一区二区三区| 999久久久精品免费观看国产| 极品教师在线免费播放| 十八禁国产超污无遮挡网站| 国产不卡一卡二| 99热6这里只有精品| 女同久久另类99精品国产91| 亚洲性夜色夜夜综合| 国产在视频线在精品| 亚洲美女视频黄频| 久久99热6这里只有精品| 亚洲一区二区三区色噜噜| 少妇猛男粗大的猛烈进出视频 | 久久久久久久久久成人| 亚洲av免费在线观看| 精品日产1卡2卡| 婷婷精品国产亚洲av在线| 午夜福利视频1000在线观看| 黄色一级大片看看| 听说在线观看完整版免费高清| 最近最新中文字幕大全电影3| 精品久久久久久,| 别揉我奶头~嗯~啊~动态视频| 99久久九九国产精品国产免费| 九色成人免费人妻av| 久久精品国产亚洲av天美| 国产一区二区三区av在线 | 国产精品亚洲一级av第二区| 校园人妻丝袜中文字幕| 亚洲中文日韩欧美视频| 亚洲av成人av| 国产aⅴ精品一区二区三区波| 国产三级在线视频| 国产人妻一区二区三区在| 精品不卡国产一区二区三区| 69av精品久久久久久| 国产毛片a区久久久久| 狠狠狠狠99中文字幕| 特级一级黄色大片| 欧美性猛交黑人性爽| 成人综合一区亚洲| 亚洲综合色惰| 亚洲国产精品久久男人天堂| eeuss影院久久| 如何舔出高潮| 3wmmmm亚洲av在线观看| xxxwww97欧美| 国产三级中文精品| 成人av在线播放网站| 色哟哟哟哟哟哟| 99久久精品国产国产毛片| 一级黄色大片毛片| 亚洲精华国产精华精| 国产精品国产高清国产av| 俄罗斯特黄特色一大片| 日日摸夜夜添夜夜添小说| 亚洲精品色激情综合| 久久久久精品国产欧美久久久| 此物有八面人人有两片| 欧美色视频一区免费| 欧美日韩综合久久久久久 | 97热精品久久久久久| 免费无遮挡裸体视频| 岛国在线免费视频观看| 内射极品少妇av片p| 欧美一区二区亚洲| 欧美激情国产日韩精品一区| 日本a在线网址| 免费观看精品视频网站| 国产精品av视频在线免费观看| 国产精品久久电影中文字幕| 久久中文看片网| 在线播放国产精品三级| 男人舔女人下体高潮全视频| 免费看日本二区| 久久国产乱子免费精品| 夜夜夜夜夜久久久久| 91久久精品电影网| 欧美色视频一区免费| 日本-黄色视频高清免费观看| 久久久久久久精品吃奶| 搞女人的毛片| 欧美日韩乱码在线| 亚洲性久久影院| 超碰av人人做人人爽久久| 久久久久免费精品人妻一区二区| 亚洲av成人av| 国产精品一区二区免费欧美| 我的女老师完整版在线观看| 级片在线观看| 91麻豆精品激情在线观看国产| 精品日产1卡2卡| av国产免费在线观看| 婷婷丁香在线五月| 日韩欧美精品v在线| 国产69精品久久久久777片| 亚洲真实伦在线观看| 午夜激情福利司机影院| 欧美+亚洲+日韩+国产| 国产成人aa在线观看| 22中文网久久字幕| 熟女人妻精品中文字幕| 欧美一区二区精品小视频在线| 桃红色精品国产亚洲av| 中文亚洲av片在线观看爽| 春色校园在线视频观看| 永久网站在线| 久久久国产成人免费| 国产高清视频在线播放一区| 久99久视频精品免费| 97碰自拍视频| 桃色一区二区三区在线观看| 少妇人妻一区二区三区视频|