nusiew 发表于 2022-9-27 19:28 给你个计算体相Fe的例子 |
2.38 KB, 下载次数 Times of downloads: 83
119.05 KB, 下载次数 Times of downloads: 18
nusiew 发表于 2022-9-27 21:13 1 显然是 2 不加UKS会当成闭壳层算，结果没意义 3 MAGNETIZATION可以是非整数。没必要非得设得那么凑巧，初始自旋设得大一点也没关系 |
本帖最后由 nusiew 于 2022-9-28 09:03 编辑 丁越 发表于 2022-9-28 08:49 非常感谢指导！我发现论坛里关于过渡金属及其氧化物SCF计算不收敛的问题特别多，但很少给出成熟的解决方案。 由于上述测试是在扩胞3*3*3的 Γ 点计算的，发现不收敛；于是我将扩胞改成使用k点，发现在相同的参数设置下实现了收敛。还有一种方法，是将Mixing中的alpha参数由0.4改成0.1，也能够在扩胞不使用k点的情况下收敛。如果振荡不收敛，在模型、smear、开壳层等设置合理的情况下，可能调整Mixing参数是一个不错的选择。 |
nusiew 发表于 2022-9-27 19:28 sob老师在论坛里解答过这个问题（http://bbs.keinsci.com/thread-22642-1-1.html） CP2K对一些电子结构复杂点的过渡金属体系计算还是存在问题的，比如之前论坛里有个人想用CP2K优化一个Ni的slab就出现SCF无法收敛的情况，这个属于程序方面的问题。这种情况下换成平面波程序如QE、VASP还是比较稳健的 |
本帖最后由 nusiew 于 2022-9-27 20:00 编辑 丁越 发表于 2022-9-26 22:42 非常感谢丁老师的指导！按照您的建议，我开了smear和UKS，并设置了MAGNETIZATION为3，同时注销了MULTIPLICITY（参考设置了smear不用指定自旋多重度），但是Fe超胞体系（3*3*3）SCF出现了不收敛的情况。于是我将注销的MULTIPLICITY 27 恢复后，但同时还保留smear、UKS、和MAGNETIZATION，发现SCF在62步后实现了收敛。现在有一些迷茫，添加smear不指定MULTIPLICITY，体系出现不收敛的情况；指定了MULTIPLICITY，但却不一定是合理的MULTIPLICITY。请问丁老师接下来该如何调整？恳请丁老师予以指导。我的inp文件见附件。 |
3.2 KB, 下载次数 Times of downloads: 23
| 常见的体相材料在material project上都有，可以直接查看上面的SCF收敛得到的磁矩信息，可以参考着那个设，就省得你自己试各种设置了 |
nusiew 发表于 2022-9-26 21:31 1.是，在Multiwfn中就有设置选项。单个Fe原子磁矩为4μ_B，表面相要小于4，体相远小于4。设置比如2，3之类的都可以，程序会自发弛豫磁矩。 2.RELAX_MULTIPLICITY设置如0.001。看看https://groups.google.com/g/cp2k ... /kQ92btddAQAJ?pli=1 |
nusiew 发表于 2022-9-26 21:37 使用了smear时可以不用管自旋多重度。或者你还可以使用RELAX_MULTIPLICITY 对，不设置磁矩算出来的就不是该物质正确的电子态，所以对那些实空间波函数分析都有影响 |
本帖最后由 nusiew 于 2022-9-26 21:43 编辑 丁越 发表于 2022-9-26 19:53 谢谢丁老师的指导！在加了UKS并将MAGNETIZATION设为3，自旋多重度设为25，发现SCF在27步后实现了收敛。请问丁老师自旋多重度需要一个一个的测试吗？如果自旋多重度设置的不是对应能量最小值的，会对输出的PDOS等产生影响吗？对于类似含有Fe、Co、Ni等原子构成的磁性体系，在未正确设置磁矩和自旋多重度的条件下即使SCF收敛，输出的ELF、电子密度等实空间函数及其衍生描述符是不是都会存在很大误差？ |
Eudaimonia 发表于 2022-9-26 16:14 非常感谢指导。 （1）请问Fe的磁矩初猜是基于&Kind 中的 MAGNETIZATION参数来设置吗？该参数一般怎么取值呢？对于我这个Fe原子体系是直接用电子数（alpha-beta）=4 来取值吗？ （2）自动优化自旋多重度是直接加上 RELAX_MULTIPLICITY 这个关键词就可以了吗？请问该关键词的取值如何根据自己的体系来设定呢？ 再次感谢您的回复！ |
| 磁矩没设，UKS没加，自旋多重度不对，不是1。 |
|
大概率是因为磁矩带来的问题，建议对Fe磁矩进行初猜并且自动优化自选多重度，不会设置的话Multiwfn里面可以手动指定原子磁矩 顺便，扩胞555跑这个晶胞远远不如加K点来的划算 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-20 23:40 , Processed in 0.230832 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites