# Publications by Katherine Blundell

## Cosmic ray acceleration in hydromagnetic flux tubes

Monthly Notices of the Royal Astronomical Society Oxford University Press 487 (2019) 4571–4579-

AR Bell, J Matthews, KM Blundell, AT Araudo

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 magnetized 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.

## Studying the H-alpha line of the B[ e ] supergiant binary GG Carinae using high-cadence optical spectroscopy

Proceedings of the International Astronomical Union (2019) 123-124

A Porter, K Blundell, S Lee

© International Astronomical Union 2019. We present a case study of GG Carinae (GG Car), a Galactic B[ e ] supergiant binary having significant eccentricity (0.28), based on Global Jet Watch spectroscopy data which has been collecting high-time-sampled optical spectra since early 2015. GG Car has so far not been observed in the X-ray band, however it is of similar phenomenology to known X-ray binaries and may therefore be an obscured X-ray source. We have discovered that the absorption component of the H-alpha line displays a '1/462-478-day period in both equivalent width and wavelength centroid indicating cycles in the dynamics of the circumstellar environment, such as precession of the circumbinary or circumprimary disk. Circumbinary disk precession is an as-of-yet underexplored origin of super-orbital variations in the X-ray flux of X-ray binaries, since the rate of precession is generally much longer than the orbital period of the inner binary.

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

Monthly Notices of the Royal Astronomical Society Oxford University Press 482 (2018) 4303–4321-

J Matthews, B Bell, K Blundell, AT Araudo

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 − 10 kpc. Approximately 10% of the jet’s energy flux is focused through a shock in the backflow of M &gt; 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 1019−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.

## Cosmic ray acceleration by shocks: spectral steepening due to turbulent magnetic field amplification

Monthly Notices of the Royal Astronomical Society Oxford University Press 488* (2019) 2466-2472

A Bell, J Matthews, K Blundell

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.

## Cosmic ray acceleration to ultrahigh energy in radio galaxies

EPJ Web of Conferences EDP Sciences 210 (2019) 04002

J Matthews, AR Bell, KM Blundell, AT Araudo

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.

## SS433's Jet Trace from ALMA Imaging and Global Jet Watch Spectroscopy: Evidence for Post-launch Particle Acceleration

ASTROPHYSICAL JOURNAL LETTERS 867 (2018) ARTN L25

KM Blundell, R Laing, S Lee, A Richards

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

Monthly Notices of the Royal Astronomical Society Oxford University Press 473 (2017) 3500–3506-

AT Araudo, AR Bell, KM Blundell, JH Matthews

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 $\mu$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 &lt; $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 $\mu$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 $\mu$G and accelerate particles up to the maximum energy $E_{e, \rm max}$ observed in the Cygnus A primary hotspot.

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

Monthly Notices of the Royal Astronomical Society: Letters Oxford University Press 479 (2018) L76–L80-

JH Matthews, AR Bell, AT Araudo, KM Blundell

The origin of ultra-high energy cosmic rays (UHECRs) is still unknown. It has recently been proposed that UHECR anisotropies can be attributed to starbust 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.

## Cosmic ray acceleration by relativistic shocks: Limits and estimates

Monthly Notices of the Royal Astronomical Society Oxford University Press (2017)

AR Bell, AT Araudo, KM Blundell, JH Matthews

We examine limits to the energy to which cosmic rays can be accelerated by relativistic shocks, showing that acceleration of light ions as high as 100 EeV is unlikely. The implication of our estimates is that if ultra-high energy cosmic rays are accelerated by shocks, then those shocks are probably not relativistic.

## Evidence that particle acceleration in hotspots of FR II galaxies is not constrained by synchrotron cooling

Nuclear and Particle Physics Proceedings 297-299 (2018) 242-248

A Araudo, AR BELL, K BLUNDELL

© 2018 We study the hotspots of powerful radiogalaxies, where electrons accelerated at the jet termination shock emit synchrotron radiation. The turnover of the synchrotron spectrum is typically observed between infrared and optical frequencies, indicating that the maximum energy of non-thermal electrons accelerated at the shock is ≲ TeV for a canonical magnetic field of ∼100 μG. We show that this maximum energy cannot be constrained by synchrotron losses as usually assumed, unless the jet density is unreasonably large and most of the jet upstream energy goes to non-thermal particles. We test this result by considering a sample of hotspots observed at radio, infrared and optical wavelengths.

## LOFAR 150-MHz observations of SS 433 and W 50

MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY 475 (2018) 5360-5377

JW Broderick, RP Fender, JCA Miller-Jones, SA Trushkin, AJ Stewart, GE Anderson, TD Staley, KM Blundell, M Pietka, S Markoff, A Rowlinson, JD Swinbank, AJ van der Horst, ME Bell, RP Breton, D Carbone, S Corbel, J Eisloeffel, H Falcke, J-M Griessmeier, JWT Hessels, VI Kondratiev, CJ Law, GJ Molenaar, M Serylak, BW Stappers, J van Leeuwen, RAMJ Wijers, R Wijnands, MW Wise, P Zarka

## Amplification of perpendicular and parallel magnetic fields by cosmic ray currents

Monthly Notices of the Royal Astronomical Society Oxford University Press 469 (2017) 1849-1860

JH Matthews, A Bell, KM Blundell, AT Araudo

Cosmic ray (CR) currents through magnetized 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 hotspots. Using magnetohydrodynamic 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.

## Deep 230-470 MHz VLA Observations of the mini-halo in the Perseus Cluster

MNRAS 469 (2017) 2017-2017

M Gendron-Marsolais, J Hlavacek-Larrondo, RJV Weeren, T Clarke, AC Fabian, HT Intema, GB Taylor, KM Blundell, JS Sanders

We present a low-frequency view of the Perseus cluster with new observations from the Karl G. Jansky Very Large Array (JVLA) at 230-470 MHz. The data reveal a multitude of new structures associated with the mini-halo. The mini-halo seems to be influenced both by the AGN activity as well as by the sloshing motion of the cool core cluster's gas. In addition, it has a filamentary structure similar to that seen in radio relics found in merging clusters. We present a detailed description of the data reduction and imaging process of the dataset. The depth and resolution of the observations allow us to conduct for the first time a detailed comparison of the mini-halo structure with the X-ray structure as seen in the Chandra X-ray images. The resulting image shows very clearly that the mini-halo emission is mostly contained behind the cold fronts, similar to that predicted by simulations of gas sloshing in galaxy clusters. However, due to the proximity of the Perseus cluster, as well as the quality of the data at low radio frequencies and at X-ray wavelengths, we also find evidence of fine structure. This structure includes several radial radio filaments extending in different directions, a concave radio structure associated with the southern X-ray bay and sharp edges that correlate with X-ray edges. Mini-halos are therefore not simply diffuse, uniform radio sources, but are rather filled with a rich variety of complex structures. These results illustrate the high-quality images that can be obtained with the new JVLA at low radio-frequencies, as well as the necessity to obtain deeper, higher-fidelity radio images of mini-halos and halos in clusters to further understand their origin.

## Fast launch speeds in radio flares, from a new determination of the intrinsic motions of SS 433's jet bolides

Monthly Notices of the Royal Astronomical Society Oxford University Press 461 (2016) 312–320-

K Blundell, RM Jeffrey, SA Trushkin, AJ Mioduszewski

We present new high-resolution, multi-epoch, VLBA radio images of the Galactic microquasar SS 433. We are able to observe plasma knots in the milliarcsecond-scale jets more than 50 days after their launch. This unprecedented baseline in time allows us to determine the bulk launch speed of the radio-emitting plasma during a radio flare, using a new method which we present here, and which is completely independent of optical spectroscopy. We also apply this method to an earlier sequence of 39 short daily VLBA observations, which cover a period in which SS 433 moved from quiescence into a flare. In both datasets we find, for the first time at radio wavebands, clear evidence that the launch speeds of the milliarcsecondscale jets rise as high as 0.32c during flaring episodes. By comparing these images of SS 433 with photometric radio monitoring from the RATAN telescope, we explore further properties of these radio flares.

## A very deep Chandra view of metals, sloshing and feedback in the Centaurus cluster of galaxies

Monthly Notices of the Royal Astronomical Society Oxford University Press 457 (2016) 82–109-

K Blundell

We examine deep Chandra X-ray observations of the Centaurus cluster of galaxies, Abell 3526. Applying a gradient magnitude filter reveals a wealth of structure, from filamentary soft emission on 100 pc (0.5 arcsec) scales close to the nucleus to features 10 s of kpc in size at larger radii. The cluster contains multiple high-metallicity regions with sharp edges. Relative to an azimuthal average, the deviations of metallicity and surface brightness are correlated, and the temperature is inversely correlated, as expected if the larger scale asymmetries in the cluster are dominated by sloshing motions. Around the western cold front are a series of ∼7 kpc ‘notches’, suggestive of Kelvin–Helmholtz instabilities. The cold front width varies from 4 kpc down to close to the electron mean free path. Inside the front are multiple metallicity blobs on scales of 5–10 kpc, which could have been uplifted by AGN activity, also explaining the central metallicity drop and flat inner metallicity profile. Close to the nucleus are multiple shocks, including a 1.9-kpc-radius inner shell-like structure and a weak 1.1–1.4 Mach number shock around the central cavities. Within a 10 kpc radius are nine depressions in surface brightness, several of which appear to be associated with radio emission. The shocks and cavities imply that the nucleus has been repeatedly active on 5–10 Myr time-scales, indicating a tight balance between heating and cooling. We confirm the presence of a series of linear quasi-periodic structures. If they are sound waves, the ∼5 kpc spacing implies a period of 6 Myr, similar to the ages of the shocks and cavities. Alternatively, these structures may be Kelvin–Helmholtz instabilities, their associated turbulence or amplified magnetic field layers.

## Particle acceleration and magnetic field amplification in hotspots of FR II galaxies: The case study 4C74.26

Argentinian Astronomical Society (2016)

AR Bell, AT Araudo, KM Blundell

It has been suggested that relativistic shocks in extragalactic sources may accelerate the most energetic cosmic rays. However, recent theoretical advances indicating that relativistic shocks are probably unable to accelerate particles to energies much larger than a PeV cast doubt on this. In the present contribution we model the radio to X-ray emission in the southern hotspot of the quasar 4C74.26. The synchrotron radio emission is resolved near the shock with the MERLIN radio-interferometer, and the rapid decay of this emission behind the shock is interpreted as the decay of the downstream magnetic field as expected for small scale turbulence. If our result is confirmed by analyses of other radiogalaxies, it provides firm observational evidence that relativistic shocks at the termination region of powerful jets in FR II radiogalaxies do not accelerate ultra high energy cosmic rays.

## Evidence that the maximum electron energy in hotspots of FR II galaxies is not determined by synchrotron cooling

Monthly Notices of the Royal Astronomical Society (2016) stw1204-stw1204

AT Araudo, AR Bell, A Crilly, KM Blundell

## A very deep Chandra view of metals, sloshing and feedback in the Centaurus cluster of galaxies

Monthly Notices of the Royal Astronomical Society Oxford University Press 457 (2016) 82-109

REA Canning, J Hlavacek-Larrondo, AC Fabian, JS Sanders, GB Taylor, K Blundell, HR Russell, SA Walker, C Grimes

We examine deep Chandra X-ray observations of the Centaurus cluster of galaxies, Abell 3526. Applying a gradient magnitude filter reveals a wealth of structure, from filamentary soft emission on 100pc (0.5 arcsec) scales close to the nucleus to features 10s of kpc in size at larger radii. The cluster contains multiple high-metallicity regions with sharp edges. Relative to an azimuthal average, the deviations of metallicity and surface brightness are correlated, and the temperature is inversely correlated, as expected if the larger scale asymmetries in the cluster are dominated by sloshing motions. Around the western cold front are a series of ~7 kpc 'notches', suggestive of Kelvin-Helmholtz instabilities. The cold front width varies from 4 kpc down to close to the electron mean free path. Inside the front are multiple metallicity blobs on scales of 5-10 kpc, which could have been uplifted by AGN activity, also explaining the central metallicity drop and flat inner metallicity profile. Close to the nucleus are multiple shocks, including a 1.9-kpc-radius inner shell-like structure and a weak 1.1-1.4 Mach number shock around the central cavities. Within a 10 kpc radius are 9 depressions in surface brightness, several of which appear to be associated with radio emission. The shocks and cavities imply that the nucleus has been repeatedly active on 5-10 Myr timescales, indicating a tight balance between heating and cooling. We confirm the presence of a series of linear quasi-periodic structures. If they are sound waves, the ~5 kpc spacing implies a period of 6 Myr, similar to the ages of the shocks and cavities. Alternatively, these structures may be Kelvin-Helmholtz instabilities, their associated turbulence or amplified magnetic field layers.

## Black Holes: A Very Short Introduction

Oxford University Press, 2015

K Blundell

In this Very Short Introduction, Katherine Blundell addresses a variety of questions, including what a black hole actually is, how they are characterised and discovered, to what happens if you get too close to one.

## Particle acceleration and magnetic field amplification in the jets of 4C74.26

Astrophysical Journal American Astronomical Society 806 (2015) ARTN 243-

A Bell, AT Araudo, K Blundell

We model the multi-wavelength emission in the southern hotspot of the radio quasar 4C74.26. The synchrotron radio emission is resolved near the shock with the MERLIN radio-interferometer, and the rapid decay of this emission behind the shock is interpreted as the decay of the amplified downstream magnetic field as expected for small scale turbulence. Electrons are accelerated to only 0.3 TeV, consistent with a diffusion coefficient many orders of magnitude greater than in the Bohm regime. If the same diffusion coefficient applies to the protons, their maximum energy is only ~100 TeV.