Nonlinear Creep in a Polymer Glass

Local dynamics and failure of inhomogeneous polymer networks

Inhomogeneous crosslinked polymers are powerful platforms for materials design, because they can be synthesized from materials that provide complimentary properties to the resulting gel. For example, a membrane with both glassy and rubbery domains will be mechanically robust while enabling transport. The dynamics, and mechanical and failure properties of rubbery/glassy conetworks are only beginning to be studied, and there is likely to be strong heterogeneities in the dynamics and mechanical response. In this study, we use coarse-grained molecular dynamics simulations to generate microphase separated rubbery/glassy polymer networks with a bicontinuous morphology via in silico crosslinking. We study the effect of phase boundary on the local mobility gradient, and our simulation results reveal an asymmetric shift in the local mobility gradient across the interface that extends deeper into the phase with a lower Tg when the system temperature is between the glass transition temperatures of the two phases. Moreover, by employing a model that allows bond breaking, we examine the microscopic mechanism for failure in these networks as a function of the molecular weight of polymer strands between crosslinks and the number fraction of the glassy domain. Under uniaxial extension, we find that the stress is initially larger in the glassy domain. As the deformation proceeds, the segmental dynamics of the two phases homogenize, and subsequently bond breaking begins.

Structural Relaxation Time of a Polymer Glass during Deformation

In order to determine the structural relaxation time of a polymer glass during deformation, a strain rate switching experiment is performed in the steady-state plastic flow regime. A lightly cross-linked poly(methylmethacrylate) glass was utilized and, simultaneously, the segmental motion in the glass was quantified using an optical probe reorientation method. After the strain rate switch, a nonmonotonic stress response is observed, consistent with previous work. The correlation time for segmental motion, in contrast, monotonically evolves toward a new steady state, providing an unambiguous measurement of the structural relaxation time during deformation, which is found to be approximately equal to the segmental correlation time. The Chen-Schweizer model qualitatively predicts the changes in the segmental correlation time and the observed nonmonotonic stress response. In addition, our experiments are reasonably consistent with the material time assumption used in polymer deformation modeling; in this approach, the response of a polymer glass to a large deformation is described by combining a linear-response model with a time-dependent segmental correlation time.

Molecular Origin of the Reinforcement Effect and Its Strain-Rate Dependence in Polymer Nanocomposite Glass

We investigate the molecular origin of mechanical reinforcement in a polymer nanocomposite (PNC) under a glass state via molecular dynamics simulations. The strength of the PNC system is found to be reinforced mainly via reduced plastic deformations of the nanoparticle neighborhood (NN). Such a reinforcement effect is found to decay with an increase in the strain rate. The Arrhenius-Eyring relation is used to analyze its origin. The amplitude of the reinforcement is found to be determined by the difference between the energy barrier (ΔE) for the activation of NN and the work (W) done by the applied stress to conquer that barrier. A larger strain rate is found to result in a larger W and, hence, a weaker reinforcement effect. Such a strain-rate dependence is verified in the experimental tensile tests of a poly(vinyl alcohol)/SiO2 composite system. These results not only provide a new understanding of the molecular origin of the reinforcement effect in the PNC system, but also pave the way for a better design of the PNC material properties.

Creep and Recovery Behavior of Vitrimers with Fast Bond Exchange Rate

Vitrimers encompass the desirable mechanical properties of thermosets with the recyclability of thermoplastics. This ability arises from the rearrangement of the vitrimer covalent network upon heating via a bond shuffling mechanism while its cross-link density remains preserved. This unique feature makes vitrimers interesting candidates for the design of materials that combine dimensional stability at high temperatures and solvent resistance with the ability to be reshaped and processed. Despite these advantages, vitrimer exhibits significant creep at operating conditions where thermosets show little or no creep. As the mechanical properties of vitrimers not only depend on their chemical composition but also on the dynamics of the polymer chains, molecular dynamics (MD) simulations can provide detailed molecular mechanisms of the system of interest under macroscopic stress-induced deformations. In this regard, the recently developed MD/Monte Carlo simulation methodology capable of capturing the bond exchange mechanics in vitrimers is used to study the creep and recovery response of a coarse-grained model thermoset and vitrimer with a fast bond exchange rate. The time-stress superposition principle is then successfully applied to the creep response. The resulting universal curves enable us to predict the long-time creep behavior of both systems extending the timescale from 4 to over 10 orders of magnitude.

Role of Local Structure in the Enhanced Dynamics of Deformed Glasses

External stress can accelerate molecular mobility of amorphous solids by several orders of magnitude. The changes in mobility are commonly interpreted through the Eyring model, which invokes an empirical activation volume. Here, we analyze constant-stress molecular dynamics simulations and propose a structure-dependent Eyring model, connecting activation volume to a machine-learned field, softness. We show that stress has a heterogeneous effect on the mobility that depends on local structure through softness. The barrier impeding relaxation reduces more for well-packed particles, which explains the narrower distribution of relaxation time observed under stress.

Mechanical reinforcement of polymer colloidal crystals by supercritical fluids

Colloidal crystals realized by self-assembled polymer nanoparticles have prominent attraction as a platform for various applications from assembling photonic and phononic crystals, acoustic metamaterials to coating applications. However, the fragility of these systems limits their application horizon. In this work the uniform mechanical reinforcement and tunability of 3D polystyrene colloidal crystals by means of cold soldering are reported. This structural strengthening is achieved by high pressure gas (N₂ or Ar) plasticization at temperatures well below the glass transition. Brillouin light scattering is employed to monitor in-situ the mechanical vibrations of the crystal and thereby determine preferential pressure, temperature and time ranges for soldering, i.e. formation of physical bonding among the nanoparticles while maintaining the shape and translational order. This low-cost method is potentially useful for fabrication and tuning of durable devices including applications in photonics, phononics, acoustic metamaterials, optomechanics, surface coatings and nanolithography.

Predictive relation for the α-relaxation time of a coarse-grained polymer melt under steady shear

We examine the influence of steady shear on structural relaxation in a simulated coarse-grained unentangled polymer melt over a wide range of temperature and shear rates. Shear is found to progressively suppress the α-relaxation process observed in the intermediate scattering function, leading ultimately to a purely inertially dominated β-relaxation at high shear rates, a trend similar to increasing temperature. On the basis of a scaling argument emphasizing dynamic heterogeneity in cooled liquids and its alteration under material deformation, we deduce and validate a parameter-free scaling relation for both the structural relaxation time τα from the intermediate scattering function and the "stretching exponent" β quantifying the extent of dynamic heterogeneity over the entire range of temperatures and shear rates that we can simulate.

Molecular mobility in driven monomeric and polymeric glasses

We show that in monomeric supercooled liquids and glasses that are plastically flowing at a constant shear stress σ while being deformed with strain rate ε[over ̇], the microscopic structural relaxation time τ_{str} is given by the universal relation σ/G_{∞}ε[over ̇] with G_{∞} a modulus. This equality holds for all rheological regimes from temperatures above the glass transition all the way to the athermal limit, and arises from the competing effects of elastic loading and viscous dissipation. In macromolecular (polymeric) glasses, however, the stress decouples from this relaxation time and τ_{str} is in fact further reduced even though σ rises during glassy strain hardening. We develop expressions to capture both effects and thus provide a framework for analyzing mobility measurements in glassy materials.

Application of activated barrier hopping theory to viscoplastic modeling of glassy polymers

An established statistical mechanical theory of amorphous polymer deformation has been incorporated as a plastic mechanism into a constitutive model and applied to a range of polymer mechanical deformations. The temperature and rate dependence of the tensile yield of PVC, as reported in early studies, has been modeled to high levels of accuracy. Tensile experiments on PET reported here are analyzed similarly and good accuracy is also achieved. The frequently observed increase in the gradient of the plot of yield stress against logarithm of strain rate is an inherent feature of the constitutive model. The form of temperature dependence of the yield that is predicted by the model is found to give an accurate representation. The constitutive model is developed in two-dimensional form and implemented as a user-defined subroutine in the finite element package ABAQUS. This analysis is applied to the tensile experiments on PET, in some of which strain is localized in the form of shear bands and necks. These deformations are modeled with partial success, though adiabatic heating of the instability causes inaccuracies for this isothermal implementation of the model. The plastic mechanism has advantages over the Eyring process, is equally tractable, and presents no particular difficulties in implementation with finite elements.

Role of Dynamical Heterogeneities on the Mechanical Response of Confined Polymer

Confinement induces various modifications in the dynamics of polymers as compared to bulk. We focus here on the role of dynamical heterogeneities on the mechanics of confined polymers. Using a simple model that allows computation of the mechanical response over 10 decades in frequency, we show that the local mechanical coupling controlling the macroscopic response in the bulk disappears in a confined geometry. The slowest domains significantly contribute to the mechanical response for increasing confinement. As a consequence, the apparent glass transition is broadened and shifted towards lower frequencies as confinement increases. We compare our numerical predictions with experiments performed on poly(ethylacrylate) chains in model filled elastomers. We suggest that the change of elastic coupling between domains induced by confinement should contribute significantly to the polymer mobility shift observed on filled systems.

Inherent structure energy is a good indicator of molecular mobility in glasses

Glasses produced via physical vapor deposition can display greater kinetic stability and lower enthalpy than glasses prepared by liquid cooling. While the reduced enthalpy has often been used as a measure of the stability, it is not obvious whether dynamic measures of stability provide the same view. Here, we study dynamics in vapor-deposited and liquid-cooled glass films using molecular simulations of a bead-spring polymer model as well as a Lennard-Jones binary mixture in two and three dimensions. We confirm that the dynamics in vapor-deposited glasses is indeed slower than in ordinary glasses. We further show that the inherent structure energy is a good reporter of local dynamics, and that aged systems and glasses prepared by cooling at progressively slower rates exhibit the same behavior as vapor-deposited materials when they both have the same inherent structure energy. These findings suggest that the stability inferred from measurements of the energy is also manifested in dynamic observables, and they strengthen the view that vapor deposition processes provide an effective strategy for creation of stable glasses.

Simulations of inorganic-bioorganic interfaces to discover new materials: insights, comparisons to experiment, challenges, and opportunities

Natural and man-made materials often rely on functional interfaces between inorganic and organic compounds. Examples include skeletal tissues and biominerals, drug delivery systems, catalysts, sensors, separation media, energy conversion devices, and polymer nanocomposites. Current laboratory techniques are limited to monitor and manipulate assembly on the 1 to 100 nm scale, time-consuming, and costly. Computational methods have become increasingly reliable to understand materials assembly and performance. This review explores the merit of simulations in comparison to experiment at the 1 to 100 nm scale, including connections to smaller length scales of quantum mechanics and larger length scales of coarse-grain models. First, current simulation methods, advances in the understanding of chemical bonding, in the development of force fields, and in the development of chemically realistic models are described. Then, the recognition mechanisms of biomolecules on nanostructured metals, semimetals, oxides, phosphates, carbonates, sulfides, and other inorganic materials are explained, including extensive comparisons between modeling and laboratory measurements. Depending on the substrate, the role of soft epitaxial binding mechanisms, ion pairing, hydrogen bonds, hydrophobic interactions, and conformation effects is described. Applications of the knowledge from simulation to predict binding of ligands and drug molecules to the inorganic surfaces, crystal growth and shape development, catalyst performance, as well as electrical properties at interfaces are examined. The quality of estimates from molecular dynamics and Monte Carlo simulations is validated in comparison to measurements and design rules described where available. The review further describes applications of simulation methods to polymer composite materials, surface modification of nanofillers, and interfacial interactions in building materials. The complexity of functional multiphase materials creates opportunities to further develop accurate force fields, including reactive force fields, and chemically realistic surface models, to enable materials discovery at a million times lower computational cost compared to quantum mechanical methods. The impact of modeling and simulation could further be increased by the advancement of a uniform simulation platform for organic and inorganic compounds across the periodic table and new simulation methods to evaluate system performance in silico.

Mapping Brittle and Ductile Behaviors of Polymeric Glasses under Large Extension

We have carried out a series of tensile extension tests on the two most common polymer glasses to describe their generic mechanical responses as a function of deformation rate at different temperatures. The essentially defect-free polystyrene and poly(methyl methacrylate) both show remarkable re-entrant failure: being ductile at intermediate rates and showing diminishing toughness at both higher and lower rates. We draw phase diagrams to map out the relationship between brittle-like and yield-like states in terms of temperature, rate, and stress. A coherent understanding of the rich phenomenology requires us to describe in more detail the interplay between the chain network and the primary structure bonded by intersegmental van der Waals forces.

Universal rescaling of flow curves for yield-stress fluids close to jamming

The experimental flow curves of four different yield-stress fluids with different interparticle interactions are studied near the jamming concentration. By appropriate scaling with the distance to jamming all rheology data can be collapsed onto master curves below and above jamming that meet in the shear-thinning regime and satisfy the Herschel-Bulkley and Cross equations, respectively. In spite of differing interactions in the different systems, master curves characterized by universal scaling exponents are found for the four systems. A two-state microscopic theory of heterogeneous dynamics is presented to rationalize the observed transition from Herschel-Bulkley to Cross behavior and to connect the rheological exponents to microscopic exponents for the divergence of the length and time scales of the heterogeneous dynamics. The experimental data and the microscopic theory are compared with much of the available literature data for yield-stress systems.

Stability of polymer glasses vitrified under stress

How stress or strain imparts mobility to glasses is a scientific issue linking ideas of jamming and the glass transition across colloids, granular materials, polymers, and molecular glasses. Here, we address for the first time how stress applied during vitrification, formation of the glassy state by a temperature quench, affects the subsequent stability of the glassy state, even after the stress has been removed. Using entangled polymers that are easily manipulated mechanically above the glass transition temperature, we find that the resulting polymer glasses become less stable, exhibiting a higher physical aging rate, when stress is applied while rapidly cooling the polymer films. The data show an initial plateau value at low stress, before transitioning rapidly to a higher aging rate at larger stress. These results are suggestive of the glassy system being left trapped in a less stable, higher energy state with faster physical aging rate when stressed above some minimum value during vitrification.

Evolution of collective motion in a model glass-forming liquid during physical aging

At temperatures moderately below their glass transition temperature, the properties of many glass-forming materials can evolve slowly with time in a process known as physical aging whereby the thermodynamic, mechanical, and dynamic properties all drift towards their equilibrium values. In this work, we study the evolution of the thermodynamic and dynamic properties during physical aging for a model polymer glass. Specifically, we test the relationship between an estimate of the size of the cooperative rearrangements taking the form of strings and the effective structural relaxation time predicted by the Adam-Gibbs relationship for both an equilibrium supercooled liquid and the same fluid undergoing physical aging towards equilibrium after a series of temperature jumps. We find that there is apparently a close correlation between a structural feature of the fluid, the size of the string-like rearrangements, and the structural relaxation time, although the relationship for the aging fluid appears to be distinct from that of the fluid at equilibrium.

Mechanical properties of thin confined polymer films close to the glass transition in the linear regime of deformation: theory and simulations

Over the past twenty years experiments performed on thin polymer films deposited on substrates have shown that the glass transition temperature T(g) can either decrease or increase depending on the strength of the interactions. Over the same period, experiments have also demonstrated that the dynamics in liquids close to the glass transition temperature is strongly heterogeneous, on the scale of a few nanometers. A model for the dynamics of non-polar polymers, based on percolation of slow subunits, has been proposed and developed over the past ten years. It proposes a unified mechanism regarding these two features. By extending this model, we have developed a 3D model, solved by numerical simulations, in order to describe and calculate the mechanical properties of polymers close to the glass transition in the linear regime of deformation, with a spatial resolution corresponding to the subunit size. We focus on the case of polymers confined between two substrates with non-negligible interactions between the polymer and the substrates, a situation which may be compared to filled elastomers. We calculate the evolution of the elastic modulus as a function of temperature, for different film thicknesses and polymer-substrate interactions. In particular, this allows to calculate the corresponding increase of glass transition temperature, up to 20 K in the considered situations. Moreover, between the bulk T(g) and T(g) + 50 K the modulus of the confined layers is found to decrease very slowly in some cases, with moduli more than ten times larger than that of the pure matrix at temperatures up to T(g) + 50 K. This is consistent with what is observed in reinforced elastomers. This slow decrease of the modulus is accompanied by huge fluctuations of the stress at the scale of a few tens of nanometers that may even be negative as compared to the solicitation, in a way that may be analogous to mechanical heterogeneities observed recently in molecular dynamics simulations. As a consequence, confinement may result not only in an increase of the glass transition temperature, but in a huge broadening of the glass transition.

Glassy dynamics and mechanical response in dense fluids of soft repulsive spheres. II. Shear modulus, relaxation-elasticity connections, and rheology

We apply the quiescent and mechanically driven versions of nonlinear Langevin equation theory to study how particle softness influences the shear modulus, the connection between shear elasticity and activated relaxation, and nonlinear rheology of the repulsive Hertzian contact model of dense soft sphere fluids. Below the soft jamming threshold, the shear modulus follows a power law dependence on volume fraction over a narrow interval with an apparent exponent that grows with particle stiffness. To a first approximation, the elastic modulus and transient localization length are controlled by a single coupling constant determined by local fluid structure. In contrast to the behavior of hard spheres, an approximately linear relation between the shear modulus and activation barrier is predicted. This connection has recently been observed for microgel suspensions and provides a microscopic realization of the elastic shoving model. Yielding, shear and stress thinning of the alpha relaxation time and viscosity, and flow curves are also studied. Yield strains are relatively weakly dependent on volume fraction and particle stiffness. Shear thinning commences at values of the effective Peclet number far less than unity, a signature of stress-assisted activated relaxation when barriers are high. Apparent power law reduction of the viscosity with shear rate is predicted with a thinning exponent less than unity. In the vicinity of the soft jamming threshold, a power law flow curve occurs over an intermediate reduced shear rate range with an apparent exponent that decreases as fluid volume fraction and/or repulsion strength increase.