计算化学公社

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

[新手求助] 关于周期性体系与小分子物质的结构优化

[复制链接 Copy URL]

17

帖子

0

威望

197

eV
积分
214

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
之前使用Gaussian对二聚体或者多聚体进行结构优化时,基本都是参照了社长的“使用molclus程序做团簇构型搜索和分子构象搜索”这个帖子,去得到一个比较合理的结构去做结构优化,确实非常好用!但是目前想通过cp2k去计算周期性体系比如Mxene与小分子物质的相互作用时,在初始的建模当中,没办法建立一个比较合理的结构,导致后续使用cp2k做结构优化时,消耗的时间巨久(目前是想探究MXene与纤维素单体即葡萄糖的相互作用)。除了参照有关的文献的做法,另外还想问问大佬们在这方面是如何做的呢?或者molclus程序是否也适用于这种体系做团簇构型搜索呢?如果适用,其使用方法是否有区别呢?

6万

帖子

99

威望

5万

eV
积分
124671

管理员

公社社长

2#
发表于 Post on 2025-11-20 00:09:52 | 只看该作者 Only view this author
当前不牵扯团簇问题,“这种体系做团簇构型搜索”描述不当

VMD里自己把小分子挪到MXene表面合适位置作为初始结构计算。如果你担心有不同吸附位点、不同吸附朝向,手动摆考虑不周全,可以CP2K里尝试用GFN1-xTB或者更便宜的GFN0-xTB跑动力学,每隔一定帧数把结构优化一次,对比能量,进行初筛。再把能量最低的几个结构进一步用DFT做优化,取能量最低的。

molclus未来可能支持结合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

17

帖子

0

威望

197

eV
积分
214

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2025-11-20 09:32:36 | 只看该作者 Only view this author
好的社长,感谢回复

17

帖子

0

威望

197

eV
积分
214

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2025-11-20 20:26:17 | 只看该作者 Only view this author
sobereva 发表于 2025-11-20 00:09
当前不牵扯团簇问题,“这种体系做团簇构型搜索”描述不当

VMD里自己把小分子挪到MXene表面合适位置作为 ...

sob老师,我按照您的方法采用CP2K用GFN1-xTB跑MD后,得到的一个xx-pos-1.xyz的文件,这里面记录了不同结构的原子坐标,以及对应的能量。就您的这句话“每隔一定帧数把结构优化一次,对比能量,进行初筛”是指在这个文件当中进行初筛是吗,后续我要对能量较低的结构用DFT做优化的时,应该如何做呢?(这个我之前在Gaussian中是直接复制能量较低的结构的原子坐标去做结构优化的,但是在CP2K中我不太清楚这个如何做)
另外,我对已经结构优化成功的体系,要去做更高精度的单点能计算时,我是使用Multiwfn载入相应的restart文件,之后在生成相应的CP2K输入文件吗?望老师解答

6万

帖子

99

威望

5万

eV
积分
124671

管理员

公社社长

5#
发表于 Post on 2025-11-21 03:49:44 | 只看该作者 Only view this author
顶顶顶顶顶 发表于 2025-11-20 20:26
sob老师,我按照您的方法采用CP2K用GFN1-xTB跑MD后,得到的一个xx-pos-1.xyz的文件,这里面记录了不 ...

跟molclus的思路一样,每隔一定帧数保存一帧结构,自己写脚本用Multiwfn批量产生几何优化的输入文件、批量让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

17

帖子

0

威望

197

eV
积分
214

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2025-11-21 09:34:53 | 只看该作者 Only view this author
sobereva 发表于 2025-11-21 03:49
跟molclus的思路一样,每隔一定帧数保存一帧结构,自己写脚本用Multiwfn批量产生几何优化的输入文件、批 ...

好的老师

本版积分规则 Credits rule

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

GMT+8, 2026-1-24 07:31 , Processed in 0.232506 second(s), 20 queries , Gzip On.

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