Next Article in Journal
Oxyethylated Fluoresceine—(thia)calix[4]arene Conjugates: Synthesis and Visible-Light Photoredox Catalysis in Water–Organic Media
Next Article in Special Issue
The Halogen Bond in Weakly Bonded Complexes and the Consequences for Aromaticity and Spin-Orbit Coupling
Previous Article in Journal
A Novel Curcumin Arginine Salt: A Solution for Poor Solubility and Potential Anticancer Activities
Previous Article in Special Issue
Dissecting Bonding Interactions in Cysteine Dimers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dihydrogen Bonding—Seen through the Eyes of Vibrational Spectroscopy

Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, TX 75275-0314, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Molecules 2023, 28(1), 263; https://doi.org/10.3390/molecules28010263
Submission received: 6 December 2022 / Revised: 20 December 2022 / Accepted: 20 December 2022 / Published: 28 December 2022
(This article belongs to the Special Issue Chemical Bond and Intermolecular Interactions)

Abstract

:
In this work, we analyzed five groups of different dihydrogen bonding interactions and hydrogen clusters with an H 3 + kernel utilizing the local vibrational mode theory, developed by our group, complemented with the Quantum Theory of Atoms–in–Molecules analysis to assess the strength and nature of the dihydrogen bonds in these systems. We could show that the intrinsic strength of the dihydrogen bonds investigated is primarily related to the protonic bond as opposed to the hydridic bond; thus, this should be the region of focus when designing dihydrogen bonded complexes with a particular strength. We could also show that the popular discussion of the blue/red shifts of dihydrogen bonding based on the normal mode frequencies is hampered from mode–mode coupling and that a blue/red shift discussion based on local mode frequencies is more meaningful. Based on the bond analysis of the H 3 + (H 2 ) n systems, we conclude that the bond strength in these crystal–like structures makes them interesting for potential hydrogen storage applications.

1. Introduction

Hydrogen bonding (HB) is one of the most prominent non–covalent chemical interactions [1,2,3,4,5], being responsible for the shape and structure of biopolymers found in nature, such as proteins and nucleic acids [6,7,8], but also being a major player in transition metal catalysis [9], solid–state chemistry and materials science [10]. Generally, HB occurs between a positively charged hydrogen atom of a proton donor, and a lone pair of an electronegative atom, π electron cloud, or a transition metal center, representing the corresponding proton acceptor. However, in recent years, in addition to this conventional HB interaction, another intriguing HB type has been discussed, where the two atoms directly involved in the weak interaction are both hydrogen atoms [11,12,13,14]. This interaction was coined by Richardson [15] dihydrogen bonding (DHB will be used for dihydrogen bonds in this study). This describes the attraction between a protonic (H δ + ) and a hydridic (H δ ) hydrogen atom, where the protonic H atom acts as the proton donor as in conventional HBs, and the hydridic H atom takes over the role of the proton acceptor [11].
Interaction energies of DHBs are generally between 1 and 7 kcal/mol, and because of their similarity to HBs, they are often treated as a subclass of HBs [16]. It has been suggested that DHBs are predominantly of electrostatic nature [16,17,18,19], accompanied in stronger DHBs with a substantial covalent contribution [16,20,21,22]; although, this has not been quantified so far. DHBs can influence molecular properties in the gas phase, in solution and in the solid state, yielding broad potential utilities in catalysis and materials sciences. In the solid state, DHBs show the unique ability to lose H 2 , which could be utilized for the design of low–cost, high–capacity hydrogen storage materials [23,24,25]. In particular, a better understanding of the transformation from a weak H−H electrostatic interaction to a strong covalent bond of H 2 type could open new routes for technologies using hydrogen as an environmentally clean and efficient fuel.
The first experimental observation of the attractive interaction between two hydrogen atoms was reported in early 1930 by Zachariasen and Mooney [26] in the crystal structure of ammonium hypophosphite ( NH 4 H 2 PO 2 ), showing a close contact between the H atoms of the hypophosphite group and those of ammonium. More than 30 years later, Burg [27] reported the direct intermolecular interaction between the H atoms of the NH and BH 3 groups in liquid ( CH 3 ) 2 NH HBH 2 , suggesting the presence of an H−H interaction comparable to that of a HB. This interaction was formally categorized as a true HB by Brown in the late 1960s [28,29,30] on the basis of the IR spectra of boron coordination compounds.
DHB is particularly attractive for molecular systems involving metal hydrides acting as proton acceptors. The first experimental evidence for the formation of a DBH involving a metal hydride was presented independently by Crabtree [31] and Morris [32]. Later, Crabtree et al. [33,34] investigated the H−H interaction in molecular complexes involving iridium hydrides, exploring the strength of this interaction with experimental and computational means. By using the Cambridge Crystallographic Database [35] and spectroscopic data of aminoboranes, Richardson et al. [15] and later Crabtree et at. [36] showed that the NH HB distances and heats of formation are in the ranges of 1.7 2.2 Å and 3–7 kcal/mol, respectively, comparable to that of conventional HBs.
Over the years, DHB has been also the subject of many theoretical studies utilizing different levels of theoretical approaches. Comprehensive overviews can be found in Refs. [1,2,3,14,37,38] and studies cited therein. Furthermore, two new data sets have been recently reported for benchmarking [39]. The focus of these theoretical studies ranges from energy decomposition and bond energies [40,41,42], geometries [1,11], normal mode frequency shifts [9,43], studies of the topological features of the total electron density ρ ( r ) [44,45,46,47], and electron localization function (ELF) investigations [48,49] to the analysis of electrostatic potential [50,51,52], which led to a number of open questions. In particular, what is missing so far is a quantitative bond strength measure coupled with the assessment of the covalent character, which allows the classification of DHB as (i) weak H⋯H interactions of the van der Waals type, (ii) moderate H⋯H interactions mostly of electrostatic type, and (iii) strong H−H interactions with substantial covalent character.
Therefore, the focus of our study was to introduce a quantitative bond strength measure based on the local vibrational mode analysis (LMA), originally developed by Konkoli and Cremer [53,54,55,56,57], paired with Weinhold’s Natural Bond Orbital (NBO) population analysis [58,59,60], and Bader’s Quantum Theory of Atoms–In–Molecules (QTAIM) analysis [61,62,63,64], in order to obtain a deeper understanding of DHBs in general, in particular their strength, and to unravel similarities/differences compared to conventional HBs. Figure 1 shows the set of DHB complexes investigated in this study (organized in six specific groups) covering a broad spectrum of DHBs, ranging from weak interactions of electrostatic origin to strong interactions with covalent character. Group I includes hydrides interacting with simple organic neutral and cationic molecules as well intramolecular DHB in boron complexes, and Group II shows DHB in ammonia–borane complexes. DHB in Fe complexes are included in Group III, while DHB in noble gas complexes are involved in Group IV. We also reassessed the H−H interaction in aromatic hydrocarbons such as phenanthrene, dibenz[a]anthracene, and biphenyl [22,65] in Group V, where both hydrogen atoms are equally charged, which has been called the H−H bond in the literature, to contrast this situation from DHB [11], leading to some controversy [66]. Group VI of our study is made up of DHB in H3+(H 2 ) n hydrogen clusters.
The manuscript is structured in the following way: In the results and discussion part, general findings for Group I–Group V are summarized, followed by a more specific discussion of each individual group and the description of H–H interactions in H3+(H 2 ) n clusters (Group VI, Figure 1). In the computational methods section, a short summary of LMA is given, and the model chemistry applied and software utilized are described. In the final part, the conclusions and future aspects of this work are presented.

2. Results

2.1. Overall Trends

In the following, general trends found for Group I–Group VI members are discussed referring to the data collected in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9 and Table 1, Table 2 and Table 3. In addition, bond lengths R in Å, local mode force constants k a in mdyn/Å, local mode vibrational frequencies ω a in cm 1 , binding energies E in kcal/mol, electron densities at the bond critical points ρ c in e/Å 3 and energy densities at bond critical points H c in Hartree/Å 3 of protonic and hydridic parts of DHBs and DHBs for all complexes and references compounds are compiled in Table S2 of the Supplementary Materials.
Strength of DHBs:Figure 2a–c shows BSO n values for the DHBs, protonic and hydridic parts of DHBs; in Figure 2d the relationship between protonic and hydridic local mode force constants k a is shown; and in Figure 2e,f the pairwise correlation between local mode force constants k a (Protonic), k a (Hydridic) and k a (DHB) is shown. BSO n (DHB) values range from 0.14 (FeF2) to 0.41 (BeO4), which is in a similar range as we previously found for HBs (BSO n = 0.135, N–H⋯F; BSO n = 0.33, N–H⋯N; see Ref. [67]), quantifying that they are both of comparable strength with the former being somewhat stronger. The protonic parts of DHBs ranging from 0.49 for BeO3 to 0.99 for CF (being close to the FH reference value of BSO n = 1) are stronger than their hydridic counter parts with BSO n values of 0.41 for NaO and 0.86 for Ar3, respectively. One could argue that a weaker protonic part of DHB caused by charge transfer to the hydridic bond should lead to the strengthening of the latter and vice versus, resulting in a correlation between the strengths of protonic and hydridic parts of DHBs. As is obvious from Figure 2d, there is no such correlation between k a (Protonic) and k a (Hydridic). On the contrary, according to Figure 2d, there is some trend that weaker protonic parts of DHBs are paired with weaker hydridic bonds, as in the case of BeO3 and that stronger protonic parts of DHBs are paired with stronger hydridic parts of DHBs, as in the case of CF. According to Figure 2e, there is also a tendency that weaker protonic parts of DHBs correspond to stronger DHBs, e.g., as found for BeO3, whereas stronger protonic parts of DHBs correspond to weaker DHBs, indicating an larger electron density transfer from a weak protonic part of DHB to DHB, that the donor bond transfers electron density to the DHB. However, as depicted in Figure 2f, this tendency is less pronounced for the hydridic parts of DHB.
Overall, these findings suggest that if one wants to design a DHB with a certain strength, the focus should be on modifying and fine–tuning the protonic rather than the hydridic part of DHB guided by the local mode force constants as a sensitive tool monitoring the changes in the bond strengths.
Figure 2. (a) BSO n values for DHB, (b) BSO n values for protonic parts of DHBs, and (c) BSO n values for hydridic parts of DHBs calculated from the corresponding local mode force constants as described in the text. (d) Correlation between k a (Protonic) and k a (Hydridic), (e) between k a (Protonic) and k a (DHB), (f) between k a (Hydridic) and k a (DHB).
Figure 2. (a) BSO n values for DHB, (b) BSO n values for protonic parts of DHBs, and (c) BSO n values for hydridic parts of DHBs calculated from the corresponding local mode force constants as described in the text. (d) Correlation between k a (Protonic) and k a (Hydridic), (e) between k a (Protonic) and k a (DHB), (f) between k a (Hydridic) and k a (DHB).
Molecules 28 00263 g002
Figure 3. (a) Correlation of k a (DHB) and normalized energy density H c / ρ c . (b) Correlation of charge differences Δ q (see Table 1 for definition) and k a (DHB).
Figure 3. (a) Correlation of k a (DHB) and normalized energy density H c / ρ c . (b) Correlation of charge differences Δ q (see Table 1 for definition) and k a (DHB).
Molecules 28 00263 g003
Covalency, protonic and hydridic charges: The covalency of the DHBs investigated in this work was assessed by the normalized energy density H c / ρ c (see Table S2 in the Supplementary Materials). In Figure 3a k a (DHB) is correlated with H c / ρ c which ranges from covalent H c / ρ c = −0.6 Hartree/electron, to electrostatic H c / ρ c = 0.2 Hartree/electron in accordance with the range of H c / ρ c values −0.6 to 0.3 Hartree/electron for the HBs investigated in our previous work [67], disproving the assumption that DHBs are generally more electrostatic than their HB counterparts. Although there is no strong correlation (R 2 = 0.860) between these two quantities, Figure 3a shows that stronger DHBs such as BeO4 or BeN2 are of covalent character whereas weak BDHs such as BN3 or BC2 are electrostatic in nature. In summary, the Cremer–Kraka criterion of covalent bond offers an efficient tool to assess the covalency of DHB interactions, which as our study shows can be quite different, depending on the electronic structure of the DHB complex and the nature of the protonic and hydridic parts of DHB.
In Table 1 and Table 2 protonic and hydridic H–atom NBO charges and their differences for Group I–V complexes are reported, and in Table S1 of Supplementary Materials the NBO hydrogen charges for the reference compounds are shown.
Table 1. Protonic and hydridic NBO charges (in e) for Group I. Protonic hydrogen charge q(H δ + ); hydridic hydrogen charge q(H δ ); charge difference Δ q = q(H δ + ) − q(H δ ).
Table 1. Protonic and hydridic NBO charges (in e) for Group I. Protonic hydrogen charge q(H δ + ); hydridic hydrogen charge q(H δ ); charge difference Δ q = q(H δ + ) − q(H δ ).
Complexq(H δ + )q(H δ ) Δ q Complexq(H δ + )q(H δ ) Δ q Complexq(H δ + )q(H δ ) Δ q
LiC10.236−0.8041.040BeN10.467−0.5561.024BN50.404−0.0780.482
LiC20.213−0.7690.982BeN20.418−0.5300.948BO10.494−0.1040.598
LiC30.189−0.7760.964BeN30.417−0.5520.969BO20.491−0.1040.595
LiC40.249−0.8151.064BeN40.439−0.5000.939BS10.166−0.0830.249
LiC50.258−0.8161.074BeO10.505−0.5371.042BS20.166−0.0920.258
LiC60.269−0.8161.086BeO20.521−0.5151.036BF0.565−0.1620.727
LiN10.326−0.7401.066BeO30.497−0.4670.964SiN0.462−0.3100.772
LiN20.428−0.7691.197BeO40.558−0.4841.042SiO0.534−0.2920.826
LiO0.493−0.7291.222BeF0.558−0.4841.042SiF0.552−0.2290.781
NaC10.238−0.8091.048BeCl0.275−0.4750.750AlO10.504−0.4450.949
NaC20.214−0.8141.028BC10.227−0.0790.306AlO20.513−0.4490.962
NaC30.190−0.8181.008BC20.218−0.0880.306AlF0.592−0.4171.009
NaC40.246−0.8121.059BC30.262−0.1400.402GaF0.575−0.3230.898
NaC50.258−0.8171.075 0.262−0.0630.325CF0.5480.1870.362
NaC60.271−0.8191.090BC40.259−0.0870.346 0.5480.2100.339
NaN10.313−0.7561.068 0.259−0.0710.330
NaO0.480−0.7591.239BN10.380−0.0930.474
BeC10.229−0.4840.713BN20.461−0.0920.554
BeC30.244−0.4650.709BN40.466−0.0850.551
As obvious from these data, the complexes investigated in this work display a large range of hydrogen charges reflecting the breadth of the different DHB scenarios covered. Protonic H–atom charges range from 0.166 e in BS1 and BS2 to 0.592 e in AlF (see Table 1) and hydridic H–atom charges from −0.061 e in AB2 (see Table 2) to −0.819 e in NaC6 (see Table 1). It is noteworthy that some of the hydridic charges in transition metal and aromatic compounds are positive, which will be further elucidated below. Figure 3b shows the correlation between Δ q and the DHB bond strength as reflected by k a (DHB). As obvious from the large scattering of data points, there is no direct correlation between these two quantities, revealing that differences in the atomic charges is only one of the components determining the strength of these complex interaction. Therefore, Δ q , although convenient, cannot serve as a direct DHB strength measure.
Table 2. Protonic and hydridic NBO charges (in e) for Group II–V. Protonic hydrogen charge q(H δ + ); hydridic hydrogen charge q(H δ ); charge difference Δ q = q(H δ + ) − q(H δ ).
Table 2. Protonic and hydridic NBO charges (in e) for Group II–V. Protonic hydrogen charge q(H δ + ); hydridic hydrogen charge q(H δ ); charge difference Δ q = q(H δ + ) − q(H δ ).
GroupComplexq(H δ + )q(H δ ) Δ q GroupComplexq(H δ + )q(H δ ) Δ q
AminoboranesAB10.448−0.0650.512TransitionFeN10.3700.1060.264
0.448−0.0650.512metalsFeN20.3500.1320.218
AB20.445−0.0610.506 FeO10.4830.0740.409
0.445−0.0640.509 FeO20.4900.1330.357
0.441(H1)(H2) − 0.0520.493 FeF10.5670.0470.519
0.441(H1)(H3) − 0.0620.503 FeF20.5920.1300.462
AB30.443−0.0580.502Noble gasesKrO0.464−0.3000.765
0.435(H1)(H2,H3) − 0.0560.491 BeAr0.270−0.5100.781
AB40.296−0.0630.359 BeKr0.147−0.5020.649
0.296−0.0630.359AromaticAr10.2040.2040.000
AB50.494−0.0740.568compoundsAr20.2040.1930.011
AB60.388−0.0650.453 Ar30.2020.2020.000
AB70.571−0.0880.659
DHB binding energies, bond lengths and bond angles and local mode force constants:Figure 4 shows how BEs and geometrical features that are frequently used to characterize the strength of a DHB, correlate with the local mode force constants k a (DHB). Because BE includes geometry relaxation and electron density reorganization effects of the fragments upon DHB dissociation, as discussed above, there is no significant correlation between BE and k a (DHB) as expected (R 2 = 0.6186; see Figure 4a), in particular not for the broad range of complexes investigated in this work, leading to a variety of different fragments formed during the breakage of the DHB and disassociation into fragments. Nevertheless, we observed some trends: e.g., larger BEs are associated with larger k a (DHB) values, reflecting stronger DHBs for the BeO complexes, whereas weak DHBs are characterized by small BE and k a (DHB) values, such as LiC1 or LiC2.
Figure 4b shows an overall power relationship between DHB lengths R(DHB) and k a (DHB) following the generalized Badger rule of Cremer and co–workers [68]. Again, no significant correlation (R 2 = 0.7014) is observed despite some trends; one of the strongest DHBs found in this study, namely for complexes BeO1, BeO2, and BeO3 are characterized by smaller R(DHB) values (1.237, 1.170 and 1.091 Å, respectively), whereas the weakest DHBs found e.g., for LiC1 and BC4 have the longer DHBs with values of 2.495 Å and 2.533 Å, respectively. However, there are also outliers, such as FeF2 which also has one of the weakest DHBs but a DHB length of 1.802 Å. Also obvious is that the Group V members do not follow this generalized Badger type rule caused by the specific topology of the bay H atoms forming these interactions. A more detailed discussion if these interaction can be considered as DHBs, can be found below.
The overall shortest DHB within Group I–IV is R(DHB) 1.091 Å (BeO3), which is still 0.348 Å longer than the R(HH) distance of 0.743 Å of the hydrogen molecule H2. Therefore, in our set of complexes we did not find any realization of a X δ + H–H⋯Y δ situation, suggested by Bakhmutov [11].
Figure 4c,d check the hypothesis that the angle X–H δ + H δ of the protonic fragment is closer to linear, whereas the hydridic fragment H δ + H δ –Y is more bent and if there is a correlation between these angles and the DHB bond strength. We do not find any such correlation reflected by the large scattering of data points in Figure 4c,d. Furthermore, we also cannot confirm that the protonic part of the complex is more linear, whereas the hydridic part is more bent.
Local mode versus normal mode frequency red/blue shifts: Frequently normal mode frequency red/blue shifts are utilized to characterize DHB and HB. Although the frequency shifts of the DHB or HB normal mode stretching frequencies can be used to detect these bonds [43,69,70,71], they are not necessarily suited to assess the bond strength of these interactions as we discussed for HBs [67]. The normal mode frequency shift Δ ω = ω ( c o m p l e x ) ω ( m o n o m e r ) can only be related to the thermochemical strength of the HB or DHB in question, if the HB or DHB normal stretching modes in both complex and monomer are localized, which is generally not the case [67,72]. Therefore, we draw upon to local mode frequency shifts instead.
Figure 4. (a) Correlation between binding energy E b i n scaled by the number of DHBs n and k a (DHB), (b) correlation between distance R(DHB) and k a (DHB), (c) correlation between protonic bond angle (X–H δ + H δ ) and k a (DHB), (d) correlation between hydridic bond angle (H δ + H δ –Y) and k a (DHB).
Figure 4. (a) Correlation between binding energy E b i n scaled by the number of DHBs n and k a (DHB), (b) correlation between distance R(DHB) and k a (DHB), (c) correlation between protonic bond angle (X–H δ + H δ ) and k a (DHB), (d) correlation between hydridic bond angle (H δ + H δ –Y) and k a (DHB).
Molecules 28 00263 g004
Figure 5a depicts the local mode frequency shift Δ ω of the protonic parts of DHB and Figure 5b that of the hydridic parts presented as a bar diagram. For both the DHB donor and acceptor, we found almost exclusively red shifts, i.e., the weakening of the bond in question upon complexation; however, the magnitude of the shifts differ significantly; protonic parts of DHBs shifts range from 139 to −2421 cm 1 , whereas hydridic parts shifts cover a considerably smaller range, namely from 18 to −406 cm 1 . It is interesting to note that largest red shifts reflecting large changes in DHB donor and acceptor upon DHB formation occur for BeO3 and BeO4, i.e., the two complexes with the strongest DHBs.
In the following, we focus on evaluating to what extent the normal vibrational modes are localized in the DHB donor and/or acceptor of these complexes, which forms the necessary prerequisite for a red/blue shift discussion based on normal vibrational modes. As assessment tool we used in this work CMN to probe the local character of both protonic and hydridic normal stretching modes in complex and references molecules. As an example the CMN for the complex BeO3 and its two references HOCL2 and BeH is shown in Figure 6a,b.
The normal vibrational stretching mode of the DHB protonic part of reference HOCL2 (O1H4, yellow bar in Figure 6a) with a frequency of 3589 cm 1 is 100% localized. In contrast, the hydridic normal vibrational stretching modes of reference BeH with frequencies of 2203 and 1995 cm 1 are composed of 50% Be2H3 (red color) and 50% Be2H1 (orange color), reflecting the symmetry of the BeH reference molecule. Upon complex formation, the protonic part of DHB is considerably weakened, and most importantly, as reflected in Figure 6a, it strongly couples with other modes in the BeO3 complex. In BeO3, the normal mode corresponding to the complex normal vibration of 1386 cm 1 has 47% O5H4 (the protonic part of DHB, yellow color) and 39% H3H4 (DHB, purple color); see Figure 6a. The normal vibrational modes corresponding to 2286 cm 1 and 2101 cm 1 have 14% and 58% B3H3 (the hydridic bond) character. This clearly shows (i) the strong influence of the protonic part of DHB on DHB and (ii) that a red/blue shift discussion based on normal vibrational mode frequencies is rather questionable, and one has to resort to local mode frequencies for a meaningful analysis, which is further elucidated in the following.
Figure 5. (a) Local mode frequency shifts Δ ω of the protonic parts of DHBs. (b) Local mode frequency shifts Δ ω of the hydridic parts of DHBs, where a positive sign denotes a blue–shift and a negative sign a red–shift.
Figure 5. (a) Local mode frequency shifts Δ ω of the protonic parts of DHBs. (b) Local mode frequency shifts Δ ω of the hydridic parts of DHBs, where a positive sign denotes a blue–shift and a negative sign a red–shift.
Molecules 28 00263 g005
Figure 7a presents the relationship between protonic and hydridic local mode frequency shifts Δ ω , and Figure 7b the relationship between the corresponding relative shifts Δ ω / ω , involving Group I–Group IV members.
Figure 6. (a) CNM analysis for the BeO3 complex. (b) CNM analysis for the two reference molecules HOCL2 (DHB donor) and BeH (DHB acceptor) forming the complex.
Figure 6. (a) CNM analysis for the BeO3 complex. (b) CNM analysis for the two reference molecules HOCL2 (DHB donor) and BeH (DHB acceptor) forming the complex.
Molecules 28 00263 g006
Figure 7. (a) Correlation between local mode frequency shifts Δ ω (Hydridic) and Δ ω (Protonic). (b) Correlation between the relative frequency shifts Δ ω / ω (Hydridic) and Δ ω / ω (Protonic). (c) Correlation between Δ ω / ω (Protonic) and k a (DHB). (d) Correlation between Δ ω / ω (Hydridic) and k a (DHB).
Figure 7. (a) Correlation between local mode frequency shifts Δ ω (Hydridic) and Δ ω (Protonic). (b) Correlation between the relative frequency shifts Δ ω / ω (Hydridic) and Δ ω / ω (Protonic). (c) Correlation between Δ ω / ω (Protonic) and k a (DHB). (d) Correlation between Δ ω / ω (Hydridic) and k a (DHB).
Molecules 28 00263 g007
There is no significant correlation between hydridic and protonic frequency shifts (R 2 = 0.6666 and 0.5592, respectively; see Figure 7a,b), although some trends are visible. Overall, larger protonic frequency shifts are in line with larger hydridic frequency shifts, as found for BeO3, BeO4 and on the other end of the spectrum smaller protonic frequency shifts are in line with smaller hydridic frequency shifts, as found for AB4. However, there are also a number of outliers, such as LiO and NaO complexes. Because of the large difference in the frequency shift values for protonic and hydridic parts of DHBs, for the following correlations only the latter were used.
Next, we explored the question if there is a correlation between the frequency shifts of the protonic and hydridic parts of DHB, and the strength of the DHB, as suggested in the literature [43,69,70]. Figure 7c shows the correlation between the relative frequency shifts Δ ω / ω (Protonic) and the local mode force constant k a (DHB), and Figure 7d shows the same situation for Δ ω / ω (Hydridic). We found some general trends for the protonic parts of DHBs shifts (R 2 = 0.7330). Stronger DHBs are characterized by larger red shifts, as found for BeO2, BeO4 and weaker DBH are characterized by small red shifts or even small blue shifts, such as in the case of AB4. Although the same overall trend can be observed for the hydridic frequency shifts scattering of data points, and the outliers LiO and NaO are more pronounced. This shows that the strength of the DHB cannot be directly deduced from the frequency shifts of the protonic or hydridic parts of DHBs upon complex formation.
Figure 8a shows the correlation between protonic frequency shifts and binding energies and Figure 8b the corresponding correlation for the hydridic frequency shifts. Again, in both cases no significant correlation could be found (R 2 = 0.7768 and 0.5933, respectively; see Figure 8a,b), except the general trend that larger binding energies corresponding to stronger DHBs are characterized by larger red shifts, whereas smaller binding energies corresponding to weaker DHBs are characterized by smaller red shifts or even a small blue shift. These findings are consistent with the observed correlation between frequency shifts and local mode DHB force constants. Figure 8c,d show the correlation between frequency shifts and scaled energy density H c / ρ c for the protonic and hydridic parts of DHB, respectively. In both cases, there is no correlation between these two quantities in contrast to H c / ρ c (DHB), which shows at least a weak correlation with k a (DHB) (see Figure 3a).
Figure 8. (a) Correlation between E b i n /n and Δ ω / ω (Protonic). (b) Correlation between E b i n /n and Δ ω / ω (Hydridic). (c) Correlation between H c / ρ c (Protonic) and Δ ω / ω (Protonic). (d) Correlation of H c / ρ c (Hydridic) and Δ ω / ω (Hydridic).
Figure 8. (a) Correlation between E b i n /n and Δ ω / ω (Protonic). (b) Correlation between E b i n /n and Δ ω / ω (Hydridic). (c) Correlation between H c / ρ c (Protonic) and Δ ω / ω (Protonic). (d) Correlation of H c / ρ c (Hydridic) and Δ ω / ω (Hydridic).
Molecules 28 00263 g008

2.2. Some Group Specific Highlights

Group I: This group represents DHBs being built from hydrides, such as LiH, NaH, BeH 2 , BH 4 , and AlH 4 interacting with simple organic molecules—for example, H 2 O , HNO, CH 4 , CFH 3 , and CF 2 H 2 —and also some cations, such as H 3 O + , NH 4 + , CFH 4 + , taken from the work of Bakhmutov [11]. With 52 members, this is the largest group investigated in this work, with DHB lengths stretching from 1.091 Å to 2.533 Å and BSOs from 0.159 to 0.415. Among them, the strongest H−H interactions were observed for the molecular BeO clusters and the weakest interaction was observed for CF sharing three DHBs. It is interesting to note that BC4 also sharing three DHBs has a longer DHB length of 2.533 Å, compared to 2.484 Å found for the DHBs of CF; however, both CF and BC4 DHBs are of comparable strengths (BSO n = 0.159 versus BSO n = 0.168, respectively), another exemption of the generalized Badger rule.
As discussed above, we observed larger changes in the local mode frequencies of the protonic parts of DHB than for the hydridic parts upon DHB complex formation, a trend which also applies to Group I representatives. An interesting example is the BeO3 complex that has a DHB binding energy of −29.61 kcal/mol. The local mode frequency of the protonic H–O bond has a value of 3584 cm 1 in the HOCL2 monomer, while in the DHB complex this frequency has a value of 1163 cm 1 , which corresponds to a large red–shift of −2421 cm 1 . On the other hand, in the similar BeF DHB complex with a DHB binding energy of −4.08 kcal/mol, the red–shift of the protonic bond H–F is with a value of −298 cm 1 almost ten times smaller. This large difference in the red–shifts can be explained by the different electron density transfer from the hydridic Be–H bond into the antibonding σ * orbital of the protonic H–O and H–F bonds. According to our calculations, the population of the H–O σ * orbital in the BeO3 complex has a value of 0.248 e, while the population of the H–F σ * orbital in the BeF has only a value of 0.019 e, which is more than ten times smaller. These results quantify the findings of a previous study, suggesting that the experimentally observed red–shifts in DHBs depend on the amount of electron density transfer into the antibonding σ * orbitals of the protonic bonds [43]. The larger this transfer is the weaker becomes the protonic part of DHB. This is also reflected in the series of charge assisted DHBs [11] BeO1, BeO2, and BeO3. Starting from H 3 O + successive chlorination, i.e., increasing the electrophilicity of the protonic group leads to a considerable decrease of the DHB length (1.237, 1.170, and 1.091 Å, respectively), however the DHB strength remains fairly unchanged (BSO n values of 0.414, 0.414, and 0.397).
Group II: This group includes dimer (AB1), trimer (AB2), and tetramer (AB3) complexes of ammonia–borane H 3 B NH 3 as well as some reference complexes (AB4AB7). Ammonia borane polymers have been suggested as storage medium for hydrogen [73], therefore we assessed these systems how, starting from AB1, the strength of the H−H interaction changes. According to the results of our study the H−H interaction between the monomeric units increases as reflected by the BSO values (0.199, 0.226, and 0.244, for AB1, AB2, and AB3, respectively), whereas the DHB distances remain fairly constant (2.003, 1.916, and 2.069 Å, for AB1, AB2, and AB3 respectively). These results suggest an increasing DHB strength trend in larger ammonia borane clusters and ammonia borane polymers. An interesting alternative is AB7, namely paring ammonia borane with HF leading to the shortest and strongest DHB in this group (BSO = 0.291 and R(DHB) = 1.582 Å). A follow–up study is planned.
Group III: Group 3 includes transition metal clusters made up of iron–phosphine complexes forming either one or two DHBs with ammonia, water, or HF. DHB length range from 1.451 to 2.198 Å wherein forming two DHB interactions (FeN2, FeO2, and FeF2) leads to longer and weaker DHB interactions than those for the complexes with one DHB (FeN1, FeO1, and FeF1) which parallels our observations for Group I members. Unsurprisingly, HF forms the shortest and strongest DHB interactions in this group (FeF1: R(DHB) = 1.451 Å, BSO n = 0.363). As compared to other groups, Group III shows smaller local mode red shifts for both protonic and hydridic parts of DHBs, with the smallest protonic and hydridic part shifts found for FeN1 (−0.017 and −0.006 cm 1 , respectively). Group III members overall show the best correlation between protonic and hydridic frequency shifts in relation to the other groups, likely because the shifts themselves are so small. We also observe a stronger correlation between frequency shifts and BEs than for any of the other groups. However, Group III members are also the most similar compared to the members in other groups, which may cause the better correlation.
The BEs of Group III range from −3.81 to −11.14 kcal/mol for FeN1 and FeF2, respectively, which are among the smallest BEs of the entire set of complexes investigated in this work. Thus, we assert that such transition metal complexes do not tend to form strong DHBs. This difference from other groups can be explained by the presence of d–orbitals, which can withdraw electron density from antibonding σ * orbitals of the hydridic parts of DHBs, leading to complexes with predominately electrostatic character, as reflected by normalized energy density H c / ρ c shown in Figure 3, with the exception of FeF1 and FeF2. This is an interesting aspect to be further explored in the future.
Group IV: This group features four noble gas complexes, whose DHB bond lengths range from 1.322 to 1.866 Å, with BSO n values between 0.246 and 0.324. Although these DHB strengths are on the lower end of the spectrum, they are predominantly covalent in nature, as reflected by the normalized energy density H c / ρ c values shown in Figure 3. For Group IV members, we find exclusively local mode red shifts indicating the weakening of both the protonic and hydridic part of DHB upon complexation. BeAr shows with −0.243 the most significant relative protonic shift in this group, as well as the largest BSO n value (0.324). This can be explained by the linearity of the DHB as well as the presence of fewer shielding electrons in argon. As with the other groups, there is a weak correlation between the frequency shifts and BE.
Group V: Naphthalene, anthracene and phenanthrene Ar1Ar3 have been added to the set of DHB complexes in order to compare the H−H interaction in these aromatic systems, where both hydrogen atoms involved are equally charged, with the DHB properties of the other groups investigated in this work. In previous work [66] we could disprove arguments in favor of DHB based on the existence of a bond path and bond critical point between these two H atoms and the stability of anthracene [11,22,65] as dubious [66]. The sole existence of a bond path and a bond critical point does not necessarily imply the existence of a chemical bond or weak chemical interaction, it may be just an artifact of a closer contact between the two atoms in question [66,74,75,76,77,78,79,80]. In addition, previous attempts to explain the higher stability of phenanthrene via a maximum electron density path between the bay H atoms are misleading in view of the properties of the electron density distribution in the bay region. The 6.8 kcal/mol larger stability of phenanthrene relative to anthracene predominantly (84%) results from its higher resonance energy, which is a direct consequence of the topology of ring annulation [66]. In summary, these observations clearly disqualify the H−H in these systems as DHBs. With the higher level of theory used in this study, we could no longer find a bond path between the bay H atoms for Ar1 and Ar2. The bond path and bond critical point found for Ar3 are an artifact of the closed contact (R(HH) = 1.998 Å). The relative large BSO n values of 0.370, 0.370 and 0.379 (see Figure 2a) are also an artifact of the topology of the bay H atoms. In addition, Figure 2d,f as well as Figure 3 and Figure 4b,c clear identify these complexes as outliers.
Group VI: This is a special group devoted to modeling of the strength of DHB in hydrogen clusters, which have been suggested in the interstellar space [81]. Figure 9a–c present results of our investigations for molecular clusters H4P, H5P, H6P, H7P, and H9P involving diatomic hydrogen H2, and triatomic hydrogen cation H3P. The BSO data as a function of the local mode force constants k a of DHBs for these clusters are shown in Figure 9a, whose values are based on the reference molecules presented in the computational details of this study. According to Figure 9a, the strongest DHB (k a = 5.835 mdyn/Å) is for diatomic hydrogen H2, and the weakest (k a = 0.101 mdyn/Å) is for DHB between the H 3 + and H 2 units of H7P. The strengths of DHB for this Group investigated in our study correlate (R 2 = 0.9218) with the DHB lengths, as shown in Figure 9b. The shortest DHB (R = 0.743 Å) is for H2, and the longest (R = 1.675 Å) for DHB between the H 3 + and H 2 units of H9P. Similarly, according to Figure 9c there is a correlation observed in our study between the strengths of DHB and the covalent characters of DHB expressed by the normalized energy densities H c / ρ c , although this correlation is less perfect (R 2 = 0.8981). The strongest DHB of H2 has the most negative normalized energy density (H c / ρ c = −1.125 Hartree/electron), while the weakest DHB between the H 3 + and H 2 units of H7P has one of the less negative normalized energy density (H c / ρ c = −0.194 Hartree/electron). Generally, we observed in our investigation of this Group, formation of weak DHBs between the H 3 + and H 2 units (light blue and orange colors in H7P and H9P in Figure 9), whose strength is smaller than the strength of DHBs within the H 3 + unit (dark blue, red, and brown colors in H7P and H9P in Figure 9), and the strength of DHB in the H 2 unit of these clusters (violet color in H7P and H9P in Figure 9). The H6P cluster represents formation of weak DHBs (yellow color in H6P in Figure 9) between the H 2 + and H 2 units (dark blue color in H6P in Figure 9), while the H5P cluster shows formation of weak DHBs between the H + and H 2 units (red color in H5P in Figure 9). However the H4P cluster presents formation of weak DHB between the H 3 + unit and the H atom (orange color in H4P in Figure 9). Generally, according to our calculations, the hydrogen clusters are formed between the positively charged unit of the cluster (H + in H5P, H 2 + in H6P, and H 3 + in H7P and H9P) and the H 2 unit. Therefore, we conclude that the results of this study are suggesting possible applications for hydrogen storage devices, where a positively charged central kernel of a molecular complex is surrounded by H 2 molecules, which can easily dissociate.
Figure 9. (a) BSO n(HH) as a function of the local mode force constant k a (HH); power relationship based on H 2 and H 2 + references see text. (b) Correlation of distance R(HH) and k a (HH). (c) Correlation of k a (HH) and normalized energy density H c / ρ c .
Figure 9. (a) BSO n(HH) as a function of the local mode force constant k a (HH); power relationship based on H 2 and H 2 + references see text. (b) Correlation of distance R(HH) and k a (HH). (c) Correlation of k a (HH) and normalized energy density H c / ρ c .
Molecules 28 00263 g009
Table 3. NBO charges (e) for the hydrogen clusters.
Table 3. NBO charges (e) for the hydrogen clusters.
ComplexAtomqBond Δ q a
Group VIH2 0.000H-H0.000
H2P 0.500H-H0.000
H3P 0.333H1-H10.000
H4PH10.161H1-H20.010
H20.171H2-H20.000
H30.496H3-H10.334
H5PH10.281H1-H20.102
H20.180H2-H20.000
H6PH10.219H1-H10.000
H30.141H1-H30.078
H7PH10.209H1-H20.075
H20.284H2-H20.000
H30.056H2-H30.228
H3-H30.000
H9PH10.262H1-H30.227
H30.035H3-H30.000

3. Methods and Computational Details

As a major assessment tool of H–H bond strength, we utilized LMA in this work. A comprehensive review of the underlying theory and recent applications can be found in Refs. [72,82]. Therefore, in the following, only some essentials are summarized.
Normal vibrational modes are generally delocalized as a result of kinematic and electronic coupling [83,84,85]. A particular vibrational stretching mode between two atoms of interest can couple with other vibrational modes such as bending or torsion, which inhibits the direct correlation between normal stretching frequency or associated normal mode force constant and bond strength, as well as the comparison between normal stretching vibrational modes in related molecules. As a consequence, the normal stretching force constant cannot be used as a direct bond strength measure. One needs to derive local counterparts that are free from any mode–mode coupling. Konkoli, Cremer and co–workers solved this problem by solving the mass–decoupled analogue of Wilson’s equation of vibrational spectroscopy [53,54,55,56,57], leading to local vibrational modes, associated local mode frequencies and local mode force constants.
There is a 1:1 relationship between the normal vibrational modes and a complete non–redundant set of local vibrational modes via an adiabatic connection scheme (ACS) [86], which forms the foundation for a new method of analyzing vibrational spectra, the characterization of normal mode (CNM) procedure [55,56,57,66,72,86,87,88,89,90,91]. CNM decomposes each normal mode into local mode contributions leading to a wealth of information about structure and bonding. Recently, CNM has been successfully applied in the investigation of pK a probes [92] and to assess the qualification of vibrational Stark effect probes [93].
Zou and Cremer showed that the local stretching force constant k a reflects the curvature of the PES in the direction of the bond stretching [94]. This important result qualifies the local stretching force constants k a as unique quantitative measures of the intrinsic strength of chemical bonds and/or weak chemical interactions based on vibrational spectroscopy. Therefore, LMA has been extensively applied for the study of chemical bonding and non–bonded interactions, as documented in Refs. [72,82] and references cited therein, including the characterization of HB ranging from conventional HBs [67], unconventional HBs [95,96], HBs in water and ice [87,97,98,99,100,101,102], and HBs in biomolecules [103,104,105,106] to HBs in catalysis [107,108]—a list which has been extended by DHB in this work. The local vibrational modes theory allows us to analyze a variety of internal coordinates in addition to bond lengths, such as angles between bonds, dihedral bond angles, puckering coordinates, and others. However, this study is focused on the bond lengths leading to the local mode force constants directly related to the strength of DHB.
The widely accepted definition of a DHB includes the following characteristics: DHBs H δ + H δ are formed from two H–atoms with opposite partial charges and a distance R(DHB) between the two H–atoms that is smaller than the sum of their van der Waals radii (2.4 Å); the angle X–H δ + H δ of the protonic fragment is closer to linear, whereas the hydridic fragment H δ + H δ –Y is more bent [11]. For a larger separation, it has been suggested to apply a criterion based on the topological electron density: The density ρ c at the bond critical point r c between the two H–atoms should be small, and the Laplacian 2 ( ρ c ) should be small and positive. Whereas bond energies are assumed to be similar for both HBs and DHBs, DHBs should have a smaller covalent character, i.e., a more pronounced electrostatic nature [11]. As discussed above, in addition to R(DHB) [1,11], the red/blue shifts of the normal mode frequency associated with the DHB [9,43], as well as components of the binding energies (BE) determined with different energy decomposition methods [40,41,42] have been applied as a measure of the DHB bond strength.
However a caveat is necessary: both the bond dissociation energies (BDE) and BE [109,110,111] as well as associated energy decomposition schemes (EDA) [112,113,114], bond distances [115,116] and normal mode frequency shifts, are not suited to measure direct bond strength. The BDE is a reaction parameter that includes all changes taking place during the dissociation process. Accordingly, it includes any (de)stabilization effects of the fragments to be formed. It reflects the energy needed for bond breaking, but also contains energy contributions due to geometry relaxation and electron density reorganization in the dissociation fragments [67]. Therefore, the BDE is not a suitable measure of the intrinsic strength of a chemical bond and its use may lead to misjudgments, as documented in the literature [117,118,119,120,121,122]. In addition, EDA schemes are not free from arbitrariness [123]. Additionally, the bond length is not always a qualified bond strength descriptor. Numerous cases have been reported, illustrating that a shorter bond is not always a stronger bond [54,55,57,124,125].
The use of normal mode stretching frequencies and/or frequency shifts as convenient bond strength measures (both are accessible experimentally and computationally) is problematic, because of the coupling of the stretching mode with other vibrational modes, which hampers a direct comparison between bonds in strongly different systems, as pointed out above. The topological analysis of the electron density can be useful to uncover possible attractive contacts between two atoms via the existence of a maximum electron density path (i.e., bond path), with a bond critical point connecting the two nuclei under consideration [126]. However, the sole existence of a bond path and a bond critical point does not necessarily imply the existence of a chemical bond; in particular, the QTAIM description of weak chemical interactions may be problematic [46,66,74,75,76,78,79]. Using the Laplacian 2 ( ρ c ) as a complementary measure is problematic too [127,128,129], because it does not reflect the complex interplay between kinetic and potential energy, which accompanies the bond forming process [130,131].
Therefore, in this work we utilized LMA to derive a quantitative DHB bond strength measure, in particular local mode force constants and related bond strength orders. It is convenient to base the comparison of the bond strength for a set of molecules on a chemically more prevalent bond strength order (BSO n) rather than on a comparison of local force constant values. Both are connected via a power relationship according to the generalized Badger rule derived by Cremer and co–workers [68]: BSO n = A ( k a ) B . Constants A and B can be determined from two reference compounds with known BSO n values and the requirement that, for a zero force constant, the BSO n must be zero. For the protonic parts of DHBs, the hydridic parts of DHBs, and the H–H interaction in DHBs investigated in this study, we use as references the FH bond in the FH molecule with BSO n = 1 and the FH bond in the [F⋯H⋯F] anion with BSO n = 0.5 [67,99,103,106]. For the ω B97X-D [132,133] /aug-cc-pVTZ [134,135] model chemistry, applied in this study, this led to k a (FH) = 9.719 mdyn/Å, k a (F⋯H) = 0.803 mdyn/Å and A = 0.531 and B = 0.278 . For the DHBs in hydrogen clusters of Group VI investigated in this work, we used as references H 2 (BSO n = 1, k a = 5.840 mdyn/Å) and H 2 + (BSO n = 0.5, k a = 1.263 mdyn/Å), leading to A = 0.450 and B = 0.453 .
LMA was complemented with QTAIM [61,62,78,126], where in this work we applied the Cremer–Kraka criterion [127,128,129] of covalent bonding to assess the covalent/electrostatic character of DHBs. The Cremer–Kraka criterion is composed of two conditions: (i) the necessary condition—existence of a bond path and bond critical point r c = c between the two atoms under consideration; (ii) the sufficient condition—the energy density H ( r c ) = H c is smaller than zero. H ( r ) is defined as H(r) = G(r) + V(r), where G ( r ) is the kinetic energy density and V ( r ) is the potential energy density. A negative V ( r ) corresponds to a stabilizing accumulation of density, whereas the positive G ( r ) corresponds to the depletion of electron density [128]. As a result, the sign of H c indicates which term is dominant [129]. If H c < 0 , the interaction is considered covalent in nature, whereas H c > 0 is indicative of electrostatic interactions. In addition to the QTAIM analysis, we calculated Natural Bond Orbital (NBO) [58,59] atomic charges and the corresponding charge transfer between the two hydrogen atoms involved in DHB, and we introduced local mode frequency shifts [67] as a more reliable parameter to distinguish different DHB situation than normal mode frequency shifts, which are contaminated by mode–mode coupling.
The geometries of all group members and references compounds (see Figure 1) were optimized using the ω B97X–D density functional [132,133] combined with Dunning’s aug–cc–pVTZ basis set [134,135]. The ω B97X–D functional was chosen as it has been proven to reliably describe weak (long–range) intermolecular interactions covering the diverse range of molecules considered in particular in combination with the augmented basis set [133,136]. Using the same model chemistry as applied in our comprehensive study on HBs [137] also allows the direct comparison of HB and DHB features. Furthermore, this model chemistry has been recently assessed and characterized as useful for the description of DHB [39]. Geometry optimizations were followed by a normal mode analysis and LMA. Geometry optimizations, normal mode analyses and NBO analyses were performed with the Gaussian09 program package [138]. For the subsequent LMA the program LModeA [139] was utilized. The QTAIM analysis was performed with the AIMAll package [140].

4. Conclusions

In this computational study, we investigated the strength and nature of dihydrogen bonding (DHB) in a variety of molecular complexes and clusters, categorized into six different groups. As an assessment tool, we used the Local Mode Analysis (LMA) based on DFT calculations involving dispersion corrections and a flexible basis set. LMA, which was used to determine the precise strength of individual DHBs in these complexes and clusters, was supported by the topological analysis of the electron density, quantifying the covalent versus electrostatic character of these DHBs. Along with the analysis of DHB, we also investigated the strength and the nature of the protonic bonds X–H δ + , as well as the hydridic bonds H δ –Y forming the DHBs (i.e., X–H δ + H δ –Y) for Group I–V members. According to our results, the strength of DHB is modulated by the strength of the protonic bonds X–H δ + rather than by the strength of the hydridic bonds H δ –Y, which was also confirmed for Group VI members involving hydrogen clusters. In accordance with our calculations, the hydrogen clusters were formed by weak DHB between the central positively charged H 3 + kernel of the cluster and the surrounding H 2 molecules, which can easily dissociate, suggesting a possible application for hydrogen storage.
We hope that our study inspires the community to add the Local Mode Analysis to their repertoire for the investigation of chemical bonds and intermolecular interactions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28010263/s1. Table S1: NBO hydrogen charges for reference compounds; Table S2: With bond lengths R in Å, local mode force constants k a in mdyn/Å, local mode vibrational frequencies ω a in cm 1 , binding energies E in kcal/mol, electron densities at the bond critical points ρ c in e/Å 3 and energy densities at critical points H c in Hartree/Å 3 of protonic and hydridic HBs and DHBs for all complexes and references compounds.

Author Contributions

Conceptualization, E.K. and M.F.; methodology, M.F. and E.K.; validation, E.K. and M.F.; formal analysis, M.F., N.B. and E.K.; investigation, M.F., N.B. and E.K.; resources, E.K.; supervision, E.K.; data curation, M.F. and M.M.; writing—original draft preparation, M.F. and E.K.; writing—review and editing, E.K., M.F. and M.M.; funding acquisition, E.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science Foundation NSF, grant CHE 2102461.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data supporting the results of this work are presented in tables and figure of the manuscript and in the Supplementary Materials.

Acknowledgments

This work was supported by the National Science Foundation, Grant CHE 2102461. We thank the Center for Research Computation at SMU for providing generous high-performance computational resources.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Not applicable.

References

  1. Karas, L.J.; Wu, C.H.; Das, R.; Wu, J.I.C. Hydrogen bond design principles. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2020, 10, e1477. [Google Scholar] [CrossRef] [PubMed]
  2. Grabowski, S.J. Understanding Hydrogen Bonds; Theoretical and Computational Chemistry Series; The Royal Society of Chemistry: London, UK, 2021. [Google Scholar]
  3. Scheiner, S. Understanding noncovalent bonds and their controlling forces. Int. J. Chem. Phys. 2020, 153, 140901. [Google Scholar] [CrossRef] [PubMed]
  4. Arunan, E. One Hundred Years After The Latimer and Rodebush Paper, Hydrogen Bonding Remains an Elephant! Indian J. Sci. 2020, 100, 249–255. [Google Scholar] [CrossRef] [Green Version]
  5. Gibb, B.C. The Centenary (maybe) of The Hydrogen Bond. Nat. Chem. 2020, 12, 665–667. [Google Scholar] [CrossRef] [PubMed]
  6. Mallamace, D.; Fazio, E.; Mallamace, F.; Corsaro, C. The role of hydrogen bonding in the folding/unfolding process of hydrated lysozyme: A review of recent NMR and FTIR results. Int. J. Mol. Sci. 2018, 19, 3825. [Google Scholar] [CrossRef] [Green Version]
  7. Bulusu, G.; Desiraju, G.R. Strong and Weak Hydrogen Bonds in Protein–Ligand Recognition. J. Indian Inst. Sci. 2019, 100, 31–41. [Google Scholar] [CrossRef]
  8. Mobika, J.; Rajkumar, M.; Sibi, S.L.; Priya, V.N. Investigation on hydrogen bonds and conformational changes in protein/polysaccharide/ceramic based tri-component system. Spectrochim. Acta A Mol. Biomol. Spectrosc. 2021, 244, 118836. [Google Scholar] [CrossRef]
  9. Belkova, N.V.; Epstein, L.M.; Filippov, O.A.; Shubina, E.S. Hydrogen and dihydrogen bonds in the reactions of metal hydrides. Chem. Rev. 2016, 116, 8545–8587. [Google Scholar] [CrossRef]
  10. Gadwal, I. A Brief Overview on Preparation of Self-Healing Polymers and Coatings via Hydrogen Bonding Interactions. Macromol 2021, 1, 18–36. [Google Scholar] [CrossRef]
  11. Bakhmutov, V.I. Dihydrogen Bonds; John Wiley & Sons: Hobroken, NJ, USA, 2008. [Google Scholar]
  12. Crabtree, R.H. Dihydrogen complexation. Chem. Rev. 2016, 116, 8750–8769. [Google Scholar] [CrossRef]
  13. Kubas, G.J. Dihydrogen complexes as prototypes for the coordination chemistry of saturated molecules. Proc. Natl. Acad. Sci. USA 2007, 104, 6901–6907. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Grabowski, S.J. Molecular Hydrogen as a Lewis Base in Hydrogen Bonds and Other Interactions. Molecules 2021, 25, 3294. [Google Scholar] [CrossRef] [PubMed]
  15. Richardson, T.; de Gala, S.; Crabtree, R.H.; Siegbahn, P.E.M. Unconventional hydrogen bonds: Intermolecular B-H⋯H-N interactions. J. Am. Chem. Soc. 1995, 117, 12875–12876. [Google Scholar] [CrossRef]
  16. Custelcean, R.; Jackson, J.E. Dihydrogen bonding: Structures, energetics, and dynamics. Chem. Rev. 2001, 101, 1963–1980. [Google Scholar] [CrossRef]
  17. Alkorta, I.; Elguero, J. Non-conventional hydrogen bonds. Chem. Soc. Rev. 1998, 27, 163–170. [Google Scholar] [CrossRef]
  18. Epstein, L.M.; Shubina, E.S. New types of hydrogen bonding in organometallic chemistry. Coord. Chem. Rev. 2002, 231, 165–181. [Google Scholar] [CrossRef]
  19. Solimannejad, M.; Amlashi, L.M.; Alkorta, I.; Elguero, J. XeH2 as a proton-accepting molecule for dihydrogen bonded systems: A theoretical study. Chem. Phys. Lett. 2006, 422, 226–229. [Google Scholar] [CrossRef]
  20. Alkorta, I.; Elguero, J.; Mo, O.; Yanez, M.; Del Bene, J.E. Ab Initio study of the structural, energetic, bonding, and IR spectroscopic properties of complexes with dihydrogen bonds. J. Phys. Chem. A 2002, 106, 9325–9330. [Google Scholar] [CrossRef]
  21. Alkorta, I.; Zborowski, K.; Elguero, J.; Solimannejad, M. Theoretical study of dihydrogen bonds between (XH)2, X = Li, Na, BeH, and MgH, and weak hydrogen bond donors (HCN, HNC, and HCCH). J. Phys. Chem. A 2006, 110, 10279–10286. [Google Scholar] [CrossRef]
  22. Grabowski, S.J.; Sokalski, W.A.; Leszczynski, J. Wide spectrum of H⋯H interactions: Van der Waals contacts, dihydrogen bonds and covalency. Chem. Phys. 2007, 337, 68–76. [Google Scholar] [CrossRef]
  23. Castilla-Martinez, C.A.; Moury, R.; Demirci, U.B. Amidoboranes and hydrazinidoboranes: State of the art, potential for hydrogen storage, and other prospects. Int. J. Hydrog. Energy 2020, 45, 30731–30755. [Google Scholar] [CrossRef]
  24. Magos-Palasyuk, E.; Litwiniuk, A.; Palasyuk, T. Experimental and theoretical evidence of dihydrogen bonds in lithium amidoborane. Sci. Rep. 2020, 10, 17431. [Google Scholar] [CrossRef] [PubMed]
  25. Zhou, W.; Jin, S.; Dai, W.; Lyon, J.T.; Lu, C. Theoretical study on the structural evolution and hydrogen storage in NbHn (n = 2–15) clusters. Int. J. Hydrog. Energy 2021, 46, 17246–17252. [Google Scholar] [CrossRef]
  26. Zachariasen, W.H.; Mooney, R.C.L. The structure of the hypophosphite group as determined from the crystal lattice of ammonium hypophosphite. J. Chem. Phys. 1934, 2, 34–38. [Google Scholar] [CrossRef]
  27. Burg, A.B. Enhancement of P-H bonding in a phosphine monoborane. Inorg. Chem. 1964, 3, 1325–1327. [Google Scholar]
  28. Brown, M.P.; Heseltine, R.W. Co-ordinated BH3 as a proton acceptor group in hydrogen bonding. Chem. Comm. 1968, 1551–1552. [Google Scholar] [CrossRef]
  29. Brown, M.P.; Heseltine, R.W.; Sutcliffe, L.H. Studies of aminoboranes. Part II. Pyrolysis of monoalkylamine-boranes: Some aminoborane trimers (1,3,5-trialkylcyclotriborazanes) and dimers. J. Chem. Soc. 1968, 612–616. [Google Scholar] [CrossRef]
  30. Brown, M.P.; Heseltine, R.W.; Smith, P.A.; Walker, P.J. An infrared study of co-ordinated BH3 and BH2 groups as proton acceptors in hydrogen bonding. J. Chem. Soc. 1970, 410–414. [Google Scholar] [CrossRef]
  31. Lee, J.C., Jr.; Rheingold, A.L.; Muller, B.; Pregosin, P.S.; Crabtree, R.H. Complexation of an amide to iridium via an iminol tautomer and evidence Ir-H⋯H-O hydrogen bond. J. Chem. Soc. 1994, 1021–1022. [Google Scholar] [CrossRef]
  32. Lough, A.J.; Park, S.; Ramachandran, R.; Morris, R.H. Switching on and off a new intramolecular hydrogen-hydrogen interaction and the heterolytic splitting of dihydrogen. Crystal and molecular structure. J. Am. Chem. Soc. 1994, 116, 8356–8357. [Google Scholar] [CrossRef]
  33. Peris, E.; Lee, J.C.; Rambo, J.R.; Eisenstein, O.; Crabtree, R.H. Factors affecting the strength of N-H⋯H-Ir hydrogen bonds. J. Am. Chem. Soc. 1995, 117, 3485–3491. [Google Scholar] [CrossRef]
  34. Clot, E.; Eisenstein, O.; Crabtree, R.H. How hydrogen bonding affects ligand binding and fluxionality in transition metal complexes: A DFT study on interligand hydrogen bonds involving HF and H2O. New J. Chem. 2001, 25, 66–72. [Google Scholar] [CrossRef]
  35. Groom, C.R.; Bruno, I.J.; Lightfoot, M.P.; Ward, S.C. The Cambridge Structural Database. Acta Cryst. B 2016, 72, 171–179. [Google Scholar] [CrossRef] [PubMed]
  36. Crabtree, R.H.; Siegbahn, P.E.M.; Eisenstein, O.; Rheingold, A.L.; Koetzle, T.F. A new intermolecular interaction: Unconventional hydrogen bonds with element hydride bonds as proton acceptor. Acc. Chem. Res. 1996, 29, 348–354. [Google Scholar] [CrossRef] [PubMed]
  37. Jayaraman, A. 100th Anniversary of Macromolecular Science Viewpoint: Modeling and Simulation of Macromolecules with Hydrogen Bonds: Challenges, Successes, and Opportunities. ACS Macro Lett. 2020, 2, 656–665. [Google Scholar] [CrossRef] [PubMed]
  38. Bauer, C.A. How to Model Inter- and Intramolecular Hydrogen Bond Strengths with Quantum Chemistry. J. Chem. Inf. Model. 2019, 59, 3735–3743. [Google Scholar] [CrossRef]
  39. Fanfrlík, J.; Pecina, A.; Řezáč, J.; Lepšík, M.; Sárosi, M.B.; Hnyk, D.; Hobza, P. Benchmark data sets of boron cluster dihydrogen bonding for the validation of approximate computational methods. ChemPhysChem 2020, 21, 2599–2604. [Google Scholar] [CrossRef]
  40. Van der Lubbe, S.C.C.; Guerra, C.F. The Nature of Hydrogen Bonds: A Delineation of the Role of Different Energy Components on Hydrogen Bond Strengths and Lengths. Chem. Asian J. 2019, 14, 2760–2769. [Google Scholar]
  41. Karimi, S.; Sanchooli, M.; Shoja-Hormozzah, F. Estimation of resonance assisted hydrogen bond (RAHB) energies using properties of ring critical points in some dihydrogen-bonded complexes. J. Mol. Struct. 2021, 1242, 130710. [Google Scholar] [CrossRef]
  42. Deshmukh, M.M.; Gadre, S.R. Molecular Tailoring Approach for the Estimation of Intramolecular Hydrogen Bond Energy. Molecules 2021, 26, 2928. [Google Scholar] [CrossRef]
  43. Uchida, M.; Shimizu, T.; Shibutani, R.; Matsumoto, Y.; Ishikawa, H. A comprehensive infrared spectroscopic and theoretical study on phenol- ethyldimethylsilane dihydrogen-bonded clusters in the S0 and S1 states. J. Chem. Phys. 2020, 153, 104305. [Google Scholar] [CrossRef] [PubMed]
  44. Matta, C.F.; Huang, L.; Massa, L. Characterization of a trihydrogen bond on the basis of the topology of the electron density. J. Phys. Chem. A 2011, 115, 12451–12458. [Google Scholar] [CrossRef] [PubMed]
  45. Nakanishi, W.; Hayashi, S.; Narahara, K. Atoms–in–Molecules Dual Parameter Analysis of Weak to Strong Interactions: Behaviors of Electronic Energy Densities versus Laplacian of Electron Densities at Bond Critical Points. J. Phys. Chem. A 2008, 112, 13593–13599. [Google Scholar] [CrossRef] [PubMed]
  46. Shahbazian, S. Why Bond Critical Points Are Not ‘Bond’ Critical Points. Chem. Eur. J. 2018, 24, 5401–5405. [Google Scholar] [CrossRef] [PubMed]
  47. Mo, Y. Can QTAIM Topological Parameters Be a Measure of Hydrogen Bonding Strength. J. Phys. Chem. A 2012, 116, 5240–5246. [Google Scholar] [CrossRef] [PubMed]
  48. Fuster, F.; Grabowski, S.J. Intramolecular Hydrogen Bonds: The QTAIM and ELF Characteristics. J. Phys. Chem. A 2011, 115, 10078–10086. [Google Scholar] [CrossRef]
  49. Alikhani, M.E.; Fuster, F.; Silvi, B. What Can Tell the Topological Analysis of ELF on Hydrogen Bonding? Struct. Chem. 2005, 16, 204–210. [Google Scholar] [CrossRef]
  50. Shahbazian, S. A Molecular Electrostatic Potential Analysis of Hydrogen, Halogen and Dihydrogen Bonds. J. Phys. Chem. A 2014, 118, 1697–1705. [Google Scholar]
  51. Murray, J.S.; Politzer, P. Hydrogen Bonding: A Coulombic σ? Hole Interaction. J. Indian Inst. Sci. 2020, 100, 21–30. [Google Scholar] [CrossRef]
  52. Murray, J.S.; Politzer, P. The Electrostatic Potential: An Overview. WIREs Comput. Mol. Sci. 2011, 1, 153–163. [Google Scholar] [CrossRef]
  53. Konkoli, Z.; Cremer, D. A New Way of Analyzing Vibrational Spectra. I. Derivation of Adiabatic Internal Modes. Int. J. Quantum Chem. 1998, 67, 1–9. [Google Scholar] [CrossRef]
  54. Konkoli, Z.; Larsson, J.A.; Cremer, D. A New Way of Analyzing Vibrational Spectra. II. Comparison of Internal Mode Frequencies. Int. J. Quantum Chem. 1998, 67, 11–27. [Google Scholar] [CrossRef]
  55. Konkoli, Z.; Cremer, D. A New Way of Analyzing Vibrational Spectra. III. Characterization of Normal Vibrational Modes in terms of Internal Vibrational Modes. Int. J. Quantum Chem. 1998, 67, 29–40. [Google Scholar] [CrossRef]
  56. Konkoli, Z.; Larsson, J.A.; Cremer, D. A New Way of Analyzing Vibrational Spectra. IV. Application and Testing of Adiabatic Modes within the Concept of the Characterization of Normal Modes. Int. J. Quantum Chem. 1998, 67, 41–55. [Google Scholar] [CrossRef]
  57. Cremer, D.; Larsson, J.A.; Kraka, E. New Developments in the Analysis of Vibrational Spectra on the Use of Adiabatic Internal Vibrational Modes. In Theoretical and Computational Chemistry; Parkanyi, C., Ed.; Elsevier: Amsterdam, The Netherlands, 1998; pp. 259–327. [Google Scholar]
  58. Reed, A.; Curtiss, L.; Weinhold, F. Intermolecular Interactions from A Natural Bond Orbital, Donor-Acceptor Viewpoint. Chem. Rev. 1988, 88, 899–926. [Google Scholar] [CrossRef]
  59. Weinhold, F.; Landis, C.R. Valency and Bonding: A Natural Bond Orbital Donor-Acceptor Perspective; Cambridge University Press: Cambridge, UK, 2005. [Google Scholar]
  60. Glendening, E.D.; Badenhoop, J.K.; Reed, A.E.; Carpenter, J.E.; Bohmann, J.A.; Morales, C.M.; Landis, C.R.; Weinhold, F. NBO6; Theoretical Chemistry Institute, University of Wisconsin: Madison, WI, USA, 2013. [Google Scholar]
  61. Bader, R. Atoms in Molecules: A Quantum Theory; Clarendron Press: Oxford, UK, 1995. [Google Scholar]
  62. Bader, R.F.W. The Quantum Mechanical Basis of Conceptual Chemistry. Monatsh. Chem. 2005, 136, 819–854. [Google Scholar] [CrossRef]
  63. Popelier, P.L.A. On Quantum Chemical Topology. In Applications of Topological Methods in Molecular Chemistry; Springer International Publishing: Berlin/Heidelberg, Germany, 2016; pp. 23–52. [Google Scholar]
  64. Matta, C.F.; Boyd, R.J. An Introduction to the Quantum Theory of Atoms in Molecules. In The Quantum Theory of Atoms in Molecules; Wiley-VCH Verlag GmbH & Co., KGaA: Weinheim, Germany, 2007; pp. 1–34. [Google Scholar]
  65. Matta, C.F.; Hernndez-Trujillo, J.; Tang, T.H.; Bader, R.F.W. Hydrogen hydrogen bonding: A stabilizing interaction in molecules and crystals. Chem. Europ. J. 2003, 9, 1940–1951. [Google Scholar] [CrossRef] [PubMed]
  66. Kalescky, R.; Kraka, E.; Cremer, D. Description of Aromaticity with the Help of Vibrational Spectroscopy: Anthracene and Phenanthrene. J. Phys. Chem. A 2013, 118, 223–237. [Google Scholar] [CrossRef]
  67. Freindorf, M.; Kraka, E.; Cremer, D. A Comprehensive Analysis of Hydrogen Bond Interactions Based on Local Vibrational Modes. Int. J. Quantum Chem. 2012, 112, 3174–3187. [Google Scholar] [CrossRef]
  68. Kalescky, R.; Kraka, E.; Cremer, D. Identification of the Strongest Bonds in Chemistry. J. Phys. Chem. A 2013, 117, 8981–8995. [Google Scholar] [CrossRef]
  69. Grabowski, S. Halogen bond and its counterparts: Bent’s rule explains the formation of nonbonding interactions. J. Phys. Chem. A 2011, 115, 12340. [Google Scholar] [CrossRef] [PubMed]
  70. Zhang, X.; Zeng, Y.; Li, X.; Meng, L.; Zheng, S. Comparison in the complexes of oxygen-containing σ-electron donor with hydrogen halide and dihalogen molecules. J. Mol. Struct. THEOCHEM 2010, 950, 27–35. [Google Scholar] [CrossRef]
  71. De Oliveira, B.G.; Ramos, M.N. Dihydrogen Bonds and Blue–Shifting Hydrogen Bonds. AH⋯CF3 and TH2⋯HCF3 Model Systems with A = Li or Na and T = Be or Mg. Int. J. Quant. Chem. 2010, 110, 307–316. [Google Scholar] [CrossRef]
  72. Kraka, E.; Zou, W.; Tao, Y. Decoding Chemical Information from Vibrational Spectroscopy Data: Local Vibrational Mode Theory. WIREs: Comput. Mol. Sci. 2020, 10, 1480. [Google Scholar] [CrossRef]
  73. Li, S.F.; Tang, Z.W.; Tan, Y.B.; Yu, X.B. Polyacrylamide Blending with Ammonia Borane: A Polymer Supported Hydrogen Storage Composite. J. Phys. Chem. C 2012, 116, 1544–1549. [Google Scholar] [CrossRef]
  74. Jabłoński, M. On the Uselessness of Bond Paths Linking Distant Atoms and on the Violation of the Concept of Privileged Exchange Channels. ChemistryOpen 2019, 8, 497–507. [Google Scholar] [CrossRef] [Green Version]
  75. Jabłoński, M. Bond paths between distant atoms do not necessarily indicate dominant interactions. J. Comput. Chem. 2018, 39, 2183–2195. [Google Scholar] [CrossRef]
  76. Wick, C.R.; Clark, T. On bond-critical points in QTAIM and weak interactions. J. Mol. Model. 2018, 24, 142. [Google Scholar] [CrossRef]
  77. Jabłoński, M. A Critical Overview of Current Theoretical Methods of Estimating the Energy of Intramolecular Interactions. Molecules 2020, 25, 5512. [Google Scholar] [CrossRef]
  78. Bader, R.F.W. Bond Paths Are Not Chemical Bonds. J. Phys. Chem. A 2009, 113, 10391–10396. [Google Scholar] [CrossRef] [Green Version]
  79. Stöhr, M.; Voorhis, T.V.; Tkatchenko, A. Theory and practice of modeling van der Waals interactions in electronic–structure calculations. Chem. Soc. Rev. 2019, 48, 4118–4154. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Hazrah, A.S.; Nanayakkara, S.; Seifert, N.A.; Kraka, E.; Jäger, W. Structural study of 1- and 2-naphthol: New insights into the non-covalent H–H interaction in cis-1-naphthol. Phys. Chem. Chem. Phys. 2022, 24, 3722–3732. [Google Scholar] [CrossRef] [PubMed]
  81. Barbatti, M.; Nascimento, M.A.C. Does the H5+ hydrogen cluster exist in dense interstellar clouds? Int. J. Quantum Chem. 2012, 112, 3169–3173. [Google Scholar] [CrossRef]
  82. Kraka, E.; Quintano, M.; Force, H.W.L.; Antonio, J.J.; Freindorf, M. The Local Vibrational Mode Theory and Its Place in the Vibrational Spectroscopy Arena. J. Phys. Chem. A 2022, 126, 8781–8798. [Google Scholar] [CrossRef] [PubMed]
  83. Wilson, E.B.; Decius, J.C.; Cross, P.C.M. Molecular Vibrations. The Theory of Infrared and Raman Vibrational Spectra; McGraw-Hill: New York, NY, USA, 1955; pp. 59–136. [Google Scholar]
  84. Califano, S. Vibrational States; Wiley: London, UK, 1976. [Google Scholar]
  85. Kelley, J.D.; Leventhal, J.J. Problems in Classical and Quantum Mechanics: Normal Modes and Coordinates; Springer: Berlin/Heidelberg, Germany, 2017; pp. 95–117. [Google Scholar]
  86. Zou, W.; Kalescky, R.; Kraka, E.; Cremer, D. Relating Normal Vibrational Modes to Local Vibrational Modes with the Help of an Adiabatic Connection Scheme. J. Chem. Phys. 2012, 137, 084114. [Google Scholar] [CrossRef] [Green Version]
  87. Kalescky, R.; Zou, W.; Kraka, E.; Cremer, D. Local Vibrational Modes of the Water Dimer—Comparison of Theory and Experiment. Chem. Phys. Lett. 2012, 554, 243–247. [Google Scholar] [CrossRef]
  88. Zou, W.; Kalescky, R.; Kraka, E.; Cremer, D. Relating Normal Vibrational Modes to Local Vibrational Modes: Benzene and Naphthalene. J. Mol. Model. 2012, 19, 2865–2877. [Google Scholar] [CrossRef]
  89. Freindorf, M.; Kraka, E. Critical Assessment of the FeC and CO Bond strength in Carboxymyoglobin—A QM/MM Local Vibrational Mode Study. J. Mol. Model. 2020, 26, 281. [Google Scholar] [CrossRef]
  90. Yannacone, S.; Freindorf, M.; Tao, Y.; Zou, W.; Kraka, E. Local Vibrational Mode Analysis of π-Hole Interactions between Aryl Donors and Small Molecule Acceptors. Crystals 2020, 10, 556. [Google Scholar] [CrossRef]
  91. Yannacone, S.; Sayala, K.D.; Freindorf, M.; Tsarevsky, N.V.; Kraka, E. Vibrational Analysis of Benziodoxoles and Benziodazolotetrazoles. PhysChem 2021, 1, 45–68. [Google Scholar] [CrossRef]
  92. Quintano, M.; Kraka, E. Theoretical insights into the linear relationship between pKa values and vibrational frequencies. Chem. Phys. Lett. 2022, 803, 139746. [Google Scholar] [CrossRef]
  93. Verma, N.; Tao, Y.; Zou, W.; Chen, X.; Chen, X.; Freindorf, M.; Kraka, E. A Critical Evaluation of Vibrational Stark Effect (VSE) Probes with the Local Vibrational Mode Theory. Sensors 2020, 20, 2358. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Zou, W.; Cremer, D. C2 in a Box: Determining its Intrinsic Bond Strength for the X1Σ+g Ground State. Chem. Eur. J. 2016, 22, 4087–4097. [Google Scholar] [CrossRef] [PubMed]
  95. Cremer, D. Stereochemistry of the Ozonolysis of Alkenes: Ozonide-versus Carbonyl Oxide-Control. Angew. Chem. Int. Ed. 1981, 20, 888–889. [Google Scholar] [CrossRef]
  96. Zou, W.; Zhang, X.; Dai, H.; Yan, H.; Cremer, D.; Kraka, E. Description of an Unusual Hydrogen Bond between Carborane and a Phenyl Group. J. Organometal. Chem. 2018, 856, 114–127. [Google Scholar] [CrossRef]
  97. Cremer, D. From Configuration Interaction to Coupled Cluster Theory: The Quadratic Configuration Interaction Approach. WIREs Comput. Mol. Sci. 2013, 3, 482–503. [Google Scholar] [CrossRef]
  98. Kalescky, R.; Zou, W.; Kraka, E.; Cremer, D. Vibrational Properties of the Isotopomers of the Water Dimer Derived from Experiment and Computations. Aust. J. Chem. 2014, 67, 426. [Google Scholar] [CrossRef] [Green Version]
  99. Tao, Y.; Zou, W.; Jia, J.; Li, W.; Cremer, D. Different Ways of Hydrogen Bonding in Water - Why Does Warm Water Freeze Faster than Cold Water? J. Chem. Theory Comput. 2017, 13, 55–76. [Google Scholar] [CrossRef]
  100. Tao, Y.; Zou, W.; Kraka, E. Strengthening of Hydrogen Bonding With the Push-Pull Effect. Chem. Phys. Lett. 2017, 685, 251–258. [Google Scholar] [CrossRef]
  101. Delgado, A.A.A.; Sethio, D.; Kraka, E. Assessing the Intrinsic Strengths of Ion-Solvent and Solvent-Solvent Interactions for Hydrated Mg2+ Clusters. Inorganics 2021, 9, 31. [Google Scholar] [CrossRef]
  102. Nanayakkara, S.; Tao, Y.; Kraka, E. Capturing Individual Hydrogen Bond Strengths in Ices via Periodic Local Vibrational Mode Theory: Beyond the Lattice Energy Picture. J. Chem. Theory Comput. 2022, 18, 562–579. [Google Scholar] [CrossRef] [PubMed]
  103. Lyu, S.; Beiranvand, N.; Freindorf, M.; Kraka, E. Interplay of Ring Puckering and Hydrogen Bonding in Deoxyribonucleosides. J. Phys. Chem. A 2019, 123, 7087–7103. [Google Scholar] [CrossRef]
  104. Yannacone, S.; Sethio, D.; Kraka, E. Quantitative Assessment of Intramolecular Hydrogen Bonds in Neutral Histidine. Theor. Chem. Acc. 2020, 139, 125. [Google Scholar] [CrossRef]
  105. Verma, N.; Tao, Y.; Kraka, E. Systematic Detection and Characterization of Hydrogen Bonding in Proteins via Local Vibrational Modes. J. Phys. Chem. B 2021, 125, 2551–2565. [Google Scholar] [CrossRef]
  106. Beiranvand, N.; Freindorf, M.; Kraka, E. Hydrogen Bonding in Natural and Unnatural Base Pairs - Explored with Vibrational Spectroscopy. Molecules 2021, 26, 2268. [Google Scholar] [CrossRef]
  107. Makoś, M.Z.; Freindorf, M.; Sethio, D.; Kraka, E. New Insights into Fe–H2 and Fe–H- Bonding of a [NiFe] Hydrogenase Mimic—A Local Vibrational Mode Study. Theor. Chem. Acc. 2019, 138, 76. [Google Scholar] [CrossRef]
  108. Freindorf, M.; Yannacone, S.; Oliveira, V.; Verma, N.; Kraka, E. Halogen Bonding Involving I2 and d8 Transition-Metal Pincer Complexes. Crystals 2021, 11, 373. [Google Scholar] [CrossRef]
  109. Luo, Y.R. Comprehensive Handbook of Chemical Bond Energies; Taylor and Francis: Boca Raton, FL, USA, 2007. [Google Scholar]
  110. Moltved, K.A.; Kepp, K.P. Chemical Bond Energies of 3d Transition Metals Studied by Density Functional Theory. J. Chem. Theory Comput. 2018, 14, 3479–3492. [Google Scholar] [CrossRef] [Green Version]
  111. Morse, M.D. Predissociation measurements of bond dissociation energies. Acc. Chem. Res. 2018, 52, 119–126. [Google Scholar] [CrossRef]
  112. Stasyuk, O.A.; Sedlak, R.; Guerra, C.F.; Hobza, P. Comparison of the DFT-SAPT and canonical EDA Schemes for the energy decomposition of various types of noncovalent interactions. J. Chem. Theory Comput. 2018, 14, 3440–3450. [Google Scholar] [CrossRef]
  113. Levine, D.S.; Head-Gordon, M. Energy decomposition analysis of single bonds within Kohn-Sham density functional theory. Proc. Natl. Acad. Sci. USA 2017, 114, 12649–12656. [Google Scholar] [CrossRef] [Green Version]
  114. Zhao, L.; Hermann, M.; Schwarz, W.H.E.; Frenking, G. The Lewis electron-pair bonding model: Modern energy decomposition analysis. Nat. Rev. Chem. 2019, 3, 49–62. [Google Scholar] [CrossRef]
  115. Zhao, L.; Pan, S.; Holzmann, N.; Schwerdtfeger, P.; Frenking, G. Chemical Bonding and Bonding Models of Main-Group Compounds. Chem. Rev. 2019, 119, 8781–8845. [Google Scholar] [CrossRef] [PubMed]
  116. Gavezzotti, A. Comparing the strength of covalent bonds, intermolecular hydrogen bonds and other intermolecular interactions for organic molecules: X-ray diffraction data and quantum chemical calculations. New J. Chem. 2016, 40, 6848–6853. [Google Scholar] [CrossRef]
  117. Cremer, D.; Kraka, E. From Molecular Vibrations to Bonding, Chemical Reactions, and Reaction Mechanism. Curr. Org. Chem. 2010, 14, 1524–1560. [Google Scholar] [CrossRef]
  118. Kalescky, R.; Kraka, E.; Cremer, D. Are Carbon-Halogen Double and Triple Bonds Possible? Int. J. Quantum Chem. 2014, 114, 1060–1072. [Google Scholar] [CrossRef]
  119. Kalescky, R.; Zou, W.; Kraka, E.; Cremer, D. Quantitative Assessment of the Multiplicity of Carbon-Halogen Bonds: Carbenium and Halonium Ions with F, Cl, Br, and I. J. Phys. Chem. A 2014, 118, 1948–1963. [Google Scholar] [CrossRef]
  120. Oliveira, V.; Kraka, E.; Cremer, D. Quantitative Assessment of Halogen Bonding Utilizing Vibrational Spectroscopy. Inorg. Chem. 2016, 56, 488–502. [Google Scholar] [CrossRef]
  121. Setiawan, D.; Sethio, D.; Cremer, D.; Kraka, E. From Strong to Weak NF Bonds: On the Design of a New Class of Fluorinating Agents. Phys. Chem. Chem. Phys. 2018, 20, 23913–23927. [Google Scholar] [CrossRef]
  122. Sethio, D.; Oliveira, V.; Kraka, E. Quantitative Assessment of Tetrel Bonding Utilizing Vibrational Spectroscopy. Molecules 2018, 23, 2763. [Google Scholar] [CrossRef] [Green Version]
  123. Andrés, J.; Ayers, P.W.; Boto, R.A.; Carbó-Dorca, R.; Chermette, H.; Cioslowski, J.; Contreras-García, J.; Cooper, D.L.; Frenking, G.; Gatti, C.; et al. Nine questions on energy decomposition analysis. J. Comput. Chem. 2019, 40, 2248–2283. [Google Scholar] [CrossRef] [PubMed]
  124. Kraka, E.; Cremer, D. Characterization of CF Bonds with Multiple-Bond Character: Bond Lengths, Stretching Force Constants, and Bond Dissociation Energies. ChemPhysChem 2009, 10, 686–698. [Google Scholar] [CrossRef] [PubMed]
  125. Kaupp, M.; Danovich, D.; Shaik, S. Chemistry is about energy and its changes: A critique of bond-length/bond-strength correlations. Coord. Chem. Rev. 2017, 344, 355–362. [Google Scholar] [CrossRef]
  126. Bader, R.F.W. A Bond Path: A Universal Indicator of Bonded Interactions. J. Phys. Chem. A 1998, 102, 7314–7323. [Google Scholar] [CrossRef]
  127. Cremer, D.; Kraka, E. Chemical Bonds without Bonding Electron Density? Does the Difference Electron-Density Analysis Suffice for a Description of the Chemical Bond? Angew. Chem. Int. Ed. 1984, 23, 627–628. [Google Scholar] [CrossRef]
  128. Cremer, D.; Kraka, E. A Description of the Chemical Bond in Terms of Local Properties of Electron Density and Energy. Croat. Chem. Acta 1984, 57, 1259–1281. [Google Scholar]
  129. Kraka, E.; Cremer, D. Chemical Implication of Local Features of the Electron Density Distribution. In Theoretical Models of Chemical Bonding. The Concept of the Chemical Bond; Maksic, Z.B., Ed.; Springer Verlag: Berlin/Heidelberg, Germany, 1990; Volume 2, pp. 453–542. [Google Scholar]
  130. Ruedenberg, K. The physical nature of the chemical bond. Rev. Mod. Phys. 1962, 44, 326–376. [Google Scholar] [CrossRef]
  131. Levine, D.; Head-Gordon, M. Clarifying the quantum mechanical origin of the covalent chemical bond. Nat. Commun. 2020, 11, 4893. [Google Scholar] [CrossRef]
  132. Chai, J.D.; Head-Gordon, M. Long-Range Corrected Hybrid Density Functionals with Damped Atom–Atom Dispersion Corrections. Phys. Chem. Chem. Phys. 2008, 10, 6615–6620. [Google Scholar] [CrossRef] [Green Version]
  133. Chai, J.D.; Head-Gordon, M. Systematic Optimization of Long-Range Corrected Hybrid Density Functionals. J. Chem. Phys. 2008, 128, 084106. [Google Scholar] [CrossRef]
  134. Kendall, R.A.; Dunning, T.H., Jr.; Harrison, R.J. Electron Affinities of the First-Row Atoms Revisited. Systematic Basis Sets and Wave Functions. J. Chem. Phys. 1992, 96, 6796–6806. [Google Scholar] [CrossRef] [Green Version]
  135. Kendall, R.A.; Dunning, T.H., Jr.; Harrison, R.J. Gaussian basis sets for use in correlated molecular calculations. IX. The atoms gallium through krypton. J. Chem. Phys. 1999, 110, 7667–7676. [Google Scholar]
  136. Sato, T.; Nakai, H. Density Functional Method Including Weak Interactions: Dispersion Coefficients Based On the Local Response Approximation. J. Chem. Phys. 2009, 131, 224104. [Google Scholar] [CrossRef] [PubMed]
  137. López, C.S.; Faza, O.N.; Freindorf, M.; Kraka, E.; Cremer, D. Solving the Pericyclic-Pseudo pericyclic Puzzle in the Ring-Closure Reactions of 1,2,4,6-Heptatetraene Derivatives. J. Org. Chem. 2015, 81, 404–414. [Google Scholar] [CrossRef]
  138. 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; Gaussian Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  139. Zou, W.; Tao, Y.; Freindorf, M.; Makoś, M.Z.; Verma, N.; Cremer, D.; Kraka, E. Local Vibrational Mode Analysis (LModeA; Computational and Theoretical Chemistry Group (CATCO), Southern Methodist University: Dallas, TX, USA, 2022. [Google Scholar]
  140. Keith, T.A. AIMAll Version 17.11.14; K Gristmill Software: Overland Park, KS, USA, 2017. [Google Scholar]
Figure 1. Sketches of Group I–Group VI members and reference compounds investigated in this study. DHB distances (in Å) are given in blue. Short names of the complexes used throughout the manuscript are given in black.
Figure 1. Sketches of Group I–Group VI members and reference compounds investigated in this study. DHB distances (in Å) are given in blue. Short names of the complexes used throughout the manuscript are given in black.
Molecules 28 00263 g001aMolecules 28 00263 g001b
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Freindorf, M.; McCutcheon, M.; Beiranvand, N.; Kraka, E. Dihydrogen Bonding—Seen through the Eyes of Vibrational Spectroscopy. Molecules 2023, 28, 263. https://doi.org/10.3390/molecules28010263

AMA Style

Freindorf M, McCutcheon M, Beiranvand N, Kraka E. Dihydrogen Bonding—Seen through the Eyes of Vibrational Spectroscopy. Molecules. 2023; 28(1):263. https://doi.org/10.3390/molecules28010263

Chicago/Turabian Style

Freindorf, Marek, Margaret McCutcheon, Nassim Beiranvand, and Elfi Kraka. 2023. "Dihydrogen Bonding—Seen through the Eyes of Vibrational Spectroscopy" Molecules 28, no. 1: 263. https://doi.org/10.3390/molecules28010263

APA Style

Freindorf, M., McCutcheon, M., Beiranvand, N., & Kraka, E. (2023). Dihydrogen Bonding—Seen through the Eyes of Vibrational Spectroscopy. Molecules, 28(1), 263. https://doi.org/10.3390/molecules28010263

Article Metrics

Back to TopTop