计算化学公社

 找回密码 Forget password
 注册 Register

CP2K XTB与PBE方法优化MOF结构不合理的问题

查看数: 5644 | 评论数: 10 | 收藏 Add to favorites 3
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2021-11-3 10:46

正文摘要:

本帖最后由 函数与激情 于 2021-11-3 10:57 编辑 大家好,我想在CP2K中使用xtb方法优化MOF分子ZIF8,但是在优化过程中发现结构完全垮掉,想询问大家有没有办法处理一下~后续也尝试了PBE-D3优化,结构依然会垮掉, ...

回复 Reply

jy1165254033 发表于 Post on 2024-8-22 14:05:40
函数与激情 发表于 2021-11-3 19:47
抱歉,问题可能找到了,配体是带单个电子的,不是电中性,整体不该是+8价,这应该是结构坍塌的原因。

请问电荷和自旋应该设置成多少,我也需要算ZIF8
renzhogn424 发表于 Post on 2021-11-4 12:04:37
函数与激情 发表于 2021-11-4 09:44
没有冻结原子,不是一个完整的晶胞,结构你可以看我的附件cc.xyz

好的,谢谢。我自己也常常拿个片段出来跑。跑着跑着 散了
函数与激情 发表于 Post on 2021-11-4 09:44:41
renzhogn424 发表于 2021-11-4 08:47
请问一下,用xtb6.3.2优化非周期的时候,你冻结边缘原子么??还是抠出一个晶胞直接优化?

没有冻结原子,不是一个完整的晶胞,结构你可以看我的附件cc.xyz
renzhogn424 发表于 Post on 2021-11-4 08:47:31
请问一下,用xtb6.3.2优化非周期的时候,你冻结边缘原子么??还是抠出一个晶胞直接优化?
喵星大佬 发表于 Post on 2021-11-4 00:18:05
函数与激情 发表于 2021-11-3 19:47
抱歉,问题可能找到了,配体是带单个电子的,不是电中性,整体不该是+8价,这应该是结构坍塌的原因。

我就想不应该有问题,去年用xtb算过ZIF8的团簇是OK的,这种问题多半是电荷/自旋没搞对
函数与激情 发表于 Post on 2021-11-3 19:47:38
抱歉,问题可能找到了,配体是带单个电子的,不是电中性,整体不该是+8价,这应该是结构坍塌的原因。
函数与激情 发表于 Post on 2021-11-3 19:31:03
丁越 发表于 2021-11-3 14:26
你到底是要优化晶胞还是要优化结构?看你的意思应该是属于优化晶胞,改成CELL_OPT。
另外优化晶胞把这个搞 ...

谢谢回复,我是优化结构,晶胞也尝试优化过,依然会散掉。
函数与激情 发表于 Post on 2021-11-3 19:30:17
wth1219 发表于 2021-11-3 12:12
根据个人使用经验提供两点建议:
1. 周期性边界条件下的计算中盒子的a、b两边似乎不足以消除镜像之间的相 ...

谢谢回复,a,b两边的盒子大小应该是没有问题的,和实验上的cif文件中的尺寸几乎相当,而且化合价0价的情况下优化出来的情况是和气相中的情况也是较为类似的,+8价就会突然坍塌。
丁越 发表于 Post on 2021-11-3 14:26:08
本帖最后由 丁越 于 2021-11-3 14:34 编辑

你到底是要优化晶胞还是要优化结构?看你的意思应该是属于优化晶胞,改成CELL_OPT。
另外优化晶胞把这个搞成类似于Slab模型干什么,ABC得是晶胞边长。MOF的话a~15埃就够了,不够的话下面的过程在多个扩胞(-9,3)。

将zif-8的cif载入Multiwfn
cp2k
-1
4
3
2
4
5
0
就ok了

202111031423329092..png (1.06 KB, 下载次数 Times of downloads: 24)

202111031423329092..png
wth1219 发表于 Post on 2021-11-3 12:12:35
根据个人使用经验提供两点建议:
1. 周期性边界条件下的计算中盒子的a、b两边似乎不足以消除镜像之间的相互作用;
2. CP2K也可以开展非周期性边界条件下的计算,如果使用POISSON_SOLVER MT的话要求盒子足够大。

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

GMT+8, 2024-11-24 04:10 , Processed in 0.187649 second(s), 26 queries , Gzip On.

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