Showing 1–50 of 141 results for author: Eldridge, J

The First Billion Years, According to JWST

Authors: Angela Adamo, Hakim Atek, Micaela B. Bagley, Eduardo Bañados, Kirk S. S. Barrow, Danielle A. Berg, Rachel Bezanson, Maruša Bradač, Gabriel Brammer, Adam C. Carnall, John Chisholm, Dan Coe, Pratika Dayal, Daniel J. Eisenstein, Jan J. Eldridge, Andrea Ferrara, Seiji Fujimoto, Anna de Graaff, Melanie Habouzit, Taylor A. Hutchison, Jeyhan S. Kartaltepe, Susan A. Kassin, Mariska Kriek, Ivo Labbé, Roberto Maiolino , et al. (24 additional authors not shown)

Abstract: With stunning clarity, JWST has revealed the Universe's first billion years. The scientific community is analyzing a wealth of JWST imaging and spectroscopic data from that era, and is in the process of rewriting the astronomy textbooks. Here, 1.5 years into the JWST science mission, we provide a snapshot of the great progress made towards understanding the initial chapters of our cosmic history.… ▽ More With stunning clarity, JWST has revealed the Universe's first billion years. The scientific community is analyzing a wealth of JWST imaging and spectroscopic data from that era, and is in the process of rewriting the astronomy textbooks. Here, 1.5 years into the JWST science mission, we provide a snapshot of the great progress made towards understanding the initial chapters of our cosmic history. We highlight discoveries and breakthroughs, topics and issues that are not yet understood, and questions that will be addressed in the coming years, as JWST continues its revolutionary observations of the Early Universe. While this compendium is written by a small number of authors, invited to ISSI Bern in March 2024 as part of the 2024 ISSI Breakthrough Workshop, we acknowledge the work of a large community that is advancing our collective understanding of the evolution of the Early Universe. △ Less

Submitted 31 May, 2024; originally announced May 2024.

Comments: review article written by the attendees of the 2024 ISSI breakthrough workshop "The first billion year of the Universe", submitted. Comments welcome

Predicting gravitational wave signals from BPASS White Dwarf Binary and Black Hole Binary populations of a Milky Way-like galaxy model for LISA

Authors: Petra Tang, Jan Eldridge, Renate Meyer, Astrid Lamberts, Guillaume Boileau, Wouter van Zeist

Abstract: Galactic white dwarf binaries (WDBs) and black hole binaries (BHBs) will be gravitational wave (GW) sources for LISA. Their detection will provide insights into binary evolution and the evolution of our Galaxy through cosmic history. Here, we make predictions of the expected WDB and BHB population within our Galaxy. We combine predictions of the compact remnant binary populations expected by stell… ▽ More Galactic white dwarf binaries (WDBs) and black hole binaries (BHBs) will be gravitational wave (GW) sources for LISA. Their detection will provide insights into binary evolution and the evolution of our Galaxy through cosmic history. Here, we make predictions of the expected WDB and BHB population within our Galaxy. We combine predictions of the compact remnant binary populations expected by stellar evolution by using the detailed Binary Population and Spectral Synthesis code (BPASS) with a Milky Way analogue galaxy model from the Feedback In Realistic Environment (FIRE) simulations. We use \textsc{PhenomA} and \textsc{LEGWORK} to simulate LISA observations. Both packages make similar predictions that on average four Galactic BHBs and 673 Galactic WDBs above the signal-to-noise ratio (SNR) threshold of 7 after a four-year mission. We compare these predictions to earlier results using the Binary Star Evolution (BSE) code with the same FIRE model galaxy. We find that BPASS predicts a few more LISA observable Galactic BHBs and a twentieth of the Galactic WDBs. The differences are due to the different physical assumptions that have gone into the binary evolution calculations. These results indicate that the expected population of compact binaries that LISA will detect depends very sensitively on the binary population synthesis models used and thus observations of the LISA population will provide tight constraints on our modelling of binary stars. Finally, from our synthetic populations we have created mock LISA signals that can be used to test and refine data processing methods of the eventual LISA observations. △ Less

Submitted 30 May, 2024; originally announced May 2024.

The Pristine Inner Galaxy Survey (PIGS) IX. The largest detailed chemical analysis of very metal-poor stars in the Sagittarius dwarf galaxy

Authors: Federico Sestito, Sara Vitali, Paula Jofre, Kim A. Venn, David S. Aguado, Claudia Aguilera-Gómez, Anke Ardern-Arentsen, Danielle de Brito Silva, Raymond Carlberg, Camilla J. L. Eldridge, Felipe Gran, Vanessa Hill, Pascale Jablonka, Georges Kordopatis, Nicolas F. Martin, Tadafumi Matsuno, Samuel Rusterucci, Else Starkenburg, Akshara Viswanathan

Abstract: The most metal-poor stars provide valuable insights into the early chemical enrichment history of a system, carrying the chemical imprints of the first generations of supernovae. The most metal-poor region of the Sagittarius dwarf galaxy remains inadequately observed and characterised. To date, only $\sim4$ stars with [Fe/H]~$<-2.0$ have been chemically analysed with high-resolution spectroscopy.… ▽ More The most metal-poor stars provide valuable insights into the early chemical enrichment history of a system, carrying the chemical imprints of the first generations of supernovae. The most metal-poor region of the Sagittarius dwarf galaxy remains inadequately observed and characterised. To date, only $\sim4$ stars with [Fe/H]~$<-2.0$ have been chemically analysed with high-resolution spectroscopy. In this study, we present the most extensive chemical abundance analysis of 12 low-metallicity stars with metallicities down to [Fe/H]~$=-3.26$ and located in the main body of Sagittarius. These targets, selected from the Pristine Inner Galaxy Survey, were observed using the MIKE high-resolution spectrograph at the {\it Magellan-Clay} telescope, which allowed us to measure up to 17 chemical species. The chemical composition of these stars reflects the imprint of a variety of type~II supernovae (SNe~II). A combination of low- to intermediate-mass high-energy SNe and hypernovae ($\sim10-70\msun$) is required to account for the abundance patterns of the lighter elements up to the Fe-peak. The trend of the heavy elements suggests the involvement of compact binary merger events and fast-rotating (up to $\sim300\kms$) intermediate-mass to massive metal-poor stars ($\sim25-120\msun$) that are the sources of rapid and slow processes, respectively. Additionally, asymptotic giant branch stars contribute to a wide dispersion of [Ba/Mg] and [Ba/Eu]. The absence of an $α-$knee in our data indicates that type Ia supernovae did not contribute in the very metal-poor region ([Fe/H]~$\leq-2.0$). However, they might have started to pollute the interstellar medium at [Fe/H]~$>-2.0$, given the relatively low [Co/Fe] in this metallicity region. △ Less

Submitted 2 July, 2024; v1 submitted 30 April, 2024; originally announced May 2024.

Comments: Accepted for publication in A&A. New plot on [Eu/Mg]. Some refs updated

The effect of stellar rotation on black hole mass and spin

Authors: Sohan Ghodla, J. J. Eldridge

Abstract: The gravitational wave signature of a binary black hole (BBH) merger is dependent on its component mass and spin. If such black holes originate from rapidly rotating progenitors, the large angular momentum reserve in the star could drive a collapsar-like supernova explosion, hence substantially impacting these characteristics of the black holes in the binary. To examine the effect of stellar rotat… ▽ More The gravitational wave signature of a binary black hole (BBH) merger is dependent on its component mass and spin. If such black holes originate from rapidly rotating progenitors, the large angular momentum reserve in the star could drive a collapsar-like supernova explosion, hence substantially impacting these characteristics of the black holes in the binary. To examine the effect of stellar rotation on the resulting black hole mass and spin, we conduct a 1D general relativistic study of the end phase of the collapse. We find that the resulting black hole mass at times differs significantly from the previously assumed values. We quantify the dependence of the black hole spin magnitude on the hydrodynamics of the accretion flow, providing analytical relations for calculating the mass and spin based on the progenitor's pre-collapse properties. Depending on the nature of the accretion flow, our findings have implications for the black hole upper mass gap resulting from pair-instability supernovae, the maximum mass of a maximally rotating stellar black hole, and the maximum effective spin of a BBH formed in tidally locked helium star - black hole binary. △ Less

Submitted 24 June, 2024; v1 submitted 16 December, 2023; originally announced December 2023.

Comments: 20 pages, 22 figures

On the evolutionary history of a simulated disc galaxy as seen by phylogenetic trees

Authors: Danielle de Brito Silva, Paula Jofré, Patricia B. Tissera, Keaghan J. Yaxley, Jenny Gonzalez Jara, Camilla J. L. Eldridge, Emanuel Sillero, Robert M. Yates, Xia Hua, Payel Das, Claudia Aguilera-Gómez, Evelyn J. Johnston, Alvaro Rojas-Arriagada, Robert Foley, Gerard Gilmore

Abstract: Phylogenetic methods have long been used in biology, and more recently have been extended to other fields - for example, linguistics and technology - to study evolutionary histories. Galaxies also have an evolutionary history, and fall within this broad phylogenetic framework. Under the hypothesis that chemical abundances can be used as a proxy for interstellar medium's DNA, phylogenetic methods a… ▽ More Phylogenetic methods have long been used in biology, and more recently have been extended to other fields - for example, linguistics and technology - to study evolutionary histories. Galaxies also have an evolutionary history, and fall within this broad phylogenetic framework. Under the hypothesis that chemical abundances can be used as a proxy for interstellar medium's DNA, phylogenetic methods allow us to reconstruct hierarchical similarities and differences among stars - essentially a tree of evolutionary relationships and thus history. In this work, we apply phylogenetic methods to a simulated disc galaxy obtained with a chemo-dynamical code to test the approach. We found that at least 100 stellar particles are required to reliably portray the evolutionary history of a selected stellar population in this simulation, and that the overall evolutionary history is reliably preserved when the typical uncertainties in the chemical abundances are smaller than 0.08 dex. The results show that the shape of the trees are strongly affected by the age-metallicity relation, as well as the star formation history of the galaxy. We found that regions with low star formation rates produce shorter trees than regions with high star formation rates. Our analysis demonstrates that phylogenetic methods can shed light on the process of galaxy evolution. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: Accepted for publication in ApJ, October 12th 2023

Reaching for the stars -- JWST/NIRSpec spectroscopy of a lensed star candidate at $z=4.76$

Authors: Lukas J. Furtak, Ashish K. Meena, Erik Zackrisson, Adi Zitrin, Gabriel B. Brammer, Dan Coe, José M. Diego, Jan J. Eldridge, Yolanda Jiménez-Teja, Vasily Kokorev, Massimo Ricotti, Brian Welch, Rogier A. Windhorst, Abdurro'uf, Felipe Andrade-Santos, Rachana Bhatawdekar, Larry D. Bradley, Tom Broadhurst, Wenlei Chen, Christopher J. Conselice, Pratika Dayal, Brenda L. Frye, Seiji Fujimoto, Tiger Y. -Y. Hsiao, Patrick L. Kelly , et al. (9 additional authors not shown)

Abstract: We present JWST/NIRSpec observations of a highly magnified star candidate at a photometric redshift of $z_{\mathrm{phot}}\simeq4.8$, previously detected in JWST/NIRCam imaging of the strong lensing (SL) cluster MACS J0647+7015 ($z=0.591$). The spectroscopic observation allows us to precisely measure the redshift of the host arc at $z_{\mathrm{spec}}=4.758\pm0.004$, and the star's spectrum displays… ▽ More We present JWST/NIRSpec observations of a highly magnified star candidate at a photometric redshift of $z_{\mathrm{phot}}\simeq4.8$, previously detected in JWST/NIRCam imaging of the strong lensing (SL) cluster MACS J0647+7015 ($z=0.591$). The spectroscopic observation allows us to precisely measure the redshift of the host arc at $z_{\mathrm{spec}}=4.758\pm0.004$, and the star's spectrum displays clear Lyman- and Balmer-breaks commensurate with this redshift. A fit to the spectrum suggests a B-type super-giant star of surface temperature $T_{\mathrm{eff,B}}\simeq15000$ K with either a redder F-type companion ($T_{\mathrm{eff,F}}\simeq6250$K) or significant dust attenuation ($A_V\simeq0.82$) along the line of sight. We also investigate the possibility that this object is a magnified young globular cluster rather than a single star. We show that the spectrum is in principle consistent with a star cluster, which could also accommodate the lack of flux variability between the two epochs. However, the lack of a counter image and the strong upper limit on the size of the object from lensing symmetry, $r\lesssim0.5$ pc, could indicate that this scenario is somewhat less likely -- albeit not completely ruled out by the current data. The presented spectrum seen at a time when the Universe was only $\sim1.2$ Gyr old showcases the ability of JWST to study early stars through extreme lensing. △ Less

Submitted 25 September, 2023; v1 submitted 31 July, 2023; originally announced August 2023.

Comments: Accepted for publication in MNRAS letters. v2 updated to match the published version

Observational implications of cosmologically coupled black holes

Authors: Sohan Ghodla, Richard Easther, M. M. Briel, J. J. Eldridge

Abstract: It was recently suggested that "cosmologically coupled" black holes with masses that increase in proportion to the volume of the Universe might constitute the physical basis of dark energy. We take this claim at face value and discuss its potential astrophysical implications. We show that the gravitational wave emission in binary systems would be significantly enhanced so that the number of black… ▽ More It was recently suggested that "cosmologically coupled" black holes with masses that increase in proportion to the volume of the Universe might constitute the physical basis of dark energy. We take this claim at face value and discuss its potential astrophysical implications. We show that the gravitational wave emission in binary systems would be significantly enhanced so that the number of black hole mergers would exceed the observed rate by orders of magnitude, with typical masses much larger than those seen by the LIGO-Virgo-KAGRA network. Separately, if the mass growth happens at fixed angular momentum, the supermassive black holes in matter-deficient elliptical galaxies should be slowly rotating. Finally, cosmological coupling would stabilize small black holes against Hawking radiation-induced evaporation. △ Less

Submitted 20 July, 2023; v1 submitted 13 June, 2023; originally announced June 2023.

Comments: 11 Pages, 5 Figures, Matches Accepted version

Sustained super-Eddington accretion around neutron stars & black holes

Authors: Sohan Ghodla, J. J. Eldridge

Abstract: Recently, it was shown that the formation of a photon-trapping surface might not be sufficient to ensure unimpeded super-Eddington (SE) accretion. In light of this finding, here we derive a condition such that sustained and unimpeded SE accretion could be achieved in optically thick slim accretion disks surrounding neutron stars (NSs) and black holes (BHs). For this, we calculate a semi-analytic a… ▽ More Recently, it was shown that the formation of a photon-trapping surface might not be sufficient to ensure unimpeded super-Eddington (SE) accretion. In light of this finding, here we derive a condition such that sustained and unimpeded SE accretion could be achieved in optically thick slim accretion disks surrounding neutron stars (NSs) and black holes (BHs). For this, we calculate a semi-analytic approximation of the self-similar global radial velocity expression for an advection-dominated flow. Neglecting the influence of relativistic jets on the accretion flow, we find that for Eddington fraction $\dot{m} \gtrsim 1.5 (ε/0.1)^{3/5}$ ($ε$ being the accretion efficiency) sustained SE accretion might be possible in slim disks around BHs irrespective of their spin. The same condition holds for NSs when $ε> 0.03$. The presence of a surface magnetic field might truncate the disk at the magnetosphere of the NS, resulting in lower efficiencies and consequently changing the condition to $\dot{m} > 0.013 ε^{-19/31}$. Our approach suggests that sustained SE accretion might almost always be possible around NSs and BHs hosting accretion disks. △ Less

Submitted 27 May, 2023; v1 submitted 16 May, 2023; originally announced May 2023.

Comments: 7 pages, 4 Figures. Matches published version

JWST NIRSpec spectroscopy of the triply-lensed $z = 10.17$ galaxy MACS0647$-$JD

Authors: Tiger Yu-Yang Hsiao, Abdurro'uf, Dan Coe, Rebecca L. Larson, Intae Jung, Matilde Mingozzi, Pratika Dayal, Nimisha Kumari, Vasily Kokorev, Anton Vikaeus, Gabriel Brammer, Lukas J. Furtak, Angela Adamo, Felipe Andrade-Santos, Jacqueline Antwi-Danso, Marusa Bradac, Larry D. Bradley, Tom Broadhurst, Adam C. Carnall, Christopher J. Conselice, Jose M. Diego, Megan Donahue, Jan J. Eldridge, Seiji Fujimoto, Alaina Henry , et al. (16 additional authors not shown)

Abstract: We present JWST/NIRSpec prism spectroscopy of MACS0647-JD, the triply-lensed $z \sim 11$ candidate discovered in HST imaging and spatially resolved by JWST imaging into two components A and B. Spectroscopy of component A yields a spectroscopic redshift $z=10.17$ based on 7 detected emission lines: CIII] $λλ$1907,1909, [OII] $λ$3727, [NeIII] $λ$3869, [NeIII] $λ$3968, H$δ$ $λ$4101, H$γ$ $λ$4340, and… ▽ More We present JWST/NIRSpec prism spectroscopy of MACS0647-JD, the triply-lensed $z \sim 11$ candidate discovered in HST imaging and spatially resolved by JWST imaging into two components A and B. Spectroscopy of component A yields a spectroscopic redshift $z=10.17$ based on 7 detected emission lines: CIII] $λλ$1907,1909, [OII] $λ$3727, [NeIII] $λ$3869, [NeIII] $λ$3968, H$δ$ $λ$4101, H$γ$ $λ$4340, and [OIII] $λ$4363. These are the second-most distant detections of these emission lines to date, in a galaxy observed just 460 million years after the Big Bang. Based on observed and extrapolated line flux ratios we derive a gas-phase metallicity $Z =$ log(O/H) = $7.5 - 8.0$, or $(0.06 - 0.2)$ $Z_\odot$, ionization parameter log($U$) $\sim -1.9\pm0.2$, and an ionizing photon production efficiency ${\rm log}(ξ_{\rm ion})=25.2\pm0.2\,$erg$^{-1}$ Hz. The spectrum has a softened Lyman-$α$ break, evidence for a strong Ly$α$ damping wing, suggesting that MACS0647-JD was unable to ionize its surroundings beyond its immediate vicinity ($R_{\text{HII}} \ll 1$ pMpc). The Ly$α$ damping wing also suppresses the F150W photometry, explaining the slightly overestimated photometric redshift $z = 10.6 \pm 0.3$. MACS0647-JD has a stellar mass log($M/M_\odot$) = $8.1 \pm 0.3$, including $\sim$ 6$\times 10^7 M_\odot$ in component A, most of which formed recently (within $\sim$ 20 Myr) with a star formation rate $2\pm1 M_\odot$ / yr, all within an effective radius $70\pm24\,$pc. The smaller component B ($r \sim 20$) pc is likely older ($\sim$100 Myr) with more dust ($A_V \sim 0.1$ mag), as found previously. Spectroscopy of a fainter companion galaxy C separated by a distance of \about\ 3$\,$kpc reveals a Lyman break consistent with $z = 10.17$. MACS0647-JD is likely the most distant galaxy merger known. △ Less

Submitted 2 July, 2024; v1 submitted 4 May, 2023; originally announced May 2023.

Comments: 22 pages, 6 figures, accepted for publication in ApJ

Exploring the impact of IMF and binary parameter stochasticity with a binary population synthesis code

Authors: E. R. Stanway, J. J. Eldridge

Abstract: Low mass star formation regions are unlikely to fully populate their initial mass functions, leading to a deficit of massive stars. In binary stellar populations, the full range of binary separations and mass ratios will also be underpopulated. To explore the effects of stochastic sampling in the integrated light of stellar clusters, we calculate models at a broad range of cluster masses, from 10^… ▽ More Low mass star formation regions are unlikely to fully populate their initial mass functions, leading to a deficit of massive stars. In binary stellar populations, the full range of binary separations and mass ratios will also be underpopulated. To explore the effects of stochastic sampling in the integrated light of stellar clusters, we calculate models at a broad range of cluster masses, from 10^2 to 10^7 M_sun, using a binary stellar population synthesis code. For clusters with stellar masses less than 10^5 M_sun, observable quantities show substantial scatter and their mean properties reflect the expected deficit of massive stars. In common with previous work, we find that purely stochastic sampling of the initial mass function appears to underestimate the mass of the most massive star in known clusters. However, even with this constraint, the majority of clusters likely inject sufficient kinetic energy to clear their birth clusters of gas. For quantities which directly measure the impact of the most massive stars, such as N_{ion}, xi_{ion} and beta_{UV}, uncertainties due to stochastic sampling dominate over those from the IMF shape or distribution of binary parameters, while stochastic sampling has a negligible effect on the stellar continuum luminosity density. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: 15 pages and 8 page appendix. Accepted for publication in MNRAS

Searching for ejected supernova companions in the era of precise proper motion and radial velocity measurements

Authors: A. A. Chrimes, A. J. Levan, J. J. Eldridge, M. Fraser, N. Gaspari, P. J. Groot, J. D. Lyman, G. Nelemans, E. R. Stanway, K. Wiersema

Abstract: The majority of massive stars are born in binaries, and most unbind upon the first supernova. With precise proper motion surveys such as Gaia, it is possible to trace back the motion of stars in the vicinity of young remnants to search for ejected companions. Establishing the fraction of remnants with an ejected companion, and the photometric and kinematic properties of these stars, offers unique… ▽ More The majority of massive stars are born in binaries, and most unbind upon the first supernova. With precise proper motion surveys such as Gaia, it is possible to trace back the motion of stars in the vicinity of young remnants to search for ejected companions. Establishing the fraction of remnants with an ejected companion, and the photometric and kinematic properties of these stars, offers unique insight into supernova progenitor systems. In this paper, we employ binary population synthesis to produce kinematic and photometric predictions for ejected secondary stars. We demonstrate that the unbound neutron star velocity distribution from supernovae in binaries closely traces the input kicks. Therefore, the observed distribution of neutron star velocities should be representative of their natal kicks. We evaluate the probability for any given filter, magnitude limit, minimum measurable proper motion (as a function of magnitude), temporal baseline, distance and extinction that an unbound companion can be associated with a remnant. We compare our predictions with results from previous companion searches, and demonstrate that the current sample of stars ejected by the supernova of their companion can be increased by a factor of 5-10 with Gaia data release 3. Further progress in this area is achievable by leveraging the absolute astrometric precision of Gaia, and by obtaining multiple epochs of deep, high resolution near-infrared imaging with the Hubble Space Telescope, JWST and next-generation wide-field near-infrared observatories such as Euclid or the Nancy Grace Roman Space Telescope. △ Less

Submitted 5 April, 2023; originally announced April 2023.

Comments: Accepted for publication in MNRAS. 19 pages, 17 figures

The IACOB project VIII. Searching for empirical signatures of binarity in fast-rotating O-type stars

Authors: N. Britavskiy, S. Simón-Díaz, G. Holgado, S. Burssens, J. Maíz Apellániz, J. J. Eldridge, Y. Nazé, M. Pantaleoni González, A. Herrero

Abstract: The empirical distribution of projected rotational velocities (vsini) in massive O-type stars is characterized by a dominant slow velocity component and a tail of fast rotators. Binary interaction has been proposed to play a dominant role in the formation of this tail. We perform a complete and homogeneous search for empirical signatures of binarity in a sample of 54 fast-rotating stars with the a… ▽ More The empirical distribution of projected rotational velocities (vsini) in massive O-type stars is characterized by a dominant slow velocity component and a tail of fast rotators. Binary interaction has been proposed to play a dominant role in the formation of this tail. We perform a complete and homogeneous search for empirical signatures of binarity in a sample of 54 fast-rotating stars with the aim of evaluating this hypothesis. This working sample has been extracted from a larger sample of 415 Galactic O-type stars which covers the full range of vsini values. We use new and archival multi-epoch spectra in order to detect spectroscopic binary systems. We complement this information with Gaia proper motions and TESS photometric data to aid in the identification of runaway stars and eclipsing binaries, respectively. The identified fraction of single-lined spectroscopic binary (SB1) systems and apparently single stars among the fast-rotating sample is $\sim$18% and $\sim$70%, respectively. When comparing these percentages with those corresponding to the slow-rotating sample we find that our sample of fast rotators is characterized by a slightly larger percentage of SB1 systems ($\sim$18% vs. $\sim$13%) and a considerably smaller fraction of clearly detected SB2 systems (8% vs. 33%). Overall, there seems to be a clear deficit of spectroscopic binaries (SB1+SB2) among fast-rotating O-type stars ($\sim$26% vs. $\sim$46%). On the contrary, the fraction of runaway stars is significantly higher in the fast-rotating domain ($\sim$33-50%) than among those stars with vsini < 200 km/s. Lastly, almost 65% of the apparently single fast-rotating stars are runaways. Our empirical results seem to be in good agreement with the idea that the tail of fast-rotating O-type stars (with vsini > 200 km/s) is mostly populated by post-interaction binary products. △ Less

Submitted 15 March, 2023; v1 submitted 2 February, 2023; originally announced February 2023.

Comments: 33 pages, 16 figures, accepted for publication in "Astronomy and Astrophysics"

Journal ref: A&A 672, A22 (2023)

A high-mass X-ray binary descended from an ultra-stripped supernova

Authors: Noel D. Richardson, Clarissa Pavao, Jan J. Eldridge, Herbert Pablo, André-Nicolas Chené, Peter Wysocki, Douglas R. Gies, George Younes, Jeremy Hare

Abstract: Ultra-stripped supernovae are different from other terminal explosions of massive stars, as they show little or no ejecta from the actual supernova event. They are thought to occur in massive binary systems after the exploding star has lost its surface through interactions with its companion. Such supernovae produce little to no kick, leading to the formation of a neutron star without loss of the… ▽ More Ultra-stripped supernovae are different from other terminal explosions of massive stars, as they show little or no ejecta from the actual supernova event. They are thought to occur in massive binary systems after the exploding star has lost its surface through interactions with its companion. Such supernovae produce little to no kick, leading to the formation of a neutron star without loss of the binary companion, which itself may also evolve into another neutron star. Here we show that a recently discovered high-mass X-ray binary, CPD -29 2176 (CD -29 5159; SGR 0755-2933), has an evolutionary history that shows the neutron star component formed during an ultra-stripped supernova. The binary has orbital elements that are similar both in period and in eccentricity to one of 14 Be X-Ray binaries that have both known orbital periods and eccentricities. The identification of the progenitors systems for ultra-stripped supernovae is necessary as their evolution pathways leads to the formation of a binary neutron star systems. Binary neutron stars, such as the system that produced the kilonova GW170817 that was observed with both electromagnetic and gravitational energy, are known to produce a large quantity of heavy elements. △ Less

Submitted 31 January, 2023; originally announced February 2023.

Comments: published in Nature, February 1, 2023

Gravitational wave spectral synthesis

Authors: Wouter G. J. van Zeist, J. J. Eldridge, Petra N. Tang

Abstract: We study the LISA sources that arise from isolated binary evolution, and how these depend on age and metallicity, using model stellar populations from BPASS. We model these as single-aged populations which are analogous to star clusters. We calculate the combined GW spectrum of all the binaries within these model clusters, including all types of compact binaries as well as those with living stars.… ▽ More We study the LISA sources that arise from isolated binary evolution, and how these depend on age and metallicity, using model stellar populations from BPASS. We model these as single-aged populations which are analogous to star clusters. We calculate the combined GW spectrum of all the binaries within these model clusters, including all types of compact binaries as well as those with living stars. These results allow us to evaluate the detectability of star clusters with LISA. We find at late times the dominant sources are WD-WD binaries by factors of 50-200, but at times between $10^8$ and $10^9$ years we find a significant population of NS-WD and BH-WD binaries (2-40 per $10^6$ M$_{\odot}$), which is related to the treatment of mass transfer and common envelope events in BPASS, wherein mass transfer is relatively likely to be stable. Metallicity also has an effect on the GW spectrum and on the relative dominance of different types of binaries. Using the information about known star clusters will aid the identification of sky locations where one could expect LISA to find GW sources. △ Less

Submitted 2 August, 2023; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: 21 pages, 12 figures, 3 tables

Journal ref: MNRAS 524 (2023) 2836-2856

End-to-end study of the home and genealogy of the first binary neutron star merger

Authors: Heloise F. Stevance, Jan J. Eldridge, Elizabeth R. Stanway, Joe Lyman, Anna F. McLeod, Andrew J. Levan

Abstract: Binary neutron star mergers are one of the ultimate events of massive binary star evolution, and our understanding of their parent system is still in its infancy. Upcoming gravitational wave detections, coupled with multi-wavelength follow-up observations, will allow us to study an increasing number of these events by characterising their neighbouring stellar populations and searching for their pr… ▽ More Binary neutron star mergers are one of the ultimate events of massive binary star evolution, and our understanding of their parent system is still in its infancy. Upcoming gravitational wave detections, coupled with multi-wavelength follow-up observations, will allow us to study an increasing number of these events by characterising their neighbouring stellar populations and searching for their progenitors. Stellar evolution simulations are essential to this work but they are also based on numerous assumptions. Additionally, the models used to study the host galaxies differ from those used to characterise the progenitors and are typically based on single star populations. Here we introduce a framework to perform an end-to-end analysis and deploy it to the first binary neutron star merger - GW170817. With the Binary Population And Spectral Synthesis (BPASS) codes we are able to retrieve the physical properties of the host galaxy NGC 4993 as well as infer progenitor candidates. In our simulations there is a >98% chance that GW170817 originated from a stellar population with Z=0.010 born between 5 and 12.5 Gyrs ago. By carefully weighing the stellar genealogies we find that GW170817 most likely came from a binary system born with a 13-24 Msol primary and 10-12 Msol secondary which underwent two or three common envelope events over their lifetime. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Comments: Main text (3 Figures, 10 pages) and Supplementary Information (23 pages, 14 Figures). Published in Nature Astronomy

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

Can neutron star mergers alone explain the r-process enrichment of the Milky Way?

Authors: Chiaki Kobayashi, Ilya Mandel, Krzysztof Belczynski, Stephane Goriely, Thomas H. Janka, Oliver Just, Ashley J. Ruiter, Dany Van Beveren, Matthias U. Kruckow, Max M. Briel, Jan J. Eldridge, Elizabeth Stanway

Abstract: Comparing Galactic chemical evolution models to the observed elemental abundances in the Milky Way, we show that neutron star mergers can be a leading r-process site only if at low metallicities such mergers have very short delay times and significant ejecta masses that are facilitated by the masses of the compact objects. Namely, black hole-neutron star mergers, depending on the black-hole spins,… ▽ More Comparing Galactic chemical evolution models to the observed elemental abundances in the Milky Way, we show that neutron star mergers can be a leading r-process site only if at low metallicities such mergers have very short delay times and significant ejecta masses that are facilitated by the masses of the compact objects. Namely, black hole-neutron star mergers, depending on the black-hole spins, can play an important role in the early chemical enrichment of the Milky Way. We also show that none of the binary population synthesis models used in this paper, i.e., COMPAS, StarTrack, Brussels, ComBinE, and BPASS, can currently reproduce the elemental abundance observations. The predictions are problematic not only for neutron star mergers, but also for Type Ia supernovae, which may point to shortcomings in binary evolution models. △ Less

Submitted 18 December, 2022; v1 submitted 8 November, 2022; originally announced November 2022.

Comments: Accepted for publication in the Astrophysical Journal Letters, 15 pages, 4 figures, 1 table

JWST reveals a possible $z \sim 11$ galaxy merger in triply-lensed MACS0647$-$JD

Authors: Tiger Yu-Yang Hsiao, Dan Coe, Abdurro'uf, Lily Whitler, Intae Jung, Gourav Khullar, Ashish Kumar Meena, Pratika Dayal, Kirk S. S. Barrow, Lillian Santos-Olmsted, Adam Casselman, Eros Vanzella, Mario Nonino, Yolanda Jimenez-Teja, Masamune Oguri, Daniel P. Stark, Lukas J. Furtak, Adi Zitrin, Angela Adamo, Gabriel Brammer, Larry Bradley, Jose M. Diego, Erik Zackrisson, Steven L. Finkelstein, Rogier A. Windhorst , et al. (41 additional authors not shown)

Abstract: MACS0647$-$JD is a triply-lensed $z\sim11$ galaxy originally discovered with the Hubble Space Telescope. Here we report new JWST imaging, which clearly resolves MACS0647$-$JD as having two components that are either merging galaxies or stellar complexes within a single galaxy. Both are very small, with stellar masses $\sim10^8\,M_\odot$ and radii $r<100\,\rm pc$. The brighter larger component "A"… ▽ More MACS0647$-$JD is a triply-lensed $z\sim11$ galaxy originally discovered with the Hubble Space Telescope. Here we report new JWST imaging, which clearly resolves MACS0647$-$JD as having two components that are either merging galaxies or stellar complexes within a single galaxy. Both are very small, with stellar masses $\sim10^8\,M_\odot$ and radii $r<100\,\rm pc$. The brighter larger component "A" is intrinsically very blue ($β\sim-2.6$), likely due to very recent star formation and no dust, and is spatially extended with an effective radius $\sim70\,\rm pc$. The smaller component "B" appears redder ($β\sim-2$), likely because it is older ($100-200\,\rm Myr$) with mild dust extinction ($A_V\sim0.1\,\rm mag$), and a smaller radius $\sim20\,\rm pc$. We identify galaxies with similar colors in a high-redshift simulation, finding their star formation histories to be out of phase. With an estimated stellar mass ratio of roughly 2:1 and physical projected separation $\sim400\,\rm pc$, we may be witnessing a galaxy merger 400 million years after the Big Bang. We also identify a candidate companion galaxy C $\sim3\,{\rm kpc}$ away, likely destined to merge with galaxies A and B. The combined light from galaxies A+B is magnified by factors of $\sim$8, 5, and 2 in three lensed images JD1, 2, and 3 with F356W fluxes $\sim322$, $203$, $86\,\rm nJy$ (AB mag 25.1, 25.6, 26.6). MACS0647$-$JD is significantly brighter than other galaxies recently discovered at similar redshifts with JWST. Without magnification, it would have AB mag 27.3 ($M_{UV}=-20.4$). With a high confidence level, we obtain a photometric redshift of $z=10.6\pm0.3$ based on photometry measured in 6 NIRCam filters spanning $1-5\rmμm$, out to $4300\,Å$ rest-frame. JWST NIRSpec observations planned for January 2023 will deliver a spectroscopic redshift and a more detailed study of the physical properties of MACS0647$-$JD. △ Less

Submitted 31 May, 2023; v1 submitted 25 October, 2022; originally announced October 2022.

Comments: 27 pages, 14 figures, accepted for publication in ApJL

Nested Dust Shells around the Wolf-Rayet Binary WR 140 observed with JWST

Authors: Ryan M. Lau, Matthew J. Hankins, Yinuo Han, Ioannis Argyriou, Michael F. Corcoran, Jan J. Eldridge, Izumi Endo, Ori D. Fox, Macarena Garcia Marin, Theodore R. Gull, Olivia C. Jones, Kenji Hamaguchi, Astrid Lamberts, David R. Law, Thomas Madura, Sergey V. Marchenko, Hideo Matsuhara, Anthony F. J. Moffat, Mark R. Morris, Patrick W. Morris, Takashi Onaka, Michael E. Ressler, Noel D. Richardson, Christopher M. P. Russell, Joel Sanchez-Bermudez , et al. (7 additional authors not shown)

Abstract: Massive colliding-wind binaries that host a Wolf-Rayet (WR) star present a potentially important source of dust and chemical enrichment in the interstellar medium (ISM). However, the chemical composition and survival of dust formed from such systems is not well understood. The carbon-rich WR (WC) binary WR~140 presents an ideal astrophysical laboratory for investigating these questions given its w… ▽ More Massive colliding-wind binaries that host a Wolf-Rayet (WR) star present a potentially important source of dust and chemical enrichment in the interstellar medium (ISM). However, the chemical composition and survival of dust formed from such systems is not well understood. The carbon-rich WR (WC) binary WR~140 presents an ideal astrophysical laboratory for investigating these questions given its well-defined orbital period and predictable dust-formation episodes every 7.93 years around periastron passage. We present observations from our Early Release Science program (ERS1349) with the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) Medium-Resolution Spectrometer (MRS) and Imager that reveal the spectral and spatial signatures of nested circumstellar dust shells around WR~140. MIRI MRS spectroscopy of the second dust shell and Imager detections of over 17 shells formed throughout the past $\gtrsim130$ years confirm the survival of carbonaceous dust grains from WR~140 that are likely carriers of "unidentified infrared" (UIR)-band features at 6.4 and 7.7 $μ$m. The observations indicate that dust-forming WC binaries can enrich the ISM with organic compounds and carbonaceous dust. △ Less

Submitted 12 October, 2022; originally announced October 2022.

Comments: Published in Nature Astronomy on Oct 12, 2022; 21 pages, 5 figures, 2 tables

Journal ref: Lau, R.M., Hankins, M.J., Han, Y. et al. Nested dust shells around the Wolf-Rayet binary WR 140 observed with JWST. Nat Astron (2022)

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

Evaluating chemically homogeneous evolution in stellar binaries: Electromagnetic implications -- Ionizing photons, SLSN-I, GRB, Ic-BL

Authors: Sohan Ghodla, J. J. Eldridge, Elizabeth R. Stanway, Héloïse F. Stevance

Abstract: We investigate the occurrence of rapid-rotation induced chemically homogeneous evolution (CHE) due to strong tides and mass accretion in binaries. To this end, we generalize the relation in Packet (1981) to calculate the minimum angular momentum (AM) accretion required by a secondary star to experience accretion-induced CHE. Contrary to traditionally assumed 5-10 percent accretion of initial mass… ▽ More We investigate the occurrence of rapid-rotation induced chemically homogeneous evolution (CHE) due to strong tides and mass accretion in binaries. To this end, we generalize the relation in Packet (1981) to calculate the minimum angular momentum (AM) accretion required by a secondary star to experience accretion-induced CHE. Contrary to traditionally assumed 5-10 percent accretion of initial mass ($Z \lesssim 0.004$, $M \gtrsim$ 20 M$_{\odot}$) for spinning up the accretor (resulting in CHE) this value can drop to $\sim$ 2 percent for efficient AM accretion while for certain systems it could be substantially larger. We conduct a population study using \textsc{bpass} by evolving stars under the influence of strong tides in short-period binaries and also account for the updated effect of accretion-induced spin-up. We find accretion CHE (compared to tidal CHE) to be the dominant means of producing homogeneous stars even at 10 percent AM accretion efficiency during mass transfer. Unlike tidal CHE, it is seen that CH stars arising due to accretion can retain a larger fraction of their AM till core collapse. Thus we show that accretion CHE could be an important formation channel for energetic electromagnetic transients like GRBs, Ic-BL (SLSN-I, Ic-BL) under the collapsar (magnetar) formalism and a single CH star could lead to both the transients under their respective formation scenario. Lastly, we show that under the current treatment of CHE, the emission rate of ionizing photons by such stars decreases more rapidly at higher metallicities than previously predicted. △ Less

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

Comments: Contents unchanged from last version. Added few references. 19 Pages, 17 Figures

New constraints on the Bray conservation-of-momentum natal kick model from multiple distinct observations

Authors: S. M. Richards, J. J. Eldridge, M. M. Briel, H. F. Stevance, R. Willcox

Abstract: Natal supernova kicks, the linear momentum compact remnants receive during their formation, are an essential part of binary population synthesis (BPS) models. Although these kicks are well-supported by evidence, their underlying distributions and incorporation into BPS models is uncertain. In this work, we investigate the nature of natal kicks using a previously proposed analytical prescription wh… ▽ More Natal supernova kicks, the linear momentum compact remnants receive during their formation, are an essential part of binary population synthesis (BPS) models. Although these kicks are well-supported by evidence, their underlying distributions and incorporation into BPS models is uncertain. In this work, we investigate the nature of natal kicks using a previously proposed analytical prescription where the strength of the kick is given by $v_\text{k}=α\frac{m_\text{ejecta}}{m_\text{remnant}}+β~\text{km s}^{-1}$, for free parameters $α$ and $β$. We vary the free parameters over large ranges of possible values, comparing these synthetic populations simultaneously against four constraints: the merger rate of compact binary neutron star (BNS) systems, the period-eccentricity distribution of galactic BNSs, the velocity distribution of single-star pulsars, and the likelihood for low-ejecta mass supernovae to produce low-velocity kicks. We find that different samples of the parameter space satisfy each tests, and only 1 per cent of the models satisfy all four constraints simultaneously. Although we cannot identify a single best kick model, we report $α=115^{+40}_{-55}~\text{km s}^{-1}, β=15^{+10}_{-15}~\text{km s}^{-1}$ as the center of the region of the parameter space that fulfils all of our constraints, and expect $β\geq0~\text{km s}^{-1}$ as a further constraint. We also suggest further observations that will enable future refinement of the kick model. A sensitive test for the kick model will be the redshift evolution of the BNS merger rate since this is effectively a direct measure of the delay-time distribution for mergers. For our best fitting values, we find that the peak of the BNS merger rate is the present-day. △ Less

Submitted 29 March, 2023; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: Accepted by MNRAS

VFTS 243 as predicted by the BPASS fiducial models

Authors: H. F. Stevance, S. Ghodla, S. Richards, J. J. Eldridge, M. M. Briel, P. Tang

Abstract: The recent discovery of an unambiguous quiescent BH and main sequence O star companion in VFTS 243 opens the door to new constraints on theoretical stellar evolution and population models looking to reproduce the progenitors of black hole - black hole binaries. Here we show that the Binary Population and Spectral Synthesis fiducial models (BPASSv2.2.1) natively predict VFTS243-like systems: We fin… ▽ More The recent discovery of an unambiguous quiescent BH and main sequence O star companion in VFTS 243 opens the door to new constraints on theoretical stellar evolution and population models looking to reproduce the progenitors of black hole - black hole binaries. Here we show that the Binary Population and Spectral Synthesis fiducial models (BPASSv2.2.1) natively predict VFTS243-like systems: We find that VFTS 243 likely originated from a binary system in a about 15 day orbit with primary mass ranging from 40 to 50 M$_{\odot}$ and secondary star with initial mass 24-25M$_{\odot}$. Additionally we find that the death of the primary star must have resulted in a low energy explosion E < 10$^{50}$ ergs.With a uniform prior we find that the kick velocity of the new-born black hole was <10 km/s. The very low eccentricity reported for VFTS 243 and the subsequent conclusion by the authors that the SN kick must have been very small is in line with the peak in the posterior distribution between 0 and 5 km/s found from our numerical simulations performed with a uniform prior. Finally, the reduced Hobbs kick distribution commonly used in black hole population synthesis is strongly disfavoured. △ Less

Submitted 10 February, 2023; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: 4 figures, 1 table, Published to MNRAS

Understanding the high-mass binary black hole population from stable mass transfer and super-Eddington accretion in BPASS

Authors: M. M. Briel, H. F. Stevance, J. J. Eldridge

Abstract: With the remarkable success of the LVK consortium in detecting binary black hole mergers, it has become possible to use the population properties to constrain our understanding of the progenitor stars' evolution. The most striking features of the observed primary black hole mass distributions are the extended tail up to 100M$_\odot$ and an excess of masses at 35M$_\odot$. Currently, isolated binar… ▽ More With the remarkable success of the LVK consortium in detecting binary black hole mergers, it has become possible to use the population properties to constrain our understanding of the progenitor stars' evolution. The most striking features of the observed primary black hole mass distributions are the extended tail up to 100M$_\odot$ and an excess of masses at 35M$_\odot$. Currently, isolated binary population synthesis have difficulty explaining these features. Using the well-tested BPASS detailed stellar binary evolution models to determine mass transfer stability, accretion rates, and remnant masses, we postulate that stable mass transfer with super-Eddington accretion is responsible for the extended tail. Furthermore, that the excess is not due to pulsation-pair instability, as previously thought, but due to stable mass transfer. These systems are able to merge within the Hubble time due to more stable mass transfer with extreme mass ratios that allows the orbits to shrink sufficiently to allow for a merger. These finding are at odds with those from other population synthesis codes but in agreement with other recent studies using detailed binary evolution models. △ Less

Submitted 3 October, 2022; v1 submitted 28 June, 2022; originally announced June 2022.

Comments: Submitted to MNRAS, comments welcome. 22 pages, 18 figures, 9 pages supplementary material

Towards an understanding of long gamma-ray burst environments through circumstellar medium population synthesis predictions

Authors: A. A. Chrimes, B. P. Gompertz, D. A. Kann, A. J. van Marle, J. J. Eldridge, P. J. Groot, T. Laskar, A. J. Levan, M. Nicholl, E. R. Stanway, K. Wiersema

Abstract: The temporal and spectral evolution of gamma-ray burst (GRB) afterglows can be used to infer the density and density profile of the medium through which the shock is propagating. In long-duration (core-collapse) GRBs, the circumstellar medium (CSM) is expected to resemble a wind-blown bubble, with a termination shock separating the stellar wind and the interstellar medium (ISM). A long standing pr… ▽ More The temporal and spectral evolution of gamma-ray burst (GRB) afterglows can be used to infer the density and density profile of the medium through which the shock is propagating. In long-duration (core-collapse) GRBs, the circumstellar medium (CSM) is expected to resemble a wind-blown bubble, with a termination shock separating the stellar wind and the interstellar medium (ISM). A long standing problem is that flat density profiles, indicative of the ISM, are often found at lower radii than expected for a massive star progenitor. Furthermore, the presence of both wind-like environments at high radii and ISM-like environments at low radii remains a mystery. In this paper, we perform a 'CSM population synthesis' with long GRB progenitor stellar evolution models. Analytic results for the evolution of wind blown bubbles are adjusted through comparison with a grid of 2D hydrodynamical simulations. Predictions for the emission radii, ratio of ISM to wind-like environments, wind and ISM densities are compared with the largest sample of afterglow-derived parameters yet compiled, which we make available for the community. We find that high ISM densities around 1000/cm3 best reproduce observations. If long GRBs instead occur in typical ISM densities of approximately 1/cm3, then the discrepancy between theory and observations is shown to persist at a population level. We discuss possible explanations for the origin of variety in long GRB afterglows, and for the overall trend of CSM modelling to over-predict the termination shock radius. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: 22 pages, 15 figures, 4 tables. Accepted for publication in MNRAS

WO-Type Wolf-Rayet Stars: the Last Hurrah of Massive Star Evolution

Authors: Erin Aadland, Philip Massey, D. John Hillier, Nidia I. Morrell, Kathryn F. Neugent, J. J. Eldridge

Abstract: Are WO-type Wolf Rayet (WR) stars in the final stage of massive star evolution before core-collapse? Although WC- and WO-type WRs have very similar spectra, WOs show a much stronger O VI $λλ$3811,34 emission-line feature. This has usually been interpreted to mean that WOs are more oxygen rich than WCs, and thus further evolved. However, previous studies have failed to model this line, leaving the… ▽ More Are WO-type Wolf Rayet (WR) stars in the final stage of massive star evolution before core-collapse? Although WC- and WO-type WRs have very similar spectra, WOs show a much stronger O VI $λλ$3811,34 emission-line feature. This has usually been interpreted to mean that WOs are more oxygen rich than WCs, and thus further evolved. However, previous studies have failed to model this line, leaving the relative abundances uncertain, and the relationship between the two types unresolved. To answer this fundamental question, we modeled six WCs and two WOs in the LMC using UV, optical, and NIR spectra with the radiative transfer code CMFGEN in order to determine their physical properties. We find that WOs are not richer in oxygen; rather, the O VI feature is insensitive to the abundance. However, the WOs have a significantly higher carbon and lower helium content than the WCs, and hence are further evolved. A comparison of our results with single-star Geneva and binary BPASS evolutionary models show that while many properties match, there is more carbon and less oxygen in the WOs than either set of evolutionary model predicts. This discrepancy may be due to the large uncertainty in the $^{12}$C+$^4$He$\rightarrow^{16}$O nuclear reaction rate; we show that if the Kunz et al. rate is decreased by a factor of 25-50%, then there would be a good match with the observations. It would also help explain the LIGO/VIRGO detection of black holes whose masses are in the theoretical upper mass gap. △ Less

Submitted 8 April, 2022; originally announced April 2022.

Comments: 55 pages, 22 figures, Accepted into ApJ

The Dependence of Theoretical Synthetic Spectra on $α$-enhancement in Young, Binary Stellar Populations

Authors: C. M. Byrne, E. R. Stanway, J. J. Eldridge, L. McSwiney, O. T. Townsend

Abstract: The enhancement of $α$ elements such as oxygen is an important phase in the chemical evolution of the early Universe, with nebular material becoming enriched in these elements sooner than iron. Here we present models which incorporate stellar spectra with $α$-enhanced compositions, focusing on the impact on the integrated light of young stellar populations, including those with large binary star f… ▽ More The enhancement of $α$ elements such as oxygen is an important phase in the chemical evolution of the early Universe, with nebular material becoming enriched in these elements sooner than iron. Here we present models which incorporate stellar spectra with $α$-enhanced compositions, focusing on the impact on the integrated light of young stellar populations, including those with large binary star fractions using the Binary Populations and Spectral Synthesis (BPASS) framework, while using Solar-scaled stellar evolution models. We find that broad spectrum outputs such as production of ionising flux, the ultraviolet spectral slope and optical colours are only weakly affected by a change in [$α$/Fe]. A number of features such as ultraviolet line indices (e.g. at 1719 and 1853Å) and optical line indices (such as MgB) are sensitive to such changes in composition for a continuously star-forming population and a single starburst population respectively. We find that at ages of more than 1Gyr, $α$-enhanced stellar populations appear bluer than their Solar-scaled counterparts, and show expected sensitivity of optical line indices to composition, in agreement with previous work. The ultraviolet stellar absorption lines are relatively insensitive to subtleties in the abundances ratios, although with sufficient measurement precision, a combination of UV line indices may enable a simultaneous measurement of total metallicity mass fraction and [$α$/Fe] in young stellar populations. The output models are designated as BPASS v2.3 and made available to the community with the aim of assisting interpretation of observations of high-redshift galaxies with the James Webb Space Telescope. △ Less

Submitted 7 April, 2022; v1 submitted 24 March, 2022; originally announced March 2022.

Comments: Accepted for publication in MNRAS 13 pages, 10 figures. Data release of BPASS v2.3 is available at http://warwick.ac.uk/bpass and http://bpass.auckland.ac.nz V2: Updated to correct the order of panels in Figure 4

New Insights into the Evolution of Massive Stars and Their Effects on Our Understanding of Early Galaxies

Authors: Jan J. Eldridge, Elizabeth R. Stanway

Abstract: The observable characteristics and subsequent evolution of young stellar populations is dominated by their massive stars. As our understanding of those massive stars and the factors affecting their evolution improves, so our interpretation of distant, unresolved stellar systems can also advance. As observations increasingly probe the distant Universe, and the rare low metallicity starbursts nearby… ▽ More The observable characteristics and subsequent evolution of young stellar populations is dominated by their massive stars. As our understanding of those massive stars and the factors affecting their evolution improves, so our interpretation of distant, unresolved stellar systems can also advance. As observations increasingly probe the distant Universe, and the rare low metallicity starbursts nearby, so the opportunity arises for these two fields to complement one another, and lead to an improved conception of both stars and galaxies. Here we review the current state of the art in modelling of massive star dominated stellar populations, and discuss their applications and implications for interpreting the distant Universe. Our principle findings include: - Binary evolutionary pathways must be included to understand the stellar populations in early galaxies. - Observations constraining the extreme ultraviolet spectrum of early galaxies are showing that current models are incomplete. The best current guess is that some form of accretion onto compact remnants is required. - The evolution and fates of very massive stars, of the order of 100Msun and above, may be key to fully understand aspects of early galaxies. △ Less

Submitted 3 February, 2022; originally announced February 2022.

Comments: Manuscript submitted to Annual Reviews of Astronomy and Astrophysics for future publication in Volume 60. This is the authors' own version. The final version will become available from https://www.annualreviews.org/journal/astro There are 43 pages, 5 figures and 2 tables

Explaining the differences in massive star models from various simulations

Authors: Poojan Agrawal, Dorottya Szécsi, Simon Stevenson, Jan J. Eldridge, Jarrod Hurley

Abstract: The evolution of massive stars is the basis of several astrophysical investigations, from predicting gravitational-wave event rates to studying star-formation and stellar populations in clusters. However, uncertainties in massive star evolution present a significant challenge when accounting for these models' behaviour in stellar population studies. In this work, we present a comparison between fi… ▽ More The evolution of massive stars is the basis of several astrophysical investigations, from predicting gravitational-wave event rates to studying star-formation and stellar populations in clusters. However, uncertainties in massive star evolution present a significant challenge when accounting for these models' behaviour in stellar population studies. In this work, we present a comparison between five published sets of stellar models from the BPASS, BoOST, Geneva, MIST, and PARSEC simulations at near-solar metallicity. The different sets of stellar models have been computed using slightly different physical inputs in terms of mass-loss rates and internal mixing properties. Moreover, these models also employ various pragmatic methods to overcome the numerical difficulties that arise due to the presence of density inversions in the outer layers of stars more massive than 40 M$_\odot$. These density inversions result from the combination of inefficient convection in the low-density envelopes of massive stars and the excess of radiative luminosity to the Eddington luminosity. We find that the ionizing radiation released by the stellar populations can change by up to 18 percent, the maximum radial expansion of a star can differ between 100-1600 R$_\odot$, and the mass of the stellar remnant can vary up to 20 M$_\odot$ between the five sets of simulations. We conclude that any attempts to explain observations that rely on the use of models of stars more massive than 40 M$_\odot$ should be made with caution. △ Less

Submitted 31 March, 2022; v1 submitted 6 December, 2021; originally announced December 2021.

Comments: added core mass information for the Geneva models, accepted for publication in MNRAS

To be or not to be a black hole: detailed binary population models as a sanity check

Authors: H. F. Stevance, S. G. Parsons, J. J. Eldridge

Abstract: We use the self-consistent stellar populations in the Binary Population A Spectral Synthesis (BPASS) models to assess whether NGC1850-BH1 is a black hole. Using search criteria based on reported physical properties in the literature we purposefully search for suitable systems with a black hole (or compact object) companion: we do not find any. Good matches to the observations are found in models w… ▽ More We use the self-consistent stellar populations in the Binary Population A Spectral Synthesis (BPASS) models to assess whether NGC1850-BH1 is a black hole. Using search criteria based on reported physical properties in the literature we purposefully search for suitable systems with a black hole (or compact object) companion: we do not find any. Good matches to the observations are found in models where the bright component is a stripped star and the companion is natively (meaning we did not impose this in our search) 1 to 2.3 magnitudes fainter than the primary in the optical bands. This alone can explain the lack of detection of the companion in the MUSE spectra without the need to invoke rapid rotation, although the conservative mass transfer exhibited by these particular models is likely to lead to rapidly rotating companions which could further smear their spectroscopic signatures. We advise that future claims of unseen black holes in binary systems would benefit from exploring detailed binary evolution models of stellar populations as a sanity check △ Less

Submitted 22 January, 2022; v1 submitted 30 November, 2021; originally announced December 2021.

Comments: 4 pages, 3 tables, 1 figure, published in MNRAS Letters - analysis code available at: https://github.com/UoA-Stars-And-Supernovae/ngc1850bh1

Estimating Transient Rates from Cosmological Simulations and BPASS

Authors: M. M. Briel, J. J. Eldridge, E. R. Stanway, H. F. Stevance, A. A. Chrimes

Abstract: The detection rate of electromagnetic (EM) and gravitational wave (GW) transients is growing exponentially. As the accuracy of the transient rates will significantly improve over the coming decades, so will our understanding of their evolution through cosmic history. To this end, we present predicted rates for EM and GW transients over the age of the Universe using Binary Population and Spectral S… ▽ More The detection rate of electromagnetic (EM) and gravitational wave (GW) transients is growing exponentially. As the accuracy of the transient rates will significantly improve over the coming decades, so will our understanding of their evolution through cosmic history. To this end, we present predicted rates for EM and GW transients over the age of the Universe using Binary Population and Spectral Synthesis (BPASS) results combined with four cosmic star formation histories (SFH). These include a widely used empirical SFH of Madau & Dickinson and those from three cosmological simulations: MilliMillennium, EAGLE and IllustrisTNG. We find that the choice of SFH significantly changes our predictions: transients with short delay times are most affected by the star formation rate, while long-delay time events tend to depend on the metallicity evolution of star formation. Importantly we find that the cosmological simulations have very different metallicity evolution that cannot be reproduced by the widely used metallicity model of Langer & Norman, which impacts the binary black hole merger and stripped-envelope supernovae rates in the local Universe most acutely. We recommend against using simple prescriptions for the metallicity evolution of the Universe when predicting the rates of events that can have long delay times and that are sensitive to metallicity evolution. △ Less

Submitted 15 November, 2021; originally announced November 2021.

Comments: Submitted to MNRAS, 19 pages 15 Figures, 4 Tables + 2 pages, 4 Tables in Appendix (Comments welcome)

Numerical experiments to help understand cause and effect in massive star evolution

Authors: Eoin Farrell, Jose Groh, Georges Meynet, JJ Eldridge

Abstract: The evolution of massive stars is affected by a variety of physical processes including convection, rotation, mass loss and binary interaction. Because these processes modify the internal chemical abundance profiles in multiple ways simultaneously, it can be challenging to determine which properties of the stellar interior are primarily driving the overall evolution. Building on previous work, we… ▽ More The evolution of massive stars is affected by a variety of physical processes including convection, rotation, mass loss and binary interaction. Because these processes modify the internal chemical abundance profiles in multiple ways simultaneously, it can be challenging to determine which properties of the stellar interior are primarily driving the overall evolution. Building on previous work, we develop a new modelling approach called SNAPSHOT that allows us to isolate the key features of the internal abundance profile that drive the evolution of massive stars. Using our approach, we compute numerical stellar structure models in thermal equilibrium covering key phases of stellar evolution. For the main sequence, we demonstrate that models with the same mass and very similar surface properties can have different internal distributions of hydrogen and convective core masses. We discuss why massive stars expand after the main sequence and the fundamental reasons for why they become red, blue or yellow supergiants. For the post-main sequence, we demonstrate that small changes in the abundance profile can cause very large effects on the surface properties. We also discuss the effects that produce blue supergiants and the cause of blue loops. Our models show that massive stars with lower metallicity tend to be more compact due to the combined effect of lower CNO abundances in the burning regions and lower opacity in the envelope. △ Less

Submitted 17 April, 2023; v1 submitted 6 September, 2021; originally announced September 2021.

Comments: Accepted in MNRAS

Binary evolution pathways of Blue Large-Amplitude Pulsators

Authors: C. M. Byrne, E. R. Stanway, J. J. Eldridge

Abstract: Blue Large-Amplitude Pulsators (BLAPs) are a recently discovered class of pulsating star, believed to be proto-white dwarfs, produced by mass stripping of a red giant when it has a small helium core. An outstanding question is why the stars in this class of pulsator seem to form two distinct groups by surface gravity, despite predictions that stars in the gap between them should also pulsate. We u… ▽ More Blue Large-Amplitude Pulsators (BLAPs) are a recently discovered class of pulsating star, believed to be proto-white dwarfs, produced by mass stripping of a red giant when it has a small helium core. An outstanding question is why the stars in this class of pulsator seem to form two distinct groups by surface gravity, despite predictions that stars in the gap between them should also pulsate. We use a binary population synthesis model to identify potential evolutionary pathways that a star can take to become a BLAP. We find that BLAPs can be produced either through common envelope evolution or Roche lobe overflow, with a Main Sequence star or an evolved compact object being responsible for the envelope stripping. The mass distribution of the inferred population indicates that fewer stars would be expected in the range of masses intermediate to the two known groups of pulsators, suggesting that the lack of observational discoveries in this region may be a result of the underlying population of pre-white dwarf stars. We also consider metallicity variation and find evidence that BLAPs at $Z = 0.010$ (half-Solar) would be pulsationally unstable and may also be more common. Based on this analysis, we expect the Milky Way to host around 12000 BLAPs and we predict the number density of sources expected in future observations such as the Legacy Survey of Space and Time at the Vera Rubin Observatory. △ Less

Submitted 30 July, 2021; originally announced July 2021.

Comments: Accepted for publication in MNRAS

Testing Evolutionary Models with Red Supergiant and Wolf-Rayet Populations

Authors: Philip Massey, Kathryn F. Neugent, Trevor Z. Dorn-Wallenstein, J. J. Eldridge, E. R. Stanway, Emily M. Levesque

Abstract: Despite the many successes that modern massive star evolutionary theory has enjoyed, reproducing the apparent trend in the relative number of red supergiants (RSGs) and Wolf-Rayet (WR) stars has remained elusive. Previous estimates show the RSG/WR ratio decreasing strongly with increasing metallicity. However, the evolutionary models have always predicted a relatively flat distribution for the RSG… ▽ More Despite the many successes that modern massive star evolutionary theory has enjoyed, reproducing the apparent trend in the relative number of red supergiants (RSGs) and Wolf-Rayet (WR) stars has remained elusive. Previous estimates show the RSG/WR ratio decreasing strongly with increasing metallicity. However, the evolutionary models have always predicted a relatively flat distribution for the RSG/WR ratio. In this paper we reexamine this issue, drawing on recent surveys for RSGs and WRs in the Magellanic Clouds, M31, and M33. The RSG surveys have used Gaia astrometry to eliminate foreground contamination, and have separated RSGs from asymptotic giant branch stars using near-infrared colors. The surveys for WRs have utilized interference filter imaging, photometry, and image subtraction techniques to identify candidates, which have then been confirmed spectroscopically. After carefully matching the observational criteria to the models, we now find good agreement in both the single-star Geneva and binary BPASS models with the new observations. The agreement is better when we shift the RSG effective temperatures derived from J-Ks photometry downwards by 200 K in order to agree with the Levesque TiO effective temperature scale. In an appendix we also present a source list of RSGs for the SMC which includes effective temperatures and luminosities derived from near-infrared 2MASS photometry, in the same manner as used for the other galaxies. △ Less

Submitted 17 July, 2021; originally announced July 2021.

Comments: Astrophysical Journal, in press

Comparing Compact Object Distributions from Mass- and Presupernova Core Structure-based Prescriptions

Authors: Rachel A. Patton, Tuguldur Sukhbold, J. J. Eldridge

Abstract: Binary population synthesis (BPS) employs prescriptions to predict final fates, explosion or implosion, and remnant masses based on one or two stellar parameters at the evolutionary cutoff imposed by the code, usually at or near central carbon ignition. In doing this, BPS disregards the integral role late-stage evolution plays in determining the final fate, remnant type, and remnant mass within th… ▽ More Binary population synthesis (BPS) employs prescriptions to predict final fates, explosion or implosion, and remnant masses based on one or two stellar parameters at the evolutionary cutoff imposed by the code, usually at or near central carbon ignition. In doing this, BPS disregards the integral role late-stage evolution plays in determining the final fate, remnant type, and remnant mass within the neutrino-driven explosion paradigm. To highlight differences between a popular prescription which relies only on the core and final stellar mass and emerging methods which rely on a star's presupernova core structure, we generate a series of compact object distributions using three different methods for a sample population of single and binary stars computed in BPASS. The first method estimates remnant mass based on a star's carbon-oxygen (CO) core mass and final total mass. The second method uses the presupernova core structure based on the bare CO-core models of \citet{Pat20} combined with a parameterized explosion criterion to first determine final fate and remnant type, then remnant mass. The third method associates presupernova helium-core masses with remnant masses determined from public explosion models which rely implicitly on core structure. We find that the core-/final mass-based prescription favors lower mass remnants, including a large population of mass gap black holes, and predicts neutron star masses which span a wide range, whereas the structure-based prescriptions favor slightly higher mass remnants, mass gap black holes only as low as 3.5 \Msun, and predict neutron star mass distributions which cluster in a narrow range. △ Less

Submitted 10 June, 2021; originally announced June 2021.

Comments: 11 pages, 9 figures Submitted to MNRAS

Forward Modelling the O3(a+b) GW transient mass distributions with BPASS by varying compact remnant mass and SNe kick prescriptions

Authors: Sohan Ghodla, Wouter G. J. van Zeist, J. J. Eldridge, Héloïse F. Stevance, Elizabeth R. Stanway

Abstract: We present forward modeling from the BPASS code suite of the population of observed gravitational wave (GW) transients reported by the LIGO/VIRGO consortium (LVC) during their third observing run, O3(a+b). Specifically, we predict the expected chirp mass and mass ratio distributions for GW transients, taking account of detector sensitivity to determine how many events should have been detected by… ▽ More We present forward modeling from the BPASS code suite of the population of observed gravitational wave (GW) transients reported by the LIGO/VIRGO consortium (LVC) during their third observing run, O3(a+b). Specifically, we predict the expected chirp mass and mass ratio distributions for GW transients, taking account of detector sensitivity to determine how many events should have been detected by the current detector network in O3(a+b). We investigate how these predictions change by alternating between four different remnant mass estimation schemes and two supernovae (SNe) kick prescriptions. We find that none of the model populations resulting from these variations accurately match the whole O3(a+b) GW transient catalog. However, agreement from some models to part of the catalog suggests ways to achieve a more complete fit. These include reducing the number of low mass black holes (BHs) close to the mass gap, while also increasing the number of higher mass BHs below the pair-instability SN limit. Finally, we find that the interaction between the value of the remnant mass from a stellar model and the choice of SN kick is complex and different kick prescriptions may be required depending on whether a neutron star or BH is formed. △ Less

Submitted 13 January, 2022; v1 submitted 10 May, 2021; originally announced May 2021.

Comments: 9 pages, 9 figures, 1 Table, Accepted by MNRAS

Binary pathways to SLSNe-I: SN 2017gci

Authors: H. F. Stevance, J. J. Eldridge

Abstract: Some hydrogen poor superluminous supernovae (SLSNe) exhibit bumps in the tails of their light-curves associated with hydrogen features in their late time spectra. Here we use the explosion parameters of one such SLSN -- SN 2017gci -- to search the stellar models of the Binary Population And Spectral Synthesis (BPASS) code for potential progenitors. We find good matches for a 30 $M_{\odot}$ progeni… ▽ More Some hydrogen poor superluminous supernovae (SLSNe) exhibit bumps in the tails of their light-curves associated with hydrogen features in their late time spectra. Here we use the explosion parameters of one such SLSN -- SN 2017gci -- to search the stellar models of the Binary Population And Spectral Synthesis (BPASS) code for potential progenitors. We find good matches for a 30 $M_{\odot}$ progenitor star in a binary system and no matches from single star models. Common envelope and mass transfer after the giant branch, combined with increased mass loss from strong stellar winds immediately before death, allow the progenitor to lose its hydrogen envelope decades before the explosion. This results in a hydrogen poor SLSN and allows for delayed interaction of the ejecta with the lost stellar material △ Less

Submitted 17 May, 2021; v1 submitted 7 April, 2021; originally announced April 2021.

Comments: Accepted MNRAS Letter - 5 pages, 2 figures, 1 table

Riroriro: Simulating gravitational waves and evaluating their detectability in Python

Authors: Wouter G. J. van Zeist, Héloïse F. Stevance, J. J. Eldridge

Abstract: Riroriro is a Python package to simulate the gravitational waveforms of binary mergers of black holes and/or neutron stars, and calculate several properties of these mergers and waveforms, specifically relating to their observability by gravitational wave detectors. The gravitational waveform simulation of Riroriro is based upon the methods of Buskirk and Babiuc-Hamilton (2019), a paper which de… ▽ More Riroriro is a Python package to simulate the gravitational waveforms of binary mergers of black holes and/or neutron stars, and calculate several properties of these mergers and waveforms, specifically relating to their observability by gravitational wave detectors. The gravitational waveform simulation of Riroriro is based upon the methods of Buskirk and Babiuc-Hamilton (2019), a paper which describes a computational implementation of an earlier theoretical gravitational waveform model by Huerta et al. (2017), using post-Newtonian expansions and an approximation called the implicit rotating source to simplify the Einstein field equations and simulate gravitational waves. Riroriro's calculation of signal-to-noise ratios (SNR) of gravitational wave events is based on the methods of Barrett et al. (2018), with the simpler gravitational wave model Findchirp (Allen et al. (2012)) being used for comparison and calibration in these calculations. △ Less

Submitted 11 March, 2021; originally announced March 2021.

Comments: 3 pages

Journal ref: Journal of Open Source Software (2021), 6(59), 2968

Progenitor, environment, and modelling of the interacting transient, AT 2016jbu (Gaia16cfr)

Authors: S. J. Brennan, M. Fraser, J. Johansson, A. Pastorello, R. Kotak, H. F. Stevance, T. -W. Chen, J. J. Eldridge, S. Bose, P. J. Brown, E. Callis, R. Cartier, M. Dennefeld, Subo Dong, P. Duffy, N. Elias-Rosa, G. Hosseinzadeh, E. Hsiao, H. Kuncarayakti, A. Martin-Carrillo, B. Monard, G. Pignata, D. Sand, B. J. Shappee, S. J. Smartt , et al. (45 additional authors not shown)

Abstract: We present the bolometric lightcurve, identification and analysis of the progenitor candidate, and preliminary modelling of AT2016jbu (Gaia16cfr). We find a progenitor consistent with a $\sim$22--25~$M_{\odot}$ yellow hypergiant surrounded by a dusty circumstellar shell, in agreement with what has been previously reported. We see evidence for significant photometric variability in the progenitor,… ▽ More We present the bolometric lightcurve, identification and analysis of the progenitor candidate, and preliminary modelling of AT2016jbu (Gaia16cfr). We find a progenitor consistent with a $\sim$22--25~$M_{\odot}$ yellow hypergiant surrounded by a dusty circumstellar shell, in agreement with what has been previously reported. We see evidence for significant photometric variability in the progenitor, as well as strong H$α$ emission consistent with pre-existing circumstellar material. The age of the environment as well as the resolved stellar population surrounding AT2016jbu, support a progenitor age of $>$10 Myr, consistent with a progenitor mass of $\sim$22~$M_{\odot}$. A joint analysis of the velocity evolution of AT2016jbu, and the photospheric radius inferred from the bolometric lightcurve shows the transient is consistent with two successive outbursts/explosions. The first outburst ejected material with velocity $\sim$650$kms^{-1}$, while the second, more energetic event, ejected material at $\sim$4500$kms^{-1}$. Whether the latter is the core-collapse of the progenitor remains uncertain. We place a limit on the ejected $^{56}$Ni mass of $<$0.016$M_{\odot}$. Using the BPASS code, we explore a wide range of possible progenitor systems, and find that the majority of these are in binaries, some of which are undergoing mass transfer or common envelope evolution immediately prior to explosion. Finally, we use the SNEC code to demonstrate that the low-energy explosion within some of these binary systems, together with sufficient CSM, can reproduce the overall morphology of the lightcurve of AT2016jbu. △ Less

Submitted 27 April, 2022; v1 submitted 18 February, 2021; originally announced February 2021.

Comments: 23 pages, 15 figures, accepted for publication in MNRAS

Photometric and spectroscopic evolution of the interacting transient AT 2016jbu (Gaia16cfr)

Authors: S. J. Brennan, M. Fraser, J. Johansson, A. Pastorello, R. Kotak, H. F. Stevance, T. -W. Chen, J. J. Eldridge, S. Bose, P. J. Brown, E. Callis, R. Cartier, M. Dennefeld, Subo Dong, P. Duffy, N. Elias-Rosa, G. Hosseinzadeh, E. Hsiao, H. Kuncarayakti, A. Martin-Carrillo, B. Monard, A. Nyholm, G. Pignata, D. Sand, B. J. Shappee , et al. (46 additional authors not shown)

Abstract: We present the results from a high cadence, multi-wavelength observation campaign of AT 2016jbu (aka Gaia16cfr), an interacting transient. This dataset complements the current literature by adding higher cadence as well as extended coverage of the lightcurve evolution and late-time spectroscopic evolution. Photometric coverage reveals that AT 2016jbu underwent significant photometric variability f… ▽ More We present the results from a high cadence, multi-wavelength observation campaign of AT 2016jbu (aka Gaia16cfr), an interacting transient. This dataset complements the current literature by adding higher cadence as well as extended coverage of the lightcurve evolution and late-time spectroscopic evolution. Photometric coverage reveals that AT 2016jbu underwent significant photometric variability followed by two luminous events, the latter of which reached an absolute magnitude of M$_V\sim$-18.5 mag. This is similar to the transient SN 2009ip whose nature is still debated. Spectra are dominated by narrow emission lines and show a blue continuum during the peak of the second event. AT 2016jbu shows signatures of a complex, non-homogeneous circumstellar material (CSM). We see slowly evolving asymmetric hydrogen line profiles, with velocities of 500km$s^{-1}$ seen in narrow emission features from a slow moving CSM, and up to 10,000km$s^{-1}$ seen in broad absorption from some high velocity material. Late-time spectra ($\sim$+1 year) show a lack of forbidden emission lines expected from a core-collapse supernova and are dominated by strong emission from H, He i and Ca ii. Strong asymmetric emission features, a bumpy lightcurve, and continually evolving spectra suggest an inhibit nebular phase. We compare the evolution of H$α$ among SN 2009ip-like transients and find possible evidence for orientation angle effects. The light-curve evolution of AT 2016jbu suggests similar, but not identical, circumstellar environments to other SN 2009ip-like transients. △ Less

Submitted 27 April, 2022; v1 submitted 18 February, 2021; originally announced February 2021.

Comments: 27 pages, 21 figures, accepted for publication in MNRAS

The orbit and stellar masses of the archetype colliding-wind binary WR 140

Authors: Joshua D. Thomas, Noel D. Richardson, J. J. Eldridge, Gail H. Schaefer, John D. Monnier, Hugues Sana, Anthony F. J. Moffat, Peredur Williams, Michael F. Corcoran, Ian R. Stevens, Gerd Weigelt, Farrah D. Zainol, Narsireddy Anugu, Jean-Baptiste Le Bouquin, Theo ten Brummelaar, Fran Campos, Andrew Couperus, Claire L. Davies, Jacob Ennis, Thomas Eversberg, Oliver Garde, Tyler Gardner, Joan Guarro Fló, Stefan Kraus, Aaron Labdon , et al. (11 additional authors not shown)

Abstract: We present updated orbital elements for the Wolf-Rayet (WR) binary WR\,140 (HD\,193793; WC7pd + O5.5fc). The new orbital elements were derived using previously published measurements along with {\color{black}160} new radial velocity measurements across the 2016 periastron passage of WR 140. Additionally, four new measurements of the orbital astrometry were collected with the CHARA Array. With thes… ▽ More We present updated orbital elements for the Wolf-Rayet (WR) binary WR\,140 (HD\,193793; WC7pd + O5.5fc). The new orbital elements were derived using previously published measurements along with {\color{black}160} new radial velocity measurements across the 2016 periastron passage of WR 140. Additionally, four new measurements of the orbital astrometry were collected with the CHARA Array. With these measurements, we derive stellar masses of $M_{\rm WR} = 10.31\pm0.45 M_\odot$ and $M_{\rm O} = 29.27\pm1.14 M_{\odot}$. We also include a discussion of the evolutionary history of this system from the Binary Population and Spectral Synthesis (BPASS) model grid to show that this WR star likely formed primarily through mass loss in the stellar winds, with only a moderate amount of mass lost or transferred through binary interactions. △ Less

Submitted 23 April, 2021; v1 submitted 26 January, 2021; originally announced January 2021.

Comments: 10 pages, 5 figures

Binary Fraction Indicators in Resolved Stellar Populations and Supernova Type Ratios

Authors: E. R. Stanway, J. J. Eldridge, A. A. Chrimes

Abstract: The binary fraction of a stellar population can have pronounced effects on its properties, and in particular the number counts of different massive star types, and the relative subtype rates of the supernovae which end their lives. Here we use binary population synthesis models with a binary fraction that varies with initial mass to test the effects on resolved stellar populations and supernovae,… ▽ More The binary fraction of a stellar population can have pronounced effects on its properties, and in particular the number counts of different massive star types, and the relative subtype rates of the supernovae which end their lives. Here we use binary population synthesis models with a binary fraction that varies with initial mass to test the effects on resolved stellar populations and supernovae, and ask whether these can constrain the poorly-known binary fraction in different mass and metallicity regimes. We show that Wolf-Rayet star subtype ratios are valuable binary diagnostics, but require large samples to distinguish by models. Uncertainties in which stellar models would be spectroscopically classified as Wolf-Rayet stars are explored. The ratio of thermonuclear, stripped envelope and other core-collapse supernovae may prove a more accessible test and upcoming surveys will be sufficient to constrain both the high mass and low mass binary fraction in the z < 1 galaxy population. △ Less

Submitted 14 July, 2020; originally announced July 2020.

Comments: 10 pages, accepted for publication in MNRAS. Author's accepted version

Revisiting the Impact of Dust Production from Carbon-Rich Wolf-Rayet Binaries

Authors: Ryan M. Lau, J. J. Eldridge, Matthew J. Hankins, Astrid Lamberts, Itsuki Sakon, Peredur M. Williams

Abstract: We present a dust spectral energy distribution (SED) and binary stellar population analysis revisiting the dust production rates (DPRs) in the winds of carbon-rich Wolf-Rayet (WC) binaries and their impact on galactic dust budgets. DustEM SED models of 19 Galactic WC ``dustars" reveal DPRs of $\dot{M}_d\sim10^{-10}-10^{-6}$ M$_\odot$ yr$^{-1}$ and carbon dust condensation fractions, $χ_C$, between… ▽ More We present a dust spectral energy distribution (SED) and binary stellar population analysis revisiting the dust production rates (DPRs) in the winds of carbon-rich Wolf-Rayet (WC) binaries and their impact on galactic dust budgets. DustEM SED models of 19 Galactic WC ``dustars" reveal DPRs of $\dot{M}_d\sim10^{-10}-10^{-6}$ M$_\odot$ yr$^{-1}$ and carbon dust condensation fractions, $χ_C$, between $0.002 - 40\%$. A large ($0.1 - 1.0$ $μ$m) dust grain size composition is favored for efficient dustars where $χ_C\gtrsim1\%$. Results for dustars with known orbital periods verify a power-law relation between $χ_C$, orbital period, WC mass-loss rate, and wind velocity consistent with predictions from theoretical models of dust formation in colliding-wind binaries. We incorporated dust production into Binary Population and Spectral Synthesis (BPASS) models to analyze dust production rates from WC dustars, asymptotic giant branch stars (AGBs), red supergiants (RSGs), and core-collapse supernovae (SNe). BPASS models assuming constant star formation (SF) and a co-eval $10^6$ M$_\odot$ stellar population were performed at low, Large Magellanic Cloud (LMC)-like, and solar metallicities (Z = 0.001, 0.008, and 0.020). Both constant SF and co-eval models show that SNe are net dust destroyers at all metallicities. Constant SF models at LMC-like metallicities show that AGB stars slightly outproduce WC binaries and RSGs by factors of $2-3$, whereas at solar metallicites WC binaries are the dominant source of dust for $\sim60$ Myr until the onset of AGBs, which match the dust input of WC binaries. Co-eval population models show that for "bursty" SF, AGB stars dominate dust production at late times ($t\gtrsim 70$ Myr). △ Less

Submitted 15 June, 2020; originally announced June 2020.

Comments: 34 pages, 12 figures, 7 tables, Accepted for publication in ApJ

Population and spectral synthesis: it doesn't work without binaries

Authors: J. J. Eldridge, Elizabeth R. Stanway

Abstract: In this chapter we discuss the population and spectral synthesis of stellar populations. We describe the method required to achieve such synthesis and discuss examples where inclusion of interacting binaries are vital to reproducing the properties of observed stellar systems. These examples include the Hertzsprung-Russell diagram, massive star number counts, core-collapse supernovae and the ionisi… ▽ More In this chapter we discuss the population and spectral synthesis of stellar populations. We describe the method required to achieve such synthesis and discuss examples where inclusion of interacting binaries are vital to reproducing the properties of observed stellar systems. These examples include the Hertzsprung-Russell diagram, massive star number counts, core-collapse supernovae and the ionising radiation from stellar populations that power both nearby HII regions and the epoch of reionization. We finally offer some speculations on the future paths of research in spectral synthesis. △ Less

Submitted 12 June, 2020; v1 submitted 24 May, 2020; originally announced May 2020.

Comments: Authors' copy of manuscript that appears in, The Impact of Binaries on Stellar Evolution, Beccari G. & Boffin H.M.J. (Eds.). Copyright, 2018 Cambridge University Press. Full text to be found here: https://www.cambridge.org/core/books/impact-of-binary-stars-on-stellar-evolution/5F7496DDDE170B26657D66D5F742906E ). 24 pages, 4 figues

Journal ref: Population and Spectral Synthesis: It Doesn't Work without Binaries. (2019). In G. Beccari & H. Boffin (Eds.), The Impact of Binary Stars on Stellar Evolution (Cambridge Astrophysics, pp. 307-328). Cambridge: Cambridge University Press

Population Synthesis of Massive Close Binary Evolution

Authors: J. J. Eldridge

Abstract: Binary population synthesis is the method by which predictions of varied observables of stellar populations can be made from theoretical models of binary stellar evolution. Binary stars have many more possible evolutionary outcomes compared to single stars and the relative rates of the different pathways, such as the rates of different supernova types, depend on uncertain or poorly constrained phy… ▽ More Binary population synthesis is the method by which predictions of varied observables of stellar populations can be made from theoretical models of binary stellar evolution. Binary stars have many more possible evolutionary outcomes compared to single stars and the relative rates of the different pathways, such as the rates of different supernova types, depend on uncertain or poorly constrained physics. In this Chapter we describe population synthesis, outline the major uncertainties and discuss the relevant predictions for core-collapse supernovae. After we overview single star evolution we outline the important physical processes that occur in binaries including Roche-lobe overflow, common-envelope evolution and supernova kicks. We also discuss how a synthetic stellar population incorporating interacting binaries can be constructed and how uncertainties, such as the strength of supernova kicks, affect any predictions. We illustrate the process by comparing predictions for the stellar populations in two young star clusters. We then discuss the important predictions from population synthesis for understanding core-collapse supernovae, their delay-time distribution and their progenitor stars. Finally we discuss how we can predict the rate of mergers of compact remnants and thus predict the initial parameters of gravitational wave sources. △ Less

Submitted 12 June, 2020; v1 submitted 24 May, 2020; originally announced May 2020.

Comments: Author version of a chapter for the 'Handbook of Supernovae', edited by A. Alsabti and P. Murdin, Springer. 22 pages, 4 figures, 1 table

SNAPSHOT: Connections between Internal and Surface Properties of Massive Stars

Authors: Eoin Farrell, Jose Groh, Georges Meynet, JJ Eldridge, Sylvia Ekström, Cyril Georgy

Abstract: We introduce SNAPSHOT, a technique to systematically compute stellar structure models in hydrostatic and thermal equilibrium based on 3 structural properties - core mass $M_{\rm core}$, envelope mass $M_{\rm env}$ and core composition. This approach allows us to connect these properties of stellar interiors to the luminosity and effective temperature $T_{\rm eff}$ in a more systematic way than wit… ▽ More We introduce SNAPSHOT, a technique to systematically compute stellar structure models in hydrostatic and thermal equilibrium based on 3 structural properties - core mass $M_{\rm core}$, envelope mass $M_{\rm env}$ and core composition. This approach allows us to connect these properties of stellar interiors to the luminosity and effective temperature $T_{\rm eff}$ in a more systematic way than with stellar evolution models. For MS models, we derive an analytical relationship between $M_{\rm core}$, $M_{\rm total}$ and central H abundance that can be used in rapid stellar evolution algorithms. Core-He burning models with $M_{\rm core}/M_{\rm total}$ from 0.2 to 0.8 have convective envelopes, low $T_{\rm eff}$ and will appear as red supergiants. For a given $M_{\rm core}$, they exhibit a small variation in luminosity (0.02 dex) and $T_{\rm eff}$ ($\sim 400\,\mathrm{K}$) over a wide range of $M_{\rm env}$ ($\sim 2 - 20\,\mathrm{M}_{\odot}$). This means that it is not possible to derive red supergiant masses from luminosities and $T_{\rm eff}$ alone. We derive the following relationship between $M_{\rm core}$ and the total luminosity of a red supergiant during core He burning: $\log M_{\rm core} \simeq 0.44 \log L/L_{\odot} - 1.38$. At $M_{\rm core}$/$M_{\rm total} \approx 0.2$, our models exhibit a bi-stability and jump from a RSG to a BSG structure. Our models with $M_{\rm core}/M_{\rm total} > 0.8$, which correspond to stripped stars produced by mass loss or binary interaction, show that $T_{\rm eff}$ has a strong dependence on $M_{\rm env}$, $M_{\rm core}$ and the core composition. We find the mass of one of these stripped stars in a binary system, HD 45166, to be less than its dynamical mass. When a large observational sample of stripped stars becomes available, our results can be used to constrain their masses and the physics of binary interaction. △ Less

Submitted 13 May, 2020; originally announced May 2020.

Comments: Accepted for publication in MNRAS

Evaluating the impact of binary parameter uncertainty on stellar population properties

Authors: E. R. Stanway, A. A. Chrimes, J. J. Eldridge, H. F. Stevance, .

Abstract: Binary stars have been shown to have a substantial impact on the integrated light of stellar populations, particularly at low metallicity and early ages - conditions prevalent in the distant Universe. But the fraction of stars in stellar multiples as a function of mass, their likely initial periods and distribution of mass ratios are all known empirically from observations only in the local Univer… ▽ More Binary stars have been shown to have a substantial impact on the integrated light of stellar populations, particularly at low metallicity and early ages - conditions prevalent in the distant Universe. But the fraction of stars in stellar multiples as a function of mass, their likely initial periods and distribution of mass ratios are all known empirically from observations only in the local Universe. Each has associated uncertainties. We explore the impact of these uncertainties in binary parameters on the properties of integrated stellar populations, considering which properties and timescales are most susceptible to uncertainty introduced by binary fractions and whether observations of the integrated light might be sufficient to determine binary parameters. We conclude that the effects of uncertainty in the empirical binary parameter distributions are likely smaller than those introduced by metallicity and stellar population age uncertainties for observational data. We identify emission in the He II 1640 Angstrom emission line and continuum colour in the ultraviolet-optical as potential indicators of a high mass binary presence, although poorly constrained metallicity, dust extinction and degeneracies in plausible star formation history are likely to swamp any measurable signal. △ Less

Submitted 24 April, 2020; originally announced April 2020.

Comments: 18 pages + appendix. Accepted for publication in MNRAS. Author's accepted version

On aging star clusters using red supergiants independent of the fraction of interacting binary stars

Authors: J. J. Eldridge, Emma R. Beasor, N. Britavskiy

Abstract: We use the Binary Population and Spectral Synthesis (BPASS) models to test the recent suggestion that red supergiants can provide an accurate age estimate of a co-eval stellar population that is unaffected by interacting binary stars. Ages are estimated by using both the minimum luminosity red supergiant and the mean luminosity of red supergiants in a cluster. We test these methods on a number of… ▽ More We use the Binary Population and Spectral Synthesis (BPASS) models to test the recent suggestion that red supergiants can provide an accurate age estimate of a co-eval stellar population that is unaffected by interacting binary stars. Ages are estimated by using both the minimum luminosity red supergiant and the mean luminosity of red supergiants in a cluster. We test these methods on a number of observed star clusters and find our results in agreement with previous estimates. Importantly we find the difference between the ages derived from stellar population models with and without a realistic population of interacting binary stars is only a few 100,000 years at most. We find that the mean luminosity of red supergiants in a cluster is the best method to determine the age of a cluster because it is based o the entire red supergiant population rather than using only the least luminous red supergiant. △ Less

Submitted 17 April, 2020; originally announced April 2020.

Comments: Accepted for publication in MNRAS, 7 pages, 3 pages, 2 tables

A systematic aging method I: HII regions D118 and D119 in NGC 300

Authors: H. F. Stevance, J. J. Eldridge, A. McLeod, E. R. Stanway, A. A. Chrimes

Abstract: Accurately determining the age of H\,{ii} regions and the stars they host is as important as it is challenging. Historically the most popular method has been isochrone fitting to Hertzsprung-Russell Diagrams or Colour-Magnitude Diagrams. Here we introduce a different method for age determination using BPASS and hoki. We infer the most likely ages of the regions D118 and D119 NGC 300 to be log(age/… ▽ More Accurately determining the age of H\,{ii} regions and the stars they host is as important as it is challenging. Historically the most popular method has been isochrone fitting to Hertzsprung-Russell Diagrams or Colour-Magnitude Diagrams. Here we introduce a different method for age determination using BPASS and hoki. We infer the most likely ages of the regions D118 and D119 NGC 300 to be log(age/years)=6.86$^{0.05}_{0.06}$ and we also deduce stellar mass and number counts by comparison with the BPASS models. We compare how our binary and single star models perform and find that the latter are unable to predict 20 per cent ($\pm$ 10 per cent) of our sample. We also discuss how results obtained from isochrone fitting would differ. We conclude that ages could be underestimated by ~0.2 dex and that the limitations of the isochrone method is not solely due to the lack of binary stars. We propose that the method presented here is more reliable and more widely applicable since it can be used on smaller samples. Alongside this study, we release new hoki features to allow easy implementation of this method. △ Less

Submitted 10 August, 2020; v1 submitted 6 April, 2020; originally announced April 2020.

Comments: Submitted to MNRAS, comments welcome, 17 pages, 14 figures, 3 Tables

Hoki: Making BPASS accessible through Python

Authors: Heloise F. Stevance, J. J. Eldridge, Elizabeth Stanway

Abstract: We now know that a large number of stars are born in multiple systems. Additionally, more than 70% of massive stars are found in close binary systems, meaning that they will interact over the course of their lifetime. This has strong implications for their evolution as well as the transients (e.g supernovae) and the potential gravitational wave progenitors they produce. Therefore, in order to unde… ▽ More We now know that a large number of stars are born in multiple systems. Additionally, more than 70% of massive stars are found in close binary systems, meaning that they will interact over the course of their lifetime. This has strong implications for their evolution as well as the transients (e.g supernovae) and the potential gravitational wave progenitors they produce. Therefore, in order to understand and correctly interpret astronomical observations of stellar populations, we must use theoretical models able to account for the effects of binary stars. This is the case of the Binary Population and Spectral Synthesis code (BPASS), which has been a staple of the field for over 10 years. As is the case for most other theoretical models, the data products of BPASS are large, varied and complex. As a result, their use requires a level of expertise that is not immediately accessible to a wider community that may hold key observational data. The goal of hoki is to bridge the gap between observation and theory, by providing a set of tools to make BPASS data easily accessible and facilitate analysis. The use of Python is deliberate as it is a ubiquitous language within Astronomy. This allows BPASS results to be used naturally within the pre-existing workflow of most astronomers. △ Less

Submitted 29 January, 2020; originally announced January 2020.

Comments: 3 pages, Published in JOSS, GitHub: https://github.com/HeloiseS/hoki