## Improving Resolution of Solid State NMR in Dense Molecular Hydrogen

Thomas Meier, 1, ∗ Saiana Khandarkhaeva, 1 Jeroen Jacobs, 2 Natalia Dubrovinskaia, 3 and Leonid Dubrovinksy 1

1 Bavarian Geoinstitute, University of Bayreuth, D-95447 Bayreuth, Germany 2 European Synchrotron Radiation Facility (ESRF), F-38043 Grenoble Cedex, France 3 Material Physics and Technology at Extreme Conditions, Laboratory of Crystallography, University of Bayreuth, D-95447 Bayreuth, Germany (Dated: August 6, 2019)

Recent advancements in radio-frequency resonator designs have led to the implementation of nuclear magnetic resonance in diamond anvil cells (DACs) at pressures well above 100 GPa. However, a relatively low resolution and the absence of decoupling sequences complicate the analysis of the results of solid state NMR in DACs. Here, we present the first application of homo-nuclear LeeGoldburg (LG) decoupling at extreme conditions on high density molecular hydrogen. Lenz lens based two-dimensional resonator structures were found to generate a homogeneous B 1 field across sample cavities as small as 12 picolitres, a prerequisite for optimal decoupling. At ideal LG conditions, the broad 1 H resonance of molecular ortho-hydrogen was narrowed 1600-fold, resulting in line-widths of 3.1 ppm.

After the first prediction of a metallic phase under high compression by Wigner and Huntington[1], and an estimation of the superconducting properties at room temperature by Ashcroft[2], H 2 became an archetypical test system in modern high pressure science[3]. In 2017, experimental[4] observations (although heavily disputed[5-7] suggested the transformation of hydrogen into a metallic state at around 500 GPa. Recent IR spectroscopic results[8] indicate a band gap closure at a pressure of 425 GPa, which may be interpreted as the onset of metallisation.

never been pursued.

Molecular hydrogen remains an elusive quantum solid[9] in high pressure research. Owing to the lack of core electrons, hydrogen nuclei are practically invisible to X-ray diffraction methods[10]. The equation of state of solid hydrogen was, however, determined up to 1 Mbar following three sets of diffraction peaks (100, 002 and 101)[11]. Other structural probes, such as neutron diffraction, cannotbe applied at the necessary pressures in excess of 100 GPa[12], and complementary spectroscopic methods such as IR or Raman provide only limited information on the atomic or molecular arrangements of dense hydrogen[13]. Nuclear Magnetic Resonance spectroscopy (NMR) is an alternative to the techniques mentioned above, and is widely regarded as one of the most versatile spectroscopic methods, probing local magnetic fields of single nuclei or whole molecular building blocks within inorganic compounds. A reason for the wide applicability of NMR lies in the development of Average-Hamiltonian-Theory[14] and subsequent decoupling pulse sequences[15, 16] as well as mechanical rotation methods[17, 18], leading to resonance line narrowing from formerly several hundred or thousands of ppm to often sub-ppm resolution[19]. While such high resolution techniques find widespread use in the fields of biochemistry, chemistry, or soft-matter physics, an application within diamond anvil cell research has

∗ thomas.meier@uni-bayreuth.de

Recent developments of high pressure NMR techniques, using three-dimensional microcoils[20] or twodimensional Lenz-lens based resonators[21-23], enabled the detection of faint nuclear induction signals at pressures well into the megabar (1 Mbar = 100 GPa) regime with record pressures of up to 202 GPa, as reported recently[24]. Most of these high pressure NMR experiments have been focused on 1 H nuclei[25], owing to the large gyromagnetic ratio of protons, and subsequently high signal-to-noise ratios[26]. At such pressures, a plethora of exotic physical phenomena can be observed, such as the elusive hydrogen-bond symmetrisation from ice VII to ice X governed by quantum-mechanical tunneling of protons within the H-bonds[22], or pressure induced hydrogen-hydrogen interactions in metal hydrides as a precursor of high temperature superconductivity[24]. Monserrat et al.[27] pointed out that NMR might be the only method available in diamond anvil cell based research to directly probe the atomic arrangement of hydrogen atoms at high densities. Using ab-initio calculations, theoretical 1 H-NMR spectra at megabar pressures have been calculated, showing a surprisingly large dispersion of chemical shielding values, yielding the possibility to distinguish between theoretical high pressure phases of dense hydrogen.However, an attempt to identify crystallographic phases of H 2 seemed unreasonable, as 1 H-NMR linewidths in excess of 10 3 ppm for ortho-hydrogen have been reported under cryogenic temperatures[28-31]. Under decreasing H-H distances, this effect would certainly be more pronounced in high density phases of H . 2

Here we present the first application of homo-nuclear LeeGoldburg decoupling[32] of molecular ortho-hydrogen in diamond anvil cells at pressures between 20 and 64 GPa.

Two DACs with 100 µ m and 250 µ m culeted diamonds have been prepared. Rhenium gaskets were pre-indented to about 15 and 25 µ m thickness and the central hole has been laser drilled into the centre of the flat indented area, forming the sample cavity. Subsequently, the diamonds have been coated with a layer of 1 µ m copper,

FIG. 1. Experimental high pressure NMR set-up. a) NMR probe with mounted DAC placed in the centre of a 1.4 T ESR magnet. b) view inside a DAC for NMR experiments. The Helmholtz coil has been made from 100 µ m PTFE insulated copper wire. A 1 µ m copper layer was deposited on the diamonds surface using physical vapor deposition. The rhenium gasket was electrically insulated using chemical vapor deposition of a 500 nm thick layer of Al 2 O . c) SEM photograph of a 100 3 µ m culet diamond with three Lenz lenses on the diamonds pavilion, the beveled region as well as the culet.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01150v1-figure-ef7b5103-04a9-4dac-a5b1-dfa0c8c336f0.jpg)

and the rhenium gaskets with a 500 nm layer of electrically insulating Al 2 O . Two dimensional Lenz lens NMR 3 resonators have been cut out of the copper layers using focused ion beam milling.

All excitation coils were manufactured from 100 µ m PTFE insulated copper wire, each consisting of 5 turns, and fixed around the diamond anvils. After loading, both coils were combined to form a Helmholtz-like excitation coil, providing a homogeneous inductive coupling into the respective Lenz lenses. Figure 1 shows three photographs of the set-up used: BX90[33] pressure cell mounted on a home-built NMR probe in the region of highest homogeneity of an ESR magnet (a); photograph of the inside of the DAC with both copper coated diamonds, the sealing rhenium gasket, as well as the Helmholtz excitation coil (b) and an SEM photograph of a prepared triple stage Lenz lens on the surface of a 100 µ m culeted diamond.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01150v1-figure-66feb323-212a-43f1-b740-8fa830a69e67.jpg)

Pressure cells were loaded with molecular hydrogen at the European Synchrotron Radiation Facility. Pressure was increased under cryogenic conditions by submerging the cells into liquid nitrogen, and Raman spectra were recorded at room temperature after every pressure increase, see figure 2. The pressure induced shift of the recorded vibron spectra at about 4250 cm -1 are in accordance with previous studies[34, 35].

All NMR experiments have been conducted at an ex-

FIG. 2. Raman spectra of molecular H 2 under increasing pressure.

ternal magnetic field of about 1063 mT, corresponding to a proton resonance frequency of 45.26 MHz. Chemical shift referencing has been performed, using a water sample as a reference standard in a BX90 pressure cell prepared in the same way as the hydrogen loaded cells. Performing 1 H-NMR nutation experiments at an average

pulse power of 8 W has shown optimal 90 ◦ pulse lengths of about 0.4 to 0.5 µ s, corresponding to RF fields of B 1 = π/ (2 γ n π/ t 2 ) = 11 to 15 mT. Supplementary figure 1 shows numerical simulations, using the FEMM software package, of the B 1 field distribution across a 40 µ m x 15 µ m sample cavity used when 100 µ m culeted anvils are used. An average value of 〈 B 1 〉 = 17 mT (green isolines in supplementary figure 1) was found, in excellent agreement with the nutation experiments. Considering a standard deviation of about σ = 3 mT (red isolines in supplementary figure 1), it was found that approximately 80 % of the sample volume lies within 〈 B 1 〉± σ . Thus, an effective volume V eff , which is subject to a homogeneous B 1 field of V eff = 0 8 . · V 0 , has been defined.

According to the equation of state of molecular hy-

FIG. 3. a) time domain data of the standard FID (blue) and of the FIDRF (orange) after optimisation of LG conditions. The yellow envelope is guide to the eye. B) corresponding 1 H-NMR spectra under standard static conditions (blue) and under LG conditions (yellow), the inset shows a zoomed-in region of the narrowed NMR signal.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01150v1-figure-3b3083b2-eb49-4c5f-8655-6fdd539b6d98.jpg)

drogen in phase I, the molar volume at 20 GPa is 5.25

cm 3 /mol [11]. Therefore, an upper limit of a total amount of 3 nmol of hydrogen loaded in the sample cavities could be estimated, corresponding to 1.7 · 10 15 H 2 molecules. Considering a normalised time domain signal-to-noise ratio of 8.8 and a limiting bandwidth of 1 MHz, a limit-ofdetection[21] of 2 10 · 11 spins/Hz 1 / 2 has been found.

Performing standard single pulse experiments, an intense and rapidly decaying 1 H signal was detected, at the respective Larmor frequency of protons at 1063 mT, figure 3a. Fourier transform results in a dipolar broadened resonance spectrum of about 5000 ppm line width, figure 3b. Using inversion recovery experiments, spin lattice relaxation times between 15 to 21 ms have been found, in accordance with other NMR studies of Phase I of H 2 at cryogenic temperatures[28, 30, 31].

Lee-Goldburg decoupling experiments have been conducted at varying offset frequencies in steps of 10 kHz around the centre of gravity of the broad hydrogen spectrum shown in figure 3b. Figure 4 demonstrates the dependence of the linewidth after a Fourier transform of the free induction decay in the rotating frame (FIDRF). As can be seen, a resonance line-width minimum is observed at about foff = 20 kHz. In these first experiments, the saturating Lee-Goldburg pulses were incremented from initially 1 s in steps of 500 ns for 8000 increments, resulting in a rotating frame time domain of 4 ms.

Figure 3a also shows the resulting FIDRF at 20 GPa.

FIG. 4. Dependence of the 1 H-NMR linewidth after FT of the free induction decay in the rotating frame for different off-resonance frequencies f off at 20 GPa.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01150v1-figure-212d1c55-0e42-46a4-9081-edd6f211c49a.jpg)

The free induction decay is stretched by a factor proportional to | λ 0 (Θ) | -1 = | 1 2 (3 cos 2 Θ - | 1 -1 where Θ denotes the angle between the direction of the irradiated LG-field and the external magnetic field B 0 acting on the high pressure chamber. Comparison of the 1 /e decay times, 900 ns for the FID (blue line in figure 3a) and about 2 ms for the FIDRF (orange trace in figure 3a) suggests an elongation of the FID by a factor of about 2000 and thus | 1 2 (3 cos 2 Θ - | 1 -1 ≈ 4 5 . · 10 -4 or Θ = 54 717 , close to . ◦ the magic angle.

As the FIDRF is not fully decayed within the predefined time frame of 4 ms, the minimal line-width of its Fourier transform can be estimated to be FWHM ≈ 1 / (4 ms ) = 250 Hz ≡ 5 5 . ppm . Experimentally, we find a line width of about 6.4 ppm. By increasing the time scale of the

LG-pulse to 16 ms (8000 increments in 2 µ s steps), the FIDRF has been found to be almost completely decayed within the time frame of the LG-pulse, resulting in a spectral resolution of 3.1 ppm (see yellow spectrum in figure 3b).

Figure 5 shows a summary of the FWHM line widths of 1 H-NMR spectra of molecular ortho-hydrogen in a pressure range of 20 to 64 GPa, after standard single pulse excitation and LG-decoupling. In almost all cases, line widths of about 3 to 5 ppm could be realised.

All resonance spectra were found to be almost fully represented by a Lorenzian line shape without any detectable signs of chemical shift anisotropy, evidencing the suitability of this approach for a purely dipolar broadened spin1/2 systems. An interpretation of the chemical shift dependencies at increasing densities will be published elsewhere.

The aim of this study was to elucidate the possibility to achieve high spectral resolutions in diamond anvil cell based NMR. It was found that even one of the most basic homo-nuclear decoupling experiments, i.e. Lee-Goldburg irradiation, can be utilised to achieve resolutions as low as 3 ppm. While this value might still be considered too high when compared to magic angle spinning techniques, it represents a major improvement in spectral resolution in diamond anvil cells at extreme conditions.

The found resolution of 3 ppm may be improved in future experiments by extending the time domain in the rotating frame, or by achieving a higher SNR per increment by using a higher number of accumulations. The appli-

- [1] E. Wigner and H. B. Huntington, 'On the Possibility of a Metallic Modification of Hydrogen', J. Chem. Phys. , 3 , 764-770 (1935).

[2] N W Ashcroft, 'Metallic hydrogen: A high-temperature superconductor?', Physical Review Letters , 21 , 17481749 (1968).

[3] H.-k. Mao and R. J. Hemley, 'Ultrahigh-pressure transitions in solid hydrogen', Reviews of Modern Physics , 66 , 671-692 (1994).

[4] R. Dias and I. F. Silvera, 'Observation of the WignerHuntington Transition to Solid Metallic Hydrogen', Science , 355 , 715-718 (2017).

- [5] M. I. Eremets and A. P. Drozdov, 'Comments on the claimed observation of the Wigner-Huntington Transition to Metallic Hydrogen', ArXiv:1702.05125 , (2017).
- [6] X.-D. Liu, P. Dalladay-Simpson, R. T. Howie, B. Li, and E. Gregoryanz, 'Comment on Observation of the Wigner-Huntington transition to metallic hydrogen', Science , 357 , eaan2286 (2017).
- [7] A. F. Goncharov and V. Struzhkin, 'Comment on Observation of the Wigner-Huntington transition to metallic hydrogen, Science , 357 , eaam9736 (2017).

FIG. 5. Comparison of FWHM line widths after standard 90 pulse excitation, and after Lee-Goldburg decoupling in a pressure range from 20 to 64 GPa.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01150v1-figure-65803519-ef95-473e-a984-60ac93cba48f.jpg)

cability of homo-nuclear Lee-Goldburg decoupling is particularly powerful in systems governed by direct dipoledipole interactions, such as H , which are almost com2 pletely averaged out during the LG-experiments leaving only isotropic chemical shift interactions present. These results open the way for experiments at even higher pressures, and may eventually resolve the predicted structural models of Phases II and III of hydrogen[36].

## I. ACKNOWLEDGEMENTS

The authors thank the German Research Foundation (Deutsche Forschungsgemeinschaft, DFG, projects DU 954/11-1, DU 393/13-1, DU 393/9-2, and ME 5206/31) and the Federal Ministry of Education and Research, Germany (BMBF, grant no. 05K19WC1) for financial support.

- [8] P. Loubeyre, F. Occelli, and P. Dumas, 'Observation of a first order phase transition to metal hydrogen near 425 GPa', arxiv:1906.05634 ,(2019).
- [9] J. Van Kranendonk, 'Solid Hydrogen' . Plenum Press, New York, NY, first edition, (1983).
- [10] N. Stojilovic, 'Why Can't We See Hydrogen in X-ray Photoelectron Spectroscopy?', J. of Chem. Edu. , 89 , (10), 1331-1332 (2012).
- [11] P. Loubeyre, R. LeToullec, D. Hausermann, M. Hanfland, R. J. Hemley, H. K. Mao, and L. W. Finger, 'X-ray diffraction and equation of state of hydrogen at megabar pressures', Nature , textbf83, 702-704 (1996).
- [12] R. Boehler, M. Guthrie, J. J. Molaison, A. M. dos Santos, S. Sinogeikin, S. Machida, N. Pradhan, and C. A. Tulk, 'Large-volume diamond cells for neutron diffraction above 90 GPa', High Pressure Research , 33 , 546554 (2013).
- [13] I. I. Mazin, Russell J. Hemley, A. F. Goncharov, M. Hanfland, and H. K. Mao, 'Quantum and classical orientational ordering in solid hydrogen', Physical Review Letters , 78 , 1066-1069 (1997).

- [14] A. Brinkmann, 'Introduction to average Hamiltonian theory. I. Basics', Concepts in Magnetic Resonance Part A: Bridging Education and Research , 45A , (2016).
- [15] M. Mehring and J. S. Waugh, 'Magic-angle NMR experiments in solids', Physical Review B , 5 , 3459-3471 (1972).
- [16] E. Vinogradov, P. K. Madhu, and S. Vega, 'Phase modulated Lee-Goldburg magic angle spinning proton nuclear magnetic resonance experiments in the solid state: A bimodal Floquet theoretical treatment'. Journal of Chemical Physics , 115 , 8983-9000 (2001).
- [17] I. J. Lowe, 'Free Induction Decays of Rotating Solids', Physical Review Letters , 2 , 285-287 (1959).
- [18] E. R. Andrew, A. Bradbury, and R. G. Eades, 'Nuclear Magnetic Resonance Spectra from a Crystal rotated at High Speed'. Nature , 182 , 1659-1659 (1958).
- [19] J. W. Hennel and J. Klinowski, 'Magic-Angle Spinning: a Historical Perspective', In Klinowski J. (eds), editor, New Techniques in Solid-State NMR. Topics in Current Chemistry , volume 2, pages 1-14. Springer, Berlin, Heidelberg, 246 edition, (2005).
- [20] T. Meier, 'At Its Extremes: NMR at Giga-Pascal Pressures', Graham Webb (eds.), Annual Reports on NMR Spectroscopy , chapter 1, pages 1-74. Elsevier, London, 93th edition, (2018).
- [21] T. Meier, N. Wang, D. Mager, J. G. Korvink, S. Petitgirard, and L. Dubrovinsky, 'Magnetic flux tailoring through Lenz lenses for ultrasmall samples: A new pathway to high-pressure nuclear magnetic resonance', Science Advances , 3 ,eaao5242 (2017).
- [22] T. Meier, S. Petitgirard, S. Khandarkhaeva, and L. Dubrovinsky, 'Observation of nuclear quantum effects and hydrogen bond symmetrisation in high pressure ice'. Nature Communications , 9 , 2766(2018).
- [23] T. Meier, S. Khandarkhaeva, S. Petitgirard, T. K¨rber, o A. Lauerer, E. R¨ssler, and Leonid Dubrovinsky, 'NMR o at pressures up to 90 GPa', Journal of Magnetic Resonance , 292 ,44-47 (2018).
- [24] T. Meier, F. Trybel, S. Khandarkhaeva, G. SteinleNeumann, S. Chariton, T. Fedotenko, S. Petitgirard, M. Hanfland, K. Glazyrin, N. Dubrovinskaia, and L. Dubrovinsky, 'Pressure-Induced Hydrogen-Hydrogen Interaction in Metallic FeH Revealed by NMR'. Physical Review X , 9 , 031008 (2019).
- [25] T. Meier. 'Journey to the centre of the Earth: Jules Vernes' dream in the laboratory from an NMR perspec-
- tive', Progress in Nuclear Magnetic Resonance Spectroscopy , 106-107 , 26-36 (2018).
- [26] M..H Levitt. Spin Dynamics: Basics of Nuclear Magnetic Resonance . (2000), ISBN 9780470511183.
- [27] B. Monserrat, S. E. Ashbrook, and C. J. Pickard, 'Nuclear Magnetic Resonance Spectroscopy as a Dynamical Structural Probe of Hydrogen under High Pressure', Physical Review Letters , 122 , 135501 (2019).
- [28] W. N. Hardy and A. J. Berlinsky, 'Anomalous NMR Line Shapes in Solid Ordered Ortho-Hydrogen and ParaDeuterium', Physical Review B , 8 , 4996-5012, (1973).
- [29] A. J. Berlinsky and W. N. Hardy, 'Theory of OrthoPara Conversion and its Effect on the NMR Spectrum of Ordered Solid Ortho-Hydrogen', Physical Review B , 8 , 5013-5027 (1973).
- [30] S. Washburn, M. Calkins, H. Meyer, and A. B. Harris. 'Pulsed NMR studies in solid H2. III. Spin-lattice relaxation times', Journal of Low Temperature Physics , 53 , 585-617 (1983).
- [31] P. L. Pedroni, R. Schweizer, and H. Meyer, 'NMR relaxation times and line shapes in solid H2 at elevated densities', Physical Review B , 14 , 896-905 (1976).
- [32] M. Lee and W. I. Goldburg, 'Nuclear-MagneticResonance Line Narrowing by a Rotating rf Field', Physical Review , 140 , A1261-A1271 (19650.
- [33] I. Kantor, V. Prakapenka, A. Kantor, P. Dera, A. Kurnosov, S. Sinogeikin, N. Dubrovinskaia, and L. Dubrovinsky, 'BX90: A new diamond anvil cell design for X-ray diffraction and optical measurements', Review of Scientific Instruments , 83 ,125102 (2012).
- [34] S. K. Sharma, H. K. Mao, and P. M. Bell, 'Raman Measurements of Hydrogen in the Pressure Range 0.2-630 kbar at Room Temperature', Physical Review Letters , 44 , 886-888 (1980).
- [35] R. J. Hemley and H. K. Mao, 'Phase transition in solid molecular hydrogen at ultrahigh pressures', Physical Review Letters , 61 , 857-860 (1988).
- [36] Y. A. Freiman and Y. Crespo, 'Isotopic and spin-nuclear effects in solid hydrogens', Low Temperature Physics , 43 , 1345-1361 (2017).

## II. SUPPLEMENTARY

FIG. 6. Finite element simulation of the RF B 1 field generated by the 3rd stage Lenz lenses across a 12 pl sample cavity. The green lines denote an average value of about 〈 B 1 〉 = 17 mT . The red lines denote the upper and lower (defined by the standard deviation) boundaries of the effective volume V eff .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01150v1-figure-1c065b0a-7a3f-45bc-af4f-534e01986a40.jpg)
Thomas Meier Saiana Khandarkhaeva Jeroen Jacobs Natalia Dubrovinskaia Leonid Dubrovinksy
Improving Resolution of Solid State NMR in Dense Molecular Hydrogen
Recent advancements in radio-frequency resonator designs have led to the
implementation of nuclear magnetic resonance in diamond anvil cells (DACs) at
pressures well above 100 GPa. However, a relatively low resolution and the
absence of decoupling sequences complicate the analysis of the results of solid
state NMR in DACs. Here, we present the first application of homo-nuclear
Lee-Goldburg (LG) decoupling at extreme conditions on high density molecular
hydrogen. Lenz lens based two-dimensional resonator structures were found to
generate a homogeneous B1 field across sample cavities as small as 12
picolitres, a prerequisite for optimal decoupling. At ideal LG conditions, the
broad 1H resonance of molecular ortho-hydrogen was narrowed 1600-fold,
resulting in line-widths of 3.1 ppm.