Westerbork Synthesis Radio Telescope Research Articles

Turbulent circumnuclear disc and cold gas outflow in the newborn radio source 4C 31.04

We present deep kiloparsec- and parsec-scale neutral atomic hydrogen (H I) absorption observations of a very young radio source (≤5000 years), 4C 31.04, using the Westerbork Synthesis Radio Telescope (WSRT) and the Global Very Long Baseline Interferometry (VLBI) array. Using z = 0.0598, derived from molecular gas observations, we detect, at both kpc and pc scales, a broad absorption feature (FWZI = 360 km s−1) centred at the systemic velocity, and narrow absorption (FWZI = 6.6 km s−1) redshifted by 220 km s−1, both previously observed. Additionally, we detect a new blueshifted, broad, shallow absorption wing. At pc scales, the broad absorption at the systemic velocity is detected across the entire radio source while the shallow wing is only seen against part of the eastern lobe. The gas has higher H I column density along the eastern lobe than along the western one. The velocity dispersion of the gas is high (≥40 km s−1) along the entire radio continuum, and is highest (≥60 km s−1) in the region including the outflow and the radio hot spot. While we detect a velocity gradient along the western lobe and parts of the eastern lobe at PA ∼ 5° −10°, most of the gas along the rest of the eastern lobe exhibits no signs of rotation. Earlier optical spectroscopy suggests that the optical AGN is very weak. We therefore conclude that the radio lobes of 4C 31.04 are expanding into a circumnuclear disc, partially disrupting it and making the gas highly turbulent. The distribution of gas is predominantly smooth at the spatial resolution of ∼4 pc studied here. However, clumps of gas are also present, particularly along the eastern lobe. This lobe appears to be strongly interacting with the clouds and driving an outflow ∼35 pc from the radio core, with a mass-outflow rate of 0.3 ≤ Ṁ ≤ 1.4 M⊙ year−1. It is likely that this interaction has caused the eastern lobe to be rebrightened, giving the source an asymmetric morphology. We compare our observations with the predictions of a recent analytical model regarding the survival of atomic gas clouds in radio-jet-driven outflows and find that the existence of a subkpc-scale outflow in this case could imply inefficient mixing of the cold gas with the hot medium and high gas density, leading to very short cooling times. Overall, our study provides further evidence of the strong impact of radio jets on the cold interstellar medium (ISM) in the early stages of their evolution and supports the predictions of numerical simulations regarding jet–ISM interactions and the nature of the circumnuclear gas into which the jets expand.

FAST observations of neutral hydrogen in the interacting galaxies NGC 3395/3396

ABSTRACT We report on high-sensitivity neutral hydrogen observations toward the gas-rich interacting galaxies NGC 3395/3396 with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Compared to previous observations carried out by the Very Large Array (VLA) and the Westerbork Synthesis Radio Telescope (WSRT), a more extended H i envelope around this system has been detected. The total H i gas mass of the NGC 3395/3396 system is estimated to be 7.8 × 109 M$_\odot$. This value is 2.7 times more than that reported based on the VLA interferometric maps. Previous observations found a large H i tail extending to the south-west and a minor tail emerging from the north of this peculiar galaxy pair. Based on the high-sensitivity observations of FAST, an extended H i plume to the north-west and a gas plume to the north-east have been detected for the first time. Neutral hydrogen of the two smaller galaxies IC 2604 and IC 2608 in the south of the system have also been detected. We discuss the origins of these extra gas and possible tidal interactions between these galaxies. The most prominent tidal feature of NGC 3395/3396, the south-west tail combined with the new detected north-west plume, behaves like a large ring. We suggest the ring might be formed by the previous fly-by interaction between NGC 3395 and 3396, which happened 500 Myr ago. Our study shows that high-sensitivity H i observations are important in revealing low column density gas, which is crucial to a deeper understanding of this interacting system.

A Candidate Period of 4.605 Days for FRB 20121102A and One Possible Implication of Its Origin

A firm establishment of the presence or the lack of periodicity in repeating fast radio bursts is crucial for determining their origins. Here, we compile 1145 radio bursts of FRB 20121102A with fluence larger than 0.15 Jy ms from observations using the Five-hundred-meter Aperture Spherical radio Telescope, the Arecibo Observatory, the Green Bank Telescope, the Effelsberg Telescope, the MeerKAT Telescope, the Lovell Telescope, the Deep Space Network 70 m radio telescopes, the Very Large Array, and the Westerbork Synthesis Radio Telescope, spanning the time interval of MJD 57175−58776. A quasi-period of 157.1−4.8+5.2 days and a candidate quasi-period of 4.605−0.010+0.003 days are found through the phase-folding probability binomial analysis. The former is consistent with previous findings and the latter is new. The 4.605 day periodicity is more obvious in high-energy bursts with fluence larger than 1 × 1038 erg. The presence of these (candidate) quasi-periods, together with the corresponding width of burst accumulation in the phase space, are consistent with the bursts originating from a binary degenerate star system with a close-by planet around the primary neutron star.

The temperature of the neutral interstellar medium in the Galaxy

ABSTRACT Atomic Hydrogen-21 cm transition (H i) is an excellent tracer to study and understand the properties of the atomic gas in the Galaxy. Using the Westerbork Synthesis Radio Telescope, we observed 12 quasar sightlines to detect galactic H i in absorption. We achieve an optical depth rms of ∼1−2 × 10−3, essential to detect the warm neutral medium (WNM). We detect H i absorption in all our sightlines except along 1006+349, for which we set a strict upper limit on the spin temperature as 〈Ts〉 > 570 K. We find around 50 per cent of our sightlines have 〈Ts〉 > 500 K, indicating a WNM dominance. Further, we calculate an upper limit of the cold neutral medium (CNM) fraction along our sightlines and find a median CNM fraction of ∼0.12. With our observations, we reconfirm the existence of a threshold column density of ∼2 × 1020$\rm atoms \, cm^{-2}$ to form CNM in the interstellar medium (ISM). Using a two-temperature model of the H i disc, we explore the distribution of spin temperature in the Galactic ISM. We find that a simple fixed axisymmetric two-temperature model could not produce either the observed column density or the integral optical depth. This indicates the existence of a more complex distribution of spin temperatures in the Galaxy.

Dark matter measurements combining stellar and H i kinematics: 30 per cent 1σ outliers with low dark matter content at 5Re

ABSTRACT We construct the Schwarzschild dynamical models for 11 early-type galaxies with the SAURON and Mitchell stellar IFUs out to 2–4Re, and construct dynamical models with combined stellar and H i kinematics for a subsample of four galaxies with H i velocity fields out to 10Re obtained from the Westerbork Synthesis Radio Telescope, thus robustly obtaining the dark matter content out to large radii for these galaxies. Adopting a generalized-NFW dark matter profile, we measure an NFW-like density cusp in the dark matter inner slopes for all sample galaxies, with a mean value of 1.00 ± 0.04 (rms scatter 0.15). The mean dark matter fraction for the sample is 0.2 within 1Re, and increases to 0.4 at 2Re, and 0.6 at 5Re. The dark matter fractions within 1Re of these galaxies are systematically lower than the predictions of both the TNG-100 and EAGLE simulations. For the dark matter fractions within 2Re and 5Re, 40 and 70 per cent galaxies are 1σ consistent with either the TNG-100 or the EAGLE predictions, while the remaining 60 and 30 per cent galaxies lie below the 1σ region. Combined with 36 galaxies with dark matter fractions measured out to 5Re in the literature, about 10 per cent of these 47 galaxies lie below the 3σ region of the TNG-100 or EAGLE predictions.

Disentangling Magnetic Fields in NGC 6946 with Wide-band Polarimetry

We present λ13 cm polarization observations of the nearby spiral galaxy NGC 6946 with the Westerbork Synthesis Radio Telescope (WSRT) to examine the nearside halo magnetic fields. Despite λ13 cm exhibiting similar two-dimensional morphology as observed at longer (λ18–22 cm) or shorter (λ3 and λ6 cm) wavelengths, more complete frequency coverage will be required to explain the gap in polarization in the southwest quadrant of the galaxy. We fit models of the turbulent and coherent line-of-sight magnetic fields to the fractional degree of linearly polarized emission at λ3, λ6, λ13, λ18, and λ22 cm from observations taken with the WSRT, Karl G. Jansky Very Large Array, and Effelsberg telescopes. The results favor a multilayer turbulent magneto-ionized medium consistent with current observations of edge-on galaxies. We constrain the physical properties of the synchrotron-emitting thin and thick disks (scale heights of 300 pc and 1.4 kpc, respectively) along with the thermal thick disk and halo (scale heights of 1 and 5 kpc, respectively). Our preferred model indicates a clumpy and highly turbulent medium within 1 kpc of the midplane, and a diffuse extraplanar layer with a substantially lower degree of Faraday depolarization. In the halo, we estimate a regular magnetic field strength of 0.4–2.2 μG and that turbulence and a total magnetic field strength of ∼6 μG result in a Faraday dispersion of σ RM = 4–48 rad m−2. This work is an example of how the advanced capabilities of modern radio telescopes are opening a new frontier for the study of cosmic magnetism.

Long baseline bistatic radar imaging of tumbling space objects for enhancing space domain awareness

AbstractLong baseline bistatic radar systems herald enhanced sensitivity and metric accuracy for space objects in geosynchronous orbits and beyond. Radio telescopes are ideal participants in such a system; in particular, they often feature large apertures with low‐noise temperatures and have stable, synchronised clocks. Pairing radio telescopes with high‐power radars creates new methodologies for Space Domain Awareness. This paper describes long baseline bistatic measurements using the Millstone Hill Radar in the USA, the Tracking and Imaging Radar in Germany, multiple receivers of the enhanced multi‐element remotely linked interferometer network array in the United Kingdom, and the Westerbork Synthesis Radio Telescope in the Netherlands. The authors, a Research Task Group formed by the NATO Science and Technology Organisation Sensors and Electronic Technology Panel (SET‐293), performed novel bistatic and monostatic radar imaging experiments with real on‐orbit tumbling rocket bodies. These experiments on tumbling objects at near‐geosynchronous orbits highlight successful demonstrations of advanced bistatic Doppler characterisation across diverse imaging geometries. Specialised Doppler processing on tumbling targets, such as the Doppler superpulse algorithm, enables high‐fidelity rotation period estimation and determination of minimum target size.

Detection of diffuse H i emission in the circumgalactic medium of NGC 891 and NGC 4565 - II

ABSTRACT We probe the neutral circumgalactic medium (CGM) along the major axes of NGC 891 and NGC 4565 in 21-cm emission out to ≳100 kpc using the Green Bank Telescope (GBT), extending our previous minor axes observations. We achieve an unprecedented 5σ sensitivity of 6.1 × 1016 cm−2 per 20 km s−1 velocity channel. We detect H i with diverse spectral shapes, velocity widths, and column densities. We compare our detections to the interferometric maps from the Westerbork Synthesis Radio Telescope (WSRT) obtained as part of the HALOGAS survey. At small impact parameters, $\gt 31\!-\!43~{{\ \rm per\ cent}}$ of the emission detected by the GBT cannot be explained by emission seen in the WSRT maps, and it increases to $\gt 64\!-\!73~{{\ \rm per\ cent}}$ at large impact parameters. This implies the presence of diffuse circumgalactic H i. The mass ratio between H i in the CGM and H i in the disc is an order of magnitude larger than previous estimates based on shallow GBT mapping. The diffuse H i along the major axes pointings is corotating with the H i disc. The velocity along the minor axes pointings is consistent with an inflow and/or fountain in NGC 891 and an inflow/outflow in NGC 4565. Including the circumgalactic H i, the depletion time and the accretion rate of NGC 4565 are sufficient to sustain its star formation. In NGC 891, most of the required accreting material is still missing.

Radio multifrequency observations of Abell 781 with the WSRT

ABSTRACT The Main galaxy cluster in the Abell 781 system is undergoing a significant merger and accretion process with peripheral emission to the north and south-eastern flanks of the merging structure. Here we present a full polarimetric study of this field, using radio interferometric data taken at 21 and 92 cm with the Westerbork Synthesis Radio Telescope (WSRT), to a sensitivity better than any 21 cm (L band) observation to date. We detect evidence of extended low-level emission of 1.9 mJy associated with the Main cluster at 21 cm, although this detection necessitates further follow-up by modern instruments due to the limited resolution of the WSRT. Our polarimetric study indicates that, most likely, the peripheral emission associated with this cluster is not a radio relic.

Apertif 1.4 GHz continuum observations of the Boötes field and their combined view with LOFAR

We present a new image of a 26.5 square degrees region in the Boötes constellation obtained at 1.4 GHz using the Aperture Tile in Focus (Apertif) system on the Westerbork Synthesis Radio Telescope. We use a newly developed processing pipeline that includes direction-dependent self-calibration, which provides a significant improvement in the quality of the images compared to those released as part of the Apertif first data release. For the Boötes region, we mosaicked 187 Apertif images and extracted a source catalog. The mosaic image has an angular resolution of 27 × 11.5″ and a median background noise of 40 μJy beam−1. The catalog has 8994 sources and is complete down to the 0.3 mJy level. We combined the Apertif image with LOFAR images of the Boötes field at 54 and 150 MHz to study the spectral properties of the sources. We find a spectral flattening toward sources with a low flux density. Using the spectral index limits from Apertif nondetections, we derive that up to 9% of the sources have ultrasteep spectra with a slope below −1.2. A steepening of the spectral index with increasing redshift is also seen in the data, which shows a different dependence for the low-and high-frequency spectral index. The explanation probably is that a population of sources has concave radio spectra with a turnover frequency of about the LOFAR band. Additionally, we discuss cases of individual extended sources with an interesting resolved spectral structure. With the improved pipeline, we aim to continue to process data from the Apertif wide-area surveys and release the improved 1.4-GHz images of several well-known fields.

Radio multifrequency observations of the galaxy cluster pair Abell 399 − 401 with WSRT

ABSTRACT Galaxy clusters are assembled via merging of smaller structures, in a process that generates shocks and turbulence in the intra cluster medium and produces radio diffuse emission in the form of halos and relics. The cluster pair A 399–A 401 represents a special case: both clusters host a radio halo. Recent Low Frequency Array (LOFAR) observations at 140 MHz revealed the presence of a radio bridge connecting the two clusters along with two relic candidates. These relics include one South of A 399 and the other in between the two clusters, in proximity of a shock front detected in X-ray observations. In this paper we present observations of the A 399–A 401 cluster pair at 1.7, 1.4, 1.2 GHz, and 346 MHz from the Westerbork Synthesis Radio Telescope (WSRT). We detect the radio halo in the A 399 cluster at 346 MHz, extending up to ∼650 kpc and with a 125 ± 6 mJy flux density. Its spectral index between 140 MHz and 346 MHz is α = 1.75 ± 0.14. The two candidate relics are also seen at 346 MHz and we determine their spectral indices to be α = 1.10 ± 0.14 and α = 1.46 ± 0.14. The low surface brightness bridge connecting the two clusters is below the noise level at 346 MHz therefore we constrain the bridge average spectral index to be steep, i.e. α > 1.5 at 95 per cent confidence level. This result favours the scenario where dynamically-induced turbulence is a viable mechanism to reaccelerate a population of mildly relativistic particles and amplify magnetic fields on scales of a few Mpcs.

Catching Tidal Dwarf Galaxies at a Later Evolutionary Stage with ALFALFA

We present deep optical imaging and photometry of four objects classified as “Almost-Dark” galaxies in the Arecibo Legacy Fast Arecibo L-band Feed Array (ALFALFA) survey because of their gas-rich nature and extremely faint or missing optical emission in existing catalogs. They have H i masses of 107–109 M ⊙ and distances of ∼9–100 Mpc. Observations with the WIYN 3.5 m telescope and One Degree Imager reveal faint stellar components with central surface brightnesses of ∼24–25 in the g band. We also present the results of H i synthesis observations with the Westerbork Synthesis Radio Telescope. These Almost-Dark galaxies have been identified as possible tidal dwarf galaxies (TDGs) based on their proximity to one or more massive galaxies. We demonstrate that AGC 229398 and AGC 333576 likely have the low dark matter content and large effective radii representative of TDGs. They are located much farther from their progenitors than previously studied TDGs, suggesting they are older and more evolved. AGC 219369 is likely dark matter dominated, while AGC 123216 has a dark matter content that is unusually high for a TDG, but low for a normal dwarf galaxy. We consider possible mechanisms for the formation of the TDG candidates such as a traditional major merger scenario and gas ejection from a high-velocity flyby. Blind H i surveys like ALFALFA enable the detection of gas-rich, optically faint TDGs that can be overlooked in other surveys, thereby providing a more complete census of the low-mass galaxy population and an opportunity to study TDGs at a more advanced stage of their life cycle.

The Apertif Radio Transient System (ARTS): Design, commissioning, data release, and detection of the first five fast radio bursts

Fast radio bursts (FRBs) must be powered by uniquely energetic emission mechanisms. This requirement has eliminated a number of possible source types, but several remain. Identifying the physical nature of FRB emitters arguably requires good localisation of more detections, as well as broad-band studies enabled by real-time alerting. In this paper, we present the Apertif Radio Transient System (ARTS), a supercomputing radio-telescope instrument that performs real-time FRB detection and localisation on the Westerbork Synthesis Radio Telescope (WSRT) interferometer. It reaches coherent-addition sensitivity over the entire field of the view of the primary-dish beam. After commissioning results verified that the system performed as planned, we initiated the Apertif FRB survey (ALERT). Over the first 5 weeks we observed at design sensitivity in 2019, we detected five new FRBs, and interferometrically localised each of them to 0.4–10 sq. arcmin. All detections are broad band, very narrow, of the order of 1 ms in duration, and unscattered. Dispersion measures are generally high. Only through the very high time and frequency resolution of ARTS are these hard-to-find FRBs detected, producing an unbiased view of the intrinsic population properties. Most localisation regions are small enough to rule out the presence of associated persistent radio sources. Three FRBs cut through the halos of M31 and M33. We demonstrate that Apertif can localise one-off FRBs with an accuracy that maps magneto-ionic material along well-defined lines of sight. The rate of one every ~7 days ensures a considerable number of new sources are detected for such a study. The combination of the detection rate and localisation accuracy exemplified by the first five ARTS FRBs thus marks a new phase in which a growing number of bursts can be used to probe our Universe.

Multi-epoch variability of AT 2000ch (SN 2000ch) in NGC 3432

Context.AT 2000ch is a highly variable massive star and supernova imposter in NGC 3432 first detected in 2000. It is similar and often compared to SN 2009ip, and it is therefore expected to undergo a core-collapse supernova (SN) – a SN imposter of similar brightness – in the near future.Aims.We characterize the long-term variability of AT 2000ch in the radio and optical regimes with archival data reaching back to the year 1984. We use these newly reduced observations in addition to observations in the literature to restrict the mass-loss rates of AT 2000ch at multiple epochs based on different approaches, and to infer the general properties of its circumstellar nebula with respect to the detected radio brightness.Methods.We extend the known optical light curve of AT 2000ch up to the beginning of 2022 by performing point spread function (PSF) photometry on archival data from the Palomar Transient Factory and the Zwicky Transient Facility. We reduced archival radio continuum observations obtained with the Very Large Array using standard calibration and imaging methods and complemented these with pre-reduced CHANG-ESobservations as well as observations obtained with the Westerbork Synthesis Radio Telescope and LOw Frequency ARray. For the analysis of AT 2000ch, we consider the optical light curve and color evolution, its radio continuum brightness at different frequencies and times, and the corresponding spectral indices. We estimated mass-loss rates and optical depths based on radio continuum brightnesses and Hαfluxes.Results.We report two newly detected outbursts of AT 2000ch similar to those found in the 2000 s and 13 re-brightening events, of which at least four are not conclusively detected because of insufficient sampling of the light curve. The dates of all outbursts and significant, well-sampled re-brightening events are consistent with a period of ∼201 ± 12 days over a total time-span of two decades. Such a behavior has never been found for any SN imposter, especially not for candidate SN 2009ip analogs. During 2010 to 2012 and 2014 to 2018, we only have a few detections, which is insufficient to come to any conclusion as to a possible less eruptive phase of the transient. We find steady dimming after the most recent re-brightening events and possible evidence of porosity in the circumstellar envelope, suggesting AT 2000ch may currently be in transition to a state of relative calm. We identified a second, unrelated source at a projected distance of ∼23 pc (∼0.5″) that has contaminated the optical measurements of AT 2000ch at its minimum luminosity over the last two decades probably on a 5%−10% level, but this does not affect our overall findings and is negligible during re-brightening. We are able to restrict the mass-loss rate of AT 2000ch to range between several 10−6M⊙yr−1and several 10−5M⊙yr−1. The fresh ejecta appear to be optically thick to radio continuum emission at least within the first ∼25 days after significant re-brightening.Conclusions.We suggest that other SN imposter and probably also candidate SN 2009ip-analogs at comparable distances emit radio continuum fluxes on the order of a few to several tens of microJanskys at 6 GHz. Deep and frequent continuum surveys in the radio range are needed to study these kinds of objects in a broader context.

A comparative study of source-finding techniques in H I emission line cubes using SoFiA, MTObjects, and supervised deep learning

Context. The 21 cm spectral line emission of atomic neutral hydrogen (H I) is one of the primary wavelengths observed in radio astronomy. However, the signal is intrinsically faint and the H I content of galaxies depends on the cosmic environment, requiring large survey volumes and survey depth to investigate the H I Universe. As the amount of data coming from these surveys continues to increase with technological improvements, so does the need for automatic techniques for identifying and characterising H I sources while considering the tradeoff between completeness and purity. Aims. This study aimed to find the optimal pipeline for finding and masking the most sources with the best mask quality and the fewest artefacts in 3D neutral hydrogen cubes. Various existing methods were explored, including the traditional statistical approaches and machine learning techniques, in an attempt to create a pipeline to optimally identify and mask the sources in 3D neutral hydrogen (H I) 21 cm spectral line data cubes. Methods. Two traditional source-finding methods were tested first: the well-established H I source-finding software SoFiA and one of the most recent, best performing optical source-finding pieces of software, MTObjects. A new supervised deep learning approach was also tested, in which a 3D convolutional neural network architecture, known as V-Net, which was originally designed for medical imaging, was used. These three source-finding methods were further improved by adding a classical machine learning classifier as a post-processing step to remove false positive detections. The pipelines were tested on H I data cubes from the Westerbork Synthesis Radio Telescope with additional inserted mock galaxies. Results. Following what has been learned from work in other fields, such as medical imaging, it was expected that the best pipeline would involve the V-Net network combined with a random forest classifier. This, however, was not the case: SoFiA combined with a random forest classifier provided the best results, with the V-Net–random forest combination a close second. We suspect this is due to the fact that there are many more mock sources in the training set than real sources. There is, therefore, room to improve the quality of the V-Net network with better-labelled data such that it can potentially outperform SoFiA.

FEASTS: IGM Cooling Triggered by Tidal Interactions through the Diffuse H i Phase around NGC 4631

We use the single-dish radio telescope Five-hundred meter Aperture Spherical Telescope (FAST) to map the H i in the tidally interacting NGC 4631 group with a resolution of 3.′24 (7 kpc), reaching a 5σ column density limit of 1017.9 cm−2 assuming a line width of 20 km s−1. Taking the existing interferometric H i image from the Hydrogen Accretion in LOcal GAlaxieS project of Westerbork Synthesis Radio Telescope as a reference, we are able to identify and characterize a significant excess of large-scale, low-density, and diffuse H i in the group. This diffuse H i extends more than 120 kpc across, and accounts for more than one-fourth of the total H i detected by FAST in and around the galaxy NGC 4631. In the region of the tidal tails, the diffuse H i has a typical column density above 1019.5 cm−2, and is highly turbulent with a velocity dispersion of around 50 km s−1. It increases in column density with the dense H i, and tends to be associated with the kinematically hotter part of the dense H i. Through simple modeling, we find that the majority of the diffuse H i in the tail region is likely to induce cooling out of the hot intragalactic medium (IGM) instead of evaporating or being radiatively ionized. Given these relations of gas in different phases, the diffuse H i may represent a condensing phase of the IGM. Ongoing and past active tidal interactions may have produced the wide-spreading H i distribution, and triggered the gas accretion to NGC 4631 through the phase of the diffuse H i.

An interference detection strategy for Apertif based on AOFlagger 3

Context. Apertif is a multi-beam receiver system for the Westerbork Synthesis Radio Telescope that operates at 1.1–1.5 GHz, which overlaps with various radio services, resulting in contamination of astronomical signals with radio-frequency interference (RFI). Aims. We analyse approaches to mitigate Apertif interference and design an automated detection procedure for its imaging mode. Using this approach, we present long-term RFI detection results of over 300 Apertif observations. Methods. Our approach is based on the AOFlagger detection approach. We introduce several new features, including ways to deal with ranges of invalid data (e.g. caused by shadowing) in both the SumThreshold and scale-invariant rank operator steps; pre-calibration bandpass calibration; auto-correlation flagging; and HI flagging avoidance. These methods have been implemented in a new framework that uses the Lua language for scripting, which is new in AOFlagger version 3. Results. Our approach removes RFI fully automatically, and it is robust and effective enough for further calibration and (continuum) imaging of these data. The analysis of 304 observations shows an average of 11.1% of lost data due to RFI with a large spread. We observe 14.6% RFI in auto-correlations. Computationally, AOFlagger achieves a throughput of 370 MB/s on a single computing node. Compared to published machine learning results, the method is one to two orders of magnitude faster.

HALOGAS: Strong constraints on the neutral gas reservoir and accretion rate in nearby spiral galaxies

Context. Galaxies in the local Universe are thought to require ongoing replenishment of their gas reservoir in order to maintain the observed star formation rates. Cosmological simulations predict that this accretion can occur in both a dynamically hot and cold mode, depending on the redshift, halo mass, and the details of the included feedback processes. However, until now, observational evidence of the accretion required to match the observed star formation histories is lacking. Aims. Within the framework of the Hydrogen Accretion in LOcal GalaxieS (HALOGAS) survey, this paper attempts to determine whether galaxies in the local Universe possess a significant reservoir of cold neutral gas and the possible accretion rates these reservoirs could provide. Additionally, with this moderately sized sample, we can start to investigate whether the observed accretion is connected to intrinsic properties of the galaxies such as Hubble type, star formation rate, or environment. Methods. We searched the vicinity of 22 nearby galaxies in a systematic and automated manner for isolated H I clouds or distinct streams that are not yet connected to the galaxy disk. The HALOGAS observations were carried out with the Westerbork Synthesis Radio Telescope and represent one of the most sensitive and detailed H I surveys to date. These observations typically reach column density sensitivities of ∼1019 cm−2 over a 20 km s−1 line width. Results. We find 14 secure H I cloud candidates without an observed optical counterpart in the entire HALOGAS sample. These cloud candidates appear to be analogs to the most massive clouds detected in the extensive cloud distributions around the Milky Way and M 31. However, their numbers seem significantly reduced on average compared to the numbers in these galaxies. Within the framework of cold accretion, we constrain upper limits for H I accretion in the local Universe. The average H I mass currently observed in a state suggestive of accretion onto the galaxies amounts to a rate of 0.05 M⊙ yr−1 with a stringent upper limit of 0.22 M⊙ yr−1, confirming previous estimates. This is much lower than the average star formation rate in this sample. Our best estimate, based on the Green Bank Telescope detection limits of several galaxies in the sample, suggests that another 0.04 M⊙ yr−1 of neutral hydrogen at most could be accreted from clouds and streams that remain undetected. Conclusions. These results show that in nearby galaxies, neutral hydrogen is not being accreted at the same rate as stars are currently being formed. Our study cannot exclude that other forms of gas accretion are at work, such as those provided by direct infall of ionized intergalactic gas or the condensation of coronal gas, triggered by galactic fountain activities. However, these observations do not reveal extended neutral gas reservoirs around most nearby spiral galaxies either.

Characterising the Apertif primary beam response

Context. Phased array feeds (PAFs) are multi-element receivers in the focal plane of a telescope that make it possible to simultaneously form multiple beams on the sky by combining the complex gains of the individual antenna elements. Recently, the Westerbork Synthesis Radio Telescope (WSRT) was upgraded with PAF receivers to carry out several observing programs, including two imaging surveys and a time-domain survey. The Apertif imaging surveys use a configuration of 40 partially overlapping compound beams (CBs) simultaneously formed on the sky and arranged in an approximately rectangular shape. Aims. This work is aimed at characterising the response of the 40 Apertif CBs to create frequency-resolved I, XX, and YY polarization empirical beam shapes. The measured CB maps can be used for the image deconvolution, primary beam correction, and mosaicking processes of Apertif imaging data. Methods. We used drift scan measurements to measure the response of each of the 40 Apertif CBs. We derived beam maps for all individual beams in I, XX, and YY polarisation in 10 or 18 frequency bins over the same bandwidth as the Apertif imaging surveys. We sampled the main lobe of the beams and the side lobes up to a radius of 0.6 degrees from the beam centres. In addition, we derived beam maps for each individual WSRT dish. Results. We present the frequency and time dependence of the beam shapes and sizes. We compared the compound beam shapes derived with the drift scan method to beam shapes derived with an independent method using a Gaussian Process Regression comparison between the Apertif continuum images and the NRAO VLA Sky Survey (NVSS) catalogue. We find a good agreement between the beam shapes derived with the two independent methods.

First release of Apertif imaging survey data

Context. Apertif is a phased-array feed system for the Westerbork Synthesis Radio Telescope, providing forty instantaneous beams over 300 MHz of bandwidth. A dedicated survey program utilizing this upgrade started on 1 July 2019, with the last observations taken on 28 February 2022. The imaging survey component provides radio continuum, polarization, and spectral line data. Aims. Public release of data is critical for maximizing the legacy of a survey. Toward that end, we describe the release of data products from the first year of survey operations, through 30 June 2020. In particular, we focus on defining quality control metrics for the processed data products. Methods. The Apertif imaging pipeline, Apercal, automatically produces non-primary beam corrected continuum images, polarization images and cubes, and uncleaned spectral line and dirty beam cubes for each beam of an Apertif imaging observation. For this release, processed data products are considered on a beam-by-beam basis within an observation. We validate the continuum images by using metrics that identify deviations from Gaussian noise in the residual images. If the continuum image passes validation, we release all processed data products for a given beam. We apply further validation to the polarization and line data products and provide flags indicating the quality of those data products. Results. We release all raw observational data from the first year of survey observations, for a total of 221 observations of 160 independent target fields, covering approximately one thousand square degrees of sky. Images and cubes are released on a per beam basis, and 3374 beams (of 7640 considered) are released. The median noise in the continuum images is 41.4 uJy beam−1, with a slightly lower median noise of 36.9 uJy beam−1 in the Stokes V polarization image. The median angular resolution is 11.6″/sin δ. The median noise for all line cubes, with a spectral resolution of 36.6 kHz, is 1.6 mJy beam−1, corresponding to a 3-σ H i column density sensitivity of 1.8 × 1020 atoms cm−2 over 20 km s−1 (for a median angular resolution of 24″ × 15″). Line cubes at lower frequency have slightly higher noise values, consistent with the global RFI environment and overall Apertif system performance. We also provide primary beam images for each individual Apertif compound beam. The data are made accessible using a Virtual Observatory interface and can be queried using a variety of standard tools.