高斯计算UV-Vis，设置nstates求助

miemiemie

老师您好
高斯TDDFT计算VU-Vis吸收光谱时，第一个分子设置nstates从20到50，但S1激发态（ 536.06 nm）没有出现吸收峰，请问这种情况是否正确。是否需要继续增大nstates？

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

wah

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

wah

好的感谢

wzkchem5

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

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

KSeGaSn

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

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

wah

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

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

KSeGaSn

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

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

wah

请问老师 我计算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

1773525076

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

miemiemie

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

谢谢老师

miemiemie

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

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

非常感谢

乐平

你这不是有 536.06 nm 的吸收峰吗……

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



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

回到计算的结果，你可以在 UV-vis 光谱上点鼠标右键——> Properties



然后在弹出的对话框里减小光谱的半高全宽的数值，也就是我前面说的提高光谱分辨率，如下图所示


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

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

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

sobereva

如果你是指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）

wzkchem5

提问方式不符合规范，按http://sobereva.com/620的要求仔细修改帖子内容，明确你问题的前因后果，我们再回答