计算化学公社

 找回密码 Forget password
 注册 Register

求助:cp2k优化Fe110晶面SCF和优化不收敛

查看数: 11605 | 评论数: 17 | 收藏 Add to favorites 7
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2021-12-6 21:48

正文摘要:

本帖最后由 zx4562430 于 2021-12-6 21:56 编辑 各位老师,想请教个问题,谢谢。最近根据社长大大的博文,练习将cp2k应用于自己实际体系中,但优化这个Fe110晶面时,总是scf不收敛,结构优化优化1-2天都距离收敛 ...

回复 Reply

zzp 发表于 Post on yesterday 23:29
请问下楼主这个最后算收敛了吗
丁越 发表于 Post on 2022-10-21 20:11:40
时光的邂逅 发表于 2022-10-21 10:38
请问一下,如果按照老师您的意思,添加smearing之后,自旋多重度值可以随便设置,那么对于铁体系的话,是 ...

需要开启UKS
时光的邂逅 发表于 Post on 2022-10-21 10:38:26
丁越 发表于 2022-4-14 23:14
添加SMEAR后默认不会保持alpha和beta电子数差值,自旋多重度值可以设置比较随意。或者你直接设RELAX_MULT ...

请问一下,如果按照老师您的意思,添加smearing之后,自旋多重度值可以随便设置,那么对于铁体系的话,是否需要打开UKS或者LSD呢?
莫雨萱 发表于 Post on 2022-4-18 08:39:44
丁越 发表于 2022-4-14 23:14
添加SMEAR后默认不会保持alpha和beta电子数差值,自旋多重度值可以设置比较随意。或者你直接设RELAX_MULT ...

谢谢老师的回答!
莫雨萱 发表于 Post on 2022-4-18 08:39:11
zx4562430 发表于 2022-4-15 16:42
利用multiwfn构建cp2k输入文件时,利用其中6 Toggle smearing electron occupation, current: Yes,开启s ...

谢谢您的回答,非常感谢!
zx4562430 发表于 Post on 2022-4-15 16:42:20
莫雨萱 发表于 2022-4-14 18:47
楼主您好,请问对于Fe表面体系,多重度的设置您是如何理解的呀?多重度设置问题解决了吗?

利用multiwfn构建cp2k输入文件时,利用其中6 Toggle smearing electron occupation, current: Yes,开启smearing电子布居选项,这样默认不会保持alpha和beta电子数差值,因此设置自旋多重度已经没有意义了,设不设都一样。重新进行计算试试。
丁越 发表于 Post on 2022-4-14 23:14:37
莫雨萱 发表于 2022-4-14 18:47
楼主您好,请问对于Fe表面体系,多重度的设置您是如何理解的呀?多重度设置问题解决了吗?

添加SMEAR后默认不会保持alpha和beta电子数差值,自旋多重度值可以设置比较随意。或者你直接设RELAX_MULTIPLICITY 0.001使其自发的弛豫自旋多重度
莫雨萱 发表于 Post on 2022-4-14 18:47:40
zx4562430 发表于 2021-12-8 17:46
您好,请教您下,根据您分享的pdf和multiwfn实际磁矩设置操作,我的理解是磁矩设置应该是针对每个原子的 ...

楼主您好,请问对于Fe表面体系,多重度的设置您是如何理解的呀?多重度设置问题解决了吗?
zx4562430 发表于 Post on 2021-12-8 20:26:02
丁越 发表于 2021-12-8 18:38
你看看sob老师的解答http://bbs.keinsci.com/thread-22642-1-1.html,说的很详细了。这里的磁矩不是每个 ...

好的谢谢您
丁越 发表于 Post on 2021-12-8 18:38:38
本帖最后由 丁越 于 2021-12-8 19:47 编辑
zx4562430 发表于 2021-12-8 17:46
您好,请教您下,根据您分享的pdf和multiwfn实际磁矩设置操作,我的理解是磁矩设置应该是针对每个原子的 ...

你看看sob老师的解答http://bbs.keinsci.com/thread-22642-1-1.html,说的很详细了。这里的磁矩不是每个原子给它设一个,这样的话体相和表面相的磁矩必然不同,如此一来太复杂了就没办法算了,我们这里只是设置原子磁矩的初猜值。另外输入文件不都给你了么?你把KIND部分Fe的磁矩设置为3,迭代过程中每个原子的磁矩会慢慢趋于合理。
zx4562430 发表于 Post on 2021-12-8 17:26:33
丁越 发表于 2021-12-8 09:04
首先你看看这个关于表面吸附模型搭建流程的帖子中卡开发发老师的解答(http://bbs.keinsci.com/thread-14 ...

好的。底下两层原子我固定了,我看看您发的关于磁矩的设定,非常谢谢您
丁越 发表于 Post on 2021-12-8 09:04:26
本帖最后由 丁越 于 2021-12-8 09:05 编辑
zx4562430 发表于 2021-12-7 16:51
您好,我这个是基于fe优化好的晶体切的110晶面,准备在真空层放吸附小分子。我看了您的inp文件,是铁金属 ...

首先你看看这个关于表面吸附模型搭建流程的帖子中卡开发发老师的解答(http://bbs.keinsci.com/thread-1447-1-1.html),注意一点是切面优化表面的时候要固定下体相的原子,一般是两层就可以了。另外这里的磁矩不是对单个原子的设定,是整体的设定。看看这个帖子以及下面的ppt有助于你理解。(http://bbs.keinsci.com/thread-22642-1-1.html

hand_on_sprkkr_05.pdf

1016.79 KB, 下载次数 Times of downloads: 148

zx4562430 发表于 Post on 2021-12-7 16:51:43
丁越 发表于 2021-12-7 11:07
首先没必要用超晶胞优化,用原胞设置k点就行。
下面是sob老师在论坛里回复的例子,参考着设置:

您好,我这个是基于fe优化好的晶体切的110晶面,准备在真空层放吸附小分子。我看了您的inp文件,是铁金属的原胞,MAGNETIZATION设为3, MULTIPLICITY设为7,但对于我这个吸附表面模型,共有45个铁原子,如果将全部铁原子MAGNETIZATION设为3, MULTIPLICITY会变为100多,不合理。所以我在想应该将部分铁原子的MAGNETIZATION进行设置,但是目前不知道该设多少,应该设哪些铁原子的磁矩。谢谢

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-15 20:37 , Processed in 0.158444 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list