通过计算两个成键原子间的COHP，可以查看成键及反键范围，如图所示：

 

计算过程如下：

一．vasp计算

1.   在自洽计算文件夹中删除原有波函数文件；

2.   自洽计算的INCAR中加入以下参数：

    ISYM=-1

    NSW=0

    NBANDS=*（大致设置为总电子数*80%或更高）总电子数可以通过自洽计算结果中：grep NELE OUTCAR获得

3.   运行vasp计算；

二．lobster计算

1.   修改lobsterin文件

    COHPstartEnergy -14 ！能量窗口，以费米能级为0

    COHPendEnergy 6

    basisfunctions Ga 4s 4p 3d  ！POTCAR中每种原子的价电子轨道

    basisfunctions As 4s 4p

    cohpbetween atom 1 and atom 2 ！计算第一个原子与第二个原子的COHP

最后一行也可以使用：

    cohpGenerator from 1.4 to 1.5 type Ga type As ! 计算所有Ga-As键长小于1.5大于1.4的键。

---

注意：

（1）怎么获得每种原子的价电子轨道？

运行：grep VRH POTCAR，再结合元素周期表，可以得到每种原子的价电子轨道；

示例：

这边轨道（s, p, d）后的数字表示电子数（此处不使用），轨道前的数字才表示相应轨道。结合以下元素周期表

可以写出正确的lobsterin文件：

（2）：以上原子顺序必须和POSCAR一致；

---

2.   运行lobster

命令：lobster-3.0.0

3.   结果：

COHPCAR.lobster文件的第一列为能量值，第四列为第一对原子的COHP值；作图-COHP，需要把第四列乘以负号；