Effect of Short Chain Branching on the Interlamellar Structure of Semicrystalline Polyethylene

Mechanisms of Shock Dissipation in Semicrystalline Polyethylene

Semicrystalline polymers are lightweight, multiphase materials that exhibit attractive shock dissipation characteristics and have potential applications as protective armor for people and equipment. For shocks of 10 GPa or less, we analyzed various mechanisms for the storage and dissipation of shock wave energy in a realistic, united atom (UA) model of semicrystalline polyethylene. Systems characterized by different levels of crystallinity were simulated using equilibrium molecular dynamics with a Hugoniostat to ensure that the resulting states conform to the Rankine-Hugoniot conditions. To determine the role of structural rearrangements, order parameters and configuration time series were collected during the course of the shock simulations. We conclude that the major mechanisms responsible for the storage and dissipation of shock energy in semicrystalline polyethylene are those associated with plastic deformation and melting of the crystalline domain. For this UA model, plastic deformation occurs primarily through fine crystallographic slip and the formation of kink bands, whose long period decreases with increasing shock pressure.

Unraveling the influential mechanism of short-chain branching on the crystallization of trimodal polyethylene by molecular dynamics simulation

Trimodal polyethylene (PE) has become the focus of research in recent years due to its excellent performance. By means of molecular dynamics simulations, we aim to expound the molecular mechanism of short-chain branching (SCB) in the nucleation process, crystallization process and chain entanglement of trimodal PE. In this study, a series of polyethylene models including different short-chain branching concentrations (SCBCs), short-chain branching lengths (SCBLs), and short-chain branching distributions (SCBDs) were considered. The increase of SCBCs greatly reduces the ability of flipping and movement of PE chains, resulting in more time for nucleation and crystallization and a significant reduction of crystallinity. In contrast, an increase in the SCBL only slightly slows down the diffusion rate of the chain, which leads to a little increase in crystallization time. Most important of all, in the study of SCBD, we find that the distribution of SCBs on a high molecular weight chain, which is the characteristic of trimodal PE, is conducive to the chain entanglement and prevents the occurrence of micro phase separation compared with the case where the SCBs are distributed on a medium molecular weight chain. The mechanism of chain entanglement is proposed to explain the effect of SCBs on tie chain entanglement.

Atomistic Modeling of Hydrogen and Oxygen Solubility in Semicrystalline PA-6 and HDPE Materials

Hydrogen is a clean and sustainable energy carrier which plays a major role in the transition of the global energy market to a less fossil fuel dependent future. Polymer-based materials are crucial in the production, storage, transportation, and energy extraction of hydrogen. More insights in the hydrogen-polymers interactions are required to guide material design and product development, especially for hydrogen solubility in polymers, which is crucial in many applications. The current study aims at rationalizing the determining factors of hydrogen solubility in two relevant polymers: polyamide-6 (PA-6) and high density polyethylene (HDPE). Based on atomistic molecular dynamics simulations and experimental data, we have reached several conclusions related to hydrogen and oxygen solubility in these two polymers. The crystal phases of PA-6 and HDPE are impenetrable to hydrogen and oxygen at elevated pressures, despite the small molecular size of hydrogen and oxygen. The practical implication for gas barrier applications is that polymer crystals act as impermeable obstacles and gas migration takes place primarily in the amorphous phase. Experimental hydrogen and oxygen solubilities in PA-6 and HDPE at elevated pressures can be predicted in a semiquantitative manner by molecular simulations. The discrepancies between experimental and predicted values could be attributed to neglect of the effect of crystal regions on the amorphous polymer domains. Although hydrogen is smaller than oxygen, it has been experimentally observed that hydrogen has a lower solubility in PA-6 and HDPE than oxygen. This observation has been confirmed by molecular simulations and attributed to the more favorable energetic interactions of oxygen with PA-6 and PE than of hydrogen. These interactions dominate the solubility behavior over the distribution of the accessible volume in the polymers.

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.

Simu-D: A Simulator-Descriptor Suite for Polymer-Based Systems under Extreme Conditions

We present Simu-D, a software suite for the simulation and successive identification of local structures of atomistic systems, based on polymers, under extreme conditions, in the bulk, on surfaces, and at interfaces. The protocol is built around various types of Monte Carlo algorithms, which include localized, chain-connectivity-altering, identity-exchange, and cluster-based moves. The approach focuses on alleviating one of the main disadvantages of Monte Carlo algorithms, which is the general applicability under a wide range of conditions. Present applications include polymer-based nanocomposites with nanofillers in the form of cylinders and spheres of varied concentration and size, extremely confined and maximally packed assemblies in two and three dimensions, and terminally grafted macromolecules. The main simulator is accompanied by a descriptor that identifies the similarity of computer-generated configurations with respect to reference crystals in two or three dimensions. The Simu-D simulator-descriptor can be an especially useful tool in the modeling studies of the entropy- and energy-driven phase transition, adsorption, and self-organization of polymer-based systems under a variety of conditions.

Critical Thickness of Free-Standing Nanothin Films Made of Melted Polyethylene Chains via Molecular Dynamics

The mechanical stability of nanothin free-standing films made of melted polyethylene chains was predicted via molecular dynamics simulations in the range of 373.15-673.15 K. The predicted critical thickness, tc, increased with the square of the temperature, T, with additional chains needed as T increased. From T = 373.15 K up to the thermal limit of stability for polyethylene, tc values were in the range of nanothin thicknesses (3.42-5.63 nm), which approximately corresponds to 44-55 chains per 100 nm2. The density at the center of the layer and the interfacial properties studied (density profiles, interfacial thickness, and radius of gyration) showed independence from the film thickness at the same T. The polyethylene layer at its tc showed a lower melting T (<373.15 K) than bulk polyethylene.

Network structure of swollen iodine-doped poly(vinyl alcohol) amorphous domain as characterized by low field NMR

The network structure in the amorphous domain of swollen iodine-doped poly(vinyl alcohol) (PVA) was systematically investigated by low-field (LF) NMR techniques to reveal the PVA-iodine complex formation mechanism. Three PVA-iodine complexes were obtained under different iodine concentrations (c_iodine) of KI/I₂ solution: (i) c_iodine < 0.1 M: PVA-I₃^-/I₅^- complex only exists in the non-crystalline region, (ii) 0.1 M < c_iodine < 1 M: formation of PVA-I₃^- complex I, and (iii) c_iodine > 1 M: formation of PVA-I₃^- complex II. It was found that there is no intermediate-magnitude chain motion of PVA under dyeing conditions to induce the substance exchange, as evidenced by the unchanged second moment M₂ (∼1.2 × 10⁴ m s^-2) at elevated temperature (<380 K). The introduction of iodine ions can affect the chain mobility of the interphase and mobile regions. With increasing c_iodine, the chain dynamics become more restricted, as detected by the faster decay of the T₂ relaxometry results, which further accelerates the complexation process. The residual dipolar coupling strength, D_res, obtained by the more quantitative double-quantum (DQ) NMR, increases abruptly at c_iodine > 1 M. This suggests more constraints form in the amorphous network for the PVA-I₃^- complex II system. The constant defects fraction further reveals that the complexation prefers to happen along the tie chains. These results supply a possible formation pathway for the PVA-iodine complexes.

Low-Temperature Mechanical Properties of High-Density and Low-Density Polyethylene and Their Blends

Low-temperature properties of high-density polyethylene (HDPE), low-density polyethylene (LDPE), linear low-density polyethylene (LLDPE), and their blends were studied. The analyzed low-temperature mechanical properties involve the deformation resistance and impact strength characteristics. HDPE is a bimodal ethylene/1-hexene copolymer; LDPE is a branched ethylene homopolymer containing short-chain branches of different length; LLDPE is a binary ethylene/1-butene copolymer and an ethylene/1-butene/1-hexene terpolymer. The samples of copolymers and their blends were studied by gel permeation chromatography (GPC), differential scanning calorimetry (DSC), ¹³С NMR spectroscopy, and dynamic mechanical analysis (DMA) using testing machines equipped with a cryochamber. It is proposed that such parameters as "relative elongation at break at -45 °C" and "Izod impact strength at -40 °C" are used instead of the ductile-to-brittle transition temperature to assess frost resistance properties because these parameters are more sensitive to deformation and impact at subzero temperatures for HDPE. LLDPE is shown to exhibit higher relative elongation at break at -45 °C and Izod impact strength at -20 ÷ 60 °C compared to those of LDPE. LLDPE terpolymer added to HDPE (at a content ≥ 25 wt.%) simultaneously increases flow properties and improves tensile properties of the blend at -45 °C. Changes in low-temperature properties as a function of molecular weight, MWD, crystallinity, and branch content were determined for HDPE, LLDPE, and their blends. The DMA data prove the resulting dependences. The reported findings allow one to understand and predict mechanical properties in the HDPE-LLDPE systems at subzero temperatures.

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.

Dominant Effects of Short-Chain Branching on the Initial Stage of Nucleation and Formation of Tie Chains for Bimodal Polyethylene as Revealed by Molecular Dynamics Simulation

The molecular mechanism of short-chain branching (SCB), especially the effects of methylene sequence length (MSL) and short-chain branching distribution (SCBD) on the initial stage of nucleation, the crystallization process, and particularly the tie chain formation process of bimodal polyethylene (BPE), were explored using molecular dynamics simulation. This work constructed two kinds of BPE models in accordance with commercial BPE pipe resins: SCB incorporated in the long chain or in the short chains. The initial stage of nucleation was determined by the MSL of the system, as the critical MSL for a branched chain to nucleate is about 60 CH2. SCB incorporated in the long chain led to a delay of the initial stage of nucleation relative to the case of SCB incorporated in the short chains. The increase of branch length could accelerate the delay to nucleation. The location of short chain relative to the long chain depended on the MSL of the short chain. As the MSL of the system decreased, the crystallinity decreased, while the tie chains concentration increased. The tie chains concentration of the BPE model with branches incorporated in the long chain was higher than that with branches incorporated in the short chain.

Steric "attraction": not by dispersion alone

Non-covalent interactions between neutral, sterically hindered organic molecules generally involve a strong stabilizing contribution from dispersion forces that in many systems turns the 'steric repulsion' into a 'steric attraction'. In addition to London dispersion, such systems benefit from electrostatic stabilization, which arises from a short-range effect of charge penetration and gets bigger with increasing steric bulk. In the present work, we quantify this contribution for a diverse set of molecular cores, ranging from unsubstituted benzene and cyclohexane to their derivatives carrying tert-butyl, phenyl, cyclohexyl and adamantyl substituents. While the importance of electrostatic interactions in the dimers of sp2-rich (e.g., π-conjugated) cores is well appreciated, less polarizable assemblies of sp3-rich systems with multiple short-range CH···HC contacts between the bulky cyclohexyl and adamantyl moieties are also significantly influenced by electrostatics. Charge penetration is drastically larger in absolute terms for the sp2-rich cores, but still has a non-negligible effect on the sp3-rich dimers, investigated herein, both in terms of their energetics and equilibrium interaction distances. These results emphasize the importance of this electrostatic effect, which has so far been less recognized in aliphatic systems compared to London dispersion, and are therefore likely to have implications for the development of force fields and methods for crystal structure prediction.

First-principle simulations of electronic structure in semicrystalline polyethylene

In order to increase our fundamental knowledge about high-voltage cable insulation materials, realistic polyethylene (PE) structures, generated with a novel molecular modeling strategy, have been analyzed using first principle electronic structure simulations. The PE structures were constructed by first generating atomistic PE configurations with an off-lattice Monte Carlo method and then equilibrating the structures at the desired temperature and pressure using molecular dynamics simulations. Semicrystalline, fully crystalline and fully amorphous PE, in some cases including crosslinks and short-chain branches, were analyzed. The modeled PE had a structure in agreement with established experimental data. Linear-scaling density functional theory (LS-DFT) was used to examine the electronic structure (e.g., spatial distribution of molecular orbitals, bandgaps and mobility edges) on all the materials, whereas conventional DFT was used to validate the LS-DFT results on small systems. When hybrid functionals were used, the simulated bandgaps were close to the experimental values. The localization of valence and conduction band states was demonstrated. The localized states in the conduction band were primarily found in the free volume (result of gauche conformations) present in the amorphous regions. For branched and crosslinked structures, the localized electronic states closest to the valence band edge were positioned at branches and crosslinks, respectively. At 0 K, the activation energy for transport was lower for holes than for electrons. However, at room temperature, the effective activation energy was very low (∼0.1 eV) for both holes and electrons, which indicates that the mobility will be relatively high even below the mobility edges and suggests that charge carriers can be hot carriers above the mobility edges in the presence of a high electrical field.