(309g) Coarse-Grained Molecular Dynamics Simulations of Polymer Melts: Crossover from the Rouse to Reptation Regime | AIChE

(309g) Coarse-Grained Molecular Dynamics Simulations of Polymer Melts: Crossover from the Rouse to Reptation Regime

Authors 

Depa, P. K. - Presenter, Pennsylvania State University
Maranas, J. K. - Presenter, Pennsylvania State University


Coarse-graining that involves eliminating uninteresting degrees of freedom to reach large length and long time scales has found tremendous importance in polymer simulations. Here, we present results from long-time coarse-grained molecular dynamics simulations (CGMD) of linear polyethylene chains, ranging in length from C50 to C250. The coarse-grained force field used in the simulations of these systems is parameterized using the distribution functions from the united atom simulation (UA) of the shortest chains (C50). While CGMD of the C50 chains reproduces the static properties, compared to UA, the dynamics are off-set by a constant factor, a, corresponding to an effective time that is larger than the actual time in the simulation. We borrow from the framework of Accelerated Molecular Dynamics method to predict this rescaling factor. This apparent rescaling of time is independent of the chain length and is related to the decrease in the inter-molecular interaction when moving from UA to CG description. Thus, with this constant speed-up factor, a, obtained by the comparison of dynamic properties from UA and CG simulations of C50 chains, we simulate chains up to a length of 250 carbons with the coarse-grained force field. Static properties, including relevant distribution functions and the chain size are comparable to theoretical predictions and experiments. Dynamic properties like scattering functions and diffusion coefficients from these CGMD are also in good agreement with those from experiments, and other atomistic and coarse-grained simulations. We observe a crossover from Rouse to Reptation dynamics near a chain length of 150 carbons. Self diffusion coefficients (D) calculated with the scaled time exhibit a change in the power law dependence on chain length (D ~ N): while chains smaller than C150 exhibit rouse-like behavior (λ = 1), longer chains display reptation-like behavior (λ = 2.3).