计算化学公社

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

[新手求助] 关于谈谈隐式溶剂模型下溶解自由能和体系自由能的计算“的例子中的问题

[复制链接 Copy URL]

111

帖子

0

威望

206

eV
积分
317

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
求助:
Sob老师在中的最后给出了一个计算的例子,其中2、3、4都是依据1中优化出来的结构来计算的,是不是因为环氧乙烷作为中性分子受到溶剂对其分子构型的影响小至可以忽略的原因呢? 对于一般的离子体系,很有可能在真空中和在溶剂中的最优构型完全不同,是不是最好在计算气相标况自由能的时候基于真空来优化分子构性;而计算溶解自由能的时候,先对于真空条件下优化好的分子计算其自由能,然后再在加了溶剂模型下重新优化分子构型之后,再来计算其电子能量呢?最后再将二者相加再加1.89kcal/mol,这样算对吗?


6万

帖子

99

威望

6万

eV
积分
125174

管理员

公社社长

2#
发表于 Post on 2019-10-2 16:51:25 | 只看该作者 Only view this author
想多了
你这样做结果更差。文中的做法是严格的
北京科音自然科学研究中心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

111

帖子

0

威望

206

eV
积分
317

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2019-10-2 20:50:00 | 只看该作者 Only view this author
sobereva 发表于 2019-10-2 16:51
想多了
你这样做结果更差。文中的做法是严格的

噢噢噢。不好意思,才看到:溶质分子在溶剂环境下的总能量Etot=溶剂化能deltaG(solv)+相同构象但是在真空中的能量Evac。

不好意思!

111

帖子

0

威望

206

eV
积分
317

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2019-10-2 21:12:37 | 只看该作者 Only view this author
sobereva 发表于 2019-10-2 16:51
想多了
你这样做结果更差。文中的做法是严格的

对了,在您给的例子里没有看到关于频率因子的关键词(scale=...),是不需要么?

111

帖子

0

威望

206

eV
积分
317

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2019-10-2 21:30:08 | 只看该作者 Only view this author
Kellyyyy 发表于 2019-10-2 20:50
噢噢噢。不好意思,才看到:溶质分子在溶剂环境下的总能量Etot=溶剂化能deltaG(solv)+相同构象但是在真空 ...

还有就是对于温度的考虑。只需要在1中加上temperature吗?
后面算真空下的电子能量和溶解自由能就不需要添加该关键词了吗?那结果也会包含温度对热力学的影响吗?

6万

帖子

99

威望

6万

eV
积分
125174

管理员

公社社长

6#
发表于 Post on 2019-10-3 03:15:08 | 只看该作者 Only view this author
Kellyyyy 发表于 2019-10-2 21:12
对了,在您给的例子里没有看到关于频率因子的关键词(scale=...),是不需要么?
影响不是特别大,就忽略了。加的话原理上有好处
北京科音自然科学研究中心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
积分
125174

管理员

公社社长

7#
发表于 Post on 2019-10-3 03:24:07 | 只看该作者 Only view this author
Kellyyyy 发表于 2019-10-2 21:30
还有就是对于温度的考虑。只需要在1中加上temperature吗?
后面算真空下的电子能量和溶解自由能就不需要 ...

默认是常温

严格来说没有办法计算非常温的溶剂下的自由能,因为这些隐式溶剂模型都是对常温参数化的,就算你气相自由能是其它温度的,溶解自由能这部分也没法严格算其它温度的

温度与电子能量的关系仔细看
谈谈温度、压力、同位素设定对量子化学计算结果产生的影响
http://sobereva.com/423http://bbs.keinsci.com/thread-10159-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

111

帖子

0

威望

206

eV
积分
317

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2019-10-3 11:36:40 | 只看该作者 Only view this author
sobereva 发表于 2019-10-3 03:15
影响不是特别大,就忽略了。加的话原理上有好处

哦哦,知道了

111

帖子

0

威望

206

eV
积分
317

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2019-10-3 11:41:14 | 只看该作者 Only view this author
sobereva 发表于 2019-10-3 03:24
默认是常温

严格来说没有办法计算非常温的溶剂下的自由能,因为这些隐式溶剂模型都是对常温参数化的, ...

噢噢噢。

仿照您的例子想做一个计算,输入文件指令如下:

#P UB3LYP/6-311++G** em=GD3BJ opt freq scrf(SMD,solvent=water)

但是一传到集群上的Linux系统上计算,立即错误终止:

QPErr --- A syntax error was detected in the input line.
#P UB3LYP/6-311++G** em=GD3BJ opt freq s

Last state="GCL"
TCursr= 1046 LCursr= 23

这是为什么呀?

6万

帖子

99

威望

6万

eV
积分
125174

管理员

公社社长

10#
发表于 Post on 2019-10-4 04:44:34 | 只看该作者 Only view this author
Kellyyyy 发表于 2019-10-3 11:41
噢噢噢。

仿照您的例子想做一个计算,输入文件指令如下:

括号打成全角了

另外如果用的是G09 D.01以前的版本,没法用D3校正

没事甭给氢用弥散,这里明确说了
谈谈弥散函数和“月份”基组
http://sobereva.com/119
北京科音自然科学研究中心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

111

帖子

0

威望

206

eV
积分
317

Level 3 能力者

11#
 楼主 Author| 发表于 Post on 2019-10-4 20:47:22 | 只看该作者 Only view this author
sobereva 发表于 2019-10-4 04:44
括号打成全角了

另外如果用的是G09 D.01以前的版本,没法用D3校正

您说得对,集群上的版本是A01,不支持D3校正。。。
那么除了用B3LYP+D3以外,还有其它A01也支持的方式吗?

111

帖子

0

威望

206

eV
积分
317

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2019-10-4 21:04:20 | 只看该作者 Only view this author
sobereva 发表于 2019-10-4 04:44
括号打成全角了

另外如果用的是G09 D.01以前的版本,没法用D3校正

那么可以直接用M06-2X代替b3lyp+d3来计算ZPE吗?

111

帖子

0

威望

206

eV
积分
317

Level 3 能力者

13#
 楼主 Author| 发表于 Post on 2019-10-4 21:16:52 | 只看该作者 Only view this author
sobereva 发表于 2019-10-4 04:44
括号打成全角了

另外如果用的是G09 D.01以前的版本,没法用D3校正

还有想请教一下:
用频率因子对热力学参数进行校正的时候,你说到最好的方法是用Gaussian自带的freqchk工具,它会从chk文件中读取储存的Hessian矩阵并作分析。启动它之后,输入chk文件路径,然后依次输入
N    //不输出Hyperchem文件
323   //温度。如果输入0则为默认的298.15K
1     //压力。如果输入0则为默认的1atm
0.8905   //校正因子
[回车]
[回车]

能否举一个简单的例子说明如何用Gaussian自带的freqchk工具(这个是什么)呢?如何在上一个频率因子校正完频率之后的基础上重新输入另一个热力学参数校正因子再算呢(即gjf的前几行包括调用等怎么写)?

139

帖子

1

威望

1436

eV
积分
1595

Level 5 (御坂)

14#
发表于 Post on 2019-10-5 05:07:39 | 只看该作者 Only view this author
Kellyyyy 发表于 2019-10-4 20:47
您说得对,集群上的版本是A01,不支持D3校正。。。
那么除了用B3LYP+D3以外,还有其它A01也支持的方式吗 ...

给自己目录下装一个高级版本就好了 m062x我记得也是晚一点的版本才有的 估计A01没有

6万

帖子

99

威望

6万

eV
积分
125174

管理员

公社社长

15#
发表于 Post on 2019-10-5 05:21:05 | 只看该作者 Only view this author
Kellyyyy 发表于 2019-10-4 21:04
那么可以直接用M06-2X代替b3lyp+d3来计算ZPE吗?

可以。或者用wB97XD
北京科音自然科学研究中心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-25 07:51 , Processed in 0.823589 second(s), 20 queries , Gzip On.

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