Langevin dynamics simulations of early-stage polymer nucleation and crystallization

Crystallization with Nodular Aggregation near the Glass Transition Temperature for Syndiotactic Polypropylene

The isothermal crystallization from the melt state of syndiotactic polypropylene (sPP) has been studied by wide-angle X-ray diffraction (WAXD), small-angle X-ray scattering (SAXS), and optical microscopy. The WAXD and SAXS results show the crystallization mechanism near the glass transition temperature in which the crystalline and mesomorphic nodules cover the entire sample with the formation of aggregation regions. For the SAXS analysis, the scattering function for the three-component system has been suggested. Furthermore, to analyze the growth kinetics of the aggregation region for sPP, the time-dependent structure factor combined with the homogeneous and inhomogeneous nucleation-and-growth kinetics has been suggested. The analysis shows that the growth kinetics of the aggregation region for sPP is the homogeneous nucleation-and-growth. The growth velocity of the aggregation region is a natural extrapolation of that of spherulite to the high supercooling region. These results might indicate that the crystallization with the nodular aggregation is a fundamental crystallization process near the glass transition temperature for polymers.

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 Chains Fold Prior to Crystallization

There are long-standing debates in crystallization mechanism of polymer chains at the molecular levels: Which comes first, chain folding or lamellae formation during crystallization? In this study, we report the local chain trajectory of ¹³C-labeled semicrystalline polymer in an extreme case of rapidly quenched glassy state as well as thermodynamically stable crystals formed via different pathways from glass and melt. Magnetically dipole interactions do not require a long-range order of molecular objects and thus enable us to trace the local chain trajectory of polymer chains even in a glassy state. To accurately characterize the local chain trajectory of polymer glass, the natural abundance effect on ¹³C-¹³C double-quantum (DQ) nuclear magnetic resonance (NMR) signal is re-examined using extended chain conformation. As results, it is found that glassy chains adopt the same adjacent re-entry structure (adjacent re-entry number, n = 1) with the melt- and cold-grown crystals. From these results, it is concluded that (i) folding occurs prior to crystallization and (ii) melt and cold crystallization do not induce additional folding but proceed with rearrangements of polymer chains in the existing templates.

Kinetics of Polymer Crystallization with Aggregating Small Crystallites

The isothermal crystallization near the glass transition temperature from the melt state of poly(trimethylene terephthalate) has been studied by wide-angle x-ray diffraction (WAXD), small-angle x-ray scattering (SAXS), and optical microscopy. The SAXS and WAXD results show the crystallization mechanism in which the crystalline nodules cover the entire sample with the formation of aggregation regions. The analysis of the SAXS results using Kolmogorov-Johnson-Mehl-Avrami theory indicates that the formation kinetics of the aggregation regions is of three-dimensional homogeneous nucleation type. The analysis of the SAXS profiles using Sekimoto's theory provides the growth velocity and the nucleation rate of the aggregation region. The temperature dependence of the growth velocity of the aggregation region is a natural extrapolation of that of spherulite to the high supercooling region. The temperature dependence of the nucleation rate of the aggregation region is also represented by the parameters of the spherulitic growth rate. The result of the growth velocities of the aggregation region and the spherulite suggests the existence of precursors at the front of the crystal growth.

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.

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.

Crystallization of semiflexible polymers in melts and solutions

We studied the crystallization of semiflexible polymer chains in melts and poor-solvent solutions with different concentrations using dissipative particle dynamics (DPD) computer simulation techniques. We used the coarse-grained polymer model to reveal the general principles and microscopic scenario of crystallization in such systems at large time and length scales. It covers both primary and secondary nucleation as well as crystallites' merging. The parameters of the DPD model were chosen appropriately to reproduce the entanglements of polymer chains. We started from an initial homogeneous disordered solution of Gaussian chains and observed the initial stages of crystallization process caused in our model by orientational ordering of polymer chains and polymer-solvent phase separation. We found that the overall crystalline fraction at the end of the crystallization process decreases with the increasing polymer volume fraction while the steady-state crystallization speed at later stages does not depend on the polymer volume fraction. The average crystallite size has a maximal value in the systems with a polymer volume fraction from 0.7 to 0.95. In our model, these polymer concentrations represent an optimal value in the sense of balance between the amount of polymer material available to increase the crystallite size and chain entanglements, that prevent crystallites' growth and merging. On large time scales, our model allows us to observe lamellar thickening linear in logarithmic time scale.

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.

Divining the shape of nascent polymer crystal nuclei

We demonstrate that nascent polymer crystals (i.e., nuclei) are anisotropic entities with neither spherical nor cylindrical geometry, in contrast to previous assumptions. In fact, cylindrical, spherical, and other high symmetry geometries are thermodynamically unfavorable. Moreover, postcritical transitions are necessary to achieve the lamellae that ultimately arise during the crystallization of semicrystalline polymers. We also highlight how inaccurate treatments of polymer nucleation can lead to substantial errors (e.g., orders of magnitude discrepancies in predicted nucleation rates). These insights are based on quantitative analysis of over four million crystal clusters from the crystallization of prototypical entangled polyethylene melts. New comprehensive bottom-up models are needed to capture polymer nucleation.

Polymer nucleation under high-driving force, long-chain conditions: Heat release and the separation of time scales

This study reveals important features of polymer crystal formation at high-driving forces in entangled polymer melts based on simulations of polyethylene. First and in contrast to small-molecule crystallization, the heat released during polymer crystallization does not appreciably influence structural details of early-stage, crystalline clusters (crystal nuclei). Second, early-stage polymer crystallization (crystal nucleation) can occur without substantial chain-level relaxation and conformational changes. This study's results indicate that local structures and environments guide crystal nucleation in entangled polymer melts under high-driving force conditions. Given that such conditions are often used to process polyethylene, local structures and the separation of time scales associated with crystallization and chain-level processes are anticipated to be of substantial importance to processing strategies. This study highlights new research directions for understanding polymer crystallization.

An accelerated united-atom molecular dynamics simulation on the fast crystallization of ring polyethylene melts

To investigate crystallinities based on trans-structures, we determined the differences in the crystallization properties of ring and linear polymers by performing united-atom-model molecular dynamics (MD) simulations of homogeneous polyethylene melts of equal length, N, which refers to the number of monomers per chain. Modified parameters based on the DREIDING force field for the CH₂ units were used in order to accelerate the crystallization process. To detect polymer crystallization, we introduced some local-order parameters that relate to trans-segments in addition to common crystallinities using neighboring bond orders. Through quenching MD simulations at 5 K/ns, we roughly determined temperature thresholds, T_th, at which crystallization is observed although it was hard to determine the precise T_th as observed in the laboratory time frame with the present computing resources. When N was relatively small (100 and 200), T_th was determined to be 320 and 350 K for the linear- and ring-polyethylene melts, respectively, while T_th was found to be 330 and 350 K, respectively, when N was 1000. Having confirmed that the crystallization of a ring-polyethylene melt occurs faster than that of the analogous linear melt, we conclude that the trans-segment-based crystallinities are effective for the analysis of local crystal behavior.

Folded alkane chains and the emergence of the lamellar crystal

The competition between chain stiffness and chain collapse gives rise to complex low temperature morphologies of single polymer chains, in our case alkanes. These structures are characterized by specific sequences of dihedral angles along the chain, i.e., dihedral angle correlations extending beyond local steric effects. To describe and classify these morphologies, one can transfer concepts from protein science, where this creation of dihedral angle correlations underlies the formation of α-helices and β-sheets. We show here by means of flat-histogram Monte Carlo simulations that, although lacking in primary structure being simple homopolymers, short alkane chains fold into non-trivial ground states (tertiary structure) consisting of chain segments of defined secondary structures. The folded lamellar crystal typical for polyethylene chains requires a minimum chain length to occur as the ground state folded structure, which we identify to be around 150 repeat units.

Chain Trajectory, Chain Packing, and Molecular Dynamics of Semicrystalline Polymers as Studied by Solid-State NMR

Chain-level structure of semicrystalline polymers in melt- and solution-grown crystals has been debated over the past half century. Recently, 13C⁻13C double quantum (DQ) Nuclear Magnetic Resonance (NMR) spectroscopy has been successfully applied to investigate chain-folding (CF) structure and packing structure of 13C enriched polymers after solution and melt crystallization. We review recent NMR studies for (i) packing structure, (ii) chain trajectory, (iii) conformation of the folded chains, (iv) nucleation mechanisms, (v) deformation mechanism, and (vi) molecular dynamics of semicrystalline polymers.

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.

Crystallization of finite-extensible nonlinear elastic Lennard-Jones coarse-grained polymers

The ability of a simple coarse-grained finite-extensible nonlinear elastic (FENE) Lennard-Jones (LJ) polymer model to be crystallized is investigated by molecular dynamics simulations. The optimal FENE Lennard-Jones parameter combinations (ratio between FENE and LJ equilibrium distances) and the optimal lattice parameters are calculated for five different perfect crystallite structures: simple tetragonal, body-centered tetragonal, body-centered orthorhombic, hexagonal primitive, and hexagonal close packed. It was found that the most energetically favorable structure is the body-centered orthorhombic. Starting with an equilibrated polymer liquid and with the optimal parameters found for the body-centered orthorhombic, an isothermal treatment led to the formation of large lamellar crystallites with a typical chain topology: folded, loop, and tie chains, and with a crystallinity of about 60%-70%, similar to real semicrystalline polymers. This simple coarse-grained Lennard-Jones model provides a qualitative tool to study semicrystalline microstructures for polymers.

Stress-Strain Relationships in Hydroxyl Substituted Polyethylene

Stress-strain relationships in semicrystalline hydroxylated polyethylene are studied using all-atom molecular dynamics simulations. Chain sizes ranging from 50 to 2000 carbons are gradually cooled from melt in order to obtain semicrystalline samples for pure, 4%, and 8% hydroxylated chains. Local orientational order of the polymer backbone and hydrogen bonding behavior is studied. The effects of -OH substitution and chain length on stress-strain relationships are examined at 300 K. The number of hydrogen bonds is found to be independent of the chain length. Stress-strain relationships are generally unaffected by 4% hydroxyl substitution in long chain polyethylene. The presence of 8% -OH tends to increase the elastic limit of the material. A method for comparing semicrystalline samples of substituted and unsubstituted polymeric chains is presented by eliminating differences in alignment, distribution, and extent of crystallization.

Chain Trajectory of Semicrystalline Polymers As Revealed by Solid-State NMR Spectroscopy

Over the last half century, a chain-folding structure of semicrystalline polymers has been debated in polymer science. Recently, ¹³C-¹³C double quantum (DQ) NMR spectroscopy combined with ¹³C selective isotope labeling has been developed to investigate re-entrance sites of the folded chains, mean values of adjacent re-entry number ⟨n⟩ and fraction ⟨F⟩ of semicrystalline polymers. This viewpoint highlights the versatile approaches of using solid-state (ss) NMR and isotope labeling for revealing (i) chain trajectory in melt- and solution-grown crystals, (ii) conformation of the folded chains in single crystals, (iii) self-folding in the early stage of crystallization, and (iv) unfolding of the folded chains under stretching.

Molecular dynamics study of the isothermal crystallization mechanism of polyethylene chain: the combined effects of chain length and temperature

A molecular level understanding of the polyethylene (PE) crystallization process was elucidated by molecular dynamics simulation of three states, with varying chain length and temperature. The process can be classified into the following three states: (1) nucleation controlled state, (2) competitive state of crystal growth process and new nuclei formation, and (3) crystal growth controlled state, which could be quantified by the evolution of nuclei number. With increasing chain length, two phenomena occur: the single crystallization mechanism changes from state (1) to (3), and the crystal size increases while the b/a axial ratio in the lateral surface decreases. These changes can be explained from a thermodynamic point of view, in that the van der Waals (vdW) interaction per CH2 unit is strengthened and more nucleation sites are generated for longer chain. Size effect (meaning different surface fractions when the chain collapses into a globule) was an important factor determining vdW energy per unit and the crystallization states of a single PE chain. On the other hand, the crystallization states were independent of chain length for short chains systems with the same size effect. In both conditions, a long chain generates multi-crystal domains, and a short chain prefers a single crystal domain. Our results not only provide molecular level evidence for crystallization states but also clarify the influence of chain length on the crystallization process.

Effect of chain stiffness on the competition between crystallization and glass-formation in model unentangled polymers

We map out the solid-state morphologies formed by model soft-pearl-necklace polymers as a function of chain stiffness, spanning the range from fully flexible to rodlike chains. The ratio of Kuhn length to bead diameter (lK/r0) increases monotonically with increasing bending stiffness kb and yields a one-parameter model that relates chain shape to bulk morphology. In the flexible limit, monomers occupy the sites of close-packed crystallites while chains retain random-walk-like order. In the rodlike limit, nematic chain ordering typical of lamellar precursors coexists with close-packing. At intermediate values of bending stiffness, the competition between random-walk-like and nematic chain ordering produces glass-formation; the range of kb over which this occurs increases with the thermal cooling rate |Ṫ| implemented in our molecular dynamics simulations. Finally, values of kb between the glass-forming and rodlike ranges produce complex ordered phases such as close-packed spirals. Our results should provide a useful initial step in a coarse-grained modeling approach to systematically determining the effect of chain stiffness on the crystallization-vs-glass-formation competition in both synthetic and colloidal polymers.

Folding of Polymer Chains in the Early Stage of Crystallization

Understanding the structure formation of an ordered domain in the early stage of crystallization is one of the long-standing issues in polymer science. In this study, we investigate the chain trajectory of isotactic polypropylene (iPP) formed via rapid and deep quenching, using solid-state NMR spectroscopy. Comparisons of experimental and simulated ¹³C-¹³C double quantum (DQ) buildup curves demonstrated that individual iPP chains adopt adjacent re-entry sequences with an average folding number ⟨n⟩ = 3-4 in the mesomorphic form, assuming an adjacent re-entry fraction ⟨F⟩ of 100%. Therefore, long flexible polymer chains naturally fold in the early stage of crystallization, and folding-initiated nucleation results in formation of mesomorphic nanodomains.

Three-Dimensional Conformation of Folded Polymers in Single Crystals

The chain-folding mechanism and structure of semicrystalline polymers have long been controversial. Solid-state NMR was applied to determine the chain trajectory of (13)C CH3-labeled isotactic poly(1-butene) (iPB1) in form III chiral single crystals blended with nonlabeled iPB1 crystallized in dilute solutions under low supercooling. An advanced (13)C-(13)C double-quantum NMR technique probing the spatial proximity pattern of labeled (13)C nuclei revealed that the chains adopt a three-dimensional (3D) conformation in single crystals. The determined results indicate a two-step crystallization process of (i) cluster formation via self-folding in the precrystallization stage and (ii) deposition of the nanoclusters as a building block at the growth front in single crystals.

Kinetic mechanism of chain folding in polymer crystallization

I develop a kinetic mechanism to explain chain folding in polymer crystallization which is based on the competition between the formation of stems, which is due to frequent occupations of trans states along the chains in the supercooled polymer melt, and the random coil structure of the polymer chains. Setting equal the average formation time of stems of length d(l) with the Rouse time of a piece of polymer of the same arc length d(l) yields a lower bound for the thickness of stems and bundles. The estimated lamellar thickness is inversely proportional to the supercooling. The present approach emphasizes the importance of repulsive interactions in polymer crystallization, which are expected to be responsible for the logarithmic lamellar thickening and the increase of lamellar thickness with pressure. An expression for the growth rate for formation and deposition of stems is derived by considering the growth as a dynamic multistage process.

Melting processes of oligomeric α and β isotactic polypropylene crystals at ultrafast heating rates

The melting behaviors of α (stable) and β (metastable) isotactic polypropylene (iPP) crystals at ultrafast heating rates are simulated with atomistic molecular dynamics method. Quantitative information about the melting processes of α- and β-iPP crystals at atomistic level is achieved. The result shows that the melting process starts from the interfaces of lamellar crystal through random dislocation of iPP chains along the perpendicular direction of lamellar crystal structure. In the melting process, the lamellar crystal gradually expands but the corresponding thickness decreases. The analysis shows that the system expansion lags behind the crystallinity decreasing and the lagging extents for α- and β-iPP are significantly different. The apparent melting points of α- and β-iPP crystals rise with the increase of the heating rate and lamellar crystal thickness. The apparent melting point of α-iPP crystal is always higher than that of β-iPP at differently heating rates. Applying the Gibbs-Thomson rule and the scaling property of the melting kinetics, the equilibrium melting points of perfect α- and β-iPP crystals are finally predicted and it shows a good agreement with experimental result.