基于Wannier函数研究赝势与结构对过渡金属氧化物电子性质的影响

《RSC Advances》：Pseudopotential and structural effects on electronic properties in transition metal oxides using Wannier functions

　　

赝势选择如何影响第一性原理计算中的电子输运性质？原子结构又如何影响所计算的性质？本研究以若干过渡金属氧化物为体系，比较了六种不同赝势库对带隙、高频介电常数及电导率等计算材料性质的影响。研究人员同时指出，Wannier函数可用于计算材料的电子输运性质。所建立的工作流表明，无论采用何种赝势库，Wannier函数都为处理复杂体系提供了一种可靠且准确的方法。研究人员将结果与已有理论工作及文献中的实验值进行了验证与对比。研究发现，不同赝势对计算的带隙、电导率及高频介电常数的影响可忽略不计。因此，Wannier函数可用于对金属、半导体及磁性过渡金属氧化物的带隙与电导率计算进行基准测试。研究人员还发现，原子结构的选择（即实验结构 vs. 弛豫/优化结构）对电导率的影响并非可忽略。与弛豫原子结构相比，实验结构通常能给出与本征电导率实验值更为吻合的结果。

论文解读：基于Wannier函数研究赝势与结构对过渡金属氧化物电子性质的影响

研究背景与意义

过渡金属氧化物（Transition metal oxides, TMOs）因其独特的电子性质和复杂结构，在能源存储（如锂离子电池正极）、电子器件（传感器、晶体管、存储器件）及电催化与光催化领域具有重要应用。其性能调控高度依赖于过渡金属d轨道与氧p轨道的强共价相互作用。然而，针对TMOs电子输运性质的从头算（ab initio）建模仍具挑战，主要源于对复杂结构中电子行为的有限理解，以及密度泛函理论（DFT）中交换关联泛函与赝势（Pseudopotential, PP）选择带来的不确定性。此前的高通量研究多未系统评估不同赝势库的影响，且常使用弛豫结构而非实验结构进行对比。为此，研究人员开展本研究，以明确赝势库与原子结构选择对TMOs电子性质（带隙、高频介电常数、电导率）计算的影响，并验证基于最大局域Wannier函数（Maximally Localized Wannier Functions, MLWFs）工作流的可靠性。该论文发表于《RSC Advances》。

主要关键技术方法

研究人员选取多种TMOs体系（如TiO₂-R、TiO₂-A、MoO₂、MoO₃、V₂O₅、V₂O₃、Cr₂O₃、ZnO），使用Quantum ESPRESSO 7.2与Wannier90 3.1.0开展DFT计算。交换关联泛函采用Perdew–Zunger局域密度近似（LDA/PZ）、Perdew–Burke–Ernzerhof广义梯度近似（GGA/PBE）及带Hubbard +U修正（PBE+U_AFLOW）。研究对比六种常用赝势库：SG15-NC（模守恒）、GBRV-US（超软）、PS-PAW（Dal Corso投影缀加平面波）、PS-US（Dal Corso超软）、SSSP-efficiency v1.3.0、PseudoDojo v0.5。电子能带通过MLWFs进行插值：金属及小带隙体系采用选定密度矩阵列（SCDM）自动投影法；绝缘体及磁性体系（E_G>1 eV）采用投影解缠Wannier函数（PDWFs）手动设置。高频介电常数（ε^∞）通过Wannier电荷中心或Berry相方法计算；电导率通过BoltzWann代码在恒定弛豫时间近似（τ=10 fs）下求解玻尔兹曼输运方程获得。原子结构效应通过对比Materials Project弛豫结构与晶体学开放数据库（COD）实验结构进行评估。

研究结果

3.1 带隙（Band gaps）

研究人员计算了半导体TMOs的带隙。结果表明，所有计算带隙较实验值普遍低估，LDA/PZ低估最显著，PBE略有改善，PBE+U_AFLOW对磁性Cr₂O₃等体系提升明显。不同赝势库对带隙的影响较小：以MoO₃为例，各PP结果几乎重合；ZnO在PZ下PS-PAW与SG15-NC偏差最大（约40%），但PBE下降至14%，PBE+U进一步降至7%。PS-PAW与PS-US虽源自同一原子数据，带隙仍可略有差异。总体结论：在所研究的多数TMOs中，带隙对赝势选择具有鲁棒性。

3.2 高频介电常数（High-frequency permittivity）

以TiO₂-R（四方）和V₂O₅（正交）为例，研究人员对比了MLWF电荷中心法与Berry相法的结果。对于近各向同性的TiO₂-R，a/b轴ε^∞~6.95–7.10，c轴~8.56–8.59，六种PP结果一致，MLWF与Berry相差<1%。对于各向异性显著的V₂O₅，a轴6.37–6.60，b轴6.45–6.66，c轴3.8–4.0；MLWF与Berry在a轴偏差~10%，b轴~3%，c轴可忽略。使用实验原胞而非弛豫原胞时，与Materials Project参考值偏差增大，主要源于结构远离PBE平衡态及k点密度差异。结论：赝势选择对ε^∞影响可忽略；MLWF可可靠重复Berry相结果，但在未弛豫实验结构上需关注Wannier函数收敛性。

3.3 电导率（Electrical conductivity）

研究人员计算了金属（MoO₂、V₂O₃）、绝缘/半导体（MoO₃、TiO₂-R、TiO₂-A、ZnO、Cr₂O₃、V₂O₅）的电导率。赝势库变化对σ及活化能E_a(σ)影响甚微。结构选择影响显著：MoO₃弛豫结构高估活化能（~1.20 eV vs. 实验~0.60 eV），实验结构更准（~0.42 eV）；金属MoO₂对结构不敏感，均与实验吻合。TiO₂-A（400–800 K本征输运区）与TiO₂-R（850–1250 K）同样显示实验结构更接近单晶实验值（~1.44 eV vs. ~1.48 eV），弛豫结构偏低（~1.06 eV）。ZnO表现反常：实验结构反而不佳，因其满壳层3d¹⁰使O 2p–Zn 3d排斥对Zn–O键长敏感，导致带隙低估与迁移率低估的偶然抵消被打破；PBE+U改善带隙却降低电导率吻合度。Cr₂O₃自旋极化计算在实验结构下与实验合理吻合。V₂O₅预测偏差较大，可能仍处极化子跳迁区、样品为退火薄膜且易含氧空位、层状结构需色散修正；但其活化能约为实验带隙一半，本征输运趋势合理。V₂O₃作为金属，绝对值在1.5倍内，结构效应不显。此外，Hubbard修正并未系统改善E_a(σ)，带隙精度与电导率在恒定弛豫时间近似下关联较弱。

讨论与结论翻译

研究人员在本工作中测试了PZ、PBE及PBE+U密度泛函近似（DFAs）下赝势与原子结构对材料电子结构与电子输运性质（带隙、高频介电常数、电导率）的影响，并评估了利用最大局域Wannier函数计算这些性质的可行性。研究发现，赝势对带隙、电导率及高频介电常数的计算结果影响通常可忽略不计。不过，为使MLWF与Berry相方法完全自洽，需保证MLWFs良好收敛、占据带恰当解缠、使用原胞及弛豫结构。此外，研究人员利用基于Wannier函数的工作流，成功计算了具绝缘与金属特性及磁性体系的TMOs电导率温度响应，可用于材料筛选。研究表明，通过对比优化/弛豫结构与实验结构计算的结果，原子结构的选择对电导率具有不可忽略的影响。最后，研究人员考察了赝势参数化与所用交换关联泛函不一致的影响：对计算带隙低于~1 eV的材料（如ZnO），这种不一致可产生显著影响；其余情况多可忽略。本研究为未来在更广材料范围内对密度泛函近似进行基准测试奠定了基础。