We present an accurate adaptive multiscale molecular dynamics method that will enable the detailed study of large molecular systems that mimic experiment. The method treats the reactive regions at the quantum mechanical level and the inactive environment regions at lower levels of accuracy, while at the same time molecules are allowed to flow across the border between active and environment regions. Among many other things, this scheme affords accurate investigation of chemical reactions in solution. A scheme like this ideally fulfills the key criteria applicable to all molecular dynamics simulations: energy conservation and computational efficiency. Approaches that fulfill both criteria can, however, result in complicated potential energy surfaces, creating rapid energy changes when the border between regions is crossed. With the difference-based adaptive solvation potential, a simple approach is introduced that meets the above requirements and reduces fast fluctuations in the potential to a minimum. In cases where none of the current adaptive QM/MM potentials are able to properly describe the system under investigation, we use a continuous force scheme instead, which, while no longer energy conserving, still retains a related conserved quantity along the trajectory. We show that this scheme does not introduce a significant temperature drift on time scales feasible for QM/MM simulations.

Recently, the outstanding properties of Cu-SSZ-13 (a zeolite in the chabazite structure) for the selective catalytic reduction of nitrous oxides were discovered. However, the true nature of the active site is still not answered satisfactorily. In this work, we identify the active site for the given reaction from first-principles simulations of the total energy of Cu(II) ions in various positions in combination with previously published catalytic activity as a function of the copper exchange level. This attribution is confirmed by the simulation of vibrational properties of CO adsorbed to the reduced Cu(I) species. The relation between energetic considerations, vibrational calculations, and experiment allows a clear statement about the distribution of active sites in the catalyst. We furthermore discuss the structural properties of the active site leading to the high stability under reaction conditions over a large temperature range. The insights from this work allow a more targeted catalyst design and represent a step toward an industrial application of copper-exchanged zeolites in cleaning car exhaust gases.