Showing 1–46 of 46 results for author: Stevance, H

Discovery and Extensive Follow-Up of SN 2024ggi, a nearby type IIP supernova in NGC 3621

Authors: Ting-Wan Chen, Sheng Yang, Shubham Srivastav, Takashi J. Moriya, Stephen J. Smartt, Sofia Rest, Armin Rest, Hsing Wen Lin, Hao-Yu Miao, Yu-Chi Cheng, Amar Aryan, Chia-Yu Cheng, Morgan Fraser, Li-Ching Huang, Meng-Han Lee, Cheng-Han Lai, Yu Hsuan Liu, Aiswarya Sankar. K, Ken W. Smith, Heloise F. Stevance, Ze-Ning Wang, Joseph P. Anderson, Charlotte R. Angus, Thomas de Boer, Kenneth Chambers , et al. (23 additional authors not shown)

Abstract: We present the discovery and early observations of the nearby Type II supernova (SN) 2024ggi in NGC 3621 at 6.64 +/- 0.3 Mpc. The SN was caught 5.8 (+1.9 -2.9) hours after its explosion by the ATLAS survey. Early-phase, high-cadence, and multi-band photometric follow-up was performed by the Kinder (Kilonova Finder) project, collecting over 1000 photometric data points within a week. The combined o… ▽ More We present the discovery and early observations of the nearby Type II supernova (SN) 2024ggi in NGC 3621 at 6.64 +/- 0.3 Mpc. The SN was caught 5.8 (+1.9 -2.9) hours after its explosion by the ATLAS survey. Early-phase, high-cadence, and multi-band photometric follow-up was performed by the Kinder (Kilonova Finder) project, collecting over 1000 photometric data points within a week. The combined o- and r-band light curves show a rapid rise of 3.3 magnitudes in 13.7 hours, much faster than SN 2023ixf (another recent, nearby, and well-observed SN II). Between 13.8 and 18.8 hours after explosion SN 2024ggi became bluer, with u-g colour dropping from 0.53 to 0.15 mag. The rapid blueward evolution indicates a wind shock breakout (SBO) scenario. No hour-long brightening expected for the SBO from a bare stellar surface was detected during our observations. The classification spectrum, taken 17 hours after the SN explosion, shows flash features of high-ionization species such as Balmer lines, He I, C III, and N III. Detailed light curve modeling reveals critical insights into the properties of the circumstellar material (CSM). Our favoured model has an explosion energy of 2 x 10^51 erg, a mass-loss rate of 10^-3 solar_mass/yr (with an assumed 10 km/s wind), and a confined CSM radius of 6 x 10^14 cm. The corresponding CSM mass is 0.4 solar_mass. Comparisons with SN 2023ixf highlight that SN 2024ggi has a smaller CSM density, resulting in a faster rise and fainter UV flux. The extensive dataset and the involvement of citizen astronomers underscore that a collaborative network is essential for SBO searches, leading to more precise and comprehensive SN characterizations. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 11 pages, 5 figures in manuscript, 6 pages in appendix, submitted to ApJL

SN 2023zaw: the low-energy explosion of an ultra-stripped star, with non-radioactive heating

Authors: Thomas Moore, James Gillanders, Matt Nicholl, Mark Huber, Stephen Smartt, Shubham Srivastav, Heloise Stevance, Ting-Wan Chen, Kenneth Chambers, Joseph Anderson, Michael Fulton, Samantha Oates, Charlotte Angus, Giuliano Pignata, Nicolas Erasmus, Hua Gao, Joanna Bulger, Chien-Cheng Lin, Thomas Lowe, Eugene Magnier, Paloma Minguez, Chow-Choong Ngeow, Xinyue Sheng, Stuart A. Sim, Ken Smith , et al. (4 additional authors not shown)

Abstract: Most stripped envelope supernova progenitors are formed through binary interaction, losing hydrogen and/or helium from their outer layers. An emerging class of supernovae with the highest degree of envelope-stripping are thought to be the product of stripping by a NS companion. However, relatively few examples are known and the outcomes of such systems can be diverse and are poorly understood at p… ▽ More Most stripped envelope supernova progenitors are formed through binary interaction, losing hydrogen and/or helium from their outer layers. An emerging class of supernovae with the highest degree of envelope-stripping are thought to be the product of stripping by a NS companion. However, relatively few examples are known and the outcomes of such systems can be diverse and are poorly understood at present. Here, we present spectroscopic observations and high cadence multi-band photometry of SN 2023zaw, a low ejecta mass and rapidly evolving supernova. SN 2023zaw was discovered in a nearby spiral galaxy at D = 39.7 Mpc, with significant Milky Way extinction, $E(B-V) = 0.21$, and significant (but uncertain) host extinction. Bayesian evidence comparison reveals that nickel is not the only power source and an additional energy source is required to explain our observations. Our models suggest an ejecta mass of $M_{\rm ej} \sim 0.07\,\rm M_\odot$ and a synthesised nickel mass of $M_{\rm ej} \sim 0.007\,\rm M_\odot$ is required to explain the explosion. However an additional heating from a magnetar or interaction with circumstellar material is required to power the early light curve. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Discovery of the optical and radio counterpart to the fast X-ray transient EP240315a

Authors: J. H. Gillanders, L. Rhodes, S. Srivastav, F. Carotenuto, J. Bright, M. E. Huber, H. F. Stevance, S. J. Smartt, K. C. Chambers, T. -W. Chen, R. Fender, A. Andersson, A. J. Cooper, P. G. Jonker, F. J. Cowie, T. deBoer, N. Erasmus, M. D. Fulton, H. Gao, J. Herman, C. -C. Lin, T. Lowe, E. A. Magnier, H. -Y. Miao, P. Minguez , et al. (14 additional authors not shown)

Abstract: Fast X-ray Transients (FXTs) are extragalactic bursts of soft X-rays first identified >10 years ago. Since then, nearly 40 events have been discovered, although almost all of these have been recovered from archival Chandra and XMM-Newton data. To date, optical sky surveys and follow-up searches have not revealed any multi-wavelength counterparts. The Einstein Probe, launched in January 2024, has s… ▽ More Fast X-ray Transients (FXTs) are extragalactic bursts of soft X-rays first identified >10 years ago. Since then, nearly 40 events have been discovered, although almost all of these have been recovered from archival Chandra and XMM-Newton data. To date, optical sky surveys and follow-up searches have not revealed any multi-wavelength counterparts. The Einstein Probe, launched in January 2024, has started surveying the sky in the soft X-ray regime (0.5-4 keV) and will rapidly increase the sample of FXTs discovered in real time. Here, we report the first discovery of both an optical and radio counterpart to a distant FXT, the fourth source publicly released by the Einstein Probe. We discovered a fast-fading optical transient within the 3 arcmin localisation radius of EP240315a with the all-sky optical survey ATLAS, and our follow-up Gemini spectrum provides a redshift, z=4.859+/-0.002. Furthermore, we uncovered a radio counterpart in the S-band (3.0 GHz) with the MeerKAT radio interferometer. The optical (rest-frame UV) and radio luminosities indicate the FXT most likely originates from either a long gamma-ray burst or a relativistic tidal disruption event. This may be a fortuitous early mission detection by the Einstein Probe or may signpost a mode of discovery for high-redshift, high-energy transients through soft X-ray surveys, combined with locating multi-wavelength counterparts. △ Less

Submitted 19 June, 2024; v1 submitted 16 April, 2024; originally announced April 2024.

Comments: Updated to match version accepted for publication in ApJL (17 pages, 4 figures, 2 tables)

The enigmatic double-peaked stripped-envelope SN 2023aew

Authors: Tuomas Kangas, Hanindyo Kuncarayakti, Takashi Nagao, Rubina Kotak, Erkki Kankare, Morgan Fraser, Heloise Stevance, Seppo Mattila, Kei'ichi Maeda, Maximilian Stritzinger, Peter Lundqvist, Nancy Elias-Rosa, Lucía Ferrari, Gastón Folatelli, Christopher Frohmaier, Lluís Galbany, Miho Kawabata, Eleni Koutsiona, Tomás E. Müller-Bravo, Lara Piscarreta, Miika Pursiainen, Avinash Singh, Kenta Taguchi, Rishabh Singh Teja, Giorgio Valerin , et al. (7 additional authors not shown)

Abstract: We present optical and near-infrared photometry and spectroscopy of SN 2023aew and our findings on its remarkable properties. This event, initially resembling a Type IIb supernova (SN), rebrightens dramatically $\sim$90 d after the first peak, at which time its spectrum transforms into that of a SN Ic. The slowly evolving spectrum specifically resembles a post-peak SN~Ic with relatively low line v… ▽ More We present optical and near-infrared photometry and spectroscopy of SN 2023aew and our findings on its remarkable properties. This event, initially resembling a Type IIb supernova (SN), rebrightens dramatically $\sim$90 d after the first peak, at which time its spectrum transforms into that of a SN Ic. The slowly evolving spectrum specifically resembles a post-peak SN~Ic with relatively low line velocities even during the second rise. The second peak, reached 119 d after the first peak, is both more luminous ($M_r = -18.75\pm0.04$ mag) and much broader than those of typical SNe Ic. Blackbody fits to SN 2023aew indicate that the photosphere shrinks almost throughout its observed evolution, and the second peak is caused by an increasing temperature. Bumps in the light curve after the second peak suggest interaction with circumstellar matter (CSM) or possibly accretion. We consider several scenarios for producing the unprecedented behavior of SN 2023aew. Two separate SNe, either unrelated or from the same binary system, require either an incredible coincidence or extreme fine-tuning. A pre-SN eruption followed by a SN requires an extremely powerful, SN-like eruption (consistent with $\sim$10$^{51}$ erg) and is also disfavored. We therefore consider only the first peak a true stellar explosion. The observed evolution is difficult to reproduce if the second peak is dominated by interaction with a distant CSM shell. A delayed internal heating mechanism is more likely, but emerging embedded interaction with a CSM disk should be accompanied by CSM lines in the spectrum, which are not observed, and is difficult to hide long enough. A magnetar central engine requires a delayed onset to explain the long time between the peaks. Delayed fallback accretion onto a black hole may present the most promising scenario, but we cannot definitively establish the power source. △ Less

Submitted 17 June, 2024; v1 submitted 30 January, 2024; originally announced January 2024.

Comments: 21 pages + 14 appendix pages, 16 figures + 6 appendix figures. Accepted by A&A

NEural Engine for Discovering Luminous Events (NEEDLE): identifying rare transient candidates in real time from host galaxy images

Authors: Xinyue Sheng, Matt Nicholl, Ken W. Smith, David R. Young, Roy D. Williams, Heloise F. Stevance, Stephen J. Smartt, Shubham Srivastav, Thomas Moore

Abstract: Known for their efficiency in analyzing large data sets, machine learning classifiers are widely used in wide-field sky surveys. The upcoming Vera C. Rubin Observatory Legacy of Time and Space Survey (LSST) will generate millions of alerts every night, enabling the discovery of large samples of rare events. Identifying such objects soon after explosion will be essential to study their evolution. T… ▽ More Known for their efficiency in analyzing large data sets, machine learning classifiers are widely used in wide-field sky surveys. The upcoming Vera C. Rubin Observatory Legacy of Time and Space Survey (LSST) will generate millions of alerts every night, enabling the discovery of large samples of rare events. Identifying such objects soon after explosion will be essential to study their evolution. This requires a machine learning framework that makes use of all available transient and contextual information. Using $\sim5400$ transients from the ZTF Bright Transient Survey as input data, we develop NEEDLE, a novel hybrid classifier to select for two rare classes with strong environmental preferences: superluminous supernovae (SLSNe) preferring dwarf galaxies, and tidal disruption events (TDEs) occurring in the centres of nucleated galaxies. The input data includes detection and reference images, photometric information from the alert packets, and host galaxy magnitudes from Pan-STARRS. Despite having only a few tens of examples of the rare classes, our average (best) completeness on an unseen test set reaches 77% (93%) for SLSNe and 72% (87%) for TDEs. This may still result in a large fraction of false positives for the rare transients, given the large class imbalance in real surveys. However, the goal of NEEDLE is to find good candidates for spectroscopic classification, rather than to select pure photometric samples. Our network is designed with LSST in mind and we expect performance to improve further with the higher resolution images and more accurate transient and host photometry that will be available from Rubin. Our system will be deployed as an annotator on the UK alert broker, Lasair, to provide predictions to the community in real time. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Comments: Submitted to MNRAS

GW190425: Pan-STARRS and ATLAS coverage of the skymap and limits on optical emission associated with FRB190425

Authors: S. J. Smartt, M. Nicholl, S. Srivastav, M. E. Huber, K. C. Chambers, K. W. Smith, D. R. Young, M. D. Fulton, J. L. Tonry, C. W. Stubbs, L. Denneau, A. J. Cooper, A. Aamer, J. P. Anderson, A. Andersson, J. Bulger, T. -W Chen, P. Clark, T. de Boer, H. Gao, J. H. Gillanders, A. Lawrence, C. C. Lin, T. B. Lowe, E. A. Magnier , et al. (10 additional authors not shown)

Abstract: GW190425 is the second of only two binary neutron star (BNS) merger events to be significantly detected by the LIGO-Virgo- Kagra gravitational wave detectors. With a detection only in LIGO Livingston, the skymap containing the source was large and no plausible electromagnetic counterpart was found in real time searching in 2019. Here we summarise our ATLAS and Pan-STARRS wide-field optical coverag… ▽ More GW190425 is the second of only two binary neutron star (BNS) merger events to be significantly detected by the LIGO-Virgo- Kagra gravitational wave detectors. With a detection only in LIGO Livingston, the skymap containing the source was large and no plausible electromagnetic counterpart was found in real time searching in 2019. Here we summarise our ATLAS and Pan-STARRS wide-field optical coverage of the skymap beginning within 1 hour and 3 hours respectively of the GW190425 merger time. More recently, a potential coincidence between GW190425 and a fast radio burst FRB 190425 has been suggested, given their spatial and temporal coincidence. The smaller sky localisation area of FRB 190425 and its dispersion measure have led to the identification of a likely host galaxy, UGC 10667 at a distance of 141 +/- 10 Mpc. Our optical imaging covered the galaxy 6.0 hrs after GW190425 was detected and 3.5 hrs after the FRB 190425. No optical emission was detected and further imaging at +1.2 and +13.2 days also revealed no emission. If the FRB 190425 and GW190425 association were real, we highlight our limits on kilonova emission from a BNS merger in UGC 10667. The model for producing FRB 190425 from a BNS merger involves a supramassive magnetised neutron star spinning down by dipole emission on the timescale of hours. We show that magnetar enhanced kilonova emission is ruled out by optical upper limits. The lack of detected optical emission from a kilonova in UGC 10667 disfavours, but does not disprove, the FRB-GW link for this source. △ Less

Submitted 20 September, 2023; originally announced September 2023.

Comments: Submitted to MNRAS, 20th Sept 2023, 9 pages

JWST detection of heavy neutron capture elements in a compact object merger

Authors: A. Levan, B. P. Gompertz, O. S. Salafia, M. Bulla, E. Burns, K. Hotokezaka, L. Izzo, G. P. Lamb, D. B. Malesani, S. R. Oates, M. E. Ravasio, A. Rouco Escorial, B. Schneider, N. Sarin, S. Schulze, N. R. Tanvir, K. Ackley, G. Anderson, G. B. Brammer, L. Christensen, V. S. Dhillon, P. A. Evans, M. Fausnaugh, W. -F. Fong, A. S. Fruchter , et al. (58 additional authors not shown)

Abstract: The mergers of binary compact objects such as neutron stars and black holes are of central interest to several areas of astrophysics, including as the progenitors of gamma-ray bursts (GRBs), sources of high-frequency gravitational waves and likely production sites for heavy element nucleosynthesis via rapid neutron capture (the r-process). These heavy elements include some of great geophysical, bi… ▽ More The mergers of binary compact objects such as neutron stars and black holes are of central interest to several areas of astrophysics, including as the progenitors of gamma-ray bursts (GRBs), sources of high-frequency gravitational waves and likely production sites for heavy element nucleosynthesis via rapid neutron capture (the r-process). These heavy elements include some of great geophysical, biological and cultural importance, such as thorium, iodine and gold. Here we present observations of the exceptionally bright gamma-ray burst GRB 230307A. We show that GRB 230307A belongs to the class of long-duration gamma-ray bursts associated with compact object mergers, and contains a kilonova similar to AT2017gfo, associated with the gravitational-wave merger GW170817. We obtained James Webb Space Telescope mid-infrared (mid-IR) imaging and spectroscopy 29 and 61 days after the burst. The spectroscopy shows an emission line at 2.15 microns which we interpret as tellurium (atomic mass A=130), and a very red source, emitting most of its light in the mid-IR due to the production of lanthanides. These observations demonstrate that nucleosynthesis in GRBs can create r-process elements across a broad atomic mass range and play a central role in heavy element nucleosynthesis across the Universe. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: Submitted. Comments welcome! Nature (2023)

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

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

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

Submitted 12 January, 2023; originally announced January 2023.

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

Panning for gold, but finding helium: discovery of the ultra-stripped supernova SN2019wxt from gravitational-wave follow-up observations

Authors: I. Agudo, L. Amati, T. An, F. E. Bauer, S. Benetti, M. G. Bernardini, R. Beswick, K. Bhirombhakdi, T. de Boer, M. Branchesi, S. J. Brennan, M. D. Caballero-García, E. Cappellaro, N. Castro Rodríguez, A. J. Castro-Tirado, K. C. Chambers, E. Chassande-Mottin, S. Chaty, T. -W. Chen, A. Coleiro, S. Covino, F. D'Ammando, P. D'Avanzo, V. D'Elia, A. Fiore , et al. (74 additional authors not shown)

Abstract: We present the results from multi-wavelength observations of a transient discovered during the follow-up of S191213g, a gravitational wave (GW) event reported by the LIGO-Virgo Collaboration as a possible binary neutron star merger in a low latency search. This search yielded SN2019wxt, a young transient in a galaxy whose sky position (in the 80\% GW contour) and distance ($\sim$150\,Mpc) were pla… ▽ More We present the results from multi-wavelength observations of a transient discovered during the follow-up of S191213g, a gravitational wave (GW) event reported by the LIGO-Virgo Collaboration as a possible binary neutron star merger in a low latency search. This search yielded SN2019wxt, a young transient in a galaxy whose sky position (in the 80\% GW contour) and distance ($\sim$150\,Mpc) were plausibly compatible with the localisation uncertainty of the GW event. Initially, the transient's tightly constrained age, its relatively faint peak magnitude ($M_i \sim -16.7$\,mag) and the $r-$band decline rate of $\sim 1$\,mag per 5\,days appeared suggestive of a compact binary merger. However, SN2019wxt spectroscopically resembled a type Ib supernova, and analysis of the optical-near-infrared evolution rapidly led to the conclusion that while it could not be associated with S191213g, it nevertheless represented an extreme outcome of stellar evolution. By modelling the light curve, we estimated an ejecta mass of $\sim 0.1\,M_\odot$, with $^{56}$Ni comprising $\sim 20\%$ of this. We were broadly able to reproduce its spectral evolution with a composition dominated by helium and oxygen, with trace amounts of calcium. We considered various progenitors that could give rise to the observed properties of SN2019wxt, and concluded that an ultra-stripped origin in a binary system is the most likely explanation. Disentangling electromagnetic counterparts to GW events from transients such as SN2019wxt is challenging: in a bid to characterise the level of contamination, we estimated the rate of events with properties comparable to those of SN2019wxt and found that $\sim 1$ such event per week can occur within the typical GW localisation area of O4 alerts out to a luminosity distance of 500\,Mpc, beyond which it would become fainter than the typical depth of current electromagnetic follow-up campaigns. △ Less

Submitted 20 June, 2023; v1 submitted 18 August, 2022; originally announced August 2022.

Comments: By the ENGRAVE collaboration (engrave-eso.org). 35 pages, 20 figures, final version accepted by A&A

Journal ref: A&A 675, A201 (2023)

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

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

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

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

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

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

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

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

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

Comments: Accepted by MNRAS

VFTS 243 as predicted by the BPASS fiducial models

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

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

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

Comments: 4 figures, 1 table, Published to MNRAS

UV Spectropolarimetry with Polstar: Massive Star Binary Colliding Winds

Authors: Nicole St-Louis, Kenneth Gayley, D. John Hillier, Richard Ignace, Carol E. Jones, Alexandre David-Uraz, Noel D. Richardson, Jorick S. Vink, Geraldine J. Peters, Jennifer L. Hoffman, Yael, Nazé, Heloise Stevance, Tomer Shenar, Andrew G. Fullard, Jaimie R. Lomax, Paul A. Scowen

Abstract: The winds of massive stars are important for their direct impact on the interstellar medium, and for their influence on the final state of a star prior to it exploding as a supernova. However, the dynamics of these winds is understood primarily via their illumination from a single central source. The Doppler shift seen in resonance lines is a useful tool for inferring these dynamics, but the mappi… ▽ More The winds of massive stars are important for their direct impact on the interstellar medium, and for their influence on the final state of a star prior to it exploding as a supernova. However, the dynamics of these winds is understood primarily via their illumination from a single central source. The Doppler shift seen in resonance lines is a useful tool for inferring these dynamics, but the mapping from that Doppler shift to the radial distance from the source is ambiguous. Binary systems can reduce this ambiguity by providing a second light source at a known radius in the wind, seen from orbitally modulated directions. From the nature of the collision between the winds, a massive companion also provides unique additional information about wind momentum fluxes. Since massive stars are strong ultraviolet (UV) sources, and UV resonance line opacity in the wind is strong, UV instruments with a high resolution spectroscopic capability are essential for extracting this dynamical information. Polarimetric capability also helps to further resolve ambiguities in aspects of the wind geometry that are not axisymmetric about the line of sight, because of its unique access to scattering direction information. We review how the proposed MIDEX-scale mission Polstar can use UV spectropolarimetric observations to critically constrain the physics of colliding winds, and hence radiatively-driven winds in general. We propose a sample of 20 binary targets, capitalizing on this unique combination of illumination by companion starlight, and collision with a companion wind, to probe wind attributes over a range in wind strengths. Of particular interest is the hypothesis that the radial distribution of the wind acceleration is altered significantly, when the radiative transfer within the winds becomes optically thick to resonance scattering in multiple overlapping UV lines. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 26 pages, 12 figures, Review in a topical collection series of Astrophysics and Space Sciences on the proposed Polstar satellite. arXiv admin note: substantial text overlap with arXiv:2111.11552

Observations of the luminous red nova AT 2021biy in the nearby galaxy NGC 4631

Authors: Y. -Z. Cai, A. Pastorello, M. Fraser, X. -F. Wang, A. V. Filippenko, A. Reguitti, K. C. Patra, V. P. Goranskij, E. A. Barsukova, T. G. Brink, N. Elias-Rosa, H. F. Stevance, W. Zheng, Y. Yang, K. E. Atapin, S. Benetti, T. J. L. de Boer, S. Bose, J. Burke, R. Byrne, E. Cappellaro, K. C. Chambers, W. -L. Chen, N. Emami, H. Gao , et al. (51 additional authors not shown)

Abstract: We present an observational study of the luminous red nova (LRN) AT\,2021biy in the nearby galaxy NGC\,4631. The field of the object was routinely imaged during the pre-eruptive stage by synoptic surveys, but the transient was detected only at a few epochs from $\sim 231$\,days before maximum brightness. The LRN outburst was monitored with unprecedented cadence both photometrically and spectroscop… ▽ More We present an observational study of the luminous red nova (LRN) AT\,2021biy in the nearby galaxy NGC\,4631. The field of the object was routinely imaged during the pre-eruptive stage by synoptic surveys, but the transient was detected only at a few epochs from $\sim 231$\,days before maximum brightness. The LRN outburst was monitored with unprecedented cadence both photometrically and spectroscopically. AT\,2021biy shows a short-duration blue peak, with a bolometric luminosity of $\sim 1.6 \times 10^{41}$\,erg\,s$^{-1}$, followed by the longest plateau among LRNe to date, with a duration of 210\,days. A late-time hump in the light curve was also observed, possibly produced by a shell-shell collision. AT\,2021biy exhibits the typical spectral evolution of LRNe. Early-time spectra are characterised by a blue continuum and prominent H emission lines. Then, the continuum becomes redder, resembling that of a K-type star with a forest of metal absorption lines during the plateau phase. Finally, late-time spectra show a very red continuum ($T_{\mathrm{BB}} \approx 2050$ K) with molecular features (e.g., TiO) resembling those of M-type stars. Spectropolarimetric analysis indicates that AT\,2021biy has local dust properties similar to those of V838\,Mon in the Milky Way Galaxy. Inspection of archival {\it Hubble Space Telescope} data taken on 2003 August 3 reveals a $\sim 20$\,\msun\ progenitor candidate with log\,$(L/{\rm L}_{\odot}) = 5.0$\,dex and $T_{\rm{eff}} = 5900$\,K at solar metallicity. The above luminosity and colour match those of a luminous yellow supergiant. Most likely, this source is a close binary, with a 17--24\,\msun\ primary component. △ Less

Submitted 27 August, 2022; v1 submitted 2 July, 2022; originally announced July 2022.

Comments: 21 pages, 14 figures. Accepted by Astronomy and Astrophysics

Journal ref: A&A 667, A4 (2022)

What can Gaussian Processes really tell us about supernova lightcurves? Consequences for Type II(b) morphologies and genealogies

Authors: H. F. Stevance, A. Lee

Abstract: Machine learning has become widely used in astronomy. Gaussian Process (GP) regression in particular has been employed a number of times to fit or re-sample supernova (SN) light-curves, however by their nature typical GP models are not suited to fit SN photometric data and they will be prone to over-fitting. Recently GP re-sampling was used in the context of studying the morphologies of type II an… ▽ More Machine learning has become widely used in astronomy. Gaussian Process (GP) regression in particular has been employed a number of times to fit or re-sample supernova (SN) light-curves, however by their nature typical GP models are not suited to fit SN photometric data and they will be prone to over-fitting. Recently GP re-sampling was used in the context of studying the morphologies of type II and IIb SNe and they were found to be clearly distinct with respect to four parameters: the rise time (t$_{\rm rise}$), the magnitude difference between 40 and 30 days post explosion ($Δm_{\rm 40-30}$), the earliest maximum (post-peak) of the first derivative (dm1) and minimum of the second derivative (dm2). Here we take a close look at GP regression and its limitations in the context of SN light-curves in general, and we also discuss the uncertainties on these specific parameters, finding that dm1 and dm2 cannot give reliable astrophysical information. We do reproduce the clustering in t$_{\rm rise}$--$Δm_{\rm 40-30}$ space although it is not as clear cut as previously presented. The best strategy to accurately populate the t$_{\rm rise}$-- $Δm_{\rm 40-30}$ space will be to use an expanded sample of high quality light-curves (such as those in the ATLAS transient survey) and analytical fitting methods. Finally, using the BPASS fiducial models, we predict that future photometric studies will reveal clear clustering of the type IIb and II light curve morphologies with a distinct continuum of transitional events. △ Less

Submitted 29 November, 2022; v1 submitted 29 June, 2022; originally announced June 2022.

Comments: 13 pages, 11 figures, 2 tables, Accepted for publication in MNRAS

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

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

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

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

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

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

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

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

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

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

Ultraviolet Spectropolarimetry:Conservative and Nonconservative Mass Transfer in OB Interacting Binaries

Authors: Geraldine J. Peters, Ken Gayley, Richard Ignace, Carol E. Jones, Yael Naze, Nicole St-Louis, Heloise Stevance, Jorick S. Vink, Noel D. Richardson, Jennifer L. Hoffman, Jamie R. Lomax, Tomer Shenar, Andrew G. Fullard, Paul A. Scowen

Abstract: One objective of the Polstar spectropolarimetry mission is to characterize the degree of nonconservative mass transfer that occurs at various stages of binary evolution, from the initial mass reversal to the late Algol phase. The proposed instrument combines spectroscopic and polarimetric capabilities, where the spectroscopy can resolve Doppler shifts in UV resonance lines with 10 km/s precision,… ▽ More One objective of the Polstar spectropolarimetry mission is to characterize the degree of nonconservative mass transfer that occurs at various stages of binary evolution, from the initial mass reversal to the late Algol phase. The proposed instrument combines spectroscopic and polarimetric capabilities, where the spectroscopy can resolve Doppler shifts in UV resonance lines with 10 km/s precision, and polarimetry can resolve linear polarization with 1e-3 precision or better. The spectroscopy will identify absorption by mass streams seen in projection against the stellar disk as a function of orbital phase, hot accretion spots, as well as scattering from extended splash structures, circumbinary disks, and other flows in and above/below the orbital plane (e.g. jets) that fail to be transferred conservatively. The polarimetry affects more the light coming from material not seen against the stellar disk, allowing the geometry of the scattering to be tracked, resolving ambiguities left by the spectroscopy and light-curve information. For example, nonconservative mass streams ejected in the polar direction will produce polarization of the opposite sign from conservative transfer accreting in the orbital plane. Also, time domain coverage over a range of phases of the binary orbit are well supported by the Polstar observing strategy. Combining these elements will significantly improve our understanding of the mass transfer process and the amount of mass that can escape from the system, an important channel for changing the final mass, and ultimate supernova, of the large number of massive stars found in binaries at close enough separation to undergo interaction. △ Less

Submitted 12 October, 2022; v1 submitted 27 November, 2021; originally announced November 2021.

Comments: Refereed paper in Ap&SS 367:Topical Collection, 23 pages, 10 figures

Ultraviolet Spectropolarimetry with Polstar: Massive Star Binary Colliding Winds

Authors: Nicole St-Louis, Kenneth Gayley, Desmond John Hillier, Richard Ignace, Carol E. Jones, Alexandre David-Uraz, Noel D. Richardson, Jorick S. Vink, Geraldine J. Peters, Jennifer L. Hoffman, Yael Naze, Heloise Stevance, Tomer Shenar, Andrew G. Fullard, Jamie R. Lomax, Paul A. Scowen

Abstract: As sources of chemical enrichment, ionizing radiation and energetic feedback, massive stars drive the ecology of their host galaxies despite their relative rarity, additionally to yielding compact remnants, which can generate gravitational waves. The evolution of massive stars is crucially informed by their detailed mass-loss history; however, wind structures on a variety of scales cause important… ▽ More As sources of chemical enrichment, ionizing radiation and energetic feedback, massive stars drive the ecology of their host galaxies despite their relative rarity, additionally to yielding compact remnants, which can generate gravitational waves. The evolution of massive stars is crucially informed by their detailed mass-loss history; however, wind structures on a variety of scales cause important uncertainties on their mass-loss rates. Binary systems can place further constraints on the mass-loss properties of massive stars, especially colliding-wind binaries. In this paper, we review how the proposed MIDEX-scale mission Polstar can critically constrain the physics of colliding winds (and hence radiatively-driven winds in general) with ultraviolet spectropolarimetric observations, providing an unprecedented improvement on the accuracy of the determination of both mass-loss rates and the velocity structure of the winds of massive stars. We propose a sample of 17 targets that will allow us to study a variety of wind-colliding systems spanning a large parameter space using the spatial information yielded by both spectroscopic and polarimetric data obtained with Polstar. △ Less

Submitted 16 December, 2021; v1 submitted 22 November, 2021; originally announced November 2021.

Comments: 21 pages, 12 figures 1 table, Whitepaper

Estimating Transient Rates from Cosmological Simulations and BPASS

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

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

Submitted 15 November, 2021; originally announced November 2021.

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

Unveiling the Nature of SN 2011fh: a Young and Massive Star Gives Rise to a Luminous SN 2009ip-like Event

Authors: Thallis Pessi, Jose L. Prieto, Berto Monard, Christopher S. Kochanek, Greg Bock, Andrew J. Drake, Ori D. Fox, Stuart Parker, Heloise F. Stevance

Abstract: In recent years, many Type IIn supernovae have been found to share striking similarities with the peculiar SN 2009ip, whose true nature is still under debate. Here, we present 10 years of observations of SN 2011fh, an interacting transient with spectroscopic and photometric similarities to SN 2009ip. SN 2011fh had a M$_r \sim -16$ mag brightening event, followed by a brighter M$_r \sim -18$ mag lu… ▽ More In recent years, many Type IIn supernovae have been found to share striking similarities with the peculiar SN 2009ip, whose true nature is still under debate. Here, we present 10 years of observations of SN 2011fh, an interacting transient with spectroscopic and photometric similarities to SN 2009ip. SN 2011fh had a M$_r \sim -16$ mag brightening event, followed by a brighter M$_r \sim -18$ mag luminous outburst in August 2011. The spectra of SN 2011fh are dominated by narrow to intermediate Balmer emission lines throughout its evolution, with P Cygni profiles indicating fast-moving material at $\sim 6400 \ \textrm{km s}^{-1}$. HST/WFC3 observations from October of 2016 revealed a bright source with M$_{F814W} \approx -13.3$ mag, indicating that we are seeing the ongoing interaction of the ejecta with the circumstellar material or that the star might be going through an eruptive phase five years after the luminous outburst of 2011. Using HST photometry of the stellar cluster around SN 2011fh, we estimated an age of $\sim 4.5$ Myr for the progenitor, which implies a stellar mass of $\sim 60$ M$_\odot$, using single-star evolution models, or a mass range of $35 - 80 \ \textrm{M}_\odot$, considering a binary system. We also show that the progenitor of SN 2011fh exceeded the classical Eddington limit by a large factor in the months preceding the luminous outburst of 2011, suggesting strong super-Eddington winds as a possible mechanism for the observed mass-loss. These findings favor an energetic outburst in a young and massive star, possibly a luminous blue variable. △ Less

Submitted 24 January, 2022; v1 submitted 18 October, 2021; originally announced October 2021.

The impact of pre-supernova feedback and its dependence on environment

Authors: Anna F. Mcleod, Ahmad A. Ali, Mélanie Chevance, Lorenza Della Bruna, Andreas Schruba, Heloise F. Stevance, Angela Adamo, J. M. Diederik Kruijssen, Steven N. Longmore, Daniel R. Weisz, Peter Zeidler

Abstract: Integral field units enable resolved studies of a large number of star-forming regions across entire nearby galaxies, providing insight on the conversion of gas into stars and the feedback from the emerging stellar populations over unprecedented dynamic ranges in terms of spatial scale, star-forming region properties, and environments. We use the VLT/MUSE legacy data set covering the central $35$… ▽ More Integral field units enable resolved studies of a large number of star-forming regions across entire nearby galaxies, providing insight on the conversion of gas into stars and the feedback from the emerging stellar populations over unprecedented dynamic ranges in terms of spatial scale, star-forming region properties, and environments. We use the VLT/MUSE legacy data set covering the central $35$ arcmin$^{2}$ (${\sim}12$ kpc$^{2}$) of the nearby galaxy NGC 300 to quantify the effect of stellar feedback as a function of the local galactic environment. We extract spectra from emission line regions identified within dendrograms, combine emission line ratios and line widths to distinguish between HII regions, planetary nebulae, and supernova remnants, and compute their ionised gas properties, gas-phase oxygen abundances, and feedback-related pressure terms. For the HII regions, we find that the direct radiation pressure ($P_\mathrm{dir}$) and the pressure of the ionised gas ($P_{HII}$) weakly increase towards larger galactocentric radii, i.e. along the galaxy's (negative) abundance and (positive) extinction gradients. While the increase of $P_{HII}$ with galactocentric radius is likely due to higher photon fluxes from lower-metallicity stellar populations, we find that the increase of $P_\mathrm{dir}$ is likely driven by the combination of higher photon fluxes and enhanced dust content at larger galactocentric radii. In light of the above, we investigate the effect of increased pre-supernova feedback at larger galactocentric distances (lower metallicities and increased dust mass surface density) on the ISM, finding that supernovae at lower metallicities expand into lower-density environments, thereby enhancing the impact of supernova feedback. △ Less

Submitted 17 September, 2021; originally announced September 2021.

Comments: Accepted for publication in MNRAS

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

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

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

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

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

Binary pathways to SLSNe-I: SN 2017gci

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

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

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

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

Using Artificial Intelligence to Shed Light on the Star of Biscuits: The Jaffa Cake

Authors: H. F. Stevance

Abstract: Before Brexit, one of the greatest causes of arguments amongst British families was the question of the nature of Jaffa Cakes. Some argue that their size and host environment (the biscuit aisle) should make them a biscuit in their own right. Others consider that their physical properties (e.g. they harden rather than soften on becoming stale) suggest that they are in fact cake. In order to finally… ▽ More Before Brexit, one of the greatest causes of arguments amongst British families was the question of the nature of Jaffa Cakes. Some argue that their size and host environment (the biscuit aisle) should make them a biscuit in their own right. Others consider that their physical properties (e.g. they harden rather than soften on becoming stale) suggest that they are in fact cake. In order to finally put this debate to rest, we re-purpose technologies used to classify transient events. We train two classifiers (a Random Forest and a Support Vector Machine) on 100 recipes of traditional cakes and biscuits. Our classifiers have 95 percent and 91 percent accuracy respectively. Finally we feed two Jaffa Cake recipes to the algorithms and find that Jaffa Cakes are, without a doubt, cakes. Finally, we suggest a new theory as to why some believe Jaffa Cakes are biscuits. △ Less

Submitted 30 March, 2021; originally announced March 2021.

Comments: April fools astro-ph submission - The topic is a joke but the numbers and the analysis are real

Riroriro: Simulating gravitational waves and evaluating their detectability in Python

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

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

Submitted 11 March, 2021; originally announced March 2021.

Comments: 3 pages

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

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

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

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

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

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

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

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

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

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

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

Evaluating the impact of binary parameter uncertainty on stellar population properties

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

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

Submitted 24 April, 2020; originally announced April 2020.

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

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

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

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

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

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

The shape of SN 1993J re-analyzed

Authors: H. F. Stevance, D. Baade, J. R. Bruten, A. Cikota, A. Clocchiatti, D. C. Hines, P. Höflich, J. R. Maund, F. Patat, P. J. Vallely, J. C. Wheeler

Abstract: 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 em… ▽ More 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β}$, He\,{\sc i} $λ5876$, $\mathrm{Hα}$, He\,{\sc i} $λ6678$, He\,{\sc i} $λ7065$, and high velocity (HV) components of He\,{\sc i} $λ5876$ and $\mathrm{Hα}$. 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. △ Less

Submitted 12 March, 2020; originally announced March 2020.

Comments: 18 pages, 12 figures, Accepted for publication in MNRAS

Hoki: Making BPASS accessible through Python

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

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

Submitted 29 January, 2020; originally announced January 2020.

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

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

Authors: J. J. Eldridge, E. R. Stanway, K. Breivik, A. R. Casey, D. T. H. Steeghs, H. F. Stevance

Abstract: 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… ▽ More 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. △ Less

Submitted 24 May, 2020; v1 submitted 7 December, 2019; originally announced December 2019.

Comments: Accepted for MNRAS, 10 pages, 5 figures, 3 tables

The Growing Importance of a Tech Savvy Astronomy and Astrophysics Workforce

Authors: Dara Norman, Kelle Cruz, Vandana Desai, Britt Lundgren, Eric Bellm, Frossie Economou, Arfon Smith, Amanda Bauer, Brian Nord, Chad Schafer, Gautham Narayan, Ting Li, Erik Tollerud, Brigitta Sipocz, Heloise Stevance, Timothy Pickering, Manodeep Sinha, Joseph Harrington, Jeyhan Kartaltepe, Dany Vohl, Adrian Price-Whelan, Brian Cherinka, Chi-kwan Chan, Benjamin Weiner, Maryam Modjaz , et al. (4 additional authors not shown)

Abstract: 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 th… ▽ More 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. △ Less

Submitted 17 October, 2019; originally announced October 2019.

Comments: Submitted as a ASTRO2020 Decadal Survey APC position paper. arXiv admin note: substantial text overlap with arXiv:1905.05116

Linear spectropolarimetry of 35 Type Ia Supernovae with VLT/FORS: An analysis of the Si II line polarization

Authors: Aleksandar Cikota, Ferdinando Patat, Lifan Wang, J. Craig Wheeler, Mattia Bulla, Dietrich Baade, Peter Höflich, Stefan Cikota, Alejandro Clocchiatti, Justyn R. Maund, Heloise F. Stevance, Yi Yang

Abstract: 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 $λ$6355… ▽ More 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 $λ$6355 $Å$ 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 $Δ$m$_{15}$-p$_{\rm Si II}$ relationship identified in a previous study, and show that subluminous and transitional objects display polarization values below the $Δ$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 $λ$6355 $Å$ line before maximum brightness and the Si II line velocity and suggest that this, along with the $Δ$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 $Δ$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. △ Less

Submitted 20 August, 2019; originally announced August 2019.

Comments: 22 pages, accepted for publication in MNRAS

Astro2020 APC White Paper: Elevating the Role of Software as a Product of the Research Enterprise

Authors: Arfon M. Smith, Dara Norman, Kelle Cruz, Vandana Desai, Eric Bellm, Britt Lundgren, Frossie Economou, Brian D. Nord, Chad Schafer, Gautham Narayan, Joseph Harrington, Erik Tollerud, Brigitta Sipőcz, Timothy Pickering, Molly S. Peeples, Bruce Berriman, Peter Teuben, David Rodriguez, Andre Gradvohl, Lior Shamir, Alice Allen, Joel R. Brownstein, Adam Ginsburg, Manodeep Sinha, Cameron Hummels , et al. (20 additional authors not shown)

Abstract: 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 o… ▽ More 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. △ Less

Submitted 14 July, 2019; originally announced July 2019.

Comments: arXiv admin note: substantial text overlap with arXiv:1905.05116

Spectropolarimetry of stripped envelope core collapse supernovae and their progenitors

Authors: H. F. Stevance

Abstract: 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 - se… ▽ More 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. △ Less

Submitted 17 June, 2019; originally announced June 2019.

Comments: 208 pages - Accepted Thesis

The Young and Nearby Normal Type Ia Supernova 2018gv: UV-Optical Observations and the Earliest Spectropolarimetry

Authors: Yi Yang, Peter A. Hoeflich, Dietrich Baade, Justyn R. Maund, Lifan Wang, Peter. J. Brown, Heloise F. Stevance, Iair Arcavi, Jamie Burke, Aleksandar Cikota, Alejandro Clocchiatti, Avishay Gal-Yam, Melissa. L. Graham, Daichi Hiramatsu, Griffin Hosseinzadeh, D. Andrew Howell, Saurabh W. Jha, Curtis McCully, Ferdinando Patat, David. J. Sand, Steve Schulze, Jason Spyromilio, Stefano Valenti, Jozsef Vinko, Xiaofeng Wang , et al. (3 additional authors not shown)

Abstract: 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… ▽ More 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} $λ$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. △ Less

Submitted 26 March, 2019; originally announced March 2019.

Comments: 55 pages, 22 figures, 6 tables, submitted to AAS journal

The 3D shape of Type IIb SN 2011hs

Authors: H. F. Stevance, J. R. Maund, D. Baade, J. Bruten, A. Cikota, P. Höflich, L. Wang, J. C. Wheeler, A. Clocchiatti, J. Spyromilio, F. Patat, Y. Yang, P. Crowther

Abstract: 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… ▽ More 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 λ5876, Hα 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. △ Less

Submitted 22 January, 2019; originally announced January 2019.

Comments: Accepted for publication in MNRAS, 16 pages, 11 figures

SPLOT: a Snapshot survey for Polarised Light in Optical Transients

Authors: A. B. Higgins, K. Wiersema, S. Covino, R. L. C. Starling, H. F. Stevance, Ł. Wyrzykowski, S. T. Hodgkin, J. R. Maund, P. T. O'Brien, N. R. Tanvir

Abstract: 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 sa… ▽ More 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. △ Less

Submitted 28 November, 2018; v1 submitted 5 November, 2018; originally announced November 2018.

Comments: Accepted in MNRAS. 20 pages, 10 figures, 3 tables. Also contains a 10 page, two part appendix with unpolarised standard star measurements used for calibration (A) and individual source information and light curves (B). Corrected more typos

Probing the rotational velocity of Galactic WO stars with spectropolarimetry

Authors: H. F. Stevance, R. Ignace, P. A. Crowther, J. R. Maund, B. Davies, G. Rate

Abstract: 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, whi… ▽ More 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. △ Less

Submitted 10 July, 2018; v1 submitted 5 July, 2018; originally announced July 2018.

Comments: Accepted for Publication in MNRAS (9 pages, 3 figures)

Mapping Circumstellar Matter with Polarized Light - The Case of Supernova 2014J in M82

Authors: Yi Yang, Lifan Wang, Dietrich Baade, Peter J. Brown, Aleksandar Cikota, Misty Cracraft, Peter A. Hoflich, Justyn Maund, Ferdinando Patat, William B. Sparks, Jason Spyromilio, Heloise F. Stevance, Xiaofeng Wang, J. Craig Wheeler

Abstract: 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 at… ▽ More 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. △ Less

Submitted 14 January, 2018; originally announced January 2018.

Comments: 8 figures, 4 tables. Accepted for publication in ApJ

The evolution of the 3D shape of the broad-lined type Ic SN 2014ad

Authors: H. F. Stevance, J. R. Maund, D. Baade, P. Höflich, S. Howerton, F. Patat, M. Rose, J. Spyromilio, J. C. Wheeler, L. Wang

Abstract: 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 spect… ▽ More 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 λ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. △ Less

Submitted 20 April, 2017; originally announced April 2017.

Comments: Accepted for publication in MNRAS - 17 pages, 16 Figures

Late-time Flattening of Type Ia Supernova Light Curves: Constraints From SN 2014J in M82

Authors: Yi Yang, Lifan Wang, Dietrich Baade, Peter J. Brown, Aleksandar Cikota, Misty Cracraft, Peter A. Hoflich, Justyn Maund, Ferdinando Patat, William B. Sparks, Jason Spyromilio, Heloise F. Stevance, Xiaofeng Wang, J. Craig Wheeler

Abstract: 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… ▽ More 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. △ Less

Submitted 28 November, 2017; v1 submitted 5 April, 2017; originally announced April 2017.

Comments: 5 figures, 3 tables, accepted for publication in ApJ

Interstellar-Medium Mapping in M82 Through Light Echoes Around Supernova 2014J

Authors: Yi Yang, Lifan Wang, Dietrich Baade, Peter J. Brown, Misty Cracraft, Peter A. Hoflich, Justyn Maund, Ferdinando Patat, William B. Sparks, Jason Spyromilio, Heloise F. Stevance, Xiaofeng Wang, J. Craig Wheeler

Abstract: 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 com… ▽ More 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. △ Less

Submitted 7 November, 2016; v1 submitted 7 October, 2016; originally announced October 2016.

Comments: 45 pages, 16 figures, accepted for publication in The Astrophysical Journal

Spectropolarimetry of the Type IIb SN 2008aq

Authors: H. F. Stevance, J. R. Maund, D. Baade, P. Hoflich, F. Patat, J. Spyromilio, J. C. Wheeler, A. Clocchiatti, L. Wang, Y. Yang, P. Zelaya

Abstract: 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… ▽ More 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. △ Less

Submitted 4 March, 2020; v1 submitted 17 June, 2016; originally announced June 2016.

Comments: 7 pages, 6 figures, Accepted for publication in MNRAS - Erratum Added