Crystallization and Melting Simulations of Oligomeric α1 Isotactic Polypropylene

Local and Global Order in Dense Packings of Semi-Flexible Polymers of Hard Spheres

The local and global order in dense packings of linear, semi-flexible polymers of tangent hard spheres are studied by employing extensive Monte Carlo simulations at increasing volume fractions. The chain stiffness is controlled by a tunable harmonic potential for the bending angle, whose intensity dictates the rigidity of the polymer backbone as a function of the bending constant and equilibrium angle. The studied angles range between acute and obtuse ones, reaching the limit of rod-like polymers. We analyze how the packing density and chain stiffness affect the chains' ability to self-organize at the local and global levels. The former corresponds to crystallinity, as quantified by the Characteristic Crystallographic Element (CCE) norm descriptor, while the latter is computed through the scalar orientational order parameter. In all cases, we identify the critical volume fraction for the phase transition and gauge the established crystal morphologies, developing a complete phase diagram as a function of packing density and equilibrium bending angle. A plethora of structures are obtained, ranging between random hexagonal closed packed morphologies of mixed character and almost perfect face centered cubic (FCC) and hexagonal close-packed (HCP) crystals at the level of monomers, and nematic mesophases, with prolate and oblate mesogens at the level of chains. For rod-like chains, a delay is observed between the establishment of the long-range nematic order and crystallization as a function of the packing density, while for right-angle chains, both transitions are synchronized. A comparison is also provided against the analogous packings of monomeric and fully flexible chains of hard spheres.

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.

A coarse-grained model for capturing the helical behavior of isotactic polypropylene

Understanding the process-property relations of helical polymers using molecular simulations has been an attractive research field over the years. Specifically, isotactic polypropylene still remains a challenge for current computational experimentation, as it exhibits phenomena such as crystallization that emerge on large spatial and temporal scales. Coarse-graining is an efficient technique for approaching such phenomena, although previous coarse-grained models lack in preserving important atomistic and structural details. In this paper we develop a new coarse-grained model, based on the popular MARTINI force field, that is able to reproduce the helical behavior of isotactic polypropylene. To test the model, the predicted statistical and structural properties (characteristic ratio, density, entanglement molecular weight, solubility parameter in the melt) are compared with previous simulation results and available experimental data. For the development of the new coarse-grained force field, a single unperturbed chain Monte Carlo algorithm has been implemented: an efficient algorithm which samples conformations representative of a melt by simulating just a single chain.

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.

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.

Modelling flow-induced crystallisation in polymers

A review of recent multiscale modelling of flow-induced crystallisation in polymers with a particular emphasis on newly emerging techniques to connect modelling and simulation techniques at different levels of coarse-graining.