The Stibium Bond or the Antimony-Centered Pnictogen Bond: The Covalently Bound Antimony Atom in Molecular Entities in Crystal Lattices as a Pnictogen Bond Donor

A stibium bond, i.e., a non-covalent interaction formed by covalently or coordinately bound antimony, occurs in chemical systems when there is evidence of a net attractive interaction between the electrophilic region associated with an antimony atom and a nucleophile in another, or the same molecular entity. This is a pnictogen bond and are likely formed by the elements of the pnictogen family, Group 15, of the periodic table, and is an inter- or intra-molecular non-covalent interaction. This overview describes a set of illustrative crystal systems that were stabilized (at least partially) by means of stibium bonds, together with other non-covalent interactions (such as hydrogen bonds and halogen bonds), retrieved from either the Cambridge Structure Database (CSD) or the Inorganic Crystal Structure Database (ICSD). We demonstrate that these databases contain hundreds of crystal structures of various dimensions in which covalently or coordinately bound antimony atoms in molecular entities feature positive sites that productively interact with various Lewis bases containing O, N, F, Cl, Br, and I atoms in the same or different molecular entities, leading to the formation of stibium bonds, and hence, being partially responsible for the stability of the crystals. The geometric features, pro-molecular charge density isosurface topologies, and extrema of the molecular electrostatic potential model were collectively examined in some instances to illustrate the presence of Sb-centered pnictogen bonding in the representative crystal systems considered.


Introduction
Naming a chemical interaction formed by an element of the periodic table and its subsequent characterization readily enable one to identify its status and visualize what its role could be in the rational design of complex systems [1,2]. After identifying, in the last century, the possible existence of "close contacts" between atomic domains in molecules, molecular complexes and crystals, considerable efforts were made to classify, characterize, and name them [2]. One such close contact that occurs between atomic domains is widely known as the hydrogen bond; it is a non-covalent (chemical) interaction caused by a force of attraction between atomic sites of opposite polarity, following the fundamental law of electrostatics, i.e., Coulomb's law, that states that opposite charges attract whereas like charges repel each other when in close proximity. As Arunan pointed out in an unpublished study [3], over fifty different proposals have been put forward to date on the characterization of hydrogen bonds, dating back to the time of Werner (1902) [4], Hantzsch (1910) [5] and Pfeifer (1913) [6], among many others [3]. However, a new definition of hydrogen bonds was adopted as recently as 2011 [7], listing many characteristics and indicators of 2 of 52 such interactions that emerged from a variety of computational and experimental measurements, including IR, Raman, and NMR spectroscopic techniques, and X-ray and neutron diffraction measurements.
First-principles [8][9][10][11][12] and density functional theory [13][14][15] methods have also been developed and applied to calculate the energy of hydrogen bonds in chemical systems and various spectroscopic and geometrical signatures, as well as the nature of the charge density profile associated with the surfaces of the interacting atomic domains in molecular entities responsible for these interactions.
Following the definition and delineation of the characteristic features of hydrogen bonding, a similar definition, together with a set of characteristic features, was recommended for halogen bonding in 2013 [16]. This appears to have been prompted by the realization that thousands of solid-state systems were synthesized during the last decade that conceived directional intermolecular interactions driven by covalently bound elements of Group 17. There is very little difference between the definitions provided for hydrogen bonds and halogen bonds (except that the word "hydrogen" replaces "halogen" in the latter, of course); both atom types in molecular entities must feature electrophilic regions on their electrostatic surfaces and be capable of accepting some fractional charge density from an electron density donor fragment in the partner molecular species when in its close proximity.
Since the definition and characteristic features of a hydrogen bond were found to be transferable to the elements of Group 17 of the periodic table, they are, in principle, transferable to other elements, for instance, those of Groups 14-16, as well. In fact, this was pointed out in an article published in 2017 that recommended a definition for and described a number of characteristic features of chalcogen bonds in chemical systems (Group 16) [17].
We believe that the same definition is also transferable to the family of pnictogen bonds formed by the pnictogen elements, i.e., Group 15. A pnictogen bond in chemical systems occurs when there is evidence of a net attractive interaction between the electrophilic region associated with a covalently or coordinately bound pnictogen atom in a molecular entity and a nucleophile in another, or the same molecular entity. It is formed by the elements of the pnictogen family, Group 15, of the periodic table, and is an inter-or intra-molecular non-covalent interaction.
Clearly, when a covalently bound nitrogen atom in a molecular entity behaves as an electrophile toward a Lewis base, one would be tempted to call such an attractive engagement a nitrogen bond, in which, the nitrogen atom is a pnictogen bond donor. The same conceptual framework is transferable to covalently bound phosphorous, arsenic, antimony, bismuth and (in principle, anyway) moscovium when they display positive sites and have the ability to donate the electrophile to a Lewis base to form phosphorous, arsenic, stibium, bismuth, and moscovium bonds, respectively. This overview is focused on the possible occurrence of antimony-centered pnictogen bonds (or simply, stibium bonds) in chemical systems that have been known since the last century. In our view, it is an overlooked non-covalent interaction yet to be fully appreciated by, among others, computational and supramolecular chemists, gas-phase spectroscopists and materials scientists. This is not surprising given that the term "pnictogen bond" was coined only recently [18], and that the majority of studies of non-covalent interactions in the last and current centuries have been focused on arriving at a basic understanding of hydrogen bonding and van der Waals interactions. Legon has argued that pnictogen bonding (and other interactions) was recognized long before it was specifically named [19].
For the reasons given in the following sections, we believe that supramolecular pnictogen bonds are one of the cornerstones for the development of novel functional materials [51][52][53] and occur simultaneously with other non-covalent interactions (viz. halogen bonding, tetrel bonding, and hydrogen bonding), although their importance may only be fully appreciated in the years to come. In this overview, we provide several illustrative examples of crystal structures deposited in the Inorganic Crystal Structure Database (ICSD) [54,55]) and the Cambridge Structural Database (CSD) [56] in which the identification and characterization of stibium bonds were apparently overlooked, even though they play a fundamental role in the assembly of these solid-state materials. This study is therefore expected to assist researchers in the future design of novel functional materials featuring stibium bonds. With this in mind, we attempt to highlight, using illustrative examples, the occurrence and modes of interaction of such bonds.

Antimony in Molecular Entities, Materials Design and Discovery
Antimony, the fourth member of the pnictogen family, has played an important role in the development of phase change and optoelectronic materials [57][58][59]. In particular, chalcogenides based on germanium-antimony-tellurides (GST-PCMs), such as GaSbTe and GaGe 2 Sb 2 Te 5 , display outstanding properties that are a prerequisite for the development of non-volatile memory (NVM) technologies due to their high write and read speeds, reversible phase transition, high degree of scalability, low power consumption, good data retention, and multi-level storage capability [60]. Kao and workers discussed the importance of some Sb-based alloys, viz. Ga 25 Te 8 Sb 67 and Ga 18 Te 12 Sb 70 , that have a crystallization temperature above 245 • C and activation energy of crystallization greater than 5 eV, and are useful for phase-change memory applications [61]. When combined with halides, they form inter alia various single and double perovskites which have been the focus of many studies because of their application in the development of photovoltaics and other areas of optoelectronics. A 3 Sb 2 X 9 and A 2 AgSbX 6 are a specific class of such perovskites. A a B b X a+3b are called rudorfites (A = methylammonium, formamidinium, Cs, Rb, etc.; B = Bi, Sb, etc.; X = I, Br, Cl), and several have been crystallized [59,62,63].
The quantum dots of the all-inorganic Sb-based lead-free double perovskite Cs 2 AgSbX 6 (X = Cl, Br or I) are examples of an air stable compound that displays strong blue emission with photoluminescence quantum yields of 31.33% [64]. The mixed metal 〈111〉-oriented layered perovskite, Cs 4 CuSb 2 Cl 12 , incorporates Cu 2+ and Sb 3+ into layers that are three octahedra thick. The compound is photo-and thermal-stable, tolerant of humidity and behaves as a semiconductor with a direct bandgap of 1.0 eV. Its conductivity is an order of magnitude greater than that of the widely studied MAPbI 3 (MA = methylammonium) perovskite [65].
Apart from the above, many crystal structures have been reported in which antimony plays a prominent role in the development of their structural stability and functionality. It appears that in many cases that covalently bound Sb can act as a pnictogen bond donor to form non-covalent interactions, a driving force for self-assembly. Our search and analysis of the structures deposited in the ICSD and CSD databases suggest that most of the interand intra-molecular bonding interactions emerge upon attraction between bound Sb in molecules and halides in the interacting partner entities. The intermolecular bond distances associated with them are often less than, but sometimes marginally greater than, the sum of the van der Waals (vdW) radii of the bound atomic basins. Primary and secondary noncovalent interactions (viz. hydrogen bonding and halogen bonding) are often present, and their interplay with pnictogen bonding drive the packing and hence the overall geometry of the crystal lattice.

Inter-or Intramolecular Bond Distance and Less Than the Sum of the van der Waals Radii Concept
Our attempts in rationalizing the possible occurrence of a stibium bond hinges on a geometric criterion: whether the inter-or intramolecular distance between covalently/coordinately bound pnictogen atom Pn (viz. Pn = Sb) and a Lewis base(s) D in a crystal lattice, Pn···D (D = an electron density donor, such as O, N, a halogen anion, etc.), is less than (or even slightly greater than) the sum of the van der Waals radii of the bound atomic basins. If the first attribute is met, we recognize the Pn···D link as a likely stibium bond, consequent on the overlapping of atomic domains takes place. In the case of the other attribute, care has to be exercised, since there are numerous occasions in which the inter-or intra-molecular distance, Pn···D, is larger than the sum of the vdW radii of the non-bonded atoms that are in close proximity. In this case, we have taken into consideration crystal systems where the Pn···D distance can exceed the vdW radii sum by several tenths of an ångstrom. This is quite justifiable according to Politzer and coworkers [66], as well as others [67,68], given that the values proposed for the vdW radii of atoms have a typical uncertainty of ±0.2 Å; hence "less than the sum of the vdW radii" concept will necessarily miss a significant number of non-covalent interactions if treated as a strict criterion to identify a non-covalent interaction in molecular entities. The reason is that hard sphere models with spherical symmetry of atoms were considered in the determination of vdW radii of atoms, which differ from one study to the other depending on the approximations invoked in their determination. The underlying reason of course is that the charge density profile of atoms in molecules is anisotropic, and hence the "vdW radius" of an atom in a molecular entity is likely to vary between molecular entities. Table 1, for instance, lists the vdW radii of some selected atoms that we invoke throughout this overview, showing that the Bondi's radius for a given atom type is different from those proposed by others. Since the specific vdW radius of an atom in different molecular entities is unknown, we use the same radius of an atomic domain to examine the "less than sum of the vdW radius concept", regardless of the lattice systems examined. While "an interatomic distance that is greater than the sum of the respective vdW radii by even several tenths of an ångstrom may still correspond to a non-covalent interaction" [66], caution needs to be exercised to validate the interaction by simultaneously examining attributes like its directional features and the nature of the polarity of the bound atomic basins. Otherwise, such an assignment could be misleading. In our discussion below, we use the vdW radii of atoms proposed recently by Alvarez [68].

Directionality of Inter-/Intra-Molecular Interactions
To maximize the integrity of our identification and subsequent characterization of pnictogen bonding in the examined crystal systems, we simultaneously examined the angle of interaction between the noncovalently bound atomic basins. We did so to see whether the angle θ of approach of the electrophile, θ = ∠R-Pn···D, associated with the bound Sb atom in a molecular entity forms either a linear, a quasilinear or a bent (non-linear) interaction with the atomic domain in the partner molecular entity. At the same time, we inspected whether the interacting atomic domain D was electrophilic or nucleophilic. As such, Type-I interactions (Scheme 1a), which can be further classified as (left) Type-Ia (trans) and (right) Type-Ib (cis), appear when the regions of interacting atomic domains in molecular entities both carry either positive or negative local polarity, with θ 1 = ∠R-Pn···D < 120 • and θ 2 = ∠R -D···Pn < 120 • for the former and 90 • < θ 1 < ∠R-Pn···D < 150 • and 90 • < θ 2 < ∠R -D···Pn < 150 • for the latter interactions, where R and R are the remaining part of the bound Sb atom in a molecular entity forms either a linear, a quasilinear or a bent (nonlinear) interaction with the atomic domain in the partner molecular entity. At the same time, we inspected whether the interacting atomic domain D was electrophilic or nucleophilic. As such, Type-I interactions (Scheme 1a), which can be further classified as (left) Type-Ia (trans) and (right) Type-Ib (cis), appear when the regions of interacting atomic domains in molecular entities both carry either positive or negative local polarity, with θ1 = ∠R-Pn···D < 120° and θ2 = ∠R′-D···Pn < 120° for the former and 90° < θ1 < ∠R-Pn···D < 150° and 90° < θ2 < ∠R′-D···Pn < 150° for the latter interactions, where R and R′ are the remaining part of the molecular entities associated with Pn and D, respectively. These are non-linear interactions, and hence, should be regarded as non-directional. Scheme 1. Illustration of (a) Type-I, (b) Type-II and (c) Type-III topologies (geometric) of noncovalent bonding interactions. θ's represents the angle of interaction between the interacting atomic domains, and Pn, D, R and R' refer to the covalently/coordinately bound pnictogen atom, the interacting atomic domain (generally nucleophilic), and the remaining part of the molecular entities associated with Pn and D atomic basins, respectively. δ± signifies the local polarity (positive or negative), and the small region on atom Pn along R-Pn σ-bond extension colored in green indicates an electron deficient region (a σ-hole). Scheme 1. Illustration of (a) Type-I, (b) Type-II and (c) Type-III topologies (geometric) of noncovalent bonding interactions. θ's represents the angle of interaction between the interacting atomic domains, and Pn, D, R and R' refer to the covalently/coordinately bound pnictogen atom, the interacting atomic domain (generally nucleophilic), and the remaining part of the molecular entities associated with Pn and D atomic basins, respectively. δ± signifies the local polarity (positive or negative), and the small region on atom Pn along R-Pn σ-bond extension colored in green indicates an electron deficient region (a σ-hole).
Directional interactions are linear or quasi-linear, and thus, generally follow a Type-IIa topology of bonding. Hence, when the angle of approach θ = ∠R-Pn···D = 180 • , we call the interaction linear, and when it deviates from linearity such that 150 • < θ < 180 • (Scheme 1b, left), we refer to the interaction as quasi-linear. In either case, the electrostatic surface of covalently/coordinately bound atom Pn conceives an electrophilic region along the extension of the R-Pn covalent or coordinate bond, and D is a Lewis base (such as N in NH 3 , O in OH 2 , and F in HF). However, when 90 • < θ < 150 • , and Pn conceives a positive site, we recognize the interaction it forms as Type-IIb (Scheme 1b, right). Both Type-IIa and Type-IIb interactions are of coulombic origin.
Type-III interactions (Scheme 1c) originate when the angle of interaction follow a Type-IIa topology of bonding, yet a portion of the interacting atomic domains are either both positive or both negative polarities. This classification scheme has been discussed elsewhere [22,49,50,72].

The σ-Hole and π-Hole Concepts, and Their Relationship with Pnictogen Bonds
A σ-hole is defined as a region of charge density deficiency, compared to the remainder of the bound atom, on the electrostatic surface of atom A that appears along the outer extension of the R-A covalent bond [73,74]. It can be positive, or negative, or even neutral [22,75]. For instance, a negative σ-hole can be observed on the electrostatic surface of atom F along the H-F bond extension of a hydrogen fluoride molecule [76], while a positive σ-hole can be seen on the surface of hydrogen halides HX (X = Cl, Br, I) [72,77]. The strength of the σ-hole on A is determined by the electronegativity and electron-withdrawing abilities of R, as well as the electronegativity and polarizability of A. Specifically, iodine in HI has a stronger σ-hole on its electrostatic surface compared to that of X in HX (X = F, Cl, Br) [72,78]. Similarly, the strength of the positive σ-hole on F in CN-F should be weaker than those of CN-X (X = Cl, Br, I), and the strength varies in the order CN-F < CN-Cl < CN-Br < CN-I. This trend is also true for systems with an arene moiety, for example, the strength of the positive σ-hole on F in C-F of C 6 F 6 is weaker than those of C-X (X = Cl, Br, I) in C 6 X 6 , and varies in the order C-F < C-Cl < C-Br < C-I [72].
A σ-hole centered inter-or intra-molecular interaction involving a positive σ-hole on the pnictogen atom along the R-Pn bond extension should be referred to simply as a pnictogen bond [50]. Hypervalent atoms in molecules have more than one σ-hole. For instance, the P in phosphorous trihalides PX 3 (P = F, Cl, Br, I) possesses three σ-holes, each along the extension of the X-P bond [50]. Similarly, C in CX 4 (X = F, Cl, Br, I) possesses four σ-holes, each along the extension of each of the four X-C bonds [22,73,76,79].
On the other hand, an inter-or intra-molecular interaction may be regarded as a π-hole centered pnictogen bond when the pnictogen atom in molecular entities features a π-hole on its electrostatic surface and has the ability to engage attractively with a negative site in a neighboring molecule, or a site that has an electron density different to that of the π-hole, thus providing stability to the geometry of the resulting structure. A π-hole may be defined as an electron deficient region on the surface of a molecular entity; it may appear on the electrostatic surface (centroid region) of a delocalized arene moiety (e.g., C 6 H 6 , C 6 X 6 ), on the surface of an atom (e.g., N in NO 3 − ) or on the central portion of a delocalized bond in a molecular entity (e.g., C≡C and C≡N bonds in HCCH and HCN, respectively). So, analogous with the σ-hole, a π-hole can be positive or negative.
A σ-hole interaction in a chemical system is generally observed to be directional, whereas a π-hole interaction is non-directional.

Model Systems and Computational Approaches
As indicated already above, the crystal systems presented in the following sections were retrieved either from the CSD [56] or the ICSD [54,55]. We show, in some cases, that the bound antimony atom in these molecular entities has a positive site, and therefore, is capable of interacting with a negative site of a partner molecular entity. This phenomenon is either wholly or partially responsible for the overall geometric architecture and stability of the resulting crystal lattice.
The characterization of the positive and negative sites on the surface of molecular entities was made possible by analyzing the electrostatic potentials. We did so by energyminimizing the geometry of some selected isolated molecular entities in the gas phase. For this, we used the Gaussian 16 program package [80]. The Density Functional Theory at the ωB97XD [81] level and Møller-Plesset's second-order perturbation theory (MP2) [82] were used. Depending on our interest, and for reasons discussed below, basis sets such as Aug-CC-pVTZ or def2-TZVPD were chosen. These were obtained from the basis set exchange library [83,84]. The calculation of the MESP [85] was done using the wavefunctions generated on the fully relaxed geometries of the monomeric entities. For this, AIMAll [86] and Multiwfn [87] codes were used. Geometry analysis, drawing of molecular entities/crystals, and isosurface plots were performed using the Mercury 4.0 [88]/Gaussview 5.0 [89] and VMD [90] suite of programs, respectively.
The isoelectron density envelope on which to calculate the electrostatic potential is quite arbitrary [23,75,77,[91][92][93] and can be computed at any contour level, such as 0.001, or 0.0015, or 0.002 a.u. (electrons bohr −3 ) of the total electronic density function ρ(r). The choice of a 0.001 or 0.002 a.u. envelope was suggested by Bader and co-workers [94,95] and others [96] because it was felt that such an envelope encompasses at least 95% of the electronic charge and should yield physically reasonable molecular "dimensions". For fluorinated systems, an isodensity envelope greater than 0.001 a.u. is recommended, since this particular envelope often fails to provide insight into the correct nature of surface potential extrema [22,49,50].
The promolecular charge density based isosurfaces were calculated within the promolecular framework of Independent Gradient Model (IGM) [97,98]. We did so because we were interested in the identification and subsequent characterization of Sb-centered pnictogen bonding within or between interacting monomer entities responsible for the crystals investigated, and to cross-check the reliability of these interactions emanated using intermolecular bond distances, directionality, and MESP model descriptions. Within the framework of IGM, promolecular atomic electron densities are summed, and the associated atomic gradients do not interfere. This can be achieved by using absolute values to sum the atomic gradients and by rejecting any electron gradient contragradience feature. The resulting total gradient |∇ρ IGM | is an upper limit of the true gradient, and the difference between them, δg, quantifies the net electron density gradient collapse due to interactions. Ultimately, this means that the δg descriptor identifies the presence of opposite signs in the components of the total electron density gradient |∇ρ(r)| due to interactions. IGM thus has the capacity to automatically separate intra-and inter-fragment interactions in a molecular entity, and that this can be plotted in 2D or in 3D (isosurfaces) to reveal the presence of inter-or intramolecular interactions between interacting atomic basins. The 3D shape of the isosurface volumes can be utilized to infer localized or delocalized interactions between interacting domains. The colors of these isosurface volumes, in blue and green, represent strong and weak attractions, respectively, while red represents a repulsive interaction. Hereafter, we refer to this approach as IGM-δg.

The Molecular Electrostatic Potential and Characterization of σand π-Holes in Molecular Entities
The extrema of potential on the surfaces of molecular entities appear in two different flavors [27][28][29][30]99,100], i.e., the local most minimum of potential, V S,min , and the local most maximum of potential, V S,max . They may be positive, or negative, or neutral in specific regions on the surface of a molecular entity depending on the extent of electron density depletion or electron density accumulation. Lone-pair regions on the surface of a covalently bound atom, or on the surface of an anion, should feature a V S,min , since such regions have high electron densities.
The positive and negative signs of V S,min or V S,max generally refer to regions that are electrophilic and nucleophilic, respectively. So an electrophilic region can be recognized when V S,min > 0 or V S,max > 0. Similarly, a nucleophilic region can be identified when V S,min < 0 or V S,max < 0. The magnitude of V S,min , or V S,max determines the strength of the potential; the larger the magnitude of V S,min , or V S,max , the stronger the nucleophilicity, or electrophilicity, of the region concerned. For instance, V S,max was observed to be positive and larger on the surface of the I than on X (X = F, Cl, Br) along the extension of the H-X and C-X covalent bonds in HX and CX 4 molecules; hence the electrophilicity of covalently bound I along the H-I and C-I bond extension is the stronger when compared to that of Br, Cl, F in each of the two families of molecular entities. Since V S,max on X in these molecules appears along the extension of the H-X/C-X σ-bonds, it describes the sign and magnitude of the σ-hole on X. The concept is transferable to other covalently bound atoms; thus a positive or negative σ-hole can be found on the surfaces of a covalently bound triel, tetrel, chalcogen, pnictogen and halogen atoms, among others, depending on whether V S,max is positive or negative along the outer extension of the covalent bond.
An attractive coulombic interaction between two atomic basins (intermolecular, or intramolecular) may be recognized when a region of an atom or fragment in a molecular entity with a positive V S,min (or V S,max ) is in close proximity to that with a negative V S,min (or V S,max ) on the same, or a neighboring, molecular entity.
A σ-hole interaction between two atomic basins (intermolecular, or intramolecular) may be recognized when a positive V S,max on atom A along the R-A bond extension is in close proximity to a negative site described by a negative V S,min or V S,max . When atom A = Pn conceives a positive V S,max , the attractive engagement is a σ-hole centered pnictogen bound interaction, or simply a pnictogen bond.
Similarly, a π-hole centered intermolecular or intramolecular interaction between two atomic basins may be recognized when a positive V S,min or V S,max centered on the p-orbital of atom, a bond, a delocalized ring, or on a fragment, in a molecular entity is able to attract a negative site described by a negative V S,min or V S,max on another similar, or different, molecular entity. When a positive π-hole is found on a pnictogen atom in a molecular entity, and if it is a capable for attracting a negative site on the same, or a different molecule, the resulting interaction between them is a consequence of π-centered pnictogen bonding.
There are many reports of σ-hole and π-hole interactions forming the basis of the stability of numerous chemical systems. Some of them, for instance, include discussions about their similarities and differences, as well as the controversies and misconceptions surrounding this concept [51,73,[101][102][103][104][105][106].

Illustrative Crystal Systems
As we show below, we also observed that Sb in molecular entities can cause the formation of pseudo-covalent bonds when it finds itself in close proximity to a negative site in another molecular species. It has the capability to accept electron density from lone-pairs centered on O, N, F, Cl, Br, and I in other molecules. Antimony trihalides SbX 3 (X = F, Cl, Br, I) are probably the most simplest, best examples that feature the ability of Sb to form non-covalent interactions with negative sites on a neighboring molecule (vide infra). The ability of covalently bound Sb to form very strong pnictogen bonds upon its attractive engagement with halide anions has been recognized, with the energy of these interactions reaching the lower bound of covalent bonding energies [107]. It was argued that such interactions comprise significant covalent character, even though these were sometimes thought to be essentially electrostatic interactions.

Antimony Trihalides
In the various examples provided in Figure 1, we mainly focus on highlighting the intermolecular bonding modes of SbX 3 and a few other compounds containing Sb. As can be seen, the SbX 3 molecules have a trigonal pyramidal molecular geometry. In order to understand the chemical reactivity of the various constituents of SbX 3 that are responsible for the intermolecular interactions in the crystals, we focused on examining the nature of the electrophilic and nucleophilic regions on the electrostatic surfaces of isolated SbX 3 molecules. We did so by exploring the MESP at the MP2(full)/def2-TZVPPD and ωB97XD/def2-TZVPPD levels of theory. Our results for the local maxima and minima of potential on the electrostatic surface of the SbX 3 molecules are shown in Table 2. Figure 2a-d shows the MESP graphs for the isolated SbX 3 molecules.
There are cases in which the 5s electron pair appears to have no stereochemical influence. Thus, in (NH 4 ) 2 Sb 2 Br 12 , Sb 3+ is O h and the Br-Sb-Br angles are insignificantly different from 90 • , while the Sb-Br bonds are significantly lengthened [125]. This may be analogus with phosphorus(III) compounds of the type Ar 2 PX (Ar = 2,5-bis(trifluoromethyl)phenyl, X = Br, Cl) that feature a stereochemically inactive lone pair, with P 3+ displaying an approximately octahedral geometry [126].
The electron density associated with the stereochemically active lone pair in NX 3 and PX 3 were previously reported from evidence of their corresponding MESP plots [49,50], but this is not apparent in the MESP of SbX 3 ( Figure 2). This does not mean that the lone pair is not stereochemically active; it is there, but what is evident is the gradual polarization of charge density from Sb to X as the softness of Sb compared to N [49] and P [50] increases (and its electronegativity decreases) down the group.
Nevertheless, from Table 2, it can be seen that F in SbF 3 is entirely negative and the σ-hole on it is completely neutralized. The trend in the stability of the σ-holes along the Sb-X bond extensions in SbX 3 increases with an increase in the polarizability of X: Sb-Cl (1.6 kcal mol −1 ) < Sb-Br (6.5 kcal mol −1 ) < Sb-I (12.3 kcal mol −1 ) (Figure 2a-d).
The lone pair regions on X in the same molecules, described by V S,min , become less negative as one passes from the lightest to the heaviest member of the halogen family. The most positive regions on the surfaces of the SbX 3 molecules are identified on Sb along the three X-Sb bond extensions and that the stability of Sb's σ-hole follows the order SbF 3 > SbCl 3 > SbBr 3 > SbI 3 . The central region of the triangular face formed by the three halogens is described by a weakly positive V S,max , and is surrounded by negative potentials. These features of the electrostatic potential, elucidated by ab initio MP2(full) calculations, are very similar to those that arise from DFT-based ωB97XD calculations (cf. Table 2).
The existence of the electrophilic sites on Sb and the nucleophilic sites on X in SbX 3 (X = Cl, Br) explains that the attractive interactions that led to the development of the Sb···X intermolecular bonding interactions that were observed in the crystal geometries shown in Figure 1 were coulombic in nature.

Tetramethyl-Antimony Iodide
Tetramethyl-antimony iodide salt [(CH 3 ) 4 Sb] + [I] − is a good example of a solid state structure featuring charge-assisted pnictogen bonding in a sterically crowded environment; see Figure 3a [127]. We did not perform a MESP calculation for the system, since the polarity of the molecular domains responsible for the crystal lattice is quite apparent. The methyl hydrogen atoms of the three nearest (CH 3 ) 4 Sb + moieties of the crystal are non-covalently bound with iodide anions through charge-assisted C-H···I hydrogen bonds. Six such prominent interactions are formed with I − . They are all equivalent and the angle of attack for the formation for each of these and the corresponding intermolecular distances are ∠C-H···I = 147.4 • and 3.222 Å, respectively. Apart from this, [(CH 3 ) 4 Sb]I features chargeassisted Type-IIa C-Sb···I pnictogen bonds (r(Sb···I) = 4.056 and 4.352 Å) and Sb-C···I tetrel bonds (r(C···I) = 4.276 Å) that are significantly longer than the hydrogen bonds, and are directional interactions (∠C-Sb···I = 180 • and ∠Sb-C···I = 180 • ). The long intermolecular distances are not unexpected, since the sum of the vdW radii for C and Sb (r(C) vdW = 1.77 Å and r(Sb) vdW = 2.47 Å) is larger than that of I and H (r(I) vdW = 2.04 Å and r(H) vdW = 1.20 Å). the polarity of the molecular domains responsible for the crystal lattice is quite apparent. The methyl hydrogen atoms of the three nearest (CH3)4Sb + moieties of the crystal are noncovalently bound with iodide anions through charge-assisted C-H···I hydrogen bonds. Six such prominent interactions are formed with I − . They are all equivalent and the angle of attack for the formation for each of these and the corresponding intermolecular distances are ∠C-H···I = 147.4° and 3.222 Å, respectively. Apart from this, [(CH3)4Sb]I features charge-assisted Type-IIa C-Sb···I pnictogen bonds (r(Sb···I) = 4.056 and 4.352 Å) and Sb-C···I tetrel bonds (r(C···I) = 4.276 Å) that are significantly longer than the hydrogen bonds, and are directional interactions (∠C-Sb···I = 180° and ∠Sb-C···I = 180°). The long intermolecular distances are not unexpected, since the sum of the vdW radii for C and Sb (r(C)vdW = 1.77 Å and r(Sb)vdW = 2.47 Å) is larger than that of I and H (r(I)vdW = 2.04 Å and r(H)vdW = 1.20 Å). Our observation is in line with that of the authors of [127], in that the coordination around antimony is distorted tetragonal pyramidal, and the tetrahedral coordination of the ions in the structure results in the formation of a wurtzite-type structure with antimony-iodine distances of 4.06 and 4.35 Å. Although the non-covalent bonding features that occur in the crystal constitute σ-hole interactions, the involvement of the σ-holes on Our observation is in line with that of the authors of [127], in that the coordination around antimony is distorted tetragonal pyramidal, and the tetrahedral coordination of the ions in the structure results in the formation of a wurtzite-type structure with antimonyiodine distances of 4.06 and 4.35 Å. Although the non-covalent bonding features that occur in the crystal constitute σ-hole interactions, the involvement of the σ-holes on H and Sb in the long-range interactions cause the formation of a layer-type molecular framework in the crystal (Figure 3b).
Our contention that C-H···I hydrogen bonds and C-Sb···I pnictogen bonds exist in the extended crystal is given credence by the IGM-δg based isosurface plots shown in Figure 3c, obtained at three different isovalues. The circular volume (green) in the graph on the left occurs between Sb and I atomic basins, whereas the isosurfaces that occur between the two molecular entities are illustrated in the central and right side of the graph. These features suggest possible attractive interactions between bound atomic basins.

The Co-Crystal of Antimony Trihalide and Molecular Sulfur
The addition compound SbI 3 :3S 8 crystallizes in the R3m space group [128]. The structure is stabilized by many non-covalent interactions. As noted in the study, each iodine atom of the SbI 3 unit is attached to a S atom of an S 8 moiety with a Sb-I···S distance of 3.60 Å and an ∠Sb-I···S angle of 169.4 • . Each I-bound Sb also links with three nearest neighbor S sites of S 8 crown-shaped molecules, forming Sb-I···S intermolecular contacts, with I-Sb···S = 3.391 Å and an ∠I-Sb···S angle of 152.5 • . These bonding features can be inferred from Figure 4. framework in the crystal (Figure 3b).
Our contention that C-H···I hydrogen bonds and C-Sb···I pnictogen bonds exist in the extended crystal is given credence by the IGM-δg based isosurface plots shown in Figure 3c, obtained at three different isovalues. The circular volume (green) in the graph on the left occurs between Sb and I atomic basins, whereas the isosurfaces that occur between the two molecular entities are illustrated in the central and right side of the graph. These features suggest possible attractive interactions between bound atomic basins.

The Co-crystal of Antimony Trihalide and Molecular Sulfur
The addition compound SbI3:3S8 crystallizes in the R3m space group [128]. The structure is stabilized by many non-covalent interactions. As noted in the study, each iodine atom of the SbI3 unit is attached to a S atom of an S8 moiety with a Sb-I···S distance of 3.60 Å and an ∠Sb-I···S angle of 169.4°. Each I-bound Sb also links with three nearest neighbor S sites of S8 crown-shaped molecules, forming Sb-I···S intermolecular contacts, with I-Sb···S = 3.391 Å and an ∠I-Sb···S angle of 152.5°. These bonding features can be inferred from Figure 4. Our calculation of the MESP of SbI3 and S8 at the ωB97XD/Aug-CC-pVTZ level of theory provides a more in-depth understanding of the bonding modes in this compound than that of the original study. Our calculation of the MESP of SbI 3 and S 8 at the ωB97XD/Aug-CC-pVTZ level of theory provides a more in-depth understanding of the bonding modes in this compound than that of the original study.
The MESP graph of S 8 , together with all selected local most maxima and minima of potential on the electrostatic surface of the molecule, are shown in Figure 5. The molecule has eight pairs of equivalent maxima of potential (V S,max = 14.0 kcal mol −1 each) that arise along each S-S bond extension, slightly off the σ-axis; each maximum is associated with a σ-hole on S (cf. Figure 5a,b). The centroid region of the molecule features two maxima of potential, one on the top and one at the bottom surface; they are equivalent (V S,max = 5.3 kcal mol −1 each) and weaker than the σ-hole on S (Figure 5b). Each of these electrophilic regions on covalently bound S has the potential for chalcogen bonding interactions when in close proximity of a negative site. By contrast, the lateral portion of each covalently bound S atom in S 8 has a negative potential (Figure 5a). These loci of negative potential toward the interior of the S 8 ring are more negative (V S,min = −5.6 kcal mol −1 ) than those on the exterior surface of the ring (V S,min = −3.4 kcal mol −1 ). potential on the electrostatic surface of the molecule, are shown in Figure 5. The molecule has eight pairs of equivalent maxima of potential (VS,max = 14.0 kcal mol −1 each) that arise along each S-S bond extension, slightly off the σ-axis; each maximum is associated with a σ-hole on S (cf. Figure 5a,b). The centroid region of the molecule features two maxima of potential, one on the top and one at the bottom surface; they are equivalent (VS,max = 5.3 kcal mol −1 each) and weaker than the σ-hole on S (Figure 5b). Each of these electrophilic regions on covalently bound S has the potential for chalcogen bonding interactions when in close proximity of a negative site. By contrast, the lateral portion of each covalently bound S atom in S8 has a negative potential (Figure 5a). These loci of negative potential toward the interior of the S8 ring are more negative (VS,min = -5.6 kcal mol −1 ) than those on the exterior surface of the ring (VS,min = -3.4 kcal mol −1 ). The authors of the original study [128] suggested the presence of intermolecular bonding between I and S sites, basing this conclusion on the intermolecular distance between I and S atoms, and within the conceptual framework of Lewis structures of chemical bonding. In particular, they observed that each iodine atom has two pairs of sulfur neighbors with I···S separation distances lying between 3.78 and 3.88 Å. The stability of the structure was attributed not only to these I···S charge transfer interactions, but also to some extent to Sb···I interactions. However, it was not clear whether these interactions are halogen bonds or chalcogen bonds in the first case, or halogen bonds or pnictogen bonds in the second case, since these bonding terminologies were not coined in 1963.
With the computational tools now available, we are able to more precisely elucidate the nature of these interactions. Our analysis indicates that each I site in SbI3 is involved in three S···I interactions. Two of them are equivalent chalcogen bonds, with r(S···I(Sb)) = 3.883 Å. They are a result of attractive coulombic interaction between the lateral negative portions of bound I atoms in SbI3 and the axial positive portions of the S atom along the S-S bond extensions in S8. They are relatively far from being linear, Type-IIa (∠S-S···I(Sb) = 154.4°). The third S···I interaction is a quasi-linear Type-IIa (Sb)I···S halogen bond, since The authors of the original study [128] suggested the presence of intermolecular bonding between I and S sites, basing this conclusion on the intermolecular distance between I and S atoms, and within the conceptual framework of Lewis structures of chemical bonding. In particular, they observed that each iodine atom has two pairs of sulfur neighbors with I···S separation distances lying between 3.78 and 3.88 Å. The stability of the structure was attributed not only to these I···S charge transfer interactions, but also to some extent to Sb···I interactions. However, it was not clear whether these interactions are halogen bonds or chalcogen bonds in the first case, or halogen bonds or pnictogen bonds in the second case, since these bonding terminologies were not coined in 1963.
With the computational tools now available, we are able to more precisely elucidate the nature of these interactions. Our analysis indicates that each I site in SbI 3 is involved in three S···I interactions. Two of them are equivalent chalcogen bonds, with r(S···I(Sb)) = 3.883 Å. They are a result of attractive coulombic interaction between the lateral negative portions of bound I atoms in SbI 3 and the axial positive portions of the S atom along the S-S bond extensions in S 8 . They are relatively far from being linear, Type-IIa (∠S-S···I(Sb) = 154.4 • ). The third S···I interaction is a quasi-linear Type-IIa (Sb)I···S halogen bond, since the σ-hole on I along the Sb-I bond extension is in an attractive engagement with the negative site on S in S 8 (Figure 5a).
Similarly, each covalently bound Sb site along the I-Sb bond extension in SbI 3 is bound to the nearest negative sites localized on three S atoms of the three nearest S 8 molecules, forming three equivalent Type-IIa I-Sb···I bonding interactions (r(Sb···S) = 3.391 Å and ∠I-Sb···S = 152.5 • for each). These interactions have the characteristics of quasi-linear pnictogen (stibium) bonds.
Other than the three predominant interactions noted above, the crystal features numerous S-S(σ)···S type interactions that are quasi-linear. The deviation from linearity is probably because the maximum of potential on the electrostatic surface of S (σ-hole) is off the S-S bond axis ( Figure 5b). Nevertheless, we observed Type-IIa S-S(σ)···S interactions between the nearest S 8 rings, featuring interaction distances of 3.469, 3.484, and 3.755 Å, corresponding to ∠S-S(σ)···S of 152.5, 161.9, and 169.5 • , respectively. Clearly, the extensive network of Type-II noncovalent interactions, together with numerous Type-I S···S interactions between the interacting units (not shown), determine the shape of the crystal, and can be rationalized from Figure 6a Similarly, each covalently bound Sb site along the I-Sb bond extension in SbI3 is bound to the nearest negative sites localized on three S atoms of the three nearest S8 molecules, forming three equivalent Type-IIa I-Sb···I bonding interactions (r(Sb···S) = 3.391 Å and ∠I-Sb···S = 152.5° for each). These interactions have the characteristics of quasi-linear pnictogen (stibium) bonds.
Other than the three predominant interactions noted above, the crystal features numerous S-S(σ)···S type interactions that are quasi-linear. The deviation from linearity is probably because the maximum of potential on the electrostatic surface of S (σ-hole) is off the S-S bond axis ( Figure 5b). Nevertheless, we observed Type-IIa S-S(σ)···S interactions between the nearest S8 rings, featuring interaction distances of 3.469, 3.484, and 3.755 Å, corresponding to ∠S-S(σ)···S of 152.5, 161.9, and 169.5°, respectively. Clearly, the extensive network of Type-II noncovalent interactions, together with numerous Type-I S···S interactions between the interacting units (not shown), determine the shape of the crystal, and can be rationalized from Figure 6a,b. Our characterization of intermolecular bonding interactions in SbI3:3S8 are supported by the IGM-δg based isosurface topologies shown in Figure 6b, in which, there are clear isosurface volumes between bound atomic basins that are bonded with each other. Another adduct between an antimony halide and molecular sulfur is SbCl 3 .S 8 ; see Figure 7. It crystallizes in the triclinic P-1 space group [129]. The unit-cell contains two SbCl 3 and two S 8 molecules, and the crystal is zero-dimensional [130]. Sb in each SbCl 3 molecule links with several S sites of the crown-shaped S 8 ring, as well as with the negative sites on Cl in close proximity (see dotted lines in Figure 7a). Compared to the formal Sb-Cl bond distances (r(Sb-Cl) = 2.60-2.34 Å; see Figure 7b) between Sb and Cl which are responsible for the trigonal pyramidal geometry of SbCl 3 , the Sb···S and Sb···Cl intermolecular interactions are much longer (r(Sb···S) in the range 3.448-4.320 Å; r(Sb···Cl) = 3.228 Å; Figure 7b). Several of them are non-linear and those appearing along the three Cl-Sb extensions are quasi-linear. All these interactions are Sb-centered pnictogen bonds. Additionally, there are Cl···S halogen bonds, and S···Cl and S···S chalcogen bonds in the crystal (not shown). The formation of these intermolecular interactions is in agreement with the electrostatic surface features of the interacting molecules revealed by the MESP model discussed above (see Figures 2b and 4 for SbCl 3 and S 8 , respectively); together they are responsible for the stability of the SbCl 3 .S 8 system in the crystalline phase. Theoretically, the SbCl 3 .S 8 crystal system is thermally stable with a bandgap of 2.88 eV [130], and is thus a semiconducting material.
ence of a σ-hole, and the intermolecular interactions are represented by dotted lines in cyan or green. Atom type is shown in each case.
Another adduct between an antimony halide and molecular sulfur is SbCl3.S8; see Figure 7. It crystallizes in the triclinic P-1 space group [129]. The unit-cell contains two SbCl3 and two S8 molecules, and the crystal is zero-dimensional [130]. Sb in each SbCl3 molecule links with several S sites of the crown-shaped S8 ring, as well as with the negative sites on Cl in close proximity (see dotted lines in Figure 7a). Compared to the formal Sb-Cl bond distances (r(Sb-Cl) = 2.60-2.34 Å; see Figure 7b) between Sb and Cl which are responsible for the trigonal pyramidal geometry of SbCl3, the Sb···S and Sb···Cl intermolecular interactions are much longer (r(Sb···S) in the range 3.448-4.320 Å; r(Sb···Cl) = 3.228 Å; Figure 7b). Several of them are non-linear and those appearing along the three Cl-Sb extensions are quasi-linear. All these interactions are Sb-centered pnictogen bonds. Additionally, there are Cl···S halogen bonds, and S···Cl and S···S chalcogen bonds in the crystal (not shown). The formation of these intermolecular interactions is in agreement with the electrostatic surface features of the interacting molecules revealed by the MESP model discussed above (see Figure 2b and Figure 4 for SbCl3 and S8, respectively); together they are responsible for the stability of the SbCl3.S8 system in the crystalline phase. Theoretically, the SbCl3.S8 crystal system is thermally stable with a bandgap of 2.88 eV [130], and is thus a semiconducting material.

Antimony Trihalide Crystals
To validate the results of the MESP model summarized in Table 2, we examined the structure of SbF3, reported in 1970 [131]. The unit cell and 2 × 2 extended crystal structure of SbF3 are shown in Figure 8a,b. The dotted lines between the SbF3 units represent intermolecular interactions. Each Sb has three close fluorine neighbors with mean Sb-F distances 1.92 Å and each SbF3 unit is linked with three fluorine bridges (r(Sb···F) = 2.61 Å) to form a pseudo three-dimensional network, with a much-distorted octahedral co-ordination around the antimony atom.

Antimony Trihalide Crystals
To validate the results of the MESP model summarized in Table 2, we examined the structure of SbF 3 , reported in 1970 [131]. The unit cell and 2 × 2 extended crystal structure of SbF 3 are shown in Figure 8a,b. The dotted lines between the SbF 3 units represent intermolecular interactions. Each Sb has three close fluorine neighbors with mean Sb-F distances 1.92 Å and each SbF 3 unit is linked with three fluorine bridges (r(Sb···F) = 2.61 Å) to form a pseudo three-dimensional network, with a much-distorted octahedral co-ordination around the antimony atom.
Our MESP calculations, discussed in Section 8.1. Antimony Trihalides above, indicate that the surface of Sb is strongly positive along the F-Sb bond extensions, while F is entirely negative. This suggests that there are three σ-holes along the extensions of the three F-Sb bonds. The nature of the electrostatic potential surfaces of Sb and F explains why is there a coulombic attraction between Sb in one SbF 3   Other than these interactions, we found a pattern of a parallel arrangement between the SbF 3 units along the crystallographic b-axis, in which the axes of neighboring SbF 3 units are aligned (Figure 8c, top). In this arrangement, Sb in an SbF 3 unit faces the center of the triangular face formed by three F atoms in another a neighboring SbF 3 unit (see the first and third columns of the cluster shown in Figure 8b). These are relatively weak and non-linear Type-Ib pnictogen bound interactions, with r(Sb···F) between 3.60 and 3.80 Å and ∠F-Sb···F in the range 131.8-136.7 • . They are also observable in the IGM-δg based isosurface plot shown in Figure 8d,e.
The extensive network of pnictogen bonding not only drives the formation of the pseudo 3D framework of the crystalline form of SbF 3 , but also the development of secondorder nonlinear optical (NLO) properties [132,133]. In fact, the SbF 3 system exhibits a powder second harmonic generation (SHG) which is about 5.8 times that of potassium dihydrogen phosphate (KDP), and is an NLO material in the infrared region. This allows excellent transparency in the range of 0.29-12 µm with high thermal stability. The band gap of the material calculated from the UV-vis-near-infrared spectrum was about 4.3 eV [132].
The 2 × 2 extended cementite crystal structures of SbX 3 (X = Cl, Br) are shown in Figure 9a-c. The unit-cell for each contains four molecular formula units SbX 3 crystallized in an orthorhombic space group-as SbI 3 , vide supra-SbCl 3 in Pnma [134] and Pbnm [135]; SbBr 3 in Pbnm [136]. There are two equivalent Br sites in the structures, forming a shorter and two longer, equivalent, Sb-X bond distances. The structures are zero-dimensional and Sb 3+ is bound in a distorted T-shaped geometry to three halogen atoms [137,138]. Kang et al. [133] have demonstrated that SbX 3 (X = Cl, Br) are good mid-IR NLO crystals materials. The DFT+U calculated bandgaps were calculated to be 3.751 and 3.119 eV for SbCl 3 [138] and SbBr 3 [137], respectively, suggesting that they are wide band gap semiconductors.
Our geometric analysis of the structure shown in Figure 9b indicates that Sb in SbCl 3 is linked to nine Cl neighbors in five surrounding SbCl 3    There are extensive Type-I and Type-II Sb-X···X and Sb-X···X-Sb (X = Cl, Br) halogen bound contacts in each of the three crystals discussed above. They are shown in Figure 9d, Figure 9e and Figure 9f for SbCl 3 (Pnma), SbCl 3 (Pbnm) and SbBr 3 (Pbnm), respectively. The Type-I Cl···Cl contacts in SbCl 3 that occur between the zig-zag sheets are shorter than those within the sheet itself, and hence responsible for the sheet's development (viz. 3.579 Å versus 3.630 Å in Figure 9d By contrast, the SbI 3 system crystallizes in two space groups, i.e., monoclinic (P2 1 /c) and trigonal (R-3), and is non-magnetic [139]. The crystal structure of the latter, depicted in Figure 10, shows that it adopts a layered geometry, and that the neighboring layers are linked to each other through a number of Type-IIa halogen-centered I···I interactions (r(I···I) = 3.958; ∠Sb-I···I = 164.5 • ) that are less than twice the vdW radius of I, 4.08 Å. However, in the monoclinic crystal, the isolated SbI 3 molecules are linked with each other by means of a 3D network of I-Sb···I intermolecular contacts. Both the monoclinic and trigonal structures of SbI 3 possess wide band gap transition energies [140] (2.217 eV for the monoclinic geometry and 2.211 eV for the trigonal geometry), indicative that they are semiconducting materials. Kang et al. have theoretically determined that SbI 3 has a strong second harmonic generation effect [133]. Onodera and coworkers have experimentally demonstrated that SbI 3 is a promising material for use in radiation detectors; it is stable under operating conditions, similar to popular semiconductor detectors CdTe and TlBr [140]. The IGM-δg results given in Figure 11 confirm the intermolecular pnictogen bonding interactions that were identified between the interacting molecules in SbX3 (X = Cl, Br, I) crystals. As discussed above, we observed that there were circular or dumbbell-shaped green isosurfaces between the Sb and X atomic basins in SbX3 (X = Cl, Br). In addition, we also identified a number of potential Type-I and Type-II Sb-X···X halogen bound interactions between the SbX3 molecules in the crystal, as marked in Figure 11a,b for SbCl3 and The IGM-δg results given in Figure 11 confirm the intermolecular pnictogen bonding interactions that were identified between the interacting molecules in SbX 3 (X = Cl, Br, I) crystals. As discussed above, we observed that there were circular or dumbbell-shaped green isosurfaces between the Sb and X atomic basins in SbX 3 (X = Cl, Br). In addition, we also identified a number of potential Type-I and Type-II Sb-X···X halogen bound interactions between the SbX 3 molecules in the crystal, as marked in Figure 11a,b for SbCl 3 and SbBr 3 , respectively. In both cases, the pnictogen bonding was stronger than the Type-I and Type-II halogen bonds, as inferred from the size and thickness of the IGM-δg based isosurfaces. The presence of Type-I and -II halogen bonding interactions between the interacting molecules was not obvious in the unit cell of the SbI 3 crystal. Their presence was rationalized when the unit-cell was expanded. The IGM-δg based isosurface analysis performed on a small portion of the expanded structure of the SbI 3 crystal is shown in Figure 11c (right), including both Type-I and Type-II I···I halogen bonds in the SbI 3 crystal.  There are many H···H and H···Br hydrogen bonding interactions in the crystal, with intermolecular bond distances around 2.374 and 2.990 Å, respectively. These are attractive interactions, and each is less than the sum of the van der Waals radii of the respective atomic basins. For instance, the H···H bond distance between the methyl groups is 2.374 Å, which is less than the sum of the vdW radii of two hydrogen atoms, 2.40 Å (rvdW (H) = 1.20 Å). Similarly, the H···Br distance of 2.990 Å is less than sum of the vdW radii of H and Br atoms, 3.06 Å (rvdW (H) = 1.20 Å and rvdW (Br) = 1.86 Å).
Bis(dimethylstibanyl)oxane, ((CH3)2Sb)2O, is an example of an organo-element species that also features pnictogen bonding [142]. As shown in Figure 13a  There are many H···H and H···Br hydrogen bonding interactions in the crystal, with intermolecular bond distances around 2.374 and 2.990 Å, respectively. These are attractive interactions, and each is less than the sum of the van der Waals radii of the respective atomic basins. For instance, the H···H bond distance between the methyl groups is 2.374 Å, which is less than the sum of the vdW radii of two hydrogen atoms, 2.40 Å (r vdW (H) = 1.20 Å). Similarly, the H···Br distance of 2.990 Å is less than sum of the vdW radii of H and Br atoms, 3.06 Å (r vdW (H) = 1.20 Å and r vdW (Br) = 1.86 Å).
Bis(dimethylstibanyl)oxane, ((CH 3 ) 2 Sb) 2 O, is an example of an organo-element species that also features pnictogen bonding [142]. As shown in Figure 13a These interactions suggest that Sb in crystalline ((CH 3 ) 2 Sb) 2 O is either tetragonal or pentagonal. The authors of [142] suggested that ((CH 3 ) 2 Sb) 2 O adopts a syn-anti conformation in the solid state with the entities arranged in zigzag chains along because of weak intermolecular Sb···O interactions. We suggest that that these interactions are actually reasonably strong, featuring a significant degree of covalency and may be best described as coordinate bonds. A similar proposition, i.e., that halogen bonds have a significant degree of covalency and are best described as analogous to coordinate bonds rather than analogous to hydrogen bonds, has been advanced [143]. We base our suggestion on the observation that there is a significant accumulation of charge density in the bonding region between the Sb and O atomic basins of a pair of interacting molecules (Figure 13c, right), indicated by blue isosurfaces. This was not so for the Sb···Sb interactions molecules (Figure 13c, left), as the isosurfaces corresponding to these were green, a signature of weakly bound interactions. These interactions suggest that Sb in crystalline ((CH3)2Sb)2O is either tetragonal or pentagonal. The authors of [142] suggested that ((CH3)2Sb)2O adopts a syn-anti conformation in the solid state with the entities arranged in zigzag chains along because of weak intermolecular Sb···O interactions. We suggest that that these interactions are actually reasonably strong, featuring a significant degree of covalency and may be best described as coordinate bonds. A similar proposition, i.e., that halogen bonds have a significant degree of covalency and are best described as analogous to coordinate bonds rather than analogous to hydrogen bonds, has been advanced [143]. We base our suggestion on the observation that there is a significant accumulation of charge density in the bonding region between the Sb and O atomic basins of a pair of interacting molecules (Figure 13c, right), indicated by blue isosurfaces. This was not so for the Sb···Sb interactions molecules ( Figure  13c, left), as the isosurfaces corresponding to these were green, a signature of weakly bound interactions.
Since the Sb-O bond distance in ((CH3)2Sb)2O is 2.009 Å, the mutual penetration between the vdW spheres of Sb and O in within a ((CH3)2Sb)2O entity should be larger than the penetration of these domains with those on neighboring entities. The Sb-O bond is, therefore, a covalent-like bond with some significant ionic characteristics (roughly 38%), which is not unexpected for coordinate bonds [144,145]. Because of its high electronegativity, O in ((CH3)2Sb)2O pulls significant electron density towards the O-Sb bonding region, leaving behind a significantly strong positive region (σ-hole) on the surface of the Sb atom along the outermost extension of the O-Sb bond. This σ-hole is responsible for formation of a strong O-Sb···O-Sb pnictogen bonding interaction in the crystal; see Figure Since the Sb-O bond distance in ((CH 3 ) 2 Sb) 2 O is 2.009 Å, the mutual penetration between the vdW spheres of Sb and O in within a ((CH 3 ) 2 Sb) 2 O entity should be larger than the penetration of these domains with those on neighboring entities. The Sb-O bond is, therefore, a covalent-like bond with some significant ionic characteristics (roughly 38%), which is not unexpected for coordinate bonds [144,145]. Because of its high electronegativity,  Figure 13c.
The crystal structure of the sulfane analog ((CH 3 ) 2 Sb) 2 S is also reported (Figure 14a) [142]. The molecules are in an approximately syn-syn conformation in the crystal structure. The Sb-S bond length is 2.449 Å, which is significantly smaller than the Sb···S intermolecular distances of 3.164 Å (Figure 14b), due to a pair of interactions that developed when three units of the ((CH 3 ) 2 Sb) 2 S came in close proximity in the crystal lattice. Both the Sb···S intermolecular contacts appear along the S-Sb bond extensions with ∠S-Sb···S is = 176.8 • , revealing the presence of antimony-centered Type-IIa σ-hole interactions. These directional interactions in the crystal are evident in the IGM-δg based bluish-green isosurfaces shown in Figure 14c (right). In addition, there are Sb···Sb long-ranged stibium bonds in the crystal, with r(Sb···Sb) = 4.208 Å and ∠S-Sb···Sb =134.8 • . These are Type-Ib pnictogen bonds, and show up in the IGM-δg based green isosurfaces shown in Figure 14c (left). There are many H···S hydrogen bonds in the extended crystal lattice (not shown).
The crystal structure of the sulfane analog ((CH3)2Sb)2S is also reported (Figure 14a) [142]. The molecules are in an approximately syn-syn conformation in the crystal structure. The Sb-S bond length is 2.449 Å, which is significantly smaller than the Sb···S intermolecular distances of 3.164 Å (Figure 14b), due to a pair of interactions that developed when three units of the ((CH3)2Sb)2S came in close proximity in the crystal lattice. Both the Sb···S intermolecular contacts appear along the S-Sb bond extensions with ∠S-Sb···S is = 176.8°, revealing the presence of antimony-centered Type-IIa σ-hole interactions. These directional interactions in the crystal are evident in the IGM-δg based bluish-green isosurfaces shown in Figure 14c (right). In addition, there are Sb···Sb long-ranged stibium bonds in the crystal, with r(Sb···Sb) = 4.208 Å and ∠S-Sb···Sb =134.8°. These are Type-Ib pnictogen bonds, and show up in the IGM-δg based green isosurfaces shown in Figure 14c (left). There are many H···S hydrogen bonds in the extended crystal lattice (not shown).
The space-filling and ball-and-stick models of the compound are shown in Figure 15a,b. The first shows that there is an overlapping of Sb and O sites in the crystal, causing the formation of three equivalent Cl-Sb···O interactions (Figure 15c). The Sb···O distance, r(Sb···O), is 2.755 Å and the angle of interaction, ∠Cl-Sb···O, is 177.2 • . These Type-IIa stibium bonds are markedly longer than the three formal Sb-Cl coordinate bonds (2.379 Å) of SbCl 3 . The Sb···O interactions can be viewed as the result of attraction between the regions of positive potential along the Cl-Sb bond extensions and the negative potential on the surface of O atoms in 1,4-dioxane. One may also expect that they have some characters expected of coordinate bonds, and are largely coulombic interactions.
interacting with the lone-pair electrons on three O-sites in three surrounding 1,4-dioxane moieties [146,147].
The space-filling and ball-and-stick models of the compound are shown in Figure  15a,b. The first shows that there is an overlapping of Sb and O sites in the crystal, causing the formation of three equivalent Cl-Sb···O interactions (Figure 15c). The Sb···O distance, r(Sb···O), is 2.755 Å and the angle of interaction, ∠Cl-Sb···O, is 177.2°. These Type-IIa stibium bonds are markedly longer than the three formal Sb-Cl coordinate bonds (2.379 Å) of SbCl3. The Sb···O interactions can be viewed as the result of attraction between the regions of positive potential along the Cl-Sb bond extensions and the negative potential on the surface of O atoms in 1,4-dioxane. One may also expect that they have some characters expected of coordinate bonds, and are largely coulombic interactions. In addition, there are also Cl···Cl halogen bonds and H···O and H···Cl hydrogen bonds between the interacting molecules in the crystal. The first developed between the Cl sites In addition, there are also Cl···Cl halogen bonds and H···O and H···Cl hydrogen bonds between the interacting molecules in the crystal. The first developed between the Cl sites on the Sb-Cl bond extensions, with r(Cl···Cl) bond distances of 3.912 Å and angle of interaction, ∠Sb-Cl···Cl, of 172.85 • . Although these interactions are directional and have a Type-IIa topology of bonding, they are not actually Type-IIa interactions (∠Sb-Cl···Cl = 172.85 • ; and ∠Cl-Cl···Cl = 172.85 • ), since the Sb-Cl···Cl-Sb interactions are formed between regions of nearly identical electrostatic potentials centered on the interacting Cl atoms. They may therefore be best characterized as Type-III halogen bound interactions [22]. On the other hand, the C-H···Cl hydrogen bonds are formed between the negative sites on the Cl atoms bonds around the Sb-Cl bonds and the nearest H sites on 1,4-dioxane, with r(H···Cl) varies in the range 2.9-3.3 Å. Each Cl site in SbCl 3 serves as an acceptor of six H···Cl hydrogen bonds, with bond distances of 2.928, 3.023, 3.106, 3.145, 3.148 and 3.198 Å. The C-H···O hydrogen bonds in the crystal are formed between the 1,4-dioxane moieties, the majority of them with a bond distance of 3.315 Å and an interaction angle of 161.1 • . There are four such bonds between a pair of two 1,4-dioxane molecules, locally forming a sandwich type 1,4-dioxane dimer (Figure 15d). Most of the hydrogen bonds are weak, and they originate from attraction between the positive site on H and the negative site(s) on Cl or O, making a local contribution to the overall stability of the crystal.

The Crystal Structure of [Co(trien)(NSC) 2 ][Sb 2 (tart)(Htart)]
Kushi and co-workers [148] reported the structure of the diastereomeric salt formed between trans-[Co(trien)(NSC) 2 ] + (trien = triethylenetetramine) and [Sb 2 (tart)(Htart)] -(tart = tartrate, C 4 H 2 O 6 2-). It features each Sb 3+ ion coordinated to four O atoms of tartrate in a distorted seesaw coordinated environment. Sb is engaged in an attractive engagement with a π-cloud around the C=N bond of the coordinated isothiocyanato anion (Figure 16a,b). The N=C and C=S bonds in the metal coordinate SCNanion form two interactions with an Sb 3+ ion of the cation, Sb(σ)···π(SCN). These are probably weak, charge assisted interactions, corresponding to r(Sb···π(midpoint of N=C)) = 3.437 Å and r(Sb···π(midpoint of C=S) = 3.607 Å). The sum of the vdW radii of Sb and S atomic basins is 4.36 Å, considerably longer than the distance between Sb and S (3.929 Å; see Figure 16b). The interaction is also directional, yet Type-IIa as ∠O-Sb···S = 158.5 • . This contact between Sb and S in the crystal is not a pnictogen bond, but rather, the result of charge assisted π-hole centered tetrel bonding, since trans-[Co(trien)(NSC) 2 ] + carries a formal charge of +1. Similarly, Sb···π(midpoint of N=C) is a π-hole centered tetrel bond, yet Sb is closer to C than to N. These two interactions, together with the four Sb-O bonds to tartrate, make the coordination environment of Sb 3+ pseudo-octahedral in this salt (Figure 16b). Apparently, the packing in the crystal is not only driven by chalcogen and tetrel bonding interactions, but also by several other non-covalent interactions formed between the interacting species (not shown). This example implies that interactions that appear to be the result of pnictogen bonding may actually correspond to something else if assigned correctly.

The Crystal of [SbCl2imR2R'2][OTf]
Henne et al. reported a series of pnictogen complexes featuring a Pn-C bond in the C2 position of a sterically hindered imidazole (imR2R'2 = 1,3-dipropyl-2,5-dimethyl imidazole) [149]. One of these, [SbCl2imR2R '2][OTf] (OTf = triflate), is shown in Figure 17a. As discussed above for other systems, the solid-state structure arises from several chargeassisted intermolecular interactions, as Sb is entirely positive. The three Sb···O contacts shown in Figure 17b are a result of coulombic attraction between it and the entirely negative O sites on two neighboring OTf anions. Two of these interactions are significantly

The crystal structures of phosphoryl isothiocyanate
Phosphoryl isothiocyanate, OP(NCS)3 can form many types of crystal structures [150]. An example is [SbCl3·OP(NCS)3]4, which crystallizes in the cubic space group I4 3m, shown in Figure 18a [150]. The three chloride ligands of SbCl3 are arranged in a primarily trigonal pyramidal geometry around Sb 3+ , which also interacts with three O sites of three neighboring OP(NCS)3 units. We also identified a long Sb···Cl bond between neighboring cations with r(Sb···Cl) = 4.123 Å and an ∠Cl-Sb···Cl = 129.4 • (Type-IIb). This is also a weak, non-linear pnictogen bond, given that it is marginally shorter than the sum of the vdW radii of the Sb and Cl atomic basins, 4.29 Å. The Sb center in [SbCl 2 imR 2 R' 2 ] + is therefore pseudo seven coordinate.

The Crystal Structures of Phosphoryl Isothiocyanate
Phosphoryl isothiocyanate, OP(NCS) 3 can form many types of crystal structures [150]. An example is [SbCl 3 ·OP(NCS) 3 ] 4 , which crystallizes in the cubic space group I43m, shown in Figure 18a [150]. The three chloride ligands of SbCl 3 are arranged in a primarily trigonal pyramidal geometry around Sb 3+ , which also interacts with three O sites of three neighboring OP(NCS) 3 units.  As shown in Figure 18b, there are three links between the Sb center in SbCl 3 and the three O sites of OP(NCS) 3 moieties. These are a result of coulombic attractions between the positive potential on Sb and the negative potential associated with the lone-pair electrons on O. They are all equivalent and are much longer (3.060 Å) than the formal Sb-Cl bonds (2.341 Å). Because the Sb center interacts with six atomic basins, this creates a pseudooctahedral around it. Therefore, the complex system exhibits a tetranuclear "cage" structure with an [Sb 4 O 4 ] core, in which four corners of the core form a regular cube of edge length 3.060 Å that are occupied by four Sb 3+ ions and the remaining four are occupied by O sites. The [Sb 4 O 4 ] core consists of µ 3 -bridging oxygen and six coordinated Sb 3+ ion such that each O in OP(NCS) 3 acts as a trifurcated center to accept electrophilic attacks. All four lone pairs of the four Sb atoms suggested in ref. [150] are located in the tetranuclear closed "cage" structure of the [Sb 4 O 4 ] core. Because the Sb···O intermolecular interactions are substantially longer than the Sb-Cl bonds, and appear along the Cl-Sb extensions (∠Cl-Sb···O = 179.3 • ), we characterize them as Type-IIa stibium bonds.
Our further analysis suggested that Cl, along the Sb-Cl bond extension, also engages attractively with the S atom in OP(NCS) 3 . Each SbCl 3 unit in the crystal donates three chlorine-centered σ-holes to three S, forming three equivalent Cl···S halogen bonds (r(Cl···S) = 3.613 Å and ∠Sb-Cl···S = 166.7 • ). These halogen-centered σ-hole interactions leads to the emergence of macrocyclic cage-like pseudo 3D structures as shown in Figure 18c.
The antimony atoms form a distorted tetrahedron with an Sb···Sb edge of 4.186(1) Å, which is much less than the sum of the vdW radii of two Sb atoms (4.94 Å). Because of this, a weak specific interaction between Sb ions was postulated [150]. Other than this, we also expect that the crystal should contain weak (CN)π···Cl interactions. These postulated intermolecular interactions that are responsible for the stability of the [SbCl 3 ·OP(NCS) 3 ] 4 crystal are elucidated by our IGM-δg analysis. The results are summarized in Figure 18d for a small portion of the crystal, consisting of three SbCl 3 units and three OP(NCS) 3 units. The IGM-δg results show that the Sb···Sb pnictogen bonds are very weak, and appear as a thin flat isosurface colored green. The Cl···S interaction is a localized interaction and the isosurface volume between Cl and S is tiny. By contrast, the isosurface volumes between the Sb and O atomic basins are thick and circular, and bluish-green, which indicate strongly bound stibium bonds.
Similar pnictogen bonding topologies were found in a crystal of 1,3,5-Trithiane antimony tribromide, SbX 3 ·(CH 2 S) 3 (X = Cl, Br) [151]. The packing between SbX 3 and (CH 2 S) 3 is illustrated in Figure 19a,b for SbCl 3 ·(CH 2 S) 3 and SbBr 3 ·(CH 2 S) 3 , respectively. The H atom positions are not provided in the CSD deposition, and thus, are missing in Figure 19. The intermolecular interactions between the interacting units in the crystal are shown in Figure 19c,d, respectively. Because of the absence of H atom positions in CH 2 S in the crystal, the probable role of H cannot be determined. H is expected to play a role in the packing by the formation of several non-covalent interactions. One such prominent interaction in SbBr 3 ·(CH 2 S) 3 could be the (HC)H···Br/(C)H···Br hydrogen bonds, for example.
As found in other structures, see above, the Sb center in SbX 3 ·(CH 2 S) 3 (X = Cl, Br) is pseudo-octahedral. In SbCl 3 ·(CH 2 S) 3 , the Sb-Cl distances are 2.374 Å, and the Sb···S distances are much longer, i.e., 3.257 Å. Similarly, in SbBr 3 ·(CH 2 S) 3 , the Sb-Br distances are 2.534 Å, while the Sb···S distances are, 3.365 Å. The presence of the exceptionally Sb···S long bonds signifies that they are not coordinate bonds, but Type-IIa stibium bonds, and that the SbX 3 units are locally trigonal pyramidal. We note further that there are many S···Cl and S···Br chalcogen bonds in the crystal, and each S site is involved in Chen and co-workers reported an organic-inorganic hybrid birefringent material, i.e., (CN4H7)SbC2O4F2(H2O)0.5, comprising π-conjugated organic groups [CN4H7] + and lin-   (Figure 20a) [152]. The compound exhibits a large birefringence (∆n = 0.126 @ 546 nm) that is almost equal to that of the well-known birefringent material, α-BaB 2 O 4 . Sb is bound to four O and two F donors, forming a distorted pentagonal pyramidal structure. There is therefore significant space around the coordination sphere of Sb, opposite the F-Sb bonds, for which the Sb center is able to accommodate interactions with negative sites, just as we observed. The positive site on Sb is involved in F-Sb···O (oxalate) and F-Sb···F pnictogen bonding interactions (Figure 20b (Figure 20a) [152]. The compound exhibits a large birefringence (Δn = 0.126 @ 546 nm) that is almost equal to that of the well-known birefringent material, α-BaB2O4. Sb is bound to four O and two F donors, forming a distorted pentagonal pyramidal structure. There is therefore significant space around the coordination sphere of Sb, opposite the F-Sb bonds, for which the Sb center is able to accommodate interactions with negative sites, just as we observed. The positive site on Sb is involved in F-Sb···O (oxalate) and F-Sb···F pnictogen bonding interactions (Figure 20b

Stibium Bonds Formed with Arene Moieties and Other Cyclic and Non-cyclic systems
Covalently bound Sb can form pnictogen bonds with π-density rich regions. For example, Bombieri and co-workers [153] reported such a case, i.e., the 2:1 complex between SbBr3 and pyrene. Figure 21a figure illustrates the prominence of pnictogen bonds between interacting SbBr3 units, as well as between SbBr3 and pyrene. Sb 3+ in each SbBr3 unit simultaneously forms three Br-Sb···Br pnictogen bonds with the Br sites in neighboring SbBr3 units, which act as pnictogen bond acceptors. At the same time, the Sb forms a pnictogen bond with the π-density on the aromatic rings of pyrene which also act as pnictogen bond acceptors. Most of the Br-Sb···Br and Br-Sb···π(C) interactions are directional, with

Stibium Bonds Formed with Arene Moieties and Other Cyclic and Non-Cyclic Systems
Covalently bound Sb can form pnictogen bonds with π-density rich regions. For example, Bombieri and co-workers [153] reported such a case, i.e., the 2:1 complex between SbBr 3 and pyrene (Figure 21a). The figure illustrates the prominence of pnictogen bonds between interacting SbBr 3 units, as well as between SbBr 3 and pyrene. Sb 3+ in each SbBr 3 unit simultaneously forms three Br-Sb···Br pnictogen bonds with the Br sites in neighboring SbBr 3 units, which act as pnictogen bond acceptors. At the same time, the Sb forms a pnictogen bond with the π-density on the aromatic rings of pyrene which also act as pnictogen bond acceptors. Most of the Br-Sb···Br and Br-Sb···π(C) interactions are directional, with different ∠Br-Sb···Br values of between 170 • and 175 • (Figure 21c), and ∠Br-Sb···π(C 6 ) = 157.8 • (Figure 21d).  The lateral negative sides around Br in SbBr 3 are also involved in an extensive number of C-H···Br hydrogen bonds. There are many such intermolecular non-covalent interactions between SbBr 3 and pyrene, with Br···H distances varying between 2.85 and 3.20 Å. They display a Type-IIa directional feature, with the angle of interaction, ∠C-H···Br in the 150.0-170.0 • range, as shown in Figure 21b. There are also Type-I hydrogen bonds; they are longer, with r(H···Br)~3.0-3.3 Å, and ∠C-H···Br in the 90.0-150 • range.
The crystal also features a significant number of Sb-Br···Br(Sb) contacts between the SbBr 3 units, between 3.60 and 3.90 Å in length. The bond distances for Type-I Sb-Br···Br(Sb) contacts occurring between the SbBr 3 units are close to 3.89 Å, with bond angles lying between 90 • and 120 • (Figure 21c). The ability of Br in SbBr 3 to engage in non-covalent interactions is remarkable. It is involved either in (i) two Type-I and one Type-II Br···Br halogen bound interactions (∠Sb-Br···Br = 160.3 • ); or (ii) in three Type-I and one Type-II Br···Br contacts; or (iii) in one hydrogen bond, one pnictogen bond, one Type-II and two Type-I Br···Br contacts; (iv) or in a set of one Sb-Br···Br pnictogen bond and two Br···H hydrogen bond(s); (v) or in a set of two halogen bonds, one pnictogen bond, and one hydrogen bond; (vi) or in three halogen bonds and one hydrogen bond; (vi) or in a set of four halogen bonds and one hydrogen bond; (vii) or in a set of one pnictogen bond, two hydrogen bonds, and one Br···π(pyrene ring) interactions.
In addition, there are also slipped-parallel π···π staking interactions between the aromatic rings of pyrene. This complicated bonding topology emerges from the variable orientation of the SbBr 3 molecule in the crystal, which plays a crucial role in packing with the pyrene entity. While all types bonding interactions responsible for the crystal are not explicitly discussed, in Figure 21e,f, we present the IGM-δg based isosurface results calculated between three SbBr 3 units, and those between SbBr 3 and pyrene, the geometries of which were extracted from the crystal. As discussed above, our IGM-δg analysis provides evidence of Type-I and Type-II halogen bonding between the SbBr 3 units-described by circular thin green volumes ( Figure 21e); Sb···Br pnictogen bonds-described by bluishgreen circular or dumbbell-shaped volumes (Figure 21e); and Sb···π pnictogen bondsdescribed by a cone-like volume (Figure 21f).  [118]; and tetrakis(trichloro-antimony) 1,2-diphenylethyne ( Figure 1l) [119]. In all cases, the intermolecular distances are less than 3.5 Å, and not unexpected for π-centered intermolecular interactions.
In all these phosphine chalcogenide complexes, the Sb-X bond distances are significantly distorted from what is expected of a trigonal pyramidal SbX 3 molecular entity. The three Sb-X bonds of SbX 3 in its complex with Cy 3 PS are not equivalent, and each of them is significantly shorter than the Sb-S bond. The Sb···X bonds responsible for the (Cy 3 PS)SbX 3 (X = Cl, Br, I) dimers are significantly longer than ordinary coordinate bonds, and are directional. They are inequivalent in (Cy 3   The structure of [S(CH2-2-C6H4SbMe3)2]I2, resulting from pnictogen bonding and extensive other primary and secondary interactions, is shown in Figures 22e and 23 [155]. A representation of the ionic species, [S(CH2-2-C6H4SbMe3)2] 2+ ·2I − , is shown in Figure 23a, and the expanded crystal system is shown in Figure 23b. Since the pnictogen bonding is between (formally) ionic species, it is charge assisted.  Figures 22e and 23 [155]. A representation of the ionic species, [S(CH 2 -2-C 6 H 4 SbMe 3 ) 2 ] 2+ ·2I − , is shown in Figure 23a, and the expanded crystal system is shown in Figure 23b. Since the pnictogen bonding is between (formally) ionic species, it is charge assisted.
While appreciating that [S(CH 2 -2-C 6 H 4 SbMe 3 ) 2 ]I 2 is an ionic species, for convenience, we shall refer to it as a "molecular unit". The detailed topology of antimony-centered pnictogen bonding cannot be inferred from the molecular unit shown in Figure 23a, and is best appreciated from Figure 22e. There are two different C-Sb···I links and two different C-Sb···S links, two intermolecular, and two intramolecular. Except for the C-Sb···S bond with an intramolecular distance of 3.555 Å, the other three are significantly longer. The C-Sb···I links have the characteristics of Type-II pnictogen bonds, whereas one of the two C-Sb···S bonds, at an intramolecular distance of 4.419 Å, is a characteristic of Type-I.
Each of the contacts is shorter than the sum of the vdW radii of the respective atomic basins, 4.51 Å. Similarly, the intramolecular distance of one of the two C-Sb···S bonds (r(Sb···S) = 3.555 Å) is shorter than the sum of the vdW radii of S and Sb, 4.36 Å, while the other is slightly longer. These observations indicate the formation of two intramolecular hypervalent C-Sb···S pnictogen bonds. The nature of the intermolecular interactions is revealed by inspecting the molecular ordering in the expanded crystal, a part of which is shown in Figure 22e. While appreciating that [S(CH2-2-C6H4SbMe3)2]I2 is an ionic species, for convenience, we shall refer to it as a "molecular unit". The detailed topology of antimony-centered pnictogen bonding cannot be inferred from the molecular unit shown in Figure 23a, and is best appreciated from Figure 22e. There are two different C-Sb···I links and two different C-Sb···S links, two intermolecular, and two intramolecular. Except for the C-Sb···S bond with an intramolecular distance of 3.555 Å, the other three are significantly longer. The C-Sb···I links have the characteristics of Type-II pnictogen bonds, whereas one of the two C-Sb···S bonds, at an intramolecular distance of 4.419 Å, is a characteristic of Type-I. Each of the contacts is shorter than the sum of the vdW radii of the respective atomic basins, 4.51 Å. Similarly, the intramolecular distance of one of the two C-Sb···S bonds (r(Sb···S) = 3.555 Å) is shorter than the sum of the vdW radii of S and Sb, 4.36 Å, while the other is slightly longer. These observations indicate the formation of two intramolecular Each of the two Sb centers of the cation is formally four-coordinate. Since each Sb center has four positive σ-holes along the four C-Sb bond extensions, they are in coulombic engagement with four nearest iodide anions. The resulting interactions, i.e., C-Sb···I, are all long-ranged and directional. We also note that the (Me)H···I hydrogen bonds are key driving forces causing the occurrence of these charge-assisted pnictogen bonds in the crystal system (r(H···I) = 3.0-3.5 Å). To confirm these observations, we performed an IGM-δg based isosurface analysis, using the geometry of the ion-pair shown in Figure 23a. The results are summarized in Figure 23c,d. The first suggests that the IGM-δg isosurface volumes in green and bluish-green representing the long and short interactions between Sb and S within the cation are weak and medium strength interactions, respectively. Similarly, the broad isosurface volumes in green appearing between I and H atomic basins in Figure 23d suggest that the (Me)H···I hydrogen bonds are the primary interactions that lead to the formation of C-Sb···I pnictogen bonds in the crystal.
Shown in Figure 25 are the complexes of SbX 3 (X = Cl, Br) with dibenzo-24-crown-8, 18-crown-6, and 12-crown-4. The large cavity of dibenzo-24-crown-8 shows the ability to accommodate two SbX 3 units (Figure 25a,d). The O lone-pairs of the crown ether face the electrophilic regions on SbBr 3 (see Figure 25d), leading to the formation of four Sb···O contacts of between 2.893(7) and 3.183(8) Å), while the hydrophobic methylene bridges face the three Br sites in the neighboring SbBr 3 unit with Sb-Br bond lengths of 2.5268(15), 2.5478 (15), and 2.5558(14) Å). The Sb···O contacts may not represent coordinate bonds, since the length the contacts (for example, the 3.147 Å contact) would suggest that these are antimony-centered pnictogen bonds. Two of them are quasi-linear, appearing along the Br-Sb bond extensions, while the other two are bent. The strength of these interactions may fall into the border region between coordinate and very strong pnictogen bonds.
The bonding pattern in the dibenzo-24-crown-8 complex with SbCl 3 is somewhat different (Figure 25a) relative to that found with SbBr 3 (Figure 25d). There are now five Sb···O close contacts; no doubt, this is a consequence of the smaller size of SbCl 3 compared to SbBr 3 . The three formal Sb-Cl bonds of SbCl 3 have bond lengths of 2.3807(8), 2.3874(8), and 2.4015(9) Å) that are predominantly shorter. Because the electrostatic surface along the X-Sb bond extensions is positive (cf. Figure 2b,c), it links to the O-sites from opposite sides of the cavity of the crown-ether. This results in a sigmoidal conformation, in which two pseudo O 5 donor cavities are created to satisfy the host-guest matching requirements. With this topology of bonding, the two O sites of the dibenzo-24-crown-8 serve as bifurcated pnictogen bond acceptors with Sb···O contacts in the range 2.852(2)-3.153(2) Å. Regardless of the nature of guest, i.e., SbCl 3 or SbBr 3 , most of the Sb···O pnictogen bonds between SbX 3 (X = Cl, Br) with dibenzo-24-crown-8 were found to be non-linear Type-II.
In order to provide further insight into the coordination chemistry of Sb 3+ ions, we performed IGM-δg analyses on their crystal geometries; our results are shown in Figure 26. The isosurface volumes in the plots shown in the top graphs were obtained between the Sb and X atoms in SbX 3 . They are colored deep blue, indicative of a reasonably high concentration of charge density in the bonding regions between Sb and X, as expected for coordinate bonds in SbX 3 . The remaining IGM-δg based isosurface plots of Figure 26 correspond to different isovalues, illustrating that intermolecular interactions are not faint. That is, the isosurfaces for the very strong interactions between Sb and O always appear, regardless of the isovalue used. They are bluish-green in nature, and demonstrate the stability anticipated from the intermolecular bond distances. However, the isosurfaces corresponding to the longer bonds formed between Sb in one unit of SbCl 3 and Cl in the second unit of the same molecule do not show up for large isovalues (0.022, 0.020 and 0.14 a.u.). This is not very surprising given that small isovalues are necessary to visualize the weak and van der Waals interactions that correspond to regions with a very low charge density gradient.
tacts of between 2.893(7) and 3.183(8) Å), while the hydrophobic methylene bridges face the three Br sites in the neighboring SbBr3 unit with Sb-Br bond lengths of 2.5268(15), 2.5478 (15), and 2.5558(14) Å). The Sb···O contacts may not represent coordinate bonds, since the length the contacts (for example, the 3.147 Å contact) would suggest that these are antimony-centered pnictogen bonds. Two of them are quasi-linear, appearing along the Br-Sb bond extensions, while the other two are bent. The strength of these interactions may fall into the border region between coordinate and very strong pnictogen bonds.  When an isovalue of 0.008 a.u. was used, most of the chemical interactions between SbCl 3 and dibenzo-24-crown-8 were revealed. In addition to several H···Cl interactions, we observed two additional Sb···Cl interactions between the two SbCl 3 units marked by circles in Figure 26a (bottom plots). They are longer (r(Sb···Cl) = 3.957 Å each) and bent (∠Cl-Sb···Sb(Cl) = 76.4 • ). These results suggest that antimony in (dibenzo-24-crown-8)(SbCl 3 ) 2 is eight-fold coordinated. At least three of the contacts (one Sb···Cl and two Sb···O crown ) are of weak-to-medium strength and are stibium bonds.
A similar result was obtained with the IGM-δg based isosurface analysis for (dibenzo-24-crown-8)(SbBr 3 ) 2 , but there are no (Br-)Sb···Br links between the two SbBr 3 units. However, the dependence of the size of the isosurface on the isovalue was notable (see Figure 26b, bottom ones); there were no IGM-δg isosurfaces found between Sb and O for an isovalue of 0.022 a.u., and they were revealed only when an isovalue < 0.020 a.u. was used, suggesting that they are genuine medium-strength pnictogen bonds. Figure 25b,c display the local bonding environment between the host and guest species in (12-crown-4)SbCl 3 [162] and (18-crown-6)SbCl 3 [161], respectively. Because the cavity of 12-crown-4 is smaller than that of 18-crown-6, the Sb 3+ in SbCl 3 forms four Sb···O contacts in its complex with the first, but six with the second, resulting in a bowl-shaped architecture with a cone-like facial stand formed by SbCl 3 .
As a means of verifying these tentative deductions concerning the nature of the bonding in these systems, the IGM-δg isosurfaces were calculated and are shown in Figure 27. There are four thick, bluish circular volumes between Sb and O atoms in Figure 27a for (12-crown-4)SbCl 3 that signify Sb-O covalent bonds [162]. In addition, we note interactions between neighboring (12-crown-4)SbCl 3 units: a Cl···Cl halogen bond; H···Cl hydrogen bonds between the methylene hydrogens of the crown ether and the negative sites on Cl of SbCl 3 ; and H···H dihydrogen type attractive interactions between the methylene groups of the crown ether (Figure 27a). In the case of (18-crown-6)SbCl 3 , there are three greenish and three bluish-green thick circular volumes between Sb and O (Figure 27b). The first three are pnictogen bonds and the others are genuine coordinate bonds. In addition, there are several H···Cl hydrogen bound engagements between SbCl 3 and 18-crown-6. However, it should be kept in mind that these are not just the interactions responsible for the stability of the infinite crystals. There are a number other interactions between the host and guest species, including H . . . O hydrogen bonds, which can be revealed when the unit-cell is expanded.
On the other hand, the intermolecular bonding modes of crown-ethers with SbF3 are somewhat different to those of the other antimony trihalides, no doubt because of the small, highly electronegative nature of fluorine. Host-guest fitting is more prominent for SbF3 with the bound F atoms of the molecule involve in additional Sb-F···O interactions when a large crown-ether is the host (Figure 28a-e). These interactions are either comparable to, or slightly longer than, the F-Sb···O interactions, but are markedly longer than the Sb-F bonds (r(Sb-F) ~ 1.9 Å). The former are Type-I and the latter are Type-IIa interactions, as deduced from our investigation of the ∠Sb-F···O and ∠F-Sb···O contact angles in the solid-state structures of SbF3 complex with a number of 18-crown-6 crown ethers (Figure 28d,e). On the other hand, the intermolecular bonding modes of crown-ethers with SbF 3 are somewhat different to those of the other antimony trihalides, no doubt because of the small, highly electronegative nature of fluorine. Host-guest fitting is more prominent for SbF 3 with the bound F atoms of the molecule involve in additional Sb-F···O interactions when a large crown-ether is the host (Figure 28a-e). These interactions are either comparable to, or slightly longer than, the F-Sb···O interactions, but are markedly longer than the Sb-F bonds (r(Sb-F)~1.9 Å). The former are Type-I and the latter are Type-IIa interactions, as deduced from our investigation of the ∠Sb-F···O and ∠F-Sb···O contact angles in the solid-state structures of SbF 3 complex with a number of 18-crown-6 crown ethers (Figure 28d,e).
We note further that in all these host-guest structures, the importance of secondary interactions cannot be overlooked. They appear either as H···X or F···O interactions between SbX 3 and crown-ether. The presence of H···X hydrogen bonds in the host-guest complex systems can be inferred from the space-filling models shown in Figure 28d; each covalently-bound F in SbF 3 kisses at least two nearest H atoms of the -CH 2 fragment of 18-crown-6, forming a set of six hydrogen bonds. Of course, the number and extent of the formation of hydrogen bonds is determined by the size of the crown-ether. That SbF 3 molecular entity fits well inside the 18-crown-6 cavity is also evident in Figure 28d (see the space filling model, extreme left).
We note further that in all these host-guest structures, the importance of secondary interactions cannot be overlooked. They appear either as H···X or F···O interactions between SbX3 and crown-ether. The presence of H···X hydrogen bonds in the host-guest complex systems can be inferred from the space-filling models shown in Figure 28d; each covalently-bound F in SbF3 kisses at least two nearest H atoms of the -CH2 fragment of 18crown-6, forming a set of six hydrogen bonds. Of course, the number and extent of the formation of hydrogen bonds is determined by the size of the crown-ether. That SbF3 molecular entity fits well inside the 18-crown-6 cavity is also evident in Figure 28d (see the space filling model, extreme left).
By contrast, in the crowded environment of (NEt4) [(η 5 -cyclopentadienyl-Fe(CO)2)(Fe(CO)4)2SbI] [173], tetra-coordinated Sb 3+ also participates in Sb···π(C≡O) interactions ( Figure 29c). There are six such interactions, with each pair developing along the non-linear extension of each Sb-X (X = I, Fe) bond. The Sb···π(mid of point of C≡O) distance is 3.399 Å, and Sb···C(≡O) distance is 3.171 Å. From an examination of the spacefilling model of the system, it was inferred that the Sb cation was also linked non-covalently with the near C sites of the Fe-bound η 5 -cyclopentadienyl moiety.
By contrast, in the crowded environment of (NEt 4 ) [(η 5 -cyclopentadienyl-Fe(CO) 2 ) (Fe(CO) 4 ) 2 SbI] [173], tetra-coordinated Sb 3+ also participates in Sb···π(C≡O) interactions ( Figure 29c). There are six such interactions, with each pair developing along the non-linear extension of each Sb-X (X = I, Fe) bond. The Sb···π(mid of point of C≡O) distance is 3.399 Å, and Sb···C(≡O) distance is 3.171 Å. From an examination of the space-filling model of the system, it was inferred that the Sb cation was also linked non-covalently with the near C sites of the Fe-bound η 5 -cyclopentadienyl moiety.
To verify these conclusions, we examined the IGM-δg isosurfaces. Two examples are shown in Figure 30. The top one, Figure 30a, is for [1.5]dibenzothia-18-crown-6)SbF 3 [170], while the bottom one, Figure 30b, is for (HDTOA) 3 (SbI 3 ) 2 (HDTOA = N,N'-dicyclohexyldithiooxamide) [171]. The IGM-δg based isosurface results are in agreement with the inferences drawn from the intermolecular bond distances between the interacting molecules (Figure 29a,b) and the MESP analysis of SbI 3 (see Figure 2d). In the case of [1.5]dibenzothia-18-crown-6)SbF 3 , of the six Sb···O contacts, the two with the shortest distance are genuine coordinate bonds, as evidenced by the bluish-green isosurfaces. The other four Sb···O contacts are associated with greenish isosurfaces and are seemingly pnictogen bonds. In addition, we observed F···O, C···F and H···F interactions between SbF 3 and the crown ether. All the Sb-I and Sb-S contacts in (HDTOA) 3 (SbI 3 ) 2 are coordinate bonds, since they feature bluish-green isosurfaces. We have also observed a number of (N)C···I/(C)C···I and H···I contacts between the interacting monomer units. These results are in agreement with what might have been speculated from the space-filling models that rely on the vdW radii of the atomic domains. (Figure 29a,b) and the MESP analysis of SbI3 (see Figure 2d). In the case of [1.5]dibenzothia-18-crown-6)SbF3, of the six Sb···O contacts, the two with the shortest distance are genuine coordinate bonds, as evidenced by the bluish-green isosurfaces. The other four Sb···O contacts are associated with greenish isosurfaces and are seemingly pnictogen bonds. In addition, we observed F···O, C···F and H···F interactions between SbF3 and the crown ether. All the Sb-I and Sb-S contacts in (HDTOA)3(SbI3)2 are coordinate bonds, since they feature bluish-green isosurfaces. We have also observed a number of (N)C···I/(C)C···I and H···I contacts between the interacting monomer units. These results are in agreement with what might have been speculated from the space-filling models that rely on the vdW radii of the atomic domains. Figure 30. IGM-δg based isosurface plots for (a) [1.5]dibenzothia-18-crown-6)SbF3 [170], and (b) (HDTOA)3(SbI3)2 (HDTOA = N,N'-dicyclohexyldithio-oxamide) [171]. Atom labeling is shown. Isosurfaces colored blue and green between the interacting molecules in the crystal signify strong and weak intermolecular interactions, respectively. Other intermolecular interactions (viz. hydrogen bonds, C···F, F···O and C···I bonds, etc.) are marked with dotted lines in (b) and circles in (a). The space-filling model is shown for each case.

Conclusions
In this overview, we have examined the nature of the electrostatic surfaces of several molecular entities in various illustrative crystal lattices to provide evidence of the occurrence of an often overlooked non-covalent interaction, the antimony-centered pnictogen Figure 30. IGM-δg based isosurface plots for (a) [1.5]dibenzothia-18-crown-6)SbF 3 [170], and (b) (HDTOA) 3 (SbI 3 ) 2 (HDTOA = N,N'-dicyclohexyldithio-oxamide) [171]. Atom labeling is shown. Isosurfaces colored blue and green between the interacting molecules in the crystal signify strong and weak intermolecular interactions, respectively. Other intermolecular interactions (viz. hydrogen bonds, C···F, F···O and C···I bonds, etc.) are marked with dotted lines in (b) and circles in (a). The space-filling model is shown for each case.

Conclusions
In this overview, we have examined the nature of the electrostatic surfaces of several molecular entities in various illustrative crystal lattices to provide evidence of the occurrence of an often overlooked non-covalent interaction, the antimony-centered pnictogen bond, or more succinctly, the stibium bond. We must emphasize that we do not claim that our survey of the CSD and ICSD databases is comprehensive; we have merely drawn illustrative examples from those databases to highlight the occurrence of overlooked non-covalent interactions.
Pnictogen bonding in crystal lattices has been known since the middle of the last century, so it is misleading to claim that it only became known to community after 2011.Our investigations of intermolecular distances, directionality, electrostatic potential and IGM-δg based isodensity surfaces in most of the selected systems have enabled to us demonstrate that hypervalent Sb atoms in molecular entities contain more than one σ-hole and have the potential for involvement in attractive engagements with the negative sites in neighboring molecular entities, thereby causing (or contributing to) the assembly and stability of crystal lattices. The MESP model is suitable for understanding of directional Type-II interactions that are Coulombic, but not for rationalizing Type-I or -III interactions; this may be overcome by applying the IGM-δg model. The "sum of the vdW radii" concept, together with other geometric and topological features, has shown to be very effective in detecting pnictogen bonding in crystals when applied carefully.
Depending on the molecular environment, the directionality of these interactions may appear in different flavors, following Type-I, -II and -III bonding topology interactions. Our investigation suggests that pnictogen bonds are not always linear, or quasi-linear; indeed, non-linear Type-II interactions were abundant in many of the crystals we examined. In most cases, other primary or secondary interactions (such as hydrogen bonding, halogen bonding, tetrel bonding and other π-centered interactions) occur simultaneously, which drive the directionality of the stibium bonds. The involvement of the same covalently bound Sb atom in multiple interactions (for example, as in crown-ether complexes) and the competition between them influence the directionality of the pnictogen bonds.
We have demonstrated that antimony in molecules can form pnictogen bonds when it finds itself in the close vicinity of Lewis bases such as O, N, F, Cl, Br, I and S in molecules with which it interacts. We have also shown that pnictogen bonds can be formed with the delocalized π-density on the fragments -C=N=S and -C≡O in molecular entities, as well as on arene moieties. These findings suggest that donors of pnictogen bonds are no different to those that have been widely examined for hydrogen bonds, halogen bonds, and other non-covalent interactions. None of the crystal systems showed evidence that stibium bonds were the only attractive interactions causing molecular assembly and therefore contributing to the entirety of crystal stability.
We have also observed that intramolecular stibium bonds in crystals are not so uncommon, even though intermolecular stibium bonds are widely present in a variety of crystals. While the majority of the illustrative systems presented in this overview have yet to be thoroughly explored using state of the art computational methods and other approaches, we expect that this work will serve as a guide for researchers who use, or intend to use, pnictogen bonding in their work. Since Sb has already played a significant role in the synthetic design of functional zero-, one-, two-and three-dimensional nano-scale materials, especially in the area of optoelectronics, the possibility of exploiting the Sb-centered pnictogen bonding interactions discussed in this study should assist future investigations into these and other materials.