计算化学公社

标题: 适合开壳层体系的方法和基组 [打印本页]

作者
Author: diaolanxinyu 时间: 2015-2-1 16:23
标题: 适合开壳层体系的方法和基组
请问适合开壳层体系的方法和基组有哪些？哪些是可以避免自旋污染的？谢谢大家！

我用b3lyp/6-31+g(d)扫描势能面的时候出现了自旋污染
Annihilation of the first spin contaminant:
S**2 before annihilation 1.5799, after 0.7574
Convergence failure -- run terminated.

加入IOP（5/13=1）虽然可以正常结束，但是不知道结果还可靠吗？

作者
Author: sobereva 时间: 2015-2-1 16:43
IOP(5/13=1)任何时候都不要用。
UDFT、耦合簇、MCSCF计算开壳层都适合。
不用管DFT的自旋污染，自旋污染这个概念对于DFT本身就是不明确的。耦合簇受自旋污染小，MCSCF或CI是基于组态函数的，没有自旋污染问题。

作者
Author: diaolanxinyu 时间: 2015-2-1 17:09

sobereva 发表于 2015-2-1 16:43
IOP(5/13=1)任何时候都不要用。
UDFT、耦合簇、MCSCF计算开壳层都适合。
不用管DFT的自旋污染，自旋污染 ...

”不用管DFT的自旋污染“，那就是我选的基组不合适唠？

那适合开壳层体系的基组有哪些呢？

谢谢sobereva老师！

作者
Author: sobereva 时间: 2015-2-1 17:26
基组和开壳层闭壳层无关。
不用管DFT的自旋污染是理论方法的缘故，和基组也无关。

作者
Author: diaolanxinyu 时间: 2015-2-1 18:53
谢谢sobereva老师！

作者
Author: ter20 时间: 2015-2-2 09:29
目测你算的是单自由基(maybe 0 2)，自旋污染程序已经annihilate到了0.75左右，没关系的
S**2 before annihilation 1.5799, after 0.7574

作者
Author: helpme 时间: 2015-2-2 09:40
这个你要试试多重度是4的情况了吧，看和多重度2相比，哪个能量低。

作者
Author: diaolanxinyu 时间: 2015-2-2 11:07

ter20 发表于 2015-2-2 09:29
目测你算的是单自由基(maybe 0 2)，自旋污染程序已经annihilate到了0.75左右，没关系的
S**2 before annih ...

我测试的是最简单的氧负离子自由基，电荷和自旋多重度（-1，2），用scan扫描时出现了自旋污染的错误，计算中断了；后来加了iop(5/13=1)可以正常结束了，但是sobereva老师说最好不要加，后来我修改了scan的步长就正常结束了，呵呵

感谢你参与讨论！

作者
Author: diaolanxinyu 时间: 2015-2-2 11:09

helpme 发表于 2015-2-2 09:40
这个你要试试多重度是4的情况了吧，看和多重度2相比，哪个能量低。

感谢你的建议，如果需要的话我会测试比较一下的。

作者
Author: liuzl 时间: 2021-10-26 09:30

diaolanxinyu 发表于 2015-2-2 11:07
我测试的是最简单的氧负离子自由基，电荷和自旋多重度（-1，2），用scan扫描时出现了自旋污染的错误，计 ...

您好我也遇到了和您一样的问题，开壳层体系计算nbo时出现了自旋污染，我是初学者，不太明白您说的修改scan步长是什么意思，请问您具体是怎么做的呢？

作者
Author: wzkchem5 时间: 2021-10-26 16:02

liuzl 发表于 2021-10-26 02:30
您好我也遇到了和您一样的问题，开壳层体系计算nbo时出现了自旋污染，我是初学者，不太明白您说的修改sca ...

不要受误导，自旋污染不会直接导致报错。有自旋污染的体系倾向于（但不必然）SCF不收敛，而SCF不收敛会导致报错，自旋污染不是导致报错的直接原因。如果只有自旋污染，但计算没有中断，那么无需仅因为有自旋污染就说结果不可靠。
修改scan步长是解决扫描任务的SCF不收敛问题（不管这个不收敛是不是因为自旋污染导致的都可以）的方法之一（但不是首选方法，也不保证奏效），不能解决其他类型任务的SCF不收敛问题，更不能解决自旋污染问题本身。
仔细阅读http://sobereva.com/264，彻底了解什么叫自旋污染再来提问

作者
Author: liuzl 时间: 2021-10-26 19:44

wzkchem5 发表于 2021-10-26 16:02
不要受误导，自旋污染不会直接导致报错。有自旋污染的体系倾向于（但不必然）SCF不收敛，而SCF不收敛会导 ...

请问如何修改scan步长？具体操作是什么，加什么关键词

作者
Author: zjxitcc 时间: 2021-10-26 19:53
本帖最后由 zjxitcc 于 2021-10-26 19:55 编辑

liuzl 发表于 2021-10-26 19:44
请问如何修改scan步长？具体操作是什么，加什么关键词

这个帖子里说的scan有两种可能的意思：刚性扫描或柔性扫描，是势能面扫描的两种常见做法。至于是哪种，上述几层没有详细指出。

看这篇sob老师的博文《详谈使用Gaussian做势能面扫描》http://sobereva.com/474学习如何做势能面扫描。当然，如果你不做势能面扫描，自然也就与scan步长无关，无需关心scan步长。

如果你碰到SCF不收敛，应该看这篇《解决SCF不收敛问题的方法》http://sobereva.com/61。如果你还有其他问题，建议另开一贴详细展示，方便大家对症下药。

作者
Author: wzkchem5 时间: 2021-10-26 20:03

liuzl 发表于 2021-10-26 12:44
请问如何修改scan步长？具体操作是什么，加什么关键词

你还是没看懂我在说什么。应该先把别人的回帖看懂了再追问

作者
Author: liuzl 时间: 2021-10-26 20:21

zjxitcc 发表于 2021-10-26 19:53
这个帖子里说的scan有两种可能的意思：刚性扫描或柔性扫描，是势能面扫描的两种常见做法。至于是哪种， ...

好的谢谢

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