International Conference on Gas Hydrates (ICGH) (6th : 2008)

INELASTIC NEUTRON SCATTERING STUDY OF HOST AND GUEST MOLECULAR MOTIONS IN METHANE HYDRATE Kamiyama, T.; Seki, N.; Iwasa, H.; Uchida, T.; Kiyanagi, Y.; Ebinuma, Takao; Narita, Hideo; Igawa, N.; Ishii, Y.; Bennington, S.M. Jul 31, 2008

Your browser doesn't seem to have a PDF viewer, please download the PDF to view this item.

Item Metadata


59278-5758.pdf [ 999.21kB ]
JSON: 59278-1.0041073.json
JSON-LD: 59278-1.0041073-ld.json
RDF/XML (Pretty): 59278-1.0041073-rdf.xml
RDF/JSON: 59278-1.0041073-rdf.json
Turtle: 59278-1.0041073-turtle.txt
N-Triples: 59278-1.0041073-rdf-ntriples.txt
Original Record: 59278-1.0041073-source.json
Full Text

Full Text

Proceedings of the 6th International Conference on Gas Hydrates (ICGH 2008), Vancouver, British Columbia, CANADA, July 6-10, 2008.  INELASTIC NEUTRON SCATTERING STUDY OF HOST AND GUEST MOLECULAR MOTIONS IN METHANE HYDRATE T. Kamiyama∗, N. Seki, H. Iwasa, T. Uchida, Y. Kiyanagi Graduate School of Engineering, Hokkaido University Kita 13 Nishi 8, Kita-ku, Sapporo 060-8628 JAPAN T. Ebinuma, H. Narita Methane Hydrate Research Laboratory National Institute of Advanced Industrial Science and Technology Sapporo 062-8517 JAPAN N. Igawa, Y. Ishii Neutron Science Research Center, Japan Atomic Energy Research Institute Tokai, Ibaraki 319-1195 JAPAN S. M. Bennington ISIS Department, Rutherford-Appleton Laboratory Chilton, Didcot OX11 OQX UK ABSTRACT Methane hydrate has a unique structure that the host water framework forms two kinds of cages, which contain one methane molecule each. Therefore, it has been expected that there may exist three kinds of translational modes of a methane molecule and also the distortion of translational mode of host water molecules compared with normal ice. We need information of the host and guest molecular dynamics over the wide momentum and energy transfer region for studying such dynamics. In this study inelastic neutron measurements were carried under 40 K with MARI spectrometer at ISIS in UK, TAS at JRR-3 and CAT at KENS in Japan. For the methane molecular motion we could confirm its freelike rotation by complementary use of MARI and TAS spectra. After the subtraction of the scattering intensity of the rotation evaluated by the free rotation model from the experimental data, three kinds of translation modes were identified at first experimentally. On the experimental spectra there still remains the excess intensity which could not explain the single mode excitation. The libration mode of the water framework shows the different momentum and energy transfer dependence with those of normal ice. The feature of the libration mode is resemble to ice-IX, that could be considered as a proton ordering of the cage structure appeared in ice-II, VIII and IX. Keywords: methane hydrates, neutron inelastic scattering, molecular motion ∗  Corresponding author: Phone: +81 11 706 6651 Fax +81 11 706 6651 E-mail:  NOMENCLATURE BT(J): Boltzman factor at T d: Bond length of C-H E: energy transfer EJJ’: Energy transfer from J to J’ J, J’: Quantum number of rotational state jl(Qd): spherical Bessel function k, k0: wave vectors for final and incident neutrons Q: momentum transfer σ: neutron cross section T: temperature [K] Ω: solid angle INTRODUCTION Methane hydrate (MH) is a kind of nonstoichiometric inclusion compound with the host framework composed of water molecules. It consists of two kinds of water cages, two pentagonal dodecahedra and six tetrakaidecahedra, which can contain one methane molecule each. Intensive studies of MH by neutron scattering have been undertaken and have made clear the structure [1] and spectral features [2,3,4] and given the structural formation [5]. In addition to these their unique structure gives interesting physical properties such as thermal conductivity [6]. The properties can be affected its isolated guest and host structures and interactions between them. As a first approximation MH can be regarded as a discrete system of guest molecular motions, because the guest molecules, methane, are  separated each other by water cages. In our previous report [7], we revealed the neutron cross sectional features on the Q = 1 - 13 Å-1 and E = 20 - 90 meV region. There was an impressive Q2 dependent peak in the region, which was caused by the methane rotational motion. The result implied the prospect of interactions between the methane molecule and the water cage was weak, and the interactions between the methane molecules were much weaker. It was said that in the cage the methane molecular rotational motion was practically free, and we had made clear its quantum rotation in detail from the lower Q-E space measurement of dynamical structure factors [8,9]. On the other hand, there are still considered the translational motions of the methane molecule, which is strongly restricted with the water cage, and the host framework motions. In this study we reveal such modes of MH based on the wide Q region measurement. EXPERIMENTAL CH4-D2O (MH-D2O) and CH4-H2O (MH-H2O) were made at the Hokkaido National Industrial Research Institute (HNIRI) for neutron scattering samples. These were prepared by grinding crushed ice at around 268 K into contact with methane gas at 5 MPa for 72 h. The methane inclusion densities were 80% for MH-D2O and 87% for MH-H2O, respectively. Crushed ice Ih samples of H2O and D2O were prepared for comparison. Samples were  excess intensity  Figure 1. Dynamical structure factor S(Q,E) of CH4-D2O at 12 K measured by MARI (right). Calculated S(Q,E) for free rotational component of methane molecule at 12 K (left).  handled under a liquid nitrogen temperature and stored in aluminum containers. An annular geometry was used for the container to minimize self-shielding effects. A double walled can with a diameter of 42 mm was used, with the thicknesses of the MH-H2O and MH-D2O samples being 0.5 and 1.5 mm, respectively. The neutron scattering experiments were performed on the MARI spectrometer at the pulsed neutron source of the ISIS at Rutherford Appleton Laboratory in UK, the triple-axis spectrometer TAS-1 at Japan Atomic Energy Research Agency in Japan and the inverted geometry spectrometer CAT at the pulsed neutron source of the KENS at High Energy Accelerator Research Organization in Japan. The sample was attached to a cryogenic system keeping the temperature under 150 K during the sample change. The neutron scattering measurements were carried around 12 K with incident energies of 15 meV, 100 meV and 500 meV on MARI and 30 meV on TAS-1. RESULTS AND DISCUSSION CH4 free rotation model On the low energy region around meV, there appear the inelastic peaks corresponding to the methane rotation and it is suggested that there C H 4-D 2O 12K  0  d 2σ σb k = dΩdE 4π k0  ∑′ JJ  J +J′ 2J ′ + 1 2 BT (J ) ∑ jl (Qd)δ (E − EJ J ′ ) l = J −J ′ 2J + 1  (1).  Figure 1 shows the obtained spectrum by MARI (right) and the calculated one using the FRM (left). The measured dynamical structure factor S(Q,E) of MH and the methane FRM shows well agreement with their characteristic peaks, but there still remains the excess intensity of the cross-section around 4 meV toward over 8 meV. Figure 2 shows the comparison between the experimental S(Q,E) and the methane FRM on the some Q axes. The experimental features are agreed with the CH4 FRM very well. On the other hand, at the energy region over 4 meV the experimental S(Q,E) and CH4 FRM has a great difference as shown in figure 3. Such difference could not be explained by CH4 FRM, then we done the detailed analysis of the E and Q dependence of the excess mode. CH4 translation mode The subtraction of the calculated CH4 FRM from the experimental S(Q,E) leaves the excess intensity (figure 4). In the molecular dynamics (MD) calculation by Tse et. al. [2] there appears three peaks due to the translational motion of CH4 molecule in the two kinds of water cages. Fig. 4 shows an example of three Gaussian functions fitting to the subtracted intensity. The fitting  2.2meV (J:1→ 2)  Intensity [arb.unit]  Intensity [arb. unit]  1.1meV (J:0→1)  exist translational modes of methane in a small and large cages. We calculated the intensity of CH4 rotational mode by assuming the free rotation (FRM). The formula of the free rotation of a spherical top molecule is  0 3.3meV (J:0→2)  0  CH4-D2O (12[K]) Free Rotation Model Excess Intensity D2O ice  400  200  5.5me (J:1→3) 0 0  0  2  4  6  8  10  12  14  Energy transfer [meV] 0  2  -1  4  Q [Å ]  Figure 2. An example of the momentum dependent scattering intensity comparison between the measured S(Q,E) and the calculated S(Q,E) for the methane free rotation model at constant Q.  Figure 3. An example of the energy dependent scattering intensity comparison between the measured S(Q,E) and the calculated S(Q,E) for the methane free rotation model at constant Q = 6 Å-1.  procedure worked well, and the fitting results show good agreements with the subtracted spectrum. Figure 5 shows the Q dependence of the peak positions of the fitting results for MARI spectrum. The decomposed peak positions appear on the nearly constant energy around E = 3.9, 4.8 and 7.8 meV. These peaks seem to correspond to the positions of methane translational peaks in the water cages calculated by MD as 3.4 and 6 meV in large cage or 8.6 meV in small cage [2]. We think that we could assign the origin of the excess peaks at the lower energy region as such hindered translational motions of the methane molecules in the water cages.  Excess intensity over 40 meV region Figure 6 shows the S(Q,E) intensity map obtained by MARI with 100 meV incident neutrons. We can recognize that there are some unknown modes between 40 to water libration regions in MH-D2O and MH-H2O compared to the D2O and H2O ice spectra. The modes are originated by CH4 existence, because they appear just on the MH spectra. Decomposition of such unknown modes is shown in figure 7. The experimental data were taken by MARI and TAS-1, and their trends were corresponded. In the figure we placed the one of decomposed peak candidate. The assumed peaks have their tops at around 38 and 46 meV. In the region the CH4 FRM makes some very weak transition modes, then the excess modes are not by Energy levels by MD  160  Q = 3 . 0 [1/Å]  8.6[meV]  8  120  Methane Translation Peak position[meV]  Intensity[arb.unit]  140  9  100 80 60 40 20  7 6.0[meV]  6 5 4  3.4[meV]  3 2  Peak 1 Peak 2 Peak 3  1  0 2  4  6  8  10  12  Energy transfer [meV]  0 0  1  2  3  4  5  Q[Å^-1]  Figure 4. An example of three Gaussian fitting to the excess intensity, which appears after subtracting the methane free rotation model from the measured S(Q,E). This spectra (Q = 3 Å-1) was obtained by MARI at 12 K. CH4-D2O (12K)  CH4-H2O (13K)  H2O (12K)  Energy Transfer [meV]  D2O (12K)  Figure 5. Q dependence of the Gaussian fitted peaks. Energy constant dotted lines correspond to the methane translational modes calculated by MD [2].  Momentum Transfer [Å-1]  Figure 6. S(Q,E) intensity map obtained by MARI with 100 meV incident neutrons  1100 1000  Intensity[arb.unit]  900 800 700 600  CH4-D2O (MARI) D2O ice (MARI) assumed mode composed (ice + assumed) CH4-D2O (TAS-1)  500  difference between them. H2O ice has flat structure, but MH-H2O has sharp peak at 70 meV. Such peak feature of the libration mode is resemble to ice-IX, then the reason would be considered as a proton ordering of the cage structure appeared in ice-II, VIII and IX.  400 300 200 100 0  36 38 40 42 44 46 48 50 52 54 Energy transfer [meV]  Figure 7. Assumed excess modes at higher energy region. Assumed modes were made by Gaussian functions.  Figure 8. Inelastic neutron scattering spectra of CH4-H2O and H2O ice. the CH4 free rotation. One of possibilities can be said as the coupled mode between CH4 molecular modes and water lattice modes, because the region is just above the ice lattice modes between 7 and 40 meV. Unfortunately, the experimental spectra have poor resolution for detailed analysis, because the Q-E region is difficult to reach by the present spectrometers. For the confirmation of such assumption, we have to carry the more precise experiment at the new neutron facility. Water lattice mode In fig. 6 we can show the difference between the intensity distributions of librational mode of water lattice. Over 70 meV for H2O and over 50 meV for D2O the scattering intensities are almost uniform. On the other hands, MH-H2O shows the characteristic intensity distribution. Figure 8 shows the neutron spectra over wide energy region for MH-H2O and normal H2O ice. The libration mode region from 60 to 120 meV shows the large  SUMMARY We have carried on the inelastic neutron scattering measurements of methane hydrate in the wide energy and momentum range. We analyzed the obtained S(Q,E) on the assumption of methane free rotation, which shows very good consistency with the S(Q,E) under 3 meV. After subtracting the rotational component from the S(Q,E), there remained excess three peaks around 4 to 8 meV. Because these peak positions are close to the methane translational modes in the water cages, we could assign the peaks as such modes. In the higher energy region around 40 to 50 meV there were exsess modes for MH-D2O. The origin of the modes were still unknown, but one possibility was coupled mode between methane molecular modes and water lattice modes. For the host lattice mode, the libration region showed the large difference between H2O ice and MH-H2O. The libration region feature of MH is resemble to ice-IX, which could be considered as a proton ordering of the cage structure.  REFERENCES [1] Gutt C, Asmussen B, Press W, Johnson MR, Handa YP, Tse JS. The structure of deuterated methane-hydrate. J. Chem. Phys. 2000; 113:47134721. [2] Tse JS, Ratcliffe CI, Powell BM, Sears VF, Handa YP. Rotational and Translational Motions of Trapped Methane. Incoherent Inelastic Neutron Scattering of Methane Hydrate. J. Phys. Chem. A 1997; 101:4491-4495. [3] Gutt C, Asmussen B, Press W, Merkl C, Casalta H, Greinert J, Bohrmann G, Tse JS, Huller A. Quantum rotations in natural methaneclathrates from the Pacific sea-floor. Europhys. Lett. 1999; 48:269-275. [4] Gutt C, Press W, Huller A, Tse JS, Casalta H. The isotope effect and orientational potentials of methane molecules in gas hydrates. J. Chem. Phys. 2001; 114:4160-4170. [5] Koh CA, Wisbey RP, Wu X, Westacott RE, Soper AK. Water ordering around methane during  hydrate formation. J. Chem. Phys. 2000; 113:6390-6397. [6] Tse JS, White MA. Origin of glassy crystalline behavior in the thermal properties of clathrate hydrates: a thermal conductivity study of tetrahydrofuran hydrate. J. Phys. Chem. 1988; 92 :5006-5011. [7] Kamiyama T, Kiyanagi Y, Horikawa T, Iwasa H, Uchida T, Ebinuma T, Narita H, Bennington SM. Dynamical features of methane hydrate around 12 K. Physica B 2004; 350(1-3, Suppl. 1) :e395-e398 [8] Kiyanagi Y, Kamiyama T, Ohnuma S, Hiraga F, Uchida T, Ebinuma T, Narrita H, Bennington SM. Dynamical Structure of Methane Hydrate Measured by Neutron Inelastic Scattering. Proc. 5th Int. Conf. on Gas Hydrates (ICGH5, Trondheim) 2005;502-506. [9] Kamiyama T, Seki N, Iwasa H, Uchida T, Ebinuma T, Narita H, Igawa N, Ishii Y, Bennington SM, Kiyanagi Y. Methane molecular motion in clathrate hydrate host framework. Physica B 2006; 385-386:202-204.  


Citation Scheme:


Citations by CSL (citeproc-js)

Usage Statistics



Customize your widget with the following options, then copy and paste the code below into the HTML of your page to embed this item in your website.
                            <div id="ubcOpenCollectionsWidgetDisplay">
                            <script id="ubcOpenCollectionsWidget"
                            async >
IIIF logo Our image viewer uses the IIIF 2.0 standard. To load this item in other compatible viewers, use this url:


Related Items