Ab Initio Study of the Structural, Energetic, Bonding, and IR Spectroscopic Properties of Complexes with Dihydrogen Bonds

Dihydrogen Bonding-Seen through the Eyes of Vibrational Spectroscopy

In this work, we analyzed five groups of different dihydrogen bonding interactions and hydrogen clusters with an H3+ 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 H3+(H2)n systems, we conclude that the bond strength in these crystal-like structures makes them interesting for potential hydrogen storage applications.

Prototypical cyclohexane dimers: spectroscopic evidence for σ stacking at low temperatures

In this work, the first unambiguous spectroscopic evidence for the existence of σ stacking interactions in cyclohexane dimers has been provided using matrix isolation infrared spectroscopy.

The borazine dimer: the case of a dihydrogen bond competing with a classical hydrogen bond

Dimers of borazine were studied using matrix isolation infrared spectroscopy and ab initio quantum chemical calculations. Computations were performed at the MP2 and M06-2X levels of theory using the 6-311++G(d,p) and aug-cc-pVDZ basis sets for the various homodimers. At both levels of theory, an aligned stacked structure was found to be the global minimum, which was nearly isoenergetic to a parallel displaced structure. A T-shaped structure, where the N-H of one borazine pointed towards the N of the second borazine, was found to be a local minimum. In addition to these structures, a dihydrogen bonded structure, where the hydrogen attached to the nitrogen of borazine interacted with the hydrogen attached to the boron atom of another borazine, was also indicated. Experimentally, we observed the T-shaped dimer and the dihydrogen bonded dimer. This is one of the rare examples of experimental evidence for a dihydrogen bond, in a system other than in a metal hydride. These results for the borazine dimer were clearly different from the benzene dimer where the parallel displaced structure was found to be the global minimum followed by the T-shaped structure at the MP2/aug-cc-pVDZ level of theory. AIM, EDA and NBO analyses were carried out for all the structures to explore the nature of interactions.

Hydrogen and Dihydrogen Bonds in the Reactions of Metal Hydrides

The dihydrogen bond-an interaction between a transition-metal or main-group hydride (M-H) and a protic hydrogen moiety (H-X)-is arguably the most intriguing type of hydrogen bond. It was discovered in the mid-1990s and has been intensively explored since then. Herein, we collate up-to-date experimental and computational studies of the structural, energetic, and spectroscopic parameters and natures of dihydrogen-bonded complexes of the form M-H···H-X, as such species are now known for a wide variety of hydrido compounds. Being a weak interaction, dihydrogen bonding entails the lengthening of the participating bonds as well as their polarization (repolarization) as a result of electron density redistribution. Thus, the formation of a dihydrogen bond allows for the activation of both the MH and XH bonds in one step, facilitating proton transfer and preparing these bonds for further transformations. The implications of dihydrogen bonding in different stoichiometric and catalytic reactions, such as hydrogen exchange, alcoholysis and aminolysis, hydrogen evolution, hydrogenation, and dehydrogenation, are discussed.

Dihydrogen bond interactions as a result of H2 cleavage at Cu, Ag and Au centres

A quantum chemical study of H2 activation at fluorides of coinage metals, MF (M = Cu, Ag and Au), and its splitting was performed. The following reaction path was analyzed: FMH2→ FHHM → HMFH, where both the molecular complexes and the corresponding transition states have been characterized at the CCSD(T)/aug-cc-pVQZ//MP2/aug-cc-pVQZ level of theory. Further single-point CASSCF/CASPT2 calculations, including spin-orbit coupling effects, were also performed to analyze the role of non-dynamic correlation. The scalar relativistic effects are included via aug-cc-pVQZ-PP basis sets used for the metals. The dihydrogen-bonded copper (FHHCu) and silver (FHHAg) complexes are observed as a result of H2 cleavage, while the corresponding FHHAu gold complex is not found but the HAuHF arrangement is observed, instead. The energetic and geometrical parameters of the complexes have been analyzed and both the Quantum Theory of Atoms in Molecules approach and the Natural Bond Orbitals method were additionally applied to analyze the intermolecular interactions.

First Example of Hydrogen Bonding to Platinum Hydride

The interaction between various proton donors (indole, CF3CH2OH, (CF3)2CHOH, (CF3)3COH) and (NNC)PtH hydrido complex (NNC-H=6-(1,1'-dimethylbenzyl)-2,2'-bipyridine) was investigated through low-temperature IR and NMR spectroscopy in combination with density functional theory calculations at the M06 level of theory. The experiment shows formation of very weak hydrogen bonded complexes (Δ HºHB ca. -1.0 kcal mol-1), which undergo subsequent proton transfer surprisingly easy. Computational analysis of the hydrogen bonded complexes geometry, electronic parameters (obtained by NBO and AIM analysis), and orbital interaction energies shows that all the complexes are better described as bonded to the metal atom. At that in case of weak alcohols (CH3OH, TFE) there is also the additional interaction with the hydride ligand. These computational results allow explaining the observed experimental trends and give the first example of hydrogen bonding to a metal atom in the presence of hydride ligand.

DENSITY FUNCTIONAL THEORY STUDY OF BINDING ENERGIES, 7Li NUCLEAR MAGNETIC SHIELDING, AND ELECTRIC FIELD GRADIENT TENSORS ON THE SMALL CLUSTERS OF LinHm (m ≤ n ≤ 4)

The binding energies, geometries, 7 Li magnetic shielding, and electric field gradient tensors of hydrogenated lithium clusters, Li n H m (m ≤ n ≤ 4), were studied via density functional theory approach. We optimized the structures using B3LYP functional and 6-311++G (2d,2p) basis set. The calculated binding energies of lithium hydride clusters indicate that hydrogenation energy of Li n H m clusters decreases as the number of hydrogen atoms within the cluster increases. Our calculations also showed that for n = 4 clusters, the three-dimensional structure is more stable than the planar one. The study of the trends in the 7 Li magnetic shielding isotropy, σiso, and anisotropies, Δσ, values are explained in terms of the interplay between the electronic and geometrical effects. The variations in the 7 Li nuclear quadrupole coupling constants, χ, and their associated asymmetry parameters, ηQ, for different isomers of the lithium hydride clusters and the influence of hydrogenation on the EFG tensors are also discussed. For n = 4, we obtained a noticeable difference in the χ value from the planar to the three-dimensional structures. The atoms in molecules (AIM) analysis at the Li–H bond critical point reveals remarkably different topographical properties of the charge density and associated Laplacian fields for the planar and three-dimensional lithium hydride clusters.

The B-H...H-P dihydrogen bonding in ion pair complexes [(CF(3))(3)BH(-)][HPH(3-n)(Me)(n)(+)] (n = 0-3) and its implication in H(2) elimination and activation reactions

The B-H(delta-)...(delta+)H-P dihydrogen bonding (DHB) in ion pair complexes [(CF(3))(3)BH(-)][HPH(3-n)(Me)(n)(+)] (n = 0-3) and its role in the combination of proton and hydride with the release of H(2) or, reversibly, the heterolytic activation of H(2) by Lewis pairs (CF(3))(3)BPH(3-n)(Me)(n) have been theoretically investigated at the MP2 and DFT levels. It is found that the B-H...H-P bonds behave similarly to those in neutral pairs and ion-molecule complexes in most respects, such as the linearity of the H...H-P moiety, the characteristics of the electron transfer and rearrangement, and the topological properties of the DHB critical point, except that in certain cases, a blue-shifting of the H-bond vibrational frequency is observed. In [(CF(3))(3)BH(-)][HPH(3-n)(Me)(n)(+)], the proton shifting within the complexes leads to the formation of the dihydrogen complex B(CF(3))(3)(eta(2)-H(2)), which is followed by a subsequent H(2) release. The stability of B(CF(3))(3)(eta(2)-H(2)) (D(e)/D(0) = 10.8/6.0 kcal/mol) makes the proton-hydride combination proceed in a fashion similar to the protonation reactions in transition-metal hydrides rather than those in group 13 hydrides EH(4)(-) (E = B, Al, Ga). As for the H(2)-splitting reaction R(3)BPR'(3) + H(2) --> [R(3)BH(-)][HPR'(3)(+)], classical Lewis pair (CLP) (CF(3))(3)BPH(3) exhibits a high barrier and results in an unstable ion pair product [(CF(3))(3)BH(-)][HPH(3)(+)] compared with the "frustrated Lewis pair" (FLP) (C(6)F(5))(3)BP(tBu)(3). A detailed analysis of the mechanistic aspects of H(2) activation by (CF(3))(3)BPH(3) and (C(6)F(5))(3)BP(tBu)(3), supported by another CLP (CF(3))(3)BP(tBu)(3) which has a binding energy comparable to (CF(3))(3)BPH(3) but a reaction exothermicity comparable to (C(6)F(5))(3)BP(tBu)(3), allows us to suggest that the low stability of FLP (C(6)F(5))(3)BP(tBu)(3) is the determining factor for the low reaction barrier. The relative stability and other properties of the ion pair products [R(3)BH(-)][HPR'(3)(+)] have also been analyzed. Results strongly support the view proposed by Rokob et al. [ Rokob , T. A. ; Hamza , A. ; Stirling , A. ; Soos , T. ; Papai , I. Angew. Chem., Int. Ed. 2008 , 47 , 2435 ] that the frustration energy lowers the energy barrier and increases the exothermicity of the reaction.

Dihydrogen bonding: donor-acceptor bonding (AH...HX) versus the H2 molecule (A-H2-X)

Dihydrogen bonds (DHBs) play a role in, among others, crystal packing, organometallic reaction mechanisms, and potential hydrogen-storage materials. In this work we have analyzed the central H-H bond in linear H(4), LiH...HX, BH(4)(-)...HX, and AlH(4)(-)...HX complexes with various X by using the quantitative molecular orbital model contained in Kohn-Sham density functional theory at the BP86/TZ2P level of theory. First, we address the questions of if and how one can distinguish, in principle, between a H...H donor-acceptor DHB and the formation of an H(2) molecule by using the simple H(4) model system. The results of these analyses have been used to gain an understanding of the bonding in more realistic model systems (some of which have been studied experimentally), and how this differs from the bonding in H(4).

Intermolecular HH vibrations of dihydrogen bonded complexes H3EH(-)...HOR in the low-frequency region: theory and IR spectra

The results of DFT calculations of harmonic and anharmonic frequencies of the dihydrogen bonded (DHB) complexes H 3EH (-)...HOR (E = B, Al, Ga and HOR = CH 3OH, CF 3CH 2OH) in gas phase and in low polar medium (by CPCM model) in comparison with the partners are presented. Normal coordinate analysis of the low-frequency modes was carried out to assign the new vibrations induced by DHB formation by the potential energy distribution values. Among them, the intermolecular H...H stretching vibrations only have individual modes. The influence of central atom mass and isotope and the strength of the proton donor effects were determined. The systems convenient for IR studies were chosen from the calculation predictions. The spectral investigation was made on the BH 4 (-)/ROH complexes (ROH = CH 2FCH 2OH (MFE), CF 3CH 2OH (TFE), (CF 3) 2CHOH (HFIP)). The results of temperature dependence, isotope substitution, and influence of the proton-donor strength studies agree with the theoretical conclusions. Combination of experimental and theoretical approaches allowed determining for the first time the intermolecular stretching mode characterizing intrinsic DHB vibrations.

Novel infrared spectra for intermolecular dihydrogen bonding of the phenol-borane-trimethylamine complex in electronically excited state

The intermolecular dihydrogen bonding in the electronically excited states of the dihydrogen-bonded phenol-BTMA complex in gas phase was theoretically investigated using the time-dependent density functional theory method for the first time. It was theoretically demonstrated that the S(1) state of the dihydrogen-bonded phenol-BTMA complex is a locally excited state, in which only the phenol moiety is electronically excited. The infrared spectra of the dihydrogen-bonded phenol-BTMA complex in ground state and the S(1) state were calculated at both the O-H and B-H stretching vibrational regions. A novel infrared spectrum of the dihydrogen-bonded phenol-BTMA complex in the electronically excited state was found. The stretching vibrational absorption bands of the dihydrogen-bonded O-H and B-H groups are very strong in the ground state, while they are disappeared in the S(1) state. At the same time, a new strong absorption band appears at the C[Double Bond]O stretching region. From the calculated bond lengths, it was found that both the O-H and B-H bonds in the dihydrogen bond O-H...H-B are significantly lengthened in the S(1) state of the dihydrogen-bonded phenol-BTMA complex. However, the C-O bond in the phenol moiety is markedly shortened in the excited state, and then has the characteristics of C[Double Bond]O group. Furthermore, it was demonstrated that the intermolecular dihydrogen bonds in the electronically excited state of the dihydrogen-bonded phenol-BTMA complex are strengthened, since calculated H...H distance is drastically shortened in the S(1) state.

Electron density topological properties are useful to assess the difference between hydrogen and dihydrogen complexes

B3LYP/6-31++G(d,p) and MP2/6-31++G(d,p) calculations for a series of hydrogen- and dihydrogen-bonded systems have been carried out in order to analyze the topology of the electron density and the energy densities at the respective energy-optimized bond critical points. Even though there are no significant differences when these properties are represented as a function of the dimerization energy, they can be separated into two well-defined sets if those properties are correlated with intermolecular distances. When analyzing the dependence of various properties with equilibrium bond lengths, the specific trends of dihydrogen bond systems consist of (a) lower electron density at the bond critical point, and (b) lower concentration/depletion of that density which can be translated in a different behavior for the Laplacian components. Furthermore, the sets of molecules form two different plots which allow for a valuable classification between hydrogen- and dihydrogen-bonded systems.

Proton-transfer and H2-elimination reactions of main-group hydrides EH4- (E = B, Al, Ga) with alcohols

The reaction of the isostructural anions of group 13 hydrides EH4- (E = B, Al, Ga) with proton donors of different strength (CH3OH, CF3CH2OH, and CF3OH) was studied with different theoretical methods [DFT/B3LYP and second-order Møller-Plesset (MP2) using the 6-311++G(d,p) basis set]. The results show the general mechanism of the reaction: the dihydrogen-bonded (DHB) adduct (EH...HO) formation leads through the activation barrier to the next concerted step of H2 elimination and alkoxo product formation. The structures, interaction energies (calculated by different approaches including the energy decomposition analysis), vibrational E-H modes, and electron-density distributions were analyzed for all of the DHB adducts. The transition state (TS) is the dihydrogen complex stabilized by a hydrogen bond with the anion [EH3(eta2-H2)...OR-]. The single exception is the reaction of BH4- with CF3OH exhibiting two TSs separated by a shallow minimum of the BH3(eta2-H2)...OR- intermediate. The structures and energies of all of the species were calculated, leading to the establishment of the potential energy profiles for the reaction. A comparison is made with the mechanism of the proton-transfer reaction to transition-metal hydrides. The solvent influence on the stability of all of the species along the reaction pathway was accounted for by means of polarizable conductor calculation model calculations in tetrahydrofuran (THF). Although in THF the DHB intermediates, the TSs, and the products are destabilized with respect to the separated reactants, the energy barriers for the proton transfer are only slightly affected by the solvent. The dependence of the energies of the DHB complexes, TSs, and products as well as the energy barriers for the H2 release on the central atom and the proton donor strength is also discussed.

Electron bridging dihydrogen bond in the imidazole-contained anion derivatives

The large contact distance of electron bridging dihydrogen bond (EBDB), which is over 2.4 A, is the most prominent characteristic for the imidazole-contained anion derivatives. The elongation of N-H bond and the shortening of H...H distance can be observed upon hydration and hydrogenation. Transformation from EBDB to dissociative H2 is convenient upon sequential hydrogenation. The H...H distance decreases with the enhancement of the electronegativity of the heavy atom which contacts directly with one of these two hydrogen atoms. NMR shielding of the bonding N varies significantly upon hydration and hydrogenation. The spin-spin coupling constants, 1J(H-H), is dominated predominantly by the paramagnetic spin-orbit and diamagnetic spin-orbit contributions instead of the Fermi-contact term. Enhancement of electronegativity of the heavy atom leads to the increase of 1J(H-H) coupling constants. The stabilization is enhanced upon hydration predominantly for the formation of O-H...N H bond, while it is reversed upon hydrogenation for the cleavage of big pi bond, Pi5(6). Enhancement of the stability is demonstrated by the increase of stabilization energy and vertical electron detachment energy with the electronegativity of the heavy atom. The dominant contributions for the formation of such electron bridging dihydrogen bond are the high polarity of each fragment, large electron density between two fragments, and strong bonding interaction of the bridging electron with H(N) atoms. The H...H interaction can be formed by X-Hdelta+ and Hdelta- -Y polar molecules in Hdelta+...Hdelta- and Hdelta+...e...Hdelta+ of two forms.

The Properties of Weak and Strong Dihydrogen-Bonded DH⋅⋅⋅HA Complexes

The properties of six dihydrogen-bonded (DHB) dimers with the BeH2 molecule as a proton acceptor were calculated by MP2, CCSD(T) and B3LYP methods. The structural, energetic and spectroscopic parameters are presented and analyzed in terms of their possible correlation with the interaction energy and the intermolecular H...H separation. The symmetry-adapted perturbation theory (SAPT) calculations were performed to gain more insight into the nature of the H...H interactions. The studied complexes are divided into three groups based on the calculated intermolecular distances and the interaction energies which range from approximately -1 to -42 kJ mol(-1). The analysis of the interaction energy components indicates that, in contrast to conventional hydrogen bonds, the induction energy is the most important term in the BeH2NH4+ complex. On the other hand, there is no sharp boundary between the DHB complexes classified as hydrogen bonded and van der Waals systems. The complexation-induced changes in vibrational frequencies and in proton shielding constants show a relationship with the interaction energy. The values of the 2hJXH and 3hJBeX coupling constants correlate well with the interaction energy and with the intermolecular distance.

Computational study of the process of hydrogen bond breaking: The case of the formamide–formic acid complex

MP2/6-311++G(d,p) and B3LYP/6-311++G(d,p) quantum calculations are used to study the formamide-formic acid complex (FFAC), a system bound by two hydrogen bonds, N--H...O and O--H...O, that forms a bond ring at equilibrium. When the intermolecular separation between monomers R increases, this ring opens at a distance for which the weaker N--H...O bond breaks remaining the stronger O--H...O bond. The computational study characterizes that process addressing changes of interaction energy DeltaE, structure and properties of the electron density rho(r) as well as spatial distributions of rho(r), the electrostatic potential U(r), and the electron localization function eta(r). It is shown that the spatial derivatives of DeltaE, the topology of rho(r), and qualitative changes noticed in U(r) = 0 isocontours allow to identify a precise distance R for which one can say the N--H...O hydrogen bond has broken. Both levels of theory predict essentially the same changes of structure and electron properties associated to the process of breaking and virtually identical distances at which it takes place.

Theoretical Investigation of the Dihydrogen Bond Linking MH2 with HCCRgF (M = Zn, Cd; Rg = Ar, Kr)

An ab initio computational study of the properties of four linear dihydrogen-bonded complexes pairing MH2 (M = Zn, Cd) with HCCRgF (Rg = Ar, Kr) was undertaken at the MP2/DGDZVP level of theory. The calculated complexation energies of the linear complexes vary between 6.5 kJ/mol for M = Zn to 8.5 kJ/mol for M = Cd. Equilibrium interatomic H...H distances are roughly 2.07 A for all four complexes. The red shifts of the H-C stretching frequency of HCCRgF correlate nicely with the interaction energies.