Researcher profile

Tommaso Treu

Tommaso Treu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
34works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

34 published item(s)

preprint2026arXiv

BEACON: JWST NIRCam Pure-parallel Imaging Survey. II. Physical Properties of $z=7-14$ Galaxies

We present photometric properties of 161 galaxy candidates at $z=7-14$ selected from the second data release (DR2) of BEACON, a JWST Cycle 2 pure-parallel NIRCam imaging program. Carefully selected from 36 independent pointings (corresponding to $\sim350$\,arcmin$^2$ sky coverage), and hence with reduced cosmic variance, our galaxy candidates provide an unbiased sample for investigating galaxy properties over a wide range of environments. We measure the physical properties, including UV continuum slope ($β_{\rm UV}$), stellar mass ($M_*$), star formation rate (SFR), and sizes. Our highest redshift galaxy candidate at $z=13.71\pm0.15$ has a remarkably bright UV luminosity of $M_{\rm UV}=-21.19\pm0.08$, making it the brightest galaxy at $z>12$ if spectroscopically confirmed. With an extremely blue UV slope, compact morphology, and high star formation rate surface density ($Σ_{\rm SFR}$), this candidate may have extremely low metallicity, high ionizing photon escape fraction, or contributions from an AGN. Among our multiple independent sightlines, we identify three fields of galaxy number overdensity with $>3σ$ significance. The properties of galaxies in various environments do not exhibit significant differences, implying either that accelerated galaxy evolution in overdense regions is not yet widespread at $z>7$, or that the current constraints are limited by sample size. Our simulations indicate that increasing the sample by an order of magnitude would allow such environmental trends to be robustly confirmed or ruled out, underscoring the importance of future pure-parallel observations.

preprint2026arXiv

TDCOSMO XIX. Measuring stellar velocity dispersion with sub-percent accuracy for cosmography

The stellar velocity dispersion ($σ$) of massive elliptical galaxies is a key ingredient in breaking the mass-sheet degeneracy and obtaining precise and accurate cosmography from gravitational time delays. The relative uncertainty on the Hubble constant H$_0$ is double the relative error on $σ$. Therefore, time-delay cosmography imposes much more demanding requirements on the precision and accuracy of $σ$ than galaxy studies. While precision can be achieved with an adequate signal-to-noise ratio (S/N), the accuracy critically depends on key factors such as the elemental abundance and temperature of stellar templates, flux calibration, and wavelength ranges. We carried out a detailed study of the problem using multiple sets of galaxy spectra of massive elliptical galaxies with S/N$\sim$30--160 Å$^{-1}$, along with state-of-the-art empirical and semi-empirical stellar libraries and stellar population synthesis templates. We show that the choice of stellar library is generally the dominant source of residual systematic errors. We propose a general recipe for mitigating and accounting for residual uncertainties. We show that a sub-percent level of accuracy can be achieved on individual spectra with our data quality, which we subsequently validated with simulated mock datasets. The covariance between velocity dispersions measured for a sample of spectra can also be reduced to sub-percent levels. We recommend this recipe for all applications that require high precision and accurate stellar kinematics. Thus, we have made all the software publicly available to facilitate its implementation. This recipe will also be used in future TDCOSMO collaboration papers.

preprint2026arXiv

TDCOSMO XXI. Accurate stellar velocity dispersions of the SL2S lens sample and the fundamental plane of the lensing mass

We reanalyzed spectra that were taken as part of the SL2S lens galaxy survey with the goal to obtain the stellar velocity dispersion with a precision and accuracy sufficient for time-delay cosmography. In order to achieve this goal, we imposed stringent cuts on the signal-to-noise ratio (S/N), and employed recently developed methods to mitigate and quantify residual systematic errors that are transferred from template libraries and fitting process. We also quantified the covariance across the sample. For galaxy spectra with S/N $>20/$Å, our new measurements have an average random uncertainty of 3-4\%, an average systematic uncertainty of 2\%, and a covariance across the sample of 1\%. We find a negligible covariance between spectra taken with different instruments. The systematic uncertainty and covariance need to be included when the sample is used as an external dataset in time-delay cosmography. We revisited empirical scaling relations of lens galaxies based on the improved kinematics. We show that the SL2S sample, the TDCOSMO time-delay lens sample, and the lower-redshift SLACS sample follow the same correlation of the effective radius, stellar velocity dispersion, and lensing mass, known as the lensing-mass fundamental plane, as the previously derived correlation that assumed isothermal mass profiles for the deflectors. We also derived for the first time the lensing-mass fundamental plane assuming free power-law mass density profiles, and we show that the three samples also follow the same correlation. This is consistent with a scenario in which massive galaxies evolve by growing their radii and mass, but stay within the plane.

preprint2026arXiv

The BoRG-$JWST$ Survey: Abundance and Mass-to-light Ratio of Luminous $z=7-9$ Galaxies from Independent Sight Lines with NIRSpec

We present new results on the rest-frame UV luminosity function (UVLF) and stellar mass-to-light (M/L) ratio of bright (M$_{\rm UV}\lesssim-20$ mag) spectroscopically-confirmed galaxies at $z=7-9$ derived from the BoRG-$JWST$ survey, a unique data set of NIRSpec prism follow up of $HST$-selected sources from random-pointing imaging. By selecting galaxies from over 200 independent sight lines, the survey minimizes cosmic variance ensuring a statistically robust sample of the bright-galaxy population during the epoch of reionization. The data is used to constrain, for the first time, the bright end of the UVLF at $z=7-9$ from spectroscopically-confirmed galaxies over eight independent fields. We find that the bright end of the UVLF is higher than found using imaging over $JWST$ legacy fields, suggesting the latter may be significantly affected by cosmic variance, and thus reducing the tension with recent findings from $JWST$ at $z>10$ and comparable to models invoking little dust attenuation and bursty star formation. Additionally, we use the galaxies' $JWST$ spectra to infer their stellar masses and M/L ratios relative to other $HST$ and $JWST$ studies. We show that the stellar mass scales almost linearly with UV luminosity (M$_* \propto L_{\rm UV}^{0.85\pm0.12}$), albeit with large ($\sim0.5$ dex) intrinsic scatter, consistent with stochastic bursts of star formation in early galaxy formation.

preprint2025arXiv

TDCOSMO 2025: Cosmological constraints from strong lensing time delays

We present cosmological constraints from 8 strongly lensed quasars (hereafter, the TDCOSMO-2025 sample). Building on previous work, our analysis incorporated new deflector stellar velocity dispersions measured from spectra obtained with the James Webb Space Telescope (JWST), the Keck Telescopes, and the Very Large Telescope (VLT), utilizing improved methods. We used integrated JWST stellar kinematics for 5 lenses, VLT-MUSE for 2, and resolved kinematics from Keck and JWST for RXJ1131-1231. We also considered two samples of non-time-delay lenses: 11 from the Sloan Lens ACS (SLACS) sample with Keck-KCWI resolved kinematics; and 4 from the Strong Lenses in the Legacy Survey (SL2S) sample. We improved our analysis of line-of-sight effects, the surface brightness profile of the lens galaxies, and orbital anisotropy, and corrected for projection effects in the dynamics. Our uncertainties are maximally conservative by accounting for the mass-sheet degeneracy in the deflectors' mass density profiles. The analysis was blinded to prevent experimenter bias. Our primary result is based on the TDCOSMO-2025 sample, in combination with $Ω_{\rm m}$ constraints from the Pantheon+ Type Ia supernovae (SN) dataset. In the flat $Λ$ cold dark matter (CDM), we find $H_0=71.6^{+3.9}_{-3.3}$ km s$^{-1}$ Mpc$^{-1}$. The SLACS and SL2S samples are in excellent agreement with the TDCOSMO-2025 sample, improving the precision on $H_0$ in flat $Λ$CDM to 4.6%. Using the Dark Energy Survey SN Year-5 dataset (DES-SN5YR) or DESI-DR2 baryonic acoustic oscillations (BAO) likelihoods instead of Pantheon+ yields very similar results. We also present constraints in the open $Λ$CDM, $w$CDM, $w_0w_a$CDM, and $w_ϕ$CDM cosmologies. The TDCOSMO $H_0$ inference is robust and consistent across all presented cosmological models, and our cosmological constraints in them agree with those from the BAO and SN.

preprint2024arXiv

AGN STORM 2: V. Anomalous Behavior of the CIV Light Curve in Mrk 817

An intensive reverberation mapping campaign on the Seyfert 1 galaxy Mrk817 using the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST) revealed significant variations in the response of the broad UV emission lines to fluctuations in the continuum emission. The response of the prominent UV emission lines changes over a $\sim$60-day duration, resulting in distinctly different time lags in the various segments of the light curve over the 14 months observing campaign. One-dimensional echo-mapping models fit these variations if a slowly varying background is included for each emission line. These variations are more evident in the CIV light curve, which is the line least affected by intrinsic absorption in Mrk817 and least blended with neighboring emission lines. We identify five temporal windows with distinct emission line response, and measure their corresponding time delays, which range from 2 to 13 days. These temporal windows are plausibly linked to changes in the UV and X-ray obscuration occurring during these same intervals. The shortest time lags occur during periods with diminishing obscuration, whereas the longest lags occur during periods with rising obscuration. We propose that the obscuring outflow shields the ultraviolet broad lines from the ionizing continuum. The resulting change in the spectral energy distribution of the ionizing continuum, as seen by clouds at a range of distances from the nucleus, is responsible for the changes in the line response.

preprint2024arXiv

Early Results from GLASS-JWST XXIII: The transmission of Lyman-alpha from UV-faint z ~ 3-6 galaxies

Lyman-alpha (Ly$α$) emission from galaxies can be used to trace neutral hydrogen in the epoch of reionization, however, there is a degeneracy between the attenuation of Ly$α$ in the intergalactic medium (IGM) and the line profile emitted from the galaxy. Large shifts of Ly$α$ redward of systemic due to scattering in the interstellar medium can boost Ly$α$ transmission in the IGM during reionization. The relationship between Ly$α$ velocity offset from systemic and other galaxy properties is not well-established at high-redshift or low luminosities, due to the difficulty of observing emission lines which trace systemic redshift. Rest-frame optical spectroscopy with JWST/NIRSpec has opened a new window into understanding of Ly$α$ at z>3. We present a sample of 12 UV-faint galaxies ($-20 \lesssim$ MUV $\lesssim -16$) at $3 \lesssim z \lesssim 6$, with Ly$α$ velocity offsets, $Δv_{\mathrm{Ly}α}$, measured from VLT/MUSE and JWST/NIRSpec from the GLASS-JWST Early Release Program. We find median $Δv_{\mathrm{Ly}α}$ of 205 km s$^{-1}$ and standard deviation 75 km s$^{-1}$, compared to 320 and 170km s$^{-1}$ for MUV < -20 galaxies in the literature. Our new sample demonstrates the previously observed trend of decreasing Ly$α$ velocity offset with decreasing UV luminosity and optical line velocity dispersion, extends to MUV $\gtrsim$ -20, consistent with a picture where the Ly$α$ profile is shaped by gas close to the systemic redshift. Our results imply that during reionization Ly$α$ from UV-faint galaxies will be preferentially attenuated, but that detecting Ly$α$ with low $Δv_{\mathrm{Ly}α}$ can be an indicator of large ionized bubbles.

preprint2023arXiv

The brightest galaxies at Cosmic Dawn

Recent JWST observations suggest an excess of $z\gtrsim10$ galaxy candidates above most theoretical models. Here, we explore how the interplay between halo formation timescales, star formation efficiency and dust attenuation affects the properties and number densities of galaxies we can detect in the early universe. We calculate the theoretical upper limit on the UV luminosity function, assuming star formation is 100% efficient and all gas in halos is converted into stars, and that galaxies are at the peak age for UV emission (~10 Myr). This upper limit is ~4 orders of magnitude greater than current observations, implying these are fully consistent with star formation in $Λ$CDM cosmology. In a more realistic model, we use the distribution of halo formation timescales derived from extended Press-Schechter theory as a proxy for star formation rate (SFR). We predict that the galaxies observed so far at $z\gtrsim10$ are dominated by those with the fastest formation timescales, and thus most extreme SFRs and young ages. These galaxies can be upscattered by ~1.5 mag compared to the median UV magnitude vs halo mass relation. This likely introduces a selection effect at high redshift whereby only the youngest ($\lesssim$10 Myr), most highly star forming galaxies (specific SFR$\gtrsim$30 Gyr$^{-1}$) have been detected so far. Furthermore, our modelling suggests that redshift evolution at the bright end of the UV luminosity function is substantially affected by the build-up of dust attenuation. We predict that deeper JWST observations (reaching m~30) will reveal more typical galaxies with relatively older ages (~100 Myr) and less extreme specific SFRs (~10 Gyr$^{-1}$ for a $M_\mathrm{UV}$ ~ -20 galaxy at z~10).

preprint2022arXiv

Concordance between observations and simulations in the evolution of the mass relation between supermassive black holes and their host galaxies

We carry out a comparative analysis of the relation between the mass of supermassive black holes (BHs) and the stellar mass of their host galaxies at $0.2<z<1.7$ using well-matched observations and multiple state-of-the-art simulations (e.g., Massive Black II, Horizon-AGN, Illustris, TNG and a semi-analytic model). The observed sample consists of 646 uniformly-selected SDSS quasars ($0.2 < z < 0.8$) and 32 broad-line active galactic nuclei (AGNs; $1.2<z<1.7$) with imaging from Hyper Suprime-Cam (HSC) for the former and Hubble Space Telescope (HST) for the latter. We first add realistic observational uncertainties to the simulation data and then construct a simulated sample in the same manner as the observations. Over the full redshift range, our analysis demonstrates that all simulations predict a level of intrinsic scatter of the scaling relations comparable to the observations which appear to agree with the dispersion of the local relation. Regarding the mean relation, Horizon-AGN and TNG are in closest agreement with the observations at low and high redshift ($z\sim$ 0.2 and 1.5, respectively) while the other simulations show subtle differences within the uncertainties. For insight into the physics involved, the scatter of the scaling relation, seen in the SAM, is reduced by a factor of two and closer to the observations after adopting a new feedback model that considers the geometry of the AGN outflow. The consistency in the dispersion with redshift in our analysis supports the importance of both quasar- and radio-mode feedback prescriptions in the simulations. Finally, we highlight the importance of increasing the sensitivity (e.g., using the James Webb Space Telescope), thereby pushing to lower masses and minimizing biases due to selection effects.

preprint2022arXiv

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

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

preprint2022arXiv

Early results from GLASS-JWST. IX: First spectroscopic confirmation of low-mass quiescent galaxies at $z>2$ with NIRISS

How passive galaxies form, and the physical mechanisms which prevent star formation over long timescales, are some of the most outstanding questions in understanding galaxy evolution. The properties of quiescent galaxies over cosmic time provide crucial information to identify the quenching mechanisms. Passive galaxies have been confirmed and studied out to $z\sim4$, but all of these studies have been limited to massive systems (mostly with $\log{(M_{\rm star}/M_{\odot})}>10.8$). Using James Webb Space Telescope (JWST) NIRISS grism slitless spectroscopic data from the GLASS JWST ERS program, we present spectroscopic confirmation of two quiescent galaxies at $z_{\rm spec}=2.650^{+0.004}_{-0.006}$ and $z_{\rm spec}=2.433^{+0.032}_{-0.016}$ (3$σ$ errors) with stellar masses of $\log{(M_{\rm star}/M_{\odot})}=10.53^{+0.18}_{-0.06}$ and $\log{(M_{\rm star}/M_{\odot})}=9.93^{+0.06}_{-0.07}$ (corrected for magnification factors of $μ=2.0$ and $μ=2.1$, respectively). The latter represents the first spectroscopic confirmation of the existence of low-mass quiescent galaxies at cosmic noon, showcasing the power of JWST to identify and characterize this enigmatic population.

preprint2022arXiv

Identifying lensed quasars and measuring their time-delays from unresolved light curves

Identifying multiply imaged quasars is challenging due to their low density in the sky and the limited angular resolution of wide field surveys. We show that multiply imaged quasars can be identified using unresolved light curves, without assuming a light curve template or any prior information. After describing our method, we show using simulations that it can attain high precision and recall when we consider high-quality data with negligible noise well below the variability of the light curves. As the noise level increases to that of the Zwicky Transient Facility (ZTF) telescope, we find that precision can remain close to $100\%$ while recall drops to $\sim 60\%$. We also consider some examples from the Time Delay Challenge 1 (TDC1) and demonstrate that the time delays can be accurately recovered from the joint light curve data in realistic observational scenarios. We further demonstrate our method by applying it to publicly available COSMOGRAIL data of the observed lensed quasar SDSS J1226-0006. We identify the system as a lensed quasar based on the unresolved light curve and estimate a time delay in good agreement with the one measured by COSMOGRAIL using the individual image light curves. The technique shows great potential to identify lensed quasars in wide field imaging surveys, especially the soon to be commissioned Vera Rubin Observatory.

preprint2022arXiv

RELICS: Small Lensed $z\geq5.5$ Galaxies Selected as Potential Lyman Continuum Leakers

We present size measurements of 78 high-redshift ($z\geq 5.5$) galaxy candidates from the Reionisation Lensing Cluster Survey (RELICS). These distant galaxies are well-resolved due to the gravitational lensing power of foreground galaxy clusters, imaged by the Hubble Space Telescope (HST) and the Spitzer Space Telescope. We compute sizes using the forward-modeling code Lenstruction and account for magnification using public lens models. The resulting size-magnitude measurements confirm the existence of many small galaxies with effective radii $R_{\rm{eff}}<200$ pc in the early universe, in agreement with previous studies. In addition, we highlight compact and highly star-forming sources with star formation rate surface densities $Σ_\text{SFR}>10M_\odot\text{yr}^{-1}\text{kpc}^{-2}$ as possible Lyman continuum leaking candidates that could be major contributors to the process of reionisation. Future spectroscopic follow-up of these compact galaxies (e.g., with the James Webb Space Telescope) will further clarify their role in reionisation and the physics of early star formation.

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Observational Facilities to Study Dark Matter

We present an overview of future observational facilities that will significantly enhance our understanding of the fundamental nature of dark matter. These facilities span a range of observational techniques including optical/near-infrared imaging and spectroscopy, measurements of the cosmic microwave background, pulsar timing, 21-cm observations of neutral hydrogen at high redshift, and the measurement of gravitational waves. Such facilities are a critical component of a multi-pronged experimental program to uncover the nature of dark matter, while often providing complementary measurements of dark energy, neutrino physics, and inflation.

preprint2022arXiv

The Broad Line Region and Black Hole Mass of NGC 4151

We present a reanalysis of reverberation-mapping data from 2005 for the Seyfert galaxy NGC 4151, supplemented with additional data from the literature to constrain the continuum variations over a significantly longer baseline than the original monitoring program. Modeling of the continuum light curve and the velocity-resolved variations across the H$β$ emission line constrains the geometry and kinematics of the broad line region (BLR). The BLR is well described by a very thick disk with similar opening angle ($θ_o \approx 57^{\circ}$) and inclination angle ($θ_i \approx 58^{\circ}$), suggesting that our sight line towards the innermost central engine skims just above the surface of the BLR. The inclination is consistent with constraints from geometric modeling of the narrow line region, and the similarity between the inclination and opening angles is intriguing given previous studies of NGC 4151 that suggest BLR gas has been observed temporarily eclipsing the X-ray source. The BLR kinematics are dominated by eccentric bound orbits, with $\sim10$% of the orbits preferring near-circular motions. With the BLR geometry and kinematics constrained, the models provide an independent and direct black hole mass measurement of $\log M_{\rm BH}/M_{\odot} = 7.22^{+0.11}_{-0.10}$ or $M_{\rm BH}=1.66^{+0.48}_{-0.34}\times10^7 M_{\odot}$, which is in good agreement with mass measurements from stellar dynamical modeling and gas dynamical modeling. NGC 4151 is one of the few nearby broad-lined Seyferts where the black hole mass may be measured via multiple independent techniques, and it provides an important test case for investigating potential systematics that could affect the black hole mass scales used in the local Universe and for high-redshift quasars.

preprint2022arXiv

The Lick AGN Monitoring Project 2016: Dynamical Modeling of Velocity-Resolved H\b{eta} Lags in Luminous Seyfert Galaxies

We have modeled the velocity-resolved reverberation response of the H\b{eta} broad emission line in nine Seyfert 1 galaxies from the Lick Active Galactic Nucleus (AGN) Monitioring Project 2016 sample, drawing inferences on the geometry and structure of the low-ionization broad-line region (BLR) and the mass of the central supermassive black hole. Overall, we find that the H\b{eta} BLR is generally a thick disk viewed at low to moderate inclination angles. We combine our sample with prior studies and investigate line-profile shape dependence, such as log10(FWHM/σ), on BLR structure and kinematics and search for any BLR luminosity-dependent trends. We find marginal evidence for an anticorrelation between the profile shape of the broad H\b{eta} emission line and the Eddington ratio, when using the root-mean-square spectrum. However, we do not find any luminosity-dependent trends, and conclude that AGNs have diverse BLR structure and kinematics, consistent with the hypothesis of transient AGN/BLR conditions rather than systematic trends.

preprint2022arXiv

The primordial matter power spectrum on sub-galactic scales

The primordial matter power spectrum quantifies fluctuations in the distribution of dark matter immediately following inflation. Over cosmic time, over-dense regions of the primordial density field grow and collapse into dark matter halos, whose abundance and density profiles retain memory of the initial conditions. By analyzing the image magnifications in eleven strongly-lensed and quadruply-imaged quasars, we infer the abundance and concentrations of low-mass halos, and cast the measurement in terms of the amplitude of the primordial matter power spectrum. We anchor the power spectrum on large scales, isolating the effect of small-scale deviations from the $Λ$CDM prediction. Assuming an analytic model for the power spectrum and accounting for several sources of potential systematic uncertainty, including three different models for the halo mass function, we obtain correlated inferences of $\log_{10}\left(P / P_{Λ\rm{CDM}}\right)$, the power spectrum amplitude relative to the predictions of the concordance cosmological model, of $0.0_{-0.4}^{+0.5}$, $0.1_{-0.6}^{+0.7}$, and $0.2_{-0.9}^{+1.0}$ at k = 10, 25 and 50 $\rm{Mpc^{-1}}$ at $68 \%$ confidence, consistent with cold dark matter and single-field slow-roll inflation.

preprint2022arXiv

The size-luminosity relation of lensed galaxies at $z=6-9$ in the Hubble Frontier Fields

We measure the size-luminosity relation of photometrically-selected galaxies within the redshift range $z\sim6-9$, using galaxies lensed by six foreground Hubble Frontier Fields (HFF) clusters. The power afforded by strong gravitational lensing allows us to observe fainter and smaller galaxies than in blank fields. We select our sample of galaxies and obtain their properties, e.g., redshift, magnitude, from the photometrically-derived ASTRODEEP catalogues. The intrinsic size is measured with the Lenstruction software, and completeness maps are created as a function of size and luminosity via the GLACiAR2 software. We perform a Bayesian analysis to estimate the intrinsic and incompleteness-corrected size-luminosity distribution, with parameterization $r_e \propto L^β$. We find slopes of $β\sim0.48\pm0.08$ at $z\sim6-7$ and $β\sim0.68\pm0.14$ at $z\sim8.5$, adopting the Bradac lens model. The slope derived by lensed galaxies is steeper than that obtained in blank fields and is consistent with other independent determinations of the size-luminosity relation from the HFF dataset. We also investigate the systematic uncertainties correlated with the choice of lens models, finding that the slopes of size-luminosity relations derived from different models are consistent with each other, i.e. the modeling errors are not a significant source of uncertainty in the size-luminosity relation.

preprint2021arXiv

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

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

preprint2021arXiv

The Lick AGN Monitoring Project 2016: Velocity-Resolved Hβ Lags in Luminous Seyfert Galaxies

We carried out spectroscopic monitoring of 21 low-redshift Seyfert 1 galaxies using the Kast double spectrograph on the 3-m Shane telescope at Lick Observatory from April 2016 to May 2017. Targeting active galactic nuclei (AGN) with luminosities of λLλ (5100 Å) = 10^44 erg/s and predicted Hβ lags of 20-30 days or black hole masses of 10^7-10^8.5 Msun, our campaign probes luminosity-dependent trends in broad-line region (BLR) structure and dynamics as well as to improve calibrations for single-epoch estimates of quasar black hole masses. Here we present the first results from the campaign, including Hβ emission-line light curves, integrated Hβ lag times (8-30 days) measured against V-band continuum light curves, velocity-resolved reverberation lags, line widths of the broad Hβ components, and virial black hole mass estimates (10^7.1-10^8.1 Msun). Our results add significantly to the number of existing velocity-resolved lag measurements and reveal a diversity of BLR gas kinematics at moderately high AGN luminosities. AGN continuum luminosity appears not to be correlated with the type of kinematics that its BLR gas may exhibit. Follow-up direct modeling of this dataset will elucidate the detailed kinematics and provide robust dynamical black hole masses for several objects in this sample.

preprint2021arXiv

The Physical Properties of Luminous $z\gtrsim8$ Galaxies and Implications for the Cosmic Star Formation Rate Density From $\sim$0.35 deg$^{2}$ of (Pure-)Parallel HST Observations

We present the largest systematic, HST-based search to date for luminous $z\gtrsim8$ galaxy candidates using $\sim$1267 arcmin$^{2}$ of (pure-)parallel observations from a compilation of 288 random sightlines with ACS and WFC3 observations, derived from the SuperBoRG data set and together representing a factor $\sim1.12\times$ larger than existing space-based data sets. Using NIR color cuts and careful photo-$z$ analyses, we find 31 $z\gtrsim8$ galaxy candidates over 29 unique sightlines, and derive global galaxy properties such as UV magnitudes and continuum slopes, sizes, and rest-frame optical properties (e.g., SFRs, stellar masses, $A_{\rm v}$). Taking advantage of the (pure-)parallel nature of our data set - making it one of the most representative thus far - and derived SFRs, we evaluate the cosmic star formation rate density for the bright end of the UV luminosity function at $z\sim8-10$ and test the validity of luminosity function-derived results using a conversion factor. We find our method yields comparable results to those derived with luminosity functions. Furthermore, we present follow up observations of 5 (Super)BoRG targets with Keck/MOSFIRE, finding no evidence of Ly$α$ in $>3$ hrs of $Y-$band observations in either, consistent with a largely neutral medium at $z\sim8$. Our results offer a definitive HST legacy on the bright end of the luminosity function and provide a valuable benchmark as well as targets for follow up with JWST.

preprint2020arXiv

A Census of Sub-kiloparsec Resolution Metallicity Gradients in Star-forming Galaxies at Cosmic Noon from HST Slitless Spectroscopy

We present hitherto the largest sample of gas-phase metallicity radial gradients measured at sub-kiloparsec resolution in star-forming galaxies in the redshift range of $z\in[1.2, 2.3]$. These measurements are enabled by the synergy of slitless spectroscopy from the Hubble Space Telescope near-infrared channels and the lensing magnification from foreground galaxy clusters. Our sample consists of 76 galaxies with stellar mass ranging from 10$^7$ to 10$^{10}$ $M_\odot$, instantaneous star-formation rate in the range of [1, 100] $M_\odot$/yr, and global metallicity [$\frac{1}{12}$, 2] solar. At 2-$σ$ confidence level, 15/76 galaxies in our sample show negative radial gradients, whereas 7/76 show inverted gradients. Combining ours and all other metallicity gradients obtained at similar resolution currently available in the literature, we measure a negative mass dependence of $Δ\log({\rm O/H})/Δr~ [\mathrm{dex~kpc^{-1}}] = \left(-0.020\pm0.007\right) + \left(-0.016\pm0.008\right) \log(M_\ast/10^{9.4} M_\odot)$ with the intrinsic scatter being $σ=0.060\pm0.006$ over four orders of magnitude in stellar mass. Our result is consistent with strong feedback, not secular processes, being the primary governor of the chemo-structural evolution of star-forming galaxies during the disk mass assembly at cosmic noon. We also find that the intrinsic scatter of metallicity gradients increases with decreasing stellar mass and increasing specific star-formation rate. This increase in the intrinsic scatter is likely caused by the combined effect of cold-mode gas accretion and merger-induced starbursts, with the latter more predominant in the dwarf mass regime of $M_\ast\lesssim10^9 M_\odot$.

preprint2020arXiv

Cosmic dissonance: new physics or systematics behind a short sound horizon?

Persistent tension between low-redshift observations and the Cosmic Microwave Background radiation (CMB), in terms of two fundamental distance scales set by the sound horizon $r_d$ and the Hubble constant $H_0$, suggests new physics beyond the Standard Model or residual systematics. We examine recently updated distance calibrations from Cepheids, gravitational lensing time-delay observations, and the Tip of the Red Giant Branch. Calibrating the Baryon Acoustic Oscillations (BAO) and Type Ia supernovae with combinations of the distance indicators, we obtain a joint and self-consistent measurement of $H_0$ and $r_d$ at low redshift, independent of cosmological models and CMB inference. In an attempt to alleviate the tension between late-time and CMB-based measurements, we consider four extensions of the standard $Λ$CDM model. The sound horizon from our different measurements is $r_d=(137\pm3^{stat.}\pm2^{syst.})$~Mpc. Depending on the adopted distance indicators, the $combined$ tension in $H_0$ and $r_d$ ranges between 2.3 and 5.1 $σ$. We find that modifications of $Λ$CDM that change the physics after recombination fail to solve the problem, for the reason that they only resolve the tension in $H_0$, while the tension in $r_d$ remains unchanged. Pre-recombination extensions (with early dark energy or the effective number of neutrinos $\rm{N}_{\rm{eff}}=3.24 \pm 0.16$) are allowed by the data, unless the calibration from Cepheids is included. Results from time-delay lenses are consistent with those from distance-ladder calibrations and point to a discrepancy between absolute distance scales measured from the CMB (assuming the standard cosmological model) and late-time observations. New proposals to resolve this tension should be examined with respect to reconciling not only the Hubble constant but also the sound horizon derived from the CMB and other cosmological probes.

preprint2020arXiv

Robotic reverberation mapping of the broad-line radio galaxy 3C 120

We carried out photometric and spectroscopic observations of the well-studied broad-line radio galaxy 3C 120 with the Las Cumbres Observatory (LCO) global robotic telescope network from 2016 December to 2018 April as part of the LCO AGN Key Project on Reverberation Mapping of Accretion Flows. Here, we present both spectroscopic and photometric reverberation mapping results. We used the interpolated cross-correlation function (ICCF) to perform multiple-line lag measurements in 3C 120. We find the H$γ$, He II $λ4686$, H$β$ and He I $λ5876$ lags of $τ_{\text{cen}} = 18.8_{-1.0}^{+1.3}$, $2.7_{-0.8}^{+0.7}$, $21.2_{-1.0}^{+1.6}$, and $16.9_{-1.1}^{+0.9}$ days respectively, relative to the V-band continuum. Using the measured lag and rms velocity width of the H$β$ emission line, we determine the mass of the black hole for 3C 120 to be $M=\left(6.3^{+0.5}_{-0.3}\right)\times10^7\,(f/5.5)$ M$_\odot$. Our black hole mass measurement is consistent with similar previous studies on 3C 120, but with small uncertainties. In addition, velocity-resolved lags in 3C 120 show a symmetric pattern across the H$β$ line, 25 days at line centre decreasing to 17 days in the line wings at $\pm4000$ km s$^{-1}$. We also investigate the inter-band continuum lags in 3C 120 and find that they are generally consistent with $τ\proptoλ^{4/3}$ as predicted from a geometrically-thin, optically-thick accretion disc. From the continuum lags, we measure the best fit value $τ_{\rm 0} = 3.5\pm 0.2$ days at $λ_{\rm 0} = 5477$A. It implies a disc size a factor of $1.6$ times larger than prediction from the standard disc model with $L/L_{\rm Edd} = 0.4$. This is consistent with previous studies in which larger than expected disc sizes were measured.

preprint2020arXiv

Spectroscopically Confirmed Lyman-Alpha Emitters from Redshift 5 to 7 Behind Ten Galaxy Cluster Lenses

We present 36 spectroscopically confirmed intrinsically UV-faint Ly$α$ emitting galaxies from follow-up observations with Keck/DEIMOS of gravitationally lensed high-redshift candidates. Candidates were selected to be between $5\lesssim z \lesssim 7$ from photometric data using \textit{HST} and \textit{Spitzer} imaging surveys. We used photometric redshift information to perform an integrated photometric redshift probability cut $>1\%$ between $5<z<7$ to construct a sample of 198 high-redshift objects. Our high-redshift sample spans intrinsic UV luminosities from a few $L^*$ down to $0.001L^*$. We identified 19 high-confidence detections of Ly$α$ and an additional 17 likely detections. We divided our sample into lower-redshift ($z\sim5.5$) and higher-redshift ($z\sim6.5$) bins and ran Monte Carlo trials, incorporating the strength of the Ly$α$ emission and the photometric redshift of the non-detections. Considering only objects where Ly$α$ could be detected at EW(Ly$α$)$>$25Å at $3σ$ at the fiducial depth of our survey, and only those galaxies with EW(Ly$α$)$>$25Å as true LAEs, and finally, only objects with $m_{AB}<26.8$, we found the LAE fraction to be flat, or modestly increase from 0.26$\pm0.04$ to 0.30$\pm0.04$. These values relative to those for lower-redshift samples are consistent with a rising LAE fraction with redshift out to $z\sim6$, but at $z\sim6.5$ there is some tension between our results and results from surveys at intrinsically brighter luminosities. We conclude intrinsically fainter galaxies have Ly$α$ emission, and there is a steep drop in the LAE fraction from our high-redshift sample at $z\sim6.5$ and from similar galaxies at $z\sim7.5$. This likely indicates we are witnessing the tail end of the epoch of reionization, as such a drop is not expected due to changes of intrinsic galaxy properties between these redshifts.

preprint2020arXiv

The Impact of Scatter in the Galaxy UV Luminosity to Halo Mass Relation on Ly$α$ Visibility During the Epoch of Reionization

The reionization of hydrogen is closely linked to the first structures in the universe, so understanding the timeline of reionization promises to shed light on the nature of these early objects. In particular, transmission of Lyman alpha (Ly$α$) from galaxies through the intergalactic medium (IGM) is sensitive to neutral hydrogen in the IGM, so can be used to probe the reionization timeline. In this work, we implement an improved model of the galaxy UV luminosity to dark matter halo mass relation to infer the volume-averaged fraction of neutral hydrogen in the IGM from Ly$α$ observations. Many models assume that UV-bright galaxies are hosted by massive dark matter haloes in overdense regions of the IGM, so reside in relatively large ionized regions. However, observations and N-body simulations indicate that scatter in the UV luminosity-halo mass relation is expected. Here, we model the scatter (though we assume the IGM topology is unaffected) and assess the impact on Ly$α$ visibility during reionization. We show that UV luminosity-halo mass scatter reduces Ly$α$ visibility compared to models without scatter, and that this is most significant for UV-bright galaxies. We then use our model with scatter to infer the neutral fraction, $\overline{x}_{\mathrm{HI}}$, at $z \sim 7$ using a sample of Lyman-break galaxies in legacy fields. We infer $\overline{x}_{\mathrm{HI}} = 0.55_{-0.13}^{+0.11}$ with scatter, compared to $\overline{x}_{\mathrm{HI}} = 0.59_{-0.14}^{+0.12}$ without scatter, a very slight decrease and consistent within the uncertainties. Finally, we place our results in the context of other constraints on the reionization timeline and discuss implications for future high-redshift galaxy studies.

preprint2020arXiv

Time Delay Lens Modeling Challenge: I. Experimental Design

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

preprint2019arXiv

Constraints on the mass-concentration relation of cold dark matter halos with 11 strong gravitational lenses

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.

preprint2019arXiv

Evolution of the Stellar Mass-Metallicity Relation. II. Constraints on Galactic Outflows from the Mg Abundances of Quiescent Galaxies

We present the stellar mass-[Fe/H] and mass-[Mg/H] relation of quiescent galaxies in two galaxy clusters at $z\sim0.39$ and $z\sim0.54$. We derive the age, [Fe/H], and [Mg/Fe] for each individual galaxy using a full-spectrum fitting technique. By comparing with the relations for $z\sim0$ SDSS galaxies, we confirm our previous finding that the mass-[Fe/H] relation evolves with redshift. The mass-[Fe/H] relation at higher redshift has lower normalization and possibly steeper slope. However, based on our sample, the mass-[Mg/H] relation does not evolve over the observed redshift range. We use a simple analytic chemical evolution model to constrain average outflow that these galaxies experience over their lifetime, via the calculation of mass-loading factor. We find that the average mass-loading factor $η$ is a power-law function of galaxy stellar mass, $η\propto M_*^{-0.21\pm0.09}$. The measured mass-loading factors are consistent with the results of other observational methods for outflow measurements and with the predictions where outflow is caused by star formation feedback in turbulent disks.

preprint2019arXiv

H0LiCOW XII. Lens mass model of WFI2033-4723 and blind measurement of its time-delay distance and $H_0$

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 $Λ$CDM cosmology with a uniform prior on $Ω_\mathrm{m}$ in the range [0.05, 0.5]. This work is part of the $H_0$ Lenses in COSMOGRAIL&#39;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).

preprint2019arXiv

H0LiCOW XIII. A 2.4% measurement of $H_{0}$ from lensed quasars: $5.3σ$ tension between early and late-Universe probes

We present a measurement of the Hubble constant ($H_{0}$) and other cosmological parameters from a joint analysis of six gravitationally lensed quasars with measured time delays. All lenses except the first are analyzed blindly with respect to the cosmological parameters. In a flat $Λ$CDM cosmology, we find $H_{0} = 73.3_{-1.8}^{+1.7}$, a 2.4% precision measurement, in agreement with local measurements of $H_{0}$ from type Ia supernovae calibrated by the distance ladder, but in $3.1σ$ tension with $Planck$ observations of the cosmic microwave background (CMB). This method is completely independent of both the supernovae and CMB analyses. A combination of time-delay cosmography and the distance ladder results is in $5.3σ$ tension with $Planck$ CMB determinations of $H_{0}$ in flat $Λ$CDM. We compute Bayes factors to verify that all lenses give statistically consistent results, showing that we are not underestimating our uncertainties and are able to control our systematics. We explore extensions to flat $Λ$CDM using constraints from time-delay cosmography alone, as well as combinations with other cosmological probes, including CMB observations from $Planck$, baryon acoustic oscillations, and type Ia supernovae. Time-delay cosmography improves the precision of the other probes, demonstrating the strong complementarity. Allowing for spatial curvature does not resolve the tension with $Planck$. Using the distance constraints from time-delay cosmography to anchor the type Ia supernova distance scale, we reduce the sensitivity of our $H_0$ inference to cosmological model assumptions. For six different cosmological models, our combined inference on $H_{0}$ ranges from $\sim73$-$78~\mathrm{km~s^{-1}~Mpc^{-1}}$, which is consistent with the local distance ladder constraints.

preprint2019arXiv

The mass relations between supermassive black holes and their host galaxies at 1<z<2 with HST-WFC3

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.

preprint2019arXiv

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

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.

preprint2019arXiv

Where do quasar hosts lie with respect to the size-mass relation of galaxies?

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.