VASP（Vienna Ab initio Simulation Package）是一款基于密度泛函理论（DFT）的量子力学模拟软件，广泛应用于材料科学、化学、物理学等领域。其核心计算方法之一是自洽场（Self-Consistent Field, SCF）算法，解决Kohn-Sham方程，该算法通过迭代计算电子波函数和密度，直到系统达到自洽状态。
自洽计算是VASP计算中最基础也是最重要的步骤，其目的是找到电子波函数和相应的总能量，使得电子在晶胞中的分布达到自洽。自洽计算通常包括电荷密度的初始化、波函数的迭代优化以及能量和力的计算。
VASP计算需要准备超算连接软件EASYCONNECT与SSH，建模软件VESTA，超算连接软件Winscp
jp-minerals.org/vesta/en/download.html
EasyConnect下载-EasyConnect最新版下载V7.6.7.0
Downloading WinSCP-6.5.3-Setup.exe :: WinSCP
ISTART=0 #开始新的任务,随机产生初始波函数
ICHARG=2 #开始新的任务,从原子电荷密度产生体系初始电荷密度
EDIFF=1E-5 #相邻两步电子迭代的能量差收敛标准
EDIFFG=-0.1 #离子弛豫的force的收敛标准
ISMEAR=0 #费米能级附近电子占据数为高斯分布，适合金属、半导体、绝缘体
Automatic generation #注释行
修改 INCAR 文件，让VASP执行自洽计算，修改IBRION=-1、NSW=0、LVHAR=.T.
华算科技朱老师今天给大家介绍 pt001 面的静电势计算，需要在自洽计算中完成，在之前已经做了 001 面的结构优化 pure文件夹这个计算，那么接下来把 pure 复制成 s ，进到自洽计算 s文件夹中，把 CONTCAR 变成POSCAR，然后修改INCAR，先改自洽计算的标签， IBRION=-1， NSW =0，再加上静电势的标签 LVHAR=.T. ，它会输出静电势文件，如果没有这个标签是不输出的，其他就不用改。保存提交任务，待这个计算完成后，会得到一个叫 LOCPOT 的文件，这个文件是三维的静电势数据，一般会把它处理成某一个方向的渐变式变化规律，那么一般会沿着这个带真空的这个方向，就是沿着 z 方向对xy平面做一个平均，这个完成之后得到了 LOCPOT 文件
接下来使用脚本vtot，这里问哪个方向保留，那么我们保留3，也就是 z 方向，它对 XY 平面做一个平均，算完之后得到一个叫 VLINE 的文件。它是一个一维的这个曲线，第一列是它的横坐标，也就是沿着 z 轴的取点个数，第二列就是它的静电势的数据，把这个文件下载下来之后，打开一个新的表格，导入刚才的 VLINE 文件，导入之后把第一行删除，然后做一个图。这就是pt 001 面的静电势分布，这里有两个极小值点，对应两层原子，这里有一个平台，正好对应它的真空区，那么这个横纵坐标可以做一个调整，横坐标 0—160，它一共是取了 160 个点，纵坐标从-20—5，然后再调整一下横纵坐标的名称。横坐标一般用距离，就是 z 轴的距离，然后纵坐标就是静电势。这个曲线稍作修饰就得到了一个静电势的结果。
🎯500+博士团队护航，累计助力5️⃣0️⃣0️⃣0️⃣0️⃣➕篇科研成果，计算数据已发表在Nature & Science正刊及大子刊、JACS、Angew、PNAS、AM系列等国际顶刊。 👏👏👏
声明：如需转载请注明出处（华算科技旗下资讯学习网站-学术资讯），并附有原文链接，谢谢！
