计算化学公社

标题: 请问AIMD计算氢气高温解离的过程体系爆掉什么原因？ [打印本页]

作者
Author: 软哨豆米 时间: 2024-12-3 22:38
标题: 请问AIMD计算氢气高温解离的过程体系爆掉什么原因？
各位老师晚上好，AIMD初学者求解，最近在用AIMD模拟氢气分子高温弛豫的过程，目的是想探索氢气分子解离的温度，我算了三组数据，分别在2000K，2500K，3000K温度下高温弛豫，2000K与3000K的体系计算正常，体系没有爆掉(散掉)，2500K的体系温度在跑了一定步数后，温度上升至上百万K，体系整个散掉(下面附2500K数据的inp文件)。
请问各位老师两个问题：
1.三组数据的inp文件仅控温字段的数字不同，其他完全相同，为什么会出现不一样的结果？
2.想算氢气分子解离温度，用AIMD高温弛豫的方式是否正确？
万分感谢

作者
Author: sobereva 时间: 2024-12-5 02:40
GFN1-xTB算这个根本不靠谱
本来氢的质量就非常小，动力学还想牵扯到断键，1fs步长根本不能用。0.1、0.2 fs还差不多
用闭壳层形式计算还想有正确的共价键解离曲线必须开smearing，参考GFN1-xTB的原文里的相关信息

作者
Author: 软哨豆米 时间: 2024-12-5 10:15

sobereva 发表于 2024-12-5 02:40
GFN1-xTB算这个根本不靠谱
本来氢的质量就非常小，动力学还想牵扯到断键，1fs步长根本不能用。0.1、0.2 fs ...

谢谢社长，请问社长该体系用什么泛函得到的结果比较可靠

作者
Author: sobereva 时间: 2024-12-5 19:38

软哨豆米发表于 2024-12-5 10:15
谢谢社长，请问社长该体系用什么泛函得到的结果比较可靠

用一个算H2的H-H键解离能相对较准的泛函
可以用Gaussian计算，对不同泛函进行测试

作者
Author: 软哨豆米 时间: 2024-12-5 20:51

sobereva 发表于 2024-12-5 19:38
用一个算H2的H-H键解离能相对较准的泛函
可以用Gaussian计算，对不同泛函进行测试

谢谢社长

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