The 30 Doradus (30\,Dor) region in the Large Magellanic Cloud (also known as the Tarantula Nebula) is the nearest massive starburst region, containing the richest sample of massive stars in the Local Group. It is the best possible laboratory to investigate aspects of the formation and evolution of massive stars. Here, we focus on rotation which is a key parameter in the evolution of these objects. We establish the projected rotational velocity, $v_{e}\sin i$, distribution of an unprecedented sample of 216 radial velocity constant ($\rm{\Delta RV\, \leq\, 20 \,km s^{-1}}$) O-type stars in 30\,Dor observed in the framework of the VLT-FLAMES Tarantula Survey (VFTS). The distribution of $v_{e}\sin i$ shows a two-component structure: a peak around 80 $\rm{km s^{-1}}$ and a high-velocity tail extending up to $\sim$600 $\rm{km s^{-1}}$. Around 75% of the sample has 0 $\leq\, v_{e}\sin i \leq$ 200 $\rm{km s^{-1}}$ with the other 25% distributed in the high-velocity tail. The presence of the low-velocity peak is consistent with that found in other studies of late-O and early-B stars. The high-velocity tail is compatible with expectations from binary interaction synthesis models and may be predominantly populated by post-binary interaction, spun-up, objects and mergers. This may have important implications for the nature of progenitors of long-duration gamma ray bursts.

Aims. Using ground based multi-object optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to establish the (projected) rotational velocity distribution for a sample of 216 presumably single O-type stars in 30 Doradus (30 Dor). Methods. We measured projected rotational velocities, \vrot, by means of a Fourier transform method and a profile fitting method applied on a set of isolated spectral lines. We also used an iterative deconvolution procedure to infer the probability density, $\rm{P(\veq)}$, of the equatorial rotational velocity, \veq. Results. The distribution of \vrot\ shows a two-component structure: a peak around 80 \kms\ and a high-velocity tail extending up to $\sim$600 \kms. This structure is also present in the inferred distribution $\rm{P(\veq)}$ with around 80% of the sample having 0 $<$ \veq\, $\leq\, 300$ \kms\ and the other 20% distributed in the high-velocity region. Conclusions. Most of the stars in our sample rotate with a rate less than 20%\ of their break-up velocity. For the bulk of the sample, mass-loss in a stellar wind and/or envelope expansion is not efficient enough to significantly spin down these stars within the first few Myr of evolution. The presence of a sizeable population of fast rotators is compatible with recent population synthesis computations that investigate the influence of binary evolution on the rotation rate of massive stars. Despite the fact that we have excluded stars that show significant radial velocity variations, our sample may have remained contaminated by post-interaction binary products. The fact that the high-velocity tail may be preferentially (and perhaps even exclusively), populated by post-binary interaction products, has important implications for the evolutionary origin of systems that produce gamma-ray bursts.