REVIEW 3 major objections 4 minor 5 references
Assessing and Improving Machine Learning Model Predictions of Polymer Glass Transition Temperatures
T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read Adding 871 polymers to a polymer informatics model's training data cuts glass-transition temperature prediction error from 57 K to 27 K.
desk verdict The 27 K RMSE for PG-1 is an in-sample training error, not an out-of-sample prediction, so the paper's headline improvement over PG-0's 57 K is not established. 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 load-bearing mechanism is the hierarchical polymer fingerprint: each repeat unit is encoded as a vector of 128 atomic-fragment counts, 39 quantitative structure-property relationship descriptors (such as estimated surface area and fraction of rotatable bonds), and 22 morphological features (such as aromatic-ring topology and side-chain lengths). This fingerprint is mapped to $T_g$ by Gaussian process regression with a sum kernel of radial-basis-function and white-noise terms, so the prediction for a new polymer is a weighted average of similar training polymers and comes with an uncertainty estimate. Retraining after data expansion changes the kernel's similarity landscape, giving accurate predictions in regions where uncertainty was previously large.
What would settle it
Retrain the model on a random 80% of the 1,321 polymers, predict the remaining 20%, and compute the RMSE; if the held-out RMSE is much larger than 27 K, for instance 50 K or more, then the reported improvement overstates what the model will do on unseen polymers.
Extended reading notes
Core claim
The central claim is that the earlier model, trained on 450 $T_g$ values, generalizes poorly to polymers outside its training domain: its root-mean-square error on the 871 newly collected polymers exceeds 50 K. Augmenting the training set with those 871 polymers and retraining the same Gaussian-process regression pipeline yields a model whose RMSE on the full 1,321-polymer dataset is 27 K, down from 57 K, with $R^2$ rising from 0.71 to 0.92 and visibly smaller uncertainty bands. The authors attribute the improvement to the new data filling chemical spaces the old model had not seen, especially the high-$T_g$ region of roughly 500–700 K, and present the retrained model as more versatile and transferable.
Load-bearing premise
The 27 K error is measured on the same 1,321 polymers used to train the new model, with no held-out test set, so the improvement assumes that the training fit reflects how well the model will predict truly new polymers.
Editorial extensions
If this is right
- The model is now accurate to roughly the scatter of experimental $T_g$ measurements, making it useful for screening candidate polymers before synthesis.
- The same data-infusion workflow can be applied to the platform's other property models, such as band gap, dielectric constant, and tensile strength, provided comparable experimental data are collected.
- The domain of applicability explicitly expands beyond the original 450 polymers; predictions for high-$T_g$ polymers in the 500–700 K range are now backed by training data.
- Future data additions should prioritize diverse chemistries, since poor predictions in this study coincided with high Gaussian-process uncertainty.
- A retrained model with calibrated uncertainties enables inverse design, such as searching for polymers whose predicted $T_g$ falls in a target window.
Reading between the lines
- Editorial inference: because the paper evaluates the retrained model on the same 1,321 polymers used to train it, the 27 K figure is a training-fit error, not a measured generalization error; a held-out test would likely show a larger RMSE.
- Editorial inference: the pattern of large errors and high uncertainty for chemically distant polymers indicates the model is an interpolation tool, so predictions for chemistry outside the expanded space remain untrustworthy.
- Editorial inference: the success of adding data in sparse regions suggests active learning—selecting the next polymers to measure by their current prediction uncertainty—could achieve similar gains with fewer experiments.
- Editorial inference: if other property models in the platform had similarly sparse training regions, comparable data expansion would likely yield similar accuracy improvements.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports on the collection of 871 additional experimental glass transition temperature (Tg) values for polymers, augmenting an earlier dataset of 450 polymers to form a combined set of 1,321 polymers. The authors evaluate the previous Polymer Genome (PG-0) model on the new data, finding a high out-of-domain RMSE of 57 K, then retrain a Gaussian process regression model (PG-1) on the full 1,321 data points and report that the RMSE 'of prediction' drops to 27 K. The paper argues that this improvement demonstrates the value of data expansion and yields a more versatile predictive model. The main contributions are the new data collection and the qualitative demonstration that PG-0 transfers poorly to out-of-domain polymers.
Significance. The new dataset of 871 experimentally measured Tg values is a valuable community resource, and the demonstration that PG-0's predictions on these new polymers carry an RMSE of 57 K is a useful, falsifiable negative result about the transferability of models trained on a narrow chemical space. If the reported improvement to 27 K RMSE were a genuine out-of-sample prediction error, the result would be significant for polymer informatics, showing that data expansion can substantially improve predictive accuracy. However, the central quantitative claim is undermined by the evaluation protocol: the 27 K figure is computed on the same data used for training and therefore cannot be interpreted as prediction error. The paper's strengths — careful data collection, use of an established fingerprinting scheme, and consistent application of GPR — are real, but they are attached to a headline result that does not support the stated conclusion of improved predictive capability. The manuscript is more of a dataset and benchmark report than a validated methodological advance.
major comments (3)
- [Results, Fig. 2(b)] The abstract and Results state that 'The RMSE of prediction for the extended dataset ... decreased to 27 K from 57 K.' This comparison is invalid because PG-1 was trained on the full set of 1,321 polymers and then evaluated on that same full set, whereas the 57 K value for PG-0 was an out-of-sample error computed on the 871 newly collected polymers not present in PG-0's training set. The 27 K figure is therefore a training-set RMSE, not a prediction error. For Gaussian process regression with a radial basis function plus white-noise kernel, the training-set RMSE can be driven toward the noise level even when held-out error is much larger, so this number carries no evidence of generalization. To support the claim of improved predictive capability, the authors must evaluate PG-1 on data not used for training — for example, treat the 871 new polymers as a held-out test set, or perform k-fold cross-validation over the full 1,321 — and report the resulting out-of-sample RMSE alongside the 57 K baseline.
- [Methods] No validation protocol is described. The Methods section explains the fingerprinting approach and the GPR kernel but does not state how the kernel hyperparameters were optimized, whether any portion of the 1,321 polymers was held out for testing, or whether cross-validation was used. In the absence of this information, the reported R2=0.92 and RMSE=27 K in Figure 2(b) cannot be interpreted as predictive performance. The manuscript should explicitly describe the train/test split or cross-validation procedure, the hyperparameter optimization scheme, and the metrics computed on held-out data.
- [Discussion] The claim that the retrained model 'is hence more versatile' is unsupported by the reported experiments. Versatility is a statement about out-of-sample and out-of-domain generalization, yet the quantitative evidence in the manuscript is an in-sample fit on the training data. Without held-out evaluation, the discussion overstates the conclusion that adding diverse data has led to 'a significantly better predictive capability of PG.'
minor comments (4)
- [References] References 7 and 8 do not support the statement about inverse design algorithms: reference 7 is a blog post on variational autoencoders and reference 8 concerns a vendor-managed inventory optimization problem; neither is a representative citation for polymer inverse design. Please cite relevant literature on generative models for materials discovery.
- [Figure 2] The axes of the parity plots in Figure 2 are not labeled, and the parity line is not identified. Units (K) should be given on both axes, and the dashed line should be clearly described as the ideal parity line.
- [Results] The statement that 'the uncertainties calculated by GPR shown by the yellow error bars have also decreased significantly' is made without quantitative support; the authors should report the average or distribution of predicted uncertainties for the new polymers under PG-0 versus PG-1.
- [Results] The manuscript says 'the RMSE ... is greater than 50 K' for PG-0 on the new 871 polymers, but the figure caption reports RMSE=57 K. Please make the value consistent or state whether 57 K is the exact value and 'greater than 50 K' is a rounding description.
Circularity Check
PG-1's 27 K RMSE is an in-sample training error on the same 1,321 polymers used to fit the model, not an independent prediction, so the claimed improvement over PG-0's 57 K out-of-sample error is not established.
-
fitted input called prediction
[Abstract and Results, Figure 2(b) caption]
"Using the full dataset of 1321 polymers, a new ML model for Tg was built following past work. The RMSE of prediction for the extended dataset, when compared to the earlier one, decreased to 27 K from 57 K."
The paper trains PG-1 on the full 1,321-polymer dataset and then reports the 27 K RMSE on that same dataset, calling it a prediction. For a GPR model with an RBF plus white-noise kernel, the posterior mean on training inputs is an interpolation of the training targets, so the training-set RMSE is a measure of fit rather than predictive accuracy. In contrast, the 57 K RMSE for PG-0 is evaluated on the 871 newly collected polymers that were not used to train PG-0. Therefore the headline comparison is between an in-sample training error (27 K) and an out-of-sample test error (57 K), and the improvement is partially forced by construction: the model was fit to the very points on which the error is reported. No held-out subset or cross-validation split is described in Results or Methods.
full rationale
The central circularity is the evaluation protocol for PG-1. The paper uses all 1,321 polymers both to train the new model and to compute the reported 27 K RMSE, then presents this in-sample error as a predictive improvement over PG-0's out-of-sample 57 K error on unseen polymers. For a Gaussian process regressor with an RBF plus white-noise kernel, predictions at training points interpolate the training targets, so a low training RMSE is expected by construction and does not demonstrate better generalization. The fingerprinting and GPR methodology is inherited from the authors' prior work, but that is a normal methodological starting point rather than a self-citation used to force the result; no uniqueness theorem or circular definition is invoked. If the authors had evaluated PG-1 on a held-out portion of the 1,321 polymers, the reported improvement would be a genuine out-of-sample result. As reported, the 27 K RMSE is a fitted-input quantity relabeled as a prediction, so the paper deserves a partial-circularity score rather than a clean bill.
Assumptions & free parameters
free parameters (1)
- GPR kernel hyperparameters (RBF length scales and white-noise variance) =
not reported
assumptions (4)
- domain assumption Experimental Tg values from the cited handbooks and database are accurate and self-consistent.
- domain assumption The hierarchical polymer fingerprint (128 atomic fragments, 39 QSPR descriptors, 22 morphological features) captures the physical features that determine Tg.
- domain assumption Gaussian process regression with a sum of RBF and white-noise kernels is an appropriate modeling choice for Tg.
- domain assumption The 1,321 polymers are representative of the polymer chemical space of interest.
Cite this review
Pith. "Pith review of Assessing and Improving Machine Learning Model Predictions of Polymer Glass Transition Temperatures." pith.science (2026). https://pith.science/paper/5C6CH4LS
@misc{pith2026190802398,
author = {Pith},
title = {Pith review of: Assessing and Improving Machine Learning Model Predictions of Polymer Glass Transition Temperatures},
year = {2026},
howpublished = {\url{https://pith.science/paper/5C6CH4LS}},
note = {Machine review of arXiv:1908.02398}
}
read the original abstract
The success of the Materials Genome Initiative has led to opportunities for data-driven approaches for materials discovery. The recent development of Polymer Genome (PG), which is a machine learning (ML) based data-driven informatics platform for polymer property prediction, has significantly increased the efficiency of polymer design. Nevertheless, continuous expansion of the "training data" is necessary to improve the robustness, versatility and accuracy of the ML predictions. In order to test the performance and transferability of the predictive models presently available in PG (which were previously trained on a dataset of 450 polymers), we have carefully collected additional experimental glass transition temperature (Tg) data for 871 polymers from multiple data sources. The Tg values predicted by the present PG models for the polymers in the newly collected dataset were compared directly with the experimental Tg to estimate the accuracy of the present model. Using the full dataset of 1321 polymers, a new ML model for Tg was built following past work. The RMSE of prediction for the extended dataset, when compared to the earlier one, decreased to 27 K from 57 K. To further improve the performance of the Tg prediction model, we are continuing to accumulate new data and exploring new ML approaches.
Figures
Reference graph
Works this paper leans on
-
[1]
Polystyrene X11, X12, X13,…, X1M 373
-
[2]
Poly(propylene) X21, X22, X23,…, X2M 264
-
[3]
Polyethyleneoxide206 … … Fingerprinting, learning and prediction Polymer FingerprintTg(K)
-
[4]
Properties of Polymers: Their Correlation with Chemical Structure
Polyethyleneoxide X31, X32, X33,…, X3M 206 … … FingerprintingLearningNew case Instant property predictionML modelf(Xi1, Xi2, Xi3, …, XiM) = Tg(i)PolymerTg(K) New polymer? (a) (b) (c) Results Although efficient, ML models are accurate and reliable only within the domain of the dataset the model was trained on. Predictions made for cases that fall outside t...
work page 1997
-
[7]
Intuitively Understanding Variational Autoencoders
“Intuitively Understanding Variational Autoencoders” (2018, Feb 4) https://towardsdatascience.com/intuitively-understanding-variational-autoencoders-1bfe67eb5daf 8. Sadeghi, J., Sadeghi, S., and Niaki, Seyed Taghi A. "Optimizing a hybrid vendor-managed inventory and transportation problem with fuzzy demand: An improved particle swarm optimization algorith...
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.