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.

Compressive failure of finite size unidirectional composite laminates with a region of fibre waviness

The compressive behaviour of finite unidirectional composites with a region of misaligned reinforcement is investigated via finite element analyses. Models with and without fibre bending stiffness are compared, confirming that compressive strength is accurately predicted without modelling fibre bending stiffness for real composite components which typically have waviness defects of several millimetres wavelength. Various defect parameters are investigated. Results confirm the well-known sensitivity of compressive strength to misalignment angle, and also show that compressive strength falls rapidly with the proportion of laminate width covered by the wavy region. A simple empirical equation is proposed to model the effect of a single patch of waviness in finite specimens. Other parameters such as length and position of the wavy region are found to have a smaller effect on compressive strength. The modelling approach is finally adapted to model distributed waviness and thus determine the compressive strength of composites with realistic waviness defects. © 2011 Elsevier Ltd. All rights reserved.

QM/MM simulation of liquid water with an adaptive quantum region.

The simulation of complex chemical systems often requires a multi-level description, in which a region of special interest is treated using a computationally expensive quantum mechanical (QM) model while its environment is described by a faster, simpler molecular mechanical (MM) model. Furthermore, studying dynamic effects in solvated systems or bio-molecules requires a variable definition of the two regions, so that atoms or molecules can be dynamically re-assigned between the QM and MM descriptions during the course of the simulation. Such reassignments pose a problem for traditional QM/MM schemes by exacerbating the errors that stem from switching the model at the boundary. Here we show that stable, long adaptive simulations can be carried out using density functional theory with the BLYP exchange-correlation functional for the QM model and a flexible TIP3P force field for the MM model without requiring adjustments of either. Using a primary benchmark system of pure water, we investigate the convergence of the liquid structure with the size of the QM region, and demonstrate that by using a sufficiently large QM region (with radius 6 Å) it is possible to obtain radial and angular distributions that, in the QM region, match the results of fully quantum mechanical calculations with periodic boundary conditions, and, after a smooth transition, also agree with fully MM calculations in the MM region. The key ingredient is the accurate evaluation of forces in the QM subsystem which we achieve by including an extended buffer region in the QM calculations. We also show that our buffered-force QM/MM scheme is transferable by simulating the solvated Cl(-) ion.

Phycobilisome-Deficient Strains of Synechocystis sp. PCC 6803 Have Reduced Size and Require Carbon-Limiting Conditions to Exhibit Enhanced Productivity.

Reducing excessive light harvesting in photosynthetic organisms may increase biomass yields by limiting photoinhibition and increasing light penetration in dense cultures. The cyanobacterium Synechocystis sp. PCC 6803 harvests light via the phycobilisome, which consists of an allophycocyanin core and six radiating rods, each with three phycocyanin (PC) discs. Via targeted gene disruption and alterations to the promoter region, three mutants with two (pcpcT→C) and one (ΔCpcC1C2:pcpcT→C) PC discs per rod or lacking PC (olive) were generated. Photoinhibition and chlorophyll levels decreased upon phycobilisome reduction, although greater penetration of white light was observed only in the PC-deficient mutant. In all strains cultured at high cell densities, most light was absorbed by the first 2 cm of the culture. Photosynthesis and respiration rates were also reduced in the ΔCpcC1C2:pcpcT→C and olive mutants. Cell size was smaller in the pcpcT→C and olive strains. Growth and biomass accumulation were similar between the wild-type and pcpcT→C under a variety of conditions. Growth and biomass accumulation of the olive mutant were poorer in carbon-saturated cultures but improved in carbon-limited cultures at higher light intensities, as they did in the ΔCpcC1C2:pcpcT→C mutant. This study shows that one PC disc per rod is sufficient for maximal light harvesting and biomass accumulation, except under conditions of high light and carbon limitation, and two or more are sufficient for maximal oxygen evolution. To our knowledge, this study is the first to measure light penetration in bulk cultures of cyanobacteria and offers important insights into photobioreactor design.