Not Only Hydrogen Bonds: Other Noncovalent Interactions

Weak Interaction Activates Esters: Reconciling Catalytic Activity and Turnover Contradiction by Tailored Chalcogen Bonding

The activation of esters by strong Lewis acids via the formation of covalent adducts is a classic strategy to give reactivity; however, this approach frequently incurs limited turnover due to the low efficiency in the dissociation of catalyst from a stable catalyst-product complex. While the use of some weak interaction catalysts that can easily dissociate from any bonding complexes in the reaction system would solve this catalyst turnover problem, the poor catalytic activity in the ester activation that can be provided by these noncovalent forces in turn sets up a formidable challenge. Herein, we describe the activation and catalytic transformation of esters by weak interactions, which provides a promising platform to reconcile the catalytic activity and turnover problems. Several tailored chalcogen-bonding catalysts were developed for the activation of esters, enabling achieving several inherently low reactive Diels-Alder reactions as well as the ring-opening polymerization of lactones through weak chalcogen bonding interactions. This supramolecular catalysis approach is particularly highlighted by its capability to promote some uncommon Diels-Alder reactions involving using dienes bearing electron-withdrawing groups coupled by α,β-unsaturated ester as dienophiles and substrate incorporating competitive Lewis basic sites, in which typical strong Lewis acids showed low catalytic efficiency, while representative hydrogen and halogen bonding catalysts were inactive.

Cationic Hypervalent Chalcogen Bond Catalysis on the Povarov Reaction: Reactivity and Stereoselectivity

Chalcogen bond catalysis, particularly cationic hypervalent chalcogen bond catalysis, is considered to be an effective strategy for organocatalysis. In this work, the cationic hypervalent chalcogen bond catalysis for the Povarov reaction between N-benzylideneaniline and ethyl vinyl ether was investigated by density functional theory (DFT). The catalytic reaction involves the cycloaddition process and the proton transfer process, and the rate-determining step is the cycloaddition process. Cationic hypervalent tellurium derivatives bearing CF3 and F groups exhibit superior catalytic activity. For the rate-determining step, the Gibbs free energy barrier decreases as the positive electrostatic potential of the chalcogen bond catalysts increases. More importantly, the Gibbs free energy barrier has a strong linear correlation with the electrostatic energy of the chalcogen bond in the catalyst-substrate complex. Furthermore, the catalytic reactions include the endo pathway and exo pathway. The C-H⋅⋅⋅π interaction between the substituent of the ethyl vinyl ether and the aryl ring of the N-benzylideneaniline contributes to the endo-selectivity of the reaction. This research contributes to a deeper understanding of chalcogen bond catalysis, providing insights for designing chalcogen bond catalysts with high performance.

Cooperativity between Intermolecular Hydrogen and Carbon Bonds in ZY···CH3CN/CH3NC···HX Trimers (ZY = H2O, H2S, HF, HCl, HBr, NH3, and H2CO; HX = HF, HCl, and HBr)

Hydrogen-bonding and carbon-bonding interactions are widespread in nature. We studied the cooperativity between these interactions in 42 trimeric complexes ZY···CH3CN/CH3NC···HX, where ZY molecules are H2O, H2S, HF, HCl, HBr, NH3, and H2CO, and HX molecules are HF, HCl, and HBr. Acetonitrile (CH3CN) and isoacetonitrile (CH3NC) act as hydrogen bond acceptors as well as carbon bond donors in these trimers. Various theoretical methods, such as electronic structure calculations, quantum theory of atoms in molecule (QTAIM), natural bond orbital (NBO), and reduced density gradient analysis, are employed to study these trimers, and the results are compared with the corresponding ZY···CH3CN/CH3NC and CH3CN/CH3NC···HX dimers. Electronic structure calculations are performed at the second-order Mo̷ller-Plesset perturbation theory using the 6-311++G(2d,2p) basis set. We show that both the interactions act synergistically in these trimers leading to an increase in their bond strength as compared to the strength in the individual dimers. The cooperative energies for these trimers are in the range of 0.69 to 3.22 kJ/mol. It is seen that the carbon bonds benefit more from the cooperativity than the hydrogen bonds. The trends of cooperativity and correlations of interaction energies and cooperative energies with relevant QTAIM and NBO parameters are reported.

Comparison of N···I and N···O Halogen Bonds in Organoiodine Cocrystals of Heterocyclic Aromatic Diazine Mono-N-oxides

A series of cocrystals of halogen bond donors 1,4-diiodotetrafluorobenzene (p-F4DIB) and tetraiodoethylene (TIE) with five aromatic heterocyclic diazine mono-N-oxides based on pyrazine, tetramethylpyrazine, quinoxaline, phenazine, and pyrimidine as halogen bonding acceptors were studied. Structural analysis of the mono-N-oxides allows comparison of the competitive occurrence of N···I vs O···I interactions and the relative strength and directionality of these two types of interactions. Of the aromatic heterocyclic diazine mono-N-oxide organoiodine cocrystals examined, six exhibited 1:1 stoichiometry, forming chains that utilized both N···I and O···I interactions. Two cocrystals presented 1:1 stoichiometry with exclusive O···I interactions. Two cocrystals displayed a 2:1 stoichiometry-one characterized solely by O···I interactions and the other solely by N···I interactions. We have also compared these interactions to those present in the corresponding diazines, some of which we report here and some which have been previously reported. In addition, a computational analysis using density functional theory (M062X/def2-SVPD) was performed on these two systems and has been compared to the experimental results. The calculated complex formation energies were, on average, 4.7 kJ/mol lower for the I···O halogen bonding interaction as compared to the corresponding N···I interaction. The average I···O interaction distances were calculated to be 0.15 Å shorter than the corresponding I···N interactions.

Computational tools to study non-covalent interactions and confinement effects in chemical systems

Confinement is a very common phenomenon in chemistry, for example, when molecules are located inside cavities. In these conditions, the electronic structure of atoms and molecules is modified. These changes could be mapped through the interaction with other molecules since non-covalent interactions between molecules are also influenced by confinement. In this work we address both topics, non-covalent interactions, and confined systems, using quantum chemistry tools with new software, emphasizing the importance of analyzing both fields simultaneously.

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.

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

Triiodide-Based Chair-Like Copper Complex Assembled by Halogen Bonding

Cocrystallization of the dimeric [Cu2(μ-I)2(CNXyl)4] (Xyl = 2,6-Me2C6H3, 1) and polymeric catena-[Cu(μ-I)(CNC6H3-2-Cl-6-Me)2] (2) complexes with I2 at different molar ratios between the reactants resulted in a series of (RNC)2CuI-based crystal polyiodides formed along with gradual accumulation of iodine, namely the cocrystals [1·I2]·[Cu(μ1,1-I3)(CNXyl)2]2 followed by the generation of [Cu(μ1,3-I3)(CNXyl)2]2·2I2 (5·2I2) or [Cu(μ1,1-I3)(CNC6H3-2-Cl-6-Me)2]2 and then [Cu(μ1,3-I3)(CNC6H3-2-Cl-6-Me)2]n·n/2I2. The polyiodide 5·2I2 exhibits a novel supramolecular motif─a purely inorganic halogen-bonded Cu2(μ1,3-I3)2 core in the chair conformation. The X-ray structure of 5·2I2 featuring I···I contacts was analyzed by a set of theoretical methods and attributed to moderately strong halogen bonding (from -3.2 to -3.9 kcal/mol); these interactions determine the supramolecular architecture of 5·2I2.

Exploring the electron donor-acceptor duality of B3N3 in noncovalent interactions

Boron nitrides are very important and are used as lubricants, insulating agents, etc. Interactions of such systems with small molecules are important. This study examined the potential of B3N3 (triboron trinitride) to act as both an electron acceptor and an electron donor in the formation of noncovalent interactions. The anisotropic electronic distribution observed in the electrostatic potential map supported the B3N3's ability to exhibit the predicted electron donor-acceptor duality. Further computational investigations on optimized gas-phase complexes of B3N3:(NH3)n=1-3, B3N3:(NCH)n=1-6, B3N3:(N2H2)n=1-3 and (B3N3)2 confirmed that the B3N3 molecule can participate in B⋯N triel bonding interactions and H···N hydrogen bonding interactions. These energetically stable complexes are primarily governed by electrostatic and polarization interactions.

Interplay of Hydrogen, Pnicogen, and Chalcogen Bonding in X(H2O)n=1-5 (X = NO, NO+, and NO-) Complexes: Energetics Insights via a Molecular Tailoring Approach

Nitric oxide (NO) and its redox congeners (NO+ and NO-), designated as X, play vital roles in various atmospheric and biological events. Understanding the interaction between X and water is inevitable to explain the different reactions that occur during these events. The present study is a unified attempt to explore the noncovalent interactions in microhydrated networks of X using the MP2/aug-cc-pVTZ//MP2/6-311++G(d,p) level of theory. The interactions between X and water have been probed by the molecular electrostatic potential (MESP) by exploiting the features of the most positive (Vmax) and most negative potential (Vmin) sites. The individual energy and cooperativity contributions of various types of noncovalent interactions present in X(H2O)n=1-5 complexes are estimated with the help of a molecular tailoring-based approach (MTA-based). The MTA-based analysis reveals that among various possible interactions in NO(H2O)n complexes, the water···water hydrogen bonds (HBs) are the strongest. Neutral NO can form hydrogen and pnicogen bonds (PBs) with water depending on the orientation; however, such HBs and PBs are the weakest. On the other hand, in the NO+(H2O)n complexes, the NO+···water interactions that occur through PBs are the strongest; the next one is the chalcogen bonding (CB), and the water···water HBs are the weakest. In the case of the NO-(H2O)n complexes, the HB interactions via both N and O atoms of NO- and water molecules are the strongest ones. The strength of water···water HB interactions is also seen to increase with the increase in the number of water molecules in NO-(H2O)n. The present study exemplifies the applicability of MTA-based calculations for quantifying various types of individual noncovalent interactions and their interplay in microhydrated networks of NO and its related ions.

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.

Color properties and non-covalent interactions in hydrated (Z)-4-(1-cyano-2-(2,4,5-trimethoxyphenyl)-vinyl)pyridin-1-ium chloride salt: Insights from experimental and theoretical studies

The optical charge-transfer (CT) property and the crystal structure of (Z)-4-(1-cyano-2-(2,4,5-trimethoxyphenyl)vinyl)pyridin-1-ium chloride monohydrate salt (I), which belongs to an acrylonitrile family, was studied. The title salt, I, was characterized using different spectroscopy techniques and a single-crystal X-ray diffraction study combined with quantum chemical computations. The results showed that the color properties of I are determined by the CT, changes in bandgap, optical absorption, and various non-covalent interactions. The HOMO-LUMO energy gaps are 5.41 eV and 5.23 eV for the precursor and salt, respectively. It was demonstrated that π-π stacking interactions lead to the formation of intercalated dimers and donor-acceptor interactions assisted by hydrogen bonds; the dimers and interactions are different between the precursor and the salt. The cation moiety is mainly stabilized by N(1)+-H···Cl, and the anion is predominantly stabilized by strong O(1W)- H⋯ Cl- bonds as well as the hydrogen bonds with the MeO group O(2W)-H⋯O(1) and O(2W)-H⋯O(1W). The charge transfer between cation and anion moieties in the structure is established through NBO analysis.

Revisiting conventional noncovalent interactions towards a complete understanding: from tetrel to pnicogen, chalcogen, and halogen bond

Typical noncovalent interactions, including tetrel (TtB), pnicogen (PniB), chalcogen (ChalB), and halogen bonds (HalB), were systematically re-investigated by modeling the N⋯Z interactions (Z = Si, P, S, Cl) between NH3 - as a nucleophilic, and SiF4, PF3, SF2, and ClF - as electrophilic components, employing highly reliable ab initio methods. The characteristics of N⋯Z interactions when Z goes from Si to Cl, were examined through their changes in stability, vibrational spectroscopy, electron density, and natural orbital analyses. The binding energies of these complexes at CCSD(T)/CBS indicate that NH3 tends to hold tightly most with ClF (-34.7 kJ mol-1) and SiF4 (-23.7 kJ mol-1) to form N⋯Cl HalB and N⋯Si TtB, respectively. Remarkably, the interaction energies obtained from various approaches imply that the strength of these noncovalent interactions follows the order: N⋯Si TtB > N⋯Cl HalB > N⋯S ChalB > N⋯P PniB, that differs the order of their corresponding complex stability. The conventional N⋯Z noncovalent interactions are characterized by the local vibrational frequencies of 351, 126, 167, and 261 cm-1 for TtB, PniB, ChalB, and HalB, respectively. The SAPT2+(3)dMP2 calculations demonstrate that the primary force controlling their strength retains the electrostatic term. Accompanied by the stronger strength of N⋯Si TtB and N⋯Cl HalB, the AIM and NBO results state that they are partly covalent in nature with amounts of 18.57% and 27.53%, respectively. Among various analysis approaches, the force constant of the local N⋯Z stretching vibration is shown to be most accurate in describing the noncovalent interactions.

Activation of alkynes by chalcogen bonding: a Se⋯π interaction catalyzed intramolecular cyclization of 1,6-diynes

The activation of the triple bond of alkynes was dominated by transition metals, while it is difficult for organocatalysts to play an effective role in this realm. Herein, we describe the activation of alkynes by chalcogen bonding, and the weak Se⋯π interaction was capable of catalyzing the intramolecular cyclization of 1,6-diynes, thus adding a new capability in the list of supramolecular catalysis.

Supramolecular catalysis with ethers enabled by dual chalcogen bonding activation

The activation of ethers by weak interactions is a long-standing objective in supramolecular catalysis, but yet it remains an underdeveloped topic. The obstacles towards solving this problem are prominent since it is difficult for a weak interaction to cleave a relatively strong C-O σ-bond and moreover, the ionic intermediate composing of an alkoxide ion and an electrophilic carbocation would deactivate weak interaction donors. Herein, we describe a distinctive activation mode, dual Se···π and Se···O bonding, that could activate benzylic as well as allylic ether C-O σ-bonds to achieve cyclization, coupling and elimination reactions. This dual Se···π and Se···O bonding catalysis approach could tolerate various alkoxide leaving groups, while the other representative weak interaction donors showed no catalytic activity.

Metastable Charged Dimers in Organometallic Species: A Look into Hydrogen Bonding between Metallocene Derivatives

Multiply charged complexes bound by noncovalent interactions have been previously described in the literature, although they were mostly focused on organic and main group inorganic systems. In this work, we show that similar complexes can also be found for organometallic systems containing transition metals and deepen in the reasons behind the existence of these species. We have studied the structures, binding energies, and dissociation profiles in the gas phase of a series of charged hydrogen-bonded dimers of metallocene (Ru, Co, Rh, and Mn) derivatives isoelectronic with the ferrocene dimer. Our results indicate that the carboxylic acid-containing dimers are more strongly bonded and present larger barriers to dissociation than the amide ones and that the cationic complexes tend to be more stable than the anionic ones. Additionally, we describe for the first time the symmetric proton transfer that can occur while in the metastable phase. Finally, we use a density-based energy decomposition analysis to shine light on the nature of the interaction between the dimers.

Production of Dihydrogen Using Ammonia Borane as Reagent and Pyrazole as Catalyst

Theoretical chemistry (DLPNO-CCSD(T)/def2-TZVP//M06-2x/aug-cc-pVDZ) was used to design a system based on ammonia boranes catalyzed by pyrazoles with the aim of producing dihydrogen, nowadays of high interest as clean fuel. The reactivity of ammonia borane and cyclotriborazane were investigated, including catalytic activation through 1H-pyrazole, 4-methoxy-1H-pyrazole, and 4-nitro-1H-pyrazole. The results point toward a catalytic cycle by which, at the same time, ammonia borane can initially store and then, through catalysis, produce dihydrogen and amino borane. Subsequently, amino borane can trimerize to form cyclotriborazane that, in presence of the same catalyst, can also produce dihydrogen. This study proposes therefore a consistent progress in using environmentally sustainable (metal free) catalysts to efficiently extract dihydrogen from small B-N bonded molecules.

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.

Halogen Bonding Involving Isomeric Isocyanide/Nitrile Groups

2,3,5,6-Tetramethyl-1,4-diisocyanobenzene (1), 1,4-diisocyanobenzene (2), and 1,4-dicyanobenzene (3) were co-crystallized with 1,3,5-triiodotrifluorobenzene (1,3,5-FIB) to give three cocrystals, 1·1,3,5-FIB, 2·2(1,3,5-FIB), and 3·2(1,3,5-FIB), which were studied by X-ray diffraction. A common feature of the three structures is the presence of I···Cisocyanide or I···Nnitrile halogen bonds (HaBs), which occurs between an iodine σ-hole and the isocyanide C-(or the nitrile N-) atom. The diisocyanide and dinitrile cocrystals 2·2(1,3,5-FIB) and 3·2(1,3,5-FIB) are isostructural, thus providing a basis for accurate comparison of the two types of noncovalent linkages of C≡N/N≡C groups in the composition of structurally similar entities and in one crystal environment. The bonding situation was studied by a set of theoretical methods. Diisocyanides are more nucleophilic than the dinitrile and they exhibit stronger binding to 1,3,5-FIB. In all structures, the HaBs are mostly determined by the electrostatic interactions, but the dispersion and induction components also provide a noticeable contribution and make the HaBs attractive. Charge transfer has a small contribution (<5%) to the HaB and it is higher for the diisocyanide than for the dinitrile systems. At the same time, diisocyanide and dinitrile structures exhibit typical electron-donor and π-acceptor properties in relation to the HaB donor.

DFT potentials from a chemical perspective: Anatomy of electron (de)localization in molecules and crystals

We introduce a fermionic potential, v f $$ {v}_f $$ , as a comprehensive measure of electron (de)localization in atomic-molecular systems. Unlike other common descriptors as ELF, LOL, etc., it characterizes all physical effects responsible for (de)localization of electrons, namely: an exchange hole depth, its tendency to change, a sensitivity of an exchange correlation hidden in a pair density and kinetic potential to local variations in electron density. Wells in the v f $$ {v}_f $$ distribution correspond to the domains of maximum electron localization, while the potential's barriers prevent delocalization of electrons through them. It also estimates bond orders and successfully reveals the impact of chemical modifications or environmental effects on the delocalization of electrons in molecules and crystals. The v f $$ {v}_f $$ components provide a unique opportunity to compare the influence of the mentioned physical effects on electron (de)localization. This merges physical and chemical views of electron delocalization using functions appearing in density functional theory.

Organic-Cation Modulated Assembly Behaviors of a Ureidopyrimidone-Grafting Cluster

Ureidopyrimidone (UPy) is an important building block for constructing functional supramolecular polymers and soft materials based on their characteristic quadruple hydrogen bonds. While the evidence from the single-crystal X-ray diffraction data for the existence of linear hydrogen bonding has still been absent up to now. To obtain the crystals of UPy-containing molecules with high quality, enhanced rigidity and crystallinity are expected. Herein, an inorganic Anderson-Evans type cluster [Mn(OH₎₆Mo₆O₁₈]^3-, which can provide suitable stiffness and charge, is used as a linker to covalently anchor two UPy units. The prepared organic-inorganic polyanion with three negative charges has a linear architecture, which is prone to form an infinite one-dimensional structure based on the supramolecular forces. The results indicate that the combination models of UPy units can be conveniently modulated by organic counter cations with different sizes, and therefore three unreported models are observed under various conditions. The present study gives a unique understanding of the intermolecular interactions in UPy-based supramolecular polymers and also provides a simple tuning method, which benefits the construction of functional materials and the adjustment of their properties.

The Tetrel Bond and Tetrel Halide Perovskite Semiconductors

The ion pairs [Cs+•TtX3−] (Tt = Pb, Sn, Ge; X = I, Br, Cl) are the building blocks of all-inorganic cesium tetrel halide perovskites in 3D, CsTtX3, that are widely regarded as blockbuster materials for optoelectronic applications such as in solar cells. The 3D structures consist of an anionic inorganic tetrel halide framework stabilized by the cesium cations (Cs+). We use computational methods to show that the geometrical connectivity between the inorganic monoanions, [TtX3−]∞, that leads to the formation of the TtX64− octahedra and the 3D inorganic perovskite architecture is the result of the joint effect of polarization and coulombic forces driven by alkali and tetrel bonds. Depending on the nature and temperature phase of these perovskite systems, the Tt···X tetrel bonds are either indistinguishable or somehow distinguishable from Tt–X coordinate bonds. The calculation of the potential on the electrostatic surface of the Tt atom in molecular [Cs+•TtX3−] provides physical insight into why the negative anions [TtX3−] attract each other when in close proximity, leading to the formation of the CsTtX3 tetrel halide perovskites in the solid state. The inter-molecular (and inter-ionic) geometries, binding energies, and charge density-based topological properties of sixteen [Cs+•TtX3−] ion pairs, as well as some selected oligomers [Cs+•PbI3−]n (n = 2, 3, 4), are discussed.

Intermolecular Covalent Interactions: Nature and Directionality

Quantum chemical methods were employed to analyze the nature and the origin of the directionality of pnictogen (PnB), chalcogen (ChB), and halogen bonds (XB) in archetypal F_m Z⋅⋅⋅F^- complexes (Z=Pn, Ch, X), using relativistic density functional theory (DFT) at ZORA-M06/QZ4P. Quantitative Kohn-Sham MO and energy decomposition analyses (EDA) show that all these intermolecular interactions have in common that covalence, that is, HOMO-LUMO interactions, provide a crucial contribution to the bond energy, besides electrostatic attraction. Strikingly, all these bonds are directional (i.e., F-Z⋅⋅⋅F^- is approximately linear) despite, and not because of, the electrostatic interactions which, in fact, favor bending. This constitutes a breakdown of the σ-hole model. It was shown how the σ-hole model fails by neglecting both, the essential physics behind the electrostatic interaction and that behind the directionality of electron-rich intermolecular interactions. Our findings are general and extend to the neutral, weaker ClI⋅⋅⋅NH₃ , HClTe⋅⋅⋅NH₃ , and H₂ ClSb⋅⋅⋅NH₃ complexes.

Chalcogen Bonding Catalysis with Phosphonium Chalcogenide (PCH)

ConspectusThe exploration of new catalysis concepts and strategies to drive chemical reactions is of vital importance for the sustainable development of organic synthesis. Recently, chalcogen bonding catalysis has emerged as a new concept for organic synthesis and has been demonstrated to be an important synthetic tool capable of addressing elusive reactivity and selectivity issues. This Account describes our progress in the research field of chalcogen bonding catalysis, including (1) the discovery of phosphonium chalcogenide (PCH) as highly efficient chalcogen bonding catalyst; (2) the development of "chalcogen-chalcogen bonding catalysis" and "chalcogen···π bonding catalysis" modes; (3) the demonstration that chalcogen bonding catalysis with PCH can activate hydrocarbons to achieve cyclization and coupling reactions of alkenes; (4) the discovery of unusual results that chalcogen bonding catalysis with PCH can solve elusive reactivity and selectivity issues that are inaccessible by classic catalysis approaches; and (5) the elucidation of chalcogen bonding mechanisms.With PCH catalysts, we systematically studied their chalcogen bonding properties, the relationship between structure and catalysis, and their application in facilitating a diverse array of reactions. Enabled by chalcogen-chalcogen bonding catalysis, an efficient assembly reaction of three molecules of β-ketoaldehyde and one indole derivative in a single operation was realized, delivering heterocycles with a newly constructed seven-membered ring. In addition, a Se···O bonding catalysis approach achieved an efficient synthesis of calix[4]pyrroles. We developed a "dual chalcogen bonding catalysis" strategy to solve reactivity and selectivity issues in the Rauhut-Currier-type reactions and related cascade cyclizations, thus shifting conventionally covalent Lewis base catalysis to a cooperative Se···O bonding catalysis approach. This strategy enables the cyanosilylation of ketones to take place in the presence of a ppm-level amount of PCH catalyst loading. Furthermore, we established chalcogen···π bonding catalysis for catalytic transformation of alkenes. In the research field of supramolecular catalysis, the activation of hydrocarbons such as alkenes by weak interactions is a highly interesting unresolved topic. We showed that the Se···π bonding catalysis approach could efficiently activate alkenes to achieve both coupling and cyclization reactions. Chalcogen···π bonding catalysis with PCH catalysts is particularly highlighted by the capability of facilitating strong Lewis-acid inaccessible transformations, such as the controlled cross coupling of triple alkenes. Overall, this Account presents a panoramic view of our research on chalcogen bonding catalysis with PCH catalysts. The works described in this Account provide a significant platform to solve synthetic problems.

Square Planar Pt(II) Ion as Electron Donor in Pnictogen Bonding Interactions

It has been proposed that late transition metals with low coordination numbers (square planar or linear) can act as nucleophiles and participate in σ-hole interactions as electron donors. This is due to the existence, in this type of metal complexes, of a pair of electrons located at high energy d-orbitals (dz2 or dx2-y2), which are adequate for interacting with antibonding σ-orbitals [σ*(X–Y)] where Y is usually an electron withdrawing element and X an element of the p-block. This type of d[M]→σ*(X–Y) interaction has been reported for metals of groups 9–11 in oxidation states +1 and +2 (d8 and d10) as electron donors and σ-holes located in halogen and chalcogen atoms as electron acceptors. To our knowledge, it has not been described for σ-holes located in pnictogen atoms. In this manuscript, evidence for the existence of pnictogen bonding involving the square planar Pt(II) metal as the electron donor and Sb as the electron acceptor is provided by using an X-ray structure retrieved from the Cambridge Structural Database (CSD) and theoretical calculations. In particular, the quantum theory of atoms in molecules (QTAIM), the noncovalent interaction plot (NCIPlot) and molecular electrostatic potential (MEP) methods were used. Moreover, to further confirm the nature of the Sb···Pt(II) contact, a recently developed method was used where the electron density (ED) and electrostatic potential (ESP) distribution were compared along the Sb···Pt(II) bond path.

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.

Interchangeability and Disorder in the Solid-State Structures of "Two Wall" Calix[4]pyrroles Equipped with Iodine and Ethynyl para-Substituents

Herein, the synthesis and X-ray structures of three α,β "two wall" aryl-extended calix[4]pyrroles having either identical (symmetrically substituted) or different (non-symmetrically substituted) meso-aryl substituents (aryl=4-ethynylphenyl and 4-iodophenyl) are reported. The X-ray structure of the co-crystal formed by the two symmetrically substituted calix[4]pyrroles is also described. In the solid state, all studied α,β-calix[4]pyrroles exhibit a 1,3-alternate conformation with two co-crystallized acetonitrile solvent molecules H-bonded to adjacent cis-pyrrole rings. Remarkably, the 1,3-conformer of the non-symmetrically substituted iodophenyl/ethynylphenyl compound is intrinsically chiral. The two enantiomers are present in the average asymmetric unit in a 65 : 35 occupancy ratio displaying a head-to-tail directional disorder. This is due to the functional complementarity and the isosteric and isoelectronic properties of the para-substituents: iodo and ethynyl. That is, the negative belt of iodine is similar to the negative π-system of the C≡C triple bond and the σ-hole in the iodine atom is similar to the positive proton at the C≡C-H group.

Factors contributing to halogen bond strength and stretch or contraction of internal covalent bond

The halogen bond formed by a series of Lewis acids TF₃X (T = C, Si, Ge, Sn, Pb; X = Cl, Br, I) with NH₃ is studied by quantum chemical calculations. The interaction energy is closely mimicked by the depth of the σ-hole on the X atom as well as the full electrostatic energy. There is a first trend by which the hole is deepened if the T atom to which X is attached becomes more electron-withdrawing: C > Si > Ge > Sn > Pb. On the other hand, larger more polarizable T atoms are better able to transmit the electron-withdrawing power of the F substituents. The combination of these two opposing factors leaves PbF₃X forming the strongest XBs, followed by CF₃X, with SiF₃X engaging in the weakest bonds. The charge transfer from the NH₃ lone pair into the σ*(TX) antibonding orbital tends to elongate the covalent TX bond, and this force is largest for the heavier X and T atoms. On the other hand, the contraction of this bond deepens the σ-hole at the X atom, which would enhance both the electrostatic component and the full interaction energy. This bond-shortening effect is greatest for the lighter X atoms. The combination of these two opposing forces leaves the T-X bond contracting for X = Cl and Br, but lengthening for I.

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.

Structure, Optical and Magnetic Properties of Two Isomeric 2-Bromomethylpyridine Cu(II) Complexes [Cu(C6H9NBr)2(NO3)2] with Very Different Binding Motives

Two isomeric 2-bromomethylpyridine Cu(II) complexes [Cu(C₆H₉NBr)₂(NO₃)₂] with 2-bromo-5-methylpyridine (L¹) and 2-bromo-4-methylpyridine (L²) were synthesized as air-stable blue materials in good yields. The crystal structures were different with [Cu(L¹)₂(NO₃)₂] (CuL¹) crystallizing in the monoclinic space group P2₁/c, while the 4-methyl derivative CuL² was solved and refined in triclinic P1¯. The orientation of the Br substituents in the molecular structure (anti (CuL¹) vs. syn (CuL²) conformations) and the geometry around Cu(II) in an overall 4 + 2 distorted coordination was very different with two secondary (axially elongated) Cu-O bonds on each side of the CuN₂O₂ basal plane in CuL¹ or both on one side in CuL². The two Br substituents in CuL² come quite close to the Cu(II) centers and to each other (Br⋯Br ~3.7 Å). Regardless of these differences, the thermal behavior (TG/DTA) of both materials is very similar with decomposition starting at around 160 °C and CuO as the final product. In contrast to this, FT-IR and Raman frequencies are markedly different for the two isomers and the UV-vis absorption spectra in solution show marked differences in the π-π* absorptions at 263 (CuL²) or 270 (CuL¹) nm and in the ligand-to-metal charge transfer bands at around 320 nm which are pronounced for CuL¹ with the higher symmetry at the Cu(II) center, but very weak for CuL². The T-dependent susceptibility measurements also show very similar results (µ_eff = 1.98 µ_B for CuL¹ and 2.00 µ_B for CuL² and very small Curie-Weiss constants of about -1. The EPR spectra of both complexes show axial symmetry, very similar averaged g values of 2.123 and 2.125, respectively, and no hyper-fine splitting.

Exploring the Dynamical Nature of Intermolecular Hydrogen Bonds in Benzamide, Quinoline and Benzoic Acid Derivatives

The hydrogen bonds properties of 2,6-difluorobenzamide, 5-hydroxyquinoline and 4-hydroxybenzoic acid were investigated by Car-Parrinello and path integral molecular dynamics (CPMD and PIMD), respectively. The computations were carried out in vacuo and in the crystalline phase. The studied complexes possess diverse networks of intermolecular hydrogen bonds (N-H…O, O-H…N and O-H…O). The time evolution of hydrogen bridges gave a deeper insight into bonds dynamics, showing that bridged protons are mostly localized on the donor side; however, the proton transfer phenomenon was registered as well. The vibrational features associated with O-H and N-H stretching were analyzed on the basis of the Fourier transform of the atomic velocity autocorrelation function. The spectroscopic effects of hydrogen bond formation were studied. The PIMD revealed quantum effects influencing the hydrogen bridges providing more accurate free energy sampling. It was found that the N…O or O…O interatomic distances decreased (reducing the length of the hydrogen bridge), while the O-H or N-H covalent bond was elongated, which led to the increase in the proton sharing. Furthermore, Quantum Theory of Atoms in Molecules (QTAIM) was used to give insight into electronic structure parameters. Finally, Symmetry-Adapted Perturbation Theory (SAPT) was employed to estimate the energy contributions to the interaction energy of the selected dimers.

Selenoxides as Excellent Chalcogen Bond Donors: Effect of Metal Coordination

The chalcogen bond has been recently defined by the IUPAC as the attractive noncovalent interaction between any element of group 16 acting as an electrophile and any atom (or group of atoms) acting as a nucleophile. Commonly used chalcogen bond donor molecules are divalent selenium and tellurium derivatives that exhibit two σ-holes. In fact, the presence of two σ-hole confers to the chalcogen bonding additional possibilities with respect to the halogen bond, the most abundant σ-hole interaction. In this manuscript, we demonstrate that selenoxides are good candidates to be used as σ-hole donor molecules. Such molecules have not been analyzed before as chalcogen bond donors, as far as our knowledge extends. The σ-hole opposite to the Se=O bond is adequate for establishing strong and directional ChBs, as demonstrated herein using the Cambridge structural database (CSD) and density functional theory (DFT) calculations. Moreover, the effect of the metal coordination of the selenoxide to transition metals on the strength of the ChB interaction has been analyzed theoretically. The existence of the ChBs has been further supported by the quantum theory of atoms in molecules (QTAIM) and the noncovalent interaction plot (NCIPlot).

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.

Halogen Bonding Involving Gold Nucleophiles in Different Oxidation States

A single-crystal X-ray diffraction (XRD) study of diaryliodonium tetrachloroaurates (or, in the recent terminology, tetrachloridoaurates), [(p-XC₆H₄)₂I][AuCl₄] (X = Cl, 1; Br, 2), was performed for 1 (the structure is denoted as 1a to show similarity with the isomorphic structure 2a) and two polymorphs─2a (obtained from MeOH) and 2b (from 1,2-C₂H₄Cl₂). Examination of the XRD data for these three structures revealed 2-center C-X···Au^III (X = Cl and Br) and 3-center bifurcated C-Br···(Cl-Au) halogen bonding (abbreviated as XB) between the p-Cl or p-Br atoms of the diaryliodonium cations and the gold(III) atom of [AuCl₄]^-. The noncovalent nature of Au^III-involving interactions, the nucleophilicity of the gold(III) atoms, and the electrophilic role of p-X atoms of the diaryliodonium cations in the XBs were studied by a set of complementary computational methods. Combined experimental and theoretical studies allowed the recognition of the d-nucleophilicity of the [d⁸Au^III] atom which, regardless of its rather substantial formal 3+ charge, can function as a d-nucleophilic partner of XB. This conclusion was also supported by theoretical calculations performed for the structures' refcodes BINXOM and ICSD 62511; the obtained data verified the nucleophilicity of Au^III toward a K⁺ ions or a σ-(Cl)-hole, respectively. All our results, together with consideration of relevant literature, indicate that gold atoms in the three oxidation states (0, I, and even III) exhibit nucleophilicity in XBs.

Unique Dispersion-Induced Tetrel Bond with Co-operative σ-hole-Induced Pnicogen Bond in the POCl3-Acetone Heterodimer: Experimental Confirmation at Low Temperatures

Both tetrel and pnicogen bonds are known to be induced through σ-/π-holes. This work reports computational and experimental evidence of the carbonyl carbon of acetone hosting a tetrel bond by dispersion rather electrostatic forces, for the first time, while phosphorus of POCl₃ sustains pnicogen bonding via the σ-hole. Heterodimers of POCl₃ with acetone (CH₃COCH₃) have been isolated within inert gas matrixes of Ar and N₂ at 12 K. Characteristic vibrational bands at P═O stretching of POCl₃ and C═O stretching of CH₃COCH₃ have been obtained in support of the computations. The potential energy surface has been traced computationally using ab initio and density functional methods. CH₃COCH₃ harboring such a tetrel bond, in itself, is quite intriguing. The interplay of these interactions has been comprehended by the quantum theory of atoms in molecules, natural bond orbital, energy decomposition, electrostatic potential mapping, and noncovalent interaction analyses.

Implications of Coexistent Halogen and Hydrogen Bonds in Amorphous Solid Dispersions on Drug Solubility, Miscibility, and Mobility

Specific noncovalent drug-polymer interactions were analytically identified using Raman and Fourier transform infrared spectroscopy for amorphous solid dispersions (ASD) formed between either chlorpropamide or tolbutamide and polyvinylpyrrolidone vinyl acetate random copolymer (PVPVA). Spectral changes in the C-Cl stretching vibrations due to changes in the electronic environment of the Cl atom confirmed halogen bond formation in chlorpropamide-PVPVA ASDs, the extent of which was established to be inversely related to the concentration of the drug using 2D correlation spectroscopy analysis. Hydrogen bonding between the secondary amide of each drug and the pyrrolidone carbonyl of the copolymer was also confirmed in all dispersions. Implications of coexistent interactions were investigated for drug-polymer solubility, mixing free energy, and molecular mobility relative to tolbutamide, which only formed hydrogen bonds with PVPVA. Chlorpropamide had a higher solubility, a larger negative mixing free energy, and lower mobility in PVPVA relative to tolbutamide. These thermodynamic and kinetic differences demonstrate the significance of halogen bond formation even when hydrogen bonding is present.

Non-covalent interactions from a Quantum Chemical Topology perspective

About half a century after its little-known beginnings, the quantum topological approach called QTAIM has grown into a widespread, but still not mainstream, methodology of interpretational quantum chemistry. Although often confused in textbooks with yet another population analysis, be it perhaps an elegant but somewhat esoteric one, QTAIM has been enriched with about a dozen other research areas sharing its main mathematical language, such as Interacting Quantum Atoms (IQA) or Electron Localisation Function (ELF), to form an overarching approach called Quantum Chemical Topology (QCT). Instead of reviewing the latter’s role in understanding non-covalent interactions, we propose a number of ideas emerging from the full consequences of the space-filling nature of topological atoms, and discuss how they (will) impact on interatomic interactions, including non-covalent ones. The architecture of a force field called FFLUX, which is based on these ideas, is outlined. A new method called Relative Energy Gradient (REG) is put forward, which is able, by computation, to detect which fragments of a given molecular assembly govern the energetic behaviour of this whole assembly. This method can offer insight into the typical balance of competing atomic energies both in covalent and non-covalent case studies. A brief discussion on so-called bond critical points is given, highlighting concerns about their meaning, mainly in the arena of non-covalent interactions.

Can We Merge the Weak and Strong Tetrel Bonds? Electronic Features of Tetrahedral Molecules Interacted with Halide Anions

Using the orbital-free quantum crystallography approach, we have disclosed the quantitative trends in electronic features for bonds of different strengths formed by tetrel (Tt) atoms in stable molecular complexes consisting of electrically neutral tetrahedral molecules and halide anions. We have revealed the role of the electrostatic and exchange-correlation components of the total one-electron static potential that are determined by the equilibrium atomic structure and by kinetic Pauli potential, which reflects the spin-dependent electron motion features of the weak and strong bonds. The gap between the extreme positions in the electrostatic and total static potentials along the line linking the Tt atom and halide anion is wide for weak bonds and narrow for strong ones. It is in very good agreement with the number of minima in the Pauli potential between the bounded atoms. This gap exponentially correlates with the exchange-correlation potential in various series with a fixed nucleophilic fragment. A criterion for categorizing the noncovalent tetrel bonds (TtB) based on the potential features is suggested.

Noncovalent interactions between benzochalcogenadiazoles and nitrogen bases

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