Source author record

I. Heywood

I. Heywood appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

54works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

54 published item(s)

preprint2025arXiv

Atomic hydrogen reservoirs in quiescent galaxies at z = 0.4

Context. Based on Local Universe observations, quiescent galaxies (QGs) host lower to no HI compared to star-forming galaxies (SFGs), but no constraints have been derived so far at higher redshift (z>0.1). Understanding whether QGs can retain significant HI reservoirs at higher z is crucial to refine quenching and gas accretion models and to constrain overall star formation efficiency at different epochs. Aims. We aim to probe HI in candidate QGs at intermediate redshifts (z=0.36) and to understand whether there exists a class of QGs retaining consistent HI reservoirs and which parameters (dust content, stellar mass, Dn4000, morphology, environment) effectively capture HI-rich QGs. Methods. We perform 21-cm spectral line stacking on MIGHTEE-HI data at z=0.36, targeting two different samples of QGs, defined by means of a color-selection criterion and a spectroscopic criterion based on Dn4000, respectively. We also perform stacking on subsamples of the spectroscopically-selected quiescent sample to investigate the correlation between the HI content and other galaxy properties. Results. We find that QGs with an IR counterpart (i.e., dusty galaxies) are found to host a substantial HI content, on average just 40% lower than SFGs. In contrast, color-selected QGs still hold HI, but lower than SFGs by a factor 3. Among dusty objects, we find morphology to have a mild impact on the atomic gas content, with spirals hosting approximately 15-30% more HI than spheroids. Environmental effects are also present, with low-density regions hosting galaxies that are HI-richer than in high-density ones, by approximately 30% for spirals and 60% for spheroids. We suggest that, in general, HI content is driven by several factors, including quenching mechanisms and ISM enrichment processes. Also, quiescent galaxies - and especially dusty systems - seem to yield HI more consistently than in the Local Universe.

preprint2022arXiv

21 new long-term variables in the GX 339-4 field: two years of MeerKAT monitoring

We present 21 new long-term variable radio sources found commensally in two years of weekly MeerKAT monitoring of the low-mass X-ray binary GX 339-4. The new sources vary on time scales of weeks to months and have a variety of light curve shapes and spectral index properties. Three of the new variable sources are coincident with multi-wavelength counterparts; and one of these is coincident with an optical source in deep MeerLICHT images. For most sources, we cannot eliminate refractive scintillation of active galactic nuclei as the cause of the variability. These new variable sources represent $2.2\pm0.5$ per cent of the unresolved sources in the field, which is consistent with the 1-2 per cent variability found in past radio variability surveys. However, we expect to find short-term variable sources in the field as well as these 21 new long-term variable sources. We present the radio light curves and spectral index variability of the new variable sources, as well as the absolute astrometry and matches to coincident sources at other wavelengths.

preprint2022arXiv

Evidence for X-ray Emission in Excess to the Jet Afterglow Decay 3.5 yrs After the Binary Neutron Star Merger GW 170817: A New Emission Component

For the first $\sim3$ years after the binary neutron star merger event GW 170817 the radio and X-ray radiation has been dominated by emission from a structured relativistic off-axis jet propagating into a low-density medium with n $< 0.01\,\rm{cm^{-3}}$. We report on observational evidence for an excess of X-ray emission at $δt>900$ days after the merger. With $L_x\approx5\times 10^{38}\,\rm{erg\,s^{-1}}$ at 1234 days, the recently detected X-ray emission represents a $\ge 3.2\,σ$ (Gaussian equivalent) deviation from the universal post jet-break model that best fits the multi-wavelength afterglow at earlier times. In the context of JetFit afterglow models, current data represent a departure with statistical significance $\ge 3.1\,σ$, depending on the fireball collimation, with the most realistic models showing excesses at the level of $\ge 3.7\,σ$. A lack of detectable 3 GHz radio emission suggests a harder broad-band spectrum than the jet afterglow. These properties are consistent with the emergence of a new emission component such as synchrotron radiation from a mildly relativistic shock generated by the expanding merger ejecta, i.e. a kilonova afterglow. In this context, we present a set of ab-initio numerical-relativity BNS merger simulations that show that an X-ray excess supports the presence of a high-velocity tail in the merger ejecta, and argues against the prompt collapse of the merger remnant into a black hole. Radiation from accretion processes on the compact-object remnant represents a viable alternative. Neither a kilonova afterglow nor accretion-powered emission have been observed before, as detections of BNS mergers at this phase of evolution are unprecedented.

preprint2022arXiv

MIGHTEE-HI: The relation between the HI gas in galaxies and the cosmic web

We study the 3D axis of rotation (3D spin) of 77 HI galaxies from the MIGHTEE-HI Early Science observations, and its relation to the filaments of the cosmic web. For this HI-selected sample, the alignment between the spin axis and the closest filament ($\lvert \cos ψ\rvert$) is higher for galaxies closer to the filaments, with $\langle\lvert \cos ψ\rvert\rangle= 0.66 \pm 0.04$ for galaxies $<5$ Mpc from their closest filament compared to $\langle\lvert \cos ψ\rvert\rangle= 0.37 \pm 0.08$ for galaxies at $5 < d <10$ Mpc. We find that galaxies with a low HI-to-stellar mass ratio ($\log_{10}(M_{\rm HI}/M_{\star}) < 0.11$) are more aligned with their closest filaments, with $\langle\lvert \cos ψ\rvert\rangle= 0.58 \pm 0.04$; whilst galaxies with ($\log_{10}(M_{\rm HI}/M_{\star}) > 0.11$) tend to be mis-aligned, with $\langle\lvert \cos ψ\rvert\rangle= 0.44 \pm 0.04$. We find tentative evidence that the spin axis of HI-selected galaxies tend to be aligned with associated filaments ($d<10$ Mpc), but this depends on the gas fractions. Galaxies that have accumulated more stellar mass compared to their gas mass tend towards stronger alignment. Our results suggest that those galaxies that have accrued high gas fraction with respect to their stellar mass may have had their spin axis alignment with the filament disrupted by a recent gas-rich merger, whereas the spin vector for those galaxies in which the neutral gas has not been strongly replenished through a recent merger tend to orientate towards alignment with the filament. We also investigate the spin transition between galaxies with a high HI content and a low HI content at a threshold of $M_{\mathrm{HI}}\approx 10^{9.5} M_{\odot}$ found in simulations, however we find no evidence for such a transition with the current data.

preprint2022arXiv

MIGHTEE: the nature of the radio-loud AGN population

We study the nature of the faint radio source population detected in the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) Early Science data in the COSMOS field, focusing on the properties of the radio-loud active galactic nuclei (AGN). Using the extensive multi-wavelength data available in the field, we are able to classify 88 per cent of the 5223 radio sources in the field with host galaxy identifications as AGN (35 per cent) or star-forming galaxies (54 per cent). We select a sample of radio-loud AGN with redshifts out to $z \sim 6$ and radio luminosities $10^{20} < \textrm{L}_{1.4~\textrm{GHz}} / \textrm{W Hz}^{-1} < 10^{27}$ and classify them as high-excitation and low-excitation radio galaxies (HERGs and LERGs). The classification catalogue is released with this work. We find no significant difference in the host galaxy properties of the HERGs and LERGs in our sample. In contrast to previous work, we find that the HERGs and LERGs have very similar Eddington-scaled accretion rates; in particular we identify a population of very slowly accreting AGN that are formally classified as HERGs at these low radio luminosities, where separating into HERGs and LERGs possibly becomes redundant. We investigate how black hole mass affects jet power, and find that a black hole mass $\gtrsim 10^{7.8}~\textrm{M}_\odot$ is required to power a jet with mechanical power greater than the radiative luminosity of the AGN ($L_\textrm{mech}/L_\textrm{bol} > 1$). We discuss that both a high black hole mass and black hole spin may be necessary to launch and sustain a dominant radio jet.

preprint2022arXiv

Radio footprints of a minor merger in the Shapley Supercluster: From supercluster down to galactic scales

The Shapley Supercluster ($\langle z \rangle\approx0.048$) contains several tens of gravitationally bound clusters and groups, making it it is an ideal subject for radio studies of cluster mergers. We used new high sensitivity radio observations to investigate the less energetic events of mass assembly in the Shapley Supercluster from supercluster down to galactic scales. We created total intensity images of the full region between A 3558 and A 3562, from $\sim 230$ to $\sim 1650$ MHz, using ASKAP, MeerKAT and the GMRT, with sensitivities ranging from $\sim 6$ to $\sim 100$ $μ$Jy beam$^{-1}$. We performed a detailed morphological and spectral study of the extended emission features, complemented with ESO-VST optical imaging and X-ray data from XMM-Newton. We report the first GHz frequency detection of extremely low brightness intercluster diffuse emission on a $\sim 1$ Mpc scale connecting a cluster and a group, namely: A 3562 and the group SC 1329--313. It is morphologically similar to the X-ray emission in the region. We also found (1) a radio tail generated by ram pressure stripping in the galaxy SOS 61086 in SC 1329-313; (2) a head-tail radio galaxy, whose tail is broken and culminates in a misaligned bar; (3) ultrasteep diffuse emission at the centre of A 3558. Finally (4), we confirm the ultra-steep spectrum nature of the radio halo in A 3562. Our study strongly supports the scenario of a flyby of SC 1329-313 north of A 3562 into the supercluster core. [abridged...]

preprint2022arXiv

Search and identification of transient and variable radio sources using MeerKAT observations: a case study on the MAXI J1820+070 field

Many transient and variable sources detected at multiple wavelengths are also observed to vary at radio frequencies. However, these samples are typically biased towards sources that are initially detected in wide-field optical, X-ray or gamma-ray surveys. Many sources that are insufficiently bright at higher frequencies are therefore missed, leading to potential gaps in our knowledge of these sources and missing populations that are not detectable in optical, X-rays or gamma-rays. Taking advantage of new state-of-the-art radio facilities that provide high quality wide-field images with fast survey speeds, we can now conduct unbiased surveys for transient and variable sources at radio frequencies. In this paper, we present an unbiased survey using observations obtained by MeerKAT, a mid-frequency ($\sim$1.4 GHz) radio array in South Africa's Karoo Desert. The observations used were obtained as part of a weekly monitoring campaign for X-ray binaries (XRBs) and we focus on the field of MAXI J1820+070. We develop methods to optimally filter transient and variable candidates that can be directly applied to other datasets. In addition to MAXI J1820+070, we identify four likely active galactic nuclei, one source that could be a Galactic source (pulsar or quiescent X-ray binary) or an AGN, and one variable pulsar. No transient sources, defined as being undetected in deep images, were identified leading to a transient surface density of $<3.7\times10^{-2}$ deg$^{-2}$ at a sensitivity of 1 mJy on timescales of one week at 1.4 GHz.

preprint2022arXiv

Statistical Properties of the Population of the Galactic Center Filaments II: The Spacing Between Filaments

We carry out a population study of magnetized radio filaments in the Galactic center using MeerKAT data by focusing on the spacing between the filaments that are grouped. The morphology of a sample of 43 groupings containing 174 magnetized radio filaments are presented. Many grouped filaments show harp-like, fragmented cometary tail-like, or loop-like structures in contrast to many straight filaments running mainly perpendicular to the Galactic plane. There are many striking examples of a single filament splitting into two prongs at a junction, suggestive of a flow of plasma along the filaments. Spatial variations in spectral index, brightness, bending and sharpening along the filaments indicate that they are evolving on a 10^{5-6}-year time scale. The mean spacings between parallel filaments in a given grouping peaks at $\sim16''$. We argue by modeling that the filaments in a grouping all lie on the same plane and that the groupings are isotropically oriented in 3D space. One candidate for the origin of filamentation is interaction with an obstacle, which could be a compact radio source, before a filament splits and bends into multiple filaments. In this picture, the obstacle or sets the length scale of the separation between the filaments. Another possibility is synchrotron cooling instability occurring in cometary tails formed as a result of the interaction of cosmic-ray driven Galactic center outflow with obstacles such as stellar winds. In this picture, the mean spacing and the mean width of the filaments are expected to be a fraction of a parsec, consistent with observed spacing.

preprint2022arXiv

Statistical Properties of the Population of the Galactic Center Filaments: The Spectral Index and Equipartition Magnetic Field

We present high-pass filtered continuum images of the inner $3.5^\circ\times2.5^\circ$ of the Galactic center at 20 cm with $6.4''$ resolution. These mosaic images are taken with MeerKAT and reveal a large number of narrow filaments, roughly an order of magnitude increase in their numbers compared to past measurements. For the first time, we carry out population studies of the spectral index and magnetic field of the entire region. The mean spectral indices of the filaments are steeper than supernova remnants (SNRs) (-0.62) with a value of $α\sim-0.83$. The variation in $α$ is much larger than for the SNRs, suggesting that these characteristics have a different origin. A large-scale cosmic-ray driven wind has recently been proposed to explain the origin of filaments and the large-scale 430 pc bipolar radio and X-ray structure. This favors the possibility that the large-scale bipolar radio/X-ray structure is produced by past activity of Sgr A* rather than coordinated burst of supernovae. A trend of steeper indices is also noted with increasing distance from the Galactic plane. This could be explained either by synchrotron cooling or weak shocks accelerating cosmic-ray particles in the context of the cosmic-ray driven wind. The mean magnetic field strengths along the filaments ranges from $\sim100$ to 400 $μ$G depending on the assumed ratio of cosmic-ray protons to electrons. Given that there is a high cosmic ray pressure in the Galactic center, the large equipartition magnetic field implies that the magnetic field is weak in most of the interstellar volume of the Galactic center.

preprint2022arXiv

The 1.28 GHz MeerKAT Galactic Center Mosaic

The inner $\sim$200 pc region of the Galaxy contains a 4 million M$_{\odot}$ supermassive black hole (SMBH), significant quantities of molecular gas, and star formation and cosmic ray energy densities that are roughly two orders of magnitude higher than the corresponding levels in the Galactic disk. At a distance of only 8.2 kpc, the region presents astronomers with a unique opportunity to study a diverse range of energetic astrophysical phenomena, from stellar objects in extreme environments, to the SMBH and star-formation driven feedback processes that are known to influence the evolution of galaxies as a whole. We present a new survey of the Galactic center conducted with the South African MeerKAT radio telescope. Radio imaging offers a view that is unaffected by the large quantities of dust that obscure the region at other wavelengths, and a scene of striking complexity is revealed. We produce total intensity and spectral index mosaics of the region from 20 pointings (144 hours on-target in total), covering 6.5 square degrees with an angular resolution of 4$"$,at a central frequency of 1.28 GHz. Many new features are revealed for the first time due to a combination of MeerKAT's high sensitivity, exceptional $u,v$-plane coverage, and geographical vantage point. We highlight some initial survey results, including new supernova remnant candidates, many new non-thermal filament complexes, and enhanced views of the Radio Arc Bubble, Sgr A and Sgr B regions. This project is a SARAO public legacy survey, and the image products are made available with this article.

preprint2022arXiv

The Curious Case of the "Heartworm" Nebula

The curious Galactic features near G357.2$-$0.2 were observed with the MeerKAT radio interferometer array in the UHF and L bands (0.56--1.68 GHz). There are two possibly related features: a newly identified faint heart-shaped partial shell (the "Heart"), and a series of previously known but now much better imaged narrow, curved features (the "Worm") interior to the heart. Polarized emission suggests that much of the emission is nonthermal and is embedded in a dense plasma. The filaments of the worm appear to be magnetic structures powered by embedded knots that are sites of particle acceleration. The morphology of the worm broadly resembles some known pulsar wind nebulae (PWNe) but there is no known pulsar or PWN which could be powering this structure. We also present eROSITA observations of the field; no part of the nebula is detected in X-rays, but the current limits do not preclude the existence of a pulsar/PWN of intermediate spin-down luminosity.

preprint2021arXiv

A MeerKAT, e-MERLIN, H.E.S.S. and Swift search for persistent and transient emission associated with three localised FRBs

We report on a search for persistent radio emission from the one-off Fast Radio Burst (FRB) 20190714A, as well as from two repeating FRBs, 20190711A and 20171019A, using the MeerKAT radio telescope. For FRB 20171019A we also conducted simultaneous observations with the High Energy Stereoscopic System (H.E.S.S.) in very high energy gamma rays and searched for signals in the ultraviolet, optical, and X-ray bands. For this FRB, we obtain a UV flux upper limit of 1.39x10^-16 erg/cm^-2/s/Amstrong, X-ray limit of ~ 6.6x10^-14 erg/cm^-2/s and a limit on the very-high-energy gamma-ray flux (Phi) (E > 120 GeV) < 1.7 x 10^-12 erg/cm^-2/s. We obtain a radio upper limit of ~15 microJy/beam for persistent emission at the locations of both FRBs 20190711A and 20171019A, but detect diffuse radio emission with a peak brightness of ~53 microJy/beam associated with FRB 20190714A at z = 0.2365. This represents the first detection of the radio continuum emission potentially associated with the host (galaxy) of FRB 20190714A, and is only the third known FRB to have such an association. Given the possible association of a faint persistent source, FRB 20190714A may potentially be a repeating FRB whose age lies between that of FRB 20121102A and FRB 20180916A. A parallel search for repeat bursts from these FRBs revealed no new detections down to a fluence of 0.08 Jy ms for a 1 ms duration burst.

preprint2021arXiv

Observations of a radio-bright, X-ray obscured GRS 1915+105

The Galactic black hole transient GRS1915+105 is famous for its markedly variable X-ray and radio behaviour, and for being the archetypal galactic source of relativistic jets. It entered an X-ray outburst in 1992 and has been active ever since. Since 2018 GRS1915+105 has declined into an extended low-flux X-ray plateau, occasionally interrupted by multi-wavelength flares. Here we report the radio and X-ray properties of GRS1915+105 collected in this new phase, and compare the recent data to historic observations. We find that while the X-ray emission remained unprecedentedly low for most of the time following the decline in 2018, the radio emission shows a clear mode change half way through the extended X-ray plateau in 2019 June: from low flux (~3mJy) and limited variability, to marked flaring with fluxes two orders of magnitude larger. GRS1915+105 appears to have entered a low-luminosity canonical hard state, and then transitioned to an unusual accretion phase, characterised by heavy X-ray absorption/obscuration. Hence, we argue that a local absorber hides from the observer the accretion processes feeding the variable jet responsible for the radio flaring. The radio-X-ray correlation suggests that the current low X-ray flux state may be a signature of a super-Eddington state akin to the X-ray binaries SS433 or V404 Cyg.

preprint2021arXiv

Radio and optical observations of the possible AE Aqr twin, LAMOST J024048.51+195226.9

Thorstensen (2020) recently argued that the cataclysmic variable (CV) LAMOST J024048.51+195226.9 may be a twin to the unique magnetic propeller system AE Aqr. If this is the case, two predictions are that it should display a short period white dwarf spin modulation, and that it should be a bright radio source. We obtained follow-up optical and radio observations of this CV, in order to see if this holds true. Our optical high-speed photometry does not reveal a white dwarf spin signal, but lacks the sensitivity to detect a modulation similar to the 33-s spin signal seen in AE Aqr. We detect the source in the radio, and measure a radio luminosity similar to that of AE Aqr and close to the highest so far reported for a CV. We also find good evidence for radio variability on a time scale of tens of minutes. Optical polarimetric observations produce no detection of linear or circular polarization. While we are not able to provide compelling evidence, our observations are all consistent with this object being a propeller system.

preprint2021arXiv

The Galactic center chimneys: The base of the multiphase outflow of the Milky Way

Outflows and feedback are key ingredients of galaxy evolution. Evidence for an outflow arising from the Galactic center (GC) has recently been discovered at different wavelength. We show that the X-ray, radio, and infrared emissions are deeply interconnected, affecting one another and forming coherent features on scales of hundreds of parsecs, therefore indicating a common physical link associated with the GC outflow. We debate the location of the northern chimney and suggest that it might be located on the front side of the GC because of a significant tilt of the chimneys toward us. We report the presence of strong shocks at the interface between the chimneys and the interstellar medium, which are traced by radio and warm dust emission. We observe entrained molecular gas outflowing within the chimneys, revealing the multiphase nature of the outflow. In particular, the molecular outflow produces a long, strong, and structured shock along the northwestern wall of the chimney. Because of the different dynamical times of the various components of the outflow, the chimneys appear to be shaped by directed large-scale winds launched at different epochs. The data support the idea that the chimneys are embedded in an (often dominant) vertical magnetic field, which likely diverges with increasing latitude. We observe that the thermal pressure associated with the hot plasma appears to be smaller than the ram pressure of the molecular outflow and the magnetic pressure. This leaves open the possibility that either the main driver of the outflow is more powerful than the observed hot plasma, or the chimneys represent a "relic" of past and more powerful activity. These multiwavelength observations corroborate the idea that the chimneys represent the channel connecting the quasi-continuous, but intermittent, activity at the GC with the base of the Fermi bubbles.

preprint2021arXiv

The MeerKAT Galaxy Cluster Legacy Survey I. Survey Overview and Highlights

MeerKAT's large number of antennas, spanning 8 km with a densely packed 1 km core, create a powerful instrument for wide-area surveys, with high sensitivity over a wide range of angular scales. The MeerKAT Galaxy Cluster Legacy Survey (MGCLS) is a programme of long-track MeerKAT L-band (900-1670 MHz) observations of 115 galaxy clusters, observed for $\sim$6-10 hours each in full polarisation. The first legacy product data release (DR1), made available with this paper, includes the MeerKAT visibilities, basic image cubes at $\sim$8" resolution, and enhanced spectral and polarisation image cubes at $\sim$8" and 15" resolutions. Typical sensitivities for the full-resolution MGCLS image products are $\sim$3-5 μJy/beam. The basic cubes are full-field and span 4 deg^2. The enhanced products consist of the inner 1.44 deg^2 field of view, corrected for the primary beam. The survey is fully sensitive to structures up to $\sim$10' scales and the wide bandwidth allows spectral and Faraday rotation mapping. HI mapping at 209 kHz resolution can be done at $0<z<0.09$ and $0.19<z<0.48$. In this paper, we provide an overview of the survey and DR1 products, including caveats for usage. We present some initial results from the survey, both for their intrinsic scientific value and to highlight the capabilities for further exploration with these data. These include a primary beam-corrected compact source catalogue of $\sim$626,000 sources for the full survey, and an optical/infrared cross-matched catalogue for compact sources in Abell 209 and Abell S295. We examine dust unbiased star-formation rates as a function of clustercentric radius in Abell 209 and present a catalogue of 99 diffuse cluster sources (56 are new), some of which have no suitable characterisation. We also highlight some of the radio galaxies which challenge current paradigms and present first results from HI studies of four targets.

preprint2020arXiv

A MeerKAT Survey of Nearby Novalike Cataclysmic Variables

We present the results of MeerKAT radio observations of eleven nearby novalike cataclysmic variables. We have detected radio emission from IM Eri, RW Sex, V3885 Sgr and V603 Aql. While RW Sex, V3885 Sgr and V603 Aql had been previously detected, this is the first reported radio detection of IM Eri. Our observations have doubled the sample of non-magnetic CVs with sensitive radio data. We observe that at our radio detection limits, a specific optical luminosity $\gtrsim 2.2\times 10^{18}$ erg/s/Hz (corresponding to $M_V \lesssim 6.0$) is required to produce a radio detection. We also observe that the X-ray and radio luminosities of our detected novalikes are on an extension of the $L_{X} \propto L_{R}^{\sim 0.7}$ power law originally proposed for non-pulsating neutron star low-mass X-ray binaries. We find no other correlations between the radio emission and emission in other wavebands or any other system parameters for the existing sample of radio-detected non-magnetic CVs. We measure in-band (0.9-1.7 GHz) radio spectral indices that are consistent with reports from earlier work. Finally, we constructed broad spectral energy distributions for our sample from published multi-wavelength data, and use them to place constraints on the mass transfer rates of these systems.

preprint2020arXiv

An extremely powerful long-lived superluminal ejection from the black hole MAXI J1820+070

Black holes in binary systems execute patterns of outburst activity where two characteristic X-ray states are associated with different behaviours observed at radio wavelengths. The hard state is associated with radio emission indicative of a continuously replenished, collimated, relativistic jet, whereas the soft state is rarely associated with radio emission, and never continuously, implying the absence of a quasi-steady jet. Here we report radio observations of the black hole transient MAXI J1820$+$070 during its 2018 outburst. As the black hole transitioned from the hard to soft state we observed an isolated radio flare, which, using high angular resolution radio observations, we connect with the launch of bi-polar relativistic ejecta. This flare occurs as the radio emission of the core jet is suppressed by a factor of over 800. We monitor the evolution of the ejecta over 200 days and to a maximum separation of 10$''$, during which period it remains detectable due to in-situ particle acceleration. Using simultaneous radio observations sensitive to different angular scales we calculate an accurate estimate of energy content of the approaching ejection. This energy estimate is far larger than that derived from state transition radio flare, suggesting a systematic underestimate of jet energetics.

preprint2020arXiv

Evidence for a jet and outflow from Sgr A*: a continuum and spectral line study

We study the environment of Sgr A* using spectral and continuum observations with the ALMA and VLA. Our analysis of sub-arcsecond H30alpha, H39alpha, H52alpha and H56alpha line emission towards Sgr A* confirm the recently published broad peak ~500 km/s~spectrum toward Sgr~A*. We also detect emission at more extreme radial velocities peaking near -2500 and 4000 km/s, within 0.2''. We then present broad band radio continuum images at multiple frequencies on scales from arcseconds to arcminutes. A number of elongated continuum structures lie parallel to the Galactic plane, extending from ~0.4'' to 10'. We note a nonthermal elongated structure on an arcminute scale emanating from Sgr A* at low frequencies between 1 and 1.4 GHz where thermal emission from the mini-spiral is depressed by optical depth effects. The position angle of this elongated structure and the sense of motion of ionized features with respect to Sgr A* suggest a symmetric, collimated jet emerging from Sgr A* with an opening angle of ~30deg and a position angle of ~60deg punching through the medium before accelerating a significant fraction of the orbiting ionized gas to high velocities. The jet with estimated mass flow rate ~1.4x10^{-5} solar mass/yr emerges perpendicular to the equatorial plane of the accretion flow near the event horizon of Sgr A* and runs along the Galactic plane. To explain a number of east-west features near Sgr A*, we also consider the possibility of an outflow component with a wider-angle launched from the accretion flow at larger radii.

preprint2020arXiv

Field sources near the southern-sky calibrator PKS B1934-638: effect on spectral line observations with SKA-MID and its precursors

Accurate instrumental bandpass corrections are essential for the reliable interpretation of spectral lines from targeted and survey-mode observations with radio interferometers. Bandpass correction is typically performed by comparing measurements of a strong calibrator source to an assumed model, typically an isolated point source. The wide field-of-view and high sensitivity of modern interferometers means that additional sources are often detected in observations of calibrators. This can introduce errors into bandpass corrections and subsequently the target data if not properly accounted for. Focusing on the standard calibrator PKS B1934-638, we perform simulations to asses this effect by constructing a wide-field sky model. The cases of ASKAP (0.7-1.9 GHz), MeerKAT (UHF: 0.58-1.05 GHz; L-band: 0.87-1.67 GHz) and Band 2 (0.95-1.76 GHz) of SKA-MID are examined. The use of a central point source model during bandpass calibration is found to impart amplitude errors into spectra measured by the precursor instruments at the ~0.2-0.5% level dropping to ~0.01% in the case of SKA-MID. This manifests itself as ripples in the source spectrum, the behaviour of which is coupled to the distribution of the array baselines, the solution interval, the primary beam size, the hour-angle of the calibration scan, as well as the weights used when imaging the target. Calibration pipelines should routinely employ complete field models for standard calibrators to remove this potentially destructive contaminant from the data, a recommendation we validate by comparing our simulation results to a MeerKAT scan of PKS B1934-638, calibrated with and without our expanded sky model.

preprint2020arXiv

MIGHTEE-HI: The HI emission project of the MeerKAT MIGHTEE survey

We present the HI emission project within the MIGHTEE survey, currently being carried out with the newly commissioned MeerKAT radio telescope. This is one of the first deep, blind, medium-wide interferometric surveys for neutral hydrogen (HI) ever undertaken, extending our knowledge of HI emission to z=0.6. The science goals of this medium-deep, medium-wide survey are extensive, including the evolution of the neutral gas content of galaxies over the past 5 billion years. Simulations predict nearly 3000 galaxies over 0<z<0.4 will be detected directly in HI, with statistical detections extending to z=0.6. The survey allows us to explore HI as a function of galaxy environment, with massive groups and galaxy clusters within the survey volume. Additionally, the area is large enough to contain as many as 50 local galaxies with HI mass $<10^8$ Msun, which allows us to study the low-mass galaxy population. The 20 deg$^2$ main survey area is centred on fields with exceptional multi-wavelength ancillary data, with photometry ranging from optical through far-infrared wavelengths, supplemented with multiple spectroscopic campaigns. We describe here the survey design and the key science goals. We also show first results from the Early Science observations, including kinematic modelling of individual sources, along with the redshift, HI, and stellar mass ranges of the sample to date.

preprint2020arXiv

Observation of inverse Compton emission from a long $γ$-ray burst

Long-duration gamma-ray bursts (GRBs) originate from ultra-relativistic jets launched from the collapsing cores of dying massive stars. They are characterised by an initial phase of bright and highly variable radiation in the keV-MeV band that is likely produced within the jet and lasts from milliseconds to minutes, known as the prompt emission. Subsequently, the interaction of the jet with the external medium generates external shock waves, responsible for the afterglow emission, which lasts from days to months, and occurs over a broad energy range, from the radio to the GeV bands. The afterglow emission is generally well explained as synchrotron radiation by electrons accelerated at the external shock. Recently, an intense, long-lasting emission between 0.2 and 1 TeV was observed from the GRB 190114C. Here we present the results of our multi-frequency observational campaign of GRB~190114C, and study the evolution in time of the GRB emission across 17 orders of magnitude in energy, from $5\times10^{-6}$ up to $10^{12}$\,eV. We find that the broadband spectral energy distribution is double-peaked, with the TeV emission constituting a distinct spectral component that has power comparable to the synchrotron component. This component is associated with the afterglow, and is satisfactorily explained by inverse Compton upscattering of synchrotron photons by high-energy electrons. We find that the conditions required to account for the observed TeV component are not atypical, supporting the possibility that inverse Compton emission is commonly produced in GRBs.

preprint2020arXiv

Radio & X-ray detections of GX 339--4 in quiescence using MeerKAT and Swift

The radio:X-ray correlation that characterises accreting black holes at all mass scales - from stellar mass black holes in binary systems to super-massive black holes powering Active Galactic Nuclei - is one of the most important pieces of observational evidence supporting the existence of a connection between the accretion process and the generation of collimated outflows - or jets - in accreting systems. Although recent studies suggest that the correlation extends down to low luminosities, only a handful of stellar mass black holes have been clearly detected, and in general only upper limits (especially at radio wavelengths) can be obtained during quiescence. We recently obtained detections of the black hole X-ray binary GX 339--4 in quiescence using the MeerKAT radio telescope and Swift X-ray Telescope instrument onboard the Neil Gehrels Swift Observatory, probing the lower end of the radio:X-ray correlation. We present the properties of accretion and of the connected generation of jets in the poorly studied low-accretion rate regime for this canonical black hole XRB system.

preprint2020arXiv

The relation between the diffuse X-ray luminosity and the radio power of the central AGN in galaxy groups

Our understanding of how AGN feedback operates in galaxy clusters has improved in recent years owing to large efforts in multi-wavelength observations and hydrodynamical simulations. However, it is much less clear how feedback operates in galaxy groups, which have shallower gravitational potentials. In this work, using very deep VLA and new MeerKAT observations from the MIGHTEE survey, we compiled a sample of 247 X-ray selected galaxy groups detected in the COSMOS field. We have studied the relation between the X-ray emission of the intra-group medium and the 1.4 GHz radio emission of the central radio galaxy. For comparison, we have also built a control sample of 142 galaxy clusters using ROSAT and NVSS data. We find that clusters and groups follow the same correlation between X-ray and radio emission. Large radio galaxies hosted in the centres of groups and merging clusters increase the scatter of the distribution. Using statistical tests and Monte-Carlo simulations, we show that the correlation is not dominated by biases or selection effects. We also find that galaxy groups are more likely than clusters to host large radio galaxies, perhaps owing to the lower ambient gas density or a more efficient accretion mode. In these groups, radiative cooling of the ICM could be less suppressed by AGN heating. We conclude that the feedback processes that operate in galaxy clusters are also effective in groups.

preprint2020arXiv

VLA imaging of the XMM-LSS / VIDEO deep field at 1-2 GHz

Modern radio telescopes are routinely reaching depths where normal starforming galaxies are the dominant observed population. Realising the potential of radio as a tracer of star formation and black hole activity over cosmic time involves achieving such depths over representative volumes, with radio forming part of a larger multiwavelength campaign. In pursuit of this we used the Karl G. Jansky Very Large Array (VLA) to image $\sim$5 deg$^{2}$ of the VIDEO/XMM-LSS extragalactic deep field at 1--2 GHz. We achieve a median depth of 16 $μ$Jy beam$^{-1}$ with an angular resolution of 4.5\arcsec. Comparisons with existing radio observations of XMM-LSS showcase the improved survey speed of the upgraded VLA: we cover 2.5 times the area and increase the depth by $\sim$20\% in 40\% of the time. Direction-dependent calibration and wide-field imaging were required to suppress the error patterns from off-axis sources of even modest brightness. We derive a catalogue containing 5,762 sources from the final mosaic. Sub-band imaging provides in-band spectral indices for 3,458 (60\%) sources, with the average spectrum becoming flatter than the canonical synchrotron slope below 1 mJy. Positional and flux-density accuracy of the observations, and the differential source counts are in excellent agreement with those of existing measurements. A public release of the images and catalogue accompanies this article.

preprint2019arXiv

The 1.28 GHz MeerKAT DEEP2 Image

We present the confusion-limited 1.28 GHz MeerKAT DEEP2 image covering one $\approx 68'$ FWHM primary beam area with $7.6''$ FWHM resolution and $0.55 \pm 0.01$ $μ$Jy/beam rms noise. Its J2000 center position $α=04^h 13^m 26.4^s$, $δ=-80^\circ 00' 00''$ was selected to minimize artifacts caused by bright sources. We introduce the new 64-element MeerKAT array and describe commissioning observations to measure the primary beam attenuation pattern, estimate telescope pointing errors, and pinpoint $(u,v)$ coordinate errors caused by offsets in frequency or time. We constructed a 1.4 GHz differential source count by combining a power-law count fit to the DEEP2 confusion $P(D)$ distribution from $0.25$ to $10$ $μ$Jy with counts of individual DEEP2 sources between $10$ $μ$Jy and $2.5$ mJy. Most sources fainter than $S \sim 100$ $μ$Jy are distant star-forming galaxies obeying the FIR/radio correlation, and sources stronger than $0.25$ $μ$Jy account for $\sim93\%$ of the radio background produced by star-forming galaxies. For the first time, the DEEP2 source count has reached the depth needed to reveal the majority of the star formation history of the universe. A pure luminosity evolution of the 1.4 GHz local luminosity function consistent with the Madau & Dickinson (2014) model for the evolution of star-forming galaxies based on UV and infrared data underpredicts our 1.4 GHz source count in the range $-5 \lesssim \log[S(\mathrm{Jy})] \lesssim -4$.

preprint2016arXiv

A deep / wide 1-2 GHz snapshot survey of SDSS Stripe 82 using the Karl G. Jansky Very Large Array in a compact hybrid configuration

We have used the Karl G. Jansky Very Large Array to image ~100 sq. deg. of SDSS Stripe 82 at 1-2 GHz. The survey consists of 1,026 snapshot observations of 2.5 minutes duration, using the hybrid CnB configuration. The survey has good sensitivity to diffuse, low surface brightness structures and extended radio emission, making it highly synergistic with existing 1.4 GHz radio observations of the region. The principal data products are continuum images, with 16 x 10 arcsecond resolution, and a catalogue containing 11,782 point and Gaussian components resulting from fits to the thresholded Stokes-I brightness distribution, forming approximately 8,948 unique radio sources. The typical effective 1σ noise level is 88 μJy / beam. Spectral index estimates are included, as derived from the 1 GHz of instantaneous bandwidth. Astrometric and photometric accuracy are in excellent agreement with existing narrowband observations. A large-scale simulation is used to investigate clean bias, which we extend into the spectral domain. Clean bias remains an issue for snapshot surveys with the VLA, affecting our total intensity measurements at the ~1σ level. Statistical spectral index measurements are in good agreement with existing measurements derived from matching separate surveys at two frequencies. At flux densities below ~35σ the median in-band spectral index measurements begin to exhibit a bias towards flatness that is dependent on both flux density and the intrinsic spectral index. In-band spectral curvature measurements are likely to be unreliable for all but the very brightest components. Image products and catalogues are publicly available via an FTP server.

preprint2016arXiv

ATCA detections of massive molecular gas reservoirs in dusty, high-z radio galaxies

Observations using the 7 mm receiver system on the Australia Telescope Compact Array have revealed large reservoirs of molecular gas in two high-redshift radio galaxies: HATLAS J090426.9+015448 (z = 2.37) and HATLAS J140930.4+003803 (z = 2.04). Optically the targets are very faint, and spectroscopy classifies them as narrow-line radio galaxies. In addition to harbouring an active galactic nucleus the targets share many characteristics of sub-mm galaxies. Far-infrared data from Herschel-ATLAS suggest high levels of dust (>10^9 M_solar) and a correspondingly large amount of obscured star formation (~1000 M_solar / yr). The molecular gas is traced via the J = 1-0 transition of 12CO, its luminosity implying total H_2 masses of (1.7 +/- 0.3) x 10^11 and (9.5 +/- 2.4) x 10^10 (alpha_CO/0.8) M_solar in HATLAS J090426.9+015448 and HATLAS J140930.4+003803 respectively. Both galaxies exhibit molecular line emission over a broad (~1000 km/s) velocity range, and feature double-peaked profiles. We interpret this as evidence of either a large rotating disk or an on-going merger. Gas depletion timescales are ~100 Myr. The 1.4 GHz radio luminosities of our targets place them close to the break in the luminosity function. As such they represent `typical' z > 2 radio sources, responsible for the bulk of the energy emitted at radio wavelengths from accretion-powered sources at high redshift, and yet they rank amongst the most massive systems in terms of molecular gas and dust content. We also detect 115 GHz rest-frame continuum emission, indicating a very steep high-radio-frequency spectrum, possibly classifying the targets as compact steep spectrum objects.

preprint2016arXiv

High-velocity OH megamasers in IRAS 20100-4156: Evidence for a Supermassive Black Hole

We report the discovery of new, high-velocity narrow-line components of the OH megamaser in IRAS 20100-4156. Results from the Australian Square Kilometre Array Pathfinder (ASKAP)'s Boolardy Engineering Test Array (BETA) and the Australia Telescope Compact Array (ATCA) provide two independent measurements of the OH megamaser spectrum. We found evidence for OH megamaser clumps at $-$409 and $-$562 km/s (blue-shifted) from the systemic velocity of the galaxy, in addition to the lines previously known. The presence of such high velocities in the molecular emission from IRAS 20100$-$4156 could be explained by a ~50 pc molecular ring enclosing an approximately 3.8 billion solar mass black hole. We also discuss two alternatives, i.e. that the narrow-line masers are dynamically coupled to the wind driven by the active galactic nucleus or they are associated with two separate galactic nuclei. The comparison between the BETA and ATCA spectra provides another scientific verification of ASKAP's BETA. Our data, combined with previous measurements of the source enabled us to study the variability of the source over a twenty-six year period. The flux density of the brightest OH maser components has reduced by more than a factor of two between 1988 and 2015, whereas a secondary narrow-line component has more than doubled in the same time. Plans for high-resolution VLBI follow-up of this source are discussed, as are prospects for discovering new OH megamasers during the ASKAP early science program.

preprint2016arXiv

Radio Galaxy Zoo: discovery of a poor cluster through a giant wide-angle tail radio galaxy

We have discovered a previously unreported poor cluster of galaxies (RGZ-CL J0823.2+0333) through an unusual giant wide-angle tail radio galaxy found in the Radio Galaxy Zoo project. We obtained a spectroscopic redshift of $z=0.0897$ for the E0-type host galaxy, 2MASX J08231289+0333016, leading to M$_r = -22.6$ and a $1.4\,$GHz radio luminosity density of $L_{\rm 1.4} = 5.5\times10^{24}$ W Hz$^{-1}$. These radio and optical luminosities are typical for wide-angle tailed radio galaxies near the borderline between Fanaroff-Riley (FR) classes I and II. The projected largest angular size of $\approx8\,$arcmin corresponds to $800\,$kpc and the full length of the source along the curved jets/trails is $1.1\,$Mpc in projection. X-ray data from the XMM-Newton archive yield an upper limit on the X-ray luminosity of the thermal emission surrounding RGZ J082312.9+033301,at $1.2-2.6\times10^{43}$ erg s$^{-1}$ for assumed intra-cluster medium temperatures of $1.0-5.0\,$keV. Our analysis of the environment surrounding RGZ J082312.9+033301 indicates that RGZ J082312.9+033301 lies within a poor cluster. The observed radio morphology suggests that (a) the host galaxy is moving at a significant velocity with respect to an ambient medium like that of at least a poor cluster, and that (b) the source may have had two ignition events of the active galactic nucleus with $10^7\,$yrs in between. This reinforces the idea that an association between RGZ J082312.9+033301, and the newly discovered poor cluster exists.

preprint2016arXiv

The Australian Square Kilometre Array Pathfinder: Performance of the Boolardy Engineering Test Array

We describe the performance of the Boolardy Engineering Test Array (BETA), the prototype for the Australian Square Kilometre Array Pathfinder telescope ASKAP. BETA is the first aperture synthesis radio telescope to use phased array feed technology, giving it the ability to electronically form up to nine dual-polarization beams. We report the methods developed for forming and measuring the beams, and the adaptations that have been made to the traditional calibration and imaging procedures in order to allow BETA to function as a multi-beam aperture synthesis telescope. We describe the commissioning of the instrument and present details of BETA's performance: sensitivity, beam characteristics, polarimetric properties and image quality. We summarise the astronomical science that it has produced and draw lessons from operating BETA that will be relevant to the commissioning and operation of the final ASKAP telescope.

preprint2016arXiv

Wide-field broadband radio imaging with phased array feeds: a pilot multi-epoch continuum survey with ASKAP-BETA

The Boolardy Engineering Test Array is a 6 x 12 m dish interferometer and the prototype of the Australian Square Kilometre Array Pathfinder (ASKAP), equipped with the first generation of ASKAP's phased array feed (PAF) receivers. These facilitate rapid wide-area imaging via the deployment of simultaneous multiple beams within a 30 square degree field of view. By cycling the array through 12 interleaved pointing positions and using 9 digitally formed beams we effectively mimic a traditional 1 hour x 108 pointing survey, covering 150 square degrees over 711 - 1015 MHz in 12 hours of observing time. Three such observations were executed over the course of a week. We verify the full bandwidth continuum imaging performance and stability of the system via self-consistency checks and comparisons to existing radio data. The combined three epoch image has arcminute resolution and a 1-sigma thermal noise level of 375 micro-Jy per beam, although the effective noise is a factor 3 higher due to residual sidelobe confusion. From this we derive a catalogue of 3,722 discrete radio components, using the 35 percent fractional bandwidth to measure in-band spectral indices for 1,037 of them. A search for transient events reveals one significantly variable source within the survey area. The survey covers approximately two-thirds of the Spitzer South Pole Telescope Deep Field. This pilot project demonstrates the viability and potential of using PAFs to rapidly and accurately survey the sky at radio wavelengths.

preprint2015arXiv

A pilot ASKAP survey of radio transient events in the region around the intermittent pulsar PSR J1107-5907

We use observations from the Boolardy Engineering Test Array (BETA) of the Australian Square Kilometre Array Pathfinder (ASKAP) telescope to search for transient radio sources in the field around the intermittent pulsar PSR J1107-5907. The pulsar is thought to switch between an "off" state in which no emission is detectable, a weak state and a strong state. We ran three independent transient detection pipelines on two-minute snapshot images from a 13 hour BETA observation in order to 1) study the emission from the pulsar, 2) search for other transient emission from elsewhere in the image and 3) to compare the results from the different transient detection pipelines. The pulsar was easily detected as a transient source and, over the course of the observations, it switched into the strong state three times giving a typical timescale between the strong emission states of 3.7 hours. After the first switch it remained in the strong state for almost 40 minutes. The other strong states lasted less than 4 minutes. The second state change was confirmed using observations with the Parkes radio telescope. No other transient events were found and we place constraints on the surface density of such events on these timescales. The high sensitivity Parkes observations enabled us to detect individual bright pulses during the weak state and to study the strong state over a wide observing band. We conclude by showing that future transient surveys with ASKAP will have the potential to probe the intermittent pulsar population.

preprint2015arXiv

ASKAP HI imaging of the galaxy group IC 1459

We present HI imaging of the galaxy group IC 1459 carried out with six antennas of the Australian SKA Pathfinder equipped with phased-array feeds. We detect and resolve HI in eleven galaxies down to a column density of $\sim10^{20}$ cm$^{-2}$ inside a ~6 deg$^2$ field and with a resolution of ~1 arcmin on the sky and ~8 km/s in velocity. We present HI images, velocity fields and integrated spectra of all detections, and highlight the discovery of three HI clouds -- two in the proximity of the galaxy IC 5270 and one close to NGC 7418. Each cloud has an HI mass of $10^9$ M$_\odot$ and accounts for ~15% of the HI associated with its host galaxy. Available images at ultraviolet, optical and infrared wavelengths do not reveal any clear stellar counterpart of any of the clouds, suggesting that they are not gas-rich dwarf neighbours of IC 5270 and NGC 7418. Using Parkes data we find evidence of additional extended, low-column-density HI emission around IC 5270, indicating that the clouds are the tip of the iceberg of a larger system of gas surrounding this galaxy. This result adds to the body of evidence on the presence of intra-group gas within the IC 1459 group. Altogether, the HI found outside galaxies in this group amounts to several times $10^9$ M$_\odot$, at least 10% of the HI contained inside galaxies. This suggests a substantial flow of gas in and out of galaxies during the several billion years of the group's evolution.

preprint2015arXiv

Discovery of HI gas in a young radio galaxy at $z = 0.44$ using the Australian Square Kilometre Array Pathfinder

We report the discovery of a new 21-cm HI absorption system using commissioning data from the Boolardy Engineering Test Array of the Australian Square Kilometre Array Pathfinder (ASKAP). Using the 711.5 - 1015.5 MHz band of ASKAP we were able to conduct a blind search for the 21-cm line in a continuous redshift range between $z = 0.4$ and 1.0, which has, until now, remained largely unexplored. The absorption line is detected at $z = 0.44$ towards the GHz-peaked spectrum radio source PKS B1740$-$517 and demonstrates ASKAP's excellent capability for performing a future wide-field survey for HI absorption at these redshifts. Optical spectroscopy and imaging using the Gemini-South telescope indicates that the HI gas is intrinsic to the host galaxy of the radio source. The narrow OIII emission lines show clear double-peaked structure, indicating either large-scale outflow or rotation of the ionized gas. Archival data from the \emph{XMM-Newton} satellite exhibit an absorbed X-ray spectrum that is consistent with a high column density obscuring medium around the active galactic nucleus. The HI absorption profile is complex, with four distinct components ranging in width from 5 to 300 km s$^{-1}$ and fractional depths from 0.2 to 20 per cent. In addition to systemic HI gas, in a circumnuclear disc or ring structure aligned with the radio jet, we find evidence for a possible broad outflow of neutral gas moving at a radial velocity of $v \sim 300$ km s$^{-1}$. We infer that the expanding young radio source ($t_{\rm age} \approx 2500$ yr) is cocooned within a dense medium and may be driving circumnuclear neutral gas in an outflow of $\sim$ 1 $\mathrm{M}_{\odot}$ yr$^{-1}$.

preprint2015arXiv

Radio Galaxy Zoo: host galaxies and radio morphologies derived from visual inspection

We present results from the first twelve months of operation of Radio Galaxy Zoo, which upon completion will enable visual inspection of over 170,000 radio sources to determine the host galaxy of the radio emission and the radio morphology. Radio Galaxy Zoo uses $1.4\,$GHz radio images from both the Faint Images of the Radio Sky at Twenty Centimeters (FIRST) and the Australia Telescope Large Area Survey (ATLAS) in combination with mid-infrared images at $3.4\,μ$m from the {\it Wide-field Infrared Survey Explorer} (WISE) and at $3.6\,μ$m from the {\it Spitzer Space Telescope}. We present the early analysis of the WISE mid-infrared colours of the host galaxies. For images in which there is $>\,75\%$ consensus among the Radio Galaxy Zoo cross-identifications, the project participants are as effective as the science experts at identifying the host galaxies. The majority of the identified host galaxies reside in the mid-infrared colour space dominated by elliptical galaxies, quasi-stellar objects (QSOs), and luminous infrared radio galaxies (LIRGs). We also find a distinct population of Radio Galaxy Zoo host galaxies residing in a redder mid-infrared colour space consisting of star-forming galaxies and/or dust-enhanced non star-forming galaxies consistent with a scenario of merger-driven active galactic nuclei (AGN) formation. The completion of the full Radio Galaxy Zoo project will measure the relative populations of these hosts as a function of radio morphology and power while providing an avenue for the identification of rare and extreme radio structures. Currently, we are investigating candidates for radio galaxies with extreme morphologies, such as giant radio galaxies, late-type host galaxies with extended radio emission, and hybrid morphology radio sources.

preprint2015arXiv

Strongly lensed neutral hydrogen emission: detection predictions with current and future radio interferometers

Strong gravitational lensing provides some of the deepest views of the Universe, enabling studies of high-redshift galaxies only possible with next-generation facilities without the lensing phenomenon. To date, 21 cm radio emission from neutral hydrogen has only been detected directly out to z~0.2, limited by the sensitivity and instantaneous bandwidth of current radio telescopes. We discuss how current and future radio interferometers such as the Square Kilometre Array (SKA) will detect lensed HI emission in individual galaxies at high redshift. Our calculations rely on a semi-analytic galaxy simulation with realistic HI disks (by size, density profile and rotation), in a cosmological context, combined with general relativistic ray tracing. Wide-field, blind HI surveys with the SKA are predicted to be efficient at discovering lensed HI systems, increasingly so at z > 2. This will be enabled by the combination of the magnification boosts, the steepness of the HI luminosity function at the high-mass end, and the fact that the HI spectral line is relatively isolated in frequency. These surveys will simultaneously provide a new technique for foreground lens selection and yield the highest redshift HI emission detections. More near term (and existing) cm-wave facilities will push the high redshift HI envelope through targeted surveys of known lenses.

preprint2015arXiv

Tracing the neutral gas environments of young radio AGN with ASKAP

At present neutral atomic hydrogen (HI) gas in galaxies at redshifts above $z \sim 0.3$ (the extent of 21-cm emission surveys in individual galaxies) and below $z \sim 1.7$ (where the Lyman-$α$ line is not observable with ground-based telescopes) has remained largely unexplored. The advent of precursor telescopes to the Square Kilometre Array will allow us to conduct the first systematic radio-selected 21-cm absorption surveys for HI over these redshifts. While HI absorption is a tracer of the reservoir of cold neutral gas in galaxies available for star formation, it can also be used to reveal the extreme kinematics associated with jet-driven neutral outflows in radio-loud active galactic nuclei. Using the six-antenna Boolardy Engineering Test Array of the Australian Square Kilometre Array Pathfinder, we have demonstrated that in a single frequency tuning we can detect HI absorption over a broad range of redshifts between $z = 0.4$ and $1.0$. As part of our early science and commissioning program, we are now carrying out a search for absorption towards a sample of the brightest GPS and CSS sources in the southern sky. These intrinsically compact sources present us with an opportunity to study the circumunuclear region of recently re-started radio galaxies, in some cases showing direct evidence of mechanical feedback through jet-driven outflows. With the sensitivity of the full ASKAP array we will be able to study the kinematics of atomic gas in a few thousand radio galaxies, testing models of radio jet feedback well beyond the nearby Universe

preprint2014arXiv

A close-pair binary in a distant triple supermassive black-hole system

Galaxies are believed to evolve through merging, which should lead to multiple supermassive black holes in some. There are four known triple black hole systems, with the closest pair being 2.4 kiloparsecs apart (the third component is more distant at 3 kiloparsecs), which is far from the gravitational sphere of influence of a black hole with mass $\sim$10$^9$ M$_\odot$ (about 100 parsecs). Previous searches for compact black hole systems concluded that they were rare, with the tightest binary system having a separation of 7 parsecs. Here we report observations of a triple black hole system at redshift z=0.39, with the closest pair separated by $\sim$140 parsecs. The presence of the tight pair is imprinted onto the properties of the large-scale radio jets, as a rotationally-symmetric helical modulation, which provides a useful way to search for other tight pairs without needing extremely high resolution observations. As we found this tight pair after searching only six galaxies, we conclude that tight pairs are more common than hitherto believed, which is an important observational constraint for low-frequency gravitational wave experiments.

preprint2014arXiv

Enabling the next generation of cm-wavelength studies of high-redshift molecular gas with the SKA

The Square Kilometre Array will be a revolutionary instrument for the study of gas in the distant Universe. SKA1 will have sufficient sensitivity to detect and image atomic 21 cm HI in individual galaxies at significant cosmological distances, complementing ongoing ALMA imaging of redshifted high-J CO line emission and far-infrared interstellar medium lines such as [CII] 157.7 um. At frequencies below ~50 GHz, observations of redshifted emission from low-J transitions of CO, HCN, HCO+, HNC, H2O and CS provide insight into the kinematics and mass budget of the cold, dense star-forming gas in galaxies. In advance of ALMA band 1 deployment (35 to 52 GHz), the most sensitive facility for high-redshift studies of molecular gas operating below 50~GHz is the Karl G. Jansky Very Large Array (VLA). Here, we present an overview of the role that the SKA could play in molecular emission line studies during SKA1 and SKA2, with an emphasis on studies of the dense gas tracers directly probing regions of active star-formation.

preprint2013arXiv

Prediction of the Cosmic Evolution of the CO-Luminosity Functions

We predict the emission line luminosity functions (LFs) of the first 10 rotational transitions of CO in galaxies at redshift z=0 to z=10. This prediction relies on a recently presented simulation of the molecular cold gas content in ~3e7 evolving galaxies based on the Millennium Simulation. We combine this simulation with a model for the conversion between molecular mass and CO-line intensities, which incorporates the following mechanisms: (i) molecular gas is heated by the CMB, starbursts (SBs), and active galactic nuclei (AGNs); (ii) molecular clouds in dense or inclined galaxies can overlap; (iii) compact gas can attain a smooth distribution in the densest part of disks; (iv) CO-luminosities scale with metallicity changes between galaxies; (v) CO-luminosities are always detected against the CMB. We analyze the relative importance of these effects and predict the cosmic evolution of the CO-LFs. The most notable conclusion is that the detection of regular galaxies (i.e. no AGN, no massive SB) at high z>7 in CO-emission will be dramatically hindered by the weak contrast against the CMB, in contradiction to earlier claims that CMB-heating will ease the detection of high-redshift CO. The full simulation of extragalactic CO-lines and the predicted CO-LFs at any redshift can be accessed online, prior registration required} and they should be useful for the modeling of CO-line surveys with future telescopes, such as ALMA, the LMT, or the SKA.

preprint2013arXiv

The preferentially magnified active nucleus in IRAS F10214+4724 - I. Lens model and spatially resolved radio emission

This is the first paper in a series that present a multi-wavelength analysis of the archetype Ultra-Luminous InfraRed Galaxy (ULIRG) IRAS FSC10214+4724, a gravitationally lensed, starburst/AGN at z=2.3. Here we present a new lens model and spatially-resolved radio data, as well as a deep HST F160W map. The lens modelling employs a Bayesian Markov Chain Monte Carlo algorithm with extended-source, forward ray-tracing. Using these high resolution HST, MERLIN and VLA maps, the algorithm allows us to constrain the level of distortion to the continuum spectral energy distribution resulting from emission components with differing magnification factors, due to their size and proximity to the caustic. Our lens model finds the narrow line region (NLR), and by proxy the active nucleus, is preferentially magnified. This supports previous claims that preferential magnification could mask the expected polycyclic aromatic hydrocarbon spectral features in the Spitzer mid-infrared spectrum which roughly trace the star-forming regions. Furthermore, we show the arc-to-counter-image flux ratio is not a good estimate of the magnification in this system, despite its common use in the IRAS FSC10214+4724 literature. Our lens modelling suggests magnifications of μ~ 15-20+-2 for the HST F814W, MERLIN 1.7 GHz and VLA 8 GHz maps, significantly lower than the canonical values of μ= 50-100 often used for this system. Systematic errors such as the dark matter density slope and co-location of stellar and dark matter centroids dominate the uncertainties in the lens model at the 40 percent level.

preprint2013arXiv

The preferentially magnified active nucleus in IRAS F10214+4724 - II. Spatially resolved cold molecular gas

We present JVLA observations of the cold (CO (1-0)) molecular gas in IRAS F10214+4724, a lensed ULIRG at z=2.3 with an obscured active nucleus. The galaxy is spatially and spectrally well-resolved in the CO (1-0) emission line. A CO (1-0) counter-image is detected at the 3-sigma level. Five of the 42 km/s channels (with >5-sigma detections) are mapped back into the source plane and their total magnification posterior PDFs sampled. This reveals a roughly linear arrangement, tentatively a rotating disk. We derive a molecular gas mass of M_gas = 1.2 +- 0.2 x 10^10 M_sun, assuming a ULIRG L_{CO}-to-M_{gas} conversion ratio of α= 0.8 M_sun / (K km/s pc^2) that agrees well with the derived range of α= 0.3 - 1.3 for separate dynamical mass estimates at assumed inclinations of i = 90 - 30 degrees. Based on the AGN and CO (1-0) peak emission positions and the lens model, we predict a distortion of the CO Spectral Line Energy Distribution (SLED) where higher order J lines that may be partially excited by AGN heating will be preferentially lensed owing to their smaller solid angles and closer proximity to the AGN and therefore the cusp of the caustic. Comparison with other lensing inversion results shows that the narrow line region and AGN radio core in IRAS F10214+4724 are preferentially lensed by a factor >~ 3 and 11 respectively, relative to the molecular gas emission. This distorts the global continuum emission Spectral Energy Distribution (SED) and suggests caution in unsophisticated uses of IRAS F10214+4724 as an archetype high-redshift ULIRG. We explore two Large Velocity Gradient (LVG) models, incorporating spatial CO (1-0) and (3-2) information and present tentative evidence for an extended, low excitation cold gas component that implies that the total molecular gas mass in IRAS F10214+4724 is a factor >~2 greater than that calculated using spatially unresolved CO observations.

preprint2013arXiv

The preferentially magnified active nucleus in IRAS F10214+4724 - III. VLBI observations of the radio core

We report 1.7 GHz Very Long Baseline Interferometry (VLBI) observations of IRAS F10214+4724, a lensed z=2.3 obscured quasar with prodigious star formation. We detect what we argue to be the obscured active nucleus with an effective angular resolution of < 50 pc at z = 2.3 . The S_{1.7} = 210 micro-Jy (9-σ) detection of this unresolved source is located within the HST rest-frame ultraviolet/optical arc, however, >~100 mas northward of the arc centre of curvature. This leads to a source plane inversion that places the European VLBI Network detection to within milli-arcseconds of the modelled cusp caustic, resulting in a very large magnification (μ~70), over an order of magnitude larger than the CO (1-0) derived magnification of a spatially resolved JVLA map, using the same lens model. We estimate the quasar bolometric luminosity from a number of independent techniques and with our X-ray modelling find evidence that the AGN may be close to Compton-thick, with an intrinsic bolometric luminosity log(L_{bol,QSO} / L_sun) = 11.34 +- 0.27 dex. We make the first black hole mass estimate of IRAS F10214+4724 and find log(M_{BH}/M_sun) = 8.36 +- 0.56 which suggests a low black hole accretion rate (λ= \dot{M} / \dot{M}_{Edd} ~ 3\pm^7_2 percent). We find evidence for a M_{BH}/M_{spheroid} ratio that is 1-2 orders of magnitude larger than that of submillimetre galaxies (SMGs) at z~2. At face value, this suggests IRAS F10214+4724 has undergone a different evolutionary path compared to SMGs at the same epoch. A primary result of this work is the demonstration that emission regions of differing size and position can undergo significantly different magnification boosts (> 1 dex) and therefore distort our view of high-redshift, gravitationally lensed galaxies.

preprint2012arXiv

Sub-millimetre source identifications and the micro-Jansky source population at 8.4 GHz in the William Herschel Deep Field

[Abridged] Sub-mm observations of the William Herschel Deep Field using LABOCA revealed possible counterparts for 2 X-ray absorbed QSOs. The aim here is to exploit EVLA imaging at 8.4 GHz to establish the QSOs as radio/sub-mm sources. The challenge in reducing the EVLA data was the presence of a strong 4C source in the field. A new calibration algorithm was applied to the data to subtract it. The resulting thermal noise limited radio map covers the 16'x16' Extended WHDF. It contains 41 sources above a 4-sigma limit, 17 of which have primary beam corrected flux. The radio observations show that the absorbed AGN with LABOCA detections are coincident with radio sources, confirming the tendency for X-ray absorbed AGN to be sub-mm bright. These sources show strong ultraviolet excess (UVX) suggesting the nuclear sightline is gas- but not dust-absorbed. Of the 3 remaining LABOCA sources within the ~5' half-power beam width, 1 is identified with a faint nuclear X-ray/radio source in a nearby galaxy, 1 with a faint radio source and 1 is unidentified in any other band. More generally, differential radio source counts are in good agreement with previous observations, showing at S<50 micro-Jy a significant excess over a pure AGN model. In the full area, of 10 sources fainter than this limit, 6 have optical counterparts of which 3 are UVX (i.e. likely QSOs) including the 2 absorbed quasar LABOCA sources. The other faint radio counterparts are not UVX but are only slightly less blue and likely to be star-forming/merging galaxies, predominantly at lower luminosities and redshifts. The 4 faint, optically unidentified radio sources may be either dust obscured QSOs or galaxies. These high-z obscured AGN and lower-z star-forming populations are thus the main candidates to explain the observed excess in faint source counts and hence the excess radio background found previously by the ARCADE2 experiment.

preprint2012arXiv

The accretion disc, jets and environment of the intermediate mass black hole candidate ESO 243-49 HLX-1

The Ultra Luminous X-ray (ULX) source HLX-1 in the galaxy ESO 243-49 has an observed maximum unabsorbed X-ray luminosity of 1.3e42 erg/s (0.2-10.0 keV). From the conservative assumption that this value exceeds the Eddington limit by at most a factor of 10, the minimum mass is then 500 solar masses. The X-ray luminosity varies by a factor of 40 with an apparent recurrence timescale of approximately one year. This X-ray variability is associated with spectral state transitions similar to those seen in black hole X-ray binaries. Here we discuss our recent modelling of all the X-ray data for HLX-1 and show that it supports the idea that this ULX is powered by sub- and near Eddington accretion onto an intermediate mass black hole. We also present evidence for transient radio emission which is consistent with a discrete jet ejection event as well as comment on the nature of the environment around HLX-1 in light of recent Hubble Space Telescope photometry.

preprint2011arXiv

Expanding e-MERLIN with the Goonhilly Earth Station

A consortium of universities has recently been formed with the goal of using the decommissioned telecommunications infrastructure at the Goonhilly Earth Station in Cornwall, UK, for astronomical purposes. One particular goal is the introduction of one or more of the ~30-metre parabolic antennas into the existing e-MERLIN radio interferometer. This article introduces this scheme and presents some simulations which quantify the improvements that would be brought to the e-MERLIN system. These include an approximate doubling of the spatial resolution of the array, an increase in its N-S extent with strong implications for imaging the most well-studied equatorial fields, accessible to ESO facilities including ALMA. It also increases the overlap between the e-MERLIN array and the European VLBI Network. We also discuss briefly some niche science areas in which an e-MERLIN array which included a receptor at Goonhilly would be potentially world-leading, in addition to enhancing the existing potential of e-MERLIN in its role as a Square Kilometer Array pathfinder instrument.

preprint2011arXiv

Goonhilly: a new site for e-MERLIN and the EVN

The benefits for the e-MERLIN and EVN arrays of using antennae at the satellite communication station at Goonhilly in Cornwall are discussed. The location of this site - new to astronomy - will provide an almost equal distribution of long baselines in the east-west- and north-south directions, and opens up the possibility to get significantly improved observations of equatorial fields with e-MERLIN. These additional baselines will improve the sensitivity on a set of critical spatial scales and will increase the angular resolution of e-MERLIN by a factor of two. e-MERLIN observations, including many allocated under the e-MERLIN Legacy programme, will benefit from the enhanced angular resolution and imaging capability especially for sources close to or below the celestial equator (where ESO facilities such as ALMA will operate) of including the Goonhilly telescopes. Furthermore, the baselines formed between Goonhilly and the existing stations will close the gap between the baselines of e-MERLIN and those of the European VLBI Network (EVN) and therefore enhance the legacy value of e-MERLIN datasets.

preprint2011arXiv

MESMER: MeerKAT Search for Molecules in the Epoch of Reionization

[Abridged] Observations of molecular gas at all redshifts are critical for measuring the cosmic evolution in molecular gas density and understanding the star-formation history of the Universe. The 12CO molecule (J=1-0 transition = 115.27 GHz) is the best proxy for extragalactic H2, which is the gas reservoir from which star formation occurs, and has been detected out to z~6. Typically, redshifted high-J lines are observed at mm-wavelengths, the most commonly targeted systems exhibiting high SFRs (e.g. submm galaxies), and far-IR-bright QSOs. While the most luminous objects are the most readily observed, detections of more typical galaxies with modest SFRs are essential for completing the picture. ALMA will be revolutionary in terms of increasing the detection rate and pushing the sensitivity limit down to include such galaxies, however the limited FoV when observing at such high frequencies makes it difficult to use ALMA for studies of the large-scale structure traced out by molecular gas in galaxies. This article introduces a strategy for a systematic search for molecular gas during the EoR (z~7 and above), capitalizing on the fact that the J=1-0 transition of 12CO enters the upper bands of cm-wave instruments at high-z. The FoV advantage gained by observing at such frequencies, coupled with modern broadband correlators allows significant cosmological volumes to be probed on reasonable timescales. In this article we present an overview of our future observing programme which has been awarded 6,500 hours as one of the Large Survey Projects for MeerKAT, the forthcoming South African SKA pathfinder instrument. Its large FoV and correlator bandwidth, and high-sensitivity provide unprecedented survey speed for such work. An existing astrophysical simulation is coupled with instrumental considerations to demonstrate the feasibility of such observations and predict detection rates.

preprint2011arXiv

The ALMA Observation Support Tool

The ALMA Observation Support Tool (OST) is an ALMA simulator which is interacted with solely via a standard web browser. It is aimed at users who may or may not be experts in interferometry, or those that do not wish to familarise themselves with the simulation components of a data reduction package. It has been designed to offer full imaging simulation capability for an arbitrary ALMA observation while maintaining the accessibility of other online tools such as the ALMA Sensitivity Calculator. Simulation jobs are defined by selecting and entering options on a standard web form. The user can specify the standard parameters that would need to be considered for an ALMA observation (e.g. pointing direction, frequency set up, duration), and there is also the option to upload arbitrary sky models in FITS format. Once submitted, jobs are sequentially processed by a remote server running a CASA-based back-end system. The user is notified by email when the job is complete, and directed to a standard web page which contains the results of the simulation and a range of downloadable data products. The system is currently hosted by the UK ALMA Regional Centre, and can be accessed by directing a web browser to http://almaost.jb.man.ac.uk.

preprint2010arXiv

Further Observations of the Intermediate Mass Black Hole Candidate ESO 243-49 HLX-1

The brightest Ultra-Luminous X-ray source HLX-1 in the galaxy ESO 243-49 currently provides strong evidence for the existence of intermediate mass black holes. Here we present the latest multi-wavelength results on this intriguing source in X-ray, UV and radio bands. We have refined the X-ray position to sub-arcsecond accuracy. We also report the detection of UV emission that could indicate ongoing star formation in the region around HLX-1. The lack of detectable radio emission at the X-ray position strengthens the argument against a background AGN.

preprint2010arXiv

Herschel-ATLAS: far-infrared properties of radio-selected galaxies

We use the Herschel-ATLAS science demonstration data to investigate the star-formation properties of radio-selected galaxies in the GAMA-9h field as a function of radio luminosity and redshift. Radio selection at the lowest radio luminosities, as expected, selects mostly starburst galaxies. At higher radio luminosities, where the population is dominated by AGN, we find that some individual objects are associated with high far-infrared luminosities. However, the far-infrared properties of the radio-loud population are statistically indistinguishable from those of a comparison population of radio-quiet galaxies matched in redshift and K-band absolute magnitude. There is thus no evidence that the host galaxies of these largely low-luminosity (Fanaroff-Riley class I), and presumably low-excitation, AGN, as a population, have particularly unusual star-formation histories. Models in which the AGN activity in higher-luminosity, high-excitation radio galaxies is triggered by major mergers would predict a luminosity-dependent effect that is not seen in our data (which only span a limited range in radio luminosity) but which may well be detectable with the full Herschel-ATLAS dataset.

preprint2010arXiv

SKA HI end2end simulation

The current status of the HI simulation efforts is presented, in which a self consistent simulation path is described and basic equations to calculate array sensitivities are given. There is a summary of the SKA Design Study (SKADS) sky simulation and a method for implementing it into the array simulator is presented. A short overview of HI sensitivity requirements is discussed and expected results for a simulated HI survey are presented.

preprint2009arXiv

Multiwavelength study of Cygnus A III. Evidence for relic lobe plasma

We study the particle energy distribution in the cocoon surrounding Cygnus A, using radio images between 151 MHz and 15 GHz and a 200 ks Chandra ACIS-I image. We show that the excess low frequency emission in the the lobe further from Earth cannot be explained by absorption or excess adiabatic expansion of the lobe or a combination of both. We show that this excess emission is consistent with emission from a relic counterlobe and a relic counterjet that are being re-energized by compression from the current lobe. We detect hints of a relic hotspot at the end of the relic X-ray jet in the more distant lobe. We do not detect relic emission in the lobe nearer to Earth as expected from light travel-time effects assuming intrinsic symmetry. We determine that the duration of the previous jet activity phase was slightly less than that of the current jet-active phase. Further, we explain some features observed at 5 and 15 GHz as due to the presence of a relic jet.