Chalcogen bond: a sister noncovalent bond to halogen bond

Triel Bonds with Methyl Groups as Electron Donors. A Pentacoordinate Carbon Atom

The triel bond (TrB) formed between Be(CH3)2/Mg(CH3)2 and TrX3 (Tr=B, Al, and Ga; X=H, F, Cl, Br, and I) is investigated via the MP2/aug-cc-pVTZ(PP) quantum chemical protocol. The C atoms of the methyl groups in M(CH3)2 are characterized by a negative electrostatic potential and act as an electron donor in a triel bond with the π-hole above the Tr atom of planar TrX3. The interaction energy spans a wide range between -2 and -69 kcal/mol. Mg(CH3)2 forms a stronger TrB than does Be(CH3)2, which comports with the more negative electrostatic potential on its methyl groups. Some of the complexes involving Mg display a high degree of transfer of the methyl group from Mg to Tr, which is accompanied by an inversion of the bridging methyl and a sizable pyramidalization of the TrX3 unit. The geometries of these complexes have the properties of the long sought pentacoordinate C which has eluded identification and characterization in the past.

Polarization Upends Convention: Halogen Bonding Propensities of Main Group Halides

The propensities for sigma hole bonding by halogen atoms bonded to central atoms below period 2 in the periodic table remain to be systematically examined. Using iodine as our reference halogen atom, a comprehensive analysis of the tendencies for halogen and other forms of significant sigma hole bonding by simple compounds of main group atoms from H to At is accomplished. An examination of the structure and bonding of complexes formed by those iodine-substituted main group compounds and sigma donating bases (ammonia and trimethylamine) is performed to probe the viability of halogen bonding by heavy main group RnM-I compounds in particular, given the historic focus on period 2. We show that propensities for halogen bonding by FnM-I systems for certain columns of the main group vary anomalously as M gets heavier due to a polarization-induced escalation of the electrostatic potential on I. In certain cases, the positive potential at the sigma hole on I is weaker than that at sigma holes on the central M or geminal R atoms. Previously unexplored cases of strong halogen bonding by the fluoroiodides of heavy group 13 atoms are identified, and prospects for other sigma hole type interactions to polarized (main group) central atoms are elucidated.

An orbital-overlap complement to σ-hole electrostatic potentials

A σ-hole is an electron-deficient region of positive electrostatic potential (ESP) opposite from a half-filled p orbital involved in forming a covalent bond. The σ-hole concept helps rationalize directional noncovalent interactions, known as σ-hole bonds, between covalently bonded group V-VII atoms and electron-pair donors. The magnitude and orientation of σ-holes are correlated with the strength and geometry of halogen bonds. However, ESP computed for isolated σ-holes are not always predictive of interaction energies. For example, the σ-holes of isolated CHFBr2 and isolated CH2FI have identical ESP on the molecule surface, but halogen bonds to these molecules generally have different strengths. Here we show that the compact/diffuse nature of the orbitals involved plays an important role. Our orbital overlap distance quantifies the compact/diffuse nature of the "test orbital" that best overlaps with a systems orbitals at each point. The overlap distance captures the response properties of σ-holes: diffuse σ-holes with large overlap distance are typically "softer" and more polarizable. This aids visualization and interpretation. A linear fit to overlap distance and ESP is predictive of the halogen bond strengths of CH3X and CF3X (X = Cl, Br and I). We suggest that the overlap distance will be a useful partner to ESP for characterizing σ-holes.

Interplay Between the N→C Dative Bond, Intramolecular Chalcogen Bond and π Conjugation in the Complexes Formed by Cyclo[18]carbon and C14 Polyyne with 1,2,5-Chalcogenadiazoles

The N→C dative bond (DB), intramolecular chalcogen bond and π conjugation play important roles in determining the structures and properties of some molecular carbon materials and organic/polymeric photovoltaic materials. In this work, the interplay between the N→C DB, intramolecular chalcogen bond and π conjugation in the complexes formed by cyclo[18]carbon and C14 polyyne with 1,2,5-chalcogenadiazoles has been investigated in detail by using reliable quantum chemical calculations. This study has made four main findings. First, only the Te-containing complexes bound by N→C DBs are much more stable than their corresponding van der Waals (vdW) complexes. Second, in addition to through-bond π conjugations, through-space π conjugations also exist in some Se/Te-containing complexes bound by N→C DBs. Third, the cooperativity between intramolecular chalcogen bond, π conjugation between two monomers and N→C DB is not very strong and can be ignored. Fourth, compared to π conjugations, intramolecular Ch⋅⋅⋅C (Ch=O, S, Se, Te) chalcogen bonds play a secondary role in stabilizing the complexes bound by N→C DBs. These findings clearly indicate that the role of "conformational lock", popular in the field of organic optoelectronic materials, may have been greatly overestimated.

Stable Antiaromatic [16]Triphyrin(2.1.1) with Core Modification: Synthesis Using a 16π Electrocyclic Reaction

Antiaromatic porphyrinoids have attracted significant attention owing to their unique electronic properties and potential applications. However, synthesis of antiaromatic contracted porphyrinoids is challenging owing to the inherent instability associated with smaller ring sizes. In this study, we report the synthesis and characterization of the first stable trioxa[16]triphyrin(2.1.1), a novel 16π antiaromatic contracted porphyrinoid. We utilized a core modification approach to stabilize the [16]triphyrin(2.1.1). X-ray crystallographic analysis revealed a nearly planar structure. Electrochemical studies demonstrated reversible oxidation behavior and a small HOMO-LUMO gap, which was consistent with its antiaromatic nature. Chemical oxidation yielded an aromatic [14]triphyrin(2.1.1) dication, highlighting the antiaromaticity-aromaticity switching capability of this system. This synthesis involved the discovery of a key intermediate, dihydrotrioxatriphyrin(2.1.1), which underwent oxidative dehydrogenation to yield the target compound. Theoretical calculations suggested that dihydrotrioxatriphyrin(2.1.1) formed via a rare 16π electrocyclic reaction. The successful synthesis and characterization of this stable trioxa[16]triphyrin(2.1.1) underscores the potential of the core modification strategies for the rational design of novel antiaromatic systems with tunable properties. Moreover, the discovery of the rare 16π electrocyclic reaction advances the understanding of high-order pericyclic processes and may inspire new synthetic strategies for complex macrocyclic compounds.

On the Versatility of the sp-, sp 2-, and sp 3-Hybridized Chalcogen-Bearing Molecules To Engage in Type I Chalcogen···Chalcogen Interactions: A Quantum Mechanical Investigation of Like···Like and Unlike Complexes

The predilection of sp-, sp 2-, and sp 3-hybridized chalcogen-bearing molecules to engage in type I chalcogen···chalcogen interactions was comparatively unveiled in like···like/unlike CY···YC, YCY···YCY, and F2Y···YF2 (where Y = O, S, and Se) complexes, respectively. Upon the optimized monomers, a potential energy surface (PES) scan was conducted to pinpoint the most favorable complexes. The energetic findings unveiled the ability of the investigated systems to engage in the interactions under study with binding energy values ranging from -0.36 to -2.33 kcal/mol. Notably, binding energies were disclosed to align in the posterior sequence; sp 2- (i.e., YCY···YCY) > sp- (i.e., CY···YC) > sp 3- (i.e., F2Y···YF2) hybridized complexes, except the like···like oxygen-bearing complexes. Instead, the highest negative binding energy values were detected for the OCO···OCO followed by those of the F2O···OF2 and CO···OC complexes. Furthermore, the like···like selenium-bearing complexes demonstrated the most considerable binding energies compared to the other investigated complexes. Remarkably, the quantum theory of atoms in molecules and noncovalent interaction index analyses revealed the highly directional and closed-shell nature of the investigated chalcogen···chalcogen interactions. Symmetry adapted-perturbation theory findings outlined the dispersion forces as the commanding forces for all the studied complexes. These observations will provide convincing justifications for the nature of chalcogens within type I chalcogen···chalcogen interactions, leading to increased progress in various domains regarding drug design and materials science.

An innovative deep eutectic solvent: chalcogen bonding as the primary driving force

Chalcogen bonding (ChB) interactions have drawn intensive attention in the last few decades as interesting alternatives to hydrogen bonding. The applications of ChB were mostly centered on the solid state and have rarely been explored in solution. In this work, a novel strategy for forming ChB-based deep eutectic solvents (DESs) was exploited. We set forth the formation, physicochemical properties, and interaction sites in detail. This work not only provides a new idea to design DES systems but also to exploit the potential application of ChB complexes.

Expanding chalcogen bonds in thiophenes to interactions with halogens

Compounds containing the thiophene moiety find several applications in physics and chemistry, such as electrical conduction, which depends on specific conformations to properly exhibiting the desired properties. In turn, chalcogen bonding has found to modulate the conformation of some N-thiophen-2-ylfomamides. Since halogens participate in a kin interaction (halogen bonding) and are abundant in agrochemicals, pharmaceuticals, and materials, we have quantum-chemically explored the interaction between organic halogen and thiophene as a conformational modulator in some model compounds. Although such interaction indeed appears, as demonstrated by atoms in molecules and natural bond orbital analysis, it is inefficient to control the conformational equilibrium. An energy decomposition analysis scheme demonstrated that halomethane and thiophene tend to move away from one another due to a core component (Pauli repulsion and exchange), which is mainly due to a deformation term. Therefore, chalcogen bonds with halogens appear weaker than with other chalcogens.

Chalcogen Noncovalent Interactions between Diazines and Sulfur Oxides in Supramolecular Circular Chains

The noncovalent chalcogen interaction between SO2/SO3 and diazines was studied through a dispersion-corrected DFT Kohn-Sham molecular orbital together with quantitative energy decomposition analyses. For this, supramolecular circular chains of up to 12 molecules were built with the aim of checking the capability of diazine molecules to detect SO2/SO3 compounds within the atmosphere. Trends in the interaction energies with the increasing number of molecules are mainly determined by the Pauli steric repulsion involved in these σ-hole/π-hole interactions. But more importantly, despite the assumed electrostatic nature of the involved interactions, the covalent component also plays a determinant role in its strength in the involved chalcogen bonds. Noticeably, π-hole interactions are supported by the charge transfer from diazines to SO2/SO3 molecules. Interaction energies in these supramolecular complexes are not only determined by the S···N bond lengths but attractive electrostatic and orbital interactions also determine the trends. These results should allow us to establish the fundamental characteristics of chalcogen bonding based on its strength and nature, which is of relevance for the capture of sulfur oxides.

Chalcogen-Bond-Assisted Formation of the N→C Dative Bonds in the Complexes between Chalcogenadiazoles/Chalcogenatriazoles and Fullerene C60

The existence of the N→C dative bonds in the complexes between N-containing molecules and fullerenes have been verified both theoretically and experimentally. However, finding stable N→C dative bonds is still a highly challenging task. In this work, we investigated computationally the N→C dative bonds in the complexes formed by fullerene C60 with 1,2,5-chalcogenadiazoles, 2,1,3-benzochalcogenadiazoles, and 1,2,4,5-chalcogenatriazoles, respectively. It was found that the N→C dative bonds are formed along with the formation of the N-Ch···C (Ch = S, Se, Te) chalcogen bonds. In the gas phase, from S-containing complexes through Se-containing complexes to Te-containing complexes, the intrinsic interaction energies become more and more negative, which indicates that the N-Ch···C chalcogen bonds can facilitate the formation of the N→C dative bonds. The intrinsic interaction energies are compensated by the large deformation energy of fullerene C60. The total interaction energies of Te-containing complexes are negative, while both total interaction energies of the S-containing complexes and Se-containing complexes are positive. This means that the N→C dative bonds in the Te-containing complexes are more easily observed in experiments in comparison with those in the S-containing complexes and Se-containing complexes. This study provides a new theoretical perspective on the experimental observation of the N→C dative bonds in complexes involving fullerenes. Further, the formation of stable N→C dative bonds in the complexes involving fullerenes can significantly change the properties of fullerenes, which will greatly simulate and expand the application range of fullerenes.

Halogen Bond-Assisted Supramolecular Dimerization of Pyridinium-Fused 1,2,4-Selenadiazoles via Four-Center Se2N2 Chalcogen Bonding

The synthesis and structural characterization of α-haloalkyl-substituted pyridinium-fused 1,2,4-selenadiazoles with various counterions is reported herein, demonstrating a strategy for directed supramolecular dimerization in the solid state. The compounds were obtained through a recently discovered 1,3-dipolar cycloaddition reaction between nitriles and bifunctional 2-pyridylselenyl reagents, and their structures were confirmed by the X-ray crystallography. α-Haloalkyl-substituted pyridinium-fused 1,2,4-selenadiazoles exclusively formed supramolecular dimers via four-center Se···N chalcogen bonding, supported by additional halogen bonding involving α-haloalkyl substituents. The introduction of halogens at the α-position of the substituent R in the selenadiazole core proved effective in promoting supramolecular dimerization, which was unaffected by variation of counterions. Additionally, the impact of cocrystallization with a classical halogen bond donor C6F3I3 on the supramolecular assembly was investigated. Non-covalent interactions were studied using density functional theory calculations and topological analysis of the electron density distribution, which indicated that all ChB, XB and HB interactions are purely non-covalent and attractive in nature. This study underscores the potential of halogen and chalcogen bonding in directing the self-assembly of functional supramolecular materials employing 1,2,4-selenadiazoles derived from recently discovered cycloaddition between nitriles and bifunctional 2-pyridylselenyl reagents.

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

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

The Missing Chalcogen Bonding Donor: Strongly Polarized Oxygen of Water

A biscyclen molecular cabin, synthesized by connecting two cyclen macrocycles with four linkages, entraps a Li+···H2O···Li+ trimer with a water molecule clamped by two Li+ ions. This configuration results in strongly polarized water, characterized by a water proton resonance shift of up to 10.00 ppm. The arrangement facilitates unprecedented O-centered chalcogen bonds between the lone pairs of pyridinyl nitrogen atoms and polarized water oxygen, as confirmed by X-ray crystallography, NMR spectroscopy, and theoretical calculations. Further observation of O-centered chalcogen bonding in a H2O·(LiCl)2 cluster suggests its widespread presence in hydrated salt systems.

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.

Self-Assembly of Wheel-Shaped Nanographdiynes and Self-Template Growth of Graphdiyne

Graphdiyne (GDY) multilayers show stacking-style-dependent physical properties; thus, controlling the stacking style of nanostructures is crucial for utilizing their electrical, optical, and transport properties in electro-optical devices. Herein, we report the assemblies of nanographdiynes decorated with substituents with different steric hindrances to adjust the stacking style. We show that the π-stacked aggregates were influenced by peripheral substituents and the substrate. Steric hexaterphenyl-substituted nanoGDY scaffolds led to dimer structures stacked in the AB-3 configuration with a twist angle of 26.01° or the AB-1 configuration with an in-plane shift along one diyne link. With the interval replacement of steric substituents with long C12 alkyl chains, nanoGDYs were stacked in the AB-2 configuration to decrease the steric congestion, eventually leading to one-dimensional (1D) nanofibrous aggregates. Self-assembly in the presence of substrates can result in ABC-stacked nanoGDYs, which endowed us with the possibility of using nanoGDY as the template for GDY growth in a homogeneous reaction. High-resolution transmission electron microscopy (HRTEM), powder X-ray diffraction (PXRD), X-ray photoelectron spectroscopy (XPS), Raman spectroscopy, and near-infrared-ultraviolet-visible (NIR-UV-vis) absorption spectroscopy indicate that the crystalline GDY prepared in this way is a 1.18 eV bandgap semiconductor.

Affinity of Telluronium Chalcogen Bond Donors for Lewis Bases in Solution: A Critical Experimental-Theoretical Joint Study

Telluronium salts [Ar2 MeTe]X were synthesized, and their Lewis acidic properties towards a number of Lewis bases were addressed in solution by physical and theoretical means. Structural X-ray diffraction analysis of 21 different salts revealed the electrophilicity of the Te centers in their interactions with anions. Telluroniums' propensity to form Lewis pairs was investigated with OPPh3 . Diffusion-ordered NMR spectroscopy suggested that telluroniums can bind up to three OPPh3 molecules. Isotherm titration calorimetry showed that the related heats of association in 1,2-dichloroethane depend on the electronic properties of the substituents of the aryl moiety and on the nature of the counterion. The enthalpies of first association of OPPh3 span -0.5 to -5 kcal mol-1 . Study of the affinity of telluroniums for OPPh3 by state-of-the-art DFT and ab-initio methods revealed the dominant Coulombic and dispersion interactions as well as an entropic effect favoring association in solution. Intermolecular orbital interactions between [Ar2 MeTe]+ cations and OPPh3 are deemed insufficient on their own to ensure the cohesion of [Ar2 MeTe ⋅ Bn ]+ complexes in solution (B=Lewis base). Comparison of Grimme's and Tkatchenko's DFT-D4/MBD-vdW thermodynamics of formation of higher [Ar2 MeTe ⋅ Bn ]+ complexes revealed significant molecular size-dependent divergence of the two methodologies, with MBD yielding better agreement with experiment.

Sigma Hole Potentials as Tools: Quantifying and Partitioning Substituent Effects

Empirical substituent constants, such as the Hammett parameters, have found important utility in organic and other areas of chemistry. They are useful both in predicting the impact of substitutions on chemical processes and in rationalizing after-the-fact observations on chemical bonding and reactivity. We assess the impact of substitutions on monoiodinated benzene rings and find that the modifications that substituents induce on the electrostatic potentials at the sigma hole on the terminal I center correlate strongly with established trends of common substituents. As an alternative to the experimental procedures involved in obtaining empirically based substituent constants, the computationally determined constants based on induced electrostatic potentials offer a model for quantifying the influence of mono- and polyatomic, neutral, and ionic substituents on their compounds. A partitioning scheme is proposed that allows us to discretely separate σ and π contributions to generate quantitative measures of substituent effects.

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

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

Influence of Internal Angular Arrangement on Pnicogen Bond Strength

The three Z-X covalent bonds of a ZX3 unit (Z = P, As, Sb, Bi) are normally arranged in a pyramidal structure. Quantum chemical calculations show that pnicogen bonds (ZBs) to the central Z are weakened if ZX3 is flattened, as in the opening of an umbrella. The partial closing of the umbrella has the opposite effect of substantially strengthening these ZBs, even amounting to a 2- or 3-fold magnification in certain cases. The strongest such bonds, wherein Sb and Bi are in a strained configuration within a ZO3CH model system, have interaction energies of 20 kcal/mol with an NH3 base. Most of these systems, whether flattened or more pyramidal, are capable of engaging in three ZBs simultaneously, despite a certain amount of negative cooperativity.

Light-Controlled Switching of Perylene Bisimide Assemblies

Light-driven changes in supramolecular interactions in perylene bisimides (PBIs) with pendant sulfur-containing functional groups at the bay position are demonstrated. In the ground state, a noncovalent S···X interaction between the σ-hole on sulfur and a heteroatom, X (X = O, N, S), of a neighboring molecule is the main driving force for intermolecular interactions, while in the excited state it is the π-π interaction between PBI scaffolds which drives assembly. The presence of heteroatoms in the solvent results in acceleration of the π-stacking process via the formation of a PBI-solvent complex. The excited-state dynamics involved in the assembly process were revealed via time-resolved fluorescence and transient absorption spectroscopies, while steady-state spectroscopy was used to evaluate the structure of the supramolecular assembly.

Quasi-atomic orbital analysis of halogen bonding interactions

A quasi-atomic orbital analysis of the halogen bonded NH3⋯XF complexes (X = F, Cl, Br, and I) is performed to gain insight into the electronic properties associated with these σ-hole interactions. It is shown that significant sharing of electrons between the nitrogen lone pair of the ammonia molecule and the XF molecule occurs, resulting in a weakening of the X-F bond. In addition, the N-X bond shows increasing covalent character as the size of the halogen atom X increases. While the Mulliken outer complex NH3⋯XF appears to be overall the main species, the strength of the covalent interaction of the N-X bond becomes increasingly similar to that of the N-X bond in the [NH3X]+ cation as the size of X increases.

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.

Elucidating mechanochemical reactivity of a ternary halogen-bonded cocrystal system by computational and calorimetric studies

Discovery of a halogen-bonded ternary cocrystal of 1,3,5-trifluoro-2,4,6-triiodobenzene with pyrazine and triphenylphosphine sulfide has revealed a complex landscape of multicomponent phases, all achievable by mechanochemical interconversion. The observed solid-state reaction pathways were explained by periodic density-functional calculations and comprehensive intermolecular interaction analysis, supported by dissolution calorimetry measurements.

Binding, Sensing, And Transporting Anions with Pnictogen Bonds: The Case of Organoantimony Lewis Acids

Motivated by the discovery of main group Lewis acids that could compete or possibly outperform the ubiquitous organoboranes, several groups, including ours, have engaged in the chemistry of Lewis acidic organoantimony compounds as new platforms for anion capture, sensing, and transport. Principal to this approach are the intrinsically elevated Lewis acidic properties of antimony, which greatly favor the addition of halide anions to this group 15 element. The introduction of organic substituents to the antimony center and its oxidation from the + III to the + V state provide for tunable Lewis acidity and a breadth of applications in supramolecular chemistry and catalysis. The performances of these antimony-based Lewis acids in the domain of anion sensing in aqueous media illustrate the favorable attributes of antimony as a central element. At the same time, recent advances in anion binding catalysis and anion transport across phospholipid membranes speak to the numerous opportunities that lie ahead in the chemistry of these unique main group compounds.

Sigma-Hole and Lone-Pair-Hole Site-Based Interactions of Seesaw Tetravalent Chalcogen-Bearing Molecules with Lewis Bases

For the first time, sigma (σ)- and lone-pair (lp)-hole site-based interactions of SF4 and SeF4 molecules in seesaw geometry with NH3 and FH Lewis bases were herein comparatively investigated. The obtained findings from the electrostatic potential analysis outlined the emergence of sundry holes on the molecular entity of the SF4 and SeF4 molecules, dubbed the σ- and lp-holes. The energetic viewpoint announced splendid negative binding energy values for σ-hole site-based interactions succeeded by lp-hole analogues, which were found to be -9.21 and -0.50 kcal/mol, respectively, for SeF4···NH3 complex as a case study. Conspicuously, a proper concurrence between the strength of chalcogen σ-hole site-based interactions and the chalcogen's atomic size was obtained, whereas a reverse pattern was proclaimed for the lp-hole counterparts. Further, a higher preference for the YF4···NH3 complexes with elevated negative binding energy was promulgated over the YF4···FH ones, indicating the eminent role of Lewis basicity. The indications of the quantum theory of atoms in molecules generally asserted the closed-shell nature of all the considered interactions. The observation of symmetry-adapted perturbation theory revealed the substantial contributing role of the electrostatic forces beyond the occurrence of σ-hole site-based interactions. In comparison, the dispersion forces were specified to govern the lp-hole counterparts. Such emerging findings would be a gate for the fruitful forthcoming applications of chalcogen bonding interactions in crystal engineering and biological systems.

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.

Competition between Binding to Various Sites of Substituted Imidazoliums

The imidazolium cation has a number of different sites that can interact with a nucleophile. Adding a halogen atom (X) or a chalcogen (YH) group introduces the possibility of an NX···nuc halogen or NY···nuc chalcogen bond, which competes against the various H-bonds (NH and CH donors) as well as the lone pair···π interaction wherein the nucleophile lies above the plane of the cation. Substituted imidazoliums are paired with the NH₃ base, and the various different complexes are evaluated by density functional theory (DFT) calculations. The strength of XB and YB increases quickly along with the size and polarizability of the X/Y atom, and this sort of bond is the strongest for the heavier Br, I, Se, and Te atoms, followed by the NH···N H-bond, but this order reverses for Cl and S. The various CH···N H-bonds are comparable to one another and to the lone pair···π bond, all with interaction energies of 10-13 kcal/mol, values which show very little dependence upon the substituent placed on the imidazolium.

Tuning of the Electrostatic Potentials on the Surface of the Sulfur Atom in Organic Molecules: Theoretical Design and Experimental Assessment

Noncovalent sulfur interactions are ubiquitous and play important roles in medicinal chemistry and organic optoelectronic materials. Quantum chemical calculations predicted that the electrostatic potentials on the surface of the sulfur atom in organic molecules could be tuned through the through-space effects of suitable substituents. This makes it possible to design different types of noncovalent sulfur interactions. The theoretical design was further confirmed by X-ray crystallographic experiments. The sulfur atom acts as the halogen atom acceptor to form the halogen bond in the cocrystal between 2,5-bis(2-pyridyl)-1,3,4-thiadiazole and 1,4-diiodotetrafluorobenzene, whereas it acts as the chalcogen atom donor to form the chalcogen bond in the cocrystal between 2,5-bis(3-pyridyl)-1,3,4-thiadiazole and 1,3,5-trifluoro-2,4,6-triiodobenzene.

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

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

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

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

Strategies for the Design of PEDOT Analogues Unraveled: the Use of Chalcogen Bonds and σ-Holes

In this theoretical study, we set out to demonstrate the substitution effect of PEDOT analogues on planarity as an intrinsic indicator for electronic performance. We perform a quantum mechanical (DFT) study of PEDOT and analogous model systems and demonstrate the usefulness of the ωB97X-V functional to simulate chalcogen bonds and other noncovalent interactions. We confirm that the chalcogen bond stabilizes the planar conformation and further visualize its presence via the electrostatic potential surface. In comparison to the prevalent B3LYP, we gain 4-fold savings in computational time and simulate model systems of up to a dodecamer. Implications for design of conductive polymers can be drawn from the results, and an example for self-doped polymers is presented where modulation of the strength of the chalcogen bond plays a significant role.

Fluorescent sensing of non-steroidal anti-inflammatory drugs naproxen and ketoprofen by dansylated squaramide-based receptors

Bis-squaramide receptors L1-L4 bearing a dansyl moiety were synthesised and their potential applications as fluorescent probes towards non steroidal anti-inflammatory drugs naproxen and ketoprofen was investigated. A detailed photophysical characterization in CH₃CN/DMSO solution (9 : 1 v/v) was conducted and demonstrated that the two macrocyclic receptors L1 and L2 show good sensitivity towards ketoprofen with an ON-OFF fluorescent response, while the two open chain receptors L3 and L4 behave similarly with the three guests considered. DFT theoretical calculations carried out on L2 and L4 as model receptors allowed to propose a possible coordination mode towards the guests. Finally, ¹H-NMR spectroscopy in DMSO-d₆/0.5% water solution demonstrated that the four receptors interact with the considered guests via H-bonds.

Se⋯π Chalcogen Bonding in 1,2,4-Selenodiazolium Tetraphenylborate Complexes

The series of substituted 1,2,4-selenodiazolium tetraphenylborate complexes were synthesized via cyclization between 2-pyridylselenylchloride, followed by the anion metathesis, and fully characterized. The utilization of tetraphenylborate anion, a strong π-electron donor via its phenyl rings, promoted the formation of assemblies exhibiting selenium–π interactions. The chalcogen bonding (ChB) interactions involving the π-systems of the tetraphenylborate anion were studied using density functional theory (DFT) calculations, where “mutated” anions were used to estimate the strength of the Se···π chalcogen bonds. Moreover, molecular electrostatic potential (MEP) surfaces were used to investigate the electron-rich and poor regions of the ion pairs. The quantum theory of atoms-in-molecules (QTAIM) and the noncovalent interaction (NCI) plot methods based on the topology of the electron density were used and combined to characterize the ChBs. The investigation reported herein disclosed that the formation of symmetrical dimers can be broken by the introduction of a stronger π-acceptor and, consequently, forming stronger Se···π contacts with selenodiazolium cations.

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.

Robust Supramolecular Dimers Derived from Benzylic-Substituted 1,2,4-Selenodiazolium Salts Featuring Selenium⋯π Chalcogen Bonding

The series of benzylic-substituted 1,2,4-selenodiazolium salts were prepared via cyclization reaction between 2-pyridylselenyl chlorides and nitriles and fully characterized. Substitution of the Cl anion by weakly binding anions promoted the formation supramolecular dimers featuring four center Se₂N₂ chalcogen bonding and two antiparallel selenium⋯π interactions. Chalcogen bonding interactions were studied using density functional theory calculations, molecular electrostatic potential (MEP) surfaces, the quantum theory of atoms-in-molecules (QTAIM), and the noncovalent interaction (NCI) plot. The investigations revealed fundamental role of the selenium⋯π contacts that are stronger than the Se⋯N interactions in supramolecular dimers. Importantly, described herein, the benzylic substitution approach can be utilized for reliable supramolecular dimerization of selenodiazolium cations in the solid state, which can be employed in supramolecular engineering.

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.

Effects of Lewis Basicity and Acidity on σ-Hole Interactions in Carbon-Bearing Complexes: A Comparative Ab Initio Study

The effects of Lewis basicity and acidity on σ-hole interactions were investigated using two sets of carbon-containing complexes. In Set I, the effect of Lewis basicity was studied by substituting the X3/X atom(s) of the NC-C6H2-X3 and NCX Lewis bases (LB) with F, Cl, Br, or I. In Set II, the W-C-F3 and F-C-X3 (where X and W = F, Cl, Br, and I) molecules were utilized as Lewis acid (LA) centers. Concerning the Lewis basicity effect, higher negative interaction energies (Eint) were observed for the F-C-F3∙∙∙NC-C6H2-X3 complexes compared with the F-C-F3∙∙∙NCX analogs. Moreover, significant Eint was recorded for Set I complexes, along with decreasing the electron-withdrawing power of the X3/X atom(s). Among Set I complexes, the highest negative Eint was ascribed to the F-C-F3∙∙∙NC-C6H2-I3 complex with a value of -1.23 kcal/mol. For Set II complexes, Eint values of F-C-X3 bearing complexes were noted within the -1.05 to -2.08 kcal/mol scope, while they ranged from -0.82 to -1.20 kcal/mol for the W-C-F3 analogs. However, Vs,max quantities exhibited higher values in the case of W-C-F3 molecules compared with F-C-X3; preferable negative Eint were ascribed to the F-C-X3 bearing complexes. These findings were delineated as a consequence of the promoted contributions of the X3 substituents. Dispersion forces (Edisp) were identified as the dominant forces for these interactions. The obtained results provide a foundation for fields such as crystal engineering and supramolecular chemistry studies that focus on understanding the characteristics of carbon-bearing complexes.

Stimulus responsive luminescence and application of rotor type 1,1'-([2,2'-bithiophene]-3,3'-diyl)bis(ethan-1-one) and 3'-acetyl-[2,2'-bithiophene]-3-carbaldehyde as molecular rotors

Two dithiophene aldehyde/ketone derivatives are designed as luminescence molecular rotors, i.e., 1,1'-([2,2'-bithiophene]-3,3'-diyl)bis(ethan-1-one) (BTBE) and 3'-acetyl-[2,2'-bithiophene]-3-carbaldehyde (BTAC). Their absorption and luminescence properties, as well as the stimulus responsive luminescence characteristics of water spikes are studied in detail. In order to further explore relationship of luminescence and molecular structure, three reference compounds are also synthesized, named 1-(2-methylthiophen-3-yl)ethanone (MTE), 2-methylthiophene-3-carbaldehyde (MTC) and 4H-cyclohepta[1,2-b:7,6-b']dithiophen-4-one (CDTO). BTBE and BTAC have two obvious absorption bands in the short wavelength region with peak wavelengths of about 212 nm and 260 nm, respectively, while there is a weak trailing type absorption band in the range of about 300-400 nm. Their fluorescence spectra show two luminescence bands in the range of 280-350 nm and 400-600 nm, respectively, and the latter is stronger. Compared with the absorption and luminescence spectra of the reference compounds, it is determined that two absorption bands of BTBE and BTAC in shorter wavelength region are derived from the single thiophene ring carbonyl planar unit, while the absorption band in longer region are derived from the integrated structure of dithiophene carbonyl. The fluorescence bands with peaks of about 300 nm and 470 nm originate respectively from the localized F-C vertical excited states (LE), i.e., the excited state from single planar thienyl-carbonyl unit, and integrated excited state (IE), i.e., the excited state from integrated di-thienyl-carbonyl rings linked covalently with less dihedral angle and greater degree of conjugation at excited state. The crystal structure data show that two thiophene rings possess larger dihedral angles in crystal states, 86.9° for BTBE and 60.8° for BTAC, respectively. However, theoretical calculation results prove the conformational stabilization energy changes little, less than 1.5 kcal/mol, as dihedral angle changes from 50° to 100°. Hydrogen bonding is sufficient to overcome the energy required for this conformational change. Therefore, both BTBE and BTAC can produce water stimulation response luminescence behavior. This stimulating response behavior of BTBE and BTAC can be applied to the preparation of water writable film materials.

Probing the Directionality of S···O/N Chalcogen Bond and Its Interplay with Weak C-H···O/N/S Hydrogen Bond Using Molecular Electrostatic Potential

The directionality of the chalcogen bond (Ch-bond) formed by S and its interplay with other weak interactions have important chemical and biological implications. Here, dimers made of CH₃-S-X and O/N containing nucleophiles are studied and found to be stabilized by coexisting S···O/N and C-H···O/N interactions. Based on experimentally accessible electron density and molecular electrostatic potentials (MESPs), we showed that reciprocity between S···O/N and C-H···O/N interactions in the stability of cumulative molecular interaction (ΔE) was dependent on the strength of the σ-hole on S (V_s,max). Direct correlation between ΔE of dimers with V_s,max of S supports the electrostatic nature of the Ch-bond. Such interplay of the Ch-bond is necessary for its directionality in complex nucleophiles (carbonyl groups) with multiple electron-rich centers, which is explained using MESP. A correlation between the MESP minima in the π-region and the strength of the S-π interaction explains the directional selectivity of the Ch-bond.

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

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

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.