Search

We present a study of a luminous, z=0.15, type-2 quasar (log [L([OIII])/(erg/s)]=42.8) from the Quasar Feedback Survey. It is classified as 'radio-quiet' (log [L(1.4 GHz)/(W/Hz)]=23.8); however, radio imaging reveals ~1 kpc low-power jets (log [Pjet/(erg/s)]=44) inclined into the plane of the galaxy disk. We combine MUSE and ALMA observations to map stellar kinematics and ionised and molecular gas properties. The jets are seen to drive galaxy-wide bi-conical turbulent outflows, reaching W80 = 1000-1300 km/s, in the ionised phase (traced via optical emission-lines), which also have increased electron densities compared to the quiescent gas. The turbulent gas is driven perpendicular to the jet axis and is escaping along the galaxy minor axis, reaching 7.5 kpc on both sides. Traced via CO(3-2) emission, the turbulent material in molecular gas phase is one-third as spatially extended and has 3 times lower velocity-dispersion as compared to ionised gas. The jets are seen to be strongly interacting with the interstellar medium (ISM) through enhanced ionised emission and disturbed/depleted molecular gas at the jet termini. We see further evidence for jet-induced feedback through significantly higher stellar velocity-dispersion aligned, and co-spatial with, the jet axis (, Accepted for publication in MNRAS. Main manuscript has 21 pages with 8 figures. Supplementary material is available for download under "Ancillary files" or by downloading the source file listed under "Other formats"

Source classifications, stellar masses and star formation rates are presented for 80,000 radio sources from the first data release of the Low Frequency Array Two-metre Sky Survey (LoTSS) Deep Fields, which represents the widest deep radio survey ever undertaken. Using deep multi-wavelength data spanning from the ultraviolet to the far-infrared, spectral energy distribution (SED) fitting is carried out for all of the LoTSS-Deep host galaxies using four different SED codes, two of which include modelling of the contributions from an active galactic nucleus (AGN). Comparing the results of the four codes, galaxies that host a radiative AGN are identified, and an optimised consensus estimate of the stellar mass and star-formation rate for each galaxy is derived. Those galaxies with an excess of radio emission over that expected from star formation are then identified, and the LoTSS-Deep sources are divided into four classes: star-forming galaxies, radio-quiet AGN, and radio-loud high-excitation and low-excitation AGN. Ninety-five per cent of the sources can be reliably classified, of which more than two-thirds are star-forming galaxies, ranging from normal galaxies in the nearby Universe to highly-starbursting systems at z>4. Star-forming galaxies become the dominant population below 150-MHz flux densities of about 1 mJy, accounting for 90 per cent of sources at a 150-MHz flux density of 100 microJy. Radio-quiet AGN comprise around 10 per cent of the overall population. Results are compared against the predictions of the SKADS and T-RECS radio sky simulations, and improvements to the simulations are suggested., Accepted for publication in MNRAS. Catalogues available at www.lofar-surveys.org/deepfields.html. 27 pages

Recent studies have suggested that red quasars are a phase in quasar evolution when feedback from black hole accretion evacuates obscuring gas from the nucleus of the host galaxy. Here, we report a direct link between dust-reddening and molecular outflows in quasars at $z\sim2.5$. By examining the dynamics of warm molecular gas in the inner region of galaxies, we detect outflows with velocities 500--1000 km s$^{-1}$ and infer timescales of $\approx0.1$ Myr that are due to ongoing quasar energy output. We observe outflows only in systems where quasar radiation pressure on dust in the vicinity of the black hole is sufficiently large to expel their obscuring gas column densities. This result is in agreement with theoretical models that predict radiative feedback regulates gas in the nuclear regions of galaxies and is a major driving mechanism of galactic-scale outflows of cold gas. Our findings suggest that radiative quasar feedback ejects star-forming gas from within nascent stellar bulges at velocities comparable to those seen on larger scales, and that molecules survive in outflows even from the most luminous quasars., Comment: Accepted by MNRAS. 18 figures and 3 tables

Identifying active galactic nuclei (AGN) and isolating their contribution to a galaxy's energy budget is crucial for studying the co-evolution of AGN and their host galaxies. Brightness temperature ($T_b$) measurements from high-resolution radio observations at GHz frequencies are widely used to identify AGN. Here we investigate using new sub-arcsecond imaging at 144 MHz with the International LOFAR Telescope to identify AGN using $T_b$ in the Lockman Hole field. We use ancillary data to validate the 940 AGN identifications, finding 83 percent of sources have AGN classifications from SED fitting and/or photometric identifications, yielding 160 new AGN identifications. Considering the multi-wavelength classifications, brightness temperature criteria select over half of radio-excess sources, 32 percent of sources classified as radio-quiet AGN, and 20 percent of sources classified as star-forming galaxies. Infrared colour-colour plots and comparison with what we would expect to detect based on peak brightness in 6 arcsec LOFAR maps, imply that the star-forming galaxies and sources at low flux densities have a mixture of star-formation and AGN activity. We separate the radio emission from star-formation and AGN in unresolved, $T_b$-identified AGN with no significant radio excess and find the AGN comprises $0.49\pm 0.16$ of the radio luminosity. Overall the non-radio excess AGN show evidence for having a variety of different radio emission mechanisms, which can provide different pathways for AGN and galaxy co-evolution. This validation of AGN identification using brightness temperature at low frequencies opens the possibility for securely selecting AGN samples where ancillary data is inadequate., Accepted for publication in MNRAS. 19 pages, 15 figures, 2 tables

Red quasi-stellar objects (QSOs) are a subset of the luminous end of the cosmic population of active galactic nuclei (AGN), most of which are reddened by intervening dust along the line-of-sight towards their central engines. In recent work from our team, we developed a systematic technique to select red QSOs from the Sloan Digital Sky Survey (SDSS), and demonstrated that they have distinctive radio properties using the Faint Images of the Radio Sky at Twenty centimeters (FIRST) radio survey. Here we expand our study using low-frequency radio data from the LOFAR Two-metre Sky Survey (LoTSS). With the improvement in depth that LoTSS offers, we confirm key results: compared to a control sample of normal "blue" QSOs matched in redshift and accretion power, red QSOs have a higher radio detection rate and a higher incidence of compact radio morphologies. For the first time, we also demonstrate that these differences arise primarily in sources of intermediate radio-loudness: radio-intermediate red QSOs are $\times 3$ more common than typical QSOs, but the excess diminishes among the most radio-loud and the most radio-quiet systems in our study. We develop Monte-Carlo simulations to explore whether differences in star formation could explain these results, and conclude that, while star formation is an important source of low-frequency emission among radio-quiet QSOs, a population of AGN-driven compact radio sources is the most likely cause for the distinct low-frequency radio properties of red QSOs. Our study substantiates the conclusion that fundamental differences must exist between the red and normal blue QSO populations., Comment: 20 pages, 16 Figures. Accepted for publication in MNRAS

The far-infrared radio correlation (FIRC) is one of the strongest correlations in astronomy, yet a model that explains this comprehensively does not exist. The new LOFAR all Sky Survey (LoTSS) deep field, ELAIS-N1, allows exploration of this relation in previously unexplored regions of parameter space of radio frequency (150 MHz), luminosity (L150 24.7), redshift (z ∼ 1), and stellar mass M* 11.4. We present accurate deblended far-infrared (FIR) flux measurements with robust errors at 24, 100, 160, 250, 350, and 500 μm from Spitzer and the Herschel Space Observatory using XID+. We find that the FIRC has a strong mass dependence, the evolution of which takes the form qTIR(M*) = (2.00 ± 0.01)+(−0.22 ± 0.02)(log(M/M*)−10.05). This matches recent findings in regards to the star formation rate–radio luminosity relation at 150 MHz and results from radio observations in COSMOS at 1.4 GHz with the Jansky Very Large Array (JVLA). Our results provide tighter constraints on the low-redshift end of the FIRC and at lower frequency than the COSMOS observations. In addition, we find a mild evolution with redshift, with a best fit relation qTIR(z) = (1.94 ± 0.01)(1 + z)−0.04 ± 0.01. This evolution is shallower than that suggested by previous results at 150 MHz with the differences explained by the fact that previous studies did not account for the mass dependence. Finally, we present deblended FIR fluxes for 79 609 galaxies across the LoTSS deep fields: Boötes, ELAIS-N1, and Lockman Hole.

Optically luminous quasars at $z > 5$ are important probes of super-massive black hole (SMBH) formation. With new and future radio facilities, the discovery of the brightest low-frequency radio sources in this epoch would be an important new probe of cosmic reionization through 21-cm absorption experiments. In this work, we systematically study the low-frequency radio properties of a sample of 115 known spectroscopically confirmed $z>5$ quasars using the second data release of the Low Frequency Array (LOFAR) Two Metre Sky survey (LoTSS-DR2), reaching noise levels of $\sim$80 $\mu$Jy beam$^{-1}$ (at 144 MHz) over an area of $\sim5720$ deg$^2$. We find that 41 sources (36%) are detected in LoTSS-DR2 at $>2 \sigma$ significance and we explore the evolution of their radio properties (power, spectral index, and radio loudness) as a function of redshift and rest-frame ultra-violet properties. We obtain a median spectral index of $-0.29^{+0.10}_{-0.09}$ by stacking 93 quasars using LoTSS-DR2 and Faint Images of the Radio Sky at Twenty Centimetres (FIRST) data at 1.4 GHz, in line with observations of quasars at $z, Comment: 11 pages, 10 figures, accepted for publication in A&A

We examine the distribution of radio emission from ~42,000 quasars from the Sloan Digital Sky Survey, as measured in the LOFAR Two-Metre Sky Survey (LoTSS). We present a model of the radio luminosity distribution of the quasars that assumes that every quasar displays a superposition of two sources of radio emission: active galactic nuclei (jets) and star-formation. Our two-component model provides an excellent match to the observed radio flux density distributions across a wide range of redshifts and quasar optical luminosities; this suggests that the jet-launching mechanism operates in all quasars but with different powering efficiency. The wide distribution of jet powers allows for a smooth transition between the 'radio-quiet' and 'radio-loud' quasar regimes, without need for any explicit bimodality. The best-fit model parameters indicate that the star-formation rate of quasar host galaxies correlates strongly with quasar luminosity and also increases with redshift at least out to z~2. For a model where star-formation rate scales as $SFR \propto L_{bol}^\alpha (1+z)^\beta$, we find $\alpha = 0.47 \pm 0.01$ and $\beta = 1.61 \pm 0.05$, in agreement with far-infrared studies. Quasars contribute ~0.15 per cent of the cosmic star-formation rate density at z=0.5, rising to 0.4 per cent by z=2. The typical radio jet power is seen to increase with both increasing optical luminosity and black hole mass independently, but does not vary with redshift, suggesting intrinsic properties govern the production of the radio jets. We discuss the implications of these results for the triggering of quasar activity and the launching of jets., Comment: Accepted for publication in MNRAS. 21 pages

Fundamental differences in the radio properties of red quasars (QSOs), as compared to blue QSOs, have been recently discovered, positioning them as a potential key population in the evolution of galaxies and black holes across cosmic time. To elucidate their nature, we exploited a rich compilation of photometry and spectroscopic data to model their spectral energy distributions (SEDs) from the UV to the FIR and characterise their emission-line properties. Following a systematic comparison approach, we infer the AGN accretion, obscuration, and host galaxy properties in a sample of ~1800 QSOs at 0.21000 km/s) in red QSOs. We find that red QSOs that exhibit evidence for high-velocity winds present a stronger signature of the infrared excess, suggesting a causal connection between reddening and the presence of hot dust in QSO winds. We propose that dusty winds at nuclear scales are potentially the physical ingredient responsible for the colours in red QSOs, as well as a key parameter for the regulation of accretion material in the nucleus., Comment: Accepted for publication in Astronomy & Astrophysics. 22 pages, 16 figures. Abstract abridged

We present the first results from the Quasar Feedback Survey, a sample of 42 z10^42.1 ergs/s) with moderate radio luminosities (i.e. L(1.4GHz)>10^23.4 W/Hz; median L(1.4GHz)=5.9x10^23 W/Hz). Using high spatial resolution (~0.3-1 arcsec), 1.5-6 GHz radio images from the Very Large Array, we find that 67 percent of the sample have spatially extended radio features, on ~1-60 kpc scales. The radio sizes and morphologies suggest that these may be lower radio luminosity versions of compact, radio-loud AGN. By combining the radio-to-infrared excess parameter, spectral index, radio morphology and brightness temperature, we find radio emission in at least 57 percent of the sample that is associated with AGN-related processes (e.g. jets, quasar-driven winds or coronal emission). This is despite only 9.5-21 percent being classified as radio-loud using traditional criteria. The origin of the radio emission in the remainder of the sample is unclear. We find that both the established anti-correlation between radio size and the width of the [O III] line, and the known trend for the most [O III] luminous AGN to be associated with spatially-extended radio emission, also hold for our sample of moderate radio luminosity quasars. These observations add to the growing evidence of a connection between the radio emission and ionised gas in quasar host galaxies. This work lays the foundation for deeper investigations into the drivers and impact of feedback in this unique sample., Accepted for publication in MNRAS. Data products from this paper and the survey pilot papers are available through our website: https://blogs.ncl.ac.uk/quasarfeedbacksurvey/. The extensive supplementary material (containing additional figures and information on individual targets) is available for download under "Ancillary files" or by downloading the source file listed under "Other formats"

We have exploited LOFAR deep observations of the Lockman Hole field at 150 MHz to investigate the relation between the radio luminosity of star-forming galaxies (SFGs) and their star formation rates (SFRs), as well as its dependence on stellar mass and redshift. The adopted source classification, SFRs and stellar masses are consensus estimates based on a combination of four different SED fitting methods. We note a flattening of radio spectra of a substantial minority of sources below $\sim 1.4 $ GHz. Such sources have thus a lower "radio-loudness" level at 150 MHz than expected from extrapolations from 1.4 GHz using the average spectral index. We found a weak trend towards a lower $\hbox{SFR}/L_{150 \rm MHz}$ ratio for higher stellar mass, $M_\star$. We argue that such a trend may account for most of the apparent redshift evolution of the $L_{150 \rm MHz}/\hbox{SFR}$ ratio, in line with previous work. Our data indicate a weaker evolution than found by some previous analyses. We also find a weaker evolution with redshift of the specific star formation rate than found by several (but not all) previous studies. Our radio selection provides a view of the distribution of galaxies in the $\hbox{SFR}$-$M_\star$ plane complementary to that of optical/near-IR selection. It suggests a higher uniformity of the star formation history of galaxies than implied by some analyses of optical/near-IR data. We have derived luminosity functions at 150 MHz of both SFGs and radio-quiet (RQ) AGN at various redshifts. Our results are in very good agreement with the T-RECS simulations and with literature estimates. We also present explicit estimates of SFR functions of SFGs and RQ AGN at several redshifts derived from our radio survey data., Comment: 18 pages, 9 figures, 2 tables, A&A accepted

We investigate the radio properties of a sample of 53 sources selected at 850 $��$m from the SCUBA-2 Cosmology Legacy Survey using new deep, low-frequency radio imaging of the Lockman Hole field from the Low Frequency Array. Combining these data with additional radio observations from the GMRT and the JVLA, we find a variety of radio spectral shapes and luminosities within our sample despite their similarly bright submillimetre flux densities. We characterise their spectral shapes in terms of multi-band radio spectral indices. Finding strong spectral flattening at low frequencies in ~20% of sources, we investigate the differences between sources with extremely flat low-frequency spectra and those with `normal' radio spectral indices. As there are no other statistically significant differences between the two subgroups of our sample as split by the radio spectral index, we suggest that any differences are undetectable in galaxy-averaged properties that we can observe with our unresolved images, and likely relate to galaxy properties that we cannot resolve, on scales $\lesssim$ 1 kpc. We attribute the observed spectral flattening in the radio to free-free absorption, proposing that those sources with significant low-frequency spectral flattening have a clumpy distribution of star-forming gas. We estimate an average spatial extent of absorbing material of at most several hundred parsecs to produce the levels of absorption observed in the radio spectra. This estimate is consistent with the highest-resolution observations of submillimetre galaxies in the literature, which find examples of non-uniform dust distributions on scales of ~100 pc, with evidence for clumps and knots in the interstellar medium. Additionally, we find two bright (> 6 mJy) submm sources undetected at all other wavelengths. We speculate that these objects may be very high redshift sources, likely residing at z > 4., 15 pages, 10 figures, accepted by A&A

We study the impact of AGN ionised outflows on star formation in high-redshift AGN hosts, by combining NIR IFS observations, mapping the H$\alpha$ emission and [OIII] outflows, with matched-resolution observations of the rest-frame FIR emission. We present high-resolution ALMA Band 7 observations of eight X-ray selected AGN at z~2 from the SUPER sample, targeting the rest-frame ~260 um continuum at ~2 kpc (0.2'') resolution. We detected 6 out of 8 targets with S/N>10 in the ALMA maps, with continuum flux densities F = 0.27-2.58 mJy and FIR half-light radii Re = 0.8-2.1 kpc. The FIR Re of our sample are comparable to other AGN and star-forming galaxies at a similar redshift from the literature. However, we find that the mean FIR size in X-ray AGN (Re = 1.16+/- 0.11 kpc) is slightly smaller than in non-AGN (Re = 1.69+/-0.13 kpc). From SED fitting, we find that the main contribution to the 260 um flux density is dust heated by star formation, with < 4% contribution from AGN-heated dust and < 1% from synchrotron emission. The majority of our sample show different morphologies for the FIR (mostly due to reprocessed stellar emission) and the ionised gas emission (H$\alpha$ and [OIII], mostly due to AGN emission). This could be due to the different locations of dust and ionised gas, the different sources of the emission (stars and AGN), or the effect of dust obscuration. We are unable to identify any residual H$\alpha$ emission, above that dominated by AGN, that could be attributed to star formation. Under the assumption that the FIR emission is a reliable tracer of obscured star formation, we find that the obscured star formation activity in these AGN host galaxies is not clearly affected by the ionised outflows. However, we cannot rule out that star formation suppression is happening on smaller spatial scales than the ones we probe with our observations (< 2 kpc) or on different timescales., Comment: 32 pages, 13 figures. Accepted for publication in A&A

The Low Frequency Array (LOFAR) is an ideal instrument to conduct deep extragalactic surveys. It has a large field of view and is sensitive to large-scale and compact emission. It is, however, very challenging to synthesize thermal noise limited maps at full resolution, mainly because of the complexity of the low-frequency sky and the direction dependent effects (phased array beams and ionosphere). In this first paper of a series, we present a new calibration and imaging pipeline that aims at producing high fidelity, high dynamic range images with LOFAR High Band Antenna data, while being computationally efficient and robust against the absorption of unmodeled radio emission. We apply this calibration and imaging strategy to synthesize deep images of the Boötes and Lockman Hole fields at ~150 MHz, totaling ~80 and ~100 h of integration, respectively, and reaching unprecedented noise levels at these low frequencies of ≲30 and ≲23μJy beam−1in the inner ~3 deg2. This approach is also being used to reduce the LOTSS-wide data for the second data release.

We use a sample of powerful z~0.1 type 2 quasars ('obscured'; log[L(AGN)/erg/s]>~45), which host kiloparsec-scale ionized outflows and jets, to identify possible signatures of AGN feedback on the total molecular gas reservoirs of their host galaxies. Specifically, we present Atacama Pathfinder EXperiment (APEX) observations of the CO(2-1) transition for nine sources and the CO(6-5) for a subset of three. We find that the majority of our sample reside in starburst galaxies (average specific star formation rates of 1.7/Gyr), with the seven CO-detected quasars also having large molecular gas reservoirs (average Mgas = 1.3x10^10Msun), even though we had no pre-selection on the star formation or molecular gas properties. Despite the presence of quasars and outflows, we find that the molecular gas fractions (Mgas/Mstar = 0.1-1.2) and depletion times (Mgas/SFR = 0.16-0.95Gyr) are consistent with those expected for the overall galaxy population with matched stellar masses and specific star formation rates. Furthermore, for at least two of the three targets with the required measurements, the CO(6-5)/CO(2-1) emission-line ratios are consistent with star formation dominating the CO excitation over this range of transitions. The targets in our study represent a gas-rich phase of galaxy evolution with simultaneously high levels of star formation and nuclear activity; furthermore, the jets and outflows do not have an immediate appreciable impact on the global molecular gas reservoirs., Accepted for publication in MNRAS. Supplementary material is available for download under "Ancillary files"

Using data from ALMA and near-infrared (NIR) integral field spectrographs including both SINFONI and KMOS on the VLT, we investigate the two-dimensional distributions of H$\alpha$ and rest-frame far-infrared (FIR) continuum in six submillimeter galaxies at $z\sim2$. At a similar spatial resolution ($\sim$0.5" FWHM; $\sim$4.5 kpc at $z=2$), we find that the half-light radius of H$\alpha$ is significantly larger than that of the FIR continuum in half of the sample, and on average H$\alpha$ is a median factor of $2.0\pm0.4$ larger. Having explored various ways to correct for the attenuation, we find that the attenuation-corrected H$\alpha$-based SFRs are systematically lower than the IR-based SFRs by at least a median factor of $3\pm1$, which cannot be explained by the difference in half-light radius alone. In addition, we find that in 40% of cases the total $V$-band attenuation ($A_V$) derived from energy balance modeling of the full ultraviolet(UV)-to-FIR spectral energy distributions (SEDs) is significantly higher than that derived from SED modeling using only the UV-to-NIR part of the SEDs, and the discrepancy appears to increase with increasing total infrared luminosity. Finally, considering all our findings along with the studies in the literature, we postulate that the dust distributions in SMGs, and possibly also in less IR luminous $z\sim2$ massive star-forming galaxies, can be decomposed into three main components; the diffuse dust heated by older stellar populations, the more obscured and extended young star-forming HII regions, and the heavily obscured central regions that have a low filling factor but dominate the infrared luminosity in which the majority of attenuation cannot be probed via UV-to-NIR emissions., Comment: 15 pages, 8 figures, A&A in press

We have recently used the Faint Images of the Radio Sky at Twenty-centimeters (FIRST) survey to show that red quasars have fundamentally different radio properties to typical blue quasars: a significant (factor $\sim3$) enhancement in the radio-detection fraction, which arises from systems around the radio-quiet threshold with compact ($, Comment: 17 pages, 17 Figures. Accepted for publication in MNRAS

We present 0.″15 (1 kpc) resolution ALMA observations of the [C II] 157.74 μm line and rest-frame 160 μm continuum emission in two z ∼ 3 dusty, star-forming galaxies—ALESS 49.1 and ALESS 57.1, combined with resolved CO (3-2) observations. In both sources, the [C II] surface brightness distribution is dominated by a compact core ≤1 kpc in radius, a factor of 2-3 smaller than the extent of the CO (3-2) emission. In ALESS 49.1, we find an additional extended (8 kpc radius), low surface brightness [C II] component. Based on an analysis of mock ALMA observations, the [C II] and 160 μm continuum surface brightness distributions are inconsistent with a single-Gaussian surface brightness distribution with the same size as the CO (3-2) emission. The [C II] rotation curves flatten at ≃2 kpc radius, suggesting that the kinematics of the central regions are dominated by a baryonic disk. Both galaxies exhibit a strong [C II]/far-IR (FIR) deficit on 1 kpc scales, with FIR surface brightness to [C II]/FIR slope steeper than in local star- forming galaxies. A comparison of the [C II]/CO (3-2) observations with photodissociation region models suggests a strong far-UV (FUV) radiation field (G 0 ∼ 104) and high gas density (n(H) ∼ 104-105 cm−3) in the central regions of ALESS 49.1 and ALESS 57.1. The most direct interpretation of the pronounced [C II]/FIR deficit is a thermal saturation of the C+ fine-structure levels at temperatures ≥500 K, driven by the strong FUV field.

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., Comment: 16 figures, 1 table and 22 pages. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available on lofar-surveys.org upon publication of the journal

Determining the star-forming properties of radio-quiet quasars is important for understanding the co-evolution of star formation and black hole accretion. Here, we present the detection of the gravitationally-lensed radio-quiet quasars SDSS J1055+4628, SDSS J1313+5151 and SBS 1520+530 at 144 MHz that fall in the HETDEX Spring Field targeted in the LOFAR Two-metre Sky Survey (LoTSS) first full data release. We compare their radio and far-infrared luminosities relative to the radio-infrared correlation and find that their radio luminosities can be explained by star formation. The implied star formation rates derived from their radio and infrared luminosities are between 20 and 300 $\rm{M_{\odot}~yr^{-1}}$. These detections represent the first study of gravitationally lensed sources with LOFAR, opening a new frequency window for investigating the star-forming properties of high-redshift quasars at radio wavelengths. We consider the implications for future data releases and estimate that many of the objects in our parent sample will be detected during LoTSS, significantly increasing the fraction of gravitationally lensed radio-quiet quasars with radio detections., 7 figures, 2 tables and 10 pages. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available on lofar-surveys.org upon publication of the journal

We present sub-kpc-scale mapping of the 870 $\mu$m ALMA continuum emission in six luminous ($L_{\rm IR}~\sim~5~\times10^{12}$ L$_{\odot}$) submillimeter galaxies (SMGs) from the ALESS survey of the Extended Chandra Deep Field South. Our high-fidelity 0.07$''$-resolution imaging ($\sim$500 pc) reveals robust evidence for structures with deconvolved sizes of $\lesssim$0.5-1 kpc embedded within (dominant) exponential dust disks. The large-scale morphologies of the structures within some of the galaxies show clear curvature and/or clump-like structures bracketing elongated nuclear emission, suggestive of bars, star-forming rings, and spiral arms. In this interpretation, the ratio of the `ring' and `bar' radii (1.9$\pm$0.3) agrees with that measured for such features in local galaxies. These potential spiral/ring/bar structures would be consistent with the idea of tidal disturbances, with their detailed properties implying flat inner rotation curves and Toomre-unstable disks (Q, Comment: 18 pages, 6 figures, Accepted to ApJ

This paper presents a study of the redshift evolution of radio-loud active galactic nuclei (AGN) as a function of the properties of their galaxy hosts in the Bo\"otes field. To achieve this we match low-frequency radio sources from deep $150$-MHz LOFAR observations to an $I$-band-selected catalogue of galaxies, for which we have derived photometric redshifts, stellar masses and rest-frame colours. We present spectral energy distribution (SED) fitting to determine the mid-infrared AGN contribution for the radio sources and use this information to classify them as High- versus Low-Excitation Radio Galaxies (HERGs and LERGs) or Star-Forming galaxies. Based on these classifications we construct luminosity functions for the separate redshift ranges going out to $z = 2$. From the matched radio-optical catalogues, we select a sub-sample of $624$ high power ($P_{150\mathrm{\,MHz}}>10^{25}$ W Hz$^{-1}$) radio sources between $0.5 \leq z < 2$. For this sample, we study the fraction of galaxies hosting HERGs and LERGs as a function of stellar mass and host galaxy colour. The fraction of HERGs increases with redshift, as does the fraction of sources in galaxies with lower stellar masses. We find that the fraction of galaxies that host LERGs is a strong function of stellar mass as it is in the local Universe. This, combined with the strong negative evolution of the LERG luminosity functions over this redshift range, is consistent with LERGs being fuelled by hot gas in quiescent galaxies., Comment: 25 pages, 21 figures, accepted for publication in MNRAS

We present a study of the low-frequency radio properties of star forming (SF) galaxies and active galactic nuclei (AGN) up to redshift $z=2.5$. The new spectral window probed by the Low Frequency Array (LOFAR) allows us to reconstruct the radio continuum emission from 150 MHz to 1.4 GHz to an unprecedented depth for a radio-selected sample of $1542$ galaxies in $\sim 7~ \rm{deg}^2$ of the LOFAR Bo\"otes field. Using the extensive multi-wavelength dataset available in Bo\"otes and detailed modelling of the FIR to UV spectral energy distribution (SED), we are able to separate the star-formation (N=758) and the AGN (N=784) dominated populations. We study the shape of the radio SEDs and their evolution across cosmic time and find significant differences in the spectral curvature between the SF galaxy and AGN populations. While the radio spectra of SF galaxies exhibit a weak but statistically significant flattening, AGN SEDs show a clear trend to become steeper towards lower frequencies. No evolution of the spectral curvature as a function of redshift is found for SF galaxies or AGN. We investigate the redshift evolution of the infrared-radio correlation (IRC) for SF galaxies and find that the ratio of total infrared to 1.4 GHz radio luminosities decreases with increasing redshift: $ q_{\rm 1.4GHz} = (2.45 \pm 0.04) \times (1+z)^{-0.15 \pm 0.03} $. Similarly, $q_{\rm 150MHz}$ shows a redshift evolution following $ q_{\rm 150GHz} = (1.72 \pm 0.04) \times (1+z)^{-0.22 \pm 0.05}$. Calibration of the 150 MHz radio luminosity as a star formation rate tracer suggests that a single power-law extrapolation from $q_{\rm 1.4GHz}$ is not an accurate approximation at all redshifts., Comment: Accepted for publication in MNRAS, 21 pages, 16 figures

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)., Comment: 23 pages, 15 figures, 3 tables, accepted for publication in A&A

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., Comment: 26 pages. Accepted for publication in MNRAS