Classical force fields have been broadly used in studies of metal–organic framework crystals. However, processes involving bond breaking or forming are prohibited due to the nonreactive nature of the potentials. Here, we present simulations of three ZIF compounds using a ReaxFF reactive force field.

Enabling Computational Design of ZIFs Using ReaxFF