求助：表面模型优化不收敛的问题

Kklove · 发表于 Post on 2021-4-14 17:04:45

请问大家有没有遇到过这种情况，同样的参数，优化Fe单胞可以，但是建立Fe的超胞表面模型之后，优化一直不收敛，能量一直震荡，不知道是什么原因。

sobereva · 发表于 Post on 2021-4-14 17:43:59

如果之前用的BFGS，试试CG没有更多信息别人不好回答

Kklove · 发表于 Post on 2021-4-14 18:01:54

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

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

sobereva · 发表于 Post on 2021-4-15 07:35:25

周期性设置不当。应该在Multiwfn里设成二维周期性体系，这样真空层还能设得小点，计算能快不少
MULTIPLICITY 5这个设置明显不合理，单个铁原子自旋多重度才是5

Kklove · 发表于 Post on 2021-4-15 08:48:12

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

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

sobereva · 发表于 Post on 2021-4-16 12:59:32

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

Kklove · 发表于 Post on 2021-4-16 19:09:00

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

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

sobereva · 发表于 Post on 2021-4-21 15:31:43

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

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

Kklove · 发表于 Post on 2021-4-22 12:12:37

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

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

sobereva · 发表于 Post on 2021-4-23 00:02:55

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

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

Kklove · 发表于 Post on 2021-4-25 09:12:09

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

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

sobereva · 发表于 Post on 2021-4-26 12:46:14

优化晶胞没必要用超胞，用原胞考虑k点效率高得多
给你我跑的Fe原胞算单点的例子

Fe_bulk.inp (2.38 KB, 下载次数 Times of downloads: 342)

Fe_bulk.out (119.05 KB, 下载次数 Times of downloads: 90)

Kklove · 发表于 Post on 2021-4-26 14:40:43

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

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

不会扣篮的后卫 · 发表于 Post on 2022-1-25 13:01:51

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，有什么用途吗?

sobereva · 发表于 Post on 2022-1-26 01:29:28

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 求助：表面模型优化不收敛的问题