We present an estimate of the optical luminosity function (OLF) of LOFAR radio-selected quasars (RSQs) at $1.4<z<5.0$ in the $9.3\:\textrm{deg}^{2}$ NOAO Deep Wide-field survey (NDWFS) of the Bo\"otes field. The selection was based on optical/mid-ir photometry used to train three different machine learning (ML) algorithms. Objects taken as quasars by the ML algorithms are required to be detected at $5\sigma$ significance in deep radio maps to be classified as candidate quasars. The optical imaging came from the SDSS and the PS1 $3\pi$ survey; mid-ir photometry was taken from the SDWFS survey; and radio data was obtained from deep LOFAR imaging of the NDWFS-Bo\"otes field. The requirement of a $5\sigma$ LOFAR detection allowed us to reduce the stellar contamination in our sample by two orders of magnitude. The sample comprises 134 objects, including both photometrically selected candidate quasars (47) and spectroscopically confirmed quasars (83). The depth of our LOFAR observations allowed us to detect the radio-emission of quasars that would be otherwise classified as radio-quiet. Around $65\%$ of the quasars in the sample are fainter than $M_{\textrm{1450}}=-24.0$, a regime where the OLF of quasars selected through their radio emission, has not been investigated in detail. It has been demonstrated that in cases where mid-ir wedge-based AGN selection is not possible due to a lack of appropriate data, the selection of quasars using ML algorithms trained with optical/mid-ir photometry in combination with LOFAR data provides an excellent approach for obtaining samples of quasars. We demonstrate that RSQs show an evolution similar to that exhibited by faint quasars $(M_{\textrm{1450}}\leq-22.0)$. Finally, we find that RSQs may compose up to $\sim20\%$ of the whole faint quasar population (radio-detected plus radio-undetected).

We have conducted a deep survey (with a central rms of $55\mu\textrm{Jy}$) with the LOw Frequency ARray (LOFAR) at 120-168 MHz of the Bo\"otes field, with an angular resolution of $3.98^{''}\times6.45^{''}$, and obtained a sample of 10091 radio sources ($5\sigma$ limit) over an area of $20\:\textrm{deg}^{2}$. The astrometry and flux scale accuracy of our source catalog is investigated. The resolution bias, incompleteness and other systematic effects that could affect our source counts are discussed and accounted for. The derived 150 MHz source counts present a flattening below sub-mJy flux densities, that is in agreement with previous results from high- and low- frequency surveys. This flattening has been argued to be due to an increasing contribution of star-forming galaxies and faint active galactic nuclei. Additionally, we use our observations to evaluate the contribution of cosmic variance to the scatter in source counts measurements. The latter is achieved by dividing our Bo\"otes mosaic into 10 non-overlapping circular sectors, each one with an approximate area of $2\:\textrm{deg}^{2}.$ The counts in each sector are computed in the same way as done for the entire mosaic. By comparing the induced scatter with that of counts obtained from depth observations scaled to 150MHz, we find that the $1\sigma$ scatter due to cosmic variance is larger than the Poissonian errors of the source counts, and it may explain the dispersion from previously reported depth source counts at flux densities $S<1\,\textrm{mJy}$. This work demonstrates the feasibility of achieving deep radio imaging at low-frequencies with LOFAR.

Atomic Carbon can be an efficient tracer of the molecular gas mass, and when combined to the detection of high-J and low-J CO lines it yields also a sensitive probe of the power sources in the molecular gas of high redshift galaxies. The recently installed SEPIA5 receiver at the focus of the APEX telescope has opened up a new window at frequencies 159 - 211 GHz allowing the exploration of the Atomic Carbon in high-z galaxies, at previously inaccessible frequencies from the ground. We have targeted three gravitationally lensed galaxies at redshift of about 3 and conducted a comparative study of the observed high-J CO/CI ~ratios with well-studied nearby galaxies. Atomic Carbon (CI(2-1)) was detected in one of the three targets and marginally in a second, while in all three targets the $J=7\to6$ CO line is detected. The CO(7-6)/CI(2-1), CO(7-6)/CO(1-0) line ratios and the CO(7-6)/(far-IR continuum) luminosity ratio are compared to those of nearby objects. A large excitation status in the ISM of these high-z objects is seen, unless differential lensing unevenly boosts the CO line fluxes from the warm and dense gas more than the CO(1-0), CI(2-1), tracing a more widely distributed cold gas phase. We provide estimates of total molecular gas masses derived from the atomic Carbon and the Carbon monoxide CO(1-0), which within the uncertainties turn out to be equal.

We investigate the clustering properties of 45441 radio-quiet quasars (RQQs) and 3493 radio-loud quasars (RLQs) drawn from a joint use of the Sloan Digital Sky Survey (SDSS) and Faint Images of the Radio Sky at 20 cm (FIRST) surveys in the range $0.3<z<2.3$. This large spectroscopic quasar sample allow us to investigate the clustering signal dependence on radio-loudness and black hole (BH) virial mass. We find that RLQs are clustered more strongly than RQQs in all the redshift bins considered. We find a real-space correlation length of $r_{0}=6.59_{-0.24}^{+0.33}\,h^{-1}\,\textrm{Mpc}$ and $r_{0}=10.95_{-1.58}^{+1.22}\,h^{-1}\,\textrm{Mpc}$ {\normalsize{}for} RQQs and RLQs, respectively, for the full redshift range. This implies that RLQs are found in more massive host haloes than RQQs in our samples, with mean host halo masses of $\sim4.9\times10^{13}\,h^{-1}\,M_{\odot}$ and $\sim1.9\times10^{12}\,h^{-1}\,M_{\odot}$, respectively. Comparison with clustering studies of different radio source samples indicates that this mass scale of $\gtrsim1\times10^{13}\,h^{-1}\,M_{\odot}$ is characteristic for the bright radio-population, which corresponds to the typical mass of galaxy groups and galaxy clusters. The similarity we find in correlation lengths and host halo masses for RLQs, radio galaxies and flat-spectrum radio quasars agrees with orientation-driven unification models. Additionally, the clustering signal shows a dependence on black hole (BH) mass, with the quasars powered by the most massive BHs clustering more strongly than quasars having less massive BHs. We suggest that the current virial BH mass estimates may be a valid BH proxies for studying quasar clustering. We compare our results to a previous theoretical model that assumes that quasar activity

The LOFAR Two-metre Sky Survey (LoTSS) is a deep 120-168 MHz imaging survey that will eventually cover the entire Northern sky. Each of the 3170 pointings will be observed for 8 hrs, which, at most declinations, is sufficient to produce ~5arcsec resolution images with a sensitivity of ~0.1mJy/beam and accomplish the main scientific aims of the survey which are to explore the formation and evolution of massive black holes, galaxies, clusters of galaxies and large-scale structure. Due to the compact core and long baselines of LOFAR, the images provide excellent sensitivity to both highly extended and compact emission. For legacy value, the data are archived at high spectral and time resolution to facilitate subarcsecond imaging and spectral line studies. In this paper we provide an overview of the LoTSS. We outline the survey strategy, the observational status, the current calibration techniques, a preliminary data release, and the anticipated scientific impact. The preliminary images that we have released were created using a fully-automated but direction-independent calibration strategy and are significantly more sensitive than those produced by any existing large-area low-frequency survey. In excess of 44,000 sources are detected in the images that have a resolution of 25arcsec, typical noise levels of less than 0.5 mJy/beam, and cover an area of over 350 square degrees in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45d00m00s to 57d00m00s).

The Lockman Hole is a well-studied extragalactic field with extensive multi-band ancillary data covering a wide range in frequency, essential for characterising the physical and evolutionary properties of the various source populations detected in deep radio fields (mainly star-forming galaxies and AGNs). In this paper we present new 150-MHz observations carried out with the LOw Frequency ARray (LOFAR), allowing us to explore a new spectral window for the faint radio source population. This 150-MHz image covers an area of 34.7 square degrees with a resolution of 18.6$\times$14.7 arcsec and reaches an rms of 160 $\mu$Jy beam$^{-1}$ at the centre of the field. As expected for a low-frequency selected sample, the vast majority of sources exhibit steep spectra, with a median spectral index of $\alpha_{150}^{1400}=-0.78\pm0.015$. The median spectral index becomes slightly flatter (increasing from $\alpha_{150}^{1400}=-0.84$ to $\alpha_{150}^{1400}=-0.75$) with decreasing flux density down to $S_{150} \sim$10 mJy before flattening out and remaining constant below this flux level. For a bright subset of the 150-MHz selected sample we can trace the spectral properties down to lower frequencies using 60-MHz LOFAR observations, finding tentative evidence for sources to become flatter in spectrum between 60 and 150 MHz. Using the deep, multi-frequency data available in the Lockman Hole, we identify a sample of 100 Ultra-steep spectrum (USS) sources and 13 peaked spectrum sources. We estimate that up to 21 percent of these could have $z>4$ and are candidate high-$z$ radio galaxies, but further follow-up observations are required to confirm the physical nature of these objects.

We present the first wide area (19 deg$^2$), deep ($\approx120-150$ {\mu}Jy beam$^{-1}$), high resolution ($5.6 \times 7.4$ arcsec) LOFAR High Band Antenna image of the Bo\"otes field made at 130-169 MHz. This image is at least an order of magnitude deeper and 3-5 times higher in angular resolution than previously achieved for this field at low frequencies. The observations and data reduction, which includes full direction-dependent calibration, are described here. We present a radio source catalogue containing 6276 sources detected over an area of $19$\,deg$^2$, with a peak flux density threshold of $5\sigma$. As the first thorough test of the facet calibration strategy, introduced by van Weeren et al., we investigate the flux and positional accuracy of the catalogue. We present differential source counts that reach an order of magnitude deeper in flux density than previously achieved at these low frequencies, and show flattening at 150 MHz flux densities below 10 mJy associated with the rise of the low flux density star-forming galaxies and radio-quiet AGN.

With Low-Frequency Array (LOFAR) observations, we have discovered a diverse assembly of steep spectrum emission that is apparently associated with the intra cluster medium (ICM) of the merging galaxy cluster Abell 2034. Such a rich variety of complex emission associated with the ICM has been observed in few other clusters. This not only indicates that Abell 2034 is a more interesting and complex system than previously thought but it also demonstrates the importance of sensitive and high-resolution, low-frequency observations. These observations can reveal emission from relativistic particles which have been accelerated to sufficient energy to produce observable emission or have had their high energy maintained by mechanisms in the ICM. The most prominent feature in our maps is a bright bulb of emission connected to two steep spectrum filamentary structures, the longest of which extends perpendicular to the merger axis for 0.5Mpc across the south of the cluster. The origin of these objects is unclear, with no shock detected in the X-ray images and no obvious connection with cluster galaxies or AGNs. We also find that the X-ray bright region of the cluster coincides with a giant radio halo with an irregular morphology and a very steep spectrum. In addition, the cluster hosts up to three possible radio relics, which are misaligned with the cluster X-ray emission. Finally, we have identified multiple regions of emission with a very steep spectral index that seem to be associated with either tailed radio galaxies or a shock.

A metric representing a slowly rotating object with quadrupole moment is obtained using a perturbation method to include rotation into the weak limit of the Erez-Rosen metric. This metric is intended to tackle relativistic astrometry and gravitational lensing problems in which a quadrupole moment has to be taken into account.

N-body simulations predict that dark matter haloes are described by specific density profiles on both galactic- and cluster-sized scales. Weak gravitational lensing through the measurements of their first and second order properties, shear and flexion, is a powerful observational tool for investigating the true shape of these profiles. One of the three-parameter density profiles recently favoured in the description of dark matter haloes is the Einasto profile. We present exact expressions for the shear and the first and second flexions of Einasto dark matter haloes derived using a Mellin-transform formalism in terms of the Fox H and Meijer G functions, that are valid for general values of the Einasto index. The resulting expressions can be written as series expansions that permit us to investigate the asymptotic behaviour of these quantities. Moreover, we compare the shear and flexion of the Einasto profile with those of different mass profiles including the singular isothermal sphere, the Navarro-Frenk-White profile, and the S\'ersic profile. We investigate the concentration and index dependences of the Einasto profile, finding that the shear and second flexion could be used to determine the halo concentration, whilst for the Einasto index the shear and first and second flexions may be employed. We also provide simplified expressions for the weak lensing properties and other lensing quantities in terms of the generalized hypergeometric function.

Recent high-resolution N-body CDM simulations indicate that nonsingular three-parameter models such as the Einasto profile perform better than the singular two-parameter models, e.g. the Navarro, Frenk and White, in fitting a wide range of dark matter haloes. While many of the basic properties of the Einasto profile have been discussed in previous studies, a number of analytical properties are still not investigated. In particular, a general analytical formula for the surface density, an important quantity that defines the lensing properties of a dark matter halo, is still lacking to date. To this aim, we used a Mellin integral transform formalism to derive a closed expression for the Einasto surface density and related properties in terms of the Fox H and Meijer G functions, which can be written as series expansions. This enables arbitrary-precision calculations of the surface density and the lensing properties of realistic dark matter halo models. Furthermore, we compared the S\'ersic and Einasto surface mass densities and found differences between them, which implies that the lensing properties for both profiles differ.

In recent high resolution N-body CDM simulations, it has been had found that nonsingular three-parameter models, e.g. the Einasto profile has a better performance better than the singular two-parameter models, e.g. the Navarro, Frenk and White in the fitting of a wide range of dark matter halos. A problem with this profile is that the surface mass density is non-analytical for general values of the Einasto index. Therefore, its other lensing properties have the same problem. We obtain an exact analytical expression for the surface mass density of the Einasto profile in terms of the Fox H-function for all values of the Einasto index. With the idea of facilitate the use of the Einasto profile in lensing studies, we calculate the surface mass density, deflection angle, lens equation, deflection potential, magnification, shear and critical curves of the Einasto profile in terms of the Meijer G-function for all rational values of the Einasto index. The Meijer G-function have been implemented in several commercial and open-source computer algebra systems, thus the use of the lensing properties of the Einasto profile in strong and weak lensing studies is straighforward. We also compare the S\'ersic and Einasto surface mass densities profiles and found differences between them. This implies that the lensing properties are not equal for both profiles.