Free Surfaces Overcome Superheating in Simulated Melting of Isotactic Polypropylene

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.

Melting of superheated molecular crystals

Melting dynamics of micrometer scale, polycrystalline samples of isobutane, dimethyl ether, methyl benzene, and 2-propanol were investigated by fast scanning calorimetry. When films are superheated with rates in excess of 10⁵ K s^-1, the melting process follows zero-order, Arrhenius-like kinetics until approximately half of the sample has transformed. Such kinetics strongly imply that melting progresses into the bulk via a rapidly moving solid-liquid interface that is likely to originate at the sample's surface. Remarkably, the apparent activation energies for the phase transformation are large; all exceed the enthalpy of vaporization of each compound and some exceed it by an order of magnitude. In fact, we find that the crystalline melting kinetics are comparable to the kinetics of dielectric α-relaxation in deeply supercooled liquids. Based on these observations, we conclude that the rate of non-isothermal melting for superheated, low-molecular-weight crystals is limited by constituent diffusion into an abnormally dense, glass-like, non-crystalline phase.

Using surface-induced ordering to probe the isotropic-to-nematic transition for semiflexible polymers

Semiflexible polymers undergo a weakly first order isotropic-to-nematic (IN) phase transition when the volume fraction ϕ is high enough that random alignment of the backbone segments is no longer viable. For semiflexible chains, the critical volume fraction ϕc is governed by the backbone stiffness Np. To locate the IN phase transition, we perform molecular dynamics (MD) simulations of bead-spring chains confined between two impenetrable parallel surfaces. We use the impenetrable surfaces to induce nematic-isotropic interfaces for semiflexible chains in the isotropic phase. By progressively increasing the backbone stiffness Np, we observe the propagation of surface-induced nematic order above a critical stiffness N for a given ϕ. Using the simulation results N(ϕ), we construct the IN phase boundry in the ϕ-Np plane, from which the scaling relation between ϕc and Np is obtained. For semiflexible chains with Np ≤ 5.78, our results suggest ϕc ∼ Np(-1), consistent with prediction by Khokhlov and Semenov. For chains with Np ≥ 5.78, we observe a new scaling regime in which ϕc ∼ Np(-2/3).