REVIEW 3 major objections 7 minor 2 references
Formation and Thermodynamic Behavior of THF-Water Hydrates in Confined Mesoporous Media
T0 review · 3 major / 7 minor · reviewed 2026-07-11 · grok-4.5
Pith's one-line read THF clathrate hydrate forms inside mesopores only when the pores are nearly full and the mixture is THF-rich enough, and its melting temperature is fixed by pore curvature, not by cooling rate or thermal history.
desk verdict Solid experimental map of when confined THF hydrate forms in SBA-15; the thermodynamics hold, the in-pore structural proof is thinner than claimed. 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 classical Gibbs–Thomson relation linking melting-point depression to inverse pore radius (ΔTₘ ∝ 1/r), applied simultaneously to confined ice and confined sII hydrate; together with the dual experimental gate of near-percolating filling fraction and stoichiometric THF content that decides whether hydrate can form inside the channels at all.
What would settle it
Prepare a sample at φ ≈ 1.1 and 1:14 THF:water, then measure only the in-pore population (for example by diffraction after removing external liquid or by a contrast method that rejects out-pore hydrate). If the −13.2 °C endotherm vanishes or the Bragg peaks are inconsistent with sII, or if reheating temperatures shift systematically with cooling rate at fixed composition and pore size, the central claim fails.
Extended reading notes
Core claim
Confined sII THF hydrate appears in SBA-15 only when near-percolating filling (φ ≈ 1.0–1.1 cm³/g) and sufficient THF (≥ 1:16 mol:mol) are both satisfied. Reheating peaks for confined ice (−14.7 °C) and confined hydrate (−13.2 °C) remain invariant under changes in cooling rate and thermal cycling, proving that capillarity and composition—not kinetic history—fix the liquidus and dissociation temperatures. Both temperature shifts follow the Gibbs–Thomson law versus inverse pore radius, and the confined product is crystallographically bulk sII.
Load-bearing premise
The endotherm near −13.2 °C is assigned to in-pore structure-II hydrate, and the X-ray patterns from overfilled samples are taken to prove that identity even though those patterns mix signals from hydrate outside the pores with hydrate inside them.
Editorial extensions
If this is right
- Dissociation temperatures of confined hydrates can be predicted from measured pore radius once the Gibbs–Thomson slope is known.
- Cooling protocols can steer the ice-to-hydrate ratio without moving the temperatures at which each phase melts.
- Pore connectivity and guest loading become the primary design switches for phase selection in promoter-assisted storage or separation media.
- The same Gibbs–Thomson framework works for both ice and hydrate, so DSC peak positions can serve as independent reporters of crystalline size.
- Ice-encapsulation anomalies typical of gas hydrates are absent here, allowing cleaner isolation of pure confinement effects.
Reading between the lines
- The dual-threshold requirement implies that continuous liquid pathways and activity buffering by excess external liquid are necessary for hydrate growth along cylindrical channels—a design rule likely to generalize to other guest–water systems in mesopores.
- Because the present X-ray patterns cannot separate in-pore from out-pore hydrate, diffraction on underfilled samples or contrast-matched neutron scattering would be the decisive test of the in-pore lattice alone.
- The roughly 10 % lower solid–liquid interfacial free energy inferred for hydrate versus ice may transfer to other sII formers and refine capillary models used for methane hydrates in sediments.
- Only a few unit cells fit across ~4 nm pores, so hydrate formation is expected to shut off below a characteristic diameter near 2 nm, offering a testable cutoff for storage-media design.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript combines DSC and WAXS to map how composition, pore filling, and cooling rate control THF–water ice versus sII clathrate formation in bulk and in SBA-15 mesopores. Bulk runs provide mass-balanced ice and hydrate enthalpies and reversible melting/dissociation temperatures. In confinement, heating endotherms separate into a Gibbs–Thomson-depressed confined-ice melt (−14.7 ± 0.2 °C) and a higher-temperature peak assigned to in-pore hydrate (−13.2 ± 0.2 °C), with an external hydrate near +4 °C when overfilled. Confined hydrate is reported only for near-percolating filling (φ ≈ 1.0–1.1 cm³/g) and THF content ≥ 1:16. Slower cooling increases the confined-hydrate enthalpy fraction without shifting reheating temperatures; thermal cycling likewise leaves peak positions invariant. ΔT_m versus inverse pore radius is linear for both ice and hydrate, and overfilled WAXS patterns are indexed to cubic sII.
Significance. If the calorimetric assignments hold, the work cleanly separates equilibrium markers (invariant reheating temperatures fixed by capillarity and composition) from kinetic phase selection (rate-dependent peak areas) in a model hydrate system that avoids high-pressure gas-hydrate complications. The pure-water SBA-15 calibration, filling- and composition-window maps, cooling-rate and cycling controls, and quantitative Gibbs–Thomson analysis for both ice and hydrate are genuine strengths and will be useful for interpreting confined hydrates and promoter-assisted storage concepts. The bulk mass-balance table and explicit f_hyd deconvolution further support reproducibility of the thermodynamic claims.
major comments (3)
- §3.6 and Fig. 9: The claim that “the phase formed in pores is crystallographically identical to bulk sII” is not established by the reported WAXS. The authors state that overfilled (φ > 1) patterns are a superposition of in-pore hydrate, out-pore hydrate, and amorphous silica, and that WAXS cannot separate the two hydrate populations. Matching sII line positions therefore proves that some sII is present in the pan—most cleanly the external hydrate already seen by DSC near +4 °C—not that the −13.2 °C endotherm is in-pore sII. Underfilled samples that show only the −14.7 °C peak are not used as a WAXS negative control. Either provide a pore-only or underfilled diffraction control, or soften the structural claim to “sII is present and no alternative crystalline phase is detected,” and rest the in-pore assignment primarily on the DSC composition/filling dependence and pure-water ice calibrat
- §3.1–3.2 and Table 1: Composition control is load-bearing for the stated hydrate-forming window (≥ 1:16; 1:17 ice-only). Bulk 1:11 shows ~30% THF shortfall relative to the hydrate enthalpy, attributed to pre-sealing evaporative loss; the same loss is then invoked to explain why nominal 1:17 never forms confined hydrate. Pre/post pan masses (Table S1) rule out in-run leakage but do not quantify the actual THF:H2O ratio after loading. Without post-loading composition assay (or a corrected effective stoichiometry series), the sharp ≥ 1:16 threshold and the 1:17 ice-reference interpretation remain partly circular. Report effective compositions or demonstrate that the window is robust to measured loading losses.
- §3.2–3.3 and Fig. 5: Assignment of the −13.2 °C endotherm as confined hydrate (versus a second ice-like or composite confined solid) rests on a ~1.5 °C separation from the pure-water-calibrated ice peak, composition/filling dependence, and peak deconvolution with R² ≥ 0.99. That inference is reasonable but should be stress-tested: show that the two-peak fit is unique (e.g., residual comparison to a single-peak model), that the −13.2 °C area tracks THF inventory rather than total water, and that pure-water overfilled controls never produce a peak at −13.2 °C. If those checks are only in SI, elevate them; if not done, add them, because the hydrate GT branch, f_hyd trends, and “competition between hydrate and ice” narrative all depend on this assignment.
minor comments (7)
- Abstract and §3.5: ΔT_m notation is garbled in the abstract (“ΔT___”); fix consistently with the main text.
- Eq. (2) writes ΔT_m ∝ 1/(r − t), but Fig. 8 fits are presented as ΔT_m = A/r + B (and the caption text “ΔT_m = 1 r⁄” is incomplete). State explicitly whether t is absorbed into the intercept B for the hydrate branch and whether the same t = 0.6 nm pure-water value is assumed.
- Fig. 8 uses “nominal manufacturer radii” for the commercial SBA-15 series while one in-house batch is characterized by N2/KJS. Prefer a single, measured r (or r_eff) scale for all GT points, or show that manufacturer vs measured radii do not change the slope ratio.
- §3.1: Specific hydrate latent heat is taken as 261.91 J/g from literature ΔH ≈ 99.5 kJ/mol and M = 378.37 g/mol. Cite the exact source values used and note uncertainty propagation into Table 1 mass fractions.
- Figs. 3–4, 6: Symbol legend (× ice, ★ confined hydrate, + external hydrate) is helpful; ensure it appears in every multi-panel figure caption, not only some.
- Pore volume is given as ρ ≈ 0.81–0.94 cm³/g and filling as φ in cm³/g; use a single symbol for pore volume throughout and clarify how φ is normalized when batches differ.
- Typographical: “a n ELGA”, “gue st–host”, “deco mposition”, “affectin g”, “confi ned-ice” spacing errors; clean for production.
Circularity Check
No circular derivation: experimental DSC/WAXS observations and standard Gibbs–Thomson fits; literature ΔH and sII structure are external inputs, not self-defined predictions.
full rationale
This is an experimental confinement study. Load-bearing claims (filling/composition window for confined hydrate, invariant reheating temperatures under rate and cycling, linear ΔT_m vs 1/r for ice and hydrate, WAXS indexing to sII) are measured outcomes or standard analyses, not quantities forced by construction from their own definitions. Bulk mass balance uses literature ΔH_fus (ice 333.55 J/g; hydrate ≈99.5 kJ/mol) as external constants to convert peak areas—not a redefinition of the claim. Gibbs–Thomson slopes A_hydrate and A_ice are fitted to measured ΔT_m(r) under the classical external form ΔT_m = A/r + B; the paper then compares the slope ratio to thermodynamic levers (v_m/ΔH_m) from literature densities and enthalpies and infers a relative γ_sl. That is ordinary parameter extraction, not a fitted input renamed as an independent prediction of the same data. Self-citations (prior SBA-15/confined-liquid methods from the same groups) support experimental protocol and pure-water GT benchmarks; they do not supply a uniqueness theorem or ansatz that forces the THF-hydrate conclusions. WAXS peak positions are compared to known bulk sII reflections. No step reduces Eq. X to Eq. Y by definition, and no central “prediction” is statistically forced by a fit to the target. Correctness concerns about peak assignment or WAXS superposition of in-pore/out-pore hydrate are separate from circularity.
Assumptions & free parameters
free parameters (4)
- A_hydrate (GT slope for confined hydrate) =
81.5 ± 2.9 K·nm
- A_ice (GT slope for confined ice) =
73.1 ± 4.2 K·nm
- ΔH_fus,hyd (specific latent heat of sII THF hydrate) =
261.91 J/g (from ~99.5 kJ/mol)
- Nonfreezing interfacial layer thickness t =
0.6 nm
assumptions (5)
- domain assumption Classical Gibbs–Thomson relation applies to both confined ice melting and confined sII hydrate dissociation in cylindrical SBA-15 pores (eq. 2).
- domain assumption sII stoichiometry is THF·17H2O with literature fusion enthalpy ≈ 99.5 kJ/mol and bulk dissociation near 277 K.
- domain assumption The pure-water confined-ice endotherm at −14.7 °C under identical matrix/filling calibrates any endotherm at that temperature in THF–water runs as confined ice, not hydrate.
- ad hoc to paper Small reproducible THF evaporative loss during loading shifts nominal 1:17 water-rich enough to suppress confined hydrate.
- domain assumption Heterogeneous nucleation under confinement is stochastic; cooling rate mainly changes residence time for phase selection (TTT-style), not equilibrium lines.
Cite this review
Pith. "Pith review of Formation and Thermodynamic Behavior of THF-Water Hydrates in Confined Mesoporous Media." pith.science (2026). https://pith.science/paper/TXXXTQLP
@misc{pith2026260704813,
author = {Pith},
title = {Pith review of: Formation and Thermodynamic Behavior of THF-Water Hydrates in Confined Mesoporous Media},
year = {2026},
howpublished = {\url{https://pith.science/paper/TXXXTQLP}},
note = {Machine review of arXiv:2607.04813}
}
abstract
Tetrahydrofuran (THF) is a benchmark guest for probing clathrate hydrate thermodynamics because a stoichiometric aqueous solution (THF.17H2O) forms structure-II (sII) hydrate at ambient pressure with a well-defined dissociation temperature. Here, we combine differential scanning calorimetry (DSC) and wide-angle X-ray scattering (WAXS) in bulk and confined media to resolve how composition, pore filling, and cooling rate govern hydrate formation in SBA-15 mesoporous silica. Bulk DSC establishes mass-balanced enthalpies for ice and sII hydrate and confirms reversible dissociation/melting temperatures. In confinement, the heating traces separate into a Gibbs-Thomson depressed ice melt (= -14.7 $\pm$ 0.2 {\textdegree}C), an in-pore hydrate dissociation (= -13.2 $\pm$ 0.2 {\textdegree}C). Confined hydrate appears only when two criteria are met: near-percolating filling ($\phi$ = 1.0 -1.1 cm3/g) and sufficient THF ($\ge$ 1:16 mol:mol). Cooling-rate experiments (1.0 vs 0.5 {\textdegree}C/min) demonstrate that slower precooling increases the confined-hydrate fraction and reduces confined ice without shifting equilibrium temperatures: at $\phi$ = 1.1, the hydrate enthalpy rises by ~60% at 1:11 and ~54% at 1:14, but by $\le$ 17% at 1:16. Temperature-cycling tests show invariant reheating peak positions, indicating that capillarity and composition, rather than kinetic history, fix the liquidus and dissociation temperatures. WAXS indicates that the phase formed in pores is crystallographically identical to bulk sII. Finally, the variation of melting points ($\Delta$T___) plotted against inverse pore radius follows the Gibbs-Thomson law for both ice melting and hydrate dissociation, quantitatively linking the observed shifts to crystalline size and clarifying how confinement, cooling rate, and composition govern the competition between hydrate formation and water crystallization.
Reference graph
Works this paper leans on
-
[1]
(36) Xue, H.; Li, L.; Wang, Y.; Lu, Y.; Cui, K.; He, Z.; Bai, G.; Liu, J.; Zhou, X.; Wang, J
DOI: https://doi.org/10.1039/D1NR00751C. (36) Xue, H.; Li, L.; Wang, Y.; Lu, Y.; Cui, K.; He, Z.; Bai, G.; Liu, J.; Zhou, X.; Wang, J. Probing the critical nucleus size in tetrahydrofuran clathrate hydrate formation using surface -anchored nanoparticles. Nature Communications 2024, 15 (1), 157. DOI: https://doi.org/10.1038/s41467- 023-44378-6. (37) Zakrze...
-
[2]
31 (47) Brodie-Linder, N.; Dosseh, G.; Alba-Simonesco, C.; Audonnet, F.; Impéror-Clerc, M
DOI: https://doi.org/10.1063/1.1762872. 31 (47) Brodie-Linder, N.; Dosseh, G.; Alba-Simonesco, C.; Audonnet, F.; Impéror-Clerc, M. SBA- 15 synthesis: Are there lasting effects of temperature change within the first 10 min of TEOS polymerization? Materials chemistry and physics 2008, 108 (1), 73 –81. DOI: https://doi.org/10.1016/j.matchemphys.2007.09.007. ...
Reviewed July 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.