<p>This paper extends the previous work [1] for the three-temperature gray radiative transfer equations to the frequency-dependent case. Since the additional frequency variable is considered, the equations are more complicated than those in the gray case. Moreover, opacity may be typically a decreasing function of the frequency variable in applications. At the same spatial location, the equations can be in the optically thick case for low frequency photons, while in the optically thin case for high frequency ones. Thus, the resulting discrete equations can significantly increase the computational cost for opacity having the multi-scale property in multiple frequency radiation. Due to the presence of the radiation-electron coupling, electron-ion coupling, and electron-ion diffusion terms, the model under consideration exhibits strong nonlinearity and strong coupling properties. In this paper, the multigroup method is used to discretize the frequency variable and the <i>H</i><Stack> _<i>N</i> ^<i>T</i> </Stack> method to discretize the angular variable first. Then, within the framework of a unified gas kinetic scheme (UGKS), a multig roup <i>H</i><Stack> _<i>N</i> ^<i>T</i> </Stack>-UGKS method is constructed to solve this complex model iteratively. Furthermore, it can be shown that as the Knudsen number tends to zero, with variations in the electron-ion coupling, absorption, and scattering coefficients, the multigroup <i>H</i><Stack> _<i>N</i> ^<i>T</i> </Stack>-UGKS scheme can converge to numerical schemes for the single-temperature, two-temperature, and the frequency-dependent three-temperature, two-temperature diffusion limit equations, respectively. Finally, several numerical examples are provided to validate the effectiveness and stability of the proposed scheme.</p>