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, $\chi_C$, between $0.002 - 40\%$. A large ($0.1 - 1.0$ $\mu$m) dust grain size composition is favored for efficient dustars where $\chi_C\gtrsim1\%$. Results for dustars with known orbital periods verify a power-law relation between $\chi_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).

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.

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.

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.

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

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.

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.

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.16}$ 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.

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.

We show that it is not possible to determine the final mass $M_{\rm fin}$ of a red supergiant (RSG) at the pre-supernova (SN) stage from its luminosity $L$ and effective temperature $T_{\rm eff}$ alone. Using a grid of stellar models, we demonstrate that for a given value of $L$ and $T_{\rm eff}$, a RSG can have a range of $M_{\rm fin}$ as wide as 3 to $45~\mathrm{M}_{\odot}$. While the probability distribution within these limits is not flat, any individual determination of $M_{\rm fin}$ for a RSG will be degenerate. This makes it difficult to determine its evolutionary history and to map $M_{\rm fin}$ to an initial mass. Single stars produce a narrower range that is difficult to accurately determine without making strong assumptions about mass loss, convection, and rotation. Binaries would produce a wider range of RSG $M_{\rm fin}$. However, the final Helium core mass M$_{\rm He-core}$ is well determined by the final luminosity and we find $\log (\mathrm{M}_{\rm He-core}/M_{\odot}) = 0.659 \log (L/\mathrm{L}_{\odot}) -2.630$ Using this relationship, we derive M$_{\rm He-core}$ for directly imaged SN progenitors and one failed SN candidate. The value of $M_{\rm fin}$ for stripped star progenitors of SNe IIb is better constrained by $L$ and $T_{\rm eff}$ due to the dependence of $T_{\rm eff}$ on the envelope mass $M_{\rm env}$ for $M_{\rm env} \lesssim 1~$M$_{\odot}$. Given the initial mass function, our results apply to the majority of progenitors of core collapse SNe, failed SNe and direct collapse black holes.

There is observational evidence of a dearth in core-collapse supernova (ccSN) explosions from stars with zero age main sequence (ZAMS) mass M_0~17-30 Msol, referred to as the 'red supergiant problem'. However, simulations now predict that above 20Msol we should indeed only expect stars within certain pockets of M_0 to produce a visible SN explosion. Validating these predictions requires large numbers of ccSNe of different types with measured M_0, which is challenging. In this paper we explore the reliability of using host galaxy emission lines and the Halpha equivalent width to constrain the age, and thus the M_0 of ccSNe progenitors. We use Binary Population and Spectral Synthesis models to infer a stellar population age from MUSE observations of the ionised gas properties and Halpha EW at the location of eleven ccSNe with reliable M_0 measurements. Comparing our results to published M_0 values, we find that models that do not consider binary systems yield stellar ages that are systematically too young (thus M_0 too large), whereas accounting for binary system interactions typically overpredict the stellar age (thus underpredict M_0). Taking into account the effects of photon leakage bring our M_0 estimates in much closer agreement with expectations. These results highlight the need for careful modelling of diffuse environments, such as are present in the vicinity of type II SNe, before ionised emission line spectra can be used as reliable tracers of progenitor stellar age.

Observations of star-forming galaxies in the distant Universe have confirmed the importance of massive stars in shaping galaxy emission and evolution. Distant stellar populations are unresolved, and the limited data available must be interpreted in the context of stellar population models. Understanding these populations, and their evolution with age and heavy element content is key to interpreting processes such as supernovae, cosmic reionization and the chemical enrichment of the Universe. With the upcoming launch of JWST and observations of galaxies within a billion years of the Big Bang, the uncertainties in modelling massive stars - particularly their interactions with binary companions - are becoming increasingly important to our interpretation of the high redshift Universe. In turn, observations of distant stellar populations provide ever stronger tests against which to gauge the success of, and flaws in, current massive star models. Here we briefly review the current status binary stellar population synthesis.

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.

Long-duration gamma-ray bursts (GRBs) are understood to be the final fate for a subset of massive, stripped envelope, rapidly rotating stars. Beyond this, our knowledge of the progenitor systems is limited. Using the BPASS (Binary Population and Spectral Synthesis) stellar evolution models, we investigate the possibility that some massive stars in binaries can maintain the angular momentum required for jet production, while still loosing their outer envelope through winds or binary interactions. We find that a total hydrogen mass of less than 0.0005 Msun and a helium ejecta mass fraction of less than 0.20 provide the best thresholds for the supernova type II/Ibc and Ib/Ic divisions respectively. Tidal interactions in binaries are accounted for by applying a tidal algorithm to post-process the stellar evolution models output by BPASS. We show that the observed volumetric gamma-ray burst rate evolution can be recreated using two distinct pathways and plausible distributions for burst parameters. In the first pathway, stars are spun up by mass accretion into a quasi-homogeneous state. In the second, tides maintain rotation where otherwise the star would spin down. Both lead to type Ic supernova progenitors, and a metallicity distribution consistent with the GRB host galaxy population. The inferred core angular momentum threshold for jet production is consistent with theoretical requirements for collapsars, given the assumptions made in our model. We can therefore reproduce several aspects of core collapse supernova/GRB observation and theory simultaneously. We discuss the predicted observable properties of GRB progenitors and their surviving companions.

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.

We find that applying a theoretical wind mass-loss rate from Monte Carlo radiative transfer models for hydrogen-deficient stars results in significantly more leftover hydrogen following stable mass transfer through Roche-lobe overflow than when we use an extrapolation of an empirical fit for Galactic Wolf-Rayet stars, for which a negligible amount of hydrogen remains in a large set of binary stellar evolution computations. These findings have implications for modelling progenitors of Type Ib and Type IIb supernovae. Most importantly, our study stresses the sensitivity of the stellar evolution models to the assumed mass-loss rates and the need to develop a better theoretical understanding of stellar winds.

The distance to the fireworks galaxy NGC 6946 is highly uncertain. Recent distance estimates using the tip of the red giant branch of 7.7 to 7.8 Mpc are larger than the distance commonly assumed for studying supernovae in this galaxy of 5.9 Mpc. Here we use the high supernova rate of the galaxy to derive the star-formation rate and predict the galaxy's FUV flux. We also account for dust extinction by different methods to derive a distance of 7.9+/-4.0 Mpc for NGC 6946. We then use the new distance to re-evaluate the understanding of the supernova progenitors 2002hh, 2004et, 2017eaw, the possible electron capture event 2008S and the possible black-hole forming event N6946-BH1. For the latter two exotic events the new distance improves the consistency between the observed progenitors and the stellar models that give rise to these events. From our findings we strongly recommend that all future studies of NGC~6946 must use the greater distance to the galaxy of 7.72+/-0.32 Mpc of Anand et al. (2018).

Observations of both galaxies in the distant Universe and local starbursts are showing increasing evidence for very hard ionizing spectra that stellar population synthesis models struggle to reproduce. Here we explore the effects of the assumed stellar initial mass function (IMF) on the ionizing photon output of young populations at wavelengths below key ionization energy thresholds. We use a custom set of binary population and spectral synthesis (BPASS) models to explore the effects of IMF assumptions as a function of metallicity, IMF slope, upper mass limit, IMF power law break mass and sampling. We find that while the flux capable of ionizing hydrogen is only weakly dependent on IMF parameters, the photon flux responsible for the He II and O VI lines is far more sensitive to assumptions. In our current models this flux arises primarily from helium and Wolf-Rayet stars which have partially or fully lost their hydrogen envelopes. The timescales for formation and evolution of both Wolf Rayet stars and helium dwarfs, and hence inferred population age, are affected by choice of model IMF. Even the most extreme IMFs cannot reproduce the He II ionizing flux observed in some high redshift galaxies, suggesting a source other than stellar photospheres. We caution that detailed interpretation of features in an individual galaxy spectrum is inevitably going to be subject to uncertainties in the IMF of its contributing starbursts. We remind the community that the initial mass function is fundamentally a statistical construct, and that stellar population synthesis models are most effective when considering entire galaxy populations rather than individual objects.

We present results of a supernova light-curve population synthesis, predicting the range of possible supernova lightcurves arising from a population of progenitor stars that include interacting binary systems. We show that the known diversity of supernova lightcurves can be interpreted as arising from binary interactions. Given detailed models of the progenitor stars, we are able to the determine what parameters within these stars determine the shape of their supernova lightcurve. The primary factors are the mass of supernova ejecta and the mass of hydrogen in the final progenitor. We find that there is a continuum of light-curve behaviour from type IIP, IIL to IIb supernovae related to the range of hydrogen and ejecta masses. Most type IIb supernovae arise from a relatively narrow range of initial masses from 10 to 15 M$_{\odot}$. We also find a few distinct lightcurves that are the result of stellar mergers.

Gravitational wave transients, resulting from the merger of two stellar remnants, are now detectable. The properties and rates of these directly relates to the stellar population which gave rise to their progenitors, and thus to other, electromagnetic transients which result from stellar death. We aim to estimate simultaneously the event rates and delay time distribution of gravitational wave-driven compact object mergers together with the rates of core collapse and thermonuclear supernovae within a single consistent stellar population synthesis paradigm. We combine event delay-time distributions at different metallicities from the Binary Population and Spectral Synthesis (BPASS) models with an analytic model of the volume-averaged cosmic star formation rate density and chemical evolution to determine the volume-averaged rates of each event rate at the current time. We estimate rates in excellent agreement with extant observational constraints on core-collapse supernovae, thermonuclear supernovae and long GRBs. We predict rates for gravitational wave mergers based on the same stellar populations, and find rates consistent with current LIGO estimates. We note that tighter constraints on the rates of these events will be required before it is possible to determine their redshift evolution, progenitor metallicity dependence or constrain uncertain aspects of stellar evolution.

A large fraction of massive stars evolve in interacting binary systems, which dramatically modifies the outcome of stellar evolution. We investigated the properties of blue supergiants in binary systems and whether they are suitable for extragalactic distance determinations using the flux-weighted gravity luminosity relationship (FGLR). This is a relationship between the absolute bolometric magnitude $M_{\rm bol}$ and the spectroscopically determined flux-weighted gravity $g_{\rm F} = g/T^4_{\rm eff}$. We computed a grid of binary stellar evolution models with MESA and use the v2.1 BPASS models to examine whether they are compatible with the relatively small scatter shown by the observed relationship. Our models have initial primary masses of 9 - 30 $M_\odot$, initial orbital periods of 10 - 2511 days, mass ratio $q$ = 0.9, and metallicity $Z$ = 0.02. We find that the majority of primary stars that produce blue supergiant stages are consistent with the observed FGLR, with a small offset towards brighter bolometric magnitudes. In between 1\% - 24\% of cases, binary evolution may produce blue supergiants after a mass transfer episode, that lie below the observed FGLR. A small number of such stars have been found in extragalactic FGLR studies, suggesting that they may have evolved through binary interaction. Some models with shorter periods could resemble blue hypergiants and LBVs. We used CMFGEN radiative transfer models to investigate the effects of unresolved secondaries on diagnostics for $T_{\rm eff}$ and $g$, and the biases on the determination of interstellar reddening and $M_{\rm bol}$. We find that the effects are small and within the observed scatter, but could lead to a small overestimate of the luminosity, of $T_{\rm eff}$ and of $g$ for extreme cases. We conclude that the observed FGLR can, in principle, be well reproduced by close binary evolution models.

Massive stars are the main objects that illuminate H II regions and they evolve quickly to end their lives in core-collapse supernovae (CCSNe). Thus it is important to investigate the association between CCSNe and H II regions. In this paper, we present emission line diagnostics of the stellar populations around nearby CCSNe, that include their host H II regions, from the PMAS/PPAK Integral-field Supernova hosts COmpilation (PISCO). We then use BPASS stellar population models to determine the age, metallicity and gas parameters for H II regions associated with CCSNe, contrasting models that consider either single star evolution alone or incorporate interacting binaries. We find binary-star models, that allow for ionizing photon loss, provide a more realistic fit to the observed CCSN hosts with metallicities that are closer to those derived from the oxygen abundance in O3N2. We also find that type II and type Ibc SNe arise from progenitor stars of similar age, mostly from 7 to 45 Myr, which corresponds to stars with masses < 20 solar mass . However these two types SNe have little preference in their host environment metallicity measured by oxygen abundance or in progenitor initial mass. We note however that at lower metallicities supernovae are more likely to be of type II.

In Bray and Eldridge (2017), we proposed a simple neutron star kick formula, v kick = alpha (M ejecta / M remnant) + beta to explain the observed 2D velocities of young single neutron stars. Using this kick we found that there is no statistically significant preference for a kick orientation nor for any of the three initial mass function (IMF) slopes tested, and that populations including binary stars reproduced the kick distribution better than single star only populations. However, recent analysis by Janka (2017), prompted us to revisit our basic assumptions and our new analysis has led to revised best-fit kick values of alpha=100 km per second and beta = -170 km per second. The reduction of beta to a negative value is due to using the 2D observed kick velocity distribution rather than the modelled 3D velocity distribution for neutron stars (NS). To further test the validity of the new kick, we have created synthetic populations of runaway star and double neutron star (DNS) binaries at solar metallicity (Z=0.02) using our best-fit kick. We find our new kick values create runaway star velocities and DNS period distributions in agreement with the comparable observational distributions with only the DNS eccentricities in tension with the observations. From our DNS and BH-BH datasets we estimate a predicted DNS merger rate at solar metallicity of 3,864 (+1,570/-2,371) per cubic Gpc per yr and a BH-BH merger rate of 5 (+40/-1) per cubic Gpc per yr.

Determining the properties of old stellar populations (those with age >1 Gyr) has long involved the comparison of their integrated light, either in the form of photometry or spectroscopic indexes, with empirical or synthetic templates. Here we reevaluate the properties of old stellar populations using a new set of stellar population synthesis models, designed to incorporate the effects of binary stellar evolution pathways as a function of stellar mass and age. We find that single-aged stellar population models incorporating binary stars, as well as new stellar evolution and atmosphere models, can reproduce the colours and spectral indices observed in both globular clusters and quiescent galaxies. The best fitting model populations are often younger than those derived from older spectral synthesis models, and may also lie at slightly higher metallicities.

We present numerical models of the nebular emission from H II regions around young stellar populations over a range of compositions and ages. The synthetic stellar pop- ulations include both single stars and interacting binary stars. We compare these models to the observed emission lines of 254 HII regions of 13 nearby spiral galax- ies and 21 dwarf galaxies drawn from archival data. The models are created using the combination of the Binary Population and Spectral Synthesis (bpass) code with the photoionization code cloudy to study the differences caused by the inclusion of interacting binary stars in the stellar population. We obtain agreement with the observed emission line ratios from the nearby star-forming regions and discuss the effect of binary star evolution pathways on the nebular ionization of H II regions. We find that at population ages above 10 Myr, single-star models rapidly decrease in flux and ionization strength, while binary-star model still produce strong flux and high O III/H\b{eta} ratios. Our models can reproduce the metallicity of H II regions from spiral galaxies but we find higher metallicities than previously estimated for the H II regions from dwarf galaxies. Comparing the equivalent width of H\b{eta} emission between mod- els and observations, we find that accounting for ionizing photon leakage can affect age estimates for HII regions. When it is included, the typical age derived for HII regions is 5Myr from single-star models, and up to 10Myr with binary-star models. This is due to the existence of binary-star evolution pathways which produce more hot WR and helium stars at older ages. For future reference, we calculate new BPASS binary maximal starburst lines as a function of metallicity, and for the total model population, and present these in Appendix A.

The first directly detected gravitational waves (GW 150914) were emitted by two coalescing black holes (BHs) with masses of ~36Msun and ~29Msun. Several scenarios have been proposed to put this detection into an astrophysical context. The evolution of an isolated massive binary system is among commonly considered models. Various groups have performed detailed binary-evolution calculations that lead to BH merger events. However, the question remains open as to whether binary systems with the predicted properties really exist. The aim of this paper is to help observers to close this gap by providing spectral characteristics of massive binary BH progenitors during a phase where at least one of the companions is still non-degenerate. Stellar evolution models predict fundamental stellar parameters. Using these as input for our stellar atmosphere code (PoWR), we compute a set of models for selected evolutionary stages of massive merging BH progenitors at different metallicities. The synthetic spectra obtained from our atmosphere calculations reveal that progenitors of massive BH merger events start their lives as O2-3V stars that evolve to early-type blue supergiants before they undergo core-collapse during the Wolf-Rayet phase. When the primary has collapsed, the remaining system will appear as a wind-fed high-mass X-ray binary. We provide feedback parameters, broad band magnitudes, and spectral templates that should help to identify such binaries in the future. Comparisons of empirically determined mass-loss rates with those assumed by evolution calculations reveal significant differences. The consideration of the empirical mass-loss rates in evolution calculations will possibly entail a shift of the maximum in the predicted binary-BH merger rate to higher metallicities, that is, more candidates should be expected in our cosmic neighborhood than previously assumed.

The progenitor systems of the class of "Ca-rich transients" is a key open issue in time domain astrophysics. These intriguing objects exhibit unusually strong calcium line emissions months after explosion, fall within an intermediate luminosity range, are often found at large projected distances from their host galaxies, and may play a vital role in enriching galaxies and the intergalactic medium. Here we present multi-wavelength observations of iPTF15eqv in NGC 3430, which exhibits a unique combination of properties that bridge those observed in Ca-rich transients and Type Ib/c supernovae. iPTF15eqv has among the highest [Ca II]/[O I] emission line ratios observed to date, yet is more luminous and decays more slowly than other Ca-rich transients. Optical and near-infrared photometry and spectroscopy reveal signatures consistent with the supernova explosion of a < 10 solar mass star that was stripped of its H-rich envelope via binary interaction. Distinct chemical abundances and ejecta kinematics suggest that the core collapse occurred through electron capture processes. Deep limits on possible radio emission made with the Jansky Very Large Array imply a clean environment ($n <$ 0.1 cm$^{-3}$) within a radius of $\sim 10^{17}$ cm. Chandra X-ray Observatory observations rule out alternative scenarios involving tidal disruption of a white dwarf by a black hole, for masses > 100 solar masses). Our results challenge the notion that spectroscopically classified Ca-rich transients only originate from white dwarf progenitor systems, complicate the view that they are all associated with large ejection velocities, and indicate that their chemical abundances may vary widely between events.

The Binary Population and Spectral Synthesis (BPASS) suite of binary stellar evolution models and synthetic stellar populations provides a framework for the physically motivated analysis of both the integrated light from distant stellar populations and the detailed properties of those nearby. We present a new version 2.1 data release of these models, detailing the methodology by which BPASS incorporates binary mass transfer and its effect on stellar evolution pathways, as well as the construction of simple stellar populations. We demonstrate key tests of the latest BPASS model suite demonstrating its ability to reproduce the colours and derived properties of resolved stellar populations, including well- constrained eclipsing binaries. We consider observational constraints on the ratio of massive star types and the distribution of stellar remnant masses. We describe the identification of supernova progenitors in our models, and demonstrate a good agreement to the properties of observed progenitors. We also test our models against photometric and spectroscopic observations of unresolved stellar populations, both in the local and distant Universe, finding that binary models provide a self-consistent explanation for observed galaxy properties across a broad redshift range. Finally, we carefully describe the limitations of our models, and areas where we expect to see significant improvement in future versions.

Superluminous supernovae (SLSNe) are found predominantly in dwarf galaxies, indicating that their progenitors have a low metallicity. However, the most nearby SLSN to date, SN 2017egm, occurred in the spiral galaxy NGC 3191, which has a relatively high stellar mass and correspondingly high metallicity. In this paper, we present detailed analysis of the nearby environment of SN 2017egm using MaNGA IFU data, which provides spectral data on kiloparsec scales. From the velocity map we find no evidence that SN 2017egm occurred within some intervening satellite galaxy, and at the SN position most metallicity diagnostics yield a solar and above solar metallicity (12 + log (O/H) = 8.8-9.1). Additionally we measure a small H-alpha equivalent width (EW) at the SN position of just 34 Angs, which is one of the lowest EWs measured at any SLSN or Gamma-Ray Burst position, and indicative of the progenitor star being comparatively old. We also compare the observed properties of NGC 3191 with other SLSN host galaxies. The solar-metallicity environment at the position of SN 2017egm presents a challenge to our theoretical understanding, and our spatially resolved spectral analysis provides further constraints on the progenitors of SLSNe.

SN 2002cx-like Type Ia supernovae (also known as SNe Iax) represent one of the most numerous peculiar SN classes. They differ from normal SNe Ia by having fainter peak magnitudes, faster decline rates and lower photospheric velocities, displaying a wide diversity in these properties. We present both integral-field and long-slit visual-wavelength spectroscopy of the host galaxies and explosion sites of SNe Iax to provide constraints on their progenitor formation scenarios. The SN Iax explosion site metallicity distribution is similar to that of core-collapse (CC) SNe and metal-poor compared to normal SNe Ia. Fainter members, speculated to form distinctly from brighter SN Iax, are found at a range of metallicities, extending to very metal-poor environments. Although the SN Iax explosion sites' ages and star-formation rates are comparatively older and less intense than the distribution of star forming regions across their host galaxies, we confirm the presence of young stellar populations (SP) at explosion environments for most SNe Iax, expanded here to a larger sample. Ages of the young SP (several $\times 10^{7}$ to $10^8$~yrs) are consistent with predictions for young thermonuclear and electron-capture SN progenitors. The lack of extremely young SP at the explosion sites disfavours very massive progenitors such as Wolf-Rayet explosions with significant fall-back. We find weak ionised gas in the only SN Iax host without obvious signs of star-formation. The source of the ionisation remains ambiguous but appears unlikely to be mainly due to young, massive stars.

Our Galaxy hosts the annihilation of a few $\times 10^{43}$ low-energy positrons every second. Radioactive isotopes capable of supplying such positrons are synthesised in stars, stellar remnants, and supernovae. For decades, however, there has been no positive identification of a main stellar positron source leading to suggestions that many positrons originate from exotic sources like the Galaxy's central super-massive black hole or dark matter annihilation. %, but such sources would not explain the recently-detected positron signal from the extended Galactic disk. Here we show that a single type of transient source, deriving from stellar populations of age 3-6 Gyr and yielding ~0.03 $M_\odot$ of the positron emitter $^{44}$Ti, can simultaneously explain the strength and morphology of the Galactic positron annihilation signal and the solar system abundance of the $^{44}$Ti decay product $^{44}$Ca. This transient is likely the merger of two low-mass white dwarfs, observed in external galaxies as the sub-luminous, thermonuclear supernova known as SN1991bg-like.

We present radio observations at 1.5 GHz of 32 local objects selected to reproduce the physical properties of $z\sim5$ star-forming galaxies. We also report non-detections of five such sources in the sub-millimetre. We find a radio-derived star formation rate which is typically half that derived from H$\alpha$ emission for the same objects. These observations support previous indications that we are observing galaxies with a young dominant stellar population, which has not yet established a strong supernova-driven synchrotron continuum. We stress caution when applying star formation rate calibrations to stellar populations younger than 100 Myr. We calibrate the conversions for younger galaxies, which are dominated by a thermal radio emission component. We improve the size constraints for these sources, compared to previous unresolved ground-based optical observations. Their physical size limits indicate very high star formation rate surface densities, several orders of magnitude higher than the local galaxy population. In typical nearby galaxies, this would imply the presence of galaxy-wide winds. Given the young stellar populations, it is unclear whether a mechanism exists in our sources that can deposit sufficient kinetic energy into the interstellar medium to drive such outflows.

We present a new model of the optical nebular emission from HII regions by combin- ing the results of the Binary Population and Spectral Synthesis (bpass) code with the photoion- ization code cloudy (Ferland et al. 1998). We explore a variety of emission-line diagnostics of these star-forming HII regions and examine the effects of metallicity and interacting binary evo- lution on the nebula emission-line production. We compare the line emission properties of HII regions with model stellar populations, and provide new constraints on their stellar populations and supernova progenitors. We find that models including massive binary stars can successfully match all the observational constraints and provide reasonable age and mass estimation of the HII regions and supernova progenitors.

The z = 6.6 Lyman-$\alpha$ emitter 'CR7' has been claimed to have a Population III-like stellar population, or alternatively, be a candidate Direct Collapse Black Hole (DCBH). In this paper we investigate the evidence for these exotic scenarios using recently available, deeper, optical, near-infrared and mid-infrared imaging. We find strong Spitzer/IRAC detections for the main component of CR7 at 3.6 and 4.5 microns, and show that it has a blue colour ([3.6] - [4.5] $= -1.2\pm 0.3$). This colour cannot be reproduced by current Pop. III or pristine DCBH models. Instead, the results suggest that the [3.6] band is contaminated by the [OIII]4959,5007 emission line with an implied rest-frame equivalent width of EW_0 (H$\beta$ + [OIII]) $\gtrsim 2000$\AA. Furthermore, we find that new near-infrared data from the UltraVISTA survey supports a weaker He II 1640 emission line than previously measured, with EW_0 $= 40 \pm 30$\AA. For the fainter components of CR7 visible in Hubble Space Telescope imaging, we find no evidence that they are particularly red as previously claimed, and show that the derived masses and ages are considerably uncertain. In light of the likely detection of strong [OIII] emission in CR7 we discuss other more standard interpretations of the system that are consistent with the data. We find that a low-mass, narrow-line AGN can reproduce the observed features of CR7, including the lack of radio and X-ray detections. Alternatively, a young, low-metallicity (~1/200 solar) starburst, modelled including binary stellar pathways, can reproduce the inferred strength of the He II line and simultaneously the strength of the observed [OIII] emission, but only if the gas shows super-solar $\alpha$-element abundances (O/Fe ~ 5 O/Fe solar).

Most types of supernovae (SNe) have yet to be connected with their progenitor stellar systems. Here, we reanalyze the 10-year SN sample collected during 1998-2008 by the Lick Observatory Supernova Search (LOSS) in order to constrain the progenitors of SNe Ia and stripped-envelope SNe (SE SNe, i.e., SNe~IIb, Ib, Ic, and broad-lined Ic). We matched the LOSS galaxy sample with spectroscopy from the Sloan Digital Sky Survey and measured SN rates as a function of galaxy stellar mass, specific star formation rate, and oxygen abundance (metallicity). We find significant correlations between the SN rates and all three galaxy properties. The SN Ia correlations are consistent with other measurements, as well as with our previous explanation of these measurements in the form of a combination of the SN Ia delay-time distribution and the correlation between galaxy mass and age. The ratio between the SE SN and SN II rates declines significantly in low-mass galaxies. This rules out single stars as SE~SN progenitors, and is consistent with predictions from binary-system progenitor models. Using well-known galaxy scaling relations, any correlation between the rates and one of the galaxy properties examined here can be expressed as a correlation with the other two. These redundant correlations preclude us from establishing causality - that is, from ascertaining which of the galaxy properties (or their combination) is the physical driver for the difference between the SE SN and SN II rates. We outline several methods that have the potential to overcome this problem in future works.

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.

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.

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}=\alpha\,(M_{\rm ejecta} / M_{\rm remnant}) + \beta\,$. 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 $\alpha$ and $\beta$ 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 $\alpha=70\, {\rm km\,s^{-1}}$ and $\beta=110\,{\rm km\,s^{-1}}$ reproduce the \cite{RN165} observed 2-dimensional velocities and $\alpha=70\, {\rm km\,s^{-1}}$ and $\beta=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.

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.

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.

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.

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.

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.

Literature data are collated for 38 stripped-envelope core-collapse supernovae (SE SNe; i.e. SNe IIb, Ib, Ic and Ic-BL) that have good light curve coverage in more than one optical band. Using bolometric corrections derived in previous work, the bolometric light curve of each SN is recovered and template bolometric light curves provided. Peak light distributions and decay rates are investigated; SNe subtypes are not cleanly distinguished in this parameter space, although some grouping of types does occur and there is a suggestion of a Phillips-like relation for most SNe Ic-BL. The bolometric light curves are modelled with a simple analytical prescription and compared to results from more detailed modelling. Distributions of the explosion parameters shows the extreme nature of SNe Ic-BL in terms of their 56Ni mass and the kinetic energy, however ejected masses are similar to other subtypes. SNe Ib and Ic have very similar distributions of explosion parameters, indicating a similarity in progenitors. SNe~IIb are the most homogeneous subtype and have the lowest average values for 56Ni mass, ejected mass, and kinetic energy. Ejecta masses for each subtype and SE SNe as a whole are inconsistent with those expected from very massive stars. The majority of the ejecta mass distribution is well described by more moderately massive progenitors in binaries, indicating these are the dominant progenitor channel for SE SNe.

The relationship between stellar populations and the ionizing flux with which they irradiate their surroundings has profound implications for the evolution of the intergalactic medium. We quantify the ionizing flux arising from synthetic stellar populations which incorporate the evolution of interacting binary stars. We determine that these show ionizing flux boosted by 60 per cent at 0.05 < Z < 0.3 Z_sun and a more modest 10-20 per cent at near-Solar metallicities relative to star-forming populations in which stars evolve in isolation. The relation of ionizing flux to observables such as 1500A continuum and ultraviolet spectral slope is sensitive to attributes of the stellar population including age, star formation history and initial mass function. For a galaxy forming 1 M_sun yr^{-1}, observed at > 100 Myr after the onset of star formation, we predict a production rate of photons capable of ionizing hydrogen, N_ion = 1.4 x 10^{53} s^{-1} at Z = Z_sun and 3.5 x 10^{53} s^{-1} at 0.1 Z_sun, assuming a Salpeter-like initial mass function. We evaluate the impact of these issues on the ionization of the intergalactic medium, finding that the known galaxy populations can maintain the ionization state of the Universe back to z ~ 9, assuming that their luminosity functions continue to M_UV = -10, and that constraints on the intergalactic medium at z ~ 2 - 5 can be satisfied with modest Lyman continuum photon escape fractions of 4 - 24 per cent depending on assumed metallicity.

We present optical observations of supernova SN 2014C, which underwent an unprecedented slow metamorphosis from H-poor type Ib to H-rich type IIn over the course of one year. The observed spectroscopic evolution is consistent with the supernova having exploded in a cavity before encountering a massive shell of the progenitor star's stripped hydrogen envelope. Possible origins for the circumstellar shell include a brief Wolf-Rayet fast wind phase that overtook a slower red supergiant wind, eruptive ejection, or confinement of circumstellar material by external influences of neighboring stars. An extended high velocity Halpha absorption feature seen in near-maximum light spectra implies that the progenitor star was not completely stripped of hydrogen at the time of core collapse. Archival pre-explosion Subaru Telescope Suprime-Cam and Hubble Space Telescope Wide Field Planetary Camera 2 images of the region obtained in 2009 show a coincident source that is most likely a compact massive star cluster in NGC 7331 that hosted the progenitor system. By comparing the emission properties of the source with stellar population models that incorporate interacting binary stars we estimate the age of the host cluster to be 30 - 300 Myr, and favor ages closer to 30 Myr in light of relatively strong Halpha emission. SN 2014C is the best-observed member of a class of core-collapse supernovae that fill the gap between events that interact strongly with dense, nearby environments immediately after explosion and those that never show signs of interaction. Better understanding of the frequency and nature of this intermediate population can contribute valuable information about the poorly understood final stages of stellar evolution.

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.

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.

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.

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.

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.

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.

The progenitors of many type II core-collapse supernovae have now been identified directly on pre-discovery imaging. Here we present an extensive search for the progenitors of type Ibc supernovae in all available pre-discovery imaging since 1998. There are 12 type Ibc supernovae with no detections of progenitors in either deep ground-based or Hubble Space Telescope archival imaging. The deepest absolute BVR magnitude limits are between -4 and -5. We compare these limits with the observed Wolf-Rayet population in the Large Magellanic Cloud and estimate a 16 per cent probability we have failed to detect such a progenitor by chance. Alternatively the progenitors evolve significantly before core-collapse or we have underestimated the extinction towards the progenitors. Reviewing the relative rates and ejecta mass estimates from lightcurve modelling of Ibc SNe, we find both incompatible with Wolf-Rayet stars with initial masses >25Msun being the only progenitors. We present binary evolution models that fit these observational constraints. Stars in binaries with initial masses <20Msun lose their hydrogen envelopes in binary interactions to become low mass helium stars. They retain a low mass hydrogen envelope until approximately 10,000 years before core-collapse; hence it is not surprising that galactic analogues have been difficult to identify.

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]

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.

We present a new study of H{\alpha}/FUV flux ratios of star forming regions within a sample of nearby spiral galaxies. We search for evidence of the existence of a cluster mass dependent truncation in the underlying stellar initial mass function (IMF). We use an automated approach to identification of extended objects based on the SExtractor algorithm to catalogue resolved Hii regions within a set of nearby spiral galaxies. Corrections due to dust attenuation effects are applied to avoid artificially boosted H{\alpha}/FUV values. We use the BPASS stellar population synthesis code of Eldridge & Stanway (2009) to create a benchmark population of star forming regions to act as a reference for our observations. Based on those models, we identify a zone of parameter space populated by regions that cannot be obtained with a cluster mass dependent truncation in the stellar IMF imposed. We find that the investigated galaxies display small subpopulations of star forming regions falling within our zone of interest, which appears to be inconsistent with imposing an IMF truncation at maximum stellar mass dependent on the total cluster mass. The presence of those regions is expected in models using both stochastic and sorted sampling of the full extent of the stellar IMF. This result persists after taking dust attenuation effects into account. We highlight the significance of stochastic effects in environments with low star formation activity and in studies describing systems associated with small physical scales.

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.

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$\beta$ 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.

The existence of bright quasars at high redshifts implies that supermassive black holes were able to form in the early Universe. Though a number of mechanisms to achieve this have been proposed, none yet stands out. A recent suggestion is the formation of quasi-stars, initially stellar-mass black holes accreting from hydrostatic giant-like envelopes of gas, formed from the monolithic collapse of pre-galactic gas clouds. In this work, we modify the Cambridge STARS stellar evolution package to construct detailed models of the evolution of these objects. We find that, in all of our models, the black hole inside the envelope is able to reach slightly more than one-tenth of the total mass of the system before hydrostatic equilibrium breaks down. This breakdown occurs after a few million years of evolution. We show that the mechanism which causes the hydrostatic evolution to end is present in polytropic models. We also show that the solutions are highly sensitive to the size of the inner boundary radius and that no physical solutions exist if the inner boundary is chosen to be less than about 0.3 of the Bondi radius.

Using our Binary Population And Spectral Synthesis (BPASS) code we explore the effects on star-formation rate indicators of stochastically sampling the stellar initial mass function, adding a cluster mass dependent stellar upper-mass limit and including binary stars. We create synthetic spectra of young clusters and star-forming galaxies and compare these to observations of H(alpha) emission from isolated clusters and the relation between H(alpha) and FUV emission from nearby galaxies. We find that observations of clusters tend to favour a purely stochastic sampling of the initial mass function for clusters less than 100M(Sun), rather than the maximum stellar mass being dependant on the total cluster mass. It is more difficult to determine whether the same is true for more massive clusters. We also find that binary stars blur some of the observational differences that occur when a cluster-mass dependent stellar upper-mass limit is imposed when filling the IMF. The effect is greatest when modelling the observed H(alpha) and FUV star-formation rate ratios in galaxies. This is because mass transfer and merging of stars owing to binary evolution creates more massive stars and stars that have greater mass than the initial maximum imposed on the stellar population.

We investigate the red supergiant problem: the apparent dearth of Type IIP supernova progenitors with masses between 16 and 30 M_sun. Although red supergiants with masses in this range have been observed, none have been identified as progenitors in pre-explosion images. We show that by failing to take into account the additional extinction resulting from the dust produced in the red supergiant winds, we risk underestimating the luminosity of the most massive red supergiants at the end of their lives. We estimate the initial masses of all Type IIP progenitors for which observations exist and analyse the resulting population. We find that the most likely maximum mass for a Type IIP progenitor is 21^{+2}_{-1} M_sun. This is in closer agreement with the limit predicted from single star evolution models.

Chandra LETG+HRC-S and XMM-Newton RGS spectra of H-like C and N formed in the corona of the lambda And primary star, a mildly metal-poor G8 III-IV first ascent giant that completed dredge-up ~50 Myr ago, have been used to make a precise measurement of its surface C/N ratio. We obtain the formal result [C/N]=0.03+/-0.07, which is typical of old disk giants and in agreement with standard dredge-up theory for stars of about 1 M_sun or lower. In contrast, these stars as a group, including lambda And, have 12C/13C ~< 20, which is much lower than standard model predictions. We show that the abundances of the old disk giants are consistent with models including thermohaline mixing that begins at the red giant branch luminosity function "bump". Instead, lambda And indicates that the 12C/13C anomaly can be present immediately following dredge-up, contrary to current models of extra mixing on the red giant branch.

The effects of rapid rotation on stellar evolution can be profound. We are now beginning to gather enough data to allow a realistic comparison between different physical models. Two key tests for any theory of stellar rotation are first whether it can match observations of the enrichment of nitrogen, and potentially other elements, in clusters containing rapid rotators and secondly whether it can reproduce the observed broadening of the main sequence in the Hertzsprung-Russel diagram. Models of stellar rotation have been steadily increasing in number and complexity over the past two decades but the lack of data makes it difficult to determine whether such additions actually give a closer reflection of reality. One of the most poorly explored features of stellar rotation models is the treatment of angular momentum transport within convective zones. If we treat the core as having uniform specific angular momentum the angular momentum distribution in the star, for a given surface rotation, is dramatically different from what it is when we assume the star rotates as a solid body. The uniform specific angular momentum also generates strong shears which can drive additional transport of chemical elements close to the boundary of a convection zone. A comparison of different models and their reproduction of observable properties with otherwise identical input physics is essential to properly distinguish between them. We compare detailed grids of stellar evolution tracks of intermediate and high-mass stars produced using several models for rotation generated with our new stellar rotation code.

Young, massive stars dominate the rest-frame ultraviolet spectra of star-forming galaxies. At high redshifts (z>2), these rest-UV features are shifted into the observed-frame optical and a combination of gravitational lensing, deep spectroscopy and spectral stacking analysis allows the stellar population characteristics of these sources to be investigated. We use our stellar population synthesis code BPASS to fit two strong rest-UV spectral features in published Lyman-break galaxy spectra, taking into account the effects of binary evolution on the stellar spectrum. In particular, we consider the effects of quasi-homogeneous evolution (arising from the rotational mixing of rapidly-rotating stars), metallicity and the relative abundance of carbon and oxygen on the observed strengths of HeII (1640 Angstroms) and CIV (1548,1551 Angstroms) spectral lines. We find that Lyman-break galaxy spectra at z=2-3 are best fit with moderately sub-solar metallicities, and with a depleted carbon-to-oxygen ratio. We also find that the spectra of the lowest metallicity sources are best fit with model spectra in which the HeII emission line is boosted by the inclusion of the effect of massive stars being spun-up during binary mass-transfer so these rapidly-rotating stars experiencing quasi-homogeneous evolution.

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}$).

Mounting discoveries of extrasolar planets orbiting post-main sequence stars motivate studies aimed at understanding the fate of these planets. In the traditional "adiabatic" approximation, a secondary's eccentricity remains constant during stellar mass loss. Here, we remove this approximation, investigate the full two-body point-mass problem with isotropic mass loss, and illustrate the resulting dynamical evolution. The magnitude and duration of a star's mass loss combined with a secondary's initial orbital characteristics might provoke ejection, modest eccentricity pumping, or even circularisation of the orbit. We conclude that Oort clouds and wide-separation planets may be dynamically ejected from 1-7 Solar-mass parent stars during AGB evolution. The vast majority of planetary material which survives a supernova from a 7-20 Solar-mass progenitor will be dynamically ejected from the system, placing limits on the existence of first-generation pulsar planets. Planets around >20 Solar-mass black hole progenitors may easily survive or readily be ejected depending on the core collapse and superwind models applied. Material ejected during stellar evolution might contribute significantly to the free-floating planetary population.

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.

When a core collapse supernova occurs in a binary system, the surviving star as well as the compact remnant emerging from the SN, may reach a substantial space velocity. With binary population synthesis modelling at solar and one fifth of solar metallicity, we predict the velocities of such runaway stars or binaries. We compile predictions for runaway OB stars, red supergiants and Wolf-Rayet stars. For those stars or binaries which undergo a second stellar explosion we compute their further evolution and the distance travelled until a Type II or Type Ibc SN or a long or short gamma-ray burst occurs. We find our predicted population of OB runaway stars broadly matches the observed population of stars but, to match the fastest observed WR runaway stars, we require that black holes receive an asymmetric kick upon formation. We find that at solar metallicity Type Ic SN progenitors travel shorter distances than the progenitors of other SN types because they are typically more massive and thus have shorter lifetimes. Those of Type IIP SN can fly farthest about 48 pc on average at solar metallicity. In considering the consequences of assuming that the progenitors of long GRBs are spun-up secondary stars that experience quasi-homogeneous evolution, we find that such evolution has a dramatic effect on the population of runaway WR stars and that some 30 per cent of GRBs could occur a hundred parsecs or more from their initial positions. We also consider mergers of double compact object binaries consisting of neutron stars and/or black holes. We find the most common type of visible mergers are neutron star--black hole mergers that are roughly ten times more common than neutron star--neutron star mergers. We also find that there may be a population of low-velocity neutron stars that are ejected from a binary rather than by their own natal kick.

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.

We study the post-main sequence stars in NGC 604, the most luminous HII region in M33. Previously, a number of Wolf-Rayet (WR) stars and one red supergiant (RSG) have been discovered. Based on broadband photometry of the region, we present evidence that is consistent with the presence of this RSG and with three more RSG candidates. Using SED fitting based on HST UVIJHK photometry we estimate the ages of the WR stars and RSGs finding that the two populations are from distinct formation episodes with ages 3.2$\pm%1.0Myrs and 12.4$\pm$2.1Myrs, respectively. The RSGs have greater extinctions towards their line of sight than the WR stars consistent with the RSGs producing large amount of dust. Using the WR and RSG populations and similar SED fits to the most massive O stars we estimate that the total stellar mass is (3.8 $\pm$ 0.6) x 10^5Msun. We find a large discrepancy between the expected H{\alpha} flux from such a massive cluster and that one observed. This suggests that 49 (+16,-19) percent of the ionizing photons produced by massive stars in NGC 604 is leaking from this HII region. We also suggest that the implications of an old RSG population mean that if NGC 604 was more distant and only observed in the infrared (IR) it would be difficult to study the youngest burst of star formation due to the contamination of RSGs.

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.

Using our new Binary Population and Spectral Synthesis (BPASS) code we explore the affect of binary populations on the integrated spectra of galaxies. We also explore the interplay of binary populations and a varying maximum stellar mass. We compare our synthetic populations to observations of H$\alpha$ emission from isolated clusters and H$\alpha$ and FUV observations of galaxies. We find that observations tend to favour a pure stochastic sampling of the initial mass function although the evidence is not significant. We also find that binaries make a stellar population less susceptible to the stochastic effects of filling the IMF. Therefore making it more difficult to determine if there is a variable maximum stellar mass.

We derive a new age for the Gamma^2 Velorum binary by comparing recent observations to our set of binary models. We find that it is very unlikely the stars have not interacted, which implies that previous estimates of the age from single-star models of 3.5+/-0.4 Myrs are incorrect. We prefer an older age of 5.5+/-1 Myrs that agrees with the age of other lower mass stars in the Vela OB association. We also find that our favoured binary model shows the components of the binary have interacted in a Case B, post main-sequence, mass-transfer event. During the mass-transfer event, the envelopes of both components where radiative and therefore the damping of tidal forces are relatively weak. This explains why the binary is still eccentric after mass-transfer.

We have constructed a new code to produce synthetic spectra of stellar populations that includes massive binaries. We have tested this code against the broadband colours of unresolved young massive stellar clusters in nearby galaxies, the equivalent widths of the Red and Blue Wolf-Rayet bumps in star-forming SDSS galaxies and the UV and optical spectra of the star forming regions Tol-A and B in NGC5398. In each case we find a good agreement between our models and observations. We find that in general binary populations are bluer and have fewer red supergiants, and thus significantly less flux in the I-band and at longer wavelengths, than single star populations. Also we find that Wolf-Rayet stars occur over a wider range of ages up to 10^7 years in a stellar population including binaries, increasing the UV flux and Wolf-Rayet spectral features at later times. In addition we find that nebula emission contributes significantly to these observed properties and must be considered when comparing stellar models with observations of unresolved stellar populations. We conclude that incorporation of massive stellar binaries can improve the agreement between observations and synthetic spectral synthesis codes, particularly for systems with young stellar populations.

We have developed a detailed stellar evolution code capable of following the simultaneous evolution of both stars in a binary system, together with their orbital properties. To demonstrate the capabilities of the code we investigate potential progenitors for the Type IIb supernova 1993J, which is believed to have been an interacting binary system prior to its primary exploding. We use our detailed binary stellar evolution code to model this system to determine the possible range of primary and secondary masses that could have produced the observed characteristics of this system, with particular reference to the secondary. Using the luminosities and temperatures for both stars (as determined by Maund et al. 2004) and the remaining mass of the hydrogen envelope of the primary at the time of explosion, we find that if mass transfer is 100 per cent efficient the observations can be reproduced by a system consisting of a 15 solar mass primary and a 14 solar mass secondary in an orbit with an initial period of 2100 days. With a mass transfer efficiency of 50 per cent, a more massive system consisting of a 17 solar mass primary and a 16 solar mass secondary in an initial orbit of 2360 days is needed. We also investigate some of the uncertainties in the evolution, including the effects of tidal interaction, convective overshooting and thermohaline mixing.

We present the results of a 10.5 yr, volume limited (28 Mpc) search for supernova (SN) progenitor stars. We compile all SNe discovered within this volume (132, of which 27% are type Ia) and determine the relative rates of each sub-type from literature studies : II-P (59%), Ib/c (29%), IIb (5%), IIn (4%) and II-L (3%). Twenty II-P SNe have high quality optical or near-IR pre-explosion images that allow a meaningful search for the progenitor stars. In five cases they are clearly red supergiants, one case is unconstrained, two fall on compact coeval star clusters and the other twelve have no progenitor detected. We review and update all the available data for the host galaxies (distance, metallicity and extinction) and determine masses and upper mass estimates using the STARS stellar evolutionary code and a single consistent homogeneous method. A maximum likelihood calculation suggests that the minimum stellar mass for a type II-P to form is m(min)=8.5 +1/-1.5 Msol and the maximum mass for II-P progenitors is m(max)=16.5 +/- 1.5 Msol, assuming a Salpeter initial mass function (in the range Gamma = -1.35 +0.3/-0.7). The minimum mass is consistent with current estimates for white dwarf progenitor masses, but the maximum mass does not appear consistent with massive star populations. Red supergiants in the Local Group have masses up to 25Msol and the minimum mass to produce a Wolf-Rayet star in single star evolution (between solar and LMC metallicity) is similarly 25-30 Msol. We term this discrepancy the "red supergiant problem" and speculate that these stars could have core masses high enough to form black holes and SNe which are too faint to have been detected. Low luminosity SNe with low 56Ni production seem to arise from explosions of low mass progenitors near the mass threshold for core-collapse. (abridged).

We report the identification of a source coincident with the position of the nearby type II-P supernova (SN) 2008bk in high quality optical and near-infrared pre-explosion images from the ESO Very Large Telescope (VLT). The SN position in the optical and near-infrared pre-explosion images is identified to within about +-70 and +-40 mas, respectively, using post-explosion Ks-band images obtained with the NAOS CONICA adaptive optics system on the VLT. The pre-explosion source detected in four different bands is precisely coincident with SN 2008bk and is consistent with being dominated by a single point source. We determine the nature of the point source using the STARS stellar evolutionary models and find that its colours and luminosity are consistent with the source being a red supergiant progenitor of SN 2008bk with an initial mass of 8.5 +- 1.0 Msun.

A source coincident with the position of the type IIb supernova (SN) 2008ax is identified in pre-explosion Hubble Space Telescope (HST) Wide Field Planetary Camera 2 observations in three optical filters. We identify and constrain two possible progenitor systems: (i) a single massive star that lost most of its hydrogen envelope through radiatively driven mass loss processes, prior to exploding as a helium-rich Wolf-Rayet star with a residual hydrogen envelope, and (ii) an interacting binary in a low mass cluster producing a stripped progenitor. Late time, high resolution observations along with detailed modelling of the SN will be required to reveal the true nature of this progenitor star.

The study of the stars that explode as supernovae used to be a forensic study, working backwards from the remnants of the star. This changed in 1987 when the first progenitor star was identified in pre-explosion images. Currently there are 8 detected progenitors with another 21 non-detections, for which only a limit on the pre-explosion luminosity can be placed. This new avenue of supernova research has led to many interesting conclusions, most importantly that the progenitors of the most common supernovae, type IIP, are red supergiants as theory has long predicted. However no progenitors have been detected thus far for the hydrogen-free type Ib/c supernovae which, given the expected progenitors, is an unlikely result. Also observations have begun to show evidence that luminous blue variables, which are among the most massive stars, may directly explode as supernovae. These results contradict current stellar evolution theory. This suggests that we may need to update our understanding.

We present new spectroscopic and photometric data of the type Ibn supernovae 2006jc, 2000er and 2002ao. We discuss the general properties of this recently proposed supernova family, which also includes SN 1999cq. The early-time monitoring of SN 2000er traces the evolution of this class of objects during the first few days after the shock breakout. An overall similarity in the photometric and spectroscopic evolution is found among the members of this group, which would be unexpected if the energy in these core-collapse events was dominated by the interaction between supernova ejecta and circumstellar medium. Type Ibn supernovae appear to be rather normal type Ib/c supernova explosions which occur within a He-rich circumstellar environment. SNe Ibn are therefore likely produced by the explosion of Wolf-Rayet progenitors still embedded in the He-rich material lost by the star in recent mass-loss episodes, which resemble known luminous blue variable eruptions. The evolved Wolf-Rayet star could either result from the evolution of a very massive star or be the more evolved member of a massive binary system. We also suggest that there are a number of arguments in favour of a type Ibn classification for the historical SN 1885A (S-Andromedae), previously considered as an anomalous type Ia event with some resemblance to SN 1991bg.

We present near- and mid-infrared (IR) photometric data of the Type Ibn supernova (SN) 2006jc obtained with the United Kingdom Infrared Telescope (UKIRT), the Gemini North Telescope, and the Spitzer Space Telescope between days 86 and 493 post-explosion. We find that the IR behaviour of SN 2006jc can be explained as a combination of IR echoes from two manifestations of circumstellar material. The bulk of the near-IR emission arises from an IR echo from newly-condensed dust in a cool dense shell (CDS) produced by the interaction of the ejecta outward shock with a dense shell of circumstellar material ejected by the progenitor in a luminous blue variable (LBV) like outburst about two years prior to the SN explosion. The CDS dust mass reaches a modest 3.0 x 10^(-4) M(solar) by day 230. While dust condensation within a CDS formed behind the ejecta inward shock has been proposed before for one event (SN 1998S), SN 2006jc is the first one showing evidence for dust condensation in a CDS formed behind the ejecta outward shock in the circumstellar material. At later epochs, a substantial and growing contribution to the IR fluxes arises from an IR echo from pre-existing dust in the progenitor wind. The mass of the pre-existing CSM dust is at least ~8 x 10^(-3) M(solar). This work therefore adds to the evidence that mass-loss from the progenitors of core-collapse supernovae could be a major source of dust in the universe. However, yet again, we see no direct evidence that the explosion of a supernova produces anything other than a very modest amount of dust.

We report our attempts to locate the progenitor of the peculiar type Ic SN 2007gr in HST pre-explosion images of the host galaxy, NGC 1058. Aligning adaptive optics Altair/NIRI imaging of SN 2007gr from the Gemini (North) Telescope with the pre-explosion HST WFPC2 images, we identify the SN position on the HST frames with an accuracy of 20 mas. Although nothing is detected at the SN position we show that it lies on the edge of a bright source, 134+/-23 mas (6.9 pc) from its nominal centre. Based on its luminosity we suggest that this object is possibly an unresolved, compact and coeval cluster and that the SN progenitor was a cluster member, although we note that model profile fitting favours a single bright star. We find two solutions for the age of this assumed cluster; 7-/+0.5 Myrs and 20-30 Myrs, with turn-off masses of 28+/-4 Msun and 12-9 Msun respectively. Pre-explosion ground-based K-band images marginally favour the younger cluster age/higher turn-off mass. Assuming the SN progenitor was a cluster member, the turn-off mass provides the best estimate for its initial mass. More detailed observations, after the SN has faded, should determine if the progenitor was indeed part of a cluster, and if so allow an age estimate to within ~2 Myrs thereby favouring either a high mass single star or lower mass interacting binary progenitor.

We compare our latest single and binary stellar model results from the Cambridge STARS code to several sets of observations. We examine four stellar population ratios, the number of blue to red supergiants, the number of Wolf-Rayet stars to O supergiants, the number of red supergiants to Wolf-Rayet stars and the relative number of Wolf-Rayet subtypes, WC to WN stars. These four ratios provide a quantitative measure of nuclear burning lifetimes and the importance of mass loss during various stages of the stars' lifetimes. In addition we compare our models to the relative rate of type Ib/c to type II supernovae to measure the amount of mass lost over the entire lives of all stars. We find reasonable agreement between the observationally inferred values and our predicted values by mixing single and binary star populations. However there is evidence that extra mass loss is required to improve the agreement further, to reduce the number of red supergiants and increase the number of Wolf-Rayet stars.

We first discuss the mass range of type IIP SN progenitors and how the upper and lower limits impose interesting constraints on stellar evolution. Then we discuss the possible implications of two SNe, 2002ap and 2006jc, for Wolf-Rayet star mass-loss rates and long Gamma-ray bursts.

We present an analysis of the faint M star population seen as foreground contaminants in deep extragalactic surveys. We use space-based data to separate such stars from high redshift galaxies in a publically-available dataset, and consider the photometric properties of the resulting sample in the optical and infrared. The inferred distances place these stars well beyond the scale height of the thick disk. We find strong similarities between this faint sample (reaching i'_{AB}=25) and the brighter disk M dwarf population studied by other authors. The optical-infrared properties of the bulk of our sources spanning 6000A-4.5microns are consistent with those 5-10 magnitudes brighter. We also present deep spectroscopy of faint M dwarf stars reaching continuum limits of i'_{AB}~26, and measure absorption line strengths in the CaH2 and TiO5 bands. Both photometrically and spectroscopically, our sources are consistent with metallicities as low as a tenth solar: metal-rich compared with halo stars at similar heliocentric distances. We comment on the possible MACHO identification of M stars at faint magnitudes.

(Abridged) We present a search for the progenitor star of the Type Ic Supernova 2002ap in deep, high quality pre-explosion observations taken with the Canada-France-Hawaii Telescope (CFHT). Aligning high-resolution Hubble Space Telescope (HST) observations of the supernova itself with the archival CFHT images allowed us to pinpoint the location of the progenitor site on the ground based observations. We find that a source visible in the B and R band pre-explosion images close to the position of the SN is (1) not coincident with the SN position within the uncertainties of our relative astrometry, and (2) is still visible ~ 4.7 yrs post-explosion in late-time observations taken with the William Herschel Telescope. We therefore conclude that it is not the progenitor of SN 2002ap. Comparing our luminosity limits with stellar models of single stars at appropriate metallicity (Z=0.008) we conclude that any single star progenitor must have experienced at least twice the standard mass loss rates during pre-Wolf-Rayet evolution, been initially > 30-40M(Sun) and exploded as a Wolf-Rayet star of final mass 10-12M(Sun). Alternatively an initially less massive progenitor may have evolved in an interacting binary system. We constrain any possible binary companion to a main sequence star of < 20M(Sun), a neutron star or a black hole. By combining the pre-explosion limits with the ejecta mass estimates and constraints from X-ray and radio observations we conclude that any binary interaction most likely occurred as Case B mass transfer, either with or without a subsequent common envelope evolution phase.

The death of massive stars produces a variety of supernovae, which are linked to the structure of the exploding stars. The detection of several precursor stars of Type II supernovae have been reported, however we do not yet have direct information on the progenitors of the hydrogen deficient Type Ib and Ic supernovae. Here we report that the peculiar Type Ib supernova SN2006jc is spatially coincident with a bright optical transient that occurred in 2004. Spectroscopic and photometric monitoring of the supernova leads us to suggest that the progenitor was a carbon-oxygen Wolf-Rayet star embedded within a helium-rich circumstellar medium. There are different possible explanations for this pre-explosion transient. It appears similar to the giant outbursts of Luminous Blue Variables (LBV) of 60-100 solar mass stars, however the progenitor of SN2006jc was helium and hydrogen deficient. An LBV-like outburst of a Wolf-Rayet star could be invoked, but this would be the first observational evidence of such a phenomenon. Alternatively a massive binary system composed of an LBV which erupted in 2004, and a Wolf-Rayet star exploding as SN2006jc, could explain the observations.

Recent observations of Wolf-Rayet (WR) binaries WR151 and WR155 infer that their stellar winds are asymmetric. We show that such asymmetries can alter the stellar-wind bubble structure, bringing the wind-termination shock closer to the WR star. If the wind asymmetry is caused by rotation, the wind density and distance to the wind-termination shock are both decreased along the rotation axis by a factor of a few for the observed equator-to-pole wind density ratio of WR151. If this asymmetry lasts until core-collapse the time taken to reach the wind-termination shock by supernova ejecta or a gamma-ray burst jet is reduced. This leads to a distorted structure of the supernova ejecta and makes it more likely a constant density environment is inferred from gamma-ray burst afterglow observations.

We calculate the predicted UBVRIJHK absolute magnitudes for models of supernova progenitors and apply the result to the case of supernova 2005cs. We agree with previous results that the initial mass of the star was of low, around 6 to 8 M(sun). However such stars are thought to go through second dredge-up to become AGB stars. We show that had this occurred to the progenitor of 2005cs it would have been observed in JHK pre-explosion images. The progenitor was not detected in these bands and therefore we conclude that it was not an AGB star. Furthermore if some AGB stars do produce supernovae they will have a clear signature in pre-explosion near-infrared images. Electron-capture supernovae are thought to occur in AGB stars, hence the implication is that 2005cs was not an electron-capture supernova but was the collapse of an iron core.

Binary stars are important for a full understanding of stellar evolution. We present a summary of how predictions of the relative supernova rates varies between single and binary stars. We also show how the parameter space of different supernova types differs between single and binary stars. We then consider an important question of how to infer a supernova progenitor's properties from pre-explosion imaging and present rescent work of producing synthe tic colours for our stellar models to make a direct comparison with any detections or limits obtained on supernova progentiors from pre-explosion images.

If Wolf-Rayet stars are the progenitors of Gamma-ray bursts (GRBs), they must rotate rapidly to produce the GRB. This rotation may effect their stellar-wind bubbles and possibly explain why so many GRB afterglows occur in a constant density medium.

Aims: Recent theoretical predictions for the winds of Wolf-Rayet stars indicate that their mass-loss rates scale with the initial stellar metallicity in the local Universe.We aim to investigate how this predicted dependence affects the models of Wolf-Rayet stars and their progeny in different chemical environments. Methods: We compute models of stellar structure and evolution for Wolf-Rayet stars for different initial metallicities, and investigate how the scaling of the Wolf-Rayet mass-loss rates affects the final masses, the lifetimes of the WN and WC subtypes, and how the ratio of the two populations vary with metallicity. Results: We find significant effects of metallicity dependent mass-loss rates for Wolf-Rayet stars. For models that include the scaling of the mass-loss rate with initial metallicity, all WR stars become neutron stars rather than black holes at twice the solar metallicity; at lower $Z$, black holes have larger masses. We also show that our models that include the mass-loss metallicity scaling closely reproduce the observed decrease of the relative population of WC over WN stars at low metallicities.

We study the evolution of the circumstellar medium of massive stars. We pay particular attention to Wolf-Rayet stars that are thought to be the progenitors of some long Gamma-Ray Bursts. We detail the mass-loss rates we use in our stellar evolution models and how we estimate the stellar wind speeds during different phases. With these details we simulate the interactions between the wind and the interstellar medium to predict the circumstellar environment around the stars at the time of core-collapse. We then investigate how the structure of the environment might affect the GRB afterglow. We find that when the afterglow jet encounters the free-wind to stalled-wind interface that rebrightening occurs and a bump is seen in the afterglow light curve. However our predicted positions of this interface are too distant from the site of the GRB to reach while the afterglow remains observable. The values of the final-wind density, A_{*}, from our stellar models are similar to the values inferred from observed afterglow lightcurves and those from observed Wolf-Rayet stars. However we do not reproduce the lowest observed A_{*} values below 0.3. For these cases we suggest that the progenitors could have been a WO type Wolf-Rayet star, be in a close binary or very low metallicity star. Finally we turn our attention to the matter of stellar wind material producing absorption lines in the afterglow spectra. We discuss the observational signatures of two Wolf-Rayet stellar types, WC and WO, in the afterglow lightcurve and spectra. We also indicate how it may be possible to constrain the initial mass and metallicity of a GRB progenitor by using the inferred wind density and wind velocity.

The progenitors of core-collapse supernovae are stars with an initial mass greater than about 8M(sun). Understanding the evolution of these stars is necessary to comprehend the evolution and differences between supernovae. We have constructed new and unique opacity tables to increase model accuracy during the latest stages of stellar evolution. We have investigated how initial mass, initial composition and mass loss affects the progenitors and their populations. There are many prescriptions for mass loss. Different research groups use their preferred rates. We have compared 12 different prescriptions and determined which provides the best fit to observations. We use our preferred mass-loss scheme to make suggestions as to the source of the differences between supernova types from our progenitor models. Binary evolution is considered in order to search for low luminosity SN progenitors and progenitor types not possible from single stars. Removal of the hydrogen envelope is more common and we find quite different hydrogen deficient SN progenitors. We discuss the implications of our binary models for ultra-luminous X-ray sources and gamma-ray bursts. We present an estimation of the mass distribution for black holes at various metallcities showing that massive black holes are not formed until very low metallicities. Finally we combine the single star and binary results to determine their relative populations and compare to observations. However it is not possible to draw many firm conclusions because of the uncertainty in observations to date.

We discuss the mass ranges over which we find AGB and super-AGB stars. The most massive super-AGB stars are candidate progenitors for type II core-collapse SNe. We discuss the two supernovae, SN1980K and SN2003gd that provide some restrictions on the upper mass limit of super-AGB stars.

We present details of our investigation of the progenitors to core-collapse supernovae. We discuss observations and the theory of the lowest-mass stars to explode as supernovae.

We present maps of the nature of single star progenitors of supernovae and their remnants in mass and metallicity space. We find our results are similar to others but we have gone further in varying the amount of mixing and using various mass-loss schemes to see how the maps change. We find that extra-mixing, in the form of convective overshooting, moves boundaries such as the minimum mass for a supernova or WR star to lower masses. We also find that the pre-WR mass-loss determines the shape of our maps. We find that different mass-loss rates lead to quite different results. We find that the rise in luminosity at 2nd dredge-up places quite tight constraints on the masses of some progenitors and in particular the progenitor of supernova 2003gd.

We have included opacity tables in our stellar evolution code that enable us to accurately model the structure of stars composed of mixtures with carbon and oxygen independently enhanced relative to solar. We present tests to demonstrate the effects of the new tables. Two of these are practical examples, the effect on the evolution of a thermally pulsing asymptotic giant branch star and a Wolf-Rayet Star. The changes are small but perceptible.

We have developed a numerical code with which we study the effects of 2D perturbations on stellar structure. We present new numerical and analytical results on the heating of a main-sequence star in a binary system by its companion.