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.

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.

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.

SN 1993J is one of the best studied Type IIb supernovae. Spectropolarimetric data analyses were published over two decades ago at a time when the field of supernova spectropolarimetry was in its infancy. Here we present a new analysis of the spectropolarimetric data of SN 1993J and an improved estimate of its interstellar polarization (ISP) as well as a critical review of ISP removal techniques employed in the field. The polarization of SN 1993J is found to show significant alignment on the $q-u$ plane, suggesting the presence of a dominant axis and therefore of continuum polarization. We also see strong line polarization features, including $\mathrm{H\beta}$, He\,{\sc i} $\lambda 5876$, $\mathrm{H\alpha}$, He\,{\sc i} $\lambda 6678$, He\,{\sc i} $\lambda 7065$, and high velocity (HV) components of He\,{\sc i} $\lambda 5876$ and $\mathrm{H\alpha}$. SN 1993J is therefore the second example of a stripped envelope supernova, alongside iPTF13bvn, with prominent HV helium polarization features, and the first to show a likely HV \halpha contribution. Overall, we determine that the observed features can be interpreted as the superposition of anisotropically distributed line forming regions over ellipsoidal ejecta. We cannot exclude the possibility of an off-axis energy source within the ejecta. These data demonstrate the rich structures that are inaccessible if solely considering the flux spectra but can be probed by spectropolarimetric observations. In future studies, the new ISP corrected data can be used in conjunction with 3D radiative transfer models to better map the geometry of the ejecta of SN 1993J.

We present optical spectroscopy and spectropolarimetry of the Type IIb SN 2008aq 16 days and 27 days post-explosion. The spectrum of SN 2008aq remained dominated by Halpha P Cygni profile at both epochs, but showed a significant increase in the strength of the helium features, which is characteristic of the transition undergone by supernovae between Type IIb and Type Ib. Comparison of the spectra of SN 2008aq to other Type IIb SNe (SN 1993J, SN 2011dh, and SN 2008ax) at similar epochs revealed that the helium lines in SN 2008aq are much weaker, suggesting that its progenitor was stripped to a lesser degree. SN 2008aq also showed significant levels of continuum polarisation at p_cont = 0.70 (+/- 0.22) % in the first epoch, increasing to p_cont = 1.21 (+/- 0.33) % by the second epoch. Moreover, the presence of loops in the q-u planes of Halpha and He I in the second epoch suggests a departure from axial symmetry.

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.

Fundamental coding and software development skills are increasingly necessary for success in nearly every aspect of astronomical and astrophysical research as large surveys and high resolution simulations become the norm. However, professional training in these skills is inaccessible or impractical for many members of our community. Students and professionals alike have been expected to acquire these skills on their own, apart from formal classroom curriculum or on-the-job training. Despite the recognized importance of these skills, there is little opportunity to develop them - even for interested researchers. To ensure a workforce capable of taking advantage of the computational resources and the large volumes of data coming in the next decade, we must identify and support ways to make software development training widely accessible to community members, regardless of affiliation or career level. To develop and sustain a technology capable astronomical and astrophysical workforce, we recommend that agencies make funding and other resources available in order to encourage, support and, in some cases, require progress on necessary training, infrastructure and policies. In this white paper, we focus on recommendations for how funding agencies can lead in the promotion of activities to support the astronomy and astrophysical workforce in the 2020s.

Spectropolarimetry enables us to measure the geometry and chemical structure of the ejecta in supernova explosions, which is fundamental for the understanding of their explosion mechanism(s) and progenitor systems. We collected archival data of 35 Type Ia Supernovae (SNe Ia), observed with FORS on the Very Large Telescope at 127 epochs in total. We examined the polarization of the Si II $\lambda$6355 $\AA$ line (p$_{\rm Si II}$) as a function of time which is seen to peak at a range of various polarization degrees and epochs relative to maximum brightness. We reproduced the $\Delta$m$_{15}$-p$_{\rm Si II}$ relationship identified in a previous study, and show that subluminous and transitional objects display polarization values below the $\Delta$m$_{15}$-p$_{\rm Si II}$ relationship for normal SNe Ia. We found a statistically significant linear relationship between the polarization of the Si II $\lambda$6355 $\AA$ line before maximum brightness and the Si II line velocity and suggest that this, along with the $\Delta$m$_{15}$-p$_{\rm Si II}$ relationship, may be explained in the context of a delayed-detonation model. In contrast, we compared our observations to numerical predictions in the $\Delta$m$_{15}$-v$_{\rm Si II}$ plane and found a dichotomy in the polarization properties between Chandrasekhar and sub-Chandrasekhar mass explosions, which supports the possibility of two distinct explosion mechanisms. A subsample of SNe display evolution of loops in the $q$-$u$ plane that suggests a more complex Si structure with depth. This insight, which could not be gleaned from total flux spectra, presents a new constraint on explosion models. Finally, we compared our statistical sample of the Si II polarization to quantitative predictions of the polarization levels for the double-detonation, delayed-detonation, and violent-merger models.

Software is a critical part of modern research, and yet there are insufficient mechanisms in the scholarly ecosystem to acknowledge, cite, and measure the impact of research software. The majority of academic fields rely on a one-dimensional credit model whereby academic articles (and their associated citations) are the dominant factor in the success of a researcher's career. In the petabyte era of astronomical science, citing software and measuring its impact enables academia to retain and reward researchers that make significant software contributions. These highly skilled researchers must be retained to maximize the scientific return from petabyte-scale datasets. Evolving beyond the one-dimensional credit model requires overcoming several key challenges, including the current scholarly ecosystem and scientific culture issues. This white paper will present these challenges and suggest practical solutions for elevating the role of software as a product of the research enterprise.

Over the past 30 years, spectropolarimetry has proven to be a great tool to probe the 3D geometry of core collapse supernovae (CCSNe). The number of high-quality multi-epoch spectropolarimetric data sets for CCSNe remains quite low, however, due to the challenging nature of such observations. In this Thesis we present and analyse the data of two Type IIb SNe (SN 2008aq - two epochs, SN 2011hs - seven epochs) and one Type Ic-bl (SN 2014ad - seven epochs). The latter is the most complete spectropolarimetric data set for a Type Ic-bl to date. We also provide a preliminary re-analysis of the spectropolarimetric data of SN 1993J. The new datasets demonstrate that the spectropolarimetry of stripped envelope supernovae is complex and varied, and although early studies of type IIb SNe showed similar behaviour to type IIP, it is now clear that this early trend is broken by a growing sample of SNe. Additionally, we identified potential signatures of a jet in SN 2014ad. Furthermore, a toy model to simulate spectropolarimetry from simple geometries (similar to Reilly et al. 2016) was created. A systematic approach to parameter space exploration was developped; this showed that unfortunately these models suffer extreme degeneracies and may not lead to accurate interpretations. Finally, spectropolarimetry can also be used to probe asymmetries caused by fast rotation in the potential progenitors of stripped envelope CCSNe. Although numerous Wolf-Rayet stars of type WN and WC have been observed with spectropolarimetry, no single WO stars have been studied to date. In our data of two Galactic WO stars, we found no line effect and could not sufficiently constrain the rotational velocity to exclude a collapsar scenario. Therefore, the absence of a line effect may not necessarily equate to a low rotational velocity.

The non-detection of companion stars in Type Ia supernova (SN) progenitor systems lends support to the notion of double-degenerate (DD) systems and explosions triggered by the merging of two white dwarfs. This very asymmetric process should lead to a conspicuous polarimetric signature. By contrast, observations consistently find very low continuum polarization as the signatures from the explosion process largely dominate over the pre-explosion configuration within several days. Critical information about the interaction of the ejecta with a companion and any circumstellar matter is encoded in the early polarization spectra. In this study, we obtain spectropolarimetry of SN\,2018gv with the ESO Very Large Telescope at $-$13.6 days relative to the $B-$band maximum light, or $\sim$5 days after the estimated explosion --- the earliest spectropolarimetric observations to date of any Type Ia SN. These early observations still show a low continuum polarization ($\lesssim$0.2\%) and moderate line polarization (0.30$\pm$0.04\% for the prominent \ion{Si}{2} $\lambda$6355 feature and 0.85$\pm$0.04\% for the high-velocity Ca component). The high degree of spherical symmetry implied by the low line and continuum polarization at this early epoch is consistent with explosion models of delayed detonations and is inconsistent with the merger-induced explosion scenario. The dense UV and optical photometry and optical spectroscopy within the first $\sim$100 days after the maximum light indicate that SN\,2018gv is a normal Type Ia SN with similar spectrophotometric behavior to SN\,2011fe.

We observed seven epochs of spectropolarimetry in optical wavelengths for the Type IIb SN 2011hs, ranging from -3 to +40 days with respect to V -band maximum. A high degree of interstellar polarization was detected (up to ~3 percent), with a peak lying blueward of 4500A. Similar behaviours have been seen in some Type Ia SNe, but had never been observed in a Type IIb. We find that it is most likely the result of a relative enhancement of small silicate grains in the vicinity of the SN. Significant intrinsic continuum polarization was recovered at -3 and +2 days (p = 0.55 +\- 0.12 percent and p = 0.75 +\- 0.11 percent, respectively). We discuss the change of the polarization angle across spectral lines and in the continuum as diagnostics for the 3D structure of the ejecta. We see a gradual rotation by about -50 degree in the continuum polarization angle between -2 and +18 days after V - band maximum. A similar rotation in He I {\lambda}5876, H{\alpha} and the Ca II infrared triplet seems to indicate a strong influence of the global geometry on the line polarization features. The differences in the evolution of their respective loops on the Stokes q - u plane suggest that line specific geometries are also being probed. Possible interpretations are discussed and placed in the context of literature. We find that the spectropolarimetry of SN 2011hs is most similar to that of SN 2011dh, albeit with notable differences.

We present SPLOT, a small scale pilot survey to test the potential of snapshot (single epoch) linear imaging polarimetry as a supplementary tool to traditional transient follow-up. Transients exist in a vast volume of observational parameter space and polarimetry has the potential to highlight sources of scientific interest and add value to near real-time transient survey streams. We observed a sample of $\sim 50$ randomly selected optical transients with the EFOSC2 and SofI instruments, on the 3.6m New Technology Telescope (NTT) to test the feasibility of the survey. Our sample contained a number of interesting individual sources: a variety of supernovae, X-ray binaries, a tidal disruption event, blazar outbursts, and, by design, numerous transients of unknown nature. We discuss the results, both for the individual sources and the survey in detail. We provide an overview on the success and limitations of SPLOT and also describe a novel calibration method for removing instrumental polarisation effects from Nasymth-mounted telescopes. We find that a SPLOT-like survey would be a benefit to the large scale future transient survey streams such as LSST. The polarimetric measurements have added scientific value to a significant number of the sources and, most importantly, has shown the potential to highlight unclassified transient sources of scientific interest for further study.

Oxygen sequence Wolf-Rayet stars (WO) are thought to be the final evolution phase of some high mass stars, as such they may be the progenitors of type Ic SNe as well as potential progenitors of broad-lined Ic and long gamma-ray bursts. We present the first spectropolarimetric observations of the Galactic WO stars WR93b and WR102 obtained with FORS1 on the VLT. We find no sign of a line effect, which could be expected if these stars were rapid rotators. We also place constraints on the amplitude of a potentially undetected line effect. This allows us to derive upper limits on the possible intrinsic continuum polarisation, and find P$_{\rm cont}$ < 0.077 percent and P$_{\rm cont}$ < 0.057 percent for WR93b and WR102, respectively. Furthermore, we derive upper limits on the rotation of our WO stars by considering our results in the context of the wind compression effect. We estimate that for an edge-on case the rotational velocity of WR93b is v$_{\rm rot}$ < 324 km/s while for WR102 v$_{\rm rot}$ < 234 km/s. These correspond to values of v$_{\rm rot}$/v$_{\rm crit}$ <19 percent and <10 percent, respectively, and values of log(j)<18.0 cm$^2$/s for WR93b and <17.6 cm^2 /s for WR102. The upper limits found on v$_{\rm rot}$/v$_{\rm crit}$ and log(j) for our WO stars are therefore similar to the estimates calculated for Galactic WR stars that do show a line effect. Therefore, although the presence of a line effect in single WR stars is indicative of fast rotation, the absence of a line effect does not rule out significant rotation, even when considering the edge-on scenario.

Optical polarimetry is an effective way of probing the environment of supernova for dust. We acquired linear HST ACS/WFC polarimetry in bands $F475W$, $F606W$, and $F775W$ of the supernova (SN) 2014J in M82 at six epochs from $\sim$277 days to $\sim$1181 days after the $B$-band maximum. The polarization measured at day 277 shows conspicuous deviations from other epochs. These differences can be attributed to at least $\sim$ 10$^{-6} M_{\odot}$ of circumstellar dust located at a distance of $\sim5\times10^{17}$ cm from the SN. The scattering dust grains revealed by these observations seem to be aligned with the dust in the interstellar medium that is responsible for the large reddening towards the supernova. The presence of this circumstellar dust sets strong constraints on the progenitor system that led to the explosion of SN\,2014J; however, it cannot discriminate between single- and double-degenerate models.

The very nearby Type Ia supernova 2014J in M82 offers a rare opportunity to study the physics of thermonuclear supernovae at extremely late phases ($\gtrsim$800 days). Using the Hubble Space Telescope (HST), we obtained six epochs of high precision photometry for SN 2014J from 277 days to 1181 days past the $B-$band maximum light. The reprocessing of electrons and X-rays emitted by the radioactive decay chain $^{57}$Co$\rightarrow ^{57}$Fe are needed to explain the significant flattening of both the $F606W$-band and the pseudo-bolometric light curves. The flattening confirms previous predictions that the late-time evolution of type Ia supernova luminosities requires additional energy input from the decay of $^{57}$Co (Seitenzahl et al. 2009). By assuming the $F606W$-band luminosity scales with the bolometric luminosity at $\sim$500 days after the $B-$band maximum light, a mass ratio $^{57}$Ni/$^{56}$Ni$\sim$0.065$_{-0.004}^{+0.005}$ is required. This mass ratio is roughly $\sim$3 times the solar ratio and favors a progenitor white dwarf with a mass near the Chandrasekhar limit. A similar fit using the constructed pseudo-bolometric luminosity gives a mass ratio $^{57}$Ni/$^{56}$Ni$\sim$0.066$_{-0.008}^{+0.009}$. Astrometric tests based on the multi-epoch HST ACS/WFC images reveal no significant circumstellar light echoes in between 0.3 pc and 100 pc (Yang et al. 2017) from the supernova.

We present optical spectropolarimetry and spectroscopy of the broad-lined Type Ic (Ic-bl) SN 2014ad. Our spectropolarimetric observations cover 7 epochs, from -2 days to 66 days after V-band maximum, and the spectroscopic data were acquired from -2 days to +107 days. The photospheric velocity estimates showed ejecta speeds similar to those of SN 1998bw and other SNe associated with GRBs. The spectropolarimetric data revealed aspherical outer ejecta and a nearly spherical interior. The polarisation associated with O I {\lambda}7774 and the Ca II infrared triplet suggests a clumpy and highly asymmetrical distribution of these two species within the ejecta. Furthermore it was shown that the two line forming regions must have been spatially distinct and oxygen was found to have higher velocities than calcium. Another oxygen line-forming region was also identified much closer to the core of the explosion and distributed in a spherical shell. It is difficult to reconcile the geometry of the deeper ejecta with a jet driven explosion, but the high ejecta velocities of SN 2014ad are typical of those observed in SNe Ic-bl with GRBs, and the behaviour of the oxygen and calcium line-forming regions is consistent with fully jet-driven models. The metallicity of the host galaxy of SN 2014ad was also calculated and compared to that of the hosts of other SNe Ic-bl with and without GRBs, but due to the overlap in the two populations no conclusion could be drawn.

We present multiple-epoch measurements of the size and surface brightness of the light echoes from supernova (SN) 2014J in the nearby starburst galaxy M82. Hubble Space Telescope (HST) ACS/WFC images were taken ~277 and ~416 days after B-band maximum in the filters F475W, F606W, and F775W. Observations with HST WFC3/UVIS images at epochs ~216 and ~365 days (Crotts 2015) are included for a more complete analysis. The images reveal the temporal evolution of at least two major light-echo components. The first one exhibits a filled ring structure with position-angle-dependent intensity. This radially extended, diffuse echo indicates the presence of an inhomogeneous interstellar dust cloud ranging from ~100 pc to ~500 pc in the foreground of the SN. The second echo component appears as an unresolved luminous quarter-circle arc centered on the SN. The wavelength dependence of scattering measured in different dust components suggests that the dust producing the luminous arc favors smaller grain sizes, while that causing the diffuse light echo may have sizes similar to those of the Milky Way dust. Smaller grains can produce an optical depth consistent with that along the supernova-Earth line of sight measured by previous studies around maximum light. Therefore, it is possible that the dust slab, from which the luminous arc arises, is also responsible for most of the extinction towards SN 2014J. The optical depths determined from the Milky Way-like dust in the scattering matters are lower than that produced by the dust slab.