计算化学公社

标题: 请教AIMD的系综选择和如何能减少SCF收敛的步数？ [打印本页]

作者
Author: 白菜 时间: 2026-3-20 01:10
标题: 请教AIMD的系综选择和如何能减少SCF收敛的步数？
本帖最后由白菜于 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收敛只要三步左右。

谢谢老师们~~~

作者
Author: sobereva 时间: 2026-3-20 04:12
1 通常可以
2 电子温度先设为预期的热分解温度跑。比如2500-3000 K一般就可以
3 高温下原子平均运动速度快，每一步位移比常温大得多，因此平均每一步SCF收敛所需迭代次数肯定比常温多得多。减小步长虽然可以减小迭代到收敛的平均圈数数，但并不能提升模拟性能。

作者
Author: 白菜 时间: 2026-3-20 08:55
谢谢伟大的sob老师~

欢迎光临计算化学公社 (http://bbs.keinsci.com/)