BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

J. J. Eldridge

J. J. Eldridge appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-ph.GAastro-ph.HEastro-ph.COastro-ph.IMgr-qc

Catalog footprint

What is connected

42works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

42 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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

0 citations0 reviewsNo code linkedOpen access
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 '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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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'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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Population and spectral synthesis: it doesn'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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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" 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).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Comprehensive Comparative Test of Seven Widely-Used Spectral Synthesis Models Against Multi-Band Photometry of Young Massive Star Clusters

We test the predictions of spectral synthesis models based on seven different massive-star prescriptions against Legacy ExtraGalactic UV Survey (LEGUS) observations of eight young massive clusters in two local galaxies, NGC 1566 and NGC 5253, chosen because predictions of all seven models are available at the published galactic metallicities. The high angular resolution, extensive cluster inventory and full near-ultraviolet to near-infrared photometric coverage make the LEGUS dataset excellent for this study. We account for both stellar and nebular emission in the models and try two different prescriptions for attenuation by dust. From Bayesian fits of model libraries to the observations, we find remarkably low dispersion in the median E(B-V) (~0.03 mag), stellar masses (~10^4 M_\odot) and ages (~1 Myr) derived for individual clusters using different models, although maximum discrepancies in these quantities can reach 0.09 mag and factors of 2.8 and 2.5, respectively. This is for ranges in median properties of 0.05-0.54 mag, 1.8-10x10^4 M_\odot and 1.6-40 Myr spanned by the clusters in our sample. In terms of best fit, the observations are slightly better reproduced by models with interacting binaries and least well reproduced by models with single rotating stars. Our study provides a first quantitative estimate of the accuracies and uncertainties of the most recent spectral synthesis models of young stellar populations, demonstrates the good progress of models in fitting high-quality observations, and highlights the needs for a larger cluster sample and more extensive tests of the model parameter space.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

BPASS predictions for Binary Black-Hole Mergers

Using the Binary Population and Spectral Synthesis code BPASS, we have calculated the rates, timescales and mass distributions for binary black hole mergers as a function of metallicity. We consider these in the context of the recently reported 1st LIGO event detection. We find that the event has a very low probability of arising from a stellar population with initial metallicity mass fraction above Z=0.010 (Z>0.5Zsun). Binary black hole merger events with the reported masses are most likely in populations below 0.008 (Z<0.4Zsun). Events of this kind can occur at all stellar population ages from ~3 Myr up to the age of the universe, but constitute only 0.1 to 0.4 per cent of binary BH mergers between metallicities of Z=0.001 to 0.008. However at metallicity Z=0.0001, 26 per cent of binary BH mergers would be expected to have the reported masses. At this metallicity the progenitor merger times can be close to ~10Gyr and rotationally-mixed stars evolving through quasi-homogeneous evolution, due to mass transfer in a binary, dominate the rate. The masses inferred for the black holes in the binary progenitor of GW\,150914 are amongst the most massive expected at anything but the lowest metallicities in our models. We discuss the implications of our analysis for the electromagnetic follow-up of future LIGO event detections.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Helium Stars: Towards an Understanding of Wolf-Rayet Evolution

Wolf-Rayet (WR) stars are massive stars that have lost most or all of their hydrogen via powerful stellar winds. Recent observations have indicated that hydrogen-free WR stars have cooler temperatures than those predicted by current evolutionary models. To investigate how varying mass-loss rate affects WR evolution, we have created a grid of pure helium star models. We compare our results with Galactic and LMC WR observations and show that the temperature ranges of observed WR stars can be reproduced by varying the mass-loss rate, which effectively determines the size of the helium envelope around the core. We also find that WN and WO stars arise from more massive stars, whereas WC stars come from lower masses. This contradicts the standard Conti scenario by which WN and WC stars evolve in an age sequence. We also predict the magnitudes of our models at core-collapse and compare with observations of nearby progenitors of Type Ib/c supernovae. We confirm the findings of previous studies that suggest WR stars are the progenitors of core-collapse supernovae: the progenitors would remain unobserved except in the cases where the progenitor is a low-mass helium giant, as is the case of iPTF13bvn.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Neutron Star Kicks and their Relationship to Supernovae Ejecta Mass

We propose a simple model to explain the velocity of young neutron stars. We attempt to confirm a relationship between the amount of mass ejected in the formation of the neutron star and the `kick' velocity imparted to the compact remnant resulting from the process. We assume the velocity is given by $v_{\rm kick}=α\,(M_{\rm ejecta} / M_{\rm remnant}) + β\,$. To test this simple relationship we use the BPASS (Binary Population and Spectral Synthesis) code to create stellar population models from both single and binary star evolutionary pathways. We then use our Remnant Ejecta and Progenitor Explosion Relationship (REAPER) code to apply different $α$ and $β$ values and three different `kick' orientations then record the resulting velocity probability distributions. We find that while a single star population provides a poor fit to the observational data, the binary population provides an excellent fit. Values of $α=70\, {\rm km\,s^{-1}}$ and $β=110\,{\rm km\,s^{-1}}$ reproduce the \cite{RN165} observed 2-dimensional velocities and $α=70\, {\rm km\,s^{-1}}$ and $β=120\,{\rm km\,s^{-1}}$ reproduce their inferred 3-dimensional velocity distribution for nearby single neutron stars with ages less than 3 Myrs. After testing isotropic, spin-axis aligned and orthogonal to spin-axis `kick' orientations, we find no statistical preference for a `kick' orientation. While ejecta mass cannot be the only factor that determines the velocity of supernovae compact remnants, we suggest it is a significant contributor and that the ejecta based `kick' should replace the Maxwell-Boltzmann velocity distribution currently used in many population synthesis codes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Progenitor constraints for core-collapse supernovae from Chandra X-ray observations

The progenitors of hydrogen-poor core-collapse supernovae (SNe) of types Ib, Ic and IIb are believed to have shed their outer hydrogen envelopes either by extremely strong stellar winds, characteristic of classical Wolf-Rayet stars, or by binary interaction with a close companion star. The exact nature of the progenitors and the relative importance of these processes are still open questions. One relatively unexplored method to constrain the progenitors is to search for high-mass X-ray binaries (HMXB) at SN locations in pre-explosion X-ray observations. In a HMXB, one star has already exploded as a core-collapse SN, producing a neutron star or a stellar-mass black hole. It is likely that the second star in the system will also explode as a supernova, which should cause a detectable long-term change in the system's X-ray luminosity. In particular, a pre-explosion detection of a HMXB coincident with a SN could be informative about the progenitor's nature. In this paper we analyze pre-explosion ACIS observations of 18 nearby type Ib, Ic and IIb supernovae from the Chandra X-ray observatory public archive. Two sources that could potentially be associated with the supernova are identified in the sample. Additionally we make similar post-explosion measurements for 46 SNe. Although our modelling indicates that progenitor systems with compact binary companions are probably quite rare, studies of this type can in the future provide more stringent constraints as the number of discovered nearby SNe and suitable pre-explosion X-ray data are both increasing.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Rapidly evolving faint transients from stripped-envelope electron-capture supernovae

We investigate the expected rates and bolometric light-curve properties of stripped-envelope electron-capture supernovae (ECSNe) using stellar models from the Binary Population and Spectral Synthesis (BPASS) code. We find that 0.8 per cent (Z = 0.020) and 1.2 per cent (Z = 0.004) of core-collapse supernovae are stripped-envelope ECSNe. Their typical ejecta masses are estimated to be about 0.3 Msun (Z = 0.020) and 0.6 Msun (Z = 0.004). Assuming ECSN explosion properties from numerical explosion simulations, an explosion energy of 1.5e50 erg and a 56Ni mass of 2.5e-3 Msun, we find that stripped-envelope ECSNe have a typical rise time of around 7 days (Z = 0.020) or 13 days (Z = 0.004) and peak luminosity of around 1e41 erg/s (-13.8 mag, Z = 0.020) or 7e40 erg/s (-13.4 mag, Z = 0.004). Their typical ejecta velocities are around 7000 km/s (Z = 0.020) or 5000 km/s (Z = 0.004). Thus, stripped-envelope ECSNe are observed as rapidly-evolving faint transients with relatively small velocities. SN 2008ha-like supernovae, which are the faintest kind of SN 2002cx-like (a.k.a. Type Iax) supernovae, may be related to stripped-envelope ECSNe.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The disappearance of the helium-giant progenitor of the type Ib supernova iPTF13bvn and constraints on its companion

We report and discuss post-explosion observations of supernova iPTF13bvn. We find that the brightness of the SN at +740 days is below the level of the pre-explosion source and thus confirm that the progenitor star has gone. We estimate that the late-time brightness is still dominated by the supernova, which constrains the magnitude and thus mass of a possible companion star to below approximately 10Msun. In turn this implies that the progenitor's initial mass is constrained to a narrow range of between 10 to 12Msun. The progenitor of iPTF13bvn would have been a helium giant rather than a Wolf-Rayet star. In addition, we suggest that sufficiently deep observations acquired in 2016 would now stand a chance to directly observe the companion star.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Wolf-Rayet stars in the Small Magellanic Cloud: II. Analysis of the binaries

Massive WR stars are evolved massive stars characterized by strong mass-loss. Hypothetically, they can form either as single stars or as mass donors in close binaries. About 40% of the known WR stars are confirmed binaries, raising the question as to the impact of binarity on the WR population. By performing a spectral analysis of all multiple WR systems in the SMC, we obtain the full set of stellar parameters for each individual component. Mass-luminosity relations are tested, and the importance of the binary evolution channel is assessed. The spectral analysis is performed with the PoWR model atmosphere code by superimposing model spectra that correspond to each component. Evolutionary channels are constrained using the BPASS evolution tool. Significant Hydrogen mass fractions (0.1 - 0.4) are detected in all WN components. A comparison with mass-luminosity relations and evolutionary tracks implies that the majority of the WR stars in our sample are not chemically homogeneous. The WR component in the binary AB 6 is found to be very luminous (Log L ~ 6.3 [Lsun]) given its orbital mass (~10 Msun), presumably because of observational contamination by a third component. Evolutionary paths derived for our objects suggest that Roche lobe overflow had occurred in most systems, affecting their evolution. However, the implied initial masses are large enough for the primaries to have entered the WR phase, regardless of binary interaction. Together with the results for the putatively single SMC WR stars, our study suggests that the binary evolution channel does not dominate the formation of WR stars at SMC metallicity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Turnoff Detached Binary Star V568 Lyr in the Kepler Field of the Oldest Open Cluster (NGC 6791) in the Galaxy

We present the Kepler photometric light-variation analysis of the late-type double-lined binary system V568 Lyr that is in the field of the high metallicity old open cluster NGC 6791. The radial velocity and the high-quality short-cadence light curve of the system are analysed simultaneously. The masses, radii and luminosities of the component stars are $M_1 = 1.0886\pm0.0031\, M{\odot}$, $M_2 = 0.8292 \pm 0.0026\, M{\odot}$, $R_1 = 1.4203\pm 0.0058\, R{\odot}$, $R_2 = 0.7997 \pm 0.0015\, R{\odot}$, $L_1 = 1.85\pm 0.15\, L{\odot}$, $L_2 = 0.292 \pm 0.018\, L{\odot}$ and their separation is $a = 31.060 \pm 0.002\, R{\odot}$. The distance to NGC 6791 is determined to be $4.260\pm 0.290\,$kpc by analysis of this binary system. We fit the components of this well-detached binary system with evolution models made with the Cambridge STARS and TWIN codes to test low-mass binary star evolution. We find a good fit with a metallicity of $Z = 0.04$ and an age of $7.704\,$Gyr. The standard tidal dissipation, included in TWIN is insufficient to arrive at the observed circular orbit unless it formed rather circular to begin with.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Core-Collapse Supernova Rate Synthesis Within 11 Mpc

The 11 Mpc H-alpha and Ultraviolet Galaxy (11HUGS) Survey traces the star formation activity of nearby galaxies. In addition within this volume the detection completeness of core-collapse supernovae (CCSNe) is high therefore by comparing these observed stellar births and deaths we can make a sensitive test of our understanding of how stars live and die. In this paper, we use the results of the Binary Population and Spectral Synthesis (BPASS) code to simulate the 11HUGS galaxies H-alpha and far-ultraviolet (FUV) star formation rate indicators (SFRIs) and simultaneously match the core-collapse supernova (CCSN) rate. We find that stellar population including interacting binary stars makes little difference to the total CCSN rate but increases the H-alpha and FUV fluxes for a constant number of stars being formed. In addition they significantly increase the predicted rate of type Ibc supernovae (SNe) relative to type II SNe to the level observed in the 11HUGS galaxies. We also find that instead of assuming a constant star formation history (SFH) for the galaxies our best fitting models have a star formation rate (SFR) that peaked more than 3 Myrs ago.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Did the progenitor of SN2011dh have a binary companion?

We present late-time Hubble Space Telescope (HST) ultraviolet (UV) and optical observations of the site of SN 2011dh in the galaxy M51, ~1164 days post-explosion. At the SN location, we observe a point source that is visible at all wavelengths, that is significantly fainter than the spectral energy distribution (SED) of the Yellow Supergiant progenitor observed prior to explosion. The previously reported photometry of the progenitor is, therefore, completely unaffected by any sources that may persist at the SN location after explosion. In comparison with the previously reported late-time photometric evolution of SN 2011dh, we find that the light curve has plateaued at all wavelengths. The SED of the late-time source is clearly inconsistent with a SED of stellar origin. Although the SED is bright at UV wavelengths, there is no strong evidence that the late-time luminosity originates solely from a stellar source corresponding to the binary companion, although a partial contribution to the observed UV flux from a companion star can not be ruled out.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the Blue Loops of Intermediate-Mass Stars

We consider the blue loops in the Hertzsprung-Russell diagram that occur when intermediate-mass stars begin core helium burning. It has long been known that the excess of helium above the burning shell, the result of the contraction of the convective core during core hydrogen burning, has the effect of making such stars redder and larger than they would be otherwise. The outward motion of the burning shell in mass removes this excess and triggers the loop. Hitherto nobody has attempted to demonstrate why the excess helium has this effect. We consider the effect of the local opacity, which is reduced by excess helium, the shell fuel supply, which is also reduced, and the local mean molecular weight, which is increased. We demonstrate that the mean molecular weight is the decisive reddening factor. The opacity has a much smaller effect and a reduced fuel supply actually favours blueward motion.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Interpreting high [O III]/Hbeta ratios with maturing starbursts

Star forming galaxies at high redshift show ubiquitously high ionization parameters, as measured by the ratio of optical emission lines. We demonstrate that local (z < 0.2) sources selected as Lyman break analogues also manifest high line ratios with a typical [O III]/H beta = 3.36(+0.14,-0.04) - comparable to all but the highest ratios seen in star forming galaxies at z ~ 2-4. We argue that the stellar population synthesis code BPASS can explain the high ionization parameters required through the ageing of rapidly formed star populations, without invoking any AGN contribution. Binary stellar evolution pathways prolong the age interval over which a starburst is likely to show elevated line ratios, relative to those predicted by single stellar evolution codes. As a result, model galaxies at near-Solar metallicities and with ages of up to ~100 Myr after a starburst typically have a line ratio [O III]/H beta~3, consistent with those seen in Lyman break galaxies and local sources with similar star formation densities. This emphasises the importance of including binary evolution pathways when simulating the nebular line emission of young or bursty stellar populations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Possible binary progenitors for the type Ib supernova iPTF13bvn

Cao et al. (2013) reported a possible progenitor detection for the type Ib supernovae iPTF13bvn for the first time. We find that the progenitor is in fact brighter than the magnitudes previously reported by approximately 0.7 to 0.2 mag with a larger error in the bluer filters. We compare our new magnitudes to our large set of binary evolution models and find that many binary models with initial masses in the range of 10 to 20M(Sun) match this new photometry and other constraints suggested from analysing the supernova. In addition these lower mass stars retain more helium at the end of the model evolution indicating that they are likely to be observed as type Ib supernovae rather than their more massive, Wolf-Rayet counter parts. We are able to rule out typical Wolf-Rayet models as the progenitor because their ejecta masses are too high and they do not fit the observed SED unless they have a massive companion which is the observed source at the supernova location. Therefore only late-time observations of the location will truly confirm if the progenitor was a helium giant and not a Wolf-Rayet star.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A new precise mass for the progenitor of the Type IIP SN 2008bk

The progenitor of the Type IIP SN 2008bk was discovered in pre-explosion g'r'i'IYJHKs images, acquired with European Southern Observatory Very Large Telescope FORS, HAWK-I and ISAAC instruments and the Gemini GMOS-S instrument. The wealth of pre-explosion observations makes the progenitor of this SN one of the best studied, since the detection of the progenitor of SN1987A. Previous analyses of the properties of the progenitor were hampered by the limited quality of the photometric calibration of the pre-explosion images and the crowded nature of the field containing the SN. We present new late-time observations of the site of SN2008bk acquired with identical instrument and filter configurations as the pre-explosion observations, and confirm that the previously identified red supergiant star was the progenitor of this SN and has now disappeared. Image subtraction techniques were used to conduct precise photometry of the now missing progenitor, independently of blending from any nearby stars. The nature of the surrounding stellar population and their contribution to the flux attributed to the progenitor in the pre-explosion images are probed using HST WFC3 UVIS/IR observations. In comparison with MARCS synthetic spectra, we find the progenitor was a highly reddened RSG with luminosity log (L/Lsun)=4.84+/-0.11, corresponding to an initial mass of Minit=12.9+/-1.7Msun. The temperature of the progenitor was hotter than previously expected for RSGs (T ~ 4330K), but consistent with new temperatures derived for RSGs using SED fitting techniques. We show that there is evidence for significant extinction of the progenitor, possibly arising in the CSM; but that this dust yields a similar reddening law to dust found in the ISM (E(B-V)=0.77 with Rv=3.1). [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A transdimensional Bayesian method to infer the star formation history of resolved stellar populations

We propose a new method to infer the star formation histories of resolved stellar populations. With photometry one may plot observed stars on a colour-magnitude diagram (CMD) and then compare with synthetic CMDs representing different star formation histories. This has been accomplished hitherto by parametrising the model star formation history as a histogram, usually with the bin widths set by fixed increases in the logarithm of time. A best fit is then found with maximum likelihood methods and we consider the different means by which a likelihood can be calculated. We then apply Bayesian methods by parametrising the star formation history as an unknown number of Gaussian bursts with unknown parameters. This parametrisation automatically provides a smooth function of time. A Reversal Jump Markov Chain Monte Carlo method is then used to find both the most appropriate number of Gaussians, thus avoiding avoid overfitting, and the posterior probability distribution of the star formation rate. We apply our method to artificial populations and to observed data. We discuss the other advantages of the method: direct comparison of different parametrisations and the ability to calculate the probability that a given star is from a given Gaussian. This allows the investigation of possible sub-populations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On how leakage can affect the Star Formation Rate estimation using Halpha luminosity

We present observational evidence that leakage of ionising photons from star-forming regions can affect the quantification of the star formation rate (SFR) in galaxies. This effect could partially explain the differences between the SFR estimates using the far ultraviolet (FUV) and the Halpha emission. We find that leakage could decrease the SFR(Ha)/SFR(FUV) ratio by up to a 25 per cent. The evidence is based on the observation that the SFR(Ha)/SFR(FUV) ratio is lower for objects showing a shell Halpha structure than for regions exhibiting a much more compact morphology. The study has been performed on three object samples: low luminosity dwarf galaxies from the Local Volume Legacy survey and star-forming regions in the Large Magellanic Cloud and the nearby Local Group galaxy M33. For the three samples we find differences (1.1-1.4sigma) between the SFR(Ha)/SFR(FUV) for compact and shell objects. Although leakage cannot entirely explain the observed trend of SFR(Ha)/SFR(FUV) ratios for systems with low SFR, we show the mechanism can lead to different SFR estimates when using Halpha and FUV luminosities. Therefore, further study is needed to constrain the contribution of leakage to the low SFR(Ha)/SFR(FUV) ratios observed in dwarf galaxies and its impact on the Halpha flux as a SFR indicator in such objects.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Red and dead: The progenitor of SN 2012aw in M95

Core-collapse supernovae (SNe) are the spectacular finale to massive stellar evolution. In this Letter, we identify a progenitor for the nearby core-collapse SN 2012aw in both ground based near-infrared, and space based optical pre-explosion imaging. The SN itself appears to be a normal Type II Plateau event, reaching a bolometric luminosity of 10$^{42}$ erg s$^{-1}$ and photospheric velocities of $\sim$11,000 \kms\ from the position of the H$β$ P-Cygni minimum in the early SN spectra. We use an adaptive optics image to show that the SN is coincident to within 27 mas with a faint, red source in pre-explosion HST+WFPC2, VLT+ISAAC and NTT+SOFI images. The source has magnitudes $F555W$=26.70$\pm$0.06, $F814W$=23.39$\pm$0.02, $J$=21.1$\pm$0.2, $K$=19.1$\pm$0.4, which when compared to a grid of stellar models best matches a red supergiant. Interestingly, the spectral energy distribution of the progenitor also implies an extinction of $A_V>$1.2 mag, whereas the SN itself does not appear to be significantly extinguished. We interpret this as evidence for the destruction of dust in the SN explosion. The progenitor candidate has a luminosity between 5.0 and 5.6 log L/\lsun, corresponding to a ZAMS mass between 14 and 26 \msun\ (depending on $A_V$), which would make this one of the most massive progenitors found for a core-collapse SN to date.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Binary progenitor models of type IIb supernovae

Massive stars that lose their hydrogen-rich envelope down to a few tenths of a solar mass explode as extended type IIb supernovae, an intriguing subtype that links the hydrogen-rich type II supernovae with the hydrogen-poor type Ib and Ic. The progenitors may be very massive single stars that lose their envelope due to their stellar wind, but mass stripping due to interaction with a companion star in a binary system is currently considered to be the dominant formation channel. We computed an extensive grid of binary models with the Eggleton binary evolution code. The predicted rate from our standard models, which assume conservative mass transfer, is about 6 times smaller than the current rate indicated by observations. It is larger but still comparable to the rate expected from single stars. To recover the observed rate we must generously allow for uncertainties and low accretion efficiencies in combination with limited angular momentum loss from the system. Motivated by the claims of detection and non-detection of companions for a few IIb supernovae, we investigate the properties of the secondary star at the moment of explosion. We identify three cases: (1) the companion appears as a hot O star in about 90% of the cases, (2) the companion becomes an over-luminous B star in about 3% of the cases, or (3) the companion will appear as a K supergiant. The second case, which applies to the well-studied case of SN 1993J and possibly to SN 2001ig, is the least common case and requires that the companion very efficiently accretes the transferred material -in contrast to what is required to recover the overall IIb rate. The fast increasing number of type IIb detections may lead to more solid constraints on model uncertainties, like stellar wind-mass loss and internal mixing.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Resolved photometry of extragalactic young massive star clusters

We present colour-magnitude diagrams (CMDs) for a sample of seven young massive clusters in the galaxies NGC 1313, NGC 1569, NGC 1705, NGC 5236 and NGC 7793. The clusters have ages in the range 5-50 million years and masses of 10^5 -10^6 Msun. Although crowding prevents us from obtaining photometry in the central regions of the clusters, we are still able to measure up to 30-100 supergiant stars in each of the richest clusters, along with the brighter main sequence stars. The resulting CMDs and luminosity functions are compared with photometry of artificially generated clusters, designed to reproduce the photometric errors and completeness as realistically as possible. In agreement with previous studies, our CMDs show no clear gap between the H-burning main sequence and the He-burning supergiant stars, contrary to predictions by common stellar isochrones. In general, the isochrones also fail to match the observed number ratios of red-to-blue supergiant stars, although the difficulty of separating blue supergiants from the main sequence complicates this comparison. In several cases we observe a large spread (1-2 mag) in the luminosities of the supergiant stars that cannot be accounted for by observational errors. This spread can be reproduced by including an age spread of 10-30 million years in the models. However, age spreads cannot fully account for the observed morphology of the CMDs and other processes, such as the evolution of interacting binary stars, may also play a role.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

SN 2009md: Another faint supernova from a low mass progenitor

We present adaptive optics imaging of the core collapse supernova (SN) 2009md, which we use together with archival \emph{Hubble Space Telescope} data to identify a coincident progenitor candidate. We find the progenitor to have an absolute magnitude of $V = -4.63^{+0.3}_{-0.4}$ mag and a colour of $V-I = 2.29^{+0.25}_{-0.39}$ mag, corresponding to a progenitor luminosity of log $L$/L$_{\odot}$ $\sim4.54\pm0.19$ dex. Using the stellar evolution code STARS, we find this to be consistent with a red supergiant progenitor with $M = 8.5_{-1.5}^{+6.5}$ M$_{\odot}$. The photometric and spectroscopic evolution of SN 2009md is similar to that of the class of sub-luminous Type IIP SNe; in this paper we compare the evolution of SN 2009md primarily to that of the sub-luminous SN 2005cs. We estimate the mass of $^{56}$Ni ejected in the explosion to be $(5.4\pm1.3) \times 10^{-3}$ M$_{\odot}$\ from the luminosity on the radioactive tail, which is in agreement with the low $^{56}$Ni masses estimated for other sub-luminous Type IIP SNe. From the lightcurve and spectra, we show the SN explosion had a lower energy and ejecta mass than the normal Type IIP SN 1999em. We discuss problems with stellar evolutionary models, and the discrepancy between low observed progenitor luminosities (log $L$/L$_{\odot}$ $\sim4.3-5$ dex) and model luminosities after the second-dredge-up for stars in this mass range, and consider an enhanced carbon burning rate as a possible solution. In conclusion, SN 2009md is a faint SN arising from the collapse of a progenitor close to the lower mass limit for core-collapse. This is now the third discovery of a low mass progenitor star producing a low energy explosion and low $^{56}$Ni ejected mass, which indicates that such events arise from the lowest end of the mass range that produces a core-collapse supernova (7-8 M$_{\odot}$).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the nature of the progenitors of three type II-P supernovae: 2004et, 2006my and 2006ov

The pre-explosion observations of the type II-P supernovae 2006my, 2006ov and 2004et, are re-analysed. In the cases of supernovae 2006my and 2006ov we argue that the published candidate progenitors are not coincident with their respective supernova sites in pre-explosion Hubble Space Telescope observations. We therefore derive upper luminosity and mass limits for the unseen progenitors of both these supernovae, assuming they are red supergiants: 2006my (log L/Lsun = 4.51; mass < 13Msun) and 2006ov (log L/Lsun = 4.29; mass < 10Msun). In the case of supernova 2004et we show that the yellow-supergiant progenitor candidate, originally identified in Canada France Hawaii Telescope images, is still visible ~3 years post-explosion in observations from the William Herschel Telescope. High-resolution Hubble Space Telescope and Gemini (North) adaptive optics late-time imagery reveal that this source is not a single yellow supergiant star, but rather is resolved into at least three distinct sources. We report the discovery of the unresolved progenitor as an excess of flux in pre-explosion Isaac Newton Telescope i'-band imaging. Accounting for the late-time contribution of the supernova using published optical spectra, we calculate the progenitor photometry as the difference between the pre- and post-explosion, ground-based observations. We find the progenitor was most likely a late K to late M-type supergiant of 8 +5/-1 Msun. In all cases we conclude that future, high-resolution observations of the supernova sites will be required to confirm these results.

0 citations0 reviewsNo code linkedOpen access