Boltzmann分布计算分子不同构象比例
更新时间:2023-06-09 09:48:01 阅读量: 实用文档 文档下载
- Boltzmann分布推荐度:
- 相关推荐
Boltzman 分布计算构象比例
[原创
]根据Boltzmann分布计算分子不同构象所占比例
根据Boltzmann分布计算分子不同构象所占比例
文/Sobereva 2012-Oct-20
一种分子往往有很多构象,每种构象能量各不相同。在平衡状态下,各种构象出现的比例也是不同的。在很多文献里都给出根据Boltzmann分布计算的分子在不同构象下的比例,也经常看到有人问怎么算。实际上计算非常简单,此文简单说明一下。用下面这张图的6-31G*的数据作为例子,来自J. Nat. Prod. 2011, 74, 1061。
Boltzmann分布的概念在一般的物理化学书上的统计热力学部分都讲过,可以写为
Boltzman 分布计算构象比例
其中p是所占比例,i是构象编号,n_i是处于第i构象的分子数,E是指构象的能量。T是温度(开尔文),R是理想气体常数。Q称作配分函数。
准确计算分子的能量是极其困难的,只有很小的体系,用很高精度量化方法才能得到定量准确的结果。好在Boltzmann公式按如下方式可以等价地写为只依赖于不同构象间相对能量的形式,而相对能量比较容易得到定量可靠的结果。
式中E_Ref代表所有构象中能量最低值(参考值),ΔE是相对值,由参考值对应的常数项C在计算p的时候同时作为分子分母而消掉了。Q下标上的Relat代表Relative的含义。
根据这个式子,在Excel里就可以很容易地计算构象分布比例了,如下图所示,Percent一列是最终结果,和前面文献图上所列的结果相吻合(微小偏差是数值精度原因或者所取温度略不同导致的)。过程是:先把文字、温度和相对能量填好。然后在C4上填=EXP(-B4/($B$1*0.0019858995)),并将此单元格复制到C5~C11上,其中0.0019858995是R在kcal/mol/K单位下的数值。然后D1填=SUM(C4:C11)。之后D4填=C4/$D$1,并将此单元格复制到D5~D11上,就得到结果了。
Boltzman 分布计算构象比例
需要说明的是,前面所谓的构象的能量是指自由能。所以,必须先优化到相应构象结构,计算电子能量,再做振动分析,把热力学校正能加上去得到自由能。在什么温度下计算分布比例,就应当在什么温度下计算热力学校正能。另外,若想根据Boltzmann公式定量精确地计算各个构象的比例必须把自由能算得很准,推荐用热力学组合方法来算,如Gn系列、CBS系列等方法。
假设体系只有相对能量为0, 0.17, 0.37, 0.66, 0.75 kcal/mol这五种构象,按照前文的方法计算,看看从低温到高温过程中它们所占比例的变化:
从图上看,温度很低时,分子基本都处于能量最低的构象,在0K时100%处于这种构象。由于几种构象能量差异并不大,在室温级别下其它构象所占比例就已经不可忽视了。当温度无穷高,能量差异已经显不出来了,各种构象所占比例就相等了。不过实际上这图并没有什么物理意义,首先,在计算时忽视了不同温度下各构象自由能的变化;其次,当温度较高时,不同构象之间的界限已经很含糊了,较高的动能可以让分子构象随意变换而几乎无视构象间势垒的存在,因此讨论哪个构象占比例多少已经没有太大意义;而温度极高时,分子就彻底解离了。
Boltzman 分布计算构象比例
正在阅读:
C++程序设计论文模板10-24
第1章三角函数测试1(苏教版必修4)03-05
半月谈中五篇文章12-30
小学六年级上册语文期中试卷及答案05-15
九年级语文上册18《杨修之死》(第二课时)教学案新人教版04-26
分段函数与复合函数01-26
荆轲刺秦王英语剧本05-30
挂企工作总结09-27
地产企业如何做好产品设计管理05-15
- 教学能力大赛决赛获奖-教学实施报告-(完整图文版)
- 互联网+数据中心行业分析报告
- 2017上海杨浦区高三一模数学试题及答案
- 招商部差旅接待管理制度(4-25)
- 学生游玩安全注意事项
- 学生信息管理系统(文档模板供参考)
- 叉车门架有限元分析及系统设计
- 2014帮助残疾人志愿者服务情况记录
- 叶绿体中色素的提取和分离实验
- 中国食物成分表2020年最新权威完整改进版
- 推动国土资源领域生态文明建设
- 给水管道冲洗和消毒记录
- 计算机软件专业自我评价
- 高中数学必修1-5知识点归纳
- 2018-2022年中国第五代移动通信技术(5G)产业深度分析及发展前景研究报告发展趋势(目录)
- 生产车间巡查制度
- 2018版中国光热发电行业深度研究报告目录
- (通用)2019年中考数学总复习 第一章 第四节 数的开方与二次根式课件
- 2017_2018学年高中语文第二单元第4课说数课件粤教版
- 上市新药Lumateperone(卢美哌隆)合成检索总结报告
- 同构
- Boltzmann
- 分布
- 分子
- 比例
- 计算
- 2021届云南、四川、贵州、西藏四省名校高三第一次大联考化学试题(解析版)_1
- 梁惠星教授“民法解释学”讲义
- 漳浦县新农村建设调查报告
- 英语与汉语的差别
- 第四章 建设工程合同管理
- 霍桑_红字_的女性主义解读
- 2012年中级经济师考试人力资源考前模拟试题3
- 广西百色市2008年初中毕业暨升学考试试卷
- 退火对Y2O3薄膜结构和电学性能的影响
- 2015最新6.3共产党领导的多党合作和政治协商制度:中国特色社会主义政党制度
- 单片机实验——秒表 (详细步骤)
- 常识知识分类汇总
- 战后对纳粹罪行的审判与德国反省历史的自觉意识
- 小学二年级第二学期数学期末试题2013
- 济南大学信号与信息处理研究生课程安排
- 隧道掘进水压爆破方案
- 小学四年级作文题目我最熟悉的人的四年级优秀作文
- 浅谈我国的耕地和基本农田保护制度7.15
- 3、帷幕注浆施工工艺工法
- 学前班数学单数和双数教学设计