Strain hardening of polymer glasses: Effect of entanglement density, temperature, and rate

Dynamical heterogeneity in periodically deformed polymer glasses

The dynamics of structural relaxation in a model polymer glass subject to spatially homogeneous, time-periodic shear deformation is investigated using molecular dynamics simulations. We study a coarse-grained bead-spring model of short polymer chains below the glass transition temperature. It is found that at small strain amplitudes, the segmental dynamics is nearly reversible over about 10^{4} cycles, while at strain amplitudes above a few percent, polymer chains become fully relaxed after a hundred cycles. At the critical strain amplitude, the transition from slow to fast relaxation dynamics is associated with the largest number of dynamically correlated monomers as indicated by the peak value of the dynamical susceptibility. The analysis of individual monomer trajectories showed that mobile monomers tend to assist their neighbors to become mobile and aggregate into relatively compact transient clusters.

Simple model for chain packing and crystallization of soft colloidal polymers

We study a simple bead-spring polymer model exhibiting competing crystallization and glass transitions. Constant-pressure molecular dynamics simulations are employed to study phase behavior and morphological order. For adequately slow quench rates, chain systems exhibit a first-order phase transition (crystallization) below a critical temperature T=T(cryst). We observe the formation of close-packed crystallites of FCC and/or HCP order, separated by domain walls, twin defects, and amorphous regions. Such crystal structures closely resemble the corresponding ordered morphologies of athermal polymer packings: fully flexible chains retain random-walk-like configurations in the crystalline state and do not form lamellae, while semiflexible chains do form lamellae. The model presented here is well suited to the modeling of granular and colloidal polymers, in particular for elucidating the factors that dictate the formation of specific ordered morphologies.

Exploring the ability of a multiscale coarse-grained potential to describe the stress-strain response of glassy polystyrene

A new particle-based bottom-up method to develop coarse-grained models of polymers is presented and applied to polystyrene. The multiscale coarse-graining (MS-CG) technique of Izvekov et al. [J. Chem. Phys. 120, 10896 (2004)] is applied to a polymer system to calculate nonbonded interactions. The inverse Boltzmann inversion method was used to parametrize the bonded and bond-angle bending interactions. Molecular dynamics simulations were performed, and the CG model exhibited a significantly lower modulus compared to the atomistic model at low temperature and high strain rate. In an attempt to improve the CG model performance, several other parametrization schemes were used to build other models from this base model. The first of these models included standard frictional forces through use of the constant-temperature dissipative particle dynamics method that improved the modulus, yet was not transferrable to higher temperatures and lower strain rates. Other models were built by increasing the attraction between CG beads through direct manipulation of the nonbonded potential, where an improvement of the stress response was found. For these models, two parametrization protocols that shifted the force to more attractive values were explored. The first protocol involved a uniform shift, while the other protocol shifted the force in a more localized region. The uniformly shifted potential greatly affected the structure of the equilibrium model as compared to the locally shifted potential, yet was more transferrable to different temperatures and strain rates. Further improvements in the coarse-graining protocol to generate models that more satisfactorily capture mechanical properties are suggested.

Disentanglement of Linear Polymer Chains Toward Unentangled Crystals

Large scale and long time molecular dynamics simulations and primitive path analysis are used to investigate the disentanglement of long linear polymer chains during their crystallization from the melt state. In general, two competitive processes, a slow decrease of average entanglement length during cooling caused by stiffening of chains and a strong increase during crystallization, can be observed. In both homogeneous and heterogeneous nucleation, disentanglement occurs via forming folds from locally unentangled segments and continues in postcrystallization processes (slow reorganization), in particular, during annealing. Re-entanglement processes after melting are slow and can lead to memory effects in heating-recooling protocols such as self-seeding.

Strain stiffening in random packings of entangled granular chains

Random packings of granular chains are presented as a model system to investigate the contribution of entanglements to strain stiffening. The chain packings are sheared in uniaxial compression experiments. For short chain lengths, these packings yield when the shear stress exceeds the scale of the confining pressure, similar to granular packings of unconnected particles. In contrast, packings of chains which are long enough to form loops exhibit strain stiffening, in which the effective stiffness of the material increases with strain, similar to many polymer materials. The latter packings can sustain stresses orders-of-magnitude greater than the confining pressure, and do not yield until the chain links break. X-ray tomography measurements reveal that the strain-stiffening packings contain system-spanning clusters of entangled chains.

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.

Evolution of entanglements during the response to a uniaxial deformation of lamellar triblock copolymers and polymer glasses

Using coarse-grained molecular-dynamics simulations, a generic styrene-block-butadiene-block-styrene triblock copolymer under lamellar conformation is used in order to investigate the mutual entanglement evolution when a structure of alternating glassy (S)/rubbery (B) layers is submitted to an imposed deformation. By varying the amount of loop chains between each phase, i.e., noncrossing chains, it is possible to generate different types of S/B interface definitions. A specific boundary driven tensile strain protocol has been developed in order to mimic "real" experiments and measure the stress-strain curve. The same protocol is also applied to a reference state consisting in a directed glassy homopolymers, as well as to an isotropic glassy polymer. The evolution of initial mutual entanglements from the undeformed samples during the whole deformation process is monitored. It is shown for all considered systems that initial entanglements mostly participate to the preyield regime of the stress-strain curve and that this network is debonded during the strain-hardening regime. For triblocks with a non-null amount of crossing chains, the lower the amount is, the longer the memory effect of the initial entanglement network in the postyield regime is. On the fly distributions of entanglements, which depart from the postyield regime, depict memory effects and long-time correlations during the strain-hardening regime. For triblocks, loop chains reinforce these effects.

Microscopic view of accelerated dynamics in deformed polymer glasses

A molecular level analysis of segmental trajectories obtained from molecular dynamics simulations is used to obtain the full relaxation time spectrum in aging polymer glasses subject to three different deformation protocols. As in experiments, dynamics can be accelerated by several orders of magnitude, and a narrowing of the distribution of relaxation times during creep is directly observed. Additionally, the acceleration factor describing the transformation of the relaxation time distributions is computed and found to obey a universal dependence on the strain, independent of age and deformation protocol.

Coarse-graining scheme for simulating uniaxial stress-strain response of glassy polymers through molecular dynamics

Simulation of the deformation of polymers below their glass transition through molecular dynamics provides an useful route to correlate their molecular architecture to deformation behavior. However, present computational capabilities severely restrict the time and length scales that can be simulated when detailed models of these macromolecules are used. Coarse-graining techniques for macromolecular structures intend to make bigger and longer simulations possible by grouping atoms into superatoms and devising ways of determining reasonable force fields for the superatoms in a manner that retains essential macromolecular features relevant to the process under study but jettisons unnecessary details. In this work we systematically develop a coarse-graining scheme aimed at simulating uniaxial stress-strain behavior of polymers below their glass transition. The scheme involves a two step process of obtaining the coarse grained force field parameters above glass transition. This seems to be enough to obtain "faithful" stress-strain responses after quenching to below the glass transition temperature. We apply the scheme developed to a commercially important polymer polystyrene, derive its complete force field parameters and thus demonstrate the effectiveness of the technique.

Combined molecular algorithms for the generation, equilibration and topological analysis of entangled polymers: methodology and performance

We review the methodology, algorithmic implementation and performance characteristics of a hierarchical modeling scheme for the generation, equilibration and topological analysis of polymer systems at various levels of molecular description: from atomistic polyethylene samples to random packings of freely-jointed chains of tangent hard spheres of uniform size. Our analysis focuses on hitherto less discussed algorithmic details of the implementation of both, the Monte Carlo (MC) procedure for the system generation and equilibration, and a postprocessing step, where we identify the underlying topological structure of the simulated systems in the form of primitive paths. In order to demonstrate our arguments, we study how molecular length and packing density (volume fraction) affect the performance of the MC scheme built around chain-connectivity altering moves. In parallel, we quantify the effect of finite system size, of polydispersity, and of the definition of the number of entanglements (and related entanglement molecular weight) on the results about the primitive path network. Along these lines we approve main concepts which had been previously proposed in the literature.

Fracture in glassy polymers: a molecular modeling perspective

Over the past 25 years, molecular modeling and simulations have provided important insights into the physics of deformation and fracture of glassy polymers. This review presents an overview of key results discussed in the context of experimentally observed polymer behavior. Both atomistic and coarse-grained polymer models have been used in different deformation protocols to study elastic properties, shear yielding, creep, physical aging, strain hardening and crazing. Simulations reproduce most of the macroscopic features of plasticity in polymer glasses such as stress-strain relations and creep response, and reveal information about the underlying atomistic processes. Trends of the shear yield stress with loading conditions, temperature and strain rate, and the atomistic dynamics under load have been systematically explored. Most polymers undergo physical aging, which leads to a history-dependent mechanical response. Simulations of strain hardening and crazing demonstrate the nature of polymer entanglements in the glassy state and the role of local plasticity and provide insight into the origin of fracture toughness of amorphous polymers.

Suppressed segmental relaxation as the origin of strain hardening in polymer glasses

A nanometer scale dynamical theory is proposed for the large amplitude strain hardening phenomenon in polymer glasses. The new physical picture is that external deformation induces anisotropic chain conformations, which modifies interchain packing, resulting in density fluctuation suppression and intensification of localizing dynamical constraints and activation barriers. The resulting stresses are of intermolecular origin and arise primarily from prolongation of segmental relaxation, not single chain entropic rubber elasticity. Theoretical predictions for the magnitude, temperature, and deformation rate dependence of the hardening modulus are consistent with experiments and simulations.

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.

Simulations of aging and plastic deformation in polymer glasses

We study the effect of physical aging on the mechanical properties of a model polymer glass using molecular dynamics simulations. The creep compliance is determined simultaneously with the structural relaxation under a constant uniaxial load below yield at constant temperature. The model successfully captures universal features found experimentally in polymer glasses, including signatures of mechanical rejuvenation. We analyze microscopic relaxation time scales and show that they exhibit the same aging characteristics as the macroscopic creep compliance. In addition, our model indicates that the entire distribution of relaxation times scales identically with age. Despite large changes in mobility, we observe comparatively little structural change except for a weak logarithmic increase in the degree of short-range order that may be correlated with an observed decrease in aging with increasing load.