计算化学公社

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

[新手求助] 变换泛函优化加频率分析出现虚频

[复制链接 Copy URL]

36

帖子

0

威望

844

eV
积分
880

Level 4 (黑子)

Sob老师好,大家好。我在优化结构和频率分析时,采用B3LYP优化正常结束且无虚频。我想做一下benchmark,在变换泛函后优化正常结束,但是出现了虚频。我将B3LYP和变换基组的关键词和虚频展示在了下图。请教大家,谢谢大家。

虚频.png (9.98 KB, 下载次数 Times of downloads: 22)

虚频.png

PBE-freq.png (13.15 KB, 下载次数 Times of downloads: 33)

PBE-freq.png

PBE-OPT.png (10.17 KB, 下载次数 Times of downloads: 22)

PBE-OPT.png

B3LYP.png (10.79 KB, 下载次数 Times of downloads: 23)

B3LYP.png

36

帖子

0

威望

844

eV
积分
880

Level 4 (黑子)

2#
 楼主 Author| 发表于 Post on 2021-11-1 10:11:50 | 只看该作者 Only view this author
就是我用B3LYP一切正常,变换泛函后出现了虚频。

4289

帖子

4

威望

9546

eV
积分
13915

Level 6 (一方通行)

MOKIT开发者

3#
发表于 Post on 2021-11-1 10:58:10 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2021-11-1 11:00 编辑

如果操作正确,正常优化完就算有虚频 也顶多1-2个虚频。你是不是在pbepbe freq那里用错了坐标?没用上pbepbe opt优化好的坐标?另外,你的分子是否有很高的对称性,这个信息也要提供
自动做多参考态计算的程序MOKIT

36

帖子

0

威望

844

eV
积分
880

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2021-11-1 11:05:23 | 只看该作者 Only view this author
zjxitcc 发表于 2021-11-1 10:58
如果操作正确,正常优化完就算有虚频 也顶多1-2个虚频。你是不是在pbepbe freq那里用错了坐标?没用上pbepb ...

非常感谢您的回复,我用的是优化正常结束的结构,拿来再去分析频率出现了这样的虚频。用B3LYP时就没有出现这种情况。

1187

帖子

5

威望

2876

eV
积分
4163

Level 6 (一方通行)

5#
发表于 Post on 2021-11-1 11:11:46 | 只看该作者 Only view this author
优化和频率计算必须在同一级别下,你传的几个输入文件看不懂你是怎么算的,哪个输入对应第一个输出。你是不是B3LYP优化的结构PBE算的频率?另外分清PBE0(Gaussian里写作PBE1PBE)和PBE(Gaussian里写作PBEPBE),这是两个完全不同的泛函。你的chk文件名起的是PBE1PBE,关键词写的PBEPBE,这就很迷惑。

36

帖子

0

威望

844

eV
积分
880

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2021-11-1 11:16:56 | 只看该作者 Only view this author
snljty 发表于 2021-11-1 11:11
优化和频率计算必须在同一级别下,你传的几个输入文件看不懂你是怎么算的,哪个输入对应第一个输出。你是不 ...

嗯嗯,谢谢您。后面部分是我自己写的问题。
还有就是我用PBE算频率之前是先用PBE优化的,用PBE优化完的结构去分析的频率。就是我的频率分析和优化是保持一致的。就是图2和图3。

1万

帖子

0

威望

9870

eV
积分
22110

Level 6 (一方通行)

7#
发表于 Post on 2021-11-1 16:00:36 | 只看该作者 Only view this author
qian04551 发表于 2021-11-1 04:16
嗯嗯,谢谢您。后面部分是我自己写的问题。
还有就是我用PBE算频率之前是先用PBE优化的,用PBE优化完的 ...

建议把opt和freq写在同一个输入文件里再跑一遍,确保你确实把结构用对了
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

36

帖子

0

威望

844

eV
积分
880

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2021-11-1 20:02:38 | 只看该作者 Only view this author
wzkchem5 发表于 2021-11-1 16:00
建议把opt和freq写在同一个输入文件里再跑一遍,确保你确实把结构用对了

好的,谢谢您

4289

帖子

4

威望

9546

eV
积分
13915

Level 6 (一方通行)

MOKIT开发者

9#
发表于 Post on 2021-11-1 20:05:03 | 只看该作者 Only view this author

是否用错结构,还有其他的检查办法,例如查看freq任务输出文件里的Maximum Force那四行,如果是正确的结构,应该是前两个YES;如果前两行是NO,说明用错了结构。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
北大-陶豫 + 2 我很赞同

查看全部评分 View all ratings

自动做多参考态计算的程序MOKIT

36

帖子

0

威望

844

eV
积分
880

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2021-11-2 09:43:47 | 只看该作者 Only view this author
zjxitcc 发表于 2021-11-1 20:05
是否用错结构,还有其他的检查办法,例如查看freq任务输出文件里的Maximum Force那四行,如果是正确的结 ...

好滴,谢谢您呀。检查后发现结构没错,就尝试去虚频,现在正常了。

本版积分规则 Credits rule

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

GMT+8, 2026-2-22 18:22 , Processed in 0.169004 second(s), 24 queries , Gzip On.

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