Molecular understanding of the adhesive interactions between silica surface and epoxy resin: Effects of interfacial water

The molecular mechanism of the adhesion between silica surface and epoxy resin under atmospheric conditions is investigated by periodic density-functional-theory (DFT) calculations. Slab models of the adhesion interface were built by integrating a fragment of epoxy resin and hydroxylated (0 0 1) surface of α-cristobalite in the presence of adsorbed water molecules. Effects of adsorbed water on the adhesion interaction are evaluated on the basis of geometry-optimized structures, adhesion energies, and forces. Calculated results demonstrate that adsorbed water molecules significantly reduce both the adhesion energies and forces of the silica surface–epoxy resin interface. The reduction of adhesion properties can be associated with structural deformation of water molecules confined in the tight space between the adhesive and adherend as well as structural flexibility of the hydrogen-bonding network in the interfacial region during detachment of the epoxy resin from the hydrophilic silica surface. © 2018 Wiley Periodicals, Inc.     

A computational study of interactions between acetic acid and water molecules

Density functional theory calculations were performed for the title reactions to elucidate the difference between the strong cyclic hydrogen bond of (Me-COOH)(2) and the electrolytic dissociation, MeCOOH <==> Me-COO(-) + H(+), as a weak acid. The association of water clusters with acetic acid dimers strengthens the cyclic hydrogen bond. A nucleophilic attack of the carboxylic carbon by a water cluster leads to a first zwitterionic intermediate, MeCOO(-) + H(3)O(+) + (HO)(3)C-Me. The intermediate is unstable and is isomerized to a neutral interacting system, MeCOOH...(HO)(3)C-Me + H(2)O. The ethanetriol, (HO)(3)-CMe is transformed to an acetic acid monomer. The monomer may be dissociated to give a second zwitterionic intermediate with reasonable proton-relay patterns and energy changes. In proton relay reaction channels, H in MeCOOH is not an acidic proton but is always a hydroxy proton.     

Comprehensive study of the interaction between hydrogen halides and methanol derivatives

The R? CH₂? HO…H? X (R = SCl, Cl, SH, NO₂, OMe, CHO, CN, C₂H₅, CH₃, H; X = F, Cl, Br) complexes are considered here as the interest sample for the consideration of different measures of H‐bond strength. The intermolecular interaction energies are predicted by using MP2/6‐31++G(d,p) and B3LYP/6‐31++G(d,p) methods with basis set superposition error and zero‐point energy corrections. The results showed that intermolecular hydrogen bonds for complexes with HF are stronger than such interactions in complexes with HCl and HBr. Quantum theory of “Atoms in Molecules” and natural bond orbitals method were applied to analyzed H‐bond interactions. The gas phase thermodynamic properties of complexes were predicted using quantum mechanical computations. The obtained results showed a strong influence of the R and X substituents on the thermodynamic properties of complexes. Numerous correlations between topological, geometrical, thermodynamic properties and energetic parameters were also found. © 2011 Wiley Periodicals, Inc.     

Coupled hydrogen-bonding interactions between beta-carboline derivatives and acetic acid

In this work, we have analysed the tendency of two beta-carboline derivatives, harmane and norharmane, in the formation of hydrogen bonds. We obtained the (1)H and (13)C NMR spectra of different mixtures of these derivatives with acetic acid (AcOH) in CDCl(3). A cyclic 1:3 complex is proposed between harmane and AcOH, while a 1:2 complex is proposed for norharmane. Chemical shifts at temperatures between 233 and 323 K were measured: lowering the temperature produces the same effect as increasing the amount of AcOH in solution. The (13)C data confirm a delocalisation of the pi electron density towards the pyridinic ring that occurs when AcOH is added.     

Dynamics of a protein and water molecules surrounding the protein: hydrogen-bonding between vibrating water molecules and a fluctuating protein

Internal and rigid-body motions of bovine pancreatic trypsin inhibitor (BPTI) and of water molecules surrounding the BPTI are studied in a vicinity of an energy minimum using a normal mode analysis proposed as the independent molecule model. Water's rigid-body motion is predominant in comparison to its internal motions. We have derived information about the relationship between the magnitude of a thermal ellipsoid of an H-bonding atom and the anisotropy of its ellipsoid, and the relationship between the magnitude of the ellipsoid and the H-bond strength. We see a relationship between vibrational frequencies (assuming rigid-body motion of the water molecules) and the H-bond strength of the water taking part in this H-bonding. Analyzing the H-bond strength, we found that a hydrogen in water is likely to H-bond to oxygen in the protein, whereas an oxygen in water has a less strong preference to H-bond to the protein. For water molecules acting as the hydrogen acceptor, strong H-bonding has longer lifetimes than weak H-bonding.     

Density functional theory and topological analysis on the hydrogen bonding interactions in cysteine‐thymine complexes

Hydrogen bonding interactions between amino acids and nucleic acid bases constitute the most important interactions responsible for the specificity of protein binding. In this study, complexes formed by hydrogen bonding interactions between cysteine and thymine have been studied by density functional theory. The relevant geometries, energies, and IR characteristics of hydrogen bonds (H‐bonds) have been systematically investigated. The quantum theory of atoms in molecule and natural bond orbital analysis have also been applied to understand the nature of the hydrogen bonding interactions in complexes. More than 10 kinds of H‐bonds including intra‐ and intermolecular H‐bonds have been found in complexes. Most of intermolecular H‐bonds involve O (or N) atom as H‐acceptor, whereas the H‐bonds involving C or S atom usually are weaker than other ones. Both the strength of H‐bonds and the structural deformation are responsible for the stability of complexes. Because of the serious deformation, the complex involving the strongest H‐bond is not the most stable structures. Relationships between H‐bond length (ΔR_X‐H), frequency shifts (Δv), and the electron density (ρ_b) and its Laplace (?²ρ_b) at bond critical points have also been investigated. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Hydrogen bonding interactions between ethylene glycol and water: density,excess molar volume,and spectral study

Studies of the density and the excess molar volume of ethylene glycol (EG)-water mixtures were carried out to illustrate the hydrogen bonding interactions of EG with water at different temperatures. The re-sults suggest that a likely complex of 3 ethylene glycol molecules bonding with 4 water molecules in an ethylene glycol-water mixture (EGW) is formed at the maximal excess molar volume,which displays stronger absorption capabilities for SO2 when the concentration of SO2 reaches 400×10?6 (volume ratio) in the gas phase. Meanwhile,FTIR and UV spectra of EGWs were recorded at various EG concentra-tions to display the hydrogen bonding interactions of EG with water. The FTIR spectra show that the stretching vibrational band of hydroxyl in the EGWs shifts to a lower frequency and the bending vibra-tional band of water shifts to a higher frequency with increasing the EG concentration,respectively. Furthermore,the UV spectra show that the electron transferring band of the hydroxyl oxygen in EG shows red shift with increasing the EG concentration. The frequency shifts in FTIR spectra and the shifts of absorption bands in UV absorption spectra of EGWs are interpreted as the strong hydrogen bonding interactions of the hydrogen atoms in water with the hydroxyl oxygen atoms of EG.     

Theoretical study of the effects of substitution,solvation, and structure on the interaction between nitriles and methanol

We report an investigation on intermolecular interactions in R? CN ··· H? OCH₃ (R = H, CH₃, F, Cl, NO₂, OH, SH, SCH₃, CHO, COCH₃, CH₂Cl, CH₂F, CH₂OH, CH₂COOH, CF₃, SCOCH₃, SCF₃, OCHF₂, CH₂CF₃, CH₂OCH₃, and CH₂CH₂OH) complexes using density functional theory. The calculations were conducted on B3LYP/6‐311++G** level of theory for optimization of geometries of complexes and monomers. An improper hydrogen bonding (HB) in the H₃CO? H ··· NC? R complexes was observed in that N atom of the nitriles functions acts as a proton acceptor. Furthermore, quantum theory of “Atoms in Molecules” (AIM) and natural bond orbital (NBO) method were applied to analyze H‐bond interactions in respective complexes. The electron density (ρ) and Laplacian (?²ρ) properties, estimated by atoms in molecules calculations, indicate that H ··· N bond possesses low ρ and positive ?²ρ values, which are in agreement with partially covalent character of the HBs, whereas O? H bonds have negative ?²ρ values. In addition, the weak intermolecular force due to dipole–dipole interaction (U) is also considered for analysis. The examination of HB in these complexes by quantum theory of NBO method fairly supports the ab initio results. Natural population analysis data, the electron density, and Laplacian properties, as well as, the ν(O? H) and γ(O? H) frequencies of complexes, calculated at the B3LYP/6‐311++G** level of theory, are used to evaluate the HB interactions. The calculated geometrical parameters and conformational analysis in water phase solution show that the H₃CO? H ··· NC? R complexes in water are more stable than that in gas phase. The obtained results demonstrated a strong influence of the R substituent on the properties of complexes. Numerous correlations between topological, geometrical, thermodynamic properties, and energetic parameters were also found. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

A new method for estimating the importance of hydrogen-bonding groups in the binding site of a protein

We introduce a new method to estimate the importance of hydrogen-bonding sitepoints in the binding site of a protein as part of a structure-based design strategy. Our method identifies hydrogen-bonding sitepoints within a binding pocket and ranks them according to both the accessibility of their hydrogen bonding regions to incoming ligands and their hydrogen-bonding strength. The combination of these components produces a prioritised list of sitepoints that are more likely to be involved in hydrogen bonding with an incoming ligand. A dataset of known protein-ligand interactions was used to compare the prioritisation of sitepoints identified by our method with those observed to be engaged in hydrogen bonding in their crystal structures. Our method was able to remove those sitepoints unable to bind the ligand due to a low accessibility or an unfavourable orientation and to award significantly higher hydrogen-bonding ranking values to those sitepoints observed to form hydrogen bonds. Our method can thus be used to identify hydrogen-bonding sitepoints that should be targeted preferentially in a drug design strategy.     

Theoretical study on interactions between thiourea S‐monoxide and water

This work investigated the hydrolyzation of thiourea monoxide using density functional theory (DFT) and Møller–Plesset second‐order (MP2) Theory. We obtained the equilibrium structures and other molecular properties of the clusters. The results show that thiourea monoxide has a good solubility in water solvent, and as indicated, the binding energies of the clusters are increased progressively by the addition of water molecules. Furthermore, the increases of the distance(C? S) and (S? O) by the addition of water molecules indicate that the strength of the C? S and S? O bonds are weakened. When n = 7 (n, the number of water molecules), the C? S bonds of Clusters VII ruptured. We conclude that thiourea monoxide can be decomposed in aqueous solution. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Different hydrogen-bonding patterns in two [Zn(ENTPP)] complexes with water or methanol as ligands

We have synthesized zinc complexes of H₂ENTPP (5-(8-ethoxycarbonyl-1-naphthyl)-10,15,20-triphenyl porphyrin) as a model to study hydrogen-bonding interactions. When water or methanol is a ligand, crystals of [Zn(ENTPP)(CH₃OH)] or [Zn(ENTPP)(H₂O)]?·?C₆H₅CH₃ were obtained. In both structures, the ligand has hydrogen-bonding interactions, but in different patterns. In [Zn(ENTPP)(CH₃OH)], the methanol oxygen and carboxylate oxygen in the naphthyl group form an intermolecular hydrogen bond. In [Zn(ENTPP)(H₂O)]?·?C₆H₅CH₃, there are two independent molecules A and B. In molecule B, there is an intramolecular hydrogen bond between the water oxygen and the carboxylate oxygen, while in molecule A, besides the intramolecular hydrogen bond, there is an intermolecular hydrogen bond between the water oxygen and the carboxylate oxygen. ¹H NMR spectra suggest the binding of methanol or water to zinc are equilibrium processes in solution. Equilibrium constant has been determined by UV-Vis measurements, and it suggests the binding affinity of zinc to methanol has been moderately increased.     

A computational NQR study on the hydrogen-bonded lattice of cytosine-5-acetic acid

A computational study at the level of density functional theory (DFT) employing 6-311++G** standard basis set was carried out to evaluate nuclear quadrupole resonance (NQR) spectroscopy parameters in cytosine-5-acetic acid (C5AA). Since the electric field gradient (EFG) tensors are very sensitive to the electrostatic environment at the sites of quadruple nuclei, the most possible interacting molecules with the target one were considered in a five-molecule model system of C5AA using X-ray coordinates transforming. The hydrogen atoms positions were optimized and two model systems of original and H-optimized C5AA were considered in NQR calculations. The calculated EFG tensors at the sites of (17)O, (14)N, and (2)H nuclei were converted to their experimentally measurable parameters, quadrupole coupling constants and asymmetry parameters. The evaluated NQR parameters reveal that the nuclei in original and H-optimized systems contribute to different hydrogen bonding (HB) interaction. The comparison of calculated parameters between optimized isolated gas-phase and crystalline monomer also shows the relationship between the structural deformation and NQR parameters in C5AA. The basis set superposition error (BSSE) calculations yielded no significant errors for employed basis set in the evaluation of NQR parameters. All the calculations were performed by Gaussian 98 package of program.     

Hydrogen bonding interaction between 1,4‐dioxane and water

This work reports an interaction of 1,4‐dioxane with one, two, and three water molecules using the density functional theory method at B3LYP/6‐311++G* level. Different conformers were studied and the most stable conformer of 1,4‐dioxane‐(water)_n (n = 1–3) complex has total energies ?384.1964038, ?460.6570694, and ?537.1032381 hartrees with one, two, and three water molecules, respectively. Corresponding binding energy (BE) for these three most stable structures is 6.23, 16.73, and 18.11 kcal/mol. The hydrogen bonding results in red shift in O? O stretching and C? C stretching modes of 1,4‐dioxane for the most stable conformer of 1,4‐dioxane with one, two, and three water molecules whereas there was a blue shift in C? O symmetric stretching and C? O asymmetric stretching modes of 1,4‐dioxane. The hydrogen bonding results in large red shift in bending mode of water and large blue shift in symmetric stretching and asymmetric stretching mode of water. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

A theoretical study of the sulfenate–sulfoxide rearrangement. Effect of the hydrogen bond complexation

A DFT study of the thermal and radical sulfenate–sulfoxide rearrangement of derivatives of 3‐propenyl sulfoxide has been carried out. The effect of the substitution and hydrogen bond complexation has been analyzed. The results show that without external factors the radical breakdown path is the one preferred by the alkyl and aromatic derivatives while the unsubstituted system proceeds preferentially through a two‐step series of [1,3]‐ and [2,3]‐sigmatropic shifts. The inclusion of a hydrogen bond donor interacting with the oxygen atom increases the stability of all the species except the radical and the final products. Thus, in the dimethyl derivative the radical and two‐step processes present similar limiting steps. The analysis of the electron density of the systems provides some relationships between the properties at the bond critical point and the interatomic distances for the S···C and H···O cases. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110:2391–2397, 2010     

Halogen bond versus hydrogen bond: The many‐body interactions approach

The analysis of interrelation between halogen bond and hydrogen bond in the (RX)(HNC)(HCN) complexes (R = CH₃, CF₃ and X = Cl, Br, I) was performed on the basis of DFT calculations. Both two‐body additive contributions and three‐body nonadditive contributions to the total interaction energy were discussed. QTAIM was used for topological analysis of electron density. Additionally, QTAIM analysis of electron density was performed for both two‐ and three‐body complexes. The electron charge transfer in trimers showed the dual character of the fragment with halogen atom involved into the investigated interactions—it acts as Lewis acid and Lewis base, depending on the type of interaction considered. The effect of cooperativity of X‐ and H‐bonding was assessed on the basis of many‐body interaction energy and electron density analysis. Additionally, an alternative two‐body model with the same situation (in the context of intermolecular interactions) is investigated. The anti‐cooperative effect was found also for this model.     

Insights into the cooperativity between multiple interactions of dimethyl sulfoxide with carbon dioxide and water

Interactions of dimethyl sulfoxide with carbon dioxide and water molecules which induce 18 significantly stable complexes are thoroughly investigated. An addition of CO₂ or H₂O molecules into the DMSO⋯1CO₂ and DMSO⋯1H₂O systems leads to an increase in the stability of the resulting complexes, in which it is larger for a H₂O addition than a CO₂. The overall stabilization energy of the DMSO⋯1,2CO₂ is mainly contributed by the S=O⋯C Lewis acid–base interaction, whereas the O − H⋯O hydrogen bond plays a significant role in stabilizing complexes of DMSO⋯1,2H₂O and DMSO⋯1CO₂⋯1H₂O. Remarkably, the complexes of DMSO⋯2H₂O are found to be more stable than DMSO⋯1CO₂⋯1H₂O and DMSO⋯2CO₂. The level of the cooperativity of multiple interactions in ternary complexes tends to decrease in going from DMSO⋯2H₂O to DMSO⋯1CO₂⋯1H₂O and finally to DMSO⋯2CO₂. It is generally found that the red shift of the O − H bond involved in an O − H⋯O hydrogen bond increases while the blue shift of a C − H bond in a C − H⋯O hydrogen bond decreases when a cooperative effect occurs in ternary complexes as compared to those of the corresponding binary complexes. © 2018 Wiley Periodicals, Inc.     

The water molecule arrangement over the side chain of some aliphatic amino acids: A quantum chemical and bottom-up investigation

The geometry of surrounding water molecules on the side chain of glycine, alanine, α-aminoisobutyric acid, α-aminobutyric acid, valine, and related hydrocarbons has been analyzed combining bottom-up and quantum chemical methodologies. To minimize the cavity size and to prevent water-water hydrogen bonding loss, the water molecules adopt a shape, resembling the one found in crystal structure of gas clathrate hydrates, with water molecules tangentially oriented to the surface of hydrophobic side chain. The cage is directly hydrogen bonded to the backbone's polar groups, thus hydration shells around hydrophobic and hydrophilic groups are folded together in amphiphilic molecules. The hydrophobe enclathration implies a substantial freedom degree reduction which makes it entropically disfavored. This disadvantageous entropic contribution is partially compensated by the favorable van der Waals interactions with guest in stabilizing clathrate hydrate formation. The water shell around the side chain relates intimately with the side-chain rotational isomerism. Present data are correlated with the experimental determined populations of the three rotamers, yielding promising results for both α-aminobutyric acid and valine.     

DFT study on the hydrogen bonds of phenol–cyclohexanone and phenol–H2O2 in the Baeyer–Villiger oxidation

Hydrogen bonds of phenol–cyclohexanone and phenol–H₂O₂ in the studied Baeyer–Villiger (B–V) oxidation have been investigated by HF, B3LYP, and MP2 methods with various basis sets. The accurate single‐point energies were performed using CCSD(T)/6‐31+G(d,p) and CCSD(T)/aug‐cc‐pVDZ on the optimized geometries of MP2/6‐31+G(d,p). It has been confirmed that B3LYP/6‐31+G(d,p) could be used to study such hydrogen bonds. Energetic analysis of complexes was carried out using the Xantheas method with BSSE corrected by CP method. Orbital energy order (ε) illuminated that phenol with good hydrogen donor‐acceptor property can interact with cyclohexanone or H₂O₂ to form hydrogen bound complexes, and the binding energies (BE) range from ?4.38 to ?14.06 kcal mol^?1. NBO analysis indicated that the redistribution of atomic charges in the complexes facilitated nucleophilic attack of H₂O₂ on cyclohexanone. The calculated results match remarkably well with the experimental phenomena. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

氯化钙甲醇溶液中离子溶剂化作用的理论研究

利用红外光谱研究不同温度下CaCl2/甲醇溶液体系的溶剂化作用,结果表明在溶液中CaCl2以离子形式与甲醇发生溶剂化作用,且溶剂化数随温度升高而降低.通过密度泛函理论（DFT）在B3LYP/6-31G＊＊水平下对CaCl2/甲醇溶液中可能存在的配位构型进行结构优化及热力学性质的计算,说明了在CaCl2/甲醇溶液中各种配位构型存在的可能性,得出温度升高热力学数据的变化规律,解释了溶剂化数随温度升高而降低的趋势.进一步对各种可能配位构型的红外吸收频率进行计算并与实验结果进行比较,推断在CaCl2/甲醇溶液中主要存在的配位构型为[CaCl（CH3OH）n]＋和[Cl（CH3OH）n]-.     

The investigation on the relative stability of different clusters of thiourea dioxide in water using gas phase quantum chemical calculations

The relative stability of different clusters of thiourea dioxide (TDO) in water is examined using gas phase quantum chemical calculations at the MP2 and B3LYP level with 6‐311++G(d,p) basis set. The possible equilibrium structures and other energetic and geometrical data of the thiourea dioxide clusters, TDO‐(H₂O)_n (n is the number of water molecules), are obtained. The calculation results show that a strong interaction exists between thiourea dioxide and water molecules, as indicated by the binding energies of the TDO clusters progressively increased by adding water molecules. PCM model is used to investigate solvent effect of TDO. We obtained a negative hydration energy of ?20.6 kcal mol^?1 and free‐energy change of ?21.0 kcal mol^?1 in hydration process. On the basis of increasing binding energies with adding water molecules and a negative hydration energy by PCM calculation, we conclude thiourea dioxide can dissolve in water molecules. Furthermore, the increases of the C? S bond distance by the addition of water molecules show that the strength of the C? S bonds is attenuated. We find that when the number of water molecules was up to 5, the C? S bonds of the clusters, TDO‐(H₂O)₅ and TDO‐(H₂O)₆ were ruptured. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009