24小时热门版块排行榜    

查看: 532  |  回复: 1

chenjingzhi

金虫 (小有名气)

[交流] monte carlo模拟 lattice gas

最近在做这个简单模型的模拟,但是发现有些问题,一直得不到正确结果。lattice gas模型要求使用grand canonical系综,因此依据《understanding of molecular simulation》这本书的伪代码,首先是判断是删除还是插入粒子,删除和插入的概率均为50%,然后随机选取粒子或者位置,进行粒子的删除或插入尝试移动,判断是否接受粒子的删除和插入,对于一般系统来说这一步无论什么情况下都是可行的,但是对于格点模型,在我们确定要进行删除和插入尝试移动之后,在我们随机选取格点之后,我们无法保证所选格点一定有粒子或无粒子,这也就导致了删除或插入尝试移动无法进行,如果我们对格点是否被粒子占有进行判断,那么我们也就无形之中改变了进行删除和插入尝试移动的概率,导致细致平衡与实际结果出现偏差。有没有人知道这个问题怎么解决?
       刚刚有个想法,是不是在删除粒子尝试移动里,我们应该随机抽取的是存在的粒子,而不是格点;而尝试插入时,应该随机选取没有粒子占有的格点?
回复此楼
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
本帖仅楼主可见
2楼2016-03-08 05:23:36
已阅   申请模拟EPI   编辑   查看我的主页
相关版块跳转 我要订阅楼主 chenjingzhi 的主题更新
信息提示
请填处理意见