您所在位置： &
&nbsp&&nbsp&nbsp&&nbsp
分子体积的计算.doc 5页
本文档一共被下载：
次 ,您可全文免费在线阅读后下载本文档。
下载提示
1.本站不保证该用户上传的文档完整性，不预览、不比对内容而直接下载产生的反悔问题本站不予受理。
2.该文档所得收入(下载+内容+预览三)归上传者、原创者。
3.登录后可充值，立即自动返金币，充值渠道很便利
需要金币：100 &&
分子体积的计算,高斯计算分子体积,圆柱体积的计算公式,圆体积的计算公式,体积的计算公式,空气的成分按体积计算,圆球的体积计算公式,分子势能与体积的关系,相对分子质量的计算,圆锥体积的计算公式
你可能关注的文档：
··········
··········
谈谈分子体积的计算文/Sobereva2011-Sep-20我曾在网上数次看到有人问怎么算分子体积，以及关于Gaussian的volume关键词使用问题，从提问以及很多人的回答上看有不少人对这问题存在错误的认识。本文就来谈谈这个问题。首先要认识到，分子体积不是一个可观测量，在计算方法上也不可能有唯一的定义，因此“怎么计算分子体积”这个问题本身就是不严格的。分子体积就是分子表面内部空间的体积，由于分子表面也没有唯一的定义，所以不同的分子表面的定义就会给出不同的分子体积定义。首先看一个简单分子的例子红色区域是每个原子的范德华球（以原子核为中心，半径为范德华半径的球体）的叠加，这片区域就是分子的范德华体积，其表面也就被称作范德华表面。图中蓝球代表作为探针的溶剂分子（显然溶剂实际形状并不是球形，所以这个蓝球半径是“等效”的，在计算程序中通常是可调参数），让这个蓝球紧贴着分子范德华表面在各处滚一遍，就产生了诸如图中绿色的轨迹，对应的表面叫做Connolly表面，由于溶剂分子不能触及到这个表面内的空间，所以也被叫做solvent-excluded表面，其内部区域的体积就叫做solvent-excluded体积。图中黄色是蓝球滚动时蓝球中心经过的表面，这个表面叫溶剂可及表面（其表面积就是所谓的SASA）。最常用的体积就是范德华体积。实际上范德华体积也有很多不同定义，上面介绍的原子范德华球叠加是比较简单的定义方式，通过解析几何的方法就能算出来，也可以用后面谈到的蒙特卡罗方法。这种定义存在两个缺点：(1)原子范德华半径没有唯一定义，不同研究者给出的范德华半径存在分歧。而且对于一些金属体系没有能用的范德华半径数据。(2)没有考虑到电子效应，成键导致的电子转移、极化效果都被忽略了。Bader提出过一种比较严格的范德华体积的定义，消除了前述定义的弊端，并已被广泛接受，也就是若分子处于气相，则将电子密度为0.001的等值面作为范德华表面，这种表面通常能够囊括分子98%以上的电子密度，这种范德华体积通常比范德华球叠加得到的范德华体积要大。而对于处于凝聚态的分子，考虑到分子间挤压，建议用电子密度为0.002的等值面作为范德华表面（显然，其体积小于0.001等值面对应的体积）。计算范德华体积最常用的办法是蒙特卡罗方法。首先，设立一个矩形盒子，将整个分子扩住，并且各个方向上都预留一定空间以避免将范德华表面截断，记这个盒子的体积为V_box。然后，在盒子里随机分布m个点，依次检验这些点是否符合条件。对于范德华球叠加方式的定义，如果当前点与任何一个原子核的距离小于相应原子的范德华半径，则认为此点符合条件；而对于Bader的定义，若当前点的电子密度大于阈值（0.001或0.002），就认为符合条件。假设最后有n个点符合条件，那么分子的范德华体积就是n/m*V_box。如果m较小，那么算出来的体积是不精确的，想要增加精度，就必须增加m。这就像人口普查，统计的人数越多结果越能反映实际国情。当然，分子越大，就需要越多的m，才能保证平均每单位体积内随机点的数目不变，即保证精度不变。由于每次用蒙特卡罗方法计算体积时随机分布的点的位置都是不同的，因此符合条件的点数也会不同，故算出来的范德华体积的数值每次肯定会不同。然而m越大，结果的波动就会越小。曾有人问什么程序计算的范德华体积更精确。这个问题太含糊，没法回答。只能概括地说，如果想精确计算范德华体积，就需要使用Bader的定义，用较高的随机点密度去做蒙特卡罗计算，并且用比较准确的电子密度。在Gaussian中，专门有个volume关键词用于通过蒙特卡罗方法计算Bader定义的分子范德华体积（也用来估算Onsager溶剂模型下要用的溶质半径），同时有两个相关的IOp。IOp(6/46=n)用来将电子密度等值面设为n*0.0001，默认算的是0.001的。IOp(6/45)可以设定每立方波尔平均有多少个随机点用于蒙特卡罗计算，默认是20。这个数值明显偏小，结果不够精确，也导致很多人发现每次计算体积结果都差异很大而质疑Gaussian计算体积的功能是不是可靠。实际上，将这个数值调大1~2个数量级后，体积计算精度会大有改善，结果波动也会明显小很多，当然，也会更耗时。我看到网上有人想通过让Gaussian计算100次体积取平均值以提高计算结果的可靠性，这是笨方法，还得写脚本去批量执行和处理。实际上直接将默认的随机点密度提高100倍，即IOp(6/45=2000)，算一次就行了，完全是等价的。用Volume关键词时有两点需要注意。第一点是对于后HF方法不要忘了加density关键词，否则密度是HF下的，但实际上对结果的影响很小，一般还不如蒙特卡罗方法的随机性的影响大。另外就是输出的结果单位问题，会有诸如这样的输出?Molarvolume
正在加载中，请稍后...小木虫 --- 600万学术达人喜爱的学术科研平台
热门搜索：
&&查看话题
请教高斯09软件如何计算分子体积
问一下高斯09软件如何计算分子的体积，把分子式画好以后选择高斯计算，然后出现Gaussian Calculation Setup，下面就不知道怎么弄了。Job Type应该选择什么？没有关于体积的单词啊？然后方法里面应该怎么设置？跪求大神指点，越详细越好，急用啊！都急出shi了~~~~~:cry::cry::cry::cry:
这年头什么都变得廉价？动辄急出翔？静下心来做学问，好好看看volume 这关键字怎用
能给一下具体的操作吗？~~~ 感谢
我又不是搞这个的只是用一下软件而已，着急就是着急，不喜欢不要回复。
学术必备与600万学术达人在线互动！
扫描下载送金币& 【求助】gaussian计算分子的表面增强拉曼光谱
需要登录并加入本群才可以回复和发新贴
【求助】gaussian计算分子的表面增强拉曼光谱
UID 131725
信誉分 100
可用分 5321
阅读权限 255
注册 状态 离线
【求助】gaussian计算分子的表面增强拉曼光谱
最近看一些文献，发现在用gaussian计算分子的表面增强拉曼光谱时，模型选择有很多方式，不知有什么考虑，本人是个菜鸟，没什么理论功底:shuai:，比较迷惑。
举例如下：比如是苯酚（C6H6O）吸附在银纳米颗粒上吧，在选择计算模型的时候主要有这样两种计算模式（1）苯酚不变外加一个银（2）苯酚中——OH中的氢消失，再加一个银。当然在此基础上还有两种计算模式1）苯酚不变外加两个银（2）苯酚中——OH中的氢消失，再加两个银。以此类推可以加3、4、5个银。
我的问题：从实验上如何解释挂1个银时，苯酚不变或者苯酚少一个氢啊，从实验角度上哪个模型更合理。增加了若干个银后，还是至少一个氢，这样的模型对吗？或者从实验角度上讲该如何设置模型？
UID 129814
信誉分 100
可用分 10775
阅读权限 255
注册 状态 离线
呵呵 这具体的我也不知道，我的方法是所有方法都用，都计算出来后，和实验对比一下，哪个符合，就是那个了。我理论基础也差所有就用这种办法啊。
UID 131600
信誉分 100
可用分 7817
阅读权限 255
注册 状态 离线
很想帮你 但我只是学习过拉曼光谱仪的使用 自己也不懂 不好意思啊
UID 128420
信誉分 102
可用分 10670
阅读权限 255
注册 状态 离线
我只做表面增强实验，没做过计算，对计算不懂，不好意思
UID 129716
信誉分 100
可用分 11045
阅读权限 255
注册 状态 离线
.我怎么记得我看的文献中都是去掉氢的呢。。。这个应该都是化学吸附吧，要是不脱氢不是物理吸附了么。涉及化学增强应该都是化学吸附的。当然要是纯粹的物理增强可能不需要，但是一般SERS中这两种增强都有吧。
2.实验上都是银纳米粒子，一般都是用银团簇模拟银纳米粒子。理论计算通常都是从最小的团簇模型开始尝试的，最简单的就是使用一个银。随着银原子数目不断增加，基于1点，还是脱掉一个氢的状态。
UID 129743
信誉分 100
可用分 8272
阅读权限 255
注册 状态 离线
谢谢了。晕，我刚看见你的回帖，谢谢了。我刚麻烦版主追加了2次金币，现在帖子又关闭了，不好意思再麻烦版主了，送你一朵鲜花吧，万分感谢！
UID 128422
信誉分 100
可用分 10762
阅读权限 255
注册 状态 离线
啊，金币啥的都是浮云。对于物理增强我也不太明白，希望对你有帮助吧。
UID 128223
信誉分 100
可用分 12135
阅读权限 255
注册 来自 内蒙古乌兰察布集宁状态 离线
谢谢你的帮助，你也是做拉曼计算的吗？
UID 128419
信誉分 100
可用分 12003
阅读权限 255
注册 状态 离线
我也想计算石墨烯的SERS，但是一直没有找到合适的方法，想请教你觉得知道声子谱的话可不可以预测一下
UID 128421
信誉分 102
可用分 9358
阅读权限 255
注册 状态 离线
不好意思啊，我完全不懂，我只看过单分子的计算。}