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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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 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.     

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.     

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.     

The triel bond: a potential force for tuning anion–π interactions

Using ab-initio calculations, the mutual influence between anion–π and B···N or B···C triel bond interactions is investigated in some model complexes. The properties of these complexes are studied by molecular electrostatic potential, noncovalent interaction index, quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analyses. According to the results, the formation of B···N or B···C triel bond interactions in the multi-component systems makes a significant shortening of anion–π distance. Such remarkable variation in the anion–π distances has not been reported previously. The strengthening of the anion–π bonding in the multi-component systems depend significantly on the nature of the anion, and it becomes larger in the order Br^? > Cl^? > F^?. The parameters derived from the QTAIM and NBO methodologies are used to study the mechanism of the cooperativity between the anion–π and triel bond interactions in the multi-component complexes.     

Effect of hydrogen bond formation on the NMR properties of microhydrated ortho-aminobenzoic acid

The in?uence of the hydrogen bond formation on the nuclear magnetic resonance parameters has been investigated in the case of microhydrated ortho-aminobenzoic acid (o-Abz) in the gas-phase. DFT-B3LYP/aug-cc-pVDZ predicted ¹H and ¹³C isotropic chemical shifts with respect to TMS of the isolated o-Abz are in reasonable agreement with available experimental data. The isotropic and anisotropic chemical shifts for all atoms of o-Abz within the o-Abz?···?(H₂O)_1-3 complexes have been calculated at the Hartree–Fock, and density functional (B3LYP) theoretical levels using the 6-31++G(2d,2p) and aug-cc-pVDZ basis sets and considering the counterpoise corrections for the basis set superposition errors. The chemical shift values of the carboxyl group atoms of microhydrated o-Abz relative to isolated o-abz do not show significant basis set dependence. Both the hydrogen and carbon atoms constituting the carboxyl group of o-Abz suffer downfield shift due to formation of hydrogen bond with water. The length of hydrogen bond formed between o-Abz and water is found to vary with the number of water molecules present around o-Abz. A direct correlation between the hydrogen bond length and isotropic chemical shift of the bridging hydrogen is observed for both C?=?O?···?H-O and O-H?···?O interactions.