计算化学公社

标题: 求助：表面模型优化不收敛的问题 [打印本页]

作者
Author: Kklove 时间: 2021-4-14 17:04
标题: 求助：表面模型优化不收敛的问题
请问大家有没有遇到过这种情况，同样的参数，优化Fe单胞可以，但是建立Fe的超胞表面模型之后，优化一直不收敛，能量一直震荡，不知道是什么原因。

作者
Author: sobereva 时间: 2021-4-14 17:43
如果之前用的BFGS，试试CG没有更多信息别人不好回答

作者
Author: Kklove 时间: 2021-4-14 18:01

sobereva 发表于 2021-4-14 17:43
如果之前用的BFGS，试试CG没有更多信息别人不好回答

这是我的输入文件，老师您看看哪有问题？谢谢您

作者
Author: sobereva 时间: 2021-4-15 07:35
周期性设置不当。应该在Multiwfn里设成二维周期性体系，这样真空层还能设得小点，计算能快不少
MULTIPLICITY 5这个设置明显不合理，单个铁原子自旋多重度才是5

作者
Author: Kklove 时间: 2021-4-15 08:48

sobereva 发表于 2021-4-15 07:35
周期性设置不当。应该在Multiwfn里设成二维周期性体系，这样真空层还能设得小点，计算能快不少
MULTIPLICI ...

太谢谢您了，就是自旋多重度这个问题很头疼，我的理论功底不够，目前我知道单个铁原子的多重度，但是对于这种建超胞之后表面模型多重度的判断还不是很懂，不知道是通过什么方式来计算的，希望老师能再指导一下，谢谢您。

作者
Author: sobereva 时间: 2021-4-16 12:59
算这个得开smearing，此时默认不会保持alpha和beta电子数的差值（除非结合FIXED_MAGNETIC_MOMENT设置，此时不兼容k点），因此设置自旋多重度已经没有意义了，设不设都一样。
关于孤立原子和体相Fe带单电子情况的差异，建议你看看https://www.fzu.cz/~sipr/posters/hand_on_sprkkr_05.pdf

作者
Author: Kklove 时间: 2021-4-16 19:09

sobereva 发表于 2021-4-16 12:59
算这个得开smearing，此时默认不会保持alpha和beta电子数的差值（除非结合FIXED_MAGNETIC_MOMENT设置，此时 ...

太好了，谢谢您的讲解。那我这个体系从最开始就不用设定自旋多重度，只需要设置LSD和smearing吗？

作者
Author: sobereva 时间: 2021-4-21 15:31

Kklove 发表于 2021-4-16 19:09
太好了，谢谢您的讲解。那我这个体系从最开始就不用设定自旋多重度，只需要设置LSD和smearing吗？

是，但还应该同时设置原子的MAGNETIZATION（大于0），使之在初猜波函数中就是带d单电子的

作者
Author: Kklove 时间: 2021-4-22 12:12

sobereva 发表于 2021-4-21 15:31
是，但还应该同时设置原子的MAGNETIZATION（大于0），使之在初猜波函数中就是带d单电子的

谢谢您的回复，但是我添加了MAGNETIZATION值为1后，超胞的计算还是不收敛，在第一步结束后会有“ WARNING: S**2 computation does not yet treat fractional occupied orbitals”的警告，希望您能再指导一下。

作者
Author: sobereva 时间: 2021-4-23 00:02

Kklove 发表于 2021-4-22 12:12
谢谢您的回复，但是我添加了MAGNETIZATION值为1后，超胞的计算还是不收敛，在第一步结束后会有“ WARNING ...

MAGNETIZATION还可以设更大。本来体相和表面的铁的磁矩就都比1大
warning不用管

作者
Author: Kklove 时间: 2021-4-25 09:12

sobereva 发表于 2021-4-23 00:02
MAGNETIZATION还可以设更大。本来体相和表面的铁的磁矩就都比1大
warning不用管

谢谢您的回复，修改了几个磁矩（2,3,4），还是不收敛，能量一直在震荡，不知道是不是别的参数有问题，麻烦老师再帮我看一下。

作者
Author: sobereva 时间: 2021-4-26 12:46
优化晶胞没必要用超胞，用原胞考虑k点效率高得多
给你我跑的Fe原胞算单点的例子
(, 下载次数 Times of downloads: 340)
(, 下载次数 Times of downloads: 90)

作者
Author: Kklove 时间: 2021-4-26 14:40

sobereva 发表于 2021-4-26 12:46
优化晶胞没必要用超胞，用原胞考虑k点效率高得多
给你我跑的Fe原胞算单点的例子

好的老师，我优化用的就是单胞，主要是建表面模型要用到超胞，可是超胞的优化一直不收敛。算单点的文件您看看有啥问题没。

作者
Author: 不会扣篮的后卫 时间: 2022-1-25 13:01

sobereva 发表于 2021-4-26 12:46
优化晶胞没必要用超胞，用原胞考虑k点效率高得多
给你我跑的Fe原胞算单点的例子

卢老师，请问在Fe_bulk.inp 中。这一行命令（来自倒数第4行），ADDED_MOS 30 30 #Number of virtual MOs to be solved for alpha and beta spins，有什么用途吗?

作者
Author: sobereva 时间: 2022-1-26 01:29

不会扣篮的后卫发表于 2022-1-25 13:01
卢老师，请问在Fe_bulk.inp 中。这一行命令（来自倒数第4行），ADDED_MOS 30 30 #Number of virtual MOs ...

用了smearing就必须求解最低的一部分空轨道，知道smearing的原理自然就明白了

作者
Author: sobereva 时间: 2022-1-26 01:30

Kklove 发表于 2021-4-26 14:40
好的老师，我优化用的就是单胞，主要是建表面模型要用到超胞，可是超胞的优化一直不收敛。算单点的文件您 ...

记得设MULTIPLICITY

作者
Author: 不会扣篮的后卫 时间: 2022-1-26 01:39

sobereva 发表于 2022-1-26 01:29
用了smearing就必须求解最低的一部分空轨道，知道smearing的原理自然就明白了

好的，谢谢老师！

作者
Author: 不会扣篮的后卫 时间: 2022-2-8 11:17

sobereva 发表于 2021-4-23 00:02
MAGNETIZATION还可以设更大。本来体相和表面的铁的磁矩就都比1大
warning不用管

卢老师，请问为什么不用管这个“ WARNING: S**2 computation does not yet treat fractional occupied orbitals”？

作者
Author: sobereva 时间: 2022-2-8 23:44

不会扣篮的后卫发表于 2022-2-8 11:17
卢老师，请问为什么不用管这个“ WARNING: S**2 computation does not yet treat fractional occupied or ...

S**2主要用来判断自旋污染，而自旋污染对DFT波函数本来就没有明确意义

作者
Author: 卡开发发 时间: 2022-2-9 05:37

不会扣篮的后卫发表于 2022-2-8 11:17
卢老师，请问为什么不用管这个“ WARNING: S**2 computation does not yet treat fractional occupied or ...

补充一点，引入展宽之后，本身就是把电子对轨道的占据按照统计来做的，这个基础上讨论自旋污染没什么价值。

作者
Author: 不会扣篮的后卫 时间: 2022-2-9 09:20

sobereva 发表于 2022-2-8 23:44
S**2主要用来判断自旋污染，而自旋污染对DFT波函数本来就没有明确意义

谢谢卢老师！

作者
Author: 不会扣篮的后卫 时间: 2022-2-9 09:21

卡开发发发表于 2022-2-9 05:37
补充一点，引入展宽之后，本身就是把电子对轨道的占据按照统计来做的，这个基础上讨论自旋污染没什么价值 ...

谢谢老师！

作者
Author: 不会扣篮的后卫 时间: 2022-2-13 12:53

sobereva 发表于 2021-4-16 12:59
算这个得开smearing，此时默认不会保持alpha和beta电子数的差值（除非结合FIXED_MAGNETIC_MOMENT设置，此时 ...

请问卢老师，您这个回复是说，如果在CP2K输入文件里加了磁矩参数和FIXED_MAGNETIC_MOMENT设置，
是不是就不用设置自旋多重度了？

作者
Author: sobereva 时间: 2022-2-13 13:42

不会扣篮的后卫发表于 2022-2-13 12:53
请问卢老师，您这个回复是说，如果在CP2K输入文件里加了磁矩参数和FIXED_MAGNETIC_MOMENT设置，
是不是 ...

应该是

作者
Author: 不会扣篮的后卫 时间: 2022-2-13 14:03

sobereva 发表于 2022-2-13 13:42
应该是

好的，谢谢卢老师！

作者
Author: wth1219 时间: 2022-4-9 21:20
抱歉借题咨询一个相关问题。我想用CP2K跑一个含有约10个Fe(2+)的体系的FPMD，已知体系中的Fe(2+)都是高自旋的、体系总体显顺磁性（除Fe(2+)外无磁性粒子。）。想向各位老师请教一下如何确定体系的自旋多重度、如何在FPMD过程中维持Fe(2+)的高自旋、以及在&KIND中高自旋的Fe(2+)的MAGNETIZATION是5吗？谢谢！

作者
Author: Yong过儿 时间: 2023-5-12 17:25
请问楼主问题解决了吗？我最近也遇到了Fe超胞无法收敛的情况，希望您指导一下

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