基态与激发态几何结构变化明显时振动分辨的电子光谱计算报错求助

wangxm · 发表于 Post on 2021-8-19 11:20:49

请教各位老师：基态与激发态几何结构变化明显，计算出的光谱输出文件报错：The Franck-Condon factor corresponding to the overlap integral between both vibrational ground states is too small:该如何解决？通过看帖子与文献《Vibrationally-resolved electronic spectra in GAUSSIAN 09 》，如果强行计算，结果只能用于测试。那么请问，针对这种情况，还能如何计算其振动分辨的电子光谱呢？

sobereva · 发表于 Post on 2021-8-19 11:28:13

如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“振动分辨的电子光谱计算”改了，以后务必注意，下次将删帖+扣分处理。

wangxm · 发表于 Post on 2021-8-19 11:36:52

sobereva 发表于 2021-8-19 11:28
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此 ...

不好意思，我下次注意

wzkchem5 · 发表于 Post on 2021-8-19 15:58:56

可以用ORCA的AHAS、VG等方法计算，相对不受基态激发态结构差别影响

wangxm · 发表于 Post on 2021-8-26 15:28:41

wzkchem5 发表于 2021-8-19 15:58
可以用ORCA的AHAS、VG等方法计算，相对不受基态激发态结构差别影响

好的，谢谢。此外，我还想请教一下，看了有关高斯计算振动分辨光谱的文献，文献提及适用条件是：非绝热耦合可以忽略以及谐波近似可靠，请问：怎么判断是否满足这两个条件，对于这种两个电子态几何结构相差较大的情况不能用FC方法是因为不满足条件吗？

wzkchem5 · 发表于 Post on 2021-8-26 16:51:14

wangxm 发表于 2021-8-26 08:28
好的，谢谢。此外，我还想请教一下，看了有关高斯计算振动分辨光谱的文献，文献提及适用条件是：非绝热耦 ...

非绝热耦合可以忽略：等价于电子波函数对结构变化不是太敏感。一般只要基态势能面和激发态势能面之间距离比较远（比如处处大于1eV，不过这不是一个严格的标准），就可以满足
谐振近似可靠（注意不是谐波是谐振）：意思是基态势能面和激发态势能面都可以近似成抛物面。大体上讲，只要基态和激发态都没有特别低的振动频率（<100 cm^-1），且无论是基态还是激发态，和附近的电子态之间的能量差都比较大，且基态和激发态的平衡位置差别不大，就可以认为满足。但是同样地，这不是一个严格的标准
注意基态和激发态的平衡位置差别大，并不直接违反了谐振近似，它只是增加了谐振近似不成立的概率。所以高斯发现基态和激发态的平衡位置差别大到一定程度，就告诉你结果不可靠，而不去直接检查谐振近似是否成立

wangxm · 发表于 Post on 2021-8-26 17:36:55

wzkchem5 发表于 2021-8-26 16:51
非绝热耦合可以忽略：等价于电子波函数对结构变化不是太敏感。一般只要基态势能面和激发态势能面之间距离 ...

首先，非常感谢您的回答。我还有个疑问：对于此类情况（由于结构变化太大造成的FC重叠积分太小），根据文献《Vibrationally-resolved electronic spectra in GAUSSIAN 09》中对于非谐振模型振动分辨光谱的计算例子，可否通过考虑非谐效应处理此类问题？当我在原有体系下进行非谐振模型振动光谱计算时，仍然存在上述问题，且重叠积分比之前谐振近似计算更小，这可以得出什么结论呢，是因为不满足第一个条件吗？

wzkchem5 · 发表于 Post on 2021-8-26 18:13:48

wangxm 发表于 2021-8-26 10:36
首先，非常感谢您的回答。我还有个疑问：对于此类情况（由于结构变化太大造成的FC重叠积分太小），根据文 ...

考虑非谐性效应的方法都是近似的，很多考虑非谐性效应的方法仍然需要基态和激发态结构不能差得太远，但是有可能对结构差别的要求比谐振近似要松一些。理论上如果能精确考虑非谐性效应的话，基态和激发态结构差多少都没问题

wangxm · 发表于 Post on 2021-8-26 19:53:38

wzkchem5 发表于 2021-8-26 18:13
考虑非谐性效应的方法都是近似的，很多考虑非谐性效应的方法仍然需要基态和激发态结构不能差得太远，但是 ...

好的，十分感谢

wangxm · 发表于 Post on 2021-8-31 15:52:35

wangxm 发表于 2021-8-26 15:28
好的，谢谢。此外，我还想请教一下，看了有关高斯计算振动分辨光谱的文献，文献提及适用条件是：非绝热耦 ...

您好，听了您的意见，最近有尝试使用ORCA来计算，有如下几个问题想要请教您：
1.VG方法中需要基态的平衡结构(.xyz文件）及Heaaian(.hess文件），这两个文件必须在ORCA中做OPT+FREQ计算才能得到吗?
2.我在高斯中进行了opt+freq的计算，得到没有虚频的几何构型，取该几何构型放入ORCA中进行OPT+FREQ计算（输入文件如下），结果出现两个虚频，该如何处理，不处理虚频能否直接调用ESD进行光谱计算，其结果是否合理？
! CAM-B3LYP 6-31+G** OPT FREQ
%maxcore 1000
%pal nprocs 24 end
* xyz 0 1
分子坐标
*
3.能否实现由高斯计算opt+freq，ORCA只用来做光谱计算呢？

zjxitcc · 发表于 Post on 2021-8-31 15:58:58

wzkchem5 发表于 2021-8-26 18:13
考虑非谐性效应的方法都是近似的，很多考虑非谐性效应的方法仍然需要基态和激发态结构不能差得太远，但是 ...

有啥方法能（比较）精确的考虑非谐性效应么？vibrational CC?

心向暖阳 · 发表于 Post on 2021-8-31 16:18:14

wangxm 发表于 2021-8-31 15:52
您好，听了您的意见，最近有尝试使用ORCA来计算，有如下几个问题想要请教您：
1.VG方法中需要基态的平衡 ...

1.VG近似只需要初始态的频率，然后用初态的结构计算末态的force就可以了
2.这个是很正常的，软件不一样，没虚频才不正常
3.可以吧，看哪个高手有空写一下接口

wzkchem5 · 发表于 Post on 2021-8-31 17:04:51

zjxitcc 发表于 2021-8-31 08:58
有啥方法能（比较）精确的考虑非谐性效应么？vibrational CC?

嗯，不过非谐性效应如果真的特别大的话，除了动力学以外没别的办法

心向暖阳 · 发表于 Post on 2021-8-31 17:24:32

wzkchem5 发表于 2021-8-31 17:04
嗯，不过非谐性效应如果真的特别大的话，除了动力学以外没别的办法

大佬能展开讲讲么

wzkchem5 · 发表于 Post on 2021-8-31 17:46:21

心向暖阳发表于 2021-8-31 10:24
大佬能展开讲讲么

举个极端的例子，假如你要算H2的磷光光谱。H2的T1态几乎是解离的，但没有完全解离，因为两个H有一点色散作用，平衡键长大概有4个埃。可以想见，在这里T1和S0几乎是简并的。你但凡用基于S0和T1的平衡结构的方法，都会有一个问题：这个方法必须正确预测出T1平衡结构处S0和T1几乎简并，但是还不能预测出S0比T1高。而如果只有平衡位置的高阶导数，这几乎必然是做不到的。尤其假如你用的是DFT，那更是不可能的，因为broken symmetry DFT在RKS instability point处不解析，所以你不可能只用S0平衡位置附近的S0势能面就预测出T1平衡位置附近的S0势能面。因此这时只能直接对T1平衡位置附近的S0势能面进行采样，而不能完全依赖两个势能面各自的平衡位置附近的信息，否则哪怕你用多高阶的势能面导数，做vibrational full CI都无济于事。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[综合交流] 基态与激发态几何结构变化明显时振动分辨的电子光谱计算报错求助