Understanding the yielding of glass-forming systems upon shearing is notoriously difficult since it is a strong non-equilibrium effect. Here we show that the concept of the potential energy landscape (PEL), developed for the quiescent state, can be extended to shearing. When introducing an appropriate coarse graining of the extended PEL for sheared systems, one can distinguish two fundamentally different types of plastic events, namely inherent structures (IS) vs. minimized structure (MS) transitions. We apply these concepts to non-cyclic and cyclic shearing of small systems, which allow us to characterize the properties of elementary plastic events. Whereas the general properties of the stress-strain curves are similar to larger systems, a closer analysis reveals significantly different properties. This allows one to identify the impact of the elastic coupling in larger systems. The concept of MS enables us to relate the stress overshoot of a single trajectory to the emergence of an MS transition. Furthermore, the occurrence of limit cycles can be characterized in great detail for the small systems and connections to the properties of the PEL can be formulated. Possible implications of our small-system results for the macroscopic limit are discussed.PACS numbers:
I. INTRODUCTIONA common feature of amorphous glass forming materials like metallic glasses, colloidal particles or foams is an overshoot in the stress-strain curve when applying a constant shear rate when starting at rest. During this overshoot period, the material performs a transition from solid-like to fluid-like behavior. It has been the focus of recent work to study the process of yielding in more detail [1][2][3][4][5][6][7][8]. There, we find many different definitions of yield points, each characterizing a different aspect of yielding.The simplest approach of defining a yield strain in glasses is taking the strain at the overshoot maximum γ peak . For the data of 3D Lennard-Jones-type systems, shown in [1], one obtains a value around 0.09, largely independent of system size. Another approach [5] uses an overlap function to determine whether there has been a larger rearrangement between two configurations. The yield is then defined as the point where half of the systems had such a transition. The same results are obtained by using susceptibilities [9]. When using oscillatory shear also the equality of the storage modulus G and the loss modulus G can be used to characterize the transition from elastic to non-elastic behavior [10], although yielding a slightly larger value. This rheological crossover strain can also be related to a sharp transition in the anisotropy of the structure factor [8].Interestingly, via oscillatory shear experiments and simulations in the low-temperature limit a critical strain γ nl can be defined, marking the onset of non-linear behavior [1,[10][11][12][13]. Whereas for strain amplitudes γ max < γ nl one observes limit cycles of the trajectory [1,6,10,11], for γ max > γ nl the system irreversibly leaves its initial region...