计算化学公社

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

[Gaussian/gview] 关于不同版本高斯版本计算差异问题

[复制链接 Copy URL]

127

帖子

0

威望

338

eV
积分
465

Level 3 能力者

本帖最后由 GHL 于 2022-10-18 09:56 编辑

采用高斯16,linux版本计算计算同一化合物的 IP值,对比实验结果与参考文献结果趋势相反(高斯03 或09)


后期对比同一化合物,相同方法基组计算,高斯09 window与高斯16计算的 Sum of electronic and zero-point Energies结果相差比较大
请问一下是什么原因,有什么办法校正吗?




G09.gjf

3.41 KB, 下载次数 Times of downloads: 3

G16.gjf

3.41 KB, 下载次数 Times of downloads: 3

out.rar

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

127

帖子

0

威望

338

eV
积分
465

Level 3 能力者

19#
 楼主 Author| 发表于 Post on 2022-10-20 09:15:40 | 只看该作者 Only view this author
是的,差异很可能是基组不统一,谢谢大佬指点!谢谢!

1万

帖子

0

威望

9871

eV
积分
22113

Level 6 (一方通行)

18#
发表于 Post on 2022-10-18 19:50:21 | 只看该作者 Only view this author
GHL 发表于 2022-10-18 09:18
我仔细对比了一下最终优化结构,您所说的糖苷上不同构象,两者在视觉上是完全一致的,不知道从哪些方面可 ...

那就说明你这里计算差异的主要原因是基组没统一,而不是优化收敛到的构象不一致。
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
本团队长期招收研究生,有意者可私信联系

127

帖子

0

威望

338

eV
积分
465

Level 3 能力者

17#
 楼主 Author| 发表于 Post on 2022-10-18 16:18:24 | 只看该作者 Only view this author
乐平 发表于 2022-10-18 10:45
首先,你的差别并不只是版本不同……

你的 g09.gjf 输入文件里

我仔细对比了一下最终优化结构,您所说的糖苷上不同构象,两者在视觉上是完全一致的,不知道从哪些方面可以看出是顺反异构

1万

帖子

0

威望

9871

eV
积分
22113

Level 6 (一方通行)

16#
发表于 Post on 2022-10-18 16:03:39 | 只看该作者 Only view this author
GHL 发表于 2022-10-18 05:40
感谢大佬的详细解答,我明白了。十分感谢!!
就是有一点,怎样确保同一个结构 在结构优化后的结构一致 ...

初始结构定性一致(不要求完全一致),且理论级别严格一致即可。
如果你的分子构象数目多到即便很相似的初始结构也会优化出不同的收敛的结构,那么需要做的是用构象搜索把所有构象找全,而不仅仅是通过严格控制计算条件来让每次都收敛到同一构象,因为每次都收敛到同一个构象并不保证结果有意义,只有收敛到最稳定构象才是有意义的,而除了把比较低能量的构象找全了然后判断哪个能量最低以外,没有办法保证收敛到最稳定构象
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
本团队长期招收研究生,有意者可私信联系

127

帖子

0

威望

338

eV
积分
465

Level 3 能力者

15#
 楼主 Author| 发表于 Post on 2022-10-18 16:02:54 | 只看该作者 Only view this author
sobereva 发表于 2022-10-18 16:00
注意G09和G16的差异,基础量化班的ppt

谢谢大佬的指点,谢谢!

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

14#
发表于 Post on 2022-10-18 16:00:00 | 只看该作者 Only view this author
注意G09和G16的差异,基础量化班的ppt



也注意看此文
谈谈不同量子化学程序计算结果的差异问题
http://sobereva.com/573http://bbs.keinsci.com/thread-20122-1-1.html

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

13#
发表于 Post on 2022-10-18 15:58:05 | 只看该作者 Only view this author
乐平 发表于 2022-10-18 10:45
首先,你的差别并不只是版本不同……

你的 g09.gjf 输入文件里
自洽场收敛标准默认是 scf=tight,这些默认值都比 G09 的默认精度高。

G09开始默认就是scf=tight,G03是scf=sleazy
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

12#
发表于 Post on 2022-10-18 15:57:32 | 只看该作者 Only view this author
GHL 发表于 2022-10-18 09:54
两个gif文件初始结构不同,结构优化后,结果完全收敛,还会产生不同构象吗?
我记得这里面要分α和β糖 ...

分清楚gif和gjf
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

127

帖子

0

威望

338

eV
积分
465

Level 3 能力者

11#
 楼主 Author| 发表于 Post on 2022-10-18 12:40:12 | 只看该作者 Only view this author
乐平 发表于 2022-10-18 10:45
首先,你的差别并不只是版本不同……

你的 g09.gjf 输入文件里

感谢大佬的详细解答,我明白了。十分感谢!!
就是有一点,怎样确保同一个结构 在结构优化后的结构一致?
保证初始结构一致就可以吗?
有时候结果不收敛,我会在此基础上再做优化,不知道会不会造成结果的偏差

1155

帖子

0

威望

4430

eV
积分
5585

Level 6 (一方通行)

10#
发表于 Post on 2022-10-18 10:45:02 | 只看该作者 Only view this author
本帖最后由 乐平 于 2022-10-18 10:48 编辑
GHL 发表于 2022-10-18 00:26
这个是同一个结构,差别就是软件版本不同

首先,你的差别并不只是版本不同……

你的 g09.gjf 输入文件里
  1. # opt freq b3lyp/6-31g(d) scrf=(solvent=water) em=gd3bj
复制代码


而你的 g16.gjf 输入文件里
  1. # opt freq b3lyp/6-31+g(d) scrf=(solvent=water)  em=gd3bj
复制代码


看出区别了吗? 基组不同


其次,你的两个结构中戊碳糖的朝向也不同。
举个最简单了例子,顺式-1,3-丁二烯,反式-1,3-丁二烯,这两个分子用相同方法相同基组计算得到的能量是不同的。

再次,G16 默认用的积分格点精度是 int=ultrafine,自洽场收敛标准默认是 scf=tight,这些默认值都比 G09 的默认精度高。
而你的 G09 输入文件里关键词里并没有指定 int 和 scf 的精度。所以即使是相同的计算参数(更何况你的基组还不同……),得到的结果也会由于默认值的差异而不同。

127

帖子

0

威望

338

eV
积分
465

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2022-10-18 09:54:33 | 只看该作者 Only view this author
wzkchem5 发表于 2022-10-18 02:02
这两个输入文件的初始结构不同,可能是分别收敛到了两个不同的构象。可以检查一下收敛的结构是否相同

两个gif文件初始结构不同,结构优化后,结果完全收敛,还会产生不同构象吗?
我记得这里面要分α和β糖苷键,怎样能检查和确保构象固定? 大佬,我也是首次遇到这些问题。请多指教!!

127

帖子

0

威望

338

eV
积分
465

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2022-10-18 09:41:25 | 只看该作者 Only view this author
sobereva 发表于 2022-10-18 00:44
置顶的新社员必读贴明确说了,较大的文本型必须压缩后再上传。这种问题没有完整的输出文件根本没法说

...

好的,大佬,记住了!

1万

帖子

0

威望

9871

eV
积分
22113

Level 6 (一方通行)

7#
发表于 Post on 2022-10-18 02:02:50 | 只看该作者 Only view this author
GHL 发表于 2022-10-17 17:26
这个是同一个结构,差别就是软件版本不同

这两个输入文件的初始结构不同,可能是分别收敛到了两个不同的构象。可以检查一下收敛的结构是否相同
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
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

6#
发表于 Post on 2022-10-18 00:44:25 | 只看该作者 Only view this author
GHL 发表于 2022-10-18 00:30
文件太大了。无法上传,只能粘贴一部分结果了
个、G09结果为
Zero-point correction=                    ...

置顶的新社员必读贴明确说了,较大的文本型必须压缩后再上传。这种问题没有完整的输出文件根本没法说

有别人回复之前若需要对帖子进行修改、补充,应直接编辑原帖,不要通过回帖进行补充,这点在置顶的新社员必读贴里明确说了。
论坛不是QQ群,要一次性把问题交代全面,不要靠发额外的回帖代替编辑帖子

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

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

GMT+8, 2026-2-23 02:13 , Processed in 0.185063 second(s), 24 queries , Gzip On.

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