Dr James Matthews

Abstract:

Tidal disruption events (TDEs) can generate non-spherical, relativistic, and optically thick outflows. Simulations show that the radiation we observe is reprocessed by these outflows. According to a unified model suggested by these simulations, the spectral energy distributions (SEDs) of TDEs depend strongly on viewing angle: low [high] optical-to-X-ray ratios (OXRs) correspond to face-on [edge-on] orientations. Post-processing with radiative transfer codes has simulated the emergent spectra but has so far been carried out only in a quasi-1D framework, with three atomic species (H, He, and O). Here, we present 2.5D Monte Carlo radiative transfer simulations which model the emission from a non-spherical outflow, including a more comprehensive set of cosmically abundant species. While the basic trend of OXR increasing with inclination is preserved, the inherently multi-D nature of photon transport through the non-spherical outflow significantly affects the emergent SEDs. Relaxing the quasi-1D approximation allows photons to preferentially escape in (polar) directions of lower optical depth, resulting in a greater variation of bolometric luminosity as a function of inclination. According to our simulations, inclination alone may not fully explain the large dynamic range of observed TDE OXRs. We also find that including metals, other than O, changes the emergent spectra significantly, resulting in stronger absorption and emission lines in the extreme ultraviolet, as well as a greater variation in the OXR as a function of inclination. Whilst our results support previously proposed unified models for TDEs, they also highlight the critical importance of multi-D ionization and radiative transfer.

Abstract:

We present numerical simulations of discrete relativistic ejecta from an X-ray binary (XRB) with initial conditions directly informed by observations. XRBs have been observed to launch powerful discrete plasma ejecta during state transitions, which can propagate up to parsec distances. Understanding these ejection events unveils new understanding of jet-launching, jet power, and jet–interstellar medium (ISM) interaction among other implications. Multifrequency quasi-simultaneous radio observations of ejecta from the black hole XRB MAXI J1820+070 produced both size and calorimetry constraints, which we use as initial conditions of a relativistic hydrodynamic simulation. We qualitatively reproduce the observed deceleration of the ejecta in a homogeneous ISM. Our simulations demonstrate that the ejecta must be denser than the ISM, the ISM be significantly low density, and the launch be extremely powerful, in order to propagate to the observed distances. The blob propagates and clears out a high-pressure low-density cavity in its wake, providing an explanation for this pre-existing low-density environment, as well as ‘bubble-like’ environments in the vicinity of XRBs inferred from other studies. As the blob decelerates, we observe the onset of instabilities and a long-lived reverse shock – these mechanisms convert kinetic to internal energy in the blob, responsible for in situ particle acceleration. We transform the outputs of our simulation into pseudo-radio images, incorporating the coverage of the MeerKAT and e-MERLIN telescopes from the original observations with real-sky background. Through this, we maximize the interpretability of the results and provide direct comparison to current data, as well as provide prediction capabilities.

Abstract:

As the closest known active galactic nucleus, Centaurus A (Cen A) provides a rich environment for astrophysical exploration. It has been observed across wavelengths from radio to gamma-rays, and indications of ongoing particle acceleration have been found on different scales. Recent measurements of very-high-energy (VHE) gamma-rays ( GeV) by the HESS observatory have inferred the presence of ultra-relativistic electrons along Cen A’s jet, yet the underlying acceleration mechanism remains uncertain. Various authors have proposed that jet substructures, known as knots, may serve as efficient particle accelerators. In this study, we investigate the hypothesis that knots are the particle acceleration sites along Cen A’s jets. We focus on stationary knots, and assume that they result from interactions between the jet and the stellar winds of powerful stars. By combining relativistic hydrodynamic simulations and shock acceleration theory with the radio and X-ray data, we compare theoretical predictions with morphological and spectral data from different knots. We estimate the maximum electron energy and the resulting VHE gamma-ray emission. Our findings suggest that electrons accelerated at the knots are responsible for the gamma-ray spectrum detected in the VHE band.

Abstract:

Blast wave models are commonly used to model relativistic outflows from ultra-relativistic gamma-ray bursts (GRBs), but are also applied to lower Lorentz factor ejections from X-ray binaries (XRBs). Here, we revisit the physics of blast waves and reverse shocks in these systems and explore the similarities and differences between the ultra-relativistic () and moderately relativistic () regimes. We first demonstrate that the evolution of the blast wave radius as a function of the observer frame time is recovered in the on-axis ultra-relativistic limit from a general energy and radius blast wave evolution, emphasizing that XRB ejections are off-axis, moderately relativistic cousins of GRB afterglows. We show that, for fixed blast wave or ejecta energy, reverse shocks cross the ejecta much later (earlier) on in the evolution for less (more) relativistic systems, and find that reverse shocks are much longer lived in XRBs and off-axis GRBs compared to on-axis GRBs. Reverse shock crossing should thus typically finish after 10–100 of days (in the observer frame) in XRB ejections. This characteristic, together with their moderate Lorentz factors and resolvable core separations, makes XRB ejections unique laboratories for shock and particle acceleration physics. We discuss the impact of geometry and lateral spreading on our results, explore how to distinguish between different shock components, and comment on the implications for GRB and XRB environments. Additionally, we argue that identification of reverse shock signatures in XRBs could provide an independent constraint on the ejecta Lorentz factor.