计算化学公社

标题: 求助：用高斯说明书中没有的溶剂优化体系不收敛 [打印本页]

作者
Author: 哈哈宝 时间: 2015-5-4 13:54
标题: 求助：用高斯说明书中没有的溶剂优化体系不收敛
在不加溶剂的时候，优化正常收敛；审稿人认为我的体系含有3个负电荷，只气态优化，计算的占据轨道是正值，这是不对的，所以，要求修改。我参看了一些类似的文章，其中加入了溶剂化效应，这样就能使得占据轨道本征值为负值，被认可。我引入乙二胺溶剂，但是该溶剂高斯说明书中没有，只查到了eps=12.9；另外查到epsinf是折射率的平方，就在网上搜到乙二胺的折射率，大致得到epsinf=2.114;并使用高斯的volume关键词，计算了波尔半径为3.74。
命令行：#p opt pop=full pbe1pbe/def2tzvp scrf(SMD,Solvent=Generic,Read)
结构行
空一行
eps=12.9
epsinf=2.114
RSolv=3.74
空两行

开始优化，结果显示：CheckU: Warning! The PCM matrix U/Ut may be ill-conditioned!而且不能正常收敛！
请问，这是我输入的溶剂不能被高斯读取么？用的是G09 B.01。谢谢高手解答！

作者
Author: 哈哈宝 时间: 2015-5-4 15:04
没有人遇到这种问题么？拼了！

加入了opt=（tight,calcall)，看看这样行不行，不知道多久才能跑一圈！

作者
Author: sobereva 时间: 2015-5-5 01:21
你只需要定义eps就够了，定义其它两个对你当前的计算没意义。epsinf仅对激发态等任务有意义。

至于解决报错，可以先尝试改成其它溶剂模型，如CPCM等，或者尝试微调一下初始结构，或者先稍微把eps改成其它数值等，需要多尝试，此报错属于数值巧合。

如果一开始就出错，tight、calcall不会解决任何问题。

作者
Author: 哈哈宝 时间: 2015-5-5 15:12
谢谢Sob老师用了CPCM优化正常结束了

作者
Author: Qing 时间: 2021-2-13 13:45

sobereva 发表于 2015-5-5 01:21
你只需要定义eps就够了，定义其它两个对你当前的计算没意义。epsinf仅对激发态等任务有意义。

至于解决 ...

老师，根据您《谈谈隐式溶剂模型下溶解自由能和体系自由能的计算》提到的观点，现在是不是不应该用CPCM模型了，而应该使用IEFPCM模型呢

作者
Author: biogon 时间: 2021-2-13 14:04

Qing 发表于 2021-2-13 13:45
老师，根据您《谈谈隐式溶剂模型下溶解自由能和体系自由能的计算》提到的观点，现在是不是不应该用CPCM模 ...

是

作者
Author: 喵星大佬 时间: 2021-2-14 07:57
不要用smd去优化啊。。。。。
日常各种莫名其妙的问题

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