| 查看: 1190 | 回复: 2 | ||
[求助]
关于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、总觉得小木虫上很少有关于分子模拟在汽液平衡这方面的交流，不知道队伍啥时候能壮大.... |
回复此楼» 猜你喜欢
基金系统什么内容也没有
已经有6人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
国自然面上复盘~欢迎讨论
已经有12人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
面上合作单位盖章
已经有5人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
申请删除本帖
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 大学物化相平衡问题
已经有4人回复
- Lammps计算油水界面，油相密度、平衡等问题求指导
已经有6人回复
- 相平衡 关于半透膜的问题
已经有10人回复
- 问一个关于聚合物体系平衡的问题
已经有25人回复
- NPT系综下判断是否平衡的标志
已经有12人回复
- 扩散系数的计算及体系平衡的判断
已经有4人回复
- Aspen Plus 软件中气液平衡问题
已经有9人回复
- 径向分布函数的问题
已经有13人回复
- 使用Towhee做NVT模拟时，气液两盒子的hmatrix 大小该怎么设置？ 有什么原则或公式么？
已经有7人回复
- 关于气液平衡的简单问题
已经有24人回复
- 相平衡部分，小问题求助
已经有4人回复
- 关于towhee 回答zhang_jaj的问题
已经有44人回复
- topol文件不识别小分子
已经有8人回复
- MS金属氧化物表面吸附有机物平衡问题
已经有4人回复
- 用热力学计算的平衡转化率只有实验测得的反应转化率的一半。。计算是否会有问题
已经有5人回复
- 使用towhee进行GEMC模拟相平衡的MC运动定义问题..
已经有15人回复
- 【求助】关于相平衡数据的资料
已经有4人回复
- 【求助】气液平衡的数据
已经有5人回复
2楼2014-07-23 09:28:37
草田露
银虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 227.1
- 散金: 24
- 帖子: 26
- 在线: 60.3小时
- 虫号: 3080499
- 注册: 2014-03-22
- 性别: MM
- 专业: 化工热力学和基础数据
3楼2015-09-07 09:04:39
