fundamental gap的值为负数的原因是什么？

Xun_L_iii · 发表于 Post on 2020-7-8 18:47:03

本帖最后由 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是负数，请问是哪一步出了问题？

sobereva · 发表于 Post on 2020-7-9 07:25:56

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

Xun_L_iii · 发表于 Post on 2020-7-9 15:10:00

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).

Xun_L_iii · 发表于 Post on 2020-7-9 15:13:03

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

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

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

[综合交流] fundamental gap的值为负数的原因是什么？

浏览过的版块