Showing 1–42 of 42 results for author: Gotberg, Y

Binary Yellow Supergiants in the Magellanic Clouds I: Photometric Candidate Identification

Authors: Anna J. G. O'Grady, Maria R. Drout, Kathryn F. Neugent, Bethany Ludwig, Ylva Gotberg, Bryan M. Gaensler

Abstract: Recent works have constrained the binary fraction of evolved populations of massive stars in local galaxies such as red supergiants and Wolf-Rayet stars, but the binary fraction of yellow supergiants (YSGs) in the Hertzsprung Gap remains unconstrained. Binary evolution theory predicts that the Hertzsprung Gap is home to multiple populations of binary systems with varied evolutionary histories. In… ▽ More Recent works have constrained the binary fraction of evolved populations of massive stars in local galaxies such as red supergiants and Wolf-Rayet stars, but the binary fraction of yellow supergiants (YSGs) in the Hertzsprung Gap remains unconstrained. Binary evolution theory predicts that the Hertzsprung Gap is home to multiple populations of binary systems with varied evolutionary histories. In this paper, we develop a method to distinguish single YSGs from YSG plus O- or B-type main sequence binaries using optical and ultraviolet photometry, and then apply this method to identify candidate YSG binaries in the Magellanic Clouds. After constructing a set of combined stellar atmosphere models, we find that optical photometry is, given typical measurement and reddening uncertainties, sufficient to discern single YSGs from YSG+OB binaries if the OB-star is at least $\sim5$M$_{\odot}$ for T$_{\mathrm{eff,YSG}}\sim$ 4000 K, but requires a $\sim$20M$_{\odot}$ OB star for YSGs up to T$_{\mathrm{eff,YSG}}\sim$ 9000 K. For these hotter YSG temperatures, ultraviolet photometry allows binaries with OB companions as small as $\sim$7M$_{\odot}$ to be identified. We use color-color spaces developed from these models to search for evidence of excess blue or ultraviolet light in a set of $\sim$1000 YSG candidates in the Magellanic Clouds. We identify hundreds of candidate YSG binary systems and report a preliminary fraction of YSGs that show a blue/UV color excess of 20-60\%. Spectroscopic follow-up is now required to confirm the true nature of this population. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: 20 pages, 10 figures, 3 tables. Submitted to The Astrophysical Journal on June 24th 2024

Tidal Disruption Events from Stripped Stars

Authors: Brenna Mockler, Monica Gallegos-Garcia, Ylva Götberg, Jon Miller, Enrico Ramirez-Ruiz

Abstract: Observations of tidal disruption events (TDEs) show signs of Nitrogen enrichment reminiscent of other astrophysical sources such as active galactic nuclei (AGN) and star-forming galaxies. Given that TDEs probe the gas from a single star, it is possible to test if the observed enrichment is consistent with expectations from the CNO cycle by looking at the observed Nitrogen/Carbon (N/C) abundance ra… ▽ More Observations of tidal disruption events (TDEs) show signs of Nitrogen enrichment reminiscent of other astrophysical sources such as active galactic nuclei (AGN) and star-forming galaxies. Given that TDEs probe the gas from a single star, it is possible to test if the observed enrichment is consistent with expectations from the CNO cycle by looking at the observed Nitrogen/Carbon (N/C) abundance ratios. Given that $\approx 20\%$ of solar mass stars (and an even larger fraction of more massive stars) live in close binaries, it is worthwhile to also consider what TDEs from stars influenced by binary evolution would look like. We show here that TDEs from stars stripped of their Hydrogen-rich (and Nitrogen-poor) envelopes through previous binary-induced mass loss can produce much higher observable N/C enhancements than even TDEs from massive stars. Additionally, we predict that the time-dependence of the N/C abundance ratio in the mass fallback rate of stripped stars will follow the inverse behavior of main-sequence stars, enabling a more accurate characterization of the disrupted star. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 11 pages, 5 figures, submitted to ApJL

X-Shooting ULLYSES: Massive Stars at Low Metallicity

Authors: Jorick S. Vink, Paul Crowther, Alex Fullerton, Miriam Garcia, Fabrice Martins, Nidia Morrell, Lida Oskinova, Nicole St. Louis, Asif ud-Doula, Andreas Sander, Hugues Sana, Jean-Claude Bouret, Brankica Kubatova, Pablo Marchant, Lucimara P. Martins, Aida Wofford, Jacco van Loon, O. Grace Telford, Ylva Götberg, Dominic Bowman, Christi Erba, Venu Kalari, The XShootU Collaboration

Abstract: The Hubble Space Telescope has devoted 500 orbits to observing 250 massive stars with low metallicity in the ultraviolet (UV) range within the framework of the ULLYSES program. The X-Shooting ULLYSES (XShootU) project enhances the legacy value of this UV dataset by providing high-quality optical and near-infrared spectra, which are acquired using the wide-wavelength-coverage X-shooter spectrograph… ▽ More The Hubble Space Telescope has devoted 500 orbits to observing 250 massive stars with low metallicity in the ultraviolet (UV) range within the framework of the ULLYSES program. The X-Shooting ULLYSES (XShootU) project enhances the legacy value of this UV dataset by providing high-quality optical and near-infrared spectra, which are acquired using the wide-wavelength-coverage X-shooter spectrograph at ESO's Very Large Telescope. XShootU emphasises the importance of combining UV with optical spectra for the consistent determination of key stellar parameters such as effective temperature, surface gravity, luminosity, abundances, and wind characteristics including mass-loss rates as a function of metallicity. Since uncertainties in these parameters have implications across various branches of astrophysics, the data and modelling generated by the XShootU project are poised to significantly advance our understanding of massive stars at low metallicity. This is particularly crucial for confidently interpreting JWST data of the earliest stellar generations, making XShootU a unique resource for comprehending individual spectra of low-metallicity stars. △ Less

Submitted 30 April, 2024; originally announced May 2024.

Comments: 6 pages, 6 figures. ESO Large Programme Overview

Journal ref: ESO Messenger, 2024

The young massive SMC cluster NGC 330 seen by MUSE III. Stellar parameters and rotational velocities

Authors: J. Bodensteiner, H. Sana, P. L. Dufton, C. Wang, N. Langer, G. Banyard, L. Mahy, A. de Koter, S. E. de Mink, C. J. Evans, Y. Götberg, V. Hénault-Brunet, L. R. Patrick, F. R. N. Schneider

Abstract: The origin of initial rotation rates of stars, and how a star's surface rotational velocity changes during the evolution, either by internal angular momentum transport or due to interactions with a binary companion, remain open questions in stellar astrophysics. Here, we aim to derive the physical parameters and study the distribution of (projected) rotational velocities of B-type stars in the 35… ▽ More The origin of initial rotation rates of stars, and how a star's surface rotational velocity changes during the evolution, either by internal angular momentum transport or due to interactions with a binary companion, remain open questions in stellar astrophysics. Here, we aim to derive the physical parameters and study the distribution of (projected) rotational velocities of B-type stars in the 35 Myr-old, massive cluster NGC 330 in the Small Magellanic Cloud. NGC 330 is in an age range where the number of post-interaction binaries is predicted to be high near the cluster turnoff (TO). We develop a simultaneous photometric and spectroscopic grid-fitting method adjusting atmosphere models on multi-band Hubble Space Telescope photometry and Multi Unit Spectroscopic Explorer spectroscopy. This allows us to homogeneously constrain the physical parameters of over 250 B and Be stars, brighter than mF814W = 18.8 mag. The rotational velocities of Be stars in NGC 330 are significantly higher than the ones of B stars. The rotational velocities vary as a function of the star's position in the color-magnitude diagram, qualitatively following predictions of binary population synthesis. A comparison to younger clusters shows that stars in NGC 330 rotate more rapidly on average. The rotational velocities of the 35 Myr old population in NGC 330 quantitatively agree with predictions for a stellar population that underwent significant binary interactions: the bulk of the B stars could be single stars or primaries in pre-interaction binaries. The rapidly spinning Be stars could be mass and angular momentum gainers in previous interactions, while those Be stars close to the TO may be spun-up single stars. The slowly rotating, apparently single stars above the TO could be merger products. The different vsini-characteristics of NGC 330 compared to younger populations can be understood in this framework. △ Less

Submitted 28 August, 2023; originally announced August 2023.

Comments: 18 pages (incl. appendix), 15 figures, 3 tables, accepted for publication in A&A

Journal ref: A&A 680, A32 (2023)

The B-type Binaries Characterisation Programme II. VFTS 291: A stripped star from a recent mass transfer phase

Authors: J. I. Villaseñor, D. J. Lennon, A. Picco, T. Shenar, P. Marchant, N. Langer, P. L. Dufton, F. Nardini, C. J. Evans, J. Bodensteiner, S. E. de Mink, Y. Götberg, I. Soszyński, W. D. Taylor, H. Sana

Abstract: Recent studies of massive binaries with putative black hole companions have uncovered a phase of binary evolution that has not been observed before, featuring a bloated stripped star that very recently ceased transferring mass to a main-sequence companion. In this study, we focus on the candidate system VFTS 291, a binary with an orbital period of 108 d and a high semi-amplitude velocity (… ▽ More Recent studies of massive binaries with putative black hole companions have uncovered a phase of binary evolution that has not been observed before, featuring a bloated stripped star that very recently ceased transferring mass to a main-sequence companion. In this study, we focus on the candidate system VFTS 291, a binary with an orbital period of 108 d and a high semi-amplitude velocity ($K_1=93.7\pm0.2$ km s$^{-1}$). Through our analysis of the disentangled spectra of the two components, together with dynamical and evolutionary arguments, we identify a narrow-lined star of ~1.5-2.5 $M_\odot$ dominating the spectrum, and an early B-type main-sequence companion of $13.2\pm1.5$ $M_\odot$. The low mass of the narrow-lined star, and the high mass ratio, suggest that VFTS 291 is a post-mass-transfer system, with the narrow-lined star being bloated and stripped of its hydrogen-rich envelope, sharing many similarities with other recently discovered stripped stars. Our finding is supported by our detailed binary evolution models, which indicate that the system can be well explained by an initial configuration consisting of an 8.1 $M_\odot$ primary with an 8 $M_\odot$ companion in a 7 d orbital period. While some open questions remain, particularly concerning the surface helium enrichment of the stripped star and the rotational velocity of the companion, we expect that high-resolution spectroscopy may help reconcile our estimates with theory. Our study highlights the importance of multi-epoch spectroscopic surveys to identify and characterize binary interaction products, and provides important insights into the evolution of massive binary stars. △ Less

Submitted 11 September, 2023; v1 submitted 15 July, 2023; originally announced July 2023.

Comments: 25 pages, 17 figures. Published by MNRAS, Advance Access

Stellar properties of observed stars stripped in binaries in the Magellanic Clouds

Authors: Y. Gotberg, M. R. Drout, A. P. Ji, J. H. Groh, B. A. Ludwig, P. A. Crowther, N. Smith, A. de Koter, S. E. de Mink

Abstract: Massive stars (~8-25Msun) stripped of their hydrogen-rich envelopes via binary interaction are thought to be the main progenitors for merging neutron stars and stripped-envelope supernovae. We recently presented the discovery of the first set of such stripped stars in a companion paper. Here, we fit the spectra of ten stars with new atmosphere models in order to constrain their stellar properties… ▽ More Massive stars (~8-25Msun) stripped of their hydrogen-rich envelopes via binary interaction are thought to be the main progenitors for merging neutron stars and stripped-envelope supernovae. We recently presented the discovery of the first set of such stripped stars in a companion paper. Here, we fit the spectra of ten stars with new atmosphere models in order to constrain their stellar properties precisely. We find that the stellar properties align well with the theoretical expectations from binary evolution models for helium-core burning envelope-stripped stars. The fits confirm that the stars have high effective temperatures (Teff~50-100kK), high surface gravities (log g ~5), and hydrogen-poor/helium-rich surfaces (X(H, surf)~0-0.4) while showing for the first time a range of bolometric luminosities (10^3-10^5 Lsun), small radii (~0.5-1Rsun), and low Eddington factors (Gamma_e~0.006-0.4). Using these properties, we derive intermediate current masses (~1-8Msun), which suggest that their progenitors were massive stars (~5-25Msun) and that a subset will reach core-collapse, leaving behind neutron stars or black holes. Using the model fits, we also estimate the emission rates of ionizing photons for these stars, which agree well with previous model expectations. Further, by computing models for a range of mass-loss rates, we find that the stellar winds are weaker than predicted by any existing scheme (Mdot(wind)<~ 1e-9 Msun/yr). The properties of this first sample of intermediate mass helium stars suggest they both contain progenitors of type Ib and IIb supernovae, and provide important benchmarks for binary evolution and population synthesis models. △ Less

Submitted 30 June, 2023; originally announced July 2023.

Comments: Accepted for publication in ApJ. Spectral models available via zenodo as described in text. Please also see a companion manuscript detailing the discovery of these systems in today's arXiv posting

Discovery of the missing intermediate-mass helium stars stripped in binaries

Authors: M. R. Drout, Y. Götberg, B. A. Ludwig, J. H. Groh, S. E. de Mink, A. J. G. O'Grady, N. Smith

Abstract: The theory of binary evolution predicts that many massive stars should lose their hydrogen-rich envelopes via interaction with a companion -- revealing hot helium stars with masses of $\sim$2--8M$_{\odot}$. However, only one candidate system had been identified, leaving a large discrepancy between theory and observation. Here, we present a new sample of stars -- identified via excess ultraviolet e… ▽ More The theory of binary evolution predicts that many massive stars should lose their hydrogen-rich envelopes via interaction with a companion -- revealing hot helium stars with masses of $\sim$2--8M$_{\odot}$. However, only one candidate system had been identified, leaving a large discrepancy between theory and observation. Here, we present a new sample of stars -- identified via excess ultraviolet emission -- whose luminosities, colors, and spectral morphologies are consistent with predictions for the missing population. We detect radial velocity variations indicative of binary motion and measure high temperatures ($T_{\rm eff}\sim60-100$kK), high surface gravities ($\log(g)\sim5$) and depleted surface hydrogen mass fractions ($X_{\rm{H,surf}}\lesssim0.3$), which match expectations for stars with initial masses between 8--25 M$_{\odot}$ that have been stripped via binary interaction. These systems fill the helium star mass gap between subdwarfs and Wolf-Rayet stars, and are thought to be of large astrophysical significance as ionizing sources, progenitors of stripped-envelope supernovae and merging double neutron stars. △ Less

Submitted 30 June, 2023; originally announced July 2023.

Comments: Original version of manuscript submitted in August 2022, per journal policies. Manuscript passed through referee process and was recommended for publication in Science in March 2023. Conclusions have not changed. Please contact authors for further information. Data and updated manuscript will be made available upon publication. Please also see a companion paper in today's archive posting

Observational predictions for Thorne-Żytkow objects

Authors: R. Farmer, M. Renzo, Y. Götberg, E. Bellinger, S. Justham, S. E de Mink

Abstract: Thorne-$Ż$ytkow objects (T$Ż$O) are potential end products of the merger of a neutron star with a non-degenerate star. In this work, we have computed the first grid of evolutionary models of T$Ż$Os with the MESA stellar evolution code. With these models, we predict several observational properties of T$Ż$Os, including their surface temperatures and luminosities, pulsation periods, and nucleosynthe… ▽ More Thorne-$Ż$ytkow objects (T$Ż$O) are potential end products of the merger of a neutron star with a non-degenerate star. In this work, we have computed the first grid of evolutionary models of T$Ż$Os with the MESA stellar evolution code. With these models, we predict several observational properties of T$Ż$Os, including their surface temperatures and luminosities, pulsation periods, and nucleosynthetic products. We expand the range of possible T$Ż$O solutions to cover $3.45 \lesssim \log \left(T/K\right) \lesssim 3.65$ and $4.85 \lesssim \log \left(L/L_{\odot}\right) \lesssim 5.5$. Due to the much higher densities our T$Ż$Os reach compared to previous models, if T$Ż$Os form we expect them to be stable over a larger mass range than previously predicted, without exhibiting a gap in their mass distribution. Using the GYRE stellar pulsation code we show that T$Ż$Os should have fundamental pulsation periods of 1000--2000 days, and period ratios of $\approx$0.2--0.3. Models computed with a large 399 isotope fully-coupled nuclear network show a nucleosynthetic signal that is different to previously predicted. We propose a new nucleosynthetic signal to determine a star's status as a T$Ż$O: the isotopologues $^{44}\rm{Ti} \rm{O}_2$ and $^{44}\rm{Ti} \rm{O}$, which will have a shift in their spectral features as compared to stable titanium-containing molecules. We find that in the local Universe (~SMC metallicities and above) T$Ż$Os show little heavy metal enrichment, potentially explaining the difficulty in finding T$Ż$Os to-date. △ Less

Submitted 14 July, 2023; v1 submitted 12 May, 2023; originally announced May 2023.

Comments: 17 pages, 16 figures, 3 Tables, Accepted MNRAS, Zenodo data available https://doi.org/10.5281/zenodo.4534425

X-Shooting ULLYSES: massive stars at low metallicity. I. Project Description

Authors: Jorick S. Vink, A. Mehner, P. A. Crowther, A. Fullerton, M. Garcia, F. Martins, N. Morrell, L. M. Oskinova, N. St-Louis, A. ud-Doula, A. A. C. Sander, H. Sana, J. -C. Bouret, B. Kubatova, P. Marchant, L. P. Martins, A. Wofford, J. Th. van Loon, O. Grace Telford, Y. Gotberg, D. M. Bowman, C. Erba, V. M. Kalari, M. Abdul-Masih, T. Alkousa , et al. (56 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 1 June, 2023; v1 submitted 10 May, 2023; originally announced May 2023.

Comments: Accepted in A&A - 35 Pages, 12 Figures, 4 Tables, 2 Large Tables

Journal ref: A&A 675, A154 (2023)

Bringing Stellar Evolution & Feedback Together: Summary of proposals from the Lorentz Center Workshop, 2022

Authors: Sam Geen, Poojan Agrawal, Paul A. Crowther, B. W. Keller, Alex de Koter, Zsolt Keszthelyi, Freeke van de Voort, Ahmad A. Ali, Frank Backs, Lars Bonne, Vittoria Brugaletta, Annelotte Derkink, Sylvia Ekström, Yvonne A. Fichtner, Luca Grassitelli, Ylva Götberg, Erin R. Higgins, Eva Laplace, Kong You Liow, Marta Lorenzo, Anna F. McLeod, Georges Meynet, Megan Newsome, G. André Oliva, Varsha Ramachandran , et al. (12 additional authors not shown)

Abstract: Stars strongly impact their environment, and shape structures on all scales throughout the universe, in a process known as ``feedback''. Due to the complexity of both stellar evolution and the physics of larger astrophysical structures, there remain many unanswered questions about how feedback operates, and what we can learn about stars by studying their imprint on the wider universe. In this whit… ▽ More Stars strongly impact their environment, and shape structures on all scales throughout the universe, in a process known as ``feedback''. Due to the complexity of both stellar evolution and the physics of larger astrophysical structures, there remain many unanswered questions about how feedback operates, and what we can learn about stars by studying their imprint on the wider universe. In this white paper, we summarize discussions from the Lorentz Center meeting `Bringing Stellar Evolution and Feedback Together' in April 2022, and identify key areas where further dialogue can bring about radical changes in how we view the relationship between stars and the universe they live in. △ Less

Submitted 31 January, 2023; originally announced January 2023.

Comments: Accepted to the Publications of the Astronomical Society of the Pacific

A Panchromatic Study of Massive Stars in the Extremely Metal-Poor Local Group Dwarf Galaxy Leo A

Authors: Maude Gull, Daniel R. Weisz, Peter Senchyna, Nathan R. Sandford, Yumi Choi, Anna F. McLeod, Kareem El-Badry, Ylva Götberg, Karoline M. Gilbert, Martha Boyer, Julianne J. Dalcanton, Puragra GuhaThakurta, Steven Goldman, Paola Marigo, Kristen B. W. McQuinn, Giada Pastorelli, Daniel P. Stark, Evan Skillman, Yuan-sen Ting, Benjamin F. Williams

Abstract: We characterize massive stars (M>8 M_sun) in the nearby (D~0.8 Mpc) extremely metal-poor (Z~5% Z_sun) galaxy Leo A using Hubble Space Telescope ultra-violet (UV), optical, and near-infrared (NIR) imaging along with Keck/LRIS and MMT/Binospec optical spectroscopy for 18 main sequence OB stars. We find that: (a) 12 of our 18 stars show emission lines, despite not being associated with an H II region… ▽ More We characterize massive stars (M>8 M_sun) in the nearby (D~0.8 Mpc) extremely metal-poor (Z~5% Z_sun) galaxy Leo A using Hubble Space Telescope ultra-violet (UV), optical, and near-infrared (NIR) imaging along with Keck/LRIS and MMT/Binospec optical spectroscopy for 18 main sequence OB stars. We find that: (a) 12 of our 18 stars show emission lines, despite not being associated with an H II region, suggestive of stellar activity (e.g., mass loss, accretion, binary star interaction), which is consistent with previous predictions of enhanced activity at low metallicity; (b) 6 are Be stars, which are the first to be spectroscopically studied at such low metallicity -- these Be stars have unusual panchromatic SEDs; (c) for stars well-fit by the TLUSTY non-local thermodynamic equilibrium (non-LTE) models, the photometric and spectroscopic values of T_eff and log(g) agree to within ~0.01 dex and ~0.18 dex, respectively, indicating that NUV/optical/NIR imaging can be used to reliably characterize massive (M ~ 8-30 M_sun) main sequence star properties relative to optical spectroscopy; (d) the properties of the most massive stars in H II regions are consistent with constraints from previous nebular emission line studies; and (e) 13 stars with M>8 M_sun are >40 pc from a known star cluster or H II region. Our sample comprises ~50% of all known massive stars at Z < 10% Z_sun with derived stellar parameters, high-quality optical spectra, and panchromatic photometry. △ Less

Submitted 28 December, 2022; v1 submitted 25 November, 2022; originally announced November 2022.

Comments: 35 pages, 18 figures

Cool, Luminous, and Highly Variable Stars in the Magellanic Clouds. II: Spectroscopic and Environmental Analysis of Thorne-Żytkow Object and Super-AGB Star Candidates

Authors: Anna J. G. O'Grady, Maria R. Drout, B. M. Gaensler, C. S. Kochanek, Kathryn F. Neugent, Carolyn L. Doherty, Joshua S. Speagle, B. J. Shappee, Michael Rauch, Ylva Götberg, Bethany Ludwig, Todd A. Thompson

Abstract: In previous work we identified a population of 38 cool and luminous variable stars in the Magellanic Clouds and examined 11 in detail in order to classify them as either Thorne-Żytkow Objects (TŻOs, red supergiants with a neutron star cores) or super-AGB stars (the most massive stars that will not undergo core collapse). This population includes HV\,2112, a peculiar star previously considered in o… ▽ More In previous work we identified a population of 38 cool and luminous variable stars in the Magellanic Clouds and examined 11 in detail in order to classify them as either Thorne-Żytkow Objects (TŻOs, red supergiants with a neutron star cores) or super-AGB stars (the most massive stars that will not undergo core collapse). This population includes HV\,2112, a peculiar star previously considered in other works to be either a TŻO or high-mass AGB star. Here we continue this investigation, using the kinematic and radio environments and local star formation history of these stars to place constraints on the age of the progenitor systems and the presence of past supernovae. These stars are not associated with regions of recent star formation, and we find no evidence of past supernovae at their locations. Finally, we also assess the presence of heavy elements and lithium in their spectra compared to red supergiants. We find strong absorption in Li and s-process elements compared to RSGs in most of the sample, consistent with super-AGB nucleosynthesis, while HV\,2112 shows additional strong lines associated with TŻO nucleosynthesis. Coupled with our previous mass estimates, the results are consistent with the stars being massive (~4-6.5M$_{\odot}$) or super-AGB (~6.5-12M$_{\odot}$) stars in the thermally pulsing phase, providing crucial observations of the transition between low- and high-mass stellar populations. HV\,2112 is more ambiguous; it could either be a maximally massive sAGB star, or a TŻO if the minimum mass for stability extends down to <13 M$_\odot$. △ Less

Submitted 22 November, 2022; originally announced November 2022.

Comments: 32 pages, 10 figures, 7 tables, spectroscopic data available at https://zenodo.org/record/7058608, accepted to The Astrophysical Journal

Spin-down and reduced mass loss in early-type stars with large-scale magnetic fields

Authors: Z. Keszthelyi, A. de Koter, Y. Götberg, G. Meynet, S. A. Brands, V. Petit, M. Carrington, A. David-Uraz, S. T. Geen, C. Georgy, R. Hirschi, J. Puls, K. J. Ramalatswa, M. E. Shultz, A. ud-Doula

Abstract: Magnetism can greatly impact the evolution of stars. In some stars with OBA spectral types there is direct evidence via the Zeeman effect for stable, large-scale magnetospheres, which lead to the spin-down of the stellar surface and reduced mass loss. So far, a comprehensive grid of stellar structure and evolution models accounting for these effects was lacking. For this reason, we computed and st… ▽ More Magnetism can greatly impact the evolution of stars. In some stars with OBA spectral types there is direct evidence via the Zeeman effect for stable, large-scale magnetospheres, which lead to the spin-down of the stellar surface and reduced mass loss. So far, a comprehensive grid of stellar structure and evolution models accounting for these effects was lacking. For this reason, we computed and studied models with two magnetic braking and two chemical mixing schemes in three metallicity environments with the MESA software instrument. We find notable differences between the subgrids, which affects the model predictions and thus the detailed characterisation of stars. We are able to quantify the impact of magnetic fields in terms of preventing quasi-chemically homogeneous evolution and producing slowly-rotating, nitrogen-enriched ("Group 2") stars. Our model grid is fully open access and open source. △ Less

Submitted 13 November, 2022; originally announced November 2022.

Comments: to appear in "Winds of Stars and Exoplanets" Proceedings IAUS 370, 2022, eds.: A.A. Vidotto, L. Fossati, J.S. Vink

The effects of surface fossil magnetic fields on massive star evolution: IV. Grids of models at Solar, LMC, and SMC metallicities

Authors: Z. Keszthelyi, A. de Koter, Y. Götberg, G. Meynet, S. A. Brands, V. Petit, M. Carrington, A. David-Uraz, S. T. Geen, C. Georgy, R. Hirschi, J. Puls, K. J. Ramalatswa, M. E. Shultz, A. ud-Doula

Abstract: Magnetic fields can drastically change predictions of evolutionary models of massive stars via mass-loss quenching, magnetic braking, and efficient angular momentum transport, which we aim to quantify in this work. We use the MESA software instrument to compute an extensive main-sequence grid of stellar structure and evolution models, as well as isochrones, accounting for the effects attributed to… ▽ More Magnetic fields can drastically change predictions of evolutionary models of massive stars via mass-loss quenching, magnetic braking, and efficient angular momentum transport, which we aim to quantify in this work. We use the MESA software instrument to compute an extensive main-sequence grid of stellar structure and evolution models, as well as isochrones, accounting for the effects attributed to a surface fossil magnetic field. The grid is densely populated in initial mass (3-60 M$_\odot$), surface equatorial magnetic field strength (0-50 kG), and metallicity (representative of the Solar neighbourhood and the Magellanic Clouds). We use two magnetic braking and two chemical mixing schemes and compare the model predictions for slowly-rotating, nitrogen-enriched ("Group 2") stars with observations in the Large Magellanic Cloud. We quantify a range of initial field strengths that allow for producing Group 2 stars and find that typical values (up to a few kG) lead to solutions. Between the subgrids, we find notable departures in surface abundances and evolutionary paths. In our magnetic models, chemical mixing is always less efficient compared to non-magnetic models due to the rapid spin-down. We identify that quasi-chemically homogeneous main sequence evolution by efficient mixing could be prevented by fossil magnetic fields. We recommend comparing this grid of evolutionary models with spectropolarimetric and spectroscopic observations with the goals of i) revisiting the derived stellar parameters of known magnetic stars, and ii) observationally constraining the uncertain magnetic braking and chemical mixing schemes. △ Less

Submitted 28 October, 2022; v1 submitted 13 September, 2022; originally announced September 2022.

Comments: Accepted for publication in MNRAS. A full reproduction package is available on Zenodo at https://doi.org/10.5281/zenodo.7069766

Science with the Ultraviolet Explorer (UVEX)

Authors: S. R. Kulkarni, Fiona A. Harrison, Brian W. Grefenstette, Hannah P. Earnshaw, Igor Andreoni, Danielle A. Berg, Joshua S. Bloom, S. Bradley Cenko, Ryan Chornock, Jessie L. Christiansen, Michael W. Coughlin, Alexander Wuollet Criswell, Behnam Darvish, Kaustav K. Das, Kishalay De, Luc Dessart, Don Dixon, Bas Dorsman, Kareem El-Badry, Christopher Evans, K. E. Saavik Ford, Christoffer Fremling, Boris T. Gansicke, Suvi Gezari, Y. Goetberg , et al. (31 additional authors not shown)

Abstract: UVEX is a proposed medium class Explorer mission designed to provide crucial missing capabilities that will address objectives central to a broad range of modern astrophysics. The UVEX design has two co-aligned wide-field imagers operating in the FUV and NUV and a powerful broadband medium resolution spectrometer. In its two-year baseline mission, UVEX will perform a multi-cadence synoptic all-sky… ▽ More UVEX is a proposed medium class Explorer mission designed to provide crucial missing capabilities that will address objectives central to a broad range of modern astrophysics. The UVEX design has two co-aligned wide-field imagers operating in the FUV and NUV and a powerful broadband medium resolution spectrometer. In its two-year baseline mission, UVEX will perform a multi-cadence synoptic all-sky survey 50/100 times deeper than GALEX in the NUV/FUV, cadenced surveys of the Large and Small Magellanic Clouds, rapid target of opportunity followup, as well as spectroscopic followup of samples of stars and galaxies. The science program is built around three pillars. First, UVEX will explore the low-mass, low-metallicity galaxy frontier through imaging and spectroscopic surveys that will probe key aspects of the evolution of galaxies by understanding how star formation and stellar evolution at low metallicities affect the growth and evolution of low-metallicity, low-mass galaxies in the local universe. Such galaxies contain half the mass in the local universe, and are analogs for the first galaxies, but observed at distances that make them accessible to detailed study. Second, UVEX will explore the dynamic universe through time-domain surveys and prompt spectroscopic followup capability will probe the environments, energetics, and emission processes in the early aftermaths of gravitational wave-discovered compact object mergers, discover hot, fast UV transients, and diagnose the early stages of stellar explosions. Finally, UVEX will become a key community resource by leaving a large all-sky legacy data set, enabling a wide range of scientific studies and filling a gap in the new generation of wide-field, sensitive optical and infrared surveys provided by the Rubin, Euclid, and Roman observatories. This paper discusses the scientific potential of UVEX, and the broad scientific program. △ Less

Submitted 17 January, 2023; v1 submitted 30 November, 2021; originally announced November 2021.

Comments: 69 pages, 43 figures

Pre-Explosion Properties of Helium Star Donors to Thermonuclear Supernovae

Authors: Tin Long Sunny Wong, Josiah Schwab, Ylva Götberg

Abstract: Helium star - carbon-oxygen white dwarf (CO WD) binaries are potential single-degenerate progenitor systems of thermonuclear supernovae. Revisiting a set of binary evolution calculations using the stellar evolution code $\texttt{MESA}$, we refine our previous predictions about which systems can lead to a thermonuclear supernova and then characterize the properties of the helium star donor at the t… ▽ More Helium star - carbon-oxygen white dwarf (CO WD) binaries are potential single-degenerate progenitor systems of thermonuclear supernovae. Revisiting a set of binary evolution calculations using the stellar evolution code $\texttt{MESA}$, we refine our previous predictions about which systems can lead to a thermonuclear supernova and then characterize the properties of the helium star donor at the time of explosion. We convert these model properties to NUV/optical magnitudes assuming a blackbody spectrum and support this approach using a matched stellar atmosphere model. These models will be valuable to compare with pre-explosion imaging for future supernovae, though we emphasize the observational difficulty of detecting extremely blue companions. The pre-explosion source detected in association with SN 2012Z has been interpreted as a helium star binary containing an initially ultra-massive WD in a multi-day orbit. However, extending our binary models to initial CO WD masses of up to $1.2\,M_{\odot}$, we find that these systems undergo off-center carbon ignitions and thus are not expected to produce thermonuclear supernovae. This tension suggests that, if SN 2012Z is associated with a helium star - WD binary, then the pre-explosion optical light from the system must be significantly modified by the binary environment and/or the WD does not have a carbon-rich interior composition. △ Less

Submitted 29 September, 2021; originally announced September 2021.

Comments: 26 pages, 12 figures, 3 tables; Accepted to ApJ

Evolution of accretor stars in massive binaries: broader implications from modeling ζ Ophiuchi

Authors: M. Renzo, Y. Gotberg

Abstract: Most massive stars are born in binaries close enough for mass transfer episodes. These modify the appearance, structure, and future evolution of both stars. We compute the evolution of a 100-day period binary consisting initially of a 25 M star and a 17 M star, which experiences stable mass transfer. We focus on the impact of mass accretion on the surface composition, internal rotation, and struct… ▽ More Most massive stars are born in binaries close enough for mass transfer episodes. These modify the appearance, structure, and future evolution of both stars. We compute the evolution of a 100-day period binary consisting initially of a 25 M star and a 17 M star, which experiences stable mass transfer. We focus on the impact of mass accretion on the surface composition, internal rotation, and structure of the accretor. To anchor our models, we show that our accretor broadly reproduces the properties of ζ Ophiuchi, which has long been proposed to have accreted mass before being ejected as a runaway star when the companion exploded. We compare our accretor to models of single rotating stars and find that the later and stronger spin-up provided by mass accretion produces significant differences. Specifically, the core of the accretor retains higher spin at the end of the main sequence, and a convective layer develops that changes its density profile. Moreover, the surface of the accretor star is polluted by CNO-processed material donated by the companion. Our models show effects of mass accretion in binaries that are not captured in single rotating stellar models. This possibly impacts the further evolution (either in a binary or as single stars), the final collapse, and the resulting spin of the compact object. △ Less

Submitted 23 September, 2021; v1 submitted 22 July, 2021; originally announced July 2021.

Comments: 19 pages+appendix, code available at https://github.com/mathren90/zeta_oph and results at https://zenodo.org/record/4701565

The young massive SMC cluster NGC 330 seen by MUSE II. Multiplicity properties of the massive-star population

Authors: J. Bodensteiner, H. Sana, C. Wang, N. Langer, L. Mahy, G. Banyard, A. de Koter, S. E. de Mink, C. J. Evans, Y. Götberg, L. R. Patrick, F. R. N. Schneider, F. Tramper

Abstract: Observations of massive stars in young open clusters (< ~8 Myr) have shown that a majority of them are in binary systems, most of which will interact during their life. Populations of massive stars older than ~20 Myr allow us to probe the outcome of such interactions after many systems have experienced mass and angular momentum transfer. Using multi-epoch integral-field spectroscopy, we investigat… ▽ More Observations of massive stars in young open clusters (< ~8 Myr) have shown that a majority of them are in binary systems, most of which will interact during their life. Populations of massive stars older than ~20 Myr allow us to probe the outcome of such interactions after many systems have experienced mass and angular momentum transfer. Using multi-epoch integral-field spectroscopy, we investigate the multiplicity properties of the massive-star population in NGC 330 (~40 Myr) in the Small Magellanic Cloud to search for imprints of stellar evolution on the multiplicity properties. From six epochs of VLT/MUSE observations supported by adaptive optics we extract spectra and measure radial velocities for stars brighter than F814W = 19. We identify single-lined spectroscopic binaries through significant RV variability as well as double-lined spectroscopic binaries, and quantify the observational biases for binary detection. The observed spectroscopic binary fraction is 13.2+/-2.0 %. Considering period and mass ratio ranges from log(P)=0.15-3.5, and q = 0.1-1.0, and a representative set of orbital parameter distributions, we find a bias-corrected close binary fraction of 34 +8 -7 %. This seems to decline for the fainter stars, which indicates either that the close binary fraction drops in the B-type domain, or that the period distribution becomes more heavily weighted towards longer orbital periods. Both fractions vary strongly in different regions of the color-magnitude diagram which probably reveals the imprint of the binary history of different groups of stars. We provide the first homogeneous RV study of a large sample of B-type stars at a low metallicity. The overall bias-corrected close binary fraction of B stars in NGC 330 is lower than the one reported for younger Galactic and LMC clusters. More data are needed to establish whether this result from an age or a metallicty effect. △ Less

Submitted 27 April, 2021; originally announced April 2021.

Comments: 19 page (incl. appendix), 10 figures, 4 tables, accepted for publication in A&A

Journal ref: A&A 652, A70 (2021)

Binary-Stripped Stars as Core-Collapse Supernovae Progenitors

Authors: David Vartanyan, Eva Laplace, Mathieu Renzo, Ylva Götberg, Adam Burrows, Selma E. de Mink

Abstract: Most massive stars experience binary interactions in their lifetimes that can alter both the surface and core structure of the stripped star with significant effects on their ultimate fate as core-collapse supernovae. However, core-collapse supernovae simulations to date have focused almost exclusively on the evolution of single stars. We present a systematic simulation study of single and binary-… ▽ More Most massive stars experience binary interactions in their lifetimes that can alter both the surface and core structure of the stripped star with significant effects on their ultimate fate as core-collapse supernovae. However, core-collapse supernovae simulations to date have focused almost exclusively on the evolution of single stars. We present a systematic simulation study of single and binary-stripped stars with the same initial mass as candidates for core-collapse supernovae (11 - 21 M$_{\odot}$). Generally, we find that binary-stripped stars core tend to be less compact, with a more prominent, deeper silicon/oxygen interface, and explode preferentially to the corresponding single stars of the same initial mass. Such a dichotomy of behavior between these two modes of evolution would have important implications for supernovae statistics, including the final neutron star masses, explosion energies, and nucleosynthetic yields. Binary-stripped remnants are also well poised to populate the possible mass gap between the heaviest neutron stars and the lightest black holes. Our work presents an improvement along two fronts, as we self-consistently account for the pre-collapse stellar evolution and the subsequent explosion outcome. Even so, our results emphasize the need for more detailed stellar evolutionary models to capture the sensitive nature of explosion outcome. △ Less

Submitted 7 April, 2021; originally announced April 2021.

Comments: submitted to ApJ

The Detection and Characterization of Be+sdO Binaries from HST/STIS FUV Spectroscopy

Authors: Luqian Wang, Douglas R. Gies, Geraldine J. Peters, Ylva Götberg, S. Drew Chojnowski, Kathryn V. Lester, Steve B. Howell

Abstract: The B-emission line stars are rapid rotators that were probably spun up by mass and angular momentum accretion through mass transfer in an interacting binary. Mass transfer will strip the donor star of its envelope to create a small and hot subdwarf remnant. Here we report on Hubble Space Telescope/STIS far-ultraviolet spectroscopy of a sample of Be stars that reveals the presence of the hot sdO c… ▽ More The B-emission line stars are rapid rotators that were probably spun up by mass and angular momentum accretion through mass transfer in an interacting binary. Mass transfer will strip the donor star of its envelope to create a small and hot subdwarf remnant. Here we report on Hubble Space Telescope/STIS far-ultraviolet spectroscopy of a sample of Be stars that reveals the presence of the hot sdO companion through the calculation of cross-correlation functions of the observed and model spectra. We clearly detect the spectral signature of the sdO star in 10 of the 13 stars in the sample, and the spectral signals indicate that the sdO stars are hot, relatively faint, and slowly rotating as predicted by models. A comparison of their temperatures and radii with evolutionary tracks indicates that the sdO stars occupy the relatively long-lived, He-core burning stage. Only one of the ten detections was a known binary prior to this investigation, which emphasizes the difficulty of finding such Be+sdO binaries through optical spectroscopy. However, these results and others indicate that many Be stars probably host hot subdwarf companions. △ Less

Submitted 25 March, 2021; originally announced March 2021.

Comments: 36 pages, 38 figures, AJ Accepted

Spectral Signatures of Population III and Envelope-stripped Stars in Galaxies at the Epoch of Reionization

Authors: Elizabeth Berzin, Amy Secunda, Renyue Cen, Alexander Menegas, Ylva Götberg

Abstract: While most simulations of the epoch of reionization have focused on single-stellar populations in star-forming dwarf galaxies, products of binary evolution are expected to significantly contribute to emissions of hydrogen-ionizing photons. Among these products are stripped stars (or helium stars), which have their envelopes stripped from interactions with binary companions, leaving an exposed heli… ▽ More While most simulations of the epoch of reionization have focused on single-stellar populations in star-forming dwarf galaxies, products of binary evolution are expected to significantly contribute to emissions of hydrogen-ionizing photons. Among these products are stripped stars (or helium stars), which have their envelopes stripped from interactions with binary companions, leaving an exposed helium core. Previous work has suggested these stripped stars can dominate the LyC photon output of high-redshift low luminosity galaxies. Other sources of hard radiation in the early universe include zero-metallicity Population III stars, which may have similar SED properties to galaxies with radiation dominated by stripped star emissions. Here, we use two metrics (the power-law exponent over wavelength intervals 240-500 Å, 600-900 Å, and 1200-2000 Å, and the ratio of total luminosity in FUV wavelengths to LyC wavelengths) to compare the SEDs of simulated galaxies with only single-stellar evolution, galaxies containing stripped stars, and galaxies containing Population III stars, with four different IMFs. We find that stripped stars significantly alter the SEDs in the LyC range of galaxies at the epoch of reionization. SEDs in galaxies with stripped stars present have lower power-law indices in the LyC range and lower FUV to LyC luminosity ratios. These differences in SEDs are present at all considered luminosities ($M_{UV} > -15$, AB system), and are most pronounced for lower luminosity galaxies. We also find that SEDs of galaxies with stripped stars and Pop III stars are distinct from each other for all tested IMFs. △ Less

Submitted 10 August, 2021; v1 submitted 16 February, 2021; originally announced February 2021.

Comments: 16 pages, 9 figures. Accepted to ApJ

Different to the core: the pre-supernova structures of massive single and binary-stripped stars

Authors: E. Laplace, S. Justham, M. Renzo, Y. Götberg, R. Farmer, D. Vartanyan, S. E. de Mink

Abstract: The majority of massive stars live in binary or multiple systems and will interact during their lifetimes, which helps to explain the observed diversity of core-collapse supernovae. Donor stars in binary systems can lose most of their hydrogen-rich envelopes through mass transfer, which not only affects the surface properties, but also the core structure. However, most calculations of the core-col… ▽ More The majority of massive stars live in binary or multiple systems and will interact during their lifetimes, which helps to explain the observed diversity of core-collapse supernovae. Donor stars in binary systems can lose most of their hydrogen-rich envelopes through mass transfer, which not only affects the surface properties, but also the core structure. However, most calculations of the core-collapse properties of massive stars rely on single-star models. We present a systematic study of the difference between the pre-supernova structures of single stars and stars of the same initial mass (11 - 21\Msun) that have been stripped due to stable post-main sequence mass transfer at solar metallicity. We present the pre-supernova core composition with novel diagrams that give an intuitive representation of the isotope distribution. As shown in previous studies, at the edge of the carbon-oxygen core, the binary-stripped star models contain an extended gradient of carbon, oxygen, and neon. This layer originates from the receding of the convective helium core during core helium burning in binary-stripped stars, which does not occur in single-star models. We find that this same evolutionary phase leads to systematic differences in the final density and nuclear energy generation profiles. Binary-stripped star models have systematically higher total masses of carbon at the moment of core collapse compared to single star models, which likely results in systematically different supernova yields. In about half of our models, the silicon-burning and oxygen-rich layers merge after core silicon burning. We discuss the implications of our findings for the explodability, supernova observations, and nucleosynthesis from these stars. Our models will be publicly available and can be readily used as input for supernova simulations. [Abridged] △ Less

Submitted 27 October, 2021; v1 submitted 9 February, 2021; originally announced February 2021.

Comments: Accepted for publication on October 5th, 2021. Reproduced with permission from Astronomy & Astrophysics, \c{opyright} ESO. The pre-supernova models are available online at https://doi.org/10.5281/zenodo.4506803

Journal ref: A&A 656, A58 (2021)

Successful Common Envelope Ejection and Binary Neutron Star Formation in 3D Hydrodynamics

Authors: Jamie A. P. Law-Smith, Rosa Wallace Everson, Enrico Ramirez-Ruiz, Selma E. de Mink, Lieke A. C. van Son, Ylva Götberg, Stefan Zellmann, Alejandro Vigna-Gómez, Mathieu Renzo, Samantha Wu, Sophie L. Schrøder, Ryan J. Foley, Tenley Hutchinson-Smith

Abstract: A binary neutron star merger has been observed in a multi-messenger detection of gravitational wave (GW) and electromagnetic (EM) radiation. Binary neutron stars that merge within a Hubble time, as well as many other compact binaries, are expected to form via common envelope evolution. Yet five decades of research on common envelope evolution have not yet resulted in a satisfactory understanding o… ▽ More A binary neutron star merger has been observed in a multi-messenger detection of gravitational wave (GW) and electromagnetic (EM) radiation. Binary neutron stars that merge within a Hubble time, as well as many other compact binaries, are expected to form via common envelope evolution. Yet five decades of research on common envelope evolution have not yet resulted in a satisfactory understanding of the multi-spatial multi-timescale evolution for the systems that lead to compact binaries. In this paper, we report on the first successful simulations of common envelope ejection leading to binary neutron star formation in 3D hydrodynamics. We simulate the dynamical inspiral phase of the interaction between a 12$M_\odot$ red supergiant and a 1.4$M_\odot$ neutron star for different initial separations and initial conditions. For all of our simulations, we find complete envelope ejection and final orbital separations of $a_{\rm f} \approx 1.3$-$5.1 R_\odot$ depending on the simulation and criterion, leading to binary neutron stars that can merge within a Hubble time. We find $α_{\rm CE}$-equivalent efficiencies of $\approx 0.1$-$2.7$ depending on the simulation and criterion, but this may be specific for these extended progenitors. We fully resolve the core of the star to $\lesssim 0.005 R_\odot$ and our 3D hydrodynamics simulations are informed by an adjusted 1D analytic energy formalism and a 2D kinematics study in order to overcome the prohibitive computational cost of simulating these systems. The framework we develop in this paper can be used to simulate a wide variety of interactions between stars, from stellar mergers to common envelope episodes leading to GW sources. △ Less

Submitted 22 July, 2022; v1 submitted 12 November, 2020; originally announced November 2020.

Comments: 28 pages, 16 figures. v2: ran some additional simulations to address comments from the community. Included new criterion for calculating a_f and alpha_CE, resulting in new range of values. Added Table 1, updated Fig 3, added Fig 4, updated Fig 12, updated Fig 13, added Fig 14, added Fig 15, added Fig 16. Now submitted to ApJ

Delayed Photons from Binary Evolution Help Reionize the Universe

Authors: Amy Secunda, Renyue Cen, Taysun Kimm, Ylva Gotberg, Selma E. de Mink

Abstract: High-resolution numerical simulations including feedback and aimed at calculating the escape fraction (fesc) of hydrogen-ionizing photons often assume stellar radiation based on single-stellar population synthesis models. However, strong evidence suggests the binary fraction of massive stars is 70%. Moreover, simulations so far yield values of fesc falling only on the lower end of the roughly 10-2… ▽ More High-resolution numerical simulations including feedback and aimed at calculating the escape fraction (fesc) of hydrogen-ionizing photons often assume stellar radiation based on single-stellar population synthesis models. However, strong evidence suggests the binary fraction of massive stars is 70%. Moreover, simulations so far yield values of fesc falling only on the lower end of the roughly 10-20% range, the amount presumed necessary to reionize the Universe. Analyzing a high-resolution (4 pc) cosmological radiation hydrodynamic simulation we study how fesc changes when we include two different products of binary stellar evolution - stars stripped of their hydrogen envelopes and massive blue stragglers. Both produce significant amounts of ionizing photons 10-200 Myr after each starburst. We find the relative importance of these photons are amplified with respect to escaped ionizing photons, because peaks in star formation rates (SFRs) and fesc are often out of phase by this 10-200 Myr. Additionally, low mass, bursty galaxies emit Lyman continuum radiation primarily from binary products when SFRs are low. Observations of these galaxies by the James Webb Space Telescope could provide crucial information on the evolution of binary stars as a function of redshift. Overall, including stripped stars and massive blue stragglers increases our photon-weighted mean escape fraction by around 13% and 10%, respectively, resulting in a mean fesc of 17%. Our results emphasize that using updated stellar population synthesis models with binary stellar evolution provides a more sound physical basis for stellar reionization. △ Less

Submitted 29 July, 2020; originally announced July 2020.

Comments: 20 pages, 7 figures, submitted to ApJ

Stars stripped in binaries -- the living gravitational wave sources

Authors: Y. Götberg, V. Korol, A. Lamberts, T. Kupfer, K. Breivik, B. Ludwig, M. R. Drout

Abstract: Binary interaction can cause stellar envelopes to be stripped, which significantly reduces the radius of the star. The orbit of a binary composed of a stripped star and a compact object can therefore be so tight that the gravitational radiation the system produces reaches frequencies accessible to the Laser Interferometer Space Antenna (LISA). Two such stripped stars in tight orbits with white dwa… ▽ More Binary interaction can cause stellar envelopes to be stripped, which significantly reduces the radius of the star. The orbit of a binary composed of a stripped star and a compact object can therefore be so tight that the gravitational radiation the system produces reaches frequencies accessible to the Laser Interferometer Space Antenna (LISA). Two such stripped stars in tight orbits with white dwarfs are known so far (ZTF J2130+4420 and CD-30 11223), but many more are expected to exist. These binaries provide important constraints for binary evolution models and may be used as LISA verification sources. We develop a Monte Carlo code that uses detailed evolutionary models to simulate the Galactic population of stripped stars in tight orbits with either neutron star or white dwarf companions. We predict 0-100 stripped star + white dwarf binaries and 0-4 stripped star + neutron star binaries with SNR>5 after 10 years of observations with LISA. More than 90% of these binaries are expected to show large radial velocity shifts of $\gtrsim 200$ km/s, which are spectroscopically detectable. Photometric variability due to tidal deformation of the stripped star is also expected and has been observed in ZTF J2130+4420 and CD-30 11223. In addition, the stripped star + neutron star binaries are expected to be X-ray bright with $L_X \gtrsim 10^{33} - 10^{36}$ erg/s. Our results show that stripped star binaries are promising multi-messenger sources for the upcoming electromagnetic and gravitational wave facilities. △ Less

Submitted 12 June, 2020; originally announced June 2020.

Comments: Submitted to ApJ

The expansion of stripped-envelope stars: consequences for supernovae and gravitational-wave progenitors

Authors: E. Laplace, Y. Götberg, S. E. de Mink, S. Justham, R. Farmer

Abstract: Massive binaries that merge as compact objects are the progenitors of gravitational-wave sources. Most of these binaries experience one or more phases of mass transfer, during which one of the stars loses part or all of its outer envelope and becomes a stripped-envelope star. The evolution of the size of these stripped stars is crucial in determining whether they experience further interactions an… ▽ More Massive binaries that merge as compact objects are the progenitors of gravitational-wave sources. Most of these binaries experience one or more phases of mass transfer, during which one of the stars loses part or all of its outer envelope and becomes a stripped-envelope star. The evolution of the size of these stripped stars is crucial in determining whether they experience further interactions and their final fate. We present new calculations of stripped-envelope stars based on binary evolution models computed with MESA. We use these to investigate their radius evolution as a function of mass and metallicity. We further discuss their pre-supernova observable characteristics and potential consequences of their evolution on the properties of supernovae from stripped stars. At high metallicity we find that practically all of the hydrogen-rich envelope is removed, in agreement with earlier findings. Only progenitors with initial masses below 10\Msun expand to large radii (up to 100\Rsun), while more massive progenitors stay compact. At low metallicity, a substantial amount of hydrogen remains and the progenitors can, in principle, expand to giant sizes (> 400\Rsun), for all masses we consider. This implies that they can fill their Roche lobe anew. We show that the prescriptions commonly used in population synthesis models underestimate the stellar radii by up to two orders of magnitude. We expect that this has consequences for the predictions for gravitational-wave sources from double neutron star mergers, in particular for their metallicity dependence. △ Less

Submitted 2 March, 2020; originally announced March 2020.

Comments: Main text 17 pages, 7 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 637, A6 (2020)

Sensitivity of the lower-edge of the pair instability black hole mass gap to the treatment of time dependent convection

Authors: M. Renzo, R. J. Farmer, S. Justham, S. E. de Mink, Y. Götberg, P. Marchant

Abstract: Gravitational-wave detections are now probing the black hole (BH) mass distribution, including the predicted pair-instability mass gap. These data require robust quantitative predictions, which are challenging to obtain. The most massive BH progenitors experience episodic mass ejections on timescales shorter than the convective turn-over timescale. This invalidates the steady-state assumption on w… ▽ More Gravitational-wave detections are now probing the black hole (BH) mass distribution, including the predicted pair-instability mass gap. These data require robust quantitative predictions, which are challenging to obtain. The most massive BH progenitors experience episodic mass ejections on timescales shorter than the convective turn-over timescale. This invalidates the steady-state assumption on which the classic mixing-length theory relies. We compare the final BH masses computed with two different versions of the stellar evolutionary code \texttt{MESA}: (i) using the default implementation of \cite{paxton:18} and (ii) solving an additional equation accounting for the timescale for convective deceleration. In the second grid, where stronger convection develops during the pulses and carries part of the energy, we find weaker pulses. This leads to lower amounts of mass being ejected and thus higher final BH masses of up to $\sim$\,$5\,M_\odot$. The differences are much smaller for the progenitors which determine the maximum mass of BHs below the gap. This prediction is robust at $M_{\rm BH, max}\simeq 48\,M_\odot$, at least within the idealized context of this study. This is an encouraging indication that current models are robust enough for comparison with the present-day gravitational-wave detections. However, the large differences between individual models emphasize the importance of improving the treatment of convection in stellar models, especially in the light of the data anticipated from the third generation of gravitational wave detectors. △ Less

Submitted 19 February, 2020; originally announced February 2020.

Comments: 7 pages + 1 appendix, accepted in MNRAS, online results at https://zenodo.org/record/3406320

Predictions for the hydrogen-free ejecta of pulsational pair-instability supernovae

Authors: M. Renzo, R. Farmer, S. Justham, Y. Götberg, S. E. de Mink, E. Zapartas, P. Marchant, N. Smith

Abstract: Present time-domain astronomy efforts will unveil a variety of rare transients. We focus here on pulsational pair-instability evolution, which can result in signatures observable with electromagnetic and gravitational waves. We simulate grids of bare helium stars to characterize the resulting black hole (BH) masses and ejecta composition, velocity, and thermal state. The stars do not react "elasti… ▽ More Present time-domain astronomy efforts will unveil a variety of rare transients. We focus here on pulsational pair-instability evolution, which can result in signatures observable with electromagnetic and gravitational waves. We simulate grids of bare helium stars to characterize the resulting black hole (BH) masses and ejecta composition, velocity, and thermal state. The stars do not react "elastically" to the thermonuclear explosion: there is not a one-to-one correspondence between pair-instability driven ignition and mass ejections, causing ambiguity in what is an observable pulse. In agreement with previous studies, we find that for carbon-oxygen core masses 28Msun< M_CO<30.5Msun the explosions are not strong enough to affect the surface. With increasing mass, they first cause large radial expansion (30.5Msun<M_CO<31.4Msun), and finally, also mass ejection episodes (M_CO>31.4Msun). The lowest mass to be fully disrupted in a pair-instability supernova is M_CO=57Msun. Models with M_CO>121Msun reach the photodisintegration regime, resulting in BHs with M_BH>125Msun. If the pulsating models produce BHs via (weak) explosions, the previously-ejected material might be hit by the blast wave. We characterize the H-free circumstellar material from the pulsational pair-instability of helium cores assuming simply that the ejecta maintain a constant velocity after ejection. Our models produce He-rich ejecta with mass 10^{-3}Msun<M_CSM<40Msun. These ejecta are typically launched at a few thousand \kms and reach distances of ~10^{12}-10^{15} cm before core-collapse. The delays between mass ejection events and the final collapse span a wide and mass-dependent range (from sub-hour to 10^4 years), and the shells ejected can also collide with each other. The range of properties we find suggests a possible connection with (some) type Ibn supernovae. △ Less

Submitted 16 June, 2020; v1 submitted 12 February, 2020; originally announced February 2020.

Comments: accepted version

Journal ref: A&A 640, A56 (2020)

The Habitable Exoplanet Observatory (HabEx) Mission Concept Study Final Report

Authors: B. Scott Gaudi, Sara Seager, Bertrand Mennesson, Alina Kiessling, Keith Warfield, Kerri Cahoy, John T. Clarke, Shawn Domagal-Goldman, Lee Feinberg, Olivier Guyon, Jeremy Kasdin, Dimitri Mawet, Peter Plavchan, Tyler Robinson, Leslie Rogers, Paul Scowen, Rachel Somerville, Karl Stapelfeldt, Christopher Stark, Daniel Stern, Margaret Turnbull, Rashied Amini, Gary Kuan, Stefan Martin, Rhonda Morgan , et al. (161 additional authors not shown)

Abstract: The Habitable Exoplanet Observatory, or HabEx, has been designed to be the Great Observatory of the 2030s. For the first time in human history, technologies have matured sufficiently to enable an affordable space-based telescope mission capable of discovering and characterizing Earthlike planets orbiting nearby bright sunlike stars in order to search for signs of habitability and biosignatures. Su… ▽ More The Habitable Exoplanet Observatory, or HabEx, has been designed to be the Great Observatory of the 2030s. For the first time in human history, technologies have matured sufficiently to enable an affordable space-based telescope mission capable of discovering and characterizing Earthlike planets orbiting nearby bright sunlike stars in order to search for signs of habitability and biosignatures. Such a mission can also be equipped with instrumentation that will enable broad and exciting general astrophysics and planetary science not possible from current or planned facilities. HabEx is a space telescope with unique imaging and multi-object spectroscopic capabilities at wavelengths ranging from ultraviolet (UV) to near-IR. These capabilities allow for a broad suite of compelling science that cuts across the entire NASA astrophysics portfolio. HabEx has three primary science goals: (1) Seek out nearby worlds and explore their habitability; (2) Map out nearby planetary systems and understand the diversity of the worlds they contain; (3) Enable new explorations of astrophysical systems from our own solar system to external galaxies by extending our reach in the UV through near-IR. This Great Observatory science will be selected through a competed GO program, and will account for about 50% of the HabEx primary mission. The preferred HabEx architecture is a 4m, monolithic, off-axis telescope that is diffraction-limited at 0.4 microns and is in an L2 orbit. HabEx employs two starlight suppression systems: a coronagraph and a starshade, each with their own dedicated instrument. △ Less

Submitted 26 January, 2020; v1 submitted 18 January, 2020; originally announced January 2020.

Comments: Full report: 498 pages. Executive Summary: 14 pages. More information about HabEx can be found here: https://www.jpl.nasa.gov/habex/

The young massive SMC cluster NGC 330 seen by MUSE. I. Observations and stellar content

Authors: J. Bodensteiner, H. Sana, L. Mahy, L. R. Patrick, A. de Koter, S. E. de Mink, C. J. Evans, Y. Götberg, N. Langer, D. J. Lennon, F. R. N. Schneider, F. Tramper

Abstract: 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 signif… ▽ More 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. △ Less

Submitted 8 November, 2019; originally announced November 2019.

Comments: 19 pages (incl. appendix), 21 figures, 3 tables, accepted for publication in A&A

Journal ref: A&A 634, A51 (2020)

The contribution from stars stripped in binaries to cosmic reionization of hydrogen and helium

Authors: Y. Gotberg, S. E. de Mink, M. McQuinn, E. Zapartas, J. H. Groh, C. Norman

Abstract: 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… ▽ More 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. △ Less

Submitted 1 November, 2019; originally announced November 2019.

Comments: Under review in A&A, suggestions welcome

Journal ref: A&A 634, A134 (2020)

The impact of stars stripped in binaries on the integrated spectra of stellar populations

Authors: Y. Götberg, S. E. de Mink, J. H. Groh, C. Leitherer, C. Norman

Abstract: 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 e… ▽ More 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. △ Less

Submitted 16 August, 2019; originally announced August 2019.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 629, A134 (2019)

The diverse lives of progenitors of hydrogen-rich core-collapse supernovae: the role of binary interaction

Authors: Emmanouil Zapartas, Selma E. de Mink, Stephen Justham, Nathan Smith, Alex de Koter, Mathieu Renzo, Iair Arcavi, Rob Farmer, Ylva Götberg, Silvia Toonen

Abstract: Hydrogen-rich supernovae, known as Type II (SNe II), are the most common class of explosions observed following the collapse of the core of massive stars. We use analytical estimates and population synthesis simulations to assess the fraction of SNe II progenitors that are expected to have exchanged mass with a companion prior to explosion. We estimate that 1/3 to 1/2 of SN II progenitors have a h… ▽ More Hydrogen-rich supernovae, known as Type II (SNe II), are the most common class of explosions observed following the collapse of the core of massive stars. We use analytical estimates and population synthesis simulations to assess the fraction of SNe II progenitors that are expected to have exchanged mass with a companion prior to explosion. We estimate that 1/3 to 1/2 of SN II progenitors have a history of mass exchange with a binary companion before exploding. The dominant binary channels leading to SN II progenitors involve the merger of binary stars. Mergers are expected to produce a diversity of SN II progenitor characteristics, depending on the evolutionary timing and properties of the merger. Alternatively, SN II progenitors from interacting binaries may have accreted mass from their companion, and subsequently been ejected from the binary system after their companion exploded. We show that the overall fraction of SN II progenitors that are predicted to have experienced binary interaction is robust against the main physical uncertainties in our models. However, the relative importance of different binary evolutionary channels is affected by changing physical assumptions. We further discuss ways in which binarity might contribute to the observed diversity of SNe II by considering potential observational signatures arising from each binary channel. For supernovae which have a substantial H-rich envelope at explosion (i.e., excluding Type IIb SNe), a surviving non-compact companion would typically indicate that the supernova progenitor star was in a wide, non-interacting binary. We argue that a significant fraction of even Type II-P SNe are expected to have gained mass from a companion prior to explosion. △ Less

Submitted 15 July, 2019; originally announced July 2019.

Comments: The paper has been accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 631, A5 (2019)

The Wolf-Rayet binaries of the nitrogen sequence in the Large Magellanic Cloud: spectroscopy, orbital analysis, formation, and evolution

Authors: T. Shenar, D. P. Sablowski, R. Hainich, H Todt, A. F. J. Moffat, L. M. Oskinova, V. Ramachandran, H. Sana, A. A. C Sander, O. Schnurr, N. St-Louis, D. Vanbeveren, Y. Goetberg, W. -R. Hamann

Abstract: Massive Wolf-Rayet (WR) stars dominate the radiative and mechanical energy budget of galaxies and probe a critical phase in the evolution of massive stars prior to core-collapse. It is not known whether core He-burning WR stars (classical WR, cWR) form predominantly through wind-stripping (w-WR) or binary stripping (b-WR). With spectroscopy of WR binaries so-far largely avoided due to its complexi… ▽ More Massive Wolf-Rayet (WR) stars dominate the radiative and mechanical energy budget of galaxies and probe a critical phase in the evolution of massive stars prior to core-collapse. It is not known whether core He-burning WR stars (classical WR, cWR) form predominantly through wind-stripping (w-WR) or binary stripping (b-WR). With spectroscopy of WR binaries so-far largely avoided due to its complexity, our study focuses on the 44 WR binaries / binary candidates of the Large Magellanic Cloud (LMC, metallicity Z~0.5 Zsun), identified on the basis of radial velocity variations, composite spectra, or high X-ray luminosities. Relying on a diverse spectroscopic database, we aim to derive the physical and orbital parameters of our targets, confronting evolution models of evolved massive stars at sub-solar metallicity, and constraining the impact of binary interaction in forming them. Spectroscopy is performed using the Potsdam Wolf-Rayet (PoWR) code and cross-correlation techniques. Disentanglement is performed using the code Spectangular or the shift-and-add algorithm. Evolutionary status is interpreted using the Binary Population and Spectral Synthesis (BPASS) code, exploring binary interaction and chemically-homogeneous evolution. No obvious dichotomy in the locations of apparently-single and binary WN stars on the Hertzsprung-Russell diagram is apparent. According to commonly used stellar evolution models (BPASS, Geneva), most apparently-single WN stars could not have formed as single stars, implying that they were stripped by an undetected companion. Otherwise, it must follow that pre-WR mass-loss/mixing (e.g., during the red supergiant phase) are strongly underestimated in standard stellar evolution models. △ Less

Submitted 8 September, 2020; v1 submitted 22 May, 2019; originally announced May 2019.

Comments: accepted to A&A on 10.05.2019; 69 pages (25 main paper + 44 appendix); Corrigendum: Shenar et al. 2020, A&A, 641, 2: An unfortunate typo in the implementation of the "transformed radius" caused errors of up to ~0.5dex in the derived mass-loss rates. This has now been corrected

Journal ref: A&A 641, C2 (2020)

Massive runaways and walkaway stars

Authors: M. Renzo, E. Zapartas, S. E. de Mink, Y. Götberg, S. Justham, R. J. Farmer, R. G. Izzard, S. Toonen, H. Sana

Abstract: Anticipating the kinematic constraints from the Gaia mission, we perform an extensive numerical study of the evolution of massive binary systems to predict the peculiar velocities that stars obtain when their companion collapses and disrupts the system. Our aim is to (1) identify which predictions are robust against model uncertainties and assess their implications, (2) investigate which physical… ▽ More Anticipating the kinematic constraints from the Gaia mission, we perform an extensive numerical study of the evolution of massive binary systems to predict the peculiar velocities that stars obtain when their companion collapses and disrupts the system. Our aim is to (1) identify which predictions are robust against model uncertainties and assess their implications, (2) investigate which physical processes leave a clear imprint and may therefore be constrained observationally and (3) provide a suite of publicly available model predictions. We find that $22_{-8}^{+26}$% of all massive binary systems merge prior to the first core collapse in the system. Of the remainder, $86_{-9}^{+11}$% become unbound because of the core-collapse. Remarkably, this rarely produce runaway stars (i.e., stars with velocities above 30 km/s). These are outnumbered by more than an order of magnitude by slower unbound companions, or "walkaway stars". This is a robust outcome of our simulations and is due to the reversal of the mass ratio prior to the explosion and widening of the orbit, as we show analytically and numerically. We estimate a $10^{+5}_{-8}$% of massive stars to be walkaways and only $0.5^{+1.0}_{-0.4}$% to be runaways, nearly all of which have accreted mass from their companion. Our findings are consistent with earlier studies, however the low runaway fraction we find is in tension with observed fractions 10%. If Gaia confirms these high fractions of massive runaway stars resulting from binaries, it would imply that we are currently missing physics in the binary models. Finally, we show that high end of the mass distributions of runaway stars is very sensitive to the assumed black hole natal kicks and propose this as a potentially stringent test for the explosion mechanism. We discuss companions remaining bound which can evolve into X-ray and gravitational wave sources. △ Less

Submitted 4 March, 2019; v1 submitted 24 April, 2018; originally announced April 2018.

Comments: accepted version: added sec. 8.1 on other possible contributions to the observed velocities of runaway stars

Journal ref: A&A 624, A66 (2019)

Clues about the scarcity of stripped-envelope stars from the evolutionary state of the sdO+Be binary system phi Persei

Authors: A. Schootemeijer, Y. Gotberg, S. E. de Mink, D. R. Gies, E. Zapartas

Abstract: Stripped-envelope stars (SESs) form in binary systems after losing mass through Roche-lobe overflow. They bear astrophysical significance as sources of UV and ionizing radiation in older stellar populations and, if sufficiently massive, as stripped supernova progenitors. Binary evolutionary models predict them to be common, but only a handful of subdwarfs (i.e., SESs) with B-type companions are kn… ▽ More Stripped-envelope stars (SESs) form in binary systems after losing mass through Roche-lobe overflow. They bear astrophysical significance as sources of UV and ionizing radiation in older stellar populations and, if sufficiently massive, as stripped supernova progenitors. Binary evolutionary models predict them to be common, but only a handful of subdwarfs (i.e., SESs) with B-type companions are known. This could be the result of observational biases hindering detection, or an incorrect understanding of binary evolution. We reanalyze the well-studied post-interaction binary phi Persei. Recently, new data improved the orbital solution of the system, which contains a ~1.2 Msun SES and a rapidly rotating ~9.6 Msun Be star. We compare with an extensive grid of evolutionary models using a Bayesian approach and find initial masses of the progenitor of 7.2+/-0.4 Msun for the SES and 3.8+/-0.4 Msun for the Be star. The system must have evolved through near-conservative mass transfer. These findings are consistent with earlier studies. The age we obtain, 57+/-9 Myr, is in excellent agreement with the age of the alpha Persei cluster. We note that neither star was initially massive enough to produce a core-collapse supernova, but mass exchange pushed the Be star above the mass threshold. We find that the subdwarf is overluminous for its mass by almost an order of magnitude, compared to the expectations for a helium core burning star. We can only reconcile this if the subdwarf is in a late phase of helium shell burning, which lasts only 2-3% of the total lifetime as a subdwarf. This could imply that up to ~50 less evolved, dimmer subdwarfs exist for each system similar to phi Persei. Our findings can be interpreted as a strong indication that a substantial population of SESs indeed exists, but has so far evaded detection because of observational biases and lack of large-scale systematic searches. △ Less

Submitted 6 March, 2018; originally announced March 2018.

Comments: 11 pages, 5 figures, accepted for publication in A&A

Spectral models for binary products: Unifying Subdwarfs and Wolf-Rayet stars as a sequence of stripped-envelope stars

Authors: Y. Götberg, S. E. de Mink, J. H. Groh, T. Kupfer, P. A. Crowther, E. Zapartas, M. Renzo

Abstract: Stars stripped of their envelope through interaction in a binary are generally not considered when accounting for ionizing radiation from stellar populations, despite the expectation that stripped stars emit hard ionizing radiation, form frequently and live 10-100 times longer than single massive stars. We compute the first grid of spectral models specially made for stars stripped in binaries for… ▽ More Stars stripped of their envelope through interaction in a binary are generally not considered when accounting for ionizing radiation from stellar populations, despite the expectation that stripped stars emit hard ionizing radiation, form frequently and live 10-100 times longer than single massive stars. We compute the first grid of spectral models specially made for stars stripped in binaries for a range of progenitor masses (2-20$M_{\odot}$) and metallicities covering a wide range. For stripped stars with masses between 0.3-7$M_{\odot}$, we find high effective temperatures (20-100 kK, increasing with mass), small radii (0.2-1$R_{\odot}$) and high bolometric luminosities, comparable to that of their progenitor before stripping. The spectra show a continuous sequence that naturally bridges subdwarf-type stars and Wolf-Rayet like spectra. For intermediate masses we find hybrid spectral classes showing a mixture of absorption and emission lines. These appear for stars with mass loss rates of $10^{-8}-10^{-6}M_{\odot}$/yr, which have semi-transparent atmospheres. At low metallicity, substantial H-rich layers are left at the surface and we predict spectra that resemble O-type stars instead. We obtain spectra undistinguishable from subdwarfs for stripped stars with masses up to 1.7$M_{\odot}$, which questions whether the widely adopted canonical value of 0.47$M_{\odot}$ is uniformly valid. Increasing the observed sample of stripped stars will provide necessary tests for the physics of interaction, internal mixing and stellar winds. We investigate the feasibility to detect stripped stars next to an optically bright companion and recommend searches for their UV excess and possible emission lines, most notably HeII4686 in the optical and HeII1640 in the UV. Our models are publicly available for further investigations or inclusion in spectral synthesis simulations. △ Less

Submitted 8 February, 2018; originally announced February 2018.

Comments: Accepted for publication in A&A

Journal ref: A&A 615, A78 (2018)

No surviving stellar companion for Cassiopeia A

Authors: Wolfgang E. Kerzendorf, Tuan Do, Selma E. de Mink, Ylva Götberg, Dan Millisaljevic, Emmanouil Zapartas, Mathieu Renzo, Stephen Justham, Philipp Podsiadlowski, Robert A. Fesen

Abstract: Massive stars in binaries can give rise to extreme phenomena such as X-ray binaries and gravitational wave sources after one or both stars end their lives as core-collapse supernovae. Stars in close orbit around a stellar or compact companion are expected to explode as "stripped-envelope supernovae", showing no (Type Ib/c) or little (Type IIb) signs of hydrogen in the spectra, because hydrogen-ric… ▽ More Massive stars in binaries can give rise to extreme phenomena such as X-ray binaries and gravitational wave sources after one or both stars end their lives as core-collapse supernovae. Stars in close orbit around a stellar or compact companion are expected to explode as "stripped-envelope supernovae", showing no (Type Ib/c) or little (Type IIb) signs of hydrogen in the spectra, because hydrogen-rich progenitors are too large to fit. The physical processes responsible for the stripping process and the fate of the companion are still very poorly understood. Aiming to find new clues, we investigate Cas~A, which is a very young ($\sim$340 \,yr) and near ($\sim$3.4\,kpc) remnant of a core collapse supernova. Cas~A has been subject to several searches for possible companions, all unsuccessfully. We present new measurements of the proper motions and photometry of stars in the vicinity based on deep HST ACS/WFC and WFC3-IR data. We identify stellar sources that are close enough in projection, but using their proper motions we show that none are compatible with being at the location of center at the time of explosion, in agreement with earlier findings. Our photometric measurements allow us to place much deeper (order of magnitude) upper limits on the brightness of possible undetected companions. We systematically compare them with model predictions for a wide variety of scenarios. We can confidently rule out the presence of any stellar companion of any reasonable mass and age (main sequence, pre main sequence or stripped). Although this finding is not in direct conflict with model predictions, it does rule out what many considered to be the most likely evolutionary scenario for type IIb supernova. (abstract abbreviated) △ Less

Submitted 5 December, 2018; v1 submitted 31 October, 2017; originally announced November 2017.

Comments: 13 pages, accepted by A&A

Journal ref: A&A 623, A34 (2019)

Extreme isolation of WN3/O3 stars and implications for their evolutionary origin as the elusive stripped binaries

Authors: Nathan Smith, Ylva Gotberg, Selma E. de Mink

Abstract: Recent surveys of the Magellanic Clouds have revealed a subtype of Wolf-Rayet (WR) star with peculiar properties. WN3/O3 spectra exhibit both WR-like emission and O3 V-like absorption - but at lower luminosity than O3 V or WN stars. We examine the projected spatial distribution of WN3/O3 stars in the LMC as compared to O-type stars. Surprisingly, WN3/O3 stars are among the most isolated of all cla… ▽ More Recent surveys of the Magellanic Clouds have revealed a subtype of Wolf-Rayet (WR) star with peculiar properties. WN3/O3 spectra exhibit both WR-like emission and O3 V-like absorption - but at lower luminosity than O3 V or WN stars. We examine the projected spatial distribution of WN3/O3 stars in the LMC as compared to O-type stars. Surprisingly, WN3/O3 stars are among the most isolated of all classes of massive stars; they have a distribution similar to red supergiants dominated by initial masses of 10-15 $M_{\odot}$, and are far more dispersed than classical WR stars or luminous blue variables (LBVs). Their lack of association with clusters of O-type stars suggests strongly that WN3/O3 stars are not the descendants of single massive stars (30 $M_{\odot}$ or above). Instead, they are likely products of interacting binaries at lower initial mass (10-18 $M_{\odot}$). Comparison with binary models suggests a probable origin with primaries in this mass range that were stripped of their H envelopes through non-conservative mass transfer by a low-mass secondary. We show that model spectra and positions on the Hertzsprung-Russell diagram for binary stripped stars are consistent with WN3/O3 stars. Monitoring radial velocities with high-resolution spectra can test for low-mass companions or runaway velocities. With lower initial mass and environments that avoid very massive stars, the WN3/O3 stars fit expectations for progenitors of Type Ib and possibly Type Ibn supernovae. △ Less

Submitted 5 December, 2017; v1 submitted 11 April, 2017; originally announced April 2017.

Comments: Accepted for publication in MNRAS

Ionizing spectra of stars that lose their envelope through interaction with a binary companion: role of metallicity

Authors: Y. Gotberg, S. E. de Mink, J. H. Groh

Abstract: Understanding ionizing fluxes of stellar populations is crucial for various astrophysical problems including the epoch of reionization. Massive short-lived Wolf-Rayet stars are generally considered as the main ionizing sources. We examine the role of less massive stars that lose their envelope through interaction with a companion. We use the evolutionary code MESA and the radiative transfer code C… ▽ More Understanding ionizing fluxes of stellar populations is crucial for various astrophysical problems including the epoch of reionization. Massive short-lived Wolf-Rayet stars are generally considered as the main ionizing sources. We examine the role of less massive stars that lose their envelope through interaction with a companion. We use the evolutionary code MESA and the radiative transfer code CMFGEN to investigate stripped stars as a function of metallicity (Z). We show that typical progenitors, initially 12$M_{\odot}$, produce hot and compact stars (~4$M_{\odot}$, 60-80 kK, ~1$R_{\odot}$) that copiously produce ionizing photons, emitting 60-85% and 30-60% of their energy as HI and HeI ionizing radiation, for Z=0.0001-0.02. This is comparable to a typical massive Wolf-Rayet star, when accounting for their longer lifetimes and the favoring initial mass function. Stripped stars further emit their ionizing photons with a delay (~20Myrs after star formation for this progenitor), allowing time for stellar feedback to disperse the birth clouds. This increases the fraction that can escape and contribute to ionization of the intergalactic medium. We further find that Roche stripping fails to fully remove the H-rich envelope at low Z. This questions the common treatment of stripped stars in rapid population synthesis simulations as pure helium stars. We expect implications for the rate of type Ib/c supernova at low Z and for the advanced evolutionary channels for type Ia supernova and gravitational wave sources. We discuss how the characteristic spectral features of stripped stars can be used to increase the observed sample, which is urgently needed to test the models. △ Less

Submitted 4 October, 2017; v1 submitted 25 January, 2017; originally announced January 2017.

Journal ref: A&A 608, A11 (2017)

Delay-time distribution of core-collapse supernovae with late events resulting from binary interaction

Authors: E. Zapartas, S. E. de Mink, R. G. Izzard, S. -C. Yoon, C. Badenes, Y. Gotberg, A. de Koter, C. J. Neijssel, M. Renzo, A. Schootemeijer, T. S. Shrotriya

Abstract: Most massive stars, the progenitors of core-collapse supernovae, are in close binary systems and may interact with their companion through mass transfer or merging. We undertake a population synthesis study to compute the delay-time distribution of core-collapse supernovae, that is, the supernova rate versus time following a starburst, taking into account binary interactions. We test the systemati… ▽ More Most massive stars, the progenitors of core-collapse supernovae, are in close binary systems and may interact with their companion through mass transfer or merging. We undertake a population synthesis study to compute the delay-time distribution of core-collapse supernovae, that is, the supernova rate versus time following a starburst, taking into account binary interactions. We test the systematic robustness of our results by running various simulations to account for the uncertainties in our standard assumptions. We find that a significant fraction, $15^{+9}_{-8}$%, of core-collapse supernovae are `late', that is, they occur 50-200 Myrs after birth, when all massive single stars have already exploded. These late events originate predominantly from binary systems with at least one, or, in most cases, with both stars initially being of intermediate mass ($4-8M_{\odot}$). The main evolutionary channels that contribute often involve either the merging of the initially more massive primary star with its companion or the engulfment of the remaining core of the primary by the expanding secondary that has accreted mass at an earlier evolutionary stage. Also, the total number of core-collapse supernovae increases by $14^{+15}_{-14}$% because of binarity for the same initial stellar mass. The high rate implies that we should have already observed such late core-collapse supernovae, but have not recognized them as such. We argue that $φ$ Persei is a likely progenitor and that eccentric neutron star - white dwarf systems are likely descendants. Late events can help explain the discrepancy in the delay-time distributions derived from supernova remnants in the Magellanic Clouds and extragalactic type Ia events, lowering the contribution of prompt Ia events. We discuss ways to test these predictions and speculate on the implications for supernova feedback in simulations of galaxy evolution. △ Less

Submitted 24 January, 2017; originally announced January 2017.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 601, A29 (2017)

Long-term stability of the HR 8799 planetary system without resonant lock

Authors: Ylva Gotberg, Melvyn B. Davies, Alexander Mustill, Anders Johansen, Ross P. Church

Abstract: HR 8799 is a star accompanied by four massive planets on wide orbits. The observed planetary configuration has been shown to be unstable on a timescale much shorter than the estimated age of the system (~ 30 Myr) unless the planets are locked into mean motion resonances. This condition is characterised by small-amplitude libration of one or more resonant angles that stabilise the system by prevent… ▽ More HR 8799 is a star accompanied by four massive planets on wide orbits. The observed planetary configuration has been shown to be unstable on a timescale much shorter than the estimated age of the system (~ 30 Myr) unless the planets are locked into mean motion resonances. This condition is characterised by small-amplitude libration of one or more resonant angles that stabilise the system by preventing close encounters. We simulate planetary systems similar to the HR 8799 planetary system, exploring the parameter space in separation between the orbits, planetary masses and distance from the Sun to the star. We find systems that look like HR 8799 and remain stable for longer than the estimated age of HR 8799. None of our systems are forced into resonances. We find, with nominal masses and in a narrow range of orbit separations, that 5 of 100 systems match the observations and lifetime. Considering a broad range of orbit separations, we find 12 of 900 similar systems. The systems survive significantly longer because of their slightly increased initial orbit separations compared to assuming circular orbits from the observed positions. A small increase in separation leads to a significant increase in survival time. The low eccentricity the orbits develop from gravitational interaction is enough for the planets to match the observations. With lower masses, but still comfortably within the estimated planet mass uncertainty, we find 18 of 100 matching and long-lived systems in a narrow orbital separation range. In the broad separation range, we find 82 of 900 matching systems. Our results imply that the planets in the HR 8799 system do not have to be in strong mean motion resonances. △ Less

Submitted 24 June, 2016; originally announced June 2016.

Comments: Accepted for publication in A&A

Journal ref: A&A 592, A147 (2016)