Aerogen Bonding Interaction: A New Supramolecular Force?

Noble Gas Anions: An Overview of Strategies and Bonding Motifs

This review article aims to provide an overview of the strategies employed to prepare noble gas anions under different environments and experimental conditions, and of the bonding motifs typically occurring in these species. Observed systems include anions fixed into synthesized salts, detected in the gas phase or in high-pressure devices. The major role of the theoretical calculations is also highlighted, not only in support of the experiments, but also as effective in predicting still unreported species. The chemistry of noble gas anions overall appears as a varied and rich paint, offering fascinating opportunities for both experimentalists and theoreticians.

σ-Hole, lone-pair-hole, and π-hole site-based interactions in aerogen-comprising complexes: a comparative study

Herein, the potential of ZO3 and ZF2 aerogen-comprising molecules (where Z = Ar, Kr, and Xe) to engage in σ-, lp-, and π-hole site-based interactions was comparatively studied using various ab initio computations. For the first time, a premier in-depth elucidation of the external electric field (EEF) influence on the strength of the σ-, lp-, and π-hole site-based interactions within the ZO3/ZF2⋯NH3 and ⋯NCH complexes was addressed using oriented EEF with disparate magnitude. Upon the energetic features, σ-hole site-based interactions were noticed with the most prominent preferability in comparison to lp- and π-hole analogs. This finding was ensured by the negative interaction energy values of -11.65, -3.50, and -2.74 kcal mol-1 in the case of σ-, lp-, and π-hole site-based interactions within the XeO3⋯ and XeF2⋯NH3 complexes, respectively. Detailedly, the strength of the σ- and lp-hole site-based interactions directly correlated with the atomic size of the aerogen atoms and the magnitude of the positively oriented EEF. Unexpectedly, an irregular correlation was noticed for the interaction energies of the π-hole site-based interactions with the size of the π-hole. Interestingly, the π-hole site-based interactions within Kr-comprising complexes exhibited higher negative interaction energies than the Ar- and Xe-comprising counterparts. Notwithstanding, a direct proportion between the interaction energies of the π-hole site-based interactions and π-hole size was obtained by employing EEF along the positive orientation with high strength. The present outcomes would be a fundamental basis for forthcoming progress in studying the σ-, lp-, and π-hole site-based interactions within aerogen-comprising complexes and their pertinent applications in materials science and crystal engineering.

Unveiling the Inverse Sandwich Complexes of XeO3: A Computational Exploration

Our study introduces the design of inverse sandwich (iSw) complexes incorporating a noble gas compound: xenon trioxide (XeO3). Through comprehensive computational analyses, we have investigated the critical factors influencing their stability by employing a variety of state-of-the-art computational tools. We demonstrated that the coordination number of xenon in the iSw complex of XeO3 with 18-crown-6 is influenced by the presence of a rare, weakly stabilizing Xe···Xe interaction between the XeO3 molecules. Furthermore, we observed that the stability of iSw complexes of 1,3,5-triphenylbenzene (TPB) and its derivatives is not solely attributed to aerogen bonding, but also involves contributions from C-H···O interactions and back-donation from the lone pair of Xe to the antibonding C-C orbitals of TPB. Additionally, the significant contributions from orbital interactions and dispersion interactions in the TPB derivatives highlight the multifaceted amphoteric properties of XeO3 and reveal that the iSw complexes of TPB and derivatives are not predominantly governed by electrostatic interactions, contrary to conventional belief.

Structural Elucidation of N2O Clusters at Low Temperatures: Exemplary Framework Stabilized by π-Hole-Driven N···O and N···N Pnicogen Bonding Interactions

N2O is a classic prototype, in which central nitrogen is sufficiently electropositive with a positive potential of 20 kcal mol-1 in magnitude to qualify it as a possible pnicogen. This was applied to a test with N2O clusters using ab initio calculations in association with various molecular topographic tools. The structure of the energetically dominant and N2O dimer was in favor of a perpendicular geometry, where the central nitrogen atom of the N2O submolecule assumed a near 90° angle with the adjacent N═O and/or N═N moiety, which provides the affirmation of central nitrogen as a possible π-hole-driven pnicogen. The terminal nitrogen and oxygen atoms of N2O continue to act as conventional electron donors (Lewis bases) with a negative potential. Overall, predominant π-hole-driven N···O and N···N pnicogen bonding interactions were observed to stabilize N2O clusters. Furthermore, N2O clusters (dimers and trimers) were synthesized at low temperatures in an Ar matrix using molecular beam (effusive and supersonic expansion) experiments. The geometries of these clusters were characterized by probing infrared spectroscopy with corroboration from ab initio computational methods. In addition to our previously investigated nitromethane and nitrobenzene systems, N2O also makes it to a pnicogen bonder's club with the central nitrogen as a π-hole-driven pnicogen.

Halogen Bond via an Electrophilic π-Hole on Halogen in Molecules: Does It Exist?

This study reveals a new non-covalent interaction called a π-hole halogen bond, which is directional and potentially non-linear compared to its sister analog (σ-hole halogen bond). A π-hole is shown here to be observed on the surface of halogen in halogenated molecules, which can be tempered to display the aptness to form a π-hole halogen bond with a series of electron density-rich sites (Lewis bases) hosted individually by 32 other partner molecules. The [MP2/aug-cc-pVTZ] level characteristics of the π-hole halogen bonds in 33 binary complexes obtained from the charge density approaches (quantum theory of intramolecular atoms, molecular electrostatic surface potential, independent gradient model (IGM-δginter)), intermolecular geometries and energies, and second-order hyperconjugative charge transfer analyses are discussed, which are similar to other non-covalent interactions. That a π-hole can be observed on halogen in halogenated molecules is substantiated by experimentally reported crystals documented in the Cambridge Crystal Structure Database. The importance of the π-hole halogen bond in the design and growth of chemical systems in synthetic chemistry, crystallography, and crystal engineering is yet to be fully explicated.

Evaluation of Concomitant Halogen and Pnictogen Bonds in Cocrystals of Imines Derived from 2-Nitrobenzaldehyde and 4-Haloaniline

Three imines have been prepared by condensation of 2-nitrobenzaldehyde and 4-haloanilines (halo = Cl, Br, and I) with functionalities that enabled them to act as both halogen and pnictogen bond donors; however, both interactions were found to be absent in the solid state. The prepared imines were further cocrystallized with 1,3-diiodotetrafluorobenzene and 1,3,5-triiodotetrafluorobenzene as halogen bond donors. Six novel cocrystals were prepared by means of liquid-assisted mechanochemical synthesis and by crystallization from solution. All six cocrystals were of 1:1 stoichiometry and comprised a N···I halogen bond between an iodine atom of the perhalogenated halogen bond donor and the imino nitrogen atom of the imine acting as an acceptor. Additionally, in all six cocrystals, the imine molecules were interconnected by NO2···NO2 pnictogen bonding interactions. Computational analysis has shown that the NO2···NO2 exhibits bond critical point electron densities in the region (4.897-8.306) × 10-3 e Å-3 and interaction energies of 23.6-27.7 kJ mol-1, whereas the N···I halogen bonds generally have higher critical point electron densities ((1.795-1.937) × 10-2 e Å-3), but the corresponding total interaction energies are lower (19.4-20.4 kJ mol-1). Statistical analysis of the appearance of NO2···NO2 contacts concomitantly with halogen or hydrogen bonds seems to indicate that there is a positive correlation between the presence of NO2···NO2 pnictogen bonding interactions and other directional interactions in crystal structures.

Solvent effects in anion recognition

Anion recognition is pertinent to a range of environmental, medicinal and industrial applications. Recent progress in the field has relied on advances in synthetic host design to afford a broad range of potent recognition motifs and novel supramolecular structures capable of effective binding both in solution and at derived molecular films. However, performance in aqueous media remains a critical challenge. Understanding the effects of bulk and local solvent on anion recognition by host scaffolds is imperative if effective and selective detection in real-world media is to be viable. This Review seeks to provide a framework within which these effects can be considered both experimentally and theoretically. We highlight proposed models for solvation effects on anion binding and discuss approaches to retain strong anion binding in highly competitive (polar) solvents. The synthetic design principles for exploiting the aforementioned solvent effects are explored.

Unraveling pnicogen bonding cooperativity: Insights from molecular electrostatic potential analysis

A theoretical investigation on the cooperativity of a series of binary, ternary, and quaternary complexes interconnected by pnicogen bonds has been conducted using calculations at the M06-2X/aug-cc-pVTZ level of density functional theory. By measuring changes in the molecular electrostatic potential (MESP) at the nucleus of interacting atoms in all of the complexes, it is possible to quantify the substantial reorganization of the electron density triggered by the formation of pnicogen bonds. The positive change in MESP, indicating a loss of electron density from the donor molecule in a dimer, facilitates the acceptance of electron density from a third molecule, resulting in the formation of a ternary complex with a stronger pnicogen bond compared to the one present in the binary complex. Similarly, the acceptor molecule in a dimer with a negative change in MESP showed an enhanced tendency to donate electron density to an electron-deficient third molecule. The MESP analysis provided valuable insights into the donor/acceptor characteristics of pnicogen bonds within the quaternary complexes. The proposed MESP hypotheses are consistent with the positive cooperativity observed in the pnicogen-bonded clusters. To quantify the changes in MESP, both at the donor atom (ΔVdonor ) and the acceptor atom (ΔVacceptor ), for all pnicogen bonds in the cluster, the total change in MESP (ΔΔVn ) was measured as ΔΔVn  = ∑(ΔVdonor )-∑(ΔVacceptor ). Remarkably, ΔΔVn exhibited a strong linear relationship with the sum of the bond energies of the pnicogen bonds in the cluster. This establishes the MESP analysis as a robust approach for understanding the strength and cooperative behavior of pnicogen-bonded clusters. Additionally, the MESP features provided clear evidence of pnicogen bond formation, further supporting the reliability of this approach.

Exploring Unconventional σ-Hole Interactions: Computational Insights into the Interaction of XeO3 with Non-Aromatic Coordinating Solvents

In order to control the explosiveness and shock sensitivity of XeO3 , we have investigated its plausible interaction with various non-aromatic coordinating solvents, serving as potential Lewis base donors, through density functional theory (DFT) calculations. Out of twenty six such solvents, the top ten were thus identified and then thoroughly examined by employing various computational tools such as the mapping of the electrostatic potential surface (MESP), Wiberg bond indices (WBIs), non-covalent interaction (NCI) plots, Bader's theory of atoms-in-molecules (AIM), natural bond orbital (NBO) analysis, and the energy decomposition analysis (EDA). The amphoteric nature of XeO3 was also explored by investigating the extent of back donation from the lone pair of Xe to the antibonding orbital of the donating atom/group of the solvent molecules. The C-H…O interactions were also found to be a contributing factor in the stabilization of these adducts. Although these aerogen-bonding interactions were found to be predominantly electrostatic, significant contributions from the orbital contributions, as well as dispersion interactions, were observed. The top three non-aromatic solvents (among the twenty six studied) which form the strongest adducts with XeO3 are proposed to be hexamethylphosphoramide (HMPA), N,N'-dimethylpropyleneurea (DMPU) and tetramethylethylenediamine (TMEDA).

σ-Hole Site-Based Interactions within Hypervalent Pnicogen, Halogen, and Aerogen-Bearing Molecules with Lewis Bases: A Comparative Study

σ-Hole site-based interactions in the trigonal bipyramidal geometrical structure of hypervalent pnicogen, halogen, and aerogen-bearing molecules with pyridine and NCH Lewis bases (LBs) were comparatively examined. In this respect, the ZF5···, XF3O2···, and AeF2O3···LB complexes (where Z = As, Sb; X = Br, I; Ae = Kr, Xe; and LB = pyridine and NCH) were investigated. The electrostatic potential (EP) analysis affirmations outlined the occurrence of σ-holes on the systems under consideration with disparate magnitudes that increased according to the following order: AeF2O3 < XF3O2 < ZF5. In line with EP outcomes, the proficiency of σ-hole site-based interactions increased as the atomic size of the central atom increased with a higher favorability for the pyridine-based complexes over NCH-based ones. The interaction energy showed the most favorable negative values of -35.97, -44.53, and -56.06 kcal/mol for the XeF2O3···, IF3O2···, and SbF5···pyridine complexes, respectively. The preferentiality pattern of the studied interactions could be explained as a consequence of (i) the dramatic rearrangement of ZF5 molecules from the trigonal bipyramid geometry to the square pyramidal one, (ii) the significant and tiny deformation energy in the case of the interaction of XF3O2 molecules with pyridine and NCH, respectively, and (iii) the absence of geometrical deformation within the AeF2O3···pyridine and ···NCH complexes other than the XeF2O3···pyridine one. Quantum theory of atoms in molecules and noncovalent interaction index findings reveal the partially covalent nature of most of the investigated interactions. Symmetry-adapted perturbation theory affirmations declared that the electrostatic component was the driving force beyond the occurrence of the considered interactions. The obtained findings will help in improving our understanding of the effect of geometrical deformation on intermolecular interactions.

Size Matters: Computational Insights into the Crowning of Noble Gas Trioxides

In pursuit of enhancing the stability of the highly explosive and shock-sensitive compound XeO3, we performed quantum chemical calculations to investigate its possible complexation with electron-rich crown ethers, including 9-crown-3, 12-crown-4, 15-crown-5, 18-crown-6, and 21-crown-7, as well as their thio analogues. Furthermore, we expanded our study to other noble gas trioxides (NgO3), namely, KrO3 and ArO3. The basis set superposition error (BSSE) corrected interaction energies for these adducts range from -13.0 kcal/mol to -48.2 kcal/mol, which is notably high for σ-hole-mediated noncovalent interactions. The formation of these adducts was observed to be more favorable with the increase in the ring size of the crowns and less favorable while going from XeO3 to ArO3. A comprehensive analysis by various computational tools such as the mapping of the electrostatic potential (ESP), Wiberg bond indices (WBIs), Bader's theory of atoms-in-molecules (AIM), natural bond orbital (NBO) analysis, noncovalent interaction (NCI) plots, and energy decomposition analysis (EDA) revealed that the C-H···O interactions, as well as dispersion interactions, play a pivotal role in stabilizing adducts involving larger crowns. A noteworthy outcome of our study is the revelation of a coordination number of 9 for xenon in the complex formed between XeO3 and the thio analogue of 18-crown-6, which is higher than the largest number reported to date.

Anisotropies in electronic densities and electrostatic potentials of Halonium Ions: focus on Chlorine, Bromine and Iodine

CONTEXT

Why are the halonium cations so effective in forming strongly-bound complexes? We directed our research to address this question and we present electrostatic potential data for the valence-state halogen atoms X and halonium cations X+, where X = Cl, Br, I. The electron densities and electrostatic potentials of the halonium cations show considerably greater anisotropy than do the valence state halogens. The distances from the electrostatic potential surface maxima to the halogen nuclei are about 0.5 Å smaller than the distances from the electrostatic potential surface minima to the nuclei, giving the halonium cations each a more disk-like shape than the corresponding neutral valence state halogens. Their surface electrostatic potentials are totally consistent with the directionalities of halonium cations in complexes and the strengths of their interactions. To add perspective to this brief report, we have included calculations of the isotropic cation K+ and noble gas Kr.

METHODS

The calculations of the electrostatic potentials of the valence states of the halogen atoms Cl, Br and I and the halonium cations Cl+, Br+ and I+, as well as K+ and Kr, on 0.001 au contours of their electronic densities were carried out with Gaussian O9 and the Wave Function Analysis - Surface Analysis Suite (WFA-SAS) at the M06-2X/6-31 + G(d,p) and M06-2X/3-21G* levels.

The Formation of σ-Hole Bonds: A Physical Interpretation

This paper discusses two quite different computational experiments relating to the formation of σ-hole bonds A···B. The first involves looking at the complex at equilibrium and finding the contour X of the electronic density which allows the iso-density envelopes of A and B to be nearly touching. This contour increases, becoming closer to the nuclei, as the strength of the interaction increases. The second experiment involves allowing A and B to approach each other, with the aim of finding the distance at which their 0.001 a.u. iso-density envelopes are nearly merging into one envelope. What is found in the second experiment may be somewhat surprising, in that the ratio of the distance between interacting atoms at this nearly merging point-divided by the sum of the van der Waals radii of these atoms-covers a narrow range, typically between 1.2 and 1.3. It is intriguing to note that for the dataset presented, approaching molecules attracted to each other appear to do so unknowing of the strength of their ultimate interaction. This second experiment also supports the notion that one should expect favorable interactions, in some instances, to have close contacts significantly greater than the sums of the van der Waals radii.

Ye Olde supramolecular chemistry, its modern rebranding and overarching trends in chemistry

We can describe current contingency of supramolecular chemistry as "post-halogen bonding", with clear reference to the success of the σ-hole model and the halogen bond concepts. This phase is characterized by a strong push towards a new nomenclature for non-covalent interactions, a group-by-group one focusing on the electrophile. As such nomenclature increasingly meets IUPAC endorsement, its proposers report resistances to such ideas, especially in the inorganic and coordination chemistry communities. The whole issue has been generating considerable debate in the last decade. Herein we fully embrace such discussion in the hope of involving a larger share of the relevant communities. Alternative descriptions are here reevaluated, novel views reconnected with older ones, and it is ultimately questioned whether the introduction of such a nomenclature and its subtending ideas would be beneficial. The themes of appreciation of general trends in chemistry, of counterintuitive interactions, of positioning of novel nomenclature with respect to existing ones, and of the extension of group-by-group naming from main block to d-block elements - as key and currently unresolved issues - are discussed. Equivalent, alternative and arguably more comprehensive descriptions are tentatively given, in the hope to overcome controversies together in the pursuit of higher rewards: a comprehensive shared view of supramolecular forces and a common language to express it.

Chalcogen bonding catalysis

Catalysts play a major role in chemical synthesis, and catalysis is considered to be a green and economic process. Catalysis is dominated by covalent interactions between the catalyst and substrate. The design of non-covalent catalysts came into limelight only recently. Hydrogen bonding (HB) catalysts are well established among non-covalent catalysts, including asymmetric HB catalysts. Though halogen bonding (XB) catalysis and its asymmetric version are gaining admiration, non-covalent chalcogen bonding catalysis (ChB) is in the budding stage. This tutorial review will focus on the recently evolved chalcogen bonding catalysis and emphasis will be given to the chalcogen bonding of chiral molecules. Since successful enantioselective chalcogen bonding catalysis is yet to be reported, this review will focus on the basics of non-covalent bonding catalysis, chalcogen bonding catalysis, chiral chalcogenide synthesis, rigidification of transition states by ChB, stabilization of cations by chiral chalcogens, details of unsuccessful asymmetric chalcogen bonding catalysis, enantioseparation of racemic molecules using ChB, and the existence of ChB in chiral biomolecules.

Hole interactions of aerogen oxides with Lewis bases: an insight into σ-hole and lone-pair-hole interactions

σ-Hole and lone-pair (lp)-hole interactions of aerogen oxides with Lewis bases (LB) were comparatively inspected in terms of quantum mechanics calculations. The ZOn ⋯ LB complexes (where Z = Kr and Xe, n = 1, 2, 3 and 4, and LB = NH3 and NCH) showed favourable negative interaction energies. The complexation features were explained in light of σ-hole and lp-hole interactions within optimum distances lower than the sum of the respective van der Waals radii. The emerging findings outlined that σ-hole interaction energies generally enhanced according to the following order: KrO4 ⋯ < KrO⋯ < KrO3⋯ < KrO2⋯LB and XeO4⋯ < XeO⋯ < XeO2⋯ < XeO3⋯LB complexes with values ranging from -2.23 to -12.84 kcal mol-1. Lp-hole interactions with values up to -5.91 kcal mol-1 were shown. Symmetry-adapted perturbation theory findings revealed the significant contributions of electrostatic forces accounting for 50-65% of the total attractive forces within most of the ZOn⋯LB complexes. The obtained observations would be useful for the understanding of hole interactions, particularly for the aerogen oxides, with application in supramolecular chemistry and crystal engineering.

Spodium Bonds Involving Methylmercury and Ethylmercury in Proteins: Insights from X-ray Analysis and Computations

In this study, the stability, directionality, and physical nature of Spodium bonds (SpBs, an attractive noncovalent force involving elements from group 12 and Lewis bases) between methylmercury (MeHg) and ethylmercury (EtHg) and amino acids (AAs) have been analyzed from both a structural (X-ray analysis) and theoretical (RI-MP2/def2-TZVP level of theory) point of view. More in detail, an inspection of the Protein Data Bank (PDB) reported evidence of noncovalent contacts between MeHg and EtHg molecules and electron-rich atoms (e.g., O atoms belonging to the protein backbone and S atoms from MET residues or the π-systems of aromatic AAs such as TYR or TRP). These results were rationalized through a computational study using MeHg coordinated to a thiolate group as a theoretical model and several neutral and charged electron-rich molecules (e.g., benzene, formamide, or chloride). The physical nature of the interaction was analyzed from electrostatics and orbital perspectives by performing molecular electrostatic potential (MEP) and natural bonding orbital (NBO) analyses. Lastly, the noncovalent interactions plot (NCIplot) technique was used to provide a qualitative view of the strength of the Hg SpBs and compare them to other ancillary interactions present in these systems as well as to shed light on the extension of the interaction in real space. We believe that the results derived from our study will be useful to those scientists devoted to protein engineering and bioinorganic chemistry as well as to expanding the current knowledge of SpBs among the chemical biology community.

Synthesis, X-ray characterization and DFT analysis of dicyanidoaurate telluronium salts: on the importance of charge assisted chalcogen bonds

In this manuscript we report the synthesis and X-ray characterization of two cyanidoaurate telluronium salts, namely (3-fluorophenyl)(methyl)(phenyl)telluronium dicyanidoaurate [(3-F-Ph)(Me)(Ph)Te][Au(CN)2] (1) and methyldiphenyltelluronium dicyanidoaurate [(Me)(Ph)2Te][Au(CN)2] (2). In the solid state, the tellurium atom establishes three concurrent and directional chalcogen bonds (ChBs) with the adjacent anions, in both compounds. These charge-assisted ChBs (CAChBs) have been analyzed using DFT calculations and several computational tools. The MEP surface analysis discloses the existence of three σ-holes at the Te-atoms capable of establishing strong CAChBs with the counter-ions. In addition, significant charge transfer from the lone pair orbital at the N-atom of the anion to the antibonding σ*(Te-C) orbital of the cation is observed in some cases.

Unified classification of non-covalent bonds formed by main group elements: a bridge to chemical bonding

Using correlation plots of binding energy and electron density at the bond critical point, we investigated the nature of intermolecular non-covalent bonds (D-X⋯A, where D = O/S/F/Cl/Br/H, mostly, X = main group elements (except noble gases), A = H2O, NH3, H2S, PH3, HCHO, C2H4, HCN, CO, CH3OH, and CH3OCH3). The binding energies were calculated at the MP2 level of theory, followed by Atoms in Molecules (AIM) analysis of the ab initio wave functions to obtain the electron density at the bond critical point (BCP). For each non-covalent bond, the slopes of the binding energy versus electron density plot have been determined. Based on their slopes, non-covalent bonds are classified as non-covalent bond closed-shell (NCB-C) or non-covalent bond shared-shell (NCB-S). Intriguingly, extrapolating the slopes of the NCB-C and NCB-S cases leads to intramolecular "ionic" and "covalent" bonding regimes, establishing a link between such intermolecular non-covalent and intramolecular chemical bonds. With this new classification, hydrogen bonds and other non-covalent bonds formed by a main-group atom in a covalent molecule are classified as NCB-S. Atoms found in ionic molecules generally form NCB-C type bonds, with the exception of carbon which also forms NCB-C type bonds. Molecules with a tetravalent carbon do behave like ions in ionic molecules such as NaCl and interact with other molecules through NCB-C type bonds. As with the chemical bonds, there are some non-covalent bonds that are intermediate cases.

N-X⋅⋅⋅O-N Halogen Bonds in Complexes of N-Haloimides and Pyridine-N-oxides: A Large Data Set Study

N-X⋅⋅⋅- O-N+ halogen-bonded systems formed by 27 pyridine N-oxides (PyNOs) as halogen-bond (XB) acceptors and two N-halosuccinimides, two N-halophthalimides, and two N-halosaccharins as XB donors are studied in silico, in solution, and in the solid state. This large set of data (132 DFT optimized structures, 75 crystal structures, and 168 1 H NMR titrations) provides a unique view to structural and bonding properties. In the computational part, a simple electrostatic model (SiElMo) for predicting XB energies using only the properties of halogen donors and oxygen acceptors is developed. The SiElMo energies are in perfect accord with energies calculated from XB complexes optimized with two high-level DFT approaches. Data from in silico bond energies and single-crystal X-ray structures correlate; however, data from solution do not. The polydentate bonding characteristic of the PyNOs' oxygen atom in solution, as revealed by solid-state structures, is attributed to the lack of correlation between DFT/solid-state and solution data. XB strength is only slightly affected by the PyNO oxygen properties [(atomic charge (Q), ionization energy (Is,min ) and local negative minima (Vs,min )], as the σ-hole (Vs,max ) of the donor halogen is the key determinant leading to the sequence N-halosaccharin>N-halosuccinimide>N-halophthalimide on the XB strength.

Halogen Bonding Assembles Anion⋅⋅⋅Anion Architectures in Non-centrosymmetric Iodate and Bromate Crystals

Single crystal X-ray diffraction of iodate and bromate salts shows that the I and Br atoms in IO3 - and BrO3 - anions form short and linear O-I/Br⋅⋅⋅O contacts with the O atoms of nearby anions. Non-centrosymmetric systems are formed wherein anions are orderly aligned into supramolecular 1D and 2D networks. Theoretical evidences, namely the outcome of QTAIM and NCIplot studies, prove the attractive nature of these contacts and the ability of iodate and bromate anions to act as robust halogen bond (HaB) donors. The HaB is proposed as a general and effective assisting tool to control the architecture of acentric iodate salts.

Tetrel bonds involving a CF3 group participate in protein-drug recognition: a combined crystallographic and computational study

In this study, the ability of CF₃ groups to bind to the electron-rich side chains and backbone groups of proteins has been investigated by combining a Protein Data Bank (PDB) survey and ab initio quantum mechanics calculations. More precisely, an inspection of the PDB involving organic ligands containing a CF₃ group and electron-rich atoms (A = N, O and S) in the vicinity revealed 419 X-ray structures exhibiting CF₃⋯A tetrel bonds (TtBs). In a posterior stage, those hits that exhibited the most relevant features in terms of directionality and intermolecular distance were selected for theoretical calculations at the RI-MP2/def2-TZVPD level of theory. Also, Hammett's regression plots of several TtB complexes involving meta- and para-substituted benzene derivatives were computed to shed light on the substituent effects. Moreover, the TtBs were characterized through several state-of-the-art computational techniques, such as the Quantum Theory of Atoms in Molecules (QTAIM) and Noncovalent Interactions plot (NCIplot) methodologies. We believe that the results gathered from our study will be useful for rational drug design and biological communities as well as for further expanding the role of this interaction to biomedical applications.

Regium-π Bonds Involving Nucleobases: Theoretical Study and Biological Implications

In this study, we provide crystallographic (Protein Data Bank (PDB) inspection) and theoretical (RI-MP2/def2-TZVP//PBE0-D3/def2-SVP level of theory) evidence of the involvement of nucleobases in Regium-π bonds (RgBs). This noncovalent interaction involves an electrophilic site located on an element of group 11 (Cu, Ag, and Au) and an electron-rich species (lone pair, LP donor, or π-system). Concretely, an initial PDB search revealed several examples where RgBs were undertaken involving DNA bases and Cu(II), Ag(I), and Au(I/III) ions. While coordination positions (mainly at the N atoms of the base) are well known, the noncovalent binding force between these counterparts has been scarcely studied in the literature. In this regard, 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, noncovalent interaction plot (NCIplot) visual index, and natural bonding orbital (NBO) analyses. As far as our knowledge extends, this is the first time that RgBs in metal-DNA complexes are systematically analyzed, and we believe the results might be useful for scientists working in the field of nucleic acid engineering and chemical biology as well as to increase the visibility of the interaction among the biological community.

Yet another perspective on hole interactions, part II: lp-hole vs. lp-hole interactions

Most of the experimental and theoretical work in hole interactions (HIs) is mainly focused on exploiting the nature and characteristics of σ and π-holes. In this perspective, we focus our attention on understanding the origin and properties of lone-pair holes. These holes are present on an atom opposite to its lone-pair region. Utilizing some new and old examples, such as X₃N/P⋯F^- (X = F/Cl/Br/I), F-Cl/Br/I⋯H₃P⋯NCH and H₃B-NBr₃ along with other molecular systems, we explored to what extent these lp-holes participate in lp-hole interactions, if they participate at all.

Why much of Chemistry may be indisputably non-bonded?

In this compendium, the wide scope of all intermolecular interactions ever known has been revisited, in particular giving emphasis the capability of much of the elements of the periodic table to form non-covalent contacts. Either hydrogen bonds, dihydrogen bonds, halogen bonds, pnictogen bonds, chalcogen bonds, triel bonds, tetrel bonds, regium bonds, spodium bonds or even the aerogen bond interactions may be cited. Obviously that experimental techniques have been used in some works, but it was through the theoretical methods that these interactions were validate, wherein the QTAIM integrations and SAPT energy partitions have been useful in this regard. Therefore, the great goal concerns to elucidate the interaction strength and if the intermolecular system shall be total, partial or non-covalently bonded, wherein this last one encompasses the most majority of the intermolecular interactions what leading to affirm that chemistry is debatably non-bonded.

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.

77Se and 125Te solid-state NMR and X-ray diffraction structural study of chalcogen-bonded 3,4-dicyano-1,2,5-chalcogenodiazole cocrystals

Three novel chalcogen-bonded cocrystals featuring 3,4-dicyano-1,2,5-selenodiazole (C₄N₄Se) or 3,4-dicyano-1,2,5-tellurodiazole (C₄N₄Te) as chalcogen-bond donors and hydroquinone (C₆H₆O₂), tetraphenylphosphonium chloride (C₂₄H₂₀P⁺·Cl^-) or tetraethylphosphonium chloride (C₈H₂₀P⁺·Cl^-) as chalcogen-bond acceptors have been prepared and characterized by single-crystal X-ray diffraction (XRD), powder X-ray diffraction and ⁷⁷Se/¹²⁵Te magic-angle spinning solid-state NMR spectroscopy. The single-crystal XRD results show that the chalcogenodiazole molecules interact with the electron donors through two σ-holes on each of the chalcogen atoms, which results in highly directional and moderately strong chalcogen bonds. Powder XRD confirms that the crystalline phases are preserved upon moderate grinding of the samples for solid-state NMR experiments. Measurement of ⁷⁷Se and ¹²⁵Te chemical shift tensors via magic-angle spinning solid-state NMR spectroscopy confirms the number of magnetically unique chalcogen sites in each asymmetric unit and reveals the impact of chalcogen-bond formation on the local electronic structure. These NMR data are further assessed in the context of analogous data for a wider range of crystalline chalcogen-bonded systems.

Structural, vibrational and electronic properties of some tetrel-bonded complexes of the fluorinated methanes methyl fluoride, difluoromethane and fluoroform: an ab initio study

A search has been conducted, by means of ab initio molecular orbital theory, for potential tetrel-bonded complexes formed between the fluorinated methanes methyl fluoride, difluoromethane and fluoroform, and the related hydrides ammonia, water, hydrogen fluoride, phosphine, hydrogen sulphide and hydrogen chloride. Eleven such complexes have been identified, six containing CH₃F and five CH₂F₂. The complexes are typically less strongly bound than their hydrogen-bonded counterparts, and the interaction energies vary in a consistent way with the periodic trend of the electron donors. The intermolecular separations and changes of the relevant intramolecular bond lengths, the wavenumber shifts of the critical vibrational modes and the extents of charge transfer correlate, by and large, with the strengths of interaction.

Stereoselective Processes Based on σ-Hole Interactions

The σ-hole interaction represents a noncovalent interaction between atoms with σ-hole(s) on their surface (such as halogens and chalcogens) and negative sites. Over the last decade, significant developments have emerged in applications where the σ-hole interaction was demonstrated to play a key role in the control over chirality. The aim of this review is to give a comprehensive overview of the current advancements in the use of σ-hole interactions in stereoselective processes, such as formation of chiral supramolecular assemblies, separation of enantiomers, enantioselective complexation and asymmetric catalysis.

Resonance-assisted intramolecular triel bonds

The possibility that the intramolecular Tr⋯S triel bond is strengthened by resonance is examined by quantum chemical calculations within the planar five-membered ring of TrH₂-CRCR-CRS (Tr = Al, Ga, In; R = NO₂, CH₃). This internal bond is found to be rather short (2.4-2.7 Å) with a large bond energy between 12 and 21 kcal mol^-1. The pattern of bond length alternation and atomic charges within the ring is consistent with resonance involving the conjugated double bonds. This resonance enhances the triel bond strength by some 25%. The electron-withdrawing NO₂ group weakens the bond, but it is strengthened by the electron-donating CH₃ substituent. NICS analysis suggests the presence of a certain degree of aromaticity within the ring.

Non-covalent interactions atlas benchmark data sets 4: σ-hole interactions

The SH250×10 dataset presented here extends the Non-Covalent Interactions Atlas database (https://www.nciatlas.org) to complexes bound by σ-hole interactions - halogen, chalcogen and pnictogen bonds. It comprises 250 complexes where Cl, Br, I, S, Se, P and As interact with diverse electron donors. An accurate CCSD(T)/CBS benchmark is provided for ten points along a dissociation curve of each complex. The SH250×10 set is used in testing a wide variety of DFT functionals and semiempirical quantum-mechanical methods. In DFT calculations, the new data set exposes large errors of some functionals related to exaggerated charge transfer. The size and diversity of the data set have also been exploited in the reparametrization of a halogen-bond correction for the PM6 semiempirical method.

Noncovalent interactions in proteins and nucleic acids: beyond hydrogen bonding and π-stacking

Understanding the noncovalent interactions (NCIs) among the residues of proteins and nucleic acids, and between drugs and proteins/nucleic acids, etc., has extraordinary relevance in biomolecular structure and function. It helps in interpreting the dynamics of complex biological systems and enzymatic activity, which is esential for new drug design and efficient drug delivery. NCIs like hydrogen bonding (H-bonding) and π-stacking have been researchers' delight for a long time. Prominent among the recently discovered NCIs are halogen, chalcogen, pnictogen, tetrel, carbo-hydrogen, and spodium bonding, and n → π* interaction. These NCIs have caught the imaginations of various research groups in recent years while explaining several chemical and biological processes. At this stage, a holistic view of these new ideas and findings lying scattered can undoubtedly trigger our minds to explore more. The present review attempts to address NCIs beyond H-bonding and π-stacking, which are mainly n → σ*, n → π* and σ → σ* type interactions. Five of the seven NCIs mentioned earlier are linked to five non-inert end groups of the modern periodic table. Halogen (group-17) bonding is one of the oldest and most explored NCIs, which finds its relevance in biomolecules due to the phase correction and inhibitory properties of halogens. Chalcogen (group 16) bonding serves as a redox-active functional group of different active sites of enzymes and acts as a nucleophile in proteases and phosphates. Pnictogen (group 15), tetrel (group 14), triel (group 13) and spodium (group 12) bonding does exist in biomolecules. The n → π* interactions are linked to backbone carbonyl groups and protein side chains. Thus, they are crucial in determining the conformational stability of the secondary structures in proteins. In addition, a more recently discovered to and fro σ → σ* type interaction, namely carbo-hydrogen bonding, is also present in protein-ligand systems. This review summarizes these grand epiphanies routinely used to elucidate the structure and dynamics of biomolecules, their enzymatic activities, and their application in drug discovery. It also briefs about the future perspectives and challenges posed to the spectroscopists and theoreticians.

Competition between Intra and Intermolecular Pnicogen Bonds. Complexes between Naphthalene Derivatives and Neutral or Anionic Bases

The PnF 2 (Pn=P,As,Sb,Bi) on a naphthalene scaffold can engage in an internal pnicogen Pn···N bond (PnB) with a NH 2 group placed close to it on the adjoining ring. An approaching neutral NH 3 molecule can engage in a second PnB with the central Pn, which tends to weaken the intramolecular bond. The presence of the latter in turn weakens the intermolecular PnB with respect to that formed in its absence. Replacement of the external NH 3 by a CN - anion causes a fundamental change in the bonding pattern, producing a fourth covalent bond with Pn, which rearranges into a trigonal bipyramidal motif. This addition disrupts the internal Pn···N pnicogen bond, recasting the PnF 2 ···NH 2 interaction into a NH···F H-bond.

Theoretical study on Xe⋯N non-covalent interactions: Three hybridization N with XeO3 and XeOF2

The interactions of complexes of XeOF2 and XeO3 with a series of different hybridization N-containing donors are studied by means of DFT and MP2 calculations. The aerogen bonding interaction energies range from 6.5 kcal/mol to 19.9 kcal/mol between XeO3 or XeOF2 and typical N-containing donors. The sequence of interaction for N-containing hybridization is sp3>sp2>sp, and XeO3 is higher than XeOF2. For some donors of sp2 and sp3 hybridization, the steric effect plays a minor role in the interaction with the evidence of reduced density gradient plots. The dominant stable part is the electrostatic interaction. In complex of XeO3, the weight of polarization is larger than dispersion, while the situation is opposite for XeOF2 complexes. Except for the sum of the maximum value of molecular electrostatic potential on Xe atom and minimum value of molecular electrostatic potential on N atom, the otherfive interaction parameters including the potential energy density at bond critical point, the equilibrium distances, interaction energies with the basis set superposition error correction, localized molecular orbital energy decomposition analysis interaction energies, and the electron charge density, show great linear correlation coefficients with each other.

The Phosphorus Bond, or the Phosphorus-Centered Pnictogen Bond: The Covalently Bound Phosphorus Atom in Molecular Entities and Crystals as a Pnictogen Bond Donor

The phosphorus bond in chemical systems, which is an inter- or intramolecular noncovalent interaction, occurs when there is evidence of a net attractive interaction between an electrophilic region associated with a covalently or coordinately bonded phosphorus atom in a molecular entity and a nucleophile in another, or the same, molecular entity. It is the second member of the family of pnictogen bonds, formed by the second member of the pnictogen family of the periodic table. In this overview, we provide the reader with a snapshot of the nature, and possible occurrences, of phosphorus-centered pnictogen bonding in illustrative chemical crystal systems drawn from the ICSD (Inorganic Crystal Structure Database) and CSD (Cambridge Structural Database) databases, some of which date back to the latter part of the last century. The illustrative systems discussed are expected to assist as a guide to researchers in rationalizing phosphorus-centered pnictogen bonding in the rational design of molecular complexes, crystals, and materials and their subsequent characterization.

Malonaldehyde-like Systems: BeF2 Clusters—A Subtle Balance between Hydrogen Bonds, Beryllium Bonds, and Resonance

The stability of malonaldehyde is governed by intramolecular hydrogen bonds (IMHBs) as well as in malonaldehyde-like systems where oxygen is replaced by N or S at any of the basic sites. As beryllium bonds have been shown to strongly cooperate with hydrogen bonds, this work explores at the high level ab initio G4 level of theory the effect of including this non-covalent interaction in the system through its association with BeF2. Although malonaldehyde follows the expected trends, where the formation of a pseudocyclic form is favored also when IMHB and Be bonds are present, the subtle balance between both non-covalent interactions leads to some surprising results when other heteroatoms are involved, to the point that interaction energies can be much larger than expected or even cyclization is not favored. A complete analysis using different computational tools gives an answer to those cases escaping the predictable trends.

Predictability of Chalcogen-Bond-Driven Crystal Engineering: An X-ray Diffraction and Selenium-77 Solid-State NMR Investigation of Benzylic Selenocyanate Cocrystals

We describe a series of new chalcogen-bonded cocrystals featuring 1,2-bis(selenocyanatomethyl)benzene (DSN) and 1,2,4,5-tetrakis(selenocyanatomethyl)-benzene (TSN) as the donor moieties and a variety of Lewis bases such as onium halides, N-oxides, and pyridine-containing heterocycles as the acceptors. Single-crystal X-ray diffraction demonstrates that, in every case, the selenocyanates consistently interact with the acceptor molecules through strong and directional Se···X chalcogen-bonds (ChBs) (X = halides, oxygen, and nitrogen). ⁷⁷Se solid-state nuclear magnetic resonance spectroscopy was applied to measure selenium chemical shift tensor magnitudes and to explore potential correlations between these tensor elements and the local ChB geometry. In every case, the isotropic ⁷⁷Se chemical shift decreases, and the chemical shift tensor span increases upon cocrystallization of DSN with the various ChB acceptors. This work contributes to a growing body of knowledge concerning the predictability and robustness of chalcogen bonds in crystal engineering as well as the NMR response to the establishment of chalcogen bonds. In particular, among the systems studied here, highly linear chalcogen bonds are formed exclusively at the stronger σ-hole of each and every selenium atom regardless of the size, charge, or denticity of the electron donor moiety.

π-Hole Tetrel Bonds—Lewis Acid Properties of Metallylenes

The MP2/aug-cc-pVTZ calculations were performed on the dihalometallylenes to indicate their Lewis acid and Lewis base sites. The results of the Cambridge Structural Database search show corresponding and related crystal structures where the tetrel center often possesses the configuration of a trigonal bipyramid or octahedron. The calculations were also carried out on dimers of dichlorogermylene and dibromogermylene and on complexes of these germylenes with one and two 1,4-dioxide molecules. The Ge⋯Cl, Ge⋯Br, and Ge⋯O interactions are analyzed. The Ge⋯O interactions in the above mentioned germylene complexes may be classified as the π-hole tetrel bonds. The MP2 calculations are supported by the results of the Quantum Theory of Atoms in Molecules (QTAIM) and the Natural Bond Orbital (NBO) approaches.

Noble-gas compounds: A general procedure of bonding analysis

This paper accounts for a general procedure of bonding analysis that is, expectedly, adequate to describe any type of interaction involving the noble-gas (Ng) atoms. Building on our recently proposed classification of the Ng-X bonds (X = binding partner) [New J. Chem. 44, 15536 (2020)], these contacts are first distinguished into three types, namely, A, B, or C, based on the topology of the electron energy density H(r) and on the shape of its plotted form. Bonds of type B or C are, then, further assigned as B-loose (B^l) or B-tight (B^t) and C-loose (C^l) or C-tight (C^t) depending on the sign that H(r) takes along the Ng-X bond path located from the topological analysis of ρ(r), particularly at around the bond critical point (BCP). Any bond of type A, B^l/B^t, or C^l/C^t is, finally, assayed in terms of contribution of covalency. This is accomplished by studying the maximum, minimum, and average value of H(r) over the volume enclosed by the low-density reduced density gradient (RDG) isosurface associated with the bond (typically, the RDG isosurface including the BCP) and the average ρ(r) over the same volume. The bond assignment is also corroborated by calculating the values of quantitative indices specifically defined for the various types of interactions (A, B, or C). The generality of our taken approach should encourage its wide application to the study of Ng compounds.

Deciphering C–H⋯O/X weak hydrogen bonding and halogen bonding interactions in aromatic peptoids

We deciphered weak interactions in aromatic peptoids, such as C–H⋯O/X, and simultaneously identified strong interactions, including N–H⋯N and N–H⋯O, in this class of foldamer.

Lead(ii) supramolecular structures formed through a cooperative influence of the hydrazinecarbothioamide derived and ancillary ligands

We report on tetrel bonding and other noncovalent interactions in the lead(ii)-derived complexes with the hydrazinecarbothioamide derived and ancillary ligands, which predominantly drive the formation of extended architectures.