Abstract
Various metal-organic frameworks (MOFs) containing trivalent cations (such as Fe(3+), Al(3+), and Cr(3+)) have been reported and have shown great potential in applications. However, the high structural diversity and strong electronic interactions between metal centers and their ligands make the molecular dynamics simulations of MOFs challenging. In this work, we developed new dummy atom models for Fe(3+), Al(3+), and Cr(3+) cations, which can be used in classical molecular dynamics simulations of MOFs. In our models, the correct solvation free energies and metal-ligand distances can be simultaneously reproduced. Furthermore, the usefulness and transferability of our models were validated using the commonly studied MIL-100(M) (M = Fe(3+), Al(3+), Cr(3+)) and MIL-88B(Fe(3+)) systems. Our developed models offer a valuable tool for simulating complex systems containing Fe(3+), Al(3+), and Cr(3+) cations with octahedral coordination structures.