| 查看: 555 | 回复: 1 | |||
chenjingzhi金虫 (小有名气)
|
[交流]
monte carlo模拟 lattice gas
|
|
最近在做这个简单模型的模拟，但是发现有些问题，一直得不到正确结果。lattice gas模型要求使用grand canonical系综，因此依据《understanding of molecular simulation》这本书的伪代码，首先是判断是删除还是插入粒子，删除和插入的概率均为50%，然后随机选取粒子或者位置，进行粒子的删除或插入尝试移动，判断是否接受粒子的删除和插入，对于一般系统来说这一步无论什么情况下都是可行的，但是对于格点模型，在我们确定要进行删除和插入尝试移动之后，在我们随机选取格点之后，我们无法保证所选格点一定有粒子或无粒子，这也就导致了删除或插入尝试移动无法进行，如果我们对格点是否被粒子占有进行判断，那么我们也就无形之中改变了进行删除和插入尝试移动的概率，导致细致平衡与实际结果出现偏差。有没有人知道这个问题怎么解决？ 刚刚有个想法，是不是在删除粒子尝试移动里，我们应该随机抽取的是存在的粒子，而不是格点；而尝试插入时，应该随机选取没有粒子占有的格点? |
回复此楼» 猜你喜欢
食品专硕 一志愿双一流 328
已经有4人回复
-
材料调剂
已经有6人回复
-
300分，材料，求调剂，英一数二
已经有5人回复
-
求调剂，一志愿:南京航空航天大学大学 ，080500材料科学与工程学硕，总分289分
已经有5人回复
-
085600材料与化工调剂
已经有13人回复
-
材料与化工考研调剂
已经有7人回复
-
0856材料专硕353求调剂
已经有6人回复
-
07化学280分求调剂
已经有6人回复
-
086003食品工程求调剂
已经有5人回复
-
08工学调剂
已经有16人回复
 本帖仅楼主可见 |
2楼2016-03-08 05:23:36
40