| 查看: 652 | 回复: 1 | |||
chenjingzhi金虫 (小有名气)
|
[交流]
monte carlo模拟 lattice gas
|
|
最近在做这个简单模型的模拟，但是发现有些问题，一直得不到正确结果。lattice gas模型要求使用grand canonical系综，因此依据《understanding of molecular simulation》这本书的伪代码，首先是判断是删除还是插入粒子，删除和插入的概率均为50%，然后随机选取粒子或者位置，进行粒子的删除或插入尝试移动，判断是否接受粒子的删除和插入，对于一般系统来说这一步无论什么情况下都是可行的，但是对于格点模型，在我们确定要进行删除和插入尝试移动之后，在我们随机选取格点之后，我们无法保证所选格点一定有粒子或无粒子，这也就导致了删除或插入尝试移动无法进行，如果我们对格点是否被粒子占有进行判断，那么我们也就无形之中改变了进行删除和插入尝试移动的概率，导致细致平衡与实际结果出现偏差。有没有人知道这个问题怎么解决？ 刚刚有个想法，是不是在删除粒子尝试移动里，我们应该随机抽取的是存在的粒子，而不是格点；而尝试插入时，应该随机选取没有粒子占有的格点? |
回复此楼» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有10人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有9人回复
-
厉娜
已经有6人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有9人回复
 本帖仅楼主可见 |
2楼2016-03-08 05:23:36
10