Modelowanie molekularne - metody Monte Carlo

Wprowadzenie Pola siłowe Klasyczna metoda Monte Carlo Kwantowe metody Monte Carlo
Kwantowa dynamika Monte Carlo - PIMC
Jednoelementowe ruchy typu cząstki swobodnej:
Wybierz losowo jeden element łańcucha.
Wygeneruj nowe położenie zgodnie z dystrybucja jednorodną:
√
N(x, σ) = N((x k−1 + x k+1 )/2, ħ 2 β/(2mP))
α(o → n) = exp{−βΦ K (n)} =
exp{−mP/2ħ 2 β[... + (x k−1 − x new
exp{−1/2[ xnew k −(x k−1 +x k+1 )/2
σ
] 2 + ...}
Akceptacja ruchu
acc(o → n) =
k
) 2 + (x new
k
− x k+1 ) 2 + ...]} ∝
N(n)α(n → o)
N(o)α(o → n) = exp{−βΦ(n)} exp{−βΦ K (o)}
exp{−βΦ(o)} exp{−βΦ K (n)}
= exp{−βΦ V (n)}
exp{−βΦ V (o)}
=
exp[−(β/P)V (xnew K )]
exp[−(β/P)V (x K )] .