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

    Mn(SiO2)3(M=Fe,Co,Ni;n=1-3)團(tuán)簇的幾何結(jié)構(gòu)、光電性質(zhì)及磁性

    2014-06-09 12:33:45智麗麗李艷青楊蓮紅趙高峰
    計算物理 2014年6期
    關(guān)鍵詞:磁矩基態(tài)鍵長

    智麗麗, 李艷青, 楊蓮紅, 趙高峰

    (1.伊犁師范學(xué)院新疆凝聚態(tài)相變與微結(jié)構(gòu)實驗室,伊寧 835000;2.新疆昌吉學(xué)院物理系,昌吉 831100;3.山東大學(xué)晶體材料國家重點實驗室,濟(jì)南 250100;4.河南大學(xué)物理與電子學(xué)院,開封 475004)

    Mn(SiO2)3(M=Fe,Co,Ni;n=1-3)團(tuán)簇的幾何結(jié)構(gòu)、光電性質(zhì)及磁性

    智麗麗1,2, 李艷青2,3, 楊蓮紅2, 趙高峰4

    (1.伊犁師范學(xué)院新疆凝聚態(tài)相變與微結(jié)構(gòu)實驗室,伊寧 835000;2.新疆昌吉學(xué)院物理系,昌吉 831100;3.山東大學(xué)晶體材料國家重點實驗室,濟(jì)南 250100;4.河南大學(xué)物理與電子學(xué)院,開封 475004)

    基于密度泛函理論中的廣義梯度近似系統(tǒng)研究Mn(SiO2)3(M=Fe,Co,Ni;n=1-3)團(tuán)簇的幾何結(jié)構(gòu)、光電性質(zhì)和磁學(xué)性質(zhì).結(jié)果表明:Fe,Co原子相對于Ni原子更易于在(SiO2)3團(tuán)簇上聚集;通過分析團(tuán)簇的分裂途徑及其產(chǎn)物,發(fā)現(xiàn)穩(wěn)定性較好的氧化硅是一種很好的用于負(fù)載過渡金屬“島膜”的載體材料;Mn(SiO2)3團(tuán)簇的能隙恰好位于近紅外光譜范圍內(nèi).通過磁性分析發(fā)現(xiàn),該復(fù)合團(tuán)簇的磁矩主要局域在過渡金屬原子周圍,而且,F(xiàn)e2(SiO2)3和Co3(SiO2)3具有相對較大的磁矩,這主要源于過渡金屬原子的d軌道間相互耦合.能隙和磁性兩方面性質(zhì)進(jìn)一步肯定了二氧化硅磁性復(fù)合材料在醫(yī)學(xué)界被用作光動力靶向治療的可觀前景.

    Mn(SiO2)3團(tuán)簇;幾何結(jié)構(gòu);光電性質(zhì);磁性

    0 引言

    隨著納米技術(shù)的發(fā)展,人們發(fā)現(xiàn)納米量級的物質(zhì)具有與塊體完全不同的性質(zhì),特別是作為團(tuán)簇分子其物理化學(xué)性質(zhì)隨團(tuán)簇的尺寸會發(fā)生明顯的變化.硅的氧化物是地球上含量和用途較廣的物質(zhì)之一,尤其是二氧化硅在微電子、光通訊、和薄膜技術(shù)等許多領(lǐng)域中有著廣泛的用途,因此,二氧化硅團(tuán)簇和納米顆粒方面的研究在近幾年成為熱點[1-3].有關(guān)報道指明(SiO2)n(n≤6)團(tuán)簇在其小尺寸時,其基態(tài)結(jié)構(gòu)是鏈狀的構(gòu)型[4-6].

    目前,金屬-氧化硅聚合物也是科學(xué)界關(guān)注的焦點,原因在于這些聚合物有著獨特的性能,如:孫強(qiáng)等人[7]研究了Au原子與二氧化硅團(tuán)簇的相互作用,他們指出該團(tuán)簇體系可吸收近紅外光,并將光轉(zhuǎn)化為熱量從而來殺死癌細(xì)胞.本課題組成員分別計算了Cu和Ag原子在二氧化硅團(tuán)簇上的吸附,結(jié)果表明貴金屬吸附二氧化硅均可以吸收近紅外光,從而被用于癌細(xì)胞的治療中[6,8].最近研究表明,(Fe\Co\Ni)/SiO2復(fù)合納米微球不僅具有高的化學(xué)穩(wěn)定性,大的比表面積,超順磁性,還具有良好的生物相容性[7].這些性能使其具有極廣闊的應(yīng)用前景,因而在核酸提取,藥物靶向等生物醫(yī)學(xué)領(lǐng)域有著廣泛的應(yīng)用.值得關(guān)注的是:李雅麗[9]等人制備的二氧化硅磁性復(fù)合顆粒類似于貴金屬二氧化硅復(fù)合材料,也可以吸收可見-近紅外光,而且該微粒還具有超順磁性,因此,該復(fù)合顆粒有可能被用于光動力靶向治療中.據(jù)有關(guān)報道指出Ni/SiO2納米復(fù)合材料也有著獨特的光學(xué)和磁學(xué)性能[10].

    鑒于(Fe\Co\Ni)/SiO2復(fù)合材料獨特的光學(xué)和磁學(xué)特性,本文采用基于密度泛函理論中的廣義梯度近似,從分子水平上系統(tǒng)地研究了Mn(SiO2)3(M=Fe,Co,Ni;n=1-3)團(tuán)簇的幾何結(jié)構(gòu)、光電性質(zhì)及磁學(xué)性質(zhì),這將有利于從微觀上理解這一復(fù)合體系所具有的光電磁特性,以及過渡金屬與二氧化硅之間的相互作用.本文以(SiO2)3作為基體,原因在于實驗上已經(jīng)成功合成了(SiO2)3團(tuán)簇分子[11],為設(shè)計新型納米M/SiO2復(fù)合材料提供有力的的理論依據(jù).

    1 計算方法

    計算采用密度泛函方法,選用廣義梯度近似(GGA)泛函,利用Dmol3計算程序包完成[12].在電子結(jié)構(gòu)計算中,采取全電子和包括d極化的雙數(shù)基組,交換關(guān)聯(lián)相互作用采用GGA-PW91方法自洽場,收斂標(biāo)準(zhǔn)為2.625 5×10-2kJ·mol-1.為了加速自洽場收斂我們使用了DIIS方法,軌道計算中使用的Smearing標(biāo)準(zhǔn)為13.125 kJ·mol-1.在幾何優(yōu)化過程中,力的收斂標(biāo)準(zhǔn)是52.5 kJ·mol-1·nm-1,位移收斂標(biāo)準(zhǔn)5.0×10-4nm,能量收斂標(biāo)準(zhǔn)為2.625 5×10-2kJ·mol-1,整個計算過程選取fine精確的計算標(biāo)準(zhǔn).

    為了驗證方法的可靠性,在相同條件下,首先計算了O2、Si2、Fe2、Co2、Ni2二聚體以及Si-O、Fe-Si、Co-Si、Ni-Si的鍵長,并與已有的實驗值或理論報道相對照,如表1所示,計算結(jié)果與已有報道吻合較好,表明我們所選計算方法是合理的.

    表1 O2、Si2、Fe2、Co2、Ni2二聚體及Si-O、Fe-Si、Co-Si、Ni-Si的鍵長(單位:?)Table 1 Bond length(in ?)of O2,Si2,F(xiàn)e2,Co2,Ni2and Si-O,F(xiàn)e-Si,Co-Si,Ni-Si

    2 結(jié)果與討論

    2.1 幾何結(jié)構(gòu)

    2.1.1 M(SiO2)3(M=Fe,Co,Ni)

    圖1給出了(SiO2)3和M(SiO2)3(M=Fe,Co,Ni)團(tuán)簇的基態(tài)結(jié)構(gòu),其中,所有幾何結(jié)構(gòu)圖中還標(biāo)出了相應(yīng)鍵長和原子的Mulliken電荷分布(括號中數(shù)值).如圖1所示,M(SiO2)3(M=Fe,Co,Ni)團(tuán)簇有著非常相似的基態(tài)結(jié)構(gòu),單個過渡金屬原子均吸附在(SiO2)3末端的Si原子上,這是由于中間的Si原子是滿配位的,而末端的兩個Si原子配位數(shù)未滿,因此,末端Si原子上的懸掛鍵為過渡金屬Fe,Co,Ni原子的吸附提供了合適的位置.有趣的是單個過渡金屬原子在(SiO2)3團(tuán)簇上吸附的方式和單個貴金屬原子在(SiO2)3團(tuán)簇上吸附方式一致[6-8].

    從圖1可以看出,F(xiàn)e-Si,Co-Si,Ni-Si的鍵長分別是2.177 ?,2.144 ?,2.130 ?,F(xiàn)e-O,Co-O,Ni-O的鍵長分別是1.818 ?,1.826 ?,1.809 ?,M(SiO2)3(M=Fe,Co,Ni)團(tuán)簇中與過渡金屬原子相鄰的末端Si-O鍵長分別是1.661 ?,1.649 ?,1.630 ?,較單純(SiO2)3團(tuán)簇末端Si-O鍵長(1.520 ?)均有所增加.在自由團(tuán)簇中,由于處在不等價空間位置的原子感受到不同的勢場,一部分原子將失去電荷,另一部分原子將得到電荷,從而出現(xiàn)電荷轉(zhuǎn)移現(xiàn)象.Fe,Co,Ni分別失去了0.378 e,0.250 e,0.261 e,這些電荷主要引起與其臨近的Si,O原子電荷的重新分布.

    2.1.2 M2(SiO2)3(M=Fe,Co,Ni)

    當(dāng)?shù)诙€過渡金屬原子被置于M(SiO2)3團(tuán)簇的不同位置,經(jīng)幾何優(yōu)化后發(fā)現(xiàn),能量最低態(tài)中該原子吸附在原有的金屬原子上,而且,三種團(tuán)簇的基態(tài)結(jié)構(gòu)非常相似(如圖1所示).其中,F(xiàn)e-Fe,Co-Co,Ni-Ni鍵長分別是2.192 ?,2.226 ?,2.239 ?,F(xiàn)e-Si,Co-Si,Ni-Si的平均鍵長為2.391 ?,2.362 ?,2.344 ?.Fe2(SiO2)3中第二個Fe原子失去0.024 e,相對于Fe(SiO2)3團(tuán)簇Fe2(SiO2)3中第一個Fe原子失去的電荷數(shù)幾乎沒變.Co2(SiO2)3中第二個Co原子得到0.033 e,Ni2(SiO2)3中第二個Ni原子得到較少的電子(0.017 e),在Ni2(SiO2)3中參與作用的主要仍是第一個Ni原子,該原子失去更多的電子,表明Fe,Co原子相對于Ni原子更易于在(SiO2)3團(tuán)簇上聚集.

    2.1.3 M3(SiO2)3(M=Fe,Co,Ni)

    如圖2所示,F(xiàn)e3(SiO2)3和Co3(SiO2)3團(tuán)簇的基態(tài)構(gòu)型相似,第三個過渡金屬原子與前兩個過渡金屬原子相結(jié)合,三個過渡金屬原子形成了一個三角形吸附在(SiO2)3團(tuán)簇上.不同的是Ni3(SiO2)3中的第三個Ni原子吸附在另一端的Si原子上,進(jìn)一步表明二氧化硅團(tuán)簇上Fe,Co原子比Ni原子更易于形成“島膜”,這是由于Ni原子的外層電子排布為3s23p63d10的閉殼層排布,不利于與更多金屬原子參加團(tuán)聚作用,而Fe,Co原子的外層電子排布均為開殼層排布,易于發(fā)生電子間的相互轉(zhuǎn)移,從而促進(jìn)了更多外來原子的聚集.

    2.1.3 混作和套種 混作和套種是指在同一個有機(jī)農(nóng)作物種植環(huán)境中多樣化種植有機(jī)農(nóng)作物,以此防治害蟲對其侵害。例如,在棉花田里種植高粱等能控制棉鈴蟲產(chǎn)卵的農(nóng)作物,以實現(xiàn)集中誘殺害蟲的目的,保護(hù)棉花作物不受棉鈴蟲的侵害。

    圖2 M3(SiO2)3(M=Fe,Co,Ni)團(tuán)簇的基態(tài)結(jié)構(gòu)(黑色球為O,灰色球為Si)Fig.2 Ground state structures of M3(SiO2)3(M=Fe,Co,Ni)clusters

    為了研究團(tuán)簇的穩(wěn)定性,我們給出Mn(SiO2)3(M=Fe,Co,Ni;n=1-3)團(tuán)簇的平均結(jié)合能,平均結(jié)合能的定義是團(tuán)簇的結(jié)合能除以團(tuán)簇的總原子數(shù),其表達(dá)式:

    由表2可以看出,三種不同過渡金屬元素和不同原子數(shù)目的原子吸附在氧化硅團(tuán)簇上構(gòu)成的復(fù)合體幾乎具有相近的平均結(jié)合能(4.810 eV~5.354 eV),這與它們具有非常相似的幾何結(jié)構(gòu)有關(guān),表明氧化硅是一種很好的用于負(fù)載過渡金屬“島膜”的載體材料.

    表2 Mn(SiO2)3(M=Fe,Co,Ni;n=1-3)團(tuán)簇的優(yōu)選分裂途徑、分裂能(Ed)、平均結(jié)合能(Eb)、能隙、團(tuán)簇總磁矩和金屬原子磁矩(Mn)Table 2 Optimal disaggregated way,disaggregated energy,average binding energy,energy gap and moment of Mn(SiO2)3(M=Fe,Co,Ni;n=1-3)clusters

    通過分析團(tuán)簇的分裂途徑及其產(chǎn)物,便于我們理解表面吸附或沉積等現(xiàn)象.表2給出了Mn(SiO2)3(M=Fe,Co,Ni;n=1-3)團(tuán)簇的優(yōu)選分裂途徑及其相應(yīng)的分裂能,分裂能定義:如果分裂途徑為AmBn→Am+Bn,則分裂能定義為Ed=E(Am)+E(Bn)-E(AmBn),其中,E(Am)、E(Bn)和E(AmBn)均代表相應(yīng)團(tuán)簇的基態(tài)能量.Mn(SiO2)3(M=Fe,Co,Ni;n=1-3)團(tuán)簇的優(yōu)選分裂產(chǎn)物中均包含有(SiO2)3團(tuán)簇,另一種產(chǎn)物是過渡金屬原子的聚集體,這表明鏈狀的(SiO2)3團(tuán)簇具有很好的穩(wěn)定性,這與之前的文獻(xiàn)報道相一致[4].而且,隨著金屬原子數(shù)目的增加,聚集體Mn內(nèi)部的團(tuán)聚力不斷增加,這將有利于過渡金屬“島膜”的生長.

    團(tuán)簇一個很重要的電子性質(zhì)就是其能隙,即電子最高占據(jù)分子軌道(HOMO)與最低未占據(jù)分子軌道(LUMO)之間的能量差,如表2所示,除Con(SiO2)3(n=1-3)之外,Mn(SiO2)3(M=Fe,Ni;n=1-3)團(tuán)簇隨著金屬原子數(shù)目的增加,能隙不斷減小,發(fā)生紅移現(xiàn)象.而且,從整體看,該復(fù)合團(tuán)簇的能隙介于0.786 eV~1.131 eV,恰好位于近紅外光譜范圍內(nèi),與已有的貴金屬吸附在(SiO2)n團(tuán)簇上結(jié)果一致[6-8],說明Mn(SiO2)3(M=Fe,Co,Ni;n=1-3)團(tuán)簇可以吸收近紅外光,并將光轉(zhuǎn)化為熱量來殺死癌細(xì)胞,可能是一種治療腫瘤和癌癥的新型功能材料,同時,F(xiàn)e,Co,Ni均是人體不可缺少的微量元素,我們的計算結(jié)果進(jìn)一步肯定了過渡金屬—二氧化硅磁性復(fù)合材料在醫(yī)學(xué)界被用作光動力靶向治療的可觀前景.

    為了進(jìn)一步分析過渡金屬在(SiO2)n團(tuán)簇上聚集的機(jī)理,圖3和圖4分別給出了Mn(SiO2)3(M=Fe,Co,Ni;n=1-3)團(tuán)簇分子的前線軌道分布,由圖可知,電子云主要密集在過渡金屬原子周圍,這些原子將有利于與外來原子發(fā)生相互作用,因此,F(xiàn)e,Co,Ni原子均在(SiO2)n團(tuán)簇上發(fā)生了聚集現(xiàn)象,與前面的基態(tài)結(jié)構(gòu)圖保持一致.但是,Ni原子在二氧化硅上的吸附狀態(tài)卻不同于Fe和Co,由圖4可知Ni2(SiO2)3的HOMO軌道發(fā)生凹陷,即Ni原子的d軌道發(fā)生了畸變,此時,第三個Ni原子主要受懸掛鍵的影響吸附到另一端的Si原子上.總體來看,過渡金屬原子在(SiO2)n團(tuán)簇上的吸附不僅受懸掛鍵的影響,還受過渡金屬原子外層電子排布和分子前線軌道分布的影響.

    為了進(jìn)一步研究該復(fù)合體系的電子性質(zhì),計算了Mn(SiO2)3(M=Fe,Co,Ni;n=1-3)團(tuán)簇的分波態(tài)密度,圖5是Con(SiO2)3(n=1-3)團(tuán)簇的自旋分波態(tài)密度,從總體來看,在費米面附近起主要貢獻(xiàn)的是p、d軌道,而且,隨著過渡金屬原子數(shù)目的增加,d軌道貢獻(xiàn)不斷增強(qiáng),p軌道卻顯著下降,表明該分子軌道主要是由過渡金屬原子的3d軌道提供的.

    圖3 Fen(SiO2)3和Con(SiO2)3(n=1-3)團(tuán)簇的HOMO和LUMO軌道Fig.3 HOMO和LUMO orbitals of Fen(SiO2)3和Con(SiO2)3(n=1-3)clusters

    圖4 Nin(SiO2)3(n=1-3)團(tuán)簇的HOMO和LUMO軌道Fig.4 HOMO和LUMO orbitals of Nin(SiO2)3(n=1-3)clusters

    2.3 磁學(xué)性質(zhì)

    對于含有過渡金屬的團(tuán)簇,我們關(guān)心的另外一個問題是團(tuán)簇的磁性.表2給出了基態(tài)Mn(SiO2)3(M=Fe,Co,Ni;n=1-3)團(tuán)簇的總磁矩和過渡金屬原子的磁矩,原來不具有磁性的氧化硅團(tuán)簇在吸附有過渡金屬后變成了很好的磁性團(tuán)簇,表明Mn(SiO2)3(M=Fe,Co,Ni;n=1-3)團(tuán)簇內(nèi)部出現(xiàn)了顯著的自旋極化現(xiàn)象,而且,由表2可知團(tuán)簇的磁矩主要是由過渡金屬原子提供的.

    Fe(SiO2)3中Fe原子磁矩近似于鐵塊體材料所具有得磁矩(2.22 μB[13]),F(xiàn)e2(SiO2)3中兩個Fe原子磁矩均大于鐵塊體材料的磁矩,使該團(tuán)簇具有6 μB的磁矩,實驗已經(jīng)表明由鐵磁性原子(Fe,Co,Ni)構(gòu)成的團(tuán)簇體系的磁矩遠(yuǎn)大于其塊體材料的磁矩[14-15].Fe3(SiO2)3的總磁矩為4 μB,第三個Fe原子出現(xiàn)了反鐵磁性(磁矩為-3.027 μB).Co2(SiO2)3和Co3(SiO2)3團(tuán)簇也具有較高的磁矩分別是4 μB、7 μB,其中,該復(fù)合團(tuán)簇中每個Co原子的磁矩均大于鈷塊體材料的磁矩(1.72 μB[13]).Fe2(SiO2)3和Co3(SiO2)3具有較大的總磁矩可能是由于過渡金屬原子的d軌道之間發(fā)生顯著耦合,加強(qiáng)了電子的自旋極化,與圖3中展示的電子云分布情況保持一致.

    值得注意的是Ni(SiO2)3團(tuán)簇出現(xiàn)磁矩淬滅現(xiàn)象,主要原因是其中的Ni原子磁矩為零,而且,該團(tuán)簇中其它原子的磁矩也均為零,Ni2(SiO2)3和Ni3(SiO2)3團(tuán)簇的磁矩都是2 μB,吸附不同數(shù)目的Ni原子卻導(dǎo)致相同的磁矩,主要是因為Ni3(SiO2)3中吸附在末端的第三個Ni原子發(fā)生了磁矩淬滅現(xiàn)象.

    圖5 Con(SiO2)3(n=1-3)團(tuán)簇的分波態(tài)密度Fig.5 Partial density of states for Con(SiO2)3(n=1-3)clusters

    為了進(jìn)一步研究過渡金屬-二氧化硅復(fù)合團(tuán)簇的磁學(xué)性質(zhì),在考慮不同構(gòu)型和不同自旋多重度的情況下,給出了Mn(SiO2)3(M=Fe,Co,Ni;n=2-3)團(tuán)簇的鐵磁(FM)和反鐵磁(AFM)能量差,對于M2(SiO2)3(M=Fe,Co,Ni)團(tuán)簇來說,分別是0.638 eV、0.490 eV、1.074 eV,M3(SiO2)3(M=Fe,Co,Ni)相應(yīng)能量差為0.249 eV、0.265 eV、0.299 eV,這表明M2(SiO2)3(M=Fe,Co,Ni)相對于M3(SiO2)3(M=Fe,Co,Ni)團(tuán)簇的鐵磁性更加穩(wěn)定.

    3 結(jié)論

    基于密度泛函理論中的廣義梯度近似系統(tǒng)地研究了Mn(SiO2)3(M=Fe,Co,Ni;n=1-3)團(tuán)簇的幾何結(jié)構(gòu)、電子性質(zhì)和磁學(xué)性質(zhì).通過對基態(tài)結(jié)構(gòu)、電子性質(zhì)和磁學(xué)性質(zhì)的分析,主要結(jié)論如下:

    1)Fe,Co原子相對于Ni原子更易于在(SiO2)3團(tuán)簇上聚集,而且,過渡金屬原子在(SiO2)n團(tuán)簇上的吸附不僅受懸掛鍵的影響,還受過渡金屬原子外層電子排布和分子前線軌道分布的影響.

    2)通過平均結(jié)合能、優(yōu)選分裂途徑及分裂能計算,從理論上證明了鏈狀的(SiO2)3團(tuán)簇具有很好的穩(wěn)定性,是一種很好的用于負(fù)載過渡金屬“島膜”的載體材料.

    3)Mn(SiO2)3(M=Fe,Ni;n=1-3)團(tuán)簇的能隙介于0.786 eV~1.131 eV,恰好位于近紅外光譜范圍內(nèi),表明該復(fù)合材料可能是一種治療腫瘤和癌癥的新型醫(yī)用功能材料.

    4)Mn(SiO2)3(M=Fe,Co,Ni;n=1-3)團(tuán)簇內(nèi)部出現(xiàn)了顯著的自旋極化現(xiàn)象,原來不具有磁性的氧化硅團(tuán)簇在吸附有過渡金屬后變成了很好的磁性團(tuán)簇,該磁矩主要局域在過渡金屬原子周圍,F(xiàn)e2(SiO2)3和Co3(SiO2)3具有相對較大的磁矩(6 μB、7 μB).

    [1] Wang L S,Nicholas J B,Dupuis M,et al.Si3Oy(y=1-6)clusters:Models for oxidation of silicon surfaces and defect sites in bulk oxide materials[J].Phys Rev Lett,1997,78:4450-4453.

    [2] Li S,Silvers S J,El-Shall M S.Surface oxidation and luminescence properties of weblike agglomeration of silicon nanocrystals produced by a laser vaporization-controlled condensation technique[J].J Phys Chem B,1997,101(10):1794-1802.

    [3] Harkless J A,Stilinger D K,Stillinger F H.Structures and energies of SiO2clusters[J].J Phys Chem,1996,100(4):1098-1103.

    [4] Nayak S K,Rao B K,Khanna S N,et al.Atomic and electronic structure of neutral and charged SinOmclusters[J].J Chem Phys,1998,109(4):1245-1250.

    [5] Zhi L L,Zhao G F,Guo L J,Jing Q.Structural,electronic,and vibrational properties of water molecules adsorbed on silica clusters[J].Phys Rev B,2008,77:235435/1-7.

    [6] Zhao G F,Zhi L L,Guo L J,Zeng Z.The structural and electronic properties of Ag-adsorbed(SiO2)n(n=1-7)clusters[J].J Chem Phys,2007,127:234705/1-7.

    [7] Sun Q,Wang Q,Rao B K,et al.Electronic structure and bonding of Au on a SiO2cluster:A nanobullet for tumors[J].Phys Rev Lett,2004,93(18):186803/1-4.

    [8] 孫建敏,趙高峰,王獻(xiàn)偉,等.Cu吸附(SiO2)n(n=1-8)團(tuán)簇幾何結(jié)構(gòu)和電子性質(zhì)的密度泛函研究[J].物理學(xué)報,2010,59(11):7830-7837.

    [9] Li Y L,Yang M L.二氧化硅磁性復(fù)合微球的制備[D].西安:西北大學(xué),2010.

    [10] Yeshchenko O A,Dmitruk I M,Alexeenko A A,Dmytruk A M.Optical properties of sol-gel fabricated Ni/SiO2glass nanocomposites[J].J Phys Chem Solids,2008,69:1615-1622.

    [11] Wang L S,Desai S R,Wu H,Nicholas J B.Small silicon oxide clusters:Chains and rings[J].Z Phys D,1997,40:36-39.

    [12] Delley B.From molecules to solids with the DMol3approach[J].J Chem Phys,2000,113(18):7756-7764.

    [13] 葛桂賢,羅有華,張建瑋.FMBen(FM=Fe,Co,Ni;n=1-12)團(tuán)簇結(jié)構(gòu)和電子性質(zhì)[J].物理化學(xué)學(xué)報,2008,24(10):1891-1896.

    [14] Billas I M L,Chatelain A,de Heer Walt A.Magnetism from the atom to the bulk for iron,cobalt and nickel clusters[J]. Science,1994,265:1682-1684.

    [15] Apsel S E,Emment J W,Deng J,Bloomfield L A.Surface-enhanced magnetism in nickel clusters[J].Phys Rev Lett,1996,76:1441-1444.

    [16] Lide D R.CRC handbook of chemistry and physics[M].CRC,Boca Raton,F(xiàn)L,2000.

    [17] Kant A,Strauss B.Dissociation energies of diatomic molecules of the transition elements.II:Titanium,chromium,manganese,and cobalt[J].J Chem Phys,1964,41(12):3806-3808.

    [18] Hales D A,Su C X,Li L,Armenstout P B.Collision-induced dissociation of Co+n(n=2-18)with Xe:Bond energies of cationic and neutral cobalt clusters,dissociation pathways,and structures[J].J Chem Phys,1994,100:1049-1057.

    [19] Purdum H,Montano P A,Shenoyand G K,Morrison T.Extended-x-ray-absorption-fine-structure study of small Fe molecules isolated in solid neon[J].Phys Rev B,1982,25:4412-4417.

    [20] Wu Z J,Su Z M.Electronic structures and chemical bonding in transition metal monosilicides MSi(M=3 d,4 d,5 d elements)[J].J Chem Phys,2006,124:184306/1-15.

    [21] Lindholm N F,Burgh D J,Rothschopf G K,Sickaf-oose S M,Morse M D.Optical spectroscopy of jet-cooled NiSi[J].J Chem Phys,2003,118:2190-2196.

    Structural,Photoelectric and Magnetic Properties of Mn(SiO2)3(M=Fe,Co,Ni;n=1-3)Clusters

    ZHI Lili1,2,LI Yanqing2,3,YANG Lianhong2,ZHAO Gaofeng4
    (1.Xinjiang Laboratory of Phase Transitions and Microstructures of Condensed Matters,Yili Normal University,Yining 835000,China;2.Department of Physics,Changji College,Xinjiang,Changji 831100,China;3.State Key Laboratory of Crystal Materials,Shandong University,Jinan 250100,China;4.School of Physics and Electronics,Henan University,Kaifeng 475004,China)

    Equilibrium geometries,electronic and magnetic properties of Mn(SiO2)3(M=Fe,Co,Ni;n=1-3)clusters are systematically studied employing density functional theory with a generalized gradient approximation.It shows that Fe and Co atoms are easier to congregate on(SiO2)3cluster than Ni atoms.It is found that stabler silica is an excellent matrix materials to carry islands of transition-metals.Energy gaps of Mn(SiO2)3(M=Fe,Co,Ni;n=1-3)clusters lie in near infrared radiation region.In analysis of magnetism,it is found that their magnetic moments are mainly located on transition-metal atoms.Fe2(SiO2)3and Co3(SiO2)3have greater magnetic moments,owing to coupling between d orbits of transition-metal atoms.Energy gap and magnetic property affirm a considerable foreground of magnetic-mulriple silica used for photodynamic target therapy in medical stage.

    Mn(SiO2)3cluster;geometrical structure;photoelectric property;magnetism

    date: 2013-12-03;Revised date: 2014-04-06

    O561

    A

    2013-12-03;

    2014-04-06

    新疆維吾爾自治區(qū)高??蒲杏媱澷Y助(XJEDU2013S42);新疆凝聚態(tài)相變與微結(jié)構(gòu)實驗室開放課題基金(XJDX0912-2010-05和XJDX0912-2010-07)及昌吉學(xué)院院級課題(2013YJYB002)資助項目

    智麗麗(1984-),女,研究生,講師,主要從事團(tuán)簇結(jié)構(gòu)和電子性質(zhì)研究,E-mail:zhilili2010@sina.com

    1001-246X(2014)06-0727-08

    猜你喜歡
    磁矩基態(tài)鍵長
    陰離子調(diào)控錳基鈣鈦礦中Mn─O的鍵長和磁性
    一類非線性Choquard方程基態(tài)解的存在性
    擬相對論薛定諤方程基態(tài)解的存在性與爆破行為
    一類反應(yīng)擴(kuò)散方程的Nehari-Pankov型基態(tài)解
    非線性臨界Kirchhoff型問題的正基態(tài)解
    密度泛函理論研究鎘的二鹵化合物分子的結(jié)構(gòu)和振動頻率
    淺議鍵能與鍵長的關(guān)系
    CoFeB/MgO磁隧道結(jié)的低電流密度磁矩翻轉(zhuǎn)特性
    兩種計算帶電輕子磁矩的嘗試
    河南科技(2014年23期)2014-02-27 14:18:52
    苯并三氧化呋咱的晶體結(jié)構(gòu)
    国产一区二区在线av高清观看| 亚洲乱码一区二区免费版| 国产99白浆流出| 一区二区三区国产精品乱码| 亚洲中文日韩欧美视频| 国产1区2区3区精品| 中文字幕av在线有码专区| 777久久人妻少妇嫩草av网站| 757午夜福利合集在线观看| 亚洲欧美日韩东京热| 在线a可以看的网站| 亚洲真实伦在线观看| 午夜日韩欧美国产| 精品国产亚洲在线| 又粗又爽又猛毛片免费看| 亚洲国产看品久久| 露出奶头的视频| 日本免费a在线| 大型av网站在线播放| 99在线视频只有这里精品首页| 亚洲av电影不卡..在线观看| 国产激情久久老熟女| 91大片在线观看| 啦啦啦韩国在线观看视频| 91在线观看av| 免费无遮挡裸体视频| 国产高清视频在线播放一区| 波多野结衣高清无吗| 精品久久久久久,| 宅男免费午夜| 一个人免费在线观看的高清视频| 亚洲狠狠婷婷综合久久图片| 亚洲九九香蕉| av视频在线观看入口| 亚洲五月天丁香| 别揉我奶头~嗯~啊~动态视频| 国产精品亚洲av一区麻豆| 日本在线视频免费播放| 日日摸夜夜添夜夜添小说| 色综合婷婷激情| 一本大道久久a久久精品| 精品一区二区三区视频在线观看免费| 2021天堂中文幕一二区在线观| 此物有八面人人有两片| 国产蜜桃级精品一区二区三区| 制服丝袜大香蕉在线| 三级毛片av免费| 国内精品久久久久精免费| 麻豆成人av在线观看| 日韩 欧美 亚洲 中文字幕| 嫩草影院精品99| 亚洲一卡2卡3卡4卡5卡精品中文| 两个人看的免费小视频| 一二三四在线观看免费中文在| 一级a爱片免费观看的视频| 天天添夜夜摸| 国产成人av教育| 啦啦啦韩国在线观看视频| 在线观看www视频免费| 在线观看一区二区三区| 亚洲人成77777在线视频| 男女之事视频高清在线观看| 成人午夜高清在线视频| 香蕉国产在线看| 久久久国产成人精品二区| 正在播放国产对白刺激| 韩国av一区二区三区四区| 精品少妇一区二区三区视频日本电影| 欧美+亚洲+日韩+国产| 久久亚洲精品不卡| 国产精品免费一区二区三区在线| 久久久久久久午夜电影| 午夜福利免费观看在线| 51午夜福利影视在线观看| 欧美人与性动交α欧美精品济南到| 免费看日本二区| 一本精品99久久精品77| 久久久久久久久中文| 黄色成人免费大全| 午夜福利在线在线| 亚洲精品美女久久av网站| 久久久久久久精品吃奶| 99在线视频只有这里精品首页| 中文在线观看免费www的网站 | 色老头精品视频在线观看| 一边摸一边做爽爽视频免费| 久久久久亚洲av毛片大全| 色噜噜av男人的天堂激情| 深夜精品福利| 日韩 欧美 亚洲 中文字幕| 12—13女人毛片做爰片一| 亚洲午夜精品一区,二区,三区| 人妻丰满熟妇av一区二区三区| 亚洲一区高清亚洲精品| 欧美人与性动交α欧美精品济南到| 久久久精品欧美日韩精品| 99精品欧美一区二区三区四区| 国产精品久久久av美女十八| 欧美日韩黄片免| 亚洲一区二区三区不卡视频| 小说图片视频综合网站| 亚洲真实伦在线观看| 99热只有精品国产| 欧美zozozo另类| 久久这里只有精品19| 国产亚洲欧美98| 亚洲avbb在线观看| 成人18禁在线播放| 亚洲天堂国产精品一区在线| 特级一级黄色大片| 精品国产亚洲在线| 亚洲国产看品久久| 国产一区二区三区在线臀色熟女| 香蕉国产在线看| 色尼玛亚洲综合影院| 成人欧美大片| 精品国产亚洲在线| 久久人人精品亚洲av| 久久久水蜜桃国产精品网| 久久天躁狠狠躁夜夜2o2o| 午夜激情福利司机影院| 制服人妻中文乱码| 欧美日韩瑟瑟在线播放| 欧美绝顶高潮抽搐喷水| 熟女电影av网| 在线观看一区二区三区| 久久中文字幕一级| 1024手机看黄色片| 亚洲精品国产一区二区精华液| 中文字幕高清在线视频| 国产一级毛片七仙女欲春2| avwww免费| 久久中文字幕一级| 久久这里只有精品19| 精品久久久久久久久久免费视频| 免费无遮挡裸体视频| 午夜福利18| 长腿黑丝高跟| 一个人观看的视频www高清免费观看 | 亚洲欧美精品综合久久99| 黑人欧美特级aaaaaa片| 亚洲一区中文字幕在线| 女人高潮潮喷娇喘18禁视频| 黄色视频不卡| 香蕉国产在线看| 色尼玛亚洲综合影院| 欧美性长视频在线观看| 亚洲精品在线观看二区| 十八禁网站免费在线| 高清毛片免费观看视频网站| 欧洲精品卡2卡3卡4卡5卡区| 国产精品亚洲美女久久久| 手机成人av网站| 此物有八面人人有两片| 国产av不卡久久| 又大又爽又粗| 精品欧美一区二区三区在线| 日日夜夜操网爽| 欧美日韩瑟瑟在线播放| 日韩精品青青久久久久久| 国产av不卡久久| 两性夫妻黄色片| 最近最新中文字幕大全电影3| 国产精品98久久久久久宅男小说| 18禁黄网站禁片免费观看直播| 成在线人永久免费视频| 黄片大片在线免费观看| www.精华液| 又黄又粗又硬又大视频| 日韩精品青青久久久久久| 麻豆国产97在线/欧美 | 免费一级毛片在线播放高清视频| 成人特级黄色片久久久久久久| 国产av一区二区精品久久| 免费看美女性在线毛片视频| 男人的好看免费观看在线视频 | 两个人免费观看高清视频| 免费观看人在逋| 欧美色欧美亚洲另类二区| 真人做人爱边吃奶动态| 亚洲精品一卡2卡三卡4卡5卡| 精品电影一区二区在线| 亚洲男人天堂网一区| 久久中文字幕人妻熟女| 三级毛片av免费| 91老司机精品| 在线观看免费视频日本深夜| 日本 av在线| 日日干狠狠操夜夜爽| 中出人妻视频一区二区| 一区二区三区国产精品乱码| 老汉色∧v一级毛片| 精品乱码久久久久久99久播| 国产精华一区二区三区| 精品久久久久久久久久久久久| 国产精品98久久久久久宅男小说| 国产高清激情床上av| 亚洲国产精品成人综合色| 老司机福利观看| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品一区二区三区四区久久| 十八禁人妻一区二区| 国产精品永久免费网站| 曰老女人黄片| 在线观看一区二区三区| 国产成人精品久久二区二区91| aaaaa片日本免费| 香蕉久久夜色| 国产亚洲精品第一综合不卡| 亚洲av电影在线进入| 国产成人aa在线观看| 91九色精品人成在线观看| 啦啦啦免费观看视频1| 久久精品成人免费网站| 日本三级黄在线观看| 亚洲 欧美 日韩 在线 免费| 美女高潮喷水抽搐中文字幕| 两个人视频免费观看高清| 欧美成人免费av一区二区三区| 超碰成人久久| 国产不卡一卡二| 在线观看美女被高潮喷水网站 | 国内精品久久久久久久电影| 人妻丰满熟妇av一区二区三区| 亚洲精品在线美女| 哪里可以看免费的av片| 午夜激情福利司机影院| 亚洲一区二区三区不卡视频| 国产精品久久久av美女十八| 岛国在线免费视频观看| 日韩欧美在线二视频| 欧美精品亚洲一区二区| 日本 av在线| 麻豆av在线久日| 国产aⅴ精品一区二区三区波| 久久精品91蜜桃| 久久精品影院6| 在线观看免费午夜福利视频| 真人做人爱边吃奶动态| 欧美又色又爽又黄视频| 黄色视频,在线免费观看| 亚洲人成伊人成综合网2020| 国产亚洲精品一区二区www| 日本一二三区视频观看| 亚洲一区高清亚洲精品| www日本在线高清视频| 婷婷六月久久综合丁香| 一边摸一边抽搐一进一小说| 最近最新中文字幕大全电影3| 黄色 视频免费看| 久久精品aⅴ一区二区三区四区| a级毛片在线看网站| 亚洲天堂国产精品一区在线| 久久婷婷成人综合色麻豆| 久久香蕉国产精品| 国产成人一区二区三区免费视频网站| 国产激情欧美一区二区| 国产爱豆传媒在线观看 | 国产精品野战在线观看| www.自偷自拍.com| 少妇裸体淫交视频免费看高清 | 亚洲最大成人中文| 久久这里只有精品中国| 午夜福利在线观看吧| 久久午夜综合久久蜜桃| 亚洲熟女毛片儿| 亚洲熟妇中文字幕五十中出| 欧美乱妇无乱码| 真人做人爱边吃奶动态| 亚洲成人精品中文字幕电影| 人成视频在线观看免费观看| 日韩欧美国产在线观看| 久久精品91蜜桃| 黄片小视频在线播放| 国产精品美女特级片免费视频播放器 | 中出人妻视频一区二区| 无限看片的www在线观看| 成人三级黄色视频| 午夜福利在线观看吧| 亚洲在线自拍视频| 国产爱豆传媒在线观看 | 国产激情偷乱视频一区二区| 亚洲全国av大片| 丝袜美腿诱惑在线| 精品日产1卡2卡| 国产精品亚洲一级av第二区| 五月玫瑰六月丁香| x7x7x7水蜜桃| 欧美日韩瑟瑟在线播放| 1024香蕉在线观看| 国产成人aa在线观看| 91九色精品人成在线观看| 男人舔女人下体高潮全视频| 嫁个100分男人电影在线观看| 午夜精品一区二区三区免费看| 日韩精品青青久久久久久| 我要搜黄色片| 国产av一区在线观看免费| 12—13女人毛片做爰片一| 91成年电影在线观看| 亚洲avbb在线观看| 亚洲 欧美 日韩 在线 免费| 欧美性猛交黑人性爽| 91麻豆av在线| 小说图片视频综合网站| 国产av不卡久久| 久久久久久亚洲精品国产蜜桃av| 精品福利观看| www.999成人在线观看| e午夜精品久久久久久久| 手机成人av网站| 久久婷婷人人爽人人干人人爱| 国产亚洲av嫩草精品影院| 国产亚洲精品av在线| 国产精品香港三级国产av潘金莲| 嫁个100分男人电影在线观看| 国产成人精品久久二区二区免费| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品精品国产色婷婷| 三级毛片av免费| 亚洲专区字幕在线| 丰满人妻一区二区三区视频av | 亚洲人成网站在线播放欧美日韩| 亚洲成人国产一区在线观看| tocl精华| 91在线观看av| 欧美黑人精品巨大| 1024手机看黄色片| 国产精品电影一区二区三区| 国产精品影院久久| 视频区欧美日本亚洲| 国产精品久久电影中文字幕| 999精品在线视频| 丝袜美腿诱惑在线| 国产又黄又爽又无遮挡在线| 99国产极品粉嫩在线观看| 午夜福利在线观看吧| 麻豆一二三区av精品| 巨乳人妻的诱惑在线观看| 在线观看免费日韩欧美大片| 国产熟女xx| 国产精品免费一区二区三区在线| 在线看三级毛片| 小说图片视频综合网站| 国产精品爽爽va在线观看网站| 亚洲国产中文字幕在线视频| 黄色 视频免费看| 亚洲七黄色美女视频| 国产精品电影一区二区三区| 美女午夜性视频免费| 18美女黄网站色大片免费观看| 看片在线看免费视频| 久久亚洲精品不卡| 99久久国产精品久久久| 亚洲精品一卡2卡三卡4卡5卡| 亚洲中文av在线| 亚洲精品国产一区二区精华液| 午夜福利高清视频| 国产精品久久久久久人妻精品电影| 日日夜夜操网爽| 国内揄拍国产精品人妻在线| 亚洲人与动物交配视频| 午夜激情av网站| 亚洲av日韩精品久久久久久密| 亚洲国产精品久久男人天堂| 久久久久亚洲av毛片大全| 国产精品 国内视频| 成人精品一区二区免费| avwww免费| 操出白浆在线播放| 欧美日韩精品网址| 视频区欧美日本亚洲| 亚洲国产中文字幕在线视频| 亚洲 国产 在线| 日韩大码丰满熟妇| 久久精品亚洲精品国产色婷小说| www.999成人在线观看| 亚洲avbb在线观看| www.自偷自拍.com| 日本一二三区视频观看| 天天一区二区日本电影三级| 黄片小视频在线播放| tocl精华| 午夜日韩欧美国产| 一进一出好大好爽视频| 久久久精品国产亚洲av高清涩受| 国语自产精品视频在线第100页| 白带黄色成豆腐渣| 亚洲精品美女久久av网站| 精品欧美一区二区三区在线| 精品福利观看| 亚洲熟妇熟女久久| 亚洲国产精品久久男人天堂| 神马国产精品三级电影在线观看 | 亚洲专区国产一区二区| 淫秽高清视频在线观看| 最近最新中文字幕大全电影3| 国产精品久久久av美女十八| 久久香蕉激情| 精品福利观看| 欧美日韩亚洲综合一区二区三区_| 九色国产91popny在线| 亚洲成人精品中文字幕电影| 真人一进一出gif抽搐免费| 日韩精品中文字幕看吧| av有码第一页| 中文字幕最新亚洲高清| 亚洲乱码一区二区免费版| 欧美成狂野欧美在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 18美女黄网站色大片免费观看| 国产伦在线观看视频一区| 床上黄色一级片| 日韩欧美国产在线观看| 国产在线观看jvid| 99热6这里只有精品| 黄色 视频免费看| 亚洲黑人精品在线| 一本综合久久免费| 国产一区二区激情短视频| 精品久久久久久久久久久久久| 国产精品久久久av美女十八| 久99久视频精品免费| 成人高潮视频无遮挡免费网站| 久久久国产精品麻豆| 国产精品一区二区免费欧美| 国产成人系列免费观看| 亚洲成人久久爱视频| 在线观看www视频免费| 色精品久久人妻99蜜桃| 亚洲一区二区三区不卡视频| 亚洲国产欧美人成| 精品久久久久久久久久免费视频| 男女下面进入的视频免费午夜| 老司机午夜十八禁免费视频| 亚洲精品中文字幕在线视频| a级毛片在线看网站| 不卡av一区二区三区| 一个人免费在线观看的高清视频| 亚洲在线自拍视频| 久久这里只有精品中国| 久久久国产成人精品二区| or卡值多少钱| 一本一本综合久久| 亚洲人成网站在线播放欧美日韩| 岛国在线免费视频观看| 亚洲国产日韩欧美精品在线观看 | 国产激情偷乱视频一区二区| avwww免费| 亚洲成a人片在线一区二区| 好男人电影高清在线观看| 制服诱惑二区| 熟女电影av网| 国产精品 国内视频| 成人高潮视频无遮挡免费网站| 中文字幕熟女人妻在线| 日日爽夜夜爽网站| 在线观看一区二区三区| av福利片在线| 亚洲熟女毛片儿| 毛片女人毛片| 亚洲熟女毛片儿| 久久精品成人免费网站| av免费在线观看网站| 国产伦人伦偷精品视频| 亚洲一区二区三区不卡视频| 天天添夜夜摸| 一区二区三区国产精品乱码| 亚洲国产欧美人成| 国产高清videossex| 男人舔女人下体高潮全视频| 欧美中文日本在线观看视频| 国产精品久久久久久亚洲av鲁大| 麻豆国产av国片精品| 国产精品99久久99久久久不卡| 亚洲九九香蕉| 在线国产一区二区在线| 国产av不卡久久| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美精品v在线| 狂野欧美白嫩少妇大欣赏| 亚洲成人免费电影在线观看| 亚洲九九香蕉| 午夜精品一区二区三区免费看| 亚洲国产精品sss在线观看| 国产熟女午夜一区二区三区| 国产乱人伦免费视频| 午夜亚洲福利在线播放| 日日夜夜操网爽| 99re在线观看精品视频| 1024手机看黄色片| 在线观看www视频免费| 国产高清有码在线观看视频 | 精品国产超薄肉色丝袜足j| 色播亚洲综合网| 香蕉久久夜色| av有码第一页| 亚洲欧洲精品一区二区精品久久久| 天天躁狠狠躁夜夜躁狠狠躁| 可以在线观看的亚洲视频| 亚洲欧美精品综合一区二区三区| 日韩欧美免费精品| 久久久久精品国产欧美久久久| 日韩欧美精品v在线| e午夜精品久久久久久久| 伊人久久大香线蕉亚洲五| 亚洲五月婷婷丁香| 香蕉久久夜色| videosex国产| 色在线成人网| 最近在线观看免费完整版| 亚洲最大成人中文| 国产一区二区激情短视频| 国产av麻豆久久久久久久| 国产精品九九99| 日韩欧美在线乱码| 久久中文字幕一级| 不卡av一区二区三区| 99国产精品一区二区蜜桃av| 成人国产综合亚洲| 一级片免费观看大全| 夜夜看夜夜爽夜夜摸| 别揉我奶头~嗯~啊~动态视频| 久久精品人妻少妇| 很黄的视频免费| 国产精品久久久久久亚洲av鲁大| 国产区一区二久久| 亚洲精品久久成人aⅴ小说| 国产成人啪精品午夜网站| 国产91精品成人一区二区三区| 精品一区二区三区av网在线观看| 久久精品影院6| 亚洲精品av麻豆狂野| 欧美性猛交╳xxx乱大交人| 亚洲av片天天在线观看| 极品教师在线免费播放| 欧美 亚洲 国产 日韩一| 啪啪无遮挡十八禁网站| 麻豆成人av在线观看| 熟女电影av网| 啪啪无遮挡十八禁网站| 国产亚洲精品久久久久久毛片| 欧美成狂野欧美在线观看| 一进一出抽搐gif免费好疼| 免费在线观看亚洲国产| 村上凉子中文字幕在线| 美女扒开内裤让男人捅视频| 国产私拍福利视频在线观看| 欧美一级a爱片免费观看看 | 夜夜夜夜夜久久久久| 国产不卡一卡二| 中文在线观看免费www的网站 | 亚洲欧美日韩高清在线视频| 国产精品一区二区三区四区久久| av在线天堂中文字幕| 三级毛片av免费| 男女午夜视频在线观看| www国产在线视频色| 精品久久久久久,| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久久久久黄片| 亚洲成av人片在线播放无| 亚洲七黄色美女视频| 国语自产精品视频在线第100页| 久久精品国产亚洲av高清一级| 国产精品免费视频内射| 国产真实乱freesex| 亚洲成人精品中文字幕电影| 国产精品 欧美亚洲| 亚洲精品美女久久久久99蜜臀| 中文字幕av在线有码专区| 欧美性猛交╳xxx乱大交人| 人人妻人人看人人澡| 欧美日韩乱码在线| 亚洲欧美精品综合一区二区三区| 青草久久国产| 久久久久久人人人人人| 久久午夜亚洲精品久久| ponron亚洲| 欧美色欧美亚洲另类二区| 日本五十路高清| 91麻豆av在线| 亚洲精品色激情综合| 久久精品国产综合久久久| 久久亚洲真实| 99久久久亚洲精品蜜臀av| 国产视频内射| 久久亚洲精品不卡| 国产探花在线观看一区二区| 中文字幕熟女人妻在线| 亚洲午夜精品一区,二区,三区| 久久天躁狠狠躁夜夜2o2o| 全区人妻精品视频| 国产精品美女特级片免费视频播放器 | 91大片在线观看| 久久精品综合一区二区三区| 一本综合久久免费| 最近最新免费中文字幕在线| 欧美一区二区精品小视频在线| 欧美一级毛片孕妇| 国产高清视频在线观看网站| 欧美乱码精品一区二区三区| 欧美zozozo另类| 三级国产精品欧美在线观看 | 级片在线观看| 一个人免费在线观看的高清视频| 成人国产综合亚洲| 1024视频免费在线观看| 亚洲成人久久性| xxxwww97欧美| 在线观看一区二区三区| 国产aⅴ精品一区二区三区波| 亚洲熟妇中文字幕五十中出| 波多野结衣巨乳人妻| 性色av乱码一区二区三区2|