Strong σ‐Hole Activation on Icosahedral Carborane Derivatives for a Directional Halide Recognition

Crystal engineering based on σ‐hole interactions is an emerging approach for realization of new materials with higher complexity. Neutral inorganic clusters derived from 1,2‐dicarba‐closo‐dodecaborane, substituted with ‐SeMe, ‐TeMe, and ‐I moieties on both skeletal carbon vertices are experimentally demonstrated herein as outstanding chalcogen‐ and halogen‐bond donors. In particular, these new molecules strongly interact with halide anions in the solid‐state. The halide ions are coordinated by one or two donor groups (μ₁‐ and μ₂‐coordinations), to stabilize a discrete monomer or dimer motifs to 1D supramolecular zig‐zag chains. Crucially, the observed chalcogen bond and halogen bond interactions feature remarkably short distances and high directionality. Electrostatic potential calculations further demonstrate the efficiency of the carborane derivatives, with V_s,max being similar or even superior to that of reference organic halogen‐bond donors, such as iodopentafluorobenzene.     

Activating Chalcogen Bonding (ChB) in Alkylseleno/Alkyltelluroacetylenes toward Chalcogen Bonding Directionality Control

Activation of a deep electron-deficient area on chalcogen atoms (Ch=Se, Te) is demonstrated in alkynyl chalcogen derivatives, in the prolongation of the (C≡)C−Ch bond. The solid-state structures of 1,4-bis(methylselenoethynyl)perfluorobenzene ( 1Se ) show the formation of recurrent chalcogen-bonded (ChB) motifs. Association of 1Se and the tellurium analogue 1Te with 4,4′-bipyridine and with the stronger Lewis base 1,4-di(4-pyridyl)piperazine gives 1:1 co-crystals with 1D extended structures linked by short and directional ChB interactions, comparable to those observed with the corresponding halogen bond (XB) donor, 1,4-bis(iodoethynyl)-perfluorobenzene. This “alkynyl” approach for chalcogen activation provides the crystal-engineering community with efficient, and neutral ChB donors for the elaboration of supramolecular 1D (and potentially 2D or 3D) architectures, with a degree of strength and predictability comparable to that of halogen bonding in iodoacetylene derivatives.     

Competition between chalcogen bond and halogen bond interactions in YOX<Subscript>4</Subscript>:NH<Subscript>3</Subscript> (Y = S,Se; X = F,Cl, Br) complexes: An ab initio investigation

Using ab initio calculations, the geometries, interaction energies and bonding properties of chalcogen bond and halogen bond interactions between YOX₄ (Y = S, Se; X = F, Cl, Br) and NH₃ molecules are studied. These binary complexes are formed through the interaction of a positive electrostatic potential region (σ-hole) on the YOX₄ with the negative region in the NH₃. The ab initio calculations are carried out at the MP2/aug-cc-pVTZ level, through analysis of molecular electrostatic potentials, quantum theory of atoms in molecules and natural bond orbital methods. Our results indicate that even though the chalcogen and halogen bonds are mainly dominated by electrostatic effects, but the polarization and dispersion effects also make important contributions to the total interaction energy of these complexes. The examination of interaction energies suggests that the chalcogen bond is always favored over the halogen bond for all of the binary YOX₄:NH₃ complexes.     

Halogen bonding between substituted chlorobenzene and trimethylamine: Decisive role of σ–hole and C?Cl bond breaking energy

Theoretical studies have been carried out on the halogen bonding interaction between para substituted chlorobenzene (Y C₆H₄Cl, Y = H, NH₂, CH₃, F, CN, NO₂) and N(CH₃)₃ using ab initio MP2/aug‐cc‐pVDZ and DFT based wB97XD/6‐311++G(d,p) methods. The positive electrostatic potential (V_S,max) on the Cl atom and the heterolytic bond breaking enthalpy of the C Cl bond have been calculated and their role on halogen bonding is discussed. The heterolytic bond breaking enthalpy of the C Cl bond is proposed as a measure of the strength of the σ‐hole on Cl atom. The binding strength of the complexes ranging between −6.13 kJ mol⁻¹ and −9.29 kJ mol⁻¹ are linearly related to the V_S,max of the Cl atom and the bond breaking enthalpy of the C Cl bond. In addition, energy decomposition analysis was performed on the halogen bonded complexes via symmetry adapted perturbation theory (SAPT) to predict the dominant energy component and the nature of the N···Cl interaction.     

Halogen Bonding： An AIM Analysis of the Weak Interactions

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

Charge-Assisted Chalcogen Bonds: CSD and DFT Analyses and Biological Implication in Glucosidase Inhibitors

This study reports a combined Cambridge Structural Database and theoretical DFT study of charge assisted chalcogen bonds involving sulfonium, selenonium, and telluronium cations. The chalcogen bond has been recently defined by IUPAC as the net attractive interaction between an electrophilic region associated with a chalcogen atom in a molecular entity and a nucleophilic region in another, or the same, molecular entity. Divalent chalcogen atoms typically have up to two σ-holes and forms up to two ChBs; the same holds for tetravalent chalcogens which adopt a seesaw arrangement. In sulfonium, selenonium, and telluronium salts chalcogen atoms form three covalent bonds, three σ-holes are located opposite to these bonds, and up to three charge assisted ChBs can be formed between these holes and the counterions. The covalent bond arrangement around these chalcogen atoms is similar to trivalent pnictogen atoms and translates into a similar pattern of noncovalent interactions. We have found and studied this type of charge-assisted chalcogen bonds in various sulfonium ion-containing inhibitors of glucosidase, for example, salacinol and kotalanol.     

Actinium coordination chemistry: A density functional theory study with monodentate and bidentate ligands

In the current study, the coordination chemistry of nine-coordinate Ac(III) complexes with 35 monodentate and bidentate ligands was investigated using density functional theory (DFT) in terms of their geometries, charges, reaction energies, and bonding interactions. The energy decomposition analysis with naturals orbitals for chemical valence (EDA-NOCV) and the quantum theory of atoms in molecules (QTAIM) were employed as analysis methods. Trivalent Ac exhibits the highest affinities toward hard acids (such as charged oxophilic donors, fluoride), so its classification as a hard acid is justified. Natural population analysis quantified the involvement of 5f orbitals on Ac to be about 30% of total valence electron natural configuration indicating that Ac is a member of the actinide series. Pearson correlation coefficients were used to study the pairwise correlations among the bond lengths, ΔG reaction energies, charges on Ac and donor atoms, and data from EDA-NOCV and QTAIM. Strong correlations and anticorrelations were found between Voronoi charges on donor atoms with ΔG, EDA-NOCV interaction energies and QTAIM bond critical point densities.     

σ-Hole bond tunability in YO<Subscript>2</Subscript>X<Subscript>2</Subscript>:NH<Subscript>3</Subscript> and YO<Subscript>2</Subscript>X<Subscript>2</Subscript>:H<Subscript>2</Subscript>O complexes (X = F,Cl, Br; Y = S,Se): trends and theoretical aspects

A σ-hole is defined as an electron-deficient region on the extension of a covalently bonded group IV–VII atoms. If the electronic density in the σ-hole is sufficiently low, then this region will have a positive electrostatic potential, which allows attractive noncovalent interactions with negative sites. SO₂X₂ and SeO₂X₂ (X = F, Cl and Br) have three Lewis acid sites of σ-hole located in the outermost of chalcogen atom and X end, participating in the chalcogen and halogen bonds with NH₃ and H₂O, respectively. MP2/aug-cc-pVTZ and M06-2X/aug-cc-pVTZ calculations reveal that for a given halogen atom, SeO₂X₂ forms stronger chalcogen bond interactions than SO₂X₂ counterpart. Almost a perfect linear relationship is evident between the interaction energies and the magnitudes of the product of most positive and negative electrostatic potentials. The interaction energies calculated by M06-2X and MP2 methods are almost consistent with each other.     

Revealing substituent effects on the concerted interaction of pnicogen, chalcogen, and halogen bonds in substituted s-triazine ring

Quantum chemical calculations have been performed to gauge the effect of substituents on concerted interactions of pnicogen, chalcogen, and halogen bonds in the X–TAZ···Y complexes (X = CN, F, Cl, Br, H, CH₃, OH, and NH₂, where TAZ and Y denote s-triazine ring and P, S, and Cl atoms, respectively) at the M06-2X/aug-cc-pVDZ level. The mutual interplay of these interactions is also investigated. The results indicate that diminutive effects are observed when the three kinds of noncovalent interactions pnicogen, chalcogen, and halogen bonds are coexisted in the complexes. These effects are studied in terms of energetic and geometric features of the complexes. In addition, Bader’s theory of “atoms in molecules” is used to analyze their strength of varying electron density at bond critical points. Natural bond orbital (NBO) theory is used to characterize the orbital interactions. The results indicate that the electron-withdrawing/donating substituents decrease/increase the magnitude of the binding energies compared to the unsubstituted X–TAZ···Y (X = H) complex. Good correlations among binding energies, Hammett constants, geometrical, atoms in molecular and NBO parameters are established in X–TAZ···Y complexes. By taking advantage of all the aforementioned computational methods, this study examines how these interactions mutually influence each other.     

Directionality of Halogen Bonds: An Interacting Quantum Atoms (IQA) and Relative Energy Gradient (REG) Study

Interacting Quantum Atoms (IQA) and Interacting Quantum Fragments (IQF) analyses are used to study (X=Cl and Br) model complexes in order to determine the origin of halogen bond directionality. IQA allows for the calculation of intra- and interatomic classical and exchange-correlation energies, which can be used to determine the energetic nature of the changes that occur when deviating from the preferred halogen bond approach. The Relative Energy Gradient (REG) method is also applied to rank the IQA energies and reveal which energy contributions best describe the total behavior of the system. Indeed, all the pairwise interactions and atomic self-energies are angularly dependent; some terms favor the linear structure and some tend toward nonlinear arrangements. For instance, when the C−X−N angle is altered, the halogen-nitrogen interaction energy behaves like the total energy of the system while the carbon-nitrogen interaction works against the total energy profile. Furthermore, the REG values reveal that the contribution of the halogen-nitrogen interaction to the total behavior of the system is small. Instead, the secondary interactions (e. g., fluorine-nitrogen and carbon-hydrogen interactions) and atomic self-energies are mainly responsible for the angular preference of these halogen bonds. Finally, IQF calculations followed by REG analysis reveal the importance of the self-energy of the fragments.     

Halogen-hydride interaction between Z-X (Z = CN, NC; X = F, Cl, Br) and H-Mg-Y (Y = H, F, Cl, Br, CH3)

Halogen-hydride interactions between Z-X (Z = CN, NC and X = F, Cl, Br) as halogen donor and H-Mg-Y (Y = H, F, Cl, Br, CH(3)) as electron donor have been investigated through the use of Becke three-parameter hybrid exchange with Lee-Yang-Parr correlation (B3LYP), second-order M?ller-Plesset perturbation theory (MP2), and coupled-cluster single and double excitation (with triple excitations) [CCSD(T)] approaches. Geometry changes during the halogen-hydride interaction are accompanied by a mutual polarization of both partners with some charge transfer occurring from the electron donor subunit. Interaction energies computed at MP2 level vary from -1.23 to -2.99 kJ/mol for Z-F···H-Mg-Y complexes, indicating that the fluorine interactions are relatively very weak but not negligible. Instead, for chlorine- and bromine-containing complexes the interaction energies span from -5.78 to a maximum of -26.42 kJ/mol, which intimate that the interactions are comparable to conventional hydrogen bonding. Moreover, the calculated interaction energy was found to increase in magnitude with increasing positive electrostatic potential on the extension of Z-X bond. Analysis of geometric, vibrational frequency shift and the interaction energies indicates that, depending on the halogen, CN-X···H interactions are about 1.3-2.0 times stronger than NC-X···H interactions in which the halogen bonds to carbon. We also identified a clear dependence of the halogen-hydride bond strength on the electron-donating or -withdrawing effect of the substituent in the H-Mg-Y subunits. Furthermore, the electronic and structural properties of the resulting complexes have been unveiled by means of the atoms in molecules (AIM) and natural bond orbital (NBO) analyses. Finally, several correlative relationships between interaction energies and various properties such as binding distance, frequency shift, molecular electrostatic potential, and intermolecular density at bond critical point have been checked for all studied systems.     

The Role of Molecular Electrostatic Potentials in the Formation of a Halogen Bond in Furan⋅⋅⋅XY and Thiophene⋅⋅⋅XY Complexes

The halogen bonding of furan???XY and thiophene???XY (X=Cl, Br; Y=F, Cl, Br), involving σ‐ and π‐type interactions, was studied by using MP2 calculations and quantum theory of “atoms in molecules” (QTAIM) studies. The negative electrostatic potentials of furan and thiophene, as well as the most positive electrostatic potential (V_S,max) on the surface of the interacting X atom determined the geometries of the complexes. Linear relationships were found between interaction energy and V_S,max of the X atom, indicating that electrostatic interactions play an important role in these halogen‐bonding interactions. The halogen‐bonding interactions in furan???XY and thiophene???XY are weak, “closed‐shell” noncovalent interactions. The linear relationship of topological properties, energy properties, and the integration of interatomic surfaces versus V_S,max of atom X demonstrate the importance of the positive σ hole, as reflected by the computed V_S,max of atom X, in determining the topological properties of the halogen bonds.     

Chalcogen Bonding Catalysis of a Nitro‐Michael Reaction

Chalcogen bonding is the non‐covalent interaction between Lewis acidic chalcogen substituents and Lewis bases. Herein, we present the first application of dicationic tellurium‐based chalcogen bond donors in the nitro‐Michael reaction between trans‐β‐nitrostyrene and indoles. This also constitutes the first activation of nitro derivatives by chalcogen bonding (and halogen bonding). The catalysts showed rate accelerations of more than a factor of 300 compared to strongly Lewis acidic hydrogen bond donors. Several comparison experiments, titrations, and DFT calculations support a chalcogen‐bonding‐based mode of activation of β‐nitrostyrene.     

A Quantitative Molecular Orbital Perspective of the Chalcogen Bond

We have quantum chemically analyzed the structure and stability of archetypal chalcogen-bonded model complexes D₂Ch⋅⋅⋅A⁻ (Ch = O, S, Se, Te; D, A = F, Cl, Br) using relativistic density functional theory at ZORA-M06/QZ4P. Our purpose is twofold: (i) to compute accurate trends in chalcogen-bond strength based on a set of consistent data; and (ii) to rationalize these trends in terms of detailed analyses of the bonding mechanism based on quantitative Kohn-Sham molecular orbital (KS-MO) theory in combination with a canonical energy decomposition analysis (EDA). At odds with the commonly accepted view of chalcogen bonding as a predominantly electrostatic phenomenon, we find that chalcogen bonds, just as hydrogen and halogen bonds, have a significant covalent character stemming from strong HOMO−LUMO interactions. Besides providing significantly to the bond strength, these orbital interactions are also manifested by the structural distortions they induce as well as the associated charge transfer from A⁻ to D₂Ch.     

Halogen Bonding Propensity in Solution: Direct Observation and Computational Prediction

Halogen-bonded complexes are often designed by consideration of electrostatic potential (ESP) predictions. ESP predictions do not capture the myriad variables associated with halogen bond (XB) donors and acceptors; thus, binding interaction cannot be quantitatively predicted. Here, a discrepancy between predictions based on ESP energy difference (ΔV_s) and computed gas phase binding energy (ΔE_bind) motivated the experimental determination of the relative strength of halogen bonding interactions in solution by Raman spectroscopic observation of complexes formed from interacting five iodobenzene-derived XB donors and four pyridine XB acceptors. Evaluation of ΔE_bind coupled with absolutely-localized molecular orbital energy decomposition analysis (ALMO-EDA) deconvolutes halogen bonding energy contributions and reveals a prominent role for charge transfer (CT) interactions. Raman spectra reveal ΔE_bind accurately predicts stronger interactions within iodopentafluorobenzene (IPFB) complexes than with 1-iodo-3,5-dinitrobenzene (IDNB) complexes even though IPFB has similar electrostatics to IDNB and contains a smaller σ-hole.     

Molecular Dynamics Simulation Study of ClF in Water: Halogen Bonding Interaction in Liquid

Does the halogen bonding interaction co-exist in liquid when it competes with the hydrogen bonding interaction? The classical molecular dynamics simulations for the solvation properties of ClF molecule in water are performed with the Lennard-Jones plus Coulomb electrostatic potential parameters that are optimized with ab initio interaction energy calculations for the pre-reactive H₂O…ClF complex. We find that the halogen bonding interactions occur between O and Cl atoms and have the comparable strength and population with respect to the hydrogen bonding interactions of Cl…H.     

The dual role of halogen,chalcogen, and pnictogen atoms as Lewis acid and base: Triangular XBr:SHX:PH2X complexes (X = F,Cl, Br,CN, NC,OH, NH2, and OCH3)

Using ab initio calculations, we have investigated the possibility of formation of triangular XBr:SHX:PH₂X complexes, where X = F, Cl, Br, CN, NC, OH, NH₂, and OCH₃. These complexes are formed through the interaction of a positive electrostatic potential region (σ‐hole) on a molecule with the negative region in another one. The results show that the combined halogen, chalcogen, and pnictogen interactions can give rise to stable cyclic structures. The interaction energies of these complexes span over a wide range, from ?3.55 to ?24.93 kcal/mol. Nice quadratic correlations are found between the interaction energies and binding distances in the trimers. To understand the nature of the interactions in these complexes, molecular electrostatic potential and quantum theory of atoms in molecule analyses are performed. © 2015 Wiley Periodicals, Inc.     

Theoretical Density Functional Theory insights into the nature of chalcogen bonding between CX2 (X = S,Se, Te) and diazine from monomer to supramolecular complexes

Chalcogen bonding is a noncovalent interaction, highly similar to halogen and hydrogen bonding, occurring between a chalcogen atom and a nucleophilic region. Two density functional theory (DFT) approaches B3LY-D3 and B97-D3 were performed on a series of complexes formed between CX₂ (X = S, Se, Te) and diazine (pyridazine, pyrimidine and pyrazine). Chalcogen atoms prefer interacting with the lone pair of a nitrogen atom rather than with the π-cloud of an aromatic ring. CTe₂ and CSe₂ form a stronger chalcogen bond than CS₂. The electrostatic potential of CX₂ (X = S, Se and Te) reveals the presence of two equivalent σ-holes, one on each chalcogen atom. These CX₂ molecules interact with diazine giving rise to supramolecular interactions. Wiberg bond index and second-order perturbation theory analysis in NBO were performed to better understand the nature of the chalcogen bond interaction.     

Theoretical investigation on charge‐assisted halogen bonding interactions in the complexes of bromocarbons with some anions

A series of dimeric complexes formed between bromocarbon molecules and two anions (Br^? and CN^?) have been investigated by using MP2 method. The quantum theory of atoms in molecules (QTAIM) and the second‐order perturbation natural bond orbital (NBO) approaches were applied to analyze the electron density distributions of these complexes and to explore the nature of charge‐assisted halogen bonding interactions. As anticipated, these interactions are significantly stronger relative to the corresponding neutral ones. The results derived from ab initio calculations described herein reveal a major contribution from the electrostatic interaction on the stability of the systems considered. Beside the electrostatic interaction, the charge‐transfer force and the second‐order orbital interaction also play an important role in the formation of the complexes, as a NBO analysis suggested. The presence of halogen bonds in the complexes has been identified in terms of the QTAIM methodology, and several linear relationships have been established to provide more insight into charge‐assisted halogen bonding interactions. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008