Direct All-Atom Molecular Dynamics Simulations of the Effects of Short Chain Branching on Polyethylene Oligomer Crystal Nucleation

Chain-Level Analysis of Reinforced Polyethylene through Stretch-Induced Crystallization

Herein, we propose a large-scale simulation approach to perform the stretch-induced crystallization of entangled polyethylene (PE) melts. Sufficiently long (1000 ns) united-atom molecular dynamics (UAMD) simulations for 16000 chains of 1000 consecutive CH2 united-atom particles under periodic boundary conditions were performed to achieve the crystallinity observed in experiments. Before the isothermal crystallization process, we applied uniaxial stretching as pre-elongation to the embedded strain memory on the entangled PE melts. We confirmed significant differences in the morphologies of crystal domains and scattering patterns for pre-elongation ratios of 400% and 800%. The obtained scattering patterns were consistent with the experimental results. Uniaxial stretching MD simulations revealed that the elastic modulus at 800% pre-elongation was stronger than that at 400% pre-elongation. From this observation, we can derive the structure-property relationship, wherein the magnitude of the pre-elongation governs the crystal domain structures and mechanical properties.

Predicting Natural Rubber Crystallinity by a Novel Machine Learning Algorithm Based on Molecular Dynamics Simulation Data

Natural rubber (NR) with excellent mechanical properties, mainly attributed to its strain-induced crystallization (SIC), has garnered significant scientific and technological interest. With the aid of molecular dynamics (MD) simulations, we can investigate the impacts of crucial structural elements on SIC on the molecular scale. Nonetheless, the computational complexity and time-consuming nature of this high-precision method constrain its widespread application. The integration of machine learning with MD represents a promising avenue for enhancing the speed of simulations while maintaining accuracy. Herein, we developed a crystallinity algorithm tailored to the SIC properties of natural rubber materials. With the data enhancement algorithm, the high evaluation value of the prediction model ensures the accuracy of the computational simulation results. In contrast to the direct utilization of small sample prediction algorithms, we propose a novel concept grounded in feature engineering. The proposed machine learning (ML) methodology consists of (1) An eXtreme Gradient Boosting (XGB) model to predict the crystallinity of NR; (2) a generative adversarial network (GAN) data augmentation algorithm to optimize the utilization of the limited training data, which is utilized to construct the XGB prediction model; (3) an elaboration of the effects induced by phospholipid and protein percentage (ω), hydrogen bond strength (εH), and non-hydrogen bond strength (εNH) of natural rubber materials with crystallinity prediction under dynamic conditions are analyzed by employing weight integration with feature importance analysis. Eventually, we succeeded in concluding that εH has the most significant effect on the strain-induced crystallinity, followed by ω and finally εNH.

Local and Global Order in Dense Packings of Semi-Flexible Polymers of Hard Spheres

The local and global order in dense packings of linear, semi-flexible polymers of tangent hard spheres are studied by employing extensive Monte Carlo simulations at increasing volume fractions. The chain stiffness is controlled by a tunable harmonic potential for the bending angle, whose intensity dictates the rigidity of the polymer backbone as a function of the bending constant and equilibrium angle. The studied angles range between acute and obtuse ones, reaching the limit of rod-like polymers. We analyze how the packing density and chain stiffness affect the chains' ability to self-organize at the local and global levels. The former corresponds to crystallinity, as quantified by the Characteristic Crystallographic Element (CCE) norm descriptor, while the latter is computed through the scalar orientational order parameter. In all cases, we identify the critical volume fraction for the phase transition and gauge the established crystal morphologies, developing a complete phase diagram as a function of packing density and equilibrium bending angle. A plethora of structures are obtained, ranging between random hexagonal closed packed morphologies of mixed character and almost perfect face centered cubic (FCC) and hexagonal close-packed (HCP) crystals at the level of monomers, and nematic mesophases, with prolate and oblate mesogens at the level of chains. For rod-like chains, a delay is observed between the establishment of the long-range nematic order and crystallization as a function of the packing density, while for right-angle chains, both transitions are synchronized. A comparison is also provided against the analogous packings of monomeric and fully flexible chains of hard spheres.

Polymer crystallization under external flow

The general aspects of polymer crystallization under external flow, i.e., flow-induced crystallization (FIC) from fundamental theoretical background to multi-scale characterization and modeling results are presented. FIC is crucial for modern polymer processing, such as blowing, casting, and injection modeling, as two-third of daily-used polymers is crystalline, and nearly all of them need to be processed before final applications. For academics, the FIC is intrinsically far from equilibrium, where the polymer crystallization behavior is different from that in quiescent conditions. The continuous investigation of crystallization contributes to a better understanding on the general non-equilibrium ordering in condensed physics. In the current review, the general theories related to polymer nucleation under flow (FIN) were summarized first as a preliminary knowledge. Various theories and models, i.e., coil-stretch transition and entropy reduction model, are briefly presented together with the modified versions. Subsequently, the multi-step ordering process of FIC is discussed in detail, including chain extension, conformational ordering, density fluctuation, and final perfection of the polymer crystalline. These achievements for a thorough understanding of the fundamental basis of FIC benefit from the development of various hyphenated rheometer, i.e., rheo-optical spectroscopy, rheo-IR, and rheo-x-ray scattering. The selected experimental results are introduced to present efforts on elucidating the multi-step and hierarchical structure transition during FIC. Then, the multi-scale modeling methods are summarized, including micro/meso scale simulation and macroscopic continuum modeling. At last, we briefly describe our personal opinions related to the future directions of this field, aiming to ultimately establish the unified theory of FIC and promote building of the more applicable models in the polymer processing.

Circular Economy Assessment in Recycling of LLDPE Bags According to European Resolution, Thermal and Structural Characterization

According to the Circular Economy Package promoted by the European directive, plastic bags companies must use in their formulations a percentage of polyethylene waste (industrial and/or domestic) greater than 70%. Following that regulation requires an understanding of its consequences in the final product from an industrial point of view. This manuscript analyzes the thermal and morphological changes related to the tear resistance of linear-low density polyethylene (LLDPE) samples from industrial waste generated by the company Sphere Spain subjected to the degradation produced by the recycling cycles. The process is analogue to the industrial, starts from samples in pellets then a film by blow extrusion is obtained (odd steps) and posteriorly this film is recycled to pellets again (even steps). The results obtained show that the LLDPE samples develop two crystalline structures (CS1 and CS2) which evolve differently through the recycling cycles with a tendency to decrease in crystallinity due to degradation that is not the same for the process of obtaining film or recycling to pellet. The molecules with a more linear structure and a longer chain break and branch. The more branched structure increases and tends to crosslinking. This leads to a decrease in tear strength in the longitudinal direction, which is not so evident in the transversal direction. The samples could admit four recycling cycles with and acceptable tear resistance. The longitudinal tear strength value decreases by 40% for each film and 20% in the case of tearing in the transverse direction. The results obtained in this research work show that the regulations included in the cited circular economy package can be applied in the manufacture of consumer bags, helping also to reduce the dependence of manufacturers on fluctuations in delivery by collapses in shipping.

Molecular dynamics simulation of extension-induced crystallization of branched bimodal HDPE: Unraveling the effects of short-chain branches

Bimodal HDPE models were designed for extension-induced crystallization imitating the architecture of industrial bimodal HDPE copolymerized with ethylene and 1-butene, 1-hexene, or 1-octene. Crystallites of bimodal HDPE experienced the emergence of precursors, shish nuclei, and lamellae. The compact conformation of branched polymers impeded the rolling-over, deposition, and folding of chains on the substrate, and thus the formation of nuclei and lamella. Moreover, this retardation was intensified with the rising branch density and length, causing a depression of crystallinity and an increment of tie-chains concentration. Besides, when branches were all located on long chains, the compact conformation enlarged the resistance to the disentanglement of main chains, thus relatively fewer branched long chains were involved in the precursors or nuclei, resulting in the attenuation of lamella formation. Furthermore, for ethyl branched polymers, the coexistent orthorhombic and monoclinic crystallites were built up, and a few expanded monoclinic cells occurred for butyl branches because of the larger butyl reeling into lamella, while hexagonal crystals were created for ethyl/1-hexyl copolymers because of cocrystallization. Additionally, relative to ethyl, larger butyl and hexyl were preferential to be repelled outside crystals to form tie-chains, and hexyl branched polymers acquired relatively fewer tie-chains because of hexagonal eutectoid.

Molecular Dynamics Simulations of Crystal Nucleation near Interfaces in Incompatible Polymer Blends

We apply molecular dynamics (MD) simulations to investigate crystal nucleation in incompatible polymer blends under deep supercooling conditions. Simulations of isothermal nucleation are performed for phase-separated blends with different degrees of incompatibility. In weakly segregated blends, slow and incompatible chains in crystallizable polymer domains can significantly hinder the crystal nucleation and growth. When a crystallizable polymer is blended with a more mobile species in interfacial regions, enhanced molecular mobility leads to the fast growth of crystalline order. However, the incubation time remains the same as that in pure samples. By inducing anisotropic alignment near the interfaces of strongly segregated blends, phase separation also promotes crystalline order to grow near interfaces between different polymer domains.

Property Decoupling across the Embryonic Nucleus-Melt Interface during Polymer Crystal Nucleation

Spatial distributions are presented that quantitatively capture how polymer properties (e.g., segment alignment, density, and potential energy) vary with distance from nascent polymer crystals (nuclei) in prototypical polyethylene melts. It is revealed that the spatial extent of nuclei and their interfaces is metric-dependent as is the extent to which nucleus interiors are solid-like. As distance from a nucleus increases, some properties, such as density, decay to melt-like behavior more rapidly than polymer segment alignment, indicating that a polymer nucleus resides in a nematic-like droplet. This nematic-like droplet region coincides with enhanced formation of ordered polymer segments that are not part of the nucleus. It is more favorable to find nonconstituent ordered polymer segments near a nucleus than in the surrounding metastable melt, pointing to the possibility of one nucleus inducing the formation of other nuclei. In this vein, there is also a second region of enhanced ordering that lies along the nematic director of a nucleus, but beyond its nematic droplet and fold regions. These results indicate that crystal stacking, a key characteristic of lamellae in semicrystalline polymeric materials, begins to emerge during the earliest stages of polymer crystallization (i.e., crystal nucleation). More generally, the findings of this study provide a conceptual bridge between polymer crystal nucleation under nonflow and flow conditions and are used to rationalize previous results.

A metastable nematic precursor accelerates polyethylene oligomer crystallization as determined by atomistic simulations and self-consistent field theory

Using PYS, TraPPE, OPLS-L, and Flexible-Williams (FW) force field models, atomistic simulations at temperatures ranging from 450 K to 600 K are performed to predict the melt density ρ, the persistence length N_p, the nematic coupling constant α, and crystallization dynamics for pentacontane (C50). The coupling constant α arises from packing entropy of rodlike Kuhn segments and increases with increasing ρ and N_p. Together with a self-consistent field theory, N_p and α are then used to predict the isotropic-to-nematic (IN) transition temperature for polyethylene (PE) oligomers as a function of chain length. The nematic phase is found to be metastable since the IN transition temperature lies below the crystal melting temperatures for C50 in simulations using different force fields. Finally, isothermal simulations of crystallization for PE C50 oligomers and C1000 polymers show that crystal nucleation may be much accelerated by quenching below the IN transition temperature, where chains in the isotropic state first rapidly form nematic ordered domains, within which crystalline order then grows. We also find that the PYS, TraPPE, and FW models overpredict the melting temperature for C50 by around 50 K, while the most flexible OPLS-L model gives a melting temperature within around 10 K of the experimental value. Although giving a more accurate melting temperature, the slow crystallization kinetics of the OPLS-L model may limit its application in direct simulations of PE crystallization.

Toward realistic computer modeling of paraffin-based composite materials: critical assessment of atomic-scale models of paraffins

Paraffin-based composites represent a promising class of materials with numerous practical applications such as e.g. heat storage. Computer modeling of these complex multicomponent systems requires a proper theoretical description of both the n-alkane matrix and the non-alkane filler molecules. The latter can be modeled with the use of a state-of-the-art general-purpose force field such as GAFF, CHARMM, OPLS-AA and GROMOS, while the paraffin matrix is traditionally described in the frame of relatively old, alkane-specific force fields (TraPPE, NERD, and PYS). In this paper we link these two types of models and evaluate the performance of several general-purpose force fields in computer modeling of paraffin by their systematic comparison with earlier alkane-specific models as well as with experimental data. To this end, we have performed molecular dynamics simulations of n-eicosane bulk samples with the use of 10 different force fields: TraPPE, NERD, PYS, OPLS-UA, GROMOS, GAFF, GAFF2, OPLS-AA, L-OPLS-AA, and CHARMM36. For each force field we calculated several thermal, structural and dynamic characteristics of n-eicosane over a wide temperature range. Overall, our findings show that the general-purpose force fields such as CHARMM36, L-OPLS-AA and GAFF/GAFF2 are able to provide a realistic description of n-eicosane samples. While alkane-specific models outperform most general-purpose force fields as far as the temperature dependence of mass density, the coefficient of volumetric thermal expansion in the liquid state, and the crystallization temperature are concerned, L-OPLS-AA, CHARMM36 and GAFF2 force fields provide a better match with experiment for the shear viscosity and the diffusion coefficient in melt. Furthermore, we show that most general-purpose force fields are able to reproduce qualitatively the experimental triclinic crystal structure of n-eicosane at low temperatures.

Atomic-scale computational models of paraffins are critically assessed and compared.