计算化学公社

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

作者
Author: miemiemie 时间: 2024-2-26 20:33
标题: 高斯计算UV-Vis，设置nstates求助
本帖最后由 miemiemie 于 2024-2-26 22:55 编辑

老师您好
高斯TDDFT计算VU-Vis吸收光谱时，第一个分子设置nstates从20到50，但S1激发态（ 536.06 nm）没有出现吸收峰，请问这种情况是否正确。是否需要继续增大nstates？
(, 下载次数 Times of downloads: 9)
Excited State 1:    Singlet-A    2.3129 eV  536.06 nm  f=1.2176  <S**2>=0.000
   142 ->143       0.69758
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-DFT) =  -2096.53653509
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State 2:    Singlet-A    2.8426 eV  436.16 nm  f=1.4960  <S**2>=0.000
   141 ->143       0.45306
   142 ->143       0.10860
   142 ->144       0.52719

Excited State 3:    Singlet-A    3.0158 eV  411.12 nm  f=0.1644  <S**2>=0.000
   141 ->143       0.53369
   142 ->144       -0.45722

Excited State 4:    Singlet-A    3.4487 eV  359.51 nm  f=0.0183  <S**2>=0.000
   139 ->143       -0.10184
   140 ->143       0.68560

Excited State 5:    Singlet-A    3.5304 eV  351.19 nm  f=0.0049  <S**2>=0.000
   138 ->143       0.10973
   141 ->144       0.67812
   142 ->145       0.10289

Excited State 6:    Singlet-A    3.6756 eV  337.31 nm  f=0.0295  <S**2>=0.000
   139 ->143       0.64415
   139 ->144       -0.22658
   142 ->145       0.10490

作者
Author: wzkchem5 时间: 2024-2-26 21:40
提问方式不符合规范，按http://sobereva.com/620的要求仔细修改帖子内容，明确你问题的前因后果，我们再回答

作者
Author: sobereva 时间: 2024-2-27 04:12
如果你是指500多nm的那个电子激发没有“独立”的峰，仔细看下文搞清楚理论光谱产生的原理，自然知道怎么回事。还可以让Multiwfn把各个电子激发对光谱的贡献单独给出，能理解得更清楚
使用Multiwfn计算特定方向的UV-Vis吸收光谱
http://sobereva.com/648（http://bbs.keinsci.com/thread-30873-1-1.html）

盲目设大nstates没有任何意义，若只关注S1，用默认的nstates=3通常就足够了，看下文了解nstates的正确设置
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314（http://bbs.keinsci.com/thread-2413-1-1.html）

作者
Author: 乐平 时间: 2024-2-27 10:50
你这不是有 536.06 nm 的吸收峰吗……

而且 536.06 nm 是第一个激发态，你增加 NStates 干嘛呢…… 而且你也算了 NStates=20 和 NStates=50，你会发现增大数值只会计算更多更高的激发态，前面的激发态是一模一样的。无论你再增大 NStates，相同体系相同计算参数得到的第一个激发态都是一样的。

(, 下载次数 Times of downloads: 11)

图中的光谱是总包的轮廓。如果你有一点点实验的知识就会知道，总包的光谱轮廓实际上是多个光谱峰叠加的结果。从仪器的角度来讲，可能是分辨率不高，或者是溶剂的影响等等诸多因素。

回到计算的结果，你可以在 UV-vis 光谱上点鼠标右键——> Properties
(, 下载次数 Times of downloads: 9)

然后在弹出的对话框里减小光谱的半高全宽的数值，也就是我前面说的提高光谱分辨率，如下图所示
(, 下载次数 Times of downloads: 11)

减小就能看到原本的大轮廓分成两个小峰。

如果你有一系列光谱计算，那么其他的 UV-vis 计算结果要将光谱的半高全宽的数值保持一致。

当然，这是粗略的方法。如果要分析谱峰的性质，看 3楼社长的帖子 http://bbs.keinsci.com/forum.php ... 53&fromuid=1532

作者
Author: miemiemie 时间: 2024-2-27 14:23

乐平发表于 2024-2-27 10:50
你这不是有 536.06 nm 的吸收峰吗……

而且 536.06 nm 是第一个激发态，你增加 NStates 干嘛呢…… 而 ...

非常感谢

作者
Author: miemiemie 时间: 2024-2-27 14:24

sobereva 发表于 2024-2-27 04:12
如果你是指500多nm的那个电子激发没有“独立”的峰，仔细看下文搞清楚理论光谱产生的原理，自然知道怎么回 ...

谢谢老师

作者
Author: 1773525076 时间: 2024-2-27 19:45
Excited State 1: Singlet-A 2.3129 eV 536.06 nm f=1.2176 <S**2>=0.000 你自己给出的信息里面都有这个吸收峰了而且就算没有你现在设置的nstates已经囊括到这个波长范围了改变nstates也没有意义了

作者
Author: wah 时间: 2024-3-27 16:04
请问老师我计算UV-Vis时，高斯没有图生成是为什么
我的部分输入文件：
%nprocshared=64
%chk=MgC60.chk
#  6-31g  m062x  TD

Title Card Required

0 1
C                -1.13358700 0.14493100 -0.00431800
C                0.32464900 0.14496800 -0.00438900
C                1.01966000 1.34905600 -0.00455900
C                0.29064900 2.61205500 -0.00435300
C                -1.09971100 2.61201900 -0.00428800
C                -1.82867600 1.34898200 -0.00437100
C                -1.58429300 -0.95710400 -0.84625800
C                -0.40446600 -1.63817200 -1.36661500
C                0.77536200 -0.95703200 -0.84636300
C                1.90021500 -0.80377400 -1.64895400
C                2.19930500 1.50986800 -0.84654800
C                1.01974400 3.55333700 -0.84635700
C                0.32451800 4.45071800 -1.64904200
C                -1.13377300 4.45068600 -1.64898000
C                -1.82890600 3.55328000 -0.84623100
C                -3.00855500 2.87198700 -1.36647900
C                -3.00837000 1.50974000 -0.84630800
C                -3.43837600 0.45893900 -1.64879100
C                -2.70920400 -0.80389600 -1.64876900
C                -0.40449100 -2.13441700 -2.66534200
C                -1.58431100 -1.97361800 -3.50719100
C                -2.70913700 -1.32428200 -3.01105400

作者
Author: KSeGaSn 时间: 2024-3-27 17:42

wah 发表于 2024-3-27 16:04
请问老师我计算UV-Vis时，高斯没有图生成是为什么
我的部分输入文件：
%nprocshared=64

在sob老师看到之前把你的计算水平至少提升到6-31G*……当下的计算恐怕毫无意义。
另外这个图上标出来了位置，之所以你没看到曲线是因为振子强度=0 所以y轴高度是0。

作者
Author: wah 时间: 2024-3-27 18:28

KSeGaSn 发表于 2024-3-27 17:42
在sob老师看到之前把你的计算水平至少提升到6-31G*……当下的计算恐怕毫无意义。
另外这个图上标出来了 ...

好的另外我想问一下振子强度为啥为0 怎样修改

作者
Author: KSeGaSn 时间: 2024-3-27 22:26

wah 发表于 2024-3-27 18:28
好的另外我想问一下振子强度为啥为0 怎样修改

显然这不是一个"修改"的问题，振子强度代表了辐射跃迁的性质，振子强度为0说明这种辐射跃迁的概率很小（乃至为暗态）。分析这些数据的提前是必须在合理的水平上进行了计算比如你现在用6-31G基组进行计算就完全谈不上合理；研究吸收光谱默认TD关键词计算三个态也是远远不够的；可能还有其他结构问题。等到把这些问题都解决了重新进行了修正，再来分析振子强度是否为0这件事也不迟。

作者
Author: wzkchem5 时间: 2024-3-27 22:50

wah 发表于 2024-3-27 11:28
好的另外我想问一下振子强度为啥为0 怎样修改

振子强度为0是正常的，说明你要算的态是暗态。
如果实验确定不是暗态，就说明你的计算不合理，例如分子结构不合理、理论级别不合理、没有考虑Herzberg-Teller效应等。

作者
Author: wah 时间: 2024-3-28 09:07
好的感谢

作者
Author: wah 时间: 2024-3-29 09:08
老师们好我，这是我对C60进行Uv-vis的计算依旧没有峰，以下是我的输入文件 6-31g* nosymm m062x TD(nstates=20) 按理说C60是近紫外区是有吸收的，求各位老师帮助。

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