激发态计算 可以完成的任务有:1.圆明园事件用CIS TD 和ZINDO 方法计算紫外吸收光谱;2.可用CIS 方法优化 退组词激发态,研究光化学反应中反应物的构型与机理的关系;3.可以计算发射光谱。 与上面的方法相比,CASSCF 方法的计算效果要更好。不过这种方法对计算者的要求很高。希望感兴趣的朋友仔细的研究高斯自带的例子和有关的讨论再做尝试。磨刀不误砍柴功。需要特别提出的是 CIS 计算对内存的要求也是比较高的。 激发态是指分子体系的稳定的,高能级的电子结构状态。 模拟激发态的理论很少,在Gaussian中,提供的是结构相关方法(Configuration Interactionapproach),采用Hartree-Fock 单置换,所以也被称为CI-Singles。这一方法被认为是很多分子激发态的零级处理的有效方法。在定性上有比较好的结果,但在定量方面,并不总是好的。和Hartree-Fock 方法一样,CI-Singles方法也不昂贵。ZINDO 是半经验方法,用的不多。TD 是含时密度泛函方法,计算结果比CIS 方法要好。据说该方法不能用于优化。 下面是激发态计算的关键词 CIS 进行CI-Silgles 激发态计算,这是方法关键词,在其前面可以加上 R 和U 作为闭壳层分子和开壳层分子的标识。 CIS=(Root=n) 确定研究哪一个激发态,默认为1 CIS=(NStates=n) 确定研究的激发态的数量。默认为3 CIS=50-50 同时计算单重态和三重态。Triplets 选项为计算单重态,默认 为单重态。当选择Triplets 时,一定要设置NStates 数量。 CIS=Read买豆腐 从临时文件中读入初始猜测。用于利用前一个工作中的结果 进行计算,与Guess=Read 和Geom=Check 一起使用。 Density=Current 要求计算激发态的布局分析 Pop=Reg 要求较详细的布局分析结果,包括分子轨道系数 例子1:乙烯的激发态 #T RCIS(NStates=2,50-50)/6-31+G(d) Test Ethylene Excited States 0,1 C C,1,R(C-C) H,1,R(C-H),2,A(CCH) H,1,R(C-H),2,A(CCH),3,180. H,2,R(C-H),1,A(CCH),3,0. H,2,R(C-H),1,A(CCH),3,180. R(C-C)=1.31696 R(C-H)=1.07601 A(CCH)=121.81335 由于即需要单重态,也需要三重态,所以设置50-50,每个状态计算两个。 注意这里使用了弥散基组,弥散基组在计算过渡态方面有好的结果。 喷泉英文Excited State 1:Triplet-B1U 3.7768eV 328.27 nm f=0.000 多重度 对称性 激发能 振荡强度 8 --> 11 0.52952 跃迁轨道数 激发态波函数系数 8 --> 17 -0.45942 This State for optimization and/or cond-order correction王思想: Total Energy,E(Cis) = -77.8969983928 激发态能量 Copying the Cisingles density for this state as the 1-particle RhoCI density。 计算结果是 1 2 3 4 对称性 3B1u 3B3u 1B3u 1B1u 能量 3.78 7.43 7.83 7.98 实验值 4.36 7.66 6.98 7.15 总的说计算的结果与实验值符合的不错,更加精确的结果就需要大的基组计算了。 |
例子2:甲醛的激发态优化 在Gaussian 中可以进行激发态的优化和频率分析。步骤是,首先进行能量计算得到激发态, 然后在此结构基础上进行结构优化和频率分析。 %Chk=c2ho_es #T RCIS/6-31+G(D) Test Formaldehyde (C2V) Vertical Excited States 0,1 C O,1,AB H,1,AH,2,HAB H,1,AH,2,HAB,3,180. AB=1.18429 AH=1.09169 HAB=122.13658 --Link1-- %Chk=c2ho_es %NoSave # RCIS(Root=1,Read)/6-31+G(D) Guess=Read Opt Freq Test Geom=AllCheck 第一个步骤计算最低的三个激发态,第二步骤是利用其结果进行结构优化和频率分析. 在本例的计算中,得到的激发态出现了虚频,其结构显示碳原子要离开原子平面。 由于原来的分子结构中,所有原子都在同一平面上,所以,也得到一个共平面的激发态结构。 这样,就要设法产生一个正确的几何结构。(激发态的结构优化目前没有包含内坐标冗余, 所以在初始结构中定义一个平面结构,得到的激发态结构就也是共平面的)。 例3 甲醛激发态优化 ,确定甲醛第一激发态的结构,比较红外光谱。 #T RCIS(Root=1)/6-31+G(D) Opt Freq Density=Current Test Formaldehyde n-->pi* Cs Optimization 昙花什么时候开花0,1 今日头条视频下载C O,1,R2 X,1,1.,2,A3 H,1,R4,3,A4,2,90. H,1,R4,3,A4,2,-90. R2=1.25546079 R4=1.08213788 A3=145. A4=60. 这里使用了虚原子,使得对于分子的描述更加简单,这个分子也不再是共平面的结构。计算采用校正因子0.8929,结果如下: 计算值 校正值 实验值 平面外弯曲 495 442 683 CH2 摇摆 978 873 898 CH2 剪式 1426 1273 1290 CO 伸缩 1647 1471 1173 对称CH伸缩 3200 2857 2847 反对称CH伸缩3295 2942 2968 基本与实验值相符 |
8个字的成语 |
本文发布于:2023-08-01 02:35:49,感谢您对本站的认可!
本文链接:https://www.wtabcd.cn/fanwen/fan/89/1103813.html
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。
留言与评论(共有 0 条评论) |