| 查看: 2930 | 回复: 6 | ||
[求助]
怎样用高斯在计算配合物的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)之类的？ |
回复此楼» 猜你喜欢
宁夏大学化学工程与技术专业招收博士研究生
已经有1人回复
-
南京大学软物质计算模拟方向硕博连读/直博生招生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有204人回复
-
微纳hub，微纳加工的宝库
已经有0人回复
-
微纳加工里的 lift-off 工艺
已经有0人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有1人回复
-
美国杰克逊州立大学周剑歌教授的化学计算与AI实验室诚聘博士后研究员
已经有4人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有1人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
Tc数据拟合
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用高斯模拟的吸收光谱中，有三个峰出现，请问怎么确定最大吸收波长呢？
已经有3人回复
- 用高斯计算磷光发射波长时，振子强度是0，实验上可以测量磷光光谱的，求交流
已经有19人回复
2楼2014-08-15 20:11:12
ZJboy
木虫 (著名写手)
- 应助: 34 (小学生)
- 金币: 5707.3
- 散金: 92
- 红花: 8
- 帖子: 1315
- 在线: 332.8小时
- 虫号: 430962
- 注册: 2007-08-11
- 性别: GG
- 专业: 配位化学
3楼2014-08-16 11:23:00
yumu0411
金虫 (著名写手)
打酱油
- 应助: 24 (小学生)
- 金币: 4052.7
- 散金: 5233
- 红花: 9
- 帖子: 1496
- 在线: 766.3小时
- 虫号: 864757
- 注册: 2009-10-07
- 性别: GG
- 专业: 有机分子功能材料化学
4楼2014-08-16 16:28:37
ZJboy
木虫 (著名写手)
- 应助: 34 (小学生)
- 金币: 5707.3
- 散金: 92
- 红花: 8
- 帖子: 1315
- 在线: 332.8小时
- 虫号: 430962
- 注册: 2007-08-11
- 性别: GG
- 专业: 配位化学
5楼2014-08-17 16:08:26
6楼2014-09-29 08:37:29
7楼2016-11-29 09:21:02
