New porous water ice metastable at atmospheric pressure obtained by emptying a hydrogen-filled ice (original) (raw)

Introduction

Water molecules in the solid state may give rise to more than 15 different forms of ices, depending on temperature and pressure1. 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 geometries2. 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 temperature3. 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 pressure4.

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 _P_>100 MPa and at 4–6 °C below 0, is a clathrate hydrate5,6 having the so-called cubic sII structure, common to other clathrates with different molecular guests7. In this non-stoichiometric compound, hydrogen molecules are trapped in two types of cages with a total hydrogen molar fraction X (_X_=mol(H2)/mol(H2O)) up to ∼35%. The hydrogen molecules in the cages perform a peculiar quantum rattling and rotational motion, which has been studied efficiently with inelastic neutron scattering8,9,10,11,12 and Raman scattering13,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 H2–H2O solid mixture at higher pressure (_P_≳700 MPa) is known since 1993 (ref. 17). These two structures, indicated with C1 and C2, do not possess the typical cage structure of clathrates and are usually and more properly named filled ices. Recently, only two experimental studies18,19 have investigated the phase diagram of the H2–H2O compounds at intermediate pressures and have demonstrated the presence of a further stable phase, named C0, at temperatures 100–270 K and pressures 360–700 MPa, which is intermediate between the stability region of the sII clathrate and C1 phase (see Fig. 1). Up to now, a definitive consensus on the structure of the C0 phase has not been reached. The structural model (later referred to as C0-I) originally proposed in ref. 18 after ex-situ X-ray diffraction measurements at room pressure and 80 K, assumes space group _P_3112 (or _P_3121) and the presence of water molecules in sites 3_a_2 with 0.5 occupancy. These water molecules would not present correct hydrogen bonds with the other molecules of the lattice. Strobel et al.19 suggest 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 Stegailov20, who considered also a variant of the C0-I structure (indicated by C0-II), and suggested C0-II and sT′ both as viable candidates. In fact, the _α_-quartz structure, with the proposed parameters20, gives rise to atypical, very short, O–O distances. After this, Oganov and colleagues21 used density functional theory-based structure searching to predict several phases of the H2–H2O mixture, but these searching methods are restricted in the number of atoms; thus, any large clathrate-sized unit cells can easily be missed. In this work we study, by means of Raman spectroscopy, the C0 phase of the water–hydrogen mixture, produced at ∼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 (that is, in the C0 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 c axis, having a diameter of ∼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.

Figure 1: Stability region of the C0 phase.

figure 1

Projection on the P_−_T plane of the phase diagram of the water–hydrogen mixture, for large H2 concentration, in the region of interest. The solid black lines summarize the available experimental data5,17,18,19,36,37,38 and do not have the same accuracy. The dashed black line shows the possible location of the C0–C1 transition19. For comparison, also the phase diagram of pure water is shown (thin cyan lines)1,39. Drawings of the molecular arrangement in the C0 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.

Full size image

Results

Sample production

For this study we produce several batches of the solid H2–H2O compound in the C0 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 C0-II structure, space group _P_3112. Lattice constant are _a_=6.3313±0.0002 Å and _c_=6.1058±0.0002 Å. The other structures considered, namely sT′ (P_42/mnm), Ih-C0 (Cc)21 and ice I_c (Fdm), would give diffraction patterns in evident disagreement with the experimental one. Details of the synthesis procedure and of the X-ray measurements are given in the Methods section.

Figure 2: X-ray diffraction pattern.

figure 2

Experimental X-ray diffraction pattern measured at ambient pressure and temperature of ∼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 C0-II structure, space group _P_3112. 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.

Full size image

Measurement and interpretation of the Raman spectra

By using our cryogenic Raman apparatus22 we measure spectra during different P_−_T cycles, observing the crystal lattice excitations and the rotations of the H2 molecule (150–650 cm−1), the OH stretching mode of the water molecule (3,000–3,400 cm−1) and the vibron of the H2 molecule (4,100–4,200 cm−1). 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 C0-II one18,20 with space group _P_3112 depicted in Fig. 1. In this structure the water molecules form spiraling channels with a free bore hole along the z axis of ∼5.26 Å and with a diameter of 6.10 Å, which can accommodate the H2 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 C0 phase might be proton ordered. This is at variance with the C1-structure filled ice, which exhibits proton order and, consequently, a lattice phonon spectrum with evident sharp lines17. The rotational spectrum of the hydrogen molecules (Fig. 3c) presents well distinct lines, proving that the H2 molecules rotate almost freely. The S0(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 configuration18, at a distance from water oxygen atoms of ∼3.1 Å and at ∼2.95 Å between each other, the same as in solid hydrogen at ∼1.5 GPa. The larger splitting of the rotational band components is accountable on the basis of the shortest distance between H2 and H2O molecules in the _P_3112 structure, and indirectly confirms it. In addition, at such a short distance the H2–H2 anisotropic interaction may sustain collective rotational excitations, as it happens in solid H2 (ref. 23) and Ar(H2)2 high-pressure compound24,25, which may contribute to the width of the rotational lines. The spectrum of the H2 vibron region (Fig. 3d) consists of only one doublet (Q1(0) and Q1(1) lines), and this is in accordance with the occupation of a single crystallographic site for the hydrogen molecules, as in the proposed structure18. The OH stretching band (Fig. 3b) has the main peak at 3,123 cm−1, 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 compound26 or, for the same structure, as a function of pressure, as reported for ice VII (ref. 27) and C1-phase filled ice28. The higher frequency of the OH stretching mode observed for the C0-structure filled ice is in accordance with the larger O–O distance in the structure _P_3112 (d _OO_≃2.83–2.85 Å ) than in sII clathrate (d _OO_≃2.77–2.79 Å) and in ice Ih (d _OO_≃2.75 Å).

Figure 3: Raman spectra measured at 30 K.

figure 3

Spectra measured for the C0 sample (black lines), compared with those measured in sII clathrate16 (red line) and ice Ih (blue line). The lattice phonon band in the range 120–310 cm−1 (a) is at lower frequencies, whereas 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, S0(0) and S0(1) of the hydrogen molecule, (c) are broader than in sII clathrates, indicating a stronger interaction with the host. A single pair of lines (Q1(0) and Q1(1)) observed for the H2 vibron attests the existence of a single crystallographic site for the H2 molecule in the C0 crystal, opposite to the sII clathrate, where the different peaks correspond to different environments for the H2 molecules (d).

Full size image

Raman band intensities and hydrogen content

We derive the hydrogen molar fraction X from the intensity ratio of the hydrogen rotational lines and the lattice phonon band, _I_rot/_I_phon. 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, as the difference of polarizability of the same molecules in different environments may lead to tiny differences in intensity. This matter has been considered in section IIIa of ref. 16 and we refer to that study 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 cages16. The spectra of these samples in the 200–900 cm−1 range (not reported in ref. 16) are used to calculate _I_rot/_I_phon. To these data we add the information we derive from new volumetric measurements performed on our C0-phase 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 _X_≃25%.

Emptying of the C0-structure filled ice

We discover that the C0 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 h at a temperature of ∼120 K. Considering the sensitivity and the signal-to-noise ratio of our detection system, we can ascertain that the H2/H2O molar fraction in the emptied sample is <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 ∼7 cm−1 for the lattice phonon bands (Fig. 4a) and a decrease of ∼20 cm−1 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 C0 ice and that of ice XVII are essentially the same. The situation for the C0 structure is different from previous observation in sII, where on cages emptying, the lattice constant and, consequently, the O–O distance increases3. The contraction of the water framework on 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 on emptying29. 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 C0 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 authors30, to explain the non-perfect refinement of the X-ray diffraction pattern with the _P_3112 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 N2 stretching mode, observing, at 20 K before the annealing, a quite evident and sharp line at a frequency 2323.7±0.3 cm−1. By virtue of their vibration frequency, sensibly lower than both that of N2 gas (2329.917 cm−1)31, and that of solid N2 in the _α_-phase (a doublet at 2327.5 and 2328.5 cm−1)32, we can establish that the N2 molecules are trapped inside the structure. The Raman line of nitrogen molecules essentially disappears (it becomes 1,000 times weaker) after the annealing process. The evident spectroscopic signatures of structural modification described above are a combined effects of the removal of both H2 and N2 from the sample. The structure of deuterated ice XVII is determined by means of neutron diffraction measurements performed on OSIRIS at ISIS, RAL (UK)[33](/articles/ncomms13394#ref-CR33 "del Rosso, L. et al. Refined structure of metastable ice XVII from neutron diffraction measurements. J. Phys. Chem. Lett Preprint at https://arxiv.org/abs/1609.04996

               (2016)."). Reitveld refinement of the data enables us to determine the structure. The empty structure is described by the space group P6122, with oxygen and deuterium atoms in 6_b_ and 12_c_ positions, respectively. The space group used for the refinement has a higher symmetry than that proposed for the C0\-phase filled ice (ref. [20](/articles/ncomms13394#ref-CR20 "Smirnov, G. S. & Stegailov, V. V. Toward determination of the new hydrogen hydrate clathrate structures. J. Phys. Chem. Lett. 4, 3560–3564 (2013).")) on the basis of the data in ref. [18](/articles/ncomms13394#ref-CR18 "Efimchenko, V. S. et al. New phase in the water-hydrogen system. J. Alloys Comp. 509, S860–S863 (2011)."). This may be due to the presence in the channels of either H2 guests and N2 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. [33](/articles/ncomms13394#ref-CR33 "del Rosso, L. et al. Refined structure of metastable ice XVII from neutron diffraction measurements. J. Phys. Chem. Lett Preprint at 
                https://arxiv.org/abs/1609.04996
                
               (2016).").

Figure 4: Effect of heating and hydrogen release.

figure 4

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 C0 ice, that is, 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,b). Red lines in c show the spectra of the hydrogen molecules for the maximum filling obtained in this work (_X_≃40%).

Full size image

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 that is pressure dependent, but may grow larger than that initially present just after synthesis. As a matter of fact, we generally obtain _X_=25% 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 it into 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 C0 sample. The shape of the rotational lines does not change significantly, demonstrating that H2 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 H2/H2O molar fraction. We notice that, in some instances, as for the most intense spectrum shown in Fig. 4c, the measured H2/H2O molar ratio in the sample exceeds 40%, reaching higher values than the established theoretical maximum (48/136≃35%) for sII hydrates.

We measure several adsorption isotherms, deriving the amount of adsorbed H2 as a function of pressure, at several different temperatures between 15 and 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 Δ_t_ from pressure increase, which is indicated in each panel. 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 Δ_t_=10 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 ∼240 mbar, after a 16.5 h-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 with sII clathrates). This is not surprising at all, as the channels in ice XVII are much wider than the free space inside both the six-membered and five-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 Van’t Hoff equation, the enthalpy of adsorption, Δ_h_, for a given hydrogen uptake X:

Figure 5: Adsorption isotherms measured at different temperatures.

figure 5

We report the hydrogen molar fraction X measured while increasing (square symbols) and decreasing (triangular symbols) hydrogen pressure in annealed ice XVII. Owing to calibration, X 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, whereas 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.

Full size image

where R is the gas constant34,35. 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 −Δ_h_ decrease with increasing molar ratio X starting from ∼5 kJ mol−1 at _X_≃10% down to ∼2 kJ mol−1 at _X_≃40%, in agreement with typical values for physisorption.

Discussion

This work is the first characterization of the dynamics of the C0-phase filled ice. Moreover, two remarkable results arise from this study. First, we discover that the C0 structure is metastable when emptied. Owing 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 _X_=50%), 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, which make ice XVII really appealing for hydrogen storage applications.

Methods

Procedure for the synthesis of the sample

We produce several instances of solid H2–H2O 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 C0-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 ∼255 K (−18 °C). We prudentially leave the sample under pressure for a few days, even though, while rising pressure during synthesis, we observe at ∼360 MPa an abrupt decrease of the rate of pressure rise, which is a clear indication of large hydrogen adsorption by the ice. Next, 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-ray diffraction, except that in this case a drop of water is frozen from the beginning inside a few common X-ray glass capillaries, having 0.5 mm diameter and a length of ∼10–15 mm. The frozen water occupies a capillary length of ∼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 cryostat22. 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 ≃10 K. Gas pressure in the cell can be varied and measured with an accuracy of ∼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 ∼30 mW on the sample, focused on a spot of about 30 μm. Scattered light is collected in an almost backscattering geometry, focused on the entrance slit of a Spex spectrograph and recorded by a cooled CCD (charge-coupled device) detector (Andor). The maximum resolution of the instrument is 0.4 cm−1. 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 (∼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 min. 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 S0(0) and S0(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 _σ_0 and _σ_1 the Raman cross-section for the S0(0) and S0(1) rotational lines, and with _I_0 and _I_1 their intensity, respectively, the quantity _I_rot=_I_0/_σ_0+_I_1/_σ_1 is proportional to the number of H2 molecules in the lowest two _J_=0 and _J_=1 rotational states, that is, at low temperature, to the number of H2 molecules giving Raman scattering. Dividing _I_rot by the intensity of the water lattice phonon band, _I_phon, we obtain a quantity which is independent of laser intensity and detector efficiency, and proportional to the H2/H2O molar ratio in the sample. This procedure needs a calibration, which we have accomplished fitting data obtained with two independent methods. For several samples of sII clathrates studied in our laboratory, the hydrogen content has been indirectly derived, by counting the molecules in the small and large cages16, whereas the spectra of the same samples in the 200–900 cm−1 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 C0 sample, while measuring, with Raman scattering, the adsorption isotherms at three temperatures, namely 20, 50 and 80 K. These results are reported as coloured 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 0.233±0.023, which is the value we have used in this work.

Figure 6: Raman intensity calibration.

figure 6

Summary of the determinations of the H2/H2O molar ratio plotted as a function of the concomitant measurements of the _I_rot/_I_phon ratio, used to determine the calibration factor. The solid black squares represent measures of several instances of sII pure H2 clathrate hydrates, synthesized at different pressures and/or undergone to different thermal cycles16, whereas 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.

Full size image

Data availability

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

Additional information

How to cite this article: del Rosso, L. et al. New porous water ice metastable at atmospheric pressure obtained by emptying a hydrogen-filled ice. Nat. Commun. 7, 13394 doi: 10.1038/ncomms13394 (2016).

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

References

  1. Bartels-Rausch, T. et al. Ice structures, patterns, and processes: A view across the icefields. Rev. Mod. Phys. 84, 885–944 (2012).
    Article CAS ADS Google Scholar
  2. Sloan, E. D. & Koh, C. A. Clathrate Hydrates of Natural Gases, 3rd Edition Taylor & Francis (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).
    Article CAS ADS Google Scholar
  4. Huang, Y. et al. A new phase diagram of water under negative pressure: The rise of the lowest-density clathrate s-III. Sci. Adv. 2, e1501010 (2016).
    Article ADS Google Scholar
  5. Dyadin, Y. A. et al. Clathrate hydrates of hydrogen and neon. Mendeleev Commun. 9, 209–210 (1999).
    Article Google Scholar
  6. Mao, W. L. et al. Hydrogen clusters in clathrate hydrate. Science 297, 2247–2249 (2002).
    Article CAS ADS Google Scholar
  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).
    Article ADS Google Scholar
  8. Ulivi, L. et al. Quantum rattling of molecular hydrogen in clathrate hydrate nanocavities. Phys. Rev. B 76, 161401 (2007).
    Article ADS Google Scholar
  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).
    Article ADS Google Scholar
  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–8 (2013).
    Article CAS ADS Google Scholar
  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).
    Article ADS Google Scholar
  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).
    Article CAS Google Scholar
  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).
    Article ADS Google Scholar
  14. Strobel, T. A., Sloan, E. D. & Koh, C. A. Raman spectroscopic studies of hydrogen clathrate hydrates. J. Phys. Chem. 130, 014506 (2009).
    Article Google Scholar
  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).
    Article ADS Google Scholar
  16. Zaghloul, M. A., Celli, M., Salem, N. M., Elsheikh, S. M. & Ulivi, L. High pressure synthesis and in situ Raman spectroscopy of H2 and HD clathrate hydrates. J. Chem. Phys. 137, 164320 (2012).
    Article ADS Google Scholar
  17. Vos, W. L., Finger, L. W., Hemley, R. J. & Mao, H. K. Novel H2+H2O clathrates at high pressures. Phys. Rev. Lett. 71, 3150–3153 (1993).
    Article CAS ADS Google Scholar
  18. Efimchenko, V. S. et al. New phase in the water-hydrogen system. J. Alloys Comp. 509, S860–S863 (2011).
    Article CAS Google Scholar
  19. Strobel, T. A., Somayazulu, M. & Hemley, R. J. Phase behavior of H2+H2O at high pressures and low temperatures. J. Phys. Chem. C 115, 4898–4903 (2011).
    Article CAS Google Scholar
  20. Smirnov, G. S. & Stegailov, V. V. Toward determination of the new hydrogen hydrate clathrate structures. J. Phys. Chem. Lett. 4, 3560–3564 (2013).
    Article CAS Google Scholar
  21. Qian, G. R., Lyakhov, A. O., Zhu, Q., Oganov, A. R. & Dong, X. Novel hydrogen hydrate structures under pressure. Sci. Rep. 4, 5606 (2014).
    Article CAS Google Scholar
  22. 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).
    Article CAS ADS Google Scholar
  23. Van Kranendonk, J. Solid Hydrogen Plenum Press (1983).
  24. Grazzi, F., Santoro, M., Moraldi, M. & Ulivi, L. Anisotropic interactions of hydrogen molecules from the pressure dependence of the rotational spectrum in the Ar(H2)2 compound. Phys. Rev. Lett. 87, 125506 (2001).
    Article CAS ADS Google Scholar
  25. Grazzi, F., Santoro, M., Moraldi, M. & Ulivi, L. Roton excitations of the hydrogen molecule in the Ar(H2)2 compound. Phys. Rev. B 66, 144303 (2002).
    Article ADS Google Scholar
  26. Pimentel, G. C. & Sederholm, C. H. Correlation of infrared stretching frequencies and hydrogen bond distances in crystals. J. Chem. Phys. 24, 639–641 (1956).
    Article CAS ADS Google Scholar
  27. Pruzan, P. Pressure effects on the hydrogen bond in ice up to 80 GPa. J. Mol. Struct. 322, 279–286 (1994).
    Article CAS ADS Google Scholar
  28. Vos, W. L., Finger, L. W., Hemley, R. J. & Mao, H. K. Pressure dependence of hydrogen bonding in a novel H2+H2O clathrate. Chem. Phys. Lett. 257, 524–530 (1996).
    Article CAS ADS Google Scholar
  29. 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).
    Article CAS ADS Google Scholar
  30. Smirnov, G. S. & Stegailov, V. V. Anomalous diffusion of guest molecules in hydrogen gas hydrates. High Temp. 53, 829–836 (2015).
    Article CAS Google Scholar
  31. Bendtsen, J. The rotational and rotation-vibrational Raman spectra of 14N2, 14N15N and 15N2 . J. Raman Spectrosc. 2, 133–145 (1974).
    Article CAS ADS Google Scholar
  32. Ouillon, R., Turc, C., Lemaistre, J. P. & Ranson, P. High resolution Raman spectroscopy in the α and β crystalline phases of N2 . J. Chem. Phys. 93, 3005–3011 (1990).
    Article CAS ADS Google Scholar
  33. del Rosso, L. et al. Refined structure of metastable ice XVII from neutron diffraction measurements. J. Phys. Chem. Lett Preprint at https://arxiv.org/abs/1609.04996 (2016).
  34. Do, D. D. Adsorption Analysis: Equilibria and Kinetics Imperial College Press (1998).
  35. Schlichtenmayer, M. & Hirscher, M. The usable capacity of porous materials for hydrogen storage. Appl. Phys. A 122, 379 (2016).
    Article ADS Google Scholar
  36. Lokshin, K. A. & Zhao, Y. Fast synthesis method and phase diagram of hydrogen clathrate hydrate. App. Phys. Lett. 88, 131909 (2006).
    Article ADS Google Scholar
  37. Efimchenko, L. et al. 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).
    Article CAS ADS Google Scholar
  38. 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).
    Article CAS Google Scholar
  39. 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).
    Article ADS Google Scholar

Download references

Acknowledgements

We 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.

Author information

Authors and Affiliations

  1. Consiglio Nazionale delle Ricerche, Istituto dei Sistemi Complessi, via Madonna del Piano 10, Sesto Fiorentino, I-50019, Italy
    Leonardo del Rosso, Milva Celli & Lorenzo Ulivi
  2. Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, via Sansone 1, Sesto Fiorentino, I-50019, Italy
    Leonardo del Rosso

Authors

  1. Leonardo del Rosso
    You can also search for this author inPubMed Google Scholar
  2. Milva Celli
    You can also search for this author inPubMed Google Scholar
  3. Lorenzo Ulivi
    You can also search for this author inPubMed Google Scholar

Contributions

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.

Corresponding author

Correspondence toLorenzo Ulivi.

Ethics declarations

Competing interests

The authors declare no competing financial interests.

Supplementary information

Rights and permissions

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Reprints and permissions

About this article

Cite this article

del Rosso, L., Celli, M. & Ulivi, L. New porous water ice metastable at atmospheric pressure obtained by emptying a hydrogen-filled ice.Nat Commun 7, 13394 (2016). https://doi.org/10.1038/ncomms13394

Download citation