Computer modeling of polymer crystallization – Toward computer-assisted materials' design

Molecular dynamics simulations of single polyethylene chain folding during fast quenching using all-atom and united-atom models

Molecular dynamics simulations have been employed to investigate the folding behavior of a single linear polyethylene (PE) chain containing 1000 backbone carbon atoms under fast quenching based on all-atom and united-atom models. The single-chain folding characteristics were studied in detail for six different force fields by analyzing the evolution of chain conformations, folded structure characterisation, free energy and crystallisation. The results show that the all-trans chain undergoes a similar two-stage chain collapse mechanism during isothermal relaxation at T = 500 K, transitioning from local collapse to global collapse into a molten globule state under different force fields. During fast quenching at 100 K ns-1, the molten globule of all-atom model transitions into a folded, significantly anisotropic ordered structure under AMBER-AA or OPLS-AA force fields, while that of the united-atom model remains unchanged in its globular structure. The chain crystallization evolution indicates that the single chain folds into ordered lamellar structures with higher crystallinity under AMBER-AA and OPLS-AA force fields. In contrast, under the other four force fields, the single chain remains in a stable amorphous state.

Prediction of directional solidification in freeze casting of biomaterial scaffolds using physics-informed neural networks

Freeze casting, a manufacturing technique widely applied in biomedical fields for fabricating biomaterial scaffolds, poses challenges for predicting directional solidification due to its highly nonlinear behavior and complex interplay of process parameters. Conventional numerical methods, such as computational fluid dynamics (CFD), require adequate and accurate boundary condition knowledge, limiting their utility in real-world transient solidification applications due to technical limitations. In this study, we address this challenge by developing a physics-informed neural networks (PINNs) model to predict directional solidification in freeze-casting processes. The PINNs model integrates physical constraints with neural network predictions, requiring significantly fewer predetermined boundary conditions compared to CFD. Through a comparison with CFD simulations, the PINNs model demonstrates comparable accuracy in predicting temperature distribution and solidification patterns. This promising model achieves such a performance with only 5000 data points in space and time, equivalent to 250,000 timesteps, showcasing its ability to predict solidification dynamics with high accuracy. The study's major contributions lie in providing insights into solidification patterns during freeze-casting scaffold fabrication, facilitating the design of biomaterial scaffolds with finely tuned microstructures essential for various tissue engineering applications. Furthermore, the reduced computational demands of the PINNs model offer potential cost and time savings in scaffold fabrication, promising advancements in biomedical engineering research and development.

Unveiling the Nanoconfinement Effect on Crystallization of Semicrystalline Polymers Using Coarse-Grained Molecular Dynamics Simulations

Semicrystalline polymers under nanoconfinement show distinct structural and thermomechanical properties compared to their bulk counterparts. Despite extensive research on semicrystalline polymers under nanoconfinement, the nanoconfinement effect on the local crystallization process and the unique structural evolution of such polymers have not been fully understood. In this study, we unveil such effects by using coarse-grained molecular dynamics simulations to study the crystallization process of a model semicrystalline polymer-polyvinyl alcohol (PVA)-under different levels of nanoconfinement induced by nanoparticles that are represented implicitly. We quantify in detail the evolution of the degree of crystallinity (XC) of PVA and examine distinct crystalline regions from simulation results. The results show that nanoconfinement can promote the crystallization process, especially at the early stage, and the interfaces between nanoparticles and polymer can function as crystallite nucleation sites. In general, the final XC of PVA increases with the levels of nanoconfinement. Further, nanoconfined cases show region-dependent XC with higher and earlier increase of XC in regions closer to the interfaces. By tracking region-dependent XC evolution, our results indicate that nanoconfinement can lead to a heterogenous crystallization process with a second-stage crystallite nucleation in regions further away from the interfaces. In addition, our results show that even under very high cooling rates, the nanoconfinement still promotes the crystallization of PVA. This study provides important insights into the underlying mechanisms for the intricate interplay between nanoconfinement and the crystallization behaviors of semicrystalline polymer, with the potential to guide the design and characterization of semicrystalline polymer-based nanocomposites.

Material Property Recovery by Controlling the Melt Memory Effects on Recrystallization and on Crystal Deformation: An Approach by the Molecular Dynamics Simulation for Polyethylene

Degradation in the mechanical properties of recycled polymer materials has been recently appearing as a big issue in polymer science. The molecular mechanism of the degradation is considered in part due to residual memories of flow experienced during molding processes, and therefore the mechanical recycling through remolding involving melting and recrystallization has been attempted in recent years. In the present paper, the molecular processes of melting and recrystallization are investigated by the molecular dynamics simulation for polyethylene with special interest in the melt memory effects. We also studied the mechanical properties of the recrystalized materials that have undergone different recrystallization processes aiming to discover better recycling strategies. A successive step-by-step approach is adopted to study the loss of the crystal memory during retention in the melt, the effects of the melt memory on the mode of recrystallization, the relation between the recrystallization mode and the resulting higher-order structure, and the mechanical properties controlled by the higher-order structures. It is shown that the melt memory clearly remains in various order parameters that persist over time scales corresponding to the Rouse time, the remaining melt memory markedly affects the crystallization mode leading to distinct crystalline morphologies, and the distinct morphologies obtained give different mechanical responses during large deformations.

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.

Comparing equilibration schemes of high-molecular-weight polymer melts with topological indicators

Recent theoretical studies have demonstrated that the behaviour of molecular knots is a sensitive indicator of polymer structure. Here, we use knots to verify the ability of two state-of-the-art algorithms-configuration assembly and hierarchical backmapping-to equilibrate high-molecular-weight (MW) polymer melts. Specifically, we consider melts with MWs equivalent to several tens of entanglement lengths and various chain flexibilities, generated with both strategies. We compare their unknotting probability, unknotting length, knot spectra, and knot length distributions. The excellent agreement between the two independent methods with respect to knotting properties provides an additional strong validation of their ability to equilibrate dense high-MW polymeric liquids. By demonstrating this consistency of knotting behaviour, our study opens the way for studying topological properties of polymer melts beyond time and length scales accessible to brute-force molecular dynamics simulations.

Coarse-grained molecular dynamics simulations of void generation and growth processes in the fracture of the lamellar structure of polyethylene

We investigate the void generation and growth process in semicrystalline polymers by large-scale coarse-grained molecular dynamics simulations. Voids are generated in the amorphous layers and grow spherically and then cylindrically, consistent with the results of previous experiments. Interestingly, the fusion of voids is indicated to enlarge the voids in the direction perpendicular to the stretching direction, but not beyond the crystalline layers. The orientational order along the stretching direction increased sharply before void generation, but the increase leveled off afterward. The simulations also clearly reveal that the crystalline layers break but do not bend in the fragmentation process. The dependence of the void growth process on stretching velocity and the stability levels of voids at constant strain are also discussed.

Insight into the Microcosm of the Human Growth Hormone and Its Interactions with Polymers and Copolymers: A Molecular Dynamics Perspective

Therapeutic proteins nowadays have increasingly been applied for their considerable potential in treating a wide variety of diseases. The effectiveness and potency of native therapeutic proteins are limited by various factors (e.g., stability, blood circulation time, specificity). Over the past years, a great deal of effort has been devoted to developing safe and efficient protein delivery systems. Entrapment of protein into polymeric and copolymeric matrices is common among the different types of protein-based drug formulation. However, despite the massive efforts toward developing therapeutic protein delivery in experimental studies and industrial applications, there is relatively little data on the influence of polymers and copolymers on therapeutic proteins at the atomic and molecular levels. Herein, molecular dynamics (MD) simulations are used to study the effects of biocompatible synthetic polymers including methoxy poly(ethylene glycol) (MPEG), poly(lactic acid) (PLA), and poly(lactic acid) copolymers (poly(lactic-co-glycolic acid)) PLGA and MPEG-PLA(PELA)) on the structure and dynamics of the human growth hormone (hGH), and the results are compared with previous experimental findings. Our results indicate that the hGH conformation remains stable both in pure water and in the presence of polymers, and these results are in good agreement with previous experimental data. It is shown that the MPEG chains are self-assembled and folded into a semicrystalline structure; therefore, only a small portion of the protein interacts with the polymer. The other three polymers, however, interact well with the protein and partially cover its surface. Our findings suggest that the use of these polymers for protein encapsulation has the advantage of reducing protein aggregation and thus increasing drug serum half-life. Eventually, we anticipate that the research results will expand the current knowledge about encapsulation mechanisms and the molecular interactions between hGH and the polymers.

Predicting Crystallization Tendency of Polymers Using Multifidelity Information Fusion and Machine Learning

The degree of crystallinity of a polymer is a critical parameter that controls a variety of polymer properties. A high degree of crystallinity is associated with excellent mechanical properties crucial for high-performing applications like composites. Low crystallinity promotes ion and gas mobility critical for battery and membrane applications. Experimental determination of the crystallinity for new polymers is time and cost intensive. A data-driven machine learning-based method capable of rapidly predicting the crystallinity could counter these disadvantages and be used to screen polymers for a myriad of applications in a fast, inexpensive fashion. In this work, we developed the first-of-its-kind, data-driven machine learning model to predict the most-likely polymer crystallinity trained on experimental data and theoretical group contribution methods. Since polymer data under consistent processing conditions are unavailable, we tackled process variability by using the "most-likely" polymer values which we refer to as the polymer's tendency to crystallize. Experimental data for polymers' tendency to crystallize is limited by number and diversity, and to tackle this, we augmented experimentation-based data with data using group contribution methods. Therefore, this work utilized two data sets, viz., a high-fidelity, experimental data set for 107 polymers and a more diverse, less accurate low-fidelity data set for 429 polymers which used group contribution methods. We used a multifidelity information fusion strategy to utilize all the information captured in the low-fidelity data set while still predicting at the high-fidelity accuracy. Although this model inherently assumed "typical" processing conditions and estimated the "most-likely" percent crystallinity value, it can help in the estimation of a polymer's tendency to crystallize in a far more cost-effective and efficient manner.

Magnetic and Dielectric Properties of Nano- and Micron-BiFeO3/LDPE Composites with Magnetization Treatments

By means of magnetization treatments at ambient temperature and elevated temperatures, the nano- and micron-bismuth ferrate/low density polyethylene (BiFeO₃/LDPE) dielectric composites are developed to explore the material processing method to modify the crystalline morphology, magnetic and dielectric properties. The magnetic field treatment can induce the dipole in the LDPE macromolecular chain which leads to preferred orientation of polyethylene crystal grains to the direction of the magnetization field. The surface morphology of the materials measured by atomic force microscope (AFM) implies that the LDPE macromolecular chains in BiFeO₃/LDPE composites have been orderly arranged and form thicker lamellae accumulated with a larger spacing after high temperature magnetization, resulting in the increased dimension and orientation of spherulites. The residual magnetization intensities of BiFeO₃/LDPE composites have been significantly improved by magnetization treatments at ambient temperature. After this magnetization at ambient temperature, the M_R of nano- and micron-BiFeO₃/LDPE composites approach to 4.415 × 10⁻³ and 0.690 × 10⁻³ emu/g, respectively. The magnetic moments of BiFeO₃ fillers are arranged parallel to the magnetic field direction, leading to appreciable enhancement of the magnetic interactions between BiFeO₃ fillers, which will inhibit the polarization of the electric dipole moments at the interface between BiFeO₃ fillers and the LDPE matrix. Therefore, magnetization treatment results in the lower dielectric constant and higher dielectric loss of BiFeO₃/LDPE composites. It is proven that the magnetic and dielectric properties of polymer dielectric composites can be effectively modified by the magnetization treatment in the melt blending process of preparing composites, which is expected to provide a technical strategy for developing magnetic polymer dielectrics.

Anomalous Ostwald Ripening Enables 2D Polymer Crystals via Fast Evaporation

We demonstrate by molecular simulations that the Ostwald ripening of crystalline polymer nuclei within the fast-evaporation-induced 2D skin layer is retarded at suitable temperatures and evaporation rates. Such an anomalous ripening can be attributed to the interplay between the thermodynamically driven diffusion of noncrystalline fragments toward the growing nuclei and the diffusive current away from the free surface caused by the densification in the nonequilibrium skin layer. The growth orientation of the nuclei inside the skin plane can be adjusted during this anomalous ripening process, which is beneficial for fabricating 2D polymer crystals.

Dominant Effects of Short-Chain Branching on the Initial Stage of Nucleation and Formation of Tie Chains for Bimodal Polyethylene as Revealed by Molecular Dynamics Simulation

The molecular mechanism of short-chain branching (SCB), especially the effects of methylene sequence length (MSL) and short-chain branching distribution (SCBD) on the initial stage of nucleation, the crystallization process, and particularly the tie chain formation process of bimodal polyethylene (BPE), were explored using molecular dynamics simulation. This work constructed two kinds of BPE models in accordance with commercial BPE pipe resins: SCB incorporated in the long chain or in the short chains. The initial stage of nucleation was determined by the MSL of the system, as the critical MSL for a branched chain to nucleate is about 60 CH2. SCB incorporated in the long chain led to a delay of the initial stage of nucleation relative to the case of SCB incorporated in the short chains. The increase of branch length could accelerate the delay to nucleation. The location of short chain relative to the long chain depended on the MSL of the short chain. As the MSL of the system decreased, the crystallinity decreased, while the tie chains concentration increased. The tie chains concentration of the BPE model with branches incorporated in the long chain was higher than that with branches incorporated in the short chain.

Effect of Mold Geometry on Pore Size in Freeze-Cast Chitosan-Alginate Scaffolds for Tissue Engineering

Freeze-casting is a popular method to produce biomaterial scaffolds with highly porous structures. The pore structure of freeze-cast biomaterial scaffolds is influenced by processing parameters but has mostly been controlled experimentally. A mathematical model integrating Computational Fluid Dynamics with Population Balance Model was developed to predict average pore size (APS) of 3D porous chitosan-alginate scaffolds and to assess the influence of the geometrical parameters of mold on scaffold pore structure. The model predicted the crystallization pattern and APS for scaffolds cast in different diameter molds and filled to different heights. The predictions demonstrated that the temperature gradient and solidification pattern affect ice crystal nucleation and growth, subsequently influencing APS homogeneity. The predicted APS compared favorably with APS measurements from a corresponding experimental dataset, validating the model. Sensitivity analysis was performed to assess the response of the APS to the three geometrical parameters of the mold: well radius; solution fill height; and spacing between wells. The pore size was most sensitive to the distance between the wells and least sensitive to solution height. This validated model demonstrates a method for optimizing the APS of freeze-cast biomaterial scaffolds that could be applied to other compositions or applications.

The Origins of Enhanced and Retarded Crystallization in Nanocomposite Polymers

Controlling the crystallinity of hybrid polymeric systems has an important impact on their properties and is essential for developing novel functional materials. The crystallization of nanocomposite polymers with gold nanoparticles is shown to be determined by free space between nanoparticles. Results of large-scale molecular dynamics simulations reveal while crystallinity is affected by the nanoparticle size and its volume fraction, their combined effects can only be measured by interparticle free space and characteristic size of the crystals. When interparticle free space becomes smaller than the characteristic extended length of the polymer molecule, nanoparticles impede the crystallization because of the confinement effects. Based on the findings from this work, equations for critical particle size or volume fraction that lead to this confinement-induced retardation of crystallization are proposed. The findings based on these equations are demonstrated to agree with the results reported in experiments for nanocomposite systems. The results of simulations also explain the origin of a two-tier crystallization regime observed in some of the hybrid polymeric systems with planar surfaces where the crystallization is initially enhanced and then retarded by the presence of nanoparticles.

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.

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.

Molecular dynamics simulation of the folding of single alkane chains with different lengths on single-walled carbon nanotubes and graphene

Chain folding is an important step during polymer crystallization. In order to study the effects of the surface on chain folding, molecular dynamics simulations of the folding of different alkane chains on three kinds of single-walled carbon nanotubes (SWCNTs) and graphene were performed. The folding behaviors of the single alkane chains on these surfaces were found to be different from their folding behaviors in vacuum. The end-to-end distances of the chains were calculated to explore the chain folding. An increasing tendency to fold into two or more stems with increasing alkane chain length was observed. This result indicates that the occurrence and the stability of chain folding are related to the surface curvature, the diameter of the SWCNT, and surface texture. In addition, the angle between the direction of the alkane chain segment and the direction of the surface texture was measured on different surfaces.

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.

Dynamics and Structure of Monolayer Polymer Crystallites on Graphene

Graphene-based nanostructured systems and van der Waals heterostructures comprise a material class of growing technological and scientific importance. Joining materials with vastly different properties, polymer/graphene heterosystems promise diverse applications in surface and nanotechnology, including photovoltaics or nanotribology. Fundamentally, molecular adsorbates are prototypical systems to study confinement-induced phase transitions exhibiting intricate dynamics, which require a comprehensive understanding of the dynamical and static properties on molecular time and length scales. Here, we investigate the dynamics and the structure of a single polyethylene chain on free-standing graphene by means of molecular dynamics simulations. In equilibrium, the adsorbed polymer is orientationally linked to the graphene as two-dimensional folded-chain crystallite or at elevated temperatures as a floating solid. The associated superstructure can be reversibly melted on a picosecond time scale upon quasi-instantaneous substrate heating, involving ultrafast heterogeneous melting via a transient floating phase. Our findings elucidate time-resolved molecular-scale ordering and disordering phenomena in individual polymers interacting with solids, yielding complementary information to collective friction and viscosity, and linking to recent experimental observables from ultrafast electron diffraction. We anticipate that the approach will help in resolving nonequilibrium phenomena of hybrid polymeric systems over a broad range of time and length scales.

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.

A Phase Field Technique for Modeling and Predicting Flow Induced Crystallization Morphology of Semi-Crystalline Polymers

Flow induced crystallization of semi-crystalline polymers is an important issue in polymer science and engineering because the changes in morphology strongly affect the properties of polymer materials. In this study, a phase field technique considering polymer characteristics was established for modeling and predicting the resulting morphologies. The considered crystallization process can be divided into two stages, which are nucleation upon the flow induced structures and subsequent crystal growth after the cessation of flow. Accordingly, the proposed technique consists of two parts which are a flow induced nucleation model based on the calculated information of molecular orientation and stretch, and a phase field crystal growth model upon the oriented nuclei. Two-dimensional simulations are carried out to predict the crystallization morphology of isotactic polystyrene under an injection molding process. The results of these simulations demonstrate that flow affects crystallization morphology mainly by producing oriented nuclei. Specifically, the typical skin-core structures along the thickness direction can be successfully predicted. More importantly, the results reveal that flow plays a dominant part in generating oriented crystal morphologies compared to other parameters, such as anisotropy strength, crystallization temperature, and physical noise.

On the Orientation-Induced Crystallization of Polymers

In order to understand orientation-induced crystallization of polymers, we introduced an intermolecular interaction between polymer chains based on quantum mechanics. We therefore considered a pair of perfectly extended chains where the intermolecular interaction is assumed to be based on the hydrogen interaction with a single chain. When two protons of each extended chain become closer together under tension, the attractive force between the extended chains is caused by the interaction between hydrogen atoms surrounding the main chains based on the hydrogen molecule ion H 2 + . The energy is split into the ground and excited states, and the spontaneous process leading to the ground state is the origin for orientation-induced crystallization.

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.