Selection and Validation of Charge and Lennard-Jones Parameters for QM/MM Simulations of Hydrocarbon Interactions with Zeolites

Quantum mechanics/molecular mechanics (QM/MM) models are an appealing method for performing zeolite simulations. In QM/MM, a small cluster chosen to encompass the active center is described by QM, while the rest of the zeolite is described by MM. In the present study, we demonstrate that the charges and Lennard-Jones parameters on Si and O must be chosen properly for QM/MM calculations of adsorption energies and activation energies to agree closely with full QM calculations. The selection of parameters for Si and O is based on using the omega B97X-D functional for DFT calculations of the QM region, which is effective in capturing the effects of van der Waals interactions. A comparison of the heats of adsorption for a variety of adsorbates and activation energies for the cracking of propane and butane reveals that energies derived from QM/MM calculation carried out with appropriately selected MM parameters agree to within an rim error of similar to 1.5 kcal/mol with QM calculations. To avoid reparametrization for new substrates, Lennard-Jones zeolite parameters are chosen to be compatible with existing CHARMM parameters. Transferability of these parameters is demonstrated by tests utilizing the B3LYP density functional and simulations of MFI and FAU zeolites. Moreover, the computational time for QM/MM calculations is considerably lower than that for QM calculations, and the ratio of computational times decreases rapidly with increasing size of the cluster used to represent the zeolite.