计算化学公社

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

[CP2K] CP2K用CI-NEB搜索过渡态时能量问题

[复制链接 Copy URL]

128

帖子

0

威望

611

eV
积分
739

Level 4 (黑子)

各位,我CP2K中用CI-NEB搜索过渡态时遇到两个问题,请教一下大家,非常感谢!
1、由于我提前优化了反应物和产物,搜索过渡态时就关闭了优化始态和末态,可是我发现搜索过渡态过程中始末态的能量和优化完的反应物和产物能量不一样,请问这是怎么回事?


2、搜索过渡态时发现如上图所示趋势不是简单的先增后减,于是我提取了能量低的2点和5点重新优化后搜索过渡态,这次还是出现差不多的趋势,请问这该怎么解决?



p.xyz

21.81 KB, 下载次数 Times of downloads: 8

r.xyz

21.81 KB, 下载次数 Times of downloads: 8

ts.inp

28.72 KB, 下载次数 Times of downloads: 33

40

帖子

0

威望

1175

eV
积分
1215

Level 4 (黑子)

2#
发表于 Post on 2023-12-25 16:30:08 | 只看该作者 Only view this author
本帖最后由 jim 于 2023-12-25 16:53 编辑

LZ,请问你的这个大分子体系用PBE泛函来计算,你的服务器能扛得住?

128

帖子

0

威望

611

eV
积分
739

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2023-12-25 17:21:07 | 只看该作者 Only view this author
jim 发表于 2023-12-25 16:30
LZ,请问你的这个大分子体系用PBE泛函来计算,你的服务器能扛得住?

体系确实有点大,不过没办法,反应物就那么大了,服务器是双路80物理核心的还可以

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

4#
发表于 Post on 2023-12-26 04:19:03 | 只看该作者 Only view this author
jim 发表于 2023-12-25 16:30
LZ,请问你的这个大分子体系用PBE泛函来计算,你的服务器能扛得住?

这种大小的体系用速度极快的CP2K做纯泛函计算毫无压力
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

5#
发表于 Post on 2023-12-26 06:16:33 | 只看该作者 Only view this author
确保优化极小点和跑CI-NEB的和能量有关的设置严格一致,如CUTOFF、色散校正、盒子参数等,输入文件对照着仔细检查
对于含Mg的体系结合GTH赝势基组,当前CUTOFF明显偏小
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

128

帖子

0

威望

611

eV
积分
739

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2023-12-26 11:36:04 | 只看该作者 Only view this author
sobereva 发表于 2023-12-26 06:16
确保优化极小点和跑CI-NEB的和能量有关的设置严格一致,如CUTOFF、色散校正、盒子参数等,输入文件对照着仔 ...

原来设置要一致啊,谢谢sob老师指点

40

帖子

0

威望

1175

eV
积分
1215

Level 4 (黑子)

7#
发表于 Post on 2023-12-26 17:19:30 | 只看该作者 Only view this author
sobereva 发表于 2023-12-26 06:16
确保优化极小点和跑CI-NEB的和能量有关的设置严格一致,如CUTOFF、色散校正、盒子参数等,输入文件对照着仔 ...

sob老师,请问这种大分子与催化材料结合的时候是否要用molclus来优化下它们结合的部位,还是直接做结构优化?

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

8#
发表于 Post on 2023-12-27 01:13:53 | 只看该作者 Only view this author
jim 发表于 2023-12-26 17:19
sob老师,请问这种大分子与催化材料结合的时候是否要用molclus来优化下它们结合的部位,还是直接做结构优 ...

拿单一初猜结构优化可能会优化到亚稳结构
有条件的话最好跑动力学,每隔一定时间取一个结构当初始结构,分别优化,取能量最低的作为最稳定的吸附结构
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

40

帖子

0

威望

1175

eV
积分
1215

Level 4 (黑子)

9#
发表于 Post on 2023-12-27 17:01:21 | 只看该作者 Only view this author
sobereva 发表于 2023-12-27 01:13
拿单一初猜结构优化可能会优化到亚稳结构
有条件的话最好跑动力学,每隔一定时间取一个结构当初始结构, ...

请问sob老师,以LZ的研究对象为例,请问是否可以按照如下步骤搜索大分子与催化材料最佳结合部位:(参照您的帖子“使用Molclus结合xtb做的动力学模拟对瑞德西韦(Remdesivir)做构象搜索”
1、在Gaussian或MS等软件建模,将大分子与催化材料摆放差不多的位置,然后生成xyz文件,命名为:A.xyz;
2、构建md.inp文件,参照您帖子中的“4第一步:用xtb在GFN0-xTB下做动力学模拟”;
3、参照您帖子中的“5 第二步:用Molclus调用xtb在GFN0-xTB下做批量优化”;
4、参照您帖子中的“6 第三步:用Molclus调用xtb在GFN2-xTB结合隐式水模型下做批量优化”;
5、还需要参照您帖子中的“7 第四步:用Molclus调用Gaussian和ORCA得到水环境下的高精度自由能”来进一步更高级别的计算吗?如果需要,但是这么大体系调用Gaussian/ORCA,成功的几率应该不大吧?
辛苦sob老师!谢谢!

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

10#
发表于 Post on 2023-12-28 10:05:57 | 只看该作者 Only view this author
jim 发表于 2023-12-27 17:01
请问sob老师,以LZ的研究对象为例,请问是否可以按照如下步骤搜索大分子与催化材料最佳结合部位:(参照 ...

当前是对周期性模型做周期性计算,跟Gaussian、ORCA、xtb根本没直接关系
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

40

帖子

0

威望

1175

eV
积分
1215

Level 4 (黑子)

11#
发表于 Post on 2023-12-28 11:25:01 | 只看该作者 Only view this author
sobereva 发表于 2023-12-28 10:05
当前是对周期性模型做周期性计算,跟Gaussian、ORCA、xtb根本没直接关系

这个帖子我知道是搜索过渡态,我是针对您回答的8#有关亚稳结构来说的,因为您的帖子“使用Molclus结合xtb做的动力学模拟对瑞德西韦(Remdesivir)做构象搜索”中没涉及到催化材料吸附大分子相关的例子,所以借楼问问您。我和LZ做得有些相似,也是催化材料催化降解大分子,所以特别关心,谢谢您。

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

12#
发表于 Post on 2023-12-29 04:38:14 | 只看该作者 Only view this author
jim 发表于 2023-12-28 11:25
这个帖子我知道是搜索过渡态,我是针对您回答的8#有关亚稳结构来说的,因为您的帖子“使用Molclus结合xtb ...

催化材料是周期性体系,和瑞德西韦那篇文章要用的程序明显不一样,瑞德西韦又不是周期性体系
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

40

帖子

0

威望

1175

eV
积分
1215

Level 4 (黑子)

13#
发表于 Post on 2023-12-30 10:17:24 | 只看该作者 Only view this author
sobereva 发表于 2023-12-29 04:38
催化材料是周期性体系,和瑞德西韦那篇文章要用的程序明显不一样,瑞德西韦又不是周期性体系

明白了,谢谢sob老师。我在论坛中查了下,好像没见molclus应用于周期性材料的例子或者帖子,请问molclus是否可以应用于这种周期性材料与大分子结合方式的搜索?

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

14#
发表于 Post on 2023-12-31 10:28:31 | 只看该作者 Only view this author
jim 发表于 2023-12-30 10:17
明白了,谢谢sob老师。我在论坛中查了下,好像没见molclus应用于周期性材料的例子或者帖子,请问molclus ...

不正常周期性的情况
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

40

帖子

0

威望

1175

eV
积分
1215

Level 4 (黑子)

15#
发表于 Post on 2024-1-2 10:23:16 | 只看该作者 Only view this author
sobereva 发表于 2023-12-31 10:28
不正常周期性的情况

sob老师新年好!非常冒昧地问一下,我搜索了“不正常周期性、Molclus”,好像没见相关的内容,请问是否能再稍微详细地指点我们一下。我现在想将周期性材料与大分子结合方式全局搜索一下,卡这儿了动不了了,我的猜想是:全局搜索后再进一步结构优化,然后再用CI-NEB来搜索过渡态,看能否突破一二,谢谢!

本版积分规则 Credits rule

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

GMT+8, 2024-11-23 23:32 , Processed in 0.181702 second(s), 24 queries , Gzip On.

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