Showing 1–50 of 82 results for author: Agnello, A

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

Emulation Techniques for Scenario and Classical Control Design of Tokamak Plasmas

Authors: A. Agnello, N. C. Amorisco, A. Keats, G. K. Holt, J. Buchanan, S. Pamela, C. Vincent, G. McArdle

Abstract: The optimisation of scenarios and design of real-time-control in tokamaks, especially for machines still in design phase, requires a comprehensive exploration of solutions to the Grad-Shafranov (GS) equation over a high-dimensional space of plasma and coil parameters. Emulators can bypass the numerical issues in the GS equation, if a large enough library of equilibria is available. We train an ens… ▽ More The optimisation of scenarios and design of real-time-control in tokamaks, especially for machines still in design phase, requires a comprehensive exploration of solutions to the Grad-Shafranov (GS) equation over a high-dimensional space of plasma and coil parameters. Emulators can bypass the numerical issues in the GS equation, if a large enough library of equilibria is available. We train an ensemble of neural networks to emulate the typical shape-control targets (separatrix at midplane, X-points, divertor strike point, flux expansion, poloidal beta) as a function of plasma parameters and active coil currents for the range of plasma configurations relevant to spherical tokamaks with a super-X divertor, with percent-level accuracy. This allows a quick calculation of the classical-control shape matrices, potentially allowing real-time calculation at any point in a shot with sub-ms latency. We devise a hyperparameter sampler to select the optimal network architectures and quantify uncertainties on the model predictions. To generate the relevant training set, we devise a Markov-Chain Monte Carlo algorithm to produce large libraries of forward Grad-Shafranov solutions without the need for user intervention. The algorithm promotes equilibria with desirable properties, while avoiding parameter combinations resulting in problematic profiles or numerical issues in the integration of the GS equation. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Comments: Physics of Plasmas in print, ICDDPS-4 special issue, 15 pages, 6 figures

Journal ref: Physics of Plasmas 2024 Vol.31, Issue 4

VarIabiLity seLection of AstrophysIcal sources iN PTF (VILLAIN) II. Supervised classification of variable sources

Authors: S. H. Bruun, J. Hjorth, A. Agnello

Abstract: Context. Large, high-dimensional astronomical surveys require efficient data analysis. Automatic fitting of lightcurve variability and machine learning may assist in identification of sources including candidate quasars. Aims. We aim to classify sources from the Palomar Transient Factory (PTF) as quasars, stars or galaxies, and to examine model performance using variability and colours. We deter… ▽ More Context. Large, high-dimensional astronomical surveys require efficient data analysis. Automatic fitting of lightcurve variability and machine learning may assist in identification of sources including candidate quasars. Aims. We aim to classify sources from the Palomar Transient Factory (PTF) as quasars, stars or galaxies, and to examine model performance using variability and colours. We determine the added value of variability information as well as quantifying the performance when colours are not available. Methods. We use supervised learning in the form of a histogram-based gradient boosting classifier to predict spectroscopic SDSS classes using photometry. For comparison, we create models with structure function variability parameters only, magnitudes only and using all parameters. Results. We achieve highly accurate predictions for 71 million sources with lightcurves in PTF. The full model correctly identifies 92.49 % of spectroscopically confirmed quasars from the SDSS with a purity of 95.64 %. With only variability, the completeness is 34.97 % and the purity is 58.71 % for quasars. The predictions and probabilities of PTF objects belonging to each class are made available in a catalogue, VILLAIN-Cat, including magnitudes and variability parameters. Conclusions. We have developed a method for automatic and effective classification of PTF sources using magnitudes and variability. For similar supervised models, we recommend using at least 100,000 labeled objects, and we show how performance scales with data volume. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: 10 pages, 5 figures

VarIabiLity seLection of AstrophysIcal sources iN PTF (VILLAIN) I. Structure function fits to 71 million objects

Authors: S. H. Bruun, A. Agnello, J. Hjorth

Abstract: Context. Lightcurve variability is well-suited for characterising objects in surveys with high cadence and long baseline. This is especially relevant in view of the large datasets to be produced by the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). Aims. We aim to determine variability parameters for objects in the Palomar Transient Factory (PTF) and explore differences betwee… ▽ More Context. Lightcurve variability is well-suited for characterising objects in surveys with high cadence and long baseline. This is especially relevant in view of the large datasets to be produced by the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). Aims. We aim to determine variability parameters for objects in the Palomar Transient Factory (PTF) and explore differences between quasars (QSOs), stars and galaxies. We will relate variability and colour information in preparation for future surveys. Methods. We fit joint likelihoods to structure functions (SFs) of 71 million PTF lightcurves with a Markov Chain Monte Carlo method. For each object, we assume a power law SF and extract two parameters: the amplitude on timescales of one year, $A$, and a power law index, $γ$. With these parameters and colours in the optical (Pan-STARRS1) and mid infrared (WISE), we identify regions of parameter space dominated by different types of spectroscopically confirmed objects from SDSS. Candidate QSOs, stars and galaxies are selected to show their parameter distributions. Results. QSOs have high amplitude variations in the $R$ band, and the strongest timescale dependence of variability. Galaxies have a broader range of amplitudes and low timescale dependency. With variability and colours, we achieve a photometric selection purity of 99.3 % for QSOs. Even though hard cuts in monochromatic variability alone are not as effective as seven-band magnitude cuts, variability is useful in characterising object sub-classes. Through variability, we also find QSOs that were erroneously classified as stars in the SDSS. We discuss perspectives and computational solutions in view of the upcoming LSST. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: Accepted by A&A on 11/04/2023, 16 pages, 14 figures

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

Hierarchical Inference of the Lensing Convergence from Photometric Catalogs with Bayesian Graph Neural Networks

Authors: Ji Won Park, Simon Birrer, Madison Ueland, Miles Cranmer, Adriano Agnello, Sebastian Wagner-Carena, Philip J. Marshall, Aaron Roodman, the LSST Dark Energy Science Collaboration

Abstract: We present a Bayesian graph neural network (BGNN) that can estimate the weak lensing convergence ($κ$) from photometric measurements of galaxies along a given line of sight. The method is of particular interest in strong gravitational time delay cosmography (TDC), where characterizing the "external convergence" ($κ_{\rm ext}$) from the lens environment and line of sight is necessary for precise in… ▽ More We present a Bayesian graph neural network (BGNN) that can estimate the weak lensing convergence ($κ$) from photometric measurements of galaxies along a given line of sight. The method is of particular interest in strong gravitational time delay cosmography (TDC), where characterizing the "external convergence" ($κ_{\rm ext}$) from the lens environment and line of sight is necessary for precise inference of the Hubble constant ($H_0$). Starting from a large-scale simulation with a $κ$ resolution of $\sim$1$'$, we introduce fluctuations on galaxy-galaxy lensing scales of $\sim$1$''$ and extract random sightlines to train our BGNN. We then evaluate the model on test sets with varying degrees of overlap with the training distribution. For each test set of 1,000 sightlines, the BGNN infers the individual $κ$ posteriors, which we combine in a hierarchical Bayesian model to yield constraints on the hyperparameters governing the population. For a test field well sampled by the training set, the BGNN recovers the population mean of $κ$ precisely and without bias, resulting in a contribution to the $H_0$ error budget well under 1\%. In the tails of the training set with sparse samples, the BGNN, which can ingest all available information about each sightline, extracts more $κ$ signal compared to a simplified version of the traditional method based on matching galaxy number counts, which is limited by sample variance. Our hierarchical inference pipeline using BGNNs promises to improve the $κ_{\rm ext}$ characterization for precision TDC. The implementation of our pipeline is available as a public Python package, Node to Joy. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: 15 pages, 8 figures (+ 6 pages, 2 figures in Appendix). Submitted to ApJ. Code at https://github.com/jiwoncpark/node-to-joy

LensWatch: I. Resolved HST Observations and Constraints on the Strongly-Lensed Type Ia Supernova 2022qmx ("SN Zwicky")

Authors: J. D. R. Pierel, N. Arendse, S. Ertl, X. Huang, L. A. Moustakas, S. Schuldt, A. J. Shajib, Y. Shu, S. Birrer, M. Bronikowski, J. Hjorth, S. H. Suyu, S. Agarwal, A. Agnello, A. S. Bolton, S. Chakrabarti, C. Cold, F. Courbin, J. M. Della Costa, S. Dhawan, M. Engesser, O. D. Fox, C. Gall, S. Gomez, A. Goobar , et al. (17 additional authors not shown)

Abstract: Supernovae (SNe) that have been multiply-imaged by gravitational lensing are rare and powerful probes for cosmology. Each detection is an opportunity to develop the critical tools and methodologies needed as the sample of lensed SNe increases by orders of magnitude with the upcoming Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope. The latest such discovery is of the quadruply-image… ▽ More Supernovae (SNe) that have been multiply-imaged by gravitational lensing are rare and powerful probes for cosmology. Each detection is an opportunity to develop the critical tools and methodologies needed as the sample of lensed SNe increases by orders of magnitude with the upcoming Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope. The latest such discovery is of the quadruply-imaged Type Ia SN 2022qmx (aka, "SN Zwicky"; Goobar et al. 2022) at z = 0.3544. SN Zwicky was discovered by the Zwicky Transient Facility (ZTF) in spatially unresolved data. Here we present follow-up Hubble Space Telescope observations of SN Zwicky, the first from the multi-cycle "LensWatch" program (www.lenswatch.org). We measure photometry for each of the four images of SN Zwicky, which are resolved in three WFC3/UVIS filters (F475W, F625W, F814W) but unresolved with WFC3/IR F160W, and produce an analysis of the lensing system using a variety of independent lens modeling methods. We find consistency between time delays estimated with the single epoch of HST photometry and the lens model predictions constrained through the multiple image positions, with both inferring time delays of <1 day. Our lens models converge to an Einstein radius of (0.168+0.009-0.005)", the smallest yet seen in a lensed SN. The "standard candle" nature of SN Zwicky provides magnification estimates independent of the lens modeling that are brighter by ~1.5 mag and ~0.8 mag for two of the four images, suggesting significant microlensing and/or additional substructure beyond the flexibility of our image-position mass models. △ Less

Submitted 8 November, 2022; v1 submitted 7 November, 2022; originally announced November 2022.

Photometric redshift estimation of strongly lensed galaxies

Authors: Danial Langeroodi, Alessandro Sonnenfeld, Henk Hoekstra, Adriano Agnello

Abstract: Around $10^5$ strongly lensed galaxies are expected to be discovered with Euclid and the LSST. Utilising these large samples to study the inner structure of lens galaxies requires source redshifts, to turn lens models into mass measurements. However, obtaining spectroscopic source redshifts for large lens samples is prohibitive with the capacity of spectroscopic facilities. Alternatively, we study… ▽ More Around $10^5$ strongly lensed galaxies are expected to be discovered with Euclid and the LSST. Utilising these large samples to study the inner structure of lens galaxies requires source redshifts, to turn lens models into mass measurements. However, obtaining spectroscopic source redshifts for large lens samples is prohibitive with the capacity of spectroscopic facilities. Alternatively, we study the possibility of obtaining source photometric redshifts (photo-zs) for large lens samples. Our strategy consists of deblending the source and lens light by simultaneously modelling the lens and background source in all available photometric bands, and feeding the derived source colours to a template-fitting photo-z algorithm. We describe the lens and source light with a Sersic profile, and the lens mass with a Singular Isothermal Ellipsoid. We test our approach on a simulated and a real sample of lenses, both in broad-band photometry of the Hyper Suprime-Cam survey. We identify the deviations of the lens light from a Sersic profile and the contrast between the lens and source image as the main drivers of the source colour measurement error. We split the real sample based on the ratio $Λ$ of the lens to source surface brightness measured at the image locations. In the $Λ<1$ regime, the photo-z outlier fraction is $20\%$, and the accuracy of photo-z estimation is limited by the performance of template-fitting process. In the opposite regime, the photo-z outlier fraction is $75\%$, and the errors from the source colour measurements dominate the photo-z uncertainty. Measuring source photo-zs for lenses with $Λ<1$ poses no particular challenges, compared to isolated galaxies. For systems with significant lens light contamination, however, improving the description of the surface brightness distribution of the lens is required: a single Sersic model is not sufficiently accurate. △ Less

Submitted 20 December, 2022; v1 submitted 1 July, 2022; originally announced July 2022.

Comments: 15 pages, 13 figures, Accepted for Publication in A&A

Journal ref: A&A 669, A154 (2023)

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

Hubble constant and nuclear equation of state from kilonova spectro-photometric light curves

Authors: M. A. Pérez-García, L. Izzo, D. Barba, M. Bulla, A. Sagués-Carracedo, E. Pérez, C. Albertus, S. Dhawan, F. Prada, A. Agnello, C. R. Angus, S. H. Bruun, C. del Burgo, C. Dominguez-Tagle, C. Gall, A. Goobar, J. Hjorth, D. Jones, A. R. López-Sánchez, J. Sollerman

Abstract: The merger of two compact objects of which at least one is a neutron star is signalled by transient electromagnetic emission in a kilonova (KN). This event is accompanied by gravitational waves and possibly other radiation messengers such as neutrinos or cosmic rays. The electromagnetic emission arises from the radioactive decay of heavy $r-$process elements synthesized in the material ejected dur… ▽ More The merger of two compact objects of which at least one is a neutron star is signalled by transient electromagnetic emission in a kilonova (KN). This event is accompanied by gravitational waves and possibly other radiation messengers such as neutrinos or cosmic rays. The electromagnetic emission arises from the radioactive decay of heavy $r-$process elements synthesized in the material ejected during and after the merger. In this paper we show that the analysis of KNe light curves can provide cosmological distance measurements and constrain the properties of the ejecta. In this respect, MAAT, the new Integral Field Unit in the OSIRIS spectrograph on the $10.4$ m Gran Telescopio CANARIAS (GTC), is well suited for the study of KNe by performing absolute spectro-photometry over the entire 3600-10000 Angstron spectral range. Here, we study the most representative cases regarding the scientific interest of KNe from binary neutron stars, and we evaluate the observational prospects and performance of MAAT on the GTC to do the following: a) study the impact of the equation of state on the KN light curve, and determine to what extent bounds on neutron star (NS) radii or compactness deriving from KN peak magnitudes can be identified and b) measure the Hubble constant, $H_0$, with precision improved by up to 40$\%$, when both gravitational wave data and photometric-light curves are used. In this context we discuss how the equation of state, the viewing angle, and the distance affect the precision and estimated value of $H_0$. △ Less

Submitted 14 December, 2022; v1 submitted 31 March, 2022; originally announced April 2022.

Comments: 21 pages, 23 figures, matches published version in A&A 666, A67 (2022)

Journal ref: A&A 666, A67 (2022)

Cosmology Intertwined: A Review of the Particle Physics, Astrophysics, and Cosmology Associated with the Cosmological Tensions and Anomalies

Authors: Elcio Abdalla, Guillermo Franco Abellán, Amin Aboubrahim, Adriano Agnello, Ozgur Akarsu, Yashar Akrami, George Alestas, Daniel Aloni, Luca Amendola, Luis A. Anchordoqui, Richard I. Anderson, Nikki Arendse, Marika Asgari, Mario Ballardini, Vernon Barger, Spyros Basilakos, Ronaldo C. Batista, Elia S. Battistelli, Richard Battye, Micol Benetti, David Benisty, Asher Berlin, Paolo de Bernardis, Emanuele Berti, Bohdan Bidenko , et al. (178 additional authors not shown)

Abstract: In this paper we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant $H_0$, the $σ_8$--$S_8$ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of system… ▽ More In this paper we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant $H_0$, the $σ_8$--$S_8$ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the $5.0\,σ$ tension between the {\it Planck} CMB estimate of the Hubble constant $H_0$ and the SH0ES collaboration measurements. After showing the $H_0$ evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the {\it Planck} CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density $Ω_m$, and the amplitude or rate of the growth of structure ($σ_8,fσ_8$). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the $H_0$--$S_8$ tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals.[Abridged] △ Less

Submitted 24 April, 2022; v1 submitted 11 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021. 224 pages, 27 figures. Accepted for publication in JHEAp

Journal ref: J. High En. Astrophys. 2204, 002 (2022)

TDCOSMO. IX. Systematic comparison between lens modelling software programs: time delay prediction for WGD 2038-4008

Authors: A. J. Shajib, K. C. Wong, S. Birrer, S. H. Suyu, T. Treu, E. Buckley-Geer, H. Lin, C. E. Rusu, J. Poh, A. Palmese, A. Agnello, M. W. Auger-Williams, A. Galan, S. Schuldt, D. Sluse, F. Courbin, J. Frieman, M. Millon

Abstract: The importance of alternative methods for measuring the Hubble constant, such as time-delay cosmography, is highlighted by the recent Hubble tension. It is paramount to thoroughly investigate and rule out systematic biases in all measurement methods before we can accept new physics as the source of this tension. In this study, we perform a check for systematic biases in the lens modelling procedur… ▽ More The importance of alternative methods for measuring the Hubble constant, such as time-delay cosmography, is highlighted by the recent Hubble tension. It is paramount to thoroughly investigate and rule out systematic biases in all measurement methods before we can accept new physics as the source of this tension. In this study, we perform a check for systematic biases in the lens modelling procedure of time-delay cosmography by comparing independent and blind time-delay predictions of the system WGD 2038-4008 from two teams using two different software programs: GLEE and LENSTRONOMY. The predicted time delays from the two teams incorporate the stellar kinematics of the deflector and the external convergence from line-of-sight structures. The un-blinded time-delay predictions from the two teams agree within $1.2σ$, implying that once the time delay is measured the inferred Hubble constant will also be mutually consistent. However, there is a $\sim$4$σ$ discrepancy between the power-law model slope and external shear, which is a significant discrepancy at the level of lens models before the stellar kinematics and the external convergence are incorporated. We identify the difference in the reconstructed point spread function (PSF) to be the source of this discrepancy. When the same reconstructed PSF was used by both teams, we achieved excellent agreement, within $\sim$0.6$σ$, indicating that potential systematics stemming from source reconstruction algorithms and investigator choices are well under control. We recommend that future studies supersample the PSF as needed and marginalize over multiple algorithms or realizations for the PSF reconstruction to mitigate the systematics associated with the PSF. A future study will measure the time delays of the system WGD 2038-4008 and infer the Hubble constant based on our mass models. △ Less

Submitted 21 September, 2022; v1 submitted 22 February, 2022; originally announced February 2022.

Comments: Accepted version by A&A. 34 pages, 26 figures, 3 tables

Journal ref: A&A 667, A123 (2022)

The halo of M105 and its group environment as traced by planetary nebula populations: II. Using kinematics of single stars to unveil the presence of intragroup light around the Leo I galaxies NGC 3384 and M105

Authors: J. Hartke, M. Arnaboldi, O. Gerhard, L. Coccato, M. Merrifield, K. Kuijken, C. Pulsoni, A. Agnello, S. Bhattacharya, C. Spiniello, A. Cortesi, K. C. Freeman, N. R. Napolitano, A. J. Romanowsky

Abstract: M105 is an early-type galaxy in the nearby Leo I group, the closest galaxy group to contain all galaxy types and therefore an excellent environment to explore the low-mass end of intra-group light (IGL) assembly. We present a new extended kinematic survey of planetary nebulae (PNe) in M105 and the surrounding 30'x30' in the Leo I group with the Planetary Nebula Spectrograph. We use PNe as kinemati… ▽ More M105 is an early-type galaxy in the nearby Leo I group, the closest galaxy group to contain all galaxy types and therefore an excellent environment to explore the low-mass end of intra-group light (IGL) assembly. We present a new extended kinematic survey of planetary nebulae (PNe) in M105 and the surrounding 30'x30' in the Leo I group with the Planetary Nebula Spectrograph. We use PNe as kinematic tracers of the diffuse stellar light in the halo and IGL and employ Gaussian mixture models to separate contributions from the companion galaxy NGC 3384 and associate PNe with halo and IGL components around M105. We present a catalogue of 314 PNe and firmly associate 93 with NGC 3384 and 169 with M105. The PNe in M105 are further associated with its halo and the surrounding exponential envelope. We construct smooth velocity and velocity dispersion fields and calculate projected rotation, velocity dispersion, and $λ_R$ profiles for each component. Halo PNe exhibit declining velocity dispersion and rotation profiles, while the velocity dispersion and rotation of the exponential envelope increase notably at large radii. We identify three regimes with distinct kinematics that are linked to distinct stellar population properties: (i) the rotating core (within $1~R_\mathrm{eff}$) formed in situ and dominated by metal-rich ([M/H]~0) stars likely formed in situ, (ii) the halo from 1 to $7.5~R_\mathrm{eff}$ consisting of intermediate-metallicity stars ([M/H]>-1), either formed in situ or brought in through major mergers, and (iii) the exponential envelope reaching beyond our farthest data point at 16 $R_\mathrm{eff}$, predominately composed of metal-poor ([M/H]<-1) stars. The high velocity dispersion and moderate rotation of the latter are consistent with that measured for dwarf satellite galaxies in the Leo I group, indicating that the exponential envelope traces the transition to the IGL. △ Less

Submitted 21 January, 2022; originally announced January 2022.

Comments: Submitted to A&A, comments welcome. 20 pages, 14 figures, abstract abridged

Journal ref: A&A 663, A12 (2022)

Detection of $^7$Be II in the Small Magellanic Cloud

Authors: Luca Izzo, Paolo Molaro, Gabriele Cescutti, Elias Aydi, Pierluigi Selvelli, Eamonn Harvey, Adriano Agnello, Piercarlo Bonifacio, Massimo Della Valle, Ernesto Guido, Margarita Hernanz

Abstract: We analyse high resolution spectra of two classical novae that exploded in the Small Magellanic Cloud. $^7$Be II resonance transitions are detected in both ASASSN-19qv and ASASSN-20ni novae. This is the first detection outside the Galaxy and confirms that thermo-nuclear runaway reactions, leading to the $^7$Be formation, are effective also in the low metallicity regime, characteristic of the SMC.… ▽ More We analyse high resolution spectra of two classical novae that exploded in the Small Magellanic Cloud. $^7$Be II resonance transitions are detected in both ASASSN-19qv and ASASSN-20ni novae. This is the first detection outside the Galaxy and confirms that thermo-nuclear runaway reactions, leading to the $^7$Be formation, are effective also in the low metallicity regime, characteristic of the SMC. Derived yields are of N($^7$Be=$^7$Li)/N(H) = (5.3 $\pm$ 0.2) $\times$ 10$^{-6}$ which are a factor 4 lower than the typical values of the Galaxy. Inspection of two historical novae in the Large Magellanic Cloud observed with IUE in 1991 and 1992 showed also the possible presence of $^7$Be and similar yields. For an ejecta of $M_{H,ej} =$ 10$^{-5}$ M$_{\odot}$, the amount of $^7$Li produced is of $M_{^7 Li} = (3.7 \pm 0.6) \times 10^{-10}$ M$_{\odot}$ per nova event. Detailed chemical evolutionary model for the SMC shows that novae could have made an amount of lithium in the SMC corresponding to a fractional abundance of A(Li) $\approx$ 2.6. Therefore, it is argued that a comparison with the abundance of Li in the SMC, as measured by its interstellar medium, could effectively constrain the amount of the initial abundance of primordial Li, which is currently controversial. △ Less

Submitted 22 December, 2021; originally announced December 2021.

Comments: 13 pages, 17 figures, 4 tables. Accepted for publication in MNRAS

J1721$+$8842: a gravitationally lensed binary quasar with a proximate damped Lyman-$α$ absorber

Authors: C. Lemon, M. Millon, D. Sluse, F. Courbin, M. Auger, J. Chan, E. Paic, A. Agnello

Abstract: High-redshift binary quasars provide key insights into mergers and quasar activity, and are useful tools for probing the spatial kinematics and chemistry of galaxies along the line-of-sight. However, only three sub-10-kpc binaries have been confirmed above $z=1$. Gravitational lensing would provide a way to easily resolve such binaries, study them in higher resolution, and provide more sightlines,… ▽ More High-redshift binary quasars provide key insights into mergers and quasar activity, and are useful tools for probing the spatial kinematics and chemistry of galaxies along the line-of-sight. However, only three sub-10-kpc binaries have been confirmed above $z=1$. Gravitational lensing would provide a way to easily resolve such binaries, study them in higher resolution, and provide more sightlines, though the required alignment with a massive foreground galaxy is rare. Through image deconvolution of StanCam Nordic Optical Telescope (NOT) monitoring data, we reveal two further point sources in the known, $z \approx 2.38$, quadruply lensed quasar (quad), J1721+8842. An ALFOSC/NOT long-slit spectrum shows that the brighter of these two sources is a quasar with $z = 2.369 \pm 0.007$ based on the C III] line, while the C III] redshift of the quad is $z = 2.364 \pm 0.003$. Lens modelling using point source positions rules out a single source model, favouring an isothermal lens mass profile with two quasar sources separated by $\sim6.0$ kpc (0.73$^{\prime \prime}$) in projection. Given the resolving ability from lensing and current lensed quasar statistics, this discovery suggests a large population of undiscovered, unlensed sub-10-kpc binaries. We also analyse spectra of two images of the quad, showing narrow Ly$α$ emission within the trough of a proximate damped Ly$α$ absorber (PDLA). An apparent mismatch between the continuum and narrow line flux ratios provides a new potential tool for simultaneously studying microlensing and the quasar host galaxy. Signs of the PDLA are also seen in the second source, however a deeper spectrum is still required to confirm this. Thanks to the multiple lines-of-sight from lensing and two quasar sources, this system offers simultaneous sub-parsec and kpc-scale probes of a PDLA. △ Less

Submitted 2 September, 2021; originally announced September 2021.

Comments: 10 pages, 7 figures, submitted to A&A

Journal ref: A&A 657, A113 (2022)

VEXAS: VISTA EXtension to Auxiliary Surveys -- Data Release 2: Machine-learning based classification of sources in the Southern Hemisphere

Authors: V. Khramtsov, C. Spiniello, A. Agnello, A. Sergeyev

Abstract: We present the second public data release (DR) of the VISTA EXtension to Auxiliary Surveys (VEXAS), where we classify objects into stars, galaxies and quasars based on an ensemble of machine learning algorithms. The aim of VEXAS is to build the widest multi-wavelength catalogue, providing reference magnitudes, colours and morphological information for a large number of scientific uses. We apply an… ▽ More We present the second public data release (DR) of the VISTA EXtension to Auxiliary Surveys (VEXAS), where we classify objects into stars, galaxies and quasars based on an ensemble of machine learning algorithms. The aim of VEXAS is to build the widest multi-wavelength catalogue, providing reference magnitudes, colours and morphological information for a large number of scientific uses. We apply an ensemble of 32 different machine learning models, based on three different algorithms and on different magnitude sets, training samples and classification problems on the three VEXAS DR1 optical+infrared (IR) tables. The tables were created in DR1 cross-matching VISTA near-IR data with WISE far-IR data and with optical magnitudes from the Dark Energy Survey (VEXAS-DESW), the Sky Mapper Survey (VEXAS-SMW), and the PanSTARRS (VEXAS-PSW). We assemble a large table of spectroscopically confirmed objects (415 628 unique objects), based on the combination of 6 different spectroscopic surveys that we use for training. We develop feature imputation to classify also objects for which magnitudes in one or more bands are missing. We classify in total ~90 million objects in the Southern Hemisphere. Among these,~62.9M (~52.6M) are classified as 'high confidence' ('secure') stars, ~920k (~750k) as 'high confidence' ('secure') quasars and ~34.8M (~34.1M) as 'high confidence' ('secure') galaxies, with probabilities $p_{\rm class}\ge 0.7$ ($p_{\rm class}\ge 0.9$). The density of high-confidence extragalactic objects varies strongly with the survey depth: at $p_{\rm class}\ge 0.7$, there are 111/deg$^2$ quasars in the VEXAS-DESW footprint and 103/deg$^2$ in the VEXAS-PSW footprint, while only 10.7/deg$^2$ in the VEXAS-SM footprint. Improved depth in the midIR and coverage in the optical and nearIR are needed for the SM footprint that is not already covered by DESW and PSW. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Comments: 25 pages, 18 figures, 8 tables. Accepted for publication on A&A. The VEXAS tables are publicly available through the ESO Phase 3 here: https://archive.eso.org/scienceportal/home?data_collection=VEXAS. The DR2 tables update the DR1 with the addition of imputed magnitudes and membership probabilities to each of the three classes

Journal ref: A&A 651, A69 (2021)

Dwarf galaxies in the MATLAS survey: Hubble Space Telescope observations of the globular cluster system in the ultra-diffuse galaxy MATLAS-2019

Authors: Oliver Müller, Patrick R. Durrell, Francine R. Marleau, Pierre-Alain Duc, Sungsoon Lim, Lorenzo Posti, Adriano Agnello, Rúben Sánchez-Janssen, Mélina Poulain, Rebecca Habas, Eric Emsellem, Sanjaya Paudel, Remco F. J. van der Burg, Jérémy Fensch

Abstract: Ultra-diffuse galaxies (UDGs) are very low-surface brightness galaxies with large effective radii. Spectroscopic measurements of a few UDGs have revealed a low dark matter content, based on the internal motion of stars or globular clusters (GCs). This is in contrast to the large number of GCs found for these systems, from which it would be expected to correspond to a large dark matter halo mass. H… ▽ More Ultra-diffuse galaxies (UDGs) are very low-surface brightness galaxies with large effective radii. Spectroscopic measurements of a few UDGs have revealed a low dark matter content, based on the internal motion of stars or globular clusters (GCs). This is in contrast to the large number of GCs found for these systems, from which it would be expected to correspond to a large dark matter halo mass. Here we present HST+ACS observations for the UDG MATLAS-2019 in the NGC5846 group. Using the F606W and F814W filters, we trace the GC population two magnitudes below the peak of the GC luminosity function (GCLF). Employing Bayesian considerations, we identify 26+-6 GCs associated with the dwarf, yielding a large specific frequency of S_N=58+-14. We use the turnover of the GCLF to derive a distance of 21+-2 Mpc, which is consistent with the NGC5846 group of galaxies. Due to the superior image quality of the HST, we are able to resolve the GCs and measure their sizes, which are consistent with the sizes of GCs around Local Group galaxies. Using the linear relation between the total mass of galaxies and of GCs, we derive a halo mass of 0.9(+-0.2) *10^11 M_solar (M_solar/L_solar>1000). The high abundance of GCs, together with the small uncertainties, make MATLAS-2019 one of the most extreme UDGs, which likely sets an upper limit of the number of GCs for UDGs. △ Less

Submitted 23 September, 2021; v1 submitted 26 January, 2021; originally announced January 2021.

Comments: 16 pages, 10 figures, 1 table, accepted for publication in APJ

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

The Young Supernova Experiment: Survey Goals, Overview, and Operations

Authors: D. O. Jones, R. J. Foley, G. Narayan, J. Hjorth, M. E. Huber, P. D. Aleo, K. D. Alexander, C. R. Angus, K. Auchettl, V. F. Baldassare, S. H. Bruun, K. C. Chambers, D. Chatterjee, D. L. Coppejans, D. A. Coulter, L. DeMarchi, G. Dimitriadis, M. R. Drout, A. Engel, K. D. French, A. Gagliano, C. Gall, T. Hung, L. Izzo, W. V. Jacobson-Galán , et al. (46 additional authors not shown)

Abstract: Time domain science has undergone a revolution over the past decade, with tens of thousands of new supernovae (SNe) discovered each year. However, several observational domains, including SNe within days or hours of explosion and faint, red transients, are just beginning to be explored. Here, we present the Young Supernova Experiment (YSE), a novel optical time-domain survey on the Pan-STARRS tele… ▽ More Time domain science has undergone a revolution over the past decade, with tens of thousands of new supernovae (SNe) discovered each year. However, several observational domains, including SNe within days or hours of explosion and faint, red transients, are just beginning to be explored. Here, we present the Young Supernova Experiment (YSE), a novel optical time-domain survey on the Pan-STARRS telescopes. Our survey is designed to obtain well-sampled $griz$ light curves for thousands of transient events up to $z \approx 0.2$. This large sample of transients with 4-band light curves will lay the foundation for the Vera C. Rubin Observatory and the Nancy Grace Roman Space Telescope, providing a critical training set in similar filters and a well-calibrated low-redshift anchor of cosmologically useful SNe Ia to benefit dark energy science. As the name suggests, YSE complements and extends other ongoing time-domain surveys by discovering fast-rising SNe within a few hours to days of explosion. YSE is the only current four-band time-domain survey and is able to discover transients as faint $\sim$21.5 mag in $gri$ and $\sim$20.5 mag in $z$, depths that allow us to probe the earliest epochs of stellar explosions. YSE is currently observing approximately 750 square degrees of sky every three days and we plan to increase the area to 1500 square degrees in the near future. When operating at full capacity, survey simulations show that YSE will find $\sim$5000 new SNe per year and at least two SNe within three days of explosion per month. To date, YSE has discovered or observed 8.3% of the transient candidates reported to the International Astronomical Union in 2020. We present an overview of YSE, including science goals, survey characteristics and a summary of our transient discoveries to date. △ Less

Submitted 5 January, 2021; v1 submitted 19 October, 2020; originally announced October 2020.

Comments: ApJ, in press; more information at https://yse.ucsc.edu/

Mixture Models for Photometric Redshifts

Authors: Z. Ansari, A. Agnello, C. Gall

Abstract: Determining photometric redshifts to high accuracy is paramount to measure distances in wide-field cosmological experiments. With only photometric information at hand, photo-zs are prone to systematic uncertainties in the intervening extinction and the unknown underlying spectral-energy distribution of different astrophysical sources. Here, we aim to resolve these model degeneracies and obtain a c… ▽ More Determining photometric redshifts to high accuracy is paramount to measure distances in wide-field cosmological experiments. With only photometric information at hand, photo-zs are prone to systematic uncertainties in the intervening extinction and the unknown underlying spectral-energy distribution of different astrophysical sources. Here, we aim to resolve these model degeneracies and obtain a clear separation between intrinsic physical properties of astrophysical sources and extrinsic systematics. We aim at estimates of the full photo-z probability distributions, and their uncertainties. We perform a probabilistic photo-z determination using Mixture Density Networks (MDN). The training data-set is composed of optical ($griz$) point-spread-function and model magnitudes and extinction measurements from the SDSS-DR15, and WISE midinfrared ($3.4 μ$m and $4.6 μ$m) model magnitudes. We use Infinite Gaussian Mixture models to classify the objects in our data-set as stars, galaxies or quasars, and to determine the number of MDN components to achieve optimal performance. The fraction of objects that are correctly split into the main classes is 94%. Our method improves the bias of photometric redshift estimation (i.e. the mean $Δz$ = (zp - zs)/(1 + zs)) by one order of magnitude compared to the SDSS photo-z, and decreases the fraction of $3 σ$ outliers (i.e. 3rms$(Δz) < Δz$). The relative, root-mean-square systematic uncertainty in our resulting photo-zs is down to 1.7% for low-redshift galaxies (zs $<$ 0.5). We have demonstrated the feasibility of machine-learning based methods that produce full probability distributions for photo-z estimates with a performance that is competitive with state-of-the art techniques. Our method can be applied to wide-field surveys where extinction can vary significantly across the sky and with sparse spectroscopic calibration samples. △ Less

Submitted 14 October, 2020; originally announced October 2020.

Comments: 14 pages, 9 figures, 7 tables, submitted to A&A 14/10/2020

Journal ref: A&A 650, A90 (2021)

Cluster Strong Lensing with Hierarchical Inference

Authors: Pietro Bergamini, Adriano Agnello, Gabriel Bartosch Caminha

Abstract: Lensing by galaxy clusters is a versatile probe of cosmology and extragalactic astrophysics, but the accuracy of some of its predictions is limited by the simplified models adopted to reduce the (otherwise untractable) number of degrees of freedom. We aim at cluster lensing models where the parameters of all cluster-member galaxies are free to vary around some common scaling relations with non-zer… ▽ More Lensing by galaxy clusters is a versatile probe of cosmology and extragalactic astrophysics, but the accuracy of some of its predictions is limited by the simplified models adopted to reduce the (otherwise untractable) number of degrees of freedom. We aim at cluster lensing models where the parameters of all cluster-member galaxies are free to vary around some common scaling relations with non-zero scatter, and deviate significantly from them if and only if the data require it. We have devised a Bayesian hierarchical inference framework, which enables the determination of all lensing parameters and of the scaling-relation hyperparameters, including intrinsic scatter, from lensing constraints and (if given) stellar kinematic measurements. We achieve this through BayesLens, a purpose-built wrapper around common parametric lensing codes for the lensing likelihood that can sample the posterior on parameters and hyperparameters, which we release with this paper. We have run functional tests of our code against simple mock cluster lensing datasets with realistic uncertainties. The parameters and hyperparameters are recovered within their 68% credibility ranges, and the positions of all the "observed" multiple images are accurately reproduced by the BayeLens best-fit model, without overfitting. We have shown that an accurate description of cluster member galaxies is attainable, despite the large number of degrees of freedom, through fast and tractable inference. This extends beyond the state-of-the-art of current cluster lensing models. The precise impact on studies of cosmography, galaxy evolution and high-redshift galaxy populations can then be quantified on real galaxy clusters. While other sources of systematics exist and may be significant in real clusters, our results show that the contribution of intrinsic scatter in cluster-member populations can now be controlled. △ Less

Submitted 2 December, 2020; v1 submitted 26 August, 2020; originally announced August 2020.

Comments: A&A in press (first subm.to A&A on 19/11/2019)

Journal ref: A&A 648, A123 (2021)

A new measurement of the Hubble constant using Type Ia supernovae calibrated with surface brightness fluctuations

Authors: N. Khetan, L. Izzo, M. Branchesi, R. Wojtak, M. Cantiello, C. Murugeshan, A. Agnello, E. Cappellaro, M. Della Valle, C. Gall, J. Hjorth, S. Benetti, E. Brocato, J. Burke, D. Hiramatsu, D. Andrew Howell, L. Tomasella, S. Valenti

Abstract: We present a new calibration of the peak absolute magnitude of SNe Type Ia based on the Surface Brightness Fluctuations (SBF) method, aimed at measuring the value of the Hubble constant. We build a sample of calibrating anchors consisting of 24 SNe hosted in galaxies having SBF distance measurements. Applying a hierarchical Bayesian approach, we calibrate the SNe luminosity and extend it into the… ▽ More We present a new calibration of the peak absolute magnitude of SNe Type Ia based on the Surface Brightness Fluctuations (SBF) method, aimed at measuring the value of the Hubble constant. We build a sample of calibrating anchors consisting of 24 SNe hosted in galaxies having SBF distance measurements. Applying a hierarchical Bayesian approach, we calibrate the SNe luminosity and extend it into the Hubble flow by using a sample of 96 SNe Ia in the redshift range $0.02 < z < 0.075$, extracted from the Combined Pantheon Sample. We estimate a value of $H_0 = 70.50 \pm 2.37(stat) \pm 3.38(sys)$ $\text{km}\ \text{s}^{-1}\ \text{Mpc}^{-1}$ (i.e. $3.4\% stat, 4.8\% sys$), which is in agreement with the value obtained using the tip of the red giant branch calibration, and consistent within the errors with the value obtained from SNe Type Ia calibrated with Cepheids and the one inferred from the analysis of the cosmic microwave background. We find that the SNe Ia distance moduli calibrated with SBF are on average larger by 0.07 mag than the ones calibrated with Cepheids. Our results point to possible differences among SNe in different types of galaxies, which could originate from different local environments and/or SNe Ia progenitor properties. Sampling different host galaxy type, SBF offers a complementary approach to Cepheids which is important in addressing possible systematics. As the SBF method has the ability to reach larger distances than Cepheids, the impending entry of LSST and JWST into operation will increase the number of SNe Ia hosted in galaxies where SBF distances can be measured, making SBF measurements attractive for improving the calibration of SNe Ia, and in the estimation of $H_0$. △ Less

Submitted 3 March, 2021; v1 submitted 18 August, 2020; originally announced August 2020.

Comments: 22 pages, 14 figures, 9 tables. Comments are very welcome

Journal ref: A&A 647, A72 (2021)

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)

White Paper on MAAT@GTC

Authors: Francisco Prada, Robert Content, Ariel Goobar, Luca Izzo, Enrique Pérez, Adriano Agnello, Carlos del Burgo, Vik Dhillon, José M. Diego, Lluis Galbany, Jorge García-Rojas, David Jones, Jon Lawrence, Eduardo Martín, Evencio Mediavilla, M. Ángeles Pérez García, Jorge Sánchez Almeida, José A. Acosta Pulido, Angel R. López-Sánchez, Santiago Arribas, Francisco J. Carrera, Amalia Corral, Inmaculada Domínguez, Silvia Mateos, Silvia Martínez Nuñez , et al. (19 additional authors not shown)

Abstract: MAAT is proposed as a visitor mirror-slicer optical system that will allow the OSIRIS spectrograph on the 10.4-m Gran telescopio CANARIAS (GTC) the capability to perform Integral Field Spectroscopy (IFS) over a seeing-limited FoV 14.20''x10'' with a slice width of 0.303''. MAAT@GTC will enhance the resolution power of OSIRIS by 1.6 times as compared to its 0.6'' wide long-slit. All the eleven OSIR… ▽ More MAAT is proposed as a visitor mirror-slicer optical system that will allow the OSIRIS spectrograph on the 10.4-m Gran telescopio CANARIAS (GTC) the capability to perform Integral Field Spectroscopy (IFS) over a seeing-limited FoV 14.20''x10'' with a slice width of 0.303''. MAAT@GTC will enhance the resolution power of OSIRIS by 1.6 times as compared to its 0.6'' wide long-slit. All the eleven OSIRIS grisms and volume-phase holographic gratings will be available to provide broad spectral coverage with moderate resolution (R=600 up to 4100) in the 3600 - 10000 Å wavelength range. MAAT unique observing capabilities will broaden its use to the needs of the GTC community to unveil the nature of most striking phenomena in the universe well beyond time-domain astronomy. The GTC equipped with OSIRIS+MAAT will also play a fundamental role in synergy with other facilities, some of them operating on the northern ORM at La Palma. This White Paper presents the different aspects of MAAT@GTC - including scientific and technical specifications, outstanding science cases, and an outline of the instrument concept. △ Less

Submitted 19 July, 2020; v1 submitted 3 July, 2020; originally announced July 2020.

Comments: 49 pages, 36 figures

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)

A spectroscopic study of MATLAS-2019 with MUSE: an ultra-diffuse galaxy with an excess of old globular clusters

Authors: Oliver Müller, Francine R. Marleau, Pierre-Alain Duc, Rebecca Habas, Jérémy Fensch, Eric Emsellem, Mélina Poulain, Sungsoon Lim, Adriano Agnello, Patrick Durrell, Sanjaya Paudel, Rubén Sánchez-Janssen, Remco F. J. van der Burg

Abstract: The MATLAS deep imaging survey has uncovered a plethora of dwarf galaxies in the low density environment it has mapped. A fraction of them are unusually extended and have a low-surface brightness. Among these so-called ultra-diffuse galaxies, a few seem to host an excess of globular clusters. With the integral-field unit spectrograph MUSE we have observed one of these galaxies - MATLAS J15052031+0… ▽ More The MATLAS deep imaging survey has uncovered a plethora of dwarf galaxies in the low density environment it has mapped. A fraction of them are unusually extended and have a low-surface brightness. Among these so-called ultra-diffuse galaxies, a few seem to host an excess of globular clusters. With the integral-field unit spectrograph MUSE we have observed one of these galaxies - MATLAS J15052031+0148447 (MATLAS-2019) - located towards the nearby group NGC 5846 and measured its systemic velocity,age, and metallicity, and that of its globular clusters candidates. For the stellar body of MATLAS-2019 we derive a metallicity of -1.33+0.19-0.01 dex and an age of 11.2+1.8-0.8 Gyr. For some of the individual GCs and the stacked GC population, we derive consistent ages and metallicities. From the 11 confirmed globular clusters and using a Markov Chain Monte Carlo approach we derived a dynamical mass-to-light ratio of 4.2+8.6-3.4M/L. This is at the lower end of the luminosity-mass scaling relation defined by the Local Group dwarf galaxies. Furthermore, we couldn't confirm nor reject the possibility of a rotational component of the GC system. If present, this would further modify the inferred mass. Follow-up observations of the globular cluster population and of the stellar body of the galaxy are needed to assess whether this galaxy is lacking dark matter like it was suggested for the pair of dwarf galaxies in the field of NGC 1052, or if this is a miss-interpretation arising from systematic uncertainties of the method commonly used for these systems and the large uncertainties of the individual globular cluster velocities. △ Less

Submitted 8 June, 2020; originally announced June 2020.

Comments: 10 pages, 9 figures, 2 tables, accepted for publication in A&A

Journal ref: A&A 640, A106 (2020)

Testing the Evolution of the Correlations between Supermassive Black Holes and their Host Galaxies using Eight Strongly Lensed Quasars

Authors: Xuheng Ding, Tommaso Treu, Simon Birrer, Adriano Agnello, Dominique Sluse, Chris Fassnacht, Matthew W. Auger, Kenneth C. Wong, Sherry H. Suyu, Takahiro Morishita, Cristian E. Rusu, Aymeric Galan

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

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

Comments: 12 pages, 4 figures, 3 tables. MNRAS in press. Comments welcome

Journal ref: MNRAS, (2021), 501, 269

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)

Observation and Confirmation of Nine Strong Lensing Systems in Dark Energy Survey Year 1 Data

Authors: B. Nord, E. Buckley-Geer, H. Lin, N. Kuropatkin, T. Collett, D. L. Tucker, H. T. Diehl, A. Agnello, A. Amara, T. M. C. Abbott, S. Allam, J. Annis, S. Avila, K. Bechtol, D. Brooks, D. L. Burke, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, C. E. Cunha, L. N. da Costa, C. Davis, J. De Vicente, P. Doel, T. F. Eifler , et al. (42 additional authors not shown)

Abstract: We describe the observation and confirmation of \nbconfirmtext\ new strong gravitational lenses discovered in Year 1 data from the Dark Energy Survey (DES). We created candidate lists based on a) galaxy group and cluster samples and b) photometrically selected galaxy samples. We selected 46 candidates through visual inspection and then used the Gemini Multi-Object Spectrograph (GMOS) at the Gemini… ▽ More We describe the observation and confirmation of \nbconfirmtext\ new strong gravitational lenses discovered in Year 1 data from the Dark Energy Survey (DES). We created candidate lists based on a) galaxy group and cluster samples and b) photometrically selected galaxy samples. We selected 46 candidates through visual inspection and then used the Gemini Multi-Object Spectrograph (GMOS) at the Gemini South telescope to acquire spectroscopic follow-up of 21 of these candidates. Through analysis of this spectroscopic follow-up data, we confirmed nine new lensing systems and rejected 2 candidates, but the analysis was inconclusive on 10 candidates. For each of the confirmed systems, we report measured spectroscopic properties, estimated \einsteinradiussub, and estimated enclosed masses. The sources that we targeted have an i-band surface brightness range of iSB ~ 22 - 24 mag arcsec^2 and a spectroscopic redshift range of zspec ~0.8 - 2.6. The lens galaxies have a photometric redshift range of zlens ~ 0.3 - 0.7. The lensing systems range in image-lens separation 2 - 9 arcsec and in enclosed mass 10^12 - 10^13 Msol. △ Less

Submitted 5 November, 2019; originally announced November 2019.

Report number: FERMILAB-PUB-17-042-PPD

Newly discovered dwarf galaxies in the MATLAS low density fields

Authors: Rebecca Habas, Francine R. Marleau, Pierre-Alain Duc, Patrick R. Durrell, Sanjaya Paudel, Mélina Poulain, Rubén Sánchez-Janssen, Sreevarsha Sreejith, Joanna Ramasawmy, Bryson Stemock, Christopher Leach, Jean-Charles Cuillandre, Stephen Gwyn, Adriano Agnello, Michal Bílek, Jérémy Fensch, Oliver Müller, Eric W. Peng, Remco F. J. van der Burg

Abstract: We present the photometric properties of 2210 newly identified dwarf galaxy candidates in the MATLAS fields. The Mass Assembly of early Type gaLAxies with their fine Structures (MATLAS) deep imaging survey mapped $\sim$142 deg$^2$ of the sky around nearby isolated early type galaxies using MegaCam on the Canada-France-Hawaii Telescope, reaching surface brightnesses of $\sim$ 28.5 - 29 in the g-ban… ▽ More We present the photometric properties of 2210 newly identified dwarf galaxy candidates in the MATLAS fields. The Mass Assembly of early Type gaLAxies with their fine Structures (MATLAS) deep imaging survey mapped $\sim$142 deg$^2$ of the sky around nearby isolated early type galaxies using MegaCam on the Canada-France-Hawaii Telescope, reaching surface brightnesses of $\sim$ 28.5 - 29 in the g-band. The dwarf candidates were identified through a direct visual inspection of the images and by visually cleaning a sample selected using a partially automated approach, and were morphologically classified at the time of identification. Approximately 75% of our candidates are dEs, indicating that a large number of early type dwarfs also populate low density environments, and 23.2% are nucleated. Distances were determined for 13.5% of our sample using pre-existing $z_{spec}$ measurements and HI detections. We confirm the dwarf nature for 99% of this sub-sample based on a magnitude cut $M_g$ = -18. Additionally, most of these ($\sim$90%) have relative velocities suggesting that they form a satellite population around nearby massive galaxies rather than an isolated field sample. Assuming that the candidates over the whole survey are satellites of the nearby galaxies, we demonstrate that the MATLAS dwarfs follow the same scaling relations as dwarfs in the Local Group as well as the Virgo and Fornax clusters. We also find that the nucleated fraction increases with $M_g$, and find evidence of a morphology-density relation for dwarfs around isolated massive galaxies. △ Less

Submitted 17 December, 2019; v1 submitted 29 October, 2019; originally announced October 2019.

Comments: 21 pages, 16 figures, 1 table. Accepted for publication in MNRAS; minor corrections to the text, and updated Figures 8 and 9 to match published version

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)

VEXAS: the VISTA EXtension to Auxiliary Surveys -- Data Release 1: the Southern Galactic Hemisphere

Authors: C. Spiniello, A. Agnello

Abstract: We present the first public data release of the VISTA EXtension to Auxiliary Surveys (VEXAS), comprising of 9 cross-matched multi-wavelength photometric catalogs where each object has a match in at least two surveys. We aim at a spatial coverage as uniform as possible in the multi-wavelength sky, with the purpose of providing the astronomical community with reference magnitudes and colours for var… ▽ More We present the first public data release of the VISTA EXtension to Auxiliary Surveys (VEXAS), comprising of 9 cross-matched multi-wavelength photometric catalogs where each object has a match in at least two surveys. We aim at a spatial coverage as uniform as possible in the multi-wavelength sky, with the purpose of providing the astronomical community with reference magnitudes and colours for various scientific uses, including: object classification (e.g. quasars, galaxies, and stars; high-z galaxies, white dwarfs, etc.); photometric redshifts of large galaxy samples; searches of exotic objects such as, for example, extremely red objects and lensed quasars. We have cross-matched the wide-field VISTA catalogs (the VISTA Hemisphere Survey and the VISTA Kilo Degree Infrared Galaxy Survey) with the AllWISE mid-infrared Survey, requiring that a match exists within 10 arcsec. We have further matched this table with X-Ray and radio data (ROSAT, XMM, SUMSS). We also performed a second cross-match between VISTA and AllWISE, with a smaller matching radius (3"), including WISE magnitudes. We have then cross-matched this resulting table ($\approx138\times10^6$ objects) with three photometric wide-sky optical deep surveys (DES, SkyMapper, PanSTARRS). We finally include matches to objects with spectroscopic follow-up by the SDSS and 6dFGS. To demonstrate the power of all-sky multi-wavelength cross-match tables, we show two examples of scientific applications of VEXAS, in particular using the publicly released tables to discover strong gravitational lenses (beyond the reach of previous searches), and to build a statistically large sample of extremely red objects. The VEXAS catalog is currently the widest and deepest, public, optical-to-IR photometric and spectroscopic database in the Southern Hemisphere. △ Less

Submitted 29 August, 2019; originally announced August 2019.

Comments: 13 pages, 10 figures, accepted for publication on Astronomy and Astrophysics

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

The Hubble-Lemaître constant and sound horizon from low-redshift probes

Authors: Radosław Wojtak, Adriano Agnello

Abstract: We revisit the claimed tension, or lack thereof, of measured values of the Hubble-Lemaître parameter $H_0$ from Cosmic Microwave Background (CMB) data and low-redshift indicators. Baryon Acoustic Oscillations (BAO) rely on the scale of the sound horizon at recombination $r_s$ to convert angular measurements into angular-diameter distances, so fixing $r_s$ from CMB measurements already constrains… ▽ More We revisit the claimed tension, or lack thereof, of measured values of the Hubble-Lemaître parameter $H_0$ from Cosmic Microwave Background (CMB) data and low-redshift indicators. Baryon Acoustic Oscillations (BAO) rely on the scale of the sound horizon at recombination $r_s$ to convert angular measurements into angular-diameter distances, so fixing $r_s$ from CMB measurements already constrains $H_0.$ If departures from concordance cosmology are to be constrained, truly independent measurements of $H_0$ are needed. We use the angular-diameter distances to three time-delay lenses from the H0LiCOW collaboration to calibrate the distance ladder, combine them with relative distances from Supernovae Ia and BAO, leaving $r_s$ completely free, and provide the inferred coefficients ($q_{0},j_{0},s_{0}$) in the polynomial expansion of H(z). We obtain $H_{0}r_{s}=(9895\pm161)$km/s and $H_0=(72\pm7)$km/s/Mpc. Combined with $H_0$ from H0LiCOW, then $r_s=(137\pm4.5)$Mpc is consistent with previous work and systematically lower than the CMB-inferred value. Our results are independent of the adopted cosmology, and removing Supernovae with z<0.1 has a negligible effect. △ Less

Submitted 6 August, 2019; originally announced August 2019.

Comments: 6 pages, 2 figures, 2 tables; published in MNRAS

Journal ref: MNRAS, 486, 5046 (2019)

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

KiDS-SQuaD II: Machine learning selection of bright extragalactic objects to search for new gravitationally lensed quasars

Authors: Vladislav Khramtsov, Alexey Sergeyev, Chiara Spiniello, Crescenzo Tortora, Nicola R. Napolitano, Adriano Agnello, Fedor Getman, Jelte T. A. de Jong, Konrad Kuijken, Mario Radovich, HuanYuan Shan, Valery Shulga

Abstract: The KiDS Strongly lensed QUAsar Detection project (KiDS-SQuaD) aims at finding as many previously undiscovered gravitational lensed quasars as possible in the Kilo Degree Survey. This is the second paper of this series where we present a new, automatic object classification method based on machine learning technique. The main goal of this paper is to build a catalogue of bright extragalactic objec… ▽ More The KiDS Strongly lensed QUAsar Detection project (KiDS-SQuaD) aims at finding as many previously undiscovered gravitational lensed quasars as possible in the Kilo Degree Survey. This is the second paper of this series where we present a new, automatic object classification method based on machine learning technique. The main goal of this paper is to build a catalogue of bright extragalactic objects (galaxies and quasars), from the KiDS Data Release 4, with a minimum stellar contamination, preserving the completeness as much as possible, to then apply morphological methods to select reliable gravitationally lensed (GL) quasar candidates. After testing some of the most used machine learning algorithms, decision trees based classifiers, we decided to use CatBoost, that was specifically trained with the aim of creating a sample of extragalactic sources as clean as possible from stars. We discuss the input data, define the training sample for the classifier, give quantitative estimates of its performances, and finally describe the validation results with Gaia DR2, AllWISE, and GAMA catalogues. We have built and make available to the scientific community the KiDS Bright EXtraGalactic Objects catalogue (KiDS-BEXGO), specifically created to find gravitational lenses. This is made of $\approx6$ millions of sources classified as quasars ($\approx 200\,000$) and galaxies ($\approx 5.7$M), up to $r<22^m$. From this catalog we selected 'Multiplets': close pairs of quasars or galaxies surrounded by at least one quasar, presenting the 12 most reliable gravitationally lensed quasar candidates, to demonstrate the potential of the catalogue, which will be further explored in a forthcoming paper. We compared our search to the previous one, presented in the first paper from this series, showing that employing a machine learning method decreases the stars-contaminators within the GL candidates. △ Less

Submitted 7 June, 2019; v1 submitted 4 June, 2019; originally announced June 2019.

Comments: Submitted to Astronomy and Astrophysics, 19 pages, 12 figures, 4 tables

Journal ref: A&A 632, A56 (2019)

On the mass mismatch between simulations and weak-lensing measurements

Authors: Jacob Svensmark, Davide Martizzi, Adriano Agnello

Abstract: The recently discovered discrepancy between galaxy mass measurements from weak lensing and predictions from abundance matching questions our understanding of cosmology, or of the galaxy-halo connection, or of both. We re-examined this tension by considering, as models, different cosmological simulations in the Illustris suite. We produced excess profiles $RΔΣ$ from subhalo snapshots at different r… ▽ More The recently discovered discrepancy between galaxy mass measurements from weak lensing and predictions from abundance matching questions our understanding of cosmology, or of the galaxy-halo connection, or of both. We re-examined this tension by considering, as models, different cosmological simulations in the Illustris suite. We produced excess profiles $RΔΣ$ from subhalo snapshots at different redshifts in Illustris-1 and IllustrisTNG (TNG100 and TNG300) simulations, enabling a direct comparison with weak-lensing measurements. We separate the individual contributions of stars, dark matter and gas within $\approx1$ Mpc (comoving length), beyond which correlated two-halo terms dominate. The mismatch between measurements and predictions is more severe than in previous studies: $RΔΣ$ profiles from IllustrisTNG are $\approx2$ times higher than the measured ones. Contrary to abundance matching results, the mismatch is mostly unchanged with increasing redshifts. The contribution of gas to the $RΔΣ$ profiles is $5-10\%$ over the scales dominated by one-halo terms. Different procedures to link stellar and halo masses (abundance matching, cosmological simulations) are still significantly discrepant with weak lensing measurements, but their trends are different. Therefore, the change in cosmological parameters advocated through abundance-matching arguments may not resolve this tension. Also, current criteria to select isolated massive galaxies in simulations are susceptible to resolution issues and may not correspond to observational criteria. The (currently subdominant) contribution of gas is non-negligible, and even if the major discrepancy within stellar and halo masses is resolved, it will be an appreciable source of systematics in the LSST era, when uncertainties on the $RΔΣ$ profiles are expected to be $\approx10$ times smaller. △ Less

Submitted 3 June, 2019; originally announced June 2019.

Comments: Submitted to A&A

Enhanced cluster lensing models with measured galaxy kinematics

Authors: P. Bergamini, P. Rosati, A. Mercurio, C. Grillo, G. B. Caminha, M. Meneghetti, A. Agnello, A. Biviano, F. Calura, C. Giocoli, M. Lombardi, G. Rodighiero, E. Vanzella

Abstract: We present an improved determination of the total mass distribution of three CLASH/HFF massive clusters, MACS J1206.2-0847 (z=0.44), MACS J0416.1-2403 (z=0.40), Abell S1063 (z=0.35). We specifically reconstruct the sub-halo mass component with robust stellar kinematics information of cluster galaxies, in combination with precise strong lensing models based on large samples of spectroscopically ide… ▽ More We present an improved determination of the total mass distribution of three CLASH/HFF massive clusters, MACS J1206.2-0847 (z=0.44), MACS J0416.1-2403 (z=0.40), Abell S1063 (z=0.35). We specifically reconstruct the sub-halo mass component with robust stellar kinematics information of cluster galaxies, in combination with precise strong lensing models based on large samples of spectroscopically identified multiple images. We use VLT/MUSE integral-field spectroscopy in the cluster cores to measure the stellar velocity dispersion, $σ$, of 40-60 member galaxies per cluster, covering 4-5 magnitudes to $m_{F160W}\simeq 21.5$. We verify the robustness and quantify the accuracy of the velocity dispersion measurements with extensive spectral simulations. With these data, we determine the normalization and slope of the galaxy $L\mbox{-}σ$ Faber-Jackson relation in each cluster and use these parameters as a prior for the scaling relations of the sub-halo population in the mass distribution modeling. When compared to our previous lens models, the inclusion of member galaxies' kinematics provides a similar precision in reproducing the positions of the multiple images. However, the inherent degeneracy between the central effective velocity dispersion, $σ_0$, and truncation radius, $r_{cut}$, of sub-halos is strongly reduced, thus significantly alleviating possible systematics in the measurements of sub-halo masses. The three independent determinations of the $σ_0\mbox{-}r_{cut}$ scaling relation in each cluster are found to be fully consistent, enabling a statistical determination of sub-halo sizes as a function of $σ_0$, or halo masses. We derive galaxy central velocity dispersion functions of the three clusters and found them in agreement with each other. Sub-halo mass functions determined with this method can be compared with those obtained from cosmological simulations. △ Less

Submitted 4 October, 2019; v1 submitted 30 May, 2019; originally announced May 2019.

Comments: v2: minor changes. Accepted for publication in A&A (18 pages, 11 figures)

Journal ref: A&A 631, A130 (2019)

Low-redshift measurement of the sound horizon through gravitational time-delays

Authors: Nikki Arendse, Adriano Agnello, Radosław Wojtak

Abstract: The matter sound horizon can be inferred from the cosmic microwave background within the Standard Model. Independent direct measurements of the sound horizon are then a probe of possible deviations from the Standard Model. We aim at measuring the sound horizon $r_s$ from low-redshift indicators, which are completely independent of CMB inference. We used the measured product $H(z)r_s$ from baryon a… ▽ More The matter sound horizon can be inferred from the cosmic microwave background within the Standard Model. Independent direct measurements of the sound horizon are then a probe of possible deviations from the Standard Model. We aim at measuring the sound horizon $r_s$ from low-redshift indicators, which are completely independent of CMB inference. We used the measured product $H(z)r_s$ from baryon acoustic oscillations (BAO) together with supernovae~\textsc{I}a to constrain $H(z)/H_{0}$ and time-delay lenses analysed by the H0LiCOW collaboration to anchor cosmological distances ($\propto H_{0}^{-1}$). {Additionally, we investigated the influence of adding a sample of quasars with higher redshift with standardisable UV-Xray luminosity distances. We adopted polynomial expansions in $H(z)$ or in comoving distances} so that our inference was completely independent of any cosmological model on which the expansion history might be based. Our measurements are independent of Cepheids and systematics from peculiar motions {to within percent-level accuracy.} The inferred sound horizon $r_s$ varies between $(133 \pm 8)$~Mpc and $(138 \pm 5)$~Mpc across different models. The discrepancy with CMB measurements is robust against model choice. Statistical uncertainties are comparable to systematics. The combination of time-delay lenses, supernovae, and BAO yields a distance ladder that is independent of cosmology (and of Cepheid calibration) and a measurement of $r_s $ that is independent of the CMB. These cosmographic measurements are then a competitive test of the Standard Model, regardless of the hypotheses on which the cosmology is based. △ Less

Submitted 19 February, 2020; v1 submitted 28 May, 2019; originally announced May 2019.

Comments: 2019A&A...632A..91A (subm. 28/05/2019), 6 pages, 3 figures

Journal ref: Astronomy & Astrophysics, 2019, Volume 632, id.A91, 6 pp. Astronomy & Astrophysics, Volume 632, id.A91, 6 pp

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

Spectroscopic confirmation and modelling of two lensed quadruple quasars in the Dark Energy Survey public footprint

Authors: C. Spiniello, A. V. Sergeyev, L. Marchetti, C. Tortora, N. R. Napolitano, V. Shalyapin, A. Agnello, F. I. Getman, M. Vaccari, S. Serjeant, L. V. E. Koopmans, A. J. Baker, T. H. Jarrett, G. Covone, G. Vernardos

Abstract: Quadruply lensed quasars are extremely rare objects, but incredibly powerful cosmological tools. Only few dozen are known in the whole sky. Here we present the spectroscopic confirmation of two new quadruplets WG0214-2105 and WG2100-4452 discovered by Agnello & Spiniello (2018) within the Dark Energy Survey (DES) public footprints. We have conducted spectroscopic follow-up of these systems with th… ▽ More Quadruply lensed quasars are extremely rare objects, but incredibly powerful cosmological tools. Only few dozen are known in the whole sky. Here we present the spectroscopic confirmation of two new quadruplets WG0214-2105 and WG2100-4452 discovered by Agnello & Spiniello (2018) within the Dark Energy Survey (DES) public footprints. We have conducted spectroscopic follow-up of these systems with the Southern African Large Telescope as part of a program that aims at confirming the largest possible number of optically selected strong gravitational lensing systems in the Equatorial and Southern Hemisphere. For both systems, we present the spectra for the sources and deflectors that allowed us to estimate the source redshifts and unambiguously confirm their lensing nature. For the brighter deflector (WG2100-4452), we measure the stellar velocity dispersion from the spectrum. We also obtain photometry for both lenses, directly from DES multi-band images, isolating the lens galaxies from the quasar images. One of the quadruplets, WG0214-2105, was also observed by Pan-STARRS, allowing us to estimate the apparent brightness of each quasar image at two different epochs, and thus to find evidence for flux variability. This result could suggest a microlensing event for the faintest components, although intrinsic variability cannot be excluded with only two epochs. Finally, we present simple lens models for both quadruplets, obtaining Einstein radii, SIE velocity dispersions, ellipticities, and position angles of the lens systems, as well as time delay predictions assuming a concordance cosmological model. △ Less

Submitted 16 January, 2019; v1 submitted 14 January, 2019; originally announced January 2019.

Comments: Submitted to MNRAS, 10 pages, 6 figures, 4 tables. Some typos have been corrected with respect to the previous version

The ultra-diffuse galaxy NGC 1052-DF2 with MUSE: II. The population of DF2: stars, clusters and planetary nebulae

Authors: Jeremy Fensch, Remco F. J. van der Burg, Tereza Jerabkova, Eric Emsellem, Anita Zanella, Adriano Agnello, Michael Hilker, Oliver Mueller, Marina Rejkuba, Pierre-Alain Duc, Patrick Durrell, Rebecca Habas, Sungsoon Lim, Francine R. Marleau, Eric Peng, Ruben Sanchez Janssen

Abstract: NGC 1052-DF2, an ultra diffuse galaxy (UDG), has been the subject of intense debate. Its alleged absence of dark matter, and the brightness and number excess of its globular clusters (GCs) at an initially assumed distance of 20Mpc, suggested a new formation channel for UDGs. We present the first systematic spectroscopic analysis of both the stellar body and the GCs (six of which were previously kn… ▽ More NGC 1052-DF2, an ultra diffuse galaxy (UDG), has been the subject of intense debate. Its alleged absence of dark matter, and the brightness and number excess of its globular clusters (GCs) at an initially assumed distance of 20Mpc, suggested a new formation channel for UDGs. We present the first systematic spectroscopic analysis of both the stellar body and the GCs (six of which were previously known, and one newly confirmed member) of this galaxy using MUSE@VLT. Even though NGC 1052-DF2 does not show any spatially extended emission lines we report the discovery of three planetary nebulae (PNe). We conduct full spectral fitting on the UDG and the stacked spectra of all GCs. The UDG's stellar population is old, 8.9$\pm$1.5 Gyr, metal-poor, with [M/H] = $-$1.07$\pm$0.12 with little or no $α$-enrichment. The stacked spectrum of all GCs indicates a similar age of 8.9$\pm$1.8 Gyr, but lower metallicity, with [M/H] = $-$1.63$\pm$0.09, and similarly low $α$-enrichment. There is no evidence for a variation of age and metallicity in the GC population with the available spectra. The significantly more metal-rich stellar body with respect to its associated GCs, the age of the population, its metallicity and alpha enrichment, are all in line with other dwarf galaxies. NGC 1052-DF2 thus falls on the same empirical mass-metallicity relation as other dwarfs, for the full distance range assumed in the literature. We find that both debated distance estimates (13 and 20 Mpc) are similarly likely, given the three discovered PNe. △ Less

Submitted 20 January, 2019; v1 submitted 18 December, 2018; originally announced December 2018.

Comments: 10 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 625, A77 (2019)

The ultra-diffuse galaxy NGC 1052-DF2 with MUSE: I. Kinematics of the stellar body

Authors: Eric Emsellem, Remco F. J. van der Burg, Jeremy Fensch, Tereza Jerabkova, Anita Zanella, Adriano Agnello, Michael Hilker, Oliver Mueller, Marina Rejkuba, Pierre-Alain Duc, Patrick Durrell, Rebecca Habas, Federico Lelli, Sungsoon Lim, Francine R. Marleau, Eric Peng, Ruben Sanchez-Janssen

Abstract: The so-called ultra-diffuse galaxy NGC~1052-DF2 was announced to be a galaxy lacking dark matter based on a spectroscopic study of its constituent globular clusters. Here we present the first spectroscopic analysis of the stellar body of this galaxy using the MUSE integral-field spectrograph at the (ESO) Very Large Telescope. The MUSE datacube simultaneously provides DF2's stellar velocity field a… ▽ More The so-called ultra-diffuse galaxy NGC~1052-DF2 was announced to be a galaxy lacking dark matter based on a spectroscopic study of its constituent globular clusters. Here we present the first spectroscopic analysis of the stellar body of this galaxy using the MUSE integral-field spectrograph at the (ESO) Very Large Telescope. The MUSE datacube simultaneously provides DF2's stellar velocity field and systemic velocities for seven globular clusters (GCs). We further discovered three planetary nebulae (PNe) that are likely part of this galaxy. While five of the clusters had velocities measured in the literature, we were able to confirm the membership of two more candidates through precise radial velocity measurements, which increases the measured specific frequency of GCs in DF2. The mean velocity of the diffuse stellar body, 1792.9$^{-1.8}_{+1.4}$~\kms, is consistent with the mean globular cluster velocity. We detect a weak but significant velocity gradient within the stellar body, with a kinematic axis close to the photometric major-axis, making it a prolate-like rotator. We estimate a velocity dispersion from the clusters and PNe of $σ_{\mathrm{int}}=10.6^{+3.9}_{-2.3}$~\kms. The velocity dispersion $σ_{\rm{DF2}\star}$(\re) for the stellar body within one effective radius is $10.8^{-4.0}_{+3.2}$~\kms. Considering various sources of systemic uncertainties this central value varies between 5 and 13~\kms, and we conservatively report a 95\% confidence upper limit to the dispersion within one \re\ of 21~\kms. We provide updated mass estimates based on these dispersions corresponding to the different distances to NGC~1052-DF2 that have been reported in the recent literature. △ Less

Submitted 1 April, 2019; v1 submitted 18 December, 2018; originally announced December 2018.

Comments: 10 pages, 5 figures, Accepted for publication in A&A. Publisher ready version: minor language editing

Bright lenses are easy to find: Spectroscopic confirmation of lensed quasars in the Southern Sky

Authors: C. Spiniello, A. Agnello, A. V. Sergeyev, T. Anguita, Ó. Rodríguez, N. R. Napolitano, C. Tortora

Abstract: Gravitationally lensed quasars are valuable, but extremely rare, probes of observational cosmology and extragalactic astrophysics. Progress in these fields has been limited just by the paucity of systems with good ancillary data. Here we present a first spectroscopic confirmation of lenses discovered in the Southern Sky from the DES and KiDS-DR3 footprints. % optical ground based surveys in the So… ▽ More Gravitationally lensed quasars are valuable, but extremely rare, probes of observational cosmology and extragalactic astrophysics. Progress in these fields has been limited just by the paucity of systems with good ancillary data. Here we present a first spectroscopic confirmation of lenses discovered in the Southern Sky from the DES and KiDS-DR3 footprints. % optical ground based surveys in the Southern Hemisphere. We have targeted 7 high-graded candidates, selected with new techniques, with NTT-EFOSC2, and confirmed 5 of them. We provide source spectroscopic redshifts, image separations, $gri$ photometry and first lens model parameters. The success rate of ~70% confirms our forecasts, based on the comparison between the number of candidate doubles and quadruplets in our searches over a ~5000 sq.deg footprint and theoretical predictions. △ Less

Submitted 7 December, 2018; v1 submitted 2 November, 2018; originally announced November 2018.

Comments: 6 pages, 5 figures. Accepted for publication to MNRAS

KiDS0239-3211: A new gravitational quadruple lens candidate

Authors: A. Sergeyev, C. Spiniello, V. Khramtsov, N. R. Napolitano, E. Bannikova, C. Tortora, F. I. Getman, A. Agnello

Abstract: We report the discovery of a candidate to quadrupole gravitationally lensed system KiDS0239-3211 based on the public data release 3 of the KiDS survey and machine learning techniques. We report the discovery of a candidate to quadrupole gravitationally lensed system KiDS0239-3211 based on the public data release 3 of the KiDS survey and machine learning techniques. △ Less

Submitted 3 October, 2018; originally announced October 2018.

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