理论计算dos

刚学理论计算dos 的时候，我对着一张总态密度图只会说”有带隙”，导师问我”哪个峰是 O 的 p 轨道”，我当场卡壳。那次之后我才明白，理论计算dos真正的用法是拆分到分波态密度（PDOS），把每个峰指认到具体原子和轨道，否则图再漂亮也只是装饰。

一、为什么理论计算dos要做 PDOS

总 DOS 把所有原子的贡献混在一起，看不出化学键本质。我算 BaTiO₃ 时，总 DOS 显示一个约 1.8 eV 的带隙，但说不清是谁贡献的；拆成 PDOS 才发现价带顶主要是 O 2p，导带底主要是 Ti 3d，Ba 几乎不参与成键。这种 p–d 杂化信息直接决定你后续怎么调控带边。做理论计算dos 时，LORBIT=11（VASP）或对应设置输出 lm 分解，是基本操作。

二、核心原理：从本征值到态密度

理论计算dos 的本质，是把 Kohn-Sham 本征值按能量直方图化，并对 k 点加权。布里渊区积分靠足够密的 k 网格（我常用 0.020.04 Å⁻¹ 间距），否则尖峰会漏采样。为了平滑，用 Gaussian 或 tetrahedron 展宽，展宽宽度 0.050.1 eV 较稳妥。我早期用太大展宽把窄 d 峰抹平，误判了轨道杂化，后来固定 0.05 eV 并交叉验证。

三、关键技术要点：带隙低估与修正

PBE 算出来的带隙普遍比实验小 30%~50%，这是理论计算dos 必须坦诚的短板。我处理有两种路：一是老实用 PBE 但定性分析带边归属，不纠结绝对数值；二是上 HSE06 或 GW，带隙能拉回实验值，代价是算力涨一个量级。曾经一个二维半导体客户非要”精确带隙”，我上了 HSE06，从 0.6 eV 修到 1.4 eV，和光吸收实测对上。选哪种，取决于结论靠不靠带隙绝对值。

四、实操流程：从自洽到画图

我做理论计算dos的步骤：结构充分弛豫 → 加大 k 网格自洽（保留 WAVECAR）→ 用更密 k 做 DOS 非自洽 → 输出 total 和 lm-decomposed → 用 p4vasp 或自己脚本画 PDOS 并分色。画图时我会把费米能级平移到 0，价带导带分上下，关键峰用文字标注。最容易被省的是非自洽这步，但直接用自洽稀疏 k 画的 DOS 会出锯齿假象。

举个 MoS₂ 单层的例子：总 DOS 看是 1.8 eV 直接带隙，PDOS 拆开才发现导带底 90% 是 Mo 的 d_z²，价带顶是 S 的 p_z，说明光吸收是层内 d–p 跃迁主导，和实验偏振依赖一致。如果不做 PDOS，我可能误判为普通带边，漏掉后面做应变调控的依据。现在凡是要解释光学或输运的 DOS，我默认都出 PDOS 分色，多花十分钟但结论稳得多，也不会在答辩时被一句”这个峰是谁的”问住。

五、常见踩坑：自旋极化忘了开

磁性材料（含过渡金属、缺陷）做理论计算dos 忘了 ISPIN=2，会得到完全错误的占据和带隙。我算过一个含 Fe 的氧化物，非磁计算显示金属性，开自旋极化后是 1.2 eV 带隙，和实验一致。另一个坑是 Fermi 面穿过窄带造成假金属态，往往是结构没弛豫稳或 k 点不够，要先排除再下结论。

六、复盘：DOS 是电子结构的”指纹”

回过头看，理论计算dos 最实用的不是报一个数，而是给你一张能指认成键、带边、磁性的指纹图。我接项目时默认先问”你想从 DOS 里看什么”——是带隙、是轨道杂化、还是磁性起源，再决定算多细。把 PDOS 拆分、带隙修正、自旋三件事做到位，图才有说服力。被证明有用的，是愿意为一个峰翻文献指认的较真。下次有人问”这个峰是什么”，我会先把 PDOS 分层，而不是盯着总态密度猜。