第六章：VASP查看赝势方法指南! | 2026新版VASP基础教程

引言

在上一章《第五章：VASP多元素赝势合并指南！| 2026新版VASP基础教程》中，我们详细介绍了VASP的POSCAR文件。VASP输入文件中的 POTCAR提供了计算所需要原子赝势。本章华算科技朱老师将正式介绍POTCAR，从价电子、芯电子、阶段能等角度来介绍VASP计算所需要原子赝势。

赝势文件的结构与重要性

在深入使用 grep 命令之前，我们首先需要了解 POTCAR 文件的内部结构。POTCAR 文件通常由多个原子块（Atomic Blocks）组成，每个块对应一个元素的赝势信息。

POTCAR 文件的基本结构

每个原子块的开头通常包含以下几行关键信息：

标题行 (Title Line)

通常以 TITEL 或 TITLE 开头，紧跟元素符号和赝势描述。

电子数 (ZVAL)

描述了该赝势中价电子的数目。

能量截断 (ENMAX)

该赝势使用时平面波基组的最大能量截断。

例如，一个典型的 POTCAR 开头可能如下所示：

TITEL = PAW_PBE H  07Sep2000 

ZVAL = 1.00 

ENMAX = 400 eV 

为什么要查看Title？

Title 行包含了赝势的核心信息：

赝势类型

是 PAW 还是 USPP（超软赝势）。

交换关联泛函

是 LDA 还是 GGA（如 PBE）。

是否包含半芯态

如 Fe_pv 表示 Fe 元素包含 p 轨道的半芯态，这对于高压计算非常重要。

基础使用：快速查看标题

查看所有元素的Title

如果你已经将多个元素的赝势文件合并成了一个 POTCAR（如用于 SiO₂ 的 POTCAR），可以使用以下命令快速查看所有元素的 Title 行：

grep TIT POTCAR 

或者：

grep TITLE POTCAR 

这条命令会输出类似以下内容：

TITEL = PAW_PBE H 

TITEL = PAW_PBE Si 

TITEL = PAW_PBE O 

只查看特定元素的Title

如果你只关心某个特定元素（如 Fe）的信息，可以结合 grep 的 -A（after）选项来提取：

grep -A 2 “Fe” POTCAR 

这会输出 Fe 元素块中 Title 以及紧接着的两行（通常包括 ZVAL 和 ENMAX）。

高级技巧：结合awk或sed提取更精准的信息

虽然 grep 可以快速定位标题，但有时我们需要提取更精确的内容（例如，仅获取 ENMAX 值）。此时可以将 grep 与 awk 或 sed 结合使用。

提取所有元素的 ENMAX

grep ENMAX POTCAR | awk ‘{print $3}’ 

这条命令会输出 POTCAR 中所有元素的 ENMAX 值列表。

提取特定元素的 ZVAL

grep -A 1 “Fe” POTCAR | grep ZVAL | awk ‘{print $3}’ 

这条命令会输出 Fe 元素的价电子数（ZVAL）。

本章要点总结

查看POTCAR信息知识点梳理

赝势文件 (POTCAR) 的结构与重要性：介绍了POTCAR中参数意义

基础使用快速查看标题 (Title)：介绍了使用grep命令查看POTCAR信息

高级技巧结合 awk 或 sed 提取更精准的信息：介绍了使用sed和awk命令查看POTCAR信息

下一步学习建议

下一章将正式引入本次教程的核心—KPOINTS文件中Gamma 与 Monkhorst-Pack 区别。我们将从KPOINTS 文件基础结构、Gamma 采样与 Monkhorst-Pack 采样的原理与区别、详细对比与选取建议方面详细介绍VASP KPOINTS文件，以及他们在VASP计算中的应用，敬请期待！