计算化学公社

标题: Gaussian 计算光谱不收敛问题 [打印本页]

作者
Author: juanshenghuang 时间: 2016-5-11 23:10
标题: Gaussian 计算光谱不收敛问题
我计算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

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

作者
Author: juanshenghuang 时间: 2016-5-13 09:12

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

多谢，我试试看

作者
Author: juanshenghuang 时间: 2016-5-13 10:59

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

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

作者
Author: jiliu482 时间: 2016-5-13 11:20

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

不客气，

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

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

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

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

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

作者
Author: juanshenghuang 时间: 2016-5-14 17:53

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

此时 ...

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

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