分析化学成键特性（COHP）的两个软件

Program1：Lobster http://www.cohp.de/

Lobster可以视为一后处理软件，它是基于（VASP计算得到的）波函数文件进行处理的。采用VASP计算时，该软件要求采用全K点进行计算（即要求设置ISYM=-1，所得K点的权重均为1）。另外，该软件对赝势有要求，比如不得使用_sv的赝势。要分析的原子/离子对通过输入文件lobsterin进行设置，有指定原子和指定距离设置两种形式。具体要求和计算流程详见手册。

Lobsterin例子：

COHPstartEnergy  -8   ！能量起始值

COHPendEnergy   12   ！能量终止值。这两个参数定义能量范围，通常通过DOS判断

includeOrbitals s p   ！价电子的轨道

cohpGenerator from 4.3 to 4.4   ！指定待分析的原子对的距离

!cohpbetween atom 1 and atom 2    ！或指定待分析的原子对的名称

Program2：STUTTGART TB-LMTO https://www2.fkf.mpg.de/andersen/

STUTTGART TB-LMTO-ASA是一独立的计算软件。输入文件（CTRL）设置基本结构信息（类似于VASP的POSCAR）即可，其它参数可通过命令进行补全。建立输入文件后，进行自洽计算，然后对要分析的原子/离子对之间的成键特性（COOP、COHP、ICOHP等）进行分析。该软件通过指定原子球的名字来定义要分析的原子/离子对。该软件的优点之一是可以对电子化合物中间隙局域电子（ISQ）和原子/ISQ之间的成键特性进行分析。具体要求和计算流程详见手册。

CTRL例子：

STRUC ALAT=11.2683 ！lattice constant （单位是bohr）

PLAT=1 0 0

0 1 0

0 0 1    ！lattice translation vectors

CLASS ATOM=NI Z=28 ！原子名字和原子序数

ATOM=Se Z=34

SITE ATOM=NI POS= .0 .0 .0    ！原子名字和原子坐标

ATOM=NI POS= .5 .0 .5

ATOM=NI POS= .0 .5 .5

ATOM=NI POS= .5 .5 .0

ATOM=Se POS= .383 .383 .383

ATOM=Se POS= .617 .617 .617

ATOM=Se POS= .117 .617 .883

ATOM=Se POS= .883 .383 .117

ATOM=Se POS= .617 .883 .117

ATOM=Se POS= .383 .117 .883

ATOM=Se POS= .883 .117 .617

ATOM=Se POS= .117 .883 .383

P.S. 两款软件都需要向开发者发送邮件进行申请。