Gamma-ray astronomy

Abstract:

DIPLODOCUS (Distribution-In-PLateaux methODOlogy for the CompUtation of transport equationS) is a framework being developed for the general transport of particle distribution functions through the seven dimensions of phase space, including forcing terms and interactions between particles. Following Paper I, which details the mathematical background, this second paper provides an overview of the numerical implementation in the form of the code package Diplodocus . jl, written in Julia, including the description of a novel Monte-Carlo sampling technique for the pre-computation of anisotropic collision integrals. In addition to the discussion of numerical implementation, a selection of test cases are presented to examine the package’s capabilities. These test cases focus on micro-scale physical effects: binary collisions, emissive interactions and external forces that are relevant to the modelling of jetted astrophysical sources, such as Active Galactic Nuclei and X-Ray Binaries.

Abstract:

DIPLODOCUS (Distribution-In-PLateaux methODOlogy for the CompUtation of transport equationS) is a novel framework being developed for the mesoscopic modelling of astrophysical systems via the transport of particle distribution functions through the seven dimensions of phase space, including continuous forces and discrete interactions between particles. This first paper in a series provides an overview of the analytical framework behind the model, consisting of an integral formulation of the relativistic transport equations (Boltzmann equations) and a discretisation procedure for the particle distribution function (Distribution-In-Plateaux). The latter allows for the evaluation of anisotropic interactions, and generates a conservative numerical scheme for a distribution function’s transport through phase space.

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.

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.

Abstract:

Abstract We present the longest and the densest quasi-simultaneous radio, X-ray and optical campaign of the black hole low mass X-ray binary GX 339–4, covering five years of weekly GX 339–4 monitoring with MeerKAT, Swift/XRT and MeerLICHT, respectively. Complementary high frequency radio data with the Australia Telescope Compact Array are presented to track in more detail the evolution of GX 339–4 and its transient ejecta. During the five years, GX 339–4 has been through two ‘hard-only’ outbursts and two ‘full’ outbursts, allowing us to densely sample the rise, quenching and re-activation of the compact jets. Strong radio flares were also observed close to the transition between the hard and the soft states. Following the radio flare, a transient optically thin ejection was spatially resolved during the 2020 outburst, and was observed for a month. We also discuss the radio/X-ray correlation of GX 339–4 during this five year period, which covers several states in detail from the rising phase to the quiescent state. This campaign allowed us to follow ejection events and provide information on the jet proper motion and its intrinsic velocity. With this work we publicly release the weekly MeerKAT L-band radio maps from data taken between September 2018 and October 2023.