Calibration-induced Systematics in SALT3 Training and Their Impact on Dark Energy Constraints from Stage IV Supernova Surveys

The Astrophysical Journal American Astronomical Society 1010:1 (2026) 50

Authors:

Kene Anumba, David O Jones, Richard Kessler, Daniel Scolnic, W D’Arcy Kenworthy, Rebecca C Chen, Bastien Carreres, Maria Vincenzi, Erik R Peterson, Maria Acevedo, Ben Rose, Dillon Brout, Jillian Paulin, Rujuta A Purohit, Rebekah Hounsell, The Roman Supernova Cosmology Project Infrastructure Team

Abstract:

In the coming years, the Vera Rubin Observatory’s Legacy Survey of Space and Time (Rubin-LSST) and the Nancy Grace Roman Space Telescope’s (Roman) High Latitude Time Domain Survey (HLTDS) are expected to discover more than a million Type Ia supernovae (SNe Ia), several orders of magnitude more than current samples and with a tighter control on systematic uncertainties. One of the largest systematic uncertainties in cosmological analyses with SNe Ia is the accuracy of the spectrophotometric model for SNe Ia time series data, which depends on the photometric calibration of the surveys. To quantify the impact of this uncertainty, we analyze simulated Rubin-LSST and HLTDS data, perturb the photometric zero-points and filter mean wavelengths, and propagate these systematics to spectral model recovery, estimated distances, and dark energy Figure of Merit (FoM) based on the w0waCDM model. Zero-point shifts of 5 mmag and filter mean wavelength shifts of 5 Å lead to a ∼50% decrease in the FoM relative to a statistical-only case when calibration uncertainties are propagated only through light-curve fitting. The same calibration shifts applied only during model training produce a smaller ∼13% degradation. Contrary to previous analyses, calibration uncertainties in light-curve fitting dominate over those from model training. Their effect during light-curve fitting varies smoothly with redshift and is nearly degenerate with cosmology, preventing mitigation through self-calibration. Finally, we show that the FoM dependence on the size of the calibration uncertainties (in the range of expected sizes) is roughly linear.

Identification of a Radio Counterpart to SN 2025ulz in the S250818k Localization Area

The Astrophysical Journal Letters American Astronomical Society 1009:2 (2026) L44-L44

Authors:

Tanner O’Dwyer, Alessandra Corsi, Deepika Yadav, Kunal P Mooley, Raphael Baer-Way, Poonam Chandra, Gregg Hallinan, Mansi M Kasliwal, Lauren Rhodes, Oleg M Smirnov, Davide Lazzati, Joeri van Leeuwen, Adam Deller, Pikky Atri, Tanazza Khanam

Abstract:

Abstract On 2025 August 18, the LIGO – Virgo – KAGRA collaboration reported S250818k, a subthreshold gravitational-wave (GW) candidate consistent with a binary neutron star (NS) merger potentially involving a subsolar-mass NS. Optical follow-up by the Zwicky Transient Facility identified a transient temporally coinciding with the GW trigger that initially resembled a kilonova but was later classified as a young stripped-envelope Type IIb supernova (SN), dubbed SN 2025ulz. A key question is whether SN 2025ulz harbors fast, possibly collimated, nonthermal ejecta indicative of a central engine, as invoked in “superkilonova” scenarios linking subsolar-mass NSs to accretion-disk fragmentation or core fission. We present early-to-late-time multiband radio observations of SN 2025ulz obtained with the Karl G. Jansky Very Large Array as part of the JAGWAR program, complemented by observations with the upgraded Giant Metrewave Radio Telescope and MeerKAT. We detect a faint but significant radio counterpart to SN 2025ulz at 6−10 GHz. The data are consistent with nonthermal emission from SN ejecta interacting with circumstellar material. Our data are also consistent with emission from an off-axis jet peaking at ∼50−100 days after the GW trigger. While we cannot distinguish between these two interpretations, our radio detection motivates future systematic multiwavelength follow-up of core-collapse events coinciding with subsolar NS GW candidates, should the association between S250818k and SN 2025ulz be supported by offline GW analyses.

The impact of flickering variability and magnetisation on the dynamics, stability and morphology of radio-loud AGN jets

Monthly Notices of the Royal Astronomical Society Oxford University Press (OUP) (2026) stag1798

Authors:

Emma L Elley, James H Matthews, Henry Whitehead, Alex J Cooper

Abstract:

Abstract The physics governing the morphology of radio-loud AGN jets is not fully understood. We investigate how magnetization, flickering jet power and their interplay affects the morphology of radio galaxies. We present a grid of relativistic magnetohydrodynamic simulations using the PLUTO code covering constant and variable jets with two levels of magnetisation. We find that the constant high magnetisation jets can lead to highly asymmetrical cocoon morphologies, whilst the variable high magnetisation jet can exhibit a broken morphology, caused by a discontinuous jet beam. Our work highlights the importance of magnetisation and variability on the stability and resulting morphology of radio-loud AGN jets, suggesting both are significant factors in addition to jet power or environment. Furthermore, we show that the interaction between magnetisation and variability can lead to the development of localised kink instabilities along the jet beam. Finally, we discuss the effects of hydrodynamic mixing in low magnetisation jets and the role of viewing angle dependence in comparisons between our simulations and observed sources. To facilitate this comparison we present a library of simulated radio images at different times in the simulations and from various viewing angles, which highlight a diverse set of complex morphologies.

Eccentric stellar-mass binary black holes: Population, detectability, and waveform analysis in the LISA and LIGO era

Physical Review D American Physical Society (APS) 114:6 (2026) 063036

Authors:

Zeyuan Xuan, Smadar Naoz, Kyle Kremer, Michael L Katz, Bence Kocsis, Erez Michaely

Abstract:

Eccentric binary black holes (BBHs) formed through dynamical interactions can significantly contribute to gravitational wave (GW) detections. In this work, we present a simulated catalog of dynamically formed, stellar-mass BBHs in the local universe, incorporating contributions from the Galactic field (flyby interactions), Galactic nucleus (eccentric Kozai-Lidov evolution), and globular clusters (N-body interactions). Our results predict a wide, highly eccentric BBH population in the Milky Way, with source counts of ∼25, 9, 3, 1, and <1 for a 4-year LISA mission (∼36, 13, 5, 2, and 1 for a 10-year extended mission), for SNR>1, 3, 8, 20, and 50, respectively. Extending this model to cosmological populations, we show that different dynamical channels can produce distinct eccentricity distributions in the LVK band and can contribute hundreds of additional low-SNR mHz sources. Specifically, our model yields a merger rate of Γ∼9 Gpc-3 yr-1 and ∼150 extragalactic mHz BBHs with SNR>1. However, due to the lower mass and weaker GW signals of stellar-mass BBHs, this number declines sharply at higher detection thresholds (e.g., ∼1 for SNR>8). We further highlight the impact of eccentric BBH signals on the LISA global fit, showing that their individual harmonics can be independently detected in the Milky Way, and may mimic circular binaries with systematically biased chirp masses. Lastly, we show that post-Newtonian waveforms converge reliably for eccentric BBHs with masses of ≲103M⊙ in the mHz band. Overall, eccentric BBHs represent a prevalent and promising target for future space-based GW observatories.

cuDART: a GPU-accelerated ray tracing code for generating synthetic observations of relativistic astrophysical sources

ArXiv 2609.09322 (2026)

Authors:

Henry Whitehead, Emma L Elley, Christoper N Everett