Symulacja MC w układzie Lennarda - Jonesa

Algorytm programu:

1.) Podajemy parametry początkowe:

  • Tytuł projektu.
  • Liczbę cykli MC
  • Co ile kroków wyniki obliczeń będą wyświetlane na ekranie.
  • Co ile kroków konfiguracja będzie zapisywana do pliku.
  • Stała potrzebna do stwierdzenia co ile kroków ma być modyfikowane maksymalne przemieszczenie cząsteczki (domyślnie = 0,15sigma).
  • Nazwa pliku z konfiguracją początkową (wcześniej przygotowaną).
  • Nazwa pliku pod który ma być zapisana konfiguracja końcowa
  • Gęstość bezwymiarowa:
klj1
  • Temperatura bezwymiarowa:
klj2
  • Zasięg potencjału:
klj3

2.) Program liczy potrzebne rzeczy (np. odwrotność temperatury, zerowanie liczników).

3.) Uruchamiana jest procedura SUMUP - sprawdzana jest konfiguracja (czy kulki na siebie nie nachodzą - jeśli tak jest to zatrzymanie programu) oraz liczona jest całkowita energia (suma każdej pary atom - atom, chyba że odległość między atomami jest większa niż Rcut).

4.) Ze wzoru:
klj4

    obliczane jest ciśnienie.

5.) Rozpoczęcie łańcucha Markowa.
NStep razy liczone jest:

    Energia I-ty atom - inne atomy (procedura ENERGY),
    Losowe przesunięcie I-tego atomu,
    Energia I-ty atom - inne atomy (po przesunięciu)
    Obliczenie różnicy energii deltaV=Vpo przesunieciu - Vprzed.
    Jeżeli deltaV<0 przesunięcie jest akceptowane (deltaV<0 ( odpowiada ruchom zmniejszającym energię)
    Jeżeli deltaV>0 dokonujemy sprawdzenia: exp(-deltaV/T*)>liczba losowa - jeśli tak jest to przesunięcie także akceptujemy;
    Jeśli : exp(-deltaV/T*)<liczba losowa to przesunięcia nie akceptujemy
boltzman

    Liczone jest ciśnienie i inne rzeczy (liczniki)
    Tak po wszystkich cząsteczkach.
    Zapisywanie każdego kroku do pliku    .

Liczona jest całkowita energia.
Zapisanie konfiguracji końcowej.

*(demonstracja programu)*

W metodzie MC ważnym elementem są generatory liczb losowych.

<< Wstecz

Dalej >>