Theoretical analysis of trends in hydrogen bonding involving halogen acceptors (F−At) covalently bonded to a group 14 atom (C−Pb)

In this work, extensive quantum-chemical calculations have been carried out to identify and elucidate trends in the hydrogen-bonding (HB) interaction involving halogen acceptors covalently bonded to a group 14 atom. A series of 25 heterodimers composed of MH₃X (where M = C?Pb and X = F?At) and HNC molecules have been selected as model complexes stabilised by the HB interaction occurring between the X atom of MH₃X and the H atom of HNC. The interaction energy (E_int) between MH₃X and HNC in the MH₃X···HNC complexes falls in the range from ?2.7 to ?10.8 kcal/mol, indicating weak or medium strength of HB in these complexes. The strength of HB in the complexes remains consistent with the well-known HB weakening as X gets heavier. Regarding the effect of M on E_int, the gradual strengthening of HB is observed while descending group 14, but only from M = Si to M = Pb. The trends in E_int are compared with various HB-related parameters obtained from vibrational analysis, the natural bond orbital (NBO) method, the symmetry-adapted perturbation theory (SAPT) and the quantum theory of atoms in molecules (QTAIM). The parameters that present clear (possibly linear) relationships with E_int have been selected to characterise the effect of M and X on the HB interaction.     

Strengthening of the halogen-bonding by an aerogen bond interaction: substitution and cooperative effects in O3Z···NCX···NCY (Z = Ar,Kr, Xe; X = Cl,Br, I; Y = H,F, OH) complexes

The geometry, interaction energy and bonding properties of ternary complexes O₃Z···NCX···NCY (Z= Ar, Kr, Xe; X = Cl, Br, I and Y = H, F, OH) are investigated with ab initio calculations at the MP2/aug-cc-pVTZ level. Two different types of intermolecular interactions are present in these complexes, namely, aerogen bond (Z···N) and halogen bond (X···N). The formation mechanism and bonding properties of these complexes are analysed with molecular electrostatic potentials, quantum theory of atoms in molecules and non-covalent interaction index. It is found that the cooperativity energies in the ternary complexes are all negative; that is, the interaction energy of the ternary complex is greater (more negative) than the sum of the interaction energies of the corresponding binary systems. Also, the cooperativity energies increase with the increase of the interaction energies. The cooperative effects in the ternary complexes make a decrease in the total spin–spin coupling constants across the aerogen bonding, J(Z–N), which can be regarded as a proof for the reinforce of Z···N interactions in the ternary complexes with respect to the binary systems.     

Theoretical studies on the nature and strength of an intermolecular non-covalent Te•••π interaction

ABSTRACT

DFT and MP2 calculations were used to determine the nature of non-covalent tellurium–π interactions in R₂Te???C₆H₆ (R = H, F, Cl, CH₃) and C₄H₄Te???C₆H₆ systems. The results showed that the strength of Te···π interaction follows the order F₂Te > Cl₂Te > tellurophene > H₂Te > Me₂Te. Also, the F₂X···π system complexes (X = Te or Se, π system = C₆F₆, C₆(CH₃)₆, Cr(C₆H₆)₂ and coronene) were studied for investigating the direction of charge flow in Te···π interaction. The obtained data expressed that the existence of electron withdrawing group on Te atom increases the strength of Te???π interaction while the fluorine atoms on benzene ring decrease it. The breakdown of ΔE_int in the R₂Te···C₆H₆ (R = H, F, Cl, CH₃) and C₄H₄Te···C₆H₆ systems using two dispersion corrected DFT methods showed that when the R group on divalent tellurium atom is an electron withdrawing substituent, the contribution of ΔE_elstat and ΔE_orb in total interaction energy increases and the value of ΔE_int is relatively large. The present data also showed that the intermolecular Te···π interactions are slightly stronger than corresponding Se···π interactions.     

Interplay and competition between the lithium bonding and halogen bonding: R3C···XCN···LiCN and R3C···LiCN···XCN as a working model (R = H,CH3; X = Cl,Br)

UMP2 calculations with aug-cc-pVDZ basis set were used to analyse intermolecular interactions in R₃C···XCN···LiCN and R₃C···LiCN···XCN triads (R = H, CH₃; X = Cl, Br) which are connected via lithium bond and halogen bond. To understand the properties of the systems better, the corresponding dyads are also studied. Molecular geometries and binding energies of dyads, and triads are investigated at the UMP2/aug-cc-pVDZ computational level. Particular attention is paid to parameters such as cooperative energies, and many-body interaction energies. All studied complexes, with the simultaneous presence of a lithium bond and a halogen bond, show cooperativity with energy values ranging between ?1.20 and ?7.71 kJ mol^?1. A linear correlation was found between the interaction energies and magnitude of the product of most positive and negative electrostatic potentials (V_S,maxV_S,min). The electronic properties of the complexes are analysed using parameters derived from the atoms in molecules (AIM) methodology. According to energy decomposition analysis, it is revealed that the electrostatic interactions are the major source of the attraction in the title complexes.     

Long distance structural consequences of H‐bonding: the case of complexes of para‐substituted phenol derivatives

H‐bonded complexes of p‐X‐PhOH/p‐X‐PhO^? with fluoride and hydrofluoric acid (X = OH, H, NO₂) were subject of optimization (by means of B3LYP/6‐311+G**) for gradually changed O···F distance from d_O···F = 4.0 Å down to (i) the distance of the proton transfer from the hydroxyl group to fluoride leading to O^?···HF interaction and (ii) fully optimized system (O^?···HF type). In this way, we simulate gradual changes of H‐bond strength estimating simultaneously the energy of interaction, E_int, energy of deformation, E_def, and the binding energy, E_tot. The obtained geometrical parameters allow us to show that H‐bond formation causes substantial changes in geometry, even at so distant parts of the system as the ring and bond length in para‐substituents (OH and NO₂). All these changes are monotonically dependent on interaction and deformation energies. Copyright © 2009 John Wiley & Sons, Ltd.     

The strengthening effect of a hydrogen or lithium bond on the Z···N aerogen bond (Z = Ar,Kr and Xe): a comparative study

An ab initio study is performed on O₃Z···NCM···NCX (Z = Ar, Kr and Xe; M = H and Li; X = H, F and CH₃) complexes to investigate cooperativity effects between aerogen and hydrogen or lithium bonding interactions in these systems. To understand the cooperative effects, a detailed analysis of the binding distances, interaction energies and bonding properties is performed on these complexes. The results indicate that all Z···N and H/Li···N binding distances in the ternary complexes are shorter than those of corresponding binary systems. For a given M or X, cooperative energies increase as Z = Xe > Kr > Ar. Moreover, O₃Z···NCLi···NCX complexes exhibit a larger cooperative energy than O₃Z···NCH···NCX ones. The non-covalent interaction (NCI) index analysis indicates that the formation of an H/Li···N interaction in the ternary complexes shifts the location of the spike associated with the Z···N interaction towards the negative λ₂ρ values. This indicates that NCI analysis can be regarded as a useful tool for the study of cooperative effects between two different non-covalent interactions. Also, cooperative effects in O₃Z···NCM···NCX complexes make a decrease in ¹⁴N nuclear quadrupole coupling constants of NCH or NCLi molecule.     

Halogen bond between hypervalent halogens YF3/YF5 (Y=Cl,Br, I) and H2X (X= O,S, Se)

ABSTRACT

The complexes of H₂X (X?=?O, S, Se) with hypervalent halogens YF₃ and YF₅ (Y?=?Cl, Br, I) have been studied. The σ-hole on the Y atom participates in a halogen bond with the lone pair on the chalcogen atom. In addition, some secondary interactions coexist with the halogen bond in most complexes. The interaction energy correlates with the nature of both X and Y atoms. In most cases, the complex is more stable for the heavier Y atom and the lighter X atom. Of course, there are some exceptions in H₂X···YF₃. YF₃ forms a more stable complex with H₂X than does YF₅. These complexes are dominated by electrostatic interaction and the halogen bond involving H₂S and H₂Se exhibits some covalent character.

Halogen bond plays an important role in chemical reactions and multivalent halogens can regulate chemical reactions by participating in a halogen bond. Thus we compare the effect of the chalcogen electron donor on the strength and nature of halogen bonding involving multivalent halogens.     

Cooperativity of tetrel bonds tuned by substituent effects

Substituent effects on cooperativity between N···X and C···X interactions are studied in the 4-Z-Py···XH₃(NC)···XH₃(NC) complexes, where X = C, Si; Z = H, F, OH, CH₃, NH₂, F, NC, CN, NO₂ and Py = pyridine. All N···X and C···X binding distances in the ternary complexes are always shorter than those in the corresponding binary complex. This indicates that the formation of the N···X interaction strengthens C···X bond in these complexes and vice versa. Our results reveal that the strength of N···X and C···X interactions in the ternary complexes considerably depends on the nature of X and Z substituents. For a given aromatic system, the shortening of N···X and C···X distances is more important for SiH₃(NC) complexes than CH₃(NC) counterparts. The mechanism of cooperative effects in the ternary complexes is unveiled by electrostatic potential analyses and natural bond theory.     

Cooperativity between the hydrogen bonding and halogen bonding in F3CX ··· NCH(CNH) ··· NCH(CNH) complexes (X=Cl,Br)

MP2 calculations with the cc-pVTZ basis set were used to analyse the intermolecular interactions in F₃CX?···?NCH(CNH)?···?NCH(CNH) triads (X=Cl, Br), which are connected via hydrogen and halogen bonds. Molecular geometries, binding energies, and infrared spectra of the dyads and triads were investigated at the MP2/cc-pVTZ computational level. Particular attention was given to parameters such as the cooperative energies, cooperative dipole moments, and many-body interaction energies. All studied complexes, with the simultaneous presence of a halogen bond and a hydrogen bond, show cooperativity with energy values ranging between ?1.32 and ?2.88?kJ?mol^?1. The electronic properties of the complexes were analysed using the Molecular Electrostatic Potential (MEP), electron density shift maps and the parameters derived from the Atoms in Molecules (AIM) methodology.     

Theoretical study of cooperativity between hydrogen bond‒hydrogen bond,halogen bond‒halogen bond and hydrogen bond‒halogen bond in ternary FX…diazine…XF (X = H and Cl) complexes

ABSTRACT

In the present work, the cooperativity between hydrogen bond?hydrogen bond, halogen bond?halogen bond and hydrogen bond?halogen bond in ternary FX…diazine…XF (X = H and Cl) complexes is theoretically investigated. The sign of cooperative energy (E_coop) obtained in all of the triads is positive which indicates that the ternary complex is less stable than the sum of the two isolated binary complexes. Moreover, our calculations show that E_coop value in triads increases as FX…pyridazine…XF > FX…pyrimidine…XF > FX…pyrazine…XF. In agreement with energetic, geometrical and topological properties, electrostatic potentials and coupling constants across ¹⁵N…X?¹⁹F (X = ¹H or ³⁵Cl) hydrogen and halogen bonds indicate that hydrogen and halogen bonds are weakened in the considered complexes where two hydrogen and halogen bonds coexist. As compared to N…H hydrogen bond, it is also observed that cooperativity has greater effect on N…Cl halogen bond.     

Tuning of carbon bonds by substituent effects: an ab initio study

ABSTRACT

An ab initio study, at the MP2/aug-cc-pVTZ level of theory, is performed to study σ-hole bond in binary XH₃C···CNY complexes, where X = CN, F, NO₂, CCH and Y = H, OH, NH₂, CH₃, C₂H₅, Li. This type of interaction is labelled as ‘carbon bond’, since a covalently bonded carbon atom acts as the Lewis acid in these systems. The geometrical and energetic parameters of the resulting complexes are analysed in details. The interaction energies of these complexes are between ?4.97 kJ/mol in (HCC)H₃C···CNH and ?23.07 kJ/mol in (O₂N)H₃C···CNLi. It is found that the electrostatic interaction plays a key role in the overall stabilisation of these carbon-bonded complexes. To deepen the understanding of the nature of the carbon-bonding, the molecular electrostatic potential, natural bond orbital, quantum theory of atoms in molecules and non-covalent interaction index analyses are also used. Our results indicate that the carbon bond is favoured over the C-H···C hydrogen bond in the all complexes considered and may suggest the possible important roles of the C···C interactions in the crystal growth and design.     

A theoretical study of substitution effects on halogen–π interactions

Ab initio calculations are performed to analyse the existence of intermolecular halogen···π interactions in NCX complexes with YC≡CY, where X = Cl, Br and Y = H, CN, F, Cl, OH, NH₂, and CH₃. Molecular geometries and interaction energies of the complexes are investigated at the MP2/aug-cc-pVTZ level of theory. Our results indicate that the interaction energies for the NCX···YC≡CY complexes lie in the range between ?0.5 and ?5.9 kcal/mol. The physical nature of the interactions is studied using symmetry-adapted perturbation theory (SAPT). The stability of the X···π interactions is predicted to be attributable mainly to electrostatic and dispersion effects.     

Multi-mode ground state interaction terms in C60-based electron donor-acceptor complexes

Quantum-chemical testing of donor-acceptor properties of binary molecular complexes, related to the singlet state, is suggested as QCh calculations of studied systems and their constituents by using both spin-nondependent (RHF) and spindependent (UHF) versions of the exploited computational tool. The avoided crossing of intermolecular interaction terms of neutral moleculesE _int (A ⁰ B ⁰) and molecular ionsE _int (A ⁺ B ⁻) causes a multi-mode character of the ground state term. The dependence of D-A complex properties on the type of the term, space positions of the term minimum, and the interrelation of the corresponding energies are discussed. The suggested approach has been applied to binary complexes C₆₀+X (X=TAE, TDAE, DMMA, COANP, 2Li, Mg).     

Prediction and characterisation of a chalcogen···π interaction with acetylene as a potential electron donor in XHS···HCCH and XHSe···HCCH (X = F,Cl, Br,CN, OH,OCH3, NH2, CH3) σ-hole complexes

It is well-known that many covalently bonded atoms of group VI have specific positive regions of electrostatic potential (σ-holes) through which they can interact with Lewis bases. This interaction is called ‘chalcogen bond’ by analogy with halogen bond and hydrogen bond. In this study, ab initio calculations are performed to predict and characterise chalcogen···π interactions in XHS···HCCH and XHSe···HCCH complexes, where X = F, Cl, Br, CN, OH, OCH₃, NH₂, CH₃. For the complexes studied here, XHS(Se) and HCCH are treated as a Lewis acid and a Lewis base, respectively. The CCSD(T)/aug-cc-pVTZ interaction energies of this type of σ-hole bonding range from ?1.18 to ?4.83 kcal/mol. The calculated interaction energies tend to increase in magnitude with increasing positive electrostatic potential on the extension of X–S(Se) bond. The stability of chalcogen···π complexes is attributed mainly to electrostatic and correlation effects. The nature of chalcogen···π interactions is unveiled by means of the atoms in molecules, natural bond orbital, and electron localisation function analyses.     

Blue and red shifts in Rg···HCN and Rg···HNC complexes (Rg=He,Ne, Ar,Kr)

The shift in the harmonic vibrational frequency of the H–C stretch of HCN on formation of the linear Rg···HCN complexes, and of the H–N stretch of HNC on the formation of Rg···HNC complexes (Rg?=?He, Ne, Ar, Kr), has been determined by ab initio computations. These shifts are in agreement with predictions from a model based on perturbation theory and involving the first and second derivatives of the interaction energy with respect to displacement of the H–C (H–N) bond length from its equilibrium value in the monomer. Small blue shifts were obtained for He···HCN, Ne···HCN and He···HNC, while red shifts were obtained for the other weakly bound complexes. These vibrational characteristics are rationalized by considering the balance between the interaction energy derivatives obtained from the perturbative model. For all complexes, the IR intensity of the H–C or H–N stretch was increased from the isolated monomer values on complexation.     

Geometry‐based analysis of azulene and azulene‐like systems with H‐ or Li‐bonding

B3LYP/6–311+G** optimization was carried out for azulene and its analogs, in which CH? CH? CH fragment was replaced with O···X···O (X = H or Li). π‐electron delocalization in four possible derivatives with H‐bonding and three possible derivatives with Li‐bonding was described by the use of HOMA index. All derivatives with Li‐bonding exhibit high π‐electron delocalization similar to that found for azulene. Among four H‐bonded systems, two exhibit lower π‐electron delocalization (HOMA < 0.39) and higher total electron energy than the other two derivatives. Copyright © 2007 John Wiley & Sons, Ltd.     

Cooperativity effects between σ-hole interactions: a theoretical evidence for mutual influence between chalcogen bond and halogen bond interactions in F2S···NCX···NCY complexes (X = F,Cl, Br,I; Y = H,F, OH)

Quantum chemical calculations are performed to study the cooperativity effects between chalcogen bond and halogen bond interactions in F₂S···NCX···NCY complexes, where X = F, Cl, Br, I and Y = H, F, OH. These effects are investigated in terms of geometric and energetic features of the complexes, which are computed by second-order Møller–Plesset perturbation theory (MP2). For each F₂S···NCX···NCY complex studied, the effect of cooperativity on the chalcogen bond is dependent on the strength of halogen bond. The results indicate that the interaction energies of chalcogen and halogen bonds in the triads are more negative relative to the respective dyads. The interaction energy of chalcogen bond is increased by 31%–49%, whereas that of halogen bond by 28%–62%. The energy decomposition analysis reveals that electrostatic force plays a main role in the cooperativity effects between the chalcogen bond and halogen bond interactions. The topological analysis, based on the quantum theory of atoms in molecules, is used to characterise the interactions and analyse their enhancement with varying electron density at bond critical points.     

Exploring hydride-π interactions and their tuning by σ-hole bonds: an ab initio study

In the present work, ab initio calculations are performed to investigate the geometry, interaction energy and bonding properties of binary complexes formed between metal-hydrides HMX (M = Be, Mg, Zn and X = H, F, CH₃) and a series of π-acidic heteroaromatic rings. In all the resulting complexes, the heteroaromatic ring acts as a Lewis acid (electron acceptor), while the H atom of the HMX molecule acts as a Lewis base (electron donor). The nature of this interaction, called ‘hydride-π’ interaction, is explored in terms of molecular electrostatic potential, non-covalent interaction, quantum theory of atoms in molecules and natural bond orbital analyses. The results show that the interaction energies of these hydride-π interactions are between ?1.24 and ?2.72 kcal/mol. Furthermore, mutual influence between the hydride-π and halogen- or pnicogen-bonding interactions is studied in complexes in which these interactions coexist. For a given π-acidic ring, the formation of the pnicogen-bonding induces a larger enhancing effect on the strength of hydride-π bond than the halogen-bonding.     

2,2,6,6-四甲基哌啶氮氧自由基与卤仿形成卤键或氢键的理论研究

利用理论计算化学研究了2,2,6,6-四甲基哌啶-N-氧自由基与卤仿形成卤键和氢键络合物的可能性. 从分子静电势、络合物分子的结构参数、络合物的作用能以及自然键轨道理论的角度着手研究. 结果表明,卤键与氢键络合物的键合能均遵循氯化物<溴化物<碘化物,氢键络合物作用强度大于相应的卤键络合物. 因此,卤仿与2,2,6,6-四甲基哌啶-N-氧自由基之间作用模式氢键为主. 需要注意的是,碘仿形成卤键的作用强度与氢键相当,因此在碘仿中,卤键与氢键两种模式应该竞争性的存在.     

Abnormal synergistic effects between Lewis acid–base interaction and halogen bond in F3B···NCX···NCM

An abnormal synergistic effect was found between the Lewis acid–base interaction and halogen bond in triads F₃B···NCX···NCM (X and M are halogen atoms), where the strong Lewis acid–base interaction between F₃B and NCX has a larger enhancement than the weak halogen bond between NCX and NCM. This is in contrast with the traditional cooperative effect. It is interesting that the alkali-metal substituent as well as the heavier halogen atom play a more remarkable role in the enhancement of the interaction F₃B···NCX than that of NCX···NCM, particularly, the alkali-metal substituent makes the abnormal synergistic effect be the traditional cooperative one.