We use molecular dynamics computer simulations to investigate the relaxation dynamics of a simple model for a colloidal gel at a low volume fraction. We find that due to the presence of the open spanning network this dynamics shows at low temperature a non-trivial dependence on the wave-vector which is very different from the one observed in dense glass-forming liquids. At high wave vectors the relaxation is due to the fast cooperative motion of the branches of the gel network, whereas at low wave vectors the overall rearrangements of the heterogeneous structure produce the relaxation process.PACS numbers: 82.70. Gg, 82.70.Dd, 64.70.Pf, 67.40.Fd In gels, the deep connection existing between their unusual dynamics and the open network characterizing their structure is still not understood [1,2]. In the dramatic slowing down of the dynamics accompanying the gel formation, the relaxation functions are often stretched and/or, most remarkably, compressed, i.e. the time correlators decay faster than an exponential [3,4,5,6,7], showing a complex dependence on length scale. These findings suggest that different relaxation mechanisms interplay at a microscopic level, which have not been elucidated yet [8,9]. In the present paper, we show that the formation of the gel network does induce a nontrivial length scale dependence of the dynamics in a simple model for colloidal gels. We use molecular dynamics computer simulations to study the gel formation from the equilibrium sol phase. Our results give evidence that in the incipient gel, the relaxation at high wave vectors is due to the fast cooperative motion of pieces of the gel structure, whereas at low wave vectors the overall rearrangements of the heterogeneous gel make the system relax via a stretched exponential decay of the time correlators. The coexistence of such diverse relaxation mechanisms is determined by the formation of the gel network (i.e. the onset of the elastic response of the system) and it is characterized by a typical crossover length which is of the order of the network mesh size. This is the first work where such a characterization of the gel dynamics in colloidal systems has been achieved, thus making important progress as compared to previous numerical studies [9,10,11,12].In colloidal suspensions at low volume fractions gelation competes and/or interplays with phase separation. As a consequence, coarsening or ordering processes due to the underlying thermodynamics often interfere with the gel dynamics. Whereas in the experiments the time scale typical of the micro or macro-phase separation is often much longer than the observation time scales [3,4], this is not the case in numerical studies using traditional models for colloidal suspensions, where the investigation of the gel dynamics has been until now severely hindered [8,9,10,11]. In order to overcome this problem, we have developed a model in which directional interactions are able to produce a persistent gel network at relatively high temperatures, where phase separation does not occur. The...