Abstract
We devise a generic strategy and simple numerical models for multi-component metallic glasses for which the swap Monte Carlo algorithm can produce highly stable equilibrium configurations equivalent to experimental systems cooled more than $10^7$ times slower than in conventional simulations. This paves the way for a deeper understanding of thermodynamic, dynamic, and mechanical properties of metallic glasses. As a first application, we extend configurational entropy measurements down to the experimental glass temperature, and demonstrate a qualitative evolution of the mechanical response of metallic glasses of increasing stability towards brittleness.