Is There a Single Ideal Parameter for Halogen-Bonding-Based Lewis Acidity?

A comparison of structure, bonding and non-covalent interactions of aryl halide and diarylhalonium halogen-bond donors

Halogen bonding permeates many areas of chemistry. A wide range of halogen-bond donors including neutral, cationic, monovalent, and hypervalent have been developed and studied. In this work we used density functional theory (DFT), natural bond orbital (NBO) theory, and quantum theory of atoms in molecules (QTAIM) to analyze aryl halogen-bond donors that are neutral, cationic, monovalent and hypervalent and in each series we include the halogens Cl, Br, I, and At. Within this diverse set of halogen-bond donors, we have found trends that relate halogen bond length with the van der Waals radii of the halogen and the non-covalent or partial covalency of the halogen bond. We have also developed a model to calculate ΔG of halogen-bond formation by the linear combination of the % p-orbital character on the halogen and energy of the σ-hole on the halogen-bond donor.

Gauging the Strength of the Molecular Halogen Bond via Experimental Electron Density and Spectroscopy

Strong and weak halogen bonds (XBs) in discrete aggregates involving the same acceptor are addressed by experiments in solution and in the solid state. Unsubstituted and perfluorinated iodobenzenes act as halogen donors of tunable strength; in all cases, quinuclidine represents the acceptor. NMR titrations reliably identify the strong intermolecular interactions in solution, with experimental binding energies of approx. 7 kJ/mol. Interaction of the σ hole at the halogen donor iodine leads to a redshift in the symmetric C-I stretching vibration; this shift reflects the interaction energy in the halogen-bonded adducts and may be assessed by Raman spectroscopy in condensed phase even for weak XBs. An experimental picture of the electronic density for the XBs is achieved by high-resolution X-ray diffraction on suitable crystals. Quantum theory of atoms in molecules (QTAIM) analysis affords the electron densities and energy densities in the bond critical points of the halogen bonds and confirms stronger interaction for the shorter contacts. For the first time, the experimental electron density shows a significant effect on the atomic volumes and Bader charges of the quinuclidine N atoms, the halogen-bond acceptor: strong and weak XBs are reflected in the nature of their acceptor atom. Our experimental findings at the acceptor atom match the discussed effects of halogen bonding and thus the proposed concepts in XB activated organocatalysis.

A Quantum-chemical Analysis on the Lewis Acidity of Diarylhalonium Ions

Cyclic diaryliodonium compounds like iodolium derivatives have increasingly found use as noncovalent Lewis acids in the last years. They are more stable toward nucleophilic substitution than acyclic systems and are markedly more Lewis acidic. Herein, this higher Lewis acidity is analyzed and explained via quantum-chemical calculations and energy decomposition analyses. Its key origin is the change in energy levels and hybridization of iodine's orbitals, leading to both more favorable electrostatic interaction and better charge transfer. Both of the latter seem to contribute in similar fashion, while hydrogen bonding as well as steric repulsion with the phenyl rings play at best a minor role. In comparison to iodolium, bromolium and chlorolium are less Lewis acidic the lighter the halogen, which is predominantly based on less favorable charge-transfer interactions.

Catalytic Activation of Imines by Chalcogen Bond Donors in a Povarov [4+2] Cycloaddition Reaction

Recently, chalcogen bonding has been investigated in more detail in organocatalysis and the scope of activated functionalities continues to increase. Herein, the activation of imines in a Povarov [4+2] cycloaddition reaction with bidentate cationic chalcogen bond donors is presented. Tellurium-based Lewis acids show superior properties compared to selenium-based catalysts and inactive sulfur-based analogues. The catalytic activity of the chalcogen bonding donors increases with weaker binding anions. Triflate, however, is not suitable due to its participation in the catalytic pathway. A solvent screening revealed a more efficient activation in less polar solvents and a pronounced effect of solvent (and catalyst) on endo : exo diastereomeric ratio. Finally, new chiral chalcogen bonding catalysts were applied but provided only racemic mixtures of the product.

Orbital analysis of bonding in diarylhalonium salts and relevance to periodic trends in structure and reactivity

Diarylhalonium compounds provide new opportunities as reagents and catalysts in the field of organic synthesis. The three center, four electron (3c-4e) bond is a center piece of their reactivity, but structural variation among the diarylhaloniums, and in comparison with other λ3-iodanes, indicates that the model needs refinement for broader applicability. We use a combination of Density Functional Theory (DFT), Natural Bond Orbital (NBO) Theory, and X-ray structure data to correlate bonding and structure for a λ3-iodane and a series of diarylchloronium, bromonium, and iodonium salts, and their isoelectronic diarylchalcogen counterparts. This analysis reveals that the s-orbital on the central halogen atom plays a greater role in the 3c-4e bond than previously considered. Finally, we show that our revised bonding model and associated structures account for both kinetic and thermodynamic reactivity for both acyclic phenyl(mesityl)halonium and cyclic dibenzohalolium salts.

Revealing Intra- and Intermolecular Interactions Determining Physico-Chemical Features of Selected Quinolone Carboxylic Acid Derivatives

The intra- and intermolecular interactions of selected quinolone carboxylic acid derivatives were studied in monomers, dimers and crystals. The investigated compounds are well-recognized as medicines or as bases for further studies in drug design. We employed density functional theory (DFT) in its classical formulation to develop gas-phase and solvent reaction field (PCM) models describing geometric, energetic and electronic structure parameters for monomers and dimers. The electronic structure was investigated based on the atoms in molecules (AIM) and natural bond orbital (NBO) theories. Special attention was devoted to the intramolecular hydrogen bonds (HB) present in the investigated compounds. The characterization of energy components was performed using symmetry-adapted perturbation theory (SAPT). Finally, the time-evolution methods of Car-Parrinello molecular dynamics (CPMD) and path integral molecular dynamics (PIMD) were employed to describe the hydrogen bond dynamics as well as the spectroscopic signatures. The vibrational features of the O-H stretching were studied using Fourier transformation of the autocorrelation function of atomic velocity. The inclusion of quantum nuclear effects provided an accurate depiction of the bridged proton delocalization. The CPMD and PIMD simulations were carried out in the gas and crystalline phases. It was found that the polar environment enhances the strength of the intramolecular hydrogen bonds. The SAPT analysis revealed that the dispersive forces are decisive factors in the intermolecular interactions. In the electronic ground state, the proton-transfer phenomena are not favourable. The CPMD results showed generally that the bridged proton is localized at the donor side, with possible proton-sharing events in the solid-phase simulation of stronger hydrogen bridges. However, the PIMD enabled the quantitative estimation of the quantum effects inclusion-the proton position was moved towards the bridge midpoint, but no qualitative changes were detected. It was found that the interatomic distance between the donor and acceptor atoms was shortened and that the bridged proton was strongly delocalized.

Theoretical treatment of IO-X (X = N2, CO, CO2, H2O) complexes

Iodine monoxide (IO) is an important component of the biogeochemical cycle of iodine. For instance, it is present in the troposphere, where it plays a crucial role in the physical chemical processes involving iodine containing compounds. Here, we present a theoretical study on a series of atmospherically relevant complexes of IO with N₂, CO, CO₂ and H₂O, where their structural and spectroscopic properties and their interaction energies are computed. Calculations are carried out by means of ab initio post Hartree-Fock (RCCSD(T) and RMP2) methods and density functional theory DFT (PBE0 and M05-2X) based approaches with and without the inclusion of dispersion correction. After comparison to RCCSD(T), we highlight the good performance of M05-2X(+D3) DFT in describing the bonding between IO and X (X = N₂, CO, CO₂, H₂O). Moreover, we found that the IO-X (X = N₂, CO, CO₂, H₂O) complexes are formed by non-covalent interactions between the two monomers. In sum, we characterized two types of complexes: I-bonded and O-bonded, where the former is more stable. The atmospheric implications of the present findings are also discussed such as in the formation of the iodine oxide particles (IOPs).

Pushing the Limits of the Hydrogen Bond Enhanced Halogen Bond —The Case of the C–H Hydrogen Bond

C–H hydrogen bonds have remarkable impacts on various chemical systems. Here we consider the influence of C–H hydrogen bonds to iodine atoms. Positioning a methyl group between two iodine halogen bond donors of the receptor engendered intramolecular C–H hydrogen bonding (HBing) to the electron-rich belt of both halogen bond donors. When coupled with control molecules, the role of the C–H hydrogen bond was evaluated. Gas-phase density functional theory studies indicated that methyl C–H hydrogen bonds help bias a bidentate binding conformation. Interaction energy analysis suggested that the charged C–H donors augment the halogen bond interaction—producing a >10 kcal mol⁻¹ enhancement over a control lacking the C–H⋯I–C interaction. X-ray crystallographic analysis demonstrated C–H hydrogen bonds and bidentate conformations with triflate and iodide anions, yet the steric bulk of the central functional group seems to impact the expected trends in halogen bond distance. In solution, anion titration data indicated elevated performance from the receptors that utilize C–H Hydrogen Bond enhanced Halogen Bonds (HBeXBs). Collectively, the results suggest that even modest hydrogen bonds between C–H donors and iodine acceptors can influence molecular structure and improve receptor performance.

C–H hydrogen bonds to iodine halogen bond donors are shown to improve halogen bonding and molecular preorganization.

Maximal occupation by bases of π-hole bands surrounding linear molecules

Linear molecules such as CO₂ contain a positive π-hole ring that surrounds C on the molecule's equator. Quantum calculations examine the question as to how many bases can simultaneously bind to this ring. Linear molecules examined are TO₂ , where T = C, Si, Ge, Sn; bases are NCH and NH₃ . CO₂ engages in the weakest of the tetrel bonds, and can bind up to three NCH and two NH₃ . Unlike σ-hole tetrel bonds, Si forms the strongest tetrel bonds, with interaction energies as high as 43 kcal/mol with NH₃ . But like GeO₂ , SiO₂ can sustain only two bases in its equatorial ring. The π-hole ring of SnO₂ can engage in up to four tetrel bonds with either NCH or NH₃ , even though these bonds are weaker than those with GeO₂ or SiO₂ . As all of these complexes cast TO₂ in the role of multiple electron acceptor, the resulting negative cooperativity makes each successive bond weaker than its predecessor as bases are added, as well as reducing the magnitude of the central molecule's π-hole.

Anatomy of π-hole bonds: Linear systems

The list of σ-hole bonds is long and growing, encompassing both H-bonds and its closely related halogen, chalcogen, etc., sisters. These bonds rely on the asymmetric distribution of electron density, whose depletion along the extension of a covalent bond leaves a positive region of electrostatic potential from which these bonds derive their name. However, the density distributions of other molecules contain analogous positive regions that lie out of the molecular plane known as π-holes, which are likewise capable of engaging in noncovalent bonds. Quantum calculations are applied to study such π-hole bonds that involve linear molecules, whose positive region is a circular belt surrounding the molecule, rather than the more restricted area of a σ-hole. These bonds are examined in terms of their most fundamental elements arising from the spatial dispositions of their relevant molecular orbitals and the π-holes in both the total electron density and the electrostatic potential to which they lead. Systems examined comprise tetrel, chalcogen, aerogen, and triel bonds, as well as those involving group II elements, with atoms drawn from various rows of the Periodic Table. The π-hole bonds established by linear molecules tend to be weaker than those of comparable planar systems.

Yet another perspective on hole interactions

Hole interactions are known by different names depending on the key atom of the bond (halogen bond, chalcogen bond, hydrogen bond, etc.), and the geometry of the interaction (σ if in line, π if perpendicular to the Lewis acid plane). However, its origin starts with the creation of a Lewis acid by an underlying covalent bond, which forms an electrostatic depletion and a virtual antibonding orbital, which can create non-covalent interactions with Lewis bases. In this (maybe subjective) perspective, we will claim that hole interactions must be defined via the molecular orbital origin of the molecule. Under this premise we can better explore the richness of such bonding patterns. For that, we will study old, recent and new systems, trying to pinpoint some misinterpretations that are often associated with them. We will use as exemplars the triel bonds, a couple of metal complexes, a discussion on convergent σ-holes, and many cases of anti-electrostatic hole interactions.

Carbonyl Hypoiodites as Extremely Strong Halogen Bond Donors

Neutral halogen-bonded O-I-N complexes were prepared from in situ formed carbonyl hypoiodites and aromatic organic bases. The carbonyl hypoiodites have a strongly polarized iodine atom with larger σ-holes than any known uncharged halogen bond donor. Modulating the Lewis basicity of the selected pyridine derivatives and carboxylates leads to halogen-bonded complexes where the classical O-I⋅⋅⋅N halogen bond transforms more into a halogen-bonded COO- ⋅⋅⋅I-N+ ion-pair (salt) with an asymmetric O-I-N moiety. X-ray analyses, NMR studies, and calculations reveal the halogen bonding geometries of the carbonyl hypoiodite-based O-I-N complexes, confirming that in the solid-state the iodine atom is much closer to the N-atom of the pyridine derivatives than its original position at the carboxylate O-atom.

Relating Bond Strength and Nature to the Thermodynamic Stability of Hypervalent Togni-Type Iodine Compounds

The bond strength and nature of a set of 32 Togni-like reagents have been investigated at the M062X/def2-TZVP(D) level of theory in acetonitrile described with the SMD continuum solvent model, to rationalize the main factors responsible for their thermodynamic stability in different conformations, and trifluoromethylation capabilities. For the assessment of bond strength, we utilized local stretching force constants and associated bond strength orders, complemented with local features of the electron density to access the nature of the bonds. Bond dissociation energies varied from 31.6 to 79.9 kcal/mol depending on the polarizing power of the ligand trans to CF3 . Based on the analysis of the Laplacian of the density, we propose that the charge-shift bond character plays an important role in the stability of the molecules studied, especially for those containing I-O bonds. New insights on the trans influence and on possible ways to fine-tune the stability of these reagents are provided.

Dissection of the Origin of π-Holes and the Noncovalent Bonds in Which They Engage

Accompanying the rapidly growing list of σ-hole bonds has come the acknowledgment of parallel sorts of noncovalent bonds which owe their stability in large part to a deficiency of electron density in the area above the molecular plane, known as a π-hole. The origins of these π-holes are probed for a wide series of molecules, comprising halogen, chalcogen, pnicogen, tetrel, aerogen, and spodium bonds. Much like in the case of their σ-hole counterparts, formation of the internal covalent π-bond in the Lewis acid molecule pulls density toward the bond midpoint and away from its extremities. This depletion of density above the central atom is amplified by an electron-withdrawing substituent. At the same time, the amplitude of the π*-orbital is enhanced in the region of the density-depleted π-hole, facilitating a better overlap with the nucleophile's lone pair orbital and a stabilizing n → π* charge transfer. The presence of lone pairs on the central atom acts to attenuate the π-hole and shift its position somewhat, resulting in an overall weakening of the π-hole bond. There is a tendency for π-hole bonds to include a higher fraction of induction energy than σ-bonds with proportionately smaller electrostatic and dispersion components, but this distinction is less a product of the σ- or π-character and more a function of the overall bond strength.

Iodine(III)-Based Halogen Bond Donors: Properties and Applications

Halogen bonding, the non-covalent interaction of Lewis bases with an electron-deficient region of halogen substituents, received increased attention recently. Consequently, the design and evaluation of numerous halogen-containing species as halogen bond donors have been subject to intense research. More recently, organoiodine compounds at the iodine(III) state have been receiving growing attention in the field. Due to their electronic and structural properties, they provide access to unique binding modes. For this reason, our groups have been involved in the development of such compounds, in the quantification of their halogen bonding strength (through the evaluation of their Lewis acidities), as well as in the evaluation of their activities as catalysts in several model reactions. This account will describe these contributions.

Multidimensional Lewis Acidity: A Consistent Data Set of Chloride, Hydride, Methide, Water and Ammonia Affinities for 183 p-Block Element Lewis Acids

The computed fluoride ion affinity (FIA) is a widely applied descriptor to gauge Lewis acidity. Like every other single-parameter Lewis acidity scale, the FIA metric suffers from the one-dimensionality, that prohibits addressing Lewis acidity by the multidimensionality it inherently requires (i. e., reference Lewis base dependency). However, a systematic screening of computed affinities other than the FIA is much less developed. Herein, we extended our CCSD(T)/CBS benchmark of different density functionals and the DLPNO-CCSD(T) method for chloride (CIA), methide (MIA), hydride (HIA), water (WA), and ammonia (AA) affinities. The best performing methods are subsequently applied to yield nearly 800 affinities for 183 p-block element compounds of group 13-16 with an estimated accuracy of <10 kJ mol-1 . The study's output serves as a consistent library for qualitative analyses and a training set for future statistical approaches. A first holistic correlation analysis underscores the need for a multidimensional description of Lewis acidity.

Halogen Bonding Involving I2 and d8 Transition-Metal Pincer Complexes

We systematically investigated iodine–metal and iodine–iodine bonding in van Koten’s pincer complex and 19 modifications changing substituents and/or the transition metal with a PBE0–D3(BJ)/aug–cc–pVTZ/PP(M,I) model chemistry. As a novel tool for the quantitative assessment of the iodine–metal and iodine–iodine bond strength in these complexes we used the local mode analysis, originally introduced by Konkoli and Cremer, complemented with NBO and Bader’s QTAIM analyses. Our study reveals the major electronic effects in the catalytic activity of the M–I–I non-classical three-center bond of the pincer complex, which is involved in the oxidative addition of molecular iodine I2 to the metal center. According to our investigations the charge transfer from the metal to the σ* antibonding orbital of the I–I bond changes the 3c–4e character of the M–I–I three-center bond, which leads to weakening of the iodine I–I bond and strengthening of the metal–iodine M–I bond, facilitating in this way the oxidative addition of I2 to the metal. The charge transfer can be systematically modified by substitution at different places of the pincer complex and by different transition metals, changing the strength of both the M–I and the I2 bonds. We also modeled for the original pincer complex how solvents with different polarity influence the 3c–4e character of the M–I–I bond. Our results provide new guidelines for the design of pincer complexes with specific iodine–metal bond strengths and introduce the local vibrational mode analysis as an efficient tool to assess the bond strength in complexes.

A Bidentate Iodine(III)-Based Halogen-Bond Donor as a Powerful Organocatalyst*

In contrast to iodine(I)-based halogen bond donors, iodine(III)-derived ones have only been used as Lewis acidic organocatalysts in a handful of examples, and in all cases they acted in a monodentate fashion. Herein, we report the first application of a bidentate bis(iodolium) salt as organocatalyst in a Michael and a nitro-Michael addition reaction as well as in a Diels-Alder reaction that had not been activated by noncovalent organocatalysts before. In all cases, the performance of this bidentate XB donor distinctly surpassed the one of arguably the currently strongest iodine(I)-based organocatalyst. Bidentate coordination to the substrate was corroborated by a structural analysis and by DFT calculations of the transition states. Overall, the catalytic activity of the bis(iodolium) system approaches that of strong Lewis acids like BF3 .

A low cost, high accuracy method for halogen bonding complexes

The ONIOM scheme M052X/[Def2TZVP+Def2TZVPD.ECP(I)]:AM1 is shown to represent halogen bond (XB) geometries nearly as well as DFT while being more than two orders of magnitude faster in systems containing >40 atoms. This finding is shown to hold for 40 XB donors, which cover most known backbones, and for a range of neutral and anionic Lewis bases. Complexation free energies can be accurately computed using these geometries and a single-point energy calculation at the DFT level. This approach circumvents the unfavorable scaling of computing time associated with modeling large systems involving halogen bonding.

Halogen Bond of Halonium Ions: Benchmarking DFT Methods for the Description of NMR Chemical Shifts

Because of their anisotropic electron distribution and electron deficiency, halonium ions are unusually strong halogen-bond donors that form strong and directional three-center, four-electron halogen bonds. These halogen bonds have received considerable attention owing to their applicability in supramolecular and synthetic chemistry and have been intensely studied using spectroscopic and crystallographic techniques over the past decade. Their computational treatment faces different challenges to those of conventional weak and neutral halogen bonds. Literature studies have used a variety of wave functions and DFT functionals for prediction of their geometries and NMR chemical shifts, however, without any systematic evaluation of the accuracy of these methods being available. In order to provide guidance for future studies, we present the assessment of the accuracy of 12 common DFT functionals along with the Hartree-Fock (HF) and the second-order Møller-Plesset perturbation theory (MP2) methods, selected from an initial set of 36 prescreened functionals, for the prediction of 1H, 13C, and 15N NMR chemical shifts of [N-X-N]+ halogen-bond complexes, where X = F, Cl, Br, and I. Using a benchmark set of 14 complexes, providing 170 high-quality experimental chemical shifts, we show that the choice of the DFT functional is more important than that of the basis set. The M06 functional in combination with the aug-cc-pVTZ basis set is demonstrated to provide the overall most accurate NMR chemical shifts, whereas LC-ωPBE, ωB97X-D, LC-TPSS, CAM-B3LYP, and B3LYP to show acceptable performance. Our results are expected to provide a guideline to facilitate future developments and applications of the [N-X-N]+ halogen bond.

Hypervalent Iodine(III) Compounds as Biaxial Halogen Bond Donors

“Hypervalent” iodine(III) derivatives have been established as powerful reagents in organic transformations, but so far only a handful of studies have addressed their potential use as halogen-bonding noncovalent Lewis acids. In contrast to “classical” halogen-bond donors based on iodine(I) compounds, iodine(III) salts feature two directional electrophilic axes perpendicular to each other. Herein we present the first systematic investigation on biaxial binding to such Lewis acids in solution. To this end, hindered and unhindered iodolium species were titrated with various substrates, including diesters and diamides, via ¹H NMR spectroscopy and isothermal titration calorimetry. Clear evidence for biaxial binding was obtained in two model systems, and the association strengths increased by 2 orders of magnitude. These findings were corroborated by density functional theory calculations (which reproduced the trend well but underestimated the absolute binding constants) and a cocrystal featuring biaxial coordination of a diamide to the unhindered iodolium compound.

Is There a Single Ideal Parameter for Halogen-Bonding-Based Lewis Acidity?

Halogen-bond donors (halogen-based Lewis acids) have now found various applications in diverse fields of chemistry. The goal of this study was to identify a parameter obtainable from a single DFT calculation that reliably describes halogen-bonding strength (Lewis acidity). First, several DFT methods were benchmarked against the CCSD(T) CBS binding data of complexes of 17 carbon-based halogen-bond donors with chloride and ammonia as representative Lewis bases, which revealed M05-2X with a partially augmented def2-TZVP(D) basis set as the best model chemistry. The best single parameter to predict halogen-bonding strengths was the static σ-hole depth, but it still provided inaccurate predictions for a series of compounds. Thus, a more reliable parameter, Ωσ* , has been developed through the linear combination of the σ-hole depth and the σ*(C-I) energy, which was further validated against neutral, cationic, halogen- and nitrogen-based halogen-bond donors with very good performance.