计算化学公社

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

[新手求助] 高斯两种方式优化几何结构能量相差很大

[复制链接 Copy URL]

33

帖子

0

威望

267

eV
积分
300

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
我用高斯优化好一个二重态的结构,想算一下四重态的结构,但是我用两种方法计算的结果不一样,一种是考取二重态的.chk文件计算四重态的能量为(HF=-2787.9442),另一种是将二重态生成的.log文件重新生成坐标计算四重态(HF=-2787.9556),能量相差很大,哪一种方法更合理呢?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2016-10-6 23:51:13 | 只看该作者 Only view this author
意义不明。你要得到四重态结构,一般过程是gview打开二重态的out文件,保存成新输入文件,写上opt关键词,自旋多重度设4,然后做优化。
无论是gview打开二重态优化后的chk还是输出文件,看到的结构都是相同的,按上述操作得到的四重态的结构和能量也是相同的。
北京科音自然科学研究中心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

2494

帖子

11

威望

7052

eV
积分
9766

Level 6 (一方通行)

3#
发表于 Post on 2016-10-7 09:18:36 | 只看该作者 Only view this author
sobereva 发表于 2016-10-6 23:51
意义不明。你要得到四重态结构,一般过程是gview打开二重态的out文件,保存成新输入文件,写上opt关键词, ...

我明白她的意思,就是一种方法是大博士你说的方法,还有一种就是用.chk文件,只是把输入文件中的多重度由2变成4,利用.chk文件读取几何和初猜进行优化,这两种方式得到的能量不一致

33

帖子

0

威望

267

eV
积分
300

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2016-10-7 09:42:45 | 只看该作者 Only view this author
sobereva 发表于 2016-10-6 23:51
意义不明。你要得到四重态结构,一般过程是gview打开二重态的out文件,保存成新输入文件,写上opt关键词, ...

老师,我是想比较一下哪个基态更稳定,但是不同的计算方法计算四重态得到的能量结果差很大,我也觉得不管用哪种方法计算得到的能量正常应该是相同的,但是我现在遇到的这种情况不知道应该怎么解释,为什么同一个结构计算四重态时能量会不一样。

33

帖子

0

威望

267

eV
积分
300

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2016-10-7 09:43:07 | 只看该作者 Only view this author
我本是个娃娃 发表于 2016-10-7 09:18
我明白她的意思,就是一种方法是大博士你说的方法,还有一种就是用.chk文件,只是把输入文件中的多重度由 ...

多谢!

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

6#
发表于 Post on 2016-10-7 11:16:56 | 只看该作者 Only view this author
嘟嘟的小木鱼 发表于 2016-10-7 09:42
老师,我是想比较一下哪个基态更稳定,但是不同的计算方法计算四重态得到的能量结果差很大,我也觉得不管 ...

说能量前你先确保两个计算最终用的结构是相同的
北京科音自然科学研究中心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
积分
125141

管理员

公社社长

7#
发表于 Post on 2016-10-7 11:17:59 | 只看该作者 Only view this author
我本是个娃娃 发表于 2016-10-7 09:18
我明白她的意思,就是一种方法是大博士你说的方法,还有一种就是用.chk文件,只是把输入文件中的多重度由 ...

结果肯定是完全一致的,肯定是最基本的操作没弄对,或者数据都没读对。
直接从chk里读结构,和gview取优化最后一帧的结构,完全是相同的初猜。
北京科音自然科学研究中心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

298

帖子

1

威望

2595

eV
积分
2913

科音成员

8#
发表于 Post on 2016-10-7 11:38:02 | 只看该作者 Only view this author
可能从多重度2的chk中读取了guess=read,而另一个方法用默认的初始guess。
这样虽然多重度都是4,但是可能电子态不一样,所以能量不一样很正常。
从目前看,第二种方法得到的能量低(这里赞一下楼主,第一次提问就把能量贴出来了,而不是简单说哪个高哪个低),第二种方法的结果应该更好。
华北电力大学数理学院,理论与计算化学,团簇、表面的结构与反应机理。(招第一性原理计算,量子化学计算方向的教师、硕士/博士研究生)

2494

帖子

11

威望

7052

eV
积分
9766

Level 6 (一方通行)

9#
发表于 Post on 2016-10-7 11:57:29 | 只看该作者 Only view this author
sobereva 发表于 2016-10-7 11:17
结果肯定是完全一致的,肯定是最基本的操作没弄对,或者数据都没读对。
直接从chk里读结构,和gview取优 ...

我的理解是这样的,在二重态优化完成之后
方法一:利用.chk,geom=check guess=read fopt 0 4
方法二:读取优化好的二重态.log文件生成.com文件后 fopt 0 4这两者之间存在能量不一致,是否应该在方法一里去掉guess=read?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

10#
发表于 Post on 2016-10-7 13:15:42 | 只看该作者 Only view this author
我本是个娃娃 发表于 2016-10-7 11:57
我的理解是这样的,在二重态优化完成之后
方法一:利用.chk,geom=check guess=read fopt 0 4
方法二: ...

fopt应当写opt,f是多余的。

guess=read纯属多余,自旋多重度都不同,拿二重态的轨道当初猜不会有任何益处。
北京科音自然科学研究中心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

33

帖子

0

威望

267

eV
积分
300

Level 3 能力者

11#
 楼主 Author| 发表于 Post on 2016-10-7 14:12:44 | 只看该作者 Only view this author
helpme 发表于 2016-10-7 11:38
可能从多重度2的chk中读取了guess=read,而另一个方法用默认的初始guess。
这样虽然多重度都是4,但是可能 ...

我知道了,谢谢您的回复

33

帖子

0

威望

267

eV
积分
300

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2016-10-7 14:13:12 | 只看该作者 Only view this author
我本是个娃娃 发表于 2016-10-7 11:57
我的理解是这样的,在二重态优化完成之后
方法一:利用.chk,geom=check guess=read fopt 0 4
方法二: ...

抱拳,多谢!!!

33

帖子

0

威望

267

eV
积分
300

Level 3 能力者

13#
 楼主 Author| 发表于 Post on 2016-10-7 14:13:27 | 只看该作者 Only view this author
sobereva 发表于 2016-10-7 13:15
fopt应当写opt,f是多余的。

guess=read纯属多余,自旋多重度都不同,拿二重态的轨道当初猜不会有任何 ...

谢谢老师,我知道了。

2494

帖子

11

威望

7052

eV
积分
9766

Level 6 (一方通行)

14#
发表于 Post on 2016-10-7 14:15:25 | 只看该作者 Only view this author
sobereva 发表于 2016-10-7 13:15
fopt应当写opt,f是多余的。

guess=read纯属多余,自旋多重度都不同,拿二重态的轨道当初猜不会有任何 ...

老师,guess=read只有自旋多重度一致的时候使用还是什么情况下都不要使用?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

15#
发表于 Post on 2016-10-7 14:29:06 | 只看该作者 Only view this author
我本是个娃娃 发表于 2016-10-7 14:15
老师,guess=read只有自旋多重度一致的时候使用还是什么情况下都不要使用?

一般情况下,电子态和之前chk的情况一样的时候才用,否则并不会比默认的初猜更好而节约时间,有害无益。

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
辣椒丛中水果糖 + 1 我很赞同

查看全部评分 View all ratings

北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 05:15 , Processed in 0.928487 second(s), 21 queries , Gzip On.

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