Strain localization in glassy polymers under cylindrical confinement

Identifying microscopic factors that influence ductility in disordered solids

There are empirical strategies for tuning the degree of strain localization in disordered solids, but they are system-specific and no theoretical framework explains their effectiveness or limitations. Here, we study three model disordered solids: a simulated atomic glass, an experimental granular packing, and a simulated polymer glass. We tune each system using a different strategy to exhibit two different degrees of strain localization. In tandem, we construct structuro-elastoplastic (StEP) models, which reduce descriptions of the systems to a few microscopic features that control strain localization, using a machine learning-based descriptor, softness, to represent the stability of the disordered local structure. The models are based on calculated correlations of softness and rearrangements. Without additional parameters, the models exhibit semiquantitative agreement with observed stress-strain curves and softness statistics for all systems studied. Moreover, the StEP models reveal that initial structure, the near-field effect of rearrangements on local structure, and rearrangement size, respectively, are responsible for the changes in ductility observed in the three systems. Thus, StEP models provide microscopic understanding of how strain localization depends on the interplay of structure, plasticity, and elasticity.

The initiation of shear band formation in deformed metallic glasses from soft localized domains

It has long been thought that shear band (SB) formation in amorphous solids initiates from relatively "soft" regions in the material in which large-scale non-affine deformations become localized. The test of this hypothesis requires an effective means of identifying "soft" regions and their evolution as the material is deformed to varying degrees, where the metric of "softness" must also account for the effect of temperature on local material stiffness. We show that the mean square atomic displacement on a caging timescale ⟨u²⟩, the "Debye-Waller factor," provides a useful method for estimating the shear modulus of the entire material and, by extension, the material stiffness at an atomic scale. Based on this "softness" metrology, we observe that SB formation indeed occurs through the strain-induced formation of localized soft regions in our deformed metallic glass free-standing films. Unexpectedly, the critical strain condition for SB formation occurs when the softness (⟨u²⟩) distribution within the emerging soft regions approaches that of the interfacial region in its undeformed state, initiating an instability with similarities to the transition to turbulence. Correspondingly, no SBs arise when the material is so thin that the entire material can be approximately described as being "interfacial" in nature. We also quantify relaxation in the glass and the nature and origin of highly non-Gaussian particle displacements in the dynamically heterogeneous SB regions at times longer than the caging time.

Confinement effects on the mechanical heterogeneity of polymer fiber glasses

Both polymer fiber glasses and bulk polymer glasses exhibit nonlinear mechanical responses under uniaxial deformation. In polymer fibers, however, polymer chains are confined strongly and the surface area is relatively large compared to their volume. The confinement and the surface may lead to the spatially heterogeneous relaxation of chains in polymer fibers. In this work we perform molecular dynamics simulations and investigate the relation between the heterogeneous dynamics and the nonlinear mechanical responses at a molecular level. Our molecular simulations capture successfully not only the nonlinear mechanical response but also the dependence of mechanical properties on the strain rate of typical polymer glasses as in experiments. We find that the local elastic modulus and the nonaffine displacement are spatially heterogeneous in the pre-yield regime, which results in a lower elastic modulus for polymer fibers than bulk polymer glasses. In the post-yield regime, those mechanical properties become relatively homogeneous. Monomers with large nonaffine displacement are localized mainly at the interfacial region in the pre-yield regime while highly nonaffine monomers are distributed throughout the fibers in the post-yield regime. We show that the nonaffine displacement during deformation relates closely to the mechanical response of the polymer fibers. We also find that in the strain-hardening regime there is a significant difference in the energetic contribution to the stress between polymer fibers and bulk polymers, for which the modulus of the strain-hardening regime of the polymer fibers is smaller than that of bulk polymers.

Effect of polymer-nanoparticle interaction on strain localization in polymer nanopillars

Polymer nanocomposites (PNCs), a class of composites consisting of typically inorganic nanoparticles (NPs) embedded in a polymer matrix, have become an emerging class of materials due to their significant potential to combine the functionality of NPs with the toughness of polymers. However, many applications are limited by their mechanical properties, and a fundamental understanding of NPs' effect on the nonlinear mechanical properties is lacking. In this study, we used molecular dynamics simulations to investigate the influence of NPs on the tendency of a polymer nanopillar to form a shear band. Even though we restrict ourselves to sufficiently low NP loadings that the elastic and yield behaviors are unaffected compared to the pure polymer, the polymer-NP interactions have a surprisingly strong effect on the location of a shear band in the sample. Different polymer-NP interactions have been used to explore their effect on the local structure of materials which is described using a recently developed machine-learned quantity, softness. Our calculations reveal a strong correlation between the strain localization pattern and the local structural signatures. Lastly, we show that weak interactions between NP and polymer matrix can form a soft region near the NP, and this leads to an attraction of the shear band to the NP surface.

Cavitation in soft matter

Cavitation is the sudden, unstable expansion of a void or bubble within a liquid or solid subjected to a negative hydrostatic stress. Cavitation rheology is a field emerging from the development of a suite of materials characterization, damage quantification, and therapeutic techniques that exploit the physical principles of cavitation. Cavitation rheology is inherently complex and broad in scope with wide-ranging applications in the biology, chemistry, materials, and mechanics communities. This perspective aims to drive collaboration among these communities and guide discussion by defining a common core of high-priority goals while highlighting emerging opportunities in the field of cavitation rheology. A brief overview of the mechanics and dynamics of cavitation in soft matter is presented. This overview is followed by a discussion of the overarching goals of cavitation rheology and an overview of common experimental techniques. The larger unmet needs and challenges of cavitation in soft matter are then presented alongside specific opportunities for researchers from different disciplines to contribute to the field.

Distinguishing failure modes in oligomeric polymer nanopillars

Brittle failure is ubiquitous in amorphous materials that are sufficiently cooled below their glass transition temperature, T_g. This catastrophic failure mode is limiting for amorphous materials in many applications, and many fundamental questions surrounding it remain poorly understood. Two challenges that prevent a more fundamental understanding of the transition between a ductile response at temperatures near T_g to brittle failure at lower temperatures are (i) a lack of computationally inexpensive molecular models that capture the failure modes observed in experiments and (ii) the lack of quantitative metrics that can distinguish various failure mechanisms. In this work, we use molecular dynamics simulations to capture ductile-to-brittle transition in glass-forming oligomeric polymer systems where we systematically vary both the temperature relative to T_g and the form of the interaction potential to induce a variety of failure modes. We characterized the effects of this new potential on macroscopic mechanical properties as well as microscopic structural and dynamical response during deformation. Finally, we develop several quantitative metrics to distinguish between different failure modes, and we find that the transition between catastrophic brittle failure, necking, and homogeneous plastic flow is gradual as the temperature is increased.

Identifying structural signatures of shear banding in model polymer nanopillars

Amorphous solids are critical in the design and production of nanoscale devices, but under strong confinement these materials exhibit changes in their mechanical properties which are not well understood. Phenomenological models explain these properties by postulating an underlying defect structure in these materials but do not detail the microscopic properties of these defects. Using machine learning methods, we identify mesoscale defects that lead to shear banding in model polymer nanopillars well below the glass transition temperature as a function of pillar diameter. Our results show that the primary structural features responsible for shear banding on this scale are fluctuations in the diameter of the pillar. Surprisingly, these fluctuations are quite small compared to the diameter of the pillar, less than half of a particle diameter in size. At intermediate pillar diameters, we find that these fluctuations tend to concentrate along the minor axis of shear band planes. We also see the importance of mean "softness" as a classifier of shear banding grow as a function of pillar diameter. Softness is a new field that characterizes local structure and is highly correlated with particle-level dynamics such that softer particles are more likely to rearrange. This demonstrates that softness, a quantity that relates particle-level structure to dynamics on short time and length scales, can predict large time and length scale phenomena related to material failure.

Structure-property relationships from universal signatures of plasticity in disordered solids

When deformed beyond their elastic limits, crystalline solids flow plastically via particle rearrangements localized around structural defects. Disordered solids also flow, but without obvious structural defects. We link structure to plasticity in disordered solids via a microscopic structural quantity, "softness," designed by machine learning to be maximally predictive of rearrangements. Experimental results and computations enabled us to measure the spatial correlations and strain response of softness, as well as two measures of plasticity: the size of rearrangements and the yield strain. All four quantities maintained remarkable commonality in their values for disordered packings of objects ranging from atoms to grains, spanning seven orders of magnitude in diameter and 13 orders of magnitude in elastic modulus. These commonalities link the spatial correlations and strain response of softness to rearrangement size and yield strain, respectively.

The glass transition and interfacial dynamics of single strand fibers of polymers

We investigate the glass transition and interfacial dynamics of single strand fibers of flexible polymers by employing molecular dynamics (MD) simulations along with a coarse grained model. While the polymer fiber has drawn significant attention due to its applicability in tissue engineering and stretchable electronics, its dynamic properties, especially the glass transition temperature (T_g), are yet to be understood at the molecular level. For example, there has been a controversy on the effect of the polymer fiber radius (R) on T_g: T_g decreased with a decrease in R for some polymer fibers, whereas T_g of other polymer fibers was not sensitive to R. In this article, we estimate the bond relaxation time of polymers and evaluate both T_g and fragility (m) as a function of R. We illustrate that T_g of the polymer fiber decreased with a decrease in R monotonically and also that the values of T_g follow faithfully the empirical equation proposed by Keddie et al. as a function of R, which was successfully employed to fit the values of T_g of both polyvinyl alcohol (PVA) fibers and polyethylene (PE) fibers. We also find that the dynamics of polymers at the interface between a polymer fiber and air is faster than that of polymers at the center. By employing Adam-Gibbs theory, we show that the fast interface dynamics of polymer fibers should influence the cooperative motion of monomers, which should be responsible for the decrease in T_g for smaller values of R. Near the interface there are more mobile monomers that participate in the cooperative motions of polymers. Interesting is that due to the curved surface (unlike flat polymer films) the cooperative motion of monomers is anisotropic in polymer fibers.

The dynamics of unentangled polymers during capillary rise infiltration into a nanoparticle packing

Although highly packed polymer nanocomposites (PNCs) are important for a wide array of applications, preparing them remains difficult because of the poor dispersion of NPs at high loading fractions. One method to successfully prepare PNCs with high loadings is through capillary rise infiltration, as previously shown by Huang et al., although the mechanism of polymer infiltration remains largely unknown. We use molecular dynamics simulations to directly simulate the process of capillary rise infiltration, and we show that the polymers follow Lucas-Washburn dynamics. We observe a wetting front that precedes bulk infiltration, and chains belonging to this front are highly adsorbed to NPs. We also investigate the viscosity of the model polymers both globally and locally in supported and free-standing films, and we find reduced viscosity near the surface of the films and increased viscosity near the supporting substrate, similar to the results of local relaxation times. The reduction in the viscosity at the free surface for short, oligomeric polymers is smaller than for higher molecular weight polymers, and the ratio of the surface viscosities is most consistent with the predictions of the Lucas-Washburn equation. Our results introduce the mechanism by which polymers infiltrate a highly packed NP film, which may shed light on better ways to prepare these materials for energy storage applications and protective coatings.

Strain fluctuations and elastic moduli in disordered solids

Recently there has been a surge in interest in using video-microscopy techniques to infer the local mechanical properties of disordered solids. One common approach is to minimize the difference between particle vibrational displacements in a local coarse-graining volume and the displacements that would result from a best-fit affine deformation. Effective moduli are then inferred under the assumption that the components of this best-fit affine deformation tensor have a Boltzmann distribution. In this paper, we combine theoretical arguments with experimental and simulation data to demonstrate that the above does not reveal information about the true elastic moduli of jammed packings and colloidal glasses.

Structural skeleton of preferentially interpenetrated clusters and correlation with shear localization in Mg-Cu-Ni ternary metallic glasses

Inherent hierarchical structure and its effect on shear localization were clarified for ternary Mg-Cu-Ni metallic glasses via molecular dynamics studies based on a newly constructed n-body potential for the system. Assisted by a proposed index to detect the medium-range correlation heterogeneity, it was found that the Cu/Ni-centered icosahedra and specific Mg-centered clusters exhibit a strong preference to interconnect, leading to the formation, over an extended scale, of a percolated network that serves as structural skeleton in the glassy matrix. In constituting the skeleton network, the clusters mainly integrate in an interpenetrating mode, while the noninterpenetrating linkages provide additional reinforcements, jointly consolidating the structural and energetic stability of the skeleton. Furthermore, by monitoring the structural evolution upon compressive deformation, it was revealed that the gradual collapse of the skeleton network is intimately correlated to the mechanical response of metallic glasses and acts as a structural signature of the initiation and propagation of shear bands.

Disordered surface vibrations in jammed sphere packings

We study the vibrational properties near a free surface of disordered spring networks derived from jammed sphere packings. In bulk systems, without surfaces, it is well understood that such systems have a plateau in the density of vibrational modes extending down to a frequency scale ω*. This frequency is controlled by ΔZ = 〈Z〉 - 2d, the difference between the average coordination of the spheres and twice the spatial dimension, d, of the system, which vanishes at the jamming transition. In the presence of a free surface we find that there is a density of disordered vibrational modes associated with the surface that extends far below ω*. The total number of these low-frequency surface modes is controlled by ΔZ, and the profile of their decay into the bulk has two characteristic length scales, which diverge as ΔZ(-1/2) and ΔZ(-1) as the jamming transition is approached.

Atomistic modeling to optimize composition and characterize structure of Ni–Zr–Mo metallic glasses

An interatomic potential was constructed for the Ni–Zr–Mo ternary metal system with the newly proposed long-range empirical formulism, which has been verified to be applicable for fcc, hcp and bcc transition metals and their alloys.