Source author record

Tommaso Treu

Tommaso Treu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

136works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

136 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'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 "Time Delay Lens Modeling Challenge" (TDLMC). The goal of this challenge is to assess the present capabilities of lens modeling codes and assumptions and test the level of accuracy of inferred cosmological parameters given realistic mock datasets. We invite scientists to model a set of simulated HST observations of 50 mock lens systems. The systems are organized in rungs, with the complexity and realism increasing going up the ladder. The goal of the challenge is to infer $H_0$ for each rung, given the HST images, the time delay, and stellar velocity dispersion of the deflector for a fixed background cosmology. The TDLMC challenge starts with the mock data release on 2018 January 8th. The deadline for blind submission is different for each rung. The deadline for Rung 0-1 is 2018 September 8; the deadline for Rung 2 is 2019 April 8 and the one for Rung 3 is 2019 September 8. This first paper gives an overview of the challenge including the data design, and a set of metrics to quantify the modeling performance and challenge details.

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'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.

preprint2016arXiv

Broad Hbeta Emission-Line Variability in a Sample of 102 Local Active Galaxies

A sample of 102 local (0.02 < z < 0.1) Seyfert galaxies with black hole masses MBH > 10^7 M_sun was selected from the Sloan Digital Sky Survey (SDSS) and observed using the Keck 10-m telescope to study the scaling relations between MBH and host galaxy properties. We study profile changes of the broad Hbeta emission line within the ~3-9 year time-frame between the two sets of spectra. The variability of the broad Hbeta emission line is of particular interest, not only since it is used to estimate MBH, but also since its strength and width is used to classify Seyfert galaxies into different types. At least some form of broad-line variability (in either width or flux) is observed in the majority (~66%) of the objects, resulting in a Seyfert-type change for ~38% of the objects, likely driven by variable accretion and/or obscuration. The broad Hbeta line virtually disappears in 3/102 (~3%) extreme cases. We discuss potential causes for these changing-look AGNs. While similar dramatic transitions have previously been reported in the literature, either on a case-by-case basis or in larger samples focusing on quasars at higher redshifts, our study provides statistical information on the frequency of H$β$ line variability in a sample of low-redshift Seyfert galaxies.

preprint2016arXiv

Detection of Lyman-Alpha Emission From a Triple Imaged z=6.85 Galaxy Behind MACS J2129.4-0741

We report the detection of Ly$α$ emission at $\sim9538$Å in the Keck/DEIMOS and \HST WFC3 G102 grism data from a triply-imaged galaxy at $z=6.846\pm0.001$ behind galaxy cluster MACS J2129.4$-$0741. Combining the emission line wavelength with broadband photometry, line ratio upper limits, and lens modeling, we rule out the scenario that this emission line is \oii at $z=1.57$. After accounting for magnification, we calculate the weighted average of the intrinsic Ly$α$ luminosity to be $\sim1.3\times10^{42}~\mathrm{erg}~\mathrm{s}^{-1}$ and Ly$α$ equivalent width to be $74\pm15$Å. Its intrinsic UV absolute magnitude at 1600Å is $-18.6\pm0.2$ mag and stellar mass $(1.5\pm0.3)\times10^{7}~M_{\odot}$, making it one of the faintest (intrinsic $L_{UV}\sim0.14~L_{UV}^*$) galaxies with Ly$α$ detection at $z\sim7$ to date. Its stellar mass is in the typical range for the galaxies thought to dominate the reionization photon budget at $z\gtrsim7$; the inferred Ly$α$ escape fraction is high ($\gtrsim 10$\%), which could be common for sub-$L^*$ $z\gtrsim7$ galaxies with Ly$α$ emission. This galaxy offers a glimpse of the galaxy population that is thought to drive reionization, and it shows that gravitational lensing is an important avenue to probe the sub-$L^*$ galaxy population.

preprint2016arXiv

Spitzer UltRa Faint SUrvey Program (SURFS UP). II. IRAC-Detected Lyman-Break Galaxies at 6 < z < 10 Behind Strong-Lensing Clusters

We study the stellar population properties of the IRAC-detected $6 \lesssim z \lesssim 10$ galaxy candidates from the Spitzer UltRa Faint SUrvey Program (SURFS UP). Using the Lyman Break selection technique, we find a total of 16 new galaxy candidates at $6 \lesssim z \lesssim 10$ with $S/N \geq 3$ in at least one of the IRAC $3.6μ$m and $4.5μ$m bands. According to the best mass models available for the surveyed galaxy clusters, these IRAC-detected galaxy candidates are magnified by factors of $\sim 1.2$--$5.5$. We find that the IRAC-detected $6 \lesssim z \lesssim 10$ sample is likely not a homogeneous galaxy population: some are relatively massive (stellar mass as high as $4 \times 10^9\,M_{\odot}$) and evolved (age $\lesssim 500$ Myr) galaxies, while others are less massive ($M_{\text{stellar}}\sim 10^8\,M_{\odot}$) and very young ($\sim 10$ Myr) galaxies with strong nebular emission lines that boost their rest-frame optical fluxes. We identify two Ly$α$ emitters in our sample from the Keck DEIMOS spectra, one at $z_{\text{Ly}α}=6.76$ (in RXJ1347) and one at $z_{\text{Ly}α}=6.32$ (in MACS0454). We show that IRAC $[3.6]-[4.5]$ color, when combined with photometric redshift, can be used to identify galaxies likely with strong nebular emission lines within certain redshift windows.

preprint2016arXiv

The Grism Lens-Amplified Survey from Space (GLASS). VI. Comparing the Mass and Light in MACSJ0416.1-2403 using Frontier Field imaging and GLASS spectroscopy

We present a strong and weak gravitational lens model of the galaxy cluster MACSJ0416.1-2403, constrained using spectroscopy from the Grism Lens-Amplified Survey from Space (GLASS) and Hubble Frontier Fields (HFF) imaging data. We search for emission lines in known multiply imaged sources in the GLASS spectra, obtaining secure spectroscopic redshifts of 31 multiple images belonging to 16 distinct source galaxies. The GLASS spectra provide the first spectroscopic measurements for 6 of the source galaxies. The weak lensing signal is acquired from 884 galaxies in the F606W HFF image. By combining the weak lensing constraints with 15 multiple image systems with spectroscopic redshifts and 9 multiple image systems with photometric redshifts, we reconstruct the gravitational potential of the cluster on an adaptive grid. The resulting total mass density map is compared with a stellar mass density map obtained from the deep Spitzer Frontier Fields imaging data to study the relative distribution of stellar and total mass in the cluster. We find that the projected stellar mass to total mass ratio, $f_{\star}$, varies considerably with the stellar surface mass density. The mean projected stellar mass to total mass ratio is $\langle f_{\star} \rangle= 0.009 \pm 0.003 $ (stat.), but with a systematic error as large as $0.004-0.005$, dominated by the choice of the IMF. We find agreement with several recent measurements of $f_{\star}$ in massive cluster environments. The lensing maps of convergence, shear, and magnification are made available to the broader community in the standard HFF format.

preprint2016arXiv

The Grism lens-amplified survey from space (GLASS). VII. The diversity of the distribution of star formation in cluster and field galaxies at 0.3<z<0.7

Exploiting the slitless spectroscopy taken as part of the Grism Lens-Amplified Survey from Space (GLASS), we present an extended analysis of the spatial distribution of star formation in 76 galaxies in 10 clusters at 0.3< z <0.7. We use 85 foreground and background galaxies in the same redshift range as a field sample. The samples are well matched in stellar mass (10^8-10^11 M_sun) and star formation rate (0.5-50 M_sun/yr). We visually classify galaxies in terms of broad-band morphology, Halpha morphology and likely physical process acting on the galaxy. Most Halpha emitters have a spiral morphology (41+/-8% in clusters, 51+/-8% in the field), followed by mergers/interactions (28+/-8%, 31+/-7%, respectively) and early-type galaxies (remarkably as high as 29+/-8% in clusters and 15+/-6% in the field). A diversity of Halpha morphologies is detected, suggesting a diversity of physical processes. In clusters, 30+/-8% of the galaxies present a regular morphology, mostly consistent with star formation diffused uniformly across the stellar population (mostly in the disk component, when present). The second most common morphology (28+/-8%) is asymmetric/jellyfish, consistent with ram pressure stripping or other non-gravitational processes in 18+/-8% of the cases. Ram pressure stripping appears significantly less prominent in the field (2+/-2%), where the most common morphology/mechanism appears to be consistent with minor gas rich mergers or clump accretion. This work demonstrates that while environment specific mechanisms affect galaxy evolution at this redshift, they are diverse and their effects subtle. A full understanding of this complexity requires larger samples and detailed and spatially resolved physical models.

preprint2016arXiv

The InfraRed Imaging Spectrograph (IRIS) for TMT: latest science cases and simulations

The Thirty Meter Telescope (TMT) first light instrument IRIS (Infrared Imaging Spectrograph) will complete its preliminary design phase in 2016. The IRIS instrument design includes a near-infrared (0.85 - 2.4 micron) integral field spectrograph (IFS) and imager that are able to conduct simultaneous diffraction-limited observations behind the advanced adaptive optics system NFIRAOS. The IRIS science cases have continued to be developed and new science studies have been investigated to aid in technical performance and design requirements. In this development phase, the IRIS science team has paid particular attention to the selection of filters, gratings, sensitivities of the entire system, and science cases that will benefit from the parallel mode of the IFS and imaging camera. We present new science cases for IRIS using the latest end-to-end data simulator on the following topics: Solar System bodies, the Galactic center, active galactic nuclei (AGN), and distant gravitationally-lensed galaxies. We then briefly discuss the necessity of an advanced data management system and data reduction pipeline.

preprint2015arXiv

A Local Baseline of the Black Hole Mass Scaling Relations for Active Galaxies. III. The BH mass - $σ$ relation

We create a baseline of the black hole (BH) mass (MBH) - stellar-velocity dispersion (sigma) relation for active galaxies, using a sample of 66 local (0.02<z<0.09) Seyfert-1 galaxies, selected from the Sloan Digital Sky Survey (SDSS). Analysis of SDSS images yields AGN luminosities free of host-galaxy contamination and morphological classification. 51/66 galaxies have spiral morphology. 28 bulges have Sersic index n<2 and are considered candidate pseudo bulges, with eight being definite pseudo bulges based on multiple classification criteria met. Only 4/66 galaxies show sign of interaction/merging. High signal-to-noise ratio Keck spectra provide the width of the broad Hbeta emission line free of FeII emission and stellar absorption. AGN luminosity and Hbeta line widths are used to estimate MBH. The Keck-based spatially-resolved kinematics is used to determine stellar-velocity dispersion within the spheroid effective radius. We find that sigma can vary on average by up to 40% across definitions commonly used in the literature, emphasizing the importance of using self-consistent definitions in comparisons and evolutionary studies. The MBH-sigma relation for our Seyfert-1 galaxies has the same intercept and scatter as that of reverberation-mapped AGNs as well as quiescent galaxies, consistent with the hypothesis that our single epoch MBH estimator and sample selection do not introduce significant biases. Barred galaxies, merging galaxies, and those hosting pseudo bulges do not represent outliers in the MBH-sigma relation. This is in contrast with previous work, although no firm conclusion can be drawn due to the small sample size and limited resolution of the SDSS images.

preprint2015arXiv

Constraints on the broad line region from regularized linear inversion: Velocity-delay maps for five nearby active galactic nuclei

Reverberation mapping probes the structure of the broad emission-line region (BLR) in active galactic nuclei (AGN). The kinematics of the BLR gas can be used to measure the mass of the central supermassive black hole. The main uncertainty affecting black hole mass determinations is the structure of the BLR. We present a new method for reverberation mapping based on regularized linear inversion (RLI) that includes modelling of the AGN continuum light curves. This enables fast calculation of velocity-resolved response maps to constrain BLR structure. RLI allows for negative response, such as when some areas of the BLR respond in inverse proportion to a change in ionizing continuum luminosity. We present time delays, integrated response functions, and velocity-delay maps for the $\rm{H}\,β$ broad emission line in five nearby AGN, as well as for $\rm{H}\,α$ and $\rm{H}\,γ$ in Arp 151, using data from the Lick AGN Monitoring Project 2008. We find indications of prompt response in three of the objects (Arp 151, NGC 5548 and SBS 1116+583A) with additional prompt response in the red wing of $\rm{H}\,β$. In SBS 1116+583A we find evidence for a multimodal broad prompt response followed by a second narrow response at 10 days. We find no clear indications of negative response. The results are complementary to, and consistent with, other methods such as cross correlation, maximum entropy and dynamical modelling. Regularized linear inversion with continuum light curve modelling provides a fast, complementary method for velocity-resolved reverberation mapping and is suitable for use on large datasets.

preprint2015arXiv

Correcting the z~8 Galaxy Luminosity Function for Gravitational Lensing Magnification Bias

We present a Bayesian framework to account for the magnification bias from both strong and weak gravitational lensing in estimates of high-redshift galaxy luminosity functions. We illustrate our method by estimating the $z\sim8$ UV luminosity function using a sample of 97 Y-band dropouts (Lyman break galaxies) found in the Brightest of Reionizing Galaxies (BoRG) survey and from the literature. We find the luminosity function is well described by a Schechter function with characteristic magnitude of $M^\star = -19.85^{+0.30}_{-0.35}$, faint-end slope of $α= -1.72^{+0.30}_{-0.29}$, and number density of $\log_{10} Ψ^\star [\textrm{Mpc}^{-3}] = -3.00^{+0.23}_{-0.31}$. These parameters are consistent within the uncertainties with those inferred from the same sample without accounting for the magnification bias, demonstrating that the effect is small for current surveys at $z\sim8$, and cannot account for the apparent overdensity of bright galaxies compared to a Schechter function found recently by Bowler et al. (2014a,b) and Finkelstein et al. (2014). We estimate that the probability of finding a strongly lensed $z\sim8$ source in our sample is in the range $\sim 3-15 \%$ depending on limiting magnitude. We identify one strongly-lensed candidate and three cases of intermediate lensing in BoRG (estimated magnification $μ>1.4$) in addition to the previously known candidate group-scale strong lens. Using a range of theoretical luminosity functions we conclude that magnification bias will dominate wide field surveys -- such as those planned for the Euclid and WFIRST missions -- especially at $z>10$. Magnification bias will need to be accounted for in order to derive accurate estimates of high-redshift luminosity functions in these surveys and to distinguish between galaxy formation models.

preprint2015arXiv

Discovery of two gravitationally lensed quasars in the Dark Energy Survey

We present spectroscopic confirmation of two new lensed quasars via data obtained at the 6.5m Magellan/Baade Telescope. The lens candidates have been selected from the Dark Energy Survey (DES) and WISE based on their multi-band photometry and extended morphology in DES images. Images of DES J0115-5244 show two blue point sources at either side of a red galaxy. Our long-slit data confirm that both point sources are images of the same quasar at $z_{s}=1.64.$ The Einstein Radius estimated from the DES images is $0.51$". DES J2200+0110 is in the area of overlap between DES and the Sloan Digital Sky Survey (SDSS). Two blue components are visible in the DES and SDSS images. The SDSS fiber spectrum shows a quasar component at $z_{s}=2.38$ and absorption compatible with Mg II and Fe II at $z_{l}=0.799$, which we tentatively associate with the foreground lens galaxy. The long-slit Magellan spectra show that the blue components are resolved images of the same quasar. The Einstein Radius is $0.68$" corresponding to an enclosed mass of $1.6\times10^{11}\,M_{\odot}.$ Three other candidates were observed and rejected, two being low-redshift pairs of starburst galaxies, and one being a quasar behind a blue star. These first confirmation results provide an important empirical validation of the data-mining and model-based selection that is being applied to the entire DES dataset.

preprint2015arXiv

High resolution imaging and spectroscopy of the gravitational lens SDSSJ1206+4332: a natural coronagraph at $z=1.789$ and a standard ruler at $z=0.745$

We present spectroscopy and laser guide star adaptive optics (LGSAO) images of the doubly imaged lensed quasar SDSS J1206+4332. We revise the deflector redshift proposed previously to $z_{d}=0.745,$ and measure for the first time its velocity dispersion $σ=(290\pm30)$ km/s. The LGSAO data show the lensed quasar host galaxy stretching over the astroid caustic thus forming an extra pair of merging images, which was previously thought to be an unrelated galaxy in seeing limited data. Owing to the peculiar geometry, the lens acts as a natural coronagraph on the broad-line region of the quasar so that only narrow [O III] emission is found in the fold arc. We use the data to reconstruct the source structure and deflector potential, including nearby perturbers. We reconstruct the point-spread function (PSF) from the quasar images themselves, since no additional point source is present in the field of view. From gravitational lensing and stellar dynamics, we find the slope of the total mass density profile to be $γ^{\prime}=-\logρ/\log r =1.93\pm0.09.$ We discuss the potential of SDSS J1206+4332 for measuring time delay distance (and thus H$_0$ and other cosmological parameters), or as a standard ruler, in combination with the time delay published by the COSMOGRAIL collaboration. We conclude that this system is very promising for cosmography. However, in order to achieve competitive precision and accuracy, an independent characterization of the PSF is needed. Spatially resolved kinematics of the deflector would reduce the uncertainties further. Both are within the reach of current observational facilities.

preprint2015arXiv

Illuminating a Dark Lens : A Type Ia Supernova Magnified by the Frontier Fields Galaxy Cluster Abell 2744

SN HFF14Tom is a Type Ia Supernova (SN) discovered at z = 1.3457 +- 0.0001 behind the galaxy cluster Abell 2744 (z = 0.308). In a cosmology-independent analysis, we find that HFF14Tom is 0.77 +- 0.15 magnitudes brighter than unlensed Type Ia SNe at similar redshift, implying a lensing magnification of mu_obs = 2.03 +- 0.29. This observed magnification provides a rare opportunity for a direct empirical test of galaxy cluster lens models. Here we test 17 lens models, 13 of which were generated before the SN magnification was known, qualifying as pure "blind tests". The models are collectively fairly accurate: 8 of the models deliver median magnifications that are consistent with the measured mu to within 1-sigma. However, there is a subtle systematic bias: the significant disagreements all involve models overpredicting the magnification. We evaluate possible causes for this mild bias, and find no single physical or methodological explanation to account for it. We do find that model accuracy can be improved to some extent with stringent quality cuts on multiply-imaged systems, such as requiring that a large fraction have spectroscopic redshifts. In addition to testing model accuracies as we have done here, Type Ia SN magnifications could also be used as inputs for future lens models of Abell 2744 and other clusters, providing valuable constraints in regions where traditional strong- and weak-lensing information is unavailable.

preprint2015arXiv

Luminous and Dark Matter Profiles from Galaxies to Clusters: Bridging the Gap with Group-scale Lenses

Observations of strong gravitational lensing, stellar kinematics, and larger-scale tracers enable accurate measures of the distribution of dark matter (DM) and baryons in massive early-type galaxies (ETGs). While such techniques have been applied to galaxy-scale and cluster-scale lenses, the paucity of intermediate-mass systems with high-quality data has precluded a uniform analysis of mass-dependent trends. With the aim of bridging this gap, we present new observations and analyses of 10 group-scale lenses at <z>=0.36 characterized by Einstein radii theta_Ein=2.5"-5.1" and a mean halo mass of M_200=10^14.0 Msol. We measure a mean concentration c_200=5.0+-0.8 consistent with unmodified cold dark matter halos. By combining our data with other lens samples, we analyze the mass structure of ETGs in 10^13 Msol-10^15 Msol halos using homogeneous techniques. We show that the slope of the total density profile gamma_tot within the effective radius depends on the stellar surface density, as demonstrated previously, but also on the halo mass. We analyze these trends using halo occupation models and resolved stellar kinematics with the goal of testing the universality of the DM profile. Whereas the central galaxies of clusters require a shallow inner DM density profile, group-scale lenses are consistent with a Navarro-Frenk-White profile or one that is slightly contracted. The largest uncertainties arise from the sample size and likely radial gradients in stellar populations. We conclude that the net effect of baryons on the DM distribution may not be universal, but more likely varies with halo mass due to underlying trends in star formation efficiency and assembly history.

preprint2015arXiv

Modeling reverberation mapping data II: dynamical modeling of the Lick AGN Monitoring Project 2008 dataset

We present dynamical modeling of the broad line region (BLR) for a sample of five Seyfert 1 galaxies using reverberation mapping data taken by the Lick AGN Monitoring Project in 2008. By modeling the AGN continuum light curve and H$β$ line profiles directly we are able to constrain the geometry and kinematics of the BLR and make a measurement of the black hole mass that does not depend upon the virial factor, $f$, needed in traditional reverberation mapping analysis. We find that the geometry of the BLR is generally a thick disk viewed close to face-on. While the H$β$ emission is found to come preferentially from the far side of the BLR, the mean size of the BLR is consistent with the lags measured with cross-correlation analysis. The BLR kinematics are found to be consistent with either inflowing motions or elliptical orbits, often with some combination of the two. We measure black hole masses of $\log_{10}(M_{\rm\,BH}/M_\odot)=6.62^{+0.10}_{-0.13}$ for Arp 151, $7.42^{+0.26}_{-0.27}$ for Mrk 1310, $7.51^{+0.23}_{-0.14}$ for NGC 5548, $6.42^{+0.24}_{-0.18}$ for NGC 6814, and $6.99^{+0.32}_{-0.25}$ for SBS 1116+583A. The $f$ factors measured individually for each AGN are found to correlate with inclination angle, although not with $M_{\rm\,BH}$, $L_{5100}$, or FWHM/$σ$ of the emission line profile.

preprint2015arXiv

Multiple Images of a Highly Magnified Supernova Formed by an Early-Type Cluster Galaxy Lens

In 1964, Refsdal hypothesized that a supernova whose light traversed multiple paths around a strong gravitational lens could be used to measure the rate of cosmic expansion. We report the discovery of such a system. In Hubble Space Telescope imaging, we have found four images of a single supernova forming an Einstein cross configuration around a redshift z=0.54 elliptical galaxy in the MACS J1149.6+2223 cluster. The cluster's gravitational potential also creates multiple images of the z=1.49 spiral supernova host galaxy, and a future appearance of the supernova elsewhere in the cluster field is expected. The magnifications and staggered arrivals of the supernova images probe the cosmic expansion rate, as well as the distribution of matter in the galaxy and cluster lenses.

preprint2015arXiv

Precision cosmology with time delay lenses: high resolution imaging requirements

Lens time delays are a powerful probe of cosmology, provided that the gravitational potential of the main deflector can be modeled with sufficient precision. Recent work has shown that this can be achieved by detailed modeling of the host galaxies of lensed quasars, which appear as "Einstein Rings" in high resolution images. We carry out a systematic exploration of the high resolution imaging required to exploit the thousands of lensed quasars that will be discovered by current and upcoming surveys with the next decade. Specifically, we simulate realistic lens systems as imaged by the Hubble Space Telescope (HST), James Webb Space Telescope (JWST), and ground based adaptive optics images taken with Keck or the Thirty Meter Telescope (TMT). We compare the performance of these pointed observations with that of images taken by the Euclid (VIS), Wide-Field Infrared Survey Telescope (WFIRST) and Large Synoptic Survey Telescope (LSST) surveys. We use as our metric the precision with which the slope $γ'$ of the total mass density profile $ρ_{tot}\propto r^{-γ'}$ for the main deflector can be measured. Ideally, we require that the statistical error on $γ'$ be less than 0.02, such that it is subdominant to other sources of random and systematic uncertainties. We find that survey data will likely have sufficient depth and resolution to meet the target only for the brighter gravitational lens systems, comparable to those discovered by the SDSS survey. For fainter systems, that will be discovered by current and future surveys, targeted follow-up will be required. However, the exposure time required with upcoming facilitites such as JWST, the Keck Next Generation Adaptive Optics System, and TMT, will only be of order a few minutes per system, thus making the follow-up of hundreds of systems a practical and efficient cosmological probe.

preprint2015arXiv

RCS2 J232727.6-020437: An Efficient Cosmic Telescope at $z=0.6986$

We present a detailed gravitational lens model of the galaxy cluster RCS2 J232727.6-020437. Due to cosmological dimming of cluster members and ICL, its high redshift ($z=0.6986$) makes it ideal for studying background galaxies. Using new ACS and WFC3/IR HST data, we identify 16 multiple images. From MOSFIRE follow up, we identify a strong emission line in the spectrum of one multiple image, likely confirming the redshift of that system to $z=2.083$. With a highly magnified ($μ\gtrsim2$) source plane area of $\sim0.7$ arcmin$^2$ at $z=7$, RCS2 J232727.6-020437 has a lensing efficiency comparable to the Hubble Frontier Fields clusters. We discover four highly magnified $z\sim7$ candidate Lyman-break galaxies behind the cluster, one of which may be multiply-imaged. Correcting for magnification, we find that all four galaxies are fainter than $0.5 L_{\star}$. One candidate is detected at ${>10σ}$ in both Spitzer/IRAC [3.6] and [4.5] channels. A spectroscopic follow-up with MOSFIRE does not result in the detection of the Lyman-alpha emission line from any of the four candidates. From the MOSFIRE spectra we place median upper limits on the Lyman-alpha flux of $5-14 \times 10^{-19}\, \mathrm{erg \,\, s^{-1} cm^{-2}}$ ($5σ$).

preprint2015arXiv

Science with KRAKENS

The Keck science community is entering an era of unprecedented change. Powerful new instrument like ZTF, JWST, LSST, and the ELTs will catalyze this change, and we must be ready to take full advantage to maintain our position of scientific leadership. The best way to do this is to continue the UC and Caltech tradition of technical excellence in instrumentation. In this whitepaper we describe a new instrument called KRAKENS to help meet these challenges. KRAKENS uses a unique detector technology (MKIDs) to enable groundbreaking science across a wide range of astrophysical research topics. This document will lay out the detailed expected science return of KRAKENS.

preprint2015arXiv

The Galaxy UV Luminosity Function Before the Epoch of Reionization

We present a model for the evolution of the galaxy ultraviolet (UV) luminosity function (LF) across cosmic time where star formation is linked to the assembly of dark matter halos under the assumption of a mass dependent, but redshift independent, efficiency. We introduce a new self-consistent treatment of the halo star formation history, which allows us to make predictions at $z>10$ (lookback time $\lesssim500$ Myr), when growth is rapid. With a calibration at a single redshift to set the stellar-to-halo mass ratio, and no further degrees of freedom, our model captures the evolution of the UV LF over all available observations ($0\lesssim z\lesssim10$). The significant drop in luminosity density of currently detectable galaxies beyond $z\sim8$ is explained by a shift of star formation toward less massive, fainter galaxies. Assuming that star formation proceeds down to atomic cooling halos, we derive a reionization optical depth $τ= 0.056^{+0.007}_{-0.010}$, fully consistent with the latest Planck measurement, implying that the universe is fully reionized at $z=7.84^{+0.65}_{-0.98}$. In addition, our model naturally produces smoothly rising star formation histories for galaxies with $L\lesssim L_*$ in agreement with observations and hydrodynamical simulations. Before the epoch of reionization at $z>10$ we predict the LF to remain well-described by a Schechter function, but with an increasingly steep faint-end slope ($α\sim-3.5$ at $z\sim16$). Finally, we construct forecasts for surveys with \JWST~and \WFIRST and predict that galaxies out to $z\sim14$ will be observed. Galaxies at $z>15$ will likely be accessible to JWST and WFIRST only through the assistance of strong lensing magnification.

preprint2015arXiv

The Grism Lens-Amplified Survey from Space (GLASS). II. Gas-phase metallicity and radial gradients in an interacting system at z~2

We present spatially resolved gas-phase metallicity for a system of three galaxies at z=1.85 detected in the Grism Lens-Amplified Survey from Space (GLASS). The combination of HST's diffraction limit and strong gravitational lensing by the cluster MACS J0717+3745 results in a spatial resolution of ~200-300 pc, enabling good spatial sampling despite the intrinsically small galaxy sizes. The galaxies in this system are separated by 50-200 kpc in projection and are likely in an early stage of interaction, evidenced by relatively high specific star formation rates. Their gas-phase metallicities are consistent with larger samples at similar redshift, star formation rate, and stellar mass. We obtain a precise measurement of the metallicity gradient for one galaxy and find a shallow slope compared to isolated galaxies at high redshift, consistent with a flattening of the gradient due to gravitational interaction. An alternative explanation for the shallow metallicity gradient and elevated star formation rate is rapid recycling of metal-enriched gas, but we find no evidence for enhanced gas-phase metallicities which should result from this effect. Notably, the measured stellar masses log(M/Msun) = 7.2-9.1 probe to an order of magnitude below previous mass-metallicity studies at this redshift. The lowest mass galaxy has properties similar to those expected for Fornax at this redshift, indicating that GLASS is able to directly study the progenitors of local group dwarf galaxies on spatially resolved scales. Larger samples from the full GLASS survey will be ideal for studying the effects of feedback, and the time evolution of metallicity gradients. These initial results demonstrate the utility of HST spectroscopy combined with gravitational lensing for characterizing resolved physical properties of galaxies at high redshift.

preprint2015arXiv

The Grism Lens-Amplified Survey from Space (GLASS). V. Extent and spatial distribution of star formation in z~0.5 cluster galaxies

We present the first study of the spatial distribution of star formation in z~0.5 cluster galaxies. The analysis is based on data taken with the Wide Field Camera 3 as part of the Grism Lens-Amplified Survey from Space (GLASS). We illustrate the methodology by focusing on two clusters (MACS0717.5+3745 and MACS1423.8+2404) with different morphologies (one relaxed and one merging) and use foreground and background galaxies as field control sample. The cluster+field sample consists of 42 galaxies with stellar masses in the range 10^8-10^11 M_sun, and star formation rates in the range 1-20 M_sun/yr. Both in clusters and in the field, Hα is more extended than the rest-frame UV continuum in 60% of the cases, consistent with diffuse star formation and inside out growth. In ~20% of the cases, the Hα emission appears more extended in cluster galaxies than in the field, pointing perhaps to ionized gas being stripped and/or star formation being enhanced at large radii. The peak of the Hα emission and that of the continuum are offset by less than 1 kpc. We investigate trends with the hot gas density as traced by the X-ray emission, and with the surface mass density as inferred from gravitational lens models and find no conclusive results. The diversity of morphologies and sizes observed in H_alpha illustrates the complexity of the environmental process that regulate star formation. Upcoming analysis of the full GLASS dataset will increase our sample size by almost an order of magnitude, verifying and strengthening the inference from this initial dataset.

preprint2015arXiv

The Sloan Lens ACS Survey. XII. Extending Strong Lensing to Lower Masses

We present observational results from a new Hubble Space Telescope (HST) Snapshot program to extend the methods of the Sloan Lens ACS (SLACS) Survey to lower lens-galaxy masses. We discover 40 new galaxy-scale strong lenses, which we supplement with 58 previously discovered SLACS lenses. In addition, we determine the posterior PDFs of the Einstein radius for 33 galaxies (18 new and 15 from legacy SLACS data) based on single lensed images. We find a less-than-unity slope of $0.64\pm0.06$ for the $\log_{10} σ_*$-$\log_{10} σ_{\rm SIE}$ relation, which corresponds to a 6-$σ$ evidence that the total mass-density profile of early-type galaxies varies systematically in the sense of being shallower at higher lens-galaxy velocity dispersions. The trend is only significant when single-image systems are considered, highlighting the importance of including both "lenses" and "non-lenses" for an unbiased treatment of the lens population when extending to lower mass ranges. By scaling simple stellar population models to the HST I-band data, we identify a strong trend of increasing dark-matter fraction at higher velocity dispersions, which can be alternatively interpreted as a trend in the stellar initial mass function (IMF) normalization. Consistent with previous findings and the suggestion of a non-universal IMF, we find that a Salpeter IMF is ruled out for galaxies with velocity dispersion less than $180$ km/s. Considered together, our mass-profile and dark-matter-fraction trends with increasing galaxy mass could both be explained by an increasing relative contribution on kiloparsec scales from a dark-matter halo with a spatial profile more extended than that of the stellar component.

preprint2014arXiv

AMUSE-Field II. Nucleation of early-type galaxies in the field vs. cluster environment

The optical light profiles of nearby early type galaxies are known to exhibit a smooth transition from nuclear light deficits to nuclear light excesses with decreasing galaxy mass, with as much as 80 per cent of the galaxies with stellar masses below 10^10 Msun hosting a massive nuclear star cluster. At the same time, while all massive galaxies are thought to harbor nuclear super-massive black holes (SMBHs), observational evidence for SMBHs is slim at the low end of the mass function. Here, we explore the environmental dependence of the nucleation fraction by comparing two homogeneous samples of nearby field vs. cluster early type galaxies with uniform Hubble Space Telescope (HST) coverage. Existing Chandra X-ray Telescope data for both samples yield complementary information on low-level accretion onto nuclear SMBHs. Specifically, we report on dual-band (F475W & F850LP) Advanced Camera for Surveys (ACS) imaging data for 28 out of the 103 field early type galaxies that compose the AMUSE-Field Chandra survey, and compare our results against the companion HST and Chandra surveys for a sample of 100 Virgo cluster early types (ACS Virgo Cluster and AMUSE-Virgo surveys, respectively). We model the two-dimensional light profiles of the field targets to identify and characterize NSCs, and find a field nucleation fraction of 26% (+17%, -11%; at the 1-sigma level), consistent with the measured Virgo nucleation fraction of 30% (+17%, -12%), across a comparable mass distribution. Coupled with the Chandra result that SMBH activity is higher for the field, our findings indicate that, since the last epoch of star formation, the funneling of gas to the nuclear regions has been inhibited more effectively for Virgo galaxies, arguably via ram pressure stripping.

preprint2014arXiv

Cosmic Evolution of Black Holes and Spheroids. V. The Relation Between Black Hole Mass and Host Galaxy Luminosity for a Sample of 79 Active Galaxies

We investigate the cosmic evolution of the black hole (BH) mass -- bulge luminosity relation using a sample of 52 active galaxies at $z \sim 0.36$ and $z \sim 0.57$ in the BH mass range of $10^{7.4-9.1} M_{\odot}$. By consistently applying multi-component spectral and structural decomposition to high-quality Keck spectra and high-resolution HST images, BH masses ($M_{\rm BH}$) are estimated using the H$β$ broad emission line combined with the 5100 Å nuclear luminosity, and bulge luminosities ($L_{\rm bul}$) are derived from surface photometry. Comparing the resulting $M_{\rm BH}-L_{\rm bul}$ relation to local active galaxies and taking into account selection effects, we find evolution of the form $M_{\rm BH} / L_{\rm bul} \propto (1+z)^γ$ with $γ=1.8\pm0.7$, consistent with BH growth preceding that of the host galaxies. Including an additional sample of 27 active galaxies with $0.5<z<1.9$ taken from the literature and measured in a consistent way, we obtain $γ=0.9\pm0.7$ for the $M_{\rm BH}-L_{\rm bul}$ relation and $γ=0.4\pm0.5$ for the $M_{\rm BH}$--total host galaxy luminosity ($L_{\rm host}$) relation. The results strengthen the findings from our previous studies and provide additional evidence for host-galaxy bulge growth being dominated by disk-to-bulge transformation via minor mergers and/or disk instabilities.

preprint2014arXiv

Data Mining for Gravitationally Lensed Quasars

Gravitationally lensed (GL) quasars are brighter than their unlensed counterparts and produce images with distinctive morphological signatures. Past searches and target selection algorithms, in particular the Sloan Quasar Lens Search (SQLS), have relied on basic morphological criteria, which were applied to samples of bright, spectroscopically confirmed quasars. The SQLS techniques are not sufficient for searching into new surveys (e.g. DES, PS1, LSST), because spectroscopic information is not readily available and the large data volume requires higher purity in target/candidate selection. We carry out a systematic exploration of machine learning techniques and demonstrate that a two step strategy can be highly effective. In the first step we use catalog-level information ($griz$+WISE magnitudes, second moments) to preselect targets, using artificial neural networks. The accepted targets are then inspected with pixel-by-pixel pattern recognition algorithms (Gradient-Boosted Trees), to form a final set of candidates. The results from this procedure can be used to further refine the simpler SQLS algorithms, with a twofold (or threefold) gain in purity and the same (or $80\%$) completeness at target-selection stage, or a purity of $70\%$ and a completeness of $60\%$ after the candidate-selection step. Simpler photometric searches in $griz$+WISE based on colour cuts would provide samples with $7\%$ purity or less. Our technique is extremely fast, as a list of candidates can be obtained from a stage III experiment (e.g. DES catalog/database) in {a few} CPU hours. The techniques are easily extendable to Stage IV experiments like LSST with the addition of time domain information.

preprint2014arXiv

Gravitational Lensing - Einstein's Unfinished Symphony

Gravitational lensing - the deflection of light rays by gravitating matter - has become a major tool in the armoury of the modern cosmologist. Proposed nearly a hundred years ago as a key feature of Einstein's theory of General Relativity, we trace the historical development since its verification at a solar eclipse in 1919. Einstein was apparently cautious about its practical utility and the subject lay dormant observationally for nearly 60 years. Nonetheless there has been rapid progress over the past twenty years. The technique allows astronomers to chart the distribution of dark matter on large and small scales thereby testing predictions of the standard cosmological model which assumes dark matter comprises a massive weakly-interacting particle. By measuring distances and tracing the growth of dark matter structure over cosmic time, gravitational lensing also holds great promise in determining whether the dark energy, postulated to explain the accelerated cosmic expansion, is a vacuum energy density or a failure of General Relativity on large scales. We illustrate the wide range of applications which harness the power of gravitational lensing, from searches for the earliest galaxies magnified by massive clusters to those for extrasolar planets which temporarily brighten a background star. We summarise the future prospects with dedicated ground and space-based facilities designed to exploit this remarkable physical phenomenon.

preprint2014arXiv

Modeling reverberation mapping data I: improved geometric and dynamical models and comparison with cross-correlation results

We present an improved and expanded simply parameterized phenomenological model of the broad line region (BLR) in active galactic nuclei (AGN) for modeling reverberation mapping data. By modeling reverberation mapping data directly, we can constrain the geometry and dynamics of the BLR and measure the black hole mass without relying on the normalization factor needed in the traditional analysis. For realistic simulated reverberation mapping datasets of high-quality, we can recover the black hole mass to $0.05-0.25$ dex uncertainty and distinguish between dynamics dominated by elliptical orbits and inflowing gas. While direct modeling of the integrated emission line light curve allows for measurement of the mean time lag, other details of the geometry of the BLR are better constrained by the full spectroscopic dataset of emission line profiles. We use this improved model of the BLR to explore possible sources of uncertainty in measurements of the time lag using cross-correlation function (CCF) analysis and in measurements of the black hole mass using the virial product. Sampling the range of geometries and dynamics in our model of the BLR suggests that the theoretical uncertainty in black hole masses measured using the virial product is on the order of 0.25 dex. These results support the use of the CCF to measure time lags and the virial product to measure black hole masses when direct modeling techniques cannot be applied, provided the uncertainties associated with the interpretation of the results are taken into account.

preprint2014arXiv

Purely dry mergers do not explain the observed evolution of massive early-type galaxies since z = 1

Several studies have suggested that the observed size evolution of massive early-type galaxies (ETGs) can be explained as a combination of dry mergers and progenitor bias, at least since z\sim1. In this paper we carry out a new test of the dry-merger scenario based on recent lensing measurements of the evolution of the mass density profile of ETGs. We construct a theoretical model for the joint evolution of the size and mass density profile slope γ' driven by dry mergers occurring at rates given by cosmological simulations. Such dry-merger model predicts a strong decrease of γ' with cosmic time, inconsistent with the almost constant γ' inferred from observations in the redshift range 0<z<1. We then show with a simple toy model that a modest amount of cold gas in the mergers -- consistent with the upper limits on recent star formation in ETGs -- is sufficient to reconcile the model with measurements of γ'. By fitting for the amount of gas accreted during mergers, we find that models with dissipation are consistent with observations of the evolution in both size and density slope, if \sim4% of the total final stellar mass arises from the gas accreted since z\sim1. Purely dry merger models are ruled out at > 99% CL. We thus suggest a scenario where the outer regions of massive ETGs grow by accretion of stars and dark matter, while small amounts of dissipation and nuclear star formation conspire to keep the mass density profile constant and approximately isothermal.

preprint2014arXiv

RingFinder: automated detection of galaxy-scale gravitational lenses in ground-based multi-filter imaging data

We present RingFinder, a tool for finding galaxy-scale strong gravitational lenses in multiband imaging data. By construction, the method is sensitive to configurations involving a massive foreground early-type galaxy and a faint, background, blue source. RingFinder detects the presence of blue residuals embedded in an otherwise smooth red light distribution by difference imaging in two bands. The method is automated for efficient application to current and future surveys, having originally been designed for the 150-deg2 Canada France Hawaii Telescope Legacy Survey (CFHTLS). We describe each of the steps of RingFinder. We then carry out extensive simulations to assess completeness and purity. For sources with magnification mu>4, RingFinder reaches 42% (resp. 25%) completeness and 29% (resp. 86%) purity before (resp. after) visual inspection. The completeness of RingFinder is substantially improved in the particular range of Einstein radii 0.8 < REin < 2. and lensed images brighter than g = 22.5, where it can be as high as 70%. RingFinder does not introduce any significant bias in the source or deflector population. We conclude by presenting the final catalog of RingFinder CFHTLS galaxy-scale strong lens candidates. Additional information obtained with Hubble Space Telescope and Keck Adaptive Optics high resolution imaging, and with Keck and Very Large Telescope spectroscopy, is used to assess the validity of our classification, and measure the redshift of the foreground and the background objects. From an initial sample of 640,000 early type galaxies, RingFinder returns 2500 candidates, which we further reduce by visual inspection to 330 candidates. We confirm 33 new gravitational lenses from the main sample of candidates, plus an additional 16 systems taken from earlier versions of RingFinder. First applications are presented in the SL2S galaxy-scale Lens Sample paper series.

preprint2014arXiv

Scaling relations between black holes and their host galaxies: comparing theoretical and observational measurements, and the impact of selection effects

We use the high-resolution simulation MassiveBlackII to examine scaling relations between black hole mass (MBH) and host galaxy properties (sigma, M*, and LV), finding good agreement with observational data, especially at the high-mass end. The simulations have less intrinsic scatter than observations, and the MBH-LV correlation has the largest scatter, suggesting it may the the least fundamental of the three relations. We find Gaussian scatter about all three relations, except among the highest mass galaxies, which host more massive black holes. Below z~2 the slopes for the full population remain roughly z-independent, and only steepen by 50% by z~4. The normalization of the sigma, LV relations evolve by 0.3, 0.43 dex, while the MBH correlation does not evolve to at least z~2. Testing for selection biases, we find samples selected by MBH or M* have steeper slopes than randomly selected samples. If unaccounted for, such a selection function would find faster evolution than inferred from a randomly selected sample, as objects at the highend of the relation tend to evolve more rapidly. We find a potential bias among high-LBH subsamples (tending to reside in higher mass galaxies), but these bright-AGN exhibit no intrinsic bias relative to fainter ones in equivalent-mass hosts, nor is there a significant difference between active- and inactive-samples. Finally we characterize the evolution of individual black holes along the scaling planes. Below the local relation, black holes grow faster than their host (72% of black holes 0.3 dex below the mean relation have a MBH-M* trajectory steeper than the local relation), while those above have shallower trajectories (only 14% are steeper than local). Black holes tend to grow faster than their hosts until surpassing the local relation, at which point their growth is suppressed while their hosts continue to grow, returning them to the mean relation.

preprint2014arXiv

Spectroscopic Confirmation of the Rich z=1.80 Galaxy Cluster JKCS 041 Using the WFC3 Grism: Environmental Trends in the Ages and Structure of Quiescent Galaxies

We present Hubble Space Telescope imaging and grism spectroscopy in the field of the distant galaxy cluster JKCS 041 using the Wide Field Camera 3. We confirm that JKCS 041 is a rich cluster and derive a redshift z=1.80 via the spectroscopic identification of 19 member galaxies, of which 15 are quiescent. These are centered upon diffuse X-ray emission seen by the Chandra observatory. As JKCS 041 is the most distant known cluster with such a large, spectroscopically confirmed quiescent population, it provides a unique opportunity to study the effect of the environment on galaxy properties at early epochs. We construct high-quality composite spectra of the quiescent cluster members that reveal prominent Balmer and metallic absorption lines. Using these, we measure the mean stellar ages in two bins of stellar mass. The quiescent cluster members' ages agree remarkably closely with that inferred by Whitaker et al. for similarly selected samples in the field, supporting the idea that the cluster environment is more efficient at truncating star formation while not having a strong effect on the mean epoch of quenching. We find some evidence (90% confidence) for a lower fraction of disk-like quiescent systems in JKCS 041 compared to a sample of coeval field galaxies drawn from the CANDELS survey. Taking this into account, we do not detect a significant difference between the mass-radius relations of the quiescent JKCS 041 members and our z~1.8 field sample. Finally, we demonstrate how differences in the morphological mixture of quenched systems can complicate measures of the environmental dependence of size growth.

preprint2014arXiv

Spitzer UltRa Faint SUrvey Program (SURFS UP) I: An Overview

SURFSUP is a joint Spitzer and HST Exploration Science program using 10 galaxy clusters as cosmic telescopes to study z >~ 7 galaxies at intrinsically lower luminosities, enabled by gravitational lensing, than blank field surveys of the same exposure time. Our main goal is to measure stellar masses and ages of these galaxies, which are the most likely sources of the ionizing photons that drive reionization. Accurate knowledge of the star formation density and star formation history at this epoch is necessary to determine whether these galaxies indeed reionized the universe. Determination of the stellar masses and ages requires measuring rest frame optical light, which only Spitzer can probe for sources at z >~ 7, for a large enough sample of typical galaxies. Our program consists of 550 hours of Spitzer/IRAC imaging covering 10 galaxy clusters with very well-known mass distributions, making them extremely precise cosmic telescopes. We combine our data with archival observations to obtain mosaics with ~30 hours exposure time in both 3.6$μ$m and 4.5$μ$m in the central 4 arcmin x 4 arcmin field and ~15 hours in the flanking fields. This results in 3-$σ$ sensitivity limits of ~26.6 and ~26.2AB magnitudes for the central field in the IRAC 3.6 and 4.5$μ$m bands, respectively. To illustrate the survey strategy and characteristics we introduce the sample, present the details of the data reduction and demonstrate that these data are sufficient for in-depth studies of z >~ 7 sources (using a z=9.5 galaxy behind MACSJ1149.5+2223 as an example). For the first cluster of the survey (the Bullet Cluster) we have released all high-level data mosaics and IRAC empirical PSF models. In the future we plan to release these data products for the entire survey.

preprint2014arXiv

Strong Lens Time Delay Challenge: I. Experimental Design

The time delays between point-like images in gravitational lens systems can be used to measure cosmological parameters. The number of lenses with measured time delays is growing rapidly; the upcoming \emph{Large Synoptic Survey Telescope} (LSST) will monitor $\sim10^3$ strongly lensed quasars. In an effort to assess the present capabilities of the community to accurately measure the time delays, and to provide input to dedicated monitoring campaigns and future LSST cosmology feasibility studies, we have invited the community to take part in a "Time Delay Challenge" (TDC). The challenge is organized as a set of "ladders," each containing a group of simulated datasets to be analyzed blindly by participating teams. Each rung on a ladder consists of a set of realistic mock observed lensed quasar light curves, with the rungs' datasets increasing in complexity and realism. The initial challenge described here has two ladders, TDC0 and TDC1. TDC0 has a small number of datasets, and is designed to be used as a practice set by the participating teams. The (non-mandatory) deadline for completion of TDC0 was the TDC1 launch date, December 1, 2013. The TDC1 deadline was July 1 2014. Here we give an overview of the challenge, we introduce a set of metrics that will be used to quantify the goodness-of-fit, efficiency, precision, and accuracy of the algorithms, and we present the results of TDC0. Thirteen teams participated in TDC0 using 47 different methods. Seven of those teams qualified for TDC1, which is described in the companion paper II.

preprint2014arXiv

Strong Lens Time Delay Challenge: II. Results of TDC1

We present the results of the first strong lens time delay challenge. The motivation, experimental design, and entry level challenge are described in a companion paper. This paper presents the main challenge, TDC1, which consisted of analyzing thousands of simulated light curves blindly. The observational properties of the light curves cover the range in quality obtained for current targeted efforts (e.g.,~COSMOGRAIL) and expected from future synoptic surveys (e.g.,~LSST), and include simulated systematic errors. \nteamsA\ teams participated in TDC1, submitting results from \nmethods\ different method variants. After a describing each method, we compute and analyze basic statistics measuring accuracy (or bias) $A$, goodness of fit $χ^2$, precision $P$, and success rate $f$. For some methods we identify outliers as an important issue. Other methods show that outliers can be controlled via visual inspection or conservative quality control. Several methods are competitive, i.e., give $|A|<0.03$, $P<0.03$, and $χ^2<1.5$, with some of the methods already reaching sub-percent accuracy. The fraction of light curves yielding a time delay measurement is typically in the range $f = $20--40\%. It depends strongly on the quality of the data: COSMOGRAIL-quality cadence and light curve lengths yield significantly higher $f$ than does sparser sampling. Taking the results of TDC1 at face value, we estimate that LSST should provide around 400 robust time-delay measurements, each with $P<0.03$ and $|A|<0.01$, comparable to current lens modeling uncertainties. In terms of observing strategies, we find that $A$ and $f$ depend mostly on season length, while P depends mostly on cadence and campaign duration.

preprint2014arXiv

The bulge-halo conspiracy in massive elliptical galaxies: implications for the stellar initial mass function and halo response to baryonic processes

Recent studies have shown that massive elliptical galaxies have total mass density profiles within an effective radius that can be approximated as ρ_{tot}\propto r^{-γ'}, with mean slope <γ'>=2.08 \pm 0.03 and scatter σ_γ'=0.16 \pm 0.02. The small scatter of the slope (known as the bulge-halo conspiracy) is not generic in LCDM based models and therefore contains information about the galaxy formation process. We compute the distribution of γ' for LCDM-based models that reproduce the observed correlations between stellar mass, velocity dispersion, and effective radius of early-type galaxies in the SDSS. The models have a range of stellar initial mass functions (IMFs) and dark halo responses to galaxy formation. The observed distribution of γ' is well reproduced by a model with cosmologically motivated but uncontracted dark matter haloes, and a Salpeter-type IMF. Other models are on average ruled out by the data, even though they may happen in individual cases. Models with adiabatic halo contraction (and lighter IMFs) predict too small values of γ'. Models with halo expansion, or mass-follows-light predict too high values of γ'. Our study shows that the non-homologous structure of massive early-type galaxies can be precisely reproduced by LCDM models if the IMF is not universal and if mechanisms such as feedback from active galactic nuclei, or dynamical friction, effectively on average counterbalance the contraction of the halo expected as a result of baryonic cooling.

preprint2014arXiv

The InfraRed Imaging Spectrograph (IRIS) for TMT: Overview of innovative science programs

IRIS (InfraRed Imaging Spectrograph) is a first light near-infrared diffraction limited imager and integral field spectrograph being designed for the future Thirty Meter Telescope (TMT). IRIS is optimized to perform astronomical studies across a significant fraction of cosmic time, from our Solar System to distant newly formed galaxies (Barton et al. [1]). We present a selection of the innovative science cases that are unique to IRIS in the era of upcoming space and ground-based telescopes. We focus on integral field spectroscopy of directly imaged exoplanet atmospheres, probing fundamental physics in the Galactic Center, measuring 10^4 to 10^10 Msun supermassive black hole masses, resolved spectroscopy of young star-forming galaxies (1 < z < 5) and first light galaxies (6 < z < 12), and resolved spectroscopy of strong gravitational lensed sources to measure dark matter substructure. For each of these science cases we use the IRIS simulator (Wright et al. [2], Do et al. [3]) to explore IRIS capabilities. To highlight the unique IRIS capabilities, we also update the point and resolved source sensitivities for the integral field spectrograph (IFS) in all five broadband filters (Z, Y, J, H, K) for the finest spatial scale of 0.004" per spaxel. We briefly discuss future development plans for the data reduction pipeline and quicklook software for the IRIS instrument suite.

preprint2014arXiv

The Luminosity Function at z~8 from 97 Y-band dropouts: Inferences About Reionization

[Abbreviated] We present the largest search to date for $z\sim8$ Lyman break galaxies (LBGs) based on 350 arcmin$^2$ of HST observations in the V-, Y-, J- and H-bands from the Brightest of Reionizing Galaxies (BoRG) survey. The BoRG dataset includes $\sim$50 arcmin$^2$ of new data and deeper observations of two previous BoRG pointings, from which we present 9 new $z\sim8$ LBG candidates, bringing the total number of BoRG LBGs to 38 with $25.5\leqslant m_{J} \leqslant 27.6$ (AB system). We introduce a new Bayesian formalism for estimating the galaxy luminosity function (LF), which does not require binning (and thus smearing) of the data and includes a likelihood based on the formally correct binomial distribution as opposed to the often used approximate Poisson distribution. We demonstrate the utility of the new method on a sample of $97$ LBGs that combines the bright BoRG galaxies with the fainter sources published in Bouwens et al. (2012) from the HUDF and ERS programs. We show that the $z\sim8$ LF is well described by a Schechter function with a characteristic magnitude $M^\star = -20.15^{+0.29}_{-0.38}$, a faint-end slope of $α= -1.87^{+0.26}_{-0.26}$, and a number density of $\log_{10} ϕ^\star [\textrm{Mpc}^{-3}] = -3.24^{+0.25}_{-0.24}$. Integrated down to $M=-17.7$ this LF yields a luminosity density, $\log_{10} ε[\textrm{erg}/\textrm{s/Hz/Mpc}^{3}] = 25.52^{+0.05}_{-0.05}$. Our LF analysis is consistent with previously published determinations within 1$σ$. We discuss the implication of our study for the physics of reionization. By assuming theoretically motivated priors on the clumping factor and the photon escape fraction we show that the UV LF from galaxy samples down to $M=-17.7$ can ionize only 10-50% of the neutral hydrogen at $z\sim8$. Full reionization would require extending the LF down to $M=-15$.

preprint2014arXiv

The SL2S Galaxy-scale Lens Sample. V. Dark Matter Halos and Stellar IMF of Massive Early-type Galaxies out to Redshift 0.8

We investigate the cosmic evolution of the internal structure of massive early-type galaxies over half of the age of the Universe. We perform a joint lensing and stellar dynamics analysis of a sample of 81 strong lenses from the SL2S and SLACS surveys and combine the results with a hierarchical Bayesian inference method to measure the distribution of dark matter mass and stellar IMF across the population of massive early-type galaxies. Lensing selection effects are taken into account. We find that the dark matter mass projected within the inner 5 kpc increases for increasing redshift, decreases for increasing stellar mass density, but is roughly constant along the evolutionary tracks of early-type galaxies. The average dark matter slope is consistent with that of an NFW profile, but is not well constrained. The stellar IMF normalization is close to a Salpeter IMF at $\log{M_*} = 11.5$ and scales strongly with increasing stellar mass. No dependence of the IMF on redshift or stellar mass density is detected. The anti-correlation between dark matter mass and stellar mass density supports the idea of mergers being more frequent in more massive dark matter halos.

preprint2014arXiv

The stellar initial mass function of early type galaxies from low to high stellar velocity dispersion: homogeneous analysis of ATLAS$^{\rm 3D}$ and Sloan Lens ACS galaxies

We present an investigation about the shape of the initial mass function (IMF) of early-type galaxies (ETGs), based on a joint lensing and dynamical analysis, and on stellar population synthesis models, for a sample of 55 lens ETGs identified by the Sloan Lens ACS (SLACS) Survey. We construct axisymmetric dynamical models based on the Jeans equations which allow for orbital anisotropy and include a dark matter halo. The models reproduce in detail the observed \textit{HST} photometry and are constrained by the total projected mass within the Einstein radius and the stellar velocity dispersion ($σ$) within the SDSS fibers. Comparing the dynamically-derived stellar mass-to-light ratios $(M_*/L)_{\rm dyn}$, obtained for an assumed halo slope $ρ_{\rm h}\propto r^{-1}$, to the stellar population ones $(M_*/L)_{\rm pop}$, derived from full-spectrum fitting and assuming a Salpeter IMF, we infer the mass normalization of the IMF. Our results confirm the previous analysis by the SLACS team that the mass normalization of the IMF of high $σ$ galaxies is consistent on average with a Salpeter slope. Our study allows for a fully consistent study of the trend between IMF and $σ$ for both the SLACS and \ATLAS samples, which explore quite different $σ$ ranges. The two samples are highly complementary, the first being essentially $σ$ selected, and the latter volume-limited and nearly mass selected. We find that the two samples merge smoothly into a single trend of the form $\logα=(0.38\pm0.04)\times\log(σ_{\rm e}/200\,\mathrm{km~s}^{-1})+(-0.06\pm0.01)$, where $α=(M_*/L)_{\rm dyn}/(M_*/L)_{\rm pop}$ and $σ_{\rm e}$ is the luminosity averaged $σ$ within one effective radius $R_{\rm e}$. This is consistent with a systematic variation of the IMF normalization from Kroupa to Salpeter in the interval $σ_{\rm e}\approx90-270\,\mathrm{km~s}^{-1}$.

preprint2014arXiv

X-ray constraints on the local supermassive black hole occupation fraction

Distinct seed formation mechanisms are imprinted upon the fraction of dwarf galaxies currently containing a central supermassive black hole. Seeding by Pop III remnants is expected to produce a higher occupation fraction than is generated with direct gas collapse precursors. Chandra observations of nearby early-type galaxies can directly detect even low-level supermassive black hole activity, and the active fraction immediately provides a firm lower limit to the occupation fraction. Here, we use the volume-limited AMUSE surveys of ~200 optically-selected early-type galaxies to characterize simultaneously, for the first time, the occupation fraction and the scaling of nuclear X-ray luminosity with stellar mass, accounting for intrinsic scatter, measurement uncertainties, and X-ray limits. For early-type galaxies with log(M_star/M_sun)<10, we obtain a lower limit to the occupation fraction of >20% (at 95% confidence), but full occupation cannot be excluded. The preferred dependence of log(L_X) upon log(M_star) has a slope of about 0.7-0.8, consistent with the "downsizing" trend previously identified from the AMUSE dataset, and a uniform Eddington efficiency is disfavored at ~2 sigma. We provide guidelines for the future precision with which these parameters may be refined with larger or more sensitive samples.

preprint2013arXiv

A low-mass cut-off near the hydrogen burning limit for Salpeter-like initial mass functions in early-type galaxies

We conduct a detailed investigation of the properties of the stellar initial mass function (IMF) in two massive early-type lens galaxies with velocity dispersions of sigma ~245 km/s and sigma ~325 km/s, for which both HST imaging and X-Shooter spectra are available. We compare the inferences obtained from two fully independent methods: (i) a combined gravitational lensing and stellar dynamics (L&D) analysis of the data sets employing self-consistent axisymmetric models, and (ii) a spectroscopic simple stellar population (SSP) analysis of optical line-strength indices, assuming single power-law IMFs. The results from the two approaches are found to be in agreement within the 1-sigma uncertainties. Both galaxies are consistent with having a Salpeter IMF (power-law slope of x = 2.35), which is strongly favoured over a Chabrier IMF (x = 1.8), with probabilities inferred from the joint analysis of 89% and 99%, respectively. Bottom-heavy IMFs significantly steeper than Salpeter (x >= 3.0) are ruled out with decisive evidence (Bayes factor B > 1000) for both galaxies, as they exceed the total mass derived from the L&D constraints. Our analysis allows, for the first time, the inference of the low-mass cut-off of the IMF (M_low). Combining the joint L&D and SSP analyses of both galaxies, we infer an IMF slope of x = 2.22 +/- 0.14, consistent with Salpeter IMF, and a low-mass limit M_low = 0.13 +/- 0.03 M_sun, just above the hydrogen burning limit.

preprint2013arXiv

An Environmental Study of the Ultraluminous X-ray Source Population in Early-type Galaxies

Ultraluminous X-ray sources (ULXs) are some of the brightest phenomena found outside of a galaxy's nucleus, and their explanation typically invokes accretion of material onto a black hole. Here, we perform the largest population study to date of ULXs in early-type galaxies, focusing on whether a galaxy's large scale environment can affect its ULX content. Using the AMUSE survey, which includes homogeneous X-ray coverage of 100 elliptical galaxies in the Virgo cluster and a similar number of elliptical galaxies in the field (spanning stellar masses of 1e8-1e12 M_sun), we identify 37.9 +/- 10.1 ULXs in Virgo and 28.1 +/- 8.7 ULXs in the field. Across both samples, we constrain the number of ULXs per unit stellar mass, i.e., the ULX specific frequency, to be 0.062 +/- 0.013 ULXs per 10^10 M_sun (or about 1 ULX per 1.6e11 M_sun of galaxy stellar mass). We find that the number of ULXs, the specific frequency of ULXs, and the average ULX spectral properties are all similar in both cluster and field environments. Contrary to late-type galaxies, we do not see any trend between specific ULX frequency and host galaxy stellar mass, and we show that dwarf ellipticals host fewer ULXs than later-type dwarf galaxies at a statistically meaningful level. Our results are consistent with ULXs in early-type galaxies probing the luminous tail of the low-mass X-ray binary population, and are briefly discussed in context of the influence of gravitational interactions on the long-term evolution of a galaxy's (older) stellar population.

preprint2013arXiv

Improving the precision of time-delay cosmography with observations of galaxies along the line of sight

In order to use strong gravitational lens time delays to measure precise and accurate cosmological parameters the effects of mass along the line of sight must be taken into account. We present a method to achieve this by constraining the probability distribution function of the effective line of sight convergence k_ext. The method is based on matching the observed overdensity in the weighted number of galaxies to that found in mock catalogs with k_ext obtained by ray-tracing through structure formation simulations. We explore weighting schemes based on projected distance, mass, luminosity, and redshift. This additional information reduces the uncertainty of k_ext from sigma_k $0.06 to ~0.04 for very overdense lines of sight like that of the system B1608+656. For more common lines of sight, sigma_k is reduced to ~<0.03, corresponding to an uncertainty of ~<3% on distance. This uncertainty has comparable effects on cosmological parameters to that arising from the mass model of the deflector and its immediate environment. Photometric redshifts based on g, r, i and K photometries are sufficient to constrain k_ext almost as well as with spectroscopic redshifts. As an illustration, we apply our method to the system B1608+656. Our most reliable k_ext estimator gives sigma_k=0.047 down from 0.065 using only galaxy counts. Although deeper multi-band observations of the field of B1608+656 are necessary to obtain a more precise estimate, we conclude that griK photometry, in addition to spectroscopy to characterize the immediate environment, is an effective way to increase the precision of time-delay cosmography.

preprint2013arXiv

Quasar Black Hole Mass Estimates in the Era of Time Domain Astronomy

We investigate the dependence of the normalization of the high-frequency part of the X-ray and optical power spectral densities (PSD) on black hole mass for a sample of 39 active galactic nuclei (AGN) with black hole masses estimated from reverberation mapping or dynamical modeling. We obtained new Swift observations of PG 1426+015, which has the largest estimated black hole mass of the AGN in our sample. We develop a novel statistical method to estimate the PSD from a lightcurve of photon counts with arbitrary sampling, eliminating the need to bin a lightcurve to achieve Gaussian statistics, and we use this technique to estimate the X-ray variability parameters for the faint AGN in our sample. We find that the normalization of the high-frequency X-ray PSD is inversely proportional to black hole mass. We discuss how to use this scaling relationship to obtain black hole mass estimates from the short time-scale X-ray variability amplitude with precision ~ 0.38 dex. The amplitude of optical variability on time scales of days is also anti-correlated with black hole mass, but with larger scatter. Instead, the optical variability amplitude exhibits the strongest anti-correlation with luminosity. We conclude with a discussion of the implications of our results for estimating black hole mass from the amplitude of AGN variability.

preprint2013arXiv

Reconstructing the Lensing Mass in the Universe from Photometric Catalogue Data

High precision cosmological distance measurements towards individual objects such as time delay gravitational lenses or type Ia supernovae are affected by weak lensing perturbations by galaxies and groups along the line of sight. In time delay gravitational lenses, "external convergence," kappa, can dominate the uncertainty in the inferred distances and hence cosmological parameters. In this paper we attempt to reconstruct kappa, due to line of sight structure, using a simple halo model. We use mock catalogues from the Millennium Simulation, and calibrate and compare our reconstructed P(kappa) to ray-traced kappa "truth" values; taking into account realistic observational uncertainties. We find that the reconstruction of kappa provides an improvement in precision of ~50% over galaxy number counts. We find that the lowest-kappa lines of sight have the best constrained P(kappa). In anticipation of large future samples of lenses, we find that selecting the third of the systems with the highest precision kappa estimates gives a sample of unbiased time delay distance measurements with just ~1% uncertainty due to line of sight external convergence effects. Photometric data are sufficient to pre-select the best-constrained lines of sight, and can be done before investment in light-curve monitoring. Conversely, we show that selecting lines of sight with high external shear could, with the reconstruction model presented, induce biases of up to 1% in time delay distance. We find that a major potential source of systematic error is uncertainty in the high mass end of the stellar mass-halo mass relation; this could introduce ~2% biases on the time-delay distance if completely ignored. We suggest areas for the improvement of this general analysis framework (including more sophisticated treatment of high mass structures) that should allow yet more accurate cosmological inferences to be made.

preprint2013arXiv

The Density Profiles of Massive, Relaxed Galaxy Clusters. I. The Total Density Over Three Decades in Radius

Clusters of galaxies are excellent locations to probe the distribution of baryons and dark matter (DM) over a wide range of scales. We study a sample of seven massive, relaxed galaxy clusters with centrally-located brightest cluster galaxies (BCGs) at z=0.2-0.3. Using the observational tools of strong and weak gravitational lensing, combined with resolved stellar kinematics within the BCG, we measure the total radial density profile, comprising both dark and baryonic matter, over scales of ~3-3000 kpc. Lensing-derived mass profiles typically agree with independent X-ray estimates within ~15%, suggesting that departures from hydrostatic equilibrium are small and that the clusters in our sample (except A383) are not strongly elongated along the line of sight. The inner logarithmic slope gamma_tot of the total density profile measured over r/r200=0.003-0.03, where rho_tot ~ r^(-gamma_tot), is found to be nearly universal, with a mean <gamma_tot> = 1.16 +- 0.05 (random) +0.05-0.07 (systematic) and an intrinsic scatter of < 0.13 (68% confidence). This is further supported by the very homogeneous shape of the observed velocity dispersion profiles, obtained via Keck spectroscopy, which are mutually consistent after a simple scaling. Remarkably, this slope agrees closely with numerical simulations that contain only dark matter, despite the significant contribution of stellar mass on the scales we probe. The Navarro-Frenk-White profile characteristic of collisionless cold dark matter is a better description of the total mass density at radii >~ 5-10 kpc than that of dark matter alone. Hydrodynamical simulations that include baryons, cooling, and feedback currently provide a poorer match. We discuss the significance of our findings for understanding the assembly of BCGs and cluster cores, particularly the influence of baryons on the inner DM halo. [abridged]

preprint2013arXiv

The Density Profiles of Massive, Relaxed Galaxy Clusters. II. Separating Luminous and Dark Matter in Cluster Cores

We present stellar and dark matter (DM) density profiles for a sample of seven massive, relaxed galaxy clusters derived from strong and weak gravitational lensing and resolved stellar kinematic observations within the centrally-located brightest cluster galaxies (BCGs). In Paper I of the series, we demonstrated that the total density profile derived from these data, which span 3 decades in radius, is consistent with numerical DM-only simulations at radii >~ 5-10 kpc, despite the significant contribution of stellar material in the core. Here we decompose the inner mass profiles of these clusters into stellar and dark components. Parametrizing the DM density profile as a power law rho_DM ~ r^{-β} on small scales, we find a mean slope <β> = 0.50 +- 0.10 (random) +0.14-0.13 (systematic). Alternatively, cored Navarro-Frenk-White (NFW) profiles with <log r_core/kpc> = 1.14 +- 0.13 (random) +0.14-0.22 (systematic) provide an equally good description. These density profiles are significantly shallower than canonical NFW models at radii <~ 30 kpc, comparable to the effective radii of the BCGs. The inner DM profile is correlated with the distribution of stars in the BCG, suggesting a connection between the inner halo and the assembly of stars in the central galaxy. The stellar mass-to-light ratio inferred from lensing and stellar dynamics is consistent with that inferred using stellar population synthesis models if a Salpeter initial mass function is adopted. We compare these results to theories describing the interaction between baryons and DM in cluster cores, including adiabatic contraction models and the possible effects of galaxy mergers and active galactic nucleus feedback, and evaluate possible signatures of alternative DM candidates.

preprint2013arXiv

The Lick AGN Monitoring Project 2011: Fe II Reverberation from the Outer Broad-Line Region

The prominent broad Fe II emission blends in the spectra of active galactic nuclei have been shown to vary in response to continuum variations, but past attempts to measure the reverberation lag time of the optical Fe II lines have met with only limited success. Here we report the detection of Fe II reverberation in two Seyfert 1 galaxies, NGC 4593 and Mrk 1511, based on data from a program carried out at Lick Observatory in Spring 2011. Light curves for emission lines including H-beta and Fe II were measured by applying a fitting routine to decompose the spectra into several continuum and emission-line components, and we use cross-correlation techniques to determine the reverberation lags of the emission lines relative to V-band light curves. In both cases the measured lag (t_cen) of Fe II is longer than that of H-beta, although the inferred lags are somewhat sensitive to the choice of Fe II template used in the fit. For spectral decompositions done using the Fe II template of Veron-Cetty et al. (2004), we find t_cen(Fe II)/t_cen(H-beta) = 1.9+-0.6 in NGC 4593 and 1.5+-0.3 in Mrk 1511. The detection of highly correlated variations between Fe II and continuum emission demonstrates that the Fe II emission in these galaxies originates in photoionized gas, located predominantly in the outer portion of the broad-line region.

preprint2013arXiv

The Low-Luminosity End of the Radius-Luminosity Relationship for Active Galactic Nuclei

We present an updated and revised analysis of the relationship between the Hbeta broad-line region (BLR) radius and the luminosity of the active galactic nucleus (AGN). Specifically, we have carried out two-dimensional surface brightness decompositions of the host galaxies of 9 new AGNs imaged with the Hubble Space Telescope Wide Field Camera 3. The surface brightness decompositions allow us to create "AGN-free" images of the galaxies, from which we measure the starlight contribution to the optical luminosity measured through the ground-based spectroscopic aperture. We also incorporate 20 new reverberation-mapping measurements of the Hbeta time lag, which is assumed to yield the average Hbeta BLR radius. The final sample includes 41 AGNs covering four orders of magnitude in luminosity. The additions and updates incorporated here primarily affect the low-luminosity end of the R-L relationship. The best fit to the relationship using a Bayesian analysis finds a slope of alpha = 0.533 (+0.035/-0.033), consistent with previous work and with simple photoionization arguments. Only two AGNs appear to be outliers from the relationship, but both of them have monitoring light curves that raise doubt regarding the accuracy of their reported time lags. The scatter around the relationship is found to be 0.19(+/-0.02) dex, but would be decreased to 0.13 dex by the removal of these two suspect measurements. A large fraction of the remaining scatter in the relationship is likely due to the inaccurate distances to the AGN host galaxies. Our results help support the possibility that the R-L relationship could potentially be used to turn the BLRs of AGNs into standardizable candles. This would allow the cosmological expansion of the Universe to be probed by a separate population of objects, and over a larger range of redshifts.

preprint2013arXiv

The SL2S Galaxy-scale Lens Sample. III. Lens Models, Surface Photometry and Stellar Masses for the final sample

We present Hubble Space Telescope (HST) imaging data and CFHT Near IR ground-based images for the final sample of 56 candidate galaxy-scale lenses uncovered in the CFHT Legacy Survey as part of the Strong Lensing in the Legacy Survey (SL2S) project. The new images are used to perform lens modeling, measure surface photometry, and estimate stellar masses of the deflector early-type galaxies. Lens modeling is performed on the HST images (or CFHT when HST is not available) by fitting the spatially extended light distribution of the lensed features assuming a singular isothermal ellipsoid mass profile and by reconstructing the intrinsic source light distribution on a pixelized grid. Based on the analysis of systematic uncertainties and comparison with inference based on different methods we estimate that our Einstein Radii are accurate to \sim3%. HST imaging provides a much higher success rate in confirming gravitational lenses and measuring their Einstein radii than CFHT imaging does. Lens modeling with ground-based images however, when successful, yields Einstein radius measurements that are competitive with spaced-based images. Information from the lens models is used together with spectroscopic information from the companion paper IV to classify the systems, resulting in a final sample of 39 confirmed (grade-A) lenses and 17 promising candidates. The redshifts of the main deflector span a range 0.3<zd< 0.8, providing an excellent sample for the study of the cosmic evolution of the mass distribution of early-type galaxies over the second half of the history of the Universe.

preprint2013arXiv

The SL2S Galaxy-scale Lens Sample. IV. The dependence of the total mass density profile of early-type galaxies on redshift, stellar mass, and size

We present optical and near infrared spectroscopy obtained at Keck, VLT, and Gemini for a sample of 36 secure strong gravitational lens systems and 17 candidates identified as part of the SL2S survey. The deflectors are massive early-type galaxies in the redshift range z_d=0.2-0.8, while the lensed sources are at z_s=1-3.5. We combine this data with photometric and lensing measurements presented in the companion paper III and with lenses from the SLACS and LSD surveys to investigate the cosmic evolution of the internal structure of massive early-type galaxies over half the age of the universe. We study the dependence of the slope of the total mass density profile γ' (ρ(r)\propto r^{-γ'}) on stellar mass, size, and redshift. We find that two parameters are sufficent to determine γ' with less than 6% residual scatter. At fixed redshift, γ' depends solely on the surface stellar mass density \partial γ'/ \partial Σ_*=0.38\pm 0.07, i.e. galaxies with denser stars also have steeper slopes. At fixed M_* and R_{eff}, γ' depends on redshift, in the sense that galaxies at a lower redshift have steeper slopes (\partial γ' / \partial z = -0.31\pm 0.10). However, the mean redshift evolution of γ' for an individual galaxy is consistent with zero dγ'/dz=-0.10\pm0.12. This result is obtained by combining our measured dependencies of γ' on z,M_*,R_{eff} with the evolution of the R_{eff}-M_* taken from the literature, and is broadly consistent with current models of the formation and evolution of massive early-type galaxies. Detailed quantitative comparisons of our results with theory will provide qualitatively new information on the detailed physical processes at work.

preprint2013arXiv

The SWELLS Survey. VI. hierarchical inference of the initial mass functions of bulges and discs

The long-standing assumption that the stellar initial mass function (IMF) is universal has recently been challenged by a number of observations. Several studies have shown that a "heavy" IMF (e.g., with a Salpeter-like abundance of low mass stars and thus normalisation) is preferred for massive early-type galaxies, while this IMF is inconsistent with the properties of less massive, later-type galaxies. These discoveries motivate the hypothesis that the IMF may vary (possibly very slightly) across galaxies and across components of individual galaxies (e.g. bulges vs discs). In this paper we use a sample of 19 late-type strong gravitational lenses from the SWELLS survey to investigate the IMFs of the bulges and discs in late-type galaxies. We perform a joint analysis of the galaxies' total masses (constrained by strong gravitational lensing) and stellar masses (constrained by optical and near-infrared colours in the context of a stellar population synthesis [SPS] model, up to an IMF normalisation parameter). Using minimal assumptions apart from the physical constraint that the total stellar mass within any aperture must be less than the total mass within the aperture, we find that the bulges of the galaxies cannot have IMFs heavier (i.e. implying high mass per unit luminosity) than Salpeter, while the disc IMFs are not well constrained by this data set. We also discuss the necessity for hierarchical modelling when combining incomplete information about multiple astronomical objects. This modelling approach allows us to place upper limits on the size of any departures from universality. More data, including spatially resolved kinematics (as in paper V) and stellar population diagnostics over a range of bulge and disc masses, are needed to robustly quantify how the IMF varies within galaxies.

preprint2012arXiv

AMUSE-Field I: Nuclear X-ray Properties of Local Field and Group Spheroids across the Stellar Mass Scale

We present the first results from AMUSE-Field, a Chandra survey designed to characterize the occurrence and intensity of low-level accretion onto supermassive black holes (SMBHs) at the center of local early-type field galaxies. This is accomplished by means of a Large Program targeting a distance-limited (<30 Mpc) sample of 103 early types spanning a wide range in stellar masses. We acquired new ACIS-S observations for 61 objects down to a limiting (0.3-10 keV) luminosity of 2.5x10^38 erg/s, and we include an additional 42 objects with archival (typically deeper) coverage. A nuclear X-ray source is detected in 52 out of the 103 galaxies. After accounting for potential contamination from low-mass X-ray binaries, we estimate that the fraction of accreting SMBHs within the sample is 45+/-7 percent, which sets a firm lower limit on the occupation fraction within the field. The measured nuclear X-ray luminosities are invariably highly sub-Eddington, with -8<log(L_X/L_Edd)<-4. As also found in a companion survey targeting Virgo early types, the active fraction increases with increasing host galaxy stellar mass, reflective of "Eddington incompleteness" within the lower-mass objects. For the Field sample, the average nuclear X-ray luminosity scales with the host stellar mass as M_star^(0.71+/-0.10), with an intrinsic scatter of 0.73+/-0.09 dex. A majority of the AMUSE-Field galaxies (78%) inhabits groups, enabling us to investigate the influence of group richness upon nuclear activity. We see no evidence for a positive correlation between nuclear X-ray luminosity, normalized to host properties, and galaxy density. Rather, while the scatter is substantial, it appears that the Eddington-scaled X-ray luminosity of group members may be slightly lower than for isolated galaxies, and that this trend continues to cluster early-types.

preprint2012arXiv

Can Minor Merging Account for the Size Growth of Quiescent Galaxies? New Results from the CANDELS Survey

The presence of extremely compact galaxies at z~2 and their subsequent growth in physical size has been the cause of much puzzlement. We revisit the question using deep infrared Wide Field Camera 3 data to probe the rest-frame optical structure of 935 host galaxies selected with 0.4<z<2.5 and stellar masses M* > 10^10.7 Msol using optical and near-infrared photometry in the UKIRT Ultra Deep Survey and GOODS-South fields of the CANDELS survey. At each redshift, the most compact sources are those with little or no star formation, and we find that the mean size of these systems grows by a factor of 3.5 +- 0.3 over this redshift interval. The new data are sufficiently deep to enable us to identify companions to these hosts whose stellar masses are ten times smaller, while still yielding suitably accurate photometric redshifts to define a likely physical association. By searching for faint companions around 404 quiescent hosts within a projected physical annulus 10 < R < 30 kpc/h, we estimate the minor merger rate over the redshift range 0.4 < z < 2. After correcting for contamination from projected pairs, we find that 13-18% of quiescent hosts have likely physical companions with stellar mass ratios of 0.1 or greater. Mergers of these companions will typically increase the host mass by 6+-2% per merger timescale. We estimate the minimum growth rate necessary to explain the declining abundance of compact galaxies. Using a simple model of merging motivated by recent numerical simulations, we then assess whether mergers of the faint companions with their hosts are sufficient to explain this minimal rate. We find that mergers with mass ratios > 0.1 may explain most of the size evolution observed at z >~ 1 if a relatively short merger timescale is assumed, but the rapid growth seen at higher redshift likely requires additional physical processes.

preprint2012arXiv

Evidence for dark matter contraction and a Salpeter IMF in a massive early-type galaxy

Stars and dark matter account for most of the mass of early-type galaxies, but uncertainties in the stellar population and the dark matter profile make it challenging to distinguish between the two components. Nevertheless, precise observations of stellar and dark matter are extremely valuable for testing the many models of structure formation and evolution. We present a measurement of the stellar mass and inner slope of the dark matter halo of a massive early-type galaxy at $z=0.222$. The galaxy is the foreground deflector of the double Einstein ring gravitational lens system SDSSJ0946+1006, also known as the "Jackpot". By combining the tools of lensing and dynamics, we first constrain the mean slope of the total mass density profile ($ρ_{\rm{tot}}\propto r^{-γ'}$) within the radius of the outer ring to be $γ' = 1.98\pm0.02\pm0.01$. Then we obtain a bulge-halo decomposition, assuming a power-law form for the dark matter halo. Our analysis yields $γ_{\rm{DM}} = 1.7\pm0.2$ for the inner slope of the dark matter profile, in agreement with theoretical findings on the distribution of dark matter in ellipticals, and a stellar mass from lensing and dynamics $M_*^{\rm{LD}} = 5.5_{-1.3}^{+0.4}\times10^{11}M_\Sun$. By comparing this measurement with stellar masses inferred from stellar population synthesis fitting we find that a Salpeter IMF provides a good description of the stellar population of the lens while the probability of the IMF being heavier than Chabrier is 95%. Our data suggest that growth by accretion of small systems from a compact red nugget is a plausible formation scenario for this object.

preprint2012arXiv

Recalibration of the virial factor and M-sigma relation for local active galaxies

Determining the virial factor of the broad-line region (BLR) gas is crucial for calibrating AGN black hole mass estimators, since the measured line-of-sight velocity needs to be converted into the intrinsic virial velocity. The average virial factor has been empirically calibrated based on the M-sigma relation of quiescent galaxies, but the claimed values differ by a factor of two in recent studies. We investigate the origin of the difference by measuring the M-sigma relation using an updated galaxy sample from the literature, and explore the dependence of the virial factor on various fitting methods. We find that the discrepancy is primarily caused by the sample selection, while the difference stemming from the various regression methods is marginal. However, we generally prefer the FITEXY and Bayesian estimators based on Monte Carlo simulations for the M-sigma relation. In addition, the choice of independent variable in the regression leads to ~0.2 dex variation in the virial factor inferred from the calibration process. Based on the determined virial factor, we present the updated M-sigma relation of local active galaxies.

preprint2012arXiv

Spectroscopic Confirmation of a z=6.740 Galaxy behind the Bullet Cluster

We present the first results of our spectroscopic follow-up of 6.5 < z < 10 candidate galaxies behind clusters of galaxies. We report the spectroscopic confirmation of an intrinsically faint Lyman break galaxy (LBG) identified as a z 850LP-band dropout behind the Bullet Cluster. We detect an emission line at λ = 9412 Å at >5σ significance using a 16 hr long exposure with FORS2 VLT. Based on the absence of flux in bluer broadband filters, the blue color of the source, and the absence of additional lines, we identify the line as Lyα at z = 6.740 \pm 0.003. The integrated line flux is f = (0.7 \pm 0.1 \pm 0.3) \times 10^{-17} erg^{-1} s^{-1} cm^{-2} (the uncertainties are due to random and flux calibration errors, respectively) making it the faintest Lyα flux detected at these redshifts. Given the magnification of μ = 3.0 \pm 0.2 the intrinsic (corrected for lensing) flux is f^int = (0.23 \pm 0.03 \pm 0.10 \pm 0.02) \times 10^{-17} erg^{-1} s^{-1} cm^{-2} (additional uncertainty due to magnification), which is ~2-3 times fainter than other such measurements in z ~ 7 galaxies. The intrinsic H 160W-band magnitude of the object is m^int(H_160W)=27.57 \pm 0.17, corresponding to 0.5 L* for LBGs at these redshifts. The galaxy is one of the two sub-L* LBG galaxies spectroscopically confirmed at these high redshifts (the other is also a lensed z = 7.045 galaxy), making it a valuable probe for the neutral hydrogen fraction in the early universe.

preprint2012arXiv

The Assembly History of Disk Galaxies: II. Probing the Emerging Tully-Fisher Relation During 1<z<1.7

Through extended integrations using the recently-installed deep depletion CCD on the red arm of the Keck I Low Resolution Imaging Spectrograph, we present new measurements of the resolved spectra of 70 morphologically-selected star-forming galaxies with i_AB<24.1 in the redshift range 1<z<1.7. Using the formalism introduced in Paper I of this series and available HST ACS images, we successfully recover rotation curves using the extended emission line distribution of [O II] 3727 A to 2.2 times the disk scale radius for a sample of 42 galaxies. Combining these measures with stellar masses derived from HST and ground-based near-infrared photometry enables us to construct the stellar mass Tully-Fisher relation in the time interval between the well-constructed relation defined at z~1 in Paper I and the growing body of resolved dynamics probed with integral field unit spectrographs at z>2. Remarkably, we find a well-defined Tully-Fisher relation with up to 60% increase in scatter and stellar mass zero-point shift constraint of 0.02+/-0.02 dex since z~1.7, compared to the local relation. Although our sample is incomplete in terms of either a fixed stellar mass or star formation rate limit, we discuss the implications that typical star-forming disk galaxies evolve to arrive on a well-defined Tully-Fisher relation within a surprisingly short period of cosmic history.

preprint2012arXiv

The Lick AGN Monitoring Project 2011: Dynamical Modeling of the Broad Line Region in Mrk 50

We present dynamical modeling of the broad line region (BLR) in the Seyfert 1 galaxy Mrk 50 using reverberation mapping data taken as part of the Lick AGN Monitoring Project (LAMP) 2011. We model the reverberation mapping data directly, constraining the geometry and kinematics of the BLR, as well as deriving a black hole mass estimate that does not depend on a normalizing factor or virial coefficient. We find that the geometry of the BLR in Mrk 50 is a nearly face-on thick disk, with a mean radius of 9.6(+1.2,-0.9) light days, a width of the BLR of 6.9(+1.2,-1.1) light days, and a disk opening angle of 25\pm10 degrees above the plane. We also constrain the inclination angle to be 9(+7,-5) degrees, close to face-on. Finally, the black hole mass of Mrk 50 is inferred to be log10(M(BH)/Msun) = 7.57(+0.44,-0.27). By comparison to the virial black hole mass estimate from traditional reverberation mapping analysis, we find the normalizing constant (virial coefficient) to be log10(f) = 0.78(+0.44,-0.27), consistent with the commonly adopted mean value of 0.74 based on aligning the M(BH)-σ* relation for AGN and quiescent galaxies. While our dynamical model includes the possibility of a net inflow or outflow in the BLR, we cannot distinguish between these two scenarios.

preprint2012arXiv

The role of environment in low-level active galactic nucleus activity: no evidence for cluster enhancement

We use the AMUSE-Virgo and AMUSE-Field surveys for nuclear X-ray emission in early-type galaxies to conduct a controlled comparison of low-level supermassive black hole activity within cluster and field spheroids. While both the Virgo and the Field samples feature highly sub-Eddington X-ray luminosities (log L_x/L_Edd between -8 and -4), we find that after accounting for the influence of host galaxy stellar mass, the field early-type galaxies tend toward marginally greater (0.38+/-0.14 dex) nuclear X-ray luminosities, at a given black hole mass, than their cluster counterparts. This trend is qualitatively consistent with the field black holes having access to a greater reservoir of fuel, plausibly in the form of cold gas located near the nucleus. We are able to rule out at high confidence the alternative of enhanced X-ray activity within clusters. Presuming nuclear X-ray emission correlates with the total energy and momentum output of these weakly accreting black holes, this indicates that low-level active galactic nucleus feedback is not generally stronger within typical cluster galaxies than in the field. These results confirm that for most cluster early-type galaxies (i.e., excluding brightest cluster galaxies) direct environmental effects, such as gas stripping, are more relevant in quenching star formation.

preprint2012arXiv

The SL2S Galaxy-scale Gravitational Lens Sample. I. The alignment of mass and light in massive early-type galaxies at z=0.2-0.9

We study the relative alignment of mass and light in a sample of 16 massive early-type galaxies at z=0.2-0.9 that act as strong gravitational lenses. The sample was identified from deep multi-band images obtained as part of the Canada France Hawaii Telescope Legacy Survey as part of the Strong Lensing Legacy Survey (SL2S). Higher resolution follow-up imaging is available for a subset of 10 systems. We construct gravitational lens models and infer total enclosed mass, elongation, and position angle of the mass distribution. By comparison with the observed distribution of light we infer that there is a substantial amount of external shear with mean value $<γ_{\rm ext}>= 0.12 \pm 0.05$, arising most likely from the environment of the SL2S lenses. In a companion paper, we combine these measurements with follow-up Keck spectroscopy to study the evolution of the stellar and dark matter content of early-type galaxies as a function of cosmic time.

preprint2012arXiv

The SWELLS survey. III. Disfavouring "heavy" initial mass functions for spiral lens galaxies

We present gravitational lens models for 20 strong gravitational lens systems observed as part of the Sloan WFC Edge-on Late-type Lens Survey (SWELLS) project. Fifteen of the lenses are taken from paper I while five are newly discovered systems. The systems are galaxy-galaxy lenses where the foreground deflector has an inclined disc, with a wide range of morphological types, from late-type spiral to lenticular. For each system, we compare the total mass inside the critical curve inferred from gravitational lens modelling to the stellar mass inferred from stellar population synthesis (SPS) models, computing the stellar mass fraction f* = M(SPS)/M(lens). We find that, for the lower mass SWELLS systems, adoption of a Salpeter stellar initial mass function (IMF) leads to estimates of f* that exceed 1. This is unphysical, and provides strong evidence against the Salpeter IMF being valid for these systems. Taking the lower mass end of the SWELLS sample sigma(SIE) < 230 km/s, we find that the IMF is lighter (in terms of stellar mass-to-light ratio) than Salpeter with 98% probability, and consistent with the Chabrier IMF and IMFs between the two. This result is consistent with previous studies of spiral galaxies based on independent techniques. In combination with recent studies of massive early-type galaxies that have favoured a heavier Salpeter-like IMF, this result strengthens the evidence against a universal stellar IMF.

preprint2012arXiv

The SWELLS survey. IV. Precision measurements of the stellar and dark matter distributions in a spiral lens galaxy

We construct a fully self-consistent mass model for the lens galaxy J2141 at z=0.14, and use it to improve on previous studies by modelling its gravitational lensing effect, gas rotation curve and stellar kinematics simultaneously. We adopt a very flexible axisymmetric mass model constituted by a generalized NFW dark matter halo and a stellar mass distribution obtained by deprojecting the MGE fit to the high-resolution K'-band LGSAO imaging data of the galaxy, with the (spatially constant) M/L ratio as a free parameter. We model the stellar kinematics by solving the anisotropic Jeans equations. We find that the inner logarithmic slope of the dark halo is weakly constrained (gamma = 0.82^{+0.65}_{-0.54}), and consistent with an unmodified NFW profile. We infer the galaxy to have (i) a dark matter fraction within 2.2 disk radii of 0.28^{+0.15}_{-0.10}, independent of the galaxy stellar population, implying a maximal disk for J2141; (ii) an apparently uncontracted dark matter halo, with concentration c_{-2} = 7.7_{-2.5}^{+4.2} and virial velocity v_{vir} = 242_{-39}^{+44} km/s, consistent with LCDM predictions; (iii) a slightly oblate halo (q_h = 0.75^{+0.27}_{-0.16}), consistent with predictions from baryon-affected models. Comparing the stellar mass inferred from the combined analysis (log_{10} Mstar/Msun = 11.12_{-0.09}^{+0.05}) with that inferred from SPS modelling of the galaxies colours, and accounting for a cold gas fraction of 20+/-10%, we determine a preference for a Chabrier IMF over Salpeter IMF by a Bayes factor of 5.7 (substantial evidence). We infer a value beta_{z} = 1 - sigma^2_{z}/sigma^2_{R} = 0.43_{-0.11}^{+0.08} for the orbital anisotropy parameter in the meridional plane, in agreement with most studies of local disk galaxies, and ruling out at 99% CL that the dynamics of this system can be described by a two-integral distribution function. [Abridged]

preprint2012arXiv

Using the Bullet Cluster as a Gravitational Telescope to Study z~7 Lyman Break Galaxies

We use imaging obtained with the Hubble Space Telescope Wide Field Camera 3 to search for z_850 dropouts at z~7 and J_110 dropouts at z~9 lensed by the Bullet Cluster. In total we find 10 z_850 dropouts in our 8.27 arcmin^2 field. Using magnification maps from a combined weak and strong lensing mass reconstruction of the Bullet Cluster and correcting for estimated completeness levels, we calculate the surface density and luminosity function of our z_850 dropouts as a function of intrinsic (accounting for magnification) magnitude. We find results consistent with published blank field surveys, despite using much shallower data, and demonstrate the effectiveness of cluster surveys in the search for z~7 galaxies.

preprint2011arXiv

A Local Baseline of the Black Hole Mass Scaling Relations for Active Galaxies. I. Methodology and Results of Pilot Study

We present high-quality Keck/LRIS longslit spectroscopy of a pilot sample of 25 local active galaxies selected from the SDSS (0.02<z<0.1; MBH>10^7 M_sun) to study the relations between black hole mass (MBH) and host-galaxy properties. We determine stellar kinematics of the host galaxy, deriving stellar-velocity dispersion profiles and rotation curves from three spectral regions (including CaH&K, MgIb triplet, and CaII triplet). In addition, we perform surface photometry on SDSS images, using a newly developed code for joint multi-band analysis. BH masses are estimated from the width of the Hbeta emission line and the host-galaxy free 5100A AGN luminosity. Combining results from spectroscopy and imaging allows us to study four MBH scaling relations: MBH-sigma, MBH-L(sph), MBH-M(sph,*), MBH-M(sph,dyn). We find the following results. First, stellar-velocity dispersions determined from aperture spectra (e.g. SDSS fiber spectra or unresolved data from distant galaxies) can be biased, depending on aperture size, AGN contamination, and host-galaxy morphology. However, such a bias cannot explain the offset seen in the MBH-sigma relation at higher redshifts. Second, while the CaT region is the cleanest to determine stellar-velocity dispersions, both the MgIb region, corrected for FeII emission, and the CaHK region, although often swamped by the AGN powerlaw continuum and emission lines, can give results accurate to within a few percent. Third, the MBH scaling relations of our pilot sample agree in slope and scatter with those of other local active and inactive galaxies. In the next papers of the series we will quantify the scaling relations, exploiting the full sample of ~100 objects.

preprint2011arXiv

AMUSE-VIRGO. III: mid-infrared photometry of early-type galaxies and limits on obscured nuclear emission

We complete our census of low-level nuclear activity in Virgo Cluster early-type galaxies by searching for obscured emission using Spitzer Space Telescope mid-infrared (MIR) imaging at 24mu. Of a total sample of 95 early-type galaxies, 53 objects are detected, including 16 showing kiloparsec-scale dust in optical images. One dimensional and two dimensional surface photometry of the 37 detections without extended dust features reveals that the MIR light is more centrally concentrated than the optical light as traced by Hubble Space Telescope F850LP-band images. No such modeling was performed for the sources with dust detected in the optical images. We explore several possible sources of the MIR excess emission, including obscured nuclear emission. We find that radial metallicity gradients in the stellar population appear to be a natural and most likely explanation for the observed behavior in a majority of the sources. Alternatively, if the concentrated MIR emission were due to nuclear activity, it would imply a MIR-to-X luminosity ratio ~5-10 for the low luminosity AGN detected in X-rays by our survey. This ratio is an order of magnitude larger than that of typical low-luminosity AGN and would imply an unusual spectral energy distribution. We conclude that the black holes found by our survey in quiescent early-type galaxies in Virgo have low bolometric Eddington ratios arising from low accretion rates and/or highly radiatively inefficient accretion.

preprint2011arXiv

Geometric and Dynamical Models of Reverberation Mapping Data

We present a general method to analyze reverberation mapping data that provides both estimates for the black hole mass and for the geometry and dynamics of the broad line region (BLR) in active galactic nuclei (AGN). Our method directly infers the spatial and velocity distribution of the BLR from the data, allowing us to easily derive a velocity-resolved transfer function and allowing for a self-consistent estimate of the black hole mass without a virial coefficient. We obtain estimates and reasonable uncertainties of the BLR model parameters by implementing a Markov Chain Monte Carlo algorithm using the formalism of Bayesian probability theory. We use Gaussian Processes to interpolate the the continuum light curve data and create mock light curves that can be fitted to the data. We test our method by creating simulated reverberation mapping data-sets with known true parameter values and by trying to recover these parameter values using our models. We are able to recover the parameters with realistic uncertainties that depend upon the variability of the AGN and the quality of the reverberation mapping campaign. With a geometry model we can recover the mean radius of the BLR to within ~0.1dex random uncertainty for simulated data with an integrated line flux uncertainty of 1.5%, while with a dynamical model we can recover the black hole mass and the mean radius to within ~0.05dex random uncertainty, for simulated data with a line profile average signal to noise ratio of 4 per spectral pixel. These uncertainties do not include modeling errors, which are likely to be present in the analysis of real data, and should therefore be considered as lower limits to the accuracy of the method.

preprint2011arXiv

The Assembly History of Disk Galaxies: I - The Tully-Fisher Relation to z~1.3 from Deep Exposures with DEIMOS

We present new measures of the evolving scaling relations between stellar mass, luminosity and rotational velocity for a morphologically-inclusive sample of 129 disk-like galaxies with z_AB<22.5 in the redshift range 0.2<z<1.3, based on spectra from DEIMOS on the Keck II telescope, multi-color HST ACS photometry, and ground-based Ks-band imaging. A unique feature of our survey is the extended spectroscopic integration times, leading to significant improvements in determining characteristic rotational velocities for each galaxy. Rotation curves are reliably traced to the radius where they begin to flatten for ~90% of our sample, and we model the HST-resolved bulge and disk components in order to accurately de-project our measured velocities, accounting for seeing and dispersion. We demonstrate the merit of these advances by recovering an intrinsic scatter on the stellar mass Tully-Fisher relation a factor of 2-3 less than in previous studies at intermediate redshift and comparable to that of locally-determined relations. With our increased precision, we find the relation is well-established by <z>~1, with no significant evolution to <z>~0.3, \DeltaM_stellar ~ 0.04+/-0.07 dex. A clearer trend of evolution is seen in the B-band Tully-Fisher relation corresponding to a decline in luminosity of \DeltaM_B ~ 0.85+/-0.28 magnitudes at fixed velocity over the same redshift range, reflecting the changes in star formation over this period. As an illustration of the opportunities possible when gas masses are available for a sample such as ours, we show how our dynamical and stellar mass data can be used to evaluate the likely contributions of baryons and dark matter to the assembly history of spiral galaxies.

preprint2011arXiv

The Dark Matter Distribution in Abell 383: Evidence for a Shallow Density Cusp from Improved Lensing, Stellar Kinematic and X-ray Data

We extend our analyses of the dark matter (DM) distribution in relaxed clusters to the case of Abell 383, a luminous X-ray cluster at z=0.189 with a dominant central galaxy and numerous strongly-lensed features. Following our earlier papers, we combine strong and weak lensing constraints secured with Hubble Space Telescope and Subaru imaging with the radial profile of the stellar velocity dispersion of the central galaxy, essential for separating the baryonic mass distribution in the cluster core. Hydrostatic mass estimates from Chandra X-ray observations further constrain the solution. These combined datasets provide nearly continuous constraints extending from 2 kpc to 1.5 Mpc in radius, allowing stringent tests of results from recent numerical simulations. Two key improvements in our data and its analysis make this the most robust case yet for a shallow slope βof the DM density profile ρ_DM ~ r^-βon small scales. First, following deep Keck spectroscopy, we have secured the stellar velocity dispersion profile to a radius of 26 kpc for the first time in a lensing cluster. Secondly, we improve our previous analysis by adopting a triaxial DM distribution and axisymmetric dynamical models. We demonstate that in this remarkably well-constrained system, the logarithmic slope of the DM density at small radii is β< 1.0 (95% confidence). An improved treatment of baryonic physics is necessary, but possibly insufficient, to reconcile our observations with the recent results of high-resolution simulations.

preprint2011arXiv

The Lick AGN Monitoring Project: Recalibrating Single-Epoch Virial Black Hole Mass Estimates

We investigate the calibration and uncertainties of black hole mass estimates based on the single-epoch (SE) method, using homogeneous and high-quality multi-epoch spectra obtained by the Lick Active Galactic Nucleus (AGN) Monitoring Project for 9 local Seyfert 1 galaxies with black hole masses < 10^8 M_sun. By decomposing the spectra into their AGN and stellar components, we study the variability of the single-epoch Hbeta line width (full width at half-maximum intensity, FWHM_Hbeta; or dispersion, sigma_Hbeta) and of the AGN continuum luminosity at 5100A (L_5100). From the distribution of the "virial products" (~ FWHM_Hbeta^2 L_5100^0.5 or sigma_Hbeta^2 L_5100^0.5) measured from SE spectra, we estimate the uncertainty due to the combined variability as ~ 0.05 dex (12%). This is subdominant with respect to the total uncertainty in SE mass estimates, which is dominated by uncertainties in the size-luminosity relation and virial coefficient, and is estimated to be ~ 0.46 dex (factor of ~ 3). By comparing the Hbeta line profile of the SE, mean, and root-mean-square (rms) spectra, we find that the Hbeta line is broader in the mean (and SE) spectra than in the rms spectra by ~ 0.1 dex (25%) for our sample with FWHM_Hbeta < 3000 km/s. This result is at variance with larger mass black holes where the difference is typically found to be much less than 0.1 dex. To correct for this systematic difference of the Hbeta line profile, we introduce a line-width dependent virial factor, resulting in a recalibration of SE black hole mass estimators for low-mass AGNs.

preprint2011arXiv

The Mass of the Black Hole in Arp 151 from Bayesian Modeling of Reverberation Mapping Data

Supermassive black holes are believed to be ubiquitous at the centers of galaxies. Measuring their masses is extremely challenging yet essential for understanding their role in the formation and evolution of cosmic structure. We present a direct measurement of the mass of a black hole in an active galactic nucleus (Arp 151) based on the motion of the gas responsible for the broad emission lines. By analyzing and modeling spectroscopic and photometric time series, we find that the gas is well described by a disk or torus with an average radius of 3.99 +- 1.25 light days and an opening angle of 68.9 (+21.4, -17.2) degrees, viewed at an inclination angle of 67.8 +- 7.8 degrees (that is, closer to face-on than edge-on). The black hole mass is inferred to be 10^(6.51 +- 0.28) solar masses. The method is fully general and can be used to determine the masses of black holes at arbitrary distances, enabling studies of their evolution over cosmic time.

preprint2011arXiv

The Relation between Black Hole Mass and Host Spheroid Stellar Mass out to z~2

We combine Hubble Space Telescope images from the Great Observatories Origins Deep Survey with archival Very Large Telescope and Keck spectra of a sample of 11 X-ray selected broad-line active galactic nuclei in the redshift range 1<z<2 to study the black hole mass - stellar mass relation out to a lookback time of 10 Gyrs. Stellar masses of the spheroidal component are derived from multi-filter surface photometry. Black hole masses are estimated from the width of the broad MgII emission line and the 3000A nuclear luminosity. Comparing with a uniformly measured local sample and taking into account selection effects, we find evolution in the form M_BH/M_spheroid ~ (1+z)^(1.96+/-0.55), in agreement with our earlier studies based on spheroid luminosity. However, this result is more accurate because it does not require a correction for luminosity evolution and therefore avoids the related and dominant systematic uncertainty. We also measure total stellar masses. Combining our sample with data from the literature, we find M_BH/M_host ~ (1+z)^(1.15+/-0.15), consistent with the hypothesis that black holes (in the range M_BH ~ 10^8-9 M_sun) predate the formation of their host galaxies. Roughly one third of our objects reside in spiral galaxies; none of the host galaxies reveal signs of interaction or major merger activity. Combined with the slower evolution in host stellar masses compared to spheroid stellar masses, our results indicate that secular evolution or minor mergers play a non-negligible role in growing both BHs and spheroids.

preprint2011arXiv

The Sloan Lens ACS Survey. XI. Beyond Hubble resolution: size, luminosity and stellar mass of compact lensed galaxies at intermediate redshift

We exploit the strong lensing effect to explore the properties of intrinsically faint and compact galaxies at intermediate redshift, at the highest possible resolution at optical wavelengths. Our sample consists of 46 strongly-lensed emission line galaxies discovered by the Sloan Lens ACS (SLACS) Survey. The galaxies have been imaged at high resolution with HST in three bands (V_HST, I_814 and H_160), allowing us to infer their size, luminosity, and stellar mass using stellar population synthesis models. Lens modeling is performed using a new fast and robust code, klens, which we test extensively on real and synthetic non-lensed galaxies, and also on simulated galaxies multiply-imaged by SLACS- like galaxy-scale lenses. Our tests show that our measurements of galaxy size, flux, and Sersic index are robust and accurate, even for objects intrinsically smaller than the HST point spread function. The median magnification is 8.8, with a long tail that extends to magnifications above 40. Modeling the SLACS sources reveals a population of galaxies with colors and Sersic indices (median n ~ 1) consistent with the objects detected in the field with HST in the GEMS survey, but that are (typically) ~ 2 magnitudes fainter and ~ 5 times smaller in apparent size. The closest analog are ultracompact emission line galaxies identified by HST grism surveys. The lowest mass galaxies in our sample are comparable to the brightest Milky Way satellites in stellar mass (10^7 solar masses) and have well-determined half light radii of 0."05 (~0.3 kpc).

preprint2011arXiv

Two-dimensional kinematics of SLACS lenses - IV. The complete VLT-VIMOS data set

This paper presents the full VLT/VIMOS-IFU data set and related data products from an ESO Large Programme with the observational goal of obtaining two-dimensional kinematic data of early-type lens galaxies, out to one effective radius. The sample consists of 17 early-type galaxies (ETG) selected from the SLACS gravitational-lens survey. The galaxies cover the redshift range from 0.08 to 0.35 and have stellar velocity dispersions between 200 and 350 km/s. This programme is complemented by a similar observational programme on Keck, using long-slit spectroscopy. In combination with multi-band imaging data, the kinematic data provide stringent constraints on the inner mass profiles of ETGs beyond the local universe. Our Large Programme thus extends studies of nearby early-type galaxies (e.g. SAURON/ATLAS3D) by an order of magnitude in distance and toward higher masses. We provide an overview of our observational strategy, the data products (luminosity-weighted spectra and Hubble Space Telescope images) and derived products (i.e. two-dimensional fields of velocity dispersions and streaming motions) that have been used in a number of published and forthcoming lensing, kinematic and stellar-population studies.These studies also pave the way for future studies of early-type galaxies at z=1 with the upcoming extremely large telescopes.

preprint2011arXiv

Two-dimensional kinematics of SLACS lenses: III. Mass structure and dynamics of early-type lens galaxies beyond z ~ 0.1

We combine in a self-consistent way the constraints from both gravitational lensing and stellar kinematics to perform a detailed investigation of the internal mass distribution, amount of dark matter, and dynamical structure of the 16 early-type lens galaxies from the SLACS Survey, at z = 0.08 - 0.33, for which both HST/ACS and NICMOS high-resolution imaging and VLT VIMOS IFU spectroscopy are available. Based on this data set, we analyze the inner regions of the galaxies, i.e. typically within one (3D) effective radius r_e, under the assumption of axial symmetry and by constructing dynamical models supported by two-integral stellar DFs. For all systems, the total mass density distribution is found to be well approximated by a simple power-law: this profile is on average slightly super-isothermal, with a logarithmic slope <gamma'> = 2.074^{+0.043}_{-0.041} (68% CL) and an intrinsic scatter 0.144^{+0.055}_{-0.014}, and is fairly round, with an average axial ratio <q> = 0.77+/-0.04. The lower limit for the dark matter fraction (fDM) inside r_e ranges, in individual systems, from nearly zero to almost a half, with a median value of 12%. By including stellar masses derived from SPS models with a Salpeter IMF, we obtain an average fDM = 31%. The fDM rises to 61% if, instead, a Chabrier IMF is assumed. For both IMFs, the dark matter fraction increases with the total mass of the galaxy (3-sigma correlation). Based on the intrinsic angular momentum parameter calculated from our models, we find that the galaxies can be divided into two dynamically distinct groups, which are shown to correspond to the usual classes of the slow and fast rotators. Overall, the SLACS systems are structurally and dynamically very similar to their nearby counterparts, indicating that the inner regions of early-type galaxies have undergone little, if any, evolution since redshift z ~ 0.35. (Abridged)

preprint2010arXiv

AMUSE-Virgo II. Down-sizing in black hole accretion

(Abridged) We complete the census of nuclear X-ray activity in 100 early type Virgo galaxies observed by the Chandra X-ray Telescope as part of the AMUSE-Virgo survey, down to a (3sigma) limiting luminosity of 3.7E+38 erg/s over 0.5-7 keV. The stellar mass distribution of the targeted sample, which is mostly composed of formally `inactive' galaxies, peaks below 1E+10 M_Sun, a regime where the very existence of nuclear super-massive black holes (SMBHs) is debated. Out of 100 objects, 32 show a nuclear X-ray source, including 6 hybrid nuclei which also host a massive nuclear cluster as visible from archival HST images. After carefully accounting for contamination from nuclear low-mass X-ray binaries based on the shape and normalization of their X-ray luminosity function, we conclude that between 24-34% of the galaxies in our sample host a X-ray active SMBH (at the 95% C.L.). This sets a firm lower limit to the black hole occupation fraction in nearby bulges within a cluster environment. At face value, the active fraction -down to our luminosity limit- is found to increase with host stellar mass. However, taking into account selection effects, we find that the average Eddington-scaled X-ray luminosity scales with black hole mass as M_BH^(-0.62^{+0.13}_{-0.12}), with an intrinsic scatter of 0.46^({+0.08}_{-0.06}) dex. This finding can be interpreted as observational evidence for `down-sizing' of black hole accretion in local early types, that is, low mass black holes shine relatively closer to their Eddington limit than higher mass objects. As a consequence, the fraction of active galaxies, defined as those above a fixed X-ray Eddington ratio, decreases with increasing black hole mass.

preprint2010arXiv

Cosmic Evolution of Virial and Stellar Mass in Massive Early-Type Galaxies

We measure the average mass properties of a sample of 41 strong gravitational lenses at moderate redshift (z ~ 0.4 - 0.9), and present the lens redshift for 6 of these galaxies for the first time. Using the techniques of strong and weak gravitational lensing on archival data obtained from the Hubble Space Telescope, we determine that the average mass overdensity profile of the lenses can be fit with a power-law profile (Delta_Sigma prop. to R^{-0.86 +/- 0.16}) that is within 1-sigma of an isothermal profile (Delta_Sigma prop. to R^{-1}) with velocity dispersion sigma_v = 260 +/- 20 km/s. Additionally, we use a two-component de Vaucouleurs+NFW model to disentangle the total mass profile into separate luminous and dark matter components, and determine the relative fraction of each component. We measure the average rest frame V-band stellar mass-to-light ratio (Upsilon_V = 4.0 +/- 0.6 h M_sol/L_sol) and virial mass-to-light ratio (tau_V = 300 +/- 90 h M_sol/L_sol) for our sample, resulting in a virial-to-stellar mass ratio of M_vir/M_* = 75 +/- 25. Finally, we compare our results to a previous study using low redshift lenses, to understand how galaxy mass profiles evolve over time. We investigate the evolution of M_vir/M_*(z) = alpha(1+z)^{beta}, and find best fit parameters of alpha = 51 +/- 36 and beta = 0.9 +/- 1.8, constraining the growth of virial to stellar mass ratio over the last ~7 Gigayears. We note that, by using a sample of strong lenses, we are able to constrain the growth of M_vir/M_*(z) without making any assumptions about the IMF of the stellar population.

preprint2010arXiv

Keck Spectroscopy of z>1 Field Spheroidals: Dynamical Constraints on the Growth Rate of Red "Nuggets"

We present deep Keck spectroscopy for 17 morphologically-selected field spheroidals in the redshift range 1.05<z<1.60 in order to investigate the continuity in physical properties between the claimed massive compact red galaxies ("nuggets") at z~2 and well-established data for massive spheroidal galaxies below z~1. By combining Keck-based stellar velocity dispersions with HST-based sizes, we find that the most massive systems (Mdyn > 10^11 Msol) grew in size over 0<z<1.6 as (1+z)^(-0.75 +- 0.10) (i.e., x2 since z=1.5) whereas intermediate mass systems (10^11 Msol > Mdyn > 10^10 Msol) did not grow significantly. These trends are consistent with a picture in which more massive spheroidals formed at higher redshift via "wetter" mergers involving greater dissipation. To examine growth under the favored "dry" merger hypothesis, we also examine size growth at a fixed velocity dispersion. This test, uniquely possible with our dynamical data, allows us to consider the effects of "progenitor bias." Above our completeness limit (sigma > 200 km/s), we find size growth consistent with that inferred for the mass-selected sample, thus ruling out strong progenitor bias. To maintain continuity in the growth of massive galaxies over the past 10 Gyr, our new results imply that size evolution over 1.3<z<2.3, a period of 1.9 Gyr, must have been even more dramatic than hitherto claimed if the red sources at z>2 are truly massive and compact.

preprint2010arXiv

The Lick AGN Monitoring Project: Alternate Routes to a Broad-line Region Radius

It is now possible to estimate black hole masses across cosmic time, using broad emission lines in active galaxies. This technique informs our views of how galaxies and their central black holes coevolve. Unfortunately, there are many outstanding uncertainties associated with these "virial" mass estimates. One of these comes from using the accretion luminosity to infer a size for the broad-line region. Incorporating the new sample of low-luminosity active galaxies from our recent monitoring campaign at Lick Observatory, we recalibrate the radius-luminosity relation with tracers of the accretion luminosity other than the optical continuum. We find that the radius of the broad-line region scales as the square root of the X-ray and Hbeta luminosities, in agreement with recent optical studies. On the other hand, the scaling appears to be marginally steeper with narrow-line luminosities. This is consistent with a previously observed decrease in the ratio of narrow-line to X-ray luminosity with increasing total luminosity. The radius of the broad-line region correlates most tightly with Hbeta luminosity, while the X-ray and narrow-line relations both have comparable scatter of a factor of two. These correlations provide useful alternative virial BH masses in objects with no detectable optical/UV continuum emission, such as high-redshift galaxies with broad emission lines, radio-loud objects, or local active galaxies with galaxy-dominated continua.

preprint2010arXiv

The Lick AGN Monitoring Project: Reverberation Mapping of Optical Hydrogen and Helium Recombination Lines

We have recently completed a 64-night spectroscopic monitoring campaign at the Lick Observatory 3-m Shane telescope with the aim of measuring the masses of the black holes in 12 nearby (z < 0.05) Seyfert 1 galaxies with expected masses in the range ~10^6-10^7M_sun and also the well-studied nearby active galactic nucleus (AGN) NGC 5548. Nine of the objects in the sample (including NGC 5548) showed optical variability of sufficient strength during the monitoring campaign to allow for a time lag to be measured between the continuum fluctuations and the response to these fluctuations in the broad Hbeta emission, which we have previously reported. We present here the light curves for the Halpha, Hgamma, HeII 4686, and HeI 5876 emission lines and the time lags for the emission-line responses relative to changes in the continuum flux. Combining each emission-line time lag with the measured width of the line in the variable part of the spectrum, we determine a virial mass of the central supermassive black hole from several independent emission lines. We find that the masses are generally consistent within the uncertainties. The time-lag response as a function of velocity across the Balmer line profiles is examined for six of the AGNs. Finally we compare several trends seen in the dataset against the predictions from photoionization calculations as presented by Korista & Goad. We confirm several of their predictions, including an increase in responsivity and a decrease in the mean time lag as the excitation and ionization level for the species increases. Further confirmation of photoionization predictions for broad-line gas behavior will require additional monitoring programs for these AGNs while they are in different luminosity states. [abridged]

preprint2010arXiv

The Lick AGN Monitoring Project: The MBH - sigma Relation For Reverberation-Mapped Active Galaxies

(Abridged) To investigate the black hole mass (MBH) vs. stellar velocity dispersion relation of active galaxies, we measured the velocity dispersions of a sample of local Seyfert 1 galaxies, for which we have recently determined MBH using reverberation mapping. For most objects, velocity dispersions were measured from high S/N ratio optical spectra centered on the Ca II triplet region (~8500 A), obtained at the Keck, Palomar, and Lick Observatories. For two objects, in which the Ca II triplet region was contaminated by nuclear emission, the measurement was based on high-quality H-band spectra obtained with the OSIRIS at the Keck-II Telescope. Combining our new measurements with data from the literature, we assemble a sample of 24 active galaxies with stellar velocity dispersions and reverberation-based MBH measurements in the range of black hole mass 10^6<MBH/M_sun<10^9. We use this sample to obtain reverberation mapping constraints on the slope and intrinsic scatter of the M-sigma relation of active galaxies. Assuming a constant virial coefficient f for the reverberation mapping black hole masses, we find a slope b=3.55+-0.60 and the intrinsic scatter sigma_int=0.43+-0.08 dex in the relation log (MBH/M_sun) = a + b log (sigma/200 km/s), which are consistent with those found for quiescent galaxies. We derive an updated value of the virial coefficient f by finding the value which places the reverberation masses in best agreement with the M-sigma relation of quiescent galaxies; using the quiescent M-sigma relation determined by Gultekin et al. we find log f=0.72^{+0.09-0.10} with an intrinsic scatter of 0.44+-0.07 dex. The uncertainty of the virial coefficient remains one of the main sources of the uncertainty in black hole mass determinations using reverberation mapping, and therefore also in single-epoch spectroscopic estimates of MBH in active galaxies.

preprint2010arXiv

The non-evolving internal structure of early-type galaxies: the case study SDSS J0728+3835 at z = 0.206

We study the internal dynamical structure of the early-type lens galaxy SDSS J0728+3835 at z = 0.206. The analysis is based on two-dimensional kinematic maps extending out to 1.7 effective radii obtained from Keck spectroscopy, on lensing geometry and on stellar mass estimates obtained from multiband Hubble Space Telescope imaging. The data are modelled under the assumptions of axial symmetry supported by a two-integral distribution function (DF), by applying the combined gravitational lensing and stellar dynamics code CAULDRON, and yielding high-quality constraints for an early-type galaxy at cosmological redshifts. Modelling the total density profile as a power-law of the form rho_tot ~ 1/r^{gamma}, we find that it is nearly isothermal (logarithmic slope gamma = 2.08^{+0.04}_{-0.02}), and quite flattened (axial ratio q = 0.60^{+0.08}_{-0.03}). The galaxy is mildly anisotropic (delta = 0.08 +/- 0.02) and shows a fair amount of rotational support, in particular towards the outer regions. We determine a dark matter fraction lower limit of 28 per cent within the effective radius. The stellar contribution to the total mass distribution is close to maximal for a Chabrier initial mass function (IMF), whereas for a Salpeter IMF the stellar mass exceeds the total mass within the galaxy inner regions. We find that the combination of a NFW dark matter halo with the maximally rescaled luminous profile provides a remarkably good fit to the total mass distribution over a broad radial range. Our results confirm and expand the findings of the SLACS survey for early-type galaxies of comparable velocity dispersion (sigma_SDSS = 214 +/- 11 km/s). The internal structure of J0728 is consistent with that of local early-type galaxies of comparable velocity dispersion as measured by the SAURON project, suggesting lack of evolution in the past two billion years.

preprint2010arXiv

The SL2S Galaxy-scale Lens Sample. II. Cosmic evolution of dark and luminous mass in early-type galaxies

We present a joint gravitational lensing and stellar-dynamical analysis of 11 early-type galaxies (median deflector redshift $\zd=0.5$) from Strong Lenses in the Legacy Survey (SL2S). Using newly measured redshifts and stellar velocity dispersions from Keck spectroscopy with lens models from Paper I, we derive the total mass density slope inside the Einstein radius for each of the 11 lenses. The average total density slope is found to be $<γ'>= 2.16^{+0.09}_{-0.09}$ ($ρ_{\rm tot}\propto r^{-γ'}$), with an intrinsic scatter of $0.25^{+0.10}_{-0.07}$. We also determine the dark matter fraction for each lens within half the effective radius, and find the average projected dark matter mass fraction to be $0.42^{+0.08}_{-0.08}$ with a scatter of $0.20^{+0.09}_{-0.07}$ for a Salpeter IMF. By combining the SL2S results with those from the Sloan Lens ACS Survey (median $\zd=0.2$) and the Lenses Structure and Dynamics survey (median $\zd=0.8$), we investigate cosmic evolution of $γ'$ and find a mild trend $\partial<γ'>/\partial\zd = -0.25^{+0.10}_{-0.12}$. This suggests that the total density profile of massive galaxies has become slightly steeper over cosmic time. If this result is confirmed by larger samples, it would indicate that dissipative processes played some role in the growth of massive galaxies since $z\sim1$.

preprint2009arXiv

Cosmic Evolution of Black Holes and Spheroids. IV. The BH Mass - Spheroid Luminosity Relation

From high-resolution images of 23 Seyfert-1 galaxies at z=0.36 and z=0.57 obtained with the Near Infrared Camera and Multi-Object Spectrometer on board the Hubble Space Telescope (HST), we determine host-galaxy morphology, nuclear luminosity, total host-galaxy luminosity and spheroid luminosity. Keck spectroscopy is used to estimate black hole mass (M_BH). We study the cosmic evolution of the M_BH-spheroid luminosity (L_sph) relation. In combination with our previous work, totaling 40 Seyfert-1 galaxies, the covered range in BH mass is substantially increased, allowing us to determine for the first time intrinsic scatter and correct evolutionary trends for selection effects. We re-analyze archival HST images of 19 local reverberation-mapped active galaxies to match the procedure adopted at intermediate redshift. Correcting spheroid luminosity for passive luminosity evolution and taking into account selection effects, we determine that at fixed present-day V-band spheroid luminosity, M_BH/L_sph \propto (1+z)^(2.8+/-1.2). When including a sample of 44 quasars out to z=4.5 taken from the literature, with luminosity and BH mass corrected to a self-consistent calibration, we extend the BH mass range to over two orders of magnitude, resulting in M_BH/L_sph \propto (1+z)^(1.4+/-0.2). The intrinsic scatter of the relation, assumed constant with redshift, is 0.3+/-0.1 dex (<0.6 dex at 95% CL). The evolutionary trend suggests that BH growth precedes spheroid assembly. Interestingly, the M_BH-total host-galaxy luminosity relation is apparently non-evolving. It hints at either a more fundamental relation or that the spheroid grows by a redistribution of stars. However, the high-z sample does not follow this relation, indicating that major mergers may play the dominant role in growing spheroids above z~1.

preprint2009arXiv

Direct Observation of Cosmic Strings via their Strong Gravitational Lensing Effect: II. Results from the HST/ACS Image Archive

We have searched 4.5 square degrees of archival HST/ACS images for cosmic strings, identifying close pairs of similar, faint galaxies and selecting groups whose alignment is consistent with gravitational lensing by a long, straight string. We find no evidence for cosmic strings in five large-area HST treasury surveys (covering a total of 2.22 square degrees), or in any of 346 multi-filter guest observer images (1.18 square degrees). Assuming that simulations ccurately predict the number of cosmic strings in the universe, this non-detection allows us to place upper limits on the unitless Universal cosmic string tension of G mu/c^2 < 2.3 x 10^-6, and cosmic string density of Omega_s < 2.1 x 10^-5 at the 95% confidence level (marginalising over the other parameter in each case). We find four dubious cosmic string candidates in 318 single filter guest observer images (1.08 square degrees), which we are unable to conclusively eliminate with existing data. The confirmation of any one of these candidates as cosmic strings would imply G mu/c^2 ~ 10^-6 and Omega_s ~ 10^-5. However, we estimate that there is at least a 92% chance that these string candidates are random alignments of galaxies. If we assume that these candidates are indeed false detections, our final limits on G mu/c^2 and Omega_s fall to 6.5 x 10^-7 and 7.3 x 10^-6. Due to the extensive sky coverage of the HST/ACS image archive, the above limits are universal. They are quite sensitive to the number of fields being searched, and could be further reduced by more than a factor of two using forthcoming HST data.

preprint2009arXiv

Dry mergers and the formation of early-type galaxies: constraints from lensing and dynamics

Dissipationless (gas-free or "dry") mergers have been suggested to play a major role in the formation and evolution of early-type galaxies, particularly in growing their mass and size without altering their stellar populations. We perform a new test of the dry merger hypothesis by comparing N-body simulations of realistic systems to empirical constraints provided by recent studies of lens early-type galaxies. We find that major and minor dry mergers: i) preserve the nearly isothermal structure of early-type galaxies within the observed scatter; ii) do not change more than the observed scatter the ratio between total mass M and "virial" mass R_e*sigma/2G (where R_e is the half-light radius and sigma the projected velocity dispersion); iii) increase strongly galaxy sizes [as M^(0.85+/-0.17)] and weakly velocity dispersions [as M^(0.06+/-0.08)] with mass, thus moving galaxies away from the local observed M-R_e and M-sigma relations; iv) introduce substantial scatter in the M-R_e and M-sigma relations. Our findings imply that, unless there is a high degree of fine tuning of the mix of progenitors and types of interactions, present-day massive early-type galaxies cannot have assembled more than ~50% of their mass, and increased their size by more than a factor ~1.8, via dry merging.

preprint2009arXiv

Environmental Effects in the Evolution of Galactic Bulges

We investigate possible environmental trends in the evolution of galactic bulges over the redshift range 0<z<0.6. For this purpose, we construct the Fundamental Plane (FP) for cluster and field samples at redshifts <z>=0.4 and <z>=0.54 using surface photometry based on HST imaging and velocity dispersions based on Keck spectroscopy. As a reference point for our study we include data for pure ellipticals, which we model as single-component Sersic profiles; whereas for multi-component galaxies we undertake decompositions using Sersic and exponential models for the bulge and disk respectively. Although the FP for both distant cluster and field samples are offset from the local relation, consistent with evolutionary trends found in earlier studies, we detect significant differences in the zero point of ~=0.2 dex between the field and cluster samples at a given redshift. For both clusters, the environmentally-dependent offset is in the sense expected for an accelerated evolution of bulges in dense environments. By matching the mass range of our samples, we confirm that this difference does not arise as a result of the mass-dependent downsizing effects seen in larger field samples. Our result is also consistent with the hypothesis that - at fixed mass and environment - the star formation histories of galactic bulges and pure spheroids are indistinguishable, and difficult to reconcile with the picture whereby the majority of large bulges form primarily via secular processes within spiral galaxies.

preprint2009arXiv

Lyman Break Galaxies at z~1.8-2.8: GALEX/NUV Imaging of the Subaru Deep Field

Abridged: A photometric sample of ~7100 V<25.3 Lyman break galaxies (LBGs) has been selected by combining Subaru/Suprime-Cam BVRci'z' data with deep GALEX/NUV imaging of the Subaru Deep Field. Follow-up spectroscopy confirmed 24 LBGs at 1.5<z<2.7. Among the optical spectra, 12 have Ly-alpha emission with rest-frame equivalent widths of ~5-60AA. The success rate for identifying LBGs as NUV-dropouts at 1.5<z<2.7 is 86%. The rest-frame UV (1700AA) luminosity function (LF) is constructed from the photometric sample with corrections for stellar contamination and z<1.5 interlopers. The LF is 1.7+/-0.1 times higher than those of z~2 BXs and z~3 LBGs. Three explanations were considered, and it is argued that significantly underestimating low-z contamination or effective comoving volume is unlikely: the former would be inconsistent with the spectroscopic sample at 93% confidence, and the second explanation would not resolve the discrepancy. The third scenario is that different photometric selection of the samples yields non-identical galaxy populations, such that some BX galaxies are LBGs and vice versa. This argument is supported by a higher surface density of LBGs at all magnitudes while the redshift distribution of the two populations is nearly identical. This study, when combined with other star-formation rate (SFR) density UV measurements from LBG surveys, indicates that there is a rise in the SFR density: a factor of 3-6 (3-10) increase from z~5 (z~6) to z~2, followed by a decrease to z~0. This result, along with past sub-mm studies that find a peak at z~2 in their redshift distribution, suggest that z~2 is the epoch of peak star-formation. Additional spectroscopy is required to characterize the complete shape of the z~2 LBG UV LF via measurements of contamination and accurate distances.

preprint2009arXiv

The Lick AGN Monitoring Project: Broad-Line Region Radii and Black Hole Masses from Reverberation Mapping of Hbeta

We have recently completed a 64-night spectroscopic monitoring campaign at the Lick Observatory 3-m Shane telescope with the aim of measuring the masses of the black holes in 12 nearby (z < 0.05) Seyfert 1 galaxies with expected masses in the range ~10^6-10^7 M_sun and also the well-studied nearby active galactic nucleus (AGN) NGC 5548. Nine of the objects in the sample (including NGC 5548) showed optical variability of sufficient strength during the monitoring campaign to allow for a time lag to be measured between the continuum fluctuations and the response to these fluctuations in the broad Hbeta emission. We present here the light curves for the objects in this sample and the subsequent Hbeta time lags for the nine objects where these measurements were possible. The Hbeta lag time is directly related to the size of the broad-line region, and by combining the lag time with the measured width of the Hbeta emission line in the variable part of the spectrum, we determine the virial mass of the central supermassive black hole in these nine AGNs. The absolute calibration of the black hole masses is based on the normalization derived by Onken et al. We also examine the time lag response as a function of velocity across the Hbeta line profile for six of the AGNs. The analysis of four leads to ambiguous results with relatively flat time lags as a function of velocity. However, SBS 1116+583A exhibits a symmetric time lag response around the line center reminiscent of simple models for circularly orbiting broad-line region (BLR) clouds, and Arp 151 shows an asymmetric profile that is most easily explained by a simple gravitational infall model. Further investigation will be necessary to fully understand the constraints placed on physical models of the BLR by the velocity-resolved response in these objects.

preprint2009arXiv

The Lick AGN Monitoring Project: Photometric Light Curves and Optical Variability Characteristics

The Lick AGN Monitoring Project targeted 13 nearby Seyfert 1 galaxies with the intent of measuring the masses of their central black holes using reverberation mapping. The sample includes 12 galaxies selected to have black holes with masses roughly in the range 10^6-10^7 solar masses, as well as the well-studied AGN NGC 5548. In conjunction with a spectroscopic monitoring campaign, we obtained broad-band B and V images on most nights from 2008 February through 2008 May. The imaging observations were carried out by four telescopes: the 0.76-m Katzman Automatic Imaging Telescope (KAIT), the 2-m Multicolor Active Galactic Nuclei Monitoring (MAGNUM) telescope, the Palomar 60-in (1.5-m) telescope, and the 0.80-m Tenagra II telescope. Having well-sampled light curves over the course of a few months is useful for obtaining the broad-line reverberation lag and black hole mass, and also allows us to examine the characteristics of the continuum variability. In this paper, we discuss the observational methods and the photometric measurements, and present the AGN continuum light curves. We measure various variability characteristics of each of the light curves. We do not detect any evidence for a time lag between the B- and V-band variations, and we do not find significant color variations for the AGNs in our sample.

preprint2009arXiv

Two-dimensional kinematics of SLACS lenses: II. Combined lensing and dynamics analysis of early-type galaxies at z = 0.08 - 0.33

We present the first detailed analysis of the mass and dynamical structure of a sample of six early-type lens galaxies, selected from the Sloan Lens ACS Survey, in the redshift range 0.08 < z < 0.33. Both Hubble Space Telescope (HST)/ACS high-resolution imaging and VLT VIMOS integral field spectroscopy are available for these systems. The galaxies are modelled - under the assumptions of axial symmetry and two-integral stellar distribution function - by making use of the CAULDRON code, which self-consistently combines gravitational lensing and stellar dynamics, and is fully embedded within the framework of Bayesian statistics. The principal results of this study are: (i) all galaxies in the sample are well described by a simple axisymmetric power-law profile for the total density, with a logarithmic slope gamma very close to isothermal (<gamma> = 1.98 +/- 0.05 and an intrinsic spread close to 5 per cent) showing no evidence of evolution over the probed range of redshift; (ii) the axial ratio of the total density distribution is rounder than 0.65 and in all cases, except for a fast rotator, does not deviate significantly from the flattening of the intrinsic stellar distribution; (iii) the dark matter fraction within the effective radius has a lower limit of about 15 to 30 per cent; (iv) the sample galaxies are only mildly anisotropic, with delta <= 0.16; (v) the physical distinction among slow and fast rotators, quantified by the v/sigma ratio and the intrinsic angular momentum, is already present at z > 0.1. Altogether, early-type galaxies at z = 0.08 - 0.33 are found to be markedly smooth and almost isothermal systems, structurally and dynamically very similar to their nearby counterparts. (Abridged)

preprint2008arXiv

Enhanced lensing rate by clustering of massive galaxies: newly discovered systems in the SLACS fields

[Abridged] We exploit the clustering of massive galaxies to perform a high efficiency imaging search for gravitational lenses. Our dataset comprises 44 fields imaged by the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS), each of which is centered on a lens discovered by the Strong Lens ACS Survey (SLACS). We compare four different search methods: 1) automated detection with the HST Archive Galaxy-scale Gravitational Lens Survey (HAGGLeS) robot, 2) examining cutout images of bright galaxies (BGs) after subtraction of a smooth galaxy light distribution, 3) examining the unsubtracted BG cutouts, and 4) performing a full-frame visual inspection of the ACS images. We compute purity and completeness and consider investigator time for the four algorithms, using the main SLACS lenses as a testbed. The first and second algorithms perform the best. We present the four new lens systems discovered during this comprehensive search, as well as one other likely candidate. For each new lens we use the fundamental plane to estimate the lens velocity dispersion and predict, from the resulting lens geometry, the redshifts of the lensed sources. Two of these new systems are found in galaxy clusters, which include the SLACS lenses in the two respective fields. Overall we find that the enhanced lens abundance (30^{+24}_{-8} lenses/degree^2) is higher than expected for random fields (12^{+4}_{-2} lenses/degree^2 for the COSMOS survey). Additionally, we find that the gravitational lenses we detect are qualitatively different from those in the parent SLACS sample: this imaging survey is largely probing higher-redshift, and lower-mass, early-type galaxies.

preprint2008arXiv

The nature of dusty starburst galaxies in a rich cluster at z=0.4: the progenitors of lenticulars?

We present the results of a Spitzer Infrared Spectrograph (IRS) survey of 24um-selected luminous infrared galaxies (LIRGs, L_IR > 10^11 L_sun) in the rich cluster Cl0024+16 at z=0.4. Optically, these LIRGs resemble unremarkable spiral galaxies with e(a)/e(c) spectral classifications and [Oii]-derived star formation rates (SFRs) of <2 M_sun/yr, generally indistinguishable from the 'quiescent' star forming population in the cluster. Our IRS spectra show that the majority of the 24um-detected galaxies exhibit polycyclic aromatic hydrocarbon (PAH) emission with implied SFRs ~30-60 M_sun/yr, with only one (<10%) of the sample displaying unambiguous evidence of an active galactic nucleus in the mid-infrared. This confirms the presence of a large population of obscured starburst galaxies in distant clusters, which comprise the bulk of the star formation occurring in these environments at z~0.5. We suggest that, although several mechanisms could be at play, these dusty starbursts could be the signature of an important evolutionary transition converting gas-rich spiral galaxies in distant clusters into the passive, bulge-dominated lenticular galaxies that become increasingly abundant in the cores of rich clusters in the ~4Gyr to the present day.

preprint2008arXiv

The Sloan Lens ACS Survey. VII. Elliptical Galaxy Scaling Laws from Direct Observational Mass Measurements

We use a sample of 53 massive early-type strong gravitational lens galaxies with well-measured redshifts (ranging from z=0.06 to 0.36) and stellar velocity dispersions (between 175 and 400 km/s) from the Sloan Lens ACS (SLACS) Survey to derive numerous empirical scaling relations. The ratio between central stellar velocity dispersion and isothermal lens-model velocity dispersion is nearly unity within errors. The SLACS lenses define a fundamental plane (FP) that is consistent with the FP of the general population of early-type galaxies. We measure the relationship between strong-lensing mass M_lens within one-half effective radius (R_e/2) and the dimensional mass variable M_dim = G^-1 sigma_e2^2 R_e/2 to be log_10 [M_lens/10^11 M_Sun] = (1.03 +/- 0.04) log_10 [M_dim/10^11 M_Sun] + (0.54 +/- 0.02) (where sigma_e2 is the projected stellar velocity dispersion within R_e/2). The near-unity slope indicates that the mass-dynamical structure of massive elliptical galaxies is independent of mass, and that the "tilt" of the SLACS FP is due entirely to variation in total (luminous plus dark) mass-to-light ratio with mass. Our results imply that dynamical masses serve as a good proxies for true masses in massive elliptical galaxies. Regarding the SLACS lenses as a homologous population, we find that the average enclosed 2D mass profile goes as log_10 [M(<R)/M_dim] = (1.10 +/- 0.09) log_10 [R/R_e] + (0.85 +/- 0.03), consistent with an isothermal (flat rotation curve) model when de-projected into 3D. This measurement is inconsistent with the slope of the average projected aperture luminosity profile at a confidence level greater than 99.9%, implying a minimum dark-matter fraction of f_DM = 0.38 +/- 0.07 within one effective radius. (abridged)

preprint2007arXiv

Comparing and calibrating black hole mass estimators for distant active galactic nuclei

Black hole mass is a fundamental property of active galactic nuclei (AGNs). In the distant universe, black hole mass is commonly estimated using the MgII, Hbeta, or Halpha emission line widths and the optical/UV continuum or line luminosities, as proxies for the characteristic velocity and size of the broad-line region. Although they all have a common calibration in the local universe, a number of different recipes are currently used in the literature. It is important to verify the relative accuracy and consistency of the recipes, as systematic changes could mimic evolutionary trends when comparing various samples. At z=0.36, all three lines can be observed at optical wavelengths, providing a unique opportunity to compare different empirical recipes. We use spectra from the Keck Telescope and the Sloan Digital Sky Survey to compare black hole mass estimators for a sample of nineteen AGNs at this redshift. We compare popular recipes available from the literature, finding that mass estimates can differ up to 0.38+-0.05 dex in the mean (or 0.13+-0.05 dex, if the same virial coefficient is adopted). Finally, we provide a set of 30 internally self consistent recipes for determining black hole mass from a variety of observables. The intrinsic scatter between cross-calibrated recipes is in the range 0.1-0.3 dex. This should be considered as a lower limit to the uncertainty of the black hole mass estimators.

preprint2007arXiv

Dynamical Evidence for Environmental Evolution of Intermediate Redshift Spiral Galaxies

Combining resolved optical spectroscopy with panoramic HST imaging, we study the dynamical properties of spiral galaxies as a function of position across two intermediate redshift clusters, and we compare the cluster population to field galaxies in the same redshift range. By modelling the observed rotation curves, we derive maximal rotation velocities for 40 cluster spirals and 37 field spirals, yielding one of the largest matched samples of cluster and field spirals at intermediate redshift. We construct the Tully-Fisher relation in both V and K_S bands, and find that the cluster Tully-Fisher relation exhibits significantly higher scatter than the field relation, in both V and K_S bands. Under the assumption that this increased scatter is due to an interaction with the cluster environment, we examine several dynamical quantities (dynamical mass, mass-to-light ratio, and central mass density) as a function of cluster environment. We find that the central mass densities of star-forming spirals exhibit a sharp break near the cluster Virial radius, with spirals in the cluster outskirts exhibiting significantly lower densities. We argue that the lower-density spirals in the cluster outskirts, combined with the high scatter in both K_S- and V-band TF relations, demonstrate that cluster spirals are kinematically disturbed by their environment, even as far as 2R_VIR from the cluster center. We propose that such disturbances may be due to a combination of galaxy merging and harassment.

preprint2007arXiv

The Dynamical Distinction between Elliptical and Lenticular Galaxies in Distant Clusters: Further Evidence for the Recent Origin of S0 Galaxies

We examine resolved spectroscopic data obtained with the Keck II telescope for 44 spheroidal galaxies in the fields of two rich clusters, Cl0024+16 (z=0.40) and MS0451-03 (z=0.54), and contrast this with similar data for 23 galaxies within the redshift interval 0.3<z<0.65 in the GOODS northern field. For each galaxy we examine the case for systemic rotation, derive central stellar velocity dispersions sigma and photometric ellipticities, epsilon. Using morphological classifications obtained via Hubble Space Telescope imaging as the basis, we explore the utility of our kinematic quantities in distinguishing between pressure-supported ellipticals and rotationally-supported lenticulars (S0s). We demonstrate the reliability of using the v/(1-epsilon) vs sigma and v/sigma vs epsilon distributions as discriminators, finding that the two criteria correctly identify 63%+-3% and 80%+-2% of S0s at z~0.5, respectively, along with 76%+8-3% and 79%+-2% of ellipticals. We test these diagnostics using equivalent local data in the Coma cluster, and find that the diagnostics are similarly accurate at z=0. Our measured accuracies are comparable to the accuracy of visual classification of morphologies, but avoid the band-shifting and surface brightness effects that hinder visual classification at high redshifts. As an example application of our kinematic discriminators, we then examine the morphology-density relation for elliptical and S0 galaxies separately at z~0.5. We confirm, from kinematic data alone, the recent growth of rotationally-supported spheroidals. We discuss the feasibility of extending the method to a more comprehensive study of cluster and field galaxies to z~1, in order to verify in detail the recent density-dependent growth of S0 galaxies.

preprint2007arXiv

The Sloan Lens ACS Survey. IV: the mass density profile of early-type galaxies out to 100 effective radii

We present a weak gravitational lensing analysis of 22 early-type strong lens galaxies, based on deep HST images obtained as part of the Sloan Lens ACS Survey. Using the most advanced techniques to control systematic uncertainties related to the variable PSF and charge transfer efficiency of the ACS, we detect weak lensing signal out to 300 kpc/h. We analyze blank control fields from the COSMOS survey in the same manner, inferring that the residual systematic uncertainty in the tangential shear is <0.3%. A joint strong and weak lensing analysis shows that the average total mass density profile is consistent with isothermal over two decades in radius (3-300 kpc/h, approximately 1-100 Reff). This finding extends by over an order of magnitude in radius previous results, based on strong lensing and/or stellar dynamics, that luminous and dark component ``conspire'' to form an isothermal mass distribution. In order to disentangle the contributions of luminous and dark matter, we fit a two-component mass model (R^1/4 + NFW) to the weak and strong lensing constraints. It provides a good fit to the data with only two free parameters; i) the average stellar mass-to-light ratio M_*/L_V=4.48 +- 0.46 hMo/Lo, in agreement with that expected for an old stellar population; ii) the average virial mass-to-light ratio M_vir/L_V = 246+101-87 hMo/Lo. [abridged]

preprint2005arXiv

The Sloan Lens ACS Survey. I. A Large Spectroscopically Selected Sample of Massive Early-Type Lens Galaxies

The Sloan Lens ACS (SLACS) Survey is an efficient Hubble Space Telescope Snapshot imaging survey for new galaxy-scale strong gravitational lenses. The targeted lens candidates are selected spectroscopically from within the Sloan Digital Sky Survey (SDSS) database of galaxy spectra for having multiple nebular emission lines at a redshift significantly higher than that of the SDSS target galaxy. In this paper, we present a catalog of 19 newly discovered gravitational lenses, along with 9 other observed candidate systems that are either possible lenses, non-lenses, or non-detections. The survey efficiency is thus >=68%. We also present Gemini and Magellan IFU data for 9 of the SLACS targets, which further support the lensing interpretation. A new method for the effective subtraction of foreground galaxy images to reveal faint background features is presented. We show that the SLACS lens galaxies have colors and ellipticities typical of the spectroscopic parent sample from which they are drawn (SDSS luminous red galaxies and quiescent main-sample galaxies), but are somewhat brighter and more centrally concentrated. Several explanations for the latter bias are suggested. The SLACS survey provides the first statistically significant and homogeneously selected sample of bright early-type lens galaxies, furnishing a powerful probe of the structure of early-type galaxies within the half-light radius. The high confirmation rate of lenses in the SLACS survey suggests consideration of spectroscopic lens discovery as an explicit science goal of future spectroscopic galaxy surveys (abridged).

preprint2004arXiv

Evolution Since z=1 of the Morphology-density Relation of Galaxies

We measure the morphology-density relation of galaxies at z=1 across the full three orders of magnitude in projected galaxy density available in low-redshift studies. Our study adopts techniques that are comparable with those applied at lower redshifts, allowing a direct investigation of how the morphological segregation of galaxies has evolved over the last 8 Gyr. Although the morphology-density relation, as described by the fraction of early-type (E+S0) galaxies, was in place at z=1, its form differs from that observed at both z=0 and z=0.5. In the highest density regions the early-type fraction has increased steadily with time from f(E+S0)=0.7+/-0.1 at z=1 to f(E+S0)=0.9+/-0.1 at the present epoch. However, in intermediate density regions corresponding to groups and the accretion regions of rich clusters, significant evolution appears to begin only after z=0.5. Finally, at the lowest densities, no evolution is observed for the early-type fraction of field galaxies which remains constant at f(E+S0)=0.4+/-0.1 at all epochs. We examine a simple picture consistent with these observations where the early-type population at z=1 is comprised largely of elliptical galaxies. Subsequent evolution in both intermediate and dense regions is attributed to the transformation of spirals into lenticulars. Further progress in verifying our hypothesis may be achieved through distinguishing ellipticals and lenticulars at these redshifts through resolved dynamical studies of representative systems.

preprint2003arXiv

Constraints on the equation of state of dark energy and the Hubble constant from stellar ages and the CMB

We place tight constraints on the redshift-averaged, effective value of the equation of state of dark energy, w, using only the absolute ages of Galactic stars and the observed position of the first peak in the angular power spectrum of the CMB. We find w<-0.8 at the 68% confidence level. If we further consider that w > -1, this finding suggests that within our uncertainties, dark energy is indistinguishable from a classical vacuum energy term. We detect a correlation between the ages of the oldest galaxies and their redshift. This opens up the possibility of measuring w(z) by computing the relative ages of the oldest galaxies in the universe as a function of redshift, dz/dt. We show that this is a realistic possibility by computing dz/dt at z~0 from SDSS galaxies and obtain an independent estimate for the Hubble constant, H_0 = 69 \pm 12 km s-1 Mpc-1. The small number of galaxies considered at z>0.2 does not yield, currently, a precise determination of w(z), but shows that the age--redshift relation is consistent with a Standard LCDM universe with $w=-1$.

preprint2003arXiv

The Dark Matter Distribution in the Central Regions of Galaxy Clusters: Implications for CDM

We have undertaken a spectroscopic survey of gravitational arcs in a carefully chosen sample of six clusters each containing a dominant brightest cluster galaxy. We use these systems to study the relative distributions of dark and baryonic material in the central regions. Three clusters present both radial and tangential arcs and provide particularly strong constraints on the mass profiles, whereas the other three display only tangential arcs and act as a control set. We analyze stellar velocity dispersion data for the brightest cluster galaxies in conjunction with the arc redshifts and lens models to constrain the dark and baryonic mass profiles jointly. For the systems containing radial arcs we find that the inner dark matter density profile is consistent with a 3-D distribution, rho_{DM} propto r^-beta, with logarithmic slope <beta>=0.52^{+0.05}_{-0.05} (68% CL). Similarly, we find that the tangential arc sample gives an upper limit, beta<0.57 (99% CL). Taking the 6 clusters together, the mean dark matter distribution is quite inconsistent with the standard Navarro, Frenk & White (1997) value, beta=1.0, at >99% confidence. In addition, we find considerable cosmic scatter in the beta (Delta beta ~0.3) values of the radial arc sample. We find no evidence that systems with radial arcs preferentially yield flatter dark matter profiles as might be expected if they were a biased subset. We discuss the validity of our 1-D mass reconstruction method and verify its conclusions by comparing with results of a more rigorous ray-tracing lensing code. Our results extend and strengthen the earlier conclusions presented by Sand et al. (2002) and suggest the relationship between dark and visible matter in the cores of clusters is more complex than anticipated from recent simulations. (abridged)