计算化学公社

标题: gaussian09中pm6几何优化260个原子体系长时间不结束的问题 [打印本页]

作者
Author: renyutong 时间: 2023-9-20 13:53
标题: gaussian09中pm6几何优化260个原子体系长时间不结束的问题
最近在gaussian09中进行pm6的初步的几何优化，一共260个原子构成的分子。在输入#p opt pm6时出现l103错误，改正输入#p opt=cartesian pm6出现l502错误，其中参考帖子也尝试了帮助收敛的方法，我猜测是 opt=cartesian会使收敛困难。于是我输入#p opt=cartesian freq pm6 SCF=conver=6运算到势能曲线发生震荡时（图一）停止运算后，保存结构改用默认opt继续运算，发现二十小时不结束，势能曲线还算平滑但在初始有很大波动并且没有报错（图二），请问这种是应该继续等呢还是尝试新的帮助收敛的方法？
(, 下载次数 Times of downloads: 4) （图1） (, 下载次数 Times of downloads: 3) （图二）

作者
Author: Freeman 时间: 2023-9-20 14:05
挑出来能量最低点，以其为初始结构，继续算即可。

作者
Author: renyutong 时间: 2023-9-20 14:23

Freeman 发表于 2023-9-20 14:05
挑出来能量最低点，以其为初始结构，继续算即可。

谢谢您的回答，我观察了一会，现在开始有轻微的震荡，也同样挑出最低点算吗？

作者
Author: Freeman 时间: 2023-9-20 19:03

renyutong 发表于 2023-9-20 14:23
谢谢您的回答，我观察了一会，现在开始有轻微的震荡，也同样挑出最低点算吗？

我的经验是，一旦开始震荡了，挑最低点算更容易收敛。也有不管它，它最后收敛的情况，但不多，且一般要等很久，不值得。

作者
Author: renyutong 时间: 2023-9-20 20:56

Freeman 发表于 2023-9-20 19:03
我的经验是，一旦开始震荡了，挑最低点算更容易收敛。也有不管它，它最后收敛的情况，但不多，且一般要等 ...

好的，我一共有三步优化，那我是继续使用pm6至完全收敛还是进行下一步#p opt wb97xd/lanlzdz继续优化呢？

作者
Author: sobereva 时间: 2023-9-20 22:23
你当前的计算最大的硬伤就是用了PM6。怎么解决不收敛那都是后事
当前的体系有非常显著的分子内色散作用，用没法描述色散作用的PM6优化的结果毫无意义。注意看下文。必须用Gaussian的话，至少用PM6D3或PM7。如果用MOPAC程序，可以用更好的PM6-D3H4X优化；如果用xtb程序，可以用更好的GFN2-xTB优化。

谈谈“计算时是否需要加DFT-D3色散校正？”
http://sobereva.com/413（http://bbs.keinsci.com/thread-9772-1-1.html）

作者
Author: renyutong 时间: 2023-9-21 16:10

sobereva 发表于 2023-9-20 22:23
你当前的计算最大的硬伤就是用了PM6。怎么解决不收敛那都是后事
当前的体系有非常显著的分子内色散作用， ...

谢谢老师，我这就修改

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