计算化学公社

标题: 请教下算一个过渡金属复合物体系，轨道能量有正值 [打印本页]

作者
Author: 413 时间: 2023-10-9 14:50
标题: 请教下算一个过渡金属复合物体系，轨道能量有正值
印象里轨道能量都是负的
最近算出来的一组数据出现了正值，请问正常吗？
1205 2.0000    -0.093922       -2.5558
......
1214 0.0000    -0.025817       -0.7025
1215 0.0000    -0.010044       -0.2733
1216 0.0000    -0.008665       -0.2358
1217 0.0000    -0.008244       -0.2243
1218 0.0000    -0.000000       -0.0000
1219 0.0000    -0.000000       -0.0000
1220 0.0000    -0.000000       -0.0000
1221 0.0000    -0.000000       -0.0000
1222 0.0000    -0.000000       -0.0000
1223 0.0000    -0.000000       -0.0000
1224 0.0000    -0.000000       -0.0000
1225 0.0000    -0.000000       -0.0000
1226 0.0000    0.000000       0.0000
1227 0.0000    0.000000       0.0000
1228 0.0000    0.000000       0.0000
1229 0.0000    0.000000       0.0000
1230 0.0000    0.000000       0.0000
1231 0.0000    0.000000       0.0000
1232 0.0000    0.000000       0.0000
1233 0.0000    0.000000       0.0000
1234 0.0000    0.000000       0.0000
1235 0.0000    0.000000       0.0000
1236 0.0000    0.000000       0.0000
1237 0.0000    0.000000       0.0000
1238 0.0000    0.006479       0.1763
1239 0.0000    0.007660       0.2084
1240 0.0000    0.009209       0.2506
......
1271 0.0000    0.095361       2.5949
1272 0.0000    0.099645       2.7115

作者
Author: wzkchem5 时间: 2023-10-9 16:17
几乎任何计算，足够高的轨道的能量都是正的。你随便打开一个计算文件（哪怕没有过渡金属），看能量最高的几个轨道就知道了。只不过这些轨道的位置特别靠后，如果你平时只看HOMO-LUMO gap附近的轨道的话，容易看不到而已。
当基组无穷大，而且体系呈电中性或带正电荷的时候，前无穷多个轨道的能量都是负的。但实际计算的基组不是无穷大，所以并不是所有轨道的能量都是负的。

作者
Author: 413 时间: 2023-10-10 09:19

wzkchem5 发表于 2023-10-9 16:17
几乎任何计算，足够高的轨道的能量都是正的。你随便打开一个计算文件（哪怕没有过渡金属），看能量最高的几 ...

好的，谢谢谢谢。

搭车想请教下，scf的过程中会出现这种warning，如何消除呢？

WARNING : small HOMO - LUMO gap : 0.013533
  490    dE -4.704463e-10    4.639610e-06                         (NR MIcro)
  490    dE -4.713915e-10    2.627936e-06                         (NR MIcro)
  490    dE -4.715548e-10    1.830265e-06                         (NR MIcro)
  490    dE -4.716427e-10    9.655585e-07                         (NR MIcro)
  491 -502222.130248927395    5.367579e-05                         (NR MAcro)
WARNING : small HOMO - LUMO gap : 0.013533
  491    dE -8.725774e-11    3.216671e-06                         (NR MIcro)
  491    dE -8.770321e-11    1.691843e-06                         (NR MIcro)
  491    dE -8.779073e-11    1.071979e-06                         (NR MIcro)
  491    dE -8.782351e-11    5.719919e-07                         (NR MIcro)

作者
Author: zjxitcc 时间: 2023-10-10 10:17
本帖最后由 zjxitcc 于 2023-10-10 10:19 编辑

413 发表于 2023-10-10 09:19
好的，谢谢谢谢。

搭车想请教下，scf的过程中会出现这种warning，如何消除呢？

这看起来像是TRAH过程中的警告，是正常信息，不必消除。适当的警告信息可以提醒用户注意检查、注意自己在算什么。要做的是检验SCF收敛后的波函数稳定性，检查SCF收敛后的HOMO-LUMO gap是否仍很小。

作者
Author: 413 时间: 2023-10-10 10:20

zjxitcc 发表于 2023-10-10 10:17
这看起来像是TRAH过程中的警告，是正常信息，不必消除。要做的是检验SCF收敛后的波函数稳定性，检查SCF收 ...

谢谢指导。对的，是trah中的。可是第二列都跑到10e-7了，还不收敛。

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