Robert S. Hoy

Topological analysis of polymeric melts: chain-length effects and fast-converging estimators for entanglement length.

Primitive path analyses of entanglements are performed over a wide range of chain lengths for both bead spring and atomistic polyethylene polymer melts. Estimators for the entanglement length N_{e} which operate on results for a single chain length N are shown to produce systematic O(1/N) errors. The mathematical roots of these errors are identified as (a) treating chain ends as entanglements and (b) neglecting non-Gaussian corrections to chain and primitive path dimensions. The prefactors for the O(1/N) errors may be large; in general their magnitude depends both on the polymer model and the method used to obtain primitive paths. We propose, derive, and test new estimators which eliminate these systematic errors using information obtainable from the variation in entanglement characteristics with chain length. The new estimators produce accurate results for N_{e} from marginally entangled systems. Formulas based on direct enumeration of entanglements appear to converge faster and are simpler to apply.

End grafted polymer nanoparticles in a polymeric matrix: Effect of coverage and curvature

It has recently been proposed that the miscibility of nanoparticles with a polymer matrix can be controlled by grafting polymer chains to the nanoparticle surface. We examine this hypothesis using molecular dynamics simulations on a single nanoparticle of radius R (4σ ≤ R ≤ 16σ, where σ is the diameter of a polymer monomer) grafted with chains of length 500 in a polymer melt of chains of length 1000. The grafting density Σ is varied between 0.04–0.32 chains/σ2. To facilitate equilibration a Monte Carlo double-bridging algorithm is applied - new bonds are formed across a pair of chains, creating two new chains each substantially different from the original. For the long brush chains studied here, the structure of the brush assumes its large particle limit even for R as small as 8σ, which is consistent with recent experimental findings and the small chain lattice simulations of Klos and Pakula. We study autophobic dewetting of the melt from the brush as a function of increasing Σ. Even these long brush and matrix chains of lengths 6 and 12 Ne, respectively, (the entanglement length is Ne ∼ 85) give somewhat ambiguous results for the interfacial width, showing that studies of two or more nanoparticles are necessary to properly understand these miscibility issues. Entanglement between the brush and melt chains is identified using path analysis. We find that the number of entanglements between the brush and melt chains scale simply with the product of the local monomer densities of brush and melt chains.

Strain hardening of polymer glasses: entanglements, energetics, and plasticity.

Simulations are used to examine the microscopic origins of strain hardening in polymer glasses. While stress-strain curves for a wide range of temperature can be fit to the functional form predicted by entropic network models, many other results are fundamentally inconsistent with the physical picture underlying these models. Stresses are too large to be entropic and have the wrong trend with temperature. The most dramatic hardening at large strains reflects increases in energy as chains are pulled taut between entanglements rather than a change in entropy. A weak entropic stress is only observed in shape recovery of deformed samples when heated above the glass transition. While short chains do not form an entangled network, they exhibit partial shape recovery, orientation, and strain hardening. Stresses for all chain lengths collapse when plotted against a microscopic measure of chain stretching rather than the macroscopic stretch. The thermal contribution to the stress is directly proportional to the rate of plasticity as measured by breaking and reforming of interchain bonds. These observations suggest that the correct microscopic theory of strain hardening should be based on glassy state physics rather than rubber elasticity.

Strain hardening in polymer glasses : Limitations of network models

Simulations are used to examine the microscopic origins of strain hardening in polymer glasses. While traditional entropic network models can be fit to the total stress, their underlying assumptions are inconsistent with simulation results. There is a substantial energetic contribution to the stress that rises rapidly as segments between entanglements are pulled taut. The thermal component of stress is less sensitive to entanglements, mostly irreversible, and directly related to the rate of local plastic rearrangements. Entangled and unentangled chains show the same strain hardening when plotted against the microscopic chain orientation rather than the macroscopic strain.

Structure of finite sphere packings via exact enumeration: implications for colloidal crystal nucleation.

We analyze the geometric structure and mechanical stability of a complete set of isostatic and hyperstatic sphere packings obtained via exact enumeration. The number of nonisomorphic isostatic packings grows exponentially with the number of spheres N, and their diversity of structure and symmetry increases with increasing N and decreases with increasing hyperstaticity H≡N_{c}-N_{ISO}, where N_{c} is the number of pair contacts and N_{ISO}=3N-6. Maximally contacting packings are in general neither the densest nor the most symmetric. Analyses of local structure show that the fraction f of nuclei with order compatible with the bulk (rhcp) crystal decreases sharply with increasing N due to a high propensity for stacking faults, five- and near-fivefold symmetric structures, and other motifs that preclude rhcp order. While f increases with increasing H, a significant fraction of hyperstatic nuclei for N as small as 11 retain non-rhcp structure. Classical theories of nucleation that consider only spherical nuclei, or only nuclei with the same ordering as the bulk crystal, cannot capture such effects. Our results provide an explanation for the failure of classical nucleation theory for hard-sphere systems of N≲10 particles; we argue that in this size regime, it is essential to consider nuclei of unconstrained geometry. Our results are also applicable to understanding kinetic arrest and jamming in systems that interact via hard-core-like repulsive and short-ranged attractive interactions.

Effect of equilibration on primitive path analyses of entangled polymers

We use recently developed primitive path analysis (PPA) methods to study the effect of equilibration on entanglement density in model polymeric systems. Values of Ne for two commonly used equilibration methods differ by a factor of 2-4 even though the methods produce similar large-scale chain statistics. We find that local chain stretching in poorly equilibrated samples increases entanglement density. The evolution of Ne with time shows that many entanglements are lost through fast processes such as chain retraction as the local stretching relaxes. Quenching a melt state into a glass has little effect on Ne. Equilibration-dependent differences in short-scale structure affect the craze extension ratio much less than expected from the differences in PPA values of Ne.

Fcc-bcc transition for Yukawa interactions determined by applied strain deformation

Calculations of the work required to transform between bcc and fcc phases yield a high-precision bcc-fcc transition line for monodisperse point Yukawa (screened-Coulomb) systems. Our results agree qualitatively but not quantitatively with recently published simulations and phenomenological criteria for the bcc-fcc transition. In particular, the bcc-fcc-fluid triple point lies at a higher inverse screening length than previously reported.

Particle-scale reversibility in athermal particulate media below jamming.

We perform numerical simulations of repulsive, frictionless athermal disks in two and three spatial dimensions undergoing cyclic quasistatic simple shear to investigate particle-scale reversible motion. We identify three classes of steady-state dynamics as a function of packing fraction φ and maximum strain amplitude per cycle γ(max). Point-reversible states, where particles do not collide and exactly retrace their intracycle trajectories, occur at low φ and γ(max). Particles in loop-reversible states undergo numerous collisions and execute complex trajectories but return to their initial positions at the end of each cycle. For sufficiently large φ and γ(max), systems display irreversible dynamics with nonzero self-diffusion. Loop-reversible dynamics enables the reliable preparation of configurations with specified structural and mechanical properties over a broad range of φ.

Viscoplasticity and large-scale chain relaxation in glassy-polymeric strain hardening

A simple theory for glassy-polymeric mechanical response that accounts for large-scale chain relaxation is presented. It captures the crossover from perfect-plastic response to Gaussian strain hardening as the degree of polymerization N increases, without invoking entanglements. By relating hardening to interactions on the scale of monomers and chain segments, we correctly predict its magnitude. Strain-activated relaxation arising from the need to maintain constant chain contour length reduces the characteristic relaxation time by a factor ~εN during active deformation at strain rate ε. This prediction is consistent with results from recent experiments and simulations, and we suggest how it may be further tested experimentally.

Effect of chain stiffness on the competition between crystallization and glass-formation in model unentangled polymers

We map out the solid-state morphologies formed by model soft-pearl-necklace polymers as a function of chain stiffness, spanning the range from fully flexible to rodlike chains. The ratio of Kuhn length to bead diameter (lK/r0) increases monotonically with increasing bending stiffness kb and yields a one-parameter model that relates chain shape to bulk morphology. In the flexible limit, monomers occupy the sites of close-packed crystallites while chains retain random-walk-like order. In the rodlike limit, nematic chain ordering typical of lamellar precursors coexists with close-packing. At intermediate values of bending stiffness, the competition between random-walk-like and nematic chain ordering produces glass-formation; the range of kb over which this occurs increases with the thermal cooling rate |Ṫ| implemented in our molecular dynamics simulations. Finally, values of kb between the glass-forming and rodlike ranges produce complex ordered phases such as close-packed spirals. Our results should provide a useful initial step in a coarse-grained modeling approach to systematically determining the effect of chain stiffness on the crystallization-vs-glass-formation competition in both synthetic and colloidal polymers.