PYTHIA, a transferable machine learning framework for simultaneous defect and adsorption site optimization

cs.LG Ioannis Kouroudis, Zhaosu Gu, Poonam, Leon Muller, Alessio Gagliardi · Oct 22, 2026
Local to this browser
What it does
Machine learning for atomistic materials discovery is often bottlenecked by the cost of generating density functional theory (DFT) training data for new chemical systems. PYTHIA attacks this by freezing a pretrained MACE foundation model,...
Why it matters
PYTHIA attacks this by freezing a pretrained MACE foundation model, adapting its atomic features through a lightweight Transformer, and feeding the result into a Gaussian Process that both predicts adsorption energies and quantifies...
Main concern
The paper presents a compelling and practically motivated pipeline. Its core contribution is combining transfer learning, expressive neural representations, and uncertainty-aware Bayesian Optimization to tackle joint defect and adsorption...
Community signal
0
0 up · 0 down
Sign in to vote with arrows
AI Review AI reviewed
Plain-language introduction

Machine learning for atomistic materials discovery is often bottlenecked by the cost of generating density functional theory (DFT) training data for new chemical systems. PYTHIA attacks this by freezing a pretrained MACE foundation model, adapting its atomic features through a lightweight Transformer, and feeding the result into a Gaussian Process that both predicts adsorption energies and quantifies uncertainty for Bayesian Optimization. The framework is designed to make complex architectures usable on small DFT datasets, enabling simultaneous optimization of adsorbate placement and substrate defects in understudied materials.

Critical review
Verdict
Bottom line

The paper presents a compelling and practically motivated pipeline. Its core contribution is combining transfer learning, expressive neural representations, and uncertainty-aware Bayesian Optimization to tackle joint defect and adsorption site optimization with limited DFT data. On two chemically distinct benchmarks the authors demonstrate strong empirical gains over their prior AUGUR model and a fully trained MACE baseline. However, the manuscript provides few quantitative implementation details and withholds data and code until acceptance, which limits independent validation.

“By combining a pre-trained, frozen foundational model (MACE) with a lightweight transformer-based translation layer and a Gaussian Process head, PYTHIA is able to achieve high accuracy with limited data, while simultaneously providing a robust uncertainty estimate.”
paper · Section 4
What holds up

The architectural choices are well reasoned. Freezing a large pretrained equivariant model and inserting a Transformer adapter avoids overfitting on small adsorption datasets while preserving learned chemical priors. The authors convincingly argue that "the typical size of datasets for adsorption application ... does not suffice for the training of such a large model," and they demonstrate successful optimization on a molecular cluster and a periodic sheet with diverse dopants. The GP head naturally supplies the uncertainty estimates required for Bayesian Optimization, and the point-generation scheme incorporates physical constraints that keep candidates feasible.

“the typical size of datasets for adsorption application, however, does not suffice for the training of such a large model.”
paper · Section 2.6
“The simplest way is to generate a number of points (in our case, 10.000) that adhere to a set of constraints.”
paper · Section 2.3
Main concerns

The baseline comparisons are arguably unfair or incomplete. The "fully trained MACE" baseline is expected to struggle on a small dataset due to overparameterization, yet the authors do not compare against other parameter-efficient transfer strategies such as linear probing or LoRA, nor against contemporary foundation-model fine-tuning protocols. Their prior AUGUR model is a simpler GNN, making it an easy target rather than a strong independent baseline. Additionally, critical experimental details are missing from the main text: exact training set sizes, the specific pretrained MACE checkpoint, Transformer hyperparameters, GP kernel specification, and the holdout protocol for unseen defect configurations are all absent. The strong claim that "AUGUR fails entirely as the high complexity of the space cannot be easily captured with a simple graph neural network" is stated without quantitative metrics, and the reported Chini prediction errors—$R^2=0.83$ and MAE of $0.33$ eV—while superior to baselines, still represent chemically significant uncertainty.

“AUGUR fails entirely as the high complexity of the space cannot be easily captured with a simple graph neural network.”
paper · Section 3.2
“limited dataset does not allow for a full MACE training.”
paper · Section 3.2
Evidence and comparison

The evidence supports the authors' central claim that their specific pipeline outperforms the two internal baselines on the selected tasks. For the Chini system, PYTHIA achieved an $R^2$ of $0.83$ and an MAE of $0.33$ eV, which they report as "a significant improvement over the next best model, ranging between $80\%-100\%$, depending on the chosen metric." The silicene case showed a discovered adsorption energy improvement of over $230\%$ relative to the initial dataset minimum. Nevertheless, the comparison space is narrow; the authors do not demonstrate advantage over other transfer-learning adapters or fully trained models given comparable regularization. The impressive percentage gains reflect the success of Bayesian Optimization in discovering better minima, not solely predictive accuracy, which conflates the contributions of the surrogate model and the acquisition strategy.

“returning an $R^2$ value of $0.83$, a mean absolute error (MAE) of $0.33$ eV, and a root mean square error (RMSE) of $0.37$ eV. This represents a significant improvement over the next best model, ranging between $80\%-100\%$, depending on the chosen metric.”
paper · Section 3.1
“improvement of more than $230\%$.”
paper · Section 3.2
Reproducibility

Reproducibility is currently poor. The authors explicitly state that "Data will be made available via github upon acceptance for publication" and "Code will be made available via github upon acceptance for publication." Beyond the withheld artifacts, the main text omits essential hyperparameters and procedural details required to rerun the study: the identity of the pretrained MACE checkpoint, the architecture and training configuration of the Transformer layer, the Gaussian Process kernel and any approximations used, the exact number of DFT training points and Bayesian Optimization iterations, and the wall-clock costs. Without this information, independent reproduction and fair benchmarking are blocked.

“Data will be made available via github upon acceptance for publication.”
paper · Data Availability
“Code will be made available via github upon acceptance for publication.”
paper · Code Availability
Abstract

Machine learning has permeated material optimization, providing high value. Nevertheless, its re- liance on a high number of samples limits its applicability to relatively explored material families. In this paper we extend the applicability of machine learning to novel and understudied materials through the PYTHIA (Probabilistic Yet Transferable Hierarchical Interaction Architecture) framework. Pythia utilizes two of the most powerful machine learning concepts, foundational models and transfer learning. We employ a complex foundational model trained on a large dataset of chemically similar molecules, on a target that is physically close to our prediction target. Then, we freeze the weights and perform subspace alignment with the use of a transformer layer and add uncertainty quantification with gaussian processes. These two layers can be trained on a limited dataset yet allow us to leverage the gained knowl- edge and the studied materials and complex architectures. The resulting model is used as a surrogate for a Bayesian Optimization scheme. To showcase the robustness of the framework, we applied it to the investigation of adsorption and defect optimization of two distinct systems. Not only were we able to correctly and efficiently identify the optimal adsorption sites of even complex a

Challenge the Review

Pick a starting point or write your own. Challenges run in the background, so you can keep reading while the AI investigates.

No challenges yet. Disagree with the review? Ask the AI to revisit a specific claim.