Search

The neutron star low-mass X-ray binary, EXO 0748--676, recently returned to outburst after a $\sim$ 16 year-long quiescence. Since its return, there has been a global effort to capture the previously unseen rise of the source and to understand its somewhat early return to outburst, as it is typical for a source to spend longer in quiescence than in outburst. Here, we report on the simultaneous optical and X-ray detection of a type I X-ray burst, captured by XMM-Newton during a DDT observation on 30th June 2024. The data show 3 X-ray eclipses consistent with the known ephemeris and one type I X-ray burst at 60492.309 MJD. The X-ray burst is reprocessed into the optical band and captured by XMM-Newton's Optical Monitor during a 4399 s exposure with the B filter in image + fast mode. We determine that the optical peak lags the X-ray peak by 4.46 $\pm$ 1.71s. The optical and X-ray rise times are similar, but the optical decay timescale is shorter than the X-ray decay timescale. The reprocessing site is likely within a few light seconds of the X-ray emitting region, so the companion star, accretion disc and ablated material are all plausible., Comment: 6 Pages, 3 Figures, Accepted for Publication in MNRAS Letters

Outflows are critical components of many astrophysical systems, including accreting compact binaries and active galactic nuclei (AGN). These outflows can significantly affect a system's evolution and alter its observational appearance by reprocessing the radiation produced by the central engine. Sirocco (Simulating Ionization and Radiation in Outflows Created by Compact Objects - or "the code formerly known as Python") is a Sobolev-based Monte Carlo ionization and radiative transfer code. It is designed to simulate the spectra produced by any system with an azimuthally-symmetric outflow, from spherical stellar winds to rotating, biconical accretion disc winds. Wind models can either be parametrized or imported, e.g. from hydrodynamical simulations. The radiation sources include an optically thick accretion disc and various central sources with flexible spectra and geometries. The code tracks the "photon packets" produced by the sources in any given simulation as they traverse and interact with the wind. The code assumes radiative near-equilibrium, so the thermal and ionization state can be determined iteratively from these interactions. Once the physical properties in the wind have converged, Sirocco can be used to generate synthetic spectra at a series of observer sightlines. Here, we describe the physical assumptions, operation, performance and limitations of the code. We validate it against tardis, cmfgen and cloudy, finding good agreement, and present illustrative synthetic spectra from disc winds in cataclysmic variables, tidal disruption events, AGN and X-ray binaries. Sirocco is publicly available on GitHub, alongside its associated data, documentation and sample input files covering a wide range of astrophysical applications., Comment: 26 pages, 16 figures. Submitted to MNRAS. This is the release paper for the SIROCCO code, which can be found at https://github.com/sirocco-rt/sirocco with links to documentation. Underlying data and figure scripts available at https://github.com/sirocco-rt/release-models -- comments on the code, paper or documentation are welcomed

Despite being operational for only a short time, the Einstein Probe mission has already significantly advanced the study of rapid variability in the soft X-ray sky. We report the discovery of luminous and variable radio emission from the Einstein Probe fast X-ray transient EP240414a, the second such source with a radio counterpart. The radio emission at $3\,\rm{GHz}$ peaks at $\sim30$ days post explosion and with a spectral luminosity $\sim2\times10^{30}\,\rm{erg}\,\rm{s}^{-1}\,\rm{Hz}^{-1}$, similar to what is seen from long gamma-ray bursts, and distinct from other extra-galactic transients including supernovae and tidal disruption events, although we cannot completely rule out emission from engine driven stellar explosions e.g. the fast blue optical transients. An equipartition analysis of our radio data reveals that an outflow with at least a moderate bulk Lorentz factor ($\Gamma\gtrsim1.6$) with a minimum energy of $\sim10^{48}\,\rm{erg}$ is required to explain our observations. The apparent lack of reported gamma-ray counterpart to EP240414a could suggest that an off-axis or choked jet could be responsible for the radio emission, although a low luminosity gamma-ray burst may have gone undetected. Our observations are consistent with the hypothesis that a significant fraction of extragalactic fast X-ray transients are associated with the deaths of massive stars., Comment: 11 pages, 6 figures, submitted to ApJL

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

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., Comment: 13 pages, 9 figures, 3 tables. Accepted for publication in MNRAS

Understanding the links between different phases of outflows from active galactic nuclei is a key goal in extragalactic astrophysics. Here we compare [OIII] $\lambda\lambda$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 $\lambda$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., Comment: 14 pages, 6 figures. Accepted for publication in MNRAS

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., Comment: 13+2 pages, 6 figures, software code available at https://github.com/sebhoof/bhsr

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., Comment: Accepted for publication in MNRAS. 14 pages, 10 figures + 3 figures in Appendix

Blueshifts - or, more accurately, blue asymmetries - in broad emission lines such as CIV $\lambda$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 $\lambda$1550 emission blueshift versus equivalent-width space and explore the implications for quasar disc wind physics., Comment: 20 pages, 12 figures + 1 appendix. Supplementary material, simulation data and plotting scripts available at: https://github.com/jhmatthews/blueshifts

We analyse the rest-frame UV spectra of 2,531 high-redshift (3.53.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 0.2 and Mbh > 10^9 Mo for the entire 1.56 are on average less massive and have higher Eddington-scaled accretion rates than their luminosity-matched counterparts at 1.5

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'., Comment: 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

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\gamma} > (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., Comment: 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"

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., Comment: Resubmitted to MNRAS. Minor edits to text; figures updated; licence updated

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 $\beta 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., Comment: 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

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\gamma} > 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\gamma}$ 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\gamma} > 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., Comment: 22 pages, 10 figures, 3 tables. Accepted for publication at the Astrophysical Journal. Comments welcome

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., Comment: 13 pages, 8 figures, 3 tables. Accepted for publication in MNRAS. Amended to include two missing references

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 \gamma}$, 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 \gamma}$ (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 $\beta_{\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 $\beta_{\rm pl}$., Comment: 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

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., Comment: 10 Pages, 20 Figures

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, Comment: Accepted to MNRAS. 19 pages, 11 figures

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 $\beta = 100$, we exclude all couplings $g_\mathrm{a\gamma} > 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., Comment: 14 pages, 7 figures. Accepted for publication in MNRAS on 24th November 2021. Comments welcome

We explore the possibility that the hotspot of ultrahigh energy cosmic rays (UHECR) detected by the Telescope Array (TA) from the approximate direction of M82 and the M81 group of galaxies might be the echo of UHECR emitted by Centaurus A in an earlier more powerful phase. Echoes from other starburst galaxies or groups of galaxies may contribute to the UHECR flux at the Earth. We use an illustrative Monte Carlo model of mono-energetic UHECR transport by small angle scattering to generate synthetic sky maps. The model informs a discussion of overall energetics and time and distance scales. We find a viable echo model for the observed UHECR hotspots if the UHECR luminosity of Centaurus A 20 Myr ago was 200 times its present luminosity and if the ordered magnetic field exceeds 10-20nG out to a distance of 400-800kpc in the circumgalactic medium of M82 and other starburst galaxies., Comment: 9 pages, 8 figures. Accepted 2022 January 4th for publication in MNRAS. Replaces the original submitted version

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_{\gamma \to a}$, is given by the one-dimensional power spectrum of the magnetic field components perpendicular to the particle trajectory. Equivalently, we express $P_{\gamma \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>\omega_{\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 $\omega_{\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., Comment: 39 pages, 14 figures. v3: matches published version

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\mu$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., Comment: 21 pages, 7 figures, submitted to AAS journals

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., Comment: 17 pages, 11 figures. Accepted to MNRAS

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$\lambda$1549-emission line that is used to infer the strength of accretion disc winds, and the HeII$\lambda$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., Comment: 13 pages, 10 figures, plus appendices. Minor changes to match published version

We solve the linearised Vlasov-Fokker-Planck (VFP) equation to show that heat flow or an electrical current in a magnetized collisional plasma is unstable to the growth of a circularly polarised transverse perturbation to a zeroth order uniform magnetic field. The Braginskii (1965) transport equations exhibit the same instability in the appropriate limit. This is relevant to laser-produced plasmas, inertial fusion energy (IFE) and to dense cold interstellar plasmas., Comment: Accepted for publication in Plasma Physics and Controlled Fusion

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~\AA~ and an enhancement of N \textsc{v}~1240~\AA~ 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., Comment: 16 pages, 9 figures

In this chapter, we review some features of particle acceleration in astrophysical jets. We begin by describing four observational results relating to the topic, with particular emphasis on jets in active galactic nuclei and parallels between different sources. We then discuss the ways in which particles can be accelerated to high energies in magnetised plasmas, focusing mainly on shock acceleration, second-order Fermi and magnetic reconnection; in the process, we attempt to shed some light on the basic conditions that must be met by any mechanism for the various observational constraints to be satisfied. We describe the limiting factors for the maximum particle energy and briefly discuss multimessenger signals from neutrinos and ultrahigh energy cosmic rays, before describing the journey of jet plasma from jet launch to cocoon with reference to the different acceleration mechanisms. We conclude with some general comments on the future outlook., Comment: To appear in New Astronomy Reviews special volume '100 Years of Astrophysical Jets' (eds. Rob Fender and Ralph Wijers). 23 pages, 8 figures

Determining accurate orbits of binary stars with powerful winds is challenging. The dense outflows increase the effective photospheric radius, precluding direct observation of the Keplerian motion; instead the observables are broad lines emitted over large radii in the stellar wind. Our analysis reveals strong, systematic discrepancies between the radial velocities extracted from different spectral lines: the more extended a line's emission region, the greater the departure from the true orbital motion. To overcome these challenges, we formulate a novel semi-analytical model which encapsulates both the star's orbital motion and the propagation of the wind. The model encodes the integrated velocity field of the out-flowing gas in terms of a convolution of past motion due to the finite flow speed of the wind. We test this model on two binary systems. (1), for the extreme case $\eta$ Carinae, in which the effects are most prominent, we are able to fit the model to 10 Balmer lines from H-alpha to H-kappa concurrently with a single set of orbital parameters: time of periastron $T_{0}=2454848$ (JD), eccentricity $e=0.91$, semi-amplitude $k=69$ km/s and longitude of periastron $\omega=241^\circ$. (2) for a more typical case, the Wolf-Rayet star in RMC 140, we demonstrate that for commonly used lines, such as He II and N III/IV/V, we expect deviations between the Keplerian orbit and the predicted radial velocities. Our study indicates that corrective modelling, such as presented here, is necessary in order to identify a consistent set of orbital parameters, independent of the emission line used, especially for future high accuracy work., Comment: 20 pages, 15 figures, accepted for publication in MNRAS

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., Comment: Accepted for publication by MNRAS

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 ($\phi_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., Comment: 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

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., Comment: MNRAS Letters in press. Updated to match published version

We show that the energy required to turbulently amplify magnetic field during cosmic ray (CR) acceleration by shocks extracts energy from the CR and steepens the CR energy spectrum., Comment: Accepted for publication in the Monthly Notices of the Royal Astronomical Society

We find that hydromagnetic flux tubes in back-flows in the lobes of radio galaxies offer a suitable environment for the acceleration of cosmic rays (CR) to ultra-high energies. We show that CR can reach the Hillas (1984) energy even if the magnetised turbulence in the flux tube is not sufficiently strong for Bohm diffusion to apply. First-order Fermi acceleration by successive weak shocks in a hydromagnetic flux tube is shown to be equivalent to second-order Fermi acceleration by strong turbulence., Comment: Accepted for publication in the Monthly Notices of the Royal Astronomical Society on 31st May 2019

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., Comment: 6 pages, 4 figures. Proceedings of UHECR 2018, 8-12 October 2018, Paris, France

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., Comment: Accepted for publication in MNRAS