利用cp2k优化a-Fe2O3表面后发生形变的问题

fanzhichao001 · 发表于 Post on 2024-7-26 16:29:19

在利用cp2k优化a-Fe2O3表面时，为了方便设置边界条件，首先用Atomsk对a-Fe2O3进行正交化处理，参考以下链接中的例子进行了设置，以下是本人的输入文件，优化后感觉会发生形变，并且与vasp结构优化后的结构不太一样，这样的形变和差异正常嘛？
参考例子：events:2015_cecam_tutorial:ling_hybrids.pdf (cp2k.org)
以下为本人cp2k的输入文件与vasp使用的INCAR及优化前后的结构。

fe2o3.inp (2.67 KB, 下载次数 Times of downloads: 11)

fanzhichao001 · 发表于 Post on 2024-7-26 16:30:16

@卡开发发卡卡老师可以帮忙瞅一眼嘛

sobereva · 发表于 Post on 2024-7-26 22:43:09

如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义和含糊性，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “cp2k计算a-Fe2O3体系” 改了，以后务必注意！下次将删帖+扣分处理。

卡开发发 · 发表于 Post on 2024-7-27 02:25:44

fanzhichao001 发表于 2024-7-26 16:30
@卡开发发卡卡老师可以帮忙瞅一眼嘛

我不确定你当中计算过程输出是怎样的情况，比如SCF是否收敛等，另外cp2k这边采用的是blyp的泛函，这和vasp应该也不一样。

fanzhichao001 · 发表于 Post on 2024-7-27 17:07:53

sobereva 发表于 2024-7-26 22:43
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此 ...

好滴，谢谢老师

fanzhichao001 · 发表于 Post on 2024-7-27 17:18:42

卡开发发发表于 2024-7-27 02:25
我不确定你当中计算过程输出是怎样的情况，比如SCF是否收敛等，另外cp2k这边采用的是blyp的泛函，这和vas ...

老师这是我输出文件的最后一步

卡开发发 · 发表于 Post on 2024-7-27 17:38:01

fanzhichao001 发表于 2024-7-27 17:18
老师这是我输出文件的最后一步

传个完整的输出看看？

fanzhichao001 · 发表于 Post on 2024-7-27 18:13:44

卡开发发发表于 2024-7-27 17:38
传个完整的输出看看？

cp2k.zip (115.28 KB, 下载次数 Times of downloads: 4)

卡开发发 · 发表于 Post on 2024-7-27 18:51:14

卡开发发发表于 2024-7-27 17:38
传个完整的输出看看？

磁矩是怎样设置的这里也看不出，因为看不到你的inp的xyz，另外存在的两个问题是：
1、CP2K在DFT+U的方案与VASP有差别，我我猜这不是导致算出来差异这么大的原因。
2、SCF前几步没收敛，这有可能是导致后续结构不合理的原因。
或许你可以在p1x1表面上先优化一下再去进行正交化的处理，会节约一些时间。

fanzhichao001 · 发表于 Post on 2024-7-28 14:41:57

卡开发发发表于 2024-7-27 18:51
磁矩是怎样设置的这里也看不出，因为看不到你的inp的xyz，另外存在的两个问题是：
1、CP2K在DFT+U的方案 ...

好滴，谢谢老师

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

[CP2K] 利用cp2k优化a-Fe2O3表面后发生形变的问题