We present a detailed numerical investigation of the tunability of a diffusive random laser when Mie resonances are excited. We solve a multimode diffusion model and calculate multiple light scattering in presence of optical gain which includes dispersion in both scattering and gain, without any assumptions about the β parameter. This allows us to investigate a realistic photonic glass made of latex spheres and rhodamine and to quantify both the lasing wavelength tunability range and the lasing threshold. Beyond what is expected by diffusive monochromatic models, the highest threshold is found when the competition between the lasing modes is strongest and not when the lasing wavelength is furthest from the maximum of the gain curve.Random lasers (RL) are mirror-less lasing systems which have attracted a lot of interest due to their structural simplicity. Nowadays they have been studied in a vast variety of scattering systems ranging from semiconductor powder to biological tissue and biocompatible materials [1]. Random lasing originates from a complex out-of-equilibrium phenomenon with rich multimodes features [2] and surprising statistical features [3]. Despite its potential for practical applications [1], random lasing technology is still in its infancy with pioneering applications such as low coherence light source [4] and biosensing [5]. One of the factors that has limited practical applications is the difficulty of controlling the frequency and directionality of the emission. In conventional lasers the lasing emission can be tuned by engineering the high finesse cavity which provides the feedback and thus defines the lasing mode. Instead, feedback in RL is provided by multiple scattering and the lasing emission properties are determined by the complex interplay between gain and losses. Recent experiments have shown lasing emission controlled by exploiting scattering dispersion via resonant scattering sustained by spherical particles [6,7] or by gain dispersion achieved by artificially increasing absorption in a spectral band [8]. Active tuning of the lasing properties has also been achieved by shaping the pump profile to selectively excite one or a few lasing modes [9][10][11].Different theoretical approaches to model random lasing action have been developed which combine multiple scattering and gain. For uncorrelated random systems in which interference between the scattered waves can be neglected, diffusive models are very accurate even in presence of optical gain [12,13] and they provide the time evolution of the lasing process and a smooth lasing spectrum with no spiking lasing behaviour [14]. The radiative transport model with gain can also be solved for * Corresponding author: michele.gaio@kcl.ac.uk instance with Monte Carlo simulations which consider a random walk of photons [15,16] and in which amplification of single paths can be important in defining the spectral properties [17], and by solving the complete radiative transfer equations [18]. These approaches allow the study of large systems (> 100s...