计算化学公社

标题: 使用Multiwfn绘制电荷转移光谱(CTS)为什么会出现波长正数到负数的情况 [打印本页]

作者
Author: ddddnight 时间: 2022-9-1 22:50
标题: 使用Multiwfn绘制电荷转移光谱(CTS)为什么会出现波长正数到负数的情况
下面是算激发态的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

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

作者
Author: ddddnight 时间: 2022-9-2 07:24

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

好的，谢谢老师

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

作者
Author: ddddnight 时间: 2022-9-2 08:25

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

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

作者
Author: wzkchem5 时间: 2022-9-2 15:17

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

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

作者
Author: ddddnight 时间: 2022-9-2 15:42

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

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

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