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

登录注册

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

24小时热门版块排行榜

>基金申请 (7214)
>论坛更新日志 (5364)
>休闲灌水 (4578)
>虫友互识 (3253)
>健康生活 (1647)
>博后之家 (1225)
>导师招生 (649)
>信息科学 (638)
>招聘信息布告栏 (504)
>硕博家园 (435)
>文献求助 (391)
>论文投稿 (353)
>考博 (266)
>找工作 (218)
>教师之家 (192)
>有奖起名 (66)

返回列表

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

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

yahoohoo

铁杆木虫 (著名写手)

模拟EPI: 10
应助: 0 (幼儿园)
贵宾: 1.55
金币: 7636.4
散金: 251
红花: 4
帖子: 1179
在线: 167.5小时
虫号: 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

» 猜你喜欢

如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》已经有7人回复
没有任何消息-是不是就凉了已经有8人回复
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有3人回复
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急已经有5人回复
人气不行了已经有10人回复
我面上完蛋了已经有10人回复
能否退出参与的面上项目解除限项已经有24人回复
filecode，4个jtjc了已经有17人回复
2026国自然函评费到账已经有19人回复
建议基金发布提前给出明确的时间点已经有16人回复

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

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的回帖

yahoohoo

铁杆木虫 (著名写手)

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

★ ★ ★
lei0736(金币+3,VIP+0):谢谢 12-22 21:44

引用回帖:

Originally posted by lei0736 at 2009-12-22 10:01:
在MD模拟里有平衡和非平衡的MD
那么有没有平衡MC和非平衡MC之分？
他们之间有哪些区别与联系？是否非平衡MC不需要满足细致平衡条件？
如果有区别他们各自应用的领域特点又是什么？

没听说过非平衡MC，所以我不确定有没有。不过有所谓的kinetic MC，有人用来模拟硬球体系的结晶成核速率。至于是否需要满足细致平衡条件，我想应该还是要满足的。

赞一下(11人)

回复此楼

7楼2009-12-22 20:52:24

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

查看全部 74 个回答

lei0736

荣誉版主 (职业作家)

模拟EPI: 5
应助: 17 (小学生)
贵宾: 10.32
金币: 19244
散金: 5577
红花: 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的回帖

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
金币: 7636.4
散金: 251
红花: 4
帖子: 1179
在线: 167.5小时
虫号: 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 个回答

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

[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3	0XLacIJUOj8D 2026-08-24	5/250	2026-08-24 21:51 by BZKMTicpDhFj
[基金申请] 我面上完蛋了 +9	且听虎啸 2026-08-20	10/500	2026-08-24 20:08 by 苦寒来香
[基金申请] 能否退出参与的面上项目解除限项 +21	koalala 2026-08-24	24/1200	2026-08-24 19:25 by 家与远方
[基金申请] filecode，4个jtjc了 +14	ziyangfang 2026-08-19	17/850	2026-08-24 18:37 by 哈哈蛤？
[基金申请] 估计是周四 +4	archvillain 2026-08-18	4/200	2026-08-24 13:53 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6	炎黄贵胄 2026-08-22	7/350	2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7	wangzilk 2026-08-18	9/450	2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧！ +13	大萍1987 2026-08-20	16/800	2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4	snowwithsea 2026-08-19	17/850	2026-08-24 10:20 by echo8914667
[基金申请] 科研孤儿太难了 +18	我4大白菜 2026-08-20	19/950	2026-08-24 09:47 by 凯恩广盛催化分�
[基金申请] 什么时候开奖？ +10	CrisMessi 2026-08-18	11/550	2026-08-24 06:50 by 开心的小狮子
[教师之家] 跳槽后在研项目怎么办？ +5	简单化xn 2026-08-22	10/500	2026-08-23 12:38 by 简单化xn
[基金申请] 今天基金会出结果吗？20260819 +16	kkkl_v 2026-08-19	17/850	2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天，20号变了 +5	archvillain 2026-08-20	5/250	2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了？ +8	chengyan1220 2026-08-21	11/550	2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13	wuchongjun 2026-08-20	19/950	2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4	哈哈蛤？ 2026-08-21	4/200	2026-08-21 10:58 by Vivilian
[基金申请] 重要消息，中午系统在维护 +11	yuleib84 2026-08-18	12/600	2026-08-20 11:09 by xskun
[基金申请] 明天放榜？ +5	Shxjjxjkx 2026-08-18	5/250	2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护祝所有人高中 +8	gjjjzhong 2026-08-18	9/450	2026-08-18 13:01 by 家与远方

24小时热门版块排行榜

yahoohoo

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

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

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

» 猜你喜欢

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

yahoohoo

lei0736

qzhaosdu

yahoohoo

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