sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Birrer, Simon

Normalized to: Birrer, S.

39 article(s) in total. 407 co-authors, from 1 to 24 common article(s). Median position in authors list is 4,0.

[1]  [pdf] - 2072774
STRIDES: Spectroscopic and photometric characterization of the environment and effects of mass along the line of sight to the gravitational lenses DES J0408-5354 and WGD 2038-4008
Comments: 40 pages, 18 figures Added missing author to the author list
Submitted: 2020-03-26, last modified: 2020-03-30
In time-delay cosmography, three of the key ingredients are 1) determining the velocity dispersion of the lensing galaxy, 2) identifying galaxies and groups along the line of sight with sufficient proximity and mass to be included in the mass model, and 3) estimating the external convergence $\kappa_\mathrm{ext}$ from less massive structures that are not included in the mass model. We present results on all three of these ingredients for two time-delay lensed quasar systems, DES J0408-5354 and WGD 2038-4008. We use the Gemini, Magellan and VLT telescopes to obtain spectra to both measure the stellar velocity dispersions of the main lensing galaxies and to identify the line-of-sight galaxies in these systems. Next, we identify 10 groups in DES J0408-5354 and 2 groups in WGD 2038-4008using a group-finding algorithm. We then identify the most significant galaxy and galaxy-group perturbers using the "flexion shift" criterion. We determine the probability distribution function of the external convergence $\kappa_\mathrm{ext}$ for both of these systems based on our spectroscopy and on the DES-only multiband wide-field observations. Using weighted galaxy counts, calibrated based on the Millennium Simulation, we find that DES J0408-5354 is located in a significantly underdense environment, leading to a tight (width $\sim3\%$), negative-value $\kappa_\mathrm{ext}$ distribution. On the other hand, WGD 2038-4008 is located in an environment of close to unit density, and its low source redshift results in a much tighter $\kappa_\mathrm{ext}$ of $\sim1\%$, as long as no external shear constraints are imposed.
[2]  [pdf] - 2060241
The first simultaneous measurement of Hubble constant and post-Newtonian parameter from Time-Delay Strong Lensing
Comments: 6 pages, 2 figures, comments welcome
Submitted: 2020-03-06
Strong gravitational lensing has been a powerful probe of cosmological models and gravity. To date, constraints in either domain have been separately. We propose a new methodology through which the cosmological model, specifically the Hubble constant, and post-Newtonian parameter can be simultaneously constrained. Using the time-delay cosmography from strong lensing combined with the stellar kinematics of the deflector lens, we demonstrate the Hubble constant and post-Newtonian parameter are incorporated in two distance ratios which reflect the lensing mass and dynamical mass, respectively. Through the reanalysis of the four publicly released lenses distance posteriors from the H0LiCOW collaboration, the simultaneous constraints of Hubble constant and post-Newtonian parameter are obtained. Our results suggests no deviation from the General Relativity, $\gamma_{\texttt{PPN}}=0.87^{+0.19}_{-0.17}$ with a Hubble constant favors the local universe value, $H_0=73.65^{+1.95}_{-2.26}$ km s$^{-1}$ Mpc$^{-1}$. Finally, we forecast the robustness of gravity tests by using the time-delay strong lensing for constraints we expect in the next few years. We find that the joint constraint from 40 lenses are able to reach the order of $7.7\%$ for the post-Newtonian parameter and $1.4\%$ for Hubble parameter.
[3]  [pdf] - 2035848
A versatile tool for cluster lensing source reconstruction. I. methodology and illustration on sources in the Hubble Frontier Field Cluster MACS J0717.5+3745
Comments: Submitted to MNRAS, 16 pages, 11 figures, 6 tables
Submitted: 2020-01-21
We describe a general purpose method to reconstruct the intrinsic properties of sources lensed by the gravitational potential of foreground clusters of galaxies. The tool Lenstruction is implemented in the publicly available multi-purpose gravitational lensing software Lenstronomy, in order to provide easy and fast solution to this common astrophysical problem. The tool is based on forward modeling the appearance of the source in the image plane, taking into account the distortion by lensing and the instrumental point spread function (PSF). For singly-imaged sources a global lens model in the format of the Hubble Frontier Fields (HFF) lensing maps is required as a starting point. For multiply-imaged sources, the tool can also fit and apply first (deflection) second (shear, convergence) and third order (flexion) corrections to the local gravitational potential to improve the overall reconstruction, depending on the quality of the data. We illustrate the performance and features of the code with two examples of multiply-imaged systems taken from the Hubble Frontier Fields, starting from five different publicly available cluster models. We find that, after our correction, the relative magnification - and other lensing properties - between the multiple images become robustly constrained. Furthermore, we find that scatter between models of the reconstructed source size and magnitude is reduced. The code and jupyter notebooks are publicly available.
[4]  [pdf] - 2065348
The STRong lensing Insights into the Dark Energy Survey (STRIDES) 2017/2018 follow-up campaign: Discovery of 10 lensed quasars and 10 quasar pairs
Comments: 22 pages, 16 figures, submitted to MNRAS
Submitted: 2019-12-19
We report the results of the STRong lensing Insights from the Dark Energy Survey (STRIDES) follow-up campaign of the late 2017/early 2018 season. We obtained spectra of 65 lensed quasar candidates either with EFOSC2 on the NTT or ESI on Keck, which confirm 10 new gravitationally lensed quasars and 10 quasar pairs with similar spectra, but which do not show a lensing galaxy in DES images. Eight lensed quasars are doubly imaged with source redshifts between 0.99 and 2.90, one is triply imaged by a group (DESJ0345-2545, $z=1.68$), and one is quadruply imaged (quad: DESJ0053-2012, $z=3.8$). Singular isothermal ellipsoid models for the doubles, based on high-resolution imaging from SAMI on SOAR or NIRC2 on Keck, give total magnifications between 3.2 and 5.6, and Einstein radii between 0.49 and 1.97 arcseconds. After spectroscopic follow-up, we extract multi-epoch $grizY$ photometry of confirmed lensed quasars and contaminant quasar+star pairs from the first 4 years of DES data using parametric multi-band modelling, and compare variability in each system's components. By measuring the reduced ${\chi}^2$ associated with fitting all epochs to the same magnitude, we find a simple cut on the less variable component that retains all confirmed lensed quasars, while removing 94 per cent of contaminant systems with stellar components. Based on our spectroscopic follow-up, this variability information can improve selection of lensed quasars and quasar pairs from 34-45 per cent to 51-70 per cent, with the majority of remaining contaminants being compact star-forming galaxies. Using mock lensed quasar lightcurves we demonstrate that selection based only on variability will over-represent the quad fraction by 10 per cent over a complete DES magnitude-limited sample (excluding microlensing differences), explained by the magnification bias and hence lower luminosity (more variable) sources in quads.
[5]  [pdf] - 2015915
TDCOSMO. I. An exploration of systematic uncertainties in the inference of $H_0$ from time-delay cosmography
Comments: 18 pages, 8 figures, 3 tables
Submitted: 2019-12-17
Time-delay cosmography of lensed quasars has achieved 2.4% precision on the measurement of the Hubble Constant, $H_0$. As part of an ongoing effort to uncover and control systematic uncertainties, we investigate three potential sources: 1- stellar kinematics, 2- line-of-sight effects, 3- deflector mass model. To meet this goal in a quantitative way, we mimic closely the H0LiCOW/SHARP/STRIDES procedures (i.e., TDCOSMO), and we find the following. First, stellar kinematics cannot be a dominant source of error or bias given current uncertainties. Second, we find no bias arising from incorrect estimation of the line-of-sight effects. Third, we show that elliptical composite (stars + dark matter halo), power-law, and cored power-law mass profiles have the flexibility to yield a broad range in $H_0$ values. However, the TDCOSMO procedures to model the data with both composite and power-law mass profiles are informative. If the models agree, as we observe in real systems owing to the "bulge-halo" conspiracy, $H_0$ is recovered precisely by both models. If the two models disagreed, as in the case of some pathological models illustrated here, the TDCOSMO procedure would either be able to discriminate between them through the goodness of fit, or account for the discrepancy in the final error bars provided by the analysis. This conclusion is consistent with a reanalysis of the TDCOSMO (real) lenses: the composite model yields $H_0=74.2^{+1.6}_{-1.6}$ ${\, \mathrm{km}\, \mathrm{s}^{-1}\, \mathrm{Mpc}^{-1}}$, while the power-law model yields $74.0^{+1.7}_{-1.8}$ ${\, \mathrm{km}\, \mathrm{s}^{-1}\, \mathrm{Mpc}^{-1}}$. In conclusion, we find no evidence of bias or errors larger than the current statistical uncertainties reported by TDCOSMO.
[6]  [pdf] - 2025707
Double dark matter vision: twice the number of compact-source lenses with narrow-line lensing and the WFC3 grism
Comments: Accepted for publication in MNRAS
Submitted: 2019-08-17, last modified: 2019-12-16
The magnifications of compact-source lenses are extremely sensitive to the presence of low mass dark matter halos along the entire sight line from the source to the observer. Traditionally, the study of dark matter structure in compact-source strong gravitational lenses has been limited to radio-loud systems, as the radio emission is extended and thus unaffected by microlensing which can mimic the signal of dark matter structure. An alternate approach is to measure quasar nuclear-narrow line emission, which is free from microlensing and present in virtually all quasar lenses. In this paper, we double the number of systems which can be used for gravitational lensing analyses by presenting measurements of narrow-line emission from a sample of 8 quadruply imaged quasar lens systems, WGD J0405-3308, HS 0810+2554, RX J0911+0551, SDSS J1330+1810, PS J1606-2333, WFI 2026-4536, WFI 2033-4723 and WGD J2038-4008. We describe our updated grism spectral modelling pipeline, which we use to measure narrow-line fluxes with uncertainties of 2-10\%, presented here. We fit the lensed image positions with smooth mass models and demonstrate that these models fail to produce the observed distribution of image fluxes over the entire sample of lenses. Furthermore, typical deviations are larger than those expected from macromodel uncertainties. This discrepancy indicates the presence of perturbations caused by small-scale dark matter structure. The interpretation of this result in terms of dark matter models is presented in a companion paper.
[7]  [pdf] - 2011870
H0LiCOW XII. Lens mass model of WFI2033-4723 and blind measurement of its time-delay distance and $H_0$
Comments: Version accepted by MNRAS. 29 pages including appendix, 17 figures, 6 tables. arXiv admin note: text overlap with arXiv:1607.01403
Submitted: 2019-05-22, last modified: 2019-12-08
We present the lens mass model of the quadruply-imaged gravitationally lensed quasar WFI2033-4723, and perform a blind cosmographical analysis based on this system. Our analysis combines (1) time-delay measurements from 14 years of data obtained by the COSmological MOnitoring of GRAvItational Lenses (COSMOGRAIL) collaboration, (2) high-resolution $\textit{Hubble Space Telescope}$ imaging, (3) a measurement of the velocity dispersion of the lens galaxy based on ESO-MUSE data, and (4) multi-band, wide-field imaging and spectroscopy characterizing the lens environment. We account for all known sources of systematics, including the influence of nearby perturbers and complex line-of-sight structure, as well as the parametrization of the light and mass profiles of the lensing galaxy. After unblinding, we determine the effective time-delay distance to be $4784_{-248}^{+399}~\mathrm{Mpc}$, an average precision of $6.6\%$. This translates to a Hubble constant $H_{0} = 71.6_{-4.9}^{+3.8}~\mathrm{km~s^{-1}~Mpc^{-1}}$, assuming a flat $\Lambda$CDM cosmology with a uniform prior on $\Omega_\mathrm{m}$ in the range [0.05, 0.5]. This work is part of the $H_0$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration, and the full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper (H0LiCOW XIII).
[8]  [pdf] - 2044602
The mass relations between supermassive black holes and their host galaxies at 1<z<2 with HST-WFC3
Comments: Accepted for publication in ApJ; 28 pages, 14 figures, 7 tables
Submitted: 2019-10-25, last modified: 2019-11-20
Correlations between the mass of a supermassive black hole and the properties of its host galaxy (e.g., total stellar mass (M*), luminosity (Lhost)) suggest an evolutionary connection. A powerful test of a co-evolution scenario is to measure the relations MBH-Lhost and MBH-M* at high redshift and compare with local estimates. For this purpose, we acquired HST imaging with WFC3 of 32 X-ray-selected broad-line AGN at 1.2<z<1.7 in deep survey fields. By applying state-of-the-art tools to decompose the HST images including available ACS data, we measured the host galaxy luminosity and stellar mass along with other properties through the 2D model fitting. The black hole mass was determined using the broad Halpha line, detected in the near-infrared with Subaru/FMOS, which potentially minimizes systematic effects using other indicators. We find that the observed ratio of MBH to total M* is 2.7 times larger at z~1.5 than in the local universe, while the scatter is equivalent between the two epochs. A non-evolving mass ratio is consistent with the data at the 2-3 sigma confidence level when accounting for selection effects and their uncertainties. The relationship between MBH-Lhost paints a similar picture. Therefore, our results cannot distinguish whether SMBHs and their total M* and Lhost proceed in lockstep or whether the growth of the former somewhat overshoots the latter, given the uncertainties. Based on a statistical estimate of the bulge-to-total mass fraction, the ratio MBH/M* is offset from the local value by a factor of ~7 which is significant even accounting for selection effects. Taken together, these observations are consistent with a scenario in which stellar mass is subsequently transferred from an angular momentum supported component of the galaxy to the pressure supported one through secular processes or minor mergers at a faster rate than mass accretion onto the SMBH.
[9]  [pdf] - 2025711
Warm dark matter chills out: constraints on the halo mass function and the free-streaming length of dark matter with 8 quadruple-image strong gravitational lenses
Submitted: 2019-08-19, last modified: 2019-11-11
The free-streaming length of dark matter depends on fundamental dark matter physics, and determines the abundance and concentration of dark matter halos on sub-galactic scales. Using the image positions and flux ratios from eight quadruply-imaged quasars, we constrain the free-streaming length of dark matter and the amplitude of the subhalo mass function (SHMF). We model both main deflector subhalos and halos along the line of sight, and account for warm dark matter (WDM) free-streaming effects on the mass function and mass-concentration relation. By calibrating the scaling of the SHMF with host halo mass and redshift using a suite of simulated halos, we infer a global normalization for the SHMF. We account for finite-size background sources, and marginalize over the mass profile of the main deflector. Parameterizing dark matter free-streaming through the half-mode mass $m_{\rm{hm}}$, we constrain the thermal relic particle mass $m_{\rm{DM}}$ corresponding to $m_{\rm{hm}}$. At $95 \%$ CI: $m_{\rm{hm}} < 10^{7.8} M_{\odot}$ ($m_{\rm{DM}} > 5.2 \ \rm{keV}$). We disfavor $m_{\rm{DM}} = 4.0 \rm{keV}$ and $ m_{\rm{DM}} = 3.0 \rm{keV}$ with likelihood ratios of 7:1 and 30:1, respectively, relative to the peak of the posterior distribution. Assuming cold dark matter, we constrain the projected mass in substructure between $10^6 - 10^{9} M_{\odot}$ near lensed images. At $68 \%$ CI, we infer $2.0 - 6.1 \times 10^{7} M_{\odot} \rm{kpc^{-2}}$, corresponding to mean projected mass fraction $\bar{f}_{\rm{sub}} = 0.035_{-0.017}^{+0.021}$. At $95 \%$ CI, we obtain a lower bound on the projected mass of $0.6 \times 10^{7} M_{\odot} \rm{kpc^{-2}}$, corresponding to $\bar{f}_{\rm{sub}} > 0.005$. These results agree with the predictions of cold dark matter.
[10]  [pdf] - 1994089
H0LiCOW XIII. A 2.4% measurement of $H_{0}$ from lensed quasars: $5.3\sigma$ tension between early and late-Universe probes
Comments: Accepted for publication in MNRAS; 23 pages, 13 figures, 8 tables
Submitted: 2019-07-10, last modified: 2019-11-05
We present a measurement of the Hubble constant ($H_{0}$) and other cosmological parameters from a joint analysis of six gravitationally lensed quasars with measured time delays. All lenses except the first are analyzed blindly with respect to the cosmological parameters. In a flat $\Lambda$CDM cosmology, we find $H_{0} = 73.3_{-1.8}^{+1.7}$, a 2.4% precision measurement, in agreement with local measurements of $H_{0}$ from type Ia supernovae calibrated by the distance ladder, but in $3.1\sigma$ tension with $Planck$ observations of the cosmic microwave background (CMB). This method is completely independent of both the supernovae and CMB analyses. A combination of time-delay cosmography and the distance ladder results is in $5.3\sigma$ tension with $Planck$ CMB determinations of $H_{0}$ in flat $\Lambda$CDM. We compute Bayes factors to verify that all lenses give statistically consistent results, showing that we are not underestimating our uncertainties and are able to control our systematics. We explore extensions to flat $\Lambda$CDM using constraints from time-delay cosmography alone, as well as combinations with other cosmological probes, including CMB observations from $Planck$, baryon acoustic oscillations, and type Ia supernovae. Time-delay cosmography improves the precision of the other probes, demonstrating the strong complementarity. Allowing for spatial curvature does not resolve the tension with $Planck$. Using the distance constraints from time-delay cosmography to anchor the type Ia supernova distance scale, we reduce the sensitivity of our $H_0$ inference to cosmological model assumptions. For six different cosmological models, our combined inference on $H_{0}$ ranges from $\sim73$-$78~\mathrm{km~s^{-1}~Mpc^{-1}}$, which is consistent with the local distance ladder constraints.
[11]  [pdf] - 1994455
Algorithms and Statistical Models for Scientific Discovery in the Petabyte Era
Comments: arXiv admin note: substantial text overlap with arXiv:1905.05116
Submitted: 2019-11-04
The field of astronomy has arrived at a turning point in terms of size and complexity of both datasets and scientific collaboration. Commensurately, algorithms and statistical models have begun to adapt --- e.g., via the onset of artificial intelligence --- which itself presents new challenges and opportunities for growth. This white paper aims to offer guidance and ideas for how we can evolve our technical and collaborative frameworks to promote efficient algorithmic development and take advantage of opportunities for scientific discovery in the petabyte era. We discuss challenges for discovery in large and complex data sets; challenges and requirements for the next stage of development of statistical methodologies and algorithmic tool sets; how we might change our paradigms of collaboration and education; and the ethical implications of scientists' contributions to widely applicable algorithms and computational modeling. We start with six distinct recommendations that are supported by the commentary following them. This white paper is related to a larger corpus of effort that has taken place within and around the Petabytes to Science Workshops (
[12]  [pdf] - 2026207
Where do quasar hosts lie with respect to the size-mass relation of galaxies?
Comments: 4 pages; 1 figure; To be resubmitted to ApJ Letters imminently; minor corrections based on comments from the referee have been included; Further comments are welcome
Submitted: 2019-10-30
The evolution of the galaxy size - stellar mass (Mstellar) relation has been a puzzle for over a decade. High redshift galaxies are significantly more compact than galaxies observed today, at an equivalent mass, but how much of this apparent growth is driven by progenitor bias, minor mergers, secular processes, or feedback from AGN is unclear. To help disentangle the physical mechanisms at work by addressing the latter, we study the galaxy size - Mstellar relation of 32 carefully-selected broad-line AGN hosts at 1.2 < z < 1.7 (7.5 < log M_BH < 8.5; L_bol/L_Edd > 0.1). Using HST with multi-band photometry and state-of-the-art modeling techniques, we measure half-light radii while accounting for uncertainties from subtracting bright central point sources. We find AGN hosts to have sizes ranging from 1 to 6 kpc at Mstellar ~ 0.3 - 1 x 10^11 Msun. Thus, many hosts have intermediate sizes as compared to equal-mass star-forming and quiescent galaxies. While inconsistent with the idea that AGN feedback may induce an increase in galaxy sizes, this finding is consistent with hypotheses in which AGNs preferentially occur in systems with prior concentrated gas reservoirs, or are involved in secular compaction processes perhaps responsible for simultaneously building bulges and shutting down star formation. If driven by minor mergers, which do not grow central black holes as fast as they do bulge-like stellar structures, such a process would explain both the galaxy size - mass relation observed here and the evolution in the black hole, bulge mass relation described in a companion paper.
[13]  [pdf] - 1981272
STRIDES: A 3.9 per cent measurement of the Hubble constant from the strong lens system DES J0408-5354
Comments: 30 pages, 17 figures, 6 tables. Submitted to MNRAS. This version: minor typo fixes
Submitted: 2019-10-14, last modified: 2019-10-15
We present a blind time-delay cosmographic analysis for the lens system DES J0408$-$5354. This system is extraordinary for the presence of two sets of multiple images at different redshifts, which provide the opportunity to obtain more information at the cost of increased modelling complexity with respect to previously analyzed systems. We perform detailed modelling of the mass distribution for this lens system using three band Hubble Space Telescope imaging. We combine the measured time delays, line-of-sight central velocity dispersion of the deflector, and statistically constrained external convergence with our lens models to estimate two cosmological distances. We measure the "effective" time-delay distance corresponding to the redshifts of the deflector and the lensed quasar $D_{\Delta t}^{\rm eff}=3382^{+146}_{-115}$ Mpc and the angular diameter distance to the deflector $D_{\rm d}=1711^{+376}_{-280}$ Mpc, with covariance between the two distances. From these constraints on the cosmological distances, we infer the Hubble constant $H_0 = 74.2^{+2.7}_{-3.0}$ km s$^{-1}$ Mpc$^{-1}$ assuming a flat $\Lambda$CDM cosmology and a uniform prior for $\Omega_{\rm m}$ as $\Omega_{\rm m} \sim \mathcal{U}(0.05, 0.5)$. This measurement gives the most precise constraint on $H_0$ to date from a single lens. Our measurement is consistent with that obtained from the previous sample of six lenses analyzed by the $H_0$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration. It is also consistent with measurements of $H_0$ based on the local distance ladder, reinforcing the tension with the inference from early universe probes, for example, with 2.2$\sigma$ discrepancy from the cosmic microwave background measurement.
[14]  [pdf] - 1963750
Cosmic dissonance: new physics or systematics behind a short sound horizon?
Comments: 9 pages, 6 figures, 3 tables; submitted to A&A
Submitted: 2019-09-17
Persistent tension between low-redshift and the Cosmic Microwave Background (CMB) measurements, in terms of the sound horizon and the Hubble-Lema\^itre constant, suggests new physics beyond the Standard Model, departures from concordance cosmology, or residual systematics. Assessing the likelihood of new physics devised to resolve the tension requires thorough consistency tests of several independent and high-precision distance calibrations. We examine recent updated distance calibrations from Cepheids, gravitational lensing time-delay observations, and the Tip of the Red Giant Branch. Through selected cosmographic methods for combining observations of the Baryon Acoustic Oscillations (BAO), Type Ia supernovae and local distance calibrators, we obtain robust measurements of the sound horizon, independent of cosmological models and the CMB. For our different models and probes, the sound horizon scale varies between $r_{\rm s}=(135\pm3)$ and $r_{\rm s}=(140\pm3)$~Mpc. The combined tension in the sound horizon and the Hubble-Lema\^itre constant ranges between $3\sigma$ and $5\sigma$, depending on whether distance ladder calibrations are used, and is independent of possible changes to the late-Universe expansion history. Early-Universe extensions with $\textrm{N}_{\textrm{eff}}=3.27\pm0.15$ are allowed by the data, unless a Cepheid distance calibration is included. Some late-dark-energy models, as constrained by the relative supernova distance moduli, do not resolve the current tension. Results from time-delay lenses are consistent with those from distance-ladder calibrations. With upcoming lens samples, the tension in the inferred sound horizon will be completely independent from local calibrations. New proposals to resolve the tension should examine CMB and late-Universe constraints separately, and the inference on both $H_0$ and $r_s.$
[15]  [pdf] - 2025776
Constraints on the mass-concentration relation of cold dark matter halos with 11 strong gravitational lenses
Comments: fix typos, edit wording
Submitted: 2019-09-05, last modified: 2019-09-09
The mass-concentration relation of dark matter halos reflects the assembly history of objects in hierarchical structure formation scenarios, and depends on fundamental quantities in cosmology such as the slope of the primordial matter power-spectrum. This relation is unconstrained by observations on sub-galactic scales. We derive the first measurement of the mass-concentration relation using the image positions and flux ratios from eleven quadruple-image strong gravitational lenses (quads) in the mass range $10^{6} - 10^{10} M_{\odot}$, assuming cold dark matter. Our analysis framework includes both subhalos and line of sight halos, marginalizes over nuisance parameters describing the lens macromodel, accounts for finite source effects on lensing observables, and simultaneously constrains the normalization and logarithmic slope of the mass-concentration relation, and the normalization of the subhalo mass function. At $z=0$, we constrain the concentration of $10^{8} M_{\odot}$ halos $c=12_{-5}^{+6}$ at $68 \%$ CI, and $c=12_{-9}^{+15}$ at $95 \%$ CI. For a $10^{7} M_{\odot}$ halo, we obtain $68 \%$ ($95 \%$) constraints $c=15_{-8}^{+9}$ ($c=15_{-11}^{+18}$), while for $10^{9} M_{\odot}$ halos $c=10_{-4}^{+7}$ ($c=10_{-7}^{+14}$). These results are consistent with the theoretical predictions from mass-concentration relations in the literature, and establish strong lensing by galaxies as a powerful probe of halo concentrations on sub-galactic scales across cosmological distance.
[16]  [pdf] - 1956501
H0LiCOW X: Spectroscopic/imaging survey and galaxy-group identification around the strong gravitational lens system WFI2033-4723
Comments: Matches the version accepted for publication by MNRAS. Note that this paper previously appeared as H0LICOW XI
Submitted: 2019-05-21, last modified: 2019-09-06
Galaxies and galaxy groups located along the line of sight towards gravitationally lensed quasars produce high-order perturbations of the gravitational potential at the lens position. When these perturbation are too large, they can induce a systematic error on $H_0$ of a few-percent if the lens system is used for cosmological inference and the perturbers are not explicitly accounted for in the lens model. In this work, we present a detailed characterization of the environment of the lens system WFI2033-4723 ($z_{\rm src} = 1.662$, $z_{\rm lens}$ = 0.6575), one of the core targets of the H0LICOW project for which we present cosmological inferences in a companion paper (Rusu et al. 2019). We use the Gemini and ESO-Very Large telescopes to measure the spectroscopic redshifts of the brightest galaxies towards the lens, and use the ESO-MUSE integral field spectrograph to measure the velocity-dispersion of the lens ($\sigma_{\rm {los}}= 250^{+15}_{-21}$ km/s) and of several nearby galaxies. In addition, we measure photometric redshifts and stellar masses of all galaxies down to $i < 23$ mag, mainly based on Dark Energy Survey imaging (DR1). Our new catalog, complemented with literature data, more than doubles the number of known galaxy spectroscopic redshifts in the direct vicinity of the lens, expanding to 116 (64) the number of spectroscopic redshifts for galaxies separated by less than 3 arcmin (2 arcmin) from the lens. Using the flexion-shift as a measure of the amplitude of the gravitational perturbation, we identify 2 galaxy groups and 3 galaxies that require specific attention in the lens models. The ESO MUSE data enable us to measure the velocity-dispersions of three of these galaxies. These results are essential for the cosmological inference analysis presented in Rusu et al. (2019).
[17]  [pdf] - 1945774
Astrometric requirements for strong lensing time-delay cosmography
Comments: 7 pages, 1 figure, 3 tables, MNRAS accepted. Notebook to perform the calculations presented in this work available here:
Submitted: 2019-04-24, last modified: 2019-08-07
The time delay between the arrival of photons of multiple images of time variable sources can be used to constrain absolute distances in the Universe (Refsdal 1964), and in turn obtain a direct estimate of the Hubble constant and other cosmological parameters. To convert the time delay into distances, it is well known that the gravitational potential of the main deflector and the contribution of the matter along the line-of-sight need to be known to a sufficient level of precision. In this paper, we discuss a new astrometric requirement that is becoming important as time-delay cosmography improves in precision and accuracy with larger samples, and better data and modelling techniques. We derive an analytic expression for the propagation of astrometric uncertainties on the multiple image positions into the inference of the Hubble constant and derive requirements depending on image separation and relative time delay. We note that this requirement applies equally to the image position measurements and to the accuracy of the model in reproducing them. To illustrate the requirement, we discuss some example lensing configurations and highlight that, especially for time delays of order 10 days or shorter, the relative astrometric requirement is of order milli-arcseconds, setting a tight requirement on both measurements and models. With current optical infrared technology, astrometric uncertainties may be the dominant limitation for strong lensing cosmography in the small image-separation regime when high-precision time-delays become accessible.
[18]  [pdf] - 1945799
The Hubble Constant determined through an inverse distance ladder including quasar time delays and Type Ia supernovae
Comments: 5 pages, 3 figures, A&A letters accepted version
Submitted: 2019-05-29, last modified: 2019-07-31
Context. The precise determination of the present-day expansion rate of the Universe, expressed through the Hubble constant $H_0$, is one of the most pressing challenges in modern cosmology. Assuming flat $\Lambda$CDM, $H_0$ inference at high redshift using cosmic-microwave-background data from Planck disagrees at the 4.4$\sigma$ level with measurements based on the local distance ladder made up of parallaxes, Cepheids and Type Ia supernovae (SNe Ia), often referred to as "Hubble tension". Independent, cosmological-model-insensitive ways to infer $H_0$ are of critical importance. Aims. We apply an inverse-distance-ladder approach, combining strong-lensing time-delay-distance measurements with SN Ia data. By themselves, SNe Ia are merely good relative distance indicators, but by anchoring them to strong gravitational lenses one can obtain an $H_0$ measurement that is relatively insensitive to other cosmological parameters. Methods. A cosmological parameter estimate is performed for different cosmological background models, both for strong-lensing data alone and for the combined lensing + SNe Ia data sets. Results. The cosmological-model dependence of strong-lensing $H_0$ measurements is significantly mitigated through the inverse distance ladder. In combination with SN Ia data, the inferred $H_0$ consistently lies around 73-74 km s$^{-1}$ Mpc$^{-1}$, regardless of the assumed cosmological background model. Our results agree nicely with those from the local distance ladder, but there is a >2$\sigma$ tension with Planck results, and a ~1.5$\sigma$ discrepancy with results from an inverse distance ladder including Planck, Baryon Acoustic Oscillations and SNe Ia. Future strong-lensing distance measurements will reduce the uncertainties in $H_0$ from our inverse distance ladder.
[19]  [pdf] - 1924076
Probing dark matter structure down to $10^7$ solar masses: flux ratio statistics in gravitational lenses with line of sight halos
Comments: matches version accepted to MNRAS
Submitted: 2019-01-30, last modified: 2019-07-26
Strong lensing provides a powerful means of investigating the nature of dark matter as it probes dark matter structure on sub-galactic scales. We present an extension of a forward modeling framework that uses flux ratios from quadruply imaged quasars (quads) to measure the shape and amplitude of the halo mass function, including line of sight (LOS) halos and main deflector subhalos. We apply this machinery to 50 mock lenses --- roughly the number of known quads --- with warm dark matter (WDM) mass functions exhibiting free-streaming cutoffs parameterized by the half-mode mass $m_{\rm{hm}}$. Assuming cold dark matter (CDM), we forecast bounds on $m_{\rm{hm}}$ and the corresponding thermal relic particle masses over a range of tidal destruction severity, assuming a particular WDM mass function and mass-concentration relation. With significant tidal destruction, at $2 \sigma$ we constrain $m_{\rm{hm}}<10^{7.9} \left(10^{8.4}\right) M_{\odot}$, or a 4.4 (3.1) keV thermal relic, with image flux uncertainties from measurements and lens modeling of $2\% \left(6\%\right)$. With less severe tidal destruction we constrain $m_{\rm{hm}}<10^{7} \left(10^{7.4}\right) M_{\odot}$, or an 8.2 (6.2) keV thermal relic. If dark matter is warm, with $m_{\rm{hm}} = 10^{7.7} M_{\odot}$ (5.1 keV), we would favor WDM with $m_{\rm{hm}} > 10^{7.7} M_{\odot}$ over CDM with relative likelihoods of 22:1 and 8:1 with flux uncertainties of $2\%$ and $6\%$, respectively. These bounds improve over those obtained by modeling only main deflector subhalos because LOS objects produce additional flux perturbations, especially for high redshift systems. These results indicate that $\sim 50$ quads can conclusively differentiate between warm and cold dark matter.
[20]  [pdf] - 1966838
A SHARP view of H0LiCOW: $H_{0}$ from three time-delay gravitational lens systems with adaptive optics imaging
Comments: 33 pages, 23 figures
Submitted: 2019-07-04
We present the measurement of the Hubble Constant, $H_0$, with three strong gravitational lens systems. We describe a blind analysis of both PG1115+080 and HE0435-1223 as well as an extension of our previous analysis of RXJ1131-1231. For each lens, we combine new adaptive optics (AO) imaging from the Keck Telescope, obtained as part of the SHARP AO effort, with Hubble Space Telescope (HST) imaging, velocity dispersion measurements, and a description of the line-of-sight mass distribution to build an accurate and precise lens mass model. This mass model is then combined with the COSMOGRAIL measured time delays in these systems to determine $H_{0}$. We do both an AO-only and an AO+HST analysis of the systems and find that AO and HST results are consistent. After unblinding, the AO-only analysis gives $H_{0}=82.8^{+9.4}_{-8.3}~\rm km\,s^{-1}\,Mpc^{-1}$ for PG1115+080, $H_{0}=70.1^{+5.3}_{-4.5}~\rm km\,s^{-1}\,Mpc^{-1}$ for HE0435-1223, and $H_{0}=77.0^{+4.0}_{-4.6}~\rm km\,s^{-1}\,Mpc^{-1}$ for RXJ1131-1231. The joint AO-only result for the three lenses is $H_{0}=75.6^{+3.2}_{-3.3}~\rm km\,s^{-1}\,Mpc^{-1}$. The joint result of the AO+HST analysis for the three lenses is $H_{0}=76.8^{+2.6}_{-2.6}~\rm km\,s^{-1}\,Mpc^{-1}$. All of the above results assume a flat $\Lambda$ cold dark matter cosmology with a uniform prior on $\Omega_{\textrm{m}}$ in [0.05, 0.5] and $H_{0}$ in [0, 150] $\rm km\,s^{-1}\,Mpc^{-1}$. This work is a collaboration of the SHARP and H0LiCOW teams, and shows that AO data can be used as the high-resolution imaging component in lens-based measurements of $H_0$. The full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper.
[21]  [pdf] - 1888804
H0LiCOW - IX. Cosmographic analysis of the doubly imaged quasar SDSS 1206+4332 and a new measurement of the Hubble constant
Comments: 30 pages, 17 figures, MNRAS published likelihood available here:, all modeling and analysis scripts available upon request
Submitted: 2018-09-04, last modified: 2019-05-24
We present a blind time-delay strong lensing (TDSL) cosmographic analysis of the doubly imaged quasar SDSS 1206+4332. We combine the relative time delay between the quasar images, Hubble Space Telescope imaging, the Keck stellar velocity dispersion of the lensing galaxy, and wide-field photometric and spectroscopic data of the field to constrain two angular diameter distance relations. The combined analysis is performed by forward modelling the individual data sets through a Bayesian hierarchical framework, and it is kept blind until the very end to prevent experimenter bias. After unblinding, the inferred distances imply a Hubble constant $H_0 = 68.8^{+5.4}_{-5.1}$ kms$^{-1}$Mpc$^{-1}$, assuming a flat Lambda cold dark matter cosmology with uniform prior on $\Omega_{\rm m}$ in [0.05, 0.5]. The precision of our cosmographic measurement with the doubly imaged quasar SDSS 1206+4332 is comparable with those of quadruply imaged quasars and opens the path to perform on selected doubles the same analysis as anticipated for quads. Our analysis is based on a completely independent lensing code than our previous three H0LiCOW systems and the new measurement is fully consistent with those. We provide the analysis scripts paired with the publicly available software to facilitate independent analysis. The consistency between blind measurements with independent codes provides an important sanity check on lens modelling systematics. By combining the likelihoods of the four systems under the same prior, we obtain $H_0 = 72.5^{+2.1}_{-2.3}$kms$^{-1}$Mpc$^{-1}$. This measurement is independent of the distance ladder and other cosmological probes.
[22]  [pdf] - 1872765
Probing the Fundamental Nature of Dark Matter with the Large Synoptic Survey Telescope
Drlica-Wagner, Alex; Mao, Yao-Yuan; Adhikari, Susmita; Armstrong, Robert; Banerjee, Arka; Banik, Nilanjan; Bechtol, Keith; Bird, Simeon; Boddy, Kimberly K.; Bonaca, Ana; Bovy, Jo; Buckley, Matthew R.; Bulbul, Esra; Chang, Chihway; Chapline, George; Cohen-Tanugi, Johann; Cuoco, Alessandro; Cyr-Racine, Francis-Yan; Dawson, William A.; Rivero, Ana Díaz; Dvorkin, Cora; Erkal, Denis; Fassnacht, Christopher D.; García-Bellido, Juan; Giannotti, Maurizio; Gluscevic, Vera; Golovich, Nathan; Hendel, David; Hezaveh, Yashar D.; Horiuchi, Shunsaku; Jee, M. James; Kaplinghat, Manoj; Keeton, Charles R.; Koposov, Sergey E.; Lam, Casey Y.; Li, Ting S.; Lu, Jessica R.; Mandelbaum, Rachel; McDermott, Samuel D.; McNanna, Mitch; Medford, Michael; Meyer, Manuel; Marc, Moniez; Murgia, Simona; Nadler, Ethan O.; Necib, Lina; Nuss, Eric; Pace, Andrew B.; Peter, Annika H. G.; Polin, Daniel A.; Prescod-Weinstein, Chanda; Read, Justin I.; Rosenfeld, Rogerio; Shipp, Nora; Simon, Joshua D.; Slatyer, Tracy R.; Straniero, Oscar; Strigari, Louis E.; Tollerud, Erik; Tyson, J. Anthony; Wang, Mei-Yu; Wechsler, Risa H.; Wittman, David; Yu, Hai-Bo; Zaharijas, Gabrijela; Ali-Haïmoud, Yacine; Annis, James; Birrer, Simon; Biswas, Rahul; Blazek, Jonathan; Brooks, Alyson M.; Buckley-Geer, Elizabeth; Caputo, Regina; Charles, Eric; Digel, Seth; Dodelson, Scott; Flaugher, Brenna; Frieman, Joshua; Gawiser, Eric; Hearin, Andrew P.; Hložek, Renee; Jain, Bhuvnesh; Jeltema, Tesla E.; Koushiappas, Savvas M.; Lisanti, Mariangela; LoVerde, Marilena; Mishra-Sharma, Siddharth; Newman, Jeffrey A.; Nord, Brian; Nourbakhsh, Erfan; Ritz, Steven; Robertson, Brant E.; Sánchez-Conde, Miguel A.; Slosar, Anže; Tait, Tim M. P.; Verma, Aprajita; Vilalta, Ricardo; Walter, Christopher W.; Yanny, Brian; Zentner, Andrew R.
Comments: 96 pages, 22 figures, 1 table
Submitted: 2019-02-04, last modified: 2019-04-24
Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at .
[23]  [pdf] - 1847070
Testing the Nature of Dark Matter with Extremely Large Telescopes
Comments: 8 pages, 2 figures. Science white paper submitted to the Astro2020 decadal survey
Submitted: 2019-03-12
For nearly 40 years, dark matter has been widely assumed to be cold and collisionless. Cold dark matter models make fundamental predictions for the behavior of dark matter on small (<10 kpc) scales. These predictions include cuspy density profiles at the centers of dark matter halos and a halo mass function that increases as dN/dM ~ M^-1.9 down to very small masses. We suggest two observational programs relying on extremely large telescopes to critically test these predictions, and thus shed new light on the nature of dark matter. (1) Combining adaptive optics-enabled imaging with deep spectroscopy to measure the three-dimensional motions of stars within a sample of Local Group dwarf galaxies that are the cleanest dark matter laboratories known in the nearby universe. From these observations the inner slope of the dark matter density profile can be determined with an accuracy of 0.20 dex, enabling a central cusp to be distinguished from a core at 5 sigma significance. (2) Diffraction-limited AO imaging and integral field spectroscopy of gravitationally lensed galaxies and quasars to quantify the abundance of dark substructures in the halos of the lens galaxies and along the line of sight. Observations of 50 lensed arcs and 50 multiply-imaged quasars will be sufficient to measure the halo mass function over the range 10^7 < M < 10^10 Msun at cosmological scales, independent of the baryonic and stellar composition of those structures. These two observational probes provide complementary information about the small scale structure, with a joint self-consistent analysis mitigating limitations of either probe. This program will produce the strongest existing constraints on the properties of dark matter on small scales, allowing conclusive tests of alternative warm, fuzzy, and self-interacting dark matter models.
[24]  [pdf] - 1847119
Measuring the Hubble Constant Near and Far in the Era of ELT's
Comments: Submitted as an Astro2020 White Paper. Please send comments to both Rachael Beaton & Simon Birrer. Development of this paper occurred as part of the The US Extremely Large Telescope Program Workshop in Oct 2018. We wish to acknowledge NOAO for bringing the co-authors together, in particular the enthusiasm and tireless leadership of Mark Dickinson
Submitted: 2019-03-12
Many of the fundamental physical constants in Physics, as a discipline, are measured to exquisite levels of precision. The fundamental constants that define Cosmology, however, are largely determined via a handful of independent techniques that are applied to even fewer datasets. The history of the measurement of the Hubble Constant (H0), which serves to anchor the expansion history of the Universe to its current value, is an exemplar of the difficulties of cosmological measurement; indeed, as we approach the centennial of its first measurement, the quest for H0 still consumes a great number of resources. In this white paper, we demonstrate how the approaching era of Extremely Large Telescopes (ELTs) will transform the astrophysical measure of H0 from the limited and few into a fundamentally new regime where (i) multiple, independent techniques are employed with modest use of large aperture facilities and (ii) 1% or better precision is readily attainable. This quantum leap in how we approach H0 is due to the unparalleled sensitivity and spatial resolution of ELT's and the ability to use integral field observations for simultaneous spectroscopy and photometry, which together permit both familiar and new techniques to effectively by-pass the conventional 'ladder' framework to minimize total uncertainty. Three independent techniques are discussed -- (i) standard candles via a two-step distance ladder applied to metal, poor stellar populations, (ii) standard clocks via gravitational lens cosmography, and (iii) standard sirens via gravitational wave sources -- each of which can reach 1% with relatively modest investment from 30-m class facilities.
[25]  [pdf] - 1846071
Dark Matter Science in the Era of LSST
Bechtol, Keith; Drlica-Wagner, Alex; Abazajian, Kevork N.; Abidi, Muntazir; Adhikari, Susmita; Ali-Haïmoud, Yacine; Annis, James; Ansarinejad, Behzad; Armstrong, Robert; Asorey, Jacobo; Baccigalupi, Carlo; Banerjee, Arka; Banik, Nilanjan; Bennett, Charles; Beutler, Florian; Bird, Simeon; Birrer, Simon; Biswas, Rahul; Biviano, Andrea; Blazek, Jonathan; Boddy, Kimberly K.; Bonaca, Ana; Borrill, Julian; Bose, Sownak; Bovy, Jo; Frye, Brenda; Brooks, Alyson M.; Buckley, Matthew R.; Buckley-Geer, Elizabeth; Bulbul, Esra; Burchat, Patricia R.; Burgess, Cliff; Calore, Francesca; Caputo, Regina; Castorina, Emanuele; Chang, Chihway; Chapline, George; Charles, Eric; Chen, Xingang; Clowe, Douglas; Cohen-Tanugi, Johann; Comparat, Johan; Croft, Rupert A. C.; Cuoco, Alessandro; Cyr-Racine, Francis-Yan; D'Amico, Guido; Davis, Tamara M; Dawson, William A.; de la Macorra, Axel; Di Valentino, Eleonora; Rivero, Ana Díaz; Digel, Seth; Dodelson, Scott; Doré, Olivier; Dvorkin, Cora; Eckner, Christopher; Ellison, John; Erkal, Denis; Farahi, Arya; Fassnacht, Christopher D.; Ferreira, Pedro G.; Flaugher, Brenna; Foreman, Simon; Friedrich, Oliver; Frieman, Joshua; García-Bellido, Juan; Gawiser, Eric; Gerbino, Martina; Giannotti, Maurizio; Gill, Mandeep S. S.; Gluscevic, Vera; Golovich, Nathan; Gontcho, Satya Gontcho A; González-Morales, Alma X.; Grin, Daniel; Gruen, Daniel; Hearin, Andrew P.; Hendel, David; Hezaveh, Yashar D.; Hirata, Christopher M.; Hložek, Renee; Horiuchi, Shunsaku; Jain, Bhuvnesh; Jee, M. James; Jeltema, Tesla E.; Kamionkowski, Marc; Kaplinghat, Manoj; Keeley, Ryan E.; Keeton, Charles R.; Khatri, Rishi; Koposov, Sergey E.; Koushiappas, Savvas M.; Kovetz, Ely D.; Lahav, Ofer; Lam, Casey; Lee, Chien-Hsiu; Li, Ting S.; Liguori, Michele; Lin, Tongyan; Lisanti, Mariangela; LoVerde, Marilena; Lu, Jessica R.; Mandelbaum, Rachel; Mao, Yao-Yuan; McDermott, Samuel D.; McNanna, Mitch; Medford, Michael; Meerburg, P. Daniel; Meyer, Manuel; Mirbabayi, Mehrdad; Mishra-Sharma, Siddharth; Marc, Moniez; More, Surhud; Moustakas, John; Muñoz, Julian B.; Murgia, Simona; Myers, Adam D.; Nadler, Ethan O.; Necib, Lina; Newburgh, Laura; Newman, Jeffrey A.; Nord, Brian; Nourbakhsh, Erfan; Nuss, Eric; O'Connor, Paul; Pace, Andrew B.; Padmanabhan, Hamsa; Palmese, Antonella; Peiris, Hiranya V.; Peter, Annika H. G.; Piacentni, Francesco; Piacentini, Francesco; Plazas, Andrés; Polin, Daniel A.; Prakash, Abhishek; Prescod-Weinstein, Chanda; Read, Justin I.; Ritz, Steven; Robertson, Brant E.; Rose, Benjamin; Rosenfeld, Rogerio; Rossi, Graziano; Samushia, Lado; Sánchez, Javier; Sánchez-Conde, Miguel A.; Schaan, Emmanuel; Sehgal, Neelima; Senatore, Leonardo; Seo, Hee-Jong; Shafieloo, Arman; Shan, Huanyuan; Shipp, Nora; Simon, Joshua D.; Simon, Sara; Slatyer, Tracy R.; Slosar, Anže; Sridhar, Srivatsan; Stebbins, Albert; Straniero, Oscar; Strigari, Louis E.; Tait, Tim M. P.; Tollerud, Erik; Troxel, M. A.; Tyson, J. Anthony; Uhlemann, Cora; Urenña-López, L. Arturo; Verma, Aprajita; Vilalta, Ricardo; Walter, Christopher W.; Wang, Mei-Yu; Watson, Scott; Wechsler, Risa H.; Wittman, David; Xu, Weishuang; Yanny, Brian; Young, Sam; Yu, Hai-Bo; Zaharijas, Gabrijela; Zentner, Andrew R.; Zuntz, Joe
Comments: 11 pages, 2 figures, Science Whitepaper for Astro 2020, more information at
Submitted: 2019-03-11
Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.
[26]  [pdf] - 1821724
Is every strong lens model unhappy in its own way? Uniform modelling of a sample of 13 quadruply+ imaged quasars
Comments: 22 pages, 9 figures, 8 tables, Accepted to MNRAS. This version: accepted version
Submitted: 2018-07-24, last modified: 2019-01-28
Strong-gravitational lens systems with quadruply-imaged quasars (quads) are unique probes to address several fundamental problems in cosmology and astrophysics. Although they are intrinsically very rare, ongoing and planned wide-field deep-sky surveys are set to discover thousands of such systems in the next decade. It is thus paramount to devise a general framework to model strong-lens systems to cope with this large influx without being limited by expert investigator time. We propose such a general modelling framework (implemented with the publicly available software Lenstronomy) and apply it to uniformly model three-band Hubble Space Telescope Wide Field Camera 3 images of 13 quads. This is the largest uniformly modelled sample of quads to date and paves the way for a variety of studies. To illustrate the scientific content of the sample, we investigate the alignment between the mass and light distribution in the deflectors. The position angles of these distributions are well-aligned, except when there is strong external shear. However, we find no correlation between the ellipticity of the light and mass distributions. We also show that the observed flux-ratios between the images depart significantly from the predictions of simple smooth models. The departures are strongest in the bluest band, consistent with microlensing being the dominant cause in addition to millilensing. Future papers will exploit this rich dataset in combination with ground based spectroscopy and time delays to determine quantities such as the Hubble constant, the free streaming length of dark matter, and the normalization of the initial stellar mass function.
[27]  [pdf] - 1783742
Lenstronomy: multi-purpose gravitational lens modelling software package
Comments: software and documentation see, matches accepted version of Physics of the Dark Universe
Submitted: 2018-03-26, last modified: 2018-11-12
We present Lenstronomy, a multi-purpose open-source gravitational lens modeling python package. Lenstronomy is able to reconstruct the lens mass and surface brightness distributions of strong lensing systems using forward modelling. Lenstronomy supports a wide range of analytic lens and light models in arbitrary combination. The software is also able to reconstruct complex extended sources (Birrer et. al 2015) as well as being able to model point sources. We designed Lenstronomy to be stable, flexible and numerically accurate, with a clear user interface that could be deployed across different platforms. Throughout its development, we have actively used Lenstronomy to make several measurements including deriving constraints on dark matter properties in strong lenses, measuring the expansion history of the universe with time-delay cosmography, measuring cosmic shear with Einstein rings and decomposing quasar and host galaxy light. The software is distributed under the MIT license. The documentation, starter guide, example notebooks, source code and installation guidelines can be found at
[28]  [pdf] - 1838190
Sparse Lens Inversion Technique (SLIT): lens and source separability from linear inversion of the source reconstruction problem
Comments: 17 pages, 10 figures, Accepted (A&A)
Submitted: 2018-09-24
Strong gravitational lensing offers a wealth of astrophysical information on the background source it affects, provided the lensed source can be reconstructed as if it was seen in the absence of lensing. In the present work, we illustrate how sparse optimisation can address the problem. As a first step towards a full free-form lens modelling technique, we consider linear inversion of the lensed source under sparse regularisation and joint deblending from the lens light profile. The method is based on morphological component analysis, assuming a known mass model. We show with numerical experiments that representing the lens and source light using an undecimated wavelet basis allows us to reconstruct the source and to separate it from the foreground lens at the same time. Both the source and lens light have a non-analytic form, allowing for the flexibility needed in the inversion to represent arbitrarily small and complex luminous structures in the lens and source. in addition, sparse regularisation avoids over-fitting the data and does not require the use of any adaptive mesh or pixel grid. As a consequence, our reconstructed sources can be represented on a grid of very small pixels. Sparse regularisation in the wavelet domain also allows for automated computation of the regularisation parameter, thus minimising the impact of arbitrary choice of initial parameters. Our inversion technique for a fixed mass distribution can be incorporated in future lens modelling technique iterating over the lens mass parameters. The python package corresponding to the algorithms described in this article can be downloaded via the github platform at
[29]  [pdf] - 1745972
The STRong lensing Insights into the Dark Energy Survey (STRIDES) 2016 follow-up campaign. I. Overview and classification of candidates selected by two techniques
Comments: (MNRAS in press; color scheme of Figures 1-4 improved with respect to v1)
Submitted: 2018-08-14, last modified: 2018-09-05
The primary goals of the STRong lensing Insights into the Dark Energy Survey (STRIDES) collaboration are to measure the dark energy equation of state parameter and the free streaming length of dark matter. To this aim, STRIDES is discovering strongly lensed quasars in the imaging data of the Dark Energy Survey and following them up to measure time delays, high resolution imaging, and spectroscopy sufficient to construct accurate lens models. In this paper, we first present forecasts for STRIDES. Then, we describe the STRIDES classification scheme, and give an overview of the Fall 2016 follow-up campaign. We continue by detailing the results of two selection methods, the Outlier Selection Technique and a morphological algorithm, and presenting lens models of a system, which could possibly be a lensed quasar in an unusual configuration. We conclude with the summary statistics of the Fall 2016 campaign. Including searches presented in companion papers (Anguita et al.; Ostrovski et al.), STRIDES followed up 117 targets identifying 7 new strongly lensed systems, and 7 nearly identical quasars (NIQs), which could be confirmed as lenses by the detection of the lens galaxy. 76 candidates were rejected and 27 remain otherwise inconclusive, for a success rate in the range 6-35\%. This rate is comparable to that of previous searches like SQLS even though the parent dataset of STRIDES is purely photometric and our selection of candidates cannot rely on spectroscopic information.
[30]  [pdf] - 1739845
Probing the nature of dark matter by forward modeling flux ratios in strong gravitational lenses
Comments: accepted to MNRAS; added a new simulation to illustrate how constraints on free streaming length scale with mass function normalization; 17 pages, 12 figures, 1 table, 5 appendices
Submitted: 2017-12-13, last modified: 2018-08-16
The free streaming length of dark matter particles determines the abundance of structure on sub-galactic scales. We present a statistical technique, amendable to any parameterization of subhalo density profile and mass function, to probe dark matter on these scales with quadrupole image lenses. We consider a warm dark matter particle with a mass function characterized by a normalization and free streaming scale $m_{\rm{hm}}$. We forecast bounds on dark matter warmth for 120-180 lenses, attainable with future surveys, at typical lens (source) redshifts of 0.5 (1.5) in early-type galaxies with velocity dispersions of 220-270 km/sec. We demonstrate that limits on $m_{\rm{hm}}$ deteriorate rapidly with increasing uncertainty in image fluxes, underscoring the importance of precise measurements and accurate lens models. For our forecasts, we assume the deflectors in the lens sample do not exhibit complex morphologies, so we neglect systematic errors in their modeling. Omitting the additional signal from line of sight halos, our constraints underestimate the true power of the method. Assuming cold dark matter, for a low normalization, corresponding the destruction of all subhalos within the host scale radius, we forecast $2\sigma$ bounds on $m_{\rm{hm}}$ (thermal relic mass) of $10^{7.5} \ (5.0)$, $10^{8} \ (3.6)$, and $10^{8.5} \ (2.7) \ M_{\odot} \left(\rm{keV}\right)$ for flux errors of $2\%$, $4\%$, and $8\%$. With a higher normalization, these constraints improve to $10^{7.2} \ (6.6)$, $10^{7.5} \ (5.3) $, and $10^{7.8} \ (4.3) \ M_{\odot} \left(\rm{keV}\right)$ with 120 systems. We are also able to measure the normalization of the mass function, which has implications for baryonic feedback models and tidal stripping.
[31]  [pdf] - 1614837
Cosmic Shear with Einstein Rings
Comments: 4 pages, 1 figure, 1 table. ApJL accepted
Submitted: 2017-10-03, last modified: 2017-12-15
We explore a new technique to measure cosmic shear using Einstein rings. In Birrer et al. (2017), we showed that the detailed modelling of Einstein rings can be used to measure external shear to high precision. In this letter, we explore how a collection of Einstein rings can be used as a statistical probe of cosmic shear. We present a forecast of the cosmic shear information available in Einstein rings for different strong lensing survey configurations. We find that, assuming that the number density of Einstein rings in the COSMOS survey is representative, future strong lensing surveys should have a cosmological precision comparable to the current ground based weak lensing surveys. We discuss how this technique is complementary to the standard cosmic shear analyses since it is sensitive to different systematic and can be used for cross-calibration.
[32]  [pdf] - 1582109
Core or cusps: The central dark matter profile of a redshift one strong lensing cluster with a bright central image
Comments: 13 Pages. Accepted for publication in ApJ
Submitted: 2017-03-24, last modified: 2017-06-02
We report on SPT-CLJ2011-5228, a giant system of arcs created by a cluster at $z=1.06$. The arc system is notable for the presence of a bright central image. The source is a Lyman Break galaxy at $z_s=2.39$ and the mass enclosed within the 14 arc second radius Einstein ring is $10^{14.2}$ solar masses. We perform a full light profile reconstruction of the lensed images to precisely infer the parameters of the mass distribution. The brightness of the central image demands that the central total density profile of the lens be shallow. By fitting the dark matter as a generalized Navarro-Frenk-White profile---with a free parameter for the inner density slope---we find that the break radius is $270^{+48}_{-76}$ kpc, and that the inner density falls with radius to the power $-0.38\pm0.04$ at 68 percent confidence. Such a shallow profile is in strong tension with our understanding of relaxed cold dark matter halos; dark matter only simulations predict the inner density should fall as $r^{-1}$. The tension can be alleviated if this cluster is in fact a merger; a two halo model can also reconstruct the data, with both clumps (density going as $r^{-0.8}$ and $r^{-1.0}$) much more consistent with predictions from dark matter only simulations. At the resolution of our Dark Energy Survey imaging, we are unable to choose between these two models, but we make predictions for forthcoming Hubble Space Telescope imaging that will decisively distinguish between them.
[33]  [pdf] - 1581334
Lensing substructure quantification in RXJ1131-1231: A 2 keV lower bound on dark matter thermal relic mass
Comments: 23 pages, JCAP published
Submitted: 2017-01-31, last modified: 2017-05-18
We study the substructure content of the strong gravitational lens RXJ1131-1231 through a forward modelling approach that relies on generating an extensive suite of realistic simulations. We use a semi-analytic merger tree prescription that allows us to stochastically generate substructure populations whose properties depend on the dark matter particle mass. These synthetic halos are then used as lenses to produce realistic mock images that have the same features, e.g. luminous arcs, quasar positions, instrumental noise and PSF, as the data. We then analyse the data and the simulations in the same way with summary statistics that are sensitive to the signal being targeted and are able to constrain models of dark matter statistically using Approximate Bayesian Computing (ABC) techniques. In this work, we focus on the thermal relic mass estimate and fix the semi-analytic descriptions of the substructure evolution based on recent literature. We are able, based on the HST data for RXJ1131-1231, to rule out a warm dark matter thermal relic mass below 2 keV at the 2$\sigma$ confidence level.
[34]  [pdf] - 1570153
Line-of-sight effects in strong lensing: Putting theory into practice
Comments: 16 pages, 9 figures, accepted for publication in JCAP
Submitted: 2016-10-05, last modified: 2017-04-28
We present a simple method to accurately infer line of sight (LOS) integrated lensing effects for galaxy scale strong lens systems through image reconstruction. Our approach enables us to separate weak lensing LOS effects from the main strong lens deflector. We test our method using mock data and show that strong lens systems can be accurate probes of cosmic shear with a precision on the shear terms of $\pm 0.003$ (statistical error) for an HST-like dataset. We apply our formalism to reconstruct the lens COSMOS 0038+4133 and its LOS. In addition, we estimate the LOS properties with a halo-rendering estimate based on the COSMOS field galaxies and a galaxy-halo connection. The two approaches are independent and complementary in their information content. We find that when estimating the convergence at the strong lens system, performing a joint analysis improves the measure by a factor of two compared to a halo model only analysis. Furthermore the constraints of the strong lens reconstruction lead to tighter constraints on the halo masses of the LOS galaxies. Joint constraints of multiple strong lens systems may add valuable information to the galaxy-halo connection and may allow independent weak lensing shear measurement calibrations.
[35]  [pdf] - 1581349
Models of the strongly lensed quasar DES J0408-5354
Comments: 12 pages, 8 figures, MNRAS subm. 30/01/2017. Based on the discovery results of H~Lin et al. (2017), ApJL subm. This paper has been approved for submission by DES, via Collaboration-Wide Review
Submitted: 2017-02-01
We present gravitational lens models of the multiply imaged quasar DES J0408-5354, recently discovered in the Dark Energy Survey (DES) footprint, with the aim of interpreting its remarkable quad-like configuration. We first model the DES single-epoch $grizY$ images as a superposition of a lens galaxy and four point-like objects, obtaining spectral energy distributions (SEDs) and relative positions for the objects. Three of the point sources (A,B,D) have SEDs compatible with the discovery quasar spectra, while the faintest point-like image (G2/C) shows significant reddening and a `grey' dimming of $\approx0.8$mag. In order to understand the lens configuration, we fit different models to the relative positions of A,B,D. Models with just a single deflector predict a fourth image at the location of G2/C but considerably brighter and bluer. The addition of a small satellite galaxy ($R_{\rm E}\approx0.2$") in the lens plane near the position of G2/C suppresses the flux of the fourth image and can explain both the reddening and grey dimming. All models predict a main deflector with Einstein radius between $1.7"$ and $2.0",$ velocity dispersion $267-280$km/s and enclosed mass $\approx 6\times10^{11}M_{\odot},$ even though higher resolution imaging data are needed to break residual degeneracies in model parameters. The longest time-delay (B-A) is estimated as $\approx 85$ (resp. $\approx125$) days by models with (resp. without) a perturber near G2/C. The configuration and predicted time-delays of J0408-5354 make it an excellent target for follow-up aimed at understanding the source quasar host galaxy and substructure in the lens, and measuring cosmological parameters. We also discuss some lessons learnt from J0408-5354 on lensed quasar finding strategies, due to its chromaticity and morphology.
[36]  [pdf] - 1668828
Quantitative Evaluation of Gender Bias in Astronomical Publications from Citation Counts
Comments: Abridged version to be submitted to Nature Astronomy. Comments welcome. For readers with very little time, the central result of the paper is covered by Figure 6 (Section 5)
Submitted: 2016-10-27
We analyze the role of first (leading) author gender on the number of citations that a paper receives, on the publishing frequency and on the self-citing tendency. We consider a complete sample of over 200,000 publications from 1950 to 2015 from five major astronomy journals. We determine the gender of the first author for over 70% of all publications. The fraction of papers which have a female first author has increased from less than 5% in the 1960s to about 25% today. We find that the increase of the fraction of papers authored by females is slowest in the most prestigious journals such as Science and Nature. Furthermore, female authors write 19$\pm$7% fewer papers in seven years following their first paper than their male colleagues. At all times papers with male first authors receive more citations than papers with female first authors. This difference has been decreasing with time and amounts to $\sim$6% measured over the last 30 years. To account for the fact that the properties of female and male first author papers differ intrinsically, we use a random forest algorithm to control for the non-gender specific properties of these papers which include seniority of the first author, number of references, total number of authors, year of publication, publication journal, field of study and region of the first author's institution. We show that papers authored by females receive 10.4$\pm$0.9% fewer citations than what would be expected if the papers with the same non-gender specific properties were written by the male authors. Finally, we also find that female authors in our sample tend to self-cite more, but that this effect disappears when controlled for non-gender specific variables.
[37]  [pdf] - 1453175
The mass-sheet degeneracy and time-delay cosmography: Analysis of the strong lens RXJ1131-1231
Comments: 25 (+8) pages, 12 figures, published in JCAP
Submitted: 2015-11-11, last modified: 2016-08-10
We present extended modeling of the strong lens system RXJ1131-1231 with archival data in two HST bands in combination with existing line-of-sight contribution and velocity dispersion estimates. Our focus is on source size and its influence on time-delay cosmography. We therefore examine the impact of mass-sheet degeneracy and especially the degeneracy pointed out by Schneider & Sluse (2013) using the source reconstruction scale. We also extend on previous work by further exploring the effects of priors on the kinematics of the lens and the external convergence in the environment of the lensing system. Our results coming from RXJ1131-1231 are given in a simple analytic form so that they can be easily combined with constraints coming from other cosmological probes. We find that the choice of priors on lens model parameters and source size are subdominant for the statistical errors for $H_0$ measurements of this systems. The choice of prior for the source is sub-dominant at present (2% uncertainty on $H_0$) but may be relevant for future studies. More importantly, we find that the priors on the kinematic anisotropy of the lens galaxy have a significant impact on our cosmological inference. When incorporating all the above modeling uncertainties, we find $H_0 = 86.6^{+6.8}_{-6.9}$ km s$^{-1}$Mpc$^{-1}$, when using kinematic priors similar to other studies. When we use a different kinematic prior motivated by Barnab\`e et al. (2012) but covering the same anisotropic range, we find $H_0 = 74.5^{+8.0}_{-7.8}$ km s$^{-1}$Mpc$^{-1}$. This means that the choice of kinematic modeling and priors have a significant impact on cosmographic inferences. The way forward is either to get better velocity dispersion measures which would down weight the impact of the priors or to construct physically motivated priors for the velocity dispersion model.
[38]  [pdf] - 1308186
Gravitational lens modeling with basis sets
Comments: 14 pages, 10 figures, accepted for publication in APJ
Submitted: 2015-04-28, last modified: 2015-10-21
We present a strong lensing modeling technique based on versatile basis sets for the lens and source planes. Our method uses high performance Monte Carlo algorithms, allows for an adaptive build up of complexity and bridges the gap between parametric and pixel based reconstruction methods. We apply our method to a HST image of the strong lens system RXJ1131-1231 and show that our method finds a reliable solution and is able to detect substructure in the lens and source planes simultaneously. Using mock data we show that our method is sensitive to sub-clumps with masses four orders of magnitude smaller than the main lens, which corresponds to about $10^8 M_{\odot}$, without prior knowledge on the position and mass of the sub-clump. The modelling approach is flexible and maximises automation to facilitate the analysis of the large number of strong lensing systems expected in upcoming wide field surveys. The resulting search for dark sub-clumps in these systems, without mass-to-light priors, offers promise for probing physics beyond the standard model in the dark matter sector.
[39]  [pdf] - 1202740
A simple model linking galaxy and dark matter evolution
Comments: 26 pages, 20 figures, 14 tables, accepted for publication in The Astrophysical Journal
Submitted: 2014-01-14, last modified: 2014-08-21
We construct a simple phenomenological model for the evolving galaxy population by incorporating pre-defined baryonic prescriptions into a dark matter hierarchical merger tree. Specifically the model is based on the simple gas-regulator model introduced by Lilly et al. 2013 coupled with the empirical quenching rules of Peng et al. 2010/12. The simplest model already does quite well in reproducing, without re-adjusting the input parameters, many observables including the Main Sequence sSFR-mass relation, the faint end slope of the galaxy mass function and the shape of the star-forming and passive mass functions. Compared with observations and/or the recent phenomenological model of Behroozi et al. 2013 based on epoch-dependent abundance-matching, our model also qualitatively reproduces the evolution of the Main Sequence sSFR(z) and SFRD(z) star formation rate density relations, the $M_s - M_h$ stellar-to-halo mass relation and also the $SFR - M_h$ relation. Quantitatively the evolution of sSFR(z) and SFRD(z) is not steep enough, the $M_s - M_h$ relation is not quite peaked enough and, surprisingly, the ratio of quenched to star-forming galaxies around M* is not quite high enough. We show that these deficiencies can simultaneously be solved by ad hoc allowing galaxies to re-ingest some of the gas previously expelled in winds, provided that this is done in a mass-dependent and epoch-dependent way. These allow the model galaxies to reduce an inherent tendency to saturate their star-formation efficiency. This emphasizes how efficient galaxies around M* are in converting baryons into stars and highlights the fact that quenching occurs just at the point when galaxies are rapidly approaching the maximum possible efficiency of converting baryons into stars.