探索具有奇異物理特性的新材料是凝聚態物理和材料科學的前沿領域。其中，過渡金屬氧化物，尤其是銅氧化物（銅酸鹽），因在高溫超導方面的突破性發現而備受矚目。然而，實現類似銅酸鹽中關鍵的Cu²⁺(3d⁹)電子構型在鎳基氧化物中卻面臨挑戰。傳統的鈣鈦礦型稀土鎳酸鹽RNiO₃中的Ni³⁺(3d⁷)離子處于電荷轉移區，其八面體配位環境常導致電荷歧化、金屬-絕緣體轉變等復雜現象，與銅酸鹽的理想環境相去甚遠。為此，科學家們將目光投向了氧缺陷型層狀鎳酸鹽。通過減少氧含量，降低鎳的配位數，有望將鎳的價態推向接近Cu²⁺的Ni¹⁺(3d⁹)。在這一系列材料中，三明治層狀的Ruddlesden-Popper (RP)相R₄Ni₃O₈(R為稀土元素)尤為引人關注。它由三層平面四方NiO₂層與巖鹽結構的稀土-氧阻擋層交替堆疊而成，其結構類似于銅酸鹽的CuO₂平面，但缺乏頂端氧原子，形成了平面四方NiO₄配位。這種結構變化如何影響其電子結構、磁性及熱力學性質？它與三維鈣鈦礦鎳酸鹽以及二維銅酸鹽體系有何異同？為回答這些問題，來自印度Jiwaji大學的Surendra Singh Sengar與Dinesh C. Gupta研究員對Sm₄Ni₃O₃和Eu₄Ni₃O₈展開了系統的第一性原理研究，相關成果發表在《Computational Condensed Matter》上。

為完成這項研究，作者主要采用了以下幾種關鍵的計算技術方法：首先，使用基于全勢線性綴加平面波(FP-LAPW)方法的WIEN2k軟件包進行第一性原理密度泛函理論(DFT)計算。在交換關聯泛函上，不僅使用了廣義梯度近似(GGA)的PBE形式，還引入了改進的Becke-Johnson (mBJ)勢以更精確地描述過渡金屬氧化物的電子能帶。為了處理Ni-3d和稀土-4f電子的強關聯效應，他們采用了DFT+U方法，對Ni-3d態施加了U_eff= 4.5 eV的有效哈伯德參數。其次，基于DFT優化得到的結構，使用BoltzTraP2代碼在半經典玻爾茲曼輸運理論框架下計算了電子的泡利磁化率。最后，利用Gibbs2軟件包中的準諧德拜模型(QHDA)研究了材料在0-900 K溫度和0-50 GPa壓力范圍內的熱力學性質，包括比熱、德拜溫度、體彈模量、熱膨脹系數和熵。同時，使用IRELAST方法計算了彈性常數，評估了材料的機械穩定性。

研究首先確定了兩種化合物的晶體結構。Sm₄Ni₃O₈和Eu₄Ni₃O₈結晶于四方晶系，空間群為I4/mmm，具有典型的層狀鈣鈦礦結構，由NiO₂平面層和稀土(Sm/Eu)-氧阻擋層交替排列構成。通過對比鐵磁(FM)、A型反鐵磁(A-AFM)和非磁(NM)構型的總能量-體積曲線，發現鐵磁構型是能量最低的基態，排除了反鐵磁有序作為最低能量構型的可能性。優化后的晶格參數與實驗報道的La₄Ni₃O₈等類似化合物吻合良好。內聚能和穩定能計算證實了這兩種化合物的熱力學穩定性。與八面體配位的鈣鈦礦SmNiO₃相比，Sm₄Ni₃O₈和Eu₄Ni₃O₈由于缺乏頂端氧原子，形成了平面四方NiO₄配位，導致更均勻的平面內Ni-O鍵長，增強了平面內的Ni-3d/O-2p雜化。

為了研究其自旋極化的電子特性，計算了能帶結構和態密度(DOS)。結果表明，在GGA、mBJ以及DFT+U框架下，兩種化合物的自旋向上和自旋向下能帶均穿過費米能級(E_F)，證實了其金屬特性。態密度分析揭示了顯著的Ni-3d與O-2p軌道雜化，而稀土(Sm/Eu)-4f態作為局域化特征出現在費米能級附近，對導電貢獻較小。DFT+U計算的加入并沒有打開能隙，費米能級處仍存在有限的態密度，表明其金屬性對電子關聯效應具有魯棒性。與三維鈣鈦礦鎳酸鹽不同，氧缺陷型RP相R₄Ni₃O₈由于形成平面四方NiO₄單元并降低了平均鎳價態（趨于Ni¹⁺/Ni²⁺混合態），抑制了典型的電荷歧化現象，促進了金屬行為。

電荷密度差分析表明，兩種化合物的成鍵性質兼具離子性和共價性特征。由于氧的電負性最高，電子從稀土和鎳向氧轉移，導致R-O (R = Sm, Eu)相互作用主要是離子性的。而Ni和O之間較小的電負性差異允許部分電子共享，從而在NiO₂平面內形成了混合的離子-共價Ni-O鍵。

磁性分析顯示，磁性主要來源于局域的稀土4f電子，而Ni原子貢獻了中等大小的磁矩。具體而言，在Sm₄Ni₃O₈中，Sm原子貢獻約5.1-5.2 μ_B/atom，Ni原子貢獻約1.0-1.1 μ_B/atom。在Eu₄Ni₃O₈中，Eu原子貢獻約6.2-6.3 μ_B/atom（與其半滿的4f⁷構型一致），而Ni的磁矩貢獻較小且在不同計算方法下變化。總能量比較證實鐵磁(FM)態是能量最低的基態，優于A型反鐵磁(A-AFM)和非磁(NM)態。與三維鈣鈦礦鎳酸鹽中由Ni-O-Ni超交換主導的強關聯反鐵磁不同，RP三明治層狀結構因氧缺陷和平面四方配位，削弱了經典超交換路徑，增強了平面內電子巡游性，從而穩定了鐵磁極化態。

基于DFT能帶結構，利用BoltzTraP2計算的溫度依賴的泡利磁化率χ_Pauli(T)顯示，兩種化合物的χ_Pauli隨溫度升高略有下降，這是金屬體系費米面附近電子態熱展寬的典型特征。Sm₄Ni₃O₈的χ略高于Eu₄Ni₃O₈，與其相對更高的費米能級態密度一致。需要強調的是，此處的χ是源自費米面附近巡游電子的泡利順磁磁化率，并非局域磁矩的居里-外斯磁化率。

通過準諧德拜模型(QHDA)研究了材料在0-900 K溫度和0-50 GPa壓力范圍內的熱力學性質。定壓熱容(C_p)和定容熱容(C_v)在低溫區(0–300 K)隨溫度迅速上升，在高溫區逐漸飽和并接近經典杜隆-珀蒂極限。壓力增加會抑制晶格振動，導致C_p和C_v在所有溫度下輕微降低。德拜溫度(θ_D)隨溫度升高略微下降，但隨壓力增加而升高，反映了晶格剛性的增強。體彈模量(B)隨溫度升高而下降（熱軟化），隨壓力增加而升高。熱膨脹系數(α)隨溫度升高而增加，但隨壓力增加而減小。熵(S)隨溫度單調增加。格林艾森參數(γ)隨溫度和壓力增加均略有下降，表明非諧晶格振動的減弱。總體而言，Sm₄Ni₃O₈展現出比Eu₄Ni₃O₈略高的C_p、θ_D和B值，表明其晶格剛性更強，而Eu₄Ni₃O₈的熱膨脹系數略高。

通過計算獨立的彈性常數并驗證其滿足四方晶系的玻恩穩定性判據，確認了Sm₄Ni₃O₈和Eu₄Ni₃O₈具有機械穩定性。采用Voigt-Reuss-Hill (VRH)近似得到的多晶模量顯示，兩者的體彈模量B約為131 GPa。Sm₄Ni₃O₈的剪切模量G和楊氏模量E略高于Eu₄Ni₃O₈，表明其剛性更強。泊松比ν在0.29–0.30之間，Pugh比B/G在1.87–2.02之間，表明材料具有韌性到臨界韌性的力學行為。彈性常數C₃₃< C₁₁反映了層狀RP結構沿c軸方向的彈性響應較軟，即各向異性特征。

本研究通過綜合運用第一性原理密度泛函理論(DFT)、改進的Becke-Johnson (mBJ)勢、DFT+U方法以及準諧德拜模型(QHDA)，系統揭示了層狀稀土鎳酸鹽Sm₄Ni₃O₈和Eu₄Ni₃O₈的結構、電子、磁性和熱力學性質。主要結論如下：首先，兩種化合物均結晶于四方層狀鈣鈦礦結構(I4/mmm空間群)，鐵磁(FM)態是其能量最低的磁基態。其次，電子結構計算明確顯示它們具有金屬性，費米能級附近的態主要源于Ni-3d和O-2p軌道的強雜化，而稀土4f態高度局域化。這種金屬性在考慮電子關聯的DFT+U計算中依然穩健，表明氧缺陷和平面四方配位抑制了鈣鈦礦鎳酸鹽中常見的電荷歧化和金屬-絕緣體轉變。第三，磁性主要由局域的稀土(Sm/Eu) 4f電子主導，Ni貢獻了中等大小的磁矩，共同穩定了鐵磁基態。第四，熱力學分析表明，比熱在高溫下趨近杜隆-珀蒂極限，德拜溫度和體彈模量隨壓力增加而升高，證實了晶格在高壓下的硬化行為。最后，彈性常數計算證實了兩種化合物具有良好的機械穩定性。

這項研究的重要意義在于，它從理論上詳細闡述了從三維鈣鈦礦鎳酸鹽到準二維銅酸鹽類似物之間的結構-性質關聯。通過闡明氧缺陷、降低的配位數（從NiO₆八面體到NiO₄平面四方）以及混合鎳價態如何協同作用，抑制電荷有序、促進金屬性并改變磁交換路徑，該工作為理解層狀過渡金屬氧化物的維度與關聯效應提供了關鍵見解。所揭示的金屬性、鐵磁序、可調的熱力學響應以及良好的機械穩定性，共同指向Sm₄Ni₃O₈和Eu₄Ni₃O₈在自旋電子學和熱電材料等領域的潛在應用價值，為后續實驗探索和材料設計提供了堅實的理論基礎和預測指南。