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

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


Introduction
A plethora of self-assembled crystal shapes, synthetically developed through the remarkable adhesion engineering of non-covalent interactions, have been widely discussed and cataloged in the Cambridge Structural Database (CSD) [1][2][3] and Inorganic Crystal Structure Database (ICSD) [4-6], among others [7].Halogen bond (or HaB for short) [8], which bears a striking resemblance to hydrogen bonding [9] and other non-covalent interactions [7,10,11], stands out as a prominent subset of non-covalent interactions.It was coined around 1978 [12], vigorously came to light after 2003 [13,14], and was integrated into the broader class of σ-hole interactions after 2007 [15].Halogen in molecular entities can form close contacts with atoms in another species with which it interacts, which can be seen in studies that have appeared since the late 1990s in the gas phase [16] and since the mideighteenth century in crystals [17][18][19].Its dynamical bonding aspect was elusive until Odd Hassel was awarded the Nobel Prize in Chemistry in 1969 for recognizing its fundamental importance in the understanding of charge transfer phenomena in complexes [20].
Halogen bonding [8], regardless of its strength, develops in a chemical system when there is evidence of a net attractive interaction between an electron-density deficient electrophilic region on the electrostatic surface of a halogen atom in a molecular entity and a close-lying electron-density rich nucleophilic region on the electrostatic surface of the same or another identical or different molecular entity [21].A number of its translucent caveats and features appeared in the non-covalent chemistry literature may be suitable for the recognition of halogen bonding in chemical systems.The IUPAC definition of the halogen bond [8] mirrors the definitions of the hydrogen bond [9], chalcogen bond [7], and pnictogen bond [10,22]; the only difference between the former and the latter three is that the terms such as "halogen" and "halogen bond" in the former were replaced by "hydrogen, or chalcogen, or pnictogen" and "hydrogen bond, or chalcogen bond, or pnictogen bond" in the latter.In other words, the definition of the halogen bond [8] is nothing but a writing variant of the definition of the hydrogen bond [9], with the term "hydrogen" in the latter being replaced by "halogen" in the former.The study that proposed the definition of the halogen bond has narrowed down a number of accompanying features and notes recommended for hydrogen and chalcogen bonds, yet all these definitions are essentially based on the same underlying concept: "An electrophile on a covalently bonded halogen/hydrogen/chalcogen/pnictogen in a molecular entity attracts a nucleophile on the same/another molecular entity".The base of the concept is clearly Coulombic.
Studies have shown that HaB is electrostatically driven [23][24][25].Contributions arising from the exchange-repulsion, polarization, and dispersion interactions meticulously play an important role in explaining the net binding energy of HaB in HaB-driven complexes [25][26][27][28][29]. Arguments in support of an orbital-based charge transfer interaction that drives the formation of halogen bonds have also been developed and appeared in journals in a timely manner [30][31][32][33].The debate continues, focusing largely on the similarities and differences between the role played by polarization and charge transfer [31,34,35], emerging from the blooming of Nobel Laureate Odd Hassel's chemistry [35].Additionally, it is commonly believed that halogen bonds are highly directional [36,37], with the R-X•••Y (X = halogen; Y = nucleophile; R = remainder part of the molecule R-X) bond angle being typically close to 180 • [35], and the linearity of the halogen bond has been argued to be better explained by charge-transfer interactions and lone pair repulsion [38,39].As shown in this study, the view is narrow.This is because halogen bonds can potentially be nonlinear depending on the mode of the interaction of covalently bonded halogen with the nucleophile; thereby, the underlying refulgent phenomena of halogen bonding are yet to be fully appreciated [40].
The key to recognizing HaBs in molecules or crystal systems is to identify the halogen bond acceptors (HaBAs) and halogen bond donors (HaBDs) that interact due to attractive forces toward the molecular assembly.Halogen bond acceptors are locally nucleophiles (Lewis bases) (viz.N-ends of N 2 ), and HaBDs are locally electrophiles (Lewis acids) (viz.X-ends of X 2 (X = F, Cl, Br, I)).Clearly, the rudimentary force that combines the HaBA and HaBD to shape a crystal is a result of the net attraction driven by their opposite charge capacities, in line with the IUPAC definition of the HaB [8] and our revisit to the definition [21].Insight into the formation of HaBs can also be revealed in terms of the conceptual theoretical framework of Pearson's hard and soft acids and bases (HSAB) [41,42], corroborating the importance of the fundamental Coulomb's law in electrostatics in that opposite charges on interacting atomic/molecular entities can attract each other when in close proximity.
In a simple chemical system, such as the hydrogen fluoride (HF) molecule, it may be immediately axiomatic that H and F have opposite capacities of charge (H δ+ and F δ− ), and the halogen-end of the molecule cannot be a HaBD [43].As the number of atoms that make up an arbitrary molecule increases, it is not very straightforward to readily capture the complicated distribution of charge density profile on the surfaces of the constituent atomic basins.This is because the distribution of charge density is generally anisotropic on the surface of atoms in a molecular entity; some parts of covalently bonded atom X in R-X may have insufficient electron density, while other positions of the same atom may have abundant electron density.This means that the electrophilicity or nucleophilicity of a region on X in R-X is determined not only by the electrons and nuclei in that region but also by the electrons and nuclei in the rest of the molecule, especially those in neighboring portions.If the electron density deficiency on atom X is sufficiently large, then atom X has what is commonly referred to as carrying a "hole" and may have a positive charge capacity (δ + ) locally and is, therefore, electrophilic.Such electrophilic regions have been identified in large numbers on the side of X opposite to the covalently bonded halogen derivatives in molecular entities.They have been called positive σ-holes (pronounced "sigma-hole") [15,44,45], or even referred to as electrophilic caps [46], and are prone to form HaBs when placed in close proximity to nucleophiles on other identical or different partner molecules.
The name "σ-hole" emerges from the notion that charge density deficiency on the covalently bonded atom X appears along the extension of the σ covalent bond, but, literally, it has nothing to do with a "hole" (the charge density deficiency is the "hole") [17].A negative σ-hole can be found on an atom X in R-X when it locally carries a negative charge capacity [47], as on F in H-F [48], H 3 C-F [49,50], and C 6 H 5 F [47]; it is not prone to HaB with a nucleophile on another identical or different interacting molecule.However, it can host itself as a nucleophile for an electrophilic σ-hole on halogen on another molecule when in close proximity, resulting in the formation of a σ-hole halogen bond.Thousands of intensive studies have been reported to identify and characterize halogen bonds formed by electrophilic σ-holes on halogen derivatives in molecules.When halogen atoms in molecules are hypervalent, they may host themselves as biaxial (σ-hole) halogen bond donors [28,51].However, the electrophilicity of a π-hole on the halogen derivative in molecular entities and the diversity of its halogen bond-forming ability have not been carefully addressed; this is likely to be the subject of intense research for the coming decade.
A π-hole is a region of abundant or deficient electron density on the surface of a covalently bonded atom, or array of atoms, that is orthogonal to the bonding direction, or plane, of the molecule.Therefore, a π-hole can be either electrophilic or nucleophilic [47,52,53].While π-holes correspond to a flattening of the electronic density surface, they are craters above and below the rings in benzene derivatives and in polyazines [52].For example, the π-holes on the electrostatic surface of either side of the centroid of the C 6 carbon ring of benzene and fluorobenzene (C 6 H 6 [54] and C 6 H 5 F [47], respectively) are nucleophilic and, thus, are referred to as negative π-holes [47,53].If all six hydrogen atoms of benzene are replaced by six fluorine atoms, positive π-holes are developed above and below the ring, and thus, C 6 F 6 has two positive π-holes [47,53].Similarly, both HCCH and N 2 have triple bonds, which were conventionally understood as electron-rich in the C≡C and N≡N internuclear regions.However, the C≡C bonding region in H-C≡C-H has a negative belt-like π-hole, while the N≡N bonding region in N 2 has a positive belt-like π-hole [55,56].
A recent study has demonstrated that π-holes in molecules can be visualized by means of Kelvin probe force microscopy [57], although the chemical binding capacity of negative π-holes has been known over several decades [29,58] and that of positive π-holes on main group elements has only begun to be studied within the current decade [59][60][61].It was argued that if both electrostatics and polarization are taken into account, counterintuitive interactions involving π-holes can be treated as Coulombic [58,62], despite the apparent repulsion between ground-state molecules [34,37].Opposite views have also been advanced [27,63].The most commonly observed chemical systems where both σand π-holes are integrated are halogenated aromatic compounds.
Does an electrophilic π-hole exist on the surface of halogen derivative in molecules?Does it form a halogen bond when it is in the vicinity of a nucleophile on the partner molecule with which it interacts?To our knowledge, this question has not been properly addressed before, and the answer to both questions is "yes, it does exist" and "capable of halogen bond".This is demonstrated in this paper with examples of crystals from CSD [1][2][3], known for many years, in which the π-hole interactions in them were neither named nor properly characterized nor assigned as π-hole halogen bonds.To this end, we have placed our focus on some simplified halogenated chemical systems that radiate π-holes on their electrostatic surfaces.We reveal this based on the application of the MP2 level of theory and the molecular electrostatic surface potential (MESP) tool to the XY 3 (X = Cl, Br; Y = F, Cl, Br) set of molecules.We then demonstrate whether the electrophilic π-hole (Lewis acid), for instance, in BrCl 3, is capable of forming π-hole halogen bonds when interacting with a series of nucleophiles (Lewis bases) hosted individually by 32 other partner molecules and the power of chemical reactiveness of the π-hole on covalently bonded halogen as π-hole halogen bond donors.We have utilized the quantum theory of atoms in molecules (QTAIM) [64][65][66] and the Independent Gradient Model (IGM-δg inter ) [67,68] approaches to characterize the π-hole halogen bonds, including a discussion of their energy strengths and geometric signatures.A higher-order symmetry-adapted perturbation theory (SAPT2+(CCD)) [69,70] was also applied to provide insight into the nature of dissected energy components that may explain the dominant factors responsible for the interaction energies of the π-hole halogen bonded complexes investigated, where the many-body treatment of dispersion was based on coupled-cluster doubles (CCD) [71].The NBO's [72] second-order perturbative estimates of 'donor-acceptor' (bond-antibond) interactions in the NBO basis [73] were discussed for some complexes to evince the possibility of hyper-conjugative charge transfer interactions between interacting monomers.

π-Hole Halogen Bond Donor Molecules
To illustrate the existence of a π-hole on halogen atoms in molecules such as ClF 3 , ClBr 3 , and BrCl 3 , two different geometries for each, planar (C 2v ) and T-shaped (D 3h ), were fully relaxed.The T-shaped configuration (left, Figure 1) is energetically favorable over the trigonal planar configuration (right, Figure 1) for each of the three cases.Both configurations represent a stationary point based on the positive sign of the eigenvalue of the Hessian matrix.The relative energy between the two configurations is ca.16.44, 1.44, and 2.11 kcal mol −1 for ClF 3 , ClBr 3 , and BrCl 3 , respectively, suggesting that the planar geometry of ClF 3 (point group, D 3h ) is unstable and lying far above the stable minimum.
Figure 1a-f shows the MESP graphs for both the configurations of ClF 3 , ClBr 3 , and BrCl 3 , each superimposed with its corresponding QTAIM's molecular graph.The nature of the prominent electron density deficit is revealed on the surface of the central halogen atom for the most stable configuration (Figure 1a,c,e).Those regions are colored blue or cyan, green, or a combination of either of the two, where the extent of electron density deficiency is appreciable, characterized by the positive sign of the potential.From Figure 1a,c,e, it can be captured that the strength of the π-hole follows the following order: ClF 3 (Tshaped) > BrCl 3 (T-shaped) > ClBr 3 (T-shaped), revealed by V S,max > 0 for the former two and V S,min > 0 for ClBr 3 .Conversely, the sign of potential is negative for a pair of σ-holes on the two halogen atoms responsible for the linearity of the Y-X-Y (X, Y = halogen) skeleton in ClF 3 (Figure 1a) and BrCl 3 (Figure 1e) and positive for ClBr 3 (Figure 1c), signifying the presence of a lack of electron density deficiency on the surface of halogen in the latter.The stronger σ-holes are located along the outermost extensions of the remaining Cl-F/F-Cl bonds in ClF 3 (Figure 1a), Cl-Br/Br-Cl bonds in ClBr 3 (Figure 1c), and Cl-Br/Br-Cl bonds in BrCl 3 (Figure 1e), characterized by V S,max >> 0. The lateral surface portions of the two linearly arranged halogen atoms in the T-shaped geometries of the three molecules are electron-density rich, which may act as nucleophiles for HaBDs.
The electrostatic surfaces of the central Br and Cl atoms of the three planar molecules have six maxima and two minimum of potential (marked by tiny red and blue circles, respectively) orthogonal to the Br-Cl and Cl-Br/Cl-F bonding axes, respectively.A set of four extrema (three maxima and one minima) lying above the molecular plane is shown (Figure 1b,d,f), while the remaining set lying below the same plane is not shown (opposite to the viewer).Three of the four in a set appear in the junction regions (as potential maxima), each between a pair of adjacent halogen atoms, and the minimum of potential appears on the centroid surface area above the central halogen atom.They are each a representative electrophilic π-hole given the charge density is sufficiently deficient on the surface of these molecules, characterized by the positive sign of V S,min or V S,max (V S,min , V S,max > 0).(electrons bohr −3 ) isoelectronic density envelope was used to map the potential.The location of σand π-holes on the surfaces of the first two molecules (a-d) are indicated by arrows in black.The regions of local most maxima and minima of potential (V S,max and V S,min ) are marked by filled tiny circles in red and blue, respectively.Values are in kcal mol −1 .The exact number of minima and maxima of potential are also on the opposite side of each of the three planar molecules (b,d,f), the side opposite to the viewer.
The energy of the strongest π-hole, characterized by V S,min , is ca.17.4 kcal mol −1 in ClF 3 (Figure 1b), 10.8 in kcal mol −1 in ClBr 3 (Figure 1d), and 15.3 kcal mol −1 BrCl 3 (Figure 1f).The strength of the π-hole in these molecules described by V S,max follows the order: ClF 3 (23.7 kcal mol −1 , Figure 1b) > BrCl 3 (17.9kcal mol −1 , Figure 1f) > ClBr 3 (11.4kcal mol −1 , Figure 1d).The trend is reasonable since F in ClF 3 is highly electronwithdrawing than Cl and Br in ClBr 3 and BrCl 3 , respectively, thereby creating a relatively strong electron density deficient region on the surface of Cl in ClF 3 .This may also be explained in terms of the increase in polarizability of the halogen derivative, which follows the order: Br > Cl > F. In all three cases, the lateral surface portions of the halogen derivative along the outer extensions Cl-F, Cl-Br, and Cl-Br covalent bonds are nucleophilic (characterized by the negative sign of V S,min ).The σ-holes located at the outer extension of the Cl-F bonds in ClF 3 are weakly nucleophilic; it is is more so for lateral portions of the same halogen atom.We did not observe a σ-hole on the surface the halogen along the F-Cl, Br-Cl, and Cl-Br covalent bond extensions in the geometries shown in Figure 1b,d,e, respectively.Nevertheless, the results above indicate that hypervalent halogen derivatives in halogenated molecules can not only be anisotropic but can host σand π-hole donors simultaneously.

Binary Complexes: Geometries and Bonding Features of π-Hole Halogen Bonds
To demonstrate the development of π-hole halogen bonding interaction between a pair of molecular entities, the T-shaped configuration of BrCl 3 as the π-hole halogen bond donor and a series of Lewis bases as π-hole halogen bond acceptors were considered.The nucleophiles in the 32 HaBAs are either anions or neutral molecules.The Lewis base appears either on the entire electrostatic surface of the atom (as in the anions X − (X = F, Cl, Br, CN) and OX − (X = F, Cl, Br)) or around the bonding region, or on a small portion of an atom that constitutes the molecular entity.For instance, the nucleophile is concentrated around the bonding regions in F 2 , Cl 2 , Br 2 , ClF, BrF, and BrCl, among others, as well as on the lateral portions of the bonded atoms forming them.The outer portions on the sides of the halogen along the halogen-halogen bond extensions in these diatomic molecules are equipped with positive σ-holes.
Figure 2 (1-33) provides the theoretical evidence that the Lewis bases on the molecular entities above, characterized by negative electrostatic potentials, can show an aptness for attracting the electrophilic π-hole on the surface of Br in BrCl 3 .This causes the formation of π-hole halogen bonding interactions in the binary complexes investigated.
The weakest 1:1 binary complex in the series is observed between Cl 3 Br and H 2 (Figure 2 (1)).Two types of intermolecular interactions may be speculated from the geometry of this complex.The outer nucleophilic part of the bonding region of the H 2 molecule may be attractively engaged with the electrophilic π-hole on Br in Cl 3 Br, forming a Br(π)•••H 2 halogen bond.Similarly, an H atom in H 2 may be engaged with the nucleophile on the nearest Cl atom in BrCl 3 , forming the H•••Cl hydrogen bond.The occurrence of these two interaction types in the complex is not very surprising since the nucleophile on H 2 is delocalized around the bonding region described by a belt-like negative potential, and the end portions of each of the two H atoms of the same molecule are reasonably charge density deficient and characterized by positive electrostatic potentials [74,75].Clearly, the attraction between the interacting regions between the two molecules leading to the formation of the Cl 3 Br(π)•••H 2 dimer is primarily a result of Coulomb forces.
The intermolecular distances in Cl 3 Br They fail the distance-based feature for hydrogen [9] and halogen bonding recommended by IUPAC [8], given that neither of the intermolecular distances is less than the sum of the van der Waals (vdW) radii of respective atomic basins.For instance, the intermolecular distance for the Br•••H contact is ca.3.391 Å and is not less than vdW radii sum of H and Br atoms, 3.06 Å (r vdW (H) = 1.20 Å; r vdW (Br) = 1.86 Å [76]).Similarly, the intermolecular distance for the H•••Cl hydrogen bond is 3.089 Å and is not less than van der Waals radii sum of H and Cl atomic basins, 3.020 (r vdW (H) = 1.20 Å; r vdW (Cl) = 1.82 Å [76]).As discussed by us and others [10,11,17,77], this type of failure of the IUPAC's recommended feature is not surprising given that the vdW radii of atoms are not exact and have involved several approximations for their determinations and are associated with an error that can be as large as ±0.20 Å [10,78].Note that the bond path between Br(π-hole) and H2 is inwardly curved near H2, which indicates possible involvement of the bonding region in H2 with Br's π-hole in developing the Br(π-hole)•••(H-H)mid-point closed-shell interaction.This is also confirmed by IGM-δg inter , revealing a flat-type distored dumbbell-shaped isosurface domain between the interacting moieties that suggests the presence of attractive interaction.The BSSE corrected interaction energy for the dimer is ca.−0.85 kcal mol −1 (see Table 1), leading to the rationale that Cl3Br•••H2 is a vdW complex (vdW complexes usually have interaction energies < −1.0 kcal mol −1 [10,11]).NBO's second-order perturbation theory revealed the H     Note that the bond path between Br(π-hole) and H 2 is inwardly curved near H 2 , which indicates possible involvement of the bonding region in H 2 with Br's π-hole in developing the Br(π-hole)•••(H-H) mid-point closed-shell interaction.This is also confirmed by IGM-δg inter , revealing a flat-type distored dumbbell-shaped isosurface domain between the interacting moieties that suggests the presence of attractive interaction.The BSSE corrected interaction energy for the dimer is ca.−0.85 kcal mol −1 (see Table 1), leading to the rationale that Cl 3 Br•••H 2 is a vdW complex (vdW complexes usually have interaction energies < −1.0 kcal mol −1 [10,11]).NBO's second-order perturbation theory revealed the H•••Cl hydrogen bond is a characteristic of np(Cl) → σ*(H-H) charge transfer delocalization, with E (2) = 0.95 kcal mol −1 , where np and σ* refer to the p-type lone-pair orbital and σtype anti-bonding orbital, respectively.A similar type of charge transfer delocalization (σ(H-H) → RY* (10) Br; E (2) = 0.12 kcal mol −1 ) was yielded to describe the Cl 3 Br(π-hole)•••H 2 halogen bond, where RY* is an anti-bonding orbital called extra valence shell Rydberg orbital [82].When N 2 and CO were used as HaBAs, we did not observe any kind of sleep parallel arrangement between either of them and the π-hole on Br in BrCl 3 .The energyminimized geometries of Cl 3 Br(π-hole) 2 (2), ( 5) and (4), respectively.The spatial arrangement between the monomers in these dimers is logical since the central portion of N 2 /CO is charge density deficient and the nucleophile on them is mainly concentrated around the end portions of the N/C/O atoms [56,83].The former causes steric repulsion between the bonding region in N 2 /CO and the electrophilic π-hole on Br in BrCl 3 , pushing them apart, aligning to an orientation so as to maximize the coulombic interaction with the electrophilic π-hole on the outer surface of the Br atom in BrCl 3 .Although QTAIM typifies the presence of the localized interaction between N/C/O in N 2 /CO and Br in BrCl 3 (see Figure 3 (2), ( 5) and ( 4)), IGM-δg inter 's isosurface plot indicates a semi-elliptical-and triangle-shaped flat-type interaction between them that vivifies the possible presence of secondary interaction (see Figure 4 (2), (5), and ( 4 The perspective above may not be very surprising since the lateral portion of P in PN [43] and S in CS are positive.This may be inferred from the MESP plot shown in Figure 5a,b for these two molecules, respectively.As such, the V S,min and V S,max are −26.3 and 20.5 kcal mol −1 on the outermost extensions of the S=C and C=S bonds in CS, respectively, and the C-end along the outermost extension of the S=C bond is entirely negative.Similarly, it is conventionally understood that the electron density is concentrated around the P≡N bond region and that N is entirely negative.This is not the case; a ring of maxima is seen around the P atom (V S,max = 25.1 kcal mol −1 each), and the negative electron density is far from the triple bond and is concentrated at the N terminus along the outermost extension of the P≡N triple bond.There is no maximum seen on the outermost extension of the P≡N bond, but there is an electron density deficit on P that appears probably to be due to a p-type orbital of the same atom and is described by a positive V S,min (V S,min = 22.5 kcal mol −1 ).The local electrophilic and nucleophilic pattern of this molecule may somehow resemble that of N 2 [55,56], given that both molecules belong to the same pnictogen family.Nevertheless, the surface features of CS and PN are consistent with the bonding environment in the Cl 3 Br•••CS and Cl 3 Br•••NP dimers, in which Coulomb forces act between monomers that are presenter sites of opposite charge capacity.These also explain the reason why primary and secondary interactions are developed between the interacting monomers in the equilibrium structures of both dimers.
Int. J. Mol.Sci.2024, 25, x FOR PEER REVIEW 13 of 27 bond, but there is an electron density deficit on P that appears probably to be due to a ptype orbital of the same atom and is described by a positive VS,min (VS,min = 22.5 kcal mol −1 ).
The results of QTAIM and IGM-δg inter justify the presence of both interaction types in the complexes involving LiX (X = H, F, Cl, Br) (( 22), ( 26), (27), and (28) in Figures 3 and 4, respectively), NaX (X = H, F, Cl, Br) (( 17), ( 25), (24), and (23) in Figures 3 and 4 The HF molecule showed different reactivity when placed near BrCl 3 compared to HCl and HBr.The latter two molecules exhibit very similar intermolecular bonding patterns toward the formation of Br(π-hole)•••X (X = Cl, Br) halogen bonds (see Figure 2 (14) and (12), respectively).Because of the strong electrophilic dispersion of the H atom due to the presence of the hemispherical σ-hole [88] and the negative hemisphere of F, the tiny HF molecule was unable to stand when its F-end was placed on the top of the Br's π-hole.It flitted toward the most positive site on Br in BrCl 3 during energy minimization, and in the fully relaxed geometry, the HF molecule is stable at a position shown in Figure 2 12), respectively.This limitation of the space partitioning approach has been reported previously, especially in the weak bonding regime [49,[89][90][91].
The anion-molecule interactions observed in Cl  29) in Figure 3 and 4, respectively.The former predicts a localized closed-shell interaction between the interacting molecules in both the dimers of CN − and misses a bond path that could describe the Cl (20)).This was not the case with the results obtained from the IGM-δg inter analysis, which suggested a flat-type distorted-dumbbell-shaped interaction between the monomers causing the stability of the two dimers (Figure 4 (19) and ( 20)) and a possible involvement of the nucleophilic π-density of CN − .
The anion-molecule complexes formed between BrCl 3 and XO − (X = H, F, Cl, Br) (Figure 2  for the corresponding complexes, respectively.The hydroxyl anion, OH − , on the other hand, shows an analogous bonding proficiency with the electrophilic σ-hole on Br in BrCl 3 , giving rise to the shortest bond distance, r(Br•••O), of 1.973 Å.Not only is the intermolecular interaction in this complex highly directional (∠Cl-Br•••O = 179.9• ), but the BSSE-corrected interaction energy of −103.2 kcal mol −1 is the largest compared to all the anion-molecule complexes examined.The large interaction energy for this complex is not solely due to the Br•••O halogen bond but also due in part to the H•••Cl hydrogen bond.
The very large interaction energies for Cl 3 Br•••OX − (X = H, F, Cl, Br) suggest that the assembly of the Cl 3 Br molecule with the anions is a result of the formation of formal covalent bonds; a similar view that is also advanced in other studies for other complexes [33,92].The argument counts on the results of the IGM-δg inter analysis (Figure 4 (30-33)), in which the isosurfaces representing the Br•••O close-contacts are thickly shaped and colored blue, originating from the appreciable accumulation of charge density at the bond critical point region.QTAIM also gave the sign of H b to be negative (H b < 0) at the Br•••O bcps of the Cl 3 Br•••OX − , with appreciable character of ionicity (∇ 2 ρ b > 0), and the charge density at the corresponding bcps (ρ b > 0.1 a.u.) was not very small (covalent bonds typically have ρ b values ∼ 1.0 eÅ −3 [93]).This result is in contrast with that unveiled for the remaining 29 complexes (see  in Figures 2-4. Comparable interaction energies of anion-molecule complexes driven by σ-hole hydrogen and halogen bonds have been reported elsewhere [92,94].It was shown [92] that the range of strongly halogen-bonded trihalides D-X•••A − and the analogous strongly hydrogen-bonded complexes D-H•••A − (D, X, A=F, Cl, Br, I) examined can be explained by the combined effect of electrostatic and covalent components stemming mainly from the HOMO-LUMO interaction between the occupied halide np atomic orbital (AO) and the D-H or D-X (X = halogen) σ* anti-bonding acceptor orbital.For the halogen-bonded complexes DX•••A − , the contribution from the orbital interaction term was large, from 43% for F-I•••F − to 97% for I-F•••F − .The larger covalent contribution is due to the lower orbital energy of the empty dihalogen σ* orbital (e.g., −0.7 eV for F-H and −6.2 eV for F-F), which results in stronger donor-acceptor orbital interactions with the halide np orbital.
It is worth noting that the proximity of the nucleophile in HaBA has induced geometric deformation into the Lewis acid.It is not just the Br-Cl bonding region that hosts the π-hole, but all three Br-Cl bonds in BrCl 3 were distorted upon the formation of halogen bonds.For example, the Br-Cl bond giving BrCl 3 its T-shape was elongated for all dimers, except for Grabowski [95][96][97], as well as others [98,99], published a number of papers considering a few halogen bond donors (viz.BrF 3 ) shown in Figure 1.In ref. [96], he labeled the σ-hole interactions as new types, even though such interactions formed by hypervalent halogen derivatives have been known for decades [99].Similarly, in ref. [98], others have examined the complexes formed of H 2 Y (X = O, S, Se) with the hypervalent halogens XF 3 and XF 5 (X = Cl, Br, I).Their focus was to uncover the role of the electrophilic σ-hole on the X atom to form halogen bonding interactions with the lone-pair on the chalcogen atom in H 2 Y. Their results have shown that the interaction energies can be correlated with the properties of the X and Y atoms, that the formation of the complexes is dominated by electrostatics, and that the halogen bonds involving charge density donors H 2 S and H 2 Se possess some covalent character.Similar arguments have been provided by others [100].

SAPT-Based Energy Decomposition Analysis
To give substance to the notions developed above, SPAT2+3(CCD) level analysis was performed for all the binary complexes; the MP2 level optimized geometries were used.The results in Table 1 show that the energy contribution resulting from electrostatics (E els ) is always attractive, meaning that Coulomb interactions play an important role in driving the monomers toward association.This, together with the induction (E ind ) and dispersion (E disp ) terms, is complementary with the exchange-repulsion energy component (E exch ), giving rise to the net interaction energy, E[SPAT2+3(CCD)], for each dimer explored.
Figure 6b-d depicts the linear dependence of E[SPAT2+3(CCD)] on E els , E ind , and E exch , corresponding to R 2 of 0.9935, 0.9836, and 0.9864, respectively.However, the dependence of E[SPAT2+3(CCD)] on E disp was reasonably poor (R 2 = 0.6198).This is because the E disp could not be properly determined for most of the complexes.The large E disp for the anion-molecule complexes (see Table 1 for E disp values) gives the impression that the halogen-bonded interactions in them comprise some high degree of covalency.Nonetheless, Sengupta et al. [101] previously plotted each of the three attractive energy components against the [SAPT0/jun-cc-pVTZ] interaction energy for 37 small molecule complexes and observed that each component was linearly correlated.The correlation was strongest with the electrostatic component (R 2 = 0.99), similar to what is observed in this study (R 2 = 0.9935, Figure 6b).This has enabled the authors to demonstrate the overriding importance of electrostatics and the possibility of a predictive model in which the electrostatic energies can be considered as a match to the gas-phase total binding energies.
the electrostatic component (R 2 = 0.99), similar to what is observed in this study (R 2 = 0.9935, Figure 6b).This has enabled the authors to demonstrate the overriding importance of electrostatics and the possibility of a predictive model in which the electrostatic energies can be considered as a match to the gas-phase total binding energies.2 for geometries and Table 1 for energies).
From the energies of the first 11 complexes in Table 1 (Figure 2 (1-13)) and Cl3Br•••P2 (Figure 2 (13)), it may be seen that the dispersion term is dominant over the electrostatic component.This convinces us that, while dispersion certainly intensifies the strength of the interaction, the origin of the interaction is electrostatically driven, and polarization is an integral part of the Coulomb interaction [102][103][104][105][106][107][108].
For the remaining complexes, the scenario is not very different, meaning that the electrostatic contribution is the main driving force for their formation.The energy due to induction contributes appreciably to E[SPAT2+3(CCD)] for all complexes.However, for the first 15 complexes, the induction term dominates over the dispersion term, while for the remaining 18 complexes, the latter dominates over the former (Figure 2 (16)(17)(18)(19)(20)(21)(22)(23)(24)(25)(26)(27)(28)(29)(30)(31)(32)(33)).For most complexes, E[SPAT2+3(CCD)] is comparable with ∆E(BSSE) [MP2/aug-cc-pVTZ].The largest difference between them is found for anion-molecule complexes, Cl 3 Br•••X − (X = F, Cl, Br, OH), with deviations ranging from 2.01 to 22.81 kcal mol −1 .It arises from the fact that SPAT2+3(CCD) relies on a "simple" pair model of each complex, and the charge polarization induced by the interacting molecule on the halide/hydronium anions could not be properly accounted for by SPAT2+3(CCD); a similar view is advanced elsewhere [109].The deviational feature may be realized if E els and E ind are compared for each Cl 3 Br•••OX − (X = H, F, Cl, Br).Both the components are reasonably competitive with each other for these complexes, except for Cl 3 Br•••OH − , in which case, the former (i.e., E els = −167.51kcal mol −1 ) is about 8.32 kcal mol −1 smaller in magnitude than the latter (E ind = −175.83kcal mol −1 ).These results enable us to conclude that the anion-molecule in-teractions are indeed dominated by electrostatics, yet reinforcing the conventional wisdom of the importance of the induction term to anion binding.

Materials and Methods
Geometry optimization and frequency calculations were performed for all 33 dimers using Gaussian 16 software [110] employing the MP2 method [111,112] in conjunction with Dunning's type triple-ζ basis set aug-cc-pVTZ.The wave functions generated after geometry optimizations were used to calculate molecular graphs, bond critical point charge density properties, and molecular electrostatic surface potential graphs using AIMAll (version 19.10.12)[113].Multiwfn (version 3.7) [114] and VMD (version 1.9.3) [115] software were used to calculate IGM-δg inter isosurfaces between interacting monomers, as well as for visualization and graph generations.
The dimers reported in the former section are identified as stationary points, evidenced by positive normal mode vibrational frequencies.Table S1 of Electronic Supplementary Information lists the fully-relaxed redundant internal coordinates of all dimers.
The σand π-hole features on the molecule's electrostatic surface were characterized based on the signs of the local maximum and minimum of potential (symbolized by V S,max and V S,min , respectively), which are the two key descriptors of the MESP model [55,[116][117][118].That is, when regions of a molecular surface are characterized by V S,max < 0 (or V S,min < 0) and V S,max > 0 (or V S,min > 0), we recognize them as electron density rich and electron density deficient regions, respectively.Accordingly, a σ-hole that is often located on the surface of atom A along the outermost extension opposite to an R-A covalent bond can be characterized by a positive or negative sign of V S,max (V S,max < 0, or V S,max > 0).When the charge density is sufficiently deficient on the side of A opposite to the covalent bond extension (V S,max > 0), we label the region as a positive σ-hole [15,17,40]; when the deficiency in the charge density is not appreciable, and characterized by a negative V S,max (V S,max < 0), we label the region on A as a negative σ-hole [47,53,119,120].Chemical instances are known in which the outer extension of the covalent bond does not always feature a V S,max , but a V S,min .This is because the charge density is shifted from other parts of the molecule and accumulated around the outer extension of the bond (as on N in N 2 [56], C/O in CO [121], and N in HCN [121], for example).
The characteristics of the π-hole can be understood based on positive or negative signs of the electrostatic potential.A region on a molecular entity can have either V S,max > 0 (or V S,min > 0) if the π-hole is electrophilic (electron density deficient) [52,122,123], or V S,max < 0 (or V S,min < 0) if nucleophilic (electron density rich) [47,55,56].When such a situation is encountered, the former and latter characteristics signify the presence of positive and negative π-holes, respectively, and appear on the outer extensions around the covalently bonded atom (or molecular plane) orthogonal to the R-A bond axis [55].It is worth noting that positive σand π-holes are not a region of electrostatic potential but correspond to a region of appreciable deficiency in electron density [17].
The reason for exploring QTAIM's bond critical point features of charge density counts on IUPAC's proposed feature put forward for hydrogen bonding [9]; this was also recommended for identifying halogen bonds [8], chalcogen bonds [7], pnictogen bonds [10,22], and tetrel bonds [11], among others [124,125].The feature suggests that the charge density topology of the bond path and critical points could restore the shape of the molecule in terms of the molecular graph and that the sign and magnitude of the Laplacian of the charge density at the bond critical point, ∇ 2 ρ b , should indicate the degree of charge density depletion (∇ 2 ρ b > 0) and concentration (∇ 2 ρ b < 0).Ourselves [10,11] and others [95] have recently suggested that the total energy density (H b ) at the bond critical should also be scrutinized because it provides insight into the ionic/covalent (or mixed) nature of the chemical interaction under investigation.Considering that the bond critical point topology of QTAIM sometimes fails to appear in weak bonding regimes [49,89,90,126], an exploration of isosurfaces of charge density emerges from the application of reduced density gradient (RDG) [127] and independent gradient model (IGM-δg inter ) [67,68] may be necessary, given these approaches have proven faithful in detecting weakly bonded interactions in chemical systems of diverse origin [10,11,56,109,128,129].The theoretical details underlying these methodologies have been repeatedly discussed in several studies.
The uncorrected interaction energy, ∆E, of a dimer is calculated using the standard procedure given by Equation (1), where E T is the electronic total energy of representative species.The ∆E was corrected upon incorporating the effect of basis set superposition error (BSSE), E(BSSE), accounted for by the counterpoise procedure of Boys and Bernardi [130], given by Equation (2).
Given that there is much ongoing debate as to whether halogen bonding is an electrostatically driven phenomenon [24,44], a combined effect of electrostatics and (orbital) charge transfer [41,131], or electrostatics, exchange, polarization, and dispersion terms [40,132,133], we attempted to rationalize the dominant energy contribution to the interaction energy of the investigated dimers.We utilized the symmetry-adapted perturbation theory (SAPT) approach implemented in the PSI4 code [134], which has several variations, with zerothorder (SAPT0) being the simplest [70].The SAPT module can perform density-fitting evaluations of SAPT2, SAPT2+, SAPT2+(3), and SAPT2+3 energies only for closed-shell systems [69,70]; for methods such as SAPT2+, the many-body treatment of dispersion can be replaced by an improved method based on coupling cluster doubling (CCD) [71].
Whether CCD dispersion provides more accurate interaction energies tends to depend on the SAPT truncation and basis set employed due to the offsetting errors.We have employed the SAPT2+(CCD) level of theory, in conjunction with the basis set aug-cc-pVDZ, to provide our level of interpretation on the interaction energy of π-hole and σ-hole halogen bonded complexes, as reported elsewhere [135].The NBO calculations were performed at the restricted Hartree-Fock level of theory, in conjunction with the basis set aug-cc-pVTZ, to estimate the second-order perturbative energy, E (2) , between the donor and acceptor orbitals responsible for non-covalent bond formation [73,136,137].

Evidence of π-Hole Halogen Bonds Formed by Halogen Derivatives in Crystals: Yet to Be Fully Rationalized by Computational Approaches
A search of vertical contacts of halogen in the CSD [1][2][3] has resulted in a number of crystals.In them, the higher derivative of halogen retains charge density-deficient regions that are probably electrophilic π-holes.Figure 7a-f provides insight into this in six such crystals, for example.The non-covalent interaction (dotted lines colored red) between the building blocks in these crystals may be readily speculated upon inspecting the intermolecular geometry.The π-holes on covalently bonded I develop I•••O and I•••F type π-hole halogen bonds, contributing to the adhesion between the building blocks to shape the crystals.They are not only predominantly non-linear but also directional in the sense that the nucleophile directs toward the π-hole on halogen.
The intermolecular distances associated with the unique halogen bonds mentioned above are shorter than the vdW sum of the respective bonded atomic basins.Some of the exemplary crystal systems in Figure 7 have been known since the 1980s.They have been called novel iodine(III) compounds.However, the I•••O and I•••F close contacts formed by covalently bonded iodine atoms in these systems were not assigned as halogen bonds.For example, the bis(methyloxalato)-phenyl-iodine (Figure 7e) and 2acetyl-1-fluoro-1,2-dihydro-3H-1,2-benziodazol-3-one (Figure 7f) crystals have been known since 1989 and 2021, respectively, but the π-hole halogen bond in them was not appreciated.The authors of the study [143]  Previous studies [146] on hypervalent halogen atoms in neutral and anionic molecular entities containing bromine and iodine atoms have shown positive and negative σ-holes in their covalent bond extensions, respectively.Some of these crystal systems may also comprise vertical interactions between electrophilic iodine and interacting nucleophiles.However, the authors have suggested that crystallographic Some of the exemplary crystal systems in Figure 7 have been known since the 1980s.They have been called novel iodine(III) compounds.However, the I•••O and I•••F close contacts formed by covalently bonded iodine atoms in these systems were not assigned as halogen bonds.For example, the bis(methyloxalato)-phenyl-iodine (Figure 7e) and 2acetyl-1-fluoro-1,2-dihydro-3H-1,2-benziodazol-3-one (Figure 7f) crystals have been known since 1989 and 2021, respectively, but the π-hole halogen bond in them was not appreciated.The authors of the study [143]  Previous studies [146] on hypervalent halogen atoms in neutral and anionic molecular entities containing bromine and iodine atoms have shown positive and negative σ-holes in their covalent bond extensions, respectively.Some of these crystal systems may also comprise vertical interactions between electrophilic iodine and interacting nucleophiles.However, the authors have suggested that crystallographic analysis can reveal close contacts, but calculated electrostatic potentials can provide insight into the extent of variations in positive potentials on the cationic surfaces and negative potentials on the anionic surfaces.Similarly, Wieske and Erdelyi [147] discussed the bonding features of hypervalent halogen (I) complexes.They argued that the halonium ions bound to two carbons, [C-X-C] + , do not possess halogen-bond donor character, as discussed by many research groups in the past.However, when nitrogen, sulfur, oxygen, or halogen-donor Lewis bases (D) are involved, hypervalent/hypercoordinate halonium ion [D•••X•••D] + halogen-bond complexes are formed.With a view to this, it is unclear to us as to why the halonium ion bonded to the two carbons in [C-X-C] + cannot be considered as halogen bond donors.It contradicts, not only to the very cationic character of halonium, but also to its versatile halogen bond donating ability.Our view is in line with Cavallo and coworkers [146], who explicitly discussed the halogen bond donating capacity of a range of iodonium and bromonium cations (viz.diphenyliodonium and dibenzo[b,d]iodonium and di-p-fluorophenylbromonium cations).

Conclusions
This study elucidates the potential existence of π-holes on halogen atoms within molecules, capable of forming π-hole halogen bonds with nucleophiles in another, representing a sister non-covalent interaction to σ-hole halogen bonds.π-hole halogen bonds exhibit many common features with σ-hole halogen bonds and other noncovalent interactions, with the primary distinction lying in their directional preference and spatial nature of occurrence.While σ-hole halogen bonds are highly directional axial interactions along the outermost extension of the halogen bond donor, π-hole halogen bonds are orthogonal interactions that occur perpendicular to the covalent bond axis or molecular plane.The intermolecular interaction distances associated with π-hole halogen bonds may not necessarily be less than the van der Waals sum of bonded atomic basins.Additionally, we observed that in some instances, particularly when atomic or diatomic anions were employed as nucleophiles, the π-hole of the halogen atom, such as that on Br in BrCl 3 , failed to form π-hole halogen bonds with interacting negative sites.This reaffirms the hypothesis of Politzer and coworkers that regions of positive or negative charge capacities on a molecular entity do not always attract regions of opposite charge capacities on another interacting molecular entity.Notably, in cases of anion-molecular complexes, the anion was unstable on the surface of the π-hole region but was able to adjust its spatial arrangement to maximize non-covalent interactions with the most positive sites, forming σ-hole halogen bonds with the interacting molecule.Furthermore, the interaction energy of π-hole halogen bonds was observed to be competitive with what has widely been reported for σ-hole halogen bonds, which can be van der Waals, weak, strong, very strong, and super-strong types depending on the nature of the nucleophile interacting with the halogen bond donor.The σand π-hole donors of halogen may be comparable in strength to aerogen, pnictogen, chalcogen, and tetrel bond donor entities that can also feature both σand π-holes, review by Zierkiewicz et al. [61].We believe that this study is poised to significantly aid researchers across diverse fields of natural science in identifying and subsequently characterizing π-hole halogen bonds in chemical systems, both theoretically and experimentally.

Figure 1 .
Figure 1.The [MP2/aug-cc-pVTZ] level potential on the electrostatic surface of (a,b) ClF3, (c,d) ClBr3, and (e,f) BrCl3 molecules, illustrating the anisotropic nature of the charge density.The 0.001 a.u.(electrons bohr −3 ) isoelectronic density envelope was used to map the potential.The location of σ-and π-holes on the surfaces of the first two molecules (a-d) are indicated by arrows in black.The regions of local most maxima and minima of potential (VS,max and VS,min) are marked by filled tiny circles in red and blue, respectively.Values are in kcal mol −1 .The exact number of minima and maxima of

Figure 1 .
Figure 1.The [MP2/aug-cc-pVTZ] level potential on the electrostatic surface of (a,b) ClF 3 , (c,d) ClBr 3, and (e,f) BrCl 3 molecules, illustrating the anisotropic nature of the charge density.The 0.001 a.u.(electrons bohr −3 ) isoelectronic density envelope was used to map the potential.The location of σand π-holes on the surfaces of the first two molecules (a-d) are indicated by arrows in black.The regions

Figure 5 .
Figure 5.The [MP2/aug-cc-pVTZ] level potential on the electrostatic surface of (a) CS, (b) PN molecules, illustrates the anisotropic nature of the charge density.The 0.001 a.u.(electrons bohr −3 ) isoelectronic density envelope was used to map the potential.The regions of local most maxima and minima of potential (VS,max and VS,min) are marked by filled tiny circles in red and blue, respectively.Values are in kcal mol −1 .Most negative and positive regions are colored as red and blue, resepectively.

2
ρb and Hb) are identical for each pair of Br(π-hole)•••F bcps in Cl3Br•••F2, Br(π-hole)•••Cl bcps in Cl3Br•••Cl2 and Br(π-hole)•••Br bcps in Cl3Br•••Br2, but not in other three.In all of these dimers, the bond paths are developed between the terminal atoms of the dihalogen molecule and the Br's π-hole, providing further indication that the attractive interaction in them is Coulombic in nature.The larger the negative end-portion of the halogen atom in the dihalogen molecule, the stronger the bend of that atom toward the Br's π-hole (as in Cl3Br•••FCl, Cl3Br•••FBr and Cl3Br•••ClBr).The secondary Br•••Br interactions between the Br-ends of the Br2 molecule and the Cl atoms of the BrCl3 molecule are developed in the

Figure 5 .
Figure 5.The [MP2/aug-cc-pVTZ] level potential on the electrostatic surface of (a) CS, (b) PN molecules, illustrates the anisotropic nature of the charge density.The 0.001 a.u.(electrons bohr −3 ) isoelectronic density envelope was used to map the potential.The regions of local most maxima and minima of potential (V S,max and V S,min ) are marked by filled tiny circles in red and blue, respectively.Values are in kcal mol −1 .Most negative and positive regions are colored as red and blue, resepectively.For Cl 3 Br•••FH (Figure 2 (16)), the secondary interaction is stronger than the primary interaction.The hydrogen bond is the primary, and the halogen bond is secondary (ρ b values 0.0242 vs. 0.019 a.u.).The importance of the secondary interactions towards the stability of this and other complexes is not only recognized by QTAIM (Figure3(2-15)) but also by IGM-δg inter 's isosurface topologies (Figure4(2-15)).NBO predicted that the strongest hyperconjugative interaction to occur from σ(F 2 ) to RY*(7)Br/RY*(8)Br with E(2)  of 0.51/0.25 kcal mol −1 for Cl 3 Br•••F 2 , σ(Cl 2 ) to RY*(6)Br with E(2) of 1.17 kcal mol −1 for Cl 3 Br•••Cl 2 , and σ(Br 2 ) to RY*(6)Br with E(2) of 0.88 kcal mol −1 for Cl 3 Br•••Br 2 .In the case of Cl 3 Br•••P 2 , the charge transfer delocalization corresponds to π(2)P 2 /π(3)P 2 → RY*(6)Br/RY*(5)Br with E(2) of 0.23/0.13kcal mol −1 and np(Cl) to σ*/π*(1/2/3) P 2 with E (2) < 1.8 kcal mol −1 ; the former is indicative of the presence of Br(π-hole)•••π(P 2 ) halogen bond.The F 2 , Cl 2 , and Br 2 molecules are symmetrically oriented on the top of the Br's π-hole while forming dimers with BrCl 3 (Figure2(3),(8), and (9), respectively).This is not the case for the dimers when ClF, BrF, and BrCl served as nucleophiles (Figure2(6),(7), and (10), respectively).Therefore, the topological properties of charge density (ρ b , ∇ 2 ρ b and H b ) are identical for each pair of Br(π-hole)•••F bcps in Cl 3 Br•••F 2 , Br(π-hole)•••Cl bcps in Cl 3 Br•••Cl 2 and Br(π-hole)•••Br bcps in Cl 3 Br•••Br 2 , but not in other three.In all of these dimers, the bond paths are developed between the terminal atoms of the dihalogen molecule and the Br's π-hole, providing further indication that the attractive interaction in them is Coulombic in nature.The larger the negative end-portion of the halogen atom in the dihalogen molecule, the stronger the bend of that atom toward the Br's π-hole (as in Cl 3 Br•••FCl, Cl 3 Br•••FBr and Cl 3 Br•••ClBr).The secondary Br•••Br interactions between the Br-ends of the Br 2 molecule and the Cl atoms of the BrCl 3 molecule are developed in the Cl 3 Br•••Br 2 dimer due to the larger size of the Br 2 molecule.They are not only weak but also the forced consequence of the primary interaction.

Figure 6 .
Figure 6.(a) Linear relationship between E[SPAT2+3(CCD)] and ∆E(BSSE) [MP2/aug-cc-pVTZ] for all the 33 binary complexes examined, described by the equation representing y.The insert marked by a black arrow in (a) corresponds to weakly to strongly bound complexes (∆E(BSSE) < −22.0 kcal mol −1 ).The square of the regression coefficient, R 2 , is shown in each case.Plot of the dependence of E[SPAT2+3(CCD)] on (b) the electrostatic energy component (E els ), (c) the exchange-repulsion energy component (E exch ), and (d) the induction energy component (E ind ) for the 33 complexes investigated (see Figure 2 for geometries andTable 1 for energies).
That they are potentially non-linear, vertical interactions, evidence of the C-I•••X (X = O, F, C) angles (∠C-I•••O/F values between 63 • and 115 • ).Our investigation shows that the halogen derivatives featuring π-holes in molecules can also form halogen bonds with charge density-rich regions hosted largely by O, F, C, and N atoms in neutral/anionicmolecular entities.density-rich regions hosted largely by O, F, C, and N atoms in neutral/anionicmolecular entities.
have labeled them as secondary I•••O bonding interactions developed between the λ 3 -iodine atom and the carbonyl oxygen of the acetyl group of the N-acetylbenziodazole skeleton.The I•••N type close contacts were also observed )iodo)-m-toluene [C13H11IN2O4S], and in crystal adducts such as [C10H11IN3O + ,C7H7O3S − ] with I•••N bond distances of 3.667 [144], and 3.587 Å [145], respectively.
have labeled them as secondary I•••O bonding interactions developed between the λ 3 -iodine atom and the carbonyl oxygen of the acetyl group of the N-acetylbenziodazole skeleton.The I•••N type close contacts were also observed in molecularly assembled crystals such as (((p-Nitrophenylsulfonyl)imino)iodo)-m-toluene [C 13 H 11 IN 2 O 4 S], and in crystal adducts such as [C 10 H 11 IN 3 O + ,C 7 H 7 O 3 S − ] with I•••N bond distances of 3.667 [144], and 3.587 Å [145], respectively.