Arcminute Microkelvin Imager Research Articles

Probing the evolution of galaxy clusters using SZ effect and non-thermal emission: First results from A1413

Mass is the most fundamental property of galaxy clusters. However, measuring it is still a challenge. Calibrating mass from intracluster medium observables such as the Sunyaev-Zel’dovich (SZ) effect is subject to uncertainty and biases because of the hydrostatic equilibrium assumption. On the other hand, merging cluster systems have been shown to exhibit radio emission which implies a link with disturbances from hydrostatic equilibrium. We present work on studying deviations of galaxy cluster gas pressure profile from the average (universal) pressure profile using an example of galaxy cluster Abell 1413 with SZ effect data from the Arcminute Microkelvin Imager and Planck. This cluster has also been observed at low radio frequency with the Murchison Widefield Array allowing the investigation of links between gas pressure profile deviations and the presence of radio emission.

Day-time-scale variability in the radio light curve of the Tidal Disruption Event AT2022cmc: confirmation of a highly relativistic outflow

ABSTRACTTidal disruption events (TDEs) are transient, multiwavelength events in which a star is ripped apart by a supermassive black hole. Observations show that in a small fraction of TDEs, a short-lived, synchrotron emitting jet is produced. We observed the newly discovered TDE AT2022cmc with a slew of radio facilities over the first 100 d after its discovery. The light curve from the Arcminute Microkelvin Imager Large Array radio interferometer shows day-time-scale variability which we attribute to a high brightness temperature emitting region as opposed to scintillation. We measure a brightness temperature of 2 × 1015 K, which is unphysical for synchrotron radiation. We suggest that the measured high brightness temperature is a result of relativistic beaming caused by a jet being launched at velocities close to the speed of light along our line of sight. We infer from day-time-scale variability that the jet associated with AT2022cmc has a relativistic Doppler factor of at least 16, which corresponds to a bulk Lorentz factor of at least 8, if we are observing the jet directly on axis. Such an inference is the first conclusive evidence that the radio emission observed from some TDEs is from relativistic jets because it does not rely on an outflow model. We also compare the first 100 d of radio evolution of AT2022cmc with that of the previous bright relativistic TDE,Swift J1644, and find a remarkable similarity in their evolution.

Radio observations of the Black Hole X-ray Binary EXO 1846−031 re-awakening from a 34-year slumber

ABSTRACT We present radio [1.3 GHz MeerKAT, 4–8 GHz Karl G. Jansky Very Large Array (VLA), and 15.5 GHz Arcminute Microkelvin Imager Large Array (AMI-LA)] and X-ray (Swift and MAXI) data from the 2019 outburst of the candidate Black Hole X-ray Binary (BHXB) EXO 1846−031. We compute a Hardness–Intensity diagram, which shows the characteristic q-shaped hysteresis of BHXBs in outburst. EXO 1846−031 was monitored weekly with MeerKAT and approximately daily with AMI-LA. The VLA observations provide sub-arcsecond-resolution images at key points in the outburst, showing moving radio components. The radio and X-ray light curves broadly follow each other, showing a peak on ∼MJD 58702, followed by a short decline before a second peak between ∼MJD 58731–58739. We estimate the minimum energy of these radio flares from equipartition, calculating values of Emin ∼ 4 × 1041 and 5 × 1042 erg, respectively. The exact date of the return to ‘quiescence’ is missed in the X-ray and radio observations, but we suggest that it likely occurred between MJD 58887 and 58905. From the Swift X-ray flux on MJD 58905 and assuming the soft-to-hard transition happened at 0.3–3 per cent Eddington, we calculate a distance range of 2.4–7.5 kpc. We computed the radio:X-ray plane for EXO 1846−031 in the ‘hard’ state, showing that it is most likely a ‘radio-quiet’ BH, preferentially at 4.5 kpc. Using this distance and a jet inclination angle of θ = 73°, the VLA data place limits on the intrinsic jet speed of βint = 0.29c, indicating subluminal jet motion.

Long-term radio monitoring of the neutron star X-ray binary Swift J1858.6−0814

ABSTRACT We present the results of our long-term radio monitoring campaign at 1.3 GHz (MeerKAT) and 15.5 GHz (Arcminute Microkelvin Imager – Large Array, AMI-LA) for the outburst of the recently discovered neutron star X-ray binary Swift J1858.6−0814. Throughout the outburst, we observe radio emission consistent with a quasi-persistent, self-absorbed jet. In addition, we see two flares at MJD 58427 and 58530. The second flare allows us to place constraints on the magnetic field and minimum energy of the jet at 0.2 G and 5 × 1037 erg, respectively. We use the multifrequency radio data in conjunction with data from Swift-BAT (Burst Alert Telescope) to place Swift J1858.6−0814 on the radio/X-ray correlation. We find that the quasi-simultaneous radio and BAT data make Swift J1858.6−0814 appear to bridge the gap in the radio/X-ray plane between atoll and Z sources. Furthermore, AMI-LA observations made while Swift J1858.6−0814 was in the soft state have allowed us to show that the radio emission during the soft state is quenched by at least a factor of 4.

A 15.5 GHz detection of the galaxy cluster minihalo in RXJ1720.1+2638

ABSTRACT RXJ1720.1+2638 is a cool-core, ‘relaxed-appearing’ cluster with a minihalo previously detected up to 8.4 GHz, confined by X-ray-detected cold fronts. We present observations of the minihalo at 13–18 GHz with the Arcminute Microkelvin Imager telescope, simultaneously modelling the Sunyaev–Zel’dovich signal of the cluster in conjunction with Planck and Chandra data in order to disentangle the non-thermal emission of the minihalo. We show that the previously reported steepening of the minihalo emission at 8.4 GHz is not supported by the AMI data and that the spectrum is consistent with a single power law up to 18 GHz. We also show the presence of a larger scale component of the minihalo extending beyond the cold fronts. Both of these observations could be explained by the ‘hadronic’ or ‘secondary’ mechanism for the production of relativistic electrons, rather than the currently favoured ‘re-acceleration’ mechanism and/or multiple episodes of jet activity from the active galactic nucleus in the brightest cluster galaxy.

Strong low-frequency radio flaring from Cygnus X-3 observed with LOFAR

ABSTRACT We present Low-Frequency Array (LOFAR) 143.5-MHz radio observations of flaring activity during 2019 May from the X-ray binary Cygnus X-3. Similar to radio observations of previous outbursts from Cygnus X-3, we find that this source was significantly variable at low frequencies, reaching a maximum flux density of about 5.8 Jy. We compare our LOFAR light curve with contemporaneous observations taken at 1.25 and 2.3 GHz with the RATAN-600 telescope, and at 15 GHz with the Arcminute Microkelvin Imager (AMI) Large Array. The initial 143.5-MHz flux density level, ∼2 Jy, is suggested to be the delayed and possibly blended emission from at least some of the flaring activity that had been detected at higher frequencies before our LOFAR observations had begun. There is also evidence of a delay of more than 4 d between a bright flare that initially peaked on May 6 at 2.3 and 15 GHz, and the corresponding peak (≳ 5.8 Jy) at 143.5 MHz. From the multifrequency light curves, we estimate the minimum energy and magnetic field required to produce this flare to be roughly 1044 erg and 40 mG, respectively, corresponding to a minimum mean power of ∼1038 erg s−1. Additionally, we show that the broadband radio spectrum evolved over the course of our observing campaign; in particular, the two-point spectral index between 143.5 MHz and 1.25 GHz transitioned from being optically thick to optically thin as the flare simultaneously brightened at 143.5 MHz and faded at GHz frequencies.

Radio/X-ray correlations and variability in the X-ray binary LS I +61°303

ABSTRACT The high-mass X-ray binary LS I +61°303 exhibits variability in its radio and X-ray emissions, ranging from minute to hour time-scales. At such short time-scales, not much is known about the possible correlations between these two emissions from this source, which might offer hints to their origin. Here, we study the relationship between these emissions using simultaneous X-ray and radio monitoring. We present new radio observations using the Arcminute Microkelvin Imager Large Array telescope at two frequency bands, 13–15.5 and 15.5–18 GHz. We also describe new X-ray observations performed using the XMM–Newton telescope. These X-ray and radio observations overlapped for five hours. We find for the first time that the radio and X-ray emission are correlated up to 81${{\ \rm per\ cent}}$ with their few per cent variability correlated up to 40 ${{\ \rm per\ cent}}$. We discuss possible physical scenarios that produces the observed correlations and variability in the radio and X-ray emission of LS I +61°303.

Radio afterglows of very high-energy gamma-ray bursts 190829A and 180720B

ABSTRACT We present high-cadence multifrequency radio observations of the long gamma-ray burst (GRB) 190829A, which was detected at photon energies above 100 GeV by the High Energy Stereoscopic System (H.E.S.S.). Observations with the Meer Karoo Array Telescope (MeerKAT, 1.3 GHz) and Arcminute Microkelvin Imager – Large Array (AMI-LA, 15.5 GHz) began one day post-burst and lasted nearly 200 d. We used complementary data from Swift X-Ray Telescope (XRT), which ran to 100 d post-burst. We detected a likely forward shock component with both MeerKAT and XRT up to over 100 d post-burst. Conversely, the AMI-LA light curve appears to be dominated by reverse shock emission until around 70 d post-burst when the afterglow flux drops below the level of the host galaxy. We also present previously unpublished observations of the other H.E.S.S.-detected GRB, GRB 180720B from AMI-LA, which shows likely forward shock emission that fades in less than 10 d. We present a comparison between the radio emission from the three GRBs with detected very high energy (VHE) gamma-ray emission and a sensitivity-limited radio afterglow sample. GRB 190829A has the lowest isotropic radio luminosity of any GRB in our sample, but the distribution of luminosities is otherwise consistent, as expected, with the VHE GRBs being drawn from the same parent distribution as the other radio-detected long GRBs.

A Rapid Change in X-Ray Variability and a Jet Ejection in the Black Hole Transient MAXI J1820+070

Abstract We present Neutron Star Interior Composition Explorer X-ray and Arcminute Microkelvin Imager Large Array radio observations of a rapid hard-to-soft state transition in the black hole X-ray transient MAXI J1820+070. During the transition from the hard state to the soft state a switch between two particular types of quasiperiodic oscillations (QPOs) was seen in the X-ray power density spectra, from type-C to type-B, along with a drop in the strength of the broadband X-ray variability and a brief flare in the 7–12 keV band. Soon after this switch (∼2–2.5 hr) a strong radio flare was observed that corresponded to the launch of superluminal ejecta. Although hints of a connection between QPO transitions and radio flares have been seen in other black hole X-ray transients, our observations constitute the strongest observational evidence to date for a link between the appearance of type-B QPOs and the launch of discrete jet ejections.

Arcminute MicroKelvin Imager observations at 15 GHz of the 2020 February Outburst of Cygnus X-3

Here we report observations of a recent giant radio flare from Cygnus X-3 in 2020 February, made with the Arcminute MicroKelvin Imager.

Late-outburst radio flaring in SS Cyg and evidence for a powerful kinetic output channel in cataclysmic variables

Abstract Accreting white dwarfs in binary systems known as cataclysmic variables (CVs) have in recent years been shown to produce radio flares during outbursts, qualitatively similar to those observed from neutron star and black hole X-ray binaries, but their ubiquity and energetic significance for the accretion flow has remained uncertain. We present new radio observations of the CV SS Cyg with Arcminute Microkelvin Imager Large Array, which show for the second time late-ouburst radio flaring, in 2016 April. This flaring occurs during the optical flux decay phase, about 10 d after the well-established early-time radio flaring. We infer that both the early- and late-outburst flares are a common feature of the radio outbursts of SS Cyg, albeit of variable amplitudes, and probably of all dwarf novae. We furthermore present new analysis of the physical conditions in the best-sampled late-outburst flare, from 2016 February, which showed clear optical depth evolution. From this we can infer that the synchrotron-emitting plasma was expanding at about 1 per cent of the speed of light, and at peak had a magnetic field of order 1 G and total energy content ≥1033 erg. While this result is independent of the geometry of the synchrotron-emitting region, the most likely origin is in a jet carrying away a significant amount of the available accretion power.

Physical modelling of galaxy cluster Sunyaev–Zel’dovich data using Einasto dark matter profiles

ABSTRACT We derive a model for Sunyaev–Zel’dovich data from a galaxy cluster that uses an Einasto profile to model the cluster’s dark matter component. This model is similar to the physical models for clusters previously used by the Arcminute Microkelvin Imager (AMI) consortium, which model the dark matter using a Navarro–Frenk–White (NFW) profile, but the Einasto profile provides an extra degree of freedom. We thus present a comparison between two physical models which differ only in the way they model dark matter: one which uses an NFW profile (PM I) and one that uses an Einasto profile (PM II). We illustrate the differences between the models by plotting physical properties of clusters as a function of cluster radius. We generate AMI simulations of clusters that are created and analysed with both models. From this we find that for 14 of the 16 simulations, the Bayesian evidence gives no preference to either of the models according to the Jeffreys scale, and for the other two simulations, weak preference in favour of the correct model. However, for the mass estimates obtained from the analyses, the values were within 1σ of the input values for 14 out of 16 of the clusters when using the correct model, but only in 6 out of 16 cases when the incorrect model was used to analyse the data. Finally, we apply the models to real data from cluster A611 obtained with AMI, and find the mass estimates to be consistent with one another except in the case of when PM II is applied using an extreme value for the Einasto shape parameter.

A detailed radio study of the energetic, nearby, and puzzling GRB 171010A

We present the results of an intensive multi-epoch radio frequency campaign on the energetic and nearby GRB 171010A with the Karl G. Janksy Very Large Array and Arcminute Microkelvin Imager Large Array. We began observing GRB 171010A a day after its initial detection, and were able to monitor the temporal and spectral evolution of the source over the following weeks. The spectra and their evolution are compared to the canonical theories for broadband GRB afterglows, with which we find a general agreement. There are, however, a number of features that are challenging to explain with a simple forward shock model, and we discuss possible reasons for these discrepancies. This includes the consideration of the existence of a reverse shock component, potential microphysical parameter evolution and the effect of scintillation.

Sunyaev–Zel’dovich profile fitting with joint AMI-Planck analysis

ABSTRACT We develop a Bayesian method of analysing Sunyaev–Zel’dovich measurements of galaxy clusters obtained from the Arcminute Microkelvin Imager (AMI) radio interferometer system and from the Planck satellite, using a joint likelihood function for the data from both instruments. Our method is applicable to any combination of Planck data with interferometric data from one or more arrays. We apply the analysis to simulated clusters and find that when the cluster pressure profile is known a priori, the joint data set provides precise and accurate constraints on the cluster parameters, removing the need for external information to reduce the parameter degeneracy. When the pressure profile deviates from that assumed for the fit, the constraints become biased. Allowing the pressure profile shape parameters to vary in the analysis allows an unbiased recovery of the integrated cluster signal and produces constraints on some shape parameters, depending on the angular size of the cluster. When applied to real data from Planck-detected cluster PSZ2 G063.80+11.42, our method resolves the discrepancy between the AMI and Planck Y-estimates and usefully constrains the gas pressure profile shape parameters at intermediate and large radii.

Physical modelling of galaxy clusters detected by thePlancksatellite

We present a comparison of mass estimates for $54$ galaxy cluster candidates from the second Planck catalogue (PSZ2) of Sunyaev-Zel'dovich sources. We compare the mass values obtained with data taken from the Arcminute Microkelvin Imager (AMI) radio interferometer system and from the Planck satellite. The former of these uses a Bayesian analysis pipeline that parameterises a cluster in terms of its physical quantities, and models the dark matter & baryonic components of a cluster using NFW and GNFW profiles respectively. Our mass estimates derived from Planck data are obtained from the results of the Bayesian detection algorithm PowellSnakes (PwS), are based on the methodology detailed in the PSZ2 paper, and produce two sets of mass estimates; one estimate is calculated directly from the angular radius $\theta$ - integrated Comptonisation parameter $Y$ posterior distributions, and the other uses a `slicing function' to provide information on $\theta$ based on X-ray measurements and previous Planck mission samples. We find that for $37$ of the clusters, the AMI mass estimates are lower than both values obtained from Planck data. However the AMI and slicing function estimates are within one combined standard deviation of each other for $31$ clusters. We also generate cluster simulations based on the slicing-function mass estimates, and analyse them in the same way as we did the real AMI data. We find that inclusion in the simulations of radio-source confusion & CMB noise and measurable radio-sources causes AMI mass estimates to be systematically low.

Tracking the variable jets of V404 Cygni during its 2015 outburst

We present multi-frequency monitoring observations of the black hole X-ray binary V404 Cygni throughout its June 2015 outburst. Our data set includes radio and mm/sub-mm photometry, taken with the Karl G. Jansky Very Large Array, Arc-Minute MicroKelvin Imager Large Array, Sub-millimeter Array, James Clerk Maxwell Telescope, and the Northern Extended Millimetre Array, combined with publicly available infrared, optical, UV, and X-ray measurements. With these data, we report detailed diagnostics of the spectral and variability properties of the jet emission observed during different stages of this outburst. These diagnostics show that emission from discrete jet ejecta dominated the jet emission during the brightest stages of the outburst. We find that the ejecta became fainter, slower, less frequent, and less energetic, before the emission transitioned (over 1--2 days) to being dominated by a compact jet, as the outburst decayed toward quiescence. While the broad-band spectrum of this compact jet showed very little evolution throughout the outburst decay (with the optically thick to thin synchrotron jet spectral break residing in the near-infrared/optical bands; $\sim2-5\times10^{14}$ Hz), the emission still remained intermittently variable at mm/sub-mm frequencies. Additionally, we present a comparison between the radio jet emission throughout the 2015 and previous 1989 outbursts, confirming that the radio emission in the 2015 outburst decayed significantly faster than in 1989. Lastly, we detail our sub-mm observations taken during the December 2015 mini-outburst of V404 Cygni, which demonstrate that, similar to the main outburst, the source was likely launching jet ejecta during this short period of renewed activity.

Free-form modelling of galaxy clusters: a Bayesian and data-driven approach

ABSTRACTA new method is presented for modelling the physical properties of galaxy clusters. Our technique moves away from the traditional approach of assuming specific parameterized functional forms for the variation of physical quantities within the cluster, and instead allows for a ‘free-form’ reconstruction, but one for which the level of complexity is determined automatically by the observational data and may depend on position within the cluster. This is achieved by representing each independent cluster property as some interpolating or approximating function that is specified by a set of control points, or ‘nodes’, for which the number of nodes, together with their positions and amplitudes, are allowed to vary and are inferred in a Bayesian manner from the data. We illustrate our nodal approach in the case of a spherical cluster by modelling the electron pressure profile Pe(r) in analyses both of simulated Sunyaev–Zel’dovich (SZ) data from the Arcminute MicroKelvin Imager (AMI) and of real AMI observations of the cluster MACS J0744+3927 in the CLASH sample. We demonstrate that one may indeed determine the complexity supported by the data in the reconstructed Pe(r), and that one may constrain two very important quantities in such an analysis: the cluster total volume integrated Comptonization parameter (Ytot) and the extent of the gas distribution in the cluster (rmax). The approach is also well-suited to detecting clusters in blind SZ surveys, in the case where the population of radio sources is known in advance.

The nature of the variable millimetre–selected AGN in the brightest cluster galaxy of Abell 851

ABSTRACT We present the detection of a bright 3 mm continuum source in the brightest cluster galaxy (BCG) in Abell 0851 (z = 0.411) with the NOrthern Extended Millimeter Array (NOEMA). When this detection is compared to other multifrequency observations across 21cm–$100 {\mu {m}}$, including new Arcminute Microkelvin Imager 15 GHz observations, we find evidence for a relatively flat, variable core source associated with the BCG. The radio power and amplitude of variability observed in this galaxy is consistent with the cores in lower redshift BCGs in X-ray–selected clusters, and the flat mm–cm spectrum is suggestive of the BCG being a low-luminosity active galactic nucleus archetype. The discovery of this system could provide a basis for a long-term study of the role of low-luminosity radio mode ‘regulatory’ feedback in massive clusters.

Investigation of the cosmic ray population and magnetic field strength in the halo of NGC 891

Context.Cosmic rays and magnetic fields play an important role for the formation and dynamics of gaseous halos of galaxies.Aims.Low-frequency radio continuum observations of edge-on galaxies are ideal to study cosmic-ray electrons (CREs) in halos via radio synchrotron emission and to measure magnetic field strengths. Spectral information can be used to test models of CRE propagation. Free–free absorption by ionized gas at low frequencies allows us to investigate the properties of the warm ionized medium in the disk.Methods.We obtained new observations of the edge-on spiral galaxy NGC 891 at 129–163 MHz with the LOw Frequency ARray (LOFAR) and at 13–18 GHz with the Arcminute Microkelvin Imager (AMI) and combine them with recent high-resolution Very Large Array (VLA) observations at 1–2 GHz, enabling us to study the radio continuum emission over two orders of magnitude in frequency.Results.The spectrum of the integrated nonthermal flux density can be fitted by a power law with a spectral steepening towards higher frequencies or by a curved polynomial. Spectral flattening at low frequencies due to free–free absorption is detected in star-forming regions of the disk. The mean magnetic field strength in the halo is 7 ± 2μG. The scale heights of the nonthermal halo emission at 146 MHz are larger than those at 1.5 GHz everywhere, with a mean ratio of 1.7 ± 0.3, indicating that spectral ageing of CREs is important and that diffusive propagation dominates. The halo scale heights at 146 MHz decrease with increasing magnetic field strengths which is a signature of dominating synchrotron losses of CREs. On the other hand, the spectral index between 146 MHz and 1.5 GHz linearly steepens from the disk to the halo, indicating that advection rather than diffusion is the dominating CRE transport process. This issue calls for refined modelling of CRE propagation.Conclusions.Free–free absorption is probably important at and below about 150 MHz in the disks of edge-on galaxies. To reliably separate the thermal and nonthermal emission components, to investigate spectral steepening due to CRE energy losses, and to measure magnetic field strengths in the disk and halo, wide frequency coverage and high spatial resolution are indispensable.

A comprehensive study of high-energy gamma-ray and radio emission from Cyg X-3

We study high-energy $\gamma$-rays observed from Cyg X-3 by the Fermi Large Area Telescope and the 15-GHz emission observed by the Ryle Telescope and the Arcminute Microkelvin Imager. We measure the $\gamma$-ray spectrum averaged over strong flares much more accurately than before, and find it well modelled by Compton scattering of stellar radiation by relativistic electrons with the power law index of $\simeq$3.5 and a low-energy cutoff at the Lorentz factor of $\sim\!10^3$. We find a weaker spectrum in the soft spectral state, but only upper limits in the hard and intermediate states. We measure strong orbital modulation during the flaring state, well modelled by anisotropic Compton scattering of blackbody photons from the donor by jet relativistic electrons. We discover a weaker orbital modulation of the 15 GHz radio emission, which is well modelled by free-free absorption by the stellar wind. We then study cross-correlations between radio, $\gamma$-ray and X-ray emissions. We find the cross-correlation between the radio and $\gamma$-ray emissions peaks at a lag less than 1 d, while we detect a distinct radio lag of $\sim$50 d with respect to the soft X-rays in the soft spectral state.