The observation of a compact object with a mass of $2.50-2.67M_{\odot}$ on August 14, 2019, by the LIGO Scientific and Virgo collaborations (LVC) has the potential to improve our understanding of the supranuclear equation of state. While the gravitational-wave analysis of the LVC suggests that GW190814 likely was a binary black hole system, the secondary component could also have been the heaviest neutron star observed to date. We use our previously derived nuclear-physics--multi-messenger astrophysics framework to address the nature of this object. Based on our findings, we determine GW190814 to be a binary black hole merger with a probability of $>99.9\%$. Even if we weaken previously employed constraints on the maximum mass of neutron stars, the probability of a binary black hole origin is still $\sim 86\%$. Furthermore, we study the impact that this observation has on our understanding of the nuclear equation of state by analyzing the allowed region in the mass-radius diagram of neutron stars for both a binary black hole or neutron star - black hole scenario. We find that the unlikely scenario in which the secondary object was a neutron star requires rather stiff equations of state with a maximum speed of sound $c_s\geq \sqrt{0.6}$ times the speed of light, while the binary black hole scenario does not offer any new insight.

Early observations of Type Ia supernovae (SNe$\,$Ia) provide essential clues for understanding the progenitor system that gave rise to the terminal thermonuclear explosion. We present exquisite observations of SN$\,$2019yvq, the second observed SN$\,$Ia, after iPTF$\,$14atg, to display an early flash of emission in the ultraviolet (UV) and optical. Our analysis finds that SN$\,$2019yvq was unusual, even when ignoring the initial flash, in that it was moderately underluminous for an SN$\,$Ia ($M_g \approx -18.5\,$mag at peak) yet featured very high absorption velocities ($v \approx 15,000\,\mathrm{km\,s}^{-1}$ for Si II $\lambda$6355 at peak). We find that many of the observational features of SN$\,$2019yvq, aside from the flash, can be explained if the explosive yield of radioactive $^{56}\mathrm{Ni}$ is relatively low (we measure $M_{^{56}\mathrm{Ni}} = 0.31 \pm 0.05\,M_\odot$) and it and other iron-group elements are concentrated in the innermost layers of the ejecta. To explain both the UV/optical flash and peak properties of SN$\,$2019yvq we consider four different models: interaction between the SN ejecta and a nondegenerate companion, extended clumps of $^{56}\mathrm{Ni}$ in the outer ejecta, a double-detonation explosion, and the violent merger of two white dwarfs. Each of these models has shortcomings when compared to the observations; it is clear additional tuning is required to better match SN$\,$2019yvq. In closing, we predict that the nebular spectra of SN$\,$2019yvq will feature either H or He emission, if the ejecta collided with a companion, strong [Ca II] emission, if it was a double detonation, or narrow [O I] emission, if it was due to a violent merger.

Joint multi-messenger observations with gravitational waves and electromagnetic data offer new insights into the astrophysical studies of compact objects. The third Advanced LIGO and Advanced Virgo observing run began on April 1, 2019; during the eleven months of observation, there have been 14 compact binary systems candidates for which at least one component is potentially a neutron star. Although intensive follow-up campaigns involving tens of ground and space-based observatories searched for counterparts, no electromagnetic counterpart has been detected. Following on a previous study of the first six months of the campaign, we present in this paper the next five months of the campaign from October 2019 to March 2020. We highlight two neutron star - black hole candidates (S191205ah, S200105ae), two binary neutron star candidates (S191213g and S200213t) and a binary merger with a possible neutron star and a "MassGap" component, S200115j. Assuming that the gravitational-wave candidates are of astrophysical origin and their location was covered by optical telescopes, we derive possible constraints on the matter ejected during the events based on the non-detection of counterparts. We find that the follow-up observations during the second half of the third observing run did not meet the necessary sensitivity to constrain the source properties of the potential gravitational-wave candidate. Consequently, we suggest that different strategies have to be used to allow a better usage of the available telescope time. We examine different choices for follow-up surveys to optimize sky localization coverage vs.\ observational depth to understand the likelihood of counterpart detection.

On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. Preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical/near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS and VINROUGE projects also carried out a search on this event. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN) possibly generated by this NS-BH merger, and for the strategy of future searches. Altogether, our observations allow us to exclude a KN with large ejecta mass $M\gtrsim 0.1\,\mathrm{M_\odot}$ to a high ($>90\%$) confidence, and we can exclude much smaller masses in a subsample of our observations. This disfavours the tidal disruption of the neutron star during the merger. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundreds Mpc will be detected only by large facilities with both high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event.

The multi-messenger analysis of neutron-star mergers is not only of astrophysical interest, it also allows studying the behavior of supranuclear dense matter, probing the principles of general relativity, and measuring the expansion rate of the Universe. We perform a multi-messenger analysis of the gravitational-wave signal GW170817 and its electromagnetic counterparts AT2017gfo and GRB170817A, as well as the gravitational-wave signal GW190425. By incorporating information from the NICER observation of PSR J0030+0451, the radio observations of PSR J0740+6620, PSR J0348+4032, PSR J1614-2230, and nuclear-theory computations using chiral effective field theory, we provide new constraints on the neutron-star equation of state. Within our analysis, we determine the radius of a $1.4M_\odot$ NS to be $R_{1.4M_\odot}=11.74^{+0.98}_{-0.79}\ \rm km$ at $90\%$ confidence and estimate the Hubble constant to be $H_0=66.2^{+4.4}_{-4.2}\ \rm km /Mpc/s$ at $1\sigma$ uncertainty. The same multi-messenger-analysis techniques will be applicable for future observations, improving our understanding of the most extreme objects in the Universe.

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)

Time domain astronomy was revolutionised with the discovery of the first kilonova, AT2017gfo, in August 2017 which was associated with the gravitational wave signal GW170817. Since this event, numerous wide-field surveys have been optimising search strategies to maximise their efficiency of detecting these fast and faint transients. With the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS), we have been conducting a volume limited survey for intrinsically faint and fast fading events to a distance of $D\simeq200$ Mpc. Two promising candidates have been identified from this archival search, with sparse data - PS15cey and PS17cke. Here we present more detailed analysis and discussion of their nature. We observe that PS15cey was a luminous, fast declining transient at 320 Mpc. Models of BH-NS mergers with a very stiff equation of state could possibly reproduce the luminosity and decline but the physical parameters are extreme. A more likely scenario is that this was a SN2018kzr-like merger event. PS17cke was a faint, but fast declining, event at 15 Mpc. We explore several explosion scenarios of this transient including models of it as a NS-NS and BH-NS merger, the outburst of a massive luminous star, and compare it against other known fast fading transients. Although there is uncertainty in the explosion scenario due to difficulty in measuring the explosion epoch, we find PS17cke to be a plausible kilonova candidate from the model comparisons.

Colors of Type Ia supernovae in the first few days after explosion provide a potential discriminant between different models. In this paper, we present $g-r$ colors of 65 Type Ia supernovae discovered within 5 days from first light by the Zwicky Transient Facility in 2018, a sample that is about three times larger than that in the literature. We find that $g-r$ colors are intrinsically rather homogeneous at early phases, with about half of the dispersion attributable to photometric uncertainties ($\sigma_\mathrm{noise}\sim\sigma_\mathrm{int}\sim$ 0.18 mag). Colors are nearly constant starting from 6 days after first light ($g-r\sim-0.15$ mag), while the time evolution at earlier epochs is characterized by a continuous range of slopes, from events rapidly transitioning from redder to bluer colors (slope of $\sim-0.25$ mag day$^{-1}$) to events with a flatter evolution. The continuum in the slope distribution is in good agreement both with models requiring some amount of $^{56}$Ni mixed in the outermost regions of the ejecta and with "double-detonation" models having thin helium layers ($M_\mathrm{He}=0.01\,M_\odot$) and varying carbon-oxygen core masses. At the same time, six events show evidence for a distinctive "red bump" signature predicted by "double-detonation" models with larger helium masses. We finally identify a significant correlation between the early-time $g-r$ slopes and supernova brightness, with brighter events associated to flatter color evolution (p-value=0.006). The distribution of slopes, however, is consistent with being drawn from a single population, with no evidence for two components as claimed in the literature based on $B-V$ colors.

The detection of the binary neutron star merger GW170817 and the associated electromagnetic (EM) counterpart, the "kilonova" (kN) AT2017gfo, opened a new era in multi-messenger astronomy. However, despite many efforts, it has been proven very difficult to find additional EM counterparts to the many LIGO/Virgo compact binary merger candidates since then. We explore the sensitivity of ongoing optical surveys searching for kNe and propose ways to optimize the choices of filters and survey depth to boost the detection efficiency for these faint and fast-evolving transients. In particular, we use kN models to explore the dependence on ejecta mass, geometry, viewing angle, wavelength coverage and source distance. Considering that the O3 LIGO/Virgo binary neutron star merger candidates were at distances about 5 times further than GW170817, we find that efficient $gri$ counterpart search would require reaching a limiting magnitude $m_{\rm lim}=$ 23 mag, given the model phase-space. We conclude that kN searches during O3 were generally too shallow to observe optical counterparts, even under optimistic assumptions, thus providing a plausible explanation for the lack of detections.

We present our spectropolarimetric observations of SN 2017egm, a Type I superluminous supernova (SLSN-I) in a nearby galaxy NGC 3191, with Subaru telescope at 185.0 days after the g-band maximum light. This is the first spectropolarimetric observation for SLSNe at late phases. We find that the degree of the polarization in the late phase significantly changes from that measured at the earlier phase. The spectrum at the late phase shows a strong Ca emission line and therefore we reliably estimate the interstellar polarization component assuming that the emission line is intrinsically unpolarized. By subtracting the estimated interstellar polarization, we find that the intrinsic polarization at the early phase is only ~ 0.2 %, which indicates an almost spherical photosphere, with an axial ratio ~ 1.05. The intrinsic polarization at the late phase increases to ~ 0.8 %, which corresponds to the photosphere with an axial ratio 1.2. A nearly constant position angle of the polarization suggests the inner ejecta are almost axisymmetric. By these observations, we conclude that the inner ejecta are more aspherical than the outer ejecta. It may suggest the presence of a central energy source producing aspherical inner ejecta.

Synoptic searches for the optical counterpart to a binary neutron star (BNS) or neutron star-black hole (NSBH) merger can pose significant challenges towards the discovery of kilonovae and performing multi-messenger science. In this work, we describe the advantage of a global multi-telescope network towards this end, with a particular focus on the key and complementary role the Dark Energy Camera (DECam) plays in multi-facility follow-up. We describe the Global Relay of Observatories Watching Transients Happen (GROWTH) Target-of-Opportunity (ToO) Marshal, a common web application we built to ingest events, plan observations, search for transient candidates, and retrieve performance summary statistics for all of the telescopes in our network. Our infrastructure enabled us to conduct observations of two events during O3a, S190426c and S190510g. Furthermore, our analysis of deep DECam observations of S190814bv conducted by the DESGW team, and access to a variety of global follow-up facilities allowed us to place meaningful constraints on the parameters of the kilonova and the merging binary. We emphasize the importance of a global telescope network in conjunction with a power telescope like DECam in performing searches for the counterparts to gravitational-wave sources.

Growing evidence suggests that Type Iax supernovae might be the result of thermonuclear deflagrations of Chandrasekhar-mass white dwarfs in binary systems. We carry out Monte Carlo radiative transfer simulations and predict spectropolarimetric features originating from the supernova explosion and subsequent ejecta interaction with the companion star. Specifically, we calculate viewing-angle dependent flux and polarisation spectra for a 3D model simulating the deflagration of a Chandrasekhar-mass white dwarf and, for a second model, simulating the ejecta interaction with a main-sequence star. We find that the intrinsic signal is weakly polarised and only mildly viewing-angle dependent, owing to the overall spherical symmetry of the explosion and the depolarising contribution of iron-group elements dominating the ejecta composition. The interaction with the companion star carves out a cavity in the ejecta and produces a detectable, but modest signal that is significant only at relatively blue wavelengths ($\lesssim$ 5000 $\unicode{x212B}$). In particular, increasingly fainter and redder spectra are predicted for observer orientations further from the cavity, while a modest polarisation signal $P\sim0.2$ per cent is found at blue wavelengths for orientations 30$^\circ$ and 45$^\circ$ away from the cavity. We find a reasonable agreement between the interaction model viewed from these orientations and spectropolarimetric data of SN 2005hk and interpret the maximum-light polarisation signal seen at blue wavelengths for this event as a possible signature of the ejecta-companion interaction. We encourage further polarimetric observations of SNe Iax to test whether our results can be extended and generalised to the whole SN Iax class.

The detection of GW170817 is revolutionizing many areas of astrophysics with the joint observation of gravitational waves and electromagnetic emissions. These multi-messenger events provide a new approach to determine the Hubble constant, thus, they are a promising candidate for mitigating the tension between measurements of Type Ia supernovae via the local distance ladder and the Cosmic Microwave Background. In addition to the "standard siren" provided by the gravitational-wave measurement, the kilonova itself has characteristics that allow to improve existing measurements or to perform yet another, independent measurement of the Hubble constant without gravitational-wave information. Here, we employ standardization techniques borrowed from the type-Ia community and apply them to kilonovae, not using any information from the gravitational-wave signal. We use two versions of this technique, one derived from direct observables measured from the lightcurve, and the other based on inferred ejecta parameters, e.g., mass, velocity, and composition, for two different models. These lead to Hubble Constant measurements of $H_0 = 109^{+49}_{-35}$\,km $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$ for the measured analysis, and $H_0 = 85^{+22}_{-17}$\,km $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$ and $H_0 = 79^{+23}_{-15}$\,km $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$ for the inferred analyses. This measurement has error bars within ~$\sim$\,2 to the gravitational-wave measurements ($H_0=74^{+16}_{-8}$\,km $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$), showing its promise as an independent constraint on $H_0$.

While it is clear that Type Ia supernovae (SNe) are the result of thermonuclear explosions in C/O white dwarfs (WDs), a great deal remains uncertain about the binary companion that facilitates the explosive disruption of the WD. Here, we present a comprehensive analysis of a unique, and large, data set of 127 SNe Ia with exquisite coverage by the Zwicky Transient Facility (ZTF). High-cadence (6 observations per night) ZTF observations allow us to measure the SN rise time and examine its initial evolution. We develop a Bayesian framework to model the early rise as a power-law in time, which enables the inclusion of priors in our model. For a volume-limited subset of normal SNe Ia, we find the mean power-law index is consistent with 2 in the $r_\mathrm{ztf}$-band ($\alpha_r = 2.01\pm0.02$), as expected in the expanding fireball model. There are, however, individual SNe that are clearly inconsistent with $\alpha_r=2$. We estimate a mean rise time of 18.5$\,$d (with a range extending from $\sim$15$-$22$\,$d), though this is subject to the adopted prior. We identify an important, previously unknown, bias whereby the rise times for higher redshift SNe within a flux-limited survey are systematically underestimated. This effect can be partially alleviated if the power-law index is fixed to $\alpha=2$, in which case we estimate a mean rise time of 21.0$\,$d (with a range from $\sim$18$-$23$\,$d). The sample includes a handful or rare and peculiar SNe Ia. Finally, we conclude with a discussion of lessons learned from the ZTF sample that can eventually be applied to Large Synoptic Survey Telescope observations.

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $\kappa < 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ \chi < 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv.

GW170817 showed that neutron star mergers not only emit gravitational waves but also can release electromagnetic signatures in multiple wavelengths. Within the first half of the third observing run of the Advanced LIGO and Virgo detectors, there have been a number of gravitational wave candidates of compact binary systems for which at least one component is potentially a neutron star. In this article, we look at the candidates S190425z, S190426c, S190510g, S190901ap, and S190910h, predicted to have potentially a non-zero remnant mass, in more detail. All these triggers have been followed up with extensive campaigns by the astronomical community doing electromagnetic searches for their optical counterparts; however, according to the released classification, there is a high probability that some of these events might not be of extraterrestrial origin. Assuming that the triggers are caused by a compact binary coalescence and that the individual source locations have been covered during the EM follow-up campaigns, we employ three different kilonova models and apply them to derive possible constraints on the matter ejection consistent with the publicly available gravitational-wave trigger information and the lack of a kilonova detection. These upper bounds on the ejecta mass can be related to limits on the maximum mass of the binary neutron star candidate S190425z and to constraints on the mass-ratio, spin, and NS compactness for the potential black hole-neutron star candidate S190426c. Our results show that deeper electromagnetic observations for future gravitational wave events near the horizon limit of the advanced detectors are essential.

We extend the range of validity of the ARTIS 3D radiative transfer code up to hundreds of days after explosion, when Type Ia supernovae are in their nebular phase. To achieve this, we add a non-local thermodynamic equilibrium (non-LTE) population and ionisation solver, a new multi-frequency radiation field model, and a new atomic dataset with forbidden transitions. We treat collisions with non-thermal leptons resulting from nuclear decays to account for their contribution to excitation, ionisation, and heating. We validate our method with a variety of tests including comparing our synthetic nebular spectra for the well-known one-dimensional W7 model with the results of other studies. As an illustrative application of the code, we present synthetic nebular spectra for the detonation of a sub-Chandrasekhar white dwarf in which the possible effects of gravitational settling of Ne22 prior to explosion have been explored. Specifically, we compare synthetic nebular spectra for a 1.06 M$_\odot$ white dwarf model obtained when 5.5 Gyr of very-efficient settling is assumed to a similar model without settling. We find that this degree of Ne22 settling has only a modest effect on the resulting nebular spectra due to increased Ni58 abundance. Due to the high ionisation in sub-Chandrasekhar models, the nebular [Ni II] emission remains negligible, while the [Ni III] line strengths are increased and the overall ionisation balance is slightly lowered in the model with Ne22 settling. In common with previous studies of sub-Chandrasekhar models at nebular epochs, these models overproduce [Fe III] emission relative to [Fe II] in comparison to observations of normal Type Ia supernovae.

There is a strong degeneracy between the luminosity distance ($D_L$) and the observer viewing angle ($\theta_{\rm obs}$; hereafter viewing angle) of the gravitational wave (GW) source with an electromagnetic counterpart, GW170817. Here, for the first time, we present independent constraints on $\theta_{\rm obs} = 32.5 ^{\circ +11.7}_{-9.7}$ from broad-band photometry of the kilonova (kN) AT2017gfo associated with GW170817. These constraints are consistent with independent results presented in the literature using the associated gamma ray burst GRB170817A. Combining the constraints on $\theta_{\rm obs}$ with the GW data, we find an improvement of 24$\%$ on $H_0$. The observer angle constraints are insensitive to other model parameters, e.g. the ejecta mass, half-opening angle of the lanthanide-rich region and the temperature. A broad wavelength coverage extending to the near infrared is helpful to robustly constrain $\theta_{\rm obs}$. While the improvement on $H_0$ presented here is smaller than the one from high angular resolution imaging of the radio counterpart of GW170817, kN observations are significantly more feasible at the typical distances of such events from current and future LIGO-Virgo Collaboration observing runs ($D_L \sim 100$ Mpc). Our results are insensitive to the assumption on the peculiar velocity of the kN host galaxy.

We present POSSIS, a time-dependent three-dimensional Monte Carlo code for modelling radiation transport in supernovae and kilonovae. The code incorporates wavelength- and time-dependent opacities and predicts viewing-angle dependent spectra, light curves and polarization for both idealized and hydrodynamical explosion models. We apply the code to a kilonova model with two distinct ejecta components, one including lanthanide elements with relatively high opacities and the other devoid of lanthanides and characterized by lower opacities. We find that a model with total ejecta mass $M_\mathrm{ej}=0.04\,M_\odot$ and half-opening angle of the lanthanide-rich component $\Phi=30^\circ$ provides a good match to GW 170817 / AT 2017gfo for orientations near the polar axis (i.e. for a system viewed close to face-on). We then show how crucial is the use of self-consistent multi-dimensional models in place of combining one-dimensional models to infer important parameters as the ejecta masses. We finally explore the impact of $M_\mathrm{ej}$ and $\Phi$ on the synthetic observables and highlight how the relatively fast computation times of POSSIS make it well-suited to perform parameter-space studies and extract key properties of supernovae and kilonovae. Spectra calculated with POSSIS in this and future studies will be made publicly available.

Early-time observations of Type Ia supernovae (SNe Ia) are essential to constrain their progenitor properties. In this paper, we present high-quality light curves of 127 SNe Ia discovered by the Zwicky Transient Facility (ZTF) in 2018. We describe our method to perform forced point spread function (PSF) photometry, which can be applied to other types of extragalactic transients. With a planned cadence of six observations per night ($3g+3r$), all of the 127 SNe Ia are detected in both $g$ and $r$ band more than 10\,d (in the rest frame) prior to the epoch of $g$-band maximum light. The redshifts of these objects range from $z=0.0181$ to 0.165; the median redshift is 0.074. Among the 127 SNe, 50 are detected at least 14\,d prior to maximum light (in the rest frame), with a subset of 9 objects being detected more than 17\,d before $g$-band peak. This is the largest sample of young SNe Ia collected to date; it can be used to study the shape and color evolution of the rising light curves in unprecedented detail. We discuss six peculiar events in this sample, including one 02cx-like event ZTF18abclfee (SN\,2018crl), one Ia-CSM SN ZTF18aaykjei (SN\,2018cxk), and four objects with possible super-Chandrasekhar mass progenitors: ZTF18abhpgje (SN\,2018eul), ZTF18abdpvnd (SN\,2018dvf), ZTF18aawpcel (SN\,2018cir) and ZTF18abddmrf (SN\,2018dsx).

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

We present a comprehensive dataset of optical and near-infrared photometry and spectroscopy of type~Ia supernova (SN) 2016hnk, combined with integral field spectroscopy (IFS) of its host galaxy, MCG -01-06-070, and nearby environment. Properties of the SN local environment are characterized by means of single stellar population synthesis applied to IFS observations taken two years after the SN exploded. SN 2016hnk spectra are compared to other 1991bg-like SNe Ia, 2002es-like SNe Ia, and Ca-rich transients. In addition, abundance stratification modelling is used to identify the various spectral features in the early phase spectral sequence and the dataset is also compared to a modified non-LTE model previously produced for the sublumnious SN 1999by. SN 2016hnk is consistent with being a sub-luminous (M$_{\rm B}=-16.7$ mag, s$_{\rm BV}$=0.43$\pm$0.03), highly reddened object. IFS of its host galaxy reveals both a significant amount of dust at the SN location, as well as residual star formation and a high proportion of old stellar populations in the local environment compared to other locations in the galaxy, which favours an old progenitor for SN 2016hnk. Inspection of a nebular spectrum obtained one year after maximum contains two narrow emission lines attributed to the forbidden [Ca II] $\lambda\lambda$7291,7324 doublet with a Doppler shift of 700 km s$^{-1}$. Based on various observational diagnostics, we argue that the progenitor of SN 2016hnk was likely a near Chandrasekhar-mass ($M_{\rm Ch}$) carbon-oxygen white dwarf that produced 0.108 $M_\odot$ of $^{56}$Ni. Our modeling suggests that the narrow [Ca II] features observed in the nebular spectrum are associated with $^{48}$Ca from electron capture during the explosion, which is expected to occur only in white dwarfs that explode near or at the $M_{\rm Ch}$ limit.

Type IIP supernovae (SNe IIP), which represent the most common class of core-collapse (CC) SNe, show a rapid increase in continuum polarization just after entering the tail phase. This feature can be explained by a highly asymmetric helium core, which is exposed when the hydrogen envelope becomes transparent. Here we report the case of a SN IIP (SN~2017gmr) that shows an unusually early rise of the polarization, $\gtrsim 30$ days before the start of the tail phase. This implies that SN~2017gmr is an SN IIP that has very extended asphericity. The asymmetries are not confined to the helium core, but reach out to a significant part of the outer hydrogen envelope, hence clearly indicating a marked intrinsic diversity in the aspherical structure of CC explosions. These observations provide new constraints on the explosion mechanism, where viable models must be able to produce such extended deviations from spherical symmetry, and account for the observed geometrical diversity.

The first two months of the third Advanced LIGO and Virgo observing run (2019 April-May) showed that distant gravitational wave (GW) events can now be readily detected. Three candidate mergers containing neutron stars (NS) were reported in a span of 15 days, all likely located more than 100 Mpc away. However, distant events such as the three new NS mergers are likely to be coarsely localized, which highlights the importance of facilities and scheduling systems that enable deep observations over hundreds to thousands of square degrees to detect the electromagnetic counterparts. On 2019-05-10 02:59:39.292 UT the GW candidate S190510g was discovered and initially classified as a BNS merger with 98% probability. The GW event was localized within an area of 3462 deg2, later refined to 1166 deg2 (90%) at a distance of 227 +- 92 Mpc. We triggered Target of Opportunity observations with the Dark Energy Camera (DECam), a wide-field optical imager mounted at the prime focus of the 4m Blanco Telescope at CTIO in Chile. This Letter describes our DECam observations and our real-time analysis results, focusing in particular on the design and implementation of the observing strategy. Within 24 hours of the merger time, we observed 65% of the total enclosed probability of the final skymap with an observing efficiency of 94%. We identified and publicly announced 13 candidate counterparts. S190510g was re-classified 1.7 days after the merger, after our observations were completed, with a "binary neutron star merger" probability reduced from 98% to 42% in favor of a "terrestrial" classification.

We present high-cadence ultraviolet (UV), optical, and near-infrared (NIR) data on the luminous Type II-P supernova SN 2017gmr from hours after discovery through the first 180 days. SN 2017gmr does not show signs of narrow, high-ionization emission lines in the early optical spectra, yet the optical lightcurve evolution suggests that an extra energy source from circumstellar medium (CSM) interaction must be present for at least 2 days after explosion. Modeling of the early lightcurve indicates a ~500R$_{\odot}$ progenitor radius, consistent with a rather compact red supergiant, and late-time luminosities indicate up to 0.130 $\pm$ 0.026 M$_{\odot}$ of $^{56}$Ni are present, if the lightcurve is solely powered by radioactive decay, although the $^{56}$Ni mass may be lower if CSM interaction contributes to the post-plateau luminosity. Prominent multi-peaked emission lines of H$\alpha$ and [O I] emerge after day 154, as a result of either an asymmetric explosion or asymmetries in the CSM. The lack of narrow lines within the first two days of explosion in the likely presence of CSM interaction may be an example of close, dense, asymmetric CSM that is quickly enveloped by the spherical supernova ejecta.

Both multi-messenger astronomy and new high-throughput wide-field surveys require flexible tools for the selection and analysis of astrophysical transients. We here introduce the Alert Management, Photometry and Evaluation of Lightcurves (AMPEL) system, an analysis framework designed for high-throughput surveys and suited for streamed data. AMPEL combines the functionality of an alert broker with a generic framework capable of hosting user-contributed code, that encourages provenance and keeps track of the varying information states that a transient displays. The latter concept includes information gathered over time and data policies such as access or calibration levels. We describe a novel ongoing real-time multi-messenger analysis using AMPEL to combine IceCube neutrino data with the alert streams of the Zwicky Transient Facility (ZTF). We also reprocess the first four months of ZTF public alerts, and compare the yields of more than 200 different transient selection functions to quantify efficiencies for selecting Type Ia supernovae that were reported to the Transient Name Server (TNS). We highlight three channels suitable for (1) the collection of a complete sample of extragalactic transients, (2) immediate follow-up of nearby transients and (3) follow-up campaigns targeting young, extragalactic transients. We confirm ZTF completeness in that all TNS supernovae positioned on active CCD regions were detected. AMPEL can assist in filtering transients in real time, running alert reaction simulations, the reprocessing of full datasets as well as in the final scientific analysis of transient data. This text introduces how users can design their own channels for inclusion in the AMPEL live instance that parses the ZTF stream and the real-time submission of high quality extragalactic supernova candidates to the TNS.

The Zwicky Transient Facility (ZTF) is a new optical time-domain survey that uses the Palomar 48-inch Schmidt telescope. A custom-built wide-field camera provides a 47 deg$^2$ field of view and 8 second readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory (PTF). We describe the design and implementation of the camera and observing system. The ZTF data system at the Infrared Processing and Analysis Center provides near-real-time reduction to identify moving and varying objects. We outline the analysis pipelines, data products, and associated archive. Finally, we present on-sky performance analysis and first scientific results from commissioning and the early survey. ZTF's public alert stream will serve as a useful precursor for that of the Large Synoptic Survey Telescope.

The Zwicky Transient Facility (ZTF), a public-private enterprise, is a new time domain survey employing a dedicated camera on the Palomar 48-inch Schmidt telescope with a 47 deg$^2$ field of view and 8 second readout time. It is well positioned in the development of time domain astronomy, offering operations at 10% of the scale and style of the Large Synoptic Survey Telescope (LSST) with a single 1-m class survey telescope. The public surveys will cover the observable northern sky every three nights in g and r filters and the visible Galactic plane every night in g and r. Alerts generated by these surveys are sent in real time to brokers. A consortium of universities which provided funding ("partnership") are undertaking several boutique surveys. The combination of these surveys producing one million alerts per night allows for exploration of transient and variable astrophysical phenomena brighter than r $\sim$ 20.5 on timescales of minutes to years. We describe the primary science objectives driving ZTF including the physics of supernovae and relativistic explosions, multi-messenger astrophysics, supernova cosmology, active galactic nuclei and tidal disruption events, stellar variability, and Solar System objects.

We present the best 265 sampled R-band light curves of spectroscopically identified Type Ia supernovae (SNe) from the Palomar Transient Factory (PTF; 2009-2012) survey and the intermediate Palomar Transient Factory (iPTF; 2013-2017). A model-independent light curve template is built from our data-set with the purpose to investigate average properties and diversity in our sample. We searched for multiple populations in the light curve properties using machine learning tools. We also utilised the long history of our light curves, up to 4000 days, to exclude any significant pre- or post- supernova flares. From the shapes of light curves we found the average rise time in the R band to be $16.8^{+0.5}_{-0.6}$ days. Although PTF/iPTF were single-band surveys, by modelling the residuals of the SNe in the Hubble-Lema\^{i}tre diagram, we estimate the average colour excess of our sample to be $<$E$($B$-$V$)> \approx 0.05(2)$ mag and thus the mean corrected peak brightness to be $M_R = -19.02\pm0.02$ $+5 \log( {\rm H}_0 [{\rm km} \cdot{\rm s}^{-1} {\rm Mpc}^{-1}]/70)$ mag with only weakly dependent on light curve shape. The intrinsic scatter is found to be $\sigma_R = 0.186 \pm 0.033$ mag for the redshift range $0.05<z<0.1$, without colour corrections of individual SNe. Our analysis shows that Malmquist bias becomes very significant at z=0.13. A similar limitation is expected for the ongoing Zwicky Transient Facility (ZTF) survey using the same telescope, but new camera expressly designed for ZTF.

The Gravitational Wave (GW) event GW 170817 was generated by the coalescence of two neutron stars (NS) and produced an electromagnetic transient, labelled AT 2017gfo, that was target of a massive observational campaign. Polarimetry, a powerful diagnostic tool for probing the geometry and emission processes of unresolved sources, was obtained for this event. The observed linear polarization was consistent with being mostly induced by intervening dust, suggesting that the intrinsic emission was weakly polarized ($P < 0.4-0.5$ %). In this paper, we present and discuss a detailed analysis of the linear polarization expected from a merging NS binary system by means of 3D Monte Carlo radiative transfer simulations assuming a range of possible configurations, wavelengths, epochs and viewing angles. We find that polarization originates from the non-homogeneous opacity distribution within the ejecta and can reach levels of $P\sim1$ % at early times (1$-$2 days after the merger) and in the optical R band. Smaller polarization signals are expected at later epochs and/or different wavelengths. From the viewing-angle dependence of the polarimetric signal, we constrain the observer orientation of AT 2017gfo within $\sim$65$^\circ$ from the polar direction. The detection of non-zero polarization in future events will unambiguously reveal the presence of a lanthanide-free ejecta component and unveil its spatial and angular distribution.

Super-luminous supernovae (SLSNe) are rare events defined as being significantly more luminous than normal terminal stellar explosions. The source of the extra powering needed to achieve such luminosities is still unclear. Discoveries in the local Universe (i.e. $z<0.1$) are scarce, but afford dense multi-wavelength observations. Additional low-redshift objects are therefore extremely valuable. We present early-time observations of the type I SLSN ASASSN-18km/SN~2018bsz. These data are used to characterise the event and compare to literature SLSNe and spectral models. Host galaxy properties are also analysed. Optical and near-IR photometry and spectroscopy were analysed. Early-time ATLAS photometry was used to constrain the rising light curve. We identified a number of spectral features in optical-wavelength spectra and tracked their time evolution. Finally, we used archival host galaxy photometry together with HII region spectra to constrain the host environment. ASASSN-18km/SN~2018bsz is found to be a type I SLSN in a galaxy at a redshift of 0.0267 (111 Mpc), making it the lowest-redshift event discovered to date. Strong CII lines are identified in the spectra. Spectral models produced by exploding a Wolf-Rayet progenitor and injecting a magnetar power source are shown to be qualitatively similar to ASASSN-18km/SN~2018bsz, contrary to most SLSNe-I that display weak/non-existent CII lines. ASASSN-18km/SN~2018bsz displays a long, slowly rising, red 'plateau' of $>$26 days, before a steeper, faster rise to maximum. The host has an absolute magnitude of --19.8 mag ($r$), a mass of M$_{*}$ = 1.5$^{+0.08}_{-0.33}$ $\times$10$^{9}$ M$_{\odot}$ , and a star formation rate of = 0.50$^{+2.22}_{-0.19}$ M$_{\odot}$ yr$^{-1}$. A nearby HII region has an oxygen abundance (O3N2) of 8.31$\pm$0.01 dex.

Key information about the progenitor system and the explosion mechanism of Type Ia supernovae (SNe~Ia) can be obtained from early observations, within a few days from explosion. iPTF16abc was discovered as a young SN~Ia with excellent early time data. Here, we present photometry and spectroscopy of the SN in the nebular phase. A comparison of the early time data with a sample of SNe~Ia shows distinct features, differing from normal SNe~Ia at early phases but similar to normal SNe~Ia at a few weeks after maximum light (i.e. the transitional phase) and well into the nebular phase. The transparency timescales ($t_0$) for this sample of SNe~Ia range between $\sim$ 25 and 41 days indicating a diversity in the ejecta masses. $t_0$ also weakly correlates with the peak bolometric luminosity, consistent with the interpretation that SNe with higher ejecta masses would produce more $^{56}$Ni. Comparing the $t_0$ and the maximum luminosity, L$_{max}$\, distribution of a sample of SNe~Ia to predictions from a wide range of explosion models we find an indication that the sub-Chandrasekhar mass models span the range of observed values. However, the bright end of the distribution can be better explained by Chandrasekhar mass delayed detonation models, hinting at multiple progenitor channels to explain the observed bolometric properties of SNe~Ia. iPTF16abc appears to be consistent with the predictions from the M$_{ch}$ models.

The colour evolution of reddened Type Ia supernovae can place strong constraints on the location of dust and help address the question of whether the observed extinction stems from the interstellar medium or from circumstellar material surrounding the progenitor. Here we analyse BV photometry of 48 reddened Type Ia supernovae from the literature and estimate the dust location from their $B-V$ colour evolution. We find a time-variable colour excess $E(B-V)$ for 15 supernovae in our sample and constrain dust to distances between 0.013 and 45 pc ($4\times10^{16}-10^{20}$ cm). For the remaining supernovae, we obtain a constant $E(B-V)$ evolution and place lower limits on the dust distance from the explosion. In all the 48 supernovae, the inferred dust location is compatible with an interstellar origin for the extinction. This is corroborated by the observation that supernovae with relatively nearby dust ($\lesssim$ 1 pc) are located close to the center of their host galaxy, in high-density dusty regions where interactions between the supernova radiation and interstellar clouds close by are likely to occur. For supernovae showing time-variable $E(B-V)$, we identify a potential preference for low $R_\mathrm{V}$ values, unusually strong sodium absorption and blue-shifted and time-variable absorption features. Within the interstellar framework, this brings evidence to a proposed scenario where cloud-cloud collisions induced by the supernova radiation pressure can shift the grain size distribution to smaller values and enhance the abundance of sodium in the gaseous phase.

Superluminous supernovae (SLSNe) are at least $\sim$5 times more luminous than common supernovae (SNe). Especially hydrogen-poor SLSN-I are difficult to explain with conventional powering mechanisms. One possible scenario that might explain such luminosities is that SLSNe-I are powered by an internal engine, such as a magnetar or an accreting black hole. Strong magnetic fields or collimated jets can circularly polarize light. In this work, we measured circular polarization of two SLSNe-I with the FOcal Reducer and low dispersion Spectrograph (FORS2) mounted at the ESO's Very Large Telescope (VLT). PS17bek, a fast evolving SLSN-I, was observed around peak, while OGLE16dmu, a slowly evolving SLSN-I, was observed 100 days after maximum. Neither SLSN shows evidence of circularly polarized light, however, these non-detections do not rule out the magnetar scenario as the powering engine for SLSNe-I. We calculate the strength of the magnetic field and the expected circular polarization as a function of distance from the magnetar, which decreases very fast. Additionally, we observed no significant linear polarization for PS17bek at four epochs, suggesting that the photosphere near peak is close to spherical symmetry.

We performed a wide-area (2000 deg$^{2}$) g and I band experiment as part of a two month extension to the Intermediate Palomar Transient Factory. We discovered 36 extragalactic transients including iPTF17lf, a highly reddened local SN Ia, iPTF17bkj, a new member of the rare class of transitional Ibn/IIn supernovae, and iPTF17be, a candidate luminous blue variable outburst. We do not detect any luminous red novae and place an upper limit on their rate. We show that adding a slow-cadence I band component to upcoming surveys such as the Zwicky Transient Facility will improve the photometric selection of cool and dusty transients.

Nearby type Ia supernovae (SNe Ia), such as SN 2017cbv, are useful events to address the question of what the elusive progenitor systems of the explosions are. Hosseinzadeh et al. (2017) suggested that the early blue excess of the lightcurve of SN 2017cbv could be due to the supernova ejecta interacting with a non-degenerate companion star. Some SN Ia progenitor models suggest the existence of circumstellar (CS) environments in which strong outflows create low density cavities of different radii. Matter deposited at the edges of the cavities, should be at distances at which photoionisation due to early ultraviolet (UV) radiation of SNe Ia causes detectable changes to the observable Na I D and Ca II H&K absorption lines. To study possible narrow absorption lines from such material, we obtained a time-series of high-resolution spectra of SN 2017cbv at phases between $-14.8$ and $+83$ days with respect to $B$-band maximum, covering the time at which photoionisation is predicted to occur. Both narrow Na I D and Ca II H&K are detected in all spectra, with no measurable changes between the epochs. We use photoionisation models to rule out the presence of Na I and Ca II gas clouds along the line-of-sight of SN 2017cbv between $\sim8\times10^{16}$--$2\times10^{19}$~cm and $\sim10^{15}$--$10^{17}$~cm, respectively. Assuming typical abundances, the mass of a homogenous spherical CS gas shell with radius $R$ must be limited to $M^{\rm CSM}_{\rm HI}<3\times10^{-4}\times(R/10^{17}[{\rm cm}])^2$ M$_{\odot}$. The bounds point to progenitor models that deposit little gas in their CS environment.

Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa <= 0.5 cm2/gm at a velocity of 0.2 +/- 0.1c. The power source is constrained to have a power law slope of beta = -1.2 +/- 0.3, consistent with radioactive powering from r-process nuclides. We identify line features in the spectra that are consistent with light r-process elements (90 < A < 140). As it fades, the transient rapidly becomes red, and emission may have contribution by a higher opacity, lanthanide-rich ejecta component. This indicates that neutron star mergers produce gravitational waves, radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements.

We present a new technique to infer dust locations towards reddened Type Ia supernovae and to help discriminate between an interstellar and a circumstellar origin for the observed extinction. Using Monte Carlo simulations, we show that the time-evolution of the light-curve shape and especially of the colour excess \ebv~places strong constraints on the distance between dust and the supernova. We apply our approach to two highly-reddened Type Ia supernovae for which dust distance estimates are available in the literature: SN 2006X and SN 2014J. For the former, we obtain a time-variable $E(B-V)$ and from this derive a distance of $27.5^{+9.0}_{-4.9}$ or $22.1^{+6.0}_{-3.8}$ pc depending on whether dust properties typical of the Large Magellanic Cloud (LMC) or the Milky Way (MW) are used. For the latter, instead, we obtain a constant $E(B-V)$ consistent with dust at distances larger than 50 and 38 pc for LMC$-$ and MW$-$type dust, respectively. Values thus extracted are in excellent agreement with previous estimates for the two supernovae. Our findings suggest that dust responsible for the extinction towards these supernovae is likely to be located within interstellar clouds. We also discuss how other properties of reddened Type Ia supernovae $-$ such as their peculiar extinction and polarization behaviour and the detection of variable, blue-shifted sodium features in some of these events $-$ might be compatible with dust and gas at interstellar-scale distances.

Context. Type Ia supernovae (SNe Ia) can be used to address numerous questions in astrophysics and cosmology. Due to their well known spectral and photometric properties, SNe Ia are well suited to study gas and dust along the lines-of-sight to the explosions. For example, narrow Na I D and Ca II H&K absorption lines can be studied easily, because of the well-defined spectral continuum of SNe Ia around these features. Aims. We study the gas and dust along the line-of-sight to iPTF16abc, which occurred in an unusual location, in a tidal arm, 80 kpc from centre of the galaxy NGC 5221. Methods. Using a time-series of high-resolution spectra, we examine narrow Na I D and Ca II H&K absorption features for variations in time, which would be indicative for circumstellar (CS) matter. Furthermore, we take advantage of the well known photometric properties of SNe Ia to determine reddening due to dust along the line-of-sight. Results. From the lack of variations in Na I D and Ca II H&K, we determine that none of the detected absorption features originate from the CS medium of iPTF16abc. While the Na I D and Ca II H&K absorption is found to be optically thick, a negligible amount of reddening points to a small column of interstellar dust. Conclusions. We find that the gas along the line-of-sight to iPTF16abc is typical of what might be found in the interstellar medium (ISM) within a galaxy. It suggests that we are observing gas that has been tidally stripped during an interaction of NGC 5221 with one of its neighbouring galaxies in the past $\sim10^9$ years. In the future, the gas clouds could become the locations of star formation. On a longer time scale, the clouds might diffuse, enriching the circum-galactic medium (CGM) with metals. The gas profile along the line-of-sight should be useful for future studies of the dynamics of the galaxy group containing NGC 5221.

The light from distant supernovae (SNe) can be magnified through gravitational lensing when a foreground galaxy is located along the line of sight. This line-up allows for detailed studies of SNe at high redshift that otherwise would not be possible. Spectroscopic observations of lensed high-redshift Type Ia supernovae (SNe Ia) are of particular interest since they can be used to test for evolution of their intrinsic properties. The use of SNe Ia for probing the cosmic expansion history has proven to be an extremely powerful method for measuring cosmological parameters. However, if systematic redshift-dependent properties are found, their usefulness for future surveys could be challenged. We investigate whether the spectroscopic properties of the strongly lensed and very distant SN Ia PS1-10afx at $z=1.4$ deviates from the well-studied populations of normal SNe Ia at nearby or intermediate distance. We created median spectra from nearby and intermediate-redshift spectroscopically normal SNe Ia from the literature at -5 and +1 days from light-curve maximum. We then compared these median spectra to those of PS1-10afx. We do not find signs of spectral evolution in PS1-10afx. The observed deviation between PS1-10afx and the median templates are within what is found for SNe at low- and intermediate-redshift. There is a noticeable broad feature centred at $\rm \lambda\sim 3500$~\AA{}, which is present only to a lesser extent in individual low and intermediate redshift SN Ia spectra. From a comparison with a recently developed explosion model, we find this feature to be dominated by iron peak elements, in particular, singly ionized cobalt and chromium.

We present the first spectropolarimetric observations of a hydrogen-free superluminous supernova at z=0.1136, namely SN 2015bn. The transient shows significant polarization at both the observed epochs: one 24 days before maximum light in the rest-frame, and the subsequent at 27 days after peak luminosity. Analysis of the Q-U plane suggests the presence of a dominant axis and no physical departure from the main axis at either epoch. The polarization spectrum along the dominant axis is characterized by a strong wavelength dependence and an increase in the signal from the first to the second epoch. We use a Monte Carlo code to demonstrate that these properties are consistent with a simple toy model that adopts an axi-symmetric ellipsoidal configuration for the ejecta. We find that the wavelength dependence of the polarisation is possibly due to a strong wavelength dependence in the line opacity, while the higher level of polarisation at the second epoch is a consequence of the increase in the asphericity of the inner layers of the ejecta or the fact that the photosphere recedes into less spherical layers. The geometry of the superluminous supernova results similar to those of stripped-envelope core collapse SNe connected with GRB, while the overall evolution of the ejecta shape could be consistent with a central engine.

Calculations of synthetic spectropolarimetry are one means to test multi-dimensional explosion models for Type Ia supernovae. In a recent paper, we demonstrated that the violent merger of a 1.1 and 0.9 M$_{\odot}$ white dwarf binary system is too asymmetric to explain the low polarization levels commonly observed in normal Type Ia supernovae. Here, we present polarization simulations for two alternative scenarios: the sub-Chandrasekhar mass double-detonation and the Chandrasekhar mass delayed-detonation model. Specifically, we study a two-dimensional double-detonation model and a three-dimensional delayed-detonation model, and calculate polarization spectra for multiple observer orientations in both cases. We find modest polarization levels ($<$ 1 per cent) for both explosion models. Polarization in the continuum peaks at $\sim$ 0.1$-$0.3 per cent and decreases after maximum light, in excellent agreement with spectropolarimetric data of normal Type Ia supernovae. Higher degrees of polarization are found across individual spectral lines. In particular, the synthetic Si ii {\lambda}6355 profiles are polarized at levels that match remarkably well the values observed in normal Type Ia supernovae, while the low degrees of polarization predicted across the O i {\lambda}7774 region are consistent with the non-detection of this feature in current data. We conclude that our models can reproduce many of the characteristics of both flux and polarization spectra for well-studied Type Ia supernovae, such as SN 2001el and SN 2012fr. However, the two models considered here cannot account for the unusually high level of polarization observed in extreme cases such as SN 2004dt.

The violent merger of two carbon-oxygen white dwarfs has been proposed as a viable progenitor for some Type Ia supernovae. However, it has been argued that the strong ejecta asymmetries produced by this model might be inconsistent with the low degree of polarisation typically observed in Type Ia supernova explosions. Here, we test this claim by carrying out a spectropolarimetric analysis for the model proposed by Pakmor et al. (2012) for an explosion triggered during the merger of a 1.1 M$_{\odot}$ and 0.9 M$_{\odot}$ carbon-oxygen white dwarf binary system. Owing to the asymmetries of the ejecta, the polarisation signal varies significantly with viewing angle. We find that polarisation levels for observers in the equatorial plane are modest ($\lesssim$ 1 per cent) and show clear evidence for a dominant axis, as a consequence of the ejecta symmetry about the orbital plane. In contrast, orientations out of the plane are associated with higher degrees of polarisation and departures from a dominant axis. While the particular model studied here gives a good match to highly-polarised events such as SN 2004dt, it has difficulties in reproducing the low polarisation levels commonly observed in normal Type Ia supernovae. Specifically, we find that significant asymmetries in the element distribution result in a wealth of strong polarisation features that are not observed in the majority of currently available spectropolarimetric data of Type Ia supernovae. Future studies will map out the parameter space of the merger scenario to investigate if alternative models can provide better agreement with observations.

The Public European Southern Observatory Spectroscopic Survey of Transient Objects (PESSTO) began as a public spectroscopic survey in April 2012. We describe the data reduction strategy and data products which are publicly available through the ESO archive as the Spectroscopic Survey Data Release 1 (SSDR1). PESSTO uses the New Technology Telescope with EFOSC2 and SOFI to provide optical and NIR spectroscopy and imaging. We target supernovae and optical transients brighter than 20.5mag for classification. Science targets are then selected for follow-up based on the PESSTO science goal of extending knowledge of the extremes of the supernova population. The EFOSC2 spectra cover 3345-9995A (at resolutions of 13-18 Angs) and SOFI spectra cover 0.935-2.53 micron (resolutions 23-33 Angs) along with JHK imaging. This data release contains spectra from the first year (April 2012 - 2013), consisting of all 814 EFOSC2 spectra and 95 SOFI spectra (covering 298 distinct objects), in standard ESO Phase 3 format. We estimate the accuracy of the absolute flux calibrations for EFOSC2 to be typically 15%, and the relative flux calibration accuracy to be about 5%. The PESSTO standard NIR reduction process does not yet produce high accuracy absolute spectrophotometry but the SOFI JHK imaging will improve this. Future data releases will focus on improving the automated flux calibration of the data products.

We present a Monte Carlo radiative transfer technique for calculating synthetic spectropolarimetry for multi-dimensional supernova explosion models. The approach utilises "virtual-packets" that are generated during the propagation of the Monte Carlo quanta and used to compute synthetic observables for specific observer orientations. Compared to extracting synthetic observables by direct binning of emergent Monte Carlo quanta, this virtual-packet approach leads to a substantial reduction in the Monte Carlo noise. This is vital for calculating synthetic spectropolarimetry (since the degree of polarisation is typically very small) but also useful for calculations of light curves and spectra. We first validate our approach via application of an idealised test code to simple geometries. We then describe its implementation in the Monte Carlo radiative transfer code ARTIS and present test calculations for simple models for Type Ia supernovae. Specifically, we use the well-known one-dimensional W7 model to verify that our scheme can accurately recover zero polarisation from a spherical model, and to demonstrate the reduction in Monte Carlo noise compared to a simple packet-binning approach. To investigate the impact of aspherical ejecta on the polarisation spectra, we then use ARTIS to calculate synthetic observables for prolate and oblate ellipsoidal models with Type Ia supernova compositions.