20.11.2014 Views

Modelowanie molekularne - metody Monte Carlo

Modelowanie molekularne - metody Monte Carlo

Modelowanie molekularne - metody Monte Carlo

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Wprowadzenie Pola siłowe Klasyczna metoda <strong>Monte</strong> <strong>Carlo</strong> Kwantowe <strong>metody</strong> <strong>Monte</strong> <strong>Carlo</strong><br />

Algorytm DMC<br />

1 Obliczenie liczby kopii, które kontynuują symulację<br />

M new = INT (η + exp{−τ[E L (R) + E L (R ′ ) − 2E T ]/2}),<br />

gdzie η - liczba losowa o rozkładzie jednorodnym na przedziale<br />

[0, 1].<br />

2 Zachować wielkości, którymi się interesujemy (np. średnie E L<br />

ze zbioru “walkerów“).<br />

3 Po początkowym rozgrzaniu symulacji kroki 2 − 7 powtarzać<br />

do uzyskania wystarczająco małego błędu. Wartość E T jest co<br />

pewien czas dostosowywana tak, aby utrzymać średnią ilość<br />

”walkerów“ według wzoru<br />

E T ← E T − C E ln(M act /M ave ),<br />

gdzie C E jest dodatnią stałą kontrolującą szybkość osiągania<br />

zakładanej liczby M ave .

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!