For a few decades, the influence of a magnetic field on the aggregation process of superparamagnetic colloids has been well known on short time scale. However, the accurate study of the equilibrium state is still challenging on some aspects. On the numerical aspect, current simulations have only access to a restricted set of experimental conditions due to the computational cost of long range interactions in many body systems. In the present paper, we numerically explore a new range of parameters thanks to sped up numerical simulations validated by a recent experimental and numerical study. We first show that our simulations reproduce results from previous study in well-established conditions. Then we show that unexpectedly long chains are observed for higher volume fractions and intermediate fields. We also present theoretical developments taking into account the interaction between the chains which are able to reproduce the data that we obtained with our simulations. We finally confirm this model thanks to experimental data.