立方与四方相NaCaF3钙钛矿的比较研究：结构、力学、电子及热电(thermoelectric)性质

《Next Materials》：Comparative study of cubic and tetragonal NaCaF3 perovskites: Structural, mechanical, electronic, and thermoelectric properties

　　

摘要：本研究旨在通过第一性原理密度泛函理论(DFT, Density Functional Theory)及玻尔兹曼输运理论(Boltzmann transport theory)，系统比较并研究立方相(c)与四方相(t) NaCaF3钙钛矿化合物的结构、力学、振动、电子及载流子依赖的热电(thermoelectric)性质，重点关注四方畸变(tetragonal distortion)的影响。研究人员采用第一性原理计算结合玻尔兹曼输运理论研究了NaCaF3-c与NaCaF3-t钙钛矿化合物的载流子依赖热电性质。计算所得负形成能(formation energy)表明NaCaF3-c与NaCaF3-t均具结构稳定性；声子能带结构中无虚频(negative frequency)，证明所研究化合物具动力学稳定性(dynamic stability)；计算所得弹性常数满足玻恩判据(Born criteria)，确保力学稳定性(mechanical stability)。由于具宽间接带隙(indirect bandgap)，研究人员强调其在高功率及高温环境下的应用潜力及作为取代掺杂(substitutional doping)基体材料(host material)的适用性。采用广义梯度近似(GGA, Generalized Gradient Approximation)与HSE(Hey–Scuseria–Ernzerhof)泛函计算，NaCaF3-c的间接电子带隙分别为4.99 eV与6.78 eV，NaCaF3-t的GGA(HSE)计算间接电子带隙为4.86 eV(6.58 eV)。材料表现绝缘行为(insulating behavior)且具间接带隙，是电子器件的有前景候选材料。NaCaF3-c与NaCaF3-t具宽带隙且适合作为缺陷调控结构及电子性质以匹配各类器件应用的候选材料。预测所研究化合物具低晶格热导率(lattice thermal conductivity, κlat)，其中NaCaF3-t在500 K时κlat最小值为2 W·m?1·K?1。此外，电子(空穴)掺杂考察揭示良好的热电性能，热电优值(figures-of-merit, ZT)接近1，表明其作为热电器件应用的高效性。

论文解读：立方与四方相NaCaF₃钙钛矿的结构、力学、电子及热电性质比较研究

研究背景与意义： ABX₃型钙钛矿因可调控的组成与性质广泛应用于光伏、传感及热电领域，氟化物钙钛矿较氧化物具更宽带隙与更低声子频率。NaCaF₃中一价Na⁺与二价Ca²⁺共占A位，Goldschmidt容忍因子(tolerance factor)接近立方与四方相边界，易受热或应变诱导发生相变，但目前缺乏对其四方相及两相对热电响应比较的系统性理论研究。该论文（发表于《Next Materials》）采用第一性原理DFT与玻尔兹曼输运方程(BTE, Boltzmann Transport Equation)，首次对立方相(Pm3m)与四方相(I4/mcm, 引入c/a≈1.05畸变) NaCaF₃进行结构—力学—电子—热电综合比较，阐明四方畸变对带隙、弹性、声子散射及热电优值(ZT, figure of merit)的影响，证实其可作宽带隙绝缘基体及p型热电候选材料。

主要关键技术方法： 研究人员使用Quantum ESPRESSO软件包，采用GGA-PBE及HSE06杂化泛函，Na/Ca/F模守恒赝势(pseudopotential)，平面波截断能70 Ry、电荷密度300 Ry，Γ中心16×16×16 k点网格做几何优化（力收敛阈≤0.003 Ry/Bohr），20×20×20做电子结构与输运计算。通过Birch–Murnaghan状态方程拟合E-V曲线获平衡晶格常数与体弹模量(bulk modulus, B_o)；用微扰法计算弹性常数并验证Born–Huang稳定性；用PHONOPY基于2×2×2超胞有限位移法(±0.01 ?)算二阶力常数(IFCs-2)得声子谱与声子态密度(PDOS)，用PHONO3PY算三阶力常数(IFCs-3)求解单模弛豫时间近似(SMRTA, single mode relaxation time approximation)下晶格热导率κ_lat(T)；用玻尔兹曼输运理论在恒定弛豫时间近似(CRTA, constant relaxation time approximation, τ=10 fs)下计算Seebeck系数(S)、电导率σ/τ、电子热导率κ_e/τ及功率因子(PF, power factor)，结合κ_lat求ZT值；带边抛物拟合求电子(m_e^*)与空穴(m_h^*)有效质量。

研究结果：

3.1 立方与四方相NaCaF₃钙钛矿的结构性质(Structural properties of NaCaF₃perovskites (cubic and tetragonal phase))： 研究人员通过vc-relax获得GGA优化晶格常数a_o(c)=4.44 ?，a_o(t)=4.41 ?、c_o(t)=4.63 ?（HSE下略小）。四方畸变使ab面a减小、c轴伸长，CaF₆八面体倾斜(octahedral tilting)，赤道Ca–F键缩短、轴向Ca–F键拉长。形成能E_f分别为?3.41 eV(NaCaF₃-c)与?3.26 eV(NaCaF₃-t)，均为负值，表明两相热力学与结构稳定。体弹模量B_o分别为51.96 GPa(立方)与44.47 GPa(四方)。

力学性质与弹性常数： 立方相应满足C₁₁>0、C₄₄>0、C₁₁?C₁₂>0、C₁₁+2C₁₂>0；四方相满足相应四端不等式。NaCaF₃-c：C₁₁=109.58 GPa，C₁₂=15.23 GPa，C₄₄=12.96 GPa；NaCaF₃-t：C₁₁=126.00 GPa，C₁₂=17.26 GPa，C₁₃=20.42 GPa，C₃₃=80.10 GPa，C₄₄=13.82 GPa。剪切模量G(Hill平均)分别为22.45 GPa与24.93 GPa；杨氏模量Y分别为82.71 GPa与88.31 GPa；泊松比ν分别为0.29与0.32；B/G比分别为2.08与2.47，均>1.75，判定两相均为延性(ductile)。德拜温度Θ_D分别为400.46 K与394.83 K。各向异性因子A由立方0.27降至四方0.25，微裂纹风险低。

声子能带与晶格热导率(Phonon band structure and lattice thermal conductivity)： 基于2×2×2超胞计算的声子能带全布里渊区无虚频，证实动力学稳定；1000 K下20 ps AIMD(Ab Initio Molecular Dynamics)总能量围绕基态涨落无异样，确证热力学稳定。四方相因CaF₆八面体倾斜致对称性降低、非谐性增强及Umklapp声子—声子散射加强，并引入低频光学模致声学模共振散射("phonon rattling")，使κ_lat低于立方相。预测室温κ_lat为3.12 W·m^?1·K^?1(立方)与约4.19 W·m^?1·K^?1(四方，原文数据有笔误此处按正文叙述四方更低，文中给出最小值在500 K为NaCaF₃-t κ_lat=2 W·m^?1·K^?1@500 K，立方相相应值略高)。

4.1 电子性质(Electronic properties)： GGA(HSE)算得NaCaF₃-c间接带隙R→Γ为4.99 eV(6.78 eV)，NaCaF₃-t为4.86 eV(6.58 eV)，均为宽带隙绝缘体。分波态密度(PDOS, Projected Density of States)显示价带顶(VBM, Valence Band Maximum)主要由F 2p轨道贡献并混有少量Ca 3p，导带底(CBM, Conduction Band Minimum)主要来自Na 3s及少量Ca 3d/4s。四方畸变轻微缩小带隙并使价带顶变平(flatter VBM)。有效质量拟合得NaCaF₃-c：m_e^*=0.45 m₀，m_h^*=2.31 m₀；NaCaF₃-t：m_e^*=0.51 m₀，m_h^*=2.98 m₀。空穴有效质量显著大于电子，m_h^*/m_e^*比值(D)分别为5.13与5.84。重空穴源自平坦价带致大Seebeck系数。

4.2 热电性质(Theroelectric properties)： 在n(p)掺杂浓度10¹⁸–10²¹cm^?3、温度300–800 K范围内，用CRTA(τ=10 fs)计算输运系数。p型掺杂因重空穴具更高Seebeck系数(S(T,p)>S(T,n))；功率因子PF=σS²随掺杂浓度增至最优值（n型~5×10²⁰cm^?3，p型~9×10²⁰cm^?3）而升高。晶格热导率κ_lat(t)< />_lat(c)，电子热导率κ_e(T,p)增长慢于κ_e(T,n)。ZT=S²σT/(κ_e+κ_lat)，p型NaCaF₃-t在800 K、空穴浓度9×10²⁰cm^?3下达峰值ZT≈0.92(文中Table 5脚注提示可达近1.0)，立方相p型ZT≈0.81@800 K。n型掺杂ZT较低。四方畸变通过增大m_h^*（提S）与降κ_lat（降总κ）共同提升ZT。

讨论与结论总结（翻译自Conclusion）： 研究人员通过GGA与HSE模拟研究了NaCaF₃钙钛矿化合物的结构、弹性、电子及热电性质，考察了四方畸变的影响。负形成能证实NaCaF₃-c与NaCaF₃-t的结构稳定性，弹性常数满足Voigt–Reuss–Hill与Born–Huang判据确保力学稳定性。优化晶格常数立方相为4.44 ?、四方相a=4.41 ?、c=4.63 ?。电子带隙立方相GGA(HSE)为4.99 eV(6.78 eV)、四方相为4.86 eV(6.58 eV)，表明绝缘行为。PDOS显示VBM主由F?2p贡献，CBM主由Na?3s影响。热电研究表明ZT强烈依赖温度与掺杂浓度，p型掺杂表现优异。NaCaF₃的超宽带隙特性暗示其在超越GaN与SiC的高功率电子器件中的应用潜力，且适合作为取代掺杂调控带隙与性质的基体材料。本研究首次详细报道四方相NaCaF₃热电性能，证明四方畸变通过增大空穴有效质量及降低晶格热导率提升p型热电响应，预测空穴掺杂NaCaF₃-t的ZT趋近1，使其成为高温废热回收有前景的候选材料，也为后续实验合成与掺杂研究奠定理论基础。