Pulsars, transients and relativistic astrophysics

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:

Supermassive black hole (SMBH) binary systems are an unavoidable outcome of galaxy mergers. Their dynamics encode valuable information about their formation and growth, the composition of their host galactic nuclei, the evolution of galaxies, and the nature of gravity. Many SMBH binaries with separations pc-kpc have been found, but closer (subparsec) binaries remain to be confirmed. Identifying these systems may elucidate how binaries evolve past the “final parsec” until gravitational radiation drives them to coalescence. Methods to discover and characterize SMBH binaries can shed light on these important questions and potentially open new multimessenger channels. Here we show that SMBH binaries in nonactive galactic nuclei can be identified and characterized by the gravitational lensing of individual bright stars, located behind them in the host galaxy. The rotation of “caustics”—regions where sources are hugely magnified due to the SMBH binary’s orbit and inspiral—leads to quasiperiodic lensing of starlight (QPLS). The extreme lensing magnification of individual bright stars produces a significant variation in the host galaxies’ luminosity; their lightcurve traces the orbit of the SMBH binary and its evolution, analogous to the waveforms recorded by gravitational-wave (GW) detectors. QPLS probes the population of sources observable by pulsar timing arrays and space detectors (LISA, TianQin), offering advance warning triggers for merging SMBHs for coincident or follow-up GW detections. SMBH population models predict 1–50 $[190–5000](n_{\star }/{\mathrm{pc}}^{-3})$ QPLS binaries with period less than 10[40] yr with comparable masses and redshift $z<0.3$ , where $n_{\star }$ is the stellar number density. Additionally, stellar and orbital motion will lead to frequent instances of single or double flares caused by SMBHBs with longer periods. This novel signature can be searched for in a wealth of existing and upcoming time-domain photometric data: identifying quasiperiodic variability in galactic lightcurves will reveal an ensemble of binary systems and illuminate outstanding questions around them.

Abstract:

We present the discovery of radio emission from the Be/X-ray binary A0538–66 with the Australian Square Kilometre Array Pathfinder (ASKAP), and results from a subsequent weekly monitoring campaign with the MeerKAT radio telescope. A0538–66, located in the Large Magellanic Cloud, hosts a neutron star with a short spin period (P ≈ 69 ms) in a highly eccentric ≈16.6-day orbit. Its rare episodes of super-Eddington accretion, rapid optical and X-ray flares, and other peculiar properties make it an interesting system among high-mass X-ray binaries. Our MeerKAT data reveal that it is also one of the most radio-luminous neutron star X-ray binaries observed to date, reaching ≈3 × 1022 erg s−1Hz−1, with radio emission that appears to be orbitally modulated. We consider several possible mechanisms for the radio emission, and place A0538–66 in context by comparing it to similar systems.

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.