Search

Luminous Blue Variables (LBVs) are massive stars that are likely to be a transitionary phase between O stars and hydrogen-free classical Wolf–Rayet stars. The variability of these stars has been an area of study for both professional and amateur astronomers for more than a century. In this paper, we present 5 yr of precision photometry of the classical LBV P Cygni taken with the BRITE-Constellation nanosatellites. We have analyzed these data with Fourier analysis to search for periodicities that could elucidate the drivers of variability for these stars. These data show some long-time-scale variability over the course of all six calendar years of observations, but the frequencies needed to reproduce the individual light curves are not consistent from 1 yr to the next. These results likely show that there is no periodic phenomenon present for P Cygni, meaning that the variability is largely stochastic. We interpret the data as being caused by internal gravity waves similar to those seen in other massive stars, with P Cygni exhibiting a larger amplitude and lower characteristic frequency than the main-sequence or blue supergiant stars previously studied. These results show evidence that LBVs may be an extrapolation of the blue supergiants, which have previously been shown to be an extension of main-sequence stars in the context of the stochastic low-frequency photometric variability.

We present a uniform analysis of six examples of embedded wind shock (EWS) O star X-ray sources observed at high resolution with the Chandra grating spectrometers. By modeling both the hot plasma emission and the continuum absorption of the soft X-rays by the cool, partially ionized bulk of the wind we derive the temperature distribution of the shock-heated plasma and the wind mass-loss rate of each star. We find a similar temperature distribution for each star's hot wind plasma, consistent with a power-law differential emission measure, $\frac{d\log EM}{d\log T}$, with a slope a little steeper than -2, up to temperatures of only about $10^7$ K. The wind mass-loss rates, which are derived from the broadband X-ray absorption signatures in the spectra, are consistent with those found from other diagnostics. The most notable conclusion of this study is that wind absorption is a very important effect, especially at longer wavelengths. More than 90 per cent of the X-rays between 18 and 25 Angstrom units produced by shocks in the wind of $\zeta$ Puppis are absorbed, for example. It appears that the empirical trend of X-ray hardness with spectral subtype among O stars is primarily an absorption effect., Comment: MNRAS accepted; 12 pages, 8 figures; this version incorporates the content of an erratum which corrects two figures, but leaves the results and conclusions unchanged

We present the results of a four-month, spectroscopic campaign of the Wolf-Rayet dust-making binary, WR137. We detect only small-amplitude, random variability in the CIII5696 emission line and its integrated quantities (radial velocity, equivalent width, skewness, kurtosis) that can be explained by stochastic clumps in the wind of the WC star. We find no evidence of large-scale, periodic variations often associated with Corotating Interaction Regions that could have explained the observed intrinsic continuum polarization of this star. Our moderately high-resolution and high signal-to-noise average Keck spectrum shows narrow double-peak emission profiles in the Halpha, Hbeta, Hgamma, HeII6678 and HeII5876 lines. These peaks have a stable blue-to-red intensity ratio with a mean of 0.997 and a root-mean-square of 0.004, commensurate with the noise level; no variability is found during the entire observing period. We suggest that these profiles arise in a decretion disk around the O9 companion, which is thus an O9e star. The characteristics of the profiles are compatible with those of other Be/Oe stars. The presence of this disk can explain the constant component of the continuum polarization of this system, for which the angle is perpendicular to the plane of the orbit, implying that the rotation axis of the O9e star is aligned with that of the orbit. It remains to be explained why the disk is so stable within the strong ultraviolet radiation field of the O star. We present a binary evolutionary scenario that is compatible with the current stellar and system parameters., 12 pages, 5 Figures

We report the results of searching pulsar-like candidates from the unidentified objects in the Third Catalog of Hard Fermi-LAT Sources (3FHL). Using a machine-learning-based classification scheme with a nominal accuracy of $\sim \!98{{\, \rm per\, cent}}$, we have selected 27 pulsar-like objects from 200 unidentified 3FHL sources for an identification campaign. Using archival data, X-ray sources are found within the γ-ray error ellipses of 10 3FHL pulsar-like candidates. Within the error circles of the much better constrained X-ray positions, we have also searched for the optical/infrared counterparts and examined their spectral energy distributions. Among our shortlisted candidates, the most secure identification is the association of 3FHL J1823.3–1339 and its X-ray counterpart with the globular cluster Mercer 5. The γ-rays from the source can be contributed by a population of millisecond pulsars residing in the cluster. This makes Mercer 5 as one of the slowly growing hard γ-ray population of globular clusters with emission >10 GeV. Very recently, another candidate picked by our classification scheme, 3FHL J1405.1–6118, has been identified as a new γ-ray binary with an orbital period of 13.7 d. Our X-ray analysis with a short Chandra observation has found a possible periodic signal candidate of ∼1.4 h and a putative extended X-ray tail of ∼20 arcsec long. Spectral energy distribution of its optical/infrared counterpart conforms with a blackbody of Tbb ∼ 40 000 K and Rbb ∼ 12 R⊙ at a distance of 7.7 kpc. This is consistent with its identification as an early O star as found by infrared spectroscopy.

We obtain distances to 383 Galactic Wolf-Rayet (WR) stars from $\textit{Gaia}$ DR2 parallaxes and Bayesian methods, with a prior based on H$\mathrm{\scriptsize{II}}$ regions and dust extinction. Distances agree with those from Bailer-Jones et al. for stars up to 2 kpc from the Sun, though deviate thereafter due to differing priors, leading to modest reductions in luminosities for recent WR spectroscopic results. We calculate visual and K-band absolute magnitudes, accounting for dust extinction contributions and binarity, and identify 187 stars with reliable absolute magnitudes. For WR and O stars within 2 kpc, we find a WR/O ratio of 0.09. The distances are used to generate absolute magnitude calibrations and obtain the $\mathrm{\textit{Gaia}}$ colour magnitude diagram for WR stars. Average $v^{\rm WR}$-band absolute magnitudes for WN stars range from -3.6 mag (WN3-4) to -7.0 mag (WN8-9ha), and -3.1 (WO2-4) to -4.6 mag (WC9), with standard deviations of $\sim$0.6 mag. Using H$\mathrm{\scriptsize{II}}$ region scale heights, we identify 31 WR stars at large (3$\sigma$, |z|$\geq$156 pc) distances from the mid-plane as potential runaways accounting for the Galactic warp, of which only 4 involve WN8-9 stars, contrary to previous claims., Comment: Amended section 3.3, table 4, table 5 and table 6 based on corrected Av to AK conversion. Minor updates to figures, WR/O ratio and WR stars with reliable results. 19 pages, 15 figures. Additional appendices updated (57 pages, 5 figures). To be published in MNRAS

Context. Luminous Blue Variables (LBVs) are thought to be in a transitory phase between O stars on the main-sequence and the Wolf-Rayet stage. Recent studies suggest that they might be formed through binary interaction. Only a few are known in binary systems but their multiplicity fraction is uncertain. Aims. This study aims at deriving the binary fraction among the Galactic (confirmed and candidate) LBV population. We combine multi-epoch spectroscopy and long-baseline interferometry. Methods. We use cross-correlation to measure their radial velocities. We identify spectroscopic binaries through significant RV variability (larger than 35 km/s). We investigate the observational biases to establish the intrinsic binary fraction. We use CANDID to detect interferometric companions, derive their parameters and positions. Results. We derive an observed spectroscopic binary fraction of 26 %. Considering period and mass ratio ranges from Porb=1 to 1000 days, and q = 0.1-1.0, and a representative set of orbital parameter distributions, we find a bias-corrected binary fraction of 62%. From interferometry, we detect 14 companions out of 18 objects, providing a binary fraction of 78% at projected separations between 1 and 120 mas. From the derived primary diameters, and the distances of these objects, we measure for the first time the exact radii of Galactic LBVs to be between 100 and 650 Rsun, making unlikely to have short-period systems. Conclusions. This analysis shows that the binary fraction among the Galactic LBV population is large. If they form through single-star evolution, their orbit must be initially large. If they form through binary channel that implies that either massive stars in short binary systems must undergo a phase of fully non-conservative mass transfer to be able to sufficiently widen the orbit or that LBVs form through merging in initially binary or triple systems., Comment: 14 pages, 9 figures, 4 tables; Accepted for publication in A&A. (Abstract shortened to respect the arXiv limit of 1920 characters)

WR22 = HD 92740 is a bright (V = 6.4 $mag$), intrinsically luminous, double-line WN7h + O9III-V binary exhibiting one sharp 8\% deep eclipse near periastron in its elliptical (e = 0.6) 80-day orbit, when the WR-star passes in front of the O star, with no secondary eclipse. We apply two models (L96, A13) to probe the optical space-based light curves from {\em BRITE-Constellation}, including three separate, complete eclipses, that show increased (o-c) scatter compared to the rest of the observations outside the eclipses, likely due to O-star light encountering WR wind-clumps. L96 is a simple atmospheric-eclipse model, often applied to close WR+O binaries, where the O-star is considered a point-source. A13 considers a finite-disk O-star and allows for atmospheric, photospheric and reflection components to the eclipse, permitting a better characterization of its shape through a more physically realistic description of the structures for both stars in WR22. Nevertheless, A13 is still susceptible to uncertainties in the luminosity of the O-star before unique values for the orbital inclination and WR mass-loss rate can be estimated. We present solutions for the two extremes of the O-star, O9V and O9III. As photometry alone cannot allow us to discriminate between these, we compared our results to the spectral models found in the literature and determined the correct solution to be O9V. Our best-fit A13 Model 1 gives $i = 83.5 \pm 0.4^{\circ}$, with $\dot M_{\rm WR} = (1.86 \pm 0.2) \times 10^{-5} \dot M_{\odot}/yr$. The flux ratio in the red {\em BRITE} band in this model is $F_{\rm O}/F_{\rm WR} = 0.064\pm 0.002$., 15 pages, 11 figures

In large scale maps of the Galactic disc, the Carina Arm stands out as a clear spiral feature, hosting prominent star clusters and associations rich in massive stars. We study the proper motions of 4199 O and early B most likely in the far Carina Arm, at distances mainly in excess of 4 kpc from the Sun, within the sky region, $282^{\circ} < \ell < 294^{\circ}$ and $-3^{\circ} < b < +1^{\circ}$ (Galactic coordinates). The sample is constructed by extending an existing blue-selected catalogue, and cross-matching with Gaia EDR3 astrometry. The observed pattern of proper motions is modulated into a saw-tooth pattern, with full amplitude approaching 1 mas yr$^{-1}$, recurring roughly every 2--3 degrees of longitude (200--300 pc at the median OB-star distance of 5.8 kpc). Kinematic perturbation of underlying circular rotation is most likely present. The data also reveal a moving group containing $>50$ OB stars at $\ell \sim 286^{\circ}$, $b \sim -1^{\circ}.4$ behind the main run of the far arm. An analysis of relative proper motions is performed that yields an incidence of runaway O stars of at least 10\% (potentially $>20$\% when full space motions become available). To map where runaways have run away from, we set up simulations for the region that assume linear trajectories and test for trajectory impact parameter in order to identify likely ejection hot spots. We find the method currently gives good results for times of flight of up to $\sim$4 Myr. It shows convincingly that only NGC 3603 and Westerlund 2 have ejected OB stars in significant numbers. Indeed, both clusters have experienced intense spells of ejection between 0.6--0.9 and 0.5--0.8 Myr ago, respectively., 19 pages, 12 figures, (accepted 5th October for publication in MNRAS, in original form 6th July)

The $2p$-$3s$ lines of Fe XVII in the X-ray spectrum of the O-type star $\zeta$ Puppis exhibit an anomalous (3G + M2)/(3F) line ratio of $\sim$1.4, in comparison with $\sim$2.4 for almost all other collisionally excited astrophysical spectra. Based on the work of Mauche et al. (2001), we conjectured that the strong UV field of $\zeta$ Puppis produces the observed ratio by depopulation of metastable $3s$ excited states, and that the ratio can potentially be used as an independent diagnostic of plasma formation radius. We used the Flexible Atomic Code (FAC) collisional-radiative model to model the effect of UV photoexcitation from O stars on the Fe XVII lines. We compared our model calculations to archival spectra of coronal and hot stars from the Chandra HETGS and XMM-Newton RGS to benchmark our calculations for various electron densities and UV field intensities. Our calculations show that UV photoexcitation does not produce a sufficiently large dynamic range in the 3F/(3F + 3G + M2) fraction to explain the difference in the observed ratio between coronal stars and $\zeta$ Pup. Thus, this effect likely cannot explain the observed line ratio of $\zeta$ Pup, and its origin is still unexplained., Comment: 21 pages, 13 figures

Context. Eta Carinae is a highly eccentric, massive binary system (semimajor axis ~15.5 au) with powerful stellar winds and a phase-dependent wind-wind collision (WWC) zone. The primary star, η Car A, is a luminous blue variable (LBV); the secondary, η Car B, is a Wolf-Rayet or O star with a faster but less dense wind. Aperture-synthesis imaging allows us to study the mass loss from the enigmatic LBV η Car. Understanding LBVs is a crucial step toward improving our knowledge about massive stars and their evolution. Aims. Our aim is to study the intensity distribution and kinematics of η Car’s WWC zone. Methods. Using the VLTI-MATISSE mid-infrared interferometry instrument, we perform Brα imaging of η Car’s distorted wind. Results. We present the first VLTI-MATISSE aperture-synthesis images of η Car A’s stellar windin several spectral channels distributed across the Brα 4.052 μm line (spectral resolving power R ~ 960). Our observations were performed close to periastron passage in February 2020 (orbital phase ~ 14.0022). The reconstructed iso-velocity images show the dependence of the primary stellar wind on wavelength or line-of-sight (LOS) velocity with a spatial resolution of 6 mas (~14 au). The radius of the faintest outer wind regions is ~26 mas (~60 au). At several negative LOS velocities, the primary stellar wind is less extended to the northwest than in other directions. This asymmetry is most likely caused by the WWC. Therefore, we see both the velocity field of the undisturbed primary wind and the WWC cavity. In continuum spectral channels, the primary star wind is more compact than in line channels. A fit of the observed continuum visibilities with the visibilities of a stellar wind CMFGEN model (CMFGEN is an atmosphere code developed to model the spectra of a variety of objects) provides a full width at half maximum fit diameter of the primary stellar wind of 2.84 ± 0.06 mas (6.54 ± 0.14 au). We comparethe derived intensity distributions with the CMFGEN stellar wind model and hydrodynamic WWC models.

We explore a simple semi-analytic model for what happens when an O star (or cluster of O stars) forms in an isolated filamentary cloud. The model is characterized by three configuration parameters: the radius of the filament, $R_{_{\rm FIL}}$, the mean density of H2 in the filament, $n_{_{\rm FIL}}$, and the rate at which the O star emits ionizing photons, $\dot{\cal N}_{_{\rm LyC}}$. We show that for a wide range of these configuration parameters, ionizing radiation from the O star rapidly erodes the filament, and the ionized gas from the filament disperses into the surroundings. Under these circumstances the distance, L, from the O star to the ionization front (IF) is given approximately by $L(t) \sim 5.2 {\rm pc} [R_{_{\rm FIL}}/0.2 {\rm pc}]^{-1/6} [n_{_{\rm FIL}}/10^4 {\rm cm^{-3}}]^{-1/3} [\dot{\cal N}_{_{\rm LyC}}/10^{49} {\rm s}^{-1}]^{1/6} [t/{\rm Myr}]^{2/3}$, and we derive similar simple power-law expressions for other quantities, for example the rate at which ionized gas boils off the filament, $\dot{M}_{_{\rm IF}}(t)$, and the mass, $M_{_{\rm SCL}}(t)$, of the shock-compressed layer that is swept up behind the IF. We show that a very small fraction of the ionizing radiation is expended locally, and a rather small amount of molecular gas is ionized and dispersed. We discuss some features of more realistic models, and the extent to which they might modify or invalidate the predictions of this idealized model. In particular we show that, for very large $R_{_{\rm FIL}}$ and/or large $n_{_{\rm FIL}}$ and/or low $\dot{\cal N}_{_{\rm LyC}}$, continuing accretion on to the filament might trap the ionizing radiation from the O star, slowing erosion of the filament even further.

In this paper, we present a sample of 27 new Galactic O-type stars found in the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) Data Release 3 (DR3). The new O-stars were identified by visual examination of spectra in a subsample of the main DR3 catalogue. The earliest of these new Galactic O-stars are two O4 stars: LS V +39 27, an O4 V((f)) star placed in the Outer Arm, probably at more than 15 kpc from the centre of the Galaxy; 2MASS J20261845+4347472, an O4 IV(f) star that has not previously been identified as OB-type. Taking into account the number of new O-stars found in the LAMOST DR3, it is very likely that more will be discovered/identified in the future. Indeed, about 63 per cent (18 of 27) of the new O-type stars are bright sources (V = 12), indicating that an important fraction of the Galactic O-stars is still unidentified, waiting to be catalogued/classified from large Galactic surveys such as LAMOST.

Hot massive stars present strong stellar winds that are driven by absorption, scattering and re\-emission of photons by the ions of the atmosphere (\textit{line-driven winds}). A better comprehension of this phenomenon, and a more accurate calculation of hydrodynamics and radiative acceleration is required to reduce the number of free parameters in spectral fitting, to determine accurate wind parameters such as mass-loss rates and velocity profiles. We use the non-LTE model-atmosphere code CMFGEN to numerically solve the radiative transfer equation in the stellar atmosphere and to calculate the radiative acceleration $g_\text{rad}(r)$. Under the assumption that the radiative acceleration depends only on the radial coordinate, we solve analytically the equation of motion by means of the Lambert $W$-function. An iterative procedure between the solution of the radiative transfer and the equation of motion is executed in order to obtain a final self-consistent velocity field that is no longer based on any $\beta$-law. We apply the Lambert-procedure to three O supergiant stars ($\zeta$-Puppis, HD~165763 and $\alpha$-Cam) and discuss the Lambert-solutions for the velocity profiles. It is found that, even without recalculation of the mass-loss rate, the Lambert-procedure allows the calculation of consistent velocity profiles that reduce the number of free parameters when a spectral fitting using CMFGEN is performed. Synthetic spectra calculated from our Lambert-solutions show significant differences compared to the initial $\beta$-law CMFGEN models. The results indicate the importance of consistent velocity profile calculation in the CMFGEN code and its usage in a fitting procedure and interpretation of observed spectra., Comment: Published in the Astrophysical Journal (ApJ), in 13 October 2021

The dispersal of protoplanetary disks sets the timescale available for planets to assemble, and thus it is one of the fundamental parameters in theories of planetary formation. Disk dispersal is determined by several properties of the central star, the disk itself, and the surrounding environment. In particular, the metallicity of disks may impact their evolution, even if to date controversial results exist: in low-metallicity clusters disks seem to rapidly disperse, while in the Magellanic Clouds some evidence supports the existence of accreting disks few tens of Myrs old. In this paper we study the dispersal timescale of disks in Dolidze~25, the young cluster in proximity of the Sun with lowest metallicity, with the aim of understanding whether disk evolution is impacted by the low-metallicity of the cluster. We have analyzed Chandra/ACIS-I observations of the cluster and combined the resulting source catalog with existing optical and infrared catalogs of the region. We selected the disk-bearing population and the disk-less population of Dolidze 25. We have derived stellar parameters from isochrones fitted to color-magnitude diagrams. We derived a disk fraction of about 34% and a median age of 1.2 Myrs. By comparing this estimate with existing estimates of the disk fraction of clusters younger than 10 Myrs, our study suggests that the disk fraction of Dolidze 25 is lower than what is expected from its age alone. Even if our results are not conclusive given the intrinsic uncertainty on stellar ages estimated from isochrones fitting to color-magnitude diagrams, we suggest that disk evolution in Dolidze 25 may be impacted by the environment. Given the poor O star population and low stellar density of the cluster, it is more likely that disks dispersal timescale is dictated more by the low metallicity of the cluster rather than external photoevaporation or dynamical encounters., Accepted for publication by Astronomy & Astrophysics

Context. Radiation-driven mass loss is key to our understanding of massive-star evolution. However, for low-luminosity O-type stars there are big discrepancies between theoretically predicted and empirically derived mass-loss rates (called the weak-wind problem). Aims. We compute radiation-line-driven wind models of a typical weak-wind star to determine its temperature structure and the corresponding impact on ultra-violet (UV) line formation. Methods. We carried out hydrodynamic simulations of the line-deshadowing instability (LDI) for a weak-wind star in the Galaxy. Subsequently, we used this LDI model as input in a short-characteristics radiative transfer code to compute synthetic UV line profiles. Results. We find that the line-driven weak wind is significantly shock heated to high temperatures and is unable to cool down efficiently. This results in a complex temperature structure where more than half of the wind volume has temperatures significantly higher than the stellar effective temperature. Therefore, a substantial portion of the weak wind will be more ionised, resulting in a reduction of the UV line opacity and therefore in weaker line profiles for a given mass-loss rate. Quantifying this, we find that weak-wind mass-loss rates derived from unsaturated UV lines could be underestimated by a factor of between 10 and 100 if the high-temperature gas is not properly taken into account in the spectroscopic analysis. This offers a tentative basic explanation for the weak-wind problem: line-driven weak winds are not really weaker than theoretically expected, but rather a large portion of their wind volume is much hotter than the stellar effective temperature.

We present the first results from a CARMA high-resolution $^{12}$CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) molecular line survey of the North America and Pelican (NAP) Nebulae. CARMA observations have been combined with single-dish data from the Purple Mountain 13.7m telescope to add short spacings and produce high-dynamic-range images. We find that the molecular gas is predominantly shaped by the W80 HII bubble that is driven by an O star. Several bright rims are probably remnant molecular clouds heated and stripped by the massive star. Matching these rims in molecular lines and optical images, we construct a model of the three-dimensional structure of the NAP complex. Two groups of molecular clumps/filaments are on the near side of the bubble, one being pushed toward us, whereas the other is moving toward the bubble. Another group is on the far side of the bubble and moving away. The young stellar objects in the Gulf region reside in three different clusters, each hosted by a cloud from one of the three molecular clump groups. Although all gas content in the NAP is impacted by feedback from the central O star, some regions show no signs of star formation, while other areas clearly exhibit star formation activity. Other molecular gas being carved by feedback includes the cometary structures in the Pelican Head region and the boomerang features at the boundary of the Gulf region. The results show that the NAP complex is an ideal place for the study of feedback effects on star formation., accepted by AJ

Dust plays a key role in the formation of planets and its emission also provides one of our most accessible views of protoplanetary discs. If set by radiative equilibrium with the central star, the temperature of dust in the disc plateaus at around $10-20$K in the outer regions. However sufficiently nearby massive stars can heat the outer disc to substantially higher temperatures. In this paper we study the radiative equilibrium temperature of discs in the presence of massive external sources and gauge the effect that it has on millimetre dust mass estimates. Since millimetre grains are not entrained in any wind we focus on geometrically simple 2D-axisymmetric disc models using radiative transfer calculations with both the host star and an external source. Recent surveys have searched for evidence of massive stars influencing disc evolution using disc properties as a function of projected separation. In assuming a disc temperature of $20$K for a disc a distance $D$ from a strong radiation source, disc masses are overestimated by a factor that scales with $D^{-1/2}$ interior to the separation that external heating becomes important. This could significantly alter dust mass estimates of discs in close proximity to $\theta^1$C in the Orion Nebular Cluster. We also make an initial assessment of the effect upon snow lines. Within a parsec of an O star like $\theta^1$C a CO snow line no longer exists, though the water snow line is virtually unaffected except for very close separations of $\leq0.01\,$pc., Comment: 11 pages, 10 figures. Accepted for publication in MNRAS

The Carina Nebula harbors a large population of high-mass stars, including at least 75 O-type and Wolf-Rayet stars, but the current census is not complete since further high-mass stars may be hidden in or behind the dense dark clouds that pervade the association. With the aim of identifying optically obscured O- and early B-type stars in the Carina Nebula, we performed the first infrared spectroscopic study of stars in the optically obscured stellar cluster Tr 16-SE, located behind a dark dust lane south of eta Car. We used the integral-field spectrograph KMOS at the ESO VLT to obtain H- and K-band spectra with a resolution of R sim 4000 (Delta lambda sim 5 A) for 45 out of the 47 possible OB candidate stars in Tr 16-SE, and we derived spectral types for these stars. We find 15 stars in Tr 16-SE with spectral types between O5 and B2 (i.e., high-mass stars with M >= 8 Msun, only two of which were known before. An additional nine stars are classified as (Ae)Be stars (i.e., intermediate-mass pre-main-sequence stars), and most of the remaining targets show clear signatures of being late-type stars and are thus most likely foreground stars or background giants unrelated to the Carina Nebula. Our estimates of the stellar luminosities suggest that nine of the 15 O- and early B-type stars are members of Tr 16-SE, whereas the other six seem to be background objects. Our study increases the number of spectroscopically identified high-mass stars (M >= 8 Msun) in Tr 16-SE from two to nine and shows that Tr 16-SE is one of the larger clusters in the Carina Nebula. Our identification of three new stars with spectral types between O5 and O7 and four new stars with spectral types O9 to B1 significantly increases the number of spectroscopically identified O-type stars in the Carina Nebula., Accepted for publication in Astronomy & Astrophysics

In the Milky Way, $\sim$18 Wolf-Rayet+O (WR+O) binaries are known with estimates of their stellar and orbital parameters. Whereas black hole+O (BH+O) binaries are thought to evolve from the former, only one such system is known in the Milky Way. To resolve this disparity, it was suggested that upon core collapse, the WR stars receive large kicks such that most of the binaries are disrupted. We reassess this issue, with emphasis on the uncertainty in the formation of an accretion disk around wind-accreting BHs in BH+O binaries, which is key to identifying such systems. We follow the methodology of previous work and apply an improved analytic criterion for the formation of an accretion disk around wind accreting BHs. We then use stellar models to predict the properties of the BH+O binaries which are expected to descend from the observed WR+O binaries, if the WR stars would form BHs without a natal kick. We find that disk formation depends sensitively on the O stars' wind velocity, the specific angular momentum carried by the wind, the efficiency of angular momentum accretion by the BH, and the spin of the BH. We show that the assumption of a low wind velocity may lead to predicting that most of the BH+O star binaries will have an extended X-ray bright period. However, this is not the case when typical wind velocities of O stars are considered. We find that a high spin of the BH can boost the duration of the X-ray active phase as well as the X-ray brightness during this phase, producing a strong bias for detecting high mass BH binaries in X-rays with high BH spin parameters. We conclude that large BH formation kicks are not required to understand the sparsity of X-ray bright BH+O stars in the Milky Way. Probing for a population of X-ray silent BH+O systems with alternative methods can inform us about BH kicks and the conditions for high energy emission from high mass BH binaries. (Abridged), Comment: 14 pages, 9 figures, accepted for publication 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.

The Hubble Space Telescope (HST) survey Measuring Young Stars in Space and Time (MYSST) entails some of the deepest photometric observations of extragalactic star formation, capturing even the lowest mass stars of the active star-forming complex N44 in the Large Magellanic Cloud. We employ the new MYSST stellar catalog to identify and characterize the content of young pre-main-sequence (PMS) stars across N44 and analyze the PMS clustering structure. To distinguish PMS stars from more evolved line of sight contaminants, a non-trivial task due to several effects that alter photometry, we utilize a machine learning classification approach. This consists of training a support vector machine (SVM) and a random forest (RF) on a carefully selected subset of the MYSST data and categorize all observed stars as PMS or non-PMS. Combining SVM and RF predictions to retrieve the most robust set of PMS sources, we find $\sim26,700$ candidates with a PMS probability above 95% across N44. Employing a clustering approach based on a nearest neighbor surface density estimate, we identify 18 prominent PMS structures at $1$ $\sigma$ significance above the mean density with sub-clusters persisting up to and beyond $3$ $\sigma$ significance. The most active star-forming center, located at the western edge of N44's bubble, is a subcluster with an effective radius of $\sim 5.6$ pc entailing more than 1,100 PMS candidates. Furthermore, we confirm that almost all identified clusters coincide with known H II regions and are close to or harbor massive young O stars or YSOs previously discovered by MUSE and Spitzer observations., Comment: 29 pages, 21 figures, accepted for publication in AJ

This paper reports 209 O-type stars found with LAMOST. All 135 new O-type stars discovered so far with LAMOST so far are given. Among them, 94 stars are firstly presented in this sample. There are 1 Iafpe star, 5 Onfp stars, 12 Oe stars, 1 Ofc stars, 3 ON stars, 16 double-lined spectroscopic binaries, and 33 single-lined spectroscopic binaries. All O-type stars are determined based on LAMOST low-resolution spectra (R ~ 1800), with their LAMOST median-resolution spectra (R~7500) as supplements., Comment: 37 pages, 22 figures, 1 table

Metal-poor massive stars dominate the light we observe from star-forming dwarf galaxies and may have produced the bulk of energetic photons that reionized the universe at high redshift. Yet, the rarity of observations of individual O stars below the $20\%$ solar metallicity ($Z_\odot$) of the Small Magellanic Cloud (SMC) hampers our ability to model the ionizing fluxes of metal-poor stellar populations. We present new Hubble Space Telescope far-ultraviolet (FUV) spectra of three O-dwarf stars in the galaxies Leo P ($3\%\,Z_\odot$), Sextans A ($6\%\,Z_\odot$), and WLM ($14\%\,Z_\odot$). We quantify equivalent widths of photospheric metal lines and strengths of wind-sensitive features, confirming that both correlate with metallicity. We infer the stars' fundamental properties by modeling their FUV through near-infrared spectral energy distributions and identify stars in the SMC with similar properties to each of our targets. Comparing to the FUV spectra of the SMC analogs suggests that (1) the star in WLM has an SMC-like metallicity, and (2) the most metal-poor star in Leo P is driving a much weaker stellar wind than its SMC counterparts. We measure projected rotation speeds and find that the two most metal-poor stars have high $v \,\mathrm{sin}(i)\,\geq\,290\,\mathrm{km}\,\mathrm{s}^{-1}$, and estimate just a $3-6\%$ probability of finding two fast rotators if the metal-poor stars are drawn from the same $v \,\mathrm{sin}(i)$ distribution observed for O dwarfs in the SMC. These observations suggest that models should include the impact of rotation and weak winds on ionizing flux to accurately interpret observations of metal-poor galaxies in both the near and distant universe., Comment: Accepted at ApJ. 9 figures, 24 pages

We report an observational study of the Galactic HII region Sh 2-305/S305 using the [CII] 158 $\mu$m line data, which are used to examine the gas dynamics and structure of photodissociation regions. The integrated [CII] emission map at [39.4, 49.5] km s$^{-1}$ spatially traces two shell-like structures (i.e., inner and outer neutral shells) having a total mass of $\sim$565 M$_\odot$. The inner neutral shell encompasses an O9.5V star at its centre and has a compact ring-like appearance. However, the outer shell is seen with more extended and diffuse [CII] emission, hosting an O8.5V star at its centre and surrounds the inner neutral shell. The velocity channel maps and position-velocity diagrams confirm the presence of a compact [CII] shell embedded in the diffuse outer shell, and both the shells seem to expand with $v_{\rm exp}\sim$1.3 km s$^{-1}$. The outer shell appears to be older than the inner shell, hinting that these shells are formed sequentially. The [CII] profiles are examined toward S305, which are either double-peaked or blue-skewed and have the brighter redshifted component. The redshifted and blueshifted components spatially trace the inner and outer neutral shell geometry, respectively. The ionized, neutral, and molecular zones in S305 are seen adjacent to one another around the O-type stars. The regularly spaced dense molecular and dust clumps (mass $\sim$10-10$^{3}$ M$_{\odot}$) are investigated around the neutral shells, which might have originated due to gravitational instability in the shell of collected materials., Comment: 16 pages, 8 figures, Accepted for publication in The Astrophysical Journal

The Bajamar Star is an early O star that ionizes the North America/Pelican Nebulae. In projection, it is near the geometric center of the H II region, but appears to lie outside any of the main stellar subgroups. Furthermore, in Gaia DR2, there were slight discrepancies between this star and the rest of the system in parallax (2$\sigma$ larger) and relative tangential velocity (~6 km/s). Using Gaia EDR3, we find that the parallax discrepancy has disappeared, but the velocity difference remains. These results are consistent with the star having escaped from a subgroup., Comment: 4 pages and 1 figure. Submitted to Research Notes of the AAS

FEEDBACK is a SOFIA (Stratospheric Observatory for Infrared Astronomy) legacy program dedicated to study the interaction of massive stars with their environment. It performs a survey of 11 galactic high mass star-forming regions in the 158 μm (1.9 THz) line of [C II] and the 63 μm (4.7 THz) line of [O I]. We employ the 14 pixel Low Frequency Array and 7 pixel High Frequency Array upGREAT heterodyne instrument to spectrally resolve (0.24 MHz) these far-infrared fine structure lines. With a total observing time of 96h, we will cover ∼6700 arcmin2 at 14 1) angular resolution for the [C II] line and 6 3 for the [O I] line. The observations started in spring 2019 (Cycle 7). Our aim is to understand the dynamics in regions dominated by different feedback processes from massive stars such as stellar winds, thermal expansion, and radiation pressure, and to quantify the mechanical energy injection and radiative heating efficiency. This is an important science topic because feedback of massive stars on their environment regulates the physical conditions and sets the emission characteristics in the interstellar medium (ISM), influences the star formation activity through molecular cloud dissolution and compression processes, and drives the evolution of the ISM in galaxies. The [C II] line provides the kinematics of the gas and is one of the dominant cooling lines of gas for low to moderate densities and UV fields. The [O I] line traces warm and high-density gas, excited in photodissociations regions with a strong UV field or by shocks. The source sample spans a broad range in stellar characteristics from single OB stars, to small groups of O stars, to rich young stellar clusters, to ministarburst complexes. It contains well-known targets such as Aquila, the Cygnus X region, M16, M17, NGC7538, NGC6334, Vela, and W43 as well as a selection of H II region bubbles, namely RCW49, RCW79, and RCW120. These [C II] maps, together with the less explored [O I] 63 μm line, provide an outstanding database for the community. They will be made publically available and will trigger further studies and follow-up observations.

We present a review of UV observations of massive stars and their analysis. We discuss O stars, luminous blue variables, and Wolf–Rayet stars. Because of their effective temperature, the UV (912−3200 Å) provides invaluable diagnostics not available at other wavebands. Enormous progress has been made in interpreting and analysing UV data, but much work remains. To facilitate the review, we provide a brief discussion on the structure of stellar winds, and on the different techniques used to model and interpret UV spectra. We discuss several important results that have arisen from UV studies including weak-wind stars and the importance of clumping and porosity. We also discuss errors in determining wind terminal velocities and mass-loss rates.

Initial results are presented from 3D MHD modelling of stellar-wind bubbles around O stars moving supersonically through the ISM. We describe algorithm updates that enable high-resolution 3D MHD simulations at reasonable computational cost. We apply the methods to the simulation of the astrosphere of a rotating massive star moving with 30 km/s through the diffuse interstellar medium, for two different stellar magnetic field strengths, 10 G and 100 G. Features in the flow are described and compared with similar models for the Heliosphere. The shocked interstellar medium becomes asymmetric with the inclusion of a magnetic field misaligned with the star's direction of motion, with observable consequences. When the Alfvenic Mach number of the wind is $\leq$10 then the stellar magnetic field begins to affect the structure of the wind bubble and features related to the magnetic axis of the star become visible at parsec scales. Prospects for predicting and measuring non-thermal radiation are discussed., Accepted for publication in Journal of Physics Conference Series, proceedings of the 19th Annual International Astrophysics Conference (Santa Fe, New Mexico, 9-13 March 2020)

With the launch of JWST and the upcoming installation of extremely large telescopes, the first galaxies in our Universe will finally be revealed. Their light will be dominated by massive stars, which peak in in the ultra-violet (UV) part of the electromagnetic spectrum. Star formation is the key driver of the evolution of our Universe. At young ages, within 10 Million years, both high and low mass stars generate complex UV emission processes which are poorly understood yet are vital for interpreting high red-shift line emission. For these reasons, the Hubble Space Telescope (HST) will devote 1000 orbits to obtaining a UV Legacy Library of Young Stars as Essential Standards (ULLYSES). The purpose of this Overview is to outline the basic physical principles driving UV emission processes from local (within 100 parsecs of) star formation, ranging from huge star-forming complexes containing hundreds of massive and very-massive stars (VMS), such as 30 Doradus (the Tarantula Nebula) in the neighbouring Magellanic Clouds (only 50 kpc away), to galaxies near and far, out to the epoch of Cosmic Reionization., 12 pages, 5 figures. Comments welcome. Accepted Editorial Overview for Special Issue in Galaxies: "Star Formation in the Ultraviolet", Ed. Jorick Vink

Context. As a building block for amino acids, formamide (NH$_2$CHO) is an important molecule in astrobiology and astrochemistry, but its formation path in the interstellar medium is not understood well. Aims. We aim to find empirical evidence to support the chemical relationships of formamide to HNCO and H$_2$CO. Methods. We examine high angular resolution (~0.2") Atacama Large Millimeter/submillimeter Array (ALMA) maps of six sources in three high-mass star-forming regions and compare the spatial extent, integrated emission peak position, and velocity structure of HNCO and H$_2$CO line emission with that of NH$_2$CHO by using moment maps. Through spectral modeling, we compare the abundances of these three species. Results. In these sources, the emission peak separation and velocity dispersion of formamide emission is most often similar to HNCO emission, while the velocity structure is generally just as similar to H$_2$CO and HNCO (within errors). From the spectral modeling, we see that the abundances between all three of our focus species are correlated, and the relationship between NH$_2$CHO and HNCO reproduces the previously demonstrated abundance relationship. Conclusions. In this first interferometric study, which compares two potential parent species to NH$_2$CHO, we find that all moment maps for HNCO are more similar to NH$_2$CHO than H$_2$CO in one of our six sources (G24 A1). For the other five sources, the relationship between NH$_2$CHO, HNCO, and H$_2$CO is unclear as the different moment maps for each source are not consistently more similar to one species as opposed to the other., Main text: 15 pages, 17 figures, 10 tables. Appendix: 21 pages, 28 figures

Some massive, merging black holes (BH) may be descendants of binary O stars. The evolution and mass transfer between these O stars determines the spins of their progeny BH. These will be measurable with future gravitational wave detectors, incentivizing the measurement of the spins of O stars in binaries. We previously measured the spins of O stars in Galactic Wolf-Rayet (WR) + O binaries. Here we measure the vsini of four LMC and two SMC O stars in WR + O binaries to determine whether lower metallicity might affect the spin rates. We find that the O stars in Galactic and Magellanic WR + O binaries display average vsini = 258 +/- 18 km/s and 270 +/- 15 km/s, respectively. Two LMC O stars measured on successive nights show significant line width variability, possibly due to differing orbital phases exhibiting different parts of the O stars illuminated differently by their WR companions. Despite this variability, the vsini are highly super-synchronous but distinctly subcritical for the O stars in all these binaries; thus we conclude that an efficient mechanism for shedding angular momentum from O stars in WR + O binaries must exist. This mechanism, probably related to Roche lobe overflow-created dynamo magnetic fields, prevents nearly 100% breakup spin rates, as expected when RLOF operates, as it must, in these stars. A Spruit-Tayler dynamo and O star wind might be that mechanism., Comment: 8 pages, 13 figures, accepted in MNRAS

Our Galactic center contains young stars, including the few million year old clockwise disk between 0.05 and 0.5 pc from the Galactic center, and the S-star cluster of B-type stars at a galactocentric distance of ~0.01 pc. Recent observations suggest the S-stars are remnants of tidally disrupted binaries from the clockwise disk. In particular, Koposov et al. 2020 discovered a hypervelocity star that was ejected from the Galactic center 5 Myr ago, with a velocity vector consistent with the disk. We perform a detailed study of this binary disruption scenario. First, we quantify the plausible range of binary semimajor axes in the disk. Dynamical evaporation of such binaries is dominated by other disk stars rather than by the isotropic stellar population. For the expected range of semimajor axes in the disk, binary tidal disruptions can reproduce the observed S-star semimajor axis distribution. Reproducing the observed thermal eccentricity distribution of the S-stars requires an additional relaxation process. The flight time of the Koposov star suggests that this process must be effective within 10 Myr. We consider three possibilities: (i) scalar resonant relaxation from the observed isotropic star cluster, (ii) torques from the clockwise disk, and (iii) an intermediate-mass black hole. We conclude that the first and third mechanisms are fast enough to reproduce the observed S-star eccentricity distribution. Finally, we show that the primary star from an unequal-mass binary would be deposited at larger semimajor axes than the secondary, possibly explaining the dearth of O stars among the S-stars., 21 pages, 16 figures + Appendix. Accepted for publication in ApJ. The new version is substantially improved with a more detailed treatment of resonant and nonresonant relaxation from the observed isotropic cluster in the Galactic center. In contrast to the first version, we now conclude that these processes alone can reproduce the S-stars' eccentricity distribution in less than 10 Myr

Mass loss through stellar winds plays a dominant role in the evolution of massive stars. In particular the mass-loss rates of very massive stars (VMSs, $> 100\,M_{\odot}$) are highly uncertain. Such stars display Wolf-Rayet spectral morphologies (WNh) whilst on the main-sequence. Metal-poor VMSs are progenitors of gamma-ray bursts and pair instability supernovae. In this study we extended the widely used stellar wind theory by Castor, Abbott & Klein from the optically thin (O star) to the optically thick main-sequence (WNh) wind regime. In particular we modify the mass-loss rate formula in a way that we are able to explain the empirical mass-loss dependence on the Eddington parameter ($\Gamma_{\rm e}$). The new mass-loss recipe is suitable for incorporation into current stellar evolution models for massive and very massive stars. It makes verifiable predictions, namely how the mass-loss rate scales with metallicity and at which Eddington parameter the transition from optically thin O star to optically thick WNh star winds occurs. In the case of the star cluster R136 in the Large Magellanic Cloud we find in the optically thin wind regime $\dot{M} \propto \Gamma_{\rm e}^{3}$ while in the optically thick wind regime $\dot{M} \propto 1/ (1 - \Gamma_{\rm e})^{3.5}$. The transition from optically thin to optically thick winds occurs at $\Gamma_{\rm e, trans} \approx 0.47$. The transition mass-loss rate is $\log \dot{M}~(M_{\odot} \mathrm{yr}^{-1}) \approx -4.76 \pm 0.18$, which is in line with the prediction by Vink & Gr\"afener assuming a volume filling factor of $f_{\rm V} = 0.23_{-0.15}^{+0.40}$., Comment: Accepted for publication in MNRAS, 9 pages, 5 figures

Microquasars are neutron star or black hole X-ray binaries with jets. These jets can create shock-ionized bubbles of hot plasma that can masquerade as peculiar supernova remnants (SNRs) in extragalactic surveys. To see if this is the case in the well-studied spiral galaxy M83, where one microquasar candidate (M83-MQ1) has already been identified, we studied the properties of nine SNR candidates, selected because of their elongated or irregular morphology, from the set of previously identified SNRs in that galaxy. Using multiwavelength data from Chandra, the Hubble Space Telescope, Gemini, and the Australia Telescope Compact Array, we found that at least six of our nine sources are best interpreted as SNRs. For one source, we do not have enough observational data to explain its linear morphology. Another source shows a nebular optical spectrum dominated by photo-ionization by O stars, but its excess [Fe II] and radio luminosity suggest a possible hidden SNR. Finally, one source (S2) shows an elongated structure of ionized gas, two radio sources along that line, and an accretion-powered X-ray source in between them (the Chandra source L14-139). While S2 could be a chance alignment of multiple SNRs and one X-ray binary, it seems more likely that it is a single physical structure powered by the jet from the accreting compact object. In the latter case, the size and luminosity of S2 suggest a kinetic power of around 10^{40} erg/s, in the same class as the most powerful microquasars in the local universe (e.g., S26 in NGC7793 and SS433 in our own Galaxy)., Comment: Accepted by ApJ. Size = 8 MB

Galactic Wolf-Rayet (WR) star membership of star forming regions can be used to constrain the formation environments of massive stars. Here, we utilise $\textit{Gaia}$ DR2 parallaxes and proper motions to reconsider WR star membership of clusters and associations in the Galactic disk, supplemented by recent near-IR studies of young massive clusters. We find that only 18$-$36% of 553 WR stars external to the Galactic Centre region are located in clusters, OB associations or obscured star-forming regions, such that at least 64% of the known disk WR population are isolated, in contrast with only 13% of O stars from the Galactic O star Catalogue. The fraction located in clusters, OB associations or star-forming regions rises to 25$-$41% from a global census of 663 WR stars including the Galactic Centre region. We use simulations to explore the formation processes of isolated WR stars. Neither runaways, nor low mass clusters, are numerous enough to account for the low cluster membership fraction. Rapid cluster dissolution is excluded as mass segregation ensures WR stars remain in dense, well populated environments. Only low density environments consistently produce WR stars that appeared to be isolated during the WR phase. We therefore conclude that a significant fraction of WR progenitors originate in low density association-like surroundings which expand over time. We provide distance estimates to clusters and associations host to WR stars, and estimate cluster ages from isochrone fitting., Comment: 19 pages, 7 figures. Accepted for publication in MNRAS

NGC 7538 IRS 1 is a very young embedded O star driving an ionized jet and accreting mass with an accretion rate > 10^-4 Msun/year, which is quenching the hypercompact HII region. We use SOFIA GREAT data, Herschel PACS and SPIRE archive data, SOFIA FORCAST archive data, Onsala 20m and CARMA data, and JCMT archive data to determine the properties of the O star and its outflow. IRS 1 appears to be a single O-star with a bolometric luminosity > 1 10^5 Lsun, i.e. spectral type O7 or earlier. We find that IRS 1 drives a large molecular outflow with the blue-shifted northern outflow lobe extending to ~ 280" or 3.6 pc from IRS 1. Near IRS 1 the outflow is well aligned with the ionized jet. The dynamical time scale of the outflow is ~ 1.3 10^5 yr. The total outflow mass is ~ 130 Msun. We determine a mass outflow rate of 1.0 10^-3 Msun/yr, roughly consistent with the observed mass accretion rate. We observe strong high velocity [CII] emission in the outflow, confirming that strong UV radiation from IRS 1 escapes into the outflow lobes and is ionizing the gas. Many O stars may form like low mass stars, but with a higher accretion rate and in a denser environment. As long as the accretion stays high enough to quench the HII region, the star will continue to grow. When the accretion rate drops, the HII region will rapidly start to expand., 21 pages, 19 figures. Accepted for publication in The Astrophysical Journal

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. The apparent lack of massive O-type stars near the zero-age main sequence, or ZAMS (at ages, With funding from the Spanish government through the "María de Maeztu Unit of Excellence" accreditation (MDM-2017-0737)

The Magellanic Bridge stretching between the SMC and LMC is the nearest tidally stripped intergalactic environment and has a low average metallicity of $Z~0.1Z_{\odot}$. Here we report the first discovery of three O-type stars in the Bridge using archival spectra collected with FLAMES at ESO/VLT. We analyze the spectra using the PoWR models, which provide the physical parameters, ionizing photon fluxes, and surface abundances. This discovery suggests that the tidally stripped low density gas is capable of producing massive O stars and their ages imply ongoing star formation in the Bridge. The multi-epoch spectra indicate that all three O stars are binaries. Despite their spatial proximity to each other, these O stars are chemically distinct. One of them is a fast-rotating giant with nearly LMC-like abundances. The other two are main-sequence stars that rotate extremely slowly and are strongly metal depleted. This includes the most nitrogen-poor O star known up to date. Taking into account the previous analyses of B stars in the Bridge, we interpret the various metal abundances as the signature of a chemically inhomogeneous interstellar medium, suggesting that the gas might have accreted during multiple episodes of tidal interaction between the Clouds. Attributing the lowest derived metal content to the primordial gas, the time of initial formation of the Bridge may date back to several Gyr. Using the Gaia and Galex color-magnitude diagrams we roughly estimate the total number of O stars in the Bridge and their total ionizing radiation. Comparing with the energetics of the diffuse ISM, we find that the contribution of the hot stars to the ionizing radiation field in the Bridge is less than 10%, and conclude that the main sources of ionizing photons are leaks from the LMC and SMC. We estimate a lower limit for the fraction of ionizing radiation that escapes from these two dwarf galaxies., Comment: to be published 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)

We report the results of our study of the high-velocity (\approx160 km/s) runaway O star CPD-64 2731 and its associated horseshoe-shaped nebula discovered with the Wide-field Infrared Survey Explorer. Spectroscopic observations with the Southern African Large Telescope and spectral analysis indicate that CPD-64 2731 is a fast-rotating main-sequence O5.5 star with enhanced surface nitrogen abundance. We derive a projected rotational velocity of \approx300 km/s which is extremely high for this spectral type. Its kinematic age of \approx6 Myr, assuming it was born near the Galactic plane, exceeds its age derived from single star models by a factor of two. These properties suggest that CPD-64 2731 is a rejuvenated and spun-up binary product. The geometry of the nebula and the almost central location of the star within it argue against a pure bow shock interpretation for the nebula. Instead, we suggest that the binary interaction happened recently, thereby creating the nebula, with a cavity blown by the current fast stellar wind. This inference is supported by our results of 2D numerical hydrodynamic modelling., Comment: 15 pages, 11 figures, accepted for publication in MNRAS

In order to test massive star evolution above 25 M⊙, we perform spectral analysis on a sample of massive stars in the Small Magellanic Cloud that includes both O stars as well as more evolved Wolf-Rayet stars. We present a grid of non-LTE stellar atmospheres that has been calculated using the cmfgen code, in order to have a systematic and homogeneous approach. We obtain stellar and wind parameters for O stars, spectral types ranging from O2 to O6, and the complete sample of known Wolf-Rayet stars. We discuss the evolutionary status of both the O and WR stars and the links between them, as well as the most likely evolutionary path towards black hole formation in a low metallicity environment, including testing theoretical predictions for mass-loss rates at low metallicities.

A small subset of Galactic O-stars possess surface magnetic fields that alter the outflowing stellar wind by magnetically confining it. Key to the magnetic confinement is that it induces rotational modulation of spectral lines over the full EM domain; this allows us to infer basic quantities, e.g., mass-loss rate and magnetic geometry. Here, we present an empirical study of the H$\alpha$ line in Galactic magnetic O-stars to constrain the mass fed from the stellar base into the magnetosphere, using realistic multi-dimensional magnetized wind models, and compare with theoretical predictions. Our results suggest that it may be reasonable to use mass-feeding rates from non-magnetic wind theory if the absolute mass-loss rate is scaled down according to the amount of wind material falling back upon the stellar surface. This provides then some empirical support to the proposal that such magnetic O-stars might evolve into heavy stellar-mass black holes (Petit et al. 2017), Comment: 4 pages, 3 figures, IAU Symposium 346 proceedings

We examine claims of strong gravity-darkening effects in the O-star components of WR+O binaries. We generate synthetic spectra for a wide range of parameters, and show that the line-width results are consistent with extensive measurements of O stars that are either single or are members of `normal' binaries. By contrast, the WR+O results are at the extremes of, or outside, the distributions of both models and other observations. Remeasurement of the WR+O spectra shows that they can be reconciled with other results by judicious choice of pseudo-continuum normalization. With this interpretation, the supersynchronous rotation previously noted for the O-star components in the WR+O binaries with the longest orbital periods appears to be unexceptional. Our investigation is therefore consistent with the aphorism that if the title of a paper ends with a question mark, the answer is probably `no'., Accepted in MNRAS

WR 25 is a colliding-wind binary star system comprised of a very massive O2.5If*/WN6 primary and an O-star secondary in a 208-day period eccentric orbit. These hot stars have strong, highly-supersonic winds which interact to form a bright X-ray source from wind-collision-shocks whose conditions change with stellar separation. Different views through the WR and O star winds are afforded with orbital phase as the stars move about their orbits, allowing for exploration of wind structure in ways not easy or even possible for single stars. We have analyzed an on-axis Chandra/HETGS spectrum of WR 25 obtained shortly before periastron when the X-rays emanating from the system are the brightest. From the on-axis observations, we constrain the line fluxes, centroids, and widths of various emission lines, including He-triplets of Si XIII and Mg XI. We have also been able to include several serendipitous off-axis HETG spectra from the archive and study their flux variation with phase. This is the first report on high-resolution spectral studies of WR 25 in X-rays.

Hot star winds are driven by the radiative force due to light absorption in lines of heavier elements. Therefore, the amount of mass lost by the star per unit of time, i.e., the mass-loss rate, is sensitive to metallicity. We provide mass-loss rate predictions for O stars with mass fraction of heavier elements 0.2 ⊙ ≤ 1. Our predictions are based on global model atmospheres. The models allow us to predict wind terminal velocity and the mass-loss rate just from basic global stellar parameters. We provide a formula that fits the mass-loss rate predicted by our models as a function of stellar luminosity and metallicity. On average, the mass-loss rate scales with metallicity as (Z/Z⊙)0.59. The predicted mass-loss rates agree with mass-loss rates derived from ultraviolet wind line profiles. At low metallicity, the rotational mixing affects the wind mass-loss rates. We study the influence of magnetic line blanketing.

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\%., Accepted for publication in A&A. The arXiv version includes appendix C, which is an online only figure in A&A

We report the identification from multi-wavelength observations of the Fermi Large Area Telescope (LAT) source 4FGL J1405.1-6119 (= 3FGL J1405.4-6119) as a high-mass gamma-ray binary. Observations with the LAT show that gamma-ray emission from the system is modulated at a period of 13.7135 +/- 0.0019 days, with the presence of two maxima per orbit with different spectral properties. X-ray observations using the Neil Gehrels Swift Observatory X-ray Telescope (XRT) show that X-ray emission is also modulated at this period, but with a single maximum that is closer to the secondary lower-energy gamma-ray maximum. A radio source, coincident with the X-ray source, is also found from Australia Telescope Compact Array (ATCA) observations, and the radio emission is modulated on the gamma-ray period with similar phasing to the X-ray emission. A large degree of interstellar obscuration severely hampers optical observations, but a near-infrared counterpart is found. Near-infrared spectroscopy indicates an O6 III spectral classification. This is the third gamma-ray binary to be discovered with the Fermi LAT from periodic modulation of the gamma-ray emission, the other two sources also have early O star, rather than Be star, counterparts. We consider at what distances we can detect such modulated gamma-ray emission with the LAT, and examine constraints on the gamma-ray binary population of the Milky Way., Comment: Accepted for publication in the Astrophysical Journal