Pith. sign in

REVIEW 2 major objections 2 minor 38 references

Short time-to-solution Quantum Monte Carlo for catalysed hydrogen synthesis. Tools give CO hydrolysis activation barriers to 1kJ/mol on Pt(111)

T0 review · 2 major / 2 minor · reviewed 2026-05-22 · grok-4.3

Challenge this review Download PDF Read on arXiv ↗

Pith's one-line read Single-determinant Quantum Monte Carlo with novel averaging computes CO hydrolysis activation barriers on Pt(111) to 0.86 kJ/mol.

desk verdict The paper claims a 0.86 kJ/mol QMC barrier for O-H dissociation on Pt(111) that matches CI, but the four-layer slab and lack of shown convergence tests leave the accuracy claim hard to judge. read the letter →

arxiv 2504.12475 v4 submitted 2025-04-16 cond-mat.mtrl-sci

Ali Bagci , Philip E Hoggan This is my paper
classification cond-mat.mtrl-sci
keywords QuantumMonteCarloPt(111)activationbarriersCOhydrolysishydrogensynthesiswater-gasshiftO-Hdissociationcatalysis
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

The paper tries to establish that Quantum Monte Carlo initialized by the ground-state Slater determinant from a simple four-layer model of Pt(111) can deliver activation barriers for the rate-limiting O-H dissociation in CO hydrolysis that match high-level configuration interaction results to within 0.16 kJ/mol. A sympathetic reader would care because this step occurs during hydrogen production via the water-gas shift reaction on platinum catalysts, and reliable barriers help predict how efficiently a surface converts reactants into hydrogen. The work focuses on an embedded active site and uses stochastic solution of the Schrödinger equation for bond breaking, which standard Hartree-Fock and DFT methods handle poorly.

What carries the argument

Embedded active site in a Pt(111) model initialized by its ground-state Slater determinant to guide single-determinant QMC for O-H dissociation energy.

What would settle it

An independent high-accuracy calculation or measurement of the O-H dissociation barrier for co-adsorbed CO and water on Pt(111) that differs from the reported QMC value by more than 1 kJ/mol.

Watch

Extended reading notes

Core claim

The author states that single-determinant work with a novel averaging procedure, using the ground-state Slater determinant of a four primitive-cell layer model of Pt(111) to initialize QMC, produces activation barriers for initial O-H bond stretch when CO is co-adsorbed with water that reach 0.86 kJ/mol accuracy, compared to 0.7 kJ/mol from the CI benchmark.

Load-bearing premise

The ground-state Slater determinant from a simple four primitive-cell layer model of Pt(111), together with the novel averaging procedure, is sufficient to produce chemically accurate barriers when used to initialize QMC for the O-H dissociation step.

Editorial extensions

If this is right

  • The barriers enable more reliable rate predictions for the water-gas shift reaction on platinum surfaces.
  • Single-determinant QMC becomes usable for heterogeneous catalysis without requiring full multi-determinant wave functions.
  • The embedded slab model plus averaging shortens the time needed to reach chemically useful accuracy.
  • The same initialization can be applied to other bond-dissociation steps that limit rates in surface reactions.

Reading between the lines

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

  • The approach could be tested on other metal surfaces such as nickel to see if the four-layer model remains adequate.
  • It may reduce the cost of screening many candidate catalysts for hydrogen synthesis by making QMC faster to run.
  • Direct comparison of the computed barriers to measured turnover frequencies on Pt(111) would check whether the numbers explain observed performance.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 2 minor

Summary. The manuscript claims that single-determinant Quantum Monte Carlo, initialized by the ground-state Slater determinant of a simple four primitive-cell layer Pt(111) slab and combined with a novel averaging procedure on an embedded active site, yields the O-H dissociation activation barrier for CO hydrolysis on Pt(111) to 0.86 kJ/mol, in close agreement with a 0.7 kJ/mol configuration-interaction benchmark.

Significance. If the numerical result and its error control hold, the work would show that QMC can deliver chemically accurate barriers for a key bond-breaking step in heterogeneous catalysis using a minimal trial wavefunction, which would be a notable technical step for surface reaction modeling where DFT and Hartree-Fock are known to struggle.

major comments (2)
  1. [Abstract / model description] Abstract and model description: the headline accuracy of 0.86 kJ/mol (versus the 0.7 kJ/mol CI value) is obtained from a ground-state Slater determinant generated on a four primitive-cell layer Pt(111) slab; no slab-thickness convergence data are shown for the work function, surface dipole, or orbital energies, yet metallic screening lengths typically demand 6–8 layers before these quantities stabilize to the 0.1 eV level that would be required to keep residual nodal errors below the claimed 0.1 kJ/mol tolerance.
  2. [Abstract] Abstract: the central numerical claim is presented without reported statistical error bars on the QMC barrier, without raw walker data, and without explicit convergence tests for the active-site embedding or the novel averaging procedure, so the stated agreement with the external CI benchmark cannot be independently verified from the information given.
minor comments (2)
  1. [Abstract] Abstract contains several typographical issues (e.g., “optimidsed”, inconsistent spacing around “Pt (111)”, and the incomplete sentence “The products are hydrogen (CO$_2$ by-product is mineralised.”).
  2. [Methods] The manuscript would benefit from an explicit statement of the number of primitive cells in the lateral directions and the k-point sampling used to obtain the trial determinant.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful review and constructive feedback on our manuscript concerning the application of single-determinant Quantum Monte Carlo to the CO hydrolysis activation barrier on Pt(111). We address each major comment below and indicate where revisions will be made to improve clarity and verifiability.

read point-by-point responses

  1. Referee: [Abstract / model description] Abstract and model description: the headline accuracy of 0.86 kJ/mol (versus the 0.7 kJ/mol CI value) is obtained from a ground-state Slater determinant generated on a four primitive-cell layer Pt(111) slab; no slab-thickness convergence data are shown for the work function, surface dipole, or orbital energies, yet metallic screening lengths typically demand 6–8 layers before these quantities stabilize to the 0.1 eV level that would be required to keep residual nodal errors below the claimed 0.1 kJ/mol tolerance.

    Authors: We agree that slab-thickness convergence is important for metallic surfaces to ensure adequate screening and control of the nodal surface in QMC. The four-layer slab was selected as the minimal model consistent with the CI benchmark setup and sufficient for the embedded active-site treatment. In the revised manuscript we will add explicit convergence tests for the work function, surface dipole, and orbital energies with increasing slab thickness, together with a discussion of how the embedding and averaging procedure limits the propagation of finite-slab effects into the reported barrier. revision: yes

  2. Referee: [Abstract] Abstract: the central numerical claim is presented without reported statistical error bars on the QMC barrier, without raw walker data, and without explicit convergence tests for the active-site embedding or the novel averaging procedure, so the stated agreement with the external CI benchmark cannot be independently verified from the information given.

    Authors: We acknowledge that the current presentation omits statistical uncertainties and detailed convergence information, which limits independent verification. In the revised manuscript we will report statistical error bars on all QMC energies and the activation barrier, include explicit tests of the active-site embedding radius and the novel averaging procedure, and provide representative raw walker statistics or a link to supplementary data to substantiate the 0.86 kJ/mol result and its agreement with the 0.7 kJ/mol CI benchmark. revision: yes

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: direct QMC computation validated against external benchmark

full rationale

The paper computes the O-H dissociation barrier on Pt(111) via QMC initialized from a four-layer slab Slater determinant plus an internal averaging procedure, then directly compares the numerical result (0.86 kJ/mol) to an independent high-level CI benchmark (0.7 kJ/mol). No load-bearing step reduces by construction to a fitted parameter, self-defined quantity, or self-citation chain; the stochastic projection solves the Schrödinger equation on the stated model, and the benchmark functions as external validation rather than input. Model-size and averaging choices are explicit assumptions whose adequacy can be tested against the benchmark or thicker-slab references, but they do not force the reported barrier value by definition.

Assumptions & free parameters 0 free parameters · 2 assumptions · 0 invented entities

The central claim rests on standard quantum-chemistry assumptions plus two modeling choices whose justification is not supplied in the abstract.

assumptions (2)
  • standard math The many-electron Schrödinger equation for the embedded active-site model is solved accurately by the chosen QMC algorithm once initialized from the Slater determinant.
    Invoked when the paper states that QMC 'solves the Schrödinger equation' for the surface reaction.
  • domain assumption A four primitive-cell layer slab with (111) orientation plus the novel averaging procedure captures the essential physics of O-H dissociation when CO and water are co-adsorbed.
    This modeling choice is required for the initialization step described in the abstract.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Short time-to-solution Quantum Monte Carlo for catalysed hydrogen synthesis. Tools give CO hydrolysis activation barriers to 1kJ/mol on Pt(111)." pith.science (2026). https://pith.science/paper/2504.12475

@misc{pith2026250412475,
  author       = {Pith},
  title        = {Pith review of: Short time-to-solution Quantum Monte Carlo for catalysed hydrogen synthesis. Tools give CO hydrolysis activation barriers to 1kJ/mol on Pt(111)},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2504.12475}},
  note         = {Machine review of arXiv:2504.12475}
}
abstract

Hydrogen synthesis is a clean, sustainable alternative to fossil fuel \cite{gals}. It has come of age: prototyping various aspects of hydrogen power are hot topics. In 9 out of 10 reactions, a solid catalyst is used. Here hydrogen production (via water-gas shift) is studied. Adsorbed reactants are optimidsed on model Pt(111). Focus is on partial O-H bond dissociation, when CO is co-adsorbed with water on this plane. hydrogen is the product. Many chemical reactions involve bond-dissociation. This process is often the key to rate-limiting reaction steps at solid surfaces. Bond-breaking is poorly described by Hartree-Fock and DFT methods, our embedded active site approach is used. We showcase Quantum Monte Carlo (QMC) methodology using the ground-state Slater Determinant of a simple four primitive-cell layer model, oriented to expose Pt (111), to initialise the QMC. This stochastic approach solves the Schr{\"o}dinger equation. It recently came of age for heterogeneous systems involving solids. During hydrolysis of carbon monoxide, initial O-H bond stretch is rate-limiting. Its dissociation energy is offset by surface Pt-H bond formation. The reactive formate (H-O-C=O) species formed by initial hydrolysis of CO, also interacting with a vicinal Pt. The products are hydrogen (CO$_2$ by-product is mineralised. A H-atom dissociates from the formate, another is desorbed from Pt(111). This yields pure hydrogen. Single-determinant work with a novel averaging procedure is compared to a high-level configuration interaction (CI) wave-function. Activation barriers are given to 0.86kJ/mol (c.f. 0.7 of the CI benchmark). Active sites embedded in metal lattice (111) faces. These trial wave-functions guide QMC.

Figures

Figures reproduced from arXiv: 2504.12475 by the authors.

Figure 1
Figure 1. Primitive unit cell for Pt, exposing (111) surfaces top and bottom. The lattice [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. Primitive unit cell for Ni, exposing (111) surfaces top and bottom. The lattice [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗

Discussion (0). Sign in to comment.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

38 extracted references · 38 canonical work pages

  1. [1]

    Green HydrogenScience et Avenir2023916, (June) p42

    work page

  2. [2]

    BBC 28.11.2022 Howarth, R.T, Jacobson, M.Z; How green is blue hy- drogen?Modelling and AnalysisDOI 10.1002/ese3.956

    work page doi:10.1002/ese3.956 2022

  3. [3]

    Quantum Monte Carlo Investigation of Two Cat- alytic Reaction Paths for Hydrogen Synthesis on Pt(111) in Recent Progress in Quantum Monte Carlo, vol 1234,2016p77-88

    Absi, N; Hoggan, P.E. Quantum Monte Carlo Investigation of Two Cat- alytic Reaction Paths for Hydrogen Synthesis on Pt(111) in Recent Progress in Quantum Monte Carlo, vol 1234,2016p77-88

    work page

  4. [4]

    Hoggan, P.E; Duchamp, N; Feschet, c; Tarragó, M. M. Quantum Monte Carlo method modeling hydrogen production using methane at Ni(111): The adsorbed formyl to CO step, Advances in Quantum Chemistry, 89, 2024, Pages 359-371

    work page 2024

  5. [5]

    Hoggan, P. E. Quantum Monte Carlo method describing supported metal catalysis: Ni(111)/alumina decomposing methane as a route to hydrogen.Adv. Quantum Chem.202188, 291-304

    work page

  6. [6]

    W; Filot, I

    Zhu, T; van Grootel, P. W; Filot, I. A. W; Sun, S-G; van Santen, R. A; Hensen, E. J. M. Microkinetics of steam methane reforming on platinum and rhodium metal surfaces.Journal of Catalysis2023297, 227-235

    work page

  7. [7]

    Gonze, X; Beuken, J.-M; Caracas, R; Detraux, F.; Fuchs, M; Rignanese, G.-M; Sindic, L; Verstraete, M; Zerah, G.; Jollet, F; Torrent, M: Roy, A; Mikami, M; Ghosez, Ph; Raty, J-Y; Allan. D.C. First-principles com- putation of material properties : the ABINIT software project. Compu- tational Materials Science 25, 478-492 (2002)

    work page 2002

  8. [8]

    Atomic and molecular adsorption on Ni(111) Bai, Y.; Kirvassilis,D.; Xu, l.; Mavrikakis, M.Surface Science2019679, 240-253

    work page

Show all 38 references
  1. [9]

    Sharma, R.O., Rantala, T.T

    Selective hydrogen production at Pt(111) investigated by Quan- tum Monte Carlo methods for metal catalysis. Sharma, R.O., Rantala, T.T. and Hoggan, P.E. Int. J. Quantum Chem. (2020) https://onlinelibrary.wiley.com/doi/abs/10.1002/qua.26198

    2020 doi

  2. [10]

    Quantum Monte Carlo approach for determining the activation barrier of water addition to carbon-monoxide adsorbed on Pt(111) within 1 kJ/mol.J

    Sharma, R.O.; Rantala, T.T.; Hoggan, P.E. Quantum Monte Carlo approach for determining the activation barrier of water addition to carbon-monoxide adsorbed on Pt(111) within 1 kJ/mol.J. Phys. Chem C.2020124, 26232-26240. 17

  3. [11]

    PRACE Promise allocation, PI Hoggan, P.E.2018

    2018

  4. [12]

    Quantum Monte Carlo calculations on dissociative chemisorption of H2 + Al(110): minimum barrier heights and their comparison to DFT values.J

    Powell, A.D.; Kroes, G.J.; Doblhoff-Dier, K. Quantum Monte Carlo calculations on dissociative chemisorption of H2 + Al(110): minimum barrier heights and their comparison to DFT values.J. Chem. Phys. 2020153, 224701

  5. [13]

    T.; Yuk, S

    Krogel, J. T.; Yuk, S. F.; Kent, P. R.C.; Cooper, V. R. Perspectives on van der Waals density functionals: the case of TiS2.J. Phys. Chem A 2020124, 47, 9867–9876

  6. [14]

    GEOMOS: semi-empirical SCF system for dealing with solvent effects and solide surface adsorption, D Rinaldi, PE Hoggan, A Cartier - QCPE Bull, 1989

    1989

  7. [15]

    Hoggan, P. E. Quantum Monte Carlo with ground-state input to investi- gate platinum doped aluminium catalyst: H2 production from adsorbed CO.Adv. Quantum Chem.202183, 155-170

  8. [16]

    Rinaldi, D: Cartier, A; Hoggan, P. E. GEOMOS: semi-empirical SCF system for dealing with solvent effects and solide surface adsorption. QCPE Bull,19899, 128

  9. [17]

    T., Energies of Adsorbed Catalytic Intermediates on Tran- sition Metal Surfaces: Calorimetric Measurements and Benchmarks for Theory

    Campbell, C. T., Energies of Adsorbed Catalytic Intermediates on Tran- sition Metal Surfaces: Calorimetric Measurements and Benchmarks for Theory. Acc. Chem. Res.2019, 52, 984–993

    2019

  10. [18]

    Fajin, J. L.C. ; Cordeiro, M. N. D.S., Probing the efficiency of platinum nanotubes for the H2production by water gas shift reaction: A DFT study. Applied Catalysis B: Environmental 2632020118301

  11. [19]

    R. J. Needs, M. D. Towler, N. D. Drummond, P. Lopez Rios, J. R. Trail. Variational and Diffusion Quantum Monte Carlo Calculations with the CASINO Code. J. Chem. Phys. 152, 1541062020. URL: https://doi.org/10.1063/1.5144288

    doi

  12. [20]

    Quantum Monte Carlo calculations on a benchmark molecule - metal surface reaction: H2 + Cu(111)

    "Quantum Monte Carlo calculations on a benchmark molecule - metal surface reaction: H2 + Cu(111)", K. Doblhoff-Dier, J. Meyer, P. E. Hoggan and G-J Kroes. J. Chem. Theory Comput.2017, 13, 7, 3208- 3219. 18

    2017

  13. [21]

    Hoggan, P. E.; Quantum Monte Carlo simulation of carbon monoxide reactivity when adsorbed at metal and oxide catalyst surfaces: Trial wave-functions with exponential type basis and quasi-exact three-body correlation. Int. J. Quantum Chem.2013, 113, 277-285

    2013

  14. [22]

    E.; Bouferguene, A

    Hoggan, P. E.; Bouferguene, A. Relative Advantages of Quantum Monte Carlo Simulation for Changing Electron Correlation: CO Reactions on Copper and Platinum Catalysts. Adv. Quantum Chem. 2014, 68, 89-104

    2014

  15. [23]

    Kroes, G. J. Towards chemically accurate simulation of molecule–surface reactions. Phys. Chem. Chem. Phys. 2012, 14, 14966-14981

    2012

  16. [24]

    A.; Busnengo, H

    Diaz, C.; Pijper, E.; Olsen, R. A.; Busnengo, H. F.; Auerbach, D. J.; Kroes, G. J. Chemically Accurate Simulation of a Prototypical

  17. [25]

    Doblhoff-Dier, J

    Diffusion Monte Carlo for accurate dissociation energies of 3d transition metal containing molecules K. Doblhoff-Dier, J. Meyer, P. E. Hoggan, G-J Kroes, and L. K Wagner. J. Chem. Theory Comput.2016, 12, 2583-2597

    2016

  18. [26]

    J.; Kinder, J

    Changlani, H. J.; Kinder, J. M.; Umrigar, C. J.; Chan, G. K-L. Ap- proximating strongly correlated wave functions with correlator product states. Phys. Rev. B.2009,<80, 245116

    2009

  19. [27]

    Beyond the locality approximation in the standard diffusion Monte Carlo method.Phys

    Casula, M. Beyond the locality approximation in the standard diffusion Monte Carlo method.Phys. Rev. B200674, 161102 R

  20. [28]

    D.; Needs, R

    Drummond, N. D.; Needs, R. J.; Sorour, A.; Foulkes, W. M. C. Finite- size errors in continuum quantum Monte Carlo calculations. Phys. Rev. B.2008, 78, 125106

    2008

  21. [29]

    J.; Towler, M

    Needs, R. J.; Towler, M. D.; Drummond, N. D.; Lopez Rios, P, J. Continuum variational and diffusion quantum Monte Carlo calculations. Phys. Condens. Matter.2010, 22, 023201

    2010

  22. [30]

    Reinhardt, P.; Hoggan, P. E. Cusps and derivatives for H2O wave- functions using Hartree-Fock Slater code: a density study. Int. J. Quan- tum Chem.2009, 109, 3191-3198

    2009

  23. [31]

    Noncovalent Interactions by the Quantum Monte Carlo Method: Strong Influence of Isotropic Jastrow Cutoff Radii

    Fanta, R.; Dubecky´, M. Noncovalent Interactions by the Quantum Monte Carlo Method: Strong Influence of Isotropic Jastrow Cutoff Radii. J. Chem. Theory Comput. 19

  24. [32]

    Farkas and H

    O. Farkas and H. B. Schlegel, Phys. Chem. Chem. Phys., 2002, 4, 11-15

    2002

  25. [33]

    P. E. Hoggan and A. Bouferguene, Quantum Monte Carlo for Activated Reactions at Solid Surfaces: Time Well Spent on Stretched Bonds, In- ternational Journal of Quantum Chemistry,2014114 p1150-1156

  26. [34]

    Growth instabilities and adsorbed impurities: a case study.Phys Rev B200877 245430

    Ben Hadj Hamouda, A.; Absi, N.; Hoggan, P.E.; Pimpinelli, A. Growth instabilities and adsorbed impurities: a case study.Phys Rev B200877 245430

  27. [35]

    MOLPRO package of ab initio programs. H.-J. Werner, P. J. Knowles, G. Knizia, F. R. Manby, M. Schütz, P. Celani, W. Györffy, D. Kats, T. Korona, R. Lindh, A. Mitrushenkov, G. Rauhut, K. R. Shamasundar, T. B. Adler, R. D. Amos, S. J. Bennie, A. Bernhardsson, A. Berning, D. L. C...

  28. [36]

    Figgen, K

    D. Figgen, K. A. Peterson, M. Dolg and H. Stoll, Energy-consistent pseudopotentialsandcorrelationconsistentbasissetsforthe5delements Hf-Pt. J. Chem. Phys.2009130 164108

  29. [37]

    G. R. Iyer, B. M. Rubenstein, Finite-Size Error Cancellation in Diffusion Monte Carlo Calculations of Surface Chemistry.Phys Chem A,2022 126, 28, 4636–4646

    2022

  30. [38]

    P.E.Hoggan, M.BensitelandJ.C.Lavalley, Anewmethodofcalculating interactions between adsorbates and metal oxide surfaces: application to the study of CO2 insertion in hydroxyl or methoxy groups on Al2O3 and TiO2. J. Mol. Struct.,1994, 320, 49-56. 20

    1994

Pith tools

Reviewed May 22, 2026 · model on record in the stance chip above.