我正在研究磁化模拟并使用 Metropolis 方法编写算法。
我正在使用能量变化和玻尔兹曼分布,但是,我的问题是关于试验步骤。磁化的试验位置应该是完全随机的,还是应该是添加到旧磁化位置上的随机微分?
我的能量基本上取决于一个角度,所以,我的范围是 0 到 Pi,为了得到接近 50% 的拒绝,我在每一步中添加或减去大约 4 的 pi - 这太大了。如果我把它减少到相当小的东西,那么新旧能量概率之间的能量差异大约是 0.999 到 1.001。(这里我创建了一个随机的小角度,并从磁化矢量的当前角度中添加或减去它。)