Showing 1–47 of 47 results for author: Matthews, J H

A disc wind origin for the optical spectra of dwarf novae in outburst

Authors: Yusuke Tampo, Christian Knigge, Knox S. Long, James H. Matthews, Noel Castro Segura

Abstract: Many high-state cataclysmic variables (CVs) exhibit blue-shifted absorption features in their ultraviolet (UV) spectra -- a smoking-gun signature of outflows. However, the impact of these outflows on {\em optical} spectra remains much more uncertain. During its recent outburst, the eclipsing dwarf nova V455 And displayed strong optical emission lines whose cores were narrower than expected from a… ▽ More Many high-state cataclysmic variables (CVs) exhibit blue-shifted absorption features in their ultraviolet (UV) spectra -- a smoking-gun signature of outflows. However, the impact of these outflows on {\em optical} spectra remains much more uncertain. During its recent outburst, the eclipsing dwarf nova V455 And displayed strong optical emission lines whose cores were narrower than expected from a Keplerian disc. Here, we explore whether disc + wind models developed for matching UV observations of CVs can also account for these optical spectra. Importantly, V455~And was extremely bright at outburst maximum: the accretion rate implied by fitting the optical continuum with a standard disc model is $\dot{M}_{\rm acc} \simeq 10^{-7}~{\rm M}_\odot~{\rm yr^{-1}}$. Allowing for continuum reprocessing in the outflow helps to relax this constraint. A disk wind can also broadly reproduce the optical emission lines, but only if the wind is (i) highly mass-loaded, with a mass-loss rate reaching $\dot{M}_{\rm wind} \simeq 0.4 \dot{M}_{\rm acc}$, and/or (ii) clumpy, with a volume filling factor $f_V \simeq 0.1$. The same models can describe the spectral evolution across the outburst, simply by lowering $\dot{M}_{\rm acc}$ and $\dot{M}_{\rm wind}$. Extending these models to lower inclinations and into the UV produces spectra consistent with those observed in face-on high-state CVs. We also find, for the first time in simulations of this type, P-Cygni-like absorption features in the Balmer series, as have been observed in both CVs and X-ray binaries. Overall, dense disc winds provide a promising framework for explaining multiple observational signatures seen in high-state CVs, but theoretical challenges persist. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 13 pages, 9 figures, 3 tables. Accepted for publication in MNRAS

[OIII] emission in z=2 quasars with and without Broad Absorption Lines

Authors: Matthew J. Temple, Amy L. Rankine, Manda Banerji, Joseph F. Hennawi, Paul C. Hewett, James H. Matthews, Riccardo Nanni, Claudio Ricci, Gordon T. Richards

Abstract: Understanding the links between different phases of outflows from active galactic nuclei is a key goal in extragalactic astrophysics. Here we compare [OIII] $λλ$4960,5008 outflow signatures in quasars with and without Broad Absorption Lines (BALs), aiming to test how the broad absorption troughs seen in the rest-frame ultraviolet are linked to the narrow line region outflows seen in the rest-frame… ▽ More Understanding the links between different phases of outflows from active galactic nuclei is a key goal in extragalactic astrophysics. Here we compare [OIII] $λλ$4960,5008 outflow signatures in quasars with and without Broad Absorption Lines (BALs), aiming to test how the broad absorption troughs seen in the rest-frame ultraviolet are linked to the narrow line region outflows seen in the rest-frame optical. We present new near-infrared spectra from Magellan/FIRE which cover [OIII] in 12 quasars with 2.1 < z < 2.3, selected to have strong outflow signatures in CIV $λ$1550. Combining with data from the literature, we build a sample of 73 BAL, 115 miniBAL and 125 non-BAL QSOs with 1.5 < z < 2.6. The strength and velocity width of [OIII] correlate strongly with the CIV emission properties, but no significant difference is seen in the [OIII] emission-line properties between the BALs, non-BALs and miniBALs once the dependence on CIV emission is taken into account. A weak correlation is observed between the velocities of CIV BALs and [OIII] emission, which is accounted for by the fact that both outflow signatures correlate with the underlying CIV emission properties. Our results add to the growing evidence that BALs and non-BALs are drawn from the same parent population and are consistent with a scenario wherein BAL troughs are intermittent tracers of persistent quasar outflows, with a part of such outflow becoming optically thick along our line-of-sight for sporadic periods of time within which BALs are observed. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 14 pages, 6 figures. Accepted for publication in MNRAS

Journal ref: MNRAS 532, 424-437 (2024)

Getting More Out of Black Hole Superradiance: a Statistically Rigorous Approach to Ultralight Boson Constraints

Authors: Sebastian Hoof, David J. E. Marsh, Júlia Sisk-Reynés, James H. Matthews, Christopher Reynolds

Abstract: Black hole (BH) superradiance can provide strong constraints on the properties of ultralight bosons (ULBs). Since most of the previous work has focused on the theoretical predictions, here we investigate the most suitable statistical framework to constrain ULB masses and self-interactions. We argue that a Bayesian approach provides a clear statistical interpretation, deals with limitations regardi… ▽ More Black hole (BH) superradiance can provide strong constraints on the properties of ultralight bosons (ULBs). Since most of the previous work has focused on the theoretical predictions, here we investigate the most suitable statistical framework to constrain ULB masses and self-interactions. We argue that a Bayesian approach provides a clear statistical interpretation, deals with limitations regarding the reproducibility of existing BH analyses, incorporates the full information from BH data, and allows us to include additional nuisance parameters or to perform hierarchical modelling with BH populations in the future. We demonstrate the feasibility of our approach using mass and spin posterior samples for the X-ray binary BH M33 X-7 and, for the first time in this context, the supermassive BH IRAS 09149-6206. We explain the differences to existing ULB constraints in the literature and illustrate the effects of various assumptions about the superradiance process (equilibrium regime vs cloud collapse, higher occupation levels). As a result, our procedure yields the most rigorous ULB constraints available in the literature, with important implications for the QCD axion and axion-like particles. We encourage all groups analysing BH data to publish likelihood functions or posterior samples as supplementary material to facilitate this type of analysis. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Comments: 13+2 pages, 6 figures, software code available at https://github.com/sebhoof/bhsr

State-of-the-art simulations of line-driven accretion disc winds: realistic radiation-hydrodynamics leads to weaker outflows

Authors: Nick Higginbottom, Nicolas Scepi, Christian Knigge, Knox S. Long, James H. Matthews, Stuart A. Sim

Abstract: Disc winds are a common feature in accreting astrophysical systems on all scales. In active galactic nuclei (AGN) and accreting white dwarfs (AWDs), specifically, radiation pressure mediated by spectral lines is a promising mechanism for driving these outflows. Previous hydrodynamical simulations have largely supported this idea, but relied on highly approximate treatments of ionization and radiat… ▽ More Disc winds are a common feature in accreting astrophysical systems on all scales. In active galactic nuclei (AGN) and accreting white dwarfs (AWDs), specifically, radiation pressure mediated by spectral lines is a promising mechanism for driving these outflows. Previous hydrodynamical simulations have largely supported this idea, but relied on highly approximate treatments of ionization and radiative transfer. Given the sensitivity of line driving to the ionization state and radiation field in the outflow, here we present a new method for carrying out 2.5D radiation-hydrodynamic simulations that takes full account of the frequency-dependent radiative transfer through the wind, the corresponding ionization state and the resulting radiative accelerations. Applying our method to AWDs, we find that it is much harder to drive a powerful line-driven outflow when the interaction between matter and radiation is treated self-consistently. This conclusion is robust to changes in the adopted system parameters. The fundamental difficulty is that discs luminous enough to drive such a wind are also hot enough to over-ionize it. As a result, the mass-loss rates in our simulations are much lower than those found in earlier, more approximate calculations. We also show that the ultraviolet spectra produced by our simulations do not match those observed in AWDs. We conclude that, unless the over-ionization problem can be mitigated (e.g. by sub-grid clumping or a softer-than-expected radiation field), line driving may not be a promising mechanism for powering the outflows from AWDs. These conclusions are likely to have significant implications for disc winds in AGN also. △ Less

Submitted 10 December, 2023; originally announced December 2023.

Comments: Accepted for publication in MNRAS. 14 pages, 10 figures + 3 figures in Appendix

Shedding far-ultraviolet light on the donor star and evolutionary state of the neutron-star LMXB Swift J1858.6-0814

Authors: N. Castro Segura, C. Knigge, J. H. Matthews, F. M. Vincentelli, P. Charles, K. S. Long, D. Altamirano, D. A. H. Buckley, D. Modiano, M. A. P. Torres, D. J. K. Buisson, S. Fijma, K. Alabarta, N. Degenaar, M. Georganti, M. C. Baglio

Abstract: The evolution of accreting X-ray binary systems is closely coupled to the properties of their donor stars. As a result, we can constrain the evolutionary track a system is by establishing the nature of its donor. Here, we present far-UV spectroscopy of the transient neutron-star low-mass X-ray binary Swift J1858 in three different accretion states (low-hard, high-hard and soft). All of these spect… ▽ More The evolution of accreting X-ray binary systems is closely coupled to the properties of their donor stars. As a result, we can constrain the evolutionary track a system is by establishing the nature of its donor. Here, we present far-UV spectroscopy of the transient neutron-star low-mass X-ray binary Swift J1858 in three different accretion states (low-hard, high-hard and soft). All of these spectra exhibit anomalous N\,{\sc v}, C\,{\sc iv}, Si\,{\sc iv} and He\,{\sc ii} lines, suggesting that its donor star has undergone CNO processing. We also determine the donor's effective temperature, $T_{d} \simeq 5700$~K, and radius, $R_d \simeq 1.7~R_{\odot}$, based on photometric observations obtained during quiescence. Lastly, we leverage the transient nature of the system to set an upper limit of $\dot{M}_{\rm acc} \lesssim 10^{-8.5}~M_{\odot}~yr^{-1}$ on the present-day mass-transfer rate. Combining all these with the orbital period of the system, $P_{\rm orb} = 21.3$~hrs, we search for viable evolution paths. The initial donor masses in the allowed solutions span the range $1~M_{\odot} \lesssim M_{d,i} \lesssim 3.5~M_{\odot}$. All but the lowest masses in this range are consistent with the strong CNO-processing signature in the UV line ratios. The present-day donor mass in the permitted tracks are $0.5~M_{\odot}\lesssim M_{d,obs} \lesssim 1.3~M_{\odot}$, higher than suggested by recent eclipse modelling. Since $P_{\rm orb}$ is close to the so-called bifurcation period, both converging and diverging binary tracks are permitted. If Swift J1858 is on a converging track, it will end its life as an ultra-compact system with a sub-stellar donor star. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: Accepted for publication in MNRAS

A disc wind model for blueshifts in quasar broad emission lines

Authors: James H. Matthews, Jago Strong-Wright, Christian Knigge, Paul Hewett, Matthew J. Temple, Knox S. Long, Amy L. Rankine, Matthew Stepney, Manda Banerji, Gordon T. Richards

Abstract: Blueshifts - or, more accurately, blue asymmetries - in broad emission lines such as CIV $λ$1550 are common in luminous quasars and correlate with fundamental properties such as Eddington ratio and broad absorption line (BAL) characteristics. However, the formation of these blueshifts is still not understood, and neither is their physical connection to the BAL phenomenon or accretion disc. In this… ▽ More Blueshifts - or, more accurately, blue asymmetries - in broad emission lines such as CIV $λ$1550 are common in luminous quasars and correlate with fundamental properties such as Eddington ratio and broad absorption line (BAL) characteristics. However, the formation of these blueshifts is still not understood, and neither is their physical connection to the BAL phenomenon or accretion disc. In this work, we present Monte Carlo radiative transfer and photoionization simulations using parametrized biconical disc-wind models. We take advantage of the azimuthal symmetry of a quasar and show that we can reproduce CIV blueshifts provided that (i) the disc-midplane is optically thick out to radii beyond the line formation region, so that the receding wind bicone is obscured; and (ii) the system is viewed from relatively low (that is, more face-on) inclinations ($\lesssim40^\circ$). We show that CIV emission line blueshifts and BALs can form in the same wind structure. The velocity profile of the wind has a significant impact on the location of the line formation region and the resulting line profile, suggesting that the shape of the emission lines can be used as a probe of wind-driving physics. While we are successful at producing blueshifts/blue asymmetries in outflows, we struggle to match the detailed shape or skew of the observed emission line profiles. In addition, our models produce redshifted emission-line asymmetries for certain viewing angles. We discuss our work in the context of the CIV $λ$1550 emission blueshift versus equivalent-width space and explore the implications for quasar disc wind physics. △ Less

Submitted 25 September, 2023; originally announced September 2023.

Comments: 20 pages, 12 figures + 1 appendix. Supplementary material, simulation data and plotting scripts available at: https://github.com/jhmatthews/blueshifts

Journal ref: MNRAS 526, pp.3967-3986 (2023)

No redshift evolution in the rest-frame UV emission line properties of quasars from z=1.5 to z=4.0

Authors: Matthew Stepney, Manda Banerji, Paul C. Hewett, Matthew J. Temple, Amy L. Rankine, James H. Matthews, Gordon T. Richards

Abstract: We analyse the rest-frame UV spectra of 2,531 high-redshift (3.5<z<4.0) quasars from the Sloan Digital Sky Survey DR16Q catalogue. In combination with previous work, we study the redshift evolution of the rest-frame UV line properties across the entire redshift range, 1.5<z<4.0. We improve the systemic redshift estimates at z>3.5 using a cross-correlation algorithm that employs high signal-to-nois… ▽ More We analyse the rest-frame UV spectra of 2,531 high-redshift (3.5<z<4.0) quasars from the Sloan Digital Sky Survey DR16Q catalogue. In combination with previous work, we study the redshift evolution of the rest-frame UV line properties across the entire redshift range, 1.5<z<4.0. We improve the systemic redshift estimates at z>3.5 using a cross-correlation algorithm that employs high signal-to-noise template spectra spanning the full range in UV emission line properties. We then quantify the evolution of C IV and He II emission line properties with redshift. The increase in C IV blueshifts with cosmological redshift can be fully explained by the higher luminosities of quasars observed at high redshifts. We recover broadly similar trends between the He II EW and C IV blueshift at both 1.5<z<2.65 and 3.5<z<4.0 suggesting that the blueshift depends systematically on the spectral energy density (SED) of the quasar and there is no evolution in the SED over the redshift range 1.5<z<4.0. C IV blueshifts are highest when L/LEdd > 0.2 and Mbh > 10^9 Mo for the entire 1.5<z<4.0 sample. We find that luminosity matching samples as a means to explore the evolution of their rest-frame UV emission line properties is only viable if the samples are also matched in the Mbh - L/LEdd plane. Quasars at z>6 are on average less massive and have higher Eddington-scaled accretion rates than their luminosity-matched counterparts at 1.5<z<4.0, which could explain the observed evolution in their UV line properties. △ Less

Submitted 6 July, 2023; originally announced July 2023.

Comments: Accepted for publication in MNRAS, July 5th 2023. Supplementary material will be available online through the journal

Studying the link between radio galaxies and AGN fuelling with relativistic hydrodynamic simulations of flickering jets

Authors: Henry W. Whitehead, James H. Matthews

Abstract: We present two- and three-dimensional hydrodynamic simulations of $\sim$kpc-scale AGN jets with mean jet powers in the range $1-7\times10^{45}\,$erg~s$^{-1}$, in which the jet power varies (through variation of the Lorentz factor) according to a flicker or pink noise power spectrum. We find the morphology and dynamics of the jet-cocoon system depends on the amplitude of the variability with a clea… ▽ More We present two- and three-dimensional hydrodynamic simulations of $\sim$kpc-scale AGN jets with mean jet powers in the range $1-7\times10^{45}\,$erg~s$^{-1}$, in which the jet power varies (through variation of the Lorentz factor) according to a flicker or pink noise power spectrum. We find the morphology and dynamics of the jet-cocoon system depends on the amplitude of the variability with a clear correspondence between the shape of the cocoon and the historical activity. The jet advances quickly during high-power states, whereas quiescent periods instead produce passive periods of inflation resembling Sedov-Taylor blast waves. Periods of high activity preferentially produce hotspots and create stronger backflow as they maximise the pressure gradient between the jet head and cocoon. The variability can also lead to propagating internal shock structures along the jet. Our work suggests that variability and flickering in the jet power has important implications, which we discuss, for observations of radio galaxies, ultrahigh energy cosmic ray acceleration and jet power to luminosity correlations. We explore the link between morphology and fuelling, and suggest that chaotic cold accretion should introduce a relatively small scatter in radio luminosity ($\sim0.2$ dex) and modest imprints on morphology; sources such as Hercules A and Fornax A, which show evidence for more dramatic variability, may therefore require redder power spectra, or be triggered by mergers or other discrete events. We suggest ways to search for jet flickering observationally and propose that radio galaxies may be an important diagnostic of Myr timescale AGN fuelling, due to their `long-term memory'. △ Less

Submitted 30 May, 2023; originally announced May 2023.

Comments: 20 pages, 15 figures, 1 table. Accepted for publication in MNRAS. Animated version of figure 3 available at https://youtu.be/YhDHokWXuKk with other movies on the same youtube channel. Full supplementary material available at https://github.com/jhmatthews/flicker-supplement

MAXI J1848-015: The First Detection of Relativistically Moving Outflows from a Globular Cluster X-ray Binary

Authors: A. Bahramian, E. Tremou, A. J. Tetarenko, J. C. A. Miller-Jones, R. P. Fender, S. Corbel, D. R. A. Williams, J. Strader, F. Carotenuto, R. Salinas, J. A. Kennea, S. E. Motta, P. A. Woudt, J. H. Matthews, T. D. Russell

Abstract: Over the past decade, observations of relativistic outflows from outbursting X-ray binaries in the Galactic field have grown significantly. In this work, we present the first detection of moving and decelerating radio-emitting outflows from an X-ray binary in a globular cluster. MAXI J1848-015 is a recently discovered transient X-ray binary in the direction of the globular cluster GLIMPSE-C01. Usi… ▽ More Over the past decade, observations of relativistic outflows from outbursting X-ray binaries in the Galactic field have grown significantly. In this work, we present the first detection of moving and decelerating radio-emitting outflows from an X-ray binary in a globular cluster. MAXI J1848-015 is a recently discovered transient X-ray binary in the direction of the globular cluster GLIMPSE-C01. Using observations from the VLA, and a monitoring campaign with the MeerKAT observatory for 500 days, we model the motion of the outflows. This represents some of the most intensive, long-term coverage of relativistically moving X-ray binary outflows to date. We use the proper motions of the outflows from MAXI J1848-015 to constrain the component of the intrinsic jet speed along the line of sight, $β_\textrm{int} \cos θ_\textrm{ejection}$, to be $=0.19\pm0.02$. Assuming it is located in GLIMPSE-C01, at 3.4 kpc, we determine the intrinsic jet speed, $β_\textrm{int}=0.79\pm0.07$, and the inclination angle to the line of sight, $θ_\textrm{ejection}=76^\circ\pm2^{\circ}$. This makes the outflows from MAXI J1848-015 somewhat slower than those seen from many other known X-ray binaries. We also constrain the maximum distance to MAXI J1848-015 to be $4.3$ kpc. Lastly, we discuss the implications of our findings for the nature of the compact object in this system, finding that a black hole primary is a viable (but as-of-yet unconfirmed) explanation for the observed properties of MAXI J1848-015. If future data and/or analysis provide more conclusive evidence that MAXI J1848-015 indeed hosts a black hole, it would be the first black hole X-ray binary in outburst identified in a Galactic globular cluster. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Comments: 14 pages, 5 figures, ApJL in press

Current and Future constraints on Very-Light Axion-Like Particles from X-ray observations of cluster-hosted Active Galaxies

Authors: Julia M. Sisk-Reynes, Christopher S. Reynolds, James H. Matthews

Abstract: We discuss our recent constraints on the coupling of Very-Light Axion-Like Particles (of masses $<$$ 10^{-12} \ \mathrm{eV}$) to electromagnetism from $Chandra$ observations of the cluster-hosted Active Galactic Nuclei (AGN) H1821+643 and NGC1275. In both cases, the inferred high-quality AGN spectra excluded all photon-ALP couplings… ▽ More We discuss our recent constraints on the coupling of Very-Light Axion-Like Particles (of masses $<$$ 10^{-12} \ \mathrm{eV}$) to electromagnetism from $Chandra$ observations of the cluster-hosted Active Galactic Nuclei (AGN) H1821+643 and NGC1275. In both cases, the inferred high-quality AGN spectra excluded all photon-ALP couplings $g_\mathrm{aγ} > (6.3 - 8.0) \times 10^{-13} \ {\mathrm{GeV}}^{-1}$ at the $99.7\%$ level, respectively, based on the non-detection of spectral distortions attributed to photon-ALP inter-conversion along the cluster line-of-sight. Finally, we present the prospects of tightening current bounds on such ALPs by up to a factor of 10 with next-generation X-ray observatories such as $Athena$, $AXIS$ and $LEM$ given their improved spectral and spatial resolution and collecting area compared to current missions. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 4 pages, 1 figure. Accepted for publication to Memorie della SAIt for the Proceedings of the European Astronomical Society 2022 (EAS 2022) Annual Meeting Symposium S3 "The Dark matter multi-messenger challenge"

The origin of optical emission lines in the soft state of X-ray binary outbursts: the case of MAXI J1820+070

Authors: K. I. I. Koljonen, K. S. Long, J. H. Matthews, C. Knigge

Abstract: The optical emission line spectra of X-ray binaries (XRBs) are thought to be produced in an irradiated atmosphere, possibly the base of a wind, located above the outer accretion disc. However, the physical nature of - and physical conditions in - the line-forming region remain poorly understood. Here, we test the idea that the optical spectrum is formed in the transition region between the cool, g… ▽ More The optical emission line spectra of X-ray binaries (XRBs) are thought to be produced in an irradiated atmosphere, possibly the base of a wind, located above the outer accretion disc. However, the physical nature of - and physical conditions in - the line-forming region remain poorly understood. Here, we test the idea that the optical spectrum is formed in the transition region between the cool, geometrically thin part of the disc near the mid-plane and a hot, vertically extended atmosphere or outflow produced by X-ray irradiation. We first present a VLT X-Shooter spectrum of XRB MAXI J1820+070 in the soft state associated with its 2018 outburst, which displays a rich set of double-peaked hydrogen and helium recombination lines. Aided by ancillary X-ray spectra and reddening estimates, we then model this spectrum with the Monte Carlo radiative transfer code Python, using a simple biconical disc wind model inspired by radiation-hydrodynamic simulations of irradiation-driven outflows from XRB discs. Such a model can qualitatively reproduce the observed features; nearly all of the optical emission arising from the transonic 'transition region' near the base of the wind. In this region, characteristic electron densities are on the order of 10$^{12-13}$ cm$^{-3}$, in line with the observed flat Balmer decrement (H$α$/H$β\approx 1.3$). We conclude that strong irradiation can naturally give rise to both the optical line-forming layer in XRB discs and an overlying outflow/atmosphere that produces X-ray absorption lines. △ Less

Submitted 16 March, 2023; originally announced March 2023.

Comments: 17 pages, 11 figures. This article has been accepted for publication in MNRAS published by Oxford University Press on behalf of the Royal Astronomical Society

UHECR Echoes from the Council of Giants

Authors: A. M. Taylor, J. H. Matthews, A. R. Bell

Abstract: Recent anisotropy studies of UHECR data at energies $\gtrsim$ 40 EeV, have disclosed a correlation of their angular distribution with the extragalactic local structure, specifically with either local starburst galaxies or AGN. Using Monte Carlo simulations taking into account photo-disintegration processes, we further explore a framework in which these UHECRs were accelerated by Centaurus A in a r… ▽ More Recent anisotropy studies of UHECR data at energies $\gtrsim$ 40 EeV, have disclosed a correlation of their angular distribution with the extragalactic local structure, specifically with either local starburst galaxies or AGN. Using Monte Carlo simulations taking into account photo-disintegration processes, we further explore a framework in which these UHECRs were accelerated by Centaurus A in a recent powerful outburst before being scattered by magnetic fields associated with local, Council of Giant, extragalactic structure. We find that the observed intermediate scale anisotropies can be accounted for by the Council of Giant structure imposing a response function on the initial outburst of UHECRs from a single source located at Centaurus A's position. The presence of these local structures creates `echoes' of UHECRs after the initial impulse, and focusing effects. The strongest echo wave has a lag of $\sim$ 20 Myr, comparable to the age of synchrotron-emitting electrons in the giant Centaurus A lobes. Through consideration of the composition of both the direct and echo wave components, we find that the distribution of the light (1$<\ln A<$1.5) component across the sky offers exciting prospects for testing the echo model using future facilities such as Auger prime. Our results demonstrate the potential that UHECR nuclei offer, as "composition clocks", for probing propagation scenarios from local sources. △ Less

Submitted 13 June, 2023; v1 submitted 13 February, 2023; originally announced February 2023.

Comments: 12 pages, 11 figures

How, where and when do cosmic rays reach ultrahigh energies?

Authors: James H. Matthews, Andrew M. Taylor

Abstract: Understanding the origins of ultrahigh energy cosmic rays (UHECRs) - which reach energies in excess of $10^{20}~{\rm eV}$ - stretches particle acceleration physics to its very limits. In this review, we discuss how such energies can be reached, using general arguments that can often be derived on the back of an envelope. We explore possible particle acceleration mechanisms, with special attention… ▽ More Understanding the origins of ultrahigh energy cosmic rays (UHECRs) - which reach energies in excess of $10^{20}~{\rm eV}$ - stretches particle acceleration physics to its very limits. In this review, we discuss how such energies can be reached, using general arguments that can often be derived on the back of an envelope. We explore possible particle acceleration mechanisms, with special attention paid to shock acceleration. Informed by the arguments derived, we discuss where UHECRs might come from and which classes of powerful astrophysical objects could be UHECR sources; generally, we favour radio galaxies, GRB afterglows and other sources which are not too compact and dissipate prodigious amounts of energy on large scales, allowing them to generate large products $βB R$ without the CRs undergoing restrictive losses. Finally, we discuss when UHECRs are accelerated by highlighting the importance of source variability, and explore the intriguing possibility that the UHECR arrival directions are partly a result of "echoes" from magnetic structures in the local Universe. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: Review contribution submitted to proceedings of science as part of the proceedings of the 27th European Cosmic Ray Symposium held 25-29 July 2022. 25 pages, 5 figures, 1 table. Comments welcome

Testing AGN outflow and accretion models with CIV and HeII emission line demographics in z=2 quasars

Authors: Matthew J. Temple, James H. Matthews, Paul C. Hewett, Amy L. Rankine, Gordon T. Richards, Manda Banerji, Gary J. Ferland, Christian Knigge, Matthew Stepney

Abstract: Using 190,000 spectra from the seventeenth data release of the Sloan Digital Sky Survey, we investigate the ultraviolet emission line properties in z=2 quasars. Specifically, we quantify how the shape of CIV 1549A and the equivalent width (EW) of HeII 1640A depend on the black hole mass and Eddington ratio inferred from MgII 2800A. Above L/L_Edd>0.2, there is a strong mass dependence in both CIV b… ▽ More Using 190,000 spectra from the seventeenth data release of the Sloan Digital Sky Survey, we investigate the ultraviolet emission line properties in z=2 quasars. Specifically, we quantify how the shape of CIV 1549A and the equivalent width (EW) of HeII 1640A depend on the black hole mass and Eddington ratio inferred from MgII 2800A. Above L/L_Edd>0.2, there is a strong mass dependence in both CIV blueshift and HeII EW. Large CIV blueshifts are observed only in regions with both high mass and high accretion rate. Including X-ray measurements for a subsample of 5,000 objects, we interpret our observations in the context of AGN accretion and outflow mechanisms. The observed trends in HeII and 2 keV strength are broadly consistent with theoretical QSOSED models of AGN spectral energy distributions (SEDs) for low spin black holes, where the ionizing SED depends on the accretion disc temperature and the strength of the soft excess. High spin models are not consistent with observations, suggesting SDSS quasars at z=2 may in general have low spins. We find a dramatic switch in behaviour at L/L_Edd<0.2: the ultraviolet emission properties show much weaker trends, and no longer agree with QSOSED predictions, hinting at changes in the structure of the broad line region. Overall the observed emission line trends are generally consistent with predictions for radiation line driving where quasar outflows are governed by the SED, which itself results from the accretion flow and hence depends on both the SMBH mass and accretion rate. △ Less

Submitted 27 April, 2023; v1 submitted 6 January, 2023; originally announced January 2023.

Comments: Resubmitted to MNRAS. Minor edits to text; figures updated; licence updated

Journal ref: MNRAS 523, 646 (2023)

Physics Beyond the Standard Model with Future X-ray Observatories: Projected Constraints on Very-Light Axion-Like Particles with $Athena$ and $AXIS$

Authors: Júlia Sisk-Reynés, Christopher S. Reynolds, Michael L. Parker, James H. Matthews, M. C. David Marsh

Abstract: Axion-Like Particles (ALPs) are well-motivated extensions of the Standard Model of Particle Physics and a generic prediction of some string theories. X-ray observations of bright Active Galactic Nuclei (AGN) hosted by rich clusters of galaxies are excellent probes of very-light ALPs, with masses $\mathrm{log}(m_\mathrm{a}/\mathrm{eV}) < -12.0$. We evaluate the potential of future X-ray observatori… ▽ More Axion-Like Particles (ALPs) are well-motivated extensions of the Standard Model of Particle Physics and a generic prediction of some string theories. X-ray observations of bright Active Galactic Nuclei (AGN) hosted by rich clusters of galaxies are excellent probes of very-light ALPs, with masses $\mathrm{log}(m_\mathrm{a}/\mathrm{eV}) < -12.0$. We evaluate the potential of future X-ray observatories, particularly $Athena$ and the proposed $AXIS$, to constrain ALPs via observations of cluster-hosted AGN, taking NGC 1275 in the Perseus cluster as our exemplar. Assuming perfect knowledge of instrument calibration, we show that a modest exposure (200-ks) of NGC 1275 by $Athena$ permits us to exclude all photon-ALP couplings $g_\mathrm{aγ} > 6.3 \times 10^{-14} \ {\mathrm{GeV}}^{-1}$ at the 95% level, as previously shown by $Conlon \ et \ al. \ (2018)$, representing a factor of 10 improvement over current limits. We then proceed to assess the impact of realistic calibration uncertainties on the $Athena$ projection by applying a standard $Cash$ likelihood procedure, showing the projected constraints on $g_\mathrm{aγ}$ weaken by a factor of 10 (back to the current most sensitive constraints). However, we show how the use of a deep neural network can disentangle the energy-dependent features induced by instrumental miscalibration and those induced by photon-ALP mixing, allowing us to recover most of the sensitivity to the ALP physics. In our explicit demonstration, the machine learning applied allows us to exclude $g_\mathrm{aγ} > 2.0 \times 10^{-13} \ {\mathrm{GeV}}^{-1}$, complementing the projected constraints of next-generation ALP dark matter birefringent cavity searches for very-light ALPs. Finally, we show that a 200-ks $AXIS$/on-axis observation of NGC 1275 will tighten the current best constraints on very-light ALPs by a factor of 3. △ Less

Submitted 25 April, 2023; v1 submitted 9 November, 2022; originally announced November 2022.

Comments: 22 pages, 10 figures, 3 tables. Accepted for publication at the Astrophysical Journal. Comments welcome

Connecting radio emission to AGN wind properties with Broad Absorption Line Quasars

Authors: J. W. Petley, L. K. Morabito, D. M. Alexander, A. L. Rankine, V. A. Fawcett, D. J. Rosario, J. H. Matthews, T. M. Shimwell, A. Drabent

Abstract: Broad Absorption Line Quasars (BALQSOs) show strong signatures of powerful outflows, with the potential to alter the cosmic history of their host galaxies. These signatures are only seen in ~10% of optically selected quasars, although the fraction significantly increases in IR and radio selected samples. A proven physical explanation for this observed fraction has yet to be found, along with a det… ▽ More Broad Absorption Line Quasars (BALQSOs) show strong signatures of powerful outflows, with the potential to alter the cosmic history of their host galaxies. These signatures are only seen in ~10% of optically selected quasars, although the fraction significantly increases in IR and radio selected samples. A proven physical explanation for this observed fraction has yet to be found, along with a determination of why this fraction increases at radio wavelengths. We present the largest sample of radio matched BALQSOs using the LOFAR Two-metre Sky Survey Data Release 2 and employ it to investigate radio properties of BALQSOs. Within the DR2 footprint, there are 3537 BALQSOs from Sloan Digital Sky Survey DR12 with continuum signal to noise >5. We find radio-detections for 1108 BALQSOs, with an important sub-population of 120 LoBALs, an unprecedented sample size for radio matched BALQSOs given the LoTSS sky coverage to date. BALQSOs are a radio-quiet population that show an increase of $\times 1.50$ radio-detection fraction compared to non-BALQSOs. LoBALs show an increase of $\times 2.22$ that of non-BALQSO quasars. We show that this detection fraction correlates with wind-strength, reddening and C_{IV} emission properties of BALQSOs and that these features may be connected, although no single property can fully explain the enhanced radio detection fraction. We create composite spectra for sub-classes of BALQSOs based on wind strength and colour, finding differences in the absorption profiles of radio-detected and radio-undetected sources, particularly for LoBALs. Overall, we favour a wind-ISM interaction explanation for the increased radio-detection fraction of BALQSOs. △ Less

Submitted 20 July, 2022; originally announced July 2022.

Comments: 16 pages, 14 figures

Evidence for a moderate spin from X-ray reflection of the high-mass supermassive black hole in the cluster-hosted quasar H1821+643

Authors: Júlia Sisk-Reynés, Christopher S. Reynolds, James H. Matthews, Robyn N. Smith

Abstract: We present an analysis of deep $Chandra$ Low-Energy and High-Energy Transmission Grating archival observations of the extraordinarily luminous radio-quiet quasar H1821+643, hosted by a rich and massive cool-core cluster at redshift $z=0.3$. These datasets provide high-resolution spectra of the AGN at two epochs, free from contamination by the intracluster medium and from the effects of photon pile… ▽ More We present an analysis of deep $Chandra$ Low-Energy and High-Energy Transmission Grating archival observations of the extraordinarily luminous radio-quiet quasar H1821+643, hosted by a rich and massive cool-core cluster at redshift $z=0.3$. These datasets provide high-resolution spectra of the AGN at two epochs, free from contamination by the intracluster medium and from the effects of photon pile-up, providing a sensitive probe of the iron-$K$ band. At both epochs, the spectrum is well described by a power-law continuum plus X-ray reflection from both the inner accretion disc and cold, slowly-moving distant matter. Adopting this framework, we proceed to examine the properties of the inner disc and the black hole spin. Using Markov chain Monte Carlo (MCMC) methods, we combine constraints from the two epochs assuming that the black hole spin, inner disc inclination, and inner disc iron abundance are invariant. The black hole spin is found to be modest, with a 90$\%$ credible range of ${a}^{*}=0.62^{+0.22}_{-0.37}$; and, with a mass $M_\mathrm{BH}$ in the range $\log (M_\mathrm{BH}/M_\odot)\sim 9.2-10.5$, this is the most massive black hole candidate for which a well-defined spin constraint has yet been obtained. The modest spin of this black hole supports previous suggestions that the most massive black holes may grow via incoherent or chaotic accretion and/or SMBH-SMBH mergers. △ Less

Submitted 30 August, 2022; v1 submitted 25 May, 2022; originally announced May 2022.

Comments: 13 pages, 8 figures, 3 tables. Accepted for publication in MNRAS. Amended to include two missing references

Journal ref: MNRAS, 514, 2568 - 2580 (2022)

The X-ray Disk/Wind Degeneracy in AGN

Authors: M. L. Parker, G. A. Matzeu, J. H. Matthews, M. J. Middleton, T. Dauser, J. Jiang, A. M. Joyce

Abstract: Relativistic Fe K emission lines from accretion disks and from disk winds encode key information about black holes, and their accretion and feedback mechanisms. We show that these two processes can in principle produce indistinguishable line profiles, such that they cannot be disentangled spectrally. We argue that it is likely that in many cases both processes contribute to the net line profile, a… ▽ More Relativistic Fe K emission lines from accretion disks and from disk winds encode key information about black holes, and their accretion and feedback mechanisms. We show that these two processes can in principle produce indistinguishable line profiles, such that they cannot be disentangled spectrally. We argue that it is likely that in many cases both processes contribute to the net line profile, and their relative contributions cannot be constrained purely by Fe K spectroscopy. In almost all studies of Fe K emission to date, a single process (either disk reflection or wind Compton scattering) is assumed to dominate the total line profile. We demonstrate that fitting a single process emission model (pure reflection or pure wind) to a hybrid line profile results in large systematic biases in the estimates of key parameters, such as mass outflow rate and spin. We discuss various strategies to mitigate this effect, such as including high energy data covering the Compton hump, and the implications for future X-ray missions. △ Less

Submitted 28 March, 2022; originally announced March 2022.

Comments: 22 pages, 17 figures, accepted for publication in MNRAS

How do Magnetic Field Models Affect Astrophysical Limits on Light Axion-like Particles? An X-ray Case Study with NGC 1275

Authors: James H. Matthews, Christopher S. Reynolds, M. C. David Marsh, Júlia Sisk-Reynés, Payton E. Rodman

Abstract: Axion-like particles (ALPs) are a well-motivated extension to the standard model of particle physics, and X-ray observations of cluster-hosted AGN currently place the most stringent constraints on the ALP coupling to electromagnetism, $g_{a γ}$, for very light ALPs ($m_a\lesssim10^{-11}$ eV). We revisit limits obtained by Reynolds et al. (2020) using Chandra X-ray grating spectroscopy of NGC 1275,… ▽ More Axion-like particles (ALPs) are a well-motivated extension to the standard model of particle physics, and X-ray observations of cluster-hosted AGN currently place the most stringent constraints on the ALP coupling to electromagnetism, $g_{a γ}$, for very light ALPs ($m_a\lesssim10^{-11}$ eV). We revisit limits obtained by Reynolds et al. (2020) using Chandra X-ray grating spectroscopy of NGC 1275, the central AGN in the Perseus cluster, examining the impact of the X-ray spectral model and magnetic field model. We also present a new publicly available code, ALPro, which we use to solve the ALP propagation problem. We discuss evidence for turbulent magnetic fields in Perseus and show that it can be important to resolve the magnetic field structure on scales below the coherence length. We re-analyse the NGC 1275 X-ray spectra using an improved data reduction and baseline spectral model. We find the limits are insensitive to whether a partially covering absorber is used in the fits. At low $m_a$ ($m_a\lesssim10^{-13}$ eV), we find marginally weaker limits on $g_{a γ}$ (by $0.1-0.3$ dex) with different magnetic field models, compared to Model B from Reynolds et al. (2020). A Gaussian random field (GRF) model designed to mimic $\sim50$ kpc scale coherent structures also results in only slightly weaker limits. We conclude that the existing Model B limits are robust assuming that $β_{\rm pl}\approx100$, and are insensitive to whether cell-based or GRF methods are used. However, astrophysical uncertainties regarding the strength and structure of cluster magnetic fields persist, motivating high sensitivity RM observations and tighter constraints on the radial profile of $β_{\rm pl}$. △ Less

Submitted 17 February, 2022; originally announced February 2022.

Comments: 17 pages, 9 figures, 2 tables. Accepted for publication in The Astrophysical Journal. ALPro code available at https://github.com/jhmatthews/alpro with accompanying documentation, as well as notebooks and data for reproducing the main results

A Decade of Black-Hole X-ray Binary Transients

Authors: Philip A Charles, David A H Buckley, Enrico Kotze, Marissa M Kotze, Jessymol K Thomas, Poshak Gandhi, John A Paice, Jean-Pierre Lasota, James H Matthews, James F Steiner

Abstract: The last decade has seen a significant gain in both space and ground-based monitoring capabilities, producing vastly better coverage of BH X-ray binaries during their (rare) transient events. This interval included two of the three brightest X-ray outbursts ever observed, namely V404 Cyg in 2015, and MAXI J1820+070 in 2018, as well as the outburst of Swift J1357.2-0933, the first such system to sh… ▽ More The last decade has seen a significant gain in both space and ground-based monitoring capabilities, producing vastly better coverage of BH X-ray binaries during their (rare) transient events. This interval included two of the three brightest X-ray outbursts ever observed, namely V404 Cyg in 2015, and MAXI J1820+070 in 2018, as well as the outburst of Swift J1357.2-0933, the first such system to show variable period optical dipping. There are now superb multi-wavelength archives of these outbursts, both photometric and spectroscopic, that show substantial outflows in the form of jets and disc winds, and X-ray spectroscopy/timing that reveals how the inner accretion disc evolves. The ground-based AAVSO optical monitoring of the MAXI J1820+070 event was the most extensive ever obtained, revealing periodic variations that evolved as it approached its state transition. These modulations were of an amplitude never seen before, and suggested the development of an irradiation-driven disc warp that persisted through the transition. All these results have demonstrated the power of extensive multi-wavelength photometric and spectroscopic monitoring on all time-scales. △ Less

Submitted 27 January, 2022; originally announced January 2022.

Comments: 10 Pages, 20 Figures

Optical line spectra of tidal disruption events from reprocessing in optically thick outflows

Authors: Edward J. Parkinson, Christian Knigge, James H. Matthews, Knox S. Long, Nick Higginbottom, Stuart A. Sim, Samuel W. Mangham

Abstract: A significant number of tidal disruption events (TDEs) radiate primarily at optical and ultraviolet (UV) wavelengths, with only weak soft X-ray components. One model for this optical excess proposes that thermal X-ray emission from a compact accretion disc is reprocessed to longer wavelengths by an optically thick envelope. Here, we explore this reprocessing scenario in the context of an optically… ▽ More A significant number of tidal disruption events (TDEs) radiate primarily at optical and ultraviolet (UV) wavelengths, with only weak soft X-ray components. One model for this optical excess proposes that thermal X-ray emission from a compact accretion disc is reprocessed to longer wavelengths by an optically thick envelope. Here, we explore this reprocessing scenario in the context of an optically thick accretion disc wind. Using state-of-the-art Monte Carlo radiative transfer and ionization software, we produce synthetic UV and optical spectra for wind and disc-hosting TDEs. Our models are inspired by observations, spanning a realistic range of accretion rates and wind kinematics. We find that such outflows can efficiently reprocess the disc emission and produce the broad Balmer and helium recombination features commonly seen in TDEs and exhibit asymmetric red wings. Moreover, the characteristic colour temperature of the reprocessed spectral energy distribution (SED) is much lower than that of the accretion disc. We show explicitly how changes in black hole mass, accretion rate and wind properties affect the observed broadband SED and line spectrum. In general, slower, denser winds tend to reprocess more radiation and produce stronger Balmer emission. Most of the outflows we consider are too highly ionized to produce UV absorption features, but this is sensitive to the input SED. For example, truncating the inner disc at just 4 $R_{ISCO}$ lowers the wind ionization state sufficiently to produce UV absorption features for sight lines looking into the wind △ Less

Submitted 5 January, 2022; originally announced January 2022.

Comments: Accepted to MNRAS. 19 pages, 11 figures

EuCAPT White Paper: Opportunities and Challenges for Theoretical Astroparticle Physics in the Next Decade

Authors: R. Alves Batista, M. A. Amin, G. Barenboim, N. Bartolo, D. Baumann, A. Bauswein, E. Bellini, D. Benisty, G. Bertone, P. Blasi, C. G. Böhmer, Ž. Bošnjak, T. Bringmann, C. Burrage, M. Bustamante, J. Calderón Bustillo, C. T. Byrnes, F. Calore, R. Catena, D. G. Cerdeño, S. S. Cerri, M. Chianese, K. Clough, A. Cole, P. Coloma , et al. (112 additional authors not shown)

Abstract: Astroparticle physics is undergoing a profound transformation, due to a series of extraordinary new results, such as the discovery of high-energy cosmic neutrinos with IceCube, the direct detection of gravitational waves with LIGO and Virgo, and many others. This white paper is the result of a collaborative effort that involved hundreds of theoretical astroparticle physicists and cosmologists, und… ▽ More Astroparticle physics is undergoing a profound transformation, due to a series of extraordinary new results, such as the discovery of high-energy cosmic neutrinos with IceCube, the direct detection of gravitational waves with LIGO and Virgo, and many others. This white paper is the result of a collaborative effort that involved hundreds of theoretical astroparticle physicists and cosmologists, under the coordination of the European Consortium for Astroparticle Theory (EuCAPT). Addressed to the whole astroparticle physics community, it explores upcoming theoretical opportunities and challenges for our field of research, with particular emphasis on the possible synergies among different subfields, and the prospects for solving the most fundamental open questions with multi-messenger observations. △ Less

Submitted 19 October, 2021; originally announced October 2021.

Comments: White paper of the European Consortium for Astroparticle Theory (EuCAPT). 135 authors, 400 endorsers, 133 pages, 1382 references

New constraints on light Axion-Like Particles using Chandra Transmission Grating Spectroscopy of the powerful cluster-hosted quasar H1821+643

Authors: Júlia Sisk Reynés, James H. Matthews, Christopher S. Reynolds, Helen R. Russell, Robyn N. Smith, M. C. David Marsh

Abstract: Axion-Like Particles (ALPs) are predicted by several Beyond the Standard Model theories, in particular, string theory. In the presence of an external magnetic field perpendicular to the direction of propagation, ALPs can couple to photons. Therefore, if an X-ray source is viewed through a magnetised plasma, such as a luminous quasar in a galaxy cluster, we may expect spectral distortions that are… ▽ More Axion-Like Particles (ALPs) are predicted by several Beyond the Standard Model theories, in particular, string theory. In the presence of an external magnetic field perpendicular to the direction of propagation, ALPs can couple to photons. Therefore, if an X-ray source is viewed through a magnetised plasma, such as a luminous quasar in a galaxy cluster, we may expect spectral distortions that are well described by photon-ALP oscillations. We present a $571 \ \mathrm{ks}$ combined High and Low Energy Transmission Grating (HETG/LETG) Chandra observation of the powerful radio-quiet quasar H1821+643, hosted by a cool-core cluster at redshift $0.3$. The spectrum is well described by a double power-law continuum and broad$+$narrow iron line emission typical of type-1 Active Galactic Nuclei (AGN), with remaining spectral features $< 2.5\%$. Using a cell-based approach to describe the turbulent cluster magnetic field, we compare our spectrum with photon-ALP mixing curves for 500 field realisations assuming that the thermal-to-magnetic pressure ratio remains constant up to the virial radius. At $99.7\%$ credibility and taking $β= 100$, we exclude all couplings $g_\mathrm{aγ} > 6.3 \times 10^{-13} \ {\mathrm{GeV}}^{-1}$ for most ALP masses $< 10^{-12} \ \mathrm{eV}$. Our results are moderately more sensitive to constraining ALPs than the best previous result from Chandra observations of the Perseus cluster, albeit with a less constrained field model. We reflect on the promising future of ALP studies with bright AGN embedded in rich clusters, especially with the upcoming Athena mission. △ Less

Submitted 25 November, 2021; v1 submitted 7 September, 2021; originally announced September 2021.

Comments: 14 pages, 7 figures. Accepted for publication in MNRAS on 24th November 2021. Comments welcome

The Fourier formalism for relativistic axion-photon conversion, with astrophysical applications

Authors: M. C. David Marsh, James H. Matthews, Christopher Reynolds, Pierluca Carenza

Abstract: We study the weak mixing of photons and relativistic axion-like particles (axions) in plasmas with background magnetic fields, ${\bf B}$. We show that, to leading order in the axion-photon coupling, the conversion probability, $P_{γ\to a}$, is given by the one-dimensional power spectrum of the magnetic field components perpendicular to the particle trajectory. Equivalently, we express… ▽ More We study the weak mixing of photons and relativistic axion-like particles (axions) in plasmas with background magnetic fields, ${\bf B}$. We show that, to leading order in the axion-photon coupling, the conversion probability, $P_{γ\to a}$, is given by the one-dimensional power spectrum of the magnetic field components perpendicular to the particle trajectory. Equivalently, we express $P_{γ\to a}$ as the Fourier transform of the magnetic field autocorrelation function, and establish a dictionary between properties of the real-space magnetic field and the energy-dependent conversion probability. For axions more massive than the plasma frequency, ($m_a>ω_{\rm pl}$), we use this formalism to analytically solve the problem of perturbative axion-photon mixing in a general magnetic field. In the general case where $ω_{\rm pl}/m_a$ varies arbitrarily along the trajectory, we show that a naive application of the standard formalism for 'resonant' conversion can give highly inaccurate results, and that a careful calculation generically gives non-resonant contributions at least as large as the resonant contribution. Furthermore, we demonstrate how techniques based on the Fast Fourier Transform provide a new, highly efficient numerical method for calculating axion-photon mixing. We briefly discuss magnetic field modelling in galaxy clusters in the light of our results and argue, in particular, that a recently proposed 'regular' model used for studying axion-photon mixing (specifically applied to the Perseus cluster) is inconsistent with observations. Our formalism suggest new methods to search for imprints of axions, and will be important for spectrographs with percent level sensitivity, which includes existing X-ray observations by Chandra as well as the upcoming Athena mission. △ Less

Submitted 3 February, 2022; v1 submitted 16 July, 2021; originally announced July 2021.

Comments: 39 pages, 14 figures. v3: matches published version

Probing the Wind Component of Radio Emission in Luminous High-Redshift Quasars

Authors: Gordon T. Richards, Trevor V. McCaffrey, Amy Kimball, Amy L. Rankine, James H. Matthews, Paul C. Hewett, Angelica B. Rivera

Abstract: We discuss a probe of the contribution of wind-related shocks to the radio emission in otherwise radio-quiet quasars. Given 1) the non-linear correlation between UV and X-ray luminosity in quasars, 2) that such correlation leads to higher likelihood of radiation-line-driven winds in more luminous quasars, and 3) that luminous quasars are more abundant at high redshift, deep radio observations of h… ▽ More We discuss a probe of the contribution of wind-related shocks to the radio emission in otherwise radio-quiet quasars. Given 1) the non-linear correlation between UV and X-ray luminosity in quasars, 2) that such correlation leads to higher likelihood of radiation-line-driven winds in more luminous quasars, and 3) that luminous quasars are more abundant at high redshift, deep radio observations of high-redshift quasars are needed to probe potential contributions from accretion disk winds. We target a sample of 50 $z\simeq 1.65$ color-selected quasars that span the range of expected accretion disk wind properties as traced by broad CIV emission. 3-GHz observations with the Very Large Array to an rms of $\approx10μ$Jy beam$^{-1}$ probe to star formation rates of $\approx400\,M_{\rm Sun}\,{\rm yr}^{-1}$, leading to 22 detections. Supplementing these pointed observations are survey data of 388 sources from the LOFAR Two-metre Sky Survey Data Release 1 that reach comparable depth (for a typical radio spectral index), where 123 sources are detected. These combined observations reveal a radio detection fraction that is a non-linear function of \civ\ emission-line properties and suggest that the data may require multiple origins of radio emission in radio-quiet quasars. We find evidence for radio emission from weak jets or coronae in radio-quiet quasars with low Eddingtion ratios, with either (or both) star formation and accretion disk winds playing an important role in optically luminous quasars and correlated with increasing Eddington ratio. Additional pointed radio observations are needed to fully establish the nature of radio emission in radio-quiet quasars. △ Less

Submitted 14 June, 2021; originally announced June 2021.

Comments: 21 pages, 7 figures, submitted to AAS journals

Particle acceleration in radio galaxies with flickering jets: GeV electrons to ultrahigh energy cosmic rays

Authors: James H. Matthews, Andrew M. Taylor

Abstract: Variability is a general property of accretion discs and their associated jets. We introduce a semi-analytic model for particle acceleration and radio jet/lobe evolution and explore the effect of Myr timescale jet variability on the particles accelerated by an AGN jet. Our work is motivated by the need for local powerful ultrahigh energy cosmic ray (UHECR) sources and evidence for variability in A… ▽ More Variability is a general property of accretion discs and their associated jets. We introduce a semi-analytic model for particle acceleration and radio jet/lobe evolution and explore the effect of Myr timescale jet variability on the particles accelerated by an AGN jet. Our work is motivated by the need for local powerful ultrahigh energy cosmic ray (UHECR) sources and evidence for variability in AGN and radio galaxies. Our main results are: i) UHECR and nonthermal radiative luminosities track the jet power but with a response set by the escape and cooling times, respectively; ii) jet variability produces structure in the electron, synchrotron and UHECR spectra that deviates from that produced for a constant jet power - in particular, spectral hardening features may be signatures of variability; iii) the cutoff in the integrated CR spectrum is stretched out due to the variation in jet power (and, consequently, maximum CR energy). The resulting spectrum is the convolution of the jet power distribution and the source term. We derive an approximate form for a log-normal distribution of powers; iv) we introduce the idea of $\sim 10$ GeV 'proxy electrons' that are cooling at the same rate that UHECRs of rigidity 10 EV are escaping from the source, and determine the corresponding photon frequencies that probe escaping UHECRs. Our results demonstrate the link between the history of an astrophysical particle accelerator and its particle contents, nonthermal emission and UHECR spectrum, with consequences for observations of radio galaxies and UHECR source models. △ Less

Submitted 11 March, 2021; originally announced March 2021.

Comments: 17 pages, 11 figures. Accepted to MNRAS

Placing LOFAR-detected quasars in CIV emission space: implications for winds, jets and star formation

Authors: Amy L. Rankine, James H. Matthews, Paul C. Hewett, Manda Banerji, Leah K. Morabito, Gordon T. Richards

Abstract: We present an investigation of the low-frequency radio and ultraviolet properties of a sample of $\simeq$10,500 quasars from the Sloan Digital Sky Survey Data Release 14, observed as part of the first data release of the Low-Frequency-Array Two-metre Sky Survey. The quasars have redshifts $1.5 < z < 3.5$ and luminosities $44.6 < \log(L_{\text{bol}}/\text{erg s}^{-1}) < 47.2$. We employ ultraviolet… ▽ More We present an investigation of the low-frequency radio and ultraviolet properties of a sample of $\simeq$10,500 quasars from the Sloan Digital Sky Survey Data Release 14, observed as part of the first data release of the Low-Frequency-Array Two-metre Sky Survey. The quasars have redshifts $1.5 < z < 3.5$ and luminosities $44.6 < \log(L_{\text{bol}}/\text{erg s}^{-1}) < 47.2$. We employ ultraviolet spectral reconstructions based on an independent component analysis to parametrize the CIV$λ$1549-emission line that is used to infer the strength of accretion disc winds, and the HeII$λ$1640 line, an indicator of the soft X-ray flux. We find that radio-detected quasars are found in the same region of CIV blueshift versus equivalent-width space as radio-undetected quasars, but that the loudest, most luminous and largest radio sources exist preferentially at low CIV blueshifts. Additionally, the radio-detection fraction increases with blueshift whereas the radio-loud fraction decreases. In the radio-quiet population, we observe a range of HeII equivalent widths as well as a Baldwin effect with bolometric luminosity, whilst the radio-loud population has mostly strong HeII, consistent with a stronger soft X-ray flux. The presence of strong HeII is a necessary but not sufficient condition to detect radio-loud emission suggesting some degree of stochasticity in jet formation. Using energetic arguments and Monte Carlo simulations, we explore the plausibility of winds, compact jets and star formation as sources of the radio quiet emission, ruling out none. The existence of quasars with similar ultraviolet properties but differing radio properties suggests, perhaps, that the radio and ultraviolet emission is tracing activity occurring on different timescales. △ Less

Submitted 2 March, 2021; v1 submitted 29 January, 2021; originally announced January 2021.

Comments: 13 pages, 10 figures, plus appendices. Minor changes to match published version

Accretion Disc Winds in Tidal Disruption Events: Ultraviolet Spectral Lines as Orientation Indicators

Authors: Edward J. Parkinson, Christian Knigge, Knox S. Long, James H. Matthews, Nick Higginbottom, Stuart A. Sim, Henrietta A. Hewitt

Abstract: Some tidal disruption events (TDEs) exhibit blueshifted broad absorption lines (BALs) in their rest-frame ultraviolet (UV) spectra, while others display broad emission lines (BELs). Similar phenomenology is observed in quasars and accreting white dwarfs, where it can be interpreted as an orientation effect associated with line formation in an accretion disc wind.We propose and explore a similar un… ▽ More Some tidal disruption events (TDEs) exhibit blueshifted broad absorption lines (BALs) in their rest-frame ultraviolet (UV) spectra, while others display broad emission lines (BELs). Similar phenomenology is observed in quasars and accreting white dwarfs, where it can be interpreted as an orientation effect associated with line formation in an accretion disc wind.We propose and explore a similar unification scheme for TDEs. We present synthetic UV spectra for disc and wind-hosting TDEs, produced by a state-of-the-art Monte Carlo ionization and radiative transfer code. Our models cover a wide range of disc wind geometries and kinematics. Such winds naturally reproduce both BALs and BELs. In general, sight lines looking into the wind cone preferentially produce BALs, while other orientations preferentially produce BELs. We also study the effect of wind clumping and CNO-processed abundances on the observed spectra. Clumpy winds tend to produce stronger UV emission and absorption lines, because clumping increases both the emission measure and the abundances of the relevant ionic species, the latter by reducing the ionization state of the outflow. The main effect of adopting CNO-processed abundances is a weakening of C~{\sc iv}~1550~Å~ and an enhancement of N \textsc{v}~1240~Å~ in the spectra. We conclude that line formation in an accretion disc wind is a promising mechanism for explaining the diverse UV spectra of TDEs. If this is correct, the relative number of BAL and BEL TDEs can be used to estimate the covering factor of the outflow. The models in this work are publicly available online and upon request. △ Less

Submitted 16 April, 2020; originally announced April 2020.

Comments: 16 pages, 9 figures

Thermal and radiation driving can produce observable disk winds in hard-state X-ray binaries

Authors: Nick Higginbottom, Christian Knigge, Stuart A. Sim, Knox S. Long, James H. Matthews, Henrietta A. Hewitt, Edward J. Parkinson, Sam W. Mangham

Abstract: X-ray signatures of outflowing gas have been detected in several accreting black-hole binaries, always in the soft state. A key question raised by these observations is whether these winds might also exist in the hard state. Here, we carry out the first full-frequency radiation hydrodynamic simulations of luminous ($\rm{L = 0.5 \, L_{\mathrm{Edd}}}$) black-hole X-ray binary systems in both the har… ▽ More X-ray signatures of outflowing gas have been detected in several accreting black-hole binaries, always in the soft state. A key question raised by these observations is whether these winds might also exist in the hard state. Here, we carry out the first full-frequency radiation hydrodynamic simulations of luminous ($\rm{L = 0.5 \, L_{\mathrm{Edd}}}$) black-hole X-ray binary systems in both the hard and the soft state, with realistic spectral energy distributions (SEDs). Our simulations are designed to describe X-ray transients near the peak of their outburst, just before and after the hard-to-soft state transition. At these luminosities, it is essential to include radiation driving, and we include not only electron scattering, but also photoelectric and line interactions. We find powerful outflows with $\rm{\dot{M}_{wind} \simeq 2 \,\dot{M}_{acc}}$ are driven by thermal and radiation pressure in both hard and soft states. The hard-state wind is significantly faster and carries approximately 20 times as much kinetic energy as the soft-state wind. However, in the hard state the wind is more ionized, and so weaker X-ray absorption lines are seen over a narrower range of viewing angles. Nevertheless, for inclinations $\gtrsim 80^{\circ}$, blue-shifted wind-formed Fe XXV and Fe XXVI features should be observable even in the hard state. Given that the data required to detect these lines currently exist for only a single system in a {\em luminous} hard state -- the peculiar GRS~1915+105 -- we urge the acquisition of new observations to test this prediction. The new generation of X-ray spectrometers should be able to resolve the velocity structure. △ Less

Submitted 27 January, 2020; v1 submitted 21 January, 2020; originally announced January 2020.

Comments: Accepted for publication by MNRAS

Stratified disc wind models for the AGN broad-line region: ultraviolet, optical and X-ray properties

Authors: James H. Matthews, Christian Knigge, Nick Higginbottom, Knox S. Long, Stuart A. Sim, Samuel W. Mangham, Edward J. Parkinson, Henrietta A. Hewitt

Abstract: The origin, geometry and kinematics of the broad line region (BLR) gas in quasars and active galactic nuclei (AGN) are uncertain. We demonstrate that clumpy biconical disc winds illuminated by an AGN continuum can produce BLR-like spectra. We first use a simple toy model to illustrate that disc winds make quite good BLR candidates, because they are self-shielded flows and can cover a large portion… ▽ More The origin, geometry and kinematics of the broad line region (BLR) gas in quasars and active galactic nuclei (AGN) are uncertain. We demonstrate that clumpy biconical disc winds illuminated by an AGN continuum can produce BLR-like spectra. We first use a simple toy model to illustrate that disc winds make quite good BLR candidates, because they are self-shielded flows and can cover a large portion of the ionizing flux-density ($φ_H$-$n_H$) plane. We then conduct Monte Carlo radiative transfer and photoionization calculations, which fully account for self-shielding and multiple scattering in a non-spherical geometry. The emergent model spectra show broad emission lines with equivalent widths and line ratios comparable to those observed in AGN, provided that the wind has a volume filling factor of $f_V\lesssim0.1$. Similar emission line spectra are produced for a variety of wind geometries (polar or equatorial) and for launch radii that differ by an order of magnitude. The line emission arises almost exclusively from plasma travelling below the escape velocity, implying that `failed winds' are important BLR candidates. The behaviour of a line-emitting wind (and possibly any `smooth flow' BLR model) is similar to that of the locally optimally-emitting cloud (LOC) model originally proposed by Baldwin et al (1995), except that the gradients in ionization state and temperature are large-scale and continuous, rather than within or between distinct clouds. Our models also produce UV absorption lines and X-ray absorption features, and the stratified ionization structure can partially explain the different classes of broad absorption line quasars. △ Less

Submitted 10 January, 2020; originally announced January 2020.

Comments: Accepted to MNRAS. 20 pages, 13 figures plus appendix. Supplementary plots, some simulation data and demo notebooks can be found at https://github.com/jhmatthews/windy-blr-2020

Hot, dense He II outflows during the 2017 outburst of the X-ray transient Swift J1357.2-0933

Authors: Phil Charles, James H. Matthews, David A. H. Buckley, Poshak Gandhi, Enrico Kotze, John Paice

Abstract: Time-resolved SALT spectra of the short-period, dipping X-ray transient, Swift J1357.2-0933, during its 2017 outburst has revealed broad Balmer and HeII4686 absorption features, blue-shifted by ~600 km/s. Remarkably these features are also variable on the ~500s dipping period, indicating their likely association with structure in the inner accretion disc. We interpret this as arising in a dense, h… ▽ More Time-resolved SALT spectra of the short-period, dipping X-ray transient, Swift J1357.2-0933, during its 2017 outburst has revealed broad Balmer and HeII4686 absorption features, blue-shifted by ~600 km/s. Remarkably these features are also variable on the ~500s dipping period, indicating their likely association with structure in the inner accretion disc. We interpret this as arising in a dense, hot (>~30,000K) outflowing wind seen at very high inclination, and draw comparisons with other accretion disc corona sources. We argue against previous distance estimates of 1.5 kpc and favour a value >~6 kpc, implying an X-ray luminosity L_X>~4x10^{36} erg/s$. Hence it is not a very faint X-ray transient. Our preliminary 1D Monte-Carlo radiative transfer and photoionization calculations support this interpretation, as they imply a high intrinsic L_X, a column density N_H>~10^{24} cm^{-2} and a low covering factor for the wind. Our study shows that Swift J1357.2-0933 is truly remarkable amongst the cohort of luminous, galactic X-ray binaries, showing the first example of He II absorption, the first (and only) variable dip period and is possibly the first black hole 'accretion disc corona' candidate. △ Less

Submitted 20 August, 2019; v1 submitted 1 August, 2019; originally announced August 2019.

Comments: MNRAS Letters in press. Updated to match published version

Do Reverberation Mapping Analyses Provide an Accurate Picture of the Broad Line Region?

Authors: S. W. Mangham, C. Knigge, P. Williams, Keith Horne, A. Pancoast, J. H. Matthews, K. S. Long, S. A. Sim, N. Higginbottom

Abstract: Reverberation mapping (RM) is a powerful approach for determining the nature of the broad-line region (BLR) in active galactic nuclei. However, inferring physical BLR properties from an observed spectroscopic time series is a difficult inverse problem. Here, we present a blind test of two widely used RM methods: MEMEcho (developed by Horne) and CARAMEL (developed by Pancoast and collaborators). Th… ▽ More Reverberation mapping (RM) is a powerful approach for determining the nature of the broad-line region (BLR) in active galactic nuclei. However, inferring physical BLR properties from an observed spectroscopic time series is a difficult inverse problem. Here, we present a blind test of two widely used RM methods: MEMEcho (developed by Horne) and CARAMEL (developed by Pancoast and collaborators). The test data are simulated spectroscopic time series that track the H$α$ emission line response to an empirical continuum light curve. The underlying BLR model is a rotating, biconical accretion disc wind, and the synthetic spectra are generated via self-consistent ionization and radiative transfer simulations. We generate two mock data sets, representing Seyfert galaxies and QSOs. The Seyfert model produces a largely *negative* response, which neither method can recover. However, both fail $``gracefully''$, neither generating spurious results. For the QSO model both CARAMEL and expert interpretation of MEMEcho's output both capture the broadly annular, rotation-dominated nature of the line-forming region, though MEMEcho analysis overestimates its size by 50%, but CARAMEL is unable to distinguish between additional inflow and outflow components. Despite fitting individual spectra well, the CARAMEL velocity-delay maps and RMS line profiles are strongly inconsistent with the input data. Finally, since the H$α$ line-forming region is rotation dominated, neither method recovers the disc wind nature of the underlying BLR model. Thus considerable care is required when interpreting the results of RM analyses in terms of physical models. △ Less

Submitted 26 June, 2019; originally announced June 2019.

Comments: 22 pages, 24 figures

Report number: MN-19-0268-MJ.R1

Cosmic ray acceleration to ultrahigh energy in radio galaxies

Authors: James H. Matthews, Anthony R. Bell, Anabella T. Araudo, Katherine M. Blundell

Abstract: The origin of ultrahigh energy cosmic rays (UHECRs) is an open question. In this proceeding, we first review the general physical requirements that a source must meet for acceleration to 10-100 EeV, including the consideration that the shock is not highly relativistic. We show that shocks in the backflows of radio galaxies can meet these requirements. We discuss a model in which giant-lobed radio… ▽ More The origin of ultrahigh energy cosmic rays (UHECRs) is an open question. In this proceeding, we first review the general physical requirements that a source must meet for acceleration to 10-100 EeV, including the consideration that the shock is not highly relativistic. We show that shocks in the backflows of radio galaxies can meet these requirements. We discuss a model in which giant-lobed radio galaxies such as Centaurus A and Fornax A act as slowly-leaking UHECR reservoirs, with the UHECRs being accelerated during a more powerful past episode. We also show that Centaurus A, Fornax A and other radio galaxies may explain the observed anisotropies in data from the Pierre Auger Observatory, before examining some of the difficulties in associating UHECR anisotropies with astrophysical sources. △ Less

Submitted 27 February, 2019; originally announced February 2019.

Comments: 6 pages, 4 figures. Proceedings of UHECR 2018, 8-12 October 2018, Paris, France

The luminosity dependence of thermally-driven disc winds in low-mass X-ray binaries

Authors: Nick Higginbottom, Christian Knigge, Knox S. Long, James H. Matthews, Edward J. Parkinson

Abstract: We have carried out radiation-hydrodynamic simulations of thermally-driven accretion disc winds in low-mass X-ray binaries. Our main goal is to study the luminosity dependence of these outflows and compare with observations. The simulations span the range $\rm{0.04 \leq L_{acc}/L_{Edd} \leq 1.0}$ and therefore cover most of the parameter space in which disc winds have been observed. Using a detail… ▽ More We have carried out radiation-hydrodynamic simulations of thermally-driven accretion disc winds in low-mass X-ray binaries. Our main goal is to study the luminosity dependence of these outflows and compare with observations. The simulations span the range $\rm{0.04 \leq L_{acc}/L_{Edd} \leq 1.0}$ and therefore cover most of the parameter space in which disc winds have been observed. Using a detailed Monte Carlo treatment of ionization and radiative transfer, we confirm two key results found in earlier simulations that were carried out in the optically thin limit: (i) the wind velocity -- and hence the maximum blueshift seen in wind-formed absorption lines -- increases with luminosity; (ii) the large-scale wind geometry is quasi-spherical, but observable absorption features are preferentially produced along high-column equatorial sightlines. In addition, we find that (iii) the wind efficiency always remains approximately constant at $\rm{\dot{M}_{wind}/\dot{M}_{acc} \simeq 2}$, a behaviour that is consistent with observations. We also present synthetic Fe XXV and Fe XXVI absorption line profiles for our simulated disc winds in order to illustrate the observational implications of our results. △ Less

Submitted 28 January, 2019; originally announced January 2019.

Comments: Accepted for publication in MNRAS

Disc wind models for FU Ori objects

Authors: Kelly Milliner, James H. Matthews, Knox S. Long, Lee Hartmann

Abstract: We present disc wind models aimed at reproducing the main features of the strong Na I resonance line P-Cygni profiles in the rapidly-accreting pre-main sequence FU Ori objects. We conducted Monte Carlo radiative transfer simulations for a standard magnetocentrifugally driven wind (MHD) model and our own "Genwind" models, which allows for a more flexible wind parameterisation. We find that the fidu… ▽ More We present disc wind models aimed at reproducing the main features of the strong Na I resonance line P-Cygni profiles in the rapidly-accreting pre-main sequence FU Ori objects. We conducted Monte Carlo radiative transfer simulations for a standard magnetocentrifugally driven wind (MHD) model and our own "Genwind" models, which allows for a more flexible wind parameterisation. We find that the fiducial MHD wind and similar Genwind models, which have flows emerging outward from the inner disc edge, and thus have polar cavities with no absorbing gas, cannot reproduce the deep, wide Na I absorption lines in FU Ori objects viewed at low inclination. We find that it is necessary to include an "inner wind" to fill this polar cavity to reproduce observations. In addition, our models assuming pure scattering source functions in the Sobolev approximation at intermediate viewing angles ($30^{\circ} \lesssim i \lesssim 60^{\circ}$) do not yield sufficiently deep line profiles. Assuming complete absorption yields better agreement with observations, but simple estimates strongly suggest that pure scattering should be a much better approximation. The discrepancy may indicate that the Sobolev approximation is not applicable, possibly due to turbulence or non-monotonic velocity fields; there is some observational evidence for the latter. Our results provide guidance for future attempts to constrain FU Ori wind properties using full MHD wind simulations, by pointing to the importance of the boundary conditions necessary to give rise to an inner wind, and by suggesting that the winds must be turbulent to produce sufficiently deep line profiles. △ Less

Submitted 29 November, 2018; originally announced November 2018.

Comments: 12 pages, 17 figures, accepted for publication in MNRAS

The origin of radio emission in broad absorption line quasars: Results from the LOFAR Two-metre Sky Survey

Authors: L. K. Morabito, J. H. Matthews, P. N. Best, G. Gürkan, M. J. Jarvis, I. Prandoni, K. J. Duncan, M. J. Hardcastle, M. Kunert-Bajraszewska, A. P. Mechev, S. Mooney, J. Sabater, H. J. A. Röttgering, T. W. Shimwell, D. J. B. Smith, C. Tasse, W. L. Williams

Abstract: We present a study of the low-frequency radio properties of broad absorption line quasars (BALQSOs) from the LOFAR Two-metre Sky-Survey Data Release 1 (LDR1). The value-added LDR1 catalogue contains Pan-STARRS counterparts, which we match with the Sloan Digital Sky Survey (SDSS) DR7 and DR12 quasar catalogues. We find that BALQSOs are twice as likely to be detected at 144$\,$MHz than their non-BAL… ▽ More We present a study of the low-frequency radio properties of broad absorption line quasars (BALQSOs) from the LOFAR Two-metre Sky-Survey Data Release 1 (LDR1). The value-added LDR1 catalogue contains Pan-STARRS counterparts, which we match with the Sloan Digital Sky Survey (SDSS) DR7 and DR12 quasar catalogues. We find that BALQSOs are twice as likely to be detected at 144$\,$MHz than their non-BAL counterparts, and BALQSOs with low-ionisation species present in their spectra are three times more likely to be detected than those with only high-ionisation species. The BALQSO fraction at 144$\,$MHz is constant with increasing radio luminosity, which is inconsistent with previous results at 1.4$\,$GHz, indicating that observations at the different frequencies may be tracing different sources of radio emission. We cross-match radio sources between the Faint Images of the Radio Sky at Twenty Centimeters (FIRST) survey and LDR1, which provides a bridge via the LDR1 Pan-STARRS counterparts to identify BALQSOs in SDSS. Consequently we expand the sample of BALQSOs detected in FIRST by a factor of three. The LDR1-detected BALQSOs in our sample are almost exclusively radio-quiet (\logr $\,<2$), with radio sizes at 144$\,$MHz typically less than $200\,$kpc; these radio sizes tend to be larger than those at 1.4$\,$GHz, suggesting more extended radio emission at low frequencies. We find that although the radio detection fraction increases with increasing balnicity index (BI), there is no correlation between BI and either low-frequency radio power or radio-loudness. This suggests that both radio emission and BI may be linked to the same underlying process, but are spatially distinct phenomena. △ Less

Submitted 19 November, 2018; originally announced November 2018.

Comments: 16 figures, 3 tables and 19 pages. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available on lofar-surveys.org upon publication of the journal

Ultra-high energy cosmic rays from shocks in the lobes of powerful radio galaxies

Authors: James H. Matthews, Anthony R. Bell, Katherine M. Blundell, Anabella T. Araudo

Abstract: The origin of ultra-high energy cosmic rays (UHECRs) has been an open question for decades. Here, we use a combination of hydrodynamic simulations and general physical arguments to demonstrate that UHECRs can in principle be produced by diffusive shock acceleration (DSA) in shocks in the backflowing material of radio galaxy lobes. These shocks occur after the jet material has passed through the re… ▽ More The origin of ultra-high energy cosmic rays (UHECRs) has been an open question for decades. Here, we use a combination of hydrodynamic simulations and general physical arguments to demonstrate that UHECRs can in principle be produced by diffusive shock acceleration (DSA) in shocks in the backflowing material of radio galaxy lobes. These shocks occur after the jet material has passed through the relativistic termination shock. Recently, several authors have demonstrated that highly relativistic shocks are not effective in accelerating UHECRs. The shocks in our proposed model have a range of non-relativistic or mildly relativistic shock velocities more conducive to UHECR acceleration, with shock sizes in the range 1-10kpc. Approximately 10% of the jet's energy flux is focused through a shock in the backflow of $M>3$. Although the shock velocities can be low enough that acceleration to high energy via DSA is still efficient, they are also high enough for the Hillas energy to approach $10^{19-20}$eV, particularly for heavier CR composition and in cases where fluid elements pass through multiple shocks. We discuss some of the more general considerations for acceleration of particles to ultra-high energy with reference to giant-lobed radio galaxies such as Centaurus A and Fornax A, a class of sources which may be responsible for the observed anisotropies from UHECR observatories. △ Less

Submitted 29 October, 2018; originally announced October 2018.

Comments: Accepted to MNRAS. 20 pages, 15 figures, 1 table. Movies in supplementary material can be viewed and downloaded at http://jhmatthews.github.io/uhecr-movies

Radiation-hydrodynamic simulations of thermally-driven disc winds in X-ray binaries: A direct comparison to GRO J1655-40

Authors: Nick Higginbottom, Christian Knigge, Knox S. Long, James H. Matthews, Stuart A. Sim, Henrietta A. Hewitt

Abstract: Essentially all low-mass X-ray binaries (LMXBs) in the soft state appear to drive powerful equatorial disc winds. A simple mechanism for driving such outflows involves X-ray heating of the top of the disc atmosphere to the Compton temperature. Beyond the Compton radius, the thermal speed exceeds the escape velocity, and mass loss is inevitable. Here, we present the first coupled radiation-hydrodyn… ▽ More Essentially all low-mass X-ray binaries (LMXBs) in the soft state appear to drive powerful equatorial disc winds. A simple mechanism for driving such outflows involves X-ray heating of the top of the disc atmosphere to the Compton temperature. Beyond the Compton radius, the thermal speed exceeds the escape velocity, and mass loss is inevitable. Here, we present the first coupled radiation-hydrodynamic simulation of such thermally-driven disc winds. The main advance over previous modelling efforts is that the frequency-dependent attenuation of the irradiating SED is taken into account. We can therefore relax the approximation that the wind is optically thin throughout which is unlikely to hold in the crucial acceleration zone of the flow. The main remaining limitations of our simulations are connected to our treatment of optically thick regions. Adopting parameters representative of the wind-driving LMXB GRO~J1655-40, our radiation-hydrodynamic model yields a mass-loss rate that is $\simeq5\times$ lower than that suggested by pure hydrodynamic, optically thin models. This outflow rate still represents more than twice the accretion rate and agrees well with the mass-loss rate inferred from Chandra/HETG observations of GRO~J1655-40 at a time when the system had a similar luminosity to that adopted in our simulations. The Fe XXV and Fe XXVI Lyman $\rmα~$ absorption line profiles observed in this state are slightly stronger than those predicted by our simulations but the qualitative agreement between observed and simulated outflow properties means that thermal driving is a viable mechanism for powering the disc winds seen in soft-state LMXBs. △ Less

Submitted 13 June, 2018; originally announced June 2018.

Fornax A, Centaurus A and other radio galaxies as sources of ultra-high energy cosmic rays

Authors: James H. Matthews, Anthony R. Bell, Katherine M. Blundell, Anabella T. Araudo

Abstract: The origin of ultra-high energy cosmic rays (UHECRs) is still unknown. It has recently been proposed that UHECR anisotropies can be attributed to starburst galaxies or active galactic nuclei. We suggest that the latter is more likely and that giant-lobed radio galaxies such as Centaurus A and Fornax A can explain the data. The origin of ultra-high energy cosmic rays (UHECRs) is still unknown. It has recently been proposed that UHECR anisotropies can be attributed to starburst galaxies or active galactic nuclei. We suggest that the latter is more likely and that giant-lobed radio galaxies such as Centaurus A and Fornax A can explain the data. △ Less

Submitted 31 May, 2018; v1 submitted 4 May, 2018; originally announced May 2018.

Comments: Accepted to MNRAS Letters on 31 May 2018. 5 pages, 2 figures

On the maximum energy of non-thermal particles in the primary hotspot of Cygnus A

Authors: Anabella T. Araudo, Anthony R. Bell, Katherine M. Blundell, James H. Matthews

Abstract: We study particle acceleration and magnetic field amplification in the primary hotspot in the northwest jet of radiogalaxy Cygnus A. By using the observed flux density at 43 GHz in a well resolved region of this hotspot, we determine the minimum value of the jet density and constrain the magnitude of the magnetic field. We find that a jet with density greater than $5\times 10^{-5}$ cm$^{-3}$ and h… ▽ More We study particle acceleration and magnetic field amplification in the primary hotspot in the northwest jet of radiogalaxy Cygnus A. By using the observed flux density at 43 GHz in a well resolved region of this hotspot, we determine the minimum value of the jet density and constrain the magnitude of the magnetic field. We find that a jet with density greater than $5\times 10^{-5}$ cm$^{-3}$ and hotspot magnetic field in the range 50-400 $μ$G are required to explain the synchrotron emission at 43 GHz. The upper-energy cut-off in the hotspot synchrotron spectrum is at a frequency < $5\times 10^{14}$ Hz, indicating that the maximum energy of non-thermal electrons accelerated at the jet reverse shock is $E_{e, \rm max} \sim 0.8$ TeV in a magnetic field of 100 $μ$G. Based on the condition that the magnetic-turbulence scale length has to be larger than the plasma skin depth, and that the energy density in non-thermal particles cannot violate the limit imposed by the jet kinetic luminosity, we show that $E_{e,\rm max}$ cannot be constrained by synchrotron losses as traditionally assumed. In addition to that, and assuming that the shock is quasi-perpendicular, we show that non-resonant hybrid instabilities generated by the streaming of cosmic rays with energy $E_{e, \rm max}$ can grow fast enough to amplify the jet magnetic field up to 50-400 $μ$G and accelerate particles up to the maximum energy $E_{e, \rm max}$ observed in the Cygnus A primary hotspot. △ Less

Submitted 26 September, 2017; originally announced September 2017.

Comments: Accepted for publication in MNRAS

The reverberation signatures of rotating disc winds in active galactic nuclei

Authors: S. W. Mangham, C. Knigge, J. H. Matthews, K. S. Long, S. A. Sim, N. Higginbottom

Abstract: The broad emission lines (BELs) in active galactic nuclei (AGN) respond to ionizing continuum variations. The time and velocity dependence of their response depends on the structure of the broad-line region: its geometry, kinematics and ionization state. Here, we predict the reverberation signatures of BELs formed in rotating accretion disc winds. We use a Monte Carlo radiative transfer and ioniza… ▽ More The broad emission lines (BELs) in active galactic nuclei (AGN) respond to ionizing continuum variations. The time and velocity dependence of their response depends on the structure of the broad-line region: its geometry, kinematics and ionization state. Here, we predict the reverberation signatures of BELs formed in rotating accretion disc winds. We use a Monte Carlo radiative transfer and ionization code to predict velocity-delay maps for representative high- (C$~IV$) and low-ionization (H$α$) emission lines in both high- and moderate-luminosity AGN. Self-shielding, multiple scattering and the ionization structure of the outflows are all self-consistently taken into account, while small-scale structure in the outflow is modelled in the micro-clumping approximation. Our main findings are: (1) The velocity-delay maps of smooth/micro-clumped outflows often contain significant negative responses. (2)~The reverberation signatures of disc wind models tend to be rotation dominated and can even resemble the classic "red-leads-blue" inflow signature. (3) Traditional "blue-leads-red" outflow signatures can usually only be observed in the long-delay limit. (4) Our models predict lag-luminosity relationships similar to those inferred from observations, but systematically underpredict the observed centroid delays. (5) The ratio between "virial product" and black hole mass predicted by our models depends on viewing angle. Our results imply that considerable care needs to be taken in interpreting data obtained by observational reverberation mapping campaigns. In particular, basic signatures such as "red-leads-blue", "blue-leads-red" and "blue and red vary jointly" are not always reliable indicators of inflow, outflow or rotation. This may help to explain the perplexing diversity of such signatures seen in observational campaigns to date. △ Less

Submitted 24 July, 2017; originally announced July 2017.

Comments: 15 pages, 17 figures, 2 tables. Accepted by MNRAS 20/7/2017

Amplification of perpendicular and parallel magnetic fields by cosmic ray currents

Authors: James H. Matthews, Anthony R. Bell, Katherine M. Blundell, Anabella T. Araudo

Abstract: Cosmic ray (CR) currents through magnetised plasma drive strong instabilities producing amplification of the magnetic field. This amplification helps explain the CR energy spectrum as well as observations of supernova remnants and radio galaxy hot spots. Using magnetohydrodynamic (MHD) simulations, we study the behaviour of the non-resonant hybrid (NRH) instability (also known as the Bell instabil… ▽ More Cosmic ray (CR) currents through magnetised plasma drive strong instabilities producing amplification of the magnetic field. This amplification helps explain the CR energy spectrum as well as observations of supernova remnants and radio galaxy hot spots. Using magnetohydrodynamic (MHD) simulations, we study the behaviour of the non-resonant hybrid (NRH) instability (also known as the Bell instability) in the case of CR currents perpendicular and parallel to the initial magnetic field. We demonstrate that extending simulations of the perpendicular case to 3D reveals a different character to the turbulence from that observed in 2D. Despite these differences, in 3D the perpendicular NRH instability still grows exponentially far into the non-linear regime with a similar growth rate to both the 2D perpendicular and 3D parallel situations. We introduce some simple analytical models to elucidate the physical behaviour, using them to demonstrate that the transition to the non-linear regime is governed by the growth of thermal pressure inside dense filaments at the edges of the expanding loops. We discuss our results in the context of supernova remnants and jets in radio galaxies. Our work shows that the NRH instability can amplify magnetic fields to many times their initial value in parallel and perpendicular shocks. △ Less

Submitted 28 July, 2017; v1 submitted 10 April, 2017; originally announced April 2017.

Comments: Published in MNRAS. 14 pages, 12 figures, 2 tables. Replacement corrects some typesetting errors

Journal ref: MNRAS, 469, 1849 (2017)

Quasar emission lines as probes of orientation: implications for disc wind geometries and unification

Authors: James H. Matthews, Christian Knigge, Knox S. Long

Abstract: The incidence of broad absorption lines (BALs) in quasar samples is often interpreted in the context of a geometric unification model consisting of an accretion disc and an associated outflow. We use the the Sloan Digital Sky Survey (SDSS) quasar sample to test this model by examining the equivalent widths (EWs) of CIV 1550Å, Mg II 2800Å, [OIII] 5007Å and C III] 1909Å. We find that the emission li… ▽ More The incidence of broad absorption lines (BALs) in quasar samples is often interpreted in the context of a geometric unification model consisting of an accretion disc and an associated outflow. We use the the Sloan Digital Sky Survey (SDSS) quasar sample to test this model by examining the equivalent widths (EWs) of CIV 1550Å, Mg II 2800Å, [OIII] 5007Å and C III] 1909Å. We find that the emission line EW distributions in BAL and non-BAL quasars are remarkably similar -- a property that is inconsistent with scenarios in which a BAL outflow rises equatorially from a geometrically thin, optically thick accretion disc. We construct simple models to predict the distributions from various geometries; these models confirm the above finding and disfavour equatorial geometries. We show that obscuration, line anisotropy and general relativistic effects on the disc continuum are unlikely to hide an EW inclination dependence. We carefully examine the radio and polarisation properties of BAL quasars. Both suggest that they are most likely viewed (on average) from intermediate inclinations, between type 1 and type 2 AGN. We also find that the low-ionization BAL quasars in our sample are not confined to one region of `Eigenvector I' parameter space. Overall, our work leads to one of the following conclusions, or some combination thereof: (i) the continuum does not emit like a geometrically thin, optically thick disc; (ii) BAL quasars are viewed from similar angles to non-BAL quasars, i.e. low inclinations; (iii) geometric unification does not explain the fraction of BALs in quasar samples. △ Less

Submitted 24 January, 2017; originally announced January 2017.

Comments: Accepted to MNRAS. 16 pages, 11 figures

Testing Quasar Unification: Radiative Transfer in Clumpy Winds

Authors: James H. Matthews, Christian Knigge, Knox S. Long, Stuart A. Sim, Nick Higginbottom, Sam W. Mangham

Abstract: Various unification schemes interpret the complex phenomenology of quasars and luminous active galactic nuclei (AGN) in terms of a simple picture involving a central black hole, an accretion disc and an associated outflow. Here, we continue our tests of this paradigm by comparing quasar spectra to synthetic spectra of biconical disc wind models, produced with our state-of-the-art Monte Carlo radia… ▽ More Various unification schemes interpret the complex phenomenology of quasars and luminous active galactic nuclei (AGN) in terms of a simple picture involving a central black hole, an accretion disc and an associated outflow. Here, we continue our tests of this paradigm by comparing quasar spectra to synthetic spectra of biconical disc wind models, produced with our state-of-the-art Monte Carlo radiative transfer code. Previously, we have shown that we could produce synthetic spectra resembling those of observed broad absorption line (BAL) quasars, but only if the X-ray luminosity was limited to $10^{43}$ erg s$^{-1}$. Here, we introduce a simple treatment of clumping, and find that a filling factor of $\sim0.01$ moderates the ionization state sufficiently for BAL features to form in the rest-frame UV at more realistic X-ray luminosities. Our fiducial model shows good agreement with AGN X-ray properties and the wind produces strong line emission in, e.g., Ly α and CIV 1550Å at low inclinations. At high inclinations, the spectra possess prominent LoBAL features. Despite these successes, we cannot reproduce all emission lines seen in quasar spectra with the correct equivalent-width ratios, and we find an angular dependence of emission-line equivalent width despite the similarities in the observed emission line properties of BAL and non-BAL quasars. Overall, our work suggests that biconical winds can reproduce much of the qualitative behaviour expected from a unified model, but we cannot yet provide quantitative matches with quasar properties at all viewing angles. Whether disc winds can successfully unify quasars is therefore still an open question. △ Less

Submitted 8 February, 2016; originally announced February 2016.

Comments: Accepted for publication in MNRAS. 14 pages, 8 figures. Comments welcome. Simulation grid output spectra available online at http://jhmatthews.github.io/quasar-wind-grid/

The Impact of Accretion Disk Winds on the Optical Spectra of Cataclysmic Variables

Authors: James H. Matthews, Christian Knigge, Knox S. Long, Stuart A. Sim, Nick Higginbottom

Abstract: Many high-state non-magnetic cataclysmic variables (CVs) exhibit blue-shifted absorption or P-Cygni profiles associated with ultraviolet (UV) resonance lines. These features imply the existence of powerful accretion disk winds in CVs. Here, we use our Monte Carlo ionization and radiative transfer code to investigate whether disk wind models that produce realistic UV line profiles are also likely t… ▽ More Many high-state non-magnetic cataclysmic variables (CVs) exhibit blue-shifted absorption or P-Cygni profiles associated with ultraviolet (UV) resonance lines. These features imply the existence of powerful accretion disk winds in CVs. Here, we use our Monte Carlo ionization and radiative transfer code to investigate whether disk wind models that produce realistic UV line profiles are also likely to generate observationally significant recombination line and continuum emission in the optical waveband. We also test whether outflows may be responsible for the single-peaked emission line profiles often seen in high-state CVs and for the weakness of the Balmer absorption edge (relative to simple models of optically thick accretion disks). We find that a standard disk wind model that is successful in reproducing the UV spectra of CVs also leaves a noticeable imprint on the optical spectrum, particularly for systems viewed at high inclination. The strongest optical wind-formed recombination lines are H$α$ and He II $\lambda4686$. We demonstrate that a higher-density outflow model produces all the expected H and He lines and produces a recombination continuum that can fill in the Balmer jump at high inclinations. This model displays reasonable verisimilitude with the optical spectrum of RW Trianguli. No single-peaked emission is seen, although we observe a narrowing of the double-peaked emission lines from the base of the wind. Finally, we show that even denser models can produce a single-peaked H$α$ line. On the basis of our results, we suggest that winds can modify, and perhaps even dominate, the line and continuum emission from CVs. △ Less

Submitted 21 April, 2015; originally announced April 2015.

Comments: 15 pages, 13 figures. Accepted to MNRAS

Line-driven Disk Winds in Active Galactic Nuclei: The Critical Importance of Ionization and Radiative Transfer

Authors: Nick Higginbottom, Daniel Proga, Christian Knigge, Knox S. Long, James H. Matthews, Stuart A. Sim

Abstract: Accretion disk winds are thought to produce many of the characteristic features seen in the spectra of active galactic nuclei (AGN) and quasi-stellar objects (QSOs). These outflows also represent a natural form of feedback between the central supermassive black hole and its host galaxy. The mechanism for driving this mass loss remains unknown, although radiation pressure mediated by spectral lines… ▽ More Accretion disk winds are thought to produce many of the characteristic features seen in the spectra of active galactic nuclei (AGN) and quasi-stellar objects (QSOs). These outflows also represent a natural form of feedback between the central supermassive black hole and its host galaxy. The mechanism for driving this mass loss remains unknown, although radiation pressure mediated by spectral lines is a leading candidate. Here, we calculate the ionization state of, and emergent spectra for, the hydrodynamic simulation of a line-driven disk wind previously presented by Proga & Kallman (2004). To achieve this, we carry out a comprehensive Monte Carlo simulation of the radiative transfer through, and energy exchange within, the predicted outflow. We find that the wind is much more ionized than originally estimated. This is in part because it is much more difficult to shield any wind regions effectively when the outflow itself is allowed to reprocess and redirect ionizing photons. As a result, the calculated spectrum that would be observed from this particular outflow solution would not contain the ultraviolet spectral lines that are observed in many AGN/QSOs. Furthermore, the wind is so highly ionized that line-driving would not actually be efficient. This does not necessarily mean that line-driven winds are not viable. However, our work does illustrate that in order to arrive at a self-consistent model of line-driven disk winds in AGN/QSO, it will be critical to include a more detailed treatment of radiative transfer and ionization in the next generation of hydrodynamic simulations. △ Less

Submitted 13 May, 2014; v1 submitted 8 February, 2014; originally announced February 2014.

Comments: 13 pages, 10 figures - Accepted for publication in ApJ

A Simple Disk Wind Model for Broad Absorption Line Quasars

Authors: N. Higginbottom, C. Knigge, K. S. Long, S. A. Sim, J. H. Matthews

Abstract: Approximately 20% of quasi-stellar objects (QSOs) exhibit broad, blue-shifted absorption lines in their ultraviolet spectra. Such features provide clear evidence for significant outflows from these systems, most likely in the form of accretion disk winds. These winds may represent the "quasar" mode of feedback that is often invoked in galaxy formation/evolution models, and they are also key to uni… ▽ More Approximately 20% of quasi-stellar objects (QSOs) exhibit broad, blue-shifted absorption lines in their ultraviolet spectra. Such features provide clear evidence for significant outflows from these systems, most likely in the form of accretion disk winds. These winds may represent the "quasar" mode of feedback that is often invoked in galaxy formation/evolution models, and they are also key to unification scenarios for active galactic nuclei (AGN) and QSOs. To test these ideas, we construct a simple benchmark model of an equatorial, biconical accretion disk wind in a QSO and use a Monte Carlo ionization/radiative transfer code to calculate the ultraviolet spectra as a function of viewing angle. We find that for plausible outflow parameters, sightlines looking directly into the wind cone do produce broad, blue-shifted absorption features in the transitions typically seen in broad absorption line QSOs. However, our benchmark model is intrinsically X-ray weak in order to prevent overionization of the outflow, and the wind does not yet produce collisionally excited line emission at the level observed in non-BAL QSOs. As a first step towards addressing these shortcomings, we discuss the sensitivity of our results to changes in the assumed X-ray luminosity and mass-loss rate, Mdot(wind). In the context of our adopted geometry, Mdot(wind) ~ Mdot(acc) is required in order to produce significant BAL features. The kinetic luminosity and momentum carried by such outflows would be sufficient to provide significant feedback. △ Less

Submitted 27 August, 2013; originally announced August 2013.

Comments: 20 pages, 13 figures (some in colour), accepted for publication in MNRAS