Theoretical Study on the Hydrogen Bonding Interactions in Complexes of 5‐Hydroxytryptamine with Water

The energies, geometries and harmonic vibrational frequencies of 1:1 5‐hydroxytryptamine‐water (5‐HT‐H₂O) complexes are studied at the MP2/6‐311++G(d,p) level. Natural bond orbital (NBO), quantum theory of atoms in molecules (QTAIM) analyses and the localized molecular orbital energy decomposition analysis (LMO‐EDA) were performed to explore the nature of the hydrogen‐bonding interactions in these complexes. Various types of hydrogen bonds (H‐bonds) are formed in these 5‐HT‐H₂O complexes. The intermolecular C4H5^5‐HT···O^w H‐bond in HTW3 is strengthened due to the cooperativity, whereas no such cooperativity is found in the other 5‐HT‐H₂O complexes. H‐bond in which nitrogen atom of amino in 5‐HT acted as proton donors was stronger than other H‐bonds. Our researches show that the hydrogen bonding interaction plays a vital role on the relative stabilities of 5‐HT‐H₂O complexes.     

Nature of the water/aromatic parallel alignment interactions

The water/aromatic parallel alignment interactions are interactions where the water molecule or one of its O? H bonds is parallel to the aromatic ring plane. The calculated energies of the interactions are significant, up to ΔE_{CCSD(T)(limit)} = ?2.45 kcal mol^?1 at large horizontal displacement, out of benzene ring and CH bond region. These interactions are stronger than CH···O water/benzene interactions, but weaker than OH···π interactions. To investigate the nature of water/aromatic parallel alignment interactions, energy decomposition methods, symmetry‐adapted perturbation theory, and extended transition state‐natural orbitals for chemical valence (NOCV), were used. The calculations have shown that, for the complexes at large horizontal displacements, major contribution to interaction energy comes from electrostatic interactions between monomers, and for the complexes at small horizontal displacements, dispersion interactions are dominant binding force. The NOCV‐based analysis has shown that in structures with strong interaction energies charge transfer of the type π → σ*(O? H) between the monomers also exists. © 2014 Wiley Periodicals, Inc.     

Supramolecular networks constructed by cationic copper(II) complexes incorporating hybrid water–anionic polymeric assemblies

Two complexes, [Cu₂(TFSA)(2,2′-bpy)₄]?·?TFSA?·?8H₂O (1) and {[Cu(4,4′-bpy)(H₂O)₂]?·?TFSA?·?6H₂O} _n (2) (H₂TFSA?=?tetrafluorosuccinic acid, 2,2′-bpy?=?2,2′-bipyridine, and 4,4′-bpy?=?4,4′-bipyridine), have been synthesized and structurally characterized by X-ray structural analyses. Complex 1 is a binuclear molecule bridged by TFSA ligands; 2 is a 1-D chain bridged by 4,4′-bpy ligands. The asymmetric units of the two complexes are composed of cationic complexes [Cu₂(TFSA)(2,2′-bpy)₄]²⁺ (1) and [Cu(4,4′-bpy)(H₂O)₂]²⁺ (2), free TFSA anion, and independent crystallization water molecules. A unique 2-D hybrid water–TFSA anionic layer by linkage of {[(H₂O)₈(TFSA)]^2?} _n fragments consisting of 1-D T6(0)A2 water tape and TFSA anionic units by hydrogen bonds in 1 was observed. Unique 2-D hybrid water–TFSA anionic layer generated by the linkage of {[(H₂O)₆(TFSA)]^2?} _n fragments consisting of cyclic water tetramers with appended water molecules and TFSA anionic units, and 1-D metal–water tape [Cu–H₂O?···?(H₂O)₆?···?H₂O^?] _n in 2 were found. 3-D supramolecular networks of the two complexes consist of cationic complexes and water–TFSA anionic assemblies connected by hydrogen bonds.     

Intermolecular complexes of nido-C<Subscript>2</Subscript>B<Subscript>3</Subscript>H<Subscript>7</Subscript> with HF and LiH molecules: the theoretical studies,bonding properties and natural bond orbital (NBO) analysis

The changes in stabilization energy upon the formation of intermolecular hydrogen, dihydrogen and lithium bond complexes between C₂B₃H₇, LiH and HF have been investigated using MP2 method with aug-cc-pVDZ basis set. The interaction of HF with nido-C₂B₃H₇ could occur through the formation of B–H^?δ···^+δH–F, C–H^+δ···F^?δ–H and B–C···H–F classical and non-classical hydrogen bonds. The B–C bonds in backbone of the C₂B₃H₇ as electron donor interact with σ* orbital of HF as electron acceptor. Also interaction of LiH with nido-C₂B₃H₇ resulted in B–C···Li–H and B–H···LiH lithium bonds as well as C–H···H–Li dihydrogen bond complexes. In some of these complexes, LiH interacts with B–C bonds. Results are indicating that more stable complexes belong to interaction of HF and LiH with backbone of the nido-C₂B₃H₇. The AIM and NBO methods were used to analyze the intermolecular interactions; also the electron density at the bond critical point and the charge transfer of obtained complexes were studied.     

Competition between H<Subscript>2</Subscript>O and CH<Subscript>3</Subscript>OH molecules in the formation of the simplest stable proton disolvates and their solvation in aqueous methanol solutions of KOH and CH<Subscript>3</Subscript>OK

IR spectroscopic and quantum chemical methods are used to study the competition between water and methanol molecules in the formation of the simplest stable proton disolvates and their subsequent solvation in the case of solutions of KOH in CH₃OH and CH₃OK in H₂O with similar stoichiometries (~1:3-3.5). The complexes found in these solutions are analysed to determine their composition and structure: they are found to be heteroions (CH₃O?H?OH)^– solvated by two similar solvent molecules. In both cases, there are virtually no complexes of the second possible type (CH₃OH·(CH₃O?H?OCH₃)^–··H₂O or CH₃OH·(HO?H?OH)^–·H₂O), which appears to be due to the stoichiometric compositions of the solutions. It is shown that a DFT calculation (B3LYP/6-31++G(d,p)) of linear complexes with strong (~15-30 kcal/mol) H bonds reproduces, with good accuracy, the IR spectra of the solutions, which consist mainly of these complexes.     

Sequential bond energies of water to sodium glycine cation

Absolute bond dissociation energies of water to sodium glycine cations and glycine to hydrated sodium cations are determined experimentally by competitive collision-induced dissociation (CID) of Na⁺Gly(H₂O)_x, x = 1–4, with xenon in a guided ion beam tandem mass spectrometer. The cross sections for CID are analyzed to account for unimolecular decay rates, internal energy of reactant ions, multiple ion–molecule collisions, and competition between reaction channels. Experimental results show that the binding energies of water and glycine to the complexes decrease monotonically with increasing number of water molecules. Ab initio calculations at four different levels show good agreement with the experimental bond energies of water to Na⁺Gly(H₂O)_x, x = 0–3, and glycine to Na⁺(H₂O), whereas the bond energies of glycine to Na⁺(H₂O)_x, x = 2–4, are systematically higher than the experimental values. These discrepancies may provide some evidence that these Na⁺Gly(H₂O)_x complexes are trapped in excited state conformers. Both experimental and theoretical results indicate that the sodiated glycine complexes are in their nonzwitterionic forms when solvated by up to four water molecules. The primary binding site for Na⁺ changes from chelation at the amino nitrogen and carbonyl oxygen of glycine for x = 0 and 1 to binding at the C terminus of glycine for x = 2–4. The present characterization of the structures upon sequential hydration indicates that the stability of the zwitterionic form of amino acids in solution is a consequence of being able to solvate all charge centers.     

Ionization of sodium in water clusters

Structures of Na(H₂O)_n and Na⁺(H₂O)_n clusters with n = 1?6, 19, and 28 are determined in the second order of the Møller-Plesset perturbation theory with the use of extended atomic basis set 6–31++G^**. It is found that when the number of molecules is sufficient for the formation of two solvation shells around sodium, a continuous hydrogen-bond network is formed in both neutral and charged clusters, and the orientation of each molecule is determined by the balance between interactions with the neighboring water molecules and with the field of the central particle. In the cations, this field is stronger, and up to the third solvation shell, molecules have a predominant orientation with respect to sodium. In the neutral clusters, with an increase in the number of water molecules, the maximum of the electron density distribution of the highest occupied molecular orbital becomes more distant from the sodium nucleus, being shifted toward the cluster surface. The energy of this orbital accordingly decreases in absolute value approaching 22 kcal/mol inmicroparticles. In the charged clusters, the distribution of the positive charge generally correlates with the character of the highest occupied orbital in the neutral systems, so that with an increase in the number of molecules, the atomic charge of sodium decreases and tends to zero as n → ∞. The ionization potential of sodium changes in inverse proportion to the linear size of the cluster, and should not exceed 1.1 eV in watermicroparticles.     

Determination of the acid dissociation constants of p-benzohydroquinone by the INDO method

The stepwise acid dissociation constants for p-benzohydroquinone (QH₂) in aqueous media have been explicitly calculated for the first time, with the INDO parametrized SCF –MO method. We have optimized the geometries of QH₂, QH^?, and Q^2? and of the QH₂ · 6H₂O, QH^? · (H₃O⁺) · 5H₂O, and Q^2? · (H₃O⁺)₂ · 4H₂O systems that model the solvated species. The presence of the associated water molecules (and hydronium ions) account for the stabilization due to hydrogen bonding as well as for a part of the effect of interaction of these molecules with the respective reaction fields in an aqueous medium. To simulate the first solvation shell in a more complete manner, four more water molecules have been considered to be placed above and below the quinonoid ring and the optimized geometries of the resulting hydrated species, QH₂ · 10H₂O, QH^? · (H₃O⁺) · 9H₂O, and QH^? · (H₃O⁺) · 8H₂O, have been determined. The standard free-energy changes calculated for the dissociation of QH₂ into QH^? and H⁺ is 0.0251 Hartree (65.9 kJ mol^?1) and that of QH^? into Q^2? and H⁺ is 0.0285 Hartree (74.8 kJ mol^?1). Experimentally observed dissociation constants for these two steps correspond to free-energy changes of 0.0214 Hartree (56.2 kJ mol^?1) and 0.0248 Hartree (65.1 kJ mol^?1), respectively. © 1995 John Wiley & Sons, Inc.     

A Monte Carlo simulation of Fe2+ aqueous solvation

A Monte Carlo simulation of Fe²⁺ aqueous solvation, at 298 K, including 100 water molecules, has been done using periodic boundary conditions under the minimum image conversion. The energy has been calculated in the pair-potential approach, employing the MCY potential for the H₂OH₂O interaction and an ab initio analytical potential generated by us for the Fe²⁺H₂O interaction. The examination of interaction energies and of the radial distribution functions clearly show that the first hydration shell is formed by eight water molecules. By classifying the generated configurations into different significant structures of the solvent, it has been found that the eight water molecules of the first hydration shell are situated in a lightly distorted D_4d structure which maximizes the water—solute stabilization and minimizes the water—water repulsion. Finally, the validity of our theoretical predictions is discussed.     

Crystal Structures and Selected Properties of CoII Containing Thiostannates prepared by a New Room Temperature Route

The room temperature reaction of Na₄Sn₂S₆ · 5H₂O with CoCl₂ · 6H₂O and 2-(aminomethyl)pyridine (2-AMP) or trans-1,2-diamino-cyclohexane (DACH) leads to crystallization of two compounds with the compositions [Co(2-(aminomethyl)pyridine)₃]₂ Sn₂S₆ · 10H₂O ( 1 ) and [Co(trans-1,2-diaminocyclohexane)₃]₂Sn₂S₆ · 8H₂O ( 2 ). In both compounds [Sn₂S₆]^4– anions are present that are charge balanced each by two Co²⁺ centered complexes. Each of the two Co^II cations are sixfold coordinated by six N atoms of three 2-AMP or DACH ligands within slightly distorted octahedra. In compound 1 , the two complexes are linked by one [Sn₂S₆]^4– anion via strong N–H ··· S hydrogen bonds into centrosymmetric charge neutral trimeric units, that are further linked by weak C–H ··· S and N–H ··· S hydrogen bonds into chains that are directed along the a axis. These chains are further joined by N–H ··· O and O–H ··· O hydrogen bonds into a 3D network, with the H₂O molecules forming chains along the b axis. The crystal structure of 2 is similar to that of 1 featuring trimeric units which are also linked into chains. Between the chains water molecules are embedded that link the chains into a 3D network. Upon heating 2 in a thermobalance the water and ligand molecules are removed in discrete steps, indicating that compounds with more condensed thiostannate networks will form.     

Charge transfer and curve crossings in the [BeH2O]2+ system

Summary Diabatic and adiabatic potential energy curves have been determined for the complexation of beryllium cation with a water molecule, by means of multi-reference perturbation CI. The quasi-diabatic states correspond to Be²⁺H₂O and to nine charge transfer states (Be⁺H₂O⁺): at short beryllium-water distances the ground state is essentially Be²⁺H₂O, but at large distances several charge transfer states have lower energies. The nature of the curve crossings of the ground and lowest excited states in the [BeH₂O]²⁺ system is clarified. The changes brought about by the presence of a second water molecule are investigated.     

Stereochemical applications of mass spectrometry. 3—Energy dependence of the fragmentation of stereoisomeric methylcyclohexanols

Breakdown graphs have been constructed from charge exchange data for the epimeric 2-methyl-, 3-methyl- and 4-methyl-cyclohexanols. Although the breakdown graphs for epimeric pairs are essentially identical above ～12 eV recombination energy, significant differences are observed for the epimeric 2-methyl- and 4-methyl-cyclohexanols at low internal energies. For the 2-methylcyclohexanols the ratio ([M? H₂O]^+·/[M]^+·)_cis/([M? H₂O]^+·/[M]^+·)_trans is 3.2 in the [C₆F₆]^+· charge exchange mass spectra. This is attributed to both energetic and conformational effects which favour the stereospecific cis-1,4-H₂O elimination for the cis epimer. The breakdown graph for trans-4-methylcyclohexanol shows a sharp peak in the abundance of the [M? H₂O]^+· ion at ～10 eV recombination energy which is absent from the breakdown graph for the cis epimer. This peak is attributed to the stereospecific cis-1,4-elimination of water from the molecular ion of the trans isomer; the reaction appears to have a low critical energy but a very unfavourable frequency factor, and alternative modes of water loss common to both epimers are observed at higher energies. As a result, in the [C₆F₆]^+· charge exchange mass spectra the ([M? H₂O]^+·/[M]^+·)_trans/([M? H₂O]^+·/[M]^+·)_cis ratio is ～24, compared to the value of 13 observed in the 70 eV EI mass spectra. No differences are observed in either the metastable ion abundances or the associated kinetic energy releases for epimeric molecules.     

Remarkable shifts of Csp2-H and O-H stretching frequencies and stability of complexes of formic acid with formaldehydes and thioformaldehydes

Thirty-six stable complexes of formic acid with formaldehydes and thioformaldehydes were determined on the potential energy surface, in which the XCHO···HCOOH complexes are found to be more stable than the XCHS···HCOOH counterparts, with X = H, F, Cl, Br, CH₃, NH₂. All complexes are stabilized by hydrogen bonds, and their contribution to the total stabilization energy of the complexes increases in going from C-H···S to C-H···O to O-H···S and finally to O-H···O. Remarkably, a significant blueshift of C_sp2-H bond by 81–96 cm⁻¹ in the C_sp2-H···O hydrogen bond has hardly ever been reported, and a considerable redshift of O-H stretching frequency by 206–544 cm⁻¹ in the O-H···O/S hydrogen bonds is also predicted. The obtained results in our present work and previous literatures support that a distance contraction and a stretching frequency blueshift of C-H bond involving hydrogen bond depend mainly on its polarity and gas phase basicity of proton acceptor, besides the rearrangement of electron density due to complex formation. Markedly, we suggest the ratio of deprotonation enthalpy to proton affinity (R _c) as an indicator to prospect for classification of hydrogen bonds. The symmetry adapted perturbation theory results show a larger role of attractive electrostatic term in XO-n as compared to that in XS-n and the electrostatic interaction is overwhelming dispersion or induction counterparts in stabilizing XO-n and XS-n , with n = 1, 2, 3. © 2019 Wiley Periodicals, Inc.     

A theoretical investigation of the interaction between substituted carbonyl derivatives and water: Open or cyclic complexes?

The structures and binding energies of complexes between substituted carbonyl bases and water are the B3LYP/6‐311++G(d,p) computational level. The calculations also include the proton affinity (PA) of the O of the C?O group, the deprotonation enthalpies (DPE) of the CH bonds along a natural bond orbital analysis. The calculations reveal that stable open C?O···H_wO_w as well as cyclic CH···O_wH_w···O?C complexes are formed. The binding energies for the open complexes are linearly related to the PAs, whereas the binding energies for the cyclic complexes depend on both the PA and DPE. Different indicators of hydrogen bonds strength such as electron charge density, intramolecular and intermolecular hyperconjugation energy, occupation of orbitals, and charge transfer show significant differences between open and cyclic complexes. The contraction of the CH bond of the formyl group and the corresponding blue shift of the ν(CH) vibration are explained by the classical trans lone pair effect. In contrast, the elongation or contraction of the CH₃ group involved in the interaction with water results from the variation of the orbital interaction energies from the σ(CH) bonding orbital to the σ* and π* antibonding orbitals of the C?O group. The resulting blue or red shifts of the ν(CH₃) vibrations are calculated in the partially deuterated isotopomers. © 2012 Wiley Periodicals, Inc.     

Hydration of GABA and of the GABA·ZN2+ complex: A Monte Carlo simulation study

With empirical and theoretical atom–atom potentials the GABA·nH₂O, n = 25, 192 and GABA·Zn²⁺ · nH₂O, n = 25, 50, 100 complexes are simulated at 298.15 K by the Monte Carlo technique. The results show that the carboxyl group of GABA coordinates six water molecules. Two geometries of the GABA·Zn²⁺ complex, corresponding to the “direct” and “through-water” interaction of Zn²⁺ with the carboxyl group of GABA were found. For the latter interaction a GABA·Zn²⁺ · 6H₂O complex was found whereas the hydration of the former interaction leads to a GABA·Zn²⁺ · 5H₂O complex. Here the carboxyl group of GABA displaces only one water molecule in the first hydration shell of Zn²⁺. Energetically the “direct” and “through-water” geometries seem to be competitive, the former being slightly favored.     

Intermolecular interactions in uracil–nitrous acid complexes: structures,binding energy,topological properties,and nuclear magnetic resonance study

Molecular interactions between uracil and nitrous acid (U–NA) [C₄N₂O₂H₄? NO₂H] have been studied using B3LYP, B3PW91, and MP2 methods with different basis sets. The optimized geometries, harmonic vibrational frequencies, charge transfer, topological properties of electron density, nucleus‐independent chemical shift (NICS), and nuclear magnetic resonance one‐ and two‐bonds spin–spin coupling constants were calculated for U–NA complexes. In interaction between U and NA, eight cyclic complexes were obtained with two intermolecular hydrogen bonds N(C)HU…N(O) and OH_NA…O_U. In these complexes, uracil (U) simultaneously acts as proton acceptor and proton donor. The most stable complexes labeled, UNA1 and UNA2, are formed via NH bond of U with highest acidity and CO group of U with lowest proton affinity. There is a relationship between hydrogen bond distances and the corresponding frequency shifts. The solvent effect on complexes stability was examined using B3LYP method with the aug‐cc‐pVDZ basis set by applying the polarizable continuum model (PCM). The binding energies in the gas phase have also been compared with solvation energies computed using the PCM. Natural bond orbital analysis shows that in all complexes, the charge transfer takes place from U to NA. The results predict that the Lone Pair (LP)(O)_U → σ*(O? H) and LP(N(O)_NA → σ*(N(C)? H)_U donor–acceptor interactions are most important interactions in these complexes. Atom in molecule analysis confirms that hydrogen bond contacts are electrostatic in nature and covalent nature of proton donor groups decreases upon complexation. The relationship between spin–spin coupling constant (^1hJ_H…_Y and ^2hJ_H…_Y) with interaction energy and electronic density at corresponding hydrogen bond critical points and H‐bonds distances are investigated. NICS used for indicating of aromaticity of U ring upon complexation. © 2013 Wiley Periodicals, Inc.     

A simple model of solvation within the molecular orbital theory

A simple model of solvation within the molecular orbital method is proposed whereby the effect of solvent molecules is simulated by the inclusion of fractional point charges at the solvent atomic centers. The method is applied to three solvation problems: the hydration of Li⁺ and F^? and the solvation effect on the interaction between NH₃ and HF. The results of the first two calculations indicate that the point charge model is capable of reliably predicting solvation energies. The calculations for H₃N···HF demonstrate that the hydration has a profound effect on the potential energy surface favoring a proton transfer structure H₃NH⁺···F^?.     

Zum Komplexbildungsverhalten der 5,6-Dihalogen-7-oxa-bicyclo[2.2.1]heptan-2,3-dicarbonsäure gegenüber 3d-Übergangsmetallen

Complex Formation of 5,6-Dihalogeno-7-oxa-bicyclo[2.2.1]heptane-2,3-dicarboxylic Acid with 3d Transition Elements Carboxylate complexes of bivalent manganese, cobalt, nickel and copper with 5,6-dichloro- and 5,6-dibromo-7-oxa-bicyclo[2.2.1]heptane-2,3-dicarboxylic acid ( 3 and 4 ) have been prepared. For cobalt and nickel two types of complexes are formed: [ML^3/4(H₂O)₃] · H₂O and [ML^3/4(H₂O)₂], the latter is thermodynamically more stable. Manganese and copper form only complexes [MnL^3/4] and [CuL^3/4(H₂O)₂], respectively. The stereochemical configuration of the compounds have been deduced from their spectroscopic and magnetic properties. The metal atoms have been found to be in an octahedral environment. The stability constants of the complexes have been determined by potentiometric measurements. The thermal decomposition of the complexes has been studied by thermogravimetry and differential thermal analysis. The complexes of 3 are thermally more stable than the corresponding ones of 4 . The X-ray structure analysis of [CoL³(H₂O)₃] · H₂O shows a monomeric structure of the complex within the crystal and an octahedral coordination of the metal ion. The dicarboxylate anion acts as a tridentate ligand, the other octahedral sites are occupied by three water molecules. The chlorine atoms are not involved in the network of hydrogen bonds within the crystal packing.     

Structures and energetics of complexation of metal ions with ammonia,water, and benzene: A computational study

Quantum chemical calculations have been performed at CCSD(T)/def2‐TZVP level to investigate the strength and nature of interactions of ammonia (NH₃), water (H₂O), and benzene (C₆H₆) with various metal ions and validated with the available experimental results. For all the considered metal ions, a preference for C₆H₆ is observed for dicationic ions whereas the monocationic ions prefer to bind with NH₃. Density Functional Theory–Symmetry Adapted Perturbation Theory (DFT‐SAPT) analysis has been employed at PBE0AC/def2‐TZVP level on these complexes (closed shell), to understand the various energy terms contributing to binding energy (BE). The DFT‐SAPT result shows that for the metal ion complexes with H₂O electrostatic component is the major contributor to the BE whereas, for C₆H₆ complexes polarization component is dominant, except in the case of alkali metal ion complexes. However, in case of NH₃ complexes, electrostatic component is dominant for s‐block metal ions, whereas, for the d and p‐block metal ion complexes both electrostatic and polarization components are important. The geometry (M⁺–N and M⁺–O distance for NH₃ and H₂O complexes respectively, and cation–π distance for C₆H₆ complexes) for the alkali and alkaline earth metal ion complexes increases down the group. Natural population analysis performed on NH₃, H₂O, and C₆H₆ complexes shows that the charge transfer to metal ions is higher in case of C₆H₆ complexes. © 2016 Wiley Periodicals, Inc.