  We study the finite-shear-rate rheology of disordered solids by means of
molecular dynamics simulations in two dimensions. By systematically varying the
damping magnitude $\zeta$ in the low-temperature limit, we identify two well
defined flow regimes, separated by a thin (temperature-dependent) crossover
region. In the overdamped regime, the athermal rheology is governed by the
competition between elastic forces and viscous forces, whose ratio gives the
Weissenberg number $Wi= \zeta \dot\gamma$ (up to elastic parameters); the
macroscopic stress $\Sigma$ follows the frequently encountered Herschel-Bulkley
law $\Sigma= \Sigma\_0 + k \sqrt{Wi}$, with yield stress
$\Sigma\_0\textgreater{}0$. In the underdamped (inertial) regime, dramatic
changes in the rheology are observed for low damping: the flow curve becomes
non-monotonic. This change is not caused by longer-lived correlations in the
particle dynamics at lower damping; instead, for weak dissipation, the sample
heats up considerably due to, and in proportion to, the driving. By suitably
thermostatting more or less underdamped systems, we show that their rheology
only depends on their kinetic temperature and the shear rate, rescaled with
Einstein's vibration frequency.
