Search

This is an Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited., Context. The Carina Nebula is one of the major massive star-forming regions in the Galaxy. Its relatively nearby distance (2.35 kpc) makes it an ideal laboratory for the study of massive star formation, structure, and evolution, both for individual stars and stellar systems. Thanks to the high-quality spectra provided by the Gaia-ESO survey and the LiLiMaRlin library, as well as Gaia EDR3 astrometry, a detailed and homogeneous spectroscopic characterization of its massive stellar content can be carried out. Aims. Our main objective is to spectroscopically characterize all massive members of the Carina Nebula in the Gaia-ESO survey footprint to provide an updated census of massive stars in the region and an updated estimate of the binary fraction of O stars. Methods. We performed accurate spectral classification using an interactive code that compares spectra with spectral libraries of OB standard stars, as well as line-based classic methods. We calculated membership using our own algorithm based on Gaia EDR3 astrometry. To check the correlation between the spectroscopic n-qualifier and the rotational velocity, we used a semi-automated tool for the line-broadening characterization of OB stars based on a combined Fourier transform and goodness-of-fit methodology. Results. The Gaia-ESO survey sample of massive OB stars in the Carina Nebula consists of 234 stars. The addition of brighter sources from the Galactic O-Star Spectroscopic Survey and additional sources from the literature allows us to create the most complete census of massive OB stars so far in the region. It contains a total of 316 stars, with 18 of them in the background and 4 in the foreground. Of the 294 stellar systems in Car OB1, 74 are of O type, 214 are of nonsupergiant B type, and 6 are of WR or nonO supergiant (II to Ia) spectral class. We identify 20 spectroscopic binary systems with an O-star primary, of which 6 are reported for the first time, and another 18 with a B-star primary, of which 13 are new detections. The average observed double-lined binary fraction of O-type stars in the surveyed region is 0.35, which represents a lower limit. We find a good correlation between the spectroscopic n-qualifier and the projected rotational velocity of the stars. The fraction of candidate runaways among the stars with and without the n-qualifier is 4.4% and 2.4%, respectively, although nonresolved double-lined binaries could be contaminating the sample of fast rotators. © The Authors 2023., This paper is based mainly on data products from spectroscopic observations made with ESO Telescopes at the Paranal Observatory under programme ID 188.B-3002. These data products have been processed by the Cambridge Astronomy Survey Unit (CASU) at the Institute of Astronomy, University of Cambridge, and by the FLAMES/UVES reduction team at INAF/Osservatorio Astrofisico di Arcetri. These data have been obtained from the Gaia-ESO Survey Data Archive, prepared and hosted by the Wide Field Astronomy Unit, Institute for Astronomy, University of Edinburgh, which is funded by the UK Science and Technology Facilities Council. This work was partly supported by the European Union FP7 programme through ERC grant number 320360 and by the Leverhulme Trust through grant RPG-2012-541. We acknowledge the support from INAF and Ministero dell’Istruzione, dell’Università e della Ricerca (MIUR) in the form of the grant ‘Premiale VLT 2012’. The results presented here benefit from discussions held during the Gaia-ESO workshops and conferences supported by the ESF (European Science Foundation) through the GREAT Research Network Programme. Additional spectra were obtained using the http://www.lco.cl/irenee-du-pont-telescope/ at the Observatorio de Las Campanas (LCO) and the https://www.eso.org/public/teles-instr/lasilla/mpg22/ at the Observatorio de La Silla (LSO). This work has made use of data from the European Space Agency (ESA) mission https://www.cosmos.esa.int/gaia, processed by the Gaia Data Processing and Analysis Consortium (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. This research is partially funded by the Spanish Government Ministerio de Ciencia e Innovación and Agencia Estatal de Investigación (MCIN/AEI/10.130 39/501 100 011 033/FEDER, UE) through grants PGC2018-093 741-B-C21/C22, PGC2018-095 049-B-C21/C22 and PID2021-122 397NB-C21/C22. SRB also acknowledges funding by MCIN under the Juan de la Cierva – Formación grant (contract FJC 2020-045 785-I) and NextGeneration EU/PRTR and MIU (UNI/551/2021) through grant Margarita Salas-ULL. AH also acknowledges support by the Severo Ochoa Program through CEX2019-000920-S. EJA also acknowledges financial support from the State Agency for Research of the Spanish MCIU through the ‘Center of Excellence Severo Ochoa’ award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709). MB is supported through the Lise Meitner grant from the Max Planck Society. We acknowledge support by the Collaborative Research centre SFB 881 (projects A5, A10), Heidelberg University, of the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation). This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (Grant agreement No. 949173)., With funding from the Spanish government through the "Severo Ochoa Centre of Excellence" accreditation (CEX2021-001131-S).

The empirical distribution of projected rotational velocities (vsini) in massive O-type stars is characterized by a dominant slow velocity component and a tail of fast rotators. Binary interaction has been proposed to play a dominant role in the formation of this tail. We perform a complete and homogeneous search for empirical signatures of binarity in a sample of 54 fast-rotating stars with the aim of evaluating this hypothesis. This working sample has been extracted from a larger sample of 415 Galactic O-type stars which covers the full range of vsini values. We use new and archival multi-epoch spectra in order to detect spectroscopic binary systems. We complement this information with Gaia proper motions and TESS photometric data to aid in the identification of runaway stars and eclipsing binaries, respectively. The identified fraction of single-lined spectroscopic binary (SB1) systems and apparently single stars among the fast-rotating sample is $\sim$18% and $\sim$70%, respectively. When comparing these percentages with those corresponding to the slow-rotating sample we find that our sample of fast rotators is characterized by a slightly larger percentage of SB1 systems ($\sim$18% vs. $\sim$13%) and a considerably smaller fraction of clearly detected SB2 systems (8% vs. 33%). Overall, there seems to be a clear deficit of spectroscopic binaries (SB1+SB2) among fast-rotating O-type stars ($\sim$26% vs. $\sim$46%). On the contrary, the fraction of runaway stars is significantly higher in the fast-rotating domain ($\sim$33-50%) than among those stars with vsini < 200 km/s. Lastly, almost 65% of the apparently single fast-rotating stars are runaways. Our empirical results seem to be in good agreement with the idea that the tail of fast-rotating O-type stars (with vsini > 200 km/s) is mostly populated by post-interaction binary products., 33 pages, 16 figures, accepted for publication in "Astronomy and Astrophysics"

We aim to hunt for massive binaries hosting a black hole companion (OB+BH) and establish the natal mass-ratio distribution of massive stars at the subsolar metallicity environment of the Large Magellanic Cloud (LMC). We use the shift-and-add grid disentangling technique to characterize the hidden companions in 51 SB1 O-type and evolved B-type binaries in the LMC monitored in the framework of the Tarantula Massive Binary Monitoring (TMBM). Out of the 51 SB1 systems, 43 (84%) are found to have non-degenerate stellar companions, of which 28 are confident detections, and 15 are less certain (SB1: or SB2:). Of these 43 targets, one is found to be a triple (VFTS 64), and two are found to be quadruples (VFTS 120, 702). The remaining eight targets (16%) retain an SB1 classification. Aside from the unambiguous case of VFTS 243, analysed in detailed in a separate paper, we identify two additional OB+BH candidates: VFTS 514 and VFTS 779. Additional black holes may be present in the sample but at a lower probability. Our study firmly establishes a virtually flat natal mass-ratio distribution for O-type stars at LMC metallicity, covering the entire mass-ratio range (0.05 < q < 1) and periods in the range 0 < log P < 3 [d]. The nature of the OB+BH candidates should be verified through future monitoring, but the frequency of OB+BH candidates is generally in line with recent predictions at LMC metallicity., 41 pages (14 main article + 27 appendix), accepted to A&A

Multiplicity is a ubiquitous characteristic of massive stars. Multiple systems offer us a unique observational constraint on the formation of high-mass systems. Herschel 36 A is a massive triple system composed of a close binary (Ab1-Ab2) and an outer component (Aa). We measured the orbital motion of the outer component of Herschel 36 A using infrared interferometry with the AMBER and PIONIER instruments of ESO’s Very Large Telescope Interferometer. Our immediate aims are to constrain the masses of all components of this system and to determine if the outer orbit is co-planar with the inner one. Reported spectroscopic data for all two components of this system and our interferometric data allow us to derive full orbital solutions for the outer orbit Aa-Ab and the inner orbit Ab1-Ab2. For the first time, we derive the absolute masses of mAa = 22.3 ± 1.7, mAb1 = 20.5 ± 1.5, and mAb2 = 12.5 ± 0.9 M⊙. Despite not being able to resolve the close binary components, we infer the inclination of their orbit by imposing the same parallax as the outer orbit. Inclinations derived from the inner and outer orbits imply a modest difference of about 22° between the orbital planes. We discuss this result and the formation of Herschel 36 A in the context of Core Accretion and Competitive Accretion models, which make different predictions regarding the statistic of the relative orbital inclinations. © 2022 The Author(s). Published by Oxford University Press on behalf of Royal Astronomical Society., RS and AA acknowledge financial support from the State Agency for Research of the Spanish MCIU through the ‘Center of Excellence Severo Ochoa’ award for the Instituto de Astrofísica de Andalucía (SEV-2017-0709). RS acknowledges financial support from national project PGC2018-095049-B-C21 (MCIU/AEI/FEDER, UE). AA acknowledges financial support from national project PID2020-117404GB-C21 (MCIU/AEI/FEDER, UE). JSB acknowledges the financial support from the ‘Visitor Scientist Program’ of the ‘Center of Excellence Severo Ochoa’ provided by the IAA-CSIC; and to the ‘ESO-Garching Visitor Program’ of the European Southern Observatory. This work presents results from the European Space Agency (ESA) space mission Gaia. Gaia data are being processed by the Gaia Data Processing and Analysis Consortium (DPAC). Funding for the DPAC is provided by national institutions, in particular the institutions participating in the Gaia MultiLateral Agreement (MLA).

Context. The young massive cluster Westerlund 1 offers the promise of a grand laboratory for the study of high-mass star evolution, but its basic parameters are still poorly known. Aims. In this paper, we aim at a better characterisation of the cluster by determining some basic kinematic properties and analysing the area surrounding the cluster and the population in its foreground. Methods. We have used Gaia early data release 3 (EDR3) data, together with spectra of a large sample of luminous stars in the field surrounding Westerlund 1, to explore the extent of the cluster. We carried out a non-parametric analysis of proper motions and membership determination. We investigated the reddening and proper motions of several dozen OB stars and red supergiants less than one degree away from Westerlund 1. Results. We identify a population of kinematic members of Westerlund 1 that largely includes the known spectroscopic members. From their EDR3 parallaxes, we derive a distance to the cluster of kpc. We analyse the extinction in this direction, finding that it increases by a large amount around 2.8 kpc, which in all likelihood is due to dark clouds associated with the Scutum-Crux arm. As a consequence, we hardly see any stars at distances comparable (or higher) than that of the cluster. The proper motions of Westerlund 1, however, are very similar to those of stars in the field surrounding it which are – almost without exception – less distant, but distinct. We find a second, astrometrically well-defined population in the foreground (d ≈ 2 kpc), centred ∼8′ away, which is likely connected to the possible open cluster BH 197. Westerlund 1 is very elongated, an effect that seems real and not driven by the very heavy extinction to the east and south. We find a low-density halo extending to distances up to 10′ from the cluster centre, mainly in the north-west quadrant. A few OB stars at larger distances from the cluster, most notably the luminous blue variable (LBV) MN48, share its proper motions, suggesting that Westerlund 1 has little or no peculiar motion with respect to the field population of the Norma arm. Despite this, we are unable to find any red supergiant that could belong to an extended population related to the cluster, although we observe several dozen such objects in the foreground, demonstrating the richness of the field population along this sightline. We find a substantial population of luminous OB members obscured by several more magnitudes of extinction than most known members. These objects, mostly located in the central region of the cluster, increase the population of OB supergiants by about 25%. © ESO 2022., This research is partially supported by the Spanish Government Ministerio de Ciencia, Innovación y Universidades under grants PGC2018-093 741-B-C21/C22 and PGC2018-095 049-B-C21/C22 (MICIU/AEI/FEDER, UE) and the Generalitat Valenciana under grant PROMETEO/2019/041. EJA acknowledges financial support from the State Agency for Research of the Spanish MCIU through the “Center of Excellence Severo Ochoa” award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709)”. RD acknowledges further support from the Spanish Government Ministerio de Ciencia e Innovación through grant SEV 2015-0548, and from the Canarian Agency for Research, Innovation and Information Society (ACIISI), of the Canary Islands Government, and the European Regional Development Fund (ERDF), under grant with reference ProID2 017 010 115. The AAOmega observations have been supported by the OPTICON project (observing proposal 2011A/014), which was funded by the European Commission under the Seventh Framework Programme (FP7). This research has made use of the Simbad, Vizier and Aladin services developed at the Centre de Données Astronomiques de Strasbourg, France. 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. The Gaia data are processed with the computer resources at Mare Nostrum and the technical support provided by BSC-CNS. In this work, we have made extensive use of TOPCAT (Taylor et al. 2005). We thank its author, and subsequent contributors, for the creation and development of this tool.

We have started an ambitious program to determine if the full diversity of extinction laws is real or if some of it is due to calibration or methodological issues. Here we start by analyzing the information on NIR extinction in a 2MASS stellar sample with good quality photometry and very red colours. We calculate the extinction at 1 $\mu$m, $A_1$, and the power-law exponent, $\alpha$ ($A_\lambda = A_1 \lambda^{-\alpha}$), for the 2MASS stars located in the extinction trajectory in the $H-K$ vs. $J-H$ plane expected for red giants with $A_1 > 5$ mag. We test the validity of the assumption about the nature of those stars, whether a single or multiple values of $\alpha$ are needed, and the spatial variations of the results. Most ($\sim$83%) of those stars can indeed be explained by high-extinction red giants and the rest is composed of extinguished AGB stars (mostly O-rich), blended sources, and smaller numbers of other objects, a contaminant fraction that can be reduced with the help of Gaia} DR2 data. Galactic red giants experience a NIR extinction with $\alpha\sim 2.27$ and an uncertainty of a few hundredths of a magnitude. There is no significant spread in $\alpha$ even though our sample is widely distributed and has a broad range of extinctions. Differences with previous results are ascribed to the treatment of non-linear photometric effects and/or the contaminant correction. Future research should concentrate in finding the correct functional form for the NIR extinction law. In the appendix we detail the treatment of non-linear photometric effects in the 2MASS bands., Comment: 13 pages, accepted for publication in MNRAS

Massive stars at low metallicity are among the main feedback agents in the early Universe and in present-day star forming galaxies. When in binaries, these stars are potential progenitors of gravitational-wave events. Knowledge of stellar masses is a prerequisite to understanding evolution and feedback of low-metallicity massive stars. Using abundant spectroscopic and photometric measurements of an outstandingly bright eclipsing binary, we compare its dynamic, spectroscopic, and evolutionary mass estimates and develop a binary evolution scenario. We comprehensively studied the eclipsing binary system, AzV 476, in the Small Magellanic Cloud. The light curve and radial velocities were analyzed to obtain the orbital parameters. The photometric and spectroscopic data in the UV and optical were analyzed using the Potsdam Wolf-Rayet model atmospheres. The obtained results are interpreted using binary-evolution tracks. AzV 476 consists of an O4IV-III((f))p primary and an O9.5:Vn secondary. Both components have similar current masses (~20 M$_{\odot}$) obtained from both the orbital and spectroscopic analysis. The wind mass-loss rate of log($\dot{M}$/(M$_{\odot}$/yr))=-6.2 of the primary is a factor of ten higher than a recent empirical prescription for single O stars in the SMC. Only close-binary evolution with mass transfer can reproduce the current stellar and orbital parameters. The binary evolutionary model reveals that the primary has lost about half of its initial mass and is already core helium burning. Our comprehensive analysis of AzV 476 yields a consistent set of parameters and suggests previous case B mass transfer. The derived stellar masses agree within their uncertainties. The moderate masses of AzV 476 underline the scarcity of bright massive stars in the SMC. The core helium burning nature of the primary indicates that stripped stars might be hidden among OB-type populations., Comment: 18 pages + 5 pages appendix, 17 figures, 8 tables. Accepted for publication in A&A

[Context] Dynamical interactions in young stellar clusters can eject massive stars early in their lives and significantly alter their mass functions. If all of the most massive stars are lost, we are left with an orphan cluster. [Aims] We study the Bermuda cluster (Villafranca O-014 NW), the most significant young stellar group in the North America and Pelican nebulae, and the massive stars that may have been ejected from it to test if it has been orphane-d. [Methods] We use Gaia EDR3 parallaxes and proper motions to search for walkaway and runaway stars in the vicinity of the North America and Pelican nebulae. The candidates are analyzed with a combination of spectroscopy and photometry to assess their nature and their trajectories are traced back in time to determine at what time they left the Bermuda cluster. [Results] We detect three ejection events, dubbed the Bajamar, Toronto, and HD 201 795 events, which expelled (a minimum of) five, two, and two systems, respectively, or six, three, and three stars if we count the individual components in spectroscopic and eclipsing binaries. The events took place 1.611 ± 0.011 Ma, 1.496 ± 0.044 Ma, and 1.905 ± 0.037 Ma ago, respectively, but our analysis is marginally consistent with the first two being simultaneous. We detect bow shocks in WISE images associated with four of the ejected systems; their orientation agrees with that of their relative proper motions with respect to the cluster. Combining the three events, the Bermuda cluster has lost over 200 M, including its three most massive stars, so it can be rightfully considered an orphan cluster. One consequence is that the present-day mass function of the cluster has been radically altered from its top-heavy initial value to one compatible with a Kroupa-like function. Another is that the cluster is currently expanding with a dynamical timescale consistent with the cause being the ejection events. A scenario in which the Bermuda cluster was formed in a conveyor-belt fashion over several hundreds of ka or even 1 Ma is consistent with all the observables., J. M. A. and M. P. G. acknowledge support from the Spanish Government Ministerio de Ciencia e Innovación through grant PGC2018-095 049-B-C22. R. H. B. acknowledges support from ANID FONDECYT Regular Project 1 211 903 and the ESAC visitors program. M. W. acknowledges support from the Spanish Government Ministerio de Ciencia e Innovación (MICI/FEDER, UE) through grant RTI2018-095 076-B-C21, and from the Instituto de Ciencias del Cosmos de la Univ. de Barcelona (ICCUB, Unidad de Excelencia “María de Maeztu”) through grant CEX2019-000 918-M.

Full list of authors: Gilmore, G.; Randich, S.; Worley, C. C.; Hourihane, A.; Gonneau, A.; Sacco, G. G.; Lewis, J. R.; Magrini, L.; Francois, P.; Jeffries, R. D.; Koposov, S. E.; Bragaglia, A.; Alfaro, E. J.; Allende Prieto, C.; Blomme, R.; Korn, A. J.; Lanzafame, A. C.; Pancino, E.; Recio-Blanco, A.; Smiljanic, R.; Van Eck, S.; Zwitter, T.; Bensby, T.; Flaccomio, E.; Irwin, M. J.; Franciosini, E.; Morbidelli, L.; Damiani, F.; Bonito, R.; Friel, E. D.; Vink, J. S.; Prisinzano, L.; Abbas, U.; Hatzidimitriou, D.; Held, E., V; Jordi, C.; Paunzen, E.; Spagna, A.; Jackson, R. J.; Maiz Apellaniz, J.; Asplund, M.; Bonifacio, P.; Feltzing, S.; Binney, J.; Drew, J.; Ferguson, A. M. N.; Micela, G.; Negueruela, I; Prusti, T.; Rix, H-W; Vallenari, A.; Bergemann, M.; Casey, A. R.; de Laverny, P.; Frasca, A.; Hill, V; Lind, K.; Sbordone, L.; Sousa, S. G.; Adibekyan, V; Caffau, E.; Daflon, S.; Feuillet, D. K.; Gebran, M.; Gonzalez Hernandez, J., I; Guiglion, G.; Herrero, A.; Lobel, A.; Montes, D.; Morel, T.; Ruchti, G.; Soubiran, C.; Tabernero, H. M.; Tautvaisiene, G.; Traven, G.; Valentini, M.; Van der Swaelmen, M.; Villanova, S.; Vazquez, C. Viscasillas; Bayo, A.; Biazzo, K.; Carraro, G.; Edvardsson, B.; Heiter, U.; Jofre, P.; Marconi, G.; Martayan, C.; Masseron, T.; Monaco, L.; Walton, N. A.; Zaggia, S.; Borsen-Koch, V. Aguirre; Alves, J.; Balaguer-Nunez, L.; Barklem, P. S.; Barrado, D.; Bellazzini, M.; Berlanas, S. R.; Binks, A. S.; Bressan, A.; Capuzzo-Dolcetta, R.; Casagrande, L.; Casamiquela, L.; Collins, R. S.; D'Orazi, V; Dantas, M. L. L.; Debattista, V. P.; Delgado-Mena, E.; Di Marcantonio, P.; Drazdauskas, A.; Evans, N. W.; Famaey, B.; Franchini, M.; Fremat, Y.; Fu, X.; Geisler, D.; Gerhard, O.; Solares, E. A. Gonzalez; Grebel, E. K.; Gutierrez Albarran, M. L.; Jimenez-Esteban, F.; Jonsson, H.; Khachaturyants, T.; Kordopatis, G.; Kos, J.; Lagarde, N.; Ludwig, H-G; Mahy, L.; Mapelli, M.; Marfil, E.; Martell, S. L.; Messina, S.; Miglio, A.; Minchev, I; Moitinho, A.; Montalban, J.; Monteiro, M. J. P. F. G.; Morossi, C.; Mowlavi, N.; Mucciarelli, A.; Murphy, D. N. A.; Nardetto, N.; Ortolani, S.; Paletou, F.; Palous, J.; Pickering, J. C.; Quirrenbach, A.; Fiorentin, P. Re; Read, J., I; Romano, D.; Ryde, N.; Sanna, N.; Santos, W.; Seabroke, G. M.; Spina, L.; Steinmetz, M.; Stonkute, E.; Sutorius, E.; Thevenin, F.; Tosi, M.; Tsantaki, M.; Wright, N.; Wyse, R. F. G.; Zoccali, M.; Zorec, J.; Zucker, D. B.-- This is an Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited., Context. The Gaia-ESO Public Spectroscopic Survey is an ambitious project designed to obtain astrophysical parameters and elemental abundances for 100 000 stars, including large representative samples of the stellar populations in the Galaxy, and a well-defined sample of 60 (plus 20 archive) open clusters. We provide internally consistent results calibrated on benchmark stars and star clusters, extending across a very wide range of abundances and ages. This provides a legacy data set of intrinsic value, and equally a large wide-ranging dataset that is of value for the homogenisation of other and future stellar surveys and Gaia's astrophysical parameters. Aims. This article provides an overview of the survey methodology, the scientific aims, and the implementation, including a description of the data processing for the GIRAFFE spectra. A companion paper introduces the survey results. Methods. Gaia-ESO aspires to quantify both random and systematic contributions to measurement uncertainties. Thus, all available spectroscopic analysis techniques are utilised, each spectrum being analysed by up to several different analysis pipelines, with considerable effort being made to homogenise and calibrate the resulting parameters. We describe here the sequence of activities up to delivery of processed data products to the ESO Science Archive Facility for open use. Results. The Gaia-ESO Survey obtained 202 000 spectra of 115 000 stars using 340 allocated VLT nights between December 2011 and January 2018 from GIRAFFE and UVES. © G. Gilmore et al. 2022. Conclusions. The full consistently reduced final data set of spectra was released through the ESO Science Archive Facility in late 2020, with the full astrophysical parameters sets following in 2022. A companion article reviews the survey implementation, scientific highlights, the open cluster survey, and data products., Based on data products from observations made with ESO Telescopes at the La Silla Paranal Observatory under programme ID 188.B-3002. These data products have been processed by the Cambridge Astronomy Survey Unit (CASU) at the Institute of Astronomy, University of Cambridge, and by the FLAMES/UVES reduction team at INAF/Osservatorio Astrofisico di Arcetri. Public access to the data products is via the ESO SAF, and the Gaia-ESO Survey Data Archive, prepared and hosted by the Wide Field Astronomy Unit, Institute for Astronomy, University of Edinburgh, which is funded by the UK Science and Technology Facilities Council. This work was partly supported by the European Union FP7 programme through ERC grant number 320360 and by the Leverhulme Trust through grant RPG-2012-541. We acknowledge the support from INAF and Ministero dell’ Istruzione, dell’ Università’ e della Ricerca (MIUR) in the form of the grant “Premiale VLT 2012”. The project presented here benefited in development from discussions held during the Gaia-ESO workshops and conferences supported by the ESF (European Science Foundation) through the GREAT Research Network Programme. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France. R.S. acknowledges support from the National Science Centre, Poland (2014/15/B/ST9/03981). This work was partly supported by the INAF grant for mainstream projects: “Enhancing the legacy of the Gaia-ESO Survey for open cluster science”. F.J.E. acknowledges financial support from the Spanish MINECO/FEDER through the grant AYA2017-84089 and MDM-2017-0737 at Centro de Astrobiología (CSIC-INTA), Unidad de Excelencia María de Maeztu, and from the European Union’s Horizon 2020 research and innovation programme under Grant Agreement no. 824064 through the ESCAPE – the European Science Cluster of Astronomy and Particle Physics ESFRI Research Infrastructures project. T.B. was funded by the “The New Milky Way” project grant from the Knut and Alice Wallenberg Foundation. S.R.B. acknowledges support by the Spanish Government under grants AYA2015-68012-C2-2-P and PGC2018-093741-B-C21/C22 (MICIU/AEI/FEDER, UE). W.J.S. acknowledges CAPES for a PhD studentship. J.M.A. acknowledges support from the Spanish Government Ministerio de Ciencia e Innovación through grants AYA2013-40611-P, AYA2016-75931-C2-2-P, and PGC2018-095049-B-C22. T.M. and others from STAR institute, Liege, Belgium are grateful to Belgian F.R.S.-FNRS for support, and are also indebted for an ESA/PRODEX Belspo contract related to the Gaia Data Processing and Analysis Consortium and for support through an ARC grant for Concerted Research Actions financed by the Federation Wallonie-Brussels.. This research has been partially supported by the following grants: MIUR Premiale “Gaia-ESO survey” (PI S. Randich), MIUR Premiale “MiTiC: Mining the Cosmos” (PI B. Garilli), the ASI-INAF contract 2014-049-R.O: “Realizzazione attività tecniche/scientifiche presso ASDC” (PI Angelo Antonelli), Fondazione Cassa di Risparmio di Firenze, progetto: “Know the star, know the planet” (PI E. Pancino), and Progetto Main Stream INAF: “Chemo-dynamics of globular clusters: the Gaia revolution” (PI E. Pancino). V.A. acknowledges the support from Fundação para a Ciência e Tecnologia (FCT) through Investigador FCT contract nr. IF/00650/2015/CP1273/CT0001. A.J.K. acknowledges support by the Swedish National Space Agency (SNSA). A.B. acknowledges support by ANID, – Millennium Science Initiative Program – NCN19_171, and FONDECYT regular 1190748. E.M. acknowledges financial support from the Spanish State Research Agency (AEI) through project MDM-2017-0737 Unidad de Excelencia “María de Maeztu” – Centro de Astrobiología (CSIC-INTA). T.Z. acknowledges financial support of the Slovenian Research Agency (research core funding no. P1-0188) and the European Space Agency (Prodex Experiment Arrangement No. C4000127986). P.J. acknowledges support FONDECYT Regular 1200703. The work of I.N. is partially supported by the Spanish Government Ministerio de Ciencia, Innovación y Universidades under grant PGC2018-093741-B-C21 (MICIU/AEI/FEDER, UE). Funding for this work has been provided by the ARC Future Fellowship FT160100402. C.A.P. acknowledges financial support from the Spanish Government through research grants MINECO AYA 2014-56359-P, MINECO AYA2017-86389-P, and MICINN PID2020-117493GB-I00. S.F. was supported by the grants 2011-5042 and 2016-03412 from the Swedish Research Council and the project grant “The New Milky Way” from the Knut and Alice Wallenberg Foundation. C.A.S.U. is supported through STFC grants: ST/H004157/1, ST/J00541X/1, ST/M007626/1, ST/N005805/1, ST/T003081/1. Work reported here benefited from support through the GREAT-ITN FP7 project Grant agreement ID: 264895. D.K.F. acknowledges funds from the Alexander von Humboldt Foundation in the framework of the Sofja Kovalevskaja Award endowed by the Federal Ministry of Education and Research and the grant 2016-03412 from the Swedish Research Council. A.H. acknowledges support from the Spanish Government Ministerio de Ciencia e Innovación and ERD Funds through grants PGC-2018-091 3741-B-C22 and CEX2019-000920-S. X.F. acknowledges the support of China Postdoctoral Science Foundation 2020M670023. M.L.L.D. acknowledges the Polish NCN grant number 2019/34/E/ST9/00133. Part of this work was funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – Project-ID 138713538 – SFB 881 (“The Milky Way System”, subproject A09). M.Z. acknowledges support from the National Agency for Research and Development (ANID) grants: FONDECYT Regular 1191505, Millennium Institute of Astrophysics ICN12-009, BASAL Center for Astrophysics and Associated Technologies AFB-170002. R.B. acknowledges support from the project PRIN-INAF 2019 “Spectroscopically Tracing the Disk Dispersal Evolution”. H.M.T. acknowledges financial support from the Agencia Estatal de Investigación of the Ministerio de Ciencia, Innovación y Universidades through projects PID2019-109522GB-C51,54/AEI/10.13039/501100011033, and the Centre of Excellence “María de Maeztu” award to Centro de Astrobiología (MDM-2017-0737). J.I.G.H. acknowledges financial support from the Spanish Ministry of Science and Innovation (MICINN) project AYA2017-86389-P, and also from the Spanish MICINN under 2013 Ramѳn y Cajal program RYC-2013-14875. V.P.D. is supported by STFC Consolidated grant ST/R000786/1. N.L. acknowledges financial support from “Programme National de Physique Stellaire” (PNPS) and the “Programme National Cosmology et Galaxies (PNCG)” of CNRS/INSU, France. A.R.C. is supported in part by the Australian Research Council through a Discovery Early Career Researcher Award (DE190100656). Parts of this research were supported by the Australian Research Council Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), through project number CE170100013. P.S.B. is supported by the Swedish Research Council through individual project grants with contract Nos. 2016-03765 and 2020-03404. A.M. acknowledges funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement no. 772293 – project ASTEROCHRONOMETRY). J.P. was supported by the project RVO: 67985815. E.D.M. acknowledges the support from FCT through the research grants UIDB/04434/2020 & UIDP/04434/2020 and through Investigador FCT contract IF/00849/2015/CP1273/CT0003. This work was (partially) supported by the Spanish Ministry of Science, Innovation and University (MICIU/FEDER, UE) through grant RTI2018-095076-B-C21, and the Institute of Cosmos Sciences University of Barcelona (ICCUB, Unidad de Excelencia 'María de Maeztu’) through grant CEX2019-000918-M. S.L.M. acknowledges the support of the UNSW Scientia Fellowship program and the Australian Research Council through Discovery Project grant DP180101791. GT acknowledges financial support of the Slovenian Research Agency (research core funding No. P1-0188) and the European Space Agency (Prodex Experiment Arrangement No. C4000127986). S.G.S acknowledges the support from FCT through Investigador FCT contract nr. CEECIND/00826/2018 and POPH/FSE (EC). H.G.L. acknowledges financial support by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – Project-ID 138713538 – SFB 881 (“The Milky Way System”, subproject A04). This work was (partially) supported by the Spanish Ministry of Science, Innovation and University (MICIU/FEDER, UE) through grant RTI2018-095076-B-C21, and the Institute of Cosmos Sciences University of Barcelona (ICCUB, Unidad de Excelencia ’María de Maeztu’) through grant CEX2019-000918-M. T.K. is supported by STFC Consolidated grant ST/R000786/1. M.V. acknowledges the support of the Deutsche Forschungsgemeinschaft (DFG, project number: 428473034). T.M. is supported by a grant from the Fondation ULB. We acknowledge financial support from the Universidad Complutense de Madrid (UCM) and by the Spanish Ministerio de Ciencia, Innovación y Universidades, Ministerio de Economía y Competitividad, from project AYA2016-79425-C3-1-P and PID2019-109522GB-C5[4]/AEI/10.13039/501100011033. U.H. acknowledges support from the Swedish National Space Agency (SNSA/Rymdstyrelsen). D.G. gratefully acknowledges support from the Chilean Centro de Excelencia en Astrofísica y Tecnologías Afines (CATA) BASAL grant AFB-170002. D.G. also acknowledges financial support from the Dirección de Investigación y Desarrollo de la Universidad de La Serena through the Programa de Incentivo a la Investigación de Académicos (PIA-DIDULS). A. Lobel acknowledges support in part by the Belgian Federal Science Policy Office under contract No. BR/143/A2/BRASS. We acknowledge financial support from the Universidad Complutense de Madrid (UCM) and by the Spanish Ministerio de Ciencia, Innovación y Universidades, Ministerio de Economía y Competitividad, from project AYA2016-79425-C3-1-P and PID2019-109522GB-C5[4]/AEI/10.13039/501100011033. A.M. acknowledges the support from the Portuguese Fundação para a Ciência e a Tecnologia (FCT) through the Portuguese Strategic Programme UID/FIS/00099/2019 for CENTRA. T.M. acknowledges financial support from the Spanish Ministry of Science and Innovation (MICINN) through the Spanish State Research Agency, under the Severo Ochoa Program 2020-2023 (CEX2019-000920-S). E.J.A. acknowledges funding from the State Agency for Research of the Spanish MCIU through the “Center of Excellence Severo Ochoa” award to the Instituto de Astrofísica de An- dalucía (SEV-2017-0709).

Full list of authors: Blomme, R.; Daflon, S.; Gebran, M.; Herrero, A.; Lobel, A.; Mahy, L.; Martins, F.; Morel, T.; Berlanas, S. R.; Blazère, A.; Frémat, Y.; Gosset, E.; Maíz Apellániz, J.; Santos, W.; Semaan, T.; Simón-Díaz, S.; Volpi, D.; Holgado, G.; Jiménez-Esteban, F.; Nieva, M. F.; Przybilla, N.; Gilmore, G.; Randich, S.; Negueruela, I.; Prusti, T.; Vallenari, A.; Alfaro, E. J.; Bensby, T.; Bragaglia, A.; Flaccomio, E.; Francois, P.; Korn, A. J.; Lanzafame, A.; Pancino, E.; Smiljanic, R.; Bergemann, M.; Carraro, G.; Franciosini, E.; Gonneau, A.; Heiter, U.; Hourihane, A.; Jofré, P.; Magrini, L.; Morbidelli, L.; Sacco, G. G.; Worley, C. C.; Zaggia, S., Context. The Gaia-ESO Survey (GES) is a large public spectroscopic survey that has collected, over a period of six years, spectra of ~105 stars. This survey provides not only the reduced spectra, but also the stellar parameters and abundances resulting from the analysis of the spectra. Aims. The GES dataflow is organised in 19 working groups. Working group 13 (WG13) is responsible for the spectral analysis of the hottest stars (O, B, and A type, with a formal cutoff of Teff > 7000 K) that were observed as part of GES. We present the procedures and techniques that have been applied to the reduced spectra in order to determine the stellar parameters and abundances of these stars. Methods. The procedure used was similar to that of other working groups in GES. A number of groups (called Nodes) each independently analyse the spectra via state-of-the-art techniques and codes. Specific for the analysis in WG13 was the large temperature range covered (Teff ≈ 7000–50 000 K), requiring the use of different analysis codes. Most Nodes could therefore only handle part of the data. Quality checks were applied to the results of these Nodes by comparing them to benchmark stars, and by comparing them to one another. For each star the Node values were then homogenised into a single result: the recommended parameters and abundances. Results. Eight Nodes each analysed part of the data. In total 17 693 spectra of 6462 stars were analysed, most of them in 37 open star clusters. The homogenisation led to stellar parameters for 5584 stars. Abundances were determined for a more limited number of stars. The elements studied are He, C, N, O, Ne, Mg, Al, Si, and Sc. Abundances for at least one of these elements were determined for 292 stars. Conclusions. The hot-star data analysed here, as well as the GES data in general, will be of considerable use in future studies of stellar evolution and open clusters. © ESO 2022., This paper is based on data products from observations made with ESO Telescopes at the La Silla Paranal Observatory under programmes ID 188.B-3002, 193.B-0936, and 197.B-1074. These data products have been processed by the Cambridge Astronomy Survey Unit (CASU) at the Institute of Astronomy, University of Cambridge, and by the FLAMES/UVES reduction team at INAF/Osservatorio Astrofisico di Arcetri. This work was partly supported by the European Union FP7 programme through ERC grant number 320360 and by the Leverhulme Trust through grant RPG-2012-541. We acknowledge the support from INAF and Ministero dell’ Istruzione, dell’ Université’ e della Ricerca (MIUR) in the form of the grant ‘Premiale VLT 2012’. The results presented here benefitted from discussions held during the Gaia-ESO workshops and conferences supported by the ESF (European Science Foundation) through the GREAT Research Network Programme. We are grateful to E. Bertone for making their high-resolution models available, as well as to R. Sordo for making the high-resolution Munari models available. This research has made use of the WEBDA database, originally developed by J.-C. Mermilliod, and now operated at the Department of Theoretical Physics and Astrophysics of the Masaryk University. This research has also made use of the SIMBAD database, operated at CDS, Strasbourg, France. A.H. acknowledges support by the Spanish Ministerio de Ciencia e Innovación through grants PGC2018-091, 3741-B-C22 and CEX-2019-000920-S. The work of S.R.B, I.N., and S.S.D. is partially supported by the Spanish Government Ministerio de Economía y Competitivad (MINECO/FEDER) under grants PGC2018-093741-B-C21/C22 (MICIU/AEI/FEDER, UE). A.L. acknowledges that his research has been funded in part by the Belgian Federal Science Policy Office under contract No. BR/143/A2/BRASS. The Liège Node is grateful to Belgian F.R.S.-FNRS for various supports. They are also indebted for an ESA/PRODEX Belspo contract related to the Gaia Data Processing and Analysis Consortium and for support through an ARC grant for Concerted Research Actions financed by the Federation Wallonie-Brussels. J.M.A. acknowledges support from the Spanish Government Ministerio de Ciencia through grants AYA2016-75 931-C2-2-P and PGC2018-095 049-B-C22. W.S. acknowledges CAPES for a Ph.D. Fellowship. G.H. acknowledges that this research has been partially funded by the Canarian Agency for Economy, Knowledge, and Employment and the European Regional Development Fund (ERDF/EU), under grant with reference ProID2020010016. F.J.E. acknowledges financial support from the Spanish MINECO/FEDER through the grant AYA2017-84089 and MDM-2017-0737 at Centro de Astrobiología (CSIC-INTA), Unidad de Excelencia María de Maeztu, and from the European Union’s Horizon 2020 research and innovation programme under Grant Agreement No. 824064 through the ESCAPE - The European Science Cluster of Astronomy & Particle Physics ESFRI Research Infrastructures project. E.J.A acknowledges funding from the State Agency for Research of the Spanish MCIU through the ‘Center of Excellence Severo Ochoa’ award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709) and from MCIU grant PGC2018-095049-B-C21. T.B. acknowledges financial support by grant No. 2018-04857 from the Swedish Research Council. A.J.K. and U.H. acknowledge support from the Swedish National Space Agency (SNSA). R.S. acknowledges support from the National Science Centre (2014/15/B/ST9/03981).

Two main scenarios have been proposed for origin of massive runaway stars -- dynamical ejection or release from a binary at the first core collapse -- but their relative contribution remains debated. Using two large spectroscopic campaigns towards massive stars in 30 Doradus, we aim to provide observational constraints on the properties of the O-type runaway population in the most massive active star-forming region in the Local group. We use RV measurements of the O-type star populations in 30 Doradus obtained by the VLT-FLAMES Tarantula Survey and the Tarantula Massive Binary Monitoring to identify single and binary O-type runaways. We discuss their rotational properties and qualitatively compare observations with expectations of ejection scenarios. We identify 23 single and one binary O-type runaway objects, most of them outside the main star-forming regions in 30 Doradus. We find an overabundance of rapid rotators (vsini > 200km/s) among the runaway population, providing an explanation of the overabundance of rapidly rotating stars in the 30 Doradus field. Considerations of the projected rotation rates and runaway line-of-sight (los) velocities reveal a conspicuous absence of rapidly rotating (vsini > 210k/ms), fast moving (v_{los} > 60km/s) runaways, and suggest the presence of two different populations of runaway stars: a population of rapidly-spinning but slowly moving runaways and a population of fast moving but slowly rotating ones. These are detected with a ratio close to 2:1 in our sample. We argue that slowly moving but rapidly spinning runaways result from binary ejections, while rapidly moving but slowly spinning runaways could result from dynamical ejections. Given that detection biases will more strongly impact the slow-moving population, our results suggest that the binary evolution scenario dominates the current massive runaway population in 30 Doradus., Accepted for publication in A&A Letters; 9 pages, 5 figures

Context. Massive stars are a key element for understanding the chemical and dynamical evolution of galaxies. Stellar evolution is conditioned by many factors: Rotation, mass loss, and interaction with other objects are the most important ones for massive stars. During the first evolutionary stages of stars with initial masses (i.e., MZAMS) in the MZAMS ∼ 18–70 M⊙ range, they are of spectral type O. Given that stars in this mass range spend roughly 90% of their lifetime as O-type stars, establishing the multiplicity frequency and binary properties of O-type stars is crucial for many fields of modern astrophysics. Aims. The aim of the MONOS project is to collect information to study northern Galactic O-type spectroscopic binaries. In this second paper, we tackle the study of the 35 single-line spectroscopic binary (SB1) systems identified in the previous paper of the series, analyze our data, and review the literature on the orbits of the systems. Methods. We have measured ∼4500 radial velocities for a selection of diagnostic lines for the ∼700 spectra of the studied systems in our database, for which we have used two different methods: a Gaussian fit for several lines per object and cross-correlation with synthetic spectra computed with the FASTWIND stellar atmospheric code. We have also explored the photometric data delivered by the TESS mission to analyze the light curve (LC) of the systems, extracting 31 of them. We have explored the possible periods with the Lomb-Scargle method and, whenever possible, calculated the orbital solutions using the SBOP and GBART codes. For those systems in which an improved solution was possible, we merged our radial velocities with those in the literature and calculated a combined solution. Results. As a result of this work, of the 35 SB1 systems identified in our first paper we have confirmed 21 systems as SB1 with good orbits, discarded the binary nature of six stars (9 Sge, HD 192 281, HDE 229 232 AB, 68 Cyg, HD 108, and α Cam), and left six stars as inconclusive due to a lack of data. The remaining two stars are 15 Mon Aa, which has been classified as SB2, and Cyg OB2-22 C, for which we find evidence that it is most likely a triple system where the O star is orbiting an eclipsing SB1. We have also recalculated 20 new orbital solutions, including the first spectroscopic orbital solution for V747 Cep. For Cyg OB2-22 C, we have obtained new ephemerides but no new orbit. E.T.P., I.N.D. and S.S.D. acknowledge support from the Spanish Government Ministerio de Ciencia through grant AYA2015-68 012-C2-1/2-P. E.T.P. and J.M.A. acknowledge support from the Spanish Government Ministerio de Ciencia through grants AYA2016-75 931-C2-2-P and PGC2018-095 049-B-C22. R.H.B. acknowledges support from the ESAC Faculty Council Visitor Program and ANID FONDECYT Project 1211903. G.H. acknowledges support from the Spanish Government Ministerio de Ciencia through grant PGC2018-95049-B-CC22 and ESA Contract No. 4000-126507/19/ES/CM. This research has made use of the SIMBAD database, operated at CDS (Strasbourg, France), NASA’s Astrophysics Data System Bibliographic Services and the Python programming language (Python Software Foundation), Lightkurve, a Python package for Kepler and TESS data analysis (Lightkurve Collaboration 2018).

[Context]: The Vela X-1 system is one of the best-studied X-ray binaries because it was detected early, has persistent X-ray emission, and a rich phenomenology at many wavelengths. The system is frequently quoted as the archetype of wind-accreting high-mass X-ray binaries, and its parameters are referred to as typical examples. Specific values for these parameters have frequently been used in subsequent studies, however, without full consideration of alternatives in the literature, even more so when results from one field of astronomy (e.g., stellar wind parameters) are used in another (e.g., X-ray astronomy). The issues and considerations discussed here for this specific, very well-known example will apply to various other X-ray binaries and to the study of their physics., [Aims]: We provide a robust compilation and synthesis of the accumulated knowledge about Vela X-1 as a solid baseline for future studies, adding new information where available. Because this overview is targeted at a broader readership, we include more background information on the physics of the system and on methods than is usually done. We also attempt to identify specific avenues of future research that could help to clarify open questions or determine certain parameters better than is currently possible., [Methods]: We explore the vast literature for Vela X-1 and on modeling efforts based on this system or close analogs. We describe the evolution of our knowledge of the system over the decades and provide overview information on the essential parameters. We also add information derived from public data or catalogs to the data taken from the literature, especially data from the Gaia EDR3 release., [Results]: We derive an updated distance to Vela X-1 and update the spectral classification for HD 77518. At least around periastron, the supergiant star may be very close to filling its Roche lobe. Constraints on the clumpiness of the stellar wind from the supergiant star have improved, but discrepancies persist. The orbit is in general very well determined, but a slight difference exists between the latest ephemerides. The orbital inclination remains the least certain factor and contributes significantly to the uncertainty in the neutron star mass. Estimates for the stellar wind terminal velocity and acceleration law have evolved strongly toward lower velocities over the years. Recent results with wind velocities at the orbital distance in the range of or lower than the orbital velocity of the neutron star support the idea of transient wind-captured disks around the neutron star magnetosphere, for which observational and theoretical indications have emerged. Hydrodynamic models and observations are consistent with an accretion wake trailing the neutron star., [Conclusions]: With its extremely rich multiwavelength observational data and wealth of related theoretical studies, Vela X-1 is an excellent laboratory for exploring the physics of accreting X-ray binaries, especially in high-mass systems. Nevertheless, much room remains to improve the accumulated knowledge. On the observational side, well-coordinated multiwavelength observations and observing campaigns addressing the intrinsic variability are required. New opportunities will arise through new instrumentation, from optical and near-infrared interferometry to the upcoming X-ray calorimeters and X-ray polarimeters. Improved models of the stellar wind and flow of matter should account for the non-negligible effect of the orbital eccentricity and the nonspherical shape of HD 77581. There is a need for realistic multidimensional models of radiative transfer in the UV and X-rays in order to better understand the wind acceleration and effect of ionization, but these models remain very challenging. Improved magnetohydrodynamic models covering a wide range of scales are required to improve our understanding of the plasma-magnetosphere coupling, and they are thus a key factor for understanding the variability of the X-ray flux and the torques applied to the neutron star. A full characterization of the X-ray emission from the accretion column remains another so far unsolved challenge., IEM has received funding from the Research Foundation Flanders (FWO), from the European Union’s Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement No 665501 and from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No 863412). SMN acknowledges funding by the Spanish Ministry MCIU under project RTI2018-096686-B-C21 (MCIU/AEI/FEDER, UE), co-funded by FEDER funds and by the Unidad de Excelencia María de Maeztu, ref. MDM-2017-0765. VG is supported through the Margarete von Wrangell fellowship by the ESF and the Ministry of Science, Research and the Arts Baden-Württemberg. IEM and PK acknowledge support from the ESA/ESAC Faculty Visiting Scientist Programme. JMA acknowledges support from the Spanish Government Ministerio de Ciencia through grant PGC2018-095049-B-C22. FJE acknowledges support from ESCAPE – The European Science Cluster of Astronomy & Particle Physics ESFRI Research Infrastructures, which received funding from the European Union’s Horizon 2020 research and innovation programme under Grant Agreement no. 824064. JvdE is supported by a Lee Hysan Junior Research Fellowship awarded by St. Hilda’s College.

LB-1 has variously been proposed as either an X-ray dim B-type star plus black hole (B+BH) binary, or a Be star plus an inflated stripped star (Be+Bstr) binary. The Space Telescope Imaging Spectrograph (STIS) on board HST was used to obtain a flux-calibrated spectrum that is compared with non-LTE spectral energy distributions (SED) and line profiles for the proposed models. The Hubble data, together with the Gaia EDR3 parallax, provide tight constraints on the properties and stellar luminosities of the system. In the case of the Be+Bstr model we adopt the published flux ratio for the Be and Bstr stars, re-determine the T$_{eff}$ of the Bstr using the silicon ionization balance, and infer Teff for the Be star from the fit to the SED. We derive stellar parameters consistent with previous results, but with greater precision enabled by the Hubble SED. While the Be+Bstr model is a better fit to the HeI lines and cores of the Balmer lines in the optical, the B+BH model provides a better fit to the Si iv resonance lines in the UV. The analysis also implies that the Bstr star has roughly twice solar silicon abundance, difficult to reconcile with a stripped star origin. The Be star on the other hand has a rather low luminosity, and a spectroscopic mass inconsistent with its possible dynamical mass. The fit to the UV can be significantly improved by reducing the T$_{eff}$ and radius of the Be star, though at the expense of leading to a different mass ratio. In the B+BH model, the single B-type spectrum is a good match to the UV spectrum. Adopting a mass ratio of 5.1$\pm$0.1 (Liu et al. 2020) implies a BH mass of $\sim$21$^{+9}_{-8}M_{\odot}$., Comment: 8 pages, 7 figures, 1 online table In press with Astronomy & Astrophysics

We cross-match the Alma catalog of OB stars with Gaia DR2 astrometry and photometry as a first step towards producing a clean sample of massive stars in the solar neighbourhood with a high degree of completeness. We analyze the resulting colour-absolute magnitude diagram to divide our sample into categories and compare extinction estimates from two sources, finding problems with both of them. The distances obtained with three different priors are found to have few differences among them, indicating that Gaia DR2 distances are robust. An analysis of the 3-D distribution of massive stars in the solar neighbourhood is presented. We show that a kinematically distinct structure we dub the Cepheus spur extends from the Orion-Cygnus spiral arm towards the Perseus arm and is located above the Galactic mid-plane, likely being related to the recently discovered Radcliffe wave. We propose that this corrugation pattern in the Galactic disk may be responsible for the recent enhanced star formation at its crests and troughs. We also discuss our plans to extend this work in the immediate future., Comment: Accepted for publication in MNRAS

The GALANTE optical photometric survey is observing the northern Galactic plane and some adjacent regions using seven narrow- and intermediate-filters, covering a total of 1618 deg2. The survey has been designed with multiple exposure times and at least two different air masses per field to maximize its photometric dynamic range, comparable to that of Gaia, and ensure the accuracy of its photometric calibration. The goal is to reach at least 1 per cent accuracy and precision in the seven bands for all stars brighter than AB magnitude 17 while detecting fainter stars with lower values of the signal-to-noise ratio. The main purposes of GALANTE are the identification and study of extinguished O+B+WR stars, the derivation of their extinction characteristics, and the cataloguing of F and G stars in the solar neighbourhood. Its data will be also used for a variety of other stellar studies and to generate a high-resolution continuum-free map of the Hα emission in the Galactic plane. We describe the techniques and the pipeline that are being used to process the data, including the basis of an innovative calibration system based on Gaia DR2 and 2MASS photometry. © 2021 The Author(s)., JMA, GH, and HGE acknowledge support from the Spanish Government Ministerio de Ciencia through grant PGC2018-095 049-B-C22. EJA and AL-G acknowledge support from the Spanish Government Ministerio de Ciencia through grant PGC2018-095 049-B-C21 and from the State Agency for Research of the Spanish Government Ministerio de Ciencia through the ‘Center of Excellence Severo Ochoa’ award for the Instituto de Astrofísica de Andalucía (SEV-2017-0709). The GALANTE team at the Centro de Astrobiología acknowledges financial support from the Science and Operations Department of the European Space Agency – Contract Number 4 000 126 507. RHB acknowledges support from the ESAC Faculty Visitor Program. Funding for OAJ, UPAD, and CEFCA has been provided by the Governments of Spain and Aragón through the Fondo de Inversiones de Teruel; the Aragonese Government through the Research Groups E96, E103, and E16_17R; the Spanish Government Ministerio de Ciencia, Innovación y Universidades (MCIU/AEI/FEDER, UE) with grantPGC2018-097 585-B-C21; the Spanish Ministerio de Economía y Competitividad (MINECO/FEDER, UE) with grants AYA2015-66 211-C2-1-P, AYA2015-66 211-C2-2, AYA2012-30 789, and ICTS-2009-14; and the European FEDER funding with grants FCDD10-4E-867 and FCDD13-4E-2685. PRTC acknowledges financial support from Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) process number 2018/05 392-8 and CNPq process number 310 041/2018-0. The main source of data for this paper is the JAST80 telescope at the Observatorio Astrofísico de Javalambre, Teruel, Spain (owned, managed, and operated by the Centro de Estudios de Física del Cosmos de Aragón). As detailed in this paper, the initial part of the reduction and astrometric calibration was done by the OAJ Data Processing and Archiving Unit (UPAD) and the main part of the reduction and calibration was done at the Centro de Astrobiología. In addition to the JAST80 data, two other types of data products have been used as sources for the catalogue generation and calibration: [1] 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. [2] From the Two Micron All Sky Survey (2MASS), which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation, where the word ‘National’ refers to the United States of America. No optically reflecting surface with a size larger than the height of the tallest of the authors was required to write this paper.

CONTEXT. Gaia EDR3 has produced parallaxes for 1.468x10^9 sources but there are calibration issues that require corrections to the published values and uncertainties. AIMS. We want to characterize the behavior of the uncertainties of the Gaia EDR3 parallaxes. We also aim to provide a procedure for the calculation of distances to stars and stellar clusters. METHODS. We reanalyze some of the data in the calibration papers for QSO and LMC parallaxes and combine those results with measurements for six bright GCs. We calculate the angular covariance of EDR3 parallaxes at small separations based on the LMC results and combine it with the results for larger angles using QSOs to obtain an analytical formula for the angular covariance over the whole sky. The results for the six GCs are used to validate the parallax bias correction as a function of magnitude, color, and ecliptic latitude and to determine the constant used to convert internal uncertainties to external ones. RESULTS. The angular covariance at zero separation is 106.2 muas^2, yielding a minimum uncertainty for EDR3 parallaxes of 10.3 muas for individual stars. That value can be only slightly reduced for GCs after considering the behavior of the angular covariance of the parallaxes for small separations. The Lindegren et al. parallax bias correction works quite well, except for the brighter magnitudes, suggesting improvements may be possible there. The value of k is 1.1-1.7 and depends on G. Stars with moderately large values of RUWE can still provide useful parallaxes albeit with larger values of k. We give accurate and precise Gaia EDR3 distances to the six GCs and for the specific case of 47 Tuc we are able to beat the angular covariance limit and derive a high-precision distance of 4.53+-0.06 kpc. Finally, a recipe for the derivation of distances to stars and stellar clusters using Gaia EDR3 parallaxes is given. [ABRIDGED], Version accepted in A&A

Context. This is the second paper of a series on Galactic OB groups that uses astrometric and photometric data from Gaia and spectral classifications from the Galactic O-Star Spectroscopic Survey (GOSSS) and the Library of Libraries of Massive-star high-Resolution spectra (LiLiMaRlin). The previous paper was based on the second Gaia data release (DR2). Since then, the early third Gaia data release (EDR3) has appeared with new astrometry and photometry. Aims. The two aims of this paper are to revise the results for the sample from Paper I using Gaia EDR3 data and to expand the sample of analyzed stellar groups to 26, from Villafranca O-001 to Villafranca O-026. Methods. We used GOSSS to select Galactic stellar groups with O stars and an updated version of the method in Paper 0 of this series, combining Gaia EDR3 G + GBP + GRP photometry, positions, proper motions, and parallaxes to assign memberships and measure distances. We present 99 spectra from GOSSS and 32 from LiLiMaRlin for stars in the analyzed groups or in their foreground. Results. We derived distances to the 26 stellar groups with unprecedented precision and accuracy, with total (random plus systematic) uncertainties lower than 1% for distances within 1 kpc and of ∼3% around 3 kpc, which are values almost four times better than for Gaia DR2. We provide homogeneous spectral types for 110 stars and correct a number of errors in the literature, especially for objects in Villafranca O-023 (Orion nebula cluster). For each group, we discuss its membership and present possible runaway and walkaway stars. At least two of the studied groups, Villafranca O-O12 S in NGC 2467 and Villafranca O-014 NW in the North America nebula, are orphan clusters in which the most massive stars have been ejected by dynamical interactions, leaving objects with a capped mass function. The existence of such clusters has important consequences for the study of the initial mass function (IMF), the distribution of supernova explosions across the Galaxy, and the population and dynamics of isolated compact objects. We fit pre-main-sequence (PMS) isochrones to the color-magnitude diagrams (CMDs) of four clusters to derive ages of 2.0 ± 0.5 Ma for Villafranca O-026 (σ Orionis cluster), 4 ± 2 Ma for Villafranca O-016 (NGC 2264), 5.0 ± 0.5 Ma for Villafranca O-021 (NGC 2362), and 8 ± 2 Ma for Villafranca O-024 (γ Velorum cluster). © ESO 2022., J.M.A., R.F.A., M.P.G., P.C.B., and A.S. acknowledge support from the Spanish Government Ministerio de Ciencia e Innovación through grant PGC2018-095 049-B-C22. R. H. B. acknowledges support from ANID FONDECYT Regular Project 1 211 903 and the ESAC visitors program. E.J.A. acknowledges support from the State Agency for Research of the Spanish Government Ministerio de Ciencia e Innovación through the “Center of Excellence Severo Ochoa” award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709) and through grant PGC2018-095 049-B-C21 This work has made use of data from the European Space Agency (ESA) mission Gaia, processed by the Gaia Data Processing and Analysis Consortium (DPAC). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. The Gaia data is processed with the computer resources at Mare Nostrum and the technical support provided by BSC-CNS. This research has made extensive use of the SIMBAD and VizieR databases, operated at CDS, Strasbourg, France.

CONTEXT. The Gaia EDR3 parallaxes constitute the most detailed and accurate dataset that can be currently used to determine stellar distances in the solar neighborhood. Nevertheless, there is still room for improvement in their calibration and systematic effects can be further reduced in some circumstances. AIMS. The aim of this paper is to determine an improved Gaia EDR3 parallax bias as a function of magnitude, color, and ecliptic latitude using a single method applied to stars in open clusters, globular clusters, the LMC, and the SMC. METHODS. I study the behavior of the residuals or differences between the individual (stellar) parallaxes and the group parallaxes, which are assumed constant for the corresponding cluster or galaxy. This is done by first applying the Lindegren et al. (2021) zero point and then calculating a new zero point from the residuals of the first analysis. RESULTS. The Lindegren zero point shows very small residuals as a function of magnitude between individual and group parallaxes for G < 13 but significant ones for brighter stars, especially blue ones. The new zero point reduces those residuals, especially in the 9.2 < G < 13 range. The k factor that is used to convert from catalog parallax uncertainties to external uncertainties is small (1.1-1.7) for 9.2 < G < 11 and G > 13, intermediate (1.7-2.0) for 11 < G < 13, and large (>2.0) for G, Accepted in A&A. New version is the accepted one. Figures 2, 6, and 7 may appear fuzzy using Preview in Mac OS, use Adobe Acrobat instead

We have detected a broad interstellar absorption band centred close to 7700 \r{A} and with a FWHM of 176.6$\pm$3.9 \r{A}. This is the first such absorption band detected in the optical range and is significantly wider than the numerous diffuse interstellar bands (DIBs). It remained undiscovered until now because it is partially hidden behind the A telluric band produced by O$_2$. The band was discovered using STIS@HST spectra and later detected in a large sample of stars of diverse type (OB stars, BA supergiants, red giants) using further STIS and ground-based spectroscopy. The EW of the band is measured and compared with our extinction and K I $\lambda\lambda$7667.021,7701.093 measurements for the same sample. The carrier is ubiquitous in the diffuse and translucent Galactic ISM but is depleted in the environment around OB stars. In particular, it appears to be absent or nearly so in sightlines rich in molecular carbon. This behaviour is similar to that of the $\sigma$-type DIBs, which originate in the low/intermediate-density UV-exposed ISM but are depleted in the high-density UV-shielded molecular clouds. We also present an update on our previous work on the relationship between $E(4405-5495)$ and $R_{5495}$ and incorporate our results into a general model of the ISM., Comment: 17 pages, accepted for publication in MNRAS

We conducted an observational campaign towards one of the most massive and luminous colliding wind binaries in the Galaxy, HD~93129A, close to its periastron passage in 2018. During this time the source was predicted to be in its maximum of high-energy emission. Here we present our data analysis from the X-ray satellites \textit{Chandra} and \textit{NuSTAR} and the $\gamma$-ray satellite \textit{AGILE}. High-energy emission coincident with HD~93129A was detected in the X-ray band up to $\sim$18~keV, whereas in the $\gamma$-ray band only upper limits were obtained. We interpret the derived fluxes using a non-thermal radiative model for the wind-collision region. We establish a conservative upper limit for the fraction of the wind kinetic power that is converted into relativistic electron acceleration, $f_\mathrm{NT,e} < 0.02$. In addition, we set a lower limit for the magnetic field in the wind-collision region as $B_\mathrm{WCR} > 0.3$~G. We also argue a putative interpretation of the emission from which we estimate $f_\mathrm{NT,e} \approx 0.006$ and $B_\mathrm{WCR} \approx 0.5$~G. We conclude that multi-wavelength, dedicated observing campaigns during carefully selected epochs are a powerful tool for characterising the relativistic particle content and magnetic field intensity in colliding wind binaries., Comment: 11 pages, 10 figures, accepted for publication in MNRAS

HD 54662 AB is one of the three O + OB binaries known so far with orbital period longer than 1000 d, offering the opportunity to test scenarios of massive star formation and models of single stellar evolution. Here, we present a detailed study of this system based on new high-resolution spectra and data. A disentangling method is used to recover the individual spectra of the primary and secondary components, which are classified as O6.5 V(n)z and O7.5 Vz, respectively. Combining radial velocity measurements and astrometric data, a new absolute orbit with a period of 2113 ± 9 d and an eccentricity of 0.062 ± 0.008 is determined, confirming previous findings. However, absolute masses of 23.8 ± 1.1 M¿ for the primary and 20.3 ± 1.1 M¿ for the secondary are obtained, differing from previous determinations but in reasonable agreement with the spectral types of the stars. Primary and secondary components show remarkably different projected rotational velocities (160 and ¿40 $\rm km\, s^{-1}$, respectively), which is probably related to the formation process of the binary. Contrary to previously interpretations, the star with broader spectral features is the most massive object in the system. Stellar and wind parameters of both stars are derived through quantitative spectroscopic analysis of the disentangled spectra using fastwind models, and they are consistent with the current calibrations for O-type stars. Evolutionary masses and ages are also computed with the bonnsai tool. Ages below 2.5 Ma are obtained, in agreement with the youth expected from their Vz nature., With funding from the Spanish government through the "María de Maeztu Unit of Excellence" accreditation (MDM-2017-0737)

This article is an erratum for Astronomy and Astrophysics 626: A20 (2019).-- http://hdl.handle.net/10261/201335, No abstract

Context. The spectral classifications of the Galactic O-Star Spectroscopic Survey (GOSSS) and the astrometric and photometric data from Gaia have significantly improved our ability to measure distances and determine memberships of stellar groups (clusters, associations, or parts thereof) with OB stars. In the near future, the situation will be further improved thanks to subsequent Gaia data releases and new photometric and spectroscopic surveys. Aims. We initiated a program to identify and determine the membership of Galactic stellar groups with OB stars and measure distances to them. Given the data currently available, we started with the identification and distance determinations of groups with O stars. In this paper, we concentrate on groups that contain stars with the earliest spectral subtypes. Methods. We used GOSSS to select Galactic stellar groups with O2-O3.5 stars and the method described in paper 0 of this series, which combines Gaia DR2 G + GBP + GRP photometry, positions, proper motions, and parallaxes to assign robust memberships and measure distances. We also included Collinder 419 and NGC 2264, the clusters cited in that paper, to generate our first list of 16 O-type Galactic stellar groups. Results. We derived distances, determined the membership, and analyzed the structure of sixteen Galactic stellar groups with O stars, Villafranca O-001 to Villafranca O-016, including the fourteen groups with the earliest-O-type optically accessible stars known in the Milky Way. We compared our distance with previous results from the literature and establish that the best consistency is with (the small number of) VLBI parallaxes and the worst is with kinematic distances. Our results indicate that very massive stars can form in relatively low-mass clusters or even in near-isolation, as is the case for the Bajamar star in the North America nebula. This lends support to the hierarchical scenario of star formation, where some stars are born in well-defined bound clusters but others are born in associations that are unbound from the beginning: groups of newborn stars come in many shapes and sizes. We propose that HD 64 568 and HD 64 315 AB could have been ejected simultaneously from Haffner 18 (Villafranca O-012 S). Our results are consistent with a difference of ≈20 μas in the Gaia DR2 parallax zero point between bright and faint stars. © ESO 2020., J.M.A. and A.S. acknowledge support from the Spanish Government Ministerio de Ciencia through grants AYA2016-75 931-C2-2-P and PGC2018-095 049-B-C22. R.H.B. acknowledges support from DIDULS Project 18 143 and the ESAC visitors program. 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. Additionally, this paper includes data obtained with the MPG/ESO 2.2m Telescope at the Observatorio de La Silla, Chile; the 2.5m du Pont Telescope at the Observatorio de Las Campanas, Chile; the 10m HobbyEberly Telescope at McDonald Observatory, Texas, USA; the 4.2m William Herschel Telescope and the 10.4m Gran Telescopio Canarias at the Observatorio del Roque de los Muchachos, La Palma, Spain; and the 2.2m Telescope at the Centro Astronomico Hispano Andaluz, Almeria, Spain. We thank the sta ff at those observatories for their support. This research has made extensive use of the SIMBAD and VizieR databases, operated at CDS, Strasbourg, France.

Context. Many massive stars have nearby companions. These hamper a characterization of massive stars through spectroscopy.Aims. We continue to obtain spatially resolved spectroscopy of close massive visual binaries to derive their spectral types.Methods. We used the lucky spectroscopy technique to obtain a large number of short long-slit spectroscopic exposures of 19 close visual binaries under good seeing conditions. We selected those with the best characteristics, extracted the spectra using multiple-profile fitting, and combined the results to derive spatially separated spectra. The results were analyzed in combination with data from lucky imaging, regular intermediate-resolution single-order spectroscopy, and echelle high-resolution spectroscopy.Results. The new application of lucky spectroscopy has allowed us (among other results) to [a] spatially disentangle two O stars (FN CMa B and 6 Cas B) with brighter BA supergiant companions for the first time; [b] determine that two B stars (alpha Sco B and HD 164 492 B) with close and more massive companions are fast rotators (in the second case, solving a case of mistaken identity); [c] extend the technique to cases with extreme magnitude differences (the previous two cases plus CS Cam A,B), shorter separations (HD 193 443 A,B), and fainter primary magnitudes down to B=11 (HD 219 460 A,B); [d] spatially disentangle the spectra of stars with companions as diverse as an A supergiant (6 Cas A), a Wolf-Rayet star (HD 219 460 B = WR 157), and an M supergiant (alpha Sco A); [e] discover the unexpected identity of some targets such as two previously unknown bright O stars (HD 51 756 B and BD +60 544) and a new member of the rare OC category (HD 8768 A); and [f] identify and classify (in some cases for the first time) which of the components of four visual binaries (sigma Ori, HD 219 460, HD 194 649, and HD 191 201) is a double-lined spectroscopic binary. For another seven systems (FN CMa, sigma Sco, HD 51 756, HD 218 195, HD 17 520, HD 24 431, and HD 164 492), we detect signs of spectroscopic binarity using high-spectral-resolution spectroscopy. We also determine the limits of the technique. © ESO 2021., J.M.A., C.F., A.S., M.P.G., and G.H. acknowledge support from the Spanish Government Ministerio de Ciencia through grant PGC2018-095049-B-C22. R.H.B. acknowledges support from the ESAC Faculty Visitor Program. I.N. and S.S.-D. acknowledge support from the Spanish Government Ministerio de Ciencia through grant PGC2018-093741-B-C21/22 (MICIU/AEI/FEDER, UE). S.S.-D. also acknowledges funding from the Spanish Government Ministerio de Ciencia through grants SEV 2015-0548 and CEX2019-000920-S, and from the Canarian Agency for Research, Innovation and Information Society (ACIISI), of the Canary Islands Government, and the European Regional Development Fund (ERDF), under grant with reference ProID2017010115. This paper is based on (a) lucky (and regular longslit) spectroscopy obtained with the 4.2m William Herschel Telescope (WHT) at the Observatorio del Roque de los Muchachos (ORM) on the island of La Palma, Spain; (b) lucky imaging obtained with the 2.2m Telescope at the Centro Astronomico Hispano en Andalucia (CAHA) in Almeria, Spain; (c) IFU spectroscopy obtained with the 2 m Liverpool Telescope (LT) at the Observatorio del Roque de los Muchachos (ORM) on the island of La Palma, Spain as part of GOSSS; (d) long-slit spectroscopy obtained with the 2.5 duPont Telescope at the Observatorio de Las Campanas (LCO) in Chile; and (e) high-resolution echelle spectroscopy from the LiLiMaRlin project obtained with a variety of spectrographs: HERMES at the 1.2 m Mercator Telescope (MT) at the Observatorio del Roque de los Muchachos (ORM) on the island of La Palma, Spain; ELODIE at the 1.93 m Observatoire de Haute-Provence (OHP) Telescope, France; FEROS at the 2.2 m Telescope of the Observatorio de La Silla in Chile; CAFE at the 2.2 m Centro Astronomico Hispano en Andalucia (CAHA) Telescope, Almeria, Spain; FIES at the 2.5 Nordic Optical Telescope (NOT) at the Observatorio del Roque de los Muchachos (ORM) on the island of La Palma, Spain; and UVES at the 8.2 m Kueyen Telescope at the Observatorio Paranal in Chile. Some of the MT and NOT data were obtained from the IACOB spectroscopic database (Simon-Diaz et al. 2011b,a, 2015b). This paper 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. This paper has also made use of the Washington Double Star (WDS) catalog (Mason et al. 2001) and the Ski ff (2014) catalog of spectral classifications. The authors would like to thank the personnel of the WHT, CAHA, LT, LCO, MT, La Silla, and NOT observatories for their support and hospitality throughout the years. We dedicate this paper to our deceased colleagues, Virpi S. Niemela and Nolan R. Walborn, who they surely would have enjoyed having access to data like the ones presented here., With funding from the Spanish government through the Severo Ochoa Centre of Excellence accreditation SEV-2017-0709.

CONTEXT. Many O-type stars have nearby companions whose presence hamper their characterization through spectroscopy. AIMS. We want to obtain spatially resolved spectroscopy of close O-type visual binaries to derive their spectral types. METHODS. We use the Space Telescope Imaging Spectrograph (STIS) of the Hubble Space Telescope (HST) to obtain long-slit blue-violet spectroscopy of eight Galactic O-type stars with nearby visual companions and use spatial-profile fitting to extract the separate spectra. We also use the ground-based Galactic O-Star Spectroscopic Survey (GOSSS) to study more distant visual components. RESULTS. We spatially resolve seven of the eight systems, present spectra for their components, and obtain their spectral types. Those seven multiple systems are iota Ori Aa,Ab,B, 15 Mon Aa,Ab,C, tau CMa Aa,Ab,B,C,D,E, HD 206 267 Aa,Ab,C,D, HD 193 443 A,B, HD 16 429 Aa,Ab, and IU Aur A,B. This is the first time that spatially resolved spectroscopy of the close visual binaries of those systems is obtained. We establish the applicability of the technique as a function of the separation and magnitude difference of the binary., Comment: Accepted for publication in A&A

Context. Cygnus OB2 provides a unique insight into the high-mass stellar content in one of the largest groups of young massive stars in our Galaxy. Although several studies of its massive population have been carried out over the last decades, an extensive spectroscopic study of the whole known O-star population in the association is still lacking. Aims. We aim to carry out a spectroscopic characterization of all the currently known O stars in Cygnus OB2, determining the distribution of rotational velocities and accurate stellar parameters to obtain an improved view of the evolutionary status of the region. Methods. Based on existing and new optical spectroscopy, we performed a detailed quantitative spectroscopic analysis of all the known O-type stars identified in the association. For this purpose, we used the user-friendly iacob-broad and iacob-gbat automatized tools, FASTWIND stellar models, and astrometry provided by the Gaia second data release. Results. We created the most complete spectroscopic census of O stars carried out so far in Cygnus OB2 using already existing and new spectroscopy. We present the spectra for 78 O-type stars, from which we identify new binary systems, obtain the distribution of rotational velocities, and determine the main stellar parameters for all the stars in the region that have not been detected as double-line spectroscopic binaries. We also derive radii, luminosities, and masses for those stars with reliable Gaia astrometry, in addition to creating the Hertzsprung-Russell Diagram to interpret the evolutionary status of the association. Finally, we inspect the dynamical state of the population and identify runaway candidates. Conclusions. Our spectroscopic analysis of the O-star population in Cygnus OB2 has led to the discovery of two new binary systems and the determination of the main stellar parameters, including rotational velocities, luminosities, masses, and radii for all identified stars. This work has shown the improvement reached when using accurate spectroscopic parameters and astrometry for the interpretation of the evolutionary status of a population, revealing, in the case of Cygnus OB2, at least two star-forming bursts at ∼3 and ∼5 Myr. We find an apparent deficit of very fast rotators in the distribution of rotational velocities. The inspection of the dynamical distribution of the sample has allowed us to identify nine O stars with peculiar proper motions and discuss a possible dynamical ejection scenario or past supernova explosions in the region. © ESO 2020, We acknowledge financial support from the Spanish Ministry of Economy and Competitiveness (MINECO) under the grants AYA2012-39364-C02-01, AYA 2015-68012-C2-01, Severo Ochoa SEV-2015-0548 and PGC2018-093741-B-C21/C22 (MICIU/AEI/FEDER, UE)

Binary systems composed of a red supergiant and a B star are useful probes of stellar evolution. We have searched the literature to create a catalog of 108 Galactic systems of such type, which is presented here.

Massive binaries (MBs) play a crucial role in the Universe. Knowing the distributions of their orbital parameters (OPs) is important for a wide range of topics, from stellar feedback to binary evolution channels, from the distribution of supernova types to gravitational wave progenitors, yet, no direct measurements exist outside the Milky Way. The Tarantula Massive Binary Monitoring was designed to help fill this gap by obtaining multi-epoch radial velocity monitoring of 102 MBs in the 30 Dor. In this paper, we analyse 32 VLT/FLAMES observations of 93 O- and 7 B-type binaries. We performed a Fourier analysis and obtained orbital solutions for 82 systems: 51 single- and 31 double-lined spectroscopic binaries. Overall, the OPs and binary fraction are remarkably similar across the 30 Dor region and compared to existing Galactic samples (GSs). This indicates that within these domains environmental effects are of second order in shaping the properties of MBs. A small difference is found in the distribution of orbital periods (OrbPs), which is slightly flatter (in log space) in 30 Dor than in the Galaxy, although this may be compatible within error estimates and differences in the fitting methodology. Also, OrbPs in 30 Dor can be as short as 1.1 d; somewhat shorter than seen in GSs. Equal mass binaries q>0.95 in 30 Dor are all found outside NGC 2070 the very young and massive cluster at 30 Dor's core. One outstanding exception however is the fact that earliest spectral types (O2-O7) tend to have shorter OrbPs than latter (O9.2-O9.7). Our results point to a relative universality of the incidence rate of MBs and their OPs in the metallicity range from solar ($Z_{\odot}$) to about $0.5Z_{\odot}$. This provides the first direct constraints on MB properties in massive star-forming galaxies at the Universes peak of star formation at redshifts z~1 to 2, which are estimated to have $Z~0.5Z_{\odot}$., 36 pages, 19 figures, 10 tables, published in A\&A

This article has an erratum: A&A 639: C1 (2020), Context. Multiplicity in massive stars is key to understanding the chemical and dynamical evolution of galaxies. Among massive stars, those of O type play a crucial role due to their high masses and short lifetimes. Aims: MONOS (Multiplicity Of Northern O-type Spectroscopic systems) is a project designed to collect information and study O-type spectroscopic binaries with ¿ > -20°. In this first paper we describe the sample and provide spectral classifications and additional information for objects with previous spectroscopic and/or eclipsing binary orbits. In future papers we will test the validity of previous solutions and calculate new spectroscopic orbits. Methods: The spectra in this paper have two sources: the Galactic O-Star Spectroscopic Survey (GOSSS), a project that obtains blue-violet R ¿ 2500 spectroscopy of thousands of massive stars, and LiLiMaRlin, a library of libraries of high-resolution spectroscopy of massive stars obtained from four different surveys (CAFÉ-BEANS, OWN, IACOB, and NoMaDS) and additional data from our own observing programs and public archives. We have also used lucky images obtained with AstraLux. Results: We present homogeneous spectral classifications for 92 O-type spectroscopic multiple systems and ten optical companions, many of them original. We discuss the visual multiplicity of each system with the support of AstraLux images and additional sources. For eleven O-type objects and for six B-type objects we present their first GOSSS spectral classifications. For two known eclipsing binaries we detect double absorption lines (SB2) or a single moving line (SB1) for the first time, to which we add a third system reported by us recently. For two previous SB1 systems we detect their SB2 nature for the first time and give their first separate spectral classifications, something we have also done for a third object just recently identified as a SB2. We also detect nine new astrometric companions and provide updated information on several others. We emphasize the results for two stars: for ¿ Ori AaAbB we provide spectral classifications for the three components with a single observation for the first time thanks to a lucky spectroscopy observation obtained close to the Aa,Ab periastron and for ¿1 Ori CaCb we add it to the class of Galactic Of?p stars, raising the number of its members to six. Our sample of O-type spectroscopic binaries contains more triple- or higher-order systems than double systems.© ESO 2019., Several authors acknowledge support from the Spanish Government Ministerio de Ciencia, Innovacion y Universidades through different grants: AYA2016-75931-C2-1/2-P (J.M.A., E.T.P., A.S., and E.J.A.), AYA2015-68 012-C2-1/2-P (E.T.P., I.N., S.S.-D., J.L., and A.M.), AYA2016-79 425-C3-2-P (J.A.C.), and SEV2015-0548 (S.S.-D.). R.H.B. acknowledges support from the ESAC Faculty Council Visitor Program. R.H.B. and J.I.A. were also supported by the Direccion de Investigacion y Desarrollo de la Universidad de La Serena through projects PR18 143 and PR16 142, respectively. S.S.-D. acknowledges support from the Gobierno de Canarias grant ProID2017010115.

This paper describes the characterization of the GALANTE photometric system, a seven intermediate-and narrow-band filter system with a wavelength coverage from 3000 Å to 9000 Å. We describe the photometric system presenting the full sensitivity curve as a product of the filter sensitivity, CCD, telescope mirror, and atmospheric transmission curves, as well as some first-and second-order moments of this sensitivity function. The GALANTE photometric system is composed of four filters from the J-PLUS photometric system, a twelve broad-to-narrow filter system, and three exclusive filters, specifically designed to measure the physical parameters of stars such as effective temperature T, log (g), metallicity, colour excess E(4405-5495), and extinction type R. Two libraries, the Next Generation Spectral Library (NGSL) and the one presented in Maíz Apellániz & Weiler (2018), have been used to determine the transformation equations between the Sloan Digital Sky Survey (SDSS) ugriz photometry and the GALANTE photometric system. We will use this transformation to calibrate the zero-points of GALANTE images. To this end, a preliminary photometric calibration of GALANTE has been made based on two different griz libraries (SDSS DR12 and ATLAS All-Sky Stellar Reference Catalog, hereinafter RefCat2). A comparison between both zero-points is performed leading us to the choice of RefCat2 as the base catalogue for this calibration, and applied to a field in the Cyg OB2 association.© 2019 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society., AL-G, EJA, and JMA acknowledge support from the Spanish Government Ministerio de Ciencia, Innovacion y Universidades though grants AYA2013-40 611-P and AYA2016-75 931-C2-1/2-P. A.L.-G and EJA also acknowledge support from the State Agency for Research of the Spanish MCIU through the 'Center of Excellence Severo Ochoa' award for the Instituto de Astrofisica de Andalucia (SEV-2017-0709).

We present the photometric calibration of the 12 optical passbands observed by the Javalambre Photometric Local Universe Survey (J-PLUS). Methods. The proposed calibration method has four steps: (i) definition of a high-quality set of calibration stars using Gaia information and available 3D dust maps; (ii) anchoring of the J-PLUS gri passbands to the Pan-STARRS photometric solution, accounting for the variation in the calibration with the position of the sources on the CCD; (iii) homogenization of the photometry in the other nine J-PLUS filters using the dust de-reddened instrumental stellar locus in (chi - r) versus (g-i) colours, where chi is the filter to calibrate. The zero point variation along the CCD in these filters was estimated with the distance to the stellar locus. Finally, (iv) the absolute colour calibration was obtained with the white dwarf locus. We performed a joint Bayesian modelling of 11 J-PLUS colour-colour diagrams using the theoretical white dwarf locus as reference. This provides the needed o ffsets to transform instrumental magnitudes to calibrated magnitudes outside the atmosphere. Results. The uncertainty of the J-PLUS photometric calibration, estimated from duplicated objects observed in adjacent pointings and accounting for the absolute colour and flux calibration errors, are similar to 19 mmag in u, J0378, and J0395; similar to 11 mmag in J0410 and J0430; and similar to 8 mmag in g, J0515, r, J0660, i, J0861, and z. Conclusions. We present an optimized calibration method for the large-area multi-filter J-PLUS project, reaching 1-2% accuracy within an area of 1022 square degrees without the need for long observing calibration campaigns or constant atmospheric monitoring. The proposed method will be adapted for the photometric calibration of J-PAS, that will observe several thousand square degrees with 56 narrow optical filters.© ESO 2019, We dedicate this paper to the memory of our six IAC colleagues and friends who met with a fatal accident in Piedra de los Cochinos, Tenerife, in February 2007, with special thanks to Maurizio Panniello, whose teachings of python were so important for this paper. We thank the anonymous referee for useful comments and suggestions. Based on observations made with the JAST/T80 telescope at the Observatorio Astrofisico de Javalambre (OAJ), in Teruel, owned, managed, and operated by the Centro de Estudios de Fisica del Cosmos de Aragon. We acknowledge the OAJ Data Processing and Archiving Unit (UPAD) for reducing and calibrating the OAJ data used in this work. Funding for the J-PLUS Project has been provided by the Governments of Spain and Aragon through the Fondo de Inversiones de Teruel; the Aragon Government through the Reseach Groups E96, E103, and E16_17R; the Spanish Ministry of Science, Innovation and Universities (MCIU/AEI/FEDER, UE) with grants PGC2018-097585-B-C21 and PGC2018-097585-B-C22; the Spanish Ministry of Economy and Competitiveness (MINECO) under AYA2015-66211-C2-1-P, AYA2015-66211-C2-2, AYA2012-30789, and ICTS-2009-14; and European FEDER funding (FCDD10-4E-867, FCDD13-4E-2685). This work was also supported by the MINECO through grant ESP2016-80079-C2-1-R and RTI2018-095076-B-C21 (MINECO/FEDER, UE) and MDM-2014-0369 of ICCUB (Unidad de Excelencia "Maria de Maeztu"). F. J. E. acknowledges financial support from the Tec2Space-CM project (P2018/NMT-4291). Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the U.S. Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council for England. The SDSS Web Site is www.sdss.org.The SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions. The Participating Institutions are the American Museum of Natural History, Astrophysical Institute Potsdam, University of Basel, University of Cambridge, CaseWestern Reserve University, University of Chicago, Drexel University, Fermilab, the Institute for Advanced Study, the Japan Participation Group, Johns Hopkins University, the Joint Institute for Nuclear Astrophysics, the Kavli Institute for Particle Astrophysics and Cosmology, the Korean Scientist Group, the Chinese Academy of Sciences (LAMOST), Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New Mexico State University, the Ohio State University, University of Pittsburgh, University of Portsmouth, Princeton University, the United States Naval Observatory, and the University of Washington. The Pan-STARRS1 Surveys (PS1) and the PS1 public science archive have been made possible through contributions by the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project O ffice, the Max-Planck Society and its participating institutes, the Max Planck Institute for Astronomy, Heidelberg, and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hopkins University, Durham University, the University of Edinburgh, the Queen's University Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incorporated, the National Central University of Taiwan, the Space Telescope Science Institute, the National Aeronautics and Space Administration under Grant No. NNX08AR22G issued through the Planetary Science Division of the NASA Science Mission Directorate, the National Science Foundation Grant No. AST-1238877, the University of Maryland, Eotvos Lorand University (ELTE), the Los Alamos National Laboratory, and the Gordon and Betty Moore Foundation. 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. This research made use of Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration 2013), and Matplotlib, a 2D graphics package used for Python for publication-quality image generation across user interfaces and operating systems (Hunter 2007).

Context. It has recently been proposed that LB-1 is a binary system at 4 kpc consisting of a B-type star of 8 M¿ and a massive stellar black hole (BH) of 70 M¿. This finding challenges our current theories of massive star evolution and formation of BHs at solar metallicity. Aims: Our objective is to derive the effective temperature, surface gravity, and chemical composition of the B-type component in order to determine its nature and evolutionary status and, indirectly, to constrain the mass of the BH. Methods: We use the non-LTE stellar atmosphere code FASTWIND to analyze new and archival high-resolution data. Results: We determine (Teff, log g) values of (14 000 ± 500 K, 3.50 ± 0.15 dex) that, combined with the Gaia parallax, imply a spectroscopic mass, from log g, of 3.2+2.1-1.9 M¿ and an evolutionary mass, assuming single star evolution, of 5.2+0.3-0.6 M¿. We determine an upper limit of 8 km s-1 for the projected rotational velocity and derive the surface abundances; we find the star to have a silicon abundance below solar, and to be significantly enhanced in nitrogen and iron and depleted in carbon and magnesium. Complementary evidence derived from a photometric extinction analysis and Gaia yields similar results for Teff and log g and a consistent distance around 2 kpc. Conclusions: We propose that the B-type star is a slightly evolved main sequence star of 3-5 M¿ with surface abundances reminiscent of diffusion in late B/A chemically peculiar stars with low rotational velocities. There is also evidence for CN-processed material in its atmosphere. These conclusions rely critically on the distance inferred from the Gaia parallax. The goodness of fit of the Gaia astrometry also favors a high-inclination orbit. If the orbit is edge-on and the B-type star has a mass of 3-5 M¿, the mass of the dark companion would be 4-5 M¿, which would be easier to explain with our current stellar evolutionary models., With funding from the Spanish government through the "María de Maeztu Unit of Excellence" accreditation (MDM-2017-0737)

The wrong versions of Figs. 1 and 2 were published in the paper. The points corresponding to Table 6 in the top panels were missing and a single “dummy” unlabeled point was shown instead. The correct versions are given here. Erratum for: A&A, 616, A149 (2018), https://doi.org/10.1051/0004-6361/201832787 With funding from the Spanish government through the "María de Maeztu Unit of Excellence" accreditation (MDM-2017-0737) Peer review

Context. The first Gaia Data Release (DR1) significantly improved the previously available proper motions for the majority of the Tycho-2 stars. Aims: We wish to detect runaway stars using Gaia DR1 proper motions and compare our results with previous searches. Methods: Runaway O stars and BA supergiants were detected using a 2D proper motion method. The sample was selected using Simbad, spectra from our GOSSS project, literature spectral types, and photometry processed using the code CHORIZOS. Results: We detect 76 runaway stars, 17 (possibly 19) of them with no prior identification as such, with an estimated detection rate of approximately one half of the real runaway fraction. An age effect appears to be present, with objects of spectral subtype B1 and later having traveled for longer distances than runaways of earlier subtypes. We also tentatively propose that the fraction of runaways is lower among BA supergiants that among O stars, but further studies using future Gaia data releases are needed to confirm this. The frequency of fast rotators is high among runaway O stars, which indicates that a significant fraction of them (and possibly the majority) is produced in supernova explosions. The spectral types in Table 2 will be added to the spectral types in the three GOSSS survey papers and in Maíz Apellániz et al. (2018 and the resulting table will be available from the GOSC web site http://gosc.cab.inta-csic.es.© ESO 2018, J.M.A., A.S., and E.T.P. acknowledge support from the Spanish Government Ministerio de Economia, Industria y Competitividad (MINECO/FEDER) through grant No. AYA2016-75931-C2-2-P. M.P.G. acknowledges support from the ESAC Trainee program. R.H.B. acknowledges support from the ESAC Faculty Council Visitor Program. S.S.-D., I.N., and E.T.P. acknowledge support from the Spanish Government Ministerio de Economia, Industria y Competitividad (MINECO/FEDER) through grant No. AYA2015-68012-C2-1/2-P.

The 30 Doradus (30 Dor) nebula in the Large Magellanic Cloud (LMC) is the brightest HII region in the Local Group and a prototype starburst similar to those found in high redshift galaxies. It is thus a stepping stone to understand the complex formation processes of stars in starburst regions across the Universe. Here, we have studied the formation history of massive stars in 30 Dor using masses and ages derived for 452 mainly OB stars from the spectroscopic VLT-FLAMES Tarantula Survey (VFTS). We find that stars of all ages and masses are scattered throughout 30 Dor. This is remarkable because it implies that massive stars either moved large distances or formed independently over the whole field of view in relative isolation. We find that both channels contribute to the 30 Dor massive star population. Massive star formation rapidly accelerated about 8 Myr ago, first forming stars in the field before giving birth to the stellar populations in NGC 2060 and NGC 2070. The R136 star cluster in NGC 2070 formed last and, since then, about 1 Myr ago, star formation seems to be diminished with some continuing in the surroundings of R136. Massive stars within a projected distance of 8 pc of R136 are not coeval but show an age range of up to 6 Myr. Our mass distributions are well populated up to $200\,\mathrm{M}_\odot$. The inferred IMF is shallower than a Salpeter-like IMF and appears to be the same across 30 Dor. By comparing our sample of stars to stellar models in the Hertzsprung-Russell diagram, we find evidence for missing physics in the models above $\log L/\mathrm{L}_\odot=6$ that is likely connected to enhanced wind mass loss for stars approaching the Eddington limit. [abridged], 20 pages (incl. appendix), 15 figures, 2 tables; accepted for publication in A&A

The 30 Doradus star-forming region in the Large Magellanic Cloud is a nearby analogue of large star-formation events in the distant Universe. We determine the recent formation history and the initial mass function (IMF) of massive stars in 30 Doradus based on spectroscopic observations of 247 stars more massive than 15 solar masses ($\mathrm{M}_\odot$). The main episode of massive star formation started about $8\,\mathrm{Myr}$ ago and the star-formation rate seems to have declined in the last $1\,\mathrm{Myr}$. The IMF is densely sampled up to $200\,\mathrm{M}_\odot$ and contains $32\pm12\%$ more stars above $30\,\mathrm{M}_\odot$ than predicted by a standard Salpeter IMF. In the mass range $15-200\,\mathrm{M}_\odot$, the IMF power-law exponent is $1.90^{+0.37}_{-0.26}$, shallower than the Salpeter value of 2.35., Authors' version of paper published in Science at http://science.sciencemag.org/content/359/6371/69 ; 15 pages main text (3 figures), 47 pages supplementary materials (10 figures, 3 tables)

Context. Many massive stars have nearby companions whose presence hamper their characterization through spectroscopy. Aims: We want to obtain spatially resolved spectroscopy of close massive visual binaries to derive their spectral types. Methods: We obtained a large number of short long-slit spectroscopic exposures of five close binaries under good seeing conditions. We selected those with the best characteristics, extracted the spectra using multiple-profile fitting, and combined the results to derive spatially separated spectra. Results: We demonstrate the usefulness of Lucky Spectroscopy by presenting the spatially resolved spectra of the components of each system, in two cases with separations of only 0.''3. Those are ¿ Ori Aa+Ab (resolved in the optical for the first time) and ¿ Ori AaAb+B (first time ever resolved). We also spatially resolve 15 Mon AaAb+B, ¿ Ori AaAb+B (both previously resolved with GOSSS, the Galactic O-Star Spectroscopic Survey), and ¿ Ori AaAb+B, a system with two spectroscopic B+B binaries and a fifth visual component. The systems have in common that they are composed of an inner pair of slow rotators orbited by one or more fast rotators, a characteristic that could have consequences for the theories of massive star formation.© 2018 EDP Sciences. All rights reserved., We acknowledge support from the Spanish Government Ministerio de Economía y Competitividad (MINECO) through grants AYA2016-75 931-C2-1/2-P (J.M.A., A.S., E.T.P., and E.J.A), AYA2015-68 012-C2-1/2-P (S.S.-D. and E.T.P.) and AYA2016-79 425-C3-2-P (J.A.C.)

In this paper, we present high-resolution HST/COS observations of the extreme magnetic O star NGC 1624-2. These represent the first ultraviolet spectra of this archetypal object. We examine the variability of its wind-sensitive resonance lines, comparing it to that of other known magnetic O stars. In particular, the observed variations in the profiles of the CIV and SiIV doublets between low state and high state are the largest observed in any magnetic O-type star, consistent with the expected properties of NGC 1624-2's magnetosphere. We also observe a redshifted absorption component in the low state, a feature not seen in most stars. We present preliminary modelling efforts based on the Analytic Dynamical Magnetosphere (ADM) formalism, demonstrating the necessity of using non-spherically symmetric models to determine wind/magnetospheric properties of magnetic O stars., Comment: 13 pages, 6 figures, accepted for publication by MNRAS (23 November 2018)

Context. On September 2016 the first data from Gaia were released (DR1). The first release included photometry for over 10^9 sources in the very broad G system. Aims. To test the correspondence between G magnitudes in DR1 and the synthetic equivalents derived using spectral energy distributions from observed and model spectrophotometry. To correct the G passband curve and to measure the zero point in the Vega system. Methods. I have computed the synthetic G and Tycho-2 BV photometry for a sample of stars using the Next Generation Spectral Library (NGSL) and the Hubble Space Telescope (HST) CALSPEC spectroscopic standards. Results. I have found that the nominal G passband curve is too blue for the DR1 photometry, as shown by the presence of a color term in the comparison between observed and synthetic magnitudes. A correction to the passband applying a power law in lambda with an exponent of 0.783 eliminates the color term. The corrected passband has a Vega zero point of 0.070$\pm$0.004 magnitudes., Four-page letter, accepted for publication in A&A

The Tarantula region in the Large Magellanic Cloud contains the richest population of spatially resolved massive O-type stars known so far. This unmatched sample offers an opportunity to test models describing their main-sequence evolution and mass-loss properties. Using ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to determine stellar, photospheric and wind properties of 72 presumably single O-type giants, bright giants and supergiants and to confront them with predictions of stellar evolution and of line-driven mass-loss theories. We apply an automated method for quantitative spectroscopic analysis of O stars combining the non-LTE stellar atmosphere model {{\sc fastwind}} with the genetic fitting algorithm {{\sc pikaia}} to determine the following stellar properties: effective temperature, surface gravity, mass-loss rate, helium abundance, and projected rotational velocity. We present empirical effective temperature versus spectral subtype calibrations at LMC-metallicity for giants and supergiants. In the spectroscopic and classical Hertzsprung-Russell diagrams, our sample O stars are found to occupy the region predicted to be the core hydrogen-burning phase by Brott et al. (2011) and K\"{o}hler et al. (2015). Except for five stars, the helium abundance of our sample stars is in agreement with the initial LMC composition. The aforementioned five stars present moderate projected rotational velocities (i.e., $v_{\mathrm{e}}\,\sin\,i\, Comment: 39 pages, 12 figures in the main text, 12 figures in the appendix. Accepted for publication in A&A. Tables C1-C5 will be available at the CDS. Appendix E has been truncated due to size limitations, the full version will be available on A&A

We have used (i) Hubble Space Telescope Advanced Camera for Surveys imaging and Space Telescope Imaging Spectrograph spectroscopy, (ii) ground-based Precision Integrated-Optics Near-infrared Imaging ExpeRiment/Very Large Telescope long-baseline interferometry, and (iii) ground-based spectroscopy from different instruments to study the orbit of the extreme multiple system HD 93 129 Aa,Ab, which is composed of (at least) two very massive stars in a long-period orbit with e > 0.92, which will pass through periastron in 2017/2018. In several ways, the system is an ¿ Car precursor. Around the time of periastron passage, the two very strong winds will collide and generate an outburst of non-thermal hard X-ray emission without precedent in an O+O binary since astronomers have been able to observe above Earth's atmosphere. A coordinated multiwavelength monitoring in the next two years will enable a breakthrough understanding of the wind interactions in such extreme close encounters. Furthermore, we have found evidence that HD 93 129 Aa may be a binary system itself. In that case, we could witness a three-body interaction which may yield a runaway star or a stellar collision close to or shortly after the periastron passage. Either of those outcomes would be unprecedented, as they are predicted to be low-frequency events in the Milky Way., JMA acknowledges support from the Spanish Government Ministerio de Economía y Competitividad (MINECO) through grant AYA2013-40 611-P and from HST GO programs 10205, 10602, and 10898. RHB acknowledges support from FONDECYT Projects1 140 076 and11 121 550. The HST data were obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555.

The IACOB and OWN surveys are two ambitious complementary observational projects which have made available a large multi-epoch spectroscopic database of optical high resolution spectra of Galactic massive O-type stars. As a first step in the study of the full sample of (more than 350) O stars surveyed by the IACOB/OWN projects, we have performed the quantitative spectroscopic analysis of a subsample of 128 stars included in the modern grid of O-type standards for spectral classification. We use semi-automatized tools to determine the set of spectroscopic parameters that can be obtained from the optical spectrum of O-type stars. We also benefit from the multi-epoch character of the surveys to perform a spectroscopic variability study of the sample, accounting for spectroscopic binarity and variability of the main wind diagnostic lines. We provide a general overview of the stellar and wind parameters of this reference sample, and updated recipes for the SpT\,--\,Teff/log g calibrations for Galactic O-type stars. We evaluate our semi-automatized analysis strategy with $\sim$40 stars from the literature, and find a good agreement. The agreement between the synthetic spectra associated with fastwind best fitting models and the observed spectra is good for most targets, but 46 stars present a particular behavior of the wind diagnostic lines that cannot be reproduced by our grid of spherically symmetric unclumped models. These are potential targets of interest for more detailed investigations of clumpy winds and/or the existence of additional circumstellar components. Last, our variability study has led to the detection of signatures of spectroscopic binarity in 27\% of the stars and small amplitude radial velocity variations in the photospheric lines of another 30\%. Additionally, 31\% of the investigated stars show variability in the wind diagnostic lines., Comment: 20 pages, 18 figures, accepted for publication in Astronomy & Astrophysics

Aims: We present the first extensive spectroscopic study of the global population in star clusters Trumpler 16, Trumpler 14, and Collinder 232 in the Carina nebula, using data from the Gaia-ESO Survey, down to solar-mass stars. Methods: In addition to the standard homogeneous survey data reduction, a special processing was applied here because of the bright nebulosity surrounding Carina stars. Results: We find about 400 good candidate members ranging from OB types down to slightly subsolar masses. About 100 heavily reddened early-type Carina members found here were previously unrecognized or poorly classified, including two candidate O stars and several candidate Herbig Ae/Be stars. Their large brightness makes them useful tracers of the obscured Carina population. The spectroscopically derived temperatures for nearly 300 low-mass members enables the inference of individual extinction values and the study of the relative placement of stars along the line of sight. Conclusions: We find a complex spatial structure with definite clustering of low-mass members around the most massive stars and spatially variable extinction. By combining the new data with existing X-ray data, we obtain a more complete picture of the three-dimensional spatial structure of the Carina clusters and of their connection to bright and dark nebulosity and UV sources. The identification of tens of background giants also enables us to determine the total optical depth of the Carina nebula along many sightlines. We are also able to put constraints on the star formation history of the region with Trumpler 14 stars found to be systematically younger than stars in other subclusters. We find a large percentage of fast-rotating stars among Carina solar-mass members, which provide new constraints on the rotational evolution of pre-main-sequence stars in this mass range.© 2017 ESO., This work was partly supported by the European Union FP7 program through ERC grant number 320360 and by the Leverhulme Trust through grant RPG-2012-541.

The VLT-FLAMES Tarantula Survey has observed hundreds of O-type stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). We study the properties of 105 apparently single O-type dwarfs. To determine stellar and wind parameters, we used the IACOB-GBAT package, an automatic procedure based on a large grid of atmospheric models calculated with the FASTWIND code. In addition to classical techniques, we applied the Bayesian BONNSAI tool to estimate evolutionary masses. We provide a new calibration of effective temperature vs. spectral type for O-type dwarfs in the LMC, based on our homogeneous analysis of the largest sample of such objects to date and including all spectral subtypes. Good agreement with previous results is found, although the sampling at the earliest subtypes could be improved. Rotation rates and helium abundances are studied in an evolutionary context. We find that most of the rapid rotators (vsini higher than 300 km/s ) in our sample have masses below 25 MSun and intermediate rotation-corrected gravities (log gc between 3.9 and 4.1). Such rapid rotators are scarce at higher gravities (i.e. younger ages) and absent at lower gravities (larger ages). This is not expected from theoretical evolutionary models, and does not appear to be due to a selection bias in our sample. We compare the estimated evolutionary and spectroscopic masses, finding a trend that the former is higher for masses below 20 MSun. This can be explained as a consequence of limiting our sample to the O-type stars, and we see no compelling evidence for a systematic mass discrepancy. For most of the stars in the sample we were unable to estimate the wind-strength parameter (hence mass-loss rates) reliably, particularly for objects with luminosity lower than logL/LSun about 5.1. Ultraviolet spectroscopy is needed to undertake a detailed investigation of the wind properties of these dwarfs., Comment: 23 pages, 16 figures. Accepted for publication in A&A

[ABRIDGED] Context. O stars are excellent tracers of the intervening ISM because of their high luminosity, blue intrinsic SED, and relatively featureless spectra. We are currently conducting GOSSS, which is generating a large sample of O stars with accurate spectral types within several kpc of the Sun. Aims. To obtain a global picture of the properties of dust extinction in the solar neighborhood based on optical-NIR photometry of O stars with accurate spectral types. Methods. We have processed a photometric set with the CHORIZOS code to measure the amount and type of extinction towards 562 O-type stellar systems. We have tested three different families of extinction laws and analyzed our results with the help of additional archival data. Results. The Ma\'iz Apell\'aniz et al. (2014) family of extinction laws provides a better description of Galactic dust that either the Cardelli et al. (1989) or Fitzpatrick (1999) families, so it should be preferentially used. In many cases O stars and late-type stars experience similar amounts of extinction at similar distances but some O stars are located close to the molecular clouds left over from their births and have larger extinctions than the average for nearby late-type populations. In qualitative terms, O stars experience a more diverse extinction than late-type stars, as some are affected by the small-grain-size, low-R_5495 effect of molecular clouds and others by the large-grain-size, high-R_5495 effect of H II regions. Late-type stars experience a narrower range of grain sizes or R_5495, as their extinction is predominantly caused by the average, diffuse ISM. We propose that the reason for the existence of large-grain-size, high-R_5495 regions in the ISM in the form of H II regions and hot-gas bubbles is the selective destruction of small dust grains by EUV photons and possibly by thermal sputtering by atoms or ions., Comment: 33 pages, accepted for publication in A&A. This second version includes language editing, other minor changes to the text, and the solution to formatting issues in the main table