Elucidating the atomistic mechanisms underpinning plasticity in Li-Si nanostructures

Amorphous lithium-silicon (a-Li-Si), especially in nanostructure form, is an attractive high-capacity anode material for next-generation Li-ion batteries. During cycles of charging and discharging, a-Li-Si undergoes substantive inelastic deformation and exhibits microcracking. The mechanical response to repeated lithiation-delithiation eventually results in the loss of electrical contact and consequent decrease of capacity, thus underscoring the importance of studying the plasticity of a-Li-Si nanostructures. In recent years, a variety of phenomenological continuum theories have been introduced that purport to model plasticity and the electro-chemo-mechanical behavior of a-Li-Si. Unfortunately, the micromechanisms and atomistic considerations underlying plasticity in Li-Si material are not yet fully understood and this impedes the development of physics-based constitutive models. Conventional molecular dynamics, although extensively used to study this material, is grossly inadequate to resolve this matter. As is well known, conventional molecular dynamics simulations can only address phenomena with characteristic time scales of (at most) a microsecond. Accordingly, in such simulations, the mechanical behavior is deduced under conditions of very high strain rates (usually, ${10}^{8}\phantom{\rule{4.pt}{0ex}}{\text{s}}^{\ensuremath{-}1}$ or even higher). This limitation severely impacts a realistic assessment of rate-dependent effects. In this work, we attempt to circumvent the time-scale bottleneck of conventional molecular dynamics and provide novel insights into the mechanisms underpinning plastic deformation of Li-Si nanostructures. We utilize an approach that allows imposition of slow strain rates and involves the employment of a new and recently developed potential energy surface sampling method---the so-called autonomous basin climbing---to identify the local minima in the potential energy surface. Combined with other techniques, such as nudged elastic band, kinetic Monte Carlo and transition state theory, we assess the behavior of a-Li-Si nanostructures under tensile strain rates ranging from ${10}^{3}$ to ${10}^{8}\phantom{\rule{4.pt}{0ex}}{\text{s}}^{\ensuremath{-}1}$. We find significant differences in the deformation behavior across the strain rates and discover that the well-known shear transformation zones (widely discussed in the context of amorphous materials) are formed by a ``diffusionlike'' process. We identify the rotation of the shear transformation zone as a key dissipation mechanism.