999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

稀磁半導(dǎo)體MnTe和CrTe的微觀結(jié)構(gòu)及電磁性能

2015-08-16 09:20:35楊光敏邢光宗
關(guān)鍵詞:結(jié)構(gòu)

徐 強(qiáng),楊光敏,邢光宗

(1.長(zhǎng)春工程學(xué)院 勘查與測(cè)繪工程學(xué)院,長(zhǎng)春 130021;2.長(zhǎng)春師范大學(xué) 物理學(xué)院,長(zhǎng)春 130032;3.吉林大學(xué) 材料科學(xué)與工程學(xué)院,長(zhǎng)春 130012)

?

稀磁半導(dǎo)體MnTe和CrTe的微觀結(jié)構(gòu)及電磁性能

徐 強(qiáng)1,楊光敏2,3,邢光宗3

(1.長(zhǎng)春工程學(xué)院 勘查與測(cè)繪工程學(xué)院,長(zhǎng)春 130021;2.長(zhǎng)春師范大學(xué) 物理學(xué)院,長(zhǎng)春 130032;3.吉林大學(xué) 材料科學(xué)與工程學(xué)院,長(zhǎng)春 130012)

基于密度泛函理論的VASP軟件包,研究CrTe和MnTe晶胞的穩(wěn)定性和電磁性能,并對(duì)比具有相同結(jié)構(gòu)的CrTe和MnTe在穩(wěn)定性和電磁性能方面的差異.結(jié)果表明,穩(wěn)態(tài)的CrTe為具有NiAs結(jié)構(gòu)的鐵磁體,其閃鋅礦結(jié)構(gòu)(ZB)、纖鋅礦結(jié)構(gòu)(W)和巖鹽結(jié)構(gòu)(RS)為亞穩(wěn)態(tài)結(jié)構(gòu),該亞穩(wěn)態(tài)結(jié)構(gòu)均為半金屬性質(zhì).

稀磁半導(dǎo)體;第一性原理;能帶結(jié)構(gòu);磁性

稀磁半導(dǎo)體(DMS)具有半導(dǎo)體和磁性的特點(diǎn),與普通磁性物質(zhì)相比,其磁性較弱[1-2].稀磁半導(dǎo)體主要包括La系稀土元素且具有NaCl結(jié)構(gòu)的硫化物和具有尖晶石結(jié)構(gòu)的三元化合物[3-6].Ⅱ-Ⅳ族的稀磁半導(dǎo)體具有良好的鐵磁性質(zhì),但其居里溫度低于室溫,Ⅲ-Ⅴ族的化合物易形成團(tuán)簇狀,不利于稀磁半導(dǎo)體生長(zhǎng).文獻(xiàn)[7]利用分子束外延技術(shù)在襯底上進(jìn)行非平衡生長(zhǎng),促進(jìn)了人們對(duì)稀磁半導(dǎo)體的研究.文獻(xiàn)[8-9]研究了純凈半金屬性稀磁半導(dǎo)體CrTe的磁性和電學(xué)性質(zhì)的第一性原理,發(fā)現(xiàn)其自旋向上具有金屬性質(zhì)而自旋向下是半導(dǎo)體性質(zhì),在Fermi能級(jí)附近具有帶隙.CrTe的基態(tài)具有NiAs型六角結(jié)構(gòu)(NA),其閃鋅礦結(jié)構(gòu)(ZB)、纖鋅礦結(jié)構(gòu)(W)和巖鹽結(jié)構(gòu)(RS)為亞穩(wěn)態(tài)結(jié)構(gòu)[10-11].本文構(gòu)建了CrTe(NiAs,RS,ZB,W)和MnTe(NiAs,ZB)的晶格結(jié)構(gòu),并分別對(duì)CrTe和MnTe的RS結(jié)構(gòu)、ZB結(jié)構(gòu)和NA結(jié)構(gòu)的穩(wěn)定性、電子結(jié)構(gòu)和磁性進(jìn)行研究,對(duì)比其結(jié)構(gòu)差異,進(jìn)而得到其電子結(jié)構(gòu)和磁性的差異.

1 理論計(jì)算

先分別建立具有閃鋅礦和纖鋅礦晶格結(jié)構(gòu)的CrTe和MnTe,再利用VASP軟件包分別計(jì)算其能帶結(jié)構(gòu)、態(tài)密度(DOS)、分波態(tài)密度(PDOS)和各原子的磁矩及總磁矩.閃鋅礦采用2×1×1晶胞(8個(gè)Cr原子和8個(gè)Te原子),纖鋅礦結(jié)構(gòu)采用2×2×1超晶胞(8個(gè)Cr原子和8個(gè)Te原子),巖鹽結(jié)構(gòu)采用2×1×1超晶胞(8個(gè)Cr原子和8個(gè)Te原子).在能帶計(jì)算過程中,閃鋅礦和纖鋅礦的能帶計(jì)算分別選用原胞和晶胞形式進(jìn)行計(jì)算,考慮到原子核和電子間的相互關(guān)系,使用基于密度泛函理論的離散Fourier變換(DFT)投影綴加平面波(PAW)贗勢(shì)處理,關(guān)聯(lián)效應(yīng)采用局域密度近似(LDA),截止能(ENCUT)為700 eV,用Monkhorst-Pack方法取Brillouin區(qū)的K點(diǎn).在優(yōu)化結(jié)構(gòu)過程中,加自旋(ISPIN)=2表示進(jìn)行磁性優(yōu)化,根據(jù)Hund規(guī)則取Cr和Te的初始磁矩分別為4,0 μB.K點(diǎn)選用線型模式研究能帶結(jié)構(gòu),高對(duì)稱點(diǎn)的積分路徑為G-K-H-A-G-M-L-A,高對(duì)稱點(diǎn)的坐標(biāo)分別為G(0,0,0),K(-0.33,0.66,0),H(-0.33,0.66,0.5),A(0,0,0.5),M(0,0.5,0),L(0,0.5,0.5).態(tài)密度計(jì)算選用的K點(diǎn)為15×15×7.

2 結(jié)果與討論

2.1CrTe和MnTe結(jié)構(gòu)的穩(wěn)定性

CrTe可能穩(wěn)定存在4種結(jié)構(gòu):閃鋅礦結(jié)構(gòu)、NiAs結(jié)構(gòu)、巖鹽結(jié)構(gòu)和纖鋅礦結(jié)構(gòu),其模型如圖1所示.由圖1(A)可見,晶胞中有8個(gè)Te原子位于Cr原子組成的四面體體心,晶格常數(shù)a=b=c=0.627 46 nm,α=β=90°,γ=120°;由圖1(B)可見,晶胞由4個(gè)原子構(gòu)成,晶格常數(shù)a=b=0.442 nm,c=0.742 nm,α=β=90°,γ=120°;由圖1(C)可見,晶格常數(shù)a=b=0.427 856 nm,c=0.645 38 nm,α=β=90°,γ=120°;圖1(D)的巖鹽結(jié)構(gòu)即NaCl結(jié)構(gòu),晶胞由4個(gè)Cr原子和4個(gè)Te原子組成,晶格常數(shù)a=b=c=0.585 58 nm,α=β=γ=90°.

圖1 CrTe的閃鋅礦結(jié)構(gòu)(A)、纖鋅礦結(jié)構(gòu)(B)、NiAs結(jié)構(gòu)(C)和巖鹽結(jié)構(gòu)(D)Fig.1 Four structures of CrTe with zinc-blende (A),wurtzite (B),NiAs (C),rock salt (D)

通過計(jì)算總能量可初步判斷其結(jié)構(gòu)的穩(wěn)定性.在常規(guī)算法中令NSW=0,即不做離子步運(yùn)算.在LDA+U方法中,U為經(jīng)驗(yàn)參數(shù),通過與電子能譜對(duì)比的方法可確定U值.通過靜態(tài)計(jì)算所得體積-能量曲線如圖2所示,其中圖2(A)和(B)分別為L(zhǎng)DA+U算法和常規(guī)算法所得CrTe的體積-能量收斂曲線.由圖2可見,其存在能量最低點(diǎn).用原胞的平均能量進(jìn)行比較,數(shù)據(jù)處理后的結(jié)果列于表1.由表1可見,交換常數(shù)對(duì)結(jié)構(gòu)總能量的影響較大,因此在計(jì)算中需考慮交換常數(shù),總能量為EZB

圖2 LDA+U算法(A)和常規(guī)算法(B)的CrTe體積-能量收斂曲線Fig.2 Volume-energy convergence curves of CrTe by LDA+U algorithm (A)and normal algorithm (B)

表1 不同算法得到不同結(jié)構(gòu)的最低能量(eV)Table 1 The lowest energies (eV)of different structures by different algorithms

通過靜態(tài)計(jì)算所得MnTe體積-能量曲線如圖3所示.由圖3可見,MnTe的ZB結(jié)構(gòu)低于NiAs結(jié)構(gòu)總能量.CrTe和MnTe的閃鋅礦結(jié)構(gòu)、巖鹽結(jié)構(gòu)和NiAs結(jié)構(gòu)能量較低.CrTe和MnTe形成能的表達(dá)式分別為E形成=ECrTe-ECr-ETe和E形成=EMnTe-EMn-ETe.若形成能小于零,則材料處于基態(tài);若形成能大于零,則材料處于亞穩(wěn)態(tài)或未形成該材料.形成能的計(jì)算結(jié)果列于表2.

Fig.3 MnTe的LDA+U算法(A)和常規(guī)算法(B)的體積-能量曲線Fig.3 Volume-energy convergence curves of MnTe by LDA+U algorithm (A)and normal algorithm (B)

表2 CrTe和MnTe各種結(jié)構(gòu)的形成能(eV)Table 2 Formation energies (eV)of different structures for CrTe and MnTe

2.2CrTe和MnTe的電子結(jié)構(gòu)

NiAs結(jié)構(gòu)CrTe的能帶結(jié)構(gòu)如圖4所示,其中圖4(A)和(B)分別表示自旋向下和自旋向上電子沿高對(duì)稱點(diǎn)的能帶結(jié)構(gòu).

圖4 NiAs結(jié)構(gòu)CrTe的能帶結(jié)構(gòu)Fig.4 Spin band structures of NiAs structure for CrTe

由圖4可見:NiAs結(jié)構(gòu)的CrTe具有半金屬性質(zhì),其自旋向下的電子呈半導(dǎo)體性質(zhì),由于價(jià)帶的最高點(diǎn)在高對(duì)稱點(diǎn)G處Fermi能級(jí)附近,導(dǎo)帶的最低點(diǎn)在高對(duì)稱點(diǎn)M處,因此其為間接帶隙半導(dǎo)體;自旋向上的電子呈金屬性質(zhì),其Fermi能級(jí)穿過能帶,不存在帶隙;自旋向下電子的半金屬帶隙為1.44 eV.NiAs結(jié)構(gòu)CrTe的態(tài)密度如圖5所示,其中圖5(A)~(C)分別對(duì)應(yīng)總態(tài)密度、Cr的分波態(tài)密度和Te的分波態(tài)密度.由圖5(A)可見,自旋向上和自旋向下電子的態(tài)密度不對(duì)稱,即該物質(zhì)具有磁性;由圖5(B)和(C)可見,Fermi能級(jí)附近的態(tài)密度主要由Te原子的4p態(tài)組成,價(jià)帶主要由Te原子的4p電子態(tài)和3d電子態(tài)組成.

圖5 NiAs結(jié)構(gòu)CrTe的態(tài)密度 Fig.5 Densities of states of NiAs structure for CrTe

RS結(jié)構(gòu)CrTe的能帶結(jié)構(gòu)如圖6所示,其中圖6(A)和(B)分別表示自旋向下和自旋向上電子沿高對(duì)稱點(diǎn)的能帶結(jié)構(gòu).由圖6可見:RS結(jié)構(gòu)的CrTe具有半金屬性質(zhì),其自旋向下的電子呈半導(dǎo)體性質(zhì),由于價(jià)帶頂和導(dǎo)帶底均在高對(duì)稱點(diǎn)G處,因此其為直接帶隙半導(dǎo)體;自旋向上電子組成的能帶呈金屬性質(zhì),其帶隙為1.96 eV.

圖6 RS結(jié)構(gòu)的CrTe的能帶結(jié)構(gòu)Fig.6 Spin band structures of RS structure for CrTe

RS結(jié)構(gòu)CrTe的態(tài)密度如圖7所示.由圖7可見:RS結(jié)構(gòu)的CrTe具有磁性;Fermi能級(jí)附近的態(tài)密度主要由Te原子的5p電子態(tài)組成,價(jià)帶主要由Te原子的5p電子態(tài)和4d電子態(tài)組成.ZB結(jié)構(gòu)CrTe的能帶結(jié)構(gòu)如圖8所示,其中圖8(A)和(B)分別表示自旋向下和自旋向上電子沿高對(duì)稱點(diǎn)的能帶結(jié)構(gòu).由圖8可見:ZB結(jié)構(gòu)的CrTe具有半金屬性質(zhì),由于自旋向下電子組成的能帶價(jià)帶頂和導(dǎo)帶底均在高對(duì)稱點(diǎn)G處,因此其為直接帶隙半導(dǎo)體;自旋向上電子組成的能帶穿過Fermi能級(jí),導(dǎo)帶和價(jià)帶重合,使其具有金屬性質(zhì);自旋向下電子組成的半金屬帶隙為3.34 eV.

ZB結(jié)構(gòu)CrTe的態(tài)密度如圖9所示.由圖9可見:ZB結(jié)構(gòu)的CrTe具有磁性;Fermi能級(jí)附近的態(tài)密度主要由Cr原子的3d態(tài)和4p態(tài)組成,價(jià)帶主要由Cr原子的4s和4p態(tài)組成.

圖7 RS結(jié)構(gòu)CrTe的態(tài)密度Fig.7 Densities of states of RS structure for CrTe

圖8 ZB結(jié)構(gòu)的CrTe能帶結(jié)構(gòu)Fig.8 Spin band structures of ZB structure for CrTe

圖9 ZB結(jié)構(gòu)CrTe的態(tài)密度Fig.9 Densities of states of ZB structure for CrTe

W結(jié)構(gòu)CrTe的能帶結(jié)構(gòu)如圖10所示,其中圖10(A)和(B)分別表示自旋向下和自旋向上電子沿高對(duì)稱點(diǎn)的能帶結(jié)構(gòu).由圖10可見:W結(jié)構(gòu)的CrTe具有半金屬性質(zhì),由于自旋向下電子組成的能帶價(jià)帶頂和導(dǎo)帶底均在高對(duì)稱點(diǎn)G處,因此其為直接帶隙半導(dǎo)體;自旋向上電子組成的能帶呈金屬性質(zhì),自旋向下電子的半金屬帶隙為3.65 eV.W結(jié)構(gòu)CrTe的態(tài)密度如圖11所示.由圖11可見:W結(jié)構(gòu)的CrTe具有磁性;Fermi能級(jí)附近的態(tài)密度主要由Te原子的5p和4d電子態(tài)組成.

圖10 W結(jié)構(gòu)CrTe的能帶結(jié)構(gòu)Fig.10 Spin band structures of W structure for CrTe

圖11 W結(jié)構(gòu)CrTe的態(tài)密度Fig.11 Densities of states of W structure for CrTe

NiAs結(jié)構(gòu)MnTe的能帶結(jié)構(gòu)如圖12所示,其中圖12(A)和(B)分別表示自旋向下和自旋向上電子沿高對(duì)稱點(diǎn)的能帶結(jié)構(gòu).由圖12可見,NiAs結(jié)構(gòu)的MnTe具有金屬性質(zhì),這是由于自旋向上和自旋向下電子均未產(chǎn)生帶隙,Fermi能級(jí)穿過能帶且價(jià)帶和導(dǎo)帶重合所致.NiAs結(jié)構(gòu)MnTe的態(tài)密度如圖13所示.由圖13可見,自旋向上和自旋向下的電子態(tài)密度不對(duì)稱,即該物質(zhì)具有磁性.

ZB結(jié)構(gòu)MnTe的能帶結(jié)構(gòu)如圖14所示,其中圖14(A)和(B)分別表示自旋向下和自旋向上電子沿高對(duì)稱點(diǎn)的能帶結(jié)構(gòu).由圖14可見:ZB結(jié)構(gòu)的MnTe具有半導(dǎo)體性質(zhì),自旋向上和自旋向下電子組成的能帶結(jié)構(gòu)均產(chǎn)生直接帶隙,且價(jià)帶頂和導(dǎo)帶底均在高對(duì)稱點(diǎn)G處;自旋向下電子的帶隙為2.73 eV,自旋向上電子的帶隙為1.08 eV.ZB結(jié)構(gòu)MnTe的態(tài)密度如圖15所示.由圖15可見:自旋向下和自旋向上電子在Fermi能級(jí)附近的態(tài)密度均為零,表明其具有半導(dǎo)體性質(zhì);在Fermi能級(jí)附近的價(jià)帶主要有Mn原子的3p電子和3d電子組成,導(dǎo)帶主要由Te原子的5s電子和4d電子組成.

2.3CrTe和MnTe結(jié)構(gòu)的磁性

進(jìn)行靜態(tài)計(jì)算時(shí),設(shè)ISPIN=2,并根據(jù)Hund規(guī)則設(shè)Mn,Cr,Te的原子初始磁矩分別為4,4,0 μB,通過計(jì)算所得CrTe和MnTe各種結(jié)構(gòu)的平均總磁矩列于表3.由表3可見,CrTe和MnTe各種結(jié)構(gòu)的磁矩接近Cr原子和Mn原子的磁矩,Te原子對(duì)磁矩貢獻(xiàn)較小.

圖12 NiAs結(jié)構(gòu)MnTe的能帶結(jié)構(gòu)Fig.12 Spin band structures of NiAs structure for MnTe

圖13 NiAs結(jié)構(gòu)MnTe的態(tài)密度Fig.13 Densities of states of NiAs structure for MnTe

圖14 ZB結(jié)構(gòu)MnTe能帶結(jié)構(gòu)Fig.14 Spin band structures of ZB structure for MnTe

圖15 ZB結(jié)構(gòu)MnTe的態(tài)密度Fig.15 Densities of states of ZB structure for MnTe

表3 CrTe和MnTe不同結(jié)構(gòu)的總磁矩(μB)Table 3 Total magnetic moments (μB)of different structures for CrTe and MnTe

綜上,本文可得如下結(jié)論:1)CrTe可能存在的4種晶格結(jié)構(gòu)為性能優(yōu)良的稀磁半導(dǎo)體,4種結(jié)構(gòu)的穩(wěn)定性為NiAs>RS>ZB>W,MnTe的穩(wěn)定性為NiAs>ZB;2)CrTe和MnTe不同結(jié)構(gòu)對(duì)應(yīng)的總磁矩均與Cr原子和Mn原子總磁矩接近,即兩種結(jié)構(gòu)的磁矩均主要由金屬原子貢獻(xiàn),非金屬原子Te對(duì)磁矩貢獻(xiàn)較小;4種不同結(jié)構(gòu)的CrTe均呈良好的半金屬性質(zhì),MnTe的NiAs結(jié)構(gòu)呈金屬性質(zhì),ZB結(jié)構(gòu)呈半導(dǎo)體性質(zhì).

[1] Ohno H.Making Nonmagnetic Semiconductors Ferromagnetic [J].Science,1998,281:951-956.

[2] Zutic I,Fabian J,Das S S.Spintronics:Fundamentals and Applications [J].Rev Mod Phys,2004,76(2):323-410.

[3] Dietl T.A Ten-Year Perspective on Dilute Magnetic Semiconductors and Oxides [J].Nature Mater,2010,9(12):965-974.

[4] Liu C,Yun F,Morkoc H.Ferromagnetism of ZnO and GaN:A Review [J].Journal of Materials Science,2005,16(9):555-597.

[5] Murakami M,Matsumoto Y,Hasegawa T,et al.Cobalt Valence States and Origins of Ferromagnetism in Co Doped TiO2Rutile Thin Films [J].J Appl Phys,2004,95(10):5330-5333.

[6] Kasuya T,Yanase A.Anomalous Transport Phenomena in Eu-Chalcogenide Alloys [J].Rev Mod Phys,1968,40(4):684-696.

[7] Munekata H,Ohno H,Molnar S,von,et al.Diluted Magnetic Ⅲ-Ⅴ Semiconductors [J].Phys Rev Lett,1989,63(17):1849-1852.

[8] Liu Y,Bose S K,KudrnovskJ.First-Principles Theoretical Studies of Half-Metallic Ferromagnetism in CrTe [J].Phys Rev B,2010,82(9):094435.

[9] Dijkstra J,Weitering H H,Bruggen C F,van,et al.Band-Structure Calculations,and Magnetic and Transport-Properties of Ferromagnetic Chromium Tellurides (CrTe,Cr3Te4,Cr2Te3)[J].J Phys:Condens Matter,1989,1(46):9141-9161.

[10] ZHANG Ming,Brück E,Boeret F R,de,et al.Half-Metallic Ferromagnetism in Hypothetical Wurtzite Structure Chromium Chalcogenides [J].J Mater Res,2004,19(9):2738-2741.

[11] ZHAO Yonghong,XIE Wenhui,ZHU Lifang,et al.Half-Metallic Ferromagnets Based on the Rock-Salt Ⅳ-Ⅵ Semiconductor GeTe [J].J Phys:Condens Matter,2006,18(45):10259-10268.

(責(zé)任編輯:王 健)

MicrostructureandElectromagnetismPropertiesofMnTeandCrTeDilutedMagneticSemiconductor

XU Qiang1,YANG Guangmin2,3,XING Guangzong3

(1.SchoolofProspectingandSurveying,ChangchunInstituteofTechnology,Changchun130021,China;2.CollegeofPhysics,ChangchunNormalUniversity,Changchun130032,China;3.CollegeofMaterialsScienceandEngineering,JilinUniversity,Changchun130012,China)

Diluted magnetic semiconductors (DMS)which are used as spin-electron devices have been researched popularly because of their excellent electronic and magnetic properties.We used VASP code based on the density functional theory to calculate the stability,electronic and magnetic properties of the unit cell of CrTe,MnTe and compared the stability,electric and magnetic properties of CrTe with those of MnTe with the same structure as that of CrTe.We found that the stable state of CrTe exibits the NiAs-type hexagonal structure.The metastable state of it displays three structures called zinc-blende structure,wurtzite structure and rock salt structure,the three structures are all half metallic properties.

diluted magnetic semiconductors;the first principle;band structure;magnesium

10.13413/j.cnki.jdxblxb.2015.03.38

2014-12-21.

徐 強(qiáng)(1981—),男,漢族,碩士,講師,從事凝聚態(tài)物理的研究,E-mail:xuqiang19810919@163.com.通信作者:楊光敏(1981—),女,漢族,博士,副教授,從事凝聚態(tài)物理的研究,E-mail:249138087@qq.com.

吉林省教育廳“十二五”科學(xué)技術(shù)研究項(xiàng)目(批準(zhǔn)號(hào):2014-250)和吉林省自然科學(xué)基金(批準(zhǔn)號(hào):2014-0101061JC).

O484

:A

:1671-5489(2015)03-0553-08

猜你喜歡
結(jié)構(gòu)
DNA結(jié)構(gòu)的發(fā)現(xiàn)
《形而上學(xué)》△卷的結(jié)構(gòu)和位置
論結(jié)構(gòu)
中華詩詞(2019年7期)2019-11-25 01:43:04
新型平衡塊結(jié)構(gòu)的應(yīng)用
模具制造(2019年3期)2019-06-06 02:10:54
循環(huán)結(jié)構(gòu)謹(jǐn)防“死循環(huán)”
論《日出》的結(jié)構(gòu)
縱向結(jié)構(gòu)
縱向結(jié)構(gòu)
我國(guó)社會(huì)結(jié)構(gòu)的重建
人間(2015年21期)2015-03-11 15:23:21
創(chuàng)新治理結(jié)構(gòu)促進(jìn)中小企業(yè)持續(xù)成長(zhǎng)
主站蜘蛛池模板: 亚洲欧美日韩高清综合678| 亚洲欧美日韩成人高清在线一区| 婷五月综合| 一级毛片免费高清视频| 国产精品一区在线麻豆| 亚洲色图另类| 激情无码视频在线看| 一级爆乳无码av| 国产人成在线观看| 高清无码一本到东京热| 婷婷伊人五月| 91欧美在线| 人妻无码中文字幕第一区| 国产一区二区免费播放| 日韩精品无码免费专网站| 91久久大香线蕉| 亚洲欧美成人| 在线a视频免费观看| 在线免费看片a| 免费国产黄线在线观看| 国产jizz| 美女潮喷出白浆在线观看视频| 久久天天躁狠狠躁夜夜躁| 九九热精品视频在线| 一级成人欧美一区在线观看| 日韩午夜片| 99久久无色码中文字幕| 成年看免费观看视频拍拍| 国产麻豆另类AV| 国产精品美女自慰喷水| 欧美福利在线观看| 精品一区二区久久久久网站| 久久黄色免费电影| 亚洲精品卡2卡3卡4卡5卡区| 中文字幕在线永久在线视频2020| 国产人碰人摸人爱免费视频| 91青青草视频| 国产波多野结衣中文在线播放| 国产H片无码不卡在线视频| 久久女人网| 日韩小视频网站hq| 国产成人精品一区二区三在线观看| 亚洲男人天堂网址| 国产成人高清在线精品| 日本成人福利视频| 欧美啪啪网| 国产精品自在自线免费观看| 自慰网址在线观看| 中文字幕人成乱码熟女免费| 天天摸天天操免费播放小视频| 亚洲成人黄色在线观看| 51国产偷自视频区视频手机观看| 国产91麻豆免费观看| 第一区免费在线观看| 男人的天堂久久精品激情| 国产成人免费观看在线视频| 区国产精品搜索视频| 日韩不卡免费视频| 亚洲中文字幕日产无码2021| 怡春院欧美一区二区三区免费| 91青青视频| 中文字幕1区2区| 国产精品毛片在线直播完整版| 国产自在线拍| 国产视频一二三区| 欧美特黄一级大黄录像| 婷婷五月在线视频| 日本午夜视频在线观看| 欧美a在线视频| 日本一区二区不卡视频| 爽爽影院十八禁在线观看| 亚洲美女久久| 久久99热这里只有精品免费看| 久久久久久午夜精品| 国产精品播放| 久久黄色免费电影| 欧美日本在线观看| 欧美福利在线| 欧美在线黄| 亚洲国产综合精品一区| 成人毛片在线播放| 啪啪国产视频|