Showing 1–50 of 199 results for author: de Mink, S

Expected insights on type Ia supernovae from LISA's gravitational wave observations

Authors: Valeriya Korol, Riccardo Buscicchio, Ruediger Pakmor, Javier Morán-Fraile, Christopher J. Moore, Selma E. de Mink

Abstract: The nature of progenitors of type Ia supernovae have long been debated, primarily due to the elusiveness of the progenitor systems to traditional electromagnetic observation methods. We argue that gravitational wave observations with the upcoming Laser Interferometer Space Antenna (LISA) offer the most promising way to test one of the leading progenitor scenarios - the double-degenerate scenario,… ▽ More The nature of progenitors of type Ia supernovae have long been debated, primarily due to the elusiveness of the progenitor systems to traditional electromagnetic observation methods. We argue that gravitational wave observations with the upcoming Laser Interferometer Space Antenna (LISA) offer the most promising way to test one of the leading progenitor scenarios - the double-degenerate scenario, which involves a binary system of two white dwarf stars. In this study, we review published results, supplementing them with additional calculations for the context of type Ia supernovae. We discuss that LISA will be able to provide a complete sample of double white dwarf type Ia supernova progenitors with orbital period less than 16-11 minutes (gravitational wave frequencies above 2-3 milli-Hertz). Such a sample will enable a statistical validation of the double-degenerate scenario by simply counting whether LISA detects enough double white dwarf binaries to account the measured type Ia merger rate in Milky Way-like galaxies. Additionally, we illustrate how LISA's capability to measure the chirp mass will set lower bounds on the primary mass, revealing whether detected double white dwarf binaries will eventually result in a type Ia supernova. We estimate that the expected LISA constraints on the type Ia merger rate for the Milky Way will be 4-9%. We also discuss the potential gravitational wave signal from a type Ia supernova assuming a double detonation mechanism and explore how multi-messenger observations could significantly advance our understanding of these transient phenomena. △ Less

Submitted 4 July, 2024; originally announced July 2024.

Comments: submitted to A&A, comments are welcome

Large-scale ordered magnetic fields generated in mergers of helium white dwarfs

Authors: Rüdiger Pakmor, Ingrid Pelisoli, Stephen Justham, Abinaya S. Rajamuthukumar, Friedrich K. Röpke, Fabian R. N. Schneider, Selma E. de Mink, Sebastian T. Ohlmann, Philipp Podsiadlowski, Javier Moran Fraile, Marco Vetter, Robert Andrassy

Abstract: Stellar mergers are one important path to highly magnetised stars. Mergers of two low-mass white dwarfs may create up to every third hot subdwarf star. The merging process is usually assumed to dramatically amplify magnetic fields. However, so far only four highly magnetised hot subdwarf stars have been found, suggesting a fraction of less than $1\%$. We present two high-resolution magnetohydrod… ▽ More Stellar mergers are one important path to highly magnetised stars. Mergers of two low-mass white dwarfs may create up to every third hot subdwarf star. The merging process is usually assumed to dramatically amplify magnetic fields. However, so far only four highly magnetised hot subdwarf stars have been found, suggesting a fraction of less than $1\%$. We present two high-resolution magnetohydrodynamical (MHD) simulations of the merger of two helium white dwarfs in a binary system with the same total mass of $0.6\,M_\odot$. We analyse one equal-mass merger with two $0.3\,M_\odot$ white dwarfs, and one unequal-mass merger with a $0.25\,M_\odot$ white dwarf and a $0.35\,M_\odot$ white dwarf. We simulate the inspiral, merger, and further evolution of the merger remnant for about $50$ rotations. We find efficient magnetic field amplification in both mergers via a small-scale dynamo, reproducing previous results of stellar merger simulations. The magnetic field saturates at similar strength for both simulations. We then identify a second phase of magnetic field amplification in both merger remnants that happens on a timescale of several tens of rotational periods of the merger remnant. This phase generates a large-scale ordered azimuthal field. We identify it as a large-scale dynamo driven by the magneto-rotational instability (MRI). Finally, we suggest that in the unequal-mass merger remnant, helium burning will eventually start in a shell around a cold core. The convection zone this generates will coincide with the region that contains most of the magnetic energy, probably erasing the strong, ordered field. The equal-mass merger remnant instead will probably ignite burning in the center, retaining its ordered field. Therefore, the mass ratio of the initial merger could be the selecting factor that decides if a merger remnant will stay highly magnetised long after the merger. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 12 pages, 9 figures, submitted to A&A, comments welcome

Close Encounters of Wide Binaries Induced by the Galactic Tide: Implications for Stellar Mergers and Gravitational-Wave Sources

Authors: Jakob Stegmann, Alejandro Vigna-Gómez, Antti Rantala, Tom Wagg, Lorenz Zwick, Mathieu Renzo, Lieke A. C. van Son, Selma E. de Mink, Simon D. M. White

Abstract: A substantial fraction of stars can be found in wide binaries with projected separations between $\sim10^2$ and $10^5\,\rm AUえーゆー$. In the standard lore of binary physics, these would evolve as effectively single stars that remotely orbit one another on stationary Keplerian ellipses. However, embedded in their Galactic environment their low binding energy makes them exceptionally prone to perturbation… ▽ More A substantial fraction of stars can be found in wide binaries with projected separations between $\sim10^2$ and $10^5\,\rm AUえーゆー$. In the standard lore of binary physics, these would evolve as effectively single stars that remotely orbit one another on stationary Keplerian ellipses. However, embedded in their Galactic environment their low binding energy makes them exceptionally prone to perturbations from the gravitational potential of the Milky Way and encounters with passing stars. Employing a fully relativistic $N$-body integration scheme, we study the impact of these perturbations on the orbital evolution of wide binaries along their trajectory through the Milky Way. Our analysis reveals that the torques exerted by the Galaxy can cause large-amplitude oscillations of the binary eccentricity to $1-e\lesssim10^{-8}$. As a consequence, the wide binary members pass close to each other at periapsis, which, depending on the type of binary, potentially leads to a mass transfer or collision of stars or to an inspiral and subsequent merger of compact remnants due to gravitational-wave radiation. Based on a simulation of $10^5$ wide binaries across the Galactic field, we find that this mechanism could significantly contribute to the rate of stellar collisions and binary black hole mergers as inferred from observations of Luminous Red Novae and gravitational-wave events by LIGO/Virgo/Kagra. We conclude that the dynamics of wide binaries, despite their large mean separation, can give rise to extreme interactions between stars and compact remnants. △ Less

Submitted 5 May, 2024; originally announced May 2024.

Comments: 20 pages, 9 figures

The Asteroseismic Imprints of Mass Transfer: A Case Study of a Binary Mass Gainer in the SPB Instability Strip

Authors: Tom Wagg, Cole Johnston, Earl P. Bellinger, Mathieu Renzo, Richard Townsend, Selma E. de Mink

Abstract: We present new simulations investigating the impact of mass transfer on the asteroseismic signals of slowly pulsating B stars. We use MESA to simulate the evolution of a binary star system and GYRE to compute the asteroseismic properties of the accretor star. We show that, compared to a single star of the same final mass, a star that has undergone accretion (of non-enriched material) has a signifi… ▽ More We present new simulations investigating the impact of mass transfer on the asteroseismic signals of slowly pulsating B stars. We use MESA to simulate the evolution of a binary star system and GYRE to compute the asteroseismic properties of the accretor star. We show that, compared to a single star of the same final mass, a star that has undergone accretion (of non-enriched material) has a significantly different internal structure, evident in both the hydrogen abundance profile and Brunt-Väisälä frequency profile. These differences result in significant changes in the observed period spacing patterns, implying that one may use this as a diagnostic to test whether a star's core has been rejuvenated as a result of accretion. We show that it is essential to consider the full multimodal posterior distributions when fitting stellar properties of mass-gainers to avoid drawing misleading conclusions. Even with these considerations, stellar ages will be significantly underestimated when assuming single star evolution for a mass-gainer. We find that future detectors with improved uncertainties would rule out single star models with the correct mass and central hydrogen fraction. Our proof of principle analysis demonstrates the need to further investigate the impact of binary interactions on stellar asteroseismic signals for a wide range of parameters, such as initial mass, amount of mass transferred and the age of the accretor star at the onset of mass transfer. △ Less

Submitted 8 March, 2024; originally announced March 2024.

Comments: 10 pages (main text), 8 figures, submitted to A&A, comments welcome. Interactive plots available at https://www.tomwagg.com/html/interact/mass-gainer-asteroseismology.html

Investigating the Chemically Homogeneous Evolution Channel and its Role in the Formation of the Enigmatic Binary Black Hole Progenitor Candidate HD 5980

Authors: K. Sharpe, L. A. C. van Son, S. E. de Mink, R. Farmer, P. Marchant, G. Koenigsberger

Abstract: Chemically homogeneous evolution (CHE) is a promising channel for forming massive binary black holes. The enigmatic, massive Wolf-Rayet (WR) binary HD 5980 A&B has been proposed to have formed through this channel. We investigate this claim by comparing its observed parameters with CHE models. Using MESA, we simulate grids of close massive binaries then use a Bayesian approach to compare them with… ▽ More Chemically homogeneous evolution (CHE) is a promising channel for forming massive binary black holes. The enigmatic, massive Wolf-Rayet (WR) binary HD 5980 A&B has been proposed to have formed through this channel. We investigate this claim by comparing its observed parameters with CHE models. Using MESA, we simulate grids of close massive binaries then use a Bayesian approach to compare them with the stars' observed orbital period, masses, luminosities, and hydrogen surface abundances. The most probable models, given the observational data, have initial periods ~3 days, widening to the present-day ~20 day orbit as a result of mass loss -- correspondingly, they have very high initial stellar masses ($\gtrsim$150 M$_\odot$). We explore variations in stellar wind-mass loss and internal mixing efficiency, and find that models assuming enhanced mass-loss are greatly favored to explain HD 5980, while enhanced mixing is only slightly favoured over our fiducial assumptions. Our most probable models slightly underpredict the hydrogen surface abundances. Regardless of its prior history, this system is a likely binary black hole progenitor. We model its further evolution under our fiducial and enhanced wind assumptions, finding that both stars produce black holes with masses ~19-37 M$_\odot$. The projected final orbit is too wide to merge within a Hubble time through gravitational waves alone. However, the system is thought to be part of a 2+2 hierarchical multiple. We speculate that secular effects with the (possible) third and fourth companions may drive the system to promptly become a gravitational-wave source. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: 23 pages, 9 figures, 2 tables, 2 appendices. Accepted for publication in ApJ

Long-term Evolution of Binary Orbits Induced by Circumbinary Disks

Authors: Ruggero Valli, Christopher Tiede, Alejandro Vigna-Gómez, Jorge Cuadra, Magdalena Siwek, Jing-Ze Ma, Daniel J. D'Orazio, Jonathan Zrake, Selma E. de Mink

Abstract: Circumbinary disks are found in a variety of astrophysical scenarios, spanning binary star formation to accreting supermassive black hole binaries. The interaction with a circumbinary disk can yield opposite effects on the binary orbit leading to circularization, or exciting the eccentricity, widening the orbit or shrinking it and facilitating mergers. We present a new formalism for the long-term… ▽ More Circumbinary disks are found in a variety of astrophysical scenarios, spanning binary star formation to accreting supermassive black hole binaries. The interaction with a circumbinary disk can yield opposite effects on the binary orbit leading to circularization, or exciting the eccentricity, widening the orbit or shrinking it and facilitating mergers. We present a new formalism for the long-term evolution of the disk-binary interaction based on the results of recent suites of hydrodynamic simulations, which resolve the complex geometry of the gas in the vicinity of the binary and fully account for the gravitational and accretion forces. We release a python package, \texttt{spindler}, that implements our model. We show that, unless the mass reservoir feeding the disk is comparable to the mass of the binary, accretion onto the binary depletes the disk mass before inducing a significant change in orbital separation or mass ratio. This finding implies that, in most scenarios, interaction with a circumbinary disk is not an efficient mechanism to shrink the orbit of the binary. However, as long as the mass of the disk is at least a few percent of the mass of the binary, the interaction can excite the eccentricity up to an equilibrium value, and induce a statistical correlation between mass ratio and eccentricity. We consider the applicability of our model to a variety of astrophysical scenarios: during star formation, in evolved stellar binaries, triples and in supermassive black hole binaries. We discuss the theoretical and observational implications of our predictions. △ Less

Submitted 30 January, 2024; originally announced January 2024.

Comments: 18 pages, 12 figures. Submitted. Comments welcome. Spindler available as a Python package on Zenodo: https://doi.org/10.5281/zenodo.10529060

Solar evolution models with a central black hole

Authors: Earl P. Bellinger, Matt E. Caplan, Taeho Ryu, Deepika Bollimpalli, Warrick H. Ball, Florian Kühnel, R. Farmer, S. E. de Mink, Jørgen Christensen-Dalsgaard

Abstract: Hawking (1971) proposed that the Sun may harbor a primordial black hole whose accretion supplies some of the solar luminosity. Such an object would have formed within the first 1 s after the Big Bang with the mass of a moon or an asteroid. These light black holes are a candidate solution to the dark matter problem, and could grow to become stellar-mass black holes (BHs) if captured by stars. Here… ▽ More Hawking (1971) proposed that the Sun may harbor a primordial black hole whose accretion supplies some of the solar luminosity. Such an object would have formed within the first 1 s after the Big Bang with the mass of a moon or an asteroid. These light black holes are a candidate solution to the dark matter problem, and could grow to become stellar-mass black holes (BHs) if captured by stars. Here we compute the evolution of stars having such a BH at their center. We find that such objects can be surprisingly long-lived, with the lightest black holes having no influence over stellar evolution, while more massive ones consume the star over time to produce a range of observable consequences. Models of the Sun born about a BH whose mass has since grown to approximately $10^{-6}~\rm{M}_\odot$ are compatible with current observations. In this scenario, the Sun would first dim to half its current luminosity over a span of 100 Myr as the accretion starts to generate enough energy to quench nuclear reactions. The Sun would then expand into a fully-convective star, where it would shine luminously for potentially several Gyr with an enriched surface helium abundance, first as a sub-subgiant star, and later as a red straggler, before becoming a sub-solar-mass BH. We also present results for a range of stellar masses and metallicities. The unique internal structures of stars harboring BHs may make it possible for asteroseismology to discover them, should they exist. We conclude with a list of open problems and predictions. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: 20 pages, 6 figures, accepted for publication in The Astrophysical Journal

Is Betelgeuse Really Rotating? Synthetic ALMA Observations of Large-scale Convection in 3D Simulations of Red Supergiants

Authors: Jing-Ze Ma, Andrea Chiavassa, Selma E. de Mink, Ruggero Valli, Stephen Justham, Bernd Freytag

Abstract: The evolved stages of massive stars are poorly understood, but invaluable constraints can be derived from spatially resolved observations of nearby red supergiants, such as Betelgeuse. Atacama Large Millimeter/submillimeter Array (ALMA) observations of Betelgeuse showing a dipolar velocity field have been interpreted as evidence for a projected rotation rate of about $5\, \mathrm{km\, s^{-1}}$. Th… ▽ More The evolved stages of massive stars are poorly understood, but invaluable constraints can be derived from spatially resolved observations of nearby red supergiants, such as Betelgeuse. Atacama Large Millimeter/submillimeter Array (ALMA) observations of Betelgeuse showing a dipolar velocity field have been interpreted as evidence for a projected rotation rate of about $5\, \mathrm{km\, s^{-1}}$. This is 2 orders of magnitude larger than predicted by single-star evolution, which led to suggestions that Betelgeuse is a binary merger. We propose instead that large-scale convective motions can mimic rotation, especially if they are only partially resolved. We support this claim with 3D CO5BOLD simulations of nonrotating red supergiants that we postprocessed to predict ALMA images and SiO spectra. We show that our synthetic radial velocity maps have a 90% chance of being falsely interpreted as evidence for a projected rotation rate of $\mathrm{2\, km\, s^{-1}}$ or larger for our fiducial simulation. We conclude that we need at least another ALMA observation to firmly establish whether Betelgeuse is indeed rapidly rotating. Such observations would also provide insight into the role of angular momentum and binary interaction in the late evolutionary stages. The data will further probe the structure and complex physical processes in the atmospheres of red supergiants, which are immediate progenitors of supernovae and are believed to be essential in the formation of gravitational-wave sources. △ Less

Submitted 17 February, 2024; v1 submitted 28 November, 2023; originally announced November 2023.

Comments: 14 pages, 10 figures, 2 tables. Published in ApJL. Comments welcome. Animations accessible to general readers and supplementary materials can be found at Zenodo: doi: 10.5281/zenodo.10199936 . Postprocessing packages and Jupiter notebook to reproduce the figures will be available in the same link soon

Journal ref: ApJL 962 L36 (2024)

The Potential of Asteroseismology to Resolve the Blue Supergiant Problem

Authors: Earl Patrick Bellinger, Selma E. de Mink, Walter E. van Rossem, Stephen Justham

Abstract: Despite major progress in our understanding of massive stars, concerning discrepancies still remain between observations and theory. Most notable are the numerous stars observed beyond the theoretical main sequence, an evolutionary phase expected to be short-lived and hence sparsely populated. This is the "Blue Supergiant Problem." Stellar models with abnormal internal structures can provide long-… ▽ More Despite major progress in our understanding of massive stars, concerning discrepancies still remain between observations and theory. Most notable are the numerous stars observed beyond the theoretical main sequence, an evolutionary phase expected to be short-lived and hence sparsely populated. This is the "Blue Supergiant Problem." Stellar models with abnormal internal structures can provide long-lived solutions for this problem: core hydrogen-burning stars with oversized cores may explain the hotter ones, and core helium-burning stars with undersized cores may explain the cooler ones. Such stars may result from enhanced or suppressed mixing in single stars, or, more likely, as the products of binary interaction and stellar mergers. Here we investigate the potential of asteroseismology to uncover the nature of blue supergiants. We construct stellar models for the above scenarios and show that they predict $g$-mode period spacings that differ by an order of magnitude: $\sim$ 200~min versus $\sim$ 20~min for long-lived core H and He burning stars, respectively. For the classical scenario of H-shell-burning stars rapidly crossing the Hertzsprung gap, we furthermore predict changes of the order of $10^{-2}\,μみゅー\rm{Hz\,yr}^{-1}$ in high-frequency modes; this effect would be in principle observable from $\sim 5$~yr of asteroseismic monitoring if these modes can be identified. This raises the possibility of revealing the internal structure of blue supergiants, and thus determining whether these stars are indeed binary merger products. These asteroseismic diagnostics may be measurable through long time-series observations from the ongoing TESS mission and upcoming PLATO mission, thereby laying a path toward resolving the blue supergiant problem. △ Less

Submitted 27 May, 2024; v1 submitted 31 October, 2023; originally announced November 2023.

Comments: Accepted for publication in the Astrophysical Journal Letters

Explosions in Roche-Lobe Distorted Stars: Relativistic Bullets in Binaries

Authors: Marcus DuPont, Andrew MacFadyen, Selma de Mink

Abstract: State-of-the-art surveys reveal that most massive stars in the universe evolve in close binaries. Massive stars in such systems are expected to develop aspherical envelopes due to tidal interactions and/or rotational effects. Recently, it was shown that point explosions in oblate stars can produce relativistic equatorial ring-like outflows. Moreover, since stripped-envelope stars in binaries can e… ▽ More State-of-the-art surveys reveal that most massive stars in the universe evolve in close binaries. Massive stars in such systems are expected to develop aspherical envelopes due to tidal interactions and/or rotational effects. Recently, it was shown that point explosions in oblate stars can produce relativistic equatorial ring-like outflows. Moreover, since stripped-envelope stars in binaries can expand enough to fill their Roche lobes anew, it is likely that these stars die with a greater degree of asphericity than the oblate spheroid geometry previously studied. We investigate the effects of this asymmetry by studying the gas dynamics of axisymmetric point explosions in stars in various stages of filling their Roche lobes. We find that point explosions in these pear-shaped stars produce trans-relativistic ejecta that coalesces into bullets pointed both toward and away from the binary companion. We present this result and comment on key morphological differences between core-collapse explosions in spherical versus distorted stars in binary systems, effects on gravitational wave sources, and observational signatures that could be used to glean these explosion geometries from current and future surveys. △ Less

Submitted 31 October, 2023; originally announced October 2023.

Variability of Blue Supergiants in the LMC with TESS

Authors: Linhao Ma, Cole Johnston, Earl P. Bellinger, Selma E. de Mink

Abstract: The blue supergiant (BSG) problem, namely, the overabundance of BSGs inconsistent with classical stellar evolution theory, remains an open question in stellar astrophysics. Several theoretical explanations have been proposed, which may be tested by their predictions for the characteristic time variability. In this work, we analyze the light curves of a sample of 20 BSGs obtained from the Transitin… ▽ More The blue supergiant (BSG) problem, namely, the overabundance of BSGs inconsistent with classical stellar evolution theory, remains an open question in stellar astrophysics. Several theoretical explanations have been proposed, which may be tested by their predictions for the characteristic time variability. In this work, we analyze the light curves of a sample of 20 BSGs obtained from the Transiting Exoplanet Survey Satellite (TESS) mission. We report a characteristic signal in the low-frequency ($f\lesssim2\;\mathrm{day}^{-1}$) range for all our targets. The amplitude spectrum has a peak frequency at $\sim0.2\;\mathrm{day}^{-1}$, and we are able to fit it by a modified Lorentzian profile. The signal itself shows strong stochasticity across different TESS sectors, suggesting its driving mechanism happens on short ($\lesssim\mathrm{months}$) timescales. Our signals resemble those obtained for a limited sample of hotter OB stars and yellow supergiants, suggesting their possible common origins. We discuss three possible physical explanations: stellar winds launched by rotation, convection motions that reach the stellar surface, and waves from the deep stellar interior. The peak frequency of the signal favors processes related to the convective zone caused by the iron opacity peak, and the shape of the spectra might be explained by the propagation of high-order, damped gravity waves excited from that zone. We discuss the uncertainties and limitations of all these scenarios. △ Less

Submitted 28 April, 2024; v1 submitted 30 October, 2023; originally announced October 2023.

Comments: Accepted to ApJ

Simulating the tidal disruption of stars by stellar-mass black holes using moving-mesh hydrodynamics

Authors: Pavan Vynatheya, Taeho Ryu, Ruediger Pakmor, Selma E. de Mink, Hagai B. Perets

Abstract: In the centers of dense star clusters, close encounters between stars and compact objects are likely to occur. We study tidal disruption events of main-sequence (MS) stars by stellar-mass black holes (termed $μみゅー$TDEs), which can shed light on the processes occurring in these clusters, including being an avenue in the mass growth of stellar-mass BHs. Using the moving-mesh hydrodynamics code \texttt{… ▽ More In the centers of dense star clusters, close encounters between stars and compact objects are likely to occur. We study tidal disruption events of main-sequence (MS) stars by stellar-mass black holes (termed $μみゅー$TDEs), which can shed light on the processes occurring in these clusters, including being an avenue in the mass growth of stellar-mass BHs. Using the moving-mesh hydrodynamics code \texttt{AREPO}, we perform a suite of hydrodynamics simulations of partial $μみゅー$TDEs of realistic, \texttt{MESA}-generated MS stars by varying the initial mass of the star ($0.5\,{\rm M}_{\rm \odot}$ and $1\,{\rm M}_{\rm \odot}$), the age of the star (zero-age, middle-age and terminal-age), the mass of the black hole ($10\,{\rm M}_{\rm \odot}$ and $40\,{\rm M}_{\rm \odot}$) and the impact parameter (yielding almost no mass loss to full disruption). We then examine the dependence of the masses, spins, and orbital parameters of the partially disrupted remnant on the initial encounter parameters. We find that the mass lost from a star decreases exponentially with increasing distance of approach and that a $1\,{\rm M}_{\rm \odot}$ star loses lesser mass than a $0.5\,{\rm M}_{\rm \odot}$. Moreover, a more evolved star is less susceptible to mass loss. Tidal torques at the closest approach spin up the remnant by factors of $10^2$--$10^4$ depending on the impact parameter. The remnant star can be bound (eccentric) or unbound (hyperbolic) to the black hole: hyperbolic orbits occur when the star's central density concentration is relatively low and the black hole-star mass ratio is high, which is the case for the disruption of a $0.5\,{\rm M}_{\rm \odot}$ star. Finally, we provide best-fit analytical formulae for a range of parameters that can be incorporated into cluster codes to model star-black hole interaction more accurately. △ Less

Submitted 4 March, 2024; v1 submitted 23 October, 2023; originally announced October 2023.

Comments: Accepted for publication by A&A (with some revisions)

Journal ref: A&A 685, A45 (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 triples with chemically homogeneously evolving inner binaries

Authors: Andris Dorozsmai, Silvia Toonen, Alejandro Vigna-Gómez, Selma E. de Mink, Floris Kummer

Abstract: Observations suggest that massive stellar triples are common. However, their evolution is not yet fully understood. We investigate the evolution of hierarchical triples in which the stars of the inner binary experience chemically homogeneous evolution (CHE), particularly to understand the role of the tertiary star in the formation of gravitational-wave (GW) sources. We use the triple-star rapid po… ▽ More Observations suggest that massive stellar triples are common. However, their evolution is not yet fully understood. We investigate the evolution of hierarchical triples in which the stars of the inner binary experience chemically homogeneous evolution (CHE), particularly to understand the role of the tertiary star in the formation of gravitational-wave (GW) sources. We use the triple-star rapid population synthesis code TRES to determine the evolution of these systems at two representative metallicities: $Z = 0.005$ and $Z = 0.0005$. About half of all triples harbouring a CHE inner binary (CHE triples) experience tertiary mass transfer (TMT) episodes, an event which is rare for classically evolving stars. In the majority of TMT episodes, the inner binary consists of two main-sequence stars (58-60 per cent) or two black holes (BHs, 24-31 per cent). Additionally, we explore the role of von Zeipel-Lidov-Kozai (ZLK) oscillations for CHE triples. ZLK oscillations can result in eccentric stellar mergers or lead to the formation of eccentric compact binaries in systems with initial outer pericenters smaller than $\sim$ 1200 $R_{\odot}$. Approximately 24-30 per cent of CHE triples form GW sources, and in 31 per cent of these, the tertiary star plays a significant role and leads to configurations that are not predicted for isolated binaries. We conclude that the evolution of CHE binaries can be affected by a close tertiary companion, resulting in astronomical transients such as BH-BH binaries that merge via GW emission orders of magnitude faster than their isolated binary counterparts and tertiary-driven massive stellar mergers. △ Less

Submitted 11 December, 2023; v1 submitted 10 July, 2023; originally announced July 2023.

Comments: 28 pages, 17 figures, 3 tables. Accepted for publication in MNRAS

Close Encounters of Star - Black Hole Binaries with Single Stars

Authors: Taeho Ryu, Selma de Mink, Rob Farmer, Ruediger Pakmor, Rosalba Perna, Volker Springel

Abstract: Multi-body dynamical interactions of binaries with other objects are one of the main driving mechanisms for the evolution of star clusters. It is thus important to bring our understanding of three-body interactions beyond the commonly employed point-particle approximation. To this end we here investigate the hydrodynamics of three-body encounters between star-black hole (BH) binaries and single st… ▽ More Multi-body dynamical interactions of binaries with other objects are one of the main driving mechanisms for the evolution of star clusters. It is thus important to bring our understanding of three-body interactions beyond the commonly employed point-particle approximation. To this end we here investigate the hydrodynamics of three-body encounters between star-black hole (BH) binaries and single stars, focusing on the identification of final outcomes and their long-term evolution and observational properties, using the moving-mesh hydrodynamics code AREPO. This type of encounters produces five types of outcomes: stellar disruption, stellar collision, weak perturbation of the original binary, binary member exchange, and triple formation. The two decisive parameters are the binary phase angle, which determines which two objects meet at the first closest approach, and the impact parameter, which sets the boundary between violent and non-violent interactions. When the impact parameter is smaller than the semimajor axis of the binary, tidal disruptions and star-BH collisions frequently occur when the BH and the incoming star first meet, while the two stars mostly merge when the two stars meet first instead. In both cases, the BHs accrete from an accretion disk at super-Eddington rates, possibly generating flares luminous enough to be observed. The stellar collision products either form a binary with the BH or remain unbound to the BH. Upon collision, the merged stars are hotter and larger than main sequence stars of the same mass at similar age. Even after recovering their thermal equilibrium state, stellar collision products, if isolated, would remain hotter and brighter than main sequence stars until becoming giants. △ Less

Submitted 10 July, 2023; v1 submitted 6 July, 2023; originally announced July 2023.

Comments: 16 pages, 10 figures, 2 tables. Submitted to MNRAS. Movie here: https://studio.youtube.com/playlist/PLxLK3qI02cQf_x96DSKQWVS251Jog7IOE/videos

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)

Close encounters of black hole - star binaries with stellar-mass black holes

Authors: Taeho Ryu, Ruggero Valli, Rudiger Pakmor, Rosalba Perna, Selma E. de Mink, Volker Springel

Abstract: Dynamical interactions involving binaries play a crucial role in the evolution of star clusters and galaxies. We continue our investigation of the hydrodynamics of three-body encounters, focusing on binary black hole (BBH) formation, stellar disruption, and electromagnetic (EM) emission in dynamical interactions between a BH-star binary and a stellar-mass BH, using the moving-mesh hydrodynamics co… ▽ More Dynamical interactions involving binaries play a crucial role in the evolution of star clusters and galaxies. We continue our investigation of the hydrodynamics of three-body encounters, focusing on binary black hole (BBH) formation, stellar disruption, and electromagnetic (EM) emission in dynamical interactions between a BH-star binary and a stellar-mass BH, using the moving-mesh hydrodynamics code {\small AREPO}. This type of encounters can be divided into two classes depending on whether the final outcome includes BBHs. This outcome is primarily determined by which two objects meet at the first closest approach. BBHs are more likely to form when the star and the incoming BH encounter first with an impact parameter smaller than the binary's semimajor axis. In this case, the star is frequently disrupted. On the other hand, when the two BHs encounter first, frequent consequences are an orbit perturbation of the original binary or a binary member exchange. For the parameters chosen in this study, BBH formation, accompanied by stellar disruption, happens in roughly 1 out of 4 encounters. The close correlation between BBH formation and stellar disruption has possible implications for EM counterparts at the binary's merger. The BH that disrupts the star is promptly surrounded by an optically and geometrically thick disk with accretion rates exceeding the Eddington limit. If the debris disk cools fast enough to become long-lived, EM counterparts can be produced at the time of the BBH merger. △ Less

Submitted 10 July, 2023; v1 submitted 4 April, 2023; originally announced April 2023.

Comments: 16 pages, 11 figures, 2 table. Accepted for publication in MNRAS, movie here: https://studio.youtube.com/playlist/PLxLK3qI02cQfXhvhRZp9gkvts_E46t0kH/videos

Nucleosynthesis of binary-stripped stars

Authors: R. Farmer, E. Laplace, Jing-ze Ma, S. E. de Mink, S. Justham

Abstract: The cosmic origin of the elements, the fundamental chemical building blocks of the Universe, is still uncertain. Binary interactions play a key role in the evolution of many massive stars, yet their impact on chemical yields is poorly understood. Using the MESA stellar evolution code we predict the chemical yields ejected in wind mass loss and the supernovae of single and binary-stripped stars. We… ▽ More The cosmic origin of the elements, the fundamental chemical building blocks of the Universe, is still uncertain. Binary interactions play a key role in the evolution of many massive stars, yet their impact on chemical yields is poorly understood. Using the MESA stellar evolution code we predict the chemical yields ejected in wind mass loss and the supernovae of single and binary-stripped stars. We do this with a large 162 isotope nuclear network at solar-metallicity. We find that binary-stripped stars are more effective producers of the elements than single stars, due to their increased mass loss and an increased chance to eject their envelopes during a supernova. This increased production by binaries varies across the periodic table, with Fluorine and Potassium being more significantly produced by binary-stripped stars than single stars. We find that the C12/C13 could be used as an indicator of the conservativeness of mass transfer, as C13 is preferentially ejected during mass transfer while C12 is preferentially ejected during wind mass loss. We identify a number of gamma-ray emitting radioactive isotopes that may be used to help constrain progenitor and explosion models of core-collapse supernovae with next-generation gamma-ray detectors. For single stars we find V44 and Mn52 are strong probes of the explosion model, while for binary-stripped stars it is Cr48. Our findings highlight that binary-stripped stars are not equivalent to two single stars and that detailed stellar modelling is needed to predict their final nucleosynthetic yields. △ Less

Submitted 8 March, 2023; originally announced March 2023.

Comments: 26 pages, 8 figures, 6 Tables, Online data is available at https://zenodo.org/record/5929871, Accepted for publication in ApJ

Spectroscopic and evolutionary analyses of the binary system AzV 14 outline paths toward the WR stage at low metallicity

Authors: D. Pauli, L. M. Oskinova, W. -R. Hamann, D. M. Bowman, H. Todt, T. Shenar, A. A. C. Sander, C. Erba, V. M. A. Gómez-González, C. Kehrig, J. Klencki, R. Kuiper, A. Mehner, S. E. de Mink, M. S. Oey, V. Ramachandran, A. Schootemeijer, S. Reyero Serantes, A. Wofford

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

Submitted 7 March, 2023; originally announced March 2023.

Comments: 21 pages (13 main body + 8 appendix), 16 figures, 9 tables

Journal ref: A&A 673, A40 (2023)

The initial spin distribution of B-type stars revealed by the split main sequences of young star clusters

Authors: Chen Wang, Ben Hastings, Abel Schootemeijer, Norbert Langer, Selma E. de Mink, Julia Bodensteiner, Antonino Milone, Stephen Justham, Pablo Marchant

Abstract: Spectroscopic observations of stars in young open clusters have revealed evidence for a dichotomous distribution of stellar rotational velocities, with 10-30% of stars rotating slowly and the remaining 70-90% rotating fairly rapidly. At the same time, high-precision multiband photometry of young star clusters shows a split main sequence band, which is again interpreted as due to a spin dichotomy.… ▽ More Spectroscopic observations of stars in young open clusters have revealed evidence for a dichotomous distribution of stellar rotational velocities, with 10-30% of stars rotating slowly and the remaining 70-90% rotating fairly rapidly. At the same time, high-precision multiband photometry of young star clusters shows a split main sequence band, which is again interpreted as due to a spin dichotomy. Recent papers suggest that extreme rotation is required to retrieve the photometric split. Our new grids of MESA models and the prevalent SYCLIST models show, however, that initial slow (0-35% of the linear Keplerian rotation velocities) and intermediate (50-65% of the Keplerian rotation velocities) rotation are adequate to explain the photometric split. These values are consistent with the recent spectroscopic measurements of cluster and field stars, and are likely to reflect the birth spin distributions of upper main-sequence stars. A fraction of the initially faster-rotating stars may be able to reach near-critical rotation at the end of their main-sequence evolution and produce Be stars in the turn-off region of young star clusters. However, we find that the presence of Be stars up to two magnitudes below the cluster turnoff advocates for a crucial role of binary interaction in creating Be stars. We argue that surface chemical composition measurements may help distinguish these two Be star formation channels. While only the most rapidly rotating, and therefore nitrogen-enriched, single stars can evolve into Be stars, slow pre-mass-transfer rotation and inefficient accretion allows for mild or no enrichment even in critically rotating accretion-induced Be stars. Our results shed new light on the origin of the spin distribution of young and evolved B-type main sequence stars. △ Less

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

Comments: 18 pages, 14 figures. Accepted for publication in A&A

Journal ref: A&A 670, A43 (2023)

The VLT-FLAMES Tarantula Survey: Observational evidence for two distinct populations of massive runaway stars in 30 Doradus

Authors: H. Sana, O. H. Ramírez-Agudelo, V. Hénault-Brunet, L. Mahy, L. A. Almeida, A. de Koter, J. M. Bestenlehner, C. J. Evans, N. Langer, F. R. N. Schneider, P. A. Crowther, S. E. de Mink, A. Herrero, D. J. Lennon, M. Gieles, J. Maíz Apellániz, M. Renzo, E. Sabbi, J. Th. van Loon, J. S. Vink

Abstract: Two main scenarios have been proposed for origin of massive runaway stars -- dynamical ejection or release from a binary at the first core collapse -- but their relative contribution remains debated. Using two large spectroscopic campaigns towards massive stars in 30 Doradus, we aim to provide observational constraints on the properties of the O-type runaway population in the most massive active… ▽ More Two main scenarios have been proposed for origin of massive runaway stars -- dynamical ejection or release from a binary at the first core collapse -- but their relative contribution remains debated. Using two large spectroscopic campaigns towards massive stars in 30 Doradus, we aim to provide observational constraints on the properties of the O-type runaway population in the most massive active star-forming region in the Local group. We use RV measurements of the O-type star populations in 30 Doradus obtained by the VLT-FLAMES Tarantula Survey and the Tarantula Massive Binary Monitoring to identify single and binary O-type runaways. We discuss their rotational properties and qualitatively compare observations with expectations of ejection scenarios. We identify 23 single and one binary O-type runaway objects, most of them outside the main star-forming regions in 30 Doradus. We find an overabundance of rapid rotators (vsini > 200km/s) among the runaway population, providing an explanation of the overabundance of rapidly rotating stars in the 30 Doradus field. Considerations of the projected rotation rates and runaway line-of-sight (los) velocities reveal a conspicuous absence of rapidly rotating (vsini > 210k/ms), fast moving (v_{los} > 60km/s) runaways, and suggest the presence of two different populations of runaway stars: a population of rapidly-spinning but slowly moving runaways and a population of fast moving but slowly rotating ones. These are detected with a ratio close to 2:1 in our sample. We argue that slowly moving but rapidly spinning runaways result from binary ejections, while rapidly moving but slowly spinning runaways could result from dynamical ejections. Given that detection biases will more strongly impact the slow-moving population, our results suggest that the binary evolution scenario dominates the current massive runaway population in 30 Doradus. △ Less

Submitted 24 November, 2022; originally announced November 2022.

Comments: Accepted for publication in A&A Letters; 9 pages, 5 figures

Journal ref: A&A 668, L5 (2022)

JWST/NIRCam Probes Young Star Clusters in the Reionization Era Sunrise Arc

Authors: E. Vanzella, A. Claeyssens, B. Welch, A. Adamo, D. Coe, J. M. Diego, G. Mahler, G. Khullar, V. Kokorev, M. Oguri, S. Ravindranath, L. J. Furtak, T. Yu-Yang Hsiao, Abdurro'uf, N. Mandelker, G. Brammer, L. D. Bradley, M. Bradac, C. J. Conselice, P. Dayal, M. Nonino, F. Andrade-Santos, R. A. Windhorst, N. Pirzkal, K. Sharon , et al. (5 additional authors not shown)

Abstract: Star cluster formation in the early universe and their contribution to reionization remains to date largely unconstrained. Here we present JWST/NIRCam imaging of the most highly magnified galaxy known at z ~ 6, the Sunrise arc. We identify six young massive star clusters (YMCs) with measured radii spanning ~ 20 pc down to ~ 1 pc (corrected for lensing magnification), estimated stellar masses of ~… ▽ More Star cluster formation in the early universe and their contribution to reionization remains to date largely unconstrained. Here we present JWST/NIRCam imaging of the most highly magnified galaxy known at z ~ 6, the Sunrise arc. We identify six young massive star clusters (YMCs) with measured radii spanning ~ 20 pc down to ~ 1 pc (corrected for lensing magnification), estimated stellar masses of ~ $10^{(6-7)}$ Msun, and with ages 1-30 Myr based on SED fitting to photometry measured in 8 filters extending to rest-frame 7000A. The resulting stellar mass surface densities are higher than 1000 Msun pc$^{-2}$ (up to a few $10^5$ Msun pc$^{-2}$) and their inferred dynamical ages qualify the majority of these systems as gravitationally-bound stellar clusters. The star cluster ages map the progression of star formation along the arc, with to evolved systems (>~ 10 Myr old) followed by very young clusters. The youngest stellar clusters (< 5 Myr) show evidence of prominent Hbeta + [OIII]4959,5007 emission, based on photometry, with equivalent widths larger than 1000 A rest-frame, and are hosted in a 200 pc sized star-forming complex. Such a region dominates the ionizing photon production, with a high efficiency log($ξくしー_{ion}$ [Hz erg$^{-1}$]) ~ 25.7. A significant fraction of the recently formed stellar mass of the galaxy (> 10-30 %) occurred in these YMCs. We speculate that such sources of ionizing radiation boost the ionizing photon production efficiency which eventually carve ionized channels that might favor the escape of Lyman continuum radiation. The survival of some of the clusters would make them the progenitors of massive and relatively metal-poor globular clusters in the local Universe. △ Less

Submitted 7 February, 2023; v1 submitted 17 November, 2022; originally announced November 2022.

Comments: 18 pages, 8 figures, 3 tables. ApJ, Accepted

Triple evolution: an important channel in the formation of type Ia supernovae

Authors: Abinaya Swaruba Rajamuthukumar, Adrian Hamers, Patrick Neunteufel, Ruediger Pakmor, Selma E. de mink

Abstract: Type Ia supernovae (SNe Ia) are thought to be the result of thermonuclear explosions in white dwarfs (WDs). Commonly considered formation pathways include two merging WDs (the double degenerate channel), and a single WD accreting material from a H or He donor (the single degenerate channel). Since the predicted SNe Ia rates from WD in binaries are thought to be insufficient to explain the observed… ▽ More Type Ia supernovae (SNe Ia) are thought to be the result of thermonuclear explosions in white dwarfs (WDs). Commonly considered formation pathways include two merging WDs (the double degenerate channel), and a single WD accreting material from a H or He donor (the single degenerate channel). Since the predicted SNe Ia rates from WD in binaries are thought to be insufficient to explain the observed SNe Ia rate, it is important to study similar interactions in higher-order multiple star systems such as triple systems. We use the evolutionary population synthesis code Multiple Stellar Evolution (MSE) to study stellar evolution, binary interactions and gravitational dynamics of the triple-star systems. Also, unlike previous studies, prescriptions are included to simultaneously take into account the single and double degenerate channels, and we consider triples across the entire parameter space (including those with tight inner binaries). We explore the impact of typically ignored or uncertain physics such as fly-bys and CE prescription parameters on our results. The majority of systems undergo circular mergers to explode as SNe Ia, while eccentric collisions contribute to $0.4-4$ per cent of SNe Ia events. The time-integrated SNe Ia rate from the triple channel is found to be $(3.60 \pm 0.04) \times {10^{-4}}\,\mathrm{M_{\odot}}^{-1}$ which is, surprisingly, similar to that of the isolated binary channel where the SNe Ia rate is $(3.2 \pm 0.1) \times {10^{-4}}\,\mathrm{M_{\odot}}^{-1}$. This implies that triples, when considering their entire parameter space, yield an important contribution to the overall SNe Ia rate. △ Less

Submitted 24 March, 2023; v1 submitted 8 November, 2022; originally announced November 2022.

Comments: accepted for Publication in APJ

Close Encounters of Tight Binary Stars with Stellar-mass Black Holes

Authors: Taeho Ryu, Rosalba Perna, Ruediger Pakmor, Jing-Ze Ma, Rob Farmer, Selma E. de Mink

Abstract: Strong dynamical interactions among stars and compact objects are expected in a variety of astrophysical settings, such as star clusters and the disks of active galactic nuclei. Via a suite of 3D hydrodynamics simulations using the moving-mesh code AREPO, we investigate the formation of transient phenomena and their properties in close encounters between an $2M_{\odot}$ or $20M_{\odot}$ equal-mass… ▽ More Strong dynamical interactions among stars and compact objects are expected in a variety of astrophysical settings, such as star clusters and the disks of active galactic nuclei. Via a suite of 3D hydrodynamics simulations using the moving-mesh code AREPO, we investigate the formation of transient phenomena and their properties in close encounters between an $2M_{\odot}$ or $20M_{\odot}$ equal-mass circular binary star and single $20M_{\odot}$ black hole (BH). Stars can be disrupted by the BH during dynamical interactions, naturally producing electromagnetic transient phenomena. Encounters with impact parameters smaller than the semimajor axis of the initial binary frequently lead to a variety of transients whose electromagnetic signatures are qualitatively different from those of ordinary disruption events involving just two bodies. These include the simultaneous or successive disruptions of both stars and one full disruption of one star accompanied by successive partial disruptions of the other star. On the other hand, when the impact parameter is larger than the semimajor axis of the initial binary, the binary is either simply tidally perturbed or dissociated into bound and unbound single stars ("micro-Hills" mechanism). The dissociation of $20M_{\odot}$ binaries can produce a runaway star and an active BH moving away from one another. Also, the binary dissociation can either produce an interacting binary with the BH, or a non-interacting, hard binary; both could be candidates of BH high- and low-mass X-ray binaries. Hence our simulations especially confirm that strong encounters can lead to the formation of the (generally difficult to form) BH low-mass X-ray binaries. △ Less

Submitted 5 January, 2023; v1 submitted 4 November, 2022; originally announced November 2022.

Comments: 13 pages, 8 figures, 2 tables. Accepted for publication in MNRAS. Movies: https://www.youtube.com/watch?v=DGz2WokfZLA&list=PLxLK3qI02cQeEc3ktlYAOLY4RHVX1V_nZ

The effects of stellar rotation along the main sequence of the 100 Myr old massive cluster NGC 1850

Authors: Sebastian Kamann, Sara Saracino, Nate Bastian, Seth Gossage, Christopher Usher, Dietrich Baade, Ivan Cabrera-Ziri, Selma E. de Mink, Sylvia Ekström, Cyril Georgy, Michael Hilker, Søren S. Larsen, Dougal Mackey, Florian Niederhofer, Imants Platais, David Yong

Abstract: Young star clusters enable us to study the effects of stellar rotation on an ensemble of stars of the same age and across a wide range in stellar mass and are therefore ideal targets for understanding the consequences of rotation on stellar evolution. We combine MUSE spectroscopy with HST photometry to measure the projected rotational velocities (Vsini) of 2,184 stars along the split main sequence… ▽ More Young star clusters enable us to study the effects of stellar rotation on an ensemble of stars of the same age and across a wide range in stellar mass and are therefore ideal targets for understanding the consequences of rotation on stellar evolution. We combine MUSE spectroscopy with HST photometry to measure the projected rotational velocities (Vsini) of 2,184 stars along the split main sequence and on the main sequence turn-off (MSTO) of the 100 Myr-old massive (10^5 M_sun) star cluster NGC 1850 in the Large Magellanic Cloud. At fixed magnitude, we observe a clear correlation between Vsini and colour, in the sense that fast rotators appear redder. The average Vsini values for stars on the blue and red branches of the split main sequence are ~100 km/s and ~200 km/s, respectively. The values correspond to about 25-30% and 50-60% of the critical rotation velocity and imply that rotation rates comparable to those observed in field stars of similar masses can explain the split main sequence. Our spectroscopic sample contains a rich population of ~200 fast rotating Be stars. The presence of shell features suggests that 23% of them are observed through their decretion disks, corresponding to a disk opening angle of 15 degrees. These shell stars can significantly alter the shape of the MSTO, hence care should be taken when interpreting this photometric feature. Overall, our findings impact our understanding of the evolution of young massive clusters and provide new observational constraints for testing stellar evolutionary models. △ Less

Submitted 12 December, 2022; v1 submitted 1 November, 2022; originally announced November 2022.

Comments: Published in MNRAS; 16 pages, 11 figures; data available via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/MNRAS/518/1505#/article

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 518, Issue 1, January 2023, Pages 1505-1521

Searching for compact objects in the single-lined spectroscopic binaries of the young Galactic cluster NGC 6231

Authors: G. Banyard, L. Mahy, H. Sana, J. Bodensteiner, J. I. Villaseńor, K. Sen, N. Langer, S. de Mink, A. Picco, T. Shenar

Abstract: Recent evolutionary computations predict that a few percent of massive OB stars in binary systems should have a dormant BH companion. Despite several reported X-ray quiet OB+BH systems over the last couple of years, finding them with certainty remains challenging. These have great importance as they can be gravitational wave (GW) source progenitors, and are landmark systems in constraining superno… ▽ More Recent evolutionary computations predict that a few percent of massive OB stars in binary systems should have a dormant BH companion. Despite several reported X-ray quiet OB+BH systems over the last couple of years, finding them with certainty remains challenging. These have great importance as they can be gravitational wave (GW) source progenitors, and are landmark systems in constraining supernova kick physics. This work aims to characterise the hidden companions to the single-lined spectroscopic binaries (SB1s) in the B star population of the young open Galactic cluster NGC 6231 to find candidate systems for harbouring compact object companions. With the orbital solutions for each SB1 previously constrained, we applied Fourier spectral disentangling to multi-epoch optical VLT/FLAMES spectra of each target to extract a potential signature of a faint companion, and to identify newly disentangled double-lined spectroscopic binaries (SB2s). For targets where the disentangling does not reveal any signature of a stellar companion, we performed atmospheric and evolutionary modelling on the primary to obtain constraints on the unseen companion. Seven newly classified SB2 systems with mass ratios down to near 0.1 were identified. From the remaining targets, for which no faint companion could be extracted from the spectra, four are found to have companion masses in the predicted mass ranges of neutron stars (NSes) and BHes. Two of these have companion masses between 1 and 3.5 $M_{\odot}$, making them potential hosts of NSes (or lower mass main sequence stars). The other two are between 2.5 to 8 $M_{\odot}$ and 1.6 and 26 $M_{\odot}$, respectively, and so are identified as candidates for harbouring BH companions. However, unambiguous identification of these systems as X-ray quiet compact object harbouring binaries requires follow up observations. △ Less

Submitted 13 October, 2022; originally announced October 2022.

Comments: Submitted to A&A

Journal ref: A&A 674, A60 (2023)

No peaks without valleys: The stable mass transfer channel for gravitational-wave sources in light of the neutron star-black hole mass gap

Authors: L. A. C. van Son, S. E. de Mink, M. Renzo, S. Justham, E. Zapartas, K. Breivik, T. Callister, W. M. Farr, C. Conroy

Abstract: Gravitational-wave (GW) detections are starting to reveal features in the mass distribution of double compact objects. The lower end of the black hole (BH) mass distribution is especially interesting as few formation channels contribute here and because it is more robust against variations in the cosmic star formation than the high mass end. In this work we explore the stable mass transfer channel… ▽ More Gravitational-wave (GW) detections are starting to reveal features in the mass distribution of double compact objects. The lower end of the black hole (BH) mass distribution is especially interesting as few formation channels contribute here and because it is more robust against variations in the cosmic star formation than the high mass end. In this work we explore the stable mass transfer channel for the formation of GW sources with a focus on the low-mass end of the mass distribution. We conduct an extensive exploration of the uncertain physical processes that impact this channel. We note that, for fiducial assumptions, this channel reproduces the peak at $\sim9 \mathrm{M_{\odot}}$ in the GW-observed binary BH mass distribution remarkably well, and predicts a cutoff mass that coincides with the upper edge of the purported neutron star BH mass gap. The peak and cutoff mass are a consequence of unique properties of this channel, namely (1) the requirement of stability during the mass transfer phases, and (2) the complex way in which the final compact object masses scale with the initial mass. We provide an analytical expression for the cutoff in the primary component mass and show that this adequately matches our numerical results. Our results imply that selection effects resulting from the formation channel alone can provide an explanation for the purported neutron star--BH mass gap in GW detections. This provides an alternative to the commonly adopted view that the gap emerges during BH formation. △ Less

Submitted 10 November, 2022; v1 submitted 27 September, 2022; originally announced September 2022.

Comments: Accepted for publication in ApJ associated code is available at https://github.com/LiekeVanSon/LowMBH_and_StableChannel

The locations of features in the mass distribution of merging binary black holes are robust against uncertainties in the metallicity-dependent cosmic star formation history

Authors: L. A. C. van Son, S. E. de Mink, M. Chruslinska, C. Conroy, R. Pakmor, L. Hernquist

Abstract: New observational facilities are probing astrophysical transients such as stellar explosions and gravitational wave (GW) sources at ever increasing redshifts, while also revealing new features in source property distributions. To interpret these observations, we need to compare them to predictions from stellar population models. Such models require the metallicity-dependent cosmic star formation h… ▽ More New observational facilities are probing astrophysical transients such as stellar explosions and gravitational wave (GW) sources at ever increasing redshifts, while also revealing new features in source property distributions. To interpret these observations, we need to compare them to predictions from stellar population models. Such models require the metallicity-dependent cosmic star formation history ($\mathcal{S}(Z,z)$) as an input. Large uncertainties remain in the shape and evolution of this function. In this work, we propose a simple analytical function for $\mathcal{S}(Z,z)$. Variations of this function can be easily interpreted, because the parameters link to its shape in an intuitive way. We fit our analytical function to the star-forming gas of the cosmological TNG100 simulation and find that it is able to capture the main behaviour well. As an example application, we investigate the effect of systematic variations in the $\mathcal{S}(Z,z)$ parameters on the predicted mass distribution of locally merging binary black holes (BBH). Our main findings are: I) the locations of features are remarkably robust against variations in the metallicity-dependent cosmic star formation history, and II) the low mass end is least affected by these variations. This is promising as it increases our chances to constrain the physics that governs the formation of these objects. △ Less

Submitted 27 February, 2023; v1 submitted 7 September, 2022; originally announced September 2022.

Comments: Accepted for publication in ApJ, 6 Figures, code available at https://github.com/LiekeVanSon/SFRD_fit, made using showyourwork

JWST Imaging of Earendel, the Extremely Magnified Star at Redshift $z=6.2$

Authors: Brian Welch, Dan Coe, Erik Zackrisson, S. E. de Mink, Swara Ravindranath, Jay Anderson, Gabriel Brammer, Larry Bradley, Jinmi Yoon, Patrick Kelly, Jose M. Diego, Rogier Windhorst, Adi Zitrin, Paola Dimauro, Yolanda Jimenez-Teja, Abdurro'uf, Mario Nonino, Ana Acebron, Felipe Andrade-Santos, Roberto J. Avila, Matthew B. Bayliss, Alex Benitez, Tom Broadhurst, Rachana Bhatawdekar, Marusa Bradac , et al. (38 additional authors not shown)

Abstract: The gravitationally lensed star WHL0137-LS, nicknamed Earendel, was identified with a photometric redshift $z_{phot} = 6.2 \pm 0.1$ based on images taken with the Hubble Space Telescope. Here we present James Webb Space Telescope (JWST) Near Infrared Camera (NIRCam) images of Earendel in 8 filters spanning 0.8--5.0$μみゅー$m. In these higher resolution images, Earendel remains a single unresolved point… ▽ More The gravitationally lensed star WHL0137-LS, nicknamed Earendel, was identified with a photometric redshift $z_{phot} = 6.2 \pm 0.1$ based on images taken with the Hubble Space Telescope. Here we present James Webb Space Telescope (JWST) Near Infrared Camera (NIRCam) images of Earendel in 8 filters spanning 0.8--5.0$μみゅー$m. In these higher resolution images, Earendel remains a single unresolved point source on the lensing critical curve, increasing the lower limit on the lensing magnification to $μみゅー> 4000$ and restricting the source plane radius further to $r < 0.02$ pc, or $\sim 4000$ AUえーゆー. These new observations strengthen the conclusion that Earendel is best explained by an individual star or multiple star system, and support the previous photometric redshift estimate. Fitting grids of stellar spectra to our photometry yields a stellar temperature of $T_{\mathrm{eff}} \simeq 13000$--16000 K assuming the light is dominated by a single star. The delensed bolometric luminosity in this case ranges from $\log(L) = 5.8$--6.6 $L_{\odot}$, which is in the range where one expects luminous blue variable stars. Follow-up observations, including JWST NIRSpec scheduled for late 2022, are needed to further unravel the nature of this object, which presents a unique opportunity to study massive stars in the first billion years of the universe. △ Less

Submitted 9 November, 2022; v1 submitted 18 August, 2022; originally announced August 2022.

Comments: Accepted to ApJL. Data products, lens models, and analysis code will be available online at https://cosmic-spring.github.io

The Present and Future of Astronomy (ASTRO2022)

Authors: Giacomo Beccari, Henri M. J. Boffin, Paola Andreani, Selma de Mink, Wendy Freedman, Michael Hill, Bruno Leibundgut, Federico Lelli, Anna Miotello, Sean Sapcariu

Abstract: Being one of the most fascinating and ancient sciences, astronomy has always played a special role in society. In 2022 ESO organised an online conference to offer the community a platform to discuss astronomical topics of sociological and philosophical relevance in a professional atmosphere. The talks touched on several crucial aspects, moving from the methodology of science to the use of metrics,… ▽ More Being one of the most fascinating and ancient sciences, astronomy has always played a special role in society. In 2022 ESO organised an online conference to offer the community a platform to discuss astronomical topics of sociological and philosophical relevance in a professional atmosphere. The talks touched on several crucial aspects, moving from the methodology of science to the use of metrics, to the importance of diversity in evaluation processes, and to the link between astronomy and society. △ Less

Submitted 8 August, 2022; originally announced August 2022.

Comments: Report on the ESO on-line conference ASTRO2022, published on the ESO Messenger 187, 2022

Horizons: Nuclear Astrophysics in the 2020s and Beyond

Authors: H. Schatz, A. D. Becerril Reyes, A. Best, E. F. Brown, K. Chatziioannou, K. A. Chipps, C. M. Deibel, R. Ezzeddine, D. K. Galloway, C. J. Hansen, F. Herwig, A. P. Ji, M. Lugaro, Z. Meisel, D. Norman, J. S. Read, L. F. Roberts, A. Spyrou, I. Tews, F. X. Timmes, C. Travaglio, N. Vassh, C. Abia, P. Adsley, S. Agarwal , et al. (140 additional authors not shown)

Abstract: Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilit… ▽ More Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 96 pages. Submitted to Journal of Physics G

Report number: LA-UR-22-23997

Probing Red Supergiant dynamics through photo-center displacements measured by Gaia

Authors: A. Chiavassa, R. Kudritzki, B. Davies, B. Freytag, S. E. de Mink

Abstract: Red supergiant (RSGs) are cool massive stars in a late phase of their evolution when the stellar envelope becomes fully convective. They are the brightest stars in the universe at infrared light and can be detected in galaxies far beyond the Local Group, allowing for accurate determination of chemical composition of galaxies. The study of their physical properties is extremely important for variou… ▽ More Red supergiant (RSGs) are cool massive stars in a late phase of their evolution when the stellar envelope becomes fully convective. They are the brightest stars in the universe at infrared light and can be detected in galaxies far beyond the Local Group, allowing for accurate determination of chemical composition of galaxies. The study of their physical properties is extremely important for various phenomena including the final fate of massive stars as type II supernovae and gravitational wave progenitors. We explore the well-studied nearby young stellar cluster chi Per. Using Gaia EDR3 data, we find the distance of the cluster (d = 2.260+-0.020 kpc). We then investigate the variability of the convection-related surface structure as a source for parallax measurement uncertainty. We use state-of-the-art 3D radiative hydrodynamics simulations with CO5BOLD and the post-processing radiative transfer code OPTIM3D to compute intensity maps in the Gaia G photometric system. We calculate the variabiltiy, as a function of time, of the intensity-weighted mean from the synthetic maps. We then select the RSG stars in the cluster and compare their uncertainty on parallaxes to the predictions of photocentre displacements. The synthetic maps of RSG show extremely irregular and temporal variable surfaces due to convection-related dynamics. Consequentially, the position of the photo-center varies during Gaia measurements between 0.033 and 0.130 AUえーゆー (up to 5% of the corresponding simulation stellar radius). We argue that the variability of the convection-related surface structures accounts for a substantial part of the Gaia EDR3 parallax error of the RSG sample. We suggest that the variation of the uncertainty on Gaia parallax could be exploited quantitatively using appropriate 3D simulations to extract, in a unique way, important information about the stellar dynamics and parameters of RSG stars. △ Less

Submitted 10 May, 2022; originally announced May 2022.

Comments: Press-release: https://www.mpa-garching.mpg.de/1054448/news20220530

Journal ref: Astronomy and Astrophysics, Volume 661, L1, May 2022

On the fate of the secondary white dwarf in double-degenerate double-detonation Type Ia supernovae

Authors: R. Pakmor, F. P. Callan, C. E. Collins, S. E. de Mink, A. Holas, W. E. Kerzendorf, M. Kromer, P. G. Neunteufel, John T. O'Brien, F. K. Roepke, A. J. Ruiter, I. R. Seitenzahl, Luke J. Shingles, S. A. Sim, S. Taubenberger

Abstract: The progenitor systems and explosion mechanism of Type Ia supernovae are still unknown. Currently favoured progenitors include double-degenerate systems consisting of two carbon-oxygen white dwarfs with thin helium shells. In the double-detonation scenario, violent accretion leads to a helium detonation on the more massive primary white dwarf that turns into a carbon detonation in its core and exp… ▽ More The progenitor systems and explosion mechanism of Type Ia supernovae are still unknown. Currently favoured progenitors include double-degenerate systems consisting of two carbon-oxygen white dwarfs with thin helium shells. In the double-detonation scenario, violent accretion leads to a helium detonation on the more massive primary white dwarf that turns into a carbon detonation in its core and explodes it. We investigate the fate of the secondary white dwarf, focusing on changes of the ejecta and observables of the explosion if the secondary explodes as well rather than survives. We simulate a binary system of a $1.05\,\mathrm{M_\odot}$ and a $0.7\,\mathrm{M_\odot}$ carbon-oxygen white dwarf with $0.03\,\mathrm{M_\odot}$ helium shells each. We follow the system self-consistently from inspiral to ignition, through the explosion, to synthetic observables. We confirm that the primary white dwarf explodes self-consistently. The helium detonation around the secondary white dwarf, however, fails to ignite a carbon detonation. We restart the simulation igniting the carbon detonation in the secondary white dwarf by hand and compare the ejecta and observables of both explosions. We find that the outer ejecta at $v~>~15\,000$\,km\,s$^{-1}$ are indistinguishable. Light curves and spectra are very similar until $\sim~40\,\mathrm{d}$ after explosion and the ejecta are much more spherical than violent merger models. The inner ejecta differ significantly slowing down the decline rate of the bolometric light curve after maximum of the model with a secondary explosion by $\sim20$ per cent. We expect future synthetic 3D nebular spectra to confirm or rule out either model. △ Less

Submitted 25 October, 2022; v1 submitted 28 March, 2022; originally announced March 2022.

Comments: 12 pages, 7 figures, accepted to MNRAS, comments welcome

Formation and fate of low metallicity stars in IllustrisTNG50

Authors: Ruediger Pakmor, Christine M. Simpson, Freeke van de Voort, Lars Hernquist, Lieke van Son, Martyna Chruślińska, Rebekka Bieri, Selma E. de Mink, Volker Springel

Abstract: Low metallicity stars give rise to unique spectacular transients and are of immense interest for understanding stellar evolution. Their importance has only grown further with the recent detections of mergers of stellar mass black holes that likely originate mainly from low metallicity progenitor systems. Moreover, the formation of low metallicity stars is intricately linked to galaxy evolution, in… ▽ More Low metallicity stars give rise to unique spectacular transients and are of immense interest for understanding stellar evolution. Their importance has only grown further with the recent detections of mergers of stellar mass black holes that likely originate mainly from low metallicity progenitor systems. Moreover, the formation of low metallicity stars is intricately linked to galaxy evolution, in particular to early enrichment and to later accretion and mixing of lower metallicity gas. Because low metallicity stars are difficult to observe directly, cosmological simulations are crucial for understanding their formation. Here we quantify the rates and locations of low metallicity star formation using the high-resolution TNG50 magnetohydrodynamical cosmological simulation, and we examine where low metallicity stars end up at $z=0$. We find that $20\%$ of stars with $Z_*<0.1\,\mathrm{Z_\odot}$ form after $z=2$, and that such stars are still forming in galaxies of all masses at $z=0$ today. Moreover, most low-metallicity stars at $z=0$ reside in massive galaxies. We analyse the radial distribution of low metallicity star formation, and discuss the curious case of seven galaxies in TNG50 that form stars from primordial gas even at $z=0$. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: 15 pages, 11 figures, accepted by MNRAS, comments welcome

Astrophysics with the Laser Interferometer Space Antenna

Authors: Pau Amaro Seoane, Jeff Andrews, Manuel Arca Sedda, Abbas Askar, Quentin Baghi, Razvan Balasov, Imre Bartos, Simone S. Bavera, Jillian Bellovary, Christopher P. L. Berry, Emanuele Berti, Stefano Bianchi, Laura Blecha, Stephane Blondin, Tamara Bogdanović, Samuel Boissier, Matteo Bonetti, Silvia Bonoli, Elisa Bortolas, Katelyn Breivik, Pedro R. Capelo, Laurentiu Caramete, Federico Cattorini, Maria Charisi, Sylvain Chaty , et al. (134 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery… ▽ More The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe. △ Less

Submitted 25 May, 2023; v1 submitted 11 March, 2022; originally announced March 2022.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

Properties of the Be-type stars in 30 Doradus

Authors: P. L. Dufton, D. J. Lennon, J. I. Villasenor, I. D. Howarth, C. J. Evans, S. E. de Mink, H. Sana, W. D. Taylor

Abstract: The evolutionary status of Be-type stars remains unclear, with both single-star and binary pathways having been proposed. Here, VFTS spectroscopy of 73 Be-type stars, in the spectral-type range, B0--B3, is analysed to estimate projected rotational velocities, radial velocities and stellar parameters. They are found to be rotating faster than the corresponding VFTS B-type sample but simulations imp… ▽ More The evolutionary status of Be-type stars remains unclear, with both single-star and binary pathways having been proposed. Here, VFTS spectroscopy of 73 Be-type stars, in the spectral-type range, B0--B3, is analysed to estimate projected rotational velocities, radial velocities and stellar parameters. They are found to be rotating faster than the corresponding VFTS B-type sample but simulations imply that their projected rotational velocities are inconsistent with them all rotating at near critical velocities. The de-convolution of the projected rotational velocities estimates leads to a mean rotational velocity estimate of 320-350 km/s, approximately 100 km/s larger than that for the corresponding B-type sample. There is a dearth of targets with rotational velocities less than 0.4 of the critical velocity, with a broad distribution reaching up to critical rotation. Our best estimate for the mean or median of the rotational velocitiy is 0.68 of the critical velocity. Rapidly-rotating B-type stars are more numerous than their Be-type counterparts, whilst the observed frequency of Be-type stars identified as binary systems is significantly lower than that for normal B-type stars, consistent with their respective radial-velocity dispersions. The semi-amplitudes for the Be-type binaries are also smaller. Similar results are found for a SMC Be-type sample centred on NGC346 with no significant differences being found between the two samples. These results are compared with the predictions of single and binary stellar evolutionary models for Be-type stars. Assuming that a single mechanism dominated the production of classical Be-type stars, our comparison would favour a binary evolutionary history. △ Less

Submitted 7 March, 2022; v1 submitted 4 March, 2022; originally announced March 2022.

Comments: Accepted for publication in MNRAS, Appendicesand Tables 2 and 4 published as supplementary material

The Candidate Progenitor Companion Star of the Type Ib/c SN 2013ge

Authors: Ori D. Fox, Schuyler D. Van Dyk, Benjamin F. Williams, Maria Drout, Emmanouil Zapartas, Nathan Smith, Dan Milisavljevic, Jennifer E. Andrews, K. Azalee Bostroem, Alexei V. Filippenko, Sebastian Gomez, Patrick L. Kelly, S. E. de Mink, Justin Pierel, Armin Rest, Stuart Ryder, Niharika Sravan, Lou Strolger, Qinan Wang, Kathryn E. Weil

Abstract: This Letter presents the detection of a source at the position of the Type Ib/c supernova (SN) 2013ge more than four years after the radioactive component is expected to have faded. This source could mark the first post-SN direct detection of a surviving companion to a stripped-envelope Type Ib/c explosion. We test this hypothesis and find the shape of the source's spectral energy distribution is… ▽ More This Letter presents the detection of a source at the position of the Type Ib/c supernova (SN) 2013ge more than four years after the radioactive component is expected to have faded. This source could mark the first post-SN direct detection of a surviving companion to a stripped-envelope Type Ib/c explosion. We test this hypothesis and find the shape of the source's spectral energy distribution is most consistent with that of a B5 I supergiant. While binary models tend to predict OB-type stars for stripped-envelope companions, the location of the source on a color-magnitude diagram (CMD) places it redward of its more likely position on the main sequence (MS). The source may be temporarily out of thermal equilibrium, or a cool and inflated non-MS companion, which is similar to the suggested companion of Type Ib SN 2019yvr that was constrained from pre-SN imaging. We also consider other possible physical scenarios for the source, including a fading SN, circumstellar shock interaction, line of site coincidence, and an unresolved host star cluster, all of which will require future observations to more definitively rule out. Ultimately, the fraction of surviving companions ("binary fraction") will provide necessary constraints on binary evolution models and the underlying physics. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Comments: Accepted to ApJL. 9 pages, 4 figures, 1 table

Stellar mergers as the origin of the blue main-sequence band in young star clusters

Authors: Chen Wang, Norbert Langer, Abel Schootemeijer, Antonino Milone, Ben Hastings, Xiao-Tian Xu, Julia Bodensteiner, Hugues Sana, Norberto Castro, D. J. Lennon, Pablo Marchant, A. de Koter, Selma E. de Mink

Abstract: Recent high-quality Hubble Space Telescope (HST) photometry shows that the main sequences (MS) stars of young star clusters form two discrete components in the color-magnitude diagram (CMD). Based on their distribution in the CMD, we show that stars of the blue MS component can be understood as slow rotators originating from stellar mergers. We derive the masses of the blue MS stars, and find that… ▽ More Recent high-quality Hubble Space Telescope (HST) photometry shows that the main sequences (MS) stars of young star clusters form two discrete components in the color-magnitude diagram (CMD). Based on their distribution in the CMD, we show that stars of the blue MS component can be understood as slow rotators originating from stellar mergers. We derive the masses of the blue MS stars, and find that they follow a nearly flat mass function, which supports their unusual formation path. Our results imply that the cluster stars gain their mass in two different ways, by disk accretion leading to rapid rotation, contributing to the red MS, or by binary merger leading to slow rotation and populating the blue MS. We also derive the approximate merger time of the individual stars of the blue MS component, and find a strong early peak in the merger rate, with a lower level merger activity prevailing for tens of Myr. This supports recent binary formation models, and explains new velocity dispersion measurements for members of young star clusters. Our findings shed new light on the origin of the bi-modal mass, spin, and magnetic field distributions of main-sequence stars. △ Less

Submitted 11 February, 2022; originally announced February 2022.

Comments: Published in Nature Astronomy online. Link to the paper:https://rdcu.be/cGLpH

Impact of Massive Binary Star and Cosmic Evolution on Gravitational Wave Observations II: Double Compact Object Rates and Properties

Authors: Floor S. Broekgaarden, Edo Berger, Simon Stevenson, Stephen Justham, Ilya Mandel, Martyna Chruślińska, Lieke A. C. van Son, Tom Wagg, Alejandro Vigna-Gómez, Selma E. de Mink, Debatri Chattopadhyay, Coenraad J. Neijssel

Abstract: Making the most of the rapidly increasing population of gravitational-wave detections of black hole (BH) and neutron star (NS) mergers requires comparing observations with population synthesis predictions. In this work we investigate the combined impact from the key uncertainties in population synthesis modelling of the isolated binary evolution channel: the physical processes in massive binary-st… ▽ More Making the most of the rapidly increasing population of gravitational-wave detections of black hole (BH) and neutron star (NS) mergers requires comparing observations with population synthesis predictions. In this work we investigate the combined impact from the key uncertainties in population synthesis modelling of the isolated binary evolution channel: the physical processes in massive binary-star evolution and the star formation history as a function of metallicity, $Z$, and redshift $z, \mathcal{S}(Z,z)$. Considering these uncertainties we create 560 different publicly available model realizations and calculate the rate and distribution characteristics of detectable BHBH, BHNS, and NSNS mergers. We find that our stellar evolution and $\mathcal{S}(Z,z)$ variations can impact the predicted intrinsic and detectable merger rates by factors $10^2$-$10^4$. We find that BHBH rates are dominantly impacted by $\mathcal{S}(Z,z)$ variations, NSNS rates by stellar evolution variations and BHNS rates by both. We then consider the combined impact from all uncertainties considered in this work on the detectable mass distribution shapes (chirp mass, individual masses and mass ratio). We find that the BHNS mass distributions are predominantly impacted by massive binary-star evolution changes. For BHBH and NSNS we find that both uncertainties are important. We also find that the shape of the delay time and birth metallicity distributions are typically dominated by the choice of $\mathcal{S}(Z,z)$ for BHBH, BHNS and NSNS. We identify several examples of robust features in the mass distributions predicted by all 560 models, such that we expect more than 95% of BHBH detections to contain a BH $\gtrsim 8\,\rm{M}_{\odot}$ and have mass ratios $\lesssim 4$. Our work demonstrates that it is essential to consider a wide range of allowed models to study double compact object merger rates and properties. △ Less

Submitted 14 December, 2021; v1 submitted 10 December, 2021; originally announced December 2021.

Comments: 25 pages of which 12 are (rainbow) figures. Main results: Fig 2, 4, 5 and 6. Code/results publicly available at https://github.com/FloorBroekgaarden/Double-Compact-Object-Mergers. Comments welcome

Gravitational wave sources in our Galactic backyard: Predictions for BHBH, BHNS and NSNS binaries detectable with LISA

Authors: Tom Wagg, Floor S. Broekgaarden, Selma E. de Mink, Lieke A. C. van Son, Neige Frankel, Stephen Justham

Abstract: Future searches for gravitational waves from space will be sensitive to double compact objects (DCOs) in our Milky Way. We present new simulations of the populations of double black holes (BHBHs), black hole neutron stars (BHNSs) and double neutron stars (NSNSs) that will be detectable by the planned space-based gravitational wave detector LISA. For our estimates, we use an empirically-informed mo… ▽ More Future searches for gravitational waves from space will be sensitive to double compact objects (DCOs) in our Milky Way. We present new simulations of the populations of double black holes (BHBHs), black hole neutron stars (BHNSs) and double neutron stars (NSNSs) that will be detectable by the planned space-based gravitational wave detector LISA. For our estimates, we use an empirically-informed model of the metallicity dependent star formation history of the Milky Way. We populate it using an extensive suite of binary population-synthesis predictions for varying assumptions relating to mass transfer, common-envelope, supernova kicks, remnant masses and wind mass loss physics. For a 4(10)-year LISA mission, we predict between 30-370(50-550) detections over these variations, out of which 6-154(9-238) are BHBHs, 2-198(3-289) are BHNSs and 3-35(4-57) are NSNSs. We discuss how the variations in the physics assumptions alter the distribution of properties of the detectable systems, even when the detection rates are unchanged. In particular we discuss the observable characteristics such as the chirp mass, eccentricity and sky localisation and how the BHBH, BHNS and NSNS populations can be distinguished, both from each other and from the more numerous double white dwarf population. We further discuss the possibility of multi-messenger observations of pulsar populations with the Square Kilometre Array (SKA) and assess the benefits of extending the LISA mission. △ Less

Submitted 26 November, 2021; originally announced November 2021.

Comments: Submitted. All code for reproducing results and figures is available https://github.com/TomWagg/detecting-DCOs-in-LISA. Data is available https://doi.org/10.5281/zenodo.4699712

LEGWORK: A python package for computing the evolution and detectability of stellar-origin gravitational-wave sources with space-based detectors

Authors: Tom Wagg, Katelyn Breivik, Selma E. de Mink

Abstract: We present LEGWORK (LISA Evolution and Gravitational Wave Orbit Kit), an open-source Python package for making predictions about stellar-origin gravitational wave sources and their detectability in LISA or other space-based gravitational wave detectors. LEGWORK can be used to evolve the orbits of sources due to gravitational wave emission, calculate gravitational wave strains (using post-Newtonian… ▽ More We present LEGWORK (LISA Evolution and Gravitational Wave Orbit Kit), an open-source Python package for making predictions about stellar-origin gravitational wave sources and their detectability in LISA or other space-based gravitational wave detectors. LEGWORK can be used to evolve the orbits of sources due to gravitational wave emission, calculate gravitational wave strains (using post-Newtonian approximations), compute signal-to-noise ratios and visualise the results. It can be applied to a variety of potential sources, including binaries consisting of white dwarfs, neutron stars and black holes. Although we focus on double compact objects, in principle LEGWORK can be used for any system with a user-specified orbital evolution, such as those affected by a third object or gas drag. We optimised the package to make it efficient for use in population studies which can contain tens-of-millions of sources. This paper describes the package and presents several potential use cases. We explain in detail the derivations of the expressions behind the package as well as identify and clarify some discrepancies currently present in the literature. We hope that LEGWORK will enable and accelerate future studies triggered by the rapidly growing interest in gravitational wave sources. △ Less

Submitted 4 August, 2022; v1 submitted 16 November, 2021; originally announced November 2021.

Comments: Published in ApJS. 17 pages, 8 figures. Install by running "pip install legwork", for documentation see https://legwork.readthedocs.io, feedback welcome

Journal ref: 2022, ApJS, 260, 52

Detailed models of interacting short-period massive binary stars

Authors: K. Sen, N. Langer, P. Marchant, A. Menon, S. E. de Mink, A. Schootemeijer, C. Schürmann, L. Mahy, B. Hastings, K. Nathaniel, H. Sana, C. Wang, X. T. Xu

Abstract: About a quarter of massive binary stars undergo mass transfer while both stars burn hydrogen at their cores, first on the thermal and then on the nuclear timescale. The nuclear timescale mass transfer leads to observable counterparts: the semi-detached so-called massive Algol binaries. However, comprehensive model predictions for these systems are sparse. We study them using a large grid of ~10,00… ▽ More About a quarter of massive binary stars undergo mass transfer while both stars burn hydrogen at their cores, first on the thermal and then on the nuclear timescale. The nuclear timescale mass transfer leads to observable counterparts: the semi-detached so-called massive Algol binaries. However, comprehensive model predictions for these systems are sparse. We study them using a large grid of ~10,000 detailed binary evolution models calculated with the stellar evolution code MESA, covering initial donor masses between 10-40 M$_{\odot}$ and initial orbital periods above 1.4 d, at a metallicity suitable for the Large Magellanic Cloud (LMC). Our models imply ~30, or ~3% of the ~1,000 core hydrogen burning O-star binaries in the LMC to be currently in the semi-detached phase. Our donor models are up to 25-times more luminous than single stars of identical mass and effective temperature, which agrees with the observed Algols. A comparison of our models with the observed orbital periods and mass ratios implies rather conservative mass transfer in some systems, while very inefficient one in others. This is generally well reproduced by our spin-dependent mass transfer algorithm, except for the lowest considered masses. The observations reflect the slow increase of the surface nitrogen enrichment of the donors during the semi-detached phase all the way to CNO equilibrium. We also investigate the properties of our models after core hydrogen depletion of the donor star, when these models correspond to Wolf-Rayet/helium+OB star binaries. A dedicated spectroscopic survey of massive Algol systems may allow to derive the dependence of the efficiency of thermal timescale mass transfer on the binary parameters, as well as the efficiency of semiconvective mixing in the stellar interior. This would be a crucial step towards reliable binary models up to the formation of supernovae and compact objects. △ Less

Submitted 9 December, 2021; v1 submitted 5 November, 2021; originally announced November 2021.

Comments: Accepted for publication in A&A (December 6, 2021)

Journal ref: A&A 659, A98 (2022)

The cosmic carbon footprint of massive stars stripped in binary systems

Authors: R. Farmer, E. Laplace, S. E. de Mink, S. Justham

Abstract: The cosmic origin of carbon, a fundamental building block of life, is still uncertain. Yield predictions for massive stars are almost exclusively based on single star models, even though a large fraction interact with a binary companion. Using the MESA stellar evolution code, we predict the carbon ejected in the winds and supernovae of single and binary-stripped stars at solar metallicity. We find… ▽ More The cosmic origin of carbon, a fundamental building block of life, is still uncertain. Yield predictions for massive stars are almost exclusively based on single star models, even though a large fraction interact with a binary companion. Using the MESA stellar evolution code, we predict the carbon ejected in the winds and supernovae of single and binary-stripped stars at solar metallicity. We find that binary-stripped stars are twice as efficient at producing carbon (1.5-2.6 times, depending on choices on the slope of the initial mass function and black hole formation). We confirm that this is because the convective helium core recedes in stars that have lost their hydrogen envelope, as noted previously. The shrinking of the core disconnects the outermost carbon-rich layers created during the early phase of helium burning from the more central burning regions. The same effect prevents carbon destruction, even when the supernova shock wave passes. The yields are sensitive to the treatment of mixing at convective boundaries, specifically during carbon-shell burning (variations up to 40%) and improving upon this should be a central priority for more reliable yield predictions. The yields are robust (variations less than 0.5%) across our range of explosion assumptions. Black hole formation assumptions are also important, implying that the stellar graveyard now explored by gravitational-wave detections may yield clues to better understand the cosmic carbon production. Our findings also highlight the importance of accounting for binary-stripped stars in chemical yield predictions and motivates further studies of other products of binary interactions. △ Less

Submitted 22 October, 2021; v1 submitted 8 October, 2021; originally announced October 2021.

Comments: 20 pages, 8 figures, 3 tables, Accepted ApJ

The redshift evolution of the binary black hole merger rate: a weighty matter

Authors: L. A. C. van Son, S. E. de Mink, T. Callister, S. Justham, M. Renzo, T. Wagg, F. S. Broekgaarden, F. Kummer, R. Pakmor, I. Mandel

Abstract: Gravitational wave detectors are starting to reveal the redshift evolution of the binary black hole (BBH) merger rate, $R_{\mathrm{BBH}}(z)$. We make predictions for $R_{\mathrm{BBH}}(z)$ as a function of black hole mass for systems originating from isolated binaries. To this end, we investigate correlations between the delay time and black hole mass by means of the suite of binary population synt… ▽ More Gravitational wave detectors are starting to reveal the redshift evolution of the binary black hole (BBH) merger rate, $R_{\mathrm{BBH}}(z)$. We make predictions for $R_{\mathrm{BBH}}(z)$ as a function of black hole mass for systems originating from isolated binaries. To this end, we investigate correlations between the delay time and black hole mass by means of the suite of binary population synthesis simulations, COMPAS. We distinguish two channels: the common envelope (CE), and the stable Roche-lobe overflow (RLOF) channel, characterised by whether the system has experienced a common envelope or not. We find that the CE channel preferentially produces BHs with masses below about $30\rm{M}_{\odot}$ and short delay times ($t_{\rm delay} \lesssim 1$Gyr), while the stable RLOF channel primarily forms systems with BH masses above $30\rm{M}_{\odot}$ and long delay times ($t_{\rm delay} \gtrsim 1$Gyr). We provide a new fit for the metallicity specific star-formation rate density based on the Illustris TNG simulations, and use this to convert the delay time distributions into a prediction of $R_{\mathrm{BBH}}(z)$. This leads to a distinct redshift evolution of $R_{\mathrm{BBH}}(z)$ for high and low primary BH masses. We furthermore find that, at high redshift, $R_{\mathrm{BBH}}(z)$ is dominated by the CE channel, while at low redshift it contains a large contribution ($\sim 40\%$) from the stable RLOF channel. Our results predict that, for increasing redshifts, BBHs with component masses above $30\rm{M}_{\odot}$ will become increasingly scarce relative to less massive BBH systems. Evidence of this distinct evolution of $R_{\mathrm{BBH}}(z)$ for different BH masses can be tested with future detectors. △ Less

Submitted 20 May, 2022; v1 submitted 4 October, 2021; originally announced October 2021.

Comments: 20 pages + appendices, data is available at https://doi.org/10.5281/zenodo.5544170

Journal ref: The Astrophysical Journal, Volume 931, (2022) 17

Rapid stellar and binary population synthesis with COMPAS

Authors: Team COMPAS, :, Jeff Riley, Poojan Agrawal, Jim W. Barrett, Kristan N. K. Boyett, Floor S. Broekgaarden, Debatri Chattopadhyay, Sebastian M. Gaebel, Fabian Gittins, Ryosuke Hirai, George Howitt, Stephen Justham, Lokesh Khandelwal, Floris Kummer, Mike Y. M. Lau, Ilya Mandel, Selma E. de Mink, Coenraad Neijssel, Tim Riley, Lieke van Son, Simon Stevenson, Alejandro Vigna-Gomez, Serena Vinciguerra, Tom Wagg , et al. (1 additional authors not shown)

Abstract: Compact Object Mergers: Population Astrophysics and Statistics (COMPAS; https://compas.science) is a public rapid binary population synthesis code. COMPAS generates populations of isolated stellar binaries under a set of parametrized assumptions in order to allow comparisons against observational data sets, such as those coming from gravitational-wave observations of merging compact remnants. It i… ▽ More Compact Object Mergers: Population Astrophysics and Statistics (COMPAS; https://compas.science) is a public rapid binary population synthesis code. COMPAS generates populations of isolated stellar binaries under a set of parametrized assumptions in order to allow comparisons against observational data sets, such as those coming from gravitational-wave observations of merging compact remnants. It includes a number of tools for population processing in addition to the core binary evolution components. COMPAS is publicly available via the github repository https://github.com/TeamCOMPAS/COMPAS/, and is designed to allow for flexible modifications as evolutionary models improve. This paper describes the methodology and implementation of COMPAS. It is a living document which will be updated as new features are added to COMPAS; the current document describes COMPAS v02.21.00. △ Less

Submitted 28 December, 2021; v1 submitted 20 September, 2021; originally announced September 2021.

Comments: Code publicly available via https://compas.science . Minor updates to match version accepted to ApJS

The B-type Binaries Characterisation Programme I. Orbital solutions for the 30 Doradus population

Authors: J. I. Villaseñor, W. D. Taylor, C. J. Evans, O. H. Ramírez-Agudelo, H. Sana, L. A. Almeida, S. E. de Mink, P. L. Dufton, N. Langer

Abstract: We present results from the B-type Binaries Characterisation (BBC) programme, a multi-epoch spectroscopic study of 88 early B-type binary candidates in the 30 Doradus region of the Large Magellanic Cloud (LMC). From radial-velocity analysis of 29 observational epochs we confirm the binary status of 64 of our targets, comprising 50 SB1 and 14 SB2 B-type binaries. A further 20 systems (classified as… ▽ More We present results from the B-type Binaries Characterisation (BBC) programme, a multi-epoch spectroscopic study of 88 early B-type binary candidates in the 30 Doradus region of the Large Magellanic Cloud (LMC). From radial-velocity analysis of 29 observational epochs we confirm the binary status of 64 of our targets, comprising 50 SB1 and 14 SB2 B-type binaries. A further 20 systems (classified as SB1*) show clear signs of periodicity but with more tentative periods. Orbital solutions are presented for these 84 systems, providing the largest homogeneous sample to date of the binary properties of early B-type stars. Our derived orbital-period distribution is generally similar to those for samples of more massive (O-type) binaries in both the LMC and the Galaxy. This similarity with the properties of the more massive O-type binaries is important as early B-type stars are expected to account for the majority of core-collapse supernovae. Differences in the period distributions of the different samples start to increase above 4 d, and are also present between the earliest (B0-0.7) and later-type (B1-2.5) systems within the BBC sample, although further study is required to understand if this is an observational bias or a real physical effect. We have examined the semi-amplitude velocities and orbital periods of our sample to identify potential candidates that could hide compact companions. Comparing with probability distributions of finding black hole companions to OB-type stars from a recent theoretical study, we have found 16 binaries in the higher probability region that warrant further study. △ Less

Submitted 21 July, 2021; originally announced July 2021.

Comments: 42 pages (main body: 28 pages, supplementary material: 14 pages), 20 figures. Accepted for publication in MNRAS