The Square Kilometre Array (SKA)

Abstract:

Abstract Synchrotron self absorption (SSA) is seen across a variety of astrophysical sources, and observation of an SSA peak in the spectrum is a powerful tool for estimating the physical conditions and the minimum energy of the emitting region. We begin with the (re)derivation of the usual SSA parameter estimates, carefully considering dependencies and assumptions, obtaining the most accurate traditional SSA minimum energy equations currently available. Traditional methods rely on the assumption that the emitting region is quasi-spherical and homogeneous. However, many observations of SSA show that the spectral index at frequencies below the peak is less than the expected +2.5 (non-thermal) or +2 (thermal). We argue that an inhomogeneous emitting region is the most likely explanation in many cases. Power law inhomogeneous cylindrical slab and broken power law inhomogeneous sphere models are used to investigate how the presence of inhomogeneity affects parameter estimates using traditional SSA methods. We find that in some cases inhomogeneity can lead to traditional SSA methods underestimating the minimum energy and the size of the emitting region by over an order of magnitude. Quantitative correction factors are found which can be applied to traditional estimates to correct for inhomogeneity, depending on the value of the observed flattened spectral index and the range in frequency over which this value is observed. Furthermore, we derive simple correction factors for non-spherical homogeneous emitting regions. Finally, we explore the effects of inhomogeneity on measurements of polarisation around the spectral peak, and on lightcurves for expanding emitting regions.

Abstract:

The detection of gravitational waves (GWs) from a binary neutron star (BNS) merger by Advanced LIGO and Advanced Virgo (GW170817), together with its electromagnetic counterpart, the short gamma-ray burst GRB 170817A, heralded the birth of multimessenger astronomy. The detection of TeV emission from GRBs motivates follow-up observations with the Cherenkov Telescope Array Observatory (CTAO), which is ideal for detecting such signals due to its unprecedented sensitivity, rapid response, and wide-field survey capabilities. The aim of this work is to evaluate GeV–TeV GW follow-up strategies for CTAO using a multistep simulation pipeline and to estimate the expected rate of joint GW–GRB detections during observing run O5. Using a simulated sample of BNS systems with corresponding GW detections, gamma-ray emission is simulated through phenomenological prescriptions based on the observed population of short GRBs, including off-axis jet scenarios. CTAO observations are simulated to account for instrument response, sky tiling strategies, integration times, and varying observing conditions. Strategies with variable and constant integration times are investigated. We find that, via an optimized follow-up strategy, about 5% of simulated GW-associated short GRBs produce GeV–TeV radiation detectable by CTAO. Detectability is strongly influenced by the jet opening angle and viewing angle, suggesting that even rough estimates of the viewing angle in GW alerts could enhance targeting. This framework motivates future follow-ups of GW-detectable events, including neutron star–black hole mergers, and further supports the development of advanced strategies incorporating galaxy distributions and synergies with future detectors such as the Einstein Telescope.

Abstract:

We present the Thousand Pulsar Array (TPA) single-pulse data release, obtained with the MeerKAT radio telescope and comprising time-series observations of 1192 pulsars, typically containing consecutive pulses per source. We describe the MeerTime Single Pulse software pipeline which calibrates the data and automatically excises interference signals to produce data products suitable for typical single-pulse studies. To demonstrate the capabilities of the data set, we carry out a population-level study of phase-averaged single-pulse energy distributions and nulling behaviour. Pulse energy distributions are modelled within a Bayesian framework choosing from a range of intrinsic energy distributions, and including an explicit nulling fraction. We find that approximately half of the pulsars require multicomponent intrinsic energy distributions, while the remainder are consistent with single-component models. Nulling is detected or constrained for most pulsars in the sample, and both the occurrence and inferred nulling fraction show systematic variation across the P– diagram. In particular, nulling fractions increase with spin period and exhibit only a weak dependence on period derivative. We also examine trends in the preferred forms of pulse energy distributions as a function of spin-down luminosity, finding modest evidence for population-level evolution. Estimates of single-pulse luminosities indicate that individual pulses can exceed the long-term average luminosity by large factors, particularly for low- pulsars. These results characterize the statistical properties of single-pulse emission across a large pulsar sample and highlight the limitations of phase-averaged energy distributions for capturing the full complexity of pulsar emission variability.