| 查看: 609 | 回复: 3 | |||
| 当前主题已经存档。 | |||
faqianliu银虫 (小有名气)
|
[交流]
【求助】关于使用RCIS和TD计算荧光的一些疑问，请高手指正
|
||
|
各位高手，经搜索论坛知，一般来说计算化合物的荧光光谱一般分四个步骤： (1)优化基态结构;一般是HF或DFT; (2)优化激发态结构: (a)单重态,可以计算荧光,优化命令 #RCIS/6-31g opt freq=Noraman (b)三重态,可以计算磷光,先优化结构 (3)在以上结构基础上算激发态的单点能 #RCIS/6-31g 或 #B3LYP/6-31g TD (4)输出文件中可以看到结构.用GV可以看到图谱!!! 现在我以讲义上计算甲醛的激发态为例进行计算， 第一步优化忽略； 第二步计算命令如下： %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 问题1： 第二步的第一小步计算#T RCIS/6-31+G(D) Test默认计算三个激发态，计算结果如下： Ground to excited state transition densities written to RWF 633 Excitation energies and oscillator strengths: Excited State 1: Singlet-A2 4.7282 eV 262.22 nm f=0.0000 8 -> 9 0.51059 8 -> 13 0.45471 8 -> 17 0.14148 This state for optimization and/or second-order correction. Copying the excited state density for this state as the 1-particle RhoCI density. Excited State 2: Singlet-B2 8.8305 eV 140.41 nm f=0.0273 8 -> 10 0.63905 8 -> 11 -0.25389 8 -> 18 0.11389 Excited State 3: Singlet-B2 9.5206 eV 130.23 nm f=0.0555 8 -> 10 0.26128 8 -> 11 0.63162 8 -> 18 -0.13313 这个是我们传统意义上的吸收光谱，也就是UV-VIS谱； 对于# RCIS(Root=1,Read)/6-31+G(D) Guess=Read Opt Freq Test Geom=AllCheck讲义上只对第一个激发态进行优化，这样可以得到第一激发态的优化坐标，为什么不对第二和第三个激发态进行优化？如果想计算甲醛的发射光谱，是否只需要计算第一激发态，还是应该计算强度f最大的态？ 以我使用TD-DFT计算金属配合物的吸收光谱为例，我计算的是150个态，其中第56个激发态的强度f最大，如下所示： Excited State 59: ?Spin -A 3.4549 eV 358.86 nm f=0.1135 这样如果要计算配合物的荧光，我是否是用# RCIS(Root=59,Read)/6-31+G(D) Guess=Read Opt Freq Test Geom=AllCheck，来对配合物的第59个激发态进行结构优化，然后再TD计算配合物的荧光呢？还是有其它的计算方法？ 问题2： 仍以上面的甲醛第一激发态优化为例，RCIS优化结束后，可以得到第一激发态的结构， 进行第三步，#B3LYP/6-31g TD，进行计算，结果如下： Ground to excited state transition densities written to RWF 633 Excitation energies and oscillator strengths: Excited State 1: Singlet-A2 3.6799 eV 336.92 nm f=0.0000 8 -> 9 0.67898 This state for optimization and/or second-order correction. Copying the excited state density for this state as the 1-particle RhoCI density. Excited State 2: Singlet-B2 6.9416 eV 178.61 nm f=0.0164 8 -> 10 0.69949 Excited State 3: Singlet-B2 7.6719 eV 161.61 nm f=0.0289 8 -> 11 0.69881 这个计算的是发射光谱吗？怎么感觉还是还是基态到激发态的跃迁，和吸收光谱一样呢，当然波长都变化了。。。。而且这个只是计算的第一激发态的光谱，如果使用# RCIS(Root=2,Read)/6-31+G(D) Guess=Read Opt Freq Test Geom=AllCheck来进行激发态结构优化的话，结果是不一样的，应该怎样计算？ 请高手指教！ |
回复此楼» 猜你喜欢
分子筛催化方向 一篇cej共一第一 可以申请上海什么学校呢
已经有4人回复
-
求助共晶分子筛的电镜图分析，有偿
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有288人回复
-
中南大学化学化工学院招收电池储能、理论计算方向博士生
已经有16人回复
-
请问哪里可以代烧带压力的炉子？
已经有6人回复
-
关于合成MIL-125的求助
已经有0人回复
-
双非本硕光催化方向想读博，请各位老师看看我！
已经有2人回复
-
瑞典林雪平大学范立洲课题组现招聘博士后、岗位制博士生
已经有0人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有22人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
yjcmwgk
禁虫 (文坛精英)
密度泛函·小卒
- QC强帖: 38
- 应助: 165 (高中生)
- 贵宾: 16.999
- 金币: 555.4
- 散金: 39880
- 红花: 395
- 沙发: 2
- 帖子: 10406
- 在线: 1528.1小时
- 虫号: 448368
- 注册: 2007-11-01
- 专业: 理论和计算化学
- 管辖: 科研工具资源
2楼2010-03-14 12:52:11
faqianliu
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1042.7
- 散金: 75
- 红花: 1
- 帖子: 165
- 在线: 180.6小时
- 虫号: 263012
- 注册: 2006-07-01
- 专业: 无机纳米化学
3楼2010-03-15 09:34:10
yjcmwgk
禁虫 (文坛精英)
密度泛函·小卒
- QC强帖: 38
- 应助: 165 (高中生)
- 贵宾: 16.999
- 金币: 555.4
- 散金: 39880
- 红花: 395
- 沙发: 2
- 帖子: 10406
- 在线: 1528.1小时
- 虫号: 448368
- 注册: 2007-11-01
- 专业: 理论和计算化学
- 管辖: 科研工具资源
4楼2010-03-15 10:41:31
10