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

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.

Wetting Behavior of a Three-Phase System in Contact with a Surface

We extend the Cahn-Landau-de Gennes mean field theory of wetting in binary mixtures to understand the wetting thermodynamics of a three phase system (e.g., polymer dispersed liquid crystals or polymer-colloid mixtures) that is in contact with an external surface, which prefers one of the phases. Using a model free-energy, which has three minima in its landscape, we show that as the central minimum becomes more stable compared to the remaining ones, the bulk phase diagram encounters a triple point and then bifurcates and we observe a novel non-monotonic dependence of the surface tension as a function of the stability of the central minimum. We show that this non-monotonicity in surface tension is associated with a complete to partial wetting transition. We obtain the complete wetting phase behavior as a function of phase stability and the surface interaction parameters when the system is close to the bulk triple point. The model free-energy that we use is qualitatively similar to that of a renormalized free energy, which arises in the context of polymer-liquid crystal mixtures. Finally, we study the thermodynamics of wetting for an explicit polymer-liquid crystal mixture and show that its thermodynamics is similar to that of our model free-energy.

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.

Semiflexible oligomers crystallize via a cooperative phase transition

Semicrystalline polymers are ubiquitous, yet despite their fundamental and industrial importance, the theory of homogeneous nucleation from a melt remains a subject of debate. A key component of the controversy is that polymer crystallization is a non-equilibrium process, making it difficult to distinguish between effects that are purely kinetic and those that arise from the underlying thermodynamics. Due to computational cost constraints, simulations of polymer crystallization typically employ non-equilibrium molecular dynamics techniques with large degrees of undercooling that further exacerbate the coupling between thermodynamics and kinetics. In a departure from this approach, in this study, we isolate the near-equilibrium nucleation behavior of a simple model of a melt of short, semiflexible oligomers. We employ several Monte Carlo methods and compute a phase diagram in the temperature-density plane along with two-dimensional free energy landscapes (FELs) that characterize the nucleation behavior. The phase diagram shows the existence of ordered nematic and crystalline phases in addition to the disordered melt phase. The minimum free energy path in the FEL for the melt-crystal transition shows a cooperative transition, where nematic order and monomer positional order move in tandem as the system crystallizes. This near-equilibrium phase transition mechanism broadly agrees with recent evidence that polymer stiffness plays an important role in crystallization but differs in the specifics of the mechanism from several recent theories. We conclude that the computation of multidimensional FELs for models that are larger and more fine-grained will be important for evaluating and refining theories of homogeneous nucleation for polymer crystallization.

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.

The Sequence of a Step-Growth Copolymer Can Be Influenced by Its Own Persistence Length

Synthetic copolymer sequences remain challenging to control, and there are features of even simple one-pot, solution-based copolymerizations that are not yet fully understood. In previous simulations on step-growth copolymerizations in solution, we demonstrated that modest variations in the attractions between type A and B monomers could significantly influence copolymer sequence through an emergent aggregation and phase separation initiated by the lengthening of nascent oligomers. Here we investigate how one aspect of a copolymer's geometry-its flexibility-can modulate those effects. Our simulations show the onset of strand alignment within the polymerization-induced aggregates as chain stiffness increases and demonstrate that this alignment can influence the resulting copolymer sequences. For less flexible copolymers, with persistence lengths ≥10 monomers, modest nonbonded attractions of ∼k_BT between monomers of the same type yield A and B blocks of a characteristic length and result in a polydispersity index that grows rapidly, peaks, and then diminishes as the reaction proceeds. These results demonstrate that for copolymer systems with modest variations in intermonomer attractions and physically realistic flexibilities a nascent copolymer's persistence length can influence its own sequence.

Direct observation of long chain enrichment in flow-induced nuclei from molecular dynamics simulations of bimodal blends

Modelling of flow-induced nucleation in polymers suggest that long chains are enriched in nuclei, relative to their melt concentration. This enrichment has important consequences for the nucleation rate and mechanism, but cannot be directly observed with current experimental techniques. Instead, we ran united atom molecular dynamics simulations of bimodal polyethylene blends, comprising linear chains at a 50 : 50 mix of long (1000 carbon) and short (500-125 carbon) chains, under shear flow. We developed a method to extract the nucleus composition during a transient start-up flow. Our simulations show significant and systematic enrichment of long-chains for all nucleus sizes up to and beyond the critical nucleus. This enrichment is quantitatively predicted by the recent polySTRAND model [Read et al. Phys. Rev. Lett. 2020, 124, 147802]. The same model parameters also correctly capture the nucleus induction time in our simulations. All parameters of the model were fitted to a small subset of our data in which long chain enhancement was absent. We conclude that long-chain enrichment is central to the mechanism of flow-induced nucleation and that this enrichment must be captured to correctly predict the nucleation rate.

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.

Monodisperse Polymer Melts Crystallize via Structurally Polydisperse Nanoscale Clusters: Insights from Polyethylene

This study demonstrates that monodisperse entangled polymer melts crystallize via the formation of nanoscale nascent polymer crystals (i.e., nuclei) that exhibit substantial variability in terms of their constituent crystalline polymer chain segments (stems). More specifically, large-scale coarse-grain molecular simulations are used to quantify the evolution of stem length distributions and their properties during the formation of polymer nuclei in supercooled prototypical polyethylene melts. Stems can adopt a range of lengths within an individual nucleus (e.g., ∼1-10 nm) while two nuclei of comparable size can have markedly different stem distributions. As such, the attainment of chemically monodisperse polymer specimens is not sufficient to achieve physical uniformity and consistency. Furthermore, stem length distributions and their evolution indicate that polymer crystal nucleation (i.e., the initial emergence of a nascent crystal) is phenomenologically distinct from crystal growth. These results highlight that the tailoring of polymeric materials requires strategies for controlling polymer crystal nucleation and growth at the nanoscale.