计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 第一性原理 (First Principle) › 求助，计算不同电场下小体系有机物分子的电子亲和能，电 ...

求助，计算不同电场下小体系有机物分子的电子亲和能，电子亲和能为什么没有稳定

查看数: 728 | 评论数: 7 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

Decade_

发布时间: 2024-6-14 21:16

正文摘要:

使用CP2K计算不同电场强度情况下环氧碎片的垂直电子亲和能以及垂直电离能，使用的公式是 EA=E(N)-E(N+1) IP=E(N-1)-E(N）其中E(N)所采取的是对分子进行单点能计算后的最终体系能量如图1所示所加电场的方 ...

回复 Reply

Decade_ 发表于 Post on 2024-7-4 15:13:08

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

横向对比需要统一

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

sobereva 发表于 Post on 2024-7-4 14:00:07

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

可以考虑wB97X

横向对比需要统一

Decade_ 发表于 Post on 2024-7-3 22:28:38

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

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

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

sobereva 发表于 Post on 2024-6-17 20:52:17

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

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

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

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

Decade_ 发表于 Post on 2024-6-17 10:45:37

本帖最后由 Decade_ 于 2024-7-3 22:29 编辑

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

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

sobereva 发表于 Post on 2024-6-15 00:44:55

当前体系应当用wavelet而非MT。盒子尺寸明显也不合适。而且当前的基组仅适合几何优化用，算能量太糙了。纯泛函算这种体系的电子亲和能误差也很大。

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

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