钛酸钡的电子结构计算

结合案例图分析，我把这个项目的完整过程、参数设置、结果解读和踩坑经验整理如下。

项目过程：钛酸钡的电子结构计算

客户做的是钛酸钡（BaTiO₃）钙钛矿，想解释它的铁电性起源和能带特征，用于一篇关于无铅压电材料的论文。审稿人要求补充电子结构图：能带结构和分波态密度（PDOS）。

软件选 VASP 6.3.2，泛函 PBE，赝势 PAW。钛酸钡的 Ti 3d 电子关联比较强，直接用 PBE 会低估带隙，所以我加了 PBE+U：U(Ti 3d)=4.0 eV，J=0.5 eV。体系用 1×1×1 原胞，晶格常数先用实验值 a=3.99 Å 做固定晶胞计算，随后也做了结构优化对照。

关键参数：ENCUT=520 eV，K 点 8×8×8，ISMEAR=-5（绝缘体/半导体用四面体法），LORBIT=11 输出 PDOS，NEDOS=3000。能带路径取 R–Γ–X–M–Γ。

案例分析：从数据到结论

第一张案例图是 BaTiO₃ 的能带结构。费米能级设在 0 eV，价带顶（VBM）位于 Γ 点，主要由 O 2p 轨道构成；导带底（CBM）也位于 Γ 点，主要由 Ti 3d 构成。因此这是一个直接带隙半导体，带隙宽度约 E_g≈1.9 eV。需要提醒的是，PBE+U 给出的 1.9 eV 仍然低于实验光学带隙（约 3.2 eV），这是 GGA 类泛函的系统性低估；但在解释价带/导带组成和铁电畸变趋势上，这个结果已经够用。

第二张案例图是 PDOS。O 2p（蓝色填充）在 -6 到 0 eV 之间形成一个宽带；Ti 3d 分成两组：t₂g（约 2–4 eV）和 e_g（约 5–7 eV）。Ti 3d 和 O 2p 之间有显著重叠，说明 Ti–O 键具有强共价/配位特征，这也是铁电畸变时 Ti 离子偏移导致自发极化的电子结构根源。

总结教训：复盘必须警惕的三件事

第一，+U 参数要合理。我试过 U=2 eV 和 U=6 eV，带隙分别是 1.4 eV 和 2.4 eV；U=4 eV 与文献报道的 PBE+U 结果最接近，而且不影响晶格常数太多（优化后 a=4.02 Å，实验 3.99 Å，偏差 0.8%）。

第二，能带路径必须覆盖高对称点。如果只取 Γ–X，会漏掉 R 点和 M 点的价带顶细节，误导带隙类型判断。

第三，BaTiO₃ 的铁电性来自 Ti 离子偏离氧八面体中心，计算时不能忽略结构优化。我对比了立方相和四方相：立方相没有带隙处的分裂，而四方相由于 Ti 位移导致价带顶轻微劈裂，这与实验观测一致。

项目经验：下次再做我会怎么做

下次做钙钛矿电子结构，我会把流程标准化：1) PBE 初算看趋势；2) PBE+U 校正 d/f 电子关联；3) 必要的话用 HSE06 算更准的带隙（成本约 PBE 的 5–10 倍）；4) 同步输出能带+PDOS+Born 有效电荷，三者一起解释铁电/压电机制。对于审稿人，我会在图注里明确标注”PBE+U(U=4 eV)”，避免被质疑泛函选择。

钛酸钡的电子结构计算

计算工程师团队：500+ ｜ 成功计算案例：5W+ ｜ 1W+ 助力顶刊发布经验 需要模拟代算服务可以联系我们，可免费评估项目的计算费用和免费出计算方案。