Halogen Bonding: A Halogen-Centered Noncovalent Interaction Yet to Be Understood

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.

On the Potentiality of X-T-X3 Compounds (T = C, Si, and Ge, and X = F, Cl, and Br) as Tetrel- and Halogen-Bond Donors

The versatility of the X-T-X3 compounds (where T = C, Si, and Ge, and X = F, Cl, and Br) to participate in tetrel- and halogen-bonding interactions was settled out, at the MP2/aug-cc-pVTZ level of theory, within a series of configurations for (X-T-X3)2 homodimers. The electrostatic potential computations ensured the remarkable ability of the investigated X-T-X3 monomers to participate in σ-hole halogen and tetrel interactions. The energetic findings significantly unveil the favorability of the tetrel···tetrel directional configuration with considerable negative binding energies over tetrel···halogen, type III halogen···halogen, and type II halogen···halogen analogs. Quantum theory of atoms in molecules and noncovalent interaction analyses were accomplished to disclose the nature of the tetrel- and halogen-bonding interactions within designed configurations, giving good correlations between the total electron densities and binding energies. Further insight into the binding energy physical meanings was invoked through using symmetry-adapted perturbation theory-based energy decomposition analysis, featuring the dispersion term as the most prominent force beyond the examined interactions. The theoretical results were supported by versatile crystal structures which were characterized by the same type of interactions. Presumably, the obtained findings would be considered as a solid underpinning for future supramolecular chemistry, materials science, and crystal engineering studies, as well as a fundamental linchpin for a better understanding of the biological activities of chemicals.

Comparative Structural Study of Three Tetrahalophthalic Anhydrides: Recognition of X···O(anhydride) Halogen Bond and πh···O(anhydride) Interaction

Structures of three tetrahalophthalic anhydrides (TXPA: halogen = Cl (TCPA), Br (TBPA), I (TIPA)) were studied by X-ray diffraction, and several types of halogen bonds (HaB) and lone pair···π-hole (lp···πh) contacts were revealed in their structures. HaBs involving the central oxygen atom of anhydride group (further X···O(anhydride) were recognized in the structures of TCPA and TBPA. In contrast, for the O(anhydride) atom of TIPA, only interactions with the π system (π-hole) of the anhydride ring (further lp(O)···πh) were observed. Computational studies by a number of theoretical methods (molecular electrostatic potentials, the quantum theory of atoms in molecules, the independent gradient model, natural bond orbital analyses, the electron density difference, and symmetry-adapted perturbation theory) demonstrated that the X···O(anhydride) contacts in TCPA and TBPA and lp(O)···πh in TIPA are caused by the packing effect. The supramolecular architecture of isostructural TCPA and TBPA was mainly affected by X···O(acyl) and X···X HaBs, and, for TIPA, the main contribution provided I···I HaBs.

Study of Beryllium, Magnesium, and Spodium Bonds to Carbenes and Carbodiphosphoranes

The aim of this article is to present results of theoretical study on the properties of C⋯M bonds, where C is either a carbene or carbodiphosphorane carbon atom and M is an acidic center of MX2 (M = Be, Mg, Zn). Due to the rarity of theoretical data regarding the C⋯Zn bond (i.e., the zinc bond), the main focus is placed on comparing the characteristics of this interaction with C⋯Be (beryllium bond) and C⋯Mg (magnesium bond). For this purpose, theoretical studies (ωB97X-D/6-311++G(2df,2p)) have been performed for a large group of dimers formed by MX2 (X = H, F, Cl, Br, Me) and either a carbene ((NH2)2C, imidazol-2-ylidene, imidazolidin-2-ylidene, tetrahydropyrymid-2-ylidene, cyclopropenylidene) or carbodiphosphorane ((PH3)2C, (NH3)2C) molecule. The investigated dimers are characterized by a very strong charge transfer effect from either the carbene or carbodiphosphorane molecule to the MX2 one. This may even be over six times as strong as in the water dimer. According to the QTAIM and NCI method, the zinc bond is not very different than the beryllium bond, with both featuring a significant covalent contribution. However, the zinc bond should be definitely stronger if delocalization index is considered.

1,12-Diiodo-Ortho-Carborane: A Classic Textbook Example of the Dihalogen Bond

The crystal structure of 1,12-diiodo-ortho-carborane 1,12-I2-1,2-C2B10H10 was determined by single crystal X-ray diffraction. In contrary to earlier studied 1,12-dibromo analogue 1,12-Br2-1,2- C2B10H10, its crystal packing is governed by the presence of the intermolecular I⋯I dihalogen bonds between the iodine atom attached to the carbon atom (acceptor) and the iodine atom attached to the antipodal boron atom (donor) of the carborane cage. The observed dihalogen bonds belong to the II type and are characterized by classical parameters: shortened I⋯I distance of 3.5687(9) Å, C–I⋯I angle of 172.61(11)° and B–I⋯I angle of 92.98(12)°.

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.

Polycentric binding in complexes of trimethylamine-N-oxide with dihalogens

Dihalogens readily interact with trimethylamine-N-oxide under ambient conditions. Accordingly, herein, stable 1 : 1 adducts were obtained in the case of iodine chloride and iodine bromide. The crystal and molecular structure of the trimethylamine-N-oxide-iodine chloride adduct was solved. Furthermore, the geometry and electronic structure of the trimethylamine-N-oxide-dihalogen complexes were studied computationally. Only molecular ensembles were found in the global minimum for the 1 : 1 stoichiometry. The O⋯X-Y halogen bond is the main factor for the thermodynamic stability of these complexes. Arguments for electrostatic interactions as the driving force for this noncovalent interaction were discussed. Also, the equilibrium structures are additionally stabilised by weak C-H⋯X hydrogen bonds. Consequently, formally monodentate ligands are bound in a polycentric manner.

Understanding the impact of anticancer halogenated inhibitors and various functional groups (X = Cl, F, CF3, CH3, NH2, OH, H) of casein kinase 2 (CK2)

Main focus of study is to understand potency of halogen (X = Br) atom that exists in tetrabromobenzotriazole (TBB) derivatives of crystal CK2 ligand along with hinge region amino acids (VAL45, PHE113, GLU114, VAL116, ASN118) through interaction energy analysis. In turn to attain profound insight on nature of stabilization of core CK2 ligands: 1ZOE-L1, 1ZOG-L2, 1ZOH-L3, 2OXX-L4, 2OXY-L5, 3KXG-L6, 3KXH-L7 -L7 and 3KXM-L8, having four bromine atoms, we attempted to mutate all bromine (X = Br) atoms by various functional groups (X = Cl, F, CF₃, CH₃, NH₂, OH, H) and binding strength along with amino acids was calculated. Most stable ligands exist in mutated NH₂ functional groups: 1ZOG-L2, 1ZOH-L3, 2OXX-L4, 3KXM-L8 having interaction energy as -5.21, -14.87, -6.69 and -11.72 kcal/mol respectively, revealing strong binding strength. Second most stable mutated Cl functional group ligands also play a major role in 1ZOH-L3, 2OXX-L4 and 3KXM-L8 having interaction energy as -6.89, -5.37, and -10.48 kcal/mol respectively. Overall, this study will pave way for crystal growth and medicinal chemist to have cleared perceptive about structural properties of CK2 halogenated ligands with new insight on CK2 mutated functional group ligands. Further, it insists us to reuse existing CK2 crystal ligand with more preferable suggested binding contacts in course of new functional groups that lead to anticancer affinity.Communicated by Ramaswamy H. Sarma.

The HOX⋯SO2 (X=F, Cl, Br, I) Binary Complexes: Implications for Atmospheric Chemistry

Sulfur dioxide and hypohalous acids (HOX, X=F, Cl, Br, I) are ubiquitous molecules in the atmosphere that are central to important processes like seasonal ozone depletion, acid rain, and cloud nucleation. We present the first theoretical examination of the HOX⋯SO₂ binary complexes and the associated trends due to halogen substitution. Reliable geometries were optimized at the CCSD(T)/aug-cc-pV(T+d)Z level of theory for HOF and HOCl complexes. The HOBr and HOI complexes were optimized at the CCSD(T)/aug-cc-pV(D+d)Z level of theory with the exception of the Br and I atoms which were modeled with an aug-cc-pwCVDZ-PP pseudopotential. 27 HOX⋯SO₂ complexes were characterized and the focal point method was employed to produce CCSDT(Q)/CBS interaction energies. Natural Bond Orbital analysis and Symmetry Adapted Perturbation Theory were used to classify the nature of each principle interaction. The interaction energies of all HOX⋯SO₂ complexes in this study ranged from 1.35 to 3.81 kcal mol^-1 . The single-interaction hydrogen bonded complexes spanned a range of 2.62 to 3.07 kcal mol^-1 , while the single-interaction halogen bonded complexes were far more sensitive to halogen substitution ranging from 1.35 to 3.06 kcal mol^-1 , indicating that the two types of interactions are extremely competitive for heavier halogens. Our results provide insight into the interactions between HOX and SO₂ which may guide further research of related systems.

Structures and Chemical Bonding in Antimony(III) Bromide Complexes with Pyridine

Weakly or "partially" bonded molecules are an important link between the chemical and van der Waals interactions. Molecular structures of six new SbBr₃ -Py complexes in the solid state have been determined by single-crystal X-ray diffraction analysis. In all complexes all Sb atoms adopt a pseudo-octahedral coordination geometry which is completed by additional Sb⋅⋅⋅Br contacts shorter than the sum of the van der Waals radii, with Br-Sb⋅⋅⋅Br angles close to 180°. To reveal the nature of Sb-Br and Sb-N interactions, the DFT calculations were performed followed by the analysis of the electrostatic potentials, the orbital interactions and the topological analysis. Based on Natural Bond Orbital (NBO) analysis, the Sb-Br interactions range from the covalent bonds to the pnictogen bonds. A simple structural parameter, non-covalence criterion (NCC) is defined as a ratio of the atom-atom distance to the linear combination of sums of covalent and van der Waals radii. NCC correlates with E(2) values for Sb-N, Sb-Cl and Sb-Br bonds, and appears to be useful criterion for a preliminary evaluation of the bonding situation.

One-Dimensional Diiodine-Iodobismuthate(III) Hybrids Cat3{[Bi2I9](I2)3}: Syntheses, Stability, and Optical Properties

One-dimensional iodine-rich iodobismuthates(III), Cat₃{[Bi₂I₉](I₂)₃} [Cat = 1,4-MePy (1) and 1-EtBMAP (2)], feature the highest amount of "trapped" diiodine units in polyhalogen-halometalates of p-block elements. Both complexes have narrow optical band gaps (1.55 and 1.63 eV, respectively) and moderate thermal stability.

Intermolecular Interactions Involving Heavy Alkenes H2Si=TH2 (T = C, Si, Ge, Sn, Pb) with H2O and HCl: Tetrel Bond and Hydrogen Bond

The intermolecular interactions between the heavy alkenes H₂Si=TH₂ (T = C, Si, Ge, Sn, Pb) and H₂O or HCl have been explored at the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVDZ level. The various hydrogen bond (HB) and tetrel bond (TB) complexes can be located on the basis of molecular electrostatic potential maps of the isolated monomers. The competition between TB and HB interactions has been investigated through the relaxed potential energy surface scan. The results indicate that the HB complexes become more and more unstable relative to the TB complexes with the increase of the T atomic number, and cannot even retain as a minimum in some cases, for H₂Si=TH2···H₂O systems. In contrast, the HB complexes are generally more stable than TB complexes, and the TB complexes exhibit rather weak binding strength, for H₂Si=TH₂···HCl systems. The majority of the TB complexes formed between H₂Si=TH₂ and H₂O possesses very strong binding strength with covalent characteristics. The noncovalent TB complexes can be divided into two types on the basis of the orbital interactions: π-hole complexes, with binding angles ranging from 91 to 111°, and hybrid σ/π-hole complexes, with binding angles ranging from 130 to 165°. The interplay between different molecular interactions has been explored, and an interesting result is that the covalent TB interaction is significantly abated and becomes noncovalent because of the competitive effect.

Phase Stability of Chloroform and Dichloromethane at High Pressure

Chloroform (CHCl3) and dichloromethane (CH2Cl2) are model systems for the study of intermolecular interactions, such as hydrogen bonds and halogen–halogen interactions. Here we report a joint computational (density-functional perturbation theory (DFPT) modelling) and experimental (Raman scattering) study on the behaviour of the crystals of these compounds up to a pressure of 32 GPa. Comparing the experimental information on the Raman band positions and intensities with the results of calculations enabled us to characterize the pressure-induced evolution of the crystal structure of both compounds. We find that the previously proposed P63 phase of CHCl3 is in fact a metastable structure, and that up to 32 GPa the ambient-pressure Pnma structure is the ground state polymorph of this compound. For CH2Cl2 we confirm the stability of the ambient-pressure Pbcn structure up to 32 GPa. We show that the high-pressure evolution of the crystal geometry of CHCl3 in the Pnma structure is a result of the subtle balance between dipole–dipole interactions, hydrogen bonds and Cl···Cl contacts. For CH2Cl2 (Pbcn structure) the dipole–dipole interactions and hydrogen bonds are the main factors influencing the pressure-induced changes in the geometry.

On single-electron magnesium bonding formation and the effect of methyl substitution

The complexes formed between MgX₂ (X = F, H) molecules and alkyl radicals Y [Y = CH₃, CH₂CH₃, CH(CH₃)₂, and C(CH₃)₃] have been characterized by using quantum chemical methods. The binding distance in all cases is less than the sum of vdW radii of Mg and C, indicating the formation of a non-covalent interaction, namely single-electron magnesium bond. Energy decomposition analysis reveals that electrostatic and polarization contributions are the major components responsible for the stability of the studied complexes. According to interaction energy, atoms in molecules, and independent gradient model analyses, methyl substitution on electron donor Y imposes a positive effect on its complexation with MgX₂. When compared with other nonbonded interactions, the single-electron magnesium bond is found to have strength comparable to those of the single-electron beryllium bond and π-magnesium bond.

The complexes formed between MgX₂ (X = F, H) molecules and alkyl radicals Y [Y = CH₃, CH₂CH₃, CH(CH₃)₂, and C(CH₃)₃] have been characterized by using quantum chemical methods.

Rising trend on the halogen and non-halogen derivatives (Br, Cl, CF3, F, CH3 and NH2) in ruminal β-d-Xylopyranose - a quantum chemical perspective

The utmost aim of the current study is to find significance of the binding affinity in the halogen and non-halogen derivatives: Br, Cl, CF₃, F, CH₃ and NH₂ of β-d-Xylopyranose with the hinge region amino acids of ruminant-β-glycosidase. The interaction energy analysis was carried out in detail through various density functional studies as M062X/def2-QZVP, M062X/LANL2DZ, B3LYP/LANL2DZ and M06HF/LANL2DZ level of theories. The total interaction energy of halogen derivatives: Br, Cl, F and CF₃ are -618.21, -599.00, -720.45 and -553.08 kcal/mol respectively, and non-halogen derivative: amine group (NH₂) is -87.96 kcal/mol at M062X/def2-QZVP level of theory, which exist with strong binding affinity. Ligand properties: dipole moment, polarizability, volume, molecular mass, electrostatic potential map was evaluated to understand its electrostatic and structural behavior. The nature of the bonds was inferred from the electrostatic potential map for all the halogen and non-halogen derivatives ligand. The stabilization energy from NBO analysis reveals the stability of single hydrogen and halogen bonds (N-H…Br, C-Br…O, N-H…Cl, C-Cl…O, O-H…F, C-H…F, N-H…F, C-F…O, N-H…O, O-H…O, N-H…N, O-H…N) in β-d-Xylopyranose and its derivatives. Overall, this study paves way for scientist and medicinal chemist in modelling new drugs. Further, it suggests mutations that increase the binding and may enhance the catalytic action and strengthen the complex diet in animals and hence recommended for experimental synthesis.Communicated by Ramaswamy H. Sarma.

Unconventional Type III Halogen···Halogen Interactions: A Quantum Mechanical Elucidation of σ-Hole···σ-Hole and Di-σ-Hole Interactions

Herein, two unconventional type III halogen···halogen interactions, namely, σ-hole···σ-hole and di-σ-hole interactions, were reported in a series of halogenated complexes. In type III, the A-halogen···halogen angles are typically equal to 180°, and the occurrence of σ-hole on halogen atoms is mandatory. Using diverse quantum mechanical calculations, it was demonstrated that the occurrence of such interactions with binding energies varied from -0.35 to -1.30 kcal/mol. Symmetry-adapted perturbation theory-based energy decomposition analysis (SAPT-EDA) revealed that type III interactions are dominated by dispersion forces, while electrostatic forces are unfavorable. Cambridge Structure Database (CSD) survey unveiled the experimental evidence for the manifestation of σ-hole···σ-hole interactions in crystal structures. This work might be deemed as a foundation for a vast number of forthcoming crystal engineering and materials science studies.

Intramolecular Nicholas Reactions in the Synthesis of Heteroenediynes Fused to Indole, Triazole, and Isocoumarin

The applicability of an intramolecular Nicholas reaction for the preparation of 10-membered O- and N-enediynes fused to indole, 1,2,3-triazole, and isocoumarin was investigated. The general approach to acyclic enediyne precursors fused to heterocycles includes inter- and intramolecular buta-1,3-diyne cyclizations with the formation of iodoethynylheterocycles, followed by Sonogashira coupling. The nature of both a heterocycle and a nucleophilic group affects the possibility of a 10-membered ring closure by the Nicholas reaction. Among oxacycles, an isocoumarin-fused enediyne was obtained. In the case of O-enediyne annulated with indole, instead of the formation of a 10-membered cycle, BF₃-promoted addition of an OH-group to the proximal triple bond at the C3 position afforded dihydrofuryl-substituted indole. For 1,2,3-triazole-fused analogues, using NH-Ts as a nucleophilic functional group allowed obtaining 10-membered azaenediyne, while the substrate with a hydroxyl group gave only traces of the desired 10-membered oxacycle. An improved method for the deprotection of Co-complexes of cyclic enediynes using tetrabutylammonium fluoride in an acetone/water mixture and the investigation of the 10-membered enediynes' reactivity in the Bergman cyclization are also reported. In the solid state, all synthesized iodoethynylheterocycles were found to be involved in halogen bond (XB) formation with either O or N atoms as XB acceptors.

Combined Molecular Dynamics and DFT Simulation Study of the Molecular and Polymer Properties of a Catechol-Based Cyclic Oligomer of Polyether Ether Ketone

The geometrical, energetic, noncovalent, and material properties of a catechol-based cyclic oligomer of Polyether Ether Ketone (PEEK) called o-PEEK were investigated using Molecular Dynamics (MD) and Density Functional Theory (DFT) simulations. The DFT (and MD) calculation performed with the PBEsol functional (and COMPASS II force field) gave a density of 1.39 (and 1.36) gcm⁻³ and a volume of 2744.5 (and 2808.5) cm³ for o-PEEK and are comparable with the corresponding experimental values of 1.328 gcm⁻³ and 2884.6 cm³, respectively. The absolute values of the glass transition temperature (T_g) MD simulated using the unit-cell and 2 × 2 × 2 supercell geometries of the o-PEEK system were 424.4 and 428.6 K, respectively. Although these values slightly differ from each other, both are close to the experiment (T_g = 418.2 K). The results of the (charge) density gradient analysis suggest that the supramolecular assembly between the o-PEEK oligomers in the experimentally observed infinite semi-crystal is driven by a wide range of noncovalent interactions. While the individual local interactions between the oligomers were recognized to be weak-to-medium in strength and are theoretically difficult to quantify, the B97-D3/cc-pVTZ level stabilization energy responsible for the formation of each of the five binary complex configurations extracted from the PBEsol relaxed 2 × 2 × 2 supercell geometry of the o-PEEK system was calculated to vary between –3.5 and –33.0 kcal mol⁻¹.

Crystal Structure and Supramolecular Architecture of Inorganic Ligand-Coordinated Salen-Type Schiff Base Complex: Insights into Halogen Bond from Theoretical Analysis and 3D Energy Framework Calculations

To identify the effects of halogen bonding in the architecture of Schiff base complex supramolecular networks, we introduced halogenated Schiff-base 3-Br-5-Cl-salen as ligand and isolated a new salen-type manganese(III) complex [MnIII(Cl)(H2O)(3-Br-5-Cl-salen)] (1) where 3-Br-5-Cl-salen = N,N’-bis(3-bromo-5-chlorosalicylidene)-1,2-diamine. The complex was investigated in the solid-state for halogen bonds (XBs) by single crystal X-ray structure analysis. Meanwhile, theoretical calculations were carried out to rationalize the formation mechanism of different types of XBs in the complex. The analysis result of electronic structure of the halogen bonds indicated that the chlorine atom coordinated to the Mn(III) center possesses the most negative potential and acts as anionic XB acceptor (electron donor) to the adjacent substituted halogens on the ligand, meanwhile the intermolecular Mn-Cl···X-C halogen bonding plays a significant role in directing the packing arrangement of adjacent molecules and linking the 2D layers into a 3D network. In order to verify the above results and acquire detailed information, the title complex was further analyzed by using the Hirshfeld surface analyses.

Halogen Bonding Interactions of Polychlorinated Biphenyls and the Potential for Thyroid Disruption

Polychlorinated biphenyl (PCB) flame retardants are persistent pollutants and inhibit neurodevelopment, particularly in the early stages of life. Halogen bonding (XB) to the iodothyronine deiodinases (Dio) that modulate thyroid hormones (THs) is a potential mechanism for endocrine disruption. Cl⋅⋅⋅Se XB interactions of PCBs with SeMe- , a small model of the Dio active site selenocysteine, are compared with previous results on polybrominated diphenylethers (PBDEs) and THs using density functional theory. PCBs generally display weaker XB interactions compared to PBDEs and THs, consistent with the dependence of XB strength on the size of the halogen (I>Br>Cl). PCBs also do not meet a proposed energy threshold for substrates to undergo dehalogenation, suggesting they may behave as competitive inhibitors of Dio in addition to other mechanisms of endocrine disruption. XB interactions in PCBs are position-dependent, with ortho interactions slightly more favorable than meta and para interactions, suggesting that PCBs may have a greater effect on certain classes of Dio. Flexibility of PCBs around the biphenyl C-C bond is limited by ortho substitutions relative to the biphenyl linkage, which may contribute to the ability to inhibit Dio and other TH-related proteins.

In Situ Assessment of Intrinsic Strength of X-I⋯OA-Type Halogen Bonds in Molecular Crystals with Periodic Local Vibrational Mode Theory

Periodic local vibrational modes were calculated with the rev-vdW-DF2 density functional to quantify the intrinsic strength of the X-I⋯OA-type halogen bonding (X = I or Cl; OA: carbonyl, ether and N-oxide groups) in 32 model systems originating from 20 molecular crystals. We found that the halogen bonding between the donor dihalogen X-I and the wide collection of acceptor molecules OA features considerable variations of the local stretching force constants (0.1–0.8 mdyn/Å) for I⋯O halogen bonds, demonstrating its powerful tunability in bond strength. Strong correlations between bond length and local stretching force constant were observed in crystals for both the donor X-I bonds and I⋯O halogen bonds, extending for the first time the generalized Badger’s rule to crystals. It is demonstrated that the halogen atom X controlling the electrostatic attraction between the σ-hole on atom I and the acceptor atom O dominates the intrinsic strength of I⋯O halogen bonds. Different oxygen-containing acceptor molecules OA and even subtle changes induced by substituents can tweak the n→σ∗(X-I) charge transfer character, which is the second important factor determining the I⋯O bond strength. In addition, the presence of the second halogen bond with atom X of the donor X-I bond in crystals can substantially weaken the target I⋯O halogen bond. In summary, this study performing the in situ measurement of halogen bonding strength in crystalline structures demonstrates the vast potential of the periodic local vibrational mode theory for characterizing and understanding non-covalent interactions in materials.

Does Chlorine in CH3Cl Behave as a Genuine Halogen Bond Donor?

The CH3Cl molecule has been used in several studies as an example purportedly to demonstrate that while Cl is weakly negative, a positive potential can be induced on its axial surface by the electric field of a reasonably strong Lewis base (such as O=CH2). The induced positive potential then has the ability to attract the negative site of the Lewis base, thus explaining the importance of polarization leading to the formation of the H3C–Cl···O=CH2 complex. By examining the nature of the chlorine’s surface in CH3Cl using the molecular electrostatic surface potential (MESP) approach, with MP2/aug-cc-pVTZ, we show that this view is not correct. The results of our calculations demonstrate that the local potential associated with the axial surface of the Cl atom is inherently positive. Therefore, it should be able to inherently act as a halogen bond donor. This is shown to be the case by examining several halogen-bonded complexes of CH3Cl with a series of negative sites. In addition, it is also shown that the lateral portions of Cl in CH3Cl features a belt of negative electrostatic potential that can participate in forming halogen-, chalcogen-, and hydrogen-bonded interactions. The results of the theoretical models used, viz. the quantum theory of atoms in molecules; the reduced density gradient noncovalent index; the natural bond orbital analysis; and the symmetry adapted perturbation theory show that Cl-centered intermolecular bonding interactions revealed in a series of 18 binary complexes do not involve a polarization-induced potential on the Cl atom.

Spin-orbit coupling is the key to unraveling intriguing features of the halogen bond involving astatine

The effect of spin-orbit coupling (SOC) on the halogen bond involving astatine has been investigated using state-of-the-art two- and four-component relativistic calculations. Adducts between Cl-X (X = Cl, Br, I and At) and ammonia have been selected to establish a trend on going down the periodic table. The SOC influence has been explored not only on the geometric and energetic features that can be used to characterize the halogen bond strength but also on the three main contributions to it that are the charge transfer, the "σ-hole" (i.e. the localized region with a net positive electrostatic potential at the halogen site) and the "polar flattening" (which is related to the effective shape of the halogen site). A surprisingly large increase of the Cl-At dipole moment, due to the inclusion of SOC, has been worked out using four-component CCSD(T) reference calculations, indicating that this bond is significantly more ionic than one may predict. Due to the SOC effect, which induces a peculiar charge accumulation on the At side in the Cl-At dimer, a weakening of the astatine-mediated halogen bond occurs arising from the (i) reduced amount of charge transfer, (ii) decrease of the polar flattening and (iii) lowering of the short-range Coulomb potential. The analysis of the electronic structure of the Cl-At moiety allows for a rationalization of the SOC effects on all the considered features of the halogen bond, including an unprecedented unsymmetrical charge back-donation from Cl-At to ammonia.

Charge Displacement Analysis-A Tool to Theoretically Characterize the Charge Transfer Contribution of Halogen Bonds

Theoretical bonding analysis is of prime importance for the deep understanding of the various chemical interactions, covalent or not. Among the various methods that have been developed in the last decades, the analysis of the Charge Displacement function (CD) demonstrated to be useful to reveal the charge transfer effects in many contexts, from weak hydrogen bonds, to the characterization of σ hole interactions, as halogen, chalcogen and pnictogen bonding or even in the decomposition of the metal-ligand bond. Quite often, the CD analysis has also been coupled with experimental techniques, in order to give a complete description of the system under study. In this review, we focus on the use of CD analysis on halogen bonded systems, describing the most relevant literature examples about gas phase and condensed phase systems. Chemical insights will be drawn about the nature of halogen bond, its cooperativity and its influence on metal-ligand bond components.

Predicting reactive sites with quantum chemical topology: carbonyl additions in multicomponent reactions

The reactivity of an atom within a molecule depends mostly on the way the electron density polarizes reflected in the quadrupole moment of the reactive atom.

Identifying intermolecular atom⋯atom interactions that are not just bonding but also competitive

This highlight criticises the QTAIM method and discusses algorithms for identifying intermolecular interactions that are both bonding and competitive.

The role of non-covalent intermolecular interactions on the diversity of crystal packing in supramolecular dihalopyridine–silver(i) nitrate complexes

The crystal structures of six novel Ag⁺ complexes with NO₃⁻ and dihalopyridines revealed intriguing differences that were interpreted by DFT calculations.

Tetrachloroplatinate(ii) anion as a square-planar tecton for crystal engineering involving halogen bonding

The tetrachloroplatinate(ii) anion behaves as a useful XB-accepting tecton toward sigma-hole-donating organohalide species.

Metal-Halogen Bonding Seen through the Eyes of Vibrational Spectroscopy

Incorporation of a metal center into halogen-bonded materials can efficiently fine-tune the strength of the halogen bonds and introduce new electronic functionalities. The metal atom can adopt two possible roles: serving as halogen acceptor or polarizing the halogen donor and acceptor groups. We investigated both scenarios for 23 metal–halogen dimers trans-M(Y₂)(NC₅H₄X-3)₂ with M = Pd(II), Pt(II); Y = F, Cl, Br; X = Cl, Br, I; and NC₅H₄X-3 = 3-halopyridine. As a new tool for the quantitative assessment of metal–halogen bonding, we introduced our local vibrational mode analysis, complemented by energy and electron density analyses and electrostatic potential studies at the density functional theory (DFT) and coupled-cluster single, double, and perturbative triple excitations (CCSD(T)) levels of theory. We could for the first time quantify the various attractive contacts and their contribution to the dimer stability and clarify the special role of halogen bonding in these systems. The largest contribution to the stability of the dimers is either due to halogen bonding or nonspecific interactions. Hydrogen bonding plays only a secondary role. The metal can only act as halogen acceptor when the monomer adopts a (quasi-)planar geometry. The best strategy to accomplish this is to substitute the halo-pyridine ring with a halo-diazole ring, which considerably strengthens halogen bonding. Our findings based on the local mode analysis provide a solid platform for fine-tuning of existing and for design of new metal–halogen-bonded materials.

Chiral Chalcogen Bond Donors Based on the 4,4'-Bipyridine Scaffold

Organocatalysis through chalcogen bonding (ChB) is in its infancy, as its proof-of-principle was only reported in 2016. Herein, we report the design and synthesis of new chiral ChB donors, as well as the catalytic activity evaluation of the 5,5′-dibromo-2,2′-dichloro-3-((perfluorophenyl)selanyl)-4,4′-bipyridine as organocatalyst. The latter is based on the use of two electron-withdrawing groups, a pentafluorophenyl ring and a tetrahalo-4,4′-bipyridine skeleton, as substituents at the selenium center. Atropisomery of the tetrahalo-4,4′-bipyridine motif provides a chiral environment to these new ChB donors. Their synthesis was achieved through either selective lithium exchange and trapping or a site-selective copper-mediated reaction. Pure enantiomers of the 3-selanyl-4,4′-bipyridine were obtained by high performance liquid chromatography enantioseparation on specific chiral stationary phase, and their absolute configuration was assigned by comparison of the measured and calculated electronic circular dichroism spectra. The capability of the selenium compound to participate in σ-hole-based interactions in solution was studied by ¹⁹F NMR. Even if no asymmetric induction has been observed so far, the new selenium motif proved to be catalytically active in the reduction of 2-phenylquinoline by Hantzsch ester.

Does Oxygen Feature Chalcogen Bonding?

Using the second-order Møller-Plesset perturbation theory (MP2), together with Dunning's all-electron correlation consistent basis set aug-cc-pVTZ, we show that the covalently bound oxygen atom present in a series of 21 prototypical monomer molecules examined does conceive a positive (or a negative) σ-hole. A σ-hole, in general, is an electron density-deficient region on a bound atom M along the outer extension of the R-M covalent bond, where R is the reminder part of the molecule, and M is the main group atom covalently bonded to R. We have also examined some exemplar 1:1 binary complexes that are formed between five randomly chosen monomers of the above series and the nitrogen- and oxygen-containing Lewis bases in N2, PN, NH3, and OH2. We show that the O-centered positive σ-hole in the selected monomers has the ability to form the chalcogen bonding interaction, and this is when the σ-hole on O is placed in the close proximity of the negative site in the partner molecule. Although the interaction energy and the various other 12 characteristics revealed from this study indicate the presence of any weakly bound interaction between the monomers in the six complexes, our result is strongly inconsistent with the general view that oxygen does not form a chalcogen-bonded interaction.

Nature of halogen-centered intermolecular interactions in crystal growth and design: Fluorine-centered interactions in dimers in crystalline hexafluoropropylene as a prototype

The wide occurrence of halogen-centered noncovalent interactions in crystal growth and design prompted this study, which includes a mini review of recent advances in the field. Particular emphasis is placed on providing compelling theoretical evidence of the formation of these interactions between sites of positive electrostatic potential, as well as between sites of negative electrostatic potential, localized on the electrostatic surfaces of the bound fluorine atoms in a prototypical system, hexafluoropropylene (C₃ F₆ ), upon its interaction with another same molecule to form (C₃ F₆ )₂ dimers. The existence of σ- and π-hole interactions is shown for the stable dimers. Even so, weakly bound interactions locally responsible in holding the molecular fragments together cannot and should not be overlooked since they are partly responsible for determining the overall geometry of the crystal. The results of combined quantum theory of atoms in molecules, molecular electrostatic surface potential, and reduced density gradient noncovalent interaction analyses showed that these latter interactions do indeed play a role in the stability and growth of crystalline C₃ F₆ itself and the (C₃ F₆ )₂ dimers. A symmetry adapted perturbation theory energy decomposition analysis leads to the conclusion that a great majority of the (C₃ F₆ )₂ dimers examined are the consequence of dispersion (and electrostatics), with nonnegligible contribution from polarization, which together competes with an exchange repulsion component to determine the equilibrium geometries. In a few structures of the (C₃ F₆ )₂ dimer, the fluorine is found to serve as a six-center five-bond donor/acceptor, as found for carbon in other systems (Malischewski and Seppelt, Angew. Chem. Int. Ed. 2017, 56, 368). © 2019 Wiley Periodicals, Inc.

The chalcogen bond: can it be formed by oxygen?

This study theoretically investigates the possibility of oxygen-centered chalcogen bonding in several complexes. Shown in the graph is such a bonding scenario formed between the electrophile on O in OF₂ and the nucleophile on O in H₂CO.