# A Comprehensive Study on the Styrene-GTR Radical Graft Polymerization: Combination of an Experimental Approach, on Different Scales, with Machine Learning Modeling

Cindy Trinh, Sandrine Hoppe, Richard Lainé, Dimitrios Meimaroglou

<figure><img src="image_1.png" /></figure>

► To cite this version:

Cindy Trinh, Sandrine Hoppe, Richard Lainé, Dimitrios Meimaroglou. A Comprehensive Study on the Styrene-GTR Radical Graft Polymerization: Combination of an Experimental Approach, on Different Scales, with Machine Learning Modeling. Macromol. 2023, 3 (1), pp.79-107. 10.3390/macromol3010007. hal-04062623

HAL Id: hal-04062623

https://hal.univ-lorraine.fr/hal-04062623

Submitted on 7 Apr 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# A Comprehensive Study on the Styrene–GTR Radical Graft Polymerization: Combination of an Experimental Approach, on Different Scales, with Machine Learning Modeling

Cindy Trinh ID, Sandrine Hoppe ID, Richard Lainé and Dimitrios Meimaroglou *ID

Laboratoire Réactions et Génie des Procédés, Université de Lorraine, CNRS UMR7274, LRGP, F-54000 Nancy, France

* Correspondence: dimitrios.meimaroglou@univ-lorraine.fr

**Abstract:** The study of the styrene–Ground Tire Rubber (GTR) graft radical polymerization is particularly challenging due to the complexity of the underlying kinetic mechanisms and nature of GTR. In this work, an experimental study on two scales (~10 mL and ~100 mL) and a machine learning (ML) modeling approach are combined to establish a quantitative relationship between operating conditions and styrene conversion. The two-scale experimental approach enables to verify the impact of upscaling on thermal and mixing effects that are particularly important in this heterogeneous system, as also evidenced in previous works. The adopted experimental setups are designed in view of multiple data production, while paying specific attention in data reliability by eliminating the uncertainty related to sampling for analyses. At the same time, all the potential sources of uncertainty, such as the mass loss along the different steps of the process and the precision of the experimental equipment, are also carefully identified and monitored. The experimental results on both scales validate previously observed effects of GTR, benzoyl peroxide initiator and temperature on styrene conversion but, at the same time, reveal the need of an efficient design of the experimental procedure in terms of mixing and of monitoring uncertainties. Subsequently, the most reliable experimental data (i.e., 69 data from the 10 mL system) are used for the screening of a series of diverse supervised-learning regression ML models and the optimization of the hyperparameters of the best-performing ones. These are gradient boosting, multilayer perceptrons and random forest with, respectively, a test $R^2$ of 0.91 ± 0.04, 0.90 ± 0.04 and 0.89 ± 0.05. Finally, the effect of additional parameters, such as the scaling method, the number of folds and the random partitioning of data in the train/test splits, as well as the integration of the experimental uncertainties in the learning procedure, are exploited as means to improve the performance of the developed models.

**Keywords:** radical graft polymerization; styrene; ground tire rubber; artificial intelligence; machine learning

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Trinh, C.; Hoppe, S.; Lainé, R.; Meimaroglou, D. A Comprehensive Study on the Styrene–GTR Radical Graft

**Polymerization:** Combination of an Experimental Approach, on Different Scales, with Machine Learning Modeling. *Macromol. 2023*, *3*, 79–107.

https://doi.org/10.3390/macromol3010007

Academic Editor: Ana María Díez-Pascual

Received: 19 December 2022

Revised: 14 February 2023

Accepted: 17 February 2023

Published: 22 February 2023

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

As the consequences of climate change have been more and more visible in recent years, it has become of prior importance for all sectors to reduce the environmental impact of their activities, which is also encouraged by the increasing number of policies going in this direction. In this sense, greener ways to deal with the increasing amount of end-of-life tires are needed to avoid their ending in landfills or burning for energy recovery [1–4].

One environmentally-friendly and popular solution consists of transforming the rubber parts of tires into a micrometric powder, commercially called Ground Tire Rubber (GTR), by mechanical grinding. This powder is subsequently used as raw material or additive in other products of interest for diverse applications. In particular, the elastomeric properties of GTR make it an interesting filler agent to improve the mechanical properties of a wide range of materials including thermoplastics, thermosets, virgin or composite
rubbers, concrete, bitumen or asphalt [1–3,5–9]. For example, the introduction of GTR fillers in polystyrene (PS), which is investigated in this work, can produce a composite material with improved stress cracking resistance and impact strength, with respect to the brittle pure PS matrix, and with reduced material costs.

However, the different structures of the two phases (i.e., the PS matrix and GTR, GTR having the form of a 3D crosslinked network) present a challenge in terms of their compatibility and the resulting poor mechanical properties of the final blend. Accordingly, several compatibilization strategies can be envisioned, such as devulcanization and reclamation or surface activation [2,3,5,6,10,11]. In this work, a chemical surface activation strategy in the form of in situ cross-linking by grafting polymerization, which is among the most prominent GTR surface modification techniques [6], is employed to produce PS-GTR composites. More specifically, PS is directly grafted onto the surface of GTR particles that are present during the radical polymerization of styrene. Similar in situ bulk radical polymerization techniques have been used to finely disperse and incorporate not only rubber particles in polymer matrices [4,12–14], but also other nanofillers, such as for the synthesis of graphene-based nanocomposite materials [15–21], carbon tube-based nanocomposite materials [22,23] and many others [24–26].

A second challenging aspect of this process lies in the complexity of the GTR itself as raw material [2]. GTR is effectively obtained from tires of diverse origins and/or with unknown compositions of rubbers and additives, depending on the vehicle type and the required performance. This limits the mastery of the process and reduces the capacity of a detailed modeling of the system to better understand the link between the process parameters (e.g., temperature/time of the reaction and initial compositions of reactants) and the final productivity and product quality indicators (e.g., styrene conversion rate and grafting efficiency). For example, a previous work demonstrated that additives and moieties found in varying compositions in GTR, such as carbon black, display a significant effect on the course of the polymerization [27]. Nevertheless, developing a model capable of handling this inherent system complexity and predicting the effects of the process conditions on the desired indicators of the produced GTR-based composite would be extremely beneficial to the associated recycling and circular economy sectors.

This work builds on previously reported results of an ongoing study on the styrene-GTR radical graft polymerization system [27–29]. More specifically, a part of this ongoing study has already focused on the modeling of the system via the introduction of a generalized kinetic scheme that has been employed to describe the evolution of the polymerization in the presence of different amounts of GTR and/or benzoyl peroxide (BPO) initiator, under different reaction temperatures [29]. The mathematical model proposed in that work, presented an initial framework for the consideration of the chemical reactions that might occur in the system, as a consequence of the presence of GTR (i.e., as opposed to the presence of pure polybutadiene rubber that is considered in most reported studies [30–33]) as well as a methodology for the quantification and the gradual “masking” of reaction sites on the surface of GTR. In the present work, a complementary approach is adopted as means to circumvent the manque of knowledge on the exact composition of the system and the associated chemical developments. Accordingly, the present study focuses exclusively on the case where GTR is present in the system, contrary to [29] where both the pure styrene homopolymerization and the GTR-grafting cases were investigated.

Machine learning (ML) is a prominent tool in polymer materials design and property prediction that allows saving significant time and effort, related to the development of purely phenomenological modeling approaches, when the underlying phenomena are not completely elucidated and/or highly complex [34,35]. Accordingly, this work presents the implementation of a data-driven modeling approach, on the basis of supervised ML techniques, under varying operating conditions. The developed ML models are implemented as standalone regression models, for the prediction of the monomer conversion developments in an attempt to assess their suitability in modeling this highly complex polymerization system, without using any prior knowledge or mechanistic description of
the underlying phenomena, as it is the case in previous works [27–29]. As a perspective for future work, the findings of this work will be used in the development of a hybrid modeling approach where both phenomenological and ML models of the system will be combined.

The proposed ML modeling framework necessitates the existence of experimental data to serve during the training of the models and the subsequent testing of their generalization capacity. However, the use of the previously reported bench-scale experimental setup [28] was prohibitive for the production of sufficient amounts of data for the ML model (i.e., under a wide range of process conditions), due to its highly time-consuming nature. In addition, this system also presented a limited capacity of controlling the homogeneity of temperature and composition within the reactor, especially at high GTR loading [28]. So in order to produce the required experimental data for this work, with the lowest possible uncertainty and within a reasonable amount of time, the grafting polymerization experiments were performed on laboratory-scale experimental setups, namely a small-sized system (~10 mL) and a medium-sized one (~100 mL). In the small-sized system, a parallel realization of grafting polymerizations in sealed test tubes, placed in a heating bath and without internal agitation, was carried out. On the other hand, the medium-sized system was based on an ad hoc experimental setup allowing the simultaneous realization of the polymerization reaction in six glass reactors, under controlled temperature and internal agitation. Both systems served to produce multiple samples rapidly and to investigate the effect of scale-up on the thermal and compositional homogeneity within the reacting mixture. In both cases, small volumes of reactants were used in order to allow using the complete reactor content in the post-reaction analyses, thus limiting the uncertainties related to the sampling of the heterogeneous mixture.

Note that, most reported studies are either limited to low GTR loadings (i.e., up to 30% wt.) or concern different polymer matrices and/or compatibilization strategies [10, 12, 36–42]. At the same time, the precision and repeatability of the results are rarely discussed. To the authors’ knowledge, this is the first time that the styrene–GTR radical graft polymerization system has been comprehensively investigated, for a wide range of GTR content, by combining a two-scale experimental approach with a ML modeling strategy. All the produced experimental data, as well as the developed mathematical models, are readily available through a GitHub repository. In addition, an extensive analysis of the uncertainty introduced in the different steps of the experimental process is also included in the discussion of the results.

## 2. Experimental Section

### 2.1. Materials

The reagents used for the polymerization reactions are styrene monomer Reagent-Plus® (with a purity ≥ 99% and stabilized with 4-tert-butylcatechol) and BPO initiator Luperox® A75. They were purchased from Sigma-Aldrich Chemie GmbH (Steinheim, Germany) and Sigma-Aldrich Co. (St. Louis, MO, USA) respectively and both were used without further purification. GTR particles were obtained from DeltaGom France and used without further purification as well. The characteristics of GTR have been discussed elsewhere [27,28]. Even if the presence of stabilizers and impurities can greatly affect polymerization kinetics, these were not removed within the perspective of implementing this process at a larger scale. Therefore, the idea in this work was to use GTR as-received without further purification with organic solvents to avoid the generation of chemical wastes and to reduce production costs.

### 2.2. Small-Sized System

On the small scale, reactions took place in glass test tubes containing a total of 1.5 g of reactants (i.e., styrene, BPO and GTR) at various proportions. The test tubes were first filled with the desired amounts of GTR. The mixtures of monomer and initiator were preliminary prepared in beakers and then introduced directly onto the GTR, in the corresponding test tubes, followed by sealing of the tubes. Once the oil bath reached the target temperature,
about ten test tubes were introduced in the bath and kept submerged during the whole reaction time. Occasionally, the sealing of the tubes was loosened to allow pressure relief in case of gas buildup. At the end of each reaction, the test tubes were immediately placed into a cold water bath for cooling. Due to the very low volumes, the heating and cooling times of the tubes was limited to a few minutes (i.e., typically less than three minutes). Finally, the test tubes were unsealed and placed in vacuum oven at 40 °C until constant weight, signifying the complete evaporation of any remaining styrene. The measured mass loss during this stage served in the calculation of styrene conversion, according to:

$$X_{\text{styrene}} = \frac{m_{i,\text{styrene}} - m_{f,\text{styrene}}}{m_{i,\text{styrene}}} \quad (1)$$

where $m_{i,styrene}$ denotes the mass of styrene that was introduced in the test tube and $m_{f,styrene}$ is the recorded mass loss during vacuum evaporation.

### 2.3. Medium-Sized System

Figure 1 depicts a schematic of the experimental system, which was designed to perform simultaneously six polymerization reactions under identical heating, stirring, inerting and refrigeration conditions. More precisely, Figure 1a corresponds to a side view of the system (i.e., showing only three of the six reactors), while the top view proposed in Figure 1b enables to better visualize all six positions. Besides these schematics, photos of the actual system are provided in Appendix A. In the following, the elements of the system are described.

* The reactors were designed with a limited volume (<100 mL) to keep the duration of sample analysis after polymerization reasonable (note that the whole reactor content is analyzed) while enabling reactions with sufficient amounts of GTR. The collars are large enough to facilitate the introduction of the reactants and the recovery of the products, the latter being more or less viscous depending on the operating conditions. The reactor diameter remains constant along the reactor height (except the bottom). This facilitates the introduction and removal of the agitator. The reactor material (glass) facilitates the observation of the mixing conditions. Finally, a lateral orifice enables the punctual introduction of a thermocouple, inside the reactor, for temperature measurement.

* A 6-cm thick aluminum plate, containing 6 drilled reactor holdings, was used for reactor temperature control (Figure 1b,c). This plate was heated via conduction by an electrical heating plate on which it was directly placed. Silicon oil was also added in the holding positions to maximize the heat-transfer between the aluminum plate and the reactors.

* A cooling system was also installed at the upper part of the setup, enabling to cool down the vapors of monomer during the reaction. It was composed of a glass tube, positioned on top of each reactor and wrapped with a transparent hose in which circulated glycerol, at a temperature of to 3.5 °C, from a cooling thermostat bath.

* The mechanical agitator (Figure A1c) was designed with double propellers and an anchor to better scrape the reacting mixture from the walls and bottom of the reactors. In fact, since the mixture of GTR with PS became quite sticky during the polymerization, it was important to make sure that reactor content would remain under mixing throughout the polymerization, without forming an inverse bell shape with the agitator spinning in void in the middle of it. The rotation speed was fixed at 30 rpm. A pulley was fixed on the top of each stirring axe and a belt system made the 6 stirring axes rotate simultaneously (Figure A1b). The system was designed to keep the 6 axes parallel, thus avoiding stirrers from scratching and damaging the reactor walls.

* Nitrogen inerting, before the reactions, was also implemented via specifically designed inlets on the seals of the glass tubes and a dedicated nitrogen feeding network.
<figure><img src="image_5.png" /><figcaption>**Figure 1.** (a) Side view of the experimental medium-sized system. (b) Top view of the aluminum plate. (c) Side view of the aluminum plate from the red arrow.</figcaption></figure>

All reactors were first filled with the desired amounts of GTR. The respective mixtures of monomer and initiator were preliminary prepared in beakers and then introduced in the corresponding reactors directly onto the GTR. An initial blend with the mechanical agitator was performed manually to ensure a good distribution of the mixture of monomer and initiator within GTR particles. The experimental system was then mounted as shown in Figure 1a. Reactors were sealed and purged with nitrogen while initiating the cooling system and the mechanical agitation. The heating plate was turned on only once the cooling thermostat had reached its set-point temperature. During the polymerization, a temperature measurement was taken inside each reactor every 15 min during the first hour, then every hour until the end of the reaction. For these measurements, agitation was momentarily paused. In general, mixtures with low GTR content displayed higher temperature variations during the initial stages of the reaction but, eventually, all reactors reached a steady reaction temperature.

At the end of the polymerization, the reactors were removed from the experimental system and cooled-down in an ice bath to room temperature. The full contents of the reactors were removed, weighed and then transferred in vacuum oven at 40 °C until constant mass. Throughout the different steps of the experimental procedure, special care was taken to measure eventual mass losses of reactants and products (i.e., due to their transfer between recipients or due to evaporation). The uncertainty related to the weighing equipment was also considered. The conversion of styrene during the polymerization reaction $X_{styrene}$ was determined on the basis of an expression similar to Equation (1):

$$X_{styrene} = \frac{m_{i,styrene} - m_{f,styrene} - \Delta}{m_{i,styrene}} \quad (2)$$

where, $\Delta$ accounts for the total loss of styrene, from the preparation steps to the final measurement, as detailed in Appendix B. Finally, the measurement uncertainty, due to the equipment used, was assessed on the basis of the principle of error propagation [43] (cf. Appendix B for details).

## 2.4. Design of Experiments

A wide range of operating conditions, for the different experimental runs, was defined during this study in order to succeed in mapping the widest domain of the feasible experimental space. These were defined in terms of three principal conditions, namely
the reaction temperature, the peroxyde initiator content and the GTR loading, as shown in Table 1. In particular, an effort was made to avoid limiting the amount of GTR to very low values, as it is often the case in other reported studies [12,36], which would also be coherent with an ultimate objective of maximizing the recycling of end-of-life tires [10] in an eventual industrial realization of the system. In this respect, a GTR loading as high as 70% was tested. On the basis of these established ranges, D-optimal designs of experiments were prepared for both experimental setups, considering a maximum reaction time of 8 h per experiment. These are depicted in Figure 2 for both systems.

However, along the experiments and in view of the obtained results, it was decided to slightly modify the initial design by introducing a few additional points to better investigate a specific experimental domain of interest. These points are shown in Figure 2, in orange for the small sized-system and in green and yellow for the medium-sized system at 90 °C and 100 °C, respectively. Inversely, some points were omitted due to extremely low reaction advancement for the adopted reaction temperature and duration (i.e., practically zero styrene conversion) or important temperature overshoot within the first minutes of the reaction. All the realized experiments, along with the measured conversion values, are grouped in Appendix C.

<table><caption>**Table 1.** Range of operating conditions.</caption><thead><tr><th>Factors</th><th>Min</th><th>Max</th></tr></thead><tbody><tr><td>Temperature</td><td>90 °C</td><td>110 °C</td></tr><tr><td>m<sub>i,BPO</sub>/m<sub>i,styrene</sub></td><td>3%</td><td>8%</td></tr><tr><td>m<sub>i,GTR</sub>/(m<sub>i,GTR</sub> + m<sub>i,styrene</sub>)</td><td>10%</td><td>70%</td></tr></tbody></table>

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Design of experiments for small- and medium-sized systems.</figcaption></figure>

## 3. ML Modeling

### 3.1. ML Algorithms

Data-driven modeling is based on the premise that patterns and trends, which are otherwise difficult to identify or extract on the basis of knowledge and/or observation, can be identified within data coming from a system or process. Among the most common data-driven approaches that are encountered in the modeling of physicochemical systems are response surface methodology (i.e., as part of design of experiments) and ML methods, the latter being particularly adapted to highly complex or multidimensional problems [35]. In the present work, ten popular supervised ML regression models were initially screened on the produced experimental data of the grafting polymerization system to identify the ones
that were more fitted to the specific problem characteristics. Indeed, each ML technique
can be more or less suitable to a modeling problem, depending on different factors such
as the quantity and nature of the data and the form of the underlying sought pattern.
During this initial screening, the parameters of the tested models were fitted to match the
experimental data, via the model training step, but without any further optimization of the
hyperparameters (HPs) of these models. The tested models belong to different categories,
including linear (linear regression (LR) and two of its regularized counterparts, namely
ridge and lasso) and non-linear ones (support vector regression (SVR), Gaussian processes
(GP) and multilayer perceptrons (MLP)), with extended application in polymer science [34].
In addition, similarity-based (k-nearest neighbors (kNN)), decision trees (DT) and ensemble
(random forest (RF) for bagging and gradient boosting (GB) for boosting) methods were
also tested.

At a second stage, the four best-performing ML models, namely GB, RF, MLP and SVR,
were subjected to HP optimization, in an attempt to further improve their performance.
Note that, it is extremely difficult to select a priori the best suited ML technique for a specific
problem, mainly due to the plethora of existing methods and to the non-deterministic
character of the outcome of the final combination of the model setup, training and HP
optimization. This is the reason why this selection is not really justified in most reported
studies. However, besides attempting a screening of different techniques and selecting
the best performing ones, as is the case here, it is also useful to attempt to identify the
characteristics of the best performing models that make them suitable to the problem
in hand. This can provide valuable conclusions and serve as future guidance in similar
problems. Accordingly, the basic principles of these four techniques are briefly explained
below and will serve in the later discussion of the results and the conclusions of the study.
Moreover, more detailed descriptions and comparative studies are also reported in the
relevant literature [44]. Finally, note that all models were developed in Python, using the
Scikit Learn library.

GB and RF are ensemble ML methods as they consist of combining the outputs of
a large number of prediction models to obtain an improved ensemble prediction. The
difference between GB and RF lies in the way the outputs are combined as well as the
construction of the constitutive prediction models. More specifically, they, respectively,
belong to boosting and bagging (or bootstrap aggregation) categories. In boosting, simple
weak estimators (in underfitting) are combined sequentially in a way that at each iteration,
a new weak learner is trained by considering the errors of the previous one (i.e., each
new learner tries to correct the errors made by its predecessor), thus reducing underfit-
ting/bias. In GB, the considered learners are DT. The principal advantages of GB are its
capacity to model complex non-linear relationships, its customization possibilities and its
generalization ability. However, it suffers from large computation cost. The main HPs of
GB are generally the number of estimators (or learners) and the learning rate. Other HPs
such as the fraction of features and samples considered in each tree can also introduce
some randomness which has been proven to ameliorate generalization performance and
calculation time. More details can be found in [45-49].

Inversely, in bagging, several well performing prediction models (in overfitting) are
trained in parallel on a bootstrap sample of the data set (i.e., a subset of the train data
set that is randomly drawn to train an individual prediction model; this subset is then
replaced so that all subsets are drawn from the same data distribution). The individual
predictions are then combined to give one final prediction (majority vote if classification
problem, average value if regression problem) which reduces the overfitting/variance of
the individual models. In RF, the prediction models are DT. However, the main particularity
of RF method lies in that each DT is build only on a random subset of features. The random
subsets of data and features reduce the possible correlations between trees, creating a
diversity of individual prediction models (specific to particular conditions), and thus
resulting in better accuracy and stability. Moreover, as the number of trees grows, the
generalization error converges which makes RF particularly robust against overfitting.
Other advantages of RF are its low sensitivity to HP values and its capacity to evaluate
the importance of input features. As RF is based on DT, it also benefits from the latter's
simplicity and interpretability. More details about RF method can be found in [44,50–54].
The main HPs of RF are the number of trees and the number of random input features and
samples considered for each tree.

MLP, which are a specific type of artificial neural networks (ANN), consist of several layers containing the so-called neurons: an input layer, containing the inputs of the problem, an output layer, containing the outputs of the problem, and one or several intermediate hidden layers. The neurons present in the different layers are interconnected and associated with a weight. The information exchange is made possible thanks to the neurons which transform their input (i.e., the sum of all the inputs arriving from the previous layer multiplied by their corresponding weights plus a bias term) to an output by means of an activation function. Weights are iteratively adjusted during the learning process by minimizing the error between the predicted and the true outputs. A great number of HPs need to be set (and eventually optimized) in MLP, such as the number of hidden layers, the number of neurons in each hidden layer, the activation function, the learning rate and the method to adjust the weights, which renders this step especially time-consuming. On the other hand, MLP presents a powerful capacity to model complex problems as it can approximate any linear or nonlinear mathematical function. On the counterpart, the training of a MLP presents several local minima and is prone to overfitting and lack of generalization, as the complexity and the size of the network architecture increases, especially for small data sets. Furthermore, the structure of the MLP itself makes it harder to interpret compared to RF. More details about ANN and MLP can be found in [55–62].

SVR is the counterpart of support vector machines (SVM) classification algorithm for regression problems and is widely used due to its good generalization ability, high prediction accuracy and robustness during the training process (convex problem), even in presence of small, high dimensional and noisy data sets [63–66]. In SVM, the goal is to find the best hyperplane that will separate the data correctly into classes, which is performed by maximizing the margin (distance between hyperplane and closest point, also called support vector) in both sides of the hyperplane and the number of data classified correctly. When data is not linearly separable, kernel functions can be employed to re-describe them in another higher dimensional space where they will be linearly separable, thus improving classification performance. More details can be found in [63,64,67–70]. The principle of SVR is also based on finding the best hyperplane. However, the goal here is to have the maximum number of points located inside the narrowest margin (or epsilon tube), the latter containing the hyperplane which represents the regression function. The main HPs are the regularization parameter C (to penalize the highest deviations outside the margin), epsilon (radius of the epsilon tube) and the kernel. Other HPs, more specific to the kernel, can also be optimized.

3.2. *ML Procedure*

In any data-driven modeling study, the data pretreatment represents the most crucial
and time-consuming step, as the quality of data (besides quantity) will largely determine
the performance of the model (i.e., widely known as “garbage in, garbage out” principle). In
this work, data pretreatment consisted essentially in keeping the most reliable experimental
measurements (i.e., eliminating the ones in which important losses or measurement errors
were observed). The final data set, that was used for the training of all ML models,
is composed of a total of 69 experiments, as shown in Table C1. These were shuffled
and partitioned into train and test sets according to a 5-fold cross-validation procedure
(Figure 3a). The latter enabled to evaluate the variability of the performance subject to
partitioning, while ensuring that each sample was present in one of the test sets. As
part of the data pretreatment process, the considered inputs need also to be standardized
(i.e., rescaled to present a common mean and standard deviation) as means to ensuring that
their eventual difference in magnitude will not bias the model. In the present work, a mean
value of 0 and a standard deviation of 1 were used for the normalization, while this was applied solely on the train set to avoid data leakage (i.e., assure the complete independence of the train and test sets).

For the HP optimization step, a grid search cross validation method was employed to exhaustively test all possible combinations. This was motivated also by the relatively low number of models/HPs to be tested. In a different case, other HP optimization strategies, such as genetic algorithms or random search, could be employed, depending on the problem characteristics and the type of ML algorithms [71]. To fine-tune the boundaries and levels of the considered HPs, a prior manual investigation of the impact of each individual HP on the model performance was initially carried out. The different types of HPs for each model, along with their respective optimization domain, are detailed in Table 2.

A nested cross-validation scheme was adopted in order to limit the introduction of bias in combined HP tuning and model selection steps. More specifically, this procedure consists of two nested cross-validation schemes, namely an external one based on a train/test split of the whole data set and an internal one based on a train/validation split of each training set of the outer loop. The outer loop is used for model training (i.e., estimation of its parameters) and model selection, while the inner one aims in HP optimization. In this work, 5-fold outer and inner cross-validations were employed (Figure 3b). Even if this procedure is more time-consuming than using a common k-fold for both model training/selection and HP optimization, it produces more robust and unbiased performance estimates against overfitting and is highly suitable for small data sets [72,73].

Among the different metrics that are typically employed to evaluate the performance of ML models, the root mean squared error (RMSE), mean absolute error (MAE), determination coefficient (R²), max and median error, are the most commonly encountered ones. In this work, RMSE was selected as the principal model evaluation metric as it penalizes specifically the largest errors, which are to be avoided if the model has to be applied for industrial use. However, in the results section, R² and MAE values are also reported in order to provide a more complete overview of the performance of the developed ML models.

<table><caption>**Table 2.** HPs tested for GB, RF, MLP and SVR models.</caption><thead><tr><th>ML Algorithm</th><th>HPs</th><th>Values</th></tr></thead><tbody><tr><td>GB</td><td>n_estimators<br>learning_rate<br>max_features<br>minSamples_Leaf<br>subsample</td><td>[50, 100, 150]<br>[0.1]<br>['sqrt', 'log2']<br>[1, 5, 10, 15]<br>[1/5, 2/5, 3/5, 4/5, 1]</td></tr><tr><td>RF</td><td>n_estimators<br>max_features<br>maxSamples</td><td>[50, 100, 150]<br>['sqrt', 'log2']<br>[None] (no bootstrap: all train samples)</td></tr><tr><td>MLP</td><td>activation<br>hiddenLayerSizes<br>solver<br>learning_rate_init<br>max_iter</td><td>['relu']<br>1 hidden layer: [(i)] with i = 25, 50, 75, 100, 125, 150<br>2 hidden layers: [(i, j)] with i, j = 5, 10, 15, 20, 25, 30, 50<br>['lbfgs', 'adam']<br>[0.001, 0.005, 0.01, 0.04, 0.07]<br>[800]</td></tr><tr><td>SVR</td><td>kernel<br>C<br>epsilon</td><td>['rbf']<br>[0.5, 1, 1.5, 2, 3, 4, 5, 6, 8, 10]<br>[0.01, 0.05, 0.1]</td></tr></tbody></table>
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Data partitioning procedures: (a) 5-fold cross-validation procedure for ML models screening without HP optimization. (b) Nested cross-validation with 5-fold outer and 5-fold inner cross-validation procedures for combined model selection and HP optimization.</figcaption></figure>

# 4. Results and Discussion

## 4.1. *Experimental Results*

All the results for the small-sized and medium-sized systems are presented in detail, respectively, in Tables C1 and C2 and are also publicly accessible in electronic format (cf. Data Availability Statement). In the following, the behavior of the two systems in terms of temperature control, polymerization time and effect of GTR and BPO is discussed, on the basis of the obtained results. In all the presented graphs, error bars represent standard deviation (i.e., whenever an experiment has been repeated several times). The average absolute uncertainty in the determination of monomer conversion was estimated (cf. Appendix C) to be equal to an average of 0.13% wt and 0.29% wt for the small-sized and the medium-sized systems, respectively. This was primarily based on the uncertainty of the equipment used (i.e., ±0.1 mg of the precision balance). At the same time, the average standard deviation of three measurement repetitions for a total of seven samples was found to be equal to 0.3 mg. Finally, other sources of introduction of random error to the measurements, such as the balance drift, influenced by air stream, temperature or humidity, were not quantified.

During the experiments, the temperature inside the reactors was investigated to ensure the reliability of the data for latter modeling. In the small-sized system, a thermocouple was kept inside three reference tubes (corresponding to GTR/(GTR + styrene) contents of 10%, 40% and 70% wt, the ratio BPO/styrene being set at 5.5% wt) during the whole reaction for the three temperatures of interest. An example of the recorded temperature evolution, for 90 °C, is shown in Figure 4 and in Table 3. Despite an initial overshoot of the reaction temperature, which is more pronounced for mixtures with lower GTR content, the set-point temperature was reached in the first 10 min of reaction, within a margin of ±2 °C in all cases. In the medium-sized system, it was impossible to follow the evolution of the internal temperature on a permanent basis due to the form of the mechanical agitator that covered the reactor radially and throughout the complete mixture height. Accordingly, at specific time-intervals, the agitation was stopped promptly and thermocouples were introduced into the reactors through the lateral orifice for temperature measurement. In this case, the temperature overshoot was significantly more pronounced, reaching even 56 °C at certain conditions (see Table C2 for all temperature measurement results). Note that, the measured temperature was also affected by the positioning of the thermocouple inside the reactor, a phenomenon that was not observed in the small-sized system. The highest temperature values were observed below the agitator, near the bottom of the reactor, which are the ones that are reported here. These observations validate the ones made during a previous work [28] in a larger scale reactor and are evidence of the scale-up difficulty of the system. They are also indicative of the necessity to accompany published experimental measurements in similar systems with details about the achieved temperature homogeneity and the eventually observed heat-transfer limitations. In view of these elements, in the
present work, 110 °C experiments were limited for the medium-sized system while in the subsequent ML modeling treatment, only small-sized system data were considered as isothermal experiments.

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Polymerization temperature for small-sized system for different GTR/(GTR + styrene) ratios (% wt) (90 °C, BPO/styrene 5.5% wt).</figcaption></figure>

<table><caption>**Table 3.** Temperature control in the small-sized system (BPO/styrene = 5.5% wt).</caption><thead><tr><th rowspan="2">T:</th><th colspan="3">90 °C</th><th colspan="3">100 °C</th><th colspan="3">110 °C</th></tr><tr><th>10%</th><th>40%</th><th>70%</th><th>10%</th><th>40%</th><th>70%</th><th>10%</th><th>40%</th><th>70%</th></tr></thead><tbody><tr><td>GTR/(GTR + Styrene):</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Mean T after 1 h 30</td><td>89</td><td>89</td><td>89</td><td>99</td><td>99</td><td>100</td><td>106</td><td>106</td><td>109</td></tr><tr><td>Min T after 1 h 30</td><td>88</td><td>87</td><td>88</td><td>98</td><td>98</td><td>98</td><td>104</td><td>105</td><td>108</td></tr><tr><td>Max T after 1 h 30</td><td>90</td><td>90</td><td>91</td><td>101</td><td>101</td><td>101</td><td>108</td><td>108</td><td>111</td></tr><tr><td>t (min) to reach T ± 2 °C</td><td>3</td><td>3</td><td>4</td><td>2</td><td>3</td><td>3</td><td>2</td><td>-</td><td>5</td></tr></tbody></table>

Another parameter of interest, for this polymerization system, is the reaction time that, depending on the conditions, may extend up to >30 h until complete styrene or BPO consumption [12,36]. In the present work, the maximum reaction time was limited to 8 h due to technical constraints. However, the experiments carried out in the small-sized system, under the different reaction conditions set in the experimental design, were performed also for shorter duration, namely for 2 h and 5 h (see Table C1). This allowed monitoring the temporal evolution of the system and including time as a factor in the subsequent ML modeling study.

The evolution of styrene conversion, with respect to the mixture composition in GTR and BPO, is presented in Figures 5 and 6, for a reaction temperature of 90 °C and 100 °C, respectively. In these figures, continuous curves correspond to the small-sized system (denoted as S) while dashed curves are for the medium-sized system (denoted as M). The overall trend, observed for both systems and in both temperatures, is completely consistent with previous observations [27–29] where an increase in the amount of GTR in the system displayed a negative effect on monomer conversion. At the same time, a rather expected positive impact of the amount of BPO on the final styrene conversion is also observed. The inhibiting effect of GTR has been attributed to the reaction between the dissociated benzoyl acid radicals and the unpaired electrons on the surface of carbon black, creating new active sites which are either inert or have lower reactivity than the BPO primary radicals, thus limiting styrene conversion. Another possible mechanism, acting synergetically towards the same direction, could be the redox reaction between BPO and the oxide groups of the
carbon black surface, which would be favored at a lower reaction temperatures due to its lower activation energy than the thermal decomposition of BPO [27].

It is also worth noting that the effect of GTR content on the value of monomer conversion is not linear, with respect to the GTR content, but is more pronounced in the range of 10–40% wt loadings. For example, at 90 °C (i.e., Figure 5), the final monomer conversion after 8 h of reaction drops from around 80% to only 10% when passing from 10% to 40% wt in GTR content. However, half of this decrease is already evident at 15% GTR content. This phenomenon had already been identified in previous studies [27,29] and a critical GTR concentration had been considered, in the range between 10% and 50% wt of GTR content, beyond which the observed inhibiting effects are gradually attenuated.

The same trend is observed in both small-sized and medium-sized systems. At the same time, there is a slight disagreement between the results of the two systems, under similar reaction conditions, with the small-sized system displaying higher styrene conversions at lower GTR/(GTR + styrene) ratios (i.e., 10% and 15%) than the medium sized one. This tendency gradually inverses as the GTR ratio increases (i.e., up to 70% wt of GTR). This difference could be attributed to the aforementioned unavoidable mass losses that occur during the different steps in the medium-sized system, especially during the transfer of the reactor content to the recipients for mass measurement and evaporation in the vacuum oven (e.g., non-recovered material from the reactor walls and agitator). Although an effort was made to meticulously monitor these losses and account for them during the calculations (cf. Δ₂ in Equations (2) and (B.1)), it was assumed that the total losses in mass is proportionally distributed among styrene and GTR, according to their initially introduced amount in the reactor. In this sense, for low GTR/(GTR + styrene) ratios, a relatively higher percentage of the measured mass loss will be attributed to styrene, thus leading to the calculation of a higher value for monomer conversion than in a case where the same mass loss would correspond to a higher GTR content. A different consideration of the distribution of losses (e.g., an equal distribution between GTR and styrene) would have artificially decreased this observed deviation between the two systems. In any case, it remains technically impossible to verify this assumption.

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Comparison of small- (S) and medium- (M) sized systems at 90 °C.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 6.** Comparison of small- (S) and medium- (M) sized systems at 100 °C.</figcaption></figure>

Figures 7 and 8 depict the combined effect of GTR content and reaction temperature, for a constant amount of initial BPO concentration, on the measured styrene conversion. Besides the negative effect of GTR content, it is seen that the reaction temperature displays an overall positive effect by increasing styrene conversion. It is also observed that the effect of temperature becomes less profound at higher GTR loadings. These observations are completely inline with the previously reported measurements [27] as well as with the proposed hypothesis of redox reaction between BPO and the oxide groups of the carbon black surface that would be favored at lower reaction temperatures.

<figure><img src="image_11.png" /><figcaption>**Figure 7.** Effect of GTR/(GTR + styrene) and temperature for small-sized system and BPO/styrene = 5.5% (dashed curves correspond to trend lines).</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Figure 8.** Effect of GTR/(GTR + styrene) and temperature for medium-sized system and BPO/styrene = 5.5% (dashed curves correspond to trend lines).</figcaption></figure>

## 4.2. ML Results

Following the pretreatment of the measured experimental data, as explained in Section 3, 10 different supervised-learning regression algorithms were screened and compared, prior to any further optimization of their HPs. In Figure 9, the parity plots of the split corresponding to the lowest RMSE, for the test data set, are presented for each model. In addition, comparative bar plots between all tested models, on the basis of the adopted performance criteria (i.e., $R^2$, RMSE, MAE), are depicted in Figures 10 and 11 for the train and test data sets, respectively. The error bars correspond to the standard deviation, calculated on the basis of the 5-fold learning procedure. The exact values of all the metrics are also provided in Table D1. The best performing models, identified in terms of their test data set performance metrics were found to be the ensemble models (GB and RF), DT, MLP and SVR. The corresponding test $R^2$ values for these models varied between 0.927 ± 0.032 and 0.762 ± 0.066.

These models, identified as the top-performing ones in this case, are among the most commonly employed in similar reported studies [35]. By their construction, ensemble models can generalize very well to new data as they are composed of several models, each trained under specific data conditions. DT models also perform particularly well, but typically show lower performance than ensemble models as they can be considered as a sub-case of RF. MLP networks are notorious for their capacity in modeling highly non-linear response surfaces but suffer from low generalization and a need for extensive data-sets. On the other hand, SVR presents a more robust training performance and perform well with limited data. It would thus seem as these models displayed better performance than the rest due to a combination of characteristics of the problem that would partially match their specific advantages. Accordingly, GB, RF, MLP and SVR were selected for the HP optimization phase, in order to investigate whether further significant improvement could be attained in their performance, differentiating one (or more) of them as a clear choice for the modeling of this system. Note that, the selection of these models was not based strictly on their performance but also on a choice to further test models that are structurally different (i.e., ensemble models vs. SVR and ANN) and quite popular in such applications. Note also that some of them, such as SVR and MLP, are typically considered to be more sensitive to their HP values than others (e.g., RF). Finally, the low values of standard deviation for the selected models are indicative of their robustness against the train/test splits, which allows generalizing the conclusions about their applicability with higher confidence.
<figure><img src="image_13.png" /><figcaption>Figure 9. Parity plots of the ten screened models (without HP optimization).</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>**Figure 10.** Train performances of the ten screened ML models (without HP optimization).</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>**Figure 11.** Test performances of the ten screened ML models (without HP optimization).</figcaption></figure>

Concerning the rather poorly performing models, it seems that the experimental data could not be described by linear models and/or by the use of a reduced number of input features (i.e., as is the case in strongly regularized linear models), thus leading to evident underfitting. This becomes obvious in the case of Lasso, which significantly reduces the number of input features (e.g., as opposed to Ridge). An optimization of the regularization parameter in Ridge and Lasso models could probably help to improve the performance results and/or to reduce underfitting but this was not further pursued in this work. On the contrary, GP model performs perfectly on the training data set but fails to generalize in the test set, displaying clear evidence of severe overfitting. Similarly, an optimization of the GP kernel could have been explored as means to improve its performance, just as an optimized value of the number of neighbors might had improved the performance of the kNN model and render it more interesting for this problem.

The four best-performing ML models (GB, RF, MLP and SVR) were subsequently subjected to HP optimization using a "nested cross-validation" procedure, instead of a simple k-fold validation (i.e., as the one applied before HP optimization). The obtained optimized HPs are presented in Table 4, for each of the 5 train/test splits of the 5-fold outer cross-validation, and are completed with the average computation time for one split. Note that certain HPs of some of the models, such as the number of estimators in RF, the hidden
layer sizes in MLP and the regularization parameter C in SVR, vary depending on the split, while the rest are not affected. It is also observed that the train/test split displays a more or less significant effect on the model performance, depending on the method. For example, Figure 12 shows that the test $R^2$ of SVR with optimized HPs varies from 0.7 to 0.9, depending on the split. However, the overall net improvement of SVR and MLP models, after HP optimization, is also clearly visible in Figure 13, where train and test data are better aligned on the $y = x$ bisector of the parity plots.

<table><caption>**Table 4.** HPs optimized for each train/test split, for GB, RF, MLP and SVR.</caption><thead><tr><th>ML Algorithm</th><th>HPs</th><th>Split 1</th><th>Split 2</th><th>Split 3</th><th>Split 4</th><th>Split 5</th><th>Time (1 Split)</th></tr></thead><tbody><tr><td rowspan="5">GB</td><td>n_estimators</td><td>150</td><td>150</td><td>150</td><td>150</td><td>150</td><td rowspan="5">3.8 s</td></tr><tr><td>learning_rate</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.1</td></tr><tr><td>max_features</td><td>'log2'</td><td>'log2'</td><td>'log2'</td><td>'log2'</td><td>'log2'</td></tr><tr><td>min Samples leaf</td><td>5</td><td>5</td><td>5</td><td>1</td><td>1</td></tr><tr><td>subsample</td><td>1</td><td>1</td><td>0.8</td><td>0.6</td><td>0.4</td></tr><tr><td rowspan="3">RF</td><td>n_estimators</td><td>50</td><td>150</td><td>150</td><td>50</td><td>50</td><td rowspan="3">1.2 s</td></tr><tr><td>max_features</td><td>'log2'</td><td>'log2'</td><td>'log2'</td><td>'log2'</td><td>'log2'</td></tr><tr><td>max Samples</td><td>None</td><td>None</td><td>None</td><td>None</td><td>None</td></tr><tr><td rowspan="5">MLP</td><td>activation</td><td>'relu'</td><td>'relu'</td><td>'relu'</td><td>'relu'</td><td>'relu'</td><td rowspan="5">27.5 s</td></tr><tr><td>hidden layer sizes</td><td>(15, 10)</td><td>(50, 5)</td><td>(5, 5)</td><td>(5, 15)</td><td>(5, 25)</td></tr><tr><td>solver</td><td>'lbfgs'</td><td>'lbfgs'</td><td>'lbfgs'</td><td>'lbfgs'</td><td>'lbfgs'</td></tr><tr><td>learning_rate_init</td><td>0.001</td><td>0.001</td><td>0.001</td><td>0.001</td><td>0.001</td></tr><tr><td>max_iter</td><td>800</td><td>800</td><td>800</td><td>800</td><td>800</td></tr><tr><td rowspan="3">SVR</td><td>kernel</td><td>'rbf'</td><td>'rbf'</td><td>'rbf'</td><td>'rbf'</td><td>'rbf'</td><td rowspan="3">0.7 s</td></tr><tr><td>C</td><td>1</td><td>4</td><td>8</td><td>1</td><td>2</td></tr><tr><td>epsilon</td><td>0.01</td><td>0.01</td><td>0.01</td><td>0.01</td><td>0.01</td></tr></tbody></table>

<figure><img src="image_16.png" /><figcaption>**Figure 12.** Comparison of test performances of GB, RF, MLP and SVR with default vs. optimized HPs.</figcaption></figure>
<figure><img src="image_17.png" /><figcaption>**Figure 13.** Comparative parity plots (with default HPs on left side, and with optimized HPs on right side).</figcaption></figure>

Once these optimized values were obtained, they were used as part of the 5-fold outer cross-validation, to train the models on train set and evaluate them on test set for each train/test split. The test performances ($R^2$, RMSE and MAE) with default (i.e., prior to HP optimization) and optimized HPs are compared in Figure 12, where full bars represent the performances with default HPs while hatched ones are for performances with optimized HPs. The error bars represent the standard deviations calculated on the basis of the five different train/test splits. The exact values of all the metrics are also provided in Table D2. Note that the final results are subject to the selected HP optimization method and metric. Accordingly, the selection of a method different than grid search cross validation or the optimization on the basis of another metric than RMSE, would have probably resulted in different values for the optimized HPs and for the model performances. A net improvement is observed in the SVR and MLP model performances on the test data set, whilst GB and RF
display no improvement but rather a slight degradation in their performance. Moreover, as discussed previously, RF models are known for showing low sensitivity to their HPs values. The improvement of SVR and MLP models after HP optimization is also visible in Figure 13, where train and test data are better aligned on the *y* = *x* bisector of the parity plots. Finally, note that the test performance of the models displays a higher variation than their respective train performance (cf. Table D2), which can be considered as further justification of the need to employ a k-fold validation procedure rather than evaluating the performance of the models on the basis of any random split.

Several options could be envisioned to further improve the prediction performances of the models. First of all, it is clear that the limited number of samples prevents the models from learning from more data and thus improving their generalization performance. Indeed, as discussed in the introduction of this work, data-driven techniques require sufficient amount and quality of data. Although it is not evident to define a priori what “sufficient” means, especially since different methods have different requirements in terms of the amount of data, relative studies have shed light to the interplay that exists between the amount of data, used for the training of ML models, and their precision [74]. Accordingly, in the case of problems of limited data, one option would be to introduce a crude estimation of the target property in the feature space to enhance the accuracy of ML models.

At the same time, it is clear that fundamental experimental studies of physicochemical systems or materials cannot compete, in terms of data availability, with problems related to the fields of computer science, medical records or industrial data analysis. In this sense, it remains crucial that the experimental procedure is rigorous, involving identification and monitoring of all the steps that introduce uncertainties to the measurements. The impact of the experimental uncertainties on ML models performance can be subsequently evaluated and taken into account in the model development phase. For example, the incorporation of known data uncertainty measures on input and/or target data has been previously investigated for ANN, by using least-squares-inspired methods, and resulted in better generalization capacity [75–78]. In the present work, this possibility was exploited via a weighing procedure of the data used in the training steps of SVR, RF and GB by the estimated uncertainties on monomer conversion measurements. In particular, the normalized inverses of the uncertainties were used to weigh the samples in the regression functions (i.e., through the option sample_weight). However, only a marginal improvement was observed solely for SVR, corresponding to an increase in $R^2$ from 0.806 to 0.817, thus this option was not further pursued. The comparison of the test performances between the default case (all samples have same weight) and the case with weighted samples, after HP optimization, is provided in Figure E1.

Other parameters related to the ML model development procedure, such as the ratio for the partitioning of data into train and test set, the data standardization method and the random data shuffling and model initialization, could also influence the final model performance. The influence of these parameters in the performances of GB, RF, MLP and SVR, prior to HP optimization, was also studied in this work. The results are provided in Figures E2–E4, respectively, for the three previously described parameters. They showed that SVR can be slightly improved by using MinMax scaling instead of standard scaling. At the same time, MLP showed a net degradation of its performance with MinMax and robust scaling while the other two methods did not display any significant difference among the three scaling techniques. The partitioning ratio 90/10 for train and test sets, corresponding to 10-fold cross-validation, despite improving the average test performance was accompanied by higher standard deviations. Indeed, more data was available to train the models in this case but, at the same time, the model was evaluated on a very small test set which resulted in a broader error distribution. Finally, the random state, which controls the sequence of the pseudo-random number generator, was also found to display an effect on the model performance. As discussed previously, this effect can only be counterbalanced by the implementation of a k-fold cross validation procedure for the model training and HP optimization steps. However, it is important to remain aware of this random fluctuation
of the model performance, especially given the limited amount of available data. Finally, another option for potentially improving the model performances would be by integrating their data-driven character with existing knowledge on the system. One way to achieve this would be to couple them with the previously developed mechanistic model [29] in an attempt to exploit simultaneously the advantages of each one of the two approaches in a unique hybrid model. This remains a perspective for future work.

## 5. Conclusions

The study of the kinetics of the styrene-GTR radical graft polymerization represents a great challenge due to the complexity of the prevailing kinetic mechanisms as well as the nature of GTR. In addition, the existing knowledge or controlling capacity of these aspects remains limited. In this respect, a rigorous experimental study and a ML modeling approach were combined in this work in an attempt to overcome these limitations and relate operating conditions with a key performance index, namely styrene conversion.

The particularity of the adopted experimental approach lied in the attention paid in ensuring data reliability and quantity, which are crucial in any ML approach. For this purpose, two laboratory-scale experimental setups were implemented to evaluate the reactive medium thermal and compositional homogeneity and investigate scale-up effects. The experimental procedure for both systems was designed so that data can be produced rapidly, while increasing data reliability by eliminating the uncertainty related to sampling for analyses. At the same time, all the potential sources of uncertainty, such as the mass loss along the different steps of the process and the precision of the experimental equipment, were carefully identified and monitored. The analysis of the results on both experimental scales revealed that the effects of GTR, BPO and temperature on styrene conversion were coherent with previous studies. However, the heterogeneous nature of the system created difficulties in the mixing and the estimation of the losses for the medium-sized system, which should not be neglected in view of an eventual further scaling-up.

The most reliable experimental data, containing 69 experiments from the small-sized system, was subsequently subjected to the screening of diverse supervised-learning regression ML models and the optimization of their HPs to identify the best-performing ones. These were identified to be GB, MLP and RF with test $R^2$ of, respectively, 0.91 ± 0.04, 0.90 ± 0.04 and 0.89 ± 0.05. In view of the small available data set, a nested cross-validation scheme was employed to account for variance in both validation and test sets, and the results revealed the impact of the random data partitioning on the performance of the developed models. Finally, the effect of additional parameters, such as the scaling method and the number of folds, as well as the integration of the experimental uncertainties in the learning procedure, were exploited as means to improve the predictions of the models.

From the perspective of upscaling the studied system towards industrial application, further research on how to enhance the homogeneity of temperature and composition is more than necessary, especially at higher GTR loadings. The latter are indeed necessary to increase the recycling of used tires. To limit the costs of the process, the use of commercially available reagents and GTR without further purification was preferred in this work. However, the presence of additives, stabilizers or impurities affects the polymerization course, and the variable GTR composition represents a major part of this added uncertainty. This complexity encourages further use of ML tools, ideally in combination with already developed phenomenological models in hybrid approaches, to quickly map and, if possible, understand the relationship between operating conditions and final product properties at different scales. Developing ML models compatible with small data sets and capable of integrating/predicting uncertainties is therefore an interesting direction for future work, as well as the extension to other molecular properties. For the latter, experimental procedures enabling the collection of sufficient data of good quality are therefore required. Finally, an effort to share experimental data and accompany each experimental result with precise and detailed conditions is crucial, as demonstrated in this work.
**Author Contributions:** Conceptualization, all; methodology, C.T. and D.M.; literature review, C.T.; experimental study, C.T.; design and implementation of the medium-sized system, R.L. and S.H.; data curation and modeling, C.T.; writing—original draft preparation, C.T.; writing—review and editing, C.T. and D.M.; materials, S.H.; supervision, D.M. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by MESRI (Ministère de l'Enseignement supérieur, de la Recherche et de l'Innovation), France.

**Data Availability Statement:** Experimental data and ML code are available at https://github.com/cindyytrinhh/Machine-Learning-modeling-of-styrene-Ground-Tire-Rubber-radical-graft-polymerization, published on 20 February 2023.

**Acknowledgments:** The authors acknowledge the "Service Etudes et Réalisation Mécanique" of LRGP as well as the "Service de soufflage du verre" of University of Lorraine who, respectively, designed and built the medium-sized system and the glass parts. The authors also acknowledge the following masters students who contributed in setting up and collecting experimental data: Redouanne Chenna, Justine Guerry, Bahia Lamari, Betty Soulet and Jérémy Lanoizelée.

**Conflicts of Interest:** The authors declare no conflict of interest.

## Abbreviations

The following abbreviations are used in this manuscript:

<table><tr><td>ANN</td><td>Artificial neural networks</td></tr><tr><td>BPO</td><td>Benzoyl peroxide</td></tr><tr><td>DT</td><td>Decision trees</td></tr><tr><td>GB</td><td>Gradient boosting</td></tr><tr><td>GP</td><td>Gaussian processes</td></tr><tr><td>GTR</td><td>Ground tire rubber</td></tr><tr><td>HP</td><td>Hyperparameter</td></tr><tr><td>kNN</td><td>k-nearest neighbors</td></tr><tr><td>LR</td><td>Linear regression</td></tr><tr><td>MAE</td><td>Mean absolute error</td></tr><tr><td>ML</td><td>Machine learning</td></tr><tr><td>MLP</td><td>Multilayer perceptrons</td></tr><tr><td>PS</td><td>Polystyrene</td></tr><tr><td>RF</td><td>Random forest</td></tr><tr><td>RMSE</td><td>Root mean squared error</td></tr><tr><td>SVM</td><td>Support vector machines</td></tr><tr><td>SVR</td><td>Support vector regression</td></tr></table>

## Appendix A. Additional Photos of the Experimental Medium-Sized System

<figure><img src="image_18.png" /><figcaption>Figure A1. Medium-sized experimental system. (a) Whole system. (b) Pulleys-belt. (c) Mechanical agitator.</figcaption></figure>
## Appendix B. Evaluation of the Experimental Uncertainties

During the different steps of experimental procedure, especially in the medium-sized setup, several losses of reactants take place. These need to be identified and quantified for the correct calculation of styrene conversion (i.e., denoted as $\Delta$ in Equation (2)). These losses are identified as follows:

$$ \Delta = \sum_{i=1}^{3} \Delta_{i} \qquad (B.1) $$

where:

* $\Delta_1$ is the loss of styrene during polymerization reaction due to evaporation;

* $\Delta_2$ is the loss of styrene during the transfer of the reactor content to the gravimetry cup at the end of the polymerization;

* $\Delta_3$ is the styrene remaining in the glass tube at the end of the polymerization reaction.

All the above losses were identified by meticulous weighing of the used equipment before and after each individual step, using a Mettler Toledo XP504 precision balance. During these measurements, the ratio of styrene to GTR was considered to be identical and constant, as defined in the recipe of each experiment and measured during the initial material introduction step.

According to [43], the uncertainty of a function $f = f(x_1, x_2, \dots, x_n)$ can be expressed in terms of the uncertainties $u_{x_i}$ of each one of its variables $x_i$ and partial derivatives $\frac{\partial f}{\partial x_i}$ (for $i = 1 \dots n$), as:

$$ u_f = \sqrt{\sum_{i=1}^{n} \left[ \frac{\partial f}{\partial x_i} u_{x_i} \right]^2} \qquad (B.2) $$

The application of Equation (B.2) to $X_{styrene}$ and $\Delta$, defined in Equations (2) and (B.1), respectively, results in:

$$ \begin{aligned} u_{X_{styrene}} &= \sqrt{\left[ \frac{\partial X_{styrene}}{\partial m_{i,styrene}} u_{m_{i,styrene}} \right]^2 + \left[ \frac{\partial X_{styrene}}{\partial m_{f,styrene}} u_{m_{f,styrene}} \right]^2 + \left[ \frac{\partial X_{styrene}}{\partial \Delta} u_{\Delta} \right]^2} \\ &= \sqrt{\left[ \frac{(m_{f,styrene} + \Delta) * u_{m_{i,styrene}}}{m_{i,styrene}^2} \right]^2 + \left[ \frac{u_{m_{f,styrene}}}{m_{i,styrene}} \right]^2 + \left[ \frac{u_{\Delta}}{m_{i,styrene}} \right]^2} \end{aligned} \qquad (B.3) $$

$$ \begin{aligned} u_{\Delta} &= \sqrt{u_{\Delta_1}^2 + u_{\Delta_2}^2 + u_{\Delta_3}^2} \end{aligned} \qquad (B.4) $$
Appendix C. Detailed Experimental Data for the Small- and Medium-Sized Systems

<table><caption>Table C1. Experimental results for the small-sized system.</caption><thead><tr><th>N°</th><th>Time, h</th><th>Temperature theo., °C</th><th>Temperature exp., °C</th><th>GTR/(GTR+Styrene) theo., %wt</th><th>GTR/(GTR+Styrene) exp., %wt</th><th>BPO/Styrene theo., %wt</th><th>BPO/Styrene exp., %wt</th><th>Conversion, %wt</th><th>Balance-Related uncertainty, %wt</th></tr></thead><tbody><tr><td>1</td><td>2</td><td>90.0</td><td>89.0</td><td>10.00%</td><td>10.23%</td><td>5.50%</td><td>5.49%</td><td>47.57%</td><td>0.13%</td></tr><tr><td>2</td><td>2</td><td>90.0</td><td>89.0</td><td>40.00%</td><td>40.35%</td><td>5.50%</td><td>5.49%</td><td>0.00%</td><td>0.13%</td></tr><tr><td>3</td><td>2</td><td>90.0</td><td>89.0</td><td>70.00%</td><td>69.94%</td><td>5.50%</td><td>5.49%</td><td>0.00%</td><td>0.07%</td></tr><tr><td>4</td><td>5</td><td>90.0</td><td>89.0</td><td>10.00%</td><td>9.79%</td><td>5.50%</td><td>5.49%</td><td>58.41%</td><td>0.13%</td></tr><tr><td>5</td><td>5</td><td>90.0</td><td>89.0</td><td>40.00%</td><td>40.06%</td><td>5.50%</td><td>5.49%</td><td>1.49%</td><td>0.11%</td></tr><tr><td>6</td><td>5</td><td>90.0</td><td>89.0</td><td>70.00%</td><td>69.87%</td><td>5.50%</td><td>5.49%</td><td>0.00%</td><td>0.08%</td></tr><tr><td>7</td><td>8</td><td>90.0</td><td>89.0</td><td>10.00%</td><td>9.79%</td><td>3.00%</td><td>2.99%</td><td>28.97%</td><td>0.08%</td></tr><tr><td>8</td><td>8</td><td>90.0</td><td>89.0</td><td>10.00%</td><td>10.28%</td><td>5.50%</td><td>5.55%</td><td>67.84%</td><td>0.10%</td></tr><tr><td>9</td><td>8</td><td>90.0</td><td>89.0</td><td>10.00%</td><td>10.16%</td><td>5.50%</td><td>5.49%</td><td>53.13%</td><td>0.08%</td></tr><tr><td>10</td><td>8</td><td>90.0</td><td>89.0</td><td>10.00%</td><td>10.00%</td><td>8.00%</td><td>8.03%</td><td>86.53%</td><td>0.08%</td></tr><tr><td>11</td><td>8</td><td>90.0</td><td>89.0</td><td>15.00%</td><td>15.03%</td><td>5.50%</td><td>5.49%</td><td>22.05%</td><td>0.08%</td></tr><tr><td>12</td><td>8</td><td>90.0</td><td>89.0</td><td>15.00%</td><td>14.98%</td><td>8.00%</td><td>8.03%</td><td>35.60%</td><td>0.08%</td></tr><tr><td>13</td><td>8</td><td>90.0</td><td>89.0</td><td>25.00%</td><td>24.99%</td><td>5.50%</td><td>5.51%</td><td>8.94%</td><td>0.09%</td></tr><tr><td>14</td><td>8</td><td>90.0</td><td>89.0</td><td>40.00%</td><td>40.01%</td><td>3.00%</td><td>2.99%</td><td>4.36%</td><td>0.15%</td></tr><tr><td>15</td><td>8</td><td>90.0</td><td>89.0</td><td>40.00%</td><td>41.65%</td><td>5.50%</td><td>5.77%</td><td>3.42%</td><td>0.11%</td></tr><tr><td>16</td><td>8</td><td>90.0</td><td>89.0</td><td>40.00%</td><td>40.13%</td><td>5.50%</td><td>5.49%</td><td>4.25%</td><td>0.11%</td></tr><tr><td>17</td><td>8</td><td>90.0</td><td>89.0</td><td>40.00%</td><td>40.11%</td><td>8.00%</td><td>8.00%</td><td>8.81%</td><td>0.11%</td></tr><tr><td>18</td><td>8</td><td>90.0</td><td>89.0</td><td>55.00%</td><td>54.96%</td><td>5.50%</td><td>5.51%</td><td>4.86%</td><td>0.14%</td></tr><tr><td>19</td><td>8</td><td>90.0</td><td>89.0</td><td>70.00%</td><td>69.97%</td><td>3.00%</td><td>2.99%</td><td>0.00%</td><td>0.31%</td></tr><tr><td>20</td><td>8</td><td>90.0</td><td>89.0</td><td>70.00%</td><td>71.68%</td><td>5.50%</td><td>5.94%</td><td>2.45%</td><td>0.13%</td></tr><tr><td>21</td><td>8</td><td>90.0</td><td>89.0</td><td>70.00%</td><td>69.95%</td><td>5.50%</td><td>5.49%</td><td>0.00%</td><td>0.13%</td></tr><tr><td>22</td><td>8</td><td>90.0</td><td>89.0</td><td>70.00%</td><td>69.95%</td><td>8.00%</td><td>8.00%</td><td>8.55%</td><td>0.21%</td></tr><tr><td>23</td><td>2</td><td>100.0</td><td>99.3</td><td>10.00%</td><td>10.42%</td><td>5.50%</td><td>5.75%</td><td>55.66%</td><td>0.21%</td></tr><tr><td>24</td><td>2</td><td>100.0</td><td>99.3</td><td>40.00%</td><td>40.20%</td><td>5.50%</td><td>5.50%</td><td>4.18%</td><td>0.11%</td></tr><tr><td>25</td><td>2</td><td>100.0</td><td>99.3</td><td>70.00%</td><td>70.15%</td><td>5.50%</td><td>5.50%</td><td>3.00%</td><td>0.11%</td></tr><tr><td>26</td><td>5</td><td>100.0</td><td>99.3</td><td>10.00%</td><td>10.02%</td><td>5.50%</td><td>5.50%</td><td>59.96%</td><td>0.21%</td></tr><tr><td>27</td><td>5</td><td>100.0</td><td>99.3</td><td>40.00%</td><td>40.20%</td><td>5.50%</td><td>5.50%</td><td>4.72%</td><td>0.11%</td></tr><tr><td>28</td><td>5</td><td>100.0</td><td>99.3</td><td>70.00%</td><td>70.05%</td><td>5.50%</td><td>5.50%</td><td>3.13%</td><td>0.08%</td></tr><tr><td>29</td><td>8</td><td>100.0</td><td>99.3</td><td>10.00%</td><td>9.98%</td><td>3.00%</td><td>3.01%</td><td>54.50%</td><td>0.29%</td></tr><tr><td>30</td><td>8</td><td>100.0</td><td>99.3</td><td>10.00%</td><td>10.44%</td><td>5.50%</td><td>5.79%</td><td>71.82%</td><td>0.15%</td></tr><tr><td>31</td><td>8</td><td>100.0</td><td>99.3</td><td>10.00%</td><td>9.99%</td><td>5.50%</td><td>5.50%</td><td>77.12%</td><td>0.11%</td></tr><tr><td>32</td><td>8</td><td>100.0</td><td>99.3</td><td>10.00%</td><td>9.98%</td><td>8.00%</td><td>7.99%</td><td>80.12%</td><td>0.13%</td></tr><tr><td>33</td><td>8</td><td>100.0</td><td>99.3</td><td>15.00%</td><td>15.01%</td><td>8.00%</td><td>7.99%</td><td>52.21%</td><td>0.08%</td></tr><tr><td>34</td><td>8</td><td>100.0</td><td>99.3</td><td>25.00%</td><td>25.08%</td><td>5.50%</td><td>5.48%</td><td>0.00%</td><td>0.08%</td></tr><tr><td>35</td><td>8</td><td>100.0</td><td>99.3</td><td>40.00%</td><td>40.05%</td><td>3.00%</td><td>3.01%</td><td>3.99%</td><td>0.08%</td></tr><tr><td>36</td><td>8</td><td>100.0</td><td>99.3</td><td>40.00%</td><td>40.59%</td><td>5.50%</td><td>5.70%</td><td>5.94%</td><td>0.08%</td></tr><tr><td>37</td><td>8</td><td>100.0</td><td>99.3</td><td>40.00%</td><td>40.25%</td><td>5.50%</td><td>5.50%</td><td>7.93%</td><td>0.08%</td></tr><tr><td>38</td><td>8</td><td>100.0</td><td>99.3</td><td>40.00%</td><td>40.08%</td><td>8.00%</td><td>8.01%</td><td>6.35%</td><td>0.13%</td></tr><tr><td>39</td><td>8</td><td>100.0</td><td>99.3</td><td>55.00%</td><td>55.01%</td><td>5.50%</td><td>5.48%</td><td>0.00%</td><td>0.11%</td></tr><tr><td>40</td><td>8</td><td>100.0</td><td>99.3</td><td>70.00%</td><td>70.02%</td><td>3.00%</td><td>3.01%</td><td>4.87%</td><td>0.11%</td></tr><tr><td>41</td><td>8</td><td>100.0</td><td>99.3</td><td>70.00%</td><td>70.41%</td><td>5.50%</td><td>5.75%</td><td>4.15%</td><td>0.13%</td></tr><tr><td>42</td><td>8</td><td>100.0</td><td>99.3</td><td>70.00%</td><td>70.53%</td><td>5.50%</td><td>5.50%</td><td>0.00%</td><td>0.21%</td></tr><tr><td>43</td><td>8</td><td>100.0</td><td>99.3</td><td>70.00%</td><td>69.99%</td><td>8.00%</td><td>8.01%</td><td>7.58%</td><td>0.21%</td></tr><tr><td>44</td><td>2</td><td>110.0</td><td>107.0</td><td>10.00%</td><td>10.06%</td><td>5.50%</td><td>5.50%</td><td>65.74%</td><td>0.21%</td></tr><tr><td>45</td><td>2</td><td>110.0</td><td>107.0</td><td>40.00%</td><td>40.06%</td><td>5.50%</td><td>5.50%</td><td>4.25%</td><td>0.11%</td></tr><tr><td>46</td><td>2</td><td>110.0</td><td>107.0</td><td>70.00%</td><td>69.96%</td><td>5.50%</td><td>5.50%</td><td>1.98%</td><td>0.08%</td></tr><tr><td>47</td><td>5</td><td>110.0</td><td>107.0</td><td>10.00%</td><td>9.96%</td><td>5.50%</td><td>5.50%</td><td>80.32%</td><td>0.21%</td></tr><tr><td>48</td><td>5</td><td>110.0</td><td>107.0</td><td>40.00%</td><td>40.07%</td><td>5.50%</td><td>5.50%</td><td>9.50%</td><td>0.11%</td></tr><tr><td>49</td><td>5</td><td>110.0</td><td>107.0</td><td>70.00%</td><td>70.12%</td><td>5.50%</td><td>5.50%</td><td>2.72%</td><td>0.08%</td></tr><tr><td>50</td><td>8</td><td>110.0</td><td>107.0</td><td>10.00%</td><td>10.33%</td><td>3.00%</td><td>3.00%</td><td>79.60%</td><td>0.29%</td></tr><tr><td>51</td><td>8</td><td>110.0</td><td>107.0</td><td>10.00%</td><td>10.10%</td><td>5.50%</td><td>5.58%</td><td>93.55%</td><td>0.14%</td></tr><tr><td>52</td><td>8</td><td>110.0</td><td>107.0</td><td>10.00%</td><td>10.01%</td><td>5.50%</td><td>5.50%</td><td>86.15%</td><td>0.11%</td></tr><tr><td>53</td><td>8</td><td>110.0</td><td>107.0</td><td>10.00%</td><td>10.02%</td><td>5.50%</td><td>5.49%</td><td>83.63%</td><td>0.21%</td></tr><tr><td>54</td><td>8</td><td>110.0</td><td>107.0</td><td>10.00%</td><td>10.72%</td><td>8.00%</td><td>8.00%</td><td>89.71%</td><td>0.11%</td></tr><tr><td>55</td><td>8</td><td>110.0</td><td>107.0</td><td>15.00%</td><td>15.04%</td><td>5.50%</td><td>5.49%</td><td>64.39%</td><td>0.08%</td></tr><tr><td>56</td><td>8</td><td>110.0</td><td>107.0</td><td>15.00%</td><td>15.07%</td><td>8.00%</td><td>8.00%</td><td>73.23%</td><td>0.08%</td></tr><tr><td>57</td><td>8</td><td>110.0</td><td>107.0</td><td>25.00%</td><td>25.08%</td><td>5.50%</td><td>5.49%</td><td>25.18%</td><td>0.13%</td></tr><tr><td>58</td><td>8</td><td>110.0</td><td>107.0</td><td>32.50%</td><td>32.56%</td><td>5.50%</td><td>5.50%</td><td>15.51%</td><td>0.10%</td></tr><tr><td>59</td><td>8</td><td>110.0</td><td>107.0</td><td>40.00%</td><td>40.01%</td><td>3.00%</td><td>3.00%</td><td>7.86%</td><td>0.08%</td></tr><tr><td>60</td><td>8</td><td>110.0</td><td>107.0</td><td>40.00%</td><td>39.59%</td><td>5.50%</td><td>5.42%</td><td>10.27%</td><td>0.08%</td></tr><tr><td>61</td><td>8</td><td>110.0</td><td>107.0</td><td>40.00%</td><td>40.07%</td><td>5.50%</td><td>5.50%</td><td>14.69%</td><td>0.09%</td></tr><tr><td>62</td><td>8</td><td>110.0</td><td>107.0</td><td>40.00%</td><td>40.07%</td><td>5.50%</td><td>5.49%</td><td>10.11%</td><td>0.09%</td></tr><tr><td>63</td><td>8</td><td>110.0</td><td>107.0</td><td>40.00%</td><td>39.87%</td><td>8.00%</td><td>8.00%</td><td>9.33%</td><td>0.11%</td></tr><tr><td>64</td><td>8</td><td>110.0</td><td>107.0</td><td>55.00%</td><td>55.04%</td><td>5.50%</td><td>5.49%</td><td>8.92%</td><td>0.11%</td></tr><tr><td>65</td><td>8</td><td>110.0</td><td>107.0</td><td>70.00%</td><td>69.96%</td><td>3.00%</td><td>3.00%</td><td>8.62%</td><td>0.11%</td></tr><tr><td>66</td><td>8</td><td>110.0</td><td>107.0</td><td>70.00%</td><td>70.32%</td><td>5.50%</td><td>5.63%</td><td>7.22%</td><td>0.14%</td></tr><tr><td>67</td><td>8</td><td>110.0</td><td>107.0</td><td>70.00%</td><td>69.99%</td><td>5.50%</td><td>5.50%</td><td>11.42%</td><td>0.21%</td></tr><tr><td>68</td><td>8</td><td>110.0</td><td>107.0</td><td>70.00%</td><td>70.00%</td><td>5.50%</td><td>5.49%</td><td>8.12%</td><td>0.21%</td></tr><tr><td>69</td><td>8</td><td>110.0</td><td>107.0</td><td>70.00%</td><td>70.08%</td><td>8.00%</td><td>8.00%</td><td>9.45%</td><td>0.21%</td></tr></tbody></table>
<table><caption>Table C2. Experimental results for the medium-sized system.</caption><thead><tr><th>N°</th><th>Time, h</th><th>Temperature theo., °C</th><th>Temperature exp., °C</th><th>GTR/(GTR+Styrene) theo., %wt</th><th>GTR/(GTR+Styrene) exp., %wt</th><th>BPO/Styrene theo., %wt</th><th>BPO/Styrene exp., %wt</th><th>Conversion, %wt</th><th>Balance-Related Uncertainty, %wt</th></tr></thead><tbody><tr><td>1</td><td>90</td><td>94.6</td><td>&gt;98.16</td><td>10.00%</td><td>10.20%</td><td>3.00%</td><td>3.04%</td><td>18.43%</td><td>0.12%</td></tr><tr><td>2</td><td>90</td><td>93.4</td><td>&gt;98.45</td><td>10.00%</td><td>10.15%</td><td>3.00%</td><td>3.01%</td><td>14.40%</td><td>0.12%</td></tr><tr><td>3</td><td>90</td><td>93.8</td><td>&gt;109.47</td><td>10.00%</td><td>10.19%</td><td>5.50%</td><td>5.58%</td><td>49.20%</td><td>0.12%</td></tr><tr><td>4</td><td>90</td><td>na</td><td>na</td><td>10.00%</td><td>10.66%</td><td>5.50%</td><td>5.89%</td><td>37.88%</td><td>0.14%</td></tr><tr><td>5</td><td>90</td><td>92.2</td><td>&gt;93</td><td>10.00%</td><td>10.15%</td><td>5.50%</td><td>5.57%</td><td>49.93%</td><td>0.13%</td></tr><tr><td>6</td><td>90</td><td>92.9</td><td>&gt;146.30</td><td>10.00%</td><td>10.19%</td><td>8.00%</td><td>8.08%</td><td>72.52%</td><td>0.12%</td></tr><tr><td>7</td><td>90</td><td>95.7</td><td>&gt;136.95</td><td>10.00%</td><td>10.28%</td><td>8.00%</td><td>8.16%</td><td>72.53%</td><td>0.12%</td></tr><tr><td>8</td><td>90</td><td>91.9</td><td>&gt;94</td><td>15.00%</td><td>15.24%</td><td>5.50%</td><td>5.54%</td><td>10.40%</td><td>0.13%</td></tr><tr><td>9</td><td>90</td><td>92.3</td><td>&gt;109.77</td><td>15.00%</td><td>15.25%</td><td>8.00%</td><td>8.07%</td><td>24.20%</td><td>0.13%</td></tr><tr><td>10</td><td>90</td><td>96.7</td><td>&gt;115.47</td><td>25.00%</td><td>25.77%</td><td>4.25%</td><td>4.37%</td><td>6.11%</td><td>0.21%</td></tr><tr><td>11</td><td>90</td><td>95</td><td>&gt;97.03</td><td>25.00%</td><td>25.60%</td><td>5.00%</td><td>5.13%</td><td>5.75%</td><td>0.21%</td></tr><tr><td>12</td><td>90</td><td>93.9</td><td>&gt;95.97</td><td>25.00%</td><td>25.38%</td><td>6.00%</td><td>6.09%</td><td>5.03%</td><td>0.21%</td></tr><tr><td>13</td><td>90</td><td>93.9</td><td>&gt;94.72</td><td>25.00%</td><td>25.35%</td><td>6.75%</td><td>6.83%</td><td>6.23%</td><td>0.21%</td></tr><tr><td>14</td><td>90</td><td>na</td><td>na</td><td>30.00%</td><td>30.47%</td><td>5.50%</td><td>5.56%</td><td>2.75%</td><td>0.22%</td></tr><tr><td>15</td><td>90</td><td>90.5</td><td>&gt;91</td><td>30.00%</td><td>30.37%</td><td>5.50%</td><td>5.51%</td><td>2.64%</td><td>0.22%</td></tr><tr><td>16</td><td>90</td><td>92.5</td><td>&gt;93.30</td><td>40.00%</td><td>40.47%</td><td>5.50%</td><td>5.54%</td><td>3.54%</td><td>0.35%</td></tr><tr><td>17</td><td>90</td><td>90.7</td><td>&gt;95.85</td><td>40.00%</td><td>40.73%</td><td>8.00%</td><td>8.16%</td><td>5.39%</td><td>0.35%</td></tr><tr><td>18</td><td>90</td><td>na</td><td>na</td><td>70.00%</td><td>72.46%</td><td>5.50%</td><td>5.88%</td><td>6.99%</td><td>1.19%</td></tr><tr><td>19</td><td>90</td><td>89.2</td><td>&gt;92</td><td>70.00%</td><td>73.21%</td><td>5.50%</td><td>6.32%</td><td>5.87%</td><td>1.22%</td></tr><tr><td>20</td><td>100</td><td>103.6</td><td>&gt;105.69</td><td>10.00%</td><td>10.18%</td><td>3.00%</td><td>3.04%</td><td>45.04%</td><td>0.12%</td></tr><tr><td>21</td><td>100</td><td>104.4</td><td>&gt;123.65</td><td>10.00%</td><td>10.20%</td><td>5.50%</td><td>5.56%</td><td>69.92%</td><td>0.12%</td></tr><tr><td>22</td><td>100</td><td>103.9</td><td>&gt;129.84</td><td>10.00%</td><td>10.62%</td><td>8.00%</td><td>8.45%</td><td>91.82%</td><td>0.13%</td></tr><tr><td>23</td><td>100</td><td>102.6</td><td>&gt;143.30</td><td>10.00%</td><td>10.32%</td><td>8.00%</td><td>8.20%</td><td>78.12%</td><td>0.13%</td></tr><tr><td>24</td><td>100</td><td>102.2</td><td>&gt;118.83</td><td>10.00%</td><td>10.16%</td><td>8.00%</td><td>8.08%</td><td>83.37%</td><td>0.13%</td></tr><tr><td>25</td><td>100</td><td>103</td><td>&gt;104.38</td><td>15.00%</td><td>15.20%</td><td>5.50%</td><td>5.49%</td><td>19.42%</td><td>0.13%</td></tr><tr><td>26</td><td>100</td><td>na</td><td>na</td><td>15.00%</td><td>15.22%</td><td>5.50%</td><td>5.53%</td><td>18.57%</td><td>0.13%</td></tr><tr><td>27</td><td>100</td><td>101.7</td><td>&gt;109.85</td><td>15.00%</td><td>15.23%</td><td>8.00%</td><td>8.08%</td><td>36.72%</td><td>0.13%</td></tr><tr><td>28</td><td>100</td><td>na</td><td>na</td><td>20.00%</td><td>21.86%</td><td>5.50%</td><td>6.11%</td><td>6.61%</td><td>0.14%</td></tr><tr><td>29</td><td>100</td><td>104</td><td>&gt;105</td><td>20.00%</td><td>20.24%</td><td>5.50%</td><td>5.57%</td><td>9.81%</td><td>0.13%</td></tr><tr><td>30</td><td>100</td><td>102</td><td>&gt;103.96</td><td>20.00%</td><td>20.28%</td><td>8.00%</td><td>8.07%</td><td>12.43%</td><td>0.14%</td></tr><tr><td>31</td><td>100</td><td>106</td><td>&gt;106.85</td><td>25.00%</td><td>25.52%</td><td>4.25%</td><td>4.31%</td><td>16.18%</td><td>0.20%</td></tr><tr><td>32</td><td>100</td><td>103.7</td><td>&gt;104.90</td><td>25.00%</td><td>25.12%</td><td>5.00%</td><td>4.99%</td><td>9.70%</td><td>0.21%</td></tr><tr><td>33</td><td>100</td><td>104.3</td><td>&gt;105.29</td><td>25.00%</td><td>25.17%</td><td>5.00%</td><td>5.03%</td><td>12.52%</td><td>0.20%</td></tr><tr><td>34</td><td>100</td><td>100.2</td><td>&gt;106.49</td><td>25.00%</td><td>25.34%</td><td>6.00%</td><td>6.05%</td><td>6.07%</td><td>0.18%</td></tr><tr><td>35</td><td>100</td><td>104.5</td><td>&gt;105.65</td><td>25.00%</td><td>25.36%</td><td>6.00%</td><td>6.06%</td><td>8.68%</td><td>0.18%</td></tr><tr><td>36</td><td>100</td><td>106.7</td><td>&gt;115.21</td><td>25.00%</td><td>26.05%</td><td>6.75%</td><td>7.06%</td><td>16.49%</td><td>0.21%</td></tr><tr><td>37</td><td>100</td><td>na</td><td>na</td><td>30.00%</td><td>30.05%</td><td>3.00%</td><td>3.01%</td><td>4.13%</td><td>0.20%</td></tr><tr><td>38</td><td>100</td><td>101.4</td><td>&gt;103</td><td>30.00%</td><td>30.33%</td><td>3.00%</td><td>3.00%</td><td>3.85%</td><td>0.21%</td></tr><tr><td>39</td><td>100</td><td>na</td><td>na</td><td>30.00%</td><td>30.30%</td><td>5.50%</td><td>5.53%</td><td>5.02%</td><td>0.21%</td></tr><tr><td>40</td><td>100</td><td>101.1</td><td>&gt;103</td><td>30.00%</td><td>30.40%</td><td>5.50%</td><td>5.56%</td><td>4.58%</td><td>0.21%</td></tr><tr><td>41</td><td>100</td><td>na</td><td>na</td><td>30.00%</td><td>30.24%</td><td>8.00%</td><td>8.05%</td><td>4.77%</td><td>0.21%</td></tr><tr><td>42</td><td>100</td><td>100.4</td><td>&gt;103</td><td>30.00%</td><td>30.72%</td><td>8.00%</td><td>8.19%</td><td>6.15%</td><td>0.21%</td></tr><tr><td>43</td><td>100</td><td>103.3</td><td>&gt;105.66</td><td>32.50%</td><td>33.32%</td><td>5.50%</td><td>5.63%</td><td>6.43%</td><td>0.27%</td></tr><tr><td>44</td><td>100</td><td>103.3</td><td>&gt;104.14</td><td>32.50%</td><td>32.73%</td><td>6.25%</td><td>6.25%</td><td>7.49%</td><td>0.26%</td></tr><tr><td>45</td><td>100</td><td>101.9</td><td>&gt;105.63</td><td>35.00%</td><td>36.84%</td><td>4.00%</td><td>4.30%</td><td>5.30%</td><td>0.28%</td></tr><tr><td>46</td><td>100</td><td>103.8</td><td>&gt;105.61</td><td>35.00%</td><td>35.92%</td><td>7.00%</td><td>7.19%</td><td>7.74%</td><td>0.28%</td></tr><tr><td>47</td><td>100</td><td>101.6</td><td>&gt;105.65</td><td>40.00%</td><td>40.43%</td><td>3.00%</td><td>3.03%</td><td>5.48%</td><td>0.35%</td></tr><tr><td>48</td><td>100</td><td>na</td><td>na</td><td>40.00%</td><td>40.75%</td><td>5.50%</td><td>5.63%</td><td>4.37%</td><td>0.34%</td></tr><tr><td>49</td><td>100</td><td>99.6</td><td>&gt;101</td><td>40.00%</td><td>40.62%</td><td>5.50%</td><td>5.56%</td><td>4.20%</td><td>0.33%</td></tr><tr><td>50</td><td>100</td><td>104.5</td><td>&gt;106.32</td><td>40.00%</td><td>41.69%</td><td>6.25%</td><td>6.55%</td><td>8.69%</td><td>0.36%</td></tr><tr><td>51</td><td>100</td><td>102</td><td>&gt;102.70</td><td>40.00%</td><td>40.90%</td><td>8.00%</td><td>8.13%</td><td>7.79%</td><td>0.36%</td></tr><tr><td>52</td><td>100</td><td>na</td><td>na</td><td>70.00%</td><td>72.44%</td><td>5.50%</td><td>6.05%</td><td>8.07%</td><td>1.14%</td></tr><tr><td>53</td><td>100</td><td>101.2</td><td>&gt;102</td><td>70.00%</td><td>72.69%</td><td>5.50%</td><td>6.00%</td><td>8.43%</td><td>1.15%</td></tr><tr><td>54</td><td>100</td><td>99.4</td><td>&gt;103.70</td><td>70.00%</td><td>73.56%</td><td>8.00%</td><td>9.16%</td><td>12.16%</td><td>1.30%</td></tr><tr><td>55</td><td>110</td><td>114.8</td><td>&gt;115.40</td><td>25.00%</td><td>25.25%</td><td>4.25%</td><td>4.18%</td><td>22.11%</td><td>0.20%</td></tr><tr><td>56</td><td>110</td><td>112.8</td><td>&gt;113.52</td><td>25.00%</td><td>25.25%</td><td>5.00%</td><td>5.04%</td><td>18.17%</td><td>0.20%</td></tr><tr><td>57</td><td>110</td><td>112.9</td><td>&gt;113.60</td><td>25.00%</td><td>25.17%</td><td>6.75%</td><td>6.78%</td><td>20.68%</td><td>0.20%</td></tr><tr><td>58</td><td>110</td><td>112.9</td><td>&gt;113.74</td><td>32.50%</td><td>32.94%</td><td>5.50%</td><td>5.52%</td><td>10.84%</td><td>0.26%</td></tr><tr><td>59</td><td>110</td><td>113.3</td><td>&gt;114.45</td><td>32.50%</td><td>33.11%</td><td>6.25%</td><td>6.29%</td><td>11.44%</td><td>0.26%</td></tr><tr><td>60</td><td>110</td><td>113.7</td><td>&gt;114.96</td><td>35.00%</td><td>35.91%</td><td>7.00%</td><td>7.16%</td><td>13.46%</td><td>0.33%</td></tr></tbody></table>
## Appendix D. Detailed ML Results

<table><caption>**Table D1.** Detailed metrics ($R^2$, RMSE and MAE) obtained for the ten screened models with default HPs.</caption><thead><tr><th>ML Algorithm</th><th>R<sup>2</sup> Train</th><th>Test</th><th>RMSE Train</th><th>Test</th><th>MAE Train</th><th>Test</th></tr></thead><tbody><tr><td>LR</td><td>0.697 ± 0.016</td><td>0.619 ± 0.107</td><td>0.168 ± 0.003</td><td>0.177 ± 0.013</td><td>0.148 ± 0.003</td><td>0.156 ± 0.011</td></tr><tr><td>Ridge</td><td>0.697 ± 0.016</td><td>0.621 ± 0.101</td><td>0.168 ± 0.003</td><td>0.177 ± 0.013</td><td>0.148 ± 0.003</td><td>0.156 ± 0.012</td></tr><tr><td>Lasso</td><td>0.000 ± 0.000</td><td>-0.096 ± 0.069</td><td>0.306 ± 0.010</td><td>0.308 ± 0.038</td><td>0.270 ± 0.011</td><td>0.273 ± 0.018</td></tr><tr><td>SVR</td><td>0.905 ± 0.013</td><td>0.762 ± 0.066</td><td>0.094 ± 0.009</td><td>0.142 ± 0.029</td><td>0.078 ± 0.007</td><td>0.111 ± 0.021</td></tr><tr><td>GP</td><td>1.000 ± 0.000</td><td>-1736 ± 1572</td><td>0.000 ± 0.000</td><td>8.644 ± 6.903</td><td>0.000 ± 0.000</td><td>3.592 ± 2.829</td></tr><tr><td>kNN</td><td>0.769 ± 0.043</td><td>0.554 ± 0.206</td><td>0.146 ± 0.012</td><td>0.193 ± 0.064</td><td>0.104 ± 0.011</td><td>0.137 ± 0.047</td></tr><tr><td>DT</td><td>1.000 ± 0.000</td><td>0.843 ± 0.050</td><td>0.000 ± 0.000</td><td>0.116 ± 0.027</td><td>0.000 ± 0.000</td><td>0.075 ± 0.017</td></tr><tr><td>RF</td><td>0.986 ± 0.002</td><td>0.899 ± 0.042</td><td>0.036 ± 0.003</td><td>0.094 ± 0.031</td><td>0.023 ± 0.002</td><td>0.063 ± 0.022</td></tr><tr><td>GB</td><td>1.000 ± 0.000</td><td>0.927 ± 0.032</td><td>0.007 ± 0.000</td><td>0.079 ± 0.026</td><td>0.005 ± 0.000</td><td>0.054 ± 0.016</td></tr><tr><td>MLP</td><td>0.925 ± 0.012</td><td>0.804 ± 0.071</td><td>0.083 ± 0.008</td><td>0.127 ± 0.030</td><td>0.061 ± 0.006</td><td>0.096 ± 0.025</td></tr></tbody></table>

<table><caption>**Table D2.** Detailed metrics ($R^2$, RMSE and MAE) obtained for SVR, RF, GB and MLP with optimized HPs.</caption><thead><tr><th></th><th>SVR</th><th>RF</th><th>GB</th><th>MLP</th></tr></thead><tbody><tr><td>R<sup>2</sup> train</td><td>0.959 ± 0.019</td><td>0.985 ± 0.002</td><td>0.996 ± 0.002</td><td>0.992 ± 0.008</td></tr><tr><td>R<sup>2</sup> test</td><td>0.806 ± 0.119</td><td>0.878 ± 0.049</td><td>0.905 ± 0.043</td><td>0.898 ± 0.045</td></tr><tr><td>RMSE train</td><td>0.061 ± 0.015</td><td>0.038 ± 0.003</td><td>0.018 ± 0.005</td><td>0.025 ± 0.012</td></tr><tr><td>RMSE test</td><td>0.124 ± 0.048</td><td>0.103 ± 0.033</td><td>0.088 ± 0.025</td><td>0.090 ± 0.017</td></tr><tr><td>MAE train</td><td>0.031 ± 0.007</td><td>0.025 ± 0.002</td><td>0.014 ± 0.003</td><td>0.016 ± 0.009</td></tr><tr><td>MAE test</td><td>0.090 ± 0.035</td><td>0.070 ± 0.023</td><td>0.059 ± 0.014</td><td>0.062 ± 0.016</td></tr></tbody></table>

## Appendix E. Ways of Improvement for the ML Model

<figure><img src="image_19.png" /><figcaption>**Figure E1.** Comparison of test performances of GB, RF and SVR with optimized HPs, without weighted vs. with weighted samples.</figcaption></figure>
<figure><img src="image_20.png" /><figcaption>**Figure E2.** Scaler impact on test RMSE, without HP optimization.</figcaption></figure>

<figure><img src="image_21.png" /><figcaption>**Figure E3.** Partitioning train/test impact on test RMSE, without HP optimization.</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>**Figure E4.** Random state impact on test RMSE, without HP optimization.</figcaption></figure>
References

1. Abbas-Abadi, M.S.; Kusenberg, M.; Shirazi, H.M.; Goshayeshi, B.; Van Geem, K.M. Towards full recyclability of end-of-life tires: Challenges and opportunities. *J. Clean. Prod.* **2022**, *374*, 134036. [CrossRef]

2. Hejna, A.; Korol, J.; Przybysz-Romatowska, M.; Zedler, Ł.; Chmielnicki, B.; Formela, K. Waste tire rubber as low-cost and environmentally-friendly modifier in thermoset polymers—A review. *Waste Manag.* **2020**, *108*, 106–118. [CrossRef] [PubMed]

3. Fazli, A.; Rodrigue, D. Recycling waste tires into ground tire rubber (Gtr)/rubber compounds: A review. *J. Compos. Sci.* **2020**, *4*, 103. [CrossRef]

4. Ramarad, S.; Khalid, M.; Ratnam, C.T.; Chuah, A.L.; Rashmi, W. Waste tire rubber in polymer blends: A review on the evolution, properties and future. *Prog. Mater. Sci.* **2015**, *72*, 100–140. [CrossRef]

5. Araujo-Morera, J.; Verdugo-Manzanares, R.; González, S.; Verdejo, R.; Lopez-Manchado, M.A.; Santana, M.H. On the use of mechano-chemically modified ground tire rubber (Gtr) as recycled and sustainable filler in styrene-butadiene rubber (sbr) composites. *J. Compos. Sci.* **2021**, *5*, 68. [CrossRef]

6. Phiri, M.M.; Phiri, M.J.; Formela, K.; Wang, S.; Hlangothi, S.P. Grafting and reactive extrusion technologies for compatibilization of ground tyre rubber composites: Compounding, properties, and applications. *J. Clean. Prod.* **2022**, *369*, 133084. [CrossRef]

7. Colom, X.; Cañavate, J.; Carrillo-Navarrete, F. Towards Circular Economy by the Valorization of Different Waste Subproducts through Their Incorporation in Composite Materials: Ground Tire Rubber and Chicken Feathers. *Polymers* **2022**, *14*, 1090. [CrossRef]

8. He, M.; Gu, K.; Wang, Y.; Li, Z.; Shen, Z.; Liu, S.; Wei, J. Development of high-performance thermoplastic composites based on polyurethane and ground tire rubber by in-situ synthesis. *Resour. Conserv. Recycl.* **2021**, *173*, 105713. [CrossRef]

9. Archibong, F.N.; Sanusi, O.M.; Médéric, P.; Hocine, N.A. An overview on the recycling of waste ground tyre rubbers in thermoplastic matrices: Effect of added fillers. *Resour. Conserv. Recycl.* **2021**, *175*, 105894. [CrossRef]

10. Liu, S.; Peng, Z.; Zhang, Y.; Rodrigue, D.; Wang, S. Compatibilized thermoplastic elastomers based on highly filled polyethylene with ground tire rubber. *J. Appl. Polym. Sci.* **2022**, *139*, e52999. [CrossRef]

11. Fan, P.; Lu, C. Surface Graft Copolymerization of Poly(methyl methacrylate) onto Waste Tire Rubber Powder through Ozonization. *J. Appl. Polym. Sci.* **2011**, *122*, 2262–2270. [CrossRef]

12. Coiai, S.; Passaglia, E.; Ciardelli, F.; Tirelli, D.; Peruzzotti, F.; Resmini, E. Modification of cross-linked rubber particles by free radical polymerization. *Macromol. Symp.* **2006**, *234*, 193–202. [CrossRef]

13. Sulcis, R.; Lotti, L.; Coiai, S.; Ciardelli, F.; Passaglia, E. Novel HDPE/ground tyre rubber composite materials obtained through in-situ polymerization and polymerization filling technique. *J. Appl. Polym. Sci.* **2014**, *131*, 1–13. [CrossRef]

14. Sulcis, R.; Vizza, F.; Oberhauser, W.; Ciardelli, F.; Spiniello, R.; Dintcheva, N.T.; Passaglia, E. Recycling ground tire rubber (GTR) scraps as high-impact filler of in situ produced polyketone matrix. *Polym. Adv. Technol.* **2014**, *25*, 1060–1068. [CrossRef]

15. Tsagkalias, I.S.; Vlachou, A.; Verros, G.D.; Achilias, D.S. Effect of graphene oxide or functionalized graphene oxide on the copolymerization kinetics of Styrene/n-butyl methacrylate. *Polymers* **2019**, *11*, 999. [CrossRef] [PubMed]

16. Potts, J.R.; Lee, S.H.; Alam, T.M.; An, J.; Stoller, M.D.; Piner, R.D.; Ruoff, R.S. Thermomechanical properties of chemically modified graphene/poly(methyl methacrylate) composites made by in situ polymerization. *Carbon* **2011**, *49*, 2615–2623. [CrossRef]

17. Tripathi, S.N.; Saini, P.; Gupta, D.; Choudhary, V. Electrical and mechanical properties of PMMA/reduced graphene oxide nanocomposites prepared via in situ polymerization. *J. Mater. Sci.* **2013**, *48*, 6223–6232. [CrossRef]

18. Wang, J.; Hu, H.; Wang, X.; Xu, C.; Zhang, M.; Shang, X. Preparation and Mechanical and Electrical Properties of Graphene Nanosheets–Poly(methyl methacrylate) Nanocomposites via In Situ Suspension Polymerization Jingchao. *J. Appl. Polym. Sci.* **2011**, *122*, 1866–1871. [CrossRef]

19. Feng, L.; Guan, G.; Li, C.; Zhang, D.; Xiao, Y.; Zheng, L.; Zhu, W. In situ synthesis of poly(methyl methacrylate)/graphene oxide nanocomposites using thermal-initiated and graphene oxide-initiated polymerization. *J. Macromol. Sci. Part A Pure Appl. Chem.* **2013**, *50*, 720–727. [CrossRef]

20. Michailidis, M.; Verros, G.D.; Deliyanni, E.A.; Andriotis, E.G.; Achilias, D.S. An experimental and theoretical study of butyl methacrylate in situ radical polymerization kinetics in the presence of graphene oxide nanoadditive. *J. Polym. Sci. Part A Polym. Chem.* **2017**, *55*, 1433–1441. [CrossRef]

21. Tsagkalias, I.S.; Papadopoulou, S.; Verros, G.D.; Achilias, D.S. Polymerization Kinetics of n-Butyl Methacrylate in the Presence of Graphene Oxide Prepared by Two Different Oxidation Methods with or without Functionalization. *Ind. Eng. Chem. Res.* **2018**, *57*, 2449–2460. [CrossRef]

22. Funck, A.; Kaminsky, W. Polypropylene carbon nanotube composites by in situ polymerization. *Compos. Sci. Technol.* **2007**, *67*, 906–915. [CrossRef]

23. Jiang, X.; Bin, Y.; Matsuo, M. Electrical and mechanical properties of polyimide-carbon nanotubes composites fabricated by in situ polymerization. *Polymer* **2005**, *46*, 7418–7424. [CrossRef]

24. Verros, G.D.; Achilias, D.S. Toward the development of a mathematical model for the bulk in situ radical polymerization of methyl methacrylate in the presence of nano-additives. *Can. J. Chem. Eng.* **2016**, *94*, 1783–1791. [CrossRef]

25. Yeh, J.M.; Liou, S.J.; Lai, M.C.; Chang, Y.W.; Huang, C.Y.; Chen, C.P.; Jaw, J.H.; Tsai, T.Y.; Yu, Y.H. Comparative studies of the properties of poly(methyl methacrylate)-clay nanocomposite materials prepared by in situ emulsion polymerization and solution dispersion. *J. Appl. Polym. Sci.* **2004**, *94*, 1936–1946. [CrossRef]
26. Meng, X.; Wu, H.; Storti, G.; Morbidelli, M. Effect of Dispersed Polymeric Nanoparticles on the Bulk Polymerization of Methyl Methacrylate. *Macromolecules* **2016**, *49*, 7758–7766. [CrossRef]

27. Florez, D.; Hoppe, S.; Hu, G.H.; Meimaroglou, D. Radical bulk polymerization of styrene in the presence of rubber particles from recycled tires: A kinetic study using DSC. *J. Therm. Anal. Calorim.* **2020**, *143*, 3073–3084. [CrossRef]

28. Florez Parra, D.C. Effects of the Presence of Recycled Tire Powders on the Kinetics of the Radical Polymerization of Styrene and the Properties of the Resulting Materials. Ph.D. Thesis, Université de Lorraine, Lorraine, France, 2019; pp. 1–197.

29. Meimaroglou, D.; Florez, D.; Hu, G.H. A kinetic modeling framework for the peroxide-initiated radical polymerization of styrene in the presence of rubber particles from recycled tires. *Chem. Eng. Sci.* **2022**, *248*, 117137. [CrossRef]

30. Cameron, G.G.; Qureshi, M.Y. Free Radical Grafting of Monomers to Polydienes-4. Kinetics and Mechanism of Methyl Methacrylate Grafting to Polybutadiene. *J. Polym. Sci. Part A-1 Polym. Chem.* **1980**, *18*, 3149–3161. [CrossRef]

31. Estenoz, D.A.; Meira, G.R.; Gomez, N.; Oliva, H.M. Mathematical model of a continuous industrial high-impact polystyrene process. *AIChE J.* **1998**, *44*, 427–441. [CrossRef]

32. Meira, G.R.; Luciani, C.V.; Estenoz, D.A. Continuous Bulk Process for the Production of High-Impact Polystyrene: Recent Developments in Modeling and Control. *Macromol. React. Eng.* **2007**, *1*, 25–39. [CrossRef]

33. Zhu, C.X.; Wu, Y.Y.; Figueira, F.L.; Van Steenberge, P.H.; D'hooge, D.R.; Zhou, Y.N.; Luo, Z.H. Sensitivity analysis of isothermal free radical induced grafting through application of the distribution—Numerical fractionation—Method of moments. *Chem. Eng. J.* **2022**, *444*, 136595. [CrossRef]

34. Xu, P.; Chen, H.; Li, M.; Lu, W. New Opportunity: Machine Learning for Polymer Materials Design and Discovery. *Adv. Theory Simul.* **2022**, *5*, 2100565. [CrossRef]

35. Trinh, C.; Meimaroglou, D.; Hoppe, S. Machine learning in chemical product engineering: The state of the art and a guide for newcomers. *Processes* **2021**, *9*, 1456. [CrossRef]

36. Zhang, J.L.; Chen, H.X.; Ke, C.M.; Zhou, Y.; Lu, H.Z.; Wang, D.L. Graft polymerization of styrene onto waste rubber powder and surface characterization of graft copolymer. *Polym. Bull.* **2012**, *68*, 789–801. [CrossRef]

37. Fazli, A.; Rodrigue, D. Effect of ground tire rubber (Gtr) particle size and content on the morphological and mechanical properties of recycled high-density polyethylene (rhdpe)/gtr blends. *Recycling* **2021**, *6*, 44. [CrossRef]

38. Zhang, J.; Chen, H.; Zhou, Y.; Ke, C.; Lu, H. Compatibility of waste rubber powder/polystyrene blends by the addition of styrene grafted styrene butadiene rubber copolymer: Effect on morphology and properties. *Polym. Bull.* **2013**, *70*, 2829–2841. [CrossRef]

39. Aggarwal, P.K.; Karmarkar, A.; Taqui, S.N. Grafting of Styrene Onto Cellulose. *Int. Res. J. Eng. Technol.* **2018**, *2*, 1081–1089.

40. Hejna, A.; Klein, M.; Saeb, M.R.; Formela, K. Towards understanding the role of peroxide initiators on compatibilization efficiency of thermoplastic elastomers highly filled with reclaimed GTR. *Polym. Test.* **2019**, *73*, 143–151. [CrossRef]

41. Liu, H.L.; Wang, X.P.; Jia, D.M. Recycling of waste rubber powder by mechano-chemical modification. *J. Clean. Prod.* **2020**, *245*, 118716. [CrossRef]

42. Fletes, R.C.; López, E.O.; Gudiño, P.O.; Mendizábal, E.; Núñez, R.G.; Rodrigue, D. Ground tire rubber/polyamide 6 thermoplastic elastomers produced by dry blending and compression molding. *Prog. Rubber Plast. Recycl. Technol.* **2022**, *38*, 38–55. [CrossRef]

43. Ku, H. Notes on the use of propagation of error formulas. *J. Res. Natl. Bur. Stand. Sect. C Eng. Instrum.* **1966**, *70C*, 263. [CrossRef]

44. Rodriguez-Galiano, V.; Sanchez-Castillo, M.; Chica-Olmo, M.; Chica-Rivas, M. Machine learning predictive models for mineral prospectivity: An evaluation of neural networks, random forest, regression trees and support vector machines. *Ore Geol. Rev.* **2015**, *71*, 804–818. [CrossRef]

45. Freund, Y.; Schapire, R.E. A Decision-Theoretic Generalization of On-Line Learning and an Application to Boosting. *J. Comput. Syst. Sci.* **1997**, *55*, 119–139. [CrossRef]

46. Friedman, J.H. Greedy Function Approximation: A Gradient Boosting Machine. *Ann. Stat.* **2001**, *29*, 1189–1232. [CrossRef]

47. Friedman, J.H. Stochastic gradient boosting. *Comput. Stat. Data Anal.* **2002**, *38*, 367–378. [CrossRef]

48. Natekin, A.; Knoll, A. Gradient boosting machines, a tutorial. *Front. Neurorobotics* **2013**, *7*, 21. [CrossRef]

49. Bentéjac, C.; Csörgő, A.; Martínez-Muñoz, G. A comparative analysis of gradient boosting algorithms. *Artif. Intell. Rev.* **2021**, *54*, 1937–1967. [CrossRef]

50. Breiman, L. Random Forests. *Mach. Learn.* **2001**, *45*, 5–32. [CrossRef]

51. Liaw, A.; Wiener, M. Classification and Regression by randomForest. *R News* **2002**, *2*, 18–22.

52. Zhang, C.; Ma, Y. *Ensemble Machine Learning*; Springer: New York, NY, USA, 2012. [CrossRef]

53. Svetnik, V.; Liaw, A.; Tong, C.; Culberson, J.C.; Sheridan, R.P.; Feuston, B.P. Random Forest: A Classification and Regression Tool for Compound Classification and QSAR Modeling. *J. Chem. Inf. Comput. Sci.* **2003**, *43*, 1947–1958. [CrossRef] [PubMed]

54. Segal, M.R. *Machine Learning Benchmarks and Random Forest Regression Publication Date Machine Learning Benchmarks and Random Forest Regression*; Center for Bioinformatics and Molecular Biostatistics: San Francisco, CA, USA, 2004; p. 15.

55. Krogh, A. What are artificial neural networks? *Nat. Biotechnol.* **2008**, *26*, 195–197. [CrossRef]

56. Murtagh, F. Multilayer perceptrons for classification and regression. *Neurocomputing* **1991**, *2*, 183–197. [CrossRef]

57. Bishop, C.M. Neural networks and their applications. *Rev. Sci. Instrum.* **1998**, *65*, 1803–1832. [CrossRef]

58. Gasteiger, J.; Zupan, J. Neural Networks in Chemistry. *Angew. Chem. Int. Ed. Engl.* **1993**, *32*, 503–527. [CrossRef]

59. Gardner, M.W.; Dorling, S.R. Artificial neural networks (the multilayer perceptron)—A review of applications in the atmospheric sciences. *Atmos. Environ.* **1998**, *32*, 2627–2636. [CrossRef]
60. Zhang, Z.; Friedrich, K. Artificial neural networks applied to polymer composites: A review. *Compos. Sci. Technol.* **2003**, *63*, 2029–2044. [CrossRef]

61. Paliwal, M.; Kumar, U.A. Neural networks and statistical techniques: A review of applications. *Expert Syst. Appl.* **2009**, *36*, 2–17. [CrossRef]

62. Abiodun, O.I.; Jantan, A.; Omolara, A.E.; Dada, K.V.; Mohamed, N.A.E.; Arshad, H. State-of-the-art in artificial neural network applications: A survey. *Heliyon* **2018**, *4*, e00938. [CrossRef] [PubMed]

63. Awad, M.; Khanna, R. *Efficient Learning Machines*; Apress: Berkeley, CA, USA, 2015.

64. Ge, Z.; Song, Z.; Ding, S.X.; Huang, B. Data Mining and Analytics in the Process Industry: The Role of Machine Learning. *IEEE Access* **2017**, *5*, 20590–20616. [CrossRef]

65. Zendehboudi, S.; Rezaei, N.; Lohi, A. Applications of hybrid models in chemical, petroleum, and energy systems: A systematic review. *Appl. Energy* **2018**, *228*, 2539–2566. [CrossRef]

66. Golkarnarenji, G.; Naebe, M.; Badii, K.; Milani, A.S.; Jazar, R.N.; Khayyam, H. A machine learning case study with limited data for prediction of carbon fiber mechanical properties. *Comput. Ind.* **2019**, *105*, 123–132. [CrossRef]

67. Vapnik, V.N. *The Nature of Statistical Learning Theory*; Springer: New York, NY, USA, 1995.

68. Drucker, H.; Surges, C.J.; Kaufman, L.; Smola, A.; Vapnik, V. Support vector regression machines. *Adv. Neural Inf. Process. Syst.* **1997**, *1*, 155–161.

69. Smola, A.J.; Schölkopf, B. A tutorial on support vector regression. *Stat. Comput.* **2004**, *14*, 199–222. [CrossRef]

70. Welling, M. Support Vector Regression. Available online: https://www.academia.edu/326409 (accessed on 1 December 2022).

71. Yang, L.; Shami, A. On hyperparameter optimization of machine learning algorithms: Theory and practice. *Neurocomputing* **2020**, *415*, 295–316. [CrossRef]

72. Vabalas, A.; Gowen, E.; Poliakoff, E.; Casson, A.J. Machine learning algorithm validation with a limited sample size. *PLoS ONE* **2019**, *14*, e0224365. [CrossRef]

73. Cawley, G.C.; Talbot, N.L. On over-fitting in model selection and subsequent selection bias in performance evaluation. *J. Mach. Learn. Res.* **2010**, *11*, 2079–2107.

74. Zhang, Y.; Ling, C. A strategy to apply machine learning to small datasets in materials science. *npj Comput. Mater.* **2018**, *4*, 28–33. [CrossRef]

75. Shahvandi, M.K.; Soja, B. Inclusion of data uncertainty in machine learning and its application in geodetic data science, with case studies for the prediction of Earth orientation parameters and GNSS station coordinate time series. *Adv. Space Res.* **2022**, *70*, 563–575. [CrossRef]

76. Czarnecki, W.M.; Podolak, I.T. Machine Learning with Known Input Data Uncertainty Measure. In Proceedings of the 12th IFIP TC 8 International Conference, CISIM 2013, Krakow, Poland, 25–27 September 2013; Volume 8104, pp. 379–388. [CrossRef]

77. Alizadehsani, R.; Roshanzamir, M.; Hussain, S.; Khosravi, A.; Koohestani, A.; Zangooei, M.H.; Abdar, M.; Beykikhoshk, A.; Shoeibi, A.; Zare, A.; et al. Handling of uncertainty in medical data using machine learning and probability theory techniques: A review of 30 years (1991–2020). *Ann. Oper. Res.* **2021**. [CrossRef] [PubMed]

78. Psaros, A.F.; Meng, X.; Zou, Z.; Guo, L.; Karniadakis, G.E. Uncertainty Quantification in Scientific Machine Learning: Methods, Metrics, and Comparisons. *J. Comput. Phys.* **2023**, *477*, 111902. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.