计算化学公社

 找回密码 Forget password
 注册 Register
Views: 8175|回复 Reply: 11
打印 Print 上一主题 Last thread 下一主题 Next thread

[CP2K] 使用cp2k对MOF几何结构优化时,发现计算出来的结构和初始结构差距太大

[复制链接 Copy URL]

37

帖子

0

威望

578

eV
积分
615

Level 4 (黑子)

有点小迷糊

各位社友们大家好,在使用cp2k对SIFSIX-2-CU进行几何结构优化时,得到的结构和初始结构差距非常大,而且跑出了晶胞外,不管是是用OT还是对角化算法,都不行,输入文件我是用Multiwfn生成的,然后修改了一下参数,采用的结构优化算法是BFGS,使用的cp2k版本是8.2。请问,会是什么原因导致的呢?初始结构不合理还是参数设置的问题呢?

202108051010132769..png (45.75 KB, 下载次数 Times of downloads: 25)

优化结果最后一帧

优化结果最后一帧

202108051010042670..png (21.57 KB, 下载次数 Times of downloads: 33)

初始结构

初始结构

OT.inp

6.78 KB, 下载次数 Times of downloads: 39

OT算法输入文件

diag.inp

8.09 KB, 下载次数 Times of downloads: 12

对角化算法输入文件

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
John_Tao + 3 谢谢

查看全部评分 View all ratings

生如逆旅,一苇以航。

691

帖子

2

威望

4250

eV
积分
4981

Level 6 (一方通行)

2#
发表于 Post on 2021-8-5 11:27:59 | 只看该作者 Only view this author
你这个体系选取的cell 多大? 使用的K点是多少?

691

帖子

2

威望

4250

eV
积分
4981

Level 6 (一方通行)

3#
发表于 Post on 2021-8-5 11:33:17 | 只看该作者 Only view this author
估计是晶胞取小了 另外可以考虑用PBE+D3试试

87

帖子

0

威望

4006

eV
积分
4093

Level 6 (一方通行)

4#
发表于 Post on 2021-8-5 12:07:26 | 只看该作者 Only view this author
结构也不大,为啥不用VASP优化,还有用OT算法的时候晶胞c轴可以乘以2,使用1*1*2超包计算

37

帖子

0

威望

578

eV
积分
615

Level 4 (黑子)

有点小迷糊

5#
 楼主 Author| 发表于 Post on 2021-8-5 12:38:06 | 只看该作者 Only view this author
jiangning198511 发表于 2021-8-5 11:27
你这个体系选取的cell 多大? 使用的K点是多少?

对角化计算时没有设置K点,  只用了Gamma点进行计算, 晶胞大小是13.63160000        13.63160000       7.96800000
生如逆旅,一苇以航。

37

帖子

0

威望

578

eV
积分
615

Level 4 (黑子)

有点小迷糊

6#
 楼主 Author| 发表于 Post on 2021-8-5 12:39:48 | 只看该作者 Only view this author
jiangning198511 发表于 2021-8-5 11:33
估计是晶胞取小了 另外可以考虑用PBE+D3试试

好的,非常感谢您。
生如逆旅,一苇以航。

37

帖子

0

威望

578

eV
积分
615

Level 4 (黑子)

有点小迷糊

7#
 楼主 Author| 发表于 Post on 2021-8-5 15:27:21 | 只看该作者 Only view this author
waitingseven 发表于 2021-8-5 12:07
结构也不大,为啥不用VASP优化,还有用OT算法的时候晶胞c轴可以乘以2,使用1*1*2超包计算

好的谢谢您,我修改了输入文件提交计算了,之所以用cp2k是因为我只会cp2k,我师姐会VASP,所以我们老师让我学cp2k
生如逆旅,一苇以航。

32

帖子

0

威望

457

eV
积分
489

Level 3 能力者

8#
发表于 Post on 2021-9-14 21:42:18 | 只看该作者 Only view this author
你好 请问你的问题解决了吗 可以相互交流一下吗 谢谢
QQ:842655799

496

帖子

11

威望

4275

eV
积分
4991

Level 6 (一方通行)

9#
发表于 Post on 2021-9-14 22:29:33 | 只看该作者 Only view this author
你自己别改那些参数,Multiwfn产生的输入文件基本都是合理的,我看你改了一大通还改的不合理,重新生成一下再算。
自由发挥,野蛮生长

5

帖子

0

威望

147

eV
积分
152

Level 3 能力者

10#
发表于 Post on 2021-9-15 09:02:25 | 只看该作者 Only view this author
请问楼主解决这个问题了吗
磨练心性,强健体魄,增长学识

37

帖子

0

威望

578

eV
积分
615

Level 4 (黑子)

有点小迷糊

11#
 楼主 Author| 发表于 Post on 2021-10-16 14:47:00 | 只看该作者 Only view this author
842655799 发表于 2021-9-14 21:42
你好 请问你的问题解决了吗 可以相互交流一下吗 谢谢
QQ:842655799

问题解决啦
生如逆旅,一苇以航。

60

帖子

0

威望

1668

eV
积分
1728

Level 5 (御坂)

12#
发表于 Post on 2021-10-16 14:55:29 | 只看该作者 Only view this author
假设结构是解单晶得到的,那是不是在优化的时候可以把非氢原子固定呢?
I AM an Amphibious Warrior.

本版积分规则 Credits rule

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

GMT+8, 2026-2-19 00:13 , Processed in 0.184685 second(s), 24 queries , Gzip On.

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