计算化学公社
标题: 求助：如何解决高斯计算振动分辨谱时输出结果含有WARNING信息，最终光谱只显示50% [打印本页]
作者Author: SunriseOkk 时间: 2024-7-31 19:51
标题: 求助：如何解决高斯计算振动分辨谱时输出结果含有WARNING信息，最终光谱只显示50%
求助各位老师，我在计算阴离子基态到中性分子基态的振动分辨光谱时，输出的结果中一直提示 WARNING: About 50% of the calculated intensity is shown in the final spectrum.
学习了sob老师之前的关于《振动分辨的电子光谱的计算》的帖子，按照里面的参数设置做了一些修改，但是问题依旧存在，想请教各位老师这个问题是什么原因导致的，该如何解决？
下面是我的计算输入文件的参数设置，恳请各位老师指正。
第一步计算：使用了中性分子的优化结构
%chk=PtNH3Cl3-a-0C-EXC.chk
# b3lyp/genecp em=gd3bj opt freq=saveNM
Title Card Required
0 2
Pt 0.00000200 0.09008700 0.04439400
Cl -2.24452700 0.21002900 -0.17510000
Cl -0.00000600 -2.18650800 0.06223500
Cl 2.24453100 0.21002900 -0.17510000
H -0.83213600 2.56477800 -0.31475700
H -0.00104200 2.50111400 1.10473100
H 0.83312400 2.56469700 -0.31299200
N -0.00000700 2.19603700 0.13652600
第二步计算：使用了阴离子的优化结构
%chk=PtNH3Cl3-FCF.chk
# opt freq=(fc,readfcht) b3lyp/genecp em=gd3bj
Title Card Required
-1 1
Pt 0.00003700 0.08565700 0.00000100
Cl -2.34319800 0.27137300 -0.00000200
Cl -0.01830100 -2.23180700 0.00000000
Cl 2.34724400 0.24116300 -0.00000200
H -0.45859100 2.50232600 -0.82210800
H -0.45853200 2.50232800 0.82214400
H 0.99625000 2.45068700 -0.00003200
N 0.02290300 2.15585900 0.00000100
.........(基组)
MaxOvr=500 SpecHwHm=30 Advanced=(ForceFCCalc) FORCEPRTSPECTRUM !PreScreening=(MaxC1=10,MaxC2=10,MaxInt=1000) RedDim=ClearLowFreq=10
PtNH3Cl3-a-0C-EXC.chk
输出结果如图所示。
作者Author: wzkchem5 时间: 2024-7-31 20:25
这是高斯计算振动分辨光谱的方法的固有缺陷，改用ORCA一般能解决。参见http://bbs.keinsci.com/forum.php ... 7&fromuid=21811
作者Author: SunriseOkk 时间: 2024-7-31 20:31
嗯嗯好的，多谢指导，我试一下
作者Author: SunriseOkk 时间: 2024-7-31 21:41
老师您好，我想改用ORCA来计算振动分辨光谱，但是查阅了orca_manual之后暂时没找到相应的案例，想请老师帮忙提示一下在哪一章节或者有没有其他相关的帖子做参考？
作者Author: wzkchem5 时间: 2024-7-31 22:59
搜vibronic，很容易搜到
作者Author: SunriseOkk 时间: 2024-8-1 11:39
老师，我按照手册里做了部分计算，在最后计算光谱时结果报错Error on ESD: Different number of vibrational modes between initial and final states!Error in calculating J and K。我是想计算阴离子基态到中性基态的振动谱，不知道输入文件里写的对不对？
作者Author: wzkchem5 时间: 2024-8-1 16:16
仔细检查结构和Hessian有没有问题，有没有多画或漏画原子
作者Author: SunriseOkk 时间: 2024-8-1 16:29
好的，多谢老师。我是直接用了高斯优化的结构，然后用orca只算了freq，现在重新用orca算opt+freq，感觉可能是这里出的问题
作者Author: SunriseOkk 时间: 2024-8-1 17:14
老师，我重新优化了结构，检查了Hess文件，里面modes数量也都是一样的，但是还是出现同样的报错
作者Author: wzkchem5 时间: 2024-8-1 17:20
可能有一个问题是不应该加doht true，因为orca没有实现电离过程的HT效应的计算，只能对电子数守恒的吸收光谱算HT效应。
此外检查ESD计算输入文件的结构是否和两个Hessian里的结构相符。注意不要光检查mode的数量，也要检查各个结构的原子数目、种类、顺序是否一致。
作者Author: SunriseOkk 时间: 2024-8-1 18:13
老师，我重新修改后这个问题解决了，但是输出结果中在显示Calculating energy of the ES at its geometry:后，直接报错ORCA finished by error termination in ORCA_ESD，现在ESD计算输入文件的结构与两个Hessian里的结构是一致的。不知道这样还有没有办法解决...麻烦老师了
作者Author: wzkchem5 时间: 2024-8-1 18:19
K*K不应该这么接近0，建议检查你的阴离子和中性态的Hessian是不是用成同一个文件了
作者Author: SunriseOkk 时间: 2024-8-1 18:47
老师，我现在有点不太清楚该如何处理或者计算中性态的hessian文件，通过几次尝试发现，对中性态只做freq计算，结构与阴离子保持一致，就会出现K*K为0的情况。但是如果对中性态做opt+freq的话，结构变化，就会出现最开始的报错Error on ESD: Different number of vibrational modes between initial and final states
作者Author: wzkchem5 时间: 2024-8-1 19:50
我才发现你用的是ORCA 5.0.4，可以试试ORCA 6.0.0行不行。如果还是不行，可以把输入文件以及两个.hess文件发给我，我们测试一下
作者Author: SunriseOkk 时间: 2024-8-1 20:37
老师，我这边计算也是用的6.0版本，不过我是参考的5.0版本手册写的计算文件...
附件里面分别是计算文件和两个Hess文件，其中中性态的文件也是计算了opt+freq之后的，请老师帮忙检查指导。
感谢老师！
作者Author: SunriseOkk 时间: 2024-8-2 15:22
老师，我今天又算了一个类似的体系，K*K值5.4左右，最后还是出现同样的问题，提示ORCA finished by error termination in ORCA ESD。是不是ORCA不太适合计算电离过程的振动分辨光谱？
作者Author: wzkchem5 时间: 2024-8-2 17:30
有可能是bug，已经和相关人员反映了。我个人的感觉是这个问题和电离无关，可能和有的键角接近180度有关
| 欢迎光临 计算化学公社 (http://bbs.keinsci.com/) |
Powered by Discuz! X3.3 |