Hard spheres interacting through a square well potential were simulated by using two different methods: Brownian cluster dynamics (BCD) and event driven Brownian dynamics (EDBD). The structure of the equilibrium states obtained by both methods was compared and found to be almost identical. Self-diffusion coefficients (D) were determined as a function of the interaction strength. The same values were found by using BCD or EDBD. Contrary to EDBD, BCD allows one to study the effect of bond rigid...