s***e
发帖数: 911 | 1 我摹拟的是一个链状分子的统计物理特性. 这个分子由系列小球相连而成. 给定
这些小球的初始位置, 我计算相应的弹性能(目前只是弯曲和拉伸). 设若相邻
小球平衡间距为1, 那么什么样的微扰幅度最优?
目前运行比较良好的微扰是这样选的: 0.05*Random[0,1]. 但是这个选择下被接受的
微扰构型很少, 程序运行时间太长; 若选择0.5*Random[0,1], 则干脆就没有被接受的
微扰构型.
请教对Metropolis 算法熟悉的同志指点一下. THX... | s*a
发帖数: 33 | 2 I assume this is a simulation using classical mechanics (i.e., not quantum
systems). Ususally one choose the acceptance ratio of the random move
to be around 50%, (in practice, 10% to 90% is acceptable)
so that the simulation can sample the phase space more
efficiently. You can change this ratio by choosing the magnitude of your
perturbance (the 0.05 or 0.5 number in your example). If even in 0.05
(the simulation time you feel already too long) and the acceptance ratio
still too low, then there
【在 s***e 的大作中提到】
: 我摹拟的是一个链状分子的统计物理特性. 这个分子由系列小球相连而成. 给定
: 这些小球的初始位置, 我计算相应的弹性能(目前只是弯曲和拉伸). 设若相邻
: 小球平衡间距为1, 那么什么样的微扰幅度最优?
: 目前运行比较良好的微扰是这样选的: 0.05*Random[0,1]. 但是这个选择下被接受的
: 微扰构型很少, 程序运行时间太长; 若选择0.5*Random[0,1], 则干脆就没有被接受的
: 微扰构型.
: 请教对Metropolis 算法熟悉的同志指点一下. THX...
| s*a
发帖数: 33 | 3
For a single chain, this might not help too much, I think.
There is a "multi-level" algorithm for polymers:
Say a segment of your chain is like this:
...abcde.... (five balls here: a, b, ... and e)
Your current algorithm is to fix the position of "a" and "c" (and
all other balls), then
try to move "b"; then fix "b" and "d", then try to move "c", right?
Take the first step for example, since "a" and "c" are fixed, so "b"
can't move too far, otherwise the chain will be streched too much.
In a mul
【在 s***e 的大作中提到】
: 我摹拟的是一个链状分子的统计物理特性. 这个分子由系列小球相连而成. 给定
: 这些小球的初始位置, 我计算相应的弹性能(目前只是弯曲和拉伸). 设若相邻
: 小球平衡间距为1, 那么什么样的微扰幅度最优?
: 目前运行比较良好的微扰是这样选的: 0.05*Random[0,1]. 但是这个选择下被接受的
: 微扰构型很少, 程序运行时间太长; 若选择0.5*Random[0,1], 则干脆就没有被接受的
: 微扰构型.
: 请教对Metropolis 算法熟悉的同志指点一下. THX...
| E*****y
发帖数: 33 | 4 accept rate 0.3左右比较好. 你可以模拟退火过程, 从高温逐渐降到你所要
的温度, 这样或许比你直接模拟你想要的温度达到热平衡快...
【在 s***e 的大作中提到】
: 我摹拟的是一个链状分子的统计物理特性. 这个分子由系列小球相连而成. 给定
: 这些小球的初始位置, 我计算相应的弹性能(目前只是弯曲和拉伸). 设若相邻
: 小球平衡间距为1, 那么什么样的微扰幅度最优?
: 目前运行比较良好的微扰是这样选的: 0.05*Random[0,1]. 但是这个选择下被接受的
: 微扰构型很少, 程序运行时间太长; 若选择0.5*Random[0,1], 则干脆就没有被接受的
: 微扰构型.
: 请教对Metropolis 算法熟悉的同志指点一下. THX...
|
|
