计算化学公社

 找回密码 Forget password
 注册 Register
Views: 693|回复 Reply: 16
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 求助:如何解决高斯计算振动分辨谱时输出结果含有WARNING信息,最终光谱只显示50%

[复制链接 Copy URL]

10

帖子

0

威望

33

eV
积分
43

Level 2 能力者

求助各位老师,我在计算阴离子基态到中性分子基态的振动分辨光谱时,输出的结果中一直提示 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

输出结果如图所示。





202407311949189901..png (40.41 KB, 下载次数 Times of downloads: 21)

202407311949189901..png

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

2#
发表于 Post on 2024-7-31 20:25:27 | 只看该作者 Only view this author
这是高斯计算振动分辨光谱的方法的固有缺陷,改用ORCA一般能解决。参见http://bbs.keinsci.com/forum.php ... 7&fromuid=21811
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

10

帖子

0

威望

33

eV
积分
43

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2024-7-31 20:31:05 | 只看该作者 Only view this author
wzkchem5 发表于 2024-7-31 20:25
这是高斯计算振动分辨光谱的方法的固有缺陷,改用ORCA一般能解决。参见http://bbs.keinsci.com/forum.php?m ...

嗯嗯好的,多谢指导,我试一下

10

帖子

0

威望

33

eV
积分
43

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2024-7-31 21:41:28 | 只看该作者 Only view this author
wzkchem5 发表于 2024-7-31 20:25
这是高斯计算振动分辨光谱的方法的固有缺陷,改用ORCA一般能解决。参见http://bbs.keinsci.com/forum.php?m ...

老师您好,我想改用ORCA来计算振动分辨光谱,但是查阅了orca_manual之后暂时没找到相应的案例,想请老师帮忙提示一下在哪一章节或者有没有其他相关的帖子做参考?

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

5#
发表于 Post on 2024-7-31 22:59:17 | 只看该作者 Only view this author
SunriseOkk 发表于 2024-7-31 14:41
老师您好,我想改用ORCA来计算振动分辨光谱,但是查阅了orca_manual之后暂时没找到相应的案例,想请老师 ...

搜vibronic,很容易搜到
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

10

帖子

0

威望

33

eV
积分
43

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2024-8-1 11:39:53 | 只看该作者 Only view this author
wzkchem5 发表于 2024-7-31 22:59
搜vibronic,很容易搜到

老师,我按照手册里做了部分计算,在最后计算光谱时结果报错Error on ESD: Different number of vibrational modes between initial and final states!Error in calculating J and K。我是想计算阴离子基态到中性基态的振动谱,不知道输入文件里写的对不对?

202408011139235842..png (35.68 KB, 下载次数 Times of downloads: 8)

202408011139235842..png

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

7#
发表于 Post on 2024-8-1 16:16:09 | 只看该作者 Only view this author
SunriseOkk 发表于 2024-8-1 04:39
老师,我按照手册里做了部分计算,在最后计算光谱时结果报错Error on ESD: Different number of vibratio ...

仔细检查结构和Hessian有没有问题,有没有多画或漏画原子
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

10

帖子

0

威望

33

eV
积分
43

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2024-8-1 16:29:10 | 只看该作者 Only view this author
wzkchem5 发表于 2024-8-1 16:16
仔细检查结构和Hessian有没有问题,有没有多画或漏画原子

好的,多谢老师。我是直接用了高斯优化的结构,然后用orca只算了freq,现在重新用orca算opt+freq,感觉可能是这里出的问题

10

帖子

0

威望

33

eV
积分
43

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2024-8-1 17:14:31 | 只看该作者 Only view this author
wzkchem5 发表于 2024-8-1 16:16
仔细检查结构和Hessian有没有问题,有没有多画或漏画原子

老师,我重新优化了结构,检查了Hess文件,里面modes数量也都是一样的,但是还是出现同样的报错

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

10#
发表于 Post on 2024-8-1 17:20:02 | 只看该作者 Only view this author
SunriseOkk 发表于 2024-8-1 10:14
老师,我重新优化了结构,检查了Hess文件,里面modes数量也都是一样的,但是还是出现同样的报错

可能有一个问题是不应该加doht true,因为orca没有实现电离过程的HT效应的计算,只能对电子数守恒的吸收光谱算HT效应。
此外检查ESD计算输入文件的结构是否和两个Hessian里的结构相符。注意不要光检查mode的数量,也要检查各个结构的原子数目、种类、顺序是否一致。
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

10

帖子

0

威望

33

eV
积分
43

Level 2 能力者

11#
 楼主 Author| 发表于 Post on 2024-8-1 18:13:06 | 只看该作者 Only view this author
wzkchem5 发表于 2024-8-1 17:20
可能有一个问题是不应该加doht true,因为orca没有实现电离过程的HT效应的计算,只能对电子数守恒的吸收 ...

老师,我重新修改后这个问题解决了,但是输出结果中在显示Calculating energy of the ES at its geometry:后,直接报错ORCA finished by error termination in ORCA_ESD,现在ESD计算输入文件的结构与两个Hessian里的结构是一致的。不知道这样还有没有办法解决...麻烦老师了

202408011812484863..png (27.65 KB, 下载次数 Times of downloads: 16)

202408011812484863..png

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

12#
发表于 Post on 2024-8-1 18:19:48 | 只看该作者 Only view this author
SunriseOkk 发表于 2024-8-1 11:13
老师,我重新修改后这个问题解决了,但是输出结果中在显示Calculating energy of the ES at its geometry ...

K*K不应该这么接近0,建议检查你的阴离子和中性态的Hessian是不是用成同一个文件了
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

10

帖子

0

威望

33

eV
积分
43

Level 2 能力者

13#
 楼主 Author| 发表于 Post on 2024-8-1 18:47:45 | 只看该作者 Only view this author
wzkchem5 发表于 2024-8-1 18:19
K*K不应该这么接近0,建议检查你的阴离子和中性态的Hessian是不是用成同一个文件了

老师,我现在有点不太清楚该如何处理或者计算中性态的hessian文件,通过几次尝试发现,对中性态只做freq计算,结构与阴离子保持一致,就会出现K*K为0的情况。但是如果对中性态做opt+freq的话,结构变化,就会出现最开始的报错Error on ESD: Different number of vibrational modes between initial and final states

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

14#
发表于 Post on 2024-8-1 19:50:26 | 只看该作者 Only view this author
SunriseOkk 发表于 2024-8-1 11:47
老师,我现在有点不太清楚该如何处理或者计算中性态的hessian文件,通过几次尝试发现,对中性态只做freq ...

我才发现你用的是ORCA 5.0.4,可以试试ORCA 6.0.0行不行。如果还是不行,可以把输入文件以及两个.hess文件发给我,我们测试一下
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

10

帖子

0

威望

33

eV
积分
43

Level 2 能力者

15#
 楼主 Author| 发表于 Post on 2024-8-1 20:37:53 | 只看该作者 Only view this author
wzkchem5 发表于 2024-8-1 19:50
我才发现你用的是ORCA 5.0.4,可以试试ORCA 6.0.0行不行。如果还是不行,可以把输入文件以及两个.hess文 ...

老师,我这边计算也是用的6.0版本,不过我是参考的5.0版本手册写的计算文件...
附件里面分别是计算文件和两个Hess文件,其中中性态的文件也是计算了opt+freq之后的,请老师帮忙检查指导。
感谢老师!


202408012035242802..png (9.17 KB, 下载次数 Times of downloads: 12)

202408012035242802..png

PtNH3Cl3-0-gs023.hess

29.2 KB, 下载次数 Times of downloads: 1

PtNH3Cl3-1-gs044.hess

29.2 KB, 下载次数 Times of downloads: 1

PtNH3Cl3-orca-FCF02.inp

753 Bytes, 下载次数 Times of downloads: 2

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-24 05:23 , Processed in 0.212317 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list