Evaluating boundary dependent errors in QM/MM simulations.

Hybrid quantum mechanics/molecular mechanics (QM/MM) simulations provide a powerful tool for studying chemical reactions, especially in complex biochemical systems. In most works to date, the quantum region is kept fixed throughout the simulation and is defined in an ad hoc way based on chemical intuition and available computational resources. The simulation errors associated with a given choice of the quantum region are, however, rarely assessed in a systematic manner. Here we study the dependence of two relevant quantities on the QM region size: the force error at the center of the QM region and the free energy of a proton transfer reaction. Taking lysozyme as our model system, we find that in an apolar region the average force error rapidly decreases with increasing QM region size. In contrast, the average force error at the polar active site is considerably higher, exhibits large oscillations and decreases more slowly, and may not fall below acceptable limits even for a quantum region radius of 9.0 A. Although computation of free energies could only be afforded until 6.0 A, results were found to change considerably within these limits. These errors demonstrate that the results of QM/MM calculations are heavily affected by the definition of the QM region (not only its size), and a convergence test is proposed to be a part of setting up QM/MM simulations.

Self-Similarity and helical symmetry in vortex generator flow simulations

According to experimental observations, the vortices generated by vortex generators have previously been observed to be self-similar for both the axial (u(z)) and azimuthal (u(circle minus)) velocity profiles. Further, the measured vortices have been observed to obey the criteria for helical symmetry. This is a powerful result, since it reduces the highly complex flow to merely four parameters. In the present work, corresponding computer simulations using Reynolds-Averaged Navier-Stokes equations have been carried out and compared to the experimental observations. The main objective of this study is to investigate how well the simulations can reproduce the physics of the flow and if the same analytical model can be applied. Using this model, parametric studies can be significantly reduced and, further, reliable simulations can substantially reduce the costs of the parametric studies themselves.