Search

At the dawn of the gamma-ray burst (GRB) afterglow era, a Cepheid-like correlation was discovered between time variability V and isotropic-equivalent peak luminosity Liso of the prompt emission of about a dozen long GRBs with measured redshift available at that time. Soon afterwards, the correlation was confirmed against a sample of about 30 GRBs, despite being affected by significant scatter. Unlike the minimum variability timescale (MVT), V measures the relative power of short-to-intermediate timescales. We aim to test the correlation using about two hundred long GRBs with spectroscopically measured redshift, detected by Swift, Fermi, and Konus/WIND, for which both observables can be accurately estimated. For all the selected GRBs, variability was calculated according to the original definition using the 64-ms background-subtracted light curves of Swift/BAT (Fermi/GBM) in the 15-150 (8-900) keV energy passband. Peak luminosities were either taken from literature or derived from modelling broad-band spectra acquired with either Konus/WIND or Fermi/GBM. The statistical significance of the correlation has weakened to <~2%, mostly due to the appearance of a number of smooth and luminous GRBs characterised by a relatively small V. At odds with most long GRBs, 3 out of 4 long-duration merger candidates have high V and low Liso. Luminosity is more tightly connected with shortest timescales measured by MVT rather than short-to-intermediate ones, measured by V. We discuss the implications on internal dissipation models and the role of the e+- photosphere. We identified a few, smooth GRBs with a single broad pulse and low V, that might have an external shock origin, in contrast with most GRBs. The combination of high variability (V>~0.1), low luminosity (Liso<~10^51 erg s^-1) and short MVT (<~ 0.1 s) could be a good indicator for a compact binary merger origin., Comment: 14 pages, 8 figures, accepted by A&A (corrected small typo in the Acknowledgments)

Context. The dissipation process responsible for the long gamma-ray burst (GRB) prompt emission and the kind of dynamics that drives the release of energy as a function of time are still key open issues. We recently found that the distribution of the number of peaks per GRB is described by a mixture of two exponentials, suggesting the existence of two behaviours that turn up as peak-rich and peak-poor time profiles. Aims. Our aims are to study the distribution of the number of peaks per GRB of the entire catalogue of about 3000 GRBs observed by the Fermi Gamma-ray Burst Monitor (GBM) and to make a comparison with previous results obtained from other catalogues. Methods. We identified GRB peaks using the MEPSA code and modelled the resulting distribution following the same procedure that was adopted in the previous analogous investigation. Results. We confirm that only a mixture of two exponentials can model the distribution satisfactorily, with model parameters that fully agree with those found from previous analyses. In particular, we confirm that (21 +- 4)% of the observed GRBs are peak-rich (8 +- 1 peaks per GRB on average), while the remaining 80% are peak-poor (2.12 +- 0.10 peaks per GRB on average). Conclusions. We confirm the existence of two different components, peak-poor and peak-rich GRBs, that make up the observed GRB populations. Together with previous analogous results from other GRB catalogues, these results provide compelling evidence that GRB prompt emission is governed by two distinct regimes., Comment: 6 pages, 3 figures, accepted as a letter by A&A

Recent observations of core-collapse supernovae (SNe) suggest aspherical explosions. Globally aspherical structures in SN explosions are regarded as the key for understanding their explosion mechanism. However, the exact explosion geometries from the inner cores to the outer envelopes are poorly understood. Here, we present photometric, spectroscopic and polarimetric observations of the Type IIP SN 2021yja and discuss its explosion geometry, in comparison to those of other Type IIP SNe that show large-scale aspherical structures in their hydrogen envelopes (SNe 2012aw, 2013ej and 2017gmr). During the plateau phase, SNe 2012aw and 2021yja exhibit high continuum polarization characterized by two components with perpendicular polarization angles. This behavior can be interpreted to be due to a bipolar explosion, composed of a polar (energetic) and an equatorial (bulk) components of the SN ejecta. In such a bipolar explosion, an aspherical axis created by the polar ejecta would be dominating at early phases, while the perpendicular axis along the equatorial ejecta would emerge at late phases after the receding of the photosphere in the polar ejecta. The interpretation of the bipolar explosions in SNe 2012aw and 2021yja is also supported by other observational properties, including the time evolution of the line velocities and the line shapes in the nebular spectra. The polarization of other Type IIP SNe that show large-scale aspherical structures in the hydrogen envelope (SNe 2013ej and 2017gmr) is also consistent with the bipolar-explosion scenario, although this is not conclusive., Comment: 16 pages, 10 figures, accepted for publication in A&A letter

The variety of long duration gamma-ray burst (LGRB) light curves (LCs) encode a wealth of information on how LGRB engines release energy following the collapse of the progenitor star. Attempts to characterise GRB LCs focused on a number of properties, such as the minimum variability timescale, power density spectra (both ensemble average and individual), or with different definitions of variability. In parallel, a characterisation as a stochastic process was pursued by studying the distributions of waiting times, peak flux, fluence of individual peaks within GRB time profiles. Yet, the question remains as to whether the diversity of profiles can be described in terms of a common stochastic process. Here we address this issue by studying for the first time the distribution of the number of peaks in a GRB profile. We used four different GRB catalogues: CGRO/BATSE, Swift/BAT, BeppoSAX/GRBM, and Insight-HXMT. The statistically significant peaks were identified by means of well tested algorithm MEPSA and further selected by applying a set of thresholds on signal-to-noise ratio. We then extracted the corresponding distributions of number of peaks per GRB. Among the different models considered (power-law, simple or stretched exponential) only a mixture of two exponentials models all the observed distributions, suggesting the existence of two distinct behaviours: (i) an average number of 2.1+-0.1 peaks per GRB ("peak poor") and accounting for about 80% of the observed population of GRBs; (ii) an average number of 8.3+-1.0 peaks per GRB ("peak rich") and accounting for the remaining 20% of the observed population. We associate the class of peak-rich GRBs with the presence of sub-second variability, which seems to be absent among peak-poor GRBs. The two classes could result from two different regimes through which GRB engines release energy or through which energy is dissipated into gamma-rays., Comment: 7 pages, 5 figures, accepted by A&A

Discovered more than 50 years ago, gamma-ray burst (GRB) prompt emission remains the most puzzling aspect of GRB physics. Its complex and irregular nature should reveal how newborn GRB engines release their energy. In this respect, the possibility that GRB engines could operate as self-organized critical (SOC) systems has been put forward. Here, we present the energy, luminosity, waiting time, and duration distributions of individual pulses of GRBs with known redshift detected by the Fermi Gamma-ray Burst Monitor (GBM). This is the first study of this kind in which selection effects are accounted for. The compatibility of our results with the framework of SOC theory is discussed. We found evidence for an intrinsic break in the power-law models that describe the energy and the luminosity distributions., Comment: Accepted by ApJ, 24 pages, 14 figures, 5 tables

GRB~230812B is a bright and relatively nearby ($z =0.36$) long gamma-ray burst (GRB) that has generated significant interest in the community and has thus been observed over the entire electromagnetic spectrum. We report over 80 observations in X-ray, ultraviolet, optical, infrared, and sub-millimeter bands from the GRANDMA (Global Rapid Advanced Network for Multi-messenger Addicts) network of observatories and from observational partners. Adding complementary data from the literature, we then derive essential physical parameters associated with the ejecta and external properties (i.e. the geometry and environment) of the GRB and compare with other analyses of this event. We spectroscopically confirm the presence of an associated supernova, SN2023pel, and we derive a photospheric expansion velocity of v $\sim$ 17$\times10^3$ km s$^{-1}$. We analyze the photometric data first using empirical fits of the flux and then with full Bayesian Inference. We again strongly establish the presence of a supernova in the data, with a maximum (pseudo-)bolometric luminosity of $5.75 \times 10^{42}$ erg/s, at $15.76^{+0.81}_{-1.21}$ days (in the observer frame) after the trigger, with a half-max time width of 22.0 days. We compare these values with those of SN1998bw, SN2006aj, and SN2013dx. Our best-fit model favours a very low density environment ($\log_{10}({n_{\rm ISM}/{\rm cm}^{-3}}) = -2.38^{+1.45}_{-1.60}$) and small values for the jet's core angle $\theta_{\rm core} = 1.54^{+1.02}_{-0.81} \ \rm{deg}$ and viewing angle $\theta_{\rm obs} = 0.76^{+1.29}_{-0.76} \ \rm{deg}$. GRB 230812B is thus one of the best observed afterglows with a distinctive supernova bump.

The explosion processes of supernovae (SNe) are imprinted in their explosion geometries. Here, we study the intrinsic polarization of 15 hydrogen-rich core-collapse SNe and explore the relation with the photometric and spectroscopic properties. Our sample shows diverse properties of the continuum polarization. The polarization of most SNe has a low degree at early phases but shows a sudden rise to $\sim 1$ \% degree at certain points during the photospheric phase as well as a slow decline during the tail phase, with a constant polarization angle. The variation in the timing of peak polarisation values implies diversity in the explosion geometry: some SNe have aspherical structures only in their helium cores, while in other SNe these reach out to a significant part of the outer hydrogen envelope with a common axis from the helium core to the hydrogen envelope. Other SNe show high polarization from early phases and a change of the polarization angle around the middle of the photospheric phase. This implies that the ejecta are significantly aspherical to the outermost layer and have multi-directional aspherical structures. Exceptionally, the Type~IIL SN~2017ahn shows low polarization at both the photospheric and tail phases. Our results show that the timing of the polarization rise in Type~IIP SNe is likely correlated with their brightness, velocity and the amount of radioactive Ni produced: brighter SNe with faster ejecta velocity and a larger $^{56}$Ni mass have more extended-aspherical explosion geometries. In particular, there is a clear correlation between the timing of the polarization rise and the explosion energy, that is, the explosion asphericity is proportional to the explosion energy. This implies that the development of a global aspherical structure, e.g., a jet, might be the key to realising an energetic SN in the mechanism of SN explosions., Comment: 27 pages, 30 figures

We consider the optimisation of the observing strategy (cadence, exposure time and filter choice) using medium size (2-m class) optical telescopes in the follow-up of kilonovae localised with arcminute accuracy to be able to distinguish among various kilonova models and viewing angles. To develop an efficient observation plan, we made use of the synthetic light curves obtained with the Monte Carlo radiative transfer code POSSIS for different kilonova models and as a function of different viewing angles and distances. By adding the appropriate photon counting noise to the synthetic light curves, we analysed four alternative sequences having the same total time exposure of 8 hours, with different time windows (0.5, 1, 2, 4 h), each with $i$, $r$, and $u$ filters, to determine the observing sequence that maximises the chance of a correct identification of the model parameters. We suggest to avoid $u$ filter and to avoid the use of colour curves. We also found that, if the error on distance is $\le$ 2%, 0.5, 1, 2-hour time window sequences are equivalent, so we suggest to use 2-hour one, because it has 1 day cadence, so it can be easily realised. When the distance of the source is unknown, 0.5 h time window sequence is preferable., Comment: 9 pages, 8 figures, published in MNRAS

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., Comment: Submitted. Comments welcome! Nature (2023)

This is the first in a collection of three papers introducing the science with an ultra-violet (UV) space telescope on an approximately 130~kg small satellite with a moderately fast re-pointing capability and a real-time alert communication system approved for a Czech national space mission. The mission, called Quick Ultra-Violet Kilonova surveyor - QUVIK, will provide key follow-up capabilities to increase the discovery potential of gravitational wave observatories and future wide-field multi-wavelength surveys. The primary objective of the mission is the measurement of the UV brightness evolution of kilonovae, resulting from mergers of neutron stars, to distinguish between different explosion scenarios. The mission, which is designed to be complementary to the Ultraviolet Transient Astronomy Satellite - ULTRASAT, will also provide unique follow-up capabilities for other transients both in the near- and far-UV bands. Between the observations of transients, the satellite will target other objects described in this collection of papers, which demonstrates that a small and relatively affordable dedicated UV-space telescope can be transformative for many fields of astrophysics., Comment: Accepted for publication in Space Science Reviews

GRB 221009A is the brightest Gamma-Ray Burst (GRB) detected in more than 50 years of study. In this paper, we present observations in the X-ray and optical domains after the GRB obtained by the GRANDMA Collaboration (which includes observations from more than 30 professional and amateur telescopes) and the Insight-HXMT Collaboration. We study the optical afterglow with empirical fitting from GRANDMA+HXMT data, augmented with data from the literature up to 60 days. We then model numerically, using a Bayesian approach, the GRANDMA and HXMT-LE afterglow observations, that we augment with Swift-XRT and additional optical/NIR observations reported in the literature. We find that the GRB afterglow, extinguished by a large dust column, is most likely behind a combination of a large Milky-Way dust column combined with moderate low-metallicity dust in the host galaxy. Using the GRANDMA+HXMT-LE+XRT dataset, we find that the simplest model, where the observed afterglow is produced by synchrotron radiation at the forward external shock during the deceleration of a top-hat relativistic jet by a uniform medium, fits the multi-wavelength observations only moderately well, with a tension between the observed temporal and spectral evolution. This tension is confirmed when using the extended dataset. We find that the consideration of a jet structure (Gaussian or power-law), the inclusion of synchrotron self-Compton emission, or the presence of an underlying supernova do not improve the predictions, showing that the modelling of GRB22109A will require going beyond the most standard GRB afterglow model. Placed in the global context of GRB optical afterglows, we find the afterglow of GRB 221009A is luminous but not extraordinarily so, highlighting that some aspects of this GRB do not deviate from the global known sample despite its extreme energetics and the peculiar afterglow evolution., Comment: Accepted to ApJL for the special issue, 37 pages, 23 pages main text, 6 tables, 13 figures

We present linear polarimetry for seven hydrogen-poor superluminous supernovae (SLSNe-I). For SN 2017gci, for which we present two epochs of spectropolarimetry at +3 d and +29 d post-peak in rest frame, accompanied by four epochs of imaging polarimetry up to +108 d. The spectropolarimetry at +3 d shows increasing polarisation degree P towards the redder wavelengths and exhibits signs of axial symmetry, but at +29 d P=0 throughout the spectrum implying that the photosphere of SN 2017gci evolved from a slightly aspherical configuration to a more spherical one in the first month post-peak. However, an increase of P to 0.5% at +55 d accompanied by a different orientation of the axial symmetry compared to +3 d implies the presence of additional sources of polarisation at this phase. The increase in polarisation is possibly caused by interaction with circumstellar matter as already suggested by a knee in the light curve and a possible detection of broad Ha emission. We also analysed the sample of all 16 SLSNe-I with polarimetry to date. The data taken during the early spectroscopic phase show consistently low P indicating spherical photospheres. No clear relation between the polarimetry and spectral phase was seen when the spectra resemble Type Ic SNe during the photospheric and nebular phases. The light curve decline rate also shows no clear relation with the polarisation properties. While only slow-evolving SLSNe-I have shown non-zero P, the fast-evolving ones have not been observed at sufficiently late times to conclude that none of them exhibit changing P. However, the four SLSNe-I with increasing polarisation degree also have irregular light curve declines. For up to half of them, the photometric, spectroscopic and polarimetric properties are affected by CSM interaction. As such CSM interaction clearly plays an important role in understanding the polarimetric evolution of SLSNe-I., Comment: Accepted to A&A on 18/03/2023

The merger of two compact objects of which at least one is a neutron star is signalled by transient electromagnetic emission in a kilonova (KN). This event is accompanied by gravitational waves and possibly other radiation messengers such as neutrinos or cosmic rays. The electromagnetic emission arises from the radioactive decay of heavy $r-$process elements synthesized in the material ejected during and after the merger. In this paper we show that the analysis of KNe light curves can provide cosmological distance measurements and constrain the properties of the ejecta. In this respect, MAAT, the new Integral Field Unit in the OSIRIS spectrograph on the $10.4$ m Gran Telescopio CANARIAS (GTC), is well suited for the study of KNe by performing absolute spectro-photometry over the entire 3600-10000 Angstron spectral range. Here, we study the most representative cases regarding the scientific interest of KNe from binary neutron stars, and we evaluate the observational prospects and performance of MAAT on the GTC to do the following: a) study the impact of the equation of state on the KN light curve, and determine to what extent bounds on neutron star (NS) radii or compactness deriving from KN peak magnitudes can be identified and b) measure the Hubble constant, $H_0$, with precision improved by up to 40$\%$, when both gravitational wave data and photometric-light curves are used. In this context we discuss how the equation of state, the viewing angle, and the distance affect the precision and estimated value of $H_0$., Comment: 21 pages, 23 figures, matches published version in A&A 666, A67 (2022)

We present a spectroscopic analysis of Type I superluminous supernova (SLSN-I), SN 2018bsz. While it closely resembles SLSNe-I, the multi-component H$\alpha$ line appearing at $\sim30$ d post-maximum is the most atypical. The H$\alpha$ is characterised by two emission components, one at $+3000$ km/s and a second at $-7500$ km/s, with a third, near-zero velocity component appearing after a delay. The blue and central components can be described by Gaussian profiles of intermediate width, but the red component is significantly broader and Lorentzian. The blue component evolves towards lower velocity before fading at $100$ d post-peak, concurrently with a light curve break. Multi-component profiles are observed in other hydrogen lines including Pa$\beta$, and in lines of Ca II and He I. Spectropolarimetry obtained before (10.2 d) and after (38.4 d) the appearance of the H lines show a large shift on the Stokes $Q$ -- $U$ plane consistent with SN 2018bsz undergoing radical changes in its geometry. Assuming the SN is almost unpolarised at 10.2 d, the continuum polarisation at 38.4 d reaches $P \sim1.8\%$ implying a highly asymmetric configuration. We propose that the observed evolution of SN 2018bsz can be explained by highly aspherical CSM. After the SN explosion, the CSM is quickly overtaken by the ejecta, but as the photosphere starts to recede, the different CSM regions re-emerge producing the peculiar line profiles. Based on the first appearance of H$\alpha$, we can constrain the distance of the CSM to be less than $430$ AU, or even lower ($<87$ AU) if the pre-peak plateau is related to an eruption that created the CSM. The presence of CSM has been inferred for other SLSNe-I. However, it is not clear whether the rare properties of SN 2018bsz can be generalised for SLSNe-I or whether they are the result of an uncommon evolutionary path, possibly involving a binary companion., Comment: Accepted to A&A on 22/06/2022

Interpreting high-energy, astrophysical phenomena, such as supernova explosions or neutron-star collisions, requires a robust understanding of matter at supranuclear densities. However, our knowledge about dense matter explored in the cores of neutron stars remains limited. Fortunately, dense matter is not only probed in astrophysical observations, but also in terrestrial heavy-ion collision experiments. In this work, we use Bayesian inference to combine data from astrophysical multi-messenger observations of neutron stars and from heavy-ion collisions of gold nuclei at relativistic energies with microscopic nuclear theory calculations to improve our understanding of dense matter. We find that the inclusion of heavy-ion collision data indicates an increase in the pressure in dense matter relative to previous analyses, shifting neutron-star radii towards larger values, consistent with recent NICER observations. Our findings show that constraints from heavy-ion collision experiments show a remarkable consistency with multi-messenger observations and provide complementary information on nuclear matter at intermediate densities. This work combines nuclear theory, nuclear experiment, and astrophysical observations, and shows how joint analyses can shed light on the properties of neutron-rich supranuclear matter over the density range probed in neutron stars., Comment: 7 pages, 2 figures, Supplemental Material, minor changes

The origin of the diverse light-curve shapes of Type II supernovae (SNe), and whether they come from similar or distinct progenitors, has been actively discussed for decades. Here we report spectropolarimetry of two fast declining Type II (Type IIL) SNe: SN 2013ej and SN 2017ahn. SN 2013ej exhibited high continuum polarization from very soon after the explosion to the radioactive tail phase with time-variable polarization angles. The origin of this polarimetric behavior can be interpreted as the combination of two different aspherical structures, namely an aspherical interaction of the SN ejecta with circumstellar matter (CSM) and an inherently aspherical explosion. Aspherical explosions are a common feature of slowly declining Type II (Type IIP) SNe. By contrast, SN 2017ahn showed low polarization not only in the photospheric phase but also in the radioactive tail phase. This low polarization in the tail phase, which has never before been observed in other Type IIP/L SNe, suggests that the explosion of SN 2017ahn was nearly spherical. These observations imply that Type IIL SNe have, at least, two different origins: they result from stars that have different explosion properties and/or different mass-loss processes. This fact might indicate that 13ej-like Type IIL SNe originate from a similar progenitor to those of Type IIP SNe accompanied by an aspherical CSM interaction, while 17ahn-like Type IIL SNe come from a more massive progenitor with less hydrogen in its envelope., Comment: 19 pages, 14 figures; accepted for publication in MNRAS

The detection of AT2017gfo proved that binary neutron star mergers are progenitors of kilonovae. Using a combination of numerical-relativity and radiative-transfer simulations, the community has developed sophisticated models for these transients for a wide portion of the expected parameter space. Using these simulations and surrogate models made from them, it has been possible to perform Bayesian inference of the observed signals to infer properties of the ejected matter. It has been pointed out that combining inclination constraints derived from the kilonova with gravitational-wave measurements increases the accuracy with which binary parameters can be measured and allows a more accurate inference of the Hubble Constant. In order to not introduce biases, constraints on the inclination angle for AT2017gfo should be insensitive to the employed models. In this work, we compare different assumptions about the ejecta and radiative reprocesses used by the community and we investigate their impact on the parameter inference. While most inferred parameters agree, we find disagreement between posteriors for the inclination angle for different geometries that have been used in the literature. According to our study, the inclusion of reprocessing of the photons between different ejecta types improves the modeling fits to AT2017gfo and in some cases affects the inferred constraints. Our study motivates the inclusion of large $\sim$ 1 mag uncertainties in the kilonova models employed for Bayesian analysis to capture yet unknown systematics, especially when inferring inclination angles, although smaller uncertainties seem appropriate to capture model systematics for other parameters. We also use this method to impose soft constraints on the ejecta geometry of the kilonova AT2017gfo.

A non-universal $U(1)_{X}$ extension to the Standard Model composed of two scalar doublets and two scalar singlets together with three additional quark singlets and two lepton singlets and three generations of right-handed and Majorana neutrinos is made to explain lepton mass hierarchy, neutrino masses via inverse seesaw mechanism and muon anomalous magnetic moment in an anomaly free framework. In the present model, exotic and Standard Model particles acquire mass thanks to vacuum expectation values at different scales, yet the electron and the lightest neutrino are tree level massless but massive at one-loop level. By considering a numerical exploration and under the constraint of the Higgs mass, neutrino mass differences and PMNS matrix, it is found that only contributions due to exotic neutrinos interacting with charged scalars are relevant to muon $g-2$, though they are negative. Thus, the SUSY extension is considered and it is found that muon $g-2$ can be explained by allowing $U(1)_{X}$ vacuum expectation values to lie in the TeV scale thanks to SUSY soft-breaking interactions for at least $\sim 10^{5}$ GeV masses. Thus, the contribution due to exotic neutrinos interacting with $W$ gauge bosons is positive and no longer negligible which added to all other contributions might explain the anomaly., Comment: 42 pages, 26 figures, 3 tables

We predict linear polarization for a radioactively-powered kilonova following the merger of a black hole and a neutron star. Specifically, we perform 3-D Monte Carlo radiative transfer simulations for two different models, both featuring a lanthanide-rich dynamical ejecta component from numerical-relativity simulations while only one including an additional lanthanide-free disk wind component. We calculate polarization spectra for nine different orientations at 1.5, 2.5 and 3.5 d after the merger and in the $0.1-2\,\mu$m wavelength range. We find that both models are polarized at a detectable level 1.5 d after the merger while show negligible levels thereafter. The polarization spectra of the two models are significantly different. The model lacking a disk wind shows no polarization in the optical, while a signal increasing at longer wavelengths and reaching $\sim1\%-6\%$ at $2\,\mu$m depending on the orientation. The model with a disk-wind component, instead, features a characteristic "double-peak" polarization spectrum with one peak in the optical and the other in the infrared. Polarimetric observations of future events will shed light on the debated neutron richness of the disk-wind component. The detection of optical polarization would unambiguously reveal the presence of a lanthanide-free disk-wind component, while polarization increasing from zero in the optical to a peak in the infrared would suggest a lanthanide-rich composition for the whole ejecta. Future polarimetric campaigns should prioritize observations in the first $\sim48$ hours and in the $0.5-2\,\mu$m range, where polarization is strongest, but also explore shorter wavelengths/later times where no signal is expected from the kilonova and the interstellar polarization can be safely estimated., Comment: 8 pages, 4 figures; accepted for publication in MNRAS

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., Comment: 27 pages, 14 figures, accepted in ApJ

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., Comment: 11 pages, 8 figures, 2 tables, Published in MNRAS

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., Comment: 7 pages, 3 figures; accepted for publication in A&A; metadata update

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., Comment: 52 pages, revised version now accepted for publication in A&A. Abstract abridged to meet arXiv requirements

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 large, unique data set of 127 SNe$\,$Ia with exquisite coverage by the Zwicky Transient Facility (ZTF). High-cadence (six 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 that 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.9$\,$d (with a range extending from $\sim$15 to 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.7$\,$d (with a range from $\sim$18 to 23$\,$d). The sample includes a handful of rare and peculiar SNe$\,$Ia. Finally, we conclude with a discussion of lessons learned from the ZTF sample that can eventually be applied to observations from the Vera C. Rubin Observatory., Comment: 31 pages, 17 figures; ApJ accepted

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., Comment: 15 pages, 12 figures. Accepted for publication by MNRAS

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., Comment: 7 pages; 2 Figures, 2 Tables, accepted to ApJ

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., Comment: 6 pages, 3 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society Letters

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., Comment: Submitted to ApJ, 25 pages, plus Appendix

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., Comment: Updated to match version accepted for publication in A&A

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

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., Comment: published in Nature Astronomy, www.nature.com/articles/s41550-018-0593-y/

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., Comment: 13 pages, 8 figures, accepted for publication in MNRAS

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., Comment: Accepted for publication in A&A after minor corrections to first arXiv version

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., Comment: 14 pages, 6 figures, 5 tables; accepted for publication in MNRAS

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., Comment: 14 pages, 17 figures. Accepted to PASP

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., Comment: Nature, in press, DOI 10.1038/nature24303. Data files will be made available at http://www.pessto.org