Pith. sign in

REVIEW 4 major objections 6 minor 18 references

Predicting outcomes of catalytic reactions using machine learning

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

Pith's one-line read The paper claims that machine learning, guided by a reaction tree, can predict the products of catalytic reactions on oxygen-covered and bare gold with up to 93% accuracy.

desk verdict A promising descriptor-and-tree idea for catalytic reaction outcome prediction, but the 93% accuracy claim is not an out-of-sample result: θ and tree depth are tuned on the full dataset, so the main quantitative claim is unverified. read the letter →

arxiv 1908.10953 v1 pith:WJULUEH5 submitted 2019-08-28 physics.chem-ph cond-mat.mtrl-sci

classification physics.chem-phcond-mat.mtrl-sci
keywords machinelearningcatalyticreactionsreactionpredictionchemicalspacebagofbondstreegoldcatalysiskernelridgeregression
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

This paper claims that machine learning can predict the products of catalytic reactions on oxygen-covered and bare gold from a small experimental dataset of 145 reactions, reaching up to 93% accuracy when three candidate products are allowed. The method has the model predict where the product lies in a two-dimensional chemical space built from oxidation-number and bond-energy descriptors, then chooses the candidate molecule from a rule-generated reaction tree that sits closest to that predicted location, with a penalty for changing the number of heavy atoms. If the claim holds, reaction-outcome prediction in heterogeneous catalysis would not require millions of examples, and a learned model could flag outcomes that ordinary chemical intuition misses, such as deoxygenation of oxiranes on bare gold. It would also give a concrete path toward computational screening of catalysts and reactants.

What carries the argument

The load-bearing machinery is the reaction tree operating as a feasibility filter on top of a chemical-space regression. Reactants are encoded as extended chemical-space vectors, and the model predicts a product point $\mu = (\mu_{\mathrm{ox}}, \mu_{\mathrm{bond}})$; the reaction tree, grown by applying SMARTS patterns that encode known reaction mechanisms, supplies a discrete set of candidate products. The final prediction is the candidate minimizing $\Lambda_t(\theta) = |\hat{\mu}_{\mathrm{prod}} - \mu_t| + \theta \sum |\Delta H_t|$, where the second term penalizes differences in heavy-atom counts between reactants and products. This coupling converts a fuzzy continuous prediction into a concrete molecular structure, and the heavy-atom penalty is what pushes accuracy from 85% to 93% at three guesses.

What would settle it

Run the pipeline on a held-out set of reaction classes whose mechanisms were deliberately omitted from the SMARTS patterns; if the true product is absent from the reaction tree in a substantial fraction of cases, accuracy will drop to the rate at which the tree happens to contain the correct answer. A simpler quantitative check is to re-run the 70/30 split many times and report the mean and variance of test accuracy, since the paper's test-set accuracy (up to 100%) and full-set accuracy (93%) would then reveal how much of the result depends on the particular split.

Watch

Extended reading notes

Core claim

The central claim is that a kernel-ridge-regression model, trained on extended chemical-space vectors of reactants, predicts the product's chemical-space location well enough that matching that location against a reaction tree selects the experimentally observed product in up to 93% of cases (for three guesses and a tuned heavy-atom penalty, $\theta = 4$). The chemical space is defined by modifying the bag-of-bonds representation so that the Coulomb-style kernel uses oxidation numbers and bond energies instead of nuclear charges, and the reaction tree is grown by applying SMARTS patterns that encode known organic reaction mechanisms up to a depth of four (extended to eight when needed). The paper further reports that the model captures counterintuitive chemistry, correctly predicting deoxygenation of oxiranes and desulfurization of thiols on bare gold where chemical intuition expects simple desorption. The authors frame the result as demonstrating that small data sets can support accurate reaction-outcome prediction when a chemical-space representation is combined with mechanism-based constraints.

Load-bearing premise

The whole approach assumes the reaction tree, grown from a fixed set of known reaction mechanisms, actually contains the experimentally observed product somewhere within its depth limit; if the true product is not generated, the machine learning step cannot select it.

Editorial extensions

If this is right

  • If the central claim is right, product prediction for gold-catalyzed reactions can be done from a few hundred experiments rather than the millions of examples usually associated with machine learning in chemistry.
  • The same pipeline could be applied to other monometallic surfaces such as silver or platinum, provided enough experimental data exist for training.
  • The heavy-atom penalty and the option to allow several guesses give a practical trade-off between confidence and precision in reporting likely products.
  • The model's success on sulfur chemistry after training on oxygen chemistry suggests that the learned descriptors transfer across related heteroatom chemistries.
  • The forward reactant-to-product mapping could be inverted to search for reactants that produce a desired product, which is the inverse-design direction the paper points to in its conclusion.

Reading between the lines

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

  • The reported 93% is computed on the entire data set, so it is not a held-out generalization estimate; a repeated random-split evaluation would likely show a lower mean accuracy with noticeable variance, and the paper's own split-dependent test-set numbers already signal this.
  • A direct way to separate learning from tree coverage would be to measure the fraction of reactions whose true product appears in the reaction tree; the model's accuracy cannot exceed that fraction, so future applications should report both rates together.
  • The carbon-ring error pattern suggests that the current descriptors do not distinguish isomers well; adding graph-structure or ring-topology features is a testable fix that follows from the paper's own diagnosis.
  • The same chemical-space plus reaction-tree structure could be ported to electrocatalytic or photocatalytic reactions, but the mechanism patterns and descriptors would need to be re-derived for those settings.
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 / 6 minor

Summary. The paper presents a machine-learning pipeline for predicting the products of catalytic reactions on oxygen-covered and bare gold surfaces. Reactants are represented by extended chemical-space descriptors; kernel ridge regression, random forest, and lasso models map these descriptors to product chemical-space coordinates. Candidate products are generated by a rule-based reaction tree constructed from SMARTS patterns up to a fixed depth, and the product whose chemical-space vector is closest to the ML prediction (with a heavy-atom penalty weighted by a hyperparameter θ) is selected. Using a dataset of 145 experimental reactions, the authors report up to 93% prediction accuracy, and they claim the model can complement or exceed chemical intuition. The central contribution is the combination of a small experimental dataset, a tailored molecular representation, and a reaction tree to constrain ML predictions to chemically feasible outcomes.

Significance. If the reported accuracy were a valid out-of-sample estimate, the work would be a useful demonstration that machine learning can guide reaction-outcome prediction in heterogeneous catalysis from small datasets, an area where large datasets are lacking. The authors also propose a sensible strategy of coupling rule-based candidate generation with ML ranking, which could be transferable to other surfaces. However, the manuscript's main quantitative claim is not supported by the evaluation protocol as reported: the final accuracy is computed on the full dataset after tuning θ, γ, and tree depth, with no stable held-out evaluation. The paper's strengths are the compilation of a curated experimental dataset and the introduction of a domain-adapted descriptor, but the evaluation must be substantially reworked before the central claim can be assessed.

major comments (4)
  1. [Sec. 3.2, Fig. 5] The headline accuracy of 93% is not a predictive estimate. The text states that accuracy is reported on the entire data set after choosing θ=4 as the value that 'optimizes the prediction accuracy.' Since θ is tuned on the same full dataset on which accuracy is reported, and γ=3 counts a reaction as correct if any of three ranked molecules matches, the reported figure conflates in-sample fitting with generalization. Please provide a proper nested cross-validation or repeated train/test split protocol in which θ, γ, and any tree-depth choices are fixed on training folds only, and report the mean and variance of the resulting held-out accuracy.
  2. [Sec. 4] The retrospective increase of reaction-tree depth from 4 to 8 for reactions 'where the observed product did not appear in the trees of depth 4' is a label-aware modification of the candidate generator. Once the observed product is inserted into the candidate set after the fact, the ML ranking is being evaluated on a set that contains the answer by construction, so the improved accuracies of 89%, 95%, and 97% do not measure predictive performance. Please either fix the tree depth a priori for all reactions, or report the fraction of cases where the observed product is absent from depth-4 trees as a recall limitation separate from the ML ranking accuracy.
  3. [Sec. 3.2] The statement that the test set can yield up to 100% accuracy but 'strongly depends on how we make the training/test set split' is not accompanied by any quantitative details. No random seed, number of splits, standard deviation, or distribution of test-set accuracies is given, and the claim 'We believe that our trained model is not overfitting' is an assertion rather than a demonstrated result. Please report the held-out accuracy over many random splits for a fixed hyperparameter configuration, including the mean and spread.
  4. [Eq. (2)] The heavy-atom penalty Λ_t(θ) uses ΔĤ_t, which is itself the output of a machine learning model, and the weight θ is tuned to maximize accuracy on the full dataset. This couples the hyperparameter tuning to the evaluation labels and further weakens the reported accuracy. The penalty weight should be selected using only training data (e.g., by inner cross-validation), and the final evaluation must be on data not used for any stage of model or hyperparameter selection.
minor comments (6)
  1. [Abstract and Sec. 1] The abstract describes the dataset as 'less than two hundred reactions' while the text specifies 145 reactions; please use a consistent number throughout.
  2. [Sec. 2.1] The exclusion of reactions forming radicals and complete-combustion products should be discussed as a potential selection bias, since these are real outcomes in some catalytic conditions.
  3. [Sec. 2.3] The paper repeatedly refers to Supplementary Information for SMARTS patterns and descriptor lists, but the arXiv version does not include the SI. Please make the SI available or include the essential details in the manuscript.
  4. [Fig. 3] The t-SNE projections would benefit from stating the perplexity parameter and the random seed, as t-SNE is stochastic and the figures are used to claim clustering of molecular classes and reaction types.
  5. [Eq. (1)] The kernel width σ of the kernel ridge regression is said to be determined by cross-validation, but no range of values or final selected values are reported; please provide these details for reproducibility.
  6. [Author list] The first author's name appears as 'T revor David Rhone' in the header; this appears to be a typographical error.

Circularity Check

2 steps flagged · score 7.0 of 10

The headline 93% accuracy is an in-sample, label-tuned figure: θ and γ are optimized on the full dataset, and tree depth is extended using observed products.

  1. fitted input called prediction [Section 3.2, Eqs. (2)-(3) and Fig. 5; also Section 2.3]
    "If we examine the test set only, we can find up to 100% prediction accuracy for the test set. However, this strongly depends on how we make the training/test set split. ... for the remainder of the discussion we report the prediction accuracy on the entire data set. ... We choose a θ value that optimizes the prediction accuracy. Imposing the heavy atom penalty with the optimal weighting, θ = 4, increases the prediction accuracy to 93% for γ = 3."

    The headline accuracy is computed on the entire data set, not on a held-out test set, and θ is chosen by optimizing accuracy on that same entire data set. The paper's own Section 2.3 describes a 70/30 split whose test accuracy should approximate expected prediction accuracy, but Section 3.2 abandons that split because test results depend strongly on the split. With γ = 3, a reaction is counted correct if any of the three closest tree molecules matches, further relaxing the metric. The 93% figure is therefore a fitted, in-sample optimum over θ and γ, not an out-of-sample predictive accuracy; no nested cross-validation or pre-selection of θ/γ is reported.

  2. self definitional [Section 4 (Discussion)]
    "We tested a larger tree depth of 8 for reactions where the observed product did not appear in the trees of depth 4, and this resulted in correct predictions from our methodology in all of these cases (for γ = 2 and 3). With this larger depth for these specific cases, the prediction accuracy increased to 89%, 95% and 97% for γ = 1, γ = 2, and γ = 3, respectively."

    The predicted product is defined by Eq. (3) as the argmin over the reaction tree, so the tree is the candidate set from which a correct prediction must be drawn. Here the tree depth is extended from 4 to 8 specifically for reactions whose observed product was missing from the depth-4 tree; the observed labels are used to decide which candidate sets to expand. Counting those label-informed expansions as correct predictions means that the improved accuracy partly measures the retrospective re-insertion of the true product into the candidate set, not prospective ML selection. The 95-97% figures are thus conditional on a candidate generator modified using the labels of the scored reactions.

full rationale

The core methodology—mapping reactants to a chemical-space value and matching to a reaction-tree candidate—is not itself circular: the BoB-derived descriptors, SMARTS patterns, and experimental reaction list are external inputs, and the t-SNE clustering demonstrations are self-contained. I find no load-bearing self-citation. The circularity burden is concentrated in the headline quantitative claim. Section 2.3 describes a 70/30 train/test split and states that the test set serves as unseen data, but Section 3.2 switches to reporting accuracy on the entire data set because the test-set result 'strongly depends on how we make the training/test set split.' θ is then chosen to optimize accuracy on that entire set, and γ = 3 is used to allow any of three guesses to count as correct; the resulting 93% is an in-sample optimum rather than a predictive accuracy. Section 4 compounds this by extending tree depth from 4 to 8 for reactions where the observed product did not appear in the depth-4 tree, using the observed product labels to decide which trees to deepen, and then reporting improved accuracy on those label-informed cases. No nested cross-validation, fixed pre-selection of θ and γ, or pre-registered tree depth is reported, so the central claim that the model predicts catalytic reaction outcomes with up to 93% accuracy is not supported by an evaluation that avoids fitting the reported metric. This is partial circularity—fitted values and label-aware candidate expansion reported as predictions—not full tautology, because the model could in principle be validated by a disciplined held-out protocol. Score 7.

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

The central claim rests on a small experimental dataset, a hand-tuned chemical space representation, and a rule-based reaction tree that must contain the true product. The free parameters are dominated by tuning of the heavy-atom penalty and tree depth against the evaluation set; no new physical entities are introduced.

free parameters (4)
  • theta (heavy atom penalty weight) = 4
    Tuned in Fig. 5 to maximize prediction accuracy on the full dataset; this tuning on the evaluation set makes the reported accuracy a fitted quantity.
  • Maximum reaction tree depth = 4 (with exceptions at 8)
    Chosen to balance tree size and inclusion of observed products; depth was extended to 8 for reactions where the observed product did not appear at depth 4, a post hoc adjustment.
  • gamma (number of guesses) = 3
    Accuracy increases with gamma; the headline 93% corresponds to gamma=3, i.e., the product is considered correct if it is among three guesses.
  • ML hyperparameters (e.g., KRR kernel width) = not specified in text
    Tuned via k-fold cross-validation on the training set; these are standard hyperparameters but were fit to data, and their values are not given.
assumptions (4)
  • domain assumption The SMARTS patterns encode all relevant reaction mechanisms for the studied Au surface chemistry.
    The reaction tree is grown by applying SMARTS patterns modeled on known organic mechanisms; if a real reaction proceeds via an unencoded mechanism, the true product will not be a candidate. Section 2.3.
  • domain assumption Euclidean distance in the modified BoB chemical space (using oxidation number and bond energy) is a valid similarity measure for selecting the correct product.
    The product is chosen as the reaction-tree molecule with chemical space vector closest to the ML prediction; no validation that this distance correlates with actual reaction outcome likelihood is provided. Eq. (2), Section 3.2.
  • domain assumption The 145 experimental reactions from UHV studies on Au(111)/Au(110) are accurate and representative of the reaction outcomes on these surfaces.
    The dataset is assembled from literature UHV experiments; no error bars or reproducibility assessment for the experimental labels is given. Section 2.1.
  • domain assumption Heavy atom count is conserved or predictable; reactions that decompose into multiple extended carbon chains are excluded.
    The method omits reactions with decomposition to multiple products with extended carbon chains, and uses a learned Delta_H penalty to enforce heavy-atom conservation. Section 2.3.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Predicting outcomes of catalytic reactions using machine learning." pith.science (2026). https://pith.science/paper/WJULUEH5

@misc{pith2026190810953,
  author       = {Pith},
  title        = {Pith review of: Predicting outcomes of catalytic reactions using machine learning},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/WJULUEH5}},
  note         = {Machine review of arXiv:1908.10953}
}
read the original abstract

Predicting the outcome of a chemical reaction using efficient computational models can be used to develop high-throughput screening techniques. This can significantly reduce the number of experiments needed to be performed in a huge search space, which saves time, effort and expense. Recently, machine learning methods have been bolstering conventional structure-activity relationships used to advance understanding of chemical reactions. We have developed a model to predict the products of catalytic reactions on the surface of oxygen-covered and bare gold using machine learning. Using experimental data, we developed a machine learning model that maps reactants to products, using a chemical space representation. This involves predicting a chemical space value for the products, and then matching this value to a molecular structure chosen from a database. The database was developed by applying a set of possible reaction outcomes using known reaction mechanisms. Our machine learning approach complements chemical intuition in predicting the outcome of several types of chemical reactions. In some cases, machine learning makes correct predictions where chemical intuition fails. We achieve up to 93% prediction accuracy for a small data set of less than two hundred reactions.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

18 extracted references · 17 canonical work pages

  1. [1]

    Tropsha, Molecular Informatics, 2010, 29, 476–488

    1 A. Tropsha, Molecular Informatics, 2010, 29, 476–488. 2 S. Szymku ´c, E. P. Gajewska, T. Klucznik, K. Molga, P. Dittwald, M. Startek, M. Bajczyk and B. A. Grzybowski, Computer-Assisted Synthetic Planning: The End of the Begin - ning, 2016, vol. 55, pp. 5904–5937. 3 J. N. Wei, D. Duvenaud and A. Aspuru-Guzik, ACS Central Science, 2016, 2, 725–732. 4 W. L...

  2. [8]

    Gasteiger , M

    36 J. Gasteiger , M. G. Hutchings, B. Christoph, L. Gann, C. Hi ller , P. Löw , M. Marsili, H. Saller and K. Yuki, Organic synthesis, reactions, and mechanisms / , Springer-Verlag, Berlin ;, 1987, pp. 19–73. 37 D. Rappoport, C. J. Galvin, D. Y. Zubarev and A. Aspuru- Guzik, Journal of Chemical Theory and Computation , 2014, 10, 897–907. 38 B. Xu, X. Liu, ...

  3. [14]

    Sanchez-Lengeling and A

    67 B. Sanchez-Lengeling and A. Aspuru-Guzik, Science, 2018, 361, 360–365. 68 K. Hansen, F. Biegler , R. Ramakrishnan, W. Pronobis, O. A. von Lilienfeld, K.-R. Müller and A. Tkatchenko, The Journal of Physical Chemistry Letters , 2015, 6, 2326–2331. 69 A. J. Chowdhury , W. Yang, E. Walker , O. Mamun, A. Heyden and G. A. Terejanu, The Journal of Physical Ch...

  4. [20]
  5. [48]

    72 S. Kim, P. A. Thiessen, E. E. Bolton, J. Chen, G. Fu, A. Gindu - lyte, L. Han, J. He, S. He, B. A. Shoemaker , J. Wang, B. Yu, J. Zhang and S. H. Bryant, Nucleic Acids Research , 2016, 44, D1202–D1213. 73 R. Tibshirani, G. James, D. Witten and T. Hastie, An introduc- tion to statistical learning-with applications in R ,

  6. [66]

    Zugic, S

    56 B. Zugic, S. Karakalos, K. J. Stowers, M. M. Biener , J. Bien er , R. J. Madix and C. M. Friend, ACS Catalysis , 2016, 6, 1833–

  7. [171]

    64 J. C. F. Rodríguez-Reyes, C. M. Friend and R. J. Madix, Catal- ysis Letters, 2018, 148, 1–5. 65 G. Skoraczy ´nski, P. Dittwald, B. Miasojedow , S. Szymku´c, E. P. Gajewska, B. A. Grzybowski and A. Gambin, Scientific Reports , 2017, 7,

  8. [199]

    van der Maaten and G

    77 L. van der Maaten and G. Hinton, JOURNAL OF MACHINE LEARNING RESEARCH, 2008, 9, 2579–2605. 78 T. Hastie, R. Tibshirani and J. Friedman, The Elements of Sta- tistical Learning, Springer New York Inc., New York, NY, USA,

Show all 18 references
  1. [306]

    15 Z. Li, S. Wang, W. S. Chin, L. E. Achenie and H. Xin, J. Mater . Chem. A, 2017, 5, 24131–24138. 16 C. W. Coley , W. H. Green and K. F. Jensen, Acc. Chem. Res. , 2018, 51, 1281–1289. 17 L. Grajciar , C. J. Heard, A. A. Bondarenko, M. V. Polynski, J. Meeprasert, E. A. Pidko a...

  2. [376]

    40 B. Xu, X. Liu, J. Haubrich and C. M. Friend, Nature Chemistry, 2010, 2, 61–65. 41 X. Liu, B. Xu, J. Haubrich, R. J. Madix and C. M. Friend, Jour- nal of the American Chemical Society , 2009, 131, 5757–5759. 42 C. G. F. Siler , T. Cremer , J. C. F. Rodriguez-Reyes, C. M. Fri...

  3. [1839]

    57 B. K. Min, X. Deng, X. Liu, C. M. Friend and A. R. Alemozafar , ChemCatChem, 2009, 1, 116–121. 58 D. A. Outka and R. J. Madix, Journal of the American Chemical Society, 1987, 109, 1708–1714. 59 X. Liu, PhD thesis, Harvard University ,

  4. [2001]

    DOI.png" is available in

    This figure "DOI.png" is available in "png" format from: http://arxiv.org/ps/1908.10953v1 This figure "dates.png" is available in "png" format from: http://arxiv.org/ps/1908.10953v1

  5. [2009]

    60 C. R. O’Connor , F. Hiebel, W. Chen, E. Kaxiras, R. J. Madix and C. M. Friend, Chemical Science, 2018, 9, 3759–3766. 61 D. M. Jaffey and R. J. Madix, Journal of the American Chemical Society, 1994, 116, 3020–3027. 62 D. M. Jaffey and R. J. Madix, Journal of the American Che...

  6. [2012]

    26 G. V. Carrera, S. Gupta and J. Aires-de Sousa, J. Comput. Aided. Mol. Des. , 2009, 23, 419–429. 27 J. N. Wei, D. Duvenaud and A. Aspuru-Guzik, ACS Cent. Sci. , 2016, 2, 725–732. 28 Q. Y. Zhang and J. Aires-De-Sousa, J. Chem. Inf. Model. , 2005, 45, 1775–1783. 29 M. D. Hughe...

  7. [2013]

    Pedregosa, G

    74 F. Pedregosa, G. Varoquaux, A. Gramfort, V. Michel, B. Thirion, O. Grisel, M. Blondel, P. Prettenhofer , R. Weiss , V. Dubourg et al., J. Mach. Learn. Res. , 2011, 12, 2825–2830. 75 A. Seko, T. Maekawa, K. Tsuda and I. Tanaka, Phys. Rev. B , 2014, 89, 054303. 76 F. Feng, L....

  8. [2017]

    34 M. Rupp, A. Tkatchenko, K.-R. Müller and O. A. von Lilienfe ld, Phys. Rev. Lett. , 2012, 108, 058301. 35 J. Gilmer , S. S. Schoenholz, P. F. Riley , O. Vinyals and G. E . Dahl, arXiv:1704.01212,

  9. [2018]

    24 M. H. Segler , M. Preuss and M. P. Waller , Nature, 2018, 555, 604–610. 25 D. M. Lowe, PhD thesis, Pembroke College,

  10. [4443]

    Takigawa, K.-i

    13 I. Takigawa, K.-i. Shimizu, K. Tsuda and S. Takakusagi, RSC Adv., 2016, 6, 52587–52595. 14 H. Li, Z. Zhang and Z. Liu, Catalysts, 2017, 7,

Pith tools

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