Lobster计算COOP和COHP(以IrO2为例) | 您所在的位置:网站首页 › him如何使用 › Lobster计算COOP和COHP(以IrO2为例) |
写在最前面COOP和COHP用来干嘛? 一句话:判断成键和反键轨道的状态,具体多看看文献家人们 Lobster的官网:http://www.cohp.de/
龙虾程序是由Richard Dronskowski开发的。
由于表示波函数的希腊字母ψ和龙虾很像,因此程序的图标设计成了龙虾的样子
极简安装,直接tar -zxvf 即可安装,cp lobster-4.1.0到./bin文件即可 lobster安装后自带多个计算例子和用户手册 以VASP为例,cp实例中的lobsterin文件到你要计算的目录下,lobsterin里面的参数基本如下 手册中都会具体写哪些是默认值,一般不需要调太多, 我一般会加上这个参数 cohpGenerator from 1.4 to 1.5 type A type B 其中A和B为要计算的原子,以金红石IrO2为例,其存在两种键长,可以判断不同键长的COHP 计算时候,先写好lobsterin文件,将已经优化好结构的CONTCAR复制为POSCAR,删除WAVECAR(不删除直接算会报错~具体因为什么我不知道),然后进行静态计算 NSW = 0 ISIF = 0 ISYM = 0 NBANDS =多少,如果不知道多少的话,可以先不写入这个参数,然后算完了以后运行lobster会自动告诉你,计算COHP最低需要多少bands. 再填一遍就好了。 生成的文件我主要看COHPCAR.lobster 打开以后里面的数据很多列,都代表什么意思呢?文档中是这样写的: 实际的数据格式如下,以IrO2为例(本文所用的晶胞为根号2的supercell,为了和goodenough模型一致,详见(原创)关于如何准确计算IrO2态密度中的分轨道分布(复现一篇PRL)) 第一列为能量 第二列为全部原子对的pCOHP 第三列为第二列的积分 第四列为No.1原子对的pCOHP 第五列为第四列的积分 ….. 以IrO2为例,计算得到的Ir-O键的COHP画图如下: 对比一下一篇nc中计算的IrO2的COHP(Iridium metallene oxide for acidic oxygen evolution catalysis): 完美~
|
CopyRight 2018-2019 实验室设备网 版权所有 |