ABSTRACT
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\, \mathrm{M}_\odot$ and half-opening angle of the lanthanide-rich component Φ = 30° 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 multidimensional models in place of combining one-dimensional models to infer important parameters, such as the ejecta masses. We finally explore the impact of Mej and Φ 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.
ABSTRACT
Neutron star mergers eject neutron-rich matter in which heavy elements are synthesized. The decay of these freshly synthesized elements powers electromagnetic transients (‘macronovae’ or ...‘kilonovae’) whose luminosity and colour strongly depend on their nuclear composition. If the ejecta are very neutron-rich (electron fraction Ye < 0.25), they contain fair amounts of lanthanides and actinides that have large opacities and therefore efficiently trap the radiation inside the ejecta so that the emission peaks in the red part of the spectrum. Even small amounts of this high-opacity material can obscure emission from lower lying material and therefore act as a ‘lanthanide curtain’. Here, we investigate how a relativistic jet that punches through the ejecta can potentially push away a significant fraction of the high opacity material before the macronova begins to shine. We use the results of detailed neutrino-driven wind studies as initial conditions and explore with 3D special relativistic hydrodynamic simulations how jets are propagating through these winds. Subsequently, we perform Monte Carlo radiative transfer calculations to explore the resulting macronova emission. We find that the hole punched by the jet makes the macronova brighter and bluer for on-axis observers during the first few days of emission, and that more powerful jets have larger impacts on the macronova.
Context. Both multi-messenger astronomy and new high-throughput wide-field surveys require flexible tools for the selection and analysis of astrophysical transients. Aims. Here we introduce the alert ...management, photometry, and evaluation of light curves (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; it 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. Methods. 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). Results. 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. Conclusions. 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 is made possible by a novel way of capturing transient information through sequences of evolving states, and interfaces that allow new code to be natively applied to a full stream of alerts. This text also introduces a method by which 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.
ABSTRACT We present the first spectropolarimetric observations of a hydrogen-free superluminous supernova (SLSN) at z = 0.1136, namely SN 2015bn. The transient shows significant polarization at both ...of the observed epochs: one 24 days before maximum light in the rest-frame, and the other 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 axisymmetric ellipsoidal configuration for the ejecta. We find that the wavelength dependence of the polarization is possibly due to a strong wavelength dependence in the line opacity, while the higher level of polarization 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 SLSN is similar to that of stripped-envelope core-collapse SNe connected to GRB, while the overall evolution of the ejecta shape could be consistent with a central engine.
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 rZTF-band ( ), as expected in the expanding fireball model. There are, however, individual SNe that are clearly inconsistent with . We estimate a mean rise time of 18.9 days (with a range extending from ∼15 to 22 days), 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 = 2, in which case we estimate a mean rise time of 21.7 days (with a range from ∼18 to 23 days). 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.
ABSTRACT
We predict linear polarization for a radioactively powered kilonova following the merger of a black hole and a neutron star. Specifically, we perform 3D 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 disc-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 disc wind shows no polarization in the optical, while a signal increasing at longer wavelengths and reaching $\sim 1\!-\!6{{\ \rm per\ cent}}$ at $2\, \mu$m depending on the orientation. The model with a disc-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 disc-wind component. The detection of optical polarization would unambiguously reveal the presence of a lanthanide-free disc-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 ∼48 h 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.
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 x 10(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 (less than or similar to 1 pc) are located close to the centre 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-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.
There is a strong degeneracy between the luminosity distance (DL) and the observer viewing angle (θ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 from broadband 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 GRB 170817A. Combining the constraints on θobs with the GW data, we find an improvement of 24% on H0. The observer angle constraints are insensitive to other model parameters, e.g., the ejecta mass, the 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 θobs. While the improvement on H0 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 (DL ∼ 100 Mpc). Our results are insensitive to the assumption of the peculiar velocity of the kN host galaxy.
ABSTRACT
The detection of the binary neutron star (BNS) merger GW170817 and the associated electromagnetic (EM) counterpart, the ‘kilonova’ (kN) AT2017gfo, opened a new era in multimessenger ...astronomy. However, despite many efforts, it has been proven very difficult to find additional kNe, even though LIGO/Virgo has reported at least one BNS event during their latest run, O3. The focus of this work is the exploration of the sensitivity of the adopted optical surveys searching for kNe during O3. We propose ways to optimize the choices of filters and survey depth to boost the detection efficiency for these faint and fast-evolving transients in the future. In particular, we use kN models to explore the dependence on ejecta mass, geometry, viewing angle, wavelength coverage, and source distance. We find that the kN detection efficiency has a strong viewing-angle dependence, especially for filters blueward of i-band. This loss of sensitivity can be mitigated by early, deep, observations. Efficient gri counterpart searches for kNe at ∼200 Mpc would require reaching a limiting magnitude mlim = 23 mag, to ensure good sensitivity over a wide range of the model phase-space. We conclude that kN searches during O3 were generally too shallow to detect BNS optical counterparts, even under optimistic assumptions.
ABSTRACT
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 behaviour 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.