Search

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., 24 pages, 15 figures, 3 tables. Accepted in MNRAS

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 $\Lambda$CDM model with fixed $\Omega_{\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 $\Delta t_{\rm BA}h\hat{\sigma}_{v}^{-2}=6.89\substack{+0.8\\-0.7}$ days, $\Delta t_{\rm CA}h\hat{\sigma}_{v}^{-2}=10.7\substack{+1.6\\-1.2}$ days, and $\Delta t_{\rm DA}h\hat{\sigma}_{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{\sigma}_{v}=\sigma^{\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., Comment: 11 pages, 7 figures

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. In this work, for the first time, we break the MSD 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 allows us to robustly estimate $H_0$, effectively implementing a maximally flexible mass model. Following a blind analysis, we estimate the angular diameter distance to the lens galaxy $D_{\rm d} = 865_{-81}^{+85}$ Mpc and the time-delay distance $D_{\Delta t} = 2180_{-271}^{+472}$ Mpc, giving $H_0 = 77.1_{-7.1}^{+7.3}$ km s$^{-1}$ Mpc$^{-1}$ - for a flat $\Lambda$ 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., Comment: To be submitted to A&A. 21 pages, 22 figures, 1 table

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., 23 pages, 12 figures, 6 tables, MNRAS accepted

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 $\Lambda$CDM cosmology with a uniform prior on $\Omega_\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)., Comment: Version accepted by MNRAS. 29 pages including appendix, 17 figures, 6 tables. arXiv admin note: text overlap with arXiv:1607.01403

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.

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., Comment: 17 pages, 11 figures, accepted for publication in MNRAS. Comments welcome

The paper ‘Is every strong lens model unhappy in its own way? Uniform modelling of a sample of 13 quadruply+ imaged quasars’ was published in MNRAS, 483, 4, 5649–5671 (2019). The coordinate values of the image positions in table 4 were wrongly printed due to a clerical error. At a later stage of writing the manuscript, we have changed the zero-point definition of the lens coordinate systems, but the relative image positions were not accounted for this change of definition while printing out table 4. We provide the updated Table 4 below. This error does not impact any other results of the paper in any way, except for the table itself. We thank Collin Werner and Paul Schechter for helping us identify this error.

We present a measurement of the Hubble constant (H0) and other cosmological parameters from a joint analysis of six gravitationally lensed quasars with measured time delays. All lenses except the first are analysed blindly with respect to the cosmological parameters. In a flat Λ cold dark matter (ΛCDM) cosmology, we find $H_{0} = 73.3_{-1.8}^{+1.7}~\mathrm{km~s^{-1}~Mpc^{-1}}$, a $2.4{{\ \rm per\ cent}}$ precision measurement, in agreement with local measurements of H0 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 H0 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 H0 inference to cosmological model assumptions. For six different cosmological models, our combined inference on H0 ranges from ∼73 to 78 km s−1 Mpc−1, which is consistent with the local distance ladder constraints.

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, Comment: 13 pages, 8 figures, submitted to MNRAS

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-$\alpha$ forest, and the number of luminous satellites in the Milky Way. We derive an upper limit of $\lambda_{\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 $\lambda_{\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., Comment: 16 pages, 3 Figures, minor updates of the results

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 (L 1.4 GHz ∼ 1025 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 the Keck/Multi-Object Spectrometer for Infra-Red Exploration 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-Very Large Telescope Deep Survey 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 ∼2.6 × 1015 M ⊙ and contains several overdensity peaks. Both RAGNs are hosted by very bright and massive galaxies, while their hosts show extreme differences in color, indicating that they are of 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 have already 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.

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., 7 pages, 3 figures

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$\sigma$. 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.

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., Comment: 19 pages, 9 figures, 4 tables, published in A&A

One of the main challenges in using high redshift active galactic nuclei to study the correlations between the mass of the supermassive Black Hole (MBH) and the properties of their active host galaxies is instrumental resolution. Strong lensing magnification effectively increases instrumental resolution and thus helps to address this challenge. In this work, we study eight strongly lensed active galactic nuclei (AGN) with deep Hubble Space Telescope imaging, using the lens modelling code Lenstronomy to reconstruct the image of the source. Using the reconstructed brightness of the host galaxy, we infer the host galaxy stellar mass based on stellar population models. MBH are estimated from broad emission lines using standard methods. Our results are in good agreement with recent work based on non-lensed AGN, demonstrating the potential of using strongly lensed AGNs to extend the study of the correlations to higher redshifts. At the moment, the sample size of lensed AGN is small and thus they provide mostly a consistency check on systematic errors related to resolution for the non-lensed AGN. However, the number of known lensed AGN is expected to increase dramatically in the next few years, through dedicated searches in ground and space based wide field surveys, and they may become a key diagnostic of black hole and galaxy co-evolution., 12 pages, 4 figures, 3 tables. MNRAS in press. Comments welcome

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 ${\chi}^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., Comment: 22 pages, 16 figures, submitted to MNRAS

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., Comment: accepted by A&A. Full analysis available at https://github.com/TDCOSMO/hierarchy_analysis_2020_public updated permanent analysis script links

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 $\Lambda$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., Comment: 11 pages, 8 figures

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., Accepted for publication in MNRAS

Full author list: 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, R McMahon, G Meylan, K C Wong, T M C Abbott, J Annis, S Avila, K Bechtol, D Brooks, D Brout, D L Burke, A Carnero Rosell, M Carrasco Kind, J Carretero, F J Castander, M Costanzi, L N da Costa, J De Vicente, S Desai, J P Dietrich, P Doel, A Drlica-Wagner, A E Evrard, D A Finley, B Flaugher, P Fosalba, J García-Bellido, D W Gerdes, D Gruen, R A Gruendl, J Gschwend, G Gutierrez, D L Hollowood, K Honscheid, D Huterer, D J James, T Jeltema, E Krause, N Kuropatkin, T S Li, M Lima, N MacCrann, M A G Maia, J L Marshall, P Melchior, R Miquel, R L C Ogando, A Palmese, F Paz-Chinchón, A A Plazas, A K Romer, A Roodman, M Sako, E Sanchez, B Santiago, V Scarpine, M Schubnell, D Scolnic, S Serrano, I Sevilla-Noarbe, M Smith, M Soares-Santos, E Suchyta, G Tarle, D Thomas, A R Walker, Y Zhang, 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 analysed 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 eff=3382-115+146 Mpc and the angular diameter distance to the deflector Dd = 1711-280+376 Mpc, with covariance between the two distances. From these constraints on the cosmological distances, we infer the Hubble constant H0= 74.2-3.0+2.7 km s-1 Mpc-1 assuming a flat ΛCDM cosmology and a uniform prior for ωm as \Omega m ∼ \mathcal {U(0.05, 0.5). This measurement gives the most precise constraint on H0 to date from a single lens. Our measurement is consistent with that obtained from the previous sample of six lenses analysed by the H0 Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration. It is also consistent with measurements of H0 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., TT acknowledges support by the Packard Foundation through a Packard Research fellowship and by the National Science Foundation through NSF grants AST-1714953 and AST-1906976. This project is partly funded by the Danish council for independent research under the project ‘Fundamentals of Dark Matter Structures’, DFF - 6108-00470. AA was supported by a grant from VILLUM FONDEN (project number 16599). JHHC, MM, CL, DS, FC, and AG acknowledge support from the Swiss National Science Foundation (SNSF) and by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (COSMICLENS: grant agreement No 787866). CS has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie actions grant agreement No 664931. GCFC acknowledges support from the Ministry of Education in Taiwan via Government Scholarship to Study Abroad (GSSA). CDF and GCFC acknowledge support for this work from the National Science Foundation under Grant No AST-1715611. TA acknowledges support from Proyecto FONDECYT N○ 1190335. KCW was supported by World Premier International Research Center Initiative (WPI Initiative), MEXT, Japan. This work used computational and storage services associated with the Hoffman2 Shared Cluster provided by UCLA Institute for Digital Research and Education’s Research Technology Group. This work used the Extreme Science and Engineering Discovery Environment (XSEDE) through allocation TG-AST190038, which is supported by National Science Foundation grant number ACI-1548562 (Towns et al. 2014). Specifically, this work used the Comet and Oasis system at the San Diego Supercomputer Center (SDSC), and the Bridges system, which is supported by NSF award number ACI-1445606, at the Pittsburgh Supercomputing Center (PSC, Nystrom et al. 2015). AJS thanks Smadar Naoz for providing access to additional computing nodes on the Hoffman2 Shared Cluster and Khalid Jawed for assisting with additional computational resources from the Structures–Computer Interaction Laboratory at UCLA. The DES data management system is supported by the National Science Foundation under Grant Numbers AST-1138766 and AST-1536171. The DES participants from Spanish institutions are partially supported by MINECO under grants AYA2015-71825, ESP2015-66861, FPA2015-68048, SEV-2016-0588, SEV-2016-0597, and MDM-2015-0509, some of which include ERDF funds from the European Union. IFAE is partially funded by the CERCA program of the Generalitat de Catalunya. Research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Program (FP7/2007-2013) including ERC grant agreements 240672, 291329, and 306478. We acknowledge support from the Australian Research Council Centre of Excellence for All-sky Astrophysics (CAASTRO), through project number CE110001020, and the Brazilian Instituto Nacional de Ciência e Tecnologia (INCT) e-Universe (CNPq grant 465376/2014-2).

The DES Collaboration: et al., 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 κ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 quad 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 two groups in WGD 2038–4008 using 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 κ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 ∼3 per cent⁠), negative-value κ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 κext of ∼1 per cent⁠, as long as no external shear constraints are imposed., JP would like to thank Gourav Khullar for his help and insightful discussions that helped improve the analysis in this paper. This work made use of computing resources and support provided by the Research Computing Center at the University of Chicago. JP is supported in part by the Kavli Institute for Cosmological Physics at the University of Chicago through grant NSF PHY-1125897 and an endowment from the Kavli Foundation and its founder Fred Kavli. AJS was supported by the National Aeronautics and Space Administration (NASA) through the Space Telescope Science Institute (STScI) grant HST-GO-15320. AJS was also supported by the Dissertation Year Fellowship from the University of California, Los Angeles (UCLA) graduate division. TA acknowledges support from Proyecto FONDECYT N: 1190335. This work was supported by World Premier International Research Center Initiative (WPI Initiative), MEXT, Japan. CDF acknowledges support for this work from the National Science Foundation under Grant No. AST-1715611. SM acknowledges the funding from the European Research Council (ERC) under the EUs Horizon 2020 research and innovation program (COSMICLENS; grant agreement No. 787886). Funding for the DES Projects has been provided by the U.S. Department of Energy, the U.S. National Science Foundation, the Ministry of Science and Education of Spain, the Science and Technology Facilities Council of the United Kingdom, the Higher Education Funding Council for England, the National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign, the Kavli Institute of Cosmological Physics at the University of Chicago, the Center for Cosmology and Astro-Particle Physics at the Ohio State University, the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University, Financiadora de Estudos e Projetos, Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Científico e Tecnológico and the Ministério da Ciência, Tecnologia e Inovação, the Deutsche Forschungsgemeinschaft and the Collaborating Institutions in the Dark Energy Survey. The Collaborating Institutions are Argonne National Laboratory, the University of California at Santa Cruz, the University of Cambridge, Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas-Madrid, the University of Chicago, University College London, the DES-Brazil Consortium, the University of Edinburgh, the Eidgenössische Technische Hochschule (ETH) Zürich, Fermi National Accelerator Laboratory, the University of Illinois at Urbana-Champaign, the Institut de Ciències de l’Espai (IEEC/CSIC), the Institut de Física d’Altes Energies, Lawrence Berkeley National Laboratory, the Ludwig-Maximilians Universität München and the associated Excellence Cluster Universe, the University of Michigan, the National Optical Astronomy Observatory, the University of Nottingham, The Ohio State University, the University of Pennsylvania, the University of Portsmouth, SLAC National Accelerator Laboratory, Stanford University, the University of Sussex, Texas A&M University, and the OzDES Membership Consortium. Based in part on observations at Cerro Tololo Inter-American Observatory, National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. The DES data management system is supported by the National Science Foundation under Grant Numbers AST-1138766 and AST-1536171. The DES participants from Spanish institutions are partially supported by MINECO under grants AYA2015-71825, ESP2015-88861, FPA2015-68048, SEV-2012-0234, SEV-2016-0597, and MDM-2015-0509, some of which include ERDF funds from the European Union. IFAE is partially funded by the CERCA program of the Generalitat de Catalunya. Research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Program (FP7/2007-2013) including ERC grant agreements 240672, 291329, and 306478. We acknowledge support from the Australian Research Council Centre of Excellence for All-sky Astrophysics (CAASTRO), through project number CE110001020. This manuscript has been authored by Fermi Research Alliance, LLC under Contract No. DE-AC02-07CH11359 with the U.S. Department of Energy, Office of Science, Office of High Energy Physics. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes.

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., Comment: 11 pages, 8 figures, 5 tables. This version: accepted to MNRAS

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 that 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, 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 of 6-35 per cent. This rate is comparable to that of previous searches like SDSS Quasar Lens Search even though the parent data set of STRIDES is purely photometric and our selection of candidates cannot rely on spectroscopic information., TT and VM acknowledge support by the David and Lucille Packard Foundation through a Packard Research Fellowship to TT. TT acknowledges support by the National Science Foundation through grants AST- 1450141 and AST-1714953. CDF and GCFC acknowledge support from the National Science Foundation through grants AST-1312329 and AST-1715611. TA and YA acknowledge support by proyecto FONDECYT 11130630 and by the Ministry for the Economy, Development, and Tourism’s Programa Inicativa Científica Milenio through grant IC 12009, awarded to The Millennium Institute of Astrophysics (MAS). FC, VB, and JC acknowledge support from the Swiss National Science Foundation. SHS thanks the Max Planck Society for support through the Max Planck Research Group

We present a weak gravitational lensing measurement of the external convergence along the line of sight to the quadruply lensed quasar HE$\,$0435$-$1223. Using deep r-band images from Subaru-Suprime-Cam we observe galaxies down to a 3$\sigma$ limiting magnitude of $\sim 26$ mags resulting in a source galaxy density of 14 galaxies / arcmin$^2$ after redshift-based cuts. Using an inpainting technique and Multi-Scale Entropy filtering algorithm, we find that the region in close proximity to the lens has an estimated external convergence of $\kappa=-0.012^{+0.020}_{-0.013}$ and is hence marginally under-dense. We also rule out the presence of any halo with a mass greater than $M_{\rm vir}=1.6\times10^{14}h^{-1}M_\odot$ (68$\%$ confidence limit). Our results, consistent with previous studies of this lens, confirm that the intervening mass along the line of sight to HE$\,$0435$-$1223 does not affect significantly the cosmological results inferred from the time delay measurements of that specific object., Comment: 14 pages, 8 figures, 2 tables

Galaxies located in the environment or on the line of sight towards gravitational lenses can significantly affect lensing observables, and can lead to systematic errors on the measurement of $H_0$ from the time-delay technique. We present the results of a systematic spectroscopic identification of the galaxies in the field of view of the lensed quasar HE0435-1223, using the W. M. Keck, Gemini and ESO-Very Large telescopes. Our new catalog triples the number of known galaxy redshifts in the vicinity of the lens, expanding to 102 the number of measured redshifts for galaxies separated by less than 3 arcmin from the lens. We complement our catalog with literature data to gather redshifts up to 15 arcmin from the lens, and search for galaxy groups or clusters projected towards HE0435-1223. We confirm that the lens is a member of a small group that includes at least 12 galaxies, and find 8 other group candidates near the line of sight of the lens. The flexion shift, namely the shift of lensed images produced by high order perturbation of the lens potential, is calculated for each galaxy/group and used to identify which objects produce the largest perturbation of the lens potential. This analysis demonstrates that i) at most three of the five brightest galaxies projected within 12 arcsec of the lens need to be explicitly used in the lens models, and ii) the groups can be treated in the lens model as an external tidal field (shear) contribution., Comment: Version revised to address referee's comments, submitted to MNRAS, 21 pages (incl. Appendix). Data associated to the paper available from the H0LICOW website www.h0licow.org

We present an exploration of the mass structure of a sample of 12 strongly lensed massive, compact early-type galaxies at redshifts $z\sim0.6$ to provide further possible evidence for their inside-out growth. We obtain new ESI/Keck spectroscopy and infer the kinematics of both lens and source galaxies, and combine these with existing photometry to construct (a) the fundamental plane (FP) of the source galaxies and (b) physical models for their dark and luminous mass structure. We find their FP to be tilted towards the virial plane relative to the local FP, and attribute this to their unusual compactness, which causes their kinematics to be totally dominated by the stellar mass as opposed to their dark matter; that their FP is nevertheless still inconsistent with the virial plane implies that both the stellar and dark structure of early-type galaxies is non-homologous. We also find the intrinsic scatter of their FP to be comparable to the local value, indicating that variations in the stellar mass structure outweight variations in the dark halo in the central regions of early-type galaxies. Finally, we show that inference on the dark halo structure -- and, in turn, the underlying physics -- is sensitive to assumptions about the stellar initial mass function (IMF), but that physically-motivated assumptions about the IMF imply haloes with sub-NFW inner density slopes, and may present further evidence for the inside-out growth of compact early-type galaxies via minor mergers and accretion., Comment: 10 pages, 3 figures, 3 tables; submitted to MNRAS

We present a new sample of strong gravitational lens systems where both the foreground lenses and background sources are early-type galaxies. Using imaging from HST/ACS and Keck/NIRC2, we model the surface brightness distributions and show that the sources form a distinct population of massive, compact galaxies at redshifts $0.4 \lesssim z \lesssim 0.7$, lying systematically below the size-mass relation of the global elliptical galaxy population at those redshifts. These may therefore represent relics of high-redshift red nuggets or their partly-evolved descendants. We exploit the magnifying effect of lensing to investigate the structural properties, stellar masses and stellar populations of these objects with a view to understanding their evolution. We model these objects parametrically and find that they generally require two S\'ersic components to properly describe their light profiles, with one more spheroidal component alongside a more envelope-like component, which is slightly more extended though still compact. This is consistent with the hypothesis of the inside-out growth of these objects via minor mergers. We also find that the sources can be characterised by red-to-blue colour gradients as a function of radius which are stronger at low redshift -- indicative of ongoing accretion -- but that their environments generally appear consistent with that of the general elliptical galaxy population, contrary to recent suggestions that these objects are predominantly associated with clusters., Comment: 21 pages; accepted for publication in MNRAS

Based on spectroscopy and multiband wide-field observations of the gravitationally lensed quasar HE 0435-1223, we determine the probability distribution function of the external convergence $\kappa_\mathrm{ext}$ for this system. We measure the under/overdensity of the line of sight towards the lens system and compare it to the average line of sight throughout the universe, determined by using the CFHTLenS as a control field. Aiming to constrain $\kappa_\mathrm{ext}$ as tightly as possible, we determine under/overdensities using various combinations of relevant informative weighing schemes for the galaxy counts, such as projected distance to the lens, redshift, and stellar mass. We then convert the measured under/overdensities into a $\kappa_\mathrm{ext}$ distribution, using ray-tracing through the Millennium Simulation. We explore several limiting magnitudes and apertures, and account for systematic and statistical uncertainties relevant to the quality of the observational data, which we further test through simulations. Our most robust estimate of $\kappa_\mathrm{ext}$ has a median value $\kappa^\mathrm{med}_\mathrm{ext} = 0.004$ and a standard deviation of $\sigma_\kappa = 0.025$. The measured $\sigma_\kappa$ corresponds to $2.5\%$ uncertainty on the time delay distance, and hence the Hubble constant $H_0$ inference from this system. The median $\kappa^\mathrm{med}_\mathrm{ext}$ value is robust to $\sim0.005$ (i.e. $\sim0.5\%$ on $H_0$) regardless of the adopted aperture radius, limiting magnitude and weighting scheme, as long as the latter incorporates galaxy number counts, the projected distance to the main lens, and a prior on the external shear obtained from mass modeling. The availability of a well-constrained $\kappa_\mathrm{ext}$ makes \hequad\ a valuable system for measuring cosmological parameters using strong gravitational lens time delays., Comment: 24 pages, 17 figures, 6 tables. Submitted to MNRAS

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., 14 pages, 7 figures, accepted for publication in MNRAS

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 ($\sigma_{\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)., Comment: Matches the version accepted for publication by MNRAS. Note that this paper previously appeared as H0LICOW XI

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

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

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$\sigma$ 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., Comment: submitted to MNRAS, 15 pages, 2 tables and 6 figures, comments welcome

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 $\Lambda$ cold dark matter cosmology with a uniform prior on $\Omega_{\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., Comment: 33 pages, 23 figures

We studied AGN activity in twelve LSSs in the ORELSE survey, at 0.65, Comment: 26 pages, 8 figures, submitted to MNRAS

We present a detection in pre-explosion Hubble Space Telescope (HST) imaging of a point source consistent with being the progenitor star of the Type IIb supernova (SN IIb) 2016gkg. Post-explosion imaging from the Keck Adaptive Optics system was used to perform relative astrometry between the Keck and HST imaging. We identify a single point source in the HST images coincident with the SN position to 0.89-sigma. The HST photometry is consistent with the progenitor star being an A0Ia star with T=9500 K and log (L/Lsun)=5.15. We find that the SN 2016gkg progenitor star appears more consistent with binary than single-star evolutionary models. In addition, early-time light curve data from SN 2016gkg revealed a rapid rise in luminosity within ~0.4 days of non-detection limits, consistent with models of the cooling phase after shock break-out. We use these data to determine an explosion date of 20.15 September 2016 and progenitor star radius of log (R/Rsun)=2.41, which agrees with photometry from the progenitor star. Our findings are also consistent with detections of other SNe IIb progenitor stars, although more luminous and bluer than most other examples., Comment: 8 pages, 5 figures, submitted to MNRAS

We present a new measurement of the Hubble Constant H0 and other cosmological parameters based on the joint analysis of three multiply imaged quasar systems with measured gravitational time delays. First, we measure the time delay of HE 0435−1223 from 13-yr light curves obtained as part of the COSMOGRAIL project. Companion papers detail the modelling of the main deflectors and line-of-sight effects, and how these data are combined to determine the time-delay distance of HE 0435−1223. Crucially, the measurements are carried out blindly with respect to cosmological parameters in order to avoid confirmation bias. We then combine the time-delay distance of HE 0435−1223 with previous measurements from systems B1608+656 and RXJ1131−1231 to create a Time Delay Strong Lensing probe (TDSL). In flat Λ cold dark matter (ΛCDM) with free matter and energy density, we find H0 =71.9+2.4-3.0 km s−1Mpc−1 andΩΛ=0.62+0.24−0.35. This measurement is completely independent of, and in agreement with, the local distance ladder measurements of H0. We explore more general cosmological models combining TDSL with other probes, illustrating its power to break degeneracies inherent to other methods. The joint constraints from TDSL and Planck are H0 = 69.2+1.4−2.2 km s−1Mpc−1, ΩΛ=0.70+0.01−0.01 and Ωk=0.003+0.004−0.006 in open ΛCDM and H0 =79.0+4.4−4.2 km s−1Mpc−1, Ωde=0.77+0.02−0.03 and w=−1.38+0.14−0.16 in flat wCDM. In combination with Planck and baryon acoustic oscillation data, when relaxing the constraints on the numbers of relativistic species we find Neff = 3.34+0.21−0.21 in NeffΛCDM and when relaxing the total mass of neutrinos we find Σmν ≤ 0.182 eV in mνΛCDM. Finally, in an open wCDM in combination with Planck and cosmic microwave background lensing, we find H0 =77.9+5.0−4.2 km s−1Mpc−1, Ωde=0.77+0.03−0.03, Ωk=−0.003+0.004−0.004 and w=−1.37+0.18−0.23.

Accurate and precise measurements of the Hubble constant are critical for testing our current standard cosmological model and revealing possibly new physics. With Hubble Space Telescope (HST) imaging, each strong gravitational lens system with measured time delays can allow one to determine the Hubble constant with an uncertainty of $\sim 7\%$. Since HST will not last forever, we explore adaptive-optics (AO) imaging as an alternative that can provide higher angular resolution than HST imaging but has a less stable point spread function (PSF) due to atmospheric distortion. To make AO imaging useful for time-delay-lens cosmography, we develop a method to extract the unknown PSF directly from the imaging of strongly lensed quasars. In a blind test with two mock data sets created with different PSFs, we are able to recover the important cosmological parameters (time-delay distance, external shear, lens mass profile slope, and total Einstein radius). Our analysis of the Keck AO image of the strong lens system RXJ1131-1231 shows that the important parameters for cosmography agree with those based on HST imaging and modeling within 1-$\sigma$ uncertainties. Most importantly, the constraint on the model time-delay distance by using AO imaging with $0.045"$resolution is tighter by $\sim 50\%$ than the constraint of time-delay distance by using HST imaging with $0.09"$when a power-law mass distribution for the lens system is adopted. Our PSF reconstruction technique is generic and applicable to data sets that have multiple nearby point sources, enabling scientific studies that require high-precision models of the PSF., Comment: 20 pages, 15 figures

We present six new time-delay measurements obtained from Rc-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.8+3.5 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.

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 $\Lambda$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 $\sigma$. We find that modifications of $\Lambda$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., Comment: Accepted for publication in A&A; 14 pages, 8 figures, 4 tables

We present time-delay estimates for the quadruply imaged quasar PG 1115+080. Our resuls are based on almost daily observations for seven months at the ESO MPIA 2.2m telescope at La Silla Observatory, reaching a signal-to-noise ratio of about 1000 per quasar image. In addition, we re-analyse existing light curves from the literature that we complete with an additional three seasons of monitoring with the Mercator telescope at La Palma Observatory. When exploring the possible source of bias we consider the so-called microlensing time delay, a potential source of systematic error so far never directly accounted for in previous time-delay publications. In fifteen years of data on PG 1115+080, we find no strong evidence of microlensing time delay. Therefore not accounting for this effect, our time-delay estimates on the individual data sets are in good agreement with each other and with the literature. Combining the data sets, we obtain the most precise time-delay estimates to date on PG 1115+080, with Dt(AB) = 8.3+1.5-1.6 days (18.7% precision), Dt(AC) = 9.9+1.1-1.1 days (11.1%) and Dt(BC) = 18.8+1.6-1.6 days (8.5%). Turning these time delays into cosmological constraints is done in a companion paper that makes use of ground-based Adaptive Optics (AO) with the Keck telescope., 15 pages, 6 figures, submitted to A&A. Comments welcome

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

We present results of a systematic search for gravitationally lensed quasars in Pan-STARRS1. Our final sample of candidates comprises of 91 systems, not including 25 rediscovered lensed quasars and quasar pairs. In the absence of spectroscopy to verify the lensing nature of the candidates, the main sources of contaminants are likely to be quasar pairs, which we consider to be a byproduct of our work, and a smaller number of quasar$+$star associations. Amongst the independently discovered quads is 2M1134$-$2103, for which we obtained spectroscopy for the first time, finding a redshift of 2.77 for the quasar. There is evidence for microlensing in at least one image. We perform detailed mass modeling of this system using archival imaging data, and find that the unusually large shear responsible for the diamond-like configuration can be attributed mainly to a faint companion $\sim 4''$ away, and to a galaxy group/cluster $\sim 30''$ away. We also set limits of $z\sim0.5-1.5$ on the redshift of the lensing galaxy, based on its brightness, the image separation of the lensed images, and an analysis of the observed photometric flux ratios., Accepted by MNRAS. Based on feedback from the referee, the final sample of candidates is significantly different from the one in the previous submission. 15 pages + 4 pages of appendix, 7 figures, 3 tables