Sposoby rozpoczęcia ab initio MD od klasycznej MD

9

Prowadzę symulacje dynamiki molekularnej wody do celów testowych. Pudełko jest dość małe, jeśli zapytasz faceta z klasyczną MD, i stosunkowo duże, jeśli zapytasz faceta z DFT: Mam 58 cząsteczek wody w okresowych warunkach brzegowych.

Aby zaoszczędzić czas procesora, optymalizuję komórkę za pomocą klasycznego pola siłowego przed uruchomieniem ab initio MD. System równoważę klasycznie w 300K przez 1 ns, a następnie robię ostatnią migawkę i używam jej jako danych wejściowych dla ab initio MD. Mój ab initio MD to zwykły oparty na DFT Born-Oppenheimer MD z zestawem podstawy fali płaskiej i potencjałem PAW (pseudo) (VASP to kod). Zarówno w symulacji klasycznej, jak i ab initio utrzymuję stałą temperaturę na poziomie 300 K za pomocą termostatu ze zmianą prędkości.

Badam dwa różne sposoby przejścia od klasycznego do ab initio:

  1. Pobierz początkowe prędkości i pozycje z klasycznej trajektorii i zaimportuj je jako początkową konfigurację do symulacji ab initio
  2. Zamroź system do zera, zachowując klasyczne pozycje, zaimportuj go do kodu DFT, a następnie szybko (obecnie robię to w 0,5 ps) rozgrzej do 300K

Miałem nadzieję, że obie strategie doprowadzą do tej samej średniej energii po krótkim (powiedzmy 10 ps) okresie równoważenia, szczególnie biorąc pod uwagę, że konfiguracja początkowa jest dokładnie taka sama (te same pozycje początkowe), z wyjątkiem wspomnianej sztuczki temperaturowej (prędkości początkowe różnią się) . Nie o to chodzi. Poniższy rysunek pokazuje, że symulacja, w której układ jest zamrażany, a następnie szybko rozgrzewany, znajduje obszar energii o 1 eV niższy pod względem energii niż drugi, gdzie prędkości są importowane z klasycznej MD.

wprowadź opis zdjęcia tutaj

Moje pytania to:

  1. czy tego należy się spodziewać;
  2. czy znane są skuteczne strategie optymalizacji przejścia z klasycznej na ab initio MD;
  3. i czy możesz wskazać mi odpowiednią literaturę na ten temat?

Edytować:

Przeprowadziłem jeszcze kilka testów i - przy ograniczonych obecnie danych - wydaje się, że może to być problem specyficzny dla systemu. Test z metanolem zamiast wody w pudełku o tym samym rozmiarze wykazał, że dwa różne schematy prędkości początkowej szybko zbiegają się do tej samej średniej energii. Jednak klasyczna konfiguracja była bardzo zbliżona do konfiguracji kwantowej w przypadku metanolu, to znaczy energia przy t = 0 była bardzo zbliżona do średniej energii po konwergencji. Woda jest niezwykle trudnym systemem, więc być może ten problem jest mniej więcej specyficzny dla wody. Jeśli nie zostaną dodane żadne odpowiedzi, spróbuję opublikować jedną na podstawie moich wyników, gdy skończę wszystkie testy.

Miguel
źródło

Odpowiedzi:

1

po krótkim (powiedzmy 10 ps) okresie równoważenia

  1. Mówisz sobie, że okres ponownej równowagi jest „krótki”. Czy próbowałeś dłużej czekać, aby zobaczyć, czy dwa ponownie zainicjowane systemy się zbiegają , a jeśli tak, to w jakim tempie?
  2. Skalowanie prędkości jest naiwnie naiwne termostatem. Może mógłbyś zastąpić to czymś bardziej realistycznym? (Berendsen, Nose-Hoover itp.)
  3. Jeśli w ogóle martwisz się próbkowaniem „dużego” systemu, możesz sprawić, by zasada ergodyczna i obliczenia równoległe działały dla Ciebie: Próbuj wiele ponownie zainicjowanych realizacji z końca swojej klasycznej trajektorii (robienie migawek klasycznie zrównoważonego systemu) . Następnie możesz uruchomić wiele wystąpień kodu kwantowego, każdy używając innej konfiguracji przestrzeni fazowej ponownej inicjalizacji i uśrednić wyniki. Ponieważ poszczególne realizacje ewoluują niezależnie od siebie, symulacje są żenująco równoległe .

Niestety nie mogę zaoferować żadnych wiarygodnych źródeł ...

Inon
źródło