计算化学公社

标题: 求助，计算不同电场下小体系有机物分子的电子亲和能，电子亲和能为什么没有稳定 [打印本页]

作者
Author: Decade_ 时间: 2024-6-14 21:16
标题: 求助，计算不同电场下小体系有机物分子的电子亲和能，电子亲和能为什么没有稳定
使用CP2K计算不同电场强度情况下环氧碎片的垂直电子亲和能以及垂直电离能，使用的公式是
EA=E(N)-E(N+1)
IP=E(N-1)-E(N）

其中E(N)所采取的是对分子进行单点能计算后的最终体系能量如图1所示

所加电场的方式如图2所示
带一个正电荷时CHARGE为1，MULTIPLICITY为2；

电中性时CHARGE为0，MULTIPLICITY为1；

所使用的泛函为PBE-D3(BJ)/DZVP-MOLOPT-SR-GTH

参照文献计算的相同的分子的这两个值，如图3.1和3.2。文献中所说电子亲和能和电离能随电场改变的变化不大，但我计算出来的两个值是呈上升趋势的（图4）

想请问一下，
1.现在的基组泛函和CP2K能不能算我的这个小体系
2.现在的体系加电场结果不合理，能不能在这个条件下算

计算出的在不加电压下，加-1和+1电压下的输入输出文件如下

(, 下载次数 Times of downloads: 14) (, 下载次数 Times of downloads: 14) (, 下载次数 Times of downloads: 14) (, 下载次数 Times of downloads: 15) (, 下载次数 Times of downloads: 16)
(, 下载次数 Times of downloads: 2)

作者
Author: sobereva 时间: 2024-6-15 00:44
当前体系应当用wavelet而非MT。盒子尺寸明显也不合适。而且当前的基组仅适合几何优化用，算能量太糙了。纯泛函算这种体系的电子亲和能误差也很大。

这种体系用Gaussian等量子化学程序远比CP2K适合得多。强烈建议改用量子化学程序算这个问题。如果有什么特殊原因必须用CP2K不可，也建议至少先用量子化学程序把文献里的结论重现出来，之后再试图用CP2K重现。

作者
Author: Decade_ 时间: 2024-6-17 10:45
本帖最后由 Decade_ 于 2024-7-3 22:29 编辑

sobereva 发表于 2024-6-15 00:44
当前体系应当用wavelet而非MT。盒子尺寸明显也不合适。而且当前的基组仅适合几何优化用，算能量太糙了。纯 ...

老师能给我推荐一个适用于cp2k计算的基组吗？我看了培训的ppt也在站内搜索了一下，只搜到适用于gaussian的基组是M062X/def2TZVPP，但是用multiwfn生成cp2k的输入文件时没有找到对应的基组

作者
Author: sobereva 时间: 2024-6-17 20:52

Decade_ 发表于 2024-6-17 10:45
老师能给我推荐一个适用于cp2k计算的基组吗？我看了培训的ppt也在站内搜索了一下，只搜到适用于gaussian ...

你上传图片方式不对，别人看不到

def2-TZVPP不带弥散函数，算电子亲和能绝对不能用

电子亲和能一般可以用aug-cc-pVTZ

作者
Author: Decade_ 时间: 2024-6-24 11:42

sobereva 发表于 2024-6-17 20:52
你上传图片方式不对，别人看不到

def2-TZVPP不带弥散函数，算电子亲和能绝对不能用

老师，请问
1.我要对这个体系（Cu晶面/环氧碎片）计算加电场下，计算金属的功函数、真空能级位移、环氧的电离能和电子亲和能，是否能有一个合适的泛函在cp2k进行这些计算？
2.我使用了aug-cc-pVTZ在cp2k中计算电离能和电子亲合能，发现随着电场变化，EA和IP还是线性上升，但环氧分子的E(N)随电场变化很小，几乎没有，主要是E(N-1)和E(N+1)在跟着电场变化，请问是我的输入文件设置的有问题吗？因为始终无法得到EA和IP随电场变化基本不变的结果
3.wavelet和MT分别是计算什么体系的，想问一下他们的区别是什么？
(, 下载次数 Times of downloads: 16) (, 下载次数 Times of downloads: 15) (, 下载次数 Times of downloads: 2)

作者
Author: Decade_ 时间: 2024-7-3 22:28

sobereva 发表于 2024-6-17 20:52
你上传图片方式不对，别人看不到

def2-TZVPP不带弥散函数，算电子亲和能绝对不能用

谢谢老师回复，我试着用aug-cc-pVTZ计算了EA，发现在加y方向的电场时，EA终于稳定！
同时我在论坛中找到计算有机分子的IP推荐使用M06-2X
但是出现了一个问题，由于还要算整个铜/环氧界面体系的费米能级，发现这个泛函满足不了铜元素的计算，想请教一下是否有一个泛函能满足铜/环氧分子界面和单独环氧分子的EAIP计算的泛函和机基组呢，因为这些值是要计算的电荷注入势垒的公式中几部分，泛函和基组应该是要统一的吗？

作者
Author: sobereva 时间: 2024-7-4 14:00

Decade_ 发表于 2024-7-3 22:28
谢谢老师回复，我试着用aug-cc-pVTZ计算了EA，发现在加y方向的电场时，EA终于稳定！
同时我在论坛中找到 ...

可以考虑wB97X

横向对比需要统一

作者
Author: Decade_ 时间: 2024-7-4 15:13

sobereva 发表于 2024-7-4 14:00
可以考虑wB97X

横向对比需要统一

好的好的，谢谢老师！我试着跑下

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