Effect of charge distribution on RDX adsorption in IRMOF-10

Quantum mechanical (QM) calculations, classical grand canonical Monte Carlo (GCMC) simulations, and classical molecular dynamics (MD) simulations are performed to test the effect of charge distribution on hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) adsorption and diffusion in IRMOF-10. Several different methods for mapping QM electron distributions onto atomic point charges are explored, including the electrostatic potential (ESP) method, Mulliken population analysis, Lowdin population analysis, and natural bond orbital analysis. Classical GCMC and MD simulations of RDX in IRMOF-10 are performed using 15 combinations of charge sources of RDX and IRMOF-10. As the charge distributions vary, interaction potential energies, the adsorption loading, and the self-diffusivities are significantly different. None of the 15 combinations are able to quantitatively capture the dependence of the energy of adsorption on local configuration of RDX as observed in the QM calculations. We observe changes in the charge distributions of RDX and IRMOF-10 with the introduction of an RDX molecule into the cage. We also observe a large dispersion contribution to the interaction energy from QM calculations that is not reproduced in the classical simulations, indicating that the source of discrepancy may not lie exclusively with the assignment of charges.