Probing Hydrogen Activation in a Dimetal Dihydride Complex by Symmetric Exchange with Parahydrogen

TL;DR

This study explores hydrogen exchange in a nickel complex using parahydrogen and NMR, revealing unexpected magnetization effects.

Contribution

The paper demonstrates spontaneous conversion of nuclear singlet order without asymmetric intermediates, confirmed through theory and simulations.

Findings

Spontaneous conversion of nuclear singlet order into in-phase longitudinal magnetization was observed.

Magnetization is retained during dihydrogen elimination and detected as an enhanced NMR signal.

Relaxation interference processes explain the phenomenon without requiring asymmetric intermediates.

Abstract

We report an investigation
into the hydrogen exchange mechanism
in a dinickel­(II) dihydride complex with C
2v symmetry. To obtain atomic level information on the low-concentration
intermediate species, we use parahydrogen in combination with nuclear
magnetic resonance (NMR). Unexpectedly, despite the chemical equivalence
of the two hydride sites imposed by the mirror symmetry of the model,
we observe spontaneous conversion of the nuclear singlet order into
in-phase longitudinal magnetization–something that normally
requires sophisticated pulse sequences. The magnetization is retained
upon reductive elimination of dihydrogen from the dihydride species
and observed as an enhanced hydrogen NMR signal. This is explained
by a chain of nuclear relaxation interference processes that do not
require asymmetric intermediates; we confirm this in two ways: using
analytical relaxation theory and with brute-force numerical simulations.

Figures

(A) Dinickel dihydride complex K­[L­(NiII–H)2] (1H 2 ) studied within this work showing pairwise H2/D2 exchange. 1H 2 is highly sensitive toward O2 forming the peroxido complexK­[LNiII 2(O2)] (1O 2 ) irreversibly, and toward water reacting to K­[LNiII 2(H)­(OH)] (1(H)­(OH)) and ultimately to K­[LNiII 2(μ–OH)] (1OH). Oxidation of 1H 2 enables reductive binding of N2 coupled to H2 release giving K­[LNiII 2(μ1,2–N2)] (1N 2 ); a scenario resembling the reactivity of nitrogenase metalloenzyme where binding of N2 to the FeMo metal cofactor is also paired with reductive elimination of H2 (see B). , The hydride arrangement as well as binding mode of N2 in nitrogenases is still under discussion. − Reproduced from ref . Copyright 2025 American Chemical Society.

(A) Pairwise exchange of dihydride ligands of complex 1H 2 with pH2 or thermally polarized H2 showing enhanced signals for the dihydride and for free H2. (B) Overlay of 1H-PHIP-NMR spectra of complex 1H 2 (13.4 mM) in THF-d 8 at 298 K with π/4 detection pulses before (black) and after 10 s bubbling with pH2 (magenta) showing pH2 induced enhancement of signals at 4.54 ppm and −24.1 ppm with pure in-phase magnetization. The inserts show subsections of the spectra showing sensitivity enhancements.

1H-PHIP-CEST profiles obtained for 1H 2 (13.4 mM) in THF-d 8. The value of the Ni–H integral at −24.1 ppm (green) or H2 integral at 4.54 ppm (orange) in a.u. is plotted vs the CW irradiation offset in ppm relative to tetramethyl silane. PHIP-CEST profiles were reproduced more than 3 times with different samples and spin-lock-frequency lists (see Figures S12–S15). NMR-invisible pH2 exchanges with the dihydride complex, leading to sensitivity enhancement of the dihydride signal at −24.1 ppm. CW irradiation at this offset leads to a sharp decrease in signal intensity at 4.54 ppm (and vice versa) indicative of polarized free H2 in solution. All PHIP-CEST spectra were acquired with π/4 flip-angle pulses as single-scan spectra at 9.4 T (400 MHz) using 2 s cw saturation after 10 s bubbling and allowing the sample to settle for 1 s according to Figure S11. Spin-lock-field amplitudes (|γH B 1|) were set to 50 Hz and spin-lock-field offsets (cw offset) were varied according to uniformly and non-uniformly distributed lists of 24 points from −25 to 7.5 ppm.

Numerical simulation, using Spinach 2.9, of the dynamics of the components of the singlet state (blue and red traces) and longitudinal magnetization (green trace) for a system initially in a pure singlet state and evolving in time under the action of the Hamiltonian and the Bloch-Redfield-Wangsness relaxation superoperator assuming a rotational correlation time of 0.5 ns and magnetic field of 18.79 T. The dipole-dipole coupling tensor was obtained from the energy minimum molecular geometry (CPCM THF) obtained using PBE0-D3­(BJ)/def2-TZVP method in ORCA. Proton chemical shielding tensors were then obtained using GIAO method and def2-TZVPP basis set in the same program.

Full text

Introduction

Understanding fundamental chemical processes observed in nature, such as hydrogen production and nitrogen fixation, is important for the development of more efficient and selective catalysts. ?,? This is particularly true for organometallic and enzymatic catalysis ?−? ? ? where transient chemical species can occur in low concentrations, but may still determine reaction pathways and the selectivity of the catalytic process. Direct observation and structural characterization of these intermediates using mainstream analytical methods is difficult and limited to special cases; there are very few generally applicable tools. ?−? ?

One popular technique is nuclear magnetic resonance (NMR) spectroscopy; it is highly chemically selective and can provide structural information. ?−? ? ? However, NMR has low sensitivity and often struggles to detect low-concentration reaction intermediates. One way to remedy this is parahydrogen-induced polarization (PHIP)–a method with nearly 40 years of history in the study of catalytic processes. ?−? ? It belongs to the broader group of spin hyperpolarization methods? and uses parahydrogen (pH_2_)–an easily obtainable nuclear spin isomer of molecular hydrogen.? The two protons in pH_2_ are 100% spin-polarized relative to each other, but the total spin of two antiparallel protons is zero and the state is therefore initially NMR-silent. Chemical manipulation is therefore necessary to break the symmetry and convert that relative spin polarization into much stronger observable magnetization than what is available in conventional NMR spectroscopy.?

Major work in PHIP methodology today is concentrated in biomedical applications,? such as cancer imaging, ?,? metabolomics, ?,? as well as in the studies of catalytic processes involving hydrogen activation. ?−? ? Even the mere presence of PHIP is informative: the mechanism must keep the two protons together to retain their spin correlation. ?,?,? This was initially limited to cis-hydrogenation but was later expanded to include trans-selective and geminal hydrogenation mediated by ruthenium carbene complexes. ?,? Those mechanistic studies later enabled hyperpolarization of fumarate? and the detection of fumarase activity in cells.? Since detectable magnetization is produced as soon as the pH_2_ symmetry is broken, the investigation of metal dihydrogen activation is possible as demonstrated by the detection of transient classical ?−? ? and nonclassical metal dihydride complexes, ?,? also by indirect detection through partially negative line (PANEL) experiments.? Other applications include frustrated Lewis pair molecular tweezers that function as metal-free hydrogenation catalysts. ?−? ?

Recent extensions of PHIP methodology have enabled mechanistic studies of hydrogenases under catalytic conditions, leading to the characterization of two previously elusive intermediates in the catalytic cycle of [Fe]-hydrogenases.? These metalloenzymes? are involved in H_2_ activation and production; they serve as blueprints for the development of eco-friendly catalysts.? Transient diamagnetic hydrogen-bound intermediates in other metalloenzymes are now within reach, for example nitrogenases wherein metal hydrides are present in the FeMo cofactor and used to store reducing equivalents, enabling reductive binding of dinitrogen to the Fe/S active site accompanied by hydrogen release. ?−? ? ?

In this communication, we study similar processes in dinickel dihydride complex K­[L­(Ni^II^–H)2] (1H _ 2 , where L is a compartmental ligand scaffold) that may be viewed as a model of the enzymatic reactions discussed above, because it is known to undergo replacement of the two hydride ligands with concomitant reductive binding of small molecules coupled to H_2 release ?−? ? ? ? ? (see Figure), and because its oxidation triggers even the H_2_-releasing reductive activation of N_2_.? The H_2_/D_2_ exchange reaction proceeds in a pairwise synchronous fashion without H/D scrambling,? and density functional theory (DFT) calculations and characterization of the proposed {LNi^I^} intermediate suggest that this may follow a dissociative stepwise mechanism. ?,? In this context, parahydrogen could help confirm the proposed mechanism through NMR signal enhancement in low-concentration intermediates, and could provide evidence for the dissociative exchange pathway.

*(A) Dinickel dihydride complex K­[L­(NiII–H)2] (1H

2 ) studied within this work showing pairwise H2/D2 exchange. 1H

2 is highly sensitive toward O2 forming the peroxido complexK­[LNiII 2(O2)] (1O

2 ) irreversibly, and toward water reacting to K­[LNiII 2(H)­(OH)] (1(H)­(OH)) and ultimately to K­[LNiII 2(μ–OH)] (1OH). Oxidation of 1H

2 enables reductive binding of N2 coupled to H2 release giving K­[LNiII 2(μ1,2–N2)] (1N

2 ); a scenario resembling the reactivity of nitrogenase metalloenzyme where binding of N2 to the FeMo metal cofactor is also paired with reductive elimination of H2 (see B). , The hydride arrangement as well as binding mode of N2 in nitrogenases is still under discussion. −

Reproduced from ref . Copyright 2025 American Chemical Society.*

The difficulty in using PHIP to elucidate the H_2_ binding mechanism of 1 is that the two hydrides in 1H _ 2 _ occupy chemically equivalent sites related by mirror symmetry. Accordingly, in the absence of structural perturbations, the two nuclear spins ought to remain in the singlet state, rendering them NMR silent under standard PHIP mechanisms. ?,?,?,? However, we demonstrate here that spontaneous polarization transfer from the initial parahydrogen singlet spin order to in-phase proton longitudinal magnetization does nonetheless occur. This is unexpected, but may be explained using a nuclear spin relaxation model taking into account chemical shift anisotropy (CSA) and dipole-dipole (DD) coupling interference effects–we demonstrate this both algebraically and numerically below. This CSA―CSA and DD―CSA interference mechanism bypasses the classical requirement for asymmetric intermediates. ?,?

PHIP-CEST ?,? (chemical exchange saturation transfer experiments) with 1H _ 2 _ indicates the absence of asymmetric intermediates in the exchange reaction and that dihydrogen acquires net longitudinal magnetization upon release from the complex. This provides further insight into the dynamics of the H_2_ exchange process at the molecular level. Our results highlight a novel pathway for the generation of longitudinal magnetization from parahydrogen in systems where the hydrogen atoms remain chemically equivalent, emphasizing the role of nuclear spin relaxation interference effects in the dynamics of parahydrogen singlet order.

Results

Studying air and moisture sensitive compounds requires manipulation under an inert atmosphere. Previous PASADENA-type experiments with 1H _ 2 _ were hampered by insufficiently inert conditions.? Consequently, a new PASADENA setup that fulfills these specifications and that is easy to operate and maintain is crucial. Within this article, we present the development and construction of such a setup operating under anoxic and anhydrous conditions (see SI for further details).

PHIP experiments under inert conditions were carried out by preparing samples under an argon atmosphere in a glovebox. The sample holder was subsequently transferred into the bore of the magnet. The atmosphere of the operational section of the setup was repeatedly evacuated and backfilled with inert gas (nitrogen) to ensure the exclusion of air and moisture. After this purging process, the valves to the sample holder were opened, and the experiments were initiated. pH_2_ was introduced to the sample either automatically via computer-controlled magnetic valves or manually by operating the valves by hand.

Given the molecular symmetry of 1H _ 2 _ we were intrigued to observe enhancement of signals in the ^1^H NMR spectrum after bubbling pH_2_ through a sample of 1H _ 2 _ under argon (pulse sequence is depicted in SI Figure 7). In the C 2v symmetric complex 1H _ 2 _ the two hydrides occupy chemically equivalent positions related by mirror symmetry (FigureA). Only minor amounts of water decomposition products (88% 1H _ 2 , 12% 1(H)­(OH)) were detected before bubbling. This shows that the sample holder is gastight on the time scale of the experiment and is well suited for use with highly sensitive compounds. Even after bubbling with pH_2 for 10 s the amounts of decomposition products increases only slightly (81% 1H _ 2 _, 12% 1­(H)­(OH), 7% 1O _ 2 ). After bubbling with pH_2 moderate enhancements of the signals at −24.1 ppm (I PHIP/I thermal = 9.2, 0.03% polarization) and 4.54 ppm are visible in the NMR spectrum of 1H _ 2 _ in THF-d 8 although the hydrogen atoms occupy chemically equivalent positions (see FigureB). The former signal is assigned to the dihydride position of 1H _ 2 _ whereas the latter signal assignment has some ambiguity that will be addressed below.

*(A) Pairwise exchange of dihydride ligands of complex 1H

2 with pH2 or thermally polarized H2 showing enhanced signals for the dihydride and for free H2. (B) Overlay of 1H-PHIP-NMR spectra of complex 1H

2 (13.4 mM) in THF-d 8 at 298 K with π/4 detection pulses before (black) and after 10 s bubbling with pH2 (magenta) showing pH2 induced enhancement of signals at 4.54 ppm and −24.1 ppm with pure in-phase magnetization. The inserts show subsections of the spectra showing sensitivity enhancements.*

Notably, both signals are purelyin-phase and not antiphase as expected for PHIP experiments under PASADENA conditions. ?,? To emphasize the absence of classical antiphase PASADENA signals in Figure we decided to use π/4 pulses for detection, even if π/2 pulses would yield the highest intensity for the in-phase PHIP signals observed. To verify that these signals originate from 1H _ 2 , control experiments were performed under identical conditions using some of the decomposition products. Accordingly, NMR samples of 1OH and 1O _ 2 _ were prepared independently under argon in THF-d 8 and exposed to pH_2 (7 bar). After 10 s of bubbling and a 2 s settling time, π/4 flip angle ^1^H NMR spectra were acquired (see Figures S9 and S10). No signal enhancement was detected in either experiment, confirming that 1H _ 2 _ is essential for the PHIP effects leading to signal enhancement.

By combining PHIP techniques with the chemical exchange saturation transfer (CEST), ?−? ? herein referred to as PHIP-CEST ?,?,? transiently formed species can be indirectly detected through their chemical exchange with more abundant species. This method was recently used to reveal H_2_-bound states in hydrogenase catalysis with extremely high sensitivity.? Herein, we use PHIP-CEST to determine the origin of the signal at 4.54 ppm and test the hypothesis of chemical exchange between the dinickel dihydride species (−24.1 ppm) and free hydrogen. The hydride region in the NMR is tested (−25 ppm as a high-field limit) to exclude the presence of possible intermediates on the exchange reaction pathway. When applying a radio frequency (RF) spin-lock pulse on resonance with a signal of a certain species, the corresponding NMR transition is saturated, and the signal is attenuated. In CEST, saturation is transferred from one species to another if both are related by chemical exchange. In the system under study, the signal at 4.54 ppm can be attributed either to the methylene proton of the β-ketiminato ligand backbone of 1H _ 2 _ or to free H_2_, presumably polarized by exchange with pH_2_ bound transiently to the nickel complex. To provide evidence that the signal at 4.54 ppm stems from freely dissolved H_2_, 24 ^1^H NMR spectra with π/4 flip-angle were acquired following the pulse sequence depicted in Figure S11. The saturation irradiation was applied for 2 s with a spin-lock amplitude of 50 Hz and central frequencies varied according to lists given in the SI. Accordingly, a reduction of the dihydrogen signal at 4.54 ppm is expected when irradiating at −24.1 ppm. By plotting the integral of these signals against the continuous wave (CW) irradiation offset, PHIP-CEST profiles are obtained (see Figure).

*1H-PHIP-CEST profiles obtained for 1H

2 (13.4 mM) in THF-d 8. The value of the Ni–H integral at −24.1 ppm (green) or H2 integral at 4.54 ppm (orange) in a.u. is plotted vs the CW irradiation offset in ppm relative to tetramethyl silane. PHIP-CEST profiles were reproduced more than 3 times with different samples and spin-lock-frequency lists (see Figures S12–S15). NMR-invisible pH2 exchanges with the dihydride complex, leading to sensitivity enhancement of the dihydride signal at −24.1 ppm. CW irradiation at this offset leads to a sharp decrease in signal intensity at 4.54 ppm (and vice versa) indicative of polarized free H2 in solution. All PHIP-CEST spectra were acquired with π/4 flip-angle pulses as single-scan spectra at 9.4 T (400 MHz) using 2 s cw saturation after 10 s bubbling and allowing the sample to settle for 1 s according to Figure S11. Spin-lock-field amplitudes (|γH B 1|) were set to 50 Hz and spin-lock-field offsets (cw offset) were varied according to uniformly and non-uniformly distributed lists of 24 points from −25 to 7.5 ppm.*

Discussion

When irradiating at the hydride position, the signal intensity at 4.54 ppm decreases dramatically and vice versa. This indicates that the signal at 4.54 ppm belongs to free hydrogen, which was released from 1H _ 2 _ by reversible reductive elimination, in accordance with the mechanism for the pairwise H_2_/D_2_ exchange.? Herein, pH_2_ presumably exchanges with the hydride ligands in 1H _ 2 , is converted to oH_2 (orthohydrogen) upon singlet–triplet transition, and is released in a polarized state by reductive elimination, leaving behind an unsaturated complex able to bind hydrogen again.

Sophisticated pulse sequences? and chemical transformations ?,? are normally required to convert nuclear singlet order into magnetization, but the process can also happen spontaneously through the dipole-dipole/chemical shift anisotropy (DD-CSA) interference mechanism? which opens cross-relaxation channels from the longitudinal correlation component I Z S Z of the singlet state to longitudinal magnetization states I Z and S Z. This mechanism is significant when the I ± S ∓ part of the singlet disappears sufficiently fast. ?,?

In this instance, we have a two-proton system with a J-coupling, a dipole-dipole coupling, and two noncollinear chemical shift tensors not aligned in the spin-spin direction

$$[eqn]$$

where B is the magnetic field vector (Tesla), I and S are vectors of three Cartesian projection operators for the two spins (dimensionless), Z = −γ­(1–σ) are Zeeman interaction tensors (rad/s·T), including the isotropic Zeeman interaction with a magnetogyric ratio γ and the chemical shielding σ; J is the scalar coupling constant (Hz), and D is the dipole-dipole coupling tensor (rad/s). In a strong vertical magnetic field B = [0 0 B 0]^T^, the isotropic part of this Hamiltonian

$$[eqn]$$

is symmetric with respect to the permutation of the two spins because here ω_I_ = ω_S_. The singlet state density matrix commutes with this Hamiltonian and therefore cannot evolve. However, the two Zeeman tensors Z I,S are only related by mirror symmetry–not by inversion or permutation symmetry, meaning that relaxation leakage is expected? from the singlet subspace.

Automatic symbolic evaluation? of Bloch-Redfield-Wangsness relaxation superoperator ?,?

$$[eqn]$$

in the case of isotropic rotational diffusion (see the Mathematica script in the Supporting Information) reveals that the cross-relaxation rate from the pure singlet state to I Z + S Z is zero, but that longitudinal and transverse components of the singlet have different relaxation rates

$$[eqn]$$

Here, the invariants of 3 × 3 Cartesian interaction tensors are defined as?

$$[eqn]$$

and J(ω) = τ_C_/(1 + ω^2^τ_C_ ^2^) is the isotropic rotational diffusion autocorrelation function.

The presence of J(0) terms in R[I + S – + I – S +] means that the transverse components of the singlet disappear faster than the longitudinal components. This breaks the symmetry lockout and allows longitudinal magnetization to accumulate via DD-CSA cross-correlation with the following rate

$$[eqn]$$

This process is easily visualized (Figure) in Spinach (see the Matlab script in the Supporting Information), where the Bloch-Redfield-Wangsness theory module is designed to account for all known and unknown Redfield type processes of all ranks automatically. ?,?

Numerical simulation, using Spinach 2.9, of the dynamics of the components of the singlet state (blue and red traces) and longitudinal magnetization (green trace) for a system initially in a pure singlet state and evolving in time under the action of the Hamiltonian and the Bloch-Redfield-Wangsness relaxation superoperator assuming a rotational correlation time of 0.5 ns and magnetic field of 18.79 T. The dipole-dipole coupling tensor was obtained from the energy minimum molecular geometry (CPCM THF) obtained using PBE0-D3­(BJ)/def2-TZVP method in ORCA. Proton chemical shielding tensors were then obtained using GIAO method and def2-TZVPP basis set in the same program.

Conclusion

We have observed parahydrogen-induced NMR signal enhancements for a C 2v symmetric dinickel dihydride complex that couples the reductive activation of challenging substrates to H_2_ release, akin to the proposed mechanism for nitrogenase enzymes. We report an experimental setup that enables PASADENA-type PHIP experiments for sensitive compounds under inert conditions.

The unusual spontaneous interconversion of the nuclear singlet order into in-phase longitudinal magnetization, apparently without symmetry breaking, was explained by a chain of nuclear spin relaxation interference effects: first CSA-CSA cross-correlation makes longitudinal and flip-flop components of the singlet relax at different rates, and then DD-CSA cross-correlation converts longitudinal two-spin order into longitudinal magnetization. The presence of longitudinally magnetized dihydrogen was corroborated by PHIP-CEST experiments; no asymmetric intermediates were found.

Although in general, contributions of paramagnetic intermediates cannot be excluded, there is no direct evidence for such contributions, since the DD-CSA cross-correlated relaxation mechanism explains our observations. In analogy to recent PHIP experiments employing phosphorus biradicaloids,? the putative {LNi^I^} intermediate likely also does not interfere with PHIP generation due to its singlet electronic ground state.?

Our findings suggest that parahydrogen techniques can be applied to symmetric compounds, improving NMR sensitivity and facilitating the study of symmetric molecules in catalysis and materials science. Further work on this methodology will enable the investigation of hydrogen-bound intermediates in hydrogenase and possibly nitrogenase enzymes, as well as their synthetic model systems.

Supplementary Material