Pith. sign in

REVIEW 4 major objections 5 minor 29 references

Multiple Diffusion-Freezing Mechanisms in Molecular Hydrogen Films

T0 review · 4 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Thin films of molecular hydrogen freeze through three detectable diffusion mechanisms, and the topmost surface layer stays mobile down to about 1 K.

desk verdict A solid experimental observation of multiple elastic anomalies in hydrogen films, undermined by an interpretation that the authors' own Table I contradicts. read the letter →

arxiv 1908.06788 v1 pith:GYHCRBWB submitted 2019-08-19 cond-mat.str-el cond-mat.mes-hall

Takahiko Makiuchi , Katsuyuki Yamashita , Michihiro Tagai , Yusuke Nago , Keiya Shirahama This is my paper · ORCID
classification cond-mat.str-elcond-mat.mes-hall
keywords molecularhydrogenisotopefilmssuperfluidityelasticanomalyquantumtunnelingofvacanciessurfacediffusiontorsionaloscillatorphasetransition
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Molecular hydrogen is a candidate quantum fluid, but bulk hydrogen solidifies before it can show superfluidity. This paper measures the elasticity of H2, HD, and D2 films in 3.9-nm pores with a torsional oscillator and finds that the films stiffen in several distinct steps as they cool, each step marking the freezing of a different kind of molecular motion. Two steps are assigned to mechanisms known from bulk solid hydrogen—classical vacancy diffusion and quantum vacancy tunneling—while a third, lowest-temperature step is assigned to diffusion within the uppermost surface layer. That surface diffusion stays active down to about 1 K, a tenth of the bulk triple point, which the authors read as evidence that hydrogen's surface is a quantum many-body state on the verge of a superfluid transition.

What carries the argument

The load-bearing object is the dissipation-peak temperature $T_p$, the temperature at which an anelastic relaxation process crosses the condition $\omega\tau(T_p)=1$ with $\tau=\tau_0\exp(E/k_B T)$. The paper fits the frequency shift and excess dissipation of each isotope by a sum of response functions $z_i$, each carrying a lognormal distribution of activation energies with median $\Delta_i$, and uses the number, position, and coverage dependence of the resulting peaks as a thermometer for each diffusion mechanism. The identity that carries the argument is the ratio $T_{p1}/T_{p2}\approx 2$, matched to the bulk hydrogen ratio $(E_v+E_b)/E_v\approx 2$, which converts the second peak into quantum vacancy tunneling. The third peak's disappearance near full-pore coverage is the operational definition of surface diffusion, since a filled pore has no free surface.

What would settle it

Measure the vacancy-formation energy $E_v$ directly in the 3.9-nm-pore films, for instance by heat-capacity or NMR; if the two activation energies recovered from the dissipation peaks do not match $E_v+E_b$ and $E_v$ for each isotope, the assignment of the middle peak to quantum vacancy tunneling is falsified.

Watch

Extended reading notes

Core claim

The central discovery is that a single adsorption system, hydrogen in porous glass, shows multiple resolvable freezing transitions in its elastic response, and that these transitions can be associated with three different diffusion mechanisms. At coverages below a monolayer only one dissipation peak appears, but at and above two layers the film shows two or three peaks whose temperatures $T_{p1} > T_{p2} > T_{p3}$ split into coverage-dependent branches. Using the bulk-solid relation $(E_v+E_b)/E_v \approx 2$, the authors identify the two higher-temperature branches with classical vacancy diffusion and quantum vacancy tunneling; the third branch, which vanishes as the pore fills and therefore tracks the free surface, is identified with surface diffusion. Fitting each anomaly with an anelastic relaxation model gives activation energies $\Delta_1$, $\Delta_2$, and $\Delta_3$, with $\Delta_3/k_B$ roughly 13, 15, and 30 K for H2, HD, and D2. The surface anomaly remains present to 1–2 K, and its ratio $\Delta_3/k_B T_{p3} \approx 12$ across isotopes matches the value $\approx 13$ found in helium films, so the authors conclude the surface layer is a quantum many-body ground state (Mott insulator or Mott glass) on the verge of a quantum phase transition to a (super)fluid state, although the transition itself is not reached.

Load-bearing premise

The claim that three separate freezing mechanisms exist depends on believing that the two known ways atoms move through solid hydrogen, with energies in a ratio of about two, still work in the tiny pores and show up as the two higher-temperature peaks.

Editorial extensions

If this is right

  • A torsional oscillator can resolve three distinct diffusion-freezing mechanisms in one adsorption system, across all three hydrogen isotopes.
  • The uppermost surface layer of hydrogen films remains mobile down to 1–2 K, about a tenth of the bulk triple point, and localizes below that with a finite surface gap.
  • The dimensionless surface-gap ratio is about 12 for all three isotopes, matching the value of about 13 found for helium films, which places the hydrogen surface in the same class of quantum many-body ground states.
  • Because the surface gap does not close, hydrogen films themselves do not reach the superfluid transition in these experiments.
  • If surface molecules are excited with phonon energy above the gap, they could in principle form a non-equilibrium superfluid, provided the relaxation from excited states is slow enough.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the surface-gap picture is right, tuning the torsional frequency over a decade should shift the lowest-temperature dissipation peak along the same Arrhenius curve; a frequency-independent peak would indicate a non-thermal relaxation instead.
  • A natural next experiment is to map the surface gap versus coverage between one and two monolayers and fit it as a power law; observing the gap close at a critical coverage would put hydrogen's surface in the same quantum-critical class as helium films.
  • The proposed non-equilibrium superfluidity could be tested by injecting phonons with frequency just above the surface gap into the film below 1 K and watching for a sudden stiffening or dissipation change in the torsional response.
  • Deliberately varying the ortho-para composition of the adsorbed hydrogen could show whether the surface gap shifts with molecular quantum statistics, separating confinement effects from intrinsic quantum behavior.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The paper reports torsional-oscillator measurements of the elastic response of H2, HD, and D2 films adsorbed in 3.9-nm porous glass, covering a range of film thicknesses up to full-pore filling. The authors observe multiple dissipation peaks and associated stiffening steps in the normalized frequency shift at temperatures well below the bulk triple point. They fit the data with a sum of thermally activated relaxation functions, each with a lognormal distribution of activation energies, and extract peak temperatures and activation energies. They attribute the three anomalies to classical thermal diffusion of vacancies, quantum tunneling of vacancies, and diffusion of molecules in the uppermost surface layer, and they argue that the surface layer is on the verge of a quantum phase transition to a (super)fluid state.

Significance. If established, the observation of a mobile surface layer persisting to about 1 K in hydrogen films would be an important step toward realizing quantum-fluid behavior in molecular hydrogen, and the extension of the elastic-anomaly method from helium and neon films to hydrogen isotopes is valuable. The experimental data show multiple reproducible dissipation peaks, and the fitting model in Eq. (4) reproduces the stepwise frequency response and the multiple loss peaks. The coverage dependence of the peak temperatures in Fig. 2 provides a useful empirical map. However, the central physical labeling of the mechanisms is not quantitatively secured by the evidence presented: the claimed Tp1/Tp2 ratio of about 2 is not satisfied by the authors' own Table I for HD and D2, and the fitted activation-energy ratios disagree even more strongly. The paper should be revised to either supply quantitative support for the three-mechanism assignment or substantially weaken the interpretive claims.

major comments (4)
  1. [Peak-identification paragraph after Fig. 2 and Table I] The central quantitative premise for assigning Tp2 to quantum vacancy tunneling is not met by the data in Table I. The text states 'The ratio Tp1/Tp2 ∼ 2 is common to all isotopes', but Table I gives Tp1/Tp2 = 6.0/2.8 ≈ 2.14 for H2, 7.0/4.1 ≈ 1.71 for HD, and 7.2/4.7 ≈ 1.53 for D2. The fitted activation-energy ratios are Δ1/Δ2 = 340/72 ≈ 4.7 for H2, 420/180 ≈ 2.3 for HD, and 420/260 ≈ 1.6 for D2. Only the HD activation ratio is close to 2; the D2 value is far below, and the H2 value is more than twice the expected ratio. Because no uncertainties are given for the tabulated entries, the deviations cannot be dismissed as statistical scatter. The identification of Tp2 with the quantum-tunneling branch of bulk self-diffusion is therefore unsupported by the evidence adduced.
  2. [Discussion of surface diffusion and QPT claim, final paragraphs] The interpretation of the third anomaly as a Mott-gap feature and the statement that the surface layer is 'on the verge of a QPT to (super)fluid state' are not established. The text itself notes that Tp3 does not decrease below 1 K and that no QPT is observed; the supporting comparison is the similarity of the dimensionless ratio Δ3/kBTp3 ≈ 12 to Δ/kBTp ≈ 13 for helium. This ratio involves fitted parameters from different models and different physical mechanisms, and no uncertainties or additional scaling tests (for example, a coverage-dependent Δ3 extrapolating to zero) are provided. The conclusion goes beyond the data and should be either reframed as a speculation or supported by additional measurements.
  3. [Eqs. (3)-(4) and fitting procedure] The three-mechanism decomposition is underdetermined at the level of the fits. Equation (3) contains, for each anomaly, the free parameters (δG/G0)i, Δi, σi, and τ0i, and Eq. (4) is a sum with N chosen to match the observed number of peaks. A good fit therefore demonstrates that the data can be represented by a sum of broadened Debye-like relaxations, but it does not by itself identify the microscopic mechanisms. The coverage and isotope trends give some support for a surface versus bulk origin, but for the Tp2 branch the quantitative tests (peak-temperature ratios, activation-energy ratios, and comparison with bulk values) fail, as noted above. Independent evidence, such as isotope-specific predictions for Δ2 or a full error analysis showing how the branch assignment is constrained by the data, is needed before the 'three different diffusion mechanisms' claim can be accepted.
  4. [Discussion of activation energies, paragraph after Table I] The comparison with bulk activation energies is internally inconsistent in a way that weakens the bulk-vacancy assignment. The authors attribute Δ1 values larger than the bulk values (Ev+Eb) to confinement in the nanopores, yet Δ2 for H2 (72 K) is smaller than the bulk vacancy formation energy Ev = 112 K cited from Ref. [23]. The text invokes larger zero-point fluctuation to reduce Δ2, but this is the opposite of the confinement argument used for Δ1 and is not quantified. Without a consistent model for how confinement and zero-point effects shift the two activation energies, the H2 value of Δ2 does not support the quantum-tunneling identification.
minor comments (5)
  1. [Text near Fig. 2] In the sentence describing the coverage dependence, 'the curve of T2 appears at the two-layer coverage' should read 'Tp2'.
  2. [Fig. 2 and Table I] No error bars are shown in Fig. 2, and Table I gives no uncertainty for the peak temperatures or activation energies. Adding uncertainties is essential for evaluating whether the branch ratios and isotope trends are meaningful.
  3. [Experimental details, Ref. [19]] The quoted gas purities are 99.99999% for H2 but only 97% for HD and 96% for D2. The possible effect of impurities on the low-temperature anomalies should be discussed, since the surface-diffusion branch is the most sensitive to contamination.
  4. [Discussion of Δ3/kBTp3] The statement that 'the ratio Δ3/kBTp3 is 12 for H2, HD, and D2 films' is only approximate: from Table I the values are 13/1.1 ≈ 11.8, 15/1.3 ≈ 11.5, and 30/2.6 ≈ 11.5. Please give the actual numbers and their uncertainties.
  5. [HD data description, Fig. 1(e)] The fourth dissipation peak in HD below about 1 K is excluded from the analysis. Because this feature lies close to the Tp3 branch, a sentence explaining why it cannot affect the surface-diffusion assignment would be helpful.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the three-mechanism assignment is an interpretive analogy to bulk hydrogen diffusion, not a derivation that reduces to its inputs.

full rationale

The paper's derivation chain is a standard anelastic-relaxation analysis: torsional-oscillator frequency and dissipation data define dissipation peaks, and Eq. (3) with lognormal activation-energy distributions is fitted to those same data. Fitting parameters (Δ_i, σ_i, τ_0i) to the peaks they describe is ordinary parameter extraction, not a prediction that reduces to an input. The assignment of Tp1, Tp2, and Tp3 to classical vacancy diffusion, quantum vacancy tunneling, and surface diffusion is an interpretive mapping based on coverage dependence and on comparison with bulk solid-hydrogen activation energies cited from Ebner and Sung [23]; the argument is analogical, and the Tp1/Tp2≈2 premise is an input from bulk theory rather than an output of the present fit. The model function and the helium/neon comparison rely on the authors' prior work [16,17], but those prior results are independent experimental studies and the anelastic model is also supported by standard references [20,21]; self-citation is therefore not load-bearing here. Any quantitative tension between the claimed Tp1/Tp2≈2 and Table I values (e.g., D2 at about 1.5, or fitted Δ1/Δ2 for H2 at about 4.7) is a correctness or robustness concern about the interpretation, not a circularity: the paper does not derive its conclusion from an equation that already contains the conclusion. No step in the paper defines a quantity in terms of the target result or renames a fitted parameter as a prediction.

Assumptions & free parameters 6 free parameters · 4 assumptions · 0 invented entities

The paper's conclusions rest on a small number of physical assumptions: the anelastic relaxation model from prior work, a geometric estimate of coverage, the transfer of bulk vacancy-diffusion mechanisms to nanopore films, and an analogy to helium for the QPT interpretation. The fitted activation energies and attempt times are free parameters of this model. No new particles, fields, or conserved quantities are introduced; the 'Mott glass' description and the proposed non-equilibrium superfluid state are interpretations or speculation, not new entities.

free parameters (6)
  • Δ1 (classical vacancy diffusion activation energy) for H2, HD, D2 = 340, 420, 420 K (Table I)
    Fitted parameter of the lognormal model; used to attribute the highest-temperature anomaly to classical vacancy diffusion.
  • Δ2 (quantum tunneling activation energy) for H2, HD, D2 = 72, 180, 260 K (Table I)
    Fitted parameter; central to the quantum-tunneling identification.
  • Δ3 (surface diffusion activation energy) for H2, HD, D2 = 13, 15, 30 K (Table I)
    Fitted parameter; used to argue the surface layer is close to a QPT.
  • τ01, τ02, τ03 (attempt relaxation times) = τ03 = 1e-9 s; τ01, τ02 in 1e-12 to 1e-30 s
    Fitted from the ratio of the low-T frequency step to the dissipation peak height.
  • σ_i (lognormal width) = ~0.3
    Fitted shape parameter for each anomaly.
  • (δG/G0)_i (relaxed modulus strength) = not reported
    Fitted amplitude of each elastic anomaly; needed to reconstruct the data but not tabulated.
assumptions (4)
  • domain assumption The elastic response of the film is governed by anelastic relaxation with a lognormal distribution of activation energies (Eq. 3), and a dissipation peak occurs at ωτ(Tp)=1.
    This model, adopted from prior work by the same group, converts measured frequency shifts and dissipation into activation energies. If the relaxation model is not valid for hydrogen films, the extracted Δ values and mechanism assignments are unsupported. Invoked in the paragraph introducing Eq. (3).
  • domain assumption Monolayer and full-pore coverages are computed from bulk molar volume via n1 = (v_m^2 N_A)^(-1/3) and nf = vp/(S v_m).
    Assumes films adopt bulk density and a simple geometric packing in the porous glass. Used to label coverage regimes in Fig. 2 and to relate Tp branches to layer numbers.
  • domain assumption The two bulk hydrogen vacancy diffusion mechanisms (classical thermal diffusion and quantum tunneling) persist in 3.9-nm pores, with activation-energy ratio (Ev+Eb)/Ev ~ 2.
    Used to assign Tp1 to classical vacancy diffusion and Tp2 to quantum vacancy tunneling. The measured Tp1/Tp2 ratios in Table I are 2.1 (H2), 1.7 (HD), and 1.5 (D2), so the claimed universal ~2 ratio is not clearly satisfied.
  • domain assumption The similarity of Δ3/kBTp3 ≈ 12 to helium's Δ/kBTp ≈ 13 implies a quantum many-body ground state (Mott glass) and proximity to a QPT.
    This is an analogy, not a derivation. The paper explicitly states that no QPT was observed, and the 'verge' claim is a speculative extrapolation from helium film behavior.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Multiple Diffusion-Freezing Mechanisms in Molecular Hydrogen Films." pith.science (2026). https://pith.science/paper/GYHCRBWB

@misc{pith2026190806788,
  author       = {Pith},
  title        = {Pith review of: Multiple Diffusion-Freezing Mechanisms in Molecular Hydrogen Films},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/GYHCRBWB}},
  note         = {Machine review of arXiv:1908.06788}
}
abstract

Molecular hydrogen is a fascinating candidate for quantum fluid showing bosonic and fermionic superfluidity. We have studied diffusion dynamics of thin films of H$_2$, HD and D$_2$ adsorbed on a glass substrate by measurements of elasticity. The elasticity shows multiple anomalies well below bulk triple point. They are attributed to three different diffusion mechanisms of admolecules and their "freezing" into localized state: classical thermal diffusion of vacancies, quantum tunneling of vacancies, and diffusion of molecules in the uppermost surface. The surface diffusion is active down to 1 K, below which the molecules become localized. This suggests that the surface layer of hydrogen films is on the verge of quantum phase transition to superfluid state.

Figures

Figures reproduced from arXiv: 1908.06788 by the authors.

Figure 1
Figure 1. FIG. 1. (a) Cross-section of the torsional oscillator. (b) S [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. The coverage dependence of the dissipation peak [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. A fitting result for (a) 2 [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

29 extracted references · 24 canonical work pages

  1. [23]

    Ebner and C

    C. Ebner and C. C. Sung, Phys. Rev. A 5, 2625 (1972)

    work page 1972

  2. [1]

    freezing

    184 cm 3, the peak value of the pore diameter dp = 3.9 nm, and the porosity p = 0 . 54. If we consider the nanopore a cylindrical pore with uniform diameter, the mean diameter is ¯d = 4 . 4 nm, which is about 11 times larger than the diameter of a hydrogen molecule. The TO was mounted on a plate thermally linked to the mixing 2 gas inlet torsional oscilla...

    work page

  3. [2]

    I. F. Silvera, Rev. Mod. Phys. 52, 393 (1980)

    work page 1980

  4. [3]

    Van Kranendonk, Solid Hydrogen: Theory of the Prop- erties of Solid H 2, HD and D 2 (Springer, New York, 1983)

    J. Van Kranendonk, Solid Hydrogen: Theory of the Prop- erties of Solid H 2, HD and D 2 (Springer, New York, 1983)

    work page 1983

  5. [4]

    M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, and E. A. Cornell, Science 269, 198 (1995)

    1995

  6. [5]

    Kasprzak, M

    J. Kasprzak, M. Richard, S. Kundermann, A. Baas, P. Jeambrun, J. M. J. Keeling, F. M. Marchetti, M. H. Szyma´ nska, R. Andr´ e, J. L. Staehli, V. Savona, P. B. Lit- tlewood, B. Deveaud, and L. S. Dang, Nature 443, 409 (2006)

    2006

  7. [6]

    A. S. Borovik-Romanov, Y. M. Bunkov, V. V. Dmitriev, and Y. M. Mukharskiy, JETP Lett. 40, 1033 (1984)

    work page 1984

  8. [7]

    S. O. Demokritov, V. E. Demidov, O. Dzyapko, G. A. Melkov, A. A. Serga, B. Hillebrands, and A. N. Slavin, Nature 443, 430 (2006)

    2006

Show all 29 references
  1. [8]

    V. L. Ginzburg and A. A. Sobyanin, Sov. Phys. JETP Lett. 15, 242 (1972)

    1972

  2. [9]

    H. J. Maris, G. M. Seidel, and T. E. Huber, J. Low Temp. Phys. 51, 471 (1983)

    1983

  3. [10]

    M. C. Gordillo and D. M. Ceperley, Physical Review Let- ters 79, 3010 (1997)

    1997

  4. [11]

    S. A. Khairallah, M. B. Sevryuk, D. M. Ceperley, and J. P. Toennies, Phys. Rev. Lett. 98, 183401 (2007)

    2007

  5. [12]

    Mezzacapo and M

    F. Mezzacapo and M. Boninsegni, Physical Review Let- ters 100, 145301 (2008)

    2008

  6. [13]

    Grebenev, B

    S. Grebenev, B. Sartakov, J. P. Toennies, and A. F. Vilesov, Science 289, 1532 (2000)

    2000

  7. [14]

    O. E. Vilches, J. Low Temp. Phys. 89, 267 (1992)

    1992

  8. [15]

    R. H. Torii, H. J. Maris, and G. M. Seidel, Phys. Rev. B 41, 7167 (1990)

    1990

  9. [16]

    Dusseault and M

    M. Dusseault and M. Boninsegni, Phys. Rev. B 97, 205403 (2018)

    2018

  10. [17]

    Makiuchi, M

    T. Makiuchi, M. Tagai, Y. Nago, D. Takahashi, and K. Shirahama, Phys. Rev. B 98, 235104 (2018)

    2018

  11. [18]

    Makiuchi, K

    T. Makiuchi, K. Yamashita, M. Tagai, Y. Nago, and K. Shirahama, J. Phys. Soc. Jpn. 88, 034601 (2019)

    2019

  12. [19]

    H. M. Roder, G. E. Childs, R. D. McCarty, and P. E. Angerhofer, Survey of the Properties of the Hydrogen Iso- topes Below Their Critical Temperatures , Technical Note 641 (National Bureau of Standards (U.S.), 1973)

    1973

  13. [20]

    The sample gas purity reported by the supplier was

  14. [21]

    R. H. Tait and J. D. Reppy, Phys. Rev. B 20, 997 (1979)

    1979

  15. [22]

    A. S. Nowick and B. S. Berry, Anelastic Relaxation in Crystalline Solids (Academic Press, New York and Lon- don, 1972)

    1972

  16. [24]

    L. H. Nosanow, J. Low Temp. Phys. 26, 613 (1977)

    1977

  17. [25]

    K. G. Sukhatme, J. E. Rutledge, and P. Taborek, J. Low Temp. Phys. 103, 301 (1996)

    1996

  18. [26]

    Maruyama, M

    M. Maruyama, M. Bienfait, F. C. Liu, O. E. Vilches, and F. Rieutord, Surf. Sci. 283, 333 (1993)

    1993

  19. [27]

    Giamarchi, P

    T. Giamarchi, P. L. Doussal, and E. Orignac, Phys. Rev. B 64, 245119 (2001)

    2001

  20. [28]

    Bloss and A

    E. Bloss and A. F. G. Wyatt, J. Low Temp. Phys. 119, 743 (2000)

    2000

  21. [99]

    99999%, 97%, and 96% for H 2, HD, and D 2, respec- tively

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.