计算化学公社

标题: 请教大家，算气体能量，丙二烯比丙炔还稳定是怎么回事呢？ [打印本页]

作者
Author: pupu1616 时间: 2023-7-10 19:25
标题: 请教大家，算气体能量，丙二烯比丙炔还稳定是怎么回事呢？
用的GGA-PBE，30*30*30的超胞，单Γ点

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

作者
Author: pupu1616 时间: 2023-7-10 20:03

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

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

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

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

作者
Author: pupu1616 时间: 2023-7-11 16:06

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

ok，谢谢谢谢

作者
Author: pupu1616 时间: 2023-7-11 16:06

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

明白了，谢谢sob老师

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