A numerical study of the onset of convection in rotating internally heated self-gravitating fluid spheres is presented. The exploration of the stability of the conduction state versus the Taylor and Prandtl numbers supplies a detailed idea of the laws that fulfil the four types of solutions obtained at low Prandtl numbers. The main result found is that axisymmetric (torsional) modes of convection are preferred at high Taylor numbers in the zero-Prandtl-number limit. This instability appears at low Rayleigh numbers and gives rise to an oscillating single vortex of very high frequency.