VASP（Vienna Ab initio Simulation Package）是一款基于密度泛函理论（DFT）的量子力学模拟软件，广泛应用于材料科学、化学、物理学等领域。其核心计算方法之一是自洽场（Self-Consistent Field, SCF）算法，解决Kohn-Sham方程，该算法通过迭代计算电子波函数和密度，直到系统达到自洽状态。
自洽计算是VASP计算中最基础也是最重要的步骤，其目的是找到电子波函数和相应的总能量，使得电子在晶胞中的分布达到自洽。自洽计算通常包括电荷密度的初始化、波函数的迭代优化以及能量和力的计算。
VASP计算需要准备超算连接软件EASYCONNECT与SSH，建模软件VESTA，以及VASP的输入文件INCAR，KPOINTS，POTCAR
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 #注释行
最后用hs-1的能量减去hs-2和hs-3的能量，即可得到结合能
今天华算科技朱老师给大家介绍Pt001面吸附H结合能计算。结合能是反映催化剂和被吸附对象之间作用强度的一个物理量，来看一下具体该如何计算。在上一次视频中已经做好了 H 的吸附，那么结合能，先要把吸附构型的能量算一遍，再把H和Pt拆开，各自能量算一遍，把整体能量减去两个拆开后结构的能量就得到了它的结合能。首先把 h 变成 hs-1 就是吸附构型的自洽计算，算一下它的总能，进到这个文件夹之后把 CONTCAR 变成POSCAR，保留下来。修改 INCAR 执行自洽计算，IBRION再改成-1， NSW 改成0，其他就不用改了，结构不用动，调整一下自洽计算的参数就行。提交任务，接下来做第二个自洽计算，就是把整体的结构中移除H，就把 hs- 1 变成 hs-2，这时候INCAR就一样了，不用动了，修改 POSCAR。把这个 H删除，以及对应的原子坐标也删除，保存，提交任务。第三个，把整体结构中的Pt移除，单独算H的能量，把 hs-1 变成 hs-3，进到这个里面修改POSCAR，把Pt删除，留下一个H原子，这时候赝势做一个调整，用H的赝势去覆盖它，这样就可以提交。
来看一下三个自洽计算做的怎么样？第一个是整体的算好了，第二个只有Pt的也算好了，第三个只有H的也算好。那么接下来把能量来计算一下，首先把三个能量输出到一起，可以抄下来，也可以用这个 tail 命令 来一起输出。那就把 s1 的减s2，减 s3 就可以。
这边来做一个减法，减下来之后是- 2.75 电子伏特，这就是它的结合能。这里简单说一下，结合能的数值一定是负值，反应之间有相互作用。如果说没有相互作用，无限远的情况下结合能是0，永远不会算出是正数来，那么这个数值越负，说明相互作用越强，这个Pt-H之间是一种共价键，作用更强的话这个数会更低。
🎯500+博士团队护航，累计助力5️⃣0️⃣0️⃣0️⃣0️⃣➕篇科研成果，计算数据已发表在Nature & Science正刊及大子刊、JACS、Angew、PNAS、AM系列等国际顶刊。 👏👏👏
声明：如需转载请注明出处（华算科技旗下资讯学习网站-学术资讯），并附有原文链接，谢谢！
