1. Introduction
Massive stars play a critical role in galactic evolution, influencing their surroundings through ionizing radiation, stellar winds, and supernova explosions. Among their fundamental parameters, rotation strongly affects stellar evolution, internal mixing, and final fate.
Observations reveal a bimodal distribution in projected rotational velocity (
$v\sin{i}$
) among O-type stars: a dominant slow-rotator population (
$40 \!-\! 100 \,\mathrm{kmp\,s}^{-1}$
) and a secondary group of fast rotators above
$200\,\mathrm{kmp\,s}^{-1}$
(Ram rez-Agudelo et al. 2013; Holgado et al. Reference Holgado, Simón-Daz, Herrero and Barbá2022). The origin of this high-velocity tail is usually associated with spin-up interaction processes in binaries, such as tidals and mass transfer (Sana et al. Reference Sana, de Mink and de Koter2012; de Mink et al. Reference de Mink, Langer, Izzard, Sana and de Koter2013).
Cygnus OB2, located within the Cygnus X complex, hosts hundreds of OB-type stars at various evolutionary stages, providing a natural laboratory for probing massive star formation and evolution.
2. Rotational velocity distribution for Cygnus OB2
Through a comprehensive spectroscopic campaign, we analyzed the projected rotational velocities of a sample of 70 O-type stars in the core of Cygnus OB2. Our dataset includes both new observations (Galán-Diéguez et al. Reference Galán-Diéguez, Berlanas and Herrero2025) and stars previously studied by Berlanas et al. (Reference Berlanas, Herrero and Comerón2020). We derive projected rotational velocities (
$v\sin{i}$
) using the iacob-broad tool, which combines the Fourier Transform and Goodness-of-Fit methods to disentangle rotational and macroturbulent broadening. The results show:
-
• Only
$\!\sim\! 9$
O-type stars exceed
$200\,\mathrm{kmp\,s}^{-1}$
. -
• None exceed
$300\,\mathrm{kmp\,s}^{-1}$
. -
• The expected
$\sim\!10\%$
of rapid rotators predicted by binary interaction models is missing.
Figure 1 compares Cygnus OB2 with other regions (Galán-Diéguez et al. Reference Galán-Diéguez, Berlanas and Herrero2025). While the Milky Way population (Holgado et al. Reference Holgado, Simón-Daz, Herrero and Barbá2022) and 30 Doradus (Ram rez-Agudelo et al. 2013) exhibit a pronounced fast-rotator tail, Cygnus OB2 emerges as an outlier. Surprisingly, the O-type population in Cygnus OB2 does not follow this trend: the expected fast-rotating tail is missing.
Bootstrapped distributions of projected rotational velocities for O-type stars from various regions within the local Universe (adapted from Galán-Diéguez et al. Reference Galán-Diéguez, Berlanas and Herrero2025). The histograms in black and green illustrate the
$v\sin{i}$
measurements for Galactic stars in the solar neighborhood (Holgado et al. Reference Holgado, Simón-Daz, Herrero and Barbá2022), and 30 Doradus (Ram rez-Agudelo et al. 2013), respectively. The orange histogram depicts the rotational velocity distribution for Cygnus OB2, including O-type stars from Berlanas et al. (Reference Berlanas, Herrero and Comerón2020) and those lately reclassified from the B0 sample (Galán-Diéguez et al. Reference Galán-Diéguez, Berlanas and Herrero2025). The error bars indicate the standard deviation in each velocity bin derived from the bootstrap analysis.

This result is especially unexpected given that massive stars are thought to form predominantly in multiple systems, with a large fraction expected to experience binary interactions (Sana et al. Reference Sana, de Mink and de Koter2012). Such processes are predicted to generate a substantial population of fast rotators (de Mink et al. Reference de Mink, Langer, Izzard, Sana and de Koter2013), yet their absence in Cygnus OB2 poses a challenge to current models of massive star formation and evolution.
3. Scenarios for the lack of fast rotators in Cygnus OB2
The lack of fast rotators in Cygnus OB2, contrary to the predictions of binary evolution theory and evidence from other massive clusters, remains a puzzling result. Several non-exclusive explanations may help to account for this anomaly: Young evolutionary stage: Binary interaction, in particular mass transfer, is the dominant channel for producing rapid rotators. Since Cygnus OB2 is very young (1-6 Myr, Berlanas et al. Reference Berlanas, Herrero and Comerón2020), many systems may not yet have undergone such episodes. However, other clusters of comparable or younger ages – such as 30 Doradus – exhibit populations of fast rotators. This suggests that youth alone cannot fully explain the deficit of a high-velocity tail in Cygnus OB2.
Runaway scenario: Stars spun up through binary interaction may have been expelled from the association as runaways, either by supernova explosions or by dynamical encounters. For instance, a star moving at
$30\,\mathrm{kmp\,s}^{-1}$
would leave the central 1
$^\circ$
of Cygnus OB2 in less than 1 Myr. Some of the missing fast rotators may therefore lie outside the surveyed region.
Molecular cloud properties: The star-forming environment in Cygnus OB2 may have favored slower initial spins or reduced the binary fraction. These would naturally leading to a lack of a fast-rotating tail in the present-day population.
4. Conclusions
Our spectroscopic analysis of 70 O-type stars in Cygnus OB2 reveals a striking absence of rapid rotators, with no stars exceeding
$300\,\mathrm{kmp\,s}^{-1}$
and only a small fraction surpassing
$200\,\mathrm{kmp\,s}^{-1}$
. This contrasts sharply with other massive star-forming regions, such as 30 Doradus or the Galactic field, where a pronounced high-velocity tail is observed.
The deficit of fast rotators in Cygnus OB2 challenges current expectations from binary evolution models, which predict that
$\sim$
10% of O-type stars should appear as rapid rotators as a consequence of mass transfer or mergers. Several scenarios may contribute to this discrepancy, including the youth of the association, the possible ejection of spun-up systems as runaways, and environmental effects related to the natal molecular cloud.
Further progress will require extending the spatial coverage of the survey to search for high-velocity runaways, and performing detailed binary population synthesis tailored to the Cygnus OB2 environment. Together, these efforts will help to clarify whether the absence of fast rotators reflects a transient evolutionary stage, an observational bias, or a genuine environmental imprint on the rotational properties of massive stars.
Acknowledgements
The authors acknowledge support funding from the Ministry of Science, Innovation and Universities (MCIU) through the Spanish State Research Agency (AEI) via grant PID2021-122397NB-C21 (cofunded by the European Regional Development Fund, FEDER) and the Severo Ochoa Program 2020-2023 (CEX2019-000920-S). SRB also thanks financial support from NextGeneration EU/PRTR and MIU (UNI/551/2021) through grant Margarita Salas-ULL. MAM project received support from the “La Caixa” Foundation (ID 100010434) under the fellowship code LCF/BQ/PI23/11970035. This work is partially based on observations made with the three following telescopes: (i) we thank the Nordic Optical Telescope, owned in collaboration by the University of Turku and Aarhus University, and operated jointly by Aarhus University, the University of Turku and the University of Oslo, representing Denmark, Finland and Norway, the University of Iceland and Stockholm University at the Observatorio del Roque de los Muchachos, La Palma, Spain, of the Instituto de Astrofisica de Canarias. The data presented in this work were obtained in part with ALFOSC, which is provided by the Instituto de Astrofisica de Andalucia (IAA) under a joint agreement with the University of Copenhagen and NOT; (ii) the Gran Telescopio Canarias is installed at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias, on the island of La Palma. GTC data was obtained with the instrument OSIRIS, built by a Consortium led by the Instituto de Astrofísica de Canarias in collaboration with the Instituto de Astronomía of the Universidad Autónoma de México. OSIRIS was funded by GRANTECAN and the National Plan of Astronomy and Astrophysics of the Spanish Government; (iii) the Isaac Newton Telescope and its service mode are operated on the island of La Palma by the Isaac Newton Group of Telescopes in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. During the spectroscopic data reduction we used PyRAF, a Python environment for IRAF. IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Association of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the National Science Foundation. In addition to the codes explicitly cited in the text, this paper makes use of the following packages: matplotlib (Hunter Reference Hunter2007), numpy (Harris et al. Reference Harris, Millman and van der Walt2020), and scipy (Virtanen et al. Reference Virtanen, Gommers and Oliphant2020).




vsini