Halogen and Chalcogen Bond Energies Evaluated Using Electron Density Properties

A novel approach for estimating the strength of argentophilic and aurophilic interactions using QTAIM parameters

Metallophilic interactions, specifically argentophilic (Ag⋯Ag) and aurophilic (Au⋯Au) interactions, play a crucial role in stabilizing various molecular and solid-state structures. In this manuscript, we present a convenient method to estimate the strength of argentophilic and aurophilic interactions based on quantum theory of atoms in molecules (QTAIM) parameters evaluated at the bond critical points connecting the metal centres. We employ density functional theory (DFT) calculations and the QTAIM parameters to develop this energy predictor. To validate the reliability and applicability of our method, we test it using a selection of X-ray crystal structures extracted from the cambridge structural database (CSD), where argentophilic and aurophilic interactions are known to be significant in their solid-state arrangements. This method offers a distinct advantage in systems where multiple interactions, beyond metallophilic interactions, contribute to the overall stability of the structure. By employing our approach, researchers can distinctly quantify the strength of argentophilic and aurophilic interactions, facilitating a deeper understanding of their impact on molecular and solid-state properties. This method fills a critical gap in the existing literature, offering a valuable tool to researchers seeking to unravel the intricate interactions in metal-containing compounds.

On the Existence of Pnictogen Bonding Interactions in As(III) S-Adenosylmethionine Methyltransferase Enzymes

As(III) S-adenosylmethionine methyltransferases, pivotal enzymes in arsenic metabolism, facilitate the methylation of arsenic up to three times. This process predominantly yields trivalent mono- and dimethylarsenite, with trimethylarsine forming in smaller amounts. While this enzyme acts as a detoxifier in microbial systems by altering As(III), in humans, it paradoxically generates more toxic and potentially carcinogenic methylated arsenic species. The strong affinity of As(III) for cysteine residues, forming As(III)-thiolate bonds, is exploited in medical treatments, notably in arsenic trioxide (Trisenox®), an FDA-approved drug for leukemia. The effectiveness of this drug is partly due to its interaction with cysteine residues, leading to the breakdown of key oncogenic fusion proteins. In this study, we extend the understanding of As(III)'s binding mechanisms, showing that, in addition to As(III)-S covalent bonds, noncovalent O⋅⋅⋅As pnictogen bonding plays a vital role. This interaction significantly contributes to the structural stability of the As(III) complexes. Our crystallographic analysis using the PDB database of As(III) S-adenosylmethionine methyltransferases, augmented by comprehensive theoretical studies including molecular electrostatic potential (MEP), quantum theory of atoms in molecules (QTAIM), and natural bond orbital (NBO) analysis, emphasizes the critical role of pnictogen bonding in these systems. We also undertake a detailed evaluation of the energy characteristics of these pnictogen bonds using various theoretical models. To our knowledge, this is the first time pnictogen bonds in As(III) derivatives have been reported in biological systems, marking a significant advancement in our understanding of arsenic's molecular interactions.

Selective recognition and extraction of iodide from pure water by a tripodal selenoimidazol(ium)-based chalcogen bonding receptor

A selenium-based tripodal chalcogen bond (ChB) donor TPI-3Se is demonstrated for the recognition and extraction of I- from 100% water medium. NMR and ITC studies with the halides reveal that the ChB donor selectively binds with the large, weakly hydrated I-. Interestingly, I- crystallizes out selectively in the presence of other halides supporting the superiority of the selective recognition of I-. The X-ray structure of the ChB-iodide complex manifests both the μ1 and μ2 coordinated interactions, which is rare in the C-Se···I chalcogen bonding. Furthermore, to validate the selective I- binding potency of TPI-3Se in pure water, comparisons are made with its hydrogen and halogen bond donor analogs. The computational analysis also provides the mode of I- recognition by TPI-3Se. Importantly, this receptor is capable of extracting I- from pure water through selenium sigma-hole and I- interaction with a high degree of efficiency (∼70%).

Chalcogen and Hydrogen Bond Team up in Driving Anion⋅⋅⋅Anion Self-Assembly

H-selenite anions (HSeO3 - ) form in the solid unprecedented anionic supramolecular chains wherein single units are assembled via alternating short Se⋅⋅⋅O and H⋅⋅⋅O contacts. Crystallographic analyses and computational studies (the quantum theory of "atoms-in-molecules", QTAIM, and the noncovalent interaction plot, NCIPlot) consistently prove the attractive nature of these chalcogen bonds (ChBs) and hydrogen honds (HBs), the Janus-type character of HSeO3 - anions which act as both donors and acceptors of ChB and HB, and the possible stability of anion dimers in solution. The effectiveness of the ChBs herein described may lead to consider the HSeO3 - moiety as a new entry in the toolbox of crystal engineering based on ChB.

An Experimental and Theoretical Insight into I2 /Br2 Oxidation of Bis(pyridin-2-yl)Diselane and Ditellane

The reactivity between bis(pyridin-2-yl)diselane o Py2 Se2 and ditellane o Py2 Te2 (L1 and L2, respectively; o Py=pyridyn-2-yl) and I2 /Br2 is discussed. Single-crystal structure analysis revealed that the reaction of L1 with I2 yielded [(HL1+ )(I- )⋅5/2I2 ]∞ (1) in which monoprotonated cations HL1+ template a self-assembled infinite pseudo-cubic polyiodide 3D-network, while the reaction with Br2 yielded the dibromide Ho PySeII Br2 (2). The oxidation of L2 with I2 and Br2 yielded the compounds Ho PyTeII I2 (3) and Ho PyTeIV Br4 (6), respectively, whose structures were elucidated by X-ray diffraction analysis. FT-Raman spectroscopy measurements are consistent with a 3c-4e description of all the X-Ch-X three-body systems (Ch=Se, Te; X=Br, I) in compounds 2, 3, Ho PyTeII Br2 (5), and 6. The structural and spectroscopic observations are supported by extensive theoretical calculations carried out at the DFT level that were employed to study the electronic structure of the investigated compounds, the thermodynamic aspects of their formation, and the role of noncovalent σ-hole halogen and chalcogen bonds in the X⋅⋅⋅X, X⋅⋅⋅Ch and Ch⋅⋅⋅Ch interactions evidenced structurally.

Computational Study of Driving Forces in ATSP, PDIQ, and P53 Peptide Binding: C═O···C═O Tetrel Bonding Interactions at Work

Understanding the molecular interactions that drive peptide folding is crucial to chemistry and biology. In this study, we analyzed the role of CO···CO tetrel bonding (TtB) interactions in the folding mechanism of three different peptides (ATSP, pDIQ, and p53), which exhibit a different propensity to fold in an α helix motif. To achieve this goal, we used both a recently developed Bayesian inference approach (MELDxMD) and Quantum Mechanics (QM) calculations at the RI-MP2/def2-TZVP level of theory. These techniques allowed us to study the folding process and to evaluate the strength of the CO···CO TtBs as well as the synergies between TtBs and hydrogen-bonding (HB) interactions. We believe that the results derived from our study will be helpful for those scientists working in computational biology, peptide chemistry, and structural biology.

σ-Hole interactions in organometallic catalysts: the case of methyltrioxorhenium(VII)

Methyltrioxorhenium(VII) (MTO) is a widely employed catalyst for metathesis, olefination, and most importantly, oxidation reactions. It is often preferred to other oxometal complexes due to its stability in air and higher efficiency. The seminal papers of K. B. Sharpless showed that when pyridine derivatives are used as co-catalysts, MTO-catalyzed olefin epoxidation with H₂O₂ as oxidant, a particularly useful reaction, is accelerated, with pyridine speeding up catalytic turnover and increasing the lifetime of MTO under the reaction conditions. In this paper, combined experimental and theoretical results show that the occurrence of σ-hole interactions in catalytic systems extends to MTO. Four crystalline adducts between MTO and aliphatic and heteroaromatic bases are obtained, and their X-ray analyses display short Re⋯N/O contacts opposite to both O-Re and C-Re covalent bonds with geometries consistent with σ-hole interactions. Computational analyses support the attractive nature of these close contacts and confirm that their features are typical of σ-hole interactions. The understanding of the nature of Re⋯N/O interactions may help to optimize the ligand-acceleration effect of pyridine in the epoxidation of olefins under MTO catalysis.

Chalcogen Bond as a Factor Stabilizing Ligand Conformation in the Binding Pocket of Carbonic Anhydrase IX Receptor Mimic

It is postulated that the overexpression of Carbonic Anhydrase isozyme IX in some cancers contributes to the acidification of the extracellular matrix. It was proved that this promotes the growth and metastasis of the tumor. These observations have made Carbonic Anhydrase IX an attractive drug target. In the light of the findings and importance of the glycoprotein in the cancer treatment, we have employed quantum-chemical approaches to study non-covalent interactions in the binding pocket. As a ligand, the acetazolamide (AZM) molecule was chosen, being known as a potential inhibitor exhibiting anticancer properties. First-Principles Molecular Dynamics was performed to study the chalcogen and other non-covalent interactions in the AZM ligand and its complexes with amino acids forming the binding site. Based on Density Functional Theory (DFT) and post-Hartree-Fock methods, the metric and electronic structure parameters were described. The Non-Covalent Interaction (NCI) index and Atoms in Molecules (AIM) methods were applied for qualitative/quantitative analyses of the non-covalent interactions. Finally, the AZM-binding pocket interaction energy decomposition was carried out. Chalcogen bonding in the AZM molecule is an important factor stabilizing the preferred conformation. Free energy mapping via metadynamics and Path Integral molecular dynamics confirmed the significance of the chalcogen bond in structuring the conformational flexibility of the systems. The developed models are useful in the design of new inhibitors with desired pharmacological properties.

Description of non-covalent interactions in benzyl chalcocyanate crystals from smoothed Cromer-Mann electron density distribution functions

A well-known method to characterize non-covalent interactions consists in the topological analysis of electron density distribution (EDD) functions, complemented by the search for minima in the reduced density gradient (RDG) distributions. Here, we characterize intermolecular interactions occurring in crystals of benzyl chalcocyanate compounds through bond critical points (BCP) of the promolecular electron density (ED) built from the crystallographic Cromer-Mann parameters, at several smoothing levelst. The trajectories formed by thet-dependent BCP locations are interpreted in terms of the intermolecular interactions occurring within the crystal arrangements. Chalcogen…nitro BCPs are clearly present in the unsmoothed EDDs but are annihilated astincreases, while chalcogen…chalcogen BCPs appear and are among the only BCPs left at the highest smoothing level. The chalcogen bonds are differentiated from the other chalcogen interactions through the linear chalcogen…BCP…nitro geometry at low smoothing level and their more negative Laplacian values. The annihilation of CPs can be followed by the apparition of a RDG minimum, associated with a very weak interaction. Along the BCP trajectories, the Laplacian shows a progressive concentration of the ED in the intermolecular space within the crystals and adopts the most negative values at the shortest atom…atom separations. At the termination point of a BCP trajectory, the drastic increase of the ellipticity value illustrates the flattening of the EDD.

Self-Assembly of Supramolecular Architectures Driven by σ-Hole Interactions: A Halogen-Bonded 2D Network Based on a Diiminedibromido Gold(III) Complex and Tribromide Building Blocks

The reaction of the complex [Au(phen)Br₂](PF₆) (phen = 1,10-phenanthroline) with molecular dibromine afforded {[Au(phen)Br₂](Br₃)}_∞ (1). Single crystal diffraction analysis showed that the [Au(phen)Br₂]⁺ complex cations were bridged by asymmetric tribromide anions to form infinite zig-zag chains featuring the motif ···Au–Br···Br–Br–Br···Au–Br···Br–Br–Br···. The complex cation played an unprecedented halogen bonding (XB) donor role engaging type-I and type-II XB noncovalent interactions of comparable strength with symmetry related [Br₃]⁻ anions. A network of hydrogen bonds connects parallel chains in an infinite 2D network, contributing to the layered supramolecular architecture. DFT calculations allowed clarification of the nature of the XB interactions, showing the interplay between orbital mixing, analyzed at the NBO level, and electrostatic contribution, explored based on the molecular potential energy (MEP) maps of the interacting synthons.

Supramolecular Chalcogen‐Bonded Semiconducting Nanoribbons at Work in Lighting Devices

This work describes the design and synthesis of a π‐conjugated telluro[3,2‐β][1]‐tellurophene‐based synthon that, embodying pyridyl and haloaryl chalcogen‐bonding acceptors, self‐assembles into nanoribbons through chalcogen bonds. The ribbons π‐stack in a multi‐layered architecture both in single crystals and thin films. Theoretical studies of the electronic states of chalcogen‐bonded material showed the presence of a local charge density between Te and N atoms. OTFT‐based charge transport measurements showed hole‐transport properties for this material. Its integration as a p‐type semiconductor in multi‐layered Cu^I‐based light‐emitting electrochemical cells (LECs) led to a 10‐fold increase in stability (38 h vs. 3 h) compared to single‐layered devices. Finally, using the reference tellurotellurophene congener bearing a C−H group instead of the pyridyl N atom, a herringbone solid‐state assembly is formed without charge transport features, resulting in LECs with poor stabilities (<1 h).

Noncovalent interactions between benzochalcogenadiazoles and nitrogen bases

A theoretical study has been carried out on the intermolecular interactions between tetrafluoro-benzochalcogenadiazoles (chalcogen = S, Se, Te) and a series of nitrogen bases (FCN, ClCN, NP, trans-N₂H₂, pyridine, pyrazole, imidazole) at the B97-D3/def2-TZVP level, to obtain a better insight into the nature and strength of Ch···N chalcogen bond and secondary interaction in the binary and 1:2 ternary complexes. The dispersion force plays a prominent role on the stability of the sulfur complexes, and the electrostatic effect enhanced for the heavier chalcogen complexes. Most of intermolecular bonds display the characters of closed-shell and noncovalent interaction. For the complexes involving pyridine and imidazole, chalcogen bond is stronger than hydrogen bond, while the strength of chalcogen bond is equivalent to the secondary interaction for other complexes. With the addition of nitrogen base in the 1:2 complexes, chalcogen bond is weakened, while the secondary interaction remains unchanged. In the 1:2 complexes formed by pyridine and imidazole, stronger chalcogen bond results in larger negative cooperativity than that of other complexes.

Towards Anion Recognition and Precipitation with Water-Soluble 1,2,4-Selenodiazolium Salts: Combined Structural and Theoretical Study

The synthesis and structural characterization of a series of supramolecular complexes of bicyclic cationic pyridine-fused 1,2,4-selenodiazoles with various anions is reported. The binding of trifluoroacetate, tetrachloroaurate, tetraphenylborate, perrhenate, and pertechnetate anions in the solid state is regarded. All the anions interact with selenodiazolium cations exclusively via a pair of “chelating” Se⋯O and H⋯O non-covalent interactions, which make them an attractive, novel, non-classical supramolecular recognition unit or a synthon. Trifluoroacetate salts were conveniently generated via novel oxidation reaction of 2,2′-dipyridyl diselenide with bis(trifluoroacetoxy)iodo)benzene in the presence of corresponding nitriles. Isolation and structural characterization of transient 2-pyridylselenyl trifluoroacetate was achieved. X-ray analysis has demonstrated that the latter forms dimers in the solid state featuring very short and strong Se⋯O and Se⋯N ChB contacts. 1,2,4-Selenodiazolium trifluoroacetates or halides show good solubility in water. In contrast, (AuCl₄)⁻, (ReO₄)⁻, or (TcO₄)⁻ derivatives immediately precipitate from aqueous solutions. Structural features of these supramolecular complexes in the solid state are discussed. The nature and energies of the non-covalent interactions in novel assembles were studied by the theoretical methods. To the best of our knowledge, this is the first study that regards perrhenate and pertechnetate as acceptors in ChB interactions. The results presented here will be useful for further developments in anion recognition and precipitation involving cationic 1,2,4-selenodiazoles.

The Pnictogen Bond: The Covalently Bound Arsenic Atom in Molecular Entities in Crystals as a Pnictogen Bond Donor

In chemical systems, the arsenic-centered pnictogen bond, or simply the arsenic bond, occurs when there is evidence of a net attractive interaction between the electrophilic region associated with a covalently or coordinately bound arsenic atom in a molecular entity and a nucleophile in another or the same molecular entity. It is the third member of the family of pnictogen bonds formed by the third atom of the pnictogen family, Group 15 of the periodic table, and is an inter- or intramolecular noncovalent interaction. In this overview, we present several illustrative crystal structures deposited into the Cambridge Structure Database (CSD) and the Inorganic Chemistry Structural Database (ICSD) during the last and current centuries to demonstrate that the arsenic atom in molecular entities has a significant ability to act as an electrophilic agent to make an attractive engagement with nucleophiles when in close vicinity, thereby forming σ-hole or π-hole interactions, and hence driving (in part, at least) the overall stability of the system’s crystalline phase. This overview does not include results from theoretical simulations reported by others as none of them address the signatory details of As-centered pnictogen bonds. Rather, we aimed at highlighting the interaction modes of arsenic-centered σ- and π-holes in the rationale design of crystal lattices to demonstrate that such interactions are abundant in crystalline materials, but care has to be taken to identify them as is usually done with the much more widely known noncovalent interactions in chemical systems, halogen bonding and hydrogen bonding. We also demonstrate that As-centered pnictogen bonds are usually accompanied by other primary and secondary interactions, which reinforce their occurrence and strength in most of the crystal structures illustrated. A statistical analysis of structures deposited into the CSD was performed for each interaction type As···D (D = N, O, S, Se, Te, F, Cl, Br, I, arene’s π system), thus providing insight into the typical nature of As···D interaction distances and ∠R–As···D bond angles of these interactions in crystals, where R is the remainder of the molecular entity.

Metal Coordination Enhances Chalcogen Bonds: CSD Survey and Theoretical Calculations

In this study the ability of metal coordinated Chalcogen (Ch) atoms to undergo Chalcogen bonding (ChB) interactions has been evaluated at the PBE0-D3/def2-TZVP level of theory. An initial CSD (Cambridge Structural Database) inspection revealed the presence of square planar Pd/Pt coordination complexes where divalent Ch atoms (Se/Te) were used as ligands. Interestingly, the coordination to the metal center enhanced the σ-hole donor ability of the Ch atom, which participates in ChBs with neighboring units present in the X-ray crystal structure, therefore dictating the solid state architecture. The X-ray analyses were complemented with a computational study (PBE0-D3/def2-TZVP level of theory), which shed light into the strength and directionality of the ChBs studied herein. Owing to the new possibilities that metal coordination offers to enhance or modulate the σ-hole donor ability of Chs, we believe that the findings presented herein are of remarkable importance for supramolecular chemists as well as for those scientists working in the field of solid state chemistry.

Computational Evaluation of Me2 TCCP as Lewis Acid

Supramolecular adducts between dimethyl-2,2,3,3-tetracyanocyclopropane (Me2 TCCP) with 21 small (polar) molecules and 10 anions were computed with DFT (B3LYP-D3/def2-TZVP). Their optimized geometries were used to obtain interaction energies, and perform energy decomposition and 'atoms-in-molecules' analyses. A set of 38 other adducts were also evaluated for comparison purposes. Selected examples were further scrutinized by inspection of the molecular electrostatic potential maps, Noncovalent Interaction index plots, the Laplacian, the orbital interactions, and by estimating the Gibbs free energy of complexation in hexane solution. These calculations divulge the thermodynamic feasibility of Me2 TCCP adducts and show that complexation is typically driven by dispersion with less polarized partners, but by orbital interactions when more polarized or anionic guests are deployed. Most Me2 TCCP adducts are more stable than simple hydrogen bonding with water, but less stable than traditional Lewis adducts involving Me3 B, or a strong halogen bond such as with Br2 . Several bonding analyses showed that the locus of interaction is found near the electron poor sp3 -hydridized (NC)2 C-C(CN)2 carbon atoms. An empty hybrid σ*/π* orbital on Me2 TCCP was identified that can be held responsible for the stability of the most stable adducts due to donor-acceptor interactions.

Unveiling the role of pyrylium frameworks on π-stacking interactions: a combined ab initio and experimental study

A multidisciplinary study is presented to shed light on how pyrylium frameworks, as π-hole donors, establish π-π interactions. The combination of CSD analysis, computational modelling (ab intitio, DFT and MD simulations) and experimental NMR spectroscopy data provides essential information on the key parameters that characterize these intereactions, opening new avenues for further applications of this versatile heterocycle.

Charge Assisted S/Se Chalcogen Bonds in SAM Riboswitches: A Combined PDB and ab Initio Study

In this study, we provide experimental (Protein Data Bank (PDB) inspection) and theoretical (RI-MP2/def2-TZVP level of theory) evidence of the involvement of charge assisted chalcogen bonding (ChB) interactions in the recognition and folding mechanisms of S-adenosylmethionine (SAM) riboswitches. Concretely, an initial PDB search revealed several examples where ChBs between S-adenosyl methionine (SAM)/adenosyl selenomethionine (EEM) molecules and uracil (U) bases belonging to RNA take place. While these interactions are usually described as a merely Coulombic attraction between the positively charged S/Se group and RNA, theoretical calculations indicated that the σ holes of S and Se are involved. Moreover, computational models shed light on the strength and directionality properties of the interaction, which was also further characterized from a charge-density perspective using Bader’s “Atoms in Molecules” (AIM) theory, Non-Covalent Interaction plot (NCIplot) visual index, and Natural Bonding Orbital (NBO) analyses. As far as our knowledge extends, this is the first time that ChBs in SAM–RNA complexes have been systematically analyzed, and we believe the results might be useful for scientists working in the field of RNA engineering and chemical biology as well as to increase the visibility of the interaction among the biological community.

Molecular Electrostatic Potential and Noncovalent Interactions in Derivatives of Group 8 Elements

This communication reports experimental and theoretical evidences of σ‐hole interactions in adducts between nitrogen or oxygen nucleophiles and tetroxides of osmium or other group 8 elements. Cocrystals between pyridine or pyridine N‐oxide derivatives and osmium tetroxide are characterized through various techniques and rationalized as σ‐hole interactions using DFT calculations and several other computational tools. We propose the term “osme bond” (OmB, Om=Fe, Ru, Os, (Hs)) for naming the noncovalent interactions wherein group 8 elements have the role of the electrophile. The word osme is the transcription of ὀσμή, the ancient Greek word for smell that was used to name the heaviest group 8 element in relation to the smoky odor of its tetroxide.

Halogen Interactions in Halogenated Oxindoles: Crystallographic and Computational Investigations of Intermolecular Interactions

X-ray structural determinations and computational studies were used to investigate halogen interactions in two halogenated oxindoles. Comparative analyses of the interaction energy and the interaction properties were carried out for Br···Br, C-H···Br, C-H···O and N-H···O interactions. Employing Møller-Plesset second-order perturbation theory (MP2) and density functional theory (DFT), the basis set superposition error (BSSE) corrected interaction energy (E_int(BSSE)) was determined using a supramolecular approach. The E_int(BSSE) results were compared with interaction energies obtained by Quantum Theory of Atoms in Molecules (QTAIM)-based methods. Reduced Density Gradient (RDG), QTAIM and Natural bond orbital (NBO) calculations provided insight into possible pathways for the intermolecular interactions examined. Comparative analysis employing the electron density at the bond critical points (BCP) and molecular electrostatic potential (MEP) showed that the interaction energies and the relative orientations of the monomers in the dimers may in part be understood in light of charge redistribution in these two compounds.

Azine Steric Hindrances Switch Halogen Bonding to N-Arylation upon Interplay with σ-Hole Donating Haloarenenitriles

An interplay between 4-bromo- and 4-iodo-5-nitrophthalonitriles (XNPN, X=Br or I) and any one of the azines (pyridine 1, 4-dimethylaminopyridine 2, isoquinoline 3, 4-cyanopyridine 4, 2-methylpyridine 5, 2-aminopyridine 6, quinoline 7, 1-methylisoquinoline 8, and 2,2'-bipyridine 9) proceeds differently depending on steric and electronic effects of the heterocycles. Sterically unhindered azines 1-3 underwent N-arylation to give the corresponding azinium salts (characterized by ¹ H and ¹³ C{H} NMR and high-resolution ESI-MS). In contrast, azines 4-9 with sterically hindered N atoms or bearing an electron-withdrawing substituent, form stable co-crystals with XNPN, where two interacting molecules are bound by halogen bonding. In all obtained co-crystals, X⋅⋅⋅N structure-directed halogen bonds were recognized and theoretically evaluated including DFT calculations (PBE0-D3/def2-TZVP level of theory), QTAIM analysis, molecular electrostatic potential surfaces, and noncovalent interaction plot index. Estimated energies of halogen bonding vary from -7.6 kcal/mol (for 6 ⋅ INPN) to -11.4 kcal/mol (5 ⋅ INPN).

Chalcogen bonds: Hierarchical ab initio benchmark and density functional theory performance study

We have performed a hierarchical ab initio benchmark and DFT performance study of D2 Ch•••A- chalcogen bonds (Ch = S, Se; D, A = F, Cl). The ab initio benchmark study is based on a series of ZORA-relativistic quantum chemical methods [HF, MP2, CCSD, CCSD(T)], and all-electron relativistically contracted variants of Karlsruhe basis sets (ZORA-def2-SVP, ZORA-def2-TZVPP, ZORA-def2-QZVPP) with and without diffuse functions. The highest-level ZORA-CCSD(T)/ma-ZORA-def2-QZVPP counterpoise-corrected complexation energies (ΔECPC ) are converged within 1.1-3.4 kcal mol-1 and 1.5-3.1 kcal mol-1 with respect to the method and basis set, respectively. Next, we used the ZORA-CCSD(T)/ma-ZORA-def2-QZVPP (ΔECPC ) as reference data for analyzing the performance of 13 different ZORA-relativistic DFT approaches in combination with the Slater-type QZ4P basis set. We find that the three-best performing functionals are M06-2X, B3LYP, and M06, with mean absolute errors (MAE) of 4.1, 4.2, and 4.3 kcal mol-1 , respectively. The MAE for BLYP-D3(BJ) and PBE amount to 8.5 and 9.3 kcal mol-1 , respectively.

Strength of the [Z-I···Hal]- and [Z-Hal···I]- Halogen Bonds: Electron Density Properties and Halogen Bond Length as Estimators of Interaction Energy

Bond energy is the main characteristic of chemical bonds in general and of non-covalent interactions in particular. Simple methods of express estimates of the interaction energy, Eint, using relationships between Eint and a property which is easily accessible from experiment is of great importance for the characterization of non-covalent interactions. In this work, practically important relationships between Eint and electron density, its Laplacian, curvature, potential, kinetic, and total energy densities at the bond critical point as well as bond length were derived for the structures of the [Z-I···Hal]- and [Z-Hal···I]- types bearing halogen bonds and involving iodine as interacting atom(s) (totally 412 structures). The mean absolute deviations for the correlations found were 2.06-4.76 kcal/mol.

Synthesis and characterization of N-heterocyclic carbene-MOEt2 complexes (M = Cu, Ag, Au). Analysis of solvated auxiliary-ligand free [(NHC)M]+ species

We report the synthesis, characterization and computational analysis of coinage metal-ether complexes supported by N-heterocyclic carbenes (NHC), SIPr and Et2CAAC. The related water adducts are also included. The [(NHC)M]+(M = Cu, Ag, Au) species show the noteworthy ability to bind Et2O and H2O. This interaction towards Et2O and H2O is partly ascribed to a σ-hole bonding with an almost linear disposition, taking advantage of the enhanced σ-hole potential evaluated for such [(NHC)M]+ species. This enhanced ability is larger than those found for non-covalent interactions involving main group species.

DFT and IsoStar Analyses to Assess the Utility of σ- and π-Hole Interactions for Crystal Engineering

The interpretation of 36 charge neutral 'contact pairs' from the IsoStar database was supported by DFT calculations of model molecules 1-12, and bimolecular adducts thereof. The 'central groups' are σ-hole donors (H2 O and aromatic C-I), π-hole donors (R-C(O)Me, R-NO2 and R-C6 F5 ) and for comparison R-C6 H5 (R=any group or atom). The 'contact groups' are hydrogen bond donors X-H (X=N, O, S, or R2 C, or R3 C) and lone-pair containing fragments (R3 C-F, R-C≡N and R2 C=O). Nearly all the IsoStar distributions follow expectations based on the electrostatic potential of the 'central-' and 'contact group'. Interaction energies (ΔEBSSE ) are dominated by electrostatics (particularly between two polarized molecules) or dispersion (especially in case of large contact area). Orbital interactions never dominate, but could be significant (∼30 %) and of the n/π→σ*/π* kind. The largest degree of directionality in the IsoStar plots was typically observed for adducts more stable than ΔEBSSE ≈-4 kcal⋅mol-1 , which can be seen as a benchmark-value for the utility of an interaction in crystal engineering. This benchmark could be met with all the σ- and π-hole donors studied.

Crystal engineering with pyrazolyl-thiazole derivatives: structure-directing role of π-stacking and σ-hole interactions

The structure-directing role of a variety of noncovalent interactions in the solid state of pyrazolyl-thiazole derivatives has been analyzed energetically using DFT calculations and by means of Hirshfeld surface analysis.

Resonance Assisted Chalcogen Bonding as a New Synthon in the Design of Dyes

Intramolecular chalcogen bonding in arylhydrazones of sulfamethizole is strengthened by conjugation in the π-system of a noncovalent five-membered ring. The S⋅⋅⋅O distance in the sulfamethizole moiety of these compounds ranges from 2.698(3) to 2.806(15) Å, which indicates its strong dependence on the attached arylhydrazone fragments. Information on the nature of the intramolecular chalcogen bond was afforded by DFT calculations.

F-Halogen Bond: Conditions for Its Existence

The question as to whether the F atom can engage in a halogen bond (XB) remains unsettled. This issue is addressed via density functional theory calculations which pair a wide range of organic and inorganic F-acids with various sorts of Lewis bases. From an energetic perspective, perfluorinated hydrocarbons with sp, sp², or sp³ C-hybridization are unable to form an XB with an N-base, but a very weak bond can be formed if electron-withdrawing C≡N substituents are added to the acid. There is little improvement for inorganic acids O₂NF, FOF, ClF, BrF, SiF₄, or GeF₄, but F₂ is capable of a stronger XB of up to 5 kcal/mol. These results are consistent with a geometric criterion, which compares the intermolecular equilibrium distance with the sum of atomic van der Waals radii. The intensity of the σ-hole on the F atom has predictive value in that a V_s,max of at least 10-15 kcal/mol is required for XB formation. Adding a positive charge to the Lewis acid enhances the strength of any XB and even more so if the base is anionic. The acid-base interaction induces a contraction of the r(AF) covalent bond in the acid in most cases and a deshielding of the NMR signal of the F nucleus.

Engineering Crystals Using sp3 -C Centred Tetrel Bonding Interactions

1,1,2,2-Tetracyanocyclopropane derivatives 1 and 2 were designed and synthesized to probe the utility of sp3 -C centred tetrel bonding interactions in crystal engineering. The crystal packing of 1 and 2 and their 1,4-dioxane cocrystals is dominated by sp3 -C(CN)2 ⋅⋅⋅O interactions, has significant C⋅⋅⋅O van der Waals overlap (≤0.266 Å) and DFT calculations indicate interaction energies of up to -11.0 kcal mol-1 . A cocrystal of 2 with 1,4-thioxane reveals that the cyclopropane synthon prefers interacting with O over S. Computational analyses revealed that the electropositive C2 (CN)4 pocket in 1 and 2 can be seen as a strongly directional 'tetrel-bond donor', similar to halogen bond or hydrogen bond donors. This disclosure is expected to have implications for the utility of such 'tetrel bond donors' in molecular disciplines such as crystal engineering, supramolecular chemistry, molecular recognition and medicinal chemistry.

Quantifying Intramolecular Halogen Bonds in Nucleic Acids: A Combined Protein Data Bank and Theoretical Study

In this study, we report experimental (Protein Data Bank (PDB) search) and theoretical (RI-MP2/def2-TZVP level of theory) evidence of the nature, stability, and directionality properties of intramolecular halogen bonding interactions (HaBs) between 5-bromo/5-iodoracil bases and backbone phosphate groups in nucleic acids (NAs). A PDB survey revealed relevant examples where intramolecular HaBs are undertaken and serve as a structural source of stability in RNA and DNA molecules. In order to develop suitable energy predictors, we started this investigation by calculating the interaction energy values and both the potential V(r) and kinetic G(r) energy densities (using Bader's "atoms in molecules" theory) of several halogen bond complexes involving 5-bromo/5-iodoracil molecules and biologically relevant electron donors. Once the energy predictors based on V(r)/G(r) energy densities were developed, we analyzed the HaBs observed in the biological examples retrieved from the PDB search in order to estimate the strength of the interaction. As far as our knowledge extends, intramolecular halogen bonds in NAs have not been previously quantified in the literature using this methodology and may be of great importance in understanding their structural properties as well as in the construction of molecular materials with DNA and other biological macromolecules.

Anion Recognition by Neutral Chalcogen Bonding Receptors: Experimental and Theoretical Investigations

The utilization of neutral receptors for the molecular recognition of anions based on chalcogen bonding (ChB) is an undeveloped area of host-guest chemistry. In this manuscript, the synthesis of two new families of sulfur, selenium, and tellurium-based ChB binding motifs are reported. The stability of the thiophene, selenophene, and tellurophene binding motifs has enabled the determination of the association constants for ChB halide anion binding in the polar aprotic solvent THF by ¹ H, ⁷⁷ Se, and ¹²⁵ Te NMR experiments. Two different aromatic cores are used and one or two Ch-binding motifs are incorporated with the purpose of encapsulating the anion, offering up to two concurrent chalcogen bonds. Theoretical calculations and NMR experiments reveal that, for S and Se receptors, hydrogen-bonding interactions involving the acidic H atom adjacent to the chalcogen atom are energetically favored over the ChB interaction. However, for the tellurophene binding motif, the σ-hole interaction is competitive and more favored than the hydrogen bond.

Intermolecular interactions in antipyrine-like derivatives 2-halo-N-(1,5-dimethyl-3-oxo-2-phenyl-2,3-dihydro-1H-pyrazol-4-yl)benzamides: X-ray structure, Hirshfeld surface analysis and DFT calculations

Based on experimental and computational data, a complex network of intermolecular interactions has been rationalized for antipyrine compounds.

Direct investigation of chalcogen bonds by multinuclear solid-state magnetic resonance and vibrational spectroscopy

We report a multifaceted experimental and computational study of three self-complementary chalcogen-bond donors as well as a series of seven chalcogen bonded cocrystals.

Diaryliodonium as a double σ-hole donor: the dichotomy of thiocyanate halogen bonding provides divergent solid state arylation by diaryliodonium cations

The reactivity of [Ar¹Ar²I](SCN) toward the solid-state arylation depends on the preorganization of halogen bond (XB)-bound SCN⁻: N-XB-bound thiocyanates, which, in contrast to N,S-XB-bound, undergoes the extremely rare N-arylation of SCN⁻.

Relevant π-hole tetrel bonding interactions in ethyl 2-triazolyl-2-oxoacetate derivatives: Hirshfeld surface analysis and DFT calculations

We report the synthesis and X-ray characterization of four triazole derivatives that include an α-ketoester functionality and two phenyl substituents. The compounds form self-assembled dimers in the solid state establishing two symmetrically equivalent O⋯π-hole interactions.

Charge assisted halogen and pnictogen bonds: insights from the Cambridge Structural Database and DFT calculations

This manuscript combines a search in the Cambridge Structural Database and DFT calculations to analyse the existence and importance of charge assisted pnictogen and halogen bonds in halophosphonium salts.