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     

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.     

Theoretical study of bifurcated hydrogen bonding effects on the 1J(N,H), 1hJ(N,H), 2hJ(N,N) couplings and 1H, 15N shieldings in model pyrroles

According to the density functional theory calculations, the X···H···N (X?N, O) intramolecular bifurcated (three‐centered) hydrogen bond with one hydrogen donor and two hydrogen acceptors causes a significant decrease of the ^1hJ(N,H) and ^2hJ(N,N) coupling constants across the N? H···N hydrogen bond and an increase of the ¹J(N,H) coupling constant across the N? H covalent bond in the 2,5‐disubsituted pyrroles. This occurs due to a weakening of the N? H···N hydrogen bridge resulting in a lengthening of the N···H distance and a decrease of the hydrogen bond angle at the bifurcated hydrogen bond formation. The gauge‐independent atomic orbital calculations of the shielding constants suggest that a weakening of the N? H···N hydrogen bridge in case of the three‐centered hydrogen bond yields a shielding of the bridge proton and deshielding of the acceptor nitrogen atom. The atoms‐in‐molecules analysis shows that an attenuation of the ^1hJ(N,H) and ^2hJ(N,N) couplings in the compounds with bifurcated hydrogen bond is connected with a decrease of the electron density ρ_H···N at the hydrogen bond critical point and Laplacian of this electron density ?²ρ_H···N. The natural bond orbital analysis suggests that the additional N? H···X interaction partly inhibits the charge transfer from the nitrogen lone pair to the σ*_{N? H} antibonding orbital across hydrogen bond weakening of the ^1hJ(N,H) and ^2hJ(N,N) trans‐hydrogen bond couplings through Fermi‐contact mechanism. An increase of the nitrogen s‐character percentage of the N? H bond in consequence of the bifurcated hydrogen bonding leads to an increase of the ¹J(N,H) coupling constant across the N? H covalent bond and deshielding of the hydrogen donor nitrogen atom. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

Hydrogen bonding interactions between <Emphasis Type="Italic">N</Emphasis>,<Emphasis Type="Italic">N</Emphasis>-dimethylformamide and cysteine: DFT studies of structures,properties, and topologies

The hydrogen bonding interactions between cysteine and N,N-dimethylformamide (DMF) were studied at the extended hybrid functional DFT-X3LYP/6-311++G(d,p) level regarding their geometries, energies, vibrational frequencies, and topological features of the electron density. The quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analyses were employed to elucidate the interaction characteristics in the complexes. The results show that two intermolecular hydrogen bonds (H-bonds) are formed in one complex except few complexes with one intermolecular H-bond. The H-bonds involving O atom of DMF as H-bond acceptor usually are red-shifting H-bonds, while the blue-shifting H-bond usually involve methyl of DMF or methenyl of cysteine moiety as H-bond donors. Both hydrogen bonding interaction and structural deformation play important roles in the relative stabilities of the complexes. Due to the π-bond cooperativity, the strongest H-bond is formed between hydroxyl of cysteine moiety and O atom of DMF, however, the serious deformation counteract the hydrogen bonding interaction to a great extent. The complex involves a stronger hydrogen bonding interaction as well as the smaller deformation is the most stable one. The electron density (ρ_b) as well as its Laplacian (∇²ρ_b) at the H-bond critical point predicted by QTAIM is strongly correlated with the H-bond structural parameter (δR _H···Y) and the second-perturbation energies E(2) in the NBO scheme.     

Synthesis and Theoretical Study of Intramolecular Hydrogen Bond at Two Possible Positions in Pyrazolo[1,2‐b]phthalazine

Properties of dimethyl 3‐(alkylamino)‐5,10‐dioxo‐5,10‐dihydro‐1H‐pyrazolo[1,2‐b]phthalazine‐1,2‐dicarboxylate and its derivatives were studied by means of ab initio method. NO₂ derivative of title compound was synthesized and the nature of its intramolecular hydrogen bond (HB) was investigated. Furthermore, the topological properties of the electron density distributions for N? H···O intramolecular bridges were analyzed in terms of the Bader theory of atoms in molecules (AIM). The electron density (ρ) and Laplacian (?²ρ) properties, estimated by AIM calculations, indicated that O···H bond possesses low ρ and positive ?²ρ values which are in agreement with electrostatic character of the HBs, whereas N? H bonds have covalent character (?²ρ<0). Moreover, steric effect of the t‐Bu group on structure and topological parameters of pyrazolo[1,2‐b]phthalazine conformers was studied. Finally, the powerful method of Espinosa was used to obtain the H‐bond energy.     

Hydrogen‐bonded complexes of nicotine with simple alcohols

Ab initio and density functional theory studies have been performed on the hydrogen‐bonded complexes of neutral and protonated nicotine with ethanol, methanol, and trifluromethanol to explore their relative stability in a systematic way. Among all the hydrogen‐bonded nicotine complexes considered here, protonated forms in nicotine–ethanol and nicotine–methanol, and neutral form in nicotine–trifluromethanol complexes have been found to be the most stable. In the former two complexes, the proton attached to the pyrrolidine nitrogen acts as a strong hydrogen bond donor, whereas the pyrrolidine nitrogen atom acts as a hydrogen bond acceptor in the latter case. Neutral complex of nicotine with trifluromethanol has been found to possess a very short hydrogen bond (1.57 Å) and basis set superposition error corrected hydrogen bond energy value of 19 kcal/mol. The nature of the various hydrogen bonds formed has been investigated through topological aspects using Bader's atoms in molecules theory. From the calculated topological results, excellent linear correlation is shown to exist among the hydrogen bond length, electron density, and its Laplacian at the bond critical points for all the complexes considered. The natural bond orbital analysis has been carried out to investigate the charge transfer in the nicotine alcohol complexes. In contrast to the blue shifting behavior that is generally exhibited by other C? H···O hydrogen bonds involving sp³ carbon atom, the C? H···O hydrogen bond in the protonated nicotine–ethanol and methanol complexes has been found to be proper with red shifting in nature. © 2011 Wiley Periodicals, Inc.     

A novel tubular hydrogen‐bond pattern in a new diazaphosphole oxide: a combination of X‐ray crystallography and theoretical study of hydrogen bonds

In the structure of 2‐(4‐chloroanilino)‐1,3,2λ⁴‐diazaphosphol‐2‐one, C₁₂H₁₁ClN₃OP, each molecule is connected with four neighbouring molecules through (N—H)₂…O hydrogen bonds. These hydrogen bonds form a tubular arrangement along the [001] direction built from R ³₃(12) and R ₄³(14) hydrogen‐bond ring motifs, combined with a C (4) chain motif. The hole constructed in the tubular architecture includes a 12‐atom arrangement (three P, three N, three O and three H atoms) belonging to three adjacent molecules hydrogen bonded to each other. One of the N—H groups of the diazaphosphole ring, not co‐operating in classical hydrogen bonding, takes part in an N—H…π interaction. This interaction occurs within the tubular array and does not change the dimension of the hydrogen‐bond pattern. The energies of the N—H…O and N—H…π hydrogen bonds were studied by NBO (natural bond orbital) analysis, using the experimental hydrogen‐bonded cluster of molecules as the input file for the chemical calculations. In the ¹H NMR experiment, the nitrogen‐bound proton of the diazaphosphole ring has a high value of 17.2 Hz for the ²J _H–P coupling constant.     

On the Correlation between the Blue Shift of Hydrogen Bonding and the Proton Donor-Proton Acceptor Distance

It is demonstrated that in all types of hydrogen bonds (X—H…Y) there is a balance between the long-range attractive orbital interactions and short-range Pauli/nucleus repulsions. When the proton acceptor approaches the proton donor from distance, the hydrogen bonding energy becomes more negative at relatively large distance, goes through a minimum, and then starts to become less negative when the short-range repulsive forces come into effect.Meanwhile, the X--H bond length increases at relatively large distances, goes through a maximum and starts to shorten when the short-range repulsive forces come into effect. Whether the hydrogen bond is red or blue shifted is dictated by the energy minimum position. If at the energy minimum position the X—H bond length is shorter than that for the free monomer, the hydrogen bond is blue shifted and vice versa. Further studies demonstrate that the recent report about the correlation of C—H bond lengths with proton donor-acceptor distance in F3C—H…OH2 and F3C—H…Cl^- is not fully correct because the authors conducted an inappropriate comparison. Furthermore, it is shown for the first time that the Pauli/nucleus repulsion theory is applicable to the blue-shifted hydrogen bonds in the X—H…π complexes and the blue-shifted lithium bonds in the X—Li…Y complexes.     

Nature of MoH···I bonds in Cp2Mo(L)H···I‐C≡C‐R Complexes (L=H,CN, PPh2, C(CH3)3; R=NO2, Cl,Br, H,OH, CH3, NH2)

The nature of the MoH···I bond in Cp₂Mo(L)H···I‐C≡C‐R (L= H, CN, PPh₂, C(CH₃)₃; R=NO₂, Cl, Br, H, OH, CH₃, NH₂) was investigated using electrostatic potential analysis, topological analysis of the electron density, energy decomposition analysis and natural bond orbital analysis. The calculated results show that MoH···I interactions in the title complexes belong to halogen‐hydride bond, which is similar to halogen bonds, not hydrogen bonds. Different to the classical halogen bonds, the directionality of MoH···I bond is low; Although electrostatic interaction is dorminant, the orbital interactions also play important roles in this kind of halogen bond, and steric interactions are weak; the strength of H···I bond can tuned by the most positive electrostatic potential of the I atom. As the electron‐withdrawing ability of the R substituent in the alkyne increases, the electrostatic potential maximum of the I atom increases, which enhances the strength of the H···I halogen bond, as well as the electron transfer.     

Comparative analysis of blue‐shifted hydrogen bond versus conventional hydrogen bond in methyl radical complexes

Methyl radical complexes H₃C…HCN and H₃C…HNC have been investigated at the UMP2(full)/aug‐cc‐pVTZ level to elucidate the nature of hydrogen bonds. To better understand the intermolecular H‐bond interactions, topological analysis of electron density at bond critical points (BCP) is executed using Bader's atoms‐in‐molecules (AIM) theory. Natural bond orbital (NBO) analysis has also been performed to study the orbital interactions and change of hybridization. Theoretical calculations show that there is no essential difference between the blue‐shift H‐bond and the conventional one. In H₃C…HNC complex, rehybridization is responsible for shortening of the N? H bond. The hyperconjugative interaction between the single electron of the methyl radical and N? H antibonding orbital is up to 7.0 kcal/mol, exceeding 3.0 kcal/mol, the upper limit of hyperconjugative n(Y)→σ*(X–H) interaction to form the blue‐shifted H‐bond according to Alabugin's theory. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

The structure, properties, and nature of HArF-HOX (X = F, Cl, Br) complex: an ab initio study and an unusual short hydrogen bond

The structure and properties (geometric, energetic, electronic, spectroscopic, and thermodynamic properties) of HArF‐HOX (X = F, Cl, Br) complex have been investigated at the MP2/aug‐cc‐pVTZ level. Three types of complexes are formed through a hydrogen bond or a halogen bond. The HArF‐HOX complex is the most stable, followed by the FArH‐OHX complex, and the HArF‐XOH complex is the most unstable. The binding distance in FArH‐OHX complex is very short (1.1–1.7 Å) and is smaller than that in HArF‐HOX complex. However, the interaction strength in the former is weaker than that in the latter. Thus, an unusual short hydrogen bond is present in FArH‐OHX complex. The associated H‐Ar bond exhibits a red shift, whereas the distant one gives a blue shift. A similar result is also found for the O? H and O? X bonds. The isotropic chemical shift is negative for the associated hydrogen atom but is positive for the associated halogen atom. However, a reverse result is found for the anisotropic chemical shift. The analyses of natural bond orbital and atoms in molecules have been performed for these complexes to understand the nature and properties of hydrogen and halogen bonds. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Biphenyl: A stress tensor and vector‐based perspective explored within the quantum theory of atoms in molecules

We use quantum theory of atoms in molecules (QTAIM) and the stress tensor topological approaches to explain the effects of the torsion φ of the C‐C bond linking the two phenyl rings of the biphenyl molecule on a bond‐by‐bond basis using both a scalar and vector‐based analysis. Using the total local energy density H( r _b), we show the favorable conditions for the formation of the controversial H–H bonding interactions for a planar biphenyl geometry. This bond‐by‐bond QTAIM analysis is found to be agreement with an earlier alternative QTAIM atom‐by‐atom approach that indicated that the H–H bonding interaction provided a locally stabilizing effect that is overwhelmed by the destabilizing role of the C‐C bond. This leads to a global destabilization of the planar biphenyl conformation compared with the twisted global minimum. In addition, the H( r _b) analysis showed that only the central torsional C‐C bond indicated a minimum for a torsion φ value coinciding with that of the conventional global energy minimum. The H–H bonding interactions are found to be topologically unstable for any torsion of the central C‐C bond away from the planar biphenyl geometry. Conversely, we demonstrate that for 0.0° < φ < 39.95° there is a resultant increase in the topological stability of the C nuclei comprising the central torsional C‐C bond. Evidence is found of the effect of the H–H bonding interactions on the torsion φ of the central C‐C bond of the biphenyl molecule in the form of the QTAIM response β of the total electronic charge density ρ( r _b). Using a vector‐based treatment of QTAIM we confirm the presence of the sharing of chemical character between adjacent bonds. In addition, we present a QTAIM interpretation of hyperconjugation and conjugation effects, the former was quantified as larger in agreement with molecular orbital (MO) theory. The stress tensor and the QTAIM H atomic basin path set areas are independently found to be new tools relevant for the incommensurate gas to solid phase transition occurring in biphenyl for a value of the torsion reaction coordinate φ ≈ 5°. © 2015 Wiley Periodicals, Inc.     

A Combined Experimental and Theoretical Electron Density Study of Intra‐ and Intermolecular Interactions in Thiourea S,S‐Dioxide

The thiourea S,S‐dioxide molecule is recognized as a zwitterion with a high dipole moment and an unusually long C? S bond. The molecule has a most interesting set of intermolecular interactions in the crystalline state—a relatively strong O???H? N hydrogen bond and very weak intermolecular C???S and N???O interactions. The molecule has C_s symmetry, and each oxygen atom is hydrogen‐bonded to two hydrogen atoms with O???H? N distances of 2.837 and 2.826 Å and angles of 176.61 and 158.38°. The electron density distribution is obtained both from Xray diffraction data at 110 K and from a periodic density functional theory (DFT) calculation. Bond characterization is made in terms of the analysis of topological properties. The covalent characters of the C? N, N? H, C? S, and S? O bonds are apparent, and the agreement on the topological properties between experiment and theory is adequate. The features of the Laplacian distributions, bond paths, and atomic domains are comparable. In a systematic approach, DFT calculations are performed based on a monomer, a dimer, a heptamer, and a crystal to see the effect on the electron density distribution due to the intermolecular interactions. The dipole moment of the molecule is enhanced in the solid state. The typical values of ρ_b and H_b of the hydrogen bonds and weak intermolecular C???S and N???O interactions are given. All the interactions are verified by the location of the bond critical point and its associated topological properties. The isovalue surface of Laplacian charge density and the detailed atomic graph around each atomic site reveal the shape of the valence‐shell charge concentration and provide a reasonable interpretation of the bonding of each atom.     

Triel–hydride triel bond between ZX3 (Z = B and Al; X = H and Me) and THMe3 (T = Si,Ge and Sn)

Complexes between THMe₃ (T = Si, Ge and Sn) and ZX₃ (Z = B and Al; X = H and Me) have been characterized using MP2/aug‐cc‐pVTZ calculations. These complexes are chiefly stabilized by a triel–hydride triel bond with the T–H bond pointing to the π‐hole on the triel atom. The triel–hydride interaction is mainly attributed to the charge transfer from the T–H bond orbital to the empty p orbital of the triel atom. These complexes are very stable with a large interaction energy (>10 kcal mol^?1) excluding THMe₃···BMe₃ (T = Si and Ge), indicating that the sp²‐hydridized triel atom has a strong affinity for the T–H bond. The formation of THMe₃···BH₃ results in proton transfer, characterized by conversion of orbital interaction and large charge transfer (ca 0.5e). The large deformation is primarily responsible for the abnormally greater interaction energy in THMe₃···BH₃ (>30 kcal mol^?1) than in the AlH₃ analogue. Methyl substitution on the triel atom weakens the triel–hydride interaction and causes a larger interaction energy in THMe₃···AlMe₃ with respect to its BMe₃ counterpart. Most of these interactions possess characteristics of covalent bonds. Polarization makes a contribution to the stability of most complexes nearly equivalent to the electrostatic term.     

Using (FH)2 and (FH)3 to Bridge the σ‐Hole and the Lone Pair at P in Complexes with H2XP,for X=CH3, OH,H, CCH,F, Cl,NC, and CN

Ab initio MP2/aug′‐cc‐pVTZ calculations are used to investigate the binary complexes H₂XP:HF, the ternary complexes H₂XP:(FH)₂, and the quaternary complexes H₂XP:(FH)₃, for X=CH₃, OH, H, CCH, F, Cl, NC, and CN. Hydrogen‐bonded (HB) binary complexes are formed between all H₂XP molecules and FH, but only H₂FP, H₂ClP, and H₂(NC)P form pnicogen‐bonded (ZB) complexes with FH. Ternary complexes with (FH)₂ are stabilized by F?H???P and F?H???F hydrogen bonds and F???P pnicogen bonds, except for H₂(CH₃)P:(FH)₂ and H₃P:(FH)₂, which do not have pnicogen bonds. All quaternary complexes H₂XP:(FH)₃ are stabilized by both F?H???P and F?H???F hydrogen bonds and P???F pnicogen bonds. Thus, (FH)₂ with two exceptions, and (FH)₃ can bridge the σ‐hole and the lone pair at P in these complexes. The binding energies of H₂XP:(FH)₃ complexes are significantly greater than the binding energies of H₂XP:(FH)₂ complexes, and nonadditivities are synergistic in both series. Charge transfer occurs across all intermolecular bonds from the lone‐pair donor atom to an antibonding σ* orbital of the acceptor molecule, and stabilizes these complexes. Charge‐transfer energies across the pnicogen bond correlate with the intermolecular P?F distance, while charge‐transfer energies across F?H???P and F?H???F hydrogen bonds correlate with the distance between the lone‐pair donor atom and the hydrogen‐bonded H atom. In binary and quaternary complexes, charge transfer energies also correlate with the distance between the electron‐donor atom and the hydrogen‐bonded F atom. EOM‐CCSD spin‐spin coupling constants ^2hJ(F–P) across F?H???P hydrogen bonds, and ^1pJ(P–F) across pnicogen bonds in binary, ternary, and quaternary complexes exhibit strong correlations with the corresponding intermolecular distances. Hydrogen bonds are better transmitters of F–P coupling data than pnicogen bonds, despite the longer F???P distances in F?H???P hydrogen bonds compared to P???F pnicogen bonds. There is a correlation between the two bond coupling constants ^2hJ(F–F) in the quaternary complexes and the corresponding intermolecular distances, but not in the ternary complexes, a reflection of the distorted geometries of the bridging dimers in ternary complexes.     

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     

Halogen Bonding： An AIM Analysis of the Weak Interactions

A series of complexes formed between halogen-containing molecules and ammonia have been investigated by means of the atoms in molecules （AIM） approach to gain a deeper insight into halogen bonding. The existence of the halogen bond critical points （XBCP） and the values of the electron density （Pb） and Laplacian of electron density （V2pb） at the XBCP reveal the closed-shell interactions in these complexes. Integrated atomic properties such as charge, energy, polarization moment, volume of the halogen bond donor atoms, and the corresponding changes （△） upon complexation have been calculated. The present calculations have demonstrated that the halogen bond represents different AIM properties as compared to the well-documented hydrogen bond. Both the electron density and the Laplacian of electron density at the XBCP have been shown to correlate well with the interaction energy, which indicates that the topological parameters at the XBCP can be treated as a good measure of the halogen bond strength In addition, an excellent linear relationship between the interatomic distance d（X…N） and the logarithm of Pb has been established.     

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     

A comparative study of some lithium and hydrogen‐bonded complexes: Ab initio and QTAIM studies

The nature of the interactions of cyanide with lithium and hydrogen halides was investigated using ab initio calculations and topological analysis of electron density. The computed properties of the lithium‐bonded complexes RCN···LiX (R = H, F, Cl, Br, C?CH, CH?CH₂, CH₃, C₂H₅; X = Cl, Br) were compared with those of corresponding hydrogen‐bonded complexes RCN···HX. The results show that both types of intermolecular interactions are “closed‐shell” noncovalent interactions. The effect of substitution on the interaction energy and electron density at the bond critical points of the lithium and hydrogen bonding interactions is similar. In comparison, the interaction energies of lithium‐bonded complexes are more negative than those of hydrogen‐bonded counterparts. The electrostatic interaction plays a more important role in the lithium bond than in the hydrogen bond. On complex formation, the net charge and energy of the Li atom decrease and the atomic volume increases, while the net charge and energy of the H atom increase and the atomic volume decreases. © 2013 Wiley Periodicals, Inc.