Spectropolarimetry of stripped envelope core collapse supernovae

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 work 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. SN 2011hs was found to have a geometry consistent with an off-axis energy source within an ellipsoidal envelope, and features similar to SN 2011dh. In SN 2014ad we showed the presence of intermediate mass elements in outer parts of the ejecta as well as significant axi-symmetry, consistent with the presence of a jet. We also provided a re-analysis of the Type IIb SN 1993J, including a novel estimate of its interstellar polarisation. Diverse asymmetries were found in all SNe, both in the global geometry of the photosphere and in the distribution of the line forming regions. They are discussed in the context of previous studies. In order to simulate the geometry of the ejecta and the resulting observations, we re- created a pre-existing toy model, improving on previous work by devising a way to explore parameter space methodically. We found that such toy models are prone to degeneracies, and warn that they should be used with great caution, if at all. 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

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

We use the self-consistent stellar populations in the Binary Population A Spectral Synthesis (BPASS) models to assess whether NGC 1850-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

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

We present optical spectropolarimetry and spectroscopy of the broad-lined Type Ic (Ic-bl) SN 2014ad. Our spectropolarimetric observations cover seven epochs, from −2 to 66 d after V-band maximum, and the spectroscopic data were acquired from −2 to +107 d. The photospheric velocity estimates showed ejecta speeds similar to those of SN 1998bw and other SNe associated with gamma-ray bursts (GRBs). The spectropolarimetric data revealed aspherical outer ejecta and a nearly spherical interior. The polarization 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

Spectropolarimetry of the Type IIb SN 2008aq

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 Hα 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 signifi- cant levels of continuum polarisation at pcont= 0.70 (± 0.22) % in the first epoch, increasing to pcont=1.21 (± 0.33) % by the second epoch. Moreover, the presence of loops in the q − u planes of Hα and He I in the second epoch suggests a departure from axial symmetry

The shape of SN 1993J re-analysed

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 H β, He I λ5876, H α, He I λ6678, He I λ7065, and high velocity (HV) components of He I λ5876 and 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 H α 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

Linear spectropolarimetry of 35 Type Ia supernovae with VLT/FORS: an analysis of the Si ii line polarization

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 Focal Reducer and Low-Dispersion Spectrograph (FORS) on the Very Large Telescope at 127 epochs in total. We examined the polarization of the Si II λ6355 Å line (⁠pSiII⁠) 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 Δm15−pSiII relationship identified in a previous study, and show that subluminous and transitional objects display polarization values below the Δm15−pSiII 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 Δm15−pSiII relationship, may be explained in the context of a delayed-detonation model. In contrast, we compared our observations to numerical predictions in the Δm15−vSiII 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 displays 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

The 3D shape of Type IIb SN 2011hs

We observed seven epochs of spectropolarimetry in optical wavelengths for the Type IIb SN 2011hs, ranging from −3 to +40 d with respect to V-band maximum. A high degree of interstellar polarization was detected (up to ∼3 per cent), with a peak lying blueward of 4500 Å. Similar behaviours have been seen in some Type Ia supernovae (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 d (p = 0.55 ± 0.12 per cent and 0.75 ± 0.11 per cent, 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° in the continuum polarization angle between −2 and +18 d 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

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

We present the results from a high-cadence, multiwavelength observation campaign of AT 2016jbu (aka Gaia16cfr), an interacting transient. This data set complements the current literature by adding higher cadence as well as extended coverage of the light-curve 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 MV ∼-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 500 km s-1 seen in narrow emission features from a slow-moving CSM, and up to 10 000 km s-1 seen in broad absorption from some high-velocity material. Late-time spectra (∼+1 yr) 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 light curve, 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

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

We present the bolometric light curve, identification and analysis of the progenitor candidate, and preliminary modelling of AT 2016jbu (Gaia16cfr). We find a progenitor consistent with a ∼22-25 M⊙ 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 AT 2016jbu, supports a progenitor age of >10 Myr, consistent with a progenitor mass of ∼22 M⊙. A joint analysis of the velocity evolution of AT 2016jbu and the photospheric radius inferred from the bolometric light curve shows the transient is consistent with two successive outbursts/explosions. The first outburst ejected material with velocity ∼650 km s-1, while the second, more energetic event ejected material at ∼4500 km s-1. Whether the latter is the core collapse of the progenitor remains uncertain. We place a limit on the ejected 56Ni mass of <0.016 M⊙. Using the Binary Population And Spectral Synthesis (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 SuperNova Explosion Code (SNEC) to demonstrate that the low-energy explosions within some of these binary systems, together with sufficient circumstellar material, can reproduce the overall morphology of the light curve of AT 2016jbu

Heavy element production in a compact object merger observed by JWST

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)1, sources of high-frequency gravitational waves (GW)2 and likely production sites for heavy element nucleosynthesis via rapid neutron capture (the r-process)3. 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 mergers4–6, and contains a kilonova similar to AT2017gfo, associated with the gravitational-wave merger GW1708177–12. 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