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.

Tests of an adaptive QM/MM calculation on free energy profiles of chemical reactions in solution.

We present reaction free energy calculations using the adaptive buffered force mixing quantum mechanics/molecular mechanics (bf-QM/MM) method. The bf-QM/MM method combines nonadaptive electrostatic embedding QM/MM calculations with extended and reduced QM regions to calculate accurate forces on all atoms, which can be used in free energy calculation methods that require only the forces and not the energy. We calculate the free energy profiles of two reactions in aqueous solution: the nucleophilic substitution reaction of methyl chloride with a chloride anion and the deprotonation reaction of the tyrosine side chain. We validate the bf-QM/MM method against a full QM simulation, and show that it correctly reproduces both geometrical properties and free energy profiles of the QM model, while the electrostatic embedding QM/MM method using a static QM region comprising only the solute is unable to do so. The bf-QM/MM method is not explicitly dependent on the details of the QM and MM methods, so long as it is possible to compute QM forces in a small region and MM forces in the rest of the system, as in a conventional QM/MM calculation. It is simple, with only a few parameters needed to control the QM calculation sizes, and allows (but does not require) a varying and adapting QM region which is necessary for simulating solutions.

Collapse of a composite beam made from ultra high molecular-weight polyethylene fibres

Hot-pressed laminates with a [0/90]48 lay-up, consisting of 83% by volume of ultra high molecular-weight polyethylene (UHMWPE) fibres, and 17% by volume of polyurethane (PU) matrix, were cut into cantilever beams and subjected to transverse end-loading. The collapse mechanisms were observed both visually and by X-ray scans. Short beams deform elastically and collapse plastically in longitudinal shear, with a shear strength comparable to that observed in double notch, interlaminar shear tests. In contrast, long cantilever beams deform in bending and collapse via a plastic hinge at the built-in end of the beam. The plastic hinge is formed by two wedge-shaped microbuckle zones that grow in size and in intensity with increasing hinge rotation. This new mode of microbuckling under macroscopic bending involves both elastic bending and shearing of the plies, and plastic shear of the interface between each ply. The double-wedge pattern contrasts with the more usual parallel-sided plastic microbuckle that occurs in uniaxial compression. Finite element simulations and analytical models give additional insight into the dominant material and geometric parameters that dictate the collapse response of the UHMWPE composite beam in bending. Detailed comparisons between the observed and predicted collapse responses are used in order to construct a constitutive model for laminated UHMWPE composites. © 2013 Elsevier Ltd.