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

登录注册

应《网络安全法》要求，自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

» 猜你喜欢

366求调剂已经有11人回复
290求调剂已经有23人回复
297，工科调剂? 已经有5人回复
290调剂生物0860 已经有33人回复
一志愿鲁东大学071000生物学学硕初试分数276求调剂已经有28人回复
279学硕食品专业求调剂院校已经有28人回复
071000生物学调剂求助已经有21人回复
279求调剂已经有12人回复
调剂已经有12人回复
271求调剂已经有41人回复

» 本主题相关价值贴推荐，对您同样有帮助:

Monte carlo模拟已经有5人回复
Monte Carlo 模拟时图是用什么软件做的已经有10人回复
【求助】是不是monte carlo 模拟只适用于高分子方向的？已经有6人回复
【求助】求如何学习monte carlo仿真已经有7人回复
【求助】用monte carlo法做分子模拟，用什么编程语言实现比较好？已经有25人回复
【资源】Introduction to Monte Carlo Method 已经有7人回复
【讨论】大家对未来GPU在monte carlo中的应用持有什么样的态度(CUDA)？已经有15人回复
【转帖】一道面试题的monte carlo实现已经有8人回复
【求助】想用Monte Carlo方法预测晶体磁性转变温度应该用什么程序已经有5人回复
【求助】有用Monte Carlo 编过有关Fluid中粒子运动的程序已经有6人回复
【求助】求Monte Carlo 关于扩散系数计算的源代码已经有5人回复

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

已阅回复此楼关注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的回帖

查看全部 74 个回答

lei0736

荣誉版主 (职业作家)

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

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

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

赞一下(1人)

回复此楼

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

已阅回复此楼关注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的回帖

yahoohoo

铁杆木虫 (著名写手)

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

★ ★ ★ ★ ★
lei0736(金币+5,VIP+0):谢谢 12-22 09:58

引用回帖:

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

不知你所说的可靠性是什么？如果采用的模型合理、算法以及程序没有问题，那得到的结果就应该是可靠的。

至于跟其他方法比较的结论大不相同，我没有看过很多这方面的文章，不过需要注意的是 SCFT一般都是平均场处理，因此涨落(fluctuation)是被忽略的，而在相分离过程中这是不对的，Marcus Mueller就是做SCFT和MC对比的，我没有看到你所谓的大不相同的结论。而DPD虽然是一种很流行的研究自组装的方法，但是它所用的相互作用参数过于简单，是否丢掉了部分物理内涵也不清楚。譬如说，DPD就不能模拟高分子融体的缠结问题。

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

赞一下(11人)

回复此楼

5楼2009-12-21 22:48:57

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

查看全部 74 个回答

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

[考研] 297，工科调剂? +5	河南农业大学-能 2026-04-14	5/250	2026-04-14 15:48 by lsb1022
[教师之家] 转长聘了 +7	简单化xn 2026-04-13	7/350	2026-04-14 14:50 by xindong
[考研] 求调剂 +16	雪逢冬 2026-04-10	16/800	2026-04-14 14:27 by 逆水乘风
[考研] 274求调剂 +11	山阿蔓 2026-04-07	11/550	2026-04-14 14:21 by 不我拉绿卡
[考研] 297求调剂 +23	ORCHID1 2026-04-10	26/1300	2026-04-14 13:52 by 陈皮皮
[考研] 药学305求调剂 +4	玛卡巴卡boom 2026-04-11	4/200	2026-04-14 10:21 by dick_runner
[考研] 332求调剂 +15	蕉蕉123 2026-04-10	15/750	2026-04-13 23:12 by pies112
[考研] 求调剂，985材料与化工348分 +9	涵竹刘 2026-04-11	14/700	2026-04-13 22:26 by 涵竹刘
[考研] 一志愿211 0703化学 346分求调剂 +26	土豆er? 2026-04-09	29/1450	2026-04-13 15:15 by 独醉梦孤城
[考研] 322求调剂，08工科 +4	今天是个小号 2026-04-08	4/200	2026-04-13 00:20 by baobaoye
[考研] 。 +10	李多米lee. 2026-04-12	11/550	2026-04-12 22:58 by yuyin1233
[考研] 344 材料专业求调剂211 无地域要求 +8	hualkop 2026-04-11	8/400	2026-04-12 22:24 by fqwang
[考研] 22408 352分求调剂 +5	努力的夏末 2026-04-09	5/250	2026-04-12 19:17 by wj165256
[考研] 291求调剂 +11	关忆北. 2026-04-09	12/600	2026-04-12 10:32 by 逆水乘风
[考研] 材料与化工300求调剂 +39	肖开文 2026-04-09	43/2150	2026-04-12 01:30 by 秋豆菜芽
[考研] 087100初试311求调剂 +4	任雅琴 2026-04-09	4/200	2026-04-11 10:33 by zhq0425
[考研] 化学工程与技术324调剂 +23	孙常华 2026-04-09	25/1250	2026-04-11 00:07 by 骑牛渡寒江
[考研] 368求调剂 +3	17385968172 2026-04-10	3/150	2026-04-10 20:12 by 电子奥德彪
[考研] 机械专368 有去处吗 +4	种大树 2026-04-10	4/200	2026-04-10 15:31 by jiajinhpu
[考研] 青岛科技大学材料学院，环境学院调剂补录4月10日以前都可以 +3	1青科大。 2026-04-09	5/250	2026-04-10 09:58 by 翩翩一书生

24小时热门版块排行榜

yahoohoo

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

» 收录本帖的淘帖专辑推荐

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

qzhaosdu

lei0736

yahoohoo

yahoohoo

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