Spodium Bonds: Noncovalent Interactions Involving Group 12 Elements

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.

Anion⋯anion self-assembly under the control of σ- and π-hole bonds

The electrostatic attraction between charges of opposite signs and the repulsion between charges of the same sign are ubiquitous and influential phenomena in recognition and self-assembly processes. However, it has been recently revealed that specific attractive forces between ions with the same sign are relatively common. These forces can be strong enough to overcome the Coulomb repulsion between ions with the same sign, leading to the formation of stable anion⋯anion and cation⋯cation adducts. Hydroden bonds (HBs) are probably the best-known interaction that can effectively direct these counterintuitive assembly processes. In this review we discuss how σ-hole and π-hole bonds can break the paradigm of electrostatic repulsion between like-charges and effectively drive the self-assembly of anions into discrete as well as one-, two-, or three-dimensional adducts. σ-Hole and π-hole bonds are the attractive forces between regions of excess electron density in molecular entities (e.g., lone pairs or π bond orbitals) and regions of depleted electron density that are localized at the outer surface of bonded atoms opposite to the σ covalent bonds formed by atoms (σ-holes) and above and below the planar portions of molecular entities (π-holes). σ- and π-holes can be present on many different elements of the p and d block of the periodic table and the self-assembly processes driven by their presence can thus involve a wide diversity of mono- and di-anions. The formed homomeric and heteromeric adducts are typically stable in the solid phase and in polar solvents but metastable or unstable in the gas phase. The pivotal role of σ- and π-hole bonds in controlling anion⋯anion self-assembly is described in key biopharmacological systems and in molecular materials endowed with useful functional properties.

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.

Investigating Recurrent Matere Bonds in Pertechnetate Compounds

In this manuscript we evaluate the X-ray structure of five new pertechnetate derivatives of general formula [M(H2O)4(TcO4)2], M=Mg, Co, Ni, Cu, Zn (compounds 1-5) and one perrhenate compound Zn(H2O)4(ReO4)2 (6). In these complexes the metal center exhibits an octahedral coordination with the pertechnetate units as axial ligands. All compounds exhibit the formation of directional Tc⋅⋅⋅O Matere bonds (MaBs) that propagate the [M(H2O)4(TcO4)2], into 1D supramolecular polymers in the solid state. Such 1D polymers are linked, generating 2D layers, by combining additional MaBs and hydrogen bonds (HBs). Such concurrent motifs have been analyzed theoretically, suggesting the noncovalent σ-hole nature of the MaBs. The interaction energies range from weak (~ -2 kcal/mol) for the MaBs to strong (~ -30 kcal/mol) for the MaB+HB assemblies, where HB dominates. In case of M=Zn, the corresponding perrhenate Zn(H2O)4(ReO4)2 complex, has been also synthesized for comparison purposes, resulting in the formation of an isostructural X-ray structure, corroborating the structure-directing role of Matere bonds.

Halogen-Decorated Metal-Organic Frameworks for Efficient and Selective CO2 Capture, Separation, and Chemical Fixation with Epoxides under Mild Conditions

In the present work, three novel halogen-appended cadmium(II) metal-organic frameworks [Cd2(L1)2(4,4'-Bipy)2]n·4n(DMF) (1), [Cd2(L2)2(4,4'-Bipy)2]n·3n(DMF) (2), and [Cd(L3)(4,4'-Bipy)]n·2n(DMF) (3) [where L1 = 5-{(4-bromobenzyl)amino}isophthalate; L2 = 5-{(4-chlorobenzyl)amino}isophthalate; L3 = 5-{(4-fluorobenzyl)amino}isophthalate; 4,4'-Bipy = 4,4'-bipyridine; and DMF = N,N'-dimethylformamide] have been synthesized under solvothermal conditions and characterized by various analytical techniques. The single-crystal X-ray diffraction analysis demonstrated that all the MOFs feature a similar type of three-dimensional structure having a binuclear [Cd2(COO)4(N)4] secondary building block unit. Moreover, MOFs 1 and 2 contain one-dimensional channels along the b-axis, whereas MOF 3 possesses a 1D channel along the a-axis. In these MOFs, the pores are decorated with multifunctional groups, i.e., halogen and amine. The gas adsorption analysis of these MOFs demonstrate that they display high uptake of CO2 (up to 5.34 mmol/g) over N2 and CH4. The isosteric heat of adsorption (Qst) value for CO2 at zero loadings is in the range of 18-26 kJ mol-1. In order to understand the mechanism behind the better adsorption of CO2 by our MOFs, we have also performed configurational bias Monte Carlo simulation studies, which confirm that the interaction between our MOFs and CO2 is stronger compared to those with N2 and CH4. Various noncovalent interactions, e.g., halogen (X)···O, Cd···O, and O···O, between CO2 and the halogen atom, the Cd(II) metal center, and the carboxylate group from the MOFs are observed, respectively, which may be a reason for the higher carbon dioxide adsorption. Ideal adsorbed solution theory (IAST) calculations of MOF 1 demonstrate that the obtained selectivity values for CO2/CH4 (50:50) and CO2/N2 (15:85) are ca. 28 and 193 at 273 K, respectively. However, upon increasing the temperature to 298 K, the selectivity value (S = 34) decreases significantly for the CO2/N2 mixture. We have also calculated the breakthrough analysis curves for all the MOFs using mixtures of CO2/CH4 (50:50) and CO2/N2 (50:50 and 15:85) at different entering gas velocities and observed larger retention times for CO2 in comparison with other gases, which also signifies the stronger interaction between our MOFs and CO2. Moreover, due to the presence of Lewis acidic metal centers, these MOFs act as heterogeneous catalysts for the CO2 fixation reactions with different epoxides in the presence of tetrabutyl ammonium bromide (TBAB), for conversion into industrially valuable cyclic carbonates. These MOFs exhibit a high conversion (96-99%) of epichlorohydrin (ECH) to the corresponding cyclic carbonate 4-(chloromethyl)-1,3-dioxolan-2-one after 12 h of reaction time at 1 bar of CO2 pressure, at 65 °C. The MOFs can be reused up to four cycles without compromising their structural integrity as well as without losing their activity significantly.

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.

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).

Spodium bonding with noble gas atoms

The nature of the bonding between a neutral group 12 member (Zn3, Cd3 and Hg3) ring and a noble gas atom was explored using quantum chemical simulations. Natural bond orbital, quantum theory of atoms in molecules, symmetry-adapted perturbation theory, and molecular electrostatic potential surface analysis were also used to investigate the type of interaction between the noble gas atom and the metal rings (Zn3, Cd3 and Hg3). The Zn3, Cd3 and Hg3 rings are bonded to the noble gas through non-covalent interactions, which was revealed by the non-covalent interaction index. Additionally, energy decomposition analysis reveals that dispersion energy is the key factor in stabilizing these systems.

Anisotropic Thermal Expansion of Structurally Related Lanthanide-Mercury(II) Cyanide Coordination Polymers

Three sets of related lanthanide-mercury(II) cyanide coordination polymers were synthesized by the reaction of LnCl3·xH2O (Ln = Ce, Nd, Eu, Gd, Tb, Dy, Ho, Tm, Yb, and Lu) with Hg(CN)2 and structurally characterized. [Ce(OH2)5][Hg(CN)2Cl]3·2H2O is a 3-D material with sheet-based architecture; its thermal expansion behavior shows uniaxial negative thermal expansion (-18.3(8), 39(2), and 68.3(16) ppm K-1 along the a, b, and c axes, respectively). This anisotropic thermal behavior is postulated to be driven elastically by weak Hg···Cl interactions: large area expansion of the sheets causes negative thermal expansion in the perpendicular direction. Using lanthanides heavier than Ce yielded 2-D sheet-based compounds with the formula [Ln(OH2)x]2[Hg(CN)2]5Cl6·2H2O (Ln = Nd and Eu, x = 7; Ln = Gd, Tb, Dy, Ho, Tm, Yb, and Lu, x = 6). Although there was also evidence for elastic behavior within these materials, both showed uniaxial zero thermal expansion (Ln = Nd: 27.9(17), 22.4(10), and 0.6(12) ppm K-1 along the I, II, and III principal axes, respectively; Ln = Tb: 39.6(12), 1.1(17), and 36.1(7) ppm K-1 along the a, b, and c axes, respectively). Despite their similar structural architecture, this zero thermal expansion was found to occur in different directions─within the plane of the 2-D sheets for [Nd(OH2)7]2[Hg(CN)2]5Cl6·2H2O but in the direction perpendicular to the 2-D sheets for [Tb(OH2)6]2[Hg(CN)2]5Cl6·2H2O. Overall, this system of compounds reveals the delicate relationship between coordination polymer structure and thermal expansion.

Spodium bonding in bis(alkynyl)mecurials

The new bis(alkynyl)mercurial Hg{CCSeCW(CO)2(Tp*)}2 (Tp* = tris(dimethylpyrazolyl)borate) forms adducts with fluoride and phenathroline, the structures of which are interpreted in the context of two-coordinate mercury presenting a σ-torroid for spodium bonding.

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.

Wolfium bonds in homodimers of MX4Y (M = Mo, W; X = F, Cl, Br; Y = O, S, Se)

The term "wolfium bond" has been recently introduced to describe the noncovalent attraction between an atom of group 6 and a nucleophile via a σ-hole binding site. Crystal structures commonly contain a motif wherein two MX4Y units are arranged in close proximity, where M represents either Mo or W, and X and Y refer to halogen and chalcogen atoms respectively. DFT calculations were thus applied to a wide range of homodimers of these molecules so as to assess their preferred arrangements, and to characterize the types of bonding that are present in each in a systematic manner. The most stable Dual-X configuration is symmetric and contains a pair of equivalent M⋯X bonds. The interaction energies range from -8 to -29 kcal mol-1, and are largest for X = F, Y = O, and M = W. The X electron donor is replaced by Y, and the two wolfium bonds are reduced to one, in the less stable Mono-Y structure, with interaction energies between -2 and -10 kcal mol-1. There is some question as to whether the weaker bonds of this type constitute true wolfium bonds.

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.

Experimental and Theoretical Survey of Intramolecular Spodium Bonds/σ/π-Holes and Noncovalent Interactions in Trinuclear Zn(II)-Salen Type Complex with OCN- Ions: A Holistic View in Crystal Engineering

In this work, one new centrosymmetric trinuclear Zn(II) complex 1, [{(OCN)Zn(L)}2Zn], using a salen-type ligand (H2L) in the presence of OCN- was synthesized and characterized via elemental, spectral, SEM-EDX, and single-crystal X-ray diffraction (SCXRD) study. In 1, SCXRD reveals two different stereochemical environments of zinc metal ions; one terminal Zn(II) center adopts square pyramidal geometries with the Addison parameter (τ) 0.095, and the central Zn(II) is tetracoordinated tetrahedral geometry. This article provides evidence of the significance and presence of spodium bonds (SpBs) in solid-state crystal structures involving a pseudotetrahedral environment of the central Zn-atom. X-ray structures reveal intramolecular Zn···O SpBs caused by the methoxy (-OCH3) substituent O-atom adjacent to the coordinated phenoxy O-atom. These noncovalent interactions have been thoroughly studied using density functional theory calculations at the RI-BP86[2]-D3[3]/def2-TZVP level of theory that characterizes the nature of SpBs, including the Baders quantum theory of atoms-in-molecules "QTAIM", molecular electrostatic potential (MEP) surface, and noncovalent index plot (NCI). In addition, a unique complex-isomer-based theoretical model has been vividly employed to estimate the SpBs energy in the complex. Natural bond orbital (NBO) analysis also tries to establish the differentiation between σ-hole and π-hole SpBs' natures more authentically. The complex energy frameworks were used to investigate noncovalent interactions. To better understand the different intermolecular interactions, we conducted a Hirshfeld surface, which revealed N···H (15.4%) and O···H (9.1%) contacts and Zn···O (5.1%) (SpBs).

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.

Biological noncovalent N/O⋯V interactions: insights from theory and protein data bank analyses

Computations at the PBE0-D3/def2-TZVP level of theory in conjunction with a Protein Data Bank (PDB) survey have provided first time evidence of favorable noncovalent interactions between ADP metavanadate (VO4) and ADP orthovanadate (VO5) and electron rich atoms. These involve a σ-hole present in the V atom and the lone pairs belonging to (i) protein residues (e.g., serine (SER), glutamate (GLU) or histidine (HIS)), (ii) backbone carbonyl groups and (iii) water molecules. A computational study has been carried out to rationalize the physical nature and directionality of the interaction in addition to its plausible biological role. The results reported herein are expected to have an impact in the fields of medicinal chemistry, bioinorganic chemistry and chemical biology.

Gold(III) derivatives as the noncovalent interaction donors: theoretical study of the π-hole regium bonds

Except for the well-known σ-hole regium bonds formed by metal nanoparticles and M(I) (M = Cu, Ag, and Au) derivatives, the existence of π-hole regions located above and below the Au atom in gold(III) derivatives suggests that gold(III) also functions as an efficient electrophilic site. In this study, a comprehensive analysis was conducted on the electrophilicity of trichloro-(p-toluonitrilo-N)-gold(III) derivatives AuL3(NCC6H4X) (L = Cl, Br, CN; X = NH2, CH3, CF3, NC, and CN) and the nature of π-hole regium bonds in the AuL3(NCC6H4X)⋯LB (LB = NH3, N(NH3)3, CH2O, C2H2, C2H4, C6H6) and (AuCl3(NCC6H4Y))n (Y = Cl, CN, NC, NO2; n = 2, 3)) complexes. The characteristics of the π-hole regium bonds were studied with respect to the influence of ligands and substituents, the strength of intermolecular interactions between Au(III) derivatives and Lewis bases, and those in the polymers. In the case of the AuL3(NCC6H4X)⋯NH3 complexes, the strength of the regium bonds increases gradually in the order of L = Cl < Br < CN and X = NH2 < CH3 < CF3 ≈ NC < CN. The ligands (L) attached to the Au atom exert a significant effect on the strength of the π-hole regium bonds in comparison to the substituents (X) on the benzene ring. The regium bonds are primarily dominated by electrostatic interaction, accompanied by moderate contribution from polarization. Linear relationships were identified between the electrostatic energies and the local most positive potentials over the Au atom, as well as between the polarization energies and the amount of charge transfer. Most of the π-hole regium bonds in the AuL3(NCC6H4X)⋯LB complexes exhibit the characters of closed shell noncovalent interactions. In the polymers (AuCl3(NCC6H4Y))n, weak face-to-face π-π stacking interactions are also present, in addition to regium bonds. The trimers displayed a slightly negative cooperativity in comparison to the dimers.

Erythronium Bonds: Noncovalent Interactions Involving Group 5 Elements as Electron-Density Acceptors

Analyses of the Cambridge Structural Database and theoretical calculations (PBE0-D3/def2-TZVP level, atoms-in-molecules, natural bond orbital studies) prove the formation of net attractive noncovalent interactions between group 5 elements and electron-rich atoms (neutral or anionic). These kinds of bonding are markedly different from coordination bonds formed by the same elements and possess the distinctive features of σ-hole interactions. The term erythronium bond is proposed to denote these bonds. X-ray structures of vanadate-dependent bromoperoxidases show that these interactions are present also in biological systems.

Ligand and Substituent Effect on Regium-π Bonding in Cu and Ag π-Conjugated Complexes: A Density Functional Study

Density functional theory investigation of regium (Rg)-π bonding using the RgL-X model system, where Rg = Cu and Ag; L = CN, NO2, and OH; X = π-conjugated system (benzene, cyanobenzene, benzoic acid, pyridine, 2-methoxy aniline, 1,4-dimethoxy benzene, and cyclophane), has been performed. Conclusive evidence of the Rg-π bond has been provided by analysis of molecular electrostatic potential surfaces, Rg-π bond length, interaction energy (ΔE), second-order perturbation energy (E2), charge transfer (Δq), quantum theory of atom in molecules, and noncovalent interaction plots for 42 structural arrangements with varying ligands and the substituted aromatic ring. The Rg-π bond length in the optimized model systems varies from 2.03 to 2.12 Å in Cu complexes (1-21) and from 2.26 to 2.38 Å in Ag complexes (22-42) at the PBE0-D3 functional. While the ligand (L) attached to the Rg metal has a bargaining effect on the strength of the Rg-π bond (in the order of -OH > -CN = -NO2), the π-conjugated systems have a diminutive effect. Two X-ray crystal structures (CUCSOI and AHIDQU) having the Rg-π bond, accessed from Cambridge Crystallographic Data Centre (CCDC), are discussed here to signify the influence of Rg-π bonding on the crystal structure.

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.

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.

TFRegNCI: Interpretable Noncovalent Interaction Correction Multimodal Based on Transformer Encoder Fusion

The interpretability is an important issue for end-to-end learning models. Motivated by computer vision algorithms, an interpretable noncovalent interaction (NCI) correction multimodal (TFRegNCI) is proposed for NCI prediction. TFRegNCI is based on RegNet feature extraction and a transformer encoder fusion strategy. RegNet is a network design paradigm that mainly focuses on local features. Meanwhile, the Vision Transformer is also leveraged for feature extraction, because it can capture global features better than RegNet while lowering the computational cost. Using a transformer encoder as the fusion strategy rather than multilayer perceptron can enhance model performance, due to its emphasis on important features with less parameters. Therefore, the proposed TFRegNCI achieved high accurate prediction (mean absolute error of ∼0.1 kcal/mol) comparing with the coupled cluster single double (triple) (CCSD(T)) benchmark. To further improve the model efficiency, TFRegNCI applies two-dimensional (2D) inputs transformed from three-dimensional (3D) electron density cubes, which saves time (30%), while the model accuracy remains. To improve model interpretability, a visualization module, Gradient-weighted Regression Activation Mapping (Grad-RAM) has been embedded. Grad-RAM is promoted from the classification algorithm, Gradient-weighted Class Activation Mapping, to perform feature visualization for the regression task. With Grad-RAM, the visual location map for features in deep learning models can be displayed. The feature map visualizations suggest that the 2D model has the similar performance as the 3D model, because of equally effective feature extractions from electron density. Moreover, the valid feature region on the location map by the 3D model is consistent with the NCIPLOT NCI isosurface. It is confirmed that the model does extract significant features related to the NCI interaction. The interpretable analyses are carried out through molecular orbital contribution on effective features. Thereby, the proposed model is likely to be a promising tool to reveal some essential information on NCIs, with regard to the level of electronic theory.

Beryllium bonding with noble gas atoms

Quantum chemical calculations were carried out to investigate the nature of the bonding between a neutral Be₃ ring and noble gas atom. Electronic structure calculation for these complexes was carried out at different computational levels in association with natural bond orbital, quantum theory of atoms in molecules, electron localization function, symmetry adapted perturbation theory, and molecular electrostatic potential surface analysis of Be₃ complexes. The Be atoms in the Be₃ moiety are chemically bonded to one another, with the BeBe bond dissociation energy being ~125 kJ mol^-1 . The Be₃ ring interacts with the noble gases through non-covalent interactions. The binding energies of the noble gas atoms with the Be₃ ring increases with increase in their atomic number. The non-covalent interaction index, density overlap region indicator and independent gradient model analyses reveal the presence of non-covalent inter-fragment interactions in the complexes. Energy decomposition analysis reveals that dispersion plays the major role towards stabilizing these systems.

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.

Anion···Anion [AuI4]-···[AuI2]- Complex Trapped in the Solid State by Tetramethylammonium Cations

A discrete π-hole···σ-hole dimer is synthesized and X-ray characterized. It presents a perfect thumbtack geometry where the σ-hole of the linear [AuI₂]^- anion points to the π-hole located above the central Au-atom of the [AuI₄]^- anion. Such discrete π-hole···σ-hole dimers are unprecedented in literature, since all mixed-valence gold(I/III) iodide compounds reported to date form infinite ···([AuI₄]^-···[AuI₂]^-) _n ·· chains in the solid state. If an excess of iodine is used for the synthesis, triiodide [I₃]^- ions are partially incorporated into the [AuI₂]^- sites, forming infinite chains. The nature of the anion···anion interaction has been studied considering two possibilities: (i) a π-hole coinage bond or (ii) σ-hole halogen bond using high-level density functional theory calculations, the quantum theory of atoms in molecules, and the noncovalent interaction plot index.

Matere Bonds vs. Multivalent Halogen and Chalcogen Bonds: Three Case Studies

The term matere bond has been recently used to refer to an attractive noncovalent interaction between any element of group 7 acting as an electrophile and any atom (or group of atoms) acting as a nucleophile. The utilization of metals such as σ-hole donors is starting to attract the attention of the scientific community. In this manuscript, a comparison between matere bonds and well-known σ-hole interactions (halogen and chalcogen bonds) is carried out using three X-ray structures, retrieved from the Cambridge structural database (CSD), and density functional theory calculations (DFT). The novelty of this work resides in the utilization of a neutral Re(VII) system as the matere bond donor and multivalent chalcogen and halogen donors. In fact, as far as our knowledge extends, the description of σ-hole interactions in Se(VI) is unprecedented in the literature. The σ-hole interactions in Re(VII), Se(VI) and Cl(VII) electron acceptors are analyzed and compared using several computational tools.

Crystallographic and Theoretical Study of Osme Bonds in Nitrido-Osmium(VI) Complexes

Osme bonds have been recently defined as the attractive interaction between an element of group 8 acting as an electrophile and any atom or group of atoms acting as a nucleophile. To date, the known examples of osme bonds in X-ray structures involve mostly the highly reactive OsO4 and amines and amine oxides. In this work, evidence supporting the existence of osme bonds in osmium(VI) derivatives is reported. In particular, nitrido-osmium(VI) complexes that present square-pyramidal geometries are well disposed to participate in osme bonds opposite to the Os≡N bond. By using a combination of experimental and theoretical results, the existence and importance of this new class of σ-hole interactions is demonstrated in the solid state of several nitrido-osmium(VI) derivatives.

Recognition in the Domain of Molecular Chirality: From Noncovalent Interactions to Separation of Enantiomers

It is not a coincidence that both chirality and noncovalent interactions are ubiquitous in nature and synthetic molecular systems. Noncovalent interactivity between chiral molecules underlies enantioselective recognition as a fundamental phenomenon regulating life and human activities. Thus, noncovalent interactions represent the narrative thread of a fascinating story which goes across several disciplines of medical, chemical, physical, biological, and other natural sciences. This review has been conceived with the awareness that a modern attitude toward molecular chirality and its consequences needs to be founded on multidisciplinary approaches to disclose the molecular basis of essential enantioselective phenomena in the domain of chemical, physical, and life sciences. With the primary aim of discussing this topic in an integrated way, a comprehensive pool of rational and systematic multidisciplinary information is provided, which concerns the fundamentals of chirality, a description of noncovalent interactions, and their implications in enantioselective processes occurring in different contexts. A specific focus is devoted to enantioselection in chromatography and electromigration techniques because of their unique feature as "multistep" processes. A second motivation for writing this review is to make a clear statement about the state of the art, the tools we have at our disposal, and what is still missing to fully understand the mechanisms underlying enantioselective recognition.

Local Charge-Displacement Analysis: Targeting Local Charge-Flows in Complex Intermolecular Interactions

Charge-displacement (CD) analysis has recently proven to be a simple and powerful scheme for quantitatively analyzing the profile of the charge redistribution occurring upon intermolecular interactions along a given interaction axis. However, when two molecular fragments bind through complex interactions involving multiple concurrent charge flows, ordinary CD analysis is capable of providing only an averaged picture of the related charge-flow profiles and no detailed information on each of them. In this article, we combine CD analysis with a Hirshfeld partitioning of the molecular charge redistribution for a local analysis on focused portions of the molecule, allowing for a detailed characterization of one charge flow at a time. The resulting scheme - the local charge-displacement (LCD) analysis - is tested on the intriguing case of the dimethyl sulfide (DMS)-sulfur dioxide (SO2) complex, characterized by concurrent charge flows relating to a sulfur-sulfur homo-chalcogen interaction and a pair of hydrogen bonds. The LCD scheme is then applied to the analysis of multiple hydrogen bonding in the acetic acid dimer, of base-pairing interactions in DNA, and of ambifunctional hydrogen bonding in the ammonia-pyridine complex.

Periodate anions as a halogen bond donor: formation of anion⋯anion dimers and other adducts

Single crystal X-ray analyses show that iodine in pyridinium periodates acts as a halogen bond (HaB) donor forming short and almost linear contacts with neutral and anionic electron donors. A combination of QTAIM and NCIplot computational tools proves the attractive nature of these contacts.

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.

Noncovalent Interactions Involving Group 6 in Biological Systems: The Case of Molybdopterin and Tungstopterin Cofactors

In this study we propose to coin the term Wolfium bond (WfB) to refer to a net attractive force (noncovalent interaction) between any element of group 6 and electron donor atoms (neutral molecules or anions) and to differentiate it from a coordination bond (metal‐ligand interaction). We provide evidence of the existence of this interaction by inspecting the X‐ray crystal structure of proteins containing Molybdopterin and Tungstopterin cofactors from the Protein Data Bank (PDB). The plausible biological role of the interaction as well as its physical nature (antibonding Wf‐Ligand orbital involved) are also analyzed by means of ab initio calculations (RI‐MP2/def2‐TZVP level of theory), Atoms in Molecules (AIM), Natural Bond Orbital (NBO) and Noncovalent Interactions plot (NCIplot) analyses.

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.

Insight into Spodium–π Bonding Characteristics of the MX2···π (M = Zn, Cd and Hg; X = Cl, Br and I) Complexes—A Theoretical Study

The spodium–π bonding between MX₂ (M = Zn, Cd, and Hg; X = Cl, Br, and I) acting as a Lewis acid, and C₂H₂/C₂H₄ acting as a Lewis base was studied by ab initio calculations. Two types of structures of cross (T) and parallel (P) forms are obtained. For the T form, the X–M–X axis adopts a cross configuration with the molecular axis of C≡C or C=C, but both of them are parallel in the P form. NCI, AIM, and electron density shifts analyses further, indicating that the spodium–π bonding exists in the binary complexes. Spodium–π bonding exhibits a partially covalent nature characterized with a negative energy density and large interaction energy. With the increase of electronegativity of the substituents on the Lewis acid or its decrease in the Lewis base, the interaction energies increase and vice versa. The spodium–π interaction is dominated by electrostatic interaction in most complexes, whereas dispersion and electrostatic energies are responsible for the stability of the MX₂⋯C₂F₂ complexes. The spodium–π bonding further complements the concept of the spodium bond and provides a wider range of research on the adjustment of the strength of spodium bond.

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.

HgBrI: a possible tecton for NLO molecular materials?

Mixed mercury(II) halogenides have been known for a long time as good NLO (non-linear optic) materials. The NLO properties are due to the halogen disposition in the solid state and the electron distribution among the bonds formed by soft elements. We investigated the possibility of using HgBrI as a asymmetric tecton in the preparation of noncentrosymmetric crystalline compounds, by exploiting the coordinating power of Hg(II) toward N-donor ligands, and seven coordination complexes have been obtained. To unravel the nature of these complex systems we combined the data from different techniques: Raman spectroscopy, SC-XRD and Second Harmonic Generation, supported by a periodic DFT computational approach. In HgBrI crystalline products with low symmetry, the presence of substitutional disorder leads to a lack of the inversion center conferring NLO activity, which is absent in analogous complexes of Hg(II) halogenides. These results indicate HgBrI as an interesting tecton to obtain metallorganic NLO materials.

Substituent effects in π-hole regium bonding interactions between Au(p-X-Py)2 complexes and Lewis bases: an ab initio study

For the first time, long range substituent effects in regium bonding interactions involving Au(I) linear complexes are investigated. The Au(I) atom is coordinated to two para -substituted pyridine ligands. The interaction energy (RI-MP2/def2-TZVP level of theory) of the π-hole regium bonding assemblies is affected by the pyridine substitution. The Hammett's plot representations for several sets of Lewis bases have been carried out and, in all cases, good regression plots have been obtained (interaction energies vs. Hammett's σ parameter). The Bader's theory of "atoms-in-molecules" has been used to evidence that the electron density computed at the bond critical point that connects the Au-atom to the electron donor can be used as a measure of bond order in regium bonding. Several X-ray structures retrieved from the Cambridge Structural Database (CSD) provide some experimental support to the existence of regium π-hole bonding in [Au(Py) 2 ] + derivatives.

Synthesis, spectroscopic findings and crystal engineering of Pb(ii)-Salen coordination polymers, and supramolecular architectures engineered by σ-hole/spodium/tetrel bonds: a combined experimental and theoretical investigation

Spontaneous self-assembly is one of the available synthetic routes to achieve structurally versatile and unique crystal complexes with selected metal-ligand combinations in the spirit of pseudohalides. In this endeavour, we designed a novel 1D coordination polymer (CP), [(Cd)(Pb)(L)(η¹-NCS)(η¹-SCN)] _n (1), using a compartmental Salen ligand (H₃L) in the presence of NaSCN. The characterization of the CP was accomplished using several spectroscopic techniques: MALDI-TOF, PXRD, SEM, EDX mapping, and single-crystal X-ray crystallography. The CP crystallizes in the monoclinic space group P2₁/c with Z = 4. SCXRD reveals Cd(ii) and Pb(ii) metal ions fulfilled distorted square pyramidal and hemi-directed coordination spheres. Cd(ii) is placed in the inner N₂O₂ and heavier Pb(ii) in the outer O₄ compartments of the de-protonated form of the ligand [L]^2-. Supramolecular interactions in the intricate crystal structure produced attractive molecular architectures of the compound. The flexible aliphatic -OH pendent group coordinates with the Pb(ii) ions. This unique binding further elevates the supramolecular crystal topographies. The supramolecular interactions were authenticated by Hirshfeld surface analysis (HSA). The observation of the recurring unconventional tetrel bonds was rationalized by DFT calculations and surface plots of molecular electrostatic potential (MEP). In the 1D polymeric chain in the complex, the O-atom of the -OH groups shows a tetrel bonding interaction with the Pb atom. We have found that the combination of QTAIM/NCI and QTAIM/ELF plots helps reveal the nature of these contacts. Moreover, the QTAIM/ELF plot determines the donor-acceptor interaction between the O-atom and the Pb atom, establishing the σ-hole. Agreeably, the σ-hole interaction also helps Pb(ii) serve as a Lewis acid in the complex. Finally, spodium and tetrel bonds are formed, possible thanks to a hemi-directional coordination sphere of the Pb atoms in the polymer described.

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.

Arsolyl-Supported Intermetallic Dative Bonding

The first examples of late transition metal η⁵-arsolyls (L = CO, P(OMe)₃; R = Ph, Me, Et, SiMe₃; R′ = Ph, H, Me, Et, Me) serve as ditopic donors to extraneous metal centres (M = Pt^II, Au^I, Hg^II) through both conventional As → M and polar-covalent (dative) Co → M interactions.

Cobalt carbonyl reacts with arsoles to provide the first late transition metal η⁵-arsolyls. These serve as ditopic donors to extraneous metal centres (M = Pt^II, Au^I, Hg^II) through both conventional AsM and polar-covalent (dative) CoM interactions.

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.

Inter-anion chalcogen bonds: Are they anti-electrostatic in nature?

Inter-anion hydrogen and halogen bonds have emerged as counterintuitive linkers and inspired us to expand the range of this unconventional bonding pattern. Here, the inter-anion chalcogen bond (IAChB) was proposed and theoretically analyzed in a series of complexes formed by negatively charged bidentate chalcogen bond donors with chloride anions. The kinetic stability of IAChB was evidenced by the minima on binding energy profiles and further supported by ab initio molecular dynamic simulations. The block-localized wave function (BLW) method and its subsequent energy decomposition (BLW-ED) approach were employed to elucidate the physical origin of IAChB. While all other energy components vary monotonically as anions get together, the electrostatic interaction behaves exceptionally as it experiences a Coulombic repulsion barrier. Before reaching the barrier, the electrostatic repulsion increases with the shortening Ch⋯Cl^- distance as expected from classical electrostatics. However, after passing the barrier, the electrostatic repulsion decreases with the Ch⋯Cl^- distance shortening and subsequently turns into the most favorable trend among all energy terms at short ranges, representing a dominating force for the kinetic stability of inter-anions. For comparison, all energy components exhibit the same trends and vary monotonically in the conventional counterparts where donors are neutral. By comparing inter-anions and their conventional counterparts, we found that only the electrostatic energy term is affected by the extra negative charge. Remarkably, the distinctive (nonmonotonic) electrostatic energy profiles were reproduced using quantum mechanical-based atomic multipoles, suggesting that the crucial electrostatic interaction in IAChB can be rationalized within the classical electrostatic theory just like conventional non-covalent interactions.

On the Importance of Pnictogen and Chalcogen Bonding Interactions in Supramolecular Catalysis

In this review, several examples of the application of pnictogen (Pn) (group 15) and chalcogen (Ch) bonding (group 16) interactions in organocatalytic processes are gathered, backed up with Molecular Electrostatic Potential surfaces of model systems. Despite the fact that the use of catalysts based on pnictogen and chalcogen bonding interactions is taking its first steps, it should be considered and used by the scientific community as a novel, promising tool in the field of organocatalysis.