Atomistic Simulations of Glassy Polystyrenes with Realistic Chain Conformations

Multiscale simulations of PS-SiO2 nanocomposites: from melt to glassy state

The interaction energetics, molecular packing, entanglement network properties, segmental dynamics, and elastic constants of atactic polystyrene-amorphous silica nanocomposites in the molten and the glassy state are studied via molecular simulations using two interconnected levels of representation: (a) a coarse-grained one, wherein each polystyrene repeat unit is mapped onto a single "superatom" and the silica nanoparticle is viewed as a solid sphere. Equilibration at all length scales at this level is achieved via connectivity-altering Monte Carlo simulations. (b) A united-atom (UA) level, wherein the polymer chains are represented in terms of a united-atom forcefield and the silica nanoparticle is represented in terms of a simplified, fully atomistic model. Initial configurations for UA molecular dynamics (MD) simulations are obtained by reverse mapping well-equilibrated coarse-grained configurations. By analysing microcanonical UA MD trajectories, the polymer density profile is studied and the polymer is found to exhibit layering in the vicinity of the nanoparticle surface. An estimate of the enthalpy of mixing between polymer and nanoparticles, derived from the UA simulations, compares favourably against available experimental values. The dynamical behaviour of polystyrene (in neat and filled melt systems) is characterized in terms of bond orientation and dihedral angle time autocorrelation functions. At low concentration in the molten polymer matrix, silica nanoparticles are found to cause a slight deceleration of the segmental dynamics close to their surface compared to the bulk polymer. Well-equilibrated coarse-grained long-chain configurations are reduced to entanglement networks via topological analysis with the CReTA algorithm, yielding a slightly lower density of entanglements in the filled than in the neat systems. UA melt configurations are glassified by MD cooling. The elastic moduli of the resulting glassy nanocomposites are computed through an analysis of strain fluctuations in the undeformed state and through explicit mechanical deformation by MD, showing a stiffening of the polymer in the presence of nanoparticles. UA simulation results for the elastic constants are compared to continuum micromechanical calculations invoked in homogenization models of the overall mechanical behaviour of heterogeneous materials. They can be interpreted in terms of the presence of an "interphase" of approximate thickness 2 nm around the nanoparticles, with elastic constants intermediate between those of the filler and the matrix.

Glass transition of polymers: atomistic simulation versus experiments

With experimental investigations and current theories, molecular modeling became an inevitable technique to study the perplexing phenomenon of glass transition. Among polymers, small variations in atomic interactions yield different values of the glass transition temperature, T{g}. To reveal the influence of differences in the atomic functionality on the value of T{g}, and thus to probe the molecular mechanisms responsible for this transition, atomistic simulations have to be undertaken. However, such simulations are argued not to accurately represent physically the glass transition due to the long relaxation times involved. Here we show the universality of the well-known Williams-Landel-Ferry equation for the experimental thermal dependence of polymer viscosities as demonstrated with atomistic simulations. Consequently, atomic aspects could be explicitly revealed. The contribution of atomistic simulation to the study of glass transition is thus confirmed. However, it has to be used complementarily with experiments and coarse-grained simulation to reveal the atomic aspects of current theories.

Molecular dynamics of polystyrene solutions in microwave fields

Equilibrium and nonequilibrium molecular dynamics simulation techniques were used to assess the influence of an applied microwave field on the dynamics of methylamine-methanol and methylamine-dimethylformamide (DMF) solutions bound within atactic polystyrene over a range of polymer densities from 35 to 96 wt % polymer. Atomistically detailed systems were studied, ranging from 3000 to 10 644 particles, using previously established potential models. Structural and dynamical properties were determined in the canonical (NVT) ensemble at 298 K. The simulated DMF self-diffusion coefficients in polystyrene solutions were compared with the zero-field experimental results established with pulsed-gradient spin-echo NMR spectrometry. A simulated external microwave field, with a rms electric field intensity of 0.1 VA, was applied to these systems and the simulated dynamical results over field frequencies up to 10(4) GHz were compared with the zero-field values. Simulated evidence of athermal effects on the diffusive characteristics of these mixtures is reported.