Dr James Matthews

Abstract:

The young neutron star X-ray binary, Cir X-1, resides within its natal supernova remnant and experiences ongoing outbursts every 16.5 d, likely due to periastron passage in an eccentric orbit. We present the deepest ever radio image of the field, which reveals relativistic jet-punched bubbles that are aligned with the mean axis of the smaller scale jets observed close to the X-ray binary core. We are able to measure the minimum energy for the bubble, which is around = erg. The nature and morphological structure of the source were investigated through spectral index mapping and numerical simulations. The spectral index map reveals a large fraction of the nebula’s radio continuum has a steep slope, associated with optically thin synchrotron emission, although there are distinct regions with flatter spectra. Our data are not sensitive enough to measure the spectral index of the protruding bubbles. We used the pluto code to run relativistic hydrodynamic simulations to try and qualitatively reproduce the observations with a combined supernova-plus-jet system. We are able to do so using a simplified model in which the asymmetrical bubbles are best represented by supernova explosion which is closely followed (within 100 yr) by a phase of very powerful jets lasting less than 1000 yr. These are the first observations revealing the initial breakout of neutron star jets from their natal supernova remnant, and further support the scenario in which Cir X-1 is a younger relation of the archetypal jet source SS433.

Abstract:

This work presents comprehensive multifrequency radio monitoring of the black hole low-mass X-ray binary (LMXB) Swift J1727.8−1613, which underwent its first recorded outburst after its discovery in 2023 August. Through a considerable community effort, we have coalesced the data from multiple, distinct observing programs; the light curves include ∼10 months and 197 epochs of monitoring from seven radio facilities with observing frequencies ranging from (approximately) 0.3–230 GHz. The primary purpose of this work is to provide the broader astronomical community with these light curves to assist with the interpretation of other observing campaigns, particularly nonradio observing frequencies. We discuss the phenomenological evolution of the source, which included (i) multiple radio flares consistent with the launching of discrete jet ejections, the brightest of which reached ∼1 Jy; (ii) temporally evolving radio spectral indices (α), reaching values steeper than expected for optically thin synchrotron emission (α < −1) and emission with significant radiative cooling (α < −1.5). We have published a digital copy of the data and intend for this work to set a precedent for the community to continue releasing comprehensive radio light curves of future LMXB outbursts.

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.