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] oai:arXiv.org:2003.12117 [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] oai:arXiv.org:2003.03277 [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] oai:arXiv.org:2001.07719 [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] oai:arXiv.org:1912.09133 [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] oai:arXiv.org:1912.08027 [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] oai:arXiv.org:1908.06344 [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] oai:arXiv.org:1905.09338 [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] oai:arXiv.org:1910.11875 [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] oai:arXiv.org:1908.06983 [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

Comments:

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] oai:arXiv.org:1907.04869 [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] oai:arXiv.org:1911.02479 [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 (https://petabytestoscience.github.io/).

[12] oai:arXiv.org:1910.14242 [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] oai:arXiv.org:1910.06306 [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] oai:arXiv.org:1909.07986 [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] oai:arXiv.org:1909.02573 [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] oai:arXiv.org:1905.08800 [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] oai:arXiv.org:1904.10965 [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:
https://github.com/sibirrer/lensing-astrometry

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] oai:arXiv.org:1905.12496 [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] oai:arXiv.org:1901.11031 [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] oai:arXiv.org:1907.02533 [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] oai:arXiv.org:1809.01274 [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:
http://shsuyu.github.io/H0LiCOW/site/notebooks/H0_from_lenses.html, 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] oai:arXiv.org:1902.01055 [pdf] - 1872765

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

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
https://lsstdarkmatter.github.io/ .

[23] oai:arXiv.org:1903.04742 [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] oai:arXiv.org:1903.05035 [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] oai:arXiv.org:1903.04425 [pdf] - 1846071

Dark Matter Science in the Era of LSST

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

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] oai:arXiv.org:1807.09278 [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] oai:arXiv.org:1803.09746 [pdf] - 1783742

Lenstronomy: multi-purpose gravitational lens modelling software package

Comments: software and documentation see https://lenstronomy.readthedocs.io,
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 https://lenstronomy.readthedocs.io.

[28] oai:arXiv.org:1809.09121 [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
https://github.com/herjy/SLIT.

[29] oai:arXiv.org:1808.04838 [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] oai:arXiv.org:1712.04945 [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] oai:arXiv.org:1710.01303 [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] oai:arXiv.org:1703.08410 [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] oai:arXiv.org:1702.00009 [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] oai:arXiv.org:1610.01599 [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] oai:arXiv.org:1702.00406 [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] oai:arXiv.org:1610.08984 [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] oai:arXiv.org:1511.03662 [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] oai:arXiv.org:1504.07629 [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] oai:arXiv.org:1401.3162 [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.