计算化学公社

标题: Gaussian计算中进行BSSE校正与scrf=pcm一起使用出现问题 [打印本页]

作者
Author: 978142355 时间: 2015-9-22 23:15
标题: Gaussian计算中进行BSSE校正与scrf=pcm一起使用出现问题
本人研究的是鸟嘌呤与胞嘧啶之间的氢键作用，打算使用PCM对G:C做BSSE校正。命令行是：#p mp2/aug-cc-pVDZ counterpoise=2 scrf=pcm

计算结果出现错误：

No SCRF with Counterpoise.
Error termination via Lnk1e in d:\software\Gaussian09\l1.exe at Tue Sep 22 23:04:41 2015.

请问有什么方法解决此问题？看错误的提示是counterpoise不能与scrf一起用，但是我在气相校正时用的是BSSE校正，如果在水相下如果不用BSSE校正，出来的结果与气相对比合理吗？计算所用的gjf和计算的out文件在附件中。

作者
Author: sobereva 时间: 2015-9-22 23:56
程序不支持这样联用，但你可以自行用Bq原子来做counterpoise过程，就没这个限制了。
见
谈谈BSSE校正与Gaussian对它的处理
http://sobereva.com/46

作者
Author: 978142355 时间: 2015-9-23 12:06

sobereva 发表于 2015-9-22 23:56
程序不支持这样联用，但你可以自行用Bq原子来做counterpoise过程，就没这个限制了。
见
谈谈BSSE校正与Ga ...

谢大神指点，我去看看。

作者
Author: 978142355 时间: 2015-9-24 12:41

sobereva 发表于 2015-9-22 23:56
程序不支持这样联用，但你可以自行用Bq原子来做counterpoise过程，就没这个限制了。
见
谈谈BSSE校正与Ga ...

scrf=pcm的BSSE校正根据老师的帖子已经完美解决。
看完帖子有一处不明，即“在计算弱相互作用时counterpoise校正能只用一半效果最佳”这个还能进行一半的操作吗？

作者
Author: sobereva 时间: 2015-9-24 19:27

978142355 发表于 2015-9-24 12:41
scrf=pcm的BSSE校正根据老师的帖子已经完美解决。
看完帖子有一处不明，即“在计算弱相互作用时counterp ...

你就把BSSE校正能除以2再加到体系原先总能量上就行了

作者
Author: 978142355 时间: 2015-9-24 20:02

sobereva 发表于 2015-9-24 19:27
你就把BSSE校正能除以2再加到体系原先总能量上就行了

哦，这样啊，谢谢Sob老师。

作者
Author: wbn 时间: 2016-12-24 06:43
最好把这个帖子的内容也附在后面... 我用ghost原子做了半天结果bsse全是负的...
http://bbs.keinsci.com/thread-3049-1-1.html

sobereva 发表于 2016-4-14 21:38
溶剂模型下计算相互作用能没有严格的办法考虑BSSE问题，相对可行的做法是使用气相下得到的BSSE校正能，加到溶剂下相互作用能上去。
溶剂模型下用Bq不合适，因为Bq原子在溶剂模型计算时也会产生相应的孔洞，而不仅仅是表现相邻原子的基组了，这会导致BSSE校正能成为无意义的负值。

作者
Author: 978142355 时间: 2016-12-24 09:07

wbn 发表于 2016-12-24 06:43
最好把这个帖子的内容也附在后面... 我用ghost原子做了半天结果bsse全是负的...
http://bbs.keinsci.com/t ...

嗯嗯，谢谢提出意见。

作者
Author: xrpeng 时间: 2023-1-1 16:59
请问如果分别在有无PCM下优化结构，这里提到的“相对可行的做法是使用气相下得到的BSSE校正能，加到溶剂下相互作用能上去”那么水相的BSSE使用哪个结构计算呢

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