Hybrid approach for generating realistic amorphous carbon structure using metropolis and reverse Monte Carlo

Abstract

An improved method for the modelling of carbon structures based on a hybrid reverse Monte Carlo (HRMC) method is presented. This algorithm incorporates an accurate environment dependent interaction potential (EDIP) in conjunction with the commonly used constraints derived from experimental data. In this work, we compare this new method with other modelling results for a small system of 2.9 g/cc amorphous carbon. We find that the new approach greatly improves the structural description, alleviating the common problem in standard reverse Monte Carlo method (RMC) of generating structures with a high proportion of unphysical small rings. The advantage of our method is that larger systems can now be modelled, allowing the incorporation of mesoscopic scale features.