计算化学公社

 找回密码 Forget password
 注册 Register

使用cp2k进行Ni(OH)2结构优化过程中SCF/几何优化无法收敛的问题

查看数: 11665 | 评论数: 20 | 收藏 Add to favorites 16
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2021-12-25 13:52

正文摘要:

本帖最后由 Eudaimonia 于 2021-12-28 21:06 编辑 各位老师好,我使用cp2k优化Ni(OH)2表面,参考卡开发发老师的回答进行计算   [综合讨论] 新手刚开始做固体表面吸附,应该怎么样入手呢1.先对下载得到 ...

回复 Reply

双木aaa 发表于 Post on 2024-6-27 21:52:02
请问一下,楼主解决了吗
Janson_Shi 发表于 Post on 2023-6-26 13:24:28
Eudaimonia 发表于 2021-12-28 20:28
手滑恢复了错误的数据。。。。重新再写一次感谢卡开发发老师的进一步指点,反铁磁的情况我之前确实没有注意 ...

请问下我什么开不了对称性K点
量化小手 发表于 Post on 2022-10-31 17:45:44
楼主解决磁性表面体系下的收敛问题了,最近我也在搞Ni(OH)2表面,试了很多方法,同发现很难收敛
卡开发发 发表于 Post on 2022-9-18 20:56:02
zmjsce 发表于 2022-9-18 20:36
谢谢卡开发发老师的回复!
那对于磁性体系的晶胞参数优化,可能就不再适合使用CP2K了。
目前CP2K的Quic ...

cp2k确实多k点方面功能支持不算完善。
zmjsce 发表于 Post on 2022-9-18 20:36:48
卡开发发 发表于 2022-9-18 18:04
如果你采用的是纯泛函(LDA/GGA),一般来说会因为这些泛函的自作用问题而导致低估带隙、磁矩等,如果要 ...

谢谢卡开发发老师的回复!
那对于磁性体系的晶胞参数优化,可能就不再适合使用CP2K了。
目前CP2K的Quickstep模块还不能在多K点优化的情况下使用DFT+U
卡开发发 发表于 Post on 2022-9-18 18:04:39
zmjsce 发表于 2022-9-18 16:27
卡开发发老师您好,请问使用CP2K优化Ni(OH)2时,在已经考虑磁性情况下,是否还需要考虑DFT+U呢?

如果你采用的是纯泛函(LDA/GGA),一般来说会因为这些泛函的自作用问题而导致低估带隙、磁矩等,如果要讨论电子结构原则上还是需要的。
zmjsce 发表于 Post on 2022-9-18 16:27:26
卡开发发 发表于 2021-12-28 21:08
针对末尾的问题:1、对称性的问题实在不行你可以用一些程序约束磁矩差异或者临时将spin down的原子换成别的 ...

卡开发发老师您好,请问使用CP2K优化Ni(OH)2时,在已经考虑磁性情况下,是否还需要考虑DFT+U呢?
卡开发发 发表于 Post on 2022-1-1 11:14:33
Eudaimonia 发表于 2022-1-1 11:04
谢谢卡开发发老师的指点,老师元旦快乐~
前几天去上海做实验了,没有及时回复和对问题进行测试,不好意 ...

1、磁性约束指的是针对不等价磁矩设置的原子应当被识别为是不同类型的原子,目前能做到的例如MS或Pymatgen,在这个基础上找对称性看看获得的原胞如何。

欢迎讨论。
Eudaimonia 发表于 Post on 2022-1-1 11:04:54
卡开发发 发表于 2021-12-28 21:08
针对末尾的问题:1、对称性的问题实在不行你可以用一些程序约束磁矩差异或者临时将spin down的原子换成别的 ...

谢谢卡开发发老师的指点,老师元旦快乐~
前几天去上海做实验了,没有及时回复和对问题进行测试,不好意思
1.对论坛中进行了检索,没有找到对磁性进行约束的讨论,请问一下老师能推荐一下相关程序吗?
4.嗯嗯
5.后续会对此处重新进行测试

谢谢老师的提示,论坛中有很多关于K点测试的讨论,10楼中的测试结果是单点的结果(不过取用的模型是已经优化过的331超胞(就是经过上百次Opt step重新收敛后的),可能和单胞扩胞直接测试有一些差别,后续可能重新测试并更新再后面楼层,最近搬家+期末事情比较杂乱
卡开发发 发表于 Post on 2021-12-28 21:08:41
本帖最后由 卡开发发 于 2021-12-28 21:25 编辑

针对末尾的问题:1、对称性的问题实在不行你可以用一些程序约束磁矩差异或者临时将spin down的原子换成别的元素,重新找一下对称性试试,我不确定是不是程序没办法正确处理AFM的对称性。
2、不明,不知道程序内部怎么设计的。
3、那个Warning其实可以不用管,对窄带系统其实可以不要管自旋污染的问题,上面帖子提到过原因。
4、按道理两个Ni磁矩相反,总的多重度应该为1。
5、先看看纠正上述问题是否还有这个情况,不行再根据输出情况进行诊断。
另外,收敛性测试最好是在单点上做,然后看能量和力的收敛情况,这样能有效避免构型优化震荡。如果这种情况下还出现震荡,要么初始结构不合理,要么可以适当调整构型优化中算法的参数(比如步长控制或Hessian更新之类的算法)。
sobereva 发表于 Post on 2021-12-28 21:06:47
以后发帖时别在贴子里出现fan qiang,过于敏感,给你改了
卡开发发 发表于 Post on 2021-12-26 22:07:25
本帖最后由 卡开发发 于 2021-12-26 22:26 编辑
Eudaimonia 发表于 2021-12-26 21:57
(3).哦哦,不好意思犯了这样的错误
我尝试了一下,对原始的NI(OH)2进行了晶胞优化,但是发现晶胞 ...

(3)反铁磁(AFM)的情况Ni有磁矩为+3的,也有为-3的,z方向的层和层之间Ni的磁矩应该是相反的。前面专门提了反铁磁的事情,你可能没注意。
后续表面也是相同的问题。
表面体系应该5 5 1就行,5 5 2是针对块体的,你可以好好测试下单点,看看在能量受力等差异不大的情况可以使用稍小一些的k点,另外&KPOINTS这里,SYMMETRY F完全没必要,如果晶体本身有对称性可以用这个对称性来使用更少的(不可约)k点进行计算,有些k的贡献是一样的;如果能确定对称性,晶胞优化也尽可能使用对称性。这体系尽量猜测原子磁矩并且自动优化磁矩。
Eudaimonia 发表于 Post on 2021-12-26 21:57:37
卡开发发 发表于 2021-12-26 14:19
0、(3)有点搞事情啊。你要看清楚研究的问题,该贴中后续是要进行一个MD的过程,如果后续是NPT,晶格在 ...

(3).哦哦,不好意思犯了这样的错误
我尝试了一下,对原始的NI(OH)2进行了晶胞优化,但是发现晶胞参数发生了非常夸张的变化,从3.16x3.16x9.35变化为7.56x7.7x4.95(还未优化完),但大体结构还在,感觉也不合理,请问老师可能是哪里可能有问题?对Ni猜测磁矩为3,增加了D3校正和Smear,其他参数使用Multiwfn默认参数
1.3.4.5 感谢老师的指点,麻烦老师了

其他:
重新对Ni猜测磁矩为3,增加了Smear进行了单胞及001切面的优化,还是可以正常完成。441晶面仍然在初始SCF计算时震荡,221、551、552的k点都进行了尝试,均震荡,相关文件更新在1楼,请问老师还有什么可以进行的尝试?
Note:k点设置为552的文件不小心被覆盖,进行了约60步的SCF计算,震荡,实际电脑不太算的动这个等级的k点
目前打算更改一下超胞大小再进行尝试

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

GMT+8, 2024-11-24 20:56 , Processed in 0.186088 second(s), 25 queries , Gzip On.

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