Search

Galaxy formation theory identifies superwinds as a key regulator of star formation rates, galaxy growth, and chemical enrichment. Thermal and radiation pressure are known to drive galactic-scale winds in dusty starbursting galaxies (e.g. M82), but modern numerical simulations have recently highlighted that cosmic-ray (CR) driven winds may be especially important in normal galaxies with modest star formation rate surface densities. However, CR-driven winds have yet to be conclusively observed -- leaving significant uncertainty in their detailed microphysics. We present MeerKAT radio continuum and HI spectral-line observations of one such normal galaxy, NGC 1532; a nearby ($D\sim15\,\mathrm{Mpc}$) and edge-on ($i \gtrsim 80^{\circ}$) spiral galaxy tidally interacting with its smaller elliptical companion, NGC 1531. We find magnetized, highly-ordered radio continuum loops extending $\sim10$ kpc above and below the disk; visibly connecting discrete star-forming regions in the disk with the nucleus. The deep MeerKAT HI observations place an upper limit on the column density of neutral gas coincident with the outflow to $N_\mathrm{HI} \lesssim 3 \times 10^{19}\,\mathrm{cm}^{-2}$. Unlike previously observed outflows -- for which ejected gas and dust can be traced across multiple wavelengths -- the loops in NGC 1532 show no detectable signs of dust or gas coincident with the radio emission far from the disk. We explore multiple possible mechanisms for driving this magnetic wind and favor an explanation where cosmic-ray pressure plays a significant role in launching these outflows., Comment: 9 pages, 5 figures, submitted to ApJL

We present a detailed study of the large-scale shock front in Stephan's Quintet, a byproduct of past and ongoing interactions. Using integral-field spectroscopy from the new William Herschel Telescope Enhanced Area Velocity Explorer (WEAVE), recent 144 MHz observations from the LOFAR Two-metre Sky Survey (LoTSS), and archival data from the Very Large Array and James Webb Space Telescope (JWST), we obtain new measurements of key shock properties and determine its impact on the system. Harnessing the WEAVE large integral field unit's (LIFU) field of view (90 $\times$ 78 arcsec$^{2}$), spectral resolution ($R\sim2500$) and continuous wavelength coverage across the optical band, we perform robust emission line modeling and dynamically locate the shock within the multi-phase intergalactic medium (IGM) with higher precision than previously possible. The shocking of the cold gas phase is hypersonic, and comparisons with shock models show that it can readily account for the observed emission line ratios. In contrast, we demonstrate that the shock is relatively weak in the hot plasma visible in X-rays (with Mach number of $\mathcal{M} \sim 2 - 4$), making it inefficient at producing the relativistic particles needed to explain the observed synchrotron emission. Instead, we propose that it has led to an adiabatic compression of the medium, which has increased the radio luminosity ten-fold. Comparison of the Balmer line-derived extinction map with the molecular gas and hot dust observed with JWST suggests that pre-existing dust may have survived the collision, allowing the condensation of H$_{2}$ - a key channel for dissipating the shock energy., Comment: 23 pages, 15 figures, accepted for publication in MNRAS

Context. Understanding the dwarf galaxy population in low density environments is crucial for testing the LCDM cosmological model. The increase in diversity towards low mass galaxies is seen as an increase in the scatter of scaling relations such as the stellar mass-size and the baryonic Tully-Fisher relation (BTFR), and is also demonstrated by recent in-depth studies of an extreme subclass of dwarf galaxies of low surface brightness, but large physical sizes, called ultra-diffuse galaxies (UDGs). Aims. We select galaxies from the Apertif HI survey, and apply a constraint on their i-band absolute magnitude to exclude high mass systems. The sample consists of 24 galaxies, and span HI mass ranges of 8.6 < log ($M_{HI}/M_{Sun}$) < 9.7 and stellar mass range of 8.0 < log ($M_*/M_{Sun}$) < 9.7 (with only three galaxies having log ($M_*/M_{Sun}$) > 9). Methods. We determine the geometrical parameters of the HI and stellar discs, build kinematic models from the HI data using 3DBarolo, and extract surface brightness profiles in g-, r- and i-band from the Pan-STARRS 1 photometric survey. Results. We find that, at fixed stellar mass, our HI selected dwarfs have larger optical effective radii than isolated, optically-selected dwarfs from the literature. We find misalignments between the optical and HI morphologies for some of our sample. For most of our galaxies, we use the HI morphology to determine their kinematics, and we stress that deep optical observations are needed to trace the underlying stellar discs. Standard dwarfs in our sample follow the same BTFR of high-mass galaxies, whereas UDGs are slightly offset towards lower rotational velocities, in qualitative agreement with results from previous studies. Finally, our sample features a fraction (25%) of dwarf galaxies in pairs that is significantly larger with respect to previous estimates based on optical spectroscopic data., Comment: 22 pages, 11 figures, 5 tables. Accepted for publication in A&A

We present the Pilot Survey Phase 2 data release for the Wide-field ASKAP L-band Legacy All-sky Blind surveY (WALLABY), carried-out using the Australian SKA Pathfinder (ASKAP). We present 1760 HI detections (with a default spatial resolution of 30") from three pilot fields including the NGC 5044 and NGC 4808 groups as well as the Vela field, covering a total of ~180 deg$^2$ of the sky and spanning a redshift up to $z \simeq 0.09$. This release also includes kinematic models for over 126 spatially resolved galaxies. The observed median rms noise in the image cubes is 1.7 mJy per 30" beam and 18.5 kHz channel. This corresponds to a 5$\sigma$ HI column density sensitivity of $\sim 9.1\times10^{19}(1 + z)^4$ cm$^{-2}$ per 30" beam and $\sim 20$ km/s channel, and a 5$\sigma$ HI mass sensitivity of $\sim 5.5\times10^8 (D/100$ Mpc)$^{2}$ M$_{\odot}$ for point sources. Furthermore, we also present for the first time 12" high-resolution images ("cut-outs") and catalogues for a sub-sample of 80 sources from the Pilot Survey Phase 2 fields. While we are able to recover sources with lower signal-to-noise ratio compared to sources in the Public Data Release 1, we do note that some data quality issues still persist, notably, flux discrepancies that are linked to the impact of side lobes associated with the dirty beams due to inadequate deconvolution. However, in spite of these limitations, the WALLABY Pilot Survey Phase 2 has already produced roughly a third of the number of HIPASS sources, making this the largest spatially resolved HI sample from a single survey to date., Comment: 25 pages, 16 figures, 6 tables, accepted for publication in Publications of the Astronomical Society of Australia (PASA)

Reinforcement learning (RL) has seen significant research and application results but often requires large amounts of training data. This paper proposes two data-efficient off-policy RL methods that use parametrized Q-learning. In these methods, the Q-function is chosen to be linear in the parameters and quadratic in selected basis functions in the state and control deviations from a base policy. A cost penalizing the $\ell_1$-norm of Bellman errors is minimized. We propose two methods: Linear Matrix Inequality Q-Learning (LMI-QL) and its iterative variant (LMI-QLi), which solve the resulting episodic optimization problem through convex optimization. LMI-QL relies on a convex relaxation that yields a semidefinite programming (SDP) problem with linear matrix inequalities (LMIs). LMI-QLi entails solving sequential iterations of an SDP problem. Both methods combine convex optimization with direct Q-function learning, significantly improving learning speed. A numerical case study demonstrates their advantages over existing parametrized Q-learning methods., Comment: Accepted for Presentation at 63rd IEEE Conference on Decision and Control, CDC 2024, Milan, Italy, 2024

We have learnt in the last decades that the majority of galaxies belong to high density regions interconnected in a sponge-like fashion. This large-scale structure is characterised by clusters, filaments, walls, where most galaxies concentrate, but also under-dense regions, called voids. The void regions and the galaxies within represent an ideal place for the study of galaxy formation and evolution as they are largely unaffected by the complex physical processes that transform galaxies in high-density environments. These void galaxies can hold the key as well to answer current challenges to the $\Lambda$CDM paradigm. The Calar Alto Void Integral-field Treasury surveY (CAVITY) is a Legacy project approved by the Calar Alto Observatory to obtain spatially resolved spectroscopic information of $\sim300$ void galaxies in the Local Universe (0.005 < z < 0.050) covering from -17.0 to -21.5 in $\rm r$ band absolute magnitude. It officially started in January 2021 and has been awarded 110 useful dark observing nights at the 3.5 m telescope using the PMAS spectrograph. Complementary follow-up projects including deep optical imaging, integrated, as well as resolved CO data, and integrated HI spectra, have joint the PMAS observations and naturally complete the scientific aim of characterising galaxies in cosmic voids. The extension data has been denominated CAVITY+. The data will be available to the whole community in different data releases, the first of which is planned for July 2024, and it will provide the community with PMAS data cubes for around 100 void galaxies through a user friendly, and well documented, database platform. We present here the survey, sample selection, data reduction, quality control schemes, science goals, and some examples of the scientific power of the CAVITY and CAVITY+ data., Comment: 19 pages, 11 figures. Accepted for publication in A&A - Replacement after A&A minor language edition

We present a new image of a 26.5 square degree region in the Bo\"otes 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 which includes direction-dependent self-calibration which provides a significant improvement of the quality of the images compared to those released as part of the Apertif first data release. For the Bo\"otes region, we mosaic 187 Apertif images and extract a source catalog. The mosaic image has an angular resolution of 27${\times}$11.5 arcseconds and a median background noise of 40 ${\mu}$Jy/beam. The catalog has 8994 sources and is complete down to the 0.3 mJy level. We combine the Apertif image with LOFAR images of the Bo\"otes field at 54 and 150 MHz to study spectral properties of the sources. We find a spectral flattening towards low flux density sources. Using the spectral index limits from Apertif non-detections we derive that up to 9 percent of the sources have ultra-steep spectra with a slope steeper than -1.2. Steepening of the spectral index with increasing redshift is also seen in the data showing a different dependency for the low-frequency spectral index and the high frequency one. This can be explained by a population of sources having concave radio spectra with a turnover frequency around 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 processing data from the Apertif wide-area surveys and release the improved 1.4 GHz images of several famous fields., Comment: 13 pages, 9 figures; to be published in A&A

We investigate the diversity in the sizes and average surface densities of the neutral atomic hydrogen (HI) gas discs in ~280 nearby galaxies detected by the Widefield ASKAP L-band Legacy All-sky Blind Survey (WALLABY). We combine the uniformly observed, interferometric HI data from pilot observations of the Hydra cluster and NGC 4636 group fields with photometry measured from ultraviolet, optical and near-infrared imaging surveys to investigate the interplay between stellar structure, star formation and HI structural parameters. We quantify the HI structure by the size of the HI relative to the optical disc and the average HI surface density measured using effective and isodensity radii. For galaxies resolved by >1.3 beams, we find that galaxies with higher stellar masses and stellar surface densities tend to have less extended HI discs and lower HI surface densities: the isodensity HI structural parameters show a weak negative dependence on stellar mass and stellar mass surface density. These trends strengthen when we limit our sample to galaxies resolved by >2 beams. We find that galaxies with higher HI surface densities and more extended HI discs tend to be more star forming: the isodensity HI structural parameters have stronger correlations with star formation. Normalising the HI disc size by the optical effective radius (instead of the isophotal radius) produces positive correlations with stellar masses and stellar surface densities and removes the correlations with star formation. This is due to the effective and isodensity HI radii increasing with mass at similar rates while, in the optical, the effective radius increases slower than the isophotal radius. Our results demonstrate that with WALLABY we can begin to bridge the gap between small galaxy samples with high spatial resolution HI data and large, statistical studies using spatially unresolved, single-dish data., Comment: 16 page, 5 figures, accepted for publication in PASA

Context: To study the effects of environment on galaxies we use HI observations of galaxies in and around the cluster A2626. The cluster can effectively be divided in three different environments: the cluster itself, a group environment in the periphery of the cluster (we call it the Swarm) and substructure in the cluster itself. We use these to study the dependence of galaxy properties on environment. Aims: We have explored the relationship between HI deficiency, HI morphology, and star formation deficiency for the galaxies in and around the A2626 galaxy cluster to investigate the environmental effects on those properties. Methods: To quantify asymmetries of the outer HI disc of a galaxy, we used 1) three visual classes based on the outermost reliable HI contour (settled, disturbed, unsettled HI discs), 2) the offset between the HI centre and the optical centre of a galaxy, and 3) the modified asymmetry parameter Amod as defined by Lelli et al. (2014). Results: The HI deficiency of a galaxy is strongly correlated with the projected distance from the centre of A2626. Furthermore, substructure galaxies tend to be more asymmetric than the isolated galaxies in A2626, plausibly because of more efficient tidal interactions within substructures than outside substructures. Moreover, asymmetric, offset, and smaller HI discs are not necessarily the result of the cluster environment, as they are also observed in substructures in A2626 and in the Swarm. This signifies that "pre-processing" of the HI discs of galaxies in groups or substructures plays an important role, together with the "processing" in the cluster environment. Finally, the galaxies in all three environments have slightly lower star-formation rates (SFRs) than the typical SFR for normal galaxies as manifested by their offset from the star formation main sequence, implying effective gas removal mechanisms in all three environments., Comment: 19 pages, 12 figures

The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed to familiarise the scientific community with SKAO data and to drive the development of new analysis techniques. We present the results from Science Data Challenge 2 (SDC2), which invited participants to find and characterise 233245 neutral hydrogen (Hi) sources in a simulated data product representing a 2000~h SKA MID spectral line observation from redshifts 0.25 to 0.5. Through the generous support of eight international supercomputing facilities, participants were able to undertake the Challenge using dedicated computational resources. Alongside the main challenge, `reproducibility awards' were made in recognition of those pipelines which demonstrated Open Science best practice. The Challenge saw over 100 participants develop a range of new and existing techniques, with results that highlight the strengths of multidisciplinary and collaborative effort. The winning strategy -- which combined predictions from two independent machine learning techniques to yield a 20 percent improvement in overall performance -- underscores one of the main Challenge outcomes: that of method complementarity. It is likely that the combination of methods in a so-called ensemble approach will be key to exploiting very large astronomical datasets., Comment: Under review by MNRAS; 28 pages, 16 figures

We present HI-based B- and R-band Tully-Fisher relations (TFRs) and the Baryonic TFR (BTFR) at z=0.2 using direct HI detections from the Blind Ultra-Deep HI Environmental Survey (BUDHIES). Deep photometry from the Isaac Newton Telescope was used for 36 out of 166 HI sources, matching the quality criteria required for a robust TFR analysis. Two velocity definitions at 20% and 50% of the peak flux were measured from the global HI profiles and adopted as proxies for the circular velocities. We compare our results with an identically constructed z=0 TFR from the Ursa Major (UMa) association of galaxies. To ensure an unbiased comparison of the TFR, all the samples were treated identically regarding sample selection and applied corrections. We provide catalogues and an atlas showcasing the properties of the galaxies. Our analysis is focused on the zero points of the TFR and BTFR with their slopes fixed to the z=0 relation. Our main results are: (1) The BUDHIES galaxies show more asymmetric HI profiles with shallower wings compared to the UMa galaxies, which is likely due to the environment in which they reside, (2) The luminosity-based z=0.2 TFRs are brighter and bluer than the z=0 TFRs, even when cluster galaxies are excluded from the BUDHIES sample, (3) The BTFR shows no evolution in its zero point over the past 2.5 billion years and does not significantly change on the inclusion of cluster galaxies, and (4) proper sample selection and consistent corrections are crucial for an unbiased analysis of the evolution of the TFR., Comment: 24 pages, 14 figures, published in MNRAS

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 analyze 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 are 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 is robust and effective enough for further calibration and (continuum) imaging of these data. Analysis of 304 observations show 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., Comment: 15 pages, 10 figures. Accepted for publication in A&A

Context. The present study addresses a key question related to our understanding of the relation between void galaxies and their environment: the relationship between luminous and dark matter in and around voids. Aims. To explore the extent to which local Universe voids are empty of matter, we study the full (dark+luminous) matter content of seven nearby cosmic voids that are fully contained within the CosmicFlows-3 volume. Methods. We obtained the matter-density profiles of seven cosmic voids using two independent methods. These were built from the galaxy redshift space two-point correlation function in conjunction with peculiar velocity gradients from the CosmicFlows-3 dataset. Results. The results are striking, because when the redshift survey is used, all voids show a radial positive gradient of galaxies, while based on the dynamical analysis, only three of these voids display a clear underdensity of matter in their center. Conclusions. This work constitutes the most detailed observational analysis of voids conducted so far, and shows that void emptiness should be derived from dynamical information. From this limited study, the Hercules void appears to be the best candidate for a local Universe pure "pristine volume", expanding in three directions with no dark matter located in that void., Comment: Submitted A\&A Nov 29, 2022 - AA/2022/45578 / Accepted March 3rd, 2023

We present the Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) Pilot Phase I HI kinematic models. This first data release consists of HI observations of three fields in the direction of the Hydra and Norma clusters, and the NGC 4636 galaxy group. In this paper, we describe how we generate and publicly release flat-disk tilted-ring kinematic models for 109/592 unique HI detections in these fields. The modelling method adopted here - which we call the WALLABY Kinematic Analysis Proto-Pipeline (WKAPP) and for which the corresponding scripts are also publicly available - consists of combining results from the homogeneous application of the FAT and 3DBAROLO algorithms to the subset of 209 detections with sufficient resolution and S/N in order to generate optimized model parameters and uncertainties. The 109 models presented here tend to be gas rich detections resolved by at least 3-4 synthesized beams across their major axes, but there is no obvious environmental bias in the modelling. The data release described here is the first step towards the derivation of similar products for thousands of spatially-resolved WALLABY detections via a dedicated kinematic pipeline. Such a large publicly available and homogeneously analyzed dataset will be a powerful legacy product that that will enable a wide range of scientific studies., Comment: Accepted to PASA

We present WALLABY pilot data release 1, the first public release of HI pilot survey data from the Wide-field ASKAP L-band Legacy All-sky Blind Survey (WALLABY) on the Australian Square Kilometre Array Pathfinder. Phase 1 of the WALLABY pilot survey targeted three $60~{\rm deg}^2$ regions on the sky in the direction of the Hydra and Norma galaxy clusters and the NGC 4636 galaxy group, covering the redshift range of z < 0.08. The source catalogue, images and spectra of nearly 600 extragalactic HI detections and kinematic models for 109 spatially resolved galaxies are available. As the pilot survey targeted regions containing nearby group and cluster environments, the median redshift of the sample of z ~ 0.014 is relatively low compared to the full WALLABY survey. The median galaxy HI mass is $2.3 \times 10^{9}~M_{\odot}$. The target noise level of 1.6 mJy per $30''$ beam and 18.5 kHz channel translates into a $5\sigma$ HI mass sensitivity for point sources of about $5.2 \times 10^{8} \, (D_{\rm L} / \mathrm{100~Mpc})^{2} \, M_{\odot}$ across 50 spectral channels (~200 km/s) and a $5\sigma$ HI column density sensitivity of about $8.6 \times 10^{19} \, (1 + z)^{4}~\mathrm{cm}^{-2}$ across 5 channels (~20 km/s) for emission filling the $30''$ beam. As expected for a pilot survey, several technical issues and artefacts are still affecting the data quality. Most notably, there are systematic flux errors of up to several 10% caused by uncertainties about the exact size and shape of each of the primary beams as well as the presence of sidelobes due to the finite deconvolution threshold. In addition, artefacts such as residual continuum emission and bandpass ripples have affected some of the data. The pilot survey has been highly successful in uncovering such technical problems, most of which are expected to be addressed and rectified before the start of the full WALLABY survey., Comment: 24 pages, 13 figures, 5 tables, accepted for publication in PASA

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 such accretion can occur in both a dynamically hot and cold mode. However, until now observational evidence of the accretion required to match the observed star formation histories is lacking. This paper attempts to determine whether galaxies in the local Universe possess a significant reservoir of HI and what would be the accretion rates derived from such reservoirs. We search the vicinity of 22 nearby galaxies for isolated HI clouds or distinct streams in a systematic and automated manner. The HALOGAS observations represent one of the most sensitive and detailed HI surveys to date. These observations typically reach column density sensitivities of 10^19 cm^-2 over a 20 km/s width. We find 14 secure HI cloud candidates without an observed optical counterpart. These cloud candidates appear to be analogues to the most massive clouds detected around the Milky Way and M31. However, on average their numbers seem significantly reduced. We constrain upper limits for HI accretion in the local Universe. The average HI mass currently observed amounts to a rate of 0.05 Msun/yr with a stringent upper limit of 0.22 Msun/yr, confirming previous estimates. This is much lower than the average star formation rate in this sample. Our best estimate, based on GBT detection limits of several galaxies, suggests that another 0.04 Msun/yr could be accreted from undetected clouds and streams. These results show that in nearby galaxies HI is not being accreted at the same rate as stars are currently being formed. Our study can not exclude that other forms of gas accretion are at work. However, these observations also do not reveal extended neutral gas reservoirs around most nearby spiral galaxies., Comment: Accepted for publication in Astronomy & Astrophysics section 4. Extragalactic astronomy. Data available at https://www.astron.nl/halogas/data.php

We present MeerKAT HI observations of six jellyfish candidate galaxies (JFCGs) in the galaxy cluster, A2626. Two of the six galaxies JW100 and JW103, that were identified as JFCGs from B-band images, are confirmed as jellyfish galaxies (JFGs). Both of the JFGs have low HI content, reside in the cluster core, and move at very high velocities ($\sim$ 3$\sigma_{cl}$). The other JFCGs, identified as non-jellyfish galaxies, are HI rich, with HI morphologies revealing warps, asymmetries, and possible tidal interactions. Both the A2626 JFGs and three other confirmed JFGs from the GASP sample show that these galaxies are HI stripped but not yet quenched. We detect HI, Halpha, and CO tails of similar extent ($\sim$ 50 kpc) in JW100. Comparing the multi-phase velocity channels, we do not detect any HI or CO emission in the northern section of the tail where Halpha emission is present, possibly due to prolonged interaction between the stripped gas and the ICM. We also observe an anti-correlation between HI and CO, which hints at an efficient conversion of HI to H2 in the southern part of the tail. We find that both RPS and HI-to-H2 conversion are significant depletion channels for atomic gas. HI-to-H2 conversion is more efficient in the disc than in the tail., Comment: 15 pages, 10 figures, accepted for publication in MNRAS

(Abridged) Apertif is a phased-array feed system for WSRT, providing forty instantaneous beams over 300 MHz of bandwidth. A dedicated survey program started on 1 July 2019, with the last observations taken on 28 February 2022. We describe the release of data products from the first year of survey operations, through 30 June 2020. We focus on defining quality control metrics for the processed data products. 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. 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 are released. The median noise in the continuum images is 41.4 uJy/bm, with a slightly lower median noise of 36.9 uJy/bm in the Stokes V polarization image. The median angular resolution is 11.6"/sin(Dec). The median noise for all line cubes, with a spectral resolution of 36.6 kHz, is 1.6 mJy/bm, corresponding to a 3-sigma HI column density sensitivity of 1.8 x 10^20 atoms cm^-2 over 20 km/s (for a median angular resolution of 24" x 15"). We also provide primary beam images for each individual Apertif compound beam. The data are made accessible using a Virtual Observatory interface., Comment: Accepted for publication in A&A, updated Figure 1

The first data release of Apertif survey contains 3074 radio continuum images covering a thousand square degrees of the sky. The observations were performed during August 2019 to July 2020. The continuum images were produced at a central frequency 1355 MHz with the bandwidth of $\sim$150 MHz and angular resolution reaching 10". In this work we introduce and apply a new method to obtain a primary beam model using a machine learning approach, Gaussian process regression. The primary beam models obtained with this method are published along with the data products for the first Apertif data release. We apply the method to the continuum images, mosaic them and extract the source catalog. The catalog contains 249672 radio sources many of which are detected for the first time at these frequencies. We cross-match the coordinates with the NVSS, LOFAR/DR1/value-added and LOFAR/DR2 catalogs resulting in 44523, 22825 and 152824 common sources respectively. The first sample provides a unique opportunity to detect long term transient sources which have significantly changed their flux density for the last 25 years. The second and the third ones combined together provide information about spectral properties of the sources as well as the redshift estimates., Comment: 12 pages, 9 figures; accepted for publication in A&A

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., Comment: Accepted for publication in A&A

Context. Phased Array Feeds (PAFs) are multi element receivers in the focal plane of a telescope that make it possible to form simultaneously 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 and carried out several observing programs including two imaging surveys and a time domain survey. The Apertif imaging surveys use a configuration, where 40 partially overlapping compound beams (CBs) are simultaneously formed on the sky and arranged in an approximately rectangular shape. Aims. This manuscript aims to characterise 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 image deconvolution, primary beam correction and mosaicing of Apertif imaging data. Methods. We use drift scan measurements to measure the response of each of the 40 CBs of Apertif. We derive 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 sample 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 derive beam maps for each individual WSRT dish as well. Results. We present the frequency and time dependence of the beam shapes and sizes. We compare 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., Comment: Accepted for publication by A&A, 14 pages, 15 figures

We investigate the effect of observational constraints such as signal-to-noise, resolution and column density level on the HI morphological asymmetry ($\mathrm{A}_\mathrm{mod}$) and the effect of noise on the HI global profile ($\mathrm{A}_\mathrm{flux}$) asymmetry indices. Using mock galaxies from the EAGLE simulations we find an optimal combination of the observational constraints that are required for robust measurement of the $\mathrm{A}_\mathrm{mod}$ value of a galaxy: a column density threshold of $5\times10^{19}cm^{-2}$ or lower at a minimal signal-to-noise of 3 and a galaxy resolved with at least 11 beams. We also use mock galaxies to investigate the effect of noise on the $\mathrm{A}_\mathrm{flux}$ values and conclude that a global profile with signal-to-noise ratio greater than 5.5 is required to achieve a robust measurement of asymmetry. We investigate the relation between $\mathrm{A}_\mathrm{mod}$ and $\mathrm{A}_\mathrm{flux}$ indices and find them to be uncorrelated which implies that $\mathrm{A}_\mathrm{flux}$ values cannot be used to predict morphological asymmetries in galaxies., Comment: 19 pages, 12 figures and 2 appendices. Accepted for publication in the Monthly Notices of the Royal Astronomical Soceity

We present new results of a 12CO(J=1-0) imaging survey using the Atacama Compact Array (ACA) for 31 HI detected galaxies in the IC 1459 and NGC 4636 groups. This is the first CO imaging survey for loose galaxy groups. We obtained well-resolved CO data (~0.7-1.5 kpc) for a total of 16 galaxies in two environments. By comparing our ACA CO data with the HI and UV data, we probe the impacts of the group environment on the cold gas components (CO and HI gas) and star formation activity. We find that CO and/or HI morphologies are disturbed in our group members, some of which show highly asymmetric CO distributions (e.g., IC 5264, NGC 7421, and NGC 7418). In comparison with isolated galaxies in the xCOLD GASS sample, our group galaxies tend to have low star formation rates and low H2 gas fractions. Our findings suggest that the group environment can change the distribution of cold gas components, including the molecular gas, and star formation properties of galaxies. This is supporting evidence that preprocessing in the group-like environment can play an important role in galaxy evolution., Comment: 42 pages, 29 figures, 6 tables, published in ApJS

We analyse five early science datasets from the APERture Tile in Focus (Apertif) phased array feed system to verify the polarisation capabilities of Apertif in view of future larger data releases. We aim to characterise the source population of the polarised sky in the L-Band using polarised source information in combination with IR and optical data. We use automatic routines to generate full field-of-view Q- and U-cubes and perform RM-Synthesis, source finding, and cross-matching with published radio, optical, and IR data to generate polarised source catalogues. SED-fitting routines were used to determine photometric redshifts, star-formation rates, and galaxy masses. IR colour information was used to classify sources as AGN or star-forming-dominated and early- or late-type. We surveyed an area of 56deg$^2$ and detected 1357 polarised source components in 1170 sources. The fraction of polarised sources is 10.57% with a median fractional polarisation of 4.70$\pm$0.14%. We confirmed the reliability of the Apertif measurements by comparing them with polarised cross-identified NVSS sources. Average RMs of the individual fields lie within the error of the best Milky Way foreground measurements. All of our polarised sources were found to be dominated by AGN activity in the radio regime with most of them being radio-loud (79%) and of the FRII class (87%). The host galaxies of our polarised source sample are dominated by intermediate disc and star-forming disc galaxies. The contribution of star formation to the radio emission is on the order of a few percent for $\approx$10% of the polarised sources while for $\approx$90% it is completely dominated by the AGN. We do not see any change in fractional polarisation for different star-formation rates of the AGN host galaxies., Comment: 24 pages, 21 figures

Fast radio bursts (FRBs) are extragalactic radio transients of extraordinary luminosity. Studying the diverse temporal and spectral behaviour recently observed in a number of FRBs may help determine the nature of the entire class. For example, a fast spinning or highly magnetised neutron star might generate the rotation-powered acceleration required to explain the bright emission. Periodic, sub-second components, suggesting such rotation, were recently reported in one FRB, and potentially in two more. Here we report the discovery of FRB 20201020A with Apertif, an FRB showing five components regularly spaced by 0.415 ms. This sub-millisecond structure in FRB 20201020A carries important clues about the progenitor of this FRB specifically, and potentially about that of FRBs in general. We thus contrast its features to the predictions of the main FRB source models. We perform a timing analysis of the FRB 20201020A components to determine the significance of the periodicity. We compare these against the timing properties of the previously reported CHIME FRBs with sub-second quasi-periodic components, and against two Apertif bursts from repeating FRB 20180916B that show complex time-frequency structure. We find the periodicity of FRB 20201020A to be marginally significant at 2.5$\sigma$. Its repeating subcomponents cannot be explained as a pulsar rotation since the required spin rate of over 2 kHz exceeds the limits set by typical neutron star equations of state and observations. The fast periodicity is also in conflict with a compact object merger scenario. These quasi-periodic components could, however, be caused by equidistant emitting regions in the magnetosphere of a magnetar. The sub-millisecond spacing of the components in FRB 20201020A, the smallest observed so far in a one-off FRB, may rule out both neutron-star rotation and binary mergers as the direct source of quasi-periodic FRBs., Comment: 13 pages, 6 figures, 3 tables, supplementary material. Submitted to A&A

Radio Astronomy is undergoing a renaissance, as the next-generation of instruments provides a massive leap forward in collecting area and therefore raw sensitivity. However, to achieve this theoretical level of sensitivity in the science data products we need to address the much more pernicious systematic effects, which are the true limitation. These become all the more significant when we consider that much of the time used by survey instruments, such as the SKA, will be dedicated to deep surveys. CHILES is a deep HI survey of the COSMOS field, with 1,000 hours of VLA time. We present our approach for creating the image cubes from the first Epoch, with discussions of the methods and quantification of the data quality from 946 to 1420MHz -- a redshift range of 0.5 to 0. We layout the problems we had to solve and describe how we tackled them. These are of importance as CHILES is the first deep wideband multi-epoch HI survey and it has relevance for ongoing and future surveys. We focus on the accumulated systematic errors in the imaging, as the goal is to deliver a high-fidelity image that is only limited by the random thermal errors. To understand and correct these systematic effects we ideally manage them in the domain in which they arise, and that is predominately the visibility domain. CHILES is a perfect test bed for many of the issues we can expect for deep imaging with the SKA or ngVLA and we discuss the lessons we have learned., Comment: Accepted in AJ

Apertif (APERture Tile In Focus) is one of the Square Kilometre Array (SKA) pathfinder facilities. The Apertif project is an upgrade to the 50-year-old Westerbork Synthesis Radio Telescope (WSRT) using phased-array feed technology. The new receivers create 40 individual beams on the sky, achieving an instantaneous sky coverage of 6.5 square degrees. The primary goal of the Apertif Imaging Survey is to perform a wide survey of 3500 square degrees (AWES) and a medium deep survey of 350 square degrees (AMES) of neutral atomic hydrogen (up to a redshift of 0.26), radio continuum emission and polarisation. Each survey pointing yields 4.6 TB of correlated data. The goal of Apercal is to process this data and fully automatically generate science ready data products for the astronomical community while keeping up with the survey observations. We make use of common astronomical software packages in combination with Python based routines and parallelisation. We use an object oriented module-based approach to ensure easy adaptation of the pipeline. A Jupyter notebook based framework allows user interaction and execution of individual modules as well as a full automatic processing of a complete survey observation. If nothing interrupts processing, we are able to reduce a single pointing survey observation on our five node cluster with 24 physical cores and 256 GB of memory each within 24h keeping up with the speed of the surveys. The quality of the generated images is sufficient for scientific usage for 44 % of the recorded data products with single images reaching dynamic ranges of several thousands. Future improvements will increase this percentage to over 80 %. Our design allowed development of the pipeline in parallel to the commissioning of the Apertif system.

We present results from our analysis of the Hydra I cluster observed in neutral atomic hydrogen (HI) as part of the Widefield ASKAP L-band Legacy All-sky Blind Survey (WALLABY). These WALLABY observations cover a 60-square-degree field of view with uniform sensitivity and a spatial resolution of 30 arcsec. We use these wide-field observations to investigate the effect of galaxy environment on HI gas removal and star formation quenching by comparing the properties of cluster, infall and field galaxies extending up to $\sim5R_{200}$ from the cluster centre. We find a sharp decrease in the HI-detected fraction of infalling galaxies at a projected distance of $\sim1.5R_{200}$ from the cluster centre from $\sim0.85\%$ to $\sim0.35\%$. We see evidence for the environment removing gas from the outskirts of HI-detected cluster and infall galaxies through the decrease in the HI to $r$-band optical disc diameter ratio. These galaxies lie on the star forming main sequence, indicating that gas removal is not yet affecting the inner star-forming discs and is limited to the galaxy outskirts. Although we do not detect galaxies undergoing galaxy-wide quenching, we do observe a reduction in recent star formation in the outer disc of cluster galaxies, which is likely due to the smaller gas reservoirs present beyond the optical radius in these galaxies. Stacking of HI non-detections with HI masses below $M_{\rm{HI}}\lesssim10^{8.4}\,\rm{M}_{\odot}$ will be required to probe the HI of galaxies undergoing quenching at distances $\gtrsim60$ Mpc with WALLABY., Comment: 18 pages, 10 figures, accepted for publication in MNRAS

We present the first molecular gas mass survey of void galaxies. We compare these new data together with data for the atomic gas mass and star formation rate ($\rm SFR$) from the literature to those of galaxies in filaments and walls in order to better understand how molecular gas and star formation are related to the large-scale environment. We observed at the IRAM 30 m telescope the CO(1-0) and CO(2-1) emission of 20 void galaxies selected from the Void Galaxy Survey (VGS), with a stellar mass range from $\rm 10^{8.5}$ to $\rm 10^{10.3}M_{\odot}$. We detected 15 objects in at least one CO line. We compared the molecular gas mass ($M_{\rm H_2}$), the star formation efficiency ($\rm SFE =SFR/M_{\rm H_2}$), the atomic gas mass, the molecular-to-atomic gas mass ratio, and the specific star formation rate (sSFR) of the void galaxies with two control samples of galaxies in filaments and walls, selected from xCOLD GASS and EDGE-CALIFA, for different stellar mass bins and taking the star formation activity into account. The results for the molecular gas mass for a sample of 20 voids galaxies allowed us to make a statistical comparison to galaxies in filaments and walls for the first time., Comment: 24 pages, 15 figures, accepted in A&A, language corrected version

We describe the APERture Tile In Focus (Apertif) system, a phased array feed (PAF) upgrade of the Westerbork Synthesis Radio Telescope which has transformed this telescope into a high-sensitivity, wide field-of-view L-band imaging and transient survey instrument. Using novel PAF technology, up to 40 partially overlapping beams can be formed on the sky simultaneously, significantly increasing the survey speed of the telescope. With this upgraded instrument, an imaging survey covering an area of 2300 deg2 is being performed which will deliver both continuum and spectral line data sets, of which the first data has been publicly released. In addition, a time domain transient and pulsar survey covering 15,000 deg2 is in progress. An overview of the Apertif science drivers, hardware and software of the upgraded telescope is presented, along with its key performance characteristics., Comment: 29 pages, 42 figures, accepted for publication by A&A

We present SoFiA 2, the fully automated 3D source finding pipeline for the WALLABY extragalactic HI survey with the Australian SKA Pathfinder (ASKAP). SoFiA 2 is a reimplementation of parts of the original SoFiA pipeline in the C programming language and makes use of OpenMP for multi-threading of the most time-critical algorithms. In addition, we have developed a parallel framework called SoFiA-X that allows the processing of large data cubes to be split across multiple computing nodes. As a result of these efforts, SoFiA 2 is substantially faster and comes with a much reduced memory footprint compared to its predecessor, thus allowing the large WALLABY data volumes of hundreds of gigabytes of imaging data per epoch to be processed in real-time. The source code has been made publicly available to the entire community under an open-source licence. Performance tests using mock galaxies injected into genuine ASKAP data suggest that in the absence of significant imaging artefacts SoFiA 2 is capable of achieving near-100% completeness and reliability above an integrated signal-to-noise ratio of about 5-6. We also demonstrate that SoFiA 2 generally recovers the location, integrated flux and w20 line width of galaxies with high accuracy. Other parameters, including the peak flux density and w50 line width, are more strongly biased due to the influence of the noise on the measurement. In addition, very faint galaxies below an integrated signal-to-noise ratio of about 10 may get broken up into multiple components, thus requiring a strategy to identify fragmented sources and ensure that they do not affect the integrity of any scientific analysis based on the SoFiA 2 output., Comment: 16 pages, 10 figures, 1 table, accepted for publication in MNRAS

This study uses HI image data from the WALLABY pilot survey with the ASKAP telescope, covering the Hydra cluster out to 2.5$r_{200}$. We present the projected phase-space distribution of HI-detected galaxies in Hydra, and identify that nearly two thirds of the galaxies within $1.25r_{200}$ may be in the early stages of ram pressure stripping. More than half of these may be only weakly stripped, with the ratio of strippable HI (i.e., where the galactic restoring force is lower than the ram pressure in the disk) mass fraction (over total HI mass) distributed uniformly below 90%. Consequently, the HI mass is expected to decrease by only a few 0.1 dex after the currently strippable portion of HI in these systems has been stripped. A more detailed look at the subset of galaxies that are spatially resolved by WALLABY observations shows that, while it typically takes less than 200 Myr for ram pressure stripping to remove the currently strippable portion of HI, it may take more than 600 Myr to significantly change the total HI mass. Our results provide new clues to understanding the different rates of HI depletion and star formation quenching in cluster galaxies., Comment: 25 pages, 10 figures, 1 table. Accepted for publication at ApJ

Detection of the electromagnetic emission from coalescing binary neutron stars (BNS) is important for understanding the merger and afterglow. We present a search for a radio counterpart to the gravitational-wave source GW190425, a BNS merger, using Apertif on the Westerbork Synthesis Radio Telescope (WSRT). We observe a field of high probability in the associated localisation region for 3 epochs at 68, 90 and 109 days post merger. We identify all sources that exhibit flux variations consistent with the expected afterglow emission of GW190425. We also look for possible transients. These are sources which are only present in one epoch. In addition, we quantify our ability to search for radio afterglows in fourth and future observing runs of the gravitational-wave detector network using Monte Carlo simulations. We found 25 afterglow candidates based on their variability. None of these could be associated with a possible host galaxy at the luminosity distance of GW190425. We also found 55 transient afterglow candidates that were only detected in one epoch. All turned out to be image artefacts. In the fourth observing run, we predict that up to three afterglows will be detectable by Apertif. While we did not find a source related to the afterglow emission of GW190425, the search validates our methods for future searches of radio afterglows., Comment: 11 pages, 7 figures, accepted for publication

We present the serendipitous detection of the two main OH maser lines at 1667 and 1665 MHz associated with IRAS 10597+5926 at z = 0.19612 in the untargeted Apertif Wide-area Extragalactic Survey (AWES), and the subsequent measurement of the OH 1612 MHz satellite line in the same source. With a total OH luminosity of log(L/L_Sun) = 3.90 +/- 0.03, IRAS 10597+5926 is the fourth brightest OH megamaser (OHM) known. We measure a lower limit for the 1667/1612 ratio of R_1612 > 45.9 which is the highest limiting ratio measured for the 1612 MHz OH satellite line to date. OH satellite line measurements provide a potentially valuable constraint by which to compare detailed models of OH maser pumping mechanisms. Optical imaging shows the galaxy is likely a late-stage merger. Based on published infrared and far ultraviolet fluxes, we find that the galaxy is an ultra luminous infrared galaxy (ULIRG) with log(L_TIR/L_Sun) = 12.24, undergoing a star burst with an estimated star formation rate of 179 +/- 40 M_Sun/yr. These host galaxy properties are consistent with the physical conditions responsible for very bright OHM emission. Finally, we provide an update on the predicted number of OH masers that may be found in AWES, and estimate the total number of OH masers that will be detected in each of the individual main and satellite OH 18 cm lines., Comment: 9 pages, 4 figures. Accepted for publication in A&A

Scientists across all disciplines increasingly rely on machine learning algorithms to analyse and sort datasets of ever increasing volume and complexity. Although trends and outliers are easily extracted, careful and close inspection will still be necessary to explore and disentangle detailed behavior, as well as identify systematics and false positives. We must therefore incorporate new technologies to facilitate scientific analysis and exploration. Astrophysical data is inherently multi-parameter, with the spatial-kinematic dimensions at the core of observations and simulations. The arrival of mainstream virtual-reality (VR) headsets and increased GPU power, as well as the availability of versatile development tools for video games, has enabled scientists to deploy such technology to effectively interrogate and interact with complex data. In this paper we present development and results from custom-built interactive VR tools, called the iDaVIE suite, that are informed and driven by research on galaxy evolution, cosmic large-scale structure, galaxy-galaxy interactions, and gas/kinematics of nearby galaxies in survey and targeted observations. In the new era of Big Data ushered in by major facilities such as the SKA and LSST that render past analysis and refinement methods highly constrained, we believe that a paradigm shift to new software, technology and methods that exploit the power of visual perception, will play an increasingly important role in bridging the gap between statistical metrics and new discovery. We have released a beta version of the iDaVIE software system that is free and open to the community., Comment: Accepted for publication in the Journal of Astronomy and Computing ; https://vislab.idia.ac.za/ ; key words: Virtual Reality, data visualization, radio astrophysics, 3D catalogues, volumetric rendering; neutral hydrogen in galaxies

We determined the HI mass function of galaxies in the Ursa Major association of galaxies using a blind VLA-D array survey, consisting of 54 pointings in a cross pattern, covering the centre as well as the outskirts of the Ursa Major volume. The calculated HI mass function has best-fitting Schechter parameters {\theta}^* = 0.19+/-0.11 Mpc^{-3}, log(M^*_{HI}/M_{\odot}) = 9.8+/-0.8 and {\alpha} = -0.92+/-0.16. The high-mass end is determined by a complementary, targeted WSRT survey, the low-mass end is determined by the blind VLA survey. The slope is significantly shallower than the slopes of the HIPASS ({\alpha} = -1.37+/-0.03+/-0.05) and ALFALFA ({\alpha} = -1.33+/-0.02) HI mass functions, which are measured over much larger volumes and cover a wider range of cosmic environments: There is a relative lack of low HI mass galaxies in the Ursa Major region. This difference in the slope strongly hints at an environmental dependence of the HI mass function slope., Comment: 19 pages, 15 figures. Accepted for publication in MNRAS

Supermassive black holes at the centres of galaxies can cycle through periods of activity and quiescence. Characterising the duty cycle of active galactic nuclei is crucial for understanding the impact of the energy they release on the host galaxy. For radio AGN, this can be done by identifying dying (remnant) and restarted radio galaxies from their radio spectral properties. Using the combination of images at 1400 MHz produced by Apertif, the new phased-array feed receiver installed on the Westerbork Synthesis Radio Telescope, and images at 150 MHz provided by LOFAR, we have derived resolved spectral index images (at a resolution of ~15 arcsec) for all the sources within ~6 deg^2 area of the Lockman Hole region. We were able to select 15 extended radio sources with emission (partly or entirely) characterised by extremely steep spectral indices (steeper than 1.2). These objects represent radio sources in the remnant or the restarted phases of their life cycle. Our findings suggest this cycle to be relatively fast. They also show a variety of properties relevant for modelling the evolution of radio galaxies. For example, the restarted activity can occur while the remnant structure from a previous phase of activity is still visible. This provides constraints on the duration of the 'off' (dying) phase. In extended remnants with ultra-steep spectra at low frequencies, the activity likely stopped a few hundred megayears ago, and they correspond to the older tail of the age distribution of radio galaxies, in agreement with simulations of radio source evolution. We find remnant radio sources with a variety of structures (from double-lobed to amorphous), suggesting different types of progenitors. The present work sets the stage for exploiting low-frequency spectral index studies of extended sources by taking advantage of the large areas common to the LOFAR and the Apertif surveys., Comment: Accepted for publication in A&A. This paper is part of the 1st data release of the LoTSS Deep Fields. 17 pages, 10 figures

Galaxy clusters are some of largest structures in the universe. These very dense environments tend to be home to higher numbers of evolved galaxies that what is found in lower density environments. It is well known that dense environments can influence the evolution of galaxies through the removal of the neutral gas (HI) reservoirs which fuel star formation. It is unclear which environment has a stronger effect: the local environment (i.e. the substructure within the cluster), or the cluster itself. Using the new HI data from the Westerbork Coma Survey, we explore the average HI content of galaxies across the cluster comparing galaxies that reside in substructure to those that do not. We apply to the Dressler-Shectman test to our newly compiled redshift catalogue of the Coma cluster to search for substructure. With so few of the Coma galaxies directly detected in HI, we use the HI stacking technique to probe average HI content below what can be directly detected. Using the Dressler-Shectman test, we find 15 substructures within the footprint of the Westerbork Coma Survey. We compare the average HI content for galaxies within substructure to those not in substructure. Using the HI stacking technique, we find that the Coma galaxies (for which are not detected in HI) are more than 10--50 times more HI deficient than expected which supports the scenario of an extremely efficient and rapid quenching mechanism. By studying the galaxies that are not directly detected in HI, we also find Coma to be more HI deficient than previously thought., Comment: 18 pages (+ 21 page appendix), 23 figures, accepted for publication in A&A

The propagation of radio waves from distant compact radio sources through turbulent interstellar plasma in our Galaxy causes these sources to twinkle, a phenomenon called interstellar scintillation. Such scintillations are a unique probe of the micro-arcsecond structure of radio sources as well as of the sub-AU-scale structure of the Galactic interstellar medium. Weak scintillations (i.e. an intensity modulation of a few percent) on timescales of a few days or longer are commonly seen at centimetre wavelengths and are thought to result from the line-of-sight integrated turbulence in the interstellar plasma of the Milky Way. So far, only three sources were known that show more extreme variations, with modulations at the level of some dozen percent on timescales shorter than an hour. This requires propagation through nearby (d <~10 pc) anomalously dense (n_e ~10^2 cm^-3) plasma clouds. Here we report the discovery with Apertif of a source (J1402+5347) showing extreme (~50%) and rapid variations on a timescale of just 6.5 minutes in the decimetre band (1.4 GHz). The spatial scintillation pattern is highly anisotropic, with a semi-minor axis of about 20,000 km. The canonical theory of refractive scintillation constrains the scattering plasma to be within the Oort cloud. The sightline to J1402+5347, however, passes unusually close to the B3 star Alkaid (eta UMa) at a distance of 32 pc. If the scintillations are associated with Alkaid, then the angular size of J1402+5347 along the minor axis of the scintels must be smaller than ~10 micro arcsec yielding an apparent brightness temperature for an isotropic source of >~ 10^ 14K. }, Comment: Accepted for Astronomy and Astrophysics Letters

The Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) is a next-generation survey of neutral hydrogen (HI) in the Local Universe. It uses the widefield, high-resolution capability of the Australian Square Kilometer Array Pathfinder (ASKAP), a radio interferometer consisting of 36 x 12-m dishes equipped with Phased-Array Feeds (PAFs), located in an extremely radio-quiet zone in Western Australia. WALLABY aims to survey three-quarters of the sky (-90 degr < Dec < +30 degr) to a redshift of z < 0.26, and generate spectral line image cubes at ~30 arcsec resolution and ~1.6 mJy/beam per 4 km/s channel sensitivity. ASKAP's instantaneous field of view at 1.4 GHz, delivered by the PAF's 36 beams, is about 30 sq deg. At an integrated signal-to-noise ratio of five, WALLABY is expected to detect over half a million galaxies with a mean redshift of z ~ 0.05 (~200 Mpc). The scientific goals of WALLABY include: (a) a census of gas-rich galaxies in the vicinity of the Local Group; (b) a study of the HI properties of galaxies, groups and clusters, in particular the influence of the environment on galaxy evolution; and (c) the refinement of cosmological parameters using the spatial and redshift distribution of low-bias gas-rich galaxies. For context we provide an overview of previous large-scale HI surveys. Combined with existing and new multi-wavelength sky surveys, WALLABY will enable an exciting new generation of panchromatic studies of the Local Universe. - First results from the WALLABY pilot survey are revealed, with initial data products publicly available in the CSIRO ASKAP Science Data Archive (CASDA)., Comment: Accepted for publication in ApSS (38 pages, 14 figures), see also https://www.atnf.csiro.au/research/WALLABY/ - Contact email: Baerbel.Koribalski@csiro.au

Repeating fast radio bursts (FRBs) present excellent opportunities to identify FRB progenitors and host environments, as well as decipher the underlying emission mechanism. Detailed studies of repeating FRBs might also hold clues to the origin of FRBs as a population. We aim to detect the first two repeating FRBs: FRB 121102 (R1) and FRB 180814.J0422+73 (R2), and characterise their repeat statistics. We also want to significantly improve the sky localisation of R2. We use the Westerbork Synthesis Radio Telescope to conduct extensive follow-up of these two repeating FRBs. The new phased-array feed system, Apertif, allows covering the entire sky position uncertainty of R2 with fine spatial resolution in one pointing. We characterise the energy distribution and the clustering of detected R1 bursts. We detected 30 bursts from R1. Our measurements indicate a dispersion measure of 563.5(2) pc cm$^{-3}$, suggesting a significant increase in DM over the past few years. We place an upper limit of 8% on the linear polarisation fraction of the brightest burst. We did not detect any bursts from R2. A single power-law might not fit the R1 burst energy distribution across the full energy range or widely separated detections. Our observations provide improved constraints on the clustering of R1 bursts. Our stringent upper limits on the linear polarisation fraction imply a significant depolarisation, either intrinsic to the emission mechanism or caused by the intervening medium, at 1400 MHz that is not observed at higher frequencies. The non-detection of any bursts from R2 implies either a highly clustered nature of the bursts, a steep spectral index, or a combination of both. Alternatively, R2 has turned off completely, either permanently or for an extended period of time., Comment: 11 pages, 7 figures, submitted to A&A

We present neutral hydrogen (HI) and ionized hydrogen (H${\alpha}$) observations of ten galaxies out to a redshift of 0.1. The HI observations are from the first epoch (178 hours) of the COSMOS HI Large Extragalactic Survey (CHILES). Our sample is HI biased and consists of ten late-type galaxies with HI masses that range from $1.8\times10^{7}$ M$_{\odot}$ to $1.1\times10^{10}$ M$_{\odot}$. We find that although the majority of galaxies show irregularities in the morphology and kinematics, they generally follow the scaling relations found in larger samples. We find that the HI and H${\alpha}$ velocities reach the flat part of the rotation curve. We identify the large-scale structure in the nearby CHILES volume using DisPerSE with the spectroscopic catalog from SDSS. We explore the gaseous properties of the galaxies as a function of location in the cosmic web. We also compare the angular momentum vector (spin) of the galaxies to the orientation of the nearest cosmic web filament. Our results show that galaxy spins tend to be aligned with cosmic web filaments and show a hint of a transition mass associated with the spin angle alignment., Comment: 24 pages, 25 figures, 6 tables, accepted for publication in MNRAS

We present a study of 16 HI-detected galaxies found in 178 hours of observations from Epoch 1 of the COSMOS HI Large Extragalactic Survey (CHILES). We focus on two redshift ranges between 0.108 <= z <= 0.127 and 0.162 <= z <= 0.183 which are among the worst affected by radio frequency interference (RFI). While this represents only 10% of the total frequency coverage and 18% of the total expected time on source compared to what will be the full CHILES survey, we demonstrate that our data reduction pipeline recovers high quality data even in regions severely impacted by RFI. We report on our in-depth testing of an automated spectral line source finder to produce HI total intensity maps which we present side-by-side with significance maps to evaluate the reliability of the morphology recovered by the source finder. We recommend that this become a common place manner of presenting data from upcoming HI surveys of resolved objects. We use the COSMOS 20k group catalogue, and we extract filamentary structure using the topological DisPerSE algorithm to evaluate the \hi\ morphology in the context of both local and large-scale environments and we discuss the shortcomings of both methods. Many of the detections show disturbed HI morphologies suggesting they have undergone a recent interaction which is not evident from deep optical imaging alone. Overall, the sample showcases the broad range of ways in which galaxies interact with their environment. This is a first look at the population of galaxies and their local and large-scale environments observed in HI by CHILES at redshifts beyond the z=0.1 Universe., Comment: 23 pages, 12 figures, 1 interactive 3D figure, accepted to MNRAS

Rotation measure (RM) synthesis is a widely used polarization processing algorithm for reconstructing polarized structures along the line of sight. Performing RM synthesis on large datasets produced by telescopes like LOFAR can be computationally intensive as the computational cost is proportional to the product of the number of input frequency channels, the number of output Faraday depth values to be evaluated and the number of lines of sight present in the data cube. The required computational cost is likely to get worse due to the planned large area sky surveys with telescopes like the Low Frequency Array (LOFAR), the Murchison Widefield Array (MWA), and eventually the Square Kilometre Array (SKA). The massively parallel General Purpose Graphical Processing Units (GPGPUs) can be used to execute some of the computationally intensive astronomical image processing algorithms including RM synthesis. In this paper, we present a GPU-accelerated code, called cuFFS or CUDA-accelerated Fast Faraday Synthesis, to perform Faraday rotation measure synthesis. Compared to a fast single-threaded and vectorized CPU implementation, depending on the structure and format of the data cubes, our code achieves an increase in speed of up to two orders of magnitude. During testing, we noticed that the disk I/O when using the Flexible Image Transport System (FITS) data format is a major bottleneck and to reduce the time spent on disk I/O, our code supports the faster HDFITS format in addition to the standard FITS format. The code is written in C with GPU-acceleration achieved using Nvidia's CUDA parallel computing platform. The code is available at https://github.com/sarrvesh/cuFFS., Comment: Accepted for publication in Astronomy & Computing

We present a 3x3 degrees, 105-pointing, high-resolution neutral hydrogen (HI) mosaic of the M81 galaxy triplet (including the galaxies M81, M82 and NGC 3077, as well as dwarf galaxy NGC 2976) obtained with the Very Large Array (VLA) C and D arrays. This uniformly covers the entire area and velocity range of the triplet with a resolution of ~20'' or ~420 pc. The data reveal many small-scale anomalous velocity features highlighting the complexity of the interacting M81 triplet. We compare our data with Green Bank Telescope (GBT) observations of the same area. This provides evidence for a substantial reservoir of low-column density gas in the northern part of the triplet, probably associated with M82. Such a reservoir is not found in the southern part. We report a number of kpc-sized low-mass HI clouds with HI masses of a few times 10^6 Msun. Their dynamical masses are much larger than their baryonic masses, which could indicate the presence of dark matter if the clouds are rotationally supported. However, due to their spatial and kinematical association with HI tidal features, it is more likely that the velocity widths indicate tidal effects or streaming motions. We do not find any clouds not associated with tidal features down to an HI mass limit of a few times 10^4 Msun. We compare the HI column densities with resolved stellar density maps and find a star formation threshold around 3-6 10^20 cm-2$. We find that extreme velocity dispersions can be explained by a superposition of multiple components along the line of sight near M81 as well as winds or outflows around M82. The velocity dispersions found are high enough that these processes could explain the linewidths of Damped-Lyman-alpha absorbers observed at high redshift., Comment: Accepted for publication in The Astronomical Journal. Full-resolution PDF available at http://www.astron.nl/~blok/M81data/m81data.pdf. Data can be downloaded from http://www.astron.nl/~blok/M81data

The detailed distribution and kinematics of the atomic and the CO-bright molecular hydrogen in the disc of the Milky Way inside the Solar circle are derived under the assumptions of axisymmetry and pure circular motions. We divide the Galactic disc into a series of rings, and assume that the gas in each ring is described by four parameters: its rotation velocity, velocity dispersion, midplane density and its scale height. We fit these parameters to the Galactic HI and CO (J=1-0) data by producing artificial HI and CO line-profiles and comparing them with the observations. Our approach allows us to fit all parameters to the data simultaneously without assuming a-priori a radial profile for one of the parameters. We present the distribution and kinematics of the HI and H2 in both the approaching (QIV) and the receding (QI) regions of the Galaxy. Our best-fit models reproduces remarkably well the observed HI and CO longitude-velocity diagrams up to a few degrees of distance from the midplane. With the exception of the innermost 2.5 kpc, QI and QIV show very similar kinematics. The rotation curves traced by the HI and H2 follow closely each other, flattening beyond R=6.5 kpc. Both the HI and the H2 surface densities show a) a deep depression at 0.5

Upcoming HI surveys will deliver such large datasets that automated processing using the full 3-D information to find and characterize HI objects is unavoidable. Full 3-D visualization is an essential tool for enabling qualitative and quantitative inspection and analysis of the 3-D data, which is often complex in nature. Here we present $\tt{SlicerAstro}$, an open-source extension of 3DSlicer, a multi-platform open source software package for visualization and medical image processing, which we developed for the inspection and analysis of HI spectral line data. We describe its initial capabilities, including 3-D filtering, 3-D selection and comparative modelling.

SKA precursors are capable of detecting hundreds of galaxies in HI in a single 12 hours pointing. In deeper surveys one will probe more easily faint HI structures, typically located in the vicinity of galaxies, such as tails, filaments, and extraplanar gas. The importance of interactive visualization has proven to be fundamental for the exploration of such data as it helps users to receive immediate feedback when manipulating the data. We have developed SlicerAstro, a 3-D interactive viewer with new analysis capabilities, based on traditional 2-D input/output hardware. These capabilities enhance the data inspection, allowing faster analysis of complex sources than with traditional tools. SlicerAstro is an open-source extension of 3DSlicer, a multi-platform open source software package for visualization and medical image processing. We demonstrate the capabilities of the current stable binary release of SlicerAstro, which offers the following features: i) handling of FITS files and astronomical coordinate systems; ii) coupled 2-D/3-D visualization; iii) interactive filtering; iv) interactive 3-D masking; v) and interactive 3-D modeling. In addition, SlicerAstro has been designed with a strong, stable and modular C++ core, and its classes are also accessible via Python scripting, allowing great flexibility for user-customized visualization and analysis tasks., Comment: 18 pages, 11 figures, Accepted by Astronomy and Computing. SlicerAstro link: https://github.com/Punzo/SlicerAstro/wiki#get-slicerastro

We analyze photometry from deep B-band images of 59 void galaxies in the Void Galaxy Survey (VGS), together with their near-infrared 3.6$\mu$m and 4.5$\mu$m Spitzer photometry. The VGS galaxies constitute a sample of void galaxies that were selected by a geometric-topological procedure from the SDSS DR7 data release, and which populate the deep interior of voids. Our void galaxies span a range of absolute B-magnitude from $\rm{M_B=-15.5}$ to $\rm{M_B=-20}$, while at the 3.6$\mu$m band their magnitudes range from $\rm{M_{3.6}=-18}$ to $\rm{M_{3.6}=-24}$. Their B-[3.6] colour and structural parameters indicate these are star forming galaxies. A good reflection of the old stellar population, the near-infrared band photometry also provide a robust estimate of the stellar mass, which for the VGS galaxies we confirm to be smaller than $3 \times 10^{10}$ M$_\odot$. In terms of the structural parameters and morphology, our findings align with other studies in that our VGS galaxy sample consists mostly of small late-type galaxies. Most of them are similar to Sd-Sm galaxies, although a few are irregularly shaped galaxies. The sample even includes two early-type galaxies, one of which is an AGN. Their S\'{e}rsic indices are nearly all smaller than $n=2$ in both bands and they also have small half-light radii. In all, we conclude that the principal impact of the void environment on the galaxies populating them mostly concerns their low stellar mass and small size., Comment: Accepted for publication in MNRAS, 16 pages, 11 figures

Soon to be operational HI survey instruments such as APERTIF and ASKAP will produce large datasets. These surveys will provide information about the HI in and around hundreds of galaxies with a typical signal-to-noise ratio of $\sim$ 10 in the inner regions and $\sim$ 1 in the outer regions. In addition, such surveys will make it possible to probe faint HI structures, typically located in the vicinity of galaxies, such as extra-planar-gas, tails and filaments. These structures are crucial for understanding galaxy evolution, particularly when they are studied in relation to the local environment. Our aim is to find optimized kernels for the discovery of faint and morphologically complex HI structures. Therefore, using HI data from a variety of galaxies, we explore state-of-the-art filtering algorithms. We show that the intensity-driven gradient filter, due to its adaptive characteristics, is the optimal choice. In fact, this filter requires only minimal tuning of the input parameters to enhance the signal-to-noise ratio of faint components. In addition, it does not degrade the resolution of the high signal-to-noise component of a source. The filtering process must be fast and be embedded in an interactive visualization tool in order to support fast inspection of a large number of sources. To achieve such interactive exploration, we implemented a multi-core CPU (OpenMP) and a GPU (OpenGL) version of this filter in a 3D visualization environment ($\tt{SlicerAstro}$)., Comment: 17 pages, 9 figures, 4 tables. Astronomy and Computing, accepted