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

    As摻雜濃度對(duì)FeS2電子結(jié)構(gòu)及光學(xué)性質(zhì)的影響

    2017-04-06 09:16:38馬萬坤陳建華張國(guó)范李玉瓊馮其明
    關(guān)鍵詞:價(jià)帶導(dǎo)帶晶胞

    馬萬坤,陳建華,張國(guó)范,李玉瓊,馮其明

    ?

    As摻雜濃度對(duì)FeS2電子結(jié)構(gòu)及光學(xué)性質(zhì)的影響

    馬萬坤1,陳建華2,張國(guó)范1,李玉瓊2,馮其明1

    (1. 中南大學(xué)資源加工與生物工程學(xué)院,長(zhǎng)沙 410083;2. 廣西大學(xué)資源與冶金學(xué)院,南寧 530004)

    采用基于密度泛函理論框架下的第一性原理平面波超軟贗勢(shì)法,對(duì)不同As摻雜濃度FeS2的幾何結(jié)構(gòu)、電子結(jié)構(gòu)和光學(xué)性質(zhì)進(jìn)行計(jì)算和討論。采用2×2×2(Fe32S63As),2×2×1(Fe16S31As)和2×1×1(Fe8S15As)的超晶胞模型,用1個(gè)As原子取代1個(gè)S原子,使摻雜濃度分別為1.93%、3.82%和7.48%(質(zhì)量分?jǐn)?shù)),然后進(jìn)行3種摻雜體系的計(jì)算,并將其與理想體系進(jìn)行對(duì)比。幾何結(jié)構(gòu)與電子結(jié)構(gòu)的計(jì)算結(jié)果表明:As摻雜使得FeS2的晶格常數(shù)和晶胞體積增大,導(dǎo)帶部分下移,禁帶寬度減小,并且在?10.4~?9.5 eV的淺部?jī)r(jià)帶產(chǎn)生了由As的p態(tài)貢獻(xiàn)的雜質(zhì)能級(jí)。光學(xué)性質(zhì)計(jì)算結(jié)果表明,摻雜后Fe(S1?xAs)2的靜態(tài)介電常數(shù)、折射率和光電導(dǎo)率在一定范圍內(nèi)均隨著摻雜量的增大而明顯增大,說明As摻雜顯著增強(qiáng)了FeS2對(duì)光的吸收以及光電轉(zhuǎn)換效率。

    As摻雜;黃鐵礦;光電性質(zhì);第一性原理

    隨著能源需求的激增和環(huán)境污染的加劇,太陽能電池已成為世界各國(guó)研究的熱點(diǎn)。硅基太陽能電池雖然具有能量轉(zhuǎn)化效率高和穩(wěn)定性好等優(yōu)點(diǎn),但高昂的制造成本限制了其大規(guī)模應(yīng)用[1?3]。因而,繼續(xù)研究并開發(fā)其他新型的太陽能電池材料及先進(jìn)的太陽能電池系統(tǒng),是目前極為重要的課題。金屬硫化物因其半導(dǎo)體屬性在太陽能電池、熱電材料以及存儲(chǔ)裝置等領(lǐng)域顯示出了廣闊的應(yīng)用前景,因而得到了廣泛關(guān)注。通過切實(shí)可行的途徑使金屬硫化物具備優(yōu)良的光學(xué)、電學(xué)或者磁學(xué)性能是其能夠得到有效應(yīng)用的先決條件。

    黃鐵礦(FeS2)是地殼中最常見的硫化礦,組成元素?zé)o毒且地球儲(chǔ)量豐富。立方晶系黃鐵礦具有優(yōu)異的半導(dǎo)體特性:合適的禁帶寬度(g≈0.95 eV)和較高的光吸收系數(shù)(≈5×105/cm,≤700 nm)[4?5],其環(huán)境相容性好,制備成本低廉,是一種較有研究?jī)r(jià)值的新型太陽能電池材料[6?12]。已有研究采用人工合成方法制備出FeS2薄膜,盡管其易于氧化降解,但仍表現(xiàn)出良好的光電性能[13]。1986年,CHATZITHEODOROU等[14]首次制備出FeS2光電化學(xué)電池和肖特基太陽能電池,但是較低的光電壓限制了FeS2在光伏器件方面的應(yīng)用。

    天然黃鐵礦物中存在不同類型的雜質(zhì)和含量,黃鐵礦的光電催化性質(zhì)會(huì)受到這些雜質(zhì)的影響。LEHNER等[15]使用化學(xué)氣相沉積來制備Ni摻雜的黃鐵礦薄膜,證明Ni的摻入在光電導(dǎo)性測(cè)量中會(huì)顯著影響光電導(dǎo)的敏感性。FERRER等[16]制備了摻Cu黃鐵礦薄膜,摻雜后電阻率從0.002 Ω·cm增加到3.28 Ω·cm,摻雜體系光吸收邊從0.8 eV增大到1.1 eV。YU等[17]制備了摻Mn的黃鐵礦薄膜,摻雜后薄膜形貌從顆粒狀變?yōu)槠瑺?,有利于光捕獲。李玉瓊等[18?21]研究了Co、Ni、As等20余種天然雜質(zhì)對(duì)黃鐵礦電子結(jié)構(gòu)和光學(xué)性質(zhì)的影響,結(jié)果表明:Co、Ni、As的摻雜均使得黃鐵礦的吸收帶邊發(fā)生了明顯的紅移。盡管國(guó)內(nèi)外研究人員對(duì)黃鐵礦薄膜的合成制備、As摻雜納米晶薄膜的研究已經(jīng)很多,但對(duì)As摻雜對(duì)黃鐵礦光學(xué)性質(zhì)的影響目前仍鮮見報(bào)道。

    本文作者以FeS2超晶胞作為基體,以As為摻雜原子,采用基于密度泛函理論的第一性原理平面波超軟贗勢(shì)法,對(duì)不同濃度As摻雜FeS2前后的能帶、電子態(tài)密度和光學(xué)性質(zhì)進(jìn)行了對(duì)比和研究,初步揭示了 通過As摻雜FeS2改變其光電性質(zhì)的內(nèi)在原因,為開發(fā)FeS2光伏材料的實(shí)驗(yàn)工作提供了理論參考。

    1 理論模型及計(jì)算方法

    半導(dǎo)體化合物FeS2屬于等軸晶系,其空間群為,每個(gè)晶胞含有4個(gè)FeS2分子單元,F(xiàn)e原子分布在晶胞的6個(gè)面心和8個(gè)頂角上,每個(gè)Fe原子與6個(gè)相鄰的硫原子配位,每個(gè)S原子與3個(gè)Fe原子和1個(gè)S原子配位,2個(gè)硫原子之間形成啞鈴狀結(jié)構(gòu),以硫二聚體形式存在,且沿著(111)方向排列。FeS2的晶格常數(shù)為===0.5379 nm,===90°。

    在本計(jì)算中,首先采用單胞模型描述理想FeS2的電子結(jié)構(gòu)及性質(zhì),然后分別采用2×2×2 (Fe32S63As)、2×2×1(Fe16S31As)和2×1×1(Fe8S15As) 3個(gè)超晶胞模型,用1個(gè)As原子取代1個(gè)S原子,使摻雜濃度分別為1.93%、3.82%和7.48%(質(zhì)量分?jǐn)?shù)),然后進(jìn)行3種摻雜體系的計(jì)算。通過計(jì)算得到3種摻雜濃度下FeS2的優(yōu)化結(jié)構(gòu)、摻雜前后的能帶結(jié)構(gòu)、電子態(tài)密度及光學(xué)性質(zhì)。計(jì)算模型如圖1所示。

    本計(jì)算采用基于密度泛函理論(DFT)框架下的第一性原理贗勢(shì)平面波法,主要的計(jì)算工作由CASTEP軟件包[22]完成。計(jì)算中采用超軟贗勢(shì)(USPP)[23]來處理離子實(shí)與電子之間的相互作用,交換關(guān)聯(lián)泛函采用廣義梯度近似(GGA)下的PW91梯度修正近似。平面波截?cái)嗄軠y(cè)試表明,截?cái)嗄苋?60 eV較為合理。計(jì)算選取的價(jià)電子為:Fe 3d64s2、S 3s23p4和As 4s24p3。采用BFGS算法對(duì)體系進(jìn)行幾何優(yōu)化,得到穩(wěn)定的結(jié)構(gòu),優(yōu)化的收斂標(biāo)準(zhǔn)為:原子位移的收斂閾值為0.0002 nm,原子間作用力的收斂閾值為0.8 eV/nm,原子間的內(nèi)應(yīng)力收斂閾值為0.1 GPa,最大能量改變的收斂閾值為2.0×10?5eV/atom,自洽計(jì)算收斂精度設(shè)置為每個(gè)原子1 meV,布里淵區(qū)的積分分別采用3×5×3、3×3×2和2×2×1的MONKHORST和PACK[24]特殊點(diǎn)對(duì)全Brillouin區(qū)求和。

    2 結(jié)果與討論

    2.1 幾何結(jié)構(gòu)優(yōu)化結(jié)果

    未摻雜的FeS2單胞優(yōu)化后及As摻雜FeS2后的超晶胞模型經(jīng)幾何優(yōu)化后折合的晶胞參數(shù)、總能量的晶格常數(shù)值如表1所列。從表1可以看出,優(yōu)化后FeS2的晶格參數(shù)與實(shí)驗(yàn)值(5.416 nm)[25]非常接近,平衡晶格常數(shù)比實(shí)驗(yàn)值小0.68%。隨著As摻雜濃度增大,F(xiàn)e(S1?xAs)2(=0.0156, 0.0312, 0.0625)的晶格常數(shù)和晶格體積均逐漸增大。引起這種變化的主要原因是As的共價(jià)半徑為0.120 nm,而S的共價(jià)半徑為0.102 nm,As的共價(jià)半徑比S的大,當(dāng)用As來置換S時(shí),就會(huì)引起晶胞體積增大。晶胞優(yōu)化的結(jié)果與文獻(xiàn)[26]的結(jié)果一致。

    圖1 FeS2單胞及As摻雜FeS2的超晶胞結(jié)構(gòu)圖

    表1 結(jié)構(gòu)優(yōu)化后Fe(S1?xAsx)2 (x=0, 0.0156, 0.03125, 0.0625)的折合晶胞參數(shù)和總能

    Total energy: Total energy of Fe(S1?xAs)2convert to energy in 2×2×2 super cell.

    2.2 電子結(jié)構(gòu)

    2.2.1 電子態(tài)密度

    圖2所示為FeS2單胞結(jié)構(gòu)優(yōu)化后的總態(tài)密度及Fe、S各亞層電子的分波態(tài)密度圖。在FeS2晶體結(jié)構(gòu)中,處于八面體中心的鐵原子由于晶體場(chǎng)作用Fe 3d軌道分裂成2g和g*兩部分,因此Fe 3d軌道間斷地分布在價(jià)帶和導(dǎo)帶中。Fe 3d的2g部分主要出現(xiàn)在上部?jī)r(jià)帶(靠近費(fèi)米能級(jí)),與S 3p態(tài)重疊較小,表明它們之間的成鍵作用較弱。由圖2(b)可知,F(xiàn)e的3d2g沒有實(shí)質(zhì)上的成鍵?反鍵劈裂,其本質(zhì)原因是軌道的對(duì)稱性和能量不匹配,因而其也被稱作非鍵軌道[27]。 Fe 3dg*和S 3p形成反鍵,反鍵峰約在2.5 eV處, Fe 3dg*與S 3p態(tài)密度重疊程度較高,表明這二者之間作用較強(qiáng)。

    在態(tài)密度圖中,軌道的貢獻(xiàn)就意味著軌道重疊,重疊度越大,軌道雜化就越強(qiáng)烈。從圖2(b)中可以看出,F(xiàn)e的3d、4s軌道以及S的3s、3p軌道在?3.2 eV的同一能量位置出現(xiàn)了態(tài)密度峰,說明這兩個(gè)軌道發(fā)生了重疊形成了一個(gè)新的軌道,從雜化峰的寬度來看Fe、S雜化峰離域性強(qiáng),軌道雜化作用較強(qiáng)。從圖2可以看出,在能量為?17.5~?9 eV范圍內(nèi),態(tài)密度幾乎全部由S 3s態(tài)貢獻(xiàn),僅有少量的Fe 4s和S 3p態(tài)貢獻(xiàn);在能量為?7~?1.5 eV范圍內(nèi),態(tài)密度主要由Fe 3d態(tài)和S 3p態(tài)貢獻(xiàn),其中S 3p態(tài)貢獻(xiàn)最大,而Fe 4s和S 3s態(tài)的貢獻(xiàn)相對(duì)較小;費(fèi)米能級(jí)附近價(jià)帶區(qū)的態(tài)密度,則主要由S的3p態(tài)和非成鍵Fe 3d態(tài)組成,且大部分由Fe 3d態(tài)貢獻(xiàn)。費(fèi)米能級(jí)附近導(dǎo)帶區(qū)0~4 eV范圍內(nèi)的態(tài)密度主要由S 3p和Fe的反鍵3dg*態(tài)貢獻(xiàn),而S 3s態(tài)貢獻(xiàn)很小。

    從圖2(b)中還可以看出,S 3s態(tài)對(duì)深部?jī)r(jià)帶的貢獻(xiàn)較大,而S 3p態(tài)則主要對(duì)淺部?jī)r(jià)帶產(chǎn)生貢獻(xiàn)。Fe 4s對(duì)整個(gè)能帶范圍態(tài)密度的貢獻(xiàn)相對(duì)較少,而費(fèi)米能級(jí)附近的態(tài)密度主要來自Fe 3d的貢獻(xiàn),這和本征態(tài)FeS2的總態(tài)密度計(jì)算結(jié)果相符合。

    圖2 FeS2單胞總態(tài)密度和各原子態(tài)密度

    As摻雜置換S后得到的Fe(S1?As)2態(tài)密度如圖3所示。從圖3中可以看出,As置換FeS2中的S原子后,由于As和S原子核外電子排布的相似性,對(duì)FeS2的電子態(tài)密度影響不是很大,價(jià)帶頂和導(dǎo)帶底仍主要由Fe 3d態(tài)決定,As的s層和p層電子分別與S的s層和p層電子混合,共同構(gòu)成Fe(S1?As)2的能態(tài)密度。

    圖3 Fe(S1?xAsx)2的態(tài)密度圖

    As的摻雜使FeS2的電子能帶在淺部?jī)r(jià)帶產(chǎn)生了由它的p態(tài)貢獻(xiàn)的雜質(zhì)能級(jí),并在更低的能量范圍內(nèi)產(chǎn)生了由其s態(tài)貢獻(xiàn)的雜質(zhì)能級(jí),并且隨著體系A(chǔ)s摻雜濃度的增大,這種貢獻(xiàn)越來越明顯。

    2.2.2 能帶結(jié)構(gòu)

    圖4所示為As摻雜前后FeS2和Fe(S1?As)2的能帶結(jié)構(gòu),圖中G、F、Q、Z、G為FeS2晶體第一布里淵區(qū)高對(duì)稱點(diǎn)。為了進(jìn)行對(duì)比,將理想FeS2的第一布里淵區(qū)沿對(duì)稱方向的點(diǎn)取樣設(shè)置為與摻雜后的體系一致。在計(jì)算中所考慮的禁帶寬度是從價(jià)帶最高點(diǎn)到導(dǎo)帶最低點(diǎn)之間的距離,計(jì)算得到理想FeS2的禁帶寬度為0.55 eV,小于實(shí)驗(yàn)值0.95 eV,一般認(rèn)為出現(xiàn)這種結(jié)果的主要原因是作為基態(tài)理論的DFT,其Kohn-Sham方程的本征值不能給出系統(tǒng)的激發(fā)態(tài)能量,導(dǎo)致位于導(dǎo)帶部分的電子態(tài)能量值比實(shí)驗(yàn)值小,從而帶隙偏小,一般情況下能帶值與理論值的誤差范圍為30%~50%[28]。

    As摻雜后,F(xiàn)eS2的帶隙為0.43 eV,這表明As摻雜后FeS2的禁帶寬度變小。從圖4(b)、(d)、(f)可以明顯看出,價(jià)帶的上方出現(xiàn)了明顯的雜質(zhì)能帶,對(duì)比圖4(e)、(f)可以看出,理想FeS2費(fèi)米能級(jí)附近的導(dǎo)帶帶寬為0.55~4.10 eV,而摻雜As后費(fèi)米能級(jí)附近的導(dǎo)帶帶寬變?yōu)?.45~4.41 eV,導(dǎo)帶不僅變寬而且還下移了0.1 eV,也說明摻雜后FeS2的禁帶寬度變小。

    2.3 光學(xué)性質(zhì)計(jì)算結(jié)果與討論

    為了研究As摻雜對(duì)FeS2的光學(xué)性質(zhì)的影響,對(duì)Fe(S1?As)2的光學(xué)躍遷特性進(jìn)行了研究。在線性響應(yīng)范圍內(nèi)固體宏觀光學(xué)響應(yīng)函數(shù)通??梢杂晒獾慕殡姾瘮?shù)() =1() + i2() 或折射率() =() +i() 來描述,體系在較小波矢下對(duì)光場(chǎng)的線性響應(yīng)由介電函數(shù)的虛部2()決定,它可以通過占據(jù)態(tài)和未占據(jù)態(tài)的動(dòng)量躍遷矩陣元計(jì)算得到。

    2.3.1 復(fù)介電函數(shù)

    介電常數(shù)的虛部2()對(duì)任何材料來說都是非常重要的,它作為溝通帶間躍遷微觀物理過程與固體電子結(jié)構(gòu)的橋梁,反映了固體能帶結(jié)構(gòu)及其他各種光譜信息。圖5所示為計(jì)算得到的未摻雜FeS2和Fe(S1?xAs)2的介電函數(shù)實(shí)部和虛部。

    由圖5(a)可知,未摻雜的FeS2的靜態(tài)介電常數(shù)1(0)=19.6,隨著As摻雜量的增加,F(xiàn)e(S1?xAs)2的靜態(tài)介電常數(shù)1(0)的值逐漸增大。當(dāng)<0.46 eV或>9.5 eV時(shí)Fe(S1?xAs)2的1值比未摻雜時(shí)的大,而當(dāng)0.46 eV<<9.5 eV時(shí),F(xiàn)e(S1?xAs)2的1值比未摻雜的FeS2的小。在<7.5 eV的范圍內(nèi),1值隨著摻雜量的增加逐漸增大,而>7.5 eV后,F(xiàn)e(S1?xAs)2的1值隨摻雜量的增加而減小并且逐漸接近未摻雜FeS2的1值。

    圖4 FeS2和Fe(S1?xAsx)2 (x=0.0156, 0.0312, 0.0625) 對(duì)應(yīng)的相同超胞的能帶結(jié)構(gòu)圖

    圖5 FeS2和Fe(S1?xAsx)2的介電函數(shù)實(shí)部和虛部

    由圖5(b)可知,對(duì)于未摻雜的FeS2在<0.28 eV的能量范圍內(nèi)2=0,2()在圖所示的能量范圍內(nèi)有3個(gè)明顯的介電峰,對(duì)應(yīng)的能量依次為2.63、7.26和9.48 eV。位于能量為2.63 eV的第一個(gè)介電峰來自S的s和p軌道發(fā)生雜化后電子到導(dǎo)帶的躍遷,能量為7.26和9.48 eV的介電峰來源于價(jià)帶中部S 3p態(tài)到導(dǎo)帶Fe 3d態(tài)的電子躍遷。當(dāng)As摻雜后,可以明顯看到2()的第一躍遷峰向低能方向移動(dòng),并且第一躍遷峰的值隨摻雜量的增加而明顯增大。當(dāng)<1.28 eV時(shí),F(xiàn)e(S1?xAs)2的光躍遷強(qiáng)度遠(yuǎn)大于未摻雜的FeS2,且光躍遷強(qiáng)度隨摻雜量的增大而增大,說明As摻雜能夠增強(qiáng)電子在<1.28 eV的低能量范圍內(nèi)的光學(xué)躍遷。在>2.5 eV的能量范圍Fe(S1?xAs)2的光躍遷強(qiáng)度小于未摻雜的FeS2,且光躍遷強(qiáng)度隨著As摻雜量的增加而增大。

    2.3.2 復(fù)折射率

    由公式1=2?2和2=2可以從復(fù)介電函數(shù)得到FeS2的復(fù)折射率。圖6所示為FeS2和Fe(S1?xAs)2的折射率和消光系數(shù)。對(duì)于FeS2理想體系,折射率(見圖6(a))在=1.66 eV處取得最大值,當(dāng)>1.66 eV時(shí),折射率快速下降,在6.3~9.4 eV范圍內(nèi)出現(xiàn)兩個(gè)變化平緩的峰,表明FeS2在這一能量范圍內(nèi)呈現(xiàn)出一定的金屬反射特性。消光系數(shù)(見圖6(b))在3.88 eV處取得最大值,隨后隨著能量的增加而減小,并在能量約為20 eV時(shí)減小到零。

    由圖6(a)可知,當(dāng)As摻雜量從本征態(tài)增加到=0.0625時(shí),折射率的值從4.43增加到5.99。當(dāng)0.56~14.1 eV時(shí),F(xiàn)e(S1-xAs)2的折射率小于FeS2的折射率,且折射率隨摻雜量的增加而有小幅增大。在<0.56 eV和>14.1 eV的范圍內(nèi),F(xiàn)e(S1?xAs)2的折射率大于FeS2的折射率,且隨摻雜量的增大而減小。消光系數(shù)(見圖6(b))在<1.42 eV范圍內(nèi)出現(xiàn)了新的峰值,并且As摻雜后消光系數(shù)較本征態(tài)均有所增大。在3.48~10.54 eV的能量范圍內(nèi)出現(xiàn)了較強(qiáng)的波動(dòng)峰,當(dāng)>10.54 eV時(shí),隨摻雜量的增加逐漸增大。

    2.3.3 光電導(dǎo)率

    光電導(dǎo)率是光電子材料的一個(gè)重要參數(shù),它描述了光照引起半導(dǎo)體電導(dǎo)率改變的現(xiàn)象。圖7所示為Fe(S1?xAs)2和FeS2的光電導(dǎo)率實(shí)部1()。從圖7可以看出,F(xiàn)eS2的光電導(dǎo)率與介電函數(shù)的虛部是對(duì)應(yīng)的,1()呈先上升后下降的變化趨勢(shì),其中在=3.60 eV處達(dá)到第一個(gè)峰值,在=7.49和9.65 eV處分別達(dá)到第二、第三個(gè)峰值,且在=3.60 eV處取得最大值。當(dāng)摻入As后,除Fe(S1?xAs)2(=0.0625)在4.8~5.5 eV范圍內(nèi)的光電導(dǎo)率大于本征態(tài)外,摻雜體系的光電導(dǎo)率均小于本征態(tài)FeS2的光電導(dǎo)率,同時(shí),在>2.66 eV時(shí)隨著摻雜量的增加,F(xiàn)e(S1?xAs)2的光電導(dǎo)率逐漸增加。Fe(S1?xAs)2的光電導(dǎo)率最大峰值出現(xiàn)在2.9~3.3 eV的范圍內(nèi),隨著As摻雜濃度的增大,摻雜體系的電導(dǎo)率逐漸增大,峰值出現(xiàn)的位置也隨摻雜量的增加向高能方向移動(dòng)。

    圖7 FeS2和Fe(S1?xAsx)2的光電導(dǎo)率

    3 結(jié)論

    1) 對(duì)不同As摻雜濃度FeS2的幾何結(jié)構(gòu)計(jì)算結(jié)果表明,As摻雜FeS2引起了晶格參數(shù)的增大,隨As摻雜濃度增大,F(xiàn)e(S1?xAs)2(=0.0156,0.0312,0.0625)的晶格常數(shù)和晶格體積均逐漸增大,這是由于As的共價(jià)半徑大于S的共價(jià)半徑所致。

    2) 理想FeS2的優(yōu)化表明其禁帶寬度為0.55 eV,費(fèi)米能級(jí)附近導(dǎo)帶區(qū)0~4 eV范圍內(nèi)的態(tài)密度主要由S 3p和Fe的反鍵3dg*態(tài)貢獻(xiàn),而S 3s態(tài)貢獻(xiàn)很小。As摻雜使得FeS2導(dǎo)帶下移,禁帶寬度減小,并且在?10.4~?9.5 eV的淺部?jī)r(jià)帶產(chǎn)生了由As的p態(tài)貢獻(xiàn)的雜質(zhì)能級(jí)。通過對(duì)單個(gè)原子態(tài)密度的分析發(fā)現(xiàn),As的摻雜對(duì)黃鐵礦費(fèi)米能級(jí)處電子態(tài)密度有貢獻(xiàn),且對(duì)Fe和S原子的態(tài)密度分布產(chǎn)生了影響。

    3) 對(duì)光學(xué)性質(zhì)計(jì)算結(jié)果分析可知,F(xiàn)e(S1?xAs)2的靜態(tài)介電常數(shù)、第一介電峰、折射率隨著As摻雜量的增大而逐漸增大,平均反射效應(yīng)的減弱有效地提高了其光電轉(zhuǎn)換效率。因此,可以通過對(duì)FeS2的As摻雜以及調(diào)整摻雜濃度來得到不同光電特性的FeS2電催化活性材料。

    [1] SONG Xiao-hui, WANG Min-qiang, ZHANG Hao, DENG Jian-ping, YANG Zhi, RAN Chen-xin, YAO Xi. Morphologically controlled electrodeposition of CdSe on mesoporous TiO2film for quantum dot-sensitized solar cells[J]. Electrochimica Acta, 2013, 108(10): 449?457.

    [2] SUEZAKI T, CHEN J I L, HATAYAMA T, FUYUKI T, OZIN G A. Electrical properties of p-type and n-type doped inverse silicon opals-towards optically amplified silicon solar cells[J]. Applied Physics Letter, 2010, 96(24): 242102-1?242102-3.

    [3] ZHAI Guang-mei, BEZRYADINA A, BREEZE A J, ZHANG Dao-li, ALERS G B, CARTER S A. Air stability of TiO2/PbS colloidal nanoparticle solar cells and its impact on power efficiency[J]. Applied Physics Letter, 2011, 99(6): 063512-1?063512-3.

    [4] ALTERMATT P P, KIESEWETTER T, ELLMER K, TRIBUTSCH H. Specifying targets of future research in photovoltaic devices containing pyrite (FeS2) by numerical modeling[J]. Solar Energy Materials and Solar Cells, 2002, 71(2): 181?195.

    [5] MIGUEL C A, FABER M S, TAN Yi-zheng, HAMERS R J, JIN Song. Synthesis and properties of semiconducting iron pyrite (FeS2) nanowires[J]. Nano Letter, 2012, 12(4): 1977?1982.

    [6] DOUGLAS A, CARTER R, OAKES L, SHARE K, COHN A P, PINT C L. Ultrafine iron pyrite (FeS2) nanocrystals improve sodium–sulfur and lithium–sulfur conversion reactions for efficient batteries[J]. ACS Nano, 2015, 9(11):11156?11165.

    [7] WALTER M, KRAVCHYK K V, IBANEZ M, KOVALENKO M V. Efficient and inexpensive sodium–magnesium hybrid battery[J]. Chemistry Material, 2015, 27(21): 7452?7458.

    [8] CAO F, PAN G X, CHEN J, ZHANG Y J, XIA X H. Synthesis of pyrite/carbon shells on cobalt nanowires forming core/branch arrays as high-performance cathode for lithium ion batteries[J]. Journal of Power Sources, 2016, 303(30): 35?40.

    [9] QIU Wen-da, XIA Jian, ZHONG Hai-min, HE Shen-xian, LAI Shu-hui, CHEN Liu-ping. L-cysteine-assisted synthesis of cubic pyrite/nitrogen-doped graphene composite as anode material for lithium-ion batteries[J]. Electrochimica Acta, 2014, 137(10): 197?205.

    [10] 黎佩珊, 孫明, 鄭育英, 徐前忠. 不同純度的正極材料黃鐵礦的表征及其性能研究[J]. 廣東化工, 2012, 39(16): 1?2. LI Pei-shan, SUN Ming, ZHENG Yu-ying, XU Qian-zhong. Characterization and electritic properties of pyrite with different purity[J]. Guangdong Chemical, 2012, 39(16): 1?2.

    [11] 王大剛, 范力仁, 王圣平, 范疇, 何明中, 栗海峰. 黃鐵礦作為鋰電池正極材料的電化學(xué)性能[J]. 材料導(dǎo)報(bào), 2012, 26(9): 93?96. WANG Da-gang, FAN Li-ren, WANG Sheng-ping, FAN Chou, HE Ming-zhong, LI Hai-feng. Electrochemical properties of pyrite as lithium battery cathode materials[J]. Materials Review, 2012, 26(9): 93?96.

    [12] 楊兆堂, 劉效疆. 水熱法合成FeS2材料及其在熱電池中的性能[J]. 材料保護(hù), 2013, 46(S1): 51?53. YANG Zhao-tang, LIU Xiao-jiang. Preparing pyrite via hydrothermal method for thermal battery applications[J]. Materials Protection, 2013, 46(S1): 51?53.

    [13] 龍 飛, 神 征, 莫淑一, 鄒正光. 黃鐵礦(FeS2)薄膜制備和摻雜改性的研究進(jìn)展[J]. 半導(dǎo)體光電, 2015, 36(6): 857?862. LONG Fei, SHEN Zheng, MO Shu-yi, ZOU Zheng-guang. Research progress on deposition and doping techniques of FeS2thin film[J]. Semiconductor Optoelectronics, 2015, 36(6): 857?862.

    [14] CHATZITHEODOROU G, FIECHTER S, KONENKAMP R, KUNST M, JAEGERMANN W, TRIBUTSCH H. Thin photoactive FeS2(pyrite) films[J]. Materials Research Bulletin, 1986, 21(12): 1481?1487.

    [15] LEHNER S W, NEWMAN N, van SCHILFGAARDE M, BANDYOPADHYAY S, SAVAGE K, BUSECK P R. Defect energy levels and electronic behavior of Ni-, Co-, and As-doped synthetic pyrite (FeS2)[J]. Journal of Applied Physics, 2012, 111(8): 083717-1?083717-7.

    [16] FERRER I J, de las HERAS C, SANCHEZ C. Physical properties of Cu-doped FeS2pyrite thin films[J]. Applied Surface Science, 1993, 70(2): 588?592.

    [17] YU Qi-kai, CAI Shu, JIN Zheng-guo, YAN Zi-peng . Evolutions of composition, microstructure and optical properties of Mn-doped pyrite (FeS2) films prepared by chemical bath deposition[J]. Materials Research Bulletin, 2013, 48(9): 3601?3606.

    [18] 李玉瓊, 陳建華, 陳 曄, 郭 進(jìn). 黃鐵礦(100)表面性質(zhì)的密度泛函理論計(jì)算及其對(duì)浮選的影響[J]. 中國(guó)有色金屬學(xué)報(bào), 2011, 21(4): 919?926. LI Yu-qiong, CHEN Jian-hua, CHEN Ye, GUO Jin. Density functional theory calculation of surface properties of pyrite (100) with implications for flotation[J]. The Chinese Journal of Nonferrous Metals, 2011, 21(4): 919?926.

    [19] CHEN Jian-hua, CHEN Ye, LI Yu-qiong. Effect of vacancy defects on electronic properties and activation of sphalerite(110) surface by first-principles[J]. Transactions of Nonferrous Metals Society of China, 2010, 21(3): 502?506.

    [20] 陳建華, 李玉瓊, 衷水平, 郭 進(jìn). 含空位缺陷黃鐵礦(100)表面吸附氫氧根和羥基鈣的量子化學(xué)研究[J]. 中國(guó)有色金屬學(xué)報(bào), 2013, 23(3): 859?865. CHEN Jian-hua, LI Yu-qiong, ZHONG Shui-ping, GUO Jin. Quantum chemical study of adsorption of hydroxyl and hydroxyl calcium on pyrite (100) surface bearing vacancy defects[J]. The Chinese Journal of Nonferrous Metals, 2013, 23(3): 859?865.

    [21] 陳建華, 鐘建蓮, 李玉瓊, 陳 曄, 郭 進(jìn). 黃鐵礦、白鐵礦和磁黃鐵礦的電子結(jié)構(gòu)及可浮性[J]. 中國(guó)有色金屬學(xué)報(bào), 2011, 21(7): 1719?1727. CHEN Jian-hua, ZHONG Jian-lian, LI Yu-qiong, CHEN Ye, GUO Jin. Electronic structures and floatability of pyrite, marcasite and pyrrhotite[J]. The Chinese Journal of Nonferrous Metals, 2011, 21(7): 1719?1727.

    [22] CLARK S J, SEGALL M D, PICKARD C J, HASNIP P J, PROBERT M J, REFSON K, PAYNE M C. First principles methods using CASTEP[J]. Zeitischrift fuer Kristallographie, 2005, 220(5/6): 567?570.

    [23] VANDERBILT D. Soft self-consistent pseudopotentials in generalized eigenvalue formalism[J]. Physical Review B, 1990, 41(11): 7892?7895.

    [24] MONKHORST H J, PACK J D. Special points for Brillouin-zone integrations[J]. Physical Review B, 1976, 13(12): 5188?5192.

    [25] PRINCE K C, MATTEUCCI M, KUEPPER K, CHIUZBAIAN S G, BARKOWSKI S, NEUMANN M. Core-level spectroscopic study of FeO and FeS2[J]. Physical Review B, 2005, 71(8): 085102-1?085102-9.

    [26] OERTZEN G U , JONES R T, GERSON A R. Electronic and optical properties of Fe, Zn and Pb sulfides[J]. Physics and Chemistry of Minerals, 2005, 32(4): 255?268.

    [27] 陳建華. 硫化礦物浮選固體物理研究[M]. 長(zhǎng)沙: 中南大學(xué)出版社, 2015. CHEN Jian-hua. The solide physics of sulphide minerals flotation[M]. Changsha: Central South University Press, 2015.

    [28] JONES R O, GUNNARSSON O. The density functional formalism, its applications and prospects[J]. Review of Modern Physics, 1989, 61(3): 689?746.

    Effect of As doping concentration on electronic structure and optical properties of FeS2

    MA Wan-kun1, CHEN Jian-hua2, ZHANG Guo-fan1, LI Yu-qiong2, FENG Qi-ming1

    (1. School of Minerals Processing and Bioengineering, Central South University, Changsha 410083, China;2. College of Resources and Metallurgy, Guangxi University, Nanning 530004, China)

    Aimed at verifying the effect of As doping concentration on FeS2, geometrical structure, electronic structure and optical properties of As-doped FeS2were calculated using pseudo-potential plane-wave method of the first principle based on the density function theory. Three supercells (2×2×2, 2×2×1 and 2×1×1) were introduced to achieve different doping concentrations (1.93%, 3.82% and 7.48%, mass fraction). The results of geometrical structure, electronic structure show that lattice constant and volume of unit cell increase with As addition, the conduction band moves negatively while the width of forbidden band decreases. Moreover, the density of states in the range of ?10.4 ? ?9.5 eV exhibits impurity level composed of As p. The calculation of optical properties reveals the static dielectric constant, refractive index and photo conductivity aggrandize with the ever-increasing doping concentration, which indicates FeS2doped As possesses stronger absorption capacity of light and higher photoelectric conversion efficiency. Therefore, the calculation lays some theory foundation for exploitation and development of FeS2photoelectric materials.

    As-doping; pyrite; optical and electrical property; first principle

    (編輯 何學(xué)鋒)

    Project(2014CB643402) supported by National Basic Research Development Program of China; Project (51304054) supported by the National Natural Science Foundation of China; Project (2013GXNSFBA019259) supported by the Guangxi Natural Science Foundation, China

    2016-02-24; Accepted date:2016-07-05

    ZHANG Guo-fan; Tel: +86-731-88830913; E-mail: zhangguofan2002@qq.com

    1004-0609(2017)-02-0414-09

    O472

    A

    國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展規(guī)劃資助項(xiàng)目(2014CB643402);國(guó)家自然科學(xué)基金資助項(xiàng)目(51304054);廣西自然科學(xué)基金資助項(xiàng)目(2013GXNSFBA019259)

    2016-02-24;

    2016-07-05

    張國(guó)范,教授,博士;電話:0731-88830913;E-mail: zhangguofan2002@qq.com

    猜你喜歡
    價(jià)帶導(dǎo)帶晶胞
    晶胞考查角度之探析
    外加電場(chǎng)和磁場(chǎng)對(duì)外爾半金屬能帶結(jié)構(gòu)的調(diào)控*
    四步法突破晶體密度的計(jì)算
    關(guān)于某型系列產(chǎn)品“上下導(dǎo)帶二次收緊”合并的工藝技術(shù)
    炮射彈藥滑動(dòng)導(dǎo)帶環(huán)結(jié)構(gòu)設(shè)計(jì)
    淺談晶胞空間利用率的計(jì)算
    導(dǎo)帶式數(shù)碼噴射印花機(jī)運(yùn)行過程中常見疵病
    打捆機(jī)導(dǎo)帶槽裝置及捆帶盤脹縮裝置的改進(jìn)
    交換場(chǎng)和非共振光對(duì)單層MoS2能帶結(jié)構(gòu)的調(diào)控*
    Modification of CaO-based sorbents prepared from calcium acetate for CO2 capture at high temperature☆
    人妻久久中文字幕网| 国产av在哪里看| 又紧又爽又黄一区二区| 三级男女做爰猛烈吃奶摸视频| 国产亚洲精品久久久久久毛片| 婷婷精品国产亚洲av| 全区人妻精品视频| 精品国产亚洲在线| 精品国产超薄肉色丝袜足j| 日韩免费av在线播放| 亚洲乱码一区二区免费版| 国产一区二区在线观看日韩 | 岛国视频午夜一区免费看| 欧美又色又爽又黄视频| 不卡一级毛片| av福利片在线观看| 国产真人三级小视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 久99久视频精品免费| 亚洲成av人片在线播放无| 老司机在亚洲福利影院| 宅男免费午夜| 久久九九热精品免费| 久久久久精品国产欧美久久久| 岛国在线观看网站| 国产真人三级小视频在线观看| 免费观看人在逋| 狂野欧美白嫩少妇大欣赏| 亚洲av成人精品一区久久| 一个人看的www免费观看视频| 国产野战对白在线观看| 69人妻影院| 首页视频小说图片口味搜索| 国产毛片a区久久久久| 亚洲成人中文字幕在线播放| 一个人观看的视频www高清免费观看| 亚洲精品456在线播放app | www.熟女人妻精品国产| 全区人妻精品视频| 国产黄a三级三级三级人| 久久午夜亚洲精品久久| 可以在线观看毛片的网站| 欧美日韩黄片免| 久久久久久久久久黄片| 国产淫片久久久久久久久 | 国产精品99久久久久久久久| 国产一区二区在线观看日韩 | 午夜免费男女啪啪视频观看 | 欧美极品一区二区三区四区| 亚洲av二区三区四区| 午夜福利欧美成人| 一边摸一边抽搐一进一小说| 怎么达到女性高潮| 波多野结衣高清无吗| 色综合站精品国产| 国产视频内射| 成人高潮视频无遮挡免费网站| 久久人人精品亚洲av| 乱人视频在线观看| 免费看日本二区| 国产免费男女视频| 国产不卡一卡二| 丝袜美腿在线中文| 成人性生交大片免费视频hd| 啪啪无遮挡十八禁网站| 90打野战视频偷拍视频| 成人特级av手机在线观看| 最近视频中文字幕2019在线8| 国产精品一区二区三区四区久久| 亚洲精品日韩av片在线观看 | 亚洲国产欧美人成| 亚洲无线在线观看| 日本a在线网址| 1024手机看黄色片| 国产精品亚洲av一区麻豆| 天堂影院成人在线观看| 人人妻人人看人人澡| 日韩精品中文字幕看吧| 亚洲精品亚洲一区二区| 午夜免费观看网址| 黄色视频,在线免费观看| 精品一区二区三区视频在线 | 免费看美女性在线毛片视频| 精品日产1卡2卡| 欧美激情久久久久久爽电影| 精品乱码久久久久久99久播| 国内精品久久久久久久电影| 手机成人av网站| 国产高清激情床上av| 免费看十八禁软件| 99热精品在线国产| 亚洲成人精品中文字幕电影| 亚洲精品久久国产高清桃花| 精品久久久久久久久久久久久| 九九在线视频观看精品| 日本精品一区二区三区蜜桃| 欧美色欧美亚洲另类二区| 好看av亚洲va欧美ⅴa在| 日韩成人在线观看一区二区三区| 久久久久国内视频| 国产成人a区在线观看| 999久久久精品免费观看国产| 国产av不卡久久| 热99re8久久精品国产| 午夜免费男女啪啪视频观看 | 最近最新中文字幕大全电影3| 亚洲 欧美 日韩 在线 免费| 亚洲狠狠婷婷综合久久图片| 国产久久久一区二区三区| 在线观看66精品国产| 在线免费观看的www视频| 真人一进一出gif抽搐免费| 日日夜夜操网爽| 大型黄色视频在线免费观看| 欧美成人性av电影在线观看| 亚洲不卡免费看| 在线观看免费视频日本深夜| 欧美性猛交黑人性爽| 最新在线观看一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 91av网一区二区| 免费人成视频x8x8入口观看| 亚洲欧美日韩高清在线视频| 国产v大片淫在线免费观看| 久久久久久久久中文| 精品人妻1区二区| 亚洲av电影在线进入| 国产成人欧美在线观看| 变态另类成人亚洲欧美熟女| 久久久国产精品麻豆| 欧美3d第一页| 国产亚洲精品久久久久久毛片| АⅤ资源中文在线天堂| 国产v大片淫在线免费观看| 变态另类成人亚洲欧美熟女| 97超视频在线观看视频| 国产精品爽爽va在线观看网站| 日本与韩国留学比较| 国产三级中文精品| 91av网一区二区| 中文字幕熟女人妻在线| 亚洲av免费高清在线观看| 日韩欧美在线二视频| 内射极品少妇av片p| 国产爱豆传媒在线观看| 男女床上黄色一级片免费看| 成人无遮挡网站| 欧美另类亚洲清纯唯美| 好看av亚洲va欧美ⅴa在| 97碰自拍视频| 成人特级av手机在线观看| 日本 欧美在线| 亚洲在线自拍视频| 亚洲国产精品999在线| 女人十人毛片免费观看3o分钟| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲av二区三区四区| 中文字幕人妻丝袜一区二区| 欧美成人免费av一区二区三区| 在线观看66精品国产| 亚洲精品粉嫩美女一区| 亚洲男人的天堂狠狠| 动漫黄色视频在线观看| 亚洲av免费在线观看| 欧美成人性av电影在线观看| 一个人看的www免费观看视频| 啦啦啦韩国在线观看视频| 人妻丰满熟妇av一区二区三区| 国产精品99久久久久久久久| 老汉色∧v一级毛片| 婷婷精品国产亚洲av在线| 午夜免费观看网址| 两个人视频免费观看高清| 日本撒尿小便嘘嘘汇集6| 窝窝影院91人妻| 精品国产超薄肉色丝袜足j| 色综合亚洲欧美另类图片| 亚洲七黄色美女视频| 首页视频小说图片口味搜索| 91九色精品人成在线观看| 神马国产精品三级电影在线观看| 国产一区二区在线av高清观看| 久久九九热精品免费| 亚洲天堂国产精品一区在线| 国产精品99久久99久久久不卡| 日韩欧美国产一区二区入口| 欧美大码av| 国产v大片淫在线免费观看| 高潮久久久久久久久久久不卡| 久久精品国产亚洲av涩爱 | 97超级碰碰碰精品色视频在线观看| 久久久久久久久中文| 成人国产一区最新在线观看| 国产一区二区在线av高清观看| 欧美色视频一区免费| 国产色婷婷99| 国产免费一级a男人的天堂| 亚洲18禁久久av| 给我免费播放毛片高清在线观看| 亚洲天堂国产精品一区在线| 在线观看66精品国产| 国产精品一及| 真人一进一出gif抽搐免费| 99热精品在线国产| 国产爱豆传媒在线观看| 国产成人福利小说| 国产麻豆成人av免费视频| 日本三级黄在线观看| 久久精品91无色码中文字幕| 最新在线观看一区二区三区| 免费av不卡在线播放| 人人妻,人人澡人人爽秒播| 波多野结衣高清作品| 禁无遮挡网站| 国产av一区在线观看免费| 亚洲精品一卡2卡三卡4卡5卡| 男人舔女人下体高潮全视频| 亚洲五月天丁香| 欧美另类亚洲清纯唯美| 成人高潮视频无遮挡免费网站| 美女免费视频网站| 欧美zozozo另类| 蜜桃久久精品国产亚洲av| 黄色日韩在线| 国内揄拍国产精品人妻在线| 在线观看午夜福利视频| 观看免费一级毛片| 午夜福利在线观看免费完整高清在 | 亚洲最大成人手机在线| 精品一区二区三区av网在线观看| 午夜两性在线视频| 久久久久久久久大av| 男插女下体视频免费在线播放| 真实男女啪啪啪动态图| 国产精品久久久久久久电影 | 一级黄色大片毛片| 精品免费久久久久久久清纯| 国产亚洲精品综合一区在线观看| АⅤ资源中文在线天堂| 日韩大尺度精品在线看网址| 欧美三级亚洲精品| 最近最新中文字幕大全免费视频| 99久久99久久久精品蜜桃| 丰满人妻一区二区三区视频av | 不卡一级毛片| 久久久成人免费电影| 久久精品91无色码中文字幕| 亚洲在线自拍视频| 一本久久中文字幕| 成人av在线播放网站| 中文字幕人妻熟人妻熟丝袜美 | 少妇的逼好多水| 日韩亚洲欧美综合| 国产亚洲精品综合一区在线观看| 亚洲自拍偷在线| 久久亚洲真实| 亚洲在线自拍视频| 91字幕亚洲| 18禁在线播放成人免费| 欧美中文日本在线观看视频| 精品福利观看| 久久精品国产清高在天天线| 好男人在线观看高清免费视频| 国产精品野战在线观看| 97人妻精品一区二区三区麻豆| 亚洲精华国产精华精| 国产探花在线观看一区二区| 色在线成人网| 亚洲黑人精品在线| 久久国产精品人妻蜜桃| 黑人欧美特级aaaaaa片| 国产精品久久久久久久久免 | 亚洲黑人精品在线| 熟女少妇亚洲综合色aaa.| 白带黄色成豆腐渣| 久久久久久久久久黄片| 最近在线观看免费完整版| 啦啦啦观看免费观看视频高清| 夜夜看夜夜爽夜夜摸| 超碰av人人做人人爽久久 | 成人国产综合亚洲| 国产一区二区在线av高清观看| 午夜激情欧美在线| 又紧又爽又黄一区二区| 久久久精品欧美日韩精品| 亚洲精品乱码久久久v下载方式 | 禁无遮挡网站| 人人妻人人看人人澡| 18禁黄网站禁片免费观看直播| 午夜免费激情av| 亚洲中文日韩欧美视频| 亚洲成人中文字幕在线播放| 亚洲五月婷婷丁香| 看片在线看免费视频| 国产一区二区三区视频了| 综合色av麻豆| 欧美一级毛片孕妇| 亚洲中文字幕一区二区三区有码在线看| 免费观看人在逋| 国产高清videossex| 亚洲中文字幕一区二区三区有码在线看| 久久久精品欧美日韩精品| 午夜久久久久精精品| 国产免费一级a男人的天堂| 香蕉丝袜av| 午夜久久久久精精品| 亚洲精品在线观看二区| 听说在线观看完整版免费高清| 国产精品 欧美亚洲| 少妇高潮的动态图| 久久精品影院6| 老司机午夜福利在线观看视频| 欧美日韩精品网址| 在线免费观看的www视频| 给我免费播放毛片高清在线观看| 亚洲国产日韩欧美精品在线观看 | 99热这里只有精品一区| 91字幕亚洲| 日本 欧美在线| 91在线观看av| eeuss影院久久| 免费观看的影片在线观看| 久久久久久久久中文| 国产一区在线观看成人免费| 国产成人av教育| 老熟妇乱子伦视频在线观看| 亚洲激情在线av| 麻豆成人午夜福利视频| 中文资源天堂在线| 精品乱码久久久久久99久播| www国产在线视频色| 国产99白浆流出| 哪里可以看免费的av片| 99久久精品国产亚洲精品| 国产精品亚洲一级av第二区| 欧美性猛交╳xxx乱大交人| 久久久精品大字幕| 午夜福利在线在线| 丝袜美腿在线中文| 日本一本二区三区精品| 亚洲电影在线观看av| 午夜福利免费观看在线| 欧美国产日韩亚洲一区| 1000部很黄的大片| 天美传媒精品一区二区| 国产成人a区在线观看| 欧美av亚洲av综合av国产av| 日韩高清综合在线| 欧美日本视频| 国产精品久久久人人做人人爽| 别揉我奶头~嗯~啊~动态视频| 无限看片的www在线观看| 99热只有精品国产| 亚洲第一欧美日韩一区二区三区| 999久久久精品免费观看国产| 制服丝袜大香蕉在线| 高潮久久久久久久久久久不卡| 亚洲国产精品sss在线观看| 1024手机看黄色片| 国产乱人伦免费视频| 舔av片在线| 国产成人a区在线观看| 丁香欧美五月| 桃色一区二区三区在线观看| 成人三级黄色视频| 99热这里只有精品一区| 色综合亚洲欧美另类图片| 国产av一区在线观看免费| 久久婷婷人人爽人人干人人爱| 亚洲avbb在线观看| 日本 av在线| 久久伊人香网站| 99久久精品热视频| 男女那种视频在线观看| 久久精品国产自在天天线| 特大巨黑吊av在线直播| avwww免费| 国产成人欧美在线观看| 午夜a级毛片| 亚洲国产精品久久男人天堂| 精品国产美女av久久久久小说| 亚洲国产精品久久男人天堂| 国产91精品成人一区二区三区| 国产又黄又爽又无遮挡在线| 长腿黑丝高跟| 亚洲国产欧美网| 嫩草影院精品99| 国产亚洲精品一区二区www| 少妇裸体淫交视频免费看高清| 一个人看视频在线观看www免费 | av在线蜜桃| 禁无遮挡网站| 网址你懂的国产日韩在线| 乱人视频在线观看| 国产亚洲av嫩草精品影院| 日韩有码中文字幕| 18+在线观看网站| 热99在线观看视频| 国产美女午夜福利| 免费无遮挡裸体视频| 日本免费a在线| 日韩国内少妇激情av| 欧美zozozo另类| 国语自产精品视频在线第100页| 一个人观看的视频www高清免费观看| 亚洲av二区三区四区| a在线观看视频网站| 999久久久精品免费观看国产| 久久久国产成人免费| 日韩免费av在线播放| 99热6这里只有精品| 久久久久亚洲av毛片大全| 国产亚洲欧美98| 婷婷精品国产亚洲av在线| 欧美丝袜亚洲另类 | 久久久久国产精品人妻aⅴ院| 夜夜爽天天搞| 久久精品91蜜桃| 真人一进一出gif抽搐免费| 日本免费一区二区三区高清不卡| 非洲黑人性xxxx精品又粗又长| 男女午夜视频在线观看| 亚洲人成网站在线播| 国产伦精品一区二区三区视频9 | 他把我摸到了高潮在线观看| 国产色婷婷99| 亚洲国产欧洲综合997久久,| www.熟女人妻精品国产| 国产精品香港三级国产av潘金莲| 欧美黑人欧美精品刺激| 久久精品91无色码中文字幕| 一级毛片女人18水好多| 成人三级黄色视频| 欧美xxxx黑人xx丫x性爽| 国产色婷婷99| eeuss影院久久| 禁无遮挡网站| 少妇的逼水好多| 久久亚洲精品不卡| 免费看日本二区| 人人妻,人人澡人人爽秒播| 亚洲国产色片| 亚洲人成网站在线播放欧美日韩| 性色av乱码一区二区三区2| 亚洲国产欧美网| 有码 亚洲区| 男女那种视频在线观看| 国产精品免费一区二区三区在线| 老司机福利观看| 精品久久久久久,| 亚洲成a人片在线一区二区| 亚洲专区中文字幕在线| 午夜激情欧美在线| 欧美bdsm另类| 亚洲欧美日韩高清在线视频| 亚洲国产欧美人成| www日本在线高清视频| 最后的刺客免费高清国语| 在线a可以看的网站| 色精品久久人妻99蜜桃| 激情在线观看视频在线高清| 一个人免费在线观看电影| 香蕉久久夜色| or卡值多少钱| 日韩欧美在线乱码| 老汉色av国产亚洲站长工具| 国产成人啪精品午夜网站| 亚洲成人中文字幕在线播放| 欧美3d第一页| av中文乱码字幕在线| 少妇人妻精品综合一区二区 | 久久久国产成人免费| 三级国产精品欧美在线观看| 欧美乱妇无乱码| 国产精品久久久久久久电影 | 亚洲国产高清在线一区二区三| 久久天躁狠狠躁夜夜2o2o| 一个人看的www免费观看视频| 18禁美女被吸乳视频| 国产探花极品一区二区| 观看美女的网站| 香蕉久久夜色| 国产精品自产拍在线观看55亚洲| 免费一级毛片在线播放高清视频| 亚洲熟妇熟女久久| 国产精品一区二区免费欧美| 黄片小视频在线播放| 观看免费一级毛片| 国产精品免费一区二区三区在线| 一本精品99久久精品77| 色综合亚洲欧美另类图片| 亚洲天堂国产精品一区在线| 香蕉av资源在线| 制服人妻中文乱码| 久久久久久久久久黄片| 亚洲男人的天堂狠狠| 五月伊人婷婷丁香| 精品午夜福利视频在线观看一区| 亚洲欧美日韩无卡精品| 51国产日韩欧美| 婷婷精品国产亚洲av在线| aaaaa片日本免费| 日韩成人在线观看一区二区三区| 欧美色欧美亚洲另类二区| 久久久久久九九精品二区国产| 亚洲国产色片| 一进一出抽搐gif免费好疼| 久9热在线精品视频| 麻豆久久精品国产亚洲av| 亚洲成av人片免费观看| 亚洲成人免费电影在线观看| 黄片大片在线免费观看| 九九在线视频观看精品| 嫁个100分男人电影在线观看| 在线观看免费视频日本深夜| 国产成人影院久久av| 亚洲国产高清在线一区二区三| 国产亚洲av嫩草精品影院| 国产99白浆流出| 亚洲国产色片| 伊人久久大香线蕉亚洲五| 亚洲人与动物交配视频| 国内精品美女久久久久久| 国产精品美女特级片免费视频播放器| 亚洲av成人不卡在线观看播放网| 日本黄色片子视频| 久久久久久九九精品二区国产| 三级毛片av免费| 国产三级在线视频| 又爽又黄无遮挡网站| 亚洲欧美日韩无卡精品| 国产69精品久久久久777片| 美女大奶头视频| www国产在线视频色| 国产不卡一卡二| 国内久久婷婷六月综合欲色啪| 欧美午夜高清在线| 99热这里只有是精品50| 亚洲熟妇熟女久久| 少妇的逼水好多| 欧美性猛交╳xxx乱大交人| 老司机深夜福利视频在线观看| 男女之事视频高清在线观看| 岛国在线观看网站| 中国美女看黄片| 精品一区二区三区视频在线观看免费| 亚洲欧美日韩高清专用| 成人av一区二区三区在线看| 午夜激情欧美在线| 亚洲欧美日韩无卡精品| 精品久久久久久成人av| 丰满的人妻完整版| 精品99又大又爽又粗少妇毛片 | 特级一级黄色大片| 女人十人毛片免费观看3o分钟| 亚洲精品一卡2卡三卡4卡5卡| 男人的好看免费观看在线视频| 三级国产精品欧美在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲成人精品中文字幕电影| 校园春色视频在线观看| 欧美午夜高清在线| 国产伦精品一区二区三区四那| 18禁在线播放成人免费| 色老头精品视频在线观看| 午夜老司机福利剧场| 3wmmmm亚洲av在线观看| 国产97色在线日韩免费| av国产免费在线观看| 国产久久久一区二区三区| 高清日韩中文字幕在线| 国内精品久久久久精免费| 亚洲av成人av| 午夜日韩欧美国产| 淫秽高清视频在线观看| 18+在线观看网站| 国产一区二区三区在线臀色熟女| 亚洲一区二区三区不卡视频| 亚洲欧美一区二区三区黑人| 色老头精品视频在线观看| 欧美黑人巨大hd| 三级国产精品欧美在线观看| 国产中年淑女户外野战色| 五月伊人婷婷丁香| 欧美大码av| 欧美日韩中文字幕国产精品一区二区三区| 亚洲无线观看免费| 日韩av在线大香蕉| 亚洲中文字幕日韩| 黄色片一级片一级黄色片| 欧美av亚洲av综合av国产av| www.熟女人妻精品国产| 亚洲18禁久久av| 亚洲午夜理论影院| 亚洲国产精品成人综合色| 国产精品久久久久久久电影 | 国产黄片美女视频| 午夜精品在线福利| 欧美又色又爽又黄视频| 精品午夜福利视频在线观看一区| 国产91精品成人一区二区三区| 丁香六月欧美| 国产黄色小视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 国产欧美日韩一区二区精品| 久久人人精品亚洲av| 国产精品亚洲av一区麻豆| 他把我摸到了高潮在线观看| 国产伦精品一区二区三区四那| 两个人视频免费观看高清| 国产午夜福利久久久久久| 欧美黄色淫秽网站| 91麻豆精品激情在线观看国产| 午夜视频国产福利| 免费人成视频x8x8入口观看|