Hydrogen bonding at the diatomics-in-molecules level: Water clusters

Hybrid diatomics-in-molecules-based quantum mechanical/molecular mechanical approach applied to the modeling of structures and spectra of mixed molecular clusters Arn(HCl)m and Arn(HF)m

A new hybrid QM/DIM approach aimed at describing equilibrium structures and spectroscopic properties of medium size mixed molecular clusters is developed. This methodology is applied to vibrational spectra of hydrogen chloride and hydrogen fluoride clusters with up to four monomer molecules embedded in argon shells Arn(H(Cl/F))m (n = 1-62, m = 1-4). The hydrogen halide complexes (QM part) are treated at the MP2/aug-cc-pVTZ level, while the interaction between HX molecules and Ar atoms (MM part) is described in terms of the semiempirical DIM methodology, based on the proper mixing between neutral and ionic states of the system [Grigorenko et al., J. Chem. Phys. 104, 5510 (1996)]. A detailed analysis of the resulting topology of the QM/DIM potential energy (hyper-)surface in the triatomic subsystem Ar-HX reveals more pronounced nonadditive atomic induction and dispersion contributions to the total interaction energy in the case of the Ar-HCl system. An extension of the original analytical DIM-based potential in the frame of the present model as well as the current limitations of the method are discussed. A modified algorithm for the gradient geometry optimization, along with partly analytical force constant matrix evaluation, is developed to treat large cages of argon atoms around molecular clusters. Calculated frequency redshifts of HX stretching vibrations in the mixed clusters relative to the isolated hydrogen-bonded complexes are in good agreement with experimental findings.

The open-shell interaction of He with the B 3Πu(0+) state of Br2: Anab initiostudy and its comparison with a diatomics-in-molecule perturbation model

The interaction of He with Br2 in electronically excited B 3Piu state is investigated using spin-unrestricted single and double coupled-cluster approach with noniterative perturbative treatment of triple excitations. Internal electrons of the Br atom are described by effective core pseudopotentials. The validity of this approach is analyzed by comparing the lowest 2Sigma+ and 2Pi electronic states of the HeBr molecule with those obtained in all electron calculations [J. Chem. Phys. 115, 10438 (2001)]. In this context, we examine the performance of different basis sets and saturation with bond functions. The comparison of theoretical blue-shifts with the experiment provides confidence about the present ab initio calculations. In addition, He-Br results of ab initio calculations at the same level are used to obtain approximate He-Br2 (3Piu) interactions in the framework of the diatomics-in-molecule first order perturbation theory (IDIM-PT1) [J. Chem. Phys. 104, 9913 (1996)]. Overall, the IDIM-PT1 model results show a good agreement with the ab initio ones, being the main difference the sensitivity to the elongation of the Br-Br bond.

Intermolecular complexes of HXeOH with water: stabilization and destabilization effects

Theoretical and matrix-isolation studies of intermolecular complexes of HXeOH with water molecules are presented. The structures and possible decomposition routes of the HXeOH-(H(2)O)(n)(n = 0, 1, 2, 3) complexes are analyzed theoretically. It is concluded that the decay of these metastable species may proceed through the bent transition states (TSs), leading to the global minima on the respective potential energy surfaces, Xe + (H(2)O)(n+1). The respective barrier heights are 39.6, 26.6, 11.2, and 0.4 kcal/mol for n = 0, 1, 2, and 3. HXeOH in larger water clusters is computationally unstable with respect to the bending coordinate, representing the destabilization effect. Another decomposition channel of HXeOH-(H(2)O)(n), via a linear TS, leads to a direct break of the H-Xe bond of HXeOH. In this case, the attached water molecules stabilize HXeOH by strengthening the H-Xe bond. Due to the stabilization, a large blue shift of the H-Xe stretching mode upon complexation of HXeOH with water molecules is featured in calculations. On the basis of this computational result, the IR absorption bands at 1681 and 1742 cm(-1) observed after UV photolysis and annealing of multimeric H(2)O/Xe matrixes are assigned to the HXeOH-H(2)O and HXeOH-(H(2)O)(2) complexes. These bands are blue-shifted by 103 and 164 cm(-1) from the known monomeric HXeOH absorption.

Water structure theory and some implications for drug design

The development of theories of water structure has been hindered in the past by the difficulty of experimental measurement. Both measurement and computer modelling studies have now reached the stage where theoretical treatments of water structure are converging to a broadly acceptable model. In current understanding, water is a mixture of randomly hydrogen-bonded molecules and larger structures comprised of tetrahedral oxygen centres which, when hydrogen-bonded to each other, lead to five-membered and other rings which can aggregate to form three-dimensional structures. Evidence is taken from studies of the ices, from clathrates and other solid solutions, as well as from liquid solutions, that certain motifs occur very frequently and have relatively high stability, such as the (H2O)20 cavity-forming structure known from studies on clathrates. The implications of recent models of water structure for an understanding of biological events, including the interactions of drugs with receptors, are profound. It is becoming clear that modelling of aqueous solutions of any molecule must consider the explicit interactions with water molecules, which should not be regarded as a continuum: water itself is not a continuum. Solute molecules which possess hydrogen-bonding groups will provoke the formation of further hydrogen-bonding chains of water molecules: if these can form rings, such ringswilltend to persist longerthan chains, givingthesolute a secondary identity of associated water which may play a role in molecular recognition. Solutes that do not have hydrogen-bonding capability, or regions of solutes which are non-polar, may also produce partial cage-like water structures that are characteristic of the solute. The classification of many solutes as structure makers or structure breakers has relevance to the interactions between ligands and large biomolecules such as proteins. While it is generally accepted that sulfate and urea, respectively structure maker and breaker, may alter protein conformation through effects on water, it has not been recognised that bioactive ligands, which also change the conformation of proteins, may do so by a related, but more selective, mechanism. Very early studies of cell contents suggested that the associated water might be different from bulk water, a concept that lost support in the mid-20th century. Current theories of water structure may invite a reappraisal of this position, given the observation that structuring may extend for many molecular diameters from an ordered surface.