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.     

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.     

An ab initio study on substituent and cooperative effects in bifurcated fluorine bonds

Bifurcated fluorine bond (BFB) interactions are studied in model binary complexes pairing N-formyl formamide derivatives and FX molecules (X = F, CN, NC, CF₃ and CCH) by means of ab initio calculations. The calculated F···O binding distances in these complexes are in the range of 2.813–3.048 Å. The corresponding interaction energies lie in a narrow range, from?2.25 to ?16.49 kJ/mol. The nature of BFBs is analysed by a vast number of methods including molecular electrostatic potential, quantum theory of atoms in molecules, non-covalent interaction index and natural bond orbital methods. According to the energy decomposition analysis, the electrostatic and dispersion effects have a dominant role in the formation of these complexes. The formation of a hydrogen- and lithium-bonding interaction tends to increase the strength of BFBs in the ternary XF:NFF-H:NH₃ and XF:NFF-Li:NH₃ complexes, respectively.     

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.     

Substituent effects in cooperativity of chalcogen bonds

In the present work, substituent effects on cooperativity of S···N chalcogen bonds are studied in XHS···NCHS···4-Z–Py (X = F, Cl; Z = H, F, OH, CH₃, NH₂, NO₂, and CN; and Py = pyridine) complexes using ab initio calculations. An increased attraction or a positive cooperativity is observed on introduction of a third molecule to the XHS···NCHS and NCHS···4-Z–Py binary systems. The shortening of each chalcogen bond distance in the ternary systems is dependent on the substituent Z and is increased in the order Z = NH₂ > OH > CH₃ > H > F > CN > NO₂. The electronic aspects of the complexes are analysed using molecular electrostatic potential, and the parameters derived from the atoms in molecules and natural bond orbital methodologies. According to interaction energy decomposition analysis, the electrostatic energies are important in the interaction energy of S···N bonds and may be regarded as being responsible for the stability of these complexes.     

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.     

Enhancement effect of lithium bonding on the strength of pnicogen bonds: XH2P···NCLi···NCY as a working model (X = F,Cl; Y = H,F, Cl,CN)

MP2 calculations with aug-cc-pVDZ basis set were used to analyse intermolecular interactions in XH₂P···NCLi···NCY triads (X = F, Cl; Y = H, F, Cl, CN) which are connected via pnicogen bond and lithium bond. To understand the properties of the systems better, the corresponding dyads are also studied. Molecular geometries and interaction energies of dyads, and triads are investigated at the MP2/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 pnicogen bond, show cooperativity with energy values ranging between ?4.73 and ?8.88 kJ mol^?1. A linear correlation was found between the interaction energies and magnitude of the product of most positive and negative electrostatic potentials. According to energy decomposition analysis, it is revealed that the electrostatic interactions are the major source of the attraction in the title complexes.     

Weakly interacting molecular clusters of CO with H2O,SO2, and NO+

Intermolecular interactions in three dimers, CO···H₂O, CO···SO₂, and CO···NO⁺, were studied at the CCSD(T) level of theory, using a series of the augmented correlation consistent polarised basis sets. Interaction energy and its components as well as vibrational spectra for local minima were computed using both harmonic and anharmonic approximations. While CO···H₂O and CO···SO₂ are weakly bound with the binding energies ?7.4 and ?6.4 kJ/mol, CO···NO⁺ is much more stable with the binding energy of ?32.8 kJ/mol corresponding to ΔG = ?4.7 kJ/mol at 254 K.     

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.     

Hydrogen bond nature in formamide (CYHNH2···XH; YO,S, Se,Te; XF,HO, NH2) complexes at their ground and low‐lying excited states

A theoretical study on the nature of hydrogen bond for formamide and its heavy complexes (CYHNH₂···XH; Y?O, S, Se, Te; X?F, HO, NH₂) was performed on the basis of density functional theory and the quantum chemistry analysis. Except for the CYHNH₂···NH₃ complexes, the substitution of O atom at formamide with less electronegative atoms (S, Se, and Te) is found to weaken the hydrogen bond (H‐bond). This substitution results in cyclic structure of hydrated and ammoniated formamide complexes by the formation of bifunctional H‐bonds (Y···H₄X; X···H₃C). Natural bond orbital analysis indicates that the H‐bond is weakened because of less charge transfer from a lone pair orbital of H‐bond acceptor to antibonding orbital of H‐bond donor. The quantum theory of atoms in molecules analysis reveals that the acyclic structure with single H‐bond stabilizes the complexes more than the cyclic structure formed by bifunctional H‐bonds. Natural energy decomposition analysis (NEDA) and block‐localized wavefunction energy decomposition (BLW‐ED) analyses show that the H‐bond stabilization energies of NEDA and BLW‐ED have good correlation with the dissociation energy of formamide complexes and charge transfer from donor to acceptor atom play an important role in H‐bonding. We have also studied the low‐lying electronic excited states (T₁, T₂, and S₁) for CYHNH₂···H₂O complexes to explore the nature of H‐bond on the basis of electronegativity and found that NEDA also establishes a good correlation with relative electronic energy (with respect to their ground state) and H‐bond strength at their excited states. Copyright © 2014 John Wiley & Sons, Ltd.     

Properties of cationic pnicogen-bonded complexes F4-nHnP+:N-base with H–P···N linear and n = 1–4

Ab initio MP2/aug'-cc-pVTZ calculations have been carried out to investigate the pnicogen-bonded complexes F_4-nH_nP⁺:N-base, for n = 1–4, each with a linear or nearly linear H_ax–P···N alignment. The sp³-hybridised nitrogen bases include NH₃, NClH₂, NFH₂, NCl₂H, NCl₃, NFCl₂, NF₂H, NF₂Cl, and NF_3, and the sp bases are NCNH₂, NCCH₃, NP, NCOH, NCCl, NCH, NCF, NCCN, and N₂. Binding energies increase as the P–N distance decreases, with an exponential curve showing this relationship when complexes with sp³ and sp hybridised bases are treated separately. However, the correlations are not as good as they are for the complexes F_4-nH_nP⁺:N-base for n = 0–3 with F–P···N linear. Different patterns are observed for the change in the binding energies of complexes with a particular base as the number of F atoms in the acid changes. Thus, the particular acid–base pair is a factor in determining the binding energies of these complexes.

Three different charge-transfer interactions stabilise these complexes, namely N_lp→σ*P–H_ax, N_lp→σ*P–F_eq, and N_lp→σ*P–H_eq. Unlike the corresponding complexes with F–P···N linear, N_lp→σ*P–H_ax is not always the dominant charge-transfer interaction, since N_lp→σ*P–F_eq is greater in some complexes. N_lp→σ*P–H_eq makes the smallest contribution to the total charge-transfer energy. The total charge-transfer energies of all complexes increase exponentially as the P–N distance decreases in a manner very similar to that observed for the series of complexes with F–P···N linear.

Equation-of-motion coupled cluster singles and doubles (EOM-CCSD) spin–spin coupling constants ^1pJ(P–N) across the pnicogen bond vary with the P–N distance, but different patterns are observed which depend on the nature of the acid, and for some acids, on the hybridisation of the nitrogen base. ^1pJ(P–N) values for complexes of F₃HP⁺ initially increase as the P–N distance decreases, reach a maximum, and then decrease with decreasing P–N distance as the P···N bond acquires increased covalent character. ^1pJ(P–N) for complexes with H–P···N linear and those with F–P···N linear exhibit similar distance dependencies depending on the number of F atoms in equatorial positions and the hybridisation of the base. Complexation may increase, decrease, or leave the P–H_ax distance unchanged, but ¹J(P–H_ax) always decreases relative to the corresponding isolated ion. Decreasing ¹J(P–H_ax) can be related to decreasing intermolecular P–N distance.     

Tuning of tetrel bonds interactions by substitution and cooperative effects in XH3Si···NCH···HM (X = H,F, Cl,Br; M = Li,Na, BeH and MgH) complexes

In this work, the interplay between the tetrel bond and the dihydrogen bond is investigated in ternary XH₃Si···NCH···HM complexes, where X = H, F, Cl, Br and M = Li, Na, BeH, MgH. The nature of Si···N and H···H interactions is studied by molecular electrostatic potential (MEP), noncovalent interaction and electron localisation function analyses. All binding distances in the ternary complexes are shorter than those of isolated XH₃Si···NCH and NCH···HM systems. That is, the two types of interactions have a cooperative effect on each other. The results of the MEP analysis indicate that the enhancement of the Si···N and H···H bonds can mainly be attributed to the electrostatic interaction. The plot of the reduced density gradient versus sign (λ₂)ρ indicates that the location of the spike associated with each interaction in the ternary systems moves slightly towards the negative (λ₂)ρ values with respect to the binary systems. This confirms that both Si···N and H···H interactions in the ternary complexes are strengthened by the presence of other. Besides, cooperative effects lead to a considerable change in the ¹⁴N nuclear quadrupole coupling constant values of the ternary complexes relative to the XH₃Si···NCH complexes.     

Cooperative effects between tetrel bond and other σ–hole bond interactions: a comparative investigation

Covalently bonded atoms of Groups IV–VII tend to have anisotropic charge distributions, the electronic densities being less on the extensions of the bonds (σ–holes) than in the intervening regions. These σ–holes often give rise to positive electrostatic potentials through which the atom can interact attractively and highly directionally with negative sites. In this work, cooperative effects between tetrel bond and halogen/chalcogen/pnicogen bond interactions are studied in multi-component YH₃M···NCX···NH₃ complexes, where Y = F, CN; M = C, Si and X = Cl, SH and PH₂. These effects are analysed in detail in terms of the structural, energetic, charge-transfer and electron density properties of the complexes. The nature of the σ–hole bonds is unveiled by quantum theory of atoms in molecules and natural bond orbital theory. A favourable cooperativity is found with values that range between ?0.34 and ?1.15 kcal/mol. Many-body decomposition of interaction energies indicate that two-body energy term is the most important source of the attraction, which its contribution accounts for 87%–96% of the total interaction energy.     

Cooperative interaction between π-hole and single-electron σ-hole interactions in O2S···NCX···CH3 and O2Se···NCX···CH3 complexes (X = F,Cl, Br and I)

Ab initio calculations are performed to analyse the cooperative effects between π-hole and single-electron σ-hole interactions in O₂S···NCX···CH₃ and O₂Se···NCX···CH₃ complexes, where X = F, Cl, Br and I. These effects are investigated in terms of geometric and energetic features of the complexes, which are computed by UMP2/aug-cc-pVTZ(-PP) method. Our results indicate that the shortening of the each π-hole bond distance in the complexes is dependent on the strength of the σ-hole interaction. The maximum and minimum energetic cooperativity values correspond to the most and least stable complexes studied in the present work. The cooperativity between both types of interaction is chiefly caused by the electrostatic effects. 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.     

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.     

Ab initio study of synergetic effects of two strong interactions of cation–π interaction and lithium bond in M+ ··· phenyl lithium ··· N (M = Li,Na, K; N = H2O and NH3) complex

Quantum chemical calculations have been performed to investigate the interplay between the cation–π interaction and lithium bonding in the M⁺?···?phenyl lithium?···?OH₂ and M⁺?···?phenyl lithium?···?NH₃ (M?=?Li, Na, K) complexes. The cation–π interaction and lithium bonding in the trimers become stronger relative to the dimers. The interaction energy of cation–π interaction is increased by about 4.4–6.3%, while that of lithium bonding is increased by about 5.2–15.9%. The cooperative energy becomes larger for the stronger cation–π interaction and lithium bond. The F atom and methyl group in the phenyl ring impose a reverse effect on the cation–π interaction and lithium bond. The interaction mechanism in the complexes has been understood with the many-body interaction analysis, electrostatic potentials, and energy decomposition.     

Tuning of pnicogen and chalcogen bonds by an aerogen-bonding interaction: a comparative ab initio study

Using high-level ab initio calculations, the cooperativity effects between an aerogen-bonding and a pnicogen- or chalcogen-bonding interactions are studied in ternary Y···PH₂CN···ZO₃ and Y···SHCN···ZO₃ complexes (Y?=?NH₃, N₂ and Z?=?Ar, Kr, Xe). A detailed analysis of the structures, interaction energies and bonding properties is performed on these systems. For each set of the complexes, a favourable cooperativity is observed between Z···N and P/S···N interactions, especially in complexes involving NH₃ and XeO₃ molecules. It is found that for a given Y or Z, the amount of cooperativity effects in Y···PH₂CN···ZO₃ complexes are important than Y···SHCN···ZO₃ ones. For each ternary complex considered, the effect of a Z···N aerogen bond on a P/S···N bond is more pronounced than that of a P/S···N bond on a Z···N bond. The mechanism of the cooperativity effects in the ternary complexes is studied by electron density difference, quantum theory of atoms in molecules and natural bond orbital analyses. The solvent effects are also studied on the interaction energy and cooperativity of Z···N and P/S···N bonds in the ternary systems.     

Electron density characteristics and charge transfer effect of hydrogen bond O–H···Pt(II): atoms in molecules study and natural bond orbital analysis

In this report, we extended the works of Rizzato et al. [Angew. Chem. Int. Ed. 49, 7440 (2010)] on the nature of O–H···Pt hydrogen bond in trans-[PtCl₂(NH₃)(N–glycine)]·H₂O(1·H₂O) complex, by computational study of O–H···Pt interaction in [NBu₄][Pt(C₆F₅)₃(8-hydroxyquinaldine)], with emphasis on charge transfer effect in this interaction of platinum(II) and hydrogen atom. According to the crystallographic geometry reported by José María Casas et al., [NBu₄][Pt(C₆F₅)₃(8-hydroxyquinaldine)] possesses one O–H···Pt hydrogen bridging interaction, similar to the case in trans-[PtCl₂(NH₃)(N–glycine)]·H₂O(1·H₂O) complex. On the basis of topological criteria of electron density, we characterised this O–H···Pt interaction. Charge transferred between platinum(II) and σ^*_O–H orbital in this complex was calculated by using NBO method. The stabilised energy associated to charge transfer was estimated using a direct proportionality, that is 2–3 eV per electron transferred. Charge transfer effects in O–H···Pt hydrogen bonds were studied for these two complexes. Our results indicate that the interaction of O–H···Pt is closed–shell in nature with significant charge transfer, and that charge transfer effect is not negligible in the interaction of O–H···Pt. The second conclusion is different from the result of Rizzato et al.     

Chalcogen bonds formed through π-holes: SO3 complexes with nitrogen and phosphorus bases

The possibility of formation of chalcogen bonds through π-holes is examined by means of ab initio calculations at MP2/aug-cc-pVTZ level. We investigate the complexes of SO₃ with a series of electron-donating nitrogen bases (ZN), including NH₃, H₂C═NH, NH₂F, NP, NCH, NCF, NF₃, and N₂, and their phosphorous analogues (ZP). The ZN:SO₃ complexes show shorter chalcogen bond distances than the ZP:SO₃ counterparts, accompanied with the more negative interaction energy in the former than in the latter. To understand the nature of interactions, molecular electrostatic potential analysis is performed. In addition, the quantum theory of atoms in molecules and electron localisation function are employed to analyse the chalcogen bond properties in these complexes.     

Chalcogen bonds tuned by an N–H···π or C–H···π interaction: investigation of substituent,cooperativity and solvent effects

In the present work, ab-initio calculations are performed to investigate cooperativity effects between chalcogen bond and H···π interactions in XHY···NCH···C₆H₆ and XHY···CNH···C₆H₆ complexes, where X = F, Cl, Br, CN, NC, and Y = S, Se. The nature of these interactions and the mechanism of cooperativity are studied by means of quantum theory of atoms in molecules, noncovalent interaction index, many-body analysis of interaction energy and electron density shift analysis. For each ternary complex, the shortening of the Y···N(C) distance is more pronounced than that of the H···π. The cooperative energies of these complexes are all negative which demonstrate a positive cooperativity between the Y···N(C) and H···π interactions. The many-body analysis of interaction energy reveals that the two-body energy term has the largest contribution to the total interaction energies of ternary complexes. A good linear correlation is established between the three-body energy and cooperative energy values in the ternary systems. The cooperative energies of XHY···CNH···C₆H₆ complexes indicate a larger sensitivity on the polarity of solvent than XHY···NCH···C₆H₆ ones.