Assessment of intermolecular interactions at three sites of the arylalkyne in phenylacetylene‐containing lithium‐bonded complexes: Ab initio and QTAIM studies

The intermolecular interactions existing at three different sites between phenylacetylene and LiX (X = OH, NH₂, F, Cl, Br, CN, NC) have been investigated by means of second‐order Møller?Plesset perturbation theory (MP2) calculations and quantum theory of “atoms in molecules” (QTAIM) studies. At each site, the lithium‐bonding interactions with electron‐withdrawing groups (? F, ? Cl, ? Br, ? CN, ? NC) were found to be stronger than those with electron‐donating groups (? OH and ? NH₂). Molecular graphs of C₆H₅C?CH···LiF and πC₆H₅C?CH···LiF show the same connectional positions, and the electron densities at the lithium bond critical points (BCPs) of the πC₆H₅C?CH···LiF complexes are distinctly higher than those of the σC₆H₅C?CH···LiF complexes, indicating that the intermolecular interactions in the C₆H₅C?CH···LiX complexes can be mainly attributed to the π‐type interaction. QTAIM studies have shown that these lithium‐bond interactions display the characteristics of “closed‐shell” noncovalent interactions, and the molecular formation density difference indicates that electron transfer plays an important role in the formation of the lithium bond. For each site, linear relationships have been found between the topological properties at the BCP (the electron density ρ_b, its Laplacian ?²ρ_b, and the eigenvalue λ₃ of the Hessian matrix) and the lithium bond length d(Li‐bond). The shorter the lithium bond length d(Li‐bond), the larger ρ_b, and the stronger the π···Li bond. The shorter d(Li‐bond), the larger ?²ρ_b, and the greater the electrostatic character of the π···Li bond. © 2012 Wiley Periodicals, Inc.     

Ab initio study of the structure,cooperativity, and vibrational properties in the mixed hydrogen‐bonded trimers of hydrogen isocyanide and water

The five trimers of H₂O···HNC···H₂O, H₂O···H₂O···HNC, HNC···H₂O···H₂O, H₂O···HNC···HNC, and HNC···HNC···H₂O have been studied with quantum chemical calculations. Their structures, harmonic vibrational frequencies and interaction energies have been calculated at the B3LYP and MP2 levels with the aug‐cc‐pVDZ and aug‐cc‐pVTZ basis sets. The cooperative effect on these properties has also been studied quantitatively. For HNC:(H₂O)₂ systems, the cyclic H₂O···H₂O···HNC trimer is most stable with an interaction energy of ?16.01 kcal/mol and a large cooperative energy of ?3.25 kcal/mol at the MP2/aug‐cc‐pVTZ level. For H₂O:(HNC)₂ systems, the interaction energy and cooperative energy in the H₂O···HNC···HNC trimer are larger than those in the HNC···HNC···H₂O trimer. The NH stretch frequency has a blue shift for the terminal HNC molecule in the HNC···H₂O···H₂O and HNC···HNC···H₂O trimers and a red shift in other cases. A many‐body analysis has also been performed to understand the interaction energies in these hydrogen‐bonded clusters. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Density functional studies on hydrogen‐bonded clusters of hydrogen halides and the interaction on halide anions

Density functional theory (DFT) calculations have been performed to study the structures and stability of X^?·(HX)_n=2–5 clusters where X = F, Cl, Br at B3LYP/6‐311++G** level of theory. The presence of halide ions in these clusters disintegrates the hydrogen halide clusters. All the hydrogen halides are then hydrogen bonded to the centrally placed halide ions, thereby forming multiple hydrogen bonds. The interaction energies have been corrected for the basis set superposition error (BSSE) using Boy's counterpoise correction method. Evidence for the destruction of hydrogen bonds in hydrogen halide clusters due to the presence of halide ions is further obtained from topological analysis and natural bond orbital analysis. The chemical hardness and chemical potential have been calculated for all the anion clusters. The above analysis reveals that hydrogen bonding in these systems is not an essentially electrostatic interaction. The nature of the stabilization interactions operative in these multiple hydrogen‐bonded clusters has been explained in terms of many‐body contribution to interaction energies. From these studies, an attempt has been made to understand the nature of the molecular properties resulting from different electronegativities of the halogens. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Ab initio and AIM studies on typical π‐type and pseudo‐π‐type halogen bonds: Comparison with hydrogen bonds

Series of typical π‐type and pseudo‐π‐type halogen‐bonded complexes B ··· ClY and B ··· BrY and hydrogen‐bonded complex B ··· HY (B = C₂H₄, C₂H₂, and C₃H₆; Y = F, Cl, and Br) have been investigated using the MP2/aug‐cc‐pVDZ method. A striking parallelism was found in the geometries, vibrational frequencies, binding energies, and topological properties between B ··· XY and B ··· HY (X = Cl and Br). It has been found that the lengths of the weak bond d(X ··· π)/d(H ··· π), the frequencies of the weak bond ν(X ··· π)/ν(H ··· π), the frequency shifts Δν(X? Y)/Δν(H? Y), the electron densities at the bond critical point of the weak bonds ρ_c(X ··· π)/ρ_c(H ··· π), and the electron density changes Δρ_c(X? Y)/Δρ_c(H? Y) could be used as measures of the strengths of typical π‐type and pseudo‐π‐type halogen/hydrogen bonds. The typical π‐type and pseudo‐π‐type halogen bond and hydrogen bond are noncovalent interactions. For the same Y, the halogen bond strengths are in the order B ··· ClY < B ··· BrY. For the same X, the halogen bond strength decreases according to the sequence F > Cl > Br that is in agreement with the hydrogen bond strengths B ··· HF > B ··· HCl > B ··· HBr. All of these typical π‐type and pseudo‐π‐type hydrogen‐bonded and halogen‐bonded complexes have the “conflict‐type” structure. Contour maps of the Laplacian of π electron density indicate that the formation of B ··· XY halogen‐bonded complex and B ··· HY hydrogen‐bonded complex is very similar. Charge transfer is observed from B to XY/HY and both the dipolar polarization and the volume of the halogen atom or hydrogen atom decrease on B ··· XY/B ··· HY complex formation. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

硝酸乙酯分子间相互作用的ab initio研究

在abinitio-HF/6-31G水平上求得硝酸乙酯二聚体势能面上的四种优化构型和电子结构。经MP2电子相关校正和基组叠加误差(BSSE)以及零点能(ZPE)校正，求得二聚体的最大结合能为11.46kJ.mol^-^1，还进行HF/6-311G和HF/6-311++G水平的总能量比较计算，发现6-31G基组对计算结合能比较适合，二子体系间的电荷转移很少，对优化构型进行振动分析，并基于统计热力学求得从单体形成二聚体的热力学性质变化。     

Ab initio studies on the polarizability of lithium clusters: Some unusual results

The polarizability of lithium clusters was studied using ab initio as well as Density Functional Theory (DFT)‐based methods. The electron correlation was found not to play any major role in the evaluation of polarizability of the lithium clusters, contrary to the case of sodium and other types of clusters. The theoretical methods used in the present study are not able to reproduce the unusual features in the size evolution of the polarizability observed experimentally. An interesting feature is that the polarizability of lithium clusters calculated by various methods at 0K temperature is found to be higher than the experimental results measured at higher temperature, which is quite unusual. The significance of these issues is highlighted and the need of highly accurate DFT exchange‐correlation functionals and ab initio methods in the study of the electronic properties of lithium clusters is emphasized. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Ab initio prediction of the vibrational spectra of the hydrogen‐bonded complexes between CO and HONO2

The vibrational characteristics (vibrational frequencies and infrared intensities) for free and complexed CO and HONO₂ have been predicted using ab initio calculations at SCF and MP2 levels with different basis sets and B3LYP/6?31G(d,p) calculations. The ab initio calculations show that the complexation between HONO₂ and CO leads to two stable structures: CO … HONO₂ (1A) and OC … HONO₂ (1B). The changes in the vibrational characteristics from free monomers to complexes have been estimated. It was established that the most sensitive to the complexation is the stretching O? H vibration. In agreement with the experiment, its vibrational frequency in the complexes is shifted to lower frequency (Δν = ?123 cm^?1). The magnitude of the wave number shift is indicative of relatively strong hydrogen‐bonded interaction. The ab initio calculations at different levels predict an increase of the infrared intensity of the stretching O? H vibration for structure 1A more than five times and for structure 1B more than nine times. The most consistent agreement between the computed values of the frequency shifts for structure 1B and those experimentally observed suggests that this structure is preferred. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Computational study of structures and proton transfer in hydrogen‐bonded ammonia complexes using semiempirical valence‐bond approach

In this study, the seGVB method was implemented for the N H bonding system, specifically for hydrogen‐bonded ammonia complexes, and the model well reproduces the MP2 geometries and energetics. A comparison between the ammonia dimer and water dimer is given from the viewpoint of valance‐bond structures in terms of the calculated bond energies and pair–pair interactions. The linear hydrogen bond is found to be stronger than the bent bonds in both cases, with the difference in energy between the linear and cyclic structures being comparable in both cases although the NH bonds are generally weaker. The energy decomposition clearly demonstrates that the changes in electronic energy are quite different in the two cases due to the presence of an additional lone pair on the water molecule, and it is this effect which leads to the net stabilization of the cyclic structure for the ammonia dimer. Proton‐transfer profiles for hydrogen‐bonded ammonia complexes [NH₂ H NH₂]⁻ and [NH₃ H NH₃]⁺ were calculated. The barrier for proton transfer in [NH₃ H NH₃]⁺ is larger than that in [NH₂ H NH₂]⁻, but smaller than that in the protonated water dimer. The different bonding structures substantially affect the barrier to proton transfer, even though they are isoelectronic systems. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 73: 357–367, 1999     

Spectroscopic and theoretical study of vibrational spectra of hydrogen‐bonded tropolone

Theoretical simulation of the bandshape and fine structure of the ν_s stretching band is presented for tropolone‐H and tropolone‐D taking into account an adiabatic coupling between the high‐frequency O–H(D) stretching and the low‐frequency intra‐ and intermolecular OO stretching modes, and linear and quadratic distortions of the potential energies for the low‐frequency vibrations in the excited state of the O–H(D) stretching vibration. In order to determine the low‐frequency vibrations, the experimental spectra of the polycrystalline tropolone in the far‐infrared and the low‐frequency Raman range have been recorded for the first time. The experimental frequencies in the low‐frequency region are compared with the results of the HF/6‐31G** and Becke3LYP/6‐31G** calculations carried out for the tropolone dimer. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 73: 275–282, 1999     

Theoretical study of the red‐ and blue‐shifted hydrogen bonds of nitroxyl and acetylene dimers

Ab initio molecular orbital and density functional theory (DFT) in conjunction with different basis sets calculations were performed to study the C? H…O red‐shifted and N? H…π blue‐shifted hydrogen bonds in HNO? C₂H₂ dimers. The geometric structures, vibrational frequencies and interaction energies were calculated by both standard and counterpoise (CP)‐corrected methods. In addition, the G3B3 method was employed to calculate the interaction energies. The topological and natural bond orbital (NBO) analysis were investigated the origin of N? H…π blue‐shifted hydrogen bond. From the NBO analysis, the electron density decrease in the σ* (N? H) is due to the significant electron density redistribution effect. The blue shifts of the N? H stretching frequency are attributed to a cooperative effect between the rehybridization and electron density redistribution. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

New basis sets for the evaluation of interaction‐induced electric properties in hydrogen‐bonded complexes

Interaction‐induced static electric properties, that is, dipole moment, polarizability, and first hyperpolarizability, of the CO? (HF)_n and N₂? (HF)_n, n = 1–9 hydrogen‐bonded complexes are evaluated within the finite field approach using the Hartree–Fock, density functional theory, Møller–Plesset second‐order perturbation theory, and coupled cluster methods, and the LPol‐n (n = ds, dl, fs, fl) basis sets. To compare the performance of the different methods with respect to the increase of the complex size, we consider as model systems linear chains of the complexes. We analyze the results in terms of the many‐body and cooperative effects. © 2012 Wiley Periodicals, Inc.     

Correlation between hydrogen‐bonding interaction and mechanical properties of polyimide fibers

Novel co‐polymerization polyimide (PI) fibers based on 4,4′‐oxydianiline (ODA)‐pyromellitic dianhydride (PMDA) were prepared. 2‐(4‐Aminophenyl)‐5‐aminobenzimidazole (PABZ) containing the N? H group was introduced into the structure of the fibers as the proton donor. The results of Fourier transform infrared (FTIR) and dynamic mechanical analysis (DMA) showed that hydrogen bonding occured between the N? H group and chains, which strongly enhanced interchain interaction. This hydrogen bonding interaction increased the tensile strength and initial modulus of the PI fibers up to 2.5 times and 26 times, respectively, compared to those of homo‐PI PMDA‐ODA fibers with no hydrogen‐bonding interaction because of the absence of proton donors after the imidization process. In the mean time, glass transition temperature (T_g) of the modified PI fibers was found to be 410–440°C, which was higher than that of the homo‐PI PMDA‐ODA fibers. From the result, a novel access to molecular design and manufacture of high performance PI fibers with good properties could be provided. Copyright © 2009 John Wiley & Sons, Ltd.     

Single‐electron halogen bond: Ab initio study

Ab initio calculations have been performed on single‐electron halogen bonds between methyl radical and bromine‐containing molecules to gain a deeper insight into the nature of such noncovalent interactions. Bader's atoms in molecules (AIM) theory have also been applied to the analysis of the linking of the single‐electron halogen bond. Various characteristics of the R? Br…CH₃ interaction, i.e., binding energies, geometrical parameters and topological properties of the electron density have been determined. The presence of the bond critical points (BCPs) between the bromine atom and methyl radical and the values of electron density and Laplacian of electron density at these BCPs indicate the closed‐shell interactions in the complexes. The single‐electron halogen bonds, which are significantly weaker than the normal halogen bonds, exhibit equally bond strength as compared to the single‐electron hydrogen bond. It has been also found that plotting of the binding energies versus topological properties of the electron density at the BCPs gives two straight lines. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

A simple model of hydrogen bonding with particular application to trends in hydrogen‐bonded dimers

A simple model has been proposed to explain trends in the computed interaction energy, bond length changes, frequency shifts and infrared intensities for the chlorofluoromethanes CF_nCl_mH, FH and FArH on complexation with the isoelectronic diatomics BF, CO, N₂ and the rare gas atoms Kr, Ar, Ne to form a series of linear or nearly linear hydrogen‐bonded complexes. The dipole moment derivative of the proton donor (with respect to the stretching coordinate) and the chemical hardness of the hydrogen‐bonded atom of the proton acceptor are identified as two useful parameters for rationalizing the changes in some of the molecular properties of the proton donor when the hydrogen bond is formed. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Density functional theory study of nitrogen‐14 nuclear quadrupole coupling parameters of L‐histidine: hydrogen‐bonded system

The calculations of nitrogen‐14 nuclear quadrupole parameters, nuclear quadrupole coupling constant, χ, and asymmetry parameter, η, of L‐His were done in two distinct environments: one as a free fully optimized molecule, an isolated molecule with the geometrical parameters taken from X‐ray, and the other in the orthorhombic and monoclinic solid states. The most probable interacting molecules with the central molecule in the crystalline phase were considered in the hexameric clusters to include hydrogen‐bonding effects in the calculations. The computations were performed with PW91P86/6‐31++G** and B3LYP6‐31++G** methods using the Gaussian 98 program. The good agreement between the nitrogen‐14 quadrupole parameters of the free His and imidazole molecules with their microwave available data demonstrates that the applied level of theory and the 6‐31++G** basis set are suitable to obtain reliable electric field gradient values. In the solid state, the shifts of quadrupole coupling parameters from the monomer to the solid phase are reasonably well reproduced for the amino and imino sites of imidazole ring in a hexameric cluster. That implies the fact that the hexameric cluster worked effectively to generate the results which are compatible with the experiment. The quadrupole coupling constant values of –N⁺H₃ group are in fair agreement with the experiment. This discrepancy is due to the absences of vibrational effects and the rotation of –N⁺H₃ group around N–C(α) bond. Copyright © 2012 John Wiley & Sons, Ltd.     

Density functional theory study on hydrogen‐bonded complexes of adenine with polyformamide molecules

The structures of hydrogen‐bonded complexes A–F_n (n = 2–7) of adenine with polyformamide molecules have been fully optimized at B3LYP/6‐31G(d) basis set level. All the formamide molecules prefer to be N? H proton donor rather than C? H proton donor and are favorably bound to the five‐numbered moiety of adenine. A displacement of formamide molecules to one side of adenine mean plane has happened with an increasing number of formamide molecules. An obvious effect of hydrogen‐bonding cooperativity can be seen during the complex process. The most interesting geometrical change of adenine upon the complex is the shortening of the bond C4? N6 resulting from the strengthening of the conjugation between the π system of the adenine ring and the lone pair of the nitrogen atom. An existence of weak N? H···π bonding interaction between the π system of adenine and N? H bond of F7 is found and further conformed by an natural bond orbital analysis specially carried out on A–F₇. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Vibrational dynamics of hydrogen‐bonded HCN complexes with OH and NH acids: Computational DFT systematic study

Vibrational properties (band position, infrared [IR], and Raman intensities) of C?N stretching mode were studied in 65 gas phase hydrogen‐bonded 1:1 complexes of HCN with OH acids and NH acids using density functional theory (DFT) calculations at the B3LYP‐6‐311++G(d,p) level. Furthermore, general characteristics of the hydrogen bonds and vibrational changes in acids OH/NH stretching bands were also considered. Experimentally observed blue shift of the C?N stretching band promoted by hydrogen bonding, which shortens the triple bond length, is very well reproduced and quantitatively depends on the hydrogen bond length. Both IR and Raman ν(C?N) band intensities are enhanced, also in good agreement with the experimental results. IR intensity increase is a direct function of the hydrogen bond energy. However, the predicted Raman intensity raise is a more complex function, depending simultaneously on characteristics of both the hydrogen bond (C?N bond length) and the H‐donating acid (polarizability). With these two parameters, ν (C?N) Raman intensities of the complexes are explained with a mean error of ±2.4%. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

C4H6O,C4H6S与XF (X＝F,Cl, Br)分子间卤键的电子密度拓扑研究

运用MP2/aug-cc-pVDZ方法对2,5-二氢呋喃, 2,5-二氢噻吩与XF (X＝F, Cl, Br)之间的卤键作用进行了理论研究. 研究发现: C₄H₆O, C₄H₆S与XF之间不仅存在O(S)…XF n型卤键, C＝C双键与XF分子亦可形成π型卤键|对于C₄H₆O与XF之间的n型和π型卤键以及C₄H₆S与XF之间的π型卤键, 卤键键能ΔE、键鞍点处的电子密度ρ(r_c)以及电子给体到受体之间的电子转移数Δq(XF)均按B…F₂＜B…ClF＜B…BrF (B＝C₄H₆O, C₄H₆S)的顺序依次增大|对于卤键键能较大的体系C₄H₆O…BrF(n), C₄H₆O…BrF(π), C₄H₆S…F₂(n), C₄H₆S…ClF(n), C₄H₆S…BrF(n), C₄H₆S…BrF(π), 卤键作用介于离子键和共价键之间|而对于其它的卤键键能较小的体系, 卤键作用为闭壳层静电作用.     

Theoretical study of hydrogen‐bonded complexes of benzene with hydrides of astrochemical interest

Post Hartree–Fock and DFT calculations have been performed for studying the possibility for a benzene support to be linked to various hydrides through a quasi Bz···H? A bond. Interaction energy of compounds, including C? H bonds (CH₄, CH₃F, CH₂O, CHN, CHN? O), N? H bonds (NH₃, NH₂F, NHC, NHCO, NH₃O), and O? H bonds (OH₂, OHF, NCOH), were evaluated, taking basis set superposition error (BSSE) and zero point vibrational energy (ZPVE) corrections into account. Numerical convergence of results with respect to the ingredients included at different steps of theory (basis set, DFT functionals, correlation treatments, geometry optimization) was tested mainly on the example of the water adduct and, for comparison, the Bz···H₃O⁺ system containing a cation instead of a neutral molecule. A rather large range of adsorption energies is obtained, from about 1 kcal/mol for methane to more than 6 kcal/mol for cyanic acid, according to the acidic character of the adsorbed species in each family of Bz···H? A bonds. Some consequences for astrophysical problems involving PAHs in the interstellar medium are pointed out. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008