计算化学公社

标题: 关于计算紫外光谱 [打印本页]

作者
Author:
lastzealot    时间: 2015-6-7 21:49
标题: 关于计算紫外光谱
以前都是优化基态结构,从没计算过激发态。最近需要计算一种络合物的紫外光谱,方法是不是:
1  优化基态结构
2  # tddft=(root, nstate) opt freq
以上两步对吗?下一步要怎么做呢?
谢谢sob老师,小弟菜鸟请多指教。
作者
Author:
小范范1989    时间: 2015-6-7 23:22
本帖最后由 小范范1989 于 2015-6-8 06:53 编辑



作者
Author:
sobereva    时间: 2015-6-8 00:27
2不对。

2的时候就写比如B3LYP/6-31G* TD(nstates=20)就行了。nstates取多大看情况,越大的体系往往要设得越大,否则计算出来的光谱范围够不到低波长范围。

PS:别信小木虫上的。
作者
Author:
lastzealot    时间: 2015-6-8 09:10
sobereva 发表于 2015-6-8 00:27
2不对。

2的时候就写比如B3LYP/6-31G* TD(nstates=20)就行了。nstates取多大看情况,越大的体系往往要设 ...

谢谢老师的指点
请问第三步中的优化激发态结构如何来进行呢?需要优化第几激发态呢?
多重度是否需要更改呢?
查看哪个输出文件可以看到光谱呢?
作者
Author:
sobereva    时间: 2015-6-8 14:05
看此贴
使用Multiwfn绘制红外、拉曼、UV-Vis、ECD和VCD光谱图
http://sobereva.com/224

我不知道你说的第三步指什么
作者
Author:
万里云    时间: 2015-6-8 15:39
TDDFT基于Mark-Casida方程计算激发能,最终转化为一个赝哈密顿量特征值问题。

求解这个赝哈密顿量特征值,一般用占据态和未占据态波函数直积作为基组,求对应矩阵本征值。理论上最低的几个本征值应该是定值,实际上与矩阵大小,也就是构造基组时所用占据态和未占据态波函数个数乘积有关,需要选定某个激发态做“激发能-基组大小”收敛测试。有的程序,比如octopus,可以直接指定用多少个占据态和多少个未占据态;如果不能直接占据态和未占据态数目,也可以指定激发态总数,因为激发态总数与基组大小相等。

物理上在求解类似的Bethe-Salpeter方程时也有类似的问题。一个常用的方法是:假如要计算0-3eV内激发态,则从最高占据态能量往上数3eV,从最低未占据态能量往下数3eV,统计在此区间内的态数量作为收敛测试起点。比如在此区间内有5个占据态,6个未占据态,则至少考虑5*6=30个激发态。
作者
Author:
lastzealot    时间: 2015-6-9 00:07
万里云 发表于 2015-6-8 15:39
TDDFT基于Mark-Casida方程计算激发能,最终转化为一个赝哈密顿量特征值问题。

求解这个赝哈密顿量特征值 ...

谢谢
作者
Author:
虎王    时间: 2015-6-9 12:31
sobereva 发表于 2015-6-8 00:27
2不对。

2的时候就写比如B3LYP/6-31G* TD(nstates=20)就行了。nstates取多大看情况,越大的体系往往要设 ...

不需要写singlet吗?
作者
Author:
sobereva    时间: 2015-6-9 14:23
虎王 发表于 2015-6-9 12:31
不需要写singlet吗?

默认就是
作者
Author:
xida    时间: 2017-2-15 09:13
sobereva 发表于 2015-6-8 14:05
看此贴
使用Multiwfn绘制红外、拉曼、UV-Vis、ECD和VCD光谱图
http://sobereva.blog.163.com/blog/static ...

sobereva老师,能不能把使用Multiwfn绘制红外、拉曼、UV-Vis、ECD和VCD光谱图的链接重新发一下,你这个打不开
作者
Author:
sobereva    时间: 2017-2-15 18:01
xida 发表于 2017-2-15 09:13
sobereva老师,能不能把使用Multiwfn绘制红外、拉曼、UV-Vis、ECD和VCD光谱图的链接重新发一下,你这个打 ...

链接已更新

如果碰见博文链接访问不了的情况,在sobereva.com中的索引里面直接搜标题即可




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3