版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

本帖产生 1 个模拟EPI ，点击这里进行查看

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yahoohoo

铁杆木虫 (著名写手)

模拟EPI: 10
应助: 0 (幼儿园)
贵宾: 1.55
金币: 7632.4
散金: 251
红花: 4
帖子: 1176
在线: 167小时
虫号: 74894
注册: 2005-06-15
性别: GG
专业: 理论和计算化学

[交流] 【在线答疑】经典粒子体系的Monte Carlo 模拟之基础篇已有29人参与

应版主lei0736之盛情相邀，特开此贴，虽名为分享，实为讨论，希望与初学者以及从事MC研究的同行一起交流。第一篇关于粒子体系的MC模拟基础，欲将涉及以下几个概念与问题：

概念:
正则系综、相空间、概率密度、细致平衡(Detailed Balance)、Metropolis重要性抽样

问题:
(1) 为什么用MC？
(2) MC能做什么？不能做什么？
(3) 怎么做MC？

讨论中打算在适当之处添加一些伪代码(c++)以供理解。

欢迎大家多提宝贵意见。

[ Last edited by lei0736 on 2009-12-21 at 18:04 ]

回复此楼

» 本帖已获得的红花（最新10朵）

zhang_jaj

» 猜你喜欢

1楼2009-12-21 16:53:16

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

tpp001

金虫 (著名写手)

应助: 0 (幼儿园)
金币: 759.3
散金: 5
红花: 1
帖子: 1039
在线: 155.9小时
虫号: 423029
注册: 2007-07-21
性别: GG
专业: 化学

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

我看的不少文献中有这样的模拟方法,用widom insertion techinique 来进行mc模拟计算溶解度,不知道楼主能不能给介绍啊，因为我下一步可能就用这个来计算溶解度!谢谢!!!另祝楼主新年快乐

赞一下(1人)

回复此楼

迷茫在知识的海洋里，需要你的指导。thankyou

22楼2010-01-04 20:32:16

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 74 个回答

lei0736

荣誉版主 (职业作家)

模拟EPI: 5
应助: 17 (小学生)
贵宾: 10.32
金币: 20552.5
散金: 5527
红花: 30
沙发: 1
帖子: 4762
在线: 1077.4小时
虫号: 62582
注册: 2005-03-17
性别: GG
专业: 凝聚态物性I:结构、力学和
管辖: 量子化学

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

祝贺开贴
顺便问下细致平衡条件的重要性如何
以及如何在模拟中注意哪些事项以保证细致平衡条件

赞一下(1人)

回复此楼

2楼2009-12-21 18:55:47

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

qzhaosdu

金虫 (著名写手)

应助: 0 (幼儿园)
金币: 1401.6
散金: 2052
红花: 7
帖子: 1969
在线: 263.5小时
虫号: 493351
注册: 2008-01-10
性别: GG
专业: 理论和计算化学
管辖: 量子化学

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+1):欢迎参与讨论 2010-01-30 17:00

请问用MC做表面活性剂自组装和高分子稀溶液自组装的可靠性？
看过很多文章，貌似用MC做同样的体系和用SCFT,DPD做出来的结论大不相同。谢谢！

赞一下(2人)

回复此楼

Anewday,anewhour,anewminute,anewpeople.

3楼2009-12-21 21:13:07

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

yahoohoo

铁杆木虫 (著名写手)

模拟EPI: 10
应助: 0 (幼儿园)
贵宾: 1.55
金币: 7632.4
散金: 251
红花: 4
帖子: 1176
在线: 167小时
虫号: 74894
注册: 2005-06-15
性别: GG
专业: 理论和计算化学

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lei0736(金币+10,VIP+0):呵呵值得推荐 12-22 09:57
zeoliters(金币+5,VIP+0):谢谢！ 12-25 16:25

引用回帖:

Originally posted by lei0736 at 2009-12-21 18:55:
祝贺开贴
顺便问下细致平衡条件的重要性如何
以及如何在模拟中注意哪些事项以保证细致平衡条件

Detailed balance 用公式表示为
$P(\mathrm{o}) \alpha(\mathrm{o}\to\mathrm{n}) P_\mathrm{acc}(\mathrm{o}\to\mathrm{n})= P(\mathrm{n}) \alpha(\mathrm{n}\to\mathrm{o}) P_\mathrm{acc}(\mathrm{n}\to\mathrm{o})$
其中， $P(\mathrm{o})$ 表示微观态或微观构象 $\mathrm{o}$ 的概率(或几率)密度分布，在NVT系综中 $P(\mathrm{o})=\frac{e^{-\beta U(\mathrm{o})}}{\int\mathrm{d}\Gamma e^{-\beta U(\Gamma)}}$ ，即Boltzmann分布。 $\alpha(\mathrm{o}\to\mathrm{n})$ 表示从构象 $\mathrm{o}$ 到构象 $\mathrm{n}$ 所设计进行的尝试运动的几率，而 $P_\mathrm{acc}(\mathrm{o}\to\mathrm{n})$ 则是接受从构象 $\mathrm{o}$ 到构象 $\mathrm{n}$ 的尝试运动的几率，因此 $\alpha(\mathrm{o}\to\mathrm{n}) P_\mathrm{acc}(\mathrm{o}\to\mathrm{n})$ 表示从构象 $\mathrm{o}$ 到构象 $\mathrm{n}$ 的转移几率。用简单的语言描述，Detailed balance 是指从构象 $\mathrm{o}$ 到构象 $\mathrm{n}$ 的流量(flow)等于从构象 $\mathrm{n}$ 到构象 $\mathrm{o}$ 的流量，这也就是所谓的微观可逆性。

Detailed balance 的重要性在于 (1) 保证在有限的MC步数内达到所需的平衡概率分布[c.f. M.E.J. Newman & G.T. Barkema, Monte Carlo Methods in Statistical Physics, Oxford]，这一点有严格的数学证明，我们只需要知道这个结论就可以了； (2) 达到平衡后平衡就不再被破坏，这意味着根据MC算法产生的构象将满足平衡概率分布。

满足Detailed balance 公式的 $P_\mathrm{acc}(\mathrm{o}\to\mathrm{n})$ 解有很多，使用最为广泛的当属Metropolis 算法：

$P_\mathrm{acc}(\mathrm{o}\to\mathrm{n})=\mathrm{min}\left\{1.0, \frac{P(\mathrm{n}) \alpha(\mathrm{n}\to\mathrm{o})}{P(\mathrm{o}) \alpha(\mathrm{o}\to\mathrm{n})} \right\}$
这里， $P(\mathrm{o})/P(\mathrm{n})$ 是已知的，对于NVT系综，
$\frac{P(\mathrm{o})}{P(\mathrm{n})}=e^{-\beta\bigl\{U(\mathrm{n})-U(\mathrm{o})\bigr\}}$ 。而 $\alpha(\mathrm{o}\to\mathrm{n})$ ( $\alpha(\mathrm{n}\to\mathrm{o})$ )取决于我们设计的Monte Carlo运动(move)是如何由构象 $\mathrm{o}$ ( $\mathrm{n}$ )到构象 $\mathrm{n}$ ( $\mathrm{o}$ )的。例如对于简单的粒子体系，我们现在设计随机移动一个粒子的位置，伪代码为

int ichosen = int(Ntot * mt19937.getReal3()); // Ntot: 粒子的数目; mt19937.getReal3()返回(0, 1)之间的随机数
double stepsize = 0.5 * sigma_LJ; //假设LJ 粒子体系,  sigma_LJ为LJ直径. 在实际模拟中一般需要调节stepsize的大小以获得一定的尝试运动成功率, 例如30%.
for (int j = 0; j < 3; ++j)
  xn[ichosen][j] = xo[ichosen][j] + stepsize * (mt19937.getReal3()-0.5); //(2.0*mt19937.getReal3()-1.0)返回(-0.5,0.5)的随机值.

这段代码的for循环实现的是以随机选取的粒子的原有位置xo[ichosen]为中心所形成的边长为 stepsize的正方体内随机产生一个新的位置xn[ichosen]，那么， $\alpha(\mathrm{o}\to\mathrm{n}) \propto \mathrm{stepsize}^{-3}$ ，因为stepsize * (mt19937.getReal3()-0.5) 保证了由旧位置到新位置的尝试运动几率等于由新位置到旧位置的尝试运动几率，所以这里 $\alpha(\mathrm{n}\to\mathrm{o})=\alpha(\mathrm{o}\to\mathrm{n})$ 。有了这样设计的尝试运动再加上利用 Metropolis 算法 $P_\mathrm{acc}(\mathrm{o}\to\mathrm{n})=\mathrm{min}\left\{1.0, P(\mathrm{n})/ P(\mathrm{o})\right\}$ 来判断是否接受这样的尝试运动，达到平衡后的构象分布就已经满足我们想要抽样的统计概率 $P(\Gamma)$ 分布了。

简单一句话， 通过Detailed balance设计MC尝试运动，保证我们能达到平衡而且平衡后的构象满足我们想要的概率分布。

仍然以NVT系综为例，任意可测量量 $\mathcal{O}$ 的统计定义为
$\langle \mathcal{O} \rangle_\mathrm{NVT} = \frac{\int\mathrm{d}\Gamma\mathcal{O}(\Gamma) e^{-\beta U(\Gamma)}}{\int\mathrm{d}\Gamma e^{-\beta U(\Gamma)}} =\sum\limits_{i=1}^{N_\mathrm{sample}} P_\mathrm{NVT}(\Gamma_i)\mathcal{O}(\Gamma_i)$
需要注意的是：我们已经知道  达到平衡后的构象分布就已经满足我们想要抽样的统计概率 $P(\Gamma)$ 分布了，这就表明在MC中，
$\langle \mathcal{O} \rangle_\mathrm{NVT} = \sum\limits_{i=1}^{N_\mathrm{sample}} \mathcal{O}(\Gamma_i)$

因为MC所产生的平衡构象 $\Gamma_i$ 已经满足平衡分布 $P_\mathrm{NVT}(\Gamma_i)$ 了。

[ Last edited by yahoohoo on 2009-12-22 at 04:21 ]

赞一下(20人)

回复此楼

4楼2009-12-21 22:28:15

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 74 个回答

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[论文投稿] 你的学校职称晋升认CSCD吗？ +3	wowsunflower 2024-09-26	3/150	2024-09-29 03:20 by lyd782
[论文投稿] applied catalysis B 投稿模版 5+4	午后夏日 2024-09-27	5/250	2024-09-28 22:59 by 莱茵润色
[基金申请] QB怎么还没消息 +21	wuxiacl 2024-09-24	57/2850	2024-09-28 14:48 by xiaoxiao270
[考博] 25有机合成申博 +3	再等你啊 2024-09-24	4/200	2024-09-28 13:09 by majun521
[基金申请] 山东省基金公示结束 +6	ishow19 2024-09-24	14/700	2024-09-28 12:20 by ishow19
[电化学] 刚装好的电池电压掉很快怎么解决？？ +3	花花你最棒 2024-09-26	6/300	2024-09-28 08:20 by cnofc
[教师之家] 压力山大 +17	暮色恋伊人 2024-09-24	25/1250	2024-09-27 16:10 by 领先于123
[硕博家园] 当前读博士还有性价比吗？ +8	苏东坡二世 2024-09-23	8/400	2024-09-27 11:39 by 半生梦君
[考博] 数学博导 +4	学术霸王 2024-09-25	6/300	2024-09-27 11:14 by 青古
[有机交流] 二氯甲烷的去除 +3	cgsa吧 2024-09-24	8/400	2024-09-27 10:55 by bear2007
[论文投稿] 期刊论文发表了但查询到sci未收录 50+5	yibuxiao 2024-09-25	11/550	2024-09-27 10:24 by 莱茵润色
[论文投稿] 为什么审稿人一直拒绝审稿啊？ +7	潇湘雪月 2024-09-24	13/650	2024-09-27 07:53 by weizhi111
[考博] 有机转材料申博难度 30+3	昕散 2024-09-22	10/500	2024-09-26 22:42 by 824282658
[有机交流] 请问胺的盐酸盐中氯化氢的氢会在核磁氢谱中出峰吗？ +3	rommel1975 2024-09-25	4/200	2024-09-26 13:07 by 091602
[基金申请] 请教2024后期资助到哪个阶段了 +3	拾光者5566 2024-09-23	8/400	2024-09-26 09:27 by 拾光者5566
[论文投稿] NC投稿多久有消息呀？ +4	Jordanblood 2024-09-24	6/300	2024-09-25 16:58 by FZX_2024
[论文投稿] laser physics期刊投稿 5+3	mengxiangcz 2024-09-23	6/300	2024-09-25 15:08 by mengxiangcz
[论文投稿] SCI投稿状态 25+4	jorden8 2024-09-23	6/300	2024-09-25 10:24 by jorden8
[论文投稿] 爱斯维尔旗下Postharvest Biology and Technology，采用什么系统查查 20+4	sdsdsxh 2024-09-22	9/450	2024-09-25 07:33 by sdsdsxh
[有机交流] 做聚酰亚胺一步法合成，机械搅拌怎么密封不漏气 +3	哼哈EN 2024-09-24	5/250	2024-09-24 14:16 by mrzhl1986