VASP 专题

VASP 体系与性质计算全景指南

从晶体、表面、缺陷到电化学与激发态,按研究目标查找 VASP 原生功能、扩展工具、计算流程和适用边界。

更新于 2026-10-24约 28 分钟阅读按任务与环境选择适用配置
阅读本篇目录

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 差不等于全部实验量

推荐先明确“需要一张图还是需要定量带边”。带隙、光学峰和激子能是不同目标,后两者见后文。显存不足时可使用分段能带工作流,但分段路径不能替代输运需要的完整均匀网格。

参考:能带、态密度、电荷密度、LPARD、功函数计算。

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. 最容易混淆的结论#

  1. 算得出数值,不等于结果适用于目标体系。 泛函、磁态、势函数和边界会产生系统误差。
  2. 结构收敛、动力学稳定、热力学稳定是三件事。 一项通过不能替代其他两项。
  3. KS 带隙、准粒子带隙、光学带隙、缺陷转变能级不是同一个量。
  4. 声子谱不等于热导率;有效质量不等于迁移率;电导除以弛豫时间不等于绝对电导。
  5. 吸附能、电子能垒、自由能垒、实验反应速率要分开。 恒电荷与恒电势也不能混用热力学量。
  6. 磁矩不等于居里温度;极化势垒不等于实验矫顽场;理想强度不等于工程强度。
  7. 经典 AIMD 不自动包含核量子效应,普通 AIMD 不等于激发态非绝热动力学。
  8. 有限真空不是免费隔离。 二维光学、带电缺陷、分子、表面必须检查镜像和单位归一。
  9. GPU 加速不增加物理功能。 某算法/库在 CPU 支持,不能据此认定 GPU 分支支持。
  10. 外部软件有某功能,不等于 VASP 接口完整提供所需数据。 对 EPW、DMFT、非线性光学、NEGF 等尤其如此。
  11. 宏观器件效率、循环寿命、吸附等温线、真实催化电流通常需要更多模型。 VASP 提供其中的原子尺度参数,不是完整设备仿真。
  12. 可复现研究应保留方法与条件。 至少记录版本、PAW、泛函/参数、结构、k 点、截断、磁态、收敛、参照、温压及后处理版本。

25. 继续阅读与资料更新#

本篇给出主要能力地图及实现边界,具体输入请沿各节的官方链接继续阅读。遇到新版本、研究补丁或尚未进入稳定发布的接口,应同时核对论文描述、公开代码、宿主版本和正式手册,以实际可用组合为准。

本站下一步阅读:VASP 版本与扩展、Libxc、分段能带、CPU 参数、GPU 参数、准确性、声子与输运软件、ATAT 专题。

上游持续更新入口:VASP 官方手册、官方版本记录、官方教程。本文不分发 VASP 源码、PAW 势或其他需要单独授权的文件。

需要根据你的体系选择资源?联系讨论计算需求。实际可用软件与权限,以当前账户和运行环境为准。

查看配图

放大后可滚动查看细节,按 Esc 关闭。

100%查看原图 ↗