计算化学公社

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

[Multiwfn使用咨询] 周期性分子晶体结构能否利用Multiwfn来计算弱相互作用

[复制链接 Copy URL]

33

帖子

0

威望

837

eV
积分
870

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
各位老师好,最近我在做周期性分子结构的分子间相互作用(氢键、π堆积及范德华力)对其结构性能的影响。我查阅了相关文献发现大部分第一性原理研究都是从电子密度差(EDD)及Hirshfeld表面来分析弱相互作用的类型及强度。其中有一部分关于独立小分子体系的研究使用了Multiwfn软件进行RDG分析来描述弱相互作用,解释的很到位。
在这想请教一下各位老师,带有周期性的分子晶体(已用MS进行结构优化),能否用该软件进行RDG分析?我的想法是把cif文件转成mol2,然后再转换成gjf和wfn文件,请问这样可行吗?以及在结构已经被优化好的前提下(能够得到相应的坐标信息),周期性是否影响对弱相互作用力的分析解释?
学生刚入此门,对相关基础知识的了解仍在努力深入中,还请各位老师不吝赐教!多谢了!

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

2#
发表于 Post on 2018-3-28 22:52:36 | 只看该作者 Only view this author
可以

一种做法是先把晶体当中抠出一个团簇,用量化程序算个单点得到波函数文件,这样就可以用Multiwfn里所有弱相互作用分析方法,见
Multiwfn支持的弱相互作用的分析方法概览
http://sobereva.com/252
诸如Multiwfn原文里对尿素晶体中的弱相互作用做的RDG分析就是这样做的。

还一种做法是直接把晶体结构保存成pdb或xyz或mol等含有结构信息的格式,这样在Multiwfn里可以用那些不依赖于波函数而只依赖于结构的弱相互作用分析方法,包括Hirshfeld surface、promolecular近似下的RDG,以及promolecular近似下的IGM。这些方法在上面的博文都有提到。



后记:以上说法已经过时,Multiwfn早已直接支持对周期性波函数做IGMH、IRI、NCI等分析,例如参考
使用Multiwfn做IGMH分析非常清晰直观地展现化学体系中的相互作用
http://sobereva.com/621http://bbs.keinsci.com/thread-28147-1-1.html
使用Multiwfn结合CP2K通过NCI和IGM方法图形化考察固体和表面的弱相互作用(注:NCI和IGM如今都强烈不推荐使用了)
http://sobereva.com/588http://bbs.keinsci.com/thread-21742-1-1.html

北京科音自然科学研究中心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

33

帖子

0

威望

837

eV
积分
870

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2018-3-29 10:47:30 | 只看该作者 Only view this author
多谢sob老师,我回去再细看一下这些博文和例子!

24

帖子

0

威望

87

eV
积分
111

Level 2 能力者

4#
发表于 Post on 2025-3-15 11:14:48 | 只看该作者 Only view this author
sobereva 发表于 2018-3-28 22:52
可以

一种做法是先把晶体当中抠出一个团簇,用量化程序算个单点得到波函数文件,这样就可以用Multiwfn里 ...

不需要结构优化吗?单点能直接得到的fchk能直接用于分析IGMH吗

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

5#
发表于 Post on 2025-3-17 11:55:34 | 只看该作者 Only view this author
zhangxinxin 发表于 2025-3-15 11:14
不需要结构优化吗?单点能直接得到的fchk能直接用于分析IGMH吗

实验测定分子结构的方法以及将实验结构用于量子化学计算需要注意的问题
http://sobereva.com/569
北京科音自然科学研究中心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

72

帖子

0

威望

387

eV
积分
459

Level 3 能力者

6#
发表于 Post on 2025-3-30 17:22:35 | 只看该作者 Only view this author
sobereva 发表于 2018-3-28 22:52
可以

一种做法是先把晶体当中抠出一个团簇,用量化程序算个单点得到波函数文件,这样就可以用Multiwfn里 ...

请教下,我从晶体结构里面扣出单胞结构保存为xyz文件,然后准备进行vasp或cp2k的周期性计算,盒子大小基于什么设置?我发现盒子太小原子跑散开了。

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

7#
发表于 Post on 2025-3-30 19:40:25 | 只看该作者 Only view this author
wilely 发表于 2025-3-30 17:22
请教下,我从晶体结构里面扣出单胞结构保存为xyz文件,然后准备进行vasp或cp2k的周期性计算,盒子大小基 ...

我不知道你算的是什么晶体、为什么要抠
通常来说,有个cif文件,直接用Multiwfn载入、按下文产生CP2K输入文件就完了,此时晶胞信息也是cif里的晶胞信息
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587http://bbs.keinsci.com/thread-21668-1-1.html

除非你要算孤立体系,否则盒子的设置没有任意性,必须和体系结构相符、满足周期性。
北京科音自然科学研究中心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

72

帖子

0

威望

387

eV
积分
459

Level 3 能力者

8#
发表于 Post on 2025-3-31 11:28:23 | 只看该作者 Only view this author
sobereva 发表于 2025-3-30 19:40
我不知道你算的是什么晶体、为什么要抠
通常来说,有个cif文件,直接用Multiwfn载入、按下文产生CP2K输 ...

我是两个有机分子组成的复合物晶体,我想把其中一个分子用类似的分子替换掉,再做周期性计算,盒子的大小怎么设置?

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

9#
发表于 Post on 2025-3-31 14:12:09 | 只看该作者 Only view this author
wilely 发表于 2025-3-31 11:28
我是两个有机分子组成的复合物晶体,我想把其中一个分子用类似的分子替换掉,再做周期性计算,盒子的大小 ...

直接替换没意义。分子的一点改变都可能造成实际晶体结构的很大变化。需要做晶体结构预测得到实际晶体结构。
北京科音自然科学研究中心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

72

帖子

0

威望

387

eV
积分
459

Level 3 能力者

10#
发表于 Post on 2025-4-1 18:18:43 | 只看该作者 Only view this author
sobereva 发表于 2025-3-31 14:12
直接替换没意义。分子的一点改变都可能造成实际晶体结构的很大变化。需要做晶体结构预测得到实际晶体结构 ...

Multiwfn有没有把cif文件导出,只含单个晶胞的盒子的形式保存?我用vasp算

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

11#
发表于 Post on 2025-4-1 22:28:18 | 只看该作者 Only view this author
wilely 发表于 2025-4-1 18:18
Multiwfn有没有把cif文件导出,只含单个晶胞的盒子的形式保存?我用vasp算

病句。语义不明
如果指把超胞约化成原胞,没这功能
北京科音自然科学研究中心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

72

帖子

0

威望

387

eV
积分
459

Level 3 能力者

12#
发表于 Post on 2025-4-4 15:55:25 | 只看该作者 Only view this author
sobereva 发表于 2025-4-1 22:28
病句。语义不明
如果指把超胞约化成原胞,没这功能

超胞约化成原胞可以用mercury直接得到,关键是后面怎么设置盒子大小用vasp来算,这个跟体系有关吧。

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

13#
发表于 Post on 2025-4-4 16:30:31 | 只看该作者 Only view this author
wilely 发表于 2025-4-4 15:55
超胞约化成原胞可以用mercury直接得到,关键是后面怎么设置盒子大小用vasp来算,这个跟体系有关吧。

约化完了自然就有现成的盒子信息,直接用于CP2K等计算就完了
北京科音自然科学研究中心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

72

帖子

0

威望

387

eV
积分
459

Level 3 能力者

14#
发表于 Post on 2025-4-6 21:35:19 | 只看该作者 Only view this author
sobereva 发表于 2025-4-4 16:30
约化完了自然就有现成的盒子信息,直接用于CP2K等计算就完了

好的,thanks

本版积分规则 Credits rule

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

GMT+8, 2025-8-12 21:01 , Processed in 0.176745 second(s), 20 queries , Gzip On.

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