  While perfect crystals may exhibit a purely elastic response to shear all the
way to yielding, the response of amorphous solids is punctuated by plastic
events. The prevalence of this plasticity depends on the number of particles
$N$ of the system, with the average strain interval before the first plastic
event, $\overline{\Delta\gamma}$, scaling like $N^\alpha$ with $\alpha$
negative: larger samples are more susceptible to plasticity due to more
numerous disorder-induced soft spots. In this paper we examine this scaling
relation in ultra-stable glasses prepared with the Swap Monte Carlo algorithm,
with regard to the possibility of protocol-dependent scaling exponent, which
would also imply a protocol dependence in the distribution of local yield
stresses in the glass. We show that, while a superficial analysis seems to
corroborate this hypothesis, this is only a pre-asymptotic effect and in fact
our data can be well explained by a simple model wherein such protocol
dependence is absent.
