Showing 1–50 of 145 results for author: McKean, J

First joint oscillation analysis of Super-Kamiokande atmospheric and T2K accelerator neutrino data

Authors: Super-Kamiokande, T2K collaborations, :, S. Abe, K. Abe, N. Akhlaq, R. Akutsu, H. Alarakia-Charles, A. Ali, Y. I. Alj Hakim, S. Alonso Monsalve, S. Amanai, C. Andreopoulos, L. H. V. Anthony, M. Antonova, S. Aoki, K. A. Apte, T. Arai, T. Arihara, S. Arimoto, Y. Asada, R. Asaka, Y. Ashida, E. T. Atkin, N. Babu , et al. (524 additional authors not shown)

Abstract: 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… ▽ More 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$σ$ exclusion of CP-conservation (defined as $J_{CP}=0$) and a preference for the normal mass ordering. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 10 pages, 3 figures

Radio imaging of gravitationally lensed radio-quiet quasars

Authors: Neal Jackson, Shruti Badole, Thomas Dugdale, Hannah R. Stacey, Philippa Hartley, J. P. McKean, .

Abstract: 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 r… ▽ More 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. △ Less

Submitted 4 April, 2024; v1 submitted 28 March, 2024; originally announced March 2024.

Comments: Accepted by MNRAS

Complex angular structure of three elliptical galaxies from high-resolution ALMA observations of strong gravitational lenses

Authors: H. R. Stacey, D. M. Powell, S. Vegetti, J. P. McKean, C. D. Fassnacht, D. Wen, C. M. O'Riordan

Abstract: 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-)Millimetr… ▽ More 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. △ Less

Submitted 27 May, 2024; v1 submitted 7 March, 2024; originally announced March 2024.

Comments: In print with A&A

At the end of cosmic noon: Short gas depletion times in unobscured quasars at $z \sim$ 1

Authors: M. Frias Castillo, M. Rybak, J. Hodge, P. van der Werf, L. J. Abbo, F. J. Ballieux, S. Ward, C. Harrison, G. Calistro Rivera, J. P. McKean, H. R. Stacey

Abstract: 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 pre… ▽ More 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 (λ_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. △ Less

Submitted 12 January, 2024; v1 submitted 9 January, 2024; originally announced January 2024.

Comments: 20 pages. Accepted for publication in A&A

Measurements of the $ν_μ$ and $\barν_μ$-induced Coherent Charged Pion Production Cross Sections on $^{12}C$ by the T2K experiment

Authors: K. Abe, N. Akhlaq, R. Akutsu, A. Ali, S. Alonso Monsalve, C. Alt, C. Andreopoulos, M. Antonova, S. Aoki, T. Arihara, Y. Asada, Y. Ashida, E. T. Atkin, M. Barbi, G. J. Barker, G. Barr, D. Barrow, M. Batkiewicz-Kwasniak, V. Berardi, L. Berns, S. Bhadra, A. Blanchet, A. Blondel, S. Bolognesi, T. Bonus , et al. (359 additional authors not shown)

Abstract: We report an updated measurement of the $ν_μ$-induced, and the first measurement of the $\barν_μ$-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_{μ,π} > 0.2$ GeV, $\cos(θ_μ) > 0.8$ and $\cos(θ_π) > 0.6$, and at a mean (anti)neutrino energy of 0.85 GeV using the T2K… ▽ More We report an updated measurement of the $ν_μ$-induced, and the first measurement of the $\barν_μ$-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_{μ,π} > 0.2$ GeV, $\cos(θ_μ) > 0.8$ and $\cos(θ_π) > 0.6$, and at a mean (anti)neutrino energy of 0.85 GeV using the T2K near detector. The measured $ν_μ$ 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ν_μ$-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. △ Less

Submitted 14 October, 2023; v1 submitted 31 August, 2023; originally announced August 2023.

Journal ref: Phys.Rev.D 108 (2023) 9, 092009

On the source counts of VLBI-detected radio sources and the prospects of all-sky surveys with current and next generation instruments

Authors: S. Rezaei, J. P. McKean, A. T. Deller, J. F. Radcliffe

Abstract: 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 aro… ▽ More 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 $η_{\rm VLBI} = -1.74\pm 0.02$, which is shallower than in the cases of the FIRST parent population ($η_{\rm FIRST} = -1.77\pm 0.02$) and for compact radio sources selected at higher frequencies ($η_{\rm JBF} = -2.06\pm 0.02$). From this, we find that all-sky ($3π~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. △ Less

Submitted 30 August, 2023; originally announced August 2023.

PRUSSIC II -- ALMA imaging of dense-gas tracers in SDP.81: Evidence for low mechanical heating and a sub-solar metallicity in a z=3.04 dusty galaxy

Authors: M. Rybak, J. van Marrewijk, J. A. Hodge, P. Andreani, G. Calistro Rivera, L. Graziani, J. P. McKean, S. Viti, P. P. van der Werf

Abstract: 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… ▽ More 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. △ Less

Submitted 25 September, 2023; v1 submitted 5 August, 2023; originally announced August 2023.

Comments: A&A accepted, corrected version

Journal ref: A&A 679, A119 (2023)

Strong gravitational lensing as a probe of dark matter

Authors: S. Vegetti, S. Birrer, G. Despali, C. D. Fassnacht, D. Gilman, Y. Hezaveh, L. Perreault Levasseur, J. P. McKean, D. M. Powell, C. M. O'Riordan, G. Vernardos

Abstract: 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-sca… ▽ More 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. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Updated T2K measurements of muon neutrino and antineutrino disappearance using 3.6 $\times$ 10$^{21}$ protons on target

Authors: K. Abe, N. Akhlaq, R. Akutsu, H. Alarakia-Charles, A. Ali, Y. I. Alj Hakim, S. Alonso Monsalve, C. Alt, C. Andreopoulos, M. Antonova, S. Aoki, T. Arihara, Y. Asada, Y. Ashida, E. T. Atkin, M. Barbi, G. J. Barker, G. Barr, D. Barrow, M. Batkiewicz-Kwasniak, F. Bench, V. Berardi, L. Berns, S. Bhadra, A. Blanchet , et al. (385 additional authors not shown)

Abstract: 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} θ_{23}$ and $Δm_{32}^2$ independently for neutrinos and antineutrinos. The aforementioned symmetry violation would manifest as an inconsis… ▽ More 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} θ_{23}$ and $Δ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. △ Less

Submitted 16 October, 2023; v1 submitted 16 May, 2023; originally announced May 2023.

First measurement of muon neutrino charged-current interactions on hydrocarbon without pions in the final state using multiple detectors with correlated energy spectra at T2K

Authors: K. Abe, N. Akhlaq, R. Akutsu, H. Alarakia-Charles, A. Ali, Y. I. Alj Hakim, S. Alonso Monsalve, C. Alt, C. Andreopoulos, M. Antonova, S. Aoki, T. Arihara, Y. Asada, Y. Ashida, E. T. Atkin, M. Barbi, G. J. Barker, G. Barr, D. Barrow, M. Batkiewicz-Kwasniak, F. Bench, V. Berardi, L. Berns, S. Bhadra, A. Blanchet , et al. (380 additional authors not shown)

Abstract: 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 corre… ▽ More 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. △ Less

Submitted 18 October, 2023; v1 submitted 24 March, 2023; originally announced March 2023.

Comments: Updated discussion in Sec. V-A; Updated author list

Measurements of neutrino oscillation parameters from the T2K experiment using $3.6\times10^{21}$ protons on target

Authors: The T2K Collaboration, K. Abe, N. Akhlaq, R. Akutsu, A. Ali, S. Alonso Monsalve, C. Alt, C. Andreopoulos, M. Antonova, S. Aoki, T. Arihara, Y. Asada, Y. Ashida, E. T. Atkin, M. Barbi, G. J. Barker, G. Barr, D. Barrow, M. Batkiewicz-Kwasniak, F. Bench, V. Berardi, L. Berns, S. Bhadra, A. Blanchet, A. Blondel , et al. (376 additional authors not shown)

Abstract: 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 introdu… ▽ More 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θ_{13}$ and the impact of priors on the $δ_\mathrm{CP}$ measurement. Both analyses prefer the normal mass ordering and upper octant of $\sin^2θ_{23}$ with a nearly maximally CP-violating phase. Assuming the normal ordering and using the constraint on $\sin^2θ_{13}$ from reactors, $\sin^2θ_{23}=0.561^{+0.021}_{-0.032}$ using Feldman--Cousins corrected intervals, and $Δm^2_{32}=2.494_{-0.058}^{+0.041}\times10^{-3}~\mathrm{eV^2}$ using constant $Δχ^{2}$ intervals. The CP-violating phase is constrained to $δ_\mathrm{CP}=-1.97_{-0.70}^{+0.97}$ using Feldman--Cousins corrected intervals, and $δ_\mathrm{CP}=0,π$ is excluded at more than 90% confidence level. A Jarlskog invariant of zero is excluded at more than $2σ$ credible level using a flat prior in $δ_\mathrm{CP}$, and just below $2σ$ using a flat prior in $\sinδ_\mathrm{CP}$. When the external constraint on $\sin^2θ_{13}$ is removed, $\sin^2θ_{13}=28.0^{+2.8}_{-6.5}\times10^{-3}$, in agreement with measurements from reactor experiments. These results are consistent with previous T2K analyses. △ Less

Submitted 10 September, 2023; v1 submitted 6 March, 2023; originally announced March 2023.

Journal ref: Eur. Phys. J. C 83, 782 (2023)

A lensed radio jet at milli-arcsecond resolution II: Constraints on fuzzy dark matter from an extended gravitational arc

Authors: Devon M. Powell, Simona Vegetti, J. P. McKean, Simon D. M. White, Elisa G. M. Ferreira, Simon May, Cristiana Spingola

Abstract: 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_χ\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, su… ▽ More 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_χ\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_{χ,\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_χ$ 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. △ Less

Submitted 7 June, 2023; v1 submitted 21 February, 2023; originally announced February 2023.

Comments: 5 pages, 2 figures. Accepted in MNRAS Letters

Probing sub-galactic mass structure with the power spectrum of surface-brightness anomalies in high-resolution observations of galaxy-galaxy strong gravitational lenses. I. Power-spectrum measurement and feasibility study

Authors: D. Bayer, L. V. E. Koopmans, J. P. McKean, S. Vegetti, T. Treu, C. D. Fassnacht, K. Glazebrook

Abstract: 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 oppor… ▽ More 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. △ Less

Submitted 1 February, 2023; originally announced February 2023.

Comments: 20 pages, 19 figures, 2 tables, submitted to MNRAS

A Highly Magnified Gravitationally Lensed Red QSO at z = 2.5 with a Significant Flux Ratio Anomaly

Authors: Eilat Glikman, Cristian E. Rusu, Geoff C. -F. Chen, James Hung-Hsu Chan, Cristiana Spingola, Hannah Stacey, John McKean, Ciprian T. Berghea, S. G. Djorgovski, Matthew J. Graham, Daniel Stern, Tanya Urrutia, Mark Lacy, Nathan J. Secrest, John M. O'Meara

Abstract: We present the discovery of a gravitationally lensed dust-reddened QSO at z = 2.517, identified in a survey for QSOs by infrared selection. Hubble Space Telescope imaging reveals a quadruply lensed system in a cusp configuration, with a maximum image separation of ~1.8\arcsec. We find that compared to the central image of the cusp, the neighboring brightest image is anomalous by a factor of ~ 7 -… ▽ More We present the discovery of a gravitationally lensed dust-reddened QSO at z = 2.517, identified in a survey for QSOs by infrared selection. Hubble Space Telescope imaging reveals a quadruply lensed system in a cusp configuration, with a maximum image separation of ~1.8\arcsec. We find that compared to the central image of the cusp, the neighboring brightest image is anomalous by a factor of ~ 7 - 10, which is the largest flux anomaly measured to date in a lensed QSO. Incorporating high-resolution Jansky Very Large Array radio imaging and sub-mm imaging with the Atacama Large (sub-)Millimetre Array, we conclude that a low-mass perturber is the most likely explanation for the anomaly. The optical through near-infrared spectrum reveals that the QSO is moderately reddened with E(B - V) = 0.7 - 0.9. We see an upturn in the ultraviolet spectrum due to ~ 1% of the intrinsic emission being leaked back into the line of sight, which suggests that the reddening is intrinsic and not due to the lens. The QSO may have an Eddington ratio as high as L/L_Edd ~ 0.2. Consistent with previous red QSO samples, this source exhibits outflows in its spectrum as well as morphological properties suggestive of it being in a merger-driven transitional phase. We find a host-galaxy stellar mass of log M_*/M_Sun = 11.4, which is higher than the local M_BH vs. M_* relation, but consistent with other high redshift QSOs. When de-magnified, this QSO is at the knee of the luminosity function, allowing for the detailed study of a more typical moderate-luminosity infrared-selected QSO at high redshift. △ Less

Submitted 7 November, 2022; originally announced November 2022.

Comments: Accepted for publication in ApJ; 29 pages, 18 figures, 8 tables. arXiv admin note: text overlap with arXiv:1807.05434

A machine learning based approach to gravitational lens identification with the International LOFAR Telescope

Authors: S. Rezaei, J. P. McKean, M. Biehl, W. de Roo1, A. Lafontaine

Abstract: 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 Networ… ▽ More 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. △ Less

Submitted 21 July, 2022; originally announced July 2022.

Comments: Accepted to be published by MNRAS

Red quasars blow out molecular gas from galaxies during the peak of cosmic star formation

Authors: H. R. Stacey, T. Costa, J. P. McKean, C. E. Sharon, G. Calistro Rivera, E. Glikman, P. P. van der Werf

Abstract: 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 50… ▽ More 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. △ Less

Submitted 26 September, 2022; v1 submitted 19 July, 2022; originally announced July 2022.

Comments: Accepted by MNRAS. 18 figures and 3 tables

A lensed radio jet at milli-arcsecond resolution I: Bayesian comparison of parametric lens models

Authors: Devon M. Powell, Simona Vegetti, J. P. McKean, Cristiana Spingola, Hannah R. Stacey, Christopher D. Fassnacht

Abstract: 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… ▽ More 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. △ Less

Submitted 18 August, 2022; v1 submitted 7 July, 2022; originally announced July 2022.

Comments: 24 pages, 15 figures, 3 tables. Accepted in MNRAS

Using strong lensing to understand the microJy radio emission in two radio quiet quasars at redshift 1.7

Authors: P. Hartley, N. Jackson, S. Badole, J. P. McKean, D. Sluse, H. Vives-Arias

Abstract: 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… ▽ More 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. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Comments: Accepted by MNRAS; 15 pages, 9 figures

DECORAS: detection and characterization of radio-astronomical sources using deep learning

Authors: S. Rezaei, J. P. McKean, M. Biehl, A. Javadpour

Abstract: 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 posit… ▽ More 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$σ$ level, and has an almost factor of two improvement in reliability at 5.5$σ$. 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. △ Less

Submitted 21 September, 2021; v1 submitted 19 September, 2021; originally announced September 2021.

Comments: submitted to MNRAS

PS J1721+8842: A gravitationally lensed dual AGN system at redshift 2.37 with two radio components

Authors: C. S. Mangat, J. P. McKean, R. Brilenkov, P. Hartley, H. R. Stacey, S. Vegetti, D. Wen

Abstract: 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-AG… ▽ More 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. △ Less

Submitted 7 September, 2021; originally announced September 2021.

Comments: 5 pages, 3 figures, accepted for publication in MNRAS Letters

High-resolution imaging with the International LOFAR Telescope: Observations of the gravitational lenses MG 0751+2716 and CLASS B1600+434

Authors: Shruti Badole, Deepika Venkattu, Neal Jackson, Sarah Wallace, Jiten Dhandha, Philippa Hartley, Christopher Riddell-Rovira, Alice Townsend, Leah K. Morabito, J. P. McKean

Abstract: 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 stud… ▽ More 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. △ Less

Submitted 16 August, 2021; originally announced August 2021.

Comments: Accepted to a special issue of A&A on sub-arcsecond imaging with LOFAR

Sub-arcsecond imaging with the International LOFAR Telescope: II. Completion of the LOFAR Long-Baseline Calibrator Survey

Authors: Neal Jackson, Shruti Badole, John Morgan, Rajan Chhetri, Kaspars Prusis, Atvars Nikolajevs, Leah Morabito, Michiel Brentjens, Frits Sweijen, Marco Iacobelli, Emanuela Orrù, J. Sluman, R. Blaauw, H. Mulder, P. van Dijk, Sean Mooney, Adam Deller, Javier Moldon, J. R. Callingham, Jeremy Harwood, Martin Hardcastle, George Heald, Alexander Drabent, J. P. McKean, A. Asgekar , et al. (47 additional authors not shown)

Abstract: 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 t… ▽ More 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. △ Less

Submitted 16 August, 2021; originally announced August 2021.

Comments: Accepted to a special issue of A&A on sub-arcsecond imaging with LOFAR

Journal ref: A&A 658, A2 (2022)

Sub-arcsecond imaging with the International LOFAR Telescope I. Foundational calibration strategy and pipeline

Authors: L. K. Morabito, N. J. Jackson, S. Mooney, F. Sweijen, S. Badole, P. Kukreti, D. Venkattu, C. Groeneveld, A. Kappes, E. Bonnassieux, A. Drabent, M. Iacobelli, J. H. Croston, P. N. Best, M. Bondi, J. R. Callingham, J. E. Conway, A. T. Deller, M. J. Hardcastle, J. P. McKean, G. K. Miley, J. Moldon, H. J. A. Röttgering, C. Tasse, T. W. Shimwell , et al. (49 additional authors not shown)

Abstract: [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… ▽ More [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 $μ$Jy/bm image rms noise, which is dependent on the distance from the phase centre; typical values are ~90 $μ$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. △ Less

Submitted 16 August, 2021; originally announced August 2021.

Comments: Accepted to a special issue of A&A on sub-arcsecond imaging with LOFAR. 24 pages, 16 figures

Journal ref: A&A 658, A1 (2022)

SHARP VIII: J0924+0219 lens mass distribution and time-delay prediction through adaptive-optics imaging

Authors: Geoff C. -F. Chen, Christopher D. Fassnacht, Sherry H. Suyu, Léon V. E. Koopmans, David J. Lagattuta, John P. McKean, Matt W. Auger, Simona Vegetti, Tommaso Treu

Abstract: Strongly lensed quasars can provide measurements of the Hubble constant ($H_{0}$) independent of any other methods. One of the key ingredients is exquisite high-resolution imaging data, such as Hubble Space Telescope (HST) imaging and adaptive-optics (AO) imaging from ground-based telescopes, which provide strong constraints on the mass distribution of the lensing galaxy. In this work, we expand o… ▽ More Strongly lensed quasars can provide measurements of the Hubble constant ($H_{0}$) independent of any other methods. One of the key ingredients is exquisite high-resolution imaging data, such as Hubble Space Telescope (HST) imaging and adaptive-optics (AO) imaging from ground-based telescopes, which provide strong constraints on the mass distribution of the lensing galaxy. In this work, we expand on the previous analysis of three time-delay lenses with AO imaging (RXJ1131-1231, HE0435-1223, and PG1115+080), and perform a joint analysis of J0924+0219 by using AO imaging from the Keck Telescope, obtained as part of the SHARP (Strong lensing at High Angular Resolution Program) AO effort, with HST imaging to constrain the mass distribution of the lensing galaxy. Under the assumption of a flat $Λ$CDM model with fixed $Ω_{\rm m}=0.3$, we show that by marginalizing over two different kinds of mass models (power-law and composite models) and their transformed mass profiles via a mass-sheet transformation, we obtain $Δt_{\rm BA}h\hatσ_{v}^{-2}=6.89\substack{+0.8\\-0.7}$ days, $Δt_{\rm CA}h\hatσ_{v}^{-2}=10.7\substack{+1.6\\-1.2}$ days, and $Δt_{\rm DA}h\hatσ_{v}^{-2}=7.70\substack{+1.0\\-0.9}$ days, where $h=H_{0}/100~\rm km\,s^{-1}\,Mpc^{-1}$ is the dimensionless Hubble constant and $\hatσ_{v}=σ^{\rm ob}_{v}/(280~\rm km\,s^{-1})$ is the scaled dimensionless velocity dispersion. Future measurements of time delays with 10% uncertainty and velocity dispersion with 5% uncertainty would yield a $H_0$ constraint of $\sim15$% precision. △ Less

Submitted 21 July, 2021; originally announced July 2021.

Comments: 11 pages, 7 figures

LOFAR imaging of Cygnus A -- Direct detection of a turnover in the hotspot radio spectra

Authors: J. P. McKean, L. E. H. Godfrey, S. Vegetti, M. W. Wise, R. Morganti, M. J. Hardcastle, D. Rafferty, J. Anderson, I. M. Avruch, R. Beck, M. E. Bell, I. van Bemmel, M. J. Bentum, G. Bernardi, P. Best, R. Blaauw, A. Bonafede, F. Breitling, J. W. Broderick, M. Bruggen, L. Cerrigone, B. Ciardi, F. de Gasperin, A. Deller, S. Duscha , et al. (53 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 31 March, 2021; originally announced March 2021.

Comments: 9 pages, 3 figures, accepted for publication in MNRAS

Journal ref: MNRAS, 463, 3143 (2016)

Gravitational lensing in LoTSS DR2 -- Extremely faint 144-MHz radio emission from two highly magnified quasars

Authors: J. P. McKean, R. Luichies, A. Drabent, G. Gurkan, P. Hartley, A. Lafontaine, I. Prandoni, H. J. A. Rottgering, T. W. Shimwell, H. R. Stacey, C. Tasse

Abstract: 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 1… ▽ More 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. △ Less

Submitted 31 March, 2021; originally announced March 2021.

Comments: 5 pages, 1 figure, accepted for publication in MNRAS Letters

Implications of the Environments of Radio-detected AGN in a Complex Protostructure at z$\sim$3.3

Authors: Lu Shen, Brian C. Lemaux, Lori M. Lubin, Olga Cucciati, Olivier Le Fevre, Guilin Liu, Wenjuan Fang, Debora Pelliccia, Adam Tomczak, John McKean, Neal A. Miller, Christopher D. Fassnacht, Roy Gal, Denise Hung, Nimish Hathi, Sandro Bardelli, Daniela Vergani, Elena Zucca

Abstract: Radio Active Galactic Nuclei (RAGNs) are mainly found in dense structures (i.e., clusters/groups) at redshifts of z$<$2 and are commonly used to detect protoclusters at higher redshift. Here, we attempt to study the host and environmental properties of two relatively faint ($\mathrm L_\mathrm{1.4GHz} \sim10^{25}$ W Hz$^{-1}$) RAGNs in a known protocluster at z=3.3 in the PCl J0227-0421 field, dete… ▽ More Radio Active Galactic Nuclei (RAGNs) are mainly found in dense structures (i.e., clusters/groups) at redshifts of z$<$2 and are commonly used to detect protoclusters at higher redshift. Here, we attempt to study the host and environmental properties of two relatively faint ($\mathrm L_\mathrm{1.4GHz} \sim10^{25}$ W Hz$^{-1}$) RAGNs in a known protocluster at z=3.3 in the PCl J0227-0421 field, detected using the latest radio observation obtained as part of the Observations of Redshift Evolution in Large-Scale Environments (ORELSE) Survey. Using new spectroscopic observations obtained from Keck/MOSFIRE as part of the Charting Cluster Construction with the VIMOS Ultra-Deep Survey (VUDS) and ORELSE (C3VO) survey and previous spectroscopic data obtained as part of the VIMOS-VLT Deep Survey (VVDS) and VUDS, we revise the three-dimensional overdensity field around this protocluster. The protocluster is embedded in a large scale overdensity protostructure. This protostructure has an estimated total mass of $\sim$2.6$\times10^{15} M_\odot$ and contains several overdensity peaks. Both RAGNs are hosted by very bright and massive galaxies, while their hosts show extreme differences color, indicating that they have different ages and are in different evolutionary stages. Furthermore, we find that they are not in the most locally dense parts of the protostructure, but are fairly close to the centers of their parent overdensity peaks. We propose a scenario where merging might already have happened in both cases, which lowered the local density of their surrounding area and boosted their stellar mass. This work is the first time that two RAGNs at low luminosity have been found and studied within a high redshift protostructure. △ Less

Submitted 7 May, 2021; v1 submitted 4 March, 2021; originally announced March 2021.

Comments: 23 pages, 11 figures, 4 tables, accepted by ApJ

Joint constraints on thermal relic dark matter from strong gravitational lensing, the Lyman-$α$ forest, and Milky Way satellites

Authors: Wolfgang Enzi, Riccardo Murgia, Oliver Newton, Simona Vegetti, Carlos Frenk, Matteo Viel, Marius Cautun, Christopher D. Fassnacht, Matt Auger, Giulia Despali, John McKean, Léon V. E. Koopmans, Mark Lovell

Abstract: We derive joint constraints on the warm dark matter (WDM) half-mode scale by combining the analyses of a selection of astrophysical probes: strong gravitational lensing with extended sources, the Lyman-$α$ forest, and the number of luminous satellites in the Milky Way. We derive an upper limit of $λ_{\rm hm}=0.089{\rm~Mpc~h^{-1} }$ at the 95 per cent confidence level, which we show to be stable fo… ▽ More We derive joint constraints on the warm dark matter (WDM) half-mode scale by combining the analyses of a selection of astrophysical probes: strong gravitational lensing with extended sources, the Lyman-$α$ forest, and the number of luminous satellites in the Milky Way. We derive an upper limit of $λ_{\rm hm}=0.089{\rm~Mpc~h^{-1} }$ at the 95 per cent confidence level, which we show to be stable for a broad range of prior choices. Assuming a Planck cosmology and that WDM particles are thermal relics, this corresponds to an upper limit on the half-mode mass of $M_{\rm hm }< 3 \times 10^{7} {\rm~M_{\odot}~h^{-1}}$, and a lower limit on the particle mass of $m_{\rm th }> 6.048 {\rm~keV}$, both at the 95 per cent confidence level. We find that models with $λ_{\rm hm}> 0.223 {\rm~Mpc~h^{-1} }$ (corresponding to $m_{\rm th }> 2.552 {\rm~keV}$ and $M_{\rm hm }< 4.8 \times 10^{8} {\rm~M_{\odot}~h^{-1}}$) are ruled out with respect to the maximum likelihood model by a factor $\leq 1/20$. For lepton asymmetries $L_6>10$, we rule out the $7.1 {\rm~keV}$ sterile neutrino dark matter model, which presents a possible explanation to the unidentified $3.55 {\rm~keV}$ line in the Milky Way and clusters of galaxies. The inferred 95 percentiles suggest that we further rule out the ETHOS-4 model of self-interacting DM. Our results highlight the importance of extending the current constraints to lower half-mode scales. We address important sources of systematic errors and provide prospects for how the constraints of these probes can be improved upon in the future. △ Less

Submitted 22 June, 2021; v1 submitted 26 October, 2020; originally announced October 2020.

Comments: 16 pages, 3 Figures, minor updates of the results

Journal ref: Mon.Not.Roy.Astron.Soc. 506 (2021), 4, 5848-5862

Extended radio AGN at z ~ 1 in the ORELSE survey: The confining effect of dense environments

Authors: Lu Shen, Guilin Liu, Mengfei Zhang, Brian C. Lemaux, Lori M. Lubin, Debora Pelliccia, Emily Moravec, Emmet G. Golden-Marx, Hongyan Zhou, Wenjuan Fang, Adam Tomczak, John McKean, Neal A. Miller, Christopher D. Fassnacht, Po-Feng Wu, Dale D. Kocevski, Roy R. Gal, Denise Hung, Gordon K Squires

Abstract: Recent hydrodynamic simulations and observations of radio jets have shown that the surrounding environment has a large effect on their resulting morphology. To investigate this we use a sample of 50 Extended Radio Active Galactic Nuclei (ERAGN) detected in the Observations of Redshift Evolution in Large Scale Environments (ORELSE) survey. These sources are all successfully cross-identified to gala… ▽ More Recent hydrodynamic simulations and observations of radio jets have shown that the surrounding environment has a large effect on their resulting morphology. To investigate this we use a sample of 50 Extended Radio Active Galactic Nuclei (ERAGN) detected in the Observations of Redshift Evolution in Large Scale Environments (ORELSE) survey. These sources are all successfully cross-identified to galaxies within a redshift range of $0.55 \leq z \leq 1.35$, either through spectroscopic redshifts or accurate photometric redshifts. We find that ERAGN are more compact in high-density environments than those in low-density environments at a significance level of 4.5$σ$. Among a series of internal properties under our scrutiny, only the radio power demonstrates a positive correlation with their spatial extent. After removing the possible radio power effect, the difference of size in low- and high-density environments persists. In the global environment analyses, the majority (86\%) of high-density ERAGN reside in the cluster/group environment. In addition, ERAGN in the cluster/group central regions are preferentially compact with a small scatter in size, compared to those in the cluster/group intermediate regions and fields. In conclusion, our data appear to support the interpretation that the dense intracluster gas in the central regions of galaxy clusters plays a major role in confining the spatial extent of radio jets. △ Less

Submitted 7 September, 2020; originally announced September 2020.

The rocky road to quiescence: compaction and quenching of quasar host galaxies at z~2

Authors: H. R. Stacey, J. P. McKean, D. M. Powell, S. Vegetti, F. Rizzo, C. Spingola, M. W. Auger, R. J. Ivison, P. P. van der Werf

Abstract: 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) galaxie… ▽ More 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 ($Σ_{\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 ($Σ_{\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. △ Less

Submitted 3 November, 2020; v1 submitted 2 September, 2020; originally announced September 2020.

Comments: Accepted by MNRAS; 22 pages

A dynamically cold disk galaxy in the early Universe

Authors: F. Rizzo, S. Vegetti, D. Powell, F. Fraternali, J. P. McKean, H. R. Stacey, S. D. M. White

Abstract: 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 pr… ▽ More 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/σ= 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. △ Less

Submitted 2 September, 2020; originally announced September 2020.

Comments: Published in Nature on 12 August 2020. The published version is available at http://www.nature.com/articles/s41586-020-2572-6

Journal ref: Nature 584, 2020, 201-204

VLBI20-30: a scientific roadmap for the next decade -- The future of the European VLBI Network

Authors: Tiziana Venturi, Zsolt Paragi, Michael Lindqvist, Anna Bartkiewicz, Rob Beswick, Tamara Bogdanović, Walter Brisken, Patrick Charlot, Francisco Colomer, John Conway, Sándor Frey, José Carlos Guirado, Leonid Gurvits, Huib van Langevelde, Andrei Lobanov, John McKean, Raffaella Morganti, Tom Muxlow, Miguel Pérez-Torres, Kazi Rygl, Robert Schulz, Arpad Szomoru, Pablo de Vicente, Tao An, Guillem Anglada , et al. (55 additional authors not shown)

Abstract: This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, dem… ▽ More This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, demonstrate the EVN's potential to generate new and exciting results that will transform our view of the cosmos. Together with e-MERLIN, the EVN provides a range of baseline lengths that permit unique studies of faint radio sources to be made over a wide range of spatial scales. The science cases are reviewed in six chapters that cover the following broad areas: cosmology, galaxy formation and evolution, innermost regions of active galactic nuclei, explosive phenomena and transients, stars and stellar masers in the Milky Way, celestial reference frames and space applications. The document concludes with identifying the synergies with other radio, as well as multi-band/multi-messenger instruments, and provide the recommendations for future improvements. The appendices briefly describe other radio VLBI arrays, the technological framework for EVN developments, and a selection of spectral lines of astrophysical interest below 100 GHz. The document includes a glossary for non-specialists, and a list of acronyms at the end. △ Less

Submitted 5 July, 2020; originally announced July 2020.

Comments: Editors: Tiziana Venturi, Zsolt Paragi, Michael Lindqvist. "EVN Vision Document (2020)", 7 chapters, 3 appendices, 191 pages, 64 figures. Chapter coordinators and all the contributors are listed in the document

A novel approach to visibility-space modelling of interferometric gravitational lens observations at high angular resolution

Authors: Devon Powell, Simona Vegetti, John P. McKean, Cristiana Spingola, Francesca Rizzo, Hannah R. Stacey

Abstract: We present a new gravitational lens modelling technique designed to model high-resolution interferometric observations with large numbers of visibilities without the need to pre-average the data in time or frequency. We demonstrate the accuracy of the method using validation tests on mock observations. Using small data sets with $\sim 10^3$ visibilities, we first compare our approach with the more… ▽ More We present a new gravitational lens modelling technique designed to model high-resolution interferometric observations with large numbers of visibilities without the need to pre-average the data in time or frequency. We demonstrate the accuracy of the method using validation tests on mock observations. Using small data sets with $\sim 10^3$ visibilities, we first compare our approach with the more traditional direct Fourier transform (DFT) implementation and direct linear solver. Our tests indicate that our source inversion is indistinguishable from that of the DFT. Our method also infers lens parameters to within 1 to 2 per cent of both the ground truth and DFT, given sufficiently high signal-to-noise ratio (SNR). When the SNR is as low as 5, both approaches lead to errors of several tens of per cent in the lens parameters and a severely disrupted source structure, indicating that this is related to the SNR and choice of priors rather than the modelling technique itself. We then analyze a large data set with $\sim 10^8$ visibilities and a SNR matching real global Very Long Baseline Interferometry observations of the gravitational lens system MG J0751+2716. The size of the data is such that it cannot be modelled with traditional implementations. Using our novel technique, we find that we can infer the lens parameters and the source brightness distribution, respectively, with an RMS error of 0.25 and 0.97 per cent relative to the ground truth. △ Less

Submitted 21 October, 2020; v1 submitted 7 May, 2020; originally announced May 2020.

Comments: Accepted in MNRAS

Cassiopeia A, Cygnus A, Taurus A, and Virgo A at ultra-low radio frequencies

Authors: F. de Gasperin, J. Vink, J. P. McKean, A. Asgekar, M. J. Bentum, R. Blaauw, A. Bonafede, M. Bruggen, F. Breitling, W. N. Brouw, H. R. Butcher, B. Ciardi, V. Cuciti, M. de Vos, S. Duscha, J. Eisloffel, D. Engels, R. A. Fallows, T. M. O. Franzen, M. A. Garrett, A. W. Gunst, J. Horandel, G. Heald, L. V. E. Koopmans, A. Krankowski , et al. (27 additional authors not shown)

Abstract: 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. Fur… ▽ More 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. △ Less

Submitted 24 February, 2020; originally announced February 2020.

Comments: 7 pages, 2 figures, accepted A&A, online data on A&A website

Smoke on the water: CO and H2O in a circumnuclear disc around a quasar at redshift 2.64

Authors: H. R. Stacey, A. Lafontaine, J. P. McKean

Abstract: 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 repor… ▽ More 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. △ Less

Submitted 17 February, 2020; originally announced February 2020.

Comments: Accepted by MNRAS. 11 pages, 9 figures

Full of Orions: a 200-pc mapping of the interstellar medium in the redshift-3 lensed dusty star-forming galaxy SDP.81

Authors: M. Rybak, J. A. Hodge, S. Vegetti, P. van der Werf, P. Andreani, L. Graziani, J. P. McKean

Abstract: 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)… ▽ More 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. △ Less

Submitted 6 April, 2020; v1 submitted 28 December, 2019; originally announced December 2019.

Comments: 27 pages, 22 figures. 24 MB in total. MNRAS accepted 2020 March 17. Data available on request

The Properties of Radio and Mid-infrared Detected Galaxies and the Effect of Environment on the Co-evolution of AGN and Star Formation at $z \sim 1$

Authors: Lu Shen, Brian C. Lemaux, Lori M. Lubin, John McKean, Neal A. Miller, Debora Pelliccia, Christopher D. Fassnacht, Adam Tomczak, Po-Feng Wu, Dale Kocevski, Roy Gal, Denise Hung, Gordon Squires

Abstract: In this study we investigate 179 radio-IR galaxies drawn from a sample of spectroscopically-confirmed galaxies that are detected in radio and mid-infrared (MIR) in the redshift range of $0.55 \leq z \leq 1.30$ in the Observations of Redshift Evolution in Large Scale Environments (ORELSE) survey. We constrain the Active Galactic Nuclei (AGN) contribution in the total IR luminosity (f… ▽ More In this study we investigate 179 radio-IR galaxies drawn from a sample of spectroscopically-confirmed galaxies that are detected in radio and mid-infrared (MIR) in the redshift range of $0.55 \leq z \leq 1.30$ in the Observations of Redshift Evolution in Large Scale Environments (ORELSE) survey. We constrain the Active Galactic Nuclei (AGN) contribution in the total IR luminosity (f$_{\text{AGN}}$), and estimate the AGN luminosity (L$_{\text{AGN}}$) and the star formation rate (SFR) using the CIGALE Spectral Energy Distribution (SED) fitting routine. Based on the f$_{\text{AGN}}$ and radio luminosity, radio-IR galaxies are split into: galaxies that host either high or low f$_{\text{AGN}}$ AGN (high-/low-f$_{\text{AGN}}$), and star forming galaxies with little to no AGN activity (SFGs). We study the colour, stellar mass, radio luminosity, L$_{\text{AGN}}$ and SFR properties of the three radio-IR sub-samples, comparing to a spec-IR sample drawn from spectroscopically-confirmed galaxies that are also detected in MIR. No significant difference between radio luminosity of these sub-samples was found, which could be due to the combined contribution of radio emission from AGN and star formation. We find a positive relationship between L$_{\text{AGN}}$and specific SFR (sSFR) for both AGN sub-samples, strongly suggesting a co-evolution scenario of AGN and SF in these galaxies. A toy model is designed to demonstrate this co-evolution scenario, where we find that, in almost all cases, a rapid quenching timescale is required, which we argue is a signature of AGN quenching. The environmental preference for intermediate/infall regions of clusters/groups remains across the co-evolution scenario, which suggests that galaxies might be in an orbital motion around the cluster/group during the scenario. △ Less

Submitted 31 October, 2019; originally announced November 2019.

Proper motion in lensed radio jets at redshift 3: a possible dual super-massive black hole system in the early Universe

Authors: C. Spingola, J. P. McKean, D. Massari, L. V. E. Koopmans

Abstract: 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 $μ$as… ▽ More 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 $μ$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 $μ$as yr$^{-1}$ for the two components at highest magnification ($μ\sim350$) and of the order of a few mas yr$^{-1}$ for the two components at lower magnification ($μ\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. △ Less

Submitted 30 August, 2019; originally announced August 2019.

Comments: 11 pages, 4 figures; accepted for publication by Astronomy and Astrophysics

Journal ref: A&A 630, A108 (2019)

A SHARP view of H0LiCOW: $H_{0}$ from three time-delay gravitational lens systems with adaptive optics imaging

Authors: Geoff C. -F. Chen, Christopher D. Fassnacht, Sherry. H. Suyu, Cristian E. Rusu, James H. H. Chan, Kenneth C. Wong, Matthew W. Auger, Stefan Hilbert, Vivien Bonvin, Simon Birrer, Martin Millon, Leon V. E. Koopmans, David J. Lagattuta, John P. McKean, Simona Vegetti, Frederic Courbin, Xuheng Ding, Aleksi Halkola, Inh Jee, Anowar J. Shajib, Dominique Sluse, Alessandro Sonnenfeld, Tommaso Treu

Abstract: We present the measurement of the Hubble Constant, $H_0$, with three strong gravitational lens systems. We describe a blind analysis of both PG1115+080 and HE0435-1223 as well as an extension of our previous analysis of RXJ1131-1231. For each lens, we combine new adaptive optics (AO) imaging from the Keck Telescope, obtained as part of the SHARP AO effort, with Hubble Space Telescope (HST) imaging… ▽ More We present the measurement of the Hubble Constant, $H_0$, with three strong gravitational lens systems. We describe a blind analysis of both PG1115+080 and HE0435-1223 as well as an extension of our previous analysis of RXJ1131-1231. For each lens, we combine new adaptive optics (AO) imaging from the Keck Telescope, obtained as part of the SHARP AO effort, with Hubble Space Telescope (HST) imaging, velocity dispersion measurements, and a description of the line-of-sight mass distribution to build an accurate and precise lens mass model. This mass model is then combined with the COSMOGRAIL measured time delays in these systems to determine $H_{0}$. We do both an AO-only and an AO+HST analysis of the systems and find that AO and HST results are consistent. After unblinding, the AO-only analysis gives $H_{0}=82.8^{+9.4}_{-8.3}~\rm km\,s^{-1}\,Mpc^{-1}$ for PG1115+080, $H_{0}=70.1^{+5.3}_{-4.5}~\rm km\,s^{-1}\,Mpc^{-1}$ for HE0435-1223, and $H_{0}=77.0^{+4.0}_{-4.6}~\rm km\,s^{-1}\,Mpc^{-1}$ for RXJ1131-1231. The joint AO-only result for the three lenses is $H_{0}=75.6^{+3.2}_{-3.3}~\rm km\,s^{-1}\,Mpc^{-1}$. The joint result of the AO+HST analysis for the three lenses is $H_{0}=76.8^{+2.6}_{-2.6}~\rm km\,s^{-1}\,Mpc^{-1}$. All of the above results assume a flat $Λ$ cold dark matter cosmology with a uniform prior on $Ω_{\textrm{m}}$ in [0.05, 0.5] and $H_{0}$ in [0, 150] $\rm km\,s^{-1}\,Mpc^{-1}$. This work is a collaboration of the SHARP and H0LiCOW teams, and shows that AO data can be used as the high-resolution imaging component in lens-based measurements of $H_0$. The full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper. △ Less

Submitted 4 July, 2019; originally announced July 2019.

Comments: 33 pages, 23 figures

SHARP -- VI. Evidence for CO (1-0) molecular gas extended on kpc-scales in AGN star forming galaxies at high redshift

Authors: C. Spingola, J. P. McKean, S. Vegetti, D. Powell, M. W. Auger, L. V. E. Koopmans, C. D. Fassnacht, D. J. Lagattuta, F. Rizzo, H. R. Stacey, F. Sweijen

Abstract: 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… ▽ More 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 $α_{\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. △ Less

Submitted 11 May, 2020; v1 submitted 15 May, 2019; originally announced May 2019.

Comments: 22 pages, 15 figures with additional 1 figure in the Appendix. Accepted for publication by MNRAS

SHARP -- VII. New constraints on the dark matter free-streaming properties and substructure abundance from gravitationally lensed quasars

Authors: Jen-Wei Hsueh, Wolfgang Enzi, Simona Vegetti, Matt Auger, Christopher D. Fassnacht, Giulia Despali, Leon V. E. Koopmans, John P. McKean

Abstract: We present an analysis of seven strongly gravitationally lensed quasars and the corresponding constraints on the properties of dark matter. Our results are derived by modelling the lensed image positions and flux-ratios using a combination of smooth macro models and a population of low-mass haloes within the mass range 10^6 to 10^9 Msun. Our lens models explicitly include higher-order complexity i… ▽ More We present an analysis of seven strongly gravitationally lensed quasars and the corresponding constraints on the properties of dark matter. Our results are derived by modelling the lensed image positions and flux-ratios using a combination of smooth macro models and a population of low-mass haloes within the mass range 10^6 to 10^9 Msun. Our lens models explicitly include higher-order complexity in the form of stellar discs and luminous satellites, as well as low-mass haloes located along the observed lines of sight for the first time. Assuming a Cold Dark Matter (CDM) cosmology, we infer an average total mass fraction in substructure of f_sub = 0.012^{+0.007}_{-0.004} (68 per cent confidence limits), which is in agreement with the predictions from CDM hydrodynamical simulations to within 1 sigma. This result is closer to the predictions than those from previous studies that did not include line-of-sight haloes. Under the assumption of a thermal relic dark matter model, we derive a lower limit on the particle relic mass of m th > 5.58 keV (95 per cent confidence limits), which is consistent with a value of m_th > 5.3 keV from the recent analysis of the Ly-alpha forest. We also identify two main sources of possible systematic errors and conclude that deeper investigations in the complex structure of lens galaxies as well as the size of the background sources should be a priority for this field. △ Less

Submitted 19 December, 2019; v1 submitted 10 May, 2019; originally announced May 2019.

Comments: 14 pages, 7 figures, accepted for publication in MNRAS

Gravitational lensing at milliarcsecond angular resolution with VLBI observations

Authors: C. Spingola, J. P. McKean, A. Deller, J. Moldon

Abstract: 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-reso… ▽ More 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. △ Less

Submitted 19 February, 2019; originally announced February 2019.

Comments: 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

A novel search for gravitationally lensed radio sources in wide-field VLBI imaging from the mJIVE-20 survey

Authors: C. Spingola, J. P. McKean, M. Lee, A. Deller, J. Moldon

Abstract: 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 compon… ▽ More 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. △ Less

Submitted 22 November, 2018; originally announced November 2018.

Comments: 31 pages, 22 figures; accepted for publication in MNRAS

LoTSS/HETDEX: Disentangling star formation and AGN activity in gravitationally-lensed radio-quiet quasars

Authors: H. R. Stacey, J. P. McKean, N. J. Jackson, P. N. Best, G. Calistro Rivera, J. R. Callingham, K. J. Duncan, G. Gürkan, M. J. Hardcastle, M. Iacobelli, A. P. Mechev, L. K. Morabito, I. Prandoni, H. J. A. Röttgering, J. Sabater, T. W. Shimwell, C. Tasse, W. L. Williams

Abstract: 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… ▽ More 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. △ Less

Submitted 19 November, 2018; originally announced November 2018.

Comments: 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) III. First Data Release: optical/IR identifications and value-added catalogue

Authors: W. L. Williams, M. J. Hardcastle, P. N. Best, J. Sabater, J. H. Croston, K. J. Duncan, T. W. Shimwell, H. J. A. Röttgering, D. Nisbet, G. Gürkan, L. Alegre, R. K. Cochrane, A. Goyal, C. L. Hale, N. Jackson, M. Jamrozy, R. Kondapally, M. Kunert-Bajraszewska, V. H. Mahatma, B. Mingo, L. K. Morabito, I. Prandoni, C. Roskowinski, A. Shulevski, D. J. B. Smith , et al. (16 additional authors not shown)

Abstract: 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 th… ▽ More 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. △ Less

Submitted 19 November, 2018; originally announced November 2018.

Comments: 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

Journal ref: A&A 622, A2 (2019)

The LOFAR Two-metre Sky Survey - II. First data release

Authors: T. W. Shimwell, C. Tasse, M. J. Hardcastle, A. P. Mechev, W. L. Williams, P. N. Best, H. J. A. Röttgering, J. R. Callingham, T. J. Dijkema, F. de Gasperin, D. N. Hoang, B. Hugo, M. Mirmont, J. B. R. Oonk, I. Prandoni, D. Rafferty, J. Sabater, O. Smirnov, R. J. van Weeren, G. J. White, M. Atemkeng, L. Bester, E. Bonnassieux, M. Brüggen, G. Brunetti , et al. (82 additional authors not shown)

Abstract: 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… ▽ More 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\,μ$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. △ Less

Submitted 19 November, 2018; originally announced November 2018.

Comments: 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

Resolving on 100 pc scales the UV-continuum in Lyman-$α$ emitters between redshift 2 to 3 with gravitational lensing

Authors: E. Ritondale, M. W. Auger, S. Vegetti, J. P. McKean

Abstract: We present a study of seventeen LAEs at redshift 2$<z<$3 gravitationally lensed by massive early-type galaxies (ETGs) at a mean redshift of approximately 0.5. Using a fully Bayesian grid-based technique, we model the gravitational lens mass distributions with elliptical power-law profiles and reconstruct the UV-continuum surface brightness distributions of the background sources using pixellated s… ▽ More We present a study of seventeen LAEs at redshift 2$<z<$3 gravitationally lensed by massive early-type galaxies (ETGs) at a mean redshift of approximately 0.5. Using a fully Bayesian grid-based technique, we model the gravitational lens mass distributions with elliptical power-law profiles and reconstruct the UV-continuum surface brightness distributions of the background sources using pixellated source models. We find that the deflectors are close to, but not consistent with isothermal models in almost all cases, at the $2σ$-level. We take advantage of the lensing magnification (typically $μ\simeq$ 20) to characterise the physical and morphological properties of these LAE galaxies. From reconstructing the ultra-violet continuum emission, we find that the star-formation rates range from 0.3 to 8.5 M$_{\odot}$ yr$^{-1}$ and that the galaxies are typically composed of several compact and diffuse components, separated by 0.4 to 4 kpc. Moreover, they have peak star-formation rate intensities that range from 2.1 to 54.1 M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$. These galaxies tend to be extended with major axis ranging from 0.2 to 1.8 kpc (median 561 pc), and with a median ellipticity of 0.49. This morphology is consistent with disk-like structures of star-formation for more than half of the sample. However, for at least two sources, we also find off-axis components that may be associated with mergers. Resolved kinematical information will be needed to confirm the disk-like nature and possible merger scenario for the LAEs in the sample. △ Less

Submitted 8 November, 2018; originally announced November 2018.

Comments: 19 pages, 7 figures, accepted for publication on MNRAS

Low-mass halo perturbations in strong gravitational lenses at redshift z$\sim$0.5 are consistent with CDM

Authors: E. Ritondale, S. Vegetti, G. Despali, M. W. Auger, L. V. E Koopmans, J. P. McKean

Abstract: 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$σ$, we report no significant detection in any of the lenses. Using the sensitivity function at the 10-$σ$ level, we… ▽ More 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$σ$, we report no significant detection in any of the lenses. Using the sensitivity function at the 10-$σ$ level, we have calculated the predicted number of detectable cold dark matter (CDM) line-of-sight haloes to be $μ_{l} = 1.17\pm1.08$, in agreement with our null detection. Assuming a detection sensitivity that improved to the level implied by a 5-$σ$ threshold, the expected number of detectable line-of-sight haloes rises to $μ_l = 9.0\pm3.0$. Whilst the current data find zero detections at this sensitivity level (which has a probability of P$^{{\rm5}σ}_{\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-$σ$ 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-$σ$ level. △ Less

Submitted 8 November, 2018; originally announced November 2018.

Comments: 14 pages, 7 figures, submitted for publication on MNRAS

Constraining the mass density of free-floating black holes using razor-thin lensing arcs

Authors: Uddipan Banik, Frank C. van den Bosch, Michael Tremmel, Anupreeta More, Giulia Despali, Surhud More, Simona Vegetti, John P. McKean

Abstract: Strong lensing of active galactic nuclei in the radio can result in razor-thin arcs, with a thickness of less than a milli-arcsecond, if observed at the resolution achievable with very long baseline interferometry (VLBI). Such razor-thin arcs provide a unique window on the coarseness of the matter distribution between source and observer. In this paper, we investigate to what extent such razor-thi… ▽ More Strong lensing of active galactic nuclei in the radio can result in razor-thin arcs, with a thickness of less than a milli-arcsecond, if observed at the resolution achievable with very long baseline interferometry (VLBI). Such razor-thin arcs provide a unique window on the coarseness of the matter distribution between source and observer. In this paper, we investigate to what extent such razor-thin arcs can constrain the number density and mass function of `free-floating' black holes, defined as black holes that do not, or no longer, reside at the centre of a galaxy. These can be either primordial in origin or arise as by-products of the evolution of super-massive black holes in galactic nuclei. When sufficiently close to the line of sight, free-floating black holes cause kink-like distortions in the arcs, which are detectable by eye in the VLBI images as long as the black hole mass exceeds $\sim 1000$ Solar masses. Using a crude estimate for the detectability of such distortions, we analytically compute constraints on the matter density of free-floating black holes resulting from null-detections of distortions along a realistic, fiducial arc, and find them to be comparable to those from quasar milli-lensing. We also use predictions from a large hydrodynamical simulation for the demographics of free-floating black holes that are not primordial in origin, and show that their predicted mass density is roughly four orders of magnitude below the constraints achievable with a single razor-thin arc. △ Less

Submitted 19 December, 2018; v1 submitted 1 November, 2018; originally announced November 2018.

Comments: 17 pages, 13 figures, 1 table, comments welcome

Journal ref: Mon. Not. R. Astron. Soc. 483, 2, 1558-1573 (2019)

A flux-ratio anomaly in the CO spectral line emission from gravitationally-lensed quasar MG J0414+0534

Authors: H. R. Stacey, J. P. McKean

Abstract: 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 imag… ▽ More 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. △ Less

Submitted 16 August, 2018; originally announced August 2018.

Comments: Accepted as MNRAS Letter