Theoretical study of hydrogen bonding interaction in nitroxyl (HNO) dimer: interrelationship of the two N-H...O blue-shifting hydrogen bonds

The hydrogen bonding interactions of the HNO dimer have been investigated using ab initio molecular orbital and density functional theory (DFT) with the 6-311++G(2d,2p) basis set. The natural bond orbital (NBO) analysis and atom in molecules (AIM) theory were applied to understand the nature of the interactions. The interrelationship between one N-H...O hydrogen bond and the other N-H...O hydrogen bond has been established by performing partial optimizations. The dimer is stabilized by the N-H...O hydrogen bonding interactions, which lead to the contractions of N-H bonds as well as the characteristic blue-shifts of the stretching vibrational frequencies nu(N-H). The NBO analysis shows that both rehybridization and electron density redistribution contribute to the large blue-shifts of the N-H stretching frequencies. A quantitative correlations of the intermolecular distance H...O (r(H...O)) with the parameters: rho at bond critical points (BCPs), s-characters of N atoms in N-H bonds, electron densities in the sigma*(N-H), the blue-shift degrees of nu(N-H) are presented. The relationship between the difference of rho (|Deltarho|) for the one hydrogen bond compared with the other one and the difference of interaction energy (DeltaE) are also illustrated. It indicates that for r(H...O) ranging from 2.05 to 2.3528 A, with increasing r(H...O), there is the descending tendency for one rho(H...O) and the ascending tendency for the other rho(H...O). r(H...O) ranging from 2.3528 to 2.85 A, there are descending tendencies for the two rho(H...O) with increasing r(H...O). On the potential energy surface of the dimer, the smaller the difference between one rho(H...O) and the other rho(H...O) is, the more stable the structure is. As r(H...O) increases, the blue-shift degrees of nu(N-H) decrease. The cooperative descending tendencies in s-characters of two N atoms with increasing r(H...O) contribute to the decreases in blue-shift degrees of nu(N-H). Ranging from 2.05 to 2.55 A, the increase of the electron density in one sigma*(N-H) with elongating r(H...O) weakens the blue-shift degrees of nu(N-H), simultaneously, the decrease of the electron density in the other sigma*(N-H) with elongating r(H...O) strengthens the blue-shift degrees of nu(N-H). Ranging from 2.55 to 2.85 A, the cooperative ascending tendencies of the electron densities in two sigma*(N-H) with increasing r(H...O) contribute to the decreases in blue-shift degrees of nu(N-H).     

Bifurcated hydrogen-bonding effect on the shielding and coupling constants in trifluoroacetyl pyrroles as studied by 1H, 13C and 15N NMR spectroscopy and DFT calculations

According to the (1)H, (13)C and (15)N NMR spectroscopic data and DFT calculations, bifurcated N--H...N and N--H...O intramolecular hydrogen bond is shown to be present in 2-trifluoroacetyl-5-(2'-pyridyl)-pyrrole. This bifurcated hydrogen bond causes an increase in the absolute size of the (1)J(N,H) coupling constant by about 6 Hz, and the deshielding of the bridge proton by 2 ppm. DFT calculations show that the influence of the N--H...N and N--H...O intramolecular hydrogen bonds on the (1)J(N,H) coupling and proton shielding is almost additive, although the components of the bifurcated hydrogen bond slightly weaken each other. In 2-trifluoroacetyl-5-(2'-pyridyl)-pyrrole, the coupling constants involving the fluorine and the N--H covalent bond nuclei depend dramatically on the spatial position of the pyridine ring. The pyridine ring rotation operates as a quantum switch controlling the spin information transfer between the (19)F and (15)N nuclei, as well as the proton.     

Quantum-chemical study of CHCl3-SO2 association

CHCl(3)-SO(2) association is studied by high-level quantum-chemical calculations of stationary points of the dimer electronic potential-energy hypersurface, including correlated second-order Moller-Plesset and CCSD(T) calculations with basis sets up to 6-311++G(d,p). During geometry optimization, frequency, and energy calculations, a self-written computer code embedding the GAMESS ab initio program suite applies counterpoise correction of the basis set superposition error. A CH...O hydrogen-bonded complex (DeltaE(0)=-8.73 kJmol) with a 2.4 A intermolecular H...O distance and two very weak van der Waals complexes (DeltaE(0)=-3.78 and -2.94 kJmol) are located on the counterpoise-corrected potential-energy surface. The intermolecular interactions are characterized by Kitaura-Morokuma interaction energy decompositions and Mulliken electron population analyses. The unusual hydrogen bond is distinguished by a CH-bond contraction, a pronounced enhancement of the IR intensity and a shift to higher frequency ("blueshift") of the CH-stretching vibration compared to the CHCl(3) monomer. Spectroscopy and association in liquid solution is also discussed; our results provide an alternative explanation for features in the CH-stretching vibration spectrum of chloroform dissolved in liquid sulfur dioxide which have been attributed previously to an intermolecular Fermi resonance.     

Magnitudes and chemical consequences of R(3)N(+)-C-H...O[double bond]C hydrogen bonding

The magnitude of the stabilizing interaction between an aliphatic C[bond]H bond attached to an ammonium nitrogen and a carbonyl oxygen was evaluated by ab initio calculations at the MP2/6-311++G** level of theory. Attractive R(3)N(+)-C-H...O[double bond]C interactions play an important role in supramolecular recognition and various types of stereoselective catalysis. Our calculations show that R(3)N(+)-C-H...O[double bond]C is the strongest hydrogen bond of the C-H...O type known to date. Such hydrogen bonds remain as stabilizing interactions even in water for amide acceptors.     

A theoretical prediction of stability in hydrogen-bonded complexes formed between oxirane and oxetane rings with HX (X=F and Cl)

The optimised geometries of heterocyclic hydrogen-bonded complexes, C2H4O...HX and C3H6O...HX, where X=F or Cl, were determined at DFT/B3LYP/6-311++G(d,p) computational level. Structural, electronic and vibrational properties of these complexes are used in order to compare the strained ring, which confer the great reactivity of these heterocyclic rings with monoprotic acids, forming a primary hydrogen bond. A secondary hydrogen bond between the hydrogen atoms of the CH2 groups and the halide species also takes place, thus causing a nonlinearity (characterized by the theta angle), in the primary hydrogen bond.     

Electronic properties of multifurcated bent hydrogen bonds CH3...Y and CH2...Y

H-bonding angle angleYHX has an important effect on the electronic properties of the H-bond Y...HX, such as intra- and intermolecular hyperconjugations and rehybridization, and topological properties of electron density. We studied the multifurcated bent H-bonds of the proton donors H3CZ (Z = F, Cl, Br), H2CO and H2CF2 with the proton acceptors Cl(-) and Br(-) at the four high levels of theory: MP2/6-311++G(d,p), MP2/6-311++G(2df,2p), MP2/6-311++G(3df,3pd) and QCISD/6-311++G(d,p), and found that they are all blue-shifted. These complexes have large interaction energies, 7-12 kcal mol(-1), and large blue shifts, delta r(HC) = -0.0025 --0.006 A and delta v(HC) = 30-90 cm(-1). The natural bond orbital analysis shows that the blue shifts of these H-bonds Y...HnCZ are mainly caused by three factors: rehybridization; indirect intermolecular hyperconjugation n(Y) -->sigma*(CZ), in that the electron density from n(Y) of the proton acceptor is transferred not to sigma*(CH), but to sigma*(CZ) of the donor; intramolecular hyperconjugation n(Z) -->sigma*(CH), in that the electron density in sigma*(CH) comes back to n(Z) of the donor such that the occupancy in sigma*(CH) decreases. The topological properties of the electron density of the bifurcated H-bonds Y...H2CZ are similar to those of the usual linear H-bonds, there is a bond critical point between Y and each hydrogen, and a ring critical point inside the tetragon YHCH. However, the topological properties of electron density of the trifurcated H-bonds Y...H3CZ are essentially different from those of linear H-bonds, in that the intermolecular bond critical point, which represents a closed-shell interaction, is not between Y and hydrogen, but between Y and carbon.     

气相中O3与HSO自由基间的氢键复合物

气相中O3与HSO自由基之间的相互作用及其反应在大气化学中非常重要. 在DFT-B3LYP/6-311++G**和MP2/6-311++G**水平上求得O3+HSO复合物势能面上的稳定构型, B3LYP方法得到了三种构型(复合物I, II和III), 而MP2方法只能得到一种构型(复合物II). 在复合物I和III中, HSO单元中的1H原子作为质子供体, 与O3分子中的端基O原子作为质子受体相互作用, 形成红移氢键复合物; 而在复合物II中, 虽与复合物I和III中具有相同的质子供体和质子受体, 却形成了蓝移氢键复合物. B3LYP/6-311++G**水平上计算的单体间相互作用能的计算考虑了基组重叠误差(BSSE)和零点振动能(ZPVE)校正, 其值在-3.37到-4.55 kJ·mol-1之间. 采用自然键轨道理论(NBO)对单体间相互作用的本质进行了考查, 并通过分子中原子理论(AIM)分析了三种复合物中氢键的电子密度拓扑性质.     

pi-H...O hydrogen bonds: multicenter covalent pi-H interaction acts as the proton-donating system

The MP2 method and the Pople-style basis sets 6-311++G(d,p), 6-311++G(2df,2pd), and 6-311++G(3df,3pd) were used to perform calculations on H3O+...C2H2 and C2H3+...C2H2 complexes and related species. Hydrogen bonds existing for the analyzed complexes were investigated as well as related pi-H...O --> pi...H-O and pi-H...pi --> pi...H-pi proton-transfer processes. For some of the complexes analyzed the multicenter pi-H interaction possessing the properties of a covalent bond acts as a proton donor; more generally it is classified as the Lewis acid. The quantum theory of "atoms in molecules" (QTAIM) was also applied to deepen the nature of these interactions in terms of characteristics of bond critical points. The pi-H...O, O-H...pi, and pi-H...pi interactions analyzed here may be classified as hydrogen bonds since their characteristics are the same as or at least similar to those of typical hydrogen bonds. H...pi interactions are common in crystal structures of organic and organometallic compounds. The analyses performed here show a continuum of such interactions since there are H...pi contacts possessing the characteristics of weak intermolecular interactions on the one hand and pi-H multicenter covalent bonds on the other. Ab initio and QTAIM results support the latter statements.     

Structural studies of metabolic products of dopamine. III. Crystal and molecular structure of (--)-adrenaline.

The crystal structure of (--)-adrenaline has been determined by X-ray methods, using 831 observed reflections collected by counter methods. The crystals are monoclinic, space group P2-1 with a=7.873(2), b=6.790(2), c=8.638(2) A and beta=98.01(2) degrees. Least-squares refinements yielded a conventional R-factor of 0.053. Standard deviations in bond lengths are 0.005-0.006 A and in bond lengths aree 0.005-0.006 A and in bond angles 0.4 degrees. The adrenaline molecules were found to exist as zwitterions in the crystals. The conformation of the adrenaline molecule corresponds closely to that usually encountered among the salts of the sympathomimetic amines. The crystals consist of molecular double layers parallel to (100). The molecules within a layer are linked through hydrogen bonds of the types N--H...O and O--H...O, whereas the layers are connected by van der Waals interactions.     

Reinvestigation of intramolecular hydrogen bond in malonaldehyde derivatives: An ab initio,AIM and NBO study

The RAHB systems in malonaldehyde and its derivatives at MP2/ 6‐311++G(d,p) level of theory were studied and their intramolecular hydrogen bond energies by using the related rotamers method was obtained. The topological properties of electron density distribution in O? H···O intramolecular hydrogen bond have been analyzed in term of quantum theory of atoms in molecules (QTAIM). Correlations between the H‐bond strength and topological parameters are probed. The results of QTAIM clearly showed that the linear correlation between the electron density distribution at HB critical point and RAHB ring critical point with the corresponding hydrogen bond energies was obtained. Moreover, it was found a linear correlation between the electronic potential energy density, V(r_cp), and hydrogen bond energy which can be used as a simple equation for evaluation of HB energy in complex RAHB systems. Finally, the similar linear treatment between the geometrical parameters, such as O···O or O? H distance, and Lp(O)→σ*_OH charge transfer energy with the intramolecular hydrogen bond energy is observed. © 2010 Wiley Periodicals, Inc., Int J Quantum Chem, 2011     

乙醇-水分子团簇C2H5OH(H2O)n (n=1-9)稳定结构的量子化学研究

采用密度泛函理论B3LYP方法, 在B3LYP/6-311++G(2d,2p)//B3LYP/6-311++G(d,p)基组水平上对乙醇-水分子团簇(C₂H₅OH(H₂O)_n (n=1-9))的各种性质进行研究, 如: 优化的几何构型、结构参数、氢键、结合能、平均氢键强度、自然键轨道(NBO)电荷分布、团簇的生长规律等. 结果表明, 从二维(2-D)环状结构到三维(3-D)笼状结构的过渡出现在n=5的乙醇-水分子团簇中. 此外, 利用团簇结合能的二阶差分、形成能、能隙等性质, 发现在n=6时乙醇-水分子团簇的最低能量结构稳定性较好, 可能为幻数结构. 最后, 为了进一步探讨氢键本质, 将C₂H₅OH(H₂O)_n (n=2-9)最低能量结构的各种性质与纯水分子团簇(H₂O)_n (n=3-10)比较, 结果表明前者与后者中的水分子之间氢键相似.     

Electron density topological analysis of hydrogen bonding in glucopyranose and hydrated glucopyranose

Topological analysis of the electron density profiles and the atomic basin integration data for the most energetically favorable (4)C(1) and (1)C(4) conformers of beta-D-glucopyranose, calculated at the B3LYP/6-31+G(d), MPWlPW91/6-311+G(2d,p), and MP2/6-31+G(d) levels, demonstrates that intramolecular hydrogen bonding between adjacent ring OH groups does not occur in glucopyranose, given the need to demonstrate a bond critical point (BCP) of correct (3,-1) topology for such an interaction to be termed a hydrogen bond. On the other hand, pyranose ring OH groups separated by three, rather than just two, carbon atoms are able to form an intramolecular hydrogen bond similar in topological properties and geometry to that found for propane-1,3-diol. Vicinal, equatorial OH groups in the (4)C(1) conformer of glucopyranose are, however, able to form strong bidentate hydrogen bonds with water molecules in a cooperative manner, each water molecule acting simultaneously as both hydrogen bond donor and acceptor, and characterized by (3,-1) bond critical points with increased values for the electron density and the Laplacian of rho(r) compared to an isolated ethane-1,2-diol/water complex.     

Bifurcate hydrogen bonds. Interaction of intramolecularly H-bonded systems with Lewis bases

The structure and the hydrogen bonding in the systems formed by the intramolecularly H-bonded systems, namely, maltol (3-hydroxy-2-methyl-4-pyrone), 5, 2,4,6-trinitrophenol, 6, acetylacetone enol, 7, with Lewis bases, phosgene, 8, dioxane, 9, and DMSO, 10, have been studied by density functional theory (B3LYP) and MP2 using the 6-311G* basis set. The continuum solvent effect was simulated by IEF-PCM model. The hydrogen bond analysis using the atoms in molecules (AIM) method was applied by using the MP2(full)/6-311++G** electron density to establish the nature of the bifurcate hydrogen bond (BHB) in these systems as well as contributory factors for its stabilization. The nature of interaction in the intermolecular H-complexes formed by compounds 5- 7 with the Lewis bases 8- 10 was shown to depend on the strength of the intramolecular hydrogen bond O...H and the strength of the base. The critical values of the CO...H and NO...H angles for which the formation of BHB is possible, have been determined.     

Topological characteristics of electron density distribution in SSXY --> XSSY (X or Y = F, Cl, Br, I) isomerization reactions

The reactions of SSXY to XSSY (X or Y = F, Cl, Br, I) have been studied at B3LYP/6-311++G(2df) and MP2/6-311++G(2df) levels based on B3LYP/6-311++G(2df) optimized geometries. There are two pathways (X or Y atom transferring) for each reaction of SSXY to XSSY(X or Y = F, Cl, Br, I). The "Atoms in Molecules" theory has been applied to analyze the topological characteristics of electron density distribution along the reaction path. (a) The S-S-X or S-S-Y ring structure transition region (STR) is contributed to S(1)-X or S(1)-Y bond formation and S(2)-X or S(2)-Y bond annihilation. The STR and structure transition state (STS) are defined. (b) Comparing the two pathways of each SSXY --> XSSY (X, Y = F, Cl, Br, I) reaction, X or Y atom transferring, the broader the STR is, the later the STS appears, and the pathway is easier. (c) When X linked to the same S site of a three-member ring (S-S-F, S-S-Cl, S-S-Br, or S-S-I), the ring STR becomes broader and the STS appears later according to the sequence of X = F, Cl, Br, and I. And in these cases, for exothermic reactions, the higher the DeltaE(P-->TS1) is, the broader the STR is and the later the STS appears. (d) When the same Y atom linked to different three-member-ring S-S-X, the STR becomes broader and the STS appears later according to the sequence of X = F, Cl, Br, and I in the ring. And in these cases, for exothermic reactions, the lower the DeltaE(P-->TS2) is, the broader the STR is and the later the STS appears.     

Weakly bound complexes of N2O: an ab initio theoretical analysis toward the design of N2O receptors

Ab initio calculations at MP2/6-311++G(2d,2p) and MP2/6-311++G(3df,3pd) computational levels have been used to analyze the interactions between nitrous oxide and a series of small and large molecules that act simultaneously as hydrogen bond donors and electron donors. The basis set superposition error (BSSE) and zero point energy (ZPE) corrected binding energies of small N2O complexes (H2O, NH3, HOOH, HOO*, HONH2, HCO2H, H2CO, HCONH2, H2CNH, HC(NH)NH2, SH2, H2CS, HCSOH, HCSNH2) vary between -0.93 and -2.90 kcal/mol at MP2/6-311++G(3df,3pd) level, and for eight large complexes of N2O they vary between -2.98 and -3.37 kcal/mol at the MP2/6-311++G(2d,2p) level. The most strongly bound among small N2O complexes (HCSNH2-N2O) contains a NH..N bond, along with S-->N interactions, and the most unstable (H2S-N2O) contains just S-->N interactions. The electron density properties have been analyzed within the atoms in molecules (AIM) methodology. Results of the present study open a window into the nature of the interactions between N2O with other molecular moieties and open the possibility to design N2O abiotic receptors.     

AIM analysis of intramolecular hydrogen bonding in O-hydroxy aryl Schiff bases

AIM analysis was applied to study the changes in such topological parameters as the electron density at critical points of all the bonds of the molecule during the so-called nonadiabatic proton transfer in intramolecular hydrogen bonding in o-hydroxy aryl Schiff bases. Proton transfer is presented by a stepwise elongation and fixing of the hydroxyl bond with complete optimization of the rest of the parameters of the molecule by the B3LYP/6-311++G(d,p) method. A more detailed study of electron density changes at the critical points of the chelate and phenol rings in the stepwise proton-transfer process is presented. It was shown that the dependency of the electron density at the critical point of the chelate ring on tautomeric equilibrium is of a complicated character, whereas it is linear for the phenol ring. A complex study of the changes in the total electron density at the hydrogen bond, the quasi-aromatic ring, and in the whole molecule has been accomlished. The calculations of the intramolecular hydrogen bond by means of conformational and topological methods are discussed.     

Estimation on the intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies in glycine and alanine peptides

Computation of accurate intramolecular hydrogen-bonding energies for peptides is of great importance in understanding the conformational stabilities of peptides and developing a more accurate force field for proteins. We have proposed a method to determine the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies in glycine and alanine peptides. In this article, the method is further applied to evaluate the intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies in peptides. The optimal structures of the intramolecular 10-membered ring N-H...O=C hydrogen bonds in glycine and alanine tripetide molecules are obtained at the MP2 level with 6-31G(d), 6-31G(d,p), and 6-31+G(d,p) basis sets. The intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies are then evaluated based on our method at the MP2/6-311++G(3df,2p) level with basis set superposition error correction. The intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies are calculated to be in the range of -6.84 to -7.66, -4.44 to -4.98, and -6.95 to -7.88 kcal/mol. The method is also applied to estimate the individual intermolecular hydrogen-bonding energies in the dimers of amino-acetaldehyde, 2-amino-acetamide, formamide, and oxalamide, each dimer having two identical intermolecular hydrogen bonds. According to our method, the individual intermolecular hydrogen-bonding energies in the four dimers are calculated to be -1.77, -1.67, -6.35, and -4.82 kcal/mol at the MP2/6-311++G(d,p) level, which are in good agreement with the values of -1.84, -1.72, -6.23, and -4.93 kcal/mol predicted by the supermolecular method.