The effect of atomic point charges on adsorption isotherms of CO2 and water in metal organic frameworks

The interactions between metal–organic frameworks (MOFs) and adsorbates have been increasingly predicted and studied by computer simulations, particularly by Grand-Canonical Monte Carlo (GCMC), as this method enables comparing the results with experimental data and also provides a degree of molecular level detail that is difficult to obtain in experiments. The assignment of atomic point charges to each atom of the framework is essential for modelling Coulombic interactions between the MOF and the adsorbate. Such interactions are important in adsorption of polar gases like water or carbon dioxide, both of which are central in carbon capture processes. The aim of this work is to systematically investigate the effect of varying atomic point charges on adsorption isotherm predictions, identify the underlying trends, and based on this knowledge to improve existing models in order to increase the accuracy of gas adsorption prediction in MOFs. Adsorption isotherms for CO2 and water in several MOFs were generated with GCMC, using the same computational parameters for each material except framework point charge sets that were obtained through a wide range of computational approaches. We carried out this work for 6 widely studied MOFs; IRMOF-1, MIL-47, UiO-66, CuBTC, Co-MOF-74 and SIFSIX-2-Cu-I. We included both MOFs with and without open metal sites (OMS), specifically to investigate whether this property affects the predicted adsorption behaviour. Our results show that point charges obtained from quantum mechanical calculations on fully periodic structures are generally more consistent and reliable than those obtained from either cluster-based QM calculations or semi-empirical approaches. Furthermore, adsorption in MOFs that contain OMS is much more sensitive to the point charge values, with particularly large variability being observed for water adsorption in such MOFs. This suggests that particular care must be taken when simulating adsorption of polar molecules in MOFs with open metal sites to ensure that accurate results are obtained.