VASP计算COHP应用案例！

【做计算 找华算】燃爆开学季！DFT理论计算全场6折特惠！催化/电池等研究领域全覆盖，限量秒杀，手慢无！

科研党必看！经费预存高至30%增值，更有8500+返利直接送，一次预存，全年无忧！

COHP是衡量原子间相互作用的重要物理量，需要通过VASP和LOBSTER软件一起计算。深圳华算科技有限公司将详细介绍晶体氧化镁COHP计算流程和数据处理方法。

结构优化

首先通过晶体结构数据库下载MgO的cif文件，并生成POSCAR文件，如下图所示。

然后编写结构优化的INCAR和KPOINTS文件。

INCAR文件：

ISTART=0    #随机产生初始波函数

ICHARG=2  #从原子电荷密度产生体系初始电荷密度

PREC=M   #计算精度为Medium

ISPIN=1    #关闭自旋极化

ALGO=N  #DAV算法优化电子波函数

NELM=60   #电子波函数最多计算60步

EDIFF=1E-5  #电子波函数能量收敛标准1E-5 eV 

ENCUT=400  #平面波截断能400 eV

IVDW=11   #考虑范德华力修正

IBRION=2   #共轭梯度法优化晶体结构和原子坐标

NSW=100   #晶体结构和原子坐标优化步数最大100步

ISIF=3   #优化原子坐标和晶体结构

EDIFFG=-0.1  #原子残余力小于0.1 eV/A

ISMEAR=0   #费米能级附近电子占据数为高斯分布

SIGMA=0.1  #高斯分布展宽0.1 eV

KPOINTS文件：

Automatic generation  #注释行

0  #自动产生K点网格

G   #布里渊区K点网格以Gamma点为中心

4 4 4  #K点网格密度

0 0 0  #K点网格中心平移矢量

自洽计算

完成结构优化计算后，保持KPOINTS，POTCAR文件不变，将CONTCAR文件复制成POSCAR文件，并对结构优化的INCAR文件作如下修改：

IBRION=-1  #固定结构和原子位置保持不变

NSW=0  #仅作一步自洽计算

NBANDS=96  #数值为结构优化时的两倍

ISYM=-1  #关闭对称性

COHP计算

VASP计算完成后得到WAVECAR，用LOBSTER软件读取WAVECAR并计算COHP，需要编写输入文件lobsterin

lobsterin：

basisset pbevaspfit2015  #选择VASP拟合基组
cohpbetween atom 1 and atom 7  #给出需要计算COHP的原子编号

数据处理

完成COHP计算后，可得到COHPCAR.lobster，用Origin软件绘制对应原子对COHP曲线即可。

#华算科技 #DFT #VASP #MS #高斯 #电子轨道 #COHP #COOP #成键态 #反键态 #化学键

【做计算 找华算】

🏅 华算科技提供专业的第一性原理、分子动力学、生物模拟、量子化学、机器学习、有限元仿真等代算服务。

🎯500+博士团队护航，累计助力5️⃣0️⃣0️⃣0️⃣0️⃣➕篇科研成果，计算数据已发表在Nature & Science正刊及大子刊、JACS、Angew、PNAS、AM系列等国际顶刊。 👏👏👏

👉 点击阅读原文，立即下单💖