Researcher profile

J. J. Eldridge

J. J. Eldridge contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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

preprint2022arXiv

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

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 &#39;CSM population synthesis&#39; 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.

preprint2022arXiv

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

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.

preprint2021arXiv

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

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&#39;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&#39;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.

preprint2021arXiv

Estimating Transient Rates from Cosmological Simulations and BPASS

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.

preprint2020arXiv

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

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.

preprint2020arXiv

Binary Fraction Indicators in Resolved Stellar Populations and Supernova Type Ratios

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.

preprint2020arXiv

Evaluating the impact of binary parameter uncertainty on stellar population properties

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.

preprint2020arXiv

Hoki: Making BPASS accessible through Python

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.

preprint2020arXiv

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

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.

preprint2020arXiv

Population and spectral synthesis: it doesn&#39;t work without binaries

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.

preprint2020arXiv

Population Synthesis of Massive Close Binary Evolution

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.

preprint2020arXiv

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

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&#34; 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 &#34;bursty&#34; SF, AGB stars dominate dust production at late times ($t\gtrsim 70$ Myr).

preprint2020arXiv

Weighing in on black hole binaries with BPASS: LB-1 does not contain a 70M$_{\odot}$ black hole

The recent identification of a candidate very massive 70 M(Sun) black hole is at odds with our current understanding of stellar winds and pair-instability supernovae. We investigate alternate explanations for this system by searching the BPASS v2.2 stellar and population synthesis models for those that match the observed properties of the system. We find binary evolution models that match the LB-1 system, at the reported Gaia distance, with more moderate black hole masses of 4 to 7 M(Sun). We also examine the suggestion that the binary motion may have led to an incorrect distance determination by Gaia. We find that the Gaia distance is accurate and that the binary system is consistent with the observation at this distance. Consequently it is highly improbable that the black hole in this system has the extreme mass originally suggested. Instead, it is more likely to be representative of the typical black hole binary population expected in our Galaxy.

preprint2019arXiv

Dependence of Gravitational Wave Transient Rates on Cosmic Star Formation and Metallicity Evolution History

We compare the impacts of uncertainties in both binary population synthesis models and the cosmic star formation history on the predicted rates of Gravitational Wave compact binary merger (GW) events. These uncertainties cause the predicted rates of GW events to vary by up to an order of magnitude. Varying the volume-averaged star formation rate density history of the Universe causes the weakest change to our predictions, while varying the metallicity evolution has the strongest effect. Double neutron-star merger rates are more sensitive to assumed neutron-star kick velocity than the cosmic star formation history. Varying certain parameters affects merger rates in different ways depending on the mass of the merging compact objects; thus some of the degeneracy may be broken by looking at all the event rates rather than restricting ourselves to one class of mergers.

preprint2019arXiv

Supernova lightCURVE POPulation Synthesis II: Validation against supernovae with an observed progenitor

We use the results of a supernova light-curve population synthesis to predict the range of possible supernova light curves arising from a population of single-star progenitors that lead to type IIP supernovae. We calculate multiple models varying the initial mass, explosion energy, nickel mass and nickel mixing and then compare these to type IIP supernovae with detailed light curve data and pre-explosion imaging progenitor constraints. Where a good fit is obtained to observations, we are able to achieve initial progenitor and nickel mass estimates from the supernova lightcurve that are comparable in precision to those obtained from progenitor imaging. For two of the eleven IIP supernovae considered our fits are poor, indicating that more progenitor models should be included in our synthesis or that our assumptions, regarding factors such as stellar mass loss rates or the rapid final stages of stellar evolution, may need to be revisited in certain cases. Using the results of our analysis we are able to show that most of the type IIP supernovae have an explosion energy of the order of log(E_exp/ergs)=50.52+/-0.10 and that both the amount of nickel in the supernovae and the amount of mixing may have a dependence on initial progenitor mass.