Modelowanie molekularne - metody Monte Carlo
Modelowanie molekularne - metody Monte Carlo
Modelowanie molekularne - metody Monte Carlo
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Wprowadzenie Pola siłowe Klasyczna metoda <strong>Monte</strong> <strong>Carlo</strong> Kwantowe <strong>metody</strong> <strong>Monte</strong> <strong>Carlo</strong><br />
Metropolis <strong>Monte</strong> <strong>Carlo</strong><br />
Decyzja czy ruch próbny zostanie zaakceptowany czy odrzucony:<br />
jeżeli wygenerowany został ruch próbny o → n taki, że<br />
U(n) > U(o) to akceptujemy go z prawdopodobieństwem równym<br />
acc(o → n) = exp{−β[U(n) − U(o)]} < 1.<br />
Następnie należy wygenerować liczbę losową o rozkładzie<br />
jednorodnym rand ∈ [0, 1]. Prawdopodobieństwo, że<br />
rand < acc(o → n) wynosi acc(o → n), zatem dany ruch<br />
wykonujemy, jeżeli<br />
rand < acc(o → n),<br />
w przeciwnym wypadku jest on odrzucany.<br />
Ważne jest żeby π(o → n) bylo ergodyczne, czyli pozwalało<br />
osiągnąć każdy punkt przestrzeni konfiguracyjnej z dowolnego<br />
innego punktu w skończonej liczbie kroków MC.