We report a molecular dynamics simulation of selenium, described by a three-body interaction. The temperatures T g and T c and the structural properties are in agreement with experiment. The mean nearest neighbor coordination number is 2.1. A small prepeak at about 1 Å Ϫ1 can be explained in terms of void correlations. In the intermediate self-scattering function, i.e., the density fluctuation correlation, classical behavior, ␣ and  regimes, is found. We also observe the plateau in the  regime below T g . In a second step, we investigated the heterogeneous and/or homogeneous behavior of the relaxations. At both short and long times the relaxations are homogeneous ͑or weakly heterogeneous͒. In the intermediate time scale, lowering the temperature increases the heterogeneity. We connect these different domains to the vibrational ͑ballistic͒,  and ␣ regimes. We have also shown that the increase in heterogeneity can be understood in terms of relaxations.