Dr James Matthews

Abstract:

Disc winds play a crucial role in many accreting astrophysical systems across all scales. In accreting white dwarfs (AWDs) and active galactic nuclei (AGNs), radiation pressure on spectral lines is a promising wind-driving mechanism. However, the efficiency of line driving is extremely sensitive to the ionization state of the flow, making it difficult to construct a reliable physical picture of these winds. Recently, we presented the first radiation-hydrodynamics simulations for AWDs that incorporated detailed, multidimensional ionization calculations via fully frequency-dependent radiative transfer, using the sirocco code coupled to pluto. These simulations produced much weaker line-driven winds ( for our adopted parameters) than earlier studies using more approximate treatments of ionization and radiative transfer (which yielded ). One remaining limitation of our work was the assumption of an isothermal outflow. Here, we relax this by adopting an ideal gas equation of state and explicitly solving for the multidimensional temperature structure of the flow. In the AWD setting, accounting for the thermal state of the wind does not change the overall conclusions drawn from the isothermal approximation. Our new simulations confirm the line-driving efficiency problem: the predicted outflows are too highly ionized, meaning they neither create optimal driving conditions nor reproduce the observed ultraviolet wind signatures. Possible solutions include wind clumping on subgrid scales, a softer-than-expected spectral energy distribution or additional driving mechanisms. With the physics now built into our simulations, we are well equipped to also explore line-driven disc winds in AGN.

Abstract:

Abstract Black hole X-ray binaries in outburst launch discrete, large-scale jet ejections which can propagate to parsec scales. The kinematics of these ejecta appear to be well described by relativistic blast wave models original devised for gamma-ray burst afterglows. In previous kinematic-only modelling, a crucial degeneracy prevented the initial ejecta energy and the interstellar medium density from being accurately determined. In this work, we present the first joint Bayesian modelling of the radiation and kinematics of a large-scale jet ejection from the X-ray binary MAXI J1535-571. We demonstrate that a reverse shock powers the bright, early ejecta emission. The joint model breaks the energetic degeneracy, and we find the ejecta has an initial energy of E0 ∼ 3 × 1043 erg, and propagates into a low density interstellar medium of nism ∼ 4 × 10−5 cm−3. The ejecta is consistent with being launched perpendicular to the disc and could be powered by an efficient conversion of available accretion power alone. This work lays the foundation for future parameter estimation studies using all available data of X-ray binary jet ejecta.

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.