VASP 体系与性质计算全景指南
从晶体、表面、缺陷到电化学与激发态,按研究目标查找 VASP 原生功能、扩展工具、计算流程和适用边界。
阅读本篇目录
VASP 能研究的内容远不止结构优化、能带和态密度。它可以作为电子结构计算程序直接给出能量、力、应力和多种响应量,也可以与声子、输运、缺陷、合金热力学、恒电势和非绝热动力学工具连接,完成更复杂的研究。
本指南按体系 → 性质 → 方法 → 扩展 → 结果解释组织,适合选题、设计计算方案和寻找下一步工具。资料核对日期为 2026年10月24日,依据 VASP 官方手册、扩展项目文档及方法论文。文中的材料名称用于说明应用场景;实际计算应根据体系选择方法,并完成相应收敛测试。
阅读时可以直接使用页面目录跳到研究方向。准备实际提交任务时,请同时查阅版本与扩展选择、CPU 参数及GPU 参数。
1. 先区分四种实现方式#
| 标记 | 含义 | 典型例子 |
|---|---|---|
| 原生 | VASP 本身提供对应算法;仍受版本、编译选项和方法组合限制 | 能带、DFT+U、SOC、介电响应、GW、AIMD |
| 组合 | 多次 VASP 计算,再按物理定义拟合、求差、求导或统计 | 形成能、吸附能、弹性模量、扩散系数、平均嵌入电压 |
| 扩展 | 需要补丁、链接库或与另一套求解程序配合 | VTST、VASPsol、CP-VASP、Libxc、i-PI、DMFT |
| 后处理 | VASP 提供结构、力或电子数据,外部程序得到最终性质 | phonopy、LOBSTER、AMSET、WannierTools、ATAT |
“原生”不等于四个输入文件一次运行就能得到全部结果。“扩展”也不等于任意插件可以叠加:VASP 版本、补丁分支、计算方法、CPU/GPU 平台必须同时匹配。VASPKIT 等程序可以简化输入和提取结果,但不会自动补上缺少的物理模型。官方功能分类可参阅 VASP 手册。
2. 可以建立哪些体系#
VASP 通常采用周期性超胞。非周期方向通过真空及适用的静电处理隔开;“放进一个大盒子”仍需检查周期镜像相互作用。元素选择也应落实到合适的 PAW 数据集,而不是只看元素符号能否输入。PAW 选择、表面与二维体系、静电修正。
| 体系 | 常见研究问题 | 建模重点 |
|---|---|---|
| 金属、金属间化合物 | 相稳定性、费米面、磁性、弹性、输运 | 密集 k 网格、占据展宽和磁态 |
| 半导体、绝缘体 | 带隙、缺陷、介电、声子、光学 | 泛函误差、带边位置及激子效应 |
| 氧化物、氮化物、碳化物、陶瓷 | 结构、强度、极化、离子迁移 | 多价态、磁序、局域电子和竞争相 |
| 过渡金属及稀土化合物 | 轨道占据、关联、SOC、磁各向异性 | d/f 电子处理、PAW 价电子及 U/J 的选择 |
| 合金、高熵材料、固溶体 | 混合焓、短程有序、偏析、相图 | 多构型采样、SQS、簇展开,不能只算一种排布 |
| 二维材料、单层与多层 | 面内力学、剥离能、带隙、拓扑 | 真空、层间色散、面量单位、库仑镜像 |
| 纳米线、纳米带、纳米管 | 边缘态、应变、吸附、电子结构 | 横向真空、端面/边缘钝化、一维采样 |
| 分子、团簇、纳米颗粒 | 构型、结合能、振动、电子亲和能 | 自旋多重度、盒子尺寸、电荷修正及异构体 |
| 表面、薄膜、重构表面 | 表面能、功函数、STM、吸附与催化 | slab 厚度、终止面、固定层、偶极与覆盖度 |
| 异质结、界面、晶界 | 带阶、黏附、界面电荷、偏析 | 晶格失配、应变分配、界面模型与厚度 |
| 空位、间隙、替位、反位、缺陷复合体 | 形成能、电荷态、迁移、俘获 | 超胞尺寸、化学势、有限尺寸修正与局域畸变 |
| 位错、层错、裂纹附近局域结构 | 缺陷核心、滑移能面、断键机制 | 长程弹性场和周期边界;大尺度常需多尺度方法 |
| 非晶、玻璃、液体、熔体 | 局域结构、扩散、密度、反应 | 制样历史、时间长度、多初态及统计误差 |
| 溶液、固液界面、电极 | 溶剂化、质子转移、双电层、恒电势 | 显式/隐式溶剂、离子浓度、电势参照和采样 |
| MOF、COF、沸石、多孔材料 | 吸附位点、结合能、振动与反应 | 色散、开放金属位点、孔道与客体构型 |
| 有机晶体、聚合物晶区、杂化钙钛矿 | 晶体堆积、氢键、光电和热性质 | 色散、构象、动态无序;无序聚合物需更大采样 |
| 高压材料、地球深部矿物、氢化物 | 状态方程、相变、振动和电子结构 | 半芯态、硬势、高截断及压力范围有效性 |
| 磁性、铁电、多铁、拓扑候选材料 | 序参量、耦合、拓扑不变量 | 对称性、SOC、磁序与合适的响应算法 |
小型生物分子片段、表面吸附的有机分子也可以建模;完整蛋白质、长时间聚合物动力学等问题通常更适合专用量子化学、经典力场或 QM/MM 工作流。巨型莫尔超胞、低浓度缺陷和极慢过程的主要困难往往是规模与采样,而不是材料名称不被识别。
3. 结构、能量与相稳定性#
| 性质或问题 | 实现方式与典型流程 | 需要怎样理解结果 |
|---|---|---|
| 平衡晶格、键长、键角、原子位置 | 原生:原子/晶胞优化,IBRION、ISIF |
优化找到附近极小值,不保证全局最低能结构 |
| 固定晶格、外延应变、部分原子弛豫 | 原生/扩展:Selective dynamics、约束晶胞;特定约束可用 optcell 分支 | 保持约束与真实实验边界一致 |
| 能量、力、应力、残余压力 | 原生:SCF 与力/应力输出 | Pulay 应力和电子收敛会影响结构与力学结果 |
| 内聚能、结合能、反应能 | 组合:体系与统一参照态能量求差 | 原子自旋、分子参照和正负号约定必须说明 |
| 形成能、混合焓、分解能 | 组合:相对于元素或竞争相计算 | 相对元素稳定不等于相对所有化合物稳定 |
| 凸包、分解产物、化学势稳定区 | 后处理:相容 DFT 数据加相图分析 | 数据库泛函、U、PAW 与能量修正需一致 |
| 状态方程、体积模量、压力导数 | 组合:多体积 E–V 数据拟合 | 不同拟合区间与 EOS 模型有差异 |
| 高压相变、焓排序 | 组合:各相随压力的 E+pV 比较 | 有限温度需进一步加入自由能 |
| 表面能、解理能、剥离能 | 组合:bulk、slab、分离层能量差 | 明确单面/双面、面积、终止与化学计量 |
| 界面黏附功、晶界能、偏析能 | 组合:界面与对应参照构型求差 | 界面数、应变参照及原子数必须配平 |
| 结构预测、未知相搜索 | 外部搜索器调用 VASP:USPEX、CALYPSO 等 | VASP 负责评价候选;搜索算法负责探索构型 |
实际开始时先用小规模收敛测试确定截断能、k 网格和力阈值,再批量生成相互可比的数据。可参阅 ISIF、VASPKIT 力学与 EOS 工具、USPEX 手册及 CALYPSO 手册。
4. 弹性、力学与应变工程#
| 性质 | 路线 | 关键边界 |
|---|---|---|
| 二阶弹性常数 Cij | 原生/组合:有限应变,读取应力或拟合能量 | 区分固定离子与内部弛豫后的弹性常数 |
| 体积、剪切、杨氏模量与泊松比 | 后处理:弹性张量、Voigt/Reuss/Hill 等平均 | 单晶方向性质与多晶平均不同 |
| 弹性各向异性、方向压缩率 | 后处理:张量旋转和柔度矩阵 | 不能用单一标量概括所有方向 |
| 力学稳定性 | 后处理:适用晶系与压力条件下的稳定性判据 | 不替代声子动力学稳定性或热力学稳定性 |
| 理想拉伸/剪切强度、应力–应变曲线 | 组合:逐步加载并弛豫 | 是指定缺陷和温度模型的理想响应,不是工程屈服强度 |
| 二维弹性模量 | 组合:应变计算后按面量归一 | 通常用 N/m;不可直接比较含不同真空厚度的 GPa |
| 声速、弹性 Debye 温度 | 后处理:弹性张量与密度、连续介质近似 | 属于模型推导,不能代替完整声子色散 |
| 层错能、广义层错能面、滑移势垒 | 组合:面内相对位移与受限弛豫 | 正确处理周期界面与位移方向 |
| 应变对带隙、磁性、吸附或极化的影响 | 组合:一系列应变结构重复性质计算 | 防止同时改变参照、磁态或晶相 |
| 硬度、脆韧趋势、断裂相关指标 | 后处理/多尺度:弹性或键合经验模型、断裂模型 | 经验硬度不等于直接计算压痕;真实韧性还涉及裂纹与塑性 |
VASP 的应变响应与静态电响应存在关联,但不能把不同张量或不同机械边界混用。线性响应总览、VASPKIT 力学功能。
5. 能带、态密度与电子分布#
| 性质 | 路线与输出 | 使用要点 |
|---|---|---|
| 能带、直接/间接带隙 | 原生:SCF 后沿路径计算,EIGENVAL/HDF5 |
高对称路径可能漏掉真正带边,需要额外网格搜索 |
| 总 DOS、原子/轨道/自旋 PDOS | 原生:均匀 k 网格、DOSCAR、投影输出 |
DOS 不能使用稀疏高对称线积分 |
| 投影能带、轨道成分、自旋纹理 | 原生+后处理:PROCAR 等 |
投影依赖球半径/投影定义,SOC 需要相应自旋处理 |
| 有效质量、带边曲率与各向异性 | 组合/后处理:带边附近细网格拟合 | 拟合范围、简并和非抛物线性影响数值 |
| 费米面、费米速度、电子口袋 | 后处理:致密三维采样或 Wannier 插值 | 不能只从一张路径能带图判断费米面 |
| 超胞能带展开、缺陷/合金谱权重 | 后处理:VASPKIT、easyunfold 等适用工具 | 展开结果是谱权重,不是新的原胞本征值 |
| 分子轨道、局域态、带分解电荷 | 原生+可视化:WAVECAR、LPARD、PARCHG |
当前官方 LPARD 路径不支持非共线计算 |
| 电荷密度、差分电荷、自旋密度 | 原生/组合:CHGCAR 及同网格差分 |
差分各分量应保持同一盒子和原子坐标定义 |
| ELF、电子局域化 | 原生+可视化:ELFCAR |
ELF、电子密度和电荷转移是不同量 |
| 功函数、真空能级与表面偶极 | 组合:LOCPOT/HDF5 平面平均加费米能 |
是指定表面的性质;真空区需出现可信平台 |
| 带阶、界面势差、电子亲和势 | 组合:共同电势参照或显式界面、带边修正 | 孤立材料的真空对齐不一定等于实际界面带阶 |
| 电离能、电子亲和能、充电能 | 组合:适用的总能差/准粒子方法 | 带电超胞需修正;KS 的 HOMO/LUMO 差不等于全部实验量 |
推荐先明确“需要一张图还是需要定量带边”。带隙、光学峰和激子能是不同目标,后两者见后文。显存不足时可使用分段能带工作流,但分段路径不能替代输运需要的完整均匀网格。
6. 电荷分析、化学键与实空间拓扑#
| 目标 | 推荐路线 | 结果解释 |
|---|---|---|
| Bader 电荷、原子盆积分 | VASP 密度加 Bader/critic2,注意 PAW 重构密度 | 是一种空间划分结果,不是唯一“真实原子电荷” |
| Hirshfeld、DDEC 等原子电荷 | 兼容密度输出加对应独立分析程序 | 与 Bader 定义不同,不直接混合比较 |
| COHP、ICOHP、COOP、COBI | VASP 波函数加 LOBSTER | 检查局域基组投影质量;成键能指标不等于孤立键解离能 |
| QTAIM 临界点、键径、电子密度拓扑 | VASP 密度加 critic2 | 网格精度和核心区域密度处理影响结论 |
| 非共价相互作用图、氢键和弱作用区 | 兼容密度加 critic2 等 | 图形是分析线索,结合能仍需明确参照计算 |
| d 带中心、带宽、轨道占据 | PDOS 积分与 VASPKIT 等分析 | 积分范围和投影方案会改变指标,不宜当通用催化活性定律 |
这类工具通常无需重编译 VASP,但对输出文件和计算设置有要求。先读工具说明再做静态计算,避免计算完才发现缺少波函数。LOBSTER 方法论文、LOBSTER 键合分析综述、critic2 功能说明。
7. 声子、振动与平衡热力学#
| 性质 | 路线 | 注意事项 |
|---|---|---|
| Γ 点振动频率、本征矢 | 原生:有限位移或适用的 DFPT | 分子平移/转动零模与真正虚频要区分 |
| 声子色散、总/分波声子 DOS | 原生力常数工作流或 VASP+phonopy | 原胞 Γ 点频率不能代替全布里渊区稳定性 |
| 动力学稳定性、软模 | 声子色散与本征矢分析 | 小虚频需查精度、尺寸和声学和规则;软模也可能反映相变 |
| LO–TO 分裂 | Born 有效电荷、介电张量加非解析修正 | 极性材料的 Γ 附近具有方向依赖 |
| 零点能、振动自由能、熵、Cv | 后处理:频率/声子 DOS 统计 | 谐近似有效范围要明确 |
| 热膨胀、Cp、温度相关体积模量 | 组合:多体积声子+准谐近似 QHA | 强非谐、熔化附近不能盲目外推 |
| 模 Grüneisen 参数、声子群速度 | 后处理:体积变化与色散导数 | 需要相同模式/晶胞定义与充分收敛 |
| 原子均方位移、热椭球 | 后处理:phonopy 等 | 可用于解释晶体热运动,与扩散 MSD 不同 |
| 气体分子热化学修正 | 频率加平动/转动/振动统计 | 温度、标准压力、对称数和低频处理需说明 |
| 吸附自由能中的振动修正 | 表面吸附体频率与参照态热化学 | 受阻转动/平动与极低频模式可能不宜简单谐振处理 |
| 中子/非弹性 X 射线散射相关 S(q,ω) | phonopy 动态结构因子等后处理 | 需要散射因子、实验几何和适用近似 |
基本流程是优化 → 位移/扰动 → VASP 力或响应 → 力常数 → 色散/热力学。不要把电子占据展宽温度直接当成晶格温度。VASP 声子、phonopy 功能与 VASP 接口、phonopy QHA。
8. 非谐振动与晶格热输运#
| 目标 | 实现方式 | 额外工作 |
|---|---|---|
| 三阶力常数、声子散射与寿命 | VASP 位移受力+phono3py/thirdorder.py | 大量构型、位移幅度与相互作用截断收敛 |
| 晶格热导率张量 | phono3py 或 ShengBTE 等求解声子 BTE | 区分 RTA 与迭代/直接求解;核对 q 网格 |
| 模分辨热导、平均自由程累积 | 上述工具后处理 | 是所包含散射机制下的分解 |
| 同位素、边界等附加散射 | 按求解器支持的模型加入 | 理想晶体结果不自动包含真实晶界、孔洞与缺陷 |
| 温度重整化声子、高温稳定相 | VASP 力样本+TDEP、ALAMODE、自洽声子等 | 采样温度、拟合阶数和统计质量重要 |
| 高阶非谐、四声子效应 | 高阶力常数加支持该方法的求解器 | 不能把普通三阶工作流称为已包含四声子 |
| 非谐自由能和相变温度 | 有限温度势能面、热力学积分或相应晶格动力学 | 需要各相一致自由能参照 |
| MD 热输运 | VASP 的 Müller–Plathe 路线,或合适势模型下热流相关方法 | 时间、长度、温度梯度及经典核近似需收敛 |
有声子谱不等于有热导率。 二阶力常数给出谐振动;定量热导还需要散射/非谐信息。声子 BTE 与 MD 是不同路线,不能直接把其输出拼在一起。phono3py、ShengBTE、TDEP、ALAMODE、Müller–Plathe 方法。
9. 极化、介电、压电与铁电#
| 性质 | 路线 | 关键边界 |
|---|---|---|
| Berry 相位极化、自发极化变化 | 原生 LCALCPOL,沿绝缘路径跟踪分支 |
极化具有量子多值性;不能直接减任意两个输出 |
| 电子介电张量 | 原生 LEPSILON 或适用有限电场方法 |
区分钳制离子与包含离子贡献 |
| 离子介电贡献、静态总介电 | Born 电荷、振动模与电子介电组合 | 软模、温度和非谐影响可能很大 |
| Born 有效电荷 | DFPT/有限电场或位移导数 | 是极化响应张量,不是 Bader 静态电荷 |
| 压电 e、d 张量 | 电响应、离子弛豫与弹性柔度组合 | e 与 d 单位不同;固定离子与弛豫离子要分开 |
| 铁电双势阱、极化翻转路径与能垒 | 受限畸变路径、NEB 和极化跟踪 | 单胞均匀翻转不等于真实畴成核与实验矫顽场 |
| 外场诱导结构/能带变化 | 适用的有限电场或 slab 外场工作流 | 周期体相与真空 slab 的电场设置不同 |
| 热释电、挠曲电、电致伸缩、磁电耦合 | 专门的温度/应变梯度/电磁响应组合研究 | 需各自理论与数值实现,不是统一的一键原生输出 |
有限电场、DFPT 和杂化泛函的可组合性应逐项查手册,不能把普通 GGA 输入直接套到所有方法。极化、LEPSILON、介电响应、slab 外场。
10. 磁性、自旋与磁激发#
| 目标 | 路线 | 解释与限制 |
|---|---|---|
| 自旋极化、总/局域磁矩 | 原生:共线自旋与磁化密度 | 初始磁矩可能收敛到不同亚稳态;局域积分依赖定义 |
| 铁磁、反铁磁、亚铁磁能量比较 | 组合:匹配的多个磁构型 | 需要容纳正确磁周期的超胞 |
| 非共线磁结构、SOC 能带 | 原生:非共线/SOC,匹配 vasp_ncl |
SOC、对称性与初态必须一致 |
| 磁各向异性能、易磁化轴 | SOC 下不同磁化方向的能量比较 | 通常能差很小,对 k 网格与收敛很敏感 |
| 约束磁矩、固定方向/大小研究 | 原生约束磁矩方法 | 罚函数项与真实物理能量需正确区分 |
| 自旋螺旋 E(q)、冻结磁振子 | 原生广义 Bloch 自旋螺旋 | 该路线不能与 SOC 同时使用;需关闭不适用对称性 |
| 交换参数 J、DMI、磁哈密顿量 | 多构型能量映射或适用的外部线性响应工具 | 映射模型、截断、磁矩刚性和 SOC 处理需验证 |
| 磁振子色散、居里/奈尔温度 | 上述参数加自旋波/Monte Carlo/自旋动力学 | 不由一次自旋极化 SCF 直接给出 |
| 轨道磁化、磁光和异常霍尔 | Berry/Wannier 方法与光学响应 | 需要相应矩阵元、SOC 和细密积分 |
| 磁致伸缩、多铁耦合 | 磁化方向、应变、极化的组合计算 | 区分本征单畴响应与实验多畴响应 |
参考:磁性功能、自旋螺旋理论与限制、WannierBerri。NMR 的轨道磁化率也不是包含所有磁性贡献的实验总磁化率。
11. 光学、激子与多体电子结构#
| 性质 | 路线 | 需要增加什么 |
|---|---|---|
| 复介电函数、光学各向异性 | 原生 LOPTICS 独立粒子近似 |
足够空带、k 点、频率范围和适当展宽 |
| 吸收系数、折射率、消光系数、反射率 | 从介电函数后处理 | 体材料公式与薄膜光学模型要区分 |
| 电子能量损失函数、等离激元特征 | 介电响应及适用的有限 q 方法 | 光学极限损失函数不等于所有动量转移 EELS |
| 光学电导、磁光 Kerr 响应 | 复介电张量、SOC 与相应后处理 | 金属低频需考虑带内贡献和散射模型 |
| 准粒子带隙、电子谱、自能修正 | 原生 GW 路线 | 空带、响应截断、频率积分与起始态收敛 |
| 激子能、激子吸收峰、束缚能 | 原生 BSE,通常结合 GW/适当能带 | 束缚能需用一致准粒子隙与激子能,区分明/暗态 |
| TDHF/TDDFT 激发与光谱 | 原生 Casida/响应方法 | 泛函核和筛选决定近似;局域核不保证正确激子 |
| 时间传播求光学响应 | 原生 ALGO=TIMEEV 的响应传播 |
不应据此声称具备任意强激光非平衡动力学 |
| RPA/ACFDT 相关能与相对能 | 原生多体相关能路线 | 是总能改进路线,与 GW 准粒子能用途不同 |
| MP2 类相关能 | 原生 MP2/LTMP2 等实现 | 规模和金属性限制严格;不等于完整分子量化软件功能集 |
| 有效相互作用 U、J | 原生 cRPA;也可选适用线性响应 U 工作流 | 所得参数依赖目标子空间与屏蔽定义 |
二维材料的介电响应随超胞真空厚度变化,应明确二维极化率/片电导等归一约定。吸收峰还可能受温度、激子、无序和实验展宽影响;仅对 GGA 能带做剪刀修正不能自动补齐这些物理。
参考:LOPTICS、多体方法总览、TDDFT 计算、时间传播算法。
12. 红外、拉曼、X 射线、NMR 与扫描探针#
| 实验对应量 | 计算路线 | 不能省略的处理 |
|---|---|---|
| 红外频率与振子强度 | 振动本征矢、Born 电荷与介电响应 | 频率不等于强度;实验峰形还需展宽/环境 |
| 非共振拉曼活性、拉曼张量 | 沿正常模位移计算介电导数,vasp_raman.py 等 |
Γ 点声子本身不会给出完整拉曼强度;共振拉曼需其他理论 |
| XANES/NEXAFS、近边吸收 | 原生超胞核空穴 SCH 或核激发 BSE | 元素/吸收边、空带和核心空穴模型需匹配 |
| XPS 核能级位移 | 初态核能级或终态核空穴总能差等路线 | 相对位移与绝对结合能不同,需要参照及电荷处理 |
| EELS/ELNES | 相应介电损失或核激发近边方法 | 区分价电子与核心激发、q 范围和实验几何 |
| NMR 屏蔽张量、化学位移 | LCHIMAG 的 GIPAW 响应;用参照化合物换算 |
官方该路线适用于有带隙的非金属体系 |
| NICS、诱导电流、轨道磁化率 | LCHIMAG 配合相关输出/空间采样 |
不等于任意金属 Knight shift 或完整 EPR 谱 |
| 超精细耦合 | 原生 LHYPERFINE |
需要自旋密度、核信息与正确单位 |
| 电场梯度 EFG、核四极耦合 | 原生 LEFG,结合核四极矩 |
EFG 本身与核同位素参数分开;重元素近似需注意 |
| 穆斯堡尔谱相关局域信息 | EFG、磁超精细等加专门标定/谱模型 | 不能把 VASP 输出直接称为完整实验谱 |
| STM/STS 近似图像、局域态分布 | 带/能窗分解电荷与 Tersoff–Hamann 等模型 | 偏压、针尖与隧穿近似决定可比范围 |
| 粉末 XRD、结构因子与 PDF/RDF | 优化结构/轨迹加衍射或统计后处理 | 衍射峰还受温度、尺寸、织构和仪器影响 |
NMR 版本提醒:LNMRSHIELD 从 VASP 6.6.0 开始控制输出屏蔽还是带负号的位移约定;跨版本比较先核对符号和参照。实验化学位移仍需参照屏蔽,不能只换一个负号。LCHIMAG、LNMRSHIELD。
**XAS 范围提醒:**VASP 的主要适用范围是近边结构;高能 EXAFS 涉及大量空态与多重散射,不能把“支持 XAS”理解成所有吸收边和能区都同样方便。XAS 官方说明。
其他入口:拉曼脚本、核能级与核空穴、超精细、EFG、critic2 的 STM 与衍射功能。
13. 电子输运、热电与电子–声子耦合#
| 性质 | 路线 | 必须区分 |
|---|---|---|
| Seebeck 系数、电子电导与电子热导 | VASP 均匀能带+BoltzTraP2 等 Boltzmann 后处理 | 常数弛豫时间近似下常得到 σ/τ、κe/τ,不能直接写成绝对值 |
| 载流子浓度、化学势随温度变化 | DOS/能带积分与电荷中性条件 | 刚性带掺杂与实际缺陷掺杂不是同一模型 |
| 有散射模型的迁移率与电导 | AMSET:能带加弹性、介电、形变势等 | 声学、杂质、极性光学、压电散射各有适用条件 |
| 电子–声子矩阵元、散射率、寿命 | 新版 VASP 扰动路线及相应后处理 | HDF5、超胞势导数与密集采样是实际需求 |
| 零点带隙重整化、带隙温度依赖 | VASP AHC 路线或随机位移采样 | 区分电子–声子贡献与热膨胀贡献 |
| 电子谱函数、声子引起的线宽 | 电子自能分析 | 结果只包含所用近似和散射机制 |
| 考虑电子–声子散射的输运 | VASP 对应 Boltzmann 工作流 | 需收敛电子/声子网格、插值与极性长程项 |
| 功率因子、热电优值 ZT | 组合:S、σ、κe 与晶格 κl | 只算电子部分不能得到完整 ZT;参数需同温度同掺杂 |
| Hall、磁阻等输运 | 支持该目标的 Wannier/Boltzmann 外部工具 | 普通、异常、自旋 Hall 是不同机制 |
VASP 的扰动电子–声子工作流从 6.5.0 起提供,带隙重整化路线要求 HDF5。官方注明 6.5.0/6.5.1 的 ISPIN=2 相关问题在 6.6.0 修复,新任务应选已包含修复且支持所需泛函的版本。带隙重整化与版本限制。
这里的“电子–声子耦合”与磁模型交换、声子–声子非谐是三类不同问题。仅从形变势估计的迁移率,也不能写成包含所有散射的严格输运结果。电子–声子理论、电子–声子功能、AMSET、BoltzTraP2。
14. 拓扑、Wannier 与非线性响应#
| 目标 | 主要工具链 | 检查重点 |
|---|---|---|
| 最大局域 Wannier 函数、紧束缚模型 | VASP 接口→Wannier90 | 投影、能窗、解纠缠与 DFT 能带复现 |
| Berry 相位、Berry 曲率、轨道磁矩 | Wannier90/postw90、WannierBerri、适用的 VASPBERRY 路线 | 不同量需要的矩阵元不完全相同 |
| Chern 数、Z2、Wilson loop、Wannier 中心 | WannierTools、Z2Pack 等对应工作流 | 检查能隙、占据子空间和采样收敛 |
| Weyl 节点、手性、节点线 | Wannier 模型与节点搜索 | 需要验证真实简并与拓扑保护条件 |
| 表面态、Fermi arc、纳米带边缘态 | WannierTools 的 slab/表面 Green 函数等 | 终止面和紧束缚模型影响结果 |
| 自旋纹理、异常 Hall、自旋 Hall | Wannier/Berry 后处理 | 自旋流定义、SOC 与额外矩阵元需要明确 |
| 不可约表示、对称性标签 | VASP 波函数+irvsp 等 | 磁性、SOC 和空间群约定必须对应 |
| 位移电流、注入电流、非线性 Hall | 支持相应响应的 Wannier 后处理,例如 WannierBerri 的 ShiftCurrent | 不能只给能带或 HR 文件就假定所有矩阵元齐全 |
| 二次谐波、三阶非线性、强场高次谐波 | 专门响应/实时动力学实现,逐项核实 VASP 接口 | 普通 LOPTICS 给线性响应,不能直接冒充 SHG/HHG |
先让 Wannier 插值在目标能窗内重现 DFT 能带,再计算精细 Berry 响应。仅“插值曲线看起来平滑”不足以证明拓扑正确。WannierTools 功能、WannierBerri 能力、响应计算器、irvsp 方法论文。
15. 缺陷、掺杂、色心与载流子复合#
| 性质 | 路线 | 关键输入 |
|---|---|---|
| 中性/带电缺陷形成能 | bulk 与缺陷能量差、化学势、电荷项与修正 | 稳定相化学势范围、费米能参照与介电信息 |
| 热力学电荷转变能级 | 各电荷态弛豫后形成能交点 | 不等于 DOS 中一条缺陷本征能级 |
| 缺陷平衡浓度、补偿与可掺杂性 | doped 等热力学后处理 | 生长/退火条件、电荷中性、冻结缺陷与简并度 |
| 缺陷复合、结合与偏聚 | 多缺陷构型的相对能 | 浓度、分离距离与电荷态一致性 |
| 缺陷迁移、离子扩散能垒 | NEB/VTST,必要时多步网络 | 单条最低能路径不等于宏观扩散系数 |
| 小极化子、自陷态 | 合适泛函、畸变初态与局域占据 | 防止电子重新离域;可能存在多个局域极小值 |
| 光学吸收/发射能、零声子线 | 适用 ΔSCF/约束占据,基态与激发态势能面 | 热力学、垂直跃迁与零声子能需要分开 |
| Franck–Condon 位移、Huang–Rhys 因子、声子边带 | 构型坐标/振动模后处理 | 单模近似与多模描述精度不同 |
| 非辐射俘获系数、复合速率 | VASP 数据+Nonrad/CarrierCapture 等 | 需电子–声子矩阵元及势能面,不能仅凭缺陷深浅定寿命 |
| 色心自旋相关信息 | SOC、磁性、超精细与专门自旋哈密顿量分析 | 超精细可计算不代表所有 g 张量/零场分裂都原生可得 |
可按生成缺陷 → 搜索畸变 → 多电荷态优化 → 修正与化学势 → 能级/浓度 → 光学或动力学组织。高对称初态有时会困住结构,ShakeNBreak 等工具可帮助寻找畸变后的低能构型。doped、ΔSCF、固定占据、Nonrad。
16. 表面催化、电化学与电池#
| 问题 | 路线 | 必须指定的条件 |
|---|---|---|
| 吸附位点、吸附能、覆盖度效应 | 多位点/取向/覆盖度优化与能量比较 | 气相或溶液参照、slab 收敛与色散 |
| 解离、脱附、表面反应势垒 | NEB/dimer 与过渡态频率 | 电子能垒和自由能垒不同 |
| HER、OER、ORR、CO₂RR、NRR 自由能图 | DFT+热化学+适用 CHE 等电化学模型 | pH、电位标尺、质子/电子计数及溶剂 |
| 限制电势、热力学过电位 | 反应各基元步自由能比较 | 是热力学指标,不能直接当实验起始电势或电流 |
| 溶剂化能、溶液中结构与吸附 | 显式溶剂或 VASPsol/VASPsol++ | 连续介质不自动给出氢键网络与特异性吸附 |
| 恒电势结构与反应、恒电势 MD | 匹配版本的 CP-VASP/适用 VASPsol++ 路线 | 固定电子数与恒电势的热力学势不同 |
| 双电层、电极电荷、零电荷电势、电容 | 电荷–电势扫描及适当界面模型 | 电解质、面积、参照和电势窗口影响显著 |
| 活性、选择性、覆盖度与周转频率 | DFT 自由能/势垒+微观动力学模型 | 需要反应网络、速率常数、传质等假设 |
| 嵌入/脱出电压、储离子容量 | 不同组分相的自由能/能量与计量 | 平均电压近似、相分离与电子数需正确 |
| 电极体积膨胀、迁移与结构退化 | 多浓度结构、NEB、MD/缺陷工作流 | 原子尺度结果不直接给电池循环寿命 |
| 固态电解质离子电导 | AIMD 扩散、相关输运或迁移网络 | 自扩散估算与集体电导之间可能有相关修正 |
| SEI/CEI、腐蚀与电化学稳定窗口 | 分解反应、界面结构、自由能与动力学组合 | 稳定性窗口依赖相、环境和动力学钝化 |
VASPsol 提供隐式连续介质路线;VASPsol++ 扩展到非线性/非局域电解质等处理。CP-VASP 通过调节电子数研究给定电极电势下的结构与动力学,并提供相应溶剂界面处理。它们的功能与版本必须分别选择,不能把普通 VASPsol 分支写成自动支持 CP-VASP 的全部算法。VASPsol、VASPsol++、CP-VASP。
做自由能图时,先写清完整反应和参照状态,再决定需要真空、隐式溶剂、显式水层还是恒电势。统一的电位标尺与热力学定义比套用一套通用 INCAR 更重要。方法入口:计算氢电极与氧还原自由能、第一性原理嵌锂电压、CatMAP 微观动力学。
17. 过渡态、反应路径与速率#
| 目标 | 路线 | 如何确认物理意义 |
|---|---|---|
| 已知初末态的最小能量路径 | VASP NEB;VTST 提供改进路径与优化工具 | 插值路径、图像数、弹簧与原子对应关系 |
| 爬坡图像与更精确鞍点 | VTST CI-NEB 等对应实现 | 需要真正收敛到鞍点,最高图像不一定就是过渡态 |
| 只知道初态时寻找鞍点 | dimer、Lanczos 等适用方法 | 初始方向影响找到的机制,需探索多个候选 |
| 过渡态确认、反应坐标 | 振动分析、沿不稳定方向检查连接 | 一级鞍点应有对应反应方向的不稳定模 |
| 尝试频率、谐过渡态理论速率 | 初态/鞍点振动+TST 后处理 | 低频、隧穿、再穿越和溶剂摩擦可能需要额外方法 |
| 多机制扩散与反应网络 | 多路径势垒+动力学 Monte Carlo/微观动力学 | 单一势垒不能概括所有温度和覆盖度 |
NEB 得到的是所选势能面上的路径;有限温度溶液反应常需要自由能采样。 VTST 的源码补丁与 VTST 后处理脚本是不同部分:装了脚本不代表正在运行的 VASP 已包含补丁。VTST 官方工具说明、脚本与输出。
18. AIMD、液体、自由能与核量子效应#
| 性质或过程 | 路线 | 采样要求 |
|---|---|---|
| 有限温度结构、密度、相行为 | 原生 AIMD,适合目标的 NVE/NVT/NpT | 正确温控/压控、平衡时间、尺寸与时间步 |
| RDF、配位数、键角分布、氢键结构 | 轨迹后处理 | 多时间块统计,避免把少数快照当平均性质 |
| 自扩散系数、迁移活化能 | MSD/速度相关与多个温度 | 正确解包周期轨迹、扩散时间区间和有限尺寸 |
| 离子电导、集体迁移 | 电荷流/位移相关或适用近似 | 集体相关不能总被 Nernst–Einstein 近似忽略 |
| 黏度、振动态密度、时间相关函数 | 压力/速度等轨迹统计 | 噪声与采样长度常比单步电子精度更关键 |
| 液体反应、溶剂重排、质子转移 | AIMD 或增强采样 | 看不到反应不说明不存在,可能只是时间太短 |
| 势均力、自由能面、自由能垒 | 约束 MD、Blue Moon、metadynamics、slow growth 等 | 选择集体变量、检查正反路径与采样重叠 |
| 相间自由能差、化学势 | 热力学积分、适当参考体系 | 积分路径、参考自由能和统计误差都要说明 |
| 熔点 | 固液共存、interface pinning 或自由能比较 | 单相升温出现熔化可能存在明显过热 |
| 量子核效应、同位素效应 | i-PI+匹配 VASP 客户端的路径积分 | 普通经典核 AIMD 不包含零点核离域;成本随珠数增加 |
VASP 自带若干增强采样方法;如果采用 PLUMED 等外部驱动,必须另行核对 VASP 接口、补丁与目标方法,不能因为集群有 PLUMED 就认定任一 VASP 分支已连接。i-PI 官方也明确要求匹配的 VASP 客户端补丁。VASP MD、统计性质、热力学积分、i-PI 文档。
19. 机器学习势、合金统计与多尺度连接#
| 方向 | VASP 的角色 | 后续可研究性质 |
|---|---|---|
| VASP 原生 MLFF | 在线获取 DFT 能量、力、应力并训练/再拟合 | 更长轨迹的结构、扩散、热力学及适用热输运 |
| DeePMD、其他机器学习原子势 | 提供统一精度的标注数据 | 大体系缺陷、相变、界面与长时间动力学 |
| 主动学习/并发学习 | 在模型不确定的构型上补充 DFT 标注 | 扩大指定温压/成分范围内的可靠采样 |
| SQS 无序近似 | VASP 评价代表性无序结构 | 混合焓、体积、弹性与电子结构趋势 |
| ATAT/CASM/icet 类簇展开 | 提供构型能量训练集 | 有序–无序转变、短程有序、相平衡 |
| 合金 Monte Carlo 与自由能 | 用簇展开替代每一步 DFT | 成分/温度相图,必要时加入振动/电子/磁自由能 |
| CALPHAD 数据生成 | 提供相容的第一性原理热力学数据 | 多尺度相图与热力学数据库拟合 |
| 连续介质、相场、器件模型 | 提供弹性、界面能、迁移或电子参数 | 微结构演化和器件行为;最终求解发生在其他程序 |
机器学习势通常预测势能面,不会因为训练数据来自 VASP 就自动输出能带、激子或任意电荷/磁态。必须检查训练覆盖、独立测试集和外推构型。单个 SQS 也不能自动给出完整相图。VASP MLFF、ATAT 功能。本网站已有ATAT 专题,包含 SQS、VASP 连接、簇展开和 Monte Carlo 的分步教程。
20. 更高阶关联、超导与超快过程#
这些方向需要单独设计完整工具链,适合作为专题研究,不能承诺普通 VASP 模块一键完成。
| 目标 | 已有研究路线 | 当前使用边界 |
|---|---|---|
| 强关联谱、质量重整化、有限温度关联 | VASP 投影/局域轨道→TRIQS DFTTools→DMFT 求解 | 还需杂质求解器、双计数、U/J 与统计收敛;不是 DFT+U 的别名 |
| 高精度波函数相关能 | VASP 数据接口→Cc4s 等耦合簇工具 | 专用接口、方法和尺寸限制;不代表标准 VASP 自带通用 CCSD(T) |
| 声子介导超导的 λ、α²F、能隙、Tc | 完整电子–声子数据加适用 Eliashberg/超导工具链 | 需确认正式可用的接口;不能从声子 DOS 或普通 DOS 直接推得 Tc |
| 激发载流子弛豫、转移、非绝热复合 | 电子轨迹/波函数→Hefei-NAMD 等 | 依赖非绝热耦合、表面跳跃/退相干模型和统计;不是普通 AIMD |
| 强激光、实时电子响应、HHG | 具有相应能力的专门实时 TDDFT/研究实现 | 逐项核对代码与接口,不把线性响应 TIMEEV 当作通用强场程序 |
| 器件 I–V、接触输运、NEGF | 将适用的电子结构模型传给专门输运求解器 | 电极、自能、开放边界与偏压自洽需另建;VASP 能带不等于 I–V |
**VASP–EPW 特别说明:**已有接口论文演示 MgB₂ 的各向异性超导能隙/Tc 及 c-BN 迁移率;EPW 文档也记录了接口开发。但是 VASP 开发者在 2026年10月9日的答复仍说明该接口尚未进入 VASP 正式发布版本。因此本指南将其列为已有研究实现、需要核对配套版本的路线,不列作本集群现成可用功能。这也不应被概括成“VASP 永远不能研究超导”。接口论文、EPW 项目说明、VASP 官方近期答复。
其他参考:TRIQS 的 VASP 接口、Cc4s 说明、Cc4s–VASP 接口、Hefei-NAMD、Hefei-NAMD 开发项目。
21. 泛函与常用扩展怎样改变计算能力#
| 方法/扩展 | 解决什么问题 | 选择建议 |
|---|---|---|
| LDA、GGA、meta-GGA | 基态结构、能量与响应的不同近似 | 按目标量验证;不同泛函能量不可直接混算 |
| DFT+U | 局域 d/f 电子、占据和关联修正 | 给出 U/J 来源、投影与磁态;不自动解决所有强关联 |
| HSE、PBE0 等杂化泛函 | 带隙、局域化与相关能量改进 | 成本较高;能带、响应与采样路线按手册选择 |
| DFT-D2/D3/D4、TS、MBD、非局域 vdW | 层间、分子晶体、吸附与弱相互作用 | 是不同色散模型,不能随意叠加重复计数 |
| Libxc | 增加可调用的交换关联泛函 | VASP 通用接口从 6.3.0 起;需编译链接,且按版本核对泛函类型 |
| libMBD、外部 DFT-D4 | 增加/改进相应色散实现 | 支持标签不代表库已链接;API 与宿主版本需匹配 |
| BEEF | 对应泛函及其集合思想相关分析 | 按实际库/接口支持使用;集合离散不等于严格总误差条 |
| VTST | 鞍点、路径、优化与 TST 辅助工具 | 补丁可执行程序与脚本分别确认 |
| optcell | 特定晶胞自由度约束 | 用于满足边界,不会生成新的物理理论 |
| VASPsol / VASPsol++ / CP-VASP | 溶剂、电解质与恒电势 | 三者功能和版本范围不同 |
| Wannier 接口 | 局域模型、插值、拓扑/输运后处理 | 独立 Wannier90 版本与 VASP 链接库版本要分清 |
| HDF5 / py4vasp | 结构化输出、分析及部分新算法依赖 | HDF5 是数据基础设施;py4vasp 不是新的电子求解器 |
Libxc 在库中列出某个泛函,不意味着每个 VASP 版本都能使用它及其全部导数。官方注明部分拉普拉斯依赖 meta-GGA、meta-GGA 杂化类型从 6.4.0 才有相应支持;meta-GGA 还涉及 Libxc 构建要求。详见本集群 Libxc 指南和官方 LIBXC1。
色散与编译参考:vdW 方法分类、IVDW 与外部库、编译选项、杂化泛函、DFT+U。
22. 在鸢算科技集群如何选择入口#
以下是本次核对的软件环境范围。文章列出的全网工具并不等于全部已安装;已安装也不等于所有方法组合都适用。 提交前通过 module avail、module show 和版本指南选择完整模块名称。
| 需求 | 当前环境选择 | 平台/配套说明 |
|---|---|---|
| 常规 VASP、GPU 加速 | 6.3.2、6.4.2、6.4.3、6.5.1、6.6.1;另保留 5.4.4 CPU | 6.x 提供 CPU 及 P100/V100 对应分支;具体算法仍需核对 |
| Libxc | 上述五个 6.x 版本的独立 CPU Libxc 分支,Libxc 7.1.2 | 不将 CPU Libxc 能力推广到 GPU 分支;5.4.4 不提供通用原生接口 |
| NEB/过渡态 | 对应版本 VTST 分支与独立脚本 | 按 CPU/GPU 对应模块选择 |
| 经典隐式溶剂 | VASPsol;GPU 为匹配的独立移植分支 | 与 VASPsol++ 分开 |
| VASPsol++ | 匹配 5.4.4 / 6.3.2 的 CPU 组合分支 | 不假定其他宿主或 GPU 同样提供 |
| 恒电势 | CP-VASP 的 6.3.2 / 6.4.2 CPU 分支 | 以相应分支和溶剂模型说明为准 |
| Wannier | VASP 5.4.4 绑定 2.1.0;6.x 绑定 3.1.0 | 独立 Wannier90/postw90 4.0.3 是另外的 CPU 工具,不自动替换绑定库 |
| 能带与常规分析 | VASPKIT、py4vasp、yskit 分段能带 | yskit 的拼接结果面向约定能带数据,不是所有输出文件的通用合并器 |
| 声子与输运 | phonopy、phono3py、ShengBTE、ALAMODE、AMSET、BoltzTraP2 等 | 主机计算与后处理的并行方式分别安排 |
| 键合、拓扑、合金 | Bader、critic2、LOBSTER、WannierTools、VASPBERRY、ATAT 等 | 使用各软件独立模块与数据接口 |
编译思路是将 VASP 宿主版本、编译器/MPI、数学库、可选库和补丁固定成独立环境。CPU 与 GPU 分支的编译工具链不同,Libxc、BEEF、HDF5、Wannier 等能力由具体构建决定;实际编译器与依赖版本集中列在版本与扩展指南及软件环境,不建议自行混接不同模块的库。
计算任务通过集群作业系统提交;登录环境用于准备输入和适量后处理。平台选择先看算法支持,再看规模、内存与速度。vasp_gam、vasp_std、vasp_ncl 也不能仅按“哪个更快”互换。
23. 按研究目标选择最短完整流程#
| 我想研究…… | 建议从哪里开始 | 最终应拿到什么 |
|---|---|---|
| 新材料是否稳定 | 多相优化→形成能/凸包→弹性→全声子→必要有限温度 | 相对稳定性、竞争相及稳定条件 |
| 半导体光电性能 | 收敛结构→能带/DOS→合适带隙方法→介电/GW-BSE→必要输运 | 明确方法的电子隙、光学谱、激子/输运数据 |
| 催化活性 | 位点/覆盖度→吸附→反应网络→热修正/溶剂/电势→势垒 | 统一参照的自由能图和关键动力学步骤 |
| 电池电极 | 多组分结构→稳定相→电压/容量→体积变化→迁移 | 电压区间、计量容量、结构变化与迁移路径 |
| 固态电解质 | 缺陷/迁移路径→NEB→必要 AIMD→扩散和相关电导 | 温度/浓度条件下的迁移和输运机制 |
| 热电材料 | 精细能带→电子输运→非谐声子热导→同条件合并 | S、σ、κe、κl、功率因子与 ZT |
| 磁性/自旋电子学 | 多磁态→SOC→各向异性/交换→必要 Wannier 和统计模型 | 基态磁序、参数、响应及模型温度依赖 |
| 拓扑候选 | SOC 能带→Wannier 校验→不变量/节点→表面态 | 拓扑证据与可解释的边界态 |
| 缺陷色心/发光 | 多电荷态→形成能→ΔSCF/势能面→振动/俘获 | 热力学能级、光学能量与复合机制 |
| 合金相图 | 构型/SQS→一致 DFT 数据→簇展开→Monte Carlo→自由能 | 成分–温度关系及模型适用范围 |
| 液体/无序材料 | 代表性制样→平衡采样→结构统计→扩散/反应 | 带统计误差的结构与动力学性质 |
| 高压/未知结构 | 搜索器→VASP 优化→焓排序→声子→必要自由能 | 候选相、相变范围与稳定性依据 |
这些流程是基于前文方法能力给出的方案建议,不是所有材料共用的固定输入模板。先确定实验或理论需要比较的可观测量,再决定是否必须上更昂贵的方法。
24. 最容易混淆的结论#
- 算得出数值,不等于结果适用于目标体系。 泛函、磁态、势函数和边界会产生系统误差。
- 结构收敛、动力学稳定、热力学稳定是三件事。 一项通过不能替代其他两项。
- KS 带隙、准粒子带隙、光学带隙、缺陷转变能级不是同一个量。
- 声子谱不等于热导率;有效质量不等于迁移率;电导除以弛豫时间不等于绝对电导。
- 吸附能、电子能垒、自由能垒、实验反应速率要分开。 恒电荷与恒电势也不能混用热力学量。
- 磁矩不等于居里温度;极化势垒不等于实验矫顽场;理想强度不等于工程强度。
- 经典 AIMD 不自动包含核量子效应,普通 AIMD 不等于激发态非绝热动力学。
- 有限真空不是免费隔离。 二维光学、带电缺陷、分子、表面必须检查镜像和单位归一。
- GPU 加速不增加物理功能。 某算法/库在 CPU 支持,不能据此认定 GPU 分支支持。
- 外部软件有某功能,不等于 VASP 接口完整提供所需数据。 对 EPW、DMFT、非线性光学、NEGF 等尤其如此。
- 宏观器件效率、循环寿命、吸附等温线、真实催化电流通常需要更多模型。 VASP 提供其中的原子尺度参数,不是完整设备仿真。
- 可复现研究应保留方法与条件。 至少记录版本、PAW、泛函/参数、结构、k 点、截断、磁态、收敛、参照、温压及后处理版本。
25. 继续阅读与资料更新#
本篇给出主要能力地图及实现边界,具体输入请沿各节的官方链接继续阅读。遇到新版本、研究补丁或尚未进入稳定发布的接口,应同时核对论文描述、公开代码、宿主版本和正式手册,以实际可用组合为准。
本站下一步阅读:VASP 版本与扩展、Libxc、分段能带、CPU 参数、GPU 参数、准确性、声子与输运软件、ATAT 专题。
上游持续更新入口:VASP 官方手册、官方版本记录、官方教程。本文不分发 VASP 源码、PAW 势或其他需要单独授权的文件。