By finding local minima of a thermodynamic-like potential, we generate jammed packings of frictionless spheres under constant shear stress σ and obtain the yield stress σy by sampling the potential energy landscape. For three-dimensional systems with harmonic repulsion, σy satisfies the finite size scaling with the limiting scaling relation σy∼ϕ-ϕc,∞, where ϕc,∞ is the critical volume fraction ...