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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » Monte Carlo » 关于GEMC模拟相平衡的一些问题。
31/1返回列表
查看: 1003  |  回复: 2
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

992648352

银虫 (小有名气)

[求助] 关于GEMC模拟相平衡的一些问题。

我使用的是MCCCS Towhee来模拟纯物质的共存曲线，每次模拟还是碰到很多类似问题，现在想就我碰到的问题请教一下大家：
1、体积的设定：两相体积的设定，文献Sample Text给的建议是总体积的设定最好保证20%的粒子数目在气相。不知道大家是怎么看这个建议的。在低温的时候20%的粒子在气相会造成气相的体积非常之大，有的时候汽液相体积之比在10以上了这样做会不会影响模拟？
2、各种概率的设置：很多文献（Ilja Siepm ann的系列文章中都有提到）上写的是体积和粒子交换（swap）这两个移动概率要保证在5~10个cycles之内分别有一次被接受；
      我是这样做的，体积概率 P 的设置：N（分子数）*nsteps（循环数）*0.5（接受概率）*P（体积变化的概率）=nsteps（循环数）/（5~10）；
       swap的概率这个起初我是随便给的0.3左右，先跑个10000cycles再看看接受的swap数目{两部分（box1到box2+box2到box1）}的大小，再根据情况调试，比如：10000cycles有3000次接受。说明swap概率大了，需要适当调小；500次被接受说明需要适当增大；剩下的概率均分给其他的移动。但是，问题出现在低温的情况，低温的swap就算给到了65%也满足不了5~10个cycles之内分别有一次被接受，并且有些文献有说，除去体积和swap的概率，剩下的概率要不要小于50%。不知道大家是怎么分配的。
3、模拟的平衡与否的判断条件：看文献写的是有明显的不变的汽液密度，说明达到了平衡。但是我模拟的时候汽液密度一直是波动的，而且波动还是挺大的，我想是不是我的步数不够，就一直接着跑，但是看起过来他就是达不到平衡（收敛）。文献上都说的一般都是10000~25000达到平衡，接着另外20000~40000产出。我都快10w步还没看到平衡。。。。。。一开始接触MC时候，不是说看体系的总能量减小吗，Metropolis抽样法不是朝着能量小的构型去构建马尔科夫链的吗?为什么结果总是看到总能量（box1+box2）在不断的减小又增大（有时差一个数量级）。
4、总觉得小木虫上很少有关于分子模拟在汽液平衡这方面的交流，不知道队伍啥时候能壮大....
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2014-07-22 21:46:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

992648352

银虫 (小有名气)
有没有觉得气相的平动和转动的最大位移经常达到最大，是的接受概率几乎达到100%与想要的50%相差甚远。。。
赞一下 回复此楼
2楼2014-07-23 09:28:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

草田露

银虫 (初入文坛)
帮顶！！！！
回复此楼
3楼2015-09-07 09:04:39
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 992648352 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭