Source author record

R. J. van Weeren

R. J. van Weeren 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

135works
4topics
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

135 published item(s)

preprint2026arXiv

Testing the correlation between bending angle and polarization properties of bent radio galaxies

The bending of radio galaxies in galaxy clusters is expected to be caused by interactions with the local environment. The physical processes responsible for jet bending, and their influence on the polarization properties of radio galaxies, remain poorly understood, leading to the question of whether jet properties in bent radio galaxies differ from those in linear radio galaxies. Using a sample of 24 polarized bent radio galaxies, observed with the Karl G. Jansky Very Large Array at 1--2 GHz, we test for correlation of bending angle with polarization parameters measuring Faraday rotation, intrinsic fractional polarization, and Faraday rotation dispersion, used here as a measure of turbulence along the line of sight. We find no statistically significant correlations. At the spatial resolution of our dataset (3--46 kpc, median 18.4 kpc), our results indicate that we are primarily probing larger-scale intracluster medium effects not related to bending angle. The absence of a statistically significant correlation suggests that bent radio galaxies are reliable probes of intracluster magnetic fields, because their intrinsic properties do not appear to introduce systematic biases into measured polarization parameters. We do detect a preference for source magnetic field vectors to align with the direction of jet bending. Finally, we estimate that the POSSUM and SKA surveys will contain $\gtrsim$300 and $\gtrsim$1000 polarized radio galaxies, respectively, providing large future samples with a range of bending angles and similar redshift distribution and number of beams per source as in our sample, enabling our results to be tested with greater statistical power.

preprint2025arXiv

LOFAR 58 MHz Legacy Survey of the 3CRR Catalog

The Low Frequency Array (LOFAR) is uniquely able to perform deep, 15" resolutions imaging at frequencies below 100 MHz. Observations in this regime, using the Low Band Antenna (LBA) system, are significantly affected by instrumental and ionospheric distortions. Recent developments in calibration techniques have enabled routine production of high-fidelity images at these challenging frequencies. The aim of this paper is to obtain images of the radio sources included in the Third Cambridge catalog, second revised version (3CRR), at an observing frequency of 58 MHz, with an angular resolution of 15"and sensitivity to both compact and diffuse radio emission. This work also aims to produce accurate flux measurements for all sources. This dataset is designed to serve as a reference for low-frequency radio galaxy studies and future spectral aging analyses. We deliver 58. MHz radio images for the complete 3CRR sample including flux density measurements. We determined that the LBA has an accurate flux density scale with an average flux uncertainty of 10%. This is an important confirmation for any future works using the LOFAR LBA system. With these results we characterize the bright radio galaxy population with new high-resolution low-frequency images. We also provide high-resolution models of these sources which will be useful for calibrating future surveys. This legacy survey significantly expands the available high-resolution data at low frequencies and is the first fully imaged high-resolution sample at ultra low frequencies (< 100 MHz). It lays the foundation for future studies of radio galaxy physics, low-energy cosmic-ray populations, and the interplay between radio jets and their environments.

preprint2025arXiv

Resurging from the ashes: A spectral study of seven candidate revived radio fossils in nearby low-mass galaxy clusters

Complex energy transfer processes in the intracluster medium (ICM) can revive fossil (with spectral ages $\gg100$ Myr) plasma initially generated by radio galaxies. This leads to the re-ignition of faint radio sources with irregular and filamentary morphologies, and ultra-steep ($α\gtrsim 1.5$) synchrotron spectra, which can be more easily detected at low frequencies ($\sim 100$ MHz). These sources offer the opportunity to investigate the microphysics of the ICM and its interplay with radio galaxies, the origin of seed relativistic electrons, the merging history of the host cluster, and the phenomenology of radio filaments. The study of revived sources has so far been hampered by the requirement of sensitive and high-resolution multi-frequency radio data at low frequencies to characterise their spatial properties and provide a proper classification. We aim to perform the analysis of a sample of candidate revived sources identified among nearby ($z\leq0.35$) and low-mass ($M_{500}\leq5\times 10^{14} M_\odot$) \textit{Planck} clusters in the footprint of LoTSS-DR2. By inspecting LoTSS-DR2 images at 144 MHz, we identified 7 targets with patchy and filamentary morphologies, which have been followed-up with the uGMRT at 400 MHz. By combining LOFAR and uGMRT data, we obtained high-resolution images and spectral index maps, which we used to interpret the nature of the sources. All targets show regions with very steep spectra, confirming the effectiveness of our morphology-based selection in identifying fossil plasma. Based on their morphology, spectral properties, and optical associations, we investigated the origin of the targets. We found a variety of promising revived fossil sources, while also showing that apparently intricate structures can be easily misclassified in the absence of high-resolution and multi-band data.

preprint2023arXiv

Extended radio emission in the galaxy cluster MS 0735.6+7421 detected with the Karl G. Jansky Very Large Array

MS 0735.6+7421 ($z = 0.216$) is a massive cool core galaxy cluster hosting one of the most powerful active galactic nuclei (AGN) outbursts known. The radio jets of the AGN have carved out an unusually large pair of X-ray cavities, each reaching a diameter of $200$ kpc. This makes MS 0735.6+7421 a unique case to investigate active galactic nuclei feedback processes, as well as other cluster astrophysics at radio wavelengths. We present new low-radio-frequency observations of MS 0735.6+7421 taken with the Karl G. Jansky Very Large Array (VLA): 5 hours of P-band ($224-480$ MHz) and 5 hours of L-band ($1-2$ GHz) observations, both in C configuration. Our VLA P-band ($224-480$ MHz) observations reveal the presence of a new diffuse radio component reaching a scale of $\sim$ $900$ kpc in the direction of the jets and of $\sim$ $500$ kpc in the direction perpendicular to the jets. This component is centered on the cluster core and has a radio power scaled at $1.4$ GHz of $P_{1.4\text{ GHz}} = (4\pm2)\times 10^{24}$ WHz$^{-1}$. Its properties are consistent with those expected from a radio mini-halo as seen in other massive cool core clusters, although it may also be associated with radio plasma that has diffused out of the X-ray cavities. Observations at higher spatial resolution are needed to fully characterize the properties and nature of this component. We also suggest that if radio mini-halos originate from jetted activity, we may be witnessing the early stages of this process.

preprint2023arXiv

The Faraday Rotation Measure Grid of the LOFAR Two-metre Sky Survey: Data Release 2

A Faraday rotation measure (RM) catalogue, or RM Grid, is a valuable resource for the study of cosmic magnetism. Using the second data release (DR2) from the LOFAR Two-metre Sky Survey (LoTSS), we have produced a catalogue of 2461 extragalactic high-precision RM values across 5720 deg$^{2}$ of sky (corresponding to a polarized source areal number density of $\sim$0.43 deg$^{-2}$). The linear polarization and RM properties were derived using RM synthesis from the Stokes $Q$ and $U$ channel images at an angular resolution of 20'' across a frequency range of 120 to 168 MHz with a channel bandwidth of 97.6 kHz. The fraction of total intensity sources ($>1$ mJy beam$^{-1}$) found to be polarized was $\sim$0.2%. The median detection threshold was 0.6 mJy beam$^{-1}$ ($8σ_{QU}$), with a median RM uncertainty of 0.06 rad m$^{-2}$ (although a systematic uncertainty of up to 0.3 rad m$^{-2}$ is possible, after the ionosphere RM correction). The median degree of polarization of the detected sources is 1.8%, with a range of 0.05% to 31%. Comparisons with cm-wavelength RMs indicate minimal amounts of Faraday complexity in the LoTSS detections, making them ideal sources for RM Grid studies. Host galaxy identifications were obtained for 88% of the sources, along with redshifts for 79% (both photometric and spectroscopic), with the median redshift being 0.6. The focus of the current catalogue was on reliability rather than completeness, and we expect future versions of the LoTSS RM Grid to have a higher areal number density. In addition, 25 pulsars were identified, mainly through their high degrees of linear polarization.

preprint2022arXiv

A MeerKAT-meets-LOFAR Study of MS 1455.0+2232: A 590 kiloparsec 'Mini'-Halo in a Sloshing Cool-Core Cluster

Radio mini-haloes are poorly-understood, moderately-extended diffuse radio sources that trace the presence of magnetic fields and relativistic electrons on scales of hundreds of kiloparsecs, predominantly in relaxed clusters. With relatively few confirmed detections to-date, many questions remain unanswered. This paper presents new radio observations of the galaxy cluster MS1455.0$+$2232 performed with MeerKAT (covering the frequency range 872$-$1712 MHz) and LOFAR (covering 120$-$168 MHz), the first results from a homogeneously selected mini-halo census. We find that this mini-halo extends for $\sim590$ kpc at 1283 MHz, significantly larger than previously believed, and has a flatter spectral index ($α= -0.97 \pm 0.05$) than typically expected. Our X-ray analysis clearly reveals a large-scale (254 kpc) sloshing spiral in the intracluster medium. We perform a point-to-point analysis, finding a tight single correlation between radio and X-ray surface brightness with a super-linear slope of $b_{\rm 1283~MHz} = 1.16^{+0.06}_{-0.07}$ and $b_{\rm 145~MHz} = 1.15^{+0.09}_{-0.08}$; this indicates a strong link between the thermal and non-thermal components of the intracluster medium. Conversely, in the spectral index/X-ray surface brightness plane, we find that regions inside and outside the sloshing spiral follow different correlations. We find compelling evidence for multiple sub-components in this mini-halo for the first time. While both the turbulent (re-)acceleration and hadronic scenarios are able to explain some observed properties of the mini-halo in MS1455.0$+$2232, neither scenario is able to account for all the evidence presented by our analysis.

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

High-resolution, High-sensitivity, Low-frequency uGMRT View of Coma Cluster of Galaxies

We present high-resolution, high-sensitivity upgraded Giant Metrewave Radio Telescope observations of the Coma cluster (A1656) at 250-500 MHz and 550-850 MHz. At 250-500 MHz, 135 sources have extensions $>$ 0.45 arcmin (with peak-to-local-noise ratio $> 4$). Of these, 24 sources are associated with Coma-member galaxies. In addition, we supplement this sample of 24 galaxies with 20 ram pressure stripped galaxies from Chen et al. (2020, eight are included in the original extended radio source sample) and an additional five are detected and extended. We present radio morphologies, radio spectra, spectral index maps, and equipartition properties for these two samples. In general, we find the equipartition properties lie within a narrow range (e.g., $P_{\rm min}$ = 1-3 dynes cm$^{-2}$). Only NGC 4874, one of the two brightest central Coma cluster galaxies, has a central energy density and pressure about five times higher and a radio source age about 50 % lower than that of the other Coma galaxies. We find a diffuse tail of radio emission trailing the dominant galaxy of the merging NGC 4839 group that coincides with the "slingshot" tail, seen in X-rays. The southwestern radio relic, B1253$+$275, has a large extent $\approx$ 32$^\prime$ $\times$ 10$^\prime$ ($\simeq$ 1.08 $\times$ 0.34 Mpc$^2$). For NGC 4789, whose long radio tails merge into the relic and may be a source of its relativistic seed electrons, and we find a transverse radio spectral gradient, a steepening from southwest to northeast across the width of the radio source. Finally, radio morphologies of the extended and RPS samples suggest that these galaxies are on their first infall into Coma on (predominantly) radial orbits.

preprint2022arXiv

Identifying active galactic nuclei via brightness temperature with sub-arcsecond International LOFAR Telescope observations

Identifying active galactic nuclei (AGN) and isolating their contribution to a galaxy's energy budget is crucial for studying the co-evolution of AGN and their host galaxies. Brightness temperature ($T_b$) measurements from high-resolution radio observations at GHz frequencies are widely used to identify AGN. Here we investigate using new sub-arcsecond imaging at 144 MHz with the International LOFAR Telescope to identify AGN using $T_b$ in the Lockman Hole field. We use ancillary data to validate the 940 AGN identifications, finding 83 percent of sources have AGN classifications from SED fitting and/or photometric identifications, yielding 160 new AGN identifications. Considering the multi-wavelength classifications, brightness temperature criteria select over half of radio-excess sources, 32 percent of sources classified as radio-quiet AGN, and 20 percent of sources classified as star-forming galaxies. Infrared colour-colour plots and comparison with what we would expect to detect based on peak brightness in 6 arcsec LOFAR maps, imply that the star-forming galaxies and sources at low flux densities have a mixture of star-formation and AGN activity. We separate the radio emission from star-formation and AGN in unresolved, $T_b$-identified AGN with no significant radio excess and find the AGN comprises $0.49\pm 0.16$ of the radio luminosity. Overall the non-radio excess AGN show evidence for having a variety of different radio emission mechanisms, which can provide different pathways for AGN and galaxy co-evolution. This validation of AGN identification using brightness temperature at low frequencies opens the possibility for securely selecting AGN samples where ancillary data is inadequate.

preprint2022arXiv

Particle re-acceleration and diffuse radio sources in the galaxy cluster Abell 1550

We study diffuse radio emission in the galaxy cluster A1550, with the aim of constraining particle re-acceleration in the intra-cluster medium. We exploit observations at four different frequencies: 54, 144, 400 and 1400 MHz. To complement our analysis, we make use of archival Chandra X-ray data. At all frequencies we detect an ultra-steep spectrum radio halo ($S_ν\propto ν^{-1.6}$) with an extent of 1.2 Mpc at 54 MHz. Its morphology follows the distribution of the thermal intra-cluster medium inferred from the Chandra observation. West of the centrally located head-tail radio galaxy, we detect a radio relic with projected extent of 500 kpc. From the relic, a 600 kpc long bridge departs and connect it to the halo. Between the relic and the radio galaxy, we observe what is most likely a radio phoenix, given its curved spectrum. The phoenix is connected to the tail of the radio galaxy through two arms, which show a nearly constant spectral index for 300 kpc. The halo could be produced by turbulence induced by a major merger, with its axis lying in the NE-SW direction. This is supported by the position of the relic, whose origin could be attributed to a shock propagating along the merger axis. It is possible that the same shock has also produced the phoenix through adiabatic compression, while the bridge could be generated by electrons which were pre-accelerated by the shock, and then re-accelerated by turbulence. Finally, we detect hints of gentle re-energisation in the two arms which depart from the tail of the radio galaxy.

preprint2022arXiv

Seeing the forest and the trees: a radio investigation of the ULIRG Mrk 273

Galaxy mergers have been observed to trigger nuclear activity by feeding gas to the central supermassive black hole. One such class of objects are Ultra Luminous InfraRed Galaxies (ULIRGs), which are mostly late stage major mergers of gas-rich galaxies. Recently, large-scale ($\sim$100 kpc) radio continuum emission has been detected in a select number of ULIRGs, all of which also harbour powerful Active Galactic Nuclei (AGN). This hints at the presence of large-scale radio emission being evidence for nuclear activity. Exploring the origin of this radio emission and its link to nuclear activity requires high sensitivity multi-frequency data. We present such an analysis of the ULIRG Mrk 273. Using the International LOFAR telescope (ILT), we detected spectacular large-scale arcs in this system. This detection includes, for the first time, a giant $\sim$190 kpc arc in the north. We propose these arcs are fuelled by a low power radio AGN triggered by the merger. We also identified a bright $\sim$45 kpc radio ridge, which is likely related to the ionised gas nebula in that region. We combined this with high sensitivity data from APERture Tile In Focus (Apertif) and archival data from the Very Large Array (VLA) to explore the spectral properties. The ILT simultaneously allowed us to probe the nucleus at a resolution of $\sim$0.3 arcsec, where we detected three components, and, for the first time, diffuse emission around these components. Combining this with archival high frequency VLA images of the nucleus allowed us to detect absorption in one component, and a steep spectrum radio AGN in another. We then extrapolate from this case study to the importance of investigating the presence of radio emission in more ULIRGs and what it can tell us about the link between mergers and the presence of radio activity.

preprint2022arXiv

The Coma cluster at LOFAR frequencies II: the halo, relic, and a new accretion relic

We present LOw Frequency ARray observations of the Coma cluster field at 144\,MHz. The cluster hosts one of the most famous radio halos, a relic, and a low surface-brightness bridge. We detect new features that allow us to make a step forward in the understanding of particle acceleration in clusters. The radio halo extends for more than 2 Mpc, which is the largest extent ever reported. To the North-East of the cluster, beyond the Coma virial radius, we discover an arc-like radio source that could trace particles accelerated by an accretion shock. To the West of the halo, coincident with a shock detected in the X-rays, we confirm the presence of a radio front, with different spectral properties with respect to the rest of the halo. We detect a radial steepening of the radio halo spectral index between 144 MHz and 342 MHz, at $\sim 30^{\prime}$ from the cluster centre, that may indicate a non constant re-acceleration time throughout the volume. We also detect a mild steepening of the spectral index towards the cluster centre. For the first time, a radial change in the slope of the radio-X-ray correlation is found, and we show that such a change could indicate an increasing fraction of cosmic ray versus thermal energy density in the cluster outskirts. Finally, we investigate the origin of the emission between the relic and the source NGC 4789, and we argue that NGC4789 could have crossed the shock originating the radio emission visible between its tail and the relic.

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

Radio halos in a mass-selected sample of 75 galaxy clusters. I. Sample selection and data analysis

Radio halos are synchrotron diffuse sources at the centre of a fraction of galaxy clusters. The study of large samples of clusters with adequate radio and X-ray data is necessary to investigate the origin of radio halos and their connection with the cluster dynamics and formation history. The aim of this paper is to compile a well-selected sample of galaxy clusters with deep radio observations to perform an unbiased statistical study of the properties of radio halos. We selected 75 clusters with M > = 6e14 Msun at z=0.08-0.33 from the Planck Sunyaev-Zel'dovich catalogue. Clusters without suitable radio data were observed with the Giant Metrewave Radio Telescope (GMRT) and/or the Jansky Very Large Array (JVLA) to complete the information about the possible presence of diffuse emission. We used archival Chandra X-ray data to derive information on the clusters' dynamical states. This observational campaign led to the detection of several cluster-scale diffuse radio sources and candidates that deserve future follow-up observations. Here we summarise their properties and add information resulting from our new observations. For the clusters where we did not detect any hint of diffuse emission, we derived new upper limits to their diffuse flux. We have built the largest mass-selected (> 80 per cent complete in mass) sample of galaxy clusters with deep radio observations available to date. The statistical analysis of the sample, which includes the connection between radio halos and cluster mergers, the radio power-mass correlation, and the occurrence of radio halos as a function of the cluster mass, will be presented in paper II.

preprint2021arXiv

Radio halos in a mass-selected sample of 75 galaxy clusters. II. Statistical analysis

Many galaxy clusters host Mpc scale diffuse radio sources called radio halos. Their origin is connected to the processes that lead to the formation of clusters themselves. In order to unveil this connection, statistical studies of radio halos are necessary. We selected a sample of galaxy clusters with M500>6e14Msun and z=0.08-0.33 from the Planck SZ catalogue. In paper I, we presented the radio and X-ray data analysis that we carried out on these clusters. Here, we study the radio properties of the sample, in connection to the mass and dynamical state of clusters. We used the dynamical information derived from the X-ray data to assess the role of mergers in the origin of radio halos. We studied the distribution of clusters in the radio power-mass diagram and the role of dynamics on the radio luminosity and emissivity of radio halos. We measured the occurrence of radio halos as a function of the cluster mass and we compared it with the expectations of turbulent acceleration models. We found that more than the 90% of radio halos are in merging clusters and that their radio power correlates with the mass of the host clusters. The correlation shows a large dispersion. Interestingly, we showed that cluster dynamics contributes significantly to this dispersion with more disturbed clusters being more radio luminous. Clusters without radio halos are generally relaxed and the upper limits to their diffuse emission lie below the correlation. We showed that the radio emissivity of clusters exhibits an apparent bimodality, with the emissivity of radio halos being at least 5 times larger than the non-emission associated with more relaxed clusters. We found that the fraction of radio halos drops from ~70% in high mass clusters to ~35% in the lower mass systems of the sample and we showed that this result is in good agreement with the expectations from turbulent re-acceleration models.

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

Sub-arcsecond imaging with the International LOFAR Telescope: II. Completion of the LOFAR Long-Baseline Calibrator Survey

The Low-Frequency Array (LOFAR) Long-Baseline Calibrator Survey (LBCS) was conducted between 2014 and 2019 in order to obtain a set of suitable calibrators for the LOFAR array. In this paper we present the complete survey, building on the preliminary analysis published in 2016 which covered approximately half the survey area. The final catalogue consists of 30006 observations of 24713 sources in the northern sky, selected for a combination of high low-frequency radio flux density and flat spectral index using existing surveys (WENSS, NVSS, VLSS, and MSSS). Approximately one calibrator per square degree, suitable for calibration of $\geq$ 200 km baselines is identified by the detection of compact flux density, for declinations north of 30 degrees and away from the Galactic plane, with a considerably lower density south of this point due to relative difficulty in selecting flat-spectrum candidate sources in this area of the sky. Use of the VLBA calibrator list, together with statistical arguments by comparison with flux densities from lower-resolution catalogues, allow us to establish a rough flux density scale for the LBCS observations, so that LBCS statistics can be used to estimate compact flux densities on scales between 300 mas and 2 arcsec, for sources observed in the survey. The LBCS can be used to assess the structures of point sources in lower-resolution surveys, with significant reductions in the degree of coherence in these sources on scales between 2 arcsec and 300 mas. The LBCS survey sources show a greater incidence of compact flux density in quasars than in radio galaxies, consistent with unified schemes of radio sources. Comparison with samples of sources from interplanetary scintillation (IPS) studies with the Murchison Widefield Array (MWA) shows consistent patterns of detection of compact structure in sources observed both interferometrically with LOFAR and using IPS.

preprint2021arXiv

The eROSITA Final Equatorial-Depth Survey (eFEDS): LOFAR view of brightest cluster galaxies and AGN feedback

During the performance verification phase of the SRG/eROSITA telescope, the eROSITA Final Equatorial-Depth Survey (eFEDS) has been carried out. It covers a 140 deg$^2$ field located at 126$^\circ <$ R.A. $< 146^\circ$ and -3$^\circ <$ Dec. $< +6^\circ$ with a nominal exposure over the field of 2.2 ks. 542 candidate clusters were detected in this field, down to a flux limit $F_X \sim 10^{-14}$ erg s$^{-1}$ cm$^{-2}$ in the 0.5-2 keV band. In order to understand radio-mode feedback in galaxy clusters, we study the radio emission of brightest cluster galaxies of eFEDS clusters, and we relate it to the X-ray properties of the host cluster. Using LOFAR we identify 227 radio galaxies hosted in the BCGs of the 542 galaxy clusters and groups detected in eFEDS. We treat non-detections as radio upper limits. We analyse the properties of radio galaxies, such as redshift and luminosity distribution, offset from the cluster centre, largest linear size and radio power. We study their relation to the intracluster medium of the host cluster. We perform statistical tests to deal with upper limits on the radio luminosities. BCGs with radio-loud AGN are more likely to lie close to the cluster centre than radio-quiet BCGs. There is a clear relation between the cluster's X-ray luminosity and the radio power of the BCG. Statistical tests indicate that this correlation is not produced by selection effects in the radio band. We see no apparent link between largest linear size of the radio galaxy and central density of the host cluster. Converting the radio luminosity to kinetic luminosity, we find that radiative losses of the intracluster medium are in an overall balance with the heating provided by the central AGN. Finally, we tentatively classify our objects into disturbed and relaxed, and we show that the link between the AGN and the ICM apparently holds regardless of the dynamical state of the cluster.

preprint2021arXiv

The intra-cluster magnetic field in the double relic galaxy cluster Abell 2345

Magnetic fields are ubiquitous in galaxy clusters, yet their radial profile, power spectrum, and connection to host cluster properties are poorly known. Merging galaxy clusters hosting diffuse polarized emission in the form of radio relics offer a unique possibility to study the magnetic fields in these complex systems. In this paper, we investigate the intra-cluster magnetic field in Abell 2345. This cluster hosts two radio relics that we detected in polarization with 1-2 GHz JVLA observations. X-ray XMM-Newton images show a very disturbed morphology. We derived the Rotation Measure (RM) of five polarized sources within $\sim$ 1 Mpc from the cluster center applying the RM synthesis. Both, the average RM and the RM dispersion radial profiles probe the presence of intra-cluster magnetic fields. Using the thermal electron density profile derived from X-ray analysis and simulating a 3D magnetic field with fluctuations following a power spectrum derived from magneto-hydrodynamical cosmological simulations, we build mock RM images of the cluster. We constrained the magnetic field profile in the eastern radio relic sector by comparing simulated and observed RM images. We find that, within the framework of our model, the data require a magnetic field scaling with thermal electron density as $B(r)\propto n_e(r)$. The best model has a central magnetic field (within a 200 kpc radius) of $2.8\pm0.1$ $μ$G. The average magnetic field at the position of the eastern relic is $\sim$0.3 $μ$G, a factor 2.7 lower than the equipartition estimate.

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

Cassiopeia A, Cygnus A, Taurus A, and Virgo A at ultra-low radio frequencies

The four persistent radio sources in the northern sky with the highest flux density at metre wavelengths are Cassiopeia A, Cygnus A, Taurus A, and Virgo A; collectively they are called the A-team. Their flux densities at ultra-low frequencies (<100 MHz) can reach several thousands of janskys, and they often contaminate observations of the low-frequency sky by interfering with image processing. Furthermore, these sources are foreground objects for all-sky observations hampering the study of faint signals, such as the cosmological 21 cm line from the epoch of reionisation. We aim to produce robust models for the surface brightness emission as a function of frequency for the A-team sources at ultra-low frequencies. These models are needed for the calibration and imaging of wide-area surveys of the sky with low-frequency interferometers. This requires obtaining images at an angular resolution better than 15 arcsec with a high dynamic range and good image fidelity. We observed the A-team with the Low Frequency Array (LOFAR) at frequencies between 30 MHz and 77 MHz using the Low Band Antenna (LBA) system. We reduced the datasets and obtained an image for each A-team source. The paper presents the best models to date for the sources Cassiopeia A, Cygnus A, Taurus A, and Virgo A between 30 MHz and 77 MHz. We were able to obtain the aimed resolution and dynamic range in all cases. Owing to its compactness and complexity, observations with the long baselines of the International LOFAR Telescope will be required to improve the source model for Cygnus A further.

preprint2020arXiv

Diffuse radio emission in the galaxy cluster SPT-CL J2031-4037: a steep spectrum intermediate radio halo?

The advent of sensitive low frequency radio observations has revealed a number of diffuse radio objects with peculiar properties that are challenging our understanding about the physics of the intracluster medium. Here, we report the discovery of a steep spectrum radio halo surrounding the central Brightest Cluster Galaxy (BCG) in the galaxy cluster SPT-CL J2031-4037. This cluster is morphologically disturbed yet has a weak cool core, an example of cool core/non-cool core transition system, which harbours a radio halo of $\sim 0.7$ Mpc in size. The halo emission detected at 1.7 GHz is less extended compared to that in the 325 MHz observation, and the spectral index of the part of the halo visible at 325 MHz to 1.7 GHz frequencies was found to be $-1.35 \pm 0.07$. Also, $P_{1.4\ \mathrm{GHz}}$ was found to be $0.77 \times 10^{24}$ W Hz$^{-1}$ which falls in the region where radio mini-halos, halo upper limits and ultra-steep spectrum (USS) halos are found in the $P_{1.4\ \mathrm{GHz}} - L_\mathrm{X}$ plane. Additionally, simulations presented in the paper provide support to the scenario of the steep spectrum. The diffuse radio emission found in this cluster may be a steep spectrum "intermediate" or "hybrid" radio halo which is transitioning into a mini-halo.

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

New mysteries and challenges from the Toothbrush relic: wideband observations from 550 MHz to 8 GHz

Radio relics are diffuse extended synchrotron sources that originate from shock fronts induced by galaxy cluster mergers. The actual particle acceleration mechanism at the shock fronts is still under debate. The galaxy cluster 1RXS J0603.3+4214 hosts one of the most intriguing examples of radio relics, known as the Toothbrush. We present new wideband radio continuum observations made with uGMRT and VLA. Our new observations, in combination with previously published data, allowed us to carry out a detailed high spatial resolution spectral and curvature analysis of the known diffuse radio emission sources, over a broad range of frequencies. The integrated spectrum of the Toothbrush follows closely a power law over close to 2 decades in frequency, with a spectral index of $-1.16\pm0.02$. We do not find any evidence of spectral steepening below 8 GHz. The subregions of the Toothbrush show an identical spectral slopes, suggesting that observed spectral index is rather set by the distribution of Mach numbers which may have a similar shape at different parts of the shock front. Indeed, numerical simulations show an intriguing similar spectral index, indicating that the radio spectrum is dominated by the average over the inhomogeneities within the shock, with most of the emission coming from the tail of the Mach number distribution. In contrast to the Toothbrush, the spectrum of the fainter relics show a high frequency steepening. The integrated spectrum of the halo follows a power law from 150 MHz to 3 GHz with a spectral index of $-1.16\pm0.04$. We do not find any evidence for spectral curvature, not even in subareas of the halo. This suggest a homogeneous acceleration throughout the cluster volume. Between the brush region of the Toothbrush and the halo, the color-color analysis revealed emission that was consistent with an overlap between the two different spectral regions.

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 BUFFALO HST Survey

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

preprint2020arXiv

The discovery of radio halos in the Frontier Fields clusters Abell S1063 and Abell 370

Massive merging galaxy clusters often host diffuse Mpc-scale radio synchrotron emission. This emission originates from relativistic electrons in the ionized intracluster medium (ICM). An important question is how these synchrotron emitting relativistic electrons are accelerated. Our aim is to search for diffuse emission in the Frontier Fields clusters Abell S1063 and Abell 370 and characterize its properties. While these clusters are very massive and well studied at some other wavelengths, no diffuse emission has been reported for these clusters so far. We obtained 325 MHz Giant Metrewave Radio Telescope (GMRT) and 1--4 GHz Jansky Very Large Array (VLA) observations of Abell S1063 and Abell 370. We complement these data with Chandra and XMM-Newton X-ray observations. In our sensitive images, we discover radio halos in both clusters. In Abell S1063, a giant radio halo is found with a size of $\sim 1.2$ Mpc. The integrated spectral index between 325 MHz and 1.5 GHz is $-0.94\pm0.08$ and it steepens to $-1.77 \pm 0.20$ between 1.5 and 3.0 GHz. This spectral steepening provides support for the turbulent re-acceleration model for radio halo formation. Abell 370 hosts a faint radio halo mostly centred on the southern part of this binary merging cluster, with a size of $\sim 500-700$ kpc. The spectral index between 325 MHz and 1.5 GHz is $-1.10\pm0.09$. Both radio halos follow the known scaling relation between the cluster mass proxy $Y_{500}$ and radio power, consistent with the idea that they are related to ongoing cluster merger events.

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

Anatomy of a Cooling Flow: The Feedback Response to Pure Cooling in the Core of the Phoenix Cluster

We present new, deep observations of the Phoenix cluster from the Chandra X-ray Observatory, the Hubble Space Telescope, and the Karl Jansky Very Large Array. These data provide an order of magnitude improvement in depth and/or angular resolution at X-ray, optical, and radio wavelengths, yielding an unprecedented view of the core of the Phoenix cluster. We find that the one-dimensional temperature and entropy profiles are consistent with expectations for pure-cooling hydrodynamic simulations and analytic descriptions of homogeneous, steady-state cooling flow models. In the inner ~10 kpc, the cooling time is shorter by an order of magnitude than any other known cluster, while the ratio of the cooling time to freefall time approaches unity, signaling that the ICM is unable to resist multiphase condensation on kpc scales. When we consider the thermodynamic profiles in two dimensions, we find that the cooling is highly asymmetric. The bulk of the cooling in the inner ~20 kpc is confined to a low-entropy filament extending northward from the central galaxy. We detect a substantial reservoir of cool (10^4 K) gas (as traced by the [OII] doublet), which is coincident with the low-entropy filament. The bulk of this cool gas is draped around and behind a pair of X-ray cavities, presumably bubbles that have been inflated by radio jets, which are detected for the first time on kpc scales. These data support a picture in which AGN feedback is promoting the formation of a multiphase medium via a combination of ordered buoyant uplift and locally enhanced turbulence. These processes ought to counteract the tendency for buoyancy to suppress condensation, leading to rapid cooling along the jet axis. The recent mechanical outburst has sufficient energy to offset cooling, and appears to be coupling to the ICM via a cocoon shock, raising the entropy in the direction orthogonal to the radio jets.

preprint2019arXiv

Probing the origin of diffuse radio emission in the cool-core of the Phoenix galaxy cluster

Cool core galaxy clusters are considered to be dynamically relaxed clusters with regular morphology and highly X-ray luminous central region. However, cool core clusters can also be sites for merging events that exhibit cold fronts in X-ray and mini-halos in radio. We present recent radio/X-ray observations of the Phoenix Cluster or SPT-CL J2344-4243 at the redshift of $z=0.596$. Using archival {\it Chandra} X-ray observations, we detect spiraling cool gas around the cluster core as well as discover two cold fronts near the core. It is perhaps the most distant galaxy cluster to date known to host cold fronts. Also, we present JVLA\footnote{Jansky Very Large Array\\ \url{https://science.nrao.edu/facilities/vla}} 1.52 GHz observations of the minihalo, previously discovered at 610 MHz with GMRT\footnote{Giant Metrewave Radio Telescope\\ \url{http://www.gmrt.ncra.tifr.res.in}} observations in the center of the Phoenix galaxy cluster. The minihalo flux density at 1.52 GHz is $9.65 \pm 0.97$ mJy with the spectral index between 610 MHz and 1.52 GHz being $-0.98 \pm 0.16$\footnote{$S_ν = ν^α$ where $S_ν$}. A possible origin of these radio sources is turbulence induced by sloshing of the gas in the cluster core.

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

The Karl G. Jansky Very Large Array Sky Survey (VLASS). Science case and survey design

The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$σ$ goal of 70 $μ$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.

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 LOFAR mini-survey for low-frequency radio emission from the nearest brown dwarfs

We have conducted a mini-survey for low-frequency radio emission from some of the closest brown dwarfs to the Sun with rapid rotation rates: SIMP J013656.5+093347, WISEPC J150649.97+702736.0, and WISEPA J174124.26+255319.5. We have placed robust 3-sigma upper limits on the flux density in the 111 - 169 MHz frequency range for these targets: WISE 1506: < 0.72 mJy; WISE 1741: < 0.87 mJy; SIMP 0136: < 0.66 mJy. At 8 hours of integration per target to achieve these limits, we find that systematic and detailed study of this class of object at LOFAR frequencies will require a substantial dedication of resources.

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

Carbon and hydrogen radio recombination lines from the cold clouds towards Cassiopeia A

We use the Low Frequency Array to perform a systematic high spectral resolution investigation of the low-frequency 33-78 MHz spectrum along the line of sight to Cassiopeia A. We complement this with a 304-386 MHz Westerbork Synthesis Radio telescope observation. In this first paper we focus on the carbon radio recombination lines. We detect Cn$α$ lines at -47 and -38 km s$^{-1}$ in absorption for quantum numbers n=438-584 and in emission for n=257-278 with high signal to noise. These lines are associated with cold clouds in the Perseus spiral arm component. Hn$α$ lines are detected in emission for n=257-278. In addition, we also detect Cn$α$ lines at 0 km s$^{-1}$ associated with the Orion arm. We analyze the optical depth of these transitions and their line width. Our models show that the carbon line components in the Perseus arm are best fit with an electron temperature 85 K and an electron density 0.04 cm$^{-3}$ and can be constrained to within 15\%. The electron pressure is constrained to within 20\%. We argue that much of these carbon radio recombination lines arise in the CO-dark surface layers of molecular clouds where most of the carbon is ionized but hydrogen has made the transition from atomic to molecular. The hydrogen lines are clearly associated with the carbon line emitting clouds, but the low-frequency upperlimits indicate that they likely do not trace the same gas. Combining the hydrogen and carbon results we arrive at a firm lower limit to the cosmic ray ionization rate of 2.5$\times$10$^{-18}$ s$^{-1}$, but the actual value is likely much larger.

preprint2016arXiv

CLASH-VLT: Dissecting the Frontier Fields Galaxy Cluster MACS J0416.1-2403 with $\sim800$ Spectra of Member Galaxies

We present VIMOS-VLT spectroscopy of the Frontier Fields cluster MACS~J0416.1-2403 (z=0.397). Taken as part of the CLASH-VLT survey, the large spectroscopic campaign provided more than 4000 reliable redshifts, including ~800 cluster member galaxies. The unprecedented sample of cluster members at this redshift allows us to perform a highly detailed dynamical and structural analysis of the cluster out to ~2.2$r_{200}$ (~4Mpc). Our analysis of substructures reveals a complex system composed of a main massive cluster ($M_{200}$~0.9$\times 10^{15} M_{\odot}$) presenting two major features: i) a bimodal velocity distribution, showing two central peaks separated by $ΔV_{rf}$~1100 km s$^{-1}$ with comparable galaxy content and velocity dispersion, ii) a projected elongation of the main substructures along the NE-SW direction, with a prominent subclump ~600 kpc SW of the center and an isolated BCG approximately halfway between the center and the SW clump. We also detect a low mass structure at z~0.390, ~10' S of the cluster center, projected at ~3Mpc, with a relative line-of-sight velocity of $ΔV_{rf}$~-1700 km s$^{-1}$. The cluster mass profile that we obtain through our dynamical analysis deviates significantly from the "universal" NFW, being best fit by a Softened Isothermal Sphere model instead. The mass profile measured from the galaxy dynamics is found to be in relatively good agreement with those obtained from strong and weak lensing, as well as with that from the X-rays, despite the clearly unrelaxed nature of the cluster. Our results reveal overall a complex dynamical state of this massive cluster and support the hypothesis that the two main subclusters are being observed in a pre-collisional phase, in line with recent findings from radio and deep X-ray data. With this article we also release the entire redshift catalog of 4386 sources in the field of this cluster.

preprint2016arXiv

Frontier Fields Clusters: Deep Chandra Observations of the Complex Merger MACS J1149.6+2223

The HST Frontier Fields cluster MACS J1149.6+2223 is one of the most complex merging clusters, believed to consist of four dark matter halos. We present results from deep (365 ks) Chandra observations of the cluster, which reveal the most distant cold front (z=0.544) discovered to date. In the cluster outskirts, we also detect hints of a surface brightness edge that could be the bow shock preceding the cold front. The substructure analysis of the cluster identified several components with large relative radial velocities, thus indicating that at least some collisions occur almost along the line of sight. The inclination of the mergers with respect to the plane of the sky poses significant observational challenges at X-ray wavelengths. MACS J1149.6+2223 possibly hosts a steep-spectrum radio halo. If the steepness of the radio halo is confirmed, then the radio spectrum, combined with the relatively regular ICM morphology, could indicate that MACS J1149.6+2223 is an old merging cluster.

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

Multi-wavelength Observations of the Dissociative Merger in the Galaxy Cluster CIZA J0107.7+5408

We present results based on X-ray, optical, and radio observations of the massive galaxy cluster CIZA J0107.7+5408. We find that this system is a post core passage, dissociative, binary merger, with the optical galaxy density peaks of each subcluster leading their associated X-ray emission peaks. This separation occurs because the diffuse gas experiences ram pressure forces while the effectively collisionless galaxies (and presumably their associated dark matter halos) do not. This system contains double peaked diffuse radio emission, possibly a double radio relic with the relics lying along the merger axis and also leading the X-ray cores. We find evidence for a temperature peak associated with the SW relic, likely created by the same merger shock that is powering the relic radio emission in this region. Thus, this system is a relatively rare clean example of a dissociative binary merger, which can in principle be used to place constraints on the self-interaction cross-section of dark matter. Low frequency radio observations reveal ultra-steep spectrum diffuse radio emission that is not correlated with the X-ray, optical, or high frequency radio emission. We suggest that these sources are radio phoenixes, which are preexisting non-thermal particle populations that have been re-energized through adiabatic compression by the same merger shocks that power the radio relics. Finally, we place upper limits on inverse Compton emission from the SW radio relic.

preprint2016arXiv

Observations and modelling of CO and [CI] in disks. First detections of [CI] and constraints on the carbon abundance

The gas-solid budget of carbon in protoplanetary disks is related to the composition of the cores and atmospheres of the planets forming in them. The key gas-phase carbon carriers CO, C$^{0}$ and C$^{+}$ can now be observed in disks. The gas-phase carbon abundance in disks has not yet been well characterized, we aim to obtain new constraints on the [C]/[H] ratio in a sample of disks, and to get an overview of the strength of [CI] and warm CO emission. We carried out a survey of the CO$\,6$--$5$ and [CI]$\,1$--$0$ and $2$--$1$ lines towards $37$ disks with APEX, and supplemented it with [CII] data from the literature. The data are interpreted using a grid of models produced with the DALI code. We also investigate how well the gas-phase carbon abundance can be determined in light of parameter uncertainties. The CO$\,6$--$5$ line is detected in $13$ out of $33$ sources, the [CI]$\,1$--$0$ in $6$ out of $12$, and the [CI]$\,2$--$1$ in $1$ out of $33$. With deep integrations, the first unambiguous detections of [CI]~$1$--$0$ in disks are obtained, in TW~Hya and HD~100546. Gas-phase carbon abundance reductions of a factor $5$--$10$ or more can be identified robustly based on CO and [CI] detections. The atomic carbon detection in TW~Hya confirms a factor $100$ reduction of [C]/[H]$_{\rm gas}$ in that disk, while the data are consistent with an ISM-like carbon abundance for HD~100546. In addition, BP~Tau, T~Cha, HD~139614, HD~141569, and HD~100453 are either carbon-depleted or gas-poor disks. The low [CI]~$2$--$1$ detection rates in the survey mostly reflect insufficient sensitivity to detect T~Tauri disks. The Herbig~Ae/Be disks with CO and [CII] upper limits below the models are debris disk like systems. A roughly order of magnitude increase in sensitivity compared to our survey is required to obtain useful constraints on the gas-phase [C]/[H] ratio in most of the targeted systems.

preprint2016arXiv

Probing Ionospheric Structures using the LOFAR radio telescope

LOFAR is the LOw Frequency Radio interferometer ARray located at mid-latitude ($52^{\circ} 53'N$). Here, we present results on ionospheric structures derived from 29 LOFAR nighttime observations during the winters of 2012/2013 and 2013/2014. We show that LOFAR is able to determine differential ionospheric TEC values with an accuracy better than 1 mTECU over distances ranging between 1 and 100 km. For all observations the power law behavior of the phase structure function is confirmed over a long range of baseline lengths, between $1$ and $80$ km, with a slope that is in general larger than the $5/3$ expected for pure Kolmogorov turbulence. The measured average slope is $1.89$ with a one standard deviation spread of $0.1$. The diffractive scale, i.e. the length scale where the phase variance is $1\, \mathrm{rad^2}$, is shown to be an easily obtained single number that represents the ionospheric quality of a radio interferometric observation. A small diffractive scale is equivalent to high phase variability over the field of view as well as a short time coherence of the signal, which limits calibration and imaging quality. For the studied observations the diffractive scales at $150$ MHz vary between $3.5$ and $30\,$ km. A diffractive scale above $5$ km, pertinent to about $90 \%$ of the observations, is considered sufficient for the high dynamic range imaging needed for the LOFAR Epoch of Reionization project. For most nights the ionospheric irregularities were anisotropic, with the structures being aligned with the Earth magnetic field in about $60\%$ of the observations.

preprint2016arXiv

The discovery of lensed radio and X-ray sources behind the Frontier Fields cluster MACS J0717.5+3745 with the JVLA and Chandra

We report on high-resolution JVLA and Chandra observations of the HST Frontier Cluster MACS J0717.5+3745. MACS J0717.5+3745 offers the largest contiguous magnified area of any known cluster, making it a promising target to search for lensed radio and X-ray sources. With the high-resolution 1.0-6.5 GHz JVLA imaging in A and B configuration, we detect a total of 51 compact radio sources within the area covered by the HST imaging. Within this sample we find 7 lensed sources with amplification factors larger than $2$. None of these sources are identified as multiply-lensed. Based on the radio luminosities, the majority of these sources are likely star forming galaxies with star formation rates of 10-50 M$_\odot$ yr$^{-1}$ located at $1 \lesssim z \lesssim 2$. Two of the lensed radio sources are also detected in the Chandra image of the cluster. These two sources are likely AGN, given their $2-10$ keV X-ray luminosities of $\sim 10^{43-44}$ erg s$^{-1}$. From the derived radio luminosity function, we find evidence for an increase in the number density of radio sources at $0.6<z<2.0$, compared to a $z < 0.3$ sample. Our observations indicate that deep radio imaging of lensing clusters can be used to study star forming galaxies, with star formation rates as low as $\sim10$ M$_{\odot}$ yr$^{-1}$, at the peak of cosmic star formation history.

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.

preprint2015arXiv

$Suzaku$ X-ray study of the double radio relic galaxy cluster CIZA J2242.8+5301

Content: We present the results from $Suzaku$ observations of the merging cluster of galaxies CIZA J2242.8+5301 at $z$=0.192. Aims. To study the physics of gas heating and particle acceleration in cluster mergers, we investigated the X-ray emission from CIZA J2242.8+5301, which hosts two giant radio relics in the northern/southern part of the cluster. Methods. We analyzed data from three-pointed Suzaku observations of CIZA J2242.8+5301 to derive the temperature distribution in four different directions. Results: The Intra-Cluster Medium (ICM) temperature shows a remarkable drop from 8.5$_{-0.6}^{+0.8}$ keV to 2.7$_{-0.4}^{+0.7}$ keV across the northern radio relic. The temperature drop is consistent with a Mach number ${\cal M}_n=2.7^{+0.7}_{-0.4}$ and a shock velocity $v_{shock:n}=2300_{-400}^{+700}\rm\,km\,s^{-1}$. We also confirm the temperature drop across the southern radio relic. However, the ICM temperature beyond this relic is much higher than beyond the northern one, which gives a Mach number ${\cal M}_s=1.7^{+0.4}_{-0.3}$ and shock velocity $v_{shock:s}=2040_{-410}^{+550}\rm \,km\,s^{-1}$. These results agree with other systems showing a relationship between the radio relics and shock fronts which are induced by merging activity. We compare the X-ray derived Mach numbers with the radio derived Mach numbers from the radio spectral index under the assumption of diffusive shock acceleration in the linear test particle regime. For the northern radio relic, the Mach numbers derived from X-ray and radio observations agree with each other. Based on the shock velocities, we estimate that CIZA J2242.8+5301 is observed approximately 0.6 Gyr after core passage. The magnetic field pressure at the northern relic is estimated to be 9% of the thermal pressure.

preprint2015arXiv

A powerful double radio relic system discovered in PSZ1 G108.18-11.53: evidence for a shock with non-uniform Mach number?

Diffuse radio emission in the form of radio halos and relics has been found in a number of merging galaxy clusters. These structures indicate that shock and turbulence associated with the merger accelerate electrons to relativistic energies. We report the discovery of a radio relic + radio halo system in PSZ1 G108.18-11.53 (z=0.335). This cluster hosts the second most powerful double radio relic system ever discovered. We observed PSZ1 G108.18-11.53 with the Giant Meterwave Radio Telescope (GMRT) and the Westerbork Synthesis Radio Telescope (WSRT). We obtained radio maps at 147, 323, 607 and 1380 MHz. We also observed the cluster with the Keck telescope, obtaining the spectroscopic redshift for 42 cluster members. From the injection index we obtained the Mach number of the shocks generating the two radio relics. For the southern shock we found M = 2.33^{+0.19}_{-0.26}, while the northern shock Mach number goes from M = 2.20^{+0.07}_{-0.14} in the north part down to M = 2.00^{+0.03}_{-0.08} in the southern region. If the relation between the injection index and the Mach number predicted by diffusive shock acceleration (DSA) theory holds, this is the first observational evidence for a gradient in the Mach number along a galaxy cluster merger shock.

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

Calibrating the absolute amplitude scale for air showers measured at LOFAR

Air showers induced by cosmic rays create nanosecond pulses detectable at radio frequencies. These pulses have been measured successfully in the past few years at the LOw Frequency ARray (LOFAR) and are used to study the properties of cosmic rays. For a complete understanding of this phenomenon and the underlying physical processes, an absolute calibration of the detecting antenna system is needed. We present three approaches that were used to check and improve the antenna model of LOFAR and to provide an absolute calibration of the whole system for air shower measurements. Two methods are based on calibrated reference sources and one on a calibration approach using the diffuse radio emission of the Galaxy, optimized for short data-sets. An accuracy of 19% in amplitude is reached. The absolute calibration is also compared to predictions from air shower simulations. These results are used to set an absolute energy scale for air shower measurements and can be used as a basis for an absolute scale for the measurement of astronomical transients with LOFAR.

preprint2015arXiv

Chandra and XMM-Newton observations of the merging cluster of galaxies PLCK G036.7+14.9

We present Chandra and XMM-Newton observations of PLCK G036.7+14.9 from the Chandra-Planck Legacy Program. The high resolution X-ray observations reveal two close subclusters, G036N and G036S, which were not resolved by previous ROSAT, optical, or recent Planck observations. We perform detailed imaging and spectral analyses and use a simplified model to study the kinematics of this system. The basic picture is that PLCK G036.7+14.9 is undergoing a major merger (mass ratio close to unity) between the two massive subclusters, with the merger largely along the line-of-sight and probably at an early stage. G036N hosts a small, moderate cool-core, while G036S has at most a very weak cool-core in the central 40 kpc region. The difference in core cooling times is unlikely to be caused by the ongoing merger disrupting a pre-existing cool-core in G036S. G036N also hosts an unresolved radio source in the center, which may be heating the gas if the radio source is extended. The Planck derived mass is higher than the X-ray measured mass of either subcluster, but is lower than the X-ray measured mass of the whole cluster, due to the fact that Planck does not resolve PLCK G036.7+14.9 into subclusters and interprets it as a single cluster. This mass discrepancy could induce significant bias to the mass function if such previously unresolved systems are common in the Planck cluster sample. High resolution X-ray observations are necessary to identify the fraction of such systems and correct such a bias for the purpose of precision cosmological studies.

preprint2015arXiv

Deep Chandra, HST-COS, and Megacam Observations of the Phoenix Cluster: Extreme Star Formation and AGN Feedback on Hundred Kiloparsec Scales

We present new ultraviolet, optical, and X-ray data on the Phoenix galaxy cluster (SPT-CLJ2344-4243). Deep optical imaging reveals previously-undetected filaments of star formation, extending to radii of ~50-100 kpc in multiple directions. Combined UV-optical spectroscopy of the central galaxy reveals a massive (2x10^9 Msun)), young (~4.5 Myr) population of stars, consistent with a time-averaged star formation rate of 610 +/- 50 Msun/yr. We report a strong detection of OVI(1032,1038) which appears to originate primarily in shock-heated gas, but may contain a substantial contribution (>1000 Msun/yr) from the cooling intracluster medium. We confirm the presence of deep X-ray cavities in the inner ~10 kpc, which are amongst the most extreme examples of radio-mode feedback detected to date, implying jet powers of 2-7 x10^45 erg/s. We provide evidence that the AGN inflating these cavities may have only recently transitioned from "quasar-mode" to "radio-mode", and may currently be insufficient to completely offset cooling. A model-subtracted residual X-ray image reveals evidence for prior episodes of strong radio-mode feedback at radii of ~100 kpc, with extended "ghost" cavities indicating a prior epoch of feedback roughly 100 Myr ago. This residual image also exhibits significant asymmetry in the inner ~200 kpc (0.15R500), reminiscent of infalling cool clouds, either due to minor mergers or fragmentation of the cooling ICM. Taken together, these data reveal a rapidly evolving cool core which is rich with structure (both spatially and in temperature), is subject to a variety of highly energetic processes, and yet is cooling rapidly and forming stars along thin, narrow filaments.

preprint2015arXiv

Imaging Jupiter's radiation belts down to 127 MHz with LOFAR

Context. Observing Jupiter's synchrotron emission from the Earth remains today the sole method to scrutinize the distribution and dynamical behavior of the ultra energetic electrons magnetically trapped around the planet (because in-situ particle data are limited in the inner magnetosphere). Aims. We perform the first resolved and low-frequency imaging of the synchrotron emission with LOFAR at 127 MHz. The radiation comes from low energy electrons (~1-30 MeV) which map a broad region of Jupiter's inner magnetosphere. Methods (see article for complete abstract) Results. The first resolved images of Jupiter's radiation belts at 127-172 MHz are obtained along with total integrated flux densities. They are compared with previous observations at higher frequencies and show a larger extent of the synchrotron emission source (>=4 $R_J$). The asymmetry and the dynamic of east-west emission peaks are measured and the presence of a hot spot at lambda_III=230 ° $\pm$ 25 °. Spectral flux density measurements are on the low side of previous (unresolved) ones, suggesting a low-frequency turnover and/or time variations of the emission spectrum. Conclusions. LOFAR is a powerful and flexible planetary imager. The observations at 127 MHz depict an extended emission up to ~4-5 planetary radii. The similarities with high frequency results reinforce the conclusion that: i) the magnetic field morphology primarily shapes the brightness distribution of the emission and ii) the radiating electrons are likely radially and latitudinally distributed inside about 2 $R_J$. Nonetheless, the larger extent of the brightness combined with the overall lower flux density, yields new information on Jupiter's electron distribution, that may shed light on the origin and mode of transport of these particles.

preprint2015arXiv

LOFAR discovery of a quiet emission mode in PSR B0823+26

PSR B0823+26, a 0.53-s radio pulsar, displays a host of emission phenomena over timescales of seconds to (at least) hours, including nulling, subpulse drifting, and mode-changing. Studying pulsars like PSR B0823+26 provides further insight into the relationship between these various emission phenomena and what they might teach us about pulsar magnetospheres. Here we report on the LOFAR discovery that PSR B0823+26 has a weak and sporadically emitting 'quiet' (Q) emission mode that is over 100 times weaker (on average) and has a nulling fraction forty-times greater than that of the more regularly-emitting 'bright' (B) mode. Previously, the pulsar has been undetected in the Q-mode, and was assumed to be nulling continuously. PSR B0823+26 shows a further decrease in average flux just before the transition into the B-mode, and perhaps truly turns off completely at these times. Furthermore, simultaneous observations taken with the LOFAR, Westerbork, Lovell, and Effelsberg telescopes between 110 MHz and 2.7 GHz demonstrate that the transition between the Q-mode and B-mode occurs within one single rotation of the neutron star, and that it is concurrent across the range of frequencies observed.

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

LOFAR Sparse Image Reconstruction

Context. The LOw Frequency ARray (LOFAR) radio telescope is a giant digital phased array interferometer with multiple antennas distributed in Europe. It provides discrete sets of Fourier components of the sky brightness. Recovering the original brightness distribution with aperture synthesis forms an inverse problem that can be solved by various deconvolution and minimization methods Aims. Recent papers have established a clear link between the discrete nature of radio interferometry measurement and the "compressed sensing" (CS) theory, which supports sparse reconstruction methods to form an image from the measured visibilities. Empowered by proximal theory, CS offers a sound framework for efficient global minimization and sparse data representation using fast algorithms. Combined with instrumental direction-dependent effects (DDE) in the scope of a real instrument, we developed and validated a new method based on this framework Methods. We implemented a sparse reconstruction method in the standard LOFAR imaging tool and compared the photometric and resolution performance of this new imager with that of CLEAN-based methods (CLEAN and MS-CLEAN) with simulated and real LOFAR data Results. We show that i) sparse reconstruction performs as well as CLEAN in recovering the flux of point sources; ii) performs much better on extended objects (the root mean square error is reduced by a factor of up to 10); and iii) provides a solution with an effective angular resolution 2-3 times better than the CLEAN images. Conclusions. Sparse recovery gives a correct photometry on high dynamic and wide-field images and improved realistic structures of extended sources (of simulated and real LOFAR datasets). This sparse reconstruction method is compatible with modern interferometric imagers that handle DDE corrections (A- and W-projections) required for current and future instruments such as LOFAR and SKA

preprint2015arXiv

Outflow forces in intermediate mass star formation

Intermediate mass protostarsprovide a bridge between theories of low- and high-mass star formation. Emerging molecular outflows can be used to determine the influence of fragmentation and multiplicity on protostellar evolution through the correlation of outflow forces of intermediate mass protostars with the luminosity. The aim of this paper is to derive outflow forces from outflows of six intermediate mass protostellar regions and validate the apparent correlation between total luminosity and outflow force seen in earlier work, as well as remove uncertainties caused by different methodology. By comparing CO 6--5 observations obtained with APEX with non-LTE radiative transfer model predictions, optical depths, temperatures, densities of the gas of the molecular outflows are derived. Outflow forces, dynamical timescales and kinetic luminosities are subsequently calculated. Outflow parameters, including the forces, were derived for all sources. Temperatures in excess of 50 K were found for all flows, in line with recent low-mass results. However, comparison with other studies could not corroborate conclusions from earlier work on intermediate mass protostars which hypothesized that fragmentation enhances outflow forces in clustered intermediate mass star formation. Any enhancement in comparison with the classical relation between outflow force and luminosity can be attributed the use of a higher excitation line and improvement in methods; They are in line with results from low-mass protostars using similar techniques. The role of fragmentation on outflows is an important ingredient to understand clustered star formation and the link between low and high-mass star formation. However, detailed information on spatial scales of a few 100 AU, covering all individual members is needed to make the necessary progress.

preprint2015arXiv

Probing Atmospheric Electric Fields in Thunderstorms through Radio Emission from Cosmic-Ray-Induced Air Showers

We present measurements of radio emission from cosmic ray air showers that took place during thunderstorms. The intensity and polarization patterns of these air showers are radically different from those measured during fair-weather conditions. With the use of a simple two-layer model for the atmospheric electric field, these patterns can be well reproduced by state-of-the-art simulation codes. This in turn provides a novel way to study atmospheric electric fields.

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

Resolving the Merging Planck Cluster PLCK G147.3-16.6 with GISMO

The $Planck$ satellite has recently completed an all-sky galaxy cluster survey exploiting the thermal Sunyaev-Zel'dovich (SZ) effect to locate some of the most massive systems observable. With a median redshift of $\left<z\right>=0.22$, the clusters found by $Planck$ at $z>0.3$ are proving to be exceptionally massive and/or disturbed systems. One notable $Planck$ discovery at $z=0.645$, PLCK G147.3-16.6, has a dual core and hosts a radio halo, indicating it is likely in the process of merging. We present a 16."5 resolution SZ observation of this high-$z$ merger using the Goddard-IRAM Superconducting 2 Millimeter Observer (GISMO), and compare it to X-ray follow-up observations with XMM-$Newton$. We find the SZ pressure substructure is offset from the core components seen in X-ray. We interpret this offset as possible line of sight temperature or density substructure due to the on-going merger.

preprint2015arXiv

The LOFAR Multifrequency Snapshot Sky Survey (MSSS) I. Survey description and first results

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.

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.

preprint2015arXiv

The Physics of the Cold Neutral Medium: Low-frequency Carbon Radio Recombination Lines with the Square Kilometre Array

The Square Kilometre Array (SKA) will transform our understanding of the role of the cold, atomic gas in galaxy evolution. The interstellar medium (ISM) is the repository of stellar ejecta and the birthsite of new stars and, hence, a key factor in the evolution of galaxies over cosmic time. Cold, diffuse, atomic clouds are a key component of the ISM, but so far this phase has been difficult to study, because its main tracer, the HI 21 cm line, does not constrain the basic physical information of the gas (e.g., temperature, density) well. The SKA opens up the opportunity to study this component of the ISM through a complementary tracer in the form of low-frequency (<350 MHz) carbon radio recombination lines (CRRL). These CRRLs provide a sensitive probe of the physical conditions in cold, diffuse clouds. The superb sensitivity, large field of view, frequency resolution and coverage of the SKA allows for efficient surveys of the sky, that will revolutionize the field of low-frequency recombination line studies. By observing these lines with the SKA we will be able determine the thermal balance, chemical enrichment, and ionization rate of the cold, atomic medium from degree-scales down to scales corresponding to individual clouds and filaments in our Galaxy, the Magellanic Clouds and beyond. Furthermore, being sensitive only to the cold, atomic gas, observations of low-frequency CRRLs with the SKA will aid in disentangling the warm and cold constituents of the HI 21 cm emission.

preprint2015arXiv

Wide-Band, Low-Frequency Pulse Profiles of 100 Radio Pulsars with LOFAR

LOFAR offers the unique capability of observing pulsars across the 10-240 MHz frequency range with a fractional bandwidth of roughly 50%. This spectral range is well-suited for studying the frequency evolution of pulse profile morphology caused by both intrinsic and extrinsic effects: such as changing emission altitude in the pulsar magnetosphere or scatter broadening by the interstellar medium, respectively. The magnitude of most of these effects increases rapidly towards low frequencies. LOFAR can thus address a number of open questions about the nature of radio pulsar emission and its propagation through the interstellar medium. We present the average pulse profiles of 100 pulsars observed in the two LOFAR frequency bands: High Band (120-167 MHz, 100 profiles) and Low Band (15-62 MHz, 26 profiles). We compare them with Westerbork Synthesis Radio Telescope (WSRT) and Lovell Telescope observations at higher frequencies (350 and1400 MHz) in order to study the profile evolution. The profiles are aligned in absolute phase by folding with a new set of timing solutions from the Lovell Telescope, which we present along with precise dispersion measures obtained with LOFAR. We find that the profile evolution with decreasing radio frequency does not follow a specific trend but, depending on the geometry of the pulsar, new components can enter into, or be hidden from, view. Nonetheless, in general our observations confirm the widening of pulsar profiles at low frequencies, as expected from radius-to-frequency mapping or birefringence theories. We offer this catalog of low-frequency pulsar profiles in a user friendly way via the EPN Database of Pulsar Profiles (http://www.epta.eu.org/epndb/).

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

A new double radio relic in PSZ1 G096.89+24.17 and a radio relic mass-luminosity relation

Radio relics are diffuse synchrotron sources in galaxy clusters that are believed to trace large-scale shock waves. We have discovered a new double radio relic system in PSZ1 G096.89+24.17 (z=0.3) and have carried out a full-polarization radio observation using the Westerbork Synthesis Radio Telescope (WSRT) at 1.4 GHz. The observation revealed the presence of two relics located on the two diametrically opposite sides of the cluster and hints of a central radio halo. The linear sizes of the relics are $\sim 0.9$ and $\sim 1.4$ Mpc. We carried out an analysis of all known double radio relics by using radio, X-ray and Sunyaev-Zeldovich (SZ) data. We find that the radio luminosity of double relics is a steep function of the cluster mass, with $L_R \propto M^{2.83\pm0.39}$. If we include single radio relics, this relation is maintained. This dependence has implications for the origin of magnetic fields at the relic's locations.

preprint2014arXiv

Abell 1033: birth of a radio phoenix

Extended steep-spectrum radio emission in a galaxy cluster is usually associated with a recent merger. However, given the complex scenario of galaxy cluster mergers, many of the discovered sources hardly fit into the strict boundaries of a precise taxonomy. This is especially true for radio phoenixes that do not have very well defined observational criteria. Radio phoenixes are aged radio galaxy lobes whose emission is reactivated by compression or other mechanisms. Here, we present the detection of a radio phoenix close to the moment of its formation. The source is located in Abell 1033, a peculiar galaxy cluster which underwent a recent merger. To support our claim, we present unpublished Westerbork Synthesis Radio Telescope and Chandra observations together with archival data from the Very Large Array and the Sloan Digital Sky Survey. We discover the presence of two sub-clusters displaced along the N-S direction. The two sub-clusters probably underwent a recent merger which is the cause of a moderately perturbed X-ray brightness distribution. A steep-spectrum extended radio source very close to an AGN is proposed to be a newly born radio phoenix: the AGN lobes have been displaced/compressed by shocks formed during the merger event. This scenario explains the source location, morphology, spectral index, and brightness. Finally, we show evidence of a density discontinuity close to the radio phoenix and discuss the consequences of its presence.

preprint2014arXiv

Comparison of algorithms for determination of rotation measure and Faraday structure I. 1100 - 1400 MHz

(abridged) We run a Faraday structure determination data challenge to benchmark the currently available algorithms including Faraday synthesis (previously called RM synthesis in the literature), wavelet, compressive sampling and $QU$-fitting. The frequency set is similar to POSSUM/GALFACTS with a 300 MHz bandwidth from 1.1 to 1.4 GHz. We define three figures of merit motivated by the underlying science: a) an average RM weighted by polarized intensity, RMwtd, b) the separation $Δϕ$ of two Faraday components and c) the reduced chi-squared. Based on the current test data of signal to noise ratio of about 32, we find that: (1) When only one Faraday thin component is present, most methods perform as expected, with occasional failures where two components are incorrectly found; (2) For two Faraday thin components, QU-fitting routines perform the best, with errors close to the theoretical ones for RMwtd, but with significantly higher errors for $Δϕ$. All other methods including standard Faraday synthesis frequently identify only one component when $Δϕ$ is below or near the width of the Faraday point spread function; (3) No methods, as currently implemented, work well for Faraday thick components due to the narrow bandwidth; (4) There exist combinations of two Faraday components which produce a large range of acceptable fits and hence large uncertainties in the derived single RMs; in these cases, different RMs lead to the same Q, U behavior, so no method can recover a unique input model.

preprint2014arXiv

Diffuse radio emission in the complex merging galaxy cluster Abell 2069

Galaxy clusters with signs for a recent merger show in many cases extended diffuse radio features. This emission originates from relativistic electrons which suffer synchrotron losses due to the intra-cluster magnetic field. The mechanisms of the particle acceleration and the properties of the magnetic field are still poorly understood. We search for diffuse radio emission in galaxy clusters. Here, we study the complex galaxy cluster Abell 2069, for which X-ray observations indicate a recent merger. We investigate the cluster's radio continuum emission by deep Westerbork Synthesis Radio Telescope (WSRT) observations at 346 MHz and a Giant Metrewave Radio Telescope (GMRT) observation at 322 MHz. We find an extended diffuse radio feature roughly coinciding with the main component of the cluster. We classify this emission as a radio halo and estimate its lower limit flux density to 25 +/- 9 mJy. Moreover, we find a second extended diffuse source located at the cluster's companion and estimate its flux density to 15 +/- 2 mJy. We speculate that this is a small halo or a mini-halo. If true, this cluster is the first example of a double-halo in a single galaxy cluster.

preprint2014arXiv

Discovery of Carbon Radio Recombination Lines in absorption towards Cygnus~A

We present the first detection of carbon radio recombination line absorption along the line of sight to Cygnus A. The observations were carried out with the LOw Frequency ARray in the 33 to 57 MHz range. These low frequency radio observations provide us with a new line of sight to study the diffuse, neutral gas in our Galaxy. To our knowledge this is the first time that foreground Milky Way recombination line absorption has been observed against a bright extragalactic background source. By stacking 48 carbon $α$ lines in the observed frequency range we detect carbon absorption with a signal-to-noise ratio of about 5. The average carbon absorption has a peak optical depth of 2$\times$10$^{-4}$, a line width of 10 km s$^{-1}$ and a velocity of +4 km s$^{-1}$ with respect to the local standard of rest. The associated gas is found to have an electron temperature $T_{e}\sim$ 110 K and density $n_{e}\sim$ 0.06 cm$^{-3}$. These properties imply that the observed carbon $α$ absorption likely arises in the cold neutral medium of the Orion arm of the Milky Way. Hydrogen and helium lines were not detected to a 3$σ$ peak optical depth limit of 1.5$\times$10$^{-4}$ for a 4 km s$^{-1}$ channel width. Radio recombination lines associated with Cygnus A itself were also searched for, but are not detected. We set a 3$σ$ upper limit of 1.5$\times$10$^{-4}$ for the peak optical depth of these lines for a 4 km s$^{-1}$ channel width.

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

Initial LOFAR observations of Epoch of Reionization windows: II. Diffuse polarized emission in the ELAIS-N1 field

This study aims to characterise the polarized foreground emission in the ELAIS-N1 field and to address its possible implications for the extraction of the cosmological 21-cm signal from the Low-Frequency Array - Epoch of Reionization (LOFAR-EoR) data. We use the high band antennas of LOFAR to image this region and RM-synthesis to unravel structures of polarized emission at high Galactic latitudes. The brightness temperature of the detected Galactic emission is on average 4 K in polarized intensity and covers the range from -10 to +13rad m^-2 in Faraday depth. The total polarized intensity and polarization angle show a wide range of morphological features. We have also used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. The LOFAR and WSRT images show a similar complex morphology, at comparable brightness levels, but their spatial correlation is very low. The fractional polarization at 150 MHz, expressed as a percentage of the total intensity, amounts to 1.5%. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies. The wide frequency range, good angular resolution and good sensitivity make LOFAR an exquisite instrument for studying Galactic polarized emission at a resolution of 1-2 rad m^-2 in Faraday depth. The different polarised patterns observed at 150 MHz and 350 MHz are consistent with different source distributions along the line of sight wring in a variety of Faraday thin regions of emission. The presence of polarised foregrounds is a serious complication for Epoch of Reionization experiments. To avoid the leakage of polarized emission into total intensity, which can depend on frequency, we need to calibrate the instrumental polarization across the field of view to a small fraction of 1%.

preprint2014arXiv

LOFAR low-band antenna observations of the 3C295 and Bootes fields: source counts and ultra-steep spectrum sources

We present LOFAR Low Band observations of the Bootes and 3C295 fields. Our images made at 34, 46, and 62 MHz reach noise levels of 12, 8, and 5 mJy beam$^{-1}$, making them the deepest images ever obtained in this frequency range. In total, we detect between 300 and 400 sources in each of these images, covering an area of 17 to 52 deg$^{2}$. From the observations we derive Euclidean-normalized differential source counts. The 62 MHz source counts agree with previous GMRT 153 MHz and VLA 74 MHz differential source counts, scaling with a spectral index of $-0.7$. We find that a spectral index scaling of $-0.5$ is required to match up the LOFAR 34 MHz source counts. This result is also in agreement with source counts from the 38 MHz 8C survey, indicating that the average spectral index of radio sources flattens towards lower frequencies. We also find evidence for spectral flattening using the individual flux measurements of sources between 34 and 1400 MHz and by calculating the spectral index averaged over the source population. To select ultra-steep spectrum ($α< -1.1$) radio sources, that could be associated with massive high redshift radio galaxies, we compute spectral indices between 62 MHz, 153 MHz and 1.4 GHz for sources in the Boötes field. We cross-correlate these radio sources with optical and infrared catalogues and fit the spectral energy distribution to obtain photometric redshifts. We find that most of these ultra-steep spectrum sources are located in the $ 0.7 \lesssim z \lesssim 2.5$ range.

preprint2014arXiv

LOFAR tied-array imaging of Type III solar radio bursts

The Sun is an active source of radio emission which is often associated with energetic phenomena such as solar flares and coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), the Sun has not been imaged extensively because of the instrumental limitations of previous radio telescopes. Here, the combined high spatial, spectral and temporal resolution of the Low Frequency Array (LOFAR) was used to study solar Type III radio bursts at 30-90 MHz and their association with CMEs. The Sun was imaged with 126 simultaneous tied-array beams within 5 solar radii of the solar centre. This method offers benefits over standard interferometric imaging since each beam produces high temporal (83 ms) and spectral resolution (12.5 kHz) dynamic spectra at an array of spatial locations centred on the Sun. LOFAR's standard interferometric output is currently limited to one image per second. Over a period of 30 minutes, multiple Type III radio bursts were observed, a number of which were found to be located at high altitudes (4 solar radii from the solar center at 30 MHz) and to have non-radial trajectories. These bursts occurred at altitudes in excess of values predicted by 1D radial electron density models. The non-radial high altitude Type III bursts were found to be associated with the expanding flank of a CME. The CME may have compressed neighbouring streamer plasma producing larger electron densities at high altitudes, while the non-radial burst trajectories can be explained by the deflection of radial magnetic fields as the CME expanded in the low corona.

preprint2014arXiv

MC^2: Galaxy Imaging and Redshift Analysis of the Merging Cluster CIZA J2242.8+5301

X-ray and radio observations of CIZA J2242.8+5301 suggest that it is a major cluster merger. Despite being well studied in the X-ray, and radio, little has been presented on the cluster structure and dynamics inferred from its galaxy population. We carried out a deep (i<25) broad band imaging survey of the system with Subaru SuprimeCam (g & i bands) and the Canada France Hawaii Telescope (r band) as well as a comprehensive spectroscopic survey of the cluster area (505 redshifts) using Keck DEIMOS. We use this data to perform a comprehensive galaxy/redshift analysis of the system, which is the first step to a proper understanding the geometry and dynamics of the merger, as well as using the merger to constrain self-interacting dark matter. We find that the system is dominated by two subclusters of comparable richness with a projected separation of 6.9'^{+0.7}_{-0.5} (1.3^{+0.13}_{-0.10} Mpc). We find that the north and south subclusters have similar redshifts of z=0.188 with a relative line-of-sight velocity difference of 69+/-190 km/s. We also find that north and south subclusters have velocity dispersions of 1160^{+100}_{-90} km/s and 1080^{+100}_{-70} km/s, respectively. These correspond to masses of 16.1^{+4.6}_{-3.3}x10^14 M_sun and 13.0^{+4.0}_{-2.5}x10^14 M_sun, respectively. While velocity dispersion measurements of merging clusters can be biased we believe the bias in this system to be minor due to the large projected separation and nearly plane-of-sky merger configuration. CIZA J2242.8+5301 is a relatively clean dissociative cluster merger with near 1:1 mass ratio, which makes it an ideal merger for studying merger associated physical phenomena.

preprint2014arXiv

Measuring a Cherenkov ring in the radio emission from air showers at 110-190 MHz with LOFAR

Measuring radio emission from air showers offers a novel way to determine properties of the primary cosmic rays such as their mass and energy. Theory predicts that relativistic time compression effects lead to a ring of amplified emission which starts to dominate the emission pattern for frequencies above ~100 MHz. In this article we present the first detailed measurements of this structure. Ring structures in the radio emission of air showers are measured with the LOFAR radio telescope in the frequency range of 110 - 190 MHz. These data are well described by CoREAS simulations. They clearly confirm the importance of including the index of refraction of air as a function of height. Furthermore, the presence of the Cherenkov ring offers the possibility for a geometrical measurement of the depth of shower maximum, which in turn depends on the mass of the primary particle.

preprint2014arXiv

Multiple density discontinuities in the merging galaxy cluster CIZA J2242.8+5301

CIZA J2242.8+5301, a merging galaxy cluster at z=0.19, hosts a double-relic system and a faint radio halo. Radio observations at frequencies ranging from a few MHz to several GHz have shown that the radio spectral index at the outer edge of the N relic corresponds to a shock of Mach number 4.6+/-1.1, under the assumptions of diffusive shock acceleration of thermal particles in the test particle regime. Here, we present results from new Chandra observations of the cluster. The Chandra surface brightness profile across the N relic only hints to a surface brightness discontinuity (<2-sigma detection). Nevertheless, our reanalysis of archival Suzaku data indicates a temperature discontinuity across the relic that is consistent with a Mach number of 2.5+/-0.5, in agreement with previously published results. This confirms that the Mach number at the shock traced by the N relic is much weaker than predicted from the radio. Puzzlingly, in the Chandra data we also identify additional inner small density discontinuities both on and off the merger axis. Temperature measurements on both sides of the discontinuities do not allow us to undoubtedly determine their nature, although a shock front interpretation seems more likely. We speculate that if the inner density discontinuities are indeed shock fronts, then they are the consequence of violent relaxation of the dark matter cores of the clusters involved in the merger.

preprint2014arXiv

No Shock Across Part of a Radio Relic in the Merging Galaxy Cluster ZwCl 2341.1+0000?

The galaxy cluster ZwCl 2341.1+0000 is a merging system at z=0.27, which hosts two radio relics and a central, faint, filamentary radio structure. The two radio relics have unusually flat integrated spectral indices of -0.49 +/- 0.18 and -0.76 +/- 0.17, values that cannot be easily reconciled with the theory of standard diffusive shock acceleration of thermal particles at weak merger shocks. We present imaging results from XMM-Newton and Chandra observations of the cluster, aimed to detect and characterise density discontinuities in the ICM. As expected, we detect a density discontinuity near each of the radio relics. However, if these discontinuities are the shock fronts that fuelled the radio emission, then their Mach numbers are surprisingly low, both <=2. We studied the aperture of the density discontinuities, and found that while the NW discontinuity spans the whole length of the NW radio relic, the arc spanned by the SE discontinuity is shorter than the arc spanned by the SE relic. This startling result is in apparent contradiction with our current understanding of the origin of radio relics. Deeper X-ray data are required to confirm our results and to determine the nature of the density discontinuities.

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 LOFAR long baseline snapshot calibrator survey

Aims. An efficient means of locating calibrator sources for International LOFAR is developed and used to determine the average density of usable calibrator sources on the sky for subarcsecond observations at 140 MHz. Methods. We used the multi-beaming capability of LOFAR to conduct a fast and computationally inexpensive survey with the full International LOFAR array. Sources were pre-selected on the basis of 325 MHz arcminute-scale flux density using existing catalogues. By observing 30 different sources in each of the 12 sets of pointings per hour, we were able to inspect 630 sources in two hours to determine if they possess a sufficiently bright compact component to be usable as LOFAR delay calibrators. Results. Over 40% of the observed sources are detected on multiple baselines between international stations and 86 are classified as satisfactory calibrators. We show that a flat low-frequency spectrum (from 74 to 325 MHz) is the best predictor of compactness at 140 MHz. We extrapolate from our sample to show that the density of calibrators on the sky that are sufficiently bright to calibrate dispersive and non-dispersive delays for the International LOFAR using existing methods is 1.0 per square degree. Conclusions. The observed density of satisfactory delay calibrator sources means that observations with International LOFAR should be possible at virtually any point in the sky, provided that a fast and efficient search using the methodology described here is conducted prior to the observation to identify the best calibrator.

preprint2014arXiv

The LOFAR Pilot Surveys for Pulsars and Fast Radio Transients

We have conducted two pilot surveys for radio pulsars and fast transients with the Low-Frequency Array (LOFAR) around 140 MHz and here report on the first low-frequency fast-radio burst limit and the discovery of two new pulsars. The first survey, the LOFAR Pilot Pulsar Survey (LPPS), observed a large fraction of the northern sky, ~1.4 x 10^4 sq. deg, with 1-hr dwell times. Each observation covered ~75 sq. deg using 7 independent fields formed by incoherently summing the high-band antenna fields. The second pilot survey, the LOFAR Tied-Array Survey (LOTAS), spanned ~600 sq. deg, with roughly a 5-fold increase in sensitivity compared with LPPS. Using a coherent sum of the 6 LOFAR "Superterp" stations, we formed 19 tied-array beams, together covering 4 sq. deg per pointing. From LPPS we derive a limit on the occurrence, at 142 MHz, of dispersed radio bursts of < 150 /day/sky, for bursts brighter than S > 107 Jy for the narrowest searched burst duration of 0.66 ms. In LPPS, we re-detected 65 previously known pulsars. LOTAS discovered two pulsars, the first with LOFAR or any digital aperture array. LOTAS also re-detected 27 previously known pulsars. These pilot studies show that LOFAR can efficiently carry out all-sky surveys for pulsars and fast transients, and they set the stage for further surveying efforts using LOFAR and the planned low-frequency component of the Square Kilometer Array.

preprint2014arXiv

The shape of the radio wavefront of extensive air showers as measured with LOFAR

Extensive air showers, induced by high energy cosmic rays impinging on the Earth's atmosphere, produce radio emission that is measured with the LOFAR radio telescope. As the emission comes from a finite distance of a few kilometers, the incident wavefront is non-planar. A spherical, conical or hyperbolic shape of the wavefront has been proposed, but measurements of individual air showers have been inconclusive so far. For a selected high-quality sample of 161 measured extensive air showers, we have reconstructed the wavefront by measuring pulse arrival times to sub-nanosecond precision in 200 to 350 individual antennas. For each measured air shower, we have fitted a conical, spherical, and hyperboloid shape to the arrival times. The fit quality and a likelihood analysis show that a hyperboloid is the best parametrization. Using a non-planar wavefront shape gives an improved angular resolution, when reconstructing the shower arrival direction. Furthermore, a dependence of the wavefront shape on the shower geometry can be seen. This suggests that it will be possible to use a wavefront shape analysis to get an additional handle on the atmospheric depth of the shower maximum, which is sensitive to the mass of the primary particle.

preprint2013arXiv

Calibrating High-Precision Faraday Rotation Measurements for LOFAR and the Next Generation of Low-Frequency Radio Telescopes

Faraday rotation measurements using the current and next generation of low-frequency radio telescopes will provide a powerful probe of astronomical magnetic fields. However, achieving the full potential of these measurements requires accurate removal of the time-variable ionospheric Faraday rotation contribution. We present ionFR, a code that calculates the amount of ionospheric Faraday rotation for a specific epoch, geographic location, and line-of-sight. ionFR uses a number of publicly available, GPS-derived total electron content maps and the most recent release of the International Geomagnetic Reference Field. We describe applications of this code for the calibration of radio polarimetric observations, and demonstrate the high accuracy of its modeled ionospheric Faraday rotations using LOFAR pulsar observations. These show that we can accurately determine some of the highest-precision pulsar rotation measures ever achieved. Precision rotation measures can be used to monitor rotation measure variations - either intrinsic or due to the changing line-of-sight through the interstellar medium. This calibration is particularly important for nearby sources, where the ionosphere can contribute a significant fraction of the observed rotation measure. We also discuss planned improvements to ionFR, as well as the importance of ionospheric Faraday rotation calibration for the emerging generation of low-frequency radio telescopes, such as the SKA and its pathfinders.

preprint2013arXiv

Challenges to our understanding of radio relics: X-ray observations of the Toothbrush cluster

The cluster 1RXS J0603.3+4214 is a merging galaxy cluster that hosts three radio relics and a giant radio halo. The northern relic, the Toothbrush, is 1.9-Mpc long and has an unusual linear morphology. According to simple diffusive shock acceleration theory, its radio spectral index indicates a Mach number of 3.3-4.6. Here, we present results from a deep XMM-Newton observation of the cluster. We observe two distinct cluster cores that have survived the merger. The presence of three shocks at or near the locations of the radio relics is confirmed by density and temperature discontinuities. However, the observation poses several puzzles that challenge our understanding of radio relics: (i) at the Toothbrush, the shock Mach number is not larger than 2, in apparent conflict with the shock strength predicted from the radio spectrum; (ii) at the Toothbrush, the shock front is, in part, spatially offset from the radio emission; (iii) at the eastern relic, we detect a temperature jump corresponding to a Mach number of approximately 2.5, but there is no associated surface brightness discontinuity. We discuss possible explanations for these findings.

preprint2013arXiv

Complex diffuse radio emission in the merging PLANCK ESZ cluster Abell 3411

We present VLA radio and Chandra X-ray observations of the merging galaxy cluster Abell 3411. For the cluster, we find an overall temperature of 6.4 keV and an X-ray luminosity of 2.8 x 10^{44} erg s^{-1} between 0.5 and 2.0 keV. The Chandra observation reveals the cluster to be undergoing a merger event. The VLA observations show the presence of large-scale diffuse emission in the central region of the cluster, which we classify as a 0.9 Mpc size radio halo. In addition, a complex region of diffuse, polarized emission is found in the southeastern outskirts of the cluster, along the projected merger axis of the system. We classify this region of diffuse emission as a radio relic. The total extent of this radio relic is 1.9 Mpc. For the combined emission in the cluster region, we find a radio spectral index of -1.0 \pm 0.1 between 74 MHz and 1.4 GHz. The morphology of the radio relic is peculiar, as the relic is broken up into five fragments. This suggests that the shock responsible for the relic has been broken up due to interaction with a large-scale galaxy filament connected to the cluster or other substructures in the ICM. Alternatively, the complex morphology reflects the presence of electrons in fossil radio bubbles that are re-accelerated by a shock.

preprint2013arXiv

Detecting cosmic rays with the LOFAR radio telescope

The low frequency array (LOFAR), is the first radio telescope designed with the capability to measure radio emission from cosmic-ray induced air showers in parallel with interferometric observations. In the first $\sim 2\,\mathrm{years}$ of observing, 405 cosmic-ray events in the energy range of $10^{16} - 10^{18}\,\mathrm{eV}$ have been detected in the band from $30 - 80\,\mathrm{MHz}$. Each of these air showers is registered with up to $\sim1000$ independent antennas resulting in measurements of the radio emission with unprecedented detail. This article describes the dataset, as well as the analysis pipeline, and serves as a reference for future papers based on these data. All steps necessary to achieve a full reconstruction of the electric field at every antenna position are explained, including removal of radio frequency interference, correcting for the antenna response and identification of the pulsed signal.

preprint2013arXiv

Differential Frequency-dependent Delay from the Pulsar Magnetosphere

Some radio pulsars show clear drifting subpulses, in which subpulses are seen to drift in pulse longitude in a systematic pattern. Here we examine how the drifting subpulses of PSR B0809+74 evolve with time and observing frequency. We show that the subpulse period (P3) is constant on timescales of days, months and years, and between 14-5100 MHz. Despite this, the shapes of the driftbands change radically with frequency. Previous studies have concluded that, while the subpulses appear to move through the pulse window approximately linearly at low frequencies (< 500 MHz), a discrete step of 180 degrees in subpulse phase is observed at higher frequencies (> 820 MHz) near to the peak of the average pulse profile. We use LOFAR, GMRT, GBT, WSRT and Effelsberg 100-m data to explore the frequency-dependence of this phase step. We show that the size of the subpulse phase step increases gradually, and is observable even at low frequencies. We attribute the subpulse phase step to the presence of two separate driftbands, whose relative arrival times vary with frequency - one driftband arriving 30 pulses earlier at 20 MHz than it does at 1380 MHz, whilst the other arrives simultaneously at all frequencies. The drifting pattern which is observed here cannot be explained by either the rotating carousel model or the surface oscillation model, and could provide new insight into the physical processes happening within the pulsar magnetosphere.

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

Initial deep LOFAR observations of Epoch of Reionization windows: I. The North Celestial Pole

The aim of the LOFAR Epoch of Reionization (EoR) project is to detect the spectral fluctuations of the redshifted HI 21cm signal. This signal is weaker by several orders of magnitude than the astrophysical foreground signals and hence, in order to achieve this, very long integrations, accurate calibration for stations and ionosphere and reliable foreground removal are essential. One of the prospective observing windows for the LOFAR EoR project will be centered at the North Celestial Pole (NCP). We present results from observations of the NCP window using the LOFAR highband antenna (HBA) array in the frequency range 115 MHz to 163 MHz. The data were obtained in April 2011 during the commissioning phase of LOFAR. We used baselines up to about 30 km. With about 3 nights, of 6 hours each, effective integration we have achieved a noise level of about 100 microJy/PSF in the NCP window. Close to the NCP, the noise level increases to about 180 microJy/PSF, mainly due to additional contamination from unsubtracted nearby sources. We estimate that in our best night, we have reached a noise level only a factor of 1.4 above the thermal limit set by the noise from our Galaxy and the receivers. Our continuum images are several times deeper than have been achieved previously using the WSRT and GMRT arrays. We derive an analytical explanation for the excess noise that we believe to be mainly due to sources at large angular separation from the NCP.

preprint2013arXiv

LOFAR detections of low-frequency radio recombination lines towards Cassiopeia A

Cassiopeia A was observed using the Low-Band Antennas of the LOw Frequency ARray (LOFAR) with high spectral resolution. This allowed a search for radio recombination lines (RRLs) along the line-of-sight to this source. Five carbon-alpha RRLs were detected in absorption between 40 and 50 MHz with a signal-to-noise ratio of > 5 from two independent LOFAR datasets. The derived line velocities (v_LSR ~ -50 km/s) and integrated optical depths (~ 13 s^-1) of the RRLs in our spectra, extracted over the whole supernova remnant, are consistent within each LOFAR dataset and with those previously reported. For the first time, we are able to extract spectra against the brightest hotspot of the remnant at frequencies below 330 MHz. These spectra show significantly higher (15-80 %) integrated optical depths, indicating that there is small-scale angular structure on the order of ~1 pc in the absorbing gas distribution over the face of the remnant. We also place an upper limit of 3 x 10^-4 on the peak optical depths of hydrogen and helium RRLs. These results demonstrate that LOFAR has the desired spectral stability and sensitivity to study faint recombination lines in the decameter band.

preprint2013arXiv

Studying Galactic interstellar turbulence through fluctuations in synchrotron emission: First LOFAR Galactic foreground detection

The characteristic outer scale of turbulence and the ratio of the random to ordered components of the magnetic field are key parameters to characterise magnetic turbulence in the interstellar gas, which affects the propagation of cosmic rays within the Galaxy. We provide new constraints to those two parameters. We use the LOw Frequency ARray (LOFAR) to image the diffuse continuum emission in the Fan region at (l,b) (137.0,+7.0) at 80"x70" resolution in the range [146,174] MHz. We detect multi-scale fluctuations in the Galactic synchrotron emission and compute their power spectrum. Applying theoretical estimates and derivations from the literature for the first time, we derive the outer scale of turbulence and the ratio of random to ordered magnetic field from the characteristics of these fluctuations . We obtain the deepest image of the Fan region to date and find diffuse continuum emission within the primary beam. The power spectrum of the foreground synchrotron fluctuations displays a power law behaviour for scales between 100 and 8 arcmin with a slope of (-1.84+/-0.19). We find an upper limit of about 20 pc for the outer scale of the magnetic interstellar turbulence toward the Fan region. We also find a variation of the ratio of random to ordered field as a function of Galactic coordinates, supporting different turbulent regimes. We use power spectra fluctuations from LOFAR as well as earlier GMRT and WSRT observations to constrain the outer scale of turbulence of the Galactic synchrotron foreground, finding a range of plausible values of 10-20 pc. Then, we use this information to deduce lower limits of the ratio of ordered to random magnetic field strength. These are found to be 0.3, 0.3, and 0.5 for the LOFAR, WSRT and GMRT fields considered respectively. Both these constraints are in agreement with previous estimates.

preprint2013arXiv

The brightness and spatial distributions of terrestrial radio sources

Faint undetected sources of radio-frequency interference (RFI) might become visible in long radio observations when they are consistently present over time. Thereby, they might obstruct the detection of the weak astronomical signals of interest. This issue is especially important for Epoch of Reionisation (EoR) projects that try to detect the faint redshifted HI signals from the time of the earliest structures in the Universe. We explore the RFI situation at 30-163 MHz by studying brightness histograms of visibility data observed with LOFAR, similar to radio-source-count analyses that are used in cosmology. An empirical RFI distribution model is derived that allows the simulation of RFI in radio observations. The brightness histograms show an RFI distribution that follows a power-law distribution with an estimated exponent around -1.5. With several assumptions, this can be explained with a uniform distribution of terrestrial radio sources whose radiation follows existing propagation models. Extrapolation of the power law implies that the current LOFAR EoR observations should be severely RFI limited if the strength of RFI sources remains strong after time integration. This is in contrast with actual observations, which almost reach the thermal noise and are thought not to be limited by RFI. Therefore, we conclude that it is unlikely that there are undetected RFI sources that will become visible in long observations. Consequently, there is no indication that RFI will prevent an EoR detection with LOFAR.

preprint2013arXiv

The nature of filamentary cold gas in the core of the Virgo Cluster

We present a multi-wavelength study of the emission-line nebulae located southeast of the nucleus of M87, the central dominant galaxy of the Virgo Cluster. We report the detection of far-infrared (FIR) [CII] line emission from the nebulae using observations made with Herschel PACS. The infrared line emission is extended and cospatial with optical Hα+[NII], far-ultraviolet CIV lines, and soft X-ray emission. The filamentary nebulae evidently contain multi-phase material spanning a temperature range of at least 5 orders of magnitude, from ~100 K to ~10^7 K. This material has most likely been uplifted by the AGN from the center of M87. The thermal pressure of the 10^4 K phase appears to be significantly lower than that of the surrounding hot intra-cluster medium (ICM) indicating the presence of additional turbulent and magnetic pressure in the filaments. If the turbulence in the filaments is subsonic then the magnetic field strength required to balance the pressure of the surrounding ICM is B~30-70 μG. The spectral properties of the soft X-ray emission from the filaments indicate that it is due to thermal plasma with kT~0.5-1 keV, which is cooling by mixing with the cold gas and/or radiatively. Charge exchange can be ruled out as a significant source of soft X-rays. Both cooling and mixing scenarios predict gas with a range of temperatures. This is at first glance inconsistent with the apparent lack of X-ray emitting gas with kT<0.5 keV. However, we show that the missing very soft X-ray emission could be absorbed by the cold gas in the filaments with an integrated hydrogen column density of ~1.6x10^21 cm^-2, providing a natural explanation for the apparent temperature floor to the X-ray emission at kT~0.5 keV. The FIR through ultra-violet line emission is most likely primarily powered by the ICM particles penetrating the cold gas following a shearing induced mixing process.

preprint2012arXiv

Diffuse radio emission in MACS J1752.0+4440

We report the discovery of large-scale diffuse radio emission in the galaxy cluster MACS J1752.0+4440 (z=0.366). Using Westerbork Synthesis Radio Telescope (WSRT) observations we find that the cluster hosts a double radio relic system as well as a 1.65 Mpc radio halo covering the region between the two relics. The relics are diametrically located on opposite sides of the cluster center. The NE and SW relics have sizes of 1.3 and 0.9 Mpc, respectively. In case of an isolated binary merger event, the relative sizes of the relics suggest a mass ratio about 2:1. We measure integrated spectra of -1.16 \pm 0.03 for the NE and -1.10 \pm 0.05 for the SW relic. We conclude that this cluster has undergone a violent binary merger event and the relics are best explained by particles (re)accelerated in outwards traveling shock waves. The spectral indices suggest the relics trace shock waves with Mach numbers (M) of around 3.5 to 4.5. These relatively high Mach numbers derived from the radio spectral index are comparable to those derived for a few other recently discovered relics. This implies that merger shocks with M > 3 are relatively common in cluster outskirts if our understanding of diffusive shock acceleration is correct.

preprint2012arXiv

First LOFAR observations at very low frequencies of cluster-scale non-thermal emission: the case of Abell 2256

Abell 2256 is one of the best known examples of a galaxy cluster hosting large-scale diffuse radio emission that is unrelated to individual galaxies. It contains both a giant radio halo and a relic, as well as a number of head-tail sources and smaller diffuse steep-spectrum radio sources. The origin of radio halos and relics is still being debated, but over the last years it has become clear that the presence of these radio sources is closely related to galaxy cluster merger events. Here we present the results from the first LOFAR Low band antenna (LBA) observations of Abell 2256 between 18 and 67 MHz. To our knowledge, the image presented in this paper at 63 MHz is the deepest ever obtained at frequencies below 100 MHz in general. Both the radio halo and the giant relic are detected in the image at 63 MHz, and the diffuse radio emission remains visible at frequencies as low as 20 MHz. The observations confirm the presence of a previously claimed ultra-steep spectrum source to the west of the cluster center with a spectral index of -2.3 \pm 0.4 between 63 and 153 MHz. The steep spectrum suggests that this source is an old part of a head-tail radio source in the cluster. For the radio relic we find an integrated spectral index of -0.81 \pm 0.03, after removing the flux contribution from the other sources. This is relatively flat which could indicate that the efficiency of particle acceleration at the shock substantially changed in the last \sim 0.1 Gyr due to an increase of the shock Mach number. In an alternative scenario, particles are re-accelerated by some mechanism in the downstream region of the shock, resulting in the relatively flat integrated radio spectrum. In the radio halo region we find indications of low-frequency spectral steepening which may suggest that relativistic particles are accelerated in a rather inhomogeneous turbulent region.

preprint2012arXiv

M87 at metre wavelengths: the LOFAR picture

M87 is a giant elliptical galaxy located in the centre of the Virgo cluster, which harbours a supermassive black hole of mass 6.4x10^9 M_sun, whose activity is responsible for the extended (80 kpc) radio lobes that surround the galaxy. The energy generated by matter falling onto the central black hole is ejected and transferred to the intra-cluster medium via a relativistic jet and morphologically complex systems of buoyant bubbles, which rise towards the edges of the extended halo. Here we present the first observations made with the new Low-Frequency Array (LOFAR) of M87 at frequencies down to 20 MHz. Images of M87 were produced at low radio frequencies never explored before at these high spatial resolution and dynamic range. To disentangle different synchrotron models and place constraints on source magnetic field, age and energetics, we also performed a detailed spectral analysis of M87 extended radio-halo using these observations together with archival data. We do not find any sign of new extended emissions; on the contrary the source appears well confined by the high pressure of the intra-cluster medium. A continuous injection of relativistic electrons is the model that best fits our data, and provides a scenario in which the lobes are still supplied by fresh relativistic particles from the active galactic nuclei. We suggest that the discrepancy between the low-frequency radio-spectral slope in the core and in the halo implies a strong adiabatic expansion of the plasma as soon as it leaves the core area. The extended halo has an equipartition magnetic field strength of ~10 uG, which increases to ~13 uG in the zones where the particle flows are more active. The continuous injection model for synchrotron ageing provides an age for the halo of ~40 Myr, which in turn provides a jet kinetic power of 6-10x10^44 erg/s.

preprint2012arXiv

Radio relics in cosmological simulations

Radio relics have been discovered in many galaxy clusters. They are believed to trace shock fronts induced by cluster mergers. Cosmological simulations allow us to study merger shocks in detail since the intra-cluster medium is heated by shock dissipation. Using high resolution cosmological simulations, identifying shock fronts and applying a parametric model for the radio emission allows us to simulate the formation of radio relics. We analyze a simulated shock front in detail. We find a rather broad Mach number distribution. The Mach number affects strongly the number density of relativistic electrons in the downstream area, hence, the radio luminosity varies significantly across the shock surface. The abundance of radio relics can be modeled with the help of the radio power probability distribution which aims at predicting radio relic number counts. Since the actual electron acceleration efficiency is not known, predictions for the number counts need to be normalized by the observed number of radio relics. For the characteristics of upcoming low frequency surveys we find that about thousand relics are awaiting discovery.

preprint2012arXiv

Simulating the Toothbrush: Evidence for a triple merger of galaxy clusters

The newly discovered galaxy cluster 1RXS J0603.3+4214 hosts a 1.9 Mpc long, bright radio relic with a peculiar linear morphology. Using hydrodynamical +N-body AMR simulations of the merger between three initially hydrostatic clusters in an idealised setup, we are able to reconstruct the morphology of the radio relic. Based on our simulation, we can constrain the merger geometry, predict lensing mass measurements and X-ray observations. Comparing such models to X-ray, redshift and lensing data will validate the geometry of this complex merger which helps to constrain the parameters for shock acceleration of electrons that produces the radio relic.

preprint2012arXiv

The "toothbrush-relic": evidence for a coherent linear 2-Mpc scale shock wave in a massive merging galaxy cluster?

Some merging galaxy clusters host diffuse extended radio emission, so-called radio halos and relics. Here we present observations between 147 MHz and 4.9 GHz of a new radio-selected galaxy cluster 1RXS J0603.3+4214 (z=0.225). The cluster is also detected as an extended X-ray source in the RASS. It hosts a large bright 1.9 Mpc radio relic, an elongated ~2 Mpc radio halo, and two smaller radio relics. The large radio relic has a peculiar linear morphology. For this relic we observe a clear spectral index gradient, in the direction towards the cluster center. We performed Rotation Measure (RM) Synthesis between 1.2 and 1.7 GHz. The results suggest that for the west part of the large relic some of the Faraday rotation is caused by ICM and is not only due to galactic foregrounds. We also carried out a detailed spectral analysis of this radio relic and created radio color-color diagrams. We find (i) an injection spectral index of -0.6 to -0.7, (ii) steepening spectral index and increasing spectral curvature in the post-shock region, and (iii) an overall power-law spectrum between 74 MHz and 4.9 GHz with α=-1.10 \pm 0.02. Mixing of emission in the beam from regions with different spectral ages is probably the dominant factor that determines the shape of the radio spectra. Changes in the magnetic field, total electron content, or adiabatic gains/losses do not play a major role. A model in which particles are (re)accelerated in a first order Fermi process at the front of the relic provides the best match to the observed spectra. We speculate that in the post-shock region particles are re-accelerated by merger induced turbulence to form the radio halo as the relic and halo are connected. The 1RXS J0603.3+4214 merger is probably more complex than the "simple'" binary merger events that are thought to give rise to symmetric double radio relics.

preprint2012arXiv

The LOFAR radio environment

Aims: This paper discusses the spectral occupancy for performing radio astronomy with the Low-Frequency Array (LOFAR), with a focus on imaging observations. Methods: We have analysed the radio-frequency interference (RFI) situation in two 24-h surveys with Dutch LOFAR stations, covering 30-78 MHz with low-band antennas and 115-163 MHz with high-band antennas. This is a subset of the full frequency range of LOFAR. The surveys have been observed with a 0.76 kHz / 1 s resolution. Results: We measured the RFI occupancy in the low and high frequency sets to be 1.8% and 3.2% respectively. These values are found to be representative values for the LOFAR radio environment. Between day and night, there is no significant difference in the radio environment. We find that lowering the current observational time and frequency resolutions of LOFAR results in a slight loss of flagging accuracy. At LOFAR's nominal resolution of 0.76 kHz and 1 s, the false-positives rate is about 0.5%. This rate increases approximately linearly when decreasing the data frequency resolution. Conclusions: Currently, by using an automated RFI detection strategy, the LOFAR radio environment poses no perceivable problems for sensitive observing. It remains to be seen if this is still true for very deep observations that integrate over tens of nights, but the situation looks promising. Reasons for the low impact of RFI are the high spectral and time resolution of LOFAR; accurate detection methods; strong filters and high receiver linearity; and the proximity of the antennas to the ground. We discuss some strategies that can be used once low-level RFI starts to become apparent. It is important that the frequency range of LOFAR remains free of broadband interference, such as DAB stations and windmills.

preprint2012arXiv

Wide-band Simultaneous Observations of Pulsars: Disentangling Dispersion Measure and Profile Variations

Dispersion in the interstellar medium is a well known phenomenon that follows a simple relationship, which has been used to predict the time delay of dispersed radio pulses since the late 1960s. We performed wide-band simultaneous observations of four pulsars with LOFAR (at 40-190 MHz), the 76-m Lovell Telescope (at 1400 MHz) and the Effelsberg 100-m Telescope (at 8000 MHz) to test the accuracy of the dispersion law over a broad frequency range. In this paper we present the results of these observations which show that the dispersion law is accurate to better than 1 part in 100000 across our observing band. We use this fact to constrain some of the properties of the ISM along the line-of-sight and use the lack of any aberration or retardation effects to determine upper limits on emission heights in the pulsar magnetosphere. We also discuss the effect of pulse profile evolution on our observations, and the implications that it could have for precision pulsar timing projects such as the detection of gravitational waves with pulsar timing arrays.

preprint2012arXiv

X-ray observations of the merging cluster CIZA J2242.8+5301

Multiwavelength studies of radio relics at merger shocks set powerful constraints on the relics origin and formation mechanism. However, for X-ray observations, a main difficulty is represented by the low X-ray surface brightness far out in the cluster outskirts, where relics are typically found. Here, we present XMM-Newton results from a 130-ks observation of CIZA J2242.8+5301, a cluster at z=0.19 that hosts a double radio relic. We focus on the well-defined northern relic. There is a difference of ~55% between the temperature we measure behind the relic, and the temperature measured with Suzaku. We analyse the reasons for this large discrepancy, and discuss the possibility of reliably measuring the temperature beyond the northern relic.

preprint2011arXiv

A deep radio and X-ray view of cluster formation at the crossroads of filaments

Deep X-ray data from CHANDRA and XMM-Newton, and GMRT radio data are presented for ZwCl 2341.1+0000, an extremely unusual and complex merging cluster of galaxies at the intersection of optical filaments. We propose that energetics of multiple mergers and accretion flows has resulted in wide-spread shocks, acceleration of cosmic ray particles and amplification of weak magnetic fields. This results in Mpc-scale peripheral radio relics and halo like non-thermal emission observed near the merging center.

preprint2011arXiv

A double radio relic in the merging galaxy cluster ZwCl 0008.8+5215

Some merging galaxy clusters host diffuse elongated radio sources, also called radio relics. It is proposed that these radio relics trace shock waves in the intracluster medium created during a cluster merger event. Within the shock waves particles are accelerated to relativistic energies, and in the presence of a magnetic field synchrotron radiation will be emitted. Here we present GMRT and WSRT observations of a new double relic in the galaxy cluster ZwCl 0008.8+5215. Optical V, R, and I band images of the cluster were taken with the INT. An optical spectrum, to determine the redshift of the cluster, was taken with the WHT. Our observations show the presence of a double radio relic in the galaxy cluster ZwCl 0008.8+5215, for which we find a spectroscopic redshift of z = 0.1032 \pm 0.0018 from an optical spectrum of one of the cD galaxies. The spectral index of the two relics steepens inwards to the cluster center. For part of the relics, we measure a polarization fraction in the range ~ 5-25%. A ROSAT X-ray image displays an elongated ICM and the large-scale distribution of galaxies reveals two cluster cores, all pointing towards a binary cluster merger event. The radio relics are located symmetrically with respect to the X-ray center of the cluster, along the proposed merger axis. The relics have a linear extent of 1.4 Mpc and 290 kpc. This factor of five difference in linear size is unlike that of previously known double relic systems, for which the sizes do not differ by more than a factor of two. We conclude that the double relics in ZwCl 0008.8+5215 are best explained by two outward moving shock waves in which particles are (re)accelerated trough the diffusive shock acceleration (DSA) mechanism. [abridged]

preprint2011arXiv

A merger mystery: no extended radio emission in the merging cluster Abell 2146

We present a new 400ks Chandra X-ray observation and a GMRT radio observation at 325MHz of the merging galaxy cluster Abell 2146. The Chandra observation reveals detailed structure associated with the major merger event including the Mach M=2.1+/-0.2 bow shock located ahead of the dense subcluster core and the first known example of an upstream shock (M=1.6+/-0.1). Surprisingly, the deep GMRT observation at 325MHz does not detect any extended radio emission associated with either shock front. All other merging galaxy clusters with X-ray detected shock fronts, including the Bullet cluster, Abell 520, Abell 754 and Abell 2744, and clusters with candidate shock fronts have detected radio relics or radio halo edges coincident with the shocks. We consider several possible factors which could affect the formation of radio relics, including the shock strength and the presence of a pre-existing electron population, but do not find a favourable explanation for this result. We calculate a 3sigma upper limit of 13mJy on extended radio emission, which is significantly below the radio power expected by the observed P_{radio}-L_{X} correlation for merging systems. The lack of an extended radio halo in Abell 2146 maybe due to the low cluster mass relative to the majority of merging galaxy clusters with detected radio halos.

preprint2011arXiv

Deep Low-Frequency Radio Observations of the NOAO Bootes Field: I. Data Reduction and Catalog Construction

In this article we present deep, high-resolution radio interferometric observations at 153 MHz to complement the extensively studied NOAO Bootes field. We provide a description of the observations, data reduction and source catalog construction. From our single pointing GMRT observation of ~12 hours we obtain a high-resolution (26" x 22") image of ~11.3 square degrees, fully covering the Bootes field region and beyond. The image has a central noise level of ~1.0 mJy/beam, which rises to 2.0-2.5 mJy/beam at the field edge, placing it amongst the deepest ~150 MHz surveys to date. The catalog of 598 extracted sources is estimated to be ~92 percent complete for >10 mJy sources, while the estimated contamination with false detections is <1 percent. The low RMS position uncertainty of 1.24" facilitates accurate matching against catalogs at optical, infrared and other wavelengths. Differential source counts are determined down to <~10 mJy. There is no evidence for flattening of the counts towards lower flux densities as observed in deep radio surveys at higher frequencies, suggesting that our catalog is dominated by the classical radio-loud AGN population that explains the counts at higher flux densities. Combination with available deep 1.4 GHz observations yields an accurate determination of spectral indices for 417 sources down to the lowest 153 MHz flux densities, of which 16 have ultra-steep spectra with spectral indices below -1.3. We confirm that flattening of the median spectral index towards low flux densities also occurs at this frequency. The detection fraction of the radio sources in NIR Ks-band is found to drop with radio spectral index, which is in agreement with the known correlation between spectral index and redshift for brighter radio sources.

preprint2011arXiv

Diffuse steep-spectrum sources from the 74 MHz VLSS survey

Galaxy clusters grow by a sequence of mergers with other clusters and galaxy groups. During these mergers, shocks and/or turbulence are created within the intracluster medium (ICM). In this process, particles could be accelerated to highly relativistic energies. The synchrotron radiation from these particles is observed in the form of radio relics and halos that are generally characterized by a steep radio spectral index. Shocks can also revive fossil radio plasma from a previous episode of AGN activity, creating a so-called radio phoenix. Here we present multi-frequency radio observations of diffuse steep-spectrum radio sources selected from the 74 MHz VLSS survey. Previous GMRT observations showed that some of these sources had filamentary and elongated morphologies, which are expected for radio relics. We carried out radio continuum observations at 325 MHz with the GMRT. Observations with the VLA and WSRT were taken at 1.4 GHz in full polarization mode. Optical images around the radio sources were taken with the WHT and INT telescopes. Most of the sources in our sample consist of old radio plasma from AGNs located in small galaxy clusters. The sources can be classified as AGN relics or radio phoenices. The spectral indices across most of the radio sources display large variations. We conclude that diffuse steep-spectrum radio sources are not only found in massive X-ray luminous galaxy clusters but also in smaller systems. Future low-frequency surveys will uncover large numbers of steep-spectrum radio relics related to previous episodes of AGN activity. [abridged]

preprint2011arXiv

Far Ultraviolet Emission in the A2597 and A2204 Brightest Cluster Galaxies

We use the Hubble Space Telescope ACS/SBC and Very Large Telescope FORS cameras to observe the Brightest Cluster Galaxies in Abell 2597 and Abell 2204 in the far-ultraviolet (FUV) F150LP and optical U, B, V, R, I Bessel filters. The FUV and U band emission is enhanced in bright, filamentary structures surrounding the BCG nuclei. These filaments can be traced out to 20 kpc from the nuclei in the FUV. Excess FUV and U band light is determined by removing emission due to the underlying old stellar population and mapped with 1 arcsec spatial resolution over the central 20 kpc regions of both galaxies. We find the FUV and U excess emission to be spatially coincident and a stellar interpretation requires the existence of a significant amount of 10000-50000 K stars. Correcting for nebular continuum emission and dust intrinsic to the BCG further increases the FUV to U band emission ratio and implies that stars alone may not suffice to explain the observations. However, lack of detailed information on the gas and dust distribution and extinction law in these systems prevents us from ruling out a purely stellar origin. Non-stellar processes, such as the central AGN, Scattering, Synchrotron and Bremsstrahlung emission are investigated and found to not be able to explain the FUV and U band measurements in A2597. Contributions from non-thermal processes not treated here should be investigated. Comparing the FUV emission to the optical H-alpha line emitting nebula shows good agreement on kpc-scales in both A2597 and A2204. In concordance with an earlier investigation by O'Dea et al. (2004) we find that O-stars can account for the ionising photons necessary to explain the observed H-alpha line emission.

preprint2011arXiv

How many radio relics await discovery?

Upcoming radio telescopes will allow to study the radio sky at low frequencies with unprecedent sensitivity and resolution. New surveys are expected to discover a large number of new radio sources. Here we investigate the abundance of radio relics, i.e. steep-spectrum diffuse radio emission coming from the periphery of galaxy clusters, which are believed to trace shock waves induced by cluster mergers. With the advent of comprehensive relic samples a framework is needed to analyze statistically the relic abundance. To this end, we introduce the probability to find a relic located in a galaxy cluster with given mass and redshift allowing us to relate the halo mass function of the Universe with the radio relic number counts. Up to date about 45 relics have been reported and we compile the resulting counts, N(>S_1.4). In principle, the parameters of the distribution could be determined using a sufficiently large relic sample. However, since the number of known relics is still small for that purpose we use the MareNostrum Universe simulation to determine the relic radio power scaling with cluster mass and redshift. Our model is able to reproduce the recently found tentative evidence for an increase in the fraction of clusters hosting relics, both with X-ray luminosity and redshift, using an X-ray flux limited cluster sample. Moreover, we find that a considerable fraction of faint relics (S_1.4 < ~10 mJy) reside in clusters with an X-ray flux below ~3e-12 erg/s/cm^2. Finally, we estimate the number of radio relics which await discovery by future low frequency surveys proposed for LOFAR and WSRT. We estimate that the WODAN survey proposed for WSRT may discover 900 relics and that the LOFAR-Tier 1-120 MHz survey may discover about 2500 relics. However, the actual number of newly discovered relics will crucially depend on the existence of sufficiently complete galaxy cluster catalogues.

preprint2011arXiv

Low-frequency radio observations of the galaxy cluster CIZA J2242.8+5301

Some disturbed galaxy clusters host diffuse elongated radio sources, also called radio relics. It is proposed that these relics trace shock waves in the intracluster medium (ICM). Within the shock waves, generated by cluster merger events, particles are accelerated to relativistic energies, and in the presence of a magnetic field synchrotron radiation will be emitted. CIZA J2242.8+5301 is a disturbed galaxy cluster hosting complex diffuse radio emission, including a so-called double radio relic. Here we present new Giant Metrewave Radio Telescope (GMRT) radio observations of CIZA J2242.8+5301 at 325 and 150 MHz. We detect the double radio relic at 150 and 325 MHz. The very deep 150 MHz image reveals the presence of large-scale diffuse emission between the two radio relics.

preprint2011arXiv

Magnetic fields and shock waves in cluster outskirts

Both, X-ray and radio observations, are beginning to probe regions in galaxy clusters close to the virial radius. These are particularly interesting regions as they are (i) crucial for the thermal evolution of the intra-cluster medium and (ii) they host plasmas under conditions that are poorly understood. Recently discovered giant radio relics in the galaxy cluster CIZA J2242.8+5301 provide unprecedented evidence for particle acceleration in cluster outskirts. We will discuss some of the intriguing puzzles that surround these relics. The polarisation of the radio emission yields information about the topology of magnetic fields in these shock waves. Here we report on MHD simulations that highlight some problems about the origin of these magnetic fields.

preprint2011arXiv

Radio continuum observations of new radio halos and relics from the NVSS and WENSS surveys

Radio halos and relics are diffuse radio sources found in galaxy clusters showing significant substructure at X-ray wavelengths. These sources provide important information about non-thermal processes taking place in the ICM. Until now only a few dozen relics and halos are known, while models predict that a much larger number of these sources exist. In this paper we present the results of an extensive observing campaign to search for new diffuse radio sources in galaxy clusters. We discovered 6 new radio relics, including a probable double relic system, and 2 radio halos. In addition, we confirm the presence of diffuse radio emission in four galaxy clusters. By constructing a sample of 35 radio relics we find that relics are mostly found along the major axis of the X-ray emission from the ICM, while their orientation is perpendicular to this axis. This is consistent with the scenario that relics trace merger shock waves. The major and minor axis ratio distribution of the ICM for clusters with relics is also broader than that of the NORAS-REFLEX sample. We compared the X-ray luminosity and redshift distributions of clusters with relics to an X-ray selected sample from the NORAS and REFLEX surveys. We find tentative evidence for an increase of the cluster's relic fraction with X-ray luminosity and redshift. [abridged]

preprint2011arXiv

Radio relics in the MareNostrum Universe

We identify shocked gas in simulated galaxy clusters extracted from the MareNostrum Universe simulation (Gottloeber et al. 2006) assuming that shock waves are regions of electron acceleration. We perform flux number counts within the framework of the non-thermal emission model developed by Hoeft et al. (2008). Results are presented at two different observing frequencies, i.e. 1.4 GHz and 120 MHz, posing interesting constraints for LOFAR and upcoming radio telescopes.

preprint2011arXiv

Relics as probes of galaxy cluster mergers

Galaxy clusters grow by mergers with other clusters and galaxy groups. These mergers create shocks within the intracluster medium (ICM). It is proposed that within the shocks particles can be accelerated to extreme energies. In the presence of a magnetic field these particles should then form large regions emitting synchrotron radiation, creating so-called radio relics. An example of a cluster with relics is CIZA J2242.8+5301. Here we present hydrodynamical simulations of idealized binary cluster collisions with the aim of constraining the merger scenario for this cluster. We conclude that by using the location, size and width of double radio relics we can set constraints on the mass ratios, impact parameters, timescales, and viewing geometries of binary cluster merger events.

preprint2011arXiv

Using double radio relics to constrain galaxy cluster mergers: A model of double radio relics in CIZA J2242.8+5301

Galaxy clusters grow by mergers with other clusters and galaxy groups. These mergers create shock waves within the intracluster medium (ICM) that can accelerate particles to extreme energies. In the presence of magnetic fields, relativistic electrons form large regions emitting synchrotron radiation, so-called radio relics. Behind the shock front, synchrotron and inverse Compton (IC) losses cause the radio spectral index to steepen away from the shock front. An example of such a cluster is CIZA J2242.8+5301, where very clear spectral steepening in the downstream region is observed. Here we present hydrodynamical simulations of idealized binary cluster mergers with the aim of constraining the merger scenario for this cluster. From our simulations, we find that CIZA J2242.8+5301 is probably undergoing a merger in the plane of the sky (less then 10 deg from edge-on) with a mass ratio of about 2:1, and an impact parameter < 400 kpc. We find that the core passage of the clusters happened about 1 Gyr ago. We conclude that double relics relics can set constraints on the mass ratios, impact parameters, timescales, and viewing geometry of binary cluster mergers, which is particularly useful when detailed X-ray observations are not available. In addition, the presence of large radio relics can be used to constrain the degree of clumping in the outskirts of the ICM, which is important to constrain the baryon fraction, density and entropy profiles, around the virial radius and beyond. We find that the amplitude of density fluctuations, with sizes of < 200 kpc, in the relic in CIZA J2242.8+5301 is not larger than 30%. [abridged]

preprint2010arXiv

Particle Acceleration on Megaparsec Scales in a Merging Galaxy Cluster

Galaxy clusters form through a sequence of mergers of smaller galaxy clusters and groups. Models of diffusive shock acceleration (DSA) suggest that in shocks that occur during cluster mergers, particles are accelerated to relativistic energies, similar to supernova remnants. Together with magnetic fields these particles emit synchrotron radiation and may form so-called radio relics. Here we report the detection of a radio relic for which we find highly aligned magnetic fields, a strong spectral index gradient, and a narrow relic width, giving a measure of the magnetic field in an unexplored site of the universe. Our observations prove that DSA also operates on scales much larger than in supernova remnants and that shocks in galaxy clusters are capable of producing extremely energetic cosmic rays.

preprint2010arXiv

The Distribution and Condition of the Warm Molecular Gas in Abell 2597 and Sersic 159-03

We have used the SINFONI integral field spectrograph to map the near-infrared K-band emission lines of molecular and ionised hydrogen in the central regions of two cool core galaxy clusters, Abell 2597 and Sersic 159-03. Gas is detected out to 20 kpc from the nuclei of the brightest cluster galaxies and found to be distributed in clumps and filaments around it. The ionised and molecular gas phases trace each other closely in extent and dynamical state. Both gas phases show signs of interaction with the active nucleus. Within the nuclear regions the kinetic luminosity of this gas is found to be somewhat smaller than the current radio luminosity. Outside the nuclear region the gas has a low velocity dispersion and shows smooth velocity gradients. There is no strong correlation between the intensity of the molecular and ionised gas emission and either the radio or X-ray emission. The molecular gas in Abell 2597 and Sersic 159-03 is well described by a gas in local thermal equilibrium (LTE) with a single excitation temperature T_exc ~ 2300 K. The emission line ratios do not vary strongly as function of position, with the exception of the nuclear regions where the ionised to molecular gas ratio is found decrease. These constant line ratios imply a single source of heating and excitation for both gas phases.

preprint2009arXiv

A search for steep spectrum radio relics and halos with the GMRT

Context: Diffuse radio emission, in the form of radio halos and relics, traces regions in clusters with shocks or turbulence, probably produced by cluster mergers. Some models of diffuse radio emission in clusters indicate that virtually all clusters should contain diffuse radio sources with a steep spectrum. External accretion shocks associated with filamentary structures of galaxies could also accelerate electrons to relativistic energies and hence produce diffuse synchrotron emitting regions. Here we report on Giant Metrewave Radio Telescope (GMRT) observations of a sample of steep spectrum sources from the 74 MHz VLSS survey. These sources are diffuse and not associated with nearby galaxies. Aims: The main aim of the observations is to search for diffuse radio emission associated with galaxy clusters or the cosmic web. Methods: We carried out GMRT 610 MHz continuum observations of unidentified diffuse steep spectrum sources. Results: We have constructed a sample of diffuse steep spectrum sources, selected from the 74 MHz VLSS survey. We identified eight diffuse radio sources probably all located in clusters. We found five radio relics, one cluster with a giant radio halo and a radio relic, and one radio mini-halo. By complementing our observations with measurements from the literature we find correlations between the physical size of relics and the spectral index, in the sense that smaller relics have steeper spectra. Furthermore, larger relics are mostly located in the outskirts of clusters while smaller relics are located closer to the cluster center.

preprint2009arXiv

Diffuse radio emission in the merging cluster MACS J0717.5+3745: the discovery of the most powerful radio halo

Hierarchical models of structure formation predict that galaxy clusters grow via mergers of smaller clusters and galaxy groups, as well as through continuous accretion of gas. MACS J0717.5+3745 is an X-ray luminous and complex merging cluster, located at a redshift of 0.55. Here we present Giant Metrewave Radio Telescope (GMRT) radio observations at 610 MHz of this cluster. The main aim of the observations is to search for diffuse radio emission within the galaxy cluster MACS J0717.5+3745 related to the ongoing merger. These GMRT observations are complemented by Very Large Array (VLA) archival observations at 1.4, 4.9 and 8.5 GHz. We have discovered a radio halo in the cluster MACS J0717.5+3745 with a size of about 1.2 Mpc. The radio power P_1.4 GHz is 5 x 10^25 W/Hz, which makes it the most powerful radio halo known till date. A 700 kpc radio structure, which we classify as a radio relic, is located in between the merging substructures of the system. The location of this relic roughly coincides with regions of the intra-cluster medium (ICM) that have a significant enhancement in temperature as shown by Chandra. The major axis of the relic is also roughly perpendicular to the merger axis. This shows that the relic might be the result of a merger-related shock wave, where particles are accelerated via the diffuse shock acceleration (DSA) mechanism. Alternatively, the relic might trace an accretion shock of a large-scale galaxy filament to the south-west. The global spectral index of radio emission within the cluster is found to be -1.24 +/-0.05 between 4.9 GHz and 610 MHz. We derive a value of 5.8 microGauss for the equipartition magnetic field strength at the location of the radio halo. [abridged].

preprint2009arXiv

Radio observations of ZwCl 2341.1+0000: a double radio relic cluster

Context: Hierarchal models of large scale structure (LSS) formation predict that galaxy clusters grow via gravitational infall and mergers of (smaller) mass concentrations, such as clusters and galaxy groups. Diffuse radio emission, in the form of radio halos and relics, is found in clusters undergoing a merger, indicating that shocks or turbulence associated with the merger are capable of accelerating electrons to highly relativistic energies. Here we report on radio observations of ZwCl 2341.1+0000, a complex merging structure of galaxies located at z=0.27, using Giant Metrewave Radio Telescope (GMRT) observations. Aims: The main aim of the observations is to study the nature of the diffuse radio emission in the galaxy cluster ZwCl 2341.1+0000. Methods: We have carried out GMRT 610, 241, and 157 MHz continuum observations of ZwCl 2341.1+0000. The radio observations are combined with X-ray and optical data of the cluster. Results: The GMRT observations show the presence of a double peripheral radio relic in the cluster ZwCl 2341.1+0000. The spectral index is -0.49 \pm 0.18 for the northern relic and -0.76 \pm 0.17 for the southern relic respectively. We have derived values of 0.48-0.93 microGauss for the equipartition magnetic field strength. The relics are probably associated with an outwards traveling merger shock waves.

preprint2009arXiv

The discovery of diffuse steep spectrum sources in Abell 2256

Context: Hierarchical galaxy formation models indicate that during their lifetime galaxy clusters undergo several mergers. Here we report on the discovery of three diffuse radio sources in the periphery of Abell 2256, using the GMRT. Aims: The aim of the observations was to search for diffuse ultra-steep spectrum radio sources within the galaxy cluster Abell 2256. Methods: We have carried out GMRT 325 MHz radio continuum observations of Abell 2256. V, R and I band images of the cluster were taken with the 4.2m WHT. Results: We have discovered three diffuse elongated radio sources located about 1 Mpc from the cluster center. Two are located to the west of the cluster center, and one to the southeast. The sources have a measured physical extent of 170, 140 and 240 kpc, respectively. The two western sources are also visible in deep low-resolution 115-165 MHz WSRT images, although they are blended into a single source. For the combined emission of the blended source we find an extreme spectral index of -2.05\pm 0.14 between 140 and 351 MHz. The extremely steep spectral index suggests these two sources are most likely the result of adiabatic compression of fossil radio plasma due to merger shocks. Conclusions: The discovery of the steep spectrum sources implies the existence of a population of faint diffuse radio sources in (merging) clusters with such steep spectra that they have gone unnoticed in higher frequency (\gtrsim 1 GHz) observations. An exciting possibility therefore is that such sources will determine the general appearance of clusters in low-frequency high resolution radio maps as will be produced by for example LOFAR or LWA.