Computational Screening of Thermoplastic Starch-Based Biopolymers Blends for Mixing Compatibility ()
1. Introduction
We are witnessing, in scientific research, a rising trend focused on the precise and systematic rationalization of various computational data into legitimate chemical investigations. Polymers are ubiquitous materials that play critical roles in various industries, from packaging [1] to biomedical devices [2]. Their properties can be tailored through molecular design, making them highly versatile materials. However, the vastness of potential polymer candidates complicates the material selection and optimization process. Computational screening emerges as a powerful tool to accelerate the discovery and development of new polymers by leveraging computational methods to predict polymer properties and behaviors before experimental validation [3]. The synthesis and characterization of polymers often require significant time and expense. Traditional methods involve trial-and-error approaches, which can be inefficient and slow. Computational screening allows researchers to explore a vast chemical space quickly and cost-effectively [4]. By simulating polymer properties and performance, computational tools can help identify promising candidates, optimize their synthesis routes, and reduce the need for extensive laboratory experimentation [5]. Various studies focusing on computational screening of Thermoplastic starch (a starch-based material obtained by plasticizing starch with agents such as glycerol) and its polymeric blends for various properties have been released lately. Ozeren et al. successfully investigated the plasticization process of starch with glycerol. Their computational approach can be used for the analysis of other polymer/plasticizer systems as well as provide a new route for the development of novel classes of plasticizers [6]. Two years later, Ozeren et al. managed to extrapolate simulation data for prediction of the tensile properties of thermoplastic starch. The accuracy with respect to the experimental data projection was acceptable in many situations; however, it depended on the extrapolation approach as well as on the particular TPS system [7]. Furthermore, Yamaguchi et al. assessed the toughening caused by the development of a core-shell molecular arrangement in starch-based biopolymer composites experimentally and computationally. The researchers performed molecular dynamics simulations for the examination of mechanisms linked to the mechanical and structural characteristics of ternary blends of poly(lactic acid) (PLA), poly(butylene succinate) (PBS), and thermoplastic starch (TPS), as well as of binary blends of PLA/TPS. The dynamics simulation results revealed that PLA/PBS/ TPS blends showed improved mechanical properties compared with PLA/TPS blends. The PLA/PBS/TPS blends with a TPS ratio of 25 - 40 wt% showed higher impact strength than the PLA/PBS blend [8]. Subsequently, TPS plasticization has been studied through both experimental and computational approaches. Results have shown that starch-based systems processed at 140˚C with 15% glycerol (as plasticizer) displayed better tensile properties and wettability. The thermal stability, however, lessened in correlation with an increase in glycerol content. [9]. Finally, Yamaguchi et al. executed a study of the phase-separated interface of four environment-friendly polyesters, such as poly(lactic acid), polybutylene succinate (PBS), polyhydroxybutyrate (PHB), and poly(butylene adipate-co-terephthalate), with TPS using molecular dynamics simulations; it was concluded that PBS, PHB, and PBAT easily spread over TPS and constitute well-matched interfaces, whereas PLA displays lesser compatibility with the TPS. Both PLA and PHB associated with TPS show poor interfacial-fracture energy due to the poor compatibility between PLA and TPS, as well as to PHB having low electrostatic interaction with TPS [10]. The aim of this investigation is to assess the mixing compatibility of TPS-based binary blends (TPS/NR and TPS/PCL) through quick computational molecular screening. This study will also examine the potential accuracy of the computational approach based on the Flory-Huggins mixing theory.
2. Theory and Material
2.1. Mixtures Compatibility
The compatibility characteristics of TPS-based polymer blends were computed with the help of the Blends module. The screening of monomers, facilitated by the Blends module, provides a simplified yet effective approach, as the interactions between monomers serve as indicators of full polymer chain compatibility within this computational framework, allowing efficient prediction of polymer blend behavior [11]. The principle behind the Blends module calculations is expressed by a combination of the Flory-Huggins model and methods of simulations applied to molecules [12]. The ultimate aim is to calculate the compatibility of binary mixtures. These binary systems are evaluated in all sizes, ranging from small systems to relatively huge molecular models. Blends module algorithms establish the distinction between the two constituents of the system by assigning individual roles: one component is assigned the base (b) role, and the screen (s) role is assigned to the second [13].
2.2. Binding Energies
When two molecules are in relative spatial proximity, a phenomenon of intermolecular interaction will take place. The quantity of interaction energy between the two chemicals is termed binding energy. Blends algorithms combine binding energy parameters with coordination numbers during the compatibility investigation to generate the mixing energy (Emix) and chi parameter [14]. Estimation of binding energies requires careful consideration and judgment of the high molecular orientation applied to the two constituents during screening [15]. An extensive number of molecular orientations is provided by the module, which subsequently estimates the energies of interaction related to every configuration. After screening for the binding energy, the Blends algorithm executes the sampling of those configurations considered to be energetically satisfactory [16]. The sampling method utilized here is an altered version of the molecular Silverware algorithm [17]: it is called excluded volume constraints. Excluded volume constraint sampling methods are employed to prevent spatial overlap between particles or molecules [18]. In this approach, a proposed configuration is generated and checked for particle overlap; configurations exhibiting overlaps are rejected, and new ones are sampled [19]. Here, the constraint is applied based on the excluded volume associated with each molecular configuration. In this sampling technique, the two molecules are joined together in order to establish a proximity of their van der Waals surfaces. Four distinct possible pairs will be generated after any base-screen combination. The binding energy value of each individual pair is expressed as:
Base-base pair (Ebb)
Screen-screen pair (Ess)
Base-screen pair (Ebs)
Screen-base pair (Esb)
Let us consider a binary mixture made of component A and component B. The binding energy of the interaction between those molecules is calculated using the excluded-volume constraint method. In this method, the calculation of the binding energy is initiated with two chemical compounds that are a representation of the condensed state [20]. After the estimation of binding energies is completed, the module automatically establishes the coordination numbers as well [21]. Once those parameters have been generated, the Blends algorithm performs the calculations of the mixing energy (Emix) based on the equations derived from the Flory-Huggins model.
2.3. Flory-Huggins Model
The efforts to understand the thermodynamics of polymeric solutions are represented by the Flory-Huggins theory. The theory provides a more precise investigation compared to the ideal solution theory. The ideal solution theory presents various conditions of limited ability [22]. Two of the most important limitations of the theory are the large size of polymeric systems compared to solvents, as well as the possibility of chemical interaction between different molecules. The Flory-Huggins model stands out as, without a doubt, the best and least complicated theory for the assessment of binary mixtures thermodynamics of phase separation and mixing [23].
According to the model, the free energy of mixing of the binary system is expressed by the general equation:
(1)
where: ΔG = free energy of mixing per mole;
Φb = volume fraction of component A;
Φs = volume fraction of component B;
Χ = interaction parameter;
T = absolute temperature;
R = gas constant.
The combinatorial entropy is expressed by the first two elements. Due to their negative mathematical contribution, they provide favorable conditions for a mixed state, which is preferred over the pure components [24]. The term at the end expresses the free energy generated due to interaction. A reliable sign of the attainment of a mixed state is the negative value of the chi interaction parameter.
If the parameter of interaction is positive, a mixed state is inferred to be disfavored [25]. Moreover, when there is equilibrium between both conditions, the algorithm will generate different phase diagrams. The following equation represents the χ interaction parameter:
(2)
where Emix = mixing energy.
Mixing energy is defined as the free energy difference due to interactivity occurring between the pure state and blended state. The classical Flory-Huggins model suggests that each constituent molecule occupies a geometric arrangement of the points in space at which the atoms or ions of a crystal occur, i.e., a lattice site. Supposing that a lattice site has a coordination number Z, the mixing energy is defined by:
. [26] (3)
where, Ebb = Binding energy of Base-base pair;
Ess = Binding energy of Screen-screen pair;
Ebs = Binding energy of Base-screen pair;
Z = Coordination number for each pair (which is the number of molecules of component B that can be packed around a single molecule of component A within the excluded-volume constraints).
A positive sign of the attainment of a mixed state is the negative value of the chi interaction parameter. If the parameter of interaction is positive, the mixed state is inferred to be disfavored. Moreover, when there is an equilibrium between both conditions, the algorithm will generate different phase diagrams. The combinatorial entropy is expressed by the first two elements. Due to their negative mathematical contribution, they provide favorable conditions for a mixed state, which is preferred over the pure components. The term at the end expresses the free energy generated due to interaction.
3. Calculations
3.1. Geometry Optimization
Three TPS-based polymer mixtures were selected for the compatibility analysis: TPS/NR and TPS/PCL. To begin, the individual structures of the constituting monomers were manually sketched on atomistic 3D windows. Subsequently, energy minimization of sketched structures was carried out through geometry optimization with 500 iteration cycles utilizing the COMPASS force field [27] and SMART algorithm. Energy minimization was conducted using the Forcite module, producing new 3D molecular structures with optimized geometries illustrated in Figures 1(a)-(c), where hydrogen, oxygen, and carbon atoms are depicted as white, red, and grey spheres, respectively.
(a) (b)
(c)
Figure 1. Geometry optimized structure: (a) Thermoplastic starch (amylose + glycerol) [28]; (b) Natural rubber monomer (isoprene); (c) Polycaprolactone monomer.
3.2. Compatibility Calculations
Calculations were executed with the help of the “blends” module. The blends module allows the screening of interaction energy between polymer and polymer, between polymer and solvent, as well as between solvent and solvent. The selection of the monomer to be studied was performed by specifying its head and tail atoms, followed by structure optimization before submission to the calculation tab.
The selected monomers were entered in the first row, and their roles—screen or base—were specified in the second. An overview of the calculation tab layout is provided in Figure 2. In this work, TPS served as the base in every blend investigation. The reference temperature was kept at 298 K. Calculation outputs of the mixing task include mixing energy (Emix), binding energies (for the base-base, base-screen, and screen-screen combinations), and χ chi values. The module task of mixing provides the user with an accurate sampling of the mixing characteristics of the selected binary blends. Blends allow the screening for interaction energy between polymer and polymer, between polymer and solvent, as well as between solvent and solvent. For further extrapolations, the blending characteristics of TPS/Natural rubber (NR) and TPS/Polycaprolactone (PCL) binary mixtures of monomers were recorded.
Figure 2. Blends calculation tab in Material Studio 5.5 (Accelrys, Inc.)
4. Results and Discussion
4.1. Polymer Compatibility
For the data analysis of mixing effectiveness of TPS blends, the following parameters have captured our attention: binding energies distributions (Ebs, Ebb, Ess), Chi parameter χ, and energy of mixing (Emix).
4.2. Binding Energies Distribution
A strong indicator of structural compatibility between the input monomer pair is the similarity among the binding energy distribution curves for base-base (Ebb), base-screen (Ebs), and screen-screen (Ess) interactions [29].
(a)
(b)
Figure 3. Plot of binding energy distribution per energy sample P(E): (a) TPS/NR blend; (b) TPS/PCL blend.
From the plots in Figure 3, a stronger overlap of binding energy distribution curves is observed for the TPS/PCL mixture (Figure 3(b)) relative to the TPS/NR system (Figure 3(a)). This suggests that TPS/PCL possesses better structural compatibility, whereas TPS/NR shows the least. This is justified by the fact that TPS/PCL forms compatible blends because of the hydrogen bonding interaction between the ester carbonyl of PCL and the OH groups on starch [30]. Furthermore, it has been observed that the addition of NR to TPS is restricted by phase separation due to poor compatibility between both constituents [31].
4.3. Chi Parameter χ and Mixing Energy (Emix)
The χ interaction parameter is a unitless property expressing the difference in energy between solvent-polymer interaction and a polymer-free solvent. Its value can be positive or negative [32]. A χ (chi) value close to zero indicates miscibility, as does an Emix value that is close to zero [33]. Lower χ (chi) and mixing energy values for the TPS/PCL system (Table 1) reflect higher compatibility. These results further confirm what has been implied from the binding energies distribution analysis.
Table 1. Calculated chi parameter and Emix of TPS/NR and TPS/PCL blends.
Polymer Blend |
χ Interaction Parameter |
Mixing Energy (kcal·mol−1) |
TPS/NR |
5.381 |
3.186 |
TPS/PCL |
4.429 |
2.623 |
Configurations with the lowest interaction energies
Through the Blends module, pair configurations (base-screen, base-base, and screen-screen) can be organized based on energy to inspect the characteristic binding and non-binding (non-contact) sites. This is demonstrated in Figure 4(a) and Figure 4(b), where non-binding sites are denoted as NON-CONTACT.
(a)
(b)
Figure 4. (a) Configurations with the lowest energy of interaction for the TPS/NR system, TPS/NR blend; (b) Configurations with the lowest energy of interaction for the TPS/PCL system, TPS/PCL blend.
5. Conclusions
A molecular screening of the interaction energy between the monomer pairs of TPS/NR and TPS/PCL was successfully performed for studying polymer structural compatibility. From the collected data, caprolactone had the best mixing capacity with TPS, whereas natural rubber had the lowest capacity. Poor compatibility of the TPS/NR blend has been reported in existing literature, where it was observed that natural rubber reinforcement of TPS is restricted by phase separation due to poor compatibility of the constituents. This study suggested that an increase in plasticizer (glycerol) content can help to overcome the problem of phase separation. The exact extent to which the use of plasticizer can ameliorate the TPS and NR compatibility is yet to be investigated. TPS/PCL displayed the best mixing potential; the hydrogen bonding interaction that is observed between the ester carbonyl of the caprolactone monomer and the hydroxyl groups present in the starch glucose monomer has explained TPS/PCL compatibility [34]. This high degree of chemical compatibility is also made possible by the elevated content of plasticizer, as well as a very efficient process of TPS plasticization.
Owing to the superior compatibility of the TPS/PCL blend, a diverse range of TPS/PCL-based materials can be developed for applications such as compostable packaging, biodegradable agricultural films, and controlled-release fertilizer coatings [35] [36]. The Flory-Huggins model of investigation for polymeric blending compatibility is a rapid and less costly computational approach compared to the traditional Molecular Dynamics (MD) simulation (although MD simulation might provide a more detailed perspective). The use of monomers as proxies for complete polymer chains, along with the limitations of the selected force field, imposes constraints on the present study. Such approximations should be considered when evaluating the results and planning subsequent studies. The Flory-Huggins approach can facilitate a preliminary comparative analysis of one or several binary blends of polymers by a quick molecular screening of the individual monomers constituting the mixtures.
Acknowledgements
My heartfelt thanks go to my wife, Kahindo Rosette, for her love, patience, and endless support during this journey. Her encouragement has been a constant source of strength. I am equally grateful to my colleagues for their teamwork, valuable input, and helpful discussions, which have played a key role in the success of this work.