Halogen bond tunability I: the effects of aromatic fluorine substitution on the strengths of halogen-bonding interactions involving chlorine,bromine, and iodine

In the past several years, halogen bonds have been shown to be relevant in crystal engineering and biomedical applications. One of the reasons for the utility of these types of noncovalent interactions in the development of, for example, pharmaceutical ligands is that their strengths and geometric properties are very tunable. That is, substitution of atoms or chemical groups in the vicinity of a halogen can have a very strong effect on the strength of the halogen bond. In this study we investigate halogen-bonding interactions involving aromatically-bound halogens (Cl, Br, and I) and a carbonyl oxygen. The properties of these halogen bonds are modulated by substitution of aromatic hydrogens with fluorines, which are very electronegative. It is found that these types of substitutions have dramatic effects on the strengths of the halogen bonds, leading to interactions that can be up to 100% stronger. Very good correlations are obtained between the interaction energies and the magnitudes of the positive electrostatic potentials (σ-holes) on the halogens. Interestingly, it is seen that the substitution of fluorines in systems containing smaller halogens results in electrostatic potentials resembling those of systems with larger halogens, with correspondingly stronger interaction energies. It is also shown that aromatic fluorine substitutions affect the optimal geometries of the halogen-bonded complexes, often as the result of secondary interactions.     

The cooperativity between hydrogen and halogen bond in the XY···HNC···XY (X,Y = F,Cl, Br) complexes

The cooperativity between hydrogen and halogen bonds in XY···HNC···XY (X, Y = F, Cl, Br) complexes was studied at the MP2/aug-cc-pVTZ level. Two hydrogen-bonded dimers, five hydrogen-bonded dimers, and ten trimers were obtained. The hydrogen- and halogen-bonded interaction energies in the trimers were larger than those in the dimers, indicating that both the hydrogen bonding interaction and the halogen bonding interaction are enhanced. The binary halogen bonding interaction plays the most important role in the ternary system. The hydrogen donor molecule influences the magnitude of the halogen bonding interaction much more than the hydrogen bonding interaction in the trimers with respect to the dimers. Our calculations are consistent with the conclusion that the stronger noncovalent interaction has a bigger effect on the weaker one. The variation in the vibrational frequency in the HNC molecule was considered. The NH antisymmetry vibration frequency has a blue shift, whereas the symmetry vibration frequency has a red shift. A dipole moment enhancement is observed upon formation of the trimers. The variation in topological properties at bond critical points was obtained using the atoms in molecules method, and was consistent with the results of the interaction energy analysis.     

Interplay between halogen bonds and hydrogen bonds in OH/SH···HOX···HY (X = Cl, Br; Y = F, Cl, Br) complexes

The character of the cooperativity between the HOX···OH/SH halogen bond (XB) and the Y―H···(H)OX hydrogen bond (HB) in OH/SH···HOX···HY (X = Cl, Br; Y = F, Cl, Br) complexes has been investigated by means of second-order Møller?Plesset perturbation theory (MP2) calculations and “quantum theory of atoms in molecules” (QTAIM) studies. The geometries of the complexes have been determined from the most negative electrostatic potentials (V _S,min) and the most positive electrostatic potentials (V _S,max) on the electron density contours of the individual species. The greater the V _S,max values of HY, the larger the interaction energies of halogen-bonded HOX···OH/SH in the termolecular complexes, indicating that the ability of cooperative effect of hydrogen bond on halogen bond are determined by V _S,max of HY. The interaction energies, binding distances, infrared vibrational frequencies, and electron densities ρ at the BCPs of the hydrogen bonds and halogen bonds prove that there is positive cooperativity between these bonds. The potentiation of hydrogen bonds on halogen bonds is greater than that of halogen bonds on hydrogen bonds. QTAIM studies have shown that the halogen bonds and hydrogen bonds are closed-shell noncovalent interactions, and both have greater electrostatic character in the termolecular species compared with the bimolecular species.

The unrestricted local properties: application in nanoelectronics and for predicting radicals reactivity

I have performed quantum chemical calculations for the CF₃X (X = Cl, Br) ???Au_n (n?=?2, 3, and 4) complexes at M05-2X/aug-cc-pVDZ(PP) level. Two types of optimized structures were obtained. Type I complexes are stabilized by the coordination force between the negative electrostatic potential of halogen atom and the gold atom, and type II complexes contain halogen bonds formed between the σ-hole of the halogen atoms and the negative electrostatic potential of Au_n. Results of the interaction energy indicate that type I complexes are more stable than type II complexes. AIM analysis reveals that type II complexes are a closed shell interaction and there is a partially covalent nature for type I complexes.     

Theoretical study of noncovalent interactions in XCN···YO<Subscript>2</Subscript>H (X = F,Cl, Br,I; Y = P,As, Sb) complexes

Noncovalent interactions in XCN···YO₂H (X = F, Cl, Br, I; Y = P, As, Sb) complexes were investigated using ab initio calculations at the MP2/aug-cc-pVDZ level of theory. There are four different configurations of these complexes, and the complexes are formed via hydrogen bonds, halogen bonds, π-hole interactions, or dual interactions. An examination of binding distances and interaction energies suggested that π-hole bonds are more stable than the other interactions. Molecular electrostatic potentials, electron densities, second-order stabilization energies, and electron density differences were computed to study the character of these interactions.     

Chalcogen- and halogen-bonds involving SX<Subscript>2</Subscript> (X = F,Cl, and Br) with formaldehyde

The capacity of SX₂ (X = F, Cl, and Br) to engage in different kinds of noncovalent bonds was investigated by ab initio calculations. SCl₂ (SBr₂) has two σ-holes upon extension of Cl (Br)?S bonds, and two σ-holes upon extension of S?Cl (Br) bonds. SF₂ contains only two σ-holes upon extension of the F?S bond. Consequently, SCl₂ and SBr₂ form chalcogen and halogen bonds with the electron donor H₂CO while SF₂ forms only a chalcogen bond, i.e., no F···O halogen bond was found in the SF₂:H₂CO complex. The S···O chalcogen bond between SF₂ and H₂CO is the strongest, while the strongest halogen bond is Br···O between SBr₂ and H₂CO. The nature of these two types of noncovalent interaction was probed by a variety of methods, including molecular electrostatic potentials, QTAIM, energy decomposition, and electron density shift maps. Termolecular complexes X₂S···H₂CO···SX′₂ (X = F, Cl, Br, and X′ = Cl, Br) were constructed to study the interplay between chalcogen bonds and halogen bonds. All these complexes contained S···O and Cl (Br)···O bonds, with longer intermolecular distances, smaller values of electron density, and more positive three-body interaction energies, indicating negative cooperativity between the chalcogen bond and the halogen bond. In addition, for all complexes studied, interactions involving chalcogen bonds were more favorable than those involving halogen bonds.

Open image in new window

Graphical Abstract Molecular electrostatic potential and contour map of the Laplacian of the electron density in Cl₂S···H₂CO···SCl₂ complex

     

Effect of superalkali substituents on the strengths and properties of hydrogen and halogen bonds

Quantum chemical calculations have been performed for the complexes Li₃OCCX–Y (X?=?Cl, Br, H; Y?=?NH₃, H₂O, H₂S) and Li₃OCN–X′Y′ (X′Y′?=?ClF, BrCl, BrF, HF) to study the role of superalkalis in hydrogen and halogen bonds. The results show that the presence of an Li₃O cluster in a Lewis acid weakens its acidity, while its presence in a Lewis base enhances its basicity. Furthermore, the latter effect is more prominent than the former one, and the presence of an Na₃O cluster causes an even greater effect than Li₃O. The strengths of hydrogen and halogen bonds were analyzed using molecular electrostatic potentials. The contributions of superalkalis to the strength of hydrogen and halogen bonds were elucidated by analyzing differences in electron density.     

Theoretical investigations of the H···π and X (X = F,Cl, Br,I)···π complexes between hypohalous acids and benzene

The H···π and X (X = F, Cl, Br, I)···π interactions between hypohalous acids and benzene are investigated at the MP2/6-311++G(2d,2p) level. Four hydrogen-bonded and three halogen-bonded complexes were obtained. Ab initio calculations indicate that the X···π interaction between HOX and C₆H₆ is mainly electrostatically driven, and there is nearly an equal contribution from both electrostatic and dispersive energies in the case of XOH–C₆H₆ complexes. Natural bond orbital (NBO) analysis reveals that there exists charge transfer from benzene to hypohalous acids. Atom in molecules (AIM) analysis locates bond critical points (BCP) linking the hydrogen or halogen atom and carbon atom in benzene.     

A theoretical study on the halogen bonding interactions of C6F5I with a series of group 10 metal monohalides

The halogen bonding interactions between C₆F₅I and a series of transition metal monohalides trans-[M(X)(2-C₅NF₄)-(PR₃)₂] (M = Ni, Pd, Pt; X = F, Cl, Br; R = Me, Cy) have been studied with quantum chemical calculations. Optimized geometries of the halogen bonding complexes indicate that angles C₁-I···X are basically linear (178–180°) and angles I···X-M mainly range from 90 to 150°. The strength of these metal-influenced halogen bonds alters with different metal centers, metal-bound halogen atoms and the substitutes on phosphine ligands. Electrostatic potential and natural bond orbital analysis show that both of the electrostatic and orbital interactions make a contribution to the formation of halogen bonds, while the electrostatic term plays a dominant role. AIM analysis suggests that, for trans-[M(F)(2-C₅NF₄)-(PR₃)₂] (M = Ni, Pd, Pt) monomers, the formed halogen bonding complexes are stabilized by local concentration of the charge of intermediate character, while for the metal monomers containing chlorine and bromine, a typical closed-shell interaction exist. These results prove that the structures and geometries of these halogen bonding complexes can be tuned by changing the halogen atoms and metal centers, which may provide useful information for the design and synthesis of new functional materials.

Theoretical studies on identity SN2 reactions of lithium halide and methyl halide: A microhydration model

Reactions of lithium halide (LiX, X = F, Cl, Br and I) and methyl halide (CH₃X, X = F, Cl, Br and I) have been investigated at the B3LYP/6-31G(d) level of theory using the microhydration model. Beginning with hydrated lithium ion, four or two water molecules have been conveniently introduced to these aqueous-phase halogen-exchange S_N2 reactions. These water molecules coordinated with the center metal lithium ion, and also interacted with entering and leaving halogen anion via hydrogen bond in complexes and transition state, which to some extent compensated hydration of halogen anion. At 298 K the reaction profiles all involve central barriers ΔE _cent which are found to decrease in the order F > Cl > Br > I. The same trend is also found for the overall barriers (ΔE _ovr ) of the title reaction. In the S_N2 reaction of sodium iodide and methyl iodide, the activation energy agrees well with the aqueous conductometric investigation.     

Competition between halogen, dihalogen and hydrogen bonds in bromo- and iodomethanol dimers

O-H…X and O-H…O H-bonds as well as C-X…X dihalogen and C-X…O halogen bonds have been investigated in halomethanol dimers (bromomethanol dimer, iodomethanol dimer, difluorobromomethanol…bromomethanol complex and difluoroiodomethanol…iodomethanol complex). Structures of all complexes were optimized at the counterpoise-corrected MP2/cc-pVTZ level and single-point energies were calculated at the CCSD(T)/aug-cc-pVTZ level. Energy decomposition for the bromomethanol dimer complex was performed using the DFT-SAPT method based on the aug-cc-pVTZ basis set. OH…O and OH…X H-bonds are systematically the strongest in all complexes investigated, with the former being the strongest bond. Halogen and dihalogen bonds, being of comparable strength, are weaker than both H-bonds but are still significant. The strongest bonds were found in the difluoroiodomethanol…iodomethanol complex, where the O-H…O H-bond exceeds 7 kcal mol^-1, and the halogen and dihalogen bonds exceed 2.5 and 2.3 kcal mol^-1, respectively. Electrostatic energy is dominant for H-bonded structures, in halogen bonded structures electrostatic and dispersion energies are comparable, and, finally, for dihalogen structures the dispersion energy is clearly dominant.

How do halogen bonds (S–O⋯I,N–O⋯I and C–O⋯I) and halogen–halogen contacts (C–I⋯I–C,C–F⋯F–C) subsist in crystal structures? A quantum chemical insight

Thirteen X-ray crystal structures containing various non-covalent interactions such as halogen bonds, halogen–halogen contacts and hydrogen bonds (I?N, I?F, I?I, F?F, I?H and F?H) were considered and investigated using the DFT-D3 method (B97D/def2-QZVP). The interaction energies were calculated at MO62X/def2-QZVP and MP2/aug-cc-pvDZ level of theories. The higher interaction and dispersion energies (2nd crystal) of ?9.58 kcal mol^?1 and ?7.10 kcal mol^?1 observed for 1,4-di-iodotetrafluorobenzene bis [bis (2-phenylethyl) sulfoxide] structure indicates the most stable geometrical arrangement in the crystal packing. The electrostatic potential values calculated for all crystal structures have a positive σ-hole, which aids understanding of the nature of σ-hole bonds. The significance of the existence of halogen bonds in crystal packing environments was authenticated by replacing iodine atoms by bromine and chlorine atoms. Nucleus independent chemical shift analysis reported on the resonance contribution to the interaction energies of halogen bonds and halogen–halogen contacts. Hirshfeld surface analysis and topological analysis (atoms in molecules) were carried out to analyze the occurrence and strength of all non-covalent interactions. These analyses revealed that halogen bond interactions were more dominant than hydrogen bonding interactions in these crystal structures.

Open image in new window

Graphical Abstract Molecluar structure of 1,4-Di-iodotetrafluorobenzene bis(thianthrene 5-oxide) moelcule and its corresponding molecular electrostatic potential map for the view of σ-hole.

     

Characteristics of beryllium bonds; a QTAIM study

The nature of beryllium bonds formed between BeX2 (X is H, F and Cl) and some Lewis bases have been investigated. The distribution of the Laplacian of electron density shows that there is a region of charge depletion around the Be atom, which, according to Laplacian complementary principal, can interact with a region of charge concentration of an atom in the base and form a beryllium bond. The molecular graphs of the investigated complexes indicate that beryllium in BeH2 and BeF2 can form “beryllium bonds” with O, N and P atoms but not with halogens. In addition, eight criteria based on QTAIM properties, including the values of electron density and its Laplacian at the BCP, penetration of beryllium and acceptor atom, charge, energy, volume and first atomic moment of beryllium atom, have been considered and compared with the corresponding ones in conventional hydrogen bonds. These bonds share many common features with very strong hydrogen bonds, however,some differences have also been observed.     

Parametrization of halogen bonds in the CHARMM general force field: Improved treatment of ligand–protein interactions

A halogen bond is a highly directional, non-covalent interaction between a halogen atom and another electronegative atom. It arises due to the formation of a small region of positive electrostatic potential opposite the covalent bond to the halogen, called the ‘sigma hole.’ Empirical force fields in which the electrostatic interactions are represented by atom-centered point charges cannot capture this effect because halogen atoms usually carry a negative charge and therefore interact unfavorably with other electronegative atoms. A strategy to overcome this problem is to attach a positively charged virtual particle to the halogen. In this work, we extend the additive CHARMM General Force Field (CGenFF) to include such interactions in model systems of phenyl-X, with X being Cl, Br or I including di- and trihalogenated species. The charges, Lennard-Jones parameters, and halogen-virtual particle distances were optimized to reproduce the orientation dependence of quantum mechanical interaction energies with water, acetone, and N-methylacetamide as well as experimental pure liquid properties and relative hydration free energies with respect to benzene. The resulting parameters were validated in molecular dynamics simulations on small-molecule crystals and on solvated protein–ligand complexes containing halogenated compounds. The inclusion of positive virtual sites leads to better agreement across experimental observables, including preservation of ligand binding poses as a direct result of the improved representation of halogen bonding.     

Synthesis and characterization of the 1:1 adducts of copper(I) halides with bidentate N,N′-bis(benzophenone)-1,2-diiminoethane Schiff base: Crystal structures of [Cu(bz2en)2][CuX2] (X = Br, I) complexes

1:1 adducts of N,N′-bis(benzophenone)-1,2-diiminoethane (bz₂en) with copper(I) chloride, bromide and iodide, [Cu(bz₂en)₂][CuX₂] (X = Cl, Br, and I), have been synthesized and the structures of the solid bromide and iodide adducts were determined by X-ray crystallography from single-crystal data. The solid-state structure reveals ionic complexes containing a cation of copper(I) ion coordinated to four nitrogen atoms of two bz₂en molecules (distorted tetrahedron) and a linear dibromocuprate(I) and a di-μ-iodo-diiododicuprate(I) anion for the bromo and iodo adducts, respectively. The bromo adduct structure contains CH?Br intermolecular hydrogen bonds. The complexes are very stable towards atmospheric oxygen in the solid state. The spectral properties of the above complexes are also discussed.     

Theoretical study on cooperative effects between X⋯N and X⋯Carbene halogen bonds (X = F,Cl,Br and I)

Quantum chemical calculations are performed to study the interplay between halogen?nitrogen and halogen?carbene interactions in NCX?NCX?CH₂ complexes, where X?=?F, Cl, Br and I. Molecular geometries and interaction energies of dyads and triads are investigated at the MP2/aug-cc-pVTZ level of theory. It is found that the X?N and X?C_carbene interaction energies in the triads are larger than those in the dyads, indicating that both the halogen bonding interactions are enhanced. The estimated values of cooperative energy E _coop are all negative with much larger E _coop in absolute value for the systems including iodine. The nature of halogen bond interactions of the complexes is analyzed using parameters derived from the quantum theory atoms in molecules methodology and energy decomposition analysis.

Theoretical study of the triangular bonding complex formed by carbon tetrabromide, a halide, and a solvent molecule in the gas phase

MP2(full)/aug-cc-pVDZ(-PP) computations predict that new triangular bonding complexes (where X^? is a halide and H–C refers to a protic solvent molecule) consist of one halogen bond and two hydrogen bonds in the gas phase. Carbon tetrabromide acts as the donor in the halogen bond, while it acts as an acceptor in the hydrogen bond. The halide (which commonly acts as an acceptor) can interact with both carbon tetrabromide and solvent molecule (CH₃CN, CH₂Cl₂, CHCl₃) to form a halogen bond and a hydrogen bond, respectively. The strength of the halogen bond obeys the order CBr₄???Cl^? > CBr₄???Br^? > CBr₄???I^?. For the hydrogen bonds formed between various halides and the same solvent molecule, the strength of the hydrogen bond obeys the order C-H???Cl^? > C-H???Br^? > C-H???I^?. For the hydrogen bonds formed between the same halide and various solvent molecules, the interaction strength is proportional to the acidity of the hydrogen in the solvent molecule. The diminutive effect is present between the hydrogen bonds and the halogen bond in chlorine and bromine triangular bonding complexes. Complexes containing iodide ion show weak cooperative effects.

Investigations into the nature of halogen- and hydrogen-bonding interactions of some heteroaromatic rings with dichlorine monoxide

We have studied the structures, properties, and nature of halogen- and hydrogen-bonding interactions between some heteroaromatic rings (C₅H₅N, C₄H₄O, and C₄H₄S) with Cl₂O at the MP2/aug-cc-pVTZ level. We also considered the solvent effect on the halogen bonds and hydrogen bonds in the C₅H₅N-Cl₂O complexes and found that the solvent has a weakening effect on the π-type halogen bond and hydrogen bond but a prominent enhancing effect on σ-type halogen bond. The complexes have also been analyzed with symmetry adapted perturbation theory method (SAPT).     

Computational investigation of a new ion-pair receptor for calix[4]pyrrole

Theoretical studies of a new ion-pair receptor, meso-octamethylcalix[4]pyrrole (OMCP), and its interactions with the halide anions F(-), Cl(-), and Br(-) and the cesium halides CsF, CsCl, and CsBr have been performed. Geometries, binding energies, and binding enthalpies were evaluated with the restricted hybrid Becke three-parameter exchange functional (B3LYP) method using the 6-31+G(d) basis set and relativistic effective core potentials. The optimized geometric structures were used to perform natural bond orbital (NBO) analysis. The two typical types of hydrogen bonds, N-H…X(-) and C-H…X(-), were investigated. The results indicate that hydrogen bonding interactions are dominant, and that the halide anions (F(-), Cl(-), and Br(-)) offer lone pair electrons to the σ*(N-H) or σ*(C-H) antibonding orbitals of OMCP. In addition, electrostatic interactions between the lone pair electrons of the halide anion and the LP* orbitals of Cs(+) as well as cation-π interactions between the metal ion and π-orbitals of the pyrrole rings have important roles to play in the Cs(+)?OMCP?X(-) complexes. The current study further demonstrates that this easy-to-make OMCP host compound functions as not only an anion receptor but also an ion-pair receptor.     

Interaction of glucosamine with uracil and thymine: a computational study

Noncovalent interactions are ubiquitous and have been well recognized in chemistry, biology and material science. Yet, there are still recurring controversies over their natures, due to the wide range of noncovalent interaction terms. In this Essay, we employed the Valence Bond (VB) methods to address two types of interactions which recently have drawn intensive attention, i.e., the halogen bonding and the CH???HC dihydrogen bonding. The VB methods have the advantage of interpreting molecular structures and properties in the term of electron-localized Lewis (resonance) states (structures), which thereby shed specific light on the alteration of the bonding patterns. Due to the electron localization nature of Lewis states, it is possible to define individually and measure both polarization and charge transfer effects which have different physical origins. We demonstrated that both the ab initio VB method and the block-localized wavefunction (BLW) method can provide consistent pictures for halogen bonding systems, where strong Lewis bases NH₃, H₂O and NMe₃ partake as the halogen bond acceptors, and the halogen bond donors include dihalogen molecules and XNO₂ (X?=?Cl, Br, I). Based on the structural, spectral, and energetic changes, we confirm the remarkable roles of charge transfer in these halogen bonding complexes. Although the weak C-H???H-C interactions in alkane dimers and graphane sheets are thought to involve dispersion only, we show that this term embeds delicate yet important charge transfer, bond reorganization and polarization interactions.