CIS方法激发态计算

更新时间:2023-08-01 02:35:49 阅读: 评论:0

CIS方法激发态计算1(带高斯中的例子)
2008-09-25 17:21
激发态计算
可以完成的任务有: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 1Triplet-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 EnergyE(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
总的说计算的结果与实验值符合的不错,更加精确的结果就需要大的基组计算了。
CIS方法激发态计算2(带高斯中的例子)
2008-09-25 17:22
例子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 条评论)
   
验证码:
推荐文章
排行榜
Copyright ©2019-2022 Comsenz Inc.Powered by © 专利检索| 网站地图