Materials Studio计算分波态密度：CASTEP的PDOS设置与轨道贡献分析

分波态密度（Projected Density of States，PDOS）是理解材料电子结构和化学键性质最直观的工具之一。与总DOS相比，PDOS能把每个原子、每种轨道（s、p、d）对电子态的贡献单独显示出来，让成键机制一目了然——哪个原子的哪个轨道在费米能级附近有贡献，哪两个原子的轨道之间有杂化，看PDOS往往比看能带更直接。

MS/CASTEP中PDOS的计算逻辑

在Materials Studio中，PDOS的计算在CASTEP的Energy或单点任务完成后，通过”Population Analysis”功能提取。CASTEP计算DOS的底层机制与VASP类似：在自洽收敛的电子密度基础上，将本征态的贡献投影到各原子的球谐函数基上，得到各轨道分量的贡献权重。

计算PDOS的前提条件：

1. 结构必须已完成优化（几何收敛，无虚频驱动）
2. 自洽计算必须充分收敛（建议能量收敛到 1×10⁻⁶ eV/atom）
3. k点网格需要比优化计算时更密，以保证DOS的分辨率

关键参数：k点密度与DOS展宽

k点密度是影响PDOS精度的第一要素。DOS是k点本征值的加权求和，k点越密，DOS曲线越光滑、细节越丰富。通常建议DOS计算的k点密度是结构优化时的2-3倍。

在Materials Studio的CASTEP设置面板中，k点可以通过”k-point separation”来控制：优化时设0.05 Å⁻¹，DOS计算时设0.03 Å⁻¹或更小。对于金属体系，费米能级附近的态密度特征往往很尖锐（范霍夫奇点），需要更密的k点才能正确描述。

展宽方法（Smearing）：CASTEP支持Gaussian、Methfessel-Paxton和冷展宽等方法。展宽参数（Smearing Width）影响DOS曲线的平滑程度——太大会模糊掉真实的峰结构，太小在k点不够密时会出现锯齿状噪声。金属体系通常用0.1-0.2 eV展宽；绝缘体和半导体的带隙区间没有态，对展宽不敏感，0.05-0.1 eV即可。

Population Analysis的设置

CASTEP的Population Analysis（在Analysis面板中访问）会输出Mulliken布居分析，同时生成PDOS数据文件（.pdos格式）。设置步骤：

1. 在CASTEP任务设置中的”Properties”标签页，勾选”Density of States”和”Population Analysis”
2. 设置DOS能量范围（通常取-20 eV到5 eV，以费米能级为零点）和步长（0.02 eV是合理选择）
3. 如需投影到特定轨道角量子数（s/p/d/f），确认”Project DOS”选项开启

计算完成后，在”Analysis → Electronic → Density of States”面板中可以直接绘制各原子、各角量子数的PDOS曲线，并叠加显示。对于包含多种元素的体系（如过渡金属氧化物），分别显示金属d轨道和氧p轨道的PDOS，能清晰地看到两者在何处发生了杂化（PDOS峰在同一能量位置重叠）。

PDOS数据的物理解读

费米能级附近的PDOS决定了材料的导电性和磁性。金属的费米能级处PDOS非零；绝缘体和半导体的费米能级落在带隙内，PDOS为零。

轨道杂化的识别：以TiO₂为例，价带主要由O 2p轨道贡献，导带底主要是Ti 3d轨道。如果在两者的PDOS峰中有部分重叠，意味着Ti-O之间存在一定的共价键成分（p-d杂化），而不是纯离子键。这种信息在纯DOS图上看不出来，但在PDOS上一目了然。

Mulliken重叠布居（Overlap Population）是成键分析的补充工具，正值表示成键作用，负值表示反键，零值表示非键。这个信息在CASTEP的Population Analysis输出文件中以文字形式列出，需要与PDOS图配合理解。

磁性体系的自旋PDOS：对ISPIN=2的计算，CASTEP会分别输出自旋向上和向下两套PDOS。两者之差（磁矩的轨道分解）可以直接指出哪个原子的哪个轨道是磁矩的主要来源，这对理解磁性机制（如Stoner判据和超交换作用）非常有价值。

科研学术网提供Materials Studio全套电子结构计算，分波态密度、能带与成键分析均可承接，欢迎访问 https://www.keyanxueshu.com 获取报价。