Pan-STARRS Follow-up of the Gravitational-wave Event S250818k and the Light Curve of SN2025ulz
The Astrophysical Journal Letters American Astronomical Society 995:1 (2025) L27
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.Upper limits on radio emission from the K2-18 system
Monthly Notices of the Royal Astronomical Society 91̽»¨ University Press 546:3 (2025) staf1998
Abstract:
Stellar and planetary magnetic fields play a crucial role in the habitability of a planet and the integrity of its atmosphere. The detection of methane and carbon dioxide, along with a tentative identification of the potential biosignature dimethyl sulfide/disulfide, in the atmosphere of K2-18 b, a sub-Neptune orbiting an M dwarf star, presents an intriguing question regarding the stellar magnetic environment and the resistance of the planet’s magnetosphere (if it exists) to erosion by magnetic activity from the host. To probe for radio emission from the system, we have conducted observations using the Karl G. Jansky Very Large Array at S, C, and X bands (2–4, 5.5–7.5, and 8–10 GHz, respectively) to search for coherent and incoherent radio emission. We detect no radio emission associated with incoherent emission mechanisms. We report Stokes I upper limits of at S band, at C band, and at X band and an upper limit of the ratio of the radio to the total bolometric luminosity of . We have also searched for short duration bursts associated with coherent emission mechanisms at C and X bands. No signals above a significance threshold are detected. Although no signals are detected, our radio observations offer constraints, albeit limited, on the stellar magnetic environment 91̽»¨ing recent X-ray observations indicating that K2-18 is a very faint emitter. Our results also contextualize any planetary transmission spectra by providing constraints on the activity level of the host.The Radio Flare and Multiwavelength Afterglow of the Short GRB 231117A: Energy Injection from a Violent Shell Collision
The Astrophysical Journal American Astronomical Society 994:1 (2025) 5
Abstract:
We present the early radio detection and multiwavelength modeling of the short gamma-ray burst (GRB) 231117A at redshift z = 0.257. The Australia Telescope Compact Array automatically triggered a 9 hr observation of GRB 231117A at 5.5 and 9 GHz following its detection by the Neil Gehrels Swift Observatory just 1.3 hr post-burst. Splitting this observation into 1 hr time bins, the early radio afterglow exhibited flaring, scintillating and plateau phases. The scintillation allowed us to place the earliest upper limit (<10 hr) on the size of a GRB blast wave to date, constraining it to <1 × 1016 cm. Multiwavelength modeling of the full afterglow required a period of significant energy injection between ∼0.02 and 1 day. The energy injection was modeled as a violent collision of two shells: a reverse shock passing through the injection shell explains the early radio plateau, while an X-ray flare is consistent with a shock passing through the leading impulsive shell. Beyond 1 day, the blast wave evolves as a classic decelerating forward shock with an electron distribution index of p = 1.66 ± 0.01. Our model also indicates a jet break at ∼2 days, and a half-opening angle of θj=16.°6±1.°1 . Following the period of injection, the total energy is ζ ∼ 18 times the initial impulsive energy, with a final collimation-corrected energy of EKf ∼ 5.7 × 1049 erg. The minimum Lorentz factors this model requires are consistent with constraints from the early radio measurements of Γ > 35 to Γ > 5 between ∼0.1 and 1 day. These results demonstrate the importance of rapid and sensitive radio follow-up of GRBs for exploring their central engines and outflow behaviour.New Metrics for Identifying Variables and Transients in Large Astronomical Surveys
The Astrophysical Journal American Astronomical Society 992:1 (2025) 109
Abstract:
A key science goal of large sky surveys such as those conducted by the Vera C. Rubin Observatory and precursors to the Square Kilometre Array is the identification of variable and transient objects. One approach is analyzing time series of the changing brightness of sources, namely, light curves. However, finding adequate statistical representations of light curves is challenging because of the sparsity of observations, irregular sampling, and nuisance factors inherent in astronomical data collection. The wide diversity of objects that a large-scale survey will observe also means that making parametric assumptions about the shape of light curves is problematic. We present a Gaussian process (GP) regression approach for characterizing light-curve variability that addresses these challenges. Our approach makes no assumptions about the shape of a light curve and, therefore, is general enough to detect a range of variable and transient source types. In particular, we propose using the joint distribution of GP amplitude hyperparameters to distinguish variable and transient candidates from nominally stable ones and apply this approach to 6394 radio light curves from the ThunderKAT survey. We compare our results with two variability metrics commonly used in radio astronomy, namely ην and Vν, and show that our approach has better discriminatory power and interpretability. Finally, we conduct a rudimentary search for transient sources in the ThunderKAT data set to demonstrate how our approach might be used as an initial screening tool. Computational notebooks in Python and R are available to help deploy this framework to other surveys.Radio emission from a nearby M dwarf binary
Monthly Notices of the Royal Astronomical Society 91̽»¨ University Press 543:3 (2025) 1935-1944