Next Article in Journal
Evaluation of Hepatoprotective Activity and Oxidative Stress Reduction of Rosmarinus officinalis L. Shoots Tincture in Rats with Experimentally Induced Hepatotoxicity
Next Article in Special Issue
Computational NMR of Carbohydrates: Theoretical Background, Applications, and Perspectives
Previous Article in Journal
Unprecedented Coordination-Induced Bright Red Emission from Group 12 Metal-Bound Triarylazoimidazoles
Previous Article in Special Issue
Selective Styrene Oxidation to Benzaldehyde over Recently Developed Heterogeneous Catalysts
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Noncovalent Bonds through Sigma and Pi-Hole Located on the Same Molecule. Guiding Principles and Comparisons

by
Wiktor Zierkiewicz
1,*,
Mariusz Michalczyk
1,* and
Steve Scheiner
2
1
Faculty of Chemistry, Wrocław University of Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
2
Department of Chemistry and Biochemistry, Utah State University Logan, Logan, UT 84322-0300, USA
*
Authors to whom correspondence should be addressed.
Molecules 2021, 26(6), 1740; https://doi.org/10.3390/molecules26061740
Submission received: 5 March 2021 / Revised: 16 March 2021 / Accepted: 17 March 2021 / Published: 20 March 2021
(This article belongs to the Special Issue Featured Reviews in Applied Chemistry)

Abstract

:
Over the last years, scientific interest in noncovalent interactions based on the presence of electron-depleted regions called σ-holes or π-holes has markedly accelerated. Their high directionality and strength, comparable to hydrogen bonds, has been documented in many fields of modern chemistry. The current review gathers and digests recent results concerning these bonds, with a focus on those systems where both σ and π-holes are present on the same molecule. The underlying principles guiding the bonding in both sorts of interactions are discussed, and the trends that emerge from recent work offer a guide as to how one might design systems that allow multiple noncovalent bonds to occur simultaneously, or that prefer one bond type over another.

Graphical Abstract

1. Introduction

The concept of the σ-hole, introduced to a wide audience in 2005 at a conference in Prague by Tim Clark [1], influenced a way of thinking about noncovalent interactions that prevails to this day. Early experimental findings [2,3,4,5] of unusual halogen·halogen contacts were explained in part by the anisotropic distribution of electronic density around a halogen atom when linked to an electron-withdrawing group. Already in 1992, it had been learned that electronegative atoms from Groups 14–17 have regions of positive molecular electrostatic potential (MEP) on their outer surfaces, along an extension of a covalent bond, which may attract an incoming Lewis base [6]. This observation led to the further computational studies which resulted in formulation of the σ-hole idea [7] which was further generalized in later works [7,8,9,10,11,12,13,14,15,16,17]. These ideas concerning the halogen bond were successfully adapted to atoms from other families of the periodic table, which were later grouped into the category of σ-hole bonds. This general sort of noncovalent bond has been subdivided by the specific family from which the bridging atom is drawn, i.e., chalcogen [18,19,20,21,22,23], pnicogen (pnictogen) [24,25,26,27,28], or tetrel bonds [29,30,31]. The former, along with the halogen bond, has been formally recognized and detailed in recent IUPAC recommendations [32,33].
Quite the opposite from representing exotic or unusual contacts, these bonds make important contributions to numerous fields of chemistry and biology. As examples, understanding the forces behind crystal engineering and supramolecular chemistry benefits from a knowledge of σ-hole interactions due to their directionality, strength, and self-organization properties which promote formation of adducts in the solid state [34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52]. The importance of σ-hole bonding has also been verified in the context of anion recognition processes [53,54,55,56,57,58,59,60,61], materials chemistry [62,63,64,65,66,67,68,69,70,71,72], or biochemistry [73,74,75,76,77,78,79,80,81]. An early work connecting σ-hole bonds with crucial concepts in chemistry occurred when Grabowski recognized that tetrel bonds can be thought of as a preliminary stage of the very important SN2 reaction [82].
As ideas concerning the σ-hole were proliferating, it was recognized that density depletion is not necessarily limited only to the extensions of covalent bonds. Depletions can also occur above planar groups as well, as for example above a carbonyl or phenyl group. Linear systems such as HCN can also suffer from low density off the molecular axis. Because of their location and association with π-electronic systems, these regions of density depletion and positive electrostatic potential have come to be called π-holes. As one specific example, tricoordinated triel atoms typically occur at the center of a planar triangle, with a π-hole located above the central triel atom [13,83,84,85,86,87], and the resulting triel bond [88,89,90,91,92,93,94,95] falls into the category of a π-hole bond. The π-hole situated above the C atom of a carbonyl group offers another common example, whose presence was manifested in early work of Burgi and Dunitz [96]. Protein structures can also involve participation of π-hole bonds [97,98], as is also true of self-assembling systems [99]. As work has progressed, there has been recognition of π-holes as providing a means by which an aerogen bond (involving rare gas atoms) can form [100,101,102,103,104]. Other newly introduced types of σ/π- hole directed interactions are alkali and alkaline earth bond (e.g., beryllium bond, magnesium bond) in which atoms of 1st and 2nd groups contribute [105,106,107,108,109] or regium and spodium bonds which employ transition metals from 11th (regium [110,111,112,113,114,115]) or 12th (spodium [105,116,117,118,119]) groups of the periodic table. The full range of these sorts of bonds, along with their designations, is summarized in Figure 1.
There have been a number of earlier reviews addressing the issue of σ-hole and π-hole bonds [13,49,51,52,81,84,105,120,121]. However, little attention has been devoted to situations where both hole types are present on a single molecule, and the competition between the two for a nucleophile. Indeed, it is also of intense interest to examine the result when both of these bonds are present at the same time. As has already been explored, the tunability of single σ or π-holes enables the construction of interesting assemblies with desired properties [122,123,124]. The possibilities multiply when both types of bonds are present and influence one another.
The driving goal for this work is to review with a critical eye what is known about systems offering both σ and π-holes to an approaching nucleophile. Are there rules that can be used to predict which will be preferred? Is there a distinction depending on whether both holes lie on the same atom or on different atoms of the same molecule? Different combinations are discussed where π-hole bonds of various types are combined with other noncovalent bonds, whether aerogen, halogen, tetrel, pnicogen, or chalcogen.

2. Origin of π-Holes

The σ-hole has been well documented in the literature, along with its explanation as emanating from the pull of electron density along the axis of a covalent bond. The origin of a π-hole is similar in some ways, as it too relies on anisotropic distribution of charge. The fundamental origin and nature of a π-hole can be understood using F2GeO as an example. This molecule is planar with C2v symmetry. The NBO localized orbitals representing the σ and π bonding orbitals of the Ge=O bond exhibited in the top half of Figure 2 both show a bias toward the more electronegative O. The region above the Ge atom and on its right, toward the two F atoms, thus suffers from a depletion of electron density. This shift is reflected by the drop in the total electron density as the point of reference moves up and out of the molecular plane, culminating in the minimum of the black ρ curve in Figure 3 for θ ~ 60°. Such a density depletion above the molecular plane is commonly referred to as a π-hole, which is in turn responsible for a maximum of the red molecular electrostatic potential (MEP) curve in Figure 3, which occurs at roughly θ = 75°. Figure 4a illustrates this π-hole in the MEP as the blue region lying above (and below) the Ge and shifted slightly toward the two F atoms. It is thus natural for a nucleophile to then approach this π-hole from this direction, as exemplified by the complex with NH3 displayed in Figure 4b. It might be noted as well that the π* orbital in Figure 2 is perfectly situated to act as electron acceptor from the lone pair of the NH3 nucleophile, another trademark of π-hole interactions. Note finally the geometrical distortion of the originally planar F2GeO in Figure 4b as the O and F atoms are shifted away from the approaching nucleophile another common characteristic of π-hole bonds.
Of course, the forgoing explanation of the origin of a π-hole can vary from one system to the next. The presence of any lone electron pairs on the central atom can influence the magnitude and actual location of any such π-hole. An example to be discussed shortly is XeF2O where the central Kr atom contains two such lone pairs. The positive region that appears above an aromatic ring of, e.g., C6F6, is not located directly above any one C atom but lies rather above the ring’s center. Such a dislocation can at times make it problematic to connect a π-hole with any particular atom, such as that above a C≡N group which lies roughly midway between the C and N atom.
It should be reiterated that the causality is as follows: depletion of electron density causes a rise in the electrostatic potential [125]. In fact, sometimes a hole is labeled as such even though the potential is not positive, just less negative than its surroundings. For practical reasons, the magnitude of a hole is typically measured on the 0.001 au electron isodensity surface and is quantified by the VS,max parameter developed for this purpose [12,126]. The electron density can be accessed not only by quantum calculations but also experimentally by diffraction methods [127].
It has been shown that the intensity of a σ-hole can be adjusted by changing the polarizability of the central atom and the electron-withdrawing power of its substituents [13,84,127], and the same considerations apply to π-holes as well [128]. However, it must be borne in mind that the strength of a given interaction is not a function solely of electrostatic considerations. Polarization and dispersion forces are important attractive forces as well [13,17,129,130,131].

3. σ-Hole and π-Hole on the Same Atom

There are only a few reports in the literature of systems where both σ and π holes appear on the same atom. One example of such a situation is furnished by AeF2O where Ae refers to an aerogen atom Kr or Xe [132]. As may be seen in Figure 5, the Ae atom of this molecule contains a σ-hole along the O=Ae bond extension, while a π-hole opens up above the molecular plane. As may be seen in the upper section of Table 1, the σ-hole is quite a bit more intense than the π-hole. The weaker nature of the latter may be attributed in part to the presence of lone pairs on the central Ae atom which share space with these holes and would dilute any density depletion. It is interesting to note that whereas the σ-hole is a bit more intense for Xe vs Kr, the reverse order is seen in their π-holes.
These two AeF2O molecules were allowed to react with a series of diazine nucleophiles (with negative MEP minima on the N atoms) [132]. It was found that the intensity of these various holes carries over into the interaction energies (Eint) of the corresponding dimers. Complexation through the σ-hole was more stable than the corresponding π-hole dimers by about 6 kcal/mol. The stability order was parallel to the hole intensities, as the Xe complexes were slightly stronger than Kr atom analogues.
Another study by Bauza and Frontera [104] of related complexes obtained a deeper σ-hole on XeOF2 by a different computational procedure. XeOF2 complexes with ammonia and CH3CN were characterized by similar interaction energies as in Reference [132]. Interaction energies were magnified by the use of an anion as nucleophile as would be expected. Unfortunately, these authors did not consider π-hole complexes for purposes of comparison.
The ability of XeOF2 to interact with a Lewis base through both hole types was considered [137] in connection with the nucleophile NCCH3, as displayed in Figure 6. The σ-hole complex on the left has the a shorter Xe···N distance than the π-hole structure by nearly 0.4 Å, and is preferred by 2–4 kcal/mol. This energetic advantage occurs despite the presence of a secondary C···O tetrel bond in the latter. The authors questioned the level of covalency in the Xe···N bond, but their topological analyses were not conclusive. Unfortunately, this work did not delve into the σ or π-hole depths as background information.
Brock et al. [138] provided experimental crystal structural evidence of the XeOF2···NCCH3 binding wherein three different nucleophiles approached Xe from three different directions simultaneously, two of π-hole character and one σ-hole, demonstrating the viability of multiple bonds to different holes of the same atom. In a different vein, the σ-hole at Xe atom in XeOF2 has also been a vehicle by which to illustrate cooperativity between regium and aerogen bond in the ternary systems of C2H2···MCN···XeOF2, C2H4···MCN···XeOF2, MCN···C2H4···XeOF2 and C2(CN)4···MCN···XeOF2 where M= Cu, Ag or Au [111]. The σ-hole was able to interact with a N lone pair as well as the unsaturated π- system of a C=C bond. It would have been particularly interesting had the authors considered similar questions with regard to the Xe π-hole.
Within a hypervalent bonding situation, halogen atoms are also capable of containing a π-hole. Within the context of the BrOF2+ cation [139] an X-ray structure indicated the central Br atom can be approached by three nucleophiles. Two KrF2 and one AsF6 molecule attack electrophilic regions on the outer surface of the Br atom. While these three Lewis bases appear to interact through σ-holes on the Br atom, it cannot be excluded that there is a π-hole site on this cation which can attract another Lewis base.
Turning to the pnicogen atom, interaction with NH3 [133] caused ZF2C6H5 (Z = P, As, Sb, Bi) to take on two different arrangements. The first contained three σ-holes in the range of 19-53 kcal/mol, the strongest for the most polarizable Bi atom. As indicated in Figure 7, approach to this σ-hole places the NH3 opposite one of the F atoms. Internal rearrangement of the ZF2C6H5 to a more planar structure opens up a π-hole above the Z atom. As delineated in Table 1, this hole has a magnitude between 36 and 61 kcal/mol, which can also attract the NH3 nucleophile, as depicted in the lower half of Figure 7. The greater depth of the π-holes leads to their larger interaction energies by a factor of 2 to 8. However, the internal rearrangement required to open up these π-holes is energetically costly, requiring between 16 and 43 kcal/mol, so that despite the greater depth of the π vs σ-holes, and their superior interaction energies, it is the set of structures in the top half of Figure 7 that are energetically preferred by a margin between 1 and 11 kcal/mol. This deformation energy diminishes with larger Z atoms, so that there is a closer competition between the σ and π-hole complex binding energies.
Another study of pnicogen bonded systems addressed the question of how many Lewis base ligands can be attached to a single ZF3 molecule [134]. The ZF3 monomer (Z = P, As, Sb, Bi) is characterized by three σ-holes, one opposite each Z-F bond, and a much shallower one that lies directly opposite the Z lone pair, amongst the three F atoms. Although ZF3 is pyramidal, the latter was considered a π-hole due to its placement. Table 1 documents the much lesser VS,max of the latter as compared to the three σ-holes on each ZF3 unit. Because of its greater electrostatic attraction, it is the σ-hole that draws in the approaching nucleophile, whether HCN, CN, or NH3. A second nucleophile of any sort occupies a second σ-hole, but such a triad is only possible for the two heavier Sb and Bi atoms for the anionic CN due to the large Coulombic repulsion required to form such a dianion, and even so, the SbF3··(CN)2 triad has a positive binding energy. HCN is too weak a nucleophile to squeeze in a third pnicogen-bonded base, while such a tetrad is possible for NH3. Addition of a fourth NH3 is possible only for Z = Bi. It seems to occupy the π-hole mentioned above, but this weak interaction is reinforced by secondary noncovalent bonding, reflecting the reluctance of these units to engage with their weak π-holes.
Shifting attention to chalcogen bonding, the YF4 (Y = S, Se, Te, Po) monomer adopts a see-saw equilibrium geometry [135] with a pair of σ-holes lying opposite each of the two equatorial Y-F bonds, with VS,max magnitudes between 42 and 76 kcal/mol, as listed in Table 1. A pair of NH3 nucleophiles can approach along these two σ-holes which would lead to an overall octahedral arrangement with the two NH3 units cis to one another, as depicted on the left side of Figure 8. An alternative trans positioning of the two nucleophiles as shown on the right side of Figure 8 would place the central YF4 in a square pyramidal configuration, closer to a square. As such, it would acquire a π-hole directly below the Y atom. It may be observed from Table 1 that these π-holes are rather intense, between 53 and 64 kcal/mol, so can easily attract one of the two NH3 bases. The other side of the nearly square YF4 unit contains a Y lone electron pair, which obstructs the second NH3 from approaching directly opposite the first, so it situates itself closer to a point opposite one of the Y-F bonds, close to its σ-hole. The energetic comparison of the cis and trans geometries must again consider both interaction and deformation energies. Whereas the interaction energies of the trans structure are much more negative than those for cis, the deformation energies required to adopt the nearly square structure are much larger as well. The net result is that the binding energies of the cis and trans conformations are comparable to one another. Cis is favored for the two smaller S and Se chalcogen atoms, while the heavier Te and Po, with their somewhat reduced deformation energy requirements, shift the equilibrium toward trans.
Similar considerations apply to tetrel bonds. TF4 (T = Si, Ge, Sn, Pb) is of course tetrahedral, so is characterized by four equivalent σ-holes opposite each T-F bond [136]. In order to accommodate a pair of bases, TF4 distorts into an octahedron. Cis approach of the two bases places the TF4 in a see-saw geometry with a pair of σ-holes, each lying opposite an equatorial T-F. An alternate deformation into a square planar structure, allowing for trans arrangement of the two bases, imparts a pair of π-holes to the TF4 unit. The lowermost section of Table 1 shows the σ-holes of the former structure are slightly deeper than the π-holes of the latter. It is worth stressing that either type of hole in these distorted arrangements, whether see-saw or square, is considerably deeper than those in the original undistorted tetrahedral geometry. The interaction energies of the pair of NCH bases with the square was considerably more negative than for the see-saw, despite the larger VS,max for the latter, by between 10 and 40 kcal/mol. On the other hand, the energy required to deform the tetrahedral TF4 into a square far exceeded that needed for the see-saw, particularly for the smaller T atoms. In quantitative terms, the deformation energies for the squares varied from 22 all the way up to 66 kcal/mol for SiF4, as compared to only 0.4–15 kcal/mol for the see-saw. The net result is that it is the cis geometry that is favored, and this preference varies from only 3 kcal/mol for PbF4 up to 23 kcal/mol for SiF4. In fact, the large deformation energies for the square structures with small T make the binding energy a positive quantity for both SiF4 and GeF4.
The forgoing highlights the importance of considering π-holes, even for Lewis acids whose undistorted monomer geometries are such that no such holes are present. The deformations which the monomer undergoes in its interaction with one or more bases can present the possibility of a π-hole whose interaction is comparable to, or even stronger than, that with the original σ-holes.

4. σ-Hole and π-Hole on Different Atoms

More common than the presence of a σ and π-hole arising on the same atom is the situation wherein these two holes are located on different atoms. Within this subgroup there are several themes that are more common than others.

4.1. Combination of σ-Hole on Halogen or Chalcogen with π-hole on Aromatic Ring

A number of works have considered situations wherein a σ-hole on a halogen (X) atom coexists with a π-hole lying above the plane of an aromatic ring. When not directly above the center of the ring, the π-hole is shifted so as to lie closer to the midpoint of a C-C or C-N bond, a topic which has been discussed at some length elsewhere [140]. With this understanding, this review subsumes all of these types into the category of π-hole. Values of the MEP positive maxima are enumerated in Table 2.
One prominent example modeled functionalization of graphene sheets by C6H5Br [141] and its physisorption on the graphene surface through either its Br σ-hole or phenyl π-hole. The adsorption energy with electron-rich graphene regions was three times larger for the π-hole interaction than for the σ-hole contact. Taking the analysis one step further, the authors noted that the adsorption seriously affects the properties of graphene. Bromopentafluorobenzene also has the option of interaction through either its σ or π-hole [142], in this case with pyridine. The latter nucleophile was able to interact either through its N lone pair or the π-electron system of its aromatic ring. As in the previous case the π-hole···π interaction yielded the more stable dimer, here by about 4 kcal/mol. As a fundamental point of distinction, energy decomposition suggested that whereas electrostatics is the dominant factor in σ-hole complexation, this role is assumed by dispersion in π-hole bonding. Another study of this type [143] involved homo-dimerization of 1,3,5-trifluoro-2,4,6-triiodobenzene (TITFB) through its σ or π holes. The σ-hole on the iodine atom was roughly three times deeper than the π-hole positioned above the benzene ring. However, the intermolecular distance was shorter of about 0.1 Å in the case of π-hole interaction.
Table 2. Comparison between σ and π holes appearing on the same molecule, where σ-hole is coupled with halogen (or chalcogen) atom and π-hole with an aromatic ring, VS,max in kcal/mol.
Table 2. Comparison between σ and π holes appearing on the same molecule, where σ-hole is coupled with halogen (or chalcogen) atom and π-hole with an aromatic ring, VS,max in kcal/mol.
MoleculeHole SourceHole TypeVS,maxReferences
1,4-DITFBIσ32.3[140]
Carbon/Benzene ringπ15.1[140]
C6H5BrBrσup to ~15[142]
Benzene Ringπup to ~15[142]
TITFBIσ30.1[143]
Carbon/Benzene Ringπ11.4[143]
Halogen Substituted 1,3,4-oxadiazol-2IJ3H)-thionesS (Chalcogen), Cl, Brσup to ~37[144]
Oxadiazole and Benzyl Ringπup to ~37[144]
HaloperfluorobenzeneCl, Br, Iσ20.9 to 32.8[145]
Benzene Ringπ12.6 to 19.8[145]
XC3H4N2+; X = F, Cl, Br, ICl, Br, Iσ105.9 to 117.5[146]
Imidazolium Ringπ111.0 to 127.9[146]
Shukla et al. [144] found that biologically active derivatives of halogen substituted 1,3,4-oxadiazol-2IJ3H)-thiones in the solid state contain a σ-hole on the S, Cl, and Br atoms and a π-hole above the oxadiazole and benzyl rings all of which participate in noncovalent interactions. The S···N σ-hole chalcogen bond manifested distinctive activity in stabilization of these amalgams. This work fits into the idea that both σ and π-hole bonding play a pivotal role in the 3D organization of crystalline structures and in various molecular scaffolds.
Li et al. provided another example [145] in the context of interactions between haloperfluorobenzene with fluoroanthene (FA) which assemble into nine luminescent cocrystals. These interactions involved the participation of the σ-hole of Cl, Br, or I atoms or the π-hole of the phenyl ring with aromatic π-electrons on the FA. The σ-holes were more intense than the π-hole by as much as 17 kcal/mol. While the balance between σ and π-hole bonding with FA varied as the σ-hole deepened and the π-hole weakened along with the growth of the halogen atom, it was nonetheless the σ-hole···π interaction that was universally favored.
Wang et al. [146] examined the interaction between protonated 2-halogenated imidazolium cation (XC3H4N2+; X = F, Cl, Br, I) and the set of (CH3)3SiY (Y = F, Cl, Br, I) Lewis bases. The acid contained both a σ-hole on its X atom and a π-hole above the imidazolium ring. The π-holes were the deeper of the two, leading to stronger interactions. In these ionic cases, both interaction types were driven mainly by electrostatic forces (the electrostatic term was well correlated with the VS,max values found on monomers) with some addition of polarization and dispersion. The dispersion term was somewhat more prominent for the weaker σ-hole complexes while the polarization component towered over dispersion for π-hole dimers.
An experimental component was contributed by Zhang et al. [147], which paired C6F5X (X = Cl, Br, I) with C6D6 in solution. As in the earlier cases, a σ-hole appeared on the X atom while a π-hole occupied the space above the aromatic ring. For the strongest σ-hole which was associated with the I atom, the interaction lay through its σ-hole, while it was the π-hole link that dominated for the smaller halogens. This pattern was attributed by the authors in large part to entropic contribution since enthalpy alone would not be sufficient to stabilize σ-π over π-π. Comparable conclusions were drawn in another paper by the same group [148], where Lewis acids participated in complexes with deuterated solvent molecules with lone-pair electrons, including CD3CN, CD3COCD3, CD3OD, and [D6]DMSO. Experiment and calculations showed that again σ-hole and π-hole bonds compete with each other for the nucleophile’s lone pair, rather than the π-electron system of C6D6 in the earlier work. The results indicated that only the iodine and a few bromine σ-hole halogen bonds are strong enough to contend successfully with the π-hole interactions in solution.
The biologically relevant ebselen derivative, 2-(2-bromophenyl)benzo[d][1,2]selenazol-3(2H)-one homodimer [149] exhibited markers of simultaneous σ-hole and π-hole bonding. The σ-hole lay along the C-Br covalent bond elongation and the π-hole was derived from the phenyl ring. Shukla et al. postulated two simultaneous interactions: σ-hole bonding from the Br σ-hole to the negative charge-concentrated C-C bond, as well as a π-hole bond between a positively charged C atom and the lone electron pair of Br, confirmed by their NBO analysis. This case represents an uncommon model wherein the same atom acts as both σ-hole bond acceptor and π-hole bond donor.
Finally, Wang et al. [140] reviewed a comparison between complexes stabilized by σ/π holes with various nucleophilic acceptors with particular emphasis on their potential application to anion recognition and transport. Among numerous examples, one of especial interest concerned 1,4-diiodoperfluorobenzene (1,4-DITFB) with both a σ and π-hole. The former was twice as deep as the latter, and the authors underscored that this molecule can act as σ/π hole donor and also by taking into account the high amount of electronic density collected on F atoms, and on the belt around the central fragments of the I atoms surface, it can serve as σ-hole or π-hole acceptor as well.

4.2. Combination of σ-Hole Halogen Bond with Other π-Hole Sources

Further exploration of the literature supplies examples where a X σ-hole is combined with π-holes derived from groups other than aromatic rings. Instances of the halogen-pnicogen combination are most numerous. Lang et al. [150] provided a textbook example in that the surface of XONO/XONO2 (X = F, Cl, Br, I) contains both a typical σ-hole on X and a π-hole localized on the N atom. The VS,max values for this set of systems are listed in the first four rows of Table 3. The authors observed that the σ-hole intensity grew in line with the increase of X atom size whereas the π-hole magnitude dropped in this same order (similar picture as at Li et al. [145]). Binary complexes of these monomers with ClO and ammonia had interaction energies correlated to the MEP maxima. The ternary ClO···XONO/XONO2···NH3 complexes revealed anti-cooperative effect between the σ-hole halogen and π-hole pnicogen bond, as expected when the central unit acts as double electron acceptor. Solimannejad et al. [151] examined adducts of NO2X (X = Cl, Br) with HCN and HNC moieties. The π-hole was stronger than the σ-hole, contrary to the Lang et al. results. Within four tested interaction schemes in binary complexes (stabilized by σ-hole, π-hole and two different hydrogen bond approaches), those bonded through the π-hole were the most stable, for the majority of the trimers displayed negative cooperativity between the σ and π-hole interactions. Subsequent work of this research group [152] deployed the NO2I monomer with a different Lewis base (ammonia). The I atom showed a deeper σ than π-hole, which led to the expected outcome that the NO2I dimer with ammonia was more strongly bound by its σ than by its π-hole. This study also supported two conclusions formulated in earlier cited works: (i) both interaction types showed a strong correlation between VS,max and interaction energy, and (ii) antagonistic effect between σ and π-hole bonds was observed in trimers with two ammonia units.
There are also examples in the literature combining a σ-halogen with a π-tetrel bond. As one example, McDowell [153], analyzed complexes of NCX (X = F, Cl, Br) with H2O. Besides the obvious appearance of a σ-hole on X, a π-hole region was noticed above and below C. As in the previous cases, the intensity of the π-hole drops along with larger X as its σ-hole deepens. Of the two potential binding modes illustrated in Figure 9, only NCF was able to form complexes with water via both o modes a and b. Whereas the π-hole on NCF was 2.5 times stronger than its σ-hole, the π-holes on NCCl and NCBr were too weak to attract the incoming nucleophile. With regard to addition of a third entity, addition of hydrogen and beryllium cations destabilized the σ-hole bond while improving binding in the π-hole complexes.
Bauza and Frontera [154] tested the ability of BH2X (X = F, Cl, Br, I) to establish various sorts of interactions in their homodimers, depending on the attack angle between them. Along with the presence of a σ-hole on X, a π-hole was also available over the central B, perpendicular to the molecular plane. The strengths of these two holes followed expected trends based on X atomic size, and the B π-hole was considerably more intense than the X σ-hole. It was therefore no surprise to find a strong interaction energy of −34.7 kcal/mol for π-hole(B)···σ-hole(I) interaction, much larger than that when I is replaced by Br, or even the halogen bond between I σ-holes in a pair of BH2I molecules.
A particularly unusual combination of holes was examined by Alikhani [109]. The π-hole on BeCl2 occurs as a narrow belt around the Be atom, with VS,max = 32.2 kcal/mol. The magnitude of the σ-hole on Cl is only 1.2 kcal/mol, too weak to expect effective attraction of a Lewis base. The author consequently focused on the nature and properties of the π-hole beryllium bond in complexes with nucleophiles Cl2, NH3, and DMF (dimethylformamide), all of which proved to be quite strong, with the Eint reaching up to −46.3 kcal/mol. One might expect that the replacement of the Cl on BeCl2 by Br or I would intensify the X σ-hole and perhaps also weaken the π-hole on Be. If that were the case, then a XB through the Br or I σ-hole might be able to successfully compete with the Be π-hole interaction.

4.3. Other Examples

The literature contains a few other sorts of combinations of σ with π-holes. Pal et al. [155] present the interplay between these binding sites in Fmoc-Leu-Ψ[CH2NCS] (Fmoc = fluorenylmethyloxycarbonyl protecting group, Leu = leucine) organic isothiocyanate which act as an intermediate in the synthesis of bioactive peptides. The N=C=S group contains a σ-hole on the S atom surface (along the extension of the C=S bond), and a π-hole perpendicular to the C-N bond. The specific location of the π-hole makes it difficult to classify as either a π-hole tetrel or pnicogen bond. For the purpose of the current review, it was classified as π-hole tetrel/pnicogen bond (see Table 4). Between these positively charged sections, there were also regions of negative potential derived from the lone electron pair of the nitrogen (σ-hole acceptor) and sulphur (π-hole acceptor) atoms. N=C=S···N=C=S contacts were observed in the crystal lattice that are associated with electrostatically driven attractions between regions of opposite charge. These simultaneous σ and π-holes interactions with the negative sites were cooperative, thereby magnifying the stabilization within the crystal.
The specific location of the π-hole was more definitive [156] in CF2=CFZH2 (Z = P, As, Sb) and CF2=CFPF2. The electron-deficient regions corresponding to π-holes were positioned directly above the C atom, while the σ-hole was localized on the pnicogen Z atom. The σ-hole became deeper in the sequence P < As < Sb attributed to rising Z polarizability, while the π-hole magnitude lessens in the same order, as documented in Table 4. Switching -PH2 moiety to -PF2 boosted the magnitude of both holes. The prediction of MEP extrema was verified by the energetic properties of the complexes formed between these Lewis acids and NH3 or NMe3. The interaction energies of the dyads were generally consistent with the MEP trends. However, some irregularity was found, as the σ-hole complexes in each case were more stable than their π-hole counterparts, as might be expected for Z = As, Sb but counter to MEP trends for P. It is worth emphasizing, however, that the differences in energy between the two kinds of complexes were less than 0.5 kcal/mol. It might be concluded then that σ-hole and π-hole sites can compete with one another in these systems. Further analyses (NBO, Eint decomposition) revealed that the formation of the complexes is grounded not only in the electrostatic term (which nevertheless governs) but also charge transfer between subunits which varied from 0.3 up to 14.4 kcal/mol, for LP(N) → σ*(C-Sb) donation. The scale of charge transfer was larger for σ-hole dimers than in the case of π-hole adducts.
Another work by this group [157] involved the F2C=CFTF3 (T = C, Si, and Ge) as Lewis acid with water, ammonia, or formaldehyde as Lewis base, comparing the σ or π-holes on the tetrel T atoms. The σ-hole was formed along the T-F bond while the π-hole appeared perpendicular to the C=C double bond, coplanar to the σ-hole. The π-hole maxima rose along with the T atom size, in contrast to the case in earlier studies. However, the magnitudes were quite limited and for T = C and Si, the π-hole intensities were very close to one another. One can assume that the π-hole connected with the C=C bond is much less sensitive to the remainder of the molecule. The MEPs of the σ-holes were mostly consistent with previous patterns and also increased with the enlargement of T. The π-hole was stronger than the σ-hole in F2C=CFCF3 while the reverse was observed for T = Si and Ge. In keeping with this MEP data, F2C=CFCF3 prefers to form a π-hole tetrel bond whereas F2C=CFSiF3 and F2C=CFGeF3 were prone to a σ-hole tetrel bond. The strong correlation between MEP and interaction energy of the dimers was confirmed.

5. Conclusions and Prospective

Molecular systems can contain both σ- and π-holes, and there is a growing literature on both of them that allow some comparisons to be drawn. The presence of both sorts of holes is reflected in an expansion of possible binding sites and an increased flexibility in the noncovalent bonding with nucleophiles. The examination of crystals displays this flexibility in an array of binding patterns [104,140,141,142,143,144,145,146,149,153,155] where the directionality of σ/π-hole bonds is integral to arrangement, stabilization, and self-organization. There are numerous examples where these interactions are important to systems with biological connections [144,149,155] or photophysical performance [141,143,145].
In addition to the crystallographic data and their supporting theoretical investigations, examinations of systems in solution or model geometries in the gas phase provide additional insights into mutual presence of both σ and π-holes. The most-commonly observed situation to date is the combination of a σ-hole on a halogen atom with a π-hole from a variety of sources, most notably an aromatic ring. However, the full list is rather extensive, as for example when both sorts of hole occur on a single tetrel atom. One interesting conclusion is that a particular sort of hole does not have to be present in the isolated monomer. For example, the approach of a nucleophile can induce geometric deformation into the Lewis acid which in turn causes the appearance of a π-hole that is not present in its isolated form.
Another pattern which has emerged from the studies is that a σ-hole typically deepens as the atom on which it occurs grows larger, e.g., Si < Ge < Sn, whereas this same enlarging atom can lead to a weakening of a π-hole. These opposing trends offer the opportunity to guide an emerging complex geometry based on atom size. When both hole types are present, they each offer an attractive site for binding by a nucleophile. However, since both bonding types utilize the Lewis acid as electron acceptor, the two bonds weaken one another in a negative cooperative manner.
Whether σ or π-hole type, the noncovalent bond appears to have a strong electrostatic component, based on rigorous calculations of this component, as well as a tight correlation of bond strength with the depth of the hole. However, these bonds rely to a large extent also on polarization and charge transfer effects, complemented by dispersive forces, although the precise mix of these different components varies from one bond to another. A slightly different perspective on some of these bonds has been offered by the Bickelhaupt group [22,23] in which emphasis is placed on the activation strain model of chemical reactivity and the energy decomposition analysis combined with molecular orbital theory.
It is hoped that this review of the current state of knowledge concerning these σ and π-hole bonds will motivate additional work to better refine our understanding of the forces that undergird and control them. In a practical sense, this enhanced knowledge base will hopefully lead to the development of new crystal packing motifs and to innovative materials with improved properties.

Author Contributions

Conceptualization, W.Z., M.M., and S.S.; software, writing—original draft preparation, W.Z. and M.M.; writing—review and editing, W.Z., M.M., and S.S; graphic materials, S. S, W.Z., and M.M.; funding acquisitions, W.Z. and S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded partially by the Polish Ministry of Science and Higher Education, grant number for the Faculty of Chemistry of Wroclaw University of Science and Technology 8211104160/K19W03D10, and by the U.S. National Science Foundation under Grant No. 1954310.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

A generous grant of computer time from the Wroclaw Supercomputer and Networking Center is acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Clark, T.; Hennemann, M.; Murray, J.S.; Politzer, P. Halogen bonding: The σ-hole: Proceedings of “Modeling interactions in biomolecules II”, Prague, September 5th–9th, 2005. J. Mol. Model. 2007, 13, 291–296. [Google Scholar] [CrossRef]
  2. Guthrie, F. On the Iodide of Iodammonium. J. Chem. Soc. 1863, 16, 239–244. [Google Scholar] [CrossRef] [Green Version]
  3. Remsen, I.; Norris, J.F. Action of the halogens on the methyl-amines. Am. Chem. J. 1896, 18, 90–96. [Google Scholar]
  4. Hope, H.; McCullough, J.D. The crystal structure of the molecular complex of iodine with tetrahydroselenophene, C4H8Se.I2. Acta Cryst. 1964, 17, 712–718. [Google Scholar] [CrossRef]
  5. Olie, K.; Mijlhoff, F.C. The crystal structure of POBr3 and intermolecular bonding. Acta Cryst. Sect. B 1969, 25, 974–977. [Google Scholar] [CrossRef]
  6. Brinck, T.; Murray, J.S.; Politzer, P. Surface electrostatic potentials of halogenated methanes as indicators of directional intermolecular interactions. Int. J. Quantum Chem. 1992, 44, 57–64. [Google Scholar] [CrossRef]
  7. Murray, J.S.; Lane, P.; Clark, T.; Politzer, P. Sigma-hole bonding: Molecules containing group VI atoms. J. Mol. Model. 2007, 13, 1033–1038. [Google Scholar] [CrossRef]
  8. Politzer, P.; Lane, P.; Concha, M.C.; Ma, Y.G.; Murray, J.S. An overview of halogen bonding. J. Mol. Model. 2007, 13, 305–311. [Google Scholar] [CrossRef] [PubMed]
  9. Politzer, P.; Murray, J.S.; Concha, M.C. Sigma-hole bonding between like atoms; a fallacy of atomic charges. J. Mol. Model. 2008, 14, 659–665. [Google Scholar] [CrossRef] [PubMed]
  10. Murray, J.S.; Lane, P.; Politzer, P. Expansion of the sigma-hole concept. J. Mol. Model. 2009, 15, 723–729. [Google Scholar] [CrossRef] [PubMed]
  11. Politzer, P.; Murray, J.S.; Clark, T. Halogen bonding: An electrostatically-driven highly directional noncovalent interaction. Phys. Chem. Chem. Phys. 2010, 12, 7748–7757. [Google Scholar] [CrossRef] [PubMed]
  12. Murray, J.S.; Politzer, P. The electrostatic potential: An overview. Wires Comput. Mol. Sci. 2011, 1, 153–163. [Google Scholar] [CrossRef]
  13. Murray, J.S.; Politzer, P. Interaction and polarization energy relationships in σ-hole and π-hole bonding. Crystals 2020, 10, 76. [Google Scholar] [CrossRef] [Green Version]
  14. Politzer, P.; Murray, J.S. Halogen Bonding: An Interim Discussion. ChemPhysChem 2013, 14, 278–294. [Google Scholar] [CrossRef] [PubMed]
  15. Politzer, P.; Murray, J.S.; Clark, T. sigma-Hole Bonding: A Physical Interpretation. Top Curr. Chem. 2015, 358, 19–42. [Google Scholar] [PubMed]
  16. Politzer, P.; Murray, J.S.; Clark, T.; Resnati, G. The sigma-hole revisited. Phys. Chem. Chem. Phys. 2017, 19, 32166–32178. [Google Scholar] [CrossRef]
  17. Politzer, P.; Murray, J.S. Electrostatics and Polarization in sigma- and pi-Hole Noncovalent Interactions: An Overview. ChemPhysChem 2020, 21, 579–588. [Google Scholar] [CrossRef]
  18. Azofra, L.M.; Alkorta, I.; Scheiner, S. Chalcogen bonds in complexes of SOXY (X, Y = F, Cl) with nitrogen bases. J. Phys. Chem. A 2015, 119, 535–541. [Google Scholar] [CrossRef] [Green Version]
  19. Vincent De Paul, N.N.; Scheiner, S. Chalcogen bonding between tetravalent SF4 and amines. J. Phys. Chem. A 2014, 118, 10849–10856. [Google Scholar]
  20. Alikhani, E.; Fuster, F.; Madebene, B.; Grabowski, S.J. Topological reaction sites-very strong chalcogen bonds. Phys. Chem. Chem. Phys. 2014, 16, 2430–2442. [Google Scholar] [CrossRef]
  21. Wang, W.; Ji, B.; Zhang, Y. Chalcogen bond: A sister noncovalent bond to halogen bond. J. Phys. Chem. A 2009, 113, 8132–8135. [Google Scholar] [CrossRef]
  22. de Azevedo Santos, L.; Ramalho, T.C.; Hamlin, T.A.; Bickelhaupt, F.M. Chalcogen bonds: Hierarchical ab initio benchmark and density functional theory performance study. J. Comput. Chem. 2021, 42, 688–698. [Google Scholar] [CrossRef]
  23. de Azevedo Santos, L.; van der Lubbe, S.C.C.; Hamlin, T.A.; Ramalho, T.C.; Matthias Bickelhaupt, F. A Quantitative Molecular Orbital Perspective of the Chalcogen Bond. ChemistryOpen 2021. [Google Scholar] [CrossRef]
  24. Alkorta, I.; Elguero, J.; Del Bene, J.E. Pnicogen Bonded Complexes of PO2X (X = F, Cl) with Nitrogen Bases. J. Phys. Chem. A 2013, 117, 10497–10503. [Google Scholar] [CrossRef] [PubMed]
  25. Del Bene, J.E.; Alkorta, I.; Elguero, J. Properties of Complexes H2C=(X)P:PXH2, for X = F, Cl, OH, CN, NC, CCH, H, CH3, and BH2: P center dot center dot center dot P Pnicogen Bonding at sigma-Holes and pi-Holes. J. Phys. Chem. A 2013, 117, 11592–11604. [Google Scholar] [CrossRef]
  26. Scheiner, S. Detailed comparison of the pnicogen bond with chalcogen, halogen, and hydrogen bonds. Int. J. Quantum Chem. 2013, 113, 1609–1620. [Google Scholar] [CrossRef] [Green Version]
  27. Scheiner, S. The pnicogen bond: Its relation to hydrogen, halogen, and other noncovalent bonds. Acc. Chem. Res. 2013, 46, 280–288. [Google Scholar] [CrossRef] [PubMed]
  28. Guan, L.Y.; Mo, Y.R. Electron Transfer in Pnicogen Bonds. J. Phys. Chem. A 2014, 118, 8911–8921. [Google Scholar] [CrossRef] [PubMed]
  29. Bauza, A.; Mooibroek, T.J.; Frontera, A. Tetrel-Bonding Interaction: Rediscovered Supramolecular Force? Angew. Chem. Int. Ed. 2013, 52, 12317–12321. [Google Scholar] [CrossRef] [PubMed]
  30. Bauza, A.; Mooibroek, T.J.; Frontera, A. Tetrel Bonding Interactions. Chem. Rec. 2016, 16, 473–487. [Google Scholar] [CrossRef]
  31. Marin-Luna, M.; Alkorta, I.; Elguero, J. Cooperativity in Tetrel Bonds. J. Phys. Chem. A 2016, 120, 648–656. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Desiraju, G.R.; Shing Ho, P.; Kloo, L.; Legon, A.C.; Marquardt, R.; Metrangolo, P.; Politzer, P.; Resnati, G.; Rissanen, K. Definition of the halogen bond (IUPAC recommendations 2013). Pure Appl. Chem. 2013, 85, 1711–1713. [Google Scholar] [CrossRef]
  33. Aakeroy, C.B.; Bryce, D.L.; Desiraju, G.; Frontera, A.; Legon, A.C.; Nicotra, F.; Rissanen, K.; Scheiner, S.; Terraneo, G.; Metrangolo, P.; et al. Definition of the chalcogen bond (IUPAC Recommendations 2019). Pure Appl. Chem. 2019, 91, 1889–1892. [Google Scholar] [CrossRef]
  34. Mahmoudi, G.; Afkhami, F.A.; Zangrando, E.; Kaminsky, W.; Frontera, A.; Safin, D.A. A supramolecular 3D structure constructed from a new metal chelate self-assembled from Sn(NCS)(2) and phenyl(pyridin-2-yl)methylenepicolinohydrazide. J. Mol. Struct. 2021, 1224, 129188. [Google Scholar] [CrossRef]
  35. Zhang, J.H.; Su, Z.Z.; Luo, J.X.; Zhao, Y.; Wang, H.G.; Ying, S.M. Synthesis, structure, and characterization of a mixed amines thiogermanate [NH4](2)[NH2(CH3)(2)](2)Ge2S6. Polyhedron 2020, 182, 114486. [Google Scholar] [CrossRef]
  36. Zelenkov, L.E.; Ivanov, D.M.; Sadykov, E.K.; Bokach, N.A.; Galmes, B.; Frontera, A.; Kukushkin, V.Y. Semicoordination Bond Breaking and Halogen Bond Making Change the Supramolecular Architecture of Metal-Containing Aggregates. Cryst. Growth Des. 2020, 20, 6956–6965. [Google Scholar] [CrossRef]
  37. Yusof, E.N.M.; Latif, M.A.M.; Tahir, M.I.M.; Sakoff, J.A.; Veerakumarasivam, A.; Page, A.J.; Tiekink, E.R.T.; Ravoof, T.B.S.A. Homoleptic tin(IV) compounds containing tridentate ONS dithiocarbazate Schiff bases: Synthesis, X-ray crystallography, DFT and cytotoxicity studies. J. Mol. Struct. 2020, 1205, 127635. [Google Scholar] [CrossRef]
  38. Walton, I.; Chen, C.; Rimsza, J.M.; Nenoff, T.M.; Walton, K.S. Enhanced Sulfur Dioxide Adsorption in UiO-66 Through Crystal Engineering and Chalcogen Bonding. Cryst. Growth Des. 2020, 20, 6139–6146. [Google Scholar] [CrossRef]
  39. Venosova, B.; Koziskova, J.; Kozisek, J.; Herich, P.; Luspai, K.; Petricek, V.; Hartung, J.; Muller, M.; Hubschle, C.B.; van Smaalen, S.; et al. Charge density of 4-methyl-3-[(tetrahydro-2H-pyran-2-yl)oxy]thiazole-2(3H)-thione. A comprehensive multipole refinement, maximum entropy method and density functional theory study. Acta Cryst. B 2020, 76, 450–468. [Google Scholar] [CrossRef]
  40. Villasenor-Granados, T.O.; Montes-Tolentino, P.; Rodriguez-Lopez, G.; Sanchez-Ruiz, S.A.; Flores-Parra, A. Structural analysis of tris(5-methyl- [1,3,5]-dithiazinan-2-yl)stibine, its reactions with chalcogens. Intramolecular chalcogen-bonding interactions. J. Mol. Struct. 2020, 1200, 127050. [Google Scholar] [CrossRef]
  41. Kumar, V.; Xu, Y.J.; Bryce, D.L. Double Chalcogen Bonds: Crystal Engineering Stratagems via Diffraction and Multinuclear Solid-State Magnetic Resonance Spectroscopy. Chem. Eur. J. 2020, 26, 3275–3286. [Google Scholar] [CrossRef]
  42. Fourmigue, M.; Dhaka, A. Chalcogen bonding in crystalline diselenides and selenocyanates: From molecules of pharmaceutical interest to conducting materials. Coord. Chem. Rev. 2020, 403, 213084. [Google Scholar] [CrossRef]
  43. Michalczyk, M.; Malik, M.; Zierkiewicz, W.; Scheiner, S. Experimental and Theoretical Studies of Dimers Stabilized by Two Chalcogen Bonds in the Presence of a N···N Pnicogen Bond. J. Phys. Chem. A 2021, 125, 657–668. [Google Scholar] [CrossRef] [PubMed]
  44. Kinzhalov, M.A.; Popova, E.A.; Petrov, M.L.; Khoroshilova, O.V.; Mahmudov, K.T.; Pombeiro, A.J.L. Pnicogen and chalcogen bonds in cyclometalated iridium(III) complexes. Inorg. Chim. Acta 2018, 477, 31–33. [Google Scholar] [CrossRef]
  45. Zapata, F.; Gonzalez, L.; Bastida, A.; Bautista, D.; Caballero, A. Formation of self-assembled supramolecular polymers by anti-electrostatic anion-anion and halogen bonding interactions. Chem. Commun. 2020, 56, 7084–7087. [Google Scholar] [CrossRef] [PubMed]
  46. Soldatova, N.S.; Postnikov, P.S.; Suslonov, V.V.; Kissler, T.Y.; Ivanov, D.M.; Yusubov, M.S.; Galmes, B.; Frontera, A.; Kukushkin, V.Y. Diaryliodonium as a double sigma-hole donor: The dichotomy of thiocyanate halogen bonding provides divergent solid state arylation by diaryliodonium cations. Org. Chem. Front. 2020, 7, 2230–2242. [Google Scholar] [CrossRef]
  47. Scilabra, P.; Terraneo, G.; Daolio, A.; Baggioli, A.; Famulari, A.; Leroy, C.; Bryce, D.L.; Resnati, G. 4,4 ‘-Dipyridyl Dioxide center dot SbF3 Cocrystal: Pnictogen Bond Prevails over Halogen and Hydrogen Bonds in Driving Self-Assembly. Cryst. Growth Des. 2020, 20, 916–922. [Google Scholar] [CrossRef] [Green Version]
  48. Nascimento, M.A.; Heyer, E.; Less, R.J.; Pask, C.M.; Arauzo, A.; Campo, J.; Rawson, J.M. An Investigation of Halogen Bonding as a Structure-Directing Interaction in Dithiadiazolyl Radicals. Cryst. Growth Des. 2020, 20, 4313–4324. [Google Scholar] [CrossRef]
  49. Pramanik, S.; Chopra, D. Unravelling the Importance of H bonds, σ–hole and π–hole-Directed Intermolecular Interactions in Nature. J. Indian I Sci. 2020, 100, 43–59. [Google Scholar] [CrossRef]
  50. Orlova, A.P.; Jasien, P.G. Halogen bonding in self-assembling systems: A comparison of intra- and interchain binding energies. Comput. Chem. 2018, 1139, 63–69. [Google Scholar] [CrossRef]
  51. Tiekink, E.R.T. A Survey of Supramolecular Aggregation Based on Main Group Element Selenium Secondary Bonding Interactions—A Survey of the Crystallographic Literature. Crystals 2020, 10, 503. [Google Scholar] [CrossRef]
  52. Seth, S.K.; Bauzá, A.; Frontera, A. Chapter 9 Quantitative Analysis of Weak Non-covalent σ-Hole and π-Hole Interactions. In Monogr Supramol Chem; The Royal Society of Chemistry: London, UK, 2019; pp. 285–333. [Google Scholar]
  53. Taylor, M.S. Anion recognition based on halogen, chalcogen, pnictogen and tetrel bonding. Coord. Chem. Rev. 2020, 413, 213270. [Google Scholar] [CrossRef]
  54. Navarro-Garcia, E.; Galmes, B.; Velasco, M.D.; Frontera, A.; Caballero, A. Anion Recognition by Neutral Chalcogen Bonding Receptors: Experimental and Theoretical Investigations. Chem. Eur. J. 2020, 26, 4706–4713. [Google Scholar] [CrossRef] [PubMed]
  55. Borissov, A.; Marques, I.; Lim, J.Y.C.; Felix, V.; Smith, M.D.; Beer, P.D. Anion Recognition in Water by Charge-Neutral Halogen and Chalcogen Bonding Foldamer Receptors. J. Am. Chem. Soc. 2019, 141, 4119–4129. [Google Scholar] [CrossRef] [PubMed]
  56. Semenov, N.A.; Gorbunov, D.E.; Shakhova, M.V.; Salnikov, G.E.; Bagryanskaya, I.Y.; Korolev, V.V.; Beckmann, J.; Gritsan, N.P.; Zibarev, A.V. Donor-Acceptor Complexes between 1,2,5-Chalcogenadiazoles (Te, Se, S) and the Pseudohalides CN- and XCN- (X = O, S, Se, Te). Chem. Eur. J. 2018, 24, 12983–12991. [Google Scholar] [CrossRef]
  57. Naseer, M.M.; Hussain, M.; Bauza, A.; Lo, K.M.; Frontera, A. Intramolecular Noncovalent Carbon Bonding Interaction Stabilizes the cis Conformation in Acylhydrazones. Chempluschem 2018, 83, 881–885. [Google Scholar] [CrossRef] [PubMed]
  58. Lim, J.Y.C.; Beer, P.D. Sigma-Hole Interactions in Anion Recognition. Chem-Us 2018, 4, 731–783. [Google Scholar] [CrossRef]
  59. Liu, Y.-Z.; Yuan, K.; Liu, L.; Yuan, Z.; Zhu, Y.-C. Anion Recognition Based on a Combination of Double-Dentate Hydrogen Bond and Double-Side Anion−π Noncovalent Interactions. J. Phys. Chem. A 2017, 121, 892–900. [Google Scholar] [CrossRef]
  60. Lim, J.Y.C.; Marques, I.; Thompson, A.L.; Christensen, K.E.; Felix, V.; Beer, P.D. Chalcogen Bonding Macrocycles and [2]Rotaxanes for Anion Recognition. J. Am. Chem. Soc. 2017, 139, 3122–3133. [Google Scholar] [CrossRef]
  61. Scheiner, S. Tetrel Bonding as a Vehicle for Strong and Selective Anion Binding. Molecules 2018, 23, 1147. [Google Scholar] [CrossRef] [Green Version]
  62. Kociok-Kohn, G.; Mahon, M.F.; Molloy, K.C.; Price, G.J.; Prior, T.J.; Smith, D.R.G. Biomimetic polyorganosiloxanes: Model compounds for new materials. Dalton Trans. 2014, 43, 7734–7746. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Politzer, P.; Murray, J.S.; Concha, M.C. Halogen bonding and the design of new materials: Organic bromides, chlorides and perhaps even fluorides as donors. J. Mol. Model. 2007, 13, 643–650. [Google Scholar] [CrossRef] [PubMed]
  64. Radovic, L.R. Probing the ‘elephant’: On the essential difference between graphenes and polycyclic aromatic hydrocarbons. Carbon 2021, 171, 798–805. [Google Scholar] [CrossRef]
  65. Maruthapandi, M.; Eswaran, L.; Cohen, R.; Perkas, N.; Luong, J.H.T.; Gedanken, A. Silica-Supported Nitrogen-Enriched Porous Benzimidazole-Linked and Triazine-Based Polymers for the Adsorption of CO2. Langmuir 2020, 36, 4280–4288. [Google Scholar] [CrossRef] [PubMed]
  66. Kar, I.; Chatterjee, J.; Harnagea, L.; Kushnirenko, Y.; Fedorov, A.V.; Shrivastava, D.; Buchner, B.; Mahadevan, P.; Thirupathaiah, S. Metal-chalcogen bond-length induced electronic phase transition from semiconductor to topological semimetal in ZrX2 (X = Se and Te). Phys. Rev. B 2020, 101, 165122. [Google Scholar] [CrossRef] [Green Version]
  67. Isaeva, A.; Ruck, M. Crystal Chemistry and Bonding Patterns of Bismuth-Based Topological Insulators. Inorg. Chem. 2020, 59, 3437–3451. [Google Scholar] [CrossRef] [PubMed]
  68. Ho, P.C.; Wang, J.Z.; Meloni, F.; Vargas-Baca, I. Chalcogen bonding in materials chemistry. Coord. Chem. Rev. 2020, 422, 213464. [Google Scholar] [CrossRef]
  69. Daniels, C.L.; Knobeloch, M.; Yox, P.; Adamson, M.A.S.; Chen, Y.H.; Dorn, R.W.; Wu, H.; Zhou, G.Q.; Fan, H.J.; Rossini, A.J.; et al. Intermetallic Nanocatalysts from Heterobimetallic Group 10-14 Pyridine-2-thiolate Precursors. Organometallics 2020, 39, 1092–1104. [Google Scholar] [CrossRef]
  70. Zhang, Y.; Wang, W.Z.; Wang, Y.B. Tetrel bonding on graphene. Comput. Chem. 2019, 1147, 8–12. [Google Scholar] [CrossRef]
  71. Sung, S.H.; Schnitzer, N.; Brown, L.; Park, J.; Hovden, R. Stacking, strain, and twist in 2D materials quantified by 3D electron diffraction. Phys. Rev. Mater. 2019, 3, 064003. [Google Scholar] [CrossRef] [Green Version]
  72. Mahmudov, K.T.; Kopylovich, M.N.; Guedes da Silva, M.F.C.; Pombeiro, A.J.L. Chalcogen bonding in synthesis, catalysis and design of materials. Dalton Trans. 2017, 46, 10121–10138. [Google Scholar] [CrossRef] [Green Version]
  73. Pina, M.D.N.; Frontera, A.; Bauza, A. Quantifying Intramolecular Halogen Bonds in Nucleic Acids: A Combined Protein Data Bank and Theoretical Study. Acs Chem. Biol. 2020, 15, 1942–1948. [Google Scholar] [CrossRef]
  74. Lundemba, A.S.; Bibelayi, D.D.; Wood, P.A.; Pradon, J.; Yav, Z.G. sigma-Hole interactions in small-molecule compounds containing divalent sulfur groups R-1-S-R-2. Acta Cryst. B 2020, 76, 707–718. [Google Scholar] [CrossRef] [PubMed]
  75. Scheiner, S. Differential Binding of Tetrel-Bonding Bipodal Receptors to Monatomic and Polyatomic Anions. Molecules 2019, 24, 227. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Riel, A.M.S.; Huynh, H.T.; Jeannin, O.; Berryman, O.; Fourmigue, M. Organic Selenocyanates as Halide Receptors: From Chelation to One-Dimensional Systems. Cryst. Growth Des. 2019, 19, 1418–1425. [Google Scholar] [CrossRef]
  77. Michalczyk, M.; Zierkiewicz, W.; Wysokinski, R.; Scheiner, S. Theoretical Studies of IR and NMR Spectral Changes Induced by Sigma-Hole Hydrogen, Halogen, Chalcogen, Pnicogen, and Tetrel Bonds in a Model Protein Environment. Molecules 2019, 24, 3329. [Google Scholar] [CrossRef] [Green Version]
  78. Wei, Y.; Li, Q.; Scheiner, S. The pi-Tetrel Bond and its Influence on Hydrogen Bonding and Proton Transfer. Chemphyschem 2018, 19, 736–743. [Google Scholar] [CrossRef]
  79. Trievel, R.C.; Scheiner, S. Crystallographic and Computational Characterization of Methyl Tetrel Bonding in S-Adenosylmethionine-Dependent Methyltransferases. Molecules 2018, 23, 2965. [Google Scholar] [CrossRef] [Green Version]
  80. Scheiner, S. Halogen Bonds Formed between Substituted Imidazoliums and N Bases of Varying N-Hybridization. Molecules 2017, 22, 1634. [Google Scholar] [CrossRef]
  81. Montaña, Á.M. The σ and π Holes. The Halogen and Tetrel Bondings: Their Nature, Importance and Chemical, Biological and Medicinal Implications. Chemistryselect 2017, 2, 9094–9112. [Google Scholar] [CrossRef]
  82. Grabowski, S.J. Tetrel bond-sigma-hole bond as a preliminary stage of the S(N)2 reaction. Phys. Chem. Chem. Phys. 2014, 16, 1824–1834. [Google Scholar] [CrossRef]
  83. Angarov, V.; Kozuch, S. On the σ, π and δ hole interactions: A molecular orbital overview. New J. Chem. 2018, 42, 1413–1422. [Google Scholar] [CrossRef]
  84. Bauza, A.; Mooibroek, T.J.; Frontera, A. The Bright Future of Unconventional sigma/-Hole Interactions. Chemphyschem 2015, 16, 2496–2517. [Google Scholar] [CrossRef]
  85. Gao, L.; Zeng, Y.; Zhang, X.; Meng, L. Comparative studies on group III sigma-hole and pi-hole interactions. J. Comput. Chem. 2016, 37, 1321–1327. [Google Scholar] [CrossRef] [PubMed]
  86. Bauza, A.; Mooibroek, T.J.; Frontera, A. Directionality of pi-holes in nitro compounds. Chem. Commun. 2015, 51, 1491–1493. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Murray, J.S.; Lane, P.; Clark, T.; Riley, K.E.; Politzer, P. sigma-Holes, pi-holes and electrostatically-driven interactions. J. Mol. Model. 2012, 18, 541–548. [Google Scholar] [CrossRef] [PubMed]
  88. Grabowski, S.J. The Nature of Triel Bonds, a Case of B and Al Centres Bonded with Electron Rich Sites. Molecules 2020, 25, 2703. [Google Scholar] [CrossRef]
  89. Grabowski, S.J. Triel bond and coordination of triel centres-Comparison with hydrogen bond interaction. Coord. Chem. Rev. 2020, 407, 213171. [Google Scholar] [CrossRef]
  90. Xu, Z.F.; Li, Y. Triel bonds in RZH(2)center dot center dot center dot NH3: Hybridization, solvation, and substitution. J. Mol. Model. 2019, 25, 219. [Google Scholar] [CrossRef]
  91. Jablonski, M. Hydride-Triel Bonds. J. Comput. Chem. 2018, 39, 1177–1191. [Google Scholar] [CrossRef]
  92. Grabowski, S.J. Two faces of triel bonds in boron trihalide complexes. J. Comput. Chem. 2018, 39, 472–480. [Google Scholar] [CrossRef]
  93. Esrafili, M.D.; Mousavian, P. The triel bond: A potential force for tuning anion-pi interactions. Mol. Phys. 2018, 116, 388–398. [Google Scholar] [CrossRef]
  94. Grabowski, S.J. Triel bonds-complexes of boron and aluminum trihalides and trihydrides with benzene. Struct. Chem. 2017, 28, 1163–1171. [Google Scholar] [CrossRef]
  95. Grabowski, S.J. pi-Hole Bonds: Boron and Aluminum Lewis Acid Centers. Chemphyschem 2015, 16, 1470–1479. [Google Scholar] [CrossRef] [PubMed]
  96. Burgi, H.B.; Dunitz, J.D.; Shefter, E. Geometrical reaction coordinates. II. Nucleophilic addition to a carbonyl group. J. Am. Chem. Soc. 1973, 95, 5065–5067. [Google Scholar] [CrossRef]
  97. Harder, M.; Kuhn, B.; Diederich, F. Efficient Stacking on Protein Amide Fragments. Chemmedchem 2013, 8, 397–404. [Google Scholar] [CrossRef]
  98. Bartlett, G.J.; Choudhary, A.; Raines, R.T.; Woolfson, D.N. n ->pi* interactions in proteins. Nat. Chem. Biol. 2010, 6, 615–620. [Google Scholar] [CrossRef]
  99. Andleeb, H.; Khan, I.; Bauza, A.; Tahir, M.N.; Simpson, J.; Hameed, S.; Frontera, A. Synthesis and supramolecular self-assembly of thioxothiazolidinone derivatives driven by pi-bonding and diverse 7 hole interactions: A combined experimental and theoretical analysis. J. Mol. Struct. 2017, 1139, 209–221. [Google Scholar] [CrossRef]
  100. Gomila, R.M.; Frontera, A. Covalent and Non-covalent Noble Gas Bonding Interactions in XeFn Derivatives (n = 2–6): A Combined Theoretical and ICSD Analysis. Front. Chem. 2020, 8, 395. [Google Scholar] [CrossRef]
  101. Bauza, A.; Frontera, A. sigma/pi-Hole noble gas bonding interactions: Insights from theory and experiment. Coord. Chem. Rev. 2019, 404, 213112. [Google Scholar] [CrossRef]
  102. Esrafili, M.D.; Vessally, E. The strengthening effect of a hydrogen or lithium bond on the Z···N aerogen bond (Z = Ar, Kr and Xe): A comparative study. Mol. Phys. 2016, 114, 3265–3276. [Google Scholar] [CrossRef]
  103. Bauza, A.; Frontera, A. pi-Hole aerogen bonding interactions. Phys. Chem. Chem. Phys. 2015, 17, 24748–24753. [Google Scholar] [CrossRef] [PubMed]
  104. Bauzá, A.; Frontera, A. Aerogen Bonding Interaction: A New Supramolecular Force? Angew. Chem. Int. Ed. 2015, 54, 7340–7343. [Google Scholar] [CrossRef] [PubMed]
  105. Alkorta, I.; Elguero, J.; Frontera, A. Not Only Hydrogen Bonds: Other Noncovalent Interactions. Crystals 2020, 10, 180. [Google Scholar] [CrossRef] [Green Version]
  106. Hou, M.C.; Zhu, Y.F.; Li, Q.Z.; Scheiner, S. Tuning the Competition between Hydrogen and Tetrel Bonds by a Magnesium Bond. Chemphyschem 2020, 21, 212–219. [Google Scholar] [CrossRef]
  107. Alkorta, I.; Legon, A.C. Non-Covalent Interactions Involving Alkaline-Earth Atoms and Lewis Bases B: An ab Initio Investigation of Beryllium and Magnesium Bonds, B center dot center dot center dot MR2 (M = Be or Mg, and R = H, F or CH3). Inorganics 2019, 7, 35. [Google Scholar] [CrossRef] [Green Version]
  108. Grabowski, S.J. Magnesium Bonds: From Divalent Mg Centres to Trigonal and Tetrahedral Coordination. Chemistryselect 2018, 3, 3147–3154. [Google Scholar] [CrossRef]
  109. Alikhani, M.E. Beryllium bonding: Insights from the sigma- and pi-hole analysis. J. Mol. Model. 2020, 26, 94. [Google Scholar] [CrossRef]
  110. Zhang, Z.; Lu, T.; Ding, L.Y.; Wang, G.Y.; Wang, Z.X.; Zheng, B.S.; Liu, Y.; Ding, X.L. Cooperativity effects between regium-bonding and pnicogen-bonding interactions in ternary MF center dot center dot center dot PH3O center dot center dot center dot MF (M = Cu, Ag, Au): An ab initio study. Mol. Phys. 2020, 118, e1784478. [Google Scholar] [CrossRef]
  111. Wang, R.J.; Wang, Z.; Yu, X.F.; Li, Q.Z. Synergistic and Diminutive Effects between Regium and Aerogen Bonds. Chemphyschem 2020, 21, 2426–2431. [Google Scholar] [CrossRef] [PubMed]
  112. Sanchez-Sanz, G.; Trujillo, C.; Alkorta, I.; Elguero, J. Rivalry between Regium and Hydrogen Bonds Established within Diatomic Coinage Molecules and Lewis Acids/Bases. Chemphyschem 2020, 21, 2557–2563. [Google Scholar] [CrossRef]
  113. Alkorta, I.; Trujillo, C.; Sanchez-Sanz, G.; Elguero, J. Regium Bonds between Silver(I) Pyrazolates Dinuclear Complexes and Lewis Bases (N-2, OH2, NCH, SH2, NH3, PH3, CO and CNH). Crystals 2020, 10, 137. [Google Scholar] [CrossRef] [Green Version]
  114. Sanchez-Sanz, G.; Trujillo, C.; Alkorta, I.; Elguero, J. Understanding Regium Bonds and their Competition with Hydrogen Bonds in Au-2:HX Complexes. Chemphyschem 2019, 20, 1572–1580. [Google Scholar] [CrossRef]
  115. Zierkiewicz, W.; Michalczyk, M.; Scheiner, S. Regium bonds between Mn clusters (M = Cu, Ag, Au and n = 2–6) and nucleophiles NH3 and HCN. Phys. Chem. Chem. Phys. 2018, 20, 22498–22509. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  116. Xia, T.; Li, D.; Cheng, L.J. Theoretical analysis of the spodium bonds in HgCl2 center dot center dot center dot L (L = ClR, SR2, and PR3) dimers. Chem. Phys. 2020, 539, 110978. [Google Scholar] [CrossRef]
  117. Mahmoudi, G.; Masoudiasl, A.; Babashkina, M.G.; Frontera, A.; Doert, T.; White, J.M.; Zangrando, E.; Zubkov, F.I.; Safin, D.A. On the importance of pi-hole spodium bonding in tricoordinated Hg-II complexes. Dalton T 2020, 49, 17547–17551. [Google Scholar] [CrossRef] [PubMed]
  118. Mahmoudi, G.; Lawrence, S.E.; Cisterna, J.; Cardenas, A.; Brito, I.; Frontera, A.; Safin, D.A. A new spodium bond driven coordination polymer constructed from mercury(ii) azide and 1,2-bis(pyridin-2-ylmethylene)hydrazine. New J. Chem. 2020, 44, 21100–21107. [Google Scholar] [CrossRef]
  119. Bauze, A.; Alkorta, I.; Elguero, J.; Mooibroek, T.J.; Frontera, A. Spodium Bonds: Noncovalent Interactions Involving Group 12 Elements. Angew. Chem. Int. Ed. 2020, 59, 17482–17487. [Google Scholar] [CrossRef]
  120. Frontera, A. σ- and π-Hole Interactions. Crystals 2020, 10, 721. [Google Scholar] [CrossRef]
  121. Cavallo, G.; Metrangolo, P.; Milani, R.; Pilati, T.; Priimagi, A.; Resnati, G.; Terraneo, G. The halogen bond. Chem. Rev. 2016, 116, 2478–2601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  122. Esrafili, M.D.; Mohammadirad, N. An ab initio study on tunability of σ-hole interactions in XHS:PH2Y and XH2P:SHY complexes (X = F, Cl, Br; Y = H, OH, OCH3, CH3, C2H5, and NH2). J. Mol. Model. 2015, 21, 2727. [Google Scholar] [CrossRef]
  123. Riley, K.E.; Murray, J.S.; Fanfrlik, J.; Rezac, J.; Sola, R.J.; Concha, M.C.; Ramos, F.M.; Politzer, P. Halogen bond tunability II: The varying roles of electrostatic and dispersion contributions to attraction in halogen bonds. J. Mol. Model. 2013, 19, 4651–4659. [Google Scholar] [CrossRef]
  124. Riley, K.E.; Murray, J.S.; Fanfrlik, J.; Rezac, J.; Sola, R.J.; Concha, M.C.; Ramos, F.M.; Politzer, P. Halogen bond tunability I: The effects of aromatic fluorine substitution on the strengths of halogen-bonding interactions involving chlorine, bromine, and iodine. J. Mol. Model. 2011, 17, 3309–3318. [Google Scholar] [CrossRef]
  125. Politzer, P.; Murray, J.S. sigma-Hole Interactions: Perspectives and Misconceptions. Crystals 2017, 7, 212. [Google Scholar] [CrossRef] [Green Version]
  126. Bulat, F.A.; Toro-Labbe, A.; Brinck, T.; Murray, J.S.; Politzer, P. Quantitative analysis of molecular surfaces: Areas, volumes, electrostatic potentials and average local ionization energies. J. Mol. Model. 2010, 16, 1679–1691. [Google Scholar] [CrossRef]
  127. Murray, J.S.; Politzer, P. Molecular electrostatic potentials and noncovalent interactions. Wires Comput. Mol. Sci. 2017, 7, e13260. [Google Scholar] [CrossRef]
  128. Politzer, P.; Murray, J.S. sigma-holes and pi-holes: Similarities and differences. J. Comput. Chem. 2018, 39, 464–471. [Google Scholar] [CrossRef]
  129. Clark, T.; Murray, J.S.; Politzer, P. Role of Polarization in Halogen Bonds. Aust. J. Chem. 2014, 67, 451–456. [Google Scholar] [CrossRef] [Green Version]
  130. Hennemann, M.; Murray, J.S.; Politzer, P.; Riley, K.E.; Clark, T. Polarization-induced sigma-holes and hydrogen bonding. J. Mol. Model. 2012, 18, 2461–2469. [Google Scholar] [CrossRef] [PubMed]
  131. Riley, K.E.; Hobza, P. The relative roles of electrostatics and dispersion in the stabilization of halogen bonds. Phys. Chem. Chem. Phys. 2013, 15, 17742–17751. [Google Scholar] [CrossRef] [PubMed]
  132. Zierkiewicz, W.; Michalczyk, M.; Scheiner, S. Aerogen bonds formed between AeOF(2) (Ae = Kr, Xe) and diazines: Comparisons between sigma-hole and pi-hole complexes. Phys. Chem. Chem. Phys. 2018, 20, 4676–4687. [Google Scholar] [CrossRef]
  133. Zierkiewicz, W.; Michalczyk, M.; Wysokinski, R.; Scheiner, S. On the ability of pnicogen atoms to engage in both sigma and pi-hole complexes. Heterodimers of ZF2C6H5 (Z = P, As, Sb, Bi) and NH3. J. Mol. Model. 2019, 25, 152. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  134. Wysokinski, R.; Zierkiewicz, W.; Michalczyk, M.; Scheiner, S. How Many Pnicogen Bonds can be Formed to a Central Atom Simultaneously? J. Phys. Chem. A 2020, 124, 2046–2056. [Google Scholar] [CrossRef] [PubMed]
  135. Zierkiewicz, W.; Wysokinski, R.; Michalczyk, M.; Scheiner, S. Chalcogen bonding of two ligands to hypervalent YF4 (Y = S, Se, Te, Po). Phys. Chem. Chem. Phys. 2019, 21, 20829–20839. [Google Scholar] [CrossRef] [PubMed]
  136. Michalczyk, M.; Zierkiewicz, W.; Wysokinski, R.; Scheiner, S. Hexacoordinated Tetrel-Bonded Complexes between TF4 (T =Si, Ge, Sn, Pb) and NCH: Competition between sigma- and pi-Holes. Chemphyschem 2019, 20, 959–966. [Google Scholar] [CrossRef] [Green Version]
  137. Makarewicz, E.; Lundell, J.; Gordon, A.J.; Berski, S. On the nature of interactions in the F2 OXe(...) NCCH3 complex: Is there the Xe(IV)N bond? J. Comput. Chem. 2016, 37, 1876–1886. [Google Scholar] [CrossRef]
  138. Brock, D.S.; Bilir, V.; Mercier, H.P.A.; Schrobilgen, G.J. XeOF2, F2OXeN equivalent to CCH3, and XeOF2 center dot nHF: Rare examples of Xe(IV) oxide fluorides. J. Am. Chem. Soc. 2007, 129, 3598–3611. [Google Scholar] [CrossRef]
  139. Brock, D.S.; de Pury, J.J.C.; Mercier, H.P.A.; Schrobilgen, G.J.; Silvi, B. A Rare Example of a Krypton Difluoride Coordination Compound: [BrOF2][AsF6]center dot 2KrF(2). J. Am. Chem. Soc. 2010, 132, 3533–3542. [Google Scholar] [CrossRef]
  140. Wang, H.; Wang, W.; Jin, W.J. σ-Hole Bond vs π-Hole Bond: A Comparison Based on Halogen Bond. Chem. Rev. 2016, 116, 5072–5104. [Google Scholar] [CrossRef]
  141. Zhang, Y.-H.; Li, Y.-L.; Yang, J.; Zhou, P.-P.; Xie, K. Noncovalent functionalization of graphene via π-hole···π and σ-hole···π interactions. Struct. Chem. 2020, 31, 97–101. [Google Scholar] [CrossRef]
  142. Yang, F.-L.; Yang, X.; Wu, R.-Z.; Yan, C.-X.; Yang, F.; Ye, W.; Zhang, L.-W.; Zhou, P.-P. Intermolecular interactions between σ- and π-holes of bromopentafluorobenzene and pyridine: Computational and experimental investigations. Phys. Chem. Chem. Phys. 2018, 20, 11386–11395. [Google Scholar] [CrossRef] [PubMed]
  143. Liu, Z.F.; Chen, X.; Wu, W.X.; Zhang, G.Q.; Li, X.; Li, Z.Z.; Jin, W.J. 1,3,5-Trifluoro-2,4,6-triiodobenzene: A neglected NIR phosphor with prolonged lifetime by σ-hole and π-hole capture. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2020, 224, 117428. [Google Scholar] [CrossRef]
  144. Shukla, R.; Khan, I.; Ibrar, A.; Simpson, J.; Chopra, D. Complex electronic interplay of σ-hole and π-hole interactions in crystals of halogen substituted 1,3,4-oxadiazol-2(3H)-thiones. Crystengcomm 2017, 19, 3485–3498. [Google Scholar] [CrossRef]
  145. Li, L.; Wang, H.; Wang, W.; Jin, W.J. Interactions between haloperfluorobenzenes and fluoranthene in luminescent cocrystals from π-hole⋯π to σ-hole⋯π bonds. Crystengcomm 2017, 19, 5058–5067. [Google Scholar] [CrossRef]
  146. Wang, J.; Mo, L.; Li, X.; Geng, Z.; Zeng, Y. The protonated 2-halogenated imidazolium cation as the noncovalent interaction donor: The σ-hole and π-hole interactions. J. Mol. Model. 2016, 22, 299. [Google Scholar] [CrossRef]
  147. Zhang, Y.; Ji, B.M.; Tian, A.M.; Wang, W.Z. Communication: Competition between pi center dot center dot center dot pi interaction and halogen bond in solution: A combined C-13 NMR and density functional theory study. J. Chem. Phys. 2012, 136, 141101. [Google Scholar] [CrossRef] [PubMed]
  148. Ma, N.; Zhang, Y.; Ji, B.M.; Tian, A.M.; Wang, W.Z. Structural Competition between Halogen Bonds and Lone-Pair center dot center dot center dot p Interactions in Solution. Chemphyschem 2012, 13, 1411–1414. [Google Scholar] [CrossRef] [PubMed]
  149. Shukla, R.; Claiser, N.; Souhassou, M.; Lecomte, C.; Balkrishna, S.J.; Kumar, S.; Chopra, D. Exploring the simultaneous [sigma]-hole/[pi]-hole bonding characteristics of a Br...[pi] interaction in an ebselen derivative via experimental and theoretical electron-density analysis. Iucrj 2018, 5, 647–653. [Google Scholar] [CrossRef] [Green Version]
  150. Lang, T.; Li, X.; Meng, L.; Zheng, S.; Zeng, Y. The cooperativity between the σ-hole and π-hole interactions in the ClO···XONO2/XONO···NH3 (X = Cl, Br, I) complexes. Struct. Chem. 2015, 26, 213–221. [Google Scholar] [CrossRef]
  151. Solimannejad, M.; Nassirinia, N.; Amani, S. A computational study of 1:1 and 1:2 complexes of nitryl halides (O2NX) with HCN and HNC. Struct. Chem. 2013, 24, 651–659. [Google Scholar] [CrossRef]
  152. Solimannejad, M.; Ramezani, V.; Trujillo, C.; Alkorta, I.; Sanchez-Sanz, G.; Elguero, J. Competition and Interplay between sigma-Hole and pi-Hole Interactions: A Computational Study of 1:1 and 1:2 Complexes of Nitryl Halides (O2NX) with Ammonia. J. Phys. Chem. A 2012, 116, 5199–5206. [Google Scholar] [CrossRef] [PubMed]
  153. McDowell, S.A.C.; Joseph, J.A. A comparative study of model halogen-bonded, π-hole-bonded and cationic complexes involving NCX AND H2O (X = F, Cl, Br). Mol. Phys. 2015, 113, 16–21. [Google Scholar] [CrossRef]
  154. Bauza, A.; Frontera, A. On the Versatility of BH2X (X = F, Cl, Br, and I) Compounds as Halogen-, Hydrogen-, and Triel-Bond Donors: An AbInitio Study. Chemphyschem 2016, 17, 3181–3186. [Google Scholar] [CrossRef] [PubMed]
  155. Pal, R.; Nagendra, G.; Samarasimhareddy, M.; Sureshbabu, V.V.; Guru Row, T.N. Observation of a reversible isomorphous phase transition and an interplay of “σ-holes” and “π-holes” in Fmoc-Leu-ψ[CH2-NCS]. Chem. Commun. 2015, 51, 933–936. [Google Scholar] [CrossRef]
  156. Dong, W.; Wang, Y.; Cheng, J.; Yang, X.; Li, Q. Competition between σ-hole pnicogen bond and π-hole tetrel bond in complexes of CF2=CFZH2 (Z = P, As, and Sb). Mol. Phys. 2019, 117, 251–259. [Google Scholar] [CrossRef]
  157. Dong, W.; Yang, X.; Cheng, J.; Li, W.; Li, Q. Comparison for σ-hole and π-hole tetrel-bonded complexes involving F2CCFTF3 (TC, Si, and Ge): Substitution, hybridization, and solvation effects. J. Fluor. Chem. 2018, 207, 38–44. [Google Scholar] [CrossRef]
Figure 1. Family of σ and π-hole interactions.
Figure 1. Family of σ and π-hole interactions.
Molecules 26 01740 g001
Figure 2. Ge-O NBO bonding and antibonding orbitals of GeF2O. Red and blue colors indicate opposite phase of the wave function.
Figure 2. Ge-O NBO bonding and antibonding orbitals of GeF2O. Red and blue colors indicate opposite phase of the wave function.
Molecules 26 01740 g002
Figure 3. Electron density (10−5 au) and molecular electrostatic potential (MEP) (10−3 au) of GeF2O as a function of the displacement from the molecular plane θ, at a distance of 2.5 Å from the Ge atom.
Figure 3. Electron density (10−5 au) and molecular electrostatic potential (MEP) (10−3 au) of GeF2O as a function of the displacement from the molecular plane θ, at a distance of 2.5 Å from the Ge atom.
Molecules 26 01740 g003
Figure 4. (a) MEP on surface of GeF2O defined as 1.5 × vdW atomic radii. Blue and red colors refer to +0.05 and −0.05 au, respectively. (b) Optimized complex of GeF2O with NH3 at the mp2/aug-cc-pVDZ level, with distance in Å.
Figure 4. (a) MEP on surface of GeF2O defined as 1.5 × vdW atomic radii. Blue and red colors refer to +0.05 and −0.05 au, respectively. (b) Optimized complex of GeF2O with NH3 at the mp2/aug-cc-pVDZ level, with distance in Å.
Molecules 26 01740 g004
Figure 5. MEP of the isolated KrOF2 and XeOF2 molecules on the 0.001 au contour of the electron isodensity, at the MP2/aug-cc-pVDZ level. Color ranges, in kcal/mol, are red, greater than 40, yellow; between 20 and 40, green; between 0 and 20, blue, less than 0 (negative). Selected surface critical points Vs,max (σ and π-holes) are indicated as black dots. Reproduced from Reference [132] with permission from the PCCP Owner Societies.
Figure 5. MEP of the isolated KrOF2 and XeOF2 molecules on the 0.001 au contour of the electron isodensity, at the MP2/aug-cc-pVDZ level. Color ranges, in kcal/mol, are red, greater than 40, yellow; between 20 and 40, green; between 0 and 20, blue, less than 0 (negative). Selected surface critical points Vs,max (σ and π-holes) are indicated as black dots. Reproduced from Reference [132] with permission from the PCCP Owner Societies.
Molecules 26 01740 g005
Figure 6. Two binding modes of XeOF2···NCCH3 complexes, from Reference [137].
Figure 6. Two binding modes of XeOF2···NCCH3 complexes, from Reference [137].
Molecules 26 01740 g006
Figure 7. Two types of the MP2/aug-cc-pVDZ optimized structures (top and side views) of complexes of NH3 with ZF2C6H5, from Reference [133].
Figure 7. Two types of the MP2/aug-cc-pVDZ optimized structures (top and side views) of complexes of NH3 with ZF2C6H5, from Reference [133].
Molecules 26 01740 g007
Figure 8. Cis and Trans arrangements of two NH3 units complexed with YF4 via chalcogen bonds. Reproduced from Reference [135] with permission from the PCCP Owner Societies.
Figure 8. Cis and Trans arrangements of two NH3 units complexed with YF4 via chalcogen bonds. Reproduced from Reference [135] with permission from the PCCP Owner Societies.
Molecules 26 01740 g008
Figure 9. Optimized structures of NCF complexes with water stabilized by (a) σ-hole and (b) π-hole, from Reference [153].
Figure 9. Optimized structures of NCF complexes with water stabilized by (a) σ-hole and (b) π-hole, from Reference [153].
Molecules 26 01740 g009
Table 1. Comparison between σ and π hole depths in molecules where they coexist on the same atom, VS,max in kcal/mol.
Table 1. Comparison between σ and π hole depths in molecules where they coexist on the same atom, VS,max in kcal/mol.
MoleculeHole SourceHole TypeVS,maxReferences
Aerogen Bond Donors
KrOF2Krσ58.7[132]
KrOF2Krπ39.1[132]
XeOF2Xeσ63.4[132]
XeOF2Xeσ90[104]
XeOF2Xeπ36.2[132]
Pnicogen Bond Donors
PF2C6H5Pσ19.4[133]
PF2C6H5Pπ36.2*
AsF2C6H5Asσ28.2[133]
AsF2C6H5Asπ44.0*
SbF2C6H5Sbσ38.4[133]
SbF2C6H5Sbπ56.6*
BiF2C6H5Biσ52.6[133]
BiF2C6H5Biπ60.9*
PF3Pσ35.6[134]
PF3Pπ9.7[134]
AsF3Asσ43.9[134]
AsF3Asπ7.1[134]
SbF3Sbσ51.6[134]
SbF3Sbπ10.6[134]
BiF3Biσ61.5[134]
BiF3Biπ12.7[134]
Chalcogen Bond Donors
SF4Sσ41.7[135]
SF4Sπ64.4[135] **
SeF4Seσ51.2[135]
SeF4Seπ61.1[135] **
TeF4Teσ59.2[135]
TeF4Teπ54.5[135] **
PoF4Poσ76.3[135]
PoF4Poπ53.2[135] **
Tetrel Bond Donors
SiF4Siσ127.3[136]
SiF4Siπ109.8[136]
GeF4Geσ120.9[136]
GeF4Geπ106.1[136]
SnF4Snσ129.4[136]
SnF4Snπ121.3[136]
PbF4Pbσ127.3[136]
PbF4Pbπ104.1[136]
* Value obtained by additional calculations for the purpose of the current review, at the same level as in Reference [136]. ** Values obtained for monomer in deformed complex geometry.
Table 3. Comparison between σ and π holes appearing on the same molecule, where σ-hole is coupled with halogen atom and π-hole with different sources (other than aromatic ring), VS,max in kcal/mol.
Table 3. Comparison between σ and π holes appearing on the same molecule, where σ-hole is coupled with halogen atom and π-hole with different sources (other than aromatic ring), VS,max in kcal/mol.
MoleculeHole SourceHole TypeVS,maxReferences
Halogen-Pnicogen
XONO (X = F, Cl, Br, I)Xσ−6.3 to 51.1[150]
XONO (X = F, Cl, Br, I)Nπ20.8 to 29.5[150]
XONO2 (X = F, Cl, Br, I)Xσ2.7 to 67.7[150]
XONO2 (X = F, Cl, Br, I)Nπ29.6 to 41.2[150]
NO2X (X = Cl, Br)Xσ13.2 and 19.0[151,152]
NO2X (X = Cl, Br)Nπ28.1 and 29.8[151,152]
NO2IIσ29.4[152]
NO2INπ23.5[152]
Halogen-Tetrel
NCX (X = F, Cl, Br)Xσ14.3 to 42.1[153]
NCX (X = F, Cl, Br)Cπ12.4 to 27.3[153]
Halogen-Triel
BH2X (X = F, Cl, Br, I)Cl, Br, Iσ3.5 to 11.3[154]
BH2X (X = F, Cl, Br, I)Bπ28.8 to 39.4[154]
Halogen-Beryllium
BeCl2Clσ1.2[109]
BeCl2Beπ32.2[109]
Table 4. Comparison between other combinations of σ and π holes appearing on the same molecule, VS,max in kcal/mol.
Table 4. Comparison between other combinations of σ and π holes appearing on the same molecule, VS,max in kcal/mol.
MoleculeHole SourceHole TypeVS,maxReferences
Chalcogen-Tetrel/Pnicogen
Fmoc-Leu-Ψ[CH2NCS]Sσ3.9 (exp.), 7.8 (theory)[155]
Fmoc-Leu-Ψ[CH2NCS]C=N bondπno data[155]
Pnicogen-Tetrel
CF2=CFZH2 (Z = P, As, Sb)Zσ19.4 to 28.9[156]
CF2=CFPF2Cπ36.4[156]
CF2=CFZH2 (Z = P, As, Sb)Zσ25.7 to 28.2[156]
CF2=CFPF2Cπ40.1[156]
Tetrel-Tetrel
F2C=CFTF3 (T = C, Si, Ge)Tσ8.2 to 46.4 *[157]
F2C=CFTF3 (T = C, Si, Ge)C=C Bondπ30.7 to 35.1 *[157]
* In cited work, these values are probably given in wrong unit (eV instead of more reliable au).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zierkiewicz, W.; Michalczyk, M.; Scheiner, S. Noncovalent Bonds through Sigma and Pi-Hole Located on the Same Molecule. Guiding Principles and Comparisons. Molecules 2021, 26, 1740. https://doi.org/10.3390/molecules26061740

AMA Style

Zierkiewicz W, Michalczyk M, Scheiner S. Noncovalent Bonds through Sigma and Pi-Hole Located on the Same Molecule. Guiding Principles and Comparisons. Molecules. 2021; 26(6):1740. https://doi.org/10.3390/molecules26061740

Chicago/Turabian Style

Zierkiewicz, Wiktor, Mariusz Michalczyk, and Steve Scheiner. 2021. "Noncovalent Bonds through Sigma and Pi-Hole Located on the Same Molecule. Guiding Principles and Comparisons" Molecules 26, no. 6: 1740. https://doi.org/10.3390/molecules26061740

APA Style

Zierkiewicz, W., Michalczyk, M., & Scheiner, S. (2021). Noncovalent Bonds through Sigma and Pi-Hole Located on the Same Molecule. Guiding Principles and Comparisons. Molecules, 26(6), 1740. https://doi.org/10.3390/molecules26061740

Article Metrics

Back to TopTop