A topology preserving method for generating equilibrated polymer melts in computer simulations

Polymorphism and Perfection in Crystallization of Hard Sphere Polymers

We present results on polymorphism and perfection, as observed in the spontaneous crystallization of freely jointed polymers of hard spheres, obtained in an unprecedentedly long Monte Carlo (MC) simulation on a system of 54 chains of 1000 monomers. Starting from a purely amorphous configuration, after an initial dominance of the hexagonal closed packed (HCP) polymorph and a transitory random hexagonal close packed (rHCP) morphology, the system crystallizes in a final, stable, face centered cubic (FCC) crystal of very high perfection. An analysis of chain conformational characteristics, of the spatial distribution of monomers and of the volume accessible to them shows that the phase transition is caused by an increase in translational entropy that is larger than the loss of conformational entropy of the chains in the crystal, compared to the amorphous state. In spite of the significant local re-arrangements, as reflected in the bending and torsion angle distributions, the average chain size remains unaltered during crystallization. Polymers in the crystal adopt ideal random walk statistics as their great length renders local conformational details, imposed by the geometry of the FCC crystal, irrelevant.

Simu-D: A Simulator-Descriptor Suite for Polymer-Based Systems under Extreme Conditions

We present Simu-D, a software suite for the simulation and successive identification of local structures of atomistic systems, based on polymers, under extreme conditions, in the bulk, on surfaces, and at interfaces. The protocol is built around various types of Monte Carlo algorithms, which include localized, chain-connectivity-altering, identity-exchange, and cluster-based moves. The approach focuses on alleviating one of the main disadvantages of Monte Carlo algorithms, which is the general applicability under a wide range of conditions. Present applications include polymer-based nanocomposites with nanofillers in the form of cylinders and spheres of varied concentration and size, extremely confined and maximally packed assemblies in two and three dimensions, and terminally grafted macromolecules. The main simulator is accompanied by a descriptor that identifies the similarity of computer-generated configurations with respect to reference crystals in two or three dimensions. The Simu-D simulator-descriptor can be an especially useful tool in the modeling studies of the entropy- and energy-driven phase transition, adsorption, and self-organization of polymer-based systems under a variety of conditions.

Accelerating the crystallization kinetics of linear polylactides by adding cyclic poly (L-lactide): Nucleation, plasticization and topological effects

Polylactide is one of the most versatile biopolymers, but its slow crystallization limits its temperature usage range. Hence finding ways to enhance it is crucial to widen its applications. Linear and cyclic poly (L-lactide) (l-PLLA and c-PLLA) of similarly low molecular weights (MW) were synthesized by ring-opening polymerization of L-lactide, and ring-expansion methodology, respectively. Two types of blends were prepared by solution mixing: (a) l-PLLA/c-PLLA, at extreme compositions (rich in linear or in cyclic chains), and (b) blends of each of these low MW materials with a commercial high MW linear PLA. The crystallization of the different blends was evaluated by polarized light optical microscopy and differential scanning calorimetry. It was found, for the first time, that in the l-PLLA rich blends, small amounts of c-PLLA (i.e., 5 and 10 wt%) increase the nucleation density, nucleation rate (1/τ₀), spherulitic growth rate (G), and overall crystallization rate (1/τ_50%), when compared to neat l-PLLA, due to a synergistic effect (i.e., nucleation plus plasticization). In contrast, the opposite effect was found in the c-PLLA rich blends. The addition of small amounts of l-PLLA to a matrix of c-PLLA chains causes a decrease in the nucleation density, 1/τ₀, G, and 1/τ_50% values, due to threading effects between cyclic and linear chains. Small amounts of l-PLLA and c-PLLA enhance the crystallization ability of a commercial high MW linear PLA without affecting its melting temperature. The l-PLLA only acts as a plasticizer for the PLA matrix, whereas c-PLLA has a synergistic effect in accelerating the crystallization of PLA that goes beyond simple plasticization. The addition of small amounts of c-PLLA affects not only PLA crystal growth but also its nucleation due to the unique cyclic chains topology.

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.

Multiscale approach to equilibrating model polymer melts

We present an effective and simple multiscale method for equilibrating Kremer Grest model polymer melts of varying stiffness. In our approach, we progressively equilibrate the melt structure above the tube scale, inside the tube and finally at the monomeric scale. We make use of models designed to be computationally effective at each scale. Density fluctuations in the melt structure above the tube scale are minimized through a Monte Carlo simulated annealing of a lattice polymer model. Subsequently the melt structure below the tube scale is equilibrated via the Rouse dynamics of a force-capped Kremer-Grest model that allows chains to partially interpenetrate. Finally the Kremer-Grest force field is introduced to freeze the topological state and enforce correct monomer packing. We generate 15 melts of 500 chains of 10.000 beads for varying chain stiffness as well as a number of melts with 1.000 chains of 15.000 monomers. To validate the equilibration process we study the time evolution of bulk, collective, and single-chain observables at the monomeric, mesoscopic, and macroscopic length scales. Extension of the present method to longer, branched, or polydisperse chains, and/or larger system sizes is straightforward.