Molecular Dynamics Simulations for the Description of Experimental Molecular Conformation, Melt Dynamics, and Phase Transitions in Polyethylene

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.

Preparation and Characterization of High-Density Polyethylene with Alternating Lamellar Stems Using Molecular Dynamics Simulations

Mechanical recycling is the most efficient way to reduce plastic pollution due to its ability to maintain the intrinsic properties of plastics as well as provide economic benefits involved in other types of recycling. On the other hand, molecular dynamics (MD) simulations provide key insights into structural deformation, lamellar crystalline axis (c-axis) orientations, and reorganization, which are essential for understanding plastic behavior during structural deformations. To simulate the influence of structural deformations in high-density polyethylene (HDPE) during mechanical recycling while paying attention to obtaining an alternate lamellar orientation, the authors examine a specific way of preparing stacked lamella-oriented HDPE united atom (UA) models, starting from a single 1000 UA (C1000) chain of crystalline conformations and then packing such chain conformations into 2-chain, 10-chain, 15-chain, and 20-chain semi-crystalline models. The 2-chain, 10-chain, and 15-chain models yielded HDPE microstructures with the desired alternating lamellar orientations and entangled amorphous segments. On the other hand, the 20-chain model displayed multi-nucleus crystal growth instead of the lamellar-stack orientation. Structural characterization using a one-dimensional density profile and local order parameter {P2(r)} analyses demonstrated lamellar-stack orientation formation. All semi-crystalline models displayed the total density (ρ) and degree of crystallinity (χ) range of 0.90-0.94 g/cm-3 and ≥42-45%, respectively. A notable stress yield (σ_yield) ≈ 100-120 MPa and a superior elongation at break (ε_break) ~250% was observed under uniaxial strain deformation along the lamellar-stack orientation. Similarly, during the MD simulations, the microstructure phase change represented the average number of entanglements per chain (). From the present study, it can be recommended that the 10-chain alternate lamellar-stack orientation model is the most reliable miniature model for HDPE that can mimic industrially relevant plastic behavior in various conditions.

Unraveling Tribochemistry and Self-Lubrication Mechanism of Polytetrafluoroethylene by Reactive Coarse-Grained Molecular Dynamics Simulations

Lubrication of polymeric materials generally involves processes of atomic-scale chemical bond forming/breaking at the interface and mesoscale chain reorientation, disentanglement, and so forth. However, it is difficult to describe the important aspects of tribochemical reactions by conventional coarse-grained molecular dynamics (CGMD) simulations. Here, reactive CGMD simulations were conducted based on the ReaxFF force field to study the tribochemical interactions between polytetrafluoroethylene (PTFE) and iron. The chemical bond forming/breaking between the molecular chain and countersurface was fitted through the bond dissociation energies of specific reaction sites from all-atom ReaxFF-MD simulations. This enabled a quantitative description of tribochemical reactions in a macromolecule system. First, the number of anchoring bonds between PTFE molecules and the countersurface showed a strong correlation with the friction coefficient. The shearing process induced breaking of the interfacial anchoring bonds as well as chain disentanglement in the matrix, which consequently led to ordering reorientation of molecular chains toward sliding direction and hence decrease of friction. Second, two competitive factors were clarified to affect polymer friction with varying temperatures. The decrease of interfacial anchoring reactivity and molecular chain mobility at low temperature prohibited reorientation of molecular chains and increased the friction coefficient. On the other hand, the hardening of PTFE and the reduction in effective contact area at low temperatures decreased the friction coefficient. This led to a turning point with a maximum friction coefficient around 100 K. These results shed light on the essential role of tribochemical reactions on polymer lubrication, especially under low temperatures, which provides design guidance of polymeric lubrication systems for engineering applications.

Cooling-Rate Computer Simulations for the Description of Crystallization of Organic Phase-Change Materials

A molecular-level insight into phase transformations is in great demand for many molecular systems. It can be gained through computer simulations in which cooling is applied to a system at a constant rate. However, the impact of the cooling rate on the crystallization process is largely unknown. To this end, here we performed atomic-scale molecular dynamics simulations of organic phase-change materials (paraffins), in which the cooling rate was varied over four orders of magnitude. Our computational results clearly show that a certain threshold (1.2 × 10¹¹ K/min) in the values of cooling rates exists. When cooling is slower than the threshold, the simulations qualitatively reproduce an experimentally observed abrupt change in the temperature dependence of the density, enthalpy, and thermal conductivity of paraffins upon crystallization. Beyond this threshold, when cooling is too fast, the paraffin's properties in simulations start to deviate considerably from experimental data: the faster the cooling, the larger part of the system is trapped in the supercooled liquid state. Thus, a proper choice of a cooling rate is of tremendous importance in computer simulations of organic phase-change materials, which are of great promise for use in domestic heat storage devices.

A combined experimental and molecular simulation study on stress generation phenomena during the Ziegler-Natta polyethylene catalyst fragmentation process

The morphology of particles obtained under different pre-polymerization conditions has been connected to the stress generation mechanism at the polymer/catalyst interface. A combination of experimental characterization techniques and atomistic molecular dynamics simulations allowed a systematic investigation of experimental conditions leading to a certain particle morphology, and hence to a final polymer with specific features. Atomistic models of nascent polymer phases in contact with magnesium dichloride surfaces have been developed and validated. Using these detailed models, in the framework of McKenna's hypothesis, the pressure increase due to the polymerization reaction has been calculated under different conditions and is in good agreement with experimental scenarios. This molecular scale knowledge and the proposed investigation strategy would allow the pre-polymerization conditions to be better defined and the properties of the nascent polymer to be tuned, ensuring proper operability along the whole polymer production process.

Interfacial Forces in Free-Standing Layers of Melted Polyethylene, from Critical to Nanoscopic Thicknesses

Molecular dynamics simulations of ultrathin free-standing layers made of melted (373.15–673.15 K) polyethylene chains, which exhibit a lower melting temperature (compared to the bulk value), were carried out to investigate the dominant pressure forces that shape the conformation of chains at the interfacial and bulk liquid regions. We investigated layer thicknesses, tL, from the critical limit of mechanical stability up to lengths of tens of nm and found a normal distribution of bonds dominated by slightly stretched chains across the entire layer, even at large temperatures. In the bulk region, the contribution of bond vibrations to pressure was one order of magnitude larger than the contributions from interchain interactions, which changed from cohesive to noncohesive at larger temperatures just at a transition temperature that was found to be close to the experimentally derived onset temperature for thermal stability. The interchain interactions produced noncohesive interfacial regions at all temperatures in both directions (normal and lateral to the surface layer). Predictions for the value of the surface tension, γ, were consistent with experimental results and were independent of tL. However, the real interfacial thickness—measured from the outermost part of the interface up to the point where γ reached its maximum value—was found to be dependent on tL, located at a distance of 62 Å from the Gibbs dividing surface in the largest layer studied (1568 chains or 313,600 bins); this was ~4 times the length of the interfacial thickness measured in the density profiles.

Comparative experimental and modelling study of the thermal and thermo-mechanical properties of LLDPE/wax blends

The interactions and morphology of molecules in a polymer blend influence the physical properties of the blend. However, little is known about the influence of molecular interaction and morphology on the thermal and mechanical properties of LLDPE/wax blends. Although cooling rate can be used to investigate blends' thermal and mechanical properties, it is inadequate to determine interactions between the molecules in the LLDPE/wax blends. However, since the morphology is related to the thermal and mechanical properties of polymer blends and could be related to the cooling rate, LLDPE/wax samples prepared by melt mixing were cooled at different rates. The thermal and mechanical properties of the LLDPE/wax blends were modelled through molecular dynamic simulations. The modelled transitions were compared to experimentally determined mechanical relaxations of LLDPE/wax blends to investigate the effect of wax addition on the blend crystallinity. The crystallization behaviour of the blends was studied by differential scanning calorimetry, dynamic mechanical behaviour by dynamic mechanical analysis, and differences in crystallinity by X-ray diffraction. There were no significant differences between the results for the slow- and quench-cooled samples, confirming the rapid crystallization of both the LLDPE and the wax. Experiments and molecular dynamics simulations confirmed the cocrystallization of wax with LLDPE.

Chiral selecting crystallization of helical polymers: A molecular dynamics simulation for the POM-like bare helix

Polymer crystallization has long been a fascinating problem and is still attracting many researchers. Most of the previous simulations are concentrated on clarifying the universal aspects of polymer crystallization using model linear polymers such as polyethylene. We are recently focusing on a nearly untouched but very interesting problem of chiral selecting crystallization in helical polymers. We previously proposed a stepwise approach using two kinds of helical polymers, simple "bare" helical polymers made of backbone atoms only such as polyoxymethylene (POM) and "general" helical polymers containing complicated side groups such as isotactic polypropylene (iPP). We have already reported on the crystallization in oligomeric POM-like helix but have observed only weak chiral selectivity during crystallization. In the present paper, we investigate the crystallization of sufficiently long POM-like polymer both from the isotropic melt and from the highly stretched melt. We find in both cases that the polymer shows a clear chiral selecting crystallization. Especially the observation of a single crystal growing from the isotropic melt is very illuminating. It shows that the crystal thickness and the crystal chirality is closely correlated; thicker crystals show definite chirality while thinner ones are mostly mixtures of the R- and the L- handed stems. The single crystal is found to have a marked lenticular shape, where the thinner growth front, since being made of the mixture, shows no chiral selectivity. Final chiral crystal is found to be completed through helix reversal processes within thicker regions.

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.

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.

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.

Investigation of the properties of nanographene in polymer nanocomposites through molecular simulations: dynamics and anisotropic Brownian motion

The dynamical behavior of nanographene sheets dispersed in polymer matrices is investigated through united-atom molecular dynamics simulations. The Brownian motion of the sheet and the anisotropy in its translational and orientational diffusion are the topics of the current study. Different polymer matrices and pristine and functionalized graphene constitute various nanocomposite systems. Interactions between the nanographene flake and the matrix determine the dynamics of the systems. The dynamics is reduced in polyethylene oxide compared to polyethylene matrix, whereas carboxylated sheets move considerably slower than the pristine nanographene in any matrix. Diffusion is anisotropic for short times, while it becomes isotropic in the long time limit. The in-plane motion of the nanographene sheet is faster than the out-of-plane component, in agreement with the diffusion of perfectly oblate ellipsoids. In functionalized graphene, the anisotropy is suppressed. By exploring the temperature effect on both the nanographene sheet and polymer close to the surface, indications for coupling in the motion of the two components are revealed. The strong effect of edge functional groups on the dynamics can be used as a way to control the Brownian motion of nanographene sheets in polymer nanocomposites and consequently tailor the properties of the materials.

Uncertainty Quantification for Mechanical Properties of Polyethylene Based on Fully Atomistic Model

This study is to assess the effect of temperature and strain rate on the mechanical properties of amorphous polyethylene (PE) based on fully atomistic model. A stochastic constitutive model using data obtained from molecular dynamics (MD) simulations for the material is constructed. Subsequently, a global sensitivity analysis approach is then employed to predict the essential parameters of the mechanical model. The sensitivity indices show that the key parameter affecting Young’s modulus and yield stress is the temperature followed by the strain rate.

Disentangling and Lamellar Thickening of Linear Polymers during Crystallization: Simulation of Bimodal and Unimodal Molecular Weight Distribution Systems

We have performed coarse-grained molecular dynamics simulations to study the isothermal crystallization of bimodal and unimodal molecular weight distribution (MWD) polymers with equivalent average molecular weight (M_w). By using primitive path analysis, we can monitor the entanglement evolution during the process of crystallization. We have discovered a quantitative correlation between the degree of disentanglement and crystallinity, indicating that chain disentanglement permits the process of crystallization. In addition, the crystalline stem length also displays a linear relation with the degree of disentanglement at different temperatures. Based on the observation in our simulations, we can build a scenario of the whole process of chain disentangling and lamellar thickening on the basis of chain sliding diffusion. Furthermore, we have enough evidence to infer that the temperature dependence of crystalline stem length is basically a result of temperature dependence of chain sliding diffusion. Our observations are also in agreement with Hikosaka's sliding diffusion theory. Compared to the unimodal system, the disentanglement degree of the bimodal system is more delayed than its crystallinity due to the slower chain sliding of the long-chain component; the bimodal system reaches a larger crystalline stem length at all temperatures due to the promotion of higher chain sliding mobility of the short-chain component.

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.

A coarse-grain model for entangled polyethylene melts and polyethylene crystallization

The Shinoda-DeVane-Klein (SDK) model is herein demonstrated to be a viable coarse-grain model for performing molecular simulations of polyethylene (PE), affording new opportunities to advance molecular-level, scientific understanding of PE materials and processes. Both structural and dynamical properties of entangled PE melts are captured by the SDK model, which also recovers important aspects of PE crystallization phenomenology. Importantly, the SDK model can be used to represent a variety of materials beyond PE and has a simple functional form, making it unique among coarse-grain PE models. This study expands the suite of tools for studying PE in silico and paves the way for future work probing PE and PE-based composites at the molecular level.

Universal localization transition accompanying glass formation: insights from efficient molecular dynamics simulations of diverse supercooled liquids

The origin of the precipitous dynamic arrest known as the glass transition is a grand open question of soft condensed matter physics. It has long been suspected that this transition is driven by an onset of particle localization and associated emergence of a glassy modulus. However, progress towards an accepted understanding of glass formation has been impeded by an inability to obtain data sufficient in chemical diversity, relaxation timescales, and spatial and temporal resolution to validate or falsify proposed theories for its physics. Here we first describe a strategy enabling facile high-throughput simulation of glass-forming liquids to nearly unprecedented relaxation times. We then perform simulations of 51 glass-forming liquids, spanning polymers, small organic molecules, inorganics, and metallic glass-formers, with longest relaxation times exceeding one microsecond. Results identify a universal particle-localization transition accompanying glass formation across all classes of glass-forming liquid. The onset temperature of non-Arrhenius dynamics is found to serve as a normalizing condition leading to a master collapse of localization data. This transition exhibits a non-universal relationship with dynamic arrest, suggesting that the nonuniversality of supercooled liquid dynamics enters via the dependence of relaxation times on local cage scale. These results suggest that a universal particle-localization transition may underpin the glass transition, and they emphasize the potential for recent theoretical developments connecting relaxation to localization and emergent elasticity to finally explain the origin of this phenomenon. More broadly, the capacity for high-throughput prediction of glass formation behavior may open the door to computational inverse design of glass-forming materials.

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.

Multiscale modeling of charge transfer in polymers with flexible backbones

In order to evaluate carrier transfer properties in polymers with flexible backbones, we have proposed a simplified multi-scale modeling approach combining molecular dynamics simulations, first-principles calculations and kinetic Monte Carlo simulations.

Atomistic modelling of entropy driven phase transitions between different crystal modifications in polymers: the case of poly(3-alkylthiophenes)

Polymorphism and related solid-state phase transitions affect the structure and morphology and hence the properties of materials, but they are not-so-well understood. Atomistic computational methods can provide molecular-level insights, but they have rarely proven successful for transitions between polymorphic forms of crystalline polymers. In this work, we report atomistic molecular dynamics (MD) simulations of poly(3-alkylthiophenes) (P3ATs), widely used organic semiconductors to explore the experimentally observed, entropy-driven transition from form II to more common form I type polymorphs, or, more precisely, to form I mesophases. The transition is followed continuously, also considering X-ray diffraction evidence, for poly(3-hexylthiophene) (P3HT) and poly(3-butylthiophene) (P3BT), evidencing three main steps: (i) loss of side chain interdigitation, (ii) partial disruption of the original stacking order and (iii) reorganization of polymer chains into new, tighter, main-chain stacks and new layers with characteristic form I periodicities, substantially larger than those in the original form II. The described approach, likely applicable to other important transitions in polymers, provides previously inaccessible insight into the structural organization and disorder features of form I structures of P3ATs, not only in their development from form II structures but also from melts or solutions.

Coarse-graining simulation approaches for polymer melts: the effect of potential range on computational efficiency

The integral equation coarse-graining (IECG) approach is a promising high-level coarse-graining (CG) method for polymer melts, with variable resolution from soft spheres to multi CG sites, which preserves the structural and thermodynamical consistencies with the related atomistic simulations. When compared to the atomistic description, the procedure of coarse-graining results in smoother free energy surfaces, longer-ranged potentials, a decrease in the number of interaction sites for a given polymer, and more. Because these changes have competing effects on the computational efficiency of the CG model, care needs to be taken when studying the effect of coarse-graining on the computational speed-up in CG molecular dynamics simulations. For instance, treatment of long-range CG interactions requires the selection of cutoff distances that include the attractive part of the effective CG potential and force. In particular, we show how the complex nature of the range and curvature of the effective CG potential, the selection of a suitable CG timestep, the choice of the cutoff distance, the molecular dynamics algorithms, and the smoothness of the CG free energy surface affect the efficiency of IECG simulations. By direct comparison with the atomistic simulations of relatively short chain polymer melts, we find that the overall computational efficiency is highest for the highest level of CG (soft spheres), with an overall improvement of the computational efficiency being about 106-108 for various CG levels/resolutions. Therefore, the IECG method can have important applications in molecular dynamics simulations of polymeric systems. Finally, making use of the standard spatial decomposition algorithm, the parallel scalability of the IECG simulations for various levels of CG is presented. Optimal parallel scaling is observed for a reasonably large number of processors. Although this study is performed using the IECG approach, its results on the relation between the level of CG and the computational efficiency are general and apply to any properly-constructed CG model.

Crystal Growth in Polyethylene by Molecular Dynamics: The Crystal Edge and Lamellar Thickness

We carried out large-scale atomistic molecular dynamics simulations to study the growth of twin lamellar crystals of polyethylene initiated by small crystal seeds. By examining the size distribution of the stems-straight crystalline polymer segments-we show that the crystal edge has a parabolic profile. At the growth front, there is a layer of stems too short to be stable, and new stable stems are formed within this layer, leading to crystal growth. Away from the edge, the lengthening of the stems is limited by a lack of available slack length in the chains. This frustration can be relieved by mobile crystal defects that allow topological relaxation by traversing through the crystal. The results shed light on the process of polymer crystal growth and help explain initial thickness selection and lamellar thickening.

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.