Showing 1–50 of 79 results for author: Vegetti, S

General Multipoles and Their Implications for Dark Matter Inference

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

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

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

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

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

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

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

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

Comments: In print with A&A

A Collection of German Science Interests in the Next Generation Very Large Array

Authors: M. Kadler, D. A. Riechers, J. Agarwal, A. -K. Baczko, H. Beuther, F. Bigiel, T. Birnstiel, B. Boccardi, D. J. Bomans, L. Boogaard, T. T. Braun, S. Britzen, M. Brüggen, A. Brunthaler, P. Caselli, D. Elsässer, S. von Fellenberg, M. Flock, C. M. Fromm, L. Fuhrmann, P. Hartogh, M. Hoeft, R. P. Keenan, Y. Kovalev, K. Kreckel , et al. (66 additional authors not shown)

Abstract: The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy com… ▽ More The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy community under the lead of the National Radio Astronomy Observatory (NRAO), and is anticipated to be built between 2027 and 2037. Two workshops have been held in 2022 and 2023 with the goal to discuss and consolidate the scientific interests in the ngVLA within the German astronomical community. This community paper constitutes a collection of 48 science ideas which the German community aims to pursue with the ngVLA in the 2030s. This is not a complete list and the ideas are not developed at the level of a "Science Book", such that the present document is mainly meant provide a basis for further discussion within the community. As such, additional contributions are welcome, and will be considered for inclusion in future revisions. △ Less

Submitted 18 June, 2024; v1 submitted 16 November, 2023; originally announced November 2023.

Comments: Version 2.0 (status June 18, 2024): 169 pages, comments and future contributions welcome [v2.0: 7 new science cases added, some minor revisions to other chapters]

Angular complexity in strong lens substructure detection

Authors: Conor M. O'Riordan, Simona Vegetti

Abstract: Strong gravitational lensing can be used to find otherwise invisible dark matter subhaloes. In such an analysis, the lens galaxy mass model is a significant source of systematic uncertainty. In this paper we analyse the effect of angular complexity in the lens model. We use multipole perturbations which introduce low-order deviations from pure ellipticity in the isodensity contours, keeping the ra… ▽ More Strong gravitational lensing can be used to find otherwise invisible dark matter subhaloes. In such an analysis, the lens galaxy mass model is a significant source of systematic uncertainty. In this paper we analyse the effect of angular complexity in the lens model. We use multipole perturbations which introduce low-order deviations from pure ellipticity in the isodensity contours, keeping the radial density profile fixed. We find that, in HST-like data, multipole perturbations consistent with those seen in galaxy isophotes are very effective at causing false positive substructure detections. We show that the effectiveness of this degeneracy depends on the deviation from a pure ellipse and the lensing configuration. We find that, when multipoles of one per cent are allowed in the lens model, the area in the observation where a subhalo could be detected drops by a factor of three. Sensitivity away from the lensed images is mostly lost. However, the mass limit of detectable objects on or close to the lensed images does not change. We do not expect the addition of multipole perturbations to lens models to have a significant effect on the ability of strong lensing to constrain the underlying dark matter model. However, given the high rate of false positive detections, angular complexity beyond the elliptical power-law should be included for such studies to be reliable. We discuss implications for previous detections and future work. △ Less

Submitted 12 January, 2024; v1 submitted 16 October, 2023; originally announced October 2023.

Comments: 11 pages, 7 figures, accepted by MNRAS

Selection functions of strong lens finding neural networks

Authors: A. Herle, C. M. O'Riordan, S. Vegetti

Abstract: Convolution Neural Networks trained for the task of lens finding with similar architecture and training data as is commonly found in the literature are biased classifiers. An understanding of the selection function of lens finding neural networks will be key to fully realising the potential of the large samples of strong gravitational lens systems that will be found in upcoming wide-field surveys.… ▽ More Convolution Neural Networks trained for the task of lens finding with similar architecture and training data as is commonly found in the literature are biased classifiers. An understanding of the selection function of lens finding neural networks will be key to fully realising the potential of the large samples of strong gravitational lens systems that will be found in upcoming wide-field surveys. We use three training datasets, representative of those used to train galaxy-galaxy and galaxy-quasar lens finding neural networks. The networks preferentially select systems with larger Einstein radii and larger sources with more concentrated source-light distributions. Increasing the detection significance threshold to 12$σ$ from 8$σ$ results in 50 per cent of the selected strong lens systems having Einstein radii $θ_\mathrm{E}$ $\ge$ 1.04 arcsec from $θ_\mathrm{E}$ $\ge$ 0.879 arcsec, source radii $R_S$ $\ge$ 0.194 arcsec from $R_S$ $\ge$ 0.178 arcsec and source Sérsic indices $n_{\mathrm{Sc}}^{\mathrm{S}}$ $\ge$ 2.62 from $n_{\mathrm{Sc}}^{\mathrm{S}}$ $\ge$ 2.55. The model trained to find lensed quasars shows a stronger preference for higher lens ellipticities than those trained to find lensed galaxies. The selection function is independent of the slope of the power-law of the mass profiles, hence measurements of this quantity will be unaffected. The lens finder selection function reinforces that of the lensing cross-section, and thus we expect our findings to be a general result for all galaxy-galaxy and galaxy-quasar lens finding neural networks. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: Submitted to MNRAS

Strong gravitational lensing as a probe of dark matter

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

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

Submitted 20 June, 2023; originally announced June 2023.

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

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

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

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

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

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

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

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

Submitted 1 February, 2023; originally announced February 2023.

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

Solving Inverse Physics Problems with Score Matching

Authors: Benjamin J. Holzschuh, Simona Vegetti, Nils Thuerey

Abstract: We propose to solve inverse problems involving the temporal evolution of physics systems by leveraging recent advances from diffusion models. Our method moves the system's current state backward in time step by step by combining an approximate inverse physics simulator and a learned correction function. A central insight of our work is that training the learned correction with a single-step loss i… ▽ More We propose to solve inverse problems involving the temporal evolution of physics systems by leveraging recent advances from diffusion models. Our method moves the system's current state backward in time step by step by combining an approximate inverse physics simulator and a learned correction function. A central insight of our work is that training the learned correction with a single-step loss is equivalent to a score matching objective, while recursively predicting longer parts of the trajectory during training relates to maximum likelihood training of a corresponding probability flow. We highlight the advantages of our algorithm compared to standard denoising score matching and implicit score matching, as well as fully learned baselines for a wide range of inverse physics problems. The resulting inverse solver has excellent accuracy and temporal stability and, in contrast to other learned inverse solvers, allows for sampling the posterior of the solutions. △ Less

Submitted 5 December, 2023; v1 submitted 24 January, 2023; originally announced January 2023.

Comments: Accepted at NeurIPS 2023; code available at https://github.com/tum-pbs/SMDP

Sensitivity of strong lensing observations to dark matter substructure: a case study with Euclid

Authors: Conor M. O'Riordan, Giulia Despali, Simona Vegetti, Mark R. Lovell, Ángeles Moliné

Abstract: We introduce a machine learning method for estimating the sensitivity of strong lens observations to dark matter subhaloes in the lens. Our training data include elliptical power-law lenses, Hubble Deep Field sources, external shear, and noise and PSF for the Euclid VIS instrument. We set the concentration of the subhaloes using a $v_\mathrm{max}$-$r_\mathrm{max}$ relation. We then estimate the da… ▽ More We introduce a machine learning method for estimating the sensitivity of strong lens observations to dark matter subhaloes in the lens. Our training data include elliptical power-law lenses, Hubble Deep Field sources, external shear, and noise and PSF for the Euclid VIS instrument. We set the concentration of the subhaloes using a $v_\mathrm{max}$-$r_\mathrm{max}$ relation. We then estimate the dark matter subhalo sensitivity in $16{,}000$ simulated strong lens observations with depth and resolution resembling Euclid VIS images. We find that, with a $3σ$ detection threshold, $2.35$ per cent of pixels inside twice the Einstein radius are sensitive to subhaloes with a mass $M_\mathrm{max}\leq 10^{10}M_\odot$, $0.03$ per cent are sensitive to $M_\mathrm{max}\leq 10^{9}M_\odot$, and, the limit of sensitivity is found to be $M_\mathrm{max}=10^{8.8\pm0.2}M_\odot$. Using our sensitivity maps and assuming CDM, we estimate that Euclid-like lenses will yield $1.43^{+0.14}_{-0.11}[f_\mathrm{sub}^{-1}]$ detectable subhaloes per lens in the entire sample, but this increases to $35.6^{+0.9}_{-0.9}[f_\mathrm{sub}^{-1}]$ per lens in the most sensitive lenses. Estimates are given in units of the inverse of the substructure mass fraction $f_\mathrm{sub}^{-1}$. Assuming $f_\mathrm{sub}=0.01$, one in every $70$ lenses in general should yield a detection, or one in every $\sim$ three lenses in the most sensitive sample. From $170,000$ new strong lenses detected by Euclid, we expect $\sim 2500$ new subhalo detections. We find that the expected number of detectable subhaloes in warm dark matter models only changes relative to cold dark matter for models which have already been ruled out, i.e., those with half-mode masses $M_\mathrm{hm}>10^8M_\odot$. △ Less

Submitted 27 February, 2023; v1 submitted 28 November, 2022; originally announced November 2022.

Comments: 15 pages, 14 figures, accepted by MNRAS

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

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

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

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

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

Constraining SIDM with halo shapes: revisited predictions from realistic simulations of early-type galaxies

Authors: Giulia Despali, Levi G. Walls, Simona Vegetti, Martin Sparre, Mark Vogelsberger, Jesús Zavala

Abstract: We study the effect of self-interacting dark matter (SIDM) and baryons on the shape of early-type galaxies (ETGs) and their dark matter haloes, comparing them to the predictions of the cold dark matter (CDM) scenario. We use five hydrodynamical zoom-in simulations of haloes hosting ETGs ($M_{\rm vir}\sim 10^{13}M_{\odot}$ and $M_{*}\sim10^{11}M_{\odot}$), simulated in CDM and a SIDM model with con… ▽ More We study the effect of self-interacting dark matter (SIDM) and baryons on the shape of early-type galaxies (ETGs) and their dark matter haloes, comparing them to the predictions of the cold dark matter (CDM) scenario. We use five hydrodynamical zoom-in simulations of haloes hosting ETGs ($M_{\rm vir}\sim 10^{13}M_{\odot}$ and $M_{*}\sim10^{11}M_{\odot}$), simulated in CDM and a SIDM model with constant cross-section of $σ_T/m_χ= 1\ \mathrm{cm}^2 \mathrm{g}^{-1}$. We measure the three-dimensional and projected shapes of the dark matter haloes and their baryonic content using the inertia tensor and compare our measurements to the results of three $HST$ samples of gravitational lenses and $Chandra$ and $XMM-Newton$ X-ray observations. We find that the inclusion of baryons greatly reduces the differences between CDM and a SIDM, together with the ability to draw constraints based on shapes. Lensing measurements reject the predictions of CDM dark-matter-only simulations and prefer one of the hydro scenarios. When we consider the total sample of lenses, observational data prefer the CDM hydro scenario. The shapes of the X-ray emitting gas are compatible with observational results in both hydro runs, with CDM predicting higher elongations only in the very centre. Contrary to previous claims at the scale of elliptical galaxies, we conclude that both CDM and our SIDM model can still explain observed shapes once we include baryons in the simulations. Our results demonstrate that this is essential to derive realistic constraints and that new simulations are needed to confirm and extend our findings. △ Less

Submitted 5 September, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: 17 pages, 13 figures, revised version accepted for publication in MNRAS

Realistic galaxy images and improved robustness in machine learning tasks from generative modelling

Authors: Benjamin J. Holzschuh, Conor M. O'Riordan, Simona Vegetti, Vicente Rodriguez-Gomez, Nils Thuerey

Abstract: We examine the capability of generative models to produce realistic galaxy images. We show that mixing generated data with the original data improves the robustness in downstream machine learning tasks. We focus on three different data sets; analytical Sérsic profiles, real galaxies from the COSMOS survey, and galaxy images produced with the SKIRT code, from the IllustrisTNG simulation. We quantif… ▽ More We examine the capability of generative models to produce realistic galaxy images. We show that mixing generated data with the original data improves the robustness in downstream machine learning tasks. We focus on three different data sets; analytical Sérsic profiles, real galaxies from the COSMOS survey, and galaxy images produced with the SKIRT code, from the IllustrisTNG simulation. We quantify the performance of each generative model using the Wasserstein distance between the distributions of morphological properties (e.g. the Gini-coefficient, the asymmetry, and ellipticity), the surface brightness distribution on various scales (as encoded by the power-spectrum), the bulge statistic and the colour for the generated and source data sets. With an average Wasserstein distance (Fréchet Inception Distance) of $7.19 \times 10^{-2}\, (0.55)$, $5.98 \times 10^{-2}\, (1.45)$ and $5.08 \times 10^{-2}\, (7.76)$ for the Sérsic, COSMOS and SKIRT data set, respectively, our best models convincingly reproduce even the most complicated galaxy properties and create images that are visually indistinguishable from the source data. We demonstrate that by supplementing the training data set with generated data, it is possible to significantly improve the robustness against domain-shifts and out-of-distribution data. In particular, we train a convolutional neural network to denoise a data set of mock observations. By mixing generated images into the original training data, we obtain an improvement of $11$ and $45$ per cent in the model performance regarding domain-shifts in the physical pixel size and background noise level, respectively. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: 33 pages, 21 figures, submitted to MNRAS, comments welcome

Astrophysical and Cosmological Probes of Dark Matter

Authors: Kimberly K. Boddy, Mariangela Lisanti, Samuel D. McDermott, Nicholas L. Rodd, Christoph Weniger, Yacine Ali-Haïmoud, Malte Buschmann, Ilias Cholis, Djuna Croon, Adrienne L. Erickcek, Vera Gluscevic, Rebecca K. Leane, Siddharth Mishra-Sharma, Julian B. Muñoz, Ethan O. Nadler, Priyamvada Natarajan, Adrian Price-Whelan, Simona Vegetti, Samuel J. Witte

Abstract: While astrophysical and cosmological probes provide a remarkably precise and consistent picture of the quantity and general properties of dark matter, its fundamental nature remains one of the most significant open questions in physics. Obtaining a more comprehensive understanding of dark matter within the next decade will require overcoming a number of theoretical challenges: the groundwork for t… ▽ More While astrophysical and cosmological probes provide a remarkably precise and consistent picture of the quantity and general properties of dark matter, its fundamental nature remains one of the most significant open questions in physics. Obtaining a more comprehensive understanding of dark matter within the next decade will require overcoming a number of theoretical challenges: the groundwork for these strides is being laid now, yet much remains to be done. Chief among the upcoming challenges is establishing the theoretical foundation needed to harness the full potential of new observables in the astrophysical and cosmological domains, spanning the early Universe to the inner portions of galaxies and the stars therein. Identifying the nature of dark matter will also entail repurposing and implementing a wide range of theoretical techniques from outside the typical toolkit of astrophysics, ranging from effective field theory to the dramatically evolving world of machine learning and artificial-intelligence-based statistical inference. Through this work, the theory frontier will be at the heart of dark matter discoveries in the upcoming decade. △ Less

Submitted 12 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021, 24 pages, 4 figures, comments welcome

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

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

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

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

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

The inner density profile of an elliptical galaxy at z=1.15 from gravitational lensing

Authors: H. R. Stacey, C. M. O'Riordan, S. Vegetti, D. M. Powell, M. W. Auger, G. Despali

Abstract: The density profiles of lensing galaxies are typically parameterised by singular power-law models with a logarithmic slope close to isothermal ($ζ=2$). This is sufficient to fit the lensed emission near the Einstein radius but may not be sufficient when extrapolated to smaller or larger radii if the large-scale density profile is more complex. Here, we consider a broken power-law model for the den… ▽ More The density profiles of lensing galaxies are typically parameterised by singular power-law models with a logarithmic slope close to isothermal ($ζ=2$). This is sufficient to fit the lensed emission near the Einstein radius but may not be sufficient when extrapolated to smaller or larger radii if the large-scale density profile is more complex. Here, we consider a broken power-law model for the density profile of an elliptical galaxy at $z=1.15$ using observations with the Atacama Large (sub-)Millimetre Array of the strong gravitational lens system SPT0532$-$50. This is the first application of such a model to real data. We find the lensed emission is best fit by a density profile that is sub-isothermal ($ζ= 1.87^{+0.02}_{-0.03}$) near the Einstein radius and steepens to super-isothermal ($ζ= 2.14^{+0.03}_{-0.02}$) at around half the Einstein radius, demonstrating that the lensing data probes the mass distribution inside the region probed by the lensed images. Assuming that a broken power-law is the underlying truth, we find that a single power-law would result in a $10\pm1$ percent underestimate of the Hubble constant from time-delay cosmography. Our results suggest that a broken power-law could be useful for precision lens modelling and probing the structural evolution of elliptical galaxies. △ Less

Submitted 12 October, 2021; v1 submitted 22 September, 2021; originally announced September 2021.

Comments: Withdrawn pending further testing

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

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

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

Submitted 7 September, 2021; originally announced September 2021.

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

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

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

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

Submitted 21 July, 2021; originally announced July 2021.

Comments: 11 pages, 7 figures

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

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

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

Submitted 31 March, 2021; originally announced March 2021.

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

Journal ref: MNRAS, 463, 3143 (2016)

Dynamical properties of z $\sim 4.5$ dusty star-forming galaxies and their connection with local early type galaxies

Authors: Francesca Rizzo, Simona Vegetti, Filippo Fraternali, Hannah Stacey, Devon Powell

Abstract: There is a large consensus that gas in high-$z$ galaxies is highly turbulent, because of a combination of stellar feedback processes and gravitational instabilities driven by mergers and gas accretion. In this paper, we present the analysis of a sample of five Dusty Star Forming Galaxies (DSFGs) at $4 \lesssim z\lesssim 5$. Taking advantage of the magnifying power of strong gravitational lensing,… ▽ More There is a large consensus that gas in high-$z$ galaxies is highly turbulent, because of a combination of stellar feedback processes and gravitational instabilities driven by mergers and gas accretion. In this paper, we present the analysis of a sample of five Dusty Star Forming Galaxies (DSFGs) at $4 \lesssim z\lesssim 5$. Taking advantage of the magnifying power of strong gravitational lensing, we quantified their kinematic and dynamical properties from ALMA observations of their [CII] emission line. We combined the dynamical measurements obtained for these galaxies with those obtained from previous studies to build the largest sample of $z \sim 4.5$ galaxies with high-quality data and sub-kpc spatial resolutions, so far. We found that all galaxies in the sample are dynamically cold, with rotation-to-random motion ratios, $V/σ$, between 7 to 15. The relation between their velocity dispersions and their star-formation rates indicates that stellar feedback is sufficient to sustain the turbulence within these galaxies and no further mechanisms are needed. In addition, we performed a rotation curve decomposition to infer the relative contribution of the baryonic (gas, stars) and dark matter components to the total gravitational potentials. This analysis allowed us to compare the structural properties of the studied DSFGs with those of their descendants, the local early type galaxies. In particular, we found that five out of six galaxies of the sample show the dynamical signature of a bulge, indicating that the spheroidal component is already in place at $z \sim 4.5$. △ Less

Submitted 11 August, 2021; v1 submitted 10 February, 2021; originally announced February 2021.

Comments: 34 pages (including appendices), 10 figures, 8 tables; accepted for publication in MNRAS

An unexpected high concentration for the dark substructure in the gravitational lens SDSSJ0946+1006

Authors: Quinn E. Minor, Sophia Gad-Nasr, Manoj Kaplinghat, Simona Vegetti

Abstract: The presence of an invisible substructure has previously been detected in the gravitational lens galaxy SDSSJ0946+1006 through its perturbation of the lensed images. Using flexible models for the main halo and the subhalo perturbation to fit the lensed images, we demonstrate that the subhalo has an extraordinarily high central density and steep density slope. The inferred concentration for the sub… ▽ More The presence of an invisible substructure has previously been detected in the gravitational lens galaxy SDSSJ0946+1006 through its perturbation of the lensed images. Using flexible models for the main halo and the subhalo perturbation to fit the lensed images, we demonstrate that the subhalo has an extraordinarily high central density and steep density slope. The inferred concentration for the subhalo is well above the expected scatter in concentrations for $Λ$CDM halos of similar mass. We robustly infer the subhalo's projected mass within 1 kpc to be $\sim 2$-$3.7\times 10^9$M$_\odot$ at $>$95% CL for all our lens models, while the average slope of the subhalo's projected density profile over the radial range 0.75-1.25 kpc is constrained to be steeper than isothermal ($γ_{2D} \lesssim -1$). By modeling the subhalo light directly, we infer a conservative upper bound on its luminosity $L_V < 1.2\times 10^8L_\odot$ at 95% CL, which shows that the perturber is dark matter dominated. To compare to $Λ$CDM expectations, we analyze subhalos within analogues of lensing galaxies in the Illustris TNG100-1 simulation over many lines of sight, and find hundreds of subhalos that achieve a projected mass within 1 kpc of $\gtrsim 2\times10^9M_\odot$. However, less than 1% of the mock observations yield a log-slope steep enough to be consistent with our lensing models, and they $all$ have stellar masses in excess of that allowed by observations by about an order of magnitude or more. We conclude that the presence of such a dark, highly concentrated subhalo is unexpected in a $Λ$CDM universe. Finally, we show that this tension with CDM is not significantly reduced if the perturber is assumed to be a line-of-sight structure, rather than a subhalo. △ Less

Submitted 20 November, 2020; originally announced November 2020.

Comments: 24 pages, 18 figures, 3 appendices

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

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

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

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

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

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

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

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

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

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

Comments: Accepted by MNRAS; 22 pages

A dynamically cold disk galaxy in the early Universe

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

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

Submitted 2 September, 2020; originally announced September 2020.

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

Journal ref: Nature 584, 2020, 201-204

Time Delay Lens Modelling Challenge

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

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

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

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

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

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

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

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

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

Comments: Accepted in MNRAS

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

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

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

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

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

Systematic errors in strong gravitational lensing reconstructions, a numerical simulation perspective

Authors: W. Enzi, S. Vegetti, G. Despali, J. -W. Hsueh, R. B. Metcalf

Abstract: We present the analysis of a sample of twenty-four SLACS-like galaxy-galaxy strong gravitational lens systems with a background source and deflectors from the Illustris-1 simulation. We study the degeneracy between the complex mass distribution of the lenses, substructures, the surface brightness distribution of the sources, and the time delays. Using a novel inference framework based on Approxima… ▽ More We present the analysis of a sample of twenty-four SLACS-like galaxy-galaxy strong gravitational lens systems with a background source and deflectors from the Illustris-1 simulation. We study the degeneracy between the complex mass distribution of the lenses, substructures, the surface brightness distribution of the sources, and the time delays. Using a novel inference framework based on Approximate Bayesian Computation, we find that for all the considered lens systems, an elliptical and cored power-law mass density distribution provides a good fit to the data. However, the presence of cores in the simulated lenses affects most reconstructions in the form of a Source Position Transformation. The latter leads to a systematic underestimation of the source sizes by 50 per cent on average, and a fractional error in $H_{0}$ of around $25_{-19}^{+37}$ per cent. The analysis of a control sample of twenty-four lens systems, for which we have perfect knowledge about the shape of the lensing potential, leads to a fractional error on $H_{0}$ of $12_{-3}^{+6}$ per cent. We find no degeneracy between complexity in the lensing potential and the inferred amount of substructures. We recover an average total projected mass fraction in substructures of $f_{\rm sub}<1.7-2.0\times10^{-3}$ at the 68 per cent confidence level in agreement with zero and the fact that all substructures had been removed from the simulation. Our work highlights the need for higher-resolution simulations to quantify the lensing effect of more realistic galactic potentials better, and that additional observational constraint may be required to break existing degeneracies. △ Less

Submitted 4 May, 2020; v1 submitted 6 November, 2019; originally announced November 2019.

Comments: Accepted by MNRAS

The lensing properties of subhaloes in massive elliptical galaxies in sterile neutrino cosmologies

Authors: Giulia Despali, Mark Lovell, Simona Vegetti, Robert A. Crain, Benjamin D. Oppenheimer

Abstract: We use high-resolution hydrodynamical simulations run with the EAGLE model of galaxy formation to study the differences between the properties of - and subsequently the lensing signal from - subhaloes of massive elliptical galaxies at redshift 0.2, in Cold and Sterile Neutrino (SN) Dark matter models. We focus on the two 7 keV SN models that bracket the range of matter power spectra compatible wit… ▽ More We use high-resolution hydrodynamical simulations run with the EAGLE model of galaxy formation to study the differences between the properties of - and subsequently the lensing signal from - subhaloes of massive elliptical galaxies at redshift 0.2, in Cold and Sterile Neutrino (SN) Dark matter models. We focus on the two 7 keV SN models that bracket the range of matter power spectra compatible with resonantly-produced SN as the source of the observed 3.5 keV line. We derive an accurate parametrisation for the subhalo mass function in these two SN models relative to CDM, as well as the subhalo spatial distribution, density profile, and projected number density and the dark matter fraction in subhaloes. We create mock lensing maps from the simulated haloes to study the differences in the lensing signal in the framework of subhalo detection. We find that subhalo convergence is well described by a log-normal distribution and that signal of subhaloes in the power spectrum is lower in SN models with respect to CDM, at a level of 10 to 80 per cent, depending on the scale. However, the scatter between different projections is large and might make the use of power-spectrum studies on the typical scales of current lensing images very difficult. Moreover, in the framework of individual detections through gravitational imaging a sample of ~30 lenses with an average sensitivity of M_sub=5x10^7M_sun would be required to discriminate between CDM and the considered sterile neutrino models. △ Less

Submitted 29 October, 2019; v1 submitted 15 July, 2019; originally announced July 2019.

Comments: 16 pages, 11 figures, accepted for publication in MNRAS

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

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

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

Submitted 4 July, 2019; originally announced July 2019.

Comments: 33 pages, 23 figures

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

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

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

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

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

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

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

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

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

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

The Detailed Science Case for the Maunakea Spectroscopic Explorer, 2019 edition

Authors: The MSE Science Team, Carine Babusiaux, Maria Bergemann, Adam Burgasser, Sara Ellison, Daryl Haggard, Daniel Huber, Manoj Kaplinghat, Ting Li, Jennifer Marshall, Sarah Martell, Alan McConnachie, Will Percival, Aaron Robotham, Yue Shen, Sivarani Thirupathi, Kim-Vy Tran, Christophe Yeche, David Yong, Vardan Adibekyan, Victor Silva Aguirre, George Angelou, Martin Asplund, Michael Balogh, Projjwal Banerjee , et al. (239 additional authors not shown)

Abstract: (Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the sc… ▽ More (Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more. △ Less

Submitted 9 April, 2019; originally announced April 2019.

Comments: 9 chapters, 301 pages, 100 figures. This version of the DSC is a comprehensive update of the original version, released in 2016, which can be downloaded at arXiv:1606.00043. A detailed summary of the design of MSE is available in the MSE Book 2018, available at arXiv:1810.08695

Astrophysical Tests of Dark Matter with Maunakea Spectroscopic Explorer

Authors: Ting S. Li, Manoj Kaplinghat, Keith Bechtol, Adam S. Bolton, Jo Bovy, Timothy Carleton, Chihway Chang, Alex Drlica-Wagner, Denis Erkal, Marla Geha, Johnny P. Greco, Carl J. Grillmair, Stacy Y. Kim, Chervin F. P. Laporte, Geraint F. Lewis, Martin Makler, Yao-Yuan Mao, Jennifer L. Marshall, Alan W. McConnachie, Lina Necib, A. M. Nierenberg, Brian Nord, Andrew B. Pace, Marcel S. Pawlowski, Annika H. G. Peter , et al. (5 additional authors not shown)

Abstract: We discuss how astrophysical observations with the Maunakea Spectroscopic Explorer (MSE), a high-multiplexity (about 4300 fibers), wide field-of-view (1.5 square degree), large telescope aperture (11.25 m) facility, can probe the particle nature of dark matter. MSE will conduct a suite of surveys that will provide critical input for determinations of the mass function, phase-space distribution, an… ▽ More We discuss how astrophysical observations with the Maunakea Spectroscopic Explorer (MSE), a high-multiplexity (about 4300 fibers), wide field-of-view (1.5 square degree), large telescope aperture (11.25 m) facility, can probe the particle nature of dark matter. MSE will conduct a suite of surveys that will provide critical input for determinations of the mass function, phase-space distribution, and internal density profiles of dark matter halos across all mass scales. N-body and hydrodynamical simulations of cold, warm, fuzzy and self-interacting dark matter suggest that non-trivial dynamics in the dark sector could have left an imprint on structure formation. Analysed within these frameworks, the extensive and unprecedented datasets produced by MSE will be used to search for deviations away from cold and collisionless dark matter model. MSE will provide an improved estimate of the local density of dark matter, critical for direct detection experiments, and will improve estimates of the J-factor for indirect searches through self-annihilation or decay into Standard Model particles. MSE will determine the impact of low mass substructures on the dynamics of Milky Way stellar streams in velocity space, and will allow for estimates of the density profiles of the dark matter halos of Milky Way dwarf galaxies using more than an order of magnitude more tracers. In the low redshift Universe, MSE will provide critical redshifts to pin down the luminosity functions of vast numbers of satellite systems, and MSE will be an essential component of future strong lensing measurements to constrain the halo mass function. Across nearly all mass scales, the improvements offered by MSE, in comparison to other facilities, are such that the relevant analyses are limited by systematics rather than statistics. △ Less

Submitted 9 April, 2019; v1 submitted 7 March, 2019; originally announced March 2019.

Comments: 44 pages, 19 figures. To appear as a chapter for "The Detailed Science Case for the Maunakea Spectroscopic Explorer, 2019"

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

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

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

Submitted 8 November, 2018; originally announced November 2018.

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

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

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

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

Submitted 8 November, 2018; originally announced November 2018.

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

The interplay of Self-Interacting Dark Matter and baryons in shaping the halo evolution

Authors: Giulia Despali, Martin Sparre, Simona Vegetti, Mark Vogelsberger, Jesús Zavala, Federico Marinacci

Abstract: We use high-resolution hydrodynamical simulation to test the difference of halo properties in cold dark matter (CDM) and a self-interacting dark matter (SIDM) scenario with a constant cross-section of $σ^\text{T}/m_χ=1\;\text{cm}^{2}\text{g}^{-1}$. We find that the interplay between dark matter self-interaction and baryonic physics induces a complex evolution of the halo properties, which depends… ▽ More We use high-resolution hydrodynamical simulation to test the difference of halo properties in cold dark matter (CDM) and a self-interacting dark matter (SIDM) scenario with a constant cross-section of $σ^\text{T}/m_χ=1\;\text{cm}^{2}\text{g}^{-1}$. We find that the interplay between dark matter self-interaction and baryonic physics induces a complex evolution of the halo properties, which depends on the halo mass and morphological type, as well as on the halo mass accretion history. While high mass haloes, selected as analogues of early-type galaxies, show cored profiles in the SIDM run, systems of intermediate mass and with a significant disk component can develop a profile that is similar or cuspier than in CDM. The final properties of SIDM haloes - measured at z=0.2 - correlate with the halo concentration and formation time, suggesting that the differences between different systems are due to the fact that we are observing the impact self-interaction. We also search for signatures of self-interacting dark matter in the lensing signal of the main haloes and hints of potential differences in the distribution of Einstein radii, which suggests that future wide-field survey might be able to distinguish between CDM and SIDM models on this basis. Finally, we find that the subhalo abundances are not altered in the adopted SIDM model with respect to CDM. △ Less

Submitted 28 January, 2019; v1 submitted 6 November, 2018; originally announced November 2018.

Comments: 11 pages, 11 figures, accepted for publication in MNRAS

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

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

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

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

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

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

A novel 3D technique to study the kinematics of lensed galaxies

Authors: Francesca Rizzo, Simona Vegetti, Filippo Fraternali, Enrico Di Teodoro

Abstract: We present a 3D Bayesian method to model the kinematics of strongly lensed galaxies from spatially-resolved emission-line observations. This technique enables us to simultaneously recover the lens-mass distribution and the source kinematics directly from the 3D data cube. We have tested this new method with simulated OSIRIS observations for nine star-forming lensed galaxies with different kinemati… ▽ More We present a 3D Bayesian method to model the kinematics of strongly lensed galaxies from spatially-resolved emission-line observations. This technique enables us to simultaneously recover the lens-mass distribution and the source kinematics directly from the 3D data cube. We have tested this new method with simulated OSIRIS observations for nine star-forming lensed galaxies with different kinematic properties. The simulated rotation curves span a range of shapes which are prototypes of different morphological galaxy types, from dwarf to massive spiral galaxies. We have found that the median relative accuracy on the inferred lens and kinematic parameters are at the level of 1 and 2 per cent, respectively. We have also tested the robustness of the technique against different inclination angles, signal-to-noise ratios, the presence of warps or non-circular motions and we have found that the accuracy stays within a few per cent in most cases. This technique represents a significant step forward with respect to the methods used until now, as the lens parameters and the kinematics of the source are derived from the same 3D data. This enables us to study the possible degeneracies between the two and estimate the uncertainties on all model parameters consistently. △ Less

Submitted 19 September, 2018; originally announced September 2018.

Comments: Accepted for publication in MNRAS

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

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

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

Submitted 15 July, 2018; originally announced July 2018.

Comments: 15 pages, 8 figures; published by MNRAS

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

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

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

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

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

Probing sub-galactic mass structure with the power spectrum of surface-brightness anomalies in high-resolution observations of galaxy-galaxy strong gravitational lenses. II.Observational constraints on the sub-galactic matter power spectrum

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

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

Submitted 1 February, 2023; v1 submitted 15 March, 2018; originally announced March 2018.

Comments: 16 pages, 10 figures, updated and resubmitted to MNRAS together with an accompanying methodology paper in response to reviewer's comments

SDSS-IV MaNGA: The Spectroscopic Discovery of Strongly Lensed Galaxies

Authors: Michael S. Talbot, Joel R. Brownstein, Adam S. Bolton, Kevin Bundy, Brett H. Andrews, Brian Cherinka, Thomas E. Collett, Anupreeta More, Surhud More, Alessandro Sonnenfeld, Simona Vegetti, David A. Wake, Anne-Marie Weijmans, Kyle B. Westfall

Abstract: We present a catalogue of 38 spectroscopically detected strong galaxy-galaxy gravitational lens candidates identified in the Sloan Digital Sky Survey IV (SDSS-IV). We were able to simulate narrow-band images for 8 of them demonstrating evidence of multiple images. Two of our systems are compound lens candidates, each with 2 background source-planes. One of these compound systems shows clear lensin… ▽ More We present a catalogue of 38 spectroscopically detected strong galaxy-galaxy gravitational lens candidates identified in the Sloan Digital Sky Survey IV (SDSS-IV). We were able to simulate narrow-band images for 8 of them demonstrating evidence of multiple images. Two of our systems are compound lens candidates, each with 2 background source-planes. One of these compound systems shows clear lensing features in the narrow-band image. Our sample is based on 2812 galaxies observed by the Mapping Nearby Galaxies at APO (MaNGA) integral field unit (IFU). This Spectroscopic Identification of Lensing Objects (SILO) survey extends the methodology of the Sloan Lens ACS Survey (SLACS) and BOSS Emission-Line Survey (BELLS) to lower redshift and multiple IFU spectra. We searched ~ 1.5 million spectra, of which 3065 contained multiple high signal-to-noise background emission-lines or a resolved [OII] doublet, that are included in this catalogue. Upon manual inspection, we discovered regions with multiple spectra containing background emission-lines at the same redshift, providing evidence of a common source-plane geometry which was not possible in previous SLACS and BELLS discovery programs. We estimate more than half of our candidates have an Einstein radius > 1.7", which is significantly greater than seen in SLACS and BELLS. These larger Einstein radii produce more extended images of the background galaxy increasing the probability that a background emission-line will enter one of the IFU spectroscopic fibres, making detection more likely. △ Less

Submitted 9 March, 2018; originally announced March 2018.

Comments: Accepted for publication in MNRAS, March 8, 2018. In press. 16 pages, 5 figures, 4 tables

Extreme conditions in the molecular gas of lensed star-forming galaxies at z~3

Authors: Paola Andreani, Edwin Retana-Montenegro, Zhi-Yu Zhang, Padelis Papadopoulos, Chentao Yang, Simona Vegetti

Abstract: Atomic Carbon can be an efficient tracer of the molecular gas mass, and when combined to the detection of high-J and low-J CO lines it yields also a sensitive probe of the power sources in the molecular gas of high redshift galaxies. The recently installed SEPIA5 receiver at the focus of the APEX telescope has opened up a new window at frequencies 159 - 211 GHz allowing the exploration of the Atom… ▽ More Atomic Carbon can be an efficient tracer of the molecular gas mass, and when combined to the detection of high-J and low-J CO lines it yields also a sensitive probe of the power sources in the molecular gas of high redshift galaxies. The recently installed SEPIA5 receiver at the focus of the APEX telescope has opened up a new window at frequencies 159 - 211 GHz allowing the exploration of the Atomic Carbon in high-z galaxies, at previously inaccessible frequencies from the ground. We have targeted three gravitationally lensed galaxies at redshift of about 3 and conducted a comparative study of the observed high-J CO/CI ~ratios with well-studied nearby galaxies. Atomic Carbon (CI(2-1)) was detected in one of the three targets and marginally in a second, while in all three targets the $J=7\to6$ CO line is detected. The CO(7-6)/CI(2-1), CO(7-6)/CO(1-0) line ratios and the CO(7-6)/(far-IR continuum) luminosity ratio are compared to those of nearby objects. A large excitation status in the ISM of these high-z objects is seen, unless differential lensing unevenly boosts the CO line fluxes from the warm and dense gas more than the CO(1-0), CI(2-1), tracing a more widely distributed cold gas phase. We provide estimates of total molecular gas masses derived from the atomic Carbon and the Carbon monoxide CO(1-0), which within the uncertainties turn out to be equal. △ Less

Submitted 10 April, 2018; v1 submitted 21 February, 2018; originally announced February 2018.

Comments: A&A, in press

Journal ref: A&A 615, A142 (2018)

Time Delay Lens Modeling Challenge: I. Experimental Design

Authors: Xuheng Ding, Tommaso Treu, Anowar J. Shajib, Dandan Xu, Geoff C. -F. Chen, Anupreeta More, Giulia Despali, Matteo Frigo, Christopher D. Fassnacht, Daniel Gilman, Stefan Hilbert, Philip J. Marshall, Dominique Sluse, Simona Vegetti

Abstract: Strong gravitational lenses with measured time delay are a powerful tool to measure cosmological parameters, especially the Hubble constant ($H_0$). Recent studies show that by combining just three multiply-imaged AGN systems, one can determine $H_0$ to 2.4% precision. Furthermore, the number of time-delay lens systems is growing rapidly, enabling the determination of $H_0$ to 1% precision in the… ▽ More Strong gravitational lenses with measured time delay are a powerful tool to measure cosmological parameters, especially the Hubble constant ($H_0$). Recent studies show that by combining just three multiply-imaged AGN systems, one can determine $H_0$ to 2.4% precision. Furthermore, the number of time-delay lens systems is growing rapidly, enabling the determination of $H_0$ to 1% precision in the near future. However, as the precision increases it is important to ensure that systematic errors and biases remain subdominant. For this purpose, challenges with simulated datasets are a key component in this process. Following the experience of the past challenge on time delay, where it was shown that time delays can indeed be measured precisely and accurately at the sub-percent level, we now present the "Time Delay Lens Modeling Challenge" (TDLMC). The goal of this challenge is to assess the present capabilities of lens modeling codes and assumptions and test the level of accuracy of inferred cosmological parameters given realistic mock datasets. We invite scientists to model a set of simulated HST observations of 50 mock lens systems. The systems are organized in rungs, with the complexity and realism increasing going up the ladder. The goal of the challenge is to infer $H_0$ for each rung, given the HST images, the time delay, and stellar velocity dispersion of the deflector for a fixed background cosmology. The TDLMC challenge starts with the mock data release on 2018 January 8th. The deadline for blind submission is different for each rung. The deadline for Rung 0-1 is 2018 September 8; the deadline for Rung 2 is 2019 April 8 and the one for Rung 3 is 2019 September 8. This first paper gives an overview of the challenge including the data design, and a set of metrics to quantify the modeling performance and challenge details. △ Less

Submitted 15 June, 2020; v1 submitted 4 January, 2018; originally announced January 2018.

Comments: 6 pages, 2 figures, submitted to MNRAS

Constraining sterile neutrino cosmologies with strong gravitational lensing observations at redshift z~0.2

Authors: S. Vegetti, G. Despali, M. R. Lovell, W. Enzi

Abstract: We use the observed amount of subhaloes and line-of-sight dark matter haloes in a sample of 11 gravitational lens systems from the Sloan Lens ACS Survey to constrain the free-streaming properties of the dark matter particles. In particular, we combine the detection of a small-mass dark matter halo by Vegetti et al. 2010 with the non-detections by Vegetti et al. 2014 and compare the derived subhalo… ▽ More We use the observed amount of subhaloes and line-of-sight dark matter haloes in a sample of 11 gravitational lens systems from the Sloan Lens ACS Survey to constrain the free-streaming properties of the dark matter particles. In particular, we combine the detection of a small-mass dark matter halo by Vegetti et al. 2010 with the non-detections by Vegetti et al. 2014 and compare the derived subhalo and halo mass functions with expectations from cold dark matter (CDM) and resonantly produced sterile neutrino models. We constrain the half-mode mass, i.e. the mass scale at which the linear matter power spectrum is reduced by 50 per cent relatively to the CDM model, to be $\log M_{\rm{hm}} \left[M_\odot\right] < 12.0$ (equivalent thermal relic mass $m_{\rm th} > 0.3$ keV) at the 2$σ$ level. This excludes sterile neutrino models with neutrino masses $m_{\rm s} < 0.8$ keV at any value of $L_{\rm 6}$. Our constraints are weaker than currently provided by the number of Milky Way satellites, observations of the 3.5 keV X-ray line, and the Lyman $α$ forest. However, they are more robust than the former as they are less affected by baryonic processes. Moreover, unlike the latter, they are not affected by assumptions on the thermal histories for the intergalactic medium. Gravitational lens systems with higher data quality and higher source and lens redshift are required to obtain tighter constraints. △ Less

Submitted 9 October, 2018; v1 submitted 4 January, 2018; originally announced January 2018.

Comments: Accepted for publication on MNRAS

Modelling the line-of-sight contribution in substructure lensing

Authors: Giulia Despali, Simona Vegetti, Simon D. M. White, Carlo Giocoli, Frank C. van den Bosch

Abstract: We investigate how Einstein rings and magnified arcs are affected by small-mass dark-matter haloes placed along the line-of-sight to gravitational lens systems. By comparing the gravitational signature of line-of-sight haloes with that of substructures within the lensing galaxy, we derive a mass-redshift relation that allows us to rescale the detection threshold (i.e. lowest detectable mass) for s… ▽ More We investigate how Einstein rings and magnified arcs are affected by small-mass dark-matter haloes placed along the line-of-sight to gravitational lens systems. By comparing the gravitational signature of line-of-sight haloes with that of substructures within the lensing galaxy, we derive a mass-redshift relation that allows us to rescale the detection threshold (i.e. lowest detectable mass) for substructures to a detection threshold for line-of-sight haloes at any redshift. We then quantify the line-of-sight contribution to the total number density of low-mass objects that can be detected through strong gravitational lensing. Finally, we assess the degeneracy between substructures and line-of-sight haloes of different mass and redshift to provide a statistical interpretation of current and future detections, with the aim of distinguishing between CDM and WDM. We find that line-of-sight haloes statistically dominate with respect to substructures, by an amount that strongly depends on the source and lens redshifts, and on the chosen dark matter model. Substructures represent about 30 percent of the total number of perturbers for low lens and source redshifts (as for the SLACS lenses), but less than 10 per cent for high redshift systems. We also find that for data with high enough signal-to-noise ratio and angular resolution, the non-linear effects arising from a double-lens-plane configuration are such that one is able to observationally recover the line-of-sight halo redshift with an absolute error precision of 0.15 at the 68 per cent confidence level. △ Less

Submitted 26 January, 2018; v1 submitted 13 October, 2017; originally announced October 2017.

Comments: 19 pages, 17 figures, accepted for publication in MNRAS

Flux-ratio anomalies from discs and other baryonic structures in the Illustris simulation

Authors: Jen-Wei Hsueh, Giulia Despali, Simona Vegetti, Dandan Xu, Christopher D. Fassnacht, R. Benton Metcalf

Abstract: The flux ratios in the multiple images of gravitationally lensed quasars can provide evidence for dark matter substructure in the halo of the lensing galaxy if the flux ratios differ from those predicted by a smooth model of the lensing galaxy mass distribution. However, it is also possible that baryonic structures in the lensing galaxy, such as edge-on discs, can produce flux-ratio anomalies. In… ▽ More The flux ratios in the multiple images of gravitationally lensed quasars can provide evidence for dark matter substructure in the halo of the lensing galaxy if the flux ratios differ from those predicted by a smooth model of the lensing galaxy mass distribution. However, it is also possible that baryonic structures in the lensing galaxy, such as edge-on discs, can produce flux-ratio anomalies. In this work, we present the first statistical analysis of flux-ratio anomalies due to baryons from a numerical simulation perspective. We select galaxies with various morphological types in the Illustris simulation and ray-trace through the simulated halos, which include baryons in the main lensing galaxies but exclude any substructures, in order to explore the pure baryonic effects. Our ray-tracing results show that the baryonic components can be a major contribution to the flux-ratio anomalies in lensed quasars and that edge-on disc lenses induce the strongest anomalies. We find that the baryonic components increase the probability of finding high flux-ratio anomalies in the early-type lenses by about 8% and by about 10 - 20% in the disc lenses. The baryonic effects also induce astrometric anomalies in 13% of the mock lenses. Our results indicate that the morphology of the lens galaxy becomes important in the analysis of flux-ratio anomalies when considering the effect of baryons, and that the presence of baryons may also partially explain the discrepancy between the observed (high) anomaly frequency and what is expected due to the presence of subhalos as predicted by the CDM simulations. △ Less

Submitted 10 January, 2018; v1 submitted 24 July, 2017; originally announced July 2017.

Comments: 16 pages, 11 figures, accepted by MNRAS

ALMA view of RX J1131-1231: Sub-kpc CO (2-1) mapping of a molecular disk in a lensed star-forming quasar host galaxy

Authors: D. Paraficz, M. Rybak, J. P. McKean, S. Vegetti, D. Sluse, F. Courbin, H. R. Stacey, S. H. Suyu, M. Dessauges-Zavadsky, C. D. Fassnacht, L. V. E. Koopmans

Abstract: We present ALMA 2-mm continuum and CO (2-1) spectral line imaging of the gravitationally lensed z=0.654 star-forming/quasar composite RX J1131-1231 at 240-400 mas angular resolution. The continuum emission is found to be compact and coincident with the optical emission, whereas the molecular gas forms a complete Einstein ring, which shows strong differential magnification. The de-lensed source str… ▽ More We present ALMA 2-mm continuum and CO (2-1) spectral line imaging of the gravitationally lensed z=0.654 star-forming/quasar composite RX J1131-1231 at 240-400 mas angular resolution. The continuum emission is found to be compact and coincident with the optical emission, whereas the molecular gas forms a complete Einstein ring, which shows strong differential magnification. The de-lensed source structure is determined on 400-pc resolution using a visibility-fitting lens modelling technique. The reconstructed molecular gas velocity-field is consistent with a rotating disk with a maximum rotational velocity of 280 km/s. From dynamical model fitting we find an enclosed mass M(r<5 kpc)=(1.46+/-0.31)*10^11 M_sol. The molecular gas distribution is highly structured, with clumps that are co-incident with higher gas velocity dispersion regions 40-50 km/s and with the intensity peaks in the optical emission, which are associated with sites of on-going turbulent star-formation. The peak in the CO (2-1) distribution is not co-incident with the AGN, where there is a paucity of molecular gas emission, possibly due to radiative feedback from the central engine. The intrinsic molecular gas luminosity is L'_CO=(1.2+/-0.3)*10^10 K km/s pc^2 and the inferred gas mass is M(H2)=(8.3+/-3.0)*10^10 M_sol, which given its dynamical mass is consistent with a CO-H2 conversion factor of alpha = 5.5+/-2.0 M_solar(K km/s pc^2)^-1. This suggests that the star-formation efficiency is dependent on the host galaxy morphology as opposed to the nature of the AGN. The far-infrared continuum spectral energy distribution shows evidence for heated dust, equivalent to an obscured star-formation rate of SFR=69^(+41)_(-25)*(7.3/u_IR)M_sol/yr, which demonstrates the composite star-forming/AGN nature of this system. RX J1131-1231 △ Less

Submitted 20 December, 2017; v1 submitted 28 May, 2017; originally announced May 2017.

Comments: 12 pages, 8 figures, accepted by A&A

Journal ref: A&A 613, A34 (2018)

SHARP - IV. An apparent flux ratio anomaly resolved by the edge-on disc in B0712+472

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

Abstract: Flux ratio anomalies in quasar lenses can be attributed to dark matter substructure surrounding the lensing galaxy and, thus, used to constrain the substructure mass fraction. Previous applications of this approach infer a substructure abundance that potentially in tension with the predictions of a $Λ$CDM cosmology. However, the assumption that all flux ratio anomalies are due to substructure is a… ▽ More Flux ratio anomalies in quasar lenses can be attributed to dark matter substructure surrounding the lensing galaxy and, thus, used to constrain the substructure mass fraction. Previous applications of this approach infer a substructure abundance that potentially in tension with the predictions of a $Λ$CDM cosmology. However, the assumption that all flux ratio anomalies are due to substructure is a strong one, and alternative explanations have not been fully investigated. Here, we use new high-resolution near-IR Keck~II adaptive optics imaging for the lens system CLASS B0712+472 to perform pixel-based lens modelling for this system and, in combination with new VLBA radio observations, show that the inclusion of the disc in the lens model can explain the flux ratio anomalies without the need for dark matter substructures. The projected disc mass comprises 16% of the total lensing mass within the Einstein radius and the total disc mass is $1.79 \times 10^{10} M_{sun}$. The case of B0712+472 adds to the evidence that not all flux ratio anomalies are due to dark subhaloes, and highlights the importance of taking the effects of baryonic structures more fully into account in order to obtain an accurate measure of the substructure mass fraction. △ Less

Submitted 27 April, 2017; v1 submitted 23 January, 2017; originally announced January 2017.

Comments: 9 pages, 5 figures, 5 tables. Accepted by MNRAS