计算化学公社

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

[新手求助] 基于玻尔兹曼分布算出来的构象比例跟文献中不一样

[复制链接 Copy URL]

588

帖子

0

威望

3310

eV
积分
3898

Level 5 (御坂)

本帖最后由 xujian 于 2016-12-12 20:37 编辑

文献中的表格如下




我根据如下公式(使用的是B3LYP的detaG):

得到298K时比例分别是68%、16%、10%和6%,478K时比例分别为52%、21%、16%、11%。跟文献中差距很大,不知道哪里错了,请高手指导一下,谢谢!

188

帖子

0

威望

4072

eV
积分
4260

Level 6 (一方通行)

2#
发表于 Post on 2016-12-12 21:30:11 | 只看该作者 Only view this author
sob老师有个专门的帖子写玻尔兹曼分布的http://sobereva.com/165,一般算能量的方法是:先计算电子的能量然后计算频率把热力学矫正加进去,如果想精确计算的话需要用热力学组合的方法,我记得在老师的培训班中讲过:庸人的做法,体面人的做法....细节不记得了
If anything is impossible? It should be done!!

6万

帖子

99

威望

6万

eV
积分
125151

管理员

公社社长

3#
发表于 Post on 2016-12-12 23:04:24 | 只看该作者 Only view this author
波尔兹曼分布对于自由能的计算准确度极为敏感,自由能的一丁点不同都会导致比例相差不少
可能你和文献在计算上在一些细节部分处理不同,导致结果不同。文献也不代表就一定是对的

正确的计算气相和溶液中自由能的方法在此文里有详细说明:
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327

对于小体系,建议直接用热力学组合方法如CBS-QB3计算自由能,又方便又准。
北京科音自然科学研究中心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

588

帖子

0

威望

3310

eV
积分
3898

Level 5 (御坂)

4#
 楼主 Author| 发表于 Post on 2016-12-13 08:34:11 | 只看该作者 Only view this author
sobereva 发表于 2016-12-12 23:04
波尔兹曼分布对于自由能的计算准确度极为敏感,自由能的一丁点不同都会导致比例相差不少
可能你和文献在计 ...

谢谢sob老师

87

帖子

0

威望

376

eV
积分
463

Level 3 能力者

5#
发表于 Post on 2017-11-25 16:03:58 | 只看该作者 Only view this author
sobereva 发表于 2016-12-12 23:04
波尔兹曼分布对于自由能的计算准确度极为敏感,自由能的一丁点不同都会导致比例相差不少
可能你和文献在计 ...

老师,如果体系很大,没办法计算自由能,是不是考虑玻尔兹曼分布就不合理了?正如您说的自由能的一丁点不同都会导致比例相差不少

6万

帖子

99

威望

6万

eV
积分
125151

管理员

公社社长

6#
发表于 Post on 2017-11-25 16:12:51 | 只看该作者 Only view this author
sofia 发表于 2017-11-25 16:03
老师,如果体系很大,没办法计算自由能,是不是考虑玻尔兹曼分布就不合理了?正如您说的自由能的一丁点不 ...

如果振动分析确实做不动,那就没辙了,或者用电子能量近似作为自由能来计算分布比例
北京科音自然科学研究中心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

87

帖子

0

威望

376

eV
积分
463

Level 3 能力者

7#
发表于 Post on 2017-11-25 20:18:33 | 只看该作者 Only view this author
sobereva 发表于 2017-11-25 16:12
如果振动分析确实做不动,那就没辙了,或者用电子能量近似作为自由能来计算分布比例

老师,您帖子中玻尔兹曼分布的公式的出处能够给出呢(是否有相应的参考文献),我在物化书中没有找到呢?
谢谢老师!

4289

帖子

4

威望

9546

eV
积分
13915

Level 6 (一方通行)

MOKIT开发者

8#
发表于 Post on 2017-11-25 21:58:15 | 只看该作者 Only view this author
sofia 发表于 2017-11-25 20:18
老师,您帖子中玻尔兹曼分布的公式的出处能够给出呢(是否有相应的参考文献),我在物化书中没有找到呢? ...

那你看的一定是假的物化书哈哈哈!随便找本物化或者统计热力学的书,上面都有写。只可能是找错书或者没认真看
自动做多参考态计算的程序MOKIT

87

帖子

0

威望

376

eV
积分
463

Level 3 能力者

9#
发表于 Post on 2017-11-26 19:35:18 | 只看该作者 Only view this author
zjxitcc 发表于 2017-11-25 21:58
那你看的一定是假的物化书哈哈哈!随便找本物化或者统计热力学的书,上面都有写。只可能是找错书或者没认 ...

在物化书上找到了一部分,没有找到一模一样的公式,不知道在文章中怎么引用,而且书上写的是KT而不是RT

204

帖子

0

威望

2733

eV
积分
2937

Level 5 (御坂)

10#
发表于 Post on 2017-11-27 20:50:39 | 只看该作者 Only view this author
sofia 发表于 2017-11-26 19:35
在物化书上找到了一部分,没有找到一模一样的公式,不知道在文章中怎么引用,而且书上写的是KT而不是RT

这个取决于E用的单位,/mol的话就用R,单个分子的话就用k

87

帖子

0

威望

376

eV
积分
463

Level 3 能力者

11#
发表于 Post on 2017-11-28 19:53:09 | 只看该作者 Only view this author
赵云跳槽 发表于 2017-11-27 20:50
这个取决于E用的单位,/mol的话就用R,单个分子的话就用k

非常感谢!

本版积分规则 Credits rule

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

GMT+8, 2026-2-23 06:37 , Processed in 0.165977 second(s), 23 queries , Gzip On.

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