Using a spectral stacking technique, we measure the neutral hydrogen (HI) properties of a sample of galaxies at $z < 0.11$ across 35 pointings of the Westerbork Synthesis Radio Telescope (WSRT). The radio data contains 1,895 galaxies with redshifts and positions known from the Sloan Digital Sky Survey (SDSS). We carefully quantified the effects of sample bias, aperture used to extract spectra, sidelobes and weighting technique and use our data to provide a new estimate for the cosmic HI mass density. We find a cosmic HI mass density of $\Omega_{\rm HI} = (4.02 \pm 0.26)\times 10^{-4} h_{70}^{-1}$ at $\langle z\rangle = 0.066$, consistent with measurements from blind HI surveys and other HI stacking experiments at low redshifts. The combination of the small interferometer beam size and the large survey volume makes our result highly robust against systematic effects due to confusion at small scales and cosmic variance at large scales. Splitting into three sub-samples with $\langle z\rangle$ = 0.038, 0.067 and 0.093 shows no significant evolution of the HI gas content at low redshift.

We present the results of a multiwavelength survey of HI-excess galaxies, an intriguing population with large HI reservoirs associated with little current star formation. These galaxies have stellar masses $M_{\star} >10^{10}$ M$_{\odot}$, and were identified as outliers in the gas fraction vs. NUV$-r$ color and stellar mass surface density scaling relations based on the GALEX Arecibo SDSS Survey (GASS). We obtained HI interferometry with the GMRT, Keck optical long-slit spectroscopy and deep optical imaging (where available) for four galaxies. Our analysis reveals multiple possible reasons for the HI excess in these systems. One galaxy, AGC 10111, shows an HI disk that is counter-rotating with respect to the stellar bulge, a clear indication of external origin of the gas. Another galaxy appears to host a Malin 1-type disk, where a large specific angular momentum has to be invoked to explain the extreme $M_{\rm HI}$/$M_{\star}$ ratio of 166$\%$. The other two galaxies have early-type morphology with very high gas fractions. The lack of merger signatures (unsettled gas, stellar shells and streams) in these systems suggests that these gas-rich disks have been built several Gyr-s ago, but it remains unclear how the gas reservoirs were assembled. Numerical simulations of large cosmological volumes are needed to gain insight into the formation of these rare and interesting systems.

We present our multiwavelength analysis of a prototype \HI-excess galaxy, GASS 3505, selected based on having a large gas content ($M_{\rm HI} = 10^{9.9}$ \msun) compared to its little associated star formation activity ($\sim$0.1 \msun\ yr$^{-1}$) in the GALEX Arecibo SDSS Survey (GASS). Very Large Array (VLA) observations show that the \HI\ in GASS 3505 is distributed in a regularly rotating, extended ($\sim$50 kpc radius) gas ring. In the SDSS optical image GASS 3505 appears as a bulge-dominated galaxy, however deep optical imaging reveals low surface brightness ($\gtrsim25$ mag arcsec$^{-2}$) stellar emission around the central bulge. Direct evidence for accretion is detected in form of an extended ($\sim$60 kpc) stellar stream, showing that GASS 3505 has experienced a minor merger in the recent past. We investigate the possibility that the \HI\ ring in GASS 3505 was accreted in such a merger event using N-body and smoothed particle hydrodynamic (SPH) simulations. The best model that reproduces the general properties (i.e., gas distribution and kinematics, stellar morphology) of the galaxy involves a merger between the central bulge and a gas-rich ($M_{\star}$ = 10$^9$ \msun\ and $M_{\rm HI}$/$M_{\star}$ = 10) disk galaxy. However, small discrepancies in the observed and modeled properties could suggest that other sources of gas have to be involved in the build-up of the gas reservoir. This work is the first step toward a larger program to investigate the physical mechanisms that drive the large scatter in the gas scaling relations of nearby galaxies.

Newly born and young radio sources are in a delicate phase of their life. Their jets are fighting their way through the surrounding gaseous medium, strongly experiencing this interaction while, at the same time, impacting and affecting the interstellar medium (ISM). Here we present the results from two studies of HI (in absorption) and molecular gas illustrating what can be learned from these phases of the gas. We first describe a statistical study with the WSRT. The study shows that the young radio sources not only have an higher detection rate of HI, but also systematically broader and more asymmetric HI profiles, most of them blueshifted. This supports the idea that we are looking at young radio jets making their way through the surrounding ISM, which also appears to be, on average, richer in gas than in evolved radio sources. Signatures of the impact of the jet are seen in the kinematics of the gas. However, even among the young sources, we identify a population that remains undetected in HI even after stacking their profiles. Orientation effects can only partly explain the result. These objects either are genuinely gas-poor or have different conditions of the medium, e.g. higher spin temperature. We further present the ALMA study of molecular gas in IC5063 to trace in detail the jet impacting the ISM. The kinematics of the cold, molecular gas co-spatial with the radio plasma shows this process in action. The ALMA data reveal a fast outflow of molecular gas extending along the entire radio jet (~1 kpc), with the highest outflow velocities at the location of the brighter hot-spot. We propose a scenario where the radio jet is expanding into a clumpy medium, interacting directly with the clouds and inflating a cocoon that drives a lateral outflow into the ISM.

We study the atomic neutral hydrogen (HI) content of $\sim$1600 galaxies up to $z \sim 0.1$ using stacking techniques. The observations were carried out with the Westerbork Synthesis Radio Telescope (WSRT) in the area of the SDSS South Galactic Cap (SSGC), where we selected a galaxy sample from the SDSS spectroscopic catalog. Multi-wavelength information is provided by SDSS, NVSS, GALEX, and WISE. We use the collected information to study HI trends with color, star-forming, and active galactic nuclei (AGN) properties. Using NUV-r colors, galaxies are divided into blue cloud, green valley and red sequence galaxies. As expected based on previous observations, we detect HI in green valley objects with lower amounts of HI than blue galaxies, while stacking only produces a 3-$\sigma$ upper limit for red galaxies with M$_{\rm HI}$ $<$ (5 $\pm$ 1.5) $\times$ 10$^{8}$ M$_{\odot}$ and M$_{\rm HI}/\rm{L}_r$ $<$ 0.02 $\pm$ 0.006 $\rm M_{\odot} / \rm L_{\odot} $. We find that the HI content is more dependent on NUV-r color, and less on ionization properties, in the sense that regardless of the presence of an optical AGN (based on optical ionization line diagnostics), green-valley galaxies always show HI, whereas red galaxies only produce an upper limit. This suggests that feedback from optical AGN is not the (main) reason for depleting large-scale gas reservoirs. Low-level radio continuum emission in our galaxies can stem either from star formation, or from AGN. We use the WISE color-color plot to separate these phenomena by dividing the sample into IR late-type and IR early-type galaxies. We find that the radio emission in IR late-type galaxies stems from enhanced star formation, and this group is detected in HI. However, IR early-type galaxies lack any sign of HI gas and star formation activity, suggesting that radio AGN are likely to be the source of radio emission in this population.

We present an analysis of the HI absorption in a sample of 101 flux-selected radio AGN (S_1.4 GHz > 50 mJy) observed with the Westerbork Synthesis Radio Telescope (WSRT). HI absorption is detected in 32 galaxies, showing a broad variety of widths, shapes and kinematical properties. We characterize the HI spectra of the individual detections using the busy function (Westmeier et al. 2014). With the goal of identifying different morphological structures of HI, we study the kinematical and radio source properties of the detections as function of their width. Narrow lines (FWHM < 100 km/s) mostly lie at the systemic velocity and are likely produced by regularly rotating HI disks or clouds. More HI disks can be present among galaxies with lines of intermediate widths (100 km/s < FWHM < 200 km/s), however the HI in these sources is more unsettled. We study the asymmetry parameter and blueshift/redshift distribution of the lines. We argue that if a profile is broad, it is also asymmetric and shifted relative to the systemic velocity because it is tracing unsettled HI gas. Besides three of the broadest (up to FW20 = 825 km/s) detections, which are associated with gas-rich mergers, we find three new cases of blueshifted broad wings (with FW20 >= 500 km/s). These detections are good candidates for being HI outflows. The detection rate of HI outflows is 5 percent in the total radio AGN sample. This fraction represents a lower limit, however it could suggests that, if outflows are a characteristic phenomenon of all radio sources, they would have a short depletion timescale compared to the lifetime of the AGN. Blueshifted and broad/asymmetric lines are more often present among young, compact sources, suggesting that in these AGN the HI is likely to be unsettled. This may arise due to jet-cloud interactions, as young radio sources clear their way through the rich ambient gaseous medium.

Using the WSRT, we carried out shallow HI absorption observations of a flux-selected (S > 50 mJy) sample of 93 radio AGN with available SDSS redshifts between 0.02 < z < 0.23. We study the gas properties of radio AGN down to fluxes not systematically explored before using, for the first time, stacking of extragalactic HI absorption. Despite the shallow observations, we obtained a direct detection rate of ~29%, comparable with deeper studies. Detections are found at every flux level, showing that HI absorption detections are not biased toward brighter sources. The stacks of detections and non-detections reveal a clear dichotomy in the presence of HI, with the 27 detections showing an average peak {\tau} = 0.02, while the 66 non-detections remain undetected with an upper limit {\tau} < 0.002. Separating the sample into compact and extended AGN increases the detection rate, {\tau}, and FWHM for the compact sample. The dichotomy for the stacked profiles of detections and non-detections still holds between these two groups. We argue that orientation effects connected to a disk-like distribution of HI can be partly responsible for the dichotomy, however some of our galaxies must be genuinely depleted of cold gas. A fraction of the compact sources are confirmed by previous studies as likely young radio sources. These show an even higher detection rate of 55%. Along with their high integrated optical depth and wider profile, this reinforces the idea that young radio AGN are particularly rich in atomic gas. Part of our motivation is to probe for the presence of HI outflows. However, the stacked profiles do not reveal any significant blueshifted wing. Our results are particularly relevant for future surveys. The lack of bias toward bright sources is encouraging for the search of HI at lower radio fluxes. The results also represent a reference point for search for HI absorption at higher redshifts.

The point source list of the Wilkinson Microwave Anisotropy Probe (WMAP) is a uniform, all-sky catalogue of bright sources with flux density measurements at high (up to 94 GHz) radio frequencies. We investigated the five-year WMAP list to compile a new catalogue of bright and compact extragalactic radio sources to be potentially studied with Very Long Baseline Interferometry at millimeter wavelengths (mm-VLBI) and Space VLBI (SVLBI). After comparing the WMAP data with the existing mm-VLBI catalogues, we sorted out the yet unexplored sources. Using the 41, 61 and 94 GHz WMAP flux densities, we calculated the spectral indices. By collecting optical identifications, lower-frequency radio flux densities and VLBI images from the literature, we created a list of objects which have not been investigated with VLBI at 86 GHz before. With total flux density at least 1 Jy and declination above -40 degree, we found 37 suitable new targets. It is a nearly 25% addition to the known mm-VLBI sources. Such objects are also potentially useful as phase-reference calibrators for the future Japanese SVLBI mission ASTRO-G at its highest observing frequency (43 GHz). The phase-referencing capability of ASTRO-G would allow long integrations and hence better sensitivity for observing faint target sources close to suitable phase calibrators in the sky.