Molecular dynamics study of polymer crystallization in the presence of a particle

We have used molecular dynamics simulations with a coarse‐grained model to study the effect of a particle on the crystallization of polymer melt. We analyzed in particular a bond order parameter to characterize the nucleation and crystallization process. Our calculations show that the presence of a particle modifies the free energy landscape of polymer melts, locally induces the ordering of polymer melts near the particle surface, and thus enhances the polymer crystallization. Because the interaction between the particle and polymers is repulsive, our results suggest that the origin of the enhancement for polymer crystallization is entropic. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 2161–2166, 2007     

Adsorption of a single polymer chain on a surface: A molecular dynamics simulation study

We performed simulations of the physical adsorption of a single globular chain on a surface of hemispherical shape by means of molecular dynamics simulations. For the chain, we took advantage of a united atom model. Interactions within the chain were limited to stretching, bending, and torsional as well as nonbonded interactions between the nonadjacent atoms. The interaction between each chain element and the surface formation are reigned by a Lennard–Jones potential. In this article, we focused on differences in the behavior of the adsorbed globule to the free unadsorbed one particularly in two different zones of the immediate vicinity of the surface. There were strong indications for a localized acceleration of the dynamics as compared with the bulk that appears in an increase of trans–gauche switches. For explanation we came up with an adsorption scenario. Special attention was given to the shift of the percentage of trans and gauche conformations within the globule in dependence on the strength of the adsorption potential that might be related to crystallization or glass transition. © 2001 John Wiley & Sons, Inc. J Polym Sci Part B: Polym Phys 39: 2333–2339, 2001     

A steered molecular dynamics simulation on the elastic behavior of adsorbed star polymer chains

<正>Elastic behavior of 4-branched star polymer chain with different chain length N adsorbed on attractive surface is investigated using steered molecular dynamics(SMD) simulation method based on the united-atom(UA) model for branched alkanes.The simulation is realized by pulling up the chain via a linear spring with a constant velocity v = 0.005 nm/ps.At the beginning,the chain lies extensionally on adsorbed surface and suffers continuous deformations during the tensile process.Statistical parameters as mean-square radii of gyration S~2_(xy),S~2_z,shape factor δ,describing the conformational changes,sectional density den which gives the states of the chain,and average surface attractive energy U_a,average total energy U,average force f probed by the spring,which characterize the thermodynamic properties, are calculated in the stimulant process.Remarkably,distinguishing from the case in linear chains that there only exists one long plateau in the curve of f,the force plateau in our study for star chains is multiple,denoting different steps of desorption,and this agrees well with the experimental results in essence.We find during the tensile process,there are three characteristic distances Z_c,Z_t and Z_0 from the attractive surface,and these values vary with N.When Z=Z_c,the chain is stripped from the surface,but due to the form of wall-monomer interaction,the surface retains weak influence on the chain till Z = Z_c.From Z=Z_t,parameters U_a,U and f respectively reach a stable value,while the shape and the size of the chain still need adjustments after Z_t till Z_0 to reach their equilibrium states.Specifically,for short chain of N= 41,Z_t and Z_0 are incorporated.These results may help us to deepen the knowledge about the elastic behavior of adsorbed star polymer chains.     

A study on the structural transition of a single polymer chain by parallel tempering molecular dynamics simulation

The structural transition of a single polymer chain with chain length of 100,200 and 300 beads was investigated by parallel tempering MD simulation.Our simulation results can capture the structural change from random coil to orientationally ordered structure with decreasing temperature.The clear transition was observed on the curves of radius of gyration and global orientational order parameter P as the function of temperature,which demonstrated structural formation of a single polymer chain.The linear relationships between three components of square radius of gyration R_gx²,R_gx²,R_gz² and global orientational order P can be obtained under the structurally transformational process.The slope of the linear relationship between x（or y-axis） component R_gx²(or R_gy²) and P is negative,while that of RL as the function of P is positive.The absolute value of slope is proportional to the chain length.Once the single polymer chain takes the random coil or ordered configuration,the linear relationship is invalid.The conformational change was also analyzed on microscopic scale.The polymer chain can be treated as the construction of rigid stems connecting by flexible loops.The deviation from exponentially decreased behavior of stem length distribution becomes prominent,indicating a stiffening of the chain arises leading to more and more segments ending up in the trans state with decreasing temperature.The stem length N_tr is about 21 bonds indicating the polymer chain is ordered with the specific fold length.So,the simulation results,which show the prototype of a liquid-crystalline polymer chain,are helpful to understand the crystallization process of crystalline polymers.     

Schmidt number effects in dissipative particle dynamics simulation of polymers

Simulation studies for dilute polymeric systems are presented using the dissipative particle dynamics method. By employing two different thermostats, the velocity-Verlet and Lowe's scheme, we show that the Schmidt number (S(c)) of the solvent strongly affects nonequilibrium polymeric quantities. The fractional extension of wormlike chains subjected to steady shear is obtained as a function of S(c). Poiseuille flow in microchannels for fixed polymer concentration and varying number of repeated units within a chain is simulated. The nonuniform concentration profiles and their dependence on S(c) are computed. We show the effect of the bounce-forward wall boundary condition on the depletion layer thickness. A power law fit of the velocity profile in stratified Poiseuille flow in a microchannel yields wall viscosities different from bulk values derived from uniform, steady plane Couette flow. The form of the velocity profiles indicates that the slip flow model is not useful for the conditions of these calculations.     

Nanospheres in phase-separating multicomponent fluids: a three-dimensional dissipative particle dynamics simulation

The dynamics of phase separation of three-dimensional fluids containing nanospheres, which interact preferentially with one of the two fluids, is studied by means of large-scale dissipative particle dynamics simulations. We systematically investigated the effect of volume fraction, radius, and mass of the nanoparticles on both kinetics and morphology of the binary mixture. We found that nanospheres lead to a reduction of domain growth which is intensified as their volume fraction is increased for a given radius of nanoparticles, or as the nanoparticles radius is decreased for a given volume fraction. Up to moderate volume fractions of nanoparticles, the growth law, however, is found to be identical to that pure binary fluids, i.e., R(t) approximately t(n), with n=1. For relatively high volume fractions of nanoparticles, a diffusive growth regime was detected. The crossover to the slower growth regime as the nanoparticles volume fraction is increased or their radius is decreased is associated with the crystallization of the nanospheres within the preferred component. These results are qualitatively in good agreement with previous two-dimensional simulations using molecular dynamics [M. Laradji and G. MacNevin, J. Chem. Phys. 119, 2275 (2003)] and a time-dependent Ginzburg-Landau model [M. Laradji, J. Chem. Phys. 120, 9330 (2004)], as well as recent experiments.     

Atomistic hybrid particle‐field molecular dynamics combined with slip‐springs: Restoring entangled dynamics to simulations of polymer melts

In hybrid particle‐field (hPF) simulations (J. Chem. Phys., 2009 130, 214106), the entangled dynamics of polymer melts is lost due to chain crossability. Chains cross, because the field‐treatment of the nonbonded interactions makes them effectively soft‐core. We introduce a multi‐chain slip‐spring model (J. Chem. Phys., 2013 138, 104907) into the hPF scheme to mimic the topological constraints of entanglements. The structure of the polymer chains is consistent with that of regular molecular dynamics simulations and is not affected by the introduction of slip‐springs. Although slight deviations are seen at short times, dynamical properties such as mean‐square displacements and reorientational relaxation times are in good agreement with traditional molecular dynamics simulations and theoretical predictions at long times.     

Mesoscopic simulation of self-assembly in surfactant oligomers by dissipative particle dynamics

Self-assembling properties of surfactant oligomers in an aqueous medium is simulated by dissipative particle dynamics (DPD). The critical micellar concentration (CMC) of dimeric (oligomerization = 2) and trimeric (oligomerization = 3) surfactant is much lower than their single-chain counterpart. All surfactants form spherical micelles at the concentration not far above their CMC. The transition from spherical to cylindrical micelles exhibits with increasing surfactant concentration. Lamellar micelles will appear with further increasing the surfactant concentration. For dimeric and trimeric surfactants, cylindrical micelles transform into extremely long “wormlike” or “threadlike” micelles before the transition to lamellar micelles. These results are in qualitative agreement with laboratory experiment. Average aggregation numbers (AN) of micelles increase with a power law of AN  c when the surfactant concentration c CMC. The self-diffusion coefficients will drop with a power law of D  c⁻ when wormlike micelles are formed.     

Influence of polymer matrix composition and architecture on polymer nanocomposite formation: Coarse-grained molecular dynamics simulation

Coarse-grained molecular dynamics simulations of stacks of two-dimensional platelets immersed in a polymer melt were performed to investigate aspects of the polymer matrix that promote the formation of intercalated or exfoliated nanocomposite structures. Such factors include temperature, copolymer architecture, and blend composition. Increasing the polymer-sheet attractive interaction led to binding of the sheets, where individual beads simultaneously attract two neighboring sheets, thus kinetically blocking intercalation by occupying the perimeter of the affected gallery. Polymers with a small polymer-sheet attraction, but having a strongly attractive chain end (end-functionalized polymers) minimized the bonding of adjacent sheets. These systems exhibited some sheet sliding because a majority of the confined polymer beads only interacted weakly with adjacent sheets; however, the number density of intercalated polymer was low. Mixtures of end-functionalized and nonfunctionalized polymers, however, yielded better intercalation efficiency. For the mixed system, the reduced number of highly attractive beads provided sufficient interaction for intercalation to occur, enabling greater intercalation rates, less sheet-bridging, and incorporation of the nonfunctionalized polymers into the galleries. © 2003 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 41: 3272–3284, 2003     

Hairpin polymer unfolding in square nanochannels

The unfolding dynamics of a flexible hairpin polymer inserted in a square nanochannel is studied using Brownian dynamics simulations of the bead‐spring model. Because the hairpin polymer is not an equilibrium configuration, the molecule starts unfolding until it reaches a stretched configuration inside the tube. We study the effect of varying the channel height and width D, and the number of monomers N in the folded arm on the unfolding times. We show that for square nanochannels, the unfolding time scales as DN2, for small values of D. The unfolding relaxation dynamics obeys similar mechanisms described in the escaping dynamics of partially inserted polymers in cylindrical nanotubes. We also show that the velocity of the polymer center of mass scales as D^?1, in agreement with DNA unfolding experiments in solid‐state nanochannels and recent computational simulations. © 2013 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2013, 51, 1411–1418     

Bead-spring models of entangled polymer melts: Comparison of hard-core and soft-core potentials

Two bead-spring models of flexible chains for generic coarse graining of entangled polymer melts, the excluded volume Kremer–Grest (KG) model and the modified segmental repulsive potential (mSRP) combined with a weakly repulsive potential, are compared. For chains containing an equivalent number of entanglements, we compare the chain characteristics of the KG and mSRP polymer models by determining the ratios of the entanglement lengths , the required total number of particles to capture comparable entanglement phenomena , and the time scaling ratios τ^mSRP/τ^KG. Our findings show that systems using the mSRP polymer model require half the number of particles and relax four times faster compared to the KG polymer model. © 2012 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys, 2012     

Multiscale Lagrangian fluid dynamics simulation for polymeric fluid

We have developed a simulation technique of multiscale Lagrangian fluid dynamics to tackle hierarchical problems relating to historical dependency of polymeric fluid. We investigate flow dynamics of dilute polymeric fluid by using the multiscale simulation approach incorporating Lagrangian particle fluid dynamics technique (the modified smoothed particle hydrodynamics) with stochastic coarse‐grained polymer simulators (the dumbbell model). We have confirmed that our approach is well in agreement with the macroscopic results obtained by a constitutive equation corresponding to the dumbbell model, and observed that microscopic thermal fluctuation appears in macroscopic fluid dynamics as dispersion phenomena. © 2010 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 48: 886–893, 2010     

Flow-induced translocation of polymers through a fluidic channel: a dissipative particle dynamics simulation study

The dynamics of flow-induced translocation of polymers through a fluidic channel has been studied by dissipative particle dynamics (DPD) approach. Unlike implicit solvent models, the many-body energetic and hydrodynamic interactions are preserved naturally by incorporating explicit solvent particles in this approach. The no-slip wall boundary and the adaptive boundary conditions have been implemented in the modified DPD approach to model the hydrodynamic flow within a specific wall structure of fluidic channel and control the particles' density fluctuations. The results show that the average translocation time versus polymer chain length satisfies a power-law scaling of τ ～N(1.152). The conformational changes and translocation dynamics of polymers through the fluidic channel have also been investigated in our simulations, and two different translocation processes, i.e., the single-file and double-folded translocation events, have been observed in detail. These findings may be helpful in understanding the conformational and dynamic behaviors of such polymer and/or DNA molecules during the translocation processes.     

Mesoscale simulation of polymer reaction equilibrium: combining dissipative particle dynamics with reaction ensemble Monte Carlo. I. Polydispersed polymer systems

We present a mesoscale simulation technique, called the reaction ensemble dissipative particle dynamics (RxDPD) method, for studying reaction equilibrium of polymer systems. The RxDPD method combines elements of dissipative particle dynamics (DPD) and reaction ensemble Monte Carlo (RxMC), allowing for the determination of both static and dynamical properties of a polymer system. The RxDPD method is demonstrated by considering several simple polydispersed homopolymer systems. RxDPD can be used to predict the polydispersity due to various effects, including solvents, additives, temperature, pressure, shear, and confinement. Extensions of the method to other polymer systems are straightforward, including grafted, cross-linked polymers, and block copolymers. To simulate polydispersity, the system contains full polymer chains and a single fractional polymer chain, i.e., a polymer chain with a single fractional DPD particle. The fractional particle is coupled to the system via a coupling parameter that varies between zero (no interaction between the fractional particle and the other particles in the system) and one (full interaction between the fractional particle and the other particles in the system). The time evolution of the system is governed by the DPD equations of motion, accompanied by changes in the coupling parameter. The coupling-parameter changes are either accepted with a probability derived from the grand canonical partition function or governed by an equation of motion derived from the extended Lagrangian. The coupling-parameter changes mimic forward and reverse reaction steps, as in RxMC simulations.     

Shape instabilities in adsorbed polymer condensates

A polymer chain in a poor solvent collapses to a spherical globule. If this globule is subsequently deformed, either by stretching it or if it has some non-zero net electric charge, it will become unstable to sinusoidal perturbations leading to novel final states. When such globules are imaged by direct methods such as atomic force microscopy the substrate on which they adsorb can affect the final morphological state. We therefore investigate strongly adsorbed polymers in poor solvents. We demonstrate that the real-space, Self-Consistent Field method is an ideal numerical tool in predicting equilibrium morphologies. New structures are predicted, which support previous explicit free energy calculations. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 3327–3337, 2007     

Effect of particle agglomeration and interphase on the glass transition temperature of polymer nanocomposites

In this article, we utilize finite element modeling to investigate the effect of nanoparticle agglomeration on the glass transition temperature of polymer nanocomposites. The case of an attractive interaction between polymer and nanofiller is considered for which an interphase domain of gradient properties is developed. This model utilizes representative volume elements that are created and analyzed with varying degrees of nanoparticle clustering and length scale of interphase domain. The viscoelastic properties of the composites are studied using a statistical approach to account for variations due to the random nature of the microstructure. Results show that a monotonic increase in nanofiller clustering not only results in the loss of interphase volume but also obstructs the formation of a percolating interphase network in the nanocomposite. The combined impacts lead to a remarkable decrease of T_g enhancement of clustering nanofillers in comparison with a well‐dispersed configuration. Our simulation results provide qualitative support for experimental observations that clustering observed at high nanofiller concentrations negatively impacts the effects of the nanofiller on overall properties. © 2011 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys, 2011     

Comparison between theoretical values and simulation results of viscosity for the dissipative particle dynamics method

In order to investigate the validity of the dissipative particle dynamics method, which is a mesoscopic simulation technique, we have derived an expression for viscosity from the equation of motion of dissipative particles. In the concrete, we have shown the Fokker-Planck equation in phase space, and macroscopic conservation equations such as the equation of continuity and the equation of momentum conservation. The basic equations of the single-particle and pair distribution functions have been derived using the Fokker-Planck equation. The solutions of these distribution functions have approximately been solved by the perturbation method under the assumption of molecular chaos. The expressions of the viscosity due to momentum and dissipative forces have been obtained using the approximate solutions of the distribution functions. Also, we have conducted nonequilibrium dynamics simulations to investigate the influence of the parameters, which have appeared in defining the equation of motion in the dissipative particle dynamics method. The theoretical values of the viscosity due to dissipative forces in the Hoogerbrugge-Koelman theory are in good agreement with the simulation results obtained by the nonequilibrium dynamics method, except in the range of small number densities. There are restriction conditions for taking appropriate values of the number density, number of particles, time interval, shear rate, etc., to obtain physically reasonable results by means of dissipative particle dynamics simulations.     

Molecular dynamics simulation on the interaction between polymer inhibitors and anhydrite surface

Investigation on the microscopic interaction between polymer inhibitors and calcium sulfate will be helpful for understanding its scale inhibition mechanism and can provide a theoretical guidance to developing new scale inhibitors. In this work, molecular dynamics simulations with COMPASS force field have been performed to simulate the interaction between hydrolyzed polymaleic anhydride (HPMA), polyaspartic acid (PASP), polyepoxysuccinic acid (PESA), polyacrylic acid (PAA) and the (001) and (020) surfaces of anhydrite (AD) crystal with and without water. The results show that the sequence of binding energies between four polymer inhibitors and AD (001) and (020) with water is PESA > PASP > HPMA > PAA. The binding energy of the same polymer inhibitor on AD (001) is smaller than that on AD (020). Water molecules weaken the deformations of HPMA and PAA but aggravate those of PASP and PESA. Natural bond orbital (NBO) charges of the repeat units of polymer inhibitors were calculated by B3LYP/6‐31G* method. The Coulomb interaction is formed between the O atoms of polymer inhibitors and the Ca atoms of AD crystal. The system of polymer–AD is mainly contributed from the non‐bonding interaction. Polymer inhibitors do not interact directly with AD crystal, but indirectly through the interactions between inhibitor–H₂O and H₂O–AD, i.e. water molecules participate in scale inhibition of polymer inhibitors to AD crystal. Water molecules cannot be ignored when the interaction models are constructed, i.e. solvent effect cannot be ignored. Copyright © 2015 John Wiley & Sons, Ltd.     

Polymer grafted nanoparticles: Effect of chemical and physical heterogeneity in polymer grafts on particle assembly and dispersion

Macroscopic properties of polymer nanocomposites depend on the microscopic composite morphology of the constituent nanoparticles and polymer matrix. One way to control the spatial arrangement of the nanoparticles in the polymer matrix is by grafting the nanoparticle surfaces with polymers that can tune the effective interparticle interactions in the polymer matrix. A fundamental understanding of how graft and matrix polymer chemistries and molecular weight, grafting density, and nanoparticle size, and chemistry affect interparticle interactions is needed to design the appropriate polymer ligands to achieve the target morphology. Theory and simulations have proven to be useful tools in this regard due to their ability to link molecular level interactions to the morphology. In this feature article, we present our recent theory and simulation studies of polymer grafted nanoparticles with chemical and physical heterogeneity in grafts to calculate the effective interactions and morphology as a function of chemistry, molecular weights, grafting densities, and so forth. © 2013 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys, 2013