Gaussian 计算光谱不收敛问题

juanshenghuang

我计算UV-vis老是出现收敛失败，不知道什么原因，本人新手，还望各位多多指教了。谢谢了
以上是我计算用gif文件。

# freq ucis=(root=1)/lanl2dz geom=connectivity density=current

Formaldehyde n-->pi* Cs Optimization

0 2
Br                 9.07500000    8.72900000    5.97800000
Br                15.12500000    8.72900000    5.97800000
Br                 9.07500000   14.77900000    5.97800000
Br                15.12500000   14.77900000    5.97800000
Br                 6.05000000    8.73000000    8.96800000
Br                 9.07500000    6.05300000    9.11000000
Br                 9.07500000    8.72900000   12.02800000
Pb                 9.07500000    9.12300000    9.05500000
Br                12.10000000    8.73000000    8.96800000
Br                15.12500000    6.05300000    9.11000000
Br                15.12500000    8.72900000   12.02800000
Pb                15.12500000    9.12300000    9.05500000
Br                18.15000000    8.73000000    8.96800000
Br                 6.05000000   14.78000000    8.96800000
Br                 9.07500000   12.10300000    9.11000000
Br                 9.07500000   14.77900000   12.02800000
Pb                 9.07500000   15.17300000    9.05500000
Br                12.10000000   14.78000000    8.96800000
Br                15.12500000   12.10300000    9.11000000
Br                15.12500000   14.77900000   12.02800000
Pb                15.12500000   15.17300000    9.05500000
Br                18.15000000   14.78000000    8.96800000
Br                 9.07500000   18.15300000    9.11000000
Br                15.12500000   18.15300000    9.11000000
Br                 6.05000000    8.73000000   15.01800000
Br                 9.07500000    6.05300000   15.16000000
Br                 9.07500000    8.72900000   18.07800000
Pb                 9.07500000    9.12300000   15.10500000
Br                12.10000000    8.73000000   15.01800000
Br                15.12500000    6.05300000   15.16000000
Br                15.12500000    8.72900000   18.07800000
Pb                15.12500000    9.12300000   15.10500000
Br                18.15000000    8.73000000   15.01800000
Br                 6.05000000   14.78000000   15.01800000
Br                 9.07500000   12.10300000   15.16000000
Br                 9.07500000   14.77900000   18.07800000
Pb                 9.07500000   15.17300000   15.10500000
Br                12.10000000   14.78000000   15.01800000
Br                15.12500000   12.10300000   15.16000000
Br                15.12500000   14.77900000   18.07800000
Pb                15.12500000   15.17300000   15.10500000
Br                18.15000000   14.78000000   15.01800000
Br                 9.07500000   18.15300000   15.16000000
Br                15.12500000   18.15300000   15.16000000
C                 12.10000000   11.43000000   12.10800000
N                 12.10000000   12.92500000   12.12600000
H                 13.01500000   11.05230000   12.61680000
H                 12.09690000   11.07110000   11.05440000
H                 11.20310000   11.05190000   12.64780000
H                 11.24860000   13.27640000   11.65260000
H                 12.92730000   13.27140000   11.64010000
H                 12.10000000   13.30700000   13.09800000

1 8 1.0
2 12 1.0
3 17 1.0
4 21 1.0
5 8 1.0
6 8 1.0
7 28 1.0 8 1.0
8 9 1.0 15 1.0
9 12 1.0
10 12 1.0
11 12 1.0 32 1.0
12 19 1.0 13 1.0
13
14 17 1.0
15 17 1.0
16 37 1.0 17 1.0
17 18 1.0 23 1.0
18 21 1.0
19 21 1.0
20 41 1.0 21 1.0
21 22 1.0 24 1.0
22
23
24
25 28 1.0
26 28 1.0
27 28 1.0
28 29 1.0 35 1.0
29 32 1.0
30 32 1.0
31 32 1.0
32 39 1.0 33 1.0
33
34 37 1.0
35 37 1.0
36 37 1.0
37 38 1.0 43 1.0
38 41 1.0
39 41 1.0
40 41 1.0
41 44 1.0 42 1.0
42
43
44
45 46 1.0 47 1.0 48 1.0 49 1.0
46 50 1.0 51 1.0 52 1.0
47
48
49
50
51
52

jiliu482

我们算吸收光谱的时候很少和频率一起算，也很少用到CIS了。一般直接TDDFT就好了
比如 #p B3LYP/6-31G(d)  scf  gfinput gfprint  td=(singlets,Nstates=40)
如果三态把singlets换成triplets，加溶剂用scrf关键词

juanshenghuang

jiliu482 发表于 2016-5-12 09:28
我们算吸收光谱的时候很少和频率一起算，也很少用到CIS了。一般直接TDDFT就好了
比如 #p B3LYP/6-31G(d)   ...

多谢，我试试看

juanshenghuang

jiliu482 发表于 2016-5-12 09:28
我们算吸收光谱的时候很少和频率一起算，也很少用到CIS了。一般直接TDDFT就好了
比如 #p B3LYP/6-31G(d)   ...

因为我得原子貌似超出了6-31G(d)所包含的原子范围，所以换成了LanL2DZ基组初步算了下，其它跟你的设置一样，得到了结果，谢谢你了

jiliu482

juanshenghuang 发表于 2016-5-13 10:59
因为我得原子貌似超出了6-31G(d)所包含的原子范围，所以换成了LanL2DZ基组初步算了下，其它跟你的设置一 ...

不客气，

sobereva

当前体系不要直接写个lanl2dz就完事，此基组对于C、H太烂，应当对C、H起码用6-31G*，其它用Lanl2DZ

此时算光谱的目的，关键词通常写诸如PBE1PBE/genecp TD就完了，默认root=1。

CIS早已过时，这种体系一定要用TDDFT，见
乱谈激发态的计算方法
http://sobereva.com/265

你若不打算得到激发态波函数就不要写density=current，否则会多花很多时间。

除非你想验证激发态优化完有没有虚频，或想得到激发态热力学数据，或要计算振动分辨的电子光谱，否则甭写freq，TDDFT做freq对这样的体系完全做不动。

juanshenghuang

sobereva 发表于 2016-5-14 15:22
当前体系不要直接写个lanl2dz就完事，此基组对于C、H太烂，应当对C、H起码用6-31G*，其它用Lanl2DZ

此时 ...

太感谢sobereva大神的回复，真没想到您能有时间回复，我刚入门不久，在不断看您的帖子学习，谢谢了