Lobster安装和环境变量设置可参考:
VASP+Lobster计算晶体轨道哈密顿布居计算示例
https://mp.weixin.qq.com/s/AvkUNRd6-DbwiaSKpzVyiQ
官网
http://www.cohp.de/
先进行结构优化
然后grep NBANDS OUTCAR
在COHP计算时VASP输入文件需要修改的参数以及提示:
NBANDS = *** #NBANDS>NELECT for COHP, 1.5 or 1.7*NBANDS
##NSW = 0 #单点计算
##LWAVE = T # delete old WAVECAR
##ISYM = -1 #不能处理k点对称性,关闭对称性
##LORBIT = 11 #DOS计算
##NEDOS = 2000 #DOS计算
##POTCAR: PAW #赝势
##VASP: no gam #不支持Gamma版VASP
Lobster输入文件lobsterin模板:
COHPstartEnergy -15 #COHP输出能量范围, eV
COHPendEnergy 10
usebasisset pbeVaspFit2015 #基组, H-Lr(1-103), J. Comput. Chem. 2016, 37, 1030-1035.
basisfunctions Fe 3d 4s #投影到的基组所包含的轨道,需要和使用的赝势一一对应, grep VRHFIN POTCAR
basisfunctions C 2s 2p #占据轨道投影误差须小于5%, 最好小于2%(Spilling误差, 投影后输出的abs. charge spilling)
basisfunctions N 2s 2p
basisfunctions O 2s 2p
basisfunctions H 1s
! cohpbetween atom 72 atom 73 orbitalwise #分析原子72与原子73间的COHP; orbitalwise: 分波COHP
! cohpGenerator from 1.4 to 1.5 ! orbitalwise
cohpGenerator from 1.5 to 2.3 type Fe type N orbitalwise #输出Fe-N原子对距离在1.5到2.3的所有数据之和。
gaussianSmearingWidth 0.05 #高斯展宽,如果在VASP里用的ISMEAR=0, 设置这个关键词和SIGMA一样的值
skipDOS #跳过DOS计算
skipCOOP #跳过COOP计算
! skipCOHP #跳过COHP计算
saveProjectionToFile #把投影计算结果保存到projectionData.lobster,由于投影这一步耗时较长,所以保存这个文件,再换原子对计算的时候节约时间。一般第一次计算的时候运行,后续计算则打loadProjectionFromFile参数
loadProjectionFromFile #读取projectionData.lobster
提交计算:
sbatch cohprun.sh
cohprun.sh包含:
#!/bin/bash
#SBATCH -p v6_384
#SBATCH -N 1
#SBATCH -n 96
lobster-4.1.0
Lobster的输出文件:lobsterout
ICOHPLIST.lobster输出文件,包含积分晶体轨道哈密顿布居
CHARGE.lobster输出文件,包含Mulliken电荷和Loewdin电荷
ICOBILIST.lobster, COBICAR.lobster输出文件,包含共价键指标: 晶体轨道键指数(Crystal Orbital Bond Index, COBI)
ICOBI: 0(离子键)~1(共价键) per bond, 自旋极化体系分spin1,2,具有加和性
MadelungEnergies.lobster, SitePotentials.lobster输出文件,包含离子键指标: 马德隆能(Madelung Energy, in eV per unit cell)
wxDragon下载地址:http://www.wxdragon.del
File - Open 打开COHPCAR.lobster文件