计算化学公社

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

[CP2K] 有关计算MOF/COF与分子之间的结合能疑问

[复制链接 Copy URL]

100

帖子

0

威望

545

eV
积分
645

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
各位老师,有关计算MOF/COF与分子之间的结合能我有两个疑问
1.我用CP2K2022.2版本计算MOF/COF吸附分子的结合能,均选择PBE、DZVP-MOLOPT-SR-GTH、DFTD3(BJ),对角化算法,开启smear,CUTOFF  500、REL_CUTOFF  60,请问是否合适?后续想做一下AIM、IGMH分析,以及请问各位老师是否还有其他分析可以做呢?2.目前我先对MOF/COF进行结构优化,把分子摆放合适位置再次结构优化,再用BSSE计算结合能。请问把分子摆放进来后做结构优化和MD哪个合适呢?
附上目前计算的结果,COF结构优化,MOF结构优化,MOF吸附分子结构优化,MOF吸附分子结合能,暂时还没有做频率分析。



MOF&COF.rar

173.82 KB, 下载次数 Times of downloads: 22

6万

帖子

99

威望

5万

eV
积分
120146

管理员

公社社长

2#
发表于 Post on 2024-7-26 04:17:09 | 只看该作者 Only view this author
基组太小,而且此时还不考虑counterpoise算物理吸附能属于开玩笑(做几何优化可以),看http://bbs.keinsci.com/thread-47399-1-1.html里我的回复
还有很多可做的分析,比如
使用Multiwfn图形化展现原子对色散能的贡献以及色散密度
http://sobereva.com/705http://bbs.keinsci.com/thread-44723-1-1.html
谈谈范德华势以及在Multiwfn中的计算、分析和绘制
http://sobereva.com/551http://bbs.keinsci.com/thread-17271-1-1.html
使用Multiwfn做Hirshfeld surface分析直观展现分子晶体和复合物中的相互作用
http://sobereva.com/701http://bbs.keinsci.com/thread-43178-1-1.html
以及小分子和固体表面静电势、密度差,等等
参考下文(不是所有都支持周期性体系)
Multiwfn支持的弱相互作用的分析方法概览
http://sobereva.com/252

如果有化学吸附出现的情况,还可以结合键级、NAdO、CDA分析,看
使用Multiwfn对周期性体系做键级分析和NAdO分析考察成键特征
http://sobereva.com/719http://bbs.keinsci.com/thread-47176-1-1.html
使用Multiwfn结合CP2K对周期性体系做电荷分解分析(CDA)
http://sobereva.com/716http://bbs.keinsci.com/thread-46878-1-1.html


不要乱用BSSE这个词,仔细看下文相关部分
计算化学中的一些常见不良写法和用词
http://sobereva.com/298http://bbs.keinsci.com/thread-1358-1-1.html
即便用counterpoise考虑BSSE校正,DZVP-MOLOPT-SR-GTH也明显偏小,没理由不用更大的TZ档次的MOLOPT基组

结合方式可能性较少或凭感觉就能摆出最理想的吸附构型的时候直接摆一个结构做优化就完了,拿不准的话跑MD,每隔一定帧数取一个结构做优化,取能量最低的用于算高精度结合能

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

100

帖子

0

威望

545

eV
积分
645

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2024-7-26 14:32:04 | 只看该作者 Only view this author
本帖最后由 红茶泡枸杞 于 2024-7-27 16:13 编辑
sobereva 发表于 2024-7-26 04:17
基组太小,而且此时还不考虑counterpoise算物理吸附能属于开玩笑(做几何优化可以),看http://bbs.keinsci ...

谢谢社长的解答,您之前提到“Multiwfn创建CP2K输入文件的界面里选择任务的时候直接就能选counterpoise”,是指选任务11BSSE吗?我现在可以在原来小基组优化结束的基础上换成TZVP-MOLOPT-GTH基组再次优化吗?把分子摆放到合适位置是否需要固定MOF/COF呢?

6万

帖子

99

威望

5万

eV
积分
120146

管理员

公社社长

4#
发表于 Post on 2024-7-26 23:41:29 | 只看该作者 Only view this author
红茶泡枸杞 发表于 2024-7-26 14:32
谢谢社长的解答,您之前提到“Multiwfn创建CP2K输入文件的界面里选择任务的时候直接就能选counterpoise” ...


可。通常来说用DZVP-MOLOPT-SR-GTH只优化一次就够了,没必要再切换成更大基组
没必要也没理由固定

那叫TZVP-MOLOPT-GTH
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

100

帖子

0

威望

545

eV
积分
645

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2024-7-27 16:17:28 | 只看该作者 Only view this author
sobereva 发表于 2024-7-26 23:41

可。通常来说用DZVP-MOLOPT-SR-GTH只优化一次就够了,没必要再切换成更大基组
没必要也没理由固定

谢谢社长,我明白了,接下来的计算我用TZVP-MOLOPT-GTH基组,用counterpoise考虑BSSE,上面写错的基组修改过来了。

100

帖子

0

威望

545

eV
积分
645

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2024-8-1 20:50:18 | 只看该作者 Only view this author
sobereva 发表于 2024-7-26 23:41

可。通常来说用DZVP-MOLOPT-SR-GTH只优化一次就够了,没必要再切换成更大基组
没必要也没理由固定

社长您好,我把分子摆放在COF中进行结构优化,优化结束后发现COF有些位置散架了,这怎么解决呢?附上输出文件,由于中途服务器关机了,我续算忘了保存之前的out文件了,故输出文件只有续算部分的。

TPA-S.out

39.29 KB, 下载次数 Times of downloads: 6

6万

帖子

99

威望

5万

eV
积分
120146

管理员

公社社长

7#
发表于 Post on 2024-8-2 03:52:14 | 只看该作者 Only view this author
红茶泡枸杞 发表于 2024-8-1 20:50
社长您好,我把分子摆放在COF中进行结构优化,优化结束后发现COF有些位置散架了,这怎么解决呢?附上输出 ...

光上传out文件毫无用处,连用了什么关键词、结构怎么变化的都不直接体现
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

100

帖子

0

威望

545

eV
积分
645

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2024-8-2 14:58:03 | 只看该作者 Only view this author
sobereva 发表于 2024-8-2 03:52
光上传out文件毫无用处,连用了什么关键词、结构怎么变化的都不直接体现

社长您好,我重新上传了输入文件,molden、xyz等文件过大,附在网盘链接中。
通过网盘分享的文件:COF链接: https://pan.baidu.com/s/1yRmuaEo4GHEKoXB5y1sAqw?pwd=6666 提取码: 6666

TPA-S.inp

31.56 KB, 下载次数 Times of downloads: 10

6万

帖子

99

威望

5万

eV
积分
120146

管理员

公社社长

9#
发表于 Post on 2024-8-4 07:43:18 | 只看该作者 Only view this author
红茶泡枸杞 发表于 2024-8-2 14:58
社长您好,我重新上传了输入文件,molden、xyz等文件过大,附在网盘链接中。
通过网盘分享的文件:COF链 ...

当前的模型严重不合理,Z方向盒子尺寸太小,导致小分子都和自己作用上了。必须先在z方向扩胞
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

100

帖子

0

威望

545

eV
积分
645

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2024-9-9 17:55:35 | 只看该作者 Only view this author
sobereva 发表于 2024-8-4 07:43
当前的模型严重不合理,Z方向盒子尺寸太小,导致小分子都和自己作用上了。必须先在z方向扩胞

社长您好,我用1*1*3超胞重新进行了计算,结果如下,优化前后变化不是很大,这是正常的吗?是否需要把z方向尺寸调大一点呢?out文件和molden文件附在百度网盘链接中了
通过百度网盘分享的文件:R输出
链接:https://pan.baidu.com/s/1chrXFcFk710EUW3wSs7cRA?pwd=z78i
提取码:z78i

R.inp

34.12 KB, 下载次数 Times of downloads: 4

R-GEO.cif

27.08 KB, 下载次数 Times of downloads: 7

R.cif

52.94 KB, 下载次数 Times of downloads: 3

6万

帖子

99

威望

5万

eV
积分
120146

管理员

公社社长

11#
发表于 Post on 2024-9-10 14:35:15 | 只看该作者 Only view this author
红茶泡枸杞 发表于 2024-9-9 17:55
社长您好,我用1*1*3超胞重新进行了计算,结果如下,优化前后变化不是很大,这是正常的吗?是否需要把z方 ...

如果你是指相对于cif里的结构变化不大,明显是好事,说明当前级别和设置对当前体系描述得合理
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

100

帖子

0

威望

545

eV
积分
645

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2024-9-10 21:10:13 | 只看该作者 Only view this author
sobereva 发表于 2024-9-10 14:35
如果你是指相对于cif里的结构变化不大,明显是好事,说明当前级别和设置对当前体系描述得合理

好的,谢谢社长解答,我可以去做分析了

100

帖子

0

威望

545

eV
积分
645

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2024-9-11 09:39:25 | 只看该作者 Only view this author
sobereva 发表于 2024-9-10 14:35
如果你是指相对于cif里的结构变化不大,明显是好事,说明当前级别和设置对当前体系描述得合理

社长,不好意思,我再跟您确认一下,我把分子放在COF中进行结构优化,COF略有一点点变形,分子的位置优化前后变化不大,这样可以继续做后续弱相互作用分析吗?

6万

帖子

99

威望

5万

eV
积分
120146

管理员

公社社长

14#
发表于 Post on 2024-9-11 11:45:19 | 只看该作者 Only view this author
红茶泡枸杞 发表于 2024-9-11 09:39
社长,不好意思,我再跟您确认一下,我把分子放在COF中进行结构优化,COF略有一点点变形,分子的位置优化 ...

COF本身有的地方就很柔,弱相互作用造成适度的变形很正常
但也有可能模型本身不合理导致的异常(比如垂直于层的方向晶胞尺寸太小而又没充分扩胞等),没具体文件没法说
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

100

帖子

0

威望

545

eV
积分
645

Level 4 (黑子)

15#
 楼主 Author| 发表于 Post on 2024-9-12 15:59:29 | 只看该作者 Only view this author
sobereva 发表于 2024-9-11 11:45
COF本身有的地方就很柔,弱相互作用造成适度的变形很正常
但也有可能模型本身不合理导致的异常(比如垂 ...

好的谢谢社长,之前做的1*1*2超胞优化后严重变形,分子内发生了相互作用,现在做1*1*3,Z方向尺寸为11.1747埃,有轻微变形,可能就是您说的弱相互作用造成的,附上优化前后的cif文件。

R.cif

52.94 KB, 下载次数 Times of downloads: 3

R-GEO.cif

27.08 KB, 下载次数 Times of downloads: 5

本版积分规则 Credits rule

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

GMT+8, 2025-8-15 22:42 , Processed in 0.234047 second(s), 23 queries , Gzip On.

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