计算化学公社

 找回密码 Forget password
 注册 Register
Views: 16858|回复 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
积分
125155

管理员

公社社长

24#
发表于 Post on 2019-10-28 01:36:32 | 只看该作者 Only view this author
Kellyyyy 发表于 2019-10-27 11:11
哦哦,谢谢你!那我再试试。

还想请问一下,模拟出来了物质的紫外吸收谱图之后,能够通过计算得到该物质 ...

模拟出的紫外谱的纵坐标本身就对应摩尔吸收系数epsilon,此文说了
使用Multiwfn绘制红外、拉曼、UV-Vis、ECD、VCD和ROA光谱图
http://sobereva.com/224
北京科音自然科学研究中心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 能力者

23#
 楼主 Author| 发表于 Post on 2019-10-27 11:11:58 | 只看该作者 Only view this author
哦哦,谢谢你!那我再试试。

还想请问一下,模拟出来了物质的紫外吸收谱图之后,能够通过计算得到该物质的摩尔吸光系数吗?

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

22#
发表于 Post on 2019-10-26 04:53:03 | 只看该作者 Only view this author
Kellyyyy 发表于 2019-10-24 23:18
您好!按照您给的算例已经把反应物和产物的能量算完了。在计算过渡态的能量时候出现了一些问题,最终计算终 ...

这不是报错
没完整输出文件看不出来
没有区别

没事别给氢加弥散函数,这里说了
谈谈弥散函数和“月份”基组
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 能力者

21#
 楼主 Author| 发表于 Post on 2019-10-24 23:18:10 | 只看该作者 Only view this author
您好!按照您给的算例已经把反应物和产物的能量算完了。在计算过渡态的能量时候出现了一些问题,最终计算终止:
No NMR shielding tensors so no spin-rotation constants.

不知道这是什么原因呢?
输入文件是:
#p opt freq ub3lyp/6-311++g(d,p) scrf=(solvent=water,smd) geom=connect
ivity em=gd3bj scale=0.9887 temperature=303.15

不知道计算过渡态的能量与计算其他稳定化合物的能量有什么区别吗?
还烦请sob大神给与解答,谢谢!

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

20#
发表于 Post on 2019-10-6 07:10:39 | 只看该作者 Only view this author
Kellyyyy 发表于 2019-10-5 23:12
抱歉,还是不大会在零点能校正和构型优化了的基础上,直接更换热力学校正因子。

%oldchk=C:%users\K\D ...

完全不知道你想做什么
用freqchk按照提示敲就完了
北京科音自然科学研究中心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 能力者

19#
 楼主 Author| 发表于 Post on 2019-10-5 23:12:27 | 只看该作者 Only view this author

抱歉,还是不大会在零点能校正和构型优化了的基础上,直接更换热力学校正因子。

%oldchk=C:\Users\K\Desktop\20191003\irc1\1narawopt-2.chk
N
303.15
1
1.0102


直接就终止了。
Route card not found.
Error termination

所以,求助

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

18#
发表于 Post on 2019-10-5 06:54:06 | 只看该作者 Only view this author
Kellyyyy 发表于 2019-10-5 06:20
这个可以在A01上使用哈?

just try
北京科音自然科学研究中心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 能力者

17#
 楼主 Author| 发表于 Post on 2019-10-5 06:20:20 | 只看该作者 Only view this author
sobereva 发表于 2019-10-5 05:21
可以。或者用wB97XD

这个可以在A01上使用哈?

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

16#
发表于 Post on 2019-10-5 05:21:44 | 只看该作者 Only view this author
Kellyyyy 发表于 2019-10-4 21:16
还有想请教一下:
用频率因子对热力学参数进行校正的时候,你说到最好的方法是用Gaussian自带的freqchk ...

没法说得更清楚了
做个freq任务显然就有含有Hessian的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

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

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

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没有

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的前几行包括调用等怎么写)?

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吗?

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

GMT+8, 2026-2-24 13:45 , Processed in 0.164335 second(s), 21 queries , Gzip On.

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