自旋-軌道耦合對(duì)雙層鈣鈦礦Ba2CoOsO6電磁性質(zhì)的影響

摘要：采用第一性原理計(jì)算方法，研究了雙層鈣鈦礦氧化物Ba2CoOsO6的晶體結(jié)構(gòu)、電子結(jié)構(gòu)及自旋-軌道耦合作用（SOC）對(duì)其電磁性質(zhì)的影響，并分析電子結(jié)構(gòu)和自旋-軌道耦合作用對(duì)其宏觀物理性質(zhì)的影響.結(jié)果表明：Ba2CoOsO6為反鐵磁窄帶隙半導(dǎo)體，存在Os1/Co1反位缺陷；電子結(jié)構(gòu)的分析結(jié)果證實(shí)了Ba2CoOsO6的電荷分布為Ba2+2Co2+Os6+O2-6，確定了2價(jià)Co離子和6價(jià)Os離子的存在；SOC的存在減小了Co離子和Os離子的自旋磁矩，同時(shí)減小了帶隙；根據(jù)第一性原理計(jì)算結(jié)果可知，自旋-軌道耦合作用對(duì)電子結(jié)構(gòu)和電磁性質(zhì)的影響不可忽視，相較于其他方法，該方法計(jì)算值與試驗(yàn)分析結(jié)果更接近，也驗(yàn)證了理論計(jì)算的準(zhǔn)確性.

關(guān)鍵詞： "Ba2CoOsO6； 雙層鈣鈦礦； 自旋-軌道耦合； 第一性原理計(jì)算； 反鐵磁耦合； 半導(dǎo)體

中圖分類號(hào)： TM 23文獻(xiàn)標(biāo)志碼： "A文章編號(hào)： ""1671-7775（2024）06-0739-06

引文格式： "李洪蘋，孫安，張瑤明，等. 自旋-軌道耦合對(duì)雙層鈣鈦礦Ba2CoOsO6電磁性質(zhì)的影響[J].江蘇大學(xué)學(xué)報(bào)（自然科學(xué)版），2024，45（6）：739-744.

Effect of spinorbit coupling on magnetic and electronic properties of

double perovskite Ba2CoOsO6

LI Hongping， SUN An， ZHANG Yaoming， GUO Baochang

（School of Materials Science and Engineering， Jiangsu University， Zhenjiang， Jiangsu 212013， China）

Abstract： The crystal structure， electronic structure and the effect of spinorbit coupling （SOC） on the magnetic and electronic properties of double perovskite Ba2CoOsO6 were systematically investigated by firstprinciples calculations. Through the coupling analysis of electron structure and spin orbit， the influence on the macroscopic physical properties was discussed. The results show that Ba2CoOsO6 is antiferromagnetic semiconductor material with narrow band gap and disordered Os1/Co1 defects. The electronic structure analysis confirms that the charge combination of Ba2CoOsO6 is Ba2+2Co2+Os6+O2-6， which verifies the existence of Co2+ and Os6+. The spin magnetic moment of Co and Os ions is decreased by SOC， and the band gap is narrowed. The results of firstprinciples calculations show that the effect of spinorbit coupling on electronic structure and electromagnetic properties can not be ignored. Compared with other methods， the calculated values are closer to the results of experimental analysis， which can verify the theoretical calculation.

Key words： "Ba2CoOsO6; double perovskite; spinorbit coupling; firstprinciples calculation;antiferromagnetic coupling; semiconductor

自從有關(guān)雙層鈣鈦礦型A2BB′O6的研究被報(bào)道以來，雙層鈣鈦礦在自旋電子材料領(lǐng)域的應(yīng)用便受到人們的關(guān)注.國內(nèi)外學(xué)者已合成了更多新型的鈣鈦礦化合物.雙層鈣鈦礦可以被看作簡單鈣鈦礦ABO3中50%的B位離子被另一種元素取代而形成的新型鈣鈦礦結(jié)構(gòu).在A2BB′O6晶體結(jié)構(gòu)中，A位通常是Ba、Sr、Ca等堿土和鑭系金屬元素，B/B′位通常被Co、Fe、Cr等過渡金屬元素占據(jù).與簡單鈣鈦礦ABO3相比，雙層鈣鈦礦可以在B-B、B′-B′和B-B′之間產(chǎn)生多種競爭性磁、電相互作用.A、B、B′位離子的不同組合使得雙層鈣鈦礦具有半金屬[1]、巨磁阻[2]和高居里溫度[3]等性質(zhì)，因此雙層鈣鈦礦具有廣闊的應(yīng)用前景.

含5d電子軌道和3d電子軌道元素的雙層鈣鈦礦材料中，電子自旋和軌道運(yùn)動(dòng)之間的相互作用十分復(fù)雜.FENG H. L.等[4]在高壓和高溫下合成具有3d軌道和5d軌道過渡金屬的雙層鈣鈦礦氧化物Ca2FeOsO6，其空間群為P21/n；Ca2FeOsO6在320 K左右的溫度下出現(xiàn)鐵磁轉(zhuǎn)變，Ca2FeOsO6的電學(xué)性質(zhì)接近半金屬狀態(tài).YAN B.H.等[5]發(fā)現(xiàn)Sr2CoOsO6從室溫冷卻到2 K的過程中，在相變溫度為108 K時(shí)，Co（3d7）和Os（5d2）的相互作用力矩是動(dòng)態(tài)的，而它們的平均作用力矩是長程有序的；在相變溫度低于67 K時(shí)，Co首先形成非共線自旋傾斜反鐵磁狀態(tài)，而動(dòng)態(tài)波動(dòng)的Os磁矩（magnetic moment）在相變溫度約為5 K時(shí)為隨機(jī)傾斜反鐵磁狀態(tài)[5].近些年，研究人員開發(fā)了含有4d和5d軌道過渡金屬的新型鈣鈦礦材料.

由于材料合成試驗(yàn)條件嚴(yán)苛，很難得到完全有序的A2BB′O6型鈣鈦礦，A位和B位離子的排列順序容易受到離子半徑、合成條件和電子構(gòu)型等因素影響而變成無序排列.在A2BB′O6型雙層鈣鈦礦中，B位和B′位陽離子更容易發(fā)生反位占據(jù)，形成反位缺陷.例如，ZHAO Q.等[6]通過高溫高壓合成方法制備了Pb2FeOsO6，結(jié)構(gòu)細(xì)化表明，F(xiàn)e和Os之間的反位占據(jù)率約為16.6%，相鄰的Fe3+和Os5+之間產(chǎn)生反鐵磁相互作用，形成了一個(gè)長程有序的鐵磁躍遷.GE Z. Z.等[7]通過試驗(yàn)研究發(fā)現(xiàn)，F(xiàn)e的半滿3d軌道和Os的t2g軌道之間存在超交換作用，理論證明Pb2FeOsO6為半導(dǎo)體.FENG H. L.等[8]通過試驗(yàn)研究發(fā)現(xiàn)，雙層鈣鈦礦氧化物Ba2Fe1.12Os0.88O6具有多種磁電特性；與Ca2FeOsO6和Sr2FeOsO6相比，該化合物具有罕見的6L型鈣鈦礦結(jié)構(gòu)，在高于室溫條件下表現(xiàn)出良好的鐵磁有序，其特殊磁性是由B位和B′位離子無序形成的不均勻磁態(tài)所致[8].

采用常壓下固相法合成的雙層鈣鈦礦氧化物Ba2CoOsO6是一種6L型鈣鈦礦結(jié)構(gòu)氧化物，其B位和B′位形成無序排列，同時(shí)具有3d電子軌道和5d電子軌道的過渡金屬元素.探討B(tài)位無序排列的鈣鈦礦氧化物Ba2CoOsO6的電磁性質(zhì)以及自旋-軌道耦合對(duì)其電磁性質(zhì)的影響，有助于將其應(yīng)用于更多的領(lǐng)域.對(duì)鈣鈦礦電磁性質(zhì)進(jìn)行探討，揭示其背后的相互耦合作用機(jī)制引發(fā)的一些特殊效應(yīng)是目前鈣鈦礦的重要研究方向.

為此，筆者采用第一性原理計(jì)算方法，研究雙層鈣鈦礦氧化物Ba2CoOsO6的晶體結(jié)構(gòu)、電子結(jié)構(gòu)及自旋-軌道耦合（SOC）對(duì)其電磁性質(zhì)的影響，并分析電子結(jié)構(gòu)和自旋-軌道耦合對(duì)其宏觀物理性質(zhì)的影響.

1計(jì)算方法

筆者采用了第一性原理計(jì)算法，利用從頭算模擬軟件包（VASP） [9]來優(yōu)化雙層鈣鈦礦氧化物Ba2CoOsO6的晶體結(jié)構(gòu)，并分析其電磁性質(zhì).計(jì)算過程中，采用廣義梯度近似（generalized gradient approximation， GGA）下的PBE交換關(guān)聯(lián)能（GGAPBE）來準(zhǔn)確地描述該化合物的物理化學(xué)性質(zhì).平面波展開的截段能量設(shè)置為400 eV.根據(jù)擴(kuò)胞后晶格參數(shù)，采用基于MonkhorstPack的3×7×1的K點(diǎn)模擬布里淵區(qū).系統(tǒng)總能量收斂標(biāo)準(zhǔn)為變化量小于1.0×10-5 eV/atom，原子間相互作用力小于0.01 eV/時(shí)晶體結(jié)構(gòu)完全優(yōu)化.Ba、 Co、Os和 O 的價(jià)電子構(gòu)型分別為6s2、3d74s2、5d66s2和 2s22p4.在此鈣鈦礦體系中，Co離子和Os離子分別具有3d電子軌道和5d電子軌道，存在電子間強(qiáng)關(guān)聯(lián).筆者采用GGA+U方法，利用有效Hubbard參數(shù)Ueff（Ueff=U-J，U為Hubbard參數(shù)，J為交換參數(shù)），研究電子強(qiáng)相關(guān)效應(yīng)對(duì)體系電磁性質(zhì)的影響.分別用UCo和UOs表示Co和Os的有效Hubbard參數(shù).此外，由于Co具有多重自旋態(tài)，以及含5d電子軌道的過渡金屬元素Os的存在，因而理論計(jì)算中考慮了SOC作用，用以準(zhǔn)確描述自旋態(tài)和電荷態(tài).

2結(jié)果和討論

2.1室溫下的晶體結(jié)構(gòu)分析

雙層鈣鈦礦氧化物Ba2CoOsO6的結(jié)構(gòu)如圖1所示，圖中藍(lán)色、灰色、紅色和綠色分別為Co、Os、O和Ba原子.其空間群為P3-m1 （No.164）[10].

中子衍射試驗(yàn)[10]表明，Co3、Co4和Os2排列有序，而Os1和Co1處于排列紊亂狀態(tài).Os1和Co1占比分別為0.88和0.12.通過雙層鈣鈦礦Ba2CoOsO6的單晶X射線衍射（SXRD）和中子粉末衍射（NPD）試驗(yàn)[10]得到晶體結(jié)構(gòu)參數(shù)試驗(yàn)值.采用GGA和GGA+U方法對(duì)雙層鈣鈦礦氧化物Ba2CoOsO6晶體結(jié)構(gòu)進(jìn)行優(yōu)化.

為了進(jìn)一步分析雙層鈣鈦礦Ba2CoOsO6晶體結(jié)構(gòu)的穩(wěn)定性，考慮Os1和Co1的反位缺陷，再次進(jìn)行基于GGA+U方法的晶格優(yōu)化.

基于3種優(yōu)化方法的雙層鈣鈦礦Ba2CoOsO6晶格常數(shù)、鍵角和鍵長計(jì)算值與試驗(yàn)值對(duì)比結(jié)果分別如表1-3[10]所示，其中L為鍵長，θ為鍵角.由表1-3可知，GGA計(jì)算得到的晶格常數(shù)（a=5.757 ，b=14.088 ）略小于GGA+U的計(jì)算值（a=5.795 ，b=14.212 ）.在合理誤差范圍內(nèi)，兩種方法計(jì)算得到的鍵角均與試驗(yàn)值接近，但采用GGA+U方法計(jì)算得到的鈣鈦礦氧化物Ba2CoOsO6晶體結(jié)構(gòu)參數(shù)更接近表格中試驗(yàn)值，可以滿足后續(xù)電磁性質(zhì)研究的要求.

由表1可知：相較于完全有序狀態(tài)下基于GGA+U計(jì)算得到的晶格常數(shù)（a=5.795 ，b=14.212 ），考慮反位缺陷計(jì)算得到的晶格常數(shù)（a=5.753 ，b=14.107 ）更接近試驗(yàn)值（a=5.742 ，b=14.079 ）.基于GGA和GGA+U的結(jié)構(gòu)優(yōu)化計(jì)算結(jié)果證明，試驗(yàn)得到的雙層鈣鈦礦Ba2CoOsO6中確實(shí)存在Os1和Co1反位缺陷.對(duì)試驗(yàn)給出的初始原子坐標(biāo)進(jìn)行計(jì)算修正，得到一個(gè)相對(duì)穩(wěn)定的基態(tài)結(jié)構(gòu)，為后續(xù)物理性質(zhì)的預(yù)測研究提供參考.

化學(xué)價(jià)是決定Ba2CoOsO6內(nèi)部電子分布的關(guān)鍵，這可以通過鍵長反映出來.例如，同為B位無序排列的雙層鈣鈦礦SrCaCoTeO6和SrCaNiTeO6中，Co2+與O2+之間距離的試驗(yàn)值為2.0～2.1 [11].而CHIN Y. Y.等[12]發(fā)現(xiàn)鈣鈦礦化合物中，CoO的鍵長為1.928 ，該鍵為Co2+高、低自旋狀態(tài)轉(zhuǎn)變的分界點(diǎn)[12].這是因?yàn)橛行Ьw場存在大小，所以低自旋狀態(tài)和高自旋狀態(tài)在一定距離上是穩(wěn)定的.采用GGA、GGA+U和考慮反位缺陷的GGA+U進(jìn)行計(jì)算時(shí)，CoO的3種鍵長（見表3）都遠(yuǎn)大于1.928 ，說明Co2+的高自旋狀態(tài)具有高穩(wěn)定性.初步研究表明，Ba2CoOsO6中的Co可以是正二價(jià)的氧化態(tài)，且為高自旋狀態(tài).

2.2磁學(xué)性質(zhì)分析

首先進(jìn)行了非自旋極化和自旋極化的計(jì)算.結(jié)果表明，雙層鈣鈦礦處于自旋極化狀態(tài)下的體系總能量更低、更穩(wěn)定.因此，后文的電、磁學(xué)計(jì)算均基于自旋極化狀態(tài).此外，Co和Os均為過渡金屬，影響鈣鈦礦整體物理性質(zhì)的主要因素通常是Co和Os的相互作用.根據(jù)文獻(xiàn)[11]提出的反鐵磁耦合構(gòu)型，對(duì)晶體結(jié)構(gòu)沿X、Z軸方向擴(kuò)胞，X、Y、Z軸方向的比為2∶1∶2.擴(kuò)胞后Ba2CoOsO6磁結(jié)構(gòu)如圖2所示.其中，紫色晶胞和藍(lán)色晶胞中Co磁矩方向分別為自旋向下和自旋向上，綠色晶胞中Os磁矩方向?yàn)樽孕蛏?，灰色晶胞中Co的磁矩方向?yàn)樽孕蛳?

不共線的磁矩在自旋方向上極化可能帶來不同分量，為了計(jì)算方便，沿同一軸向設(shè)置初始磁矩.由于材料中存在3d電子軌道和5d電子軌道的過渡金屬Co和Os，電子關(guān)聯(lián)效應(yīng)對(duì)材料的電、磁性能有重要影響，因此采用了GGA和GGA+U兩種計(jì)算方法.根據(jù)La2CoMnO6[13]和Sr2FeOsO6[14]鈣鈦礦體系中采用的UCo=4 eV， UOs=2 eV的Hubbard參數(shù)組合，進(jìn)行基于GGA+U方法的計(jì)算和分析結(jié)果的討論.Eg-為自旋向下的帶隙，Eg+為自旋向上的帶隙.帶隙及每個(gè)化學(xué)單位的總磁矩（mt）、Co磁矩（mCo）和Os磁矩（mOs）數(shù)據(jù)如表4所示，磁矩以玻爾磁子μB來表示.

由表4可知，Co、Os和單胞整體的磁矩對(duì)電子相關(guān)效應(yīng)非常敏感.采用添加U值的GGA+U方法后，Co的自旋磁矩（spim magnetic moment）從1.81μB增加到2.79μB，Os自旋磁矩從0.78μB減小到0.46μB，與試驗(yàn)值（0.28μB）較為接近.可見，添加U值后的磁矩更加合理.同時(shí)，由于存在3d電子軌道和5d電子軌道的過渡金屬Co離子和Os離子，因此進(jìn)一步考慮體系內(nèi)可能存在的SOC作用，以及SOC作用對(duì)電磁性質(zhì)的影響.

2.3電子結(jié)構(gòu)分析

圖3為基于GGA方法得到的態(tài)密度（density of states，DOS）圖，其中紅色虛線為費(fèi)米能級(jí)，設(shè)為能量零點(diǎn).由圖3可知：采用GGA計(jì)算方法時(shí)，雙層鈣鈦礦Ba2CoOsO6的上、下通道中態(tài)密度都穿過費(fèi)米能級(jí)，表現(xiàn)為金屬行為；費(fèi)米能級(jí)上，Co 的3d電子軌道和Os的5d電子軌道發(fā)生強(qiáng)烈的雜化作用，對(duì)費(fèi)米能級(jí)附近的總態(tài)密度有重要貢獻(xiàn).

圖4為基于GGA+U方法的Ba2CoOsO6能帶結(jié)構(gòu)（band structures）和態(tài)密度（projected density of state， PDOS）圖，其中G、F、Q、Z皆為高對(duì)稱點(diǎn).

加U后，未占據(jù)態(tài)被推到更高的能量區(qū)域，造成自旋通道的改變.在GGA+U的計(jì)算中，Ba2CoOsO6表現(xiàn)為窄帶隙的半導(dǎo)體.利用GGA+U方法，在自旋向上和自旋向下軌道中，價(jià)帶（valence band）最大值Evb，max位于高對(duì)稱點(diǎn)，與導(dǎo)帶（conduction band）最小值Ecb，min所處位置并非同一個(gè)點(diǎn)，從而證明Ba2CoOsO6是間接帶隙的半導(dǎo)體.通過GGA+U計(jì)算得到的Co、Os、O分波態(tài)密度圖如圖5所示.

圖5a中，Co的3d電子軌道除了自旋向下的dx2-y2、dz2、dxy未被占據(jù)，其他自旋向上及自旋向下的軌道全被占據(jù)，表明Co的化合價(jià)為正二價(jià)，處于3d7： t52ge2g電子構(gòu)型，Ba2CoOsO6高自旋極化.圖5b中，扭曲的OsO6八面體將5個(gè)5d電子軌道分解為非簡并的z2、雙簡并的{x2-y2，xy}.自旋向下通道z2、{x2-y2，xy}和{xz，yz}的電子占據(jù)數(shù)分別為1.00、0.96和0，自旋向下通道全部被占據(jù)，表明Os的氧化態(tài)化合價(jià)為正六價(jià)，電子構(gòu)型為5d2.綜上可知，雙層鈣鈦礦Ba2CoOsO6的電荷分布為Ba2+2Co2+Os6+O2-6.

進(jìn)一步考慮可能存在自旋-軌道耦合作用，基于GGA+U+SOC計(jì)算方法，根據(jù)第一性原理計(jì)算方法，計(jì)算得到的磁矩和帶隙匯總?cè)缦拢篛s和Co自旋磁矩分別為0.39μB和2.64μB；Os和Co軌道磁矩（orbital magnetic moment）分別為0.10μB和0.21μB；帶隙為0.16 eV.

由表4以及計(jì)算得到的磁矩和帶隙匯總可知，采用GGA+U方法計(jì)算得到的自旋向下和自旋向上的帶隙分別0.37 eV和0.40 eV，而采用GGA+U+SOC方法得到的兩個(gè)帶隙急劇減小，均只有0.16 eV，表明自旋-軌道耦合作用對(duì)電子結(jié)構(gòu)和電磁性質(zhì)的影響不可忽視.以上計(jì)算結(jié)果與文獻(xiàn)[11]的試驗(yàn)結(jié)果更接近，證實(shí)了本研究中的理論模擬結(jié)果是準(zhǔn)確的.

圖6為基于GGA+U+SOC方法的雙層鈣鈦礦Ba2CoOsO6的態(tài)密度圖.由圖6可知，Co的3d電子軌道和O的2p電子軌道在費(fèi)米能級(jí)附近產(chǎn)生強(qiáng)烈的雜化作用，表明角共享的Co6八面體對(duì)體系的電子輸運(yùn)性質(zhì)起重要作用.

3結(jié)論

1） 利用第一性原理計(jì)算，系統(tǒng)地研究了B位和B′位的Os1和Co1無序排列的雙層鈣鈦礦Ba2CoOsO6晶體結(jié)構(gòu)、反鐵磁耦合構(gòu)型和電子結(jié)構(gòu)，證實(shí)了其存在Os1和Co1反位缺陷.

2） 進(jìn)一步研究了摻雜多重自旋態(tài)的3d電子軌道過渡金屬Co和5d電子軌道過渡金屬Os的雙層鈣鈦礦Ba2CoOsO6電子結(jié)構(gòu)，證實(shí)了Ba2CoOsO6的電荷分布為Ba2+2Co2+Os6+O2-6.

3） 自旋-軌道耦合的存在對(duì)Ba2CoOsO6的電磁性質(zhì)具有一定的影響，降低了Co離子和Os離子的自旋磁矩.更重要的是Ba2CoOsO6具有獨(dú)特的反鐵磁耦合構(gòu)型，為研究其他雙層鈣鈦礦的電磁性能提供了新方法.此外，具有非常規(guī)的電子結(jié)構(gòu)和磁性性質(zhì)的雙層鈣鈦礦Ba2CoOsO6作為一種新型材料值得進(jìn)一步研究與開發(fā)應(yīng)用.

參考文獻(xiàn)（References）

［1］PHILIPP J B， MAJEWSKKI P， ALFF L， et al. Structural and doping effects in the halfmetallic double perovskite A2CrWO6 （A=Sr， Ba， and Ca）[J]. Physical Review B， DOI：10.1103/PhysRevB.68.144431.

［2］TOKURA Y， TOMIOKA Y， KUWAHARA H， et al. Origins of colossal magnetoresistance in perovskitetype manganese oxides[J]. Journal of Applied Physics， 1996，79（8）：5288-5291.

［3］KATO H， OKUDA T， OKKIMOTO Y， et al. Metallic ordered doubleperovskite Sr2CrReO6 with maximal Curie temperature of 635 K[J]. Applied Physics Letters， 2002，81（2）：328-330.

［4］FENG H L， ARAI M， MATSUSHITA Y， et al. Hightemperature ferrimagnetism driven by lattice distortion in double perovskite Ca2FeOsO6[J]. Journal of the American Chemical Society， 2014，136（9）：3326-3329.

［5］YAN B H， PAUL A K， KANUNGO S， et al. Latticesitespecific spin dynamics in double perovskite Sr2CoOsO6[J]. Physical Review Letters， DOI：10.1103/PhysRevLett.112.147202.

［6］ZHAO Q， LIU M， DAI J H， et al. Highpressure synthesis and ferrimagnetic ordering of the Bsiteordered cubic perovskite Pb2FeOsO6[J]. Inorganic Chemistry， 2016，55（19）：9816-9821.

［7］GE Z Z， TIAN Y， SUN A， et al. The impact of Bsite antisite defects on the magnetic and electronic properties in double perovskite Pb2FeOsO6[J]. Ceramics International， 2020，47（1）：992-1001.

［8］FENG H L， ADLER P， REEHUIS M， et al. Hightemperature ferrimagnetism with large coercivity and exchange bias in the partially ordered 3d/5d hexagonal perovskite Ba2Fe1.12Os0.88O6[J]. Chemistry of Materials， 2017，29（2）：886-895.

［9］HAFNER J.Abinitio simulations of materials using VASP： densityfunctional theory and beyond[J]. Journal of Computational Chemistry， 2008，29（13）：2044-2078.

［10］KAYSER P，ORLANDI F， AVDEEV M， et al. Structural chemistry and magnetic properties of the hexagonal double perovskite Ba2CoOsO6[J].Inorganic Chemistry， 2020，59（9）：6613-6622.

［11］ZARAQ A， ORAYECH B， FAIK A， et al. High temperature induced phase transitions in SrCaCoTeO6 and SrCaNiTeO6 ordered double perovskites[J]. Polyhedron， 2016，110：119-124.

［12］CHIN Y Y， HU Z W， SHIMAKAWA Y， et al. Charge and spin degrees of freedom in Asite ordered YCu3Co4O12 and CaCu3Co4O12[J]. Physical Review B， DOI：10.1103/PhysRevB.103.115149.

［13］SAHADASGUPTA T， BAIDYA S. Electronic structure and phonons in La2CoMnO6： a ferromagnetic insulator driven by Coulombassisted spinorbit coupling[J]. Physical Review B， DOI：10.1103/PhysRevB.84.035131.

［14］WANG J， ZU N N， HAO X F， et al. Firstprinciple investigations of the peculiar magnetic behavior of Sr2FeOsO6[J]. Physica Status SolidiRapid Research Letters， 2014，8（9）：776-780.

（責(zé)任編輯趙鷗）

收稿日期： ""2022-02-25

基金項(xiàng)目： "國家自然科學(xué)基金資助項(xiàng)目（21301075，21901086）； 江蘇大學(xué)高級(jí)人才基金資助項(xiàng)目（12JDG096，19JDG017）

作者簡介： ""李洪蘋（1984—），女，江蘇鎮(zhèn)江人，副研究員（hpli@ujs.edu.cn），主要從事多功能磁電材料的研究.

孫安（1996—），女，江蘇鎮(zhèn)江人，碩士研究生（sun_an1996@126.com），主要從事鈣鈦礦電磁性質(zhì)第一性原理計(jì)算的研究.