Search

With GW170817 being the only multimessenger gravitational wave (GW) event with an associated kilonova (KN) detected so far, there exists a pressing need for realistic estimation of the GW localization uncertainties and rates, as well as optimization of available telescope time to enable the detection of new KNe. For this purpose, we simulate GW events assuming a data-driven, GW-motivated distribution of binary parameters for the LIGO/Virgo/KAGRA (LVK) fourth and fifth observing runs (O4 and O5). For our particular case of estimating KN detection rates in the O4 and O5 runs, we map the binary neutron star (BNS) and neutron star-black hole (NSBH) properties to the optical light curves arising from r-process nucleosynthesis in the ejecta. We use the simulated population of KNe to generate follow-up observing plans, with the primary goal of optimizing detection with the Gravitational Wave Multi-Messenger Astronomy DECam Survey (GW-MMADS). We explore KN detectability as a function of mass, distance, and spin of the binaries. Out of the mergers that produce a KN in our simulations, we expect detectable KNe for DECam-like instruments at a per-year rate of: $1-16$ ($0-1$) for BNS (NSBH) in O4, and $16-165$ ($1-12$) for BNS (NSBH) in O5, using a fiducial exposure time and conditional on the uncertainty on the equation of state (EOS) and volumetric rates of the mergers. Taking into account observability constraints, our scheduler covers the location of the KN $\sim 30-38\%$ of the times for our fiducial EOS. We provide the depths needed to detect a significant fraction of our simulated mergers for the astronomical community to use in their follow-up campaigns., Comment: 18 pages, 10 figures, 10 tables

Kilonovae, possible electromagnetic counterparts to neutron star mergers, provide important information about high-energy transient phenomena and, in principle, also allow us to obtain information about the source properties responsible for powering the kilonova. Unfortunately, numerous uncertainties exist in kilonova modeling that, at the current stage, hinder accurate predictions. Hence, one has to account for possible systematic modeling uncertainties when interpreting the observed transients. In this work, we provide a data-driven approach to account for time-dependent and frequency-dependent uncertainties in kilonova models. Through a suite of tests, we find that the most reliable recovery of the source parameters and description of the observational data can be obtained through a combination of kilonova models with time- and frequency-dependent systematic uncertainties. We apply our new method to analyze AT2017gfo. While recovering a total ejecta mass consistent with previous studies, our approach gives insights into the temporal and spectral evolution of the systematic uncertainties of this kilonova. We consistently find a systematic error below $1$ mag between $1$ to $5$ days after the merger. Our work addresses the need for early follow-up of kilonovae at earlier times, and improved modeling of the kilonova at later times, to reduce the uncertainties outside of this time window.

The detection of the gravitational-wave event GW230529, presumably a neutron star-black hole (NSBH) merger, by the LIGO-Virgo-KAGRA (LVK) Collaboration is an exciting discovery for multimessenger astronomy. The black hole (BH) has a high probability of falling within the ''mass gap'' between the peaks of the neutron star (NS) and the BH mass distributions. Because of the low primary mass, the binary is more likely to produce an electromagnetic counterpart than previously detected NSBH mergers. We investigate the possible kilonova (KN) emission from GW230529, and find that if it was an NSBH, there is a $\sim$ 2-41% probability (depending on the assumed equation of state) that GW230925 produced a KN with magnitude peaking at $\sim 1-2$ day post merger at $g \lesssim 23.5$, $i<23$. Hence, it could have been detected by ground-based telescopes. If it was a binary neutron star (BNS) merger, we find $\sim$ 0-12% probability that it produced a KN. Motivated by these numbers, we simulated a broader population of mgNSBH mergers that may be detected in O4, and we obtained a 9-21% chance of producing a KN, which would be detectable with $g\lesssim 25$ and $ i \lesssim 24$, typically fainter than what is expected from GW230529. Based on these findings, DECam-like instruments may be able to detect up to 80% of future mgNSBH KNe, thus up to $\sim1$ multimessenger mgNSBH per year may be discoverable at the current level of sensitivity (O4)., Comment: 13 pages, 4 figures, 2 tables

The cores of stars are the cosmic furnaces where light elements are fused into heavier nuclei. The fusion of hydrogen to helium initially powers all stars. The ashes of the fusion reactions are then predicted to serve as fuel in a series of stages, eventually transforming massive stars into a structure of concentric shells. These are composed of natal hydrogen on the outside, and consecutively heavier compositions inside, predicted to be dominated by helium, carbon/oxygen, oxygen/neon/magnesium, and oxygen/silicon/sulphur. Silicon and sulphur are fused into inert iron, leading to the collapse of the core and either a supernova explosion or the direct formation of a black hole. Stripped stars, where the outer hydrogen layer has been removed and the internal He-rich layer (in Wolf-Rayet WN stars) or even the C/O layer below it (in Wolf-Rayet WC/WO stars) are exposed, provide evidence for this shell structure, and the cosmic element production mechanism it reflects. The types of supernova explosions that arise from stripped stars embedded in shells of circumstellar material (most notably Type Ibn supernovae from stars with outer He layers, and Type Icn supernovae from stars with outer C/O layers) confirm this scenario. However, direct evidence for the most interior shells, which are responsible for the production of elements heavier than oxygen, is lacking. Here, we report the discovery of the first-of-its-kind supernova arising from a star peculiarly stripped all the way to the silicon and sulphur-rich internal layer. Whereas the concentric shell structure of massive stars is not under debate, it is the first time that such a thick, massive silicon and sulphur-rich shell, expelled by the progenitor shortly before the SN explosion, has been directly revealed., Comment: 48 pages, 12 figures and 10 tables. Submitted to a high-impact journal. The reduced spectra and photometry will be made available via the journal webpage and the WISeREP archive after the acceptance of the paper

During the first half of the fourth observing run (O4a) of the International Gravitational Wave Network (IGWN), the Zwicky Transient Facility (ZTF) conducted a systematic search for kilonova (KN) counterparts to binary neutron star (BNS) and neutron star-black hole (NSBH) merger candidates. Here, we present a comprehensive study of the five high-significance (FAR < 1 per year) BNS and NSBH candidates in O4a. Our follow-up campaigns relied on both target-of-opportunity observations (ToO) and re-weighting of the nominal survey schedule to maximize coverage. We describe the toolkit we have been developing, Fritz, an instance of SkyPortal, instrumental in coordinating and managing our telescope scheduling, candidate vetting, and follow-up observations through a user-friendly interface. ZTF covered a total of 2841 deg$^2$ within the skymaps of the high-significance GW events, reaching a median depth of g~20.2 mag. We circulated 15 candidates, but found no viable KN counterpart to any of the GW events. Based on the ZTF non-detections of the high-significance events in O4a, we used a Bayesian approach, nimbus, to quantify the posterior probability of KN model parameters that are consistent with our non-detections. Our analysis favors KNe with initial absolute magnitude fainter than -16 mag. The joint posterior probability of a GW170817-like KN associated with all our O4a follow-ups was 64%. Additionally, we use a survey simulation software, simsurvey, to determine that our combined filtered efficiency to detect a GW170817-like KN is 36%, when considering the 5 confirmed astrophysical events in O3 (1 BNS and 4 NSBH), along with our O4a follow-ups. Following Kasliwal et al. (2020), we derived joint constraints on the underlying KN luminosity function based on our O3 and O4a follow-ups, determining that no more than 76% of KNe fading at 1 mag/day can peak at a magnitude brighter than -17.5 mag., Comment: submitted

Context. The complexity and variety exhibited by the light curves of long gamma-ray bursts (GRBs) enclose a wealth of information that still awaits being fully deciphered. Despite the tremendous advance in the knowledge of the energetics, structure, and composition of the relativistic jet that results from the core collapse of the progenitor star, the nature of the inner engine, how it powers the relativistic outflow, and the dissipation mechanisms remain open issues. Aims. A promising way to gain insights is describing GRB light curves as the result of a common stochastic process. In the Burst And Transient Source Experiment (BATSE) era, a stochastic pulse avalanche model was proposed and tested through the comparison of ensemble-average properties of simulated and real light curves. Here we aim to revive and further test this model. Methods. We apply it to two independent data sets, BATSE and Swift/BAT, through a machine learning approach: the model parameters are optimised using a genetic algorithm. Results. The average properties are successfully reproduced. Notwithstanding the different populations and passbands of both data sets, the corresponding optimal parameters are interestingly similar. In particular, for both sets the dynamics appears to be close to a critical state, which is key to reproduce the observed variety of time profiles. Conclusions. Our results propel the avalanche character in a critical regime as a key trait of the energy release in GRB engines, which underpins some kind of instability., Comment: Accepted by A&A, 9 pages, 3 figures. Code available at: https://github.com/LBasz/geneticgrbs/tree/arxiv-v1

The merging of a binary system involving two neutron stars (NSs), or a black hole (BH) and a NS, often results in the emission of an electromagnetic (EM) transient. One component of this EM transient is the epic explosion known as a kilonova (KN). The characteristics of the KN emission can be used to probe the equation of state (EoS) of NS matter responsible for its formation. We predict KN light curves from computationally simulated BH-NS mergers, by using the 3D radiative transfer code \texttt{POSSIS}. We investigate two EoSs spanning most of the allowed range of the mass-radius diagram. We also consider a soft EoS compatible with the observational data within the so-called 2-families scenario in which hadronic stars coexist with strange stars. Computed results show that the 2-families scenario, characterized by a soft EoS, should not produce a KN unless the mass of the binary components are small ($M_{\rm BH} \leq 6M_{\odot}$, $M_{\rm NS} \leq 1.4M_{\odot}$) and the BH is rapidly spinning ($\chi_{\rm BH} \geq 0.3$). In contrast, a strong KN signal potentially observable from future surveys (e.g. VRO/LSST) is produced in the 1-family scenario for a wider region of the parameter space, and even for non-rotating BHs ($\chi_{\rm BH} = 0$) when $M_{\rm BH} = 4M_{\odot}$ and $M_{\rm NS} = 1.2M_{\odot}$. We also provide a fit that allows for the calculation of the unbound mass from the observed KN magnitude, without running timely and costly radiative transfer simulations. Findings presented in this paper will be used to interpret light curves anticipated during the fourth observing run (O4), of the advanced LIGO, advanced Virgo and KAGRA interferometers and thus to constrain the EoS of NS matter., Comment: 15 pages, 17 figures, 2 tables; accepted for publication in MNRAS

GW170817 and its associated electromagnetic counterpart AT2017gfo continue to be a treasure trove as observations and modeling continue. Recent precision astrometry of AT2017gfo with the Hubble Space Telescope combined with previous constraints from Very Long Baseline Interferometry (VLBI) constraints narrowed down the inclination angle to 19-25 deg (90\% confidence). This paper explores how the inclusion of precise inclination information can reveal new insights about the ejecta properties, in particular, about the composition of the dynamical ejecta of AT2017gfo. Our analysis relies on updated kilonova modeling, which includes state-of-the-art heating rates, thermalization efficiencies, and opacities and is parameterized by $\bar{Y}_{\rm e,dyn}$, the average electron fraction of the dynamical ejecta component. Using this model, we incorporate the latest inclination angle constraint of AT2017gfo into a light curve fitting framework to derive updated parameter estimates. Our results suggest that the viewing angle of the observer is pointed towards the lanthanide-poor ($Y_{\rm e,dyn}\gtrsim0.25$), squeezed polar dynamical ejecta component, which can explain the early blue emission observed in the light curve of AT2017gfo. In contrast to a recent claim of spherical ejecta powering AT2017gfo, our study indicates that the composition of the dynamical ejecta has a strong angular dependence, with a lanthanide-rich ($Y_{\rm e,dyn}\lesssim0.25$), tidal component distributed around the merger plane with a half-opening angle of $35^\circ$. The inclination angle constraint reduces $\bar{Y}_{\rm e,dyn}$ from $0.24$ to $0.22$, with values $0.17\lesssim Y_{\rm e, dyn} \lesssim0.41$ enabling the robust production of $r$-process elements up to the $3^{\rm rd}$ peak in the tidal dynamical ejecta., Comment: 8 pages, 6 figures

Binary neutron star mergers and neutron star-black hole mergers are multi-messenger sources that can be detected in gravitational waves and in electromagnetic radiation. The low electron fraction of neutron star merger ejecta favors the production of heavy elements such as lanthanides and actinides via rapid neutron capture (r-process). The decay of these unstable nuclei powers an infrared-bright transient called a "kilonova". The discovery of a population of kilonovae will allow us to determine if neutron star mergers are the dominant sites for r-process element nucleosynthesis, constrain the equation of state of nuclear matter, and make independent measurements of the Hubble constant. The Nancy Grace Roman Space Telescope (Roman) will have a unique combination of depth, near-infrared sensitivity, and wide field of view. These characteristics will enable Roman's discovery of GW counterparts that will be missed by optical telescopes, such as kilonova that are associated with large distances, high lanthanide fractions, high binary mass-ratios, large dust extinction in the line of sight, or that are observed from equatorial viewing angles. Our analysis suggests to (i) make available a rapid (about 1 week) Target of Opportunity mode for GW follow-up; (ii) include observations of the High Latitude Time-Domain survey footprint in at least two filters (preferably the F158 and F213 filters) with a cadence of < 8 days; (iii) operate in synergy with Rubin Observatory. Following these recommendations, we expect that 1-6 kilonovae can be identified by Roman via ToO observations of well localized (A < 10 sq. deg., 90% C.I.) neutron star mergers during 1.5 years of the LIGO-Virgo-KAGRA fifth (or about 4-21 in during the sixth) observing run. A sample of 5-40 serendipitously discovered kilonovae can be collected in a 5-year high latitude survey., Comment: accepted for publication in Astroparticle Physics

Neutrino interactions are essential for an accurate understanding of the binary neutron star merger process. In this article, we extend the code infrastructure of the well-established numerical-relativity code BAM that until recently neglected neutrino-driven interactions. In fact, while previous work allowed already the usage of nuclear-tabulated equations of state and employing a neutrino leakage scheme, we are moving forward by implementing a first-order multipolar radiation transport scheme (M1) for the advection of neutrinos. After testing our implementation on a set of standard scenarios, we apply it to the evolution of four low-mass binary systems, and we perform an analysis of ejecta properties. We also show that our new ejecta analysis infrastructure is able to provide numerical relativity-informed inputs for the codes $\texttt{POSSIS}$ and $\texttt{Skynet}$, for the computation of kilonova lightcurves and nucleosynthesis yields, respectively., Comment: 29 pages, 27 figures

Early observations of transient explosions can provide vital clues to their progenitor origins. In this paper we present the nearby Type Iax (02cx-like) supernova (SN), SN 2020udy that was discovered within hours ($\sim$7 hr) of estimated first light. An extensive dataset of ultra-violet, optical, and near-infrared observations was obtained, covering out to $\sim$150 d after explosion. SN 2020udy peaked at -17.86$\pm$0.43 mag in the r band and evolved similarly to other 'luminous' SNe Iax, such as SNe 2005hk and 2012Z. Its well-sampled early light curve allows strict limits on companion interaction to be placed. Main-sequence companion stars with masses of 2 and 6 M$_\odot$ are ruled out at all viewing angles, while a helium-star companion is allowed from a narrow range of angles (140-180$^\circ$ away from the companion). The spectra and light curves of SN2020udy are in good agreement with those of the 'N5def' deflagration model of a near Chandrasekhar-mass carbon-oxygen white dwarf. However, as has been seen in previous studies of similar luminosity events, SN 2020udy evolves slower than the model. Broad-band linear polarisation measurements taken at and after peak are consistent with no polarisation, in agreement with the predictions of the companion-star configuration from the early light curve measurements. The host galaxy environment is low metallicity and is consistent with a young stellar population. Overall, we find the most plausible explosion scenario to be the incomplete disruption of a CO white dwarf near the Chandrasekhar-mass limit, with a helium-star companion., Comment: 18 pages, 14 figures, submitted to MNRAS

Over the last few years, there has been an increasing interest in sub-solar mass black holes due to their potential to provide valuable information about cosmology or the black hole population. Motivated by this, we study observable phenomena connected to the merger of a sub-solar mass black hole with a neutron star. For this purpose, we perform new numerical-relativity simulations of a binary system composed of a black hole with mass $0.5M_\odot$ and a neutron star with mass $1.4 M_\odot$. We investigate the merger dynamics of this exotic system and provide information about the connected gravitational-wave and kilonova signals. Our study indicates that current gravitational-waveform models are unable to adequately describe such systems and that phenomenological relations connecting the binary parameters with the ejecta and remnant properties are not applicable to our system. Furthermore, we find a dependence of the kilonova signal on the azimuthal viewing angle due to the asymmetric mass ejection. This first-of-its-kind simulation opens the door for the study of sub-solar mass black hole - neutron star mergers and could serve as a testing ground for future model development., Comment: 16 pages, 14 figures, to be submitted to PRD, comments welcome

An open question in SN Ia research is where the boundary lies between 'normal' Type Ia supernovae (SNe Ia) that are used in cosmological measurements and those that sit off the Phillips relation. We present the spectroscopic modelling of one such '86G-like' transitional SN Ia, SN 2021rhu, that has recently been employed as a local Hubble Constant calibrator using a tip of the red-giant branch measurement. We detail its modelling from -12 d until maximum brightness using the radiative-transfer spectral-synthesis code tardis. We base our modelling on literature delayed-detonation and deflagration models of Chandrasekhar mass white dwarfs, as well as the double-detonation models of sub-Chandrasekhar mass white dwarfs. We present a new method for 'projecting' abundance profiles to different density profiles for ease of computation. Due to the small velocity extent and low outer densities of the W7 profile, we find it inadequate to reproduce the evolution of SN 2021rhu as it fails to match the high-velocity calcium components. The host extinction of SN 2021rhu is uncertain but we use modelling with and without an extinction correction to set lower and upper limits on the abundances of individual species. Comparing these limits to literature models we conclude that the spectral evolution of SN 2021rhu is also incompatible with double-detonation scenarios, lying more in line with those resulting from the delayed detonation mechanism (although there are some discrepancies, in particular a larger titanium abundance in SN 2021rhu compared to the literature). This suggests that SN 2021rhu is likely a lower luminosity, and hence lower temperature, version of a normal SN Ia., Comment: 25 pages, 22 figures, accepted for publication in MNRAS

A merger of binary neutron stars creates heavy unstable elements whose radioactive decay produces a thermal emission known as a kilonova. In this paper, we predict the photometric and polarimetric behaviour of this emission by performing 3-D Monte Carlo radiative transfer simulations. In particular, we choose three hydrodynamical models for merger ejecta, two including jets with different luminosities and one without a jet structure, to help decipher the impact of jets on the light curve and polarimetric behaviour. In terms of photometry, we find distinct color evolutions across the three models. Models without a jet show the highest variation in light curves for different viewing angles. In contrast, to previous studies, we find models with a jet to produce fainter kilonovae when viewed from orientations close to the jet axis, compared to a model without a jet. In terms of polarimetry, we predict relatively low levels (<~0.3-0.4%) at all orientations that, however, remain non-negligible until a few days after the merger and longer than previously found. Despite the low levels, we find that the presence of a jet enhances the degree of polarization at wavelengths ranging from 0.25 to 2.5\micron, an effect that is found to increase with the jet luminosity. Thus, future photometric and polarimetric campaigns should observe kilonovae in blue and red filters for a few days after the merger to help constrain the properties of the ejecta (e.g. composition) and jet., Comment: 11 pages, 7 figures, accepted for publication and in press

Although GRB 211211A is one of the closest gamma-ray bursts (GRBs), its classification is challenging because of its partially inconclusive electromagnetic signatures. In this paper, we investigate four different astrophysical scenarios as possible progenitors for GRB~211211A: a binary neutron-star merger, a black-hole--neutron-star merger, a core-collapse supernova, and an r-process enriched core collapse of a rapidly rotating massive star (a collapsar). We perform a large set of Bayesian multi-wavelength analyses based on different models describing these scenarios and priors to investigate which astrophysical scenarios and processes might be related to GRB~211211A. Our analysis supports previous studies in which the presence of an additional component, likely related to $r$-process nucleosynthesis, is required to explain the observed light curves of GRB~211211A, as it can not solely be explained as a GRB afterglow. Fixing the distance to about $350~\rm Mpc$, namely the distance of the possible host galaxy SDSS J140910.47+275320.8, we find a statistical preference for a binary neutron-star merger scenario., Comment: 14 pages, 6 figures

Tidal disruption events (TDEs) occur when a star orbiting a massive black hole is sufficiently close to be tidally ripped apart by the black hole. AT 2022cmc is the first relativistic TDE that was observed (and discovered) as an optically bright and fast transient, showing signatures of non-thermal radiation induced by a jet which is oriented towards the Earth. In this work, we present optical linear and circular polarization measurements, observed with VLT/FORS2 in the $R$-band (which corresponds to the blue/UV part of the spectrum in rest frame), $\sim$ 7.2 and $\sim$ 12.2 rest-frame days after the first detection, respectively, when the light curve of the transient had settled in a bright blue plateau. Both linear and circular polarization are consistent with zero, $p_{lin}$ = 0.14 $\pm$ 0.73 % and $p_{cir}$ = $-$0.30 $\pm$ 0.53 %. This is the highest S/N linear polarization measurement obtained for a relativistic TDE and the first circular polarimetry for such a transient. The non detection of the linear and circular polarization is consistent with the scenario of AT 2022cmc being a TDE where the thermal component (disk+outflows) is viewed pole-on, assuming an axially symmetric geometry. The presence and effect of a jet and/or external shocks are, however, difficult to disentangle., Comment: Accepted for publication in ApJL

TDEs have been observed in the optical and UV for more than a decade but the underlying emission mechanism still remains a puzzle. It has been suggested that viewing angle effects could potentially explain their large photometric and spectroscopic diversity. Polarization is indeed sensitive to the viewing angle and the first polarimetry studies of TDEs are now available, calling for a theoretical interpretation. In this study, we model the continuum polarization levels of TDEs using the radiative transfer code POSSIS and the collision-induced outflow (CIO) TDE emission scenario where unbound shocked gas originating from a debris stream intersection point offset from the black hole, reprocesses the hard emission from the accretion flow into UV and optical bands. We explore two different cases of peak mass fallback rates M'p (~3 and ~0.3 Msol/yr) while varying the following geometrical parameters: the distance R_int from the black hole (BH) to the intersection point, the radius of the photosphere around the BH R_ph, on the surface of which the photons are generated, and the opening angle Deltheta (anisotropic emission). For the high mass fallback rate case, we find for every viewing angle polarization levels below one (P<1%) and P<0.5% for 10/12 simulations. The absolute value of polarization reaches its maximum (P_max) for equatorial viewing angles. For the low mass fallback rate case, the maximum value predicted is P~8.8% and P_max is reached for intermediate viewing angles. We find that the polarization depends strongly on i) the optical depths at the central regions set by the different M'p values and ii) the viewing angle. Finally, by comparing our model predictions to polarization observations of a few TDEs, we attempt to constrain their observed viewing angles and we show that multi-epoch polarimetric observations can become a key factor in constraining the viewing angle of TDEs., Comment: Accepted for publication in Astronomy and Astrophysics journal; 20 pages

Narrow transient emission lines (flash-ionization features) in early supernova (SN) spectra trace the presence of circumstellar material (CSM) around the massive progenitor stars of core-collapse SNe. The lines disappear within days after the SN explosion, suggesting that this material is spatially confined, and originates from enhanced mass loss shortly (months to a few years) prior to explosion. We performed a systematic survey of H-rich (Type II) SNe discovered within less than two days from explosion during the first phase of the Zwicky Transient Facility (ZTF) survey (2018-2020), finding thirty events for which a first spectrum was obtained within $< 2$ days from explosion. The measured fraction of events showing flash ionisation features ($>36\%$ at $95\%$ confidence level) confirms that elevated mass loss in massive stars prior to SN explosion is common. We find that SNe II showing flash ionisation features are not significantly brighter, nor bluer, nor more slowly rising than those without. This implies that CSM interaction does not contribute significantly to their early continuum emission, and that the CSM is likely optically thin. We measured the persistence duration of flash ionisation emission and find that most SNe show flash features for $\approx 5 $ days. Rarer events, with persistence timescales $>10$ days, are brighter and rise longer, suggesting these may be intermediate between regular SNe II and strongly-interacting SNe IIn.

Tidal disruption events (TDEs) are bursts of electromagnetic energy released when supermassive black holes (SMBHs) at the centers of galaxies violently disrupt a star that passes too close. TDEs provide a new window to study accretion onto SMBHs; in some rare cases, this accretion leads to launching of a relativistic jet, but the necessary conditions are not fully understood. The best studied jetted TDE to date is Swift J1644+57, which was discovered in gamma-rays, but was too obscured by dust to be seen at optical wavelengths. Here we report the optical discovery of AT2022cmc, a rapidly fading source at cosmological distance (redshift z=1.19325) whose unique lightcurve transitioned into a luminous plateau within days. Observations of a bright counterpart at other wavelengths, including X-rays, sub-millimeter, and radio, supports the interpretation of AT2022cmc as a jetted TDE containing a synchrotron "afterglow", likely launched by a SMBH with spin $a \gtrsim 0.3$. Using 4 years of Zwicky Transient Facility (ZTF) survey data, we calculate a rate of $0.02 ^{+ 0.04 }_{- 0.01 }$ Gpc$^{-3}$ yr$^{-1}$ for on-axis jetted TDEs based on the luminous, fast-fading red component, thus providing a measurement complementary to the rates derived from X-ray and radio observations. Correcting for the beaming angle effects, this rate confirms that about 1% of TDEs have relativistic jets. Optical surveys can use AT2022cmc as a prototype to unveil a population of jetted TDEs., Comment: Accepted for publication in Nature

We present an improved version of the 3D Monte Carlo radiative transfer code POSSIS to model kilonovae from neutron star mergers, wherein nuclear heating rates, thermalization efficiencies and wavelength-dependent opacities depend on local properties of the ejecta and time. Using an axially-symmetric two-component ejecta model, we explore how simplistic assumptions on heating rates, thermalization efficiencies and opacities often found in the literature affect kilonova spectra and light curves. Specifically, we compute five models: one ($\texttt{FIDUCIAL}$) with an appropriate treatment of these three quantities, one ($\texttt{SIMPLE-HEAT}$) with uniform heating rates throughout the ejecta, one ($\texttt{SIMPLE-THERM}$) with a constant and uniform thermalization efficiency, one ($\texttt{SIMPLE-OPAC}$) with grey opacities and one ($\texttt{SIMPLE-ALL}$) with all these three simplistic assumptions combined. We find that deviations from the $\texttt{FIDUCIAL}$ model are of several ($\sim1-10$) magnitudes and are generally larger for the $\texttt{SIMPLE-OPAC}$ and $\texttt{SIMPLE-ALL}$ compared to the $\texttt{SIMPLE-THERM}$ and $\texttt{SIMPLE-HEAT}$ models. The discrepancies generally increase from a face-on to an edge-on view of the system, from early to late epochs and from infrared to ultraviolet/optical wavelengths. Our work indicates that kilonova studies using either of these simplistic assumptions ought to be treated with caution and that appropriate systematic uncertainties ought to be added to kilonova light curves when performing inference on ejecta parameters., Comment: 13 pages, 11 figures, 1 table; accepted for publication in MNRAS; added missing citations

Investigating interstellar (IS) dust properties in external galaxies is important not only to infer the intrinsic properties of astronomical objects but also to understand the star/planet formation in the galaxies. From the non-Milky-Way-like extinction and interstellar polarization (ISP) observed in reddened Type Ia supernovae (SNe), it has been suggested that their host galaxies contain dust grains whose properties are substantially different from the Milky-Way (MW) dust. It is important to investigate the universality of such non-MW-like dust in the universe. Here we report spectropolarimetry of two highly-extinguished Type II SNe (SN 2022aau and SN 2022ame). SN 2022aau shows a polarization maximum at a shorter wavelength than MW stars, which is also observed in some Type Ia SNe. This is clear evidence for the existence of non-MW-like dust in its host galaxy (i.e., NGC 1672). This fact implies that such non-MW-like dust might be more common in some environments than expected, and thus it might affect the picture of the star/planet formation. On the other hand, SN 2022ame shows MW-like ISP, implying the presence of MW-like dust in its host galaxy (i.e., NGC 1255). Our findings confirm that dust properties of galaxies are diverse, either locally or globally. The present work demonstrates that further investigation of IS dust properties in external galaxies using polarimetry of highly-reddened SNe is promising, providing a great opportunity to study the universality of such non-MW-like dust grains in the universe., Comment: 11 pages, 5 figures, accepted for publication in ApJ Letters

A rare class of supernovae (SNe) is characterized by strong interaction between the ejecta and several solar masses of circumstellar matter (CSM) as evidenced by strong Balmer-line emission. Within the first few weeks after the explosion, they may display spectral features similar to overluminous Type Ia SNe, while at later phase their observation properties exhibit remarkable similarities with some extreme case of Type IIn SNe that show strong Balmer lines years after the explosion. We present polarimetric observations of SN2018evt obtained by the ESO Very Large Telescope from 172 to 219 days after the estimated time of peak luminosity to study the geometry of the CSM. The nonzero continuum polarization decreases over time, suggesting that the mass loss of the progenitor star is aspherical. The prominent H$\alpha$ emission can be decomposed into a broad, time-evolving component and an intermediate-width, static component. The former shows polarized signals, and it is likely to arise from a cold dense shell (CDS) within the region between the forward and reverse shocks. The latter is significantly unpolarized, and it is likely to arise from shocked, fragmented gas clouds in the H-rich CSM. We infer that SN2018evt exploded inside a massive and aspherical circumstellar cloud. The symmetry axes of the CSM and the SN appear to be similar. SN\,2018evt shows observational properties common to events that display strong interaction between the ejecta and CSM, implying that they share similar circumstellar configurations. Our preliminary estimate also suggests that the circumstellar environment of SN2018evt has been significantly enriched at a rate of $\sim0.1$ M$_\odot$ yr$^{-1}$ over a period of $>100$ yr., Comment: 29 pages, 20 figures, 10 tables, accepted for publication in MNRAS

Accurate numerical-relativity simulations are essential to study the rich phenomenology of binary neutron star systems. In this work, we focus on the material that is dynamically ejected during the merger process and on the kilonova transient it produces. Typically, radiative transfer simulations of kilonova light curves from ejecta make the assumption of homologous expansion, but this condition might not always be met at the end of usually very short numerical-relativity simulations. In this article, we adjust the infrastructure of the BAM code to enable longer simulations of the dynamical ejecta with the aim of investigating when the condition of homologous expansion is satisfied. In fact, we observe that the deviations from a perfect homologous expansion are about 30% at roughly 100ms after the merger. To determine how these deviations might affect the calculation of kilonova light curves, we extract the ejecta data for different reference times and use them as input for radiative transfer simulations. Our results show that the light curves for extraction times later than 80ms after the merger deviate by less than 0.4mag and are mostly consistent with numerical noise. Accordingly, deviations from the homologous expansion for the dynamical ejecta component are negligible for the purpose of kilonova modelling.

A star crossing the tidal radius of a supermassive black hole will be spectacularly ripped apart with an accompanying burst of radiation. A few tens of such tidal disruption events (TDEs) have now been identified in the optical wavelengths, but the exact origin of the strong optical emission remains inconclusive. Here we report polarimetric observations of three TDEs. The continuum polarization is independent of wavelength, while emission lines are partially depolarized. These signatures are consistent with optical photons being scattered and polarized in an envelope of free electrons. An almost axisymmetric photosphere viewed from different angles is in broad agreement with the data, but there is also evidence for deviations from axial symmetry before the peak of the flare and significant time evolution at early times, compatible with the rapid formation of an accretion disk. By combining a super-Eddington accretion model with a radiative transfer code we generate predictions for the degree of polarization as a function of disk mass and viewing angle, and we show that the predicted levels are compatible with the observations, for extended reprocessing envelopes of $\sim$1000 gravitational radii. Spectropolarimetry therefore constitutes a new observational test for TDE models, and opens an important new line of exploration in the study of TDEs., Comment: Author's version of paper to appear in Nature Astronomy. In the journal version the detailed discussion on the ISP determination will be moved from the Methods section to a Supplementary Information section. 58 pages in double spacing format, including 5 Figures, 10 Extended Data Figures and 2 Tables

The simultaneous detection of gravitational waves and light from the binary neutron star merger GW170817 led to independent measurements of distance and redshift, providing a direct estimate of the Hubble constant $H_0$ that does not rely on a cosmic distance ladder nor assumes a specific cosmological model. By using gravitational waves as ''standard sirens'', this approach holds promise to arbitrate the existing tension between the $H_0$ value inferred from the cosmic microwave background and those obtained from local measurements. However, the known degeneracy in the gravitational-wave analysis between distance and inclination of the source lead to a $H_0$ value from GW170817 that was not precise enough to resolve the existing tension. In this review, we summarize recent works exploiting the viewing-angle dependence of the electromagnetic signal, namely the associated short gamma-ray burst and kilonova, to constrain the system inclination and improve on $H_0$. We outline the key ingredients of the different methods, summarize the results obtained in the aftermath of GW170817 and discuss the possible systematics introduced by each of these methods., Comment: invited review accepted for publication in the journal Universe as part of the special issue 'Gamma-Ray Bursts: Observational and Theoretical Prospects in the Era of Multi-Messenger Astronomy', 23 pages, 3 figures, 1 table

The multi-messenger detection of the gravitational-wave signal GW170817, the corresponding kilonova AT2017gfo and the short gamma-ray burst GRB170817A, as well as the observed afterglow has delivered a scientific breakthrough. For an accurate interpretation of all these different messengers, one requires robust theoretical models that describe the emitted gravitational-wave, the electromagnetic emission, and dense matter reliably. In addition, one needs efficient and accurate computational tools to ensure a correct cross-correlation between the models and the observational data. For this purpose, we have developed the Nuclear-physics and Multi-Messenger Astrophysics framework NMMA. The code allows incorporation of nuclear-physics constraints at low densities as well as X-ray and radio observations of isolated neutron stars. In previous works, the NMMA code has allowed us to constrain the equation of state of supranuclear dense matter, to measure the Hubble constant, and to compare dense-matter physics probed in neutron-star mergers and in heavy-ion collisions, and to classify electromagnetic observations and perform model selection. Here, we show an extension of the NMMA code as a first attempt of analyzing the gravitational-wave signal, the kilonova, and the gamma-ray burst afterglow simultaneously. Incorporating all available information, we estimate the radius of a $1.4M_\odot$ neutron star to be $R=11.98^{+0.35}_{-0.40}$km., Comment: code available at https://github.com/nuclear-multimessenger-astronomy

The Fermi Gamma-ray Burst Monitor (GBM) triggers on-board in response to $\sim$ 40 short gamma-ray bursts (SGRBs) per year; however, their large localization regions have made the search for optical counterparts a challenging endeavour. We have developed and executed an extensive program with the wide field of view of the Zwicky Transient Facility (ZTF) camera, mounted on the Palomar 48 inch Oschin telescope (P48), to perform target-of-opportunity (ToO) observations on 10 Fermi-GBM SGRBs during 2018 and 2020-2021. Bridging the large sky areas with small field of view optical telescopes in order to track the evolution of potential candidates, we look for the elusive SGRB afterglows and kilonovae (KNe) associated with these high-energy events. No counterpart has yet been found, even though more than 10 ground based telescopes, part of the Global Relay of Observatories Watching Transients Happen (GROWTH) network, have taken part in these efforts. The candidate selection procedure and the follow-up strategy have shown that ZTF is an efficient instrument for searching for poorly localized SGRBs, retrieving a reasonable number of candidates to follow-up and showing promising capabilities as the community approaches the multi-messenger era. Based on the median limiting magnitude of ZTF, our searches would have been able to retrieve a GW170817-like event up to $\sim$ 200 Mpc and SGRB afterglows to z = 0.16 or 0.4, depending on the assumed underlying energy model. Future ToOs will expand the horizon to z = 0.2 and 0.7 respectively.

The last decade has brought about a profound transformation in multimessenger science. Ten years ago, facilities had been built or were under construction that would eventually discover the nature of objects in our universe could be detected through multiple messengers. Nonetheless, multimessenger science was hardly more than a dream. The rewards for our foresight were finally realized through IceCube's discovery of the diffuse astrophysical neutrino flux, the first observation of gravitational waves by LIGO, and the first joint detections in gravitational waves and photons and in neutrinos and photons. Today we live in the dawn of the multimessenger era. The successes of the multimessenger campaigns of the last decade have pushed multimessenger science to the forefront of priority science areas in both the particle physics and the astrophysics communities. Multimessenger science provides new methods of testing fundamental theories about the nature of matter and energy, particularly in conditions that are not reproducible on Earth. This white paper will present the science and facilities that will provide opportunities for the particle physics community renew its commitment and maintain its leadership in multimessenger science., Comment: 174 pages, 12 figures. Contribution to Snowmass 2021. Solicited white paper from CF07. Comments and endorsers welcome. Still accepting contributions (contact editors)

On the 17th of August, 2017 came the simultaneous detections of GW170817, a gravitational wave that originated from the coalescence of two neutron stars, along with the gamma-ray burst GRB170817A, and the kilonova counterpart AT2017gfo. Since then, there has been much excitement surrounding the study of neutron star mergers, both observationally, using a variety of tools, and theoretically, with the development of complex models describing the gravitational-wave and electromagnetic signals. In this work, we improve upon our pipeline to infer kilonova properties from observed light-curves by employing a Neural-Network framework that reduces execution time and handles much larger simulation sets than previously possible. In particular, we use the radiative transfer code POSSIS to construct 5-dimensional kilonova grids where we employ different functional forms for the angular dependence of the dynamical ejecta component. We find that incorporating an angular dependence improves the fit to the AT2017gfo light-curves by up to ~50% when quantified in terms of the weighted Mean Square Error.

Some highly reddened Type Ia supernovae (SNe Ia) display low total-to-selective extinction ratios ($R_V \lesssim 2$) in comparison to that of typical Milky Way dust ($R_V \approx 3.3$), and polarization curves that rise steeply to blue wavelengths, with peak polarization values at short wavelengths ($\lambda_{\rm max} < 0.4$ ${\mu}$m) in comparison to the typical Galactic values ($\lambda_{\rm max} \approx 0.55$ ${\mu}$m). Understanding the source of these properties could provide insight into the progenitor systems of SNe Ia. We aim to determine whether they are the result of the host galaxy's interstellar dust or circumstellar dust. This is accomplished by analysing the continuum polarization of 66 SNe Ia in dust-rich spiral galaxies and 13 SNe Ia in dust-poor elliptical galaxies as a function of normalised galactocentric distance. We find that there is a general trend of SNe Ia in spiral galaxies displaying increased polarization values when located closer to the host galaxies' centre, while SNe Ia in elliptical host galaxies display low polarization. Furthermore, all highly polarized SNe Ia in spiral host galaxies display polarization curves rising toward blue wavelengths, while no evidence of such polarization properties is shown in elliptical host galaxies. This indicates that the source of the peculiar polarization curves is likely the result of interstellar material as opposed to circumstellar material. The peculiar polarization and extinction properties observed toward some SNe Ia may be explained by the radiative torque disruption mechanism induced by the SN or the interstellar radiation field., Comment: Accepted for publication in MNRAS

The discovery of the electromagnetic counterpart to the binary neutron star merger GW170817 has opened the era of gravitational-wave multi-messenger astronomy. Rapid identification of the optical/infrared kilonova enabled a precise localization of the source, which paved the way to deep multi-wavelength follow-up and its myriad of related science results. Fully exploiting this new territory of exploration requires the acquisition of electromagnetic data from samples of neutron star mergers and other gravitational wave sources. After GW170817, the frontier is now to map the diversity of kilonova properties and provide more stringent constraints on the Hubble constant, and enable new tests of fundamental physics. The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) can play a key role in this field in the 2020s, when an improved network of gravitational-wave detectors is expected to reach a sensitivity that will enable the discovery of a high rate of merger events involving neutron stars (about tens per year) out to distances of several hundred Mpc. We design comprehensive target-of-opportunity observing strategies for follow-up of gravitational-wave triggers that will make the Rubin Observatory the premier instrument for discovery and early characterization of neutron star and other compact object mergers, and yet unknown classes of gravitational wave events., Comment: Accepted for publication in ApJS. arXiv admin note: text overlap with arXiv:1812.04051

Detailed spectropolarimetric studies may hold the key to probing the explosion mechanisms and the progenitor scenarios of Type Ia supernovae (SNe Ia). We present multi-epoch spectropolarimetry and imaging polarimetry of SN 2019ein, an SN Ia showing high expansion velocities at early phases. The spectropolarimetry sequence spans from $\sim -11$ to $+$10 days relative to peak brightness in the $B$-band. We find that the level of the continuum polarization of SN 2019ein, after subtracting estimated interstellar polarization, is in the range $0.0-0.3\%$, typical for SNe Ia. The polarization position angle remains roughly constant before and after the SN light-curve peak, implying that the inner regions share the same axisymmetry as the outer layers. We observe high polarization ($\sim 1\%$) across both the Si II $\lambda6355$ and Ca II near-infrared triplet features. These two lines also display complex polarization modulations. The spectropolarimetric properties of SN 2019ein rule out a significant departure from spherical symmetry of the ejecta for up to a month after the explosion. These observations disfavour merger-induced and double-detonation models for SN 2019ein. The imaging polarimetry shows weak evidence for a modest increase in polarization after $\sim 20$ days since the $B$-band maximum. If this rise is real and is observed in other SNe Ia at similar phases, we may have seen, for the first time, an aspherical interior similar to what has been previously observed for SNe IIP. Future polarization observations of SNe Ia extending to post-peak epochs will help to examine the inner structure of the explosion., Comment: Accepted by MNRAS. 13 pages, 7 figures. Updated Figure 6

Searches for electromagnetic counterparts of gravitational-wave signals have redoubled since the first detection in 2017 of a binary neutron star merger with a gamma-ray burst, optical/infrared kilonova, and panchromatic afterglow. Yet, one LIGO/Virgo observing run later, there has not yet been a second, secure identification of an electromagnetic counterpart. This is not surprising given that the localization uncertainties of events in LIGO and Virgo's third observing run, O3, were much larger than predicted. We explain this by showing that improvements in data analysis that now allow LIGO/Virgo to detect weaker and hence more poorly localized events have increased the overall number of detections, of which well-localized, gold-plated events make up a smaller proportion overall. We present simulations of the next two LIGO/Virgo/KAGRA observing runs, O4 and O5, that are grounded in the statistics of O3 public alerts. To illustrate the significant impact that the updated predictions can have, we study the follow-up strategy for the Zwicky Transient Facility. Realistic and timely forecasting of gravitational-wave localization accuracy is paramount given the large commitments of telescope time and the need to prioritize which events are followed up. We include a data release of our simulated localizations as a public proposal planning resource for astronomers., Comment: ApJ, in press

This work aims to study different probes of Type Ia supernova progenitors that have been suggested to be linked to the presence of circumstellar material (CSM). In particular, we have investigated, for the first time, the link between narrow blueshifted NaID absorption profiles and the presence and strength of the broad high-velocity CaII near infrared triplet absorption features seen in Type Ia supernovae around maximum light. With the probes exploring different distances from the supernova; NaID > 10$^{17}$cm, high-velocity CaII features < 10$^{15}$cm. For this, we have used a new intermediate-resolution X-shooter spectral sample of 15 Type Ia supernovae. We do not identify a link between these two probes, implying either that, one (or both) is not physically related to the presence of CSM or that the occurrence of CSM at the distance explored by one probe is not linked to its presence at the distance probed by the other. However, the previously identified statistical excess in the presence of blueshifted (over redshifted) NaID absorption is confirmed in this sample at high significance and is found to be stronger in Type Ia supernovae hosted by late-type galaxies. This excess is difficult to explain as being from an interstellar-medium origin as has been suggested by some recent modelling, as such an origin is not expected to show a bias for blueshifted absorption. However, a circumstellar origin for these features also appears unsatisfactory based on our new results given the lack of link between the two probes of CSM investigated., Comment: This is a pre-copyedited, author-produced PDF of an article accepted for publication in Monthly Notices of the Royal Astronomical Society following peer review. 25 pages, 16 Figures

We present ${\tt nimbus}$ : a hierarchical Bayesian framework to infer the intrinsic luminosity parameters of kilonovae (KNe) associated with gravitational-wave (GW) events, based purely on non-detections. This framework makes use of GW 3-D distance information and electromagnetic upper limits from multiple surveys for multiple events, and self-consistently accounts for finite sky-coverage and probability of astrophysical origin. The framework is agnostic to the brightness evolution assumed and can account for multiple electromagnetic passbands simultaneously. Our analyses highlight the importance of accounting for model selection effects, especially in the context of non-detections. We show our methodology using a simple, two-parameter linear brightness model, taking the follow-up of GW190425 with the Zwicky Transient Facility (ZTF) as a single-event test case for two different prior choices of model parameters -- (i) uniform/uninformative priors and (ii) astrophysical priors based on surrogate models of Monte Carlo radiative transfer simulations of KNe. We present results under the assumption that the KN is within the searched region to demonstrate functionality and the importance of prior choice. Our results show consistency with ${\tt simsurvey}$ -- an astronomical survey simulation tool used previously in the literature to constrain the population of KNe. While our results based on uniform priors strongly constrain the parameter space, those based on astrophysical priors are largely uninformative, highlighting the need for deeper constraints. Future studies with multiple events having electromagnetic follow-up from multiple surveys should make it possible to constrain the KN population further., Comment: 18 pages, 5 figures, accepted for publication in ApJ

Current and future optical and near-infrared wide-field surveys have the potential of finding kilonovae, the optical and infrared counterparts to neutron star mergers, independently of gravitational-wave or high-energy gamma-ray burst triggers. The ability to discover fast and faint transients such as kilonovae largely depends on the area observed, the depth of those observations, the number of re-visits per field in a given time frame, and the filters adopted by the survey; it also depends on the ability to perform rapid follow-up observations to confirm the nature of the transients. In this work, we assess kilonova detectability in existing simulations of the LSST strategy for the Vera C. Rubin Wide Fast Deep survey, with focus on comparing rolling to baseline cadences. Although currently available cadences can enable the detection of more than 300 kilonovae out to 1400 Mpc over the ten-year survey, we can expect only 3-32 kilonovae similar to GW170817 to be recognizable as fast-evolving transients. We also explore the detectability of kilonovae over the plausible parameter space, focusing on viewing angle and ejecta masses. We find that observations in redder izy bands are crucial for identification of nearby (within 300 Mpc) kilonovae that could be spectroscopically classified more easily than more distant sources. Rubin's potential for serendipitous kilonova discovery could be increased by gain of efficiency with the employment of individual 30s exposures (as opposed to 2x15s snap pairs), with the addition of red-band observations coupled with same-night observations in g- or r-bands, and possibly with further development of a new rolling-cadence strategy., Comment: Submitted to ApJS

In the past few years, new observations of neutron stars and neutron-star mergers have provided a wealth of data that allow one to constrain the equation of state of nuclear matter at densities above nuclear saturation density. However, most observations were based on neutron stars with masses of about 1.4 solar masses, probing densities up to $\sim$ 3-4 times the nuclear saturation density. Even higher densities are probed inside massive neutron stars such as PSR J0740+6620. Very recently, new radio observations provided an update to the mass estimate for PSR J0740+6620 and X-ray observations by the NICER and XMM telescopes constrained its radius. Based on these new measurements, we revisit our previous nuclear-physics multi-messenger astrophysics constraints and derive updated constraints on the equation of state describing the neutron-star interior. By combining astrophysical observations of two radio pulsars, two NICER measurements, the two gravitational-wave detections GW170817 and GW190425, detailed modeling of the kilonova AT2017gfo, as well as the gamma-ray burst GRB170817A, we are able to estimate the radius of a typical 1.4-solar mass neutron star to be $11.94^{+0.76}_{-0.87} \rm{km}$ at $90\%$ confidence. Our analysis allows us to revisit the upper bound on the maximum mass of neutron stars and disfavors the presence of a strong first-order phase transition from nuclear matter to exotic forms of matter, such as quark matter, inside neutron stars., Comment: 13 pages, 5 figures

Gamma-ray bursts (GRBs) are among the brightest and most energetic events in the universe. The duration and hardness distribution of GRBs has two clusters, now understood to reflect (at least) two different progenitors. Short-hard GRBs (SGRBs; T90 <2 s) arise from compact binary mergers, while long-soft GRBs (LGRBs; T90 >2 s) have been attributed to the collapse of peculiar massive stars (collapsars). The discovery of SN 1998bw/GRB 980425 marked the first association of a LGRB with a collapsar and AT 2017gfo/GRB 170817A/GW170817 marked the first association of a SGRB with a binary neutron star merger, producing also gravitational wave (GW). Here, we present the discovery of ZTF20abwysqy (AT2020scz), a fast-fading optical transient in the Fermi Satellite and the InterPlanetary Network (IPN) localization regions of GRB 200826A; X-ray and radio emission further confirm that this is the afterglow. Follow-up imaging (at rest-frame 16.5 days) reveals excess emission above the afterglow that cannot be explained as an underlying kilonova (KN), but is consistent with being the supernova (SN). Despite the GRB duration being short (rest-frame T90 of 0.65 s), our panchromatic follow-up data confirms a collapsar origin. GRB 200826A is the shortest LGRB found with an associated collapsar; it appears to sit on the brink between a successful and a failed collapsar. Our discovery is consistent with the hypothesis that most collapsars fail to produce ultra-relativistic jets., Comment: submitted to Nature Astronomy

While optical surveys regularly discover slow transients like supernovae on their own, the most common way to discover extragalactic fast transients, fading away in a few nights, is via follow-up observations of gamma-ray burst and gravitational-wave triggers. However, wide-field surveys have the potential to also identify rapidly fading transients independently of such external triggers. The volumetric survey speed of the Zwicky Transient Facility (ZTF) makes it sensitive to faint and fast-fading objects as kilonovae, the optical counterparts to binary neutron stars and neutron star-black hole mergers, out to almost 200Mpc. We introduce an open-source software infrastructure, the ZTF REaltime Search and Triggering, ZTFReST, designed to identify kilonovae and fast optical transients in ZTF data. Using the ZTF alert stream combined with forced photometry, we have implemented automated candidate ranking based on their photometric evolution and fitting to kilonova models. Automated triggering of follow-up systems, such as Las Cumbres Observatory, has also been implemented. In 13 months of science validation, we found several extragalactic fast transients independent of any external trigger (though some counterparts were identified later), including at least one supernova with post-shock cooling emission, two known afterglows with an associated gamma-ray burst, two known afterglows without any known gamma-ray counterpart, and three new fast-declining sources (ZTF20abtxwfx, ZTF20acozryr, and ZTF21aagwbjr) that are likely associated with GRB200817A, GRB201103B, and GRB210204A. However, we have not found any objects which appear to be kilonovae; therefore, we constrain the rate of GW170817-like kilonovae to $R < 900$Gpc$^{-3}$yr$^{-1}$. A framework such as ZTFReST could become a prime tool for kilonova and fast transient discovery with the Vera C. Rubin Observatory., Comment: Submitted for publication

Searches for gravitational-wave counterparts have been going in earnest since GW170817 and the discovery of AT2017gfo. Since then, the lack of detection of other optical counterparts connected to binary neutron star or black hole - neutron star candidates has highlighted the need for a better discrimination criterion to support this effort. At the moment, the low-latency gravitational-wave alerts contain preliminary information about the binary properties and, hence, on whether a detected binary might have an electromagnetic counterpart. The current alert method is a classifier that estimates the probability that there is a debris disc outside the black hole created during the merger as well as the probability of a signal being a binary neutron star, a black hole - neutron star, a binary black hole or of terrestrial origin. In this work, we expand upon this approach to predict both the ejecta properties and provide contours of potential lightcurves for these events in order to improve follow-up observation strategy. The various sources of uncertainty are discussed, and we conclude that our ignorance about the ejecta composition and the insufficient constraint of the binary parameters, by the low-latency pipelines, represent the main limitations. To validate the method, we test our approach on real events from the second and third Advanced LIGO-Virgo observing runs.

In recent years, there have been significant advances in multi-messenger astronomy due to the discovery of the first, and so far only confirmed, gravitational wave event with a simultaneous electromagnetic (EM) counterpart, as well as improvements in numerical simulations, gravitational wave (GW) detectors, and transient astronomy. This has led to the exciting possibility of performing joint analyses of the GW and EM data, providing additional constraints on fundamental properties of the binary progenitor and merger remnant. Here, we present a new Bayesian framework that allows inference of these properties, while taking into account the systematic modeling uncertainties that arise when mapping from GW binary progenitor properties to photometric light curves. We extend the relative binning method presented in Zackay et al. (2018) to include extrinsic GW parameters for fast analysis of the GW signal. The focus of our EM framework is on light curves arising from r-process nucleosynthesis in the ejected material during and after merger, the so called kilonova, and particularly on black hole - neutron star systems. As a case study, we examine the recent detection of GW190425, where the primary object is consistent with being either a black hole (BH) or a neutron star (NS). We show quantitatively how improved mapping between binary progenitor and outflow properties, and/or an increase in EM data quantity and quality are required in order to break degeneracies in the fundamental source parameters., Comment: 22 pages, 12 figures, comments are welcome

Massive stars die an explosive death as a core-collapse supernova (CCSN). The exact physical processes that cause the collapsing star to rebound into an explosion are not well-understood, and the key in resolving this issue may lie in the measurement of the shape of CCSNe ejecta. Spectropolarimetry is the only way to perform this measurement for CCSNe outside of the Milky Way and Magellanic Clouds. We present an infrared (IR) spectropolarimetric detection of a CCSN, enabled by the new highly sensitive WIRC+Pol instrument at Palomar Observatory, that can observe CCSNe (M = -17 mags) out to 20 Mpc to ~0.1% polarimetric precision. IR spectropolarimetry is less affected than optical by dust scattering in the circumstellar and interstellar media, thereby providing a more unbiased probe of the intrinsic geometry of the SN ejecta. SN 2018hna, a SN 1987A-like explosion, shows 2.0+-0.3% continuum polarization in the J band oriented at ~160 degree on-sky at 182 d after the explosion. Assuming prolate geometry like in SN 1987A, we infer an ejecta axis ratio of <0.48 with the axis of symmetry pointing at 70 degree position angle. The axis ratio is similar to that of SN 1987A suggesting that they may share intrinsic geometry and inclination angle. Our data do not rule out oblate ejecta. We also observe one other core-collapse and two thermonuclear SNe in the J band. SN 2020oi, a stripped-envelope Type Ic SN in Messier 100 has p = 0.37+-0.09% at peak light, indicative of either a 10% asymmetry or host interstellar polarization. The SNe Ia, 2019ein and 2020ue have p < 0.33% and < 1.08% near peak light, indicative of asymmetries of less than 10% and 20%, respectively.

The rise of multi-messenger astronomy has brought with it the need to exploit all available data streams and learn more about the astrophysical objects that fall within its breadth. One possible avenue is the search for serendipitous optical/near-infrared counterparts of gamma-ray bursts (GRBs) and gravitational-wave (GW) signals, known as kilonovae. With surveys such as the Zwicky Transient Facility (ZTF), which observes the sky with a cadence of ~ three days, the existing counterpart locations are likely to be observed; however, due to the significant amount of sky to explore, it is difficult to search for these fast-evolving candidates. Thus, it is beneficial to optimize the survey cadence for realtime kilonova identification and enable further photometric and spectroscopic observations. We explore how the cadence of wide field-of-view surveys like ZTF can be improved to facilitate such identifications. We show that with improved observational choices, e.g., the adoption of three epochs per night on a ~ nightly basis, and the prioritization of redder photometric bands, detection efficiencies improve by about a factor of two relative to the nominal cadence. We also provide realistic hypothetical constraints on the kilonova rate as a form of comparison between strategies, assuming that no kilonovae are detected throughout the long-term execution of the respective observing plan. These results demonstrate how an optimal use of ZTF increases the likelihood of kilonova discovery independent of GWs or GRBs, thereby allowing for a sensitive search with less interruption of its nominal cadence through Target of Opportunity programs., Comment: 10 pages, 8 figures

We study a sample of 16 Type Ia supernovae (SNe Ia) having both spectroscopic and photometric observations within 2 $-$ 3 days after the first light. The early $B-V$ colors of such a sample tends to show a continuous distribution. For objects with normal ejecta velocity (NV), the C~II $\lambda$6580 feature is always visible in the early spectra while it is absent or very weak in the high-velocity (HV) counterpart. Moreover, the velocities of the detached high-velocity features (HVFs) of Ca~II NIR triplet (CaIR3) above the photosphere are found to be much higher in HV objects than in NV objects, with typical values exceeding 30,000 km~s$^{-1}$ at 2 $-$ 3 days. We further analyze the relation between %velocities of Si~II~$\lambda$6355 at maximum, $v_{\rm Si,max}$, the velocity shift of late-time [Fe~II] lines ($v_{\rm [Fe~II]}$) and host galaxy mass. We find that all HV objects have redshifted $v_{\rm [Fe~II]}$ while NV objects have both blue- and redshifted $v_{\rm [Fe~II]}$. It is interesting to point out that the objects with redshifted $v_{\rm [Fe~II]}$ are all located in massive galaxies, implying that HV and a portion of NV objects may have similar progenitor metallicities and explosion mechanisms. We propose that, with a geometric/projected effect, the He-detonation model may account for the similarity in birthplace environment and the differences seen in some SNe Ia, including $B-V$ colors, C~II feature, CaIR3 HVFs at early time and $v_{\rm [Fe~II]}$ in the nebular phase. Nevertheless, some features predicted by He-detonation simulation, such as the rapidly decreasing light curve, deviate from the observations, and some NV objects with blueshifted nebular $v_{\rm [Fe~II]}$ may involve other explosion mechanisms., Comment: accepted by the ApJ

Neutron star mergers eject neutron-rich matter in which heavy elements are synthesised. The decay of these freshly synthesised 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 $Y_\mathrm{e} < 0.25$), they contain fair amounts of lanthanides and actinides which 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., Comment: 13 pages, 6 figures, 1 table, accepted by MNRAS

LIGO and Virgo's third observing run (O3) revealed the first neutron star-black hole (NSBH) merger candidates in gravitational waves. These events are predicted to synthesize r-process elements creating optical/near-IR "kilonova" (KN) emission. The joint gravitational-wave (GW) and electromagnetic detection of an NSBH merger could be used to constrain the equation of state of dense nuclear matter, and independently measure the local expansion rate of the universe. Here, we present the optical follow-up and analysis of two of the only three high-significance NSBH merger candidates detected to date, S200105ae and S200115j, with the Zwicky Transient Facility (ZTF). ZTF observed $\sim$\,48\% of S200105ae and $\sim$\,22\% of S200115j's localization probabilities, with observations sensitive to KNe brighter than $-$17.5\,mag fading at 0.5\,mag/day in g- and r-bands; extensive searches and systematic follow-up of candidates did not yield a viable counterpart. We present state-of-the-art KN models tailored to NSBH systems that place constraints on the ejecta properties of these NSBH mergers. We show that with depths of $\rm m_{\rm AB}\approx 22$ mag, attainable in meter-class, wide field-of-view survey instruments, strong constraints on ejecta mass are possible, with the potential to rule out low mass ratios, high BH spins, and large neutron star radii.

Spectroscopic detection of narrow emission lines traces the presence of circumstellar mass distributions around massive stars exploding as core-collapse supernovae. Transient emission lines disappearing shortly after the supernova explosion suggest that the spatial extent of such material is compact, and hence imply an increased mass loss shortly prior to explosion. Here, we present a systematic survey for such transient emission lines (Flash Spectroscopy) among Type II supernovae detected in the first year of the Zwicky Transient Facility (ZTF) survey. We find that at least six out of ten events for which a spectrum was obtained within two days of estimated explosion time show evidence for such transient flash lines. Our measured flash event fraction ($>30\%$ at $95\%$ confidence level) indicates that elevated mass loss is a common process occurring in massive stars that are about to explode as supernovae.

Kilonovae produced by the coalescence of compact binaries with at least one neutron star are promising standard sirens for an independent measurement of the Hubble constant ($H_0$). Through their detection via follow-up of gravitational-wave (GW), short gamma-ray bursts (sGRBs) or optical surveys, a large sample of kilonovae (even without GW data) can be used for $H_0$ contraints. Here, we show measurement of $H_0$ using light curves associated with four sGRBs, assuming these are attributable to kilonovae, combined with GW170817. Including a systematic uncertainty on the models that is as large as the statistical ones, we find $H_0 = 73.8^{+6.3}_{-5.8}$\,$\mathrm{km}$ $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$ and $H_0 = 71.2^{+3.2}_{-3.1}$\,$\mathrm{km}$ $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$ for two different kilonova models that are consistent with the local and inverse-distance ladder measurements. For a given model, this measurement is about a factor of 2-3 more precise than the standard-siren measurement for GW170817 using only GWs., Comment: Updating GRB citations