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

    Cu、Ni摻雜FeS2電子結構與光學性質(zhì)的第一性原理計算

    2017-04-26 09:30:15武佳佳馬萬坤覃文慶
    中國有色金屬學報 2017年3期

    武佳佳,馬萬坤,焦 芬,覃文慶

    ?

    Cu、Ni摻雜FeS2電子結構與光學性質(zhì)的第一性原理計算

    武佳佳,馬萬坤,焦 芬,覃文慶

    (中南大學資源加工與生物工程學院,長沙 410083)

    采用基于密度泛函理論框架下的第一性原理平面波超軟贗勢法,研究Cu、Ni單摻雜與共摻雜體系中FeS2的晶體結構、電子結構和光學性質(zhì)。結果表明:摻雜后黃鐵礦發(fā)生晶格畸變,晶格常數(shù)變大,摻雜在FeS2禁帶中引入雜質(zhì)能級,使禁帶變窄,費米能級上移進入導帶,摻雜黃鐵礦的態(tài)密度穿過費米能級,形成簡并半導體,體系的導電率增強。光學性質(zhì)計算表明:摻雜后介電函數(shù)虛部主峰、吸收系數(shù)和光電導率均出現(xiàn)紅移、峰值減小。共摻雜后的光躍遷強度明顯增強,可見光區(qū)范圍內(nèi)的光吸收系數(shù)和光電導率均增大,說明Cu-Ni共摻雜顯著增強FeS2對光的吸收以及光電轉換效率。

    Cu-Ni摻雜;FeS2;光學性質(zhì);第一性原理計算

    隨著能源需求的激增和環(huán)境污染的加劇,長期以來,太陽能電池備受關注并發(fā)展迅速。但是發(fā)展至今,其年產(chǎn)量與世界能源需求量相比還相差甚遠,制約其發(fā)展的關鍵問題是高昂的制造成本[1?3]。金屬硫化物因其半導體屬性在太陽能電池、熱電材料以及存儲裝置等領域顯示出了廣闊的應用前景,因而得到了廣泛 關注。

    黃鐵礦(FeS2)是地殼中最常見的硫化礦,組成元素無毒且地球儲量豐富。立方晶系黃鐵礦具有優(yōu)異的半導體特性[4?5],并且其環(huán)境相容性好,制備成本低廉,是一種較有研究價值的新型太陽能電池材料[6?10]。已有研究[11]表明:采用人工合成方法制備出FeS2薄膜,盡管其易于氧化降解,但仍表現(xiàn)出良好的光電性能。

    摻雜是調(diào)制材料的電子結構并改變其光電性能的常用手段,對光電材料摻雜改性的理論研究具有重大意義[12]。對于黃鐵礦來說,天然黃鐵礦物中存在不同類型的雜質(zhì)和含量,李玉瓊等[13?16]的研究結果表明,鈷、鎳、砷的摻雜均使得黃鐵礦的吸收帶邊發(fā)生明顯的紅移。FERRER等[17]制備了摻Cu黃鐵礦薄膜,發(fā)現(xiàn)摻雜后電阻率和摻雜體系的光吸收邊均顯著增大。LEHNER等[18]使用化學氣相沉積制備了Ni摻雜的黃鐵礦薄膜,證明Ni的摻入在光電導性測量中會顯著影響其光電導的敏感性。龍飛等[19]研究Ni摻雜FeS2的可見光催化活性,發(fā)現(xiàn)Ni摻雜有利于提高Fe1?yNiS2的可見光催化活性,當Ni的摻入量為0.125時,其光催化活性最好。盡管對黃鐵礦光電性能、Ni摻雜[20]的研究已有報道,但用第一性原理研究Cu、Ni共摻雜對黃鐵礦電子結構及光學性質(zhì)的影響尚未見報道。

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

    1 理論模型和計算方法

    半導體化合物FeS2屬于等軸晶系,其空間群為,每個晶胞含有4個FeS2分子單元,F(xiàn)e原子分布在晶胞的6個面心和8個頂角上,每個Fe原子與6個相鄰的硫原子配位,形成空間八面體構造,而每個S原子與3個Fe原子和1個S原子配位,形成四面體構造,實驗測得其晶格常數(shù)為===0.5379 nm,===90°。計算采用(2×2×2)的FeS2超晶胞,晶胞中包含32個FeS2分子單元共96個原子。摻雜時分別用Cu和Ni原子取代Fe原子建立Fe29Cu3S64和Fe29Ni3S64超晶胞,用Cu和Ni同時取代Fe建立Fe26Cu3Ni3S64超晶胞。摻雜后形成Fe1?xCuS2、Fe1?yNiS2和Fe1?x?yCuNiS2(、分別代表Cu和Ni在FeS2超晶胞中的原子含量),其中==0.094。圖1所示為FeS2單胞及Cu、Ni取代Fe原子的(2×2×2) FeS2超晶胞模型。

    本計算中采用基于密度泛函理論(DFT)框架下的第一性原理贗勢平面波法,主要的計算工作由CASTEP軟件包[21]完成。計算中采用超軟贗勢(USPP)[22]來處理離子實與電子之間的相互作用,交換關聯(lián)泛函采用廣義梯度近似GGA下的PW91梯度修正近似。平面波截斷能測試表明,截斷能取350 eV較為合理。采用BFGS算法對體系進行幾何優(yōu)化,得到穩(wěn)定的結構,優(yōu)化的收斂標準為:原子位移的收斂閾值為0.02 pm,原子間作用力的收斂閾值為0.8 eV/nm,原子間的內(nèi)應力收斂標準為0.1 GPa,最大能量改變的收斂標準為2.0×10?5eV/atom,自洽計算收斂精度設置為每個原子1×10?3eV。為保證晶胞體系能量和構型在準完備平面波基水平上的收斂,布里淵區(qū)的積分計算采用4×4×2的Monkhorst Pack[23]特殊點。參與計算的價態(tài)電子為Fe 3d64s2、S 3s23p4、Cu 3d104s2和Ni 3d84s2。

    2 結果與討論

    2.1 幾何結構優(yōu)化結果

    Cu、Ni摻雜FeS2前后的超晶胞模型經(jīng)幾何優(yōu)化后,所得晶格常數(shù)及總能量見表1。從表1可以看出,計算結果與實驗和其他計算結果符合得很好,平衡晶格常數(shù)比實驗值小0.68%。實驗誤差在1%左右,說明計算方法是非??煽康?,可以進行其他性質(zhì)的計算。另外,摻雜后的晶胞體積相比原FeS2晶體稍有增大,原因是摻雜原子半徑大于Fe原子半徑一定程度上破壞了晶格周期性,造成晶格畸變,導致鍵長變長。晶胞優(yōu)化的結果與文獻[24]的結果一致。

    圖1 FeS2單胞及Cu、Ni取代Fe原子的(2×2×2)FeS2超晶胞結構

    表1 幾何優(yōu)化后FeS2和Cu、Ni摻雜后FeS2的晶胞參數(shù)和總能

    2.2 電子結構

    2.2.1 電子態(tài)密度

    為了分析摻雜對FeS2電子結構的影響,先對FeS2的幾何結構進行了優(yōu)化并在此基礎上計算了摻雜Ni、Cu后FeS2的態(tài)密度圖和能帶結構。圖2所示為摻雜前后FeS2的總態(tài)密度以及分態(tài)密度。從圖2(a)可以看出,本征態(tài)FeS2的價帶由3部分組成:?8.19~0 eV的上價帶主要由非成鍵Fe 3d態(tài)和少量S的3p態(tài)耦合而成;其次,?13.3~?10.1 eV和?17.5~?13.5 eV范圍內(nèi)的能帶幾乎全部由S 3s態(tài)貢獻,僅有少量的Fe 4s和S 3p態(tài)貢獻;費米能級附近價帶區(qū)的態(tài)密度,主要由S的3p態(tài)和非成鍵Fe 3d態(tài)組成,且大部分由Fe 3d態(tài)貢獻。而對于導帶部分,其主要由S 3p和Fe的反鍵3d eg*態(tài)貢獻,S 3s態(tài)貢獻很小。從Fe和S的分態(tài)密度可以看出,S 3s態(tài)對深部價帶的貢獻較大,S 3p態(tài)則主要對淺部價帶產(chǎn)生貢獻。Fe 4s態(tài)對整個能帶范圍態(tài)密度的貢獻相對較少,而費米能級附近的態(tài)密度主要來自Fe 3d的貢獻,這和FeS2的總態(tài)密度相符合。當Cu取代Fe原子后(見圖2(b)),費米能級附近出現(xiàn)了雜質(zhì)能級,這些雜質(zhì)能級是由Cu的3d態(tài)所形成。摻雜后導帶整體向左移動,導帶寬度由4.15 eV減小到3.16 eV,其局域性增強;費米能級與總態(tài)密度的截距明顯減小,費米能級的這種對導帶的相對運動增強了FeS2的金屬性質(zhì)。當Ni摻雜后(見圖2(c)),由于Ni和Cu原子核外電子排布的相似性,兩者單獨摻雜對FeS2的電子態(tài)密度影響基本相同,但是Ni摻雜Fe的3d軌道態(tài)密度較Cu摻雜向深部價帶有一個較大的偏移,而價帶頂和導帶底仍主要由Fe 3d態(tài)決定。圖2(d)所示為Cu-Ni共摻雜時的態(tài)密度圖,可以看到費米能級附近依然存在雜質(zhì)能帶而強度有所不同。費米能級附近的態(tài)密度主要還是由Fe 3d和S 3s貢獻,Ni和Cu的d軌道電子態(tài)密度對FeS2晶體的價帶與導帶的貢獻很小,并且只在?3.18 eV~?1.05 eV的淺部價帶處有一定作用。Cu-Ni共摻雜后FeS2電子態(tài)發(fā)生了簡并,費米能級上移進入導帶,價帶部分離域性增強。這表明Cu-Ni摻雜后FeS2體系呈現(xiàn)金屬性,Ni摻雜比Cu摻雜表現(xiàn)出更強的金屬性,共摻時金屬性最強。

    2.2.2 能帶結構

    圖3所示為FeS2摻雜前后的能帶結構。圖3中、、、、為FeS2晶體第一布里淵區(qū)高對稱點。為進行對比,將理想FeS2的第一布里淵區(qū)沿對稱方向的點取樣設置為與摻雜后的體系一致。在計算中所考慮的禁帶寬度是從價帶最高點到導帶最低點之間的距離,從圖3(a)中本征態(tài)FeS2的能帶圖可以看出,價帶頂和導帶底都位于點處,說明FeS2是直接帶隙半導體,帶隙為0.55 eV,小于實驗值0.95 eV,這是局域密度泛函的結果。一般認為與局域密度泛函理論中Kohn-Sham方程的本征值不能給出系統(tǒng)的激發(fā)態(tài)能量有較大關系,使位于導帶的電子態(tài)能量值比實驗值偏小,從而帶隙偏小[25],一般能帶值與理論值相差30%~50%。對于本研究中本征FeS2結構而言,由于Fe 3d態(tài)的能量被過高估計,使得其與S 3p態(tài)之間的相互作用增強,價帶帶寬增大,因此帶隙偏低。但理論計算與實驗之間帶隙的偏差并不影響對FeS2電子結構及相關性質(zhì)的理論分析。摻雜后,3種結構的導帶部分圖譜比較相似,與未摻雜FeS2相比,導帶部分變化較為平緩,摻雜后FeS2的費米能級F上移進入導帶內(nèi),這是由于雜質(zhì)元素的摻入使得導電載流子數(shù)目增多,導帶底附近的量子態(tài)已基本被電子所占據(jù),發(fā)生載流子的簡并化,形成簡并半導體。從圖3(b)~(d)可以看出,Cu、Ni單獨摻雜和共摻雜后,F(xiàn)eS2的帶隙分別為0.44、0.61和0.35 eV,這表明Ni摻雜后FeS2的禁帶寬度增大而其他兩種摻雜禁帶寬度減小。對比圖4中各圖可以發(fā)現(xiàn),摻雜后FeS2導帶帶寬均減小而且均發(fā)生了下移,禁帶寬度均變小,其中Cu-Ni共摻雜的變化最顯著,這也與態(tài)密度中分析的能帶變化的結果相一致。

    圖2 FeS2和Cu、Ni摻雜后FeS2的態(tài)密度

    圖3 FeS2和Cu、Ni摻雜后FeS2的能帶結構圖

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

    為了研究Cu、Ni摻雜FeS2的光學躍遷特性,必須對其光學性質(zhì)進行研究。本研究中利用第一性原理對摻雜體系進行了光譜計算和分析。在線性響應范圍內(nèi)固體宏觀光學響應函數(shù)通??梢杂晒獾膹徒殡姾瘮?shù)()=1()+i2()或折射率()=()+i()來描述,體系在較小波矢下對光場的線性響應由介電函數(shù)的虛部2()決定,它可以通過占據(jù)態(tài)和未占據(jù)態(tài)的動量躍遷矩陣元計算得到。根據(jù)Kramers-Kronig色散關系和電子躍遷概率的定義可推導出晶體的介電函數(shù)、吸收系數(shù)。介電函數(shù)實部1()和虛部2()的計算公式如下所示[26]:

    (1)

    (2)

    2.3.1 復介電函數(shù)

    介電常數(shù)的虛部2()對任何材料來說都是非常重要的,它作為溝通帶間躍遷微觀物理過程與固體電子結構的橋梁,主要反映了電子在能級間的躍遷機理。圖4所示為FeS2摻雜前后的介電函數(shù)虛部與能量的關系。

    圖4 Cu、Ni摻雜FeS2前后的介電函數(shù)虛部與能量的關系

    本征態(tài)介電函數(shù)虛部譜在=2.57 eV附近出現(xiàn)波峰,其主要是由價帶頂?shù)腟 2p態(tài)和導帶底的Fe 3d態(tài)電子躍遷引起的。與本征態(tài)相比,摻雜后峰值強度都有所降低,共摻的峰值降幅較大,且峰譜向低能方向發(fā)生了明顯的紅移,在0.22 eV附近出現(xiàn)了新的介電峰,這主要是Cu、Ni摻入的雜質(zhì)能級與相應的價帶頂和導帶底之間的電子躍遷引起。Cu-Ni共摻雜之后,譜線主峰明顯紅移,在=0.3 eV處出現(xiàn)了一個明顯的小波峰,這說明相比本征態(tài)共摻雜后由于雜質(zhì)能級的作用,體系中導帶和價帶之間的電子躍遷呈現(xiàn)臺階特性,并促使更多深部價帶的電子躍遷到導帶。在<2 eV的可見光區(qū),Cu-Ni共摻雜FeS2的光躍遷強度遠大于未摻雜FeS2的,并且共摻雜的光躍遷強度明顯高于Cu、Ni單摻雜的光躍遷強度,說明Cu-Ni共摻雜能夠很好的改善FeS2電子在<2eV的低能量范圍的光學躍遷特性。

    2.3.2吸收系數(shù)

    圖5 Cu、Ni摻雜FeS2前后的吸收系數(shù)與能量的關系

    2.3.4 光電導率

    光電導率是光電子材料的一個重要參數(shù),它描述的是半導體材料電導率隨著光照能量變化的物理參量。在它與介電函數(shù)的關系中其實部1和介電函數(shù)虛部2互相對應,也能夠反應光吸收和能隙的情況。

    圖6所示為FeS2摻雜前后光電導率隨能量變化的關系。由圖6可知,F(xiàn)eS2的光電導率與介電函數(shù)的虛部是對應的,1()總體呈先上升后下降的變化趨勢。本征態(tài)FeS2在<10 eV的能量范圍,光電導率實部1有2個峰值點,各峰值點對應的能量分別為2.93和5.23 eV。進行Cu、Ni摻雜后,這兩個峰值位置出現(xiàn)小幅紅移,并且Cu-Ni共摻雜體系中出現(xiàn)了位于0.63 eV處的峰值。在<2 eV的可見光區(qū)范圍內(nèi),摻雜體系的光電導率均大于本征態(tài)的光電導率,綜合來看,Cu、Ni共摻雜體系的光電導率明顯優(yōu)于單摻雜體系的光電導率,摻雜后材料表現(xiàn)出顯著的半導體特性,有利于成為電子或空穴的俘獲中心,減少電子和空穴的復合,提高量子產(chǎn)量,表現(xiàn)出明顯的協(xié)同效應。

    圖6 FeS2摻雜前后的光電導率與能量的關系

    3 結論

    1) 對Cu、Ni摻雜FeS2的幾何結構優(yōu)化表明,單摻雜和共摻雜均使FeS2的晶格常數(shù)增大,摻雜FeS2發(fā)生晶格畸變,鍵長變長。

    2) 本征態(tài)FeS2的電子結構優(yōu)化表明,其禁帶寬度為0.55 eV,費米能級附近導帶區(qū)0~4 eV范圍內(nèi)的態(tài)密度主要由S 3p和Fe的反鍵3d eg*態(tài)貢獻,而S 3s態(tài)貢獻很小。Cu、Ni摻雜使FeS2的電子態(tài)發(fā)生了簡并,費米能級上移進入導帶,價帶部分離域性增強,并且在?3.18~?1.05 eV的淺部價帶產(chǎn)生了由Cu和Ni的3d態(tài)貢獻的雜質(zhì)能級。

    3) 光學性質(zhì)計算結果表明:單元素摻雜后,F(xiàn)eS2的介電函數(shù)虛部、光學吸收譜以及光電導率各主峰均出現(xiàn)小幅紅移,而共摻后介電函數(shù)譜較Cu、Ni單摻雜均出現(xiàn)顯著紅移,價帶頂和導帶之間的光躍遷強度明顯增強,可見光區(qū)范圍內(nèi)的光吸收系數(shù)和光電導率均增大,共摻雜顯著增強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] JAGANNATH PANIGRAHI, VANDANA, RAJBIR SINGH, NEHA BATRA, JHUMA GOPE, MUKUL SHARMA, PATHIB P, SRIVASTAVA S K, RAUTHAN C M S, SINGH P K. Impedance spectroscopy of crystalline silicon solar cell: Observation of negative capacitance[J]. Solar Energy, 2016, 136: 412?420.

    [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: 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 S. 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] LIU S T, WU J, YU P, DING Q H, ZHOU Z H, LI H D, LAI C C, CHUEH Y L, WANG Z M. Phase-pure iron pyrite nanocrystals for low-cost photodetectors[J]. Nanoscale Research Letters, 2014, 9: 549?556.

    [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, 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.

    [11] 張 輝, 張仁剛, 萬冬云, 王寶義, 魏 龍. 摻雜對FeS2薄膜光電性能的影響[J]. 太陽能學報, 2006, 27(5): 423?427. ZHANG Hui, ZHANG Ren-gang, WAN Dong-yun, WANG Bao-yi, WEI Long. Effects of doping impurities in pyrite films on the optical and electrical properties[J]. Acta Energiae Solaris Sinica, 2006, 27(5): 423?427.

    [12] 荊 濤, 張 蘋, 闞 偉, 田景芝, 鄧啟剛. N-Al共摻雜TiO2電子結構及光學性質(zhì)的理論研究[J]. 中國有色金屬學報, 2015, 25(4): 1018?1024. JING Tao, ZHANG Ping, KAN Wei, TIAN Jing-zhi, DENG Qi-gang. Theory studies on electronic structure and optical properties of N-Al co-doped anatase TiO2[J]. The Chinese Journal of Nonferrous Metals, 2015, 25(4): 1018?1024.

    [13] 李玉瓊, 陳建華, 陳 曄, 郭 進. 黃鐵礦(100)表面性質(zhì)的密度泛函理論計算及其對浮選的影響[J]. 中國有色金屬學報, 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.

    [14] 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.

    [15] 陳建華, 鐘建蓮, 李玉瓊, 陳 曄, 郭 進. 黃鐵礦、白鐵礦和磁黃鐵礦的電子結構及可浮性[J]. 中國有色金屬學報, 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.

    [16] 李玉瓊, 陳建華, 郭 進. 天然雜質(zhì)對黃鐵礦的電子結構及催化活性的影響[J]. 物理學報, 2011, 60(9): 097801-1?097801-8. LI Yu-qiong, CHEN Jian-hua, GUO Jin. Influence of natural impurity on electronic structure and catalytic activity of pyrite[J]. Acta Phys Sin, 2011, 60(9): 097801-1?097801-8.

    [17] 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.

    [18] 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.

    [19] 龍 飛, 張 勁, 張明月, 何金云, 武曉鸝, 鄒正光. Ni摻雜黃鐵礦型FeS2的溶劑熱合成及其可見光催化活性[J]. 無機化學學報, 2015, 31(6): 1119?1124.LONG Fei, ZHANG Jin, ZHANG Ming-yue, HE Jin-yun, WU Xiao-li, ZOU Zheng-guang. Ni-doped FeS2: Solvothermal synthesis and the visible-light photocatalytic properties[J]. Chinese Journal of Inorganic Chemistry, 2015, 31(6): 1119?1124.

    [20] 熊志華, 饒建平, 江風益. CdS摻Mg和Ni電子結構和光學性質(zhì)的密度泛函理論研究[J]. 光學學報, 2007, 27(12): 2225?2228.XIONG Zhi-hua, RAO Jian-ping, JIANG Feng-yi. Density functional calculations of electronic structure and optical properties on Mg and Ni-doped CdS[J]. Acta Optical Sinica, 2007, 27(12): 2225?2228.

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

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

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

    [24] 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.

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

    [26] 沈?qū)W礎. 半導體光學性質(zhì)[M]. 北京: 科學出版社, 1992: 24. SHEN Xue-chu. The optical properties of semiconductor[M]. Beijing: Science Press, 1992: 24.

    (編輯 龍懷中)

    First principle calculation of electronic structures and optical properties of copper and nickel doped FeS2

    WU Jia-jia, MA Wan-kun, JIAO Fen, QIN Wen-qing

    (School of Minerals Processing and Bioengineering, Central South University, Changsha 410083, China)

    The crystal and electronic structures, optical properties of copper and nickel mono- and co-doped FeS2were studied by the ?rst-principles plane-wave ultrasoft pseudopotential method with the generalized gradient approximation. The results show that, the lattice deformation of FeS2occurs and the lattice parameters increase after doping. And the impurity energy level is introduced by Cu and Ni doping, which also leads to decrease of the band gap of FeS2. Fermi levels moves upward and into conduction band and density of states traverses the Fermi level, indicating that the doped FeS2is degenerate semiconductor which shows stronger metallic characteristic and electrical conductivity. In addition, the doping results in peaks redshift and amplitudes decrease of imaginary part of dielectric function, absorption coefficient and conductivity. Optical transition is significantly enhanced and so are absorption coefficient and conductivity within visible region in the co-doping system. The photoabsorption capacity and photoelectrical efficiency of FeS2are highly strengthened in Cu-Ni co-doping system.

    Cu-Ni doping; FeS2; optical property; ?rst-principles calculation

    Project (51274255) supported by the National Natural Science Foundation of China; Project (2016RS2016) supported by Provincial Science and Technology Leader Innovation Team of Interface Chemistry of Efficient and Clean Utilization of Complex Mineral Resources, Hunan Province, China; Project (2015CX005) supported by Innovation Driven Plan of Central South University, China; Project supported by the Co-Innovation Centre for Clean and Efficient Utilization of Strategic Metal Mineral Resources, China

    2016-04-06; Accepted date: 2016-08-29

    JIAO Fen, Tel: +86-731-88830884; E-mail: jfen0601@126.com

    10.19476/j.ysxb.1004.0609.2017.03.019

    1004-0609(2017)-03-0605-08

    O474

    A

    國家自然科學基金資助項目(51274255);湖南省科技領軍人才(復雜礦產(chǎn)資源高效清潔利用的界面化學創(chuàng)新團隊)(2016RS2016);中南大學“創(chuàng)新驅(qū)動計劃”資助項目(2015CX005);中南大學戰(zhàn)略金屬礦產(chǎn)資源清潔高效利用協(xié)同創(chuàng)新中心資助項目

    2016-04-06;

    2016-08-29

    焦 芬,副教授,博士;電話:0731-88830884;E-mail:jfen0601@126.com

    男女高潮啪啪啪动态图| 黑人猛操日本美女一级片| 18禁裸乳无遮挡动漫免费视频| 少妇的逼好多水| 国产亚洲午夜精品一区二区久久| 国产精品久久久久久av不卡| 999精品在线视频| av黄色大香蕉| 免费人成在线观看视频色| 亚洲久久久国产精品| 免费av不卡在线播放| 男女边摸边吃奶| 久久精品久久久久久噜噜老黄| 黑人猛操日本美女一级片| 国产成人免费无遮挡视频| 久久人人爽人人片av| 香蕉精品网在线| 日韩电影二区| 精品人妻偷拍中文字幕| 少妇人妻久久综合中文| 成人漫画全彩无遮挡| 日本欧美国产在线视频| 18禁动态无遮挡网站| 一级毛片aaaaaa免费看小| 伊人久久精品亚洲午夜| 免费高清在线观看视频在线观看| 九色成人免费人妻av| 亚洲精品av麻豆狂野| 亚洲精品日本国产第一区| 亚洲精品久久午夜乱码| 亚洲精品乱久久久久久| 国产国语露脸激情在线看| 最新的欧美精品一区二区| 亚洲国产精品一区二区三区在线| 中国美白少妇内射xxxbb| 性色av一级| 熟女人妻精品中文字幕| 大香蕉97超碰在线| 日韩强制内射视频| 亚洲精品色激情综合| 亚洲av日韩在线播放| 日韩成人伦理影院| 精品一品国产午夜福利视频| 九九爱精品视频在线观看| 亚洲欧美日韩另类电影网站| 夜夜骑夜夜射夜夜干| 在线观看美女被高潮喷水网站| 精品一区在线观看国产| 一级毛片 在线播放| 日本av免费视频播放| 三级国产精品欧美在线观看| 国产精品.久久久| 97超碰精品成人国产| av在线老鸭窝| 新久久久久国产一级毛片| 97超碰精品成人国产| 日韩电影二区| 免费人成在线观看视频色| 国产精品熟女久久久久浪| xxxhd国产人妻xxx| 欧美少妇被猛烈插入视频| 国产精品一区二区在线不卡| 菩萨蛮人人尽说江南好唐韦庄| 草草在线视频免费看| 久久99一区二区三区| 夫妻性生交免费视频一级片| 丁香六月天网| 国产欧美日韩一区二区三区在线 | 国产视频内射| 一本一本综合久久| 亚洲精品国产色婷婷电影| 国产伦精品一区二区三区视频9| 中文欧美无线码| 日韩人妻高清精品专区| 永久网站在线| 国产精品.久久久| 高清视频免费观看一区二区| 免费人成在线观看视频色| 亚洲欧美色中文字幕在线| videos熟女内射| 在线观看人妻少妇| av线在线观看网站| 黑人巨大精品欧美一区二区蜜桃 | 国产一区二区在线观看av| 午夜福利网站1000一区二区三区| 亚洲精品美女久久av网站| 成年美女黄网站色视频大全免费 | 天天躁夜夜躁狠狠久久av| 熟女电影av网| 久久鲁丝午夜福利片| 亚洲无线观看免费| a级片在线免费高清观看视频| 两个人免费观看高清视频| 国产又色又爽无遮挡免| 国产成人精品福利久久| 国产精品一区二区在线观看99| 国产一区二区三区综合在线观看 | 久久精品国产亚洲av天美| 久久久久久久久久久丰满| 亚洲少妇的诱惑av| 99热6这里只有精品| 精品少妇久久久久久888优播| 黄色怎么调成土黄色| 啦啦啦中文免费视频观看日本| 在线观看国产h片| 国产高清国产精品国产三级| 久久久国产一区二区| 成人影院久久| 激情五月婷婷亚洲| 中国国产av一级| 国产亚洲欧美精品永久| 日韩av免费高清视频| 成人国产麻豆网| 国产伦精品一区二区三区视频9| 婷婷色综合大香蕉| 欧美bdsm另类| 亚洲精品乱码久久久久久按摩| 永久网站在线| 中文天堂在线官网| 日韩一本色道免费dvd| 亚洲精品亚洲一区二区| 国产成人精品福利久久| 伦精品一区二区三区| 欧美日韩综合久久久久久| 日本-黄色视频高清免费观看| 永久网站在线| a级毛片在线看网站| 视频区图区小说| 国国产精品蜜臀av免费| av不卡在线播放| 欧美97在线视频| 99精国产麻豆久久婷婷| 99热网站在线观看| 一级毛片我不卡| 国产在线一区二区三区精| 大香蕉久久网| 午夜日本视频在线| 一区二区三区四区激情视频| 亚洲av不卡在线观看| 热99久久久久精品小说推荐| 欧美精品亚洲一区二区| 欧美成人精品欧美一级黄| 精品亚洲乱码少妇综合久久| 五月开心婷婷网| 日本av免费视频播放| 国产av国产精品国产| 99热全是精品| 国产精品99久久99久久久不卡 | 亚洲情色 制服丝袜| 多毛熟女@视频| 国产成人a∨麻豆精品| 我要看黄色一级片免费的| 内地一区二区视频在线| 亚洲久久久国产精品| 永久免费av网站大全| 国产精品国产三级专区第一集| 中文字幕人妻丝袜制服| av网站免费在线观看视频| 狠狠婷婷综合久久久久久88av| 九九久久精品国产亚洲av麻豆| 国产亚洲最大av| 18禁动态无遮挡网站| 日本黄色日本黄色录像| 中国三级夫妇交换| 91久久精品国产一区二区成人| 最后的刺客免费高清国语| 久久久久网色| 国产乱来视频区| 国产精品蜜桃在线观看| 99久久人妻综合| 18禁动态无遮挡网站| 内地一区二区视频在线| 欧美亚洲日本最大视频资源| 中文字幕免费在线视频6| 久热久热在线精品观看| 国产一区有黄有色的免费视频| 免费黄频网站在线观看国产| a级片在线免费高清观看视频| 香蕉精品网在线| 最近中文字幕2019免费版| 波野结衣二区三区在线| 一区在线观看完整版| 3wmmmm亚洲av在线观看| av卡一久久| 成人二区视频| 欧美日韩亚洲高清精品| 亚洲第一av免费看| 黄色配什么色好看| 亚洲欧美精品自产自拍| 国产黄频视频在线观看| 亚洲av.av天堂| 欧美国产精品一级二级三级| 黄色怎么调成土黄色| 99九九在线精品视频| 日韩 亚洲 欧美在线| 熟妇人妻不卡中文字幕| 亚洲久久久国产精品| 人妻夜夜爽99麻豆av| 五月伊人婷婷丁香| 又粗又硬又长又爽又黄的视频| 久久ye,这里只有精品| 精品久久久精品久久久| 99久国产av精品国产电影| 99热6这里只有精品| 蜜桃国产av成人99| 精品99又大又爽又粗少妇毛片| 丁香六月天网| 日韩视频在线欧美| 亚洲在久久综合| 日日摸夜夜添夜夜添av毛片| 国产毛片在线视频| h视频一区二区三区| 免费看不卡的av| 免费高清在线观看日韩| 国产女主播在线喷水免费视频网站| 日本午夜av视频| 亚洲伊人久久精品综合| 国产精品欧美亚洲77777| 久久97久久精品| 欧美xxxx性猛交bbbb| 一个人免费看片子| 国产精品一区二区在线观看99| 国产成人精品福利久久| 亚洲欧美成人综合另类久久久| 亚洲高清免费不卡视频| 亚洲少妇的诱惑av| 久久99精品国语久久久| 亚洲精品国产av成人精品| 精品视频人人做人人爽| 久久国产精品大桥未久av| 一级毛片电影观看| a级毛色黄片| 欧美激情国产日韩精品一区| 日本猛色少妇xxxxx猛交久久| 麻豆乱淫一区二区| 高清视频免费观看一区二区| 七月丁香在线播放| 免费黄网站久久成人精品| 人人妻人人澡人人爽人人夜夜| 午夜av观看不卡| 男男h啪啪无遮挡| 久热这里只有精品99| 日本免费在线观看一区| 一边亲一边摸免费视频| 超碰97精品在线观看| 看十八女毛片水多多多| 国产国拍精品亚洲av在线观看| 王馨瑶露胸无遮挡在线观看| 亚洲精品国产av蜜桃| 另类亚洲欧美激情| av有码第一页| 亚洲精品久久久久久婷婷小说| 亚洲伊人久久精品综合| 亚洲欧洲精品一区二区精品久久久 | 日日摸夜夜添夜夜爱| 极品少妇高潮喷水抽搐| 精品国产露脸久久av麻豆| 亚洲综合色网址| 午夜激情福利司机影院| 日韩成人伦理影院| 婷婷色综合大香蕉| 搡老乐熟女国产| 午夜日本视频在线| 少妇被粗大猛烈的视频| 亚州av有码| 天堂8中文在线网| 国产老妇伦熟女老妇高清| 欧美xxxx性猛交bbbb| 久久亚洲国产成人精品v| 欧美日韩成人在线一区二区| 国产成人a∨麻豆精品| 天天影视国产精品| 国产精品成人在线| 九色亚洲精品在线播放| 日本欧美视频一区| 高清黄色对白视频在线免费看| 七月丁香在线播放| 热99国产精品久久久久久7| 免费看不卡的av| 亚洲精品乱久久久久久| 超碰97精品在线观看| 久久av网站| 2018国产大陆天天弄谢| 欧美 日韩 精品 国产| 91午夜精品亚洲一区二区三区| 久久99热这里只频精品6学生| 亚洲少妇的诱惑av| 国产免费一级a男人的天堂| 王馨瑶露胸无遮挡在线观看| 另类亚洲欧美激情| 国产精品一国产av| 午夜激情福利司机影院| 成人亚洲精品一区在线观看| 男女国产视频网站| av视频免费观看在线观看| 最黄视频免费看| 亚洲欧美清纯卡通| 亚洲av二区三区四区| 国国产精品蜜臀av免费| 伊人久久国产一区二区| 少妇人妻 视频| 久久久久精品性色| 女性生殖器流出的白浆| 一级二级三级毛片免费看| 涩涩av久久男人的天堂| 亚洲精品久久久久久婷婷小说| 在线 av 中文字幕| 日本黄色日本黄色录像| xxx大片免费视频| 日韩欧美精品免费久久| 最近2019中文字幕mv第一页| 在线观看一区二区三区激情| 女性生殖器流出的白浆| 亚洲美女黄色视频免费看| 97超视频在线观看视频| 最近2019中文字幕mv第一页| 日韩精品免费视频一区二区三区 | 国产综合精华液| 999精品在线视频| 高清视频免费观看一区二区| 91久久精品电影网| 男女边吃奶边做爰视频| 国产精品嫩草影院av在线观看| 日韩视频在线欧美| 99国产综合亚洲精品| 日本vs欧美在线观看视频| 王馨瑶露胸无遮挡在线观看| 亚洲人成网站在线播| 九草在线视频观看| 少妇丰满av| 热99国产精品久久久久久7| 成人亚洲欧美一区二区av| 国产高清国产精品国产三级| 国产日韩欧美在线精品| 亚洲av男天堂| 在线亚洲精品国产二区图片欧美 | 欧美人与性动交α欧美精品济南到 | 久久毛片免费看一区二区三区| 国产欧美日韩综合在线一区二区| 韩国av在线不卡| 男女边摸边吃奶| 一区二区三区免费毛片| 91精品三级在线观看| 精品久久久久久久久av| 精品人妻一区二区三区麻豆| 国产熟女午夜一区二区三区 | 国产有黄有色有爽视频| 欧美日韩视频高清一区二区三区二| 亚洲欧美色中文字幕在线| 色婷婷av一区二区三区视频| 国产精品一二三区在线看| 国产av精品麻豆| 黄片播放在线免费| 精品久久久噜噜| 国产欧美另类精品又又久久亚洲欧美| 国产69精品久久久久777片| 国产精品一区二区在线观看99| 男女高潮啪啪啪动态图| 美女视频免费永久观看网站| 最近最新中文字幕免费大全7| 国产精品一区二区在线不卡| 中文乱码字字幕精品一区二区三区| h视频一区二区三区| 日本午夜av视频| 久久精品国产亚洲网站| 国产精品.久久久| 色吧在线观看| 水蜜桃什么品种好| 久久毛片免费看一区二区三区| 亚洲国产精品一区二区三区在线| 夫妻午夜视频| 国产亚洲av片在线观看秒播厂| 欧美日韩成人在线一区二区| av不卡在线播放| 18+在线观看网站| 亚洲怡红院男人天堂| 一边摸一边做爽爽视频免费| 国产免费一级a男人的天堂| 久久青草综合色| 九九爱精品视频在线观看| 成人毛片60女人毛片免费| 青青草视频在线视频观看| 在线观看一区二区三区激情| 街头女战士在线观看网站| 午夜老司机福利剧场| 日韩成人伦理影院| 中文精品一卡2卡3卡4更新| 在线免费观看不下载黄p国产| 肉色欧美久久久久久久蜜桃| 日本欧美国产在线视频| 交换朋友夫妻互换小说| 51国产日韩欧美| 人体艺术视频欧美日本| 99久国产av精品国产电影| 熟女电影av网| 欧美日韩视频高清一区二区三区二| 国产亚洲av片在线观看秒播厂| 国产精品秋霞免费鲁丝片| 亚洲精品中文字幕在线视频| 汤姆久久久久久久影院中文字幕| 久久久久国产精品人妻一区二区| 成人毛片a级毛片在线播放| 大片免费播放器 马上看| 在线播放无遮挡| 中文精品一卡2卡3卡4更新| 高清黄色对白视频在线免费看| 日韩免费高清中文字幕av| 亚洲人成77777在线视频| 久久午夜综合久久蜜桃| 乱人伦中国视频| 欧美日韩精品成人综合77777| 免费观看性生交大片5| 国产探花极品一区二区| 青青草视频在线视频观看| 天天影视国产精品| 国产亚洲午夜精品一区二区久久| 久久久久久久亚洲中文字幕| 国产黄频视频在线观看| 国产日韩欧美视频二区| 人妻夜夜爽99麻豆av| 久久久久精品性色| 大话2 男鬼变身卡| 国产 精品1| 男女免费视频国产| 国产亚洲精品久久久com| 男女免费视频国产| 2021少妇久久久久久久久久久| 国产成人精品在线电影| 一本大道久久a久久精品| 欧美+日韩+精品| 黑人猛操日本美女一级片| 久久精品国产亚洲av天美| 国产精品人妻久久久影院| 最近的中文字幕免费完整| 久久久久国产网址| 亚洲国产精品一区三区| 久久精品国产亚洲av天美| 亚洲久久久国产精品| 亚洲精品色激情综合| 国产精品99久久99久久久不卡 | 亚洲精品色激情综合| 一级黄片播放器| 伊人久久精品亚洲午夜| 桃花免费在线播放| 国产免费视频播放在线视频| av在线老鸭窝| 如日韩欧美国产精品一区二区三区 | av播播在线观看一区| 亚洲精品,欧美精品| 成人国产麻豆网| 在线观看美女被高潮喷水网站| 一边摸一边做爽爽视频免费| av在线app专区| 母亲3免费完整高清在线观看 | 国产免费现黄频在线看| 成人影院久久| 国产色婷婷99| 777米奇影视久久| 免费少妇av软件| 最近最新中文字幕免费大全7| 国产在线一区二区三区精| 一级毛片 在线播放| av专区在线播放| 国产综合精华液| 国产成人精品在线电影| av女优亚洲男人天堂| 2022亚洲国产成人精品| 在线亚洲精品国产二区图片欧美 | 亚洲欧美成人精品一区二区| 国精品久久久久久国模美| 欧美日韩一区二区视频在线观看视频在线| 人妻少妇偷人精品九色| 国产午夜精品久久久久久一区二区三区| 性色av一级| 国产成人精品无人区| av电影中文网址| 欧美激情国产日韩精品一区| 久久99精品国语久久久| 亚洲,一卡二卡三卡| 夜夜爽夜夜爽视频| 欧美人与性动交α欧美精品济南到 | 九九在线视频观看精品| av天堂久久9| 成人二区视频| av播播在线观看一区| 能在线免费看毛片的网站| 亚洲人成网站在线观看播放| 日韩不卡一区二区三区视频在线| 五月天丁香电影| 在线播放无遮挡| 桃花免费在线播放| 天堂俺去俺来也www色官网| 最近中文字幕2019免费版| 国产精品一区二区在线不卡| 蜜桃久久精品国产亚洲av| 天堂中文最新版在线下载| 好男人视频免费观看在线| h视频一区二区三区| 国产一区二区在线观看av| 麻豆精品久久久久久蜜桃| 免费av不卡在线播放| 91久久精品国产一区二区三区| 亚洲美女搞黄在线观看| 美女主播在线视频| 一本久久精品| 亚洲,欧美,日韩| 精品久久蜜臀av无| 如何舔出高潮| 国产成人精品无人区| 久久久国产一区二区| tube8黄色片| 五月伊人婷婷丁香| 老司机影院成人| 国产在线免费精品| 多毛熟女@视频| 欧美成人午夜免费资源| 22中文网久久字幕| 国产成人免费无遮挡视频| 尾随美女入室| av又黄又爽大尺度在线免费看| 亚洲国产欧美日韩在线播放| 国产成人精品久久久久久| 色视频在线一区二区三区| videos熟女内射| 亚洲熟女精品中文字幕| 在现免费观看毛片| 国产精品一区二区三区四区免费观看| 免费观看性生交大片5| 80岁老熟妇乱子伦牲交| 人人妻人人澡人人看| 午夜激情av网站| 亚洲av二区三区四区| 久久影院123| 国产免费又黄又爽又色| 婷婷色综合大香蕉| 精品熟女少妇av免费看| 狂野欧美激情性xxxx在线观看| 久久久久久久久久久丰满| 日产精品乱码卡一卡2卡三| 99久久精品国产国产毛片| 久久婷婷青草| 中文精品一卡2卡3卡4更新| 蜜桃久久精品国产亚洲av| 欧美亚洲 丝袜 人妻 在线| 一本大道久久a久久精品| 不卡视频在线观看欧美| 黑人猛操日本美女一级片| 国产日韩欧美视频二区| 国产精品.久久久| 成人漫画全彩无遮挡| 精品视频人人做人人爽| 七月丁香在线播放| 一本久久精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 中文字幕亚洲精品专区| 精品国产一区二区三区久久久樱花| 午夜激情久久久久久久| 亚洲人成网站在线播| 久久99精品国语久久久| 国产老妇伦熟女老妇高清| 国产白丝娇喘喷水9色精品| 亚洲av电影在线观看一区二区三区| 在线观看一区二区三区激情| 少妇人妻久久综合中文| 欧美 日韩 精品 国产| 国产成人精品久久久久久| 久久久久人妻精品一区果冻| 久久久精品区二区三区| 欧美精品一区二区大全| 最后的刺客免费高清国语| 久久鲁丝午夜福利片| 成人影院久久| 精品视频人人做人人爽| 国产欧美另类精品又又久久亚洲欧美| 国产深夜福利视频在线观看| 亚洲熟女精品中文字幕| av女优亚洲男人天堂| 人人妻人人澡人人爽人人夜夜| 国产日韩一区二区三区精品不卡 | 最近手机中文字幕大全| 欧美人与善性xxx| av一本久久久久| 日韩一区二区三区影片| 久久精品人人爽人人爽视色| 又粗又硬又长又爽又黄的视频| 精品一区二区三卡| 日韩大片免费观看网站| 乱人伦中国视频| 美女主播在线视频| 国产亚洲一区二区精品| 看十八女毛片水多多多| 少妇猛男粗大的猛烈进出视频| 免费人成在线观看视频色| 日韩av在线免费看完整版不卡| 国产精品久久久久久av不卡| 国产伦理片在线播放av一区| 桃花免费在线播放| 亚洲国产色片| 插逼视频在线观看| 国产精品久久久久久av不卡| 精品少妇黑人巨大在线播放| 99九九线精品视频在线观看视频| 免费大片黄手机在线观看| 最近手机中文字幕大全| videos熟女内射| 永久免费av网站大全| 久久人人爽人人片av| 五月伊人婷婷丁香| 91成人精品电影| 丰满乱子伦码专区| 国产在线免费精品| 日本欧美国产在线视频| 日产精品乱码卡一卡2卡三| 高清毛片免费看| 日韩av免费高清视频|