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

    單軸應(yīng)變下CaO和SrO電子結(jié)構(gòu)和光學(xué)性質(zhì)的第一性原理研究

    2016-09-06 01:16:22邢美靜姚倩倩
    關(guān)鍵詞:導(dǎo)帶能帶單軸

    邢美靜,梁 凱,姚倩倩,李 實(shí),申 潔,趙 輝

    (天津師范大學(xué)物理與材料科學(xué)學(xué)院,天津300387)

    單軸應(yīng)變下CaO和SrO電子結(jié)構(gòu)和光學(xué)性質(zhì)的第一性原理研究

    邢美靜,梁 凱,姚倩倩,李 實(shí),申 潔,趙 輝

    (天津師范大學(xué)物理與材料科學(xué)學(xué)院,天津300387)

    為探究不同單軸應(yīng)變條件對(duì)堿土金屬氧化物CaO和SrO晶體材料電子結(jié)構(gòu)和光學(xué)性質(zhì)的影響,采用第一性原理計(jì)算方法,利用密度泛函近似理論,對(duì)晶體的能帶結(jié)構(gòu)、態(tài)密度和介電函數(shù)進(jìn)行計(jì)算,并利用克喇末-克朗尼格關(guān)系得出晶體的光學(xué)常數(shù).研究結(jié)果表明:無(wú)應(yīng)變時(shí),CaO和SrO晶體是寬禁帶且具有直接帶隙的絕緣體,施加單軸應(yīng)變使晶體帶隙的寬度變窄,導(dǎo)電性能增強(qiáng),且單軸壓應(yīng)變對(duì)晶體電子結(jié)構(gòu)的影響比拉應(yīng)變對(duì)晶體電子結(jié)構(gòu)影響明顯.2種晶體均為各向同性材料,單軸拉應(yīng)變使靜介電常數(shù)變大,單軸壓應(yīng)變使靜介電常數(shù)變小.光子能量在0~15eV時(shí),2種材料折射率、吸收系數(shù)和反射率等光學(xué)譜的峰位和峰值均受單軸應(yīng)變的影響較大,單軸拉應(yīng)變使得各光學(xué)譜的峰位向能量低的方向移動(dòng),峰值升高;單軸壓應(yīng)變使得光學(xué)譜的峰位向能量高的方向移動(dòng),峰值降低.同時(shí)單軸拉應(yīng)變使晶體靜態(tài)介電常數(shù)和靜折射率增大,壓應(yīng)變使之降低.研究表明單軸應(yīng)變可以有效調(diào)節(jié)CaO和SrO的電子結(jié)構(gòu)和光學(xué)性質(zhì).

    CaO晶體;SrO晶體;第一性原理;單軸拉應(yīng)變;單軸壓應(yīng)變;電子結(jié)構(gòu);光學(xué)性質(zhì)

    金屬氧化物是一類在催化作用、微電子和光電設(shè)備等領(lǐng)域應(yīng)用廣泛的材料[1-4],其中堿土金屬氧化物是地球地下幔的重要成分,是一種廣泛應(yīng)用的工業(yè)科學(xué)技術(shù)基礎(chǔ)材料[5-8].CaO是一種寬帶隙的氧化物,可用于漂白和自動(dòng)加熱包等方面;SrO可應(yīng)用于燃料電池,根據(jù)它的光學(xué)性質(zhì),SrO還可以用來(lái)提高特殊玻璃的質(zhì)量[9].2種材料均為寬帶隙金屬氧化物,在實(shí)驗(yàn)和理論兩方面被廣泛研究.CaO在靜水壓為61 GPa的條件下可以實(shí)現(xiàn)從B1(NaCl)相到B2(CsCl)相的轉(zhuǎn)變[10],且至少在135 GPa條件下才能使B2相處于穩(wěn)定狀態(tài)[11],Kalpana[12]和Baltachea[8]等已對(duì)CaO的相變問(wèn)題進(jìn)行了具體研究;另一方面Mehl等[13]預(yù)測(cè)了CaO的彈性性質(zhì)受壓強(qiáng)的影響.與此同時(shí),SrO的彈性性質(zhì)、熱力學(xué)性質(zhì)和相變等問(wèn)題也有大量研究[2-3,8,14].

    綜上所述可知,現(xiàn)有的大部分研究多集中在晶體材料的彈性性質(zhì)和相變等問(wèn)題上,部分研究涉及靜水壓條件下CaO和SrO的光學(xué)性,對(duì)單軸應(yīng)變條件下二者的電子結(jié)構(gòu)和光學(xué)性質(zhì)研究較少,本研究基于密度泛函理論的第一性原理,討論單軸應(yīng)變條件對(duì)CaO和SrO的電子結(jié)構(gòu)和光學(xué)性質(zhì)的影響.

    1 計(jì)算方法與模型

    本研究從第一性原理出發(fā),采用基于密度泛函理論的平面波贗勢(shì)方法的CASTEP軟件包[15-16]進(jìn)行計(jì)算,利用超軟贗勢(shì)[15,17]描述離子勢(shì)與價(jià)電子間的相互作用勢(shì),其中Ca和Sr各原子的價(jià)電子組態(tài)分別為3s23p64s2和4s24p65s2,O原子的電子組態(tài)為2 s22 p4.電子交換相關(guān)項(xiàng)采用廣義梯度近似(GGA)[18].在倒易空間中,平面波截止能為380.0eV,布里淵區(qū)網(wǎng)格[19]選取4×4×4的格點(diǎn).在對(duì)模型結(jié)構(gòu)優(yōu)化時(shí)采用BFGS[20]算法,自洽循環(huán)計(jì)算收斂精度為每個(gè)原子平均能量小于5×10-6eV,原子間的相互作用力低于0.1eV/nm,最大內(nèi)應(yīng)力為0.02GPa,原子的最大位移為5×10-5nm.

    CaO和SrO屬于立方晶系,空間群為Fm-3m,與NaCl晶體結(jié)構(gòu)相似.CaO晶體原子坐標(biāo)分別為Ca(0,0,0)和O(0.5,0.5,0.5),優(yōu)化后晶格常數(shù)a=b=c=0.4828nm,與其他文獻(xiàn)計(jì)算值0.4819nm符合良好[21],且與實(shí)驗(yàn)測(cè)量值a=b=c=0.481 1 nm非常相近;SrO晶體中的原子坐標(biāo)分別為Sr(0,0,0)、O(0.5,0.5,0.5),優(yōu)化后晶格常數(shù)a= b=c=0.5190nm,與其他研究計(jì)算值0.5181nm和實(shí)驗(yàn)值0.516 0 nm[2]基本一致.單軸應(yīng)變施加在z軸方向,施加單軸應(yīng)變的晶格常數(shù)的相對(duì)變化定義為η=(c-c0)/c0×100%,其中c為施加應(yīng)變后z軸方向晶體的晶格常數(shù),c0為不受應(yīng)變作用時(shí)z軸方向晶體的晶格常數(shù).本研究通過(guò)對(duì)CaO和SrO施加15%單軸拉應(yīng)變和-15%的單軸壓應(yīng)變來(lái)討論單軸應(yīng)變對(duì)晶體電子結(jié)構(gòu)和光學(xué)性質(zhì)的影響.

    2 計(jì)算結(jié)果與討論

    2.1 單軸應(yīng)變對(duì)晶體電子結(jié)構(gòu)的影響

    優(yōu)化后未施加單軸應(yīng)變時(shí),CaO和SrO晶體沿布里淵區(qū)對(duì)稱點(diǎn)方向的能帶結(jié)構(gòu)、總態(tài)密度和分波態(tài)密度如圖1所示.

    圖1 無(wú)應(yīng)變時(shí),CaO和SrO的能帶結(jié)構(gòu)圖和態(tài)密度圖Fig.1 Electronic band structure and corresponding DOS of CaO and SrO without strain

    由圖1(a)可知,在未發(fā)生應(yīng)變的情況下,CaO晶體導(dǎo)帶底和價(jià)帶頂均在布里淵區(qū)的G點(diǎn),帶隙寬度為3.658 eV,說(shuō)明此種材料為直接帶隙的絕緣體.該結(jié)果與文獻(xiàn)[22]中計(jì)算的結(jié)果相近,但比實(shí)驗(yàn)值(7.1 eV)明顯偏小,這是因?yàn)闊o(wú)論采用廣義梯度近似(GGA)還是局域密度近似(LDA)都不能完全描述多電子相互作用,因此均存在禁帶寬度計(jì)算值偏低的問(wèn)題,對(duì)比圖1(b)可知,SrO與CaO的能帶結(jié)構(gòu)存在導(dǎo)帶底和價(jià)帶頂均出現(xiàn)在G點(diǎn)的共同點(diǎn),說(shuō)明2種材料都是直接帶隙的絕緣體(SrO禁帶寬度為3.337 eV).此外,由圖1(a)中CaO總態(tài)密度和各原子分態(tài)密度圖可以看出,CaO的價(jià)帶分為4個(gè)部分:-38.53~-36.58 eV處的態(tài)密度由Ca的s態(tài)貢獻(xiàn);-20.3~-17.45 eV處的峰值主要由Ca 的p態(tài)貢獻(xiàn),同時(shí)O的s態(tài)也有少量貢獻(xiàn);-16.19 eV~-13.37 eV處的峰值主要由O的p態(tài)貢獻(xiàn);價(jià)帶頂部主要由O的p態(tài)貢獻(xiàn),并摻雜Ca的d態(tài).導(dǎo)帶底部的態(tài)密度由Ca的d態(tài)貢獻(xiàn),并有少量的O的p態(tài)貢獻(xiàn).

    施加單軸拉應(yīng)變和壓應(yīng)變后,CaO和SrO晶體沿布里淵區(qū)對(duì)稱點(diǎn)方向的能帶結(jié)構(gòu)和態(tài)密度分別如圖2和圖3所示.

    圖2 施加應(yīng)變后,CaO的能帶結(jié)構(gòu)圖和態(tài)密度圖Fig.2 Electronic band structure and corresponding DOS of CaO under uniaxial strains

    圖2(a)和圖2(b)分別給出單軸拉應(yīng)變和單軸壓應(yīng)變條件下,CaO的能帶結(jié)構(gòu)和態(tài)密度.由能帶結(jié)構(gòu)圖可以看出,無(wú)論是施加拉應(yīng)變(η=15%)還是施加壓應(yīng)變(η=-15%),導(dǎo)帶底和價(jià)帶頂均轉(zhuǎn)移到X點(diǎn),仍然是直接帶隙.2種條件下,CaO的帶隙均變窄,禁帶寬度分別減小為3.254 eV和2.474 eV,η=-15%的單軸壓應(yīng)變對(duì)帶隙寬度的影響明顯比η=15%單軸應(yīng)變影響大.2種應(yīng)變條件對(duì)CaO能帶結(jié)構(gòu)的影響不同,單軸拉應(yīng)變使能帶曲線更平滑,單軸壓應(yīng)變使曲線波動(dòng)起伏更大,能帶寬度變寬.由態(tài)密度圖可以看出,由Ca的d態(tài)電子占據(jù)的導(dǎo)帶底均向低能量方向移動(dòng),從而使帶隙變窄,這與能帶結(jié)構(gòu)的變化一致,因?yàn)槟軒г狡?,態(tài)密度峰值越尖銳,由于拉應(yīng)變使能帶曲線變得更加平緩,所以η=15%單軸拉應(yīng)變情況下,CaO態(tài)密度的峰值升高;而η=-15%單軸壓應(yīng)變條件下CaO態(tài)密度的峰值降低;且在η=-15%單軸壓應(yīng)變條件下CaO態(tài)密度所在區(qū)域變大,造成體系離域性增強(qiáng),這與單軸壓應(yīng)變條件下能帶變寬一致.

    對(duì)比圖3(a)和圖3(b)中SrO在15%單軸拉應(yīng)變和-15%單軸壓應(yīng)變條件下的電子結(jié)構(gòu)可以發(fā)現(xiàn),SrO 和CaO電子能帶結(jié)構(gòu)和態(tài)密度的變化趨勢(shì)基本相同,如單軸應(yīng)變條件下SrO還是直接帶隙,其導(dǎo)帶底和價(jià)帶頂均出現(xiàn)在X點(diǎn).區(qū)別是單軸應(yīng)變對(duì)SrO晶體帶隙影響比較明顯,單軸拉應(yīng)變和壓應(yīng)變的帶隙分別為件下晶體的導(dǎo)電性有所增強(qiáng).2.645 eV和1.697 eV.由態(tài)密度圖(圖3(b))可知,SrO價(jià)帶頂?shù)碾娮釉趬簯?yīng)變條件下向能量高的方向移動(dòng),跨越費(fèi)米面,這是O的p態(tài)電子向高能量方向移動(dòng)造成的.此外,SrO態(tài)密度圖的變化說(shuō)明在單軸壓應(yīng)變條

    圖3 施加應(yīng)變后,SrO的能帶結(jié)構(gòu)圖和態(tài)密度圖Fig.3 Electronic band structure and corresponding DOS of SrO under different uniaxial strains

    2.2 單軸應(yīng)變對(duì)晶體光學(xué)性質(zhì)的影響

    在線性光學(xué)響應(yīng)范圍內(nèi),通常用復(fù)介電函數(shù)

    描述固體宏觀光學(xué)響應(yīng)函數(shù),式(1)中:ε1(ω)為介電實(shí)部,ε2(ω)為介電虛部[22].復(fù)介電函數(shù)ε(ω)與電子和光子的相互作用密切相關(guān),是連接固體材料微觀電子結(jié)構(gòu)和宏觀光學(xué)性質(zhì)的橋梁,通過(guò)ε(ω)可以方便地得到其他各種光譜信息.介電函數(shù)虛部ε2(ω)主要表征電子占有態(tài)與非占有態(tài)間的躍遷,通過(guò)占據(jù)態(tài)和非占據(jù)態(tài)波函數(shù)的動(dòng)量矩陣元可得介電函數(shù)的虛部ε2(ω),再根據(jù)克喇末-克朗尼格(Kramer-Kronig)色散關(guān)系,由ε2(ω)可以得到介電函數(shù)的實(shí)部[23-24].ε1(ω)和ε2(ω)的表示方法分別為

    式(2)和式(3)中:e為電子電量;m為自由電子質(zhì)量;ω為入射光子的頻率;BZ為第一布里淵區(qū);c和v分別為導(dǎo)帶和價(jià)帶;k為倒格失;Ec(k)和Ev(k)為導(dǎo)帶和價(jià)帶的本征能量;Mcv(k)為動(dòng)量躍遷矩陣元.其他光學(xué)常數(shù),如折射率n(ω)、反射率R(ω)、吸收系數(shù)α(ω)和能量損失函數(shù)L(ω)[24]均可由介電函數(shù)求出[25],具體關(guān)系式分別為

    2.2.1 介電函數(shù)

    圖4和圖5分別給出了不同應(yīng)變條件下CaO和SrO復(fù)介電函數(shù)實(shí)部ε1(ω)和虛部ε2(ω)的變化曲線.

    圖4 不同應(yīng)變條件下,CaO晶體復(fù)介電函數(shù)的變化情況Fig.4 Complex dielectric function of CaO under different uniaxial strains

    由圖4(a)可知,無(wú)應(yīng)變時(shí)CaO的靜介電常數(shù)ε1(0)為2.85,施加單軸拉應(yīng)變時(shí),靜介電常數(shù)增加,在η=-15%單軸應(yīng)變情況下靜介電常數(shù)降低,分別變?yōu)?.227和2.696.在8.33~11.57 eV和26.65~28.58 eV這2個(gè)能量范圍內(nèi),CaO復(fù)介電函數(shù)的實(shí)部出現(xiàn)負(fù)值,在這一區(qū)域內(nèi)光無(wú)法傳播,被稱為金屬反射區(qū),說(shuō)明CaO在這些能量范圍內(nèi)表現(xiàn)出金屬性質(zhì).由圖4(b)可以看出,不受應(yīng)變作用時(shí),ε2(ω)有2個(gè)主峰,峰的位置分別在8.21 eV和26.31 eV處.第1個(gè)峰值由價(jià)帶O的p態(tài)與導(dǎo)帶Ca的d態(tài)的電子間躍遷產(chǎn)生,第2個(gè)峰值主要源于價(jià)帶O的s態(tài)和導(dǎo)帶Ca的d態(tài)的電子間躍遷.施加拉應(yīng)變后2個(gè)峰位均稍微向低能量方向移動(dòng)(產(chǎn)生紅移),峰值增大;但在施加壓應(yīng)變條件下,2個(gè)峰均分裂成多個(gè)峰,且峰值變小.這是由于在單軸壓應(yīng)變條件下,體系的離域性特征增強(qiáng)造成的.從計(jì)算結(jié)果可以看出,單軸壓應(yīng)變對(duì)CaO復(fù)介電函數(shù)虛部峰值大小的影響比拉應(yīng)變的影響明顯,且第1個(gè)峰值附近的曲線受單軸應(yīng)變影響比較明顯.

    對(duì)比圖5可知,在無(wú)應(yīng)變情況下SrO的主峰個(gè)數(shù)比CaO的主峰個(gè)數(shù)多,但SrO晶體復(fù)介電函數(shù)隨應(yīng)變條件下的變化趨勢(shì)與CaO的基本相同.

    圖5 不同應(yīng)變條件下,SrO晶體復(fù)介電函數(shù)的變化情況Fig.5 Complex dielectric function of SrO under different uniaxial strains

    2.2.2 光學(xué)常數(shù)

    圖6和圖7為不同應(yīng)變條件下CaO和SrO晶體折射率、吸收系數(shù)、反射率和能量損失函數(shù)隨光子能量的變化曲線.

    由圖6(a)和圖7(a)可知,在無(wú)應(yīng)變、η=15%單軸拉應(yīng)變和η=-15%單軸壓應(yīng)變條件下,CaO的靜態(tài)折射率分別為1.689、1.796和1.641,SrO的靜態(tài)折射率分別為1.697、1.802和1.623,該結(jié)果表明在單軸拉應(yīng)變情況下靜態(tài)折射率增大,單軸壓應(yīng)變條件下靜態(tài)折射率降低,折射率隨光子能量變化趨勢(shì)與復(fù)介電函數(shù)實(shí)部的變化規(guī)律基本一致,符合麥克斯韋關(guān)系ε1(ω)= n2.在8.33~11.57 eV和26.65~28.58 eV這2個(gè)區(qū)域內(nèi),折射率趨于0,與復(fù)介電函數(shù)實(shí)部金屬反射區(qū)相對(duì)應(yīng).

    由圖6(b)和圖7(b)可以看出,不同應(yīng)變條件下CaO和SrO吸收系數(shù)α(ω)的吸收范圍均分為2個(gè)區(qū)域.CaO在未發(fā)生應(yīng)變時(shí)吸收范圍為2.52~12.46 eV和17.41~29.69 eV,在二者中間的區(qū)域無(wú)吸收,表明在此能量區(qū)域CaO是透明的.無(wú)應(yīng)變時(shí),在CaO第1個(gè)吸收范圍內(nèi),能量為8.94 eV處存在一個(gè)峰值0.24.當(dāng)施加單軸應(yīng)變后,峰值變化幅度變小,且在拉應(yīng)變情況下峰值向能量低的方向移動(dòng),產(chǎn)生紅移,峰值出現(xiàn)在能量為7.86 eV處,這是由于導(dǎo)帶中Ca的d態(tài)向低能量方向移動(dòng),使電子躍遷需要吸收的能量降低.而在施加單軸壓應(yīng)變的條件下,7.73 eV處出現(xiàn)1個(gè)弱峰,隨著能量的增加在10.16 eV位置上出現(xiàn)1個(gè)主峰,峰值為0.19,吸收譜產(chǎn)生藍(lán)移,原因是壓應(yīng)變條件下態(tài)密度離域性增強(qiáng),導(dǎo)帶能級(jí)展寬,且峰值向能量高的方向移動(dòng).在第2個(gè)吸收范圍內(nèi),應(yīng)變條件對(duì)峰位影響很小,單軸拉應(yīng)變使吸收峰升高,壓應(yīng)變使峰值降低,變化比拉應(yīng)變情況明顯.SrO吸收系數(shù)隨能量的變化趨勢(shì)與CaO基本上相同.

    圖6 不同應(yīng)變條件下,CaO光學(xué)常數(shù)的變化情況Fig.6 Optical constants of CaO under different uniaxial strains

    圖7 不同應(yīng)變條件下,SrO光學(xué)常數(shù)的變化情況Fig.7 Optical constants of SrO under different uniaxial strains

    由圖6(c)可知,無(wú)應(yīng)變時(shí),CaO的反射率R(0)= 6.6%,光子能量在11.47 eV位置處反射率達(dá)到峰值96%,趨近于100%,因?yàn)樵擖c(diǎn)位于金屬反射區(qū),表現(xiàn)為金屬性質(zhì),在該點(diǎn)光無(wú)法傳播,因此反射率接近于100%;同時(shí)在28.48 eV位置上另一個(gè)反射率峰的峰值為87%.施加單軸拉應(yīng)變后,第1個(gè)反射峰值向能量低的方向移動(dòng),且比無(wú)應(yīng)變時(shí)峰值增大,同時(shí)第2個(gè)峰位置變化很小,峰值稍微增大.而在壓應(yīng)變條件下,第1個(gè)峰向能量高的方向移動(dòng),且峰值增大,第2個(gè)峰值顯著降低.SrO反射率隨能量的變化趨勢(shì)受單軸應(yīng)變的影響與CaO的變化趨勢(shì)基本相同.

    L(ω)是重要的光學(xué)參數(shù),可以表征電子通過(guò)均勻電子介質(zhì)時(shí)能量的損失情形,而能量損失峰則可以描繪等離子的共振頻率.由圖6(d)和圖7(d)可以看出,單軸應(yīng)變對(duì)2種材料的影響存在很大差別.CaO在無(wú)應(yīng)變時(shí)的最大能量損失峰達(dá)到399.98,施加單軸應(yīng)變后,損失峰明顯減小,說(shuō)明等離子共振頻率減?。粏屋S拉應(yīng)變和單軸壓應(yīng)變使損失峰位分別向能量低和能量高的方向移動(dòng).而SrO在單軸拉應(yīng)變下?lián)p失峰明顯升高,峰值向能量低的方向移動(dòng);在單軸壓應(yīng)變條件下,損失峰降低且向能量高的方向移動(dòng).

    3 結(jié)論

    采用基于密度泛函理論的第一性原理計(jì)算方法,對(duì)單軸應(yīng)變條件下CaO和SrO的電子結(jié)構(gòu)及光學(xué)性質(zhì)進(jìn)行系統(tǒng)研究.計(jì)算結(jié)果表明:

    (1)單軸應(yīng)變可以調(diào)節(jié)CaO和SrO的帶隙,使帶隙變窄,從而改變材料的導(dǎo)電性,其主要原因是Ca和Sr導(dǎo)帶底中的d態(tài)電子向低能量方向移動(dòng),且SrO價(jià)帶頂O的p態(tài)電子向高能量方向移動(dòng).

    (2)單軸壓應(yīng)變使電子體系的離域性增強(qiáng),且單軸壓應(yīng)變對(duì)電子結(jié)構(gòu)的影響比拉應(yīng)變的影響大.

    (3)光子能量在0~15 eV能量范圍內(nèi),2種材料的復(fù)介電函數(shù)、折射率、吸收系數(shù)和反射率的峰位和峰值受單軸應(yīng)變的影響明顯,單軸拉應(yīng)變使光學(xué)譜的峰位向能量低的方向移動(dòng),峰值升高;單軸壓應(yīng)變使光學(xué)譜的峰位向能量高的方向移動(dòng),峰值降低.

    (4)CaO和SrO光學(xué)性質(zhì)受單軸應(yīng)變影響的變化趨勢(shì)基本相同,但能量損失情況差別較大,單軸應(yīng)變使CaO的能量損失明顯減小,但單軸拉應(yīng)變使SrO的能量損失增大.

    [1]DAI W,SHUI Z H,LI K,et al.First-principle investigations of CaO (100)surface and adsorption of H2O on CaO(100)[J].Computational and Theoretical Chemistry,2011,967(1):185-190.

    [2]GHEBOULI B,GHEBOULI M A,F(xiàn)ATMI M,et al.First-principles calculations of structural,elastic,electronic and optical properties of XO(X=Ca,Sr and Ba)compounds under pressure effect[J].Materials Science in Semiconductor Processing,2010,13(2):92-101.

    [3]DUAN Y,QIN L,TANG G,et al.First-principles study of ground-and metastable-state properties of XO(X=Be,Mg,Ca,Sr,Ba,Zn and Cd)[J].Eur Phys J B,2008,66(2):201-209.

    [4]ANDREW J L,DAVID M F,DAVID O S,et al.Structural,energetic and electronic properties of(100)surfaces for alkaline earth metal oxides as calculated with hybrid density functional theory[J].Surface Science,2015,642:58-65.

    [5]SU M X,YANG R,LI M.Biodiesel production from hempseed oil using alkaline earth metal oxides supporting copper oxide as bifunctional catalysts for transesterification and selective hydrogenation[J].Fuel,2013,103(1):398-407.

    [6]WANG Z Y,LI X H,ZHANG Y,et al.Effect of alkaline-earth metal oxides on Cu/SiO2-Al2O3catalyst for vapor-phase synthesis of 3-methylindole from glycerol and aniline[J].Chinese Journal of Catalysis,2012,33(7):1139-1145.

    [7]SKORODUMOVA N V,HERMANSSON K,JOHANSSON B.Structural and electronic properties of the(100)surface and bulk of alkaline-earth metal oxides[J].Physical Review B,2005,72(12):125414.

    [8]BALTACHEA H,KHENATA R,SAHNOUN M,et al.Full potential calculation of structural,electronic and elastic properties of alkaline earth oxides MgO,CaO and SrO[J].Physica B,2004,344(1/2/3/4):334-342.

    [9]FRANCISCO G C,JUAN B M.Combustion synthesis process for the rapid preparation of high-purity SrO powders[J].Materials Science-Poland,2014,32(4):682-687.

    [10]JEANLOZ R,AHRENS T,MAO H K,et al.B1-B2 transition in calcium oxide from shock-wave and diamond-cell experiments[J].Science,1979,206(4420):829.

    [11]DENG Y,JIA O H,CHEN X R,et al.Phase transition and elastic constants of CaO from first-principle calculations[J].Physica B,2007,392(1/2):229-232.

    [12]KALPANA G,PALANIVEL B,RAJAGOPALAN M.Electronic arid structural properties of alkaline-earth oxides under high pressure[J].Phys Rev B,1995,52:4-7.

    [13]MEHL M J,HEMLEY R J,BOYER L L.Potential-induced breathing model for the elastic moduli and high-pressure behavior of the cubic alkaline-earth oxides[J].Phys Rev B,1986,33(12):8685.

    [14]TAURIAN O E,SPRINGBORG M,CHRISTENSEN N E.Self-consistent electronic structures of MgO and SrO[J].Solid State Commun,1985,55(4):351-355.

    [15]VANDERBILT D.Soft self-consistent pseudopotentials in a generalized eigenvalue formalism[J].Phys Rev B,1990,41(11):7892-7895.

    [16]SEGALL M D,LINDAN P J D,PROBERT M J,et al.First-principles simulation:ideas,illustrations and the CASTE Pcode[J].Condens Matter,2002,14(11):2717-2744.

    [17]MILMAN V,WINKLER B,WHITE J A,t al.Electronic structure,properties,and phase stability of inorganic crystals:a pseudopotential plane-wave study[J].Int J Quantun Chem,2000,77(5):895-910.

    [18]PERDEW J P,BURKE K,ERNZERHOF M.Generalized gradient approximations made simple[J].Phys Rev Lett,1996,77(18):3865-3868.

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

    [20]FISCHER T H,ALMLOF J.General methods for geometry and wave function optimization[J].The Journal of Physical Chemistry,1992,96 (24):9768-9774.

    [21]MEDEIROSA S K,ALBUQUERQUEA E L,MAIA J F F,et al.Structural and optical properties of CaO[J].Microelectronics Journal,2005,36(11):1058-1061.

    [22]AMBROSCH C,SOFO J O.Linear optical properties of solids within the full potential linearized augmented plane wave method[J].Comput Phys Commun,2006,1751:1-14.

    [23]FOX M.Optical Properties of Solids[M].New York:Oxford University Press,2001.

    [24]RAVINDRAN P,DELIN A,AHUJA R,et al.Optical properties of monoclinic SnI2from relativistic first-principles theory[J].Phys Rev B,1997,56:6851-6861.

    (責(zé)任編校 亢原彬)

    First-principle study on the electronic structures and optical properties of CaO and SrO under uniaxial strains

    XING Meijing,LIANG Kai,YAO Qianqian,LI Shi,SHEN Jie,ZHAO Hui
    (College of Physics and Materials Science,Tianjin Normal University,Tianjin 300387,China)

    In order to investigate the impact of uniaxial strain on the electronic structure and optical properties of alkaline earth metal oxides CaO and SrO,the energy band structure,density of states and the dielectric function of the crystals were calculated using the first-principle calculation and density-functional theory(DFT)method,and the optical constants were obtained based on Kramers-Kroning relations.The calculated results show that CaO and SrO crystals are both wide and direct band gap insulators in the case of no strain.The band gaps widths of two materials decrease with both uniaxial tensile and compressive strain,in these two cases their conductivities are enhanced,and the influence of compressive strain on the electronic structure of CaO(SrO)is more obvious than tensile strain.The two crystals are both isotropic materials,the uniaxial tensile strain makes the static dielectric constant larger,and the uniaxial compressive strain makes the static dielectric constant smaller.The peaks of positions and values of refractive index of the two material refractive index,absorption coefficient and reflectivity are influenced obviously by uniaxial strain in the energy range from 0 to 15 eV.The uniaxial tensile strain makes each optical spectrum peak moving towards the direction of lower energy and the peak value becoming higher,and the uniaxial compressive strain makes optical spectrum peak position moving towards the direction of higher energy and the values of peaks becoming lower.The static dielectric constant and refractive index of two crystals are both increase with uniaxial tensile strain and decrease with uniaxial compressive strain.The study shows that uniaxial strain can effectively adjust the electronic structure and optical properties of CaO and SrO.

    CaO crystal;SrO crystal;first-principle;uniaxial tensile strain;uniaxial compressive strain;electronic structure;optical properties

    O482.31

    A

    1671-1114(2016)03-0021-07

    2015-12-02

    邢美靜(1991—),女,碩士研究生.

    趙 輝(1965—),男,教授,主要從事凝聚態(tài)理論和計(jì)算物理方面的研究.

    猜你喜歡
    導(dǎo)帶能帶單軸
    關(guān)于某型系列產(chǎn)品“上下導(dǎo)帶二次收緊”合并的工藝技術(shù)
    吃東西時(shí)注意多
    炮射彈藥滑動(dòng)導(dǎo)帶環(huán)結(jié)構(gòu)設(shè)計(jì)
    單軸壓縮條件下巖石峰后第Ⅱ種類型應(yīng)力——應(yīng)變曲線的新解釋
    汽車轉(zhuǎn)向管柱吸能帶變形研究和仿真優(yōu)化
    Life OR Death Decision
    CFRP-鋼復(fù)合板的單軸拉伸力學(xué)性能
    導(dǎo)帶式數(shù)碼噴射印花機(jī)運(yùn)行過(guò)程中常見(jiàn)疵病
    打捆機(jī)導(dǎo)帶槽裝置及捆帶盤脹縮裝置的改進(jìn)
    單軸應(yīng)變Si NMOS電流模型研究
    国产精品久久久久久久久免 | 桃红色精品国产亚洲av| 中文字幕熟女人妻在线| 婷婷六月久久综合丁香| 不卡一级毛片| 99久久九九国产精品国产免费| 亚洲人成网站在线播| 国产伦人伦偷精品视频| 日韩精品青青久久久久久| 久久精品91蜜桃| 天天一区二区日本电影三级| 夜夜躁狠狠躁天天躁| 99久久精品国产亚洲精品| 亚洲精品影视一区二区三区av| 免费看美女性在线毛片视频| 九色成人免费人妻av| 亚洲精品色激情综合| 国内少妇人妻偷人精品xxx网站| 一个人看的www免费观看视频| 中文字幕高清在线视频| 国产精品一区二区免费欧美| 国产精品伦人一区二区| 久久精品国产清高在天天线| 日本熟妇午夜| 99久国产av精品| 欧美一区二区亚洲| 不卡一级毛片| 国产成人啪精品午夜网站| 在线观看美女被高潮喷水网站 | 久久国产乱子伦精品免费另类| 国产亚洲精品久久久久久毛片| 天天一区二区日本电影三级| 国产亚洲精品综合一区在线观看| 怎么达到女性高潮| 久久久久精品国产欧美久久久| 精品午夜福利视频在线观看一区| 日韩有码中文字幕| 欧美黑人巨大hd| 国内毛片毛片毛片毛片毛片| 久久精品国产亚洲av香蕉五月| 国产麻豆成人av免费视频| 亚洲自偷自拍三级| 精品人妻偷拍中文字幕| 内射极品少妇av片p| www.熟女人妻精品国产| 亚洲最大成人手机在线| 亚洲精品粉嫩美女一区| 成人毛片a级毛片在线播放| 欧美日本视频| 欧美国产日韩亚洲一区| 麻豆成人午夜福利视频| 最近最新免费中文字幕在线| 亚洲自偷自拍三级| 亚洲精品日韩av片在线观看| 亚洲经典国产精华液单 | 每晚都被弄得嗷嗷叫到高潮| 人人妻人人看人人澡| 日本三级黄在线观看| 国产单亲对白刺激| 极品教师在线视频| 国产亚洲欧美98| 成年人黄色毛片网站| 一区二区三区免费毛片| 91久久精品国产一区二区成人| 夜夜爽天天搞| 亚洲av日韩精品久久久久久密| 少妇丰满av| 熟女电影av网| 网址你懂的国产日韩在线| 黄色一级大片看看| 亚洲va日本ⅴa欧美va伊人久久| 人人妻人人看人人澡| 日韩欧美在线二视频| 久久精品国产99精品国产亚洲性色| 亚洲五月婷婷丁香| 国产成人影院久久av| 性插视频无遮挡在线免费观看| 亚洲狠狠婷婷综合久久图片| 国产亚洲欧美在线一区二区| 亚洲,欧美精品.| 国产免费av片在线观看野外av| 人人妻人人澡欧美一区二区| 久久精品国产亚洲av天美| 久久6这里有精品| www.色视频.com| 欧美日韩黄片免| 精品人妻一区二区三区麻豆 | 黄色女人牲交| 此物有八面人人有两片| 日本黄色视频三级网站网址| 中文亚洲av片在线观看爽| 日韩中文字幕欧美一区二区| 999久久久精品免费观看国产| 91久久精品电影网| 99久久久亚洲精品蜜臀av| 国产乱人视频| 人人妻人人澡欧美一区二区| 日韩成人在线观看一区二区三区| 久久国产乱子免费精品| 日韩欧美国产一区二区入口| 免费观看精品视频网站| 成人性生交大片免费视频hd| 亚洲最大成人av| 91麻豆精品激情在线观看国产| 免费看a级黄色片| 成年女人看的毛片在线观看| 又黄又爽又免费观看的视频| 国产伦精品一区二区三区四那| 在线播放国产精品三级| 首页视频小说图片口味搜索| 亚洲国产精品999在线| 欧美黄色淫秽网站| 中文字幕精品亚洲无线码一区| 男插女下体视频免费在线播放| 91在线精品国自产拍蜜月| 国产国拍精品亚洲av在线观看| 欧美精品国产亚洲| 成年女人永久免费观看视频| 琪琪午夜伦伦电影理论片6080| 久久久久久大精品| 亚洲七黄色美女视频| 亚洲欧美日韩卡通动漫| 我要看日韩黄色一级片| 国产一区二区亚洲精品在线观看| 亚洲av成人av| 免费高清视频大片| 色综合站精品国产| 国产黄色小视频在线观看| 性色avwww在线观看| 亚洲av第一区精品v没综合| 欧美成人性av电影在线观看| 午夜老司机福利剧场| 国产精品美女特级片免费视频播放器| xxxwww97欧美| 欧美高清性xxxxhd video| 十八禁国产超污无遮挡网站| 午夜免费男女啪啪视频观看 | 亚洲18禁久久av| 久久久久久久午夜电影| 老鸭窝网址在线观看| 精品一区二区免费观看| 亚洲欧美日韩无卡精品| 国产精品久久久久久亚洲av鲁大| 欧美乱色亚洲激情| 久久热精品热| 亚洲激情在线av| 亚洲av免费在线观看| 人人妻,人人澡人人爽秒播| 女同久久另类99精品国产91| 999久久久精品免费观看国产| 免费av观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 免费av观看视频| 免费无遮挡裸体视频| 12—13女人毛片做爰片一| 国产男靠女视频免费网站| av欧美777| 宅男免费午夜| 少妇的逼好多水| 俄罗斯特黄特色一大片| 色精品久久人妻99蜜桃| 国产精品久久视频播放| 日韩高清综合在线| а√天堂www在线а√下载| 亚洲av成人av| 黄片小视频在线播放| 久久亚洲精品不卡| 人人妻人人澡欧美一区二区| 成人毛片a级毛片在线播放| 不卡一级毛片| 国产午夜福利久久久久久| 午夜两性在线视频| 九九在线视频观看精品| 日本成人三级电影网站| 国产亚洲精品久久久久久毛片| 国语自产精品视频在线第100页| 国产伦在线观看视频一区| av福利片在线观看| 亚洲中文字幕日韩| 韩国av一区二区三区四区| 在线国产一区二区在线| 欧美+日韩+精品| 国产白丝娇喘喷水9色精品| 色哟哟哟哟哟哟| 97人妻精品一区二区三区麻豆| 亚洲av免费在线观看| 亚洲aⅴ乱码一区二区在线播放| 午夜精品在线福利| 性欧美人与动物交配| 可以在线观看的亚洲视频| 国产精品美女特级片免费视频播放器| 精品一区二区三区视频在线| 又爽又黄无遮挡网站| 观看免费一级毛片| 中文资源天堂在线| 精品不卡国产一区二区三区| 国产亚洲精品av在线| 国产精品亚洲一级av第二区| av在线观看视频网站免费| 国产蜜桃级精品一区二区三区| 中文字幕熟女人妻在线| 男人的好看免费观看在线视频| 亚洲成a人片在线一区二区| 真实男女啪啪啪动态图| 好看av亚洲va欧美ⅴa在| .国产精品久久| 免费看a级黄色片| 五月玫瑰六月丁香| 91久久精品国产一区二区成人| 久久久久久久午夜电影| 波多野结衣高清作品| www.999成人在线观看| 最新中文字幕久久久久| 久久国产乱子伦精品免费另类| 亚洲国产精品sss在线观看| 国内少妇人妻偷人精品xxx网站| 欧美三级亚洲精品| 国产一级毛片七仙女欲春2| 色哟哟哟哟哟哟| 亚洲美女黄片视频| www.999成人在线观看| 国产伦精品一区二区三区四那| www.www免费av| 午夜激情福利司机影院| 夜夜夜夜夜久久久久| 又黄又爽又刺激的免费视频.| 欧美最黄视频在线播放免费| 禁无遮挡网站| 1024手机看黄色片| 久9热在线精品视频| 99在线视频只有这里精品首页| 51国产日韩欧美| 99久久精品国产亚洲精品| 一级黄片播放器| 色视频www国产| 欧美另类亚洲清纯唯美| 高清在线国产一区| 99热精品在线国产| 搞女人的毛片| 亚洲色图av天堂| 脱女人内裤的视频| 天堂网av新在线| 日韩欧美在线二视频| 久久久久久久亚洲中文字幕 | 成人av一区二区三区在线看| av国产免费在线观看| 日本与韩国留学比较| 国产免费av片在线观看野外av| 久久久久久国产a免费观看| 色精品久久人妻99蜜桃| 99热6这里只有精品| 99热精品在线国产| 变态另类成人亚洲欧美熟女| 亚洲精品一区av在线观看| 88av欧美| 欧美zozozo另类| 亚洲成av人片免费观看| 熟女人妻精品中文字幕| 又爽又黄无遮挡网站| 免费人成在线观看视频色| 亚洲 欧美 日韩 在线 免费| 欧美午夜高清在线| 亚洲av电影在线进入| 久久午夜福利片| 直男gayav资源| 亚洲av五月六月丁香网| 欧美在线黄色| 日日干狠狠操夜夜爽| 亚洲精品456在线播放app | av黄色大香蕉| 国产老妇女一区| 国产蜜桃级精品一区二区三区| 99国产精品一区二区蜜桃av| 色综合欧美亚洲国产小说| 国产熟女xx| 午夜免费成人在线视频| 亚洲美女黄片视频| 亚洲av.av天堂| 国产精品日韩av在线免费观看| 久久天躁狠狠躁夜夜2o2o| 毛片女人毛片| 尤物成人国产欧美一区二区三区| 观看免费一级毛片| 亚洲国产高清在线一区二区三| 欧美日韩黄片免| 精品乱码久久久久久99久播| 久久精品影院6| www.www免费av| 精品人妻熟女av久视频| 欧美国产日韩亚洲一区| 欧美色欧美亚洲另类二区| 可以在线观看毛片的网站| 网址你懂的国产日韩在线| 91字幕亚洲| 久久亚洲精品不卡| 亚洲av日韩精品久久久久久密| 88av欧美| 丁香欧美五月| 十八禁网站免费在线| а√天堂www在线а√下载| 69av精品久久久久久| 在线观看美女被高潮喷水网站 | 日本精品一区二区三区蜜桃| 十八禁网站免费在线| 久久久久久久久大av| 亚洲欧美日韩无卡精品| 久久久成人免费电影| 免费高清视频大片| 偷拍熟女少妇极品色| 亚洲五月天丁香| 国产人妻一区二区三区在| 亚洲欧美精品综合久久99| 级片在线观看| 亚洲av中文字字幕乱码综合| 亚洲欧美日韩高清在线视频| 窝窝影院91人妻| 内射极品少妇av片p| 人妻制服诱惑在线中文字幕| 99久久久亚洲精品蜜臀av| 色综合亚洲欧美另类图片| 偷拍熟女少妇极品色| 好男人在线观看高清免费视频| 国产精品久久电影中文字幕| 18+在线观看网站| 综合色av麻豆| 国产欧美日韩精品亚洲av| 乱人视频在线观看| www.www免费av| 国产成人啪精品午夜网站| 91在线观看av| 成人特级黄色片久久久久久久| 村上凉子中文字幕在线| 久久人人爽人人爽人人片va | 嫩草影院精品99| 69av精品久久久久久| aaaaa片日本免费| 免费在线观看日本一区| 在线天堂最新版资源| 国产主播在线观看一区二区| 免费看美女性在线毛片视频| 国产亚洲欧美98| 麻豆一二三区av精品| 12—13女人毛片做爰片一| 日韩欧美三级三区| 国产精品免费一区二区三区在线| 性色av乱码一区二区三区2| 国产黄色小视频在线观看| 很黄的视频免费| 狂野欧美白嫩少妇大欣赏| 欧美黄色淫秽网站| 亚洲美女搞黄在线观看 | 午夜福利在线观看免费完整高清在 | 99热这里只有是精品在线观看 | 久久久久久九九精品二区国产| 国产高清激情床上av| 嫩草影院新地址| 亚洲一区高清亚洲精品| 国产在线男女| 日本免费一区二区三区高清不卡| 国产成年人精品一区二区| 久久久久国内视频| 麻豆av噜噜一区二区三区| netflix在线观看网站| 亚洲国产精品sss在线观看| 深爱激情五月婷婷| 一级作爱视频免费观看| 日本三级黄在线观看| 免费看光身美女| 国产男靠女视频免费网站| 亚洲激情在线av| 免费观看人在逋| 亚洲人成电影免费在线| 亚洲男人的天堂狠狠| 亚洲在线自拍视频| 国产在线男女| 可以在线观看的亚洲视频| 老司机午夜十八禁免费视频| 久久国产精品人妻蜜桃| 99久久九九国产精品国产免费| 婷婷精品国产亚洲av在线| 亚洲av免费在线观看| 国产精品久久视频播放| 如何舔出高潮| 搞女人的毛片| 亚洲真实伦在线观看| 欧美精品国产亚洲| 欧美成人一区二区免费高清观看| 国产单亲对白刺激| 国产精品不卡视频一区二区 | 国产午夜福利久久久久久| 国产成年人精品一区二区| 青草久久国产| 亚洲在线自拍视频| 日本免费一区二区三区高清不卡| 88av欧美| 久久性视频一级片| 精品久久久久久成人av| 亚洲va日本ⅴa欧美va伊人久久| 日韩欧美 国产精品| 久久久久久久久久黄片| 亚洲欧美日韩高清专用| 国产精品一区二区免费欧美| 亚洲av二区三区四区| 亚洲成人精品中文字幕电影| 99在线人妻在线中文字幕| 久9热在线精品视频| 亚洲在线自拍视频| 国产免费一级a男人的天堂| 日本精品一区二区三区蜜桃| 成人午夜高清在线视频| 人妻久久中文字幕网| 欧美日韩瑟瑟在线播放| 日韩中字成人| 免费人成视频x8x8入口观看| 久久草成人影院| 在线十欧美十亚洲十日本专区| 国产野战对白在线观看| 欧美性猛交╳xxx乱大交人| 在线看三级毛片| 99热这里只有是精品50| 成人高潮视频无遮挡免费网站| 成人av一区二区三区在线看| www.色视频.com| 超碰av人人做人人爽久久| a在线观看视频网站| 桃色一区二区三区在线观看| 波多野结衣巨乳人妻| 国语自产精品视频在线第100页| 2021天堂中文幕一二区在线观| av在线蜜桃| 亚洲av免费在线观看| 禁无遮挡网站| 亚洲在线自拍视频| 国产精品一区二区性色av| 精品午夜福利在线看| 国产欧美日韩一区二区精品| 亚洲精品久久国产高清桃花| 亚洲人成网站在线播| 99久久九九国产精品国产免费| 草草在线视频免费看| 国产成人影院久久av| 国产精品久久久久久久久免 | 亚洲欧美日韩东京热| 亚洲乱码一区二区免费版| av视频在线观看入口| 级片在线观看| 欧美最新免费一区二区三区 | 欧美日本亚洲视频在线播放| 国产精品一区二区三区四区久久| 禁无遮挡网站| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 精品久久久久久成人av| 亚洲最大成人中文| 五月伊人婷婷丁香| 性色av乱码一区二区三区2| 最新在线观看一区二区三区| 最近最新中文字幕大全电影3| .国产精品久久| 国产精品,欧美在线| 熟女人妻精品中文字幕| 日日夜夜操网爽| 男人舔奶头视频| 久久久精品欧美日韩精品| 很黄的视频免费| 亚洲av五月六月丁香网| 国产不卡一卡二| 精品欧美国产一区二区三| 亚洲熟妇熟女久久| 亚洲成人精品中文字幕电影| 成人特级av手机在线观看| 国产美女午夜福利| 18+在线观看网站| 97超视频在线观看视频| 午夜激情欧美在线| 欧美日韩瑟瑟在线播放| 久久精品国产清高在天天线| 久久久国产成人免费| 色5月婷婷丁香| 午夜a级毛片| 欧美成狂野欧美在线观看| 18+在线观看网站| 午夜福利欧美成人| 免费观看的影片在线观看| 免费av观看视频| 男女做爰动态图高潮gif福利片| 青草久久国产| 一本综合久久免费| 丰满人妻一区二区三区视频av| 欧美激情在线99| 桃色一区二区三区在线观看| 搞女人的毛片| 亚洲精品日韩av片在线观看| 成年人黄色毛片网站| 乱码一卡2卡4卡精品| 国产精品一区二区性色av| 亚洲精品色激情综合| 男女视频在线观看网站免费| 简卡轻食公司| 天堂av国产一区二区熟女人妻| 日韩欧美 国产精品| 变态另类丝袜制服| 国产精品久久久久久久久免 | 欧美潮喷喷水| 国产野战对白在线观看| 亚洲精华国产精华精| 日本在线视频免费播放| 我的老师免费观看完整版| 久久久精品欧美日韩精品| 99久国产av精品| 精品一区二区三区视频在线观看免费| 久久久久久大精品| 午夜福利成人在线免费观看| 1000部很黄的大片| 欧美激情在线99| 国产伦一二天堂av在线观看| 国产av不卡久久| 真实男女啪啪啪动态图| 精品久久久久久,| 在线观看午夜福利视频| 美女被艹到高潮喷水动态| 少妇人妻一区二区三区视频| 性色avwww在线观看| av专区在线播放| 每晚都被弄得嗷嗷叫到高潮| 嫩草影院入口| 亚洲在线观看片| 亚洲真实伦在线观看| a在线观看视频网站| 日本熟妇午夜| 日韩欧美三级三区| 精品无人区乱码1区二区| 亚洲色图av天堂| 色综合站精品国产| 日本一本二区三区精品| 男人舔奶头视频| 毛片一级片免费看久久久久 | 黄片小视频在线播放| 一级作爱视频免费观看| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久久久久久黄片| 国产视频内射| 久久香蕉精品热| 亚洲狠狠婷婷综合久久图片| 内地一区二区视频在线| 一区二区三区高清视频在线| 国产伦人伦偷精品视频| 国产成人a区在线观看| 最好的美女福利视频网| www.www免费av| 亚洲,欧美,日韩| 深夜精品福利| 国产精品亚洲av一区麻豆| 色哟哟哟哟哟哟| 毛片一级片免费看久久久久 | 成年人黄色毛片网站| a级一级毛片免费在线观看| 欧美日本亚洲视频在线播放| av视频在线观看入口| 乱码一卡2卡4卡精品| 日本五十路高清| 高清毛片免费观看视频网站| 欧美激情国产日韩精品一区| 欧美黑人欧美精品刺激| 亚洲av不卡在线观看| 国产精品一区二区性色av| av视频在线观看入口| 99久久精品一区二区三区| 丰满的人妻完整版| 亚洲精品456在线播放app | 国产成人欧美在线观看| 身体一侧抽搐| 男女视频在线观看网站免费| 可以在线观看毛片的网站| 黄片小视频在线播放| 亚洲av不卡在线观看| 九九在线视频观看精品| 少妇裸体淫交视频免费看高清| 男女那种视频在线观看| 日韩欧美国产一区二区入口| 两个人视频免费观看高清| 中文字幕久久专区| 成人av一区二区三区在线看| 午夜亚洲福利在线播放| 男人狂女人下面高潮的视频| 亚洲一区二区三区不卡视频| 久久欧美精品欧美久久欧美| 成人av在线播放网站| 欧美不卡视频在线免费观看| 十八禁网站免费在线| 国产精品亚洲av一区麻豆| 悠悠久久av| 69av精品久久久久久| 精品国内亚洲2022精品成人| 少妇被粗大猛烈的视频| 好男人在线观看高清免费视频| 国产精品亚洲美女久久久| 国产成年人精品一区二区| 久久午夜福利片| 好看av亚洲va欧美ⅴa在| 欧美午夜高清在线| 精品无人区乱码1区二区| 亚洲 欧美 日韩 在线 免费| 免费大片18禁| 人妻丰满熟妇av一区二区三区| 国产乱人伦免费视频| 俄罗斯特黄特色一大片| 丁香六月欧美| 亚洲av美国av| 日韩欧美在线二视频| 久久人人精品亚洲av| 欧美成狂野欧美在线观看| 婷婷精品国产亚洲av| а√天堂www在线а√下载| 搡老熟女国产l中国老女人| or卡值多少钱|