在上一章《第二十一VASP计算实战：能带结构计算与绘图步骤讲解！| 2026新版VASP基础教程》中，华算科技朱老师详细介绍了VASP示例任务。VASP计算一般从结构优化开始，然后对优化后的结构进行性质计算。本章将介绍VASP结构优化计算实例，让大家更好的了解计算流程，具体包括结构优化、半导体性质案例、表面与催化案例。
态密度D(E)描述了在能量范围E到E+dE内，每单位能量范围内的电子态数目。
总态密度（TDOS）
反映材料整体的电子分布情况。它是分析材料是否为金属（EF 处D(E)≠0）还是绝缘体/半导体（EF处D（F）＝0）的直接依据。
投影态密度（PDOS）
通过将总态密度投影到特定原子、特定轨道（如s、p、d）上，揭示了哪些原子或轨道主导了导电性或光学跃迁。
费米能级（EF）
通常将费米能级设为零点（E=0），便于观察价带（Valence Band）和导带（Conduction Band）的相对位置。
VASP 计算态密度的标准流程严谨且高效，必须严格执行以下三个阶段：
目的
获得材料的基态晶体结构，消除内应力。
文件准备
INCAR（设置 IBRION=2 进行离子弛豫，NSW=60 设定最大步数），POSCAR（初始结构），POTCAR（赝势），KPOINTS（k点网格）。
关键点
确保力收敛（EDIFFG）足够小。优化结束后会生成 CONTCAR（最终结构）、CHGCAR（电荷密度）和 WAVECAR（波函数）文件。
目的
在固定几何结构下获得高精度的电子密度分布和费米能级。
文件准备
复制第一步的 CONTCAR 为新的 POSCAR，保留 CHGCAR 和 WAVECAR。
INCAR 参数
ISTART = 1 或 3（从 CHGCAR/WAVECAR 读取上一步数据）。
ICHARG = 11（从电荷密度文件读取）。
ISMEAR = 0（高斯展宽，适用于非金属）。
SIGMA = 0.02（展宽参数，通常较小）。
关键点
此步至关重要，它决定了最终 DOS 的精度。如果直接使用结构优化的结果（第一步），可能导致电子结构不收敛或费米能级不准确。
目的
基于已知电荷密度，计算密度矩阵并输出 DOS 数据。
文件准备
复制第二步的 CHGCAR 为 CHGCAR（覆盖），复制 WAVECAR。
INCAR 参数
ISTART = 1 或 3。
ICHARG = 11。
LORBIT = 11（打开投影态密度输出，生成 PROCAR）。
NEDOS = 2000（DOS 点数，点越多图越光滑）。
EMIN / EMAX（设置能量范围，如 -10 到 10 eV）。
运行：生成 DOSCAR（总态密度）和 PROCAR（投影态密度）文件。
态密度的质量在很大程度上取决于 INCAR 中的细节设置：
NEDOS：决定了能量轴的分辨率。一般设为 1000-2000 足以绘制光滑曲线。如果分析细节（如 Van Hove 奇异点），建议设置更高。
LORBIT：控制投影输出。
LORBIT = 11：输出 PROCAR，可分析原子分波态密度（s/p/d）。
LORBIT = 10：仅输出总态密度。
ISMEAR 与 SIGMA：金属材料通常设 ISMEAR = 1（高斯展宽），非金属设 0。SIGMA 越小，费米能级越精确，但计算收敛难度增加。
NBANDS：如果分析高能区的导带，需要增加能带数，确保空带数足够（通常设为原子数的 2 倍以上）。
VASP 计算结束后，最关键的工作是将 DOSCAR 和 PROCAR 文件转换为可视化图形。常用方法包括：
对于 Python 开发者，pymatgen 的 Dos 对象非常强大。
代码示例：
from pymatgen.io.vasp import Vasprun
from pymatgen.electronic_structure.plotter import DosPlotter
vr = Vasprun(‘vasprun.xml’)
dos = vr.complete_dos
plotter = DosPlotter()
plotter.add_dos(‘Total DOS’, dos)
plotter.get_plot().show()
DOSCAR 前 5 行是头信息（包括能量范围、NEDOS、是否投影）。
第 6 行起，每行数据格式为：
能量 总态密度 电子数
如果是投影态密度文件（PROCAR），则包含每个原子、每个轨道的分布。
什么是态密度：介绍了总态密度、投影态密度、费米能级
计算步骤：介绍了结构优化，自洽计算，非自洽计算
DOS数据提取与可视化：介绍了pymatgen + sumo工具与手动处理方法
下一章将正式引入本次教程的核心—吸附位点建模。我们将从模型选择、吸附位点选取与建模、输入文件方面详细介绍VASP表面与催化案例，敬请期待！
