Deep galaxy surveys have revealed that the global star formation rate (SFR) density in the Universe peaks at 1 < z < 2 and sharply declines towards z = 0. But a clear picture of the underlying processes, in particular the evolution of cold atomic (~100 K) and molecular gas phases, that drive such a strong evolution is yet to emerge. MALS is designed to use MeerKAT's L- and UHF-band receivers to carry out the most sensitive (N(HI)>10$^{19}$ cm$^{-2}$) dust-unbiased search of intervening HI 21-cm and OH 18-cm absorption lines at 0 < z < 2. This will provide reliable measurements of the evolution of cold atomic and molecular gas cross-sections of galaxies, and unravel the processes driving the steep evolution in the SFR density. The large sample of HI and OH absorbers obtained from the survey will (i) lead to tightest constraints on the fundamental constants of physics, and (ii) be ideally suited to probe the evolution of magnetic fields in disks of galaxies via Zeeman Splitting or Rotation Measure synthesis. The survey will also provide an unbiased census of HI and OH absorbers, i.e. cold gas associated with powerful AGNs (>10$^{24}$ W Hz$^{-1}$) at 0 < z < 2, and will simultaneously deliver a blind HI and OH emission line survey, and radio continuum survey. Here, we describe the MALS survey design, observing plan and the science issues to be addressed under various science themes.

We built a multi-wavelength dataset for galaxies from the Local Volume HI Survey (LVHIS), which comprises 82 galaxies. We also select a sub-sample of ten large galaxies for investigating properties in the galactic outskirts. The LVHIS sample covers nearly four orders of magnitude in stellar mass and two orders of magnitude in HI mass fraction (fHI). The radial distribution of HI gas with respect to the stellar disc is correlated with fHI but with a large scatter. We confirm the previously found correlations between the total HI mass and star formation rate (SFR), and between HI surface densities and SFR surface densities beyond R25. However, the former correlation becomes much weaker when the average surface densities rather than total mass or rate are considered, and the latter correlation also becomes much weaker when the effect of stellar mass is removed or controlled. Hence the link between SFR and HI is intrinsically weak in these regions, consistent with what was found on kpc scales in the galactic inner regions. We find a strong correlation between the SFR surface density and the stellar mass surface density, which is consistent with the star formation models where the gas is in quasi-equilibrium with the mid-plane pressure. We find no evidence for HI warps to be linked with decreasing star forming efficiencies.

We examine the properties of atomic hydrogen (HI) associated with galaxies in the EAGLE simulations of galaxy formation. EAGLE's feedback parameters were calibrated to reproduce the stellar mass function and galaxy sizes at $z=0.1$, and we assess whether this calibration also yields realistic HI properties. We estimate the self-shielding density with a fitting function calibrated using radiation transport simulations, and correct for molecular hydrogen with empirical or theoretical relations. The `standard-resolution' simulations systematically underestimate HI column densities, leading to an HI deficiency in low-mass ($M_\star < 10^{10}M_\odot$) galaxies and poor reproduction of the observed HI mass function. These shortcomings are largely absent from EAGLE simulations featuring a factor of 8 (2) better mass (spatial) resolution, within which the HI mass of galaxies evolves more mildly from $z=1$ to $0$ than in the standard-resolution simulations. The largest-volume simulation reproduces the observed clustering of HI systems, and its dependence on HI-richness. At fixed $M_\star$, galaxies acquire more HI in simulations with stronger feedback, as they become associated with more massive haloes and higher infall rates. They acquire less HI in simulations with a greater star formation efficiency, since the star formation and feedback necessary to balance the infall rate is produced by smaller gas reservoirs. The simulations indicate that the HI of present-day galaxies was acquired primarily by the smooth accretion of ionized, intergalactic gas at $z\simeq1$, which later self-shields, and that only a small fraction is contributed by the reincorporation of gas previously heated strongly by feedback. HI reservoirs are highly dynamic: over $40$ percent of HI associated with $z=0.1$ galaxies is converted to stars or ejected by $z=0$.

Asymmetries in the neutral hydrogen gas distribution and kinematics of galaxies are thought to be indicators for both gas accretion and gas removal processes. These are of fundamental importance for galaxy formation and evolution. Upcoming large blind HI surveys will provide tens of thousands of galaxies for a study of these asymmetries in a proper statistical way. Due to the large number of expected sources and the limited resolution of the majority of objects, detailed modelling is not feasible for most detections. We need fast, automatic and sensitive methods to classify these objects in an objective way. Existing non-parametric methods suffer from effects like the dependence on signal to noise, resolution and inclination. Here we show how to correctly take these effects into account and show ways to estimate the precision of the methods. We will use existing and modelled data to give an outlook on the performance expected for galaxies observed in the various sky surveys planned for e.g. WSRT/APERTIF and ASKAP.

We use the EAGLE suite of cosmological hydrodynamical simulations to study how the HI content of present-day galaxies depends on their environment. We show that EAGLE reproduces observed HI mass-environment trends very well, while semi-analytic models typically overpredict the average HI masses in dense environments. The environmental processes act primarily as an on/off switch for the HI content of satellites with stellar mass Mstar>10^9 Msun. At a fixed Mstar, the fraction of HI-depleted satellites increases with increasing host halo mass M200 in response to stronger environmental effects, while at a fixed M200 it decreases with increasing satellite Mstar as the gas is confined by deeper gravitational potentials. HI-depleted satellites reside mostly, but not exclusively, within the virial radius r200 of their host halo. We investigate the origin of these trends by focussing on three environmental mechanisms: ram pressure stripping by the intra-group medium, tidal stripping by the host halo, and satellite-satellite encounters. By tracking back in time the evolution of the HI-depleted satellites, we find that the most common cause of HI removal is satellite encounters. The timescale for HI removal is typically less than 0.5 Gyr. Tidal stripping occurs in halos of M200<10^{14} Msun within half r200, while the other processes act also in more massive halos, generally within r200. Conversely, we find that ram pressure stripping is the most common mechanism that disturbs the HI morphology of galaxies at redshift z=0. This implies that HI removal due to satellite-satellite interactions occurs on shorter timescales than the other processes.

We revisit the HI size-mass (D$_{\rm HI}$-M$_{\rm HI}$) relation of galaxies with a sample of more than 500 nearby galaxies covering over five orders of magnitude in HI mass and more than ten $B$-band magnitudes. The relation is remarkably tight with a scatter $\sigma \sim$0.06 dex, or 14%. The scatter does not change as a function of galaxy luminosity, HI richness or morphological type. The relation is linked to the fact that dwarf and spiral galaxies have a homogenous radial profile of HI surface density in the outer regions when the radius is normalised by D$_{\rm HI}$. The early-type disk galaxies typically have shallower HI radial profiles, indicating a different gas accretion history. We argue that the process of atomic-to-molecular gas conversion or star formation cannot explain the tightness of the D$_{\rm HI}$-M$_{\rm HI}$ relation. This simple relation puts strong constraints on simulation models for galaxy formation.

Using data taken as part of the Bluedisk project we study the connection between neutral hydrogen (HI) in the environment of spiral galaxies and that in the galaxies themselves. We measure the total HI mass present in the environment in a statistical way by studying the distribution of noise peaks in the HI data cubes obtained for 40 galaxies observed with WSRT. We find that galaxies whose HI mass fraction is high relative to standard scaling relations have an excess HI mass in the surrounding environment as well. Gas in the environment consists of gas clumps which are individually below the detection limit of our HI data. These clumps may be hosted by small satellite galaxies and\or be the high-density peaks of a more diffuse gas distribution in the inter-galactic medium. We interpret this result as an indication for a picture in which the HI-rich central galaxies accrete gas from an extended gas reservoir present in their environment.

We present the study of a set of N-body+SPH simulations of a Milky Way-like system produced by the radiative cooling of hot gas embedded in a dark matter halo. The galaxy and its gaseous halo evolve for 10 Gyr in isolation, which allows us to study how internal processes affect the evolution of the system. We show how the morphology, the kinematics and the evolution of the galaxy are affected by the input supernova feedback energy E$_{\rm SN}$, and we compare its properties with those of the Milky Way. Different values of E$_{\rm SN}$ do not significantly affect the star formation history of the system, but the disc of cold gas gets thicker and more turbulent as feedback increases. Our main result is that, for the highest value of E$_{\rm SN}$ considered, the galaxy shows a prominent layer of extra-planar cold (log(T)<4.3) gas extended up to a few kpc above the disc at column densities of $10^{19}$ cm$^{-2}$. The kinematics of this material is in agreement with that inferred for the HI halos of our Galaxy and NGC 891, although its mass is lower. Also, the location, the kinematics and the typical column densities of the hot (5.3<log(T)<5.7) gas are in good agreement with those determined from the O$_{\rm VI}$ absorption systems in the halo of the Milky Way and external galaxies. In contrast with the observations, however, gas at log(T)<5.3 is lacking in the circumgalactic region of our systems.

As part of the Bluedisk survey we analyse the radial gas-phase metallicity profiles of 50 late-type galaxies We compare the metallicity profiles of a sample of HI-rich galaxies against a control sample of HI-'normal' galaxies. We find the metallicity gradient of a galaxy to be strongly correlated with its HI mass fraction (M(HI) / Mstar). We note that some galaxies exhibit a steeper metallicity profile in the outer disc than in the inner disc. These galaxies are found in both the HI-rich and control samples. This contradicts a previous indication that these outer drops are exclusive to HI-rich galaxies. These effects are not driven by bars, although we do find some indication that barred galaxies have flatter metallicity profiles. By applying a simple analytical model we are able to account for the variety of metallicity profiles that the two samples present. The success of this model implies that the metallicity in these isolated galaxies may be in a local equilibrium, regulated by star formation. This insight could provide an explanation of the observed local mass-metallicity relation.

We introduce SoFiA, a flexible software application for the detection and parameterization of sources in 3D spectral-line datasets. SoFiA combines for the first time in a single piece of software a set of new source-finding and parameterization algorithms developed on the way to future HI surveys with ASKAP (WALLABY, DINGO) and APERTIF. It is designed to enable the general use of these new algorithms by the community on a broad range of datasets. The key advantages of SoFiA are the ability to: search for line emission on multiple scales to detect 3D sources in a complete and reliable way, taking into account noise level variations and the presence of artefacts in a data cube; estimate the reliability of individual detections; look for signal in arbitrarily large data cubes using a catalogue of 3D coordinates as a prior; provide a wide range of source parameters and output products which facilitate further analysis by the user. We highlight the modularity of SoFiA, which makes it a flexible package allowing users to select and apply only the algorithms useful for their data and science questions. This modularity makes it also possible to easily expand SoFiA in order to include additional methods as they become available. The full SoFiA distribution, including a dedicated graphical user interface, is publicly available for download.

We analyze the radial distribution of HI gas for 23 disk galaxies with unusually high HI content from the Bluedisk sample, along with a similar-sized sample of "normal" galaxies. We propose an empirical model to fit the radial profile of the HI surface density, an exponential function with a depression near the center. The radial HI surface density profiles are very homogeneous in the outer regions of the galaxy; the exponentially declining part of the profile has a scale-length of $\sim 0.18$ R1, where R1 is the radius where the column density of the HI is 1 M$_{\odot}$ pc$^{-2}$. This holds for all galaxies, independent of their stellar or HI mass. The homogenous outer profiles, combined with the limited range in HI surface density in the non-exponential inner disk, results in the well-known tight relation between HI size and HI mass. By comparing the radial profiles of the HI-rich galaxies with those of the control systems, we deduce that in about half the galaxies, most of the excess gas lies outside the stellar disk, in the exponentially declining outer regions of the HI disk. In the other half, the excess is more centrally peaked. We compare our results with existing smoothed-particle hydrodynamical simulations and semi-analytic models of disk galaxy formation in a $\Lambda$ Cold Dark Matter universe. Both the hydro simulations and the semi-analytic models reproduce the HI surface density profiles and the HI size-mass relation without further tuning of the simulation and model inputs. In the semi-analytic models, the universal shape of the outer HI radial profiles is a consequence of the {\em assumption} that infalling gas is always distributed exponentially. The conversion of atomic gas to molecular form explains the limited range of HI surface densities in the inner disk. These two factors produce the tight HI mass-size relation.

We measure the neutral atomic hydrogen (HI) gas content of field galaxies at intermediate redshifts of z ~ 0.1 and z ~ 0.2 using hydrogen 21-cm emission lines observed with the Westerbork Synthesis Radio Telescope (WSRT). In order to make high signal-to-noise ratio detections, an HI signal stacking technique is applied: HI emission spectra from multiple galaxies, optically selected by the CNOC2 redshift survey project, are co-added to measure the average HI mass of galaxies in the two redshift bins. We calculate the cosmic HI gas densities ({\Omega}_{HI}) at the two redshift regimes and compare those with measurements at other redshifts to investigate the global evolution of the HI gas density over cosmic time. From a total of 59 galaxies at z ~ 0.1 we find {\Omega}_{HI} = (0.33 $\pm$ 0.05) ~ $\times$ 10$^{-3}$, and at z ~ 0.2 we find {\Omega}_{HI} = (0.34 $\pm$ 0.09) ~ $\times$ 10$^{-3}$, based on 96 galaxies. These measurements help bridge the gap between high-z damped Lyman-$\alpha$ observations and blind 21-cm surveys at $z=$ 0. We find that our measurements of {\Omega}_{HI} at z ~ 0.1 and 0.2 are consistent with the HI gas density at z ~ 0 and that all measurements of {\Omega}_{HI} from 21-cm emission observations at $z \la$ ~ 0.2 are in agreement with no evolution of the HI gas content in galaxies during the last 2.4 Gyr.

We introduce the "Bluedisk" project, a large program at the Westerbork Synthesis Radio Telescope (WSRT) that has mapped the HI in a sample of 23 nearby galaxies with unusually high HI mass fractions, along with a similar-sized sample of control galaxies matched in stellar mass, size, inclination and redshift. This paper presents the sample selection, observational set-up, data reduction strategy, and a first analysis of the sizes and structural properties of the HI disks. We find that the HI-rich galaxies lie on the same HI mass versus HI size relation as normal spiral galaxies, extending it to total HI masses of $2 \times 10^{10} M_{\odot}$ and radii R1 of $\sim 100$ kpc (where R1 is defined as the radius where the HI column density reaches 1 $M_{\odot}$ pc$^{-2}$). HI-rich galaxies have significantly larger values of HI-to-optical size ratio at fixed stellar mass, concentration index, stellar and star formation rate surface density compared to the control sample. The disks of HI-rich galaxies are also significantly more clumpy (i.e. have higher HI Gini and $\Delta$Area coefficient) than those of normal spirals. There is no evidence that the disks of HI-rich galaxies are more disturbed: HI-rich galaxies exhibit no difference with respect to control samples in their distributions of HI asymmetry indices or optical/HI disk position angle differences. In fact, the center of the HI distribution corresponds more closely with the center of the optical light in the HI-rich galaxies than in the controls. All these results argue against a scenario in which new gas has been brought in by mergers. It is possible that they may be more consistent with cooling from a surrounding quasi-static halo of warm/hot gas.

Our aim is to assess the benefits and limitations of using the redundant visibility information in regular phased array systems for improving the calibration. Regular arrays offer the possibility to use redundant visibility information to constrain the calibration of the array independent of a sky model and a beam models of the station elements. It requires a regular arrangement in the configuration of array elements and identical beam patterns. We revised a calibration method for phased array stations using the redundant visibility information in the system and applied it successfully to a LOFAR station. The performance and limitations of the method were demonstrated by comparing its use on real and simulated data. The main limitation is the mutual coupling between the station elements, which leads to non-identical beams and stronger baseline dependent noise. Comparing the variance of the estimated complex gains with the Cramer-Rao Bound (CRB) indicates that redundancy is a stable and optimum method for calibrating the complex gains of the system. Our study shows that the use of the redundant visibility does improve the quality of the calibration in phased array systems. In addition it provides a powerful tool for system diagnostics. Our results demonstrate that designing redundancy in both the station layout and the array configuration of future aperture arrays is strongly recommended. In particular in the case of the Square Kilometre Array with its dynamic range requirement which surpasses any existing array by an order of magnitude.

We review and discuss aspects of Cosmic Voids that form the background for our Void Galaxy Survey (see accompanying paper by Stanonik et al.). Following a sketch of the general characteristics of void formation and evolution, we describe the influence of the environment on their development and structure and the characteristic hierarchical buildup of the cosmic void population. In order to be able to study the resulting tenuous void substructure and the galaxies populating the interior of voids, we subsequently set out to describe our parameter free tessellation-based watershed void finding technique. It allows us to trace the outline, shape and size of voids in galaxy redshift surveys. The application of this technique enables us to find galaxies in the deepest troughs of the cosmic galaxy distribution, and has formed the basis of our void galaxy program.

Evidence for the accretion of material in spiral galaxies has grown over the past years and clear signatures can be found in HI observations of galaxies. We describe here new detailed and sensitive HI synthesis observations of a few nearby galaxies (NGC 3359, NGC 4565 and NGC 6946) which show that indeed accretion of small amounts of gas is taking place. These should be regarded as examples illustrating a general phenomenon of gas infall in galaxies. Such accretion may also be at the origin of the gaseous halos which are being found around spirals. Probably it is the same kind of phenomenon of material infall as observed in the stellar streams in the halo and outer parts of our galaxy and M 31

Observations of the nearby starburst galaxy NGC 253 in the 21-cm line reveal the presence of neutral hydrogen in the halo, up to 12 kpc from the galactic plane. This extra-planar HI is found only in one half of the galaxy and is concentrated in a half-ring structure and plumes which are lagging in rotation with respect to the disk. The HI plumes are seen bordering the bright Halpha and X-ray halo emission. It is likely that, as proposed earlier for the Halpha and the X-rays, also the origin of the extra-planar HI is related to the central starburst and to the active star formation in the disk. A minor merger and gas accretion are also discussed as possible explanations. The HI disk is less extended than the stellar disk. This may be the result of ionization of its outer parts or, alternatively, of tidal or ram pressure stripping.

Multi-wavelength observations of nearby spiral galaxies have shown that neutral and ionized gas are present up to a few kpc from the disk and that star formation and supernovae probably play an important role in bringing gas into the halo. We have obtained very sensitive HI observations of the face-on galaxy NGC 6946 and of the nearly edge-on starburst galaxy NGC 253. We find high velocity HI clouds in NGC 6946 and extra-planar gas with anomalous velocities in NGC 253. In both galaxies there seems to be a close connection between the star-forming disk and the halo HI. In the outer parts of NGC 6946 there is also evidence for recent gas accretion.

Studies of the HI in galaxies have clearly shown that subtle details of the HI distribution and kinematics often harbour key information for understanding the structure and evolution of galaxies. Evidence for the accretion of material has grown over the past many years and clear signatures can be found in HI observations of galaxies. We have obtained new detailed and sensitive HI synthesis observations of three nearby galaxies which are suspected of capturing small amounts of HI and show that indeed accretion of small amounts of gas is taking place in these galaxies. This could be the same kind of phenomenon of material infall as observed in the stellar streams in the halo and outer parts of our galaxy and M31

We are studying the properties of the holes and the high velocity gas in NGC 6946. Here we present some puzzling results.

We present a deep multi-colour CCD mosaic of the Coma cluster (Abell 1656), covering 5.2 deg^2 in the B and r bands, and 1.3 deg^2 in the U band. This large, homogeneous data set provides a valuable low redshift comparison sample for studies of galaxies in distant clusters. In this paper we present our survey, and study the dependence of the galaxy luminosity function (LF) on passband and radial distance from the cluster centre. The U, B and r band LFs of the complete sample cannot be represented by single Schechter functions. For the central area, r<245 h_{100}^{-1} kpc, we find best-fitting Schechter parameters of M^{*}_U=-18.60^{+0.13}_{-0.18} and \alpha_U=-1.32^{+0.018}_{-0.028}, M^{*}_B=-19.79^{+0.18}_{-0.17} and \alpha_B=-1.37^{+0.024}_{-0.016} and M^{*}_r=-20.87^{+0.12}_{-0.17} and \alpha_r=-1.16^{+0.012}_{-0.019}. The LF becomes steeper at larger radial distance from the cluster centre. The effect is most pronounced in the U band. This result is consistent with the presence of a star forming dwarf population at large distance from the cluster centre, which may be in the process of being accreted by the cluster. The shapes of the LFs of the NGC 4839 group support a scenario in which the group has already passed through the centre.