Simulation study of semi-crystalline polymer interphases

Molecular Simulations of Hydrogen Sorption in Semicrystalline High-Density Polyethylene: The Impact of the Surface Fraction of Tie-Chains

The modeling of the barrier properties of semicrystalline polymers has gained interest following the possible application of such materials as protective liners for the safe supply of pressurized hydrogen. The mass transport in such systems is intimately related to the complex intercalation between the crystal and amorphous phases, which was approached in this work through an all-atom representation of high-density polyethylene structures with a tailored fraction of amorphous-crystalline connections (tie-chains). Simulations of the polymer pressure-volume-temperature data and hydrogen sorption were performed by means of molecular dynamics and the Widom test particle insertion method. The discretization of the simulation domains of the semicrystalline structures allowed us to obtain profiles of density, degree of order, and gas solubility. The results indicated that the gas sorption in the crystalline regions is negligible and that the confinement of the amorphous phase between crystals induces a significant increase in density and a drop in the sorption capacity, even in the absence of tie-chains. Adding ties between the crystal and the amorphous phase results in further densification, an increase of the lamella tilt angle, and a decrease in the degree of crystallinity and hydrogen sorption coefficient, in agreement with several literature references.

Quantitative Model of Multiple Crystal Growth Rate Minima in Polymers with Regularly Spaced Substituent Groups

A simple theory has been developed to explain quantitatively the multiple crystal growth rate minima observed experimentally in polyethylene brassylates (PEBs), polymers with regularly spaced "chemical defects", in this case, diester groups separated by 11 methylenes. The minima occur at the transitions where the fold length drops from 4 to 3 repeat units and from 3 to 2 units. An analytical rate-equation model was developed with elementary attachment and detachment steps of individual monomer repeat units, also including postattachment stem lengthening (stem conversion). The model produced a good fit to experimental crystallization rate curves for PEBs of three different molecular weights. The fits confirm in a quantitative way that the anomalies are caused by the self-poisoning effect, as proposed in the original experimental report on PEBs, based on the ideas developed in previous studies on long-chain n-alkanes. It is concluded that the rate minima in PEBs are the result of temporary attachment to the growth surface of stems that are too short to be stable yet long enough and close to stability to obstruct productive growth by stems of sufficient length. The results confirm the ubiquitous presence of self-poisoning at the growth front of polymer crystals in general and will help to achieve a better understanding of the complex process of crystallization of polymers. It will also allow the determination of more realistic parameters controlling their lamellar growth kinetics.

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.

Quasi-continuous melting of model polymer monolayers prompts reinterpretation of polymer melting

Condensed matter textbooks teach us that melting cannot be continuous and indeed experience, including with polymers and other long-chain compounds, tells us that it is a strongly first-order transition. However, here we report nearly continuous melting of monolayers of ultralong n-alkane C₃₉₀H₇₈₂ on graphite, observed by AFM and reproduced by mean-field theory and MD simulation. On heating, the crystal-melt interface moves steadily and reversibly from chain ends inward. Remarkably, the final melting point is 80 K above that of the bulk, and equilibrium crystallinity decreases continuously from ~100% to <50% prior to final melting. We show that the similarity in melting behavior of polymers and non-polymers is coincidental. In the bulk, the intermediate melting stages of long-chain crystals are forbidden by steric overcrowding at the crystal-liquid interface. However, there is no crowding in a monolayer as chain segments can escape to the third dimension.

A Local-Exchange Model of Folding Chain Surface of Polymer Crystal Based on Worm-Like Chain Model within Single-Chain in Mean-Field Theory

The structure of amorphous layer of folding surface controls the properties of the polymer lamellar crystal, which consists of chains with a loop conformation. The surface tension depends on the length and the distance between two injection points of the loop which involving the reptation motion and lateral exchange motion of the stems. In the present work, a local-exchange motion model based on the worm-like chain model is developed to investigate the effects of lateral motion of stems on the release the surface tension. The optimal distance between two injection points is determined by the balance of chain bending energy and conformational entropy. The numerical results provide evidences to the adjacent re-entry model for various loop lengths. A possible explanation involving density of injection points is proposed to interpret the mechanism.

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.

Self-consistent field model simulations for statistics of amorphous polymer chains in crystalline lamellar structures

We calculate statistical properties of amorphous polymer chains between crystalline lamellae by self-consistent field model simulations. In our model, an amorphous subchain is modelled as a polymer chain of which ends are grafted onto the crystal-amorphous interfaces. The crystal-amorphous interfaces are expressed as impenetrable surfaces. We incorporate the interaction between segments to satisfy the incompressible condition for the segment density field. The simulation results show that amorphous polymer chains feel thin potential layers, which are mainly repulsive, near the crystal-amorphous interfaces. The impenetrable and incompressible conditions affect the statistics of polymer chains and the chain statistics becomes qualitatively different from the ideal Gaussian chain statistics without any constraints. We show the effects of the system size and the graft density to statistical quantities. We also show that the tie subchain statistics obey rather simple statistics.