NMR and INS Line Shapes of Transition Metal Hydrides in the Presence of Coherent and Incoherent Dihydrogen Exchange

Quantum Dynamical Investigation of Dihydrogen-Hydride Exchange in a Transition-Metal Polyhydride Complex

The ongoing shift toward clean, sustainable energy is a primary driving force behind hydrogen fuel research. Safe and effective storage of hydrogen is a major challenge (particularly for mobile applications) and requires a detailed understanding of the atomic level interactions of hydrogen with its host materials. The light mass of hydrogen, however, implies that quantum effects are important, so a quantum dynamical treatment is required to properly account for these effects in computational simulations. As one such example, we describe herein the hydrogen exchange dynamics between a hydride and a dihydrogen ligand in the [FeH(H2)(PH3)4]+ model complex. A global three-dimensional (3D) potential energy surface (PES) was constructed by fitting to and interpolating from a discrete set of grid points computed using density functional theory; exact quantum dynamical calculations were then carried out on the 3D PES using discrete variable representation basis sets. Energy levels and their quantum tunneling splittings were computed up to 3000 cm-1 above the ground state. Within that energy range, all three fundamentals have been identified using wave function plots, as well as the first three overtones of the exchange (reaction coordinate) motion and several of its combination bands. From the tunneling splittings, the Boltzmann-averaged tunneling rates were computed. The Arrhenius plot of the total exchange rate shows a clear transition around 150 K, below which the activation energy is essentially zero and above which it is less than half of the electronic structure barrier. This indicates that exchange rates are governed by quantum tunneling throughout the relevant temperature range with the low-temperature regime dominated by a single quantum (ground) state. This work is the first-ever fully quantum dynamical study to investigate the hydrogen exchange dynamics between hydride and dihydrogen ligands coordinated to a transition-metal complex.

Zero-field nuclear magnetic resonance of chemically exchanging systems

Zero- to ultralow-field (ZULF) nuclear magnetic resonance (NMR) is an emerging tool for precision chemical analysis. In this work, we study dynamic processes and investigate the influence of chemical exchange on ZULF NMR J-spectra. We develop a computational approach that allows quantitative calculation of J-spectra in the presence of chemical exchange and apply it to study aqueous solutions of [¹⁵N]ammonium (¹⁵N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{H}}_4^ +$$\end{document}H4+) as a model system. We show that pH-dependent chemical exchange substantially affects the J-spectra and, in some cases, can lead to degradation and complete disappearance of the spectral features. To demonstrate potential applications of ZULF NMR for chemistry and biomedicine, we show a ZULF NMR spectrum of [2-¹³C]pyruvic acid hyperpolarized via dissolution dynamic nuclear polarization (dDNP). We foresee applications of affordable and scalable ZULF NMR coupled with hyperpolarization to study chemical exchange phenomena in vivo and in situations where high-field NMR detection is not possible to implement.

Zero-field nuclear magnetic resonance can identify species and collective behaviors in mixtures without applied magnetic fields. Here the authors demonstrate its use for resolving proton exchange in ammonium and for the detection of hyperpolarized pyruvic acid, an important imaging biomarker.

Cross-term Splittings Due to the Orientational Inequivalence of Proton Magnetic Shielding Tensors: Do Water Molecules Trapped in Crystals Hop or Tunnel?

Water molecules trapped in crystals of barium chlorate monohydrate have been investigated by magic-angle spinning (MAS) proton NMR spectroscopy in the temperature range 110-300 K. At high temperatures, a single spinning sideband pattern is observed. Below 150 K, however, two interleaved spinning sideband manifolds appear, with distinct centerbands that do not coincide with the average isotropic chemical shift seen at high temperatures. This hitherto unknown "cross-term splitting" results from the interplay of the homonuclear dipole-dipole coupling and two anisotropic proton shielding tensors that have identical principal components but nonequivalent orientations. The resulting cross terms cannot be averaged out by rotation about the magic angle. The analysis of the exchange-induced broadening, coalescence, and narrowing of the cross-term splitting in MAS spectra allows one to estimate the rate of exchange of the two protons between 140 and 190 K. The experimental data is compared with ²H and ¹H NMR studies of the same sample reported in the literature. Density functional theory methods are utilized to estimate the thermal activation energy for a 2-fold hopping process of proton exchange about the H-O-H bisector. The Bell-Limbach model allows one to take into account contributions due to incoherent quantum tunneling in the low-temperature regime.

A quantum dynamical study of the rotation of the dihydrogen ligand in the Fe(H)2(H2)(PEtPh2)3 coordination complex

Progress in the hydrogen fuel field requires a clear understanding and characterization of how materials of interest interact with hydrogen. Due to the inherently quantum mechanical nature of hydrogen nuclei, any theoretical studies of these systems must be treated quantum dynamically. One class of material that has been examined in this context are dihydrogen complexes. Since their discovery by Kubas in 1984, many such complexes have been studied both experimentally and theoretically. This particular study examines the rotational dynamics of the dihydrogen ligand in the Fe(H)₂(H₂)(PEtPh₂)₃ complex, allowing for full motion in both the rotational degrees of freedom and treating the quantum dynamics (QD) explicitly. A "gas-phase" global potential energy surface is first constructed using density functional theory with the Becke, 3-parameter, Lee-Yang-Parr functional; this is followed by an exact QD calculation of the corresponding rotation/libration states. The results provide insight into the dynamical correlation of the two rotation angles as well as a comprehensive analysis of both ground- and excited-state librational tunneling splittings. The latter was computed to be 6.914 cm^-1-in excellent agreement with the experimental value of 6.4 cm^-1. This work represents the first full-dimensional ab initio exact QD calculation ever performed for dihydrogen ligand rotation in a coordination complex.

Gas phase 1H NMR studies and kinetic modeling of dihydrogen isotope equilibration catalyzed by Ru-nanoparticles under normal conditions: dissociative vs. associative exchange

Exposure of surface H-containing Ru-nanoparticles to D₂ gas produces HD via associative adsorption, surface H-transfer and associative desorption.

Spin polarization transfer mechanisms of SABRE: A magnetic field dependent study

We have investigated the magnetic field dependence of Signal Amplification By Reversible Exchange (SABRE) arising from binding of para-hydrogen (p-H2) and a substrate to a suitable transition metal complex. The magnetic field dependence of the amplification of the (1)H Nuclear Magnetic Resonance (NMR) signals of the released substrates and dihydrogen, and the transient transition metal dihydride species shows characteristic patterns, which is explained using the theory presented here. The generation of SABRE is most efficient at low magnetic fields due to coherent spin mixing at nuclear spin Level Anti-Crossings (LACs) in the SABRE complexes. We studied two Ir-complexes and have shown that the presence of a (31)P atom in the SABRE complex doubles the number of LACs and, consequently, the number of peaks in the SABRE field dependence. Interestingly, the polarization of SABRE substrates is always accompanied by the para-to-ortho conversion in dihydride species that results in enhancement of the NMR signal of free (H2) and catalyst-bound H2 (Ir-HH). The field dependences of hyperpolarized H2 and Ir-HH by means of SABRE are studied here, for the first time, in detail. The field dependences depend on the chemical shifts and coupling constants of Ir-HH, in which the polarization transfer takes place. A negative coupling constant of -7Hz between the two chemically equivalent but magnetically inequivalent hydride nuclei is determined, which indicates that Ir-HH is a dihydride with an HH distance larger than 2Å. Finally, the field dependence of SABRE at high fields as found earlier has been investigated and attributed to polarization transfer to the substrate by cross-relaxation. The present study provides further evidence for the key role of LACs in the formation of SABRE-derived polarization. Understanding the spin dynamics behind the SABRE method opens the way to optimizing its performance and overcoming the main limitation of NMR, its notoriously low sensitivity.

Quantum rotation and translation of hydrogen molecules encapsulated inside C₆₀: temperature dependence of inelastic neutron scattering spectra

The quantum dynamics of a hydrogen molecule encapsulated inside the cage of a C60 fullerene molecule is investigated using inelastic neutron scattering (INS). The emphasis is on the temperature dependence of the INS spectra which were recorded using time-of-flight spectrometers. The hydrogen endofullerene system is highly quantum mechanical, exhibiting both translational and rotational quantization. The profound influence of the Pauli exclusion principle is revealed through nuclear spin isomerism. INS is shown to be exceptionally able to drive transitions between ortho-hydrogen and para-hydrogen which are spin-forbidden to photon spectroscopies. Spectra in the temperature range 1.6≤T≤280 K are presented, and examples are given which demonstrate how the temperature dependence of the INS peak amplitudes can provide an effective tool for assigning the transitions. It is also shown in a preliminary investigation how the temperature dependence may conceivably be used to probe crystal field effects and inter-fullerene interactions.

Molecular flexibility and structural instabilities in crystalline l-methionine

We have investigated the dynamics in polycrystalline samples of l-methionine related to the structural transition at about 307K by incoherent inelastic and quasielastic neutron scattering, X-ray powder diffraction as well as ab-initio calculations. l-Methionine is a sulfur amino acid which can be considered a derivative of alanine with the alanine R-group CH3 exchanged by CH3S(CH2)2. Using X-ray powder diffraction we have observed at ~190K an anomalous drop of the c-lattice parameter and an abrupt change of the β-monoclinic angle that could be correlated to the anomalies observed in previous specific heat measurements. Distinct changes in the quasielastic region of the neutron spectra are interpreted as being due to the onset and slowing-down of reorientational motions of the CH3-S group, are clearly distinguished above 130K in crystalline l-methionine. Large-amplitude motions observed at low frequencies are also activated above 275K, while other well-defined vibrations are damped. The ensemble of our results suggests that the crystalline structure of l-methionine is dynamically highly disordered above 275K, and such disorder can be linked to the flexibility of the molecular thiol-ether group.

Immobilization and Characterization of RuCl2(PPh3)3 Mesoporous Silica SBA-3

The five-fold coordinated ruthenium (II) catalyst RuCl2(PPh3)3 is heterogenized on the surface of amine functionalized mesoporous silica material SBA-3 and studied by BET, XRD and 31P CP-MAS solid-state NMR. The spectra indicate the replacement of one or two of the triphenyl-phosphine groups in the course of the grafting. To distinguish between both possibilities two-dimensional J-resolved 31P MAS solid-state NMR combined with quantum chemical calculations is employed. The changes of the scalar coupling pattern between the neat catalyst and the immobilized catalyst suggest the replacement of two of the PPh3-groups by coordinative bonds to the amine-functions of the linker. DFT calculations reveal the absence of two chemically and magnetically equivalent phosphines and confirm the replacement of two PPh3-groups. This finding is similar to the results observed for the binding of the rhodium in the Wilkinson's catalyst, despite the different coordination of the metals (four-fold versus five-fold), the presence of two chlorine ligands and the different transition metal.

High-resolution 2H MAS NMR applied to deuterium analogs of hydrido eta2-dihydrogen complexes

(2)H solid-state, variable temperature magic angle spinning (MAS) NMR spectra of precipitated samples of the deutero dideuterium complexes Ru(D)(2)(eta(2)-D(2))(2)(PCy(3))(2) and RuD(eta(2)-D(2))I(PCy(3))(2) [Cy=cyclohexyl] are presented. They show that even at moderate MAS speed, high resolution is achieved at 7 and 14T allowing (2)H chemical shifts and quadrupole couplings to be obtained and assigned to different solid and gaseous (2)H species. These two parameters allow identifying chemically different hydrogen species in the material. The analysis of these parameters in this study reveals the presence of three different species in the sample, namely the complexes RuD(eta(2)-D(2))I(PCy(3))(2) and RuD(eta(2)-D(2))(2)I(PCy(3))(2), and highly mobile HD/D(2). These assignments are supported by (2)H T(1) relaxation times and (31)P MAS NMR spectra. Moreover, variable temperature MAS NMR spectra reveal temperature-dependent line-shape changes, which are clear indications of intramolecular hydrogen exchange of the deutero and the dideuterium ligands and which give an estimate for the activation energy of this process.

Mechanisms of Dipolar Ortho/Para-H2O Conversion in Ice

In this paper a possible explanation for an unexpected ortho/para-water ratio in the gas clouds of comets is given. The description is based on the quantum-mechanical density matrix formalism and the spin temperature concept. Only the nuclear spin system is treated quantum-mechanically. Employing the model of a four spin system, created by two nearest neighbour water molecules, spin eigenstates and their dynamics under the influence of their mutual dipolar interactions are studied. It is shown that a fast conversion between ortho- and para-states occurs on a msec time scale, caused by the intermolecular homonuclear magnetic dipolar interaction. Moreover the spin eigenstates of water in an ice crystal are determined by magnetic dipolar interactions and are not given by normal ortho- and para-H2O states of gaseous water. As a result of this the spin temperature of gaseous water evaporated from ice depends strongly on its evaporation history and the ortho/para-ratio of water molecules are only an indirect measure of the temperature of ice crystals from where they descend. This result could explain the unexpected experimentally observed ortho/para-ratios in the clouds of comets.

The sigma-CAM Mechanism: sigma complexes as the basis of sigma-bond metathesis at late-transition-metal centers

Complexes in which a sigma-H--E bond (E=H, B, Si, C) acts as a two-electron donor to the metal center are called sigma complexes. Clues that it is possible to interconvert sigma ligands without a change in oxidation state derive from C--H activation reactions effecting isotope exchange and from dynamic rearrangements of sigma complexes (see Frontispiece). Through these pathways, metathesis of M--E bonds can occur at late transition metals. We call this process sigma-complex-assisted metathesis, or sigma-CAM, which is distinct from the familiar sigma-bond metathesis (typical for d(0) metals and requiring no intermediate) and from oxidative-reductive elimination mechanisms (inherently requiring intermediates with changed oxidation states and sometimes involving sigma complexes). There are examples of sigma-CAM mechanisms in catalysis, especially for alkane borylation and isotope exchange of alkanes. It may also occur in silylation and alkene hydrogenation.

Novel Insights into the Mechanism of the Ortho/Para Spin Conversion of Hydrogen Pairs: Implications for Catalysis and Interstellar Water

The phenomenon of exchange coupling is taken into account in the description of the magnetic nuclear spin conversion between bound ortho- and para-dihydrogen. This conversion occurs without bond breaking, in contrast to the chemical spin conversion. It is shown that the exchange coupling needs to be reduced so that the corresponding exchange barrier can increase and the given magnetic interaction can effectively induce a spin conversion. The implications for related molecules such as water are discussed. For ice, a dipolar magnetic conversion and for liquid water a chemical conversion are predicted to occur within the millisecond timescale. It follows that a separation of water into its spin isomers, as proposed by Tikhonov and Volkov (Science 2002, 296, 2363), is not feasible. Nuclear spin temperatures of water vapor in comets, which are smaller than the gas-phase equilibrium temperatures, are proposed to be diagnostic for the temperature of the ice or the dust surface from which the water was released.

Mechanism of nuclear spin initiated para-H2 to ortho-H2 conversion

In this paper a quantitative explanation for a diamagnetic ortho/para H2 conversion is given. The description is based on the quantum-mechanical density matrix formalism originally developed by Alexander and Binsch for studies of exchange processes in NMR spectra. Only the nuclear spin system is treated quantum-mechanically. Employing the model of a three spin system, the reactions of the hydrogen gas with the catalysts are treated as a phenomenological rate process, described by a rate constant. Numerical calculations reveal that for nearly all possible geometrical arrangements of the three spin system an efficient spin conversion is obtained. Only in the chemically improbable case of a linear group H-X-H no spin conversion is obtained. The efficiency of the spin conversion depends strongly on the lifetime of the H-X-H complex and on the presence of exchange interactions between the two hydrogens. Even moderate exchange couplings cause a quench of the spin conversion. Thus a sufficiently strong binding of the dihydrogen to the S spin is necessary to render the quenching by the exchange interaction ineffective.

Bridging the gap between homogeneous and heterogeneous catalysis: ortho/para H(2) conversion, hydrogen isotope scrambling, and hydrogenation of olefins by Ir(CO)Cl(PPh(3))(2)

Some transition metal complexes are known to catalyze ortho/para hydrogen conversion, hydrogen isotope scrambling, and hydrogenation reactions in liquid solution. Using the example of Vaska's complex, we present here evidence by NMR that the solvent is not necessary for these reactions to occur. Thus, solid frozen solutions or polycrystalline powdered samples of homogeneous catalysts may become heterogeneous catalysts. Comparative liquid- and solid-state studies provide novel insight into the reaction mechanisms.

Structural distortions in mer-M(H)3(NO)L2 (M = Ru, Os) and their influence on intramolecular fluxionality and quantum exchange coupling

Molecules of the type mer-M(H)3(NO)L2 [M = Ru (1), Os (2); L = PR3] are characterized on the basis of 1H NMR T1min values and IR spectra as pseudo-octahedral trihydrides significantly distorted by compression of the cis H-M-H angles to approximately 75 degrees. The distortion, uncharacteristic of six-coordinate d6 complexes, is rationalized with DFT (B3LYP) calculations as being driven by increased H-to-M sigma donation and by the exceptional pi-accepting ability of linear NO+. In both 1 and 2, hydrides undergo intramolecular site exchange with delta HHH++(1) = 10-11 kcal/mol and delta HHH++(2) = 16-20 kcal/mol, depending on L, whereas for mer-Ru(H)3(NO)(PtBu2Me)2 (1b), moderate exchange couplings (up to 77 Hz) are featured in the low-temperature 1H NMR spectra, in addition to chemical exchange. On the basis of experimental and theoretical results, a dihydrogen intermediate is suggested to mediate hydride site exchange in 1. The cis H-M-H distortion shortens the tunneling path for the exchanging hydrides in 1, thereby increasing the tunneling rate; diminishes the "conflict" between trans hydrides in the mer geometry; and decreases the nucleophilicity of the hydrides. The generality of the observed structural distortion and its dependence on the ligand environment in late transition metal tri- and dihydrides are discussed. A less reducing metal center is generally characterized by greater distortion.

Probing the motion of the η2-dideuterium ligand by solution and solid-state 2H NMR spectroscopy

Variable temperature 2H NMR is used to measure the T1min values of the η2-D2 and D ligands in trans-[M(η2-D2)(D)(dppe)2]+, M = Ru (1) and Os (2) in solution. The rapid spinning motion of the η2-D2 ligand results in a much longer T1min than that of the terminal deuteride. The quadrupole coupling constant (CQ) for the terminal deuteride is calculated to be 79 kHz for 1 and 81 kHz for 2 while motion-reduced coupling constant (CQeff) for the D in the η2-D2 ligands in 1 and 2 are in the range 19-22 and 27-31 kHz, respectively. The actual CQ for these ligands with short D-D distances (<1 Å) should be at least 2CQeff and probably greater than that of the terminal deuteride but less than that of D2(g), 227 kHz. A fast spinning (>>61 MHz) and tilting of the primary electric field gradient component of between 90° and 60° or between 50° and 40° with respect to the axis of D2 rotation is an explanation for the small CQeff. Therefore neither D-D nor M-D bonding dominates the electric field gradient direction in these M(η2-D2) bonds. The complex [RuD2(C5Me5)(dppm)]+ 3, which exists in solution as a 3.3:1 mixture of Ru(D···D) (3a) and Ru(D)2 (3b) tautomers, has T1min(D) values that provide CQeff values of 66 kHz for 3a and 71 kHz for 3b. The elongated D···D ligand in 3a with d(D···D) ~ 1.10 Å may be "static" compared to the 61.4 MHz 2H spectrometer frequency and therefore have "compressed dihydride" character and display similar Ru-D bonds as in 3b. However, it is more likely that the static CQ for the D of 3a is much larger than that of a terminal Ru-D but averaged to 66 kHz by a 180° flip of the D2 as observed in 7. The complexes trans-[M(D···D)(Cl)(dppe)2]PF6, M = Ru (4), Os (5), and Ru(η2-D2)(dppb)(µ-Cl)3RuCl(dppb) (6) also have motion-reduced CQ values. Some 2H NMR quadrupole echo wide-line spectra of 1, 4, and Os(D···D)(Cl)2(CO)(PiPr3)2 (7) (contaminated with 20% Os(H···D)) were recorded in the solid phase from 293 to 123 K. These also indicate that the CQeff of the D in the D2 ligands are motion-reduced. Simulation of the spectra of 7 are suggestive of a twofold reorientation of a D2 ligand with a static CQ value in the range of 120-167 kHz, a tilt angle of the electric field gradient of about 50°, and an asymmetry parameter near to zero. Significantly, the MAS 2H spectrum of 7 has a broad doublet, possibly due to a non-averaged D-D dipolar/quadrupolar interference phenomenon; this is first time this has been observed. Variable temperature T1 data for solid 7 are also reported, which allow the evaluation of the activation barrier to the twofold flipping motion of the D2; the only other way of quantifying such an energy barrier is by use of an inelastic neutron scattering method.Key words: dihydrogen, deuterium, NMR, ruthenium, osmium, complexes, dynamics, hydride, bonding.