计算化学公社

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

[Gaussian/gview] 求助高斯的freqchk命令计算不同温度下的吉布斯自由能的问题

[复制链接 Copy URL]

24

帖子

0

威望

101

eV
积分
125

Level 2 能力者

各位老师,本人还是量化计算中的小白,有个问题想要请教。就是在利用高斯的freqchk命令计算得到的不同温度下的吉布斯自由能结果,为什么和直接利用temperature=(目标温度)关键词得到的吉布斯自由能的结果不一样呢?

4114

帖子

4

威望

8875

eV
积分
13069

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2020-11-12 13:05:00 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2020-11-12 13:07 编辑

最好能晒举出具体的数据,以及你的命令行参数。有的新手看到小数点后第6位有变化,也认为是“不一样”
如果还能晒出你操作时用的.chk或.fchk文件,那大家便可迅速找到问题所在

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

524

帖子

3

威望

6224

eV
积分
6808

Level 6 (一方通行)

3#
发表于 Post on 2020-11-12 16:27:19 | 只看该作者 Only view this author
zjxitcc 发表于 2020-11-12 13:05
最好能晒举出具体的数据,以及你的命令行参数。有的新手看到小数点后第6位有变化,也认为是“不一样”
如 ...

新手看到小数点后第6位有变化,也认为是“不一样”真的太真实了
恍惚月余,深谙人与人之间的差距。以后还应努力学习,才能与强者比肩。

24

帖子

0

威望

101

eV
积分
125

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2020-11-12 16:35:06 | 只看该作者 Only view this author
本帖最后由 wuyanzhou 于 2020-11-12 16:37 编辑
zjxitcc 发表于 2020-11-12 13:05
最好能晒举出具体的数据,以及你的命令行参数。有的新手看到小数点后第6位有变化,也认为是“不一样”
如 ...

file:///C:/Users/2019620960/Desktop/%E5%BE%AE%E4%BF%A1%E5%9B%BE%E7%89%87_20201112163045.png

微信图片_20201112163053.png (169.73 KB, 下载次数 Times of downloads: 50)

直接设置温度关键词得到的G,H

直接设置温度关键词得到的G,H

微信图片_20201112163045.png (109.21 KB, 下载次数 Times of downloads: 32)

利用freqchk关键词得到的G,H

利用freqchk关键词得到的G,H

微信图片_20201112163004.png (67.91 KB, 下载次数 Times of downloads: 23)

没有设置温度,直接得到的298k下的G,H值

没有设置温度,直接得到的298k下的G,H值

24

帖子

0

威望

101

eV
积分
125

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2020-11-12 16:39:12 | 只看该作者 Only view this author
ABetaCarw 发表于 2020-11-12 16:27
新手看到小数点后第6位有变化,也认为是“不一样”真的太真实了

老师,好像不只是小数点后几位的问题,直接差了很多

24

帖子

0

威望

101

eV
积分
125

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2020-11-12 16:40:54 | 只看该作者 Only view this author
zjxitcc 发表于 2020-11-12 13:05
最好能晒举出具体的数据,以及你的命令行参数。有的新手看到小数点后第6位有变化,也认为是“不一样”
如 ...

老师,下面的三张图片,分别对应直接设置关键词temperature=313得到的。第二张对应freqchk得到的,第三张是常温298k的数据

1万

帖子

0

威望

9004

eV
积分
20764

Level 6 (一方通行)

7#
发表于 Post on 2020-11-12 19:11:23 | 只看该作者 Only view this author
wuyanzhou 发表于 2020-11-12 16:40
老师,下面的三张图片,分别对应直接设置关键词temperature=313得到的。第二张对应freqchk得到的,第三张 ...

你用freqchk的时候读Hessian信息没有读成功,因此计算结果没有包括振动的贡献。
仔细重新研究一下freqchk的用法,确保该读的文件的格式都是对的,放在了正确的位置,再试
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?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

4114

帖子

4

威望

8875

eV
积分
13069

Level 6 (一方通行)

MOKIT开发者

8#
发表于 Post on 2020-11-12 22:00:17 | 只看该作者 Only view this author
wuyanzhou 发表于 2020-11-12 16:40
老师,下面的三张图片,分别对应直接设置关键词temperature=313得到的。第二张对应freqchk得到的,第三张 ...

freqchk的操作对象是 加了关键词Freq算出的chk或fchk文件,你不会是随便算了个单点->就扔给freqchk操作了吧
自动做多参考态计算的程序MOKIT

本版积分规则 Credits rule

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

GMT+8, 2025-8-17 23:31 , Processed in 0.185802 second(s), 23 queries , Gzip On.

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