计算化学公社

标题: 求助MoS2表面吸附能的分解，及SAPT(DFT)电离能的计算问题 [打印本页]

作者
Author: yufeng 时间: 2024-10-8 15:33
标题: 求助MoS2表面吸附能的分解，及SAPT(DFT)电离能的计算问题
各位老师好，

1、我想算一个水分子在MoS2表面吸附能的分解，大约90~100个原子，用PSI4里的SAPT0计算很难收敛，而且计算速度非常慢。学校没有高斯软件，用不了sob老师说的sobEDA。请问有什么解决的办法吗？
2、另外想问下PSI4里的SAPT(DFT)电离能如何计算呢？我看官网教程写着用中性单体减去少一个电子的阳离子单体的能量就可以，这个能量用VASP软件GGA-PBE来算可以吗？官网给出的例子中H2O的电离能是0.1307 Hartree，但是查找NIST Chemistry Webbook中H2O的电离能是0.46 Hartree，这是为什么呢？

作者
Author: sobereva 时间: 2024-10-8 16:41
置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义和含糊性，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “SAPT(DFT)电离能计算” 改了，以后务必注意

作者
Author: sobereva 时间: 2024-10-8 16:42
最便宜的SAPT0都算不动，SAPT(DFT)算这个也根本算不动
电离能按照常规方式，少一个电子状态能量和原状态能量求差就完了

SAPT不支持周期性体系。倘若你用簇模型来算，又毫无理由用VASP，没买Gaussian可以用ORCA或PSI4自己

作者
Author: yufeng 时间: 2024-10-9 10:25

sobereva 发表于 2024-10-8 16:42
最便宜的SAPT0都算不动，SAPT(DFT)算这个也根本算不动
电离能按照常规方式，少一个电子状态能量和原状态能 ...

谢谢sob老师。那为什么PSI4给出水的电离能是0.1307 Hartee，而查NIST表是0.46 Hartree，相差这么大呢？
这个是PSI4说明链接https://psicode.org/psi4manual/1.4.0/sapt#a-first-example

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