24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

szhshuan

银虫 (正式写手)

应助: 3 (幼儿园)
金币: 44.1
散金: 527
红花: 2
帖子: 403
在线: 169.5小时
虫号: 930547
注册: 2009-12-19
性别: GG
专业: 有机合成

[求助] 怎样用高斯在计算配合物的UV光谱时，只计算某一波长范围内的紫外光谱跃迁已有2人参与

我要计算一个Eu的配合物的UV光谱，初始输入文件如下：
%chk=Eu-2p1.chk
%mem=2000MW
%nproclinda=1
%nprocshared=12
#p b3lyp/genecp geom=check guess=read scf=(maxcyc=128,xqc) td(nstates=10,solvent=acetonitrile)

title

3 1

C N H 0
6-311G(d,p)
****
Eu 0
S 5 1.00
70059.4200000 0.000097
10776.2350000 0.000707
2482.4901000 0.003362
702.1526000 0.010206
216.7926000 0.016597
S 1 1.00
41.4950000 1.
S 1 1.00
29.4969000 1.
S 1 1.00
15.0309000 1.
S 1 1.00
3.7772000 1.
S 1 1.00
1.9219000 1.
S 1 1.00
0.7131000 1.
S 1 1.00
0.3097000 1.
S 1 1.00
0.0549000 1.
S 1 1.00
0.0227000 1.
P 6 1.00
4058.0207000 0.000048
1023.3259000 0.000357
339.3852000 0.001687
122.3485000 0.005693
29.0558000 0.072398
20.6269000 -0.046807
P 1 1.00
14.6179000 1.
P 1 1.00
5.4444000 1.
P 1 1.00
2.7414000 1.
P 1 1.00
1.2283000 1.
P 1 1.00
0.5769000 1.
P 1 1.00
0.2489000 1.
P 1 1.00
0.0800000 1.
D 6 1.00
389.0841000 0.000534
117.7602000 0.004258
44.9869000 0.017160
22.8599000 0.012427
7.8331000 0.270416
4.1134000 0.445272
D 1 1.00
2.0999000 1.
D 1 1.00
1.0043000 1.
D 1 1.00
0.3422000 1.
D 1 1.00
0.1051000 1.
F 5 1.00
128.2341000 0.003424
46.2685000 0.034313
20.7366000 0.123194
9.4479000 0.251824
4.3149000 0.351292
F 1 1.00
1.9121000 1.
F 1 1.00
0.7897000 1.
F 1 1.00
0.2859000 1.
G 4 1.00
20.7366000 0.026516
9.4479000 0.099066
4.3149000 0.213378
1.9121000 0.453713
G 1 1.00
0.7897000 1.
G 1 1.00
0.2859000 1.
****

Eu 0
Eu-ECP 5 28
H-K
1
2 1.000000 0.000000
S-H
1
2 23.471384 607.659331
P-H
1
2 16.772479 264.385476
D-H
1
2 13.981343 115.381375
F-H
1
2 23.962888 -49.400794
G-H
1
2 21.232458 -26.748273


而输出文件中UV光谱部分的各阶段波长对应的峰强度都为零：
Excitation energies and oscillator strengths:

Excited State 1:    Singlet-A    -0.4834 eV -2565.04 nm  f=-0.0000  <S**2>=0.000
   138 -> 164       0.10362
   138 -> 165    -0.47261
   139 -> 164    -0.47290
   139 -> 165    -0.10402
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) =  -2432.12358821
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State 2:    Singlet-A    -0.3962 eV -3129.05 nm  f=-0.0000  <S**2>=0.000
   138 -> 164    -0.48398
   138 -> 165    -0.10551
   139 -> 164    -0.10641
   139 -> 165       0.48077

Excited State 3:    Singlet-A    -0.0647 eV -19163.44 nm  f=-0.0000  <S**2>=0.000
   122 -> 166    -0.63717
   124 -> 166       0.11421
   131 -> 166    -0.57690
   137 -> 166    -0.94432
   151 -> 166       0.26579
   154 -> 166       0.16082
   157 -> 166    -0.16119
   122 <- 166    -0.54377
   131 <- 166    -0.49209
   137 <- 166    -0.80440
   151 <- 166       0.22368
   154 <- 166       0.13521
   157 <- 166    -0.13525

Excited State 4:    Singlet-A    0.1869 eV 6634.23 nm  f=0.0000  <S**2>=0.000
   122 -> 165       0.11670
   131 -> 165       0.10632
   137 -> 165       0.17467
   138 -> 164    -0.18373
   138 -> 165       0.34442
   138 -> 167       0.34843
   139 -> 164    -0.34327
   139 -> 165    -0.18731
   139 -> 167       0.22262
   138 <- 167       0.19899
   139 <- 167       0.12713

Excited State 5:    Singlet-A    0.1888 eV 6566.04 nm  f=0.0000  <S**2>=0.000
   122 -> 164       0.11676
   131 -> 164       0.10541
   137 -> 164       0.17338
   138 -> 164    -0.34230
   138 -> 165    -0.18490
   138 -> 167    -0.22200
   139 -> 164       0.18624
   139 -> 165    -0.34538
   139 -> 167       0.34731
   138 <- 167    -0.12599
   139 <- 167       0.19710

Excited State 6:    Singlet-A    1.2663 eV  979.12 nm  f=0.0000  <S**2>=0.000
   137 -> 165       0.12582
   138 -> 166    -0.10799
   138 -> 167    -0.15740
   139 -> 166    -0.64730
   139 <- 166       0.13076

Excited State 7:    Singlet-A    1.2668 eV  978.72 nm  f=0.0000  <S**2>=0.000
   137 -> 164    -0.12680
   138 -> 166    -0.64625
   139 -> 166       0.10759
   139 -> 167       0.15929
   138 <- 166       0.13045

Excited State 8:    Singlet-A    1.3734 eV  902.75 nm  f=0.0000  <S**2>=0.000
   122 -> 164    -0.22921
   131 -> 164    -0.20943
   137 -> 164    -0.34930
   138 -> 166       0.25995
   139 -> 167       0.41007
   151 -> 164       0.11683

Excited State 9:    Singlet-A    1.3741 eV  902.28 nm  f=0.0000  <S**2>=0.000
   122 -> 165    -0.22961
   131 -> 165    -0.20987
   137 -> 165    -0.34995
   138 -> 167       0.41001
   139 -> 166    -0.25781
   151 -> 165       0.11713

Excited State  10:    Singlet-A    2.3613 eV  525.06 nm  f=0.0001  <S**2>=0.000
   137 -> 164    -0.11054
   138 -> 164    -0.11933
   139 -> 165    -0.11984
   139 -> 167    -0.13904
   162 -> 164       0.38668
   162 -> 165       0.23650
   163 -> 164    -0.26629
   163 -> 165       0.38025
SavETr:  write IOETrn= 770 NScale= 10 NData=  16 NLR=1 LETran=    190.
Leave Link  914 at Sun Aug  3 15:12:28 2014, MaxMem= 2097152000 cpu:  160552.9
(Enter /public/software/gaussian/g09c01/g09/l601.exe)
Copying SCF densities to generalized density rwf, IOpCl= 0 IROHF=0.

而该配合物实验数据中吸收在200-350nm（22.png）。
请问我是否可以只计算200-350nm之间的紫外光谱，请大家帮我编辑一下输入文件，谢谢。或者是在原有chk基础上改为td(add=90,full)之类的？

回复此楼