A highly coarse-grained model to simulate entangled polymer melts

We introduce a highly coarse-grained model to simulate the entangled polymer melts. In this model, a polymer chain is taken as a single coarse-grained particle, and the creation and annihilation of entanglements are regarded as stochastic events in proper time intervals according to certain rules and possibilities. We build the relationship between the probability of appearance of an entanglement between any pair of neighboring chains at a given time interval and the rate of variation of entanglements which describes the concurrence of birth and death of entanglements. The probability of disappearance of entanglements is tuned to keep the total entanglement number around the target value. This useful model can reflect many characteristics of entanglements and macroscopic properties of polymer melts. As an illustration, we apply this model to simulate the polyethylene melt of C(1000)H(2002) at 450 K and further validate this model by comparing to experimental data and other simulation results.     

Multicentered Gaussian-based potentials for coarse-grained polymer simulations: Linking atomistic and mesoscopic scales

A new analytical form for bond and angle potentials suitable for obtaining mesoscale effective force fields from target distributions is reported. Applications to realistic coarse-grained models of atactic polystyrene and polyamide-6,6 are described. The potential optimization procedure, despite its simplicity, allows the accurate reproduction of the target atomistic distributions. The procedure has been validated for both bond and angle potentials. Nonbonded numerical potentials have been optimized by pressure-corrected iterative Boltzmann inversion. Thus, the proposed coarse-graining strategy consists of hybrid analytical and numerical contributions to the mesoscale polymer force field. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 43: 871–885, 2005     

A molecular model for cohesive slip at polymer melt/solid interfaces

A molecular model is proposed which predicts wall slip by disentanglement of polymer chains adsorbed on a wall from those in the polymer bulk. The dynamics of the near-wall boundary layer is found to be governed by a nonlinear equation of motion, which accounts for such mechanisms on surface chains as convection, retraction, constraint release, and thermal fluctuations. This equation is valid over a wide range of grafting regimes, including those in which interactions between neighboring adsorbed molecules become essential. It is not closed since the dynamics of adsorbed chains is shown to be coupled to that of polymer chains in the bulk via constraint release. The constitutive equations for the layer and bulk, together with continuity of stress and velocity, are found to form a closed system of equations which governs the dynamics of the whole "bulk+boundary layer" ensemble. Its solution provides a stick-slip law in terms of the molecular parameters and extruder geometry. The model is quantitative and contains only those parameters that can be measured directly, or extracted from independent rheological measurements. The model predictions show a good agreement with available experimental data.     

Numerical simulation of polymer nanocomposites using self-consistent mean-field model

Clay-containing polymeric nanocomposites (PNC) are mixtures of dispersed clay platelets in a polymeric matrix. These materials show enhancement of physical properties, such as modulus, strength, and dimensional stability, as well as a reduction of gas permeability and flammability. The performance is related to the degree of clay dispersion (i.e., intercalation or exfoliation) and the bonding between the clay and the matrix. The main goal of this work has been to map the degree of dispersion as a function of independent variables (viz., magnitude of the interaction parameters, molecular weights, composition, etc.). In this paper, we present the results of the numerical analysis of the equilibrium thermodynamic miscibility using one- and two-dimensional (1D and 2D) models based on the self-consistent mean-field theory. In the limit, the 2D model reproduced the 1D model published results. The adopted 2D model considers the presence of four PNC components: solid clay platelets, low molecular weight intercalant, polymeric matrix, and end-functionalized compatibilizer. The simulations, with realistic values of the binary interaction parameters, were analyzed for potential exfoliation of PNC with a polyolefin as the matrix. The simulation results show that intercalation and exfoliation is expected within limited ranges of the independent variables. The presence of a bare clay surface (e.g., generated by thermal decomposition of intercalant or extraction by molten polymer) has a strong negative effect on the dispersion process. The simulation successfully identified the most influential factors, e.g., optimum ranges of the compatibilizer and the intercalant concentration.     

A molecular approach to the spurt effect in polymer melt flow

The Doi-Edwards theory of polymer melts, extended to include relaxation processes associated with chain-length equilibration, is used to make quantitative predictions of a discontinuity in the flow curve of a monodisperse melt in a capillary. A fluid interface between regions of high and low deformation rates is found to propagate from the former into the latter. Our results for the “spurt” and its hysteresis compare favorably with experiment using a molecular weight dependence of the ratio of “reptation time” (T_d) to “equilibration time” (T_eq) in agreement with that determined from nonlinear stress relaxation.     

Electronic shell structure and metal clusters: the self-consistent spheroidal jellium model

The electronic properties of small metal particles within the recently proposed self-consistent spheroidal jellium model [1] are further explored and compared to recent experimental data. Physical properties investigated include ionization potentials, electron affinities and the binding energy of neutral monomers to cationic clusters. The formalism is applied within the size-range 2≦N≦41, but could easily be extended beyondN=41. Finally, we discuss briefly the implications for the study of the dynamical response of open-shell clusters. In sharp contrast to earlier studies the functional is now corrected for self-interaction error, in a way first proposed by Pedew and Zunger [2]. This enables us to calculate reliable values for the electron affinities within ajellium-based model. This has the advantage, that we can calculate the affinities for Cu for all particle numbers for which experimental data are available. In all cases investigated we obtain excellent agreement with experiment, with pronounced shell-effects both for the electron affinities and for the binding energies, confirming in this way that the abundances map the relative stability of (Me) _N clusters, with Me being a sp-metal atom (Na, K, Li, Cu, Ag, Au etc.).     

Transferability of coarse-grained force fields: the polymer case

A key question for all coarse-graining methodologies is the degree of transferability of the resulting force field between various systems and thermodynamic conditions. Here we present a detailed study of the transferability over different thermodynamic states of a coarse-grained (CG) force field developed using the iterative Boltzmann inversion method. The force field is optimized against distribution functions obtained from atomistic simulations. We analyze the polymer case by investigating the bulk of polystyrene and polyamide-6,6 whose coarse-grained models differ in the chain length and in the number of atoms lumped in one bead. The effect of temperature and pressure on static, dynamic, and thermodynamic properties is tested by comparing systematically the coarse-grain results with the atomistic ones. We find that the CG model describing the polystyrene is transferable only in a narrow range of temperature and it fails in describing the change of the bulk density when temperature is 80 K lower than the optimization one. Moreover the calculation of the self-diffusion coefficient shows that the CG model is characterized by a faster dynamics than the atomistic one and that it overestimates the isothermal compressibility. On the contrary, the polyamide-6,6 CG model turns out to be fully transferable between different thermodynamic conditions. The transferability is checked by changing either the temperature or the pressure of the simulation. We find that, in this case, the CG model is able to follow all the intra- and interstructural rearrangements caused by the temperature changes. In addition, while at low temperature the difference between the CG and atomistic dynamics is remarkable due to the presence of hydrogen bonds in the atomistic systems, for high temperatures, the speedup of the CG dynamics is strongly reduced, leading to a CG diffusion coefficient only six times bigger than the atomistic one. Moreover, the isothermal compressibility calculated at different temperatures agrees very well with the experimental one. We find that the polymer chain length does not affect the transferability of the force field and we attribute such transferability mainly to the finer model used in describing the polyamide-6,6 than the polystyrene.     

The dynamics of single chains within a model polymer melt

Discontinuous molecular dynamics simulations are performed on a system containing 32 hard chains of length 192 at a volume fraction of phi = 0.45 to explore the idea that localized entanglements have a significant effect on the dynamics of the individual chains within an entangled polymer melt. Anomalous behavior can still be observed when studying the dynamics of the individual chains, although increased time averaging causes the anomalous relaxation-memory-release behavior that was observed previously in the system to smooth out. First, the individual chain mean squared displacements and apparent diffusion coefficients are calculated, and a wide distribution of diffusive behavior is found. Although the apparent diffusion coefficient curve averaged over all chains displays the predicted long-time diffusive behavior, the curves for the individual chains differ both qualitatively and quantitatively. They display superdiffusive, diffusive, and subdiffusive behavior, with the largest percentage of chains exhibiting superdiffusive behavior and the smallest percentage exhibiting the predicted diffusive behavior. Next, the individual chain end-to-end vector autocorrelation functions and relaxation times are determined, and a wide distribution of stress relaxation behavior is found. The times when the end-to-end vector autocorrelation functions relax completely span almost an order of magnitude in reduced time. For some chains, the end-to-end vector autocorrelation function relaxes smoothly toward zero similar to the system average; however, for other chains the relaxation is slowed greatly, indicating the presence of additional entanglements. Almost half of the chains exhibit the anomalous behavior in the end-to-end vector autocorrelation function. Finally, the dynamic properties are displayed for a single chain exhibiting anomalous relaxation-memory-release behavior, supporting the idea that the relaxation-memory-release behavior is a single-chain property.     

Statistical temperature molecular dynamics: application to coarse-grained beta-barrel-forming protein models

Recently the authors proposed a novel sampling algorithm, "statistical temperature molecular dynamics" (STMD) [J. Kim et al., Phys. Rev. Lett. 97, 050601 (2006)], which combines ingredients of multicanonical molecular dynamics and Wang-Landau sampling. Exploiting the relation between the statistical temperature and the density of states, STMD generates a flat energy distribution and efficient sampling with a dynamic update of the statistical temperature, transforming an initial constant estimate to the true statistical temperature T(U), with U being the potential energy. Here, the performance of STMD is examined in the Lennard-Jones fluid with diverse simulation conditions, and in the coarse-grained, off-lattice BLN 46-mer and 69-mer protein models, exhibiting rugged potential energy landscapes with a high degree of frustration. STMD simulations combined with inherent structure (IS) analysis allow an accurate determination of protein thermodynamics down to very low temperatures, overcoming quasiergodicity, and illuminate the transitions occurring in folding in terms of the energy landscape. It is found that a thermodynamic signature of folding is significantly suppressed by accurate sampling, due to an incoherent contribution from low-lying non-native IS in multifunneled landscapes. It is also shown that preferred accessibility to such IS during the collapse transition is intimately related to misfolding or poor foldability.     

Kinetics of molecular reactions in polymer solutions: Numerical model

The numerical simulation of the kinetics of the molecular interaction in a system containing mobile polymer clusters is performed on the basis of the bond-fluctuation model (BFM). Power-law dependences of the rates of molecular reactions on the concentration of these monomers are observed above the critical concentration. These dependences are related to the features of macromolecular structures in solutions.     

Internal structure of dendrimers in the melt under shear: a molecular dynamics study

The molecular structure of fluids composed of dendrimers of different generations is studied using nonequilibrium molecular dynamics (NEMD). NEMD results for dendrimer melts undergoing planar Couette flow are reported and analyzed with particular attention paid to the shear-induced changes in the internal structure of dendrimers. The radii of gyration, pair distribution functions and the fractal dimensionality of the dendrimers are determined at different strain rates. The location of the terminal groups is analyzed and found to be uniformly distributed throughout the space occupied by the molecules. The fractal dimension as a function of strain rate displays crossover behavior analogous to the Newtonian/non-Newtonian transition of shear viscosity.     

Application of mean-field model of polymer melt intercalation in organo-silicates for nanocomposites

The mean-field, lattice-based model of polymer melt intercalation in organically-modified layered silicates (OLS) originally developed by Vaia and Giannelis was applied for different polymers such as poly(methyl methacrylate) (PMMA), polypropylene (PP), and poly(ethylene oxide) (PEO). The nature of each polymer controls significantly the intercalation of the system. The internal energy change caused by the interaction of polymer, surfactant and clay is the strongest factor in determining the equilibrium structure of the nanocomposite system.     

Structure of an associating polymer melt in a narrow slit by molecular dynamics simulation

Molecular dynamics simulation has been used to study the equilibrium properties of a generic coarse-grained polymer melt with associating terminal groups, confined in a narrow slit by two atomically smooth walls. Simulations were carried out as a function of wall separation and attracting strength as well as polymer end-end interaction strength. We find that confinement has an important effect on the melt properties. In particular, strongly attracting walls can produce radical changes in chain conformation, the nature of the transient network, and the structure of the aggregates formed by the associating terminals.     

The directional contact distance of two ellipsoids: coarse-grained potentials for anisotropic interactions

We obtain the distance of closest approach of the surfaces of two arbitrary ellipsoids valid at any orientation and separation measured along their intercenter vector. This directional distance is derived from the elliptic contact function. The geometric meaning behind this approach is clarified. An elliptic pair potential for modeling arbitrary mixtures of elliptic particles, whether hard or soft, is proposed based on this distance. Comparisons with Gay-Berne potentials are discussed. Analytic expressions for the forces and torques acting on the elliptic particles are given.