计算化学公社

 找回密码 Forget password
 注册 Register
Views: 18506|回复 Reply: 9
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] 关于激发态的几个小问题

[复制链接 Copy URL]

177

帖子

0

威望

1567

eV
积分
1744

Level 5 (御坂)

最近开始做一些跟激发态有关的计算,有几个问题一直搞不太明白,请大家指点一下:
1. 激发态的概念
问题:激发态跟分子轨道的关系?
个人理解,激发态是分子在受到各种激发源(e.g., 光,电,热等)激发后,从基态势能面进入某一个激发态势能面。一直有个搞不清的概念,就是激发态的电子态是不是基态和一系列激发态电子组态的线性组合(类似于组态相互作用的公式)?
比如紫外-可见光谱中的轨道跃迁归属中,常常出现n->π*跃迁等说法,这是不是指分子受到一定波长光激发后,电子从某一占据轨道(n)跃迁到某一空轨道(π*)这样产生的电子组态占主导?其余的只占很小部分?
高斯的TD光谱计算结果文件中会出现如下图所示的归属,其中的27A ->106A我想就是指这样的电子跃迁,后面的0.00119是不是跃迁概率?可以理解为各种组态的系数吗?


2. 在高斯的TD或者CIS计算中,选项NStates是指激发组态的个数吗?
问题可能比较粗浅,求解答,不胜感激!

1.png (12.92 KB, 下载次数 Times of downloads: 65)

UV轨道

UV轨道

6万

帖子

99

威望

5万

eV
积分
120132

管理员

公社社长

2#
发表于 Post on 2015-7-15 11:27:06 | 只看该作者 Only view this author
1 理解没什么问题。
那些数值就是相应轨道跃迁对应的单激发组态的系数,对应的贡献值的计算见此文
电子激发任务中轨道跃迁贡献的计算
http://sobereva.com/230

另外建议参看
图解电子激发的分类
http://sobereva.com/284
跃迁密度分析方法-自然跃迁轨道(NTO)简介
http://sobereva.com/91

2 nstates是指算多少激发态,并非是激发组态的数目。对于TDDFT、CIS这些计算,激发态波函数由单激发组态函数展开,单激发组态函数的数目是占据轨道数目(不计冻芯的轨道)乘以空轨道数目。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

177

帖子

0

威望

1567

eV
积分
1744

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2015-7-15 11:54:19 | 只看该作者 Only view this author
sobereva 发表于 2015-7-15 11:27
1 理解没什么问题。
那些数值就是相应轨道跃迁对应的单激发组态的系数,对应的贡献值的计算见此文
电子激 ...

明白了,感谢Sob!
感觉涉及到激发态的计算都比较复杂,失败率较高,CASSCF看起来比较主流和高端,不过高斯似乎做的不怎么好,而且操作也比较复杂,期待Sob大神的《实用量子化学》里能够多些相关的介绍!

6万

帖子

99

威望

5万

eV
积分
120132

管理员

公社社长

4#
发表于 Post on 2015-7-15 13:35:50 | 只看该作者 Only view this author
激发态计算并不复杂。CASSCF、CASPT2等绝大多数情况不需要,TDDFT还是主流,能够应付大多问题。如果研究小分子时要求精度更高点,就用EOM-CCSD就很好,也很省事,高斯能直接做。只有静态相关强、涉及势能面交叉、双激发等等情况,才必须得用多参考方法处理。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

177

帖子

0

威望

1567

eV
积分
1744

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2015-7-15 14:29:48 | 只看该作者 Only view this author
sobereva 发表于 2015-7-15 13:35
激发态计算并不复杂。CASSCF、CASPT2等绝大多数情况不需要,TDDFT还是主流,能够应付大多问题。如果研究小 ...

谢谢Sob,努力学习中

177

帖子

0

威望

1567

eV
积分
1744

Level 5 (御坂)

6#
 楼主 Author| 发表于 Post on 2015-7-17 11:02:05 | 只看该作者 Only view this author
sobereva 发表于 2015-7-15 13:35
激发态计算并不复杂。CASSCF、CASPT2等绝大多数情况不需要,TDDFT还是主流,能够应付大多问题。如果研究小 ...

最近在算荧光谱,在优化得到激发态Minimum后,用此结构做基态的TD计算得到荧光谱,听一个朋友说这一步需要将NStates设为1以得到S1态,我不清楚是不是该这样?求解

177

帖子

0

威望

1567

eV
积分
1744

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2015-7-17 11:29:39 | 只看该作者 Only view this author
sobereva 发表于 2015-7-15 13:35
激发态计算并不复杂。CASSCF、CASPT2等绝大多数情况不需要,TDDFT还是主流,能够应付大多问题。如果研究小 ...

我仔细想想朋友的意思应该是指要设定Root=1吧,NStates=1我总觉得不合理,求Sob释疑

6万

帖子

99

威望

5万

eV
积分
120132

管理员

公社社长

8#
发表于 Post on 2015-7-17 11:30:31 | 只看该作者 Only view this author
ter20 发表于 2015-7-17 11:02
最近在算荧光谱,在优化得到激发态Minimum后,用此结构做基态的TD计算得到荧光谱,听一个朋友说这一步需 ...


不对。
首先“基态的TD计算”这句话有问题。
而且nstates根本不用变,默认的nstates=3就正合适。虽然只需要计算S0-S1的激发能和振子强度,但不要把nstates正好只设到1,这样算得容易不太准。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

177

帖子

0

威望

1567

eV
积分
1744

Level 5 (御坂)

9#
 楼主 Author| 发表于 Post on 2015-7-17 15:14:44 | 只看该作者 Only view this author
sobereva 发表于 2015-7-17 11:30
不对。
首先“基态的TD计算”这句话有问题。
而且nstates根本不用变,默认的nstates=3就正合适。虽然 ...

1. “基态的TD计算”我是指就用激发态优化得到的结构,做TD计算(其实就是算UV谱),关键字比如:
# PBE1PBE/6-31G(d) TD(Singlet,Direct) geom=allcheck guess=read
我想应该是对的吧,我之前的表达有些问题

2. 过去我算紫外光谱的时候一般NStates都设到了40,因为怕丢掉一些态,其实平时我们只关心最大吸收波长附近的三四个吸收峰,低波长部分(比如<300nm)都不怎么关心,所以按照您的意思,其实NStates根本不需要设那么大对吧?

6万

帖子

99

威望

5万

eV
积分
120132

管理员

公社社长

10#
发表于 Post on 2015-7-17 19:32:59 | 只看该作者 Only view this author
ter20 发表于 2015-7-17 15:14
1. “基态的TD计算”我是指就用激发态优化得到的结构,做TD计算(其实就是算UV谱),关键字比如:
# PBE ...

1 对。direct没必要写
2 设多大看体系特征、原子数和研究目的。比如有机大共轭体系即便算到300nm,nstates也得非常多,好几百个原子的大分子体系要想得到完整的光谱甚至要算好几百乃至上千个态。如果对于你的体系,nstates设几就足够满足你的研究目的那就没必要设好几十。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-15 02:13 , Processed in 0.169744 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list