A recent analysis of cosmic ray air showers observed at the Pierre Auger Observatory indicates that nearby starburst galaxies (SBGs) might be the cause of ~10% of the Ultra-High-Energy Cosmic Ray flux at energies E > 39 EeV. Since high energy neutrinos are a direct product of cosmic ray interactions, we investigate SBGs as a possible source of some of the ~ 10^-2 - 1 PeV neutrinos observed at IceCube. A statistical analysis is performed to establish the degree of positional correlation between the observed neutrinos and a set of 45 nearby radio- and infrared-bright SBGs. Our results are consistent with no causal correlation. However, a scenario where ~10% of the High Energy Starting Events (HESE) in the detector are coming from the candidate SBGs is not excluded. The same conclusion is reached for different data subsets, as well as two different subsets of SBGs motivated by the Pierre Auger Observatory analysis.

(abridged) Radio recombination lines (RRLs) at frequencies $\nu$ < 250 MHz trace the cold, diffuse phase of the ISM. Next generation low frequency interferometers, such as LOFAR, MWA and the future SKA, with unprecedented sensitivity, resolution, and large fractional bandwidths, are enabling the exploration of the extragalactic RRL universe. We observed the radio quasar 3C 190 (z~1.2) with the LOFAR HBA. In reducing this data for spectroscopic analysis, we have placed special emphasis on bandpass calibration. We devised cross-correlation techniques to significantly identify the presence of RRLs in a low frequency spectrum. We demonstrate the utility of this method by applying it to existing low-frequency spectra of Cassiopeia A and M 82, and to the new observations of 3C 190. RRLs have been detected in the foreground of 3C 190 at z = 1.12355 (assuming a carbon origin), owing to the first detection of RRLs outside of the local universe (first reported in Emig et al. 2019). Towards the Galactic supernova remnant Cas A, we uncover three new detections: (1) C$\epsilon$-transitions ($\Delta$n = 5) for the first time at low radio frequencies, (2) H$\alpha$-transitions at 64 MHz with a FWHM of 3.1 km/s, the most narrow and one of the lowest frequency detections of hydrogen to date, and (3) C$\alpha$ at v$_{LSR}$ = 0 km/s in the frequency range 55-78 MHz for the first time. Additionally we recover C$\alpha$, C$\beta$, C$\gamma$, and C$\delta$ from the -47 km/s and -38 km/s components. In the nearby starburst galaxy, M 82, we do not find a significant feature. Our current searches for RRLs in LOFAR observations are limited to narrow (< 100 km/s) features, owing to the relatively small number of channels available for continuum estimation. Future strategies making use of larger contiguous frequency coverage would aid calibration to deeper sensitivities and broader features.

Context. The combined use of carbon radio recombination lines (CRRLs) and the 158 $\mu$m-[CII] line is a powerful tool for the study of the energetics and physical conditions (e.g., temperature and density) of photodissociation regions (PDRs). However, there are few observational studies that exploit this synergy. Aims. Here we explore the relation between CRRLs and the 158 $\mu$m-[CII] line in light of new observations and models. Methods. We present new and existing observations of CRRLs in the frequency range 0.15--230 GHz with ALMA, VLA, the GBT, Effelsberg 100m, and LOFAR towards Orion~A (M42). We complement these observations with SOFIA observations of the 158 $\mu$m-[CII] line. We studied two PDRs: the foreground atomic gas, known as the Veil, and the dense PDR between the HII region and the background molecular cloud. Results. In the Veil we are able to determine the gas temperature and electron density, which we use to measure the ionization parameter and the photoelectric heating efficiency. In the dense PDR, we are able to identify a layered PDR structure at the surface of the molecular cloud to the south of the Trapezium cluster. There we find that the radio lines trace the colder portion of the ionized carbon layer, the C$^{+}$/C/CO interface. By modeling the emission of the $158$~$\mu$m-[CII] line and CRRLs as arising from a PDR we derive a thermal pressure $>5\times10^{7}$ K cm$^{-3}$ and a radiation field $G_{0}\approx10^{5}$ close to the Trapezium. Conclusions. This work provides additional observational support for the use of CRRLs and the 158 $\mu$m-[CII] line as complementary tools to study dense and diffuse PDRs, and highlights the usefulness of CRRLs as probes of the C$^{+}$/C/CO interface.

Recombination lines involving high principal quantum numbers populate the radio spectrum in large numbers. Low-frequency (<1 GHz) observations of radio recombination lines (RRLs) primarily from carbon and hydrogen offer a new, if not unique, way to probe cold, largely atomic gas and warm, ionised gas in other galaxies. Furthermore, RRLs can be used to determine the physical state of the emitting regions, such as temperature and density. These properties make RRLs, potentially, a powerful tool of extragalactic ISM physics. At low radio frequencies, its conceivable to detect RRLs out to cosmological distances when illuminated by a strong radio continuum. However, they are extremely faint (tau ~ 1e-3 -- 1e-4) and have so far eluded detection outside of the local universe. With LOFAR observations of the radio quasar 3C 190 (z=1.1946), we aim to demonstrate that the ISM can be explored out to great distances through low-frequency RRLs. We report the detection of RRLs in the frequency range 112--163 MHz in the spectrum of 3C 190. Stacking 13 a-transitions with principal quantum numbers n=266-301, a peak 6sigma feature of optical depth, tau(peak) = (1.0 +\- 0.2) x 1e-3 and FWHM = 31.2 +/- 8.3 km/s was found at z=1.124. This corresponds to a velocity offset of -9965 km/s with respect to the systemic redshift of 3C 190. We consider three interpretations of the origin of the RRL emission: an intervening dwarf-like galaxy, an AGN-driven outflow, and the inter-galactic medium. We argue that the RRLs most likely originate in a dwarf-like galaxy (M ~ 1e9 Msun) along the line of sight, although we cannot rule out an AGN-driven outflow. We do find the RRLs to be inconsistent with an inter-galactic medium origin. With this detection, we have opened up a new way to study the physical properties of cool, diffuse gas out to cosmological distances.

Context: New generation low-frequency telescopes are exploring a new parameter space in terms of depth and resolution. The data taken with these interferometers, for example with the LOw Frequency ARray (LOFAR), are often calibrated in a low signal-to-noise ratio regime and the removal of critical systematic effects is challenging. The process requires an understanding of their origin and properties. Aim: In this paper we describe the major systematic effects inherent to next generation low-frequency telescopes, such as LOFAR. With this knowledge, we introduce a data processing pipeline that is able to isolate and correct these systematic effects. The pipeline will be used to calibrate calibrator observations as the first step of a full data reduction process. Methods: We processed two LOFAR observations of the calibrator 3C196: the first using the Low Band Antenna (LBA) system at 42-66 MHz and the second using the High Band Antenna (HBA) system at 115-189 MHz. Results: We were able to isolate and correct for the effects of clock drift, polarisation misalignment, ionospheric delay, Faraday rotation, ionospheric scintillation, beam shape, and bandpass. The designed calibration strategy produced the deepest image to date at 54 MHz. The image has been used to confirm that the spectral energy distribution of the average radio source population tends to flatten at low frequencies. Conclusions: We prove that LOFAR systematic effects can be described by a relatively small number of parameters. Furthermore, the identification of these parameters is fundamental to reducing the degrees of freedom when the calibration is carried out on fields that are not dominated by a strong calibrator.

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,\mu$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS.

Quantitative understanding of the interstellar medium requires knowledge of its physical conditions. Low frequency carbon radio recombination lines (CRRLs) trace cold interstellar gas, and can be used to determine its physical conditions (e.g., electron temperature and density). In this work we present spatially resolved observations of the low frequency ($\leq390$ MHz) CRRLs centered around C$268\alpha$, C$357\alpha$, C$494\alpha$ and C$539\alpha$ towards Cassiopeia A on scales of $\leq1.2$ pc. We compare the spatial distribution of CRRLs with other ISM tracers. This comparison reveals a spatial offset between the peak of the CRRLs and other tracers, which is very characteristic for photodissociation regions and that we take as evidence for CRRLs being preferentially detected from the surfaces of molecular clouds. Using the CRRLs we constrain the gas electron temperature and density. These constraints on the gas conditions suggest variations of less than a factor of two in pressure over $\sim1$ pc scales, and an average hydrogen density of $200$-$470$ cm$^{-3}$. From the electron temperature and density maps we also constrain the ionized carbon emission measure, column density and path length. Based on these, the hydrogen column density is larger than $10^{22}$ cm$^{-2}$, with a peak of $\sim4\times10^{22}$ cm$^{-2}$ towards the South of Cassiopeia A. Towards the southern peak the line of sight length is $\sim40$ pc over a $\sim2$ pc wide structure, which implies that the gas is a thin surface layer on a large (molecular) cloud that is only partially intersected by Cassiopeia A. These observations highlight the utility of CRRLs as tracers of low density extended HI and CO-dark gas halo's around molecular clouds.

We present a study of carbon radio recombination lines towards Cassiopeia A using LOFAR observations in the frequency range 10-33 MHz. Individual carbon $\alpha$ lines are detected in absorption against the continuum at frequencies as low as 16 MHz. Stacking several C$\alpha$ lines we obtain detections in the 11-16 MHz range. These are the highest signal-to-noise measurements at these frequencies. The peak optical depth of the C$\alpha$ lines changes considerably over the 11-33 MHz range with the peak optical depth decreasing from 4$\times10^{-3}$ at 33 MHz to 2$\times10^{-3}$ at 11 MHz, while the line width increases from 20 km s$^{-1}$ to 150 km s$^{-1}$. The combined change in peak optical depth and line width results in a roughly constant integrated optical depth. We interpret this as carbon atoms close to local thermodynamic equilibrium. In this work we focus on how the 11-33 MHz carbon radio recombination lines can be used to determine the gas physical conditions. We find that the ratio of the carbon radio recombination lines to that of the 158 $\mu$m [CII] fine-structure line is a good thermometer, while the ratio between low frequency carbon radio recombination lines provides a good barometer. By combining the temperature and pressure constraints with those derived from the line width we are able to constrain the gas properties (electron temperature and density) and radiation field intensity. Given the 1$\sigma$ uncertainties in our measurements these are; $T_{e}\approx68$-$98$ K, $n_{e}\approx0.02$-$0.035$ cm$^{-3}$ and $T_{r,100}\approx1500$-$1650$ K. Despite challenging RFI and ionospheric conditions, our work demonstrates that observations of carbon radio recombination lines in the 10-33 MHz range can provide insight into the gas conditions.

We combine wide and deep galaxy number-count data from GAMA, COSMOS/G10, HST ERS, HST UVUDF and various near-, mid- and far- IR datasets from ESO, Spitzer and Herschel. The combined data range from the far-UV (0.15microns) to far-IR (500microns), and in all cases the contribution to the integrated galaxy light (IGL) of successively fainter galaxies converges. Using a simple spline fit, we derive the IGL and the extrapolated-IGL in all bands. We argue undetected low surface brightness galaxies and intra-cluster/group light is modest, and that our extrapolated-IGL measurements are an accurate representation of the extra-galactic background light. Our data agree with most earlier IGL estimates and with direct measurements in the far-IR, but disagree strongly with direct estimates in the optical. Close agreement between our results and recent very high-energy experiments (H.E.S.S. and MAGIC), suggest that there may be an additional foreground affecting the direct estimates. The most likely culprit could be the adopted Zodiacal light model. Finally we use a modified version of the two-component model to integrate the EBL and obtain measurements of the Cosmic Optical Background (COB) and Cosmic Infrared Background (CIB) of: $24^{+4}_{-4}$nWm$^{-2}$sr$^{-1}$ and $26^{+5}_{-5}$nWm$^{-2}$sr$^{-1}$ respectively (48:52%). Over the next decade, upcoming space missions such as Euclid and WFIRST, have the capacity to reduce the COB error to $<1\%$, at which point comparisons to the very high energy data could have the potential to provide a direct detection and measurement of the reionisation field.

The IceCube Neutrino Observatory has provided the first map of the high energy (~ 0.01 -- 1 PeV) sky in neutrinos. Since neutrinos propagate undeflected, their arrival direction is an important identifier for sources of high energy particle acceleration. Reconstructed arrival directions are consistent with an extragalactic origin, with possibly a galactic component, of the neutrino flux. We present a statistical analysis of positional coincidences of the IceCube neutrinos with known astrophysical objects from several catalogs. When considering starburst galaxies with the highest flux in gamma-rays and infrared radiation, up to n=8 coincidences are found, representing an excess over the ~4 predicted for the randomized, or "null" distribution. The probability that this excess is realized in the null case, the p-value, is p=0.042. This value falls to p=0.003 for a partial subset of gamma-ray-detected starburst galaxies and superbubble regions in the galactic neighborhood. Therefore, it is possible that starburst galaxies, and the typically hundreds of superbubble regions within them, might account for a portion of IceCube neutrinos. The physical plausibility of such correlation is discussed briefly.