请教大家，算气体能量，丙二烯比丙炔还稳定是怎么回事呢？

pupu1616 · 发表于 Post on 2023-7-10 19:25:05

用的GGA-PBE，30*30*30的超胞，单Γ点

pal · 发表于 Post on 2023-7-10 19:53:09

算的是单点能，还是优化呢，如果是单点能，怎么构建的两个分子的结构，如果是优化，优化之后还是丙二烯么

pupu1616 · 发表于 Post on 2023-7-10 20:03:19

pal 发表于 2023-7-10 19:53
算的是单点能，还是优化呢，如果是单点能，怎么构建的两个分子的结构，如果是优化，优化之后还是丙二烯么

是优化，分子结构从库里扒的，优化后还是丙二烯

agent99 · 发表于 Post on 2023-7-11 00:29:45

我查了一下NIST数据库，这两个分子的能量顺序和方法有关，就连热力学组合方法给出的顺序都不一样。对于你用的PBE泛函，在TZVP, cc-PVTZ这些基组下确实都是丙二烯能量更低。只能说反应活性不能和能量高低划等号（虽然大多数情况可以）

sobereva · 发表于 Post on 2023-7-11 12:32:08

我不用VASP不对输入文件进行评价。只是强调一点，对于能量相近的同分异构分子，PBE的精度经常不足以正确判断能量最低结构，如果必须用纯泛函的话，建议用其中算热力学性质几乎最好的B97M-rV或者MN15L，如果程序都不支持，用r2SCAN也明显好于PBE。如果能用杂化泛函，首选wB97M-V，至少用更流行、支持得很广泛的M06-2X。以如今的视角来看，PBE算有机体系能量问题的精度实在太烂了。

pupu1616 · 发表于 Post on 2023-7-11 16:06:39

agent99 发表于 2023-7-11 00:29
我查了一下NIST数据库，这两个分子的能量顺序和方法有关，就连热力学组合方法给出的顺序都不一样。对于你用 ...

ok，谢谢谢谢

pupu1616 · 发表于 Post on 2023-7-11 16:06:50

sobereva 发表于 2023-7-11 12:32
我不用VASP不对输入文件进行评价。只是强调一点，对于能量相近的同分异构分子，PBE的精度经常不足以正确判 ...

明白了，谢谢sob老师

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

[VASP] 请教大家，算气体能量，丙二烯比丙炔还稳定是怎么回事呢？