Search

Context. Current implementations of mass loss for hot, massive stars in stellar evolution models include a sharp increase in mass loss when blue supergiants become cooler than Teff 20-22kK. This drastic mass-loss jump has been motivated by the potential presence of a so-called bistability ionisation effect, which may occur for line-driven winds in this temperature region due to recombination of important line-driving ions. Aims. We perform quantitative spectroscopy using UV (ULLYSES program) and optical (XShootU collaboration) data for 17 OB-supergiant stars in the LMC (covering the range Teff 14-32kK), deriving absolute constraints on global stellar, wind, and clumping parameters. We examine whether there are any empirical signs of a mass-loss jump in the investigated region, and we study the clumped nature of the wind. Methods. We use a combination of the model atmosphere code fastwind and the genetic algorithm code Kiwi-GA to fit synthetic spectra of a multitude of diagnostic spectral lines in the optical and UV. Results. We find no signs of any upward mass loss jump anywhere in the examined region. Standard theoretical comparison models, which include a strong bistability jump thus severely over predict the empirical mass-loss rates on the cool side of the predicted jump. Additionally, we find that on average about 40% of the total wind mass seems to reside in the diluted medium in between dense clumps. Conclusions. Our derived mass-loss rates suggest that for applications like stellar evolution one should not include a drastic bistability jump in mass loss for stars in the temperature and luminosity region investigated here. The derived high values of interclump density further suggest that the common assumption of an effectively void interclump medium (applied in the vast majority of spectroscopic studies of hot star winds) is not generally valid in this parameter regime.

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

Surveys in the Milky Way and Large Magellanic Cloud revealed that the majority of massive stars will interact with companions during their lives. However, knowledge of the binary properties of massive stars at low metallicity, which approaches the conditions of the Early Universe, remains sparse. We present the Binarity at LOw Metallicity (BLOeM) campaign - an ESO large programme designed to obtain 25 epochs of spectroscopy for 929 massive stars in the SMC - the lowest metallicity conditions in which multiplicity is probed to date (Z = 0.2 Zsun). BLOeM will provide (i) the binary fraction, (ii) the orbital configurations of systems with periods P < 3 yr, (iii) dormant OB+BH binaries, and (iv) a legacy database of physical parameters of massive stars at low metallicity. The stars are observed with the LR02 setup of the giraffe instrument of the Very Large Telescope (3960-4570A, resolving power R=6200; typical signal-to-noise ratio S/N=70-100). This paper utilises the first 9 epochs obtained over a three-month time. We describe the survey and data reduction, perform a spectral classification of the stacked spectra, and construct a Hertzsprung-Russell diagram of the sample via spectral-type and photometric calibrations. The sample covers spectral types from O4 to F5, spanning the effective temperature and luminosity ranges 6.5

CONTEXT: The spectral analysis of hot, massive stars is a fundamental astrophysical method to obtain their intrinsic properties and their feedback. Quantitative spectroscopy for hot, massive stars requires detailed numerical modeling of the atmosphere and an iterative treatment to obtain the best solution within a given framework. AIMS: We present an overview of different techniques for the quantitative spectroscopy of hot stars employed within the X-Shooting ULLYSES collaboration, from grid-based approaches to tailored fits. By performing a blind test, we gain an overview about the similarities and differences of the resulting parameters. Our study aims to provide an overview of the parameter spread caused by different approaches. METHODS: For three different stars from the sample (SMC O5 star AzV 377, LMC O7 star Sk -69 50, and LMC O9 star Sk -66 171), we employ different atmosphere codes (CMFGEN, Fastwind, PoWR) and strategies to determine their best-fitting model. For our analyses, UV and optical spectra are used to derive the properties with some methods relying purely on optical data for comparison. To determine the overall spectral energy distribution, we further employ additional photometry from the literature. RESULTS: Effective temperatures for each of three sample stars agree within 3 kK while the differences in log g can be up to 0.2 dex. Luminosity differences of up to 0.1 dex result from different reddening assumptions, which seem to be larger for the methods employing a genetic algorithm. All sample stars are nitrogen-enriched. CONCLUSIONS: We find a reasonable agreement between the different methods. Tailored fitting tends to be able to minimize discrepancies obtained with more course or automatized treatments. UV spectral data is essential for the determination of realistic wind parameters. For one target (Sk -69 50), we find clear indications of an evolved status., Comment: 19+15 pages, 21+4 figures, accepted version (A&A 689, A30) including language editing, condensed abstract

On the route towards merging neutron stars and stripped-envelope supernovae, binary population synthesis predicts a large number of post-interaction systems with massive stars that have stripped off their outer layers. Yet, observations of such stars in the intermediate-mass regime below the Wolf-Rayet masses are rare. Using X-Shooting ULLYSES (XShootU) data, we discovered three partially stripped star + Be/Oe binaries in the Magellanic Clouds. We analyzed the UV and optical spectra using the PoWR model atmosphere code by superimposing model spectra corresponding to each component. The estimated current masses of the partially stripped stars fall within the intermediate mass range of 4-8 $M_{\odot}$. These objects are overluminous for their stellar masses, matching core He-burning luminosities. Their Be/Oe secondaries have much higher masses than their stripped primaries (mass ratio > 2). All three partially stripped stars show significant nitrogen enrichment and carbon and oxygen depletion on their surfaces. Additionally, one of our sample stars exhibits significant helium enrichment. Our study provides the first comprehensive determination of the wind parameters of partially stripped stars in the intermediate mass range. The wind mass-loss rates of these stars are found to be on the order of $10^{-7} M_\odot$ yr$^{-1}$, which is over ten times higher than that of OB stars of the same luminosity. Current evolutionary models characterizing this phase typically employ OB or WR mass-loss rates, which underestimate or overestimate stripped stars' mass-loss rates by an order of magnitude. Binary evolution models indicate that the observed primaries had initial masses of 12-17 $M_{\odot}$, making them potential candidates for stripped-envelope supernovae that form neutron stars. If they survive the explosion, these systems may become Be X-ray binaries and later double neutron stars., Comment: Accepted for publication in Astronomy & Astrophysics

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 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