Pulsars, transients and relativistic astrophysics

Abstract:

Abstract We compute and test analytic models for the plunging region dynamics, thermodynamics, and magnetic fields against dedicated 3D global general relativistic magnetohydrodynamics (MHD) simulations of thin accretion discs around black holes across the spin range, using the code ATHENAK. We find that the dynamics of the plunging fluid closely resembles that of a gravity-dominated geodesic plunge, with the best agreement at low spins. Additionally, we find good agreement between the thermodynamic framework and the simulated quantities across the spin range. Finally, we develop a new model for the magnetic fields in the plunging region that assumes a fixed geodesic inflow, into which the magnetic fields are frozen. Overall, our simulations are in good concordance with this model, albeit with some discrepancies that suggest a degree of non-ideal MHD dissipation. In addition, we investigate how the MHD stresses in the plunging region depend on the black hole spin, interpreting our results through the lens of our flux-freezing model. We find that the magnitude of the stress increases as the black hole spin is increased in the prograde direction. This question is of particular importance for observers who wish to determine the black hole spin from X-ray measurements of the inner accretion disc, since a low-stress, high-spin solution is degenerate with a high-stress, low-spin solution. The spin-stress relationship that we report is approximately orthogonal to the contour of degenerate spin-stress pairings, indicating that the degeneracy is not fundamental. We show this explicitly for the case of M33 X-7.

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.

Abstract:

The propagation of very-high-energy (VHE; E γ ≥ 100 GeV) gamma rays from extragalactic sources is affected by interactions with photons of the extragalactic background light (EBL), resulting in pair production that attenuates the intrinsic gamma-ray flux. This interaction renders the Universe increasingly opaque to VHE photons at high energies and redshifts. New physics scenarios involving axion-like particles (ALPs) could modify this expected optical depth. In particular, ALPs with masses ma ∼ 10 eV can decay into two photons over cosmological timescales, thereby contributing to the diffuse EBL. If such ALPs constitute a significant fraction of the dark matter density, their decay would enhance the EBL intensity and consequently increase the gamma-ray optical depth. In this study, we investigate this scenario using a large sample of gamma-ray spectra observed with the High Energy Stereoscopic System (H.E.S.S.) and the Fermi Large Area Telescope. We model the contribution of decaying ALPs to the EBL and assess their impact on the spectra of blazars across redshifts. By comparing these observations with standard EBL models, we place constraints on the properties of heavy ALPs — specifically their mass and photon coupling — and evaluate their viability as a dark matter candidate capable of modifying the gamma-ray transparency of the Universe. From the combined analysis, and under the assumption that ALPs constitute the entire dark matter density, we derive 95% confidence exclusion limits on the photon-ALP coupling down to maγ ∼ 7×10-12 GeV-1 for masses ma ∼ 15 eV. These constraints are competitive with existing astrophysical bounds and provide complementary sensitivity to other techniques, closing a previously unconstrained region of parameter space in the ma ∼ 2.5–20 eV range.

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.