使用Multiwfn绘制电荷转移光谱(CTS)为什么会出现波长正数到负数的情况

ddddnight

下面是算激发态的gif文件内容，请问各位老师怎么回事？%nprocshared=18
%mem=24GB
%chk=D:\gaussian\calculation\mix base\02state.chk
#p pbe1pbe/genecp IOp(9/40=4) TD(nstates=50) SCF=QC

02

3 2
C                 -3.03705300    0.97340900    0.86841300
C                 -1.93886200    1.08942200    0.05456400
C                 -1.81626200    0.38065500   -1.17419700
C                 -2.82014800   -0.51526700   -1.54915800
C                 -3.94954800   -0.66343200   -0.76319900
C                 -4.12417800    0.10991800    0.46698300
H                 -1.15260600    1.77674500    0.35440100
H                 -2.75636600   -1.08489800   -2.46863900
C                  0.67064900    0.55143700   -1.35207600
C                  1.66213800    1.53693900   -1.34119100
C                  0.94674700   -0.70367300   -0.70393700
C                  2.91121000    1.29745300   -0.78814000
H                  1.46429700    2.52085500   -1.75153300
C                  2.17902300   -0.96181500   -0.17125800
H                  0.17625100   -1.47042500   -0.67942100
C                  3.23851400   -0.01492000   -0.26913900
C                 -0.64297200    0.71828500   -2.09382300
C                 -0.55161300   -0.22763300   -3.32184700
H                  0.32564400    0.02655200   -3.92004800
H                 -0.48860900   -1.28315700   -3.04365700
H                 -1.43043000   -0.07927700   -3.95553100
C                 -0.86068900    2.15804600   -2.60948700
H                 -0.84730300    2.89931700   -1.80586000
H                 -0.10239700    2.41619200   -3.35400600
H                 -1.83223100    2.21817400   -3.10381500
O                  4.44136600   -0.38271600    0.15457600
O                 -5.25272600   -0.08115200    1.04198300
Br                -5.31494400   -1.82839700   -1.28552600
Br                -3.05393000    1.96126800    2.49220800
Br                 4.10112100    2.73835100   -0.59829900
Br                 2.51976900   -2.69269800    0.60694100
C                 -5.88443900    0.55204700    2.19719700
H                 -6.88875600    0.13746500    2.19093500
H                 -5.34296300    0.26535900    3.09792300
H                 -5.89592500    1.63156800    2.05308300
C                  5.70636000    0.00099900   -0.53070100
H                  5.48156300    0.36193700   -1.53243500
H                  6.20839400    0.75217600    0.07515400
H                  6.27750500   -0.92639400   -0.58385700
Fe                 4.78234700   -1.93810800    1.72324900

Fe Br
SDD
****
C H O
6-311G**
****

Fe Br
SDD

wzkchem5

波长是负的是因为SCF波函数不稳定。原因是（1）自旋多重度错误；（2）结构不合理。
自旋多重度的问题，建议找本无机化学（配位化学）的教材看看，学会判断这个配合物是高自旋的还是低自旋的，然后设定正确的自旋多重度。
结构的问题，必须对Fe做显式溶剂化，并且做构型搜索、优化结构，得到没有虚频的结构以后才能算。优化结构、算频率的时候建议加溶剂模型，必须加色散校正，建议做stable=opt

ddddnight

wzkchem5 发表于 2022-9-1 23:43
波长是负的是因为SCF波函数不稳定。原因是（1）自旋多重度错误；（2）结构不合理。
自旋多重度的问题，建 ...

好的，谢谢老师

sobereva

不建议轻易用SCF=QC，收敛到不稳定波函数的概率会大为提高
遇到SCF不收敛的正确考虑方式看
解决SCF不收敛问题的方法
http://sobereva.com/61

ddddnight

sobereva 发表于 2022-9-2 08:11
不建议轻易用SCF=QC，收敛到不稳定波函数的概率会大为提高
遇到SCF不收敛的正确考虑方式看
解决SCF不收敛 ...

老师您好。我所有的方式都试过了，都不行，最后选用的scf＝qc来解决收敛的，结果就成了这样了

wzkchem5

ddddnight 发表于 2022-9-2 01:25
老师您好。我所有的方式都试过了，都不行，最后选用的scf＝qc来解决收敛的，结果就成了这样了

如果很难收敛，有可能是结构不合理或者自旋多重度不合理，一定要反复检查结构、自旋多重度，不能盲目加收敛选项

ddddnight

wzkchem5 发表于 2022-9-2 15:17
如果很难收敛，有可能是结构不合理或者自旋多重度不合理，一定要反复检查结构、自旋多重度，不能盲目加收 ...

谢谢老师，应该是结构有些不合理，因为我换了几个自旋多重度都不能收敛。