Showing 1–50 of 186 results for author: Fassnacht, C D

TDCOSMO. XVII. Measurement of the Hubble Constant from the Lensed Quasar WGD$\,$2038$-$4008

Authors: Kenneth C. Wong, Frédéric Dux, Anowar J. Shajib, Sherry H. Suyu, Martin Millon, Pritom Mozumdar, Patrick R. Wells, Adriano Agnello, Simon Birrer, Elizabeth J. Buckley-Geer, Frédéric Courbin, Christopher D. Fassnacht, Joshua Frieman, Aymeric Galan, Huan Lin, Philip J. Marshall, Jason Poh, Stefan Schuldt, Dominique Sluse, Tommaso Treu

Abstract: Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant ($H_{0}$). The TDCOSMO collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar~WGD$\,$2038$-$4008~using new time-delay measurements and previous mass models by TDCOSMO. Th… ▽ More Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant ($H_{0}$). The TDCOSMO collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar~WGD$\,$2038$-$4008~using new time-delay measurements and previous mass models by TDCOSMO. This is the first TDCOSMO lens to incorporate multiple lens modeling codes and the full time-delay covariance matrix into the cosmological inference. The models are fixed before the time delay is measured, and the analysis is performed blinded with respect to the cosmological parameters to prevent unconscious experimenter bias. We obtain $D_{Δt} = 1.68^{+0.40}_{-0.38}$ Gpc using two families of mass models, a power-law describing the total mass distribution, and a composite model of baryons and dark matter, although the composite model is disfavored due to kinematics constraints. In a flat $Λ$CDM cosmology, we constrain the Hubble constant to be $H_{0} = 65^{+23}_{-14}\, \rm km\ s^{-1}\,Mpc^{-1}$. The dominant source of uncertainty comes from the time delays, due to the low variability of the quasar. Future long-term monitoring, especially in the era of the Vera C. Rubin Observatory's Legacy Survey of Space and Time, could catch stronger quasar variability and further reduce the uncertainties. This system will be incorporated into an upcoming hierarchical analysis of the entire TDCOSMO sample, and improved time delays and spatially-resolved stellar kinematics could strengthen the constraints from this system in the future. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: 9 pages, 5 figures, 3 tables; submitted for publication in Astronomy & Astrophysics

TDCOSMO XV: Population Analysis of Lines of Sight of 25 Strong Galaxy-Galaxy Lenses with Extreme Value Statistics

Authors: Patrick R. Wells, Christopher D. Fassnacht, Simon Birrer, Devon Williams

Abstract: Time-Delay Cosmography is a technique for measuring $H_0$ with strong gravitational lensing. It requires a correction for line of sight perturbations, and it is necessary to build tools to assess populations of these lines of sight efficiently. We aim demonstrate the techniques necessary to analyze line of sight effects at a population level, and investigate whether strong lenses fall in preferabl… ▽ More Time-Delay Cosmography is a technique for measuring $H_0$ with strong gravitational lensing. It requires a correction for line of sight perturbations, and it is necessary to build tools to assess populations of these lines of sight efficiently. We aim demonstrate the techniques necessary to analyze line of sight effects at a population level, and investigate whether strong lenses fall in preferably overdense environments. We analyze a set of 25 galaxy-galaxy lens lines of sight in the Strong Lensing Legacy Survey sample using standard techniques, then perform a hierarchical analysis to constrain the population-level parameters. We introduce a new statistical model for these posteriors that may provide insight into the underlying physics of the system. We find the median value of $κ_{\rm{ext}}$ in the population model to be $0.033 \pm 0.010$. The median value of $κ_{\rm{ext}}$ for the individual lens posteriors is $0.008 \pm 0.015$. Both approaches demostrate that our systems are drawn from an overdense sample. The different results from these two approaches show the importance of population models that do not multiply the effect of our priors. △ Less

Submitted 17 June, 2024; v1 submitted 15 March, 2024; originally announced March 2024.

Comments: 12 pages, 4 figures, 1 table. Accepted for Publication in A&A

General Multipoles and Their Implications for Dark Matter Inference

Authors: Jacob S. Cohen, Christopher D. Fassnacht, Conor M. O'Riordan, Simona Vegetti

Abstract: The flux ratios of strongly lensed quasars have previously been used to infer the properties of dark matter. In these analyses it is crucial to separate the effect of the main lensing galaxy and the low-mass dark matter halo population. In this work, we investigate flux-ratio perturbations resulting from general third- and fourth-order multipole perturbations to the main lensing galaxy's mass prof… ▽ More The flux ratios of strongly lensed quasars have previously been used to infer the properties of dark matter. In these analyses it is crucial to separate the effect of the main lensing galaxy and the low-mass dark matter halo population. In this work, we investigate flux-ratio perturbations resulting from general third- and fourth-order multipole perturbations to the main lensing galaxy's mass profile. We simulate four lens systems, each with a different lensing configuration, without multipoles. The simulated flux ratios are perturbed by 10-40 per cent by a population of low-mass haloes consistent with CDM and, in one case, also a satellite galaxy. This level of perturbation is comparable to the magnitude of flux-ratio anomalies in real data that has been previously analyzed. We then attempt to fit the simulated systems using multipoles instead of low-mass haloes. We find that multipoles with amplitudes of 0.01 or less can produce flux-ratio perturbations in excess of 40 per cent. In all cases, third- or fourth-order multipoles can individually reduce the magnitude of, if not eliminate, flux-ratio anomalies. When both multipole orders are jointly included, all simulated flux ratios can be fit to within the observational uncertainty. Our results indicate that low-mass haloes and multipoles are highly degenerate when modelling quadruply-imaged quasars based just on image positions and flux ratios. In the presence of this degeneracy, flux-ratio anomalies in lensed quasars alone cannot be used to place strong constraints on the properties of dark matter without additional information that can inform our priors. △ Less

Submitted 15 March, 2024; v1 submitted 13 March, 2024; originally announced March 2024.

Comments: 13 pages, 9 figures, to be submitted to MNRAS; figures 2-7 updated for clarity

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

JWST lensed quasar dark matter survey I: Description and First Results

Authors: A. M. Nierenberg, R. E. Keeley, D. Sluse, D. Gilman, S. Birrer, T. Treu, K. N. Abazajian, T. Anguita, A. J. Benson, V. N. Bennert, S. G. Djorgovski, X. Du, C. D. Fassnacht, S. F. Hoenig, A. Kusenko, C. Lemon, M. Malkan, V. Motta, L. A. Moustakas, D. Stern, R. H. Wechsler

Abstract: The flux ratios of gravitationally lensed quasars provide a powerful probe of the nature of dark matter. Importantly, these ratios are sensitive to small-scale structure, irrespective of the presence of baryons. This sensitivity may allow us to study the halo mass function even below the scales where galaxies form observable stars. For accurate measurements, it is essential that the quasar's light… ▽ More The flux ratios of gravitationally lensed quasars provide a powerful probe of the nature of dark matter. Importantly, these ratios are sensitive to small-scale structure, irrespective of the presence of baryons. This sensitivity may allow us to study the halo mass function even below the scales where galaxies form observable stars. For accurate measurements, it is essential that the quasar's light is emitted from a physical region of the quasar with an angular scale of milli-arcseconds or larger; this minimizes microlensing effects by stars within the deflector. The warm dust region of quasars fits this criterion, as it has parsec-size physical scales and dominates the spectral energy distribution of quasars at wavelengths greater than 10$μ$m. The JWST Mid-Infrared Instrument (MIRI) is adept at detecting redshifted light in this wavelength range, offering both the spatial resolution and sensitivity required for accurate gravitational lensing flux ratio measurements. Here, we introduce our survey designed to measure the warm dust flux ratios of 31 lensed quasars. We discuss the flux-ratio measurement technique and present results for the first target, DES J0405-3308. We find that we can measure the quasar warm dust flux ratios with 3% precision. Our simulations suggest that this precision makes it feasible to detect the presence of 10$^7$ M$_\odot$ dark matter halos at cosmological distances. Such halos are expected to be completely dark in Cold Dark Matter models. △ Less

Submitted 18 September, 2023; originally announced September 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.

TDCOSMO XIV: Practical Techniques for Estimating External Convergence of Strong Gravitational Lens Systems and Applications to the SDSS J0924+0219 System

Authors: Patrick Wells, Christopher D. Fassnacht, C. E. Rusu

Abstract: Time-delay cosmography uses strong gravitational lensing of a time-variable source to infer the Hubble Constant. The measurement is independent from both traditional distance ladder and CMB measurements. An accurate measurement with this technique requires considering the effects of objects along the line of sight outside the primary lens, which is quantified by the external convergence (… ▽ More Time-delay cosmography uses strong gravitational lensing of a time-variable source to infer the Hubble Constant. The measurement is independent from both traditional distance ladder and CMB measurements. An accurate measurement with this technique requires considering the effects of objects along the line of sight outside the primary lens, which is quantified by the external convergence ($κ_{\rm{ext}}$). In absence of such corrections, $H_0$ will be biased towards higher values in overdense fields and lower values in underdense fields. We discuss the current state of the methods used to account for environment effects. We present a new software package built for this kind of analysis and others that can leverage large astronomical survey datasets. We apply these techniques to the SDSS J0924+0219 strong lens field. We infer the relative density of the SDSS J0924+0219 field by computing weighted number counts for all galaxies in the field, and comparing to weighted number counts computed for a large number of fields in a reference survey. We then compute weighted number counts in the Millennium Simulation and compare these results to infer the external convergence of the lens field.Results. Our results show the SDSS J0924+0219 field is a fairly typical line of sight, with median $κ_{\rm{ext}} = -0.012$ and standard deviation $σ_κ = 0.028$. △ Less

Submitted 31 May, 2023; v1 submitted 6 February, 2023; originally announced February 2023.

Comments: Submitted to A&A. 10 pages, 5 figures

Journal ref: A&A 676, A95 (2023)

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

TDCOSMO. XII. Improved Hubble constant measurement from lensing time delays using spatially resolved stellar kinematics of the lens galaxy

Authors: Anowar J. Shajib, Pritom Mozumdar, Geoff C. -F. Chen, Tommaso Treu, Michele Cappellari, Shawn Knabel, Sherry H. Suyu, Vardha N. Bennert, Joshua A. Frieman, Dominique Sluse, Simon Birrer, Frederic Courbin, Christopher D. Fassnacht, Lizvette Villafaña, Peter R. Williams

Abstract: Strong-lensing time delays enable measurement of the Hubble constant ($H_{0}$) independently of other traditional methods. The main limitation to the precision of time-delay cosmography is mass-sheet degeneracy (MSD). Some of the previous TDCOSMO analyses broke the MSD by making standard assumptions about the mass density profile of the lens galaxy, reaching 2% precision from seven lenses. However… ▽ More Strong-lensing time delays enable measurement of the Hubble constant ($H_{0}$) independently of other traditional methods. The main limitation to the precision of time-delay cosmography is mass-sheet degeneracy (MSD). Some of the previous TDCOSMO analyses broke the MSD by making standard assumptions about the mass density profile of the lens galaxy, reaching 2% precision from seven lenses. However, this approach could potentially bias the $H_0$ measurement or underestimate the errors. For this work, we broke the MSD for the first time using spatially resolved kinematics of the lens galaxy in RXJ1131$-$1231 obtained from the Keck Cosmic Web Imager spectroscopy, in combination with previously published time delay and lens models derived from Hubble Space Telescope imaging. This approach allowed us to robustly estimate $H_0$, effectively implementing a maximally flexible mass model. Following a blind analysis, we estimated the angular diameter distance to the lens galaxy $D_{\rm d} = 865_{-81}^{+85}$ Mpc and the time-delay distance $D_{Δt} = 2180_{-271}^{+472}$ Mpc, giving $H_0 = 77.1_{-7.1}^{+7.3}$ km s$^{-1}$ Mpc$^{-1}$ - for a flat $Λ$ cold dark matter cosmology. The error budget accounts for all uncertainties, including the MSD inherent to the lens mass profile and the line-of-sight effects, and those related to the mass-anisotropy degeneracy and projection effects. Our new measurement is in excellent agreement with those obtained in the past using standard simply parametrized mass profiles for this single system ($H_0 = 78.3^{+3.4}_{-3.3}$ km s$^{-1}$ Mpc$^{-1}$) and for seven lenses ($H_0 = 74.2_{-1.6}^{+1.6}$ km s$^{-1}$ Mpc$^{-1}$), or for seven lenses using single-aperture kinematics and the same maximally flexible models used by us ($H_0 = 73.3^{+5.8}_{-5.8}$ km s$^{-1}$ Mpc$^{-1}$). This agreement corroborates the methodology of time-delay cosmography. △ Less

Submitted 19 November, 2023; v1 submitted 6 January, 2023; originally announced January 2023.

Comments: 21 pages, 22 figures, 1 table. Accepted by A&A (this version: language improvements only)

Journal ref: A&A 673, A9 (2023)

TDCOSMO. XI. New lensing galaxy redshift and velocity dispersion measurements from Keck spectroscopy of eight lensed quasar systems

Authors: P. Mozumdar, C. D. Fassnacht, T. Treu, C. Spiniello, A. J. Shajib

Abstract: We have measured the redshifts and single-aperture velocity dispersions of eight lens galaxies using the data collected by the Echellette Spectrograph and Imager (ESI) and Low Resolution Imaging Spectrometer (LRIS) at W.M. Keck observatory on different observing nights spread over three years (2018-2020). These results, combined with other ancillary data, such as high-resolution images of the lens… ▽ More We have measured the redshifts and single-aperture velocity dispersions of eight lens galaxies using the data collected by the Echellette Spectrograph and Imager (ESI) and Low Resolution Imaging Spectrometer (LRIS) at W.M. Keck observatory on different observing nights spread over three years (2018-2020). These results, combined with other ancillary data, such as high-resolution images of the lens systems, and time delays, are necessary to increase the sample size of the quasar-galaxy lens systems for which the Hubble constant can be measured, using the time-delay strong lensing method, hence increasing the precision of its inference. Typically, the 2D spectra of the quasar-galaxy lens systems get spatially blended due to seeing by ground-based observations. As a result, the extracted lensing galaxy (deflector) spectra become significantly contaminated by quasar light, which affects the ability to extract meaningful information about the deflector. To account for spatial blending and extract less contaminated and higher signal-to-noise ratio (S/N) 1D spectra of the deflectors, a forward modeling method has been implemented. From the extracted spectra, we have measured redshifts using prominent absorption lines and single aperture velocity dispersions using the penalized pixel fitting code pPXF. In this paper, we report the redshifts and single aperture velocity dispersions of eight lens galaxies - J0147+4630, B0445+123, B0631+519, J0659+1629, J0818-2613, J0924+0219, J1433+6007, and J1817+2729. Among these systems, six do not have previously measured velocity dispersions; for the other two, our measurements are consistent with previously reported values. Additionally, we have measured the previously unknown redshifts of the deflectors in J0818-2613 and J1817+2729 to be $0.866 \pm 0.002$ and $0.408 \pm 0.002$, respectively. △ Less

Submitted 16 February, 2023; v1 submitted 28 September, 2022; originally announced September 2022.

Comments: 13 pages, 6 figures, 3 tables; accepted in A&A

Journal ref: A&A 672, A20 (2023)

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

STRIDES: Automated uniform models for 30 quadruply imaged quasars

Authors: T. Schmidt, T. Treu, S. Birrer, A. J. Shajib, C. Lemon, M. Millon, D. Sluse, A. Agnello, T. Anguita, M. W. Auger-Williams, R. G. McMahon, V. Motta, P. Schechter, C. Spiniello, I. Kayo, F. Courbin, S. Ertl, C. D. Fassnacht, J. A. Frieman, A. More, S. Schuldt, S. H. Suyu, M. Aguena, F. Andrade-Oliveira, J. Annis , et al. (46 additional authors not shown)

Abstract: Gravitational time delays provide a powerful one step measurement of $H_0$, independent of all other probes. One key ingredient in time delay cosmography are high accuracy lens models. Those are currently expensive to obtain, both, in terms of computing and investigator time (10$^{5-6}$ CPU hours and $\sim$ 0.5-1 year, respectively). Major improvements in modeling speed are therefore necessary to… ▽ More Gravitational time delays provide a powerful one step measurement of $H_0$, independent of all other probes. One key ingredient in time delay cosmography are high accuracy lens models. Those are currently expensive to obtain, both, in terms of computing and investigator time (10$^{5-6}$ CPU hours and $\sim$ 0.5-1 year, respectively). Major improvements in modeling speed are therefore necessary to exploit the large number of lenses that are forecast to be discovered over the current decade. In order to bypass this roadblock, building on the work by Shajib et al. (2019), we develop an automated modeling pipeline and apply it to a sample of 30 quadruply imaged quasars and one lensed compact galaxy, observed by the Hubble Space Telescope in multiple bands. Our automated pipeline can derive models for 30/31 lenses with few hours of human time and <100 CPU hours of computing time for a typical system. For each lens, we provide measurements of key parameters and predictions of magnification as well as time delays for the multiple images. We characterize the cosmography-readiness of our models using the stability of differences in Fermat potential (proportional to time delay) w.r.t. modeling choices. We find that for 10/30 lenses our models are cosmography or nearly cosmography grade (<3% and 3-5% variations). For 6/30 lenses the models are close to cosmography grade (5-10%). These results are based on informative priors and will need to be confirmed by further analysis. However, they are also likely to improve by extending the pipeline modeling sequence and options. In conclusion, we show that uniform cosmography grade modeling of large strong lens samples is within reach. △ Less

Submitted 9 June, 2022; originally announced June 2022.

Comments: 40 pages, 24 figures, 11 tables. Submitted to MNRAS

Detecting low-mass haloes with strong gravitational lensing I: the effect of data quality and lensing configuration

Authors: Giulia Despali, Simona Vegetti, Simon D. M. White, Devon M. Powell, Hannah R. Stacey, Christopher D. Fassnacht, Francesca Rizzo, Wolfgang Enzi

Abstract: This paper aims to quantify how the lowest halo mass that can be detected with galaxy-galaxy strong gravitational lensing depends on the quality of the observations and the characteristics of the observed lens systems. Using simulated data, we measure the lowest detectable NFW mass at each location of the lens plane, in the form of detailed \emph{sensitivity maps}. In summary, we find that: (i) th… ▽ More This paper aims to quantify how the lowest halo mass that can be detected with galaxy-galaxy strong gravitational lensing depends on the quality of the observations and the characteristics of the observed lens systems. Using simulated data, we measure the lowest detectable NFW mass at each location of the lens plane, in the form of detailed \emph{sensitivity maps}. In summary, we find that: (i) the lowest detectable mass $M_{\rm low}$ decreases linearly as the signal-to-noise ratio (SNR) increases and the sensitive area is larger when we decrease the noise; (ii) a moderate increase in angular resolution (0.07" vs 0.09") and pixel scale (0.01" vs 0.04") improves the sensitivity by on average 0.25 dex in halo mass, with more significant improvement around the most sensitive regions; (iii) the sensitivity to low-mass objects is largest for bright and complex lensed galaxies located inside the caustic curves and lensed into larger Einstein rings (i.e $r_{E}\geq1.0"$). We find that for the sensitive mock images considered in this work, the minimum mass that we can detect at the redshift of the lens lies between $1.5\times10^{8}$ and $3\times10^{9}M_{\odot}$. We derive analytic relations between $M_{\rm low}$, the SNR and resolution and discuss the impact of the lensing configuration and source structure. Our results start to fill the gap between approximate predictions and real data and demonstrate the challenging nature of calculating precise forecasts for gravitational imaging. In light of our findings, we discuss possible strategies for designing strong lensing surveys and the prospects for HST, Keck, ALMA, Euclid and other future observations. △ Less

Submitted 4 December, 2021; v1 submitted 16 November, 2021; originally announced November 2021.

Comments: 17 pages, 11 figures, accepted for publication in MNRAS. Comments welcome

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

Point-spread function reconstruction of adaptive-optics imaging: Meeting the astrometric requirements for time-delay cosmography

Authors: Geoff C. -F. Chen, Tommaso Treu, Christopher D. Fassnacht, Sam Ragland, Thomas Schmidt, Sherry H. Suyu

Abstract: Astrometric precision and knowledge of the point spread function are key ingredients for a wide range of astrophysical studies including time-delay cosmography in which strongly lensed quasar systems are used to determine the Hubble constant and other cosmological parameters. Astrometric uncertainty on the positions of the multiply-imaged point sources contributes to the overall uncertainty in inf… ▽ More Astrometric precision and knowledge of the point spread function are key ingredients for a wide range of astrophysical studies including time-delay cosmography in which strongly lensed quasar systems are used to determine the Hubble constant and other cosmological parameters. Astrometric uncertainty on the positions of the multiply-imaged point sources contributes to the overall uncertainty in inferred distances and therefore the Hubble constant. Similarly, knowledge of the wings of the points spread function (PSF) is necessary to disentangle light from the background sources and the foreground deflector. We analyze adaptive optics (AO) images of the strong lens system J0659+1629 obtained with the W. M. Keck Observatory using the laser guide star AO system. We show that by using a reconstructed point spread function we can i) obtain astrometric precision of $< 1$ milliarcsecond (mas), which is more than sufficient for time-delay cosmography; and ii) subtract all point-like images resulting in residuals consistent with the noise level. The method we have developed is not limited to strong lensing, and is generally applicable to a wide range of scientific cases that have multiple point sources nearby. △ Less

Submitted 24 June, 2021; v1 submitted 17 June, 2021; originally announced June 2021.

Comments: 7 pages, 3 figures

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

An Optical Observational Cluster Mass Function at $z\sim1$ with the ORELSE Survey

Authors: D. Hung, B. C. Lemaux, R. R. Gal, A. R. Tomczak, L. M. Lubin, O. Cucciati, D. Pelliccia, L. Shen, O. Le Fèvre, G. Zamorani, P-F. Wu, D. D. Kocevski, C. D. Fassnacht, G. K. Squires

Abstract: We present a new mass function of galaxy clusters and groups using optical/near-infrared wavelength spectroscopic and photometric data from the Observations of Redshift Evolution in Large-Scale Environments (ORELSE) survey. At $z\sim1$, cluster mass function studies are rare regardless of wavelength and have never been attempted from an optical/near-infrared perspective. This work serves as a proo… ▽ More We present a new mass function of galaxy clusters and groups using optical/near-infrared wavelength spectroscopic and photometric data from the Observations of Redshift Evolution in Large-Scale Environments (ORELSE) survey. At $z\sim1$, cluster mass function studies are rare regardless of wavelength and have never been attempted from an optical/near-infrared perspective. This work serves as a proof of concept that $z\sim1$ cluster mass functions are achievable without supplemental X-ray or Sunyaev-Zel'dovich (SZ) data. Measurements of the cluster mass function provide important contraints on cosmological parameters and are complementary to other probes. With ORELSE, a new cluster finding technique based on Voronoi tessellation Monte-Carlo (VMC) mapping, and rigorous purity and completeness testing, we have obtained $\sim$240 galaxy overdensity candidates in the redshift range $0.55<z<1.37$ at a mass range of $13.6<\log(M/M_{\odot})<14.8$. This mass range is comparable to existing optical cluster mass function studies for the local universe. Our candidate numbers vary based on the choice of multiple input parameters related to detection and characterization in our cluster finding algorithm, which we incorporated into the mass function analysis through a Monte-Carlo scheme. We find cosmological constraints on the matter density of $Ω_{m} = 0.250^{+0.104}_{-0.099}$ and on the amplitude of fluctuations of $σ_{8} = 1.150^{+0.260}_{-0.163}$. While our $Ω_{m}$ value is close to concordance, our $σ_{8}$ value is $\sim2σ$ higher because of the inflated observed number densities compared to theoretical mass function models owing to how our survey targeted overdense regions. With Euclid and several other large, unbiased optical surveys on the horizon, VMC mapping will enable optical/NIR cluster cosmology at redshifts much higher than what has been possible before. △ Less

Submitted 6 January, 2021; originally announced January 2021.

Comments: 13 pages, 8 figures, submitted to MNRAS

TDCOSMO VI: Distance Measurements in Time-delay Cosmography under the Mass-sheet transformation

Authors: Geoff C. -F. Chen, Christopher D. Fassnacht, Sherry H. Suyu, Akın Yıldırım, Eiichiro Komatsu, Jose Luis Bernal

Abstract: Time-delay cosmography with gravitationally lensed quasars plays an important role in anchoring the absolute distance scale and hence measuring the Hubble constant, $H_{0}$, independent of traditional distance ladder methodology. A current potential limitation of time delay distance measurements is the "mass-sheet transformation" (MST) which leaves the lensed imaging unchanged but changes the dist… ▽ More Time-delay cosmography with gravitationally lensed quasars plays an important role in anchoring the absolute distance scale and hence measuring the Hubble constant, $H_{0}$, independent of traditional distance ladder methodology. A current potential limitation of time delay distance measurements is the "mass-sheet transformation" (MST) which leaves the lensed imaging unchanged but changes the distance measurements and the derived value of $H_0$. In this work we show that the standard method of addressing the MST in time delay cosmography, through a combination of high-resolution imaging and the measurement of the stellar velocity dispersion of the lensing galaxy, depends on the assumption that the ratio, $D_{\rm s}/D_{\rm ds}$, of angular diameter distances to the background quasar and between the lensing galaxy and the quasar can be constrained. This is typically achieved through the assumption of a particular cosmological model. Previous work (TDCOSMO IV) addressed the mass-sheet degeneracy and derived $H_{0}$ under the assumption of $Λ$CDM model. In this paper we show that the mass sheet degeneracy can be broken without relying on a specific cosmological model by combining lensing with relative distance indicators such as supernovae type Ia and baryon acoustic oscillations, which constrain the shape of the expansion history and hence $D_{\rm s}/D_{\rm ds}$. With this approach, we demonstrate that the mass-sheet degeneracy can be constrained in a cosmological-model-independent way, and hence model-independent distance measurements in time-delay cosmography under mass-sheet transformations can be obtained. △ Less

Submitted 11 November, 2020; originally announced November 2020.

Comments: 11 pages, 8 figures

Journal ref: A&A 652, A7 (2021)

High-resolution imaging follow-up of doubly imaged quasars

Authors: Anowar J. Shajib, Eden Molina, Adriano Agnello, Peter R. Williams, Simon Birrer, Tommaso Treu, Christopher D. Fassnacht, Takahiro Morishita, Louis Abramson, Paul L. Schechter, Lutz Wisotzki

Abstract: We report upon three years of follow-up and confirmation of doubly imaged quasar lenses through imaging campaigns from 2016-2018 with the Near-Infrared Camera2 (NIRC2) on the W. M. Keck Observatory. A sample of 57 quasar lens candidates are imaged in adaptive-optics-assisted or seeing-limited $K^\prime$-band observations. Out of these 57 candidates, 15 are confirmed as lenses. We form a sample of… ▽ More We report upon three years of follow-up and confirmation of doubly imaged quasar lenses through imaging campaigns from 2016-2018 with the Near-Infrared Camera2 (NIRC2) on the W. M. Keck Observatory. A sample of 57 quasar lens candidates are imaged in adaptive-optics-assisted or seeing-limited $K^\prime$-band observations. Out of these 57 candidates, 15 are confirmed as lenses. We form a sample of 20 lenses adding in a number of previously-known lenses that were imaged with NIRC2 in 2013-14 as part of a pilot study. By modelling these 20 lenses, we obtain $K^\prime$-band relative photometry and astrometry of the quasar images and the lens galaxy. We also provide the lens properties and predicted time delays to aid planning of follow-up observations necessary for various astrophysical applications, e.g., spectroscopic follow-up to obtain the deflector redshifts for the newly confirmed systems. We compare the departure of the observed flux ratios from the smooth-model predictions between doubly and quadruply imaged quasar systems. We find that the departure is consistent between these two types of lenses if the modelling uncertainty is comparable. △ Less

Submitted 27 May, 2021; v1 submitted 3 November, 2020; originally announced November 2020.

Comments: 11 pages, 8 figures, 5 tables. This version: accepted to MNRAS

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.

TDCOSMO IV: Hierarchical time-delay cosmography -- joint inference of the Hubble constant and galaxy density profiles

Authors: S. Birrer, A. J. Shajib, A. Galan, M. Millon, T. Treu, A. Agnello, M. Auger, G. C. -F. Chen, L. Christensen, T. Collett, F. Courbin, C. D. Fassnacht, L. V. E. Koopmans, P. J. Marshall, J. -W. Park, C. E. Rusu, D. Sluse, C. Spiniello, S. H. Suyu, S. Wagner-Carena, K. C. Wong, M. Barnabè, A. S. Bolton, O. Czoske, X. Ding , et al. (2 additional authors not shown)

Abstract: The H0LiCOW collaboration inferred via gravitational lensing time delays a Hubble constant $H_0=73.3^{+1.7}_{-1.8}$ km s$^{-1}{\rm Mpc}^{-1}$, describing deflector mass density profiles by either a power-law or stars plus standard dark matter halos. The mass-sheet transform (MST) that leaves the lensing observables unchanged is considered the dominant source of residual uncertainty in $H_0$. We qu… ▽ More The H0LiCOW collaboration inferred via gravitational lensing time delays a Hubble constant $H_0=73.3^{+1.7}_{-1.8}$ km s$^{-1}{\rm Mpc}^{-1}$, describing deflector mass density profiles by either a power-law or stars plus standard dark matter halos. The mass-sheet transform (MST) that leaves the lensing observables unchanged is considered the dominant source of residual uncertainty in $H_0$. We quantify any potential effect of the MST with flexible mass models that are maximally degenerate with H0. Our calculation is based on a new hierarchical approach in which the MST is only constrained by stellar kinematics. The approach is validated on hydrodynamically simulated lenses. We apply the method to the TDCOSMO sample of 7 lenses (6 from H0LiCOW) and measure $H_0=74.5^{+5.6}_{-6.1}$ km s$^{-1}{\rm Mpc}^{-1}$. In order to further constrain the deflector mass profiles, we then add imaging and spectroscopy for 33 strong gravitational lenses from the SLACS sample. For 9 of the SLAC lenses we use resolved kinematics to constrain the stellar anisotropy. From the joint analysis of the TDCOSMO+SLACS sample, we measure $H_0=67.4^{+4.1}_{-3.2}$ km s$^{-1}{\rm Mpc}^{-1}$, assuming that the TDCOSMO and SLACS galaxies are drawn from the same parent population. The blind H0LiCOW, TDCOSMO-only and TDCOSMO+SLACS analyses are in mutual statistical agreement. The TDCOSMO+SLACS analysis prefers marginally shallower mass profiles than H0LiCOW or TDCOSMO-only. While our new analysis does not statistically invalidate the mass profile assumptions by H0LiCOW, and thus their $H_0$ measurement relying on those, it demonstrates the importance of understanding the mass density profile of elliptical galaxies. The uncertainties on $H_0$ derived in this paper can be reduced by physical or observational priors on the form of the mass profile, or by additional data, chiefly spatially resolved kinematics of lens galaxies. △ Less

Submitted 19 December, 2020; v1 submitted 6 July, 2020; originally announced July 2020.

Comments: accepted by A&A. Full analysis available at https://github.com/TDCOSMO/hierarchy_analysis_2020_public updated permanent analysis script links

Journal ref: A&A 643, A165 (2020)

TDCOSMO II: 6 new time delays in lensed quasars from high-cadence monitoring at the MPIA 2.2m telescope

Authors: M. Millon, F. Courbin, V. Bonvin, E. Buckley-Geer, C. D. Fassnacht, J. Frieman, P. J. Marshall, S. H. Suyu, T. Treu, T. Anguita, V. Motta, A. Agnello, J. H. H. Chan, D. C. -Y Chao, M. Chijani, D. Gilman, K. Gilmore, C. Lemon, J. R. Lucey, A. Melo, E. Paic, K. Rojas, D. Sluse, P. R. Williams, A. Hempel , et al. (3 additional authors not shown)

Abstract: We present six new time-delay measurements obtained from $R_c$-band monitoring data acquired at the Max Planck Institute for Astrophysics (MPIA) 2.2 m telescope at La Silla observatory between October 2016 and February 2020. The lensed quasars HE 0047-1756, WG 0214-2105, DES 0407-5006, 2M 1134-2103, PSJ 1606-2333 and DES 2325-5229 were observed almost daily at high signal-to-noise ratio to obtain… ▽ More We present six new time-delay measurements obtained from $R_c$-band monitoring data acquired at the Max Planck Institute for Astrophysics (MPIA) 2.2 m telescope at La Silla observatory between October 2016 and February 2020. The lensed quasars HE 0047-1756, WG 0214-2105, DES 0407-5006, 2M 1134-2103, PSJ 1606-2333 and DES 2325-5229 were observed almost daily at high signal-to-noise ratio to obtain high-quality light curves where we can record fast and small-amplitude variations of the quasars. We measured time delays between all pairs of multiple images with only one or two seasons of monitoring with the exception of the time delays relative to image D of PSJ 1606-2333. The most precise estimate was obtained for the delay between image A and image B of DES 0407-5006, where $τ_{AB} = -128.4^{+3.5}_{-3.8}$ d (2.8% precision) including systematics due to extrinsic variability in the light curves. For HE 0047-1756, we combined our high-cadence data with measurements from decade-long light curves from previous COSMOGRAIL campaigns, and reach a precision of 0.9 d on the final measurement. The present work demonstrates the feasibility of measuring time delays in lensed quasars in only one or two seasons, provided high signal-to-noise ratio data are obtained at a cadence close to daily. △ Less

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

Comments: 14 pages, 9 figures, 4 Tables, published in A&A

Journal ref: A&A 642, A193 (2020)

Time Delay Lens Modelling Challenge

Authors: X. Ding, T. Treu, S. Birrer, G. C. -F. Chen, J. Coles, P. Denzel, M. Frigo A. Galan, P. J. Marshall, M. Millon, A. More, A. J. Shajib, D. Sluse, H. Tak, D. Xu, M. W. Auger, V. Bonvin, H. Chand, F. Courbin, G. Despali, C. D. Fassnacht, D. Gilman, S. Hilbert, S. R. Kumar, Y. -Y. Lin, J. W. Park , et al. (4 additional authors not shown)

Abstract: In recent years, breakthroughs in methods and data have enabled gravitational time delays to emerge as a very powerful tool to measure the Hubble constant $H_0$. However, published state-of-the-art analyses require of order 1 year of expert investigator time and up to a million hours of computing time per system. Furthermore, as precision improves, it is crucial to identify and mitigate systematic… ▽ More In recent years, breakthroughs in methods and data have enabled gravitational time delays to emerge as a very powerful tool to measure the Hubble constant $H_0$. However, published state-of-the-art analyses require of order 1 year of expert investigator time and up to a million hours of computing time per system. Furthermore, as precision improves, it is crucial to identify and mitigate systematic uncertainties. With this time delay lens modelling challenge we aim to assess the level of precision and accuracy of the modelling techniques that are currently fast enough to handle of order 50 lenses, via the blind analysis of simulated datasets. The results in Rung 1 and Rung 2 show that methods that use only the point source positions tend to have lower precision ($10 - 20\%$) while remaining accurate. In Rung 2, the methods that exploit the full information of the imaging and kinematic datasets can recover $H_0$ within the target accuracy ($ |A| < 2\%$) and precision ($< 6\%$ per system), even in the presence of poorly known point spread function and complex source morphology. A post-unblinding analysis of Rung 3 showed the numerical precision of the ray-traced cosmological simulations to be insufficient to test lens modelling methodology at the percent level, making the results difficult to interpret. A new challenge with improved simulations is needed to make further progress in the investigation of systematic uncertainties. For completeness, we present the Rung 3 results in an appendix, and use them to discuss various approaches to mitigating against similar subtle data generation effects in future blind challenges. △ Less

Submitted 17 February, 2021; v1 submitted 15 June, 2020; originally announced June 2020.

Comments: 23 pages, 12 figures, 6 tables, MNRAS accepted

Journal ref: MNRAS, (2021), 503, 1096

H0LiCOW XI. A weak lensing measurement of the external convergence in the field of the lensed quasar B1608+656 using HST and Subaru deep imaging

Authors: O. Tihhonova, F. Courbin, D. Harvey, S. Hilbert, A. Peel, C. E. Rusu, C. D. Fassnacht, V. Bonvin, P. J. Marshall, G. Meylan, D. Sluse, S. H. Suyu, T. Treu, K. C. Wong

Abstract: We investigate the environment and line of sight of the H0LiCOW lens B1608+656 using Subaru Suprime-Cam and the Hubble Space Telescope (HST) to perform a weak lensing analysis. We compare three different methods to reconstruct the mass map of the field, i.e. the standard Kaiser-Squires inversion coupled with inpainting and Gaussian or wavelet filtering, and $\tt{Glimpse}$ a method based on sparse… ▽ More We investigate the environment and line of sight of the H0LiCOW lens B1608+656 using Subaru Suprime-Cam and the Hubble Space Telescope (HST) to perform a weak lensing analysis. We compare three different methods to reconstruct the mass map of the field, i.e. the standard Kaiser-Squires inversion coupled with inpainting and Gaussian or wavelet filtering, and $\tt{Glimpse}$ a method based on sparse regularization of the shear field. We find no substantial difference between the 2D mass reconstructions, but we find that the ground-based data is less sensitive to small-scale structures than the space-based observations. Marginalising over the results obtained with all the reconstruction techniques applied to the two available HST filters F606W and F814W, we estimate the external convergence, $κ_{\rm ext}$ at the position of B1608+656 is $κ_{\rm ext} = 0.11^{+0.06}_{-0.04}$, where the error bars corresponds respectively to the 16th and 84th quartiles. This result is compatible with previous estimates using the number-counts technique, suggesting that B1608+656 resides in an over-dense line of sight, but with a completely different technique. Using our mass reconstructions, we also compare the convergence at the position of several groups of galaxies in the field of B1608+656 with the mass measurements using various analytical mass profiles, and find that the weak lensing results favor truncated halo models. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: Accepted MNRAS

STRIDES: Spectroscopic and photometric characterization of the environment and effects of mass along the line of sight to the gravitational lenses DES J0408-5354 and WGD 2038-4008

Authors: E. J. Buckley-Geer, H. Lin, C. Rusu, J. Poh, A. Palmese, A. Agnello, L. Christensen, J. Frieman, A. J. Shajib, T. Treu, T. Collett, S. Birrer, T. Anguita, C. D. Fassnacht, G. Meylan, S. Mukherjee, K. C. Wong, M. Aguena, S. Allam, S. Avila, E. Bertin, S. Bhargava, D. Brooks, A. Carnero Rosell, M. Carrasco Kind , et al. (43 additional authors not shown)

Abstract: In time-delay cosmography, three of the key ingredients are 1) determining the velocity dispersion of the lensing galaxy, 2) identifying galaxies and groups along the line of sight with sufficient proximity and mass to be included in the mass model, and 3) estimating the external convergence $κ_\mathrm{ext}$ from less massive structures that are not included in the mass model. We present results o… ▽ More In time-delay cosmography, three of the key ingredients are 1) determining the velocity dispersion of the lensing galaxy, 2) identifying galaxies and groups along the line of sight with sufficient proximity and mass to be included in the mass model, and 3) estimating the external convergence $κ_\mathrm{ext}$ from less massive structures that are not included in the mass model. We present results on all three of these ingredients for two time-delay lensed quasar systems, DES J0408-5354 and WGD 2038-4008. We use the Gemini, Magellan and VLT telescopes to obtain spectra to both measure the stellar velocity dispersions of the main lensing galaxies and to identify the line-of-sight galaxies in these systems. Next, we identify 10 groups in DES J0408-5354 and 2 groups in WGD 2038-4008using a group-finding algorithm. We then identify the most significant galaxy and galaxy-group perturbers using the "flexion shift" criterion. We determine the probability distribution function of the external convergence $κ_\mathrm{ext}$ for both of these systems based on our spectroscopy and on the DES-only multiband wide-field observations. Using weighted galaxy counts, calibrated based on the Millennium Simulation, we find that DES J0408-5354 is located in a significantly underdense environment, leading to a tight (width $\sim3\%$), negative-value $κ_\mathrm{ext}$ distribution. On the other hand, WGD 2038-4008 is located in an environment of close to unit density, and its low source redshift results in a much tighter $κ_\mathrm{ext}$ of $\sim1\%$, as long as no external shear constraints are imposed. △ Less

Submitted 30 March, 2020; v1 submitted 26 March, 2020; originally announced March 2020.

Comments: 40 pages, 18 figures Added missing author to the author list

The STRong lensing Insights into the Dark Energy Survey (STRIDES) 2017/2018 follow-up campaign: Discovery of 10 lensed quasars and 10 quasar pairs

Authors: C. Lemon, M. W. Auger, R. McMahon, T. Anguita, Y. Apostolovski, G. C. -F. Chen, C. D. Fassnacht, A. Melo, V. Motta, A. Shajib, T. Treu, A. Agnello, E. Buckley-Geer, P. L. Schechter, S. Birrer, T. Collett, F. Courbin, C. E. Rusu, T. M. C. Abbott, S. Allam, J. Annis, S. Avila, E. Bertin, D. Brooks, D. L. Burke , et al. (43 additional authors not shown)

Abstract: We report the results of the STRong lensing Insights from the Dark Energy Survey (STRIDES) follow-up campaign of the late 2017/early 2018 season. We obtained spectra of 65 lensed quasar candidates either with EFOSC2 on the NTT or ESI on Keck, which confirm 10 new gravitationally lensed quasars and 10 quasar pairs with similar spectra, but which do not show a lensing galaxy in DES images. Eight len… ▽ More We report the results of the STRong lensing Insights from the Dark Energy Survey (STRIDES) follow-up campaign of the late 2017/early 2018 season. We obtained spectra of 65 lensed quasar candidates either with EFOSC2 on the NTT or ESI on Keck, which confirm 10 new gravitationally lensed quasars and 10 quasar pairs with similar spectra, but which do not show a lensing galaxy in DES images. Eight lensed quasars are doubly imaged with source redshifts between 0.99 and 2.90, one is triply imaged by a group (DESJ0345-2545, $z=1.68$), and one is quadruply imaged (quad: DESJ0053-2012, $z=3.8$). Singular isothermal ellipsoid models for the doubles, based on high-resolution imaging from SAMI on SOAR or NIRC2 on Keck, give total magnifications between 3.2 and 5.6, and Einstein radii between 0.49 and 1.97 arcseconds. After spectroscopic follow-up, we extract multi-epoch $grizY$ photometry of confirmed lensed quasars and contaminant quasar+star pairs from the first 4 years of DES data using parametric multi-band modelling, and compare variability in each system's components. By measuring the reduced $χ^2$ associated with fitting all epochs to the same magnitude, we find a simple cut on the less variable component that retains all confirmed lensed quasars, while removing 94 per cent of contaminant systems with stellar components. Based on our spectroscopic follow-up, this variability information can improve selection of lensed quasars and quasar pairs from 34-45 per cent to 51-70 per cent, with the majority of remaining contaminants being compact star-forming galaxies. Using mock lensed quasar lightcurves we demonstrate that selection based only on variability will over-represent the quad fraction by 10 per cent over a complete DES magnitude-limited sample (excluding microlensing differences), explained by the magnification bias and hence lower luminosity (more variable) sources in quads. △ Less

Submitted 19 December, 2019; originally announced December 2019.

Comments: 22 pages, 16 figures, submitted to MNRAS

TDCOSMO. I. An exploration of systematic uncertainties in the inference of $H_0$ from time-delay cosmography

Authors: M. Millon, A. Galan, F. Courbin, T. Treu, S. H. Suyu, X. Ding, S. Birrer, G. C. -F. Chen, A. J. Shajib, D. Sluse, K. C. Wong, A. Agnello, M. W. Auger, E. J. Buckley-Geer, J. H. H. Chan, T. Collett, C. D. Fassnacht, S. Hilbert, L. V. E. Koopmans, V. Motta, S. Mukherjee, C. E. Rusu, A. Sonnenfeld, C. Spiniello, L. Van de Vyvere

Abstract: Time-delay cosmography of lensed quasars has achieved 2.4% precision on the measurement of the Hubble constant, $H_0$. As part of an ongoing effort to uncover and control systematic uncertainties, we investigate three potential sources: 1- stellar kinematics, 2- line-of-sight effects, and 3- the deflector mass model. To meet this goal in a quantitative way, we reproduced the H0LiCOW/SHARP/STRIDES… ▽ More Time-delay cosmography of lensed quasars has achieved 2.4% precision on the measurement of the Hubble constant, $H_0$. As part of an ongoing effort to uncover and control systematic uncertainties, we investigate three potential sources: 1- stellar kinematics, 2- line-of-sight effects, and 3- the deflector mass model. To meet this goal in a quantitative way, we reproduced the H0LiCOW/SHARP/STRIDES (hereafter TDCOSMO) procedures on a set of real and simulated data, and we find the following. First, stellar kinematics cannot be a dominant source of error or bias since we find that a systematic change of 10% of measured velocity dispersion leads to only a 0.7% shift on $H_0$ from the seven lenses analyzed by TDCOSMO. Second, we find no bias to arise from incorrect estimation of the line-of-sight effects. Third, we show that elliptical composite (stars + dark matter halo), power-law, and cored power-law mass profiles have the flexibility to yield a broad range in $H_0$ values. However, the TDCOSMO procedures that model the data with both composite and power-law mass profiles are informative. If the models agree, as we observe in real systems owing to the "bulge-halo" conspiracy, $H_0$ is recovered precisely and accurately by both models. If the two models disagree, as in the case of some pathological models illustrated here, the TDCOSMO procedure either discriminates between them through the goodness of fit, or it accounts for the discrepancy in the final error bars provided by the analysis. This conclusion is consistent with a reanalysis of six of the TDCOSMO (real) lenses: the composite model yields $74.0^{+1.7}_{-1.8}$ $km.s^{-1}.Mpc^{-1}$, while the power-law model yields $H_0=74.2^{+1.6}_{-1.6}$ $km.s^{-1}.Mpc^{-1}$. In conclusion, we find no evidence of bias or errors larger than the current statistical uncertainties reported by TDCOSMO. △ Less

Submitted 22 July, 2020; v1 submitted 17 December, 2019; originally announced December 2019.

Comments: 19 pages, 9 figures, 4 tables, published in A&A

Journal ref: A&A 639, A101 (2020)

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.

STRIDES: a 3.9 per cent measurement of the Hubble constant from the strong lens system DES J0408-5354

Authors: A. J. Shajib, S. Birrer, T. Treu, A. Agnello, E. J. Buckley-Geer, J. H. H. Chan, L. Christensen, C. Lemon, H. Lin, M. Millon, J. Poh, C. E. Rusu, D. Sluse, C. Spiniello, G. C. -F. Chen, T. Collett, F. Courbin, C. D. Fassnacht, J. Frieman, A. Galan, D. Gilman, A. More, T. Anguita, M. W. Auger, V. Bonvin , et al. (66 additional authors not shown)

Abstract: We present a blind time-delay cosmographic analysis for the lens system DES J0408$-$5354. This system is extraordinary for the presence of two sets of multiple images at different redshifts, which provide the opportunity to obtain more information at the cost of increased modelling complexity with respect to previously analyzed systems. We perform detailed modelling of the mass distribution for th… ▽ More We present a blind time-delay cosmographic analysis for the lens system DES J0408$-$5354. This system is extraordinary for the presence of two sets of multiple images at different redshifts, which provide the opportunity to obtain more information at the cost of increased modelling complexity with respect to previously analyzed systems. We perform detailed modelling of the mass distribution for this lens system using three band Hubble Space Telescope imaging. We combine the measured time delays, line-of-sight central velocity dispersion of the deflector, and statistically constrained external convergence with our lens models to estimate two cosmological distances. We measure the "effective" time-delay distance corresponding to the redshifts of the deflector and the lensed quasar $D_{Δt}^{\rm eff}=3382^{+146}_{-115}$ Mpc and the angular diameter distance to the deflector $D_{\rm d}=1711^{+376}_{-280}$ Mpc, with covariance between the two distances. From these constraints on the cosmological distances, we infer the Hubble constant $H_0 = 74.2^{+2.7}_{-3.0}$ km s$^{-1}$ Mpc$^{-1}$ assuming a flat $Λ$CDM cosmology and a uniform prior for $Ω_{\rm m}$ as $Ω_{\rm m} \sim \mathcal{U}(0.05, 0.5)$. This measurement gives the most precise constraint on $H_0$ to date from a single lens. Our measurement is consistent with that obtained from the previous sample of six lenses analyzed by the $H_0$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration. It is also consistent with measurements of $H_0$ based on the local distance ladder, reinforcing the tension with the inference from early Universe probes, for example, with 2.2$σ$ discrepancy from the cosmic microwave background measurement. △ Less

Submitted 14 April, 2020; v1 submitted 14 October, 2019; originally announced October 2019.

Comments: 32 pages, 17 figures, 10 tables. Accepted by MNRAS

Cosmic dissonance: new physics or systematics behind a short sound horizon?

Authors: Nikki Arendse, Radosław J. Wojtak, Adriano Agnello, Geoff C. -F. Chen, Christopher D. Fassnacht, Dominique Sluse, Stefan Hilbert, Martin Millon, Vivien Bonvin, Kenneth C. Wong, Frédéric Courbin, Sherry H. Suyu, Simon Birrer, Tommaso Treu, Leon V. E. Koopmans

Abstract: Persistent tension between low-redshift observations and the Cosmic Microwave Background radiation (CMB), in terms of two fundamental distance scales set by the sound horizon $r_d$ and the Hubble constant $H_0$, suggests new physics beyond the Standard Model or residual systematics. We examine recently updated distance calibrations from Cepheids, gravitational lensing time-delay observations, and… ▽ More Persistent tension between low-redshift observations and the Cosmic Microwave Background radiation (CMB), in terms of two fundamental distance scales set by the sound horizon $r_d$ and the Hubble constant $H_0$, suggests new physics beyond the Standard Model or residual systematics. We examine recently updated distance calibrations from Cepheids, gravitational lensing time-delay observations, and the Tip of the Red Giant Branch. Calibrating the Baryon Acoustic Oscillations (BAO) and Type Ia supernovae with combinations of the distance indicators, we obtain a joint and self-consistent measurement of $H_0$ and $r_d$ at low redshift, independent of cosmological models and CMB inference. In an attempt to alleviate the tension between late-time and CMB-based measurements, we consider four extensions of the standard $Λ$CDM model. The sound horizon from our different measurements is $r_d=(137\pm3^{stat.}\pm2^{syst.})$~Mpc. Depending on the adopted distance indicators, the $combined$ tension in $H_0$ and $r_d$ ranges between 2.3 and 5.1 $σ$. We find that modifications of $Λ$CDM that change the physics after recombination fail to solve the problem, for the reason that they only resolve the tension in $H_0$, while the tension in $r_d$ remains unchanged. Pre-recombination extensions (with early dark energy or the effective number of neutrinos $\rm{N}_{\rm{eff}}=3.24 \pm 0.16$) are allowed by the data, unless the calibration from Cepheids is included. Results from time-delay lenses are consistent with those from distance-ladder calibrations and point to a discrepancy between absolute distance scales measured from the CMB (assuming the standard cosmological model) and late-time observations. New proposals to resolve this tension should be examined with respect to reconciling not only the Hubble constant but also the sound horizon derived from the CMB and other cosmological probes. △ Less

Submitted 13 May, 2020; v1 submitted 17 September, 2019; originally announced September 2019.

Comments: Accepted for publication in A&A; 14 pages, 8 figures, 4 tables

Journal ref: A&A 639, A57 (2020)

A measurement of the Hubble constant from angular diameter distances to two gravitational lenses

Authors: Inh Jee, Sherry Suyu, Eiichiro Komatsu, Christopher D. Fassnacht, Stefan Hilbert, Léon V. E. Koopmans

Abstract: The local expansion rate of the Universe is parametrized by the Hubble constant, $H_0$, the ratio between recession velocity and distance. Different techniques lead to inconsistent estimates of $H_0$. Observations of Type Ia supernovae (SNe) can be used to measure $H_0$, but this requires an external calibrator to convert relative distances to absolute ones. We use the angular diameter distance to… ▽ More The local expansion rate of the Universe is parametrized by the Hubble constant, $H_0$, the ratio between recession velocity and distance. Different techniques lead to inconsistent estimates of $H_0$. Observations of Type Ia supernovae (SNe) can be used to measure $H_0$, but this requires an external calibrator to convert relative distances to absolute ones. We use the angular diameter distance to strong gravitational lenses as a suitable calibrator, which is only weakly sensitive to cosmological assumptions. We determine the angular diameter distances to two gravitational lenses, $810^{+160}_{-130}$ and $1230^{+180}_{-150}$~Mpc, at redshifts of $z=0.295$ and $0.6304$. Using these absolute distances to calibrate 740 previously-measured relative distances to SNe, we measure the Hubble constant to be $H_0=82.4^{+8.4}_{-8.3} ~{\rm km\,s^{-1}\,Mpc^{-1}}$. △ Less

Submitted 14 September, 2019; originally announced September 2019.

Comments: This paper presents the measurements of angular diameter distances to two time-delay lenses, and the Hubble constant derived only from these two distances and the JLA supernova sample. One of the distance measurements is further used for the cosmological inference in the H0LiCOW XIII paper (arxiv:1907.04869). Published in Science

Journal ref: Science 13 Sep 2019:Vol. 365, Issue 6458, pp. 1134-1138

Double dark matter vision: twice the number of compact-source lenses with narrow-line lensing and the WFC3 grism

Authors: A. M. Nierenberg, D. Gilman, T. Treu, G. Brammer, S. Birrer, L. Moustakas, A. Agnello, T. Anguita, C. D. Fassnacht, V. Motta, A. H. G. Peter, D. Sluse

Abstract: The magnifications of compact-source lenses are extremely sensitive to the presence of low mass dark matter halos along the entire sight line from the source to the observer. Traditionally, the study of dark matter structure in compact-source strong gravitational lenses has been limited to radio-loud systems, as the radio emission is extended and thus unaffected by microlensing which can mimic the… ▽ More The magnifications of compact-source lenses are extremely sensitive to the presence of low mass dark matter halos along the entire sight line from the source to the observer. Traditionally, the study of dark matter structure in compact-source strong gravitational lenses has been limited to radio-loud systems, as the radio emission is extended and thus unaffected by microlensing which can mimic the signal of dark matter structure. An alternate approach is to measure quasar nuclear-narrow line emission, which is free from microlensing and present in virtually all quasar lenses. In this paper, we double the number of systems which can be used for gravitational lensing analyses by presenting measurements of narrow-line emission from a sample of 8 quadruply imaged quasar lens systems, WGD J0405-3308, HS 0810+2554, RX J0911+0551, SDSS J1330+1810, PS J1606-2333, WFI 2026-4536, WFI 2033-4723 and WGD J2038-4008. We describe our updated grism spectral modelling pipeline, which we use to measure narrow-line fluxes with uncertainties of 2-10\%, presented here. We fit the lensed image positions with smooth mass models and demonstrate that these models fail to produce the observed distribution of image fluxes over the entire sample of lenses. Furthermore, typical deviations are larger than those expected from macromodel uncertainties. This discrepancy indicates the presence of perturbations caused by small-scale dark matter structure. The interpretation of this result in terms of dark matter models is presented in a companion paper. △ Less

Submitted 16 December, 2019; v1 submitted 17 August, 2019; originally announced August 2019.

Comments: Accepted for publication in MNRAS

H0LiCOW XIII. A 2.4% measurement of $H_{0}$ from lensed quasars: $5.3σ$ tension between early and late-Universe probes

Authors: Kenneth C. Wong, Sherry H. Suyu, Geoff C. -F. Chen, Cristian E. Rusu, Martin Millon, Dominique Sluse, Vivien Bonvin, Christopher D. Fassnacht, Stefan Taubenberger, Matthew W. Auger, Simon Birrer, James H. H. Chan, Frederic Courbin, Stefan Hilbert, Olga Tihhonova, Tommaso Treu, Adriano Agnello, Xuheng Ding, Inh Jee, Eiichiro Komatsu, Anowar J. Shajib, Alessandro Sonnenfeld, Roger D. Blandford, Leon V. E. Koopmans, Philip J. Marshall , et al. (1 additional authors not shown)

Abstract: We present a measurement of the Hubble constant ($H_{0}$) and other cosmological parameters from a joint analysis of six gravitationally lensed quasars with measured time delays. All lenses except the first are analyzed blindly with respect to the cosmological parameters. In a flat $Λ$CDM cosmology, we find $H_{0} = 73.3_{-1.8}^{+1.7}$, a 2.4% precision measurement, in agreement with local measure… ▽ More We present a measurement of the Hubble constant ($H_{0}$) and other cosmological parameters from a joint analysis of six gravitationally lensed quasars with measured time delays. All lenses except the first are analyzed blindly with respect to the cosmological parameters. In a flat $Λ$CDM cosmology, we find $H_{0} = 73.3_{-1.8}^{+1.7}$, a 2.4% precision measurement, in agreement with local measurements of $H_{0}$ from type Ia supernovae calibrated by the distance ladder, but in $3.1σ$ tension with $Planck$ observations of the cosmic microwave background (CMB). This method is completely independent of both the supernovae and CMB analyses. A combination of time-delay cosmography and the distance ladder results is in $5.3σ$ tension with $Planck$ CMB determinations of $H_{0}$ in flat $Λ$CDM. We compute Bayes factors to verify that all lenses give statistically consistent results, showing that we are not underestimating our uncertainties and are able to control our systematics. We explore extensions to flat $Λ$CDM using constraints from time-delay cosmography alone, as well as combinations with other cosmological probes, including CMB observations from $Planck$, baryon acoustic oscillations, and type Ia supernovae. Time-delay cosmography improves the precision of the other probes, demonstrating the strong complementarity. Allowing for spatial curvature does not resolve the tension with $Planck$. Using the distance constraints from time-delay cosmography to anchor the type Ia supernova distance scale, we reduce the sensitivity of our $H_0$ inference to cosmological model assumptions. For six different cosmological models, our combined inference on $H_{0}$ ranges from $\sim73$-$78~\mathrm{km~s^{-1}~Mpc^{-1}}$, which is consistent with the local distance ladder constraints. △ Less

Submitted 5 November, 2019; v1 submitted 10 July, 2019; originally announced July 2019.

Comments: Accepted for publication in MNRAS; 23 pages, 13 figures, 8 tables

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

H0LiCOW XII. Lens mass model of WFI2033-4723 and blind measurement of its time-delay distance and $H_0$

Authors: Cristian E. Rusu, Kenneth C. Wong, Vivien Bonvin, Dominique Sluse, Sherry H. Suyu, Christopher D. Fassnacht, James H. H. Chan, Stefan Hilbert, Matthew W. Auger, Alessandro Sonnenfeld, Simon Birrer, Frederic Courbin, Tommaso Treu, Geoff C. -F. Chen, Aleksi Halkola, Leon V. E. Koopmans, Philip J. Marshall, Anowar J. Shajib

Abstract: We present the lens mass model of the quadruply-imaged gravitationally lensed quasar WFI2033-4723, and perform a blind cosmographical analysis based on this system. Our analysis combines (1) time-delay measurements from 14 years of data obtained by the COSmological MOnitoring of GRAvItational Lenses (COSMOGRAIL) collaboration, (2) high-resolution $\textit{Hubble Space Telescope}$ imaging, (3) a me… ▽ More We present the lens mass model of the quadruply-imaged gravitationally lensed quasar WFI2033-4723, and perform a blind cosmographical analysis based on this system. Our analysis combines (1) time-delay measurements from 14 years of data obtained by the COSmological MOnitoring of GRAvItational Lenses (COSMOGRAIL) collaboration, (2) high-resolution $\textit{Hubble Space Telescope}$ imaging, (3) a measurement of the velocity dispersion of the lens galaxy based on ESO-MUSE data, and (4) multi-band, wide-field imaging and spectroscopy characterizing the lens environment. We account for all known sources of systematics, including the influence of nearby perturbers and complex line-of-sight structure, as well as the parametrization of the light and mass profiles of the lensing galaxy. After unblinding, we determine the effective time-delay distance to be $4784_{-248}^{+399}~\mathrm{Mpc}$, an average precision of $6.6\%$. This translates to a Hubble constant $H_{0} = 71.6_{-4.9}^{+3.8}~\mathrm{km~s^{-1}~Mpc^{-1}}$, assuming a flat $Λ$CDM cosmology with a uniform prior on $Ω_\mathrm{m}$ in the range [0.05, 0.5]. This work is part of the $H_0$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration, and the full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper (H0LiCOW XIII). △ Less

Submitted 8 December, 2019; v1 submitted 22 May, 2019; originally announced May 2019.

Comments: Version accepted by MNRAS. 29 pages including appendix, 17 figures, 6 tables. arXiv admin note: text overlap with arXiv:1607.01403

H0LiCOW X: Spectroscopic/imaging survey and galaxy-group identification around the strong gravitational lens system WFI2033-4723

Authors: D. Sluse, C. E. Rusu, C. D. Fassnacht, A. Sonnenfeld, J. Richard, M. W. Auger, L. Coccato, K. C. Wong, S. H. Suyu, T. Treu, A. Agnello, S. Birrer, V. Bonvin, T. Collett, F. Courbin, S. Hilbert, L. V. E. Koopmans, O. Tihhanova, P. J. Marshall, G. Meylan, A. J. Shajib, J. Annis, S. Avila, E. Bertin, D. Brooks , et al. (43 additional authors not shown)

Abstract: Galaxies and galaxy groups located along the line of sight towards gravitationally lensed quasars produce high-order perturbations of the gravitational potential at the lens position. When these perturbation are too large, they can induce a systematic error on $H_0$ of a few-percent if the lens system is used for cosmological inference and the perturbers are not explicitly accounted for in the len… ▽ More Galaxies and galaxy groups located along the line of sight towards gravitationally lensed quasars produce high-order perturbations of the gravitational potential at the lens position. When these perturbation are too large, they can induce a systematic error on $H_0$ of a few-percent if the lens system is used for cosmological inference and the perturbers are not explicitly accounted for in the lens model. In this work, we present a detailed characterization of the environment of the lens system WFI2033-4723 ($z_{\rm src} = 1.662$, $z_{\rm lens}$ = 0.6575), one of the core targets of the H0LICOW project for which we present cosmological inferences in a companion paper (Rusu et al. 2019). We use the Gemini and ESO-Very Large telescopes to measure the spectroscopic redshifts of the brightest galaxies towards the lens, and use the ESO-MUSE integral field spectrograph to measure the velocity-dispersion of the lens ($σ_{\rm {los}}= 250^{+15}_{-21}$ km/s) and of several nearby galaxies. In addition, we measure photometric redshifts and stellar masses of all galaxies down to $i < 23$ mag, mainly based on Dark Energy Survey imaging (DR1). Our new catalog, complemented with literature data, more than doubles the number of known galaxy spectroscopic redshifts in the direct vicinity of the lens, expanding to 116 (64) the number of spectroscopic redshifts for galaxies separated by less than 3 arcmin (2 arcmin) from the lens. Using the flexion-shift as a measure of the amplitude of the gravitational perturbation, we identify 2 galaxy groups and 3 galaxies that require specific attention in the lens models. The ESO MUSE data enable us to measure the velocity-dispersions of three of these galaxies. These results are essential for the cosmological inference analysis presented in Rusu et al. (2019). △ Less

Submitted 6 September, 2019; v1 submitted 21 May, 2019; originally announced May 2019.

Comments: Matches the version accepted for publication by MNRAS. Note that this paper previously appeared as H0LICOW XI

Report number: FERMILAB-PUB-19-168-AE

COSMOGRAIL XVIII: time delays of the quadruply lensed quasar WFI2033-4723

Authors: V. Bonvin, M. Millon, J. H. H. Chan, F. Courbin, C. E. Rusu, D. Sluse, S. H. Suyu, K. C. Wong, C. D. Fassnacht, P. J. Marshall, T. Treu, E. Buckley-Geer, J. Frieman, A. Hempel, S. Kim, R. Lachaume, M. Rabus, D. C. -Y. Chao, M. Chijani, D. Gilman, K. Gilmore, K. Rojas, P. Williams, T. Anguita, C. S. Kochanek , et al. (4 additional authors not shown)

Abstract: We present new measurements of the time delays of WFI2033-4723. The data sets used in this work include 14 years of data taken at the 1.2m Leonhard Euler Swiss telescope, 13 years of data from the SMARTS 1.3m telescope at Las Campanas Observatory and a single year of high-cadence and high-precision monitoring at the MPIA 2.2m telescope. The time delays measured from these different data sets, all… ▽ More We present new measurements of the time delays of WFI2033-4723. The data sets used in this work include 14 years of data taken at the 1.2m Leonhard Euler Swiss telescope, 13 years of data from the SMARTS 1.3m telescope at Las Campanas Observatory and a single year of high-cadence and high-precision monitoring at the MPIA 2.2m telescope. The time delays measured from these different data sets, all taken in the R-band, are in good agreement with each other and with previous measurements from the literature. Combining all the time-delay estimates from our data sets results in Dt_AB = 36.2-0.8+0.7 days (2.1% precision), Dt_AC = -23.3-1.4+1.2 days (5.6%) and Dt_BC = -59.4-1.3+1.3 days (2.2%). In addition, the close image pair A1-A2 of the lensed quasars can be resolved in the MPIA 2.2m data. We measure a time delay consistent with zero in this pair of images. We also explore the prior distributions of microlensing time-delay potentially affecting the cosmological time-delay measurements of WFI2033-4723. There is however no strong indication in our measurements that microlensing time delay is neither present nor absent. This work is part of a H0LiCOW series focusing on measuring the Hubble constant from WFI2033-4723. △ Less

Submitted 20 May, 2019; originally announced May 2019.

Comments: Submitted to Astronomy and Astrophysics

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

Dark Matter Science in the Era of LSST

Authors: Keith Bechtol, Alex Drlica-Wagner, Kevork N. Abazajian, Muntazir Abidi, Susmita Adhikari, Yacine Ali-Haïmoud, James Annis, Behzad Ansarinejad, Robert Armstrong, Jacobo Asorey, Carlo Baccigalupi, Arka Banerjee, Nilanjan Banik, Charles Bennett, Florian Beutler, Simeon Bird, Simon Birrer, Rahul Biswas, Andrea Biviano, Jonathan Blazek, Kimberly K. Boddy, Ana Bonaca, Julian Borrill, Sownak Bose, Jo Bovy , et al. (155 additional authors not shown)

Abstract: Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We… ▽ More Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 11 pages, 2 figures, Science Whitepaper for Astro 2020, more information at https://lsstdarkmatter.github.io

Probing the Fundamental Nature of Dark Matter with the Large Synoptic Survey Telescope

Authors: Alex Drlica-Wagner, Yao-Yuan Mao, Susmita Adhikari, Robert Armstrong, Arka Banerjee, Nilanjan Banik, Keith Bechtol, Simeon Bird, Kimberly K. Boddy, Ana Bonaca, Jo Bovy, Matthew R. Buckley, Esra Bulbul, Chihway Chang, George Chapline, Johann Cohen-Tanugi, Alessandro Cuoco, Francis-Yan Cyr-Racine, William A. Dawson, Ana Díaz Rivero, Cora Dvorkin, Denis Erkal, Christopher D. Fassnacht, Juan García-Bellido, Maurizio Giannotti , et al. (75 additional authors not shown)

Abstract: Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We d… ▽ More Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ . △ Less

Submitted 24 April, 2019; v1 submitted 4 February, 2019; originally announced February 2019.

Comments: 96 pages, 22 figures, 1 table

Report number: FERMILAB-PUB-19-048-A-AE

Possible Evidence of the Radio AGN Quenching of Neighboring Galaxies at z $\sim$ 1

Authors: Lu Shen, Adam R. Tomczak, Brian C. Lemaux, Debora Pelliccia, Lori M. Lubin, Neal A. Miller, Serena Perrotta, Christopher D. Fassnacht, Robert H. Becker, Roy R. Gal, Po-Feng Wu, Gordon Squires

Abstract: Using 57 Radio Active Galactic nuclei (RAGN) at 0.55 $\leq$ z $\leq$ 1.3 drawn from five fields of the Observations of Redshift Evolution in Large Scale Environments (ORELSE) survey, we study the effect of injection of energy from outbursts of RAGN on their spectroscopically-confirmed neighboring galaxies (SNGs). We observe an elevated fraction of quenched neighbors (fq) within 500 kpc projected r… ▽ More Using 57 Radio Active Galactic nuclei (RAGN) at 0.55 $\leq$ z $\leq$ 1.3 drawn from five fields of the Observations of Redshift Evolution in Large Scale Environments (ORELSE) survey, we study the effect of injection of energy from outbursts of RAGN on their spectroscopically-confirmed neighboring galaxies (SNGs). We observe an elevated fraction of quenched neighbors (fq) within 500 kpc projected radius of RAGN in the most dense local environments compared to those of non-RAGN control samples matched to the RAGN population in colour, stellar mass, and local environment at 2$σ$ significance. Further analyses show that there are offsets at similar significance between fqs of RAGN-SNGs and the appropriate control samples for galaxies specifically in cluster environments and those hosted by most massive cluster galaxies, which tentatively suggests that some negative feedback from the RAGN is occurring in these dense environments. In addition, we find that the median radio power of RAGN increases with increasing local overdensity, an effect which may lend itself to the quenching of neighboring galaxies. Furthermore, we find that, in the highest local overdensities, the fq of the sub-sample of lower stellar mass RAGN-SNGs is larger than that of the higher stellar mass RAGN-SNGs sub-sample, which indicates a more pronounced effect from RAGN on lower stellar mass galaxies. We propose a scenario in which RAGN residing within clusters might heat the intracluster medium (ICM) affecting both in situ star formation and any inflowing gas that remains in their neighboring galaxies. △ Less

Submitted 11 December, 2018; originally announced December 2018.

Comments: submitted to MNRAS, 15 pages, 2 tables and 6 figures, comments welcome

H0LiCOW - IX. Cosmographic analysis of the doubly imaged quasar SDSS 1206+4332 and a new measurement of the Hubble constant

Authors: S. Birrer, T. Treu, C. E. Rusu, V. Bonvin, C. D. Fassnacht, J. H. H. Chan, A. Agnello, A. J. Shajib, G. C. -F. Chen, M. Auger, F. Courbin, S. Hilbert, D. Sluse, S. H. Suyu, K. C. Wong, P. Marshall, B. C. Lemaux, G. Meylan

Abstract: We present a blind time-delay strong lensing (TDSL) cosmographic analysis of the doubly imaged quasar SDSS 1206+4332. We combine the relative time delay between the quasar images, Hubble Space Telescope imaging, the Keck stellar velocity dispersion of the lensing galaxy, and wide-field photometric and spectroscopic data of the field to constrain two angular diameter distance relations. The combine… ▽ More We present a blind time-delay strong lensing (TDSL) cosmographic analysis of the doubly imaged quasar SDSS 1206+4332. We combine the relative time delay between the quasar images, Hubble Space Telescope imaging, the Keck stellar velocity dispersion of the lensing galaxy, and wide-field photometric and spectroscopic data of the field to constrain two angular diameter distance relations. The combined analysis is performed by forward modelling the individual data sets through a Bayesian hierarchical framework, and it is kept blind until the very end to prevent experimenter bias. After unblinding, the inferred distances imply a Hubble constant $H_0 = 68.8^{+5.4}_{-5.1}$ kms$^{-1}$Mpc$^{-1}$, assuming a flat Lambda cold dark matter cosmology with uniform prior on $Ω_{\rm m}$ in [0.05, 0.5]. The precision of our cosmographic measurement with the doubly imaged quasar SDSS 1206+4332 is comparable with those of quadruply imaged quasars and opens the path to perform on selected doubles the same analysis as anticipated for quads. Our analysis is based on a completely independent lensing code than our previous three H0LiCOW systems and the new measurement is fully consistent with those. We provide the analysis scripts paired with the publicly available software to facilitate independent analysis. The consistency between blind measurements with independent codes provides an important sanity check on lens modelling systematics. By combining the likelihoods of the four systems under the same prior, we obtain $H_0 = 72.5^{+2.1}_{-2.3}$kms$^{-1}$Mpc$^{-1}$. This measurement is independent of the distance ladder and other cosmological probes. △ Less

Submitted 24 May, 2019; v1 submitted 4 September, 2018; originally announced September 2018.

Comments: 30 pages, 17 figures, MNRAS published likelihood available here: http://shsuyu.github.io/H0LiCOW/site/notebooks/H0_from_lenses.html, all modeling and analysis scripts available upon request

The STRong lensing Insights into the Dark Energy Survey (STRIDES) 2016 follow-up campaign. I. Overview and classification of candidates selected by two techniques

Authors: T. Treu, A. Agnello, M. A. Baumer, S. Birrer, E. J. Buckley-Geer, F. Courbin, Y. J. Kim, H. Lim, P. J. Marshall, B. Nord, P. L. Schechter, P. R. Sivakumar, L. E. Abramson, T. Anguita, Y. Apostolovski, M. W. Auger, J. H. H. Chan, G. C. F. Chen, T. E. Collett, C. D. Fassnacht, J. -W. Hsueh, C. Lemon, R. G. McMahon, V. Motta, F. Ostrovski , et al. (57 additional authors not shown)

Abstract: The primary goals of the STRong lensing Insights into the Dark Energy Survey (STRIDES) collaboration are to measure the dark energy equation of state parameter and the free streaming length of dark matter. To this aim, STRIDES is discovering strongly lensed quasars in the imaging data of the Dark Energy Survey and following them up to measure time delays, high resolution imaging, and spectroscopy… ▽ More The primary goals of the STRong lensing Insights into the Dark Energy Survey (STRIDES) collaboration are to measure the dark energy equation of state parameter and the free streaming length of dark matter. To this aim, STRIDES is discovering strongly lensed quasars in the imaging data of the Dark Energy Survey and following them up to measure time delays, high resolution imaging, and spectroscopy sufficient to construct accurate lens models. In this paper, we first present forecasts for STRIDES. Then, we describe the STRIDES classification scheme, and give an overview of the Fall 2016 follow-up campaign. We continue by detailing the results of two selection methods, the Outlier Selection Technique and a morphological algorithm, and presenting lens models of a system, which could possibly be a lensed quasar in an unusual configuration. We conclude with the summary statistics of the Fall 2016 campaign. Including searches presented in companion papers (Anguita et al.; Ostrovski et al.), STRIDES followed up 117 targets identifying 7 new strongly lensed systems, and 7 nearly identical quasars (NIQs), which could be confirmed as lenses by the detection of the lens galaxy. 76 candidates were rejected and 27 remain otherwise inconclusive, for a success rate in the range 6-35\%. This rate is comparable to that of previous searches like SQLS even though the parent dataset of STRIDES is purely photometric and our selection of candidates cannot rely on spectroscopic information. △ Less

Submitted 5 September, 2018; v1 submitted 14 August, 2018; originally announced August 2018.

Comments: (MNRAS in press; color scheme of Figures 1-4 improved with respect to v1)

Is every strong lens model unhappy in its own way? Uniform modelling of a sample of 13 quadruply+ imaged quasars

Authors: A. J. Shajib, S. Birrer, T. Treu, M. W. Auger, A. Agnello, T. Anguita, E. J. Buckley-Geer, J. H. H. Chan, T. E. Collett, F. Courbin, C. D. Fassnacht, J. Frieman, I. Kayo, C. Lemon, H. Lin, P. J. Marshall, R. McMahon, A. More, N. D. Morgan, V. Motta, M. Oguri, F. Ostrovski, C. E. Rusu, P. L. Schechter, T. Shanks , et al. (48 additional authors not shown)

Abstract: Strong-gravitational lens systems with quadruply-imaged quasars (quads) are unique probes to address several fundamental problems in cosmology and astrophysics. Although they are intrinsically very rare, ongoing and planned wide-field deep-sky surveys are set to discover thousands of such systems in the next decade. It is thus paramount to devise a general framework to model strong-lens systems to… ▽ More Strong-gravitational lens systems with quadruply-imaged quasars (quads) are unique probes to address several fundamental problems in cosmology and astrophysics. Although they are intrinsically very rare, ongoing and planned wide-field deep-sky surveys are set to discover thousands of such systems in the next decade. It is thus paramount to devise a general framework to model strong-lens systems to cope with this large influx without being limited by expert investigator time. We propose such a general modelling framework (implemented with the publicly available software Lenstronomy) and apply it to uniformly model three-band Hubble Space Telescope Wide Field Camera 3 images of 13 quads. This is the largest uniformly modelled sample of quads to date and paves the way for a variety of studies. To illustrate the scientific content of the sample, we investigate the alignment between the mass and light distribution in the deflectors. The position angles of these distributions are well-aligned, except when there is strong external shear. However, we find no correlation between the ellipticity of the light and mass distributions. We also show that the observed flux-ratios between the images depart significantly from the predictions of simple smooth models. The departures are strongest in the bluest band, consistent with microlensing being the dominant cause in addition to millilensing. Future papers will exploit this rich dataset in combination with ground based spectroscopy and time delays to determine quantities such as the Hubble constant, the free streaming length of dark matter, and the normalization of the initial stellar mass function. △ Less

Submitted 3 September, 2020; v1 submitted 24 July, 2018; originally announced July 2018.

Comments: 22 pages, 9 figures, 8 tables, Accepted to MNRAS. This version: updated Table 4 from the erratum

Journal ref: MNRAS 483 (2019) 5649

SHARP - V. Modelling gravitationally-lensed radio arcs imaged with global VLBI observations

Authors: C. Spingola, J. P. McKean, M. W. Auger, C. D. Fassnacht, L. V. E. Koopmans, D. J. Lagattuta, S. Vegetti

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

Submitted 15 July, 2018; originally announced July 2018.

Comments: 15 pages, 8 figures; published by MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Vol. 478, Issue 4, 2018

The STRong lensing Insights into the Dark Energy Survey (STRIDES) 2016 follow-up campaign. II. New quasar lenses from double component fitting

Authors: T. Anguita, P. L. Schechter, N. Kuropatkin, N. D. Morgan, F. Ostrovski, L. E. Abramson, A. Agnello, Y. Apostolovski, C. D. Fassnacht, J. W. Hsueh, V. Motta, K. Rojas, C. E. Rusu, T. Treu, P. Williams, M. Auger, E. Buckley-Geer, H. Lin, R. McMahon, T. M. C. Abbott, S. Allam, J. Annis, R. A. Bernstein, E. Bertin, D. Brooks , et al. (34 additional authors not shown)

Abstract: We report upon the follow up of 34 candidate lensed quasars found in the Dark Energy Survey using NTT-EFOSC, Magellan-IMACS, KECK-ESI and SOAR-SAMI. These candidates were selected by a combination of double component fitting, morphological assessment and color analysis. Most systems followed up are indeed composed of at least one quasar image and 13 with two or more quasar images: two lenses, four… ▽ More We report upon the follow up of 34 candidate lensed quasars found in the Dark Energy Survey using NTT-EFOSC, Magellan-IMACS, KECK-ESI and SOAR-SAMI. These candidates were selected by a combination of double component fitting, morphological assessment and color analysis. Most systems followed up are indeed composed of at least one quasar image and 13 with two or more quasar images: two lenses, four projected binaries and seven Nearly Identical Quasar Pairs (NIQs). The two systems confirmed as genuine gravitationally lensed quasars are one quadruple at $z_s=1.713$ and one double at $z_s=1.515$. Lens modeling of these two systems reveals that both systems require very little contribution from the environment to reproduce the image configuration. Nevertheless, small flux anomalies can be observed in one of the images of the quad. Further observations of 9 inconclusive systems (including 7 NIQs) will allow to confirm (or not) their gravitational lens nature. △ Less

Submitted 30 May, 2018; originally announced May 2018.

Comments: 13 pages, 8 figures. MNRAS submitted (05/24). Approved by the DES Collaboration Wide Review

Evaluating Tests of Virialization and Substructure Using Galaxy Clusters in the ORELSE Survey

Authors: N. Rumbaugh, B. C. Lemaux, A. R. Tomczak, L. Shen, D. Pelliccia, L. M. Lubin, D. D. Kocevski, P. -F. Wu, R. R. Gal, S. Mei, C. D. Fassnacht, G. K. Squires

Abstract: We evaluated the effectiveness of different indicators of cluster virialization using 12 large-scale structures in the ORELSE survey spanning from $0.7<z<1.3$. We located diffuse X-ray emission from 16 galaxy clusters using Chandra observations. We studied the properties of these clusters and their members, using Chandra data in conjunction with optical and near-IR imaging and spectroscopy. We mea… ▽ More We evaluated the effectiveness of different indicators of cluster virialization using 12 large-scale structures in the ORELSE survey spanning from $0.7<z<1.3$. We located diffuse X-ray emission from 16 galaxy clusters using Chandra observations. We studied the properties of these clusters and their members, using Chandra data in conjunction with optical and near-IR imaging and spectroscopy. We measured X-ray luminosities and gas temperatures of each cluster, as well as velocity dispersions of their member galaxies. We compared these results to scaling relations derived from virialized clusters, finding significant offsets of up to 3-4$σ$ for some clusters, which could indicate they are disturbed or still forming. We explored if other properties of the clusters correlated with these offsets by performing a set of tests of virialization and substructure on our sample, including Dressler-Schectman tests, power ratios, analyses of the velocity distributions of galaxy populations, and centroiding differences. For comparison to a wide range of studies, we used two sets of tests: ones that did and did not use spectral energy distribution fitting to obtain rest-frame colours, stellar masses, and photometric redshifts of galaxies. Our results indicated that the difference between the stellar mass or light mean-weighted center and the X-ray center, as well as the projected offset of the most-massive/brightest cluster galaxy from other cluster centroids had the strongest correlations with scaling relation offsets, implying they are the most robust indicators of cluster virialization and can be used for this purpose when X-ray data is insufficiently deep for reliable $L_X$ and $T_X$ measurements. △ Less

Submitted 31 May, 2018; v1 submitted 1 May, 2018; originally announced May 2018.

Comments: 24 pages, 6 figures, 10 tables, accepted to MNRAS, v2 fixes a rendering issue and typos in Table 9, v3 corrects further typos in Table 1 and elsewhere

Constraining the microlensing effect on time delays with new time-delay prediction model in $H_{0}$ measurements

Authors: Geoff C. -F. Chen, James H. H. Chan, Vivien Bonvin, Christopher D. Fassnacht, Karina Rojas, Martin Millon, Fred Courbin, Sherry H. Suyu, Kenneth C. Wong, Dominique Sluse, Tommaso Treu, Anowar J. Shajib, Jen-Wei Hsueh, David J. Lagattuta, Leon V. E. Koopmans, Simona Vegetti, John P. McKean

Abstract: Time-delay strong lensing provides a unique way to directly measure the Hubble constant ($H_{0}$). The precision of the $H_{0}$ measurement depends on the uncertainties in the time-delay measurements, the mass distribution of the main deflector(s), and the mass distribution along the line of sight. Tie and Kochanek (2018) have proposed a new microlensing effect on time delays based on differential… ▽ More Time-delay strong lensing provides a unique way to directly measure the Hubble constant ($H_{0}$). The precision of the $H_{0}$ measurement depends on the uncertainties in the time-delay measurements, the mass distribution of the main deflector(s), and the mass distribution along the line of sight. Tie and Kochanek (2018) have proposed a new microlensing effect on time delays based on differential magnification of the coherent accretion disc variability of the lensed quasar. If real, this effect could significantly broaden the uncertainty on the time delay measurements by up to $30\%$ for lens systems such as PG1115+080, which have relatively short time delays and monitoring over several different epochs. In this paper we develop a new technique that uses the time-delay ratios and simulated microlensing maps within a Bayesian framework in order to limit the allowed combinations of microlensing delays and thus to lessen the uncertainties due to the proposed effect. We show that, under the assumption of Tie and Kochanek (2018), the uncertainty on the time-delay distance ($D_{Δt}$, which is proportional to 1/$H_{0}$) of short time-delay ($\sim18$ days) lens, PG1115+080, increases from $\sim7\%$ to $\sim10\%$ by simultaneously fitting the three time-delay measurements from the three different datasets across twenty years, while in the case of long time-delay ($\sim90$ days) lens, the microlensing effect on time delays is negligible as the uncertainty on $D_{Δt}$ of RXJ1131-1231 only increases from $\sim2.5\%$ to $\sim2.6\%$. △ Less

Submitted 28 August, 2018; v1 submitted 25 April, 2018; originally announced April 2018.