An Overview of Strengths and Directionalities of Noncovalent Interactions: σ-Holes and π-Holes

Insight into the Metal-Involving Chalcogen Bond in the PdII/PtII-Based Complexes: Comparison with the Conventional Chalcogen Bond

The metal-involving Ch···M chalcogen bond and the conventional Ch···O chalcogen bond between ChX2 (Ch = Se, Te; X = CCH, CN) acting as a Lewis acid and M(acac)2 (M = Pd, Pt; Hacac = acetylacetone) acting as a Lewis base were studied by density functional theory calculations. It has been observed that the nucleophilicity of the PtII complexes is higher than that of the corresponding PdII complexes. As a result, the PtII complexes tend to exhibit a more negative interaction energy and larger orbital interaction. The strength of the chalcogen bond increases with the increase of the chalcogen atom and the electronegativity of the substituent on the Lewis acid and vice versa. The metal-involving chalcogen bond shows a typical weak closed-shell noncovalent interaction in the (HCC)2Ch···M(acac)2 complexes, while it exhibits a partially covalent nature in the (NC)2Ch···M(acac)2 complexes. The conventional Ch···O chalcogen bond displays the character of a weak noncovalent interaction, and its strength is generally weaker than that of metal-involving Ch···M interactions. It could be argued that the metal-involving chalcogen bond is primarily determined by the correlation term, whereas the conventional chalcogen bond is mainly governed by the electrostatic interaction.

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

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

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.

Disorder and Sorption Preferences in a Highly Stable Fluoride-Containing Rare-Earth fcu-Type Metal-Organic Framework

Rare-earth (RE) metal-organic frameworks (MOFs) synthesized in the presence of fluorine-donating modulators or linkers are an important new subset of functional MOFs. However, the exact nature of the REaXb core of the molecular building block (MBB) of the MOF, where X is a μ2 or 3-bridging group, remains unclear. Investigation of one of the archetypal members of this family with the stable fcu framework topology, Y-fum-fcu-MOF (1), using a combination of experimental techniques, including high-field (20 T) solid-state nuclear magnetic resonance spectroscopy, has determined two sources of framework disorder involving the μ3-X face-capping group of the MBB and the fumarate (fum) linker. The core of the MBB of 1 is shown to contain a mixture of μ3-F- and (OH)- groups with preferential occupation at the crystallographically different face-capping sites that result in different internally lined framework tetrahedral cages. The fum linker is also found to display a disordered arrangement involving bridging- or chelating-bridging bis-bidentate modes over the fum linker positions without influencing the MBB orientation. This linker disorder will, upon activation, result in the creation of Y3+ ions with potentially one or two additional uncoordinated sites possessing differing degrees of Lewis acidity. Crystallographically determined host-guest relationships for simple sorbates demonstrate the favored sorption sites for N2, CO2, and CS2 molecules that reflect the chemical nature of both the framework and the sorbate species with the structural partitioning of the μ3-groups apparent in determining the favored sorption site of CS2. The two types of disorder found within 1 demonstrate the complexity of fluoride-containing RE-MOFs and highlight the possibility to tune this and other frameworks to contain different proportions and segregations of μ3-face-capping groups and degrees of linker disorder for specifically tailored applications.

Computational Insight into the Nature and Strength of the π-Hole Type Chalcogen∙∙∙Chalcogen Interactions in the XO2∙∙∙CH3YCH3 Complexes (X = S, Se, Te; Y = O, S, Se, Te)

In recent years, the non-covalent interactions between chalcogen centers have aroused substantial research interest because of their potential applications in organocatalysis, materials science, drug design, biological systems, crystal engineering, and molecular recognition. However, studies on π-hole-type chalcogen∙∙∙chalcogen interactions are scarcely reported in the literature. Herein, the π-hole-type intermolecular chalcogen∙∙∙chalcogen interactions in the model complexes formed between XO2 (X = S, Se, Te) and CH3YCH3 (Y = O, S, Se, Te) were systematically studied by using quantum chemical computations. The model complexes are stabilized via one primary X∙∙∙Y chalcogen bond (ChB) and the secondary C-H∙∙∙O hydrogen bonds. The binding energies of the studied complexes are in the range of -21.6~-60.4 kJ/mol. The X∙∙∙Y distances are significantly smaller than the sum of the van der Waals radii of the corresponding two atoms. The X∙∙∙Y ChBs in all the studied complexes except for the SO2∙∙∙CH3OCH3 complex are strong in strength and display a partial covalent character revealed by conducting the quantum theory of atoms in molecules (QTAIM), a non-covalent interaction plot (NCIplot), and natural bond orbital (NBO) analyses. The symmetry-adapted perturbation theory (SAPT) analysis discloses that the X∙∙∙Y ChBs are primarily dominated by the electrostatic component.

Halide Complexes of 5,6-Dicyano-2,1,3-Benzoselenadiazole with 1 : 4 Stoichiometry: Cooperativity between Chalcogen and Hydrogen Bonding

The [M4 -Hal]- (M=the title compound; Hal=Cl, Br, and I) complexes were isolated in the form of salts of [Et4 N]+ cation and characterized by XRD, NMR, UV-Vis, DFT, QTAIM, EDD, and EDA. Their stoichiometry is caused by a cooperative interplay of σ-hole-driven chalcogen (ChB) and hydrogen (HB) bondings. In the crystal, [M4 -Hal]- are connected by the π-hole-driven ChB; overall, each [Hal]- is six-coordinated. In the ChB, the electrostatic interaction dominates over orbital and dispersion interactions. In UV-Vis spectra of the M+[Hal]- solutions, ChB-typical and [Hal]- -dependent charge-transfer bands are present; they reflect orbital interactions and allow identification of the individual [Hal]- . However, the structural situation in the solutions is not entirely clear. Particularly, the UV-Vis spectra of the solutions are different from the solid-state spectra of the [Et4 N]+ [M4 -Hal]- ; very tentatively, species in the solutions are assigned [M-Hal]- . It is supposed that the formation of the [M4 -Hal]- proceeds during the crystallization of the [Et4 N]+ [M4 -Hal]- . Overall, M can be considered as a chromogenic receptor and prototype sensor of [Hal]- . The findings are also useful for crystal engineering and supramolecular chemistry.

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

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

About Perfluoropolyhedranes, Their Electron-Accepting Ability and Questionable Supramolecular Hosting Capacity

Polyhedral molecules are appealing for their eye-catching architecture and distinctive chemistry. Perfluorination of such, often greatly strained, compounds is a momentous challenge. It drastically changes the electron distribution, structure and properties. Notably, small high-symmetry perfluoropolyhedranes feature a centrally located, star-shaped low-energy unoccupied molecular orbital that can host an extra electron within the polyhedral frame, thus producing a radical anion, without loss of symmetry. This predicted electron-hosting capacity was definitively established for perfluorocubane, the first perfluorinated Platonic polyhedrane to be isolated pure. Hosting atoms, molecules, or ions in such "cage" structures is, however, all but forthright, if not illusionary, offering no easy access to supramolecular constructs. While adamantane and cubane have fostered numerous applications in materials science, medicine, and biology, specific uses for their perfluorinated counterparts remain to be established. Some aspects of highly fluorinated carbon allotropes, such as fullerenes and graphite, are briefly mentioned for context.

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

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

Substituent Effects in Tetrel Bonds Involving Aromatic Silane Derivatives: An ab initio Study

In this manuscript substituent effects in several silicon tetrel bonding (TtB) complexes were investigated at the RI-MP2/def2-TZVP level of theory. Particularly, we have analysed how the interaction energy is influenced by the electronic nature of the substituent in both donor and acceptor moieties. To achieve that, several tetrafluorophenyl silane derivatives have been substituted at the meta and para positions by several electron donating and electron withdrawing groups (EDG and EWG, respectively), such as -NH₂, -OCH₃, -CH₃, -H, -CF₃ and -CN substituents. As electron donor molecules, we have used a series of hydrogen cyanide derivatives using the same EDGs and EWGs. We have obtained the Hammett's plots for different combinations of donors and acceptors and in all cases we have obtained good regression plots (interaction energies vs. Hammet's σ parameter). In addition, we have used the electrostatic potential (ESP) surface analysis as well as the Bader's theory of atoms in molecules (AIM) and noncovalent interaction plot (NCI plot) techniques to further characterize the TtBs studied herein. Finally, a Cambridge Structural Database (CSD) inspection was carried out, retrieving several structures where halogenated aromatic silanes participate in tetrel bonding interactions, being an additional stabilization force of their supramolecular architectures.

Xenon Derivatives as Aerogen Bond-Donating Catalysts for Organic Transformations: A Theoretical Study on the Metaphorical "Spherical Cow in a Vacuum" Provides Insights into Noncovalent Organocatalysis

Computations indicate that cationic and noncharged xenon derivatives should exhibit higher catalytic activity than their iodine-based noncovalent organocatalytic congeners. Perfluorophenyl xenonium(II) is expected to demonstrate the best balance between catalytic activity and chemical stability for use in organocatalysis. Comparing its catalytic activity with that of isoelectronic perfluoroiodobenzene indicates that the high catalytic activity of cationic noncovalent organocatalysts is predominantly attributed to the electrostatic interactions with the reaction substrates, which cause the polarization of ligated species during the reaction progress. In contrast, the electron transfer and covalent contributions to the bonding between the catalyst and substrate have negligible effects. The dominant effect of electrostatic interactions results in a strong negative correlation between the calculated Gibbs free energies of activation for the modeled reactions and the highest potentials of the σ-holes on the central atoms of the catalysts. No such correlation is observed for noncharged catalysts.

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.

Stacking Interactions: A Supramolecular Approach to Upgrade Weak Halogen Bond Donors

The co-crystallization of tetracyanobenzene (TCB) with haloarenes ArX provided six new co-crystals TCB ⋅ ArX (ArX=PhCl, PhBr, 4-MeC6 H4 Cl, 4-MeC6 H4 Br, 4-MeOC6 H4 Cl, 1,2-Br2 C6 H4 ) which were studied by X-ray diffraction. In these systems, the strong collective effect of π⋅⋅⋅π stacking interactions and lone pair-(X)⋅⋅⋅π-hole-(C) bondings between TCB and ArX promote the strength of X⋅⋅⋅Ncyano halogen bonding (HaB). Theoretical studies showed that the stacking interactions affect the σ-hole depth of the haloarenes, thus significantly boosting their ability to function as HaB donors. According to the molecular electrostatic potential calculations, the σ- hole-(Cl) value (1.5 kcal/mol) in the haloarene 4-MeOC6 H4 Cl (featuring an electron-rich arene moiety and exhibiting very poor σ-hole-(Cl) ability) increases significantly in the stacked trimer (TCB)2  ⋅ 4-MeOC6 H4 Cl (12.5 kcal/mol). Theoretical DFT calculations demonstrate the dramatic increase of X⋅⋅⋅Ncyano HaB strength for stacked trimers in comparison with parent unstacked haloarenes.

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.

Halonium, chalconium, and pnictonium salts as noncovalent organocatalysts: a computational study on relative catalytic activity

This theoretical study sheds light on the relative catalytic activity of pnictonium, chalconium, and halonium salts in reactions involving elimination of chloride and electrophilic activation of a carbonyl group. DFT calculations indicate that for cationic aromatic onium salts, values of the electrostatic potential on heteroatom σ-holes gradually increase from pnictogen- to halogen-containing species. The higher values of the potential on the halogen atoms of halonium salts result in the overall higher catalytic activity of these species, but in the case of pnictonium and chalconium cations, weak interactions from the side groups provide an additional stabilization effect on the reaction transition states. Based upon quantum-chemical calculations, the catalytic activity of phosphonium(V) and arsenonium(V) salts is expected to be too low to obtain effective noncovalent organocatalytic compounds, whereas stibonium(V), telluronium(IV) and iodonium(III) salts exhibit higher potential in application as noncovalent organocatalysts.

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

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

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

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

Metal Centers as Nucleophiles: Oxymoron of Halogen Bond-Involving Crystal Engineering

This review highlights recent studies discovering unconventional halogen bonding (HaB) that involves positively charged metal centers. These centers provide their filled d‐orbitals for HaB, and thus behave as nucleophilic components toward the noncovalent interaction. This role of some electron‐rich transition metal centers can be considered an oxymoron in the sense that the metal is, in most cases, formally cationic; consequently, its electron donor function is unexpected. The importance of Ha⋅⋅⋅d‐[M] (Ha=halogen; M is Group 9 (Rh, Ir), 10 (Ni, Pd, Pt), or 11 (Cu, Au)) interactions in crystal engineering is emphasized by showing remarkable examples (reported and uncovered by our processing of the Cambridge Structural Database), where this Ha⋅⋅⋅d‐[M] directional interaction guides the formation of solid supramolecular assemblies of different dimensionalities.

On the Importance of σ–Hole Interactions in Crystal Structures

Elements from groups 14–18 and periods 3–6 commonly behave as Lewis acids, which are involved in directional noncovalent interactions (NCI) with electron-rich species (lone pair donors), π systems (aromatic rings, triple and double bonds) as well as nonnucleophilic anions (BF4−, PF6−, ClO4−, etc.). Moreover, elements of groups 15 to 17 are also able to act as Lewis bases (from one to three available lone pairs, respectively), thus presenting a dual character. These emerging NCIs where the main group element behaves as Lewis base, belong to the σ–hole family of interactions. Particularly (i) tetrel bonding for elements belonging to group 14, (ii) pnictogen bonding for group 15, (iii) chalcogen bonding for group 16, (iv) halogen bonding for group 17, and (v) noble gas bondings for group 18. In general, σ–hole interactions exhibit different features when moving along the same group (offering larger and more positive σ–holes) or the same row (presenting a different number of available σ–holes and directionality) of the periodic table. This is illustrated in this review by using several examples retrieved from the Cambridge Structural Database (CSD), especially focused on σ–hole interactions, complemented with molecular electrostatic potential surfaces of model systems.

Predicting the catalytic activity of azolium-based halogen bond donors: an experimentally-verified theoretical study

This report demonstrates the successful application of electrostatic surface potential distribution analysis for evaluating the relative catalytic activity of a series of azolium-based halogen bond donors. A strong correlation (R² > 0.97) was observed between the positive electrostatic potential of the σ-hole on the halogen atom and the Gibbs free energy of activation of the model reactions (i.e., halogen abstraction and carbonyl activation). The predictive ability of the applied approach was confirmed experimentally. It was also determined that the catalytic activity of azolium-based halogen bond donors was generally governed by the structure of the azolium cycle, whereas the substituents on the heterocycle had a limited impact on the activity. Ultimately, this study highlighted four of the most promising azolium halogen bond donors, which are expected to exhibit high catalytic activity.

Quantifying magnetic anisotropy using X-ray and neutron diffraction

In this work, the magnetic anisotropy in two iso-structural distorted tetrahedral Co(II) complexes, CoX 2tmtu2 [X = Cl(1) and Br(2), tmtu = tetra-methyl-thio-urea] is investigated, using a combination of polarized neutron diffraction (PND), very low-temperature high-resolution synchrotron X-ray diffraction and CASSCF/NEVPT2 ab initio calculations. Here, it was found consistently among all methods that the compounds have an easy axis of magnetization pointing nearly along the bis-ector of the compression angle, with minute deviations between PND and theory. Importantly, this work represents the first derivation of the atomic susceptibility tensor based on powder PND for a single-molecule magnet and the comparison thereof with ab initio calculations and high-resolution X-ray diffraction. Theoretical ab initio ligand field theory (AILFT) analysis finds the d xy orbital to be stabilized relative to the d xz and d yz orbitals, thus providing the intuitive explanation for the presence of a negative zero-field splitting parameter, D, from coupling and thus mixing of d xy and . Experimental d-orbital populations support this interpretation, showing in addition that the metal-ligand covalency is larger for Br-ligated 2 than for Cl-ligated 1.

Acid-Base and Anion Binding Properties of Tetrafluorinated 1,3-Benzodiazole, 1,2,3-Benzotriazole and 2,1,3-Benzoselenadiazole

The influence of fluorination on the acid-base properties and the capacity of structurally related 6-5 bicyclic compounds - 1,3-benzodiazole 1, 1,2,3-benzotriazole 2 and 2,1,3-benzoselenadiazole 3 to σ-hole interactions, i. e. hydrogen (1 and 2) and chalcogen (3) bondings, is studied experimentally and computationally. The tetrafluorination increases the Brønsted acidity of the diazole and triazole scaffolds and the Lewis acidity of selenadiazole scaffold decreases the basicity. Increased Brønsted acidity facilitates anion binding via the formation of hydrogen bonds; particularly, tetrafluorinated derivative of 1 (compound 4) binds Cl^- . Increased Lewis acidity of tetrafluorinated derivative of 3 (compound 10), however, is not enough for binding with Cl^- and F^- via chalcogen bonds in contrast to previously studied Te analog of 10. It is suggested that the maximum positive values of molecular electrostatic potential at the σ-holes, V_S,max , can be a reasonable metric for design and synthesis of new anion receptors with selenadiazole-diazole/triazole hybrids as a special target. Related chlorinated compounds are also discussed.

The π-hole revisited

It follows from the Schrödinger equation that the forces operating within molecules and molecular complexes are Coulombic, which necessarily entails both electrostatics and polarization. A common and important class of molecular complexes is due to π-holes. These are molecular regions of low electronic density that are perpendicular to planar portions of the molecular frameworks. π-Holes often have positive electrostatic potentials associated with them, which result in mutually polarizing attractive forces with negative sites such as lone pairs, π electrons or anions. In many molecules, π-holes correspond to a flattening of the electronic density surface but in benzene derivatives and in polyazines the π-holes are craters above and below the rings. The interaction energies of π-hole complexes can be expressed quite well in terms of regression relationships that account for both the electrostatics and the polarization. There is a marked gradation in the interaction energies, from quite weak (about -2 kcal mol-1) to relatively strong (about -40 kcal mol-1). Gradations are also evident in the ratios of the intermolecular separations to the sums of the respective van der Waals radii and in the gradual transition of the π-hole atoms from trigonal to quasi-tetrahedral configurations. These trends are consistent with the concept that chemical interactions form a continuum, from very weak to very strong.

Theoretical study on the noncovalent interactions involving triplet diphenylcarbene

The properties of some types of noncovalent interactions formed by triplet diphenylcarbene (DPC³) have been investigated by means of density functional theory (DFT) calculations and quantum theory of atoms in molecule (QTAIM) studies. The DPC³···LA (LA = AlF₃, SiF₄, PF₅, SF₂, ClF) complexes have been analyzed from their equilibrium geometries, binding energies, and properties of electron density. The triel bond in the DPC³···AlF₃ complex exhibits a partially covalent nature, with the binding energy - 65.7 kJ/mol. The tetrel bond, pnicogen bond, chalcogen bond, and halogen bond in the DPC³···LA (LA = SiF₄, PF₅, SF₂, ClF) complexes show the character of a weak closed-shell noncovalent interaction. Polarization plays an important role in the formation of the studied complexes. The strength of intermolecular interaction decreases in the order LA = AlF₃ > ClF > SF₂ > SiF₄ > PF₅. The electron spin density transfers from the radical DPC³ to ClF and SF₂ in the formation of halogen bond and chalcogen bond, but for the DPC³···AlF₃/SiF₄/PF₅ complexes, the transfer of electron spin density is minimal.

The competition and cooperativity of hydrogen/halogen bond and π-hole bond involving the heteronuclear ethylene analogues

The noncovalent interactions involving heteronuclear ethylene analogues H₂ CEH₂ (E = Si, Ge and Sn) have been studied by the Møller-Plesset perturbation theory to investigate the competition and cooperativity between the hydrogen/halogen bond and π-hole bond. H₂ CEH₂ has a dual role of being a Lewis base and acid with the region of π-electron accumulation above the carbon atom and the region of π-electron depletion (π-hole) above the E atom to participate in the NCX···CE (X = H and Cl) hydrogen/halogen bond and CE···NCY (Y = H, Cl, Li and Na) π-hole bond, respectively. When HCN/ClCN interacts with H₂ CEH₂ by two sites, the strength of hydrogen bond/halogen bond is stronger than that of π-hole bond. The π-hole bond becomes obviously stronger when the metal substituent of YCN (Y = Li and Na) interacting with H₂ CEH₂ , showing the character of partial covalent, its strength is much greater than that of hydrogen/halogen bond. In the ternary complexes, both hydrogen/halogen bond and π-hole bond are simultaneously strengthened compared to those in the binary complexes, especially in the systems containing alkali metal.

Hydrogen versus tetrel bonds in complexes of 3-oxetanone with water and formaldehyde

The ability and preference of 3-oxetanone to form hydrogen or tetrel bonds have been investigated in its complexes with water and formaldehyde by using Fourier transform microwave spectroscopy complemented with quantum chemical calculations. Different types of interactions and internal dynamics have been observed in the targeted complexes. With water, the ether oxygen of 3-oxetanone is the favoured interaction site forming a classical O-HO hydrogen bond. Quite differently, the carbonyl group of 3-oxetanone plays the dual role as a tetrel donor and a proton acceptor in the 3-oxetanone-formaldehyde complex, featuring the CO tetrel bond and C-HO weak hydrogen bond interactions. Splittings originated from the internal rotation of formaldehyde around its C₂ axis were also observed. The V₂ barrier was estimated to be 375(10) cm^-1 based on Meyer's one-dimensional flexible model. The changes in geometries and electronic densities upon complexation would shed light on the impact of archetype solvent and organic substrate molecules on the reactivity of 3-oxetanone.

Bifurcated Halogen Bonding Involving Two Rhodium(I) Centers as an Integrated σ-Hole Acceptor

The complexes [RhX(COD)]2 (X = Cl, Br; COD = 1,5-cyclooctadiene) form cocrystals with σ-hole iodine donors. X-ray diffraction studies and extensive theoretical considerations indicate that the d z 2-orbitals of two positively charged rhodium(I) centers provide sufficient nucleophilicity to form a three-center halogen bond (XB) with the σ-hole donors. The two metal centers function as an integrated XB acceptor, providing assembly via a metal-involving XB.

Enantioseparation of 5,5'-Dibromo-2,2'-Dichloro-3-Selanyl-4,4'-Bipyridines on Polysaccharide-Based Chiral Stationary Phases: Exploring Chalcogen Bonds in Liquid-Phase Chromatography

The chalcogen bond (ChB) is a noncovalent interaction based on electrophilic features of regions of electron charge density depletion (σ-holes) located on bound atoms of group VI. The σ-holes of sulfur and heavy chalcogen atoms (Se, Te) (donors) can interact through their positive electrostatic potential (V) with nucleophilic partners such as lone pairs, π-clouds, and anions (acceptors). In the last few years, promising applications of ChBs in catalysis, crystal engineering, molecular biology, and supramolecular chemistry have been reported. Recently, we explored the high-performance liquid chromatography (HPLC) enantioseparation of fluorinated 3-arylthio-4,4′-bipyridines containing sulfur atoms as ChB donors. Following this study, herein we describe the comparative enantioseparation of three 5,5′-dibromo-2,2′-dichloro-3-selanyl-4,4′-bipyridines on polysaccharide-based chiral stationary phases (CSPs) aiming to understand function and potentialities of selenium σ-holes in the enantiodiscrimination process. The impact of the chalcogen substituent on enantioseparation was explored by using sulfur and non-chalcogen derivatives as reference substances for comparison. Our investigation also focused on the function of the perfluorinated aromatic ring as a π-hole donor recognition site. Thermodynamic quantities associated with the enantioseparation were derived from van’t Hoff plots and local electron charge density of specific molecular regions of the interacting partners were inspected in terms of calculated V. On this basis, by correlating theoretical data and experimental results, the participation of ChBs and π-hole bonds in the enantiodiscrimination process was reasonably confirmed.

Supramolecular architectures sustained by delocalised C–I⋯π(arene) interactions in molecular crystals and the propensity of their formation

A survey of delocalised C–I⋯π(chelate ring) interactions is presented.

On the importance of RH3C⋯N tetrel bonding interactions in the solid state of a dinuclear zinc complex with a tetradentate Schiff base ligand

The tetrel bonding and π-stacking interactions in a new dinuclear zinc complex using a tetradentate N₂O₂ donor Schiff base have been analysed energetically using DFT calculations and several computational tools.

Anion-anion and anion-neutral triel bonds

The ability of a TrCl4- anion (Tr = Al, Ga, In, Tl) to engage in a triel bond with both a neutral NH3 and CN- anion is assessed by ab initio quantum calculations in both the gas phase and in aqueous medium. Despite the absence of a positive σ or π-hole on the Lewis acid, strong triel bonds can be formed with either base. The complexation involves an internal restructuring of the tetrahedral TrCl4- monomer into a trigonal bipyramid shape, where the base can occupy either an axial or equatorial position. Although this rearrangement requires a substantial investment of energy, it aids the complexation by imparting a much more positive MEP to the site that is to be occupied by the base. Complexation with the neutral base is exothermic in the gas phase and even more so in water where interaction energies can exceed 30 kcal mol-1. Despite the long-range coulombic repulsion between any pair of anions, CN- can also engage in a strong triel bond with TrCl4-. In the gas phase, complexation is endothermic, but dissociation of the metastable dimer is obstructed by an energy barrier. The situation is entirely different in solution, with large negative interaction energies of as much as -50 kcal mol-1. The complexation remains an exothermic process even after the large monomer deformation energy is factored in.

Spodium bonding in five coordinated Zn(ii): a new player in crystal engineering?

This highlight evidences the existence and importance of spodium bonds (SpB) in solid state structures involving five-coordinated square-pyramidal Zn(ii) spodium atom.

Copper (II)-based halogen-substituted chromone antitumor drug entities: Studying biomolecular interactions with ct-DNA mediated by sigma hole formation and cytotoxicity activity

Copper-based antitumor drug entities 1-3 derived from substituted (F^-, Br^-, -CH₃) 3-formylchromone pharmacophore were synthesized and thoroughly characterized by spectroscopic and single X-ray crystallographic studies. These complexes show structural novelty due to presence of the X-bonds in chromone scaffold which could facilitate higher propensity for nucleic acids via sigma σ-hole interactions. Therefore, structure-activity relationship of 1-3 was studied by performing ct-DNA binding, pBR322 cleavage and cytotoxicity activity to validate their potential to act as chemotherapeutic drug entities. The binding studies of 1-3 with ct- DNA were carried out employing many biophysical techniques and the corroborative results of these experiments showed intercalation mode of binding and the order of binding was found to be 2 > 1 > 3. The structure of drug entities could facilitated strong halogen bonding interaction (in case of 1 &2) and stability of X bond was rationalized by sigma hole region of positive electrostatic potential on the surface of C-X covalent bond, as determined by gas phase B3LYP computational DFT studies. Interestingly, 2 exhibited most avid binding affinity due to presence of Br^- electron withdrawing and polarizable group. Further, cleavage studies of 1-3 with pBR322 plasmid DNA were performed which demonstrated significant cleavage activity, the supercoiled form (Form I) of plasmid DNA was converted to nicked form (Form II) with the appearance of linearized form (Form III) in between two, implicating lethal double strand breaks of DNA. 2 showed predominantly higher cleavage activity following the similar trend as observed for binding studies. The cytotoxicity of the complexes 1-3 was evaluated by MTT assay against the human liver carcinoma (Huh-7) and prostate cancer (DU-145) cell lines; complex 2 exhibited specific and selective cytotoxicity for the DU-145 cancer cell line with LC₅₀ value of 1.6 μM.

Tetrel Bonding Interactions Involving Carbon at Work: Recent Advances in Crystal Engineering and Catalysis

The σ- and π-hole interactions are used to define attractive forces involving elements of groups 12–18 of the periodic table acting as Lewis acids and any electron rich site (Lewis base, anion, and π-system). When the electrophilic atom belongs to group 14, the resulting interaction is termed a tetrel bond. In the first part of this feature paper, tetrel bonds formed in crystalline solids involving sp3-hybridized carbon atom are described and discussed by using selected structures retrieved from the Cambridge Structural Database. The interaction is characterized by a strong directionality (close to linearity) due to the small size of the σ-hole in the C-atom opposite the covalently bonded electron withdrawing group. The second part describes the utilization of two allotropic forms of carbon (C60 and carbon nanotubes) as supramolecular catalysts based on anion–π interactions (π-hole tetrel bonding). This part emphasizes that the π-hole, which is considerably more accessible by nucleophiles than the σ-hole, can be conveniently used in supramolecular catalysis.

Factors Impacting σ- and π-Hole Regions as Revealed by the Electrostatic Potential and Its Source Function Reconstruction: The Case of 4,4'-Bipyridine Derivatives

Positive electrostatic potential (V) values are often associated with σ- and π-holes, regions of lower electron density which can interact with electron-rich sites to form noncovalent interactions. Factors impacting σ- and π-holes may thus be monitored in terms of the shape and values of the resulting V. Further precious insights into such factors are obtained through a rigorous decomposition of the V values in atomic or atomic group contributions, a task here achieved by extending the Bader-Gatti source function (SF) for the electron density to V. In this article, this general methodology is applied to a series of 4,4'-bipyridine derivatives containing atoms from Groups VI (S, Se) and VII (Cl, Br), and the pentafluorophenyl group acting as a π-hole. As these molecules are characterized by a certain degree of conformational freedom due to the possibility of rotation around the two C-Ch bonds, from two to four conformational motifs could be identified for each structure through conformational search. On this basis, the impact of chemical and conformational features on σ- and π-hole regions could be systematically evaluated by computing the V values on electron density isosurfaces (V_S) and by comparing and dissecting in atomic/atomic group contributions the V_S maxima (V_S,max) values calculated for different molecular patterns. The results of this study confirm that both chemical and conformational features may seriously impact σ- and π-hole regions and provide a clear analysis and a rationale of why and how this influence is realized. Hence, the proposed methodology might offer precious clues for designing changes in the σ- and π-hole regions, aimed at affecting their potential involvement in noncovalent interactions in a desired way.

Dominance of unique Pπ phosphorus bonding with π donors: evidence using matrix isolation infrared spectroscopy and computational methodology

Albeit the first account of hypervalentπ interactions has been reported with halogenπ interactions, the feasibility of their extension to other hypervalent atoms as possible Lewis acids is still open. In this work, the role of phosphorus as an acceptor from the π electron cloud (Pπ pnicogen or phosphorus bonding) in PCl3-C2H2 and PCl3-C2H4 heterodimers is explored, by combining matrix isolation infrared spectroscopy with ab initio and DFT computational methodologies. The respective potential energy surfaces of the PCl3-C2H2 and PCl3-C2H4 heterodimers reveal unique minima stabilized by a concert of reasonably strong to weak interactions, of which Pπ phosphorus bonding was energetically dominant. Heterodimers, trimers and tetramers bound primarily by this unique phosphorus bond were generated at low temperatures. The dominance of phosphorus bonding in the PCl3-C2H2 and PCl3-C2H4 heterodimers over other interactions (such as Hπ, HCl, HP, Clπ and lone pair-π interactions) was confirmed and substantiated using extended quantum theory of atoms in molecules, natural bond orbital, electrostatic potential mapping and energy decomposition analyses. The following inferences in correlation with results from non-covalent-interaction analysis offer a complete understanding of the nature of the Pπ phosphorus bonding interactions. The significance of electrostatic forces kinetically favoring the formation of phosphorus bonded heterodimers, in addition to thermodynamic stabilization, is demonstrated experimentally.

Theoretical Study of Intermolecular Interactions between Critical Residues of Membrane Protein MraYAA and Promising Antibiotic Muraymycin D2

Phospho-N-acetylmuramoyl-pentapeptide translocase (MraYAA) from Aquifex aeolicus is the binding target for the nucleotide antibiotic muraymycin D2 (MD2). MraYAA in the presence of the MD2 ligand has been crystallized and released, while the interactions between the ligand and active-site residues remain less quantitatively and qualitatively defined. We characterized theoretically the key residues involved in noncovalent interactions with MD2 in the MraYAA active site. We applied the quantum theory of atoms in molecules and natural bond orbital analyses based on the density functional theory method on the solved crystal structure of MraY with the MD2 to quantitatively estimate the intermolecular interactions. The obtained results revealed the presence of multiple hydrogen bonds in the investigated active site with strength ranging from van der Waals to covalent limits. Lys70, Asp193, Gly194, Asp196, Gly264, Ala321, Gln305, and His325 are key active-site residues interacting with MD2. Conventional and unconventional hydrogen bonds in addition with charge-dipole and dipole-dipole interactions contribute significantly to stabilize the MD2 binding to the MraYAA active site. It was also found that water molecules inside the active site have substantial effects on its structure stability through hydrogen-bonding interactions with MD2 and the interacting residues.

Variations of quantum electronic pressure under the external compression in crystals with halogen bonds assembled in Cl3-, Br3-, I3-synthons

The inner-crystal quantum electronic pressure was estimated for unstrained C6Cl6, C6Br6, and C6I6 crystals and for those under external compression simulated from 1 to 20 GPa. The changes in its distribution were analyzed for the main structural elements in considered crystals: for triangles of the typical halogen bonds assembled in Hal3-synthons, where Hal = Cl, Br, I; for Hal...Hal stacking interactions, as well as for covalent bonds. Under simulated external compression, the quantum electronic pressure in the intermolecular space reduces as the electron density increases, indicating spatial areas of relatively less crystal resistance to external compression. The most compliant C6Cl6 crystal shows the largest changes of quantum electronic pressure in the centre of Cl3-synthon while the deformation of rigid I3-synthon under external compression depends only on the features of I...I halogen bonds.

Radicalradical chalcogen bonds: CSD analysis and DFT calculations

This manuscript reports a combination of crystallographic analysis (Cambridge Structural Database) and theoretical DFT calculations in chalcogen bonding interactions involving radicals in both the Ch bond (ChB) donor and acceptor. As a radical ChB acceptor (nucleophile) we have used benzodithiazolyl radical (BDTA) and as Ch bond donors (electrophile) we have used dithiadiazolyl and diselenadiazolyl radicals of the general formula p-X-C6F4-CNChChN (Ch = S, and Se). We have evaluated how the para substituent (X) affects the interaction energy, spin density and charge/spin transfer from the electron rich BDTA radical to the electron poor dichalcogenadiazolyl ring. The ability of the latter rings to form ChBs in the solid state has been examined by a comprehensive search in the CSD; several cases are used to exemplify the preferred geometric features of the complexes and they are compared with the theory. The molecular surface electrostatic potentials calculated for these ChB donors allow for a very precise rationalization of the self-assembly motifs most frequently adopted in the crystalline state and of their relative robustness.