| 查看: 542 | 回复: 1 | |||
chenjingzhi金虫 (小有名气)
|
[交流]
monte carlo模拟 lattice gas
|
|
最近在做这个简单模型的模拟，但是发现有些问题，一直得不到正确结果。lattice gas模型要求使用grand canonical系综，因此依据《understanding of molecular simulation》这本书的伪代码，首先是判断是删除还是插入粒子，删除和插入的概率均为50%，然后随机选取粒子或者位置，进行粒子的删除或插入尝试移动，判断是否接受粒子的删除和插入，对于一般系统来说这一步无论什么情况下都是可行的，但是对于格点模型，在我们确定要进行删除和插入尝试移动之后，在我们随机选取格点之后，我们无法保证所选格点一定有粒子或无粒子，这也就导致了删除或插入尝试移动无法进行，如果我们对格点是否被粒子占有进行判断，那么我们也就无形之中改变了进行删除和插入尝试移动的概率，导致细致平衡与实际结果出现偏差。有没有人知道这个问题怎么解决？ 刚刚有个想法，是不是在删除粒子尝试移动里，我们应该随机抽取的是存在的粒子，而不是格点；而尝试插入时，应该随机选取没有粒子占有的格点? |
回复此楼» 猜你喜欢
退学或坚持读
已经有26人回复
-
面上基金申报没有其他的参与者成吗
已经有5人回复
-
有70后还继续奋斗在职场上的吗？
已经有5人回复
-
遇见不省心的家人很难过
已经有22人回复
-
博士延得我，科研能力直往上蹿
已经有5人回复
-
多组分精馏求助
已经有6人回复
 本帖仅楼主可见 |
2楼2016-03-08 05:23:36