Search

The gravitationally lensed star WHL 0137-LS, nicknamed Earendel, was identified with a photometric redshift z (phot) = 6.2 +/- 0.1 based on images taken with the Hubble Space Telescope. Here we present James Webb Space Telescope (JWST) Near Infrared Camera images of Earendel in eight filters spanning 0.8-5.0 mu m. In these higher-resolution images, Earendel remains a single unresolved point source on the lensing critical curve, increasing the lower limit on the lensing magnification to mu > 4000 and restricting the source plane radius further to r < 0.02 pc, or similar to 4000 au. These new observations strengthen the conclusion that Earendel is best explained by an individual star or multiple star system and support the previous photometric redshift estimate. Fitting grids of stellar spectra to our photometry yields a stellar temperature of T (eff) similar to 13,000-16,000 K, assuming the light is dominated by a single star. The delensed bolometric luminosity in this case ranges from log(L)=5.8 L-theta, which is in the range where one expects luminous blue variable stars. Follow-up observations, including JWST NIRSpec scheduled for late 2022, are needed to further unravel the nature of this object, which presents a unique opportunity to study massive stars in the first billion years of the universe.

We present a study of the dust-to-gas ratios in five nearby galaxies: NGC 628 (M74), NGC 6503, NGC 7793, UGC 5139 (Holmberg I), and UGC 4305 (Holmberg II). Using Hubble Space Telescope broadband WFC3/UVIS UV and optical images from the Treasury program Legacy ExtraGalactic UV Survey (LEGUS) combined with archival HST/Advanced Camera for Surveys data, we correct thousands of individual stars for extinction across these five galaxies using an isochrone-matching (reddening-free Q) method. We generate extinction maps for each galaxy from the individual stellar extinctions using both adaptive and fixed resolution techniques and correlate these maps with neutral H i and CO gas maps from the literature, including the H i Nearby Galaxy Survey and the HERA CO-Line Extragalactic Survey. We calculate dust-to-gas ratios and investigate variations in the dust-to-gas ratio with galaxy metallicity. We find a power-law relationship between dust-to-gas ratio and metallicity, consistent with other studies of dust-to-gas ratio compared to metallicity. We find a change in the relation when H2 is not included. This implies that underestimation of in low-metallicity dwarfs from a too-low CO-to-H2 conversion factor X CO could have produced too low a slope in the derived relationship between dust-to-gas ratio and metallicity. We also compare our extinctions to those derived from fitting the spectral energy distribution (SED) using the Bayesian Extinction and Stellar Tool for NGC 7793 and find systematically lower extinctions from SED fitting as compared to isochrone matching.

The Legacy ExtraGalactic UV Survey (LEGUS) is a multiwavelength Cycle 21 Treasury program on the Hubble Space Telescope. It studied 50 nearby star-forming galaxies in 5 bands from the near-UV to the I-band, combining new Wide Field Camera 3 observations with archival Advanced Camera for Surveys data. LEGUS was designed to investigate how star formation occurs and develops on both small and large scales, and how it relates to the galactic environments. In this paper we present the photometric catalogs for all the apparently single stars identified in the 50 LEGUS galaxies. Photometric catalogs and mosaicked images for all filters are available for download. We present optical and near-UV color-magnitude diagrams for all the galaxies. For each galaxy we derived the distance from the tip of the red giant branch. We then used the NUV color-magnitude diagrams to identify stars more massive than 14 M o, and compared their number with the number of massive stars expected from the GALEX FUV luminosity. Our analysis shows that the fraction of massive stars forming in star clusters and stellar associations is about constant with the star formation rate. This lack of a relation suggests that the timescale for evaporation of unbound structures is comparable or longer than 10 Myr. At low star formation rates this translates to an excess of mass in clustered environments as compared to model predictions of cluster evolution, suggesting that a significant fraction of stars form in unbound systems

We report the large effort that is producing comprehensive high-level young star cluster (YSC) catalogs for a significant fraction of galaxies observed with the Legacy ExtraGalactic UV Survey (LEGUS) Hubble treasury program. We present the methodology developed to extract cluster positions, verify their genuine nature, produce multiband photometry (from NUV to NIR), and derive their physical properties via spectral energy distribution fitting analyses. We use the nearby spiral galaxy NGC 628 as a test case for demonstrating the impact that LEGUS will have on our understanding of the formation and evolution of YSCs and compact stellar associations within their host galaxy. Our analysis of the cluster luminosity function from the UV to the NIR finds a steepening at the bright end and at all wavelengths suggesting a dearth of luminous clusters. The cluster mass function of NGC 628 is consistent with a power-law distribution of slopes and a truncation of a few times 105. After their formation, YSCs and compact associations follow different evolutionary paths. YSCs survive for a longer time frame, confirming their being potentially bound systems. Associations disappear on timescales comparable to hierarchically organized star-forming regions, suggesting that they are expanding systems. We find mass-independent cluster disruption in the inner region of NGC 628, while in the outer part of the galaxy there is little or no disruption. We observe faster disruption rates for low mass (≤104) clusters, suggesting that a mass-dependent component is necessary to fully describe the YSC disruption process in NGC 628.

We present the star cluster catalogues for 17 dwarf and irregular galaxies in the HST Treasury Program 'Legacy ExtraGalactic UV Survey' (LEGUS). Cluster identification and photometry in this sub-sample are similar to that of the entire LEGUS sample, but special methods were developed to provide robust catalogues with accurate fluxes due to low cluster statistics. The colours and ages are largely consistent for two widely used aperture corrections, but a significant fraction of the clusters are more compact than the average training cluster. However, the ensemble luminosity, mass, and age distributions are consistent suggesting that the systematics between the two methods are less than the random errors. When compared with the clusters from previous dwarf galaxy samples, we find that the LEGUS catalogues are more complete and provide more accurate total fluxes. Combining all clusters into a composite dwarf galaxy, we find that the luminosity and mass functions can be described by a power law with the canonical index of -2 independent of age and global SFR binning. The age distribution declines as a power law, with an index of approximate to -0.80 +/- 0.15, independent of cluster mass and global SFR binning. This decline of clusters is dominated by cluster disruption since the combined star formation histories and integrated-light SFRs are both approximately constant over the last few hundred Myr. Finally, we find little evidence for an upper-mass cut-off (< 2 sigma) in the composite cluster mass function, and can rule out a truncation mass below approximate to 10(4.5)M(circle dot) but cannot rule out the existence of a truncation at higher masses.

The recent gravitational wave measurements have demonstrated the existence of stellar mass black hole binaries. It is essential for our understanding of massive star evolution to identify the contribution of binary evolution to the formation of double black holes. A promising way to progress is investigating the progenitors of double black hole systems and comparing predictions with local massive star samples such as the population in 30 Doradus in the Large Magellanic Cloud (LMC). Methods. To this purpose, we analyse a large grid of detailed binary evolution models at LMC metallicity with initial primary masses between 10 and 40 Msun, and identify which model systems potentially evolve into a binary consisting of a black hole and a massive main sequence star. We then derive the observable properties of such systems, as well as peculiarities of the OB star component. We find that about 3% of the LMC late O and early B stars in binaries are expected to possess a black hole companion, when assuming stars with a final helium core mass above 6.6 M to form black holes. While the vast majority of them may be X-ray quiet, our models suggest that these may be identified in spectroscopic binaries, either by large amplitude radial velocity variations ( > 50 km s ) and simultaneous nitrogen surface enrichment, or through a moderate radial velocity ( > 10 km/s ) and simultaneously rapid rotation of the OB star. The predicted mass ratios are such that main sequence companions could be excluded in most cases. A comparison to the observed OB+WR binaries in the LMC, Be/X-ray binaries, and known massive BH binaries supports our conclusion. We expect spectroscopic observations to be able to test key assumptions in our models, with important implications for massive star evolution in general, and for the formation of double-black hole mergers in particular., Comment: 19 pages, 14 figures, Astronomy and Astrophysics, in press

A high fraction of massive stars are found to be binaries but only a few of them are reported as photometrically variable. By studying the populations of SB2 in the 30 Doradus region, we found a subset of them that have photometry from the OGLE project and that display variations in their light curves related to orbital motions. The goal of this study is to determine the dynamical masses and radii of the 26 binary components to investigate the mass-discrepancy problem and to provide an empirical mass-luminosity relation for the LMC. We use the PHOEBE programme to perform a systematic analysis of the OGLE V and I light curves obtained for 13 binary systems in 30 Dor. We adopt Teff, and orbital parameters derived previously to obtain the inclinations of the systems and the parameters of the individual components. Three systems display eclipses in their light curves, while the others only display ellipsoidal variations. We classify two systems as over-contact, five as semi-detached, and four as detached. The two remaining systems have uncertain configurations due to large uncertainties on their inclinations. The fact that systems display ellipsoidal variations has a significant impact on the inclination errors. From the dynamical masses, luminosities, and radii, we provide LMC-based empirical mass-luminosity and mass-radius relations, and we compare them to other relations given for the Galaxy, the LMC, and the SMC. These relations differ for different mass ranges, but do not seem to depend on the metallicity regimes. We also compare the dynamical, spectroscopic, and evolutionary masses of the stars in our sample. While the dynamical and spectroscopic masses agree with each other, the evolutionary masses are systematically higher, at least for stars in semi-detached systems. This suggests that the mass discrepancy can be partly explained by past or ongoing interactions between the stars., Comment: Accepted to A&A The core of the paper 9 pages, 5 pages of Appendix

Accurate stellar parameters of individual objects in binary systems are essential to constrain the effects of binarity on stellar evolution. These parameters serve as a prerequisite to probing existing and future theoretical evolutionary models. We aim to derive the atmospheric parameters of the 31 SB2s in the TMBM sample. This sample, composed of detached, semi-detached and contact systems with at least one of the components classified as an O star, is an excellent test-bed to study how binarity can impact our knowledge of the evolution of massive stars. 32 epochs of FLAMES/GIRAFFE spectra are analysed using spectral disentangling to construct the individual spectra of 62 components. We apply the CMFGEN atmosphere code to determine their stellar parameters and their He, C and N surface abundances. From these properties, we show that the effects of tides on chemical mixing are limited. Components on longer-period orbits show higher nitrogen enrichment at their surface than those on shorter-period orbits, in contrast to expectations of rotational or tidal mixing, implying that other mechanisms play a role in this process. Components filling their Roche lobe are mass donors. They exhibit higher nitrogen content at their surface and rotate more slowly than their companions. By accreting new material, their companions spin faster and are rejuvenated. Their locations in the N-vsini diagram tend to show that binary products are good candidates to populate the two groups of stars (slowly rotating, nitrogen-enriched and rapidly rotating non-enriched) that cannot be reproduced through single-star population synthesis. This sample is the largest sample of binaries to be studied in such a homogeneous way. The study of these objects gives us strong observational constraints to test theoretical binary evolutionary tracks., Comment: accepted to A&A Core of the paper 15 pages, Appendix: 41 pages

A majority of massive stars are part of binary systems, a large fraction of which will inevitably interact during their lives. Binary-interaction products (BiPs), i.e. stars affected by such interaction, are expected to be commonly present in stellar populations. BiPs are thus a crucial ingredient in the understanding of stellar evolution. We aim to identify and characterize a statistically significant sample of BiPs by studying clusters of 10-40 Myr, an age at which binary population models predict the abundance of BiPs to be highest. One example of such a cluster is NGC 330 in the Small Magellanic Cloud. Using MUSE WFM-AO observations of NGC 330, we resolve the dense cluster core for the first time and are able to extract spectra of its entire massive star population. We develop an automated spectral classification scheme based on the equivalent widths of spectral lines in the red part of the spectrum. We characterize the massive star content of the core of NGC 330 which contains more than 200 B stars, 2 O stars, 6 A-type supergiants and 11 red supergiants. We find a lower limit on the Be star fraction of $32 \pm 3\%$ in the whole sample. It increases to at least $46 \pm 10\%$ when only considering stars brighter than V=17 mag. We estimate an age of the cluster core between 35 and 40 Myr and a total cluster mass of $88^{+17}_{-18} \times 10^3 M_{\odot}$. We find that the population in the cluster core is different than the population in the outskirts: while the stellar content in the core appears to be older than the stars in the outskirts, the Be star fraction and the observed binary fraction are significantly higher. Furthermore, we detect several BiP candidates that will be subject of future studies., Comment: 19 pages (incl. appendix), 21 figures, 3 tables, accepted for publication in A&A

Massive stars are often found in binary systems and it has been argued that binary products boost the ionizing radiation of stellar populations. Accurate predictions for binary products are needed to understand and quantify their contribution to Cosmic Reionization. We investigate the contribution of stars stripped in binaries since (1) they are, arguably, the best-understood products of binary evolution, (2) we recently produced the first non-LTE radiative transfer calculations for the atmospheres of these stripped stars that predict their ionizing spectra, and (3) they are very promising sources since they boost the ionizing emission of stellar populations at late times. This allows stellar feedback to clear the surroundings such that a higher fraction of their photons can escape and ionize the intergalactic medium. Combining our detailed predictions for the ionizing spectra with a simple cosmic reionization model, we estimate that stripped stars contributed tens of percent of the photons that caused cosmic reionization of hydrogen, depending on the assumed escape fractions. More importantly, stripped stars harden the ionizing emission. We estimate that the spectral index for the ionizing part of the spectrum can increase to -1 compared to <-2 for single stars. At high redshift, stripped stars and massive single stars combined dominate the HeII-ionizing emission, but we expect active galactic nuclei drive cosmic helium reionization. Further observational consequences we expect are (1) high ionization states for the intergalactic gas surrounding stellar systems, such as CIV and SiIV and (2) additional heating of the intergalactic medium of up to a few thousand Kelvin. Quantifying these warrants the inclusion of accurate models for stripped stars and other binary products in full cosmological simulations., Comment: Under review in A&A, suggestions welcome

Gravitational-wave detections are now starting to probe the mass distribution of stellar-mass black holes (BHs). Robust predictions from stellar models are needed to interpret these. Theory predicts the existence of a gap in the BH mass distribution because of pair-instability supernova. The maximum BH mass below the gap is the result of pulsational mass loss. We evolve massive helium stars through their late hydrodynamical phases of evolution using the open-source MESA stellar evolution code. We find that the location of the lower edge of the mass gap at 45$M_\odot$ is remarkably robust against variations in the metallicity ($\approx 3M_\odot$), the treatment of internal mixing ($\approx 1M_\odot$), stellar wind mass loss ($\approx 4M_\odot$), making it the most robust predictions for the final stages of massive star evolution. The reason is that the onset of the instability is dictated by the near-final core mass, which in turn sets the resulting BH mass. However, varying $^{12}C\left(\alpha,\gamma\right)^{16}O$ reaction rate within its $1\sigma$ uncertainties shifts the location of the gap between $40M_\odot$ and $56M_\odot$. We provide updated analytic fits for population synthesis simulations. Our results imply that the detection of merging BHs can provide constraints on nuclear astrophysics. Furthermore, the robustness against metallicity suggests that there is a universal maximum for the location of the lower edge of the gap, which is insensitive to the formation environment and redshift for first-generation BHs. This is promising for the possibility to use the location of the gap as a "standard siren" across the Universe., Comment: 17 pages, 5 figures, 1 Table, Accepted ApJ

Stars stripped of their envelopes from interaction with a binary companion emit a significant fraction of their radiation as ionizing photons. They are potentially important stellar sources of ionizing radiation, however, they are still often neglected in spectral synthesis simulations or simulations of stellar feedback. We modeled the radiative contribution from stripped stars by using detailed evolutionary and spectral models. We estimated their impact on the integrated spectra and specifically on the emission rates of HI-, HeI-, and HeII-ionizing photons from stellar populations. We find that stripped stars have the largest impact on the ionizing spectrum of a population in which star formation halted several Myr ago. In such stellar populations, stripped stars dominate the emission of ionizing photons, mimicking a younger stellar population in which massive stars are still present. Our models also suggest that stripped stars have harder ionizing spectra than massive stars. The additional ionizing radiation affects observable properties that are related to the ionizing emission from stellar populations. In co-eval stellar populations, the ionizing radiation from stripped stars increases the ionization parameter and the production efficiency of HI-ionizing photons. They also cause high values for these parameters for about ten times longer than what is predicted for massive stars. The hard ionizing radiation from stripped stars likely introduces a characteristic ionization structure of the nebula, which leads to the emission of highly ionized elements such as O$^{2+}$ and C$^{3+}$. We, therefore, expect that the presence of stripped stars affects the location in the BPT diagram and the diagnostic ratio of OIII to OII nebular emission lines. Our models are publicly available through CDS database and on the Starburst99 website., Comment: Accepted for publication in Astronomy & Astrophysics

Cosmic History has witnessed the lives and deaths of multiple generations of massive stars, all of them invigorating their host galaxies with ionizing photons, kinetic energy, fresh material and stellar-mass black holes. Ubiquitous engines as they are, Astrophysics needs a good understanding of their formation, evolution, properties and yields throughout the history of the Universe, and with decreasing metal content mimicking the environment at the earliest epochs. Ultimately, a physical model that could be extrapolated to zero metallicity would enable tackling long-standing questions such as "What did the First, very massive stars of the Universe look like?" or "What was their role in the re-ionization of the Universe?". Yet, most our knowledge of metal-poor massive stars is drawn from one single point in metallicity. Massive stars in the Small Magellanic Cloud (SMC, $\sim 1/5 Z_{\odot}$) currently serve as templates for low-metallicity objects in the early Universe, even though significant differences with respect to massive stars with poorer metal content have been reported. This White Paper summarizes the current knowledge on extremely (sub-SMC) metal poor massive stars, highlighting the most outstanding open questions and the need to supersede the SMC as standard. A new paradigm can be built from nearby extremely metal-poor galaxies that make a new metallicity ladder, but massive stars in these galaxies are out of reach to current observational facilities. Such task would require an L-size mission, consisting of a 10m-class space telescope operating in the optical and the ultraviolet ranges. Alternatively, we propose that ESA unites efforts with NASA to make the LUVOIR mission concept a reality, thus continuing the successful partnership that made Hubble Space Telescope one of the greatest observatories of all time., Comment: A white paper submitted for the Voyage 2050 long-term plan in the ESA Science Programme. 21 pages, 1 table, 6 figures. arXiv admin note: substantial text overlap with arXiv:1903.05235