“第10届量子化学波函数分析与Multiwfn程序培训班将于5月4-8日于北京举办,这是一次性完整、系统学习波函数分析的各种理论知识和全面掌握强大的Multiwfn波函数分析程序使用的最不可错过的机会!请点击此链接查看详情和报名方式,欢迎参加!

“第18届北京科音分子动力学与GROMACS培训班” 将于5月23-26日于北京举办。这是一次性全面、系统学习分子动力学模拟知识和最流行的分子动力学程序GROMACS的关键机会!报名正在进行中,请点击此链接查看详情,欢迎参加!

计算化学公社

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

[CP2K] 请教AIMD的系综选择和如何能减少SCF收敛的步数?

[复制链接 Copy URL]

53

帖子

0

威望

598

eV
积分
651

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 白菜 于 2026-3-23 08:21 编辑

老师们好~~想请教老师们关于AIMD的一些设置问题:
我在做一个铵盐离子晶体A的AIMD,想得到它的热分解机理(文献报道的后期涉及自由基)。下面是目前的计算考量:
1. 建立了某个晶面的slab模型,没有用体相是因为要和A上吸附别的分子的热分解做对比。
2. 取优化后的结构作为初始结构,全程用NVT做模拟,因为文献里报道的不多,所以想用1000K、1500K、2000K、2500K、3000K分别跑2ps观察一下,再决定用什么温度。
3. 在群里请教过sob老师后,"RKS并开较高温度的smearing能勉强给出正确的解离曲线,参考GFN1-xTB原文里的讨论。",所以用RKS结合smearing(电子温度分别设为和NVT设置的体系温度一样)
遇到的问题是:
1. 全程用NVT可以吗?我看见有一篇文献是先用NVT模拟“点燃”,然后再切换到NVE。(我认为用全程用NVT就行了,但也不确定。)
2. 电子温度分别设为和1000K、1500K、2000K、2500K、3000K可以吗?(我设置的原因:GFN1-xTB原文里是对比了电子温度300K和4000K的解离曲线,平衡区和键长完全不变;结合搜到cp2k论坛上电子温度设成和离子温度一样,理由是希望体系达到热力学平衡。)
3. 现在几个温度下的模拟的每一step的SCF收敛步数都是十几(还没有成断键),导致整体AIMD花费时间很久,请问老师们这个SCF收敛步数还有降低的空间吗?之前做了同A的slab上+特定分子的AIMD,298K,0.5fs,开OT,SCF收敛只要三步左右。

谢谢老师们~~~



6万

帖子

99

威望

6万

eV
积分
126147

管理员

公社社长

2#
发表于 Post on 2026-3-20 04:12:39 | 只看该作者 Only view this author
1 通常可以
2 电子温度先设为预期的热分解温度跑。比如2500-3000 K一般就可以
3 高温下原子平均运动速度快,每一步位移比常温大得多,因此平均每一步SCF收敛所需迭代次数肯定比常温多得多。减小步长虽然可以减小迭代到收敛的平均圈数数,但并不能提升模拟性能。

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
白菜 + 5 谢谢伟大的sob老师~

查看全部评分 View all ratings

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

53

帖子

0

威望

598

eV
积分
651

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2026-3-20 08:55:16 | 只看该作者 Only view this author
谢谢伟大的sob老师~

本版积分规则 Credits rule

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

GMT+8, 2026-4-13 23:53 , Processed in 0.156887 second(s), 21 queries , Gzip On.

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