Hintze Centre for Astrophysical Surveys

Abstract:

Kilonovae are the scientifically rich—but observationally elusive—optical transient phenomena associated with compact binary mergers. Only a handful of events have been discovered to date, all through multiwavelength (gamma-ray) and multimessenger (gravitational-wave) signals. Given their scarcity, it is important to maximise the discovery possibility of new kilonova events. To this end, we present our follow-up observations of the gravitational-wave signal S250818k—a plausible binary neutron star merger at a distance of 237 ± 62 Mpc. Pan-STARRS tiled 286 and 318 deg2 (32% and 34% of the 90% sky localisation region) within 3 and 7 days of the GW signal, respectively. ATLAS covered 65% of the sky map within 3 days, but with lower sensitivity. These observations uncovered 47 new transients; however, none were deemed to be linked to S250818k. We undertook an expansive follow-up campaign of AT2025ulz, the purported counterpart to S250818k. The griz-band light curve, combined with our redshift measurement (z = 0.0849 ± 0.0003), all indicate that SN2025ulz is a type IIb supernova and thus not the counterpart to S250818k. We rule out the presence of an AT2017gfo-like kilonova within ≈27% of the distance posterior sampled by our Pan-STARRS pointings (≈9.1% across the total 90% 3D sky localisation). We demonstrate that early observations are optimal for probing the distance posterior of the 3D gravitational-wave sky map, and that SN2025ulz was a plausible kilonova candidate for ≲5 days, before ultimately being ruled out.

Abstract:

We present optical observations of the type IIb supernova (SN) 2017ckj, covering approximately 180 days after the explosion. Its early-time multi-band light curves display no clear evidence of a shock-cooling tail, resembling the behaviour of SN 2008ax. The V -band light curve exhibits a short rise time of about 5 days and reaches an absolute fitted peak magnitude of M V  = −18.49 ± 0.18 mag. The late-time multi-band light curves reveal a linear decline. We modelled the bolometric light curve of SN 2017ckj to constrain the progenitor and the explosion parameters. We estimated a total mass of 56 Ni synthesised by SN 2017ckj of M Ni = 0.21 +0.05 −0.03 M ⊙ , with a massive H-rich envelope of M env = 0.4 +0.1 −0.1 M ⊙ . Both the 56 Ni mass and the envelope mass of SN 2017ckj are higher than those of typical SNe IIb, in agreement with its peculiar light curve evolution. The early-time spectra of SN 2017ckj are dominated by a blue continuum, accompanied by narrow H α and He  II emission lines. The earliest spectrum exhibits flash ionisation features, from which we estimated a progenitor mass-loss rate of ∼3 × 10 −4  M ⊙  yr −1 . At later epochs, the spectra develop broad P-Cygni profiles and become increasingly similar to those of SNe IIb, especially SN 2018gk. The late-time spectrum at around 139 days does not show a distinct decline in the strength of its H α emission profile, also indicating a relatively massive envelope of its progenitor. Aside from the H α feature, the nebular spectrum exhibits prominent emission lines of O  I , Ca  II , [Ca  II ], and Mg  I ], which are consistent with the prototypical SN 1993J.

Abstract:

Context. The H I mass function (HIMF) is a crucial tool for understanding the evolution of the H I content in galaxies over cosmic time and, hence, to constraining both the baryon cycle in galaxy evolution and the reionization history of the Universe. Aims. We aim to derive semiempirical constraints at z  ∼ 0.37 by combining literature results on the stellar mass function from optical surveys with recent findings on the M HI  −  M ⋆ scaling relation derived via spectral stacking analysis applied to 21 cm line interferometric data from the MIGHTEE and CHILES surveys, conducted with the MeerKAT and VLA radio telescopes, respectively. Methods. We drew synthetic stellar mass samples directly from the publicly available results underlying the analysis of the COSMOS2020 galaxy photometric sample. We then converted M ⋆ into M HI using analytical fitting functions to the data points from H I stacking. We next fit a Schechter function to the median HIMF from all the samples via Monte Carlo Markov chains. We finally derived the posterior distribution for Ω HI by integrating the models for the HIMF built from the posteriors samples of the Schechter parameters. Results. We find a deviation of the HIMF at z  ∼ 0.37 from the results at z  ∼ 0 from the ALFALFA survey and at z  ∼ 1 from uGMRT data. Our results for Ω HI are in broad agreement with other literature results and follow the overall trend on Ω HI as a function of redshift. The derived value Ω HI = (7.02 +0.59 −0.52 ) × 10 −4 at z  ∼ 0.37 from the combined analysis deviates by ∼2.9 σ from the ALFALFA result at z  ∼ 0. Conclusions. Our findings regarding the HIMF and Ω HI derived from deep, state-of-the-art interferometric surveys differ from previous literature results at z  ∼ 0 and z  ∼ 1. We are unable to confirm at this stage whether these differences are due to cosmic evolution consistent with a smooth transition of the H I content of galaxies over the last 8 Gyr or due to selection biases and systematics.