计算化学公社

标题: 高斯16扫描氢分子多个激发态势能面问题 [打印本页]

作者
Author: nuke 时间: 2019-8-9 03:05
标题: 高斯16扫描氢分子多个激发态势能面问题
本帖最后由 nuke 于 2019-8-9 03:07 编辑

用下面的h2-excited.gjf扫描激发态势能曲线输出了（Summary of the potential surface scan:）两条曲线SCF应该是基态的，CIS应该是第一激发态的。但是其他的呢？如何找到啊谢谢。
如何在这个Summary of the potential surface scan 中多数出几个势能曲线谢谢
---------------------------------------------------------------
h2-excited.gjf
%nprocs=2
%mem=60MW
# scan td=(nstates=4,root=3) b3lyp/6-311++g(d,p) nosymm

Title Card Required

0 1
H
H                1          B1

B1          0.60000000 60 0.010000
-------------------------------------------------------------------
Output
Summary of the potential surface scan:
N    B1       SCF       CIS
----  ---------  -----------  -----------
1    0.6000    -1.16032    -0.54408
2    0.6100    -1.16331    -0.55068
3    0.6200    -1.16598    -0.55684
4    0.6300    -1.16836    -0.56269
5    0.6400    -1.17046    -0.56880
6    0.6500    -1.17230    -0.57626
7    0.6600    -1.17390    -0.57833
8    0.6700    -1.17527    -0.58267
9    0.6800    -1.17642    -0.58697
10    0.6900    -1.17738    -0.59112
11    0.7000    -1.17815    -0.59509
12    0.7100    -1.17874    -0.59891
13    0.7200    -1.17916    -0.60266
14    0.7300    -1.17943    -0.60598
15    0.7400    -1.17956    -0.61037
16    0.7500    -1.17955    -0.61263
17    0.7600    -1.17941    -0.61552
18    0.7700    -1.17915    -0.61828
19    0.7800    -1.17878    -0.62092
20    0.7900    -1.17830    -0.62367
21    0.8000    -1.17772    -0.62607
22    0.8100    -1.17705    -0.62836
23    0.8200    -1.17630    -0.63054
24    0.8300    -1.17546    -0.63264
25    0.8400    -1.17454    -0.63467
26    0.8500    -1.17356    -0.63656
27    0.8600    -1.17251    -0.63836
28    0.8700    -1.17139    -0.64008
29    0.8800    -1.17021    -0.64173
30    0.8900    -1.16899    -0.64330
31    0.9000    -1.16771    -0.64480
32    0.9100    -1.16638    -0.64623
33    0.9200    -1.16501    -0.64759
34    0.9300    -1.16359    -0.64889
35    0.9400    -1.16214    -0.65014
36    0.9500    -1.16065    -0.65133
37    0.9600    -1.15913    -0.65246
38    0.9700    -1.15758    -0.65355
39    0.9800    -1.15600    -0.65458
40    0.9900    -1.15440    -0.65557
41    1.0000    -1.15277    -0.65651
42    1.0100    -1.15112    -0.65740
43    1.0200    -1.14945    -0.65826
44    1.0300    -1.14776    -0.65907
45    1.0400    -1.14606    -0.65985
46    1.0500    -1.14434    -0.66058
47    1.0600    -1.14260    -0.66129
48    1.0700    -1.14086    -0.66195
49    1.0800    -1.13910    -0.66258
50    1.0900    -1.13733    -0.66318
51    1.1000    -1.13556    -0.66375
52    1.1100    -1.13378    -0.66429
53    1.1200    -1.13199    -0.66480
54    1.1300    -1.13020    -0.66528
55    1.1400    -1.12840    -0.66574
56    1.1500    -1.12660    -0.66616
57    1.1600    -1.12480    -0.66656
58    1.1700    -1.12299    -0.66694
59    1.1800    -1.12119    -0.66729
60    1.1900    -1.11938    -0.66761
61    1.2000    -1.11758    -0.66792
----  ---------  -----------  -----------

作者
Author: sobereva 时间: 2019-8-9 03:40
由于你用的是root=3，所以CIS是第3激发态的，而非第1激发态
其它的态的自行用激发能加到基态能量上就得到了

作者
Author: nuke 时间: 2019-8-9 05:37

sobereva 发表于 2019-8-9 03:40
由于你用的是root=3，所以CIS是第3激发态的，而非第1激发态
其它的态的自行用激发能加到基态能量上就得到 ...

感谢回复。
所以说如果要扫描三个激发态的势能曲线，只需改变input文件root=1,2,3，并分别运行。然后手工把summary里面的cis势能面和一起就可以了吧
---------------------------------------------------------------
h2-excited.gjf
%nprocs=2
%mem=60MW
# scan td=(nstates=4,root=1) b3lyp/6-311++g(d,p) nosymm

Title Card Required

0 1
H
H                1          B1

B1          0.60000000 60 0.010000

作者
Author: sobereva 时间: 2019-8-9 06:20

nuke 发表于 2019-8-9 05:37
感谢回复。
所以说如果要扫描三个激发态的势能曲线，只需改变input文件root=1,2,3，并分别运行。然后手 ...

你不嫌麻烦你就这么算

作者
Author: nuke 时间: 2019-8-9 12:37
本帖最后由 nuke 于 2019-8-9 14:17 编辑

sobereva 发表于 2019-8-9 06:20
你不嫌麻烦你就这么算

如大师所说感谢指点迷津。不过不知能否在Summary再加一列能量输出之类，这样感觉更方便些

作者
Author: sobereva 时间: 2019-8-10 08:18

nuke 发表于 2019-8-9 12:37
如大师所说感谢指点迷津。不过不知能否在Summary再加一列能量输出之类，这样感觉更方便些

自行用grep命令把其它激发态的激发能批量提取出来，加上基态能量就行了，非常容易

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