立方钙钛矿BaPdO3的结构、热物理、电子和光学性质的第一性原理研究

《Materials Today Communications》：First principles investigation of the structural, thermophysical, electronic and optical properties of cubic perovskite BaPdO 3

　　

在本工作中，研究人员首次利用密度泛函理论（DFT）计算研究了立方钙钛矿BaPdO₃，以推测其弹性、热学、电子和光学性质。BaPdO₃的负形成能证实了其热力学稳定性。对计算得到的弹性常数和模量的详细分析表明，BaPdO₃具有机械稳定性、延展性、可加工性、中等硬度和弹性各向异性。研究人员根据弹性常数预测了BaPdO₃的德拜温度和熔化温度。对电子电荷密度和Mulliken键布居的分析表明，其具有离子-共价混合键合特征，其中离子相互作用起主导作用。计算的电子能带结构（包括和不包括自旋-轨道耦合（SOC））以及态密度均证实了BaPdO₃的金属本性。费米面由类电子和类空穴片层组成。电导率曲线表明其具有高导电性。在有限温度下，非谐效应稳定了在0 K时虚设的声子，从而导致稳定相。计算的光学参数与能带结构和态密度的特征一致。BaPdO₃在可见光范围内表现出高反射率，在紫外区域表现出强吸收，并在低光子能量下具有高折射率。这些电子和光学性质表明，BaPdO₃可能有望用于需要金属行为的应用中，例如电极和反射涂层。

**论文解读：立方钙钛矿BaPdO₃结构与物性的第一性原理研究**

**1. 研究背景与问题**

近年来，第一性原理计算方法的发展极大地拓展了凝聚态物理研究的范围，能够以高精度计算固体的电子结构和结构性质。先进功能材料对于当代和未来技术应用至关重要。ABO₃钙钛矿材料因其多样的电学特征和固态现象（如金属性、超导性、半导体性和绝缘性）而受到广泛关注，并在波导、激光倍频、衬底、传感器、燃料电池、存储单元、光电子和自旋电子学等领域具有应用前景。其中，金属性钙钛矿材料因其高电导率、高电流负载能力、高硬度和独特的晶体结构，有望用于电极材料。然而，寻找与钙钛矿结构及化学性质兼容、且具有极低电阻率的高质量电极材料仍是发展功能氧化物电子器件的关键挑战。立方钙钛矿SrMoO₃、BaMoO₃、PbMoO₃、BaRuO₃、BaNbO₃等已被研究或用作电极材料。尽管Li等人曾报道过立方BaPdO₃的容差因子和稳定性，但该化合物的详细物理特性尚属未知。因此，研究人员首次通过第一性原理计算，系统研究BaPdO₃的结构、弹性、电子、动力学和光学性质，并评估其潜在应用。

**2. 主要技术方法**

研究人员使用基于密度泛函理论（DFT）的CASTEP代码进行计算。采用局域密度近似（LDA）下的CA-PZ形式描述交换关联效应，使用模守恒赝势和Koelling-Harmon相对论方案处理离子-电子相互作用。通过BFGS最小化技术优化晶体结构，采用Monkhorst-Pack方案生成12×12×12的k点网格，平面波截断能为1100 eV。对于费米面和自旋-轨道耦合（SOC）计算，分别使用18×18×18和6×6×6的k点网格。弹性常数通过应力-应变方法计算，Mulliken和Hirshfeld布居分析用于研究键合性质。声子计算采用有限位移法结合密度泛函微扰理论，使用12×12×12的q点网格；有限温度声子色散通过ALAMODE包中的自洽声子（SCP）形式获得，并采用2×2×2超胞。电导率使用BoltzTraP2计算。

**3. 研究结果**

**3.1 结构分析**
BaPdO₃结晶为立方结构（空间群221，Pm-3m），Ba占据角位置，Pd占据体心，O占据面心。优化晶格常数为4.154 ?，体积为71.67 ?³。形成能为-1.407 eV/atom，证实其热力学稳定性。

**3.2 弹性与力学性质**
计算得到三个独立弹性常数C₁₁=242.7 GPa，C₁₂=116.03 GPa，C₄₄=77.54 GPa，满足Born-Huang力学稳定性判据。Kleinman参数为0.606，表明键弯曲贡献大于键拉伸。正的Cauchy压力（38.49 GPa）和B/G比（2.213）表明材料具有延展性。体积模量B=158.25 GPa，剪切模量G=71.51 GPa，杨氏模量Y=186.45 GPa。泊松比ν=0.30，显示离子-金属混合键合。可加工指数μ_M=2.04，表明材料具有良好的可加工性。硬度H=9.562 GPa，为中等硬度。各向异性因子A=1.22，表现为弹性各向异性。

**3.3 热力学性质**
通过声速计算得到德拜温度θ_D=445.57 K，熔化温度T_m=1446.24 K，均低于SrMoO₃和BaMoO₃等类似化合物。

**3.4 电子能带结构与态密度**
能带结构表明，费米能级处存在能带穿越，导带和价带重叠，无带隙，呈现金属性。包括SOC后，重元素Ba和Pd引起的能带分裂出现，但整体特征相似。总态密度（TDOS）在费米能级处非零，主要贡献来自O-2p和Pd-4d态。价带-11~-9 eV区域由Ba-5p态主导，-8~-3 eV区域由Pd-4d和O-2p杂化贡献，导带6~12 eV区域主要由Ba-6s和Ba-4d态占据。

**3.5 费米面**
费米面由五个能带（21-25）构成，包含类电子和类空穴片层。能带21和22在Γ-R方向出现类电子片层，能带23在Γ-R方向及X点附近有类电子片层，能带24在Γ-M方向有类空穴片层，能带25在Γ点附近有类电子片层。

**3.6 电导率**
计算的电导率随温度升高呈指数下降，但在800 K时仍保持显著导电性。50 K时，电导率约为8×10²⁰ (Ω·m)^-1·s^-1，表明BaPdO₃是良好的导体，有望用于燃料电池电极。

**3.7 布居分析**
Mulliken和Hirshfeld电荷分析显示Ba和Pd带正电，O带负电，存在离子键合。O-Pd键长最短（1.988 ?），重叠布居最大（0.57），表明共价键主要来自O-Pd相互作用。O-Ba重叠布居小（0.06），主要为离子键。O-O和Pd-Ba具有负重叠布居，呈反键特征。

**3.8 电荷密度图**
(001)和(002)晶面电荷密度图显示，Pd和O周围电荷积累较大，Ba周围较小。Pd和O之间存在少量电荷共享，表明弱共价键。Ba和O之间无明显电荷共享，主要为离子键。整体均匀的低电荷密度区域表明存在金属键。因此，BaPdO₃具有离子、金属和共价混合键合，以离子键为主。

**3.9 声子色散**
0 K下声子谱出现负频率，表明动态不稳定性，这与氧八面体倾斜和面内氧原子旋转畸变有关。然而，采用自洽声子（SCP）方法在300 K、600 K和900 K下计算声子色散后，虚频消失，表明非谐效应在有限温度下稳定了立方相。

**3.10 光学性质**
光学参数计算使用Drude阻尼0.05 eV和等离子体能量18 eV。介电函数实部ε₁(ω)在低能量下为大的负值，呈Drude响应；虚部ε₂(ω)在低能量下为大的正值，表明强吸收。吸收谱从零光子能量开始，证实金属性；在紫外区（15-25 eV）出现强吸收。光导率在低能区（1-5 eV）最高，并在31 eV附近趋近于零。反射率在零光子能量处最大，在红外和可见光区超过90%，表明可用作热障涂层和太阳热反射表面。折射率在低能下较高，随能量增加而降低。损失函数在31 eV处出现峰值，对应等离子体频率，超过此能量后材料变为透明，从金属性转变为介电性。

**4. 讨论与结论**

研究首次通过第一性原理计算全面揭示了BaPdO₃的结构、弹性、热学、电子、动力学和光学性质。负形成能确认了热力学稳定性，弹性常数满足Born-Huang判据，表明机械稳定性。材料表现为延性、中等硬度和良好的可加工性。电子能带结构和态密度证实其金属性，电导率计算显示高导电性，因此有望作为固体氧化物燃料电池的电极材料。0 K下的虚声子在有限温度下通过非谐相互作用被稳定，消除了动态不稳定性。光学性质进一步支持其金属行为，在可见光区具有高反射率，在紫外区具有有效吸收，潜在应用包括太阳能反射器、热障涂层和深紫外探测器。这些新确定的性质为未来的实验和理论研究奠定了坚实基础。