计算化学公社

标题: 激发态计算问题 [打印本页]

作者
Author: 风飞 时间: 2016-12-23 12:24
标题: 激发态计算问题
老师：您好！我想算激发态，我的体系是As和Cu的一个-3价的阴离子体系，有文献说Cu的簇化物，且具有较短的Cu-Cu键长，可以有发光性质，是发冷光，因此我想计算一下激发态，看我的这种材料是否具有这种性质，看了您的帖子：http://sobereva.com/265，，还存在如下问题，1. 如果我要算激发态，研究其发光性质。根据您的回复帖子：http://bbs.keinsci.com/forum.php ... =%BC%A4%B7%A2%CC%AC，中研究激发态的关键词
   1）nstates=N：表明总共算能量最低的N个激发态的信息，如激发能、跃迁电/磁偶极矩、振子强度、组态系数等。默认为3。（因此我的也设置为3，可以吗？）
   2）root=i：选择感兴趣的态。如果当前任务是与某个激发态有关的，比如几何优化、振动分析、获得偶极矩等，则表明是对第i个激发态做的。只是计算激发态信息则不需要管此设定。默认             为1。（对于我这个发光，请问感兴趣的态指的是？）
   3）singlet：要求计算的N个激发态都为单重态，此为默认。
   4）triplet：要求计算的N个激发态都为三重态。
2.老师对于我的这样计算是属于价层激发吗？
3，我之前的体系是在PBE1PBE/Aug-cc-pvtz（-pp）条件下优化的，并且加了GD3. 现在要算激发态，如果是价层激发，使用的方法：若不使用耗时的wB97XD，采用b3ly进行计算，差别会很大吗？
4. 溶剂化效应就采用老师您推荐的PCM模型您觉得考虑不足的有哪些？谢谢！

作者
Author: sobereva 时间: 2016-12-23 13:37
1 (1)可 (2)我不知道你说的发冷光的具体机理是怎么回事，反正肯定得优化激发态，你认为最终从第几个激发态发射的就优化哪个态

2 应该属于

3 泛函不同结果差异很大。二者耗时根本差不了太多，都在一个数量级，反倒是你把弥散去了，或者降到某个月份基组耗时能巨幅降低

4 PCM就够

作者
Author: 冰释之川 时间: 2016-12-23 13:47
本帖最后由冰释之川于 2016-12-23 13:50 编辑

1. nstates=N这一设置看你研究的是第几个激发态了，如果是low-lying激发态，设置成5-10没问题，如果你要研究的激发态（用root来指定）较高，则相应的nstates也要调大。有一点规则是：nstates设置的大小必须比root指定的要大，最好大5-10个态以上。感兴趣的态就是指你对第i个激发态感兴趣，举个例子，我做完TD以后，我想要继续研究第二激发态的偶极矩大小，则需要设置root=2 来申明我接下去的计算是针对第2个激发态。

2&3.你可以查阅相关文献；如果是价层激发，B3LYP没问题，参考《简谈量子化学计算中DFT泛函的选择》http://www.sobereva.com/272

4. PCM等隐式溶剂模型存在的最大问题是，没法描述溶剂和溶质分子之间的强相互作用（例如氢键），参考《谈谈隐式溶剂模型下溶解自由能和体系自由能的计算》http://www.sobereva.com/327

作者
Author: 风飞 时间: 2016-12-23 13:58

冰释之川发表于 2016-12-23 13:47
1. nstates=N这一设置看你研究的是第几个激发态了，如果是low-lying激发态，设置成5-10没问题，如果你要研 ...

恩，好的，谢谢老师！我现在去试试

作者
Author: 风飞 时间: 2016-12-23 14:48
本帖最后由风飞于 2016-12-23 15:01 编辑

sobereva 发表于 2016-12-23 13:37
1 (1)可 (2)我不知道你说的发冷光的具体机理是怎么回事，反正肯定得优化激发态，你认为最终从第几个激发态 ...

好的，谢谢老师! 请问怎样确定是价层激发还是电荷激发呢？价层激发是我猜测的，您帖子中说可以用multiwfn判断电荷转移、里德堡激发，请问怎样判断价层激发呢？是否也可以用multiwfn？，在不知道是什么激发的情况下，我要做激发态，您的建议是？谢谢，因为您说不同形式的激发选择的泛函不一样，

作者
Author: sobereva 时间: 2016-12-24 06:30

风飞发表于 2016-12-23 14:48
好的，谢谢老师! 请问怎样确定是价层激发还是电荷激发呢？价层激发是我猜测的，您帖子中说可以用multiwfn ...

没有“电荷激发”这个词

算出来激发态后，参考里面的图，看空轨道等值面图就知道算不算里德堡激发
图解电子激发的分类
http://sobereva.com/284

作者
Author: 风飞 时间: 2016-12-26 14:25

sobereva 发表于 2016-12-24 06:30
没有“电荷激发”这个词

算出来激发态后，参考里面的图，看空轨道等值面图就知道算不算里德堡激发

恩好的谢谢老师！我马上去看看

作者
Author: 风飞 时间: 2016-12-26 16:46
本帖最后由风飞于 2016-12-26 16:50 编辑

sobereva 发表于 2016-12-24 06:30
没有“电荷激发”这个词

算出来激发态后，参考里面的图，看空轨道等值面图就知道算不算里德堡激发

老师：您好！我现在的输入文件和部分输出文件如下：C:\Users\T\Desktop请问出现了对称性不识别这种情况，对吗？谢谢！

作者
Author: sobereva 时间: 2016-12-26 17:40

风飞发表于 2016-12-26 16:46
老师：您好！我现在的输入文件和部分输出文件如下：请问出现了对称性不识别这种情况，对吗？谢谢！

去掉弥散，加上int=ultrafine再试
数值精度问题容易造成判断不出来

作者
Author: 风飞 时间: 2016-12-26 22:27
本帖最后由风飞于 2016-12-26 22:40 编辑

sobereva 发表于 2016-12-26 17:40
去掉弥散，加上int=ultrafine再试
数值精度问题容易造成判断不出来

老师：您好！这个问题解决了，谢谢！请问去掉弥散的原因是？之前在您的帖子中看到说，对于里德堡激发要确定有充足的弥散，所以我就加了，去掉之后，会不会就会减少里德堡激发呢？我要计算荧光和磷光，我的想法是，想计算激发态，看他是属于什么激发，根据激发的类型在确定泛函，用于荧光和磷光的计算，您觉得这样可以吗？谢谢！

作者
Author: sobereva 时间: 2016-12-26 22:36

风飞发表于 2016-12-26 22:27
老师：您好！这个问题解决了，谢谢！我要计算荧光和磷光，我的想法是，想计算激发态，看他是属于什么激发 ...

可以

作者
Author: 风飞 时间: 2016-12-26 22:43

sobereva 发表于 2016-12-26 22:36
可以

老师：您好！抱歉，我的帖子更新慢了，请问之前去掉弥散的原因是？之前在您的帖子中看到说，对于里德堡激发要确定有充足的弥散，所以我就加了，去掉之后，会不会就会减少里德堡激发呢？谢谢!

作者
Author: sobereva 时间: 2016-12-27 07:38

风飞发表于 2016-12-26 22:43
老师：您好！抱歉，我的帖子更新慢了，请问之前去掉弥散的原因是？之前在您的帖子中看到说，对于里德堡激 ...

弥散容易因为造成积分不精确而导致破坏对称性或对称性难以判断
如果你不在意能否显示不可约表示，或者经过测试没法判断不可约表示不是因为这个造成的，对于怀疑有里德堡激发的时候要带弥散函数。

作者
Author: 风飞 时间: 2016-12-27 09:30

sobereva 发表于 2016-12-27 07:38
弥散容易因为造成积分不精确而导致破坏对称性或对称性难以判断
如果你不在意能否显示不可约表示，或者经 ...

恩好的，谢谢老师！

作者
Author: 风飞 时间: 2016-12-27 11:58

sobereva 发表于 2016-12-27 07:38
弥散容易因为造成积分不精确而导致破坏对称性或对称性难以判断
如果你不在意能否显示不可约表示，或者经 ...

老师：您好！我的图解激发态如下，因此判断为电荷转移激发，您看觉得怎样？对吗？谢谢！

作者
Author: 风飞 时间: 2016-12-27 12:31

sobereva 发表于 2016-12-27 07:38
弥散容易因为造成积分不精确而导致破坏对称性或对称性难以判断
如果你不在意能否显示不可约表示，或者经 ...

老师：您好！我现在开始算荧光，关键词如下，根据您的帖子http://sobereva.com/314，中您的TD用的都是默认值，而我的是改了nstates=3为nstates=10. 这样是否有影响？此外，您的帖子中说“第二步最后的结构就是S1结构了。每一步优化都会做一次电子激发计算”，那么我的最后结构也是否S1结构，请问对于S2到的结构呢？是保存S1的结构，然后在这个基础上按照计算S1的关键词再次计算吗？谢谢！前一条回复是我该结构计算的垂直激发能，根据我之前的垂直激发能的计算，请问您觉得我这里计算的S1合适吗？谢谢！

作者
Author: sobereva 时间: 2016-12-27 16:46

风飞发表于 2016-12-27 12:31
老师：您好！我现在开始算荧光，关键词如下，根据您的帖子http://sobereva.com/314，中您的TD用的都是默 ...

什么叫“有影响”？10的话显然多算出7个态来，根据你的目的来设
默认root=1，显然最后优化出来S1结构。要S2结构就root=2，用基态结构还是S1结构当初始结构无所谓

作者
Author: sobereva 时间: 2016-12-27 16:47

风飞发表于 2016-12-27 11:58
老师：您好！我的图解激发态如下，因此判断为电荷转移激发，您看觉得怎样？对吗？谢谢！

对

作者
Author: 风飞 时间: 2016-12-27 20:17

sobereva 发表于 2016-12-27 16:46
什么叫“有影响”？10的话显然多算出7个态来，根据你的目的来设
默认root=1，显然最后优化出来S1结构。 ...

明白了，谢谢老师！请问一般算荧光和磷光一般需要算到Sn,n等于多少合适呢？

作者
Author: sobereva 时间: 2016-12-28 08:30

风飞发表于 2016-12-27 20:17
明白了，谢谢老师！请问一般算荧光和磷光一般需要算到Sn,n等于多少合适呢？

荧光默认的nstates=3就够了
磷光哪来的Sn，磷光是从三重态激发态发射的

作者
Author: 风飞 时间: 2016-12-28 09:45

sobereva 发表于 2016-12-28 08:30
荧光默认的nstates=3就够了
磷光哪来的Sn，磷光是从三重态激发态发射的

恩，谢谢老师!还有就是：
1)磷光是Tn,这个n一般要算到n等于几呢？
2)nstates=3是否也是足够的？
3)对于磷光，算T2，是否也是写root=2？算T3，就写root=3. 以此类推？
4)荧光，算S2，写root=2.算S3，写root=3,以此类推？
谢谢！

作者
Author: sobereva 时间: 2016-12-28 13:44

风飞发表于 2016-12-28 09:45
恩，谢谢老师!还有就是：
1)磷光是Tn,这个n一般要算到n等于几呢？
2)nstates=3是否也是足够的？

用默认的

再仔细阅读和揣摩我的博文，里面该说的都说了

作者
Author: 风飞 时间: 2016-12-30 09:31

sobereva 发表于 2016-12-27 07:38
弥散容易因为造成积分不精确而导致破坏对称性或对称性难以判断
如果你不在意能否显示不可约表示，或者经 ...

老师：您好！这次激发态计算又出现了不能识别不可约表示，输入文件和部分输出文件如下，谢谢！

作者
Author: sobereva 时间: 2016-12-30 10:14

风飞发表于 2016-12-30 09:31
老师：您好！这次激发态计算又出现了不能识别不可约表示，输入文件和部分输出文件如下，谢谢！

自行根据轨道的不可约表示，通过直积手动计算电子态不可约表示
如果轨道不可约表示都没判断出来，通过图形根据特征表标自行判断

作者
Author: 虎王 时间: 2017-1-16 10:16

sobereva 发表于 2016-12-23 13:37
1 (1)可 (2)我不知道你说的发冷光的具体机理是怎么回事，反正肯定得优化激发态，你认为最终从第几个激发态 ...

发冷光就是luminescence，即普通的磷光或荧光，应该不包括红外光，因为红外光可以加热。OLED发的都是冷光，所以开久了也不会变热。白炽灯则明显发热光。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)