在上一章《第四章：VASP赝势选择原则！| 2026新版VASP基础教程》中，华算科技朱老师介绍了VASP的POSCAR文件。VASP输入文件中的 POTCAR提供了计算所需要原子赝势。本章将正式介绍POTCAR，从价电子、芯电子、阶段能等角度来介绍VASP计算所需要原子赝势。
VASP的计算核心在于POTCAR文件，它是由每个元素的赝势文件（POTCAR或POTCAR.Z）按原子顺序拼接而成的。
在VASP中，POTCAR文件负责描述系统中各元素的原子核与价电子相互作用。VASP推荐使用投影缀加波（PAW）方法的赝势，因为它兼具全电子精度和高效性，尤其适合处理过渡金属和高压体系。PAW方法实际上蕴含了传统赝势的理论，是更一般的理论框架。
合并POTCAR的标准流程如下：
步骤一：选择合适的赝势目录
VASP自带多个赝势目录，最常用的是potpaw_PBE（GGA-PBE泛函的PAW赝势）。如果你的体系中包含4d、5d或f电子，需要额外注意使用_pv（半芯态）或_sv（软阱）赝势。
步骤二：解压单元素赝势
进入赝势目录后，你会看到类似POTCAR.Z的压缩文件。使用zcat POTCAR.Z > Element（将Element替换为具体的元素符号）解压。例如：
zcat POTCAR.Z > Ga
zcat POTCAR.Z > As
这一步是生成可供合并的中间文件。
步骤三：合并并检查顺序
最关键的一步是合并且顺序必须与POSCAR文件中的原子顺序严格一致。例如，对于GaAs体系：
cat Ga As > POTCAR
如果顺序反了（如cat As Ga > POTCAR），计算结果会完全错误。
步骤四：确认能量截断（ENCUT）
合并完成后，建议查看POTCAR文件头部的ENMAX值。VASP计算时的平面波截断能必须不低于POTCAR中所有元素的最大ENMAX。如果你的计算需要较高精度，建议将ENCUT设置为ENMAX * 1.3左右。
4d/5d元素
强烈建议使用X_sv（软阱）或X_pv（半芯）赝势。例如，对于钼（Mo）和钨（W），Mo_sv和W_sv通常比默认赝势更稳定。
稀土与f电子
如果处理稀土元素（如Ce、Pr），通常只需要保留4f电子价态（使用_pv或_sv），全电子计算非常耗时。
对于大型体系，建议使用_POTCAR作为后缀（如Ga_POTCAR），这样可以避免覆盖默认的POTCAR文件。在脚本自动化中，通常使用如下伪代码：
for element in $(cat POSCAR | grep -v ‘#’); do
cp /path/to/potpaw_PBE/$element/POTCAR.Z ./$element.POTCAR.Z
zcat $element.POTCAR.Z >> POTCAR
done
幽灵态（Ghost State）
如果计算结果出现非物理的高能态，通常是赝势选择不当。解决办法是更换硬度更高的赝势（如使用_sv代替默认）。
压力计算错误
如果计算体系在高压下（如>100 GPa），需要检查赝势的适用压力范围。部分软赝势（_s）在高压下失效，需要更换硬赝势。
ASE（Atomic Simulation Environment）
对于多元素体系，手动合并容易出错。建议使用ASE的VaspIO插件，它可以自动根据POSCAR中的元素顺序读取并合并POTCAR，极大提高效率。
POTCAR生成脚本
VASP官方通常提供了一个名为pseudopotential_select.py的脚本，输入元素列表即可自动生成POTCAR。
核心原理与标准流程：介绍了zcat和cat命令解压和合并赝势的步骤
高级技巧：介绍了for循环与ASE命令合并赝势方法
下一章将正式引入本次教程的核心—vasp查看赝势方法。我们将从赝势文件 (POTCAR) 的结构与重要性、基础使用快速查看标题 (Title)、高级技巧结合 awk 或 sed 提取更精准的信息方面详细介绍VASP赝势选择原则，以及他们在VASP计算中的应用，敬请期待！
