计算化学公社

 找回密码 Forget password
 注册 Register

ORCA中SMD的CPCM疑似存在bug

查看数: 9864 | 评论数: 8 | 收藏 Add to favorites 2
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2020-4-28 12:04

正文摘要:

本帖最后由 linq5203 于 2020-4-30 14:28 编辑 我对一批同系物使用M062X/6-31G*级别靠SMD计算溶解自由能时,发现ORCA对其中一个物质在正辛醇中溶解自由能的计算有问题。(ORCA好像不支持M052X)对于同系物,溶解 ...

回复 Reply

LiHuaYu 发表于 Post on 2023-9-13 09:17:14
Orca 存在数百万个错误,从溶剂化模型实现(SMD、CPCM 等)到激发态(TDDFT、CASSCF 等)以及与 DLPNO-CCSD 方法相关的其他内容。即使是简单的 HF 实现在非平衡距离上也存在一些错误,因此 CCSD(T) 。我们对 ORCA 进行了很多基准测试,并且始终避免使用 ORCA 进行实际计算。到目前为止只使用了EDS模块,但似乎也不好。

所以:使用 GAUSSIAN 进行研究,使用 MOLCAS 和 MOLPRO 进行多参考计算,远离 ORCA

Google Translate may be not good, sorry.
sobereva 发表于 Post on 2023-9-12 22:50:24
JCenter 发表于 2023-9-12 22:40
老师,想请教您几个疑问:1.ORCA中SMD溶剂化模型及算完单点后,输出文件中的SMD CDS free energy correct ...

随便找个简单体系,用HF/def2-SVP对比一下两个程序SMD的输出便知,如果能找到和Gaussian输出的非极性部分的数值差不多的量自然就是那个了
JCenter 发表于 Post on 2023-9-12 22:40:44
sobereva 发表于 2020-4-30 07:14
如果你的计算方式合理,那么按理说应当是ORCA的SMD的非极性部分比较诡异。
按理说ORCA和Gaussian的SMD的非 ...

老师,想请教您几个疑问:1.ORCA中SMD溶剂化模型及算完单点后,输出文件中的SMD CDS free energy correction energy是不是对应着溶剂对溶质的非静电或者叫非极性部分的能量贡献。2,极性部分的贡献是不是通过用SMD溶剂模型下的单点能减去真空下的单点能得到一个总的溶剂化能量,再用这个能量再减去上述的非极性部分,即可得到极性部分(静电部分),我看了ORCA中的手册,貌似输出文件没有直接给出极性部分的能量。3.溶剂化能,极性部分和非极性部分的值均为正代表什么物理意义呢
linq5203 发表于 Post on 2020-4-30 14:26:25
sobereva 发表于 2020-4-30 07:14
如果你的计算方式合理,那么按理说应当是ORCA的SMD的非极性部分比较诡异。
按理说ORCA和Gaussian的SMD的非 ...

加入了非极性部分的对比,确实是就这一个结构的SMD中的CPCM的极性部分出现问题
sobereva 发表于 Post on 2020-4-30 07:14:09
如果你的计算方式合理,那么按理说应当是ORCA的SMD的非极性部分比较诡异。
按理说ORCA和Gaussian的SMD的非极性部分应该是一样,应当对比一下看看是否有显著差别
linq5203 发表于 Post on 2020-4-29 09:49:34
sobereva 发表于 2020-4-29 02:09
表格里每一项是什么含义都应当详细说明,否则理解起来很累

CPCM本来就和Gaussian默认的IEFPCM形式差异很 ...

感谢sob老师的提醒,已经改正
sobereva 发表于 Post on 2020-4-29 02:09:35
表格里每一项是什么含义都应当详细说明,否则理解起来很累

CPCM本来就和Gaussian默认的IEFPCM形式差异很大,所以不能指望能够很好重现Gaussian的溶解自由能的结果。并且还需要把情况具体分析,看到底是极性部分带来的差异还是非极性部分就有所不同。
plus 发表于 Post on 2020-4-28 13:48:10
建议发到ORCA论坛或者给Frank Neese写信

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-26 17:48 , Processed in 2.296056 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list