Search

Star formation has often been studied by separating the low- and high-mass regimes with an approximate boundary at 8M_sun. While some of the outcomes of the star-formation process are different between the two regimes, it is less clear whether the physical processes leading to these outcomes are that different at all. Here, we systematically compare low- and high-mass star formation by reviewing the most important processes and quantities from an observational and theoretical point of view. We identify three regimes where processes are either similar, quantitatively or qualitatively different between low- and high-mass star formation. Similar characteristics can be identified for the turbulent gas properties and density structures of the star-forming regions. Many of the observational characteristics also do not depend that strongly on the environment. Quantitative differences can be found for outflow, infall and accretion rates as well as mean column and volume densities. Also the multiplicity significantly rises from low- to high-mass stars. The importance of the magnetic field for the formation processes appears still less well constrained. Qualitative differences between low- and high-mass star formation relate mainly to the radiative and ionizing feedback that occurs almost exclusively in regions forming high-mass stars. Nevertheless, accretion apparently can continue via disk structures in ionized accretion flows. Finally, we discuss to what extent a unified picture of star formation over all masses is possible and which issues need to be addressed in the future., Comment: 45 pages, 6 figures, accepted for Annual Reviews of Astronomy and Astrophysics

Massive stars lead the ionization and mechanical feedback within young star-forming regions. The Large Magellanic Cloud (LMC) is an ideal galaxy for studying individual massive stars and quantifying their feedback contribution to the environment. We analyze eight exemplary targets in LMC N11 B from the Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) program, using novel spectra from HST (COS and STIS) in the UV, and from VLT (X-shooter) in the optical. We model the spectra of early to late O-type stars by using state-of-the-art PoWR atmosphere models. We determine the stellar and wind parameters (e.g., $T_\star$, $\log g$, $L_{\star}$, $\dot{M}$, $v_\infty$) of the analyzed objects, chemical abundances (C, N, O), ionizing and mechanical feedback ($Q_\mathrm{H}$, $Q_\mathrm{He{\small{I}}}$, $Q_\mathrm{He{\small{II}}}$, $L_\mathrm{mec}$) and X-rays. We report ages of $2-4.5$ Myr and masses of $30-60$ $M_\odot$ for the analyzed stars in N11 B, consistent with a scenario of sequential star formation. We note that the observed wind-momentum luminosity relation is consistent with theoretical predictions. We detect nitrogen enrichment in most of the stars, up to a factor of seven. However, we do not find a correlation between nitrogen enrichment and projected rotational velocity. Finally, based on their spectral type, we estimate the total ionizing photons injected from the O-type stars in N11 B into its environment. We report $\log$ ($\sum$ $Q_\mathrm{H}$)$=50.5$ ph s$^{-1}$, $\log$ ($\sum$ $Q_\mathrm{He{\small{I}}}$)$=49.6$ ph s$^{-1}$ and $\log$ ($\sum$ $Q_\mathrm{He{\small{II}}}$)$=44.4$ ph s$^{-1}$, consistent with the total ionizing budget in N11., Comment: Accepted in Astronomy and Astrophysics

Mass loss through a stellar wind is an important physical process that steers the evolution of massive stars and controls the properties of their end-of-life products, such as the supernova type and the mass of compact remnants. For an accurate mass loss determination, the inhomogeneities in the wind, known as clumping, needs to be taking into account. We aim to improve empirical estimates of mass loss and wind clumping for hot main-sequence massive stars, study the dependence of both properties on the metallicity, and compare the theoretical predictions to our findings. We analyzed the optical and UV spectra of 13 O-type stars in the Small Magellanic Cloud galaxy, which has a metallicity of $\sim 0.2\,Z_\odot$. We quantified the stellar atmosphere, outflow, and wind-clumping properties. To probe the role of metallicity, we compared our findings to studies of Galactic and Large Magellanic Cloud samples that were analyzed with similar methods. We find significant variations in the wind-clumping properties of the target stars, with clumping starting at flow velocities $0.01 - 0.23$ of the terminal wind velocity and reaching clumping factors $f_{\rm cl} = 2 - 30$. In the luminosity ($\log L / L_{\odot} = 5.0 - 6.0$) and metallicity ($Z/Z_{\odot} = 0.2 - 1$) range we considered, we find that the scaling of the mass loss $\dot{M}$ with metallicity $Z$ varies with luminosity. At $\log L/L_{\odot} = 5.75$, we find $\dot{M} \propto Z^m$ with $m = 1.02 \pm 0.30$, in agreement with pioneering work in the field. For lower luminosities, however, we obtain a significantly steeper scaling of $m > 2$. The monotonically decreasing $m(L)$ behavior adds a complexity to the functional description of the mass-loss rate of hot massive stars. Although the trend is present in the predictions, it is much weaker than we found here., Comment: 37 pages, 29 figures. Accepted for publication in Astronomy and Astrophysics

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.

We use data from the ALMA Evolutionary Study of High Mass Protocluster Formation in the Galaxy (ALMAGAL) survey to study 100 ALMAGAL regions at $\sim$ 1 arsecond resolution located between $\sim$ 2 and 6 kpc distance. Using ALMAGAL $\sim$ 1.3mm line and continuum data we estimate flow rates onto individual cores. We focus specifically on flow rates along filamentary structures associated with these cores. Our primary analysis is centered around position velocity cuts in H$_2$CO (3$_{0,3}$ - 2$_{0,2}$) which allow us to measure the velocity fields, surrounding these cores. Combining this work with column density estimates we derive the flow rates along the extended filamentary structures associated with cores in these regions. We select a sample of 100 ALMAGAL regions covering four evolutionary stages from quiescent to protostellar, Young Stellar Objects (YSOs), and HII regions (25 each). Using dendrogram and line analysis, we identify a final sample of 182 cores in 87 regions. In this paper, we present 728 flow rates for our sample (4 per core), analysed in the context of evolutionary stage, distance from the core, and core mass. On average, for the whole sample, we derive flow rates on the order of $\sim$10$^{-4}$ M$_{sun}$yr$^{-1}$ with estimated uncertainties of $\pm$50%. We see increasing differences in the values among evolutionary stages, most notably between the less evolved (quiescent/protostellar) and more evolved (YSO/HII region) sources. We also see an increasing trend as we move further away from the centre of these cores. We also find a clear relationship between the flow rates and core masses $\sim$M$^{2/3}$ which is in line with the result expected from the tidal-lobe accretion mechanism. Overall, we see increasing trends in the relationships between the flow rate and the three investigated parameters; evolutionary stage, distance from the core, and core mass., Comment: 11 pages, 11 figures, accepted for publication in A&A

Context. B supergiants (BSGs) represent an important connection between the main sequence and more extreme evolutionary stages of massive stars. Additionally, lying toward the cool end of the hot star regime, determining their wind properties is crucial to constrain the evolution and feedback of massive stars as, for instance, they might manifest the bi-stability jump phenomenon. Aims. We undertake a detailed analysis of a representative sample of 18 Small Magellanic Cloud (SMC) BSGs within the ULLYSES and XShootU datasets. Our UV and optical analysis spans BSGs from B0 to B8 - covering the bi-stability jump region. We aim to evaluate their evolutionary status and verify what their wind properties say about the bi-stability jump in a low-metallicity environment. Methods. We used the CMFGEN to model the spectra and photometry (from UV to infrared) of our sample. We compare our results with different evolutionary models, with previous determinations in the literature of OB stars, and with diverging mass-loss recipes at the bi-stability jump. Additionally, we provide the first BSG models in the SMC including X-rays. Results. (i) Within a single-stellar evolution framework, the evolutionary status of early BSGs seem less clear than that of late BSGs, which agree with H-shell burning models. (ii) UV analysis shows evidence that BSGs contain X-rays in their atmospheres, for which we provide constraints. In general, we find higher X-ray luminosity (close to the standard log(L_X/L) ~ -7) for early BSGs. For cooler BSGs, lower values are preferred, log(L_X/L) ~ -8.5. (iii) The obtained mass-loss rates suggest neither a jump nor a monotonic decrease with temperature. Instead, a rather constant trend is observed, which is at odds with the increase found for Galactic BSGs. (iv) The wind velocity behavior with temperature shows a sharp drop at ~19 kK, similar to what is observed for Galactic BSGs., Comment: 33 pages (23+10)+(22 at Zenodo), 34 figures (21+13)+(21 at Zenodo), 7 tables (3+4)+(1 at Zenodo), accepted for publication

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

Methods: Observing the large pc-scale Stokes I mm dust continuum emission with the IRAM 30m telescope and the intermediate-scale (<0.1pc) polarized submm dust emission with the Submillimeter Array toward a sample of 20 high-mass star-forming regions allows us to quantify the dependence of the fragmentation behaviour of these regions depending on the density and magnetic field structures. Results: We infer density distributions n~r^{-p} of the regions with typical power-law slopes p around ~1.5. There is no obvious correlation between the power-law slopes of the density structures on larger clump scales (~1pc) and the number of fragments on smaller core scales (<0.1pc). Comparing the large-scale single-dish density profiles to those derived earlier from interferometric observations at smaller spatial scales, we find that the smaller-scale power-law slopes are steeper, typically around ~2.0. The flattening toward larger scales is consistent with the star-forming regions being embedded in larger cloud structures that do not decrease in density away from a particular core. Regarding the magnetic field, for several regions it appears aligned with filamentary structures leading toward the densest central cores. Furthermore, we find different polarization structures with some regions exhibiting central polarization holes whereas other regions show polarized emission also toward the central peak positions. Nevertheless, the polarized intensities are inversely related to the Stokes I intensities. We estimate magnetic field strengths between ~0.2 and ~4.5mG, and we find no clear correlation between magnetic field strength and the fragmentation level of the regions. Comparison of the turbulent to magnetic energies shows that they are of roughly equal importance in this sample. The mass-to-flux ratios range between ~2 and ~7, consistent with collapsing star-forming regions., Comment: Accepted for Astronomy & Astrophysics, 14 pages, 14 figures plus appendices, also download option at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

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

The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in the 1.37 mm wavelength regime in both continuum and line emission, reaching $\sim$0.4" resolution (800 au at 2 kpc). Using the dense gas tracer CH$_3$CN, we find velocity gradients across 13 cores perpendicular to the directions of bipolar molecular outflows, making them excellent disk candidates. Specific angular momentum ($j$) radial profiles are on average $\sim10^{-3}$ km /s pc and follow $j \propto r^{1.7}$, consistent with a poorly resolved rotating and infalling envelope/disk model. Fitting the velocity profiles with a Keplerian model, we find protostellar masses in the range of $\sim 10-25$ $M_\odot$. Modelling the level population of CH$_3$CN lines, we present temperature maps and find median gas temperatures in the range $70-210$ K. We create Toomre $Q$ maps to study the stability of the disks and find almost all (11 of 13) disk candidates to be prone to fragmentation due to gravitational instabilities at the scales probed by our observations. In particular, disks with masses greater than $\sim10-20\%$ of the mass of their host (proto)stars are Toomre unstable, and more luminous protostellar objects tend to have disks that are more massive and hence more prone to fragmentation. Our finings show that most disks around high-mass protostars are prone to disk fragmentation early in their formation due to their high disk to stellar mass ratio. This impacts the accretion evolution of high-mass protostars which will have significant implications for the formation of the most massive stars., Comment: 27 pages, 12 figures, 6 appendices - accepted for publication in Astronomy and Astrophysics

The origin of the observed population of Wolf-Rayet (WR) stars in low-metallicity (low-Z) galaxies, such as the Small Magellanic Cloud (SMC), is not yet understood. Standard, single-star evolutionary models predict that WR stars should stem from very massive O-type star progenitors, but these are very rare. On the other hand, binary evolutionary models predict that WR stars could originate from primary stars in close binaries. We conduct an analysis of the massive O star, AzV 14, to spectroscopically determine its fundamental and stellar wind parameters, which are then used to investigate evolutionary paths from the O-type to the WR stage with stellar evolutionary models. Multi-epoch UV and optical spectra of AzV 14 are analyzed using the non-LTE stellar atmosphere code PoWR. An optical TESS light curve was extracted and analyzed using the PHOEBE code. The obtained parameters are put into an evolutionary context, using the MESA code. AzV 14 is a close binary system consisting of two similar main sequence stars with masses of 32 Msol. Both stars have weak stellar winds with mass-loss rates of log $\dot{M}$ = -7.7. Binary evolutionary models can explain the empirically derived stellar and orbital parameters. The model predicts that the primary will evolve into a WR star with T = 100 kK, while the secondary, which will accrete significant amounts of mass during the first mass transfer phase, will become a cooler WR star with T = 50 kK and are predicted to have compared to other WR stars increased oxygen abundances. This model prediction is supported by a spectroscopic analysis of a WR star in the SMC. We hypothesize that the populations of WR stars in low-Z galaxies may have bimodal temperature distributions. Hotter WR stars might originate from primary stars, while cooler WR stars are the evolutionary descendants of the secondary stars if they accreted a significant amount of mass., Comment: 21 pages (13 main body + 8 appendix), 16 figures, 9 tables

Recent observations found close-in planets with significant atmospheres of hydrogen and helium in great abundance. These are the so-called super-Earths and mini-Neptunes. Their atmospheric composition suggests that they formed early during the gas-rich phase of the circumstellar disk and were able to avoid becoming hot Jupiters. As a possible explanation, recent studies explored the recycling hypothesis and showed that atmosphere-disk recycling is able to fully compensate for radiative cooling and thereby halt Kelvin-Helmholtz contraction to prevent runaway gas accretion. To understand the parameters that determine the efficiency of atmospheric recycling, we extend our earlier studies by exploring the effects of the core mass, the effect of circumstellar gas on sub-Keplerian orbits (headwind), and the optical depth of the surrounding gas on the recycling timescale. Additionally, we analyze their effects on the size and mass of the forming atmosphere. For the explored parameter space, all simulations eventually reach an equilibrium where heating due to hydrodynamic recycling fully compensates radiative cooling. In this equilibrium, the atmosphere-to-core mass ratio stays well below $10 \, \%$, preventing the atmosphere from becoming self-gravitating and entering runaway gas accretion. Higher core masses cause the atmosphere to become turbulent, which further enhances recycling. Even for our highest core mass of $10 \, M_\mathrm{Earth}$, atmosphere-disk recycling is efficient enough to fully compensate for radiative cooling and prevent the atmosphere from becoming self-gravitating. Hence, in-situ formation of hot Jupiters is very unlikely, and migration of gas giants is a key process required to explain their existence. Our findings imply that atmosphere-disk recycling is the most natural explanation for the prevalence of close-in super-Earths and mini-Neptunes., Comment: 18 pages, 13 figures, Accepted for publication by A&A on 2/7/2022