计算化学公社

标题: 使用VASP计算天然气水合物的客体分子相互作用能是否合理 [打印本页]

作者
Author: D-opamine 时间: 2023-6-7 16:41
标题: 使用VASP计算天然气水合物的客体分子相互作用能是否合理
本帖最后由 D-opamine 于 2023-6-7 16:45 编辑

各位老师同学好，我是做天然气水合物的，想问一下，使用VASP计算天然气水合物中的客体分子与笼子之间的相互作用能，（使用VASP）是否合理？下图是lammps跑的结构，我想通过取这个超胞中的一个水笼子，计算与客体分子的相互作用能。

起因是，因为看到了一个帖子 http://bbs.keinsci.com/thread-23358-1-1.html 这个大佬说：孤立的、不是太小的有机分子，适合用高斯、ORCA等软件来算。VASP基本上只适合算固体和表面。
需要用VASP算分子只有一种情况，就是你研究一个分子在表面上的吸附，因为计算吸附能之类的，需要减去孤立分子的能量，又因为理论级别统一的要求，孤立分子必须和表面用同一个软件、同样的计算方法，这时用VASP算分子是合理的。或者比如说你算cohesive energy，需要减去原子的能量，这个原子用VASP算，是合理的。其他情况基本不要用VASP算孤立的分子、原子。

sob老师也说，VASP算分子体系太烂了

所以对我这样使用VASP计算相互作用能的合理性产生了怀疑，是否合理，或者有更好的软件可以代替VASP计算？

作者
Author: Eudaimonia 时间: 2023-6-7 17:09
实际上vasp更适合的场景是周期性体系中的性质，对于孤立体系的分子不合适，并不是说只适用于固体和表面

其他软件也可以考虑QE，Cp2k等

作者
Author: sobereva 时间: 2023-6-8 06:59
VASP不是不能算，但明显有更好的选择，用Gaussian或ORCA能用更高级别理论方法结合全电子基组计算，而且还不用考虑结果随真空层厚度收敛性问题。
如果你当周期性体系而不取团簇来算，VASP也很难算得动动，而算大体系周期性计算极快的CP2K则没问题（哪怕上千原子做纯泛函+D3的计算在主流双路服务器上也没问题）

作者
Author: D-opamine 时间: 2023-6-8 11:36

sobereva 发表于 2023-6-8 06:59
VASP不是不能算，但明显有更好的选择，用Gaussian或ORCA能用更高级别理论方法结合全电子基组计算，而且还不 ...

明白了，感谢sob老师

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