Simulations of the organic-molecular ordering on the solid surfaces using pair potentials

Abstract

Organic molecules on a flat, unreactive, periodical surface self-assemble to ordered superstructures. The superstructures can be used to selectively tune the properties of the surface. This problem can be studied using quantum mechanical models, such as Density functional theory (DFT). The use of those models for problems with hundreds of atoms and thousand of iterations is, however, due to their complexity, computationally impractical. In this thesis we extend and test an existing algorithm developed by the advisor of this thesis. The algorithm sacrifices accuracy for a quicker simulation time by using the method of paired potentials. This method simplifies the quantum mechanical description of an atom by considering a point source of change and replaces the complex interaction between atoms with pair parametrized forces. This method is incorporated into a Monte Carlo simulation which, from a given superstructure obtained from exper- imental data, predicts the position and deformation of the molecules that minimize the total potential energy. We test the algorithm on two real-world problem sets. The algo- rithm gives results comparable to that of DFT when tested on already studied problems, yet converges to the solution much faster. It can be used as a standalone program or as an initialization to the DFT...

Organické molekuly rozptýlené na plochý, nereaktivní, periodický povrch se spontálně uspořádají do superstruktur, které mohou být použity k ladění vlastností povrchu. Tento problém může být studovám za použití velmi přesných kvantově mechanických výpočet- ních modelů. Použití těchto modelů na problémy obsahující stovky atomů a tisíce iter- ací je, nicméně, kvůli jejich složitosti, výpočetně nepraktické. V této práci rozšiřujeme a testujeme existující algoritmus používaný vedoucím této práce. Algoritmus používá metodu párových potenciálů pro zrychlení výpočtů, avšak na úkor přesnosti. Tato metoda zjednodušuje kvantově mechanický popis atomu a uvažuje atom jako bodový náboj a mění komplexní interakci mezi atomy na párově parametrizovatelné funkce. Tato metoda je in- tegrovaná do Monte Carlo simulace, která pro zadanou superstrukturu, získanou z exper- imentálních dat, předpoví konfiguraci molekul v superstruktuře minimalizující celkovou potenciální energii. Algoritmus dává výsledky srovnatelné s přesnějšími, kvantově me- chanickými modely a konverguje k řešení mnohem rychleji. Lze ho použít samostatně nebo jako inicializaci výpočtů založených na kvantové mechanice. 1