计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › 求助：用ORCA对一系列连续变化的结构计算单点能出现能量 ...

求助：用ORCA对一系列连续变化的结构计算单点能出现能量跳跃点

查看数: 13452 | 评论数: 15 | 收藏 Add to favorites 6

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

enthalpy

发布时间: 2020-8-18 20:38

正文摘要:

本帖最后由 enthalpy 于 2020-8-20 10:38 编辑大家好：我先用gaussian09E 对核苷酸类似物进行二面角进行柔性扫描（10度一步）, 其中计算关键词为# opt=modredundant m062x/def2tzvpp em=GD3 scrf=(smd,solve ...

回复 Reply

enthalpy 发表于 Post on 2020-8-20 10:59:59

zjxitcc 发表于 2020-8-20 10:44
我是随便旋转了一个就是出来了，没有尝试其他的方位朝向。我用的技巧就是GV打开gjf，另存为内坐标，再打 ...

谢谢你的解答。我又学到一些技巧。

zjxitcc 发表于 Post on 2020-8-20 10:44:37

enthalpy 发表于 2020-8-20 10:35
谢谢你的帮助！如果不用SMD模型，ORCA得到结果与Gaussian是下上一致的。将分子坐标旋转有什么讲究吗？还 ...

我是随便旋转了一个就是出来了，没有尝试其他的方位朝向。我用的技巧就是GV打开gjf，另存为内坐标，再打开内坐标gjf，另存为直角坐标，这样分子朝向一般就会改变了。

ORCA的隐式溶剂模型选择性比较低吧，就那么几个，你这碰到的也是特殊情况，。。。结构优化我还是选择用高斯，不会去用ORCA

enthalpy 发表于 Post on 2020-8-20 10:35:23

本帖最后由 enthalpy 于 2020-8-20 10:58 编辑

zjxitcc 发表于 2020-8-19 21:04
我做了些试算：
（1）换用更高精度的格点grid7是解决不了问题的，能量还是相差很多。
（2）ORCA不支持 ...

谢谢你的帮助！如果不用SMD模型，ORCA得到结果与Gaussian是一致的。将分子坐标旋转有什么讲究吗？还是说某种取向得到结果不对，在随机调整试试看？还有一个问题是ORCA4.2.1还有什么可靠的隐式模型可以用？

hebrewsnabla 发表于 Post on 2020-8-19 22:22:01

zjxitcc 发表于 2020-8-19 21:04
我做了些试算：
（1）换用更高精度的格点grid7是解决不了问题的，能量还是相差很多。
（2）ORCA不支持 ...

可以看看CDS能量有无异常，判断一下是极性部分还是非极性部分出了问题。
但是按理说这俩都是旋转不变的啊……

zjxitcc 发表于 Post on 2020-8-19 21:04:21

本帖最后由 zjxitcc 于 2020-8-19 21:05 编辑

enthalpy 发表于 2020-8-19 15:05
谢谢sob老师的解答。我尝试加上kdiis 重算了这三个结构的单点能。发现结构29的能量还是要比另外两个高出 ...

我做了些试算：
（1）换用更高精度的格点grid7是解决不了问题的，能量还是相差很多。
（2）ORCA不支持隐式溶剂模型下检验波函数稳定性，因此这条没法做。但是基本可以排除波函数稳定性问题，因为我从高斯传轨道过去，收敛性很好，但收敛的能量还是很离谱。
（3）将分子坐标旋转过再算（见附件），结果好像就很正常了

n29_m062x_rot.inp (2.42 KB, 下载次数 Times of downloads: 7)

建议你到ORCA论坛上反映这点，督促他们找原因。

喵星大佬 发表于 Post on 2020-8-19 18:13:20

enthalpy 发表于 2020-8-19 18:03
你说到SMD溶剂模型提醒我重新计算在它们气相下的能量，发现之前ORCA计算中出现的跳跃点消失。此时三个结 ...

那看来orca的SMD也有这个问题，我在用Gaussian带SMD溶剂模型优化结构的时候碰到过几次。

enthalpy 发表于 Post on 2020-8-19 18:03:23

喵星大佬发表于 2020-8-19 16:01
这个还算好了，重新算可解决，其实还有个就是SMD溶剂模型在Gaussian里有可能出现很奇葩的情况，好像是因为 ...

你说到SMD溶剂模型提醒我重新计算在它们气相下的能量，发现之前ORCA计算中出现的跳跃点消失。此时三个结构的相对能量变化与gaussian计算的一致。看来时ORCA对某些结构SMD溶剂化处理出问题。

喵星大佬 发表于 Post on 2020-8-19 16:01:14

这个还算好了，重新算可解决，其实还有个就是SMD溶剂模型在Gaussian里有可能出现很奇葩的情况，好像是因为溶剂表面有溶剂进不去的坑，然后SCF的时候能量会巨幅变化然后报错。那个解决起来很麻烦麻烦程度相当高。

enthalpy 发表于 Post on 2020-8-19 15:05:29

本帖最后由 enthalpy 于 2020-8-19 15:08 编辑

sobereva 发表于 2020-8-19 11:28
KDIIS是SCF收敛方法，和几何优化的GDIIS差异极大

谢谢sob老师的解答。我尝试加上kdiis 重算了这三个结构的单点能。发现结构29的能量还是要比另外两个高出~60kJ/mol.
@zjxitcc， ORCA4.2.1计算得到的电子能量分别为 -996.647693959302， -996.624219156508， -996.647318397864。
gaussian09E 得到为 -996.652346174，-996.652659893， -996.652691888

ps 附件是ORCA4.2.1 计算时的输入文件。