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

    鈣鈦礦APbI3結(jié)構(gòu)穩(wěn)定性及光電性質(zhì)的理論研究?

    2017-08-01 17:16:48劉娜危陽馬新國祝林徐國旺楚亮黃楚云
    物理學報 2017年5期
    關(guān)鍵詞:光吸收價帶導帶

    劉娜 危陽 馬新國 祝林 徐國旺 楚亮 黃楚云

    1)(湖北工業(yè)大學理學院,武漢 430068)

    2)(湖北工業(yè)大學,太陽能高效利用湖北省協(xié)同創(chuàng)新中心,武漢 430068)

    3)(南京郵電大學理學院,南京 210046)

    鈣鈦礦APbI3結(jié)構(gòu)穩(wěn)定性及光電性質(zhì)的理論研究?

    劉娜1)危陽1)馬新國1)2)?祝林1)徐國旺1)2)楚亮3)黃楚云1)2)?

    1)(湖北工業(yè)大學理學院,武漢 430068)

    2)(湖北工業(yè)大學,太陽能高效利用湖北省協(xié)同創(chuàng)新中心,武漢 430068)

    3)(南京郵電大學理學院,南京 210046)

    (2016年7月29日收到;2016年12月4日收到修改稿)

    采用基于色散修正的平面波超軟贗勢方法研究了鈣鈦礦材料APbI3結(jié)構(gòu)中四種陽離子Cs+,NH+4,MA+,FA+分別處于A位時,其結(jié)構(gòu)的穩(wěn)定性、電子結(jié)構(gòu)及光學性質(zhì).研究結(jié)果顯示受陽離子種類和尺寸的影響,PbI基體骨架發(fā)生不同程度的扭曲,A位置陽離子(除Cs+外)半徑越大,其與PbI基體骨架之間的作用力越強,在MAPbI3和FAPbI3中PbI6八面體顯示出較大的電偶極矩.計算得出的能帶結(jié)構(gòu)表明,四種體系在費米能級附近的能帶相似,即價帶頂均由I 5p軌道組成,導帶底由Pb 6p軌道和部分I 5p軌道雜化而成.電子結(jié)構(gòu)和光學性質(zhì)的差異主要源于PbI6八面體結(jié)構(gòu)的扭曲.在四種結(jié)構(gòu)體系中,CsPbI3顯示出最窄的直接帶隙、最小載流子有效質(zhì)量和較強的光吸收能力.這些結(jié)果可為進一步深入研究鈣鈦礦材料在太陽能電池領(lǐng)域的應用提供理論指導.

    鈣鈦礦結(jié)構(gòu)材料,結(jié)構(gòu)穩(wěn)定性,電子結(jié)構(gòu),第一性原理

    1 引 言

    鈣鈦礦太陽能電池是基于染料敏化太陽能電池發(fā)展起來的一種新型太陽能電池,具有光電轉(zhuǎn)換效率較高、制備簡單和易于規(guī)?;a(chǎn)等優(yōu)點,掀起了這類太陽能電池的研究熱潮[1?3].2009年,Kojima等[4]在染料敏化太陽能電池中首次將有機金屬鹵化物CH3NH3PbX3(X為Br,I)作為光吸收劑,制作的太陽能電池器件獲得了3.8%的光電轉(zhuǎn)換效率.之后,有機金屬鹵化物作為鈣鈦礦太陽能電池的關(guān)鍵材料引起了人們極大的關(guān)注.盡管在提高鈣鈦礦太陽能電池器件光電轉(zhuǎn)換效率方面的研究進展十分迅速[5],但是要實現(xiàn)鈣鈦礦太陽能電池的商業(yè)化,除了提高光電轉(zhuǎn)換效率外,更大的挑戰(zhàn)在于解決電池的穩(wěn)定性問題[6?8],尤其是作為光吸收層的鈣鈦礦材料的穩(wěn)定性.

    鈣鈦礦材料表現(xiàn)出較強的光吸收能力和較高的載流子遷移率,源于其特殊的有機無機雜化結(jié)構(gòu),分子式為ABX3,其中BX6構(gòu)成八面體,并相互接觸組成了三維結(jié)構(gòu)框架,帶電的A離子嵌入其中[9].其穩(wěn)定性能通過容差因子t進行判斷,其中RA,RB和RX分別表示A,B和X原子的有效半徑.目前對鈣鈦礦材料CH3NH3PbX3體系的研究較多,為了提高其結(jié)構(gòu)穩(wěn)定性和光電轉(zhuǎn)換性能,基于該結(jié)構(gòu)中A位置的陽離子替代的改性研究已經(jīng)展開[10?14].如選取離子半徑更大的替代形成的結(jié)構(gòu)晶格常數(shù)更大,帶隙稍小[10?12].盡管利用該鈣鈦礦材料制作的太陽能電池的效率達到20%以上,但是其結(jié)構(gòu)穩(wěn)定性較差,嚴重影響了實際應用.最近,部分實驗采用Cs+替代A位置MA+,發(fā)現(xiàn)其量子效率并沒有降低,同時鈣鈦礦材料的結(jié)構(gòu)穩(wěn)定性更高,顯示出極大的應用潛力[13,14].

    目前大量的研究主要集中在鈣鈦礦材料的合成和摻雜改性等工藝上[15],但是因其結(jié)構(gòu)的特殊性和A位陽離子與PbI6骨架之間較弱的相互作用,很難通過實驗直接探測,其微結(jié)構(gòu)與電子結(jié)構(gòu)之間的關(guān)聯(lián)性仍然需要進一步深入研究[16].采用第一性原理的方法可以獲得實驗研究無法獲取的有用信息,從材料設(shè)計的角度分析微結(jié)構(gòu)對光電性能的影響.最近有研究人員采用該方法研究了這類材料的電學性質(zhì)[17?21],但是有機陽離子的種類對結(jié)構(gòu)穩(wěn)定性和光電性質(zhì)影響的系統(tǒng)研究仍然沒有展開,其較強的光吸收能力和較高的載流子遷移率的物理機理仍然不清楚.為此,我們在A位選取不同的陽離子建立鈣鈦礦物理模型,即A位為Cs+,NH+4,MA+和FA+,采用基于密度泛函理論的第一性原理方法深入研究陽離子對鈣鈦礦結(jié)構(gòu)穩(wěn)定性、電子結(jié)構(gòu)和光學性質(zhì)的影響,為進一步提高鈣鈦礦材料的光吸收能力和光電轉(zhuǎn)換效率提供理論指導.

    2 物理模型與計算方法

    選取MAPbI3室溫下的四方相(空間群為I4/mcm)為研究對象,一個晶胞由4個MAPbI3分子組成,共48個原子,初始晶格常數(shù)采用實驗數(shù)據(jù)a=b=0.880 nm,c=1.269 nm[22].在此結(jié)構(gòu)的基礎(chǔ)上,分別用取代MAPbI3結(jié)構(gòu)中A位置的有機陽離子MA+,獲得了CsPbI3,NH4PbI3和FAPbI3的初始晶體結(jié)構(gòu),再分別對其進行幾何結(jié)構(gòu)優(yōu)化.A位置的四種陽離子分別為和FA+,其分子構(gòu)型如圖1所示.在MAPbI3和FAPbI3結(jié)構(gòu)中,考慮了MA+和FA+的不同取向?qū)w系總能量的影響,而電子結(jié)構(gòu)計算僅研究了對應結(jié)構(gòu)中體系總能量最低的結(jié)構(gòu)取向情形.

    采用基于密度泛函理論的平面波超軟贗勢方法[23],考慮陽離子基團和PbI骨架之間的相互作用,分別在廣義梯度近似(GGA)的PBE方案中考慮了TS色散修正[24],在局域密度近似(LDA)的CAPZ方案中考慮了OBS色散修正[25],獲得了相應的晶格常數(shù)、Pb—I鍵長、能帶結(jié)構(gòu)及光學性質(zhì)等.在描述離子實與價電子之間的相互作用時,選取的價電子組態(tài)分別為C:2s22p2,H:1s1,N:2s22p3,Pb:5d106s26p2,I:5s25p5,Cs:5s25p66s1.布里淵區(qū)k點網(wǎng)格均選為4×4×3[26]. 采用BFGS算法對所有模型進行幾何結(jié)構(gòu)優(yōu)化,平面波截斷能設(shè)置為400 eV,自洽收斂精度設(shè)置為5.0×10?5eV/atom,原子間的力場收斂精度設(shè)置為0.01 eV/nm,最大應力設(shè)置為0.2 GPa,最大位移設(shè)置為5×10?4nm.所有計算均由CASTEP軟件完成[27].

    圖1 (網(wǎng)刊彩色)鈣鈦礦APbI3結(jié)構(gòu)中A位陽離子的分子構(gòu)型(紫色球為Cs原子,藍色球為N原子,灰色球為C原子,白色球為H原子) (a)Cs+;(b)NH;(c)CH3NH;(d)為CH(NH2)Fig.1. (color online)Four molecule models of A cations in perovskite APbI3(purple ball is Cs atom,blue balls are N atoms,gray balls are C atoms,and white balls are H atoms):(a)Cs+;(b)NH;(c)CH3NH;(d)CH(NH2).

    3 結(jié)果與討論

    3.1 結(jié)構(gòu)穩(wěn)定性

    為了獲得結(jié)構(gòu)合理且總能量較低的穩(wěn)定晶體結(jié)構(gòu),對建立的所有結(jié)構(gòu)進行全局幾何結(jié)構(gòu)優(yōu)化,即不固定原子位置、晶格常數(shù)以及各棱之間夾角α,β,γ,優(yōu)化后的晶體結(jié)構(gòu)如圖2所示.可以看出,結(jié)構(gòu)優(yōu)化后NH4PbI3,MAPbI3,FAPbI3中PbI6八面體在一定程度上發(fā)生了較明顯的傾斜和旋轉(zhuǎn),而在CsPbI3中PbI6八面體的傾斜和旋轉(zhuǎn)不明顯.表1為采用基于TS色散修正的GGA(PBE)方法和OBS色散修正的LDA(CAPZ)方法對4種鈣鈦礦結(jié)構(gòu)進行全局結(jié)構(gòu)優(yōu)化后,獲得的晶體結(jié)構(gòu)參數(shù)、禁帶寬度和鍵長等結(jié)果.

    由于Cs原子與周圍PbI骨架之間是共價鍵作用,不存在較弱的范德瓦耳斯力,因此對CsPbI3體系進行幾何結(jié)構(gòu)優(yōu)化時沒有考慮色散修正.眾所周知,LDA計算的晶格常數(shù)比實驗值偏小,而GGA計算的晶格常數(shù)比實驗值偏大.例如,在MAPbI3結(jié)構(gòu)中分別采用GGA和LDA方法,計算出的晶胞體積分別比實驗值大7.54%和小7.05%.由此可見計算方法的選擇在一定程度上影響晶體結(jié)構(gòu)和電子結(jié)構(gòu)的結(jié)果.總體來說,除了Cs+外,晶胞體積會隨著有機陽離子FA+半徑RFA+=0.279 nm)的增大而增大.由于Pb2+和I?離子的半徑分別為0.133 nm和0.203 nm,穩(wěn)定的鈣鈦礦結(jié)構(gòu)中占據(jù)A位置的粒子半徑RA應該滿足0.130 nm

    表1列出了采用GGA(PBE)和LDA(CAPZ)兩種方法獲得的Pb—I鍵長以及H和I之間的距離,這里的“數(shù)量”為H和I間距小于0.28 nm的近鄰數(shù)目.可以看出,與其他體系相比,在CsPbI3中Pb—I鍵長的平均值最小,且Pb—I鍵長的變化范圍最小,這表明PbI6構(gòu)成的八面體扭曲及旋轉(zhuǎn)程度最小.隨著A位陽離子(除Cs+)有效半徑的增大,在GGA算法中,Pb—I的平均鍵長從0.325 nm增大至0.329 nm,H和I之間的近鄰數(shù)目從7增大至11;而在LDA算法中,Pb—I的平均鍵長從0.316 nm增大到0.320 nm,H和I之間的近鄰數(shù)目從17增大到21.這說明陽離子A的半徑越大,與無機骨架之間的作用力將越強.在MAPbI3中Pb—I鍵長的變化范圍分別為0.046 nm(GGA)和0.064 nm(LDA),其對應的PbI6八面體電偶極矩分別為0.23D和0.28D;在FAPbI3體系中Pb—I鍵長的變化范圍分別為0.086 nm(GGA)和0.072 nm(LDA),其對應的PbI6八面體電偶極矩分別為0.32D和0.29D,較高的電偶極矩有利于光生電子和空穴的分離.

    隨后對MAPbI3和FAPbI3結(jié)構(gòu)中MA+和FA+陽離子的取向展開了研究.發(fā)現(xiàn)當MA+陽離子中C—N鍵的方向平行于xoy平面,取向如圖2(c)所示時,MAPbI3具有最低的能量.當FA+陽離子中N—C—N鍵分子所在平面垂直于xoy平面,取向如圖2(d)所示時,FAPbI3具有最低的能量.MAPbI3和FAPbI3結(jié)構(gòu)中MA+和FA+陽離子取向及轉(zhuǎn)動引起體系總能量的變化小于0.82 eV和0.97 eV,這表明有機陽離子與周圍骨架之間的靜電作用并不是很強,因此有機陽離子的取向容易受到制備條件和外界環(huán)境的影響,導致其結(jié)構(gòu)不穩(wěn)定而容易發(fā)生相變.

    表1 采用GGA和LDA方法對APbI3(A為Cs+,NH+4,MA+,FA+)進行結(jié)構(gòu)優(yōu)化后的晶格常數(shù)、禁帶寬度、Pb—I鍵長,以及H和I的近鄰數(shù)目Table 1.The calculated crystal parameters,band gap,Pb—I bond length and the nearest neighbor number of H and I atoms for perovskite APbI3(A denotes Cs+,NH+4,MA+,FA+)optimized with GGA and LDA.

    圖2 (網(wǎng)刊彩色)采用GGA方法進行結(jié)構(gòu)優(yōu)化后APbI3(A為Cs+,NH,MA+,FA+)的晶體結(jié)構(gòu),其中紫色球為Cs原子,藍色球為N原子,白色球為H原子,淺灰色球為C原子,深灰色為Pb原子,褐色球為I原子Fig.2.(color online)The geometrical structures of perovskite APbI3(A denotes Cs+,NH,MA+,FA+)optimized with GGA.Purple balls are Cs atoms,blue balls are N atoms,white balls are H atoms,light gray balls are C atoms,dark gray balls are Pb atoms,and brown balls are I atoms.

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

    材料的電子結(jié)構(gòu)決定光生載流子的產(chǎn)生、復合及輸運等性質(zhì).采用APbI3(A為Cs+,NH+4,MA+,FA+)全局結(jié)構(gòu)優(yōu)化后的幾何結(jié)構(gòu)進行能帶結(jié)構(gòu)和態(tài)密度計算.由于GGA方法計算的晶格常數(shù)更接近實驗值,后續(xù)所有的能帶結(jié)構(gòu)和光學性質(zhì)計算均采用GGA方法.表1列出了采用GGA和LDA方法獲得的禁帶寬度,總體上加入了OBS色散修正的LDA方法獲得的禁帶寬度小于加入了TS色散修正的GGA方法,兩種方法計算得到的禁帶寬度大小順序是一致的.由此可見,A位置陽離子種類對其電子結(jié)構(gòu)的影響較大,計算結(jié)果表明晶胞體積越小,禁帶寬度越大.例如兩種方法優(yōu)化后的NH4PbI3具有最小的晶胞尺寸,卻具有最大的禁帶寬度,比FAPbI3大約0.1 eV,因而其可見光的吸收范圍稍小.而CsPbI3的禁帶寬度分別為1.52 eV(GGA)和0.96 eV(LDA),與其他體系相比均為最小.

    四種鈣鈦礦材料的能帶結(jié)構(gòu)如圖3所示,從圖中可以看出,四種鈣鈦礦材料具有相似的能帶結(jié)構(gòu),其中CsPbI3和FAPbI3均為高對稱G(0,0,0)點的直接帶隙半導體,禁帶寬度分別為1.52 eV和1.72 eV,CsPbI3的禁帶寬度計算值與實驗值相當符合,僅相差0.05 eV[28],FAPbI3的禁帶寬度計算值與實驗值相差0.25 eV[32].而NH4PbI3和MAPbI3表現(xiàn)為間接帶隙半導體,禁帶寬度分別為1.83和1.79 eV,他們的導帶底在G點左側(cè),能量值均比G點能量值高0.03 eV,價帶頂分別在高對稱點Z(0,0,0.5)和G(0,0,0)處,其中MAPbI3禁帶寬度計算值與實驗值比較符合(表1).由于這些能量值相差不大,可近似地認為NH4PbI3和MAPbI3是在G點的直接帶隙半導體.這里計算的MAPbI3和FAPbI3禁帶寬度與Lee等[33]獲得的結(jié)果非常接近,差值分別為0.098和0.099 eV.這可能是由于在計算過程中沒有考慮自旋耦合等因素,但是這并不影響對其電子結(jié)構(gòu)的分析.

    為了進一步深入分析APbI3(A為Cs+,NH+4,MA+,FA+)的電子結(jié)構(gòu),分別繪出了他們的態(tài)密度、最高占據(jù)分子軌道(HOMO)和最低未占據(jù)分子軌道(LUMO)圖,如圖4所示.可以看出,所有體系的導帶底由Pb 6p軌道和部分I 5p軌道雜化而成,價帶頂由I 5p軌道組成.圖4中LUMO和HOMO圖更直觀地顯示了這些導帶底和價帶頂附近的軌道分布.可以看出鈣鈦礦結(jié)構(gòu)中Cs+和FA+對價帶頂附近的貢獻比NH+4和MA+的貢獻稍大,其中CsPbI3在其價帶頂附近(?1.13 eV)存在一個較弱的態(tài)密度峰,而FAPbI3在其價帶頂和導帶底附近(?0.60 eV和3.31 eV)有明顯的態(tài)密度峰.總體上來說,A位陽離子對費米能級附近能帶的直接貢獻比較弱,僅起到電荷平衡的作用.

    鈣鈦礦材料具有優(yōu)異的電荷傳輸能力,為此計算了APbI3(A為Cs+,NH+4,MA+,FA+)光生載流子的有效質(zhì)量.載流子的有效質(zhì)量為

    式中a=?2E(k)/?k2,為能帶結(jié)構(gòu)中E(k)的二次函數(shù)方程中的二次項常系數(shù),h為普朗克常數(shù).通過第一性原理計算,獲得導帶底和價帶頂?shù)腅(k)和k的關(guān)系,之后在其極值點附近用二次函數(shù)方程進行擬合,即可得到二次項常系數(shù)a.擬合曲線(虛線)如圖3所示,其中,Eg為禁帶寬度,ae和ah分別為導帶底最低點和價帶頂最高點二次函數(shù)方程擬合的二次項系數(shù).再通過(1)式求出光生載流子的有效質(zhì)量.從圖4可以看出,CsPbI3和FAPbI3體系中導帶底和價帶頂擬合得到的二次項常系數(shù)a比NH4PbI3和MAPbI3大很多,即電子和空穴的有效質(zhì)量較小,受外場的影響較大.

    圖3 鈣鈦礦材料APbI3(A為Cs+,NH,MA+,FA+)導帶底和價帶頂?shù)哪軒€(實線),以及導帶底的最低位置和價帶頂最高位置附近的二次函數(shù)擬合曲線(虛線) (a)CsPbI3,Eg=1.52 eV,μ=0.08;(b)NH4PbI3,Eg=1.83 eV,μ=0.22;(c)CH3NH3PbI3,Eg=1.79 eV,μ=0.18;(d)CH(NH2)2PbI3,Eg=1.72 eV,μ=0.11Fig.3.The top curves of valence band and the bottom curves of conduction band for perovskite APbI3(A denotes Cs+,NH,MA+,FA+)(solid lines),and quadratic fitting curves(dash lines)at bottom point of conduction band and top point of valence band:(a)CsPbI3,Eg=1.52 eV,μ=0.08;(b)NH4PbI3,Eg=1.83 eV,μ=0.22;(c)CH3NH3PbI3,Eg=1.79 eV,μ=0.18;(d)CH(NH2)2PbI3,Eg=1.72 eV,μ=0.11.

    圖4 (網(wǎng)刊彩色)采用GGA方法計算得到的APbI3(A為Cs+,NH,MA+,FA+)的總態(tài)密度和分態(tài)密度,內(nèi)嵌圖為LUMO圖(綠色)和HOMO圖(藍色) (a)CsPbI3;(b)NH4PbI3;(c)CH3NH3PbI3;(d)CH(NH2)2PbI3Fig.4.(color online)The total and partial density of states of perovskite APbI3(A denotes Cs+,NH,MA+,FA+).The insets are LUMO(green part)and HOMO(blue part):(a)CsPbI3;(b)NH4PbI3;(c)CH3NH3PbI3;(d)CH(NH2)2PbI3.

    3.3 光學性質(zhì)

    光學性質(zhì)是半導體物理性質(zhì)最重要的方面之一.計算電子結(jié)構(gòu)時無論是帶間還是帶內(nèi)躍遷頻率都遠超過聲子頻率,而且使用的方法是單電子近似方法,故僅考慮電子激發(fā).從量子力學的觀點看,帶間躍遷光吸收過程是電子在輻射電磁場微擾的作用下從低能占據(jù)態(tài)到高能未占據(jù)態(tài)之間的躍遷過程.根據(jù)費米黃金定律,介電函數(shù)虛部的計算可利用電偶極子近似進行,從直接躍遷概率的定義推導出介電函數(shù)的虛部ε2為

    式中u為入射電場的極化方向量;e為電子電量;?為倒空間體積;ε0為真空介電常數(shù);C,V分別表示導帶和價帶;K為倒格矢;為動量躍遷矩陣;分別為導帶和價帶上的本征能級.根據(jù)Kramers-Kronig變換,可以利用介電函數(shù)的虛部得到其實部.利用吸收系數(shù)η與復折射率虛部k之間的關(guān)系可以得到吸收光譜η=2kω/c,其中c為真空中的光速,ω為圓頻率.后文將采用第一性原理的計算方法獲得晶體材料的吸收光譜等信息.

    APbI3特殊的晶體結(jié)構(gòu)使其在可見光區(qū)域內(nèi)可能具有優(yōu)異的光吸收特性.圖5為采用GGA方法計算得到的APbI3(A為Cs+,MA+,FA+)的光吸收譜.可以看出,這四種鈣鈦礦材料在E//x和E//y方向存在兩個較明顯的光吸收峰,而且這兩個方向的光吸收曲線幾乎重合,但是與E//z方向光吸收峰明顯不同.在CsPbI3和FAPbI3體系中,E//x,E//y和E//z的光吸收峰均處于相同的波長位置,僅光吸收系數(shù)的峰高不相同.而在NH4PbI3和MAPbI3體系中光吸收峰的位置有明顯的差異.總體說來,這四種鈣鈦礦材料的光吸收都表現(xiàn)出一定的各向異性,這是由晶體結(jié)構(gòu)的不對稱性決定的.

    圖5 (網(wǎng)刊彩色)鈣鈦礦材料APbI3(A為Cs+,NH,MA+,FA+)的光吸收譜 (a)CsPbI3;(b)NH4PbI3;(c)MAPbI3;(d)FAPbI3Fig.5. (color online)Absorption spectra of perovskite APbI3(A denotes Cs+,NH,MA+,FA+):(a)CsPbI3;(b)NH4PbI3;(c)MAPbI3;(d)FAPbI3.

    在CsPbI3,NH4 PbI3,MAPbI3和FAPbI3四種體系中最大波長的吸收峰分別位于373 nm(E//x,E//y和E//z),369 nm(E//x,E//y),376 nm(E//x,E//y)和319 nm(E//x,E//y和E//z).以MAPbI3為例,計算出的光吸收峰與通過紫外-可見光譜儀測出的吸收峰390 nm情形非常接近[21].進一步觀察光吸收譜可以發(fā)現(xiàn),在CsPbI3和FAPbI3體系中E//z方向的光吸收峰非常高,與E//x和E//y方向相比,其光吸收邊向長波方向有較大的擴展.CsPbI3的吸收帶邊在850 nm附近,而NH4PbI3,MAPbI3和FAPbI3的吸收帶邊分別位于670,690,720 nm附近,這都是由其自身的禁帶寬度決定的,對應著吸收光子的電子從價帶頂躍遷到導帶底的電子躍遷過程.總體來說,CsPbI3體系在長波范圍內(nèi)的可見光和近紅外區(qū)域均有較強的光吸收能力.

    4 結(jié) 論

    采用平面波超軟贗勢方法研究了APbI3(A為Cs+,NH,MA+,FA+)結(jié)構(gòu)的穩(wěn)定性、電子結(jié)構(gòu)及光學性質(zhì).考慮到A位有機陽離子與PbI基質(zhì)之間的相互作用較弱,計算時引入了相應的色散修正方案.結(jié)果顯示受陽離子種類和尺寸的影響,PbI基體骨架發(fā)生不同程度的扭曲,MAPbI3和FAPbI3中PbI6八面體顯示出較大的電偶極矩,這有利于光生電子和空穴的分離.四種體系中費米能級附近的能帶比較相似,即價帶頂均由I 5p軌道組成,導帶底由Pb 6p軌道和部分I 5p軌道雜化而成.電子結(jié)構(gòu)和光學性質(zhì)的差異與PbI6八面體結(jié)構(gòu)的扭曲程度密切相關(guān).通過比較可以發(fā)現(xiàn),CsPbI3顯示出最窄的直接帶隙、最小載流子有效質(zhì)量和較強的光吸收能力.由于CsPbI3表現(xiàn)出優(yōu)異的性質(zhì),將在隨后的研究工作中進一步探討混合鈣鈦礦CsxA1?xPbI3(A為NH+4,MA+,FA+)體系的穩(wěn)定性和物理性質(zhì).

    [1]Mei A Y,Li X,Liu L F,Ku Z L,Liu T F,Rong Y G,Xu M,Hu M,Chen J Z,Yang Y,Gr?tzel M,Han H W 2014Science345 295

    [2]Kim H S,Lee C R,Im J H,Lee K B,Moehl T,Marchioro A,Moon S J,Baker R H,Yum J H,Moser J E,Gr?tzel M,Park N G 2012Sci.Rep.2 591

    [3]Wang F Z,Tan Z A,Dai S Y,Li Y F 2015Acta Phys.Sin.64 038401(in Chinese)[王福芝,譚占鰲,戴松元,李永舫2015物理學報64 038401]

    [4]Kojima A,Teshima K,Shirai Y,Miyasaka T 2009J.Am.Chem.Soc.131 6050

    [5]Yang W S,Noh J H,Jeon N J,Kim Y C,Ryu S,Seo J,Seok S II 2015Science348 1234

    [6]Zhang D F,Zheng L L,Ma Y Z,Wang S F,Bian Z Q,Huang C H,Gong Q H,Xiao L X 2015Acta Phys.Sin.64 038803(in Chinese)[張丹霏,鄭靈靈,馬英壯,王樹峰,卞祖強,黃春輝,龔旗煌,肖立新2015物理學報64 038803]

    [7]Cappel U B,Daeneke T,Bach U 2012Nano Lett.12 4925

    [8]Liu M Z,Johnston M B,Snaith H J 2013Nature501 395

    [9]Knop O,Wasylishen R E,White M A,Oort M J M V 1990Can.J.Chem.68 412

    [10]Lee J W,Seol D J,Cho A N 2014Adv.Mater.26 4991

    [11]Zhou Y Y,Yang M J,Pang S P,Zhu K,Padture N P 2016J.Am.Chem.Soc.138 5535

    [12]Pang S P,Hu H,Zhang J L,Lv S L,Yu Y M,Wei F,Qin T S,Xu H X,Liu Z L,Cui G L 2014Chem.Mater.26 1485

    [13]Choi H,Jeong J,Kim H B,Kim S,Walker B,Kim G H,Kim J Y 2014Nano Energy7 80

    [14]Saliba M,Matsui T,Seo J Y,Domanski K,Correa-Baena J P,Nazeeruddin M K,Zakeeruddin S M,Tress W,Abate A,Hagfeldt A,Gr?tzel M 2016Energy Environ.Sci.9 1989

    [15]Baikie T,Fang Y A,Kadro J M,Schreyer M,Wei F X,Mhaisalkar S G,Gr?tzel M,White T J 2013J.Mater.Chem.A1 5628

    [16]Motta C,Mellouhi F E,Kais S,Tabet N,Alharbi F,Sanvito S 2015Nat.Commun.6 7026

    [17]Filippetti A,Mattoni A 2014Phys.Rev.B89 12503

    [18]Mosconi E,Amat A,Nazeeruddin M K,Gr?tzel M,De Angelis F 2013J.Phys.Chem.C117 13902

    [19]Geng W,Zhang L,Zhang Y N,Lau W M,Liu L M 2014J.Phys.Chem.C118 19565

    [20]Wang Y,Gould T,Dobson J F,Zhang H M,Yang H G,Yao X D,Zhao H J 2014J.Phys.Chem.Chem.Phys.16 1424

    [21]Umari P,Mosconi E,De Angelis F 2014Sci.Rep.4 4467

    [22]Kawamura Y,Mashiyama H,Hasebe K 2002J.Phys.Soc.Jpn.71 1694

    [23]Vanderbilt D 1990Phys.Rev.B41 7892

    [24]Tkatchenko A,Scheffler M 2009Phys.Rev.Lett.102 073005

    [25]Ortmann F,Bechstedt F,Schmidt W G 2006Phys.Rev.B73 205101

    [26]Monkhorst H J,Pack J D 1976Phys.Rev.B13 5188

    [27]Segall M D,Lindan P J D,Probert M J,Pickard C J,Hasnip P J,Clark S J,Payne M C 2002J.Phys:Condens.Matter.14 2717

    [28]Chung L,Lee B,He J Q,Chang R P H,Kanatzidis M G 2012Nature485 486

    [29]Gao X,Uehara K,Klug D D,Patchkovskii S,Tse J S,Tritt T M 2005Phys.Rev.B72 125202

    [30]Tanaka K,Takahashi T,Ban T,Kondo T,Uchida K,Miura N 2003Solid State Commun.127 619

    [31]Schulz P E,Edri E,Kirmayer S,Hodes G,Cahen D,Kahn A 2014Energy Environ.Sci.7 1377

    [32]Jeon N J,Noh J H,Yang W S,Kim Y C,Ryu S 2015Nature517 476

    [33]Lee C,Hong J,Stroppa A,Whangbo M H,Shim J H 2015RSC Advances5 78701

    PACS:71.15.Mb,71.20.–b,81.07.Pr,88.40.H– DOI:10.7498/aps.66.057103

    Theoretical study on the stability and photoelectric properties of APbI3perovskite?

    Liu Na1)Wei Yang1)Ma Xin-Guo1)2)?Zhu Lin1)Xu Guo-Wang1)2)Chu Liang3)Huang Chu-Yun1)2)?

    1)(School of Science,Hubei University of Technology,Wuhan 430068,China)
    2)(Hubei Collaborative Innovation Center for High-Efficiency Utilization of Solar Energy,Hubei University of Technology,Wuhan 430068,China)
    3)(School of Science,Nanjing University of Posts and Telecommunications,Nanjing 210046,China)

    29 July 2016;revised manuscript

    4 December 2016)

    The rapid development of organic-inorganic hybrid perovskite solar cells has recently attracted the worldwide attention because their power conversion efficiency has risen from 4%to higher than 20%within just six years.It is well known that the perovskite materials with APbI3crystal structure have a 3D framework of corner-sharing PbI6octahedra,in which each Pb atom bonds with six I atoms,and the A cations fill in the octahedral interstices.At present,a lot of researches have focused on the synthesis and doping modification of perovskite materials.However,it is hard to detect directly the weak interactions between A cations and PbI6skeleton in the APbI3crystal structure through experiments,which have effect on the structural stability and electronic properties.To provide a full understanding of the interplay among size,structure,and organic/inorganic interactions,the stability,electronic structures and optical properties of APbI3were investigated by the plane-wave ultra soft pseudo potentials.Two dispersion corrections were taken into account in the weak interactions between A cations and PbI6skeleton in the APbI3crystal structure,respectively.The results show that the type and size of cations affect the distortion of PbI framework,indicating that the larger the radius of the A cation is,the stronger the interaction between the A cation and the PbI framework is.Further,it is identified that after geometry relaxation,the orientation of A cations(A denotesis easy to change,and the PbI frameworks present structural distortion.CsPbI3is more stable energetically than other three kinds of perovskite materials.For the PbI6octahedra,the large dipole moments of 0.23D and 0.32D for the generalized-gradient approximation method or 0.28D and 0.29D for the local-density approximation method are also present in MAPbI3and FAPbI3,respectively.In addition,the energy band structures,which affect the generation and migration of photon-generated carriers and optical properties,will alter with the structural distortion of PbI frameworks.By analyzing the energy band structures and corresponding density of states,we find that four systems have similar band structures near the Fermi energy,namely,the top of valance band is mainly contributed by I 5p orbitals,while the bottom of conduction band is dominated by Pb 6p orbitals and partly contributed by I 5p orbitals.A little difference of their electronic structures and optical absorption spectra originates from the distortion of PbI6octahedra in APbI3crystal structures.It is noted that the contribution of the ions Cs+and FA+on the top of valance band is slightly larger than that of the ions NH+4and MA+.Compared with other three kinds of perovskite materials,CsPbI3presents the narrowest direct band gap,the lowest effective carrier mass and excellent visible-light and infrared absorption.The results may provide some theoretical guidance for further research on perovskite materials in the application of solar cells.

    perovskite material,structural stability,electronic structure,first-principles

    PACS:71.15.Mb,71.20.–b,81.07.Pr,88.40.H–

    10.7498/aps.66.057103

    ?國家自然科學基金(批準號:51472081)、湖北工業(yè)大學高層次人才啟動基金(批準號:GCRC13014)、綠色工業(yè)引領(lǐng)計劃(批準號:YXQN2016005)和湖北省協(xié)同創(chuàng)新中心開放基金(批準號:HBSKFZD2014003,HBSKFZD2014011,HBSKFZD2015004)資助的課題.

    ?通信作者.E-mail:maxg2013@sohu.com

    ?通信作者.E-mail:chuyunh@163.com

    *Project supported by the National Natural Science Foundation of China(Grant No.51472081),the Foundation of Hubei University of Technology for High-Level Talents(Grant No.GCRC13014),the Leading Plan of Green Industry(Grant No.YXQN2016005),and the Development Founds of Hubei Collaborative Innovation Center(Grant Nos.HBSKFZD2014003,HBSKFZD2014011,HBSKFZD2015004).

    ?Corresponding author.E-mail:maxg2013@sohu.com

    ?Corresponding author.E-mail:chuyunh@163.com

    猜你喜歡
    光吸收價帶導帶
    外加電場和磁場對外爾半金屬能帶結(jié)構(gòu)的調(diào)控*
    關(guān)于某型系列產(chǎn)品“上下導帶二次收緊”合并的工藝技術(shù)
    炮射彈藥滑動導帶環(huán)結(jié)構(gòu)設(shè)計
    導帶式數(shù)碼噴射印花機運行過程中常見疵病
    打捆機導帶槽裝置及捆帶盤脹縮裝置的改進
    重型機械(2019年2期)2019-04-28 11:52:04
    交換場和非共振光對單層MoS2能帶結(jié)構(gòu)的調(diào)控*
    Modification of CaO-based sorbents prepared from calcium acetate for CO2 capture at high temperature☆
    多元稀土硼化物Ce1-xNdxB6的制備及性能研究*
    功能材料(2016年1期)2016-05-17 03:38:24
    鑲嵌納米晶硅的氧化硅薄膜微觀結(jié)構(gòu)調(diào)整及其光吸收特性
    Photoabsorp tion of Na n C lusters and Vo lum e P lasm on:Theory Rev iew
    午夜福利在线观看吧| 欧美中文综合在线视频| 成人国产综合亚洲| av欧美777| 久久久久久大精品| 国产极品粉嫩免费观看在线| 亚洲色图综合在线观看| 天天躁夜夜躁狠狠躁躁| 长腿黑丝高跟| 给我免费播放毛片高清在线观看| 精品日产1卡2卡| 真人做人爱边吃奶动态| 黑人巨大精品欧美一区二区蜜桃| 变态另类丝袜制服| 国产精华一区二区三区| 校园春色视频在线观看| 777久久人妻少妇嫩草av网站| 国产av在哪里看| av有码第一页| 波多野结衣一区麻豆| 国产精品 欧美亚洲| 亚洲最大成人中文| 一a级毛片在线观看| 香蕉丝袜av| 日本一区二区免费在线视频| 欧美最黄视频在线播放免费| 女人被躁到高潮嗷嗷叫费观| 91精品三级在线观看| 18美女黄网站色大片免费观看| 欧美一级毛片孕妇| 美女国产高潮福利片在线看| 亚洲片人在线观看| 久久中文字幕一级| 在线观看免费午夜福利视频| 人人澡人人妻人| 变态另类成人亚洲欧美熟女 | 丝袜人妻中文字幕| 欧美另类亚洲清纯唯美| 女人爽到高潮嗷嗷叫在线视频| 18禁观看日本| 日本精品一区二区三区蜜桃| 中亚洲国语对白在线视频| 大码成人一级视频| 91九色精品人成在线观看| 动漫黄色视频在线观看| 日韩欧美一区二区三区在线观看| 久久精品国产亚洲av高清一级| 两人在一起打扑克的视频| 丁香六月欧美| 国产一级毛片七仙女欲春2 | 国产成人影院久久av| 欧美日韩亚洲国产一区二区在线观看| 91精品国产国语对白视频| 999久久久精品免费观看国产| 亚洲视频免费观看视频| 久久亚洲真实| 成熟少妇高潮喷水视频| 亚洲最大成人中文| 在线视频色国产色| 国产精品精品国产色婷婷| 国产极品粉嫩免费观看在线| 啦啦啦韩国在线观看视频| 精品一区二区三区视频在线观看免费| 精品国产美女av久久久久小说| 亚洲精华国产精华精| 成人亚洲精品av一区二区| 99国产精品一区二区三区| 国产三级在线视频| 免费高清在线观看日韩| 免费在线观看黄色视频的| 可以在线观看毛片的网站| 久久精品91蜜桃| 亚洲视频免费观看视频| 一边摸一边抽搐一进一小说| 国产精品久久电影中文字幕| 少妇熟女aⅴ在线视频| 久久国产乱子伦精品免费另类| 精品久久久久久久毛片微露脸| 一进一出抽搐动态| 黄色毛片三级朝国网站| 少妇粗大呻吟视频| 久久精品影院6| tocl精华| 国产99久久九九免费精品| 免费在线观看影片大全网站| 深夜精品福利| 成人手机av| 制服人妻中文乱码| 一级a爱视频在线免费观看| 男女下面进入的视频免费午夜 | 欧美人与性动交α欧美精品济南到| 我的亚洲天堂| 啦啦啦观看免费观看视频高清 | 欧美色欧美亚洲另类二区 | 18美女黄网站色大片免费观看| 国产精品永久免费网站| 变态另类成人亚洲欧美熟女 | 久久久久国产一级毛片高清牌| 黄色丝袜av网址大全| 亚洲国产精品合色在线| 国产又色又爽无遮挡免费看| 久久久久久久精品吃奶| www.精华液| 国语自产精品视频在线第100页| 午夜a级毛片| a级毛片在线看网站| 欧洲精品卡2卡3卡4卡5卡区| 97人妻天天添夜夜摸| 黄色女人牲交| 亚洲人成电影观看| 欧美久久黑人一区二区| 妹子高潮喷水视频| 久久久水蜜桃国产精品网| av福利片在线| 亚洲av成人一区二区三| 亚洲美女黄片视频| 日日干狠狠操夜夜爽| 亚洲男人的天堂狠狠| 亚洲专区国产一区二区| 少妇被粗大的猛进出69影院| 亚洲精品久久国产高清桃花| 最新美女视频免费是黄的| 不卡一级毛片| 亚洲精品美女久久av网站| 国产欧美日韩综合在线一区二区| 中文字幕高清在线视频| 亚洲精品一区av在线观看| 久久中文字幕人妻熟女| 久久精品国产99精品国产亚洲性色 | 免费一级毛片在线播放高清视频 | 亚洲av电影不卡..在线观看| 777久久人妻少妇嫩草av网站| 欧美午夜高清在线| 国产成人影院久久av| 亚洲欧美日韩高清在线视频| 91字幕亚洲| 久久久精品欧美日韩精品| 国产精品亚洲av一区麻豆| 又黄又爽又免费观看的视频| 少妇被粗大的猛进出69影院| 精品欧美国产一区二区三| 99国产精品99久久久久| 亚洲精品在线观看二区| 国产成人精品久久二区二区91| 免费久久久久久久精品成人欧美视频| 黄色成人免费大全| 久久中文看片网| 欧美日韩福利视频一区二区| 亚洲av成人一区二区三| 美女午夜性视频免费| 麻豆av在线久日| 精品国内亚洲2022精品成人| 乱人伦中国视频| 成在线人永久免费视频| 99精品久久久久人妻精品| 国产片内射在线| 欧美激情久久久久久爽电影 | 丁香六月欧美| 91成人精品电影| 国产高清有码在线观看视频 | 欧美成人一区二区免费高清观看 | 一夜夜www| 久久伊人香网站| 亚洲精品中文字幕在线视频| 亚洲av电影在线进入| 大码成人一级视频| 国产亚洲精品久久久久久毛片| www.999成人在线观看| 非洲黑人性xxxx精品又粗又长| 成人永久免费在线观看视频| av在线播放免费不卡| 欧美乱妇无乱码| 在线国产一区二区在线| 午夜影院日韩av| 国产成人精品无人区| 国产精品永久免费网站| ponron亚洲| 精品国产乱码久久久久久男人| 多毛熟女@视频| 亚洲精品中文字幕在线视频| 制服丝袜大香蕉在线| 国产精品 欧美亚洲| 少妇粗大呻吟视频| 免费看美女性在线毛片视频| 国产精品久久电影中文字幕| 黄色 视频免费看| 嫩草影院精品99| 色综合欧美亚洲国产小说| 国产亚洲av高清不卡| 亚洲精品国产精品久久久不卡| 国产单亲对白刺激| 在线观看免费视频日本深夜| 久久久久久国产a免费观看| 狂野欧美激情性xxxx| 亚洲成人免费电影在线观看| or卡值多少钱| 亚洲男人的天堂狠狠| 丝袜在线中文字幕| 亚洲九九香蕉| 国产精品一区二区在线不卡| 999久久久国产精品视频| 国产高清有码在线观看视频 | 禁无遮挡网站| 一区二区三区国产精品乱码| 国产xxxxx性猛交| 久久久久国内视频| 淫妇啪啪啪对白视频| 免费看a级黄色片| 亚洲精华国产精华精| 久久 成人 亚洲| 99re在线观看精品视频| 两性夫妻黄色片| 久久精品国产综合久久久| 男人舔女人下体高潮全视频| 最近最新免费中文字幕在线| 99国产精品99久久久久| 一进一出抽搐gif免费好疼| 国产一区二区激情短视频| 精品久久久久久成人av| 99久久久亚洲精品蜜臀av| 丝袜美腿诱惑在线| 91大片在线观看| 波多野结衣高清无吗| 老鸭窝网址在线观看| 国产精品野战在线观看| 男女下面进入的视频免费午夜 | 香蕉丝袜av| 国产片内射在线| 国产男靠女视频免费网站| 亚洲精品久久国产高清桃花| 91在线观看av| 91麻豆精品激情在线观看国产| aaaaa片日本免费| 在线观看66精品国产| or卡值多少钱| av天堂久久9| 涩涩av久久男人的天堂| 高清黄色对白视频在线免费看| 国产在线精品亚洲第一网站| 国产三级在线视频| 午夜福利视频1000在线观看 | 国产精华一区二区三区| 亚洲专区字幕在线| 在线免费观看的www视频| 成人国语在线视频| 亚洲 欧美 日韩 在线 免费| 日韩欧美免费精品| 精品第一国产精品| 色在线成人网| 久久九九热精品免费| 搞女人的毛片| 国产黄a三级三级三级人| 欧美成人免费av一区二区三区| 不卡一级毛片| 大型av网站在线播放| 国产成人免费无遮挡视频| 男女做爰动态图高潮gif福利片 | 亚洲国产欧美日韩在线播放| 亚洲欧美一区二区三区黑人| 久久青草综合色| 久久精品aⅴ一区二区三区四区| 亚洲欧美激情综合另类| 久久精品91蜜桃| 性色av乱码一区二区三区2| 精品国产乱子伦一区二区三区| 欧美日本亚洲视频在线播放| 午夜精品在线福利| 亚洲国产精品合色在线| 亚洲av电影在线进入| 久久人妻av系列| 搡老岳熟女国产| 操出白浆在线播放| 侵犯人妻中文字幕一二三四区| 国产av精品麻豆| 国产成人av激情在线播放| 一a级毛片在线观看| 国产精品秋霞免费鲁丝片| 又紧又爽又黄一区二区| 国产精品av久久久久免费| 人人澡人人妻人| 在线观看一区二区三区| aaaaa片日本免费| 国产一区二区三区视频了| 久久久久久久精品吃奶| 精品国产国语对白av| 国产精品一区二区三区四区久久 | 怎么达到女性高潮| 男女做爰动态图高潮gif福利片 | www.自偷自拍.com| 禁无遮挡网站| 午夜激情av网站| 黄色视频,在线免费观看| xxx96com| 久久亚洲精品不卡| 97超级碰碰碰精品色视频在线观看| 亚洲成人免费电影在线观看| av视频在线观看入口| 国产1区2区3区精品| 久久中文看片网| 欧美日韩亚洲国产一区二区在线观看| 国产av在哪里看| 精品欧美国产一区二区三| 精品福利观看| 丝袜美腿诱惑在线| 看黄色毛片网站| 美女高潮喷水抽搐中文字幕| 亚洲av电影不卡..在线观看| 国产野战对白在线观看| 一级作爱视频免费观看| 中文亚洲av片在线观看爽| 搡老妇女老女人老熟妇| 亚洲国产欧美网| 亚洲av成人av| 日韩欧美国产在线观看| 亚洲国产看品久久| 日本vs欧美在线观看视频| 91精品国产国语对白视频| 午夜免费成人在线视频| 搡老妇女老女人老熟妇| 首页视频小说图片口味搜索| 国产精品99久久99久久久不卡| 国产主播在线观看一区二区| 亚洲第一青青草原| 给我免费播放毛片高清在线观看| 在线播放国产精品三级| 亚洲专区国产一区二区| 午夜福利一区二区在线看| 欧美乱妇无乱码| 日韩大尺度精品在线看网址 | 日韩有码中文字幕| 黑人操中国人逼视频| 精品人妻1区二区| 亚洲欧美一区二区三区黑人| 成人亚洲精品一区在线观看| 欧美成人免费av一区二区三区| 免费在线观看影片大全网站| 一区二区日韩欧美中文字幕| 99在线视频只有这里精品首页| 亚洲国产日韩欧美精品在线观看 | 国产精品 欧美亚洲| 国产一级毛片七仙女欲春2 | 欧美成人性av电影在线观看| 国产精品一区二区在线不卡| 91字幕亚洲| 亚洲黑人精品在线| 啦啦啦观看免费观看视频高清 | 亚洲成国产人片在线观看| 一边摸一边做爽爽视频免费| 亚洲av日韩精品久久久久久密| 欧美成人免费av一区二区三区| 精品国产亚洲在线| 久久中文字幕人妻熟女| 国产高清videossex| 美女午夜性视频免费| 久久久国产成人精品二区| 久久久久亚洲av毛片大全| 欧美日韩福利视频一区二区| 纯流量卡能插随身wifi吗| 一区在线观看完整版| 久久久久久免费高清国产稀缺| 午夜福利视频1000在线观看 | 久久这里只有精品19| 老汉色∧v一级毛片| 成人手机av| 久久久久久免费高清国产稀缺| a在线观看视频网站| 一二三四在线观看免费中文在| 亚洲aⅴ乱码一区二区在线播放 | 9191精品国产免费久久| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲欧美日韩另类电影网站| av超薄肉色丝袜交足视频| 亚洲国产精品合色在线| 日日干狠狠操夜夜爽| 搞女人的毛片| 国产伦人伦偷精品视频| 熟妇人妻久久中文字幕3abv| 一本久久中文字幕| 午夜两性在线视频| 国产成人影院久久av| 久热爱精品视频在线9| 人成视频在线观看免费观看| 欧美激情久久久久久爽电影 | 欧美成人午夜精品| 国产精品久久视频播放| 国产精品99久久99久久久不卡| 男女下面插进去视频免费观看| 亚洲国产精品合色在线| 亚洲av熟女| 欧美激情久久久久久爽电影 | 不卡一级毛片| 波多野结衣一区麻豆| 欧美日韩黄片免| 欧美成狂野欧美在线观看| 亚洲一区二区三区色噜噜| 久久午夜亚洲精品久久| 国产一级毛片七仙女欲春2 | 国产精品野战在线观看| 国产亚洲欧美在线一区二区| 国产亚洲精品综合一区在线观看 | 老汉色av国产亚洲站长工具| 中文字幕久久专区| 亚洲精品粉嫩美女一区| 亚洲色图综合在线观看| 天堂动漫精品| 欧美成人免费av一区二区三区| 国产午夜福利久久久久久| 久久久久久久精品吃奶| 欧美激情久久久久久爽电影 | 亚洲成人久久性| 中文字幕人妻熟女乱码| 天天躁夜夜躁狠狠躁躁| av片东京热男人的天堂| 男女午夜视频在线观看| 狂野欧美激情性xxxx| 日本欧美视频一区| 午夜成年电影在线免费观看| 久久婷婷人人爽人人干人人爱 | 国产一区二区激情短视频| 91精品国产国语对白视频| 露出奶头的视频| 在线国产一区二区在线| 免费看美女性在线毛片视频| 亚洲va日本ⅴa欧美va伊人久久| 两个人视频免费观看高清| 曰老女人黄片| 亚洲精华国产精华精| 黑人巨大精品欧美一区二区mp4| 男人的好看免费观看在线视频 | 侵犯人妻中文字幕一二三四区| 国产成人影院久久av| 亚洲精品在线美女| 欧美日韩一级在线毛片| 黑人巨大精品欧美一区二区蜜桃| 日本vs欧美在线观看视频| 91成人精品电影| 可以免费在线观看a视频的电影网站| 久久国产精品男人的天堂亚洲| 最新在线观看一区二区三区| 97人妻精品一区二区三区麻豆 | 免费看a级黄色片| 又大又爽又粗| 老汉色∧v一级毛片| 国产欧美日韩综合在线一区二区| 国产男靠女视频免费网站| 中文字幕色久视频| 精品国产一区二区三区四区第35| 在线播放国产精品三级| 国产精品二区激情视频| 在线观看舔阴道视频| 大型av网站在线播放| 中文亚洲av片在线观看爽| 麻豆国产av国片精品| 亚洲国产高清在线一区二区三 | 老司机福利观看| 中文亚洲av片在线观看爽| 国产一区二区三区在线臀色熟女| 国内久久婷婷六月综合欲色啪| 亚洲美女黄片视频| 国产欧美日韩综合在线一区二区| 亚洲熟女毛片儿| 久久精品国产亚洲av高清一级| 可以在线观看毛片的网站| 国产私拍福利视频在线观看| 91麻豆av在线| 久久人人精品亚洲av| 操出白浆在线播放| 久久精品91蜜桃| 满18在线观看网站| 一区二区三区国产精品乱码| 宅男免费午夜| 精品午夜福利视频在线观看一区| 国产高清视频在线播放一区| 国产成人影院久久av| 777久久人妻少妇嫩草av网站| 日韩欧美一区视频在线观看| 日本vs欧美在线观看视频| 悠悠久久av| 亚洲精品中文字幕一二三四区| 一本综合久久免费| 国产亚洲欧美在线一区二区| 91老司机精品| 男女午夜视频在线观看| 99久久国产精品久久久| 老汉色∧v一级毛片| 国产亚洲av嫩草精品影院| 国产精品免费视频内射| 18美女黄网站色大片免费观看| 日韩精品中文字幕看吧| 国产精品爽爽va在线观看网站 | 69精品国产乱码久久久| 国产亚洲精品久久久久久毛片| 一级片免费观看大全| 啦啦啦韩国在线观看视频| av欧美777| 首页视频小说图片口味搜索| 97碰自拍视频| 亚洲精品久久国产高清桃花| 免费无遮挡裸体视频| 国产精品自产拍在线观看55亚洲| 丁香六月欧美| 久久婷婷人人爽人人干人人爱 | 日本免费一区二区三区高清不卡 | 亚洲成人免费电影在线观看| 日韩大码丰满熟妇| 亚洲最大成人中文| 国产三级黄色录像| 大码成人一级视频| 欧美日本亚洲视频在线播放| 成年人黄色毛片网站| 一二三四在线观看免费中文在| 精品久久久久久久毛片微露脸| 夜夜躁狠狠躁天天躁| 日本精品一区二区三区蜜桃| 级片在线观看| 亚洲国产毛片av蜜桃av| 精品一区二区三区四区五区乱码| 精品电影一区二区在线| 后天国语完整版免费观看| 亚洲中文字幕日韩| 中亚洲国语对白在线视频| 在线观看免费视频日本深夜| 人人妻人人澡人人看| 中文字幕人成人乱码亚洲影| 欧美久久黑人一区二区| 亚洲人成77777在线视频| 久久精品亚洲熟妇少妇任你| 十八禁人妻一区二区| www.自偷自拍.com| 亚洲国产中文字幕在线视频| 天堂动漫精品| 最近最新中文字幕大全电影3 | 妹子高潮喷水视频| aaaaa片日本免费| 国产xxxxx性猛交| 国产日韩一区二区三区精品不卡| 久久久久久大精品| 国产乱人伦免费视频| 久久久国产成人免费| 人人妻人人澡欧美一区二区 | 欧美激情极品国产一区二区三区| 日韩国内少妇激情av| 久久久久国内视频| 91字幕亚洲| √禁漫天堂资源中文www| 国产欧美日韩精品亚洲av| 老司机深夜福利视频在线观看| 男人的好看免费观看在线视频 | 精品国产亚洲在线| 国产亚洲欧美98| 精品国产美女av久久久久小说| 波多野结衣巨乳人妻| 高清黄色对白视频在线免费看| 国产一区二区三区视频了| 女人被狂操c到高潮| 国产私拍福利视频在线观看| 91九色精品人成在线观看| 日韩大码丰满熟妇| 国产亚洲精品第一综合不卡| 成人18禁高潮啪啪吃奶动态图| 香蕉丝袜av| 不卡一级毛片| 美女扒开内裤让男人捅视频| 90打野战视频偷拍视频| 成人特级黄色片久久久久久久| e午夜精品久久久久久久| 日韩欧美三级三区| 无限看片的www在线观看| 波多野结衣一区麻豆| 美女午夜性视频免费| 91精品国产国语对白视频| 中国美女看黄片| 午夜影院日韩av| 国产在线精品亚洲第一网站| 久久人妻熟女aⅴ| 亚洲国产欧美日韩在线播放| √禁漫天堂资源中文www| 日韩视频一区二区在线观看| 亚洲三区欧美一区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲中文av在线| 精品少妇一区二区三区视频日本电影| 久久久久久亚洲精品国产蜜桃av| 天堂动漫精品| 午夜激情av网站| 桃色一区二区三区在线观看| 十八禁网站免费在线| 老司机在亚洲福利影院| 亚洲成a人片在线一区二区| 日本精品一区二区三区蜜桃| 美女国产高潮福利片在线看| 免费在线观看亚洲国产| 国产精品二区激情视频| 成年版毛片免费区| 亚洲欧美激情综合另类| 亚洲久久久国产精品| 97人妻天天添夜夜摸| 桃色一区二区三区在线观看| АⅤ资源中文在线天堂| 变态另类成人亚洲欧美熟女 | 欧美黄色片欧美黄色片| 亚洲一区中文字幕在线| 国产麻豆69| 国产精品免费一区二区三区在线| 欧美精品啪啪一区二区三区| 美女 人体艺术 gogo| 亚洲男人天堂网一区| 成年版毛片免费区| 国产亚洲欧美98| 久久久水蜜桃国产精品网| 国产精品98久久久久久宅男小说| av片东京热男人的天堂| 变态另类丝袜制服| 一边摸一边做爽爽视频免费| 人妻久久中文字幕网|