求助关于电子亲合能输出文件能量读取问题

Decade_ · 发表于 Post on 2024-7-31 10:33:17

最近在使用cp2k计算小分子的垂直电子亲合能和电离能，使用的公式是
EA=E(N)-E(N+1)
IP=E(N-1)-E(N）

想请教一下输出文件中得这么一大串能量，应该选取哪个呢？
看了老师的http://bbs.keinsci.com/thread-13450-1-1.html这篇帖子中，请问是要在Hartree energy和Total energy中选择Hartree energy吗？

1.inp (6.17 KB, 下载次数 Times of downloads: 3)

1.out (147.5 KB, 下载次数 Times of downloads: 1)

zjxitcc · 发表于 Post on 2024-7-31 19:57:17

1.out文件中第1543行Total energy。提醒一下，这类问题最好采用量化程序进行孤立体系计算，各有各的擅长之处和应用场景。先不说你的盒子可能不够大，纯泛函PBE精度是不够的，在量子化学程序中有一箩筐更准确的方法。

sobereva · 发表于 Post on 2024-8-1 06:33:36

用纯泛函算电子亲和能属于开玩笑，离域性误差极大。对于DFT，需要用长程HF成分较高的泛函，如wB97XD，才能算得不错
而且DZVP-MOLOPT-SR-GTH仅仅适合优化、振动分析、MD等对基组要求不高但又耗时的情况用，拿来算能量相关问题太烂了。

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

[CP2K] 求助关于电子亲合能输出文件能量读取问题