d带中心定义为d轨道态密度（d-DOS）的一阶矩：
其中：
：d带中心能量位置（相对于费米能级）
E：能量（eV）
：d轨道态密度
：费米能级
d带中心位置 |
物理意义 |
催化活性 |
|
较高（接近E_F） |
d轨道与吸附物种轨道重叠强 |
吸附强，但可能过强 |
|
较低（远离E_F） |
d轨道与吸附物种轨道重叠弱 |
吸附弱 |
|
最佳位置 |
适中的吸附强度 |
Sabatier原理 |
金属表面催化：Pt、Pd、Rh、Ru等贵金属
合金催化：Pt-Ni、Pt-Co等合金表面
氧化物表面：过渡金属氧化物催化剂
单原子催化：金属中心原子电子结构
建立合适的表面slab模型是d带中心计算的前提：
slab模型要求：
足够的层数（通常4-6层）
底部层固定，顶部几层优化
适当的真空层（≥15 Å）
合理的面外晶格常数
示例：Pt(111)表面模型
Pt(111)
3层slab + 12Å真空
4×4超胞
底部2层固定
顶部1层允许弛豫
POSCAR文件：
Pt(111) Surface Optimization
1.0
26.76000 0.00000 -15.45800
-13.38000 23.16931 -15.45800
0.00000 0.00000 25.00000
Pt
36
Direct
0.00000 0.00000 0.00000
0.33333 0.66667 0.00000
0.66667 0.33333 0.00000
0.50000 0.50000 0.33333
…
优化INCAR关键参数：
SYSTEM = Pt(111) Optimization
PREC = Accurate
ENCUT = 400
EDIFF = 1E-4
EDIFFG = -0.02 # 力收敛标准
IBRION = 2
NSW = 50
ISIF = 2 # 固定晶胞，优化原子位置
ISYM = 0 # 关闭对称性
LWAVE = .FALSE.
LCHARG = .TRUE. # 保存电荷密度
优化完成后，检查：
力收敛（所有力）
能量收敛（能量变化）
晶胞参数合理性
原子位置合理无重叠
d带中心计算需要高精度的态密度数据：
|
参数 |
推荐值 |
说明 |
|
LORBIT |
11 |
输出轨道投影态密度 |
|
ICHARG |
1 |
从CHGCAR读取电荷密度 |
|
NEDOS |
3000-5000 |
能量点数，越高越精确 |
|
EMIN |
-10–15 |
能量下限（d带通常位于费米能级以下） |
|
EMAX |
5-10 |
能量上限 |
|
ISMEAR |
-5 |
四面体方法，适合绝缘体/半导体 |
|
ISMEAR |
0 |
Gaussian展宽，适合金属 |
|
SIGMA |
0.05-0.1 |
展宽宽度（金属体系） |
金属体系（Pt、Pd等）：
ISMEAR = 0
SIGMA = 0.1
半导体/绝缘体体系：
ISMEAR = -5
SIGMA = 0.05
# ========== 基本参数 ==========
SYSTEM = Pt(111) d-band Calculation
PREC = Accurate
ENCUT = 400
EDIFF = 1E-6
# ========== 态密度参数 ==========
ISMEAR = 0
SIGMA = 0.1
NEDOS = 4000
EMIN = -10.0
EMAX = 5.0
LORBIT = 11
ICHARG = 1
# ========== 计算控制 ==========
ISYM = 0
LWAVE = .FALSE.
LCHARG = .FALSE.
ISPIN = 2 # 如需要磁性计算
# ========== K点设置 ==========
KPPRA = 15000
LORBIT=11：
输出总态密度（DOS）
输出各原子分态密度（PDOS）
输出轨道角动量投影（s、p、d、f）
NEDOS选择：
NEDOS=1000：快速预览
NEDOS=2000-3000：常规计算
NEDOS=4000-5000：高精度d带中心
能量范围：
EMIN：一般设为费米能级以下8-10 eV
EMAX：一般设为费米能级以上2-3 eV
d带中心的物理定义：d带中心的数学表达式、d带中心的物理意义、应用范围
结构优化：表面模型构建、结构优化计算，优化后处理
INCAR参数设置详解：态密度计算参数、金属体系vs绝缘体体系、完整INCAR示例
下一章将正式引入本次教程的核心—VASP电荷密度计算。我们将从电荷密度的物理意义与重要性、VASP电荷密度计算设置、CHGCAR文件格式详解方面详细介绍VASP电荷密度计算。
