Role of local structure on motions on the potential energy landscape for a model supercooled polymer.

We have conducted detailed Monte Carlo and molecular dynamics simulations of a model glass forming polymeric system near its apparent glass transition temperature. We have characterized the local structure of the glass using a Voronoi-Delaunay analysis of local particle arrangements. After a perturbative face elimination, we find that a significant fraction of Voronoi polyhedra consist of 12 pentagonal faces, a sign of icosahedral ordering. Further, we have identified metabasins of particle vibrations on the potential energy landscape on the basis of persistence of particle positions and neighbors over a simulated trajectory. We find that the residence times for vibrations are correlated with a particular Voronoi volume and number of neighbors of a particle; the largest metabasins correspond to particles whose average Voronoi volume is close to the value expected on the basis of the density, and whose approximate number of neighbors is close to 12. The local distortion around a particle, measured in terms of the tetrahedricity of the Delaunay simplices, reveals that the particles with a higher degree of local distortion are likely to transition faster to a neighboring metabasin. In addition to the transition between metabasins, we have also examined the influence of vibrations at inherent structures (IS) on the local structure, and find that the the low frequency modes at the IS exhibit the greatest curvature with respect to the local structure. We believe that these results establish an important connection between the local structure of glass formers and the activated dynamics, thereby providing insights into the origins of dynamic heterogeneities.

Local dynamic mechanical properties in model free-standing polymer thin films.

High-frequency sinusoidal oscillations of a coarse-grained polymer model are used to calculate the local dynamic mechanical properties (DMPs) of free-standing polymer thin films. The storage modulus G(') and loss modulus G(") are examined as a function of position normal to the free surfaces. It is found that mechanically soft layers arise near the free surfaces of glassy thin films, and that their thickness becomes comparable to the entire film thickness as the temperature approaches the glass transition T(g). As a result, the overall stiffness of glassy thin films decreases with film thickness. It is also shown that two regions coexist in thin films just at the bulk T(g); a melt-like region (G(')G(")) in the middle of the film. Our findings on the existence of a heterogeneous distribution of DMPs in free-standing polymer thin films provide insights into recent experimental measurements of the mechanical properties of glassy polymer thin films.

Mechanical heterogeneities in model polymer glasses at small length scales.

Molecular simulations of a model, deeply quenched polymeric glass show that the elastic moduli become strongly inhomogeneous at length scales comprising several tens of monomers; these calculations reveal a broad distribution of local moduli, with regions of negative moduli coexisting within a matrix of positive moduli. It is shown that local moduli have the same physical meaning as that traditionally ascribed to moduli obtained from direct measurements of local constitutive behaviors of macroscopic samples.

Influence of confinement on the vibrational density of states and the Boson peak in a polymer glass.

We have performed a normal-mode analysis on a glass forming polymer system for bulk and free-standing film geometries prepared under identical conditions. It is found that for free-standing film glasses, the normal-mode spectrum exhibits significant differences from the bulk glass with the appearance of an additional low-frequency peak and a higher intensity at the Boson peak frequency. A detailed eigenvector analysis shows that the low-frequency peak corresponds to a shear-horizontal mode which is predicted by continuum theory. The peak at higher frequency (Boson peak) corresponds to motions that are correlated over a length scale of approximately twice the interaction site diameter. These observations shed some light on the microscopic dynamics of glass formers, and help explain decreasing fragility that arises with decreasing thickness in thin films.

Density-of-states Monte Carlo simulation of a binary glass.

A newly proposed Monte Carlo formalism has been used to simulate a glass-forming liquid above and below the glass transition temperature. The heat capacity exhibits a sharp peak at a temperature lower than that reported from extensive molecular dynamics simulations. Its height is larger than that reported earlier. At temperatures below mode coupling, the average inherent-structure energy of the configurations generated in this work is significantly lower than that reported in the literature. The entropy of the supercooled liquid is calculated directly from our simulations and that of a disordered solid is calculated by a normal-mode analysis. We find that at low temperatures these two entropy curves become essentially parallel. They do not intersect each other, raising questions about the existence of a Kauzmann temperature in this glass-forming mixture.

Investigation of transition States in bulk and freestanding film polymer glasses.

We have performed transition state searches on the potential energy landscape for bulk and freestanding film polymer glasses and identified connected minima. An analysis of the displacements between minima shows that the sites that undergo the greatest displacement are highly localized in space for both the bulk and the thin-film systems studied. In the case of the thin film, the clusters originate at the surface and penetrate into the center of the film thereby coupling the relaxation in the center of the film to the mobile surface layer. Furthermore, the energy barriers between minima are lower in the thin film than in the bulk system. These findings can rationalize the experimentally observed depression of the glass transition temperature in freestanding polymer films.