Search

The magnetised near detector (ND280) of the T2K long-baseline neutrino oscillation experiment has been recently upgraded aiming to satisfy the requirement of reducing the systematic uncertainty from measuring the neutrinonucleus interaction cross section, which is the largest systematic uncertainty in the search for leptonic charge-parity symmetry violation. A key component of the upgrade is SuperFGD, a 3D segmented plastic scintillator detector made of approximately 2,000,000 optically-isolated 1 cm3 cubes. It will provide a 3D image of GeV neutrino interactions by combining tracking and stopping power measurements of final state particles with sub-nanosecond time resolution. The performance of SuperFGD is characterized by the precision of its response to charged particles as well as the systematic effects that might affect the physics measurements. Hence, a detailed Geant4 based optical simulation of the SuperFGD building block, i.e. a plastic scintillating cube read out by three wavelength shifting fibers, has been developed and validated with the different datasets collected in various beam tests. In this manuscript the description of the optical model as well as the comparison with data are reported., Comment: 31 pages, 15 figures

We introduce a unified approach that, given a strong gravitationally lensed polarised source, self-consistently infers its complex surface brightness distribution and the lens galaxy mass-density profile, magnetic field and electron density from interferometric data. The method is fully Bayesian, pixellated and three-dimensional: the source light is reconstructed in each frequency channel on a Delaunay tessellation with a magnification-adaptive resolution. We tested this technique using simulated interferometric observations with a realistic model of the lens, for two different levels of source polarisation and two different lensing configurations. For all data sets, the presence of a Faraday rotating screen in the lens is supported by the data with strong statistical significance. In the region probed by the lensed images, we can recover the Rotation Measure and the parallel component of the magnetic field with an average error between 0.6 and 11 rad m$^{-2}$ and 0.3 and 3 nG, respectively. Given our choice of model, we find the electron density is the least well-constrained component due to a degeneracy with the magnetic field and disk inclination. The background source total intensity, polarisation fraction, and polarisation angle are inferred with an error between 4 and 10 per cent, 15 and 50 per cent, and 1 to 12 degrees, respectively. Our analysis shows that both the lensing configuration and the intrinsic model degeneracies play a role in the quality of the constraints that can be obtained., Comment: submitted for review to MNRAS

The Super-Kamiokande and T2K collaborations present a joint measurement of neutrino oscillation parameters from their atmospheric and beam neutrino data. It uses a common interaction model for events overlapping in neutrino energy and correlated detector systematic uncertainties between the two datasets, which are found to be compatible. Using 3244.4 days of atmospheric data and a beam exposure of $19.7(16.3) \times 10^{20}$ protons on target in (anti)neutrino mode, the analysis finds a 1.9$\sigma$ exclusion of CP-conservation (defined as $J_{CP}=0$) and a preference for the normal mass ordering., Comment: 12 pages, 4 figures

We present 6-GHz Very Large Array radio images of 70 gravitational lens systems at 300-mas resolution, in which the source is an optically-selected quasar, and nearly all of which have two lensed images. We find that about in half of the systems (40/70, with 33/70 secure), one or more lensed images are detected down to our detection limit of 20microJy/beam, similar to previous investigations and reinforcing the conclusion that typical optically-selected quasars have intrinsic GHz radio flux densities of a few microJy ($\sim10^{23}$WHz$^{-1}$ at redshifts of 1--2). In addition, for ten cases it is likely that the lensing galaxies are detected in the radio. Available detections of, and limits on the far-infrared luminosities from the literature, suggest that nearly all of the sample lie on the radio-FIR correlation typical of star-forming galaxies, and that their radio luminosities are at least compatible with the radio emission being produced by star formation processes. One object, WISE2329$-$1258, has an extra radio component that is not present in optical images, and is difficult to explain using simple lens models. In-band spectral indices, where these can be determined, are generally moderately steep and consistent with synchrotron processes either from star-formation/supernovae or AGN. Comparison of the A/B image flux ratios at radio and optical wavelengths suggests a 10 per cent level contribution from finite source effects or optical extinction to the optical flux ratios, together with sporadic larger discrepancies that are likely to be due to optical microlensing., Comment: Accepted by MNRAS

The large-scale mass distributions of galaxy-scale strong lenses have long been assumed to be well-described by a singular ellipsoidal power-law density profile with external shear. However, the inflexibility of this model could lead to systematic errors in astrophysical parameters inferred with gravitational lensing observables. Here, we present observations with the Atacama Large (sub-)Millimetre Array (ALMA) of three strongly lensed dusty star-forming galaxies at $\simeq30$ mas angular resolution and investigate the sensitivity of these data to angular structure in the lensing galaxies. We jointly infer the lensing mass distribution and the full surface brightness of the lensed sources with multipole expansions of the power-law density profile up to fourth order using a technique developed for interferometric data. All three data sets strongly favour third and fourth-order multipole amplitudes of $\approx1$ percent of the convergence. While the infrared stellar isophotes and isodensity shapes agree for one lens system, for the other two the isophotes disagree to varying extents, suggesting contributions to the angular structure from dark matter intrinsic or extrinsic to the lensing galaxy., Comment: In print with A&A

Unobscured quasars (QSOs) are predicted to be the final stage in the evolutionary sequence from gas-rich mergers to gas-depleted, quenched galaxies. Studies of this population, however, find a high incidence of far-infrared-luminous sources -suggesting significant dust-obscured star formation-but direct observations of the cold molecular gas fuelling this star formation are still necessary. We present a NOEMA study of CO(2-1) emission, tracing the cold molecular gas, in ten lensed z=1-1.5 unobscured QSOs. We detected CO(2-1) in seven of our targets, four of which also show continuum emission (\lambda_rest = 1.3mm). After subtracting the foreground galaxy contribution to the photometry, spectral energy distribution fitting yielded stellar masses of 10^9-11 M_\odot, with star formation rates of 25-160 M_\odot yr^-1 for the host galaxies. These QSOs have lower $L'_\mathrm{CO}$ than star-forming galaxies with the same L_IR, and show depletion times spanning a large range (50-900 Myr), but with a median of just 90 Myr. We find molecular gas masses in the range 2-40 x 10^9(alpha_CO/4) M_\odot, which suggest gas fractions above ~50% for most of the targets. Despite the presence of an unobscured QSO, the host galaxies are able to retain significant amounts of cold gas. However, with a median depletion time of ~90 Myr, the intense burst of star formation taking place in these targets will quickly deplete their molecular gas reservoirs in the absence of gas replenishment, resulting in a quiescent host galaxy. The non-detected QSOs are three of the four radio-loud QSOs in the sample, and their properties indicate that they are likely already transitioning into quiescence. Recent cosmological simulations tend to overestimate the depletion times expected for these z~1 QSO-host galaxies, which is likely linked to their difficulty producing starbursts across the general high-redshift galaxy population., Comment: 20 pages. Accepted for publication in A&A

We report an updated measurement of the $\nu_{\mu}$-induced, and the first measurement of the $\bar{\nu}_{\mu}$-induced coherent charged pion production cross section on $^{12}C$ nuclei in the T2K experiment. This is measured in a restricted region of the final-state phase space for which $p_{\mu,\pi} > 0.2$ GeV, $\cos(\theta_{\mu}) > 0.8$ and $\cos(\theta_{\pi}) > 0.6$, and at a mean (anti)neutrino energy of 0.85 GeV using the T2K near detector. The measured $\nu_{\mu}$ CC coherent pion production flux-averaged cross section on $^{12}C$ is $(2.98 \pm 0.37 (stat.) \pm 0.31 (syst.) \substack{ +0.49 \\ -0.00 } \mathrm{ (Q^2\,model)}) \times 10^{-40}~\mathrm{cm}^{2}$. The new measurement of the $\bar{\nu}_{\mu}$-induced cross section on $^{12}{C}$ is $(3.05 \pm 0.71 (stat.) \pm 0.39 (syst.) \substack{ +0.74 \\ -0.00 } \mathrm{(Q^2\,model)}) \times 10^{-40}~\mathrm{cm}^{2}$. The results are compatible with both the NEUT 5.4.0 Berger-Sehgal (2009) and GENIE 2.8.0 Rein-Sehgal (2007) model predictions.

We present an analysis of the detection fraction and the number counts of radio sources imaged with Very Long Baseline Interferometry (VLBI) at 1.4 GHz as part of the mJIVE-20 survey. From a sample of 24,903 radio sources identified by FIRST, 4,965 are detected on VLBI-scales, giving an overall detection fraction of $19.9\pm2.9~$per cent. However, we find that the detection fraction falls from around 50 per cent at a peak surface brightness of $80~mJy~beam^{-1}$ in FIRST to around 8 per cent at the detection limit, which is likely dominated by the surface brightness sensitivity of the VLBI observations, with some contribution from a change in the radio source population. We also find that compactness at arcsec-scales is the dominant factor in determining whether a radio source is detected with VLBI, and that the median size of the VLBI-detected radio sources is 7.7 mas. After correcting for the survey completeness and effective sky area, we determine the slope of the differential number counts of VLBI-detected radio sources with flux densities $S_{\rm 1.4~GHz} > 1~mJy$ to be $\eta_{\rm VLBI} = -1.74\pm 0.02$, which is shallower than in the cases of the FIRST parent population ($\eta_{\rm FIRST} = -1.77\pm 0.02$) and for compact radio sources selected at higher frequencies ($\eta_{\rm JBF} = -2.06\pm 0.02$). From this, we find that all-sky ($3\pi~sr$) surveys with the EVN and the VLBA have the potential to detect $(7.2\pm0.9)\times10^{5}$ radio sources at mas-resolution, and that the density of compact radio sources is sufficient (5.3~deg$^{-2}$) for in-beam phase referencing with multiple sources (3.9 per primary beam) in the case of a hypothetical SKA-VLBI array.

We present deep ALMA Band 3 observations of the HCN, HCO+, and HNC (4-3) emission in SDP.81, a well-studied z = 3.042 strongly lensed galaxy. These lines trace the high-density gas, which remains almost entirely unexplored in z$\geq$1 galaxies. Additionally, these dense-gas tracers are potentially powerful diagnostics of the mechanical heating of the interstellar medium. While the HCN(4-3) and HNC(4-3) lines are not detected, the HCO+(4-3) emission is clearly detected and resolved. This is the third detection of this line in a high-redshift star-forming galaxy. We find an unusually high HCO+/HCN intensity ratio of $\geq$2.2. Based on the photodissociation region modelling, the most likely explanation for the elevated HCO+/HCN ratio is that SDP.81 has low mechanical heating - less than 10% of the total energy budget - and a sub-solar metallicity, Z=0.5 Z$_\odot$. While such conditions might not be representative of the general population of high-redshift dusty galaxies, lower-than-solar metallicity might have a significant impact on gas masses inferred from CO observations. In addition, we report the detection of CO(0-1) absorption from the foreground lensing galaxy and CO(1-0) emission from a massive companion to the lensing galaxy, approximately 50 kpc to the southeast., Comment: A&A accepted, corrected version

Dark matter structures within strong gravitational lens galaxies and along their line of sight leave a gravitational imprint on the multiple images of lensed sources. Strong gravitational lensing provides, therefore, a key test of different dark matter models in a way that is independent of the baryonic content of matter structures on subgalactic scales. In this chapter, we describe how galaxy-scale strong gravitational lensing observations are sensitive to the physical nature of dark matter. We provide a historical perspective of the field, and review its current status. We discuss the challenges and advances in terms of data, treatment of systematic errors and theoretical predictions, that will enable one to deliver a stringent and robust test of different dark matter models in the near future. With the advent of the next generation of sky surveys, the number of known strong gravitational lens systems is expected to increase by several orders of magnitude. Coupled with high-resolution follow-up observations, these data will provide a key opportunity to constrain the properties of dark matter with strong gravitational lensing.

Muon neutrino and antineutrino disappearance probabilities are identical in the standard three-flavor neutrino oscillation framework, but CPT violation and non-standard interactions can violate this symmetry. In this work we report the measurements of $\sin^{2} \theta_{23}$ and $\Delta m_{32}^2$ independently for neutrinos and antineutrinos. The aforementioned symmetry violation would manifest as an inconsistency in the neutrino and antineutrino oscillation parameters. The analysis discussed here uses a total of 1.97$\times$10$^{21}$ and 1.63$\times$10$^{21}$ protons on target taken with a neutrino and antineutrino beam respectively, and benefits from improved flux and cross-section models, new near detector samples and more than double the data reducing the overall uncertainty of the result. No significant deviation is observed, consistent with the standard neutrino oscillation picture.

This paper reports the first measurement of muon neutrino charged-current interactions without pions in the final state using multiple detectors with correlated energy spectra at T2K. The data was collected on hydrocarbon targets using the off-axis T2K near detector (ND280) and the on-axis T2K near detector (INGRID) with neutrino energy spectra peaked at 0.6 GeV and 1.1 GeV respectively. The correlated neutrino flux presents an opportunity to reduce the impact of the flux uncertainty and to study the energy dependence of neutrino interactions. The extracted double-differential cross sections are compared to several Monte Carlo neutrino-nucleus interaction event generators showing the agreement between both detectors individually and with the correlated result., Comment: Updated discussion in Sec. V-A; Updated author list

The T2K experiment presents new measurements of neutrino oscillation parameters using $19.7(16.3)\times10^{20}$ protons on target (POT) in (anti-)neutrino mode at the far detector (FD). Compared to the previous analysis, an additional $4.7\times10^{20}$ POT neutrino data was collected at the FD. Significant improvements were made to the analysis methodology, with the near-detector analysis introducing new selections and using more than double the data. Additionally, this is the first T2K oscillation analysis to use NA61/SHINE data on a replica of the T2K target to tune the neutrino flux model, and the neutrino interaction model was improved to include new nuclear effects and calculations. Frequentist and Bayesian analyses are presented, including results on $\sin^2\theta_{13}$ and the impact of priors on the $\delta_\mathrm{CP}$ measurement. Both analyses prefer the normal mass ordering and upper octant of $\sin^2\theta_{23}$ with a nearly maximally CP-violating phase. Assuming the normal ordering and using the constraint on $\sin^2\theta_{13}$ from reactors, $\sin^2\theta_{23}=0.561^{+0.021}_{-0.032}$ using Feldman--Cousins corrected intervals, and $\Delta{}m^2_{32}=2.494_{-0.058}^{+0.041}\times10^{-3}~\mathrm{eV^2}$ using constant $\Delta\chi^{2}$ intervals. The CP-violating phase is constrained to $\delta_\mathrm{CP}=-1.97_{-0.70}^{+0.97}$ using Feldman--Cousins corrected intervals, and $\delta_\mathrm{CP}=0,\pi$ is excluded at more than 90% confidence level. A Jarlskog invariant of zero is excluded at more than $2\sigma$ credible level using a flat prior in $\delta_\mathrm{CP}$, and just below $2\sigma$ using a flat prior in $\sin\delta_\mathrm{CP}$. When the external constraint on $\sin^2\theta_{13}$ is removed, $\sin^2\theta_{13}=28.0^{+2.8}_{-6.5}\times10^{-3}$, in agreement with measurements from reactor experiments. These results are consistent with previous T2K analyses.

Using a single gravitational lens system observed at $\lesssim5$ milli-arcsecond resolution with very long baseline interferometry (VLBI), we place a lower bound on the mass of the fuzzy dark matter (FDM) particle, ruling out $m_\chi \leq 4.4\times10^{-21}~\mathrm{eV}$ with a 20:1 posterior odds ratio relative to a smooth lens model. We generalize our result to non-scalar and multiple-field models, such as vector FDM, with $m_{\chi,\mathrm{vec}} > 1.4 \times 10^{-21}~\mathrm{eV}$. Due to the extended source structure and high angular resolution of the observation, our analysis is directly sensitive to the presence of granule structures in the main dark matter halo of the lens, which is the most generic prediction of FDM theories. A model based on well-understood physics of ultra-light dark matter fields in a gravitational potential well makes our result robust to a wide range of assumed dark matter fractions and velocity dispersions in the lens galaxy. Our result is competitive with other lower bounds on $m_\chi$ from past analyses, which rely on intermediate modelling of structure formation and/or baryonic effects. Higher resolution observations taken at 10 to 100 GHz could improve our constraints by up to 2 orders of magnitude in the future., Comment: 5 pages, 2 figures. Accepted in MNRAS Letters

While the direct detection of the dark-matter particle remains very challenging, the nature of dark matter could be possibly constrained by comparing the observed abundance and properties of small-scale sub-galactic mass structures with predictions from the phenomenological dark-matter models, such as cold, warm or hot dark matter. Galaxy-galaxy strong gravitational lensing provides a unique opportunity to search for tiny surface-brightness anomalies in the extended lensed images (i.e. Einstein rings or gravitational arcs), induced by possible small-scale mass structures in the foreground lens galaxy. In this paper, the first in a series, we introduce and test a methodology to measure the power spectrum of such surface-brightness anomalies from high-resolution Hubble Space Telescope (HST) imaging. In particular, we focus on the observational aspects of this statistical approach, such as the most suitable observational strategy and sample selection, the choice of modelling techniques and the noise correction. We test the feasibility of the power-spectrum measurement by applying it to a sample of galaxy-galaxy strong gravitational lens systems from the Sloan Lens ACS Survey, with the most extended, bright, high-signal-to-noise-ratio lensed images, observed in the rest frame ultraviolet. In the companion paper, we present the methodology to relate the measured power spectrum to the statistical properties of the underlying small-scale mass structures in the lens galaxy and infer the first observational constraints on the sub-galactic matter power spectrum in a massive elliptical (lens) galaxy., Comment: 20 pages, 19 figures, 2 tables, submitted to MNRAS

We present a novel machine learning based approach for detecting galaxy-scale gravitational lenses from interferometric data, specifically those taken with the International LOFAR Telescope (ILT), which is observing the northern radio sky at a frequency of 150 MHz, an angular resolution of 350 mas and a sensitivity of 90 uJy beam-1 (1 sigma). We develop and test several Convolutional Neural Networks to determine the probability and uncertainty of a given sample being classified as a lensed or non-lensed event. By training and testing on a simulated interferometric imaging data set that includes realistic lensed and non-lensed radio sources, we find that it is possible to recover 95.3 per cent of the lensed samples (true positive rate), with a contamination of just 0.008 per cent from non-lensed samples (false positive rate). Taking the expected lensing probability into account results in a predicted sample purity for lensed events of 92.2 per cent. We find that the network structure is most robust when the maximum image separation between the lensed images is greater than 3 times the synthesized beam size, and the lensed images have a total flux density that is equivalent to at least a 20 sigma (point-source) detection. For the ILT, this corresponds to a lens sample with Einstein radii greater than 0.5 arcsec and a radio source population with 150 MHz flux densities more than 2 mJy. By applying these criteria and our lens detection algorithm we expect to discover the vast majority of galaxy-scale gravitational lens systems contained within the LOFAR Two Metre Sky Survey., Comment: Accepted to be published by MNRAS

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

We investigate the mass structure of a strong gravitational lens galaxy at $z=0.350$, taking advantage of the milli-arcsecond (mas) angular resolution of very long baseline interferometric (VLBI) observations. In the first analysis of its kind at this resolution, we jointly infer the lens model parameters and pixellated radio source surface brightness. We consider several lens models of increasing complexity, starting from an elliptical power-law density profile. We extend this model to include angular multipole structures, a separate stellar mass component, additional nearby field galaxies, and/or a generic external potential. We compare these models using their relative Bayesian log-evidence (Bayes factor). We find strong evidence for angular structure in the lens; our best model is comprised of a power-law profile plus multipole perturbations and external potential, with a Bayes factor of $+14984$ relative to the elliptical power-law model. It is noteworthy that the elliptical power-law mass distribution is a remarkably good fit on its own, with additional model complexity correcting the deflection angles only at the $\sim5$ mas level. We also consider the effects of added complexity in the lens model on time-delay cosmography and flux-ratio analyses. We find that an overly simplistic power-law ellipsoid lens model can bias the measurement of $H_0$ by $\sim3$ per cent and mimic flux ratio anomalies of $\sim8$ per cent. Our results demonstrate the power of high-resolution VLBI observations to provide strong constraints on the inner density profiles of lens galaxies., Comment: 24 pages, 15 figures, 3 tables. Accepted in MNRAS

The radio quasar luminosity function exhibits an upturn around $L_{6\rm\:GHz}=10^{23}$ W Hz$^{-1}$ that is well-modelled by a star-forming host galaxy population. This distribution leads some authors to cite star formation as the main radio emission mechanism in so-called radio-quiet quasars (RQQs). Understanding the origin of RQQ radio emission is crucial for our understanding of quasar feedback mechanisms -- responsible for the regulation of star-formation in the host galaxy -- and for understanding galaxy evolution as a whole. By observing RQQs that have been magnified by strong gravitational lensing, we have direct access to the RQQ population out to cosmic noon, where evidence for twin mini-jets has recently been found in a sub-\textmu Jy RQQ. Here we present radio observations of two lensed RQQs using the VLA at 5~GHz, the latest objects to be observed in a sample of quadruply-imaged RQQs above -30$^{\circ}$. In SDSS~J1004+4112 we find strong evidence for AGN-related radio emission in the variability of the source. In PG~1115+080 we find tentative evidence for AGN-related emission, determined by comparing the radio luminosity with modelled dust components. If confirmed in the case of PG~1115+080, which lies on the radio--FIR correlation, the result would reinforce the need for caution when applying the correlation to rule out jet activity and when assuming no AGN heating of FIR-emitting dust when calculating star formation rates. Our programme so far has shown that two of the faintest radio sources ever imaged show strong evidence for AGN-dominated radio emission., Comment: Accepted by MNRAS; 15 pages, 9 figures

We present DECORAS, a deep learning based approach to detect both point and extended sources from Very Long Baseline Interferometry (VLBI) observations. Our approach is based on an encoder-decoder neural network architecture that uses a low number of convolutional layers to provide a scalable solution for source detection. In addition, DECORAS performs source characterization in terms of the position, effective radius and peak brightness of the detected sources. We have trained and tested the network with images that are based on realistic Very Long Baseline Array (VLBA) observations at 20 cm. Also, these images have not gone through any prior de-convolution step and are directly related to the visibility data via a Fourier transform. We find that the source catalog generated by DECORAS has a better overall completeness and purity, when compared to a traditional source detection algorithm. DECORAS is complete at the 7.5$\sigma$ level, and has an almost factor of two improvement in reliability at 5.5$\sigma$. We find that DECORAS can recover the position of the detected sources to within 0.61 $\pm$ 0.69 mas, and the effective radius and peak surface brightness are recovered to within 20 per cent for 98 and 94 per cent of the sources, respectively. Overall, we find that DECORAS provides a reliable source detection and characterization solution for future wide-field VLBI surveys., Comment: submitted to MNRAS

Dual-Active Galactic Nuclei (AGN) are a natural consequence of the hierarchical structure formation scenario, and can provide an important test of various models for black hole growth. However, due to their rarity and difficulty to find at high redshift, very few confirmed dual-AGN are known at the epoch where galaxy formation peaks. Here we report the discovery of a gravitationally lensed dual-AGN system at redshift 2.37 comprising two optical/IR quasars separated by 6.5+/-0.6 kpc, and a third compact (R_eff = 0.45+/-0.02 kpc) red galaxy that is offset from one of the quasars by 1.7+/-0.1 kpc. From Very Large Array imaging at 3 GHz, we detect 600 and 340 pc-scale radio emission that is associated with both quasars. The 1.4 GHz luminosity densities of the radio sources are about 10^24.35 W / Hz, which is consistent with weak jets. However, the low brightness temperature of the emission is also consistent with star-formation at the level of 850 to 1150 M_sun / yr. Although this supports the scenario where the AGN and/or star-formation is being triggered through an ongoing triple-merger, a post-merger scenario where two black holes are recoiling is also possible, given that neither has a detected host galaxy., Comment: 5 pages, 3 figures, accepted for publication in MNRAS Letters

We present Low-Frequency Array (LOFAR) telescope observations of the radio-loud gravitational lens systems MG 0751+2716 and CLASS B1600+434. These observations produce images at 300 milliarcseconds (mas) resolution at 150 MHz. In the case of MG 0751+2716, lens modelling is used to derive a size estimate of around 2 kpc for the low-frequency source, which is consistent with a previous 27.4 GHz study in the radio continuum with Karl G. Jansky Very Large Array (VLA). This consistency implies that the low-frequency radio source is cospatial with the core-jet structure that forms the radio structure at higher frequencies, and no significant lobe emission or further components associated with star formation are detected within the magnified region of the lens. CLASS B1600+434 is a two-image lens where one of the images passes through the edge-on spiral lensing galaxy, and the low radio frequency allows us to derive limits on propagation effects, namely scattering, in the lensing galaxy. The observed flux density ratio of the two lensed images is 1.19 +/- 0.04 at an observed frequency of 150 MHz. The widths of the two images give an upper limit of 0.035 kpc m^-20/3 on the integrated scattering column through the galaxy at a distance approximately 1 kpc above its plane, under the assumption that image A is not affected by scattering. This is relatively small compared to limits derived through very long baseline interferometry (VLBI) studies of differential scattering in lens systems. These observations demonstrate that LOFAR is an excellent instrument for studying gravitational lenses. We also report on the inability to calibrate three further lens observations: two from early observations that have less well determined station calibration, and a third observation impacted by phase transfer problems., Comment: Accepted to a special issue of A&A on sub-arcsecond imaging with LOFAR

[abridged] The International LOFAR Telescope is an interferometer with stations spread across Europe. With baselines of up to ~2,000 km, LOFAR has the unique capability of achieving sub-arcsecond resolution at frequencies below 200 MHz, although this is technically and logistically challenging. Here we present a calibration strategy that builds on previous high-resolution work with LOFAR. We give an overview of the calibration strategy and discuss the special challenges inherent to enacting high-resolution imaging with LOFAR, and describe the pipeline, which is publicly available, in detail. We demonstrate the calibration strategy by using the pipeline on P205+55, a typical LOFAR Two-metre Sky Survey (LoTSS) pointing. We perform in-field delay calibration, solution referencing to other calibrators, self-calibration, and imaging of example directions of interest in the field. For this specific field and these ionospheric conditions, dispersive delay solutions can be transferred between calibrators up to ~1.5 degrees away, while phase solution transferral works well over 1 degree. We demonstrate a check of the astrometry and flux density scale. Imaging in 17 directions, the restoring beam is typically 0.3" x 0.2" although this varies slightly over the entire 5 square degree field of view. We achieve ~80 to 300 $\mu$Jy/bm image rms noise, which is dependent on the distance from the phase centre; typical values are ~90 $\mu$Jy/bm for the 8 hour observation with 48 MHz of bandwidth. Seventy percent of processed sources are detected, and from this we estimate that we should be able to image ~900 sources per LoTSS pointing. This equates to ~3 million sources in the northern sky, which LoTSS will entirely cover in the next several years. Future optimisation of the calibration strategy for efficient post-processing of LoTSS at high resolution (LoTSS-HR) makes this estimate a lower limit., Comment: Accepted to a special issue of A&A on sub-arcsecond imaging with LOFAR. 24 pages, 16 figures

The Low-Frequency Array (LOFAR) Long-Baseline Calibrator Survey (LBCS) was conducted between 2014 and 2019 in order to obtain a set of suitable calibrators for the LOFAR array. In this paper we present the complete survey, building on the preliminary analysis published in 2016 which covered approximately half the survey area. The final catalogue consists of 30006 observations of 24713 sources in the northern sky, selected for a combination of high low-frequency radio flux density and flat spectral index using existing surveys (WENSS, NVSS, VLSS, and MSSS). Approximately one calibrator per square degree, suitable for calibration of $\geq$ 200 km baselines is identified by the detection of compact flux density, for declinations north of 30 degrees and away from the Galactic plane, with a considerably lower density south of this point due to relative difficulty in selecting flat-spectrum candidate sources in this area of the sky. Use of the VLBA calibrator list, together with statistical arguments by comparison with flux densities from lower-resolution catalogues, allow us to establish a rough flux density scale for the LBCS observations, so that LBCS statistics can be used to estimate compact flux densities on scales between 300 mas and 2 arcsec, for sources observed in the survey. The LBCS can be used to assess the structures of point sources in lower-resolution surveys, with significant reductions in the degree of coherence in these sources on scales between 2 arcsec and 300 mas. The LBCS survey sources show a greater incidence of compact flux density in quasars than in radio galaxies, consistent with unified schemes of radio sources. Comparison with samples of sources from interplanetary scintillation (IPS) studies with the Murchison Widefield Array (MWA) shows consistent patterns of detection of compact structure in sources observed both interferometrically with LOFAR and using IPS., Comment: Accepted to a special issue of A&A on sub-arcsecond imaging with LOFAR

The low-frequency radio spectra of the hotspots within powerful radio galaxies can provide valuable information about the physical processes operating at the site of the jet termination. These processes are responsible for the dissipation of jet kinetic energy, particle acceleration, and magnetic-field generation. Here we report new observations of the powerful radio galaxy Cygnus A using the Low Frequency Array (LOFAR) between 109 and 183 MHz, at an angular resolution of ~3.5 arcsec. The radio emission of the lobes is found to have a complex spectral index distribution, with a spectral steepening found towards the centre of the source. For the first time, a turnover in the radio spectrum of the two main hotspots of Cygnus A has been directly observed. By combining our LOFAR imaging with data from the Very Large Array at higher frequencies, we show that the very rapid turnover in the hotspot spectra cannot be explained by a low-energy cut-off in the electron energy distribution, as has been previously suggested. Thermal (free-free) absorption or synchrotron self absorption models are able to describe the low-frequency spectral shape of the hotspots, however, as with previous studies, we find that the implied model parameters are unlikely, and interpreting the spectra of the hotspots remains problematic., Comment: 9 pages, 3 figures, accepted for publication in MNRAS

We report extremely faint 144 MHz radio emission from two gravitationally lensed quasars, SDSS J1004+4112 (z = 1.730) and SDSS J2222+2745 (z = 2.803), using the LOFAR Two Metre Sky Survey (LoTSS) data release 2. After correcting for the lensing magnifications, the two objects have intrinsic flux-densities of 13+/-2 and 58+/-6 uJy, respectively, corresponding to 144 MHz rest-frame luminosities of 10^(23.2+/-0.2) and 10^(24.42+/-0.05) W / Hz, respectively. In the case of SDSS J1004+4112, the intrinsic flux density is close to the confusion limit of LoTSS, making this radio source the faintest to be detected thus far at low frequencies, and the lowest luminosity known at z > 0.65. Under the assumption that all of the radio emission is due to star-formation processes, the quasar host galaxies are predicted to have star-formation rates of 5.5^(+1.8)_(-1.4) and 73^(+34)_(-22) M / yr, respectively. Further multi-wavelength observations at higher angular resolution will be needed to determine if any of the detected radio emission is due to weak jets associated with the quasars., Comment: 5 pages, 1 figure, accepted for publication in MNRAS Letters

The extreme astrophysical processes and conditions that characterize the early Universe are expected to result in young galaxies that are dynamically different from those observed today. This is because the strong effects associated with galaxy mergers and supernova explosions would lead to most young star-forming galaxies being dynamically hot, chaotic and strongly unstable. Here we report the presence of a dynamically cold, but highly star-forming, rotating disk in a galaxy at redshift ($z$) 4.2, when the Universe was just 1.4 billion years old. Galaxy SPT-S J041839-4751.9 is strongly gravitationally lensed by a foreground galaxy at $z = 0.263$, and it is a typical dusty starburst, with global star-forming and dust properties that are in agreement with current numerical simulations and observations of its galaxy population. Interferometric imaging at a spatial resolution of about 60 pc reveals a ratio of rotational-to-random motions of $V/\sigma = 9.7\pm 0.4$, which is at least four times larger than expected from any galaxy evolution model at this epoch, but similar to the ratios of spiral galaxies in the local Universe. We derive a rotation curve with the typical shape of nearby massive spiral galaxies, which demonstrates that at least some young galaxies are dynamically akin to those observed in the local Universe, and only weakly affected by extreme physical processes., Comment: Published in Nature on 12 August 2020. The published version is available at http://www.nature.com/articles/s41586-020-2572-6

We resolve the host galaxies of seven gravitationally lensed quasars at redshift 1.5 to 2.8 using observations with the Atacama Large (sub-)Millimetre Array. Using a visibility-plane lens modelling technique, we create pixellated reconstructions of the dust morphology, and CO line morphology and kinematics. We find that the quasar hosts in our sample can be distinguished into two types: 1) galaxies characterised by clumpy, extended dust distributions ($R_{\rm eff}\sim2$ kpc) and mean star formation rate surface densities comparable to sub-mm-selected dusty star-forming galaxies ($\Sigma_{\rm SFR}\sim3$ M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$); 2) galaxies that have sizes in dust emission similar to coeval passive galaxies and compact starbursts ($R_{\rm eff}\sim0.5$ kpc), with high mean star formation rate surface densities ($\Sigma_{\rm SFR}=$ 400$-$4500 M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$) that may be Eddington-limited or super-Eddington. The small size of some quasar hosts suggests that we observe them at a stage in their transformation into compact spheroids, where a high density of dynamically unstable gas leads to efficient star formation and black hole accretion. For the one system where we probe the mass of the gas reservoir, we find a gas fraction of just $0.06 \pm 0.04$ and a depletion timescale of $50 \pm 40$ Myr, suggesting it is transitioning into quiescence. In general, we expect that the extreme level of star formation in the compact quasar host galaxies will rapidly exhaust their gas reservoirs and could quench with or without help from active galactic nuclei feedback., Comment: Accepted by MNRAS; 22 pages

The four persistent radio sources in the northern sky with the highest flux density at metre wavelengths are Cassiopeia A, Cygnus A, Taurus A, and Virgo A; collectively they are called the A-team. Their flux densities at ultra-low frequencies (<100 MHz) can reach several thousands of janskys, and they often contaminate observations of the low-frequency sky by interfering with image processing. Furthermore, these sources are foreground objects for all-sky observations hampering the study of faint signals, such as the cosmological 21 cm line from the epoch of reionisation. We aim to produce robust models for the surface brightness emission as a function of frequency for the A-team sources at ultra-low frequencies. These models are needed for the calibration and imaging of wide-area surveys of the sky with low-frequency interferometers. This requires obtaining images at an angular resolution better than 15 arcsec with a high dynamic range and good image fidelity. We observed the A-team with the Low Frequency Array (LOFAR) at frequencies between 30 MHz and 77 MHz using the Low Band Antenna (LBA) system. We reduced the datasets and obtained an image for each A-team source. The paper presents the best models to date for the sources Cassiopeia A, Cygnus A, Taurus A, and Virgo A between 30 MHz and 77 MHz. We were able to obtain the aimed resolution and dynamic range in all cases. Owing to its compactness and complexity, observations with the long baselines of the International LOFAR Telescope will be required to improve the source model for Cygnus A further., Comment: 7 pages, 2 figures, accepted A&A, online data on A&A website

We present an analysis of observations with the Atacama Large (sub-)Millimetre Array (ALMA) of the 380 GHz water emission line and CO (11-10) emission line from MG J0414+0534, a gravitationally lensed dusty star-forming galaxy that hosts a type 1 quasar. We also present observations at 1.6 GHz with global very long baseline interferometry (VLBI) of the radio source. We confirm the previously reported detection of the 380 GHz water line and that the flux density ratio between the two merging lensed images is reversed with respect to the radio/mm continuum. We further find tentative evidence of variability in the integrated line intensity on timescales of days. We show that the 380 GHz water line has two components of emission: a disc that is around 35 pc in diameter around the quasar, and another component of emission that is offset $\sim$600 pc perpendicular to the disc that lies close to the lensing caustic. With lens modelling of the multi-wavelength data sets, we construct a model for the quasar system consisting of a circumnuclear disc of molecular gas with a size of about 60 pc bisected by radio jets extending to a distance of about 200 pc from the radio core. Our findings suggest that observations with ALMA of high-excitation molecular lines from strongly lensed quasars could allow detailed studies of AGN accretion and feedback at the cosmic peak of black hole and galaxy growth., Comment: Accepted by MNRAS. 11 pages, 9 figures

We present a sub-kpc resolved study of the interstellar medium properties in SDP.81, a z=3.042 strongly gravitationally lensed dusty star-forming galaxy, based on high-resolution, multi-band ALMA observations of the FIR continuum, CO ladder and the [CII] line. Using a visibility-plane lens modelling code, we achieve a median source-plane resolution of ~200 pc. We use photon-dominated region (PDR) models to infer the physical conditions - far-UV field strength, density, and PDR surface temperature - of the star-forming gas on 200-pc scales, finding a FUV field strength of ~10^3-10^4 G0, gas density of ~10^5 cm^-3 and cloud surface temperatures up to 1500 K, similar to those in the Orion Trapezium region. The [CII] emission is significantly more extended than that FIR continuum: ~50 per cent of [CII] emission arises outside the FIR-bright region. The resolved [CII]/FIR ratio varies by almost 2 dex across the source, down to ~2x10^-4 in the star-forming clumps. The observed [CII]/FIR deficit trend is consistent with thermal saturation of the C+ fine-structure level occupancy at high gas temperatures. We make the source-plane reconstructions of all emission lines and continuum data publicly available., Comment: 27 pages, 22 figures. 24 MB in total. MNRAS accepted 2020 March 17. Data available on request

In this paper, we exploit the gravitational lensing effect to detect proper motion in the highly magnified gravitationally lensed source MG B2016+112. We find positional shifts up to 6 mas in the lensed images by comparing two Very Long Baseline Interferometric (VLBI) radio observations at 1.7 GHz that are separated by 14.359 years, and provide an astrometric accuracy of the order of tens of $\mu$as. From lens modelling, we exclude a shift in the lensing galaxy as the cause of the positional change of the lensed images, and we assign it to the background source. The source consists of four sub-components separated by $\sim 175$ pc, with proper motion of the order of tens $\mu$as yr$^{-1}$ for the two components at highest magnification ($\mu\sim350$) and of the order of a few mas yr$^{-1}$ for the two components at lower magnification ($\mu\sim2$). We propose single AGN and dual AGN scenarios to explain the source plane. Although, the latter interpretation is supported by the archival multi-wavelength properties of the object. In this case, MG B2016+112 would represent the highest redshift dual radio-loud AGN system discovered thus far, and would support the merger interpretation for such systems. Also, given the low probability ($\sim10^{-5}$) of detecting a dual AGN system that is also gravitationally lensed, if confirmed, this would suggest that such dual AGN systems must be more abundant in the early Universe than currently thought., Comment: 11 pages, 4 figures; accepted for publication by Astronomy and Astrophysics

We present a study of the stellar host galaxy, CO (1$-$0) molecular gas distribution and AGN emission on 50 to 500 pc-scales of the gravitationally lensed dust-obscured AGN MG J0751+2716 and JVAS B1938+666 at redshifts 3.200 and 2.059, respectively. By correcting for the lensing distortion using a grid-based lens modelling technique, we spatially locate the different emitting regions in the source plane for the first time. Both AGN host galaxies have 300 to 500 pc-scale size and surface brightness consistent with a bulge/pseudo-bulge, and 2 kpc-scale AGN radio jets that are embedded in extended molecular gas reservoirs that are 5 to 20 kpc in size. The CO (1$-$0) velocity fields show structures possibly associated with discs (elongated velocity gradients) and interacting objects (off-axis velocity components). There is evidence for a decrement in the CO (1$-$0) surface brightness at the location of the host galaxy, which may indicate radiative feedback from the AGN, or offset star formation.We find CO-H$_2$ conversion factors of around $\alpha_{\rm CO} = 1.5\pm0.5$ (K km s$^{-1}$ pc$^2$)$^{-1}$, molecular gas masses of $> 3\times10^{10}$ M$_{\odot}$, dynamical masses of $\sim 10^{11}$ M$_{\odot}$ and gas fractions of around 60 per cent. The intrinsic CO line luminosities are comparable to those of unobscured AGN and dusty star-forming galaxies at similar redshifts, but the infrared luminosities are lower, suggesting that the targets are less efficient at forming stars. Therefore, they may belong to the AGN feedback phase predicted by galaxy formation models, because they are not efficiently forming stars considering their large amount of molecular gas., Comment: 22 pages, 15 figures with additional 1 figure in the Appendix. Accepted for publication by MNRAS

Gravitational lensing is a powerful tool for quantifying the mass content and distribution in distant galaxies. By using milliarcsecond angular resolution observations of radio-loud gravitationally lensed sources it is also possible to detect and quantify small deviations from a smooth mass density distribution, which can be due to low mass substructures in the lensing galaxy. We present high-resolution global VLBI observations of the gravitationally lensed radio source MG J0751+2716 (at z = 3.2), that shows evidence of both compact and extended structure (core-jet morphology) across several gravitational arcs. These data provide a wealth of observational constraints that are used to determine the inner (baryonic and dark matter) mass profile of a group of galaxies and also investigate the smoothness of the dark matter distribution on mas-scales, which is sensitive to possible structures of $10^{6-7}$ M$_{\odot}$ within the lensing halo or along the line-of-sight. Our lens modelling finds evidence for astrometric anomalies in this system, which suggest presence of extra mass structure in the lens model. To date this kind of detailed studies of gravitational lensing systems like MG J0751+2716 has been limited by the currently small sample of radio-loud gravitational lenses. In this context, we also present a new pilot gravitational lens search in the VLBI survey mJIVE-20, in perspective of future surveys with the next generation of radio interferometers., Comment: 7 pages, 3 figures; accepted for the proceedings of the 14th European VLBI Network Symposium & Users Meeting, held on 8-11 October 2018 in Granada, Spain. It will appear in Proceedings of Science PoS(EVN2018)033

We present a novel pilot search for gravitational lenses in the mJIVE-20 survey, which observed $24\,903$ radio sources selected from FIRST with the VLBA at an angular resolution of 5 mas. We have taken the visibility data for an initial $3\,640$ sources that were detected by the mJIVE-20 observations and re-mapped them to make wide-field images, selecting fourteen sources that had multiple components separated by $\geq100$ mas, with a flux-ratio of $\leq15$:$1$ and a surface brightness consistent with gravitational lensing. Two of these candidates are re-discoveries of gravitational lenses found as part of CLASS. The remaining twelve candidates were then re-observed at 1.4 GHz and then simultaneously at 4.1 and 7.1 GHz with the VLBA to measure the spectral index and surface brightness of the individual components as a function of frequency. Ten were rejected as core-jet or core-hotspot(s) systems, with surface brightness distributions and/or spectral indices inconsistent with gravitational lensing, and one was rejected after lens modelling demonstrated that the candidate lensed images failed the parity test. The final lens candidate has an image configuration that is consistent with a simple lens mass model, although further observations are required to confirm the lensing nature. Given the two confirmed gravitational lenses in the mJIVE-20 sample, we find a robust lensing-rate of $1$:($318\pm225$) for a statistical sample of 635 radio sources detected on mas-scales, which is consistent with that found for CLASS., Comment: 31 pages, 22 figures; accepted for publication in MNRAS

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., Comment: 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

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168 MHz survey of the Northern sky with diverse and ambitious science goals. Many of the scientific objectives of LoTSS rely upon, or are enhanced by, the association or separation of the sometimes incorrectly catalogued radio components into distinct radio sources, and the identification and characterisation of the optical counterparts to these sources. Here we present the source associations and optical and/or IR identifications for sources in the first data release, which are made using a combination of statistical techniques and visual association and identification. We document in detail the colour- and magnitude-dependent likelihood ratio method used for statistical identification as well as the Zooniverse project, called LOFAR Galaxy Zoo, used for the visual classification. We describe the process used to select which of these two different methods is most appropriate for each LoTSS source. The final LoTSS-DR1-IDs value-added catalogue presented contains 318,520 radio sources, of which 231,716 (73%) have optical and/or IR identifications in Pan-STARRS and WISE. The value-added catalogue is available online at https://lofar-surveys.org/, as part of this data release., Comment: 8 figures, 5 tables and 23 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 at https://lofar-surveys.org/ upon publication in the journal

We use a sample of 17 strong gravitational lens systems from the BELLS GALLERY survey to quantify the amount of low-mass dark matter haloes within the lensing galaxies and along their lines of sight, and to constrain the properties of dark matter. Based on a detection criterion of 10$\sigma$, we report no significant detection in any of the lenses. Using the sensitivity function at the 10-$\sigma$ level, we have calculated the predicted number of detectable cold dark matter (CDM) line-of-sight haloes to be $\mu_{l} = 1.17\pm1.08$, in agreement with our null detection. Assuming a detection sensitivity that improved to the level implied by a 5-$\sigma$ threshold, the expected number of detectable line-of-sight haloes rises to $\mu_l = 9.0\pm3.0$. Whilst the current data find zero detections at this sensitivity level (which has a probability of P$^{{\rm5}\sigma}_{{\rm CDM}}(n_{\rm det}=0)$=0.0001 and would be in strong tension with the CDM framework), we find that such a low detection threshold leads to many spurious detections and non-detections and therefore the current lack of detections is unreliable and requires data with improved sensitivity. Combining this sample with a subsample of 11 SLACS lenses, we constrain the half-mode mass to be $\log$(M$_{\rm hm}) < 12.26$ at the 2-$\sigma$ level. The latter is consistent with resonantly produced sterile neutrino masses m$_{\rm s} < 0.8$ keV at any value of the lepton asymmetry at the 2-$\sigma$ level., Comment: 14 pages, 7 figures, submitted for publication on MNRAS

We present a study of seventeen LAEs at redshift 2$

We present an analysis of archival observations with the Atacama Large (sub-)Millimetre Array (ALMA) of the gravitationally lensed quasar MG J0414+0534, which show four compact images of the quasar and an Einstein ring from the dust associated with the quasar host galaxy. We confirm that the flux-ratio anomalies observed in the mid-infrared and radio persists into the sub-mm for the continuum images of the quasar. We report the detection of CO (11-10) spectral line emission, which traces a region of compact gas around the quasar nucleus. This line emission also shows evidence of a flux-ratio anomaly between the merging lensed images that is consistent with those observed at other wavelengths, suggesting high-excitation CO can also provide a useful probe of substructures that is unaffected by microlensing or dust extinction. However, we do not detect the candidate dusty dwarf galaxy that was previously reported with this dataset, which we conclude is due to a noise artefact. Thus, the cause of the flux-ratio anomaly between the merging lensed images is still unknown. The composite compact and diffuse emission in this system suggest lensed quasar-starbursts will make excellent targets for detecting dark sub-haloes and testing models for dark matter., Comment: Accepted as MNRAS Letter

We present milliarcsecond (mas) angular resolution observations of the gravitationally lensed radio source MG J0751+2716 (at z=3.2) obtained with global Very Long Baseline Interferometry (VLBI) at 1.65 GHz. The background object is highly resolved in the tangential and radial directions, showing evidence of both compact and extended structure across several gravitational arcs that are 200 to 600~mas in size. By identifying compact sub-components in the multiple images, we constrain the mass distribution of the foreground z=0.35 gravitational lens using analytic models for the main deflector [power-law elliptical mass model; $\rho(r) \propto r^{-\gamma}$, where $\gamma=2$ corresponds to isothermal] and for the members of the galaxy group. Moreover, our mass models with and without the group find an inner mass-density slope steeper than isothermal for the main lensing galaxy, with $\gamma_1 = 2.08 \pm 0.02$ and $\gamma_2 = 2.16 \pm 0.02$ at the 4.2$\sigma$ level and 6.8$\sigma$ level, respectively, at the Einstein radius ($b_1 = 0.4025 \pm 0.0008$ and $b_2 = 0.307 \pm 0.002$ arcsec, respectively). We find randomly distributed image position residuals of about 3 mas, which are much larger that the measurement errors ($40$ $\mu$as on average). This suggests that at the mas level, the assumption of a smooth mass distribution fails, requiring additional structure in the model. However, given the environment of the lensing galaxy, it is not clear whether this extra mass is in the form of sub-haloes within the lens or along the line of sight, or from a more complex halo for the galaxy group., Comment: 15 pages, 8 figures; published by MNRAS

Stringent observational constraints on the sub-galactic matter power spectrum would allow one to distinguish between the concordance $\Lambda$CDM and the various alternative dark-matter models that predict significantly different properties of mass structure in galactic haloes. Galaxy-galaxy strong gravitational lensing provides a unique opportunity to probe the sub-galactic mass structure in lens galaxies beyond the Local Group. Here, we demonstrate the first application of a novel methodology to observationally constrain the sub-galactic matter power spectrum in the inner regions of massive elliptical lens galaxies on 1-10 kpc scales from the power spectrum of surface-brightness anomalies in highly magnified galaxy-scale Einstein rings and gravitational arcs. The pilot application of our approach to Hubble Space Telescope (HST/WFC3/F390W) observations of the SLACS lens system SDSS J0252+0039 allows us to place the following observational constraints (at the 99 per cent confidence level) on the dimensionless convergence power spectrum $\Delta^{2}_{\delta\kappa}$ and the standard deviation in the aperture mass $\sigma_{\rm AM}$: $\Delta^{2}_{\delta\kappa}<1$ ($\sigma_{\rm AM}< 0.8 \times 10^8 M_\odot$) on 0.5-kpc scale, $\Delta^{2}_{\delta\kappa}<0.1$ ($\sigma_{\rm AM}< 1 \times 10^8 M_\odot$) on 1-kpc scale and $\Delta^{2}_{\delta\kappa}<0.01$ ($\sigma_{\rm AM}< 3 \times 10^8 M_\odot$) on 3-kpc scale. These first upper-limit constraints still considerably exceed the estimated effect of CDM subhaloes. However, future analysis of a larger sample of galaxy-galaxy strong lens systems can substantially narrow down these limits and possibly rule out dark-matter models that predict a significantly higher level of density fluctuations on the critical sub-galactic scales., Comment: 16 pages, 10 figures, updated and resubmitted to MNRAS together with an accompanying methodology paper in response to reviewer's comments

We use new and archival Chandra observations of Cygnus A, totalling $\sim$1.9 Ms, to investigate the distribution and temperature structure of gas lying within the projected extent of the cocoon shock and exhibiting a rib-like structure. We confirm that the X-rays are dominated by thermal emission with an average temperature of around 4 keV, and have discovered an asymmetry in the temperature gradient, with the southwestern part of the gas cooler than the rest by up to 2 keV. Pressure estimates suggest that the gas is a coherent structure of single origin located inside the cocoon, with a mass of roughly $2\times10^{10} M_{\odot}$. We conclude that the gas is debris resulting from disintegration of the cool core of the Cygnus A cluster after the passage of the jet during the early stages of the current epoch of activity. The 4 keV gas now lies on the central inside surface of the hotter cocoon rim. The temperature gradient could result from an offset between the centre of the cluster core and the Cygnus A host galaxy at the switch-on of current radio activity., Comment: 15 pages, 12 figures, accepted by MNRAS

We present a deep, low-frequency radio continuum study of the nearby Fanaroff--Riley class I (FR I) radio galaxy 3C 31 using a combination of LOw Frequency ARray (LOFAR; 30--85 and 115--178 MHz), Very Large Array (VLA; 290--420 MHz), Westerbork Synthesis Radio Telescope (WSRT; 609 MHz) and Giant Metre Radio Telescope (GMRT; 615 MHz) observations. Our new LOFAR 145-MHz map shows that 3C 31 has a largest physical size of $1.1$ Mpc in projection, which means 3C 31 now falls in the class of giant radio galaxies. We model the radio continuum intensities with advective cosmic-ray transport, evolving the cosmic-ray electron population and magnetic field strength in the tails as functions of distance to the nucleus. We find that if there is no in-situ particle acceleration in the tails, then decelerating flows are required that depend on radius $r$ as $v\propto r^{\beta}$ ($\beta\approx -1$). This then compensates for the strong adiabatic losses due to the lateral expansion of the tails. We are able to find self-consistent solutions in agreement with the entrainment model of Croston & Hardcastle, where the magnetic field provides $\approx$$1/3$ of the pressure needed for equilibrium with the surrounding intra-cluster medium (ICM). We obtain an advective time-scale of $\approx$$190$ Myr, which, if equated to the source age, would require an average expansion Mach number ${\cal M} \approx 5$ over the source lifetime. Dynamical arguments suggest that instead, either the outer tail material does not represent the oldest jet plasma or else the particle ages are underestimated due to the effects of particle acceleration on large scales., Comment: MNRAS accepted, 17 pages, 9 figures, 2 tables + 3 pages appendix with 3 figures and 2 tables

Context. The Sun is an active star that produces large-scale energetic events such as solar flares and coronal mass ejections and numerous smaller-scale events such as solar jets. These events are often associated with accelerated particles that can cause emission at radio wavelengths. The reconfiguration of the solar magnetic field in the corona is believed to be the cause of the majority of solar energetic events and accelerated particles. Aims. Here, we investigate a bright J-burst that was associated with a solar jet and the possible emission mechanism causing these two phenomena. Methods. We used data from the Solar Dynamics Observatory (SDO) to observe a solar jet, and radio data from the Low Frequency Array (LOFAR) and the Nan\c{c}ay Radioheliograph (NRH) to observe a J-burst over a broad frequency range (33-173 MHz) on 9 July 2013 at ~11:06 UT. Results. The J-burst showed fundamental and harmonic components and it was associated with a solar jet observed at extreme ultraviolet wavelengths with SDO. The solar jet occurred at a time and location coincident with the radio burst, in the northern hemisphere, and not inside a group of complex active regions in the southern hemisphere. The jet occurred in the negative polarity region of an area of bipolar plage. Newly emerged positive flux in this region appeared to be the trigger of the jet. Conclusions. Magnetic reconnection between the overlying coronal field lines and the newly emerged positive field lines is most likely the cause of the solar jet. Radio imaging provides a clear association between the jet and the J-burst which shows the path of the accelerated electrons., Comment: 11 pages, 8 figures

We present a new analysis of the widely used relation between cavity power and radio luminosity in clusters of galaxies with evidence for strong AGN feedback. We study the correlation at low radio frequencies using two new surveys - the First Alternative Data Release of the TIFR GMRT Sky Survey (TGSS ADR1) at 148 MHz and LOFAR's first all-sky survey, the Multifrequency Snapshot Sky Survey (MSSS) at 140 MHz. We find a scaling relation $P_{\rm cav} \propto L_{148}^{\beta}$, with a logarithmic slope of $\beta = 0.51 \pm 0.14$, which is in good agreement with previous results based on data at 327 MHz. The large scatter present in this correlation confirms the conclusion reached at higher frequencies that the total radio luminosity at a single frequency is a poor predictor of the total jet power. We show that including measurements at 148 MHz alone is insufficient to reliably compute the bolometric radio luminosity and reduce the scatter in the correlation. For a subset of four well-resolved sources, we examine the detected extended structures at low frequencies and compare with the morphology known from higher frequency images and Chandra X-ray maps. In Perseus we discuss details in the structures of the radio mini-halo, while in the 2A 0335+096 cluster we observe new diffuse emission associated with multiple X-ray cavities and likely originating from past activity. For A2199 and MS 0735.6+7421, we confirm that the observed low-frequency radio lobes are confined to the extents known from higher frequencies. This new low-frequency analysis highlights the fact that existing cavity power to radio luminosity relations are based on a relatively narrow range of AGN outburst ages. We discuss how the correlation could be extended using low frequency data from the LOFAR Two-metre Sky Survey (LoTSS) in combination with future, complementary deeper X-ray observations.

We have observed 104 gravitationally-lensed quasars at $z\sim1-4$ with Herschel/SPIRE, the largest such sample ever studied. By targeting gravitational lenses, we probe intrinsic far-infrared (FIR) luminosities and star formation rates (SFRs) more typical of the population than the extremely luminous sources that are otherwise accessible. We detect 72 objects with Herschel/SPIRE and find 66 percent (69 sources) of the sample have spectral energy distributions (SEDs) characteristic of dust emission. For 53 objects with sufficiently constrained SEDs, we find a median effective dust temperature of $38^{+12}_{-5}$ K. By applying the radio-infrared correlation, we find no evidence for an FIR excess which is consistent with star-formation-heated dust. We derive a median magnification-corrected FIR luminosity of $3.6^{+4.8}_{-2.4}~\times 10^{11}~{\rm L_{\odot}}$ and median SFR of $120^{+160}_{-80}~{\rm M_{\odot}~yr^{-1}}$ for 94 quasars with redshifts. We find $\sim10$ percent of our sample have FIR properties similar to typical dusty star-forming galaxies at $z\sim2-3$ and a range of SFRs $<20-10000~{\rm M_{\odot}~yr^{-1}}$ for our sample as a whole. These results are in line with current models of quasar evolution and suggests a coexistence of dust-obscured star formation and AGN activity is typical of most quasars. We do not find a statistically-significant difference in the FIR luminosities of quasars in our sample with a radio excess relative to the radio-infrared correlation. Synchrotron emission is found to dominate at FIR wavelengths for $<15$ percent of those sources classified as powerful radio galaxies., Comment: 47 pages, 89 figures, accepted for publication in MNRAS