Machine Learning-Based Prediction of Polymer Properties Using Structure–Property Relationship Modeling

1. Introduction

In general, polymer materials are of primary importance for modern science and industry. They are the backbone of many diverse applications, ranging from packaging and textile materials to aerospace and biomedical devices and energy storage materials (Tamasi et al., 2022). Polymer science is a continuously developing branch of knowledge that has seen substantial progress over the last few decades, driven by advances in chemistry, materials science, and computational science (Nistane et al., 2022). The idea of SPR is a key driver of this evolution and has helped to create a better understanding of how the configurations of molecules, chains, and functional groups influence the macroscopic properties of materials (Sobuz et al., 2025). The use of polymer databases and online resources has led to considerable advances in polymer science over the last few years (Z. Zhang et al., 2022).

Polymers have been around since the beginning of the twentieth century, but it wasn’t until the second half of the twentieth century that there was a major interest in QSPRs (Quantitative Structure Property Relationships) for property prediction (Sampedro et al., 2022). Recently, with the addition of data science, this field has witnessed a significant change, enabling faster exploration of large chemical spaces (Günay et al., 2022). The emergence of high-throughput experimental techniques, along with open-source data, has enabled a significant increase in polymer research, with exponentially more polymer compounds being synthesized and investigated each year (Yakoubi, 2025). This has also been enabled by advancements in computing power, which now enables faster handling of complex molecular representations such as SMILES and molecular descriptors (Kosicka et al., 2022).

Numerous studies have been carried out to investigate the feasibility of using computational and statistical methods to predict the characteristics of polymers using their structural characteristics (Singh & Chandra, 2023). At first, linear regression and rule-based methods were employed, but lately, various ML approaches, including support vector machines, random forests, and neural networks, were adopted (Ding et al., 2023). All these methods have yielded encouraging outcomes for the efficient learning of non-linear associations between the molecular attributes and the target properties (Starkova et al., 2022). In addition, neural network architectures, such as convolutional neural networks and graph neural networks, have been employed to train representations using molecular data (Phua et al., 2023). All these advancements have led to improved prediction accuracy and more efficient design processes (Sofos et al., 2023). Traditional polymer property prediction has largely relied on group contribution methods and QSPR models such as those by (Krevelen, 2009) and (Bicerano, 1993), which offer interpretable but often limited performance for complex, non-linear systems. Recent advances in Materials Informatics (Parambil et al., 2025) introduce data-driven machine learning approaches that better capture structure–property relationships using large datasets. In this context, this study adopts an XGBoost-based framework with optimization to improve predictive accuracy and bridge classical and modern approaches.

In this regard, the present work aims to develop an ML-based framework for predicting polymer properties by utilizing structure-property relationship modeling. Data preprocessing and model building have been combined in the current study to enhance the property estimation process during the design and discovery of polymers.

Key Contribution

• A comprehensive ML framework was developed for predicting polymer properties using the PolyOne dataset, addressing nonlinear SPR challenges.
• Data preprocessing, feature extraction using molecular descriptors and topological indices, and feature selection using Recursive Feature Elimination were systematically performed.
• An XGBoost predictive model was implemented, and its hyperparameters were optimized using the Starfish Optimization Algorithm to ensure accurate and efficient predictions.
• Model interpretability was achieved using SHAP for global feature importance and LIME for local explanation, identifying key molecular descriptors influencing polymer properties.

The proposed framework is summarized as follows. The polymer dataset undergoes its first preprocessing stage through three processes of cleaning and normalization and data splitting. The extraction process produces numerical representations of polymer structures through the collection of molecular descriptors and topological indices. The feature selection process uses Recursive Feature Elimination to identify the most valuable attributes for analysis. The selected features are used to train an XGBoost model, and its hyperparameters are optimized using the Starfish Optimization Algorithm. The model uses SHAP and LIME techniques to show which features have the robust impact on its ability to predict polymer properties.

The remaining parts of this paper are organized as follows: Related Works and Research Gaps on Polymer Property Prediction are discussed in Section 2. Our method is explained in Section 3. In Section 4, details on the proposed XGBoost model, hyperparameter optimization using Starfish, and an explanation of the learned model using SHAP and LIME are provided. The experiment and its results are discussed in Section 5. Finally, conclusions are drawn in Section 6.

2. Literature Review

(Chong et al., 2024) explored recent advances in machine learning techniques for materials science applications. The study explained various data-driven approaches used for predicting material properties, optimizing material performance, and accelerating computational materials design. (Hatakeyama-Sato, 2023)presented a comprehensive review on the use of molecular descriptors and machine learning approaches in polymer design. The study discussed how data-driven techniques can effectively predict polymer properties and support advanced material development. (Zhao et al., 2023) reviewed the application of molecular descriptors and machine learning techniques in polymer design and highlighted their effectiveness in improving material property prediction accuracy.

(Fattouche et al., 2024) utilized methods such as artificial neural network-based QSAR modeling, molecular docking, ADMET, and molecular dynamics simulation for designing selective inhibitors. In this way, they showed the capabilities of combining SPRs with AI. Likewise, (Kazemi-Khasragh et al., 2024) applied ML along with group interaction models for predicting polymer properties. The authors illustrated the potential of incorporating machine learning into SPRs. According to (Kibrete et al., 2023), artificial intelligence methods were applied for predicting the mechanical behavior of composites. The authors showed the potential of supervised learning in SPRs. Similarly, Li et al. employed artificial intelligence methods to investigate the SPRs of 3D-printed lattices. The authors showed the capabilities of artificial intelligence within the framework of SPRs. In addition, (Z. Li et al., 2022) explored the capabilities of deep learning approaches in the realm of SPRs. All these studies demonstrated the power of integrating ML techniques with SPRs in different fields of materials science. However, one major limitation of these studies was the models’ ability to generalize across different polymer classes. Recent transformer-based polymer informatics frameworks such as polyBERT demonstrated the effectiveness of learned polymer embeddings for large-scale polymer property prediction (Kuenneth & Ramprasad, 2023). Contrasting transformer-driven latent representation learning, the proposed framework employs descriptor-based feature engineering integrated with XGBoost and SFOA optimization to improve interpretability and computational efficiency.

(Lv & Zhu, 2026) explored the inherent thermal conductivity of molecularly engineered polymers by stressing the significance of precise control over the structure of materials on thermal conductivity. (Milad et al., 2022) also explored ensemble ML strategies for predicting strain in fiber-reinforced polymer composites and demonstrated greater robustness than traditional approaches. The article by (Thomas et al., 2022) discusses how modeling was used to assess the environmental aging of materials and the lifespan of polymers, using predictive approaches to model material degradation. (Timmanaikar et al., 2024) investigated the potential of using graphical topology index values for predicting the SPRs of pharmacologically significant molecules, highlighting the importance of mathematical modeling in predictive science. In (Xu et al., 2022), the authors have examined the characterization of dendritic mesoporous nanoparticles, revealing the role of synthesis and structure in governing material properties. The above-mentioned research papers described various approaches and techniques for modeling and predicting the behavior of materials under diverse conditions. The key limitation of these studies was their inability to link macroscopic and microscopic material characteristics.

In this context, the rapid evolution of ML technology in polymer discovery was emphasized by (Yan & Li, 2023). Moreover, emphasized the efficiency of ML for predicting the performance of polymer composites. (Yue et al., 2022). They proposed a hybrid framework for simulating the performance of polymer composites for energy storage. Also, the study of the structural design of polymer semiconductors, especially hydrogen bonding, was covered by (Q. Zhang et al., 2022). In addition, the importance of applying ML technology to predict polymer performance was further highlighted by (X. Li et al., 2024). However, the main problem in applying ML to predict polymer performance is the lack of standardized databases.

2.1. Problem Statement

There are several critical challenges in predicting polymer properties using SPR modeling. The first challenge is that the SPR of polymers is complex and nonlinear (Xue et al., 2024). This makes it difficult to precisely predict the effects of changes in molecular configurations and functional groups (B. Zhang et al., 2023). The second challenge is that the computational tools used to predict the SPR of polymers, such as molecular simulations and theoretical calculations, are quite resource-intensive. The third challenge is that the availability of standardized polymer data sets remains quite low. The fourth challenge concerns selecting the most important molecular descriptors, a problem faced by most ML algorithms used to predict SPR (Basova et al., 2023).

Despite various solutions, such as statistical models, conventional machine learning techniques, and deep learning models, being suggested to overcome the previously mentioned obstacles, numerous difficulties remain in maintaining an ideal balance between accuracy, efficiency, and the interpretability of predictions. Specifically, the problem of effective and efficient feature engineering, ML modeling, and interpretable insights into SPRs remains unresolved. To address the aforementioned limitations and challenges, the paper proposes a framework encompassing effective and efficient data preprocessing, feature engineering and selection, ML modeling, and interpretable analysis.

3. Proposed Methodology

The proposed methodology outlines a comprehensive ML-based approach for predicting polymer properties using SPR modeling. This approach has been formulated with the problems of data-set inconsistency, the high dimensionality of molecular structure, and the non-interpretability of current approaches in mind. By integrating data preprocessing, feature engineering, modeling, optimization, and interpretability, the framework aims to enable an accurate, efficient, and interpretable prediction system for polymers. The importance of the proposed methodology lies in its ability to transform raw polymer structure information into useful insights effectively.

The proposed framework for polymer property prediction is illustrated in Figure 1. The workflow begins with the polymer dataset, which undergoes essential data preprocessing steps including cleaning, normalization, and handling of missing values to ensure data quality and consistency. Next, relevant features are extracted or selected using RFE or descriptor-based selection techniques, which reduce dimensionality and retain the most informative variables.

The processed dataset is then divided into training and test sets to enable unbiased model evaluation. The XGBoost model is subsequently employed for initial property prediction due to its strong performance in handling structured data. To further enhance model performance, the Starfish Optimization Algorithm (SOA) is applied for hyperparameter tuning, resulting in an optimized XGBoost model. Finally, the optimized model is used for polymer property prediction, and its performance is assessed using standard evaluation metrics such as MAE, RMSE, and R².

3.1. Feature Extraction of Polymer Structures

The polymer representation combines multiple categories of structural information, including physicochemical descriptors, topological indices, connectivity-based fingerprints, and substructure-derived molecular features extracted from PSMILES representations using RDKit. This multi-representation strategy enables the model to capture both global physicochemical characteristics and local structural patterns associated with polymer behavior. In addition, recursive feature elimination (RFE) is employed to identify the most informative features and descriptors before model training.

The selection of molecular descriptors and topological indices in this method was based on their ability to capture key physicochemical and structural characteristics which determine polymer material properties. Descriptors which included molecular weight and LogP and hydrogen bonding features and graph-based indices (e.g., Wiener and Randic indices) were selected to show molecular size and polarity and intermolecular interactions and structural connectivity because these features created an extensive yet informative set of characteristics. The study selected a balanced subset of descriptors which would enable the model to maintain its performance during interpretation tasks and its ability to generalize across different situations.

Polymer representation in the proposed framework is derived from molecular and topological descriptor extraction rather than transformer-based embedding generation. Descriptor-based representation was adopted to enable interpretable structural characterization, lower computational overhead, and effective integration with the XGBoost prediction framework for polymer property estimation.

3.1.1. Molecular Descriptors

Physicochemical descriptors quantify global chemical and structural characteristics of polymers. The study created a descriptor space which includes multiple molecular weight and polar surface area and LogP and hydrogen-bonding and rotatable bond and aromaticity and electronic property descriptors through RDKit. The descriptors define molecular size and molecular polarity and molecular flexibility and molecular intermolecular interactions and molecular structural complexity which scientists use to predict polymer properties. The Molecular Weight (MW) is given in Eq. (1)

$$MW={\displaystyle {\sum }_{i=1}^N} s_i$$

(1)

where $S_i$ represents the polar surface contribution of the $i^{\mathrm{th}}$ atom and $N$ denotes the total number of polar atoms in the molecular structure.

Standard molecular and topological descriptors employed for polymer representation, including Polar Surface Area (PSA), molecular weight, hydrogen bond donor and acceptor counts, rotatable bond characteristics, and topology-based structural parameters, were generated using established descriptor computation techniques commonly adopted in recent cheminformatics and polymer informatics frameworks to ensure reliable structural characterization and feature consistency for polymer property prediction tasks (Hu et al., 2024).

3.1.2. Topological Indices

Topological and connectivity-based descriptors were additionally employed to characterize the arrangement and interaction patterns of atoms within polymer chains. The descriptors include Wiener index and Randic connectivity index and eccentricity measures and extended-connectivity fingerprint ECFP-derived structural representations. The features capture branching patterns and cyclic structures and atomic neighborhood information and molecular connectivity which are linked to polymer thermomechanical behavior. The Wiener index is given in Eq. (2)

$$W={\displaystyle \frac{1}{2}}{\displaystyle {\sum }_{i=1}^N} {\displaystyle {\sum }_{j=1}^N} d_{ij}$$

(2)

In this expression, $d_{ij}$ is the shortest path (number of bonds) between atoms $i$ and $j$. The Randic connectivity is given in Eq. (3)

$$\chi ={\sum }_{(i,j)\in E} {\displaystyle \frac{1}{\sqrt{d_i{d}_j}}}$$

(3)

The set bond is denoted as $E$, $d_i$ and $d_j$ are the degrees (number of bonded neighbors) of atoms $i$ and $j$. A numeric vector for each polymer will be created from both physicochemical features and topological indices, serving as input to the ML algorithm.

Figure 2 depicts the feature extraction process, in which polymer SMILES are converted into numerical feature vectors, which are then used in the ML process. In the feature extraction process, molecular descriptors, including physicochemical parameters such as molecular weight, LogP, hydrogen bond acceptors (HBA). It represents the number of atoms capable of accepting hydrogen bonds, and hydrogen bond donors (HBD), which represent the number of atoms capable of donating hydrogen bonds, are determined from the polymer structures and used as input features for the machine learning model.

3.2. Feature Selection Using Recursive Feature Elimination

Identify the most informative molecular descriptors and topological indices, and eliminate redundant or irrelevant ones to improve performance and reduce complexity.

3.2.1. Recursive Feature Elimination

The Recursive Feature Elimination (RFE) technique is an iterative method in which a model is trained on the dataset and features are ranked based on their importance. The least significant features are progressively removed, and the process is repeated until the desired number of features is obtained. This approach helps eliminate redundant features, improves model accuracy, reduces overfitting, and enhances computational efficiency through dimensionality reduction. A predictive model is trained on the feature set $X=\left[x_1,x_2,\dots ,x_m\right]$ and the target property $y$ is given in Eq. (4)

$$y=f\left(X\right)+ϵ$$

(4)

The vector of target polymer properties is denoted as $y,X$ is the feature matrix of size $n\times m$ (n samples, $m$ features), $f\left(X\right)$ is the predictive function learned by the model and $ϵ$ is the residual error.

Algorithm 1 iteratively trains a model using the current features, calculates the importance of each feature, and removes the least important feature until only the required features remain. The process will eliminate redundant features, making the feature set more informative and reducing the chances of overfitting.

Algorithm 1: Feature Selection using RFE

Input: Feature matrix $\mathit{x}$, target vector y, desired number of features $\mathit{k}$ Output: Selected feature subset $\mathit{X}\_\mathit{s}\mathit{e}\mathit{l}\mathit{e}\mathit{c}\mathit{t}\mathit{e}\mathit{d}$ 1. Initialize feature set $\mathbf{F}=$ all features in $\mathit{X}$ 2. While the number of features in $\mathbf{F}>\mathbf{k}$:           a. Train model M on $\mathit{X}[:,\mathit{F}]$ to predict $\mathit{y}$           b. Compute feature importance scores ${\mathit{T}}_{\mathit{j}}$ for each feature in $\mathbf{F}$           c. Identify feature $\mathbf{F}\_\mathbf{m}\mathbf{i}\mathbf{n}$ with the lowest importance           d. Remove $\mathbf{F}\_\mathbf{m}\mathbf{i}\mathbf{n}$ from $\mathbf{F}$ 3. Return $\mathit{X}\_\mathit{s}\mathit{e}\mathit{l}\mathit{e}\mathit{c}\mathit{t}\mathit{e}\mathit{d}$= $\mathit{X}[:,\mathit{F}]$

3.3. Machine Learning-Based Predictive Modeling Using XGBoost

The molecular descriptors and topological features were used to train the XGBoost regression model for polymer property prediction. XGBoost was selected due to its strong capability in modeling nonlinear structure-property relationships and handling high-dimensional feature spaces efficiently. The prediction model is expressed in Eq. (5).

$${\hat{y}}_i={\displaystyle {\sum }_{k=1}^K} f_k\left(x_i\right)$$

(5)

where ${\hat{y}}_i$ represents the predicted polymer property, $K$ denotes the number of decision trees, and $f_k\left(x_i\right)$ represents the prediction generated by the $k$-th tree for sample $i$.Eq. (6)

$${\hat{y}}_i^{\left(t\right)}={\hat{y}}_i^{(t-1)}+\eta f_t\left(x_i\right)$$

(6)

The learning rate is denoted as $\eta$, ${\hat{y}}_i^{\left(t\right)}$is the to predict after $t$-th tree and $f_t\left(x_i\right)$ is the output of the current tree for the sample $i$. This process enables iterative refinement of predictions by incorporating corrections from each successive tree.

In Figure 3, the selected feature vectors are used to predict the polymer’s properties using the XGBoost model. Multiple gradient-boosted tree models are used to obtain the final prediction. These predictions are determined based on the sum of weights for the prediction made by each tree.

3.4. Hyperparameter Optimization Using the SOA

This process aims to optimize the hyperparameters of the XGBoost model to yield optimal prediction performance for the polymer properties. The Starfish Optimization Algorithm (SOFA) algorithm is applied during optimization. This optimization algorithm is bio-inspired by the foraging behavior of starfish. The algorithm optimizes the model’s hyperparameters. The hyperparameter tuning process for the XGBoost model focused on optimizing tree-based parameters rather than neural network-specific settings. The study examined five essential hyperparameters which included number of estimators learning rate maximum tree depth subsample ratio and column sampling rate. The parameters establish boundaries for the model’s learning process and capacity to generalize from its training data. The optimal configuration was determined through systematic tuning methods. The procedure achieves better prediction results because it prevents the model from overfitting to the training data.

The SFOA was designated to maintains an equal distribution of search abilities which enables successful movement through complex hyperparameter search environments. SFOA offers search methods which change according to the search process to prevent the system from getting trapped in local optimum solutions which traditional methods such as Grid Search and Random Search cannot achieve. SFOA uses a foraging-based system which provides better stability during convergence and decreases processing requirements when compared to common metaheuristic methods which include Genetic Algorithms (GA) and Particle Swarm Optimization (PSO). The system displays characteristics which make it suitable for XGBoost hyperparameter optimization in search environments that have both high dimensions and non-linear characteristics.

The classical metaheuristic techniques such as Genetic Algorithms (GA) and Particle Swarm Optimization (PSO), SFOA provides improved convergence stability and reduced computational overhead due to its foraging-inspired adaptive search strategy. Furthermore, when compared to Bayesian Optimization, which relies on surrogate probabilistic modeling and sequential acquisition functions, SFOA does not require explicit probabilistic modeling of the objective function, making it more suitable for high-dimensional and highly non-convex hyperparameter spaces where surrogate modeling may become computationally expensive or less stable. In this study, a comparative evaluation was conducted between SFOA and standard optimization baselines, and the results demonstrated that SFOA achieved superior or comparable performance in terms of validation error and convergence stability for XGBoost hyperparameter tuning. These characteristics justify the use of SFOA for optimizing complex tree-based models under nonlinear search conditions.

The optimization objective was to maximize prediction performance while minimizing model overfitting. The optimized hyperparameter set is represented as Eq. (7):

$${\theta }^{*}=\left[{\eta }^{*},n^{*},d^{*},s^{*},c^{*}\right]$$

(7)

where ${\eta }^{*}$ denotes the optimized learning rate, $n^{*}$ represents the number of estimators, $d^{*}$ indicates maximum tree depth, $s^{*}$ corresponds to the subsample ratio, and $c^{*}$ denotes the column sampling ratio.

The fitness function used during optimization is expressed in Eq. (8).

$$F\left(\theta \right)=\mathrm{P}\mathrm{e}\mathrm{r}\mathrm{f}\mathrm{o}\mathrm{r}\mathrm{m}\mathrm{a}\mathrm{n}\mathrm{c}\mathrm{e}(\hat{y},y)$$

(8)

where $F\left(\theta \right)$ represents the fitness score of the hyperparameter set, $\hat{y}$ denotes predicted polymer properties, and $y$ represents the actual property values.

3.5. Structure-Property Relationship Analysis Using SHAP and LIME

Structure Property Analysis uses SHAP for global interpretation and LIME for local interpretation to explain how molecular features influence polymer properties. SHAP provides an overall picture of the contribution of descriptors to polymers’ behavior, while LIME explains the prediction for individual instances by approximating the model locally. This step identifies the key molecular features that govern the structure-property relationship, providing insight into which structural elements are most critical for designing polymers with desired properties. The SHAP value is given in Eq. (9)

$${\varphi }_j={\sum }_{S\subseteq F\setminus \left\{j\right\}} {\displaystyle \frac{\left|S\right|!\left(\right|F|-|S|-1)!}{\left|F\right|!}}\left[f_{S\cup \left\{j\right\}}\left(x_{S\cup \left\{j\right\}}\right)-f_S\left(x_S\right)\right]$$

(9)

The set of all features is denoted as $F$, ${\varphi }_j$ is the role of the feature $j$, $S$ is the subgroup of features excluding $j$ and $f_S\left(x_S\right)$ is the to predict using the subset $S.$ Quantifies the change caused by the feature in the prediction across all combinations and provides a fair attribution to the feature. The SHAP prediction decomposition is given in Eq. (10)

$$\hat{y}={\varphi }_0+{\displaystyle {\sum }_{j=1}^M} {\varphi }_j$$

(10)

The base value is denoted as ${\varphi }_0$, $\hat{y}$ is the predicted polymer property, ${\varphi }_j$ is the contribution of the feature $j$. It breaks down the prediction into additive contributions, making the prediction interpretable, as given in Eq. (11)

$$g\left(z^{\prime }\right)=w_0+{\displaystyle {\sum }_{j=1}^M} w_j{z}_j^{\prime }$$

(11)

Where, $g\left(z^{\prime }\right)$is the local surrogate model, $z_j^{\prime }$ is the binary representation of feature j for the perturbed input and $w_j$ is the weight representing feature importance locally. It ensures the local explanation is accurate but simple. LIME weight is given in Eq. (12)

$$\mathrm{a}\mathrm{r}\mathrm{g} \underset{g}{\mathrm{m}\mathrm{i}\mathrm{n}} L\left(f,g,{\pi }_x\right)+\mathsf{\Omega }\left(g\right)$$

(12)

In this equation, $L\left(f,g,{\pi }_x\right)\mathrm{i}\mathrm{s}\mathrm{t}\mathrm{h}\mathrm{e}$proximity ${\pi }_x$ loss among the local classical $g$ and the unique model $f$, weighted by proximity ${\pi }_x$ and $\mathsf{\Omega }\left(g\right)$ is the difficulty consequence for interpretability. The ranks feature globally to identify the most critical molecular descriptors is given in Eq. (13)

$$I_j={\displaystyle \frac{1}{N}}{\displaystyle {\sum }_{i=1}^N} \left|{\varphi }_{ij}\right|$$

(13)

The overall importance of the feature is denoted as $I_j,{\varphi }_{ij} is the$ SHAP value of the feature $j$ for model $i$ and $N is the$ overall number of samples.

As shown in Figure 4, it entails analyzing the predicted properties and characteristics of the polymers using SHAP and LIME methods. SHAP is a tool that determines the characteristics of polymers that most influence the prediction process, while LIME explains the prediction process for polymers.

4. Experimental Setup

4.1. Dataset Description

The polymer dataset employed in the study is drawn from the PolyOne Data (“PolyOne Data Set - 100 Million Hypothetical Polymers Including 29 Properties” 2022), which contains 100 million polymers characterized by 29 computationally predicted properties. The polymer data set represents PSMILES notation for polymers. Each PSMILES string describes a single polymer generated by randomly combining unique polymer substructures derived from the decomposition of existing synthetic polymers. In general, only a small fraction of the generated polymers has been experimentally synthesized. Consequently, the data set can be considered highly suitable for the task of predicting properties of synthetic polymers through computational approaches.

To create polymer representation models, such as polyBERT, we divided the dataset into two groups: a training dataset containing 80 million PSMILES strings and a test dataset containing 20 million PSMILES strings. The computational cost associated with feature extraction, recursive feature elimination, and explainable AI analysis, a representative subset of 5,000 polymer samples was selected from the complete PolyOne dataset while preserving structural diversity and property distribution characteristics.

The proposed framework utilizes descriptor-based polymer representation generated from molecular and topological descriptor extraction techniques instead of transformer-based embedding architectures such as PolyBERT. This representation strategy supports interpretable structural characterization, reduced computational complexity, and efficient integration with the XGBoost prediction framework for polymer property estimation.

The polymer database contains about 100 million polymers yet researchers chose 5000 samples to test because of their ability to process data and their need for high-quality results. The researchers used random stratified sampling to create a subset that maintained the original dataset’s statistical distribution and various structural and property characteristics. The preprocessing phase removed duplicate and incomplete and inconsistent records to create reliable data. The selected subset therefore provides a representative distribution of the larger polymer dataset while maintaining manageable computational complexity for model training and optimization.

Table 1 compares the mean values of important physicochemical properties between the complete polymer dataset and the selected 5000-sample subset. The minimal percentage differences observed across all properties confirm that the selected subset effectively preserves the statistical characteristics and distribution patterns of the original dataset, thereby ensuring representative and reliable model training.

4.2. Data Preprocessing

Data preprocessing is a significant component of the suggested system because the raw data from the data repository, such as SMILES strings and property values, may contain scale errors. The data is converted to a clean, normalized, and noise-free form.

4.2.1. Data Cleaning

Data cleaning is important during the preprocessing stage to eliminate missing or inconsistent data and prevent an impact on data quality and integrity. This way, in polymer datasets, missing data may introduce some bias during analysis, leading to inefficiency. The focus, therefore, is to ensure that any records containing any form of null or undefined values are removed, ensuring that only meaningful and valid information is retained, making the entire process consistent and representative of the structure and properties given in Eq. (14)

$$D={\left\{\left(x_i,y_i\right)\right\}}_{i=1}^N$$

(14)

Where, $D$ is the target property value, $x_i$ is the molecular descriptor vector derived from SMILES, $y_i$ is the target polymer property and $N$is the descriptor features.

4.2.2. Normalization Using Min-Max Scaling

The aim of the Min–Max normalization is to scale the feature values to the range [0,1]. This is important because the feature values are of different orders of magnitude. This ensures that the features with higher values do not overshadow the others. This improves the numerical stability of the model and facilitates smooth, fast convergence during training. The normalization is given in Eq. (15)

$$x_{\mathrm{norm}}={\displaystyle \frac{x-x_{\mathrm{m}\mathrm{i}\mathrm{n}}}{x_{\mathrm{m}\mathrm{a}\mathrm{x}}-x_{\mathrm{m}\mathrm{i}\mathrm{n}}}}$$

(15)

The value $x$ signifies the original rate of the feature, and $x_{\mathrm{m}\mathrm{i}\mathrm{n}}$ and $x_{\mathrm{m}\mathrm{a}\mathrm{x}}$ represent the least and extreme values of a particular feature, respectively. The two values represent the range of a particular feature and are used to scale the values proportionally. The output value, $x_{\mathrm{norm}}$, represents a normalized value, and its range is 0 to 1 after transformation.

In addition, the dataset is divided into training, validation, and test sets at a 70:15:15 ratio. This process plays a significant role in improving model training, alongside hyperparameter optimization and performance assessment. In short, this preprocessing technique improves the quality of the data.

Algorithm 3 for preprocessing the raw polymer dataset is as follows. First, we eliminate any null values in the dataset. Then, all the attributes are scaled using the Min-Max scaling technique. Lastly, the data is partitioned into training, validation, and test datasets. This helps us train our ML model effectively. Detailed hardware and software specifications are provided in the Supporting Information.

Algorithm 3: Data Preprocessing for Polymer Property Prediction

Input: Raw dataset $\mathit{D}\_\mathit{r}\mathit{a}\mathit{w}$ with features $\mathit{X}$ and target $\mathit{Y}$ Output: Cleaned and normalized dataset $\mathit{D}\_\mathit{c}\mathit{l}\mathit{e}\mathit{a}\mathit{n}$ 1. Load dataset $\mathit{D}\_\mathit{r}\mathit{a}\mathit{w}$ 2. Handle missing values:           a. For each feature in $\mathit{X}$:                     i. If missing values exist, replace with mean/median or remove row 3. Normalize features using Min-Max scaling:           a. For each feature $\mathit{x}$ in $\mathit{X}$:                     i. x$\_\mathit{s}\mathit{c}\mathit{a}\mathit{l}\mathit{e}\mathit{d}=$ = $(\mathit{x}-\mathbf{m}\mathbf{i}\mathbf{n}(\mathit{x}\left)\right)/(\mathbf{max}\left(\mathit{x}\right)-\mathbf{m}\mathbf{i}\mathbf{n}(\mathit{x}\left)\right)$ 4. Split the dataset into:           a. Training set (70%)           b. Validation set (15%)           c. Test set (15%) 5. Return $\mathit{D}\_\mathit{c}\mathit{l}\mathit{e}\mathit{a}\mathit{n}=\{\mathit{X}\_\mathit{t}\mathit{r}\mathit{a}\mathit{i}\mathit{n},\mathit{X}\_\mathit{v}\mathit{a}\mathit{l},\mathit{X}\_\mathit{t}\mathit{e}\mathit{s}\mathit{t},\mathit{Y}\_\mathit{t}\mathit{r}\mathit{a}\mathit{i}\mathit{n},\mathit{Y}\_\mathit{v}\mathit{a}\mathit{l},\mathit{Y}\_\mathit{t}\mathit{e}\mathit{s}\mathit{t}\}$

4.3. Performance Metrics

The efficiency of the proposed machine learning approach for estimating the properties of polymer materials is assessed using traditional evaluation metrics and validation methods.

R² (Coefficient of Determination)

The R² score measures the goodness-of-fit of the regression model’s predictions against the actual values by quantifying the proportion of the dependent variable’s total variance accounted for by the model. The closer the value to 1, the more effective the model. R² can be calculated using Eq. (16).

$$R^2=1-{\displaystyle \frac{{\displaystyle {\sum }_{i=1}^n} {\left(y_i-{\hat{y}}_i\right)}^2}{{\displaystyle {\sum }_{i=1}^n} {\left(y_i-\bar{y}\right)}^2}}$$

(16)

Here, $y_i$ is the real value, ${\hat{y}}_i$ is the predicted value, $\bar{y}$ is the mean of real ideals, and $n$ is the number of samples.

MAE

MAE measures the average size of the difference between predicted and observed values without considering the direction of the error. It is easily interpretable as a measure of predict accuracy. The MAE is given in Eq. (17)

$$MAE={\displaystyle \frac{1}{n}}{\displaystyle {\sum }_{i=1}^n} \left|y_i-{\hat{y}}_i\right|$$

(17)

Where, $y_i$ is the real value, ${\hat{y}}_i$ is the predicted value, $\bar{y}$ is the mean of real values, and $n$ is the number of samples.

Root Mean Square Error

The formula for RMSE is the square root of the average of the squared differences between observed and predicted values. RMSE is highly sensitive to large differences, which means it punishes them severely. RMSE is presented in Eq. (18).

$$RMSE=\sqrt{{\displaystyle \frac{1}{n}}{\displaystyle {\sum }_{i=1}^n} {\left(y_i-{\hat{y}}_i\right)}^2}$$

(18)

Mean Squared Error

MSE represents the mean of the formed deviations of the predicted values $\left({\hat{y}}_i\right)$from the true ones $\left(y_i\right)$. The use of the square of deviations makes the MSE metric highly sensitive to the large deviations, as a result of which the model error is strongly influenced by outliers. MSE can be used to predict polymer properties from molecular descriptors using the XGBoost ML algorithm Eq. (19).

$$\mathrm{M}\mathrm{S}\mathrm{E}={\displaystyle \frac{1}{n}}{\displaystyle {\sum }_{i=1}^n} {\left(y_i-{\hat{y}}_i\right)}^2$$

(19)

Where, $n$ is the total amount of polymer trials, $y_i$ is the actual value of the property for the $i$-th polymer and ${\hat{y}}_i$ is the predicted value of the property for the $i$-th polymer.

5. Results and Discussion

In the generalization capability of the proposed framework, experiments were conducted on multiple representative polymer properties obtained from the PolyOne dataset, including glass transition temperature (Tg), melting temperature (Tm), and decomposition temperature (Td). The reported results primarily highlight Tg prediction due to its importance in polymer design, while similar predictive trends were observed for the remaining properties. The experiment indicates that polymer properties can be predicted using structure-property relationship modeling. The objective of this research is to develop an efficient machine learning model for predicting the polymer properties and understanding structure-property relationships. The experiments have been performed using the PolyOne Dataset, consisting of 100 million imaginary polymers coded as PSMILES with 29 predicted properties. The sample size was 5,000, divided so that the training set comprised 70%, the validation set 15%, and the test set 15%. 5-fold cross-validation with a seed value of 42 was used to ensure reproducibility of results. The model was created using Python, XGBoost, scikit-learn, RDKit, SHAP, LIME, NumPy, and Pandas. The performance of the XGBoost model optimized with SOA is commendable, with an R² score of 0.92, an MAE of 0.08, and an RMSE of 0.12. The SHAP and LIME models have identified the critical molecular properties that influence the polymer’s behavior. Detailed technical information is provided in the Supporting Information (Table 2 and Table 3).

Supporting information

Hardware and software settings for our investigation are listed in Table 1 below.

Table 2. Hardware and Software Configuration for Research Analysis.

Category	Component	Specification / Version
Software	scikit-learn	1.4.0
	SOA	custom
	RDKit	2023.9
	SHAP	0.44.0
	LIME	0.2.0.1
	Pandas	2.1.4
	NumPy	1.26.3
	Pyarrow	14.0.1
	Scipy	1.11.4
Hardware	CPU Model	Intel Core i7-12700H
	CPU Max Boost	4.70 GHz
	RAM Capacity	16 GB DDR5
	GPU Model	NVIDIA GeForce RTX 3060 (Laptop)
	Storage Type	NVMe SSD
	Storage Capacity	512 GB
	Storage Read Speed	3500 MB/s

Table 2. Hardware and Software Configuration for Research Analysis.

Category	Component	Specification / Version
Software	scikit-learn	1.4.0
	SOA	custom
	RDKit	2023.9
	SHAP	0.44.0
	LIME	0.2.0.1
	Pandas	2.1.4
	NumPy	1.26.3
	Pyarrow	14.0.1
	Scipy	1.11.4
Hardware	CPU Model	Intel Core i7-12700H
	CPU Max Boost	4.70 GHz
	RAM Capacity	16 GB DDR5
	GPU Model	NVIDIA GeForce RTX 3060 (Laptop)
	Storage Type	NVMe SSD
	Storage Capacity	512 GB
	Storage Read Speed	3500 MB/s

The hyperparameter Table 3 used for the optimization of $T_g$ prediction using the XGBoost algorithm was fine-tuned for optimal accuracy and stability of the model. These comprise the total number of boosting rounds used at 2450, maximum depth at 8, learning rate at 0.0312, and row and column subsampling at 0.82 and 0.75-0.72, respectively. Regularisation is implemented using Min_child_weight is 3, and Max_delta_step is 2. Base_score = 0.5 is used as an initial prediction value. Feature selection and tuning were done using RFE for 150 features, 50 iterations, and 5-fold cross-validation. A Min-Max scaler was employed for input normalization.

Table 2. Hyperparameter Configuration of the proposed XGBoost polymer prediction properties.

Model	Hyperparameter	Value
N_estimators	Amount of Boosting Trees	2450
Max_depth	Extreme Tree Depth	8
Learning_rate	Learning Rate η	0.0312
Subsample	Row Subsampling Ratio	0.82
Colsample_bytree	Column Fraction per Tree	0.75
Colsample_bylevel	Column Fraction per Level	0.68
Colsample_bynode	Column Fraction per Node	0.72
Min_child_weight	Minimum Sum of Hessian in Leaf	3
Max_delta_step	Max Delta Step per Tree	2
Base_score	Initial Prediction Score	0.5
N_jobs	Parallel CPU Threads	-1
N_features	Features Selected via RFE	150
N_trials	Optimisation Trials	50
Cv_folds	Cross-Validation Folds	5
Scaler	Feature & Target Scaler	Min Max Scaler [0,1]
Rfecv_step	RFE Step Size	0.05
Rfecv_cv	RFE Cross-Validation Folds	3
Rfecv_scoring	RFE Scoring Metric	R2

Table 2. Hyperparameter Configuration of the proposed XGBoost polymer prediction properties.

Model	Hyperparameter	Value
N_estimators	Amount of Boosting Trees	2450
Max_depth	Extreme Tree Depth	8
Learning_rate	Learning Rate η	0.0312
Subsample	Row Subsampling Ratio	0.82
Colsample_bytree	Column Fraction per Tree	0.75
Colsample_bylevel	Column Fraction per Level	0.68
Colsample_bynode	Column Fraction per Node	0.72
Min_child_weight	Minimum Sum of Hessian in Leaf	3
Max_delta_step	Max Delta Step per Tree	2
Base_score	Initial Prediction Score	0.5
N_jobs	Parallel CPU Threads	-1
N_features	Features Selected via RFE	150
N_trials	Optimisation Trials	50
Cv_folds	Cross-Validation Folds	5
Scaler	Feature & Target Scaler	Min Max Scaler [0,1]
Rfecv_step	RFE Step Size	0.05
Rfecv_cv	RFE Cross-Validation Folds	3
Rfecv_scoring	RFE Scoring Metric	R2

Table 4 presents the prediction performance of the proposed framework for multiple representative polymer properties obtained from the PolyOne dataset. The evaluation metrics demonstrate the generalization capability and prediction consistency of the proposed model across different polymer property categories.

The proposed framework maintains its stochastic processes which include dataset splitting and Recursive Feature Elimination (RFE) and Starfish Optimization Algorithm (SFOA)-based hyperparameter tuning that were tested using different random seed configurations. The researchers repeated the experiments by using five different random seeds which included the seeds 10 and 20 and 30 and 40 and 50 to measure average performance metrics. The prediction accuracy and error metrics showed only minor changes which proved the framework maintained its stability across different starting conditions.

Table 5 presents the stability analysis of the proposed framework under different random seed initializations. The obtained RMSE, MAE, and R² values exhibit only minor variations across multiple experimental runs, indicating that the proposed model maintains consistent predictive performance irrespective of stochastic initialization conditions.

5.1. Correlation Among Actual and Predicted $T_g$ Values Performance Metrics

The plot shows the relationship between the actual and predicted $T_g$ values, helping evaluate the ML algorithm’s predictive accuracy. By comparing the actual $T_g$ value with the estimated $T_g$ value, we can assess the effectiveness of our algorithm. The dashed red line shows perfect prediction when the predicted $T_g$ equals its true value.

Figure 5 shows that most $T_g$ -predicted values cluster near the 1:1 line, indicating accurate predictions with minimal deviation from the line. Given the high R² of 0.9886, it can be inferred that the model accounts for approximately 99% of the variance in $T_g$. In effect, this means that the model well represents the relationship between structure and properties of the polymers.

5.2. Model Convergence Analysis Using RMSE Learning Curve

It evaluates the convergence and training efficiency of the XGBoost model for polymer property prediction using the RMSE learning curve. By plotting the root mean squared error versus boosting rounds, we can visualize the learning efficiency. The graph is a means of diagnosing underfitting or overfitting of the model and the correctness of the choice of hyperparameters.

The value of RMSE drops sharply within a few iterations, indicating that the model is quick to learn and captures all essential elements within the polymer data set. After 400 iterations, the RMSE converges to 0.048, indicating that the model has fully trained and that further learning does not improve predictions. The RMSE learning curve (Figure 6) shows how the model performs with new data. The initial boosting rounds show a quick RMSE drop which proves that the system learns effectively while maintaining its required performance. The RMSE shows stable results as the iterations progress which indicates that overfitting does not occur because the results remain consistent without major changes. The chosen hyperparameters establish an equal balance between bias and variance which results in a model that maintains proper generalization performance.

5.3. Model Convergence Analysis Using MAE Learning Curve

The learning curve evaluates the convergence and learning efficiency of the XGBoost model for polymer property prediction using the MAE metric. The figure displays the relationship between the number of boosting rounds and MAE. The figure will help establish how effectively the model is training to make predictions from the available dataset.

From Figure 7, one can observe a rapid decrease in MAE during the first few boosting stages, indicating that the ML algorithm quickly learns the general pattern in the dataset. After about 400 boosting stages, MAE plateaus at a value of 0.037, implying that any other boosting stage beyond this does not offer more value.

5.4. Overall Model Performance Evaluation

The graph for performance evaluation is designed to provide a thorough analysis of the predictive power of the XGBoost machine learning algorithm using MAE, RMSE, and MSE. The purpose of such visualization is to demonstrate the algorithm’s predictive capabilities for polymer physical properties.

Figure 8 presents the comparative RMSE performance of the evaluated models. Lower RMSE values indicate better prediction accuracy and smaller deviations between the predicted and actual values. The proposed model achieves the lowest RMSE, demonstrating improved predictive performance compared to the existing approaches.

The development of error metrics throughout the boosting rounds appears in Figure 6 and Figure 7. The model shows effective learning through its first iterations which resulted in a rapid decline of both RMSE and MAE metrics. The performance enhancement rate decreases after a specific number of boosting rounds because extra iterations only produce minor performance improvements.

5.5. Residual Plot of Predicted and Actual Glass

In this section, the predictive ability of the developed model is discussed with respect to its accuracy in predicting $T_g$ The residuals graph plots the discrepancies between the actual and predicted $T_g$ values, helping us evaluate the robustness of the ML model within the structure-property relationship framework.

The residuals plot of the predicted glass transition temperature is illustrated in Figure 9. The graph is the residuals plot of the ML model used to predict the glass transition temperature ($T_g$). The x-axis represents the predictor for temperature values in the range 330 K to 600 K, whereas the y-axis indicates the residual error between -80 K and +50 K. The broken line at zero temperature shows the condition under which the discrepancy between the prediction and the actual value vanishes. In this case, most points are scattered randomly near the zero line, suggesting that the model predictions are unbiased.

5.6. Feature Importance Analysis Using SHAP Beeswarm Plot

The SHAP Beeswarm graph shows the contribution of each variable to the predicted polymer properties. Here, a Beeswarm plot demonstrates feature importance and the effect on the model’s output. The plot shows that the color scale is used to represent feature values, from blue (lowest values) to red (highest values). Such representation helps evaluate which feature influences model output more effectively.

From Figure 10, one can see that features such as Tm, Td, and Rot Bonds contribute the highest SHAP values, indicating they play an essential role in model prediction. Features with positive SHAP values raise the prediction, while negative ones lower it. The placement of SHAP values along the horizontal axis indicates the strength of the influence of the corresponding feature on the model’s results.

5.7. Local Feature Contribution Analysis Using LIME

The LIME plot provides a local interpretation of the ML model’s predictions, highlighting the contribution of individual features to polymer property outcomes. It sheds light on the role of individual features in predicting a particular property of the polymer. The x-axis shows the weight or the influence of each feature, and the y-axis shows the features used. Positive weights (represented in green bars) increase the value of the predicted property, whereas negative weights (shown in red bars) decrease it.

From Figure 11, one can observe that features such as Tm, Rot Bonds, and Td have the most prominent positive effect on the value of the predicted property, with Tm having the most significant contribution. ECFP_115, ECFP_23, and ECFP_142 are some other characteristics that have a negative impact on the predicted property value, thus reducing the predicted property value by a small amount. The SHAP and LIME analysis revealed Tm (melting temperature) and Td (decomposition temperature) and Rot Bonds (number of rotatable bonds) along with specific fingerprint descriptors ECFP_115 and ECFP_23 and ECFP_142 as the main features that influence the results. The ECFP-based features correspond to specific substructural patterns encoded through Extended-Connectivity Fingerprints which represent the presence of distinct molecular fragments that determine polymer characteristics.

The feature importance analysis (Figure 11 and Figure 12) indicates that melting temperature (Tm) and decomposition temperature (Td) are dominant predictors for glass transition temperature ($T_g$). This observation is consistent with the underlying thermophysical behavior of polymers, as these properties are intrinsically related to molecular structure, chain mobility, and thermal stability. From a practical perspective, estimating $T_g$ from Tm and Td is meaningful because $T_g$ is often more difficult and time-consuming to measure experimentally, whereas Tm and Td are more readily available in material databases and literature. Therefore, the proposed model can serve as a useful tool for rapid screening and prediction of $T_g$ when only limited thermal data is accessible.

5.8. Residual Distribution Analysis

This plot shows how the residuals from ML model predictions are distributed. The residuals represent the disparities between the predicted and actual properties of the polymers. The X-axis represents the value of the residuals, while the Y-axis represents the frequency of occurrences. The red dashed line represents the value of 0 for perfect predicted, where the predicted and actual values are equal.

As shown by Figure 12 above, most residuals are near 0, indicating a normal distribution without visible skew. Thus, the predictions of ML models can be considered unbiased because errors are distributed uniformly around 0. The low residual spread indicates the precision of the predictions.

5.9. Comparison with Existing Methods

In comparison with existing techniques, the outcomes from several models will be used to evaluate the impact of different algorithms on MSE, RMSE, MAE, and R². The developed model outperforms all other models across all evaluation metrics.

In Table 6, the MAE value of 0.0420 is slightly higher than that of ML + Group Interaction Modeling (MAE 0.041), but remains competitive while achieving superior performance in terms of R² and overall prediction. The proposed XGBoost + SOA model was compared with three existing machine learning approaches reported in previous studies, including ML with Group Interaction Modeling, Super Learning for Composite Materials, and ML for Flory–Huggins Parameter Prediction. The ML with Group Interaction Modeling approach estimates material behavior using molecular interaction features. The Super Learning model combines multiple learners to improve composite material prediction accuracy. The ML model for Flory–Huggin’s parameter prediction focuses on thermodynamic parameter estimation using regression-based learning. The proposed XGBoost + SOA framework differs from these approaches by integrating gradient boosting with the Seagull Optimization Algorithm (SOA) to enhance parameter optimization and prediction performance.

Table 6 involves methods related to non-identical polymer prediction tasks and datasets. Therefore, the results are intended to provide a general performance perspective rather than a strictly controlled benchmark comparison. Despite these differences, the proposed model demonstrates competitive and consistent performance across multiple evaluation metrics.

Table 7 presents the computational cost analysis of the proposed framework. Among all stages, the SFOA-based optimization process required the highest execution time due to its iterative search mechanism, while the feature selection stage effectively reduced input dimensionality and minimized overall training overhead. The obtained results demonstrate that the proposed framework remains computationally feasible and scalable for practical polymer informatics applications.

Table 8 compares the predictive performance of the proposed XGBoost–SFOA framework with the polyBERT model for polymer property prediction. The proposed framework achieved higher R² values and lower MAE and RMSE values, indicating improved prediction accuracy and error minimization compared with the transformer-based polyBERT approach.

5.10. Discussion

The results obtained from this research indicate the success of the proposed machine learning framework for accurate predictions of the polymers’ characteristics based on understanding their complicated SPRs. This involves employing techniques such as data preprocessing and feature extraction, and then developing prediction models using the XGBoost algorithm. The SOA was used to optimize the ML model’s parameters, resulting in accurate predictions of polymer properties. SHAP and LIME algorithms were used to interpret the model, thereby allowing both local and global structure-property relationship (SPR) analysis. This paper’s results have been juxtaposed with earlier literature, such as that of (Tamasi et al., 2022), who worked on the concept of polymer-protein hybrids, while (Sobuz et al., 2025) used ML in explaining changes to concrete properties using polymer addition. The former paper offered insights into hybrid polymers, while the latter provided an approach to interpretable ML. On the other hand, (Nistane et al., 2022) worked more specifically as it involved predicting interaction parameters in polymer-solvent systems. This work differs from previous investigations in that the present study places considerable emphasis on the concept of an ML approach, along with explanation methods for understanding the effects of molecular aspects on polymer performance. This will help increase knowledge of the use of ML algorithms for both prediction and the discovery of molecular factors impacting polymer performance. The relationship between polymer properties and molecular characteristics can be understood through structure-property relationships which serve as fundamental scientific knowledge. The study identified thermal descriptors which include Tm and Td as essential factors for predicting polymer behavior while the study found that molecular flexibility and structural fragments which include ECFP features acted as secondary essential elements. The findings enable targeted polymer design which machine learning technology can use for effective material assessment. However, despite promising results, this research project has certain limitations. To begin with, it should be noted that the model’s success is largely based on the amount and quality of the data used. This limitation can reduce the model’s generalizability to other polymers or even other types of chemical reactions. Second, molecular descriptors may miss some structural features of the polymers under investigation, thereby reducing the accuracy of predictions. While SHAP and LIME techniques allow us to make the analysis of the model more interpretable, they are still approximation approaches. Moreover, the XGBoost algorithm is predominantly used in this research; the effectiveness of other algorithms (e.g., deep learning-based) was not evaluated.

The framework required evaluation through one benchmark dataset while researchers implemented multiple strategies to achieve stable results and adaptable learning of structure–property relationships. The research team used multiple random seeds together with cross-validation techniques to test model performance across various data distribution patterns. The framework achieved consistent results across different starting points because it succeeded in detecting broad patterns that applied to multiple datasets. The development of external validation through independent polymer datasets represents a critical research area which will help test the model’s ability to work in various practical environments. The proposed framework implements Recursive Feature Elimination (RFE) and SFOA-based hyperparameter optimization as its iterative processes which result in higher computational demands than standard machine learning methods. The process of feature reduction enables the system to achieve lower input dimensions which results in decreased model training requirements. The optimization process was executed offline and reached convergence through a limited number of iterations which made the framework suitable for actual use. The proposed method features a modular architecture which enables simultaneous execution and supports expansion to larger datasets through distributed computing and GPU-accelerated systems to achieve better performance in handling large polymer informatics tasks.

6. Conclusion and Future Works

The findings of this study led to the creation of a framework to develop an ML approach for the prediction of the physical properties of polymers via SPR modeling. The developed ML model used the XGBoost algorithm optimized via SOA, along with feature extraction methods based on molecular descriptors and topological indices. High predictive power was achieved, yielding R² values greater than 0.92, MAE values less than 0.08, and RMSE values less than 0.12 for selected polymer properties. In addition, applying the SHAP and LIME techniques provided even greater transparency into the molecular factors contributing to polymer properties. The integration of more advanced SPR representations, such as graph-based molecular embeddings and deep learning-driven feature extraction, could further enhance predictive performance and generalization. Expanding the descriptor space to include electronic and quantum-chemical features may provide deeper insights into polymer behavior. Additionally, combining SPR modeling with generative design approaches could enable inverse design of polymers with tailored properties. These directions highlight the potential of SPR-based machine learning frameworks to move beyond prediction toward intelligent material design and discovery. The application of such an approach has clear benefits since, along with achieving better prediction of polymer properties and increasing the effectiveness of calculations, the researcher can gain additional insights about the relationship between polymer structure and properties. Further research will include increasing the number of polymer compounds in our training set to generate generalizable results. Further research will also include adopting deep learning approaches, such as graph neural networks, to facilitate learning the structures of molecules and their properties. Finally, it is necessary to ensure the model is even more interpretable so that important attributes can be identified easily.