计算化学公社

标题: fundamental gap的值为负数的原因是什么？ [打印本页]

作者
Author: Xun_L_iii 时间: 2020-7-8 18:47
标题: fundamental gap的值为负数的原因是什么？
本帖最后由 Xun_L_iii 于 2020-7-8 18:47 编辑

各位老师，我想要计算Co2+与取代酞菁形成的配合物的fundamental gap，下面是我的计算步骤（各步骤输出文件已上传）。
1. 使用PBE0/genecp进行结构优化（中心原子及配位原子用def2-TZVP，其他原子用6-31G*，中心原子加SDD赝势），检查为4 YES；
2. freq检查没有虚频；
3. 使用PBE0/def2-TZVP pseudo=read进行单点计算(金属原子加SDD赝势)；
4. 比较G确定是4重态。
5. N-1态自旋多重度是1、3、5，N+1态自旋多重度是1、3。计算各自旋多重度结构的单点能，比较数据。（加关键词“stable”检查波函数稳定性，部分结构波函数不稳定而加关键词“stable=opt”）
6. 取最稳定结构的单点能数据与N态单点能数据计算得出VIP和VEA，从而得出fundamental gap。
数据如附件excel中所示，但fundamental gap是负数，请问是哪一步出了问题？

作者
Author: sobereva 时间: 2020-7-9 07:25
我目前没时间仔细看
要上传就把gjf文件一起上传，否则别人不好检查你赝势使用到底对不对
没必要计算G，就拿电子能量对比就够了
算VEA必须带弥散函数，否则根本不准。当前中性、阳离子、阴离子体系都应当带弥散函数算

作者
Author: Xun_L_iii 时间: 2020-7-9 15:10

sobereva 发表于 2020-7-9 07:25
我目前没时间仔细看
要上传就把gjf文件一起上传，否则别人不好检查你赝势使用到底对不对
没必要计算G，就 ...

老师，gjf文件被覆盖掉了，这是我使用的赝势：（从http://www.tc.uni-koeln.de/PP/clickpse.en.html上下载）
Co 0
ECP10MDF 3 10
F-Komponente
1
2 1.000000 0.000000
S-F
2
2 23.660000 283.960566
2 10.610000 47.156846
P-F
2
2 25.040000 182.212236
2 10.440000 35.233352
D-F
2
2 29.540000 -26.475333
2 10.180000 -1.825787
! References:
! [7] M. Dolg, U. Wedig, H. Stoll, H. Preuss, J. Chem. Phys. 86, 866 (1987).

作者
Author: Xun_L_iii 时间: 2020-7-9 15:13

sobereva 发表于 2020-7-9 07:25
我目前没时间仔细看
要上传就把gjf文件一起上传，否则别人不好检查你赝势使用到底对不对
没必要计算G，就 ...

谢谢老师，我用ma-TZVP试一试

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