A Continuum from Halogen Bonds to Covalent Bonds : Where Do l 3 Iodanes Fit ?

The intrinsic bonding nature of l3-iodanes was investigated to determine where its hypervalent bonds fit along the spectrum between halogen bonding and covalent bonding. Density functional theory with an augmented Dunning valence triple zeta basis set (wB97X-D/aug-cc-pVTZ) coupled with vibrational spectroscopy was utilized to study a diverse set of 34 hypervalent iodine compounds. This level of theory was rationalized by comparing computational and experimental data for a small set of closely-related and well-studied iodine molecules and by a comparison with CCSD(T)/aug-cc-pVTZ results for a subset of the investigated iodine compounds. Axial bonds in l3-iodanes fit between the three-center four-electron bond, as observed for the trihalide species IF 2 and the covalent FI molecule. The equatorial bonds in l 3-iodanes are of a covalent nature. We explored how the equatorial ligand and axial substituents affect the chemical properties of l3-iodanes by analyzing natural bond orbital charges, local vibrational modes, the covalent/electrostatic character, and the three-center four-electron bonding character. In summary, our results show for the first time that there is a smooth transition between halogen bonding ! 3c–4e bonding in trihalides ! 3c–4e bonding in hypervalent iodine compounds ! covalent bonding, opening a manifold of new avenues for the design of hypervalent iodine compounds with specific properties.


Introduction
Hypervalent iodine compounds (HVI) are important alternatives to transition metal reagents because of their reactivity, synthetic utility, low cost, abundance, and non-toxic nature [1][2][3][4][5][6].HVIs are involved in a multitude of reactions such as: reductive elimination, ligand exchange, oxidative addition, and ligand coupling [7,8].The three-center four-electron bonds (3c-4e) in HVI are weak and polarizable, which is valuable in synthetic organic chemistry, as they can exchange leaving groups or accept electrophilic/nucleophilic ligands depending on their surroundings [9].Despite such utility, there are still unknowns regarding the intrinsic bonding nature of HVIs and hypervalency in general.Though iodine is a halogen, it behaves like a metal; it is the heaviest non-radioactive element of the periodic table and is the most polarizable halogen [10,11].Because of its diffuse electron density (van der Waals (vdW) radius of ca 2 Å), iodine is a good electron donor, but can also serve as an electron acceptor [12,13].Iodine is not known to participate in d-orbital or p-interactions, though this could be further investigated [14,15].
HVIs commonly exist in the oxidation states 3, 5, and 7, which support 10, 12, and 14 valence electrons, respectively [16].Most common are the oxidation states 3 and 5, which are referred to as l 3 and l 5 -iodanes [17].l 3 -iodanes form distorted T-shaped molecular geometries, while l 5 -iodanes generally prefer square pyramidal geometries, as confirmed through both X-ray crystallography and computational studies [18,19].These somewhat unusual molecular geometries are the result of the pseudo-Jahn-Teller effect [20].The atoms that make up the "T" in l 3 -iodanes form improper dihedrals and non-ideal bond angles.The causes of these angular and dihedral deviations are unknown, but have been related to the anisotropic nature of the electronic density distribution in iodine [21][22][23][24][25]. Bader et al. showed that the vdW radius in iodine is larger at the equatorial position than at the axial position [26][27][28].This supports the observation that electronegative ligands favor the axial positions in iodine [29].
Hypervalency has been defined in several ways: Musher characterized main group elements in higher oxidation states as hypervalent [30].A successful concept employed to explain hypervalency without involving d-orbitals is the formation of multi-centered electron-deficient bonds [31][32][33].In this context, the 3c-4e bond model of Pimentel-Rundle is especially useful.According to this model, three atoms linearly align, each of which contributes an atomic orbital to form three molecular orbitals; a bonding orbital, a non-bonding orbital, and an antibonding orbital.Since only four electrons are involved, the antibonding orbital is unoccupied.As a result, two bonds share a single bonding electron pair (i.e., they have a fractional bond order of 0.5).The formation of two or more electron-deficient bonds allows hypervalent compounds to have higher oxidation states without necessarily expanding their octets.A direct consequence of this model is that the 3c-4e bond is expected to be substantially weaker than the two-center two-electron (2c-2e) bond in a given hypervalent molecule.Even though there are various works showing that d-orbital contributions to hypervalent bonds (HVB) are minimal, many chemistry text books still make use of the idea of an extended octet and the formation of spd-hybrid orbitals to explain HVB [34].There is a strong overlap between the concepts of fractional bond order, the 3c-4e bond, and the halogen bond (XB) in trihalides, which are considered prime examples of 3c-4e bonding, but also strong XB [35,36].A formal definition of XB is given in the following paragraph.
3c-4e HVI bonding (3c-4e HVIB) draws comparisons to the secondary bonding interaction due to the weak bond strength, high reactivity, and long internuclear distances exceeding covalent bond lengths [37].3c-4e HVIB also shares similarities with non-covalent interactions, along with hydrogen bonding [38][39][40], XB [35,[41][42][43], pnicogen bonding [43][44][45], chalcogen bonding [43,46], and tetrel bonding [47].XB is a non-covalent interaction between an electrophilic halogen (X) and a nucleophile with a lone pair (lp(A)) of donating electrons [42,43].For the remainder of this work, we will express lone pairs as (lp).The nucleophile/halogen acceptor (A), donates electrons to the antibonding (s*(XY)) orbital of the halogen donor (Y) [48,49].XB is also known to have an X-A distance that is shorter than the sum of the vdW radii with Y-X-A angles close to 180 degrees [35,41,50,51].Because of the obvious similarities between 3c-4e HVIB and XB, it has been argued that HVB should not be considered as a special bonding class [52,53].On the other hand, the term hypervalency has been widely accepted by the chemistry community, and therefore, its continuous use has been advocated [54].Based on this controversy, we decided to delve deeper into the bonding nature of l 3 -iodanes.
In this work, we investigated the intrinsic nature of HVIB in l 3 -iodanes and its relation to XB, 3c-4e bonding in trihalides, and covalent bonding to determine if there is a smooth transition between these interactions.Additionally, we studied the role equatorial ligands play in strengthening the 3c-4e bond in l 3 -iodanes, as well as substituent effects in the axial ligands.We utilized density functional theory (DFT), vibrational spectroscopy, quantum theory of atoms in molecules (QTAIM) combined with the Cremer-Kraka criterion for covalent bonding [55][56][57], and the natural bond orbital (NBO) analysis to characterize the nature and the intrinsic strength of HVIB.This investigation was rationalized by studying a diverse group of 34 HVI compounds shown in Figure 1, including known chemical compounds complemented by some model compounds.The remainder of this work is presented as follows: data, results, and discussion are presented in Section 2; Section 3 gives a description of the computational methods utilized; and Section 4 gives conclusions, the outlook, and future goals.

Results and Discussion
Figure 1 shows the 34 HVI compounds with selected NBO charges investigated in this work.Note that the abbreviation (PhI) is used to refer to iodobenzenes.The compounds are organized into four groups.Group 1 contains four reference compounds; the covalent complex FI (1.1), the 3c-4e trihalide IF 2 (1.2), fluorophenyl iodide (1.3), and the l 3 -iodane, IF 3 (1.4), which is used to study the effects of electronegative equatorial ligands.Group 2 and Group 3 compounds, iodobenzenes (2.1-2.10) and (3.1-3.11),respectively, represent l 3 -iodanes with a phenyl group equatorially bound to I.While both axial ligands in Group 2 molecules are halogen atoms, in Group 3, the axial ligands consist of a halogen and a non-halogen lone-pair-bearing functional group (CN, NH 2 , and OH).Group 4 consists of four halobenziodazoles (4.1-4.4) and four halobenziodoxoles (4.5-4.8).

NBO Charge Analysis
XB is a non-covalent interaction formed between a halogen donor molecule (YX) (e.g., a dihalogen, interhalogen, or halogenated molecule) and a halogen acceptor atom A, where A is an electron-rich atom; e.g., a nucleophilic heteroatom with lone-pair (lp) electrons [9,35,[41][42][43].The general charge transfer picture in XB describes a transfer of charge from lp (A) into the empty s ⇤ (YX) orbital [36,58,59].Applied to our set of HVI compounds, there are two possibilities to define the YI and the IA part of Group 2 molecules 2.5-2.10.In these cases, we chose IA to be the weaker of the two axial bonds, in analogy to XB.The same definition was also applied for Group 3 and Group 4 compounds.For all 34 HVI compounds investigated in this work, the central iodine atom I holds a positive charge (see purple numbers in Figure 1), ranging from +0.339 e in 1.3 to +1.803 e in 1.4.There is a significant difference in the central iodine charge when comparing the l 3 -iodanes with the three reference compounds 1.1-1.3 for which the iodine charge ranges from +0.339 e to +0.521 e, whereas iodine charges in l 3 -iodanes range from +0.779 e in 2.4 to +1.803 e in 1.4.The obvious difference between the l 3 HVIs and the non-HVIs is the equatorial ligand, which is absent in 1.1-1.3.In every case, except 1.4, the equatorial ligand is a phenyl (Ph) group (Groups 2-3) or a benzene ring (Group 4).Group charges for Ph/benzene ligands are negative in every case, as are charges on the Ph/benzene C atom bound to I, as shown in Figure 1.Increased positive charge on I indicates that the equatorial ligands are pulling charge from the central I.In 1.4, the equatorial ligand is F; this is an extreme case of a strong electron-withdrawing ligand in the equatorial position, which polarizes the central I atom.As a result, each I-F interaction involved in the 3c-4e bond becomes more polar.
2.1-2.4 comprise l 3 -iodanes of the type PhIA 2 , for A = F, Cl, Br, and I. Caused by the increasing electronegativity from I to F, the bond polarity increases from I-I < Br-I < Cl-I < F-I, with the positive charge on the central I atom increasing in the same order, i.e., +0.779 e in 2.4 to +1.466 e in 2.1.The same trend is observed in 2.5-2.7.As the second substituent changes from Cl to I, charge on the F ligand remains almost unchanged (Dq = 0.012 e), while the positive charge on I decreases.The same pattern continues for compounds 2.8-2.10.The charge on the equatorial ligand appears to be independent of the axial ligands for Group 2, with the exception of 2.1, 2.5-2.7 with F as the common axial ligand.The charge on the equatorial ligand in these cases tends to be most negative (ranging between 0.089 e and 0.120 e) compared to the other Group 2 members (ranging between 0.069 e and 0.079 e).
In Group 3, one of the axial halogen atoms is replaced with a functional group (NH 2 , OH, and CN).For comparison, we refer to OCO and OC(CH 3 ) in the halobenziodazoles 4.1-4.4 and halobenziodoxoles 4.5-4.8 of Group 4 as functional groups; see Figure 1.In 3.2-3.4,charge on NH 2 remains consistent.3.1 is the exception, but the charge difference between these four molecules is 0.034 e.This trend is observed in 3.5-4.8as well, where Dq (OH) = 0.008 e, Dq (OCO) = 0.016 e, Dq (OC(CH 3 ) 2 ) = 0.007 e, and Dq (CN) = 0.014 e.An important trend is revealed: although the charge at the central I atom is dependent on both axial ligands, charge on the functional group is insensitive to the halogen substituents at the opposite side, particularly for Cl, Br, and I.The behavior of the axial halogens in Group 3 is the same as described for Group 2. The functional group with the largest net negative charge is the OCO group, followed by OH, OC(CH 3 ) 2 , CN, and NH 2 , respectively.The data reveal yet another important trend: the more electronegative the axial ligands are, the more negative charge they collect and the more negative charge they impose on the equatorial ligand.Considering the charge on the Ph group for Group 3 and the benzene ring for Group 4, almost the same functional group trend emerges: there is an increase of negative charge from NH 2 < CN < OH < OC(CH 3 ) 2 < OCO, with OH and OC(CH 3 ) 2 groups being interconverted.That means, regarding the equatorial ligands, the same trends are observed in Group 2 molecules and in Groups 3-4 molecules.For all compounds with an axial F atom, charge on the equatorial ligand becomes more negative.The benzene rings of Group 4 molecules have a substantial negative charge ranging from 0.172 e to 0.207 e, however less than the Ph group of 1.3, being 0.538 e.

Bond Strength Order
Interatomic distances (r), r(r b ), H(r b ), k a , BSO n values, and local mode frequencies (w a ) are summarized in Tables 1 and 2. Figures 2-4 1, the F• • • I interactions in 3.9 and 4.1 are the stronger bonds.These are the only two cases in Groups 3-4 where the halogen• • • I interaction is the stronger interaction.Therefore, CN and OCO will be considered A and F will be considered Y in Table 1 for these two cases only.However, this convention is not used for the BSO n plots in Figures 3 and 5.In these two figures, the Y• • • IA interaction is the bond between the non-halogen and I, and YI• • • A is the bond between the halogen and I. Figure 2a is a BSO n plot for interactions in 1.1-2.1.Comparing bond strength in axial I-F interactions in FI, IF 3 , PhIF 2 , IF 2 , and PhIF reveals the hypothesized trend: FI > IF 3 (3c-4e HVIB) > PhIF 2 (3c-4e HVIB) > IF 2 (3c-4e) > PhIF (XB).As expected, the 2c-2e FI bonds in 1.1 and 1.4 are stronger than 3c-4e in 1.2 and 1.4.The latter are stronger than the I• • • F XB in 1.3, as shown in Table 1 and Figure 2. It is notable that the hypervalent IF 3 forms a shorter and stronger 2c-2e bond compared to the FI molecule.Axial F atoms pull charge from iodine, resulting in a more polar 2c-2e I-F bond, and also contract the I orbital, improving I-F orbital overlap.The 3c-4e bond in IF 3 is stronger than that of PhIF 2 because the equatorial F polarizes the central I, causing more polar and stronger interactions at the axial positions.The equatorial Ph and benzene groups do not have the polarizing ability of F, but they do bind strongly to (BSO n = 0.648-1.046)and pull charge from the central I.The effect of the equatorial ligand causes the difference in bond strength between PhIF 2 and IF 2 .IF 2 has no equatorial ligand and less polar bonds than the 3c-4e bond in PhIF 2 .When replacing a F with a Ph group (1.3), the negative charge becomes localized in the fluorine due to its higher electronegativity; as a result, two different types of bond are formed: one is the 2c-2e C-I bond, and the other is an XB between I and F (the lower polarizing power of Ph results in a less positive charge at the iodine).
Figure 2b shows Y• • • I• • • A (where A = F, Cl, Br, and I in this case) in 2.1-2.10.In the case of PhIY 2 (2.1-2.4),there is a correlation between bond strength and bond polarity.Charge on the central I atom increases in the series: PhI 3 < PhIBr 2 < PhICl 2 < PhIF 2 .This matches the trend in 3c-4e bond strength.Charge on the axial ligand also matches this order, but with charge becoming more negative.In 2.5-2.10,there is a marked difference in I• • • F bond strength (BSO n > 0.562) and all other axial bonds (BSO n = 0.272-0.423).I• • • Cl, I• • • Br, and I• • • I interactions are similar in bond strength, but vary slightly depending on the atom on the opposite side of the 3c-4e bond.This result is in accord with observed bond polarity and electronegativity trends of halogens.
Figure 3a shows BSO n plots for Y The same trend emerges again when replacing one halogen with an electron-donating functional group; the bond strength of I• • • A increases when A changes from I to F. Keeping the Y constant and substituting A again reproduce the trend that bonds become stronger when going up the periodic table from I to F for all five functional groups.When comparing the functional groups, bond strength follows this order: OCO > OC(CH 3 ) 2 > OH > CN > NH 2 .This order holds regardless of the axial halogen.This pattern nearly matches the order observed in group charges where the more negatively-charged the group, the higher the BSO n.The exception is OH and OC(CH 3 ) 2 .OH groups are more negatively charged, but do not bind as strongly as OC(CH 3 ) 2 groups.This is because the benzene in 4.5-4.8binds more strongly on average in Group 4 (BSO n = 0.898-1.046)than Ph in Group 2 (BSO n = 0.691-0.914)and Ph in Group 3 (BSO n = 0.648-0.893).The stronger equatorial bond correlates to a more positive charge on the central I, which allows for stronger 3c-4e bonds.The key difference between Group 4 and Groups 2-3 is that all of Group 4 has functional groups bound directly to benzene and Groups 2-3 do not.In this case, it is justifiable to state that the C(CH 3 ) 2 group in 4.5-4.8 will be an electron donor to benzene, which accounts for the lower group charge.2.1 and 3.8 are the only exceptions to this trend.5a), and YI• • • A in Groups 3-4 (Figure 5b).The vertical dashed line through the origin separates the covalent region from the electrostatic region according to the Cremer-Kraka criterion.H(r b ) < 0 for every I-A, I-Y, and I-C equatorial (Ph, F, or benzene) interaction, putting them in the covalent bonding region or very close to the electrostatic region in some cases.There is significant covalent contribution for the axial bonding interactions in 2.1-4.8,indicating that charge accumulation in the bonding region produces a net stabilizing effect.For the plot in Figure 4b, there is a good linear correlation between BSO n and H(r b ), as indicated by a value of R 2 = 0.930.These data correlate higher bond strength to an increase in covalent character of the interaction.The weaker the bond, the closer to the electrostatic region.The plot is sectioned off into regions to show agreement with Figure 2a.The XB interaction in PhIF is at the bottom of the plot, closest to the electrostatic region.The 3c-4e region is next, where IF 2 is found, along with all of the weakly-bound l 3 -iodanes containing axial Cl, Br, and I atoms.It is necessary to note that we expect the 3c-4e bond in IF 2 to be the strongest of all trihalides and to be at the very top of the spectrum.Therefore, if considering other trihalide systems, one would expect to see better separation between 3c-4e bonds in HVI and 3c-4e bonds in trihalides, as is observed for the more closely-related IF 2 , PhIF 2 , and IF 3 .All of the l 3 -iodanes containing F are at the very top of this region bordering the next region or in the next region, which is 3c-4e HVI.IF 3 and PhIF 2 give prime examples of the 3c-4e HVIB.At the very top right corner is the covalent F-I complex.As we follow the linear data from the weak electrostatic region to the strong covalent region, we once again reproduce the smooth continuum: partially-covalent XB < 3c-4e bond in trihalides < 3c-4e bond in HVI < covalent bond.Now, the trend holds in terms of covalent/electrostatic character and H(r b ).Note that in Figures 4 and 5, r(r b ) could be plotted against BSO n in place of H(r b ), and the same correlation would occur, but with a positive slope instead of a negative one.

Covalent/Electrostatic Contributions
The same general trend is observed in Figure 5 for Groups 3-4.As BSO n increases, H(r b ) becomes more negative (deeper into the covalent region).In Figure 5a, points for Y• • • I are scattered, and the correlation weakens when taking the data as a whole.However, if considering each functional group individually, a strong linear correlation once again occurs.The periodic trend emerges that as A, the halogen homolog becomes smaller, the bond strengthens and becomes more covalent in nature.This is not a direct result of the axial ligand, rather it is the result of the polarizing effect the axial halogen has on the central I atom.The Y• • • IA interactions (H(r b ) < -0.237 Hartree/Å 3 ) sit significantly farther into the covalent region compared to the YI• • • A interactions (H(r b ) < -0.055 Hartree/Å 3 ).Figure 5b again shows a reasonable linear correlation with R 2 = 0.917.The trend amongst functional groups previously noted in 3.1-3.2 is once again evident here: OCO > OC(CH 3 ) 2 > OH > CN > NH 2 in terms of pulling charge from the central I, which results in strengthening and to some degree increasing H(r b ) of the I• • • A bond.Another important point here is that one must not assume certain functional groups will behave the same way in all situations as they behave when bound to benzene.A prime example is CN: a strong electron withdrawing group (when bound to benzene) is the second weakest withdrawing group in this study.

Computational Methods
DFT was utilized to optimize molecular geometries and to calculate for each stationary point molecular vibrational frequencies including the L-modes of Konkoli and Cremer [60][61][62] and the determination of local mode force constants (k a ), NBO charges, electron densities r(r b ), and energy densities H(r b ); where r b is a bond critical point.Each stationary point was confirmed as a minimum by absence of imaginary normal mode frequencies.Available experimental geometries for the ICl 3 dimer, IF 3 , IF 5 , dichloroiodobenzene (PhICl 2 ), and diacetoxyiodobenzene (PhI(OAc) 2 ) [18,[63][64][65][66] were used to gauge the accuracy of the DFT calculations.Experimental and calculated geometries using different model chemistries for this set of compounds are compared in Tables A1 and A2 (See Appendix A).We initially employed Grimme's Rung 5 double hybrid density functional B2PLYP [67] and Dunning's cc-pVDZ basis set [68][69][70][71] with a tight convergence criterion and an ultra-fine integration grid.The B2PLYP functional combines the generalized gradient approximation exchange functional of Becke [72,73] and the Lee-Yang-Parr correlation functional [74] with exact Hartree-Fock exchange and Møller-Plesset perturbation theory [75][76][77][78] of second order (MP2) [79][80][81].This functional has shown close agreement between calculated and experimental geometries and vibrational frequencies for heavy atoms [82,83].However, for our set of molecules, the cc-pVDZ basis set did not produce the desired accuracy (Tables A1 and A2), and the B2PLYP/aug-cc-pVTZ level of theory became computationally expensive.The combination of MP2 and a relatively small double-zeta basis set is known to provide a fortuitous cancellation of error [84,85].MP2 overestimates correlation energy, but this is compensated by the cc-pVDZ basis set [86].Therefore, we tested MP2/cc-pVDZ for reducing the computer time.However, results calculated at this level of theory gave less accurate results than wB97X-D/aug-cc-pVTZ [87,88], while calculations at the B2PLYP/Def2TZP level of theory led to inaccurate results in several cases.For Br and I, scalar relativistic effects were assessed by using effective core potentials (ECPs) in combination with the Dunning basis sets [89,90].
Although geometries are first order properties and therefore less sensitive to the level of theory, B2PLYP/aug-cc-pVTZ and wB97X-D/aug-cc-pVTZ calculations turned out to be in closest agreement with experimental data, while for a small subset of compounds, close agreement between wB97X-D/aug-cc-pVTZ and CCSD(T)/aug-cc-pVTZ was obtained (Table A3).To further rationalize these results, gauge-independant atomic orbital (GIAO) magnetic shielding tensors [91][92][93][94][95] were calculated and isotropic shielding constants were converted into chemical shifts utilizing the linear regression method of Tantillo et al. for PhICl 2 , PhI(OAc) 2 , and 1-Hydroxy-1,2-benziodoxol-3(1H)-one [96][97][98][99][100][101].This method requires the calculation of isotropic magnetic shielding tensors for a test set of molecules at a given level of theory (in our case, wB97X-D/aug-cc-pVTZ and B3LYP/aug-cc-pVTZ), plotting the raw calculated isotropic value against experimental NMR chemical shifts, and using the following relationship to develop an equation for calculating chemical shifts (Figure A1): where d is the derived chemical shift and s is the calculated isotropic magnetic shielding tensor.The margin of error for proton-NMR chemical shifts turned out to be 0.24-6.91%for the B3LYP functional and 0.19-5.81%for the wB97X-D functional (Table A4) [18,63,102].Although both wB97X-D and B3LYP gave satisfactory and similar calculated chemical shifts, wB97X-D gave more accurate geometries and frequencies.
Based on these findings, the wB97X-D/aug-cc-pVTZ level of theory was chosen for this study due to its displayed ability to predict accurate first and second order experimental properties in HVI molecules in addition to the previous findings of Oliveira et al. that this level of theory is suitable for the detailed analysis of XB [41].
Vibrational spectroscopy was applied to quantify the intrinsic strength of HVIBs.Chemists have utilized vibrational spectroscopy to obtain information about the electronic structure of molecules and their framework of bonds.However, normal vibrational modes cannot be used as a direct bond strength measure because they are delocalized due to electronic and mass coupling, a fact that often has been overlooked [103,104].The electronic coupling is eliminated by solving the Wilson equation of spectroscopy [105] and transforming to normal coordinates.Konkoli and Cremer showed that the remaining mass (kinematic) coupling can be eliminated by solving a mass-decoupled equivalent of the Wilson equation, leading to local vibrational modes, which are associated with internal coordinates q n such as bond lengths, bond angles, and dihedral angles [60,106].Zou and Cremer verified that there is a one-to-one relationship between local and normal vibrational modes through an adiabatic connection scheme (ACS) [107][108][109], allowing a normal mode decomposition into local mode contributions [44,110,111] and, as such, the detailed analysis of a vibrational spectrum.This is of particular value, given the fact that L-modes can be applied to both calculated and measured spectra [61,112].
Another important feature of L-modes is the direct relationship between the local stretching force constant (k a ) of a chemical bond and its intrinsic strength [113].This has enhanced our knowledge about chemical bonding and the often overlooked, but highly important weak intermolecular interactions, providing a wealth of new insight into: (i) covalent bonding [113], stretching from peculiar cases of reversed bond length-bond strength relationships [114,115], to a new design recipe for fluorinating agents [116]; (ii) weak chemical interactions including hydrogen bonding [117,118], XB [35,41,42], pnicogen bonding [43], chalcogen bonding [50], weak interactions in gold clusters [119], as well as non-classical hydrogen bonds in boron-hydrogen • • • p interactions [120,121].In addition, new electronic parameters and rules were derived [122][123][124].
When comparing a larger set of k a , the use of a relative bond strength order (BSO n) is convenient [103,104].The BSO n of a bond is obtained by utilizing the extended Badger rule [103,104,125] according to which BSO n is related to k a by a power relationship, which is fully determined by two reference values and the requirement that for a zero-force constant, the BSO n value becomes zero: The constants a and b are calculated from k a values of two reference compounds with known BSO n values n 1 and n 2 via: and, b = ln(n 2 /n 1 )/ln(k a 2 /k a 1 ) In this work, we chose as reference compounds FI and IF 2 representing BSO n values of one and 0.5, respectively, guided by the corresponding Mayer bond orders [126] of 0.940 and 0.543 evaluated at the wB97X-D/cc-pVTZ level of theory.More than 50% of iodine bonds in this work include an atom from the second period.This renders the FI/IF 2 reference system ideal (a second period atom bound to iodine), in addition to providing a spectrum with a full 3c-4e bond from a trihalide (IF 2 ) on the one end and a full covalent bond (FI) on the other end.
Using the k a 1 of 1.913 mdyn/Å for IF 2 and 3.953 mdyn/Å for FI (wB97X-D/cc-pVTZ level of theory), the constants a and b in the power relationship Equation (3) were determined to be a = 0.269 and b = 0.955, leading to: BSO n = 0.269(k a ) 0.955 (6) Because the chosen reference system was designed for 3c-4e interactions particular to this study, a scaling procedure was used to obtain appropriate BSO n values for covalent I-C interactions between the equatorial ligands and the central iodine.The equatorial bonds are fully-formed single bonds, but the C-I bond is much less polar and weaker than the I-F bond used as a reference.We calculated k a = 2.557 mdyn/Å for the I-C bond in iodobenzene.From Equation ( 6), we calculated BSO n = 0.659.The scaling factor was obtained by setting n = 1 for this I-C bond.The scaling factor is 1/0.659= 1.517, which was applied to BSO n of all equatorial I-C bonds.Multiplying the scaling factor through Equation ( 6) provided a new BSO n equation for assessing the strength of the equatorial I-C bonds in this study: BSO n(scaled) = 0.408(k a ) 0.955 The Cremer-Kraka criterion was applied to assess the covalent nature of HVIB [42,55,56,77,127].According to this criterion, a covalent bond between two atoms A and B is defined by (1) the existence of a zero-flux surface and bond critical point (r b ) between atoms A and B (necessary condition) and (2) a negative and thereby stabilizing local energy density H(r b ) (sufficient condition).H(r b ) will be close to zero or positive if the interaction between A and B is non-covalent, that is electrostatic or of the dispersion type.H(r) is defined as: where G(r) is the kinetic energy density (always positive, destabilizing) and V(r) is the potential energy density (always negative, stabilizing).In addition to the established Cremer-Kraka criterion, a molecular fragmentation scheme for estimating electron density shifts has recently emerged as a potential tool for the qualitative investigation of non-covalent interactions at low computational cost [128].L-modes was carried out with the program COLOGNE2018 [129], and Mayer bond orders were determined with the program ORCA [130].NBO populations were computed using NBO 6 [131][132][133][134].The electron density analysis, in particular the calculation of electron density at the bond critical point (r(r b )) and H(r b ), was performed with the program AIMAll [135,136].All DFT calculations were carried out with GAUSSIAN16 [137].

Conclusions
In this work, we quantified the intrinsic bond strength and bonding nature of a series of HVI compounds through vibrational spectroscopy.Use of DFT in this work was rationalized by testing several levels of theory against first and second order experimental properties of a small set of known HVI reagents.The computed set of 34 HVI molecules was then compared to XB, 3c-4e bonding, and covalent bonding in terms of BSO n, k a , r(r b ), H(r b ), and NBO charges.Recently, Politzer and coworkers [138] showed that by substituting a ligand in trihalides with a negative point charge, the positive electrostatic potential at the polarized s-hole collinear to the point charge correlates qualitatively well with the interaction energy; substantiating the key role played by electrostatics, which is also reflected in the atomic charge distribution (see Figure 1) and can be rationalized in terms of the 3c-4e model.The more negative charge at the ligands Y and A compared to the central iodine is due to the presence of a node at the center of the occupied non-bonding orbital [139].This charge separation is responsible for the lower covalent character of 3c-4e bonds compared to a classical 2c-2e bond.Politzer and coworkers proposed the existence of a continuum between non-covalent and covalent bonds, the latter being a result of an increased degree of polarization [138].Our results do also suggest the existence of such a continuum, but whether covalency can be seen as a degree of polarization is still disputable, especially in view of Ruedenberg's description of covalent bonding, where energy lowering is a result of the complex interplay of kinetic and potential energy contributions [140,141].The 3c-4e bonds in HVI share properties with XB, but are more closely related to the 3c-4e bonds in trihalides or covalent bonding in extreme cases.The equatorial 2c-2e HVI bond is stronger than comparable 3c-4e bonds (bonds involving the same ligands like in IF 3 ) and is more closely related to a covalent bond.Our results support the following transition: XB < 3c-4e bond in trihalides < 3c-4e bond in HVI < 2c-2e bond in PhIF 2 < covalent bond.When comparing the difference (equatorial ligands) between trihalides and l 3 -iodanes, we found that the 3c-4e HVIB is strengthened by the equatorial ligand by comparing IF 2 , PhIF 2 , and IF 3 .The equatorial ligand contributes significantly in pulling electron density from the central I, allowing for more polar interactions.Thus, highly electronegative ligands at the equatorial position will form strong interactions, as will axial ligands in such a case.We also found that axial ligands in HVIs have a minimal direct effect on one another in terms of NBO charge analysis, but do play a role in altering charge on the central I. Substituent effects in HVI can alter bond strength in both axial ligands and the equatorial ligand, particularly when F atoms are involved as ligands.The five functional groups studied here play a bond-strengthening and -polarizing role in the following order: OCO > OC(CH 3 ) 2 > OH > CN > NH 2 , with OH and OC(CH 3 ) 2 being partially interchangeable.In terms of H(r b ), we found a strong linear correlation with BSO n.H(r b ) becoming more negative correlates to an increase in bond strength.Furthermore, large V(r) stabilization in the bonding region correlates to the increased covalent character of a bond.Finally, we found the 3c-4e bond concept to be a valuable descriptor in terms of the linear portion of l 3 -iodanes.
Future goals are to utilize L-modes and the analysis of the electrostatic potential to explain why the T-shaped molecular geometry in l 3 -iodanes contains improper dihedrals and non-ideal bond angles.We also plan to investigate 3c-4e bonding and intramolecular HB in a series of HVI reagents utilizing L-modes and to explore the chemical reactivity of HVI compounds utilizing the unified reaction valley approach developed in our group [103,[142][143][144].In addition, we will perform a conformational and geometrical study of a series of novel HVI monomeric materials with a strong potential of forming useful polymers [145,146].B2PLYP/aug-cc-pVTZ and wB97X-D/aug-cc-pVTZ levels of theory give superior results for the given geometry parameters.Table A2 shows bond lengths, bond angles, and % error compared to experimentally measured data for PhICl 2 and PhI(OAc) 2 computed at all of thee aforementioned levels of theory.These two molecules are similar, or the same in the case of PhICl 2 as the majority of the molecules in this work.The wB97X-D/aug-cc-pVTZ gave remarkable accuracy in calculating geometry parameters for these two molecules.at the wB97X-D/CCSD(T)/aug-cc-pVTZ level of theory.Once again, wB97X-D/aug-cc-pVTZ performs remarkably well compared to the gold standard CCSD(T).Table A4 shows calculated and experimental NMR shifts for PhICl 2 , PhI(OAc) 2 , and 1-hydroxy-1,2-benziodoxol-3(1H)-one using the B3LYP and w97X-D functionals with the aug-cc-pVTZ basis set. Figure A1 shows a strong linear correlation between calculated isotropic magnetic stretching tensors and experimentally measured chemical shifts.The calculations done at the wB97X-D/aug-cc-pVTZ level of theory are slightly more in agreement with experimental measurements than the B3LYP/aug-cc-pVTZ.
show the power relationship between BSO n and k a for Group 1, Y• • • IA and YI• • • A in Groups 2-4, where YI• • • A is weaker than Y• • • IA, and I-C equatorial (Ph, F, benzene) bonds in Groups 1-4.Note: in Table

Figure 3 .
Figure 3. (a) BSO n versus for k a of Y• • • IA (3.1-4.8)(b) and BSO n versus k a of YI• • • A for complexes 3.1-4.8according to eq 2, where A = the axial halogen atom and Y = the non-halogen atom bound axial to I.

Figure 4 .
Figure 4. (a) BSO n versus (k a ) of I-Ph (2.1-3.12) and I-benzene (Group 4) according to eq 2. BSO n for the equatorial I-C bonds has been scaled by a factor of 1.517.(b) Comparison of BSO n with H(r b ) of all halogen-iodine interactions in complexes 1.1-2.10.The vertical dashed line separates the electrostatic region from the covalent region.

Figures 4b and 5
Figures 4b and 5 contain three plots correlating BSO n with H(r b ) of I• • • axial halogens in Group 2 (Figure 4b), Y• • • IA in Groups 3-4 (Figure5a), and YI• • • A in Groups 3-4 (Figure5b).The vertical dashed line through the origin separates the covalent region from the electrostatic region according to the Cremer-Kraka criterion.H(r b ) < 0 for every I-A, I-Y, and I-C equatorial (Ph, F, or benzene) interaction, putting them in the covalent bonding region or very close to the electrostatic region in some cases.There is significant covalent contribution for the axial bonding interactions in 2.1-4.8,indicating that charge accumulation in the bonding region produces a net stabilizing effect.For the plot in Figure4b, there is a good linear correlation between BSO n and H(r b ), as indicated by a value of R 2 = 0.930.These data correlate higher bond strength to an increase in covalent character of the interaction.The weaker the bond, the closer to the electrostatic region.The plot is sectioned off into regions to show agreement with Figure2a.The XB interaction in PhIF is at the bottom of the plot, closest to the electrostatic region.The 3c-4e region is next, where IF 2 is found, along with all of the weakly-bound l 3 -iodanes containing axial Cl, Br, and I atoms.It is necessary to note that we expect the 3c-4e bond in IF 2 to be the strongest of all trihalides and to be at the very top of the spectrum.Therefore, if considering other trihalide systems, one would expect to see better separation between 3c-4e bonds in HVI and 3c-4e bonds in trihalides, as is observed for the more closely-related IF 2 ,

Figure 5 . 1 - 4 . 8 .
Figure 5. (a) Comparison of BSO n with H(r b ) of iodine-non halogen axial atom (Y• • • IA) interactions in complexes 3.1-4.8and (b) a comparison of BSO n with H(r b ) of axial halogen-iodine (YI• • • A) interactions in complexes 3.1-4.8.The vertical dashed line separates the electrostatic region from the covalent region.

Figure A1 .
Figure A1.(a) Computed at the B3LYP/aug-cc-pVTZ level of theory, isotropic magnetic shielding tensors plotted against experimental NMR chemical shifts showing strong linear correlation and a slope close to 1 which is indicative of minimization of systematic error.(b) Computed at the wB97X-D/aug-cc-pVTZ level of theory, isotropic magnetic shielding tensors plotted against experimental NMR chemical shifts again showing strong linear correlation and slope close to 1.
Schematic of the 34 molecules investigated showing the numbering system (group.molecule,e.g., 1.1-4.8)given in bold face, and natural bond orbital (NBO) charges calculated at the wB97X-D/aug-cc-pVTZ level of theory.Note: charges in grey, blue, red, and orange are Ph/benzene, NH 2 , OH, OCO (red), OC(CH 3 ) 2 (red), and CN group charges, respectively, and not atomic charges.
Figure 2. (a) Power relationship between bond strength order (BSO) n and k a of Group 1 bonds (axial bonds only in the case of PhIF , PhIF 2 , and IF 3 ) according to eq 2 and (b) BSO n versus k a of Y• • • I• • • A bonds for complexes 2.1-2.10,where Y and A are halogen atoms.

Table 2 .
r, r(r b ), H(r b ), and local vibrational data for I-equatorial ligand (I-F, I-C) interactions of complexes 1.4-4.8,where BSO n is scaled by a factor of 1.517.
(Y• • • I) > BSO n (A• • • I).A• • • I isthe weaker, less covalent bond, and Y• • • I is the stronger, more covalent bond.In 2.1-2.4,A=Y; therefore, 3c-4e is 100%.In 2.5-2.7,%3c-4edecreasesfrom 68% in 2.5 to 57% in 2.7 as DBSO n becomes larger.2.8-2.10 contain weakly-bound halogens, which promote high 3c-4e bonding character (88-97%).In 3.1-3.12,there is a large range of 3c-4e contributions to the Y• • • I• • • A interactions (45-94%).The highest percentage is in 3.5, where the 3c-4e interaction is HO• • • I• • • F. Both substituents have lp electrons and are highly electronegative.The bonds formed are strong and polar, as the central I is the most polarized of all Group 3 molecules with an NBO charge of +1.401 e. OH and F have similar BSO n and NBO charges: n = 0.625, 0.565 e and n = 0.590, 0.703 e, respectively.The 3.1 has high 3c-4e character for the same reason as 3.5, but with NH 2 involved instead of OH.N is slightly less electronegative than O, and NH 2 has a more positive charge than OH, and thus forms a slightly weaker, less polar bond.3.3-3.4and 3.7-3.8have the lowest 3c-4e character in Group 3.These species contain mostly I-Br or I-I bonds, which bind weakly, while on the other side of the Y• • • I• • • A, we have polar functional groups OH and NH 2 .There is a strong polar interaction on one side of I and a weak non-polar interaction on the other side, which decreases the 3c-4e character.4.1-4.4have high 3c-4e character (81-85%).The I-O oxygen is part of an ester group which carries a large negative charge and contributes resonance stabilization.In 4.5-4.8,I is bound to the O on a T-butoxy group, which is slightly less electron rich and does not have the benefit of resonance.The T-butoxy-O binds strongly to I compared to Cl, Br, and I.

Table A3 .
Calculated r and k a of all FI bonds in FI, IF 2 , IF3, and IF5 computed at the wB97X-D/CCSD(T)/aug-cc-pVTZ level of theory.