计算化学公社

标题: 求助：vasp算IR光谱 [打印本页]

作者
Author: 不明嚼栗子 时间: 2024-7-16 22:57
标题: 求助：vasp算IR光谱
老师同学们大家好，我参考论坛里的一个帖子用vasp+phonopy算晶体的红外，但是我算完DFPT输出的vasprun.xml用phonopy --fc vasprun.xml --hdf5时报错'vasprun.xml' dones not contain necessary information. 请问这个是怎么回事啊？？？
PREC = Accurate
ICHARG =    2
ENCUT = 450
EDIFF = 1E-05
LREAL = Auto
LWAVE = .FALSE.
LCHARG = .FALSE.
IALGO = 38
IBRION = 8
LPEAD = .TRUE.
LEPSILON = .TRUE.
LMAXTAU = 0
EDIFFG =  -0.02
NSW =    0
ISMEAR =    0
SIGMA=0.05
谢谢大家！拜托了！

作者
Author: archer 时间: 2024-7-17 18:05
phononpy不能算红外吧，https://phonopy.github.io/phonopy/external-tools.html。你确定任务成功结束了吗？

作者
Author: logzzz 时间: 2024-7-17 22:34
本帖最后由 logzzz 于 2024-7-17 22:36 编辑

NSW = 1， IBRION =8算一遍，IBRION =-1，NSW =0，LEPSILON = .TRUE. 算一遍。

作者
Author: 不明嚼栗子 时间: 2024-7-19 10:05

archer 发表于 2024-7-17 18:05
phononpy不能算红外吧，https://phonopy.github.io/phonopy/external-tools.html。你确定任务成功结束了吗 ...

成功了，但是好像vasprun.xml中没有hessian，只有born电荷。。

作者
Author: 不明嚼栗子 时间: 2024-7-19 10:06

logzzz 发表于 2024-7-17 22:34
NSW = 1， IBRION =8算一遍，IBRION =-1，NSW =0，LEPSILON = .TRUE. 算一遍。

请问这样做是什么原理呀，我算了两遍但是感觉好像没有变化，还是同样的报错。。

作者
Author: 卡开发发 时间: 2024-7-19 10:28

不明嚼栗子发表于 2024-7-19 10:06
请问这样做是什么原理呀，我算了两遍但是感觉好像没有变化，还是同样的报错。。

你要把两步分到两个任务去，算声子做完用phonopy提取Hessian信息，然后再进行一个任务单独算Born电荷，然后按照2#提供的工具去处理，当然更具体的程序和例子你可以看这里。

如果还不行你可能要检查一下OUTCAR当中声子计算的部分是不是有问题。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)