Search

The outcome of the formation of massive stars is an important anchor point in their evolution. It provides insight into the physics of the assembly process, and sets the conditions for stellar evolution. We characterize a population of 18 highly reddened O4.5 to B9 stars in the very young massive star-forming region M17. Their properties allow us to identify the empirical location of the ZAMS, and rotation and mass-loss rate of stars there. We performed quantitative spectroscopic modeling of VLT/X-shooter spectra using NLTE atmosphere code Fastwind and fitting approach Kiwi-GA. The observed SEDs were used to determine the line-of-sight extinction. From a comparison of their positions in the HRD with MIST evolutionary tracks, we inferred the stellar masses and ages. We find an age of $0.4_{-0.2}^{+0.6}$ Myr for our sample, however we also identify a strong relation between the age and the mass of the stars. The extinction towards the sources ranges from $A_V = 3.6$ to 10.6. Stars more massive than 10 M$_{\odot}$ have reached the ZAMS. Their projected ZAMS spin rate distribution extends to 0.3 of the critical velocity; their mass-loss rates agree with those of other main-sequence OB stars. Stars with a mass in the range $3 < M/$M$_{\odot} < 7$ are still on the pre-main sequence (PMS). Evolving their $v \sin i$ to the ZAMS yields values up to $\sim 0.6 v_{\rm crit}$. For PMS stars without disks, we find tentative mass-loss rates up to $10^{-8.5}\,$M$_{\odot}$\,yr$^{-1}$. We constrain the empirical location of the ZAMS for massive ($10 < M/$M$_{\odot} < 50$) stars. The ZAMS rotation rates for intermediate-mass stars are twice as high as for massive stars, suggesting that the angular momentum gain processes differ between the two groups. The relation between the age and mass of the stars suggests a lag in the formation of more massive stars relative to lower mass stars., Comment: 37 pages, 31 figures, accepted for publication in A&A

Significant progress has been made toward understanding the formation of massive ($M > 8~$M$_{\odot}$) binaries in close orbits. For example, the detection of a very low velocity dispersion among the massive stars in the young region M17 and the measurement of a positive trend of velocity dispersion with age in Galactic clusters. The velocity dispersion observed in M17 could be explained either by the lack of binaries among the stars in this region or by larger binary separations than typically observed, but with a binary fraction similar to other young Galactic clusters. The latter implies that over time, the binary components migrate toward each other. We aim to determine the origin of the strikingly low velocity dispersion by determining the observed and intrinsic binary fraction of massive stars in M17 through multi-epoch spectroscopy. We performed a multi-epoch spectroscopic survey consisting of three epochs separated by days and months. We determine the radial velocity of each star at each epoch by fitting the stellar absorption profiles. We determine an observed binary fraction of 27% and an intrinsic binary fraction of 87%, consistent with that of other Galactic clusters. We conclude that the low velocity dispersion is due to a large separation among the young massive binaries in M17. Our result is in agreement with a migration scenario in which massive stars are born in binaries or higher order systems at large separation and harden within the first million years of evolution. Such an inward migration may either be driven by interaction with a remnant accretion disk, with other young stellar objects present in the system or by dynamical interactions within the cluster. Our results imply that possibly both dynamical interactions and binary evolution are key processes in the formation of gravitational wave sources., Comment: Accepted for publication in A&A. 18 pages, 8 figures. Appendices A and C will be available in electronic form at the CDS

Using the medium resolution spectrograph X-shooter, spectra of 235 OB and Wolf-Rayet (WR) stars in sub-solar metallicity environments have been secured. [...]This second paper focuses on the optical observations of 232 Magellanic Clouds targets. It describes the uniform reduction of the UVB (300 - 560 nm) and VIS (550 - 1020 nm) XShootU data as well as the preparation of advanced data products [...] . The data reduction of the raw data is based on the ESO CPL X-shooter pipeline. We paid particular attention to the determination of the response curves [...] We implemented slit-loss correction, absolute flux calibration, (semi-)automatic rectification to the continuum, and a correction for telluric lines. The spectra of individual epochs were corrected for the barycentric motion, re-sampled and co-added, and the spectra from the two arms were merged into a single flux calibrated spectrum covering the entire optical range with maximum signal-to-noise ratio. [...] We provide three types of data products: (i) two-dimensional spectra for each UVB and VIS exposure; (ii) one-dimensional UVB and VIS spectra before and after response-correction, as well as after applying various processing, including absolute flux calibration, telluric removal, normalisation and barycentric correction; and (iii) co-added flux-calibrated and rectified spectra over the full optical range, for which all available XShootU exposures were combined. For many of the targets, the final signal-to-noise ratio per resolution element is above 200 in both the UVB and the VIS co-added spectra. The reduced data and advanced scientific data products will be made available to the community upon publication of this paper. [...], Comment: 22 pages, 19 figures; accepted for publication in Astronomy & Astrophysics. Links to online tables and databases will be included upon publication

It is a challenge to study the formation process of massive stars: their formation time is short, they are few, often deeply embedded, and at relatively large distances. Our strategy is to study the outcome of the star formation process and to look for signatures remnant of the formation. We have access to a unique sample of (massive) pre-main-sequence (PMS) stars in the giant HII region M17, showing a photosphere and circumstellar disk. The aim is to determine the variability properties of the hot gaseous disks to understand the physical origin of the emission lines and identify dominant physical processes in these disks. We have obtained multiple-epoch (4-5 epochs) VLT/X-shooter spectra of six young stars in M17 covering about a decade. Using stacked spectra we update the spectral classification and identify circumstellar features. With the temporal variance method (TVS) we determine the extent and amplitude of the spectral line variations. The double-peaked emission lines in the PMS stars with gaseous disks are used to determine peak-to-peak velocities, V/R-ratios and the radial velocity of the systems. We identify many disk features, under which a new detection of CO bandhead and CI emission. In three of the stars we detect spectral variability, mainly in lines originating in the circumstellar disk, in a velocity range up to 320 km/s. In two PMS stars the ratio between the blue and red peaks shows a correlation with the peak-to-peak velocity, possibly explained by a spiral-arm structure in the disk. The PMS stars with variability are at similar positions in the HRD but show significant differences in disk lines and variability. The extent and timescale of the variability differs for each star and per line (sets). We find indications for an accretion flow, slow disk winds and/or disk structures in the hot gaseous inner disk as the cause of the variability in these PMS stars., Comment: 27 pages, 24 figures, accepted for publication in Astronomy and Astrophysics, abstract abbreviated

Observations of individual massive stars, super-luminous supernovae, gamma-ray bursts, and gravitational-wave events involving spectacular black-hole mergers, indicate that the low-metallicity Universe is fundamentally different from our own Galaxy. Many transient phenomena will remain enigmatic until we achieve a firm understanding of the physics and evolution of massive stars at low metallicity (Z). The Hubble Space Telescope has devoted 500 orbits to observe 250 massive stars at low Z in the ultraviolet (UV) with the COS and STIS spectrographs under the ULLYSES program. The complementary ``X-Shooting ULLYSES'' (XShootU) project provides enhanced legacy value with high-quality optical and near-infrared spectra obtained with the wide-wavelength coverage X-shooter spectrograph at ESO's Very Large Telescope. We present an overview of the XShootU project, showing that combining ULLYSES UV and XShootU optical spectra is critical for the uniform determination of stellar parameters such as effective temperature, surface gravity, luminosity, and abundances, as well as wind properties such as mass-loss rates in function of Z. As uncertainties in stellar and wind parameters percolate into many adjacent areas of Astrophysics, the data and modelling of the XShootU project is expected to be a game-changer for our physical understanding of massive stars at low Z. To be able to confidently interpret James Webb Space Telescope (JWST) spectra of the first stellar generations, the individual spectra of low Z stars need to be understood, which is exactly where XShootU can deliver., Comment: Accepted in A&A - 35 Pages, 12 Figures, 4 Tables, 2 Large Tables

Recently much progress has been made in probing the embedded stages of massive star formation, pointing to formation scenarios akin to a scaled up version of low-mass star formation. However, the latest stages of massive star formation have rarely been observed. Using 1st and 2nd overtone CO bandhead emission and near- to mid-infrared photometry we aim to characterize the remnant formation disks around 5 unique pre-main-sequence (PMS) stars with masses $6-12~\rm M_{\odot}$, that have constrained stellar parameters thanks to their detectable photospheres. We seek to understand this emission and the disks it originates from in the context of the evolutionary stage of the studied sources. We use an analytic LTE disk model to fit the CO bandhead and the dust emission, found to originate in different disk regions. For the first time we modeled the 2nd overtone emission. Furthermore, we fit continuum normalized bandheads and show the importance of this in constraining the emission region. We also include $^{13}\rm CO$ in our models as an additional probe of the young nature of the studied objects. We find that the CO emission originates in a narrow region close to the star (<1 AU) and under very similar disk conditions (temperatures and densities) for the different objects. This is consistent with previous modeling of this emission in a diverse range of young stellar objects. We discuss these results in the context of the positions of these PMS stars in the Hertzsprung-Russel diagram and the CO emission's association with early age and high accretion rates in (massive) young stellar objects. We conclude that, considering their mass range and for the fact that their photospheres are detected, the M17 PMS stars are observed in a relatively early formation stage. They are therefore excellent candidates for longer wavelength studies to further constrain the end stages of massive star formation., Comment: 21 pages, 12 figures

The young massive-star-forming region M17 contains optically visible massive pre-main-sequence stars that are surrounded by circumstellar disks. Such disks are expected to disappear when these stars reach the main sequence. The physical and dynamical structure of these remnant disks are poorly constrained, especially the inner regions where accretion, photo-evaporation, and companion formation and migration may be ongoing. We aim to constrain the physical properties of the inner parts of the circumstellar disks of massive young stellar objects B243 (6 Msun) and B331 (12 Msun), two systems for which the central star has been detected and characterized previously despite strong dust extinction. Two-dimensional radiation thermo-chemical modelling with ProDiMo of double-peaked hydrogen lines of the Paschen and Brackett series observed with X-shooter was used to probe the properties of the inner disks. Additionally, the dust structure was studied by fitting the optical and near-infrared spectral energy distribution. B243 features a hot gaseous inner disk with dust at the sublimation radius at 3 AU. The disk appears truncated at roughly 6.5 AU; a cool outer disk of gas and dust may be present, but it cannot be detected with our data. B331 also has a hot gaseous inner disk. A gap separates the inner disk from a colder dusty outer disk starting at up to 100 AU. In both sources the inner disk extends to almost the stellar surface. Chemistry is essential for the ionization of hydrogen in these disks. The lack of a gap between the central objects and these disks suggests that they accrete through boundary-layer accretion. This would exclude the stars having a strong magnetic field. Their structures suggest that both disks are transitional in nature, that is to say they are in the process of being cleared, either through boundary-layer accretion, photo-evaporation, or through companion activity., Comment: 18 pages, 16 figures. Accepted for publication in A&A

Spectroscopic multiplicity surveys of O stars in young clusters and OB associations have revealed that a large portion ($\sim$ 70%) of these massive stars (M$_{i}$ $\gt$ 15 $M_{\odot}$) belong to close and short-period binaries (physical separation d $\lt$few au). Despite the recent and significant progress, the formation mechanisms leading to such close massive multiple systems remain to be elucidated. As a result, young massive close binaries (or higher-order multiple systems) are unique laboratories to figure out the pairing mechanism of high-mass stars. We present the first VLTI/GRAVITY observations of six young O stars in the M17 star-forming region ($\lesssim$ 1 Myr) and two additional foreground stars. From the interferometric model fitting of visibility amplitudes and closure phases, we search for companions and measure their positions and flux ratios. Combining the resulting magnitude difference with atmosphere models and evolutionary tracks, we further constrain the masses of the individual components. All of the six high-mass stars are in multiple systems, leading to a multiplicity fraction (MF) of 100%, yielding a 68% confidence interval of 94-100%. We detect a total number of 9 companions with separations up to 120 au. Including previously identified spectroscopic companions, the companion fraction of the young O stars in our sample reaches 2.3$\pm$0.6. The derived masses span a wide range from 2.5 to 50 $M_{\odot}$, with a great tendency towards high-mass companions. While based on a modest sample, our results clearly indicate that the origin of the high degree of multiplicity is rooted in their star formation mechanism. No clear evidence for one of the competing concepts of massive star formation (core accretion or competitive accretion) could be found. However, our results are compatible with migration as a scenario for the formation of close massive binaries., Comment: 20 pages, 14 figures, 4 tables, 3 appendices

The majority of massive stars ($>8$ $\rm{M_{\odot}}$) in OB associations are found in close binary systems. Nonetheless, the formation mechanism of these close massive binaries is not understood yet. Using literature data, we measured the radial-velocity dispersion ($\sigma_\mathrm{RV}$) as a proxy for the close binary fraction in ten OB associations in the Galaxy and the Large Magellanic Cloud, spanning an age range from 1 to 6 Myrs. We find a positive trend of this dispersion with the cluster's age, which is consistent with binary hardening. Assuming a universal binary fraction of $f_\mathrm{bin}$ = 0.7, we converted the $\sigma_\mathrm{RV}$ behavior to an evolution of the minimum orbital period $P_\mathrm{cutoff}$ from $\sim$9.5 years at 1 Myr to $\sim$1.4 days for the oldest clusters in our sample at $\sim$6 Myr. Our results suggest that binaries are formed at larger separations, and they harden in around 1 to 2 Myrs to produce the period distribution observed in few million year-old OB binaries. Such an inward migration may either be driven by an interaction with a remnant accretion disk or with other young stellar objects present in the system. Our findings constitute the first empirical evidence in favor of migration as a scenario for the formation of massive close binaries., Comment: 7 pages. 5 figures. Accepted for publication in Astronomy & Astrophysics

51 Oph is a Herbig Ae/Be star that exhibits strong near-infrared CO ro-vibrational emission at 2.3 micron, most likely originating in the innermost regions of a circumstellar disc. We aim to obtain the physical and geometrical properties of the system by spatially resolving the circumstellar environment of the inner gaseous disc. We used the second-generation VLTI/GRAVITY to spatially resolve the continuum and the CO overtone emission. We obtained data over 12 baselines with the auxiliary telescopes and derive visibilities, and the differential and closure phases as a function of wavelength. We used a simple LTE ring model of the CO emission to reproduce the spectrum and CO line displacements. Our interferometric data show that the star is marginally resolved at our spatial resolution, with a radius of 10.58+-2.65 Rsun.The K-band continuum emission from the disc is inclined by 63+-1 deg, with a position angle of 116+-1 deg, and 4+-0.8 mas (0.5+-0.1 au) across. The visibilities increase within the CO line emission, indicating that the CO is emitted within the dust-sublimation radius.By modelling the CO bandhead spectrum, we derive that the CO is emitted from a hot (T=1900-2800 K) and dense (NCO=(0.9-9)x10^21 cm^-2) gas. The analysis of the CO line displacement with respect to the continuum allows us to infer that the CO is emitted from a region 0.10+-0.02 au across, well within the dust-sublimation radius. The inclination and position angle of the CO line emitting region is consistent with that of the dusty disc. Our spatially resolved interferometric observations confirm the CO ro-vibrational emission within the dust-free region of the inner disc. Conventional disc models exclude the presence of CO in the dust-depleted regions of Herbig AeBe stars. Ad hoc models of the innermost disc regions, that can compute the properties of the dust-free inner disc, are therefore required., Comment: Accepted in A&A

The process of massive star ($M\geq8~M_\odot$) formation is still poorly understood. Observations of massive young stellar objects (MYSOs) are challenging due to their rarity, short formation timescale, large distances, and high circumstellar extinction. Here, we present the results of a spectroscopic analysis of a population of MYSOs in the Large Magellanic Cloud (LMC). We took advantage of the spectral resolution and wavelength coverage of X-shooter (300-2500 nm), mounted on the European Southern Observatory Very Large Telescope, to detect characteristic spectral features in a dozen MYSO candidates near 30 Doradus, the largest starburst region in the Local Group hosting the most massive stars known. The X-shooter spectra are strongly contaminated by nebular emission. We used a scaling method to subtract the nebular contamination from our objects. We detect H$\alpha,\beta$, [O I] 630.0 nm, Ca II infrared triplet, [Fe II] 1643.5 nm, fluorescent Fe II 1687.8 nm, H$_2$ 2121.8 nm, Br$\gamma$, and CO bandhead emission in the spectra of multiple candidates. This leads to the spectroscopic confirmation of 10 candidates as bona fide MYSOs. We compare our observations with photometric observations from the literature and find all MYSOs to have a strong near-infrared excess. We compute lower limits to the brightness and luminosity of the MYSO candidates, confirming the near-infrared excess and the massive nature of the objects. No clear correlation is seen between the Br$\gamma$ luminosity and metallicity. Combining our sample with other LMC samples results in a combined detection rate of disk features such as fluorescent Fe II and CO bandheads which is consistent with the Galactic rate (40\%). Most of our MYSOs show outflow features., Comment: Accepted to A&A, 26 pages, 19 figures, 4 tables

Context: The identification and characterisation of populations of young massive stars in (giant) HII regions provides important constraints on i) the formation process of massive stars and their early feedback on the environment, and ii) the initial conditions for population synthesis models predicting the evolution of ensembles of stars. Aims: We identify and characterise the stellar populations of the following young giant HII regions: M8, G333.6-0.2, and NGC6357. Methods: We acquired H- and K-band spectra of around 200 stars using The K-band KMOS on the ESO Very Large Telescope. The targets for M8 and NGC6357 were selected from the MYStIX project, which combines X-ray observations with near-infrared and mid-infrared data. For G333.6-0.2, the sample selection is based on the near-infrared colours combined with X-ray data. We introduce an automatic spectral classification method in order to obtain temperatures and luminosities for the observed stars. We analyse the stellar populations using their photometric, astrometric, and spectroscopic properties and compared the position of the stars in the Hertzprung-Russell diagram with stellar evolution models to constrain their ages and mass ranges. Results: We confirm the presence of candidate ionising sources in the three regions and report new ones, including the first spectroscopically identified O stars in G333.6-0.2. In M8 and NGC6357, two populations are identified: i) OB main-sequence stars ($M > 5~\rm{M_{\odot}}$) and ii) pre-main sequence stars ($M\approx0.5-5~\rm{M_{\odot}}$). The ages of the clusters are $\sim$1-3~Myr, $< 3$~Myr, and $\sim$0.5-3~Myr for M8, G333.6-0.2, and NGC6357, respectively. We show that MYStIX selected targets have $>$ 90\% probability of being members of the HII region, whereas a selection based on near infrared (NIR) colours leads to a membership probability of only $\sim$70\%., Comment: Accepted for publication in A&A. The arXiv version includes appendix C, which is an online only figure in A&A

Diffuse interstellar bands (DIBs) are broad absorption features measured in sightlines probing the diffuse interstellar medium. Although large carbon-bearing molecules have been proposed as the carriers producing DIBs, their identity remains unknown. The sight line to the young massive star-forming region M17 shows anomalous extinction in the sense that the total-to-selective extinction parameter differs significantly from the average Galactic value and may reach values $R_{V} > 4$. Given the high $R_V$ values, we investigate whether the DIBs in sight lines towards young OB stars in M17 show a peculiar behaviour. We measure the properties of the most prominent DIBs in M17 and study these as a function of $E(B-V)$ and $R_{V}$. The DIB strengths in M17 concur with the observed relations between DIB equivalent width and reddening $E(B-V)$ in Galactic sight lines. For several DIBs we discover a linear relation between the normalised DIB strength EW/$A_{V}$ and $R_{V}^{-1}$. These trends suggest two groups: (i) a group of ten moderately strong DIBs that show a sensitivity to changes in $R_{V}$ that is modest and proportional to DIB strength, and (ii) a group of four very strong DIBs that react sensitively and to a similar degree to changes in $R_{V}$, but in a way that does not appear to depend on DIB strength. The DIB behaviour as a function of reddening is not peculiar in sight lines to M17. Also, we do not detect anomalous DIB profiles as seen in Her 36. DIBs are stronger, per unit visual extinction, in sight lines characterised by a smaller value of $R_{V}$ (large fraction of small dust particles). New relations between extinction normalised DIB strengths, EW/$A_V$, and $R_V$ support the idea that DIB carriers and interstellar dust are connected. Given the distinct behaviour of two groups of DIBs, different types of carriers do not necessarily relate to the dust grains in a similar way., Comment: Abstract has been shortened. Accepted for publication in A&A. 14 pages, 7 pages of appendix, 28 figures