Na Liu1
Qingzhong Li
Sean A. C. McDowell- 1The Laboratory of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Yantai University, Yantai, China
- 2Department of Biological and Chemical Sciences, The University of the West Indies, Cave Hill Campus, Cave Hill, Barbados
The pnicogen, chalcogen, and halogen bonds between 6-OXF2-fulvene (X = As, Sb, Se, Te, Br, and I) and three nitrogen-containing bases (FCN, HCN, and NH3) are compared. For each nitrogen base, the halogen bond is strongest, followed by the pnicogen bond, and the chalcogen bond is weakest. For each type of bond, the binding increases in the FCN < HCN < NH3 pattern. Both FCN and HCN engage in a bond with comparable strengths and the interaction energies of most bonds are < −6 kcal/mol. However, the strongest base NH3 forms a much more stable complex, particularly for the halogen bond with the interaction energy going up to −18 kcal/mol. For the same type of interaction, its strength increases as the mass of the central X atom increases. These bonds are different in strength, but all of them are dominated by the electrostatic interaction, with the polarization contribution important for the stronger interaction. The presence of these bonds changes the geometries of 6-OXF2-fulvene, particularly for the halogen bond formed by NH3, where the F-X-F arrangement is almost vertical to the fulvene ring.
Introduction
Intermolecular interactions can regulate many chemical and biological processes (Oshovsky et al., 2007; Schneider, 2009; Zayed et al., 2010). Insight into these interactions is helpful in promoting development of supramolecular chemistry (Smith, 2005; Uhlenheuer et al., 2010; Bauzá et al., 2014), materials science (Müller-Dethlefs and Hobza, 2000; Vickaryous et al., 2004), and the rational design of new drugs and biochemistry (Xu et al., 2011; Lu et al., 2012). Among intermolecular interactions, the hydrogen bond (HB) is still widely investigated because of its importance in aqueous and biological systems. Of course, interest in other intermolecular interactions has been growing rapidly. These new types of intermolecular interactions were named for the type of atom that replaces the bridging proton in the HB, and these atoms include aerogen (Bauzá and Frontera, 2015a,b,c), halogen (Desiraju et al., 2013; Cavallo et al., 2016; Terraneo and Resnati, 2017), chalcogen (Azofra et al., 2015; Nziko and Scheiner, 2015; Nayak et al., 2017), pnicogen (Bauzá et al., 2016; Joy et al., 2016; Sánchez-Sanz et al., 2016), and tetrel (Grabowski, 2014a; Mani and Arunan, 2014; Li et al., 2015) atoms. Scientists have tried to provide a consistent explanation for the origin of these different intermolecular interactions. By means of molecular electrostatic potential (MEP) analyses, it was found that the surfaces of these covalently-bonded atoms all have positive MEPs (holes). Thus, these holes produce an attractive force with an approaching nucleophile. Usually, this hole is called σ-hole if this positive MEP region is located at the end of the covalent bond (Murray et al., 2007). On the other hand, the stability of all these bonds depends in part upon charge transfer from the electron-donor atom into an antibonding orbital of the acceptor (Scheiner, 2013a).
In studying these new intermolecular interactions, it is natural to compare them with HBs. Besides the consistency in the formation mechanism, they have similar applications in crystal materials (Gilday et al., 2015; Wang et al., 2016; Scilabra et al., 2019), molecular recognition (Cavallo et al., 2010; Ariga et al., 2012; Wenzel et al., 2012), chemical reactions (Walter et al., 2011; Liu et al., 2016; Wonner et al., 2017), and biological systems (Parisini et al., 2011; Lange et al., 2015; García-LLinás et al., 2017). For example, Scheiner compared the ability of the HB, halogen bond (XB), and tetrel bond (TB) in recognizing halides and found that TBs may be an effective halide receptor which displays good selectivity for both F− and Cl− anions (Scheiner, 2017). A systematic comparison of the angular geometries for B···ClF XB and B···HCl HB, as B is varied, was performed by Legon and parallels between angular geometries were observed for both XB and HB; the author claimed that the empirical rules for predicting angular geometries of HB can be extended to XB (Legon, 2008). Scheiner made a detailed comparison of the pnicogen bond (ZB) with the chalcogen bond (ChB), XB, and HB with regards to energetics, geometries, electron density shifts, and energy decomposition, and concluded that ZB has certain parallels with other bonds (Scheiner, 2013b). On the other hand, some differences were found among these bonds. Because pnicogen, chalcogen, and halogen atoms are strongly anisotropic, ZB, ChB, and XB are more sensitive to angular distortions than HB (Adhikari and Scheiner, 2012). Unlike HB and XB, the tetrel and pnicogen donor molecules are often distorted in order to accommodate the approaching base (Zierkiewicz et al., 2018).
The applications of intermolecular interactions are largely based upon their strengths; consequently, more attention was focused on their strengths and the competition between these bonds (Alkorta et al., 2008; Li et al., 2010; Solimannejad et al., 2012; An et al., 2013; Nagels et al., 2014; Guo et al., 2015; Dong et al., 2018; Lee et al., 2019). Hypohalous acids engage in a stronger HB, but a weaker XB, with various bases such as formaldehyde (Li et al., 2010), formamidine (An et al., 2013), NH3 and HCN (Alkorta et al., 2008). Dong et al. (2018) assessed the relative strengths of TB, ZB, ChB, and XB, which are represented by the third-row atoms Ge, As, Se, and Br, respectively. The H atoms around these central atoms were substituted firstly by a methyl group and then by F substitution in various locations (Dong et al., 2018). In the context of unsubstituted acids, the strengths of these bonds vary in the order ChB > XB > ZB ≈ TB, while F substitution at a position directly opposite the base results in a XB > ChB > ZB > TB ordering (Dong et al., 2018). In comparing ZB, ChB, and XB, the lone pairs around the central atoms are respectively three, two, and one. It is known that substituents have an important effect on the strength of these non-covalent bonds (Geboes et al., 2017). It is therefore necessary to bind equivalent atoms/groups to these pnicogen, chalcogen, and halogen atoms when comparing their strengths.
Fulvene is a non-aromatic molecule, but its aromaticity can be realized by substituents exo (6-position) to the five-membered ring (Krygowski et al., 2010). Fulvene and its derivatives, and particularly this 6-position substitution, have been studied due to their function in synthesis (Stone and Little, 1984; Strohfeldt and Tacke, 2008; Peloquin et al., 2012). The acidity of 6-OH-fulvene can be greatly increased by strong electron-withdrawing groups adjoined to the fulvene ring (Maksić and Vianello, 2004), which prompted us to explore the intermolecular interactions of the substituent at the 6-position of fulvene (Hou et al., 2019a,b). The presence of the fulvene skeleton strengthens the carbon bonding between 6-OCH3-fulvene and NH3 since its interaction energy is double that of the complex of methanol and NH3 (Hou et al., 2019a). The relative strengths of various bonds are also related to the nature and basicity of the base. Thioformaldehyde can reduce the difference in the binding strength between XB and HB in the HOBr complex (Li et al., 2011). If TeH2 acts as the electron donor, the XB is more favorable than the HB in complexes with 6-OX-fulvene (X = H, Cl, Br, I) (Hou et al., 2019b). This shows that introduction of a fulvene skeleton to an atom/group may cause an interesting result when this atom/group participates in an intermolecular interaction.
In this work, 6-OXF2-fulvene (X = As, Sb, Se, Te, Be, I) was the Lewis acid chosen to compare the strengths of ZB, ChB, and XB. All pnicogen, chalcogen, halogen atoms have the same number of substituents in their molecules. Moreover, these substituents are also the same for different X atoms. Thus, this molecule allows for a reliable comparison of these interactions since only X is varied. If X is a halogen atom, the corresponding 6-OXF2-fulvene molecule contains a hypervalent halogen. Hypervalent halogens have been shown to form a XB with various electron donors (Grabowski, 2014b; 2017; Kirshenboim and Kozuch, 2016). For X = Se and Te, there is similar bonding environment (Nordheider et al., 2015). Although 6-OXF2-fulvene (X = Se and Te) is neutral, it is a doublet, which is different from the halogen and pnicogen analogs. The electron donors selected range in electron-donating ability from FCN to HCN to NH3. The electron donors PH3 and AsH3 are compared with NH3. These complexes are also compared with the ZB, ChB, and XB of more conventional molecules found in the literature. We explore their similarities and differences in the present study by means of natural bond orbital (NBO), atoms in molecules (AIM), MEP, and energy decomposition (ED) analyses.
Computational Methods
All calculations were performed using the Gaussian 09 suite of programs (Frisch et al., 2009). Geometries were optimized at the MP2 computational level with the aug-cc-pVDZ basis set for all atoms except I, Sb, and Te, for which the aug-cc-pVDZ-PP basis set with relativistic corrections was adopted. Vibrational frequency calculations at the same level confirmed that the structures obtained correspond to energetic minima. The interaction energy was calculated by the supermolecular method involving the energies of the monomers at the geometries they adopt within the complex. This quantity was corrected for the basis set superposition error (BSSE) by the counterpoise protocol proposed by Boys and Bernardi (1970). The interaction energy was also decomposed into five terms including electrostatic, exchange, repulsion, polarization, and dispersion energies using the GAMESS program (Schmidt et al., 1993) and localized molecular orbital-energy decomposition analysis (LMOEDA) method (Su and Li, 2009) at the MP2/aug-cc-pVDZ(PP) level. Using the natural bond orbital (NBO) method (Reed et al., 1988) implemented in the Gaussian 09 program, the charge transfer and second-order perturbation energies were obtained. The AIM2000 package (Bader, 1990) was used to assess the topological parameters at each bond critical point (BCP), including the electron density, its Laplacian, and the energy density. Molecular electrostatic potentials (MEPs) were calculated on the 0.001 au isodensity surface at the MP2/aug-cc-pVDZ level using the WFA-SAS program (Bulat et al., 2010).
Results and Discussion
Structures and MEPs of Monomers
Figure 1 shows the structures of six Lewis acid molecules. This figure shows that the structures of the Lewis acid molecules formed by different X atoms are not the same. For the As and Sb monomers, the two F atoms are oriented toward (and are about the same distance away from) the H atom of the -C=C-H group. This indicates that the molecular structures of the As and Sb monomers possess Cs symmetry. The F atoms of both molecules are negatively charged, so they attract the H atoms. The F···H distance is shorter in the As monomer, which suggests that the attraction is stronger in this molecule. The Cs symmetry can also be measured by the C=C-O-X dihedral angle (Table 1). For As and Sb monomers, this dihedral angle is 180°, indicating that the X atom is coplanar with the plane of the fulvene ring.
Figure 1. The structures of monomers and distances in Å.
Table 1. Angles of C-O-X (α) and F-X-F (β) as well as C=C-O-X dihedral angle (θ) in the monomers, all in deg.
For the Se and Te monomers, one of the F atoms is also close to the H atom of the -C=C-H group, and its separation (~2.1 Å) is much shorter than the corresponding separation in the As and Sb monomers. This value is smaller than the sum of the van der Waals (vdW) radius of the corresponding atoms (2.45 Å), suggesting that the F···H interaction is not weak. This F···H interaction makes the F-X-F arrangement non-linear, and the corresponding F-X-F angle is smaller in the Te monomer. The C=C-O-X dihedral angle shows that the Te atom is coplanar with the fulvene skeleton, while the Se atom deviates slightly from this plane.
For the Br and I monomers, the halogen atom is completely turned away from the plane of the fulvene skeleton since the C=C-O-X dihedral angle is about 130°. The Br atom deviates more than the I atom, with a smaller C=C-O-X dihedral angle for the former. Similarly, the F-X-F group is not arranged in a straight line and the F-I-F arrangement is even more non-linear. The F···H interaction also exists in these two molecules, but it is weaker than it is in the Se and Te monomers.
For 6-OXF2-fulvene, the C-O-X angle increases slightly when X goes from a Group V to VI to VII atom, but it decreases by about 10° only in the Br monomer. The C=C-O-X dihedral angle is significantly reduced in the molecule with the Group VII atom, especially in the Br monomer. The F···H interaction is strongest for the Group VI atom and weakest for the Group V atom.
Figure 2 shows the MEP diagrams of all monomers with their MEP values. In 6-OXF2-fulvene, there is a positive MEP region (σ-hole) at the O-X bond end, and the F atom has a negative MEP due to its large electronegativity. For the same group, the σ-hole increases as the atom size increases (going down the group). For the same period, the σ-hole first decreases, then increases, with the maximum appearing for the Group V atom. Atomic radius decreases and electronegativity increases going from left to right along a particular period. A Group VI atom is difficult to polarize and thus has a smaller σ-hole. The F···H interaction for the Group VI and VII atoms reduces the F electron-withdrawing ability and thus enlarges the σ-hole. The Group VII atom has a weaker F···H interaction and thus a smaller σ-hole.
Figure 2. MEP maps of monomers. Color ranges are: red, larger than 0.02; yellow, between 0.02 and 0; green, between 0 and −0.002; blue, <-0.002. All are in a.u.
The last two rows in Figure 2 are the MEP diagrams of the Lewis base molecules. For the nitrogen-containing bases, HCN and NH3 are often used as electron donors in studying non-covalent interactions, and HCN is less basic than NH3. To further reduce its electron-donating ability, the H atom of HCN is replaced by an F atom. The negative MEP on the N atom confirms their relative basicities. For different types of ZH3 Lewis bases, as Z increases, the negative MEP gradually decreases, consistent with the electronegativity of Z.
Geometries of Complexes
Figure 3 shows the structures of the complexes. We first focus on the linear arrangement involving the N and X atoms, which can be measured by the O-X···N angle in Table 2. This angle ranges between 154 and 178°, with a larger angle indicating a more linear arrangement. For the same type of interaction with a heavier X atom, the O-X···N angle decreases in almost all complexes. The O-X···N angle in the NH3 complex is not the largest, contrary to our expectations: i.e., the stronger the interaction, the more linear the arrangement. The linear HCN and FCN molecules are obviously not on the same line as the O-X bond. For the ZB complex, HCN and FCN are generally biased toward the two F atoms of 6-OXF2-fulvene. For the ChB complex, HCN and FCN are generally in the same plane as the fulvene ring and biased to one F atom of 6-OXF2-fulvene, particularly for FCN. For the XB complex, neither HCN nor FCN is in the same plane as the fulvene ring because the O-X bond has deviated from the plane of the fulvene ring. In the two ZB complexes with NH3, the configurations involving the F atoms of -AsF2/-SbF2 and the H atoms of NH3 are different; in the former a crossed arrangement is evident, while in the latter an overlapping arrangement is evident. For the ChB and XB complexes with NH3, one N-H bond is coplanar with one of the X-F bonds.
Figure 3. The optimized structures of complexes, distances in Å.
Table 2. Angles of C-O-X (α), F-X-F (β), and O-X···N (γ), C=C-O-X dihedral angle (θ) in the complexes as well as their difference relative to the monomers (Δ), and ratio (χ) of the X···N distance relative to the sum of vdW radii of both atoms, all in deg.
When the X atom is fixed, the X···N distance shortens for the stronger nitrogen electron donor in most cases. However, for the ChB complex, the binding distance increases from FCN to HCN. For all interactions, when the nitrogen base goes from FCN to HCN, the binding distance hardly changes, rarely exceeding 0.03 Å; in going from HCN to NH3, the binding distance changes greatly, by up to 0.4 Å. The heavier X atom generally results in a longer X···N distance, however, for ZB, the Sb···N distance is shorter than the As···N distance.
Complex formation affects the structure of 6-OXF2-fulvene. For most complexes formed by the weak bases HCN and FCN, the C-O-X angle is almost unchanged, and it is reduced by ~8° only in the I-NCF and I-NCH complexes where the O-I bond deviates from the fulvene ring. The change of the C-O-X angle can be observed in the NH3 complex although its change is not consistent. Similarly, the largest reduction of the C-O-X angle is found in I-NH3.
In most complexes, the F-X-F angle hardly changes, especially for the weak bases HCN and FCN. Only in some complexes formed by NH3, this angle has an observed change, such as a reduction of 6° in Se-NH3. The configuration of the -AsF2 and -SbF2 groups in the complex are almost the same as in the monomer since the C-O-X and F-X-F angles, as well as the C=C-O-X dihedral angle, show little change. In the ZB complex, the F···H distance is shortened and its shortening is larger for the stronger Lewis base. In the ChB complex, the relative configuration of the -XF2 group is similar to that in the monomer since the change of the C=C-O-X dihedral angle, as well as the C-O-X and F-X-F angles, is small. In most ChB complexes, the F···H distance is also shortened, but it is increased in the NH3 complexes. Although the F-X-F angle varies little in the XB complex, it has a large change with respect to the fulvene skeleton since the change in the C=C-O-X dihedral angle is substantial. This dihedral angle has a larger change for the stronger nitrogen base. In the Br···NH3 and I···NH3 complexes, the F-X-F group is perpendicular to the plane containing the fulvene ring. Similarly, the F···H distance in most XB complexes has a similar change to that in the ChB complexes.
Interaction Energies of Complexes
Table 3 presents the interaction energy (Eint), binding energy (Eb), and deformation energy (DE) for all complexes. The binding energy is the difference in the energy of the complex relative to the sum of the energies for the isolated monomers (in their optimized geometries). In general, the trend for Eb is the same as for Eint, but Eb is smaller than Eint, and their difference is DE. The FCN and HCN complexes have a small DE, not exceeding 0.5 kcal/mol due to the weak interaction. The NH3 complex shows a large DE, especially for XB where DE is equivalent to 30% of Eint, consistent with the structural change described above. Such a large contribution of DE to Eint is seldom found for XBs with conventional halogen donors.
Table 3. Interaction energy (Eint), binding energy (Eb), and deformation energy (DE), all in kcal/mol.
Figure 4 plots the variation of the interaction energy for different systems. For each type of interaction, Eint gradually increases for the stronger nitrogen base, which is consistent with the negative MEP on the N atom. From FCN to HCN, the increase of Eint is small (not exceeding 1 kcal/mol), while from HCN to NH3, the increase of Eint is large, up to 11 kcal/mol. The Eint of the NH3 complex is increased by 65–185% relative to the HCN analog. This percentage increase is smallest for ChB and largest for XB. For the same type of interaction, Eint increases with the increase of X atomic mass, consistent with the σ-hole on the X atom. Regardless of the Lewis base, the interaction becomes stronger in the order ChB < ZB < XB, with some inconsistencies regarding the magnitude of the σ-hole on the X atom, due mainly to large structural changes in the halogen-containing molecule. This shows that the strength of the non-covalent bond is determined not only by the electrostatic potential, but also by other factors.
Figure 4. Plots of the interaction energy (Eint) for different bases.
In Figure 4, the interaction energy for the complexes of PH3/AsH3 is also plotted. The interactions are weaker for the heavier pnicogen atom, excluding Te-NH3, consistent with the negative MEP on the pnicogen atom. We compare the interaction strength of FCN/HCN with PH3/AsH3. The former has the more negative MEP, thus it forms a stronger ZB with As and ChB with Se. However, the opposite is found for the XB with Br. This abnormal result is also found for all three types of interactions when X = Sb, Te, and I, and is pronounced since the interaction energy with PH3/AsH3 is about twice as large as that for FCN/HCN.
In most complexes, the interaction energy is a net result of more than one interaction. To estimate the change of X···N interaction strength, the ratio of X···N distance relative to the sum of the van der Waals radii of both atoms is calculated (Table 2). In most cases, the smaller ratio corresponds to the larger interaction energy and thus a strong X···N interaction. However, some exceptions are present. For example, both Br-NCH and I-NCH have an equivalent ratio, but their interaction energies show an obvious difference.
There are some studies available comparing the strengths of ZB, ChB, and XB, which allow us to compare our results with those from previous studies. Dong and coauthors performed a detailed comparison for the XB, ChB, ZB, and tetrel bond (TB) in complexes of BrH, SeH2, AsH3, and GeH4 with NH3, respectively. The ChB is the strongest, followed by XB, ZB, and TB. The replacement of the atom opposite the base by F changes the order to XB > ChB > ZB > TB (Dong et al., 2018). Both orderings are different from our results. This indicates that the number and type of substituent adjoined to the acidic center has an important effect on the strength of the interaction.
Origin of Interactions
Figure 5 shows the AIM diagrams for the three types of complexes formed by NCF; the types of complexes for the other bases are similar. For each type of interaction, an intermolecular bond critical point (BCP) is found along the X···N path, confirming the formation of the corresponding non-covalent interaction. There is also a F···H BCP in the ChB complex, indicating the existence of a F···H interaction. This BCP is not found in the ZB and XB complexes since the F···H interaction is very weak. Actually, such weak interaction can be confirmed by the green region between the F and H atoms in the NCI analyses (Figure 6).
Figure 5. AIM diagrams for three different types of complexes.
Figure 6. NCI diagrams for ZB and XB types of complexes.
The topological data including electron density (ρ), Laplacian (∇2ρ), and total energy density (H) at the X···N BCP are collected in Table 4. For the complexes with the two weaker bases NCF and NCH, both ∇2ρ and H are positive, which indicates that the corresponding non-covalent bond belongs to a completely closed shell interaction, consistent with its relatively weak nature. For the complexes with NH3, except Se-NH3, ∇2ρ is still positive but H becomes negative, which indicates that the corresponding interaction has a partially covalent character (Arnold and Oldfield, 2000). As the base goes from FCN to HCN to NH3, the positive value for H gradually decreases and even becomes negative, while ∇2ρ generally becomes larger. For the same type of interaction with a heavier X atom and a fixed base, the positive H decreases, while the negative H becomes more negative. Therefore, ∇2ρ not only characterizes the type of interaction (by its sign) but can also gauge the interaction strength by its magnitude.
Table 4. Electron density (ρ), Laplacian (Δ2ρ), and total energy density (H) in the complexes, all in a.u.
For the same type of interaction, the electron density is almost the same in the complexes formed by FCN and HCN. When the nitrogen base changes from HCN to NH3, the electron density significantly increases and the maximal increase is 0.038 a.u in Br-NH3. For the different types of interactions, the electron density becomes larger in the order XB > ZB> ChB, although the intermolecular BCP is different.
Table 5 presents charge transfer (CT) and second-order perturbation energy data for the three types of complexes. For the same type of interaction, the CT value for the complex formed by the two weaker bases FCN and HCN is almost the same. When the Lewis base is NH3, CT increases significantly and this increase varies with the type of complex. The smallest increase is found for the ChB complex, while the largest is for the XB complex. For a fixed base, CT is smallest for ChB, while the relative CT for XB and ZB, is less straightforward but even so, XB has the largest CT for NH3. For the same type of interaction with the heavier X atom, CT also increases, but an unusual result occurs for Br-NH3, which has the largest CT (0.142e). Generally, there is correlation between CT and Eint.
Table 5. Charge transfer (CT, e) and second-order perturbation energy (E2, kcal/mol) in the complexes.
For XB, the main orbital interaction is LpN → BD*O-X, where LpN represents the lone pair orbital on the N atom and BD*O-X is the anti-bonding orbital of O-X bond. For ZB and ChB, the main orbital interaction is also LpN → BD*O-X when X = As and Se, but it is changed to LpN → Lp*X (Lp*X the empty orbital of the X atom) when X = Sb and Te. If the base is stronger, the corresponding orbital interaction is also stronger except the ChB formed by HCN. For the same type of interaction with the heavier X atom, the perturbation energy is larger excluding Br-NH3. When NH3 acts as a base, the orbital interaction is strongest in XB and weakest in ChB. When the Lewis base is FCN and HCN, ChB also has the weakest orbital interaction, while the relative strength of orbital interaction between XB and ZB is a little complicated.
To assess the similarities and differences between ZB, ChB and XB, the interaction energy was decomposed into electrostatic (Eele), exchange (Eex), repulsion (Erep), polarization (Epol) and dispersion energies (Edisp), as shown in Table 6. The three attractive terms (Eele, Epol, and Edisp) in the NH3 complexes are plotted in Figure 7. The pattern is basically the same for the other two Lewis bases. Each term has a similar trend with the interaction energy. Whether XB or ZB and ChB, the electrostatic energy is dominant, indicating that the three types of interactions are electrostatic in nature. When the base is FCN or HCN, both Epol and Edisp are comparable in magnitude. When the base is NH3, Epol is larger than Edisp. This indicates that the increase in the interaction energy is partly due to the increase in the polarization energy. Each term becomes more negative going from ChB to ZB to XB. For some complexes, the interaction energy is not consistent with the σ-hole MEP on the X atom, but the electrostatic interaction is still dominant in all complexes. This inconsistency is mainly attributable to the structural change of 6-OXF2-fulvene in the complex relative to the monomer.
Table 6. Electrostatic (Eele), exchange (Eex), repulsion (Erep), polarization (Epol), and dispersion energies (Edisp) in the complexes, all in kcal/mol.
Figure 7. Three attractive energy components.
Conclusions
Three different types of complexes formed between 6-OXF2-fulvene (X = As, Sb, Se, Te, Br, and I) and the three N bases FCN, HCN and NH3 were studied. The three types of interactions operative in the complexes have a predominantly electrostatic nature with differing strengths. The relative ordering for these interactions is halogen bond > pnicogen bond > chalcogen bond. Each type of interaction becomes stronger when the N base varies from FCN to HCN to NH3, consistent with the negative MEP on the N atom. The strength of each type of interaction is also consistent with electron density, charge transfer and orbital interaction analyses. The heavier X atom engages in a stronger interaction than its lighter analog, and the relative degree of enhancement depends on the type of interaction. For the stronger interactions, the polarization contribution becomes more dominant, as confirmed by the negative energy density at the corresponding intermolecular BCP, and this relatively large polarization causes a structural change of 6-OXF2-fulvene, particularly for X = halogen.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author Contributions
QL: conceptualization, software, resources, data curation, supervision, project administration, and funding acquisition. NL: methodology, formal analysis, investigation, writing—original draft preparation, and visualization. NL and QL: validation. QL and SM: writing—review and editing. All authors: contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (21573188).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2020.608486/full#supplementary-material
References
Adhikari, U., and Scheiner, S. (2012). Sensitivity of pnicogen, chalcogen, halogen and H-bonds to angular distortions. Chem. Phys. Lett. 532, 31–35. doi: 10.1016/j.cplett.2012.02.064
Alkorta, I., Blanco, F., Solimannejad, M., and Elguero, J. (2008). Competition of hydrogen bonds and halogen bonds in complexes of hypohalous acids with nitrogenated bases. J. Phys. Chem. A 112, 10856–10863. doi: 10.1021/jp806101t
An, X., Zhuo, H., Wang, Y., and Li, Q. (2013). Competition between hydrogen bonds and halogen bonds in complexes of formamidine and hypohalous acids. J. Mol. Model. 19, 4529–4535. doi: 10.1007/s00894-013-1969-7
Ariga, K., Ito, H., Hill, J. P., and Tsukube, H. (2012). Molecular recognition: from solution science to nano/materials technology. Chem. Soc. Rev. 41, 5800–5835. doi: 10.1039/c2cs35162e
Arnold, W. D., and Oldfield, E. (2000). The chemical nature of hydrogen bonding in proteins via NMR: J-couplings, chemical shifts, and AIM theory. J. Am. Chem. Soc. 122, 12835–12841. doi: 10.1021/ja0025705
Azofra, L. M., Alkorta, I., and Scheiner, S. (2015). Chalcogen bonds in complexes of SOXY (X, Y= F, Cl) with nitrogen bases. J. Phys. Chem. A 119, 535–541. doi: 10.1021/jp511828h
Bauzá, A., and Frontera, A. (2015a). Aerogen bonding interaction: a new supramolecular force? Angew. Chem. Int. Ed. 54, 7340–7343. doi: 10.1002/anie.201502571
Bauzá, A., and Frontera, A. (2015b). Theoretical study on the dual behavior of XeO3 and XeF4 toward aromatic rings: lone pair–π versus aerogen–π interactions. Chemphyschem 16, 3625–3630. doi: 10.1002/cphc.201500757
Bauzá, A., and Frontera, A. (2015c). π-Hole aerogen bonding interactions. Phys. Chem. Chem. Phys. 17, 24748–24753. doi: 10.1039/C5CP03272E
Bauzá, A., Mooibroek, T. J., and Frontera, A. (2014). Non-covalent sp3 carbon bonding with ArCF3 is analogous to CH–π interactions. Chem. Commun. 50, 12626–12629. doi: 10.1039/C4CC05602G
Bauzá, A., Mooibroek, T. J., and Frontera, A. (2016). σ-Hole opposite to a lone pair: unconventional pnicogen bonding interactions between ZF3 (Z = N, P, As, and Sb) compounds and several donors. Chemphyschem 17, 1608–1614. doi: 10.1002/cphc.201600073
Boys, S. F., and Bernardi, F. (1970). The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors. Mol. Phys. 19, 553–566. doi: 10.1080/00268977000101561
Bulat, F. A., Toro-Labbe, A., Brinck, T., Murray, J. S., and Politzer, P. (2010). Quantitative analysis of molecular surfaces: areas, volumes, electrostatic potentials and average local ionization energies. J. Mol. Model. 16, 1679–1691. doi: 10.1007/s00894-010-0692-x
Cavallo, G., Metrangolo, P., Milani, R., Pilati, T., Priimagi, A., Resnati, G., et al. (2016). The halogen bond. Chem. Rev. 116, 2478–2601. doi: 10.1021/acs.chemrev.5b00484
Cavallo, G., Metrangolo, P., Pilati, T., Resnati, G., Sansotera, M., and Terraneo, G. (2010). Halogen bonding: a general route in anion recognition and coordination. Chem. Soc. Rev. 39, 3772–3783. doi: 10.1039/b926232f
Desiraju, G. R., Ho, P. S., Kloo, L., Legon, A. C., Marquardt, R., Metrangolo, P., et al. (2013). Definition of the halogen bond (IUPAC Recommendations 2013). Pure Appl. Chem. 85, 1711–1713. doi: 10.1351/PAC-REC-12-05-10
Dong, W., Li, Q., and Scheiner, S. (2018). Comparative strengths of tetrel, pnicogen, chalcogen, and halogen bonds and contributing factors. Molecules 23:1681. doi: 10.3390/molecules23071681
Frisch, M. J., Trucks, G. W., Schlegel, H. B., Scuseria, G. E., Robb, M. A., Cheeseman, J. R., et al. (2009). Gaussian 09, Revision A. 02. Wallingford, CT: Gaussian, Inc.
García-LLinás, X., Bauzá, A., Seth, S. K., and Frontera, A. (2017). Importance of R–CF3···O tetrel bonding interactions in biological systems. J. Phys. Chem. A 121, 5371–5376. doi: 10.1021/acs.jpca.7b06052
Geboes, Y., De Proft, F., and Herrebout, W. A. (2017). Effect of fluorination on the competition of halogen bonding and hydrogen bonding: complexes of fluoroiodomethane with dimethyl ether and trimethylamine. J. Phys. Chem. A 121, 4180–4188. doi: 10.1021/acs.jpca.7b03206
Gilday, L. C., Robinson, S. W., Barendt, T. A., Langton, M. J., Mullaney, B. R., and Beer, P. D. (2015). Halogen bonding in supramolecular chemistry. Chem. Rev. 115, 7118–7195. doi: 10.1021/cr500674c
Grabowski, S. J. (2014b). Halogen bond with the multivalent halogen acting as the Lewis acid center. Chem. Phys. Lett. 605, 131–136. doi: 10.1016/j.cplett.2014.05.029
Grabowski, S. J. (2014a). Tetrel bond–σ-hole bond as a preliminary stage of the SN2 reaction. Phys. Chem. Chem. Phys. 16, 1824–1834. doi: 10.1039/C3CP53369G
Grabowski, S. J. (2017). New type of halogen bond: multivalent halogen interacting with π-and σ-electrons. Molecules 22:2150. doi: 10.3390/molecules22122150
Guo, X., Liu, Y. W., Li, Q. Z., Li, W. Z., and Cheng, J. B. (2015). Competition and cooperativity between tetrel bond and chalcogen bond in complexes involving F2CX (X = Se and Te). Chem. Phys. Lett. 620, 7–12. doi: 10.1016/j.cplett.2014.12.015
Hou, M. C., Li, Q. Z., and Scheiner, S. (2019b). Comparison between hydrogen and halogen bonds in complexes of 6-OX-fulvene with pnicogen and chalcogen electron donors. Chemphyschem 20, 1978–1984. doi: 10.1002/cphc.201900340
Hou, M. C., Yang, S. B., Li, Q. Z., Cheng, J. B., Li, H. B., and Liu, S. F. (2019a). Tetrel bond between 6-OTX3-fulvene and NH3: substituents and aromaticity. Molecules 24:10. doi: 10.3390/molecules24010010
Joy, J., Jose, A., and Jemmis, E. D. (2016). Continuum in the X-Z—Y weak bonds: Z = main group elements. J. Comput. Chem. 37, 270–279. doi: 10.1002/jcc.24036
Kirshenboim, O., and Kozuch, S. (2016). How to twist, split and warp a σ-hole with hypervalent halogens. J. Phys. Chem. A 120, 9431–9445. doi: 10.1021/acs.jpca.6b07894
Krygowski, T. M., Oziminski, W. P., Palusiak, M., Fowler, P. W., and McKenzie, A. D. (2010). Aromaticity of substituted fulvene derivatives: substituent-dependent ring currents. Phys. Chem. Chem. Phys. 12, 10740–10745. doi: 10.1039/c003686b
Lange, A., Günther, M., Büttner, F. M., Zimmermann, M. O., Heidrich, J., Hennig, S., et al. (2015). Targeting the gatekeeper MET146 of C-Jun N-terminal kinase 3 (JNK3) induces a bivalent halogen/chalcogen bond. J. Am. Chem. Soc. 137, 14640–14652. doi: 10.1021/jacs.5b07090
Lee, L. M., Tsemperouli, M., Poblador-Bahamonde, A. I., Benz, S., Sakai, N., Sugihara, K., et al. (2019). Anion transport with pnictogen bonds in direct comparison with chalcogen and halogen bonds. J. Am. Chem. Soc. 141, 810–814. doi: 10.1021/jacs.8b12554
Legon, A. C. (2008). The interaction of dihalogens and hydrogen halides with Lewis bases in the gas phase: an experimental comparison of the halogen bond and the hydrogen bond. Halogen Bonding 126, 17–64. doi: 10.1007/430_2007_063
Li, Q. Z., Jing, B., Li, R., Liu, Z. B., Li, W. Z., Luan, F., et al. (2011). Some measures for making halogen bonds stronger than hydrogen bonds in H2CS–HOX (X = F, Cl, and Br) complexes. Phys. Chem. Chem. Phys. 13, 2266–2271. doi: 10.1039/C0CP01543A
Li, Q. Z., Xu, X. S., Liu, T., Jing, B., Li, W. Z., Cheng, J. Z., et al. (2010). Competition between hydrogen bond and halogen bond in complexes of formaldehyde with hypohalous acids. Phys. Chem. Chem. Phys. 12, 6837–6843. doi: 10.1039/b926355a
Li, Q. Z., Zhuo, H. Y., Li, H. B., Liu, Z. B., Li, W. Z., and Cheng, J. B. (2015). Tetrel–hydride interaction between XH3F (X = C, Si, Ge, Sn) and HM (M = Li, Na, BeH, MgH). J. Phys. Chem. A 119, 2217–2224. doi: 10.1021/jp503735u
Liu, M., Li, Q., Cheng, J., Li, W., and Li, H. B. (2016). Tetrel bond of pseudohalide anions with XH3F (X = C, Si, Ge, and Sn) and its role in SN2 reaction. J. Chem. Phys. 145:224310. doi: 10.1063/1.4971855
Lu, Y., Liu, Y., Xu, Z., Li, H., Liu, H., and Zhu, W. (2012). Halogen bonding for rational drug design and new drug discovery. Expert Opin. Drug Disc. 7, 375–383. doi: 10.1517/17460441.2012.678829
Maksić, Z. B., and Vianello, R. (2004). Tailoring of strong neutral organic superacids: DFT-B3LYP calculations on some fulvene derivatives. New J. Chem. 28, 843–846. doi: 10.1039/b316290g
Mani, D., and Arunan, E. (2014). The X–C···π (X = F, Cl, Br, Cn) carbon bond. J. Phys. Chem. A 118, 10081–10089. doi: 10.1021/jp507849g
Müller-Dethlefs, K., and Hobza, P. (2000). Noncovalent interactions: a challenge for experiment and theory. Chem. Rev. 100, 143–168. doi: 10.1021/cr9900331
Murray, J. S., Lane, P., and Politzer, P. (2007). A predicted new type of directional noncovalent interaction. Int. J. Quantum Chem. 107, 2286–2292. doi: 10.1002/qua.21352
Nagels, N., Geboes, Y., Pinter, B., De Proft, F., and Herrebout, W. A. (2014). Tuning the halogen/hydrogen bond competition: a spectroscopic and conceptual DFT study of some model complexes involving CHF2I. Chem. Eur. J. 20, 8433–8443. doi: 10.1002/chem.201402116
Nayak, S. K., Kumar, V., Murray, J. S., Politzer, P., Terraneo, G., Pilati, T., et al. (2017). Fluorination promotes chalcogen bonding in crystalline solids. CrystEngComm 19, 4955–4959. doi: 10.1039/C7CE01070B
Nordheider, A., Woollins, J. D., and Chivers, T. (2015). Organophosphorus–tellurium chemistry: from fundamentals to applications. Chem. Rev. 115, 10378–10406. doi: 10.1021/acs.chemrev.5b00279
Nziko, V. D. P. N., and Scheiner, S. (2015). Intramolecular S···O chalcogen bond as stabilizing factor in geometry of substituted phenyl-SF3 molecules. J. Org. Chem. 80, 2356–2363. doi: 10.1021/acs.joc.5b00012
Oshovsky, G. V., Reinhoudt, D. N., and Verboom, W. (2007). Supramolecular chemistry in water. Angew. Chem. Int. Ed. 46, 2366–2393. doi: 10.1002/anie.200602815
Parisini, E., Metrangolo, P., Pilati, T., Resnati, G., and Terraneo, G. (2011). Halogen bonding in halocarbon–protein complexes: a structural survey. Chem. Soc. Rev. 40, 2267–2278. doi: 10.1039/c0cs00177e
Peloquin, A. J., Stone, R. L., Avila, S. E., Rudico, E. R., Horn, C. B., Gardner, K. A., et al. (2012). Synthesis of 1, 3-diphenyl-6-alkyl/aryl-substituted fulvene chromophores: observation of π-π interactions in a 6-pyrene-substituted 1,3-diphenylfulvene. J. Org. Chem. 77, 6371–6376. doi: 10.1021/jo301101x
Reed, A. E., Curtiss, L. A., and Weinhold, F. (1988). Intermolecular interactions from a natural bond orbital, donor-acceptor viewpoint. Chem. Rev. 88, 899–926. doi: 10.1021/cr00088a005
Sánchez-Sanz, G., Trujillo, C., Alkorta, I., and Elguero, J. (2016). Modulating intramolecular P···N pnictogen interactions. Phys. Chem. Chem. Phys. 18, 9148–9160. doi: 10.1039/C6CP00227G
Scheiner, S. (2013b). Detailed comparison of the pnicogen bond with chalcogen, halogen, and hydrogen bonds. Int. J. Quantum Chem. 113, 1609–1620. doi: 10.1002/qua.24357
Scheiner, S. (2013a). The pnicogen bond: its relation to hydrogen, halogen, and other noncovalent bonds. Acc. Chem. Res. 46, 280–288. doi: 10.1021/ar3001316
Scheiner, S. (2017). Assembly of effective halide receptors from components. Comparing hydrogen, halogen, and tetrel bonds. J. Phys. Chem. A 121, 3606–3615. doi: 10.1021/acs.jpca.7b02305
Schmidt, M. W., Baldridge, K. K., Boatz, J. A., Elbert, S. T., and Montgomery, J. A. (1993). General atomic and molecular electronic structure system. J. Comput. Chem. 14, 1347–1363. doi: 10.1002/jcc.540141112
Schneider, H. J. (2009). Binding mechanisms in supramolecular complexes. Angew. Chem. Int. Ed. 48, 3924–3977. doi: 10.1002/anie.200802947
Scilabra, P., Terraneo, G., and Resnati, G. (2019). The chalcogen bond in crystalline solids: a world parallel to halogen bond. Acc. Chem. Res. 52, 1313–1324. doi: 10.1021/acs.accounts.9b00037
Smith, D. K. (2005). A supramolecular approach to medicinal chemistry: medicine beyond the molecule. J. Chem. Edu. 82, 393–400. doi: 10.1021/ed082p393
Solimannejad, M., Ramezani, V., Trujillo, C., Alkorta, I., Sánchez-Sanz, G., and Elguero, J. (2012). Competition and interplay between σ-hole and π-hole interactions: a computational study of 1:1 and 1:2 complexes of nitryl halides (O2NX) with ammonia. J. Phys. Chem. A 116, 5199–5206. doi: 10.1021/jp300540z
Stone, K. J., and Little, R. D. (1984). An exceptionally simple and efficient method for the preparation of a wide variety of fulvenes. J. Org. Chem. 49, 1849–1853. doi: 10.1021/jo00185a001
Strohfeldt, K., and Tacke, M. (2008). Bioorganometallic fulvene-derived titanocene anti-cancer drugs. Chem. Soc. Rev. 37, 1174–1187. doi: 10.1039/b707310k
Su, P., and Li, H. (2009). Energy decomposition analysis of covalent bonds and intermolecular interactions. J. Chem. Phys. 131:014102. doi: 10.1063/1.3159673
Terraneo, G., and Resnati, G. (2017). Bonding matters. Cryst. Growth Des. 17, 1439–1440. doi: 10.1021/acs.cgd.7b00309
Uhlenheuer, D. A., Petkau, K., and Brunsveld, L. (2010). Combining supramolecular chemistry with biology. Chem. Soc. Rev. 39, 2817–2826. doi: 10.1039/b820283b
Vickaryous, W. J., Herges, R., and Johnson, D. W. (2004). Arsenic–π interactions stabilize a self-assembled As2L3 supramolecular complex. Angew. Chem. Int. Ed. 116, 5955–5957. doi: 10.1002/ange.200461011
Walter, S. M., Kniep, F., Herdtweck, E., and Huber, S. M. (2011). Halogen-bond-induced activation of a carbon–heteroatom bond. Angew. Chem. Int. Ed. 50, 7187–7191. doi: 10.1002/anie.201101672
Wang, H., Wang, W., and Jin, W. J. (2016). σ-Hole bond vs π-hole bond: a comparison based on halogen bond. Chem. Rev. 116, 5072–5104. doi: 10.1021/acs.chemrev.5b00527
Wenzel, M., Hiscock, J. R., and Gale, P. A. (2012). Anion receptor chemistry: highlights from 2010. Chem. Soc. Rev. 41, 480–520. doi: 10.1039/C1CS15257B
Wonner, P., Vogel, L., Düser, M., Gomes, L., Kniep, F., Mallick, B., et al. (2017). Carbon–halogen bond activation by selenium-based chalcogen bonding. Angew. Chem. Int. Ed. 56, 12009–12012. doi: 10.1002/anie.201704816
Xu, Z., Liu, Z., Chen, T., Chen, T., Wang, Z., Tian, G., et al. (2011). Utilization of halogen bond in lead optimization: a case study of rational design of potent phosphodiesterase type 5 (PDE5) inhibitors. J. Med. Chem. 54, 5607–5611. doi: 10.1021/jm200644r
Zayed, J. M., Nouvel, N., Rauwald, U., and Scherman, O. A. (2010). Chemical complexity—supramolecular self-assembly of synthetic and biological building blocks in water. Chem. Soc. Rev. 39, 2806–2816. doi: 10.1039/b922348g
Keywords: halogen bond, pnicogen bond, NBO, AIM, chalcogen bond
Citation: Liu N, Li Q and McDowell SAC (2020) Reliable Comparison of Pnicogen, Chalcogen, and Halogen Bonds in Complexes of 6-OXF2-Fulvene (X = As, Sb, Se, Te, Be, I) With Three Electron Donors. Front. Chem. 8:608486. doi: 10.3389/fchem.2020.608486
Received: 20 September 2020; Accepted: 27 October 2020;
Published: 09 December 2020.
Edited by:
Soumyajit Roy, Indian Institute of Science Education and Research Kolkata, IndiaReviewed by:
Antonio Frontera, University of the Balearic Islands, SpainFangfang Pan, Central China Normal University, China
Copyright © 2020 Liu, Li and McDowell. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qingzhong Li, liqingzhong1990@sina.com; Sean A. C. McDowell, sean.mcdowell@cavehill.uwi.edu