Next Article in Journal
Ultrasonic-Assisted Efficient Extraction of Coumarins from Peucedanum decursivum (Miq.) Maxim Using Deep Eutectic Solvents Combined with an Enzyme Pretreatment
Next Article in Special Issue
A QCT View of the Interplay between Hydrogen Bonds and Aromaticity in Small CHON Derivatives
Previous Article in Journal
A Branched and Double Alpha-Gal-Bearing Synthetic Neoglycoprotein as a Biomarker for Chagas Disease
Previous Article in Special Issue
High-Density “Windowpane” Coordination Patterns of Water Clusters and Their NBO/NRT Characterization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Coexistence of the Carbene⋯H-D Hydrogen Bond and Other Accompanying Interactions in Forty Dimers of N-Heterocyclic-Carbenes (I, IMe2, IiPr2, ItBu2, IMes2, IDipp2, IAd2; I = imidazol-2-ylidene) and Some Fundamental Proton Donors (HF, HCN, H2O, MeOH, NH3)

by
Mirosław Jabłoński
Faculty of Chemistry, Nicolaus Copernicus University in Toruń, ul. Gagarina 7, 87-100 Toruń, Poland
Molecules 2022, 27(17), 5712; https://doi.org/10.3390/molecules27175712
Submission received: 26 July 2022 / Revised: 25 August 2022 / Accepted: 29 August 2022 / Published: 5 September 2022
(This article belongs to the Collection Hydrogen Bonds)

Abstract

:
The subject of research is forty dimers formed by imidazol-2-ylidene (I) or its derivative (IR 2 ) obtained by replacing the hydrogen atoms in both N-H bonds with larger important and popular substituents of increasing complexity (methyl = Me, iso-propyl = i Pr, tert-butyl = t Bu, phenyl = Ph, mesityl = Mes, 2,6-diisopropylphenyl = Dipp, 1-adamantyl = Ad) and fundamental proton donor (HD) molecules (HF, HCN, H 2 O, MeOH, NH 3 ). While the main goal is to characterize the generally dominant C⋯H-D hydrogen bond engaging a carbene carbon atom, an equally important issue is the often omitted analysis of the role of accompanying secondary interactions. Despite the often completely different binding possibilities of the considered carbenes, and especially HD molecules, several general trends are found. Namely, for a given carbene, the dissociation energy values of the IR 2 HD dimers increase in the following order: NH 3 < H 2 O < HCN ≤ MeOH ≪ HF. Importantly, it is found that, for a given HD molecule, IDipp 2 forms the strongest dimers. This is attributed to the multiplicity of various interactions accompanying the dominant C⋯H-D hydrogen bond. It is shown that substitution of hydrogen atoms in both N-H bonds of the imidazol-2-ylidene molecule by the investigated groups leads to stronger dimers with HF, HCN, H 2 O or MeOH. The presented results should contribute to increasing the knowledge about the carbene chemistry and the role of intermolecular interactions, including secondary ones.

Graphical Abstract

1. Introduction

In carbenes [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80], the carbon atom is s p 2 -hybridized. However, unlike the vast majority of organic compounds, in carbenes, it is merely divalent, consequently forming only one (R=C) or at most two (R 1 R 2 C) covalent bonds. Thus, two valence electrons remain available, giving the possibility of two spin states, singlet or triplet [5,6,8,70]. In the spin singlet state, both electrons are paired to form a lone electron pair. Additionally, the carbene carbon atom has an unfilled p-orbital perpendicular to the s p 2 hybrids. While the triplet state is also possible in some carbenes, a number of factors are known to increase singlet stability. The close presence of a strongly electronegative atom, such as, e.g., N, is important here, which stabilizes the singlet state through both the inductive σ -electron withdrawing effect and the π -electron charge donation from the lone electron pair of this atom to the vacant p-orbital of the carbene carbon atom [5,6,8,18,70]. In the case of cyclic carbenes, another factor stabilizing the singlet state is the π -electron delocalization within the entire ring, which may be related to its aromaticity, as is the case, for example, in the imidazol-2-ylidene molecule [38]. Another factor that stabilizes the singlet state of a carbene is a low value of the R 1 –C–R 2 angle [3,4,5]. Naturally, this condition is obviously present in cyclic carbenes. While all of the conditions mentioned so far exist in imidazol-2-ylidene, they need not necessarily occur simultaneously in other N-heterocyclic carbenes (NHCs) [61]. In their case, an important stabilizing factor is the size of the substituents attached to both nitrogen atoms, as larger substituents prevent carbene dimerization [1].
The presence of a lone electron pair on the carbene carbon atom makes carbenes good Lewis bases, showing strongly nucleophilic properties. Carbenes are mainly known for their association with transition metal atoms (see, e.g., refs. [11,30,31,36,55,63]), making them extremely useful compounds in organic, organometallic and materials chemistry as well as homogeneous catalysis [21,30,33,57,61,62,69]. In addition to bonds with metals, however, carbenes quite willingly also form other intermolecular connections, such as hydrogen bonds [9,10,16,19,38,51,56,58], lithium bonds [39,45,46,62,69], beryllium bonds [17,53,54,62,69,73,79], magnesium bonds [14,22,40,41,62,69,73,79], triel bonds [13,20,23,27,69], tetrel bonds [42,65,66,69], pnictogen bonds [47,48,50,69], chalcogen bonds [52,69], halogen bonds [12,15,43,59,69] (in particular to iodine [12,15]), and aerogen bonds [67]. Moreover, in addition to these possibilities resulting from the presence of a lone electron pair, the presence of an empty p orbital in singlet carbenes also gives them electrophilic properties that seem to be much less studied [7,24,29,68,71,74,75].
Considering the fundamental role of hydrogen bonds [81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96], it is somewhat surprising that hydrogen bonds involving carbenes are studied only very sporadically [9,10,16,19,38,51,56,58]. Obviously, due to considerable methodological limitations, the first theoretical reports (in 1983 by Pople et al. [9] and then in 1986 by Pople [10]) concerned a simple H 2 C⋯HF dimer. Then, in 1996, Alkorta and Elguero investigated dimers between H 2 C or F 2 C and a few simple proton donors [19]. More recently, Jabłoński and Palusiak studied hydrogen bonds between carbenes CF 2 , CCl 2 and imidazol-2-ylidene and such proton donors as H 2 O and HCF 3 [38]. In a wonderful extensive work on the application of theoretical methods in the study of carbene chemistry, Gerbig and Ley [56] also mentioned the imidazol-2-ylidene⋯HCF 3 dimer in which there is C⋯H-C hydrogen bonding and the accompanying N-H⋯F (both described earlier in ref. [38]). Samanta et al. [58] showed the possibility of the formation of either C⋯H-N or C⋯H-O hydrogen bonds between the simple heterocyclic derivative of 1,3-di(methyl)imidazol-2-ylidene (IMe 2 ) and MeNH 2 or, in particular, MeOH, respectively. In the latter case, the NHC acts as an esterification catalyst activating the alcohol molecule. On the other hand, from the experimental point of view, hydrogen bonding with carbene was first announced by Arduengo et al. [16] in 1995 with the report essentially relating to the C–H–C bridge in a bis(carbene)–proton complex formed by 1,3-di(2,4,6-trimethylphenyl)-imidazol-2-ylidene (IMes 2 ). Much later, in 2011, it was shown that IMes 2 and 1,3-di(2,6-diisopropylphenyl)imidazolidin-2-ylidene (SI i Pr 2 ) can form a C⋯H-O type hydrogen bond to 1-hydroxy-2,2,6,6-tetramethyl-piperidine (TEMPO-H) [51]. In addition to crystallographic and NMR studies, results of theoretical calculations were also reported.
As shown above, the theoretical studies for hydrogen bonding involving carbenes have generally been down to very small carbenes, at most imidazol-2-ylidene (I) or 1,3-di(methyl)imidazol-2-ylidene (IMe 2 ). However, there are no theoretical studies in which the imidazol-2-ylidene derivative would contain even larger substituents on both nitrogen atoms. Moreover, the influence of the presence of these substituents on the possible interactions accompanying the leading C⋯H-D hydrogen bond has not been investigated yet. This article aims to fill this gap. Namely, this article examines the hydrogen bonds between imidazol-2-ylidene (I) and its seven popular derivatives containing gradually more bulky substituents (methyl = Me, iso-propyl = i Pr, tert-butyl = t Bu, phenyl = Ph, mesityl = Mes, 2,6-diisopropylphenyl = Dipp, 1-adamantyl = Ad) on both nitrogen atoms and five fundamental proton donor molecules (HF, HCN, H 2 O, MeOH, NH 3 ). The possible combinations form forty dimers, which can be briefly designated as IR 2 HD, where R is one of the substituents mentioned previously. The general scheme of the dimers in question is shown in Figure 1.
It should be emphasized that the C⋯H-D hydrogen bonds studied here engage the carbene species in their singlet spin states and thus a lone electron pair on the carbene carbon atom (Figure 1) and are therefore considerably distant from hydrogen bonds involving radicals [97,98]. In addition to the description of C⋯HD hydrogen bonds, an equally important goal is to analyze the possibility of the emergence of various types of accompanying interactions and their impact on the structure of the obtained dimers. This concerns the important issue of the coexistence of various types of interactions in molecular systems. Recently, the author of the present paper has shown [80] that various types of secondary interactions have a significant effect on the mutual orientation of the ZnX 2 molecular plane relative to the imidazol-2-ylidene ring in various types of IR 2 carbenes. This result is important because the possible torsion of the planes is usually attributed to the steric effects resulting from the presence of bulky R substituents.

2. Theoretical Methods

Geometries of all the systems were fully optimized on the ω B97X-D/6-311++G(d,p) level of theory, i.e., using the ω B97X-D range-separated gradient- and dispersion-corrected hybrid exchange-correlation functional [99] of density functional theory [100,101,102] and the 6-311++G(d,p) basis set being of the triple-zeta type and containing both polarization and diffuse functions on all atoms [103,104,105,106,107]. It is worth noting that ω B97X-D was one of the best functionals out of 200 tested [108]. There were no imaginary frequencies showing that equilibrium structures were obtained each time. Both the geometry optimization and the frequency analysis were performed using the Gaussian 16 (Revision C.01) program [109].
Values of the electron density at the bond critical points (bcp) [110,111,112] of the C⋯H-D hydrogen bonds and other accompanying interactions were computed using the AIMAll program [113]. Indeed, the C⋯H-D hydrogen bond should also be investigated through quantum chemical topology, which, over the years, was shown to be an efficient theoretical approach [114,115].
The dissociation energy was calculated as the difference between the total energy of a dimer and the sum of total energies of isolated subsystems in their own fully optimized structures. Total energies were corrected for the zero-point vibrational energies (ZPVE). Dissociation energies are given as positive values. In order to determine the binding energy (in fact, the interaction energy [116]) of the individual interaction of interest, the formula suggested recently by Emamian et al. [116] was used:
E b [ kcal / mol ] = 223.08 · ρ bcp [ a . u . ] + 0.7423
This formula is based on the electron density value determined at the bond critical point of this interaction, which is the parameter that, as has been shown [116], is best correlated with binding energy among many wave function-based HB descriptors. As a result, Emamian et al. proposed using this equation for a quick estimate of the energy of hydrogen-bond-forming networks. It is worth emphasizing that the values of the dissociation energies and the (sum) of the binding energies determined by Equation (1) are not comparable, since the former is a global quantity relating to the entire dimer, while the latter is a local parameter relating to an individual interaction. Besides, the dissociation energy takes into account deformation energies and ZPVEs.

3. Results and Discussion

3.1. Hydrogen Bonds to Carbenes Observed in the Solid State

It should be noted that despite the high reactivity of carbene compounds, hydrogen bonds with the participation of the carbene carbon atom as a proton acceptor quite often occur in the crystals. Of course, due to the abundance of C-H bonds, the most common hydrogen bond of this type is C⋯H-C, but C⋯H-N or C⋯H-O can also be found in the Cambridge Structural Database [117]. Some more interesting examples are shown in Figure 2.
There are two C⋯H-C type hydrogen bonds in the crystal form of 1,3-diisopropyl-4,5-dimethylimidazol-2-ylidene (GUXJAK02) [118]. Both the dimer between 1,3-dimesityl-imidazol-2-ylidene and diphenylamine (MODVEG) [119] and between 1,3-bis[2,6- diisopropylphenyl]imidazol-2-ylidene and N 1 ,N 4 -diphenylbenzene-1,4-diamine (LAGYEB) [120] have one C⋯H-N hydrogen bond, while the dimer of N-methyl-2-(3-phenyl-2,3-dihydro-1H-imidazol-1-yl)ethanamine (AWIBOY) [121] has two. 1,3-dimesitylimidazol-2-ylidene in methanol (JAPDEK) [122], on the other hand, is the simplest example of the occurrence of a C⋯H-O type hydrogen bond and was reported as the first X-ray structure of a carbene–alkohol hydrogen-bonded complex. The -OH group donors are somewhat more complex in QIGPEE (i.e., the dimer of 1-(2-hydroxy-2-methylpropyl)-3-(2,4,6-trimethylphenyl)-imidazol-2-ylidene) [123], USINAM (i.e., the dimer between 1,3-dimesitylimidazol-2-ylidene and 2,2,6,6-tetramethylpiperidin-1-ol) [51], and XORMUP (the dimer between 1,3-bis[2,6-bis(propan-2-yl)phenyl]-imidazol-2-ylidene and (2,4,6,2 ,4 ,6 -hexamethyl[1,1 :3 ,1 -terphenyl]-2 -yl)boronic acid) [124]. Interestingly, in the last system, there is the -B(OH) 2 group, of which, however, only one -O-H bond is involved in the formation of the C⋯H-O bridge with the carbene molecule.
The examples presented show clearly that the carbene carbon can actually engage in hydrogen bonding quite easily by acting as a proton acceptor. It can also be seen that these carbenes have large substituents (most often in the form of mesityl groups) in the 1 and 3 positions of the imidazol-2-ylidene parent molecule. Therefore, the study of these types of systems using theoretical methods is also very important [80]. However, the presence of bulky substituents significantly increases the size of the systems and thus significantly increases the computational cost of theoretical research. Moreover, their presence increases the probability of the occurrence of interactions other than C⋯H-D, which additionally complicates the research. Probably for these reasons the amount of theoretical studies on this type of systems is small. As stated in the Introduction, the purpose of this article is to fill this gap to some extent. The investigated systems will be described in the next subsection.

3.2. Investigated Systems

The subject of the research in this article is the hydrogen bonds between imidazol-2-ylidene (I) and its seven derivatives having a pair of R substituents (methyl = Me, iso-propyl = i Pr, tert-butyl = t Bu, phenyl = Ph, mesityl = Mes, 2,6-diisopropylphenyl = Dipp, 1-adamantyl = Ad) attached to both nitrogens and five popular and important H-D proton donors (HF, HCN, H 2 O, MeOH, NH 3 ). Dimers (complexes) having this bond will be denoted as IR 2 HD for simplicity. Apart from the description of C⋯H-D hydrogen bonds in these systems, an important topic will be the possible presence of other accompanying bonds as they can sometimes significantly influence the structure of a complex [80]. Values of some selected parameters characterizing the C⋯H-D hydrogen bond or the proton-donor molecule in the IR 2 HD dimers are listed in Table 1. Because imidazol-2-ylidene is somewhat different in terms of its binding capacity, the imidazol-2-ylidene complexes will be discussed first.

3.3. Imidazol-2-ylidene Complexes

Compared to all other molecules, imidazol-2-ylidene is distinguished by the presence of highly polar N-H bonds (Figure 3), which are good proton donors and therefore can eagerly form hydrogen bonds as long as there is an atom with good proton-acceptor properties, such as O, N or a halogen atom, in their presence.
The stable dimers between imidazol-2-ylidene and HF, HCN, H 2 O, MeOH or NH 3 are shown in Figure 4.
In the case of HF and HCN, the obtained hydrogen bonds are linear, the former of which is very strong, which is expressed by a large value of the determined dissociation (16.1 kcal/mol) and binding ( 13.7 kcal/mol) energies, a short distance C⋯H (1.676 Å) and a fairly high value of the electron density determined at the critical point ( ρ bcp ) of this bond (0.065 a.u.). In addition, the formation of the C⋯H-F hydrogen bond leads to a significant elongation of the H-F bond (+0.061 Å) and a very large red-shift of its stretching vibration frequency (−1317 cm 1 ). On the other hand, the interaction in the I⋯HCN dimer is almost half as weak (8.6 kcal/mol), which is also related to the much longer C⋯H distance (2.138 Å) and the much lower ρ BCP value (0.025 a.u.). The elongation of the C-H bond in HCN and the red-shift of its stretching vibration frequency are only +0.029 Å and −397 cm 1 , respectively.
As can be seen in Figure 4, unlike I⋯HF and I⋯HCN, in the remaining systems, i.e., I⋯NH 3 , I⋯H 2 O and I⋯MeOH, the C⋯H-D hydrogen bond is not linear due to the presence of another hydrogen bond. The presence of this accompanying hydrogen bond results not only from the presence of a strongly polar N-H bond in the imidazol-2-ylidene molecule, but above all from the presence of the D atom having an easily accessible lone electron pair. This possibility leads to the formation of the N-H⋯O hydrogen bond in the dimers I⋯H 2 O and I⋯MeOH and N-H⋯N in the I⋯NH 3 dimer. Importantly, this situation, i.e., the simultaneous presence of two strong hydrogen bonds, means that the determined values of the dissociation energy do not represent strengths of the C⋯HD hydrogen bonds only, because, as representing intermolecular interactions globally, they should rather be assigned to both the hydrogen bonds, i.e., C⋯H-D and the accompanying one. Moreover, it is not easy to extract energy values for individual bonds. Nevertheless, as already mentioned in the Theoretical Methods section, Emamian et al. [116] have recently shown that among many wave function-based HB descriptors (including those based on QTAIM), the electron density determined at the bond critical point of HB best correlates with the binding (in fact, interaction) energy of this HB, giving Equation (1). As a result, they proposed using this equation for a quick estimate of the energy of hydrogen bonds forming networks [116]. Using this proposal, in the case of the I⋯NH 3 dimer, the value is −4.6 kcal/mol for the N-H⋯N hydrogen bond and only −1.7 kcal/mol for C⋯H-N. Thus, the N-H⋯N bond is definitely stronger than N-H⋯C, which is also expressed in a definitely shorter distance H⋯N (2.066 Å) than H⋯C (2.521 Å). However, in the case of I⋯H 2 O and I⋯MeOH dimers, the O-H⋯C bond is definitely stronger (−5.3 and −5.5 kcal/mol, respectively) than the accompanying N-H⋯O bond (−1.9 and −2.4 kcal/mol, respectively). This is also reflected in shorter C⋯H distance (2.035 Å) than H⋯O (2.319 and 2.313 Å, respectively).
The dimers I⋯NH 3 , I⋯H 2 O and I⋯MeOH are good simple examples clearly showing that imidazol-2-ylidene willingly forms accompanying hydrogen bonds by the strongly polar N-H bond. The formation of such bonds is particularly easy when the D atom has lone electron pairs readily available, as is the case, for example, of N and O atoms. The formation of these bonds is prevented after replacing the H atoms in both N-H bonds with non-polar R substituents. However, in addition to the frequently mentioned steric effects, this type of modification of the imidazol-2-ylidene molecule may lead to other accompanying secondary interactions, which is discussed in the other subsections.

3.4. The IR 2 HD Dimers

3.4.1. The IR 2 HF and IR 2 HCN Dimers

Due to the fairly large similarity in the values of the C-H-D angle, i.e., the local geometry of the C⋯H-D hydrogen bond, in the HF and HCN dimers, it is convenient to consider these dimers together. The fully optimized structures of the dimers IR 2 HF and IR 2 HCN are shown in Figure 5 along with the obtained values of dissociation energy, distance C⋯H, angle C-H-D and value of the electron density at the bond critical point of the C⋯H hydrogen bond.
First, let us note that, as with the I⋯HF and I⋯HCN dimers (Figure 4), the HF dimers are much stronger than their HCN counterparts. The dissociation energies are in the 17.5–19.9 kcal/mol range while the HCN dimer values are only 9.4–11.4 kcal/mol. In the former case, the strongest complex is IDipp 2 HF, and the weakest is IPh 2 HF, while in the latter, the strongest complex is also the one involving IDipp 2 , but the weakest dimer is formed by IMe 2 instead. Much greater strength of hydrogen bonds C⋯H in IR 2 HF dimers than IR 2 HCN is also visible in much smaller distances H⋯C (1.599–1.632 Å vs 2.055–2.105 Å), much higher values of ρ bcp (0.073–0.079 a.u. vs 0.027–0.030 a.u.), much greater extensions of the H-D bond (0.075–0.088 Å vs 0.031–0.040 Å) and much higher values of the ν HD stretching vibration frequency red-shift (from −1591 up to −1814 cm 1 vs. from −462 up to −548 cm 1 ). It is worth noting that although the dimers involving Ad (i.e., IAd 2 HF and IAd 2 HCN) are not the strongest (although the difference compared to their counterparts involving Dipp is small, especially in the case of the HCN dimer), the effect of extending the proton-donor bond and red-shift values are in these systems greatest. Therefore, it can be concluded that together with IDipp 2 , IAd 2 also forms strong C⋯H hydrogen bonds, which have the greatest impact on the characteristics of HF and HCN proton-donor molecules. It is also worth noting that for both HF and HCN, the dissociation energies obtained for the dimers shown in Figure 5 are clearly greater than the energies obtained for I⋯HF and I⋯HCN (16.1 and 8.6 kcal/mol, respectively, as shown in Table 1 and Figure 4), and thus the substitution of hydrogen atoms in both N-H bonds in imidazol-2-ylidene even by methyl groups increases the strength of the dimer. This can most likely be explained by the weak inductive effect of these groups, which increases the charge on lone electron pairs.
It may be asked why the systems with IDipp 2 and IAd 2 are characterized by the highest values of dissociation energies (of course also the shortest H⋯C distances and the highest ρ bcp , Δ d HD and red-shift values). Very often, the increase in dimer strength is due to the presence of additional competing interactions, the presence of which is often suggested by the significant non-linearity of the dominant hydrogen bond, as was the case with the dimers I⋯H 2 O, I⋯MeOH and I⋯NH 3 shown in Figure 4 and discussed earlier. However, the reference to the C-H-D angle value can be very deceptive because these values for both IDipp 2 and IAd 2 are exactly 180 , so the C⋯H-D hydrogen bonds in the dimers of these carbenes with HF and HCN are linear. What is more, interestingly, carbenes IPh 2 , and especially IMes 2 , are characterized by a distinct nonlinearity of C⋯H-D hydrogen bonds. In these cases, the non-linearity may indeed result from the presence of other competing interactions which affect the alignment of the proton-donor molecule and thus also the geometry of the C⋯H-D hydrogen bond.
In order to search for such competitive interactions, the determination of a molecular graph defined by QTAIM [110,111,112] may be a particularly helpful tool. Figure 6 presents molecular graphs of several selected IR 2 HF and IR 2 HCN dimers.
As can be clearly seen, in each of the examples, the molecular graph shows the presence of accompanying interatomic interactions, which are indicated by color-coded arrows. IPh 2 HF is a simple example in which the molecular graph suggests the presence of two accompanying C-H⋯F type hydrogen bonds. A similar situation occurs in the slightly larger IMes 2 HF, where the proton-donating C-H bond is derived from one of the methyl groups in the mesityl group. Both these bonds should be very weak, which is suggested by low electron density values (ca. 0.008 a.u.). Using Equation (1) gives the value of −1.0 kcal/mol. In contrast, in the IMes 2 HCN dimer, a similar pair of bond paths indicates the presence of CH⋯C-type hydrogen bonds, which should be even weaker ( ρ bcp amounts to ca. 0.0034 a.u. only, and therefore, the bonding effect is negligible) than C-H⋯F. The IAd 2 HF dimer is a very simple example of a system featuring a C⋯F tetrel bonding pair. This time the binding energy value is significant (ca. −1.4 kcal/mol), suggesting that the pair of these interactions may contribute significantly to the overall binding effect of the dimer.
The presence of 2,6-diisopropylphenyl (Dipp) substituent allows a particularly large number of accompanying interactions [80]. Apart from pairs of hydrogen bonds of the type C⋯H-F (or C⋯H-C in IDipp 2 HCN), both these dimers, i.e., IDipp 2 HF and IDipp 2 HCN, experience the presence of bond paths between the hydrogen atom of a C-H bond and the carbene carbon atom. Therefore, this result suggests that apart from the dominant C⋯H-D hydrogen bond (D = F or C), the carbene carbon atom engages in two additional C⋯H-C hydrogen bonds, which, however, should be very weak (ca. −0.9 and −1.1 kcal/mol in IDipp 2 HF and IDipp 2 HCN, respectively). In the case of IDipp 2 HCN, the molecular graph also suggests the presence of a pair of intramolecular hydrogen bonds of the C-H⋯N type, as nitrogen atoms belong to the imidazol-2-ylidene ring. Since the values of ρ bcp are significant (ca. 0.016 a.u.), the use of Equation (1) yields the value of ca. −2.8 kcal/mol. Apart from the aforementioned accompanying hydrogen bonds, the molecular graphs of dimers IDipp 2 HF and IDipp 2 HCN also show the presence of numerous C-H⋯H-C interactions, which, however, seems to be a fairly common feature in the case of crowded systems with many C-H bonds [80,125,126,127].
As shown (Figure 6), molecular graphs suggest the presence of many different interactions, including intermolecular ones. Their role in the overall binding of the system and their influence on the overall structure is not easy to quantify, although there are many descriptors for individual, i.e., local, interactions [116]. One such descriptor is the electron density value determined at the bond critical point of a given interaction, which, as already mentioned, best correlates with the binding energy [116]. In the case of IR 2 HF and IR 2 HCN dimers, the estimates based on Equation (1) show that these interactions should be much weaker than the dominant C⋯H-D hydrogen bond. Nevertheless, also taking into account their large number, it is highly likely that their presence can influence the overall structure of the dimer. In particular, their presence can significantly influence the angle values as they are generally associated with small force constants.

3.4.2. The IR 2 H 2 O and IR 2 MeOH Dimers

The significant non-linearity of the hydrogen bonds in the dimers I⋯H 2 O and I⋯MeOH already indicates the presence of additional intermolecular interactions, and indeed, as shown in Figure 4 and discussed previously, there is an additional hydrogen bond of the N-H⋯O type in both of these dimers. The structures of these dimers also prove that lone electron pairs on oxygen are readily available, and therefore, both H 2 O and MeOH will be willing to engage in the formation of accompanying hydrogen bonds. The resulting IR 2 H 2 O and IR 2 MeOH dimer structures are shown in Figure 7.
Among the dimers with water, IDipp 2 H 2 O and IMes 2 H 2 O are characterized by the highest dissociation energy value (11.2 and 10.7 kcal/mol, respectively), and the lowest is for IMe 2 H 2 O (9.0 kcal/mol). It is similar in the case of dimers with MeOH: IDipp 2 MeOH (14.4 kcal/mol), IMes 2 MeOH (13.1 kcal/mol), IMe 2 MeOH (10.4 kcal/mol). As was the case with dimers involving HF or HCN, the presence of any R substituent clearly increases the value of the dissociation energy. Clearly (Table 1, Figure 7), for a given carbene IR 2 , MeOH forms a stronger complex than H 2 O. This is also reflected in the shorter distances C⋯H. For example, for IPh 2 the distances are 1.927 and 1.954 Å, respectively, and for IAd 2 , they are 1.931 and 1.991 Å, respectively. The dissociation energies for H 2 O are rather similar and generally slightly lower than those for HCN. It is worth noting that in the case of dimers with H 2 O and MeOH, the highest (i.e., most negative) values of E b , i.e., the binding energy obtained by Equation (1), have been obtained for the I t Bu 2 carbene, characterized by the formation of an almost linear (ca. 174 ) C⋯H-O bond. On the contrary, the smallest values of E b have been obtained for IMe 2 (−6.8 and −7.2 kcal/mol for H 2 O and MeOH, respectively). However, they are much larger than in the unsubstituted imidazol-2-ylidene (−5.3 and −5.5 kcal/mol, respectively).
Taking into account that the molecular graphs obtained for the IR 2 HF and IR 2 HCN dimers (Figure 6) in many cases indicated the presence of two, four or even more additional coexisting interactions, it is now worth analyzing the molecular graphs for IR 2 H 2 O and IR 2 MeOH dimers. The most interesting examples are shown in Figure 8.
In the dimers shown in Figure 8, generally one (e.g., I i Pr 2 H 2 O) or two (e.g., I t Bu 2 H 2 O) bond paths are present for the accompanying C-H⋯O hydrogen bonds (indicated by red arrows). It is clear, however, that for some carbenes, the MeOH dimers contain additional bond paths for C-H⋯H-C contacts (indicated by yellow arrows). Moreover, the presence of a methyl group in MeOH allows in some cases (see the dimers with IMes 2 and IDipp 2 ) an additional (i.e., compared with H 2 O) C-H π -type hydrogen bond (green arrow). As was the case with the HF and HCN dimers, IDipp 2 produces the most intricate molecular graph, suggesting the existence of many interactions accompanying the main C⋯H-O hydrogen bond involving the carbene carbon atom. In addition to three C-H⋯H-C interactions, numerous accompanying hydrogen bonds are present. Taking this into account, this fact may explain the exceptionally high value of dissociation energy in dimers with IDipp 2 , and especially in IDipp 2 MeOH (Table 1 and Figure 7). Nevertheless, all these accompanying interactions should be much weaker than the leading C⋯H-O hydrogen bond involving the carbene carbon atom. For example, C-H⋯O hydrogen bonds are generally in the range of −0.6 to −1.5 kcal/mol (e.g., −0.8 kcal/mol in I t Bu 2 H 2 O, −1.0 kcal/mol in I t Bu 2 MeOH, −1.1 kcal/mol in I i Pr 2 H 2 O and I i Pr 2 MeOH, and −1.4 kcal/mol in IMes 2 H 2 O), however they are slightly stronger in IAd 2 H 2 O (−1.6 kcal/mol) and especially in IAd 2 MeOH (−1.9 kcal/mol). The C-H π bonds in IMes 2 MeOH and IDipp 2 MeOH are much weaker (−0.6 kcal/mol). On the other hand, the energy of C-H⋯H-C interactions is negligible, below −0.5 kcal/mol (e.g., −0.4 kcal/mol in IAd 2 MeOH). On the other hand, the intramolecular C-H⋯N hydrogen bonds in IDipp 2 H 2 O and IDipp 2 MeOH are significantly strong (−2.8 kcal/mol). All these energy values should be compared with the energies of the dominant C⋯H-O hydrogen bond, which in the case of IR 2 H 2 O and IR 2 MeOH dimers range from −6.6 to −8.4 kcal/mol (Table 1).

3.4.3. The IR 2 NH 3 Dimers

The values of dissociation ( D 0 ) and binding ( E b ) energies listed in Table 1 show that the dimers with NH 3 as well as the C⋯H-N hydrogen bonds in them should be the weakest. Namely, the D 0 values are from 4.5 kcal/mol in I t Bu 2 NH 3 to 7.3 kcal/mol in IDipp 2 NH 3 . The E b values form a fairly narrow range from −3.2 kcal/mol (not including I⋯NH 3 with a much lower value of −1.6 kcal/mol) to −3.7 for IDipp 2 NH 3 . Such small values of D 0 and E b may result from the fact that the N-H bond in ammonia is a worse proton donor than the O-H bond in water or methanol. The C-H-D angle ( α CHD ) values show that in none of the dimers, the C⋯H-N hydrogen bond is linear (although in I t Bu 2 NH 3 and IAd 2 NH 3 this angle is ca. 174 ); quite the opposite, the deviation from linearity is significant (157 –166 ), although much smaller than in I⋯NH 3 (123 ). Such large deviations from linearity may indicate the presence of additional coexisting interactions. The structures of the IR 2 NH 3 dimers are shown in Figure 9, while their molecular graphs can be found in Figure 10.
A careful comparison of the molecular graphs obtained for the IR 2 H 2 O and IR 2 MeOH dimers (Figure 8) with the molecular graphs obtained for the IR 2 NH 3 dimers (Figure 10) shows some impoverishment in the second case. In addition, there is a visible change in the nature of some accompanying interactions in some systems with the same IR 2 carbene. For example, the molecular graph of the I i Pr 2 NH 3 dimer contains only one bond path for the interaction (C-H⋯N) accompanying the dominant C⋯H-N hydrogen bond, while I i Pr 2 MeOH in addition to the additional C-H⋯O also contains a bond path for C-H⋯H-C. The I t Bu 2 H 2 O dimer contains two bond paths for C-H⋯O hydrogen bonds, whereas I t Bu 2 NH 3 also has two additional bond paths, but for C-H⋯H-C interactions. The situation is similar for the molecular graph involving Ad; in the case of dimer with H 2 O, there are two bond paths for C-H⋯O hydrogen bonds, whereas in the case of IAd 2 NH 3 , there are two bond paths for C-H⋯H-N. Keeping in mind that C-H⋯H-C interactions are weaker than C-H⋯O/N hydrogen bonds, together with a smaller number of accompanying interactions, this finding may explain the weaker strength of complexes with NH 3 . Similarly, it can be seen that IDipp 2 creates many bond paths, of which as many as three (except three for N-H bonds, of course) lead to the nitrogen atom of the ammonia molecule. Two of them represent very weak C-H⋯N hydrogen bonds (−0.1 and −1.3 kcal/mol), while the other, interestingly, determines the N π contact (−0.5 kcal/mol). Of course, all these interactions are much weaker than the dominant C⋯H-N hydrogen bond (−3.2 kcal/mol). Nevertheless, their multitude makes the dissociation energy for IDipp 2 NH 3 clearly the highest (7.3 kcal/mol) among all the IR 2 NH 3 dimers (see Table 1 or Figure 9). Of course, it should be noted that a similar dissociation energy value was obtained for a simple I⋯NH 3 dimer (7.2 kcal/mol); however, as already shown (Figure 4), such a large value in this case results from the presence of a strong N-H⋯N hydrogen bond (−4.6 kcal/mol), which overwhelms the weaker C⋯H-N bond (−1.6 kcal/mol).

3.5. Relationship between Dimer Strength and the C⋯H Distance

As has been shown in the previous subsections, many of the dimers considered here have some secondary interactions accompanying the leading C⋯H-D hydrogen bond. The strength of at least some of them may be considerable. Therefore, one should not expect a good overall (i.e., for all the IR 2 HD dimers) correlation between the dissociation energy, D 0 , and the distance C⋯H. Moreover, good correlations are also not to be expected even in cases where individual proton-donor molecules are considered. Of course, the exception here is HF, and perhaps HCN, which more often forms linear or nearly linear C⋯H hydrogen bonds (Table 1). This is confirmed in Figure 11 (left), where it is shown that the R 2 (coefficient of determination) values for the linear correlation between D 0 and d C H are very poor for NH 3 (even with the rejection of the clearly outlier point for I⋯NH 3 ), H 2 O and MeOH (0.217, 0.319, 0.550, respectively) and much better for HCN (0.886) and especially for HF (0.930).
This result clearly shows that in the case of more complex dimers, i.e., having more significant intermolecular contacts, one should not expect a good relationship between the dimer dissociation energy and the length of any single intermolecular contact. Rather, a good correlation can be expected in the case of a clearly local parameter, which is the binding (interaction) energy determined from Equation (1), and therefore based on the electron density value at the critical point of the C⋯H-D hydrogen bond. Indeed, as shown in Figure 11 (right), for HF, MeOH and NH 3 , the linear correlations between E b and d C H are very good ( R 2 > 0.98) and only slightly worse for HCN and H 2 O (0.968 and 0.962, respectively). Importantly, for all the IR 2 HD dimers, an excellent exponential relationship, E b = 755.8 e 2.38 d C H , has been found (black curve in Figure 11). This result is not unexpected as a good common exponential relationship between ρ bcp and the hydrogen bond length has been found for both different X and Y atoms in the X-H⋯Y hydrogen bond [128,129,130,131].

4. Conclusions

In this article, forty dimers formed by imidazol-2-ylidene or its derivatives, in which the hydrogen atoms of both N-H bonds were replaced by important and popular substituents of increasing complexity (methyl = Me, iso-propyl = i Pr, tert-butyl = t Bu, phenyl = Ph, mesityl = Mes, 2,6-diisopropylphenyl = Dipp, 1-adamantyl = Ad), and five fundamental proton donor (HD) molecules (HF, HCN, H 2 O, MeOH, NH 3 ), have been studied. Although the most important goal was to describe the hydrogen bond formed by the carbene carbon atom, C⋯H-D, an equally important issue, which is, in the author’s opinion, insufficiently described in the literature, was the possibility of the formation of various accompanying interactions and their influence on the structure of the considered dimers. Despite different interaction abilities represented by carbenes and especially HD molecules, the following general trends have been found, in addition to many specific results:
  • For a given carbene, dissociation energies of the IR 2 HD dimers increase in the following order: NH 3 < H 2 O < HCN ≤ MeOH ≪ HF.
  • For a given HD molecule (HF, HCN, H 2 O, MeOH, or NH 3 ), IDipp 2 , i.e., 1,3-bis[2,6-diisopropylphenyl]imidazol-2-ylidene, has been found to form the strongest dimers. This has been attributed to the multiplicity of various interactions accompanying the dominant C⋯H-D hydrogen bond.
  • The substitution of hydrogen atoms in both N-H bonds of the imidazol-2-ylidene molecule by Me, i Pr, t Bu, Ph, Mes, Dipp or Ad groups leads to stronger dimers with HF, HCN, H 2 O or MeOH.
The article clearly shows that various types of secondary interactions, thus often omitted in the analysis, can have a significant impact on the structure of a molecular system and its strength. The results should improve our understanding of carbene chemistry and the role of intermolecular interactions.

Funding

This research received no external funding.

Data Availability Statement

Data available from the author on reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this article:
NHCN-heterocyclic carbene
Iimidazol-2-ylidene
IR 2 R-substituted imidazol-2-ylidene
Memethyl group
i Prisopropyl group
t Butert-butyl group
Phphenyl group
Mesmesityl group
Dipp2,6-diisopropylphenyl group
Adadamantyl group
QTAIMquantum theory of atoms in molecules

References

  1. Wanzlick, H.W. Aspects of Nucleophilic Carbene Chemistry. Angew. Chem. Int. Ed. Engl. 1962, 1, 75–80. [Google Scholar] [CrossRef]
  2. Kirmse, W. Carbene Chemistry; Academic Press: Cambridge, MA, USA, 1964. [Google Scholar]
  3. Hoffmann, R.; Zeiss, G.D.; Van Dine, G.W. The Electronic Structure of Methylenes. J. Am. Chem. Soc. 1968, 90, 1485–1499. [Google Scholar] [CrossRef]
  4. Gleiter, R.; Hoffmann, R. On Stabilizing a Singlet Methylene. J. Am. Chem. Soc. 1968, 90, 5457–5460. [Google Scholar] [CrossRef]
  5. Baird, N.C.; Taylor, K.F. Multiplicity of the Ground State and Magnitude of the T1–S0 gap in Substituted Carbenes. J. Am. Chem. Soc. 1978, 100, 1333–1338. [Google Scholar] [CrossRef]
  6. Harrison, J.F.; Liedtke, R.C.; Liebman, J.F. The Multiplicity of Substituted Acyclic Carbenes and Related Molecules. J. Am. Chem. Soc. 1979, 101, 7162–7168. [Google Scholar] [CrossRef]
  7. Schoeller, W.W. Electrophilicity and nucleophilicity in singlet carbenes. II. Electrophilic selectivity. Tetrahedron Lett. 1980, 21, 1509–1510. [Google Scholar] [CrossRef]
  8. Mueller, P.H.; Rondan, N.G.; Houk, K.N.; Harrison, J.F.; Hooper, D.; Willen, B.H.; Liebman, J.F. Carbene Singlet–Triplet Gaps. Linear Correlations with Substituent π Donation. J. Am. Chem. Soc. 1981, 103, 5049–5052. [Google Scholar] [CrossRef]
  9. Pople, J.A.; Raghavachari, K.; Frisch, M.J.; Binkley, J.S.; Schleyer, P.v.R. Comprehensive Theoretical Study of Isomers and Rearrangement Barriers of Even-Electron Polyatomic Molecules HmABHn (A, B = C, N, O, and F). J. Am. Chem. Soc. 1983, 105, 6389–6398. [Google Scholar] [CrossRef]
  10. Pople, J.A. A theoretical search for the methylenefluoronium ylide. Chem. Phys. Lett. 1986, 132, 144–146. [Google Scholar] [CrossRef]
  11. Schubert, U. Advances in Metal Carbene Chemistry; Springer: Dordrecht, The Netherlands, 1989. [Google Scholar]
  12. Arduengo, A.J., III; Kline, M.; Calabrese, J.C.; Davidson, F. Synthesis of a Reverse Ylide from a Nucleophilic Carbene. J. Am. Chem. Soc. 1991, 113, 9704–9705. [Google Scholar] [CrossRef]
  13. Arduengo, A.J., III; Dias, H.V.R.; Calabrese, J.C.; Davidson, F. A Stable Carbene-Alane Adduct. J. Am. Chem. Soc. 1992, 114, 9724–9725. [Google Scholar] [CrossRef]
  14. Arduengo, A.J., III; Rasika Dias, H.V.; Davidson, F.; Harlow, R.L. Carbene adducts of magnesium and zinc. J. Organomet. Chem. 1993, 462, 13–18. [Google Scholar] [CrossRef]
  15. Kuhn, N.; Kratz, T.; Henkel, G. A Stable Carbene Iodine Adduct: Secondary Bonding in 1,3-Diethyl-2-iodo-4,5-dimethylimidazolium Iodide. J. Chem. Soc. Chem. Commun. 1993, 1778–1779. [Google Scholar] [CrossRef]
  16. Arduengo, A.J., III; Gamper, S.F.; Tamm, M.; Calabrese, J.C.; Davidson, F.; Craig, H.A. A Bis(carbene)–Proton Complex: Structure of a C–H–C Hydrogen Bond. J. Am. Chem. Soc. 1995, 117, 572–573. [Google Scholar] [CrossRef]
  17. Herrmann, W.A.; Runte, O.; Artus, G. Synthesis and structure of an ionic beryllium–“carbene” complex. J. Organomet. Chem. 1995, 501, Cl–C4. [Google Scholar] [CrossRef]
  18. Boehme, C.; Frenking, G. Electronic Structure of Stable Carbenes, Silylenes, and Germylenes. J. Am. Chem. Soc. 1996, 118, 2039–2046. [Google Scholar] [CrossRef]
  19. Alkorta, I.; Elguero, J. Carbenes and Silylenes as Hydrogen Bond Acceptors. J. Phys. Chem. 1996, 100, 19367–19370. [Google Scholar] [CrossRef]
  20. Li, X.-W.; Su, J.; Robinson, G.H. Syntheses and molecular structure of organo-group 13 metal carbene complexes. Chem. Commun. 1996, 2683–2684. [Google Scholar] [CrossRef]
  21. Herrmann, W.A.; Köcher, C. N-Heterocyclic Carbenes. Angew. Chem. Int. Ed. Engl. 1997, 36, 2162–2187. [Google Scholar] [CrossRef]
  22. Arduengo, A.J., III; Davidson, F.; Krafczyk, R.; Marshall, W.J.; Tamm, M. Adducts of Carbenes with Group II and XII Metallocenes. Organometallics 1998, 17, 3375–3382. [Google Scholar] [CrossRef]
  23. Hibbs, D.E.; Hursthouse, M.B.; Jones, C.; Smithies, N.A. Synthesis, crystal and molecular structure of the first indium trihydride complex, [InH3CN(Pri)C2Me2N(Pri)]. Chem. Commun. 1998, 869–870. [Google Scholar] [CrossRef]
  24. Goumri-Magnet, S.; Polishchuck, O.; Gornitzka, H.; Marsden, C.J.; Baceiredo, A.; Bertrand, G. The Electrophilic Behavior of Stable Phosphanylcarbenes Towards Phosphorus Lone Pairs. Angew. Chem. Int. Ed. 1999, 38, 3727–3729. [Google Scholar] [CrossRef]
  25. Bourissou, D.; Guerret, O.; Gabbaï, F.P.; Bertrand, G. Stable Carbenes. Chem. Rev. 2000, 100, 39–91. [Google Scholar] [CrossRef] [PubMed]
  26. Bertrande, G. Carbene Chemistry: From Fleeting Intermediates to Powerful Reagents; FontisMedia S.A.: Lausanne, Switzerland; Marcel Dekker, Inc.: New York, NY, USA, 2002. [Google Scholar]
  27. Merceron, N.; Miqueu, K.; Baceiredo, A.; Bertrand, G. Stable (Amino)(phosphino)carbenes: Difunctional Molecules. J. Am. Chem. Soc. 2002, 124, 6806–6807. [Google Scholar] [CrossRef] [PubMed]
  28. Wang, D.; Wurst, K.; Buchmeiser, M.R. N-heterocyclic carbene complexes of Zn(II): Synthesis, X-ray structures and reactivity. J. Organomet. Chem. 2004, 689, 2123–2130. [Google Scholar] [CrossRef]
  29. Moss, R.A.; Platz, M.S.; Jones, M., Jr. (Eds.) Reactive Intermediate Chemistry; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2005. [Google Scholar]
  30. Scott, N.M.; Nolan, S.P. Stabilization of Organometallic Species Achieved by the Use of N-Heterocyclic Carbene (NHC) Ligands. Eur. J. Inorg. Chem. 2005, 1815–1828. [Google Scholar] [CrossRef]
  31. Frenking, G.; Solà, M.; Vyboishchikov, S.F. Chemical bonding in transition metal carbene complexes. J. Organomet. Chem. 2005, 690, 6178–6204. [Google Scholar] [CrossRef]
  32. Garrison, J.C.; Youngs, W.J. Ag(I) N-Heterocyclic Carbene Complexes: Synthesis, Structure, and Application. Chem. Rev. 2005, 105, 3978–4008. [Google Scholar] [CrossRef]
  33. Nolan, S.P. N-Heterocyclic Carbenes in Synthesis; Wiley-VCH: Weinheim, Germany, 2006. [Google Scholar]
  34. Carey, F.A.; Sundberg, R.J. Carbenes, Part B: Reactions and Synthesis. Advanced Organic Chemistry; Springer: New York, NY, USA, 2007. [Google Scholar]
  35. Díez-González, S.; Nolan, S.P. Stereoelectronic parameters associated with N-heterocyclic carbene (NHC) ligands: A quest for understanding. Coord. Chem. Rev. 2007, 251, 874–883. [Google Scholar] [CrossRef]
  36. Jacobsen, H.; Correa, A.; Poater, A.; Costabile, C.; Cavallo, L. Understanding the M–(NHC) (NHC = N-heterocyclic carbene) bond. Coord. Chem. Rev. 2009, 253, 687–703. [Google Scholar] [CrossRef]
  37. de Frémont, P.; Marion, N.; Nolan, S.P. Carbenes: Synthesis, properties, and organometallic chemistry. Coord. Chem. Rev. 2009, 253, 862–892. [Google Scholar] [CrossRef]
  38. Jabłoński, M.; Palusiak, M. Divalent carbon atom as the proton acceptor in hydrogen bonding. Phys. Chem. Chem. Phys. 2009, 11, 5711–5719. [Google Scholar] [CrossRef] [PubMed]
  39. Li, Q.; Wang, H.; Liu, Z.; Li, W.; Cheng, J.; Gong, B.; Sun, J. Ab Initio Study of Lithium-Bonded Complexes with Carbene as an Electron Donor. J. Phys. Chem. A 2009, 113, 14156–14160. [Google Scholar] [CrossRef] [PubMed]
  40. Arrowsmith, M.; Hill, M.S.; MacDougall, D.J.; Mahon, M.F. A Hydride-Rich Magnesium Cluster. Angew. Chem. Int. Ed. 2009, 48, 4013–4016. [Google Scholar] [CrossRef] [PubMed]
  41. Arnold, P.L.; Casely, I.J.; Turner, Z.R.; Bellabarba, R.; Tooze, R.B. Magnesium and zinc complexes of functionalised, saturated N-heterocyclic carbene ligands: Carbene lability and functionalisation, and lactide polymerisation catalysis. Dalton Trans. 2009, 35, 7236–7247. [Google Scholar] [CrossRef] [PubMed]
  42. Wang, Y.; Robinson, G.H. Unique homonuclear multiple bonding in main group compounds. Chem. Commun. 2009, 5201–5213. [Google Scholar] [CrossRef] [PubMed]
  43. Li, Q.; Wang, Y.; Liu, Z.; Li, W.; Cheng, J.; Gong, B.; Sun, J. An unconventional halogen bond with carbene as an electron donor: An ab initio study. Chem. Phys. Lett. 2009, 469, 48–51. [Google Scholar] [CrossRef]
  44. Hindi, K.M.; Panzner, M.J.; Tessier, C.A.; Cannon, C.L.; Youngs, W.J. The Medicinal Applications of Imidazolium Carbene–Metal Complexes. Chem. Rev. 2009, 109, 3859–3884. [Google Scholar] [CrossRef]
  45. Wang, Y.; Xie, Y.; Abraham, M.Y.; Wei, P.; Schaefer, H.F., III; Schleyer, P.v.R.; Robinson, G.H. A Viable Anionic N-Heterocyclic Dicarbene. J. Am. Chem. Soc. 2010, 132, 14370–14372. [Google Scholar] [CrossRef]
  46. Li, Z.-F.; Yang, S.; Li, H.-X. Theoretical prediction characters of unconventional weak bond with carbene as electron donors and Li-Y (Y = OH, H, F, NC and CN) as electron acceptors. J. Mol. Struct. THEOCHEM 2010, 952, 56–60. [Google Scholar]
  47. Wang, Y.; Xie, Y.; Abraham, M.Y.; Gilliard, R.J., Jr.; Wei, P.; Schaefer, H.F., III; Schleyer, P.v.R.; Robinson, G.H. Carbene-Stabilized Parent Phosphinidene. Organometallics 2010, 29, 4778–4780. [Google Scholar] [CrossRef]
  48. Abraham, M.Y.; Wang, Y.; Xie, Y.; Wei, P.; Schaefer, H.F., III; Schleyer, P.v.R.; Robinson, G.H. Carbene Stabilization of Diarsenic: From Hypervalency to Allotropy. Chem. Eur. J. 2010, 16, 432–435. [Google Scholar] [CrossRef]
  49. Mercs, L.; Albrecht, M. Beyond catalysis: N-Heterocyclic Carbene Complexes as Components for Medicinal, Luminescent, and Functional Materials Applications. Chem. Soc. Rev. 2010, 39, 1903–1912. [Google Scholar] [CrossRef] [PubMed]
  50. Patel, D.S.; Bharatam, P.V. Divalent N(I) Compounds with Two Lone Pairs on Nitrogen. J. Phys. Chem. A 2011, 115, 7645–7655. [Google Scholar] [CrossRef] [PubMed]
  51. Giffin, N.A.; Makramalla, M.; Hendsbee, A.D.; Robertson, K.N.; Sherren, C.; Pye, C.C.; Masuda, J.D.; Clyburne, J.A.C. Anhydrous TEMPO-H: Reactions of a good hydrogen atom donor with low-valent carbon centres. Org. Biomol. Chem. 2011, 9, 3672–3680. [Google Scholar] [CrossRef] [PubMed]
  52. Zhao, Q.; Feng, D.; Sun, Y.; Hao, J.; Cai, Z. Theoretical Investigations on the Weak Nonbonded C=S⋯CH2 Interactions: Chalcogen-Bonded Complexes With Singlet Carbene as an Electron Donor. Int. J. Quant. Chem. 2011, 111, 3881–3887. [Google Scholar] [CrossRef]
  53. Gilliard, R.J., Jr.; Abraham, M.Y.; Wang, Y.; Wei, P.; Xie, Y.; Quillian, B.; Schaefer, H.F., III; Schleyer, P.v.R.; Robinson, G.H. Carbene-Stabilized Beryllium Borohydride. J. Am. Chem. Soc. 2012, 134, 9953–9955. [Google Scholar] [CrossRef]
  54. Arrowsmith, M.; Hill, M.S.; Kociok-Köhn, G.; MacDougall, D.J.; Mahon, M.F. Beryllium-Induced C–N Bond Activation and Ring Opening of an N-Heterocyclic Carbene. Angew. Chem. Int. Ed. 2012, 51, 2098–2100. [Google Scholar] [CrossRef]
  55. Moss, R.A.; Doyle, M.P. Contemporary Carbene Chemistry; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2013. [Google Scholar]
  56. Gerbig, D.; Ley, D. Computational methods for contemporary carbene chemistry. WIREs Comput. Mol. Sci. 2013, 3, 242–272. [Google Scholar] [CrossRef]
  57. Nelson, D.J.; Nolan, S.P. Quantifying and understanding the electronic properties of N-heterocyclic carbenes. Chem. Soc. Rev. 2013, 42, 6723–6753. [Google Scholar] [CrossRef]
  58. Samanta, R.C.; De Sarkar, S.; Fröhlich, R.; Grimme, S.; Studer, A. N-Heterocyclic carbene (NHC) catalyzed chemoselective acylation of alcohols in the presence of amines with various acylating reagents. Chem. Sci. 2013, 4, 2177–2184. [Google Scholar] [CrossRef]
  59. Esrafili, M.D.; Mohammadirad, N. Insights into the strength and nature of carbene⋯halogen bond interactions: A theoretical perspective. J. Mol. Model. 2013, 19, 2559–2566. [Google Scholar] [CrossRef] [PubMed]
  60. Budagumpi, S.; Endud, S. Group XII Metal–N-Heterocyclic Carbene Complexes: Synthesis, Structural Diversity, Intramolecular Interactions, and Applications. Organometallics 2013, 32, 1537–1562. [Google Scholar] [CrossRef]
  61. Hopkinson, M.N.; Richter, C.; Schedler, M.; Glorius, F. An overview of N-heterocyclic carbenes. Nature 2014, 510, 485–496. [Google Scholar] [CrossRef] [PubMed]
  62. Bellemin-Laponnaz, S.; Dagorne, S. Group 1 and 2 and Early Transition Metal Complexes Bearing N-Heterocyclic Carbene Ligands: Coordination Chemistry, Reactivity, and Applications. Chem. Rev. 2014, 114, 8747–8774. [Google Scholar] [CrossRef]
  63. Visbal, R.; Concepción Gimeno, M. N-Heterocyclic Carbene Metal Complexes: Photoluminescence and Applications. Chem. Soc. Rev. 2014, 43, 3551–3574. [Google Scholar] [CrossRef]
  64. Santoro, O.; Nahra, F.; Cordes, D.B.; Slawin, A.M.Z.; Nolan, S.P.; Cazin, C.S.J. Synthesis, characterization and catalytic activity of stable [(NHC)H][ZnXY2] (NHC = N-Heterocyclic carbene, X, Y = Cl, Br) species. J. Mol. Catal. 2016, 423, 85–91. [Google Scholar] [CrossRef]
  65. Del Bene, J.E.; Alkorta, I.; Elguero, J. Carbon–Carbon Bonding between Nitrogen Heterocyclic Carbenes and CO2. J. Phys. Chem. A 2017, 121, 8136–8146. [Google Scholar] [CrossRef]
  66. Liu, M.; Li, Q.; Li, W.; Cheng, J. Carbene tetrel-bonded complexes. Struct. Chem. 2017, 28, 823–831. [Google Scholar] [CrossRef]
  67. Esrafili, M.D.; Sabouri, A. Carbene–aerogen bonds: An ab initio study. Mol. Phys. 2017, 115, 971–980. [Google Scholar] [CrossRef]
  68. Moss, R.A.; Wang, L.; Cang, H.; Krogh-Jespersen, K. Extremely reactive carbenes: Electrophiles and nucleophiles. J. Phys. Org. Chem. 2017, 30, e3555. [Google Scholar] [CrossRef]
  69. Nesterov, V.; Reiter, D.; Bag, P.; Frisch, P.; Holzner, R.; Porzelt, A.; Inoue, S. NHCs in Main Group Chemistry. Chem. Rev. 2018, 118, 9678–9842. [Google Scholar] [CrossRef] [PubMed]
  70. Alkorta, I.; Elguero, J. A LFER analysis of the singlet-triplet gap in a series of sixty-six carbenes. Chem. Phys. Lett. 2018, 691, 33–36. [Google Scholar] [CrossRef]
  71. Jabłoński, M. The first theoretical proof of the existence of a hydride-carbene bond. Chem. Phys. Lett. 2018, 710, 78–83. [Google Scholar] [CrossRef]
  72. Dagorne, S. Recent Developments on N-Heterocyclic Carbene Supported Zinc Complexes: Synthesis and Use in Catalysis. Synthesis 2018, 50, 3662–3670. [Google Scholar] [CrossRef]
  73. Walley, J.E.; Wong, Y.-O.; Freeman, L.A.; Dickie, D.A.; Gilliard, R.J., Jr. N-Heterocyclic Carbene-Supported Aryl- and Alk- oxides of Beryllium and Magnesium. Catalysts 2019, 9, 934. [Google Scholar] [CrossRef]
  74. Jabłoński, M. In search for a hydride-carbene bond. J. Phys. Org. Chem. 2019, 32, e3949. [Google Scholar] [CrossRef]
  75. Yourdkhani, S.; Jabłoński, M. Physical nature of silane⋯carbene dimers revealed by state-of-the-art ab initio calculations. J. Comput. Chem. 2019, 40, 2643–2652. [Google Scholar] [CrossRef]
  76. Procter, R.J.; Uzelac, M.; Cid, J.; Rushworth, P.J.; Ingleson, M.J. Low-Coordinate NHC–Zinc Hydride Complexes Catalyze Alkyne C–H Borylation and Hydroboration Using Pinacolborane. ACS Catal. 2019, 9, 5760–5771. [Google Scholar] [CrossRef]
  77. Dzieszkowski, K.; Barańska, I.; Mroczyńska, K.; Słotwiński, M.; Rafiński, Z. Organocatalytic Name Reactions Enabled by NHCs. Materials 2020, 13, 3574. [Google Scholar] [CrossRef]
  78. Specklin, D.; Fliedel, C.; Dagorne, S. Recent Representative Advances on the Synthesis and Reactivity of N-Heterocyclic-Carbene-Supported Zinc Complexes. Chem. Rec. 2021, 21, 1130–1143. [Google Scholar] [CrossRef] [PubMed]
  79. Jabłoński, M. Study of Beryllium, Magnesium, and Spodium Bonds to Carbenes and Carbodiphosphoranes. Molecules 2021, 26, 2275. [Google Scholar] [CrossRef]
  80. Jabłoński, M. Theoretical Study of N-Heterocyclic-Carbene–ZnX2 (X = H, Me, Et) Complexes. Materials 2021, 14, 6147. [Google Scholar] [CrossRef]
  81. Pauling, L. The Nature of the Chemical Bond; Cornell University Press: New York, NY, USA, 1960. [Google Scholar]
  82. Pimentel, G.C.; McClellan, A.L. The Hydrogen Bond; W.H. Freeman & Co.: San Francisco, CA, USA, 1960. [Google Scholar]
  83. Hamilton, W.C.; Ibers, J.A. Hydrogen Bonding in Solids; W. A. Benjamin: New York, NY, USA, 1968. [Google Scholar]
  84. Vinogradov, S.N.; Linnell, R.H. Hydrogen Bonding; Van Nostrand-Reinhold: Princeton, NJ, USA, 1971. [Google Scholar]
  85. Schuster, P.; Zundel, G.; Sandorfy, C. (Eds.) The Hydrogen Bond. Recent Developments in Theory and Experiments; North Holland: Amsterdam, The Netherlands, 1976; Volumes I–III. [Google Scholar]
  86. Schuster, P. Intermolecular Interactions: From Diatomics to Biopolymers; Pullman, B., Ed.; John Wiley: New York, NY, USA, 1978. [Google Scholar]
  87. Hobza, P.; Zahradník, R. Weak Intermolecular Interactions in Chemistry and Biology; Academia: Prague, Czech Republic, 1980. [Google Scholar]
  88. Jeffrey, G.A.; Saenger, W. Hydrogen Bonding in Biological Structures; Springer: Berlin/Heidelberg, Germany, 1991. [Google Scholar]
  89. Hadži, D. (Ed.) Theoretical Treatments of Hydrogen Bonding; John Wiley: Chichester, UK, 1997. [Google Scholar]
  90. Jeffrey, G.A. An Introduction to Hydrogen Bonding; Oxford University Press: New York, NY, USA, 1997. [Google Scholar]
  91. Scheiner, S. (Ed.) Molecular Interactions. From van der Waals to Strongly Bound Complexes; Wiley: Chichester, UK, 1997. [Google Scholar]
  92. Scheiner, S. Hydrogen Bonding: A Theoretical Perspective; Oxford University Press: New York, NY, USA, 1997. [Google Scholar]
  93. Desiraju, G.R.; Steiner, T. The Weak Hydrogen Bond in Structural Chemistry and Biology; Oxford University Press: New York, NY, USA, 1999. [Google Scholar]
  94. Grabowski, S.J. (Ed.) Hydrogen Bonding—New Insights; Springer: Dordrecht, The Netherlands, 2006. [Google Scholar]
  95. Maréchal, Y. The Hydrogen Bond and the Water Molecule; Elsevier: Amsterdam, The Netherlands, 2007. [Google Scholar]
  96. Gilli, G.; Gilli, P. The Nature of the Hydrogen Bond. Outline of a Comprehensive Hydrogen Bond Theory; Oxford University Press: Oxford, UK, 2009. [Google Scholar]
  97. Scheiner, S. Evaluation of DFT Methods to Study Reactions of Benzene With OH Radical. Int. J. Quant. Chem. 2012, 112, 1879–1886. [Google Scholar] [CrossRef]
  98. Li, Q.-Z.; Li, H.-B. Hydrogen Bonds Involving Radical Species. In Noncovalent Forces; Scheiner, S., Ed.; Springer: Cham, Switzerland, 2015. [Google Scholar]
  99. Chai, J.-D.; Head-Gordon, M. Long-range corrected hybrid density functionals with damped atom–atom dispersion correction. Phys. Chem. Chem. Phys. 2008, 10, 6615–6620. [Google Scholar] [CrossRef]
  100. Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. Phys. Rev. 1964, 136, B864–B871. [Google Scholar] [CrossRef]
  101. Kohn, W.; Sham, L.J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133–A1138. [Google Scholar] [CrossRef]
  102. Parr, R.G.; Yang, W. Density-Functional Theory of Atoms and Molecules; Oxford University Press: New York, NY, USA, 1989. [Google Scholar]
  103. Krishnan, R.; Binkley, J.S.; Seeger, R.; Pople, J.A. Self-consistent molecular orbital methods. XX. A basis set for correlated wave functions. J. Chem. Phys. 1980, 72, 650–654. [Google Scholar] [CrossRef]
  104. McLean, A.D.; Chandler, G.S. Contracted Gaussian basis sets for molecular calculations. I. second row atoms, Z=11–18. J. Chem. Phys. 1980, 72, 5639–5648. [Google Scholar] [CrossRef]
  105. Curtiss, L.A.; McGrath, M.P.; Blandeau, J.-P.; Davis, N.E.; Binning, R.C., Jr.; Radom, L. Extension of Gaussian-2 theory to molecules containing third-row atoms Ga–Kr. J. Chem. Phys. 1995, 103, 6104–6113. [Google Scholar] [CrossRef]
  106. Frisch, M.J.; Pople, J.A.; Binkley, J.S. Self-consistent molecular orbital methods 25. Supplementary functions for Gaussian basis sets. J. Chem. Phys. 1984, 80, 3265–3269. [Google Scholar] [CrossRef]
  107. Clark, T.; Chandrasekhar, J.; Spitznagel, G.W.; Schleyer, P.v.R. Efficient Diffuse Function-Augmented Basis Sets for Anion Calculations. III. The 3-21+G Basis Set for First-Row Elements, Li–F. J. Comput. Chem. 1983, 4, 294–301. [Google Scholar] [CrossRef]
  108. Mardirossian, N.; Head-Gordon, M. Thirty years of density functional theory in computational chemistry: An overview and extensive assessment of 200 density functionals. Mol. Phys. 2017, 19, 2315–2372. [Google Scholar] [CrossRef]
  109. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16, Revision C.01; Gaussian, Inc.: Wallingford, CT, USA, 2019. [Google Scholar]
  110. Bader, R.F.W. Atoms in Molecules: A Quantum Theory; Oxford University Press: New York, NY, USA, 1990. [Google Scholar]
  111. Popelier, P.L.A. Atoms in Molecules. An Introduction; Longman: Singapore, 2000. [Google Scholar]
  112. Matta, C.F.; Boyd, R.J. The Quantum Theory of Atoms in Molecules; Wiley-VCH: Weinheim, Germany, 2007. [Google Scholar]
  113. Keith, T.A. AIMAll (Version 15.05.18); TK Gristmill Software: Overland Park, KS, USA, 2015; Available online: aim.tkgristmill.com (accessed on 26 July 2022).
  114. Gatti, C. Chemical bonding in crystals: New directions. Z. Kristallogr. 2005, 220, 399–457. [Google Scholar] [CrossRef]
  115. Gatti, C.; Macchi, P. (Eds.) A Guided Tour Through Modern Charge-Density Analysis; Springer: Dordrecht, The Netherlands, 2011. [Google Scholar]
  116. Emamian, S.; Lu, T.; Kruse, H.; Emamian, H. Exploring Nature and Predicting Strength of Hydrogen Bonds: A Correlation Analysis Between Atoms-in-Molecules Descriptors, Binding Energies, and Energy Components of Symmetry-Adapted Perturbation Theory. J. Comput. Chem. 2019, 40, 2868–2881. [Google Scholar] [CrossRef]
  117. Groom, C.R.; Bruno, I.J.; Lightfoot, M.P.; Ward, S.C. The Cambridge Structural Database. Acta Cryst. 2016, B72, 171–179. [Google Scholar] [CrossRef]
  118. Bläser, D.; Boese, R.; Göhner, M.; Herrmann, F.; Kuhn, N.; Ströbele, M. 2,3-Dihydro-1,3,4,5-tetraisopropylimidazol-2-yliden / 2,3-Dihydro-1,3,4,5-tetraisopropylimidazol-2-ylidene. Z. Naturforsch. B 2014, 69, 71–76. [Google Scholar] [CrossRef]
  119. Cowan, J.A.; Clyburne, J.A.C.; Davidson, M.G.; Harris, R.L.W.; Howard, J.A.K.; Küpper, P.; Leech, M.A.; Richards, S.P. On the Interaction between N-Heterocyclic Carbenes and Organic Acids: Structural Authentication of the First N–H⋯C Hydrogen Bond and Remarkably Short C–H⋯O Interactions. Angew. Chem., Int. Ed. 2002, 41, 1432–1434. [Google Scholar] [CrossRef]
  120. Kinney, Z.J.; Rheingold, A.L.; Protasiewicz, J.D. Preferential N–H⋯:C< hydrogen bonding involving ditopic NH-containing systems and N-heterocyclic carbenes. RSC Adv. 2020, 10, 42164–42171. [Google Scholar]
  121. Li, C.-Y.; Kuo, Y.-Y.; Tsai, J.-H.; Yap, G.P.A.; Ong, T.-G. Amine-Linked N-Heterocyclic Carbenes: The Importance of an Pendant Free-Amine Auxiliary in Assisting the Catalytic Reaction. Chem. Asian J. 2011, 6, 1520–1524. [Google Scholar] [CrossRef]
  122. Movassaghi, M.; Schmidt, M.A. N-Heterocyclic Carbene-Catalyzed Amidation of Unactivated Esters with Amino Alcohols. Org. Lett. 2005, 7, 2453–2456. [Google Scholar] [CrossRef] [PubMed]
  123. Srivastava, R.; Moneuse, R.; Petit, J.; Pavard, P.-A.; Dardun, V.; Rivat, M.; Schiltz, P.; Solari, M.; Jeanneau, E.; Veyre, L.; et al. Early/Late Heterobimetallic Tantalum/Rhodium Species Assembled Through a Novel Bifunctional NHC-OH Ligand. Chem. Eur. J. 2018, 24, 4361–4370. [Google Scholar] [CrossRef]
  124. Guo, R.; Huang, X.; Zhao, M.; Lei, Y.; Ke, Z.; Kong, L. Bifurcated Hydrogen-Bond-Stabilized Boron Analogues of Carboxylic Acids. Inorg. Chem. 2019, 58, 13370–13375. [Google Scholar] [CrossRef] [PubMed]
  125. Jabłoński, M. Does the Presence of a Bond Path Really Mean Interatomic Stabilization? The Case of the Ng@Superphane (Ng = He, Ne, Ar, and Kr) Endohedral Complexes. Symmetry 2021, 13, 2241. [Google Scholar] [CrossRef]
  126. Jabłoński, M. Endo- and exohedral complexes of superphane with cations. J. Comput. Chem. 2022, 43, 1120–1133. [Google Scholar] [CrossRef] [PubMed]
  127. Jabłoński, M. The physical nature of the ultrashort spike–ring interaction in iron maiden molecules. J. Comput. Chem. 2022, 43, 1206–1220. [Google Scholar] [CrossRef] [PubMed]
  128. Mallinson, P.R.; Smith, G.T.; Wilson, C.C.; Grech, E.; Wozniak, K. From Weak Interactions to Covalent Bonds: A Continuum in the Complexes of 1,8-Bis(dimethylamino)naphthalene. J. Am. Chem. Soc. 2003, 125, 4259–4270. [Google Scholar] [CrossRef]
  129. Dominiak, P.M.; Makal, A.; Mallinson, P.R.; Trzcinska, K.; Eilmes, J.; Grech, E.; Chruszcz, M.; Minor, W.; Woźniak, K. Continua of Interactions between Pairs of Atoms in Molecular Crystals. Chem. Eur. J. 2006, 12, 1941–1949. [Google Scholar] [CrossRef]
  130. Jabłoński, M.; Palusiak, M. Basis Set and Method Dependence in Atoms in Molecules Calculations. J. Phys. Chem. A 2010, 114, 2240–2244. [Google Scholar] [CrossRef]
  131. Jabłoński, M.; Solà, M. Influence of Confinement on Hydrogen Bond Energy. The Case of the FH⋯NCH Dimer. J. Phys. Chem. A 2010, 114, 10253–10260. [Google Scholar] [CrossRef]
Figure 1. General scheme of the IR 2 HD dimers (R = H, Me, i Pr, t Bu, Ph, Mes, Dipp, Ad). The colon on the carbene carbon atom represents a lone electron pair.
Figure 1. General scheme of the IR 2 HD dimers (R = H, Me, i Pr, t Bu, Ph, Mes, Dipp, Ad). The colon on the carbene carbon atom represents a lone electron pair.
Molecules 27 05712 g001
Figure 2. Examples of the occurrence of hydrogen bonds (marked with a dashed green line) involving a carbene carbon atom in crystallographic structures. The atoms are labeled as follows: carbon—gray, hydrogen—white, nitrogen—blue, oxygen—red, and boron—pink.
Figure 2. Examples of the occurrence of hydrogen bonds (marked with a dashed green line) involving a carbene carbon atom in crystallographic structures. The atoms are labeled as follows: carbon—gray, hydrogen—white, nitrogen—blue, oxygen—red, and boron—pink.
Molecules 27 05712 g002
Figure 3. Map of the electrostatic potential (in a.u.) of imidazol-2-ylidene: −0.06, red; −0.03, yellow; 0.00, green; 0.03, cyan; 0.06, blue.
Figure 3. Map of the electrostatic potential (in a.u.) of imidazol-2-ylidene: −0.06, red; −0.03, yellow; 0.00, green; 0.03, cyan; 0.06, blue.
Molecules 27 05712 g003
Figure 4. Imidazol-2-ylidene dimers with HF, HCN, NH 3 , H 2 O and MeOH. The values of the dissociation energy are given in bold, while the three values next to the interaction shown refer to the length of the hydrogen bond, its angle and the value of the electron density at the bond critical point.
Figure 4. Imidazol-2-ylidene dimers with HF, HCN, NH 3 , H 2 O and MeOH. The values of the dissociation energy are given in bold, while the three values next to the interaction shown refer to the length of the hydrogen bond, its angle and the value of the electron density at the bond critical point.
Molecules 27 05712 g004
Figure 5. The IR 2 HF and IR 2 HCN (R = Me, i Pr, t Bu, Ph, Mes, Dipp, Ad) dimers. The values of the dissociation energy are given in bold, while the three values refer to the length of the hydrogen bond, the C-H-D angle and the value of the electron density at the bond critical point of the C⋯H-D hydrogen bond.
Figure 5. The IR 2 HF and IR 2 HCN (R = Me, i Pr, t Bu, Ph, Mes, Dipp, Ad) dimers. The values of the dissociation energy are given in bold, while the three values refer to the length of the hydrogen bond, the C-H-D angle and the value of the electron density at the bond critical point of the C⋯H-D hydrogen bond.
Molecules 27 05712 g005
Figure 6. Molecular graphs of selected IR 2 HF and IR 2 HCN dimers. Arrows show the presence of bond paths for some accompanying interatomic interactions (hydrogen bonds—red, tetrel bonds—brown, and C-H⋯H-C contacts—yellow). Large balls represent atoms (hydrogen—white, carbon—gray, nitrogen—blue, and fluorine—light green) and small balls represent critical points (bond critical points—light green, ring critical points—red, and cage critical points—blue).
Figure 6. Molecular graphs of selected IR 2 HF and IR 2 HCN dimers. Arrows show the presence of bond paths for some accompanying interatomic interactions (hydrogen bonds—red, tetrel bonds—brown, and C-H⋯H-C contacts—yellow). Large balls represent atoms (hydrogen—white, carbon—gray, nitrogen—blue, and fluorine—light green) and small balls represent critical points (bond critical points—light green, ring critical points—red, and cage critical points—blue).
Molecules 27 05712 g006
Figure 7. The IR 2 H 2 O and IR 2 MeOH (R = Me, i Pr, t Bu, Ph, Mes, Dipp, Ad) dimers. The values of the dissociation energy are given in bold, while the three values refer to the length of the hydrogen bond, the C-H-D angle and the value of the electron density at the bond critical point of the C⋯H-D hydrogen bond.
Figure 7. The IR 2 H 2 O and IR 2 MeOH (R = Me, i Pr, t Bu, Ph, Mes, Dipp, Ad) dimers. The values of the dissociation energy are given in bold, while the three values refer to the length of the hydrogen bond, the C-H-D angle and the value of the electron density at the bond critical point of the C⋯H-D hydrogen bond.
Molecules 27 05712 g007
Figure 8. Molecular graphs of selected IR 2 H 2 O and IR 2 MeOH dimers. Arrows show the presence of bond paths for some accompanying interatomic interactions (hydrogen bonds—red, C-H⋯H-C contacts—yellow, and C-H π hydrogen bonds—light green). Large balls represent atoms (hydrogen—white, carbon—gray, nitrogen—blue, and oxygen—red) and small balls represent critical points (bond critical points—light green, ring critical points—red, and cage critical points—blue).
Figure 8. Molecular graphs of selected IR 2 H 2 O and IR 2 MeOH dimers. Arrows show the presence of bond paths for some accompanying interatomic interactions (hydrogen bonds—red, C-H⋯H-C contacts—yellow, and C-H π hydrogen bonds—light green). Large balls represent atoms (hydrogen—white, carbon—gray, nitrogen—blue, and oxygen—red) and small balls represent critical points (bond critical points—light green, ring critical points—red, and cage critical points—blue).
Molecules 27 05712 g008
Figure 9. The IR 2 NH 3 (R = Me, i Pr, t Bu, Ph, Mes, Dipp, and Ad) dimers. The values of the dissociation energy are given in bold, while the three values refer to the length of the hydrogen bond, the C-H-D angle and the value of the electron density at the bond critical point of the C⋯H-D hydrogen bond.
Figure 9. The IR 2 NH 3 (R = Me, i Pr, t Bu, Ph, Mes, Dipp, and Ad) dimers. The values of the dissociation energy are given in bold, while the three values refer to the length of the hydrogen bond, the C-H-D angle and the value of the electron density at the bond critical point of the C⋯H-D hydrogen bond.
Molecules 27 05712 g009
Figure 10. Molecular graphs of the IR 2 NH 3 dimers. Arrows show the presence of bond paths for some accompanying interatomic interactions (hydrogen bonds—red, C-H⋯H-C or C-H⋯H-N contacts—yellow, and N π contact—dark green). Large balls represent atoms (hydrogen—white, carbon—gray, and nitrogen—blue) and small balls represent critical points (bond critical points—light green, ring critical points—red, and cage critical points—blue).
Figure 10. Molecular graphs of the IR 2 NH 3 dimers. Arrows show the presence of bond paths for some accompanying interatomic interactions (hydrogen bonds—red, C-H⋯H-C or C-H⋯H-N contacts—yellow, and N π contact—dark green). Large balls represent atoms (hydrogen—white, carbon—gray, and nitrogen—blue) and small balls represent critical points (bond critical points—light green, ring critical points—red, and cage critical points—blue).
Molecules 27 05712 g010
Figure 11. The dependence of either the dissociation energy (left) or the binding energy according to Equation (1) (right) on the distance C⋯H in the IR 2 HD dimers.
Figure 11. The dependence of either the dissociation energy (left) or the binding energy according to Equation (1) (right) on the distance C⋯H in the IR 2 HD dimers.
Molecules 27 05712 g011
Table 1. Values of selected parameters characterizing the C⋯H-D hydrogen bond or the proton-donor molecule in the IR 2 HD dimers (I = imidazol-2-ylidene; R = Me, i Pr, t Bu, Ph, Mes, Dipp, Ad; HD = HF, HCN, H 2 O, MeOH, NH 3 ): D 0 –dissociation energy (in kcal/mol), ρ C H –the electron density at the bond critical point of C⋯H-D (in a.u.), E b –binding energy of the C⋯H-D hydrogen bond obtained by Equation (1) (in kcal/mol), d C H –length of the C⋯H distance (in Å), α CHD –the C-H-D angle (in degrees), Δ d HD –elongation of the H-D proton-donor bond (in Å) relative to isolated HD (HF–0.917 Å, HCN–1.068 Å, H 2 O–0.958 Å, MeOH–0.957 Å, and NH 3 –1.013 Å).
Table 1. Values of selected parameters characterizing the C⋯H-D hydrogen bond or the proton-donor molecule in the IR 2 HD dimers (I = imidazol-2-ylidene; R = Me, i Pr, t Bu, Ph, Mes, Dipp, Ad; HD = HF, HCN, H 2 O, MeOH, NH 3 ): D 0 –dissociation energy (in kcal/mol), ρ C H –the electron density at the bond critical point of C⋯H-D (in a.u.), E b –binding energy of the C⋯H-D hydrogen bond obtained by Equation (1) (in kcal/mol), d C H –length of the C⋯H distance (in Å), α CHD –the C-H-D angle (in degrees), Δ d HD –elongation of the H-D proton-donor bond (in Å) relative to isolated HD (HF–0.917 Å, HCN–1.068 Å, H 2 O–0.958 Å, MeOH–0.957 Å, and NH 3 –1.013 Å).
HD Property Carbene
I IMe 2 I i Pr 2 I t Bu 2 IPh 2 IMes 2 IDipp 2 IAd 2
HF D 0 16.118.119.018.317.518.919.919.2
ρ C H 0.0650.0740.0770.0750.0730.0760.0790.076
E b −13.7−15.7−16.3−15.9−15.5−16.1−16.8−16.3
d C H 1.6761.6261.6121.6261.6321.6141.5991.618
α CHD 180.0180.0180.0179.4171.5169.8180.0180.0
Δ d HD 0.0610.0750.0800.0840.0750.0800.0830.088
HCN D 0 8.69.49.910.09.811.011.411.3
ρ C H 0.0250.0270.0280.0290.0270.0290.0300.030
E b −4.8−5.3−5.5−5.8−5.2−5.6−5.9−6.0
d C H 2.1382.1052.0872.0712.1032.0702.0582.055
α CHD 180.0180.0180.0180.0173.7167.2179.7180.0
Δ d HD 0.0290.0330.0350.0380.0310.0340.0350.040
H 2 O D 0 8.49.09.510.39.710.711.29.9
ρ C H 0.0270.0340.0350.0400.0350.0360.0330.033
E b −5.3−6.8−7.0−8.3−7.1−7.3−6.7−6.6
d C H 2.0351.9681.9591.9011.9541.9301.9761.991
α CHD 140.4163.8166.4174.4166.6163.4168.5172.9
Δ d HD 0.0190.0250.0260.0330.0260.0280.0230.027
MeOH D 0 9.310.411.212.211.513.114.412.2
ρ C H 0.0280.0360.0370.0410.0380.0370.0380.038
E b −5.5−7.2−7.4−8.4−7.6−7.6−7.7−7.6
d C H 2.0351.9561.9451.8981.9271.9231.9271.931
α CHD 140.8162.1163.6173.9172.3165.2173.2176.5
Δ d HD 0.0180.0250.0260.0320.0260.0260.0260.030
NH 3 D 0 7.25.15.34.55.26.77.35.4
ρ C H 0.0110.0180.0190.0180.0190.0200.0180.018
E b −1.6−3.3−3.4−3.2−3.5−3.7−3.2−3.2
d C H 2.5212.2802.2762.3142.2712.2222.2912.307
α CHD 123.1157.0158.2173.6165.9158.5161.1174.4
Δ d HD 0.0070.0110.0110.0110.0100.0120.0090.011
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jabłoński, M. On the Coexistence of the Carbene⋯H-D Hydrogen Bond and Other Accompanying Interactions in Forty Dimers of N-Heterocyclic-Carbenes (I, IMe2, IiPr2, ItBu2, IMes2, IDipp2, IAd2; I = imidazol-2-ylidene) and Some Fundamental Proton Donors (HF, HCN, H2O, MeOH, NH3). Molecules 2022, 27, 5712. https://doi.org/10.3390/molecules27175712

AMA Style

Jabłoński M. On the Coexistence of the Carbene⋯H-D Hydrogen Bond and Other Accompanying Interactions in Forty Dimers of N-Heterocyclic-Carbenes (I, IMe2, IiPr2, ItBu2, IMes2, IDipp2, IAd2; I = imidazol-2-ylidene) and Some Fundamental Proton Donors (HF, HCN, H2O, MeOH, NH3). Molecules. 2022; 27(17):5712. https://doi.org/10.3390/molecules27175712

Chicago/Turabian Style

Jabłoński, Mirosław. 2022. "On the Coexistence of the Carbene⋯H-D Hydrogen Bond and Other Accompanying Interactions in Forty Dimers of N-Heterocyclic-Carbenes (I, IMe2, IiPr2, ItBu2, IMes2, IDipp2, IAd2; I = imidazol-2-ylidene) and Some Fundamental Proton Donors (HF, HCN, H2O, MeOH, NH3)" Molecules 27, no. 17: 5712. https://doi.org/10.3390/molecules27175712

APA Style

Jabłoński, M. (2022). On the Coexistence of the Carbene⋯H-D Hydrogen Bond and Other Accompanying Interactions in Forty Dimers of N-Heterocyclic-Carbenes (I, IMe2, IiPr2, ItBu2, IMes2, IDipp2, IAd2; I = imidazol-2-ylidene) and Some Fundamental Proton Donors (HF, HCN, H2O, MeOH, NH3). Molecules, 27(17), 5712. https://doi.org/10.3390/molecules27175712

Article Metrics

Back to TopTop