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

    典型磁性材料價電子結構研究面臨的機遇與挑戰(zhàn)*

    2020-02-18 03:17:32唐貴德李壯志馬麗吳光恒胡鳳霞
    物理學報 2020年2期
    關鍵詞:電離能巡游氧化物

    唐貴德 李壯志 馬麗 吳光恒 胡鳳霞

    1) (中國科學院物理研究所,磁學國家重點實驗室,北京 100190)

    2) (河北師范大學物理學院,河北省新型薄膜材料實驗室,石家莊 050024)

    目前在磁性材料磁有序現(xiàn)象研究中廣泛使用的交換作用、超交換作用和雙交換作用模型形成于1950年代及其以前,這些模型都涉及材料中的價電子狀態(tài),但那時還沒有充分的價電子狀態(tài)實驗依據(jù).1970年代以來,有關價電子結構實驗結果的報道越來越多,這些實驗結果表明傳統(tǒng)的磁有序模型需要改進.首先,大量電子譜實驗表明,在氧化物中除存在負二價氧離子之外,還存在負一價氧離子,并且負一價氧離子的含量可達30%或更多.這說明以所有氧離子都是負二價離子為基本假設的超交換和雙交換作用模型需要改進.其次,一些實驗證明,鐵、鈷、鎳自由原子的一部分4s電子在形成鐵磁性金屬的過程中變成了3d電子,這為探討金屬磁性與電輸運性質的關系提供了依據(jù).此外,即使在現(xiàn)代的密度泛函計算中,仍不能給出磁性交換作用能的函數(shù)表達式,只能采取各種不同模型進行模擬計算,從而使磁性材料的模擬計算遇到嚴重困難.尋求一個磁有序能的函數(shù)表達式可能是解決這個困難的途徑.這些研究表明磁性材料價電子結構研究面臨著重大的機遇與挑戰(zhàn).本文首先介紹一些典型的實驗例證,然后介紹了基于這些實驗結果的一套典型磁性材料的磁有序新模型,隨后介紹了基于新模型的磁性材料價電子結構與舊模型的主要區(qū)別,最后指出了未來研究工作面臨的挑戰(zhàn).

    1 引 言

    對于鐵磁性材料磁疇中原子磁矩的有序排列,外斯首先提出存在一個分子場,這個分子場能量是磁有序能的來源[1].為解釋這個分子場的物理機制,逐漸形成了關于金屬磁性的直接交換作用模型[2,3]、關于氧化物反鐵磁性的超交換相互作用(superexchange,SE)模型[4]以及關于氧化物中陽離子磁矩鐵磁耦合的雙交換相互作用(double exchange,DE)模型[5].這些形成于1950年代及以前的模型都涉及到磁性材料的價電子狀態(tài).然而,那時還沒有充分的價電子狀態(tài)實驗依據(jù).1970年代以來,有關價電子結構實驗結果的報道越來越多.這些實驗結果表明傳統(tǒng)的磁有序模型需要改進.這是21世紀磁學和磁性材料研究面臨的重要機遇.

    本研究組在改進磁有序模型方面進行了一些嘗試,提出一個關于氧化物磁有序的O 2p巡游電子模型,一個關于金屬磁有序的新巡游電子模型和一個關于外斯分子場能量的外斯電子對模型,發(fā)表了一系列研究論文,相應的綜述文章于2018年發(fā)表在《Physics Reports》上[6].此外,基于這些磁有序模型,還對典型磁性材料的電輸運性質進行了探討[7-9],認為NiCu合金與典型鈣鈦礦結構錳氧化物在居里溫度以下都存在巡游電子的自旋相關導電.這些新模型中的價電子結構與傳統(tǒng)模型存在顯著區(qū)別.本文首先簡單介紹文獻報道的一些相關典型實驗結果,然后介紹本研究組的有關工作,最后指出這方面研究面臨的挑戰(zhàn).

    2 價電子狀態(tài)的典型實驗例證

    2.1 O 2p空穴和負一價氧離子的實驗研究

    在包括磁性材料在內(nèi)的氧化物材料研究中,傳統(tǒng)觀點認為所有的氧離子都是負二價離子,其最外電子殼層是8個電子(2s22p6)的滿殼層結構.然而實驗表明在氧化物中存在O 2p空穴,即存在負一價氧離子(2s22p5),其比例隨材料不同而有很大差別,可達30%或更高.說明基于氧離子全部為負二價離子假設的傳統(tǒng)超交換作用和雙交換作用模型需要改進.

    2.1.1 關于O 2p空穴的實驗研究

    1988 年,Nücker等[10]報道了關于超導材料YBa2Cu3O7—y的電子能量損失譜(electron energy loss spectroscopy,EELS)研究,發(fā)現(xiàn)在費米能級處的電子態(tài)密度具有顯著的O 2p空穴特征,并且O 2p空穴含量隨材料中氧離子空位含量y而變化,所以認為這種超導化合物的載流子是O 2p空穴.

    Ju等[11]利用聚合物溶膠-凝膠技術,在(100)面的LaAlO3基片上制備了鈣鈦礦錳氧化物薄膜材料 La1—xSrxMnO3(0≤x≤0.7).樣品最后在空氣中700 ℃退火1 h.通過X射線衍射分析,顯示所有樣品都具有單相ABO3鈣鈦礦結構.測量了樣品電阻率隨測試溫度的變化關系,得到的結果如圖1(a)所示.曲線的變化規(guī)律與Urushibara等[12]研究單晶樣品 La1—xSrxMnO3(0≤x≤0.4)時得到的結果十分相似.圖1(b)為La1—xSrxMnO3(x=0,0.3,0.5) 樣品薄膜和LaAlO3基片的O K邊電子能量損失譜,其橫軸已經(jīng)換算成樣品發(fā)射光電子相對于費米能級的束縛能.可以看到3個薄膜樣品在529 eV 附近有一個很強的譜峰,而LaAlO3基片沒有這個譜峰.這個譜峰反映出樣品中存在O 2p的空態(tài),說明La1—xSrxMnO3系列樣品的導電是由于O 2p空穴造成的p型導電.

    Mizoroki等[13]制備了多晶樣品 La1—xSrxMnO3(x=0.1,0.2,0.3,04,0.5),利用磁性康普頓散射實驗證明,當Sr摻雜量為0.1和0.2時,摻雜的空穴擇優(yōu)進入O 2p 態(tài).Grenier等[14]利用O K邊X射線衍射,在La7/8Sr1/8MnO3中觀察到2p 電荷的有序變化,即存在“富空穴”和“貧空穴”MnO平面的交替變化.Ibrahim等[15]制備了龐磁電阻材料Pr1—xSrxMnO3(x=0,0.3),研究了材料的X射線吸收譜和俄歇電子譜,認為材料中存在O 2p空穴,且其濃度隨Sr摻雜量的增加而增大.Papavassiliou等[16]比較了材料 La1—xCaxMnO3+δ的核磁共振結果與X射線吸收數(shù)據(jù),發(fā)現(xiàn)在O 2p軌道上形成了自旋極化的空穴.

    圖1 由Ju等[11]報道的 La1—xSrxMnO3系列樣品 (a)電阻率ρ隨測試溫度T的變化關系;(b)電子能量損失譜Fig.1.La1—xSrxMnO3 reported by Ju et al.[11]:(a) Curves of the resistivity ρ versus the test temperature T;(b) electron energy loss spectra.

    2.1.2 關于負一價氧離子的實驗和理論研究

    Dupin等[17]研究了一些氧化物的X射線光電子譜(X-ray photoelectron spectroscopy,XPS),提出一個觀點:在氧化物的O 1s譜中出現(xiàn)在527.7—530.6 eV范圍內(nèi)的譜峰對應于O2—離子的光電子,在531.1—532.0 eV 范圍內(nèi)的譜峰對應于O—離子的光電子;與材料表面化學吸附氧(O0)對應的譜峰出現(xiàn)在531.1—533.5 eV 范圍內(nèi),如圖2所示.

    圖2 Dupin等[17]提出的O 1s譜峰所對應的氧離子價態(tài)示意圖Fig.2.A binding energy scale for valence state of oxygen at the O 1s peaks,proposed by Dupin et al.[17].

    Cohen等[18,19]利用密度泛函理論(density functional theory,DFT)計算了BaTiO3的價電子態(tài)密度,計算結果給出Ba,Ti和O 的平均化合價分別是2.00,2.89 和—1.63.這個計算結果得到了XPS實驗結果的證明.分析一些氧化物的XPS,包括BaTiO3,一氧化物CaO,MnO,CoO,ZnO,NiO,CuO粉末樣品[20],一個SrTiO3多晶和一個SrTiO3單晶片體樣品[21],得到 BaTiO3粉末、SrTiO3多晶和SrTiO3單晶樣品中氧離子的平均化合價為—1.55,—1.62 和—1.76,與 Cohen 的計算結果—1.63接近.此外,發(fā)現(xiàn)一氧化物中氧的平均化合價絕對值隨其中陽離子第二電離能(V2)的增大而減小,即其中負二價氧離子的含量比例隨V2的增大而減小[20].其原因在于氧的第二電子親和能(8.08 eV)小于陽離子的第二電離能,所以氧離子要從陽離子得到第二個電子存在不同程度的困難.

    2.2 金屬的價電子結構研究

    實驗表明,在3d過渡金屬Fe和Co從自由原子形成金屬晶體的過程中,一部分4s電子進入3d軌道,變成了3d電子;絕大部分3d電子是局域電子.說明用能帶論處理金屬磁有序問題時,把全部3d和4s都看成巡游電子的假設[1]需要改進.

    1995 年,美國學者Chen等[22]在超高真空條件下制備了厚度為50—70 ?的Fe和Co金屬薄膜,原位測量了樣品的L2,3吸收邊的X射線吸收譜(X-ray absorption spectra,XAS)和X射線圓二色譜(X-ray magnetic circular dichroism,XMCD).Chen等的分析給出兩個值得注意的信息:1) Fe和Co薄膜中平均每個原子的自旋磁矩分別為1.98μB和 1.55μB,軌道磁矩分別為 0.086μB和0.153μB.軌道磁矩分別為自旋磁矩的4.3%和9.9%,說明軌道磁矩的貢獻很小;2) 分析過程中,對于Fe和Co的3d電子平均數(shù)目,采用了理論計算結果6.61和7.51[23,24],明顯多于自由原子的3d電子數(shù)目6和7,說明在從自由原子形成金屬的過程中,一部分4s電子進入3d軌道,變成了3d電子.

    2007 年,德國學者 Jauch和 Reehuis[25]使用316.5 keV的γ射線,在295 K下測量了α-Fe單晶的高精度結構因子.通過對這個結構因子的分析,認為在α-Fe中的價電子組態(tài)應為3d74s1,而不是自由Fe原子的3d64s2.

    2015 年,瑞士學者Pacchioni等[26]研究了吸附在Cu(111)面上的Fe單原子和小團簇的磁學性質.在一個真空腔中,利用氬離子濺射制備了(111)面的Cu單晶.然后,利用電子束在原位蒸發(fā)Fe.在蒸發(fā)過程中,Cu基片的溫度為3.5 K,以保證Fe原子是吸附在Cu(111)面上,而不存在與Cu基片的較強相互作用.與制備樣品的真空腔連接有掃描隧道顯微鏡,用以原位檢測Fe層的厚度,并且可在不破壞真空的條件下,把樣品轉移至XAS和XMCD測量室進行測量.他們制備了系列樣品,Fe的覆蓋范圍在0.007—0.145單層(ML)之間.其中ML定義為每個Cu 晶胞的(111)面上有一個Fe原子.他們首先研究了0.007 ML的樣品,這相當于Fe原子為吸附在Cu(111)面上的孤立原子.XAS測量在溫度T=2.5 K 和磁場 B=6.8 T 條件下進行.通過對XAS和XMCD實驗結果的分析,得到吸附在Cu(111)面上的孤立Fe原子3d次殼層的3d空穴數(shù),hd=3.03,價電子結構近似為3d74s1.相對于自由的Fe原子,相當于有一個4s電子進入了3d軌道.對于吸附在Cu(111)面上的Fe原子團簇,隨著團簇中平均原子數(shù)目逐漸增加,平均每個原子的軌道磁矩迅速減小,從孤立Fe原子時的0.66μB減小到 5個Fe原子時的0.2μB.而由電子自旋形成的平均原子磁矩在5個Fe原子時為2.4μB,與公認的塊體金屬Fe的平均原子磁矩2.22μB非常接近.平均軌道磁矩與自旋磁矩比值為8.3%.

    3 改進磁有序模型的探索及相關研究

    上述實驗和理論結果為改進磁有序模型提供了重要依據(jù),而磁有序模型的改進是磁性材料價電子結構研究邁向新發(fā)展階段的重要機遇,在這方面已經(jīng)有一些嘗試,但是還遠沒有得到廣泛的重視.

    3.1 O 2p空穴在高溫超導和鈣鈦礦結構錳氧化物研究中的應用

    基于Nücker等[10]的實驗結果以及相應的其他一些實驗結果,在高溫超導性質研究方面,已經(jīng)考慮到O 2p空穴的影響.1998年韓汝珊[27]所著《高溫超導物理》一書在介紹用固體能帶論研究超導體氧化物時,給出一個Hubbard模型三帶形式的能量表達式,其中就包含了與O 2p空穴相應的能量項.

    對于鈣鈦礦錳氧化物 La1—xSrxMnO3(0≤x≤0.4),傳統(tǒng)的磁有序理論用雙交換作用(DE)模型和超交換作用(SE)模型解釋其中Sr含量從0到0.4變化過程中電輸運性質的變化,認為其中所有氧離子都處于負二價,在三價和四價Mn離子之間躍遷的Mn 3d電子是其中電流載流子的來源[5,12].然而基于上述電子能量損失譜等實驗結果[10-16],Alexandrov等[28,29]指出,傳統(tǒng)的磁耦合DE模型與這些實驗結果相沖突,這些實驗結果清楚地表明鐵磁錳氧化物中的電流載流子是氧的p空穴而不是d電子.從而提出一個完全不同于DE模型的O 2p空穴載流子模型,用于解釋鈣鈦礦錳氧化物的導電性質.

    3.2 磁有序新模型探索

    基于上述氧化物和金屬磁性材料的實驗和理論研究結果,針對傳統(tǒng)鐵磁學在解釋一些典型磁性材料磁有序問題時遇到的困難,提出了一個關于氧化物磁有序的O 2p巡游電子(itinerant electron model for magnetic oxides,IEO)模型[6,30,31]、一個解釋金屬中平均原子磁矩與電阻率關系的金屬磁有序巡游電子(itinerant electron model for magnetic metals,IEM)模型[6,32]、以及一個關于解釋氧化物和金屬磁有序能來源的外斯電子對(WEP)模型[6,33].

    3.2.1 IEO模型

    1) 磁性氧化物中同時存在O2—(2s22p6)離子和O—(2s22p5)離子,O2—離子的外電子殼層為滿殼層結構,O—離子的外電子殼層存在一個2p空穴.O2—離子的2p電子有一定的幾率以陽離子為媒介躍遷到鄰近O—離子外層軌道的2p空穴上,成為巡游電子,巡游電子在躍遷過程中自旋方向保持不變.

    2) 由于每個O2—離子的外層軌道存在自旋方向相反的兩個2p電子,造成一個O2—離子周圍的陽離子分成兩個磁性子晶格,在每個子晶格中巡游電子的自旋方向相同,但是兩個磁性子晶格中巡游電子的自旋方向相反.例如尖晶石結構鐵氧體中八面體位和四面體位的陽離子分別處于(A)和[B]子晶格.在每個子晶格中巡游電子的自旋方向相同,但是(A)和[B]子晶格中巡游電子的自旋方向相反.

    3) 在同一子晶格中,由于巡游電子在躍遷過程中自旋的方向保持不變,每個離子的電子自旋方向(包括局域電子和巡游電子)都必須遵守洪特(Hund)定則.因此,在同一子晶格中,如果兩個3d過渡金屬離子的3d電子數(shù)目(nd)同時滿足nd≥5(或它們都滿足nd≤4),則它們的磁矩平行排列;如果一個陽離子的nd≥5,另一個陽離子的nd≤4,則它們的磁矩反平行排列.

    應用IEO模型替代傳統(tǒng)的超交換和雙交換作用模型,研究了多系列尖晶石結構[30,34-45]和鈣鈦礦結構[8,9,31,46-49]磁性氧化物的磁有序問題,包括應用傳統(tǒng)模型難以給出合理解釋的Cr摻雜尖晶石鐵氧體的磁有序問題[34-39]、CrFe2O4的反常紅外光譜[40]、Ti摻雜導致尖晶石鐵氧體出現(xiàn)附加反鐵磁相[41,42]、鈣鈦礦結構錳氧化物La1—xSrxMnO3的磁矩隨Sr摻雜量變化的關系[46]等典型問題.

    3.2.2 IEM模型

    1) 基于γ射線衍射等研究結果,認為在3d過渡族原子(除Cu和Zn外)結合成金屬的過程中,由于受到原子間電子的泡利排斥力的擠壓作用,原子的大部分4s電子進入3d軌道,變成3d電子,剩余的4s電子作為自由電子.

    2) 處于費米能級附近的3d電子有一定幾率在鄰近原子的外層軌道間發(fā)生躍遷,形成巡游電子,其余的3d電子都是局域電子.

    3) 金屬的電阻率隨自由電子含量的增加而減小.自由電子的遷移過程受到晶格的弱周期性勢場影響,不受任何離子外層軌道的束縛,所以自由電子的自旋對材料磁矩沒有貢獻.在居里溫度以下巡游電子躍遷屬于自旋相關躍遷,當溫度接近居里溫度時,其躍遷幾率迅速減小;溫度越低,巡游電子躍遷對電導率的貢獻越大.

    應用IEM模型,解釋了Fe,Ni,Co,Cu金屬[32]以及NiCu合金[7]的電阻率與磁矩的關系.

    3.2.3 WEP模型

    1) 假設電子在一個離子外殼層中高速運動時自旋方向不變,由于相鄰離子的最外層軌道十分接近,其電子分別有一定幾率處于圖3(a),(b),(c)所示的狀態(tài).

    圖3 近鄰離子外層電子軌道的(a)外斯電子對和(b),(c)巡游電子示意圖[33]Fig.3.Illustrations of (a) a Weiss electron pair and (b) and(c) itinerant electrons in the outer orbits of adjacent ions[33].

    2) 當兩離子處于圖3(a)的狀態(tài)時,離子間的兩電子自旋磁矩反平行.由于最外層軌道不能同時容納自旋方向相同的兩個電子,所以電子不能在兩離子間交換,產(chǎn)生靜磁吸引能,同時也存在泡利排斥能,從而可處于吸引能和排斥能的短暫平衡態(tài),具有確定的平衡間距和壽命,稱為外斯電子對(WEP).本研究小組認為WEP的能量是磁性材料磁有序能的來源.

    3) 當兩個離子處于圖3(b)所示狀態(tài)時,兩個電子的自旋磁矩平行,容易發(fā)生互相交換,交換前后電子的自旋方向保持不變,并且這種狀態(tài)的兩個電子間存在磁性排斥能;當一個離子外層軌道有兩個電子,其相鄰的離子外層軌道只有一個電子,且處于圖3(c)的狀態(tài)時,中間的電子可以躍遷到右側離子上,并且保持自旋方向不變.圖3(b)和(c)所示的躍遷統(tǒng)稱為巡游電子的躍遷.

    上面分析解釋了在磁性材料中巡游電子的自旋方向保持不變的原因,如果自旋方向不同,就不能發(fā)生圖3(b)所示的交換或圖3(c)所示的躍遷.

    綜合應用IEO,IEM和WEP模型,可以解釋為什么Co,Fe,Ni,NiCu合金,Fe3O4和La1—xSrxMnO3具有不同的居里溫度[7,50].

    3.3 新模型中價電子結構與傳統(tǒng)模型的主要區(qū)別

    3.3.1 是否考慮氧化物中的負一價氧離子,導致不同的磁有序機制

    在傳統(tǒng)的超交換和雙交換作用模型中,假設磁性氧化物中的所有氧離子都是負二價離子,具有2s22p6的滿殼層結構,因而陽離子也都具有傳統(tǒng)的化合價數(shù)值.例如,在鈣鈦礦結構錳氧化物La1—xSrxMnO3中,當Sr含量為0時,所有Mn離子為Mn3+,認為樣品的反鐵磁性源于Mn3+與Mn3+離子間的超交換作用;當摻雜Sr時,出現(xiàn)Mn4+,并且Mn4+含量與Sr含量x相同,認為在Mn3+和Mn4+間存在雙交換作用,導致鐵磁耦合[5,12].

    基于前述電子譜的實驗結果,氧化物中同時存在負二價和負一價氧離子,其平均化合價的絕對值低于傳統(tǒng)觀點的數(shù)值,隨陰/陽離子含量比的增大而減小.例如,二氧化錳中氧離子平均化合價的絕對值小于一氧化錳中的數(shù)值.這是由于錳離子的第三、四電離能分別為33.67和51.2 eV,即使在二氧化錳中,錳離子也幾乎不會失去第四個電子[6].這直接導致利用IEO模型對于鈣鈦礦結構錳氧化物La1—xSrxMnO3磁有序問題的解釋不同于傳統(tǒng)模型[46]:當Sr含量x小于0.15時,樣品的反鐵磁性源于Mn2+與Mn3+離子間的反鐵磁耦合,樣品磁矩隨x的增加而增大,是由于Mn3+/Mn2+離子含量比隨 x的增加而增大;當 0.15<x<0.40時,所有Mn離子為Mn3+離子,樣品磁矩隨x的增加而減小,源于Mn3+離子間為傾角鐵磁耦合,并且傾角隨x的增加而增大.

    3.3.2 新舊模型的巡游電子定義不同

    在傳統(tǒng)理論中,用巡游電子模型解釋關于金屬磁性的直接交換作用.例如,對于磁性金屬Fe,Co,Ni,把3d和4s電子都看成可以在晶格中自由巡游的巡游電子,處于自旋向上和自旋向下的兩個能帶中[1].

    改進后的模型認為,在Fe,Co,Ni從自由原子形成磁性金屬的過程中,一部分4s電子進入3d軌道,變?yōu)?d電子,因而價電子分為自由電子、巡游電子和局域電子.自由電子不受離子實外層電子軌道的束縛,對樣品磁矩沒有貢獻;巡游電子指的是離子實的最外層軌道電子有一定幾率在最近鄰離子間躍遷,因而對樣品磁矩有貢獻;除巡游電子外的3d電子都是局域電子.由于價電子的軌道磁矩對樣品磁矩的貢獻很小[22,26],故可以把金屬和合金中平均每個原子磁矩的實驗值作為估算3d電子數(shù)目的依據(jù)[6,7,32].

    3.3.3 新舊模型中對鈣鈦礦錳氧化物導電機制的理解不同

    傳統(tǒng)的雙交換作用模型認為La1—xSrxMnO3的巡游電子源于Mn離子的3d電子,在居里溫度以下為金屬導電性,在居里溫度以上為半導體導電性,在居里溫度附近發(fā)生導電性的“金屬-半導體轉變”[12].

    根據(jù)IEO模型,對于鈣鈦礦結構錳氧化物La1—xSrxMnO3,在居里溫度以下巡游電子的電輸運行為屬于自旋相關電輸運,在居里溫度以上的半導體導電性屬于自旋無關電輸運[8],在居里溫度附近,巡游電子發(fā)生從自旋相關輸運到自旋無關輸運的轉變.

    3.3.4 新舊模型中對磁性金屬和合金導電機制的理解不同

    在傳統(tǒng)模型中,沒有考慮金屬中巡游電子和自由電子導電機制的區(qū)別.根據(jù)IEM模型,對于NiCu合金,在居里溫度以上,為自由電子導電,巡游電子的電導遠小于自由電子的電導;在居里溫度以下,是巡游電子與自由電子共同導電,溫度越低,巡游電子的電導越大[7].

    可見,在NiCu合金和鈣鈦礦結構錳氧化物La1—xSrxMnO3中,巡游電子的自旋相關躍遷對電輸運性質具有相似的影響[7,8]:在居里溫度以下,巡游電子的電輸運行為都屬于自旋相關輸運,在居里溫度附近,巡游電子發(fā)生從自旋相關輸運到自旋無關輸運的轉變.

    3.3.5 新舊模型對磁有序能來源的解釋不同

    對于外斯分子場能量的來源,傳統(tǒng)觀點認為完全是一種量子效應[1],在密度泛函理論中歸結在交換關聯(lián)能中.而交換關聯(lián)能不僅包括與磁有序相關的能量,還包括在泛函表達式中沒有明確給出的所有能量[51],并且迄今為止,還一直沒能給出交換關聯(lián)能的函數(shù)表達式,對于不同的材料,只能用不同模型進行擬合.

    根據(jù)外斯電子對模型,外斯分子場能量來源于在相鄰離子間形成的自旋方向相反的外斯電子對,這種處于激發(fā)態(tài)的電子對可在一定條件下形成,具有一定的形成幾率和一定的壽命[6,7,33,50].利用這個觀點,可成功解釋幾種典型的磁性材料為什么具有不同的居里溫度[6,7].

    3.4 自由原子電離能和電子親和能對晶體中價電子結構的影響

    對于上述改進的價電子結構模型的合理性,可利用自由原子的電離能和電子親和能進行分析.

    3.4.1 晶體中電子的束縛能與自由原子電離能的關系

    自由原子的電離能[52]是自由原子中電子離開所在能級到達(相對于離子實作用場)無窮遠處所需要的能量.晶體中離子對其電子的束縛能數(shù)據(jù)可以從XPS手冊[53]查到,這種束縛能是電子從所在能級到達費米能級所需的能量.圖4給出晶體中離子對1s電子的束縛能Eb和自由原子中1s電子的電離能VN隨原子序數(shù)的變化關系.可見,盡管Eb小于VN,但二者都隨原子序數(shù)的增大而迅速增大,也可理解為Eb隨VN的增大而增大.Eb與VN的這種關系是利用自由原子電離能分析材料中電子結構的依據(jù).

    圖4 晶體中離子對1s電子的束縛能Eb和自由原子中1s電子的電離能VN隨原子序數(shù)N的變化[52,53]Fig.4.Dependences on the atom number (N) of the binding energy (Eb) of 1 s electron in a crystal and the ionization energy (VN) of 1 s electron in an free atom[52,53].

    利用XPS還可得到晶體的價電子譜.Ley 等[54]給出了從Sc到Fe金屬的價帶譜,表明這些金屬的價電子分布在費米能級以下約12 eV的范圍內(nèi).研究CaO,ZnO,MnFe2O4,ZnFe2O4的價帶譜,這些氧化物的價電子也分布在費米能級以下約12 eV的范圍內(nèi)(圖5)[55].從圖5還可以看到一個有趣的結果,在ZnO和ZnFe2O4中,Zn 3d 電子集中分布在費米能級以下約8.2—11.5 eV范圍內(nèi),這是因為其中的Zn2+離子有10個 3d 電子,具有3d滿殼層結構,自由Zn原子的第三電離能為39.72 eV,極難失去第三個電子.這些價電子譜實驗結果說明絕大多數(shù)價電子是局域電子,只有費米能級附近很小范圍的電子形成巡游電子.這也說明上述關于巡游電子和局域電子的定義是合理的.

    圖5 CaO,ZnO,MnFe2O4,ZnFe2O4的價帶光電子譜[55]Fig.5.Valence band photoelectron spectra of samples CaO,ZnO,MnFe2O4 and ZnFe2O4[55].

    3.4.2 晶體中的價電子狀態(tài)與自由原子電離能和電子親和能的關系

    1) 鹵堿化合物是典型的離子晶體

    Li,Na,K,Rb,Cs的第一電離能在5.39—3.89 eV之間,F,Cl,Br,I的電子親和能在3.614—3.061 eV 之間.由于這些堿金屬的第一電離能與鹵素的電子親和能比較接近,當它們結合成鹵堿化合物時,鹵素原子容易從堿金屬得到一個電子,形成負一價離子,同時堿金屬原子形成正一價離子,所以鹵堿化合物是典型的離子晶體.

    2) 氧化物的平均化合價絕對值小于傳統(tǒng)觀點的數(shù)值

    首先考慮一氧化物,如果按照傳統(tǒng)觀點,在氧化物中每個氧離子都得到2個電子,形成負二價離子.注意到氧的第二電子親和能為8.08 eV,堿土離子Mg,Ca,Sr,Ba的第二電離能在15.04—10.00 eV 之間,3d過渡金屬離子的第二電離能在12.80—20.29 eV 之間.如果這些金屬離子與氧離子形成氧化物,氧離子要得到第二個電子,存在不同程度的困難,因而有一部分氧離子不能得到第二個電子,只能形成負一價氧離子.如前所述,陽離子的第二電離能越大,形成負一價氧離子的幾率越大,平均化合價的絕對值越小.

    當從Ti到Cu的3d過渡金屬與氧形成比例為3/4,2/3,1/2的氧化物時,由于氧離子的比例逐漸增大,按傳統(tǒng)觀點應出現(xiàn)三價和四價陽離子.但是,因為這些金屬離子的第三電離能在27.49—36.83 eV 之間,第四電離能在 43.27—55.2 eV 之間,難于形成四價陽離子,出現(xiàn)三價陽離子的比例也小于傳統(tǒng)觀點的值.所以氧化物的平均化合價絕對值隨其中氧離子含量比的增加而減小.例如,按照傳統(tǒng)觀點,在BaTiO3和SrTiO3中Ti離子應為正四價,而按照Cohen的計算[18,19]和本研究組的實驗結果[20,21],在BaTiO3和SrTiO3中絕大部分Ti離子為正三價,少量為正二價,沒有正四價的Ti離子.按照 Dupin等[17]的 XPS分析結果,在TiO2和ZrO2中,Ti和Zr離子的平均化合價分別為+2.3和+2.36,而不是+4.0.

    3) 金屬中的價電子絕大部分是局域電子

    在金屬中存在自由電子,這是因為金屬原子的第一電離能一般為幾個電子伏特,例如Fe,Co,Ni,Cu的第一電離能在7.87—7.73 eV 之間.對于作為良導體的金屬Cu,理論和實驗研究表明,其自由電子濃度為平均每個Cu原子貢獻一個自由電子.除自由電子外,其余的價電子受到原子核的束縛,形成離子實.這種離子實與離子晶體中的離子非常相似,只有最外層軌道的電子具有一定幾率在相鄰離子實間躍遷形成巡游電子(金屬Cu的離子實具有3d10價電子態(tài),不易形成巡游電子),其他電子都是局域電子,所以金屬中的絕大部分價電子是局域電子.

    4 相關研究面臨的挑戰(zhàn)

    4.1 密度泛函理論研究方法的改進

    利用密度泛函做材料模擬計算,對許多材料的預測取得了很好的效果,但是對于磁性材料的模擬和預測還存在很大困難,其主要原因在于把磁有序能與所有未知能量都包含在交換關聯(lián)能中,而交換關聯(lián)能至今還是未知函數(shù),對于不同材料交換關聯(lián)能的計算只能采用不同的方法進行模擬.如果能夠找到磁有序能函數(shù)的表達式,使其從交換關聯(lián)能中獨立出來,作為能量泛函表達式的一項,將推動磁性材料的模擬取得顯著進展.本研究組提出的WEP模型和對典型磁性材料居里溫度的解釋為此提供了一定的線索[6,7,50],但還需要深入研究.

    此外,目前利用密度泛函進行材料模擬,把每個原子的所有價電子都在程序中計算,計算量非常大,對于稍復雜的體系,計算量大到無法完成.然而,基于上述實驗結果和初步改進的價電子結構模型,大量的價電子都是局域電子,對于晶體的結合能并無影響.如果對軟件進行改進,把能夠確定的局域電子都當做內(nèi)層電子,只計算自由電子和巡游電子,可大幅度減少計算工作量.

    4.2 基于氧化物中存在負一價氧離子,對XPS分析方法進行系統(tǒng)研究

    目前關于氧化物的XPS分析,多數(shù)作者還沒有注意到負一價氧離子的存在,對于其中陽離子譜峰的分析還建立在所有氧離子都是負二價離子的基礎上.例如,對鈣鈦礦結構錳氧化物La1—xSrxMnO3中Mn離子的XPS主峰進行分峰,從中得到正三價和四價錳離子的含量比,而實際上其主峰并沒有雙峰的特點,這樣的分峰具有非常大的任意性,只有譜峰本身具有比較明顯的雙峰特點時才適于做分峰處理,例如O 1s峰,二、三價Ti離子峰[6,21].因此,應組織團隊,在這方面開展系統(tǒng)研究,建立權威的數(shù)據(jù)庫.

    4.3 多系列鈣鈦礦結構錳氧化物單晶薄膜的原位磁性、電輸運和XPS系統(tǒng)研究

    以La1—xSrxMnO3為代表的ABO3型鈣鈦礦結構錳氧化物是一類磁性和電輸運性質可多因素調節(jié)的材料[6,9,46-49]:調節(jié)A位的Sr含量,可以使Mn離子的磁矩從反鐵磁耦合變?yōu)殍F磁耦合[46];把A位的La換成Pr,可在A位引入弱的磁矩[47];以Cr替代部分Mn,可引起傾角鐵磁耦合[47];以Fe,Co,Ni替代部分Mn,可引起傾角反鐵磁耦合[9,48].對于這些替代,在樣品磁性發(fā)生變化的同時,電輸運性質也發(fā)生巨大變化[9,47,48].基于改進的價電子結構模型,對這類材料磁性和電輸運性質的精細而系統(tǒng)的分析,必將加深對于氧化物磁性和電輸運物理機制的理解.

    在這方面,雖然過去三十多年中有大量研究工作[56-58],但是,絕大部分工作都是基于其中氧離子全部為負二價離子,用雙交換和超交換作用模型分析其磁性和電輸運性質,對于存在的問題又用多種不同的模型進行修補,形成許多不同的觀點,實際上造成了對這類材料物理機制理解的混亂現(xiàn)象.

    本研究組在這方面做了多個系列樣品的研究工作,用O 2p巡游電子模型給出了系統(tǒng)的解釋[6,8,9,46-49].但都是建立在假設其中Mn離子和其他3d過渡金屬離子都不存在正四價離子的基礎上,沒有對其進行精細的XPS分析.這是因為所用的樣品都是粉末樣品,由于XPS探測深度較小,表面效應影響較大,對于其中的負一價和負二價氧離子含量比不能給出精細的分析結果.

    所以,組織力量對這類材料進行單晶薄膜的原位磁性、電輸運和XPS系統(tǒng)研究,澄清這類氧化物磁性和電輸運物理機制,將推動對于其他種類晶體結構氧化物磁耦合機制的理解.

    4.4 多系列金屬、合金磁性和電輸運性質的XPS系統(tǒng)研究

    利用提出的磁性金屬巡游電子模型(IEM模型)和磁有序能來源的外斯電子對模型(WEP模型),可擬合方俊鑫和陸棟所著《固體物理學》[59]中給出的NiCu合金在不同Cu含量下電阻率隨測試溫度變化的曲線,并能成功解釋Co,Fe,Ni和這些NiCu合金為什么具有不同的磁矩和居里溫度[7].迄今為止,還沒有發(fā)現(xiàn)其他研究組關于同時系統(tǒng)解釋金屬與合金磁性和電輸運性質的類似報道.

    如本文開頭所指出,由于材料的磁性和電輸運性質都是由其價電子結構所決定的,只有能夠同時解釋磁性和電輸運性質的磁有序模型,并且有足夠的實驗依據(jù),才能真正反映材料的實際價電子結構.

    然而,對于磁性金屬和合金,由于其電阻率非常小,一般的設備難于測量.可以查到一些早期的電阻率數(shù)據(jù),缺乏近年來的系統(tǒng)研究數(shù)據(jù),更難于找到同時給出磁性和電性的研究報告.所以,組織力量開展這方面的系統(tǒng)研究,非常必要.

    5 結 語

    總之,1950年代及其以前形成的價電子結構理論,雖然極大地推動了凝聚態(tài)物理學的發(fā)展,但是仍有大量的實驗結果不能得到合理的解釋.1970年代以來關于價電子結構的大量實驗結果為改進價電子結構理論提供了重要依據(jù),但是缺乏站在凝聚態(tài)物理的角度進行系統(tǒng)地梳理和研究.在居里溫度以下,磁性材料存在磁有序能和自旋相關輸運,可以對價電子結構研究提供更多的信息,所以,基于大量電子譜實驗結果,從鐵磁學基礎研究的角度,對典型磁性材料的價電子結構開展系統(tǒng)而深入的研究,是本世紀鐵磁學研究乃至凝聚態(tài)物理研究面臨的重大機遇與挑戰(zhàn).

    猜你喜歡
    電離能巡游氧化物
    學科融合視野下的探究教學實踐*
    ——以物質結構與性質模塊“元素周期律”教學為例
    化學教與學(2023年3期)2023-02-09 08:33:22
    藝術巡游
    家居廊(2022年12期)2023-01-05 01:52:11
    “龍馬”巡游
    相轉化法在固體氧化物燃料電池中的應用
    陶瓷學報(2020年6期)2021-01-26 00:37:56
    細說『碳和碳的氧化物』
    氧化物的分類及其中的“不一定”
    淺談A Level 化學中電離能的影響因素及變化趨勢與普高教學的差異
    乾隆巡游
    寶藏(2018年1期)2018-01-31 02:05:09
    看巡游踩街
    小主人報(2016年9期)2016-12-01 06:23:23
    基于圖表分析的課堂知識建構與能力提升
    化學教與學(2014年6期)2014-07-03 10:04:09
    一个人看的www免费观看视频| 亚洲精品久久国产高清桃花| 国产精品1区2区在线观看.| 啦啦啦观看免费观看视频高清| 久久伊人香网站| 蜜桃亚洲精品一区二区三区| 能在线免费观看的黄片| av视频在线观看入口| 亚洲av电影在线进入| 亚洲自偷自拍三级| 久久亚洲精品不卡| 亚洲 国产 在线| 成人特级黄色片久久久久久久| 国产午夜精品论理片| 赤兔流量卡办理| 国产大屁股一区二区在线视频| 欧美中文日本在线观看视频| 久99久视频精品免费| 国产精品久久久久久久久免 | 美女xxoo啪啪120秒动态图 | 极品教师在线视频| 在线免费观看不下载黄p国产 | 日本一二三区视频观看| 久久天躁狠狠躁夜夜2o2o| 免费搜索国产男女视频| 制服丝袜大香蕉在线| 九九热线精品视视频播放| 国产黄a三级三级三级人| 99在线人妻在线中文字幕| 在线免费观看的www视频| 久久精品久久久久久噜噜老黄 | 成人亚洲精品av一区二区| 18禁黄网站禁片午夜丰满| 最新在线观看一区二区三区| 精品免费久久久久久久清纯| 波多野结衣高清作品| 国产私拍福利视频在线观看| 97人妻精品一区二区三区麻豆| 亚洲精品在线观看二区| 国产极品精品免费视频能看的| 久久香蕉精品热| 免费无遮挡裸体视频| 亚洲,欧美,日韩| 亚洲成av人片在线播放无| 看黄色毛片网站| 精品一区二区三区av网在线观看| 蜜桃久久精品国产亚洲av| 十八禁人妻一区二区| 久久久久九九精品影院| 成人无遮挡网站| 成人无遮挡网站| 51国产日韩欧美| 麻豆成人午夜福利视频| 青草久久国产| 18禁裸乳无遮挡免费网站照片| 欧美黑人巨大hd| 一个人观看的视频www高清免费观看| 国产伦人伦偷精品视频| 一个人免费在线观看电影| 制服丝袜大香蕉在线| 久久6这里有精品| www.www免费av| 精品国产三级普通话版| 男人舔奶头视频| 大型黄色视频在线免费观看| ponron亚洲| 91午夜精品亚洲一区二区三区 | 日本黄色视频三级网站网址| 亚洲 欧美 日韩 在线 免费| 国产午夜福利久久久久久| 黄色日韩在线| 成人国产综合亚洲| 国产午夜精品论理片| 色5月婷婷丁香| 日韩精品中文字幕看吧| 国产午夜精品论理片| 亚洲午夜理论影院| 中文资源天堂在线| 午夜日韩欧美国产| 日本一本二区三区精品| 九色国产91popny在线| 99久久成人亚洲精品观看| 尤物成人国产欧美一区二区三区| 99精品久久久久人妻精品| 中国美女看黄片| 性色avwww在线观看| 国产欧美日韩一区二区精品| 在现免费观看毛片| 青草久久国产| 18禁在线播放成人免费| 日本精品一区二区三区蜜桃| 国产免费av片在线观看野外av| 国产蜜桃级精品一区二区三区| 久久亚洲精品不卡| 久久精品国产清高在天天线| 亚洲av美国av| 久久精品影院6| 亚洲av日韩精品久久久久久密| 天美传媒精品一区二区| 天美传媒精品一区二区| 欧美在线一区亚洲| 成人午夜高清在线视频| 日韩av在线大香蕉| 免费搜索国产男女视频| 一边摸一边抽搐一进一小说| 最后的刺客免费高清国语| 国产精品嫩草影院av在线观看 | 色综合亚洲欧美另类图片| 热99在线观看视频| www.999成人在线观看| 精品午夜福利在线看| 变态另类丝袜制服| 久久精品影院6| АⅤ资源中文在线天堂| 国产欧美日韩精品亚洲av| 九九热线精品视视频播放| 深爱激情五月婷婷| 成人午夜高清在线视频| 色播亚洲综合网| 色播亚洲综合网| 最后的刺客免费高清国语| 男插女下体视频免费在线播放| 午夜福利视频1000在线观看| 男人和女人高潮做爰伦理| 真人一进一出gif抽搐免费| 日本一二三区视频观看| 国产真实乱freesex| a在线观看视频网站| 淫秽高清视频在线观看| 久久久精品欧美日韩精品| 乱人视频在线观看| 亚洲av第一区精品v没综合| 国产高清有码在线观看视频| 在线观看午夜福利视频| 国产真实乱freesex| 国产成人啪精品午夜网站| 国产 一区 欧美 日韩| 免费观看精品视频网站| 久久久久免费精品人妻一区二区| 黄色配什么色好看| 熟女电影av网| 日韩免费av在线播放| 亚洲欧美日韩无卡精品| 国产精品av视频在线免费观看| 日本精品一区二区三区蜜桃| a级毛片a级免费在线| 少妇裸体淫交视频免费看高清| 欧美成人性av电影在线观看| 夜夜爽天天搞| 国产精品嫩草影院av在线观看 | 久久中文看片网| 国产精品三级大全| 国产精品女同一区二区软件 | 91在线精品国自产拍蜜月| 精品熟女少妇八av免费久了| 51国产日韩欧美| 久久久国产成人精品二区| 中文字幕免费在线视频6| 中文字幕久久专区| 国产伦一二天堂av在线观看| 波野结衣二区三区在线| 亚洲电影在线观看av| 成人毛片a级毛片在线播放| 日本五十路高清| 日韩欧美国产一区二区入口| 亚洲成人精品中文字幕电影| 色综合亚洲欧美另类图片| 欧美日韩中文字幕国产精品一区二区三区| 少妇裸体淫交视频免费看高清| 国模一区二区三区四区视频| 麻豆成人午夜福利视频| 蜜桃亚洲精品一区二区三区| 高清在线国产一区| 欧美潮喷喷水| 国产伦一二天堂av在线观看| 日韩精品中文字幕看吧| 午夜福利在线观看免费完整高清在 | 69人妻影院| 国产亚洲精品久久久com| 日韩有码中文字幕| 一级av片app| 亚洲avbb在线观看| 人妻丰满熟妇av一区二区三区| 国产免费男女视频| 亚洲国产高清在线一区二区三| 又紧又爽又黄一区二区| 国产视频一区二区在线看| 国产精品人妻久久久久久| 欧美日韩国产亚洲二区| 欧美+亚洲+日韩+国产| 亚洲熟妇熟女久久| 一区二区三区免费毛片| 亚洲不卡免费看| 99久久精品一区二区三区| 人妻制服诱惑在线中文字幕| 亚洲第一区二区三区不卡| 精品人妻1区二区| 天堂动漫精品| 性色av乱码一区二区三区2| 美女cb高潮喷水在线观看| 极品教师在线视频| 免费人成在线观看视频色| 黄色配什么色好看| 悠悠久久av| 久久精品夜夜夜夜夜久久蜜豆| 久久久久久大精品| a在线观看视频网站| 老司机午夜福利在线观看视频| 窝窝影院91人妻| 在线观看免费视频日本深夜| 亚洲精品一区av在线观看| 国产精品,欧美在线| 亚洲狠狠婷婷综合久久图片| a在线观看视频网站| 极品教师在线视频| 99热精品在线国产| 欧美激情久久久久久爽电影| 悠悠久久av| 亚洲片人在线观看| 亚洲狠狠婷婷综合久久图片| 欧美成人一区二区免费高清观看| 狂野欧美白嫩少妇大欣赏| 天堂影院成人在线观看| 欧美激情国产日韩精品一区| 欧美日韩黄片免| 亚洲精品成人久久久久久| 亚洲熟妇中文字幕五十中出| 一级毛片久久久久久久久女| 成年人黄色毛片网站| 午夜福利成人在线免费观看| 欧美最黄视频在线播放免费| 国产探花极品一区二区| 熟女电影av网| 久久久久九九精品影院| 99在线视频只有这里精品首页| 国产精品久久电影中文字幕| 人人妻人人看人人澡| 日本黄色视频三级网站网址| 五月伊人婷婷丁香| 亚洲国产欧美人成| a在线观看视频网站| 国产av一区在线观看免费| 日本黄大片高清| www.色视频.com| 久久精品综合一区二区三区| 精品久久久久久久末码| 两人在一起打扑克的视频| 精品国内亚洲2022精品成人| 国语自产精品视频在线第100页| 日韩欧美国产一区二区入口| 天堂√8在线中文| 亚洲七黄色美女视频| 亚洲成a人片在线一区二区| 怎么达到女性高潮| 亚洲专区国产一区二区| 我的女老师完整版在线观看| 午夜免费激情av| 怎么达到女性高潮| 久久精品91蜜桃| 久久亚洲精品不卡| 日日干狠狠操夜夜爽| 欧美3d第一页| 9191精品国产免费久久| 91字幕亚洲| 变态另类丝袜制服| 欧美乱妇无乱码| 十八禁国产超污无遮挡网站| 日本成人三级电影网站| 久久亚洲精品不卡| 69人妻影院| 免费看日本二区| 日日摸夜夜添夜夜添av毛片 | 三级男女做爰猛烈吃奶摸视频| 亚洲精品日韩av片在线观看| av在线天堂中文字幕| 国产成人aa在线观看| 夜夜躁狠狠躁天天躁| 亚洲精品在线美女| 久久99热这里只有精品18| 欧美成人性av电影在线观看| 蜜桃亚洲精品一区二区三区| 亚洲无线观看免费| 久久国产乱子伦精品免费另类| 99国产综合亚洲精品| 精品无人区乱码1区二区| 看免费av毛片| 成年女人看的毛片在线观看| 国产精品av视频在线免费观看| 岛国在线免费视频观看| 蜜桃亚洲精品一区二区三区| 成人永久免费在线观看视频| 国产野战对白在线观看| 亚洲无线在线观看| 熟女人妻精品中文字幕| 国产av麻豆久久久久久久| 国产乱人视频| 嫩草影院精品99| 很黄的视频免费| 午夜免费男女啪啪视频观看 | 熟女人妻精品中文字幕| 国产探花在线观看一区二区| 精品人妻偷拍中文字幕| 国产精品三级大全| 国产麻豆成人av免费视频| 麻豆一二三区av精品| 99久久久亚洲精品蜜臀av| 亚洲精品粉嫩美女一区| 午夜精品久久久久久毛片777| 真实男女啪啪啪动态图| 亚洲中文字幕日韩| 精品久久久久久久久av| 亚洲人与动物交配视频| 丁香六月欧美| 亚洲av成人不卡在线观看播放网| 亚洲最大成人av| 欧美黄色淫秽网站| 在线十欧美十亚洲十日本专区| 搡女人真爽免费视频火全软件 | 免费观看精品视频网站| 久久久久久九九精品二区国产| 亚洲无线在线观看| 欧美3d第一页| 亚洲中文字幕日韩| 噜噜噜噜噜久久久久久91| 亚洲美女搞黄在线观看 | 97超级碰碰碰精品色视频在线观看| 好看av亚洲va欧美ⅴa在| 亚洲成人中文字幕在线播放| 中文字幕熟女人妻在线| 亚洲成a人片在线一区二区| 国产精品亚洲美女久久久| 美女免费视频网站| 国产高清有码在线观看视频| 亚洲avbb在线观看| 在线观看免费视频日本深夜| 九色成人免费人妻av| 老司机午夜十八禁免费视频| 特级一级黄色大片| 久久精品夜夜夜夜夜久久蜜豆| 亚洲精品久久国产高清桃花| 国产高清三级在线| 91久久精品电影网| 欧美日韩中文字幕国产精品一区二区三区| 97超级碰碰碰精品色视频在线观看| 国产熟女xx| 村上凉子中文字幕在线| 国产黄色小视频在线观看| 国产精品久久久久久久久免 | 亚洲真实伦在线观看| 日本a在线网址| 亚洲avbb在线观看| 小说图片视频综合网站| 亚洲国产欧洲综合997久久,| 亚洲欧美激情综合另类| 在线播放无遮挡| 午夜福利欧美成人| 午夜福利视频1000在线观看| 午夜影院日韩av| 一级毛片久久久久久久久女| 欧美高清性xxxxhd video| 一区二区三区激情视频| 99在线人妻在线中文字幕| 午夜激情欧美在线| 日本熟妇午夜| 国产私拍福利视频在线观看| 成年人黄色毛片网站| 国产熟女xx| 亚洲自偷自拍三级| 国产视频内射| 亚洲av五月六月丁香网| 欧美黄色片欧美黄色片| 51国产日韩欧美| 伦理电影大哥的女人| 亚洲国产欧洲综合997久久,| 免费看日本二区| 亚洲成人久久性| 日本免费a在线| 乱码一卡2卡4卡精品| 成人毛片a级毛片在线播放| 欧美日韩中文字幕国产精品一区二区三区| 97超级碰碰碰精品色视频在线观看| 国产探花极品一区二区| 高清日韩中文字幕在线| 国产真实乱freesex| 亚洲一区二区三区不卡视频| 热99在线观看视频| 亚洲欧美日韩高清专用| 好看av亚洲va欧美ⅴa在| 日本 av在线| 在现免费观看毛片| 中文字幕高清在线视频| 无遮挡黄片免费观看| 一个人观看的视频www高清免费观看| 又黄又爽又刺激的免费视频.| 大型黄色视频在线免费观看| 欧美日韩乱码在线| 一本精品99久久精品77| 国内少妇人妻偷人精品xxx网站| 免费看美女性在线毛片视频| 丰满人妻熟妇乱又伦精品不卡| 身体一侧抽搐| 国内毛片毛片毛片毛片毛片| 久久伊人香网站| 免费在线观看影片大全网站| 美女xxoo啪啪120秒动态图 | av国产免费在线观看| 日本一二三区视频观看| 国产熟女xx| 亚洲综合色惰| 亚洲人与动物交配视频| 老熟妇乱子伦视频在线观看| 天堂动漫精品| 婷婷丁香在线五月| 国产一级毛片七仙女欲春2| 国产精品精品国产色婷婷| 两人在一起打扑克的视频| 一个人看的www免费观看视频| 少妇的逼好多水| 亚洲av电影不卡..在线观看| 亚洲电影在线观看av| 黄色视频,在线免费观看| 99热这里只有精品一区| 国产午夜福利久久久久久| 麻豆国产97在线/欧美| 91麻豆av在线| 亚洲第一区二区三区不卡| 国产黄色小视频在线观看| 亚洲三级黄色毛片| 久久国产乱子伦精品免费另类| 毛片女人毛片| 国产免费男女视频| 长腿黑丝高跟| 中文字幕免费在线视频6| 好男人电影高清在线观看| 亚洲av熟女| 99久久精品一区二区三区| 婷婷色综合大香蕉| 日本一二三区视频观看| 搡老熟女国产l中国老女人| 少妇被粗大猛烈的视频| 99在线视频只有这里精品首页| 伊人久久精品亚洲午夜| 欧美黑人欧美精品刺激| x7x7x7水蜜桃| 一区二区三区激情视频| 欧美一区二区国产精品久久精品| 亚洲av免费高清在线观看| 亚洲成av人片在线播放无| 久久人人精品亚洲av| 90打野战视频偷拍视频| 丰满人妻熟妇乱又伦精品不卡| 久久精品人妻少妇| 久久午夜福利片| 国产精品不卡视频一区二区 | 午夜福利成人在线免费观看| 国产一级毛片七仙女欲春2| 搡老妇女老女人老熟妇| 中文字幕久久专区| 少妇丰满av| 狂野欧美白嫩少妇大欣赏| av福利片在线观看| 国产精品精品国产色婷婷| 久久午夜亚洲精品久久| 日本黄色片子视频| 国产一区二区三区视频了| 亚洲七黄色美女视频| av专区在线播放| 香蕉av资源在线| 国产精品伦人一区二区| 1000部很黄的大片| 亚洲五月天丁香| 欧美日韩国产亚洲二区| 欧美精品国产亚洲| www.熟女人妻精品国产| 亚洲成人精品中文字幕电影| 可以在线观看的亚洲视频| 亚洲人成网站在线播放欧美日韩| 一夜夜www| 精品国内亚洲2022精品成人| 搡女人真爽免费视频火全软件 | 国产激情偷乱视频一区二区| 亚洲美女视频黄频| 亚洲欧美精品综合久久99| 亚洲精品在线观看二区| 麻豆国产av国片精品| 久久香蕉精品热| 一个人观看的视频www高清免费观看| 91狼人影院| 国产三级黄色录像| 久久久久久久亚洲中文字幕 | 亚洲av五月六月丁香网| 有码 亚洲区| 久久人妻av系列| 中文字幕人成人乱码亚洲影| 很黄的视频免费| 欧美日韩黄片免| 毛片一级片免费看久久久久 | 一卡2卡三卡四卡精品乱码亚洲| 婷婷色综合大香蕉| 亚洲一区高清亚洲精品| 精品人妻1区二区| 国产亚洲欧美98| 嫩草影视91久久| 色视频www国产| 日韩成人在线观看一区二区三区| 国产成人av教育| 午夜视频国产福利| 嫩草影院精品99| 一本综合久久免费| 美女xxoo啪啪120秒动态图 | 国产黄a三级三级三级人| 十八禁国产超污无遮挡网站| 88av欧美| 国产精品野战在线观看| 露出奶头的视频| 国产成+人综合+亚洲专区| 国产欧美日韩精品亚洲av| 国产成人a区在线观看| 国产国拍精品亚洲av在线观看| 内地一区二区视频在线| 宅男免费午夜| 精品99又大又爽又粗少妇毛片 | 3wmmmm亚洲av在线观看| 亚洲欧美激情综合另类| 91午夜精品亚洲一区二区三区 | 亚洲成人久久爱视频| 嫩草影视91久久| 在线看三级毛片| 亚洲自拍偷在线| 人妻久久中文字幕网| 亚洲美女黄片视频| 国产精品自产拍在线观看55亚洲| av专区在线播放| 国产一区二区三区在线臀色熟女| 小蜜桃在线观看免费完整版高清| 久久久精品欧美日韩精品| 看片在线看免费视频| 色哟哟·www| 成人一区二区视频在线观看| 亚洲国产欧洲综合997久久,| 精品久久久久久久久久免费视频| 他把我摸到了高潮在线观看| 亚洲国产精品久久男人天堂| 五月伊人婷婷丁香| 亚洲欧美精品综合久久99| 国内揄拍国产精品人妻在线| 99热这里只有精品一区| 久久精品夜夜夜夜夜久久蜜豆| 免费黄网站久久成人精品 | 成年人黄色毛片网站| 每晚都被弄得嗷嗷叫到高潮| 1024手机看黄色片| 两个人的视频大全免费| 国产人妻一区二区三区在| 亚洲精品亚洲一区二区| 午夜激情福利司机影院| 又粗又爽又猛毛片免费看| 99热这里只有是精品在线观看 | 亚洲美女搞黄在线观看 | 三级男女做爰猛烈吃奶摸视频| 十八禁人妻一区二区| 在线观看免费视频日本深夜| 1024手机看黄色片| 男人舔奶头视频| 精品欧美国产一区二区三| 久久热精品热| 特大巨黑吊av在线直播| 国产高清激情床上av| 亚洲在线观看片| a级毛片a级免费在线| 国产欧美日韩精品一区二区| 偷拍熟女少妇极品色| 99久久九九国产精品国产免费| 久久久久国内视频| 身体一侧抽搐| 久99久视频精品免费| 国产真实伦视频高清在线观看 | 色尼玛亚洲综合影院| 久久中文看片网| 成人欧美大片| 久久久久久久午夜电影| 久久久久久久久中文| 深爱激情五月婷婷| 老司机深夜福利视频在线观看| 亚洲av熟女| 亚洲电影在线观看av| 日韩有码中文字幕| 亚洲自拍偷在线| 亚洲av不卡在线观看| 精品一区二区免费观看| 欧美一级a爱片免费观看看| 综合色av麻豆| 欧美精品啪啪一区二区三区| 毛片一级片免费看久久久久 | 午夜激情福利司机影院| 97碰自拍视频| 欧美色视频一区免费| 欧美日韩综合久久久久久 | 亚洲国产欧美人成| 色5月婷婷丁香| 亚洲aⅴ乱码一区二区在线播放| 99久久九九国产精品国产免费| 18美女黄网站色大片免费观看| 两个人的视频大全免费| 精品无人区乱码1区二区| 亚洲成人中文字幕在线播放| 国产久久久一区二区三区| 亚洲av美国av| av专区在线播放| 深爱激情五月婷婷| 欧美不卡视频在线免费观看| 亚洲五月婷婷丁香| 久久欧美精品欧美久久欧美| 亚洲avbb在线观看| 天堂动漫精品|