计算化学公社

标题: dimer计算过渡态时输出文件的诡异现象 [打印本页]

作者
Author:
夜夜衣    时间: 2018-2-17 21:21
标题: dimer计算过渡态时输出文件的诡异现象
本帖最后由 夜夜衣 于 2018-2-17 21:23 编辑

各位老师前辈,我在用CINEB计算过渡态后得到了一个能量最高的像点,想着拿来继续跑下dimer,但是一跑就把我给郁闷的不轻,因为过渡态的能量总会直接升高18eV,检查了很久也没出现问题。然后今天在研究输出文件OUTCAR和OSZICAR事后发现了个诡异的现象。先看OSZICAR文件(图一)可以看到一直到17步能量差不多都稳定在-668.8eV左右,但是我用命令grep 'energy without entropy' OUTCAR 却得到了图四。。。,然后呢百思不得其解,又跑去检查OUTCAR,结果发现了OUTCAR中有一次电子步的能量直接从-668跳到了651eV,(图三到图二,抱歉图的顺序有误),(可以看行数)然后以为自己对这个不熟悉误解了,然后去检查其他的文件,发现OSZICAR和OUTCAR中的能量都匹配的很好,就这个有问题,那么现在问题来了,我是不是应该用OSZICAR文件中的数呢还是OUTCAR中的数呢

作者
Author:
archer    时间: 2018-2-20 10:10
用你算出的结构,跑下频率验证下是否是过渡态。
作者
Author:
夜夜衣    时间: 2018-2-20 15:56
archer 发表于 2018-2-20 10:10
用你算出的结构,跑下频率验证下是否是过渡态。

我现在算的结构只能在gamma点计算频率,优化用的551,在我关注的几个原子上还是有虚频的,没法准确判断是否是过渡态
作者
Author:
shaddock    时间: 2020-10-9 17:26
您好,请问MODECAR这个文件是怎么生成的呢,非常感谢!!!
作者
Author:
Stars    时间: 2023-12-6 00:40
虽然不知道是什么原因,但我发现在不同的k点下跑出来的同一结构的磁矩都是不一样的,可能导致很大的能量变化




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3