计算化学公社

标题: 高斯计算UV-Vis中，nstates设置求助。 [打印本页]

作者
Author: Ethan_WY 时间: 2018-9-28 09:10
标题: 高斯计算UV-Vis中，nstates设置求助。
大家好，老师们好：
我想请问一下，看了”思想家公社“计算UV-Vis那篇文章后，有个疑问，现在我有一个化合物，分子式如图所示，感觉属于大共轭系统，我将nstates设置的20，不知合理不合理？

想要验证，是不是需要我再设置不同的nstates，然后对比结果？只是计算一个化合物时间太长了。感觉这样比较下来太浪费时间，请老师们指点一二。

我的输入文件大概是：
1：优化结构；# opt b3lyp/6-31g(d,p) scrf=(solvent=dichloromethane) geom=connectivity
2：计算激发；# b3lyp/6-31g(d,p) TD(nstates=20) scrf=(solvent=dichloromethane)
3:取1得到的log文件，进行S1优化，得到荧光；# b3lyp/6-31g(d,p) TD scrf=(solvent=dichloromethane) opt

作者
Author: 让你变成回忆 时间: 2018-9-28 09:13
（1）大的共轭体系，利用B3LYP泛函是否合适？是否应该考虑其他的泛函？我觉得这个问题比nstates的问题更重要；这是一个明显的D-A-D体系，因此泛函的问题肯定更加重要；
（2）1,3都没有问题，2的话要看你想要研究多少个态，而nstates一般比你想研究的态的数目多一点就行（比如多算5个）；

作者
Author: sobereva 时间: 2018-9-28 09:28
给氢加极化完全多余，浪费时间
这体系没多大，有个像样的机子算100个态都不费劲。20个态足够覆盖可见光区域，如果想算得完整点，一次算50个态。没必要反复测试，一次算50个态的耗时比你依次测试20、30、40低得多。

作者
Author: Ethan_WY 时间: 2018-9-28 09:29

让你变成回忆发表于 2018-9-28 09:13
（1）大的共轭体系，利用B3LYP泛函是否合适？是否应该考虑其他的泛函？我觉得这个问题比nstates的问题更重 ...

你好：
我看了下sob老师的贴子，sob老师有这样一些选择建议：
计算垂直局域价层激发能：PBE0（误差约在0.2~0.3eV，平均好于B3LYP一丝）
计算大共轭体系激发能、电荷转移、高阶价层、里德堡垂直激发：wB97XD、CAM-B3LYP
计算各类激发能，只用一个泛函的情况：wB97XD、M06-2X
有机分子结构优化：PBE0（弱相互作用明显时应加DFT-D3校正）、M06-2X
计算HOMO-LUMO gap：HSE06、B3PW91（整体来说这个好，HSE06的杂化参数有点依赖于体系）

有几个疑问：
1.我优化分子结构，选择PBE0（弱相互作用明显时应加DFT-D3校正）、M06-2X；
2.计算激发选择wB97XD、CAM-B3LYP；
3.如果我还想计算HOMO-LUMO gap，是不是又要换HSE06、B3PW91重新计算？

哎，真是难为住我了。还请老师指点指点。

作者
Author: Ethan_WY 时间: 2018-9-28 09:32

sobereva 发表于 2018-9-28 09:28
给氢加极化完全多余，浪费时间
这体系没多大，有个像样的机子算100个态都不费劲。20个态足够覆盖可见光区 ...

恩恩，谢谢老师指点。
那这个B3LYP泛函是否合适呢？我是看文献上，和我的分子差不多，文献用的也是b3lyp。

作者
Author: sobereva 时间: 2018-9-28 09:36

Ethan_WY 发表于 2018-9-28 09:29
你好：
我看了下sob老师的贴子，sob老师有这样一些选择建议：
计算垂直局域价层激发能：PBE0（误 ...

都说得很明确，一点也不难为
PBE0优化
CAM-B3LYP做TDDFT
gap还用PBE0，优化完了直接就有了
把PBE0换成更常用的B3LYP也没问题

作者
Author: Ethan_WY 时间: 2018-9-28 10:36

sobereva 发表于 2018-9-28 09:36
都说得很明确，一点也不难为
PBE0优化
CAM-B3LYP做TDDFT

好的，感谢老师指点

作者
Author: Ethan_WY 时间: 2018-9-30 17:26

sobereva 发表于 2018-9-28 09:36
都说得很明确，一点也不难为
PBE0优化
CAM-B3LYP做TDDFT

老师，我想了想，最终这样输入，来计算我的紫外和荧光，化合物还是上面发帖时带的图，您看这个选择怎么样？

1.优化结构,得到HOMO-LUMO gap：# opt B3LYP/6-31g(d) scrf=(solvent=dichloromethane) geom=connectivity
2.计算50个激发态，得到UV-Vis图：# CAM-B3LYP/6-31g(d) TD(nstates=50) scrf=(solvent=dichloromethane)
3:取1优化的结果log文件，优化S1，计算荧光：# opt td B3LYP/6-31g(d) scrf=(solvent=dichloromethane)

长假来临，祝老师国庆节快乐！！

作者
Author: sobereva 时间: 2018-9-30 20:50

Ethan_WY 发表于 2018-9-30 17:26
老师，我想了想，最终这样输入，来计算我的紫外和荧光，化合物还是上面发帖时带的图，您看这个选择怎么样 ...

激发态相关任务泛函应统一，优化S1应当也用CAM-B3LYP

作者
Author: Ethan_WY 时间: 2018-10-1 10:50

sobereva 发表于 2018-9-30 20:50
激发态相关任务泛函应统一，优化S1应当也用CAM-B3LYP

好的，我知道了，感谢老师！！

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