Pulsars, transients and relativistic astrophysics

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.