计算化学公社

标题: Gaussian分子筛结构优化，不能收敛 [打印本页]

作者
Author: youfish 时间: 2017-6-13 23:09
标题: Gaussian分子筛结构优化，不能收敛
新手求助！~
参照文献建立了一个Fe-zsm-5的分子筛分子模型，采用Lanl2DZ描述Fe，其余原子采用b3lyp描述。计算过程中出现如下报错，请各位多多指教~

Annihilation of the first spin contaminant:
S**2 beforeannihilation    3.7803, after 3.7503
Convergence failure -- runterminated.
Error termination via Lnk1e in/home/jiang337/g09/l502.exe at Tue Jun 13 20:22:58 2017.
Job cpu time:  0 days 12 hours 17 minutes 54.6seconds.
File lengths (MBytes):  RWF= 144 Int=    0 D2E=    0Chk=    5 Scr=    1

作者
Author: youfish 时间: 2017-6-13 23:31
补充构建的模型，其中粉色是Al，蓝色是Fe，灰色是Si，白色是H
输入文件的语句是# opt freq b3lyp/genecp geom=connectivity

作者
Author: sobereva 时间: 2017-6-14 00:19
解决SCF不收敛问题的方法
http://sobereva.com/61

“采用Lanl2DZ描述Fe，其余原子采用b3lyp描述”是错误的描述，一个说的基组一个说的理论方法
这种奇形怪状的体系，再加上过渡金属，SCF收敛肯定不容易

作者
Author: youfish 时间: 2017-6-14 10:12

sobereva 发表于 2017-6-14 00:19
解决SCF不收敛问题的方法
http://sobereva.com/61

感谢sob老师解答！
Fe使用Lanl2DZ，其余原子使用6-31+g(d,p)，今后会注意表述的严谨性~
至于结构的话，是从zsm-5的cif中截取出的一个团簇，奇怪也没办法T T。所有的H以及与它直接相连的Si和O都固定，保证分子筛定性，其余保持松弛。
我再参照您说的办法试试。

作者
Author: zyj19831206 时间: 2017-6-15 12:43

youfish 发表于 2017-6-14 10:12
感谢sob老师解答！
Fe使用Lanl2DZ，其余原子使用6-31+g(d,p)，今后会注意表述的严谨性~
至于结构的话， ...

gaussian收敛效果应该比dmol要好吧？

作者
Author: sobereva 时间: 2017-6-16 01:08

zyj19831206 发表于 2017-6-15 12:43
gaussian收敛效果应该比dmol要好吧？

dmol3我不常用，应该是的，至少gaussian里对解决不收敛可调选项更多

作者
Author: woshiblake 时间: 2023-2-20 21:29
您好，我也碰到这个问题，后来怎么解决的？

作者
Author: 绫波レイ 时间: 2023-12-14 16:34

woshiblake 发表于 2023-2-20 21:29
您好，我也碰到这个问题，后来怎么解决的？

不是自旋污染造成的，就是SCF不收敛而已，看《解决SCF不收敛问题的方法》http://sobereva.com/61

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