计算化学公社

标题: 求助关于电子亲合能输出文件能量读取问题 [打印本页]

作者
Author: Decade_ 时间: 2024-7-31 10:33
标题: 求助关于电子亲合能输出文件能量读取问题
最近在使用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吗？
(, 下载次数 Times of downloads: 19) (, 下载次数 Times of downloads: 19)
(, 下载次数 Times of downloads: 3) (, 下载次数 Times of downloads: 1)

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

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

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