Predicting glass transition of amorphous polymers by application of cheminformatics and molecular dynamics simulations

Phenolic Acid-β-Cyclodextrin Complexation Study to Mask Bitterness in Wheat Bran: A Machine Learning-Based QSAR Study

The need to solvate and encapsulate hydro-sensitive molecules drives noticeable trends in the applications of cyclodextrins in the pharmaceutical industry, in foods, polymers, materials, and in agricultural science. Among them, β-cyclodextrin is one of the most used for the entrapment of phenolic acid compounds to mask the bitterness of wheat bran. In this regard, there is still a need for good data and especially for a robust predictive model that assesses the bitterness masking capabilities of β-cyclodextrin for various phenolic compounds. This study uses a dataset of 20 phenolic acids docked into the β-cyclodextrin cavity to generate three different binding constants. The data from the docking study were combined with topological, topographical, and quantum-chemical features from the ligands in a machine learning-based structure-activity relationship study. Three different models for each binding constant were computed using a combination of the genetic algorithm (GA) and multiple linear regression (MLR) approaches. The developed ML/QSAR models showed a very good performance, with high predictive ability and correlation coefficients of 0.969 and 0.984 for the training and test sets, respectively. The models revealed several factors responsible for binding with cyclodextrin, showing positive contributions toward the binding affinity values, including such features as the presence of six-membered rings in the molecule, branching, electronegativity values, and polar surface area.

Interpretable Machine Learning Framework to Predict the Glass Transition Temperature of Polymers

The glass transition temperature of polymers is a key parameter in meeting the application requirements for energy absorption. Previous studies have provided some data from slow, expensive trial-and-error procedures. By recognizing these data, machine learning algorithms are able to extract valuable knowledge and disclose essential insights. In this study, a dataset of 7174 samples was utilized. The polymers were numerically represented using two methods: Morgan fingerprint and molecular descriptor. During preprocessing, the dataset was scaled using a standard scaler technique. We removed the features with small variance from the dataset and used the Pearson correlation technique to exclude the features that were highly connected. Then, the most significant features were selected using the recursive feature elimination method. Nine machine learning techniques were employed to predict the glass transition temperature and tune their hyperparameters. The models were compared using the performance metrics of mean absolute error (MAE), root mean square error (RMSE), and coefficient of determination (R2). We observed that the extra tree regressor provided the best results. Significant features were also identified using statistical machine learning methods. The SHAP method was also employed to demonstrate the influence of each feature on the model's output. This framework can be adaptable to other properties at a low computational expense.

Quantitative Structure-Permittivity Relationship Study of a Series of Polymers

Dielectric constant is an important property which is widely utilized in many scientific fields and characterizes the degree of polarization of substances under the external electric field. In this work, a structure-property relationship of the dielectric constants (ε) for a diverse set of polymers was investigated. A transparent mechanistic model was developed with the application of a machine learning approach that combines genetic algorithm and multiple linear regression analysis, to obtain a mechanistically explainable and transparent model. Based on the evaluation conducted using various validation criteria, four- and eight-variable models were proposed. The best model showed a high predictive performance for training and test sets, with R2 values of 0.905 and 0.812, respectively. Obtained statistical performance results and selected descriptors in the best models were analyzed and discussed. With the validation procedures applied, the models were proven to have a good predictive ability and robustness for further applications in polymer permittivity prediction.

Relationships between the Photodegradation Reaction Rate and Structural Properties of Polymer Systems

The development of reusable polymeric materials inspires an attempt to combine renewable biomass with upcycling to form a biorenewable closed system. It has been reported that 2,5-furandicarboxylic acid (FDCA) can be recovered for recycling when incorporated as monomers into photodegradable polymeric systems. Here, we develop a procedure to better understand the photodegradation reactions combining density functional theory (DFT) based time-dependent excited-state molecular dynamics (TDESMD) studies with machine learning-based quantitative structure-activity relationships (QSAR) methodology. This procedure allows for the unveiling of hidden structural features between active orbitals that affect the rate of photodegradation and is coined InfoTDESMD. Findings show that electrotopological features are influential factors affecting the rate of photodegradation in differing environments. Additionally, statistical validations and knowledge-based analysis of descriptors are conducted to further understand the structural features' influence on the rate of photodegradation of polymeric materials.

Combined Machine Learning, Computational, and Experimental Analysis of the Iridium(III) Complexes with Red to Near-Infrared Emission

Various coordination complexes have been the subject of experimental and theoretical studies in recent decades because of their fascinating photophysical properties. In this work, a combined experimental and computational approach was applied to investigate the optical properties of monocationic Ir(III) complexes. An interpretative machine learning-based quantitative structure-property relationship (ML/QSPR) model was successfully developed that could reliably predict the emission wavelength of the Ir(III) complexes and provide a foundation for the theoretical evaluation of the optical properties of Ir(III) complexes. A hypothesis was proposed to explain the differences in the emission wavelengths between structurally different individual Ir(III) complexes. The efficacy of the developed model was demonstrated by high R2 values of 0.84 and 0.87 for the training and test sets, respectively. It is worth noting that a relationship between the N-N distance in the diimine ligands of the Ir(III) complexes and emission wavelengths is detected. This effect is most probably associated with a degree of distortion in the octahedral geometry of the complexes, resulting in a perturbed ligand field. This combined experimental and computational approach shows great potential for the rational design of new Ir(III) complexes with the desired optical properties. Moreover, the developed methodology could be extended to other transition-metal complexes.

A Molecular Dynamics Simulation Study on Enhancement of Mechanical and Tribological Properties of Nitrile-Butadiene Rubber with Varied Contents of Acrylonitrile

The molecular models of nitrile-butadiene rubber (NBR) with varied contents of acrylonitrile (ACN) were developed and investigated to provide an understanding of the enhancement mechanisms of ACN. The investigation was conducted using molecular dynamics (MD) simulations to calculate and predict the mechanical and tribological properties of NBR through the constant strain method and the shearing model. The MD simulation results showed that the mechanical properties of NBR showed an increasing trend until the content of ACN reached 40%. The mechanism to enhance the strength of the rubber by ACN was investigated and analyzed by assessing the binding energy, radius of gyration, mean square displacement, and free volume. The abrasion rate (AR) of NBR was calculated using Fe-NBR-Fe models during the friction processes. The wear results of atomistic simulations indicated that the NBR with 40% ACN content had the best tribological properties due to the synergy among appropriate polarity, rigidity, and chain length of the NBR molecules. In addition, the random forest regression model of predicted AR, based on the dataset of feature parameters extracted by the MD models, was developed to obtain the variable importance for identifying the highly correlated parameters of AR. The torsion-bend-bend energy was obtained and used to successfully predict the AR trend on the new NBR models with other acrylonitrile contents.

Ring Repeating Unit: An Upgraded Structure Representation of Linear Condensation Polymers for Property Prediction

Unique structure representation of polymers plays a crucial role in developing models for polymer property prediction and polymer design by data-centric approaches. Currently, monomer and repeating unit (RU) approximations are widely used to represent polymer structures for generating feature descriptors in the modeling of quantitative structure-property relationships (QSPR). However, such conventional structure representations may not uniquely approximate heterochain polymers due to the diversity of monomer combinations and the potential multi-RUs. In this study, the so-called ring repeating unit (RRU) method that can uniquely represent polymers with a broad range of structure diversity is proposed for the first time. As a proof of concept, an RRU-based QSPR model was developed to predict the associated glass transition temperature (T_g) of polyimides (PIs) with deterministic values. Comprehensive model validations including external, internal, and Y-random validations were performed. Also, an RU-based QSPR model developed based on the same large database of 1321 PIs provides nonunique prediction results, which further prove the necessity of RRU-based structure representation. Promising results obtained by the application of the RRU-based model confirm that the as-developed RRU method provides an effective representation that accurately captures the sequence of repeat units and thus realizes reliable polymer property prediction by data-driven approaches.

Understanding the role of cross-link density in the segmental dynamics and elastic properties of cross-linked thermosets

Cross-linking is known to play a pivotal role in the relaxation dynamics and mechanical properties of thermoset polymers, which are commonly used in structural applications because of their light weight and inherently strong nature. Here, we employ a coarse-grained (CG) polymer model to systematically explore the effect of cross-link density on basic thermodynamic properties as well as corresponding changes in the segmental dynamics and elastic properties of these network materials upon approaching their glass transition temperatures (T_g). Increasing the cross-link density unsurprisingly leads to a significant slowing down of the segmental dynamics, and the fragility K of glass formation shifts in lockstep with T_g, as often found in linear polymer melts when the polymer mass is varied. As a consequence, the segmental relaxation time τ_α becomes almost a universal function of reduced temperature, (T - T_g)/T_g, a phenomenon that underlies the applicability of the "universal" Williams-Landel-Ferry (WLF) relation to many polymer materials. We also test a mathematical model of the temperature dependence of the linear elastic moduli based on a simple rigidity percolation theory and quantify the fluctuations in the local stiffness of the network material. The moduli and distribution of the local stiffness likewise exhibit a universal scaling behavior for materials having different cross-link densities but fixed (T - T_g)/T_g. Evidently, T_g dominates both τ_α and the mechanical properties of our model cross-linked polymer materials. Our work provides physical insights into how the cross-link density affects glass formation, aiding in the design of cross-linked thermosets and other structurally complex glass-forming materials.

Synergistic Behavior of Plant Proteins and Biobased Latexes in Bioplastic Food Packaging Materials: Experimental and Machine Learning Study

Plant-based proteins are attractive components which may serve as sustainable alternatives to current petrochemical products. Both soy protein and major corn protein, zein, are of interest in food packaging applications due to their sustainability, biodegradation properties, and inherent physicochemical properties. This study discusses the development of bioplastic materials, where it explores the effects of combining zein, soy protein, and plasticizing latexes derived from plant oil-based monomers (POBMs) on properties of resulting bioplastic films. By looking for synergistic effects of soy protein's inherent film formation ability and zein's higher strength, we prepare strong yet flexible soy-zein films as materials, called proteoposites. Incorporation of natural additive POBM-latexes helps to plasticize and hydrophobize the bioplastic films and thus to improve mechanical and barrier properties. Variation of the POBM-latexes' particle size further aims to enhance the performance of resulting bioplastic films. As a result, modified soy-zein proteoposite films with improved moisture resistance, enhanced mechanical behavior, and greater barrier properties were developed. Machine learning-based computational models were utilized in order to find main structural factors affecting the bioplastic's properties and develop a quantitative structure-property relationship model between the physicochemical properties of the film components and the resulted bioplastics' properties and performance. The developed model effectively predicts experimental outcomes with >85% (R²: 0.85) accuracy. The newly synthesized proteoposites confirmed the machine learning model predictions. As a result, proteoposite films made of two plant proteins and modified with POBM-latexes can be considered as an attractive and viable replacement for petrochemical food packaging products.

M SB, George JJ, Rasulev B, Raghavan P. Experimental and Simulation Studies on Nonwoven Polypropylene-Nitrile Rubber Blend: Recycling of Medical Face Masks to an Engineering Product

The battle against the COVID-19 pandemic counters the waste management system, as billions of single-use face masks are used per day all over the world. Proper disposal of used face masks without jeopardizing the health and the environment is a challenge. Herein, a novel method for recycling of medical face masks has been studied. This method incorporates the nonwoven polypropylene (PP) fiber, which is taken off from the mask after disinfecting it, with acrylonitrile butadiene rubber (NBR) using maleic anhydride as the compatibilizer, which results in a PP-NBR blend with a high percentage economy. The PP-NBR blends show enhanced thermomechanical properties among which, 70 wt % PP content shows superior properties compared to other composites with 40, 50, and 60 wt % of PP. The fully Atomistic simulation of PP-NBR blend with compatibilizer shows an improved tensile and barrier properties, which is in good agreement with the experimental studies. The molecular dynamics simulation confirms that the compatibility between non-polar PP and polar NBR phases are vitally important for increasing the interfacial adhesion and impeding the phase separation.

An Enhanced Scheme for Multiscale Modeling of Thermomechanical Properties of Polymer Bulks

While multiscale modeling significantly enhances the capability of molecular simulations of polymer systems, it is well realized that the systematically derived coarse-grained (CG) models generally underestimate the thermomechanical properties. In this work, a charge-based mapping scheme has been adopted to include explicit electrostatic interactions and benchmarked against two typical polymers, atactic poly(methyl methacrylate) (PMMA) and polystyrene (PS). The CG potentials are parameterized against the oligomer bulks of nine monomers per chain to match the essential structural features and the two basic pressure-volume-temperature (PVT) properties, which are obtained from the all-atomistic (AA) molecular dynamics (MD) simulations at a single elevated temperature. The so-parameterized CG potentials are extended with the MD method to simulate the two polymer bulks of one hundred monomers per chain over a wide temperature range. Without any scaling, all the simulated results, including mass densities and bulk moduli at room temperature, thermal expansion coefficients at rubbery and glassy states, and glass transition temperatures (T_g), compare well with the corresponding experimental data. The proposed scheme not only contributes to realistically simulating various thermomechanical properties of both apolar and polar polymers but also allows for directly simulating their electrical properties.