计算化学公社

标题: CuO表面模型SCF一直不收敛什么原因？ [打印本页]

作者
Author: yaoyuan0711 时间: yesterday 10:15
标题: CuO表面模型SCF一直不收敛什么原因？
各位老师好，我构建了一个CuO(OO2)的表面模型，用PBE/DZVP-MOLOPT-SR-GTH优化，截断能450/55，电荷0，自旋多重度1，k点3 2 1，固定最下面两层原子，SCF change一直在个位数徘徊，没有下降的趋势，Mixing method改为Pulay_mixing也没用，基组改为pob-TZVP更差了，SCF change一直在两位数了。 (, 下载次数 Times of downloads: 8)

作者
Author: sobereva 时间: yesterday 11:08
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

作者
Author: yaoyuan0711 时间: yesterday 13:35

sobereva 发表于 2024-11-26 11:08
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.ht ...

“即便体系没有磁性，利用&BS使初猜波函数中的原子带电状态使之尽可能符合实际也往往是有益的。比如对于NaCl这样高度离子性体系计算时SCF不收敛的话，可以尝试令初猜中Na和Cl的原子组态分别对应Na+和Cl-的状态（默认情况下分别对应中性的Na和Cl）。曾经我用PBEsol/pob-TZVP-rev2计算10层的NaCl板时通过这种做法成功解决了SCF不收敛问题。”

请问sob老师，针对我这个CuO该如何设置磁矩？Cu2+最外层3d9，O2-最外层2p6。

作者
Author: sobereva 时间: yesterday 13:39

yaoyuan0711 发表于 2024-11-26 13:35
“即便体系没有磁性，利用&BS使初猜波函数中的原子带电状态使之尽可能符合实际也往往是有益的。比如对于N ...

CuO早就有无数文章算过（不限于体相，slab应该也有），看文中报道的自旋磁矩，恰当用MAGNETIZATION或&BS设置。

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