Next Article in Journal
Catalysis by Hybrid Nanomaterials
Next Article in Special Issue
A Valence-Bond-Based Multiconfigurational Density Functional Theory: The λ-DFVB Method Revisited
Previous Article in Journal
Valorization of a By-Product from the Production of Mechanically Deboned Chicken Meat for Preparation of Gelatins
Previous Article in Special Issue
Metal–Metal Bond in the Light of Pauling’s Rules
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 26
Issue 2
10.3390/molecules26020350
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of Bifurcated Halogen with Hydrogen Bonds

by
Steve Scheiner
Department of Chemistry and Biochemistry, Utah State University, Logan, UT 84322-0300, USA
Molecules 2021, 26(2), 350; https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26020350
Submission received: 28 December 2020 / Revised: 8 January 2021 / Accepted: 9 January 2021 / Published: 12 January 2021
(This article belongs to the Special Issue Chemical Bonding: A Commemorative Special Issue Honoring Professor Linus Pauling)
Browse Figures
Versions Notes

Abstract

:
Bifurcated halogen bonds are constructed with FBr and FI as Lewis acids, paired with NH3 and NCH bases. The first type considered places two bases together with a single acid, while the reverse case of two acids sharing a single base constitutes the second type. These bifurcated systems are compared with the analogous H-bonds wherein FH serves as the acid. In most cases, a bifurcated system is energetically inferior to a single linear bond. There is a larger energetic cost to forcing the single σ-hole of an acid to interact with a pair of bases, than the other way around where two acids engage with the lone pair of a single base. In comparison to FBr and FI, the H-bonding FH acid is better able to participate in a bifurcated sharing with two bases. This behavior is traced to the properties of the monomers, in particular the specific shape of the molecular electrostatic potential, the anisotropy of the orbitals of the acid and base that interact directly with one another, and the angular extent of the total electron density of the two molecules.
Keywords:
cooperativity; σ-hole; AIM; NBO; IR; NMR

Graphical Abstract

1. Introduction

Linus Pauling’s early explanation [1] of what has come to be called the hydrogen bond (HB) ushered in a long period of its investigation, including follow-up work by himself and coworkers [2,3,4]. Studies on this topic continue to this day, and are summarized periodically by books dedicated to descriptions of developments up to that point in time [5,6,7,8,9,10,11,12,13,14,15,16]. One of the chief advances of recent years has been an expansion of the list of eligible atoms that participate in such bonds [17]. While the early thinking focused on highly electronegative atoms like N, O, and F, this field has been greatly generalized to much less electronegative atoms like C, P, Se, and even metals [18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]. Along with this expansion in terms of atoms, has also come a broader concept of the source of electrons from the nucleophile that extends well beyond a lone pair, to π-systems, σ-bonds, and the half-filled orbitals of radicals [33,34,35,36,37,38,39].
Another extension of the HB concept arises in connection with the number of groups participating. A classical HB places the bridging proton squarely between a proton-donor and acceptor atom, as might occur for example between the two O atoms of the water dimer. While such a geometry might in fact represent a preferred structure, there have been numerous observations of what have come to be termed bifurcated HBs. What is usually meant by this term is the placement of the bridging proton between two different acceptors, D1 and D2 as depicted in Scheme 1A. Neither of the two alignments in such a proton-shared configuration are linear, and the position of the H is not necessarily precisely symmetric with respect to D1 and D2. Another sort of bifurcated HB shares a common electron donor between two different proton donors, as indicated in Scheme 1B.
Bifurcated H-bonds of either type are quite common as surveys of crystal structures have revealed [40]. In general, these bifurcated H-bonds have served as fertile ground for study over the years [41,42,43,44,45,46]. The 2-fluoroethanol trimer [47] furnishes one such particular example, where a single H atom can engage with O and N atoms simultaneously [48]. An example drawn from biology [49] occurs within the α-helix N-Caps of ankyrin repeat proteins. Bifurcated H-bonds are also necessary to stabilize fibrils of poly(L-glutamic) acid [50]. Another example [51] places a halide between two H atoms of an alkenic=CH2 group.
Due to their similarity to H-bonds it is not surprising that bifurcated halogen bonds are a distinct possibility, and they have undergone some scrutiny as well [52,53,54]. These sorts of bonds are equivalent to those depicted in Scheme 1, but with the H atoms replaced by a halogen atom. There is recent evidence of a bifurcated halogen bond [55] with a C-Br group interacting with both Cl and Pt atoms. In another example, an aryl I atom participated [56] in a bifurcated halogen bond with a pair of O electron donors of a PtO4 system. An I atom of a PtI2Br2 moiety [57] can involve itself in several halogen bonds simultaneously. A halogen atom placed between two N atoms in a bidentate diazaheterocyclic compound [58] provides another example. Even the first-row F halogen appears capable of participating in a bifurcated halogen bond [59] under certain conditions.
Like halogen atoms, chalcogen atoms are also capable of engaging in bifurcated interactions, as for example [60] wherein a single chalcogen atom binds symmetrically to a pair of electron donor atoms. Another example [61] pairs one of the chalcogen (Y) atoms of a Y=C=Y molecule with the two O atoms of 1,2-dihydroxybenzene. Chalcogen bonds can mix with H-bonds in a hybrid sort of bifurcation as found recently in a bifurcated supramolecular synthon in ebselen [62]. Pnicogen atoms are also capable of engaging in bifurcated interactions [63], as can triel atoms, for example when two triels interact with a common halide [64]. The same is true of rare gas atoms [65] which can accept electrons from a symmetrically disposed pair of O atoms.
A great deal has been learned about the halogen bond (XB) in recent years, and its parallels to the HB. Both sorts of bonds involve a strong electrostatic attraction between the electron donor and acceptor groups. While the bridging proton takes on an overall positive partial charge within its A-H subunit, the opposite is true of the X atom. However the electrostatic attraction to the nucleophile is made possible by a small positive region lying along the extension of the A-X bond, commonly referred to as a σ-hole. Both sorts of bonds are stabilized further by charge transfer from the nucleophile into the σ *(A-H/X) antibonding orbital. Also in common, electron-withdrawing substituents on the Lewis acid strengthen either bond by enhancing the positive charge residing on the H or X atom.
While much is now understood about the fundamental properties of standard halogen bonds, there is far less information about the corresponding bifurcated XBs. Under what conditions might they occur, and which or both of types A and B in Scheme 1 are possible? Is there a particular strength or type of base that is required in order for a bifurcated XB to occur? How does the strength of a bifurcated XB compare with that of a single linear bond, and how does bifurcation affect the lengths of the bonds involved? How close together can the pair of bases in Scheme 1A approach one another, and likewise for the two acids in Scheme 1B. Tying all these questions together is a comparison of the properties of bifurcated halogen bonds with the analogous H-bonds.
In order to answer these questions, a number of systems are devised here and evaluated via quantum chemical calculations. FBr and FI are both capable of forming halogen bonds by virtue of fairly deep σ-holes on the Br and I atoms, respectively. The larger size of the I atom imbues it with the ability to form somewhat stronger XBs of the two. FH is employed as the corresponding H-bonding acid for purposes of comparison. Two different bases are considered. NH3 is a stronger base than NCH, in part due to its sp3 hybridization in comparison to sp for the latter. Each of the three Lewis acids is paired with both NH3 and NCH, and for each case both the acid-shared Scheme 1A type and the base-shared type 1B are evaluated.

2. Results

In the case of simple dimers, whether halogen or H-bonded, the geometric preference is for a linear arrangement, as delineated in Figure 1a,b for sample systems FI··NH3 and FH··NCH, respectively. There are two types of bifurcated bonds that might be imagined. In the first place, two separate bases can compete for the same σ-hole on the halogen atom (or the H atom in the case of H-bonds) as delineated in Scheme 1A. An alternate scheme in Scheme 1B would place two Lewis acids in competition for a single lone pair on a base molecule.

2.1. Shared Lone Pair of Single Base

The systems representing Scheme 1B are illustrated in Figure 2 for both NH3 and NCH as bases, and each combined with a pair of FBr, FI, or FH molecules. None of the triads pictured there represent a true minimum. When two acid units were initially placed in positions similar to those in the figure, a full optimization led to large-scale motion of one acid, leaving a single linear dyad as in Figure 1. The remaining acid molecule would completely reposition itself so as to engage in another interaction altogether. For example, one of the two FI molecules would engage in a I··I halogen bond when paired with NH3. In the case of the NCH base, the second FI molecule maintained a I··N XB with N of the base, but a highly distorted one. The geometries illustrated in Figure 2 were the result of imposing a restriction that the two acid molecules are symmetrically placed around the base. In the case of 2 FI molecules combined with NCH, for example, both θ(CN··I) angles were forced to be equal to one another. In other words, each of the bifurcated bonding configurations of Figure 2 is less stable than one that contains a single, strong, and linear XB or HB.
A quantitative assessment of the bonding in these triads may be gleaned from Table 1. For each combination of acid and base, the first row provides the parameters for the pure, fully optimized and linear dyads. As an example, the first row shows that the R(Br··N) halogen bond length in FBr··NH3 is equal to 2.339 Å. The next row indicates this bond is stretched to 2.760 Å in the triad when a second FBr must share the N lone pair on NH3, as in Figure 2a. Note that the triad is not fully symmetric in that the two R(Br··N) distances are not quite equal, so it is the shorter of the two that is catalogued in Table 2. The two Br atoms are situated such that the θ(Br··N··Br) angle is 79.7° in the triad.
The interaction energies in the next column were evaluated by comparing the energy of each complex with that of the sum of monomers which were left in their geometries within the confines of the dimer. These quantities verify the preference of the single XB within the dimer (16.30 kcal/mol) to that of the triad (15.60 kcal/mol) where the two FBr units must share a single N lone pair. The electron density at the Br··N bond critical points in the next column amplify this energetic preference, in that the single linear Br··N density of 0.056 au greatly exceeds the 0.021 au in each of the two bifurcated XBs within the trimer. The same pattern is evident in the NBO values of E(2) for the charge transfer from the N lone pair into the σ *(FBr) antibonding orbitals. This perturbation energy of 51.1 kcal/mol in the linear dimer is much higher than the 9.0 kcal/mol for each of the two XBs in the bifurcated trimer.
It is anticipated that there are two factors contributing to the bond weakening in the triads. In the first place, each bond is distorted from its preferred linear arrangement which can only be achieved within the dimer. The second factor arises since there are two Lewis acids, both competing for the same N lone pair. The fact that the base must serve as electron donor to two acids simultaneously would lead to a negative cooperativity. The third row of each section of Table 1 allows a separation of these two factors. The quantities reported there represent those of a dimer containing only a single acid unit, but it is placed in the same position it occupies within the trimer. In other words, these values document the result of the geometric distortion from the optimized dimer, but absent the negative cooperativity that would arise were the second acid unit physically present.
Again taking the (FBr)2NH3 system as our example, the placement of a single FBr molecule in the position it occupies in the trimer reduces the interaction energy of the fully optimized dimer (16.30 kcal/mol) down to 9.43 kcal/mol. So, the bifurcation configuration cuts the interaction energy nearly in half. The reduction in ρBCP is even greater, and E(2) is reduced to only about 1/5 its optimized value. Adding in the second FBr molecule does not double the interaction energy, as would be the case in the absence of any cooperative effects. Instead of 18.86 kcal/mol (twice the value of the individual distorted Br··N), the interaction energy of the trimer is only 15.60, so one might assess the negative cooperativity as their difference, i.e., 3.26 kcal/mol.
Very similar trends can be observed when FBr is replaced by FI. The halogen bonds in the bifurcated triad are stretched relative to the fully optimized linear dimer. The interaction energy within the trimer is less than in the linear dimer, and there is a substantial negative cooperativity. The bifurcated H-bonds associated with a pair of HF molecules are slightly different. There is negligible weakening of the interaction energy on going from dimer to trimer, and the two HF molecules are placed a bit closer together, with θ(H··N··H) only 66°. However, the negative cooperativity remains.
There is a qualitative change when the sp3 lone pair of NH3 is replaced by the sp-hybridization of NCH. In the first place, the latter is a weaker base, so interaction energies, and their AIM and NBO markers are similarly reduced, and there is a concomitant stretching of the corresponding XBs. However, there is also the change that the interaction energies of the triads are larger than those within the fully optimized linear dimers. The 7.2 kcal/mol interaction energy of the linear FBr··NCH dimer rises to 10.5 kcal/mol when a second FBr molecule is added, and the system is forced into a bifurcated arrangement. In other words, there is an energetic advantage to bifurcation about the NCH lone pair. Part of the reason for this change is that the energetic cost of moving each acid from its fully optimized position to that which it occupies in the bifurcated trimer is less stringent for NCH. For example, the interaction energy of the FBr··NCH dimer is 7.20 kcal/mol, which is reduced only slightly to 6.07 kcal/mol within the context of the trimer structure. The negative cooperativity observed for NH3 occurs as well for NCH, so that is a common feature of these bifurcated trimers.

2.2. Shared σ-Hole of Single Acid

An alternate form of bifurcation would have two separate bases sharing a single σ-hole on a Lewis acid molecule. Systems of this sort are depicted in Figure 3 where two NCH bases are placed accordingly around the FBr, FI, and H-bonding FH molecules. As in the earlier cases, these structures are not the result of a full geometry optimization which aligns one HCN along the FX axis, while the second moves around to bind to the first. The bifurcated structures were derived by enforcing a degree of symmetry wherein the two FX··N angles were restricted to be equal to one another. Bifurcated systems involving NH3 rather than NCH were not possible even with such a restriction, since one NH3 would rotate so as to engage in NH··N H-bonding with the other.
The geometric and other properties of these bifurcated triads are reported in Table 2, along with the fully optimized dimers containing a linear XB or HB. Unlike the addition of a second FX molecule to the FX··NCH dyads which raised the overall interaction energy, the addition instead of a second NCH base reduces this quantity. Taking FI···NCH as an example, its XB energy is 10.2 kcal/mol. The addition of a second FI raises the total interaction energy of the triad up to 13.5 kcal/mol, whereas this quantity drops down to 6.8 kcal/mol upon addition of a second NCH. Part of the reason for this behavior is the more precipitous drop in the individual interaction energies when the deformation from optimized linear dimer to the bent structure is taken into account. As may be observed in Table 1, this distortion results in FX··NCH interaction energies of some 6–8 kcal/mol, while this quantity drops down to 3–5 kcal/mol in Table 2 when the geometry is influenced by a pair of NCH bases.
This opposite behavior is reflected to some extent in the geometries. The addition of a second Lewis acid lengthens the intermolecular distances by a fairly small amount, between 0.10 and 0.25 Å. These stretches are much longer when it is a second NCH that is added, in the range between 0.33 and 0.46 Å. The latter longer intermolecular distances are a major factor in the reduction of the distorted dimer geometries mentioned above. The pattern continues in the other measures of bonding as well. Both the bond critical point densities and the NBO E(2) quantities are reduced much more by the addition of a second NCH than by a second FX unit. The smaller θ angles in Table 2 as compared to Table 1 show that the two NCH bases can approach one another more closely than the pair of FX molecules. The θ(N··X··N) angles of 59–61° are considerably smaller than the θ(X··N··X) angles of 88–92°. However, there is no such distinction for FH where the corresponding angles are comparable to one another.

2.3. Effects of Complexation on Monomers

It is understood that the formation of a HB or XB will affect the internal properties of each subunit. For example, the internal bond length of a FH acid will stretch when engaged in a H-bond, and its stretching frequency shift to the red. These changes will be accompanied by a downfield shift of the bridging proton’s NMR signal, caused in part by its deshielding. The effects of both dimerization and trimerization on the individual monomers are reported in Table 3. The Δr quantities listed there show that complexation results in a stretching of each FX covalent bond. This stretching is accentuated for NH3 as compared to NCH, consistent with the stronger basicity of the former. The elongations within the halogen-containing FX are more pronounced than within FH. The degree of stretching is reduced when the fully linear optimized dimer is distorted into the bent bifurcated structures. It might be noted as well that the reductions imposed in the trimers relative to the dyads are smaller when it is a second NCH that is added to FX··NCH, as compared to a second FX.
The next column of Table 3 displays the change in the associated ν(FX) stretching frequency, which is a negative red shift in the majority of cases. These shifts are particularly pronounced for FH with the much lighter H as compared to the halogenated FX units. Like the FX elongations, the shifts are larger for NH3 as compared to NCH. Also consistent with bond length changes, there is less of an effect arising from addition of a second NCH to the dimer than occurs from addition of FX.
The last three columns of Table 3 contain the changes in chemical shielding arising on the X, F, and N atoms that accompany the formation of the indicated dyads and triads. The shieldings of the halogen atoms are increased by the complexation, particularly Br, while H is deshielded. The F atom is deshielded, as is the basic N atom. In the majority of cases, these NMR properties share the common behavior that the change in shielding occurring within the dimer is attenuated within the trimer. With respect to the Lewis acid X and F atoms, one sees again the pattern noted earlier for the energetics, that the addition of a second NCH unit has a greater effect than does a second FX.

2.4. Understanding the Trends

The data show that enforcing a symmetric sort of bifurcated pair of bonds is energetically inferior to a single linear such bond in most cases. The exceptions arise when a pair of H or X bonding acids interact with a single weak NCH base, wherein the total interaction energy of the triad is slightly larger than that of a single linear bond. A refined examination of the numerical data indicates that there is less of an energetic cost when two acids share a single lone pair, than the reverse situation wherein a single acid σ-hole must accommodate two bases. It is not only the energy itself, but the other indicators of noncovalent bond strength, AIM, and NBO markers, as well as geometric and spectroscopic data, that affirm the angular deformations necessary to align two bases with a single σ-hole is more damaging than when two acids share a single base lone pair.
It is possible to understand the underlying reasons for this distinction by examination of the molecular electrostatic potential (MEP) that surrounds each molecule. Figure 4a–c presents the MEP of each of the three Lewis acid molecules in the form of a contour plot. The blue curve indicates the van der Waals surface of each, as measured by an isodensity surface with ρ = 0.001 au. Focusing attention on the right side of each which encompasses its σ-hole, one can see that this vdW surface crosses several different contours, an indication of its anisotropy. In other words, the MEP drops in value quickly as the reference point moves away from the F-X axis. This anisotropy appears to be more dramatic for FBr and HI than for FH in Figure 4c. This behavior can be contrasted with the N lone pair region of the NH3 and NCH bases, whose negative MEP with its excess density might be termed a σ-lump. Their MEPs, displayed as the broken contours on the right sides of Figure 4d,e, respectively, undergo less crossing of the blue vdW surface, particularly for NCH.
An alternate manner of viewing the comparative anisotropies is via Figure 5 which plots the MEP on a circle of constant distance from the Br, I, and H atoms of the acids, and the N atom of the two bases. (Note that the sign of the MEP has been reversed for the two bases for purposes of comparison with the acids.) The red (FBr) and purple (FI) curves display the greatest drop as θ increases away from the F-X axis, while the black FH curve is a bit flatter. The blue curve associated with the anisotropy of the σ-lump of NH3 is much steeper than that for NCH.
These visualizations conform with the numerical data. The greater anisotropy of the NH3 σ-lump region as compared to NCH accounts for the observation that the angular deformations required to add two acids is more energetically costly for the former. Moreover, the high anisotropy of the FBr and FI σ-holes explains the large drop in interaction energy when two bases are placed around them that are forced to lie off of the F-X axis.
Of course, the electrostatic interaction represents only one chapter of the story. As represented through the NBO formalism, the bonding is partly due also to charge transfer from the N lone pair orbital to the unoccupied σ *(FX) orbital of the acid. This component too ought to have an angular dependence that diminishes as the bond is forced from linearity. The angular characteristics of the relevant NBO orbitals are provided in Figure 6 in the form of contour plots, where the blue curve again indicates the vdW surface of each unit. The diagrams are not unlike the MEPs of Figure 4 in that the σ * orbitals of FBr and FI cross the most contour lines. The σ *(FH) orbital is a bit more isotropic. Most isotropic of all are the N lone pair orbitals in Figure 6d,e where the blue vdW surface nearly traces out a single contour line for much of its path.
The near isotropy of the NCH lone pair is particularly obvious by the nearly flat nature of the green curve in Figure 7. The blue NH3 curve is a bit more anisotropic, reminiscent of the MEP in Figure 5. Also in common with the MEPs, it is again the FBr and FI curves which diminish the most as one shifts further from the F-X axis. So, both the MEP and individual orbital diagrams demonstrate that much of the behavior of the bifurcated systems is consistent with both an electrostatic or orbital overlap interpretation.
The anisotropies of the MEP and lone pair orbital of the bases fit squarely with the geometries of the bifurcated systems in Figure 2. In all cases, whether the two acids are FBr, FI, or FH, these acids are situated further apart from one another when paired with NCH than with NH3. Figure 5 and Figure 7 are consistent with this wider spacing as the MEP and lone pair orbital around the N atom both deteriorate more slowly as one shifts away from the linear arrangement of NCH vs. NH3, so can tolerate a wider displacement, which in turn minimizes repulsions between the two acids. The same can be said for the structures exhibited in Figure 3 where two bases are paired with a single acid. The two bases are spaced most closely together for FI, which also has the quickest deterioration of the σ-hole and σ * orbital. These two properties diminish much more slowly for FH, which allows 30° greater separation of the two bases in Figure 3c, again minimizing base–base repulsions.
Still another factor in the ability of a partner molecule to locate itself off of the primary axis is the anisotropy of the electron density. A high electron density would tend to increase the exchange repulsion with its partner. An example may be drawn from the purple FI plot of the electron density around FI, at a constant distance of 3 Å from the I atom in Figure 8. This density rises rather dramatically with θ which would impede a base from shifting off of the FI axis, thus hindering the ability of FI to engage with two bases, each of which would have to move in this direction. The minimum in the density along the FI axis is consistent with the basic idea of a depleted σ-hole that sucks density out of this region and the concept of polar flattening. A similar rise, albeit not quite as steep, is observed in the red FBr curve, and the increase is even more gradual for FH. The bases behave quite differently. The change in ρ as θ increases is either barely perceptible as in the case of NCH, or even a decrease for NH3. This maximum at θ = 0 comports with the idea of a σ-lump that aligns with the N lone pair. The absence of a higher electron density resulting from a displacement off the lone pair direction is another factor accounting for the greater ability of the bases to share a lone pair with two Lewis acids, as compared to the diminished capability of a Lewis acid to share a single σ-hole with two bases. This same feature also helps explain the greater separation of the two FX acids in Figure 2a,b (~90°) as compared to the 60° separation of the two bases in Figure 3a,b.
With regard to the electrostatic potential surrounding each monomer, an alternate view of this quantity is known as the potential acting on one electron in a molecule (PAEM) [66] which contains exchange as well as Coulomb elements. This quantity has seen some applications in systems containing halogen and other sorts of noncovalent bonds [67,68]. The PAEMs of the five monomers are displayed in Figure 9 as a function of the displacement from their symmetry line, all with constant distance of 3 Å from the central atom. It might be noted first that unlike the MEP which is positive in the σ-hole regions of the Lewis acids, and negative along the lone pairs, the PAEM is negative for all. Like the MEP, the PAEM of FI and FBr have maxima along the F-X axis, but they do not decline monotonously, instead reaching minima at about 30°. The PAEM of FH is very flat until finally tailing off near 60°. Bases NCH and NH3 behave oppositely from one another. While the latter has its least negative value along the lone pair direction, this quantity is most negative along the lone pair of NCH.

3. Discussion

According to the calculations discussed above, a bifurcated arrangement which encompasses a pair of distorted noncovalent bonds is less favorable than is a single, fully optimized linear arrangement even though the former has two such bonds to the single bond in the latter. This situation applies to H-bonds as well as halogen bonds. In fact, a geometry encompassing a bifurcated pair of noncovalent bonds to either a common lone pair or a common σ-hole does not represent a true minimum on the potential energy surface of these triads. The exceptions occur when a pair of H or X bonding acids interact with a single weak NCH base, wherein the total interaction energy of the triad is larger than that of a single linear bond, even if only by a small amount.
This situation arises due to several principal factors. Angular deformations implicit to bifurcated bonds cause a significant weakening. Coupled to this issue is the negative cooperativity that arises when a single molecule acts as either double electron donor or acceptor. There are also repulsive interactions between the pair of molecules competing to bond with the central unit. These repulsions are less than 2 kcal/mol for a pair of FBr or FI units but are a bit larger for two FH or NCH molecules. Were these repulsions to be removed from the triad interaction energies, the remaining quantity representing only the attractions between the central NH3 and the two FX acids is no longer less than that within a fully optimized linear FX··NH3 dimer. However, even removal of the repulsion between the two NCH bases in the FX(NCH)2 triads does not raise the remaining pair of FX··N attractions up to the level of a single linear FX··NCH dimer, underscoring the lesser ability of a σ-hole to accommodate two separate bases.
All energetics have been corrected for basis set superposition error by the counterpoise protocol. Failure to make this correction enlarges each of the interaction energies in Table 1 and Table 2 by a small amount, roughly 1 kcal/mol. Moreover, the uncorrected energies obey the same trends as the properly corrected values. As a second point, the interaction energies which are the focus here refer to individual monomers in the geometries they adopt within the complex. If instead, the referrals are to the monomers in their fully optimized geometries, the resulting binding energies differ very little from the interaction energies, by less than 1 kcal/mol. So, again the trends displayed in the tables here would be unchanged if binding energies were considered instead.
Although there is an innate energetic preference in many cases for a single noncovalent XB or HB as opposed to a bifurcated arrangement, the latter obviously does arise rather commonly within the confines of larger macromolecular systems, or when crystal packing forces are present. Configurations found there reflect the fact that external forces preclude the optimal linear arrangement. In such a situation, it makes sense that a second Lewis acid can take advantage of the unoccupied side of the base’s lone pair, which is preferable to the former lying dormant. After all, a bent bond is preferable to no bond at all. The data in Table 1 and Table 2 support this contention, in that there is a strong energetic advantage in adding a second bond to an already bent bond of the same type. Taking the XB formed between FBr and NH3 as an example, the interaction energy of a bent XB is 9.4 kcal/mol, which grows to 15.6 kcal/mol upon formation of a second FBr··N XB.
The bifurcated bonds considered here are idealized in the sense that they are symmetric. That is, the two CN··X angles when a pair of FX acids interact with the NCH lone pair were held to be equal to one another, as was true for the NH3 base, as presented in Figure 2. Likewise, the two FX··N angles in Figure 3 were restricted to be equivalent. This idealized situation obviously does not represent the general bifurcated case where there is no restriction to a symmetric structure. What the work presented here has demonstrated is that there will be a natural tendency for one of the two bonds to lean toward linearity, slowly squeezing the other bond out of its way if conditions permit. However, if this approach toward linearity is not precluded by circumstances beyond its control, a bifurcated arrangement comprising two noncovalent bonds would be the second-best scenario.
The work described here considers the case where the nucleophile contains only a single lone pair that must be shared between two Lewis acids. However, when this restriction is lifted, allotting more than one lone pair to the base, it becomes much easier for the latter to participate in multiple noncovalent bonds simultaneously. Recent examples can be drawn from the rapidly developing field of selective anion binding. The base in question is typically a halide anion containing not only several available lone pairs, but also a full negative charge. Bulfield and Huber [69] have reviewed the manner in which catalysis may be facilitated by the bifurcated halogen bonding of an anion to a bipodal receptor with two symmetrically disposed halogen atoms. Lim et al. [70] have extended these same ideas to chalcogen bonds, where again an anion is bound by a bifurcated pair of noncovalent bonds to a bipodal receptor. Quantum calculations have explored this sort of bonding, and extended the ideas to bifurcated pnicogen, chalcogen, and tetrel bonds as well [71,72]. Other studies have extended these ideas to more than two noncovalent bonds to a given anion, again containing multiple lone pairs, in a general sense where the two electron-accepting atoms are located on different molecules [73,74,75].
There is precedent for the finding that is confirmed here that the halogen bond tends more toward linearity than does the H-bond [76]. For purposes of interpretation, the calculations presented here have focused on the directionality of the MEP and the electron density, expressed as either its total or individual MOs, to understand the ability of different sorts of bonds to accommodate bifurcated interactions. Indeed, there is precedent in the literature for the idea that both electrostatics and charge transfer/induction are important components of the anisotropy of noncovalent bonds [77,78,79,80,81,82]. There is likewise accumulating evidence that electron repulsion also plays an important role [77,78,83,84,85,86,87,88].
The work presented here has focused on halogen bonds, and their potential to share their single σ-hole with a pair of nucleophiles. Specifically, the halogen atoms considered here are all univalent in the sense that they are involved in only one covalent bond, in this case to F. However, halogen atoms can participate in hypervalent situations where they may be involved in more than one covalent bond. In such a situation one can anticipate that there will be one σ-hole that corresponds to each such internal bond. Kirshenboim and Kozuch [89] explored this scenario, and the manner in which the multiple σ-holes might affect the halogen bonding. Heinen et al. [90] later documented these ideas in the particular case of the hypervalent I atom through both experimental and computational means. The ability of a single Lewis acid atom with multiple σ-holes to engage in several noncovalent bonds has been explored also in the framework of chalcogen, pnicogen, and tetrel bonds [91,92,93,94,95,96,97]. These sorts of hypervalent bonding situations, with multiple σ-holes, are thus much more amenable to bifurcated arrangements than is the univalent scenario.

4. Computational Methods

All quantum calculations were carried out within the framework of the Gaussian-09 set of codes [98]. DFT procedures were used to include electron correlation via the M06-2X functional along with the aug-cc-pVDZ basis set. Relativistic effects which may be important for the heavy I atom were incorporated through the aug-cc-pVDZ-PP pseudopotential [99,100]. There exists ample evidence of the ability of this level of theory to accurately reproduce results from much higher levels when applied to noncovalent interactions of this type [30,101,102,103,104,105,106,107,108,109]. The interaction energy, Eint, of each complex was computed as the difference between the energy of the complex and the sum of the energies of monomers in the geometries they adopt within the complex. Basis set superposition error was corrected via the standard counterpoise [110,111,112] protocol. Molecular electrostatic potentials (MEPs) were computed in the framework of the Multiwfn program [113]. The topology of the electron density was analyzed with the aid of AIM procedures, applying the AIMALL program [114]. The natural bond orbital (NBO) prescription [115,116] was used to elucidate charge transfers and associated energetics between individual orbitals.

5. Conclusions

With some exceptions, a bifurcated pair of either H-bonds or halogen bonds is energetically inferior to a single such bond within a dimer. There is a larger energetic cost to forcing the single σ-hole of an acid to interact with a pair of bases, than the other way around where two acids engage with the lone pair of a single base. This distinction can be traced to details of the MEP that surrounds the various molecules, as well as to the shape of the pertinent molecular orbitals, and the total electron density. The acid’s σ-hole is more anisotropic than is the negative MEP at the site of the base’s lone pair. The same is true of the σ *(FX) antibonding orbital which is more directional than is the N lone pair orbital. On both of these scores, the weaker NCH base is less anisotropic than is NH3, which allows the former to better accommodate a bifurcated arrangement with two acid molecules. Likewise, the lesser anisotropy of FH as compared to FBr and FI allows H-bonds to better share a single proton than can XBs share their X atom. A shift of a base off of the bond axis of the Lewis acid would encounter a higher electron density, and associated exchange repulsion, another factor inhibiting the formation of a bifurcated bond with two bases, particularly for FI. In contrast, there is no such exchange repulsion factor inhibiting the sort of nonlinear geometry required to link two acids to a single base.

Funding

This material is based upon work supported by the National Science Foundation under Grant No. 1954310.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original data presented in this study are available by request to the author.

Conflicts of Interest

The author declares no conflict of interest.

Sample Availability

Not available.

References

  1. Pauling, L. The Nature of the Chemical Bond; Cornell University Press: Ithaca, NY, USA, 1940; p. 450. [Google Scholar]
  2. Pauling, L.; Corey, R.B. Configurations of polypeptide chains with favored orientations around single bonds: Two new pleated sheets. Proc. Nat. Acad. Sci. USA 1951, 37, 729–740. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Pauling, L.; Corey, R.B. Two hydrogen-bonded spiral configurations of the polypeptide chain. J. Am. Chem. Soc. 1950, 72, 5349. [Google Scholar] [CrossRef]
  4. Pauling, L.; Corey, R.B.; Branson, H.R. The structure of proteins: Two hydrogen-bonded helical configurations of the polypeptide chain. Proc. Nat. Acad. Sci. USA 1951, 37, 205–211. [Google Scholar] [CrossRef] [Green Version]
  5. Pimentel, G.C.; McClellan, A.L. The Hydrogen Bond; Freeman: San Francisco, CA, USA, 1960. [Google Scholar]
  6. Sutor, D.J. The C-H...O hydrogen bond in crystals. Nature 1962, 195, 68–69. [Google Scholar] [CrossRef]
  7. Hamilton, W.C.; Ibers, J.A. Hydrogen Bonding in Solids; W. A. Benjamin: New York, NY, USA, 1968; p. 284. [Google Scholar]
  8. Vinogradov, S.N.; Linnell, R.H. Hydrogen Bonding; Van Nostrand-Reinhold: New York, NY, USA, 1971. [Google Scholar]
  9. Joesten, M.D.; Schaad, L.J. Hydrogen Bonding; Marcel Dekker: New York, NY, USA, 1974; p. 622. [Google Scholar]
  10. Schuster, P.; Zundel, G.; Sandorfy, C. (Eds.) The Hydrogen Bond. Recent Developments in Theory and Experiments; North-Holland Publishing Co.: Amsterdam, The Netherlands, 1976. [Google Scholar]
  11. Jeffrey, G.A.; Saenger, W. Hydrogen Bonding in Biological Structures; Springer: Berlin, Germany, 1991. [Google Scholar]
  12. Scheiner, S. Hydrogen Bonding: A Theoretical Perspective; Oxford University Press: New York, NY, USA, 1997; p. 375. [Google Scholar]
  13. Scheiner, S. Fundamental Features of Hydrogen Bonds. In Pauling’s Legacy—Modern Modelling of the Chemical Bond; Maksic, Z.B., Orville-Thomas, W.J., Eds.; Elsevier: Amsterdam, The Netherlands, 1997; Volume 6, pp. 571–591. [Google Scholar]
  14. Scheiner, S. Pauling’s Legacy—Modern Modelling of the Chemical Bond; Maksic, Z.B., Orville-Thomas, W.J., Eds.; Elsevier: Amsterdam, The Netherlands, 1997. [Google Scholar]
  15. Desiraju, G.R.; Steiner, T. The Weak Hydrogen Bond in Structural Chemistry and Biology; Oxford: New York, NY, USA, 1999; p. 507. [Google Scholar]
  16. Gilli, G.; Gilli, P. The Nature of the Hydrogen Bond; Oxford University Press: Oxford, UK, 2009; p. 313. [Google Scholar]
  17. Arunan, E.; Desiraju, G.R.; Klein, R.A.; Sadlej, J.; Scheiner, S.; Alkorta, I.; Clary, D.C.; Crabtree, R.H.; Dannenberg, J.J.; Hobza, P.; et al. Definition of the Hydrogen Bond. Pure Appl. Chem. 2011, 83, 1637–1641. [Google Scholar] [CrossRef]
  18. Chand, A.; Sahoo, D.K.; Rana, A.; Jena, S.; Biswal, H.S. The Prodigious Hydrogen Bonds with Sulfur and Selenium in Molecular Assemblies, Structural Biology, and Functional Materials. Acc. Chem. Res. 2020, 53, 1580–1592. [Google Scholar] [CrossRef]
  19. Kumar, M.; Francisco, J.S. Evidence of the Elusive Gold-Induced Non-classical Hydrogen Bonding in Aqueous Environments. J. Am. Chem. Soc. 2020, 142, 6001–6006. [Google Scholar] [CrossRef]
  20. Chand, A.; Biswal, H.S. Hydrogen Bonds with Chalcogens: Looking Beyond the Second Row of the Periodic Table. J. Indian Inst. Sci. 2020, 100, 77–100. [Google Scholar] [CrossRef]
  21. Mishra, K.K.; Singh, S.K.; Kumar, S.; Singh, G.; Sarkar, B.; Madhusudhan, M.S.; Das, A. Water-Mediated Selenium Hydrogen-Bonding in Proteins: PDB Analysis and Gas-Phase Spectroscopy of Model Complexes. J. Phys. Chem. A 2019, 123, 5995–6002. [Google Scholar] [CrossRef]
  22. Schmidbaur, H. Proof of Concept for Hydrogen Bonding to Gold, Au⋅⋅⋅H−X. Angew. Chem. Int. Ed. 2019, 58, 5806–5809. [Google Scholar] [CrossRef]
  23. Wang, P.; Xu, H.-G.; Cao, G.-J.; Zhang, W.-J.; Xu, X.-L.; Zheng, W.-J. Nonconventional Hydrogen Bonds between Silver Anion and Nucleobases: Size-Selected Anion Photoelectron Spectroscopy and Density Functional Calculations. J. Phys. Chem. A 2017, 121, 8973–8981. [Google Scholar] [CrossRef] [PubMed]
  24. Mundlapati, V.R.; Sahoo, D.K.; Ghosh, S.; Purame, U.K.; Pandey, S.; Acharya, R.; Pal, N.; Tiwari, P.; Biswal, H.S. Spectroscopic Evidences for Strong Hydrogen Bonds with Selenomethionine in Proteins. J. Phys. Chem. Lett. 2017, 8, 794–800. [Google Scholar] [CrossRef] [PubMed]
  25. Scheiner, S. Comparison of CH···O, SH···O, Chalcogen, and Tetrel Bonds Formed by Neutral and Cationic Sulfur-Containing Compounds. J. Phys. Chem. A 2015, 119, 9189–9199. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Zabardasti, A.; Talebi, N.; Kakanejadifard, A.; Saki, Z. The B–C and C–C bonds as preferred electron source for H-bond and Li-bond interactions in complex pairing of C4B2H6 with HF and LiH molecules. Struct. Chem. 2016, 27, 573–581. [Google Scholar] [CrossRef]
  27. Møller, K.H.; Hansen, A.S.; Kjaergaard, H.G. Gas Phase Detection of the NH–P Hydrogen Bond and Importance of Secondary Interactions. J. Phys. Chem. A 2015, 119, 10988–10998. [Google Scholar] [CrossRef]
  28. Viana, R.B.; da Silva, A.B.F. Interaction between PH3 and small water clusters: Understanding the electronic and spectroscopic properties. Comput. Theor. Chem. 2015, 1059, 35–44. [Google Scholar] [CrossRef]
  29. Sanchez-de-Armas, R.; Ahlquist, M.S.G. On the nature of hydrogen bonds to platinum(II)—Which interaction can predict their strength? Phys. Chem. Chem. Phys. 2015, 17, 812–816. [Google Scholar] [CrossRef]
  30. Gu, Y.; Kar, T.; Scheiner, S. Comparison of the CH...N and CH...O interactions involving substituted alkanes. J. Mol. Struct. 2000, 552, 17–31. [Google Scholar] [CrossRef]
  31. Orlova, G.; Scheiner, S. Intermolecular MH...HF bonding in monohydride Mo and W complexes. J. Phys. Chem. A 1998, 102, 260–269. [Google Scholar] [CrossRef]
  32. Adhikari, U.; Scheiner, S. Comparison of P...D (D=P,N) with other noncovalent bonds in molecular aggregates. J. Chem. Phys. 2011, 135, 184306. [Google Scholar] [CrossRef]
  33. Grabowski, S.J. A−H…σ Hydrogen Bonds: Dihydrogen and Cycloalkanes as Proton Acceptors. ChemPhysChem 2019, 20, 565–574. [Google Scholar] [CrossRef] [PubMed]
  34. Alkorta, I.; Martín-Fernández, C.; Montero-Campillo, M.M.; Elguero, J. Hydrogen-Bonding Acceptor Character of Be3, the Beryllium Three-Membered Ring. J. Phys. Chem. A 2018, 122, 1472–1478. [Google Scholar] [CrossRef] [PubMed]
  35. Lu, N.; Chung, W.-C.; Ley, R.M.; Lin, K.-Y.; Francisco, J.S.; Negishi, E.-I. Molecularly Tuning the Radicaloid N–H···O═C Hydrogen Bond. J. Phys. Chem. A 2016, 120, 1307–1315. [Google Scholar] [CrossRef] [PubMed]
  36. Mirzaei, M.; Nikpour, M.; Bauzá, A.; Frontera, A. On the Importance of C-H/π and C-H...H-C Interactions in the Solid State Structure of 15-Lipoxygenase Inhibitors Based on Eugenol Derivatives. ChemPhysChem 2015, 16, 2260–2266. [Google Scholar] [CrossRef]
  37. Grabowski, S.J. Dihydrogen bond and X–H…σ interaction as sub-classes of hydrogen bond. J. Phys. Org. Chem. 2013, 26, 452–459. [Google Scholar] [CrossRef]
  38. Li, Q.-Z.; Kou, H.; Li, R.; Li, W.-Z.; Cheng, J.-B. Is a Hm-1X (X = O, S, Se, m = 1; X = N, m = 2; X = C, m = 3) radical a better proton donor than HmX–H in hydrogen bonding? Comput. Theor. Chem. 2011, 976, 83–87. [Google Scholar] [CrossRef]
  39. Kar, T.; Bettinger, H.F.; Scheiner, S.; Roy, A.K. Noncovalent π−π Stacking and CH---π Interactions of Aromatics on the Surface of Single-Wall Carbon Nanotubes: An MP2 Study. J. Phys. Chem. C 2008, 112, 20070–20075. [Google Scholar] [CrossRef]
  40. Taylor, R.; Kennard, O.; Versichel, W. Geometry of the N-H...O=C hydrogen bond. 2. Three-center (“bifurcated”) and four-center (“trifurcated”) bonds. J. Am. Chem. Soc. 1984, 106, 244–248. [Google Scholar] [CrossRef]
  41. Li, A.Y. Theoretical study of linear and bifurcated H-bonds in the systems Y...H2CZn (n = 1, 2; Z = O, S, Se, F, Cl, Br; Y = Cl−, Br−). J. Mol. Struct. (Theochem.) 2008, 862, 21–27. [Google Scholar] [CrossRef]
  42. Blanco, S.; Lopez, J.C.; Lesarri, A.; Caminati, W.; Alonso, J.L. Bifurcated CH2...O and (C-H)2...F-C weak hydrogen bonds: The oxirane-difluoromethane complex. ChemPhysChem 2004, 5, 1779–1782. [Google Scholar] [CrossRef]
  43. Cense, J.M.; Agafonov, V.; Ceolin, R.; Ladure, P.; Rodier, N. Crystal and molecular structure analysis of flutamide. Bifurcated helicoidal C-H...O hydrogen bonds. Struct. Chem. 1994, 5, 79–84. [Google Scholar] [CrossRef]
  44. Duan, X.; Scheiner, S. Behavior of interaction energy and intramolecular bond stretch in linear and bifurcated hydrogen bonds. Int. J. Quantum Chem. QBS 1993, 20, 181–190. [Google Scholar] [CrossRef]
  45. van Hensbergen, B.; Block, R.; Jansen, L. Effect of direct and indirect exchange interactions on geometries and relative stabilities of H2O and H2S dimers in bifurcated, cyclic, and linear configurations. J. Chem. Phys. 1982, 76, 3161–3168. [Google Scholar] [CrossRef]
  46. Marsden, C.J.; Smith, B.J.; Pople, J.A.; Schaefer, H.F.; Radom, L. Characterization of the bifurcated structure of the water dimer. J. Chem. Phys. 1991, 95, 1825–1828. [Google Scholar] [CrossRef]
  47. Thomas, J.; Liu, X.; Jäger, W.; Xu, Y. Unusual H-Bond Topology and Bifurcated H-bonds in the 2-Fluoroethanol Trimer. Angew. Chem. Int. Ed. 2015, 54, 11711–11715. [Google Scholar] [CrossRef]
  48. Mishra, S.K.; Suryaprakash, N. Intramolecular hydrogen bonds involving organic fluorine in the derivatives of hydrazides: An NMR investigation substantiated by DFT based theoretical calculations. Phys. Chem. Chem. Phys. 2015, 17, 15226–15235. [Google Scholar] [CrossRef]
  49. Preimesberger, M.R.; Majumdar, A.; Aksel, T.; Sforza, K.; Lectka, T.; Barrick, D.; Lecomte, J.T.J. Direct NMR Detection of Bifurcated Hydrogen Bonding in the α-Helix N-Caps of Ankyrin Repeat Proteins. J. Am. Chem. Soc. 2015, 137, 1008–1011. [Google Scholar] [CrossRef] [Green Version]
  50. Fulara, A.; Dzwolak, W. Bifurcated hydrogen bonds stabilize fibrils of poyl(L-glutamic) acid. J. Phys. Chem. B 2010, 114, 8278–8283. [Google Scholar] [CrossRef]
  51. Quinonero, D. Sigma-hole carbon-bonding interactions in carbon-carbon double bonds: An unnoticed contact. Phys. Chem. Chem. Phys. 2017, 19, 15530–15540. [Google Scholar] [CrossRef]
  52. Esrafili, M.D.; Vessally, E.; Solimannejad, M. Symmetric bifurcated halogen bonds: Substituent and cooperative effects. Mol. Phys. 2016, 114, 3610–3619. [Google Scholar] [CrossRef]
  53. Novak, M.; Foroutan-Nejad, C.; Marek, R. Asymmetric bifurcated halogen bonds. Phys. Chem. Chem. Phys. 2015, 17, 6440–6450. [Google Scholar] [CrossRef]
  54. Oliveira, B.G.d.; Araújo, R.d.C.M.U.d.; Leite, E.S.; Ramos, M.N. A theoretical analysis of topography and molecular parameters of the CFCl3···O3 complex: Linear and bifurcate halogen-oxygen bonding interactions. Int. J. Quantum Chem. 2011, 111, 111–116. [Google Scholar] [CrossRef]
  55. Dabranskaya, U.; Ivanov, D.M.; Novikov, A.S.; Matveychuk, Y.V.; Bokach, N.A.; Kukushkin, V.Y. Metal-Involving Bifurcated Halogen Bonding C–Br···η2(Cl–Pt). Cryst. Growth Des. 2019, 19, 1364–1376. [Google Scholar] [CrossRef]
  56. Rozhkov, A.V.; Ivanov, D.M.; Novikov, A.S.; Ananyev, I.V.; Bokach, N.A.; Kukushkin, V.Y. Metal-involving halogen bond Ar–I⋯[dz2PtII] in a platinum acetylacetonate complex. CrystEngComm 2020, 22, 554–563. [Google Scholar] [CrossRef]
  57. Ivanov, D.M.; Novikov, A.S.; Ananyev, I.V.; Kirina, Y.V.; Kukushkin, V.Y. Halogen bonding between metal centers and halocarbons. Chem. Commun. 2016, 52, 5565–5568. [Google Scholar] [CrossRef] [Green Version]
  58. Bartashevich, E.; Troitskaya, E.; Pendás, Á.M.; Tsirelson, V. Understanding the bifurcated halogen bonding N...Hal...N in bidentate diazaheterocyclic compounds. Comput. Theor. Chem. 2015, 1053, 229–237. [Google Scholar] [CrossRef]
  59. Esrafili, M.D.; Mohammadian-Sabet, F.; Vessally, E. An ab initio study on substituent and cooperative effects in bifurcated fluorine bonds. Mol. Phys. 2017, 115, 278–287. [Google Scholar] [CrossRef]
  60. Massahi, S.; Ghobadi, M.; Nikoorazm, M. Exceptional bifurcated chalcogen bonding interaction between Ph2N2O2 and only one σ–hole on XCY (X = S, Se, Te and Y = O, S, Se, Te): A DFT study. Theor. Chem. Acc. 2020, 139, 162. [Google Scholar] [CrossRef]
  61. Esrafili, M.D.; Mohammadian-Sabet, F. Bifurcated chalcogen bonds: A theoretical study on the structure, strength and bonding properties. Chem. Phys. Lett. 2015, 634, 210–215. [Google Scholar] [CrossRef]
  62. Daolio, A.; Scilabra, P.; Di Pietro, M.E.; Resnati, C.; Rissanen, K.; Resnati, G. Binding motif of ebselen in solution: Chalcogen and hydrogen bonds team up. New J. Chem. 2020, 44, 20697–20703. [Google Scholar] [CrossRef]
  63. Esrafili, M.D.; Mohammadian-Sabet, F. Substituent effects on geometry and bonding properties of asymmetric bifurcated pnicogen bonds: A theoretical study. Chem. Phys. Lett. 2016, 650, 52–56. [Google Scholar] [CrossRef]
  64. Grabowski, S.J. Bifurcated Triel Bonds—Hydrides and Halides of 1,2-Bis(Dichloroboryl)Benzene and 1,8-Bis(Dichloroboryl)Naphthalene. Crystal 2019, 9, 503. [Google Scholar] [CrossRef] [Green Version]
  65. Esrafili, M.D.; Sadr-Mousavi, A. A computational study on the strength and nature of bifurcated aerogen bonds. Chem. Phys. Lett. 2018, 698, 1–6. [Google Scholar] [CrossRef]
  66. Zhao, D.-X.; Yang, Z.-Z. Investigation of the distinction between van der Waals interaction and chemical bonding based on the PAEM-MO diagram. J. Comput. Chem. 2014, 35, 965–977. [Google Scholar] [CrossRef] [PubMed]
  67. Bartashevich, E.; Matveychuk, Y.; Tsirelson, V. Identification of the Tetrel Bonds between Halide Anions and Carbon Atom of Methyl Groups Using Electronic Criterion. Molecules 2019, 24, 1083. [Google Scholar] [CrossRef] [Green Version]
  68. Bartashevich, E.V.; Yushina, I.D.; Stash, A.I.; Tsirelson, V.G. Halogen Bonding and Other Iodine Interactions in Crystals of Dihydrothiazolo(oxazino)quinolinium Oligoiodides from the Electron-Density Viewpoint. Cryst. Growth Des. 2014, 14, 5674–5684. [Google Scholar] [CrossRef]
  69. Bulfield, D.; Huber, S.M. Halogen Bonding in Organic Synthesis and Organocatalysis. Chem. Eur. J. 2016, 22, 14434–14450. [Google Scholar] [CrossRef]
  70. Lim, J.Y.C.; Liew, J.Y.; Beer, P.D. Thermodynamics of Anion Binding by Chalcogen Bonding Receptors. Chem. Eur. J. 2018, 24, 14560–14566. [Google Scholar] [CrossRef]
  71. Scheiner, S. Differential Binding of Tetrel-Bonding Bipodal Receptors to Monatomic and Polyatomic Anions. Molecules 2019, 24, 227. [Google Scholar] [CrossRef] [Green Version]
  72. Scheiner, S. Tetrel Bonding as a Vehicle for Strong and Selective Anion Binding. Molecules 2018, 23, 1147. [Google Scholar] [CrossRef] [Green Version]
  73. Scheiner, S.; Michalczyk, M.; Zierkiewicz, W. Coordination of anions by noncovalently bonded σ-hole ligands. Coord. Chem. Rev. 2020, 405, 213136. [Google Scholar] [CrossRef]
  74. Scheiner, S.; Michalczyk, M.; Wysokiński, R.; Zierkiewicz, W. Structures and energetics of clusters surrounding diatomic anions stabilized by hydrogen, halogen, and other noncovalent bonds. Chem. Phys. 2020, 530, 110590. [Google Scholar] [CrossRef]
  75. Scheiner, S.; Michalczyk, M.; Zierkiewicz, W. Structures of clusters surrounding ions stabilized by hydrogen, halogen, chalcogen, and pnicogen bonds. Chem. Phys. 2019, 524, 55–62. [Google Scholar] [CrossRef]
  76. Zheng, Y.-Z.; Deng, G.; Zhou, Y.; Sun, H.-Y.; Yu, Z.-W. Comparative Study of Halogen- and Hydrogen-Bond Interactions between Benzene Derivatives and Dimethyl Sulfoxide. ChemPhysChem 2015, 16, 2594–2601. [Google Scholar] [CrossRef]
  77. Tafipolsky, M. Challenging Dogmas: Hydrogen Bond Revisited. J. Phys. Chem. A 2016, 120, 4550–4559. [Google Scholar] [CrossRef]
  78. Kolár, M.; Hostaš, J.; Hobza, P. The strength and directionality of a halogen bond are co-determined by the magnitude and size of the s-hole. Phys. Chem. Chem. Phys. 2014, 16, 9987–9996. [Google Scholar] [CrossRef]
  79. Tsuzuki, S.; Wakisaka, A.; Ono, T.; Sonoda, T. Magnitude and origin of the attraction and directionality of the halogen bonds of the complexes of C6F5X and C6H5X (X = I, Br, Cl and F) with Pyridine. Chem. Eur. J. 2012, 18, 951–960. [Google Scholar] [CrossRef]
  80. Shields, Z.P.; Murray, J.S.; Politzer, P. Directional tendencies of halogen and hydrogen bonds. Int. J. Quantum Chem. 2010, 110, 2823–2832. [Google Scholar] [CrossRef]
  81. Riley, K.E.; Merz, K.M. Insights into the strength and origin of halogen bonding: The halobenzene-formaldehyde dimer. J. Phys. Chem. A 2007, 111, 1688–1694. [Google Scholar] [CrossRef]
  82. Riley, K.E.; Hobza, P. The relative roles of electrostatics and dispersion in the stabilization of halogen bonds. Phys. Chem. Chem. Phys. 2013, 15, 17742–17751. [Google Scholar] [CrossRef]
  83. Riley, K.E. Critical comparison of RX...Y and RH...Y directionality in halogen and hydrogen bonds using modern computational chemistry methods. Chem. Phys. Lett. 2020, 744, 137221. [Google Scholar] [CrossRef]
  84. Hill, J.G.; Legon, A.C. On the directionality and non-linearity of halogen and hydrogen bonds. Phys. Chem. Chem. Phys. 2015, 17, 858–867. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Kerdawy, A.E.; Murray, J.S.; Politzer, P.; Bleiziffer, P.; Heßelmann, A.; Görling, A.; Clark, T. Directional noncovalent interactions: Repulsion and dispersion. J. Chem. Theory Comput. 2013, 9, 2264–2275. [Google Scholar] [CrossRef]
  86. Stone, A.J. Are halogen bonded structures electrostatically driven? J. Am. Chem. Soc. 2013, 135, 7005–7009. [Google Scholar] [CrossRef]
  87. Adhikari, U.; Scheiner, S. Sensitivity of pnicogen, chalcogen, halogen and H-bonds to angular distortions. Chem. Phys. Lett. 2012, 532, 31–35. [Google Scholar] [CrossRef]
  88. Riley, K.E.; Vazquez, M.; Umemura, C.; Miller, C.; Tran, K.-A. Exploring the (Very Flat) Potential Energy Landscape of R−Br⋅⋅⋅π Interactions with Accurate CCSD(T) and SAPT Techniques. Chem. Eur. J. 2016, 22, 17690–17695. [Google Scholar] [CrossRef] [Green Version]
  89. Kirshenboim, O.; Kozuch, S. How to Twist, Split and Warp a σ-Hole with Hypervalent Halogens. J. Phys. Chem. A 2016, 120, 9431–9445. [Google Scholar] [CrossRef]
  90. Heinen, F.; Engelage, E.; Cramer, C.J.; Huber, S.M. Hypervalent Iodine(III) Compounds as Biaxial Halogen Bond Donors. J. Am. Chem. Soc. 2020, 142, 8633–8640. [Google Scholar] [CrossRef]
  91. Zierkiewicz, W.; Wysokiński, R.; Michalczyk, M.; Scheiner, S. Chalcogen bonding of two ligands to hypervalent YF4 (Y = S, Se, Te, Po). Phys. Chem. Chem. Phys. 2019, 21, 20829–20839. [Google Scholar] [CrossRef]
  92. Michalczyk, M.; Zierkiewicz, W.; Wysokiński, R.; Scheiner, S. Hexacoordinated Tetrel-Bonded Complexes between TF4 (T = Si, Ge, Sn, Pb) and NCH: Competition between σ- and π-Holes. ChemPhysChem 2019, 20, 959–966. [Google Scholar] [CrossRef] [Green Version]
  93. Scheiner, S.; Lu, J. Halogen, Chalcogen, and Pnicogen Bonding Involving Hypervalent Atoms. Chem. Eur. J. 2018, 24, 8167–8177. [Google Scholar] [CrossRef] [PubMed]
  94. Scheiner, S. Comparison of halide receptors based on H, halogen, chalcogen, pnicogen, and tetrel bonds. Faraday Disc. 2017, 203, 213–226. [Google Scholar] [CrossRef] [PubMed]
  95. Scheiner, S. Assembly of Effective Halide Receptors from Components. Comparing Hydrogen, Halogen, and Tetrel Bonds. J. Phys. Chem. A 2017, 121, 3606–3615. [Google Scholar] [CrossRef] [PubMed]
  96. Scheiner, S. Highly Selective Halide Receptors Based on Chalcogen, Pnicogen, and Tetrel Bonds. Chem. Eur. J. 2016, 22, 18850–18858. [Google Scholar] [CrossRef] [PubMed]
  97. Nepal, B.; Scheiner, S. Building a Better Halide Receptor: Optimum Choice of Spacer, Binding Unit, and Halosubstitution. ChemPhysChem 2016, 17, 836–844. [Google Scholar] [CrossRef] [Green Version]
  98. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09, Revision B.01; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  99. Feller, D. The role of databases in support of computational chemistry calculations. J. Comput. Chem. 1996, 17, 1571–1586. [Google Scholar] [CrossRef]
  100. Schuchardt, K.L.; Didier, B.T.; Elsethagen, T.; Sun, L.; Gurumoorthi, V.; Chase, J.; Li, J.; Windus, T.L. Basis Set Exchange:  A Community Database for Computational Sciences. J. Chem. Infor. Model. 2007, 47, 1045–1052. [Google Scholar] [CrossRef] [Green Version]
  101. Orlova, A.P.; Jasien, P.G. Halogen bonding in self-assembling systems: A comparison of intra- and interchain binding energies. Comput. Theor. Chem. 2018, 1139, 63–69. [Google Scholar] [CrossRef]
  102. Forni, A.; Pieraccini, S.; Franchini, D.; Sironi, M. Assessment of DFT Functionals for QTAIM Topological Analysis of Halogen Bonds with Benzene. J. Phys. Chem. A 2016, 120, 9071–9080. [Google Scholar] [CrossRef]
  103. Nziko, V.d.P.N.; Scheiner, S. Chalcogen Bonding between Tetravalent SF4 and Amines. J. Phys. Chem. A 2014, 118, 10849–10856. [Google Scholar] [CrossRef]
  104. Esrafili, M.D.; Vessally, E. A theoretical evidence for cooperative enhancement in aerogen-bonding interactions: Open-chain clusters of KrOF2 and XeOF2. Chem. Phys. Lett. 2016, 662, 80–85. [Google Scholar] [CrossRef]
  105. Liu, M.; Li, Q.; Scheiner, S. Comparison of tetrel bonds in neutral and protonated complexes of pyridineTF3 and furanTF3 (T = C, Si, and Ge) with NH3. Phys. Chem. Chem. Phys. 2017, 19, 5550–5559. [Google Scholar] [CrossRef] [Green Version]
  106. Scheiner, S. Systematic Elucidation of Factors That Influence the Strength of Tetrel Bonds. J. Phys. Chem. A 2017, 121, 5561–5568. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  107. Bauzá, A.; Alkorta, I.; Frontera, A.; Elguero, J. On the Reliability of Pure and Hybrid DFT Methods for the Evaluation of Halogen, Chalcogen, and Pnicogen Bonds Involving Anionic and Neutral Electron Donors. J. Chem. Theory Comput. 2013, 9, 5201–5210. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  108. Georg, H.C.; Fileti, E.E.; Malaspina, T. Ab initio study of weakly bound halogen complexes: RX...PH3. J. Mol. Model. 2013, 19, 329–336. [Google Scholar] [CrossRef] [PubMed]
  109. Lu, J.; Scheiner, S. Effects of Halogen, Chalcogen, Pnicogen, and Tetrel Bonds on IR and NMR Spectra. Molecules 2019, 24, 2822. [Google Scholar] [CrossRef] [Green Version]
  110. Boys, S.F.; Bernardi, F. The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors. Mol. Phys. 1970, 19, 553–566. [Google Scholar] [CrossRef]
  111. Brauer, B.; Kesharwani, M.K.; Martin, J.M.L. Some Observations on Counterpoise Corrections for Explicitly Correlated Calculations on Noncovalent Interactions. J. Chem. Theory Comput. 2014, 10, 3791–3799. [Google Scholar] [CrossRef]
  112. Latajka, Z.; Scheiner, S. Primary and secondary basis set superposition error at the SCF and MP2 levels: H3N--Li+ and H2O--Li+. J. Chem. Phys. 1987, 87, 1194–1204. [Google Scholar] [CrossRef] [Green Version]
  113. Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef]
  114. Keith, T.A. AIMALL; TK Gristmill Software: Overland Park, KS, USA, 2013. [Google Scholar]
  115. Reed, A.E.; Curtiss, L.A.; Weinhold, F. Natural population analysis. J. Chem. Phys. 1985, 83, 735–746. [Google Scholar] [CrossRef]
  116. Reed, A.E.; Curtiss, L.A.; Weinhold, F. Intermolecular interactions from a natural bond orbital, donor-acceptor viewpoint. Chem. Rev. 1988, 88, 899–926. [Google Scholar] [CrossRef]
Molecules 26 00350 sch001 550
Scheme 1. Diagrams of two sorts of bifurcated HBs. (A) shares a single proton, and (B) shares a single electron donor.
Scheme 1. Diagrams of two sorts of bifurcated HBs. (A) shares a single proton, and (B) shares a single electron donor.
Molecules 26 00350 sch001
Molecules 26 00350 g001 550
Figure 1. Optimized geometries of (a) FI···NH3 and (b) FH··NCH dyads. Distances in Å.
Figure 1. Optimized geometries of (a) FI···NH3 and (b) FH··NCH dyads. Distances in Å.
Molecules 26 00350 g001
Molecules 26 00350 g002 550
Figure 2. Geometries of triads wherein a pair of Lewis acids share a single lone pair of NH3 or NCH, optimized with constraint of equal θ(X··N··X) angles. Distances in Å and angles in degs. NH3 is shared by two (a) FBr, (b) FI and (c) FH. NCH is shared by two (d) FBr, (e) FI, and (f) FH.
Figure 2. Geometries of triads wherein a pair of Lewis acids share a single lone pair of NH3 or NCH, optimized with constraint of equal θ(X··N··X) angles. Distances in Å and angles in degs. NH3 is shared by two (a) FBr, (b) FI and (c) FH. NCH is shared by two (d) FBr, (e) FI, and (f) FH.
Molecules 26 00350 g002
Molecules 26 00350 g003 550
Figure 3. Geometries of triads wherein a pair of NCH bases share a single σ-hole of (a) FBr, (b) FI, and (c) FH, optimized with constraint of equal θ(N··X··N) angles. Distances in Å and angles in degs.
Figure 3. Geometries of triads wherein a pair of NCH bases share a single σ-hole of (a) FBr, (b) FI, and (c) FH, optimized with constraint of equal θ(N··X··N) angles. Distances in Å and angles in degs.
Molecules 26 00350 g003
Molecules 26 00350 g004 550
Figure 4. Contour plots of MEP of (a) FI, (b) FBr, (c) FH, (d) NH3 and (e) NCH. Blue line represents a vdW surface defined as isodensity ρ = 0.001 au.
Figure 4. Contour plots of MEP of (a) FI, (b) FBr, (c) FH, (d) NH3 and (e) NCH. Blue line represents a vdW surface defined as isodensity ρ = 0.001 au.
Molecules 26 00350 g004
Molecules 26 00350 g005 550
Figure 5. MEP of indicated molecule as function of θ defined as deviation from FX axis for acids, NCH axis for NCH, and C3 symmetry axis of NH3. Potentials were evaluated along a circle centered about the X or N atom, with a radius of 3 Å. The negative of the potential is shown for NCH and NH3.
Figure 5. MEP of indicated molecule as function of θ defined as deviation from FX axis for acids, NCH axis for NCH, and C3 symmetry axis of NH3. Potentials were evaluated along a circle centered about the X or N atom, with a radius of 3 Å. The negative of the potential is shown for NCH and NH3.
Molecules 26 00350 g005
Molecules 26 00350 g006 550
Figure 6. Contour plots of the σ *(FX) antibonding orbital (a) FI, (b) FBr, (c) FH, and the N lone pair orbitals of (d) NH3 and (e) NCH. Blue line represents a vdW surface defined as isodensity ρ = 0.001 au.
Figure 6. Contour plots of the σ *(FX) antibonding orbital (a) FI, (b) FBr, (c) FH, and the N lone pair orbitals of (d) NH3 and (e) NCH. Blue line represents a vdW surface defined as isodensity ρ = 0.001 au.
Molecules 26 00350 g006
Molecules 26 00350 g007 550
Figure 7. Amplitudes of the σ *(FX) antibonding orbital of FI, FBr, FH, and the N lone pair orbitals of NH3 and NC, plotted against deviation angle θ from FX axis for acids, NCH axis for NCH, and C3 symmetry axis of NH3. Wave functions were evaluated along a circle centered about the X or N atom, with a radius of 2 Å for all except 1.5 Å for FH.
Figure 7. Amplitudes of the σ *(FX) antibonding orbital of FI, FBr, FH, and the N lone pair orbitals of NH3 and NC, plotted against deviation angle θ from FX axis for acids, NCH axis for NCH, and C3 symmetry axis of NH3. Wave functions were evaluated along a circle centered about the X or N atom, with a radius of 2 Å for all except 1.5 Å for FH.
Molecules 26 00350 g007
Molecules 26 00350 g008 550
Figure 8. Total electron densities of indicated monomers, plotted against deviation angle θ from FX axis for acids, NCH axis for NCH, and C3 symmetry axis of NH3. Potentials were evaluated along a circle centered about the X or N atom, with a radius of 3 Å for all except 2 Å for FH.
Figure 8. Total electron densities of indicated monomers, plotted against deviation angle θ from FX axis for acids, NCH axis for NCH, and C3 symmetry axis of NH3. Potentials were evaluated along a circle centered about the X or N atom, with a radius of 3 Å for all except 2 Å for FH.
Molecules 26 00350 g008
Molecules 26 00350 g009 550
Figure 9. PAEM of indicated monomers, plotted against deviation angle θ from FX axis for acids, NCH axis for NCH, and C3 symmetry axis of NH3. Quantities were evaluated along a circle centered about the X or N atom, with a radius of 3 Å.
Figure 9. PAEM of indicated monomers, plotted against deviation angle θ from FX axis for acids, NCH axis for NCH, and C3 symmetry axis of NH3. Quantities were evaluated along a circle centered about the X or N atom, with a radius of 3 Å.
Molecules 26 00350 g009
Table
Table 1. Intermolecular distances (R, Å), interaction energies (kcal/mol), bond critical point density (ρBCP, 10−4 au), and NBO Nlp→σ *(FX) perturbation energy E(2) (kcal/mol) for complexes where one or two Lewis acids share a single lone pair on the N-base, and θ(X··N··X) angle (degs) for bifurcated bonds.
Table 1. Intermolecular distances (R, Å), interaction energies (kcal/mol), bond critical point density (ρBCP, 10−4 au), and NBO Nlp→σ *(FX) perturbation energy E(2) (kcal/mol) for complexes where one or two Lewis acids share a single lone pair on the N-base, and θ(X··N··X) angle (degs) for bifurcated bonds.
BaseAcidRθ-EintρBCPE(2)
NH31 FBr2.339 16.3056051.10
2 FBr2.76079.715.602078.95
1 FBr(def) a 9.4321310.71
NH31 FI2.493 19.1746543.61
2 FI2.73085.717.5126314.65
1 FI(def) 11.6326716.37
NH31 FH1.672 13.8853440.09
2 FH1.89566.213.8329013.85
1 FH(def) 9.7329015.44
NCH1 FBr2.629 7.2024712.13
2 FBr2.76688.210.511785.68
1 FBr(def) 6.071796.33
NCH1 FI2.621 10.2030120.09
2 FI2.86891.813.481736.66
1 FI(def) 8.051757.73
NCH1 FH1.875 7.3227112.29
2 FH1.97485.49.372055.89
1 FH(def) 6.042056.51
a single acid unit in geometry adopted within triad.
Table
Table 2. Intermolecular distances (R, Å), interaction energies (kcal/mol), bond critical point density (ρBCP, 10−4 au), and NBO Nlp→σ *(FX) perturbation energy E(2) (kcal/mol) for complexes where one or two NCH bases share a single σ-hole on the acid, and θ(N··X··N) angle (degs) for bifurcated bonds.
Table 2. Intermolecular distances (R, Å), interaction energies (kcal/mol), bond critical point density (ρBCP, 10−4 au), and NBO Nlp→σ *(FX) perturbation energy E(2) (kcal/mol) for complexes where one or two NCH bases share a single σ-hole on the acid, and θ(N··X··N) angle (degs) for bifurcated bonds.
AcidBaseRθ-EintρBCPE(2)
FBr1 NCH2.629 7.2024712.13
2 NCH3.09060.84.311111.86
1 NCH (def) a 3.331121.93
FI1 NCH2.621 10.2030120.09
2 NCH3.08159.06.831393.74
1 NCH (def) 5.171393.96
FH1 NCH1.875 7.3227112.29
2 NCH2.20790.36.461442.87
1 NCH (def) 4.661472.32
a single base unit in geometry adopted within triad.
Table
Table 3. Changes relative to monomer of internal bond length (Å), stretching frequency (cm−1), and NMR chemical shieldings (ppm) in dyads and triads a.
Table 3. Changes relative to monomer of internal bond length (Å), stretching frequency (cm−1), and NMR chemical shieldings (ppm) in dyads and triads a.
BaseAcidΔrΔν(FX)Δσ(X)Δσ(F)Δσ(N)
NH31 FBr0.061−107.21670.5−790.8−17.6
2 FBr0.012−27.0687.5−168.0−4.3
NH31 FI0.055−48.571.5−421.3−7.0
2 FI0.013+8.137.7−214.1+10.3
NH31 FH0.036−884.8−7.4−35.9−8.8
2 FH0.010−275.6−3.0−21.7−10.6
NCH1 FBr0.013−29.1680.6−183.8−21.6
2 FBr0.006−8.7484.5−122.0−12.8
NCH1 FI0.026−26.248.8−297.0−12.7
2 FI0.008+15.430.4−173.1−2.6
NCH1 FH0.008−235.5−2.2−6.9−19.4
2 FH0.005−130.1−1.6−11.9−21.6
1 NCHFBr0.013−29.1680.6−183.8−21.6
2NCH 0.013−20.578.2−66.5−45.2
1 NCHFI0.026−26.248.8−297.0−12.7
2NCH 0.020−20.716.3−136.1−41.1
1 NCHFH0.008−235.5−2.2−6.9−19.4
2NCH 0.009−182.7−1.3−16.3−41.4
a average of two values when there are two equivalent properties.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Scheiner, S. Comparison of Bifurcated Halogen with Hydrogen Bonds. Molecules 2021, 26, 350. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26020350

AMA Style

Scheiner S. Comparison of Bifurcated Halogen with Hydrogen Bonds. Molecules. 2021; 26(2):350. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26020350

Chicago/Turabian Style

Scheiner, Steve. 2021. "Comparison of Bifurcated Halogen with Hydrogen Bonds" Molecules 26, no. 2: 350. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26020350

Article Metrics

Back to TopTop