计算化学公社

标题: 开smearing后UKS与RKS对于同一结构的处理的问题？ [打印本页]

作者
Author: LIJINWEI 时间: 2025-9-10 14:35
标题: 开smearing后UKS与RKS对于同一结构的处理的问题？
以下均开smearing；
例如，现在有一个srvo3结构，该结构中V是+4价，电子组态是3d1，最外层有一个电子，按照正常的做法，自旋多重度设置为2，UKS开自旋极化计算，在开smearing，可以看见scf迭代收敛的很快，最后得出的能量：-150.68471951565894。

此时，在同样开smearing时，将自旋多重度设置为1，不在开启UKS计算，可以看见scf迭代的更快，最后得出的能量：-150.68400289626967。
我想问的问题是：
1.在不开smearing时，只有自旋多重度为2时，才不会报错，但是开启smearing后，会发现自旋多重度无论是2或者1均不会报错，这是什么原因。
2.针对第一个问题在开smearing时，使用UKS计算和不使用UKS计算两者差别大吗？其实按照结构分析来看，自旋多重度设置为2更合理。
3.我在最初自旋多重度设置为2时，同时开启smearing后，得到了优化后的结构，针对这个结构进行了3*3*3的扩胞，再次同样的进行了结构优化，对于这个优化的结构往里面的八面体掺杂了一个中性原子，对于这个掺杂后的结构出现了scf震荡的结果，我在进行多方面调整参数和测试之后，发现其中一条可行的路。
对于这个扩胞后掺杂的结构，即便分析结构后，得出的结论是含有磁性的，即自旋多重度是大于1的，此时应该开UKS计算，最后是出现scf迭代震荡，当我把自旋多重度设置为1时，不开UKS计算，突然发现scf就收敛了，而且还非常快。
所以此时，我面临一个问题，也有这样一个想法，就是即便这个结构分析之后应该是含有磁性的，也把这个自旋多重度设置为1，同时关闭UKS计算，这样做得出的结果会不会是错误的。
如果按照上面所做的，那么我将从原胞开始重新按照没有磁性开始结构优化，这样做是否可行？（可能是我自身的问题，在开磁性之后，scf迭代问题就一直存在，即便我分析了多次，认为磁性应该是没有错误的，但在cp2k计算中一直碰到scf迭代震荡问题，从原胞开始到扩胞，虽然碰到了这样的问题，但靠着社长的《CP2K中遇到SCF难收敛时的解决方法http://bbs.keinsci.com/forum.php ... 37196&fromuid=69037(出处: 计算化学公社)》，也大多调整过来，直到碰到掺杂之后，即便怎么调整也是出现scf迭代震荡问题，我将掺杂的输入文件放在附件上面了，十分麻烦大家。）

作者
Author: sobereva 时间: 2025-9-10 19:42
不要一味地管报错不报错，要搞清楚计算原理，不当的设置下就算不报错结果也毫无意义。解决SCF不收敛问题必须在正确的设置的前提下进行。磁性体系用限制性闭壳层计算根本没意义，显然必须写UKS

作者
Author: LIJINWEI 时间: 2025-9-10 20:13
本帖最后由 LIJINWEI 于 2025-9-10 20:22 编辑

sobereva 发表于 2025-9-10 19:42
不要一味地管报错不报错，要搞清楚计算原理，不当的设置下就算不报错结果也毫无意义。解决SCF不收敛问题必 ...

社长，明白了，多谢您。

作者
Author: LIJINWEI 时间: 2025-9-10 20:30

sobereva 发表于 2025-9-10 19:42
不要一味地管报错不报错，要搞清楚计算原理，不当的设置下就算不报错结果也毫无意义。解决SCF不收敛问题必 ...

社长，在电子组态中可以认定v得到0.5个电子吗？例如zn的两个电子，给了周围四个v每个0.5个电子，还是只需要将其中两个定义为+3价就行。

作者
Author: sobereva 时间: 2025-9-10 22:28

LIJINWEI 发表于 2025-9-10 20:30
社长，在电子组态中可以认定v得到0.5个电子吗？例如zn的两个电子，给了周围四个v每个0.5个电子，还是只需 ...

&BS或magnetization设置只影响初猜波函数，SCF过程中该怎么变还怎么变。关注这个没意义

作者
Author: LIJINWEI 时间: 2025-9-11 09:46

sobereva 发表于 2025-9-10 22:28
&BS或magnetization设置只影响初猜波函数，SCF过程中该怎么变还怎么变。关注这个没意义

社长，我参考了您的CP2K中《遇到SCF难收敛时的解决方法http://bbs.keinsci.com/forum.php ... 37196&fromuid=69037(出处: 计算化学公社)》，能调的大概都调了，我试着smearing加温度，发现scf迭代依然震荡，其他调整完之后，也变化不大（除开时间变长之外），在这些不行后，我又决定重新来看结构，最初我是认定zn是中性原子，不与周围反应，我现在想着，zn失去两个电子，同时有两个v从+4价变成+3价，按照现在重新认定的结构试试看，可是，社长按照您的说法，这样不行吗？

作者
Author: sobereva 时间: 2025-9-11 14:44

LIJINWEI 发表于 2025-9-11 09:46
社长，我参考了您的CP2K中《遇到SCF难收敛时的解决方法http://bbs.keinsci.com/forum.php ... 37196&from ...

先弄清楚你改的到底是什么，别光一个劲盲目地说价态什么的。并且看下文了解价态到底是个什么
使用Multiwfn结合CP2K计算晶体中原子的氧化态
http://sobereva.com/711（http://bbs.keinsci.com/thread-45876-1-1.html）

我前面已经说了，&BS只影响初猜波函数。如果你是试图让初猜波函数的电子分布方式和收敛的波函数更接近从而试图让SCF收敛更容易点，那你就试

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