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

    合金納米粒子生長形態(tài)的分子動力學(xué)模擬

    2015-12-20 06:57:48余藝程慶軍郭詩奇楊劍瑜
    關(guān)鍵詞:基底原子合金

    余藝, 程慶軍, 郭詩奇, 楊劍瑜

    ?

    合金納米粒子生長形態(tài)的分子動力學(xué)模擬

    余藝, 程慶軍, 郭詩奇, 楊劍瑜

    (湖南工程學(xué)院理學(xué)院, 湖南湘潭, 411104)

    選用Cu–Ag納米粒子作為研究對象, 采用分子動力學(xué)方法和嵌入原子模型模擬了將Cu原子注入到由Ag原子構(gòu)成的納米粒子的過程, 研究了注入能量(10, 30, 50 eV)與納米粒子結(jié)構(gòu)的關(guān)系。采用共近鄰分析技術(shù)對合金納米的結(jié)構(gòu)演變進(jìn)行了研究。模擬結(jié)果表明: 當(dāng)入射原子動能為50 eV時, 合金納米粒子從截角八面體轉(zhuǎn)變?yōu)槎骟w結(jié)構(gòu)。這種轉(zhuǎn)變是因?yàn)槿肷湓优c基底原子碰撞, 促使系統(tǒng)溫度升高。

    Cu–Ag合金納米粒子; 分子動力學(xué); 近鄰分析技術(shù); 截角八面體; 二十面體

    納米合金是由2種以上納米粒子構(gòu)成的金屬[1–2], 因其獨(dú)特的物理、化學(xué)性質(zhì), 在材料科學(xué)中有著重要的地位及廣泛的應(yīng)用。納米合金的性質(zhì)與其結(jié)構(gòu)(包括構(gòu)型、尺寸、金屬組元及分布、成分比例等)有密切關(guān)系。根據(jù)應(yīng)用需要調(diào)控其結(jié)構(gòu)是納米合金制備技術(shù)發(fā)展的方向。實(shí)驗(yàn)制備納米合金, 一般是讓其在一定基底上生長成形, 其結(jié)構(gòu)隨生長條件(溫度、壓強(qiáng)、組元成分比例、基底構(gòu)型與尺寸等)的變化而演變[3]。如何實(shí)現(xiàn)生長條件可控、協(xié)調(diào)變化, 進(jìn)而對納米合金結(jié)構(gòu)進(jìn)行調(diào)控, 是有重要理論價值與實(shí)踐意義的課題。合金納米粒子的結(jié)構(gòu)源于形成合金納米粒子的組元表面能、原子半徑的相對大小、合金形成能力的強(qiáng)弱、環(huán)境溫度、粒子尺寸與合金成分比例等多種因素。目前國內(nèi)外的研究幾乎都是從靜態(tài)角度, 利用能量優(yōu)化的方法研究合金納米粒子的結(jié)構(gòu)。Ferrando課題組應(yīng)用這一方法研究了同為面心立方結(jié)構(gòu)的2種金屬組成的合金納米粒子結(jié)構(gòu)[4–6], 發(fā)現(xiàn)由于表面聚集, 這些合金納米粒子都有形成核殼結(jié)構(gòu)的趨勢, 但其具體構(gòu)型與粒子尺寸、成分比例有密切關(guān)系; Aguado課題組研究了2種體心立方結(jié)構(gòu)的金屬元素組成的合金納米粒子的結(jié)構(gòu), 發(fā)現(xiàn)尺寸較小時, 合金納米粒子不是核殼結(jié)構(gòu), 而是與純元素納米粒子相似的多面體結(jié)構(gòu)[7]。如上所述, 實(shí)驗(yàn)與理論都表明, 合金納米粒子制備過程中, 其結(jié)構(gòu)與吸附原子在納米粒子表面的擴(kuò)散與生長行為有密切聯(lián)系。在原子尺度上對其結(jié)構(gòu)演變、生長過程進(jìn)行研究有助于理解其性質(zhì), 并達(dá)到人為控制其生長形態(tài)的目的。這些研究結(jié)果說明, 合金納米粒子結(jié)構(gòu)與合金成分比例、組元與溫度等有密切關(guān)系。雖然能量優(yōu)化的方法被廣泛應(yīng)用, 但這種方法有2個缺陷: 第一, 能量局部最小的數(shù)量隨納米粒子包含原子數(shù)的增加呈指數(shù)遞增, 因此要找出全局最小需要極大的計(jì)算量, 這就使得此方法只適用于小尺度系統(tǒng)[6]; 第二, 合金納米粒子制備是一個動力學(xué)過程, 制備條件對納米粒子結(jié)構(gòu)影響非常大, 能量最小對應(yīng)的合金納米粒子結(jié)構(gòu)僅是理想組態(tài), 因此用能量最優(yōu)化方法得到的合金納米粒子結(jié)構(gòu)與實(shí)驗(yàn)結(jié)果有較大差別。由于實(shí)驗(yàn)只是調(diào)節(jié)某一個生長條件, 也沒有考慮到制備過程中生長條件是動態(tài)變化的, 這些實(shí)驗(yàn)也就不能實(shí)現(xiàn)生長條件的可控、協(xié)調(diào)變化, 而且生長條件影響納米合金結(jié)構(gòu)的物理機(jī)制也還有待深入研究。因此, 已有實(shí)驗(yàn)還只是一些零星的、非系統(tǒng)性的研究結(jié)果[3, 8–9]。理論研究方面, 分子動力學(xué)能跟蹤原子運(yùn)動軌跡[10–12], 已經(jīng)被廣泛應(yīng)用于研究合金納米粒子的生長形態(tài)。為研究Cu–Ag合金納米粒子生長形態(tài)的調(diào)控和合金納米粒子結(jié)構(gòu)隨溫度、尺寸、合金成份比例與組元結(jié)構(gòu)等生長條件的變化規(guī)律, 將Cu沉積原子注入到Ag原子納米尺度的基底中, 基底的動能設(shè)置為0, 用高能粒子(Cu原子)的碰撞使基底轉(zhuǎn)變生長方式。為實(shí)現(xiàn)這一目標(biāo), 本文從理論上闡明影響合金納米粒子結(jié)構(gòu)的因素及演變規(guī)律,著重探討注入能量這一因素對截角八面體結(jié)構(gòu)的201個Ag原子生長形態(tài)的影響。

    1 模型與方法

    利用分子動力學(xué)進(jìn)行計(jì)算機(jī)模擬時, 準(zhǔn)確描述原子間作用力勢函數(shù)是模擬成功與否的關(guān)鍵[12–14]。嵌入原子方法(EAM)能夠很好地描述金屬原子間多體相互作用。表面能和晶格常數(shù)2個物理量是納米粒子生長模擬的重要因素, 其中晶格常數(shù)的實(shí)驗(yàn)值被作為輸入?yún)?shù)來確定EAM模型參數(shù)。因此, 應(yīng)用EAM模型計(jì)算了2個元素的表面能()可以驗(yàn)證模型的正確性。在一般情況下, 表面能隨著表面原子的配位數(shù)減少而相應(yīng)減少, 因此對于面心立方結(jié)構(gòu)的金屬, 其表面能的順序應(yīng)該是110>100>111(表1)。EAM模型準(zhǔn)確地預(yù)測了表面能的這一趨勢, 計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果和第一原理結(jié)果也比較接近。因此, 本模型能準(zhǔn)確描述Cu–Ag納米粒子。

    表1 利用第一原理(FP)與嵌入原子方法(EAM)計(jì)算的Cu和Ag表面能及其實(shí)驗(yàn)值(Exp) /( mJ·m-2)

    研究發(fā)現(xiàn)[1–2], 對于面心立方結(jié)構(gòu)(FCC)金屬, 截角八面體(TO)與二十面體(ICO)結(jié)構(gòu)的納米粒子是最穩(wěn)定的結(jié)構(gòu)。因此, 選擇201個原子的截角八面體(TO201)和147個原子的二十面體(ICO147)作為研究對象, 這2種納米粒子的構(gòu)型如圖1所示。截角八面體由8個{111}晶面和6個{100}晶面構(gòu)成, 需要指出的是{111}與{100}面是面心立方結(jié)構(gòu)的密排面, 也是最穩(wěn)定的表面, 截角八面體構(gòu)型與球形比較接近, 這也使得它的結(jié)構(gòu)更為穩(wěn)定。二十面體則是以2個五角形金字塔為基礎(chǔ), 它的表面只有排列最為緊密的{111}晶面, 從這點(diǎn)來說, 它也是一種穩(wěn)定結(jié)構(gòu)。另外, 為便于說明這2種納米粒子的結(jié)構(gòu)特點(diǎn), 分別在TO201與ICO147納米粒子表面上描繪出一個正六邊形與正五邊形, 如圖1中的白色多邊形所示。

    圖1 201個原子的截角八面體和147個原子的二十面體結(jié)構(gòu)

    在生長模擬中, 利用分子動力學(xué)方法模擬將具有一定動能的Cu原子沉積到截角八面體的Ag納米粒子上, 沉積原子的能量范圍是10~50 eV。本文以注入能量為10, 30, 50 eV為例來說明模擬結(jié)構(gòu)。生長模擬在等容等能(NVE)與等容等溫(NVT)的正則系統(tǒng)中進(jìn)行。首先在NVE系統(tǒng)中, 模擬具有一定動能的Cu原子的注入沉積過程, 經(jīng)過0.3 ns的模擬, 系統(tǒng)基本達(dá)到平衡狀態(tài)。然后經(jīng)過2.7 ns的NVT系統(tǒng)模擬, 將系統(tǒng)平衡在溫度為100 K的狀態(tài)。最后記錄所有原子的位置和速度, 以便于做結(jié)構(gòu)和動力學(xué)特性分析。在分子動力學(xué)模擬中, 利用預(yù)測–校正算法求解原子的牛頓運(yùn)動方程, 時間步長取2.0 fs。為保證系統(tǒng)的動量為0, 3個Cu原子被均等地沉積到Ag基底上, 其中每個Cu原子的運(yùn)動方向都指向Ag納米粒子中心。

    如圖2所示, 當(dāng)入射的Cu原子動能為50 eV時, 在開始階段(NVE系統(tǒng)), 系統(tǒng)溫度(約800 K)遠(yuǎn)高于設(shè)定的溫度(100 K), 這是因?yàn)锳g基底原子獲得了入射Cu原子的動能。然后在NVT系統(tǒng)下, Cu–Ag系統(tǒng)溫度逐漸降低, 最后平衡在設(shè)定的100 K左右。所有物理量都在100 K下統(tǒng)計(jì)獲得。

    圖2 Cu原子入射能量為50 eV時基底納米粒子溫度隨時間變化

    2 結(jié)果與討論

    當(dāng)把Cu原子注入到Ag納米粒子TO201, 納米合金首先在NVE中模擬沉積過程(見圖2), 入射能量等于50 eV時, Cu原子的入射能量傳遞到基底, 基底原子之間劇烈碰撞, 使系統(tǒng)溫度變?yōu)榧s850 K。NVE中模擬的850 K高溫會導(dǎo)致Ag納米粒子基底開始快速扭曲變形, 然后經(jīng)過2.7 ns的NVT模擬, 系統(tǒng)平衡溫度冷卻至100 K, 納米粒子由高溫(850 K)時的液態(tài)凝固成一定結(jié)構(gòu)的固態(tài)。通常的分析技術(shù), 如徑向分布函數(shù)和靜態(tài)結(jié)構(gòu)因子, 只能從整體上對結(jié)構(gòu)作有序與無序的判斷。為更細(xì)致地研究納米合金的結(jié)構(gòu)變化過程, 采用共近鄰分析技術(shù)(CNA)確定凝固的Cu–Ag合金納米粒子結(jié)構(gòu)。共近鄰分析技術(shù)通過3位數(shù)(,,)描述原子所屬的狀態(tài),表示成鍵的2原子的共有最近鄰原子數(shù),表示共有最近鄰原子之間的成鍵數(shù),表示共有最近鄰原子所成鍵鏈中最長鏈上的鍵數(shù)。當(dāng)將Cu原子入射到Ag基底中, 入射能量為10, 30, 50 eV時, 共近鄰分析的結(jié)果列于表2。

    表2 Cu注入到Ag的TO201中(3, 1, 1), (4, 2, 1),(4, 2, 2)和(5, 5, 5)等鍵對占總鍵對的百分比 %

    根據(jù)共近鄰分析方法, 每一種鍵對都對應(yīng)著特定的原子空間幾何構(gòu)型。其中FCC型晶體中以(4, 2, 1)鍵為特征鍵對; HCP型晶體則同時以(4, 2, 1)和(4, 2, 2)為特征鍵對。(3, 1, 1)是納米粒子的表面原子, (5, 5, 5)是二十面體5重對稱性的晶序特性。因此僅需要(3, 1, 1), (4, 2, 1), (4, 2, 2), (5, 5, 5)這幾種晶序就能區(qū)分TO201, ICO147結(jié)構(gòu)。

    如表2所示, 當(dāng)Cu原子以10 eV和30 eV的能量入射TO201時, 新的Cu–Ag合金納米粒子中沒有(4, 2, 2)和(5, 5, 5)構(gòu)型, 僅有(3, 1, 1)與(4, 2, 1)鍵, 且這2種鍵型所占百分比與Ag的TO201類似, TO201的Ag基底具有{111}和{100}面, 經(jīng)CNA分析, (4, 2, 2)和(5, 5, 5)鍵對數(shù)為0。這表明當(dāng)Cu原子以較低能量入射時, Ag基底仍保持著原來的結(jié)構(gòu)而沒有發(fā)生變化。當(dāng)入射的Cu原子動能增加大到50 eV時, (4, 2, 2)鍵所占百分比提高到了18%, (4, 2, 1)鍵的比例減少到41%。需要指出的是, 在此情形下, 有少量的(5, 5, 5)鍵型。這種鍵型特征與ICO147類似。共近鄰分析方法說明, 由于高能量的Cu原子注入, Ag納米粒子由TO轉(zhuǎn)變?yōu)镮CO結(jié)構(gòu)。用Cu原子注入后, 合金納米粒子結(jié)構(gòu)演變過程圖更直接地說明這種轉(zhuǎn)變。

    如圖3所示, 動能為10 eV的Cu原子入射到Ag的TO基底上, Ag基底幾乎保持了原來的TO結(jié)構(gòu), 入射的Cu只是嵌入到Ag基底中。最終平衡后, 合金納米粒子保持了原來的六邊對稱性。

    當(dāng)入射的Cu原子動能為50 eV時(圖4), 與10 eV相比, 50 eV的Cu原子與基底的Ag原子強(qiáng)烈作用, 粒子的運(yùn)動過程要更加復(fù)雜, 溫度急劇升高促使結(jié)構(gòu)發(fā)生變化, 由6重對稱演變成5重對稱。Ag納米粒子出現(xiàn)了5重對稱結(jié)構(gòu), 并出現(xiàn)了幾個{100}面, 這樣的結(jié)構(gòu)更類似于ICO構(gòu)型。另外, 因Ag原子表面能較低、原子半徑較大, 最終導(dǎo)致50 eV入射的Cu沉積原子在進(jìn)入TO后偏析于Ag納米粒子表面。

    如圖5所示, 對尺寸較小的Cu納米粒子, ICO比TO穩(wěn)定, 因此TO比ICO更容易演變成其它結(jié)構(gòu)。對所有尺寸的Ag納米粒子, TO比ICO更穩(wěn)定, 這一結(jié)果與其它理論的計(jì)算結(jié)果一致。當(dāng)Cu注入到Ag的TO201基底中, 因TO很穩(wěn)定, 其很難演變成ICO。因此, 只有當(dāng)以較大動能的Cu原子注入到Ag基底中, 使基底原子發(fā)生較大變化, 才能使Ag基底從TO轉(zhuǎn)化為ICO。

    圖3 注入能量為10 eV的Cu原子沉積到結(jié)構(gòu)為TO的Ag基底中納米粒子構(gòu)型演變

    圖4 注入能量為50 eV的Cu原子沉積到結(jié)構(gòu)為TO的Ag基底中納米粒子構(gòu)型演變

    圖5 (Cu, Ag)TO和ICO結(jié)構(gòu)結(jié)合能隨原子尺寸的變化

    為了解釋這一模擬結(jié)果, 進(jìn)行了一系列不同尺寸的TO和ICO納米粒子的MD模擬實(shí)驗(yàn), 證明每個原子結(jié)合能的大小與原子尺寸相關(guān)。

    3 結(jié)論

    用分子動力學(xué)和原子嵌入模型研究了能量在10~50 eV的Cu原子沉積到TO201構(gòu)型的Ag納米粒子基底的過程, 分析了生長出的Cu–Ag納米合金結(jié)構(gòu)的早期演變過程。應(yīng)用共近鄰分析技術(shù), 發(fā)現(xiàn)當(dāng)入射Cu原子的能量為10 eV和30 eV時, 生長出的納米合金粒子基本保持了原有TO形態(tài); 在50 eV時, 納米粒子由六重對稱性變?yōu)槲逯貙ΨQ性。因沉積原子的能量被均勻地分配到基底每個原子, 基底原子不斷擠壓碰撞導(dǎo)致了系統(tǒng)溫度的升高, 溫度的升高進(jìn)而導(dǎo)致了納米粒子結(jié)構(gòu)的轉(zhuǎn)變。

    參考文獻(xiàn):

    [1] Baletto F, Ferrando R. Structural properties of nanoclusters: Energetic, thermodynamic, and kinetic effects [J]. Rev Mod Phys, 2005, 77(1): 371–423.

    [2] Ferrando R, Jellinek J, Johnston R L. Nanoalloys: from theory to applications of alloy clusters and nanoparticles [J]. Chem Rev, 2008, 108(3): 845–910.

    [3] Chaudhuri R G, Paria S. Core/shell nanoparticles: classes, properties, synthesis mechanisms, characterization, and applications [J]. Chem Rev, 2012, 112(4): 2 373–2 433.

    [4] Rapallo A, Rossi G, Ferrando R, et al. Global optimization of bimetallic cluster structures. I. Size-mismatched Ag-Cu, Ag-Ni, and Au-Cu systems [J]. J Chem Phys, 2005, 122(19): 194308(1–13).

    [5] Rapallo A, Rossi G, Ferrando R, et al. Global optimization of bimetallic cluster structures. II. Size-matched Ag-Pd, Ag-Au, and Pd-Pt systems [J]. J Chem Phys, 2005, 122(19): 194309(1–9).

    [6] Negreiros F R, Barcaro G, Kuntova Z. et al. Structures of gas-phase Ag-Pd nanoclusters: A computational study [J]. J Chem Phys, 2010, 132(23): 234703(1–8).

    [7] Aguado A, Lopez J M. Structure determination in 55-atom Li-Na and Na-K nanoalloys [J]. J Chem Phys, 2010, 133(9): 094302(1–10).

    [8] Serpell C J, Cookson J, Ozkaya D, et al. Core@shell bimetal-lic nanoparticle synthesis via anion coordination [J]. Nat Chem, 2011, 3(6): 478–483.

    [9] Radilovic V, Ophus C, Marquis E A, et al. Highly monodisperse core-shell particles created by solid-state reactions [J]. Nat Mater, 2011, 10(9): 710–715.

    [10] 文玉華, 朱如曾, 周富信. 分子動力學(xué)模擬的主要技術(shù)[J]. 力學(xué)進(jìn)展, 2003, 33(1): 67–71.

    [11] 蔡錫年. 分子動力學(xué)和物理力學(xué)[M]. 北京: 科學(xué)出版社, 1986: 100–110.

    [12] 周健, 陸小華, 王延儒, 等. 分子動力學(xué)模擬[J].南京化工大學(xué)學(xué)報(bào), 1997, 19(2): 20–25.

    [13] Vitos L, Ruban A V, Skriver H L. et al. The surface energy of metals [J]. Surf Sci, 1998, 411(1–2): 186–202.

    [14] Tyson W R, Miller W A. A Comprehensive Search for Stable Pt-Pd Nanoalloy Configurations and Their Use as Tunable Catalysts [J]. Surf Sci, 1977, 62(1): 267–276.

    (責(zé)任編校:劉剛毅)

    The growth configuration of nanoalloy: a molecular dynamics simulation

    Yu Yi, Cheng Qingjun, Guo Shiqi, Yang Jianyu

    (Science College, Hunan Institute of Engineering, Xiangtan 411104, China)

    The impact depositions of Cu atoms onto the Ag truncated octahedron are simulated by molecular dynamics along with the embedded atom method. The dependence of incident energies (10, 30 and 50 eV) and the nanoalloy configuration is found. The common neighbor analysis is used to identify the structural evolutions. A structural transformation from truncated octahedral to icosahedral arrangement is observed at incident energy of 50 eV. The structural transformation can be induced by the increase of system temperature, because the incident energy of the deposited atom is delivered to the substrate atoms as a result of collisions.

    Cu-Ag nanoalloy; molecular dynamics; common neighbor analysis; truncated octahedron; icosahedron

    10.3969/j.issn.1672–6146.2015.01.005

    O 781

    1672–6146(2015)01–0017–04

    楊劍瑜, hnieyjy@aliyun.com; 余藝, 1847997096@qq.com。

    2014–11–15

    國家自然科學(xué)基金(50971058), 2014年湖南省大學(xué)生科技創(chuàng)新項(xiàng)目。

    猜你喜歡
    基底原子合金
    《我要我們在一起》主打現(xiàn)實(shí)基底 務(wù)必更接地氣
    中國銀幕(2022年4期)2022-04-07 21:28:24
    原子可以結(jié)合嗎?
    原子究竟有多?。?/a>
    帶你認(rèn)識原子
    鉬鎢合金燒結(jié)致密化行為
    Zr-Nb-Fe合金在LiOH溶液中的腐蝕行為
    可溶巖隧道基底巖溶水處理方案探討
    Hastelloy C-2000合金的焊接工藝
    焊接(2016年8期)2016-02-27 13:05:16
    磁共振顯像對老年椎基底動脈缺血的診斷價值
    為什么錫鉍合金被稱為最突出的焊接材料?
    久久天躁狠狠躁夜夜2o2o| 精品亚洲成国产av| 少妇 在线观看| 777久久人妻少妇嫩草av网站| 精品高清国产在线一区| 丁香六月欧美| 女人爽到高潮嗷嗷叫在线视频| 在线观看免费高清a一片| 久久 成人 亚洲| 高清欧美精品videossex| 变态另类成人亚洲欧美熟女 | 国产欧美日韩精品亚洲av| 在线永久观看黄色视频| 涩涩av久久男人的天堂| 天天躁狠狠躁夜夜躁狠狠躁| 男女午夜视频在线观看| 99re在线观看精品视频| 啦啦啦免费观看视频1| 精品人妻在线不人妻| 十八禁网站免费在线| 亚洲九九香蕉| 高清毛片免费观看视频网站 | 两性午夜刺激爽爽歪歪视频在线观看 | 午夜91福利影院| 久久久久网色| 12—13女人毛片做爰片一| 亚洲欧美色中文字幕在线| 国产精品国产av在线观看| 亚洲精品粉嫩美女一区| 性色av乱码一区二区三区2| 高清毛片免费观看视频网站 | 国产深夜福利视频在线观看| 久久精品aⅴ一区二区三区四区| 久热这里只有精品99| 国产精品二区激情视频| 国产免费现黄频在线看| 亚洲国产欧美一区二区综合| 欧美人与性动交α欧美精品济南到| 别揉我奶头~嗯~啊~动态视频| 菩萨蛮人人尽说江南好唐韦庄| 国产精品亚洲av一区麻豆| 成人18禁高潮啪啪吃奶动态图| 91精品国产国语对白视频| 成年人午夜在线观看视频| 成人av一区二区三区在线看| 精品国产一区二区三区四区第35| 亚洲,欧美精品.| 久久中文字幕人妻熟女| 欧美亚洲 丝袜 人妻 在线| aaaaa片日本免费| 亚洲精品国产色婷婷电影| 亚洲av成人一区二区三| 女同久久另类99精品国产91| 国产深夜福利视频在线观看| 另类亚洲欧美激情| 免费高清在线观看日韩| 99riav亚洲国产免费| 十分钟在线观看高清视频www| √禁漫天堂资源中文www| 91精品国产国语对白视频| 老司机深夜福利视频在线观看| 在线观看免费视频日本深夜| 涩涩av久久男人的天堂| 人人妻人人爽人人添夜夜欢视频| 十八禁网站网址无遮挡| 国产亚洲精品久久久久5区| 视频在线观看一区二区三区| 亚洲国产av影院在线观看| 久久久久精品国产欧美久久久| 国产1区2区3区精品| 在线观看66精品国产| 丁香六月天网| 精品少妇黑人巨大在线播放| 亚洲国产精品一区二区三区在线| 久久九九热精品免费| 2018国产大陆天天弄谢| 在线观看免费高清a一片| 麻豆乱淫一区二区| 国产真人三级小视频在线观看| 高清视频免费观看一区二区| 亚洲七黄色美女视频| 欧美日韩一级在线毛片| 黄色丝袜av网址大全| 久久精品亚洲精品国产色婷小说| 成人特级黄色片久久久久久久 | 黄频高清免费视频| 一区在线观看完整版| 亚洲精品乱久久久久久| 久久99热这里只频精品6学生| 欧美性长视频在线观看| 中文字幕制服av| 久久毛片免费看一区二区三区| 欧美激情极品国产一区二区三区| 久久人妻熟女aⅴ| 精品国产乱码久久久久久小说| 中文字幕色久视频| 亚洲精品av麻豆狂野| 9热在线视频观看99| 嫩草影视91久久| 国产精品免费一区二区三区在线 | 日韩欧美三级三区| 亚洲欧美一区二区三区黑人| 久久久精品国产亚洲av高清涩受| svipshipincom国产片| 午夜91福利影院| 欧美 日韩 精品 国产| 一区福利在线观看| 国产男女超爽视频在线观看| 香蕉国产在线看| 国产成人免费观看mmmm| 啦啦啦中文免费视频观看日本| 欧美国产精品va在线观看不卡| 侵犯人妻中文字幕一二三四区| 99久久人妻综合| 人人妻人人爽人人添夜夜欢视频| 女同久久另类99精品国产91| 嫩草影视91久久| 高清视频免费观看一区二区| 一级片免费观看大全| 可以免费在线观看a视频的电影网站| 久久青草综合色| 一级黄色大片毛片| 精品一区二区三卡| 大陆偷拍与自拍| 69精品国产乱码久久久| www.精华液| 两性午夜刺激爽爽歪歪视频在线观看 | 久久天堂一区二区三区四区| 18禁美女被吸乳视频| 国产麻豆69| 99精品在免费线老司机午夜| 国产亚洲一区二区精品| 精品人妻在线不人妻| 久久精品国产亚洲av高清一级| 国产欧美日韩综合在线一区二区| 国产精品欧美亚洲77777| 欧美 亚洲 国产 日韩一| 99久久人妻综合| 国产一区二区激情短视频| 国产黄频视频在线观看| 国产精品 欧美亚洲| 不卡av一区二区三区| av天堂在线播放| 成年人免费黄色播放视频| 桃红色精品国产亚洲av| 多毛熟女@视频| 曰老女人黄片| 精品国产一区二区三区四区第35| 国产一区二区三区在线臀色熟女 | 久久国产精品影院| 黄网站色视频无遮挡免费观看| 99久久人妻综合| 国产一区二区三区综合在线观看| 亚洲精品美女久久久久99蜜臀| 欧美日韩亚洲综合一区二区三区_| 亚洲精品美女久久av网站| 亚洲欧美色中文字幕在线| tube8黄色片| 人人妻人人爽人人添夜夜欢视频| 亚洲,欧美精品.| 人成视频在线观看免费观看| 欧美黄色淫秽网站| 亚洲综合色网址| 咕卡用的链子| 人人妻,人人澡人人爽秒播| 女人久久www免费人成看片| 9色porny在线观看| 老司机福利观看| 脱女人内裤的视频| 欧美日韩精品网址| 精品一品国产午夜福利视频| 欧美精品一区二区大全| 欧美中文综合在线视频| 国产精品久久久久成人av| 国产野战对白在线观看| 亚洲av电影在线进入| 欧美人与性动交α欧美精品济南到| 纵有疾风起免费观看全集完整版| 日本av免费视频播放| 菩萨蛮人人尽说江南好唐韦庄| 国产成人av教育| 亚洲伊人久久精品综合| 如日韩欧美国产精品一区二区三区| 超色免费av| 黑人猛操日本美女一级片| 国产亚洲午夜精品一区二区久久| 亚洲精品乱久久久久久| 一本综合久久免费| 亚洲精品一卡2卡三卡4卡5卡| 不卡一级毛片| 国产精品久久久av美女十八| www.999成人在线观看| 人人妻人人澡人人爽人人夜夜| 国产欧美日韩一区二区精品| 国产精品秋霞免费鲁丝片| 久久久久久亚洲精品国产蜜桃av| 黄片播放在线免费| 飞空精品影院首页| 女性生殖器流出的白浆| 久久人妻熟女aⅴ| 色精品久久人妻99蜜桃| 久9热在线精品视频| 50天的宝宝边吃奶边哭怎么回事| 建设人人有责人人尽责人人享有的| 午夜福利乱码中文字幕| 亚洲国产欧美一区二区综合| 亚洲精品国产色婷婷电影| 午夜两性在线视频| 一本一本久久a久久精品综合妖精| 最新的欧美精品一区二区| 最近最新中文字幕大全电影3 | 最新的欧美精品一区二区| 欧美久久黑人一区二区| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品国产精品久久久不卡| 久9热在线精品视频| 少妇猛男粗大的猛烈进出视频| 飞空精品影院首页| 色94色欧美一区二区| 一进一出抽搐动态| 精品第一国产精品| 757午夜福利合集在线观看| 国产又爽黄色视频| 一级片'在线观看视频| 亚洲av成人一区二区三| 午夜精品国产一区二区电影| 一二三四社区在线视频社区8| 黄色视频,在线免费观看| 69精品国产乱码久久久| 2018国产大陆天天弄谢| 午夜日韩欧美国产| 动漫黄色视频在线观看| 精品国产一区二区三区四区第35| 1024香蕉在线观看| 99国产综合亚洲精品| 中国美女看黄片| 18禁裸乳无遮挡动漫免费视频| 午夜免费鲁丝| 久久精品亚洲精品国产色婷小说| 亚洲精品中文字幕在线视频| 欧美av亚洲av综合av国产av| 成年人免费黄色播放视频| 国产成人精品久久二区二区91| 亚洲国产欧美在线一区| 757午夜福利合集在线观看| www.自偷自拍.com| 亚洲精品久久成人aⅴ小说| 性高湖久久久久久久久免费观看| 国产熟女午夜一区二区三区| 亚洲中文日韩欧美视频| 又紧又爽又黄一区二区| 成人国产一区最新在线观看| 91精品三级在线观看| 高清毛片免费观看视频网站 | 男人舔女人的私密视频| 18禁国产床啪视频网站| 国产福利在线免费观看视频| 99久久国产精品久久久| 无限看片的www在线观看| tube8黄色片| 国产免费现黄频在线看| 1024视频免费在线观看| 搡老乐熟女国产| 欧美性长视频在线观看| 欧美日韩亚洲高清精品| 成人亚洲精品一区在线观看| 热99re8久久精品国产| 国精品久久久久久国模美| 又紧又爽又黄一区二区| 亚洲成人国产一区在线观看| 人成视频在线观看免费观看| 中文欧美无线码| 国产精品久久电影中文字幕 | 免费在线观看黄色视频的| 久久99热这里只频精品6学生| 国产精品免费视频内射| 欧美黑人欧美精品刺激| 日本a在线网址| 十分钟在线观看高清视频www| 啦啦啦视频在线资源免费观看| 精品国产超薄肉色丝袜足j| 亚洲第一青青草原| 久久人妻熟女aⅴ| 黄色丝袜av网址大全| 成年人午夜在线观看视频| 国产日韩欧美亚洲二区| 美女高潮到喷水免费观看| 成年女人毛片免费观看观看9 | 99久久99久久久精品蜜桃| 久久国产精品男人的天堂亚洲| 91精品三级在线观看| av福利片在线| 国产成人av激情在线播放| 一级黄色大片毛片| 人人妻人人爽人人添夜夜欢视频| 黄片大片在线免费观看| 夫妻午夜视频| 亚洲精品国产区一区二| 国产精品.久久久| 黄片大片在线免费观看| 精品久久蜜臀av无| 一本综合久久免费| 欧美在线黄色| 国产男女超爽视频在线观看| 久久久国产成人免费| 在线观看免费高清a一片| 999精品在线视频| 精品国产乱码久久久久久男人| 亚洲国产欧美在线一区| 另类精品久久| 久久久水蜜桃国产精品网| 中文字幕人妻熟女乱码| 国产精品熟女久久久久浪| 亚洲男人天堂网一区| 亚洲人成伊人成综合网2020| 免费在线观看日本一区| 国产一区有黄有色的免费视频| 夜夜爽天天搞| 怎么达到女性高潮| 国产亚洲av高清不卡| 人妻一区二区av| 欧美精品一区二区免费开放| 天天躁日日躁夜夜躁夜夜| 久久亚洲精品不卡| 少妇裸体淫交视频免费看高清 | 午夜91福利影院| 一进一出好大好爽视频| 国产成人一区二区三区免费视频网站| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲五月色婷婷综合| 国产在线免费精品| 国产精品免费大片| 老司机深夜福利视频在线观看| 人妻 亚洲 视频| 一进一出抽搐动态| 一二三四在线观看免费中文在| 老汉色av国产亚洲站长工具| 午夜激情久久久久久久| 不卡一级毛片| cao死你这个sao货| 色综合婷婷激情| 欧美老熟妇乱子伦牲交| 国产99久久九九免费精品| 婷婷成人精品国产| 国产成人精品久久二区二区91| 操美女的视频在线观看| 麻豆成人av在线观看| 国产精品av久久久久免费| 国产在线精品亚洲第一网站| 丝袜喷水一区| 91国产中文字幕| 欧美国产精品va在线观看不卡| 日韩欧美三级三区| 久久久久国产一级毛片高清牌| 日韩中文字幕欧美一区二区| 久久中文字幕一级| 9色porny在线观看| 精品国产一区二区三区四区第35| 欧美成狂野欧美在线观看| 老司机深夜福利视频在线观看| 欧美激情 高清一区二区三区| 精品久久久久久电影网| 女同久久另类99精品国产91| 国内毛片毛片毛片毛片毛片| 黄色怎么调成土黄色| 欧美黑人精品巨大| 免费不卡黄色视频| 久久久精品区二区三区| 18禁美女被吸乳视频| 国产1区2区3区精品| 一进一出好大好爽视频| 国产单亲对白刺激| 久久久久久久国产电影| 国产黄色免费在线视频| 欧美 亚洲 国产 日韩一| 欧美激情极品国产一区二区三区| 欧美精品一区二区免费开放| 老司机深夜福利视频在线观看| 高清视频免费观看一区二区| 国产亚洲欧美精品永久| 亚洲 国产 在线| 丝袜人妻中文字幕| 99国产精品一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 欧美日韩视频精品一区| 9191精品国产免费久久| 精品一区二区三卡| 一区在线观看完整版| 色婷婷久久久亚洲欧美| 久久久国产精品麻豆| 黑人欧美特级aaaaaa片| 国产不卡一卡二| 国产av又大| av国产精品久久久久影院| 欧美变态另类bdsm刘玥| 亚洲少妇的诱惑av| 欧美精品人与动牲交sv欧美| 搡老熟女国产l中国老女人| 国产欧美亚洲国产| 久久热在线av| 一区福利在线观看| 首页视频小说图片口味搜索| 日韩免费av在线播放| 国产精品成人在线| 亚洲精品一卡2卡三卡4卡5卡| 国产欧美日韩精品亚洲av| 久久久久精品国产欧美久久久| 国产男靠女视频免费网站| 最黄视频免费看| 午夜两性在线视频| 久久亚洲真实| 无限看片的www在线观看| 久久久国产精品麻豆| 2018国产大陆天天弄谢| 成人国产一区最新在线观看| 制服诱惑二区| 亚洲熟女毛片儿| 纵有疾风起免费观看全集完整版| 亚洲中文日韩欧美视频| 99在线人妻在线中文字幕 | 精品少妇久久久久久888优播| 午夜福利视频在线观看免费| 国产黄频视频在线观看| 精品福利观看| 久久天躁狠狠躁夜夜2o2o| 欧美激情久久久久久爽电影 | 一级片免费观看大全| 中文字幕高清在线视频| 日韩欧美一区视频在线观看| 精品国产亚洲在线| 久久中文字幕人妻熟女| 午夜福利视频精品| 91精品三级在线观看| 国产精品 欧美亚洲| 韩国精品一区二区三区| 香蕉国产在线看| av免费在线观看网站| 国产av精品麻豆| 丰满迷人的少妇在线观看| avwww免费| 日韩中文字幕视频在线看片| 嫩草影视91久久| av超薄肉色丝袜交足视频| 午夜日韩欧美国产| 十八禁网站网址无遮挡| 精品第一国产精品| 国产欧美日韩一区二区三| 欧美另类亚洲清纯唯美| 亚洲第一av免费看| 一本—道久久a久久精品蜜桃钙片| 国产精品久久久久久精品古装| 欧美精品人与动牲交sv欧美| 免费日韩欧美在线观看| 久久中文字幕人妻熟女| 一区二区三区国产精品乱码| 亚洲视频免费观看视频| 亚洲国产欧美在线一区| 在线观看www视频免费| 亚洲国产欧美一区二区综合| 激情在线观看视频在线高清 | 日韩人妻精品一区2区三区| 久久精品亚洲熟妇少妇任你| 国内毛片毛片毛片毛片毛片| 一区二区三区国产精品乱码| 精品少妇黑人巨大在线播放| 国产激情久久老熟女| 久久久久久久国产电影| 中文字幕人妻熟女乱码| 12—13女人毛片做爰片一| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲中文字幕日韩| 在线观看免费日韩欧美大片| 亚洲成人免费电影在线观看| 桃红色精品国产亚洲av| 老司机在亚洲福利影院| videosex国产| 亚洲色图综合在线观看| 一二三四社区在线视频社区8| 久久午夜综合久久蜜桃| 国产不卡av网站在线观看| 搡老乐熟女国产| 宅男免费午夜| 91大片在线观看| 欧美激情 高清一区二区三区| 法律面前人人平等表现在哪些方面| 色在线成人网| 在线av久久热| 国产精品久久久久久精品古装| 亚洲av日韩精品久久久久久密| 最新美女视频免费是黄的| 国产伦理片在线播放av一区| 国产一区二区 视频在线| 人人妻人人爽人人添夜夜欢视频| 国产片内射在线| 亚洲成人国产一区在线观看| 制服人妻中文乱码| a级片在线免费高清观看视频| 狂野欧美激情性xxxx| 免费少妇av软件| 久久久久国产一级毛片高清牌| tube8黄色片| 可以免费在线观看a视频的电影网站| 精品亚洲乱码少妇综合久久| aaaaa片日本免费| 在线观看人妻少妇| av又黄又爽大尺度在线免费看| 国产男靠女视频免费网站| 亚洲免费av在线视频| 精品亚洲成a人片在线观看| 十分钟在线观看高清视频www| 在线亚洲精品国产二区图片欧美| 亚洲精品自拍成人| 国产精品免费视频内射| 午夜精品久久久久久毛片777| 午夜免费成人在线视频| 国产人伦9x9x在线观看| 天天添夜夜摸| 中文亚洲av片在线观看爽 | 看免费av毛片| 狠狠狠狠99中文字幕| 国产无遮挡羞羞视频在线观看| 国产麻豆69| 国产av一区二区精品久久| 午夜福利视频精品| 日韩人妻精品一区2区三区| 亚洲国产欧美日韩在线播放| 国产国语露脸激情在线看| 亚洲伊人色综图| 免费不卡黄色视频| 国产在线一区二区三区精| 极品人妻少妇av视频| 无遮挡黄片免费观看| 欧美另类亚洲清纯唯美| 国产不卡av网站在线观看| 日本av手机在线免费观看| 欧美日韩中文字幕国产精品一区二区三区 | 色综合欧美亚洲国产小说| 国产主播在线观看一区二区| 国产aⅴ精品一区二区三区波| 国产av一区二区精品久久| 黑人欧美特级aaaaaa片| 亚洲精品国产精品久久久不卡| 中文字幕av电影在线播放| 天天添夜夜摸| 女性被躁到高潮视频| 麻豆av在线久日| 亚洲精品中文字幕在线视频| 无人区码免费观看不卡 | 久久久久久亚洲精品国产蜜桃av| 精品亚洲成a人片在线观看| 成人三级做爰电影| 国产精品国产av在线观看| 天堂8中文在线网| 国产免费现黄频在线看| 欧美人与性动交α欧美软件| 亚洲国产欧美一区二区综合| 日本黄色视频三级网站网址 | 成人国产av品久久久| 亚洲天堂av无毛| av天堂在线播放| 欧美日韩亚洲国产一区二区在线观看 | 一本久久精品| 一级片免费观看大全| 亚洲精品中文字幕在线视频| 菩萨蛮人人尽说江南好唐韦庄| 在线观看人妻少妇| 国产亚洲午夜精品一区二区久久| 一进一出好大好爽视频| 国产av国产精品国产| 俄罗斯特黄特色一大片| 欧美精品亚洲一区二区| 男人操女人黄网站| 亚洲少妇的诱惑av| 久久久久国产一级毛片高清牌| 无遮挡黄片免费观看| 欧美日韩亚洲国产一区二区在线观看 | 久久影院123| 99在线人妻在线中文字幕 | 超碰成人久久| 久久天躁狠狠躁夜夜2o2o| 亚洲精品中文字幕在线视频| 别揉我奶头~嗯~啊~动态视频| 精品国产亚洲在线| 一进一出抽搐动态| 搡老岳熟女国产| 丰满迷人的少妇在线观看| av有码第一页| 在线十欧美十亚洲十日本专区| av超薄肉色丝袜交足视频| 老司机靠b影院| 国产精品亚洲av一区麻豆| 久久毛片免费看一区二区三区| 久久久久久久久免费视频了| 老司机影院毛片| 老司机在亚洲福利影院| 免费观看a级毛片全部| 一个人免费看片子| 久久午夜综合久久蜜桃| 91成年电影在线观看| 曰老女人黄片| 国产熟女午夜一区二区三区| tube8黄色片| 丰满迷人的少妇在线观看| 久久国产精品大桥未久av| 欧美精品亚洲一区二区| 欧美日韩一级在线毛片| 日韩视频在线欧美| 日本av手机在线免费观看| 少妇裸体淫交视频免费看高清 | www.999成人在线观看| 9191精品国产免费久久| 久久亚洲真实| 亚洲专区中文字幕在线| 精品午夜福利视频在线观看一区 | av福利片在线|