Penghui Cao

Strain-rate and temperature-driven transition in the shear transformation zone for two-dimensional amorphous solids

We couple the recently developed self-learning metabasin escape algorithm, which enables efficient exploration of the potential energy surface (PES), with shear deformation to elucidate strain-rate and temperature effects on the shear transformation zone (STZ) characteristics in two-dimensional amorphous solids. In doing so, we report a transition in the STZ characteristics that can be obtained through either increasing the temperature or decreasing the strain rate. The transition separates regions having two distinct STZ characteristics. Specifically, at high temperatures and high strain rates, we show that the STZs have characteristics identical to those that emerge from purely strain-driven, athermal quasistatic atomistic calculations. At lower temperatures and experimentally relevant strain rates, we use the newly coupled PES + shear deformation method to show that the STZs have characteristics identical to those that emerge from a purely thermally activated state. The specific changes in STZ characteristics that occur in moving from the strain-driven to thermally activated STZ regime include a 33% increase in STZ size, faster spatial decay of the displacement field, a change in the deformation mechanism inside the STZ from shear to tension, a reduction in the stress needed to nucleate the first STZ, and finally a notable loss in characteristic quadrupolar symmetry of the surrounding elastic matrix that has previously been seen in athermal, quasistatic shear studies of STZs.

Density functional theory calculation of edge stresses in monolayer MoS2

We utilize density functional theory to calculate the edge energy and edge stress for monolayer MoS2 nanoribbons. In contrast to previous reports for graphene, for both armchair and zigzag chiralities, the edge stresses for MoS2 nanoribbons are found to be tensile, indicating that their lowest energy configuration is one of compression in which Mo-S bond lengths are shorter than those in a bulk, periodic MoS2 monolayer. The edge energy and edge stress is found to converge for both chiralities for nanoribbon widths larger than about 1 nm.

Surface shear-transformation zones in amorphous solids.

We perform a systematic study of the characteristics of shear transformation zones (STZs) that nucleate at free surfaces of two-dimensional amorphous solids subject to tensile loading using two different atomistic simulation methods, the standard athermal, quasistatic (AQ) approach and our recently developed self-learning metabasin escape (SLME) method, to account for the finite temperature and strain-rate effects. In the AQ, or strain-driven limit, the nonaffine displacement fields of surface STZs decay exponentially away from their centers at similar decay rates as their bulk counterparts, though the direction of maximum nonaffine displacement is tilted away from the tensile axis due to surface effects. Using the SLME method at room temperature and at the high strain rates that are seen in classical molecular dynamics simulations, the characteristics for both bulk and surface STZs are found to be identical to those seen in the AQ simulations. However, using the SLME method at room temperature and experimentally relevant strain rates, we find a transition in the surface STZ characteristics where a loss in the characteristic angular tensile-compression symmetry is observed. Finally, the thermally activated surface STZs exhibit a slower decay rate in the nonaffine displacement field than do strain-driven surface STZs, which is characterized by a larger drop in potential energy resulting from STZ nucleation that is enabled by the relative compliance of the surface as compared to the bulk.

Atomistic modeling at experimental strain rates and timescales

Modeling physical phenomena with atomistic fidelity and at laboratory timescales is one of the holy grails of computational materials science. Conventional molecular dynamics (MD) simulations enable the elucidation of an astonishing array of phenomena inherent in the mechanical and chemical behavior of materials. However, conventional MD, with our current computational modalities, is incapable of resolving timescales longer than microseconds (at best). In this short review article, we briefly review a recently proposed approach—the so-called autonomous basin climbing (ABC) method—that in certain instances can provide valuable information on slow timescale processes. We provide a general summary of the principles underlying the ABC approach, with emphasis on recent methodological developments enabling the study of mechanically-driven processes at slow (experimental) strain rates and timescales. Specifically, we show that by combining a strong physical understanding of the underlying phenomena, kinetic Monte Carlo, transition state theory and minimum energy pathway methods, the ABC method has been found to be useful in a variety of mechanically-driven problems ranging from the prediction of creep-behavior in metals, constitutive laws for grain boundary sliding, void nucleation rates, diffusion in amorphous materials to protein unfolding. Aside from reviewing the basic ideas underlying this approach, we emphasize some of the key challenges encountered in our own personal research work and suggest future research avenues for exploration.

Atomistic Simulation of the Rate-Dependent Ductile-to-Brittle Failure Transition in Bicrystalline Metal Nanowires

The mechanical properties and plastic deformation mechanisms of metal nanowires have been studied intensely for many years. One of the important yet unresolved challenges in this field is to bridge the gap in properties and deformation mechanisms reported for slow strain rate experiments (∼10-2 s-1), and high strain rate molecular dynamics (MD) simulations (∼108 s-1) such that a complete understanding of strain rate effects on mechanical deformation and plasticity can be obtained. In this work, we use long time scale atomistic modeling based on potential energy surface exploration to elucidate the atomistic mechanisms governing a strain-rate-dependent incipient plasticity and yielding transition for face centered cubic (FCC) copper and silver nanowires. The transition occurs for both metals with both pristine and rough surfaces for all computationally accessible diameters (<10 nm). We find that the yield transition is induced by a transition in the incipient plastic event from Shockley partials nucleated on primary slip systems at MD strain rates to the nucleation of planar defects on non-Schmid slip planes at experimental strain rates, where multiple twin boundaries and planar stacking faults appear in copper and silver, respectively. Finally, we demonstrate that, at experimental strain rates, a ductile-to-brittle transition in failure mode similar to previous experimental studies on bicrystalline silver nanowires is observed, which is driven by differences in dislocation activity and grain boundary mobility as compared to the high strain rate case.

The Role of Binding Site on the Mechanical Unfolding Mechanism of Ubiquitin

We apply novel atomistic simulations based on potential energy surface exploration to investigate the constant force-induced unfolding of ubiquitin. At the experimentally-studied force clamping level of 100 pN, we find a new unfolding mechanism starting with the detachment between β5 and β3 involving the binding site of ubiquitin, the Ile44 residue. This new unfolding pathway leads to the discovery of new intermediate configurations, which correspond to the end-to-end extensions previously seen experimentally. More importantly, it demonstrates the novel finding that the binding site of ubiquitin can be responsible not only for its biological functions, but also its unfolding dynamics. We also report in contrast to previous single molecule constant force experiments that when the clamping force becomes smaller than about 300 pN, the number of intermediate configurations increases dramatically, where almost all unfolding events at 100 pN involve an intermediate configuration. By directly calculating the life times of the intermediate configurations from the height of the barriers that were crossed on the potential energy surface, we demonstrate that these intermediate states were likely not observed experimentally due to their lifetimes typically being about two orders of magnitude smaller than the experimental temporal resolution.

Understanding the mechanisms of amorphous creep through molecular simulation

Significance The individual and collective molecular displacements in an amorphous solid undergoing plastic deformation are simulated by an atomistic method that allows incremental motions to be observed over a time window of fractions of seconds. Because the timescale matches well with the conditions of experimental measurements, simulation details provide dynamical evidence for the fundamental mechanisms of amorphous creep. In particular, knowledge of the interplay between diffusion (flow) and mechanical deformation processes enables us to explain the stress and temperature behavior of the experimental data as well as the validity of model descriptions of molecular mechanisms in terms of spatially and temporally heterogeneous fluctuations. Molecular processes of creep in metallic glass thin films are simulated at experimental timescales using a metadynamics-based atomistic method. Space–time evolutions of the atomic strains and nonaffine atom displacements are analyzed to reveal details of the atomic-level deformation and flow processes of amorphous creep in response to stress and thermal activations. From the simulation results, resolved spatially on the nanoscale and temporally over time increments of fractions of a second, we derive a mechanistic explanation of the well-known variation of creep rate with stress. We also construct a deformation map delineating the predominant regimes of diffusional creep at low stress and high temperature and deformational creep at high stress. Our findings validate the relevance of two original models of the mechanisms of amorphous plasticity: one focusing on atomic diffusion via free volume and the other focusing on stress-induced shear deformation. These processes are found to be nonlinearly coupled through dynamically heterogeneous fluctuations that characterize the slow dynamics of systems out of equilibrium.

Superplastic Creep of Metal Nanowires from Rate-Dependent Plasticity Transition

Understanding the time-dependent mechanical behavior of nanomaterials such as nanowires is essential to predict their reliability in nanomechanical devices. This understanding is typically obtained using creep tests, which are the most fundamental loading mechanism by which the time-dependent deformation of materials is characterized. However, due to existing challenges facing both experimentalists and theorists, the time-dependent mechanical response of nanowires is not well-understood. Here, we use atomistic simulations that can access experimental time scales to examine the creep of single-crystal face-centered cubic metal (Cu, Ag, Pt) nanowires. We report that both Cu and Ag nanowires show significantly increased ductility and superplasticity under low creep stresses, where the superplasticity is driven by a rate-dependent transition in defect nucleation from twinning to trailing partial dislocations at the micro- or millisecond time scale. The transition in the deformation mechanism also governs a corresponding transition in the stress-dependent creep time at the microsecond (Ag) and millisecond (Cu) time scales. Overall, this work demonstrates the necessity of accessing time scales that far exceed those seen in conventional atomistic modeling for accurate insights into the time-dependent mechanical behavior and properties of nanomaterials.