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

    團(tuán)簇Co3NiB2極化率、偶極矩及態(tài)密度研究

    2022-01-26 04:48方志剛王智瑤鄭新喜秦渝毛智龍曾鑫漁朱依文王倩

    方志剛 王智瑤 鄭新喜 秦渝 毛智龍 曾鑫漁 朱依文 王倩

    摘 要:為探究團(tuán)簇Co3NiB2內(nèi)部結(jié)構(gòu)的相關(guān)狀況及極性強(qiáng)弱,基于拓?fù)鋵W(xué)原理和密度泛函理論,在B3LYP/lanl2dz水平下,對(duì)團(tuán)簇12種優(yōu)化構(gòu)型的極化率、偶極矩及態(tài)密度進(jìn)行深入研究,最終得出以下結(jié)論:通過(guò)對(duì)團(tuán)簇形變程度的排序分析可得,極化率張量對(duì)團(tuán)簇的幾何結(jié)構(gòu)具有很強(qiáng)的依賴性;分析團(tuán)簇對(duì)外場(chǎng)響應(yīng)程度的排序可以發(fā)現(xiàn),極化率各向異性不變量對(duì)團(tuán)簇的幾何結(jié)構(gòu)具有較強(qiáng)的依賴性,但依賴性不及極化率張量;從偶極矩角度分析發(fā)現(xiàn),團(tuán)簇所有優(yōu)化構(gòu)型均為極性分子,其中構(gòu)型4(4)的分子極性最強(qiáng),構(gòu)型4(2)的分子極性最弱;對(duì)態(tài)密度進(jìn)行分析發(fā)現(xiàn),團(tuán)簇的成鍵主要存在p-p、d-d-p、d-p-p-p、s-s-p-p-p和d-d-p-p-p 5種雜化方式,其中,p-p、d-d-p雜化作用較強(qiáng),d-p-p-p、s-s-p-p-p、d-d-p-p-p雜化作用較弱。此外,團(tuán)簇所有優(yōu)化構(gòu)型的A、B、D、E峰處的雜化情況完全一致,僅有C峰的雜化情況因構(gòu)型重態(tài)的不同而存在差異,說(shuō)明C峰處的雜化情況會(huì)在一定程度上受團(tuán)簇多重度的影響。

    關(guān)鍵詞:團(tuán)簇Co3NiB2;極化率;偶極矩;態(tài)密度

    中圖分類(lèi)號(hào):O64112

    文獻(xiàn)標(biāo)志碼:A

    物質(zhì)根據(jù)其內(nèi)在的原子排列方式可以劃分為兩種,即原子呈規(guī)則排列的晶體材料和原子呈無(wú)規(guī)則排列的非晶體材料。自從1934年德國(guó)物理學(xué)家Kramer等采用蒸發(fā)沉積法制備出非晶態(tài)薄膜[1],對(duì)于非晶態(tài)合金的合成與應(yīng)用研究便迅速開(kāi)展起來(lái)。而隨著現(xiàn)代工業(yè)技術(shù)的飛速發(fā)展和人們環(huán)保意識(shí)的不斷提高,如何得到清潔、可持續(xù)的能源載體來(lái)替代傳統(tǒng)化石燃料,以解決一次能源的匱乏和環(huán)境污染嚴(yán)重等問(wèn)題,已經(jīng)引起了大眾的密切關(guān)注。現(xiàn)如今,通過(guò)借助各種催化劑來(lái)電解水以獲得可持續(xù)發(fā)展的清潔能源(氫氣),是科學(xué)家們慣用的方法之一[2-3]。在經(jīng)過(guò)大量實(shí)驗(yàn)驗(yàn)證后發(fā)現(xiàn),處于熱力學(xué)亞穩(wěn)態(tài)的非晶態(tài)合金,因其長(zhǎng)程無(wú)序、短程有序的“玻璃態(tài)”各向同性排列的結(jié)構(gòu)特點(diǎn),而具有不同于傳統(tǒng)金屬材料的優(yōu)異性能[4-7]。例如,CHENG等[8]基于雙金屬間的協(xié)同作用發(fā)現(xiàn),非晶態(tài)Co-Ni體系摻雜非金屬原子的加氫催化活性明顯優(yōu)于單金屬原子與非金屬原子組合成分的活性。文獻(xiàn)[9-11]通過(guò)研究發(fā)現(xiàn),過(guò)渡金屬原子摻雜非金屬原子所形成的非晶態(tài)合金,因其內(nèi)在的協(xié)同效應(yīng),在一定程度上提高了物質(zhì)的催化性能;而非晶態(tài)合金“玻璃態(tài)”各向同性排列的特點(diǎn),也使其表面分布更多的活性位,以達(dá)到吸附更多反應(yīng)物分子的目的,從而增強(qiáng)對(duì)反應(yīng)物的活化能力[12]。查閱大量文獻(xiàn)后發(fā)現(xiàn),作為催化劑領(lǐng)域的佼佼者,非晶態(tài)合金體系一直是科學(xué)家們研究的熱點(diǎn)對(duì)象,其中非貴金屬化合物Co-Ni-B三元體系,更是因被證實(shí)在催化方面具有優(yōu)異性能而備受青睞[13-15]。鑒于非晶態(tài)Co-Ni-B三元體系具有優(yōu)異的催化性質(zhì),其常被用于不同類(lèi)型反應(yīng)的催化劑,這也就意味著,非晶態(tài)合金Co-Ni-B三元體系具有較強(qiáng)的研究意義和廣闊的發(fā)展前景。但截至目前,對(duì)于該體系的研究大多停留在宏觀實(shí)驗(yàn)角度,微觀理論方面的研究鮮少,因此,本文以已構(gòu)建出的Co3NiB2團(tuán)簇模型為基礎(chǔ)[16],對(duì)其極化率、偶極矩及態(tài)密度方面進(jìn)行微觀理論的探究,以期能為今后有關(guān)非晶態(tài)合金Co-Ni-B三元體系的研究提供有價(jià)值的理論參考。

    1 計(jì)算方法

    基于拓?fù)鋵W(xué)原理[17]并依托于HP-Z440計(jì)算機(jī)上的Gaussian09程序,在B3LYP/Lan12dz水平下,對(duì)團(tuán)簇Co3NiB2所設(shè)計(jì)出的全部初始構(gòu)型,采用密度泛函理論[18-20],在自旋多重度為二、四的情況下進(jìn)行全參數(shù)優(yōu)化和相關(guān)頻率計(jì)算。對(duì)計(jì)算所得12種優(yōu)化構(gòu)型的過(guò)渡金屬Co、Ni原子采用18-eECP雙ξ基組(3s,3p,3d/2s,2p,2d),對(duì)類(lèi)金屬B原子采用Dunning/Huzinaga雙ξ基組(9s,5p/3s,2p)[21-24],其中與態(tài)密度相關(guān)的各項(xiàng)參數(shù)均使用Multiwfn程序[25]完成。

    2 實(shí)驗(yàn)結(jié)果與數(shù)據(jù)分析

    2.1 團(tuán)簇Co3NiB2的優(yōu)化構(gòu)型

    依據(jù)上述計(jì)算方法,對(duì)優(yōu)化所得的穩(wěn)定構(gòu)型進(jìn)行整理和篩選之后,可將團(tuán)簇Co3NiB2中的相同構(gòu)型和含虛頻的不穩(wěn)定構(gòu)型排除,最終得到如圖1所示的12種優(yōu)化構(gòu)型,其分別以五棱錐、單帽四棱錐、單帽三角雙錐、四角雙錐4種空間結(jié)構(gòu)類(lèi)型存在。

    為方便進(jìn)一步探究,設(shè)能量最低的五棱錐構(gòu)型1(4)的能量為0 kJ/mol,將其余優(yōu)化構(gòu)型的相對(duì)能量由低到高進(jìn)行排序。其中,能量越低,穩(wěn)定性越強(qiáng)。結(jié)合構(gòu)型能量與空間結(jié)構(gòu)分析可得,團(tuán)簇Co3NiB2 12種優(yōu)化構(gòu)型的穩(wěn)定性由高到低依次為:五棱錐構(gòu)型>單帽四棱錐構(gòu)型>單帽三角雙錐構(gòu)型>四角雙錐構(gòu)型。

    22 團(tuán)簇Co3NiB2的極化率分析

    電子云是表征原子或分子中電子出現(xiàn)的概率密度。由于電子易受到外電場(chǎng)中磁場(chǎng)的作用,因此電子云形狀也會(huì)隨著外場(chǎng)作用的變化而發(fā)生改變。極化率作為衡量電子云形狀改變程度的物理量,不僅在一定程度上反映了極化過(guò)程的難易程度,而且因其對(duì)團(tuán)簇結(jié)構(gòu)極為敏感的特性可用于分析團(tuán)簇的變形程度。與此同時(shí),極化率還可表征物質(zhì)對(duì)外界電場(chǎng)的響應(yīng),是描述物質(zhì)與光非線性作用的物理量之一,在一定程度上反映出了分子間碰撞過(guò)程中的色散力、誘導(dǎo)力和取向作用力等力之間的聯(lián)系和作用,是分析物質(zhì)變形性的常用手段。

    極化率張量平均值是極化率在xx、yy和zz軸上張量的平均值。極化率張量平均值(〈α〉)越小,說(shuō)明團(tuán)簇結(jié)構(gòu)結(jié)合得越緊密,團(tuán)簇原子間的相互作用就越強(qiáng),其在外場(chǎng)作用下就越易保持原有的電子云形狀,團(tuán)簇形變程度就越小;反之,極化率張量平均值越大,團(tuán)簇形變程度就越大。極化率的各向異性不變量(Δα)可表征團(tuán)簇在外加電場(chǎng)中的響應(yīng)程度,是反映在外場(chǎng)作用下團(tuán)簇極化率各向異性的物理量。Δα值越大,團(tuán)簇構(gòu)型就越易受到外場(chǎng)影響,構(gòu)型極化率的各向異性就越大;反之,Δα值越小,團(tuán)簇構(gòu)型就越難受到外場(chǎng)影響,構(gòu)型極化率的各向異性就越小。團(tuán)簇Co3NiB2各優(yōu)化構(gòu)型的極化率張量平均值與極化率的各向異性不變量,通過(guò)以下公式進(jìn)行計(jì)算,計(jì)算結(jié)果如表1所示。

    〈α〉=(αxx+αyy+αzz)/3

    △α=(αxx-αyy)2+(αyy-αzz)2+(αzz-αxx)2+6(α2xy+α2xz+2yz)2

    由表1數(shù)據(jù)可知:該團(tuán)簇體系的極化率張量主要集中在xx(173625~295870 a.u.)、yy(158509~181113 a.u.)及zz(121631~166432 a.u.)3條主軸方向上,且主軸方向上的極化率張量數(shù)值遠(yuǎn)遠(yuǎn)大于xy(-10240~10380 a.u.)、xz(-9087~7585 a.u.)、yz(-3306~4886 a.u.)3條非主軸方向上的極化率張量。其中,極化率張量最大xx分量為構(gòu)型2(4)(295870 a.u.),最大yy分量為構(gòu)型7(4)(181113 a.u.),最大zz分量為構(gòu)型5(4)(166432 a.u.)。

    分析團(tuán)簇Co3NiB2各優(yōu)化構(gòu)型的極化率張量平均值可知:構(gòu)型4(4)、7(4)的極化率張量平均值最小且相近,說(shuō)明這兩者的團(tuán)簇結(jié)構(gòu)結(jié)合得最緊密,團(tuán)簇原子間的相互作用最強(qiáng),其在外場(chǎng)作用下最易保持原有的電子云形狀,團(tuán)簇形變程度最小;而構(gòu)型1(4)、1(2)、2(4)則與之相反,它們的極化率張量平均值最大且相近,說(shuō)明這三者在外場(chǎng)作用下最難保持原有的電子云形狀,團(tuán)簇形變程度最大。由此,根據(jù)表1數(shù)據(jù)可判斷,團(tuán)簇Co3NiB2的形變程度由大到小依次為:1(4)>1(2)>2(4)>5(4)>3(2)>2(2)>3(4)>5(2)>4(2)>6(4)>4(4)>7(4)。結(jié)合圖1各優(yōu)化構(gòu)型的空間結(jié)構(gòu)分析可得,該團(tuán)簇的形變程度由大到小依次為:五棱錐構(gòu)型>單帽四棱錐構(gòu)型>單帽三角雙錐構(gòu)型>四角雙錐構(gòu)型,這與21節(jié)中團(tuán)簇Co3NiB2穩(wěn)定性排序一致,說(shuō)明極化率張量對(duì)團(tuán)簇的幾何結(jié)構(gòu)具有很強(qiáng)的依賴性。

    分析團(tuán)簇Co3NiB2各優(yōu)化構(gòu)型的極化率各向異性不變量可知:構(gòu)型4(4)、6(4)的極化率各向異性不變量最小且相近(均小于15 au),說(shuō)明這兩者對(duì)外場(chǎng)的響應(yīng)最弱,最不易受到外場(chǎng)影響;而構(gòu)型2(4)、1(2)、1(4)的極化率各向異性不變量最大且均大于130 au,說(shuō)明這三者對(duì)外場(chǎng)的響應(yīng)最強(qiáng),最易受到外場(chǎng)的影響。此外,根據(jù)表1數(shù)據(jù)還可判斷,團(tuán)簇Co3NiB2對(duì)外場(chǎng)的響應(yīng)程度由強(qiáng)到弱依次為:2(4)>1(2)>1(4)>4(2)>2(2)>3(4)>3(2)>5(4)>7(4)>5(2)>6(4)>4(4)。同樣,結(jié)合圖1各優(yōu)化構(gòu)型的空間結(jié)構(gòu)分析可得,該團(tuán)簇受外場(chǎng)影響的難易程度由強(qiáng)到弱依次為:?jiǎn)蚊彼睦忮F構(gòu)型>五棱錐構(gòu)型>單帽三角雙錐構(gòu)型>四角雙錐構(gòu)型,這與21節(jié)中團(tuán)簇Co3NiB2穩(wěn)定性排序相比,僅五棱錐構(gòu)型與單帽四棱錐構(gòu)型的順序進(jìn)行了互換,且二者空間構(gòu)型的能量值相差不大,說(shuō)明極化率各向異性不變量對(duì)團(tuán)簇的幾何結(jié)構(gòu)具有較強(qiáng)的依賴性,但依賴性不及極化率張量。

    23 團(tuán)簇Co3NiB2的偶極矩分析

    由一定種類(lèi)和數(shù)目的原子通過(guò)錯(cuò)綜復(fù)雜的空間排列可以形成團(tuán)簇,若團(tuán)簇中某個(gè)原子的位置發(fā)生變化,那么其所對(duì)應(yīng)的電子分布也會(huì)隨之改變。偶極矩可以表征物質(zhì)分子中正負(fù)電荷的分布情況及運(yùn)行規(guī)律,也是判斷物質(zhì)是否存在極性及量度團(tuán)簇極性大小的常用方法之一。團(tuán)簇的偶極矩與其極性密切相關(guān),偶極矩越大,構(gòu)型的極性越強(qiáng)。偶極矩和X、Y、Z軸3個(gè)方向上的偶極矩分量關(guān)系式為μ=(μ2X+μ2Y+μ2Z),式中:μX、μY、μZ分別為X、Y、Z軸3個(gè)方向上的偶極矩分量;μ為物質(zhì)的總偶極矩。表2即為團(tuán)簇Co3NiB2各優(yōu)化構(gòu)型的偶極矩分量及總偶極矩。

    分析表2數(shù)據(jù)可知:團(tuán)簇Co3NiB2的12個(gè)優(yōu)化構(gòu)型在X、Y、Z方向上的偶極矩均不為0,說(shuō)明團(tuán)簇Co3NiB2所有優(yōu)化構(gòu)型均為極性分子。其中,構(gòu)型4(4)的總偶極矩最大(3549 3 D),說(shuō)明該構(gòu)型的分子極性最強(qiáng);而與之相反的構(gòu)型4(2)總偶極矩最小(0356 9 D),表明該構(gòu)型在團(tuán)簇所有的優(yōu)化構(gòu)型中極性最弱。不難發(fā)現(xiàn),團(tuán)簇Co3NiB2所有優(yōu)化構(gòu)型的極性大小排序?yàn)椋?(4)>3(2)>3(4)>2(4)>2(2)>5(4)>1(2)>1(4)>6(4)>5(2)>7(4)>4(2)。整體來(lái)看,X、Y、Z軸三個(gè)方向上的偶極矩?cái)?shù)據(jù)絕對(duì)值的平均值滿足Y軸(1507 5)>X軸(1234 3)>Z軸(1041 9),說(shuō)明在Y軸方向上正、負(fù)電中心偏移距離最大,Z軸正、負(fù)電中心偏移距離最小,這也就意味著Y軸對(duì)團(tuán)簇Co3NiB2偶極矩的貢獻(xiàn)最大,Z軸對(duì)團(tuán)簇Co3NiB2偶極矩的貢獻(xiàn)最小。

    此外,偶極矩屬于分子的靜態(tài)性質(zhì),即在分子的點(diǎn)群對(duì)稱下,偶極矩的大小和方向均不變,這也就意味著團(tuán)簇的結(jié)構(gòu)對(duì)稱性也可通過(guò)偶極矩間接反映。將Gaussian 09程序計(jì)算所得的out文件輸入Gauss View軟件后可知,團(tuán)簇Co3NiB2的12種優(yōu)化構(gòu)型的點(diǎn)群均為C1,這表明該團(tuán)簇所有優(yōu)化構(gòu)型的偶極矩都位于各自構(gòu)型的一重對(duì)稱軸上。

    24 團(tuán)簇Co3NiB2軌道雜化分析

    態(tài)密度(density of state,DOS)是能量介于E與E+dE之間的電子狀態(tài)數(shù)目,一般用E表示,可用于表征電子結(jié)構(gòu)的微觀分布。分析態(tài)密度可以得到團(tuán)簇的成鍵性質(zhì)及團(tuán)簇軌道的具體雜化情況。其中,總態(tài)密度(total density of state,TDOS)可用于表征各優(yōu)化構(gòu)型的總體電子分布情況,分波態(tài)密度(partial-wave density of state,PDOS)則可表征各優(yōu)化構(gòu)型單個(gè)原子中不同軌道的電子狀態(tài)及貢獻(xiàn)電子的能力。因此,通過(guò)分析圖2中的總態(tài)密度圖與分波態(tài)密度圖可以得到團(tuán)簇Co3NiB2軌道的雜化情況。圖中的豎直虛線EF表示Fermi費(fèi)米能級(jí)。

    觀察圖2中各優(yōu)化構(gòu)型的總態(tài)密度圖發(fā)現(xiàn),團(tuán)簇Co3NiB2所有優(yōu)化構(gòu)型均具有A、B、C、D、E 5個(gè)主要高峰,其中A、B均位于Fermi能級(jí)左側(cè),C、D、E均位于Fermi能級(jí)右側(cè)。對(duì)各個(gè)峰進(jìn)行具體分析發(fā)現(xiàn):所有優(yōu)化構(gòu)型的A峰均由B原子的2s軌道貢獻(xiàn);B峰均由Co-3d、Ni-3d、B-2p軌道貢獻(xiàn),說(shuō)明該峰存在d-d-p雜化,其中Co-3d、Ni-3d二者起主要貢獻(xiàn)作用;D峰均由Co-4p與Ni-4p軌道貢獻(xiàn),說(shuō)明所有優(yōu)化構(gòu)型的D峰均存在p-p雜化;而E峰則是均由Co-4s、Co-4p、Ni-4s、Ni-4p以及B-2p軌道貢獻(xiàn),說(shuō)明該峰處存在s-s-p-p-p雜化,其中Co-4s、Ni-4s及B-2p起主要貢獻(xiàn)作用。再對(duì)唯一不同的C峰進(jìn)行分析發(fā)現(xiàn):在四重態(tài)中,構(gòu)型1(4)—4(4)均由Co-3d、Co-4p、Ni-4p及B-2p軌道貢獻(xiàn),說(shuō)明它們存在d-p-p-p雜化,而構(gòu)型5(4)—7(4)均由Co-4p、Ni-3d、Ni-4p及B-2p軌道貢獻(xiàn),說(shuō)明這3種構(gòu)型也存在d-p-p-p雜化,這也就意味著所有四重態(tài)構(gòu)型的雜化方式均為d-p-p-p雜化;在二重態(tài)中,除構(gòu)型5(2)與四重態(tài)構(gòu)型1(4)—4(4)的雜化方式相同外,其余二重態(tài)構(gòu)型1(2)—4(2)的雜化方式均一致,即所有二重態(tài)構(gòu)型(除構(gòu)型5(2)外)均是由Co-3d、Co-4p、Ni-3d、Ni-4p及B-2p軌道主要貢獻(xiàn),存在d-d-p-p-p雜化。

    綜上分析可知:團(tuán)簇Co3NiB2的成鍵主要存在5種雜化方式,即p-p、d-d-p、d-p-p-p、s-s-p-p-p和d-d-p-p-p雜化,其中,p-p、d-d-p雜化作用較強(qiáng),d-p-p-p、s-s-p-p-p、d-d-p-p-p雜化作用較弱。因此,可推斷出p-p、d-d-p較強(qiáng)的雜化作用在促進(jìn)團(tuán)簇Co3NiB2成鍵形成的同時(shí)也增強(qiáng)了各優(yōu)化構(gòu)型的穩(wěn)定性。與此同時(shí),所有優(yōu)化構(gòu)型的A、B、D、E峰處的雜化情況完全一致,僅有C峰處的雜化情況因構(gòu)型重態(tài)的不同而存在差異。其中,所有四重態(tài)構(gòu)型的C峰處均是d-p-p-p雜化,所有二重態(tài)構(gòu)型(除構(gòu)型5(2)外)的C峰處均是d-d-p-p-p雜化,這說(shuō)明C峰處的雜化情況會(huì)在一定程度上受團(tuán)簇多重度的影響。

    3 結(jié)論

    本文以已構(gòu)建出的Co3NiB2團(tuán)簇模型為基礎(chǔ),通過(guò)結(jié)合構(gòu)型能量與空間結(jié)構(gòu)所得的團(tuán)簇穩(wěn)定性對(duì)其極化率、偶極矩及態(tài)密度方面進(jìn)行微觀理論的探究,最終得到以下結(jié)論:

    1)從極化率張量平均值角度分析:團(tuán)簇Co3NiB2的形變程度由大到小依次為:1(4)>1(2)>2(4)>5(4)>3(2)>2(2)>3(4)>5(2)>4(2)>6(4)>4(4)>7(4),即五棱錐構(gòu)型>單帽四棱錐構(gòu)型>單帽三角雙錐構(gòu)型>四角雙錐構(gòu)型,這與團(tuán)簇Co3NiB2的穩(wěn)定性排序一致,說(shuō)明極化率對(duì)于團(tuán)簇的幾何結(jié)構(gòu)具有很強(qiáng)的依賴性。

    2)從極化率各向異性不變量角度分析:團(tuán)簇Co3NiB2對(duì)外場(chǎng)的響應(yīng)程度由強(qiáng)到弱依次為:2(4)>1(2)>1(4)>4(2)>2(2)>3(4)>3(2)>5(4)>7(4)>5(2)>6(4)>4(4),即單帽四棱錐構(gòu)型>五棱錐構(gòu)型>單帽三角雙錐構(gòu)型>四角雙錐構(gòu)型,這與21節(jié)中團(tuán)簇Co3NiB2的穩(wěn)定性排序相差不大,說(shuō)明極化率各向異性不變量對(duì)于團(tuán)簇的幾何結(jié)構(gòu)具有較強(qiáng)的依賴性,但依賴性不及極化率張量。

    3)從偶極矩角度分析:團(tuán)簇Co3NiB2所有優(yōu)化構(gòu)型均為極性分子,其中構(gòu)型4(4)的分子極性最強(qiáng),構(gòu)型4(2)的分子極性最弱,且團(tuán)簇Co3NiB2優(yōu)化構(gòu)型的極性大小排序?yàn)椋?(4)>3(2)>3(4)>2(4)>2(2)>5(4)>1(2)>1(4)>6(4)>5(2)>7(4)>4(2)。其中,Y軸對(duì)團(tuán)簇偶極矩的貢獻(xiàn)最大,Z軸對(duì)團(tuán)簇偶極矩的貢獻(xiàn)最小。

    4)從態(tài)密度角度分析:團(tuán)簇Co3NiB2的成鍵主要存在p-p、d-d-p、d-p-p-p、s-s-p-p-p和d-d-p-p-p 5種雜化方式,其中,p-p、d-d-p雜化作用較強(qiáng),d-p-p-p、s-s-p-p-p、d-d-p-p-p雜化作用較弱。與此同時(shí),所有優(yōu)化構(gòu)型的A、B、D、E峰處的雜化情況完全一致,僅有C峰處的雜化情況因構(gòu)型重態(tài)的不同而存在差異。其中,所有四重態(tài)構(gòu)型的C峰處均是d-p-p-p雜化,所有二重態(tài)構(gòu)型(除構(gòu)型5(2)外)的C峰處均是d-d-p-p-p雜化,這說(shuō)明C峰處的雜化情況會(huì)在一定程度上受團(tuán)簇多重度的影響。

    參考文獻(xiàn):

    [1]YIN X L, WANG Q, DUAN D H, et al. Amorphous NiB alloy decorated by Cu as the anode catalyst for a direct borohydride fuel cell[J]. International Journal of Hydrogen Energy, 2019, 44(21): 10971-10981.

    [2] LI J Z, WEI G D, ZHU Y K, et al. Hierarchical NiCoP nanocone arrays supported on Ni foam as an efficient and stable bifunctional electrocatalyst for overall water splitting[J]. Journal of Materials Chemistry A, 2017, 5(28): 14828-14837.

    [3] HU E L, FENG Y F, NAI J W, et al. Construction of hierarchical Ni-Co-P hollow nanobricks with oriented nanosheets for efficient overall water splitting[J]. Energy & Environmental Science, 2018, 11(4): 872-880.

    [4] JIN H Y, LIU X, CHEN S M, et al. Heteroatom-doped transition metal electrocatalysts for hydrogen evolution reaction[J]. ACS Energy Letters, 2019, 4(4): 805-810.

    [5] LU W T, LI X W, WEI F, et al. In-situ transformed Ni, S-codoped CoO from amorphous Co-Ni sulfide as efficient electrocatalyst for hyd rogen evolution in alkaline media[J]. ACS Sustainable Chemistry & Engineering, 2019, 7(14): 12501-12509.

    [6] WU Y Y, TARIQ M, ZAMAN W Q, et al. Ni-Co codoped RuO2 with outstanding oxygen evolution reaction performance[J]. ACS Applied Energy Materials, 2019, 2(6):4105-4110.

    [7] PIAO J Y, LIU X C, WU J, et al. Construction of uniform Cobalt-based nanoshells and its potential for improving Li-ion battery performance[J]. ACS Applied Materials & Interfaces, 2018, 10(27): 22896-22901.

    [8] CHENG Q Y, LIU D J, WANG M M, et al. Study on catalytic performance of Ni-Co-P amorphous alloy for HDO of vanillin[J]. Journal of Fuel Chemistry Technology, 2019, 47(10): 1205-1213.

    [9] ZHOU L M,ZHANG T R,TAO Z L. Ni nano particles supported on carbon as efficient catalysts for the hydrolysis of ammonia borane[J]. Nano Research, 2014, 7(5): 774 -781.

    [10]FIGEN K A, PISKIN M B, COSKUNER B, et al. Synthesis, structural characterization, and hydrolysis of ammonia borane (NH3BH3) as a hydrogen storage carrier[J]. International Journal of Hydrogen Energy, 2013, 38(36): 16215-16228.

    [11]WANG W J, LI H, LI H X, et al. Regeneration of the amorphous NiB/SiO2 catalyst poisoned by carbon disulfide in cyclopentadiene hydrogenation[J]. Applied Catalysis A: General, 2000, 203(2): 301-306.

    [12]YUAN Z Z,ZHANG D P. Preparation and magnetic properties of amorphous Co-Cu-B alloy nano-powders[J]. Intermetallics, 2009, 17(4): 281-284.

    [13]FERNANDES R, PATEL N, MIOTELLO A, et al. Studies on catalytic behavior of Co-Ni-B in hydrogen production by hydrolysis of NaBH4[J]. Journal of Molecular Catalysis A: Chemical, 2009, 298(1/2): 1-6.

    [14]GUO H J, ZHANG H R, ZHANG L Q, et al. Selective hydrogenation of furfural to furfuryl alcohol over acid-activated attapulgite-supported NiCoB amorphous alloy catalyst[J]. Industrial & Engineering Chemistry Research, 2018, 57(2): 498-551.

    [15]ZHANG J, LI X X, LIU Y T, et al. Bi-metallic boride electrocatalysts with enhanced activity for the oxygen evolution reaction[J]. Nanoscale, 2018, 10(25): 11997-12002.

    [16]秦渝, 方志剛, 張偉, 等. 團(tuán)簇Co3NiB催化析氫活性研究[J]. 江西師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2020, 44(1): 56-62.

    [17]YU X G. Hyperbolic multi-topology and the basic principle in quantum mechanics[J]. Advances in Applied Clifford Algebras, 1999, 9(1): 109-118.

    [18]譚興鳳, 蔡雪梅, 隆正文. 電場(chǎng)作用下MgS分子電子激發(fā)特性的理論研究[J]. 貴州大學(xué)學(xué)報(bào)(自然科學(xué)版), 2011, 28(3): 15-20.

    [19]HAY P J. Gaussian basis sets for molecular calculations. The representation of 3D orbitals in transition-metal atoms[J]. The Journal of Chemical Physics, 1977, 66(10): 4377-4384.

    [20]李歷紅, 方志剛, 趙振寧, 等. 團(tuán)簇Ni3CoP電子性質(zhì)與磁性研究[J]. 江西師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2019, 43(2): 54-60.

    [21]秦渝, 方志剛, 趙璐璐, 等.團(tuán)簇Co3NiB2異構(gòu)化反應(yīng)的動(dòng)力學(xué)與熱力學(xué)研究[J]. 江西師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2021, 45(1): 67-74.

    [22]徐詩(shī)浩, 方志剛, 韓建銘, 等. 團(tuán)簇V3B2成鍵及磁學(xué)性質(zhì)研究[J]. 廣西師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2017, 35(3): 89-96.

    [23]李雯博, 方志剛, 趙振寧, 等. 團(tuán)簇Co5B2反應(yīng)活性的DFT研究[J]. 廣西師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2017, 35(4): 76-83.

    [24]HAY P J, WADT W R. Ab initio effective core potentials for molecular calculations. Potentials for the transition metal atoms Sc to Hg[J]. The Journal of Chemical Physics, 1985, 82(1): 270-283.

    [25]LU T, CHEN F W. Multiwfn: a multifunctional wavefunction analyzer[J]. Journal of Computational Chemistry, 2012, 33(5): 580-592.

    (責(zé)任編輯:周曉南)

    Study on the Dolarizability, Dipole Moment and

    Density of States of Cluster Co3NiB2

    FANG Zhigang*, WANG Zhiyao, ZHENG Xinxi, QIN Yu, MAO Zhilong,

    ZENG Xinyu, ZHU Yiwen, WANG Qian

    (School of Chemical Engineering,University of Science and Technology Liaoning,Anshan 114051,China)

    Abstract:

    In order to explore the internal structure stability and polarity of cluster Co3NiB2, based on density functional theory and topology theory, the polarizability, dipole moment and density of states of 12 optimized configurations of cluster Co3NiB2 were studied at B3LYP/Lanl2dz level, and the following conclusions were obtained. The order analysis of cluster deformation shows that the polarizability is strongly dependent on the cluster geometry. By analyzing the order of the cluster response to the external field, it is found that the polarizability anisotropy invariant has a strong dependence on the cluster geometry, but the dependence is less than that of polarizability. From the perspective of dipole moment, it is found that all the optimized configurations of the cluster are polar molecules, in which the polarity of configuration 4(4) is the strongest and that of configuration 4 (2) is the weakest. By analyzing the density of states, it is found that there are mainly five hybrid modes of p-p, d-d-p, d-p-p, s-s-p-p-p and d-d-p-p-p hybrid modes in the cluster bonding, among which p-p and d-d-p-p hybrid modes are stronger than d-p-p, s-s-p-p-p and d-d-p-p hybrid modes. In addition, the hybridization at A, B, D and E peaks of all optimized configurations of the cluster are identical, and only the hybridization at C peak is different due to the different configuration heavy states, which indicates that the hybridization at C peak is affected by the cluster multiplicity to a certain extent.

    Key words:

    cluster Co3NiB2; polarizability; dipole moment; density of states

    91久久精品国产一区二区三区| 黄色怎么调成土黄色| 禁无遮挡网站| 久久久久久久久久久丰满| 国产精品99久久久久久久久| 只有这里有精品99| 亚洲内射少妇av| 国产免费一区二区三区四区乱码| 九九爱精品视频在线观看| 99re6热这里在线精品视频| 国产极品天堂在线| 狂野欧美激情性xxxx在线观看| 麻豆成人午夜福利视频| 国产淫语在线视频| 午夜福利高清视频| 一级毛片电影观看| 亚洲成人久久爱视频| 亚洲在线观看片| 大陆偷拍与自拍| 日韩伦理黄色片| 午夜亚洲福利在线播放| 久久99热6这里只有精品| 成人午夜精彩视频在线观看| 国产亚洲5aaaaa淫片| 国产亚洲午夜精品一区二区久久 | 在线播放无遮挡| 国产爱豆传媒在线观看| 日韩av免费高清视频| 春色校园在线视频观看| 亚洲精品乱码久久久v下载方式| 精品一区二区三卡| 欧美xxxx黑人xx丫x性爽| 国产黄色免费在线视频| 亚洲成人精品中文字幕电影| 一级片'在线观看视频| 深爱激情五月婷婷| 2021少妇久久久久久久久久久| 超碰97精品在线观看| 涩涩av久久男人的天堂| 亚洲三级黄色毛片| 18+在线观看网站| 久久精品国产亚洲av涩爱| 老司机影院毛片| 国产有黄有色有爽视频| 国产女主播在线喷水免费视频网站| 午夜福利高清视频| 欧美日韩精品成人综合77777| 丰满人妻一区二区三区视频av| 国产亚洲午夜精品一区二区久久 | 欧美xxxx性猛交bbbb| 国产久久久一区二区三区| 久久国内精品自在自线图片| 国产亚洲av嫩草精品影院| 国产精品国产三级国产专区5o| 免费黄网站久久成人精品| 亚洲av在线观看美女高潮| 午夜精品国产一区二区电影 | 欧美日本视频| 少妇高潮的动态图| 国产真实伦视频高清在线观看| 久久久久久久久大av| 七月丁香在线播放| 亚洲欧洲日产国产| 有码 亚洲区| 色播亚洲综合网| 亚洲婷婷狠狠爱综合网| 亚洲最大成人手机在线| 麻豆乱淫一区二区| 国产综合懂色| 免费av观看视频| 99久久九九国产精品国产免费| 欧美精品人与动牲交sv欧美| 天天躁夜夜躁狠狠久久av| 国产精品无大码| 五月伊人婷婷丁香| 国产伦精品一区二区三区视频9| 免费不卡的大黄色大毛片视频在线观看| 国产亚洲91精品色在线| 少妇裸体淫交视频免费看高清| 一个人观看的视频www高清免费观看| 嫩草影院入口| 一级黄片播放器| 内地一区二区视频在线| 涩涩av久久男人的天堂| 国产午夜精品一二区理论片| 老司机影院毛片| 蜜臀久久99精品久久宅男| 国产精品无大码| 久久国产乱子免费精品| 高清日韩中文字幕在线| 激情五月婷婷亚洲| 久久6这里有精品| 高清欧美精品videossex| 男人爽女人下面视频在线观看| 97精品久久久久久久久久精品| a级毛色黄片| 中文字幕免费在线视频6| 免费少妇av软件| 欧美+日韩+精品| 成人亚洲欧美一区二区av| 精品人妻偷拍中文字幕| av国产免费在线观看| 久久久久久国产a免费观看| 亚洲婷婷狠狠爱综合网| 久久久久国产精品人妻一区二区| 别揉我奶头 嗯啊视频| 小蜜桃在线观看免费完整版高清| 99精国产麻豆久久婷婷| 制服丝袜香蕉在线| 一级片'在线观看视频| 亚洲av成人精品一二三区| 我要看日韩黄色一级片| 丰满少妇做爰视频| 欧美丝袜亚洲另类| 超碰97精品在线观看| 午夜视频国产福利| 一级毛片黄色毛片免费观看视频| 麻豆久久精品国产亚洲av| 久久久久久久久久久免费av| 黄色日韩在线| 免费观看在线日韩| 日本-黄色视频高清免费观看| 亚洲色图av天堂| 成年版毛片免费区| 亚洲欧美精品自产自拍| 男女啪啪激烈高潮av片| 日韩伦理黄色片| a级毛片免费高清观看在线播放| 一本久久精品| 特大巨黑吊av在线直播| 精品视频人人做人人爽| 插逼视频在线观看| 亚洲欧洲国产日韩| 欧美性猛交╳xxx乱大交人| 熟女人妻精品中文字幕| 蜜桃久久精品国产亚洲av| 久久久精品免费免费高清| 欧美xxⅹ黑人| 国产成人免费观看mmmm| 欧美zozozo另类| 在线免费观看不下载黄p国产| 亚洲天堂国产精品一区在线| 中国美白少妇内射xxxbb| 夜夜看夜夜爽夜夜摸| 男人舔奶头视频| 成人黄色视频免费在线看| 久久久久久伊人网av| 2018国产大陆天天弄谢| 亚洲国产精品成人综合色| 久久久色成人| 啦啦啦啦在线视频资源| 成年av动漫网址| 一级爰片在线观看| 成人高潮视频无遮挡免费网站| 中文字幕av成人在线电影| 七月丁香在线播放| 亚洲av电影在线观看一区二区三区 | 日日啪夜夜爽| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品乱久久久久久| 国产视频首页在线观看| 嫩草影院入口| 天堂俺去俺来也www色官网| 亚洲精品一区蜜桃| 国产成人一区二区在线| 国产有黄有色有爽视频| 日本色播在线视频| 成人特级av手机在线观看| 毛片女人毛片| 尤物成人国产欧美一区二区三区| 国产精品不卡视频一区二区| 老司机影院毛片| 日韩精品有码人妻一区| 成人漫画全彩无遮挡| 一级二级三级毛片免费看| 国产精品嫩草影院av在线观看| 国产日韩欧美在线精品| 精品少妇黑人巨大在线播放| 久久久亚洲精品成人影院| 亚洲自拍偷在线| 色哟哟·www| 王馨瑶露胸无遮挡在线观看| 干丝袜人妻中文字幕| 91精品一卡2卡3卡4卡| 在线观看av片永久免费下载| av又黄又爽大尺度在线免费看| 午夜老司机福利剧场| 午夜精品国产一区二区电影 | 免费观看av网站的网址| 午夜日本视频在线| 免费电影在线观看免费观看| 王馨瑶露胸无遮挡在线观看| 欧美老熟妇乱子伦牲交| 亚洲国产精品专区欧美| 婷婷色av中文字幕| 九九在线视频观看精品| 亚洲精品久久午夜乱码| 国产色婷婷99| 九九久久精品国产亚洲av麻豆| 下体分泌物呈黄色| 国产老妇女一区| 中国国产av一级| 中文字幕人妻熟人妻熟丝袜美| 久久久亚洲精品成人影院| 91精品伊人久久大香线蕉| 美女被艹到高潮喷水动态| 精品久久久久久久久亚洲| 久久久色成人| 色网站视频免费| 国产视频首页在线观看| 国产色婷婷99| 婷婷色综合大香蕉| 日韩 亚洲 欧美在线| 欧美老熟妇乱子伦牲交| 蜜臀久久99精品久久宅男| 黑人高潮一二区| 久久久久久久精品精品| 麻豆精品久久久久久蜜桃| 美女被艹到高潮喷水动态| 人妻夜夜爽99麻豆av| 久久精品国产自在天天线| 噜噜噜噜噜久久久久久91| 日韩亚洲欧美综合| 尤物成人国产欧美一区二区三区| 大陆偷拍与自拍| 色吧在线观看| 最近手机中文字幕大全| 国产亚洲av片在线观看秒播厂| 全区人妻精品视频| 国产成人免费观看mmmm| 插阴视频在线观看视频| 亚洲av在线观看美女高潮| 国产精品.久久久| 亚洲av成人精品一二三区| 好男人视频免费观看在线| 伊人久久国产一区二区| 日本与韩国留学比较| 插阴视频在线观看视频| 亚洲色图av天堂| 男女国产视频网站| 免费少妇av软件| 美女主播在线视频| 成年免费大片在线观看| 男的添女的下面高潮视频| 午夜福利视频精品| av在线亚洲专区| 久久久久精品久久久久真实原创| 亚洲成人一二三区av| 国产午夜福利久久久久久| 亚洲精品久久午夜乱码| 日日摸夜夜添夜夜添av毛片| 亚洲内射少妇av| 久久精品国产a三级三级三级| 午夜福利在线在线| 水蜜桃什么品种好| 国产午夜福利久久久久久| 亚洲欧美日韩卡通动漫| 伦精品一区二区三区| 国产国拍精品亚洲av在线观看| 内地一区二区视频在线| 精品久久久久久久久亚洲| 亚洲成人久久爱视频| 99视频精品全部免费 在线| 久久久久久国产a免费观看| 精品人妻一区二区三区麻豆| 黄片wwwwww| 亚洲成人久久爱视频| 国产av码专区亚洲av| 一区二区三区四区激情视频| 尾随美女入室| 又粗又硬又长又爽又黄的视频| 熟妇人妻不卡中文字幕| 亚洲欧美日韩无卡精品| 国产精品国产三级国产专区5o| 寂寞人妻少妇视频99o| 丝袜美腿在线中文| 特级一级黄色大片| 激情 狠狠 欧美| 九九爱精品视频在线观看| 2021天堂中文幕一二区在线观| 日本一二三区视频观看| 插逼视频在线观看| 日日摸夜夜添夜夜爱| 内地一区二区视频在线| 22中文网久久字幕| av女优亚洲男人天堂| 熟女av电影| 欧美日韩精品成人综合77777| 秋霞在线观看毛片| 1000部很黄的大片| 热99国产精品久久久久久7| 久久久久久久久久久免费av| 国产黄片视频在线免费观看| 日韩视频在线欧美| 在线观看三级黄色| av在线老鸭窝| 少妇人妻久久综合中文| 久久久久久久久久久丰满| 亚洲成人中文字幕在线播放| 国产69精品久久久久777片| 在线观看av片永久免费下载| 久久国内精品自在自线图片| 国产黄片美女视频| a级毛片免费高清观看在线播放| 一级爰片在线观看| 亚洲精品色激情综合| 蜜桃久久精品国产亚洲av| 国产 一区精品| 日韩一区二区三区影片| 午夜日本视频在线| 久久人人爽人人片av| 国产免费福利视频在线观看| 美女cb高潮喷水在线观看| 白带黄色成豆腐渣| 女人十人毛片免费观看3o分钟| 国产午夜精品一二区理论片| 激情五月婷婷亚洲| 国产精品国产av在线观看| 高清av免费在线| 又黄又爽又刺激的免费视频.| 九九久久精品国产亚洲av麻豆| 极品教师在线视频| 亚洲天堂av无毛| 麻豆乱淫一区二区| 久久99精品国语久久久| 天堂俺去俺来也www色官网| 男女边摸边吃奶| 久久久久精品性色| 久久热精品热| 少妇人妻精品综合一区二区| 国产高潮美女av| 欧美xxxx黑人xx丫x性爽| 中文资源天堂在线| 久久综合国产亚洲精品| 在线精品无人区一区二区三 | 高清欧美精品videossex| 五月开心婷婷网| 一区二区三区四区激情视频| 大香蕉久久网| 91久久精品国产一区二区成人| 精品少妇黑人巨大在线播放| 91在线精品国自产拍蜜月| 中文资源天堂在线| 黄色怎么调成土黄色| 人人妻人人爽人人添夜夜欢视频 | av播播在线观看一区| 麻豆乱淫一区二区| 亚洲成人精品中文字幕电影| 久热久热在线精品观看| 少妇丰满av| 成人亚洲精品av一区二区| 国产精品国产三级专区第一集| 99热6这里只有精品| 亚洲激情五月婷婷啪啪| 日韩一区二区三区影片| 国产黄色免费在线视频| 建设人人有责人人尽责人人享有的 | 久久久久久久久久久免费av| 嫩草影院新地址| 久久久久久久久久成人| 亚洲,一卡二卡三卡| 日本熟妇午夜| 亚洲国产精品成人综合色| www.av在线官网国产| 嫩草影院精品99| 久久这里有精品视频免费| 精品久久久久久久久av| 久热这里只有精品99| 精品久久久久久久久av| 国产成人午夜福利电影在线观看| 国产成人精品婷婷| 久久这里有精品视频免费| 精品久久国产蜜桃| 国产午夜精品一二区理论片| 国产精品女同一区二区软件| 国产一区二区三区综合在线观看 | 成年av动漫网址| 少妇人妻 视频| 久久精品国产a三级三级三级| 久久99热这里只有精品18| 99久久人妻综合| 精品一区二区免费观看| tube8黄色片| 在线观看三级黄色| 国产日韩欧美亚洲二区| 九九久久精品国产亚洲av麻豆| 久久久久久久久久久丰满| 少妇人妻久久综合中文| 久久久久久久亚洲中文字幕| 麻豆国产97在线/欧美| 精品久久久噜噜| 久久久国产一区二区| 99热这里只有是精品50| 亚洲国产欧美人成| 国产毛片在线视频| 久久精品久久精品一区二区三区| 国产精品一区二区在线观看99| 18禁在线无遮挡免费观看视频| 亚洲成色77777| 最后的刺客免费高清国语| 国产精品久久久久久久电影| 好男人在线观看高清免费视频| 欧美精品一区二区大全| 三级国产精品片| 国产精品成人在线| 国产爽快片一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 新久久久久国产一级毛片| 99久久中文字幕三级久久日本| 久久精品国产鲁丝片午夜精品| 免费电影在线观看免费观看| 亚洲丝袜综合中文字幕| 在线观看国产h片| 伦精品一区二区三区| 一区二区三区四区激情视频| 大陆偷拍与自拍| 热99国产精品久久久久久7| 亚洲av欧美aⅴ国产| 51国产日韩欧美| 2021天堂中文幕一二区在线观| 美女视频免费永久观看网站| 菩萨蛮人人尽说江南好唐韦庄| 国产色爽女视频免费观看| 国产av不卡久久| 人人妻人人看人人澡| 成人漫画全彩无遮挡| 国产成人freesex在线| 国产精品久久久久久av不卡| 亚洲欧美精品自产自拍| 亚洲精品成人av观看孕妇| 欧美日韩在线观看h| 中文字幕免费在线视频6| 久久久久久久久大av| 欧美一级a爱片免费观看看| 久久久久久国产a免费观看| 国产伦理片在线播放av一区| 在线观看免费高清a一片| 国产av码专区亚洲av| 中文乱码字字幕精品一区二区三区| 欧美精品人与动牲交sv欧美| a级毛色黄片| 黄片wwwwww| 中文字幕人妻熟人妻熟丝袜美| 国产亚洲91精品色在线| 中国美白少妇内射xxxbb| 69av精品久久久久久| 在线天堂最新版资源| 国产精品一区www在线观看| 国产精品偷伦视频观看了| 中文字幕制服av| 日韩制服骚丝袜av| 国产淫语在线视频| 观看美女的网站| 精品久久久久久久久av| 99久久精品一区二区三区| 久久影院123| 九草在线视频观看| 色视频在线一区二区三区| 欧美变态另类bdsm刘玥| 亚洲av一区综合| 国产伦理片在线播放av一区| 久久久精品欧美日韩精品| 成人一区二区视频在线观看| 久久人人爽av亚洲精品天堂 | 亚洲欧洲国产日韩| 日韩视频在线欧美| 青青草视频在线视频观看| 亚洲成色77777| 国产老妇伦熟女老妇高清| 国产极品天堂在线| 亚洲美女搞黄在线观看| xxx大片免费视频| 国产精品不卡视频一区二区| 嫩草影院入口| 国产爱豆传媒在线观看| 国产精品蜜桃在线观看| 国产一区亚洲一区在线观看| 国产精品不卡视频一区二区| 亚洲婷婷狠狠爱综合网| 国产永久视频网站| 18禁动态无遮挡网站| 婷婷色麻豆天堂久久| 男人狂女人下面高潮的视频| 国产伦精品一区二区三区视频9| 精品一区二区三卡| 少妇裸体淫交视频免费看高清| 久久久久久久久久人人人人人人| 国产成人aa在线观看| 国产真实伦视频高清在线观看| 尤物成人国产欧美一区二区三区| 一级毛片电影观看| 三级经典国产精品| 黄色视频在线播放观看不卡| 少妇被粗大猛烈的视频| 欧美97在线视频| 美女xxoo啪啪120秒动态图| 免费av不卡在线播放| 国产 精品1| 欧美日韩精品成人综合77777| 国产精品麻豆人妻色哟哟久久| 国产亚洲最大av| 日韩国内少妇激情av| 日本一二三区视频观看| 国产精品偷伦视频观看了| 久久精品国产鲁丝片午夜精品| 色哟哟·www| 九草在线视频观看| av女优亚洲男人天堂| 免费观看a级毛片全部| videos熟女内射| 搡老乐熟女国产| 看十八女毛片水多多多| 水蜜桃什么品种好| 亚洲欧美成人综合另类久久久| 精品一区二区三卡| 国产高清有码在线观看视频| 三级国产精品欧美在线观看| 91精品国产九色| 国产一区亚洲一区在线观看| 特级一级黄色大片| 99久国产av精品国产电影| 亚洲自拍偷在线| 真实男女啪啪啪动态图| 欧美日韩视频高清一区二区三区二| 一级毛片黄色毛片免费观看视频| 日本黄色片子视频| 中国三级夫妇交换| 亚洲无线观看免费| 99久久精品国产国产毛片| 亚洲av中文av极速乱| 丝瓜视频免费看黄片| 夫妻午夜视频| 国产白丝娇喘喷水9色精品| 熟女人妻精品中文字幕| 国产精品av视频在线免费观看| 午夜免费观看性视频| 久久韩国三级中文字幕| 亚洲激情五月婷婷啪啪| 日日啪夜夜爽| 美女xxoo啪啪120秒动态图| 高清午夜精品一区二区三区| 欧美3d第一页| 在线观看三级黄色| 九草在线视频观看| 亚洲欧美一区二区三区国产| 九九久久精品国产亚洲av麻豆| 国产又色又爽无遮挡免| 久久午夜福利片| 汤姆久久久久久久影院中文字幕| 久久精品国产鲁丝片午夜精品| 久久ye,这里只有精品| 精品国产一区二区三区久久久樱花 | 亚洲四区av| 欧美日韩视频精品一区| 日韩电影二区| 天天躁日日操中文字幕| 2022亚洲国产成人精品| 真实男女啪啪啪动态图| 久久久久久久精品精品| 日韩欧美 国产精品| 免费黄网站久久成人精品| 日本猛色少妇xxxxx猛交久久| 成年版毛片免费区| 在线免费观看不下载黄p国产| 国产成人91sexporn| 人妻一区二区av| 毛片女人毛片| 黄片wwwwww| 国产在线男女| 亚洲va在线va天堂va国产| 国产黄频视频在线观看| 午夜精品一区二区三区免费看| 欧美zozozo另类| 亚洲精品乱码久久久久久按摩| 直男gayav资源| 18禁在线播放成人免费| av在线蜜桃| 熟妇人妻不卡中文字幕| 噜噜噜噜噜久久久久久91| 欧美日韩一区二区视频在线观看视频在线 | 久久久欧美国产精品| 亚洲综合精品二区| 欧美日韩在线观看h| 一级毛片aaaaaa免费看小| 好男人视频免费观看在线| 2021天堂中文幕一二区在线观| 国产免费又黄又爽又色| 久久韩国三级中文字幕| 只有这里有精品99| 精品久久久噜噜| 毛片一级片免费看久久久久| 免费看a级黄色片| 天堂网av新在线| av卡一久久| 日本与韩国留学比较| 国产探花在线观看一区二区| 18禁在线无遮挡免费观看视频| 国产精品福利在线免费观看| 精品久久久噜噜| 国产成人a∨麻豆精品| 国产精品偷伦视频观看了| 精品一区二区三卡| 18禁在线无遮挡免费观看视频| 国产午夜精品久久久久久一区二区三区| 男人添女人高潮全过程视频| 久久韩国三级中文字幕| 国产黄频视频在线观看| 精品视频人人做人人爽| 最后的刺客免费高清国语| av免费观看日本| 久久久精品免费免费高清| 亚洲欧美清纯卡通| 欧美bdsm另类| kizo精华| 麻豆成人午夜福利视频|