原子振动频率的计算通常基于对系统能量对原子坐标的二阶导数（Hessian 矩阵）的计算。Hessian 矩阵描述了系统在原子位置变化时的能量变化，其特征值对应于振动频率的平方。通过计算 Hessian 矩阵，可以得到系统的振动模式和对应的振动频率。
首先，通过几何优化（如共轭梯度算法）找到系统的最低能量结构，以确保计算的稳定性。在优化后的结构基础上，通过计算 Hessian 矩阵或使用有限位移法（finite difference method）或密度泛函微扰理论（DFPT）来计算振动频率。最后通过分析输出文件（如 OUTCAR）中的振动频率数据，可以得到系统的振动模式和对应的频率值。
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 #费米能级附近电子占据数为高斯分布，适合金属、半导体、绝缘体
NFREE = 2 #中心差分法计算能量对位移的微分
Automatic generation #注释行
修改 INCAR 文件，修改IBRION=5、NSW=1、NFREE=2、POTIM=0.015。
今天华算科技朱老师给大家介绍Pt001面吸附H零点振动能与熵的计算。在之前的视频中，我们已经完成了吸附H的结构优化计算，也就是 h 文件夹。接下来 把 h 文件夹复制成 hfreq 振动计算，进入hfreq，首先把 CONTCAR 变成POSCAR，结构保留，然后修改 INCAR。振动计算的话，这里 IBRION改成5， NSW 改成1，增加 NFREE 等于2， POTIM 等于0.015，这样就把 INCAR 修改好了。然后修改 POSCAR，催化剂原子全部要固定，就是把这里的 T和F 全都改成F，而吸附的H原子还是保留T，也就是催化剂原子不参与振动计算。改完之后提交计算任务，等这个计算完成之后，我们需要用一个脚本去提取它的振动和熵的数值。
这边一共只有一个H原子参与振动计算，所以自由度是 3， X、 Y、 Z 方向各一个自由度，那么它是用有限差分法来算这个能量的，因此有三个自由度。有限差分法会有前向、后向各一个点，就是每个方向有两个计算，总共有 6 个计算。六个方向的能量都已经算好了，算完之后，用到一个脚本叫 zpe来提取它的振动能量 0.14 eV，再用一个脚本叫ts，那么这里要输入温度298.15K，得到它的振动的熵，这个熵是包含温度的这个结果的，这样我们就把零点振动能和熵给大家算完了，这个在催化中是非常常见的。
🎯500+博士团队护航，累计助力5️⃣0️⃣0️⃣0️⃣0️⃣➕篇科研成果，计算数据已发表在Nature & Science正刊及大子刊、JACS、Angew、PNAS、AM系列等国际顶刊。 👏👏👏
声明：如需转载请注明出处（华算科技旗下资讯学习网站-学术资讯），并附有原文链接，谢谢！
