求助铜/环氧界面体系费米能级及电子亲和能电离能计算泛函与基组选择问题

Decade_

最近在进行铜/环氧界面电荷注入势垒的研究，体系包含200个铜原子和三十多个C、H、O、N的环氧分子

目前对于电荷注入势垒计划用图1的公式进行计算，即分别要计算：①体系的真空能级位移、②铜salb的功函数、③环氧分子的垂直电子亲和能和垂直电离能
故而需要进行计算的体系分别为：具有236个分子的界面模型的费米能级和真空能级、铜slab模型的费米能级和真空能级、环氧小分子的E(N)E(N-1)E(N+1)

之前的计算中，在对整个体系进行几何优化时使用的是PBE/DZVP-MOLOPT-SR-GTH，之后依旧用这个进行①和②的计算，都是没有问题的
后面在进行环氧的EAIP计算的时候发现加电场后EAIP出线线性增长，在论坛中提问后采用aug-cc-pVTZ进行计算，发现加Y方向的电场时，EA稳定但出现加正电场时难收敛且IP计算中趋势不收敛的情况

有几个问题想请教一下各位老师
1.能否有一个基组和泛函能够适合进行以上这些全部的计算？如果没有的话，是否可以分部分用不同的泛函和基组进行计算，最后直接用电荷注入势垒公式进行带入计算，就是这个基组和泛函时要统一吗？
2.我使用aug-cc-pVTZ进行EA计算的时候，加负电场收敛情况良好但加正电场收敛很难，请问出线这种情况是什么原因呢？

EAIP.rar (402.62 KB, 下载次数 Times of downloads: 1)

2024-7-1 23:18 上传 Uploaded

点击下载
Click to download
EAIP计算文件