# Fast methane diffusion at the interface of two clathrate structures

Umbertoluca Ranieri, Michael Marek Koza, Werner F Kuhs, Stefan Klotz, Andrzej Falenty, Philippe Gillet, Livia E Bove

► To cite this version:

Umbertoluca Ranieri, Michael Marek Koza, Werner F Kuhs, Stefan Klotz, Andrzej Falenty, et al. Fast methane diffusion at the interface of two clathrate structures. Nature Communications, 2017, 8, pp.1076. 10.1038/s41467-017-01167-2. hal-01629537

HAL Id: hal-01629537

https://hal.sorbonne-universite.fr/hal-01629537

Submitted on 6 Nov 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ARTICLE

DOI: 10.1038/s41467-017-01167-2

OPEN

# Fast methane diffusion at the interface of two clathrate structures

Umbertoluca Ranieri ID <fnref n="1" />,², Michael Marek Koza ID ², Werner F. Kuhs ID ³, Stefan Klotz⁴, Andrzej Falenty ID ³, Philippe Gillet<fnref n="1" /> & Livia E. Bove ID <fnref n="1" />,⁴

Methane hydrates naturally form on Earth and in the interiors of some icy bodies of the Universe, and are also expected to play a paramount role in future energy and environmental technologies. Here we report experimental observation of an extremely fast methane diffusion at the interface of the two most common clathrate hydrate structures, namely clathrate structures I and II. Methane translational diffusion—measured by quasielastic neutron scattering at 0.8 GPa—is faster than that expected in pure supercritical methane at comparable pressure and temperature. This phenomenon could be an effect of strong confinement or of methane aggregation in the form of micro-nanobubbles at the interface of the two structures. Our results could have implications for understanding the replacement kinetics during sl-sll conversion in gas exchange experiments and for establishing the methane mobility in methane hydrates embedded in the cryosphere of large icy bodies in the Universe.

<footnote n="1">¹EPSL, ICMP, École polytechnique fédérale de Lausanne (EPFL), Station 3, CH-1015 Lausanne, Switzerland. ²Institut Laue-Langevin, 71 avenue des Martyrs, CS 20156, 38042 Grenoble cedex 9, France. ³GZG Abt. Kristallographie, Universität Göttingen, Goldschmidtstrasse 1, 37077 Göttingen, Germany. ⁴Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Université Pierre et Marie Curie Paris 06, CNRS Unité Mixte de Recherche 7590, Sorbonne Universités, F-75252 Paris, France. Correspondence and requests for materials should be addressed to U.R. (email: ranieriu@ill.fr) or to L.E.B. (email: livia.bove@impmc.upmc.fr)</footnote>
Gas clathrate hydrates are a general class of compounds composed of water (hosts) molecules and gas (guests) atoms or molecules¹. Among them, clathrate hydrates of methane are probably the most extensively studied. The natural occurrence of methane hydrate-containing sediments in permafrost areas and ocean shelves, the risk due to their formation in oil and gas pipelines, as well as their potential application as gas transportation media in soft conditions (i.e., close to atmospheric pressure and room temperature) explain the wide interest shown for these materials¹, ². Exchanging the guests in natural gas hydrate deposits with CO₂ has been suggested as a two-in-one approach of energy recovery and concomitant CO₂ mitigation³. As they are believed to be the dominant methane-bearing phase in the nebula from which the outer planets and satellites are formed, the properties of methane hydrates are also crucial to models of bodies in the outer solar system⁴. From a physical–chemical point of view, methane hydrates represent prototypical examples of hydrates of hydrophobic guests: the combination of low temperature, high pressure, a weak gas–water interaction between the guest molecules and the host lattice, and the relatively strong hydrogen bonds between host molecules allow for the formation of open crystalline water networks encaging CH₄ molecules. The topology of the water cages and the number of gas molecules trapped in these cages critically depend on the specific thermodynamic conditions of formation of the clathrate hydrate and on its formation kinetics¹, ⁵.

The most common structures formed by clathrate hydrates at relatively moderate pressures are the clathrate structures I and II (noted sI and sII). The unit cell of clathrate sI (space group Pm$\overline{3}$n) contains two small dodecahedral (5$^{12}$) water cages and six bigger tetrakaidecahedral (5$^{12}$6$^{2}$) cages. The unit cell of sII (space group Fd$\overline{3}$m) contains 16 5$^{12}$ cages and eight large hexadecahedral (5$^{12}$6$^{4}$) cages¹ (Fig. 1). It is well accepted that methane hydrates preferentially crystallise into sI¹. However, cages characteristic of sII have been transiently detected in the initial stages of the formation of methane hydrates in both experiments⁶, ⁷ and simulations⁸⁻¹². This is not surprising since (i) the difference in free energy between sI and sII is small¹³ and (ii) appearance of metastable polymorphs or transient non-equilibrium states is commonly observed during nucleation of hydrates⁵, ¹³⁻¹⁶. It is noteworthy that sI and sII are topologically incompatible without

the intercalation of pentakaidecahedral (5$^{12}$6$^{3}$) cages;⁸ the interplay between kinetic factors and thermodynamic stability during sI–sII cross-nucleation has been discussed in details¹⁷. In methane hydrates at room temperature and pressures up to 0.6 GPa, sII has been reported to persistently coexist with sI$^{18–20}$. Therefore, the resulting coexistence of structures in high-pressure samples can be seen as a frozen form on laboratory timescales of the metastable sI–sII polymorphs usually encountered during nucleation of methane hydrates.

Low-temperature translational and rotational excitations, as well as cage-to-cage hopping of CH₄ molecules trapped in clathrate sI were previously investigated at ambient and low pressures$^{21–24}$. However, no information is available on the extra-cage diffusivity of the guest molecules in methane hydrates; this information could be highly relevant for the modelling of the subcrustal layers of methane clathrates embedded in the cryosphere of icy planets and large icy satellites$^{25, 26}$. Recently, a study based on molecular dynamics simulations reported diffusion coefficient values in the nanosecond time scale for methane diffusion at grain boundary-like structures of defective clathrates$^{27}$.

In this work, we probe the microscopic diffusion of methane in a methane hydrate (CH₄–D₂O) sample exhibiting coexistence of clathrate sI and sII by quasielastic neutron scattering (QENS) measurements. Coexistence of structures is promoted by applying high pressure. QENS is a well-suited technique to study dynamics on the picosecond time and Å length scales$^{28}$. Spectra of the sI–sII clathrate show a clear quasielastic signal whose analysis reveals a very fast extra-cage translational diffusion of methane molecules on the picosecond time scale. For comparison, we also measure methane hydrates in pure sI clathrate, in pure hexagonal clathrate structure H (space group P6/mmm)$^{26}$ and during transformation from sI to structure H (noted sH); the spectra of sI and sH do not exhibit any visible quasielastic signal, and the spectra of sI–sH show a very weak signal, orders of magnitude smaller than the signal from sI–sII.

## Results

**QENS experiments and elastic scattering.** The experiments were performed at the time-of-flight spectrometer IN6 at the Institut Laue-Langevin in Grenoble (France) using a Paris-Edinburgh

<figure><img src="image_2.png" /><figcaption>**Fig. 1** Neutron diffraction patterns. Powder diffraction patterns of methane hydrate in pure sI clathrate at 0.4 GPa and 290 K and in the sI–sII clathrate at 0.8 GPa and 282 K. Breaks correspond to noisy detectors and to the strong Bragg peak of alumina from the anvils at 95°. The positions of the Bragg peaks for sI (cell parameter 11.7 Å) and for sII (cell parameter 17.0 Å) are reported. On the right, we present views of the unit cells of sI and sII (5$^{12}$ cages in cyan, 5$^{12}$6$^{2}$ cages in purple, 5$^{12}$6$^{4}$ cages in blue)</figcaption></figure>
<figure><img src="image_3.png" /><figcaption>**Fig. 2** Examples of measured QENS spectra. QENS spectra of methane hydrate in the sI-sII clathrate at 0.8 GPa and selected temperature *T* and momentum transfer *Q* values. Experimental data (empty circles) are compared to their best fits (black lines). Error bars were calculated by the square root of absolute neutron count combined with the law of propagation of errors. Quasielastic Lorentzian (solid green lines) and elastic (dashed green lines) components are also shown (upshifted by the value of the flat background for clarity). In the inset, a Lorentzian fit is compared to the 2D diffusion fit (blue line) of the same spectrum</figcaption></figure>

press and recently developed ceramics anvils²⁹. The wavelength of the incoming neutrons was 5.12 Å, resulting in an instrumental resolution of 0.08 meV. This corresponds to an observation time of ~8 ps. The sample exhibiting coexistence of clathrate sI and sII was prepared according to the following procedure: we compressed methane hydrate (originally in sI) to 0.8 GPa at liquid nitrogen temperature and then warmed it up to 282 K. The neutron powder diffraction pattern of the sample at 0.8 GPa and

282 K is presented in red in Fig. 1. It indicates that the sample contained about half as much sII than sI, on the basis of peak heights. The pattern was obtained directly on IN6 by comparing the intensities of the elastic peaks, at each scattering angle, with those measured on a vanadium standard which gives isotropic elastic scattering. All the Bragg peaks of the sample can be indexed within the space groups of sI and sII. Figure 1 also depicts the diffraction pattern of pure sI clathrate. The diffraction
patterns of pure sH clathrate and of the sample transforming from sI to sH are presented in Supplementary Fig. 1. We recorded QENS spectra of the sI-sII clathrate at the constant temperature T of 282 K and pressure of 0.8 GPa during 6 h. The amount of sII was constant during this time. Then we continuously decreased the temperature to 200 K over 15 h to characterize the T dependence of the probed diffusion at 0.8 GPa. Spectra measured between 282 and 200 K were merged into three groups of 5 h of acquisition time each, corresponding to the following average temperature values: 267, 238 and 212 K. Upon cooling down, the relative amount of sI and sII remained approximately constant. The diffraction patterns recorded at 267, 238 and 212 K are reported in Supplementary Fig. 2. More details about the experiments are given in the Methods section.

**Fitting of the QENS spectra.** Figure 2 depicts typical QENS spectra of the methane hydrate sample exhibiting coexistence of clathrate sI and sII at 0.8 GPa. Examples of spectra of methane hydrate in pure sI clathrate, in pure sH clathrate and in the sI-sH clathrate are shown in Supplementary Fig. 3. Spectra in Fig. 2 show a clear quasielastic signal, i.e., a broadening of the elastic line produced by interactions of the neutrons with diffusing atoms of the sample. Since the incoherent cross-section of hydrogen is almost two orders of magnitude larger than that of other atoms, the measured signal is essentially due to the dynamics of protons in the guest molecules³⁰, ³¹. We first applied the most common model used to fit quasielastic data, i.e., a Lorentzian function (whose half-width-half-maximum is noted Γ). A delta function was used to fit the elastic line of the spectra, which is due to the contribution of the water network and of the slowly-diffusing or non-diffusing methane molecules trapped in the clathrate structures. Total best fits to the experimental data are presented in Fig. 2 and can be seen to accurately describe the spectra. Based on the integrated areas of the quasielastic and elastic lines (after subtraction of the empty cell measurement), we roughly estimate that about one third of the methane molecules in the sample contribute to the fitted quasielastic signal, at each investigated T. More details on this estimation are given in Supplementary Note 1. Since cage occupancies in the newly formed sII clathrate might be lower compared to the starting sI clathrate hydrate, part of the methane molecules in the starting sI clathrate could indeed have been released from the starting sI hydrate into the grain boundary network during transformation from sI to sII and would be available to perform extra-cage translational diffusion. However, a minimum level of occupancies is required to ensure stability of sII and one can estimate that no more than 10% of the methane in the sample could have been released without destabilization of the water matrix. The existence of a fraction of fast diffusing methane molecules as high as one third strongly suggests that an appreciable fraction of water molecules in the sample are in a disordered state. Such disordered regions would form at the front line of the transformation between clathrate sI and sII, and their sizes are most likely far below the typical size of the crystallites (that is a few micrometres³²). This point is further discussed in Supplementary Note 2. Moreover, the absence of a prominent quasielastic signal in the spectra of the sI-sH methane clathrate hydrate highlights the very particular nature of the interfaces between coexisting sI and sII, compared to the temperature-induced or pressure-induced structural transition taking place at high driving forces between two stable forms of methane hydrates such as sI and sH. The micro-structural properties of sI and sII coexisting assemblies certainly deserve to be further investigated.

**Momentum transfer Q dependence of the QENS signal.** Figure 3 depicts the parameter Γ as a function of Q². The Q dependence of

<figure><img src="image_4.png" /><figcaption>**Fig. 3** Momentum transfer Q dependence of the QENS signal. Half-width-half-maximum Γ of the Lorentzian quasielastic component of the fits (Fig. 2) as a function of Q² at 0.8 GPa and the investigated temperatures. Error bars correspond to one standard deviation. The best fits to the data according to Eq. (1) are shown as dashed lines</figcaption></figure>

Γ provides information about the characteristic time and nature of the probed motion. The monotonic increase of Γ rules out that the measured quasielastic signal is due to a localised (rotational) dynamics of methane, which would be indicated by a Q-independent Γ. Instead, it clearly highlights that a translational diffusion process is at the origin of the signal²⁸. It must be also noted that the rotational quasielastic contribution of methane molecules trapped in clathrate sI is very large (half width above 5 meV) at 150 K²² and thus only contributes to the background of the spectra here. Similar rotational behaviour can be reasonably expected for CH₄ molecules in a clathrate sII, as no indication of inequivalent environments for the guest molecule emerged from the low-temperature rotational spectra of sI methane clathrate²³ (although methane occupies the two types of cages of sI). As can be seen in Fig. 3, Γ extrapolates to 0 for Q → 0. Hence, the measured quasielastic signal is not associated with an intra-cage diffusive motion of CH₄ molecules, since for a particle restricted to move in a confined geometry Γ shows³⁰, ³¹ a plateau at small Q. For example, for a particle moving within a sphere of radius R, Γ shows³³ a plateau for Q < π/R. The Q dependence of Γ is best approximated within the random jump diffusion model of Singwi and Sjolander by:

$$ \Gamma(Q) = \frac{\hbar DQ^2}{1 + DQ^2\tau}, \qquad (1) $$

with D representing the isotropic translational diffusion coefficient and τ the mean residence time between jumps²⁸. The corresponding formula for a continuous free translational diffusion process would be Γ(Q) = ħDQ². Fits of Γ(Q²) according to Eq. (1) are presented in Fig. 3; the values obtained for D and τ are reported in green in Fig. 4. The translational diffusion coefficient turns out to be of the order of 10⁻⁴ cm² s⁻¹ and its temperature dependence is rather weak (25% over the investigated T range). An Arrhenius fit of D provides an activation energy of 0.48 ± 0.11 kcal mol⁻¹. This value is small compared to the activation energies reported in literature for the cage-to-cage hopping of CH₄ in sI clathrates (for example, 12.4 kcal mol⁻¹ in ref. ²¹) and points at van der Waals interactions as main rate-limiting interactions for the observed methane diffusion. The parameter τ is a fraction of picosecond and does not show any temperature dependence within the error bars over the investigated T range.

**2D diffusion model.** The choice of a Lorentzian fit function for the quasielastic signal implicitly assumes that the probed motion is three-dimensional (3D)²⁸. An other possibility is that the
<figure><img src="image_5.png" /><figcaption>**Fig. 4** Translational diffusion coefficient *D* and mean residence time *τ*. Temperature dependencies of *D* (a) and *τ* (b) for methane hydrate in the sI-sII clathrate at 0.8 GPa, as obtained in the 3D and 2D diffusion models employed in this work. Error bars correspond to one standard deviation. Arrhenius fits of *D* are shown as dashed lines. Estimated value for *D* in pure methane from literature<sup>40, 41</sup> is also reported. Legend of **a** also applies to **b**</figcaption></figure>

methane diffuses essentially bi-dimensionally on the length scale probed by the instrument, if the grain boundary network or the intercalated disordered regions between crystals of sI and sII are sufficiently thin. In such a case the fit function for the quasielastic signal is no longer Lorentzian and has a logarithmic singularity at $\omega=0$<sup>34</sup> (see Supplementary Note 3 for its expression). Nevertheless, the singularity is suppressed by the convolution with the instrumental resolution and the convoluted fit function differs from the convoluted Lorentzian only near $\omega=0$ where it is more peaked<sup>35</sup>. The inset of Fig. 2 depicts an example of fit using this 2D diffusion fit function and compares it to the Lorentzian fit of the same QENS spectrum. The two fits are actually indistinguishable outside the instrumental resolution-dominated region close to $\omega=0$ and this is true for all other measured spectra. Therefore it is not possible to unequivocally establish if bulk or planar diffusion is taking place. The values for the translational diffusion coefficient and the mean residence time obtained in the 2D diffusion model are reported in blue in Fig. 4. The activation energy (0.57 ± 0.12 kcal mol<sup>-1</sup>) is comparable to that obtained in the 3D diffusion model. More details about the data analysis are given in the Methods section.

## Discussion

The methane diffusion probed in the present study is much faster than that reported in the literature for the cage-to-cage hopping of CH₄ molecules through clathrate sI. Cage-to-cage hopping is a rare event that requires distortion of the host network<sup>36</sup> and the associated diffusion coefficient is of the order of 10<sup>-11</sup> to 10<sup>-12</sup> cm² s<sup>-1</sup> at 250 K, as revealed by experimental<sup>21</sup> and computational<sup>24</sup> studies. Similar conclusions have been reported for the cage-to-cage hopping of other guest molecules<sup>37, 38</sup>, including molecules forming clathrate sII<sup>37</sup>. It is also interesting to compare the present results to the translational diffusion coefficients of CH₄ in bulk water-methane mixtures and bulk pure methane. At 0.02 GPa and 273 K, the diffusion coefficient of methane in water was found to be 0.16 × 10<sup>-4</sup> cm² s<sup>-1</sup>, i.e., an order of magnitude smaller than those measured here<sup>39</sup>. This value was obtained for the moderate methane-saturated concentration that is possible at low pressures<sup>39</sup>. In pure methane the experimental diffusion coefficient is 2.08 × 10<sup>-4</sup> cm² s<sup>-1</sup> at 0.164 GPa and 298 K<sup>40</sup>. Its temperature dependence is rather strong, with an activation energy of ~1.0 kcal mol<sup>-1</sup> between 223 and 323 K at 0.15 GPa. The diffusion coefficient at 0.8 GPa can be estimated based on the assumption that its product with the shear viscosity is constant along isotherms (Stokes-Einstein relation). The pressure dependence of the viscosity in methane at 298 K is known<sup>41</sup> and one gets a value of 0.5 × 10<sup>-4</sup> cm² s<sup>-1</sup> at 0.8 GPa and 298 K. This value is reported in

Fig. 4 and is a factor of 2–3 smaller than our results extrapolated at the same *T*. Based on the same assumption, it is possible to estimate that pure methane at about 0.2–0.3 GPa should show a diffusion coefficient comparable to that measured here.

To summarise, we observed a remarkably fast mobility of methane molecules at the interface of two clathrate structures (I and II) and measured the associated translational diffusion coefficient *D* at 0.8 GPa and temperatures between 212 and 282 K. The obtained coefficients are 7–8 orders of magnitude higher than those reported in literature for cage-to-cage hopping of methane through clathrate sI at low pressure, one order of magnitude higher than that of methane dissolved in water at low pressure and a factor of 2–3 higher than that expected for pure bulk supercritical methane at comparable pressure and temperature. The activation energy (of about 0.5 kcal mol<sup>-1</sup>) is a factor of two smaller than that of pure methane at 0.15 GPa and more than one order of magnitude smaller than that of the hydrogen bond in the water network and of the cage-to-cage hopping process as reported in literature<sup>21</sup>. This fast mobility involves a sizable fraction of the methane in the sample (approximately one third, as rough estimation), does not induce destabilization of the clathrate structures and is probably observable for times much longer than the duration of our experiment (~21 h).

We infer that the rapidity of the methane diffusion probed here could be an effect of confinement in the extensive grain boundary network<sup>32</sup> or intercalated disordered regions between crystals of clathrate sI and sII. Similar behaviour was already reported in literature. For example, the diffusivity of CH₄ is only 4 × 10<sup>-11</sup> cm² s<sup>-1</sup> in zeolite 4A<sup>42</sup>, ~10<sup>-4</sup> cm² s<sup>-1</sup> in silicalite<sup>43</sup> and is predicted to be of the same order of magnitude as that of the gas phase (10<sup>-1</sup> cm² s<sup>-1</sup>) in infinitely long single-walled carbon nanotubes<sup>44</sup>. Alternatively, the observed fast diffusion could also well be explained by the spontaneous formation of micro-scale to nano-scale gas bubbles from a supersaturated water-methane mixture. Micro-nanobubbles formation was suggested to occur after decomposition of hydrates in different works<sup>45-47</sup>. The diffusion properties of methane inside these bubbles can be considerably different from the bulk fluid and a first study of CH₄-mobility in nanobubbles suggested indeed an enhanced diffusion<sup>48</sup>. Further investigation including large-scale molecular dynamics simulations of the guest diffusivity at the structures interface are needed to shed light on these points.

In the context of energy recovery from natural gas hydrate deposits with CO₂ injections, gas replacement rates are key parameters to assess its technological viability. Earlier experimental evidences underlined greatly enhanced replacement rates during sI-sII conversion<sup>49, 50</sup> in comparison to the case of isostructural sI-sI replacement<sup>51</sup>. If extended to moderate pressures, our results might provide an explanation for that. Likewise, our results should be taken into account in the modelling of methane clathrates layers existing at depth in the interiors of large icy bodies in both solar and extra-solar systems<sup>25, 26</sup>, for which the steady-states depend on the diffusion timescales as compared to the formation and dissociation rates. As an example, the observed fast mobility of methane could be relevant to understand the phenomenon of methane release into the atmosphere of Titan, which is likely to originate from methane clathrates embedded in its crust and mantle<sup>25, 26</sup>.

## Methods

**Sample production.** The procedure followed to prepare the CH₄-D₂O methane hydrate sample was described in refs. 21, 52. It basically consists in keeping D₂O ice under an atmosphere of 6 MPa of CH₄ gas at a temperature close to the melting during 4 weeks. The starting deuterated ice was a powder of ice Ih of spherical shape (typical diameter of several tens of micrometres<sup>21</sup>) previously produced by a shock-freezing method through spraying liquid D₂O (99.9% deuterated) into liquid nitrogen. The spraying was done in a glove box under dry nitrogen atmosphere to
avoid contamination with atmospheric water. The quality of the prepared methane hydrate sample was checked by X-ray diffraction. We found that the sample was in clathrate SI with a negligible amount of water ice impurity (below 2%). The size of the crystallites is typically a few micrometres³². Typical methane occupation is 86% for the small cages and 99% for the large cages³².

**Experimental details.** The QENS experiments were carried out using the VX5 Paris-Edinburgh press. The procedure of loading the methane hydrate sample in the clamp module⁵³ of the press was done under liquid nitrogen. The sample was first compacted to a spherical pellet (of ~40 mm³) using a dedicated press operating under liquid nitrogen. The pellet was subsequently loaded into a precooled type-25 copper-beryllium encapsulating gasket and the sample-gasket assembly was placed in an aluminium ring between precooled ceramics anvils. We used recently developed zirconia-toughened alumina ceramics anvils which are highly transparent to neutrons. Their performances are described in ref. ²⁹.

To prepare the sample exhibiting coexistence of clathrate SI and SII, the gasket was sealed by applying a load of 100 kN on the anvils under liquid nitrogen. This corresponds to a pressure of about (0.8 ± 0.1) GPa in the sample, on the basis of our calibration of the used anvils. The assembled clamp was then warmed up from liquid nitrogen temperature to room temperature out of the beam before insertion (~12 h later) in the Paris-Edinburgh press. During the experiment, temperature was decreased by cooling down the whole Paris-Edinburgh press in a liquid nitrogen cryostat. It is known that the cooling of samples in such a pressure cell is approximately isochoric and this leads to a small pressure drop (typically below 5% for a change in temperature between 282 and 200 K). The measured Bragg peaks did not shift with temperature within the angular resolution of the instrument.

During a different sample loading, the gasket was sealed by applying a smaller load (50 kN), corresponding to a sample pressure of 0.4 GPa. After being warmed up to 290 K, this sample was still in pure structure I (see Fig. 1 for the diffraction pattern). We compressed this sample isothermally at 290 K and observed transformation to clathrate sH, in agreement with previous studies²⁰, ²⁶.

During another sample loading, the gasket was sealed by applying a higher load (120 kN), corresponding to a sample pressure of 1.0 GPa. After being warmed up to 295 K, this sample was found to contain a mixture of structure I and structure H. The relative amount of structure H was found to slowly increase over time and the transformation was completed within ~12 h.

The instrumental energy resolution was estimated by measuring a sphere of vanadium of the same size as the sample, which was loaded into the gasket and the Paris-Edinburgh press in the same set-up as the sample at ambient pressure and ambient temperature.

**Data analysis details.** The scattering angles $2\theta$ covered by the detectors of IN6 are in the range $10^{\circ}-115^{\circ}$. Spectra measured by several detectors were grouped together into constant-Q spectra with $0.2 \AA^{-1}$ steps, from 0.4 to $1.8 \AA^{-1}$. For the data analysis, we did not consider the two highest Q values (1.6 and $1.8 \AA^{-1}$) for which competition of the quasielastic signal with the flat background gives rise to large uncertainties for the free-fitting parameters. Six free-fitting parameters were used in the data fitting with the 3D diffusion model: intensities and half-width-half-maximum of the Lorentzian and delta functions, flat background and zero-shift of the energy-transfer axis. Six free-fitting parameters were used as well in the data fitting with the 2D diffusion model, $D_{2D}(Q)$ substituting the half-width-half-maximum of the Lorentzian. Stokes/anti-Stokes detailed balance and convolution with the instrumental energy resolution function were taken into account. Multiple scattering contribution to the spectra was neglected as the estimated sample transmission is about 89% of the incident beam.

**Data availability.** Raw data were generated at the Institut Laue-Langevin large-scale facility. Derived data supporting the findings of this study are available from the corresponding authors upon request.

Received: 10 May 2017 Accepted: 23 August 2017

Published online: 20 October 2017

## References

1. Sloan, E. D. & Koh, C. A. *Clathrate Hydrates of Natural Gases* 3rd edn (CRC Press: Boca Raton, Florida, 2008).

2. Koh, C. A., Sum, A. K. & Sloan, E. D. Gas hydrates: Unlocking the energy from icy cages. *J. Appl. Phys.* **106**, 061101 (2009).

3. Brewer, P. G., Friederich, G., Peltzer, E. T. & Orr, F. M. Jr. Direct experiments on the ocean disposal of fossil fuel $CO_2$. *Science* **284**, 943–945 (1999).

4. Luspay-Kuti, A. et al. The presence of clathrates in comet 67P/Churyumov-Gerasimenko. *Sci. Adv.* **2**, e1501781 (2016).

5. Staykova, D. K., Kuhs, W. F., Salamatin, A. N. & Hansen, T. Formation of porous gas hydrates from ice powders: Diffraction experiments and multistage model. *J. Phys. Chem. B* **107**, 10299–10311 (2003).

6. Schicks, J. M. & Ripmeester, J. A. The coexistence of two different methane hydrate phases under moderate pressure and temperature conditions: Kinetic versus thermodynamic products. *Angew. Chem. Int. Ed.* **43**, 3310–3313 (2004).

7. Schicks, J. M. & Luzi-Helbing, M. Kinetic and thermodynamic aspects of clathrate hydrate nucleation and growth. *J. Chem. Eng. Data* **60**, 269–277 (2015).

8. Vatamanu, J. & Kusalik, P. G. Unusual crystalline and polycrystalline structures in methane hydrates. *J. Am. Chem. Soc.* **128**, 15588–15589 (2006).

9. Walsh, M. R., Koh, C. A., Sloan, E. D., Sum, A. K. & Wu, D. T. Microsecond simulations of spontaneous methane hydrate nucleation and growth. *Science* **326**, 1095–1098 (2009).

10. Jiménez-Ángeles, F. & Firoozabadi, A. Nucleation of methane hydrates at moderate subcooling by molecular dynamics simulations. *J. Phys. Chem. C* **118**, 11310–11318 (2014).

11. Malolepsza, E. & Keyes, T. Pathways through equilibrated states with coexisting phases for gas hydrate formation. *J. Phys. Chem. B* **119**, 15857–15865 (2015).

12. Lauricella, M. et al. Clathrate structure-type recognition: Application to hydrate nucleation and crystallisation. *J. Chem. Phys.* **142**, 244503 (2015).

13. Jacobson, L. C. & Molinero, V. Can amorphous nuclei grow crystalline clathrates? The size and crystallinity of critical clathrate nuclei. *J. Am. Chem. Soc.* **133**, 6458–6463 (2011).

14. Moudrakovski, I. L., Sanchez, A. A., Ratcliffe, C. I. & Ripmeester, J. A. Nucleation and growth of hydrates on ice surfaces: New insights from $^{129}$Xe NMR experiments with hyperpolarized xenon. *J. Phys. Chem. B* **105**, 12338–12347 (2001).

15. Miguez, J. M. et al. Molecular dynamics simulation of $CO_2$ hydrates: Prediction of three phase coexistence line. *J. Chem. Phys.* **142**, 124505 (2015).

16. Sarupria, S. & Debenedetti, P. G. Homogeneous nucleation of methane hydrate in microsecond molecular dynamics simulations. *J. Phys. Chem. Lett.* **3**, 2942–2947 (2012).

17. Nguyen, A. H. & Molinero, V. Cross-nucleation between clathrate hydrate polymorphs: Assessing the role of stability, growth rate, and structure matching. *J. Chem. Phys.* **140**, 084506 (2014).

18. Choukroun, M., Morizet, Y. & Grasset, O. Raman study of methane clathrate hydrates under pressure: New evidence for the metastability of structure II. *J. Raman Spectrosc.* **38**, 440–451 (2007).

19. Chou, I.-M. et al. Transformations in methane hydrates. *Proc. Natl Acad. Sci. USA* **97**, 13484–13487 (2000).

20. Shu, J. et al. Structural stability of methane hydrate at high pressures. *Geosci. Front.* **2**, 93–100 (2011).

21. Kuhs, W. F., Staykova, D. K. & Salamatin, A. N. Formation of methane hydrate from polydisperse ice powders. *J. Phys. Chem. B* **110**, 13283–13295 (2006).

22. Baumert, J. et al. Lattice dynamics of methane and xenon hydrate: Observation of symmetry-avoided crossing by experiment and theory. *Phys. Rev. B* **68**, 174301 (2003).

23. Tse, J. S., Ratcliffe, C. I., Powell, B. M., Sears, V. F. & Handa, Y. P. Rotational and translational motions of trapped methane. Incoherent inelastic neutron scattering of methane hydrate. *J. Phys. Chem. A* **101**, 4491–4495 (1997).

24. Peters, B., Zimmermann, N. E. R., Beckham, G. T., Tester, J. W. & Trout, B. L. Path sampling calculation of methane diffusivity in natural gas hydrates from a water-vacancy assisted mechanism. *J. Am. Chem. Soc.* **130**, 17342–17350 (2008).

25. Choukroun, M., Grasset, O., Tobie, G. & Sotin, C. Stability of methane clathrate hydrates under pressure: Influence on outgassing processes of methane on Titan. *Icarus* **205**, 581–593 (2010).

26. Loveday, J. S. et al. Stable methane hydrate above 2 GPa and the source of Titan's atmospheric methane. *Nature* **410**, 661–663 (2001).

27. Lo, H., Lee, M.-T. & Lin, S.-T. Water vacancy driven diffusion in clathrate hydrates: Molecular dynamics simulation study. *J. Phys. Chem. C* **121**, 8280–8289 (2017).

28. Bée, M. *Quasielastic Neutron Scattering* (Adam Hilger: Bristol, 1988).

29. Komatsu, K. et al. Performance of ceramic anvils for high pressure neutron scattering. *High Press Res.* **34**, 494–499 (2014).

30. Desmedt, A., Soetens, J. C., Prager, M., Russina, M. & Ollivier, J. Dynamics of methyl iodide clathrate hydrate, investigated by MD simulations and QENS experiments. *J. Phys. Chem. C* **115**, 12689–12701 (2011).

31. Pefoute, E., Kemner, E., Soetens, J. C., Russina, M. & Desmedt, A. Diffusive motions of molecular hydrogen confined in THF clathrate hydrate. *J. Phys. Chem. C* **116**, 16823–16829 (2012).

32. Chaouachi, M., Neher, S. H., Falenty, A. & Kuhs, W. F. Time resolved coarsening of clathrate crystals: The case of gas hydrates. *Cryst. Growth Des.* **17**, 2458–2472 (2017).

33. Volino, F. & Dianoux, A. J. Neutron incoherent scattering law for diffusion in a potential of spherical symmetry: General formalism and application to diffusion inside a sphere. *Mol. Phys.* **41**, 271–279 (1980).

34. Dianoux, A. J., Volino, F. & Hervet, H. Incoherent scattering law for neutron quasi-elastic scattering in liquid crystals. *Mol. Phys.* **30**, 1181–1194 (1975).
35. Lechner, R. E. Effects of low-dimensionality in solid-state protonic conductors. *Solid State Ion.* **77**, 280–286 (1995).

36. Vidal-Vidal, À., Pérez-Rodríguez, M. & Piñeiro, M. M. Direct transition mechanism for molecular diffusion in gas hydrates. *RSC Adv.* **6**, 1966–1972 (2016).

37. Burnham, C. J. & English, N. J. Free-energy calculations of the intercage hopping barriers of hydrogen molecules in clathrate hydrates. *J. Phys. Chem. C* **120**, 16561–16567 (2016).

38. Liang, S., Liang, D., Wu, N., Yi, L. & Hu, G. Molecular mechanisms of gas diffusion in CO₂ hydrates. *J. Phys. Chem. C* **120**, 16298–16304 (2016).

39. Guo, H., Chen, Y., Lu, W., Li, L. & Wang, M. In situ Raman spectroscopic study of diffusion coefficients of methane in liquid water under high pressure and wide temperatures. *Fluid Phase Equilib.* **360**, 274–278 (2013).

40. Harris, K. R. & Trappeniers, N. J. The density dependence of the self-diffusion coefficient of liquid methane. *Physica 104A*, 262–280 (1980).

41. Van der Gulik, P. S., Mostert, R. & Van den Berg, H. R. The viscosity of methane at 25 °C up to 10 kbar. *Physica. A* **151**, 153–166 (1988).

42. Kärger, J. & Ruthven, D. M. *Diffusion in Zeolites and other Microporous Solids* (John Wiley and Sons: New York, 1992).

43. Skoulidas, A. I. & Sholl, D. S. Transport diffusivities of CH₄, CF₄, He, Ne, Ar, Xe, and SF₆ in silicalite from atomistic simulations. *J. Phys. Chem. B* **106**, 5058–5067 (2002).

44. Skoulidas, A. I., Ackerman, D. M., Johnson, J. K. & Sholl, D. S. Rapid transport of gases in carbon nanotubes. *Phys. Rev. Lett.* **89**, 185901 (2002).

45. Bagherzadeh, S. A., Alavi, S., Ripmeester, J. & Englezos, P. Formation of methane nano-bubbles during hydrate decomposition and their effect on hydrate growth. *J. Chem. Phys.* **142**, 214701 (2015).

46. Uchida, T., Yamazaki, K. & Gohara, K. Generation of micro- and nano-bubbles in water by dissociation of gas hydrates. *Korean J. Chem. Eng.* **33**, 1749–1755 (2016).

47. Yang, L., Falenty, A., Chaouachi, M., Haberthür, D. & Kuhs, W. F. Synchrotron X-ray computed microtomography study on gas hydrate decomposition in a sedimentary matrix. *Geochem. Geophys. Geosyst.* **17**, 3717–3732 (2016).

48. Uddin, M. & Coombe, D. Kinetics of CH₄ and CO₂ hydrate dissociation and gas bubble evolution via MD simulation. *J. Phys. Chem. A* **118**, 1971–1988 (2014).

49. Halpern, Y., Thieu, V., Henning, R. W., Wang, X. & Schultz, A. J. Time-resolved in situ neutron diffraction studies of gas hydrate: Transformation of structure II (sII) to structure I (sI). *J. Am. Chem. Soc.* **123**, 12826–12831 (2001).

50. Murshed, M. M., Schmidt, B. C. & Kuhs, W. F. Kinetics of methane-ethane gas replacement in clathrate-hydrates studied by time-resolved neutron diffraction and Raman spectroscopy. *J. Phys. Chem. A* **114**, 247–255 (2010).

51. Falenty, A., Qin, J., Salamatin, A. N., Yang, L. & Kuhs, W. F. Fluid composition and kinetics of the in situ replacement in CH₄–CO₂ hydrate system. *J. Phys. Chem. C* **120**, 27159–27172 (2016).

52. Hansen, T. C., Falenty, A. & Kuhs, W. F. Lattice constants and expansivities of gas hydrates from 10 K up to the stability limit. *J. Chem. Phys.* **144**, 054301 (2016).

53. Klotz, S., Philippe, J., Bull, C. L., Loveday, J. S. & Nelmes, R. J. A 3 kbar hydrogen-compatible gas loader for Paris-Edinburgh presses. *High Press Res.* **33**, 214–220 (2013).

## Acknowledgements

This work was supported by the Swiss National Science Foundation through FNS Grant 200021-149847, and by the French state funds managed by ANR within the Blanc International programme PACS (reference ANR-13-IS04-0006-01) and the Investissements d'Avenir programme (reference ANR-11-IDEX-0004-02) and more specifically within the framework of the Cluster of Excellence MATériaux Interfaces Surfaces Environnement (MATISSE) led by Sorbonne Universités. We acknowledge the Institut Laue-Langevin for provision of beam time through LTP 6-6, and Claude Payre and James Maurice for technical assistance during the experiments. We thank José Teixeira (LLB) and Robert Pick (IMPMC) for a critical reading of the manuscript.

## Author contributions

U.R., W.F.K. and A.F. prepared the sample. U.R., M.M.K., W.F.K., S.K. and L.E.B. performed the experiments. U.R., M.M.K. and L.E.B. analysed the data. All authors discussed the results and contributed to writing the manuscript.

## Additional information

Supplementary Information accompanies this paper at doi:10.1038/s41467-017-01167-2.

**Competing interests:** The authors declare no competing financial interests.

**Reprints and permission** information is available online at http://npg.nature.com/ reprintsandpermissions/

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2017