Source author record

H. J. A. Röttgering

H. J. A. Röttgering 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

70works
7topics
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

70 published item(s)

preprint2022arXiv

Accretion mode versus radio morphology in the LOFAR Deep Fields

Radio-loud active galaxies have two accretion modes [radiatively inefficient (RI) and radiatively efficient (RE)], with distinct optical and infrared signatures, and two jet dynamical behaviours, which in arcsec- to arcmin-resolution radio surveys manifest primarily as centre- or edge-brightened structures [Fanaroff-Riley (FR) class I and II]. The nature of the relationship between accretion mode and radio morphology (FR class) has been the subject of long debate. We present a comprehensive investigation of this relationship for a sample of 286 well-resolved radio galaxies in the LOFAR Two-metre Sky Survey Deep Fields (LoTSS-Deep) first data release, for which robust morphological and accretion mode classifications have been made. We find that two-thirds of luminous FRII radio galaxies are RI, and identify no significant differences in the visual appearance or source dynamic range (peak/mean surface brightness) of the RI and RE FRIIs, demonstrating that both RI and RE systems can produce FRII structures. We also find a significant population of low-luminosity FRIIs (predominantly RI), supporting our earlier conclusion that FRII radio structures can be produced at all radio luminosities. We demonstrate that in the luminosity range where both morphologies are present, the probability of producing FRI or FRII radio morphology is directly linked to stellar mass, while across all morphologies and luminosities, RE accretion occurs in systems with high specific star formation rate, presumably because this traces fuel availability. In summary, the relationship between accretion mode and radio morphology is very indirect, with host-galaxy environment controlling these two key parameters in different ways.

preprint2022arXiv

Cosmic evolution of low-excitation radio galaxies in the LOFAR Two-meter Sky Survey Deep Fields

Feedback from low-excitation radio galaxies (LERGs) plays a key role in the lifecycle of massive galaxies in the local Universe; their evolution, and the impact of these active galactic nuclei on early galaxy evolution, however, remain poorly understood. We use a sample of 10481 LERGs from the first data release of the LOFAR Two-meter Sky Survey Deep Fields, covering $\sim$ 25 deg$^2$, to present the first measurement of the evolution of the radio luminosity function (LF) of LERGs out to $z\sim2.5$; this shows relatively mild evolution. We split the LERGs into those hosted by quiescent and star-forming galaxies, finding a new dominant population of LERGs hosted by star-forming galaxies at high redshifts. The incidence of LERGs in quiescent galaxies shows a steep dependence on stellar-mass out to $z \sim1.5$, consistent with local Universe measurements of accretion occurring from cooling of hot gas haloes. The quiescent-LERGs dominate the LFs at $z<1$, showing a strong decline in space density with redshift, tracing that of the available host galaxies, while there is an increase in the characteristic luminosity. The star-forming LERG LF increases with redshift, such that this population dominates the space densities at most radio-luminosities by $z \sim 1$. The incidence of LERGs in star-forming galaxies shows a much weaker stellar-mass dependence, and increases with redshift, suggesting a different fuelling mechanism compared to their quiescent counterparts, potentially associated with the cold gas supply present in the star-forming galaxies.

preprint2022arXiv

Deep low-frequency radio observations of Abell 2256 I: The filamentary radio relic

We present deep and high fidelity images of the merging galaxy cluster Abell 2256 at low frequencies, using the upgraded Giant Metrewave Radio Telescope (uGMRT) and LOw-Frequency ARray (LOFAR). This cluster hosts one of the most prominent known relics, with a remarkably spectacular network of filamentary substructures. The new uGMRT (300-850 MHz) and LOFAR (120-169 MHz) observations, combined with the archival Karl G. Jansky Very Large Array (VLA; 1-4 GHz) data, allowed us to carry out the first spatially resolved spectral analysis of the exceptional relic emission down to 6 arcsec resolution over a broad range of frequencies. Our new sensitive radio images confirm the presence of complex filaments of magnetized relativistic plasma also at low frequencies. We find that the integrated spectrum of the relic is consistent with a single power law, without any sign of spectral steepening, at least below 3 GHz. Unlike previous claims, the relic shows an integrated spectral index of $-1.07\pm0.02$ between 144 MHz and 3 GHz, which is consistent with the (quasi)stationary shock approximation. The spatially resolved spectral analysis suggests that the relic surface very likely traces the complex shock front, with a broad distribution of Mach numbers propagating through a turbulent and dynamically active intracluster medium. Our results show that the northern part of the relic is seen edge-on and the southern part close to face-on. We suggest that the complex filaments are regions where higher Mach numbers dominate the (re-)acceleration of electrons that are responsible for the observed radio emission.

preprint2022arXiv

Deep sub-arcsecond widefield imaging of the Lockman Hole field at 144 MHz

High quality low-frequency radio surveys have the promise of advancing our understanding of many important topics in astrophysics, including the life cycle of active galactic nuclei (AGN), particle acceleration processes in jets, the history of star formation, and exoplanet magnetospheres. Currently leading low-frequency surveys reach an angular resolution of a few arcseconds. However, this resolution is not yet sufficient to study the more compact and distant sources in detail. Sub-arcsecond resolution is therefore the next milestone in advancing these fields. The biggest challenge at low radio frequencies is the ionosphere. If not adequately corrected for, ionospheric seeing blurs the images to arcsecond or even arcminute scales. Additionally, the required image size to map the degree-scale field of view of low-frequency radio telescopes at this resolution is far greater than what typical soft- and hardware is currently capable of handling. Here we present for the first time (to the best of our knowledge) widefield sub-arcsecond imaging at low radio frequencies. We derive ionospheric corrections in a few dozen individual directions and apply those during imaging efficiently using a recently developed imaging algorithm (arXiv:1407.1943, arXiv:1909.07226). We demonstrate our method by applying it to an eight hour observation of the International LOw Frequency ARray (LOFAR) Telescope (ILT) (arXiv:1305.3550). Doing so we have made a sensitive $7.4\ \mathrm{deg}^2$ $144\ \mathrm{MHz}$ map at a resolution of $0.3''$ reaching $25\ μ\mathrm{Jy\ beam}^{-1}$ near the phase centre. The estimated $250,000$ core hours used to produce this image, fit comfortably in the budget of available computing facilities. This result will enable future mapping of the entire northern low-frequency sky at sub-arcsecond resolution.

preprint2022arXiv

Diffuse radio emission from non-Planck galaxy clusters in the LoTSS-DR2 fields

The presence of large-scale magnetic fields and ultra-relativistic electrons in the intra-cluster medium (ICM) is confirmed through the detection of diffuse radio synchrotron sources, so-called radio halos and relics. Due to their steep-spectrum nature, these sources are rarely detected at frequencies above a few GHz, especially in low-mass systems. The aim of this study is to discover and characterise diffuse radio sources in low-mass galaxy clusters in order to understand their origin and their scaling with host cluster properties. We searched for cluster-scale radio emission from low-mass galaxy clusters in the Low Frequency Array (LOFAR) Two-metre Sky Survey - Data Release 2 (LoTSS-DR2) fields. We made use of existing optical (Abell, DESI, WHL) and X-ray (comPRASS, MCXC) catalogues. The LoTSS-DR2 data were processed further to improve the quality of the images that are used to detect and characterize diffuse sources. We have detected diffuse radio emission in 28 galaxy clusters. The number of confirmed (candidates) halos and relics are six (seven) and 10 (three), respectively. Among these, 11 halos and 10 relics, including candidates, are newly discovered by LOFAR. Beside these, five diffuse sources are detected in tailed radio galaxies and are probably associated with mergers during the formation of the host clusters. We are unable to classify other 13 diffuse sources. We compare our newly detected, diffuse sources to known sources by placing them on the scaling relation between the radio power and the mass of the host clusters.

preprint2022arXiv

Filamentary structures of ionized gas in Cygnus X

Ionized gas probes the influence of massive stars on their environment. The Cygnus X region (d~1.5 kpc) is one of the most massive star forming complexes in our Galaxy, in which the Cyg OB2 association (age of 3-5 Myr and stellar mass $2 \times 10^{4}$ M$_{\odot}$) has a dominant influence. We observe the Cygnus X region at 148 MHz using the Low Frequency Array (LOFAR) and take into account short-spacing information during image deconvolution. Together with data from the Canadian Galactic Plane Survey, we investigate the morphology, distribution, and physical conditions of low-density ionized gas in a $4^{\circ} \times 4^{\circ}$ (100 pc $\times$ 100 pc) region at a resolution of 2' (0.9 pc). The Galactic radio emission in the region analyzed is almost entirely thermal (free-free) at 148 MHz, with emission measures of $10^3 < EM~{\rm[pc~cm^{-6}]} < 10^6$. As filamentary structure is a prominent feature of the emission, we use DisPerSE and FilChap to identify filamentary ridges and characterize their radial ($EM$) profiles. The distribution of radial profiles has a characteristic width of 4.3 pc and a power-law distribution ($β= -1.8 \pm 0.1$) in peak $EM$ down to our completeness limit of 4200 pc cm$^{-6}$. The electron densities of the filamentary structure range from $10 < n_e~{\rm[cm^{-3}]} < 400$ with a median value of 35 cm$^{-3}$, remarkably similar to [N II] surveys of ionized gas. Cyg OB2 may ionize at most two-thirds of the total ionized gas and the ionized gas in filaments. More than half of the filamentary structures are likely photoevaporating surfaces flowing into a surrounding diffuse (~5 cm$^{-3}$) medium. However, this is likely not the case for all ionized gas ridges. A characteristic width in the distribution of ionized gas points to the stellar winds of Cyg OB2 creating a fraction of the ionized filaments through swept-up ionized gas or dissipated turbulence.

preprint2022arXiv

Nearby galaxies in LoTSS-DR2: insights into the non-linearity of the radio-SFR relation

Context. Cosmic rays and magnetic fields are key ingredients in galaxy evolution, regulating both stellar feedback and star formation. Their properties can be studied with low-frequency radio continuum observations, free from thermal contamination. Aims. We define a sample of 76 nearby (< 30 Mpc) galaxies, with rich ancillary data in the radio continuum and infrared from the CHANG-ES and KINGFISH surveys, which will be observed with the LOFAR Two-metre Sky Survey (LoTSS) at 144 MHz. Methods. We present maps for 45 of them as part of the LoTSS data release 2 (LoTSS-DR2), where we measure integrated flux densities and study integrated and spatially resolved radio spectral indices. We investigate the radio-SFR relation, using star-formation rates (SFR) from total infrared and H $α$ + 24-$μ$m emission. Results. The radio-SFR relation at 144 MHz is clearly super-linear with $L_{144} \propto SFR^{1.4-1.5}$. The mean integrated radio spectral index between 144 and $\approx$1400 MHz is $\langle α\rangle = -0.56 \pm 0.14$, in agreement with the injection spectral index for cosmic ray electrons (CRE). However, the radio spectral index maps show a variation of spectral indices with flatter spectra associated with star-forming regions and steeper spectra in galaxy outskirts and, in particular, in extra-planar regions. We found that galaxies with high star-formation rates (SFR) have steeper radio spectra; we find similar correlations with galaxy size, mass, and rotation speed. Conclusions. Galaxies that are larger and more massive are better electron calorimeters, meaning that the CRE lose a higher fraction of their energy within the galaxies. This explains the super-linear radio-SFR relation, with more massive, star-forming galaxies being radio bright. We propose a semi-calorimetric radio-SFR relation, which employs the galaxy mass as a proxy for the calorimetric efficiency.

preprint2022arXiv

Revealing new high redshift quasar populations through Gaussian mixture model selection

We present a novel method to identify candidate high redshift quasars (HzQs; ($z\gtrsim5.5$), which are unique probes of supermassive black hole growth in the early Universe, from large area optical/infrared photometric surveys. Using Gaussian Mixture Models to construct likelihoods and incorporate informed priors based on population statistics, our method uses a Bayesian framework to assign posterior probabilities that differentiate between HzQs and contaminating sources. We additionally include deep radio data to obtain informed priors. Using existing HzQ data in the literature, we set a posterior threshold that accepts ${\sim}90\%$ of known HzQs while rejecting $>99\%$ of contaminants such as dwarf stars or lower redshift galaxies. Running the probability selection on test samples of simulated HzQs and contaminants, we find that the efficacy of the probability method is higher than traditional colour cuts, decreasing the fraction of accepted contaminants by 86% while retaining a similar fraction of HzQs. As a test, we apply our method to the Pan-STARRS Data Release 1 (PS1) source catalogue within the HETDEX Spring field area on the sky, covering 400 sq. deg. and coinciding with deep radio data from the LOFAR Two-metre Sky Survey Data Release 1 (LoTSS DR1). From an initial sample of ${\sim}5\times10^5$ sources in PS1, our selection shortlists 251 candidate HzQs, which are further reduced to 63 after visual inspection. Shallow spectroscopic follow-up of 13 high probability HzQs resulted in the confirmation of a previously undiscovered quasar at $z=5.66$ with photometric colours $i-z = 1.4$, lying outside the typically probed regions when selecting HzQs based on colours. This discovery demonstrates the efficacy of our probabilistic HzQ selection method in selecting more complete HzQ samples, which holds promise when employed on large existing and upcoming photometric data sets.

preprint2022arXiv

The galaxy group NGC 507: newly detected AGN remnant plasma transported by sloshing

Jets from active galactic nuclei (AGN) are known to recurrently enrich their surrounding medium with mildly-relativistic particles and magnetic fields. Here, we present a detailed multi-frequency analysis of the nearby (z=0.01646) galaxy group NGC 507. In particular, we present new high-sensitivity and high spatial resolution radio images in the frequency range 144-675 MHz obtained using LOFAR and uGMRT observations. These reveal the presence of previously undetected diffuse radio emission with complex, filamentary morphology, likely related to a previous outburst of the central galaxy. Based on spectral ageing considerations, we derived that the plasma was first injected by the AGN 240-380 Myr ago and is now cooling. Our analysis of deep archival XMM-Newton data confirms that the system is dynamically disturbed, as previously suggested. We detect two discontinuities in the X-ray surface brightness distribution (in East and South direction) tracing a spiral pattern, which we interpret as cold fronts produced by sloshing motions. The remarkable spatial coincidence observed between the newly-detected arc-like radio filament and the southern concave X-ray discontinuity strongly suggests that the remnant plasma has been displaced by the sloshing motions on large scales. Overall, NGC 507 represents one of the clearest examples known to date in which a direct interaction between old AGN remnant plasma and the external medium is observed in a galaxy group. Our results are consistent with simulations, which suggest that filamentary emission can be created by the cluster/group weather disrupting AGN lobes and spreading their relativistic content into the surrounding medium.

preprint2022arXiv

The LOFAR Two-metre Sky Survey -- V. Second data release

In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for the direction-independent instrumental properties as well as direction-dependent ionospheric distortions during extensive, but fully automated, data processing. A catalogue of 4,396,228 radio sources is derived from our total intensity (Stokes I) maps, where the majority of these have never been detected at radio wavelengths before. At 6" resolution, our full bandwidth Stokes I continuum maps with a central frequency of 144MHz have: a median rms sensitivity of 83$μ$Jy/beam; a flux density scale accuracy of approximately 10%; an astrometric accuracy of 0.2"; and we estimate the point-source completeness to be 90% at a peak brightness of 0.8mJy/beam. By creating three 16MHz bandwidth images across the band we are able to measure the in-band spectral index of many sources, albeit with an error on the derived spectral index of +/-0.2 which is a consequence of our flux-density scale accuracy and small fractional bandwidth. Our circular polarisation (Stokes V) 20" resolution 120-168MHz continuum images have a median rms sensitivity of 95$μ$Jy/beam, and we estimate a Stokes I to Stokes V leakage of 0.056%. Our linear polarisation (Stokes Q and Stokes U) image cubes consist of 480 x 97.6 kHz wide planes and have a median rms sensitivity per plane of 10.8mJy/beam at 4' and 2.2mJy/beam at 20"; we estimate the Stokes I to Stokes Q/U leakage to be approximately 0.2%. Here we characterise and publicly release our Stokes I, Q, U and V images in addition to the calibrated uv-data.

preprint2022arXiv

The Planck clusters in the LOFAR sky. I. LoTSS-DR2: new detections and sample overview

Relativistic electrons and magnetic fields permeate the intra-cluster medium (ICM) and manifest themselves as diffuse sources of synchrotron emission observable at radio wavelengths, namely radio halos and radio relics. Although there is broad consensus that the formation of these sources is connected to turbulence and shocks in the ICM, the details of the required particle acceleration, the strength and morphology of the magnetic field in the cluster volume, and the influence of other sources of high-energy particles are poorly known. Sufficiently large samples of radio halos and relics, which would allow us to examine the variation among the source population and pinpoint their commonalities and differences, are still missing. At present, large numbers of these sources are easiest to detect at low radio frequencies, where they shine brightly. We examined the low-frequency radio emission from all 309 clusters in the second catalog of Planck Sunyaev Zel'dovich detected sources that lie within the 5634 deg$^2$ covered by the Second Data Release of the LOFAR Two-meter Sky Survey (LoTSS-DR2). We produced LOFAR images at different resolutions, with and without discrete sources subtracted, and created overlays with optical and X-ray images before classifying the diffuse sources in the ICM, guided by a decision tree. Overall, we found 83 clusters that host a radio halo and 26 that host one or more radio relics (including candidates). About half of them are new discoveries. The detection rate of clusters hosting a radio halo and one or more relics in our sample is $30\pm11$% and $10\pm6$%, respectively. Extrapolating these numbers, we anticipate that once LoTSS covers the entire northern sky it will provide the detection of $251\pm92$ clusters hosting a halo and $83\pm50$ clusters hosting at least one relic from Planck clusters alone.

preprint2021arXiv

A LOFAR view into the stormy environment of the galaxy cluster 2A0335+096

Radio observations represent a powerful probe of the physics occurring in the intracluster medium (ICM) because they trace the relativistic cosmic rays in the cluster magnetic fields, or within galaxies themselves. By probing the low-energy cosmic rays, low-frequency radio observations are especially interesting because they unveil emission powered by low-efficiency particle acceleration processes, which are believed to play a crucial role in the origin of diffuse radio emission. We investigate the origin of the radio mini-halo at the center of the galaxy cluster 2A0335+096 and its connection to the central galaxy and the sloshing cool core. We also study the properties of the head-tail galaxy GB6 B0335+096 hosted in the cluster to explore the lifecycle of the relativistic electrons in its radio tails. We use new LOw Frequency ARray (LOFAR) observations from the LOFAR Two-meter Sky Survey at 144 MHz to map the low-frequency emission with a high level of detail. The new data were combined with archival Giant Metrewave Radio Telescope and Chandra observations to carry out a multi-wavelength study. We have made the first measurement of the spectral index of the mini-halo ($α=-1.2\pm0.1$ between 144 MHz and 1.4 GHz) and the lobes of the central source ($α\simeq-1.5\pm0.1$ between 144 and 610 MHz). Based on the low-frequency radio emission morphology with respect to the thermal ICM, we propose that the origin of the diffuse radio emission is linked to the sloshing of the cool core. The new data revealed the presence of a Mpc-long radio tail associated with GB6 B0335+096. The observed projected length is a factor 3 longer than the expected cooling length, with evidence of flattening in the spectral index trend along the tail. Therefore, we suggest that the electrons toward the end of the tail are kept alive by the ICM gentle re-acceleration.

preprint2021arXiv

Origin of the ring structures in Hercules A -- Sub-arcsecond 144 MHz to 7 GHz observations

The prominent radio source Hercules A features complex structures in its radio lobes. Although it is one of the most comprehensively studied sources in the radio sky, the origin of the ring structures in the Hercules A radio lobes remains an open question. We present the first sub-arcsecond angular resolution images at low frequencies (<300 MHz) of Hercules A, made with the International LOFAR Telescope. With the addition of data from the Karl G. Jansky Very Large Array, we mapped the structure of the lobes from 144 MHz to 7 GHz. We explore the origin of the rings within the lobes of Hercules A, and test whether their properties are best described by a shock model, where shock waves are produced by the jet propagating in the radio lobe, or by an inner-lobe model, where the rings are formed by decelerated jetted plasma. From spectral index mapping our large frequency coverage reveals that the curvature of the different ring spectra increases with distance away from the central active galactic nucleus. We demonstrate that the spectral shape of the rings is consistent with synchrotron aging, which speaks in favor of an inner-lobe model where the rings are formed from the deposition of material from past periods of intermittent core activity.

preprint2021arXiv

Pushing subarcsecond resolution imaging down to 30 MHz with the trans-European International LOFAR Telescope

Relatively little information is available about the Universe at ultra-low radio frequencies, i.e. below 50 MHz (ULF), although the ULF spectral window contains a wealth of unique diagnostics for studying galactic and extragalactic phenomena. Sub-arcsecond resolution imaging at these frequencies is extremely difficult, due to the long baselines (>1000 km) required and large ionospheric perturbations. We have conducted a pilot project to investigate the ULF performance and potential of the International LOFAR Telescope (ILT), a trans-European interferometric array with baselines up to ~2000 km and observing frequencies down to 10 MHz. We have successfully produced images with sub-arcsecond resolution for 6 radio sources at frequencies down to 30 MHz. This is more than an order of magnitude better resolution than pre-ILT observations at similar frequencies. The six targets that we have imaged (3C 196, 3C 225, 3C 273, 3C 295, 3C 298 and 3C 380) are bright radio sources with compact structures. By comparing our data of 3C 196 and 3C 273 with observations at higher frequencies, we investigate their spatially resolved radio spectral properties. Our success shows that at frequencies down to 30 MHz, sub-arcsecond imaging with the ILT is possible. Further analysis is needed to determine the feasibility of observations of fainter sources or sources with less compact emission.

preprint2021arXiv

Sub-arcsecond imaging with the International LOFAR Telescope I. Foundational calibration strategy and pipeline

[abridged] The International LOFAR Telescope is an interferometer with stations spread across Europe. With baselines of up to ~2,000 km, LOFAR has the unique capability of achieving sub-arcsecond resolution at frequencies below 200 MHz, although this is technically and logistically challenging. Here we present a calibration strategy that builds on previous high-resolution work with LOFAR. We give an overview of the calibration strategy and discuss the special challenges inherent to enacting high-resolution imaging with LOFAR, and describe the pipeline, which is publicly available, in detail. We demonstrate the calibration strategy by using the pipeline on P205+55, a typical LOFAR Two-metre Sky Survey (LoTSS) pointing. We perform in-field delay calibration, solution referencing to other calibrators, self-calibration, and imaging of example directions of interest in the field. For this specific field and these ionospheric conditions, dispersive delay solutions can be transferred between calibrators up to ~1.5 degrees away, while phase solution transferral works well over 1 degree. We demonstrate a check of the astrometry and flux density scale. Imaging in 17 directions, the restoring beam is typically 0.3" x 0.2" although this varies slightly over the entire 5 square degree field of view. We achieve ~80 to 300 $μ$Jy/bm image rms noise, which is dependent on the distance from the phase centre; typical values are ~90 $μ$Jy/bm for the 8 hour observation with 48 MHz of bandwidth. Seventy percent of processed sources are detected, and from this we estimate that we should be able to image ~900 sources per LoTSS pointing. This equates to ~3 million sources in the northern sky, which LoTSS will entirely cover in the next several years. Future optimisation of the calibration strategy for efficient post-processing of LoTSS at high resolution (LoTSS-HR) makes this estimate a lower limit.

preprint2021arXiv

The LOFAR LBA Sky Survey: Deep Fields I. The Boötes Field

We present the first sub-mJy ($\approx0.7$ mJy beam$^{-1}$) survey to be completed below 100 MHz, which is over an order of magnitude deeper than previously achieved for widefield imaging of any field at these low frequencies. The high resolution ($15 \times 15$ arcsec) image of the Boötes field at 34-75 MHz is made from 56 hours of observation with the LOw Frequency ARray (LOFAR) Low Band Antenna (LBA) system. The observations and data reduction, including direction-dependent calibration, are described here. We present a radio source catalogue containing 1,948 sources detected over an area of $23.6$ deg$^2$, with a peak flux density threshold of $5σ$. Using existing datasets, we characterise the astrometric and flux density uncertainties, finding a positional uncertainty of $\sim1.2$ arcsec and a flux density scale uncertainty of about 5 per cent. Using the available deep 144-MHz data, we identified 144-MHz counterparts to all the 54-MHz sources, and produced a matched catalogue within the deep optical coverage area containing 829 sources. We calculate the Euclidean-normalised differential source counts and investigate the low-frequency radio source spectral indices between 54 and 144 MHz, both of which show a general flattening in the radio spectral indices for lower flux density sources, from $\sim-0.75$ at 144-MHz flux densities between 100-1000 mJy to $\sim-0.5$ at 144-MHz flux densities between 5-10 mJy, due to a growing population of star forming galaxies and compact core-dominated AGN.

preprint2020arXiv

A giant radio bridge connecting two clusters in Abell 1758

Collisions between galaxy clusters dissipate enormous amounts of energy in the intra-cluster medium (ICM) through turbulence and shocks. In the process, Mpc-scale diffuse synchrotron emission in form of radio halos and relics can form. However, little is known about the very early phase of the collision. We used deep radio observations from 53 MHz to 1.5 GHz to study the pre-merging galaxy clusters A1758N and A1758S that are $\sim2$ Mpc apart. We confirm the presence of a giant bridge of radio emission connecting the two systems that was reported only tentatively in our earlier work. This is the second large-scale radio bridge observed to date in a cluster pair. The bridge is clearly visible in the LOFAR image at 144 MHz and tentatively detected at 53 MHz. Its mean radio emissivity is more than one order of magnitude lower than that of the radio halos in A1758N and A1758S. Interestingly, the radio and X-ray emissions of the bridge are correlated. Our results indicate that non-thermal phenomena in the ICM can be generated also in the region of compressed gas in-between infalling systems.

preprint2020arXiv

Coherent radio emission from a quiescent red dwarf indicative of star-planet interaction

Low frequency ($ν\lesssim 150\,$MHz) stellar radio emission is expected to originate in the outer corona at heights comparable to and larger than the stellar radius. Such emission from the Sun has been used to study coronal structure, mass ejections, space-weather conditions around the planets (Schwenn 2006). Searches for low-frequency emission from other stars have only detected a single active flare-star (Lynch et al. 2017) that is not representative of the wider stellar population. Here we report the detection of low-frequency radio emission from a quiescent star, GJ 1151--- a member of the most common stellar type (red dwarf or spectral class M) in the Galaxy. The characteristics of the emission are similar to those of planetary auroral emissions (Zarka 1998) (e.g. Jupiter's decametric emission), suggesting a coronal structure dominated by a global magnetosphere with low plasma density. Our results show that large-scale currents that power radio aurorae operate over a vast range of mass and atmospheric composition, ranging from terrestrial planets to main-sequence stars. The Poynting flux required to produce the observed radio emission cannot be generated by GJ 1151's slow rotation, but can originate in a sub-Alfvénic interaction of its magnetospheric plasma with a short-period exoplanet. The emission properties are consistent with theoretical expectations (Zarka 2007; Lanza 2009; Saur et al. 2013; Turnpenney et al. 2018) for interaction with an Earth-size planet in a $\sim 1-5$ day-long orbit.

preprint2020arXiv

LOFAR 144-MHz follow-up observations of GW170817

We present low-radio-frequency follow-up observations of AT 2017gfo, the electromagnetic counterpart of GW170817, which was the first binary neutron star merger to be detected by Advanced LIGO-Virgo. These data, with a central frequency of 144 MHz, were obtained with LOFAR, the Low-Frequency Array. The maximum elevation of the target is just 13.7 degrees when observed with LOFAR, making our observations particularly challenging to calibrate and significantly limiting the achievable sensitivity. On time-scales of 130-138 and 371-374 days after the merger event, we obtain 3$σ$ upper limits for the afterglow component of 6.6 and 19.5 mJy beam$^{-1}$, respectively. Using our best upper limit and previously published, contemporaneous higher-frequency radio data, we place a limit on any potential steepening of the radio spectrum between 610 and 144 MHz: the two-point spectral index $α^{610}_{144} \gtrsim -2.5$. We also show that LOFAR can detect the afterglows of future binary neutron star merger events occurring at more favourable elevations.

preprint2020arXiv

LOFAR Observations of X-ray Cavity Systems

We present LOFAR observations at 120-168 MHz of 42 systems with possible X-ray cavities in their hot atmosphere, of which 17 are groups or ellipticals, 19 are nearby clusters (z<0.3), and six are higher-redshift clusters (z>0.3). The X-ray cavities, formed by the radio lobes of the central active galactic nucleus (AGN), are evidence of radio-mode AGN feedback. In the groups and ellipticals sample, more than half of the systems have X-ray cavities for which no associated lobe emission was detected. Conversely, we report the discovery of large radio lobes in NGC 6338, extending far beyond the emission reported previously in the literature. In the case of the nearby clusters, our observations show that there is little low-frequency radio emission that extends beyond the cavities (e.g., MS 0735.6+7421 and A2052). For the first time, we report secure cavity-radio associations in 2A 0335+096, ZwCl 2701, and ZwCl 8276 that strengthens their interpretation as AGN-created cavities. However, in some known cavity systems (e.g., A1795 and ZwCl 3146) we report the lack of detectable low-frequency radio emission associated with the cavities. Our sample of higher-redshift systems is small, and unfortunately the present LOFAR observations are not able to resolve the lobes in many of them. Nevertheless, our sample represents one of the best available for investigating the connection between radio and jet power in radio-mode AGN feedback.

preprint2020arXiv

Physical insights from the spectrum of the radio halo in MACS J0717.5+3745

We present new LOFAR observations of the massive merging galaxy cluster MACS J0717.5+3745. The cluster hosts the most powerful radio halo known to date. These new observations, in combination with published uGMRT (300$-$850 MHz) and VLA (1$-$6.5 GHz) data, reveal that the halo is more extended than previously thought, with a largest linear size of $\sim2.2 \rm Mpc$. The halo shows a steep spectrum ($α_{144\,\text{MHz}}^{1.5\,\text{GHz}}\sim-1.4$) and a steepening ($α_{1.5 \text{GHz}}^{5.5 \text{GHz}}\sim-1.9$) above 1.5 GHz. We find a strong scattering in spectral index maps on scales of 50$-$100 kpc. We suggest that such a strong scattering may be a consequence of the regime where inverse Compton dominate the energy losses of electrons. The spectral index becomes steeper and shows an increased curvature in the outermost regions of the halo. We combined the radio data with \textit{Chandra} observations to investigate the connection between the thermal and non-thermal components of the intracluster medium (ICM). Despite a significant substructure in the halo emission, the radio brightness correlates strongly with the X-ray brightness at all observed frequencies. The radio-versus-X-ray brightness correlation slope steepens at a higher radio frequency (from $b_{144 \text{MHz}}=0.67\pm0.05$ to $b_{3.0 \text{GHz}}=0.98\pm0.09$) and the spectral index shows a significant anti correlation with the X-ray brightness. Both pieces of evidence further support a spectral steepening in the external regions. The compelling evidence for a steep spectral index, the existence of a spectral break above 1.5 GHz, and the dependence of radio and X-ray surface brightness correlation on frequency are interpreted in the context of turbulent reacceleration models. Under this scenario, our results allowed us to constrain that the turbulent kinetic pressure of the ICM is up to 10%.

preprint2020arXiv

Probabilistic direction-dependent ionospheric calibration for LOFAR-HBA

Direction dependent calibration and imaging is a vital part of producing deep, high fidelity, high-dynamic range radio images with a wide-field low-frequency array like LOFAR. Currently, state-of-the-art facet-based direction dependent calibration algorithms rely on the assumption that the isoplanatic-patch size is much larger than the separation between bright in-field calibrators. This assumption is often violated due to the dynamic nature of the ionosphere, and as a result direction dependent errors affect image quality between calibrators. In this paper we propose a probabilistic physics-informed model for inferring ionospheric phase screens, providing a calibration for all sources in the field of view. We apply our method to a randomly selected observation from the LOFAR Two-Metre Sky Survey archive, and show that almost all direction dependent effects between bright calibrators are corrected and that the root-mean-squared residuals around bright sources is reduced by 32\% on average.

preprint2020arXiv

The beautiful mess in Abell 2255

We present LOFAR observations of one of the most spectacular objects in the radio sky: Abell 2255. This is a nearby ($z = 0.0806$) merging galaxy cluster hosting one of the first radio halos ever detected in the intra-cluster medium (ICM). The deep LOFAR images at 144 MHz of the central $\sim10$ Mpc$^2$ region show a plethora of emission on different scales, from tens of kpc to above Mpc sizes. In this work, we focus on the innermost region of the cluster. Among the numerous interesting features observed, we discover remarkable bright and filamentary structures embedded in the radio halo. We incorporate archival WSRT 1.2 GHz data to study the spectral properties of the diffuse synchrotron emission and find a very complex spectral index distribution in the halo spanning a wide range of values. We combine the radio data with Chandra observations to investigate the connection between the thermal and non-thermal components by quantitatively comparing the radio and X-ray surface brightness and the spectral index of the radio emission with the thermodynamical quantities of the ICM. Despite the multitude of structures observed in the radio halo, we find that the X-ray and radio emission are overall well correlated. The fact that the steepest spectrum emission is located in the cluster center and traces regions with high entropy possibly suggests the presence of seed particles injected by radio galaxies that are spread in the ICM by turbulence generating the extended radio halo.

preprint2020arXiv

The life cycle of radio galaxies in the LOFAR Lockman Hole field

Radio galaxies are known to go through cycles of activity, where phases of apparent quiescence can be followed by repeated activity of the central supermassive black hole. A better understanding of this cycle is crucial for ascertaining the energetic impact that the jets have on the host galaxy, but little is known about it. We used deep LOFAR images at 150 MHz of the Lockman Hole extragalactic field to select a sample of 158 radio sources with sizes $> 60^{\prime\prime}$ in different phases of their jet life cycle. Using a variety of criteria (e.g. core prominence combined with low-surface brightness of the extended emission and steep spectrum of the central region) we selected a subsample of candidate restarted radio galaxies representing between 13% and 15% of the 158 sources of the main sample. We compare their properties to the rest of the sample, which consists of remnant candidates and active radio galaxies. Optical identifications and characterisations of the host galaxies indicate similar properties for candidate restarted, remnant, and active radio galaxies, suggesting that they all come from the same parent population. The fraction of restarted radio galaxies is slightly higher with respect to remnants, suggesting that the restarted phase can often follow after a relatively short remnant phase (the duration of the remnant phase being a few times 10$^{7}$ years). This confirms that the remnant and restarted phases are integral parts of the life cycle of massive elliptical galaxies. A preliminary investigation does not suggest a strong dependence of this cycle on the environment surrounding any given galaxy.

preprint2020arXiv

The optical luminosity function of LOFAR radio-selected quasars at $1.4\leq z\leq5.0$ in the NDWFS-Boötes field

We present an estimate of the optical luminosity function (OLF) of LOFAR radio-selected quasars (RSQs) at $1.4<z<5.0$ in the $9.3\:\textrm{deg}^{2}$ NOAO Deep Wide-field survey (NDWFS) of the Boötes field. The selection was based on optical/mid-ir photometry used to train three different machine learning (ML) algorithms. Objects taken as quasars by the ML algorithms are required to be detected at $5σ$ significance in deep radio maps to be classified as candidate quasars. The optical imaging came from the SDSS and the PS1 $3π$ survey; mid-ir photometry was taken from the SDWFS survey; and radio data was obtained from deep LOFAR imaging of the NDWFS-Boötes field. The requirement of a $5σ$ LOFAR detection allowed us to reduce the stellar contamination in our sample by two orders of magnitude. The sample comprises 134 objects, including both photometrically selected candidate quasars (47) and spectroscopically confirmed quasars (83). The depth of our LOFAR observations allowed us to detect the radio-emission of quasars that would be otherwise classified as radio-quiet. Around $65\%$ of the quasars in the sample are fainter than $M_{\textrm{1450}}=-24.0$, a regime where the OLF of quasars selected through their radio emission, has not been investigated in detail. It has been demonstrated that in cases where mid-ir wedge-based AGN selection is not possible due to a lack of appropriate data, the selection of quasars using ML algorithms trained with optical/mid-ir photometry in combination with LOFAR data provides an excellent approach for obtaining samples of quasars. We demonstrate that RSQs show an evolution similar to that exhibited by faint quasars $(M_{\textrm{1450}}\leq-22.0)$. Finally, we find that RSQs may compose up to $\sim20\%$ of the whole faint quasar population (radio-detected plus radio-undetected).

preprint2020arXiv

Tied-array holography with LOFAR

A radio interferometer uses time delays to maximize its response to radiation coming from a particular direction. These time delays compensate for differences in the time of arrival of the wavefront at the different elements of the interferometer, and for delays in the instrument's signal chain. If the radio interferometer is operated as a phased array (tied array), the time delays cannot be accounted for after an observation, so they must be determined in advance. Our aim is to characterize the time delays between the stations in the core of the LOw Frequency ARray (LOFAR). We used radio holography to determine the time delays for the core stations of LOFAR (innermost 3.5 km). Using the multibeaming capability of LOFAR we map the voltage beam faster than with a raster scan, while simultaneously calibrating the observed beam continuously. For short radio holographic observations (60 s and 600 s) of 3C196, 3C147, and 3C48 we are able to derive time delays with errors of less than one nanosecond. After applying the derived time delays to the beamformer, the beam shows residuals of less than $20\%$ with respect to the theoretical beam shape. Tied-array holography could be a way towards semi-real-time beam calibration for the Square Kilometer Array.

preprint2020arXiv

Understanding the radio relic emission in the galaxy cluster MACS J0717.5+3745: spectral analysis

Radio relics are diffuse, extended synchrotron sources that originate from shock fronts generated during cluster mergers. The massive merging galaxy cluster MACS J0717.5+3745 hosts one of the more complex relics known to date. We present upgraded Giant Metrewave Radio Telescope band 3 (300-500 MHz) and band 4 (550-850 MHz) observations. These new observations, combined with published VLA and the new LOFAR HBA data, allow us to carry out a detailed, high spatial resolution spectral analysis of the relic over a broad range of frequencies. The integrated spectrum of the relic closely follows a power-law between 144 MHz and 5.5 GHz with a mean spectral slope $α=-1.16\pm0.03$. Despite its complex morphology, the subregions of the relic and the other isolated filaments also follow power-law behaviors, and show similar spectral slopes. Assuming Diffusive Shock Acceleration, we estimate a dominant Mach number of $\sim 3.7$ for the shocks that make up the relic. Comparison with recent numerical simulations suggests that in the case of radio relics, the slopes of the integrated radio spectra are determined by the Mach number of the accelerating shock, with $α$ nearly constant, namely between $-1.13$ and $-1.17$, for Mach numbers $3.5 - 4.0$. The spectral shapes inferred from spatially resolved regions show curvature, we speculate that the relic is inclined along the line-of-sight. The locus of points in the simulated color-color plots changes significantly with the relic viewing angle. We conclude that projection effects and inhomogeneities in the shock Mach number dominate the observed spectral properties of the relic in this complex system. Based on the new observations we raise the possibility that the relic and a narrow-angle-tailed radio galaxy are two different structures projected along the same line-of-sight.

preprint2020arXiv

Very Large Array observations of the mini-halo and AGN feedback in the Phoenix cluster

(Abridged) The relaxed cool-core Phoenix cluster (SPT-CL J2344-4243) features an extremely strong cooling flow, as well as a mini-halo. Strong star-formation in the brightest cluster galaxy indicates that AGN feedback has been unable to inhibit this cooling flow. We have studied the strong cooling flow in the Phoenix cluster by determining the radio properties of the AGN and its lobes. In addition, we use spatially resolved observations to investigate the origin of the mini-halo. We present new Very Large Array 1-12 GHz observations of the Phoenix cluster which resolve the AGN and its lobes in all four frequency bands, and resolve the mini-halo in L- and S-band. Using our L-band observations, we measure the total flux density of the radio lobes at 1.5 GHz to be $7.6\pm0.8$ mJy, and the flux density of the mini-halo to be $8.5\pm0.9$ mJy. Using L- and X-band images, we produce the first spectral index maps of the lobes from the AGN and measure the spectral indices of the northern and southern lobes to be $-1.35\pm0.07$ and $-1.30\pm0.12$, respectively. Similarly, using L- and S-band data, we map the spectral index of the mini-halo, and obtain an integrated spectral index of $α=-0.95 \pm 0.10$. We find that the mini-halo is most likely formed by turbulent re-acceleration powered by sloshing in the cool core due to a recent merger. In addition, we find that the feedback in the Phoenix cluster is consistent with the picture that stronger cooling flows are to be expected for massive clusters like the Phoenix cluster, as these may feature an underweight supermassive black hole due to their merging history. Strong time variability of the AGN on Myr-timescales may help explain the disconnection between the radio and the X-ray properties of the system. Finally, a small amount of jet precession likely contributes to the relatively low ICM re-heating efficiency of the mechanical feedback.

preprint2019arXiv

A probabilistic approach to direction-dependent ionospheric calibration

Calibrating for direction-dependent ionospheric distortions in visibility data is one of the main technical challenges that must be overcome to advance low-frequency radio astronomy. In this paper, we propose a novel probabilistic, tomographic approach that utilises Gaussian processes to calibrate direction-dependent ionospheric phase distortions in low-frequency interferometric data. We suggest that the ionospheric free electron density can be modelled to good approximation by a Gaussian process restricted to a thick single layer, and show that under this assumption the differential total electron content must also be a Gaussian process. We perform a comparison with a number of other widely successful Gaussian processes on simulated differential total electron contents over a wide range of experimental conditions, and find that, in all experimental conditions, our model is better able to represent observed data and generalise to unseen data. The mean equivalent source shift imposed by our predictive errors are half as large as the best competitor model's. We find that it is possible to partially constrain the ionosphere's hyperparameters from sparse-and-noisy observed data. Our model provides an alternative explanation for observed phase structure functions deviating from Kolmogorov's 5/3 turbulence, turnover at high baselines, and diffractive scale anisotropy. We show that our model implicitly cheaply performs tomography of the free electron density. Moreover, we find that even a fast, low-resolution approximation to our model yields better results than the best alternative Gaussian process, implying that the geometric coupling between directions and antennae is a powerful prior that should not be ignored.

preprint2019arXiv

Revived Fossil Plasma Sources in Galaxy Clusters

It is well established that particle acceleration by shocks and turbulence in the intra-cluster medium can produce cluster-scale synchrotron emitting sources. However, the detailed physics of these particle acceleration processes is still not well understood. One of the main open questions is the role of fossil relativistic electrons that have been deposited in the intra-cluster medium by radio galaxies. These synchrotron-emitting electrons are very difficult to study, as their radiative life time is only tens of Myrs at GHz frequencies, and are therefore a relatively unexplored population. Despite the typical steep radio spectrum due to synchrotron losses, these fossil electrons are barely visible even at radio frequencies well below a GHz. However, when a pocket of fossil radio plasma is compressed, it boosts the visibility at sub-GHz frequencies, creating so-called radio phoenices. This compression can be the result of bulk motion and shocks in the ICM due to merger activity. In this paper, we demonstrate the discovery potential of low frequency radio sky surveys to find and study revived fossil plasma sources in galaxy clusters. We used the 150~MHz TGSS and 1.4 GHz NVSS sky surveys to identify candidate radio phoenices. A subset of three candidates were studied in detail using deep multi-band radio observations (LOFAR and GMRT), X-ray (\textit{Chandra} or \textit{XMM-Newton}) and archival optical observations. Two of the three sources are new discoveries. Using these observations, we identified common observational properties (radio morphology, ultra-steep spectrum, X-ray luminosity, dynamical state) that will enable us to identify this class of sources more easily, and help to understand the physical origin of these sources.

preprint2019arXiv

VLT/SINFONI study of black hole growth in high redshift radio-loud quasars from the CARLA survey

We present VLT/SINFONI observations of 35 quasars at 2.1 < z < 3.2, the majority of which were selected from the Clusters Around Radio-Loud AGN (CARLA) survey. CARLA quasars have large CIV-based black hole (BH) masses (M(BH) > 10^9 Msun) and powerful radio emission (P(500MHz) > 27.5 W/Hz). We estimate Ha-based M(BH), finding a scatter of 0.35 dex compared to CIV. We evaluate several recipes for correcting CIV-based masses, which reduce the scatter to 0.24 dex. The radio power of the radio-loud quasars is at most weakly correlated with the interconnected quantities Ha-width, L(5100A) and M(BH), suggesting that it is governed by different physical processes. However, we do find a strong inverse correlation between CIV blueshift and radio power linked to higher Eddington ratios and L(5100A). Under standard assumptions, the BH growth time is longer than the cosmic age for many CARLA quasars, suggesting that they must have experienced more efficient growth in the past. If these BHs were growing from seeds since the epoch of reionization, it is possible that they grew at the Eddington limit like the quasars at z ~ 6-7, and then continued to grow at the reduced rates observed until z ~ 2-3. Finally, we study the relation between M(BH) and environment, finding a weak positive correlation between M(BH) and galaxy density measured by CARLA.

preprint2016arXiv

$Suzaku$ and $XMM-Newton$ observations of a newly-discovered early-stage cluster merger 1E2216.0-0401 and 1E2215.7-0404

We present the results of $Suzaku$ and $XMM-Newton$ X-ray observations of the cluster pair 1E2216.0-0401 and 1E2215.7-0404. We discover an X-ray bridge between the clusters. $Suzaku$ and $XMM-Newton$ observations revealed that each cluster hosts gas with moderate temperature of $kT_{1E2216.0-0401}=$4.8$\pm$0.1 keV and $kT_{1E2215.7-0404}=$5.8$\pm$0.2 keV, respectively. On the other hand, the bridge region shows a remarkably high temperature ({\it kT}=6.6$\pm$0.5 keV). Furthermore, at the position of the bridge, we detected an enhancement in the wavelet-decomposed soft-band (0.5-4.0 keV) $XMM-Newton$ image at 3 sigma significance, this is most likely due to a compression of the intracluster medium (ICM) as a consequence of the merging activity. This X-ray intensity and temperature enhancement are not consistent with those expected from a late phase, but are in agreement with the predictions by numerical simulations of an early phase merger. From the temperature jump at the location of the bridge, the Mach number is estimated to be ${\cal M}=1.4\pm0.1$, which corresponds to a shock propagation velocity of about 1570 km/s. From the shock properties, we estimate that core-passage will occur in 0.3-0.6 Gyr and that the age of the shock structure is 50--100 Myr. Based on the measured properties of the ICM at the bridge and estimation of timescales, we find indications for non-equilibrium ionization. We also discover possible diffuse radio emission located between the merging clusters. Combining the radio, X-ray, and optical image data, we speculate that the detected radio sources are most likely related to the merger event. Thus, 1E2216.0-0401 and 1E2215.7-0404 is a new example of an early phase cluster merger with remarkable characteristics.

preprint2016arXiv

A large light-mass component of cosmic rays at 10^{17} - 10^{17.5} eV from radio observations

Cosmic rays are the highest energy particles found in nature. Measurements of the mass composition of cosmic rays between 10^{17} eV and 10^{18} eV are essential to understand whether this energy range is dominated by Galactic or extragalactic sources. It has also been proposed that the astrophysical neutrino signal comes from accelerators capable of producing cosmic rays of these energies. Cosmic rays initiate cascades of secondary particles (air showers) in the atmosphere and their masses are inferred from measurements of the atmospheric depth of the shower maximum, Xmax, or the composition of shower particles reaching the ground. Current measurements suffer from either low precision, or a low duty cycle and a high energy threshold. Radio detection of cosmic rays is a rapidly developing technique, suitable for determination of Xmax with a duty cycle of in principle nearly 100%. The radiation is generated by the separation of relativistic charged particles in the geomagnetic field and a negative charge excess in the shower front. Here we report radio measurements of Xmax with a mean precision of 16 g/cm^2 between 10^{17}-10^{17.5} eV. Because of the high resolution in $Xmax we can determine the mass spectrum and find a mixed composition, containing a light mass fraction of ~80%. Unless the extragalactic component becomes significant already below 10^{17.5} eV, our measurements indicate an additional Galactic component dominating at this energy range.

preprint2016arXiv

A plethora of diffuse steep spectrum radio sources in Abell 2034 revealed by LOFAR

With Low-Frequency Array (LOFAR) observations, we have discovered a diverse assembly of steep spectrum emission that is apparently associated with the intra cluster medium (ICM) of the merging galaxy cluster Abell 2034. Such a rich variety of complex emission associated with the ICM has been observed in few other clusters. This not only indicates that Abell 2034 is a more interesting and complex system than previously thought but it also demonstrates the importance of sensitive and high-resolution, low-frequency observations. These observations can reveal emission from relativistic particles which have been accelerated to sufficient energy to produce observable emission or have had their high energy maintained by mechanisms in the ICM. The most prominent feature in our maps is a bright bulb of emission connected to two steep spectrum filamentary structures, the longest of which extends perpendicular to the merger axis for 0.5Mpc across the south of the cluster. The origin of these objects is unclear, with no shock detected in the X-ray images and no obvious connection with cluster galaxies or AGNs. We also find that the X-ray bright region of the cluster coincides with a giant radio halo with an irregular morphology and a very steep spectrum. In addition, the cluster hosts up to three possible radio relics, which are misaligned with the cluster X-ray emission. Finally, we have identified multiple regions of emission with a very steep spectral index that seem to be associated with either tailed radio galaxies or a shock.

preprint2016arXiv

LOFAR 150-MHz observations of the Boötes field: Catalogue and Source Counts

We present the first wide area (19 deg$^2$), deep ($\approx120-150$ μJy beam$^{-1}$), high resolution ($5.6 \times 7.4$ arcsec) LOFAR High Band Antenna image of the Boötes field made at 130-169 MHz. This image is at least an order of magnitude deeper and 3-5 times higher in angular resolution than previously achieved for this field at low frequencies. The observations and data reduction, which includes full direction-dependent calibration, are described here. We present a radio source catalogue containing 6276 sources detected over an area of $19$\,deg$^2$, with a peak flux density threshold of $5σ$. As the first thorough test of the facet calibration strategy, introduced by van Weeren et al., we investigate the flux and positional accuracy of the catalogue. We present differential source counts that reach an order of magnitude deeper in flux density than previously achieved at these low frequencies, and show flattening at 150 MHz flux densities below 10 mJy associated with the rise of the low flux density star-forming galaxies and radio-quiet AGN.

preprint2016arXiv

LOFAR facet calibration

LOFAR, the Low-Frequency Array, is a powerful new radio telescope operating between 10 and 240 MHz. LOFAR allows detailed sensitive high-resolution studies of the low-frequency radio sky. At the same time LOFAR also provides excellent short baseline coverage to map diffuse extended emission. However, producing high-quality deep images is challenging due to the presence of direction dependent calibration errors, caused by imperfect knowledge of the station beam shapes and the ionosphere. Furthermore, the large data volume and presence of station clock errors present additional difficulties. In this paper we present a new calibration scheme, which we name facet calibration, to obtain deep high-resolution LOFAR High Band Antenna images using the Dutch part of the array. This scheme solves and corrects the direction dependent errors in a number of facets that cover the observed field of view. Facet calibration provides close to thermal noise limited images for a typical 8 hr observing run at $\sim$ 5arcsec resolution, meeting the specifications of the LOFAR Tier-1 northern survey.

preprint2016arXiv

LOFAR, VLA, and Chandra observations of the Toothbrush galaxy cluster

We present deep LOFAR observations between 120-181 MHz of the "Toothbrush" (RX J0603.3+4214), a cluster that contains one of the brightest radio relic sources known. Our LOFAR observations exploit a new and novel calibration scheme to probe 10 times deeper than any previous study in this relatively unexplored part of the spectrum. The LOFAR observations, when combined with VLA, GMRT, and Chandra X-ray data, provide new information about the nature of cluster merger shocks and their role in re-accelerating relativistic particles. We derive a spectral index of $α= -0.8 \pm 0.1$ at the northern edge of the main radio relic, steepening towards the south to $α\approx - 2$. The spectral index of the radio halo is remarkably uniform ($α= -1.16$, with an intrinsic scatter of $\leq 0.04$). The observed radio relic spectral index gives a Mach number of $\mathcal{M} = 2.8^{+0.5}_{-0.3}$, assuming diffusive shock acceleration (DSA). However, the gas density jump at the northern edge of the large radio relic implies a much weaker shock ($\mathcal{M} \approx 1.2$, with an upper limit of $\mathcal{M} \approx 1.5$). The discrepancy between the Mach numbers calculated from the radio and X-rays can be explained if either (i) the relic traces a complex shock surface along the line of sight, or (ii) if the radio relic emission is produced by a re-accelerated population of fossil particles from a radio galaxy. Our results highlight the need for additional theoretical work and numerical simulations of particle acceleration and re-acceleration at cluster merger shocks.

preprint2016arXiv

LOFAR/H-ATLAS: A deep low-frequency survey of the Herschel-ATLAS North Galactic Pole field

We present LOFAR High-Band Array (HBA) observations of the Herschel-ATLAS North Galactic Pole survey area. The survey we have carried out, consisting of four pointings covering around 142 square degrees of sky in the frequency range 126--173 MHz, does not provide uniform noise coverage but otherwise is representative of the quality of data to be expected in the planned LOFAR wide-area surveys, and has been reduced using recently developed `facet calibration' methods at a resolution approaching the full resolution of the datasets ($\sim 10 \times 6$ arcsec) and an rms off-source noise that ranges from 100 $μ$Jy beam$^{-1}$ in the centre of the best fields to around 2 mJy beam$^{-1}$ at the furthest extent of our imaging. We describe the imaging, cataloguing and source identification processes, and present some initial science results based on a 5-$σ$ source catalogue. These include (i) an initial look at the radio/far-infrared correlation at 150 MHz, showing that many Herschel sources are not yet detected by LOFAR; (ii) number counts at 150 MHz, including, for the first time, observational constraints on the numbers of star-forming galaxies; (iii) the 150-MHz luminosity functions for active and star-forming galaxies, which agree well with determinations at higher frequencies at low redshift, and show strong redshift evolution of the star-forming population; and (iv) some discussion of the implications of our observations for studies of radio galaxy life cycles.

preprint2016arXiv

The Lockman Hole project: LOFAR observations and spectral index properties of low-frequency radio sources

The Lockman Hole is a well-studied extragalactic field with extensive multi-band ancillary data covering a wide range in frequency, essential for characterising the physical and evolutionary properties of the various source populations detected in deep radio fields (mainly star-forming galaxies and AGNs). In this paper we present new 150-MHz observations carried out with the LOw Frequency ARray (LOFAR), allowing us to explore a new spectral window for the faint radio source population. This 150-MHz image covers an area of 34.7 square degrees with a resolution of 18.6$\times$14.7 arcsec and reaches an rms of 160 $μ$Jy beam$^{-1}$ at the centre of the field. As expected for a low-frequency selected sample, the vast majority of sources exhibit steep spectra, with a median spectral index of $α_{150}^{1400}=-0.78\pm0.015$. The median spectral index becomes slightly flatter (increasing from $α_{150}^{1400}=-0.84$ to $α_{150}^{1400}=-0.75$) with decreasing flux density down to $S_{150} \sim$10 mJy before flattening out and remaining constant below this flux level. For a bright subset of the 150-MHz selected sample we can trace the spectral properties down to lower frequencies using 60-MHz LOFAR observations, finding tentative evidence for sources to become flatter in spectrum between 60 and 150 MHz. Using the deep, multi-frequency data available in the Lockman Hole, we identify a sample of 100 Ultra-steep spectrum (USS) sources and 13 peaked spectrum sources. We estimate that up to 21 percent of these could have $z>4$ and are candidate high-$z$ radio galaxies, but further follow-up observations are required to confirm the physical nature of these objects.

preprint2016arXiv

The WEAVE-LOFAR Survey

In these proceedings we highlight the primary scientific goals and design of the WEAVE-LOFAR survey, which will use the new WEAVE spectrograph on the 4.2m William Herschel Telescope to provide the primary source of spectroscopic information for the LOFAR Surveys Key Science Project. Beginning in 2018, WEAVE-LOFAR will generate more than 10$^6$ R=5000 365-960 nm spectra of low-frequency selected radio sources, across three tiers designed to efficiently sample the redshift-luminosity plane, and produce a data set of enormous legacy value. The radio frequency selection, combined with the high multiplex and throughput of the WEAVE spectrograph, make obtaining redshifts in this way very efficient, and we expect that the redshift success rate will approach 100 per cent at $z < 1$. This unprecedented spectroscopic sample - which will be complemented by an integral field component - will be transformational in key areas, including studying the star formation history of the Universe, the role of accretion and AGN-driven feedback, properties of the epoch of reionisation, cosmology, cluster haloes and relics, as well as the nature of radio galaxies and protoclusters. Each topic will be addressed in unprecedented detail, and with the most reliable source classifications and redshift information in existence.

preprint2015arXiv

AGN duty cycle estimates for the ultra-steep spectrum radio relic VLSS J1431.8+1331

Steep spectrum radio sources associated with active galactic nuclei (AGN) may contain remnants of past AGN activity episodes. Novel instruments like the LOw Frequency ARray (LOFAR) are enabling studies of these fascinating structures to be made at tens to hundreds of MHz with sufficient resolution to analyse their complex morphology. Our goal is to characterize the integrated and resolved spectral properties of VLSS J1431+1331 and estimate source ages based on synchrotron radio emission models, thus putting constraints on the AGN duty cycle. Using a broad spectral coverage, we have derived spectral and curvature maps, and used synchrotron ageing models to determine the time elapsed from the last time the source plasma was energized. We used LOFAR, Giant Metrewave Radio Telescope (GMRT) and Jansky Very Large Array (VLA) data. Based on our ageing analysis, we infer that the AGN that created this source currently has very low levels of activity or that it is switched off. The derived ages for the larger source component range from around 60 to 130 Myr, hinting that the AGN activity decreased or stopped around 60 Myr ago. Our analysis suggests that VLSS J1431.8+1331 is an intriguing, two-component source. The larger component seems to host a faint radio core, suggesting that the source may be an AGN radio relic. The spectral index we observe from the smaller component is distinctly flatter at lower frequencies than the spectral index of the larger component, suggesting the possibility that the smaller component may be a shocked plasma bubble. From the integrated source spectrum, we deduce that its shape and slope can be used as tracers of the activity history of this type of steep spectrum radio source.

preprint2015arXiv

Cosmic downsizing of powerful radio galaxies to low radio luminosities

At bright radio powers ($P_{\rm 1.4 GHz} > 10^{25}$ W/Hz) the space density of the most powerful sources peaks at higher redshift than that of their weaker counterparts. This paper establishes whether this luminosity-dependent evolution persists for sources an order of magnitude fainter than those previously studied, by measuring the steep--spectrum radio luminosity function (RLF) across the range $10^{24} < P_{\rm 1.4 GHz} < 10^{28}$ W/Hz, out to high redshift. A grid-based modelling method is used, in which no assumptions are made about the RLF shape and high-redshift behaviour. The inputs to the model are the same as in Rigby et al. (2011): redshift distributions from radio source samples, together with source counts and determinations of the local luminosity function. However, to improve coverage of the radio power vs. redshift plane at the lowest radio powers, a new faint radio sample is introduced. This covers 0.8 sq. deg., in the Subaru/XMM-Newton Deep Field, to a 1.4 GHz flux density limit of $S_{\rm 1.4 GHz} \geq 100~μ$Jy, with 99% redshift completeness. The modelling results show that the previously seen high-redshift declines in space density persist to $P_{\rm 1.4 GHz} < 10^{25}$ W/Hz. At $P_{\rm 1.4 GHz} > 10^{26}$ W/Hz the redshift of the peak space density increases with luminosity, whilst at lower radio luminosities the position of the peak remains constant within the uncertainties. This `cosmic downsizing' behaviour is found to be similar to that seen at optical wavelengths for quasars, and is interpreted as representing the transition from radiatively efficient to inefficient accretion modes in the steep-spectrum population. This conclusion is supported by constructing simple models for the space density evolution of these two different radio galaxy classes; these are able to successfully reproduce the observed variation in peak redshift.

preprint2015arXiv

LOFAR discovery of a 700-kpc remnant radio galaxy at low redshift

Remnant radio galaxies represent the final dying phase of radio galaxy evolution, in which the jets are no longer active. Due to their rarity in flux limited samples and the difficulty of identification, this dying phase remains poorly understood and the luminosity evolution largely unconstrained. Here we present the discovery and detailed analysis of a large (700 kpc) remnant radio galaxy with a low surface brightness that has been identified in LOFAR images at 150 MHz. By combining LOFAR data with new follow-up Westerbork observations and archival data at higher frequencies, we investigated the source morphology and spectral properties from 116 to 4850 MHz. By modelling the radio spectrum we probed characteristic timescales of the radio activity. The source has a relatively smooth, diffuse, amorphous appearance together with a very weak central compact core which is associated with the host galaxy located at z=0.051. From our ageing and morphological analysis it is clear that the nuclear engine is currently switched off or, at most, active at a very low power state. The host galaxy is currently interacting with another galaxy located at a projected separation of 15 kpc and a radial velocity offset of 300 km/s. This interaction may have played a role in the triggering and/or shut down of the radio jets. The spectral shape of this remnant radio galaxy differs from the majority of the previously identified remnant sources, which show steep or curved spectra at low to intermediate frequencies. In light of this finding and in preparation for new-generation deep low-frequency surveys, we discuss the selection criteria to be used to select representative samples of these sources.

preprint2015arXiv

LOFAR MSSS: Detection of a low-frequency radio transient in 400 hrs of monitoring of the North Celestial Pole

We present the results of a four-month campaign searching for low-frequency radio transients near the North Celestial Pole with the Low-Frequency Array (LOFAR), as part of the Multifrequency Snapshot Sky Survey (MSSS). The data were recorded between 2011 December and 2012 April and comprised 2149 11-minute snapshots, each covering 175 deg^2. We have found one convincing candidate astrophysical transient, with a duration of a few minutes and a flux density at 60 MHz of 15-25 Jy. The transient does not repeat and has no obvious optical or high-energy counterpart, as a result of which its nature is unclear. The detection of this event implies a transient rate at 60 MHz of 3.9 (+14.7, -3.7) x 10^-4 day^-1 deg^-2, and a transient surface density of 1.5 x 10^-5 deg^-2, at a 7.9-Jy limiting flux density and ~10-minute time-scale. The campaign data were also searched for transients at a range of other time-scales, from 0.5 to 297 min, which allowed us to place a range of limits on transient rates at 60 MHz as a function of observation duration.

preprint2015arXiv

Radio halos in SZ-selected clusters of galaxies: the making of a halo?

Radio halos are synchrotron radio sources detected in some massive galaxy clusters. Their Mpc-size indicates that (re)acceleration processes are taking place in the host cluster. X-ray catalogues of galaxy clusters have been used in the past to search for radio halos and to understand their connection with cluster-cluster mergers and with the thermal component of the intra-cluster medium. More recently, the Sunyaev-Zel'dovich effect has been proven to be a better route to search for massive clusters in a wider redshift range. With the aim of discovering new radio halos and understanding their connection with cluster-cluster mergers, we have selected from the Planck Early source catalog the most massive clusters, and we have observed with the Giant Metrewave Radio Telescope at 323 MHz those objects for which deep observations were not available. We have discovered new peculiar radio emission in three of the observed clusters finding: (i) a radio halo in the cluster RXCJ0949.8+1708; (ii) extended emission in Abell 1443 that we classify as a radio halo plus a radio relic, with a bright filament embedded in the radio halo; (iii) low-power radio emission is found in CIZA J1938.3+5409 which is ten times below the radio - X-ray correlation, and represents the first direct detection of the radio emission in the "upper-limit" region of the radio - X-ray diagram. We discuss the properties of these new radio halos in the framework of theoretical models for the radio emission.

preprint2015arXiv

Radio-AGN Feedback: When the Little Ones were Monsters

We present a study of the evolution of the fraction of radio-loud active galactic nuclei (AGN) as a function of their host stellar mass. We make use of two samples of radio galaxies: one in the local universe, $0.01 < z < 0.3$, using a combined SDSS-NVSS sample and one at higher redshifts, $0.5 < z < 2$, constructed from the VLA-COSMOS_DEEP Radio Survey at 1.4 GHz and a K$_s$-selected catalogue of the COSMOS/UltraVISTA field. We observe an increase of more than an order of magnitude in the fraction of lower mass galaxies ($M_* < 10^{10.75}$ M$_{\odot}$) which host Radio-Loud AGN with radio powers $P_{1.4GHz} > 10^{24}$ W/Hz at z ~ 1-2 while the radio-loud fraction for higher mass galaxies ($M_* > 10^{11.25}$ M$_{\odot}$) remains the same. We argue that this increase is driven largely by the increase in cold or radiative mode accretion with increasing cold gas supply at earlier epochs. The increasing population of low mass Radio-Loud AGN can thus explain the upturn in the Radio Luminosity Function (RLF) at high redshift which is important for understanding the impact of AGN feedback in galaxy evolution.

preprint2015arXiv

The peculiar radio galaxy 4C 35.06: a case for recurrent AGN activity?

Using observations obtained with the LOw Fequency ARray (LOFAR), the Westerbork Synthesis Radio Telescope (WSRT) and archival Very Large Array (VLA) data, we have traced the radio emission to large scales in the complex source 4C 35.06 located in the core of the galaxy cluster Abell 407. At higher spatial resolution (~4"), the source was known to have two inner radio lobes spanning 31 kpc and a diffuse, low-brightness extension running parallel to them, offset by about 11 kpc (in projection). At 62 MHz, we detect the radio emission of this structure extending out to 210 kpc. At 1.4 GHz and intermediate spatial resolution (~30"), the structure appears to have a helical morphology. We have derived the characteristics of the radio spectral index across the source. We show that the source morphology is most likely the result of at least two episodes of AGN activity separated by a dormant period of around 35 Myr. The AGN is hosted by one of the galaxies located in the cluster core of Abell 407. We propose that it is intermittently active as it moves in the dense environment in the cluster core. Using LOFAR, we can trace the relic plasma from that episode of activity out to greater distances from the core than ever before. Using the the WSRT, we detect HI in absorption against the center of the radio source. The absorption profile is relatively broad (FWHM of 288 km/s), similar to what is found in other clusters. Understanding the duty cycle of the radio emission as well as the triggering mechanism for starting (or restarting) the radio-loud activity can provide important constraints to quantify the impact of AGN feedback on galaxy evolution. The study of these mechanisms at low frequencies using morphological and spectral information promises to bring new important insights in this field.

preprint2014arXiv

A distant radio mini-halo in the Phoenix galaxy cluster

We report the discovery of extended radio emission in the Phoenix cluster (SPT-CL J2344-4243, z=0.596) with the GMRT at 610 MHz. The diffuse emission extends over a region of at least 400-500 kpc and surrounds the central radio source of the Brightest Cluster Galaxy, but does not appear to be directly associated with it. We classify the diffuse emission as a radio mini-halo, making it the currently most distant mini-halo known. Radio mini-halos have been explained by synchrotron emitting particles re-accelerated via turbulence, possibly induced by gas sloshing generated from a minor merger event. Chandra observations show a non-concentric X-ray surface brightness distribution, which is consistent with this sloshing interpretation. The mini-halo has a flux density of $17\pm5$ mJy, resulting in a 1.4 GHz radio power of ($10.4\pm3.5) \times 10^{24}$ W Hz$^{-1}$. The combined cluster emission, which includes the central compact radio source, is also detected in a shallow GMRT 156 MHz observation and together with the 610 MHz data we compute a spectral index of $-0.84\pm0.12$ for the overall cluster radio emission. Given that mini-halos typically have steeper radio spectra than cluster radio galaxies, this spectral index should be taken as an upper limit for the mini-halo.

preprint2014arXiv

A giant radio halo in the cool core cluster CL1821+643

Giant radio halos are Mpc-size sources found in some merging galaxy clusters. The synchrotron emitting electrons are thought to be (re)accelerated by plasma turbulence induced by the merging of two massive clusters. Cool core galaxy clusters have a low temperature core, likely an indication that a major merger has not recently occurred. CL1821+643 is one of the strongest cool core clusters known so far. Surprisingly, we detect a giant radio halo with a largest linear size of $\sim$ 1.1 Mpc. We discuss the radio and X-ray properties of the cluster in the framework of the proposed models for giant radio halos. If a merger is causing the radio emission, despite the presence of a cool-core, we suggest that it should be off-axis, or in an early phase, or a minor one.

preprint2014arXiv

Evidence for particle re-acceleration in the radio relic in the galaxy cluster PLCKG287.0 +32.9

Radio relics are diffuse radio sources observed in galaxy clusters, probably produced by shock acceleration during cluster-cluster mergers. Their large size, of the order of 1 Mpc, indicates that the emitting electrons need to be (re)accelerated locally. The usually invoked Diffusive Shock Acceleration models have been challenged by recent observations and theory. We report the discovery of complex radio emission in the Galaxy cluster PLCKG287.0 +32.9, which hosts two relics, a radio halo, and several radio filamentary emission. Optical observations suggest that the cluster is elongated, likely along an intergalactic filament, and displays a significant amount of substructure. The peculiar features of this radio relic are that (i) it appears to be connected to the lobes of a radio galaxy and (ii) the radio spectrum steepens on either side of the radio relic. We discuss the origins of these features in the context of particle re-acceleration.

preprint2014arXiv

The discovery of a radio halo in PLCK G147.3-16.6 at z=0.65

Recent X-ray and Sunyaev-Zel'dovich (SZ) observations have dramatically increased the number of known distant galaxy clusters. In some merging, low-redshift (z<0.4) clusters, centrally-located, diffuse, extended radio emission (called radio halos) has been found. Using the Giant Metrewave Radio Telescope (GMRT), we report the detection of diffuse radio emission in the binary-merging cluster PLCK G147.3-16.6 located at z = 0.65. We classify the emission as a giant radio halo due to the large physical extent of about 0.9 Mpc and low-surface brightness. We measure an integrated flux density of $7.3\pm1.1$ mJy at 610 MHz for the radio halo, resulting in a 1.4 GHz radio power of $5.1 \times 10^{24}$ W/Hz. The radio halo power is consistent with that expected from the known correlation between X-ray luminosity or the cluster integrated SZ signal and radio power. Our observations also suggest that more of these distant radio halos could be discovered with the GMRT.

preprint2014arXiv

The nature of the low-frequency emission of M51: First observations of a nearby galaxy with LOFAR

The grand-design spiral galaxy M51 was observed with the LOFAR High Frequency Antennas (HBA) and imaged in total intensity and polarisation. This observation covered the frequencies between 115 MHz and 175 MHz. We produced an image of total emission of M51 at the mean frequency of 151 MHz with 20 arcsec resolution and 0.3 mJy rms noise, which is the most sensitive image of a galaxy at frequencies below 300 MHz so far. The integrated spectrum of total radio emission is described well by a power law, while flat spectral indices in the central region indicate thermal absorption. We observe that the disk extends out to 16 kpc and see a break in the radial profile near the optical radius of the disk. Our main results, the scale lengths of the inner and outer disks at 151 MHz and 1.4 GHz, arm--interarm contrast, and the break scales of the radio--far-infrared correlations, can be explained consistently by CRE diffusion, leading to a longer propagation length of CRE of lower energy. The distribution of CRE sources drops sharply at about 10 kpc radius, where the star formation rate also decreases sharply. We find evidence that thermal absorption is primarily caused by HII regions. The non-detection of polarisation from M51 at 151 MHz is consistent with the estimates of Faraday depolarisation. Future searches for polarised emission in this frequency range should concentrate on regions with low star formation rates.

preprint2013arXiv

Discovery of spectral curvature in the shock downstream region: CIZA J2242.8+5301

Giant cluster radio relics are thought to form at shock fronts in the course of collisions between galaxy clusters. Via processes that are still poorly understood, these shocks accelerate or re-accelerate cosmic-ray electrons and might amplify magnetic fields. The best object to study this phenomenon is the galaxy cluster CIZA J2242.8+5301 as it shows the most undisturbed relic. By means of Giant Metrewave Radio Telescope (GMRT) and Westerbork Synthesis Radio Telescope (WSRT) data at seven frequencies spanning from 153 MHz to 2272 MHz, we study the synchrotron emission in this cluster. We aim at distinguishing between theoretical injection and acceleration models proposed for the formation of radio relics. We also study the head-tail radio sources to reveal the interplay between the merger and the cluster galaxies. We produced spectral index, curvature maps and radio colour-colour plots and compared our data with predictions from models. We present one of the deepest 153 MHz maps of a cluster ever produced, reaching a noise level of 1.5 mJy/beam. We derive integrated spectra for four relics in the cluster, discovering extremely steep spectrum diffuse emission concentrated in multiple patches. We find a possible radio phoenix embedded in the relic to the south of the cluster. The spectral index of the northern relic retains signs of steepening from the front towards the back of the shock also at the radio frequencies below 600 MHz. The spectral curvature in the same relic also increases in the downstream area. The data is consistent with the Komissarov-Gubanov injection models, meaning that the emission we observe is produced by a single burst of spectrally-aged accelerated radio electrons.

preprint2013arXiv

Searching for large-scale structures around high-redshift radio galaxies with Herschel

This paper presents the first results of a far-infrared search for protocluster-associated galaxy overdensities using the SPIRE instrument on-board the {\it Herschel} Space Observatory. Large ($\sim$400 arcmin$^{2}$) fields surrounding 26 powerful high-redshift radio galaxies ($2.0 < z < 4.1$; $L_{\rm 500 MHz} > 10^{28.5}$ WHz$^{-1}$) are mapped at 250, 350 and 500\mic to give a unique wide-field sample. On average the fields have a higher than expected, compared to blank fields, surface density of 500\mic sources within 6 comoving Mpc of the radio galaxy. The analysis is then restricted to potential protocluster members only, which are identified using a far-infrared colour selection; this reveals significant overdensities of galaxies in 2 fields, neither of which are previously known protoclusters. The probability of finding 2 overdensities of this size by chance, given the number of fields observed is $5 \times 10^{-4}$. Overdensities here exist around radio galaxies with $L_{\rm 500 MHz} \gtrsim 10^{29}$ WHz$^{-1}$ and $z < 3$. The radial extent of the average far-infrared overdensity is found to be $\sim$6 comoving Mpc. Comparison with predictions from numerical simulations shows that the overdensities are consistent with having masses $> 10^{14}$Msolar. However, the large uncertainty in the redshift estimation means that it is possible that these far-infrared overdensities consist of several structures across the redshift range searched.

preprint2013arXiv

The host galaxy of the z=2.4 radio-loud AGN MRC 0406-244 as seen by HST

We present multicolour Hubble Space Telescope images of the powerful z=2.4 radio galaxy MRC 0406-244 and model its complex morphology with several components including a host galaxy, a point source, and extended nebular and continuum emission. We suggest that the main progenitor of this radio galaxy was a normal, albeit massive (M ~10^{11} solar masses), star-forming galaxy. The optical stellar disc of the host galaxy is smooth and well described by a Sérsic profile, which argues against a recent major merger, however there is also a point-source component which may be the remnant of a minor merger. The half-light radius of the optical disc is constrained to lie in the range 3.5 to 8.2kpc, which is of similar size to coeval star forming galaxies. Biconical shells of nebular emission and UV-bright continuum extend out from the host galaxy along the radio jet axis, which is also the minor axis of the host galaxy. The origin of the continuum emission is uncertain, but it is most likely to be young stars or dust-scattered light from the AGN, and it is possible that stars are forming from this material at a rate of 200^{+1420}_{-110} solar masses per year.

preprint2012arXiv

A Deep Chandra Observation of the AGN Outburst and Merger in Hickson Compact Group 62

We report on an analysis of new Chandra data of the galaxy group HCG 62, well known for possessing cavities in its intragroup medium (IGM) that were inflated by the radio lobes of its central active galactic nucleus (AGN). With the new data, a factor of three deeper than previous Chandra data, we re-examine the energetics of the cavities and determine new constraints on their contents. We confirm that the ratio of radiative to mechanical power of the AGN outburst that created the cavities is less than 10^-4, among the lowest of any known cavity system, implying that the relativistic electrons in the lobes can supply only a tiny fraction of the pressure required to support the cavities. This finding implies additional pressure support in the lobes from heavy particles (e.g., protons) or thermal gas. Using spectral fits to emission in the cavities, we constrain any such volume-filling thermal gas to have a temperature kT > 4.3 keV. For the first time, we detect X-ray emission from the central AGN, with a luminosity of L(2-10 keV) = (1.1 +/- 0.4) x 10^39 erg s^-1 and properties typical of a low-luminosity AGN. Lastly, we report evidence for a recent merger from the surface brightness, temperature, and metallicity structure of the IGM.

preprint2012arXiv

T-RaMiSu: The Two-meter Radio Mini Survey I. The Bootes Field

We present wide area, deep, high-resolution 153 MHz GMRT observations of the NOAO Bootes field, adding to the extensive, multi-wavelength data of this region. The observations, data reduction, and catalogue construction and description are described here. The seven pointings produced a final mosaic covering 30 square degrees with a resolution of 25". The rms noise is 2 mJy/beam in the centre of the image, rising to 4-5 mJy/beam on the edges, with an average of 3 mJy/beam. Seventy-five per cent of the area has an rms < 4 mJy/beam. The extracted source catalogue contains 1289 sources detected at 5σ, of which 453 are resolved. We estimate the catalogue to be 92 per cent reliable and 95 per cent complete at an integrated flux density limit of 14 mJy. The flux densities and astrometry have been corrected for systematic errors. We calculate the differential source counts {which are in good agreement with those in the literature and provide an important step forward in quantifying the source counts at these low frequencies and low flux densities}. The GMRT 153 MHz sources have been matched to the 1.4 GHz NVSS and 327 MHz WENSS catalogues and spectral indices were derived.

preprint2012arXiv

The duty cycle of radio-mode feedback in complete samples of clusters

The Chandra X-ray Observatory has revealed X-ray bubbles in the intracluster medium (ICM) of many nearby cooling flow clusters. The bubbles trace feedback that is thought to couple the central active galactic nucleus (AGN) to the ICM, helping to stabilize cooling flows and govern the evolution of massive galaxies. However, the prevalence and duty cycle of such AGN outbursts is not well understood. To this end, we study how cooling is balanced by bubble heating for complete samples of clusters (the Brightest 55 clusters of galaxies, hereafter B55, and the HIghest X-ray FLUx Galaxy Cluster Sample, HIFLUGCS). We find that the radio luminosity of the central galaxy only exceeds 2.5 x 10^30 erg s^-1 Hz^-1 in cooling flow clusters. This result implies a connection between the central radio source and the ICM, as expected if AGN feedback is operating. Additionally, we find a duty cycle for radio mode feedback, the fraction of time that a system possesses bubbles inflated by its central radio source, of > 69 per cent for B55 and > 63 per cent for HIFLUGCS. These duty cycles are lower limits since some bubbles are likely missed in existing images. We used simulations to constrain the bubble power that might be present and remain undetected in the cooling flow systems without detected bubbles. Among theses systems, almost all could have significant bubble power. Therefore, our results imply that the duty cycle of AGN outbursts with the potential to heat the gas significantly in cooling flow clusters is at least 60 per cent and could approach 100 per cent.

preprint2012arXiv

The faint source population at 15.7 GHz - I. The radio properties

We have studied a sample of 296 faint (> 0.5 mJy) radio sources selected from an area of the Tenth Cambridge (10C) survey at 15.7 GHz in the Lockman Hole. By matching this catalogue to several lower frequency surveys (e.g. including a deep GMRT survey at 610 MHz, a WSRT survey at 1.4 GHz, NVSS, FIRST and WENSS) we have investigated the radio spectral properties of the sources in this sample; all but 30 of the 10C sources are matched to one or more of these surveys. We have found a significant increase in the proportion of flat spectrum sources at flux densities below approximately 1 mJy - the median spectral index between 15.7 GHz and 610 MHz changes from 0.75 for flux densities greater than 1.5 mJy to 0.08 for flux densities less than 0.8 mJy. This suggests that a population of faint, flat spectrum sources is emerging at flux densities below 1 mJy. The spectral index distribution of this sample of sources selected at 15.7 GHz is compared to those of two samples selected at 1.4 GHz from FIRST and NVSS. We find that there is a significant flat spectrum population present in the 10C sample which is missing from the samples selected at 1.4 GHz. The 10C sample is compared to a sample of sources selected from the SKADS Simulated Sky by Wilman et al. and we find that this simulation fails to reproduce the observed spectral index distribution and significantly underpredicts the number of sources in the faintest flux density bin. It is likely that the observed faint, flat spectrum sources are a result of the cores of FRI sources becoming dominant at high frequencies. These results highlight the importance of studying this faint, high frequency population.

preprint2012arXiv

The Herschel view of the environment of the radio galaxy 4C+41.17 at z = 3.8

We present Herschel observations at 70, 160, 250, 350 and 500 micron of the environment of the radio galaxy 4C+41.17 at z = 3.792. About 65% of the extracted sources are securely identified with mid-IR sources observed with the Spitzer Space Telescope at 3.6, 4.5, 5.8, 8 and 24 micron. We derive simple photometric redshifts, also including existing 850 micron and 1200 micron data, using templates of AGN, starburst-dominated systems and evolved stellar populations. We find that most of the Herschel sources are foreground to the radio galaxy and therefore do not belong to a structure associated with 4C+41.17. We do, however, find that the SED of the closest (~ 25" offset) source to the radio galaxy is fully consistent with being at the same redshift as 4C+41.17. We show that finding such a bright source that close to the radio galaxy at the same redshift is a very unlikely event, making the environment of 4C+41.17 a special case. We demonstrate that multi-wavelength data, in particular on the Rayleigh-Jeans side of the spectral energy distribution, allow us to confirm or rule out the presence of protocluster candidates that were previously selected by single wavelength data sets.

preprint2011arXiv

H-alpha emitters in z~2 proto-clusters: evidence for faster evolution in dense environments

This is a study of H-alpha emitters in two dense galaxy proto-clusters surrounding radio galaxies at z~2. We show that the proto-cluster surrounding MRC 1138-262 contains 14+/-2 times more H-alpha candidates than the average field (9 sigma significance), and the z=2.35 radio galaxy 4C+10.48 is surrounded by 12+/-2 times more emitters than the field (5 sigma), so it is also likely to reside in a dense proto-cluster environment. We compared these H-alpha emitters, situated in dense environments, to a control field sample selected from 3 separate fields forming a total area of 172 arcmin^2. We constructed and compared H-alpha and rest-frame R continuum luminosity functions of the emitters in both environments. The star formation density is on average 13 times greater in the proto-clusters than the field at z~2, so the total star formation rate within the central 1.5Mpc of the proto-clusters exceeds 3000Msun/yr. However, we found no significant difference in the shape of the H-alpha luminosity functions, implying that environment does not substantially affect the strength of the H-alpha line from strongly star forming galaxies. The proto-cluster emitters are typically 0.8mag brighter in rest-frame R continuum than field emitters, implying they are twice as massive as their field counterparts at the same redshift. We also show the proto-cluster galaxies have lower specific star formation rates than field galaxies, meaning the emitters in the dense environments formed more of their stars earlier than the field galaxies. We conclude that galaxy growth in the early Universe was accelerated in dense environments, and that cluster galaxies differed from field galaxies even before the cluster had fully formed.

preprint2011arXiv

New Insights on the z-alpha Correlation from Complete Radio Samples

The existence of a correlation between observed radio spectral index (alpha) and redshift (z) has long been used as a method for locating high-z radio galaxies. We use 9 highly spectroscopically complete radio samples, selected at different frequencies and flux limits, to determine the efficiency of this method, and compare consistently observed correlations between alpha, luminosity, linear size, and redshift. We observe a weak correlation between z and alpha which remains even when Malmquist bias is removed. The strength of this correlation depends on both the k-correction and sample selection frequency, in addition to the frequency at which alpha is measured, and consistent results for both high and low frequency selected samples are only seen if analysis is restricted to just extended radio galaxies. Many of the highest redshift radio galaxies are very compact and often display a negatively curved or peaked spectrum, and therefore the low-frequency radio spectrum as a whole should be studied; this is something for which the LOFAR will be crucial. We quantify both the efficiency and the completeness of various techniques used to select high-z radio galaxies. A steep-spectrum cut applied to low-frequency selected samples can more than double the fraction of high-z sources, but at a cost of excluding over half of the high-z sources present in the original sample. An angular size cut is an almost as equally effective method as a steep-spectrum cut, and works for both high and low frequency selected samples. In multi-wavelength data, selection first of infrared-faint radio sources remains by far the most efficient method of selecting high-z sources. We present a simple method for selecting high-z radio galaxies, based purely on combining their radio properties of alpha and angular size, with the addition of the K-band magnitude if available.[abridged]

preprint2011arXiv

The luminosity-dependent high-redshift turnover in the steep spectrum radio luminosity function: clear evidence for downsizing in the radio-AGN population

This paper presents a new grid-based method for investigating the evolution of the steep-spectrum radio luminosity function, with the aim of quantifying the high-redshift cut-off suggested by previous work. To achieve this, the Combined EIS-NVSS Survey of Radio Sources (CENSORS) has been developed; this is a 1.4 GHz radio survey, containing 135 sources complete to a flux density of 7.2 mJy, selected from the NRAO VLA Sky Survey (NVSS) over 6 deg^2 of the ESO Imaging Survey (EIS) Patch D. The sample is currently 73% spectroscopically complete, with the remaining redshifts estimated via the K-z or I-z magnitude-redshift relation. CENSORS is combined with additional radio data from the Parkes All-Sky, Parkes Selected Regions, Hercules and VLA COSMOS samples to provide comprehensive coverage of the radio power vs. redshift plane. The redshift distributions of these samples, together with radio source count determinations, and measurements of the local luminosity function, provide the input to the fitting process. The modelling reveals clear declines, at > 3sigma significance, in comoving density at z > 0.7 for lower luminosity sources (log P = 25-26); these turnovers are still present at log P > 27, but move to z > 3, suggesting a luminosity-dependent evolution of the redshift turnover, similar to the `cosmic downsizing' seen for other AGN populations. These results are shown to be robust to the estimated redshift errors and to increases in the spectral index for the highest redshift sources. Analytic fits to the best-fitting steep spectrum grid are provided so that the results presented here can be easily accessed by the reader, as well as allowing plausible extrapolations outside of the regions covered by the input datasets

preprint2010arXiv

Galaxy protocluster candidates around z ~ 2.4 radio galaxies

We study the environments of 6 radio galaxies at 2.2 < z < 2.6 using wide-field near-infrared images. We use colour cuts to identify galaxies in this redshift range, and find that three of the radio galaxies are surrounded by significant surface overdensities of such galaxies. The excess galaxies that comprise these overdensities are strongly clustered, suggesting they are physically associated. The colour distribution of the galaxies responsible for the overdensity are consistent with those of galaxies that lie within a narrow redshift range at z ~ 2.4. Thus the excess galaxies are consistent with being companions of the radio galaxies. The overdensities have estimated masses in excess of 10^14 solar masses, and are dense enough to collapse into virizalised structures by the present day: these structures may evolve into groups or clusters of galaxies. A flux-limited sample of protocluster galaxies with K < 20.6 mag is derived by statistically subtracting the fore- and background galaxies. The colour distribution of the protocluster galaxies is bimodal, consisting of a dominant blue sequence, comprising 77 +/- 10% of the galaxies, and a poorly populated red sequence. The blue protocluster galaxies have similar colours to local star-forming irregular galaxies (U -V ~ 0.6), suggesting most protocluster galaxies are still forming stars at the observed epoch. The blue colours and lack of a dominant protocluster red sequence implies that these cluster galaxies form the bulk of their stars at z < 3.

preprint2009arXiv

Probing Fine-Scale Ionospheric Structure with the Very Large Array Radio Telescope

High resolution (~1 arcminute) astronomical imaging at low frequency (below 150 MHz) has only recently become practical with the development of new calibration algorithms for removing ionospheric distortions. In addition to opening a new window in observational astronomy, the process of calibrating the ionospheric distortions also probes ionospheric structure in an unprecedented way. Here we explore one aspect of this new type of ionospheric measurement, the differential refraction of celestial source pairs as a function of their angular separation. This measurement probes variations in the spatial gradient of the line-of-sight total electron content (TEC) to 0.001 TECU/km accuracy over spatial scales of under 10 km to over 100 km. We use data from the VLA Low-frequency Sky Survey (VLSS; Cohen et al. 2007, AJ 134, 1245), a nearly complete 74 MHz survey of the entire sky visible to the Very Large Array (VLA) telescope in Socorro, New Mexico. These data comprise over 500 hours of observations, all calibrated in a standard way. While ionospheric spatial structure varies greatly from one observation to the next, when analyzed over hundreds of hours, statistical patterns become apparent. We present a detailed characterization of how the median differential refraction depends on source pair separation, elevation and time of day. We find that elevation effects are large, but geometrically predictable and can be "removed" analytically using a "thin-shell" model of the ionosphere. We find significantly greater ionospheric spatial variations during the day than at night. These diurnal variations appear to affect the larger angular scales to a greater degree indicating that they come from disturbances on relatively larger spatial scales (100s of km, rather than 10s of km).

preprint2007arXiv

The XMM-LSS catalogue: X-ray sources and associated optical data. Version I

Following the presentation of the XMM-LSS X-ray source detection package by Pacaud et al., we provide the source lists for the first 5.5 surveyed square degrees. The catalogues pertain to the [0.5-2] and [2-10] keV bands and contain in total 3385 point-like or extended sources above a detection likelihood of 15 in either band. The agreement with deep logN-logS is excellent. The main parameters considered are position, countrate, source extent with associated likelihood values. A set of additional quantities such as astrometric corrections and fluxes are further calculated while errors on the position and countrate are deduced from simulations. We describe the construction of the band-merged catalogue allowing rapid sub-sample selection and easy cross-correlation with external multi-wavelength catalogues. A small optical CFHTLS multi-band subset of objects is associated wich each source along with an X-ray/optical overlay. We make the full X-ray images available in FITS format. The data are available at CDS and, in a more extended form, at the Milan XMM-LSS database.

preprint2007arXiv

The XMM-LSS survey: the Class 1 cluster sample over the initial 5 square degrees and its cosmological modelling

We present a sample of 29 galaxy clusters from the XMM-LSS survey over an area of some 5deg2 out to a redshift of z=1.05. The sample clusters, which represent about half of the X-ray clusters identified in the region, follow well defined X-ray selection criteria and are all spectroscopically confirmed. For all clusters, we provide X-ray luminosities and temperatures as well as masses. The cluster distribution peaks around z=0.3 and T =1.5 keV, half of the objects being groups with a temperature below 2 keV. Our L-T(z) relation points toward self-similar evolution, but does not exclude other physically plausible models. Assuming that cluster scaling laws follow self-similar evolution, our number density estimates up to z=1 are compatible with the predictions of the concordance cosmology and with the findings of previous ROSAT surveys. Our well monitored selection function allowed us to demonstrate that the inclusion of selection effects is essential for the correct determination of the evolution of the L-T relation, which may explain the contradictory results from previous studies. Extensive simulations show that extending the survey area to 10deg2 has the potential to exclude the non-evolution hypothesis, but that constraints on more refined ICM models will probably be limited by the large intrinsic dispersion of the L-T relation. We further demonstrate that increasing the dispersion in the scaling laws increases the number of detectable clusters, hence generating further degeneracy [in addition to sigma8, Omega_m, L(M,z) and T(M,z)] in the cosmological interpretation of the cluster number counts. We provide useful empirical formulae for the cluster mass-flux and mass-count-rate relations as well as a comparison between the XMM-LSS mass sensitivity and that of forthcoming SZ surveys.

preprint1997arXiv

A jet-cloud interaction in 3C34 at redshift z = 0.69

We report the detection of a strong jet-cloud interaction at a distance of 120 kpc from the nucleus of the radio galaxy 3C34, which has redshift z=0.69. Hubble Space Telescope images of the radio galaxy show a long narrow region of blue emission orientated along the radio axis and directed towards a radio hotspot. The William Herschel Telescope has been used to provide long-slit spectroscopic data of this object, and infrared observations made with the United Kingdom InfraRed Telescope have enabled its spectral energy distribution to be modelled. We propose that the aligned emission is associated with a region of massive star-formation, induced by the passage of the radio jet through a galaxy within the cluster surrounding 3C34. A star-formation rate of about 100 solar masses per year is required, similar to the values necessary to produce the alignment effect in high-redshift radio galaxies. The consequences of this result for models of star formation in distant radio galaxies are discussed.

preprint1996arXiv

A search for molecular gas in high redshift radio galaxies

We present results of an extensive search for molecular gas in 14 high redshift ($z>2$) radio galaxies. Radio galaxies with redshifts between 2.0 and 4.3 were observed during several sessions on two different single-dish telescopes (IRAM and JCMT) and two interferometric arrays (VLA and OVRO). No significant CO emission to a 2$σ$ limit of a few times 10^{10} K km s^{-1} pc^2 was detected in any of these objects. The limits on the CO emission achieved indicate that, assuming a CO-H$_2$ conversion factor similar to the Galactic value, the mass of enriched molecular gas in these galaxies is less than $10^{11} $M$_{\odot}$. This suggests that distant radio galaxies may not be forming stars at the extremely high rates that have earlier been proposed, although they may still be forming stars at rates comparable to those of local starburst galaxies.

preprint1996arXiv

WN1626+5153: A Giant Radio Galaxy from the WENSS survey

An FRII radio source with an angular extension of 18.4 arcminutes has been discovered in the Westerbork Northern Sky Survey that is being carried out at 325 MHz. This source is identified with Mrk 1498 at $z=0.056$ indicating that its projected extension is 1.6 \h50\ Mpc. The optical spectrum of Mrk 1498 has the typical emission lines of a narrow line radio galaxy; H$α$ is the only permitted line that has a faint broad component. The discovery shows that the good sensitivity of WENSS at low frequencies to low surface brightness features at scales of a few minutes, is uniquely suited to provide well defined and relatively large samples of giant radio galaxies. The future definition of such samples will be important for our understanding of giant radio galaxies.