REVIEW 2 major objections 5 minor 44 references
Hydrogen-rich hydrate at high pressures up to 104 GPa
T0 review · 2 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read A water–hydrogen ice doubles its molecular H2 content at high pressure.
desk verdict A solid experimental determination of the C3 filled-ice structure and its room-temperature formation path, with the composition resting on a fixed H2 occupancy that deserves a free-refinement check. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing identity is the crystallographic site switch between the two filled ices: in the cubic space group Fd-3m, C2 places molecular hydrogen on the 8b sites while C3 places it on the 16c sites, which are twice as numerous, so occupancy of the same water network doubles the H2 content to (H2O)(H2)2. The argument is carried by single-crystal X-ray diffraction on multiply twinned crystals, which returned 41 independent reflections and structure solution from difference Fourier maps (R1 = 4.06%), together with Raman measurements of the H2 intramolecular vibron, whose downward shift marks the denser hydrogen packing. First-principles calculations using hydrogen-ordered Pna21 and P41 proxy structures supply the enthalpy ordering and the theoretical volume and vibron curves that the experiment is compared against. The cubic experimental C3 structure is interpreted as a disordered average of the lower-symmetry tetragonal ordered form, analogous to ice VII versus ice VIII.
What would settle it
Measure the single-crystal structure of the cold-compressed sample near 100 GPa and refine the hydrogen positions directly; if the H2 molecules do not sit at 16c with the (H2O)(H2)2 composition, or if the C2 reflections persist at higher pressures than reported, the central claim fails.
Extended reading notes
Core claim
At 69 GPa, laser heating of the C2 filled ice produced single crystals of a new phase whose structure the paper solves completely by single-crystal X-ray diffraction. The C3 phase is cubic, space group Fd-3m, with oxygen atoms on the 8b sites of the same hydrogen-bonded water network as C2; the difference is where molecular hydrogen sits. In C2 the H2 molecules occupy the 8b positions (one H2 per formula unit), while in C3 they occupy the 16c positions, twice as many sites, refined as a single scatterer of occupancy 2, yielding (H2O)(H2)2. Refinement gives R1 = 4.06% from 41 independent reflections, with hydrogen positions taken directly from difference Fourier maps. The paper further claims that cold compression under hydrogen-rich conditions converts C2 into the same C3 phase gradually over 47–104 GPa, with the C2 diffraction peaks disappearing by 103 GPa, and that laser heating above about 1500 K makes the transition abrupt at pressures as low as 47 GPa. Density functional calculations reproduce the volume relations and the softening of the H2 vibron, and place thermodynamic stability of C3 above about 23 GPa.
Load-bearing premise
The identification of the cold-compressed material at 103–104 GPa with the laser-heated C3 structure solved at 69 GPa rests on matching Fd-3m diffraction peak positions, volume extrapolation, and vibron convergence rather than on a full structure solution at the highest pressure, and the laser-heating temperatures were estimated visually rather than measured.
Editorial extensions
If this is right
- C3 doubles the molecular-hydrogen capacity of the filled ice relative to C2, so hydrogen storage in hydrate form can hold more H2 at a given pressure.
- The C2-to-C3 transformation at 300 K proceeds slowly and only reaches completion after weeks at 103 GPa, implying a kinetically limited reaction that can be completed by laser heating above about 1500 K.
- Once formed, C3 survives decompression to 40 GPa, which defines a practical window for recovering hydrogen-rich hydrate to moderate pressures.
- Hydrogen-bond symmetrization in the water sublattice is delayed from about 32–40 GPa in C2 to about 78–87 GPa in C3, showing that extra H2 stiffens the water framework against symmetrization.
- The band gap of C3 narrows rapidly with pressure up to about 100 GPa and then levels off, keeping the hydrate insulating at the highest pressures studied.
Reading between the lines
- If C3 is representative of deep ice mantles, planetary models may need to include a hydrogen reservoir that is twice as concentrated as previously assumed at the same depths; the paper's 104 GPa stability supports that possibility but does not itself model planetary interiors.
- The cubic Fd-3m symmetry may be a time and space average of ordered tetragonal domains; low-temperature experiments on the same composition could reveal an ordered analog analogous to ice VIII.
- The gradual pressure-driven uptake suggests that intermediate pressures below 47 GPa might also convert C2 to C3 given longer times or higher H2 fugacity, a testable kinetic prediction.
- Comparing the 1.94 Å H2–H2 distance at 100 GPa with bulk hydrogen suggests the hydrate cavities keep H2 near the onset of metallization physics, and the saturating band-gap trend could be probed by optical conductivity measurements.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports a combined synchrotron single-crystal X-ray diffraction, Raman spectroscopy, and first-principles DFT study of the hydrogen hydrate system (H2O)(H2)n up to 104 GPa. The authors find that the C2 filled ice transforms to the C3 filled ice at room temperature gradually over 47–104 GPa in hydrogen-rich conditions, and that laser heating induces the transformation at lower pressures. At 69 GPa, a single-crystal structure solution assigns the C3 phase to cubic Fd3m, with H2 molecules in the 16c positions and a nominal composition (H2O)(H2)2, twice the H2 content of C2. The paper also reports Raman vibron softening across the transition, theoretical stability calculations placing C2 stable to 14 GPa and C3 stable above 23 GPa, and band-gap calculations for both phases.
Significance. If the structural assignment and composition hold, this is a significant result: it provides the first direct crystallographic evidence for the structure of the most hydrogen-dense filled ice and updates the pressure-composition stability landscape relevant to giant-planet interiors and hydrogen storage. The SC XRD work is technically impressive: 41 independent reflections, a complete structure solution with R1 = 4.06%, twinning handled by simultaneous integration, and the structure deposited in the CSD. The DFT calculations are largely independent of the experimental identification and support the phase stability picture; the hybrid-functional band-gap results are not load-bearing. The main caveat, discussed below, is that the doubled H2 content relies on a fixed, unrefined site occupancy.
major comments (2)
- [Table 1 and Section III (SC XRD refinement at 69 GPa)] The central compositional claim, (H2O)(H2)2, is stated as determined by single-crystal XRD, but the H2 molecule is modeled as a single H atom on the 16c site with the occupancy fixed at 2, and no free refinement or estimated standard deviation for this occupancy is reported. Because hydrogen is a weak X-ray scatterer and the refinement uses only 41 independent reflections, the diffraction data alone do not independently establish the full 16c occupancy. The agreement of the measured volume per formula unit with the theoretical P41 C3 volume is supporting evidence, but it is not equivalent to a diffraction determination. The authors should either report a free refinement of the 16c occupancy (with its esd) or explicitly state that the composition is inferred from the site assignment plus volume/Raman evidence rather than determined by refinement.
- [Section III, cold-compression results at 103–104 GPa] The identification of the cold-compressed phase at 103–104 GPa as the C3 phase is based on matching Fd3m peak positions, extrapolated unit-cell volumes, and Raman vibron convergence with the 69 GPa laser-heated structure, not on a full structure solution at the highest pressure. The text should state this explicitly and discuss the possibility that the cold-compressed material may be compositionally intermediate, since the room-temperature transformation is described as gradual and the H2 uptake may be incomplete. This is load-bearing for the abstract claim that the C2-to-C3 transformation occurs at room temperature up to 104 GPa.
minor comments (5)
- [Section II.A (Methods)] Please correct 'spectradiometric' to 'spectroradiometric' and remove the doubled 'by' in 'by via absorption'.
- [Figures 5 and 7] Pressure, volume, and Raman frequency data are presented without error bars or quantitative uncertainties; adding error estimates from the Au/ruby pressure markers and peak fits would strengthen the comparison with theory.
- [Section II.A (Laser heating)] The maximum temperature of approximately 1500 K is stated as a visual estimate; since the lower-pressure laser-heating pathway to C3 depends on this temperature, the authors should clearly mark it as approximate and give an uncertainty range.
- [Reference 12] The Supplemental Material reference contains 'Supplemental Figures S1-S10' twice in a single sentence; the duplicate should be removed.
- [Figure 6 caption] The caption lists an excitation wavelength of 630 nm, while the text reports 660 nm; please make these consistent.
Circularity Check
No significant circularity: the experimental structure solution and C2→C3 transition are independent of the DFT calculations; the fixed H2 occupancy is a refinement limitation rather than a circular input.
full rationale
The paper's central claims—the room-temperature C2→C3 transformation, the cubic Fd3̅m structure of C3, and the doubling of the H2-site multiplicity (8b→16c)—rest on synchrotron single-crystal XRD, which is an independent experimental measurement, together with Raman and volume-pressure data. The DFT calculations (PBE/Grimme and hybrid BLYP) are first-principles; the BLYP mixing parameters are fit to experimental energies, but the band-gap results are not load-bearing for the phase-transition or composition claims. The self-citations are to the authors' own method papers (e.g., Refs. [15]–[17]) and do not supply any load-bearing premise. One caveat is that the absolute composition (H2O)(H2)2 in Table 1 uses H2 modeled as a single H atom with occupancy fixed at 2 rather than freely refined; this means the absolute stoichiometry is partly assumed, not independently determined from the 41-reflection refinement. This is a crystallographic/correctness limitation, not a circularity: the 16c position and its multiplicity are experimental outputs, and the composition is corroborated by the volume/EOS match to the predicted P41 C3 phase and by prior work. The paper does not define its inputs in terms of its conclusions, so no circular step is established.
Assumptions & free parameters
free parameters (1)
- Hybrid BLYP exchange-correlation mixing parameters =
not specified
assumptions (3)
- domain assumption PBE+DFT-D with ordered Pna21 and P41 structures faithfully represents the energetics of the disordered experimental C2 and C3 phases.
- domain assumption The C3 phase synthesized by laser heating at 69 GPa is the same phase as that obtained by cold compression at 104 GPa.
- domain assumption H2 molecules can be represented as a single atom at the 16c site with occupancy 2 due to rotational disorder.
Cite this review
Pith. "Pith review of Hydrogen-rich hydrate at high pressures up to 104 GPa." pith.science (2026). https://pith.science/paper/HZ43237V
@misc{pith2026250507091,
author = {Pith},
title = {Pith review of: Hydrogen-rich hydrate at high pressures up to 104 GPa},
year = {2026},
howpublished = {\url{https://pith.science/paper/HZ43237V}},
note = {Machine review of arXiv:2505.07091}
}
read the original abstract
Gas hydrates are considered fundamental building blocks of giant icy planets like Neptune and similar exoplanets. The existence of these materials in the interiors of giant icy planets, which are subject to high pressures and temperatures, depends on their stability relative to their constituent components. In this study, we reexamine the structural stability and hydrogen content of hydrogen hydrates, (H2O)(H2)n, up to 104 GPa, focusing on hydrogen-rich materials. Using synchrotron single-crystal X-ray diffraction, Raman spectroscopy, and first-principles theoretical calculations, we find that the C2-filled ice phase undergoes a transformation to C3-filled ice phase over a broad pressure range of 47 - 104 GPa at room temperature. The C3 phase contains twice as much molecular H2 as the C2 phase. Heating the C2-filled ice above approximately 1500 K induces the transition to the C3 phase at pressures as low as 47 GPa. Upon decompression, this phase remains metastable down to 40 GPa. These findings establish new stability limits for hydrates, with implications for hydrogen storage and the interiors of planetary bodies.
Reference graph
Works this paper leans on
-
[1]
We have demonstrated that the C2 to C3 transition can be achieved through cold compression under H₂-rich conditions at pressures as high as 103 GPa. This transition occurs gradually and is marked by both structural and vibrational changes, offering a clearer understanding of the pressure- induced phase behavior of these hydrogen-water inclusion compounds
-
[2]
Single-crystal XRD analysis reveals that the C 3 phase crystallizes in the cubic space group Fd3̅m, which is isostructural with the C2 phase with respect to the positions of the oxygen atoms. The key difference lies in the arrangement of hydrogen molecules: in the C 3 phase, H₂ molecules occupy the 16c positions, whereas in the C 2 phase, they are positio...
-
[3]
Raman spectroscopy across the C2 to C3 transition reveals a notable downward shift in the main vibron mode frequency. This frequency shift supports the structural changes identified via XRD, particularly the increased H₂ density in the C3 phase. The observed broadening and changes in the spectral features confirm intensified molecular interaction during t...
-
[4]
Our theoretical calculations of the electronic structure of the C2 and C3 phases under high pressure reveal a pressure induced narrowing of the bandgap of the C 3 phase primarily driven by enhanced intermolecular coupling between H2 molecules. However, this effect appears to level off at pressures near 100 GPa, indicating that the system reaches a saturat...
-
[5]
T.A. Strobel, M. Somayazulu, S.V . Sinogeikin, P. Dera, R.J. Hemley, Hydrogen-Stuffed, Quartz-like Water Ice, Journal of the American Chemical Society, 138 (2016) 13786-13789
work page 2016
-
[6]
Y . Wang, K. Glazyrin, V . Roizen, A.R. Oganov, I. Chernyshov, X. Zhang, E. Greenberg, V .B. Prakapenka, X. Yang, S.-q. Jiang, A.F. Goncharov, Novel Hydrogen Clathrate Hydrate, Physical Review Letters, 125 (2020) 255702
work page 2020
- [7]
-
[8]
T.A. Strobel, M. Somayazulu, R.J. Hemley, Phase Behavior of H2 + H2O at High Pressures and Low Temperatures, The Journal of Physical Chemistry C, 115 (2011) 4898-4903
work page 2011
Show all 44 references
-
[9]
Efimchenko, M.A
V .S. Efimchenko, M.A. Kuzovnikov, V .K. Fedotov, M.K. Sakharov, S.V . Simonov, M. Tkacz, New phase in the water–hydrogen system, Journal of Alloys and Compounds, 509 (2011) S860- S863
2011
-
[10]
V os, L.W
W.L. V os, L.W. Finger, R.J. Hemley, H.-k. Mao, Novel H2-H2O clathrates at high pressures, Physical Review Letters, 71 (1993) 3150-3153
1993
-
[11]
Hirai, S
H. Hirai, S. Kagawa, T. Tanaka, T. Matsuoka, T. Yagi, Y . Ohishi, S. Nakano, Y . Yamamoto, T. Irifune, Structural changes of filled ice Ic hydrogen hydrate under low temperatures and high pressures from 5 to 50 GPa, The Journal of Chemical Physics, 137 (2012)
2012
-
[12]
V os, L.W
W.L. V os, L.W. Finger, R.J. Hemley, H.-k. Mao, Pressure dependence of hydrogen bonding in a novel H2O-H2 clathrate, Chemical Physics Letters, 257 (1996) 524-530
1996
-
[13]
Machida, H
S.-i. Machida, H. Hirai, T. Kawamura, Y . Yamamoto, T. Yagi, Structural changes of filled ice Ic structure for hydrogen hydrate under high pressure, The Journal of Chemical Physics, 129 (2008) 224505
2008
-
[14]
Qian, A.O
G.-R. Qian, A.O. Lyakhov, Q. Zhu, A.R. Oganov, X. Dong, Novel Hydrogen Hydrate Structures under Pressure, Scientific Reports, 4 (2014) 5606
2014
-
[15]
Ranieri, S
U. Ranieri, S. Di Cataldo, M. Rescigno, L. Monacelli, R. Gaal, M. Santoro, L. Andriambariarijaona, P. Parisiades, C. De Michele, L.E. Bove, Observation of the most H2-dense filled ice under high pressure, Proceedings of the National Academy of Sciences, 120 (2023) e2312665120
2023
-
[16]
See Supplemental Material at http://link.aps.org/supplemental/ 10.1103/PhysRevB
T.S.-S. See Supplemental Material at http://link.aps.org/supplemental/ 10.1103/PhysRevB. xxx This pdf file contains Supplemental Figures S1-S10, Supplemental Figures S1-S10, Tables S1-S2
-
[17]
Mezouar, R
M. Mezouar, R. Giampaoli, G. Garbarino, I. Kantor, A. Dewaele, G. Weck, S. Boccato, V . Svitlyk, A.D. Rosa, R. Torchio, O. Mathon, O. Hignette, S. Bauchau, Methodology for in situ synchrotron X-ray studies in the laser-heated diamond anvil cell, High Pressure Research, 37 (201...
2017
-
[18]
Holtgrewe, E
N. Holtgrewe, E. Greenberg, C. Prescher, V .B. Prakapenka, A.F. Goncharov, Advanced integrated optical spectroscopy system for diamond anvil cell studies at GSECARS, High Pressure Research, 39 (2019) 457-470
2019
-
[19]
Goncharov, N
A.F. Goncharov, N. Holtgrewe, G. Qian, C. Hu, A.R. Oganov, M. Somayazulu, E. Stavrou, C.J. Pickard, A. Berlie, F. Yen, M. Mahmood, S.S. Lobanov, Z. Konôpková, V .B. Prakapenka, Backbone NxH compounds at high pressures, The Journal of Chemical Physics, 142 (2015) 214308
2015
-
[20]
Bykova, I.G
E. Bykova, I.G. Batyrev, M. Bykov, E. Edmund, S. Chariton, V .B. Prakapenka, A.F. Goncharov, Structural evolution of iodine on approach to the monatomic state, Physical Review B, 108 (2023) 024104. 23
2023
-
[21]
H. Chen, M. Bykov, I.G. Batyrev, L. Brüning, E. Bykova, M.F. Mahmood, S. Chariton, V .B. Prakapenka, T. Fedotenko, H.-P. Liermann, K. Glazyrin, A. Steele, A.F. Goncharov, High- pressure Synthesis of Cobalt Polynitrides: Unveiling Intriguing Crystal Structures and Nitridation B...
2024 doi
-
[22]
CrysAlisPro Software System (Rigaku Oxford Diffraction) (Oxford, Oxford, UK, 2014)
2014
-
[23]
Aslandukov, M
A. Aslandukov, M. Aslandukov, N. Dubrovinskaia, L. Dubrovinsky, Domain Auto Finder (DAFi) program: the analysis of single-crystal X-ray diffraction data from polycrystalline samples, Journal of Applied Crystallography, 55 (2022) 1383-1391
2022
-
[24]
Hohenberg, W
P. Hohenberg, W. Kohn, Inhomogeneous Electron Gas, Physical Review, 136 (1964) B864- B871
1964
-
[25]
Kohn, L.J
W. Kohn, L.J. Sham, Self-Consistent Equations Including Exchange and Correlation Effects, Physical Review, 140 (1965) A1133-A1138
1965
-
[26]
Perdew, K
J.P. Perdew, K. Burke, M. Ernzerhof, Generalized Gradient Approximation Made Simple, Physical Review Letters, 77 (1996) 3865-3868
1996
-
[27]
Segall, J.D.L
M.D. Segall, J.D.L. Philip, M.J. Probert, C.J. Pickard, P.J. Hasnip, S.J. Clark, M.C. Payne, First-principles simulation: ideas, illustrations and the CASTEP code, Journal of Physics: Condensed Matter, 14 (2002) 2717
2002
-
[28]
Becke, A new mixing of Hartree–Fock and local density‐functional theories, The Journal of Chemical Physics, 98 (1993) 1372-1377
A.D. Becke, A new mixing of Hartree–Fock and local density‐functional theories, The Journal of Chemical Physics, 98 (1993) 1372-1377
1993
-
[29]
Grimme, Semiempirical GGA-type density functional constructed with a long-range dispersion correction, Journal of Computational Chemistry, 27 (2006) 1787-1799
S. Grimme, Semiempirical GGA-type density functional constructed with a long-range dispersion correction, Journal of Computational Chemistry, 27 (2006) 1787-1799
2006
-
[30]
Refson, P.R
K. Refson, P.R. Tulip, S.J. Clark, Variational density-functional perturbation theory for dielectrics and lattice dynamics, Physical Review B, 73 (2006) 155114
2006
-
[31]
Sheldrick, SHELXT - Integrated space-group and crystal-structure determination, Acta Crystallographica Section A, 71 (2015) 3-8
G. Sheldrick, SHELXT - Integrated space-group and crystal-structure determination, Acta Crystallographica Section A, 71 (2015) 3-8
2015
-
[32]
Dolomanov, L.J
O.V . Dolomanov, L.J. Bourhis, R.J. Gildea, J.A.K. Howard, H. Puschmann, OLEX2: a complete structure solution, refinement and analysis program, Journal of Applied Crystallography, 42 (2009) 339-341
2009
-
[33]
Sheldrick, Crystal structure refinement with SHELXL, Acta Crystallographica Section C, 71 (2015) 3-8
G. Sheldrick, Crystal structure refinement with SHELXL, Acta Crystallographica Section C, 71 (2015) 3-8
2015
-
[34]
www.ccdc.cam.ac.uk/structures
-
[35]
Loubeyre, R
P. Loubeyre, R. LeToullec, D. Hausermann, M. Hanfland, R.J. Hemley, H.K. Mao, L.W. Finger, X-ray diffraction and equation of state of hydrogen at megabar pressures, Nature, 383 (1996) 702-704
1996
-
[36]
Goncharov, V .V
A.F. Goncharov, V .V . Struzhkin, H.K. Mao, R.J. Hemley, Raman spectroscopy of dense H2O and the transition to symmetric hydrogen bonds, Physical Review Letters, 83 (1999) 1998- 2001
1999
-
[37]
Silvera, The solid molecular hydrogens in the condensed phase: Fundamentals and static properties, Reviews of Modern Physics, 52 (1980) 393-452
I.F. Silvera, The solid molecular hydrogens in the condensed phase: Fundamentals and static properties, Reviews of Modern Physics, 52 (1980) 393-452
1980
-
[38]
B. Li, Y . Ding, D.Y . Kim, L. Wang, T.-C. Weng, W. Yang, Z. Yu, C. Ji, J. Wang, J. Shu, J. Chen, K. Yang, Y . Xiao, P. Chow, G. Shen, W.L. Mao, H.-K. Mao, Probing the Electronic Band Gap of Solid Hydrogen by Inelastic X-Ray Scattering up to 90 GPa, Physical Review Letters, 12...
2021
-
[39]
Mao, R.J
H.-k. Mao, R.J. Hemley, Ultrahigh-pressure transitions in solid hydrogen, Reviews of Modern Physics, 66 (1994) 671-692. 24
1994
-
[40]
Kranendonk, Solid Hydrogen; Theory of the Properties of Solid H, HD and D
J.v. Kranendonk, Solid Hydrogen; Theory of the Properties of Solid H, HD and D. Plenum Press New York and London 1982. , 1983
1982
-
[41]
Machida, H
S.-i. Machida, H. Hirai, T. Kawamura, Y . Yamamoto, T. Yagi, Raman spectra for hydrogen hydrate under high pressure: Intermolecular interactions in filled ice Ic structure, Journal of Physics and Chemistry of Solids, 71 (2010) 1324-1328
2010
-
[42]
Benoit, A.H
M. Benoit, A.H. Romero, D. Marx, Reassigning Hydrogen-Bond Centering in Dense Ice, Physical Review Letters, 89 (2002) 145501
2002
-
[43]
Meier, S
T. Meier, S. Petitgirard, S. Khandarkhaeva, L. Dubrovinsky, Observation of nuclear quantum effects and hydrogen bond symmetrisation in high pressure ice, Nature Communications, 9 (2018) 2766
2018
-
[44]
Guthrie, R
M. Guthrie, R. Boehler, J.J. Molaison, B. Haberl, A.M. dos Santos, C. Tulk, Structure and disorder in ice VII on the approach to hydrogen-bond symmetrization, Physical Review B, 99 (2019) 184112. 1 Supplementary Information Hydrogen-rich hydrate at high pressures up to 104 GPa...
2019
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.