Dynamic shocks powered by a wide, relativistic, super-Eddington outflow launched by an accreting neutron star in the mid-20th century

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

Authors:

FJ Cowie, RP Fender, I Heywood, F Carotenuto, JH Matthews, B Reville, L Olivera-Nieto, AJ Cooper, AK Hughes, K Savard, PA Woudt, J van den Eijnden, N Grollimund, P Saikia

Abstract:

Abstract Accreting systems can launch powerful outflows which interact with the surrounding medium. We combine new radio observations of the accreting neutron star X-ray binary (XRB) Circinus X-1 (Cir X-1) with archival radio observations going back 24 years. The ∼3 pc scale wide-angle radio and X-ray emitting caps found around Cir X-1 are identified as synchrotron emitting shocks with significant proper motion and morphological evolution on decade timescales. Proper motion measurements of the shocks reveal they are mildly relativistic and decelerating, with apparent velocity of 0.14c ± 0.03c at a propagation distance of 2 pc. We demonstrate that these shocks are likely powered by a hidden relativistic (≳ 0.3c) wide-angle conical outflow launched in 1972 ± 3, in stark contrast to known structures around other XRBs formed by collimated jets over 1000s of years. The minimum time-averaged power of the outflow required to produce the observed synchrotron emission is ∼0.1LEdd, while the time-averaged power required for the kinetic energy of the shocks is $\sim 40 \left(\frac{n}{10^{-2} \textrm{cm}^{-3}}\right)L_\textrm{Edd}$, where n is the average ambient medium number density. This reveals the outflow powering the shocks is likely significantly super-Eddington. We measure significant linear polarisation up to 52 ± 6% in the shocks demonstrating the presence of an ordered magnetic field of strength ∼200 μG. We show that the shocks are potential PeVatrons, capable of accelerating electrons to ∼0.7 PeV and protons to ∼20 PeV, and we estimate the injection and energetic efficiencies of electron acceleration in the shocks. Finally, we predict that next generation gamma-ray facilities may be able to detect hadronic signatures from the shocks.

Simulating radio emission from flickering AGN jets: travelling shocks and hotspot brightening

Monthly Notices of the Royal Astronomical Society Oxford University Press (OUP) 546:2 (2026) stag131

Authors:

Emma L Elley, James H Matthews, Dipanjan Mukherjee, Bhargav Vaidya

Abstract:

ABSTRACT We investigate the impact of flickering variability in jet power on the luminosity and morphology of radio galaxies. We use a Lagrangian particle method together with relativistic hydrodynamics simulations using the pluto code to track the evolution of electron spectra through particle acceleration at shocks and cooling processes. We introduce an adapted version of this method which improves tracking of adiabatic cooling in regimes where low density jet material mixes with high density from the ambient medium in the lobes. We find that rapid increases in jet power can lead to large increases in hotspot luminosity due to the interaction of a travelling shock structure with the pre-existing shock structure at the jet head. We show that in some cases it may be possible to identify a bright region of emission corresponding to a shock travelling along the jet axis. We find that the time-averaged radiative efficiency of variable jets is similar to their steady counterparts, but find significant departures from this on an instantaneous basis. We suggest that, together with environmental effects and differences in the average powers of jets, variable jet powers may have a significant impact on how we understand the diversity of radio jets seen in observations and have significant implications for interpretations of jet powers, energy budgets, and luminosity-linear size diagrams.

Calibrating Mid-Infrared Emission Features As Diagnostics of Star Formation in Infrared-Luminous Galaxies via Radiative Transfer Modeling

(2026)

Authors:

L Robinson, D Farrah, A Efstathiou, A Engholm, E Hatziminaoglou, M Joyce, V Lebouteiller, S Petty, LK Pitchford, J Afonso, D Clements, M Lacy, C Pearson, D Rigopoulou, M Rowan-Robinson, L Wang

Cosmic rays, gamma rays and neutrinos from discrete black hole X-ray binary ejecta

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

Authors:

Nicolas J Bacon, Alex J Cooper, Dimitrios Kantzas, James H Matthews, Rob Fender

Abstract:

Abstract The origin of cosmic rays from outside the Solar system are unknown, as they are deflected by the interstellar magnetic field. Supernova remnants are the main candidate for cosmic rays up to PeV energies but due to lack of evidence, they cannot be concluded as the sources of the most energetic Galactic CRs. We investigate discrete ejecta produced in state transitions of black hole X-ray binary systems as a potential source of cosmic rays, motivated by recent >100 TeV γ-ray detections by LHAASO. Starting from MAXI J1820+070, we examine the multi-wavelength observations and find that efficient particle acceleration may take place (i.e. into a robust power-law), up to ∼2 × 1016μ−1/2 eV, where μ is the ratio of particle energy to magnetic energy. From these calculations, we estimate the global contribution of ejecta to the entire Galactic spectrum to be $\sim 1~{{\ \rm per\ cent}}$, with the cosmic ray contribution rising to $\sim 5~{{\ \rm per\ cent}}$ at PeV energies, assuming roughly equal energy in non-thermal protons, non-thermal electrons and magnetic fields. In addition, we calculate associated γ-ray and neutrino spectra of the MAXI J1820+070 ejecta to investigate new detection methods with CTAO, which provide strong constraints on initial ejecta size of order 107 Schwarzschild radii (10−5 pc) assuming a period of adiabatic expansion.

Exploring the quasar disc-wind-jet connection with LoTSS and SDSS

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

Authors:

Charlotte L Jackson, James H Matthews, Imogen H Whittam, Matt J Jarvis, Matthew J Temple, Amy L Rankine, Paul C Hewett

Abstract:

Abstract We investigate the relationship between disc winds, radio jets, accretion rates and black hole masses of a sample of ∼100k quasars at z ≈ 2. Combining spectra from the 17th data release of the Sloan Digital Sky Survey (SDSS) with radio fluxes from the 2nd data release of the Low Frequency ARray (LOFAR) Two-Meter Sky Survey (LoTSS), we statistically characterise a radio loud and radio quiet population using a two-component Gaussian Mixture model, and perform population matching in black hole mass and Eddington fraction. We determine how the fraction of radio loud sources changes across this parameter space, finding that jets are most efficiently produced in quasars with either a very massive central black hole (MBH > 109M⊙) or one that is rapidly accreting (λEdd > 0.3). We also show that there are differences in the blueshift of the $\textrm {C}\, \rm \small {IV}$ λ1549Å line and the equivalent width of the $\rm {He}\, \rm \small {II}$ λ1640Å line in radio loud and radio quiet quasars that persist even after accounting for differences in the mass and accretion rate of the central black hole. Generally, we find an anti-correlation between the inferred presence of disc winds and jets, which we suggest is mediated by differences in the quasars’ spectral energy distributions. The latter result is shown through the close coupling between tracers of wind kinematics and the ionising flux– which holds for both radio loud and radio quiet sources, despite differences between their emission line properties– and is hinted at by a different Baldwin effect in the two populations.