计算化学公社

标题: 508错误 [打印本页]

作者
Author: lei234 时间: 2016-3-8 02:31
标题: 508错误
错误信息如下：
Density matrix breaks symmetry, PCut= 1.00D-04
Rerun with SCF=IntRep.
Error termination via Lnk1e in /opt/g09.d01/g09/l508.exe at Sat Mar 5 01:12:07 2016.
输入文件已附上，主要关键词如下：
#p opt=(cartesian,gdiis) freq PBEPBE/3-21G scf=(xqc,maxconventional=100)

望指点一二~

作者
Author: liyuanhe211 时间: 2016-3-8 03:06
3-21G算出来的没有意义，最少6-31G(d)。
常规体系一般用不着xqc，用其他方法经常都能解决。
另外有个不太确定的建议，听人说类似结构可能有单电子存在（自旋极化），不知道提供一个对称性破缺的初猜有没有好处？

作者
Author: sobereva 时间: 2016-3-8 04:30
提示已经写了，加上SCF=IntRep重算
没事别用xqc
opt里cartesian若无必要应当去掉，否则会需要更多步数才能收敛

作者
Author: 卡开发发 时间: 2016-3-8 06:39
体系确实有点大，Gaussian用8个核心跑估计也就3-21G了，如果没啥特定的针对性，考虑PBC可能会容易一些。

SCF难收敛的原因我觉得可能会是前线轨道太密集导致频繁交换，如果是这样可以考虑vshift，把PBE换成HSE也有可能能够解决。

作者
Author: lei234 时间: 2016-3-8 08:14

liyuanhe211 发表于 2016-3-8 03:06
3-21G算出来的没有意义，最少6-31G(d)。
常规体系一般用不着xqc，用其他方法经常都能解决。
另外有个不太 ...

您说的提供一个对称性破缺的初猜要怎么做？加nosymm吗？

作者
Author: liyuanhe211 时间: 2016-3-8 09:57
3-21G算出来的还没有某些力场或者半经验的效果好，几乎没有任何意义。
先试Sob说的。

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