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

    鐵磁性過(guò)渡離子摻雜對(duì)纖鋅礦相硫化鋅電子結(jié)構(gòu)的影響

    2012-09-15 11:45:36董會(huì)寧張博陳沛達(dá)張翠玲劉
    關(guān)鍵詞:金屬性鐵磁性價(jià)帶

    董會(huì)寧張 博陳沛達(dá)張翠玲劉 俊

    (1重慶郵電大學(xué)數(shù)理學(xué)院,重慶 400065)

    (2重慶郵電大學(xué)通信學(xué)院,重慶 400065)

    (3重慶工商大學(xué)計(jì)算機(jī)科學(xué)與信息工程學(xué)院,重慶 400067)

    鐵磁性過(guò)渡離子摻雜對(duì)纖鋅礦相硫化鋅電子結(jié)構(gòu)的影響

    董會(huì)寧*,1張 博2陳沛達(dá)1張翠玲3劉 ?。?1

    (1重慶郵電大學(xué)數(shù)理學(xué)院,重慶 400065)

    (2重慶郵電大學(xué)通信學(xué)院,重慶 400065)

    (3重慶工商大學(xué)計(jì)算機(jī)科學(xué)與信息工程學(xué)院,重慶 400067)

    通過(guò)第一性原理計(jì)算,優(yōu)化了鐵磁性過(guò)渡離子摻雜的纖鋅礦相硫化鋅Fm0.125Zn0.875S(Fm=Fe、Co、Ni)的幾何結(jié)構(gòu),計(jì)算了其電子結(jié)構(gòu),分析了其半金屬性及其微觀機(jī)制。結(jié)果表明:對(duì)不同的鐵磁性雜質(zhì)離子,F(xiàn)m0.125Zn0.875S在費(fèi)米面處的自旋極化率均為-100%,具有半金屬性,是潛在的優(yōu)質(zhì)自旋注入材料。Fm0.125Zn0.875S具有較寬的自旋帶隙,從而具有較高的居里溫度和廣泛的應(yīng)用前景。Fe0.125Zn0.875S、Co0.125Zn0.875S 和 Ni0.125Zn0.875S的2×2×1 超胞的磁矩分別為 3.96μB、2.90μB和 2.00μB,主要來(lái)自于鐵磁性過(guò)渡離子Fe、Co 和 Ni離子。這 3 種離子的電子結(jié)構(gòu)分別為。

    半金屬鐵磁體;電子結(jié)構(gòu);超胞磁矩

    1983年de Groot等從理論上發(fā)現(xiàn)了三元合金NiMnSb和PtMnSb具有與普通鐵磁性材料不同的輸運(yùn)性質(zhì)[1],它們的一種自旋取向的電子的能帶結(jié)構(gòu)呈現(xiàn)金屬性,即費(fèi)米面處于導(dǎo)帶中,具有金屬的行為;而另一種自旋取向的電子的能帶則呈現(xiàn)絕緣體性質(zhì)或半導(dǎo)體性質(zhì),即費(fèi)米面處于帶隙中。這種材料因具有單自旋輸運(yùn)而被稱為半金屬鐵磁體(HMFs)。由于半金屬材料具有±100%的自旋極化率,且電阻率較普通鐵磁金屬更接近于半導(dǎo)體材料,因此半金屬鐵磁體是一種理想的自旋半導(dǎo)體的自旋注入源(自旋電極材料)[2]。高自旋極化率使得半金屬鐵磁體在自旋閥、隧道結(jié)等磁電子器件中具有廣泛的應(yīng)用前景,它的使用可以明顯改善磁電阻器件的質(zhì)量和物理性能[3-8],故受到眾多研究者的密切關(guān)注。已有多種半金屬材料被發(fā)現(xiàn),如Heusler合金Co2MnSi[9-10]和 Co2FeSi[11];金屬氧化物CrO2[12]和Fe3O4[13];鈣鈦礦化合物L(fēng)a0.7Sr0.3MnO3[14]和 Sr2FeMo O6[15];閃鋅礦化合物CrAs和CrSb[16-17]等。越來(lái)越多理論預(yù)測(cè)的半金屬鐵磁體已被實(shí)驗(yàn)證實(shí)。另外,人們?cè)谶^(guò)渡金屬摻雜的半導(dǎo)體中也發(fā)現(xiàn)了半金屬性,如Mn 和 Cr摻雜的 ZnTe[18-19]、Cr摻雜的 ZnSe[20]及 Zn(Tm)O2(Tm=Cr、Mn、Fe、Co 和 Ni)[21]。ZnS 是一種常見(jiàn)半導(dǎo)體材料,其結(jié)構(gòu)有纖鋅礦和閃鋅礦結(jié)構(gòu)兩種。近年來(lái)人們對(duì)摻雜的ZnS材料及其納米顆粒進(jìn)行了大量的實(shí)驗(yàn)研究。例如,ZnS∶Ag[22],ZnS∶Cu[23],ZnS∶Mn[24],ZnS∶Co[25-26]等納米顆粒等在實(shí)驗(yàn)上也被成功合成,F(xiàn)e0.44Zn0.56S[27],Co0.25Zn0.75S[28],Ni0.03Zn0.97S[29]等 塊 體 材料也已被成功合成。在理論方面,人們已通過(guò)第一原理計(jì)算對(duì)閃鋅礦結(jié)構(gòu)的Mn,C和Fe等摻雜的ZnS進(jìn)行了研究[24,30],但對(duì)纖鋅礦結(jié)構(gòu)的過(guò)渡離子摻雜的ZnS的研究較少,尤其是對(duì)鐵磁性元素?fù)诫s的ZnS的磁性能和半金屬性的研究更少見(jiàn)報(bào)道?;诖?,本文擬通過(guò)第一原理計(jì)算來(lái)詳細(xì)研究鐵磁性過(guò)渡元素?fù)诫s的 ZnS 化合物,即 Fm0.125Zn0.875S(Fm=Fe、Co和Ni)的電子結(jié)構(gòu)、磁學(xué)性質(zhì),分析摻雜對(duì)體系電學(xué)和磁學(xué)性質(zhì)的影響,并根據(jù)配位場(chǎng)理論,探討相關(guān)性質(zhì)的微觀機(jī)理。

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

    本文選用纖鋅礦結(jié)構(gòu)ZnS為基體,其晶格常數(shù)的實(shí)驗(yàn)值為 a=b=0.382 nm、c=0.626 nm[31],空間群為P63mc,其晶胞由S的六角密堆積和Zn的六角密堆積反向套構(gòu)而成。為與實(shí)驗(yàn)摻雜含量(10%~30%)一致,我們選取包含8個(gè)Zn和8個(gè)S原子的2×2×1的超胞,分別用鐵磁性過(guò)渡元素Fe、Co和Ni取代超胞中的一個(gè)Zn原子實(shí)現(xiàn)摻雜,故雜質(zhì)含量為12.5%。摻雜后的超胞見(jiàn)圖1。

    計(jì)算軟件采用Materials studio 5.5的Castep模塊[32-33]。該模塊基于密度泛函理論,由劍橋大學(xué)凝聚態(tài)理論研究組開(kāi)發(fā)而成。它采用周期性邊界條件,利用平面波贗勢(shì)的方法,將離子勢(shì)用贗勢(shì)替代。電子波函數(shù)通過(guò)平面波基組展開(kāi),電子-電子交換相互作用和相關(guān)勢(shì)由廣義梯度近似進(jìn)行校正。孤立的Fe、Co、Ni、S和Zn原子的電子組態(tài)分 別為 Fe:3s23p63d64s2,Co:3s23p63d74s2,Ni:3s23p63d84s2,Zn:3d104s2和S:3s23p4。計(jì)算精度采用Fine,相應(yīng)的截?cái)鄤?dòng)能為340 eV,k 網(wǎng)格點(diǎn)為 9×4×5。

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

    2.1 優(yōu)化后的幾何結(jié)構(gòu)參數(shù)

    首先對(duì)ZnS的2×2×1超胞的幾何結(jié)構(gòu)進(jìn)行優(yōu)化并與實(shí)驗(yàn)值進(jìn)行比較以檢驗(yàn)本文計(jì)算方法的可靠性。本文計(jì)算所得的晶格常數(shù)為a=b=0.7858 nm、c=0.6469 nm,c/a=0.8232,折合為 1×1×1 超胞后,c/a=1.647 6,與文獻(xiàn)[34]中的實(shí)驗(yàn)值 1.638 7 偏差不到1.00%,表明了本文計(jì)算方法的可靠性。

    優(yōu)化后 Fm0.125Zn0.875S 的幾何結(jié)構(gòu)參數(shù)見(jiàn)表 1。由表1可以看出,由于Fe、Co、Ni核外電子數(shù)比Zn原子少,摻入鐵磁性過(guò)渡金屬后材料的晶格常數(shù)變小。但是,由于Fe、Co、Ni與Zn原子半徑相差不大,因此摻雜后的Fm0.125Zn0.875S晶格常數(shù)變化不大,c/a值變化不超過(guò)1.00%。由于晶格常數(shù)會(huì)影響半導(dǎo)體的禁帶寬度,因此,摻雜后Fm0.125Zn0.875S的帶隙會(huì)有所改變。

    表1 優(yōu)化后Fm0.125Zn0.875S(Fm=Fe、Co和Ni)的2×2×1超胞的幾何結(jié)構(gòu)參數(shù)Table 1 Geometric parameters of the optimized 2×2×1 supercells of Fm0.125Zn0.875S,where Fm=Fe,Co and Ni

    2.2 自旋極化的總態(tài)密度

    圖2給出了Fm0.125Zn0.875S的自旋極化的總態(tài)密度(TDOS),圖中實(shí)線和虛線分別表示自旋向上和向下的總態(tài)密度,圖中的能量相對(duì)于費(fèi)米能而定義的,即令費(fèi)米能為零。從圖2,TDOS處于-70到6 eV之間,被分割成7個(gè)部分,為描述方便,定義為7個(gè)子能帶。在-70到-10 eV之間的3個(gè)子能帶中,前2個(gè)子能帶關(guān)于態(tài)密度為零的橫軸呈現(xiàn)錯(cuò)位對(duì)稱分布,第3個(gè)子能帶則呈現(xiàn)軸對(duì)稱分布,因此這些子能帶中的自旋向上和自旋向下電子數(shù)相同,即電子的總自旋為零。同時(shí),它們距離費(fèi)米面較遠(yuǎn),因此對(duì)材料的磁學(xué)和電學(xué)性能均沒(méi)有影響,因此未在圖上畫出。在能量-10到6 eV范圍內(nèi),F(xiàn)m0.125Zn0.875S分布著4個(gè)子能帶,其中第1個(gè)和第2個(gè)子能帶以橫軸(DOS=0)為軸呈現(xiàn)軸對(duì)稱分布。對(duì)不同的雜質(zhì)元素Fe、Co和Ni,F(xiàn)m0.125Zn0.875S在此能量范圍內(nèi)有相似的態(tài)密度分布,主要原因在于Fe、Co和Ni原子此能量范圍內(nèi)存在相似的電子軌道結(jié)構(gòu)。第3個(gè)子能帶中自旋向下的子帶橫跨費(fèi)米面而自旋向上的子帶則體現(xiàn)出半導(dǎo)體性,即只有自旋向下的電子可以在材料中傳輸,因此具有半金屬性。該子能帶對(duì)Fm0.125Zn0.875S的磁學(xué)和電學(xué)性能有重要影響。此外,與Fe摻雜相比,Co和Ni摻雜后費(fèi)米面整體相對(duì)向低能方向移動(dòng)。這可能是因?yàn)殡S著摻雜元素的電子數(shù)的增加,超胞內(nèi)電子的庫(kù)侖排斥相互作用及電子交換作用變強(qiáng),從而使得體系的總能量有增加的趨勢(shì),即能帶整體向高能方向移動(dòng)。同時(shí),F(xiàn)e0.125Zn0.875S,Co0.125Zn0.875S和Ni0.125Zn0.875S的費(fèi)米能分別為3.528,3.230和3.332 eV,后兩者明顯小于前者。上述2個(gè)因素使得后兩者的費(fèi)米面相對(duì)它們的價(jià)帶整體向低能方向移動(dòng)。因此,F(xiàn)e0.125Zn0.875S的費(fèi)米面與自旋向上價(jià)帶頂?shù)木嚯x明顯大于Co0.125Zn0.875S和Ni0.125Zn0.875S。這表明Fe0.125Zn0.875S的半金屬性最穩(wěn)定。第4個(gè)子帶的能量大于費(fèi)米能,且呈現(xiàn)近似軸對(duì)稱分布,對(duì)材料的性能影響也較小。由于Co0.125Zn0.875S和Ni0.125Zn0.875S的自旋向上子價(jià)帶的帶頂較接近于費(fèi)米能,當(dāng)系統(tǒng)受到外界的微小擾動(dòng)如熱擾動(dòng)或電磁擾動(dòng)等,極可能使該價(jià)帶頂?shù)碾娮幼孕D(zhuǎn)到自旋向下的導(dǎo)帶中,因此在該子價(jià)帶中可能存在著自旋向上的空穴。因此Co0.125Zn0.875S和Ni0.125Zn0.875S的導(dǎo)電性主要來(lái)自自旋向下的電子,但可能部分來(lái)源于自旋向上的空穴。相反,F(xiàn)e0.125Zn0.875S的自旋向上價(jià)帶頂距離費(fèi)米面較遠(yuǎn),因此微小的外界擾動(dòng)很難使得該價(jià)帶頂?shù)碾娮臃D(zhuǎn)到自旋向下的子帶中。因此,F(xiàn)e0.125Zn0.875S中一般不存在自旋向上的空穴,其導(dǎo)電性主要來(lái)自自旋向下的傳導(dǎo)電子。因此Fe0.125Zn0.875S的半金屬性明顯比Co0.125Zn0.875S和Ni0.125Zn0.875S穩(wěn)定。

    從圖2的放大圖,摻雜后Fm0.125Zn0.875S的自旋帶隙,即自旋向上子價(jià)帶頂?shù)搅硗?,由圖2的放大圖可看出Fe0.125Zn0.875S,Co0.125Zn0.875S和Ni0.125Zn0.875S的自旋向上子價(jià)帶頂?shù)阶訉?dǎo)帶底的距離分別為2.12、1.98和1.88 eV。因此,摻Fe元素后的自旋帶隙明顯大于摻Co和Ni元素,表明Fe0.125Zn0.875S的半金屬性明顯比Co0.125Zn0.875S和Ni0.125Zn0.875S穩(wěn)定。由于Castep軟件對(duì)激發(fā)態(tài)電子的關(guān)聯(lián)效應(yīng)的低估,因此理論上得出的帯隙明顯小于實(shí)驗(yàn)值[35],這表明摻雜后材料帶隙的實(shí)驗(yàn)值極可能大于計(jì)算值。因此,可預(yù)測(cè)這些半金屬材料有較寬的自旋帶隙,從而具有較高的居里溫度和廣泛的應(yīng)用前景。

    2.3 自旋極化的局域態(tài)密度及分波態(tài)密度

    圖3為Fm0.125Zn0.875S的自旋極化的局域態(tài)密度(左圖)和分波態(tài)密度(右圖),圖中位于橫軸(DOS=0)上和下的曲線分別表示自旋向上和自旋向下的子帶。分波態(tài)密度包括s,p和d軌道上的態(tài)密度,局域態(tài)密度分別給出1個(gè)過(guò)渡離子及其最近鄰的1個(gè)S離子和Zn離子的電子態(tài)密度。因只有費(fèi)米面附近的子帶內(nèi)電子會(huì)對(duì)材料的電學(xué)和磁學(xué)性質(zhì)產(chǎn)生重要影響,本文只畫出了費(fèi)米面附近的子帶的態(tài)密度。由圖3可知,F(xiàn)m0.125Zn0.875S的第2、第3和第4子能帶中態(tài)密度的來(lái)源非常相似,第2子能帶 (在能量-8到-6 eV之間)主要來(lái)自S-3p和Fm-3p軌道,部分來(lái)自Zn-3d和Zn-4s軌道;第3子能帶(位于費(fèi)米面附近)主要由Fm-3d和S-3p軌道雜化而成,即材料中存在p-d雜化軌道,原因在于S-3p和Fm-3d軌道的能量相近,在強(qiáng)晶體場(chǎng)作用下易產(chǎn)生雜化。第四子能帶(能量遠(yuǎn)高于費(fèi)米面的空帶)主要來(lái)自于S-3p和Zn-4s軌道。不過(guò),由于Fe0.125Zn0.875S的晶體場(chǎng)比Co0.125Zn0.875S和Ni0.125Zn0.875S弱,所以Fe0.125Zn0.875S的第2和第3子能帶沒(méi)有重疊。而Co0.125Zn0.875S和Ni0.125Zn0.875S的第2和第3子帶發(fā)生重疊,且重疊部分隨過(guò)渡元素的序號(hào)增大而依次增大。

    為定量研究Fm0.125Zn0.875S的磁電性能,表2給出了Fm0.125Zn0.875S的2×2×1超胞中各離子的磁電參數(shù)。這些參數(shù)包括s、p、d軌道上帶的電子數(shù)、總電子數(shù)、離子電荷和離子磁矩。所有參數(shù)相同的離子列在同一行中。由表2,F(xiàn)e、Co和Ni離子磁矩分別為 3.62μB、2.46μB、1.32μB, 遠(yuǎn)大于 Zn 和 S 的離子磁矩。同時(shí),由圖3左邊3個(gè)圖可看出,Zn和S離子在費(fèi)米面處自旋向下子帶的態(tài)密度為零。這表明材料的總磁矩及其半金屬性主要來(lái)源于鐵磁性過(guò)渡金屬離子。將超胞中所有離子磁矩相加,可以得到超胞的總磁矩。Fe0.125Zn0.875S、Co0.125Zn0.875S和Ni0.125Zn0.875S超胞的總磁矩分別為 3.96μB、2.90μB和 2.00μB,近似等于整數(shù)磁矩 4.00μB、3.00μB、2.00μB。具有整數(shù)磁矩是半金屬鐵磁體一個(gè)非常重要的特征,這再次表明Fm0.125Zn0.875S具有半金屬性,這與由圖2和圖3的態(tài)密度給出的結(jié)論是一致的。

    表2 Fm0.125Zn0.875S(Fm=Fe、Co 和 Ni)中各離子的 s、p、d 軌道上的電子數(shù)(Ns,Np,Nd)、總電子數(shù)(Nt)、電荷數(shù)及離子磁矩Table 2 Ionic parameters of Fm0.125Zn0.875S(Fm=Fe,Co and Ni),including the number of s,p,d and total electrons(Ns,Np,Nd),the charges,ion magnetic moments

    本文以Fe0.125Zn0.875S為例來(lái)分析Fm0.125Zn0.875S出現(xiàn)半金屬性的機(jī)理。根據(jù)配位場(chǎng)理論[35-36],摻雜Fe原子后,中心Fe離子的周圍存在4個(gè)配體S離子,這4個(gè)配體離子與中心離子形成四面體配合物ML4。ML4內(nèi)存在強(qiáng)烈的晶體場(chǎng),從而使得過(guò)渡元素Fe的3d軌道發(fā)生能帶劈裂,形成能量較低的二重簡(jiǎn)并的eg成鍵軌道和三重簡(jiǎn)并的t2g成鍵軌道,以及它們對(duì)應(yīng)的能量較高的二重反鍵軌道e*g和三重反鍵軌道t*2g。Fe離子有6個(gè)3d電子,根據(jù)能量最低原理,它們將進(jìn)入eg成鍵軌道和三重簡(jiǎn)并的t2g成鍵軌道。由于Fe0.125Zn0.875S的晶體場(chǎng)不夠強(qiáng),使得電子的配對(duì)能高于eg和t2g的分離能,電子將按高自旋方式排列,因此Fe離子的電子結(jié)構(gòu)為eg2↑eg1↓t2g3↑。類似地,1個(gè)Co離子有7個(gè)3d電子,它們分別占據(jù)eg和t2g成鍵軌道,其電子結(jié)構(gòu)為eg2↑eg2↓t2g3↑。Ni離子3d軌道上有8個(gè)電子,其電子組態(tài)為eg2↑eg2↓t2g3↑t2g1↓。在這樣的電子結(jié)構(gòu)下,F(xiàn)e0.125Zn0.875S、Co0.125Zn0.875S和Ni0.125Zn0.875S超胞總磁矩分別為 4.00μB、3.00μB和 2.00μB,與表 2 給出的超胞磁矩一致。

    3 結(jié) 論

    摻雜后體系均在費(fèi)米面處產(chǎn)生-100%的自旋極化,具有半金屬性,是潛在的優(yōu)質(zhì)自旋注入材料。Fm0.125Zn0.875S可能具有較寬的自旋帶隙,從而具有較高的居里溫度和廣泛的應(yīng)用前景。Fe0.125Zn0.875S、Co0.125Zn0.875S和Ni0.125Zn0.875S的超胞磁矩分別為3.96、2.90和2.00 μB,主要來(lái)自于過(guò)渡金屬離子Fe、Co和Ni離子,它們的電子結(jié)構(gòu)分別為eg2↑eg1↓t2g3↑,eg2↑eg2↓t2g3↑和 eg2↑eg2↓t2g3↑t2g1↓。

    [1]de Groot R A,Mueller F M,van Engen P G,et al.Phys.Rev.Lett.,1983,50(25):2024-2027

    [2]Schmidt G,Molenkamp L W.Semicond.Sci.Tech.,2002,17(4):320-321

    [3]Baibich M N,Broto J M,Fert A,et al.Phys.Rev.Lett.,1988,61(21):2472-2475

    [4]Mathieu R,Akahoshi D,Asamitsu A,et al.Phys.Rev.Lett.,2004,93(22):227202(4 pages)

    [5]Aziz A,Wessely O P,Ali M,et al.Phys.Rev.Lett.,2009,103(23):237203(4 pages)

    [6]Miao G X,Müller M T,Moodera J S,et al.Phys.Rev.Lett.,2009,102(7):076601(4 pages)

    [7]Shi L J,Liu B G.Phys.Rev.B,2007,76(11):115201(9 pagers)

    [8]Hu X.Adv.Mater.,2012,24(2):294-298

    [9]Picozzi S,Continenza A,Freeman A.J.Phys.Rev.B,2002,66(9):094421(9 pagers)

    [10]Galanakis I.Phys.Rev.B,2005,71(1):012413(4 pagers)

    [11]Wurmehl S,Fecher G H,Kandpalet H C,et al.Appl.Phys.Lett.,2006,88(3):032503-032505

    [12]Watts S M,Wirth S,Molnár S V,et al.Phys.Rev.B,2000,61(14):9621-9628

    [13]Soeya S,Hayakawa J,Takahashi H,et al.Appl.Phys.Lett.,2002,80(5),823-845

    [14]Soulen Jr R J,Byers J M,Osofsky M S,et al.Science,1998,282(5386):85-88

    [15]Kobayashi K I,Kimura T,Sawada H,et al.Nature,1998,395:677-680

    [16]Zhang Z Y,Guo W.J.Phys.Chem.C,2012,116(1):893-900

    [17]Liu J,Chen L,Dong H N,et al.Appl.Phys.Lett.,2009,95(13):132502-132504

    [18]Giebultowicz T M,K?osowski P,Samarthet N,et al.Phys.Rev.B,1993,48(17):12817-12833

    [19]Saito H,Zayets V,Yamagata S,et al.Phys.Rev.Lett.,2003,90(20):207202(1-4)

    [20]Dong S J,Zhao H.Appl.Phys.Lett.,2011,98(18):182501-182503

    [21]Thangavel R,Rajagopalan M,Kumar J.J.Mag.Mag.Mater.,2008,320(5):774-778

    [22]GUO Guang-Sheng(郭廣生),LIU Ying-Rong(劉穎榮),WANG Zhi-Hua(王志華),et al.Chinese J.Inorg.Chem.(Wuji Huaxue Xuebao),2000,16(3):492-496

    [23]SUN Yuan-Guang(孫遠(yuǎn)光),CAO Li-Xin(曹立新),LIU Wei(柳偉),et al.J.Funct.Mater.(Gongneng Cailiao),2009,40(4):689-692

    [24]LI Sheng-Zhi(李勝旨),LIU Jin-Chao(劉錦超),YANG Xiang-Dong(楊向東),et al.Chinese J.High Press.Phys.(Gaoya Wuli Xuebao),2010,24(6):449-454

    [25]Bi C,Pan L Q,Xu M,et al.Mater.Chem.Phys.,2009,116(2/3):363-367

    [26]YANG Xu(楊 旭),ZHOU Hong(周 宏),SHEN Bin(沈 彬),et al.Acta Phys.-Chim.Sin.(Wuli Huaxue Xuebao),2010,26(1):244-248

    [27]Kullerud G.Norsk Geologisk Tidsskrift,1953,32:61-147

    [28]Becker W,Lutz H D.Mater.Res.Bull.,1978,13(13):907-911

    [29]Wu P,Kershaw R,Dwight K,et al.Mater.Res.Bull.,1989,24(1):49-53

    [30]Fan S W,Yao K L,Liu Z L.Appl.Phys.Lett.,2009,94(15):152506-152508

    [31]Xu Y N,Ching W Y.Phys.Rev.B,1993,48(7):4335-4351

    [32]Segall M D,Lindan P L D,Probert M J,et al.J.Phys.:Cond.Matt.,2002,14(11):2717-2743

    [33]Clark S J,Segall M D,Pickard C J,et al.Zeitschrift für Kristallographie,2005,220(5/6):567-570

    [34]Kisi E H,Elcombe M M.Acta Cryst.,1989,C45(12):1867-1870

    [35]Liu J,Chen X M,Liu Y,et al.Solid State Ionics,2008,179(21/22/23/24/25/26):881-886

    [36]Liu Y,Bose S K,Kudrnovsk J.Phys.Rev.B,2010,82(9):094435:1-8

    Influence of Ferromagnetic Transition Dopants on the Electronic Structures of Wurtzite ZnS

    DONG Hui-Ning*,1ZHANG Bo2CHEN Pei-Da1ZHANG Cui-Lin3LIU Jun*,1
    (1School of Mathematics and Physics,Chongqing University of Posts and Telecommunications,Chongqing 400065,China)
    (2School of Communication and Information Engineering,Chongqing University of Posts and Telecommunications,Chongqing 400065,China)
    (3Computer Science and Information Engineering College,Chongqing Technology and Business University,Chongqing 400067,China)

    The geometrical structures of Fm0.125Zn0.875S (Fm=Fe,Co and Ni)were optimized,and their electronic structures were investigated by the first principles calculations based on the density functional theory.The negative complete spin polarization at the Fermi level was found for these compounds.Therefore,they are potential half-metallic ferromagnets and spin injectors.They have wide spin gaps,so they have high Curie temperatures possibly and wide application prospects.The supercell magnetic moments of Fe0.125Zn0.875S,Co0.125Zn0.875S and Ni0.125Zn0.875S are 3.96μB,2.90μB,2.00μB,respectively,which arise mainly from doped ferromagnetic ions.The electronic structures of Fe,Co and Ni ions are,respectively.

    half-metallic ferromagnets;electronic structures;supercell magnetic moments

    O641

    A

    1001-4861(2012)07-1447-06

    2011-11-15。收修改稿日期:2012-03-01。

    重慶市自然科學(xué)基金(No.CSCT2010BB4405和CSTC2008BB4083)資助項(xiàng)目。

    *通訊聯(lián)系人。E-mail:liujun@cqupt.edu.cn,donghn@cqupt.edu.cn

    猜你喜歡
    金屬性鐵磁性價(jià)帶
    外加電場(chǎng)和磁場(chǎng)對(duì)外爾半金屬能帶結(jié)構(gòu)的調(diào)控*
    鐵磁性物質(zhì)對(duì)地磁觀測(cè)影響的野外測(cè)試
    大提離下脈沖渦流對(duì)鐵磁性材料測(cè)厚研究
    交換場(chǎng)和非共振光對(duì)單層MoS2能帶結(jié)構(gòu)的調(diào)控*
    淺談金屬性與金屬活動(dòng)性
    Modification of CaO-based sorbents prepared from calcium acetate for CO2 capture at high temperature☆
    元素金屬性、非金屬性強(qiáng)弱的判斷依據(jù)
    新課程(下)(2016年5期)2016-08-15 02:57:43
    核電站鐵磁性高加管漏磁檢測(cè)技術(shù)淺析
    科技視界(2015年30期)2015-10-22 11:26:44
    日本發(fā)現(xiàn)鐵系超導(dǎo)材料的超導(dǎo)狀態(tài)和反鐵磁性狀態(tài)可以同時(shí)存在
    浙江電力(2014年6期)2014-01-27 03:31:37
    關(guān)于元素金屬性與非金屬性的討論
    国产真实伦视频高清在线观看| 日本撒尿小便嘘嘘汇集6| 国产一区二区在线av高清观看| 久久热精品热| 欧美区成人在线视频| 国产亚洲精品av在线| 久久久久久久久久久丰满| 99热全是精品| 可以在线观看的亚洲视频| 精品久久国产蜜桃| 国产一区二区在线观看日韩| 精品一区二区三区视频在线观看免费| 国产探花极品一区二区| 看黄色毛片网站| 成人高潮视频无遮挡免费网站| 麻豆久久精品国产亚洲av| 22中文网久久字幕| 美女大奶头视频| 久久这里只有精品中国| 亚洲美女视频黄频| ponron亚洲| 在线观看66精品国产| 精品福利观看| 两性午夜刺激爽爽歪歪视频在线观看| 午夜精品一区二区三区免费看| 亚洲人成网站在线播| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品乱码久久久v下载方式| 成人性生交大片免费视频hd| 男女视频在线观看网站免费| 一级毛片电影观看 | 听说在线观看完整版免费高清| 亚洲高清免费不卡视频| 亚洲精品久久国产高清桃花| 少妇人妻精品综合一区二区 | 国产一区二区亚洲精品在线观看| 亚洲精品国产av成人精品 | 九九热线精品视视频播放| 丰满的人妻完整版| 亚洲精品日韩av片在线观看| 色综合亚洲欧美另类图片| 久久精品综合一区二区三区| 毛片一级片免费看久久久久| 国产一区二区激情短视频| 亚洲一区二区三区色噜噜| www日本黄色视频网| 狂野欧美激情性xxxx在线观看| 一本久久中文字幕| 亚洲av不卡在线观看| 色哟哟·www| 少妇人妻一区二区三区视频| www日本黄色视频网| 精品国产三级普通话版| 亚洲欧美日韩无卡精品| 97人妻精品一区二区三区麻豆| 在线观看av片永久免费下载| 天天躁夜夜躁狠狠久久av| 国产高潮美女av| av在线蜜桃| 亚洲欧美精品自产自拍| 深夜a级毛片| 国产精品一及| 国产一区二区激情短视频| 国产精品爽爽va在线观看网站| 免费看美女性在线毛片视频| 国产精品爽爽va在线观看网站| 精华霜和精华液先用哪个| 韩国av在线不卡| 国产精品久久久久久av不卡| 国产高清有码在线观看视频| 少妇裸体淫交视频免费看高清| 国产一区二区在线av高清观看| 变态另类丝袜制服| 美女xxoo啪啪120秒动态图| 欧美性猛交╳xxx乱大交人| 国产黄色小视频在线观看| 亚洲精品国产成人久久av| 日韩欧美三级三区| 丰满的人妻完整版| 成人三级黄色视频| 日韩av不卡免费在线播放| 麻豆久久精品国产亚洲av| 亚洲一区二区三区色噜噜| 久久久久国产精品人妻aⅴ院| 欧美激情国产日韩精品一区| 亚洲国产日韩欧美精品在线观看| 欧美一级a爱片免费观看看| 国产综合懂色| 别揉我奶头 嗯啊视频| 啦啦啦观看免费观看视频高清| 可以在线观看毛片的网站| 国产欧美日韩精品一区二区| 人妻丰满熟妇av一区二区三区| 日本在线视频免费播放| 亚洲天堂国产精品一区在线| 日本三级黄在线观看| 欧美成人精品欧美一级黄| 欧美3d第一页| 99热这里只有是精品在线观看| 秋霞在线观看毛片| 日本撒尿小便嘘嘘汇集6| 国产成人aa在线观看| 国产av不卡久久| 国产精品人妻久久久影院| 久久久成人免费电影| 草草在线视频免费看| 国产成人福利小说| 精品久久国产蜜桃| 插阴视频在线观看视频| 女人被狂操c到高潮| 真实男女啪啪啪动态图| 午夜久久久久精精品| 人妻丰满熟妇av一区二区三区| 亚洲丝袜综合中文字幕| 亚洲av二区三区四区| 波野结衣二区三区在线| 一级黄色大片毛片| 国产熟女欧美一区二区| 亚洲人成网站在线播放欧美日韩| 女同久久另类99精品国产91| 国产成人精品久久久久久| 男人舔奶头视频| 久久亚洲国产成人精品v| 精品一区二区三区av网在线观看| 一级毛片aaaaaa免费看小| av专区在线播放| 欧美绝顶高潮抽搐喷水| 国产精品一区二区性色av| 五月伊人婷婷丁香| 俄罗斯特黄特色一大片| 久99久视频精品免费| 欧美成人精品欧美一级黄| 联通29元200g的流量卡| 伦精品一区二区三区| 欧美日韩一区二区视频在线观看视频在线 | 国产亚洲精品久久久com| 欧美丝袜亚洲另类| 真实男女啪啪啪动态图| 熟妇人妻久久中文字幕3abv| 99久久久亚洲精品蜜臀av| 香蕉av资源在线| 久久综合国产亚洲精品| 免费观看的影片在线观看| 欧美国产日韩亚洲一区| 久久热精品热| videossex国产| 99热精品在线国产| 国产三级在线视频| 国产精品久久久久久av不卡| 亚洲中文日韩欧美视频| 国产精品av视频在线免费观看| 性色avwww在线观看| 欧美潮喷喷水| 97人妻精品一区二区三区麻豆| 成人无遮挡网站| 精品国内亚洲2022精品成人| 在线观看一区二区三区| 亚洲人与动物交配视频| av福利片在线观看| 最近视频中文字幕2019在线8| 欧美另类亚洲清纯唯美| 欧美+亚洲+日韩+国产| 日韩成人av中文字幕在线观看 | 亚洲av成人av| 欧美成人精品欧美一级黄| 插逼视频在线观看| 成人漫画全彩无遮挡| 欧美高清性xxxxhd video| av.在线天堂| 九九热线精品视视频播放| 色av中文字幕| 国产高清视频在线播放一区| 丝袜美腿在线中文| 欧美另类亚洲清纯唯美| 精品少妇黑人巨大在线播放 | 婷婷亚洲欧美| 91久久精品国产一区二区三区| 成年女人毛片免费观看观看9| 欧美区成人在线视频| 久久久精品94久久精品| 老女人水多毛片| 日日啪夜夜撸| 欧美xxxx黑人xx丫x性爽| 可以在线观看毛片的网站| 18禁黄网站禁片免费观看直播| 久久天躁狠狠躁夜夜2o2o| 免费一级毛片在线播放高清视频| 色综合亚洲欧美另类图片| 国产麻豆成人av免费视频| 婷婷亚洲欧美| 欧美性感艳星| 一级毛片aaaaaa免费看小| 3wmmmm亚洲av在线观看| 秋霞在线观看毛片| 亚洲中文字幕一区二区三区有码在线看| 99九九线精品视频在线观看视频| 日本一本二区三区精品| 国产精品无大码| 午夜福利视频1000在线观看| 国产毛片a区久久久久| 国产精品人妻久久久影院| 免费看av在线观看网站| 精品一区二区三区视频在线观看免费| 99九九线精品视频在线观看视频| 国内精品久久久久精免费| 日韩欧美精品v在线| 亚洲精品影视一区二区三区av| 欧美国产日韩亚洲一区| 久久久久国产精品人妻aⅴ院| www日本黄色视频网| 偷拍熟女少妇极品色| 亚洲国产精品成人久久小说 | 男人狂女人下面高潮的视频| 日本在线视频免费播放| 国产亚洲精品久久久久久毛片| 国产精品,欧美在线| 美女xxoo啪啪120秒动态图| 在线看三级毛片| 97碰自拍视频| 热99re8久久精品国产| 国产精品永久免费网站| 在线国产一区二区在线| 精品人妻偷拍中文字幕| 色综合亚洲欧美另类图片| 我的女老师完整版在线观看| 亚洲最大成人av| 中国美女看黄片| 久久人人爽人人片av| 国产久久久一区二区三区| 熟妇人妻久久中文字幕3abv| 欧美又色又爽又黄视频| 免费高清视频大片| 天美传媒精品一区二区| 国产真实乱freesex| 色av中文字幕| 午夜视频国产福利| 高清午夜精品一区二区三区 | 国产精品乱码一区二三区的特点| 国产精品亚洲一级av第二区| 日本-黄色视频高清免费观看| 日韩欧美一区二区三区在线观看| 最近视频中文字幕2019在线8| 亚洲欧美清纯卡通| 桃色一区二区三区在线观看| 级片在线观看| 熟妇人妻久久中文字幕3abv| 国产三级在线视频| 中国国产av一级| 高清毛片免费观看视频网站| 国产单亲对白刺激| 99久久精品一区二区三区| 看十八女毛片水多多多| 亚洲国产精品sss在线观看| 久久久久久久久久黄片| 男女视频在线观看网站免费| 一级黄片播放器| 日本一本二区三区精品| 一个人看的www免费观看视频| av天堂在线播放| 亚州av有码| 插阴视频在线观看视频| 日本三级黄在线观看| 精品少妇黑人巨大在线播放 | 国产单亲对白刺激| 天堂av国产一区二区熟女人妻| 小说图片视频综合网站| 久久精品国产自在天天线| 插阴视频在线观看视频| 成人精品一区二区免费| 欧美日韩精品成人综合77777| 久久6这里有精品| 午夜爱爱视频在线播放| 秋霞在线观看毛片| 国产aⅴ精品一区二区三区波| 人人妻人人澡人人爽人人夜夜 | 国产精品人妻久久久久久| 99视频精品全部免费 在线| 亚洲成人中文字幕在线播放| 变态另类丝袜制服| 97热精品久久久久久| 亚洲av电影不卡..在线观看| 欧美人与善性xxx| 精品人妻偷拍中文字幕| 九色成人免费人妻av| 亚洲国产精品成人综合色| 精品午夜福利视频在线观看一区| 在线播放无遮挡| 国产高潮美女av| 日韩三级伦理在线观看| av在线播放精品| 精品午夜福利视频在线观看一区| av在线亚洲专区| 国产高清视频在线播放一区| 在线免费观看的www视频| 亚洲婷婷狠狠爱综合网| 国产精品永久免费网站| 国产高清不卡午夜福利| 色噜噜av男人的天堂激情| 午夜亚洲福利在线播放| 麻豆乱淫一区二区| 成人一区二区视频在线观看| 波野结衣二区三区在线| av国产免费在线观看| 国产一区二区激情短视频| 毛片女人毛片| 真实男女啪啪啪动态图| 色综合色国产| 18禁黄网站禁片免费观看直播| 亚洲精品粉嫩美女一区| 在线免费十八禁| 免费观看精品视频网站| 日本黄色视频三级网站网址| 久久国产乱子免费精品| 亚洲最大成人av| 成人av在线播放网站| 欧美一级a爱片免费观看看| .国产精品久久| 婷婷色综合大香蕉| 国产黄色小视频在线观看| 国产精品嫩草影院av在线观看| 99热这里只有是精品50| 日韩人妻高清精品专区| 最近视频中文字幕2019在线8| 51国产日韩欧美| 国产精品乱码一区二三区的特点| 国产精品av视频在线免费观看| 日本一二三区视频观看| 啦啦啦观看免费观看视频高清| 婷婷精品国产亚洲av在线| 亚洲av熟女| 高清毛片免费看| 一个人看视频在线观看www免费| 日韩精品青青久久久久久| 国产精品国产三级国产av玫瑰| 久久午夜亚洲精品久久| 亚洲av中文字字幕乱码综合| 国产爱豆传媒在线观看| 亚洲熟妇熟女久久| 国产亚洲av嫩草精品影院| 观看免费一级毛片| 婷婷精品国产亚洲av| 最新中文字幕久久久久| 一个人免费在线观看电影| 久久精品国产亚洲网站| 欧美潮喷喷水| 精品人妻偷拍中文字幕| 国产高清视频在线观看网站| 久久久久国内视频| 九九爱精品视频在线观看| 免费看光身美女| 色在线成人网| 乱人视频在线观看| 欧美极品一区二区三区四区| 日本黄色视频三级网站网址| 在线观看美女被高潮喷水网站| 99久久成人亚洲精品观看| 自拍偷自拍亚洲精品老妇| 国产成人影院久久av| 啦啦啦韩国在线观看视频| 日本五十路高清| 69人妻影院| 久久人妻av系列| 秋霞在线观看毛片| 亚洲三级黄色毛片| 嫩草影视91久久| 狠狠狠狠99中文字幕| 内地一区二区视频在线| 波野结衣二区三区在线| 国产精品国产三级国产av玫瑰| 男插女下体视频免费在线播放| 国产在线精品亚洲第一网站| 国产高清视频在线播放一区| 亚洲aⅴ乱码一区二区在线播放| 啦啦啦韩国在线观看视频| 看非洲黑人一级黄片| 最后的刺客免费高清国语| 成人美女网站在线观看视频| 男女啪啪激烈高潮av片| 午夜亚洲福利在线播放| 乱码一卡2卡4卡精品| 久久久久久九九精品二区国产| 91久久精品国产一区二区三区| 啦啦啦啦在线视频资源| 日韩精品中文字幕看吧| 亚洲欧美成人综合另类久久久 | 欧美潮喷喷水| 午夜福利在线在线| 岛国在线免费视频观看| 中文在线观看免费www的网站| 亚洲第一区二区三区不卡| 国产精品一区二区三区四区久久| 简卡轻食公司| 亚洲无线观看免费| 亚洲美女黄片视频| 丝袜美腿在线中文| av在线老鸭窝| ponron亚洲| 日韩欧美免费精品| 午夜精品国产一区二区电影 | 国产蜜桃级精品一区二区三区| 国产成人福利小说| 99久国产av精品国产电影| 亚洲精品乱码久久久v下载方式| 国产在线男女| 国产真实伦视频高清在线观看| 欧美zozozo另类| 中文字幕人妻熟人妻熟丝袜美| 12—13女人毛片做爰片一| 亚洲性夜色夜夜综合| 午夜福利成人在线免费观看| 高清午夜精品一区二区三区 | 波多野结衣高清作品| 午夜激情欧美在线| 精品少妇黑人巨大在线播放 | 精品一区二区三区av网在线观看| 国产三级中文精品| 丝袜美腿在线中文| 老女人水多毛片| 男人舔奶头视频| 一区二区三区四区激情视频 | 女的被弄到高潮叫床怎么办| 国产精品久久久久久久电影| 国产精品三级大全| 99久国产av精品| 99久久久亚洲精品蜜臀av| 偷拍熟女少妇极品色| 亚洲久久久久久中文字幕| 你懂的网址亚洲精品在线观看 | 精品国内亚洲2022精品成人| 婷婷精品国产亚洲av在线| 欧美成人一区二区免费高清观看| 欧美一级a爱片免费观看看| 欧美极品一区二区三区四区| 在线a可以看的网站| 两个人视频免费观看高清| 天堂av国产一区二区熟女人妻| 亚洲一级一片aⅴ在线观看| 老司机福利观看| av中文乱码字幕在线| 国产精品美女特级片免费视频播放器| 亚洲电影在线观看av| 免费电影在线观看免费观看| 亚洲国产色片| 网址你懂的国产日韩在线| 国产成人福利小说| 久久精品国产亚洲av天美| 蜜臀久久99精品久久宅男| 一本一本综合久久| 亚洲av五月六月丁香网| or卡值多少钱| 国产在视频线在精品| 此物有八面人人有两片| 性色avwww在线观看| 国产成人一区二区在线| 国产精品一区二区三区四区久久| 老司机午夜福利在线观看视频| 日韩欧美 国产精品| 成人午夜高清在线视频| 啦啦啦韩国在线观看视频| 久久精品国产自在天天线| 久久久国产成人免费| 成年版毛片免费区| 久久久久久国产a免费观看| 中文在线观看免费www的网站| 欧美日韩在线观看h| 亚洲四区av| 看免费成人av毛片| 国产国拍精品亚洲av在线观看| 亚洲久久久久久中文字幕| 嫩草影视91久久| 成人无遮挡网站| 国产成人freesex在线 | 内地一区二区视频在线| 婷婷精品国产亚洲av在线| 亚洲精品国产av成人精品 | 久久久久久大精品| 亚洲国产精品国产精品| 日本五十路高清| 日本欧美国产在线视频| 国内精品一区二区在线观看| 中文字幕熟女人妻在线| 国产精品一区二区免费欧美| 深夜a级毛片| 99久国产av精品国产电影| 伦理电影大哥的女人| 插逼视频在线观看| 精品久久国产蜜桃| 国产aⅴ精品一区二区三区波| 18禁黄网站禁片免费观看直播| 乱系列少妇在线播放| 精品不卡国产一区二区三区| 一夜夜www| 国产欧美日韩精品一区二区| 亚洲av.av天堂| 精品日产1卡2卡| 欧美一区二区精品小视频在线| 高清毛片免费观看视频网站| 看免费成人av毛片| 狂野欧美激情性xxxx在线观看| 99精品在免费线老司机午夜| 99久久精品热视频| 男人狂女人下面高潮的视频| 久久久午夜欧美精品| 午夜爱爱视频在线播放| 亚洲中文字幕日韩| 亚洲婷婷狠狠爱综合网| 国产高清不卡午夜福利| 禁无遮挡网站| 成人av一区二区三区在线看| 久99久视频精品免费| ponron亚洲| 精品人妻熟女av久视频| 精品一区二区三区人妻视频| 高清毛片免费观看视频网站| 日本成人三级电影网站| 一夜夜www| 久久久久免费精品人妻一区二区| 少妇的逼好多水| 精品久久久久久久久久久久久| 一区福利在线观看| 别揉我奶头 嗯啊视频| 国产精品av视频在线免费观看| 国产精品久久久久久久久免| 国产三级中文精品| 亚洲成人久久爱视频| 国产女主播在线喷水免费视频网站 | 看非洲黑人一级黄片| 国产日本99.免费观看| 精品一区二区三区人妻视频| 久久天躁狠狠躁夜夜2o2o| 熟妇人妻久久中文字幕3abv| 欧美性猛交╳xxx乱大交人| 91在线观看av| 色综合站精品国产| 成人美女网站在线观看视频| 丰满的人妻完整版| 晚上一个人看的免费电影| 日产精品乱码卡一卡2卡三| 国产亚洲精品久久久com| 久久精品91蜜桃| 最后的刺客免费高清国语| 国产高清不卡午夜福利| 欧美zozozo另类| 国产免费男女视频| 我要看日韩黄色一级片| 看片在线看免费视频| 99久久中文字幕三级久久日本| 成年av动漫网址| av国产免费在线观看| 男人狂女人下面高潮的视频| 亚洲精品日韩在线中文字幕 | 一级毛片我不卡| 免费av不卡在线播放| 亚洲精品久久国产高清桃花| 欧美最新免费一区二区三区| 国产精品久久视频播放| 国语自产精品视频在线第100页| av卡一久久| 婷婷亚洲欧美| 日韩av在线大香蕉| 欧美最新免费一区二区三区| 日韩欧美精品免费久久| 亚洲第一区二区三区不卡| av女优亚洲男人天堂| 日韩人妻高清精品专区| 国产黄色视频一区二区在线观看 | 午夜老司机福利剧场| av天堂中文字幕网| 六月丁香七月| 精品福利观看| 99国产精品一区二区蜜桃av| 国产成人freesex在线 | 国产亚洲精品av在线| 欧美一区二区亚洲| 亚洲,欧美,日韩| 观看美女的网站| 日日摸夜夜添夜夜添av毛片| 国产一区二区三区在线臀色熟女| 波多野结衣巨乳人妻| 亚洲国产欧美人成| 啦啦啦韩国在线观看视频| 久久人人爽人人爽人人片va| 亚洲三级黄色毛片| АⅤ资源中文在线天堂| 天堂av国产一区二区熟女人妻| 日韩强制内射视频| 中国美白少妇内射xxxbb| 黄色视频,在线免费观看| 亚洲国产日韩欧美精品在线观看| 狂野欧美激情性xxxx在线观看| 真人做人爱边吃奶动态| 国产成人福利小说| a级一级毛片免费在线观看| 国产一区二区亚洲精品在线观看| 成人欧美大片| 丰满乱子伦码专区| 尾随美女入室| 国产成人a∨麻豆精品| 国产精品久久久久久久电影| 国产黄色视频一区二区在线观看 | 精品午夜福利视频在线观看一区| www日本黄色视频网| 成年av动漫网址| 精品一区二区免费观看| 国产91av在线免费观看| 国产高清激情床上av| 观看免费一级毛片| 国产高清视频在线观看网站| 精品国内亚洲2022精品成人| 国产一区二区激情短视频| 欧美高清成人免费视频www| 一夜夜www| 一级黄色大片毛片| 别揉我奶头~嗯~啊~动态视频| 日本与韩国留学比较|