Vitrimers represent an emerging class of sustainable polymers with self-healing capabilities enabled by dynamic covalent adaptive networks. However, their limited molecular diversity constrains their property space and potential applications. Recent developments in machine learning (ML) techniques accelerate polymer design by predicting properties and virtually screening candidates, yet the scarcity of available experimental vitrimer data poses challenges in training ML models. To address this, we leverage molecular dynamics (MD) data generated by our previous work to train and benchmark seven ML models covering six feature representations for glass transition temperature (Tg) prediction. By averaging predicted Tg from different models, the model ensemble approach outperforms individual models, allowing for accurate and efficient property prediction on unlabeled datasets. Two novel vitrimers are identified and synthesized, exhibiting experimentally validated higher Tg than existing bifunctional transesterification vitrimers, along with demonstrated healability. This work explores the possibility of using MD data to train ML models in the absence of sufficient experimental data, enabling the discovery of novel, synthesizable polymer chemistries with a wide range of desirable properties. The integrated MD–ML approach offers polymer chemists an efficient tool for accurate property prediction and designing polymers tailored to diverse applications.

Vitrimer is a new, exciting class of sustainable polymers with healing abilities due to their dynamic covalent adaptive networks. However, a limited choice of constituent molecules restricts their property space and potential applications. To overcome this challenge, an innovative approach coupling molecular dynamics (MD) simulations and a novel graph variational autoencoder (VAE) model for inverse design of vitrimer chemistries with desired glass transition temperature (T_g) is presented. The first diverse vitrimer dataset of one million chemistries is curated and T_g for 8,424 of them is calculated by high-throughput MD simulations calibrated by a Gaussian process model. The proposed VAE employs dual graph encoders and a latent dimension overlapping scheme which allows for individual representation of multi-component vitrimers. High accuracy and efficiency of the framework are demonstrated by discovering novel vitrimers with desirable T_g beyond the training regime. To validate the effectiveness of the framework in experiments, vitrimer chemistries are generated with a target T_g = 323 K. By incorporating chemical intuition, a novel vitrimer with T_g of 311–317 K is synthesized, experimentally demonstrating healability and flowability. The proposed framework offers an exciting tool for polymer chemists to design and synthesize novel, sustainable polymers for various applications.