PairPotMCinator: a tool for fast simulations of the organic-molecular ordering on the solid surfaces using pair...

Organic molecules on nonreactive surfaces self-assemble to ordered superstructures. This problem is typically studied using quantum mechanical models, such as Density functional theory (DFT). However, the use of those models for systems with hundreds of atoms is computationally impractical if at least approximate structures are not known. Here, we demonstrate an algorithm employing a Monte Carlo simulation and a method of pair potentials - thus simplifying the quantum mechanical description by using pair-parameterized forces. We incorporate this model into a Monte Carlo simulation which, from a given approximate molecular superstructure periodicity obtained from experimental data, predicts the position and deformation of the ordered molecules to minimize the potential energy of the system. We test the algorithm on published results of self-ordering of selected pigments in different environments: 1) on highly oriented pyrolytic graphite in solution, and 2) on the Si(111)-In[Formula: see text] surface prepared under ultra-high vacuum. The generated structures can be used as close-to-relaxed candidates for inputs to more accurate DFT calculations.