New porous water ice metastable at atmospheric pressure obtained by emptying a hydrogen filled ice

New porous water ice metastable at atmospheric pressure obtained by emptying a hydrogen filled ice


The properties of some forms of water ice reserve still intriguing surprises. Besides the several stable or metastable phases of pure ice, solid mixtures of water with gases are precursors of other ices, since in some cases they may be emptied, leaving a metastable hydrogen bound water structure. We present here the first characterization of a new form of ice, obtained from the crystalline solid compound of water and molecular hydrogen called C filled ice. By means of Raman spectroscopy, we measure the hydrogen release at different temperatures, and succeed in rapidly removing all the hydrogen molecules, obtaining a new form of ice (ice XVII). Its structure is determined by means of neutron diffraction measurements. Of paramount interest is that the emptied crystal can adsorb again hydrogen and release it repeatedly, showing a temperature dependent hysteresis.

  • Consiglio Nazionale delle Ricerche, Istituto dei Sistemi Complessi, via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, via Sansone 1, I-50019 Sesto Fiorentino, Italy.
    *Corresponding author.

  • Introduction

    Water molecules in the solid state may give rise to more than fifteen different forms of ices, depending on temperature and pressure[(1)]. In addition, when water freezes in the presence of some other molecular substance that does not bind chemically to water, it may form peculiar crystal structures, known as clathrate-hydrates, trapping guest molecules inside cages of different geometries[(2)]. Until recently, it has been believed that caged guest molecules are essential for the stability of the clathrate-hydrate crystals, so that the water skeleton would collapse without them. One experimental study has recently demonstrated, however, that at least one of these clathrate structures can be emptied of its guests, by letting neon atoms diffuse out of the solid, and persist in a metastable state if preserved at low temperature[(3)]. These and other similar low-density lattices of water molecules have raised recently large interest from a theoretical point of view because they are believed to be the stable phase of solid water at negative pressure[(4)].

    Water does form crystalline compounds with molecular hydrogen as well. Four ordered structures of this binary mixture are known. The first in order of increasing pressure, stable at MPa and at 4-6 C below zero, is a clathrate-hydrate[(5), (6)] having the so called cubic sII structure, common to other clathrates with different molecular guests[(7)]. In this non-stoichiometric compound hydrogen molecules are trapped in two types of cages, with a total hydrogen molar fraction ( mol(H)/mol(HO)) up to about 35 %. The hydrogen molecules in the cages perform a peculiar quantum rattling and rotational motion, which has been studied efficiently with inelastic neutron scattering[(8), (9), (10), (11), (12)] and Raman scattering[(13), (14), (15), (16)]. The latter technique is a very powerful and convenient one to identify the interactions of the molecules with the water environment, and to measure the composition of the sample. Besides sII clathrates, the existence of two other stable phases for the H-HO solid mixture at higher pressure ( 700 MPa) is known since 1993[(17)]. These two structures, indicated with C and C, do not possess the typical cage structure of clathrates, and are usually and more properly named filled ices. Recently, only two experimental studies[(18), (19)] have investigated the phase diagram of the H-HO compounds at intermediate pressures and have demonstrated the presence of a further stable phase, named C, at temperatures 100-270 K and pressures 360-700 MPa, that is intermediate between the stability region of the sII clathrate and C phase (see Fig. 1). Up to now, a definitive consensus on the structure of the C phase has not been reached. The structural model (later referred to as C-I) originally proposed in ref.  (18) after ex-situ x-ray diffraction measurements at room pressure and 80 K, assumes space group (or ) and the presence of water molecules in sites with 0.5 occupancy. These water molecules would not present correct hydrogen bonds with the other molecules of the lattice. Strobel(19) suggests two possible structures (indicated as -quartz and sT) after the analysis of x-ray diffraction collected from one sample pressurized in a diamond anvil cell. These models were later examined theoretically by Smirnov and Stegailov(24), who considered also a variant of the C-I structure (indicated by C-II), and suggested C-II and sT both as viable candidates. Actually, the -quartz structure, with the proposed parameters(24) gives rise to atypical, very short, O-O distances. After this, Oganov and coworkers(25) used DFT-based structure searching to predict several phases of the H-HO mixture, but these searching methods are restricted in the number of atoms, so any large clathrate-sized unit cells can easily be missed. In this work we study, by means of Raman spectroscopy, the C phase of the water-hydrogen mixture, produced at about 400 MPa and recovered at room pressure and liquid nitrogen temperature, still containing a large fraction of molecular hydrogen. We demonstrate that, by means of a thermal treatment under vacuum, all the hydrogens can be removed and that the new form of ice so obtained, named ice XVII, is metastable at room pressure below 120 K. The structure of the sample is checked before the removal of the hydrogen (i.e. in the C phase) and after it, showing that minimal changes, if any, are produced by the outgassing. The arrangement of water molecules in ice XVII gives rise to spiraling channels parallel to the crystallographic axis, having a diameter of about 6.10 Å making this ice a porous material. As a matter of fact, we show in this work that ice XVII can adsorb and release hydrogen gas over and over, without evident change of structure.


    Sample production
    For this study we produce several batches of the solid H-HO compound in the C phase and recover the samples at room pressure and liquid nitrogen temperature, by application of a standardized and reproducible procedure. The structure of some specimens is examined by x-ray diffraction at the CRIST laboratory of the University of Firenze. The sample examined presents a considerable texture, which prevents a Rietvield analysis, but the measured diffraction pattern can be used to discriminate among several proposed structures. A Le Bail fit of the diffraction pattern is presented in Fig. 2, assuming the C-II structure, space group . Lattice constant are Å and Å. The other structures considered, namely sT (), Ih-C ()(25) and ice I (), would give diffraction patterns in evident disagreement with the experimental one. Details of the synthesis procedure and of the x-rays measurements are given in the Methods section.

    Measurement and interpretation of the Raman spectra
    By using our cryogenic Raman apparatus(26) we measure spectra during different cycles, observing the crystal lattice excitations and the rotations of the H molecule (150-650 cm), the OH stretching mode of the water molecule (3000-3400 cm), and the vibron of the H molecule (4100-4200 cm). Our Raman spectra (see Fig. 3) and overall results add important information on the structure, which are consistent with only one of the proposed structures, namely the C-II one(18); (24) with space group depicted in Fig. 1. In this structure the water molecules form spiraling channels with a free bore hole along the axis of about 5.26 Å and with a diameter of 6.10 Å, which can accommodate the H molecules. Observing the spectra (Fig. 3a), we notice that the lattice phonon band (black line) has a broad smooth shape, showing similarities with the same band in both ice Ih (blue line) and sII clathrates (red line). The absence of sharp lines rules out the possibility that the water lattice in the C phase might be proton ordered. This is at variance with the C filled ice, which exhibits proton order and, consequently, a lattice phonon spectrum with evident sharp lines(17). The rotational spectrum of the hydrogen molecules (Fig. 3c) presents well distinct lines, proving that the H molecules rotate almost freely. The S(0) rotational band is split in probably three components, with a larger splitting than in clathrate-hydrates. This is an indication of a more intense interaction of the hydrogen molecule with the environment and a larger perturbation of the rotational motion. According to the proposed structure, the hydrogen molecules are arranged in the channels, probably in a spiraling configuration(18), at a distance from water oxygen atoms of about 3.1 Å and at about 2.95 Å between each other, the same as in solid hydrogen at about 1.5 GPa. The larger splitting of the rotational band components is accountable on the basis of the shortest distance between H and HO molecules in the structure, and indirectly confirms it. In addition, at such a short distance the H-H anisotropic interaction may sustain collective rotational excitations, as it happens in solid H(27) and Ar(H) high pressure compound(28); (29), which may contribute to the width of the rotational lines. The spectrum of the H vibron region (Fig. 3d) consists of only one doublet (Q(0) and Q(1) lines), and this is in accordance with the occupation of a single crystallographic site for the hydrogen molecules, as in the proposed structure(18). The OH stretching band (Fig. 3b) has the main peak at 3123 cm, a frequency higher than that of both sII clathrate and ice Ih. The OH vibrational frequency is known to decrease while decreasing the O-O distance in different ices. This is observed both as a function of the considered compound(30) or, for the same structure, as a function of pressure, as reported for ice VII(31) and C filled ice(32). The higher frequency of the OH stretching mode observed for the C filled ice is in accordance with the larger O-O distance in the structure ( Å ) than in sII clathrate ( Å) and in ice Ih ( Å).

    Raman band intensities and hydrogen content
    We derive the hydrogen molar fraction from the intensity ratio of the hydrogen rotational lines and the lattice phonon band, . This analysis is done thanks to the assumption, that sounds obvious, that the Raman intensity of a band is proportional to the number of molecules giving rise to it. This is anyhow an approximation, since difference of the polarisability of the same molecules in different environments may lead to tiny differences in intensity. This matter has been considered in sec. IIIa of Ref.  (16), and we refer to that paper for a more thorough discussion. We calculate the calibration factor by fitting data arising from two independent experimental methods. A first set of data is obtained by using similar spectra of sII clathrates, for which the hydrogen content is calculated by counting the number of molecules in the large and small cages(16). The spectra of these samples in the 200-900 cm range (not reported in ref.  (16)) are used to calculate . To these data we add the information we derive from new volumetric measurements performed on our C sample at three temperatures, namely 20, 50 and 80 K. The procedure is explained in detail in the Methods section. For the samples examined just after synthesis, whose typical spectra are represented in Fig. 3, we obtain %.

    Emptying of C filled ice
    We discover that the C structure is metastable at room pressure even when emptied, and is a new form of ice, namely ice XVII. The mechanical stability of the sample is tested by increasing the temperature while keeping the sample under dynamic vacuum. We then observe the gradual release of the hydrogen from the sample, up to the complete undetectability of the hydrogen rotational lines obtained after pumping for 1-2 hours at a temperature of about 120 K. Considering the sensitivity and the signal to noise ratio of our detection system, we can ascertain that the H/HO molar fraction in the emptied sample is less then 0.5%. During the heating process, we do not observe any abrupt change of the lattice phonon and OH stretching bands, demonstrating that no structural phase transitions have occurred. Decreasing again the temperature of ice XVII and comparing the spectra measured at the same low temperature we notice after the annealing an increase in the frequency of about 7 cm for the lattice phonon bands (Fig. 4a) and a decrease of about 20 cm for the OH mode (Fig. 4b). This indicates a decrease of the average O-O distance. The same indication we obtain from the concurrent increase of the lattice modes, indicating a stronger binding. However, the overall similarity of the band shape before and after annealing leads us to believe that the structure of the filled C ice and that of ice XVII are essentially the same. The situation for the C structure is different from previous observation in sII, where upon cages emptying, the lattice constant, and consequently the O-O distance increases(3). The contraction of the water framework upon inclusion of Ne in the cages of sII clathrate may be imputed to an attractive interaction between this atom and water. It is interesting to compare also with the situation of the He-hydrate, which become ice II upon emptying(33). Here the phenomenology is more complex, but the changes in the structure can be related to the mainly repulsive interaction between the He and O atoms. The precise nature of the effect of the annealing on the structure of the C samples is not clear yet. The increased sharpness of the spectral features measured after the annealing suggests the decrease of defects in the structure. We recall that recently some authors(34), in order to explain the non perfect refinement of the x-rays diffraction pattern with the structure, have hypothesized that some nitrogen molecules, instead of water molecules, may have been trapped in the channels during the recovering of the sample in liquid nitrogen. A possible effect of the annealing might be the removal of these nitrogen molecules. To check this hypothesis, we measure Raman spectra of our samples, before and after this annealing process, also in the region of the N stretching mode, observing, at 20 K before the annealing, a quite evident and sharp line at a frequency cm. By virtue of their vibration frequency, sensibly lower than both that of N gas (2329.917 cm)(35), and that of solid N in the -phase (a doublet at 2327.5 and 2328.5 cm)(36), we can establish that the N molecules are trapped inside the structure. The Raman line of nitrogen molecules essentially disappears (it becomes 1000 times weaker) after the annealing process. The evident spectroscopic signatures of structural modification described above are a combined effects of the removal of both H and N from the sample. The structure of deuterated ice XVII is determined by means of neutron diffraction measurements performed on OSIRIS at ISIS, RAL (UK)(37). Reitveld refinement of the data enables us to determine the structure. The empty structure is described by the space group P622, with oxygen and deuterium atoms in and positions respectively. The space group used for the refinement has a higher symmetry than that proposed for the filled ice C(24) on the basis of the data in ref. (18). This may be due to the presence in the channels of either H guests and N impurities. A complete description of the neutron diffraction experiment and fit procedure, together with the fitted structural parameters and water molecules geometry are reported in ref. (37).

    Refilling of ice XVII
    The unexpected property of ice XVII is that, when exposed to even a low pressure of hydrogen gas, it adsorbs the molecules up to an amount which is pressure dependent, but may grows larger than that initially present just after synthesis. As a matter of fact, we generally obtain % for pristine samples, possibly because of the presence of adsorbed nitrogen in the channel and/or of the handling of the sample for the insertion of the Raman apparatus. The adsorption process has initially a fast kinetics, even at 15 K. The rotational Raman spectra measured at 40 K after subsequent partial refilling steps are shown in Fig. 4c. Maximum gas pressure is in this case only 250 mbar, but the rotational intensity grows higher than that in the pristine C sample. The shape of the rotational lines does not change significantly, demonstrating that H is penetrating again in the same positions in the channels as those initially occupied at the time of synthesis. Normalizing the rotational Raman intensity with respect to the lattice phonon band and applying the calculated conversion constant, we estimate the H/HO molar fraction. We notice that, in some instances, as for the most intense spectrum shown in Fig 4c, the measured H/HO molar ratio in the sample exceeds 40 %, reaching higher values than the established theoretical maximum () for sII hydrates.

    We measure then several adsorption isotherms, deriving the amount of adsorbed H as a function of pressure, at several different temperatures between 15 to 80 K. Results for four temperatures are presented in Fig. 5. Pressure is increased gradually (squares in Fig. 5) and, when possible, measurements are taken after the same delay from pressure increase, which is indicated in each figure. The pressure at which saturation is reached is strongly dependent on temperature and ranges between a few millibar at 15 K to several bar at 80 K. By measuring with the same method, we also examine the gas desorption while releasing pressure. A large hysteresis and other kinetics effects are evident at the lower temperatures. At 15 K, after decreasing pressure down to zero, the sample does not sensibly release hydrogen in min. The same effect is observed at 40 K, the filling remaining almost constant for at least 36 min. At this temperature, the increase in the adsorbed gas observed at about 240 mbar, after a 16.5 hour long exposure to hydrogen, indicates that at least two time scales dominate the phenomenon: one fast adsorption in the outermost sites of the channels and a slower one, limited by molecular diffusion inside the channels. At higher temperatures, 65 and 80 K, we observe a very small hysteresis, indicating that the time needed to hydrogen to diffuse for distance of the order of the linear dimension of the sample grains is quite shorter than the time interval between measurements. At temperature above 50 K, the desorption process is quick (compared to sII clathrates). This is not surprising at all, since the channels in ice XVII are much wider than the free space inside both the 6-membered and 5-membered rings. From the isotherm hydrogen uptake as a function of the pressure presented in fig. 5, it is possible to estimate, by means of the Vant Hoff equation, the enthalpy of adsorption, , for a given hydrogen uptake :


    where is the gas constant.(38); (39) For this calculation we use the data collected at 50, 65 and 80 K, estimating the derivative with a straight line fit. The resulting values of decrease with increasing molar ratio starting from about 5 kJ/mol at down to about 2 kJ/mol at , in agreement with typical values for physisorption.


    This work is the first characterization of the dynamics of the C filled ice. Moreover, two remarkable results arise from this study. First, we discover that the C structure is metastable when emptied. Due to its large range of stability (up to 120 K at least), it should be counted as a new form of ice and named ice XVII. Second, we find that ice XVII is capable to adsorb hydrogen gas in large quantities (presumably, on the basis of its structure up to ), without substantial change in its structure. We believe that, more than the maximum estimated hydrogen uptake, it is the fast reversibility of gas adsorption, the theoretically infinite number of possible cycles (there is no chemical change of the substance involved) and the relatively modest pressures at which the process occurs, at temperature close to that of liquid nitrogen, that make ice XVII really appealing for hydrogen storage applications.


    Procedure for the synthesis of the sample
    We produce several instances of solid H-HO compound by application of a standardized procedure. The check of the recovered samples with Raman scattering confirms that we have always obtain the sample in the same C phase. To synthesize the samples, we introduce a few grams of finely ground ice into a beryllium-copper autoclave, and we expose the sample to hydrogen gas at a pressure in the range or above 430 MPa, at a temperature of about 255 K (-18C). We prudentially leave the sample under pressure for a few days, even though, while rising pressure during synthesis, we observe at about 360 MPa an abrupt decrease of the rate of pressure rise, which is a clear indication of large hydrogen adsorption by the ice. Then we quench the autoclave in liquid nitrogen when still under pressure, and finally, after releasing the pressure, we recover the sample at liquid nitrogen temperature, in the form of a fine powder. The same procedure is applied for the synthesis of the samples to be investigated by x-rays diffraction, except that in this case a drop of water is frozen from the beginning inside a few common x-rays glass capillaries, having 0.5 mm diameter and a length of about 10-15 mm. The frozen water occupies a capillary length of about 2 mm. The capillaries are inserted in the autoclave where they undergo the same standardized procedure, and are then recovered and handled at liquid nitrogen temperature.

    For Raman measurements, a small amount of sample is inserted in our optical cell, in contact with the cold finger of a closed-cycle He cryostat(26). The transfer process has to be accomplished at almost liquid nitrogen temperature and in a dry-nitrogen atmosphere, exerting particular care to avoid sample heating. Once the Raman cell is filled with the sample, it is sealed with the optical window and purged with helium gas. The sample cell can then be set to any temperature between room temperature and K. Gas pressure in the cell can be varied and measured with an accuracy of about 1 %.

    Raman and x-ray measurements
    Raman spectra are excited by means of an Ar ion laser at 514.5 nm, with a power of about 30 mW on the sample, focused on a spot of about 30 m. Scattered light is collected in an almost back-scattering geometry, focused on the entrance slit of a Spex spectrograph, and recorded by a cooled CCD detector (Andor). The maximum resolution of the instrument is 0.4 cm. Similar care has to be exerted to mount the capillary on the x-ray diffractometer (four circle Oxford Diffraction XcaliburPX) The capillary is maintained at low temperature (about 90-100 K) by a cold nitrogen blow provided by a Oxford Cryosteam system. Diffractograms are collected using a copper x-ray source, and recorded by a 165 mm diameter CCD detector. Total collection time is of the order of 7 minutes. During collection, the capillary is turning rapidly.

    Calibration of the Raman intensity
    We measure the molar fraction of hydrogen molecules contained in the samples of filled ice by computing the ratio of the intensity of the S(0) and S(1) rotational lines, weighted by the inverse of the respective cross section with the intensity of the lattice phonon band. Specifically, if we indicate with and the Raman cross section for the S(0) and S(1) rotational lines, and with and their intensity, respectively, the quantity is proportional to the number of H molecules in the lowest two and rotational states, that is, at low temperature, to the number of H molecules giving Raman scattering. Dividing by the intensity of the water lattice phonon band, , we obtain a quantity which is independent of laser intensity and detector efficiency, and proportional to the H/HO molar ratio in the sample. This procedure needs a calibration, that we have accomplished fitting data obtained with two independent methods. For several samples of sII clathates studied in our laboratory, the hydrogen content has been indirectly derived, by counting the molecules in the small and large cages(16), while the spectra of the same samples in the 200-900 cm range are in our records, even if not reported in ref.  (16). This allows us to draw some calibration points (black squares) in Fig. 6. We have performed additional volumetric measurements on our C sample while measuring, with Raman scattering, the adsorption isotherms at three temperatures, namely 20, 50 and 80 K. These results are reported as colored dots in Fig. 6. The estimated uncertainty on these experimental points depends critically on the precision of the different calibrated volumes, of the pressure transducer used and on the pressure of the measurement. The cumulative linear fit gives a calibration factor of , that is the value we have used in this work.

    Data availability
    The data that support the findings of this study are available from the corresponding author upon reasonable request.

    Figure 1: Stability region of the C phase. Projection on the plane of the phase diagram of the water-hydrogen mixture, for large H concentration, in the region of interest. The solid black lines summarize the available experimental data(17); (20); (21); (22); (5); (19); (18) and do not have the same accuracy. The dashed black line shows the possible location of the C-C transition(19). For comparison, also the phase diagram of pure water is shown (thin cyan lines)(1); (23). Drawings of the molecular arrangement in the C phase are pictured according to Ref.  (18) and to the results discussed later on, both with and without filling of hydrogen molecules, which are represented as dumbbell in a random orientation. The black stars represent the typical thermodynamic conditions for the synthesis of our samples.
    Figure 2: X-ray diffraction pattern. Experimental x-ray diffraction pattern measured at ambient pressure and temperature of about 100 K using Cu K radiation (black crosses). The red line represents the fit, performed with Le Bail method, to the experimental data using the C-II structure, space group . The value of the lattice parameters obtained from the fit procedure are reported in the figure. The blue curve is the difference between the experimental pattern and the fit.
    Figure 3: Raman spectra measured at 30 K. Spectra measured for the C sample (black lines), compared to those measured in sII clathrate(16) (red line) and ice Ih (blue line). The lattice phonon band in the range 120-310 cm (a) is at lower frequencies, while the OH stretching band (b) has the main peak at higher frequency with respect to ice and clathrates, as a manifestation of a larger O-O distance in the structure. The two rotational excitations, S(0) and S(1) of the hydrogen molecule, (c) are broader than in sII clathrates, indicating a stronger interaction with the host. A single pair of lines (Q(0) and Q(1)) observed for the H vibron attests the existence of a single crystallographic site for the H molecule in the C crystal, oppositely to the sII clathrate, where the different peaks correspond to different environments for the H molecules (d).
    Figure 4: Effect of heating and hydrogen release. The spectra of two different pristine samples measured at 30 K (black and grey lines in all panels) show some different filling (see c), probably due to unintentional heating during handling, but hardly any difference, or a very small one, in the lattice phonon (a) and OH stretching band (b) respectively. The empty C ice, i.e. ice XVII, cooled again at 20 K after annealing, (blue lines) show shifted and narrower spectra, because of increased bond strength between molecules and decrease of defects. Filling again ice XVII at 40 K does not alter much the spectra of the water molecules (green lines in a and b). Red lines in c show the spectra of the hydrogen molecules for the maximum filling obtained in this work ().
    Figure 5: Adsorption isotherms measured at different temperatures. We report the hydrogen molar fraction measured while increasing (square symbols) and decreasing (triangular symbols) hydrogen pressure in annealed ice XVII. Due to calibration, has a systematic uncertainty of 10 %. The time lag between pressure increase in the cell and Raman measurement in same series (same symbol) is indicated as t in each panel, while we have joined with a thin dotted line measurements performed with larger delays. The kinetic effects that produce the large hysteresis observed at the two lower temperatures are absent or less evident increasing the temperature.
    Figure 6: Raman intensity calibration. Summary of the determinations of the H/HO molar ratio plotted as a function of the concomitant measurements of the / ratio, used to determine the calibration factor. The solid black squares represent measures of several instances of sII pure H clathrate-hydrates, synthesized at different pressures and/or undergone to different thermal cycles(16), while the other points joined by straight lines are volumetric measures performed in this work during the measurements of the 20, 50 and 80 K isotherms (red, green and blue respectively). Different symbols (solid dots or squares) indicate the use of a different standard volume and/or a different transducer. The solid black line is the linear fit through zero, with slope 0.233 and the two dashed lines, with slopes 0.233 0.023, represent the confidence interval. Error bars are estimates of maximum uncertainty as explained in the text.
  • Acknowledgments

    The authors thank all the personnel at the Centro di Cristallografia of the Università degli Studi di Firenze (CRIST) for their assistance in the x-ray measurements. We are also grateful for the support of Francesco Grazzi (CNR-ISC, Firenze) in the analysis of the x-ray and neutron diffraction data. Competing Interests The authors declare that they have no competing financial interests. Author Contribution L. d. R. performed the high pressure synthesis of the sample, the low temperature Raman measurements and the x-rays and neutron diffraction measurements. M. C. set up the high pressure apparatus for the sample synthesis, performed the synthesis of the samples for Raman scattering, x-ray and neutron diffraction, and executed the x-ray and neutron diffraction measurements. L.U. designed the study, discussed data collection and analysis and wrote the paper. Correspondence Correspondence and requests for materials should be addressed to L.U. (email:

  • References

    1. Bartels-Rausch, T. et al. Ice structures, patterns, and processes: A view across the icefields. Rev. Mod. Phys. 84, 885-944 (2012).
    2. Sloan, E. D. & Koh, C. A. Clathrate Hydrates of Natural Gases, 3rd Edition (Taylor & Francis, New York, 2008).
    3. Falenty, A., Hansen, T. C. & Kuhs, W. F. Formation and properties of ice XVI obtained by emptying a type sII clathrate hydrate. Nature 516, 231-233 (2014).
    4. Huang, Y. et al. A new phase diagram of water under negative pressure: The rise of the lowest-density clathrate s-III. Science Advances 2, e1501010 (2016).
    5. Dyadin, Y. A. et al. Clathrate hydrates of hydrogen and neon. Mendeleev Commun. 9, 209-210 (1999).
    6. Mao, W. L. et al. Hydrogen clusters in clathrate hydrate. Science 297, 2247-2249 (2002).
    7. Mak, T. C. & McMullan, R. K. Polyhedral clathrate hydrates. X. Structure of the double hydrate of tetrahydrofuran and hydrogen sulfide. J. Chem. Phys. 42, 2732-2737 (1965).
    8. Ulivi, L. et al. Quantum rattling of molecular hydrogen in clathrate hydrate nanocavities. Phys. Rev. B 76, 161401 (2007).
    9. Xu, M., Ulivi, L., Celli, M., Colognesi, D. & Bačić, Z. Quantum calculation of inelastic neutron scattering spectra of a hydrogen molecule inside a nanoscale cavity based on rigorous treatment of the coupled translation-rotation dynamics. Phys. Rev. B 83, 241403 (2011).
    10. Xu, M., Ulivi, L., Celli, M., Colognesi, D. & Bačić, Z. Rigorous quantum treatment of inelastic neutron scattering spectra of a heteronuclear diatomic molecule in a nanocavity: HD in the small cage of structure II clathrate hydrate. Chem. Phys. Lett. 563, 1 (2013).
    11. Celli, M. et al. Experimental inelastic neutron scattering spectrum of hydrogen hexagonal clathrate-hydrate compared with rigorous quantum simulations. J. Chem. Phys. 139, 164507 (2013).
    12. Colognesi, D., Celli, M., Ulivi, L., Xu, M. & Bačić, Z. Neutron scattering measurements and computation of the quantum dynamics of hydrogen molecules trapped in the small and large cages of clathrate hydrates. J. Phys. Chem. A 117, 7314-7326 (2013).
    13. Giannasi, A., Celli, M., Ulivi, L. & Zoppi, M. Low temperature Raman spectra of hydrogen in simple and binary clathrate hydrates. J. Chem. Phys. 129, 084705 (2008).
    14. Strobel, T. A., Sloan, E. D. & Koh, C. A. Raman spectroscopic studies of hydrogen clathrate hydrates. J. Phys. Chem. 130, 014506 1-10 (2009).
    15. Giannasi, A., Celli, M., Zoppi, M., Moraldi, M. & Ulivi, L. Experimental and theoretical analysis of the rotational Raman spectrum of hydrogen molecules in clathrate hydrates. J. Chem. Phys. 135, 054506 (2011).
    16. Zaghloul, M. A., Celli, M., Salem, N. M., Elsheikh, S. M. & Ulivi, L. High pressure synthesis and in situ Raman spectroscopy of H and HD clathrate hydrates. J. Chem. Phys. 137, 164320 (2012).
    17. Vos, W. L., Finger, L. W., Hemley, R. J. & Mao, H. K. Novel H+HO clathrates at high pressures. Phys. Rev. Lett. 71, 3150-3153 (1993).
    18. Efimchenko, V. S. et al. New phase in the water-hydrogen system. J. Alloys Comp. 509, S860-S863 (2011).
    19. Strobel, T. A., Somayazulu, M. & Hemley, R. J. Phase behavior of H+HO at high pressures and low temperatures. J. Phys. Chem. C 115, 4898-4903 (2011).
    20. Lokshin, K. A. & Zhao, Y. Fast synthesis method and phase diagram of hydrogen clathrate hydrate. App. Phys. Lett. 88 131909 (2006).
    21. Efimchenko, V. S., Antonov, V. E., Barkalov, O. I., Beskrovnyy, A. I., Fedotov, V. K. & Klyamkin, S. N. Phase transitions and equilibrium hydrogen content of phases in the water–hydrogen system at pressures to 1.8 kbar. High Press. Res., 26 439-443, (2006).
    22. Antonov, V. E., Efimchenko, V. S. & Tkacz, M. Phase Transitions in the Water-Hydrogen System at Pressures up to 4.7 kbar. J. Phys. Chem. B 113, 779-785 (2009).
    23. Choukroun, M. & Grasset, O. Thermodynamic model for water and high-pressure ices up to 2.2 GPa and down to the metastable domain. J. Chem. Phys. 127, 124506 (2007).
    24. Smirnov, G. S. & Stegailov, V. V. Toward determination of the new hydrogen hydrate clathrate structures. J. Phys. Chem. Lett. 4, 3560-3564 (2013).
    25. Qian, G. R., Lyakhov, A. O., Zhu, Q., Oganov, A. R. & Dong, X. Novel Hydrogen Hydrate Structures under Pressure. Sci. Reports 4, 5606 (2014).
    26. Giannasi, A., Celli, M., Grazzi, F., Ulivi, L. & Zoppi, M. An apparatus for simultaneous thermodynamic and optical measurements, with large temperature excursions. Rev. Sci. Instrum. 79, 013105 (2008).
    27. Van Kranendonk, J. Solid Hydrogen (Plenum Press, New York, 1983).
    28. Grazzi, F., Santoro, M., Moraldi, M., & Ulivi, L. Anisotropic interactions of hydrogen molecules from the pressure dependence of the rotational spectrum in the Ar(H) compound. Phys. Rev. Lett. 87, 125506 (2001).
    29. Grazzi, F., Santoro, M., Moraldi, M. & Ulivi, L. Roton excitations of the hydrogen molecule in the Ar(H) compound. Phys. Rev. B 66, 144303 (2002).
    30. Pimentel, G. C. & Sederholm, C. H. Correlation of infrared stretching frequencies and hydrogen bond distances in crystals. J. Chem. Phys. 24, 639-641 (1956).
    31. Pruzan, P. Pressure effects on the hydrogen bond in ice up to 80 GPa. J. Mol. Struct. 322, 279-286 (1994).
    32. Vos, W. L., Finger, L. W., Hemley, R. J. & Mao, H. K. Pressure dependence of hydrogen bonding in a novel H+HO clathrate. Chem. Phys. Lett. 257, 524-530 (1996).
    33. Lobban, C., Finney, J. L. & Kuhs, W. F. The p-T dependency of the ice II crystal structure and the effect of helium inclusion. J. Chem. Phys. 117, 3928-3934 (2002).
    34. Smirnov G.S. & Stegailov V. V. Anomalous diffusion of guest molecules in hydrogen gas hydrates. High Temperature 53, 829-836 (2015).
    35. Bendtsen, J. The rotational and rotation-vibrational Raman spectra of N, NN and N. J. Raman Spectr. 2, 133-145 (1974).
    36. Ouillon, R., Turc, C., Lemaistre, J. P. & Ranson, P. High resolution Raman spectroscopy in the and crystalline phases of N. J. Chem. Phys., 93, 3005-3011 (1990).
    37. del Rosso, L., Grazzi, F., Celli, M., Colognesi, D. & Ulivi, L. Refined structure of metastable ice XVII from neutron diffraction measurements. J. Phys. Chem. Lett., (submitted) 2016. Preprint at (2016).
    38. Do, D. D. Adsorption Analysis: Equilibria and Kinetics. (Imperial College press, UK, 1998).
    39. Schlichtenmayer, M. Hirscher, M. The usable capacity of porous materials for hydrogen storage. Appl. Phys. A 122, 379 (2016).
    This is a comment super asjknd jkasnjk adsnkj
    The feedback cannot be empty
    Comments 0
    The feedback cannot be empty
    Add comment

    You’re adding your first comment!
    How to quickly get a good reply:
    • Offer a constructive comment on the author work.
    • Add helpful links to code implementation or project page.