We compute the transient dynamics of phonons in contact with high energy "hot" charge carriers in 12 polar and non-polar semiconductors, using a first-principles Boltzmann transport framework. For most materials, we find that the decay in electronic temperature departs significantly from a single-exponential model at times ranging from 1 ps to 15 ps after electronic excitation, a phenomenon concomitant with the appearance of non-thermal vibrational modes. We demonstrate that these effects result from the slow thermalization within the phonon subsystem, caused by the large heterogeneity in the timescales of electron-phonon and phonon-phonon interactions in these materials. We propose a generalized 2-temperature model accounting for the phonon thermalization as a limiting step of electron-phonon thermalization, which captures the full thermal relaxation of hot electrons and holes in semiconductors. A direct consequence of our findings is that, for semiconductors, information about the spectral distribution of electron-phonon and phonon-phonon coupling can be extracted from the multi-exponential behavior of the electronic temperature.Following the seminal works of Kaganov et al. [1] and Allen [2], the thermalization of a system of highly energetic charge carriers with a lattice is frequently understood as an electron-phonon mediated, temperature equilibration process with a single characteristic timescale τ el-ph . Such description, referred to as the two temperature (2T) model, relies on the central assumption that both electrons and phonons remain in distinct thermal equilibria and can therefore be described by timedependent temperatures T el (t) and T ph (t) during the thermal equilibration process. In metals, due to the relative homogeneity of the electron-phonon interactions and the rates of thermalization within the electronic and phononic subsystems, the hypothesis of subsystem-wide thermal equilibrium is generally accurate, and the 2T model has been successful in modeling ultra-fast laser heating [3][4][5], despite some notable deviations from the 2T predictions in graphene and aluminum [6][7][8]. In semiconductors, the highly heterogeneous electron-phonon interactions (e.g. in polar semiconductors with Fröhlich interactions [9]) and, in some cases, the higher lattice thermal conductivity in comparison to metals weaken the hypothesis of a thermalized phononic subsystem [10,11], hence calling for the reexamination of the 2T physical picture in semiconductors.In this context, the advent of first-principles techniques able to predict the mode-and energy-resolved electronphonon [12][13][14] and phonon-phonon interactions [15,16] provides an important opportunity: In their modern implementations [13,16,17], these methods have been able to predict lattice thermal conductivities [18][19][20][21], the temperature-and pressure-dependence of the electronic bandgap [22][23][24][25][26][27][28], electrical conductivities [29,30], and hot carrier dynamics [31,32]. However, to the best of our knowledge and despite these ...