计算化学公社

标题: 利用cp2k优化a-Fe2O3表面后发生形变的问题 [打印本页]

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

(, 下载次数 Times of downloads: 44) (, 下载次数 Times of downloads: 48) (, 下载次数 Times of downloads: 47) (, 下载次数 Times of downloads: 44)

作者
Author: fanzhichao001 时间: 2024-7-26 16:30
@卡开发发卡卡老师可以帮忙瞅一眼嘛

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

作者
Author: 卡开发发 时间: 2024-7-27 02:25

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

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

作者
Author: fanzhichao001 时间: 2024-7-27 17:07

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

好滴，谢谢老师

作者
Author: fanzhichao001 时间: 2024-7-27 17:18

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

老师这是我输出文件的最后一步
(, 下载次数 Times of downloads: 41)

作者
Author: 卡开发发 时间: 2024-7-27 17:38

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

传个完整的输出看看？

作者
Author: fanzhichao001 时间: 2024-7-27 18:13

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

(, 下载次数 Times of downloads: 5)

作者
Author: 卡开发发 时间: 2024-7-27 18:51

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

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

作者
Author: fanzhichao001 时间: 2024-7-28 14:41

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

好滴，谢谢老师

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