Simulation of water cluster assembly on a graphite surface

The assembly of small water clusters (H2O)n, n = 1-6, on a graphite surface is studied using a density functional tight-binding method complemented with an empirical van der Waals force correction, with confirmation using second-order M?ller-Plesset perturbation theory. It is shown that the optimized geometry of the water hexamer may change its original structure to an isoenergy one when interacting with a graphite surface in some specific orientation, while the smaller water cluster will maintain its cyclic or linear configurations (for the water dimer). The binding energy of water clusters interacting with graphite is dependent on the number of water molecules that form hydrogen bonds, but is independent of the water cluster size. These physically adsorbed water clusters show little change in their IR peak position and leave an almost perfect graphite surface.     

Validation of the density-functional based tight-binding approximation method for the calculation of reaction energies and other data

We investigated the performance of the approximative density functional method DFTB versus BLYP and G2 with respect to zero-point corrected reaction energies, vibrational frequencies, and geometry parameters for a set of 28 reactions and 22 representative molecules containing C, H, N, and O (DFTB--density-functional based tight-binding approximation). The DFTB reaction energies show a mean absolute deviation versus the G2 reference of 4.3 kcalmol only. The corresponding value for the vibrational frequencies amounts to 75 cm(-1) versus BLYP/cc-pVTZ. With very few exceptions bond lengths and angles are in excellent agreement with the results of higher-level methods.     

Rotamer decomposition and protein dynamics: Efficiently analyzing dihedral populations from molecular dynamics

Molecular mechanics methods have matured into powerful methods to understand the dynamics and flexibility of macromolecules and especially proteins. As multinanosecond to microsecond length molecular dynamics (MD) simulations become commonplace, advanced analysis tools are required to generate scientifically useful information from large amounts of data. Some of the key degrees of freedom to understand protein flexibility and dynamics are the amino acid residue side chain dihedral angles. In this work, we present an easily automated way to summarize and understand the relevant dihedral populations. A tremendous reduction in complexity is achieved by describing dihedral timeseries in terms of histograms decomposed into Gaussians. Using the familiar and widely studied protein lysozyme, it is demonstrated that our approach captures essential properties of protein structure and dynamics. A simple classification scheme is proposed that indicates the rotational state population for each dihedral angle of interest and allows a decision if a given side chain or peptide backbone fragment remains rigid during the course of an MD simulation, adopts a converged distribution between conformational substates or has not reached convergence yet. © 2012 Wiley Periodicals, Inc.     

Host-Guest Chemistry of Truncated Tetrahedral Imine Cages with Ammonium Ions

Three shape-persistent [4+4] imine cages with truncated tetrahedral geometry with different window sizes were studied as hosts for the encapsulation of tetra-n-alkylammonium salts of various bulkiness. In various solvents the cages behave differently. For instance, in dichloromethane the cage with smallest window size takes up NEt₄⁺ but not NMe₄ ^{+ ,} which is in contrast to the two cages with larger windows hosting both ions. To find out the reason for this, kinetic experiments were carried out to determine the velocity of uptake but also to deduce the activation barriers for these processes. To support the experimental results, calculations for the guest uptakes have been performed by molecular mechanics’ simulations. Finally, the complexation of pharmaceutical interested compounds, such as acetylcholine, muscarine or denatonium have been determined by NMR experiments.     

A global investigation of excited state surfaces within time-dependent density-functional response theory

This work investigates the capability of time-dependent density functional response theory to describe excited state potential energy surfaces of conjugated organic molecules. Applications to linear polyenes, aromatic systems, and the protonated Schiff base of retinal demonstrate the scope of currently used exchange-correlation functionals as local, adiabatic approximations to time-dependent Kohn-Sham theory. The results are compared to experimental and ab initio data of various kinds to attain a critical analysis of common problems concerning charge transfer and long range (nondynamic) correlation effects. This analysis goes beyond a local investigation of electronic properties and incorporates a global view of the excited state potential energy surfaces.     

Effect of nitrogen adsorption on the mid-infrared spectrum of water clusters

Experimental Fourier-transform infrared spectra and DFT calculated infrared spectra are compared to investigate the effect of adsorbed nitrogen on the OH-stretch band complex of water clusters. Using a collisional cooling experiment, pure as well as partially and completely N(2)-covered water clusters consisting of 20-200 water molecules have been generated in thermal equilibrium in the aerosol phase within the temperature range of 5-80 K. Computational IR-spectra simulations have been performed for discrete pure and N(2)-covered water clusters including 10, 15, 20, and 30 water molecules. The adsorbed N(2) molecules especially affect the three-coordinated water molecules at the cluster surface which could be observed as a blue shift of the companion O-H band at 2900 cm(-1) and a red shift of the dangling O-H band at 3700 cm(-1) by about 20 cm(-1) in both cases. The most striking effect of the N(2) adsorbate is an intensity increase of the dangling O-H band by a factor of 3-5. Furthermore, the onset temperature of nitrogen adsorption at the water cluster surface was experimentally found to be roughly 30 K for cluster sizes of about 100 water molecules. Experimental and computational results are in good agreement. The presented results are based on and support the work of V. Buch, J. P. Devlin, and co-workers (e.g., J. Phys. Chem. B, 1997; J. Phys. Chem. A, 2003; Int. Rev. Phys. Chem., 2004).     

Di- and Tetracyano-Substituted Pyrene-Fused Pyrazaacenes: Aggregation in the Solid State

Five di- and tetracyano-substituted pyrene-fused pyrazaacenes were synthesized and studied as potential electron acceptors in the solid state. Single crystals of all compounds were grown and the crystal packing studied by DFT calculations (transfer integrals and reorganization energies) to get insight into possible use for semiconducting charge transport.     

Effect of polarization on the opsin shift in rhodopsins. 1. A combined QM/QM/MM model for bacteriorhodopsin and pharaonis sensory rhodopsin II

The optical and IR-spectroscopic properties of the protonated Schiff base of retinal are highly sensitive to the electrostatic environment. This feature makes retinal a useful probe to study structural differences and changes in rhodopsins. It also raises an interest to theoretically predict the spectroscopic response to mutation and structural evolution. Computational models appropriate for this purpose usually combine sophisticated quantum mechanical (QM) methods with molecular mechanics (MM) force fields. In an effort to test and improve the accuracy of these QM/MM models, we consider in this article the effects of polarization and inter-residual charge transfer within the binding pocket of bacteriorhodopsin (bR) and pharaonis sensory rhodopsin II (psRII, also called pharaonis phoborhodopsin, ppR) on the excitation energy using an ab initio QM/QM/MM approach. The results will serve as reference for assessing empirical polarization models in a consecutive article.