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 CH2 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, Tth, at which crystallization is observed although it was hard to determine the precise Tth as observed in the laboratory time frame with the present computing resources. When N was relatively small (100 and 200), Tth was determined to be 320 and 350 K for the linear- and ring-polyethylene melts, respectively, while Tth 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.

Multipoint segmental repulsive potential for entangled polymer simulations with dissipative particle dynamics.

A model is developed for simulating entangled polymers by dissipative particle dynamics (DPD) using the segmental repulsive potential (SRP). In contrast to previous SRP models that define a single-point interaction on each bond, the proposed SRP model applies a dynamically adjustable multipoint on the bond. Previous SRP models could not reproduce the equilibrium properties of Groot and Warren's original DPD model [R. D. Groot and P. B. Warren, J. Chem. Phys. 107, 4423 (1997)] because the introduction of a single SRP induces a large excluded volume, whereas, the proposed multipoint SRP (MP-SRP) introduces a cylindrical effective excluded bond volume. We demonstrate that our MP-SRP model exhibits equilibrium properties similar to those of the original DPD polymers. The MP-SRP model parameters are determined by monitoring the number of topology violations, thermodynamic properties, and the polymer internal structure. We examine two typical DPD polymers with different bond-length distributions; one of them was used in the modified SRP model by Sirk et al. [J. Chem. Phys. 136, 134903 (2012)], whereas the other was used in the original DPD model. We demonstrate that for both polymers, the proposed MP-SRP model captures the entangled behaviors of a polymer melt naturally, by calculating the slowest relaxation time of a chain in the melt and the shear relaxation modulus. The results indicate that the proposed MP-SRP model can be applied to a variety of DPD polymers.

Thinning Approximation for Calculating Two-Dimensional Scattering Patterns in Dissipative Particle Dynamics Simulations under Shear Flow.

Modifications to improve thinning approximation (TA) were considered in order to calculate two-dimensional scattering patterns (2DSPs) for dissipative particle dynamics (DPD) simulations of polymer melts under a shear flow. We proposed multipoint TA and adaptive TA because the bond lengths in DPD chains vary widely when compared to those in Kremer⁻Grest (KG) chains, and the effectiveness of these two types of TA for the two major DPD parameter sets were investigated. In this paper, we report our findings on the original DPD model with soft bonds and that with rigid bonds. Based on the behavior of the 2DSPs and the distribution of orientations of the bond vectors, two spot patterns originating from the oriented chain correlations were observed when distinct distributions of the highly oriented bond vectors in the shear direction were obtained. For multipoint TA, we concluded that at least two additional midpoints ( n mid ≥ 2 ) are required to clearly observe the two spot patterns. For adaptive TA, a dividing distance of l ATA ≤ 0.4 is sufficient for clear observation, which is consistent with the requirement of n mid ≥ 2 for multipoint TA.