| 查看: 2589 | 回复: 18 | ||||
ys322铜虫 (正式写手)
|
[求助]
微扰理论硬球模型 已有1人参与
|
|||
| 请问一下 微扰理论硬球模型 ，分别怎么解释吖，小白不懂，求大神解释一下，看了半天理论也不懂，硬球半径，链节什么的 |
回复此楼» 猜你喜欢
0831生医工第一轮调剂失败求助
已经有17人回复
-
291求调剂
已经有11人回复
-
297，工科调剂?河南农业大学本科
已经有14人回复
-
0854求调剂
已经有20人回复
-
26药学专硕105500求调剂
已经有7人回复
-
一志愿华中农业071010，320求调剂
已经有15人回复
-
2026 WR青拔
已经有7人回复
-
一志愿中科大材料与化工，353分还有调剂学校吗
已经有9人回复
-
271求调剂
已经有37人回复
-
申博自荐
已经有3人回复
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
abinitio: 金币+5, 感谢精彩应助！ 2014-12-22 08:09:23
感谢参与，应助指数 +1
abinitio: 金币+5, 感谢精彩应助！ 2014-12-22 08:09:23
|
微扰理论简单说来就是你有一个reference state，对于reference state的热力学性质都已知，如果对reference state的哈密顿量做一个微小的扰动得到一个新的系统，那么新系统的热力学性质可以利用统计力学的理论从reference state近似推出。如果想获得动力学性质的话需要更复杂的微扰理论。 硬球模型(hard sphere）就是半径为r的球，当球与球之间距离大于2r时无相互作用，距离小于2r时作用为无穷大。硬球模型是最简单的描述相互作用的模型之一，其性质被研究得很透彻，状态方程也有解析表达式。所以经常被拿来做reference state。 比较经典的例子就是Weeks-Chandler-Andersen (WCA) Perturbation Theory http://scitation.aip.org/content ... 2/10.1063/1.1674820 这个实际上就是拿硬球模型作为reference state（虽然这里用的是WCA势，但WCA势和硬球模型差别不大），通过微扰理论可以求出Lennard-Jones流体的热力学性质。 |
2楼2014-12-18 23:42:19
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
6楼2014-12-19 09:31:16
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
【答案】应助回帖
★ ★ ★ ★ ★
abinitio: 金币+5, 感谢精彩应助！ 2014-12-22 08:09:33
abinitio: 金币+5, 感谢精彩应助！ 2014-12-22 08:09:33
|
每条链由若干个球组成，segment number就是球的数量，segment diameter 就是球的直径，segment energy就是球与球之间相互吸引的能阱深度(Eq. 1). 这几个参数不是直接体现在状态方程里的，而是在推导状态方程的过程中在不同的地方出现。segment number和segment diameter出现在硬球模型的compressibility和radial distribution function (Eq. 7 and 8) 里，因为硬球模型的参数和square well模型比其实就是少了个segment energy。 segment energy在这里基本是以微扰出现的，所以在Eq. 11 and 12中出现。 另外纠正一个我之前回答中的错误，我误以为Eq 4是配分函数，其实这里是compressibility。 |
11楼2014-12-20 11:31:51
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
13楼2014-12-21 12:25:16
ys322
铜虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 5325.3
- 红花: 1
- 帖子: 638
- 在线: 204.7小时
- 虫号: 916773
- 注册: 2009-11-30
- 性别: GG
- 专业: 化工系统工程
3楼2014-12-18 23:47:48
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
4楼2014-12-19 00:33:29
ys322
铜虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 5325.3
- 红花: 1
- 帖子: 638
- 在线: 204.7小时
- 虫号: 916773
- 注册: 2009-11-30
- 性别: GG
- 专业: 化工系统工程
5楼2014-12-19 09:10:54
ys322
铜虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 5325.3
- 红花: 1
- 帖子: 638
- 在线: 204.7小时
- 虫号: 916773
- 注册: 2009-11-30
- 性别: GG
- 专业: 化工系统工程
7楼2014-12-19 10:09:03
ys322
铜虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 5325.3
- 红花: 1
- 帖子: 638
- 在线: 204.7小时
- 虫号: 916773
- 注册: 2009-11-30
- 性别: GG
- 专业: 化工系统工程
8楼2014-12-19 10:38:03
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
【答案】应助回帖
★ ★ ★ ★ ★
abinitio: 金币+5, 感谢精彩应助！ 2014-12-22 08:09:44
abinitio: 金币+5, 感谢精彩应助！ 2014-12-22 08:09:44
|
我把文章前半部分粗略过了一遍，计算细节太繁琐，我也不是太懂，但大致可以理解这篇文章讲的是什么意思，在这里简单讲一下。 这篇文章的目的是为了求出链状分子的状态方程，状态方程就是液体的压力和温度还有密度（在这里用的是packing fraction)的关系。 微扰理论的基本思想是我已知一个reference state的状态方程，想求出另外一个state的状态方程，前提是另外一个state和reference state的哈密顿量相差不能太大，否则误差会比较大。在这里reference state就是hard chain，其状态方程已知（Eq 5）。那么为什么不能直接用hard chain的状态方程来表征体系呢？因为hard chain模型分子间只有排斥力，而我们都知道分子间是短程排斥，长程相吸引的，所以要对hard chain的状态方程做出修正来得到更符合实际情况的状态方程。 为了得到Eq 1 对应的体系（square well potential）的状态方程，我们需要得到体系的配分函数，也就是Eq 4. 跟hard chain的配分函数相比，square well体系的配分函数多了一项Z_disp. 因为配分函数和自由能有着一一对应的关系，所以我们就是要求出square well体系和hard chain体系自由能之差。 为了求出自由能差值，在这里使用了Barker-Henderson (BH)的微扰理论（ref 23，24）。BH微扰理论跟我在之前帖子提到的WCA微扰理论在一些技术处理上有差别（我个人觉得WCA更好点，这里不细说），但基本思想是差不多的。BH的处理方法就是Eq 10. Eq 10中的A1和A2两项都可以写成含有hard chain体系的径向分布函数的积分，而hard chain体系的径向分布函数我们是已知的（有解析形式）。 至于为什么A1和A2两项都可以写成含有hard chain体系的径向分布函数的积分。这是有缘由的。在液体物理里，一个很重要的结论是对于简单流体（比如Lennard－Jones流体），其结构更多的是由粒子间排斥力决定的，而吸引力对结构的贡献不大。也就是说对于Lennard－Jones流体，在高密度下，它的径向分布函数和具有相同半径的硬球模型（hard sphere，无排斥力）是非常相似的。那么对于这篇文献所提到的体系，也是类似的情况，在高密度下，我们可以认为square well体系的径向分布函数和hard chain体系非常相似。 知道了体系的自由能后，自由能对体积或者packing fraction求导数，就可以得到状态方程了。这里还进一步把状态方程简化成只含几个参数的方程。 液体的微扰理论在上个世纪六七十年代比较热门的研究方向，虽然现在用的人不多，但还是有不少物理意义的。如果想进一步了解，可以看一下这篇文章后面关于BH微扰的参考文献，或者参阅以下书籍 1. Statistical Mechanics by Donald McQuarrie, Chapter 14. 2. Theory of Simple Liquids by Hansen & McDonald, Chapter 5 以上仅仅是我快速扫了一眼的理解，有些细节可能讲的不对，可以指出。 |
9楼2014-12-19 11:41:31
ys322
铜虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 5325.3
- 红花: 1
- 帖子: 638
- 在线: 204.7小时
- 虫号: 916773
- 注册: 2009-11-30
- 性别: GG
- 专业: 化工系统工程
10楼2014-12-19 15:16:01
5