Modelowanie molekularne - metody Monte Carlo

Wprowadzenie Pola siłowe Klasyczna metoda Monte Carlo Kwantowe metody Monte Carlo
Kwantowa dynamika Monte Carlo - PIMC
Akceptacja ruchu
N(n)α(n → o)
acc(o → n) =
N(o)α(o → n) = exp{−βΦ(n)} exp{−βΦ K (o)}
exp{−βΦ(o)} exp{−βΦ K (n)}
= exp{−βΦ V (n)}
exp{−βΦ V (o)} = exp{−βΦ V (x 1 , x2 new , ..., xj new , ..., x P )}
exp{−βΦ V (x 1 , x 2 , ..., x j , ..., x P )}
exp{−(β/P)[V (xnew 2 ) + V (x new
3 ) + ... + V (x new
j )]}
= .
exp{−(β/P)[V (x 2 ) + V (x 3 ) + ... + V (x j )]}
Długość segmentu j wybiera się tak aby prawdopodobieństwo
zaakceptowania ruchu wynosiło 40 − 50%.
Powtórzyć poprzednie kroki P razy i obliczyć estymatory.