508错误

lei234

错误信息如下：
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)

望指点一二~

liyuanhe211

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

sobereva

提示已经写了，加上SCF=IntRep重算
没事别用xqc
opt里cartesian若无必要应当去掉，否则会需要更多步数才能收敛

卡开发发

体系确实有点大，Gaussian用8个核心跑估计也就3-21G了，如果没啥特定的针对性，考虑PBC可能会容易一些。

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

lei234

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

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

liyuanhe211

3-21G算出来的还没有某些力场或者半经验的效果好，几乎没有任何意义。
先试Sob说的。