VASPKIT专题

VASPKIT：压电与介电张量

压电应力张量e、压电应变张量d、介电常数和弹性矩阵需要使用同一坐标约定。此处用官方示例的VASP结果在集群重新提取和组合。

更新于 2026-10-26约 4 分钟阅读按任务与环境选择适用配置

阅读本篇目录

压电应力张量e、压电应变张量d、介电常数和弹性矩阵需要使用同一坐标约定。此处用官方示例的VASP结果在集群重新提取和组合。

先准备示例#

unzip vaspkit-examples-*.zip
cd vaspkit-examples
python3 prepare.py
module load vaspkit/1.5.1

prepare.py从集群安装的官方examples复制需要的原始示例；附件本身不转发这些分发文件。每个功能的目录独立，下面的cd cases/...均从附件根目录开始。完成一个功能后可用cd ../..回到根目录。

主要例子的文件#

安装包示例 piezoelectric，在集群重新执行本功能；不是本轮重新进行全部上游VASP计算。

完整的配套文件通过附件和prepare.py准备。若本功能需要WAVECAR/WAVEDER，先按新计算与波函数准备运行自己的VASP，不使用其他体系的波函数代替。

各功能操作#

442 · 提取介电与压电张量#

按下面顺序输入；原子编号从1开始，元素名与当前POSCAR一致。更换体系时先核对所需的文件和编号。

实例来源：安装包示例 piezoelectric，在集群重新执行本功能；不是本轮重新进行全部上游VASP计算。

在附件根目录执行：

cd cases/442
vaspkit < answers.in

answers.in实际内容如下，每行对应一次输入；也可运行vaspkit后逐行手动键入。

目录中的主要配套文件：KPOINTS、OUTCAR、POSCAR。

本次得到：DIELECTRIC_TENSOR、DIELECTRIC_TENSOR_ELECTRONIC、DIELECTRIC_TENSOR_IONIC、PIEZOELECTRIC_STRESS_TENSOR、PIEZOELECTRIC_STRESS_TENSOR_ELECTRONIC、PIEZOELECTRIC_STRESS_TENSOR_IONIC。原始输入与结果保持在同一例子目录；表格文件先查看首行列名，再选择横纵轴。

可用head 文件名查看表头，结构结果可在VESTA中打开，曲面数据按对应网格绘制。

443 · 压电应变与机电耦合#

按下面顺序输入；原子编号从1开始，元素名与当前POSCAR一致。更换体系时先核对所需的文件和编号。

实例来源：安装包示例 piezoelectric，在集群重新执行本功能；不是本轮重新进行全部上游VASP计算。

在附件根目录执行：

cd cases/443
vaspkit < answers.in

answers.in实际内容如下，每行对应一次输入；也可运行vaspkit后逐行手动键入。

目录中的主要配套文件：KPOINTS、OUTCAR、POSCAR。

本次得到：ELECTROMECHANICAL_COUPLING_COEFFICIENTS、PIEZOELECTRIC_STRAIN_TENSOR。原始输入与结果保持在同一例子目录；表格文件先查看首行列名，再选择横纵轴。

本例输出摘要（数值仅对应上述输入）：

0.000000   0.000000   0.000000  -0.320540   0.000000   0.000000
-0.640800  -0.640800   1.602830   0.000000   0.000000   0.000000
0.000000   0.000000   0.000000  -0.339200   0.000000   0.000000
0.000000   0.000000  -0.320540   0.000000   0.000000   0.000000
-0.640800   1.602830   0.000000   0.000000   0.000000   0.000000
0.000000   0.000000   0.000000   0.000000  -2.691916   0.000000
0.000000   0.000000   0.000000  -2.543829   0.000000   0.000000
-2.130027  -2.130027   5.115841   0.000000   0.000000   0.000000
0.000000   0.000000   0.000000   0.000000  -0.113359   0.000000
0.000000   0.000000   0.000000  -0.107123   0.000000   0.000000
-0.139249  -0.139249   0.336315   0.000000   0.000000   0.000000
[1] Wu X.F., Vanderbilt D., et al. Phys. Rev. B 72, 035105 (2005)

换成自己的体系#

先复制一份计算目录，确认POSCAR、INCAR、KPOINTS和各结果来自同一计算，再运行所需菜单。涉及静态、优化或声子等前置VASP任务，先完成计算再后处理；VASPKIT不会替代这些计算。保留结构、输入参数、工具版本和提取后的数据，后续改图时不必反复重算。

方法参考：VASPKIT官方教程；菜单与交互输入以本专题实际使用的安装版本为准。使用结果开展研究时引用VASPKIT论文。

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