请问AIMD计算氢气高温解离的过程体系爆掉什么原因？

查看数: 638 | 评论数: 4 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

发布时间: 2024-12-3 22:38

各位老师晚上好，AIMD初学者求解，最近在用AIMD模拟氢气分子高温弛豫的过程，目的是想探索氢气分子解离的温度，我算了三组数据，分别在2000K，2500K，3000K温度下高温弛豫，2000K与3000K的体系计算正常，体系没有爆 ...

回复 Reply

软哨豆米 发表于 Post on 2024-12-5 20:51:53

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

谢谢社长

sobereva 发表于 Post on 2024-12-5 19:38:13

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

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

软哨豆米 发表于 Post on 2024-12-5 10:15:13

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

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

sobereva 发表于 Post on 2024-12-5 02:40:40

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register