The H0LiCOW collaboration inferred via gravitational lensing time delays a Hubble constant $H_0=73.3^{+1.7}_{-1.8}$ km s$^{-1}{\rm Mpc}^{-1}$, describing deflector mass density profiles by either a power-law or stars plus standard dark matter halos. The mass-sheet transform (MST) that leaves the lensing observables unchanged is considered the dominant source of residual uncertainty in $H_0$. We quantify any potential effect of the MST with flexible mass models that are maximally degenerate with H0. Our calculation is based on a new hierarchical approach in which the MST is only constrained by stellar kinematics. The approach is validated on hydrodynamically simulated lenses. We apply the method to the TDCOSMO sample of 7 lenses (6 from H0LiCOW) and measure $H_0=74.5^{+5.6}_{-6.1}$ km s$^{-1}{\rm Mpc}^{-1}$. In order to further constrain the deflector mass profiles, we then add imaging and spectroscopy for 33 strong gravitational lenses from the SLACS sample. For 9 of the SLAC lenses we use resolved kinematics to constrain the stellar anisotropy. From the joint analysis of the TDCOSMO+SLACS sample, we measure $H_0=67.4^{+4.1}_{-3.2}$ km s$^{-1}{\rm Mpc}^{-1}$, assuming that the TDCOSMO and SLACS galaxies are drawn from the same parent population. The blind H0LiCOW, TDCOSMO-only and TDCOSMO+SLACS analyses are in mutual statistical agreement. The TDCOSMO+SLACS analysis prefers marginally shallower mass profiles than H0LiCOW or TDCOSMO-only. While our new analysis does not statistically invalidate the mass profile assumptions by H0LiCOW, and thus their $H_0$ measurement relying on those, it demonstrates the importance of understanding the mass density profile of elliptical galaxies. The uncertainties on $H_0$ derived in this paper can be reduced by physical or observational priors on the form of the mass profile, or by additional data, chiefly spatially resolved kinematics of lens galaxies.

Time delay cosmography uses the arrival time delays between images in strong gravitational lenses to measure cosmological parameters, in particular the Hubble constant $H_0$. The lens models used in time delay cosmography omit dark matter subhalos and line-of-sight halos because their effects are assumed to be negligible. We explicitly quantify this assumption by analyzing realistic mock lens systems that include full populations of dark matter subhalos and line-of-sight halos, applying the same modeling assumptions used in the literature to infer $H_0$. We base the mock lenses on six quadruply-imaged quasars that have delivered measurements of the Hubble constant, and quantify the additional uncertainties and/or bias on a lens-by-lens basis. We show that omitting dark substructure does not bias inferences of $H_0$. However, perturbations from substructure contribute an additional source of random uncertainty in the inferred value of $H_0$ that scales as the square root of the lensing volume divided by the longest time delay. This additional source of uncertainty, for which we provide a fitting function, ranges from $0.6 - 2.4\%$. It may need to be incorporated in the error budget as the precision of cosmographic inferences from single lenses improves, and sets a precision limit on inferences from single lenses.

We study the projected spatial offset between the ultraviolet continuum and Ly$\alpha$ emission for a sample of 65 lensed and unlensed galaxies in the Epoch of Reionization (5<z<7), the first such study at these redshifts, in order to understand the potential for these offsets to confuse estimates of the Ly$\alpha$ properties of galaxies observed in slit spectroscopy. While we find that ~40% of galaxies in our sample show significant projected spatial offsets ($|\Delta_{Ly\alpha-UV}|$), we find a modest average offset of 0.61$\pm$0.08 kpc. A small fraction of our sample, ~10%, exhibits offsets of 2-4 kpc, sizes that are larger than the effective radii of typical galaxies at these redshifts. An internal comparison and a comparison to studies at lower redshift yielded no significant evidence of evolution of $|\Delta_{Ly\alpha-UV}|$ with redshift. In our own sample, UV-bright galaxies showed offsets a factor of three greater than their fainter counterparts, 0.89$\pm$0.18 vs. 0.27$\pm$0.05 kpc, respectively. We argue that offsets are likely not the result of merging processes, but are rather due to anisotropic internal processes internal to the galaxies as a result of stellar feedback, which facilitates Ly$\alpha$ fluorescence and/or backscattering from nearby or outflowing gas. The reduction in the Ly$\alpha$ flux due to offset effects for various observational setups was quantified through mock observations of simple simulations. It was found that the loss of Ly$\alpha$ photons for galaxies with average offsets is not, if corrected for, a limiting factor for all but the narrowest slit widths (<0.4''). However, for the largest offsets, if such offsets are mostly perpendicular to the slit major axis, slit losses were found to be extremely severe in cases where slit widths of $\le$1'' were employed, such as those planned for James Webb Space Telescope/NIRSpec observations. (abridged)

Time-delay cosmography of lensed quasars has achieved 2.4% precision on the measurement of the Hubble constant, $H_0$. As part of an ongoing effort to uncover and control systematic uncertainties, we investigate three potential sources: 1- stellar kinematics, 2- line-of-sight effects, and 3- the deflector mass model. To meet this goal in a quantitative way, we reproduced the H0LiCOW/SHARP/STRIDES (hereafter TDCOSMO) procedures on a set of real and simulated data, and we find the following. First, stellar kinematics cannot be a dominant source of error or bias since we find that a systematic change of 10% of measured velocity dispersion leads to only a 0.7% shift on $H_0$ from the seven lenses analyzed by TDCOSMO. Second, we find no bias to arise from incorrect estimation of the line-of-sight effects. Third, we show that elliptical composite (stars + dark matter halo), power-law, and cored power-law mass profiles have the flexibility to yield a broad range in $H_0$ values. However, the TDCOSMO procedures that model the data with both composite and power-law mass profiles are informative. If the models agree, as we observe in real systems owing to the "bulge-halo" conspiracy, $H_0$ is recovered precisely and accurately by both models. If the two models disagree, as in the case of some pathological models illustrated here, the TDCOSMO procedure either discriminates between them through the goodness of fit, or it accounts for the discrepancy in the final error bars provided by the analysis. This conclusion is consistent with a reanalysis of six of the TDCOSMO (real) lenses: the composite model yields $H_0=74.2^{+1.6}_{-1.6}$ $km.s^{-1}.Mpc^{-1}$, while the power-law model yields $74.0^{+1.7}_{-1.8}$ $km.s^{-1}.Mpc^{-1}$. In conclusion, we find no evidence of bias or errors larger than the current statistical uncertainties reported by TDCOSMO.

We present results from the KMOS Lens-Amplified Spectroscopic Survey (KLASS), an ESO Very Large Telescope (VLT) large program using gravitational lensing to study the spatially resolved kinematics of 44 star-forming galaxies at 0.6<z<2.3 with a stellar mass of 8.1<log(M$_\star$/M$_{\odot}$)<11.0. These galaxies are located behind six galaxy clusters selected from the HST Grism Lens-Amplified Survey from Space (GLASS). We find that the majority of the galaxies show a rotating disk, but most of the rotation-dominated galaxies only have a low $\upsilon_{rot}/\sigma_0$ ratio (median of $\upsilon_{rot}/\sigma_0\sim2.5$). We explore the Tully-Fisher relation by adopting the circular velocity, $V_{circ}=(\upsilon_{rot}^2+3.4\sigma_0^2)^{1/2}$, to account for pressure support. We find that our sample follows a Tully-Fisher relation with a positive zero-point offset of +0.18 dex compared to the local relation, consistent with more gas-rich galaxies that still have to convert most of their gas into stars. We find a strong correlation between the velocity dispersion and stellar mass in the KLASS sample. When combining our data to other surveys from the literature, we also see an increase of the velocity dispersion with stellar mass at all redshift. We obtain an increase of $\upsilon_{rot}/\sigma_0$ with stellar mass at 0.5<z<1.0. This could indicate that massive galaxies settle into regular rotating disks before the low-mass galaxies. For higher redshift (z>1), we find a weak increase or flat trend. We investigate the relation between the rest-frame UV clumpiness of galaxies and their global kinematic properties. We find no clear trend between the clumpiness and the velocity dispersion and $\upsilon_{rot}/\sigma_0$. This could suggest that the kinematic properties of galaxies evolve after the clumps formed in the galaxy disk or that the clumps can form in different physical conditions.

We present six new time-delay measurements obtained from $R_c$-band monitoring data acquired at the MPIA 2.2m telescope at La Silla observatory between October 2016 and February 2020. The lensed quasars HE 0047-1756, WG 0214-2105, DES 0407-5006, 2M 1134-2103, PSJ 1606-2333 and DES 2325-5229 were observed almost daily at high signal-to-noise ratio to obtain high-quality light curves where we can record fast and small-amplitude variations of the quasars. We measure time delays between all pairs of multiple images with only one or two seasons of monitoring at the exception of the time delays relative to image D of PSJ 1606-2333. The most precise estimate is obtained for the delay between image A and image B of DES 0407-5006, with $\tau_{AB} = -128.4^{+3.5}_{-3.8}$ days (2.8% precision) including systematics due to extrinsic variability in the light curves. For HE 0047-1756, we combine our high-cadence data with measurements from decade-long light-curves from previous COSMOGRAIL campaigns, and reach a precision of 0.9 day on the final measurement. The present work demonstrates the feasibility of measuring time delays in lensed quasars in only one or two seasons, provided high signal-to-noise ratio data are obtained at a cadence close to daily.

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

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.

Physical properties of galaxies at z>7 are of interest for understanding both the early phases of star formation and the process of cosmic reionization. Chemical abundance measurements offer valuable information on the integrated star formation history, and hence ionizing photon production, as well as the rapid gas accretion expected at such high redshifts. We use reported measurements of [O III] 88$\mu$m emission and star formation rate to estimate gas-phase oxygen abundances in five galaxies at z=7.1-9.1 using the direct T_e method. We find typical abundances 12+log(O/H) = 7.9 ($\sim$0.2 times the solar value) and an evolution of 0.9$\pm$0.5 dex in oxygen abundance at fixed stellar mass from z$\simeq$8 to 0. These results are compatible with theoretical predictions, albeit with large (conservative) uncertainties in both mass and metallicity. We assess both statistical and systematic uncertainties to identify promising means of improvement with the Atacama Large Millimeter Array (ALMA) and the James Webb Space Telescope (JWST). In particular we highlight [O III] 52$\mu$m as a valuable feature for robust metallicity measurements. Precision of 0.1-0.2 dex in T_e-based O/H abundance can be reasonably achieved for galaxies at z$\approx$5-8 by combining [O III] 52$\mu$m with rest-frame optical strong lines. It will also be possible to probe gas mixing and mergers via resolved T_e-based abundances on kpc scales. With ALMA and JWST, direct metallicity measurements will thus be remarkably accessible in the reionization epoch.

One of the main challenges in using high redshift active galactic nuclei to study the correlations between the mass of the supermassive Black Hole (MBH) and the properties of their active host galaxies is instrumental resolution. Strong lensing magnification effectively increases instrumental resolution and thus helps to address this challenge. In this work, we study eight strongly lensed active galactic nuclei (AGN) with deep Hubble Space Telescope imaging, using the lens modelling code Lenstronomy to reconstruct the image of the source. Using the reconstructed brightness of the host galaxy, we infer the host galaxy stellar mass based on stellar population models. MBH are estimated from broad emission lines using standard methods. Our results are in good agreement with recent work based on non-lensed AGN, providing additional evidence that the correlation evolves over cosmic time. At the moment, the sample size of lensed AGN is small and thus they provide mostly a consistency check on systematic errors related to resolution for the non-lensed AGN. However, the number of known lensed AGN is expected to increase dramatically in the next few years, through dedicated searches in ground and space based wide field surveys, and they may become a key diagnostic of black hole and galaxy co-evolution.

We investigate the environment and line of sight of the H0LiCOW lens B1608+656 using Subaru Suprime-Cam and the Hubble Space Telescope (HST) to perform a weak lensing analysis. We compare three different methods to reconstruct the mass map of the field, i.e. the standard Kaiser-Squires inversion coupled with inpainting and Gaussian or wavelet filtering, and $\tt{Glimpse}$ a method based on sparse regularization of the shear field. We find no substantial difference between the 2D mass reconstructions, but we find that the ground-based data is less sensitive to small-scale structures than the space-based observations. Marginalising over the results obtained with all the reconstruction techniques applied to the two available HST filters F606W and F814W, we estimate the external convergence, $\kappa_{\rm ext}$ at the position of B1608+656 is $\kappa_{\rm ext} = 0.11^{+0.06}_{-0.04}$, where the error bars corresponds respectively to the 16th and 84th quartiles. This result is compatible with previous estimates using the number-counts technique, suggesting that B1608+656 resides in an over-dense line of sight, but with a completely different technique. Using our mass reconstructions, we also compare the convergence at the position of several groups of galaxies in the field of B1608+656 with the mass measurements using various analytical mass profiles, and find that the weak lensing results favor truncated halo models.

We examine the scaling relations between the mass of a supermassive black hole (SMBH) and its host galaxy properties at $1.2<z<1.7$ using both observational data and simulations. Recent measurements of 32 X-ray-selected broad-line Active Galactic Nucleus (AGNs) are compared with two independent state-of-the-art efforts, including the hydrodynamic simulation MassiveBlackII (MBII) and a semi-analytic model (SAM). After applying an observational selection function to the simulations, we find that both MBII and SAM agree well with the data, in terms of the central distribution. However, the dispersion in the mass ratio between black hole mass and stellar mass is significantly more consistent with the MBII prediction ($\sim0.3~$dex), than with the SAM ($\sim0.7~$dex), even when accounting for observational uncertainties. Hence, our observations can distinguish between the different recipes adopted in the models. The mass relations in the MBII are highly dependent on AGN feedback while the relations in the SAM are more sensitive to galaxy merger events triggering nuclear activity. Moreover, the intrinsic scatter in the mass ratio of our high-$z$ sample is comparable to that observed in the local sample, all but ruling out the proposed scenario the correlations are purely stochastic in nature arising from some sort of cosmic central limit theorem. Our results support the hypothesis of AGN feedback being responsible for a causal link between the SMBH and its host galaxy, resulting in a tight correlation between their respective masses.

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 $\Lambda$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 $\sigma$. We find that modifications of $\Lambda$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.

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$\alpha$) 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$\alpha$ 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$\alpha$ visibility during reionization. We show that UV luminosity-halo mass scatter reduces Ly$\alpha$ 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.

We present a blind time-delay cosmographic analysis for the lens system DES J0408$-$5354. This system is extraordinary for the presence of two sets of multiple images at different redshifts, which provide the opportunity to obtain more information at the cost of increased modelling complexity with respect to previously analyzed systems. We perform detailed modelling of the mass distribution for this lens system using three band Hubble Space Telescope imaging. We combine the measured time delays, line-of-sight central velocity dispersion of the deflector, and statistically constrained external convergence with our lens models to estimate two cosmological distances. We measure the "effective" time-delay distance corresponding to the redshifts of the deflector and the lensed quasar $D_{\Delta t}^{\rm eff}=3382^{+146}_{-115}$ Mpc and the angular diameter distance to the deflector $D_{\rm d}=1711^{+376}_{-280}$ Mpc, with covariance between the two distances. From these constraints on the cosmological distances, we infer the Hubble constant $H_0 = 74.2^{+2.7}_{-3.0}$ km s$^{-1}$ Mpc$^{-1}$ assuming a flat $\Lambda$CDM cosmology and a uniform prior for $\Omega_{\rm m}$ as $\Omega_{\rm m} \sim \mathcal{U}(0.05, 0.5)$. This measurement gives the most precise constraint on $H_0$ to date from a single lens. Our measurement is consistent with that obtained from the previous sample of six lenses analyzed by the $H_0$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration. It is also consistent with measurements of $H_0$ based on the local distance ladder, reinforcing the tension with the inference from early Universe probes, for example, with 2.2$\sigma$ discrepancy from the cosmic microwave background measurement.

In time-delay cosmography, three of the key ingredients are 1) determining the velocity dispersion of the lensing galaxy, 2) identifying galaxies and groups along the line of sight with sufficient proximity and mass to be included in the mass model, and 3) estimating the external convergence $\kappa_\mathrm{ext}$ from less massive structures that are not included in the mass model. We present results on all three of these ingredients for two time-delay lensed quasar systems, DES J0408-5354 and WGD 2038-4008. We use the Gemini, Magellan and VLT telescopes to obtain spectra to both measure the stellar velocity dispersions of the main lensing galaxies and to identify the line-of-sight galaxies in these systems. Next, we identify 10 groups in DES J0408-5354 and 2 groups in WGD 2038-4008using a group-finding algorithm. We then identify the most significant galaxy and galaxy-group perturbers using the "flexion shift" criterion. We determine the probability distribution function of the external convergence $\kappa_\mathrm{ext}$ for both of these systems based on our spectroscopy and on the DES-only multiband wide-field observations. Using weighted galaxy counts, calibrated based on the Millennium Simulation, we find that DES J0408-5354 is located in a significantly underdense environment, leading to a tight (width $\sim3\%$), negative-value $\kappa_\mathrm{ext}$ distribution. On the other hand, WGD 2038-4008 is located in an environment of close to unit density, and its low source redshift results in a much tighter $\kappa_\mathrm{ext}$ of $\sim1\%$, as long as no external shear constraints are imposed.

We discovered a strongly lensed (\mu >40) Lya emission at z=6.629 (S/N~18) in the MUSE Deep Lensed Field (MDLF) targeting the Hubble Frontier Field galaxy cluster MACS~J0416. Dedicated lensing simulations imply that the Lya emitting region necessarily crosses the caustic. The arc-like shape of the Lya extends 3 arcsec on the observed plane and is the result of two merged multiple images, each one with a de-lensed Lya luminosity L<~2.8 x 10^(40) erg/s arising from a confined region (< 150 pc effective radius). A spatially unresolved HST counterpart is barely detected at S/N~2 after stacking the near-infrared bands, corresponding to an observed(intrinsic) magnitude m_(1500)>~30.8(>~35.0). The inferred rest-frame Lya equivalent width is EWo > 1120 A if the IGM transmission is T(IGM)<0.5. The low luminosities and the extremely large Lya EWo match the case of a Population~III star complex made of several dozens stars (~ 10^4 Msun) which irradiate a HII region crossing the caustic. While the Lya and stellar continuum are among the faintest ever observed at this redshift, the continuum and the Lya emissions could be affected by differential magnification, possibly biasing the EWo estimate. The aforementioned tentative HST detection tend to favor a large EWo, making such a faint Pop~III candidate a key target for the James Webb Space Telescope and Extremely Large Telescopes.

We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon.

We present 36 spectroscopically confirmed intrinsically UV-faint Ly$\alpha$ 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$\alpha$ 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$\alpha$ emission and the photometric redshift of the non-detections. Considering only objects where Ly$\alpha$ could be detected at EW(Ly$\alpha$)$>$25{\AA} at $3\sigma$ at the fiducial depth of our survey, and only those galaxies with EW(Ly$\alpha$)$>$25{\AA} 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$\alpha$ 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.

The structure of the broad line region (BLR) is an essential ingredient in the determination of active galactic nuclei (AGN) virial black hole masses, which in turn are important to study the role of black holes in galaxy evolution. Constraints on the BLR geometry and dynamics can be obtained from velocity-resolved studies using reverberation mapping data (i.e. monitoring data). However, monitoring data are observationally expensive and only available for a limited sample of AGN, mostly confined to the local Universe. Here we explore a new version of a Bayesian inference, physical model of the BLR which uses an individual spectrum and prior information on the BLR size from the radius-luminosity relation, to model the AGN BLR geometry and dynamics. We apply our model to a sample of 11 AGN, which have been previously modelled using monitoring data. Our single-epoch BLR model is able to constrain some of the BLR parameters with inferred parameter values that agree within the uncertainties with those determined from the modelling of monitoring data. We find that our model is able to derive stronger constraints on the BLR for AGN with broad emission lines that qualitatively have more substructure and more asymmetry, presumably as they contain more information to constrain the physical model. The performance of this model makes it a practical and cost-effective tool to determine some of the BLR properties of a large sample of low and high redshift AGN, for which monitoring data are not available.

We describe a general purpose method to reconstruct the intrinsic properties of sources lensed by the gravitational potential of foreground clusters of galaxies. The tool Lenstruction is implemented in the publicly available multi-purpose gravitational lensing software Lenstronomy, in order to provide easy and fast solution to this common astrophysical problem. The tool is based on forward modeling the appearance of the source in the image plane, taking into account the distortion by lensing and the instrumental point spread function (PSF). For singly-imaged sources a global lens model in the format of the Hubble Frontier Fields (HFF) lensing maps is required as a starting point. For multiply-imaged sources, the tool can also fit and apply first (deflection) second (shear, convergence) and third order (flexion) corrections to the local gravitational potential to improve the overall reconstruction, depending on the quality of the data. We illustrate the performance and features of the code with two examples of multiply-imaged systems taken from the Hubble Frontier Fields, starting from five different publicly available cluster models. We find that, after our correction, the relative magnification - and other lensing properties - between the multiple images become robustly constrained. Furthermore, we find that scatter between models of the reconstructed source size and magnitude is reduced. The code and jupyter notebooks are publicly available.

We report the results of the STRong lensing Insights from the Dark Energy Survey (STRIDES) follow-up campaign of the late 2017/early 2018 season. We obtained spectra of 65 lensed quasar candidates either with EFOSC2 on the NTT or ESI on Keck, which confirm 10 new gravitationally lensed quasars and 10 quasar pairs with similar spectra, but which do not show a lensing galaxy in DES images. Eight lensed quasars are doubly imaged with source redshifts between 0.99 and 2.90, one is triply imaged by a group (DESJ0345-2545, $z=1.68$), and one is quadruply imaged (quad: DESJ0053-2012, $z=3.8$). Singular isothermal ellipsoid models for the doubles, based on high-resolution imaging from SAMI on SOAR or NIRC2 on Keck, give total magnifications between 3.2 and 5.6, and Einstein radii between 0.49 and 1.97 arcseconds. After spectroscopic follow-up, we extract multi-epoch $grizY$ photometry of confirmed lensed quasars and contaminant quasar+star pairs from the first 4 years of DES data using parametric multi-band modelling, and compare variability in each system's components. By measuring the reduced ${\chi}^2$ associated with fitting all epochs to the same magnitude, we find a simple cut on the less variable component that retains all confirmed lensed quasars, while removing 94 per cent of contaminant systems with stellar components. Based on our spectroscopic follow-up, this variability information can improve selection of lensed quasars and quasar pairs from 34-45 per cent to 51-70 per cent, with the majority of remaining contaminants being compact star-forming galaxies. Using mock lensed quasar lightcurves we demonstrate that selection based only on variability will over-represent the quad fraction by 10 per cent over a complete DES magnitude-limited sample (excluding microlensing differences), explained by the magnification bias and hence lower luminosity (more variable) sources in quads.

The magnifications of compact-source lenses are extremely sensitive to the presence of low mass dark matter halos along the entire sight line from the source to the observer. Traditionally, the study of dark matter structure in compact-source strong gravitational lenses has been limited to radio-loud systems, as the radio emission is extended and thus unaffected by microlensing which can mimic the signal of dark matter structure. An alternate approach is to measure quasar nuclear-narrow line emission, which is free from microlensing and present in virtually all quasar lenses. In this paper, we double the number of systems which can be used for gravitational lensing analyses by presenting measurements of narrow-line emission from a sample of 8 quadruply imaged quasar lens systems, WGD J0405-3308, HS 0810+2554, RX J0911+0551, SDSS J1330+1810, PS J1606-2333, WFI 2026-4536, WFI 2033-4723 and WGD J2038-4008. We describe our updated grism spectral modelling pipeline, which we use to measure narrow-line fluxes with uncertainties of 2-10\%, presented here. We fit the lensed image positions with smooth mass models and demonstrate that these models fail to produce the observed distribution of image fluxes over the entire sample of lenses. Furthermore, typical deviations are larger than those expected from macromodel uncertainties. This discrepancy indicates the presence of perturbations caused by small-scale dark matter structure. The interpretation of this result in terms of dark matter models is presented in a companion paper.

We present the lens mass model of the quadruply-imaged gravitationally lensed quasar WFI2033-4723, and perform a blind cosmographical analysis based on this system. Our analysis combines (1) time-delay measurements from 14 years of data obtained by the COSmological MOnitoring of GRAvItational Lenses (COSMOGRAIL) collaboration, (2) high-resolution $\textit{Hubble Space Telescope}$ imaging, (3) a measurement of the velocity dispersion of the lens galaxy based on ESO-MUSE data, and (4) multi-band, wide-field imaging and spectroscopy characterizing the lens environment. We account for all known sources of systematics, including the influence of nearby perturbers and complex line-of-sight structure, as well as the parametrization of the light and mass profiles of the lensing galaxy. After unblinding, we determine the effective time-delay distance to be $4784_{-248}^{+399}~\mathrm{Mpc}$, an average precision of $6.6\%$. This translates to a Hubble constant $H_{0} = 71.6_{-4.9}^{+3.8}~\mathrm{km~s^{-1}~Mpc^{-1}}$, assuming a flat $\Lambda$CDM cosmology with a uniform prior on $\Omega_\mathrm{m}$ in the range [0.05, 0.5]. This work is part of the $H_0$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration, and the full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper (H0LiCOW XIII).

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 79 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-$\sigma$ confidence level, 14/79 galaxies in our sample show negative radial gradients, whereas 7/79 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 $\Delta\log({\rm O/H})/\Delta r~ [\mathrm{dex~kpc^{-1}}] = (-0.020\pm0.007) + (-0.014\pm0.008) \log(M_\ast/10^{9.4}M_\odot)$ with the intrinsic scatter being $\sigma=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$.

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.

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.

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

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.

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 $\eta$ is a power-law function of galaxy stellar mass, $\eta \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.

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.

Galaxies and galaxy groups located along the line of sight towards gravitationally lensed quasars produce high-order perturbations of the gravitational potential at the lens position. When these perturbation are too large, they can induce a systematic error on $H_0$ of a few-percent if the lens system is used for cosmological inference and the perturbers are not explicitly accounted for in the lens model. In this work, we present a detailed characterization of the environment of the lens system WFI2033-4723 ($z_{\rm src} = 1.662$, $z_{\rm lens}$ = 0.6575), one of the core targets of the H0LICOW project for which we present cosmological inferences in a companion paper (Rusu et al. 2019). We use the Gemini and ESO-Very Large telescopes to measure the spectroscopic redshifts of the brightest galaxies towards the lens, and use the ESO-MUSE integral field spectrograph to measure the velocity-dispersion of the lens ($\sigma_{\rm {los}}= 250^{+15}_{-21}$ km/s) and of several nearby galaxies. In addition, we measure photometric redshifts and stellar masses of all galaxies down to $i < 23$ mag, mainly based on Dark Energy Survey imaging (DR1). Our new catalog, complemented with literature data, more than doubles the number of known galaxy spectroscopic redshifts in the direct vicinity of the lens, expanding to 116 (64) the number of spectroscopic redshifts for galaxies separated by less than 3 arcmin (2 arcmin) from the lens. Using the flexion-shift as a measure of the amplitude of the gravitational perturbation, we identify 2 galaxy groups and 3 galaxies that require specific attention in the lens models. The ESO MUSE data enable us to measure the velocity-dispersions of three of these galaxies. These results are essential for the cosmological inference analysis presented in Rusu et al. (2019).

The time delay between the arrival of photons of multiple images of time variable sources can be used to constrain absolute distances in the Universe (Refsdal 1964), and in turn obtain a direct estimate of the Hubble constant and other cosmological parameters. To convert the time delay into distances, it is well known that the gravitational potential of the main deflector and the contribution of the matter along the line-of-sight need to be known to a sufficient level of precision. In this paper, we discuss a new astrometric requirement that is becoming important as time-delay cosmography improves in precision and accuracy with larger samples, and better data and modelling techniques. We derive an analytic expression for the propagation of astrometric uncertainties on the multiple image positions into the inference of the Hubble constant and derive requirements depending on image separation and relative time delay. We note that this requirement applies equally to the image position measurements and to the accuracy of the model in reproducing them. To illustrate the requirement, we discuss some example lensing configurations and highlight that, especially for time delays of order 10 days or shorter, the relative astrometric requirement is of order milli-arcseconds, setting a tight requirement on both measurements and models. With current optical infrared technology, astrometric uncertainties may be the dominant limitation for strong lensing cosmography in the small image-separation regime when high-precision time-delays become accessible.

We report the first sub-kiloparsec spatial resolution measurements of strongly inverted gas-phase metallicity gradients in two dwarf galaxies at $z$$\sim$2. The galaxies have stellar masses $\sim$$10^9M_\odot$, specific star-formation rate $\sim$20 Gyr$^{-1}$, and global metallicity $12+\log({\rm O/H})\sim8.1$ (1/4 solar), assuming the Maiolino et al. (2008) strong line calibrations of OIII/Hb and OII/Hb. Their metallicity radial gradients are measured to be highly inverted, i.e., 0.122$\pm$0.008 and 0.111$\pm$0.017 dex/kpc, which is hitherto unseen at such small masses in similar redshift ranges. From the Hubble Space Telescope observations of the source nebular emission and stellar continuum, we present the 2-dimensional spatial maps of star-formation rate surface density, stellar population age, and gas fraction, which show that our galaxies are currently undergoing rapid mass assembly via disk inside-out growth. More importantly, using a simple chemical evolution model, we find that the gas fractions for different metallicity regions cannot be explained by pure gas accretion. Our spatially resolved analysis based on a more advanced gas regulator model results in a spatial map of net gaseous outflows, triggered by active central starbursts, that potentially play a significant role in shaping the spatial distribution of metallicity by effectively transporting stellar nucleosynthesis yields outwards. The relation between wind mass loading factors and stellar surface densities measured in different regions of our galaxies shows that a single type of wind mechanism, driven by either energy or momentum conservation, cannot explain the entire galaxy. These sources present a unique constraint on the effects of gas flows on the early phase of disk growth from the perspective of spatially resolved chemical evolution within individual systems.

This paper presents multiwavelength photometric catalogues of the last two Hubble Frontier Fields (HFF), the massive galaxy clusters Abell 370 and RXC J2248.7-4431. The photometry ranges from imaging performed on the Hubble Space Telescope (HST) to ground based Very Large Telescope (VLT) and Spitzer/IRAC, in collaboration with the ASTRODEEP team, and using the ASTRODEEP pipeline. While the main purpose of this paper is to release the catalogues, we also perform, as a proof of concept, a brief analysis of z > 6 objects selected using drop-out method, as well as spectroscopically confirmed sources and multiple images in both clusters. While dropout methods yield a sample of high-z galaxies, the addition of longer wavelength data reveals that as expected the samples have substantial contamination at the ~30-45% level by dusty galaxies at lower redshifts. Furthermore, we show that spectroscopic redshifts are still required to unambiguously determine redshifts of multiply imaged systems. Finally, the now publicly available ASTRODEEP catalogues were combined for all HFFs and used to explore stellar properties of a large sample of 20,000 galaxies across a large photometric redshift range. The powerful magnification provided by the HFF clusters allows us an exploration of the properties of galaxies with intrinsic stellar masses as low as $M_* \gtrsim 10^7M_{\odot}$ and intrinsic star formation rates $\mbox{SFRs}\sim 0.1\mbox{-}1M_\odot/\mbox yr$ at z > 6.

Strong lensing provides a powerful means of investigating the nature of dark matter as it probes dark matter structure on sub-galactic scales. We present an extension of a forward modeling framework that uses flux ratios from quadruply imaged quasars (quads) to measure the shape and amplitude of the halo mass function, including line of sight (LOS) halos and main deflector subhalos. We apply this machinery to 50 mock lenses --- roughly the number of known quads --- with warm dark matter (WDM) mass functions exhibiting free-streaming cutoffs parameterized by the half-mode mass $m_{\rm{hm}}$. Assuming cold dark matter (CDM), we forecast bounds on $m_{\rm{hm}}$ and the corresponding thermal relic particle masses over a range of tidal destruction severity, assuming a particular WDM mass function and mass-concentration relation. With significant tidal destruction, at $2 \sigma$ we constrain $m_{\rm{hm}}<10^{7.9} \left(10^{8.4}\right) M_{\odot}$, or a 4.4 (3.1) keV thermal relic, with image flux uncertainties from measurements and lens modeling of $2\% \left(6\%\right)$. With less severe tidal destruction we constrain $m_{\rm{hm}}<10^{7} \left(10^{7.4}\right) M_{\odot}$, or an 8.2 (6.2) keV thermal relic. If dark matter is warm, with $m_{\rm{hm}} = 10^{7.7} M_{\odot}$ (5.1 keV), we would favor WDM with $m_{\rm{hm}} > 10^{7.7} M_{\odot}$ over CDM with relative likelihoods of 22:1 and 8:1 with flux uncertainties of $2\%$ and $6\%$, respectively. These bounds improve over those obtained by modeling only main deflector subhalos because LOS objects produce additional flux perturbations, especially for high redshift systems. These results indicate that $\sim 50$ quads can conclusively differentiate between warm and cold dark matter.

The standard cosmological model successfully describes many observations from widely different epochs of the Universe, from primordial nucleosynthesis to the accelerating expansion of the present day. However, as the basic cosmological parameters of the model are being determined with increasing and unprecedented precision, it is not guaranteed that the same model will fit more precise observations from widely different cosmic epochs. Discrepancies developing between observations at early and late cosmological time may require an expansion of the standard model, and may lead to the discovery of new physics. The workshop "Tensions between the Early and the Late Universe" was held at the Kavli Institute for Theoretical Physics on July 15-17 2019 (More details of the workshop (including on-line presentations) are given at the website: https://www.kitp.ucsb.edu/activities/enervac-c19) to evaluate increasing evidence for these discrepancies, primarily in the value of the Hubble constant as well as ideas recently proposed to explain this tension. Multiple new observational results for the Hubble constant were presented in the time frame of the workshop using different probes: Cepheids, strong lensing time delays, tip of the red giant branch (TRGB), megamasers, Oxygen-rich Miras and surface brightness fluctuations (SBF) resulting in a set of six new ones in the last several months. Here we present the summary plot of the meeting that shows combining any three independent approaches to measure H$_0$ in the late universe yields tension with the early Universe values between 4.0$\sigma$ and 5.8$\sigma$. This shows that the discrepancy does not appear to be dependent on the use of any one method, team, or source. Theoretical ideas to explain the discrepancy focused on new physics in the decade of expansion preceding recombination as the most plausible. This is a brief summary of the workshop.

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state.

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

Fast Radio Bursts (FRBs) are brief radio emissions from distant astronomical sources. Some are known to repeat, but most are single bursts. Non-repeating FRB observations have had insufficient positional accuracy to localize them to an individual host galaxy. We report the interferometric localization of the single pulse FRB 180924 to a position 4 kpc from the center of a luminous galaxy at redshift 0.3214. The burst has not been observed to repeat. The properties of the burst and its host are markedly different from the only other accurately localized FRB source. The integrated electron column density along the line of sight closely matches models of the intergalactic medium, indicating that some FRBs are clean probes of the baryonic component of the cosmic web.

Observations have revealed massive (logM*/Msun>11) galaxies that were already dead when the universe was only ~2 Gyr. Given the short time before these galaxies were quenched, their past histories and quenching mechanism(s) are of particular interest. In this paper, we study star formation histories (SFHs) of 24 massive galaxies at 1.6<z<2.5. A deep slitless spectroscopy + imaging data set collected from multiple Hubble Space Telescope surveys allows robust determination of their spectral energy distributions and SFHs with no functional assumption on their forms. We find that most of our massive galaxies had formed > 50% of their extant masses by ~1.5 Gyr before the time of observed redshifts, with a trend where more massive galaxies form earlier. Their stellar-phase metallicities are already compatible with those of local early-type galaxies, with a median value of logZ*/Zsun=0.25 and scatter of ~0.15dex. In combination with the reconstructed SFHs, we reveal their rapid metallicity evolution from z~5.5 to ~2.2 at a rate of ~0.2dex/Gyr in log Z*/Zsun. Interestingly, the inferred stellar-phase metallicities are, when compared at half-mass time, ~0.25dex higher than observed gas-phase metallicities of star forming galaxies. While systematic uncertainties remain, this may imply that these quenched galaxies have continued low-level star formation, rather than abruptly terminating their star formation activity, and kept enhancing their metallicity until recently.

We present a catalogue of 22755 objects with slitless, optical, Hubble Space Telescope (HST) spectroscopy from the Grism Lens-Amplified Survey from Space (GLASS). The data cover $\sim$220 sq. arcmin to 7-orbit ($\sim$10 ks) depth in 20 parallel pointings of the Advanced Camera for Survey's G800L grism. The fields are located 6' away from 10 massive galaxy clusters in the HFF and CLASH footprints. Thirteen of the fields have ancillary HST imaging from these or other programs to facilitate a large number of applications, from studying metal distributions at $z\sim0.5$, to quasars at $z\sim4$, to the star formation histories of hundreds of galaxies in between. The spectroscopic catalogue has a median redshift of $\langle z\rangle=0.6$ with a median uncertainty of $\Delta z / (1+z)\lesssim2\%$ at $\rm F814W\lesssim23$ AB. Robust continuum detections reach a magnitude fainter. The 5 $\sigma$ limiting line flux is $f_{\rm lim}\approx5\times10^{-17}\rm~erg~s^{-1}~cm^{-2}$ and half of all sources have 50% of pixels contaminated at $\lesssim$1%. All sources have 1- and 2-D spectra, line fluxes/uncertainties and identifications, redshift probability distributions, spectral models, and derived narrow-band emission line maps from the Grism Redshift and Line Analysis tool (GRIZLI). We provide other basic sample characterisations, show data examples, and describe sources and potential investigations of interest. All data and products will be available online along with software to facilitate their use.

We present a blind time-delay strong lensing (TDSL) cosmographic analysis of the doubly imaged quasar SDSS 1206+4332. We combine the relative time delay between the quasar images, Hubble Space Telescope imaging, the Keck stellar velocity dispersion of the lensing galaxy, and wide-field photometric and spectroscopic data of the field to constrain two angular diameter distance relations. The combined analysis is performed by forward modelling the individual data sets through a Bayesian hierarchical framework, and it is kept blind until the very end to prevent experimenter bias. After unblinding, the inferred distances imply a Hubble constant $H_0 = 68.8^{+5.4}_{-5.1}$ kms$^{-1}$Mpc$^{-1}$, assuming a flat Lambda cold dark matter cosmology with uniform prior on $\Omega_{\rm m}$ in [0.05, 0.5]. The precision of our cosmographic measurement with the doubly imaged quasar SDSS 1206+4332 is comparable with those of quadruply imaged quasars and opens the path to perform on selected doubles the same analysis as anticipated for quads. Our analysis is based on a completely independent lensing code than our previous three H0LiCOW systems and the new measurement is fully consistent with those. We provide the analysis scripts paired with the publicly available software to facilitate independent analysis. The consistency between blind measurements with independent codes provides an important sanity check on lens modelling systematics. By combining the likelihoods of the four systems under the same prior, we obtain $H_0 = 72.5^{+2.1}_{-2.3}$kms$^{-1}$Mpc$^{-1}$. This measurement is independent of the distance ladder and other cosmological probes.

We constrain the distribution of spatially offset Lyman-alpha emission (Ly$\alpha$) relative to rest-frame ultraviolet emission in $\sim300$ high redshift ($3<z<5.5$) Lyman-break galaxies (LBGs) exhibiting Ly$\alpha$ emission from VANDELS, a VLT/VIMOS slit-spectroscopic survey of the CANDELS Ultra Deep Survey and Chandra Deep Field South fields (${\simeq0.2}~\mathrm{deg}^2$ total). Because slit spectroscopy compresses two-dimensional spatial information into one spatial dimension, we use Bayesian inference to recover the underlying Ly$\alpha$ spatial offset distribution. We model the distribution using a 2D circular Gaussian, defined by a single parameter $\sigma_{r,\mathrm{Ly}\alpha}$, the standard deviation expressed in polar coordinates. Over the entire redshift range of our sample ($3<z<5.5$), we find $\sigma_{r,\mathrm{Ly}\alpha}=1.70^{+0.09}_{-0.08}$ kpc ($68\%$ conf.), corresponding to $\sim0.25$ arcsec at $\langle z\rangle=4.5$. We also find that $\sigma_{r,\mathrm{Ly}\alpha}$ decreases significantly with redshift. Because Ly$\alpha$ spatial offsets can cause slit-losses, the decrease in $\sigma_{r,\mathrm{Ly}\alpha}$ with redshift can partially explain the increase in the fraction of Ly$\alpha$ emitters observed in the literature over this same interval, although uncertainties are still too large to reach a strong conclusion. If $\sigma_{r,\mathrm{Ly}\alpha}$ continues to decrease into the reionization epoch, then the decrease in Ly$\alpha$ transmission from galaxies observed during this epoch might require an even higher neutral hydrogen fraction than what is currently inferred. Conversely, if spatial offsets increase with the increasing opacity of the IGM, slit losses may explain some of the drop in Ly$\alpha$ transmission observed at $z>6$. Spatially resolved observations of Ly$\alpha$ and UV continuum at $6<z<8$ are needed to settle the issue.

We present new measurements of the time delays of WFI2033-4723. The data sets used in this work include 14 years of data taken at the 1.2m Leonhard Euler Swiss telescope, 13 years of data from the SMARTS 1.3m telescope at Las Campanas Observatory and a single year of high-cadence and high-precision monitoring at the MPIA 2.2m telescope. The time delays measured from these different data sets, all taken in the R-band, are in good agreement with each other and with previous measurements from the literature. Combining all the time-delay estimates from our data sets results in Dt_AB = 36.2-0.8+0.7 days (2.1% precision), Dt_AC = -23.3-1.4+1.2 days (5.6%) and Dt_BC = -59.4-1.3+1.3 days (2.2%). In addition, the close image pair A1-A2 of the lensed quasars can be resolved in the MPIA 2.2m data. We measure a time delay consistent with zero in this pair of images. We also explore the prior distributions of microlensing time-delay potentially affecting the cosmological time-delay measurements of WFI2033-4723. There is however no strong indication in our measurements that microlensing time delay is neither present nor absent. This work is part of a H0LiCOW series focusing on measuring the Hubble constant from WFI2033-4723.

We introduce the OSIRIS Lens-Amplified Survey (OLAS), a kinematic survey of gravitationally lensed galaxies at cosmic noon taken with Keck adaptive optics. In this paper we present spatially resolved spectroscopy and nebular emission kinematic maps for 17 star forming galaxies with stellar masses 8 < log($M_*$/$M_{\odot}$) < 9.8 and redshifts 1.2 < z < 2.3. OLAS is designed to probe the stellar mass ($M_*$) and specific star formation rate (sSFR) range where simulations suggest that stellar feedback is most effective at driving gaseous outflows that create galaxy-wide potential fluctuations which can generate dark matter cores. We compare our kinematic data with the trend between sSFR, $M_*$ and H$\alpha$ velocity dispersion, $\sigma$, from the Feedback In Realistic Environments (FIRE) simulations. Our observations reveal a correlation between sSFR and sigma at fixed $M_*$ that is similar to the trend predicted by simulations: feedback from star formation drives star-forming gas and newly formed stars into more dispersion dominated orbits. The observed magnitude of this effect is in good agreement with the FIRE simulations, in which feedback alters the central density profiles of low mass galaxies, converting dark matter cusps into cores over time. Our data support the scenario that stellar feedback drives gaseous outflows and potential fluctuations, which in turn drive dark matter core formation in dwarf galaxies.

We have compelling evidence for stellar-mass black holes (BHs) of ~5-80 M_sun that form through the death of massive stars. We also have compelling evidence for so-called supermassive BHs (10^5-10^10 M_sun) that are predominantly found in the centers of galaxies. We have very good reason to believe there must be BHs with masses in the gap between these ranges: the first ~10^9 M_sun BHs are observed only hundreds of millions of years after the Big Bang, and all theoretically viable paths to making supermassive BHs require a stage of "intermediate" mass. However, no BHs have yet been reliably detected in the 100-10}^5 M_sun mass range. Uncovering these intermediate-mass BHs of 10^3-10^5 M_sun is within reach in the coming decade. In this white paper we highlight the crucial role that 30-m class telescopes will play in dynamically detecting intermediate-mass black holes, should they exist.

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

This document illustrates the feasibility of a few per cent level measurement of the local black hole occupation fraction in low mass galaxies through wide-field, high angular resolution X-ray imaging observations of local volume galaxies. The occupation fraction, particularly at the low end of the galaxy luminosity function, is a key benchmark for any model which aims to reproduce the formation and growth of super-massive black holes and their host galaxies. Our proposed measurement will complement orthogonal efforts that are planned in X-rays at high red-shifts, as well as in the local Universe with ground-based facilities.

Many of the fundamental physical constants in Physics, as a discipline, are measured to exquisite levels of precision. The fundamental constants that define Cosmology, however, are largely determined via a handful of independent techniques that are applied to even fewer datasets. The history of the measurement of the Hubble Constant (H0), which serves to anchor the expansion history of the Universe to its current value, is an exemplar of the difficulties of cosmological measurement; indeed, as we approach the centennial of its first measurement, the quest for H0 still consumes a great number of resources. In this white paper, we demonstrate how the approaching era of Extremely Large Telescopes (ELTs) will transform the astrophysical measure of H0 from the limited and few into a fundamentally new regime where (i) multiple, independent techniques are employed with modest use of large aperture facilities and (ii) 1% or better precision is readily attainable. This quantum leap in how we approach H0 is due to the unparalleled sensitivity and spatial resolution of ELT's and the ability to use integral field observations for simultaneous spectroscopy and photometry, which together permit both familiar and new techniques to effectively by-pass the conventional 'ladder' framework to minimize total uncertainty. Three independent techniques are discussed -- (i) standard candles via a two-step distance ladder applied to metal, poor stellar populations, (ii) standard clocks via gravitational lens cosmography, and (iii) standard sirens via gravitational wave sources -- each of which can reach 1% with relatively modest investment from 30-m class facilities.

For nearly 40 years, dark matter has been widely assumed to be cold and collisionless. Cold dark matter models make fundamental predictions for the behavior of dark matter on small (<10 kpc) scales. These predictions include cuspy density profiles at the centers of dark matter halos and a halo mass function that increases as dN/dM ~ M^-1.9 down to very small masses. We suggest two observational programs relying on extremely large telescopes to critically test these predictions, and thus shed new light on the nature of dark matter. (1) Combining adaptive optics-enabled imaging with deep spectroscopy to measure the three-dimensional motions of stars within a sample of Local Group dwarf galaxies that are the cleanest dark matter laboratories known in the nearby universe. From these observations the inner slope of the dark matter density profile can be determined with an accuracy of 0.20 dex, enabling a central cusp to be distinguished from a core at 5 sigma significance. (2) Diffraction-limited AO imaging and integral field spectroscopy of gravitationally lensed galaxies and quasars to quantify the abundance of dark substructures in the halos of the lens galaxies and along the line of sight. Observations of 50 lensed arcs and 50 multiply-imaged quasars will be sufficient to measure the halo mass function over the range 10^7 < M < 10^10 Msun at cosmological scales, independent of the baryonic and stellar composition of those structures. These two observational probes provide complementary information about the small scale structure, with a joint self-consistent analysis mitigating limitations of either probe. This program will produce the strongest existing constraints on the properties of dark matter on small scales, allowing conclusive tests of alternative warm, fuzzy, and self-interacting dark matter models.

The rest-frame UV emission from massive stars contains a wealth of information about the physical nature and conditions of star formation in galaxies. Using studies of the rest-frame UV, the past decade has witnessed the beginning of knowledge about the existence and properties of galaxies during the first few billion years after the Big Bang. This period of history corresponds to the formation of the first stars, the rapid formation of galaxy stellar populations, the reionization of the IGM, the production and dissemination of heavy elements, and the formation of the first black holes. Massive stars in these galaxies drive all of these events, and their light dominates the spectral energy distributions of galaxies. As we look to the 2020s, fundamental questions remain about the nature of these stellar populations and their evolution, from just before the peak of the cosmic star formation density (z~3), up to the epoch of reionization (z > 6). This next decade will provide transformative gains both in our ability to identify star-forming galaxies and accreting supermassive black holes at these early epochs with imaging surveys in the rest-frame UV (e.g., LSST, WFIRST). Ground-based, rest-frame UV spectroscopy on >20 m-class telescopes (e.g., GMT/TMT) offers the ability to investigate the astrophysical conditions in galaxies at the earliest cosmic times. This includes studies of the evolution in galaxy stellar populations, gas ionization (temperature, pressure), metallicity, and interstellar (and circumgalactic) gas kinematics and covering fractions. In this white paper, we describe the scientific prospects and the requirements for research in this area.

The discovery of luminous quasars at redshifts up to 7.5 demonstrates the existence of several billion M_sun supermassive black holes (SMBHs) less than a billion years after the Big Bang. They are accompanied by intense star formation in their host galaxies, pinpointing sites of massive galaxy assembly in the early universe, while their absorption spectra reveal an increasing neutral intergalactic medium (IGM) at the epoch of reionization. Extrapolating from the rapid evolution of the quasar density at z=5-7, we expect that there is only one luminous quasar powered by a billion M_sun SMBH in the entire observable universe at z~9. In the next decade, new wide-field, deep near-infrared (NIR) sky surveys will push the redshift frontier to the first luminous quasars at z~9-10; the combination with new deep X-ray surveys will probe fainter quasar populations that trace earlier phases of SMBH growth. The identification of these record-breaking quasars, and the measurements of their BH masses and accretion properties require sensitive spectroscopic observations with next generation of ground-based and space telescopes at NIR wavelengths. High-resolution integral-field spectroscopy at NIR, and observations at millimeter and radio wavelengths, will together provide a panchromatic view of the quasar host galaxies and their galactic environment at cosmic dawn, connecting SMBH growth with the rise of the earliest massive galaxies. Systematic surveys and multiwavelength follow-up observations of the earliest luminous quasars will strongly constrain the seeding and growth of the first SMBHs in the universe, and provide the best lines of sight to study the history of reionization.

Detections and non-detections of Lyman alpha (Ly$\alpha$) emission from $z>6$ galaxies ($<1$ Gyr after the Big Bang) can be used to measure the timeline of cosmic reionization. Of key interest to measuring reionization's mid-stages, but also increasing observational challenge, are observations at z > 7, where Ly$\alpha$ redshifts to near infra-red wavelengths. Here we present a search for z > 7.2 Ly$\alpha$ emission in 53 intrinsically faint Lyman Break Galaxy candidates, gravitationally lensed by massive galaxy clusters, in the KMOS Lens-Amplified Spectroscopic Survey (KLASS). With integration times of ~7-10 hours, we detect no Ly$\alpha$ emission with S/N>5 in our sample. We determine our observations to be 80% complete for 5$\sigma$ spatially and spectrally unresolved emission lines with integrated line flux $>5.7\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$. We define a photometrically selected sub-sample of 29 targets at $z=7.9\pm0.6$, with a median 5$\sigma$ Ly$\alpha$ EW limit of 58A. We perform a Bayesian inference of the average intergalactic medium (IGM) neutral hydrogen fraction using their spectra. Our inference accounts for the wavelength sensitivity and incomplete redshift coverage of our observations, and the photometric redshift probability distribution of each target. These observations, combined with samples from the literature, enable us to place a lower limit on the average IGM neutral hydrogen fraction of $> 0.76 \; (68\%), \; > 0.46 \; (95\%)$ at z ~ 8, providing further evidence of rapid reionization at z~6-8. We show that this is consistent with reionization history models extending the galaxy luminosity function to $M_\textrm{UV} \lesssim -12$, with low ionizing photon escape fractions, $f_\textrm{esc} \lesssim 15\%$.

Individual highly magnified stars have been recently discovered at lookback times of more than half the age of the Universe, in lensed galaxies that straddle the critical curves of massive galaxy clusters. Having confirmed their detectability, it is now important to carry out systematic searches for them in order to establish their frequency, and in turn learn about the statistical properties of high-redshift stars and of the granularity of matter in the foreground deflector. Here we report the discovery of a highly magnified star at redshift $z = 0.94$ in a strongly lensed arc behind a Hubble Frontier Field galaxy cluster, MACS J0416.1-2403, discovered as part of a systematic archival search. The bright transient (dubbed "Warhol") was discovered in Hubble Space Telescope data taken on 2014 September 15 and 16. This single image faded over a period of two weeks, and observations taken on 2014 September 1 show that the duration of the microlensing event was at most four weeks in total. The light curve may also exhibit slow changes over a period of years consistent with the level of microlensing expected from stars responsible for the intracluster light (ICL) of the cluster. Optical and infrared observations taken near peak brightness can be fit by a stellar spectrum with moderate host-galaxy extinction. A blue supergiant matches the measured spectral energy distribution near peak, implying a temporary magnification of at least several thousand. While the spectrum of an O-type star would also fit the transient's spectral energy distribution, extremely luminous O-type stars are much less common than blue supergiants. The short timescale of the event and the estimated effective temperature indicate that the lensed source is an extremely magnified star.

We characterize the incidence and intensity of low-level super-massive black hole activity within the Fornax cluster, through X-ray observations of the nuclei of 29 quiescent early-type galaxies. Using the \textit{Chandra X-ray Telescope}, we target 17 galaxies from the HST Fornax Cluster Survey, down to a uniform (3$\sigma$) limiting X-ray luminosity threshold of $5\cdot10^{38}$ergs$^{-1}$, which we combine with deeper, archival observations for an additional 12 galaxies. A nuclear X-ray point-source is detected in 11 out of 29 targets. After accounting for the low mass X-ray binary contamination to the nuclear X-ray signal, the X-ray active fraction is measured at $26.6\% \pm 9.6\%$. The results from this analysis are compared to similar investigations targeting quiescent early types in the Virgo cluster, as well as the field. After correcting for the different mass distributions, the measured Fornax active fraction is less than the field fraction, at more than 3$\sigma$, confirming that the funneling of gas to the nuclear regions of cluster members is inhibited compared to those galaxies in the field. At the same time, we find no statistically significant difference between Fornax and Virgo galaxies, with only marginal evidence for a lower active fraction in Fornax (1 $\sigma$); if real, owing to Fornax's higher galaxy number density, this could indicate that galaxy-galaxy interactions are more effective at gas removal than galaxy-gas effects.

Strong-gravitational lens systems with quadruply-imaged quasars (quads) are unique probes to address several fundamental problems in cosmology and astrophysics. Although they are intrinsically very rare, ongoing and planned wide-field deep-sky surveys are set to discover thousands of such systems in the next decade. It is thus paramount to devise a general framework to model strong-lens systems to cope with this large influx without being limited by expert investigator time. We propose such a general modelling framework (implemented with the publicly available software Lenstronomy) and apply it to uniformly model three-band Hubble Space Telescope Wide Field Camera 3 images of 13 quads. This is the largest uniformly modelled sample of quads to date and paves the way for a variety of studies. To illustrate the scientific content of the sample, we investigate the alignment between the mass and light distribution in the deflectors. The position angles of these distributions are well-aligned, except when there is strong external shear. However, we find no correlation between the ellipticity of the light and mass distributions. We also show that the observed flux-ratios between the images depart significantly from the predictions of simple smooth models. The departures are strongest in the bluest band, consistent with microlensing being the dominant cause in addition to millilensing. Future papers will exploit this rich dataset in combination with ground based spectroscopy and time delays to determine quantities such as the Hubble constant, the free streaming length of dark matter, and the normalization of the initial stellar mass function.

We present a large spectroscopic campaign with Keck/MOSFIRE targeting Lyman-alpha emission (Ly$\alpha$) from intrinsically faint Lyman-break Galaxies (LBGs) behind 12 efficient galaxy cluster lenses. Gravitational lensing allows us to probe the more abundant faint galaxy population to sensitive Ly$\alpha$ equivalent width limits. During the campaign we targeted 70 LBG candidates with MOSFIRE Y-band, selected photometrically to cover Ly$\alpha$ over the range $7<z<8.2$. We detect $S/N>5$ emission lines in 2 of these galaxies and find that they are likely Ly$\alpha$ at $z=7.148\pm0.001$ and $z=7.161\pm0.001$. We present new lens models for 4 of the galaxy clusters, using our previously published lens models for the remaining clusters to determine the magnification factors for the source galaxies. Using a Bayesian framework that employs large scale reionization simulations of the intergalactic medium (IGM) as well as realistic properties of the interstellar medium and circumgalactic medium, we infer the volume-averaged neutral hydrogen fraction, $\overline{x}_{\mathrm{HI}}$, in the IGM during reionization to be $\overline{x}_{\mathrm{HI}}=0.88^{+0.05}_{-0.10}$ at $z=7.6\pm0.6$. Our result is consistent with a late and rapid reionization scenario inferred by Planck.

We present a full data analysis of the pure-parallel Hubble Space Telescope (HST) imaging observations in the Brightest of Reionizing Galaxies Survey (BoRG[z9]) in Cycle 22. The medium-deep exposures with five HST/WFC3IR+UVIS filter bands from 79 independent sightlines (~370 arcmin^2) provide the least biased determination of number density for z>9 bright galaxies against cosmic variance. After a strict two-step selection for candidate galaxies, including dropout color and photometric redshift analyses, and revision of previous BoRG candidates, we identify one source at z~10 and two sources at z~9. The z~10 candidate shows evidence of line-of-sight lens magnification (mu~1.5), yet it appears surprisingly luminous (MUV ~ -22.6\pm0.3 mag), making it one of the brightest candidates at z > 8 known (~ 0.3 mag brighter than the z = 8.68 galaxy EGSY8p7, spectroscopically confirmed by Zitrin and collaborators). For z ~ 9 candidates, we include previous data points at fainter magnitudes and find that the data are well fitted by a Schechter luminosity function with alpha ~ -2.1, MUV ~ -21.5 mag, and log phi ~ -4.5 Mpc^-3mag^-1, for the first time without fixing any parameters. The inferred cosmic star formation rate density is consistent with unaccelerated evolution from lower redshift.

We present a new gravitational lens model of the Hubble Frontier Fields cluster Abell 370 ($z = 0.375$) using imaging and spectroscopy from Hubble Space Telescope and ground-based spectroscopy. We combine constraints from a catalog of 1344 weakly lensed galaxies and 39 multiply-imaged sources comprised of 114 multiple images, including a system of multiply-imaged candidates at $z=7.93 \pm 0.02$, to obtain a best-fit mass distribution using the cluster lens modeling code Strong and Weak Lensing United. As the only analysis of A370 using strong and weak lensing constraints from Hubble Frontier Fields data, our method provides an independent check on assumptions in other methods on the mass distribution. Convergence, shear, and magnification maps are made publicly available through the HFF website. We find that the model we produce is similar to models produced by other groups, with some exceptions due to the differences in lensing code methodology. In an effort to study how our total projected mass distribution traces light, we measure the stellar mass density distribution using Spitzer/Infrared Array Camera imaging. Comparing our total mass density to our stellar mass density in a radius of 0.3 Mpc, we find a mean projected stellar to total mass ratio of $\langle f* \rangle = 0.011 \pm 0.003$ (stat.) using the diet Salpeter initial mass function. This value is in general agreement with independent measurements of $\langle f* \rangle$ in clusters of similar total mass and redshift.

We present models of the H$\beta$-emitting broad-line region (BLR) in seven Seyfert 1 galaxies from the Lick AGN (Active Galactic Nucleus) Monitoring Project 2011 sample, drawing inferences on the BLR structure and dynamics as well as the mass of the central supermassive black hole. We find that the BLR is generally a thick disk, viewed close to face-on, with preferential emission back toward the ionizing source. The dynamics in our sample range from near-circular elliptical orbits to inflowing or outflowing trajectories. We measure black hole masses of $\log_{10}(M_{\rm BH}/M_\odot) = 6.48^{+0.21}_{-0.18}$ for PG 1310$-$108, $7.50^{+0.25}_{-0.18}$ for Mrk 50, $7.46^{+0.15}_{-0.21}$ for Mrk 141, $7.58^{+0.08}_{-0.08}$ for Mrk 279, $7.11^{+0.20}_{-0.17}$ for Mrk 1511, $6.65^{+0.27}_{-0.15}$ for NGC 4593, and $6.94^{+0.14}_{-0.14}$ for Zw 229$-$015. We use these black hole mass measurements along with cross-correlation time lags and line widths to recover the scale factor $f$ used in traditional reverberation mapping measurements. Combining our results with other studies that use this modeling technique, bringing our sample size to 16, we calculate a scale factor that can be used for measuring black hole masses in other reverberation mapping campaigns. When using the root-mean-square (rms) spectrum and using the line dispersion to measure the line width, we find $\log_{10}(f_{{\rm rms},\sigma})_{\rm pred} = 0.57 \pm 0.19$. Finally, we search for correlations between $f$ and other AGN and BLR parameters and find marginal evidence that $f$ is correlated with $M_{\rm BH}$ and the BLR inclination angle, but no significant evidence of a correlation with the AGN luminosity or Eddington ratio.

The Initial Mass Function (IMF) for massive galaxies can be constrained by combining stellar dynamics with strong gravitational lensing. However, this method is limited by degeneracies between the density profile of dark matter and the stellar mass-to-light ratio. In this work we reduce this degeneracy by combining weak lensing together with strong lensing and stellar kinematics. Our analysis is based on two galaxy samples: 45 strong lenses from the SLACS survey and 1,700 massive quiescent galaxies from the SDSS main spectroscopic sample with weak lensing measurements from the Hyper Suprime-Cam survey. We use a Bayesian hierarchical approach to jointly model all three observables. We fit the data with models of varying complexity and show that a model with a radial gradient in the stellar mass-to-light ratio is required to simultaneously describe both galaxy samples. This result is driven by a subset of strong lenses with very steep total density profile, that cannot be fitted by models with no gradient. Our measurements are unable to determine whether $M_*/L$ gradients are due to variations in stellar population parameters at fixed IMF, or to gradients in the IMF itself. The inclusion of $M_*/L$ gradients decreases dramatically the inferred IMF normalisation, compared to previous lensing-based studies, with the exact value depending on the assumed dark matter profile. The main effect of strong lensing selection is to shift the stellar mass distribution towards the high mass end, while the halo mass and stellar IMF distribution at fixed stellar mass are not significantly affected.

The primary goals of the STRong lensing Insights into the Dark Energy Survey (STRIDES) collaboration are to measure the dark energy equation of state parameter and the free streaming length of dark matter. To this aim, STRIDES is discovering strongly lensed quasars in the imaging data of the Dark Energy Survey and following them up to measure time delays, high resolution imaging, and spectroscopy sufficient to construct accurate lens models. In this paper, we first present forecasts for STRIDES. Then, we describe the STRIDES classification scheme, and give an overview of the Fall 2016 follow-up campaign. We continue by detailing the results of two selection methods, the Outlier Selection Technique and a morphological algorithm, and presenting lens models of a system, which could possibly be a lensed quasar in an unusual configuration. We conclude with the summary statistics of the Fall 2016 campaign. Including searches presented in companion papers (Anguita et al.; Ostrovski et al.), STRIDES followed up 117 targets identifying 7 new strongly lensed systems, and 7 nearly identical quasars (NIQs), which could be confirmed as lenses by the detection of the lens galaxy. 76 candidates were rejected and 27 remain otherwise inconclusive, for a success rate in the range 6-35\%. This rate is comparable to that of previous searches like SQLS even though the parent dataset of STRIDES is purely photometric and our selection of candidates cannot rely on spectroscopic information.

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

The free streaming length of dark matter particles determines the abundance of structure on sub-galactic scales. We present a statistical technique, amendable to any parameterization of subhalo density profile and mass function, to probe dark matter on these scales with quadrupole image lenses. We consider a warm dark matter particle with a mass function characterized by a normalization and free streaming scale $m_{\rm{hm}}$. We forecast bounds on dark matter warmth for 120-180 lenses, attainable with future surveys, at typical lens (source) redshifts of 0.5 (1.5) in early-type galaxies with velocity dispersions of 220-270 km/sec. We demonstrate that limits on $m_{\rm{hm}}$ deteriorate rapidly with increasing uncertainty in image fluxes, underscoring the importance of precise measurements and accurate lens models. For our forecasts, we assume the deflectors in the lens sample do not exhibit complex morphologies, so we neglect systematic errors in their modeling. Omitting the additional signal from line of sight halos, our constraints underestimate the true power of the method. Assuming cold dark matter, for a low normalization, corresponding the destruction of all subhalos within the host scale radius, we forecast $2\sigma$ bounds on $m_{\rm{hm}}$ (thermal relic mass) of $10^{7.5} \ (5.0)$, $10^{8} \ (3.6)$, and $10^{8.5} \ (2.7) \ M_{\odot} \left(\rm{keV}\right)$ for flux errors of $2\%$, $4\%$, and $8\%$. With a higher normalization, these constraints improve to $10^{7.2} \ (6.6)$, $10^{7.5} \ (5.3) $, and $10^{7.8} \ (4.3) \ M_{\odot} \left(\rm{keV}\right)$ with 120 systems. We are also able to measure the normalization of the mass function, which has implications for baryonic feedback models and tidal stripping.

For a sample of $\sim$80 local ($0.02 \leq z \leq 0.1$) Seyfert-1 galaxies with high-quality long-slit Keck spectra and spatially-resolved stellar-velocity dispersion ($\sigma_{\star}$) measurements, we study the profile of the [OIII]$\lambda$5007A emission line to test the validity of using its width as a surrogate for $\sigma_{\star}$. Such an approach has often been used in the literature, since it is difficult to measure $\sigma_{\star}$ for type-1 active galactic nuclei (AGNs) due to the AGN continuum outshining the stellar-absorption lines. Fitting the [OIII] line with a single Gaussian or Gauss-Hermite polynomials overestimates $\sigma_{\star}$ by 50-100%. When line asymmetries from non-gravitational gas motion are excluded in a double Gaussian fit, the average ratio between the core [OIII] width ($\sigma_{\rm {[OIII],D}}$) and $\sigma_{\star}$ is $\sim$1, but with individual data points off by up to a factor of two. The resulting black-hole-mass-$\sigma_{\rm {[OIII],D}}$ relation scatters around that of quiescent galaxies and reverberation-mapped AGNs. However, a direct comparison between $\sigma_{\star}$ and $\sigma_{\rm {[OIII],D}}$ shows no close correlation, only that both quantities have the same range, average and standard deviation, probably because they feel the same gravitational potential. The large scatter is likely due to the fact that line profiles are a luminosity-weighted average, dependent on the light distribution and underlying kinematic field. Within the range probed by our sample (80-260 km s$^{-1}$), our results strongly caution against the use of [OIII] width as a surrogate for $\sigma_{\star}$ on an individual basis. Even though our sample consists of radio-quiet AGNs, FIRST radio-detected objects have, on average, a $\sim$10% larger [OIII] core width.

We present followup imaging of two bright (L > L*) galaxy candidates at z > 8 from the Brightest of Reionizing Galaxies (BoRG) survey with the F098M filter on HST/WFC3. The F098M filter provides an additional constraint on the flux blueward of the spectral break, and the observations are designed to discriminate between low- and high-z photometric redshift solutions for these galaxies. Our results confirm one galaxy, BoRG 0116+1425 747, as a highly probable z ~ 8 source, but reveal that BoRG 0116+1425 630 - previously the brightest known z > 8 candidate (mAB = 24.5) - is likely to be a z ~ 2 interloper. As this source was substantially brighter than any other z > 8 candidate, removing it from the sample has a significant impact on the derived UV luminosity function in this epoch. We show that while previous BoRG results favored a shallow power-law decline in the bright end of the luminosity function prior to reionization, there is now no evidence for departure from a Schechter function form and therefore no evidence for a difference in galaxy formation processes before and after reionization.

We present a gravitational lensing model of MACS J1149.5+2223 using ultra-deep Hubble Frontier Fields imaging data and spectroscopic redshifts from HST grism and VLT/MUSE spectroscopic data. We create total mass maps using 38 multiple images (13 sources) and 608 weak lensing galaxies, as well as 100 multiple images of 31 star-forming regions in the galaxy that hosts Supernova Refsdal. We find good agreement with a range of recent models within the HST field of view. We present a map of the ratio of projected stellar mass to total mass ($f_{\star}$), and find that the stellar mass fraction for this cluster peaks on the primary BCG. Averaging within a radius of 0.3 Mpc, we obtain a value of $\langle f_{\star} \rangle = 0.012^{+0.004}_{-0.003}$, consistent with other recent results for this ratio in cluster environments, though with a large global error (up to $\delta f_{\star} = 0.005$) primarily due to the choice of an IMF. We compare values of $f_{\star}$ and measures of star formation efficiency for this cluster to other Hubble Frontier Fields clusters studied in the literature, finding that MACS1149 has a higher stellar mass fraction than these other clusters, but a star formation efficiency typical of massive clusters.

We report upon the follow up of 34 candidate lensed quasars found in the Dark Energy Survey using NTT-EFOSC, Magellan-IMACS, KECK-ESI and SOAR-SAMI. These candidates were selected by a combination of double component fitting, morphological assessment and color analysis. Most systems followed up are indeed composed of at least one quasar image and 13 with two or more quasar images: two lenses, four projected binaries and seven Nearly Identical Quasar Pairs (NIQs). The two systems confirmed as genuine gravitationally lensed quasars are one quadruple at $z_s=1.713$ and one double at $z_s=1.515$. Lens modeling of these two systems reveals that both systems require very little contribution from the environment to reproduce the image configuration. Nevertheless, small flux anomalies can be observed in one of the images of the quad. Further observations of 9 inconclusive systems (including 7 NIQs) will allow to confirm (or not) their gravitational lens nature.

We present the single-epoch black hole mass (M$_{\rm BH}$) calibrations based on the rest-frame UV and optical measurements of Mg II 2798\AA\ and H$\beta$ 4861\AA\ lines and AGN continuum, using a sample of 52 moderate-luminosity AGNs at z$\sim$0.4 and z$\sim$0.6 with high-quality Keck spectra. We combine this sample with a large number of luminous AGNs from the Sloan Digital Sky Survey to increase the dynamic range for a better comparison of UV and optical velocity and luminosity measurements. With respect to the reference M$_{\rm BH}$ based on the line dispersion of H$\beta$ and continuum luminosity at 5100\AA, we calibrate the UV and optical mass estimators, by determining the best-fit values of the coefficients in the mass equation. By investigating whether the UV estimators show systematic trend with Eddington ratio, FWHM of H$\beta$, the Fe II strength, and the UV/optical slope, we find no significant bias except for the slope. By fitting the systematic difference of Mg II-based and H$\beta$-based masses with the L$_{3000}$/L$_{5100}$ ratio, we provide a correction term as a function of the spectral index as $\Delta$C = 0.24 (1+$\alpha_{\lambda}$) + 0.17, which can be added to the Mg II-based mass estimators if the spectral slope can be well determined. The derived UV mass estimators typically show $>$$\sim$0.2 dex intrinsic scatter with respect to H$\beta$-based M$_{\rm BH}$, suggesting that the UV-based mass has an additional uncertainty of $\sim$0.2 dex, even if high quality rest-frame UV spectra are available.

VANDELS is a uniquely-deep spectroscopic survey of high-redshift galaxies with the VIMOS spectrograph on ESO's Very Large Telescope (VLT). The survey has obtained ultra-deep optical (0.48 < lambda < 1.0 micron) spectroscopy of ~2100 galaxies within the redshift interval 1.0 < z < 7.0, over a total area of ~0.2 sq. degrees centred on the CANDELS UDS and CDFS fields. Based on accurate photometric redshift pre-selection, 85% of the galaxies targeted by VANDELS were selected to be at z>=3. Exploiting the red sensitivity of the refurbished VIMOS spectrograph, the fundamental aim of the survey is to provide the high signal-to-noise ratio spectra necessary to measure key physical properties such as stellar population ages, masses, metallicities and outflow velocities from detailed absorption-line studies. Using integration times calculated to produce an approximately constant signal-to-noise ratio (20 < t_int < 80 hours), the VANDELS survey targeted: a) bright star-forming galaxies at 2.4 < z < 5.5, b) massive quiescent galaxies at 1.0 < z < 2.5, c) fainter star-forming galaxies at 3.0 < z < 7.0 and d) X-ray/Spitzer-selected active galactic nuclei and Herschel-detected galaxies. By targeting two extragalactic survey fields with superb multi-wavelength imaging data, VANDELS will produce a unique legacy data set for exploring the physics underpinning high-redshift galaxy evolution. In this paper we provide an overview of the VANDELS survey designed to support the science exploitation of the first ESO public data release, focusing on the scientific motivation, survey design and target selection.

Realizing Refsdal's original idea from 1964, we present estimates of the Hubble constant that are complementary to and potentially competitive with those of other cosmological probes. We use the observed positions of 89 multiple images, with extensive spectroscopic information, from 28 background sources and the measured time delays between the images S1-S4 and SX of supernova "Refsdal" ($z = 1.489$), which were obtained thanks to Hubble Space Telescope (HST) deep imaging and Multi Unit Spectroscopic Explorer (MUSE) data. We extend the strong lensing modeling of the Hubble Frontier Fields (HFF) galaxy cluster MACS J1149.5$+$2223 ($z = 0.542$), published by Grillo et al. (2016), and explore different $\Lambda$CDM models. Taking advantage of the lensing information associated to the presence of very close pairs of multiple images at various redshifts and to the extended surface brightness distribution of the SN Refsdal host, we can reconstruct the total mass density profile of the cluster very precisely. The combined dependence of the multiple image positions and time delays on the cosmological parameters allows us to infer the values of $H_{0}$ and $\Omega_{\rm m}$ with relative (1$\sigma$) statistical errors of, respectively, 6% (7%) and 31% (26%) in flat (general) cosmological models, assuming a conservative 3% uncertainty on the final time delay of image SX and, remarkably, no priors from other cosmological experiments. Our best estimate of $H_{0}$, based on the model described in this work, will be presented when the final time-delay measurement becomes available. Our results show that it is possible to utilize time delays in lens galaxy clusters as an important alternative tool for measuring the expansion rate and the geometry of the Universe.

We present time-delay estimates for the quadruply imaged quasar PG 1115+080. Our resuls are based on almost daily observations for seven months at the ESO MPIA 2.2m telescope at La Silla Observatory, reaching a signal-to-noise ratio of about 1000 per quasar image. In addition, we re-analyse existing light curves from the literature that we complete with an additional three seasons of monitoring with the Mercator telescope at La Palma Observatory. When exploring the possible source of bias we consider the so-called microlensing time delay, a potential source of systematic error so far never directly accounted for in previous time-delay publications. In fifteen years of data on PG 1115+080, we find no strong evidence of microlensing time delay. Therefore not accounting for this effect, our time-delay estimates on the individual data sets are in good agreement with each other and with the literature. Combining the data sets, we obtain the most precise time-delay estimates to date on PG 1115+080, with Dt(AB) = 8.3+1.5-1.6 days (18.7% precision), Dt(AC) = 9.9+1.1-1.1 days (11.1%) and Dt(BC) = 18.8+1.6-1.6 days (8.5%). Turning these time delays into cosmological constraints is done in a companion paper that makes use of ground-based Adaptive Optics (AO) with the Keck telescope.

Strongly gravitational lensed quasars can be used to measure the so-called time-delay distance $D_{\Delta t}$, and thus the Hubble constant $H_0$ and other cosmological parameters. Stellar kinematics of the deflector galaxy play an essential role in this measurement by: (i) helping break the mass-sheet degeneracy; (ii) determining in principle the angular diameter distance $D_{\rm d}$ to the deflector and thus further improving the cosmological constraints. In this paper we simulate observations of lensed quasars with integral field spectrographs and show that spatially resolved kinematics of the deflector enable further progress by helping break the mass-anisotropy degeneracy. Furthermore, we use our simulations to obtain realistic error estimates with current/upcoming instruments like OSIRIS on Keck and NIRSPEC on the James Webb Space Telescope for both distances (typically $\sim6$ per cent on $D_{\Delta t}$ and $\sim10$ per cent on $D_{\rm d}$). We use the error estimates to compute cosmological forecasts for the sample of nine lenses that currently have well measured time delays and deep Hubble Space Telescope images and for a sample of 40 lenses that is projected to be available in a few years through follow-up of candidates found in ongoing wide field surveys. We find that $H_0$ can be measured with 2 per cent (1 per cent) precision from nine (40) lenses in a flat $\Lambda$cold dark matter cosmology. We study several other cosmological models beyond the flat $\Lambda$cold dark matter model and find that time-delay lenses with spatially resolved kinematics can greatly improve the precision of the cosmological parameters measured by cosmic microwave background data.

Recent detections of Lyman alpha (Ly$\alpha$) emission from $z>7.5$ galaxies were somewhat unexpected given a dearth of previous non-detections in this era when the intergalactic medium (IGM) is still highly neutral. But these detections were from UV bright galaxies, which preferentially live in overdensities which reionize early, and have significantly Doppler-shifted Ly$\alpha$ line profiles emerging from their interstellar media (ISM), making them less affected by the global IGM state. Using a combination of reionization simulations and empirical ISM models we show, as a result of these two effects, UV bright galaxies in overdensities have $>2\times$ higher transmission through the $z\sim7$ IGM than typical field galaxies, and this boosted transmission is enhanced as the neutral fraction increases. The boosted transmission is not sufficient to explain the observed high Ly$\alpha$ fraction of $M_\mathrm{UV} \lesssim -22$ galaxies (Stark et al. 2017), suggesting Ly$\alpha$ emitted by these galaxies must be stronger than expected due to enhanced production and/or selection effects. Despite the bias of UV bright galaxies to reside in overdensities we show Ly$\alpha$ observations of such galaxies can accurately measure the global neutral hydrogen fraction, particularly when Ly$\alpha$ from UV faint galaxies is extinguished, making them ideal candidates for spectroscopic follow-up into the cosmic Dark Ages.

Galaxy-cluster gravitational lenses can magnify background galaxies by a total factor of up to ~50. Here we report an image of an individual star at redshift z=1.49 (dubbed "MACS J1149 Lensed Star 1 (LS1)") magnified by >2000. A separate image, detected briefly 0.26 arcseconds from LS1, is likely a counterimage of the first star demagnified for multiple years by a >~3 solar-mass object in the cluster. For reasonable assumptions about the lensing system, microlensing fluctuations in the stars' light curves can yield evidence about the mass function of intracluster stars and compact objects, including binary fractions and specific stellar evolution and supernova models. Dark-matter subhalos or massive compact objects may help to account for the two images' long-term brightness ratio.

This paper describes the observations and the first data release (DR1) of the ESO public spectroscopic survey "VANDELS, a deep VIMOS survey of the CANDELS CDFS and UDS fields". VANDELS' main targets are star-forming galaxies at 2.4<z<5.5 and massive passive galaxies at 1<z<2.5. By adopting a strategy of ultra-long exposure times, from 20 to 80 hours per source, VANDELS is designed to be the deepest ever spectroscopic survey of the high-redshift Universe. Exploiting the red sensitivity of the VIMOS spectrograph, the survey has obtained ultra-deep spectra covering the wavelength 4800-10000 A with sufficient signal-to-noise to investigate the astrophysics of high-redshift galaxy evolution via detailed absorption line studies. The VANDELS-DR1 is the release of all spectra obtained during the first season of observations and includes data for galaxies for which the total (or half of the total) scheduled integration time was completed. The release contains 879 individual objects with a measured redshift and includes fully wavelength and flux-calibrated 1D spectra, the associated error spectra, sky spectra and wavelength-calibrated 2D spectra. We also provide a catalog with the essential galaxy parameters, including spectroscopic redshifts and redshift quality flags. In this paper we present the survey layout and observations, the data reduction and redshift measurement procedure and the general properties of the VANDELS-DR1 sample. We also discuss the spectroscopic redshift distribution, the accuracy of the photometric redshifts and we provide some examples of data products. All VANDELS-DR1 data are publicly available and can be retrieved from the ESO archive. Two further data releases are foreseen in the next 2 years with a final release scheduled for June 2020 which will include improved re-reduction of the entire spectroscopic data set. (abridged)

Constraining the sub-galactic matter-power spectrum on 1-10 kpc scales would make it possible to distinguish between the concordance $\Lambda$CDM model and various alternative dark-matter models due to the significantly different levels of predicted mass structure. Here, we demonstrate a novel approach to observationally constrain the population of overall law-mass density fluctuations in the inner regions of massive elliptical lens galaxies, based on the power spectrum of the associated surface-brightness perturbations observable in highly magnified galaxy-scale Einstein rings and gravitational arcs. The application of our method to the SLACS lens system SDSS J0252+0039 results in the following limits (at the 99 per cent confidence level) on the dimensionless convergence-power spectrum (and the associated standard deviation in aperture mass): $\Delta^{2}_{\delta\kappa}<1$ ($\sigma_{AM}< 0.8 \times 10^8 M_\odot$) on 0.5-kpc scale, $\Delta^{2}_{\delta\kappa}<0.1$ ($\sigma_{AM}< 1 \times 10^8 M_\odot$) on 1-kpc scale and $\Delta^{2}_{\delta\kappa}<0.01$ ($\sigma_{AM}< 3 \times 10^8 M_\odot$) on 3-kpc scale. The estimated effect of CDM sub-haloes lies considerably below these first observational upper-limit constraints on the level of inhomogeneities in the projected total mass distribution of galactic haloes. Future analysis for a larger sample of galaxy-galaxy strong lens systems will narrow down these constraints and rule out all cosmological models predicting a significantly larger level of clumpiness on these critical sub-galactic scales.

The Seyfert 1 galaxy Arp 151 was monitored as part of three reverberation mapping campaigns spanning $2008-2015$. We present modeling of these velocity-resolved reverberation mapping datasets using a geometric and dynamical model for the broad line region (BLR). By modeling each of the three datasets independently, we infer the evolution of the BLR structure in Arp 151 over a total of seven years and constrain the systematic uncertainties in non-varying parameters such as the black hole mass. We find that the BLR geometry of a thick disk viewed close to face-on is stable over this time, although the size of the BLR grows by a factor of $\sim 2$. The dynamics of the BLR are dominated by inflow and the inferred black hole mass is consistent for the three datasets, despite the increase in BLR size. Combining the inference for the three datasets yields a black hole mass and statistical uncertainty of $\log_{10}($M$_{\rm BH}/\rm{M}_{\odot})=6.82^{+0.09}_{-0.09}$ with a standard deviation in individual measurements of 0.13 dex.

Local massive early-type galaxies are believed to have completed most of their star formation $\sim10$Gyr ago and evolved without having substantial star formation since. If so, their progenitors should have roughly solar stellar metallicities ($Z_*$), comparable to their values today. We report the discovery of two lensed massive ($\log M_*/M_\odot\sim11$), $z\sim2.2$ dead galaxies, that appear markedly metal deficient given this scenario. Using 17-band $HST$+$K_{s}$+$Spitzer$ photometry and deep $HST$ grism spectra from the GLASS and SN Refsdal follow-up campaigns covering features near $\lambda_{\rm rest}\sim4000$\AA, we find these systems to be dominated by A-type stars with $\log Z_*/Z_\odot=-0.40\pm0.02$ and $-0.49\pm0.03$ ($30$-$40\%$ solar) under standard assumptions. The second system's lower metallicity is robust to isochrone changes, though this choice can drive the first system's from $\log Z_*/Z_\odot=-0.6$ to 0.1. If these two galaxies are representative of larger samples, this finding suggests that evolutionary paths other than dry minor-merging are required for these massive galaxies. Future analyses with direct metallicity measurements-e.g., by the $James\ Webb\ Space\ Telescope$-will provide critical insight into the nature of such phenomena.

We present the stellar mass-stellar metallicity relationship (MZR) in the Cl0024+1654 galaxy cluster at z~0.4 using full spectrum stellar population synthesis modeling of individual quiescent galaxies. The lower limit of our stellar mass range is $M_*=10^{9.7}M_\odot$, the lowest galaxy mass at which individual stellar metallicity has been measured beyond the local universe. We report a detection of an evolution of the stellar MZR with observed redshift at $0.037\pm0.007$ dex per Gyr, consistent with the predictions from hydrodynamical simulations. Additionally, we find that the evolution of the stellar MZR with observed redshift can be explained by an evolution of the stellar MZR with their formation time, i.e., when the single stellar population (SSP)-equivalent ages of galaxies are taken into account. This behavior is consistent with stars forming out of gas that also has an MZR with a normalization that decreases with redshift. Lastly, we find that over the observed mass range, the MZR can be described by a linear function with a shallow slope, ($[Fe/H] \propto (0.16 \pm 0.03) \log M_*$). The slope suggests that galaxy feedback, in terms of mass-loading factor, might be mass-independent over the observed mass and redshift range.

We present a new flexible Bayesian framework for directly inferring the fraction of neutral hydrogen in the intergalactic medium (IGM) during the Epoch of Reionization (EoR, z~6-10) from detections and non-detections of Lyman Alpha (Ly$\alpha$) emission from Lyman break galaxies (LBGs). Our framework combines sophisticated reionization simulations with empirical models of the interstellar medium (ISM) radiative transfer effects on Ly$\alpha$. We assert that the Ly$\alpha$ line profile emerging from the ISM has an important impact on the resulting transmission of photons through the IGM, and that these line profiles depend on galaxy properties. We model this effect by considering the peak velocity offset of Ly$\alpha$ lines from host galaxies' systemic redshifts, which are empirically correlated with UV luminosity and redshift (or halo mass at fixed redshift). We use our framework on the sample of LBGs presented in Pentericci et al. (2014) and infer a global neutral fraction at z~7 of $\overline{x}_\mathrm{HI} = 0.59_{-0.15}^{+0.11}$, consistent with other robust probes of the EoR and confirming reionization is on-going ~700 Myr after the Big Bang. We show that using the full distribution of Ly$\alpha$ equivalent width detections and upper limits from LBGs places tighter constraints on the evolving IGM than the standard Ly$\alpha$ emitter fraction, and that larger samples are within reach of deep spectroscopic surveys of gravitationally lensed fields and JWST NIRSpec.

We use Multi Unit Spectroscopic Explorer (MUSE) observations of the galaxy cluster MACS J1149.5+2223 to explore the kinematics of the grand-design spiral galaxy Sp1149 hosting the SN Refsdal. Sp1149 lies at $z\simeq1.49$, has a stellar mass $M_*\simeq5\times10^9 \, \mathrm{M_\odot}$, a star-formation rate $\mathrm{SFR} \simeq1-6 \, \mathrm{M_\odot/yr}$ and represents a likely progenitor of a Milky-Way-like galaxy. All the four multiple images of Sp1149 in our data show strong OII-line emissions pointing to a clear rotation pattern. We take advantage of the gravitational lensing magnification effect ($\simeq 4 \times$) on the OII emission of the least distorted image to fit 3D kinematic models to the MUSE data-cube and derive the rotation curve and the velocity dispersion profile of Sp1149. We find that the rotation curve steeply rises, peaks at $R\simeq1$ kpc and then (initially) declines and flattens to an average $V_\mathrm{flat} = 128^{+29}_{-19}$ km/s. The shape of the rotation curve is well determined but the actual value of $V_\mathrm{flat}$ is quite uncertain because of the nearly face-on configuration of the galaxy. The intrinsic velocity dispersion due to gas turbulence is almost constant across the entire disc with an average of $27\pm5$ km/s. This value is consistent with $z=0$ measurements in the ionized gas component and a factor 2-4 lower than other estimates in different galaxies at similar redshifts. The average stellar-to-total mass fraction is of the order of one fifth. Our kinematic analysis returns the picture of a regular star-forming, mildly turbulent, rotation-dominated ($V / \sigma\simeq5$) spiral galaxy in a 4 Gyr old Universe.

We present deep spectroscopic observations of a Lyman-break galaxy candidate (hereafter MACS1149-JD) at $z\sim9.5$ with the $\textit{Hubble}$ Space Telescope ($\textit{HST}$) WFC3/IR grisms. The grism observations were taken at 4 distinct position angles, totaling 34 orbits with the G141 grism, although only 19 of the orbits are relatively uncontaminated along the trace of MACS1149-JD. We fit a 3-parameter ($z$, F160W mag, and Ly$\alpha$ equivalent width) Lyman-break galaxy template to the three least contaminated grism position angles using an MCMC approach. The grism data alone are best fit with a redshift of $z_{\mathrm{grism}}=9.53^{+0.39}_{-0.60}$ ($68\%$ confidence), in good agreement with our photometric estimate of $z_{\mathrm{phot}}=9.51^{+0.06}_{-0.12}$ ($68\%$ confidence). Our analysis rules out Lyman-alpha emission from MACS1149-JD above a $3\sigma$ equivalent width of 21 \AA{}, consistent with a highly neutral IGM. We explore a scenario where the red $\textit{Spitzer}$/IRAC $[3.6] - [4.5]$ color of the galaxy previously pointed out in the literature is due to strong rest-frame optical emission lines from a very young stellar population rather than a 4000 \AA{} break. We find that while this can provide an explanation for the observed IRAC color, it requires a lower redshift ($z\lesssim9.1$), which is less preferred by the $\textit{HST}$ imaging data. The grism data are consistent with both scenarios, indicating that the red IRAC color can still be explained by a 4000 \AA{} break, characteristic of a relatively evolved stellar population. In this interpretation, the photometry indicate that a $340^{+29}_{-35}$ Myr stellar population is already present in this galaxy only $\sim500~\mathrm{Myr}$ after the Big Bang.

We present optical continuum lags for two Seyfert 1 galaxies, MCG+08-11-011 and NGC 2617, using monitoring data from a reverberation mapping campaign carried out in 2014. Our light curves span the ugriz filters over four months, with median cadences of 1.0 and 0.6 days for MCG+08-11-011 and NGC\,2617, respectively, combined with roughly daily X-ray and near-UV data from Swift for NGC 2617. We find lags consistent with geometrically thin accretion-disk models that predict a lag-wavelength relation of $\tau \propto \lambda^{4/3}$. However, the observed lags are larger than predictions based on standard thin-disk theory by factors of 3.3 for MCG+08-11-011 and 2.3 for NGC\,2617. These differences can be explained if the mass accretion rates are larger than inferred from the optical luminosity by a factor of 4.3 in MCG+08-11-011 and a factor of 1.3 in NGC\,2617, although uncertainty in the SMBH masses determines the significance of this result. While the X-ray variability in NGC\,2617 precedes the UV/optical variability, the long 2.6 day lag is problematic for coronal reprocessing models.

We present the full sample of 118 galaxy-scale strong-lens candidates in the Sloan Lens ACS (SLACS) Survey for the Masses (S4TM) Survey, which are spectroscopically selected from the final data release of the Sloan Digital Sky Survey. Follow-up Hubble Space Telescope (HST) imaging observations confirm that 40 candidates are definite strong lenses with multiple lensed images. The foreground lens galaxies are found to be early-type galaxies (ETGs) at redshifts 0.06 to 0.44, and background sources are emission-line galaxies at redshifts 0.22 to 1.29. As an extension of the SLACS Survey, the S4TM Survey is the first attempt to preferentially search for strong-lens systems with relatively lower lens masses than those in the pre-existing strong-lens samples. By fitting HST data with a singular isothermal ellipsoid model, we find total projected mass within the Einstein radius of the S4TM strong-lens sample ranges from $3 \times10^{10} M_{\odot}$ to $2 \times10^{11} M_{\odot}$. In [Shu15], we have derived the total stellar mass of the S4TM lenses to be $5 \times10^{10} M_{\odot}$ to $1 \times10^{12} M_{\odot}$. Both total enclosed mass and stellar mass of the S4TM lenses are on average almost a factor of 2 smaller than those of the SLACS lenses, which also represent typical mass scales of the current strong-lens samples. The extended mass coverage provided by the S4TM sample can enable a direct test, with the aid of strong lensing, for transitions in scaling relations, kinematic properties, mass structure, and dark-matter content trends of ETGs at intermediate-mass scales as noted in previous studies.

We present the first results from an optical reverberation mapping campaign executed in 2014, targeting the active galactic nuclei (AGN) MCG+08-11-011, NGC 2617, NGC 4051, 3C 382, and Mrk 374. Our targets have diverse and interesting observational properties, including a "changing look" AGN and a broad-line radio galaxy. Based on continuum-H$\beta$ lags, we measure black hole masses for all five targets. We also obtain H$\gamma$ and He{\sc ii}\,$\lambda 4686$ lags for all objects except 3C 382. The He{\sc ii}\,$\lambda 4686$ lags indicate radial stratification of the BLR, and the masses derived from different emission lines are in general agreement. The relative responsivities of these lines are also in qualitative agreement with photoionization models. These spectra have extremely high signal-to-noise ratios (100--300 per pixel) and there are excellent prospects for obtaining velocity-resolved reverberation signatures.

We present a weak gravitational lensing measurement of the external convergence along the line of sight to the quadruply lensed quasar HE$\,$0435$-$1223. Using deep r-band images from Subaru-Suprime-Cam we observe galaxies down to a 3$\sigma$ limiting magnitude of $\sim 26$ mags resulting in a source galaxy density of 14 galaxies / arcmin$^2$ after redshift-based cuts. Using an inpainting technique and Multi-Scale Entropy filtering algorithm, we find that the region in close proximity to the lens has an estimated external convergence of $\kappa=-0.012^{+0.020}_{-0.013}$ and is hence marginally under-dense. We also rule out the presence of any halo with a mass greater than $M_{\rm vir}=1.6\times10^{14}h^{-1}M_\odot$ (68$\%$ confidence limit). Our results, consistent with previous studies of this lens, confirm that the intervening mass along the line of sight to HE$\,$0435$-$1223 does not affect significantly the cosmological results inferred from the time delay measurements of that specific object.

We report on discovery results from a quasar lens search in the ATLAS public footprint, extending quasar lens searches to a regime without $u-$band or fiber-spectroscopic information, using a combination of data mining techniques on multi-band catalog magnitudes and image-cutout modelling. Spectroscopic follow-up campaigns, conducted at the 2.6m Nordic Optical Telescope (La Palma) and 3.6m New Technology Telescope (La Silla) in 2016, yielded seven pairs of quasars exhibiting the same lines at the same redshift and monotonic flux-ratios with wavelength (hereafter NIQs, Nearly Identical Quasar pairs). The quasar redshifts range between $\approx1.2$ and $\approx 2.7;$ contaminants are typically pairs of bright blue stars, quasar-star alignments along the line of sight, and narrow-line galaxies at $0.3<z<0.7.$ Magellan data of A0140-1152 (01$^h$40$^m$03.0$^s$-11$^d$52$^m$19.0$^s$, $z_{s}=1.807$) confirm it as a lens with deflector at $z_{l}=0.277$ and Einstein radius $\theta_{\rm E}=(0.73\pm0.02)^\ase$. We show the use of spatial resolution from the Gaia mission to select lenses and list additional systems from a WISE-Gaia-ATLAS search, yielding three additional lenses (02$^h$35$^m$27.4$^s$-24$^d$33$^m$13.2$^s$, 02$^h$59$^m$33.$^s$-23$^d$38$^m$01.8$^s$, 01$^h$46$^m$32.9$^s$-11$^d$33$^m$39.0$^s$). The overall sample consists of 11 lenses/NIQs, plus three lenses known before 2016, over the ATLAS-DR3 footprint ($\approx3500$~deg$^2$). Finally, we discuss future prospects for objective classification of pair/NIQ/contaminant spectra.

We report the discovery, spectroscopic confirmation and first lens models of the first two, strongly lensed quasars from a combined search in WISE and Gaia over the DES footprint. The four-image lensWGD2038-4008 (r.a.=20:38:02.65, dec.=-40:08:14.64) has source- and lens-redshifts $z_{s}=0.777 \pm 0.001$ and $z_l = 0.230 \pm 0.002$ respectively. Its deflector has effective radius $R_{\rm eff} \approx 3.4^{\prime\prime}$, stellar mass $\log(M_{\star}/M_{\odot}) = 11.64^{+0.20}_{-0.43}$, and shows extended isophotal shape variation. Simple lens models yield Einstein radii $R_{\rm E}=(1.30\pm0.04)^{\prime\prime},$ axis ratio $q=0.75\pm0.1$ (compatible with that of the starlight) and considerable shear-ellipticity degeneracies. The two-image lensWGD2021-4115 (r.a.=20:21:39.45, dec.=--41:15:57.11) has $z_{s}=1.390\pm0.001$ and $z_l = 0.335 \pm 0.002$, and Einstein radius $R_{\rm E} = (1.1\pm0.1)^{\prime\prime},$ but higher-resolution imaging is needed to accurately separate the deflector and faint quasar image. We also show high-rank candidate doubles selected this way, some of which have been independently identified with different techniques, and discuss a DES+WISE quasar multiplet selection.

We present time-delay measurements for the new quadruply imaged quasar DES J0408-5354, the first quadruply imaged quasar found in the Dark Energy Survey (DES). Our result is made possible by implementing a new observational strategy using almost daily observations with the MPIA 2.2m telescope at La Silla observatory and deep exposures reaching a signal-to-noise ratio of about 1000 per quasar image. This data quality allows us to catch small photometric variations (a few mmag rms) of the quasar, acting on temporal scales much shorter than microlensing, hence making the time delay measurement very robust against microlensing. In only 7 months we measure very accurately one of the time delays in DES J0408-5354: Dt(AB) = -112.1 +- 2.1 days (1.8%) using only the MPIA 2.2m data. In combination with data taken with the 1.2m Euler Swiss telescope, we also measure two delays involving the D component of the system Dt(AD) = -155.5 +- 12.8 days (8.2%) and Dt(BD) = -42.4 +- 17.6 days (41%), where all the error bars include systematics. Turning these time delays into cosmological constraints will require deep HST imaging or ground-based Adaptive Optics (AO), and information on the velocity field of the lensing galaxy.

Modern data empower observers to describe galaxies as the spatially and biographically complex objects they are. We illustrate this through case studies of four, $z\sim1.3$ systems based on deep, spatially resolved, 17-band + G102 + G141 Hubble Space Telescope grism spectrophotometry. Using full spectrum rest-UV/-optical continuum fitting, we characterize these galaxies' observed $\sim$kpc-scale structures and star formation rates (SFRs) and reconstruct their history over the age of the universe. The sample's diversity---passive to vigorously starforming; stellar masses $\log M_*/M_\odot=10.5$ to $11.2$---enables us to draw spatio-temporal inferences relevant to key areas of parameter space (Milky Way- to super-Andromeda-mass progenitors). Specifically, we find signs that bulge mass-fractions ($B/T$) and SF history shapes/spatial uniformity are linked, such that higher $B/T$s correlate with "inside-out growth" and central specific SFRs that peaked above the global average for all starforming galaxies at that epoch. Conversely, the system with the lowest $B/T$ had a flat, spatially uniform SFH with normal peak activity. Both findings are consistent with models positing a feedback-driven connection between bulge formation and the switch from rising to falling SFRs ("quenching"). While sample size forces this conclusion to remain tentative, this work provides a proof-of-concept for future efforts to refine or refute it: JWST, WFIRST, and the 30-m class telescopes will routinely produce data amenable to this and more sophisticated analyses. These samples---spanning representative mass, redshift, SFR, and environmental regimes---will be ripe for converting into thousands of sub-galactic-scale empirical windows on what individual systems actually looked like in the past, ushering in a new dialog between observation and theory.

We present inferences on the geometry and kinematics of the broad-Hbeta line-emitting region in four active galactic nuclei monitored as a part of the fall 2010 reverberation mapping campaign at MDM Observatory led by the Ohio State University. From modeling the continuum variability and response in emission-line profile changes as a function of time, we infer the geometry of the Hbeta- emitting broad line regions to be thick disks that are close to face-on to the observer with kinematics that are well-described by either elliptical orbits or inflowing gas. We measure the black hole mass to be log (MBH) = 7.25 (+/-0.10) for Mrk 335, 7.86 (+0.20, -0.17) for Mrk 1501, 7.84 (+0.14, -0.19) for 3C 120, and 6.92 (+0.24, -0.23) for PG 2130+099. These black hole mass measurements are not based on a particular assumed value of the virial scale factor f, allowing us to compute individual f factors for each target. Our results nearly double the number of targets that have been modeled in this manner, and investigate the properties of a more diverse sample by including previously modeled objects. We measure an average scale factor f in the entire sample to be log10(f) = 0.54 +/- 0.17 when the line dispersion is used to characterize the line width, which is consistent with values derived using the normalization of the MBH-sigma relation. We find that the scale factor f for individual targets is likely correlated with the black hole mass, inclination angle, and opening angle of the broad line region but we do not find any correlation with the luminosity.

Strongly lensed active galactic nuclei (AGN) provide a unique opportunity to make progress in the study of the evolution of the correlation between the mass of supermassive black holes ($\mathcal M_{BH}$) and their host galaxy luminosity ($L_{host}$). We demonstrate the power of lensing by analyzing two systems for which state-of-the-art lens modelling techniques have been applied to Hubble Space Telescope imaging data. We use i) the reconstructed images to infer the total and bulge luminosity of the host and ii) published broad-line spectroscopy to estimate $\mathcal M_{BH}$ using the so-called virial method. We then enlarge our sample with new calibration of previously published measurements to study the evolution of the correlation out to z~4.5. Consistent with previous work, we find that without taking into account passive luminosity evolution, the data points lie on the local relation. Once passive luminosity evolution is taken into account, we find that BHs in the more distant Universe reside in less luminous galaxies than today. Fitting this offset as $\mathcal M_{BH}$/$L_{host}$ $\propto$ (1+z)$^{\gamma}$, and taking into account selection effects, we obtain $\gamma$ = 0.6 $\pm$ 0.1 and 0.8$\pm$ 0.1 for the case of $\mathcal M_{BH}$-$L_{bulge}$ and $\mathcal M_{BH}$-$L_{total}$, respectively. To test for systematic uncertainties and selection effects we also consider a reduced sample that is homogeneous in data quality. We find consistent results but with considerably larger uncertainty due to the more limited sample size and redshift coverage ($\gamma$ = 0.7 $\pm$ 0.4 and 0.2$\pm$ 0.5 for $\mathcal M_{BH}$-$L_{bulge}$ and $\mathcal M_{BH}$-$L_{total}$, respectively), highlighting the need to gather more high-quality data for high-redshift lensed quasar hosts. Our result is consistent with a scenario where the growth of the black hole predates that of the host galaxy.

We investigate the intra-cluster light (ICL) in the 6 Hubble Frontier Field clusters at $0.3<z<0.6$. We employ a new method, which is free from any functional form of the ICL profile, and exploit the unprecedented depth of this Hubble Space Telescope imaging to map the ICL's diffuse light out to clustrocentric radii $R\sim300$kpc ($\mu_{\rm ICL}\sim27$mag arcsec$^{-2}$). From these maps, we construct radial color and stellar mass profiles via SED fitting and find clear negative color gradients in all systems with increasing distance from the Brightest Cluster Galaxy (BCG). While this implies older/more metal rich stellar components in the inner part of the ICL, we find the ICL mostly consists of a $<2$Gyr population, and plausibly originated with $\log M_*/M_\odot<10$ cluster galaxies. Further, we find 10-15% of the ICL's mass at large radii ($>150$kpc) lies in a younger/bluer stellar population ($\sim1$Gyr), a phenomenon not seen in local samples. We attribute this light to the higher fraction of starforming/(post-)starburst galaxies in clusters at $z\sim0.5$. Ultimately, we find the ICL's total mass to be $\log M_{\rm *}^{\rm ICL}/M_\odot\sim11$-12, constituting 5%-20% of the clusters' total stellar mass, or about a half of the value at $z\sim0$. The above implies distinct formation histories for the ICL and BCGs/other massive cluster galaxies; i.e. the ICL at this epoch is still being constructed rapidly ($\sim40M_\odot$yr$^{-1}$), while the BCGs have mostly completed their evolution. To be consistent with the ICL measurements of local massive clusters, such as the Virgo, our data suggest mass acquisition mainly from quiescent cluster galaxies is the principal source of ICL material in the subsequent $\sim$5 Gyr of cosmic time.

A massive galaxy cluster can serve as a magnifying glass for distant stellar populations, with strong gravitational lensing exposing details in the lensed background galaxies that would otherwise be undetectable. The MACS J0416.1-2403 cluster (hereafter MACS0416) is one of the most efficient lenses in the sky, and in 2014 it was observed with high-cadence imaging from the Hubble Space Telescope (HST). Here we describe two unusual transient events that appeared behind MACS0416 in a strongly lensed galaxy at redshift $z = 1.0054 \pm 0.0002$. These transients---designated HFF14Spo-NW and HFF14Spo-SE and collectively nicknamed "Spock"---were faster and fainter than any supernova (SN), but significantly more luminous than a classical nova. They reached peak luminosities of $\sim10^{41}$ erg s$^{-1}$ ($M_{\rm AB} < -14$ mag) in 5 rest-frame days, then faded below detectability in roughly the same time span. Models of the cluster lens suggest that these events may be spatially coincident at the source plane, but are most likely not temporally coincident. We find that HFF14Spo can be explained as a luminous blue variable (LBV), a recurrent nova (RN), or a pair of stellar microlensing events. To distinguish between these hypotheses will require a clarification of the positions of nearby critical curves, along with high-cadence monitoring of the field that could detect new transient episodes in the host galaxy.

Strong gravitational lensing provides a powerful test of Cold Dark Matter (CDM) as it enables the detection and mass measurement of low mass haloes even if they do not contain baryons. Compact lensed sources such as Active Galactic Nuclei (AGN) are particularly sensitive to perturbing subhalos, but their use as a test of CDM has been limited by the small number of systems which have significant radio emission which is extended enough avoid significant lensing by stars in the plane of the lens galaxy, and red enough to be minimally affected by differential dust extinction. Narrow-line emission is a promising alternative as it is also extended and, unlike radio, detectable in virtually all optically selected AGN lenses. We present first results from a WFC3 grism narrow-line survey of lensed quasars, for the quadruply lensed AGN HE0435-1223. Using a forward modelling pipeline which enables us to robustly account for spatial blending, we measure the [OIII] 5007 \AA~ flux ratios of the four images. We find that the [OIII] fluxes and positions are well fit by a simple smooth mass model for the main lens. Our data rule out a $M_{600}>10^{8} (10^{7.2}) M_\odot$ NFW perturber projected within $\sim$1\farcs0 (0\farcs1) arcseconds of each of the lensed images, where $M_{600}$ is the perturber mass within its central 600 pc. The non-detection is broadly consistent with the expectations of $\Lambda$CDM for a single system. The sensitivity achieved demonstrates that powerful limits on the nature of dark matter can be obtained with the analysis of $\sim20$ narrow-line lenses.

We present the discovery of 3 quasar lenses in the Sloan Digital Sky Survey (SDSS), selected using two novel photometry-based selection techniques. The J0941+0518 system, with two point sources separated by 5.46" on either side of a galaxy, has source and lens redshifts $z_s = 1.54$ and $z_l = 0.343$. The AO-assisted images of J2211+1929 show two point sources separated by 1.04", corresponding to the same quasar at $z_s = 1.07,$ besides the lens galaxy and Einstein ring. Images of J2257+2349 show two point sources separated by 1.67" on either side of an E/S0 galaxy. The extracted spectra show two images of the same quasar at redshift $z_s = 2.10$. In total, the two selection techniques identified 309 lens candidates, including 47 known lenses, and 6 previously ruled out candidates. 55 of the remaining candidates were observed using NIRC2 and ESI at Keck Observatory, EFOSC2 at the ESO-NTT (La Silla), and SAM and the Goodman spectrograph at SOAR. Of the candidates observed, 3 were confirmed as lenses, 36 were ruled out, and 16 remain inconclusive. Taking into account that we recovered known lenses, this gives us a success rate of at least 50/309 (16%). This initial campaign demonstrates the power of purely photometric selection techniques in finding lensed quasars. Developing and refining these techniques is essential for efficient identification of these rare lenses in ongoing and future photometric surveys.

Within one billion years of the Big Bang, intergalactic hydrogen was ionized by sources emitting ultraviolet and higher energy photons. This was the final phenomenon to globally affect all the baryons (visible matter) in the Universe. It is referred to as cosmic reionization and is an integral component of cosmology. It is broadly expected that intrinsically faint galaxies were the primary ionizing sources due to their abundance in this epoch. However, at the highest redshifts ($z>7.5$; lookback time 13.1 Gyr), all galaxies with spectroscopic confirmations to date are intrinsically bright and, therefore, not necessarily representative of the general population. Here, we report the unequivocal spectroscopic detection of a low luminosity galaxy at $z>7.5$. We detected the Lyman-$\alpha$ emission line at $\sim 10504$ {\AA} in two separate observations with MOSFIRE on the Keck I Telescope and independently with the Hubble Space Telescope's slit-less grism spectrograph, implying a source redshift of $z = 7.640 \pm 0.001$. The galaxy is gravitationally magnified by the massive galaxy cluster MACS J1423.8+2404 ($z = 0.545$), with an estimated intrinsic luminosity of $M_{AB} = -19.6 \pm 0.2$ mag and a stellar mass of $M_{\star} = 3.0^{+1.5}_{-0.8} \times 10^8$ solar masses. Both are an order of magnitude lower than the four other Lyman-$\alpha$ emitters currently known at $z > 7.5$, making it probably the most distant representative source of reionization found to date.

Swift monitoring of NGC 4151 with ~6 hr sampling over a total of 69 days in early 2016 is used to construct light curves covering five bands in the X-rays (0.3-50 keV) and six in the ultraviolet (UV)/optical (1900-5500 A). The three hardest X-ray bands (>2.5 keV) are all strongly correlated with no measurable interband lag while the two softer bands show lower variability and weaker correlations. The UV/optical bands are significantly correlated with the X-rays, lagging ~3-4 days behind the hard X-rays. The variability within the UV/optical bands is also strongly correlated, with the UV appearing to lead the optical by ~0.5-1 day. This combination of >~3 day lags between the X-rays and UV and <~1 day lags within the UV/optical appears to rule out the "lamp-post" reprocessing model in which a hot, X-ray emitting corona directly illuminates the accretion disk, which then reprocesses the energy in the UV/optical. Instead, these results appear consistent with the Gardner & Done picture in which two separate reprocessings occur: first, emission from the corona illuminates an extreme-UV-emitting toroidal component that shields the disk from the corona; this then heats the extreme-UV component which illuminates the disk and drives its variability.

We provide an updated calibration of CIV $\lambda1549$ broad emission line-based single-epoch (SE) black hole (BH) mass estimators for active galactic nuclei (AGNs) using new data for six reverberation-mapped AGNs at redshift $z=0.005-0.028$ with BH masses (bolometric luminosities) in the range $10^{6.5}-10^{7.5}$ $M_{\odot}$ ($10^{41.7}-10^{43.8}$ erg s$^{\rm -1}$). New rest-frame UV-to-optical spectra covering 1150-5700 \AA\ for the six AGNs were obtained with the Hubble Space Telescope (HST). Multi-component spectral decompositions of the HST spectra were used to measure SE emission-line widths for the CIV, MgII, and H$\beta$ lines as well as continuum luminosities in the spectral region around each line. We combine the new data with similar measurements for a previous archival sample of 25 AGNs to derive the most consistent and accurate calibrations of the CIV-based SE BH mass estimators against the H$\beta$ reverberation-based masses, using three different measures of broad-line width: full-width at half maximum (FWHM), line dispersion ($\sigma_{\rm line}$) and mean absolute deviation (MAD). The newly expanded sample at redshift $z=0.005-0.234$ covers a dynamic range in BH mass (bolometric luminosity) of $\log\ M_{\rm BH}/M_{\odot} = 6.5-9.1$ ($\log\ L_{\rm bol}/$erg s$^{\rm -1}=41.7-46.9$), and we derive the new CIV-based mass estimators using a Bayesian linear regression analysis over this range. We generally recommend the use of $\sigma_{\rm line}$ or MAD rather than FWHM to obtain a less biased velocity measurement of the CIV emission line, because its narrow-line component contribution is difficult to decompose from the broad-line profile.

We carry out a systematic investigation of the total mass density profile of massive (Mstar>2e11 Msun) early-type galaxies and its dependence on galactic properties and host halo mass with the aid of a variety of lensing/dynamical data and large mock galaxy catalogs. The latter are produced via semi-empirical models that, by design, are based on just a few basic input assumptions. Galaxies, with measured stellar masses, effective radii and S\'{e}rsic indices, are assigned, via abundance matching relations, host dark matter halos characterized by a typical LCDM profile. Our main results are as follows: (i) In line with observational evidence, our semi-empirical models naturally predict that the total, mass-weighted density slope at the effective radius gamma' is not universal, steepening for more compact and/or massive galaxies, but flattening with increasing host halo mass. (ii) Models characterized by a Salpeter or variable initial mass function and uncontracted dark matter profiles are in good agreement with the data, while a Chabrier initial mass function and/or adiabatic contractions/expansions of the dark matter halos are highly disfavored. (iii) Currently available data on the mass density profiles of very massive galaxies (Mstar>1e12 Msun), with Mhalo>3e14 Msun, favor instead models with a stellar profile flatter than a S\'{e}rsic one in the very inner regions (r<3-5 kpc), and a cored NFW or Einasto dark matter profile with median halo concentration a factor of ~2 or <1.3, respectively, higher than those typically predicted by N-body numerical simulations.

The CIII] and CIV rest-frame UV emission lines are powerful probes of the ionizations states of galaxies. They have furthermore been suggested as alternatives for spectroscopic redshift confirmation of objects at the epoch of reionization ($z>6$), where the most frequently used redshift indicator, Ly$\alpha$, is attenuated by the high fraction of neutral hydrogen in the inter-galactic medium. However, currently only very few confirmations of carbon UV lines at these high redshifts exist, making it challenging to quantify these claims. Here, we present the detection of CIV$\lambda\lambda$1548,1551\AA\ in \HST\ slitless grism spectroscopy obtained by GLASS of a Ly$\alpha$ emitter at $z=6.11$ multiply imaged by the massive foreground galaxy cluster RXJ2248. The CIV emission is detected at the 3--5$\sigma$ level in two images of the source, with marginal detection in two other images. We do not detect significant CIII]$\lambda\lambda$1907,1909\AA\ emission implying an equivalent width EW$_\textrm{CIII]}<20$\AA\ (1$\sigma$) and $\textrm{CIV/CIII}>0.7$ (2$\sigma$). Combined with limits on the rest-frame UV flux from the HeII$\lambda$1640\AA\ emission line and the OIII]$\lambda\lambda$1661,1666\AA\ doublet, we put constraints on the metallicity and the ionization state of the galaxy. The estimated line ratios and equivalent widths do not support a scenario where an AGN is responsible for ionizing the carbon atoms. SED fits including nebular emission lines imply a source with a mass of log(M/M$_\odot)\sim9$, SFR of around 10M$_\odot$/yr, and a young stellar population $<50$Myr old. The source shows a stronger ionizing radiation field than objects with detected CIV emission at $z<2$ and adds to the growing sample of low-mass (log(M/M$_\odot)\lesssim9$) galaxies at the epoch of reionization with strong radiation fields from star formation.

Exploiting the data of the Grism Lens-Amplified Survey from Space (GLASS), we characterize the spatial distribution of star formation in 76 high star forming galaxies in 10 clusters at 0.3< z <0.7. All these galaxies are likely restricted to first infall. In a companion paper we contrast the properties of field and cluster galaxies, whereas here we correlate the properties of H{\alpha} emitters to a number of tracers of the cluster environment to investigate its role in driving galaxy transformations. H{\alpha} emitters are found in the clusters out to 0.5 virial radii, the maximum radius covered by GLASS. The peak of the H{\alpha} emission is offset with respect to the peak of the UV-continuum. We decompose this offsets into a radial and tangential component. The radial compo- nent points away from the cluster center in 60% of the cases, with 95% confidence. The decompositions agree with cosmological simulations, i.e. the H{\alpha} emission offset correlates with galaxy velocity and ram-pressure stripping signatures. Trends between H{\alpha} emitter properties and surface mass density distributions and X-ray emissions emerge only for unrelaxed clusters. The lack of strong correlations with the global environment does not allow us to identify a unique environmental effect originating from the cluster center. In contrast, correla- tions between H{\alpha} morphology and local number density emerge. We conclude that local effects, uncorrelated to the cluster-centric radius, play a more important role in shaping galaxy properties.

We present the first results of the KMOS Lens-Amplified Spectroscopic Survey (KLASS), a new ESO Very Large Telescope (VLT) large program, doing multi-object integral field spectroscopy of galaxies gravitationally lensed behind seven galaxy clusters selected from the HST Grism Lens-Amplified Survey from Space (GLASS). Using the power of the cluster magnification we are able to reveal the kinematic structure of 25 galaxies at $0.7 \lesssim z \lesssim 2.3$, in four cluster fields, with stellar masses $8 \lesssim \log{(M_\star/M_\odot)} \lesssim 11$. This sample includes 5 sources at $z>1$ with lower stellar masses than in any previous kinematic IFU surveys. Our sample displays a diversity in kinematic structure over this mass and redshift range. The majority of our kinematically resolved sample is rotationally supported, but with a lower ratio of rotational velocity to velocity dispersion than in the local universe, indicating the fraction of dynamically hot disks changes with cosmic time. We find no galaxies with stellar mass $<3 \times 10^9 M_\odot$ in our sample display regular ordered rotation. Using the enhanced spatial resolution from lensing, we resolve a lower number of dispersion dominated systems compared to field surveys, competitive with findings from surveys using adaptive optics. We find that the KMOS IFUs recover emission line flux from HST grism-selected objects more faithfully than slit spectrographs. With artificial slits we estimate slit spectrographs miss on average 60% of the total flux of emission lines, which decreases rapidly if the emission line is spatially offset from the continuum.

We report the discovery of the quadruply lensed quasar J1433+6007, mined in the SDSS DR12 photometric catalogues using a novel outlier-selection technique, without prior spectroscopic or UV excess information. Discovery data obtained at the Nordic Optical telescope (NOT, La Palma) show nearly identical quasar spectra at $z_s=2.74$ and four quasar images in a fold configuration, one of which sits on a blue arc. The deflector redshift is $z_{l}=0.407,$ from Keck-ESI spectra. We describe the selection procedure, discovery and follow-up, image positions and $BVRi$ magnitudes, and first results and forecasts from simple lens models.

We present the results of an optical spectroscopic monitoring program targeting NGC 5548 as part of a larger multi-wavelength reverberation mapping campaign. The campaign spanned six months and achieved an almost daily cadence with observations from five ground-based telescopes. The H$\beta$ and He II $\lambda$4686 broad emission-line light curves lag that of the 5100 $\AA$ optical continuum by $4.17^{+0.36}_{-0.36}$ days and $0.79^{+0.35}_{-0.34}$ days, respectively. The H$\beta$ lag relative to the 1158 $\AA$ ultraviolet continuum light curve measured by the Hubble Space Telescope is roughly $\sim$50% longer than that measured against the optical continuum, and the lag difference is consistent with the observed lag between the optical and ultraviolet continua. This suggests that the characteristic radius of the broad-line region is $\sim$50% larger than the value inferred from optical data alone. We also measured velocity-resolved emission-line lags for H$\beta$ and found a complex velocity-lag structure with shorter lags in the line wings, indicative of a broad-line region dominated by Keplerian motion. The responses of both the H$\beta$ and He II $\lambda$4686 emission lines to the driving continuum changed significantly halfway through the campaign, a phenomenon also observed for C IV, Ly $\alpha$, He II(+O III]), and Si IV(+O IV]) during the same monitoring period. Finally, given the optical luminosity of NGC 5548 during our campaign, the measured H$\beta$ lag is a factor of five shorter than the expected value implied by the $R_\mathrm{BLR} - L_\mathrm{AGN}$ relation based on the past behavior of NGC 5548.

We present gravitational lens models of the multiply imaged quasar DES J0408-5354, recently discovered in the Dark Energy Survey (DES) footprint, with the aim of interpreting its remarkable quad-like configuration. We first model the DES single-epoch $grizY$ images as a superposition of a lens galaxy and four point-like objects, obtaining spectral energy distributions (SEDs) and relative positions for the objects. Three of the point sources (A,B,D) have SEDs compatible with the discovery quasar spectra, while the faintest point-like image (G2/C) shows significant reddening and a `grey' dimming of $\approx0.8$mag. In order to understand the lens configuration, we fit different models to the relative positions of A,B,D. Models with just a single deflector predict a fourth image at the location of G2/C but considerably brighter and bluer. The addition of a small satellite galaxy ($R_{\rm E}\approx0.2$") in the lens plane near the position of G2/C suppresses the flux of the fourth image and can explain both the reddening and grey dimming. All models predict a main deflector with Einstein radius between $1.7"$ and $2.0",$ velocity dispersion $267-280$km/s and enclosed mass $\approx 6\times10^{11}M_{\odot},$ even though higher resolution imaging data are needed to break residual degeneracies in model parameters. The longest time-delay (B-A) is estimated as $\approx 85$ (resp. $\approx125$) days by models with (resp. without) a perturber near G2/C. The configuration and predicted time-delays of J0408-5354 make it an excellent target for follow-up aimed at understanding the source quasar host galaxy and substructure in the lens, and measuring cosmological parameters. We also discuss some lessons learnt from J0408-5354 on lensed quasar finding strategies, due to its chromaticity and morphology.

We report the discovery and spectroscopic confirmation of the quad-like lensed quasar system DES J0408-5354 found in the Dark Energy Survey (DES) Year 1 (Y1) data. This system was discovered during a search for DES Y1 strong lensing systems using a method that identified candidates as red galaxies with multiple blue neighbors. DES J0408-5354 consists of a central red galaxy surrounded by three bright (i < 20) blue objects and a fourth red object. Subsequent spectroscopic observations using the Gemini South telescope confirmed that the three blue objects are indeed the lensed images of a quasar with redshift z = 2.375, and that the central red object is an early-type lensing galaxy with redshift z = 0.597. DES J0408-5354 is the first quad lensed quasar system to be found in DES and begins to demonstrate the potential of DES to discover and dramatically increase the sample size of these very rare objects.

We present a new measurement of the Hubble Constant H0 and other cosmological parameters based on the joint analysis of three multiply-imaged quasar systems with measured gravitational time delays. First, we measure the time delay of HE0435-1223 from 13-year light curves obtained as part of the COSMOGRAIL project. Companion papers detail the modeling of the main deflectors and line of sight effects, and how these data are combined to determine the time-delay distance of HE 0435-1223. Crucially, the measurements are carried out blindly with respect to cosmological parameters in order to avoid confirmation bias. We then combine the time-delay distance of HE0435-1223 with previous measurements from systems B1608+656 and RXJ1131-1231 to create a Time Delay Strong Lensing probe (TDSL). In flat $\Lambda$CDM with free matter and energy density, we find $H_0$ = 71.9 +2.4 -3.0 km/s/Mpc and $\Omega_{\Lambda}$ = 0.62 +0.24 -0.35 . This measurement is completely independent of, and in agreement with, the local distance ladder measurements of H0. We explore more general cosmological models combining TDSL with other probes, illustrating its power to break degeneracies inherent to other methods. The TDSL and Planck joint constraints are $H_0$ = 69.2 +1.4 -2.2 km/s/Mpc, $\Omega_{\Lambda}$ = 0.70 +0.01 -0.01 and $\Omega_k$ = 0.003 +0.004 -0.006 in open $\Lambda$CDM and $H_0$ = 79.0 +4.4 -4.2 km/s/Mpc, $\Omega_{de}$ = 0.77 +0.02 -0.03 and $w$ = -1.38 +0.14 -0.16 in flat $w$CDM. Combined with Planck and Baryon Acoustic Oscillation data, when relaxing the constraints on the numbers of relativistic species we find $N_{eff}$ = 3.34 +0.21 -0.21 and when relaxing the total mass of neutrinos we find 0.182 eV. In an open $w$CDM in combination with Planck and CMB lensing we find $H_0$ = 77.9 +5.0 -4.2 km/s/Mpc, $\Omega_{de}$ = 0.77 +0.03 -0.03, $\Omega_k$ = -0.003 +0.004 -0.004 and $w$ = -1.37 +0.18 -0.23.

Galaxies located in the environment or on the line of sight towards gravitational lenses can significantly affect lensing observables, and can lead to systematic errors on the measurement of $H_0$ from the time-delay technique. We present the results of a systematic spectroscopic identification of the galaxies in the field of view of the lensed quasar HE0435-1223, using the W. M. Keck, Gemini and ESO-Very Large telescopes. Our new catalog triples the number of known galaxy redshifts in the vicinity of the lens, expanding to 102 the number of measured redshifts for galaxies separated by less than 3 arcmin from the lens. We complement our catalog with literature data to gather redshifts up to 15 arcmin from the lens, and search for galaxy groups or clusters projected towards HE0435-1223. We confirm that the lens is a member of a small group that includes at least 12 galaxies, and find 8 other group candidates near the line of sight of the lens. The flexion shift, namely the shift of lensed images produced by high order perturbation of the lens potential, is calculated for each galaxy/group and used to identify which objects produce the largest perturbation of the lens potential. This analysis demonstrates that i) at most three of the five brightest galaxies projected within 12 arcsec of the lens need to be explicitly used in the lens models, and ii) the groups can be treated in the lens model as an external tidal field (shear) contribution.

Strong gravitational lens systems with time delays between the multiple images allow measurements of time-delay distances, which are primarily sensitive to the Hubble constant that is key to probing dark energy, neutrino physics, and the spatial curvature of the Universe, as well as discovering new physics. We present H0LiCOW ($H_0$ Lenses in COSMOGRAIL's Wellspring), a program that aims to measure $H_0$ with $<3.5\%$ uncertainty from five lens systems (B1608+656, RXJ1131-1231, HE0435-1223, WFI2033-4723 and HE1104-1805). We have been acquiring (1) time delays through COSMOGRAIL and Very Large Array monitoring, (2) high-resolution Hubble Space Telescope imaging for the lens mass modeling, (3) wide-field imaging and spectroscopy to characterize the lens environment, and (4) moderate-resolution spectroscopy to obtain the stellar velocity dispersion of the lenses for mass modeling. In cosmological models with one-parameter extension to flat $\Lambda$CDM, we expect to measure $H_0$ to $<3.5\%$ in most models, spatial curvature $\Omega_{\rm k}$ to 0.004, $w$ to 0.14, and the effective number of neutrino species to 0.2 (1$\sigma$ uncertainties) when combined with current CMB experiments. These are, respectively, a factor of $\sim15$, $\sim2$, and $\sim1.5$ tighter than CMB alone. Our data set will further enable us to study the stellar initial mass function of the lens galaxies, and the co-evolution of supermassive black holes and their host galaxies. This program will provide a foundation for extracting cosmological distances from the hundreds of time-delay lenses that are expected to be discovered in current and future surveys.

We present the results of ALMA spectroscopic follow-up of a $z=6.765$ Lyman-$\alpha$ emitting galaxy behind the cluster RXJ1347-1145. We report the detection of [CII]158$\mu$m line fully consistent with the Lyman-$\alpha$ redshift and with the peak of the optical emission. Given the magnification of $\mu=5.0 \pm 0.3$ the intrinsic (corrected for lensing) luminosity of the [CII] line is $L_{[CII]} =1.4^{+0.2}_{-0.3} \times 10^7L_{\odot}$, which is ${\sim}5$ times fainter than other detections of $z\sim 7$ galaxies. The result indicates that low $L_{[CII]}$ in $z\sim 7$ galaxies compared to the local counterparts might be caused by their low metallicities and/or feedback. The small velocity off-set ($\Delta v = 20_{-40}^{+140} \rm km/s$) between the Lyman-$\alpha$ and [CII] line is unusual, and may be indicative of ionizing photons escaping.

We present an exploration of the mass structure of a sample of 12 strongly lensed massive, compact early-type galaxies at redshifts $z\sim0.6$ to provide further possible evidence for their inside-out growth. We obtain new ESI/Keck spectroscopy and infer the kinematics of both lens and source galaxies, and combine these with existing photometry to construct (a) the fundamental plane (FP) of the source galaxies and (b) physical models for their dark and luminous mass structure. We find their FP to be tilted towards the virial plane relative to the local FP, and attribute this to their unusual compactness, which causes their kinematics to be totally dominated by the stellar mass as opposed to their dark matter; that their FP is nevertheless still inconsistent with the virial plane implies that both the stellar and dark structure of early-type galaxies is non-homologous. We also find the intrinsic scatter of their FP to be comparable to the local value, indicating that variations in the stellar mass structure outweight variations in the dark halo in the central regions of early-type galaxies. Finally, we show that inference on the dark halo structure -- and, in turn, the underlying physics -- is sensitive to assumptions about the stellar initial mass function (IMF), but that physically-motivated assumptions about the IMF imply haloes with sub-NFW inner density slopes, and may present further evidence for the inside-out growth of compact early-type galaxies via minor mergers and accretion.

Using deep Hubble Frontier Fields imaging and slitless spectroscopy from the Grism Lens-Amplified Survey from Space, we analyze 2200 cluster and 1748 field galaxies at $0.2\leq z\leq0.7$ to determine the impact of environment on galaxy size and structure at $\log M_*/M_\odot>7.8$, an unprecedented limit at these redshifts. Based on simple assumptions-$r_e=f(M_*)$-we find no significant differences in half-light radii ($r_e$) between equal-mass cluster or field systems. More complex analyses-$r_e=f(M_*,U-V,n,z,\Sigma$)-reveal local density $(\Sigma$) to induce only a $7\% \pm 3\%$ ($95\%$ confidence) reduction in $r_e$ beyond what can be accounted for by $U-V$ color, Sersic index ($n$), and redshift ($z$) effects.Almost any size difference between galaxies in high- and low-density regions is thus attributable to their different distributions in properties other than environment. Indeed, we find a clear color-$r_e$ correlation in low-mass passive cluster galaxies ($\log M_*/M_\odot<9.8$) such that bluer systems have larger radii, with the bluest having sizes consistent with equal-mass star-forming galaxies. We take this as evidence that large-$r_e$ low-mass passive cluster galaxies are recently acquired systems that have been environmentally quenched without significant structural transformation (e.g., by ram pressure stripping or starvation).Conversely, $\sim20\%$ of small-$r_e$ low-mass passive cluster galaxies appear to have been in place since $z\sim3$. Given the consistency of the small-$r_e$ galaxies' stellar surface densities (and even colors) with those of systems more than ten times as massive, our findings suggest that clusters mark places where galaxy evolution is accelerated for an ancient base population spanning most masses, with late-time additions quenched by environment-specific mechanisms are mainly restricted to the lowest masses.

Strong gravitational lenses with measured time delays between the multiple images allow a direct measurement of the time-delay distance to the lens, and thus a measure of cosmological parameters, particularly the Hubble constant, $H_{0}$. We present a blind lens model analysis of the quadruply-imaged quasar lens HE 0435-1223 using deep Hubble Space Telescope imaging, updated time-delay measurements from the COSmological MOnitoring of GRAvItational Lenses (COSMOGRAIL), a measurement of the velocity dispersion of the lens galaxy based on Keck data, and a characterization of the mass distribution along the line of sight. HE 0435-1223 is the third lens analyzed as a part of the $H_{0}$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) project. We account for various sources of systematic uncertainty, including the detailed treatment of nearby perturbers, the parameterization of the galaxy light and mass profile, and the regions used for lens modeling. We constrain the effective time-delay distance to be $D_{\Delta t} = 2612_{-191}^{+208}~\mathrm{Mpc}$, a precision of 7.6%. From HE 0435-1223 alone, we infer a Hubble constant of $H_{0} = 73.1_{-6.0}^{+5.7}~\mathrm{km~s^{-1}~Mpc^{-1}}$ assuming a flat $\Lambda$CDM cosmology. The cosmographic inference based on the three lenses analyzed by H0LiCOW to date is presented in a companion paper (H0LiCOW Paper V).

Wide-field photometric surveys enable searches of rare yet interesting objects, such as strongly lensed quasars or quasars with a bright host galaxy. Past searches for lensed quasars based on their optical and near infrared properties have relied on photometric cuts and spectroscopic pre-selection (as in the Sloan Quasar Lens Search), or neural networks applied to photometric samples. These methods rely on cuts in morphology and colours, with the risk of losing many interesting objects due to scatter in their population properties, restrictive training sets, systematic uncertainties in catalog-based magnitudes, and survey-to-survey photometric variations. Here, we explore the performance of a Gaussian Mixture Model to separate point-like quasars, quasars with an extended host, and strongly lensed quasars using griz psf and model magnitudes and WISE W1, W2. The choice of optical magnitudes is due to their presence in all current and upcoming releases of wide-field surveys, whereas UV information is not always available. We then assess the contamination from blue galaxies and the role of additional features such as W3 magnitudes or psf-model terms as morphological information. As a demonstration, we conduct a search in a random 10% of the SDSS footprint, and we provide the catalog of the 43 SDSS object with the highest `lens' score in our selection that survive visual inspection, and are spectroscopically confirmed to host active nuclei. We inspect archival data and find images of 5/43 objects in the Hubble Legacy Archive, including 2 known lenses. The code and materials are available to facilitate follow-up.

The Type~Ia supernova (SN~Ia) 2016coj in NGC 4125 (redshift $z=0.004523$) was discovered by the Lick Observatory Supernova Search 4.9 days after the fitted first-light time (FFLT; 11.1 days before $B$-band maximum). Our first detection (pre-discovery) is merely $0.6\pm0.5$ day after the FFLT, making SN 2016coj one of the earliest known detections of a SN Ia. A spectrum was taken only 3.7 hr after discovery (5.0 days after the FFLT) and classified as a normal SN Ia. We performed high-quality photometry, low- and high-resolution spectroscopy, and spectropolarimetry, finding that SN 2016coj is a spectroscopically normal SN Ia, but with a high velocity of \ion{Si}{2} $\lambda$6355 ($\sim 12,600$\,\kms\ around peak brightness). The \ion{Si}{2} $\lambda$6355 velocity evolution can be well fit by a broken-power-law function for up to a month after the FFLT. SN 2016coj has a normal peak luminosity ($M_B \approx -18.9 \pm 0.2$ mag), and it reaches a $B$-band maximum \about16.0~d after the FFLT. We estimate there to be low host-galaxy extinction based on the absence of Na~I~D absorption lines in our low- and high-resolution spectra. The spectropolarimetric data exhibit weak polarization in the continuum, but the \ion{Si}{2} line polarization is quite strong ($\sim 0.9\% \pm 0.1\%$) at peak brightness.

We present a new sample of strong gravitational lens systems where both the foreground lenses and background sources are early-type galaxies. Using imaging from HST/ACS and Keck/NIRC2, we model the surface brightness distributions and show that the sources form a distinct population of massive, compact galaxies at redshifts $0.4 \lesssim z \lesssim 0.7$, lying systematically below the size-mass relation of the global elliptical galaxy population at those redshifts. These may therefore represent relics of high-redshift red nuggets or their partly-evolved descendants. We exploit the magnifying effect of lensing to investigate the structural properties, stellar masses and stellar populations of these objects with a view to understanding their evolution. We model these objects parametrically and find that they generally require two S\'ersic components to properly describe their light profiles, with one more spheroidal component alongside a more envelope-like component, which is slightly more extended though still compact. This is consistent with the hypothesis of the inside-out growth of these objects via minor mergers. We also find that the sources can be characterised by red-to-blue colour gradients as a function of radius which are stronger at low redshift -- indicative of ongoing accretion -- but that their environments generally appear consistent with that of the general elliptical galaxy population, contrary to recent suggestions that these objects are predominantly associated with clusters.

The arrival times, positions, and fluxes of multiple images in strong lens systems can be used to infer the presence of dark subhalos in the deflector, and thus test predictions of cold dark matter models. However, gravitational lensing does not distinguish between perturbations to a smooth gravitational potential arising from baryonic and non-baryonic mass. In this work, we quantify the extent to which the stellar mass distribution of a deflector can reproduce flux ratio and astrometric anomalies typically associated with the presence of a dark matter subhalo. Using Hubble Space Telescope images of nearby galaxies, we simulate strong lens systems with real distributions of stellar mass as they would be observed at redshift $z_d=0.5$. We add a dark matter halo and external shear to account for the smooth dark matter field, omitting dark substructure, and use a Monte Carlo procedure to characterize the distributions of image positions, time delays, and flux ratios for a compact background source of diameter 5 pc. By convolving high-resolution images of real galaxies with a Gaussian PSF, we simulate the most detailed smooth potential one could construct given high quality data, and find scatter in flux ratios of $\approx 10\%$, which we interpret as a typical deviation from a smooth potential caused by large and small scale structure in the lensing galaxy. We demonstrate that the flux ratio anomalies arising from galaxy-scale baryonic structure can be minimized by selecting the most massive and round deflectors, and by simultaneously modeling flux ratio and astrometric data.

The empirical correlation between the mass of a super-massive black hole (MBH) and its host galaxy properties is widely considered to be evidence of their co-evolution. A powerful way to test the co-evolution scenario and learn about the feedback processes linking galaxies and nuclear activity is to measure these correlations as a function of redshift. Unfortunately, currently MBH can only be estimated in active galaxies at cosmological distances. At these distances, bright active galactic nuclei (AGN) can outshine the host galaxy, making it extremely difficult to measure the host's luminosity. Strongly lensed AGNs provide in principle a great opportunity to improve the sensitivity and accuracy of the host galaxy luminosity measurements as the host galaxy is magnified and more easily separated from the point source, provided the lens model is sufficiently accurate. In order to measure the MBH-L correlation with strong lensing, it is necessary to ensure that the lens modelling is accurate, and that the host galaxy luminosity can be recovered to at least a precision and accuracy better than that of the typical MBH measurement. We carry out extensive and realistic simulations of deep Hubble Space Telescope observations of lensed AGNs obtained by our collaboration. We show that the host galaxy luminosity can be recovered with better accuracy and precision than the typical uncertainty on MBH(~ 0.5 dex) for hosts as faint as 2-4 magnitudes dimmer than the AGN itself. Our simulations will be used to estimate bias and uncertainties on the actual measurements to be presented in a future paper.

(Abridged) We combine deep HST grism spectroscopy with a new Bayesian method to derive maps of gas-phase metallicity, nebular dust extinction, and star-formation rate for 10 star-forming galaxies at high redshift ($1.2<z<2.3$). Exploiting lensing magnification by the foreground cluster MACS1149.6+2223, we reach sub-kpc spatial resolution and push the stellar mass limit associated with such high-z spatially resolved measurements below $10^8M_\odot$ for the first time. Our maps exhibit diverse morphologies, indicative of various effects such as efficient radial mixing from tidal torques, rapid accretion of low-metallicity gas, etc., which can affect the gas and metallicity distributions in individual galaxies. Based upon an exhaustive sample of all existing sub-kpc metallicity gradients at high-z, we find that predictions given by analytical chemical evolution models assuming a relatively extended star-formation profile in the early disk formation phase can explain the majority of observed gradients, without involving galactic feedback or radial outflows. We observe a tentative correlation between stellar mass and metallicity gradient, consistent with the downsizing galaxy formation picture that more massive galaxies are more evolved into a later phase of disk growth, where they experience more coherent mass assembly at all radii and thus show shallower metallicity gradients. In addition, we compile a sample of homogeneously cross-calibrated integrated metallicity measurements spanning three orders of magnitude in stellar mass at $z\sim1.8$. We use this sample to study the mass-metallicity relation (MZR) and test the fundamental metallicity relation (FMR). The slope of the observed MZR can rule out the momentum-driven wind model at 3-$\sigma$ confidence level. We find no significant offset with respect to the FMR, taking into account the intrinsic scatter and measurement uncertainties.

We present a detection in pre-explosion Hubble Space Telescope (HST) imaging of a point source consistent with being the progenitor star of the Type IIb supernova (SN IIb) 2016gkg. Post-explosion imaging from the Keck Adaptive Optics system was used to perform relative astrometry between the Keck and HST imaging. We identify a single point source in the HST images coincident with the SN position to 0.89-sigma. The HST photometry is consistent with the progenitor star being an A0Ia star with T=9500 K and log (L/Lsun)=5.15. We find that the SN 2016gkg progenitor star appears more consistent with binary than single-star evolutionary models. In addition, early-time light curve data from SN 2016gkg revealed a rapid rise in luminosity within ~0.4 days of non-detection limits, consistent with models of the cooling phase after shock break-out. We use these data to determine an explosion date of 20.15 September 2016 and progenitor star radius of log (R/Rsun)=2.41, which agrees with photometry from the progenitor star. Our findings are also consistent with detections of other SNe IIb progenitor stars, although more luminous and bluer than most other examples.

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.

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.

The stellar initial mass function (IMF) of early-type galaxies is the combination of the IMF of the stellar population formed in-situ and that of accreted stellar populations. Using as an observable the effective IMF $\alpha_{IMF}$, defined as the ratio between the true stellar mass of a galaxy and the stellar mass inferred assuming a Salpeter IMF, we present a theoretical model for its evolution as a result of dry mergers. We use a simple dry merger evolution model, based on cosmological $N$-body simulations, together with empirically motivated prescriptions for the IMF to make predictions for how the effective IMF of massive early-type galaxies changes from $z=2$ to $z=0$. We find that the IMF normalization of individual galaxies becomes lighter with time. At fixed velocity dispersion, $\alpha_{IMF}$ is predicted to be constant with redshift. Current constraints on the evolution of the IMF are in slight tension with this prediction, even though systematic uncertainties prevent a conclusive statement. The correlation of $\alpha_{IMF}$ with stellar mass becomes shallower with time, while the correlation between $\alpha_{IMF}$ and velocity dispersion is mostly preserved by dry mergers. We also find that dry mergers can mix the dependence of the IMF on stellar mass and velocity dispersion, making it challenging to infer, from $z=0$ observations of global galactic properties, what is the quantity that is originally coupled with the IMF.

Based on spectroscopy and multiband wide-field observations of the gravitationally lensed quasar HE 0435-1223, we determine the probability distribution function of the external convergence $\kappa_\mathrm{ext}$ for this system. We measure the under/overdensity of the line of sight towards the lens system and compare it to the average line of sight throughout the universe, determined by using the CFHTLenS as a control field. Aiming to constrain $\kappa_\mathrm{ext}$ as tightly as possible, we determine under/overdensities using various combinations of relevant informative weighing schemes for the galaxy counts, such as projected distance to the lens, redshift, and stellar mass. We then convert the measured under/overdensities into a $\kappa_\mathrm{ext}$ distribution, using ray-tracing through the Millennium Simulation. We explore several limiting magnitudes and apertures, and account for systematic and statistical uncertainties relevant to the quality of the observational data, which we further test through simulations. Our most robust estimate of $\kappa_\mathrm{ext}$ has a median value $\kappa^\mathrm{med}_\mathrm{ext} = 0.004$ and a standard deviation of $\sigma_\kappa = 0.025$. The measured $\sigma_\kappa$ corresponds to $2.5\%$ uncertainty on the time delay distance, and hence the Hubble constant $H_0$ inference from this system. The median $\kappa^\mathrm{med}_\mathrm{ext}$ value is robust to $\sim0.005$ (i.e. $\sim0.5\%$ on $H_0$) regardless of the adopted aperture radius, limiting magnitude and weighting scheme, as long as the latter incorporates galaxy number counts, the projected distance to the main lens, and a prior on the external shear obtained from mass modeling. The availability of a well-constrained $\kappa_\mathrm{ext}$ makes \hequad\ a valuable system for measuring cosmological parameters using strong gravitational lens time delays.

The Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST) enabled the search for the first galaxies observed at z ~ 8 - 11 (500 - 700 Myr after the Big Bang). To continue quantifying the number density of the most luminous galaxies (M_AB ~ -22.0) at the earliest epoch observable with HST, we search for z ~ 10 galaxies (F125W-dropouts) in archival data from the Brightest of Reionizing Galaxies (BoRG[z8]) survey, originally designed for detection of z ~ 8 galaxies (F098M-dropouts). By focusing on the deepest 293 arcmin^2 of the data along 62 independent lines of sight, we identify six z ~ 10 candidates satisfying the color selection criteria, detected at S/N > 8 in F160W with M_AB = -22.8 to -21.1 if at z = 10. Three of the six sources, including the two brightest, are in a single WFC3 pointing (~ 4 arcmin^2), suggestive of significant clustering, which is expected from bright galaxies at z ~ 10. However, the two brightest galaxies are too extended to be likely at z ~ 10, and one additional source is unresolved and possibly a brown dwarf. The remaining three candidates have m_AB ~ 26, and given the area and completeness of our search, our best estimate is a number density of sources that is marginally higher but consistent at 2{\sigma} with searches in legacy fields. Our study highlights that z ~ 10 searches can yield a small number of candidates, making tailored follow-ups of HST pure-parallel observations viable and effective.

We report the detection of Ly$\alpha$ emission at $\sim9538$\AA{} 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$\alpha$ luminosity to be $\sim1.3\times10^{42}~\mathrm{erg}~\mathrm{s}^{-1}$ and Ly$\alpha$ equivalent width to be $74\pm15$\AA{}. Its intrinsic UV absolute magnitude at 1600\AA{} 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$\alpha$ 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$\alpha$ escape fraction is high ($\gtrsim 10$\%), which could be common for sub-$L^*$ $z\gtrsim7$ galaxies with Ly$\alpha$ 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.

Gravitational time delays, observed in strong lens systems where the variable background source is multiply-imaged by a massive galaxy in the foreground, provide direct measurements of cosmological distance that are very complementary to other cosmographic probes. The success of the technique depends on the availability and size of a suitable sample of lensed quasars or supernovae, precise measurements of the time delays, accurate modeling of the gravitational potential of the main deflector, and our ability to characterize the distribution of mass along the line of sight to the source. We review the progress made during the last 15 years, during which the first competitive cosmological inferences with time delays were made, and look ahead to the potential of significantly larger lens samples in the near future.

In Hubble Space Telescope (HST) imaging taken on 10 November 2014, four images of supernova (SN) "Refsdal" (redshift z=1.49) appeared in an Einstein-cross--like configuration (images S1-S4) around an early-type galaxy in the cluster MACS J1149.5+2223 (z=0.54). Almost all lens models of the cluster have predicted that the SN should reappear within a year in a second host-galaxy image created by the cluster's potential. In HST observations taken on 11 December 2015, we find a new source at the predicted position of the new image of SN Refsdal approximately 8" from the previous images S1-S4. This marks the first time the appearance of a SN at a particular time and location in the sky was successfully predicted in advance! We use these data and the light curve from the first four observed images of SN Refsdal to place constraints on the relative time delay and magnification of the new image (SX), compared to images S1-S4. This enables us, for the first time, to test "blind" lens model predictions of both magnifications and time delays for a lensed SN. We find that the timing and brightness of the new image are consistent with the blind predictions of a fraction of the models. The reappearance illustrates the discriminatory power of this blind test and its utility to uncover sources of systematic uncertainty. From planned HST photometry, we expect to reach a precision of 1-2% on the time delay between S1-S4 and SX.

During an intensive Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS) UV monitoring campaign of the Seyfert~1 galaxy NGC 5548 performed from 2014 February to July, the normally highly correlated far-UV continuum and broad emission-line variations decorrelated for ~60 to 70 days, starting ~75 days after the first HST/COS observation. Following this anomalous state, the flux and variability of the broad emission lines returned to a more normal state. This transient behavior, characterised by significant deficits in flux and equivalent width of the strong broad UV emission lines, is the first of its kind to be unambiguously identified in an active galactic nucleus reverberation mapping campaign. The largest corresponding emission-line flux deficits occurred for the high-ionization collisionally excited lines, C IV and Si IV(+O IV]), and also He II(+O III]), while the anomaly in Ly-alpha was substantially smaller. This pattern of behavior indicates a depletion in the flux of photons with E_{\rm ph} > 54 eV, relative to those near 13.6 eV. We suggest two plausible mechanisms for the observed behavior: (i) temporary obscuration of the ionizing continuum incident upon BLR clouds by a moving veil of material lying between the inner accretion disk and inner BLR, perhaps resulting from an episodic ejection of material from the disk, or (ii) a temporary change in the intrinsic ionizing continuum spectral energy distribution resulting in a deficit of ionizing photons with energies > 54 eV, possibly due to a transient restructuring of the Comptonizing atmosphere above the disk. Current evidence appears to favor the latter explanation.

We present MUSE observations in the core of the HFF galaxy cluster MACS J1149.5+2223, where the first magnified and spatially-resolved multiple images of SN 'Refsdal' at redshift 1.489 were detected. Thanks to a DDT program with the VLT and the extraordinary efficiency of MUSE, we measure 117 secure redshifts with just 4.8 hours of total integration time on a single target pointing. We spectroscopically confirm 68 galaxy cluster members, with redshift values ranging from 0.5272 to 0.5660, and 18 multiple images belonging to 7 background, lensed sources distributed in redshifts between 1.240 and 3.703. Starting from the combination of our catalog with those obtained from extensive spectroscopic and photometric campaigns using the HST, we select a sample of 300 (164 spectroscopic and 136 photometric) cluster members, within approximately 500 kpc from the BCG, and a set of 88 reliable multiple images associated to 10 different background source galaxies and 18 distinct knots in the spiral galaxy hosting SN 'Refsdal'. We exploit this valuable information to build 6 detailed strong lensing models, the best of which reproduces the observed positions of the multiple images with a rms offset of only 0.26". We use these models to quantify the statistical and systematic errors on the predicted values of magnification and time delay of the next emerging image of SN 'Refsdal'. We find that its peak luminosity should should occur between March and June 2016, and should be approximately 20% fainter than the dimmest (S4) of the previously detected images but above the detection limit of the planned HST/WFC3 follow-up. We present our two-dimensional reconstruction of the cluster mass density distribution and of the SN 'Refsdal' host galaxy surface brightness distribution. We outline the roadmap towards even better strong lensing models with a synergetic MUSE and HST effort.

It is widely believed that the large discrepancy between the observed number of satellite galaxies and the predicted number of dark subhalos can be resolved via a variety of baryonic effects which suppress star formation in low mass halos.Supporting this hypothesis, numerous high resolution simulations with star formation, and associated feedback have been shown to reproduce the satellite luminosity function around Milky Way-mass simulated galaxies at redshift zero. However, a more stringent test of these models is their ability to simultaneously match the satellite luminosity functions of a range of host halo masses and redshifts. In this work we measure the luminosity function of faint (sub-Small Magellanic Cloud luminosity) satellites around hosts with stellar masses 10.5$<\log_{10}$M$_*$/M$_\odot<11.5$ to an unprecedented redshift of 1.5. This new measurement of the satellite luminosity function provides powerful new constraining power; we compare these results with predictions from four different simulations and show that although the models perform similarly over-all, no one model reproduces the satellite luminosity function reliably at all redshifts and host stellar masses. This result highlights the continued need for improvement in understanding the fundamental physics that governs satellite galaxy evolution.

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.

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$\beta$ line variability in a sample of low-redshift Seyfert galaxies.

We present ground-based optical photometric monitoring data for NGC 5548, part of an extended multi-wavelength reverberation mapping campaign. The light curves have nearly daily cadence from 2014 January to July in nine filters (\emph{BVRI} and \emph{ugriz}). Combined with ultraviolet data from the \emph{Hubble Space Telescope} and \emph{Swift}, we confirm significant time delays between the continuum bands as a function of wavelength, extending the wavelength coverage from 1158\,\AA\ to the $z$ band ($\sim\!9160$\,\AA). We find that the lags at wavelengths longer than the {\it V} band are equal to or greater than the lags of high-ionization-state emission lines (such as He\,{\sc ii}\,$\lambda 1640$ and $\lambda 4686$), suggesting that the continuum-emitting source is of a physical size comparable to the inner broad-line region (BLR). The trend of lag with wavelength is broadly consistent with the prediction for continuum reprocessing by an accretion disk with $\tau \propto \lambda^{4/3}$. However, the lags also imply a disk radius that is 3 times larger than the prediction from standard thin-disk theory, assuming that the bolometric luminosity is 10\% of the Eddington luminosity ($L = 0.1L_{\rm Edd}$). Using optical spectra from the Large Binocular Telescope, we estimate the bias of the interband continuum lags due to BLR emission observed in the filters. We find that the bias for filters with high levels of BLR contamination ($\sim\! 20\%$) can be important for the shortest continuum lags, and likely has a significant impact on the {\it u} and {\it U} bands owing to Balmer continuum emission.

We present the first year of Hubble Space Telescope imaging of the unique supernova (SN) 'Refsdal', a gravitationally lensed SN at z=1.488$\pm$0.001 with multiple images behind the galaxy cluster MACS J1149.6+2223. The first four observed images of SN Refsdal (images S1-S4) exhibited a slow rise (over ~150 days) to reach a broad peak brightness around 20 April, 2015. Using a set of light curve templates constructed from SN 1987A-like peculiar Type II SNe, we measure time delays for the four images relative to S1 of 4$\pm$4 (for S2), 2$\pm$5 (S3), and 24$\pm$7 days (S4). The measured magnification ratios relative to S1 are 1.15$\pm$0.05 (S2), 1.01$\pm$0.04 (S3), and 0.34$\pm$0.02 (S4). None of the template light curves fully captures the photometric behavior of SN Refsdal, so we also derive complementary measurements for these parameters using polynomials to represent the intrinsic light curve shape. These more flexible fits deliver fully consistent time delays of 7$\pm$2 (S2), 0.6$\pm$3 (S3), and 27$\pm$8 days (S4). The lensing magnification ratios are similarly consistent, measured as 1.17$\pm$0.02 (S2), 1.00$\pm$0.01 (S3), and 0.38$\pm$0.02 (S4). We compare these measurements against published predictions from lens models, and find that the majority of model predictions are in very good agreement with our measurements. Finally, we discuss avenues for future improvement of time delay measurements -- both for SN Refsdal and for other strongly lensed SNe yet to come.

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\mu$m and $4.5\mu$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$\alpha$ emitters in our sample from the Keck DEIMOS spectra, one at $z_{\text{Ly}\alpha}=6.76$ (in RXJ1347) and one at $z_{\text{Ly}\alpha}=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.

We have acquired Hubble Space Telescope (HST) and Very Large Telescope near-infrared spectra and images of supernova (SN) Refsdal after its discovery as an Einstein cross in Fall 2014. The HST light curve of SN Refsdal matches the distinctive, slowly rising light curves of SN 1987A-like supernovae (SNe), and we find strong evidence for a broad H-alpha P-Cygni profile in the HST grism spectrum at the redshift (z = 1.49) of the spiral host galaxy. SNe IIn, powered by circumstellar interaction, could provide a good match to the light curve of SN Refsdal, but the spectrum of a SN IIn would not show broad and strong H-alpha absorption. From the grism spectrum, we measure an H-alpha expansion velocity consistent with those of SN 1987A-like SNe at a similar phase. The luminosity, evolution, and Gaussian profile of the H-alpha emission of the WFC3 and X-shooter spectra, separated by ~2.5 months in the rest frame, provide additional evidence that supports the SN 1987A-like classification. In comparison with other examples of SN 1987A-like SNe, SN Refsdal has a blue B-V color and a high luminosity for the assumed range of potential magnifications. If SN Refsdal can be modeled as a scaled version of SN 1987A, we estimate it would have an ejecta mass of 20+-5 solar masses. The evolution of the light curve at late times will provide additional evidence about the potential existence of any substantial circumstellar material (CSM). Using MOSFIRE and X-shooter spectra, we estimate a subsolar host-galaxy metallicity (8.3+-0.1 dex and <8.4 dex, respectively) near the explosion site.

We present the first results and design from the redshift z~9-10 Brightest of the Reionizing Galaxies {\it Hubble Space Telescope} survey BoRG[z9-10], aimed at searching for intrinsically luminous unlensed galaxies during the first 700 Myr after the Big Bang. BoRG[z9-10] is the continuation of a multi-year pure-parallel near-IR and optical imaging campaign with the Wide Field Camera 3. The ongoing survey uses five filters, optimized for detecting the most distant objects and offering continuous wavelength coverage from {\lambda}=0.35{\mu}m to {\lambda}=1.7{\mu}m. We analyze the initial ~130 arcmin$^2$ of area over 28 independent lines of sight (~25% of the total planned) to search for z>7 galaxies using a combination of Lyman break and photometric redshift selections. From an effective comoving volume of (5-25) $times 10^5$ Mpc$^3$ for magnitudes brighter than $m_{AB}=26.5-24.0$ in the $H_{160}$-band respectively, we find five galaxy candidates at z~8.3-10 detected at high confidence (S/N>8), including a source at z~8.4 with mAB=24.5 (S/N~22), which, if confirmed, would be the brightest galaxy identified at such early times (z>8). In addition, BoRG[z9-10] data yield four galaxies with $7.3 \lesssim z \lesssim 8$. These new Lyman break galaxies with m$\lesssim26.5$ are ideal targets for follow-up observations from ground and space based observatories to help investigate the complex interplay between dark matter growth, galaxy assembly, and reionization.

Supernova "Refsdal," multiply imaged by cluster MACSJ1149.5+2223, represents a rare opportunity to make a true blind test of model predictions in extragalactic astronomy, on a time scale that is short compared to a human lifetime. In order to take advantage of this event, we produced seven gravitational lens models with five independent methods, based on Hubble Space Telescope (HST) Hubble Frontier Field images, along with extensive spectroscopic follow-up observations by HST, the Very Large and the Keck Telescopes. We compare the model predictions and show that they agree reasonably well with the measured time delays and magnification ratios between the known images, even though these quantities were not used as input. This agreement is encouraging, considering that the models only provide statistical uncertainties, and do not include additional sources of uncertainties such as structure along the line of sight, cosmology, and the mass sheet degeneracy. We then present the model predictions for the other appearances of SN "Refsdal." A future image will reach its peak in the first half of 2016, while another image appeared between 1994 and 2004. The past image would have been too faint to be detected in existing archival images. The future image should be approximately one third as bright as the brightest known image (i.e., H_AB~25.7 mag at peak and H_AB~26.7 mag six months before peak), and thus detectable in single-orbit HST images. We will find out soon whether our predictions are correct.

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 $\tau = 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 ($\alpha\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.

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.

[abbreviated] We present a census of Ly\alpha\ emission at $z\gtrsim7$ utilizing deep near infrared HST grism spectroscopy from the first six completed clusters of the Grism Lens-Amplified Survey from Space (GLASS). In 24/159 photometrically selected galaxies we detect emission lines consistent with Ly\alpha\ in the GLASS spectra. Based on the distribution of signal-to-noise ratios and on simulations we expect the completeness and the purity of the sample to be 40-100% and 60-90%, respectively. For the objects without detected emission lines we show that the observed (not corrected for lensing magnification) 1$\sigma$ flux limits reaches $5\times10^{-18}$erg/s/cm$^{2}$ per position angle over the full wavelength range of GLASS (0.8-1.7$\mu$m). Based on the conditional probability of Ly\alpha\ emission measured from the ground at $z\sim7$ we would have expected 12-18 Ly\alpha\ emitters. This is consistent with the number of detections, within the uncertainties, confirming the drop in Ly\alpha\ emission with respect to $z\sim6$. These candidates include a promising source at $z=8.1$. The spatial extent of Ly\alpha\ in a deep stack of the most convincing Ly\alpha\ emitters with $\langle z\rangle=7.2$ is consistent with that of the rest-frame UV continuum. Extended Ly$\alpha$ emission, if present, has a surface brightness below our detection limit, consistent with the properties of lower redshift comparison samples. From the stack we estimate upper limits on rest-frame UV emission line ratios and find $f_\textrm{CIV} / f_\textrm{Ly${\alpha}$} \lesssim 0.32$ and $f_\textrm{CIII]} / f_\textrm{Ly$\alpha$} \lesssim 0.23$ in good agreement with other values published in the literature.

Supernova (SN) Refsdal is the first multiply-imaged, highly-magnified, and spatially-resolved SN ever observed. The SN exploded in a highly-magnified spiral galaxy at z=1.49 behind the Frontier Fields Cluster MACS1149, and provides a unique opportunity to study the environment of SNe at high z. We exploit the time delay between multiple images to determine the properties of the SN and its environment, before, during, and after the SN exploded. We use the integral-field spectrograph MUSE on the VLT to simultaneously target all observed and model-predicted positions of SN Refsdal. We find MgII emission at all positions of SN Refsdal, accompanied by weak FeII* emission at two positions. The measured ratios of [OII] to MgII emission of 10-20 indicate a high degree of ionization with low metallicity. Because the same high degree of ionization is found in all images, and our spatial resolution is too coarse to resolve the region of influence of SN Refsdal, we conclude that this high degree of ionization has been produced by previous SNe or a young and hot stellar population. We find no variability of the [OII] line over a period of 57 days. This suggests that there is no variation in the [OII] luminosity of the SN over this period, or that the SN has a small contribution to the integrated [OII] emission over the scale resolved by our observations.

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{\alpha} 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{\alpha} 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{\alpha} 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.

We present the first results from the Las Cumbres Observatory Global Telescope (LCOGT) Network's Active Galactic Nuclei Key Project, a large program devoted to using the robotic resources of LCOGT to perform time domain studies of active galaxies. We monitored the Seyfert 1 galaxy Arp~151 (Mrk~40) for $\sim$200 days with robotic imagers and with the FLOYDS robotic spectrograph at Faulkes Telescope North. Arp~151 was highly variable during this campaign, with $V$-band light curve variations of $\sim$0.3 mag and H$\beta$ flux changing by a factor of $\sim$3. We measure robust time lags between the $V$-band continuum and the H$\alpha$, H$\beta$ and H$\gamma$ emission lines, with $\tau_\mathrm{cen} = 13.89^{+1.39}_{-1.41}$, 7.52$^{+1.43}_{-1.06}$ and 7.40$^{+1.50}_{-1.32}$ days, respectively. The lag for the \ion{He}{2} $\lambda4686$ emission line is unresolved. We measure a velocity-resolved lag for the H$\beta$ line, which is clearly asymmetric with higher lags on the blue wing of the line which decline to the red, possibly indicative of radial inflow, and is similar in morphology to past observations of the H$\beta$ transfer function shape. Assuming a virialization factor of $f$=5.5, we estimate a black hole mass of $M_\mathrm{BH}=6.2^{+1.4}_{-1.2}\times$10$^{6}$~$M_{\odot}$, also consistent with past measurements for this object. These results represent the first step to demonstrate the powerful robotic capabilities of LCOGT for long-term, AGN time domain campaigns that human intensive programs cannot easily accomplish. Arp 151 is now one of just a few AGN where the virial product is known to remain constant against substantial changes in H$\beta$ lag and luminosity.

We give an overview of the Grism Lens Amplified Survey from Space (GLASS), a large Hubble Space Telescope program aimed at obtaining grism spectroscopy of the fields of ten massive clusters of galaxies at redshift z=0.308-0.686, including the Hubble Frontier Fields (HFF). The Wide Field Camera 3 yields near infrared spectra of the cluster cores, covering the wavelength range 0.81-1.69mum through grisms G102 and G141, while the Advanced Camera for Surveys in parallel mode provides G800L spectra of the infall regions of the clusters. The WFC3 spectra are taken at two almost orthogonal position angles in order to minimize the effects of confusion. After summarizing the scientific drivers of GLASS, we describe the sample selection as well as the observing strategy and data processing pipeline. We then utilize MACSJ0717.5+3745, a HFF cluster and the first one observed by GLASS, to illustrate the data quality and the high-level data products. Each spectrum brighter than H_AB=23 is visually inspected by at least two co-authors and a redshift is measured when sufficient information is present in the spectra. Furthermore, we conducted a thorough search for emission lines through all the GLASS WFC3 spectra with the aim of measuring redshifts for sources with continuum fainter than H_AB=23. We provide a catalog of 139 emission-line based spectroscopic redshifts for extragalactic sources, including three new redshifts of multiple image systems (one probable, two tentative). In addition to the data itself we also release software tools that are helpful to navigate the data.

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 $\gamma'$ of the total mass density profile $\rho_{tot}\propto r^{-\gamma'}$ for the main deflector can be measured. Ideally, we require that the statistical error on $\gamma'$ 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.

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}\,\beta$ broad emission line in five nearby AGN, as well as for $\rm{H}\,\alpha$ and $\rm{H}\,\gamma$ 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}\,\beta$. 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.

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.

We present a strong and weak lensing reconstruction of the massive cluster Abell 2744, the first cluster for which deep Hubble Frontier Field (HFF) images and spectroscopy from the Grism Lens-Amplified Survey from Space (GLASS) are available. By performing a targeted search for emission lines in multiply imaged sources using the GLASS spectra, we obtain 5 high-confidence spectroscopic redshifts and 2 tentative ones. We confirm 1 strongly lensed system by detecting the same emission lines in all 3 multiple images. We also search for additional line emitters blindly and use the full GLASS spectroscopic catalog to test reliability of photometric redshifts for faint line emitters. We see a reasonable agreement between our photometric and spectroscopic redshift measurements, when including nebular emission in photometric redshift estimations. We introduce a stringent procedure to identify only secure multiple image sets based on colors, morphology, and spectroscopy. By combining 7 multiple image systems with secure spectroscopic redshifts (at 5 distinct redshift planes) with 18 multiple image systems with secure photometric redshifts, we reconstruct the gravitational potential of the cluster pixellated on an adaptive grid, using a total of 72 images. The resulting mass map is compared with a stellar mass map obtained from the deep Spitzer Frontier Fields data to study the relative distribution of stars and dark matter in the cluster. We find that the stellar to total mass ratio varies substantially across the cluster field, suggesting that stars do not trace exactly the total mass in this interacting system. The maps of convergence, shear, and magnification are made available in the standard HFF format.

We present a detailed analysis of an individual case of gravitational lensing of a $z\sim8$ Lyman-Break galaxy (LBG) in a blank field, identified in Hubble Space Telescope imaging obtained as part of the Brightest of Reionizing Galaxies survey. To investigate the close proximity of the bright ($m_{AB}=25.8$) $Y_{098}$-dropout to a small group of foreground galaxies, we obtained deep spectroscopy of the dropout and two foreground galaxies using VLT/X-Shooter. We detect H-$\alpha$, H-$\beta$, [OIII] and [OII] emission in the brightest two foreground galaxies (unresolved at the natural seeing of $0.8$ arcsec), placing the pair at $z=1.327$. We can rule out emission lines contributing all of the observed broadband flux in $H_{160}$ band at $70\sigma$, allowing us to exclude the $z\sim8$ candidate as a low redshift interloper with broadband photometry dominated by strong emission lines. The foreground galaxy pair lies at the peak of the luminosity, redshift and separation distributions for deflectors of strongly lensed $z\sim8$ objects, and we make a marginal detection of a demagnified secondary image in the deepest ($J_{125}$) filter. We show that the configuration can be accurately modelled by a singular isothermal ellipsoidal deflector and a S\'{e}rsic source magnified by a factor of $\mu=4.3\pm0.2$. The reconstructed source in the best-fitting model is consistent with luminosities and morphologies of $z\sim8$ LBGs in the literature. The lens model yields a group mass of $9.62\pm0.31\times10^{11} M_{\odot}$ and a stellar mass-to-light ratio for the brightest deflector galaxy of $M_{\star}/L_{B}=2.3^{+0.8}_{-0.6} M_{\odot}/L_{\odot}$ within its effective radius. The foreground galaxies' redshifts would make this one of the few strong lensing deflectors discovered at $z>1$.

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.

We present the X-shooter Lens Survey (XLENS) data. The main goal of XLENS is to disentangle the stellar and dark matter content of massive early-type galaxies (ETGs), through combined strong gravitational lensing, dynamics and spectroscopic stellar population studies. The sample consists of 11 lens galaxies covering the redshift range from $0.1$ to $0.45$ and having stellar velocity dispersions between $250$ and $380\,\mathrm{km}\,\mathrm{s}^{-1}$. All galaxies have multi-band, high-quality HST imaging. We have obtained long-slit spectra of the lens galaxies with X-shooter on the VLT. We are able to disentangle the dark and luminous mass components by combining lensing and extended kinematics data-sets, and we are also able to precisely constrain stellar mass-to-light ratios and infer the value of the low-mass cut-off of the IMF, by adding spectroscopic stellar population information. Our goal is to correlate these IMF parameters with ETG masses and investigate the relation between baryonic and non-baryonic matter during the mass assembly and structure formation processes. In this paper we provide an overview of the survey, highlighting its scientific motivations, main goals and techniques. We present the current sample, briefly describing the data reduction and analysis process, and we present the first results on spatially resolved kinematics.

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.

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.

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 $\sigma=(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 $\gamma^{\prime}=-\log\rho/\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.

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

We describe the first results from a six-month long reverberation-mapping experiment in the ultraviolet based on 170 observations of the Seyfert 1 galaxy NGC 5548 with the Cosmic Origins Spectrograph on the Hubble Space Telescope. Significant correlated variability is found in the continuum and broad emission lines, with amplitudes ranging from ~30% to a factor of two in the emission lines and a factor of three in the continuum. The variations of all the strong emission lines lag behind those of the continuum, with He II 1640 lagging behind the continuum by ~2.5 days and Lyman alpha 1215, C IV 1550, and Si IV 1400 lagging by ~5-6 days. The relationship between the continuum and emission lines is complex. In particular, during the second half of the campaign, all emission-line lags increased by a factor of 1.3-2 and differences appear in the detailed structure of the continuum and emission-line light curves. Velocity-resolved cross-correlation analysis shows coherent structure in lag versus line-of-sight velocity for the emission lines; the high-velocity wings of C IV respond to continuum variations more rapidly than the line core, probably indicating higher velocity BLR clouds at smaller distances from the central engine. The velocity-dependent response of Lyman alpha, however, is more complex and will require further analysis.

Recent intensive Swift monitoring of the Seyfert 1 galaxy NGC 5548 yielded 282 usable epochs over 125 days across six UV/optical bands and the X-rays. This is the densest extended AGN UV/optical continuum sampling ever obtained, with a mean sampling rate <0.5 day. Approximately daily HST UV sampling was also obtained. The UV/optical light curves show strong correlations (r_max = 0.57 - 0.90) and the clearest measurement to date of interband lags. These lags are well-fit by a \tau propto \lambda^4/3 wavelength dependence, with a normalization that indicates an unexpectedly large disk radius of 0.35 +/- 0.05 lt-day at 1367 A, assuming a simple face-on model. The U-band shows a marginally larger lag than expected from the fit and surrounding bands, which could be due to Balmer continuum emission from the broad-line region as suggested by Korista and Goad. The UV/X-ray correlation is weaker (r_max < 0.45) and less consistent over time. This indicates that while Swift is beginning to measure UV/optical lags in general agreement with accretion disk theory (although the derived size is larger than predicted), the relationship with X-ray variability is less well understood. Combining this accretion disk size estimate with those from quasar microlensing studies suggests that AGN disk sizes scale approximately linearly with central black hole mass over a wide range of masses.

We present new results from our search for z~7 galaxies from deep spectroscopic observations of candidate z-dropouts in the CANDELS fields. Despite the extremely low flux limits achieved by our sensitive observations, only 2 galaxies have robust redshift identifications, one from its Lyalpha emission line at z=6.65, the other from its Lyman-break, i.e. the continuum discontinuity at the Lyalpha wavelength consistent with a redshift 6.42, but with no emission line. In addition, for 23 galaxies we present deep limits in the Lyalpha EW derived from the non detections in ultra-deep observations. Using this new data as well as previous samples, we assemble a total of 68 candidate z~7 galaxies with deep spectroscopic observations, of which 12 have a line detection. With this much enlarged sample we can place solid constraints on the declining fraction of Ly$\alpha$ emission in z~7 Lyman break galaxies compared to z~6, both for bright and faint galaxies. Applying a simple analytical model, we show that the present data favor a patchy reionization process rather than a smooth one.

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 $\alpha = -1.72^{+0.30}_{-0.29}$, and number density of $\log_{10} \Psi^\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 $\mu>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.

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 ($\mu\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\sigma}$ 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\sigma$).

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.

In the Spring of 2011 we carried out a 2.5 month reverberation mapping campaign using the 3 m Shane telescope at Lick Observatory, monitoring 15 low-redshift Seyfert 1 galaxies. This paper describes the observations, reductions and measurements, and data products from the spectroscopic campaign. The reduced spectra were fitted with a multicomponent model in order to isolate the contributions of various continuum and emission-line components. We present light curves of broad emission lines and the AGN continuum, and measurements of the broad H-beta line widths in mean and root-mean square (rms) spectra. For the most highly variable AGNs we also measured broad H-beta line widths and velocity centroids from the nightly spectra. In four AGNs exhibiting the highest variability amplitudes, we detect anticorrelations between broad H-beta width and luminosity, demonstrating that the broad-line region "breathes" on short timescales of days to weeks in response to continuum variations. We also find that broad H-beta velocity centroids can undergo substantial changes in response to continuum variations; in NGC 4593 the broad H-beta velocity shifted by ~250 km/s over a one-month duration. This reverberation-induced velocity shift effect is likely to contribute a significant source of confusion noise to binary black hole searches that use multi-epoch quasar spectroscopy to detect binary orbital motion. We also present results from simulations that examine biases that can occur in measurement of broad-line widths from rms spectra due to the contributions of continuum variations and photon-counting noise.

Detection of Lyman-Break Galaxies (LBGs) at high-redshift can be affected by gravitational lensing induced by foreground deflectors not only in galaxy clusters, but also in blank fields. We quantify the impact of strong magnification in the samples of $B$, $V$, $i$, $z$ $\&$ $Y$ LBGs ($4\lesssim z \lesssim8$) observed in the XDF and GOODS/CANDELS fields, by investigating the proximity of dropouts to foreground objects. We find that $\sim6\%$ of bright LBGs ($m_{H_{160}}<26$) at $z\sim7$ have been strongly lensed ($\mu>2$) by foreground objects. This fraction decreases from $\sim 3.5\%$ at $z\sim6$ to $\sim1.5\%$ at $z\sim4$. Since the observed fraction of strongly lensed galaxies is a function of the shape of the luminosity function (LF), it can be used to derive Schechter parameters, $\alpha$ and $M_{\star}$, independently from galaxy number counts. Our magnification bias analysis yields Schechter-function parameters in close agreement with those determined from galaxy counts albeit with larger uncertainties. Extrapolation of our analysis to $z\gtrsim 8$ suggests that future surveys with JSWT, WFIRST and EUCLID should find excess LBGs at the bright-end, even if there is an intrinsic exponential cutoff of number counts. Finally, we highlight how the magnification bias measurement near the detection limit can be used as probe of the population of galaxies too faint to be detected. Preliminary results using this novel idea suggest that the magnification bias at $M_{UV}\sim -18$ is not as strong as expected if $\alpha\lesssim -1.7$ extends well below the current detection limits in the XDF. At face value this implies a flattening of the LF at $M_{UV}\gtrsim-16.5$. However, selection effects and completeness estimates are difficult to quantify precisely. Thus, we do not rule out a steep LF extending to $M_{UV}\gtrsim -15$.

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} {\sigma}_*$-$\log_{10} {\sigma}_{\rm SIE}$ relation, which corresponds to a 6-$\sigma$ 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.

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.

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$\beta$ 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$\beta$ 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/$\sigma$ of the emission line profile.

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.

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.

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.

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.

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 $\chi^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 $\chi^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.

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.

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$\beta$ broad emission line combined with the 5100 \AA\ 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)^{\gamma}$ with $\gamma=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 $\gamma=0.9\pm0.7$ for the $M_{\rm BH}-L_{\rm bul}$ relation and $\gamma=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.

Central Galaxies (CGs) in massive halos live in unique environments with formation histories closely linked to that of the host halo. In local clusters they have larger sizes ($R_e$) and lower velocity dispersions (sigma) at fixed stellar mass M_star, and much larger R_e at a fixed $\sigma$ than field and satellite galaxies (non-CGs). Using spectroscopic observations of group galaxies selected from the COSMOS survey, we compare the dynamical scaling relations of early-type CGs and non-CGs at z~0.6, to distinguish possible mechanisms that produce the required evolution. CGs are systematically offset towards larger R_e at fixed $\sigma$ compared to non-CGs with similar M_star. The CG R_e-M_star relation also shows differences, primarily driven by a sub-population (~15%) of galaxies with large $R_e$, while the M_star-sigma relations are indistinguishable. These results are accentuated when double Sersic profiles, which better fit light in the outer regions of galaxies, are adopted. They suggest that even group-scale CGs can develop extended components by these redshifts that can increase total $R_e$ and M_star estimates by factors of ~2. To probe the evolutionary link between our sample and cluster CGs, we also analyze two cluster samples at z~0.6 and z~0. We find similar results for the more massive halos at comparable z, but much more distinct CG scaling relations at low-z. Thus, the rapid, late-time accretion of outer components, perhaps via the stripping and accretion of satellites, would appear to be a key feature that distinguishes the evolutionary history of CGs.

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.

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.

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 ($\sigma$) within the SDSS fibers. Comparing the dynamically-derived stellar mass-to-light ratios $(M_*/L)_{\rm dyn}$, obtained for an assumed halo slope $\rho_{\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 $\sigma$ galaxies is consistent on average with a Salpeter slope. Our study allows for a fully consistent study of the trend between IMF and $\sigma$ for both the SLACS and \ATLAS samples, which explore quite different $\sigma$ ranges. The two samples are highly complementary, the first being essentially $\sigma$ 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\alpha =(0.38\pm0.04)\times\log(\sigma_{\rm e}/200\,\mathrm{km~s}^{-1})+(-0.06\pm0.01)$, where $\alpha=(M_*/L)_{\rm dyn}/(M_*/L)_{\rm pop}$ and $\sigma_{\rm e}$ is the luminosity averaged $\sigma$ 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 $\sigma_{\rm e}\approx90-270\,\mathrm{km~s}^{-1}$.

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.

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.

We present the results of a search for galaxy substructures in a sample of 11 gravitational lens galaxies from the Sloan Lens ACS Survey. We find no significant detection of mass clumps, except for a luminous satellite in the system SDSS J0956+5110. We use these non-detections, in combination with a previous detection in the system SDSS J0946+1006, to derive constraints on the substructure mass function in massive early-type host galaxies with an average redshift z ~ 0.2 and an average velocity dispersion of 270 km/s. We perform a Bayesian inference on the substructure mass function, within a median region of about 32 kpc squared around the Einstein radius (~4.2 kpc). We infer a mean projected substructure mass fraction $f = 0.0076^{+0.0208}_{-0.0052}$ at the 68 percent confidence level and a substructure mass function slope $\alpha$ < 2.93 at the 95 percent confidence level for a uniform prior probability density on alpha. For a Gaussian prior based on Cold Dark Matter (CDM) simulations, we infer $f = 0 .0064^{+0.0080}_{-0.0042}$ and a slope of $\alpha$ = 1.90$^{+0.098}_{-0.098}$ at the 68 percent confidence level. Since only one substructure was detected in the full sample, we have little information on the mass function slope, which is therefore poorly constrained (i.e. the Bayes factor shows no positive preference for any of the two models).The inferred fraction is consistent with the expectations from CDM simulations and with inference from flux ratio anomalies at the 68 percent confidence level.

Under the assumption of a flat Lambda-CDM cosmology, recent data from the Planck satellite point toward a Hubble constant that is in tension with that measured by gravitational lens time delays and by the local distance ladder. Prosaically, this difference could arise from unknown systematic uncertainties in some of the measurements. More interestingly -- if systematics were ruled out -- resolving the tension would require a departure from the flat Lambda-CDM cosmology, introducing for example a modest amount of spatial curvature, or a non-trivial dark energy equation of state. To begin to address these issues, we present here an analysis of the gravitational lens RXJ1131-1231 that is improved in one particular regard: we examine the issue of systematic error introduced by an assumed lens model density profile. We use more flexible gravitational lens models with baryonic and dark matter components, and find that the exquisite Hubble Space Telescope image with thousands of intensity pixels in the Einstein ring and the stellar velocity dispersion of the lens contain sufficient information to constrain these more flexible models. The total uncertainty on the time-delay distance is 6.6% for a single system. We proceed to combine our improved time-delay distance measurements with the WMAP9 and Planck posteriors. In an open Lambda-CDM model, the data for RXJ1131-1231 in combination with Planck favor a flat universe with Omega_k = 0.00+0.01/-0.02 (68% CI). In a flat wCDM model, the combination of RXJ1131-1231 and Planck yields w = -1.52+0.19/-0.20 (68% CI).

Strong gravitational lenses can be used to detect low mass subhalos, based on deviations in image fluxes and positions from what can be achieved with a smooth mass distribution. So far, this method has been limited by the small number of (radio-loud, microlensing free) systems which can be analysed for the presence of substructure. Using the gravitational lens B1422+231, we demonstrate that adaptive optics integral field spectroscopy can also be used to detect dark substructures. We analyse data obtained with OSIRIS on the Keck I Telescope, using a Bayesian method that accounts for uncertainties relating to the point spread function and image positions in the separate exposures. The narrow-line [OIII] fluxes measured for the lensed images are consistent with those measured in the radio, and show a significant deviation from what would be expected in a smooth mass distribution, consistent with the presence of a perturbing low mass halo. Detailed lens modelling shows that image fluxes and positions are fit significantly better when the lens is modelled as a system containing a single perturbing subhalo in addition to the main halo, rather than by the main halo on its own, indicating the significant detection of substructure.The inferred mass of the subhalo depends on the subhalo mass density profile: the 68% confidence interval for the perturber mass within 600 pc are: 8.2$^{+0.6}_{-0.8}$, 8.2$^{+0.6}_{-1}$ and 7.6$\pm0.3$ $\log_{10}$[M$_{\rm{sub}}$/M$_\odot$] respectively for a singular isothermal sphere, a pseudo-Jaffe, and an NFW mass profile. This method can extend the study of flux ratio anomalies to virtually all quadruply imaged quasars, and therefore offers great potential to improve the determination of the subhalo mass function in the near future.

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.

Using deep two-color near-infrared HST imaging and unbiased grism spectroscopy we present a detailed study of the z=1.803 JKCS041 cluster. Uniquely, for a high redshift cluster, we confirm a mass of $\log M=14.2$ in solar units using four different techniques based on the X-ray temperature, the X-ray luminosity and the cluster richness. JKCS041 is thus a progenitor of a local system like the Coma cluster. Our rich dataset and the abundant population of 14 spectroscopically-confirmed red sequence galaxies allows us to explore the past star formation history of this system in unprecedented detail. Remarkably, we find a prominent red sequence down to stellar masses as low as $\log M=9.8$, corresponding to a mass range of 2 dex. These quiescent galaxies are concentrated around the cluster center with a core radius of 330 kpc. Blue members are few and avoid the cluster center. In JKCS041 quenching was therefore largely completed by a look-back time of 10 Gyr and we can constrain the epoch at which this occurred via spectroscopic age-dating of the individual galaxies. Most galaxies were quenched about 1.1 Gyr prior to the epoch of observation. The less massive quiescent galaxies are somewhat younger, corresponding to a decrease in age of 650 Myr per mass dex, but the scatter in age at fixed mass is only 380 Myr (at $\log M=11$). The size-mass relation of quiescent galaxies in JKCS041 is consistent with that observed for local clusters within our uncertainties. Comparing our data on JKCS041 with 41 clusters at lower redshift, we find that the form of the mass function of red sequence galaxies has hardly evolved in the past 10 Gyr, both in terms of its faint end slope and characteristic mass. Despite observing JKCS041 soon after its quenching and the three fold expected increase in mass in the next 10 Gyr, it is already remarkably similar to present-day clusters.

Galaxy masses play a fundamental role in our understanding of structure formation models. This review addresses the variety and reliability of mass estimators that pertain to stars, gas, and dark matter. The different sections on masses from stellar populations, dynamical masses of gas-rich and gas-poor galaxies, with some attention paid to our Milky Way, and masses from weak and strong lensing methods, all provide review material on galaxy masses in a self-consistent manner.

We present Spitzer/IRAC observations of nine $z'$-band dropouts highly magnified (2<mu<12) by the Bullet Cluster. We combine archival imaging with our Exploratory program (SURFS UP), which results in a total integration time of ~30 hr per IRAC band. We detect (>3sigma) in both IRAC bands the brightest of these high-redshift galaxies, with [3.6]=23.80+-0.28 mag, [4.5]=23.78+-0.25 mag, and (H-[3.6])=1.17+-0.32 mag. The remaining eight galaxies are undetected to [3.6]~26.4 mag and [4.5]~26.0 mag with stellar masses of ~5x10^7 M_sol. The detected galaxy has an estimated magnification of mu=12+-4, which implies this galaxy has an ultraviolet luminosity of L_1500~0.3 L*_{z=7} --- the lowest luminosity individual source detected in IRAC at z>7. By modeling the broadband photometry, we estimate the galaxy has an intrinsic star-formation rate of SFR~1.3 M_sol/yr and stellar mass of M~2x10^9 M_sol, which gives a specific star-formation rate of sSFR~0.7 Gyr^-1. If this galaxy had sustained this star-formation rate since z~20, it could have formed the observed stellar mass (to within a factor of ~2), we also discuss alternate star-formation histories and argue the exponentially-increasing model is unlikely. Finally, based on the intrinsic star-formation rate, we estimate this galaxy has a likely [C II] flux of <f_[C II]> = 10^{-17} erg/s/cm2.

[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 $\alpha = -1.87^{+0.26}_{-0.26}$, and a number density of $\log_{10} \phi^\star [\textrm{Mpc}^{-3}] = -3.24^{+0.25}_{-0.24}$. Integrated down to $M=-17.7$ this LF yields a luminosity density, $\log_{10} \epsilon [\textrm{erg}/\textrm{s/Hz/Mpc}^{3}] = 25.52^{+0.05}_{-0.05}$. Our LF analysis is consistent with previously published determinations within 1$\sigma$. 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$.

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 \gamma' driven by dry mergers occurring at rates given by cosmological simulations. Such dry-merger model predicts a strong decrease of \gamma' with cosmic time, inconsistent with the almost constant \gamma' 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 \gamma'. 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.

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$\mu$m and 4.5$\mu$m in the central 4 arcmin x 4 arcmin field and ~15 hours in the flanking fields. This results in 3-$\sigma$ sensitivity limits of ~26.6 and ~26.2AB magnitudes for the central field in the IRAC 3.6 and 4.5$\mu$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.

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.

Recent studies have shown that massive elliptical galaxies have total mass density profiles within an effective radius that can be approximated as \rho_{tot}\propto r^{-\gamma'}, with mean slope <\gamma'>=2.08 \pm 0.03 and scatter \sigma_\gamma'=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 \gamma' 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 \gamma' 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 \gamma'. Models with halo expansion, or mass-follows-light predict too high values of \gamma'. 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.

The Grism Lens-Amplified Survey from Space (GLASS) is a Hubble Space Telescope (HST) Large Program, which will obtain 140 orbits of grism spectroscopy of the core and infall regions of 10 galaxy clusters, selected to be among the very best cosmic telescopes. Extensive HST imaging is available from many sources including the CLASH and Frontier Fields programs. We introduce the survey by analyzing spectra of faint multiply-imaged galaxies and $z\gtrsim6$ galaxy candidates obtained from the first seven orbits out of fourteen targeting the core of the Frontier Fields cluster MACS0717.5+3745. Using the G102 and G141 grisms to cover the wavelength range 0.8-1.7$\mu$m, we confirm 4 strongly lensed systems by detecting emission lines in each of the images. For the 9 $z\gtrsim6$ galaxy candidates clear from contamination, we do not detect any emission lines down to a seven-orbit 1$\sigma$ noise level of $\sim$5$\times$10$^{-18}$erg s$^{-1}$cm$^{-2}$. Taking lensing magnification into account, our flux sensitivity reaches $\sim$0.2-5$\times$10$^{-18}$erg s$^{-1}$cm$^{-2}$. These limits over an uninterrupted wavelength range rule out the possibility that the high-$z$ galaxy candidates are instead strong line emitters at lower redshift. These results show that by means of careful modeling of the background - and with the assistance of lensing magnification - interesting flux limits can be reached for large numbers of objects, avoiding pre-selection and the wavelength restrictions inherent to ground-based multi-slit spectroscopy. These observations confirm the power of slitless HST spectroscopy even in fields as crowded as a cluster core.

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.

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.

Strong gravitational lensing can be used to directly measure the mass function of their satellites, thus testing one of the fundamental predictions of cold dark matter cosmological models. Given the importance of this test it is essential to ensure that galaxies acting as strong lenses have dark and luminous satellites which are representative of the overall galaxy population. We address this issue by measuring the number and spatial distribution of luminous satellites in ACS imaging around lens galaxies from the Sloan Lens Advanced Camera for Surveys (SLACS) lenses, and comparing them with the satellite population in ACS imaging of non lens galaxies selected from COSMOS, which has similar depth and resolution to the ACS images of SLACS lenses. In order to compare the samples of lens and non lens galaxies, which have intrinsically different stellar mass distributions, we measure, for the first time, the number of satellites per host as a continuous function of host stellar mass for both populations. We find that the number of satellites as a function of host stellar mass, as well as the spatial distribution are consistent between the samples. Using these results, we predict the number of satellites we would expect to find around a subset of the Cosmic Lens All Sky Survey (CLASS) lenses, and find a result consistent with the the number observed by Jackson et al. 2010. Thus we conclude that within our measurement uncertainties there is no significant difference in the satellite populations of lens and non lens galaxies.

We present MOSFIRE spectroscopy of 13 candidate z~8 galaxies selected as Y-dropouts as part of the BoRG pure parallel survey. We detect no significant lya emission (our median 1-sigma rest frame equivalent width sensitivity is in the range 2-16 AA). Using the Bayesian framework derived in a previous paper, we perform a rigorous analysis of a statistical subsample of non-detections for ten Y-dropouts, including data from the literature, to study the cosmic evolution of the lya emission of Lyman Break Galaxies. We find that lya emission is suppressed at z~8 by at least a factor of three with respect to z~6 continuing the downward trend found by previous studies of z-dropouts at z~7. This finding suggests a dramatic evolution in the conditions of the intergalactic or circumgalactic media in just 300 Myrs, consistent with the onset of reionization or changes in the physical conditions of the first generations of starforming regions.

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.

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.

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 \gamma' (\rho(r)\propto r^{-\gamma'}) on stellar mass, size, and redshift. We find that two parameters are sufficent to determine \gamma' with less than 6% residual scatter. At fixed redshift, \gamma' depends solely on the surface stellar mass density \partial \gamma'/ \partial \Sigma_*=0.38\pm 0.07, i.e. galaxies with denser stars also have steeper slopes. At fixed M_* and R_{eff}, \gamma' depends on redshift, in the sense that galaxies at a lower redshift have steeper slopes (\partial \gamma' / \partial z = -0.31\pm 0.10). However, the mean redshift evolution of \gamma' for an individual galaxy is consistent with zero d\gamma'/dz=-0.10\pm0.12. This result is obtained by combining our measured dependencies of \gamma' 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.

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.

Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays are a low redshift (z~0-2) probe and thus allow one to break degeneracies in the interpretation of data from higher-redshift probes like the cosmic microwave background in terms of the dark energy equation of state. Current studies are limited by the size of the sample of known lensed quasars, but this situation is about to change. Even in this decade, wide field imaging surveys are likely to discover thousands of lensed quasars, enabling the targeted study of ~100 of these systems and resulting in substantial gains in the dark energy figure of merit. In the next decade, a further order of magnitude improvement will be possible with the 10000 systems expected to be detected and measured with LSST and Euclid. To fully exploit these gains, we identify three priorities. First, support for the development of software required for the analysis of the data. Second, in this decade, small robotic telescopes (1-4m in diameter) dedicated to monitoring of lensed quasars will transform the field by delivering accurate time delays for ~100 systems. Third, in the 2020's, LSST will deliver 1000's of time delays; the bottleneck will instead be the aquisition and analysis of high resolution imaging follow-up. Thus, the top priority for the next decade is to support fast high resolution imaging capabilities, such as those enabled by the James Webb Space Telescope and next generation adaptive optics systems on large ground based telescopes.

The standard cosmological model based on cold dark matter (CDM) predicts a large number of subhalos for each galaxy-size halo. It is well known that matching the subhalos to the observed properties of luminous satellites of galaxies in the local universe poses a significant challenge to our understanding of the astrophysics of galaxy formation. We show that the cosmic evolution and host mass dependence of the luminosity function of satellites provides a powerful new diagnostic to disentangle astrophysical effects from variations in the underlying dark matter mass function. We illustrate this by comparing the results of recent observations of satellites out to $z=0.8$ based on Hubble Space Telescope images with the predictions of three different sets of state-of-the art semi-analytic models with underlying CDM power spectra and one semi-analytic model with an underlying Warm Dark Matter (WDM) power spectrum. We find that even though CDM models provide a reasonable fit to the local luminosity function of satellites around galaxies comparable or slightly larger than the Milky Way, they do not reproduce the data as well for different redshift and host galaxy stellar mass. This tension indicates that further improvements are likely to be needed in the description of star formation if the models are to be reconciled with the data. The WDM model matches the observed mass dependence and redshift evolution of satellite galaxies more closely than any of the CDM models, indicating that a modification of the underlying power spectrum may offer an alternative solution to this tension. We conclude by presenting predictions for the color magnitude relation of satellite galaxies to demonstrate how future observations will be able to further distinguish between these models and help constrain baryonic and non-baryonic physics.

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.

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.

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.

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.

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]

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^{-\beta} on small scales, we find a mean slope <\beta> = 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.

Strong gravitational lenses with measured time delays between the multiple images and models of the lens mass distribution allow a one-step determination of the time-delay distance, and thus a measure of cosmological parameters. We present a blind analysis of the gravitational lens RXJ1131-1231 incorporating (1) the newly measured time delays from COSMOGRAIL, (2) archival HST imaging of the lens system, (3) a new velocity-dispersion measurement of the lens galaxy of 323+/-20km/s based on Keck spectroscopy, and (4) a characterization of the line-of-sight structures via observations of the lens' environment and ray tracing through the Millennium Simulation. Our blind analysis is designed to prevent experimenter bias. The joint analysis of the data sets allows a time-delay distance measurement to 6% precision that takes into account all known systematic uncertainties. In combination with the WMAP7 data set in flat wCDM cosmology, our unblinded cosmological constraints for RXJ1131-1231 are: H_0=80.0+5.8/-5.7km/s/Mpc, OmegaDE=0.79+/-0.03 and w=-1.25+0.17/-0.21. We find the results to be statistically consistent with those from the analysis of the gravitational lens B1608+656. The joint constraints from the two lenses and WMAP7 are H_0=75.2+4.4/-4.2km/s/Mpc, OmegaDE=0.76+0.02/-0.03 and w=-1.14+0.17/-0.20 in flat wCDM, and H_0=73.1+2.4/-3.6km/s/Mpc, OmegaL=0.75+0.01/-0.02 and OmegaK=0.003+0.005/-0.006 in open LCDM. Time-delay lenses constrain especially tightly the Hubble constant (5.7% and 4.0% respectively in wCDM and open LCDM) and curvature of the universe. They complement well other cosmological probes, and provide an independent check of unknown systematics. Our measurement of the Hubble constant is completely independent of those based on the local distance ladder method, providing an important consistency check of the standard cosmological model and of general relativity.

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 $<\gamma_{\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.

Recent work has suggested that the stellar initial mass function (IMF) is not universal, but rather is correlated with galaxy stellar mass, stellar velocity dispersion, or morphological type. In this paper, we investigate variations of the IMF within individual galaxies. For this purpose, we use strong lensing and gas kinematics to measure independently the normalisation of the IMF of the bulge and disk components of a sample of 5 massive spiral galaxies with substantial bulge components taken from the SWELLS survey. We find that the stellar mass of the bulges are tightly constrained by the lensing and kinematic data. A comparison with masses based on stellar population synthesis models fitted to optical and near infrared photometry favors a Salpeter-like normalisation of the IMF. Conversely, the disk masses are less well constrained due to degeneracies with the dark matter halo, but are consistent with Milky Way type IMFs in agreement with previous studies. The disks are submaximal at 2.2 disk scale lengths, but due to the contribution of the bulges, the galaxies are baryon dominated at 2.2 disk scale lengths. Globally, our inferred IMF normalisation is consistent with that found for early-type galaxies of comparable stellar mass (> 10^11 M_sun). Our results suggest a non-universal IMF within the different components of spiral galaxies, adding to the well-known differences in stellar populations between disks and bulges.

We report the discovery of 33 Lyman break galaxy (LBG) candidates at z~8 detected in Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) imaging as part of the Brightest of Reionizing Galaxies (BoRG) pure-parallel survey. The ongoing BoRG survey currently has the largest area (274 arcmin^2) with Y_098 (or Y_105), J_125, and H_160 band coverage needed to search for z~8 galaxies, about three times the current CANDELS area, and slightly larger than what will be the final CANDELS wide component with Y_105 data. Our sample of 33 relatively bright Y_098-dropout galaxies have J_125 band magnitudes between 25.5 and 27.4 mag. This is the largest sample of bright (J_125 <~ 27.4) z~8 galaxy candidates presented to date. Combining our dataset with the Hubble Ultra-Deep Field (HUDF09) dataset, we constrain the rest-frame ultraviolet galaxy luminosity function at z~8 over the widest dynamic range currently available. The combined datasets are well fitted by a Schechter function, i.e. \phi(L) = \phi_* (L/L_*)^{\alpha}\ e^{-(L/L_*)}, without evidence for an excess of sources at the bright end. At 68% confidence, we derive \phi_* = (4.3^{+3.5}_{-2.1}) \times 10^{-4} Mpc^{-3}, M_* = -20.26^{+0.29}_{-0.34}, and a very steep faint-end slope \alpha = -1.98^{+0.23}_{-0.22}. While the best-fit parameters still have a strong degeneracy, especially between \phi_* and M_*, our improved coverage at the bright end has reduced the uncertainty of the faint-end power-law slope at z~8 compared to the best previous determination at +/-0.4. With a future expansion of the BoRG survey, combined with planned ultradeep WFC3/IR observations, it will be possible to further reduce this uncertainty and clearly demonstrate the steepening of the faint-end slope compared to measurements at lower redshift, thereby confirming the key role played by small galaxies in the reionization of the universe.

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.

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 {\lambda} = 9412 {\AA} at >5{\sigma} 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{\alpha} 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{\alpha} flux detected at these redshifts. Given the magnification of {\mu} = 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.

This Resource Letter provides a guide to a selection of the literature on gravitational lensing and its applications. Journal articles, books, popular articles, and websites are cited for the following topics: foundations of gravitational lensing, foundations of cosmology, history of gravitational lensing, strong lensing, weak lensing, and microlensing.

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)-{\sigma}* 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.

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.

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 ($\rho_{\rm{tot}}\propto r^{-\gamma'}$) within the radius of the outer ring to be $\gamma' = 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 $\gamma_{\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.

We infer the normalization and the radial and angular distributions of the number density of satellites of massive galaxies ($\log_{10}[M_{h}^*/M\odot]>10.5$) between redshifts 0.1 and 0.8 as a function of host stellar mass, redshift, morphology and satellite luminosity. Exploiting the depth and resolution of the COSMOS HST images, we detect satellites up to eight magnitudes fainter than the host galaxies and as close as 0.3 (1.4) arcseconds (kpc). Describing the number density profile of satellite galaxies to be a projected power law such that $P(R)\propto R^{\rpower}$, we find $\rpower=-1.1\pm 0.3$. We find no dependency of $\rpower$ on host stellar mass, redshift, morphology or satellite luminosity. Satellites of early-type hosts have angular distributions that are more flattened than the host light profile and are aligned with its major axis. No significant average alignment is detected for satellites of late-type hosts. The number of satellites within a fixed magnitude contrast from a host galaxy is dependent on its stellar mass, with more massive galaxies hosting significantly more satellites. Furthermore, high-mass late-type hosts have significantly fewer satellites than early-type galaxies of the same stellar mass, likely a result of environmental differences. No significant evolution in the number of satellites per host is detected. The cumulative luminosity function of satellites is qualitatively in good agreement with that predicted using subhalo abundance matching techniques. However, there are significant residual discrepancies in the absolute normalization, suggesting that properties other than the host galaxy luminosity or stellar mass determine the number of satellites.

Early-type galaxies (ETGs) are observed to be more compact at z>2 than in the local Universe. Remarkably, much of this size evolution appears to take place in a short (1.8 Gyr) time span between z=2.2 and z=1.3, which poses a serious challenge to hierarchical galaxy formation models where mergers occurring on a similar timescale are the main mechanism for galaxy growth. We compute the merger-driven redshift evolution of stellar mass Mstar\propto(1+z)^aM, half-mass radius Re\propto(1+z)^aR and velocity-dispersion sigma0\propto(1+z)^asigma predicted by concordance Lambda cold dark matter for a typical massive ETG in the redshift range z=1.3-2.2. Neglecting dissipative processes, and thus maximizing evolution in surface density, we find -1.5<aM<-0.6, -1.9<aR<-0.7 and 0.06<asigma<0.22, under the assumption that the accreted satellites are spheroids. It follows that the predicted z=2.2 progenitors of z=1.3 ETGs are significantly less compact (on average a factor of 2 larger Re at given Mstar) than the quiescent galaxies observed at z>2. Furthermore, we find that the scatter introduced in the size-mass correlation by the predicted merger-driven growth is difficult to reconcile with the tightness of the observed scaling law. We conclude that - barring unknown systematics or selection biases in the current measurements - minor and major mergers with spheroids are not sufficient to explain the observed size growth of ETGs within the standard model.

Measuring the relative mass contributions of luminous and dark matter in spiral galaxies is important for understanding their formation and evolution. The combination of a galaxy rotation curve and strong lensing is a powerful way to break the disk-halo degeneracy that is inherent in each of the methods individually. We present an analysis of the 10-image radio spiral lens B1933+503 at z_l=0.755, incorporating (1) new global VLBI observations, (2) new adaptive-optics assisted K-band imaging, (3) new spectroscopic observations for the lens galaxy rotation curve and the source redshift. We construct a three-dimensionally axisymmetric mass distribution with 3 components: an exponential profile for the disk, a point mass for the bulge, and an NFW profile for the halo. The mass model is simultaneously fitted to the kinematics and the lensing data. The NFW halo needs to be oblate with a flattening of a/c=0.33^{+0.07}_{-0.05} to be consistent with the radio data. This suggests that baryons are effective at making the halos oblate near the center. The lensing and kinematics analysis probe the inner ~10 kpc of the galaxy, and we obtain a lower limit on the halo scale radius of 16 kpc (95% CI). The dark matter mass fraction inside a sphere with a radius of 2.2 disk scale lengths is f_{DM,2.2}=0.43^{+0.10}_{-0.09}. The contribution of the disk to the total circular velocity at 2.2 disk scale lengths is 0.76^{+0.05}_{-0.06}, suggesting that the disk is marginally submaximal. The stellar mass of the disk from our modeling is log_{10}(M_{*}/M_{sun}) = 11.06^{+0.09}_{-0.11} assuming that the cold gas contributes ~20% to the total disk mass. In comparison to the stellar masses estimated from stellar population synthesis models, the stellar initial mass function of Chabrier is preferred to that of Salpeter by a probability factor of 7.2.

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]

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.

We report the outcome of a 3-day workshop on the Hubble constant (H_0) that took place during February 6-8 2012 at the Kavli Institute for Particle Astrophysics and Cosmology, on the campus of Stanford University. The participants met to address the following questions. Are there compelling scientific reasons to obtain more precise and more accurate measurements of H_0 than currently available? If there are, how can we achieve this goal? The answers that emerged from the workshop are (1) better measurements of H_0 provide critical independent constraints on dark energy, spatial curvature of the Universe, neutrino physics, and validity of general relativity, (2) a measurement of H_0 to 1% in both precision and accuracy, supported by rigorous error budgets, is within reach for several methods, and (3) multiple paths to independent determinations of H_0 are needed in order to access and control systematics.

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.

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.

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.

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.

Spectroscopic confirmation of galaxies at z~7 and above has been extremely difficult, owing to a drop in intensity of Ly-alpha emission in comparison with samples at z~6. This crucial finding could potentially signal the ending of cosmic reionization. However it is based on small datasets, often incomplete and heterogeneous in nature. We introduce a flexible Bayesian framework, useful to interpret such evidence. Within this framework, we implement two simple phenomenological models: a smooth one, where the distribution of Ly-alpha is attenuated by a factor \es with respect to z~6; a patchy one where a fraction \ep is absorbed/non-emitted while the rest is unabsorbed. From a compilation of 39 observed z~7 galaxies we find \es=0.69+-0.12 and \ep=0.66+-0.16. The models can be used to compute fractions of emitters above any equivalent width W. For W>25\AA, we find X^{25}_{z=7}=0.37+-0.11 (0.14+-0.06) for galaxies fainter (brighter) than M_{UV}=-20.25 for the patchy model, consistent with previous work, but with smaller uncertainties by virtue of our full use of the data. At z~8 we combine new deep (5-\sigma flux limit 10^{-17}ergs^{-1}cm^{-2}) Keck-NIRSPEC observations of a bright Y-dropout identified by our BoRG Survey, with those of three objects from the literature and find that the inference is inconclusive. We compute predictions for future near-infrared spectroscopic surveys and show that it is challenging but feasible to constrain the distribution of Ly-alpha emitters at z~8 and distinguish between models.

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.

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.

Theoretical and numerical modeling of dark-matter halo assembly predicts that the most luminous galaxies at high redshift are surrounded by overdensities of fainter companions. We test this prediction with HST observations acquired by our Brightest of Reionizing Galaxies (BoRG) survey, which identified four very bright z~8 candidates as Y-dropout sources in four of the 23 non-contiguous WFC3 fields observed. We extend here the search for Y-dropouts to fainter luminosities (M_* galaxies with M_AB\sim-20), with detections at >5sigma confidence (compared to >8sigma confidence adopted earlier) identifying 17 new candidates. We demonstrate that there is a correlation between number counts of faint and bright Y-dropouts at >99.84% confidence. Field BoRG58, which contains the best bright z\sim8 candidate (M_AB=-21.3), has the most significant overdensity of faint Y-dropouts. Four new sources are located within 70arcsec (corresponding to 3.1 comoving Mpc at z=8) from the previously known brighter z\sim8 candidate. The overdensity of Y-dropouts in this field has a physical origin to high confidence (p>99.975%), independent of completeness and contamination rate of the Y-dropout selection. We modeled the overdensity by means of cosmological simulations and estimate that the principal dark matter halo has mass M_h\sim(4-7)x10^11Msun (\sim5sigma density peak) and is surrounded by several M_h\sim10^11Msun halos which could host the fainter dropouts. In this scenario, we predict that all halos will eventually merge into a M_h>2x10^14Msun galaxy cluster by z=0. Follow-up observations with ground and space based telescopes are required to secure the z\sim8 nature of the overdensity, discover new members, and measure their precise redshift.

The Lick AGN Monitoring Project 2011 observing campaign was carried out over the course of 11 weeks in Spring 2011. Here we present the first results from this program, a measurement of the broad-line reverberation lag in the Seyfert 1 galaxy Mrk 50. Combining our data with supplemental observations obtained prior to the start of the main observing campaign, our dataset covers a total duration of 4.5 months. During this time, Mrk 50 was highly variable, exhibiting a maximum variability amplitude of a factor of 4 in the U-band continuum and a factor of 2 in the H-beta line. Using standard cross-correlation techniques, we find that H-beta and H-gamma lag the V-band continuum by tau_cen = 10.64(-0.93,+0.82) and 8.43(-1.28,+1.30) days, respectively, while the lag of He II 4686 is unresolved. The H-beta line exhibits a symmetric velocity-resolved reverberation signature with shorter lags in the high-velocity wings than in the line core, consistent with an origin in a broad-line region dominated by orbital motion rather than infall or outflow. Assuming a virial normalization factor of f=5.25, the virial estimate of the black hole mass is (3.2+-0.5)*10^7 solar masses. These observations demonstrate that Mrk 50 is among the most promising nearby active galaxies for detailed investigations of broad-line region structure and dynamics.

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.

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.

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.

We present the first results from the X-shooter Lens Survey (XLENS): an analysis of the massive early-type galaxy SDSS J1148+1930 at redshift z=0.444. We combine its extended kinematic profile -derived from spectra obtained with X-shooter on the ESO VLT- with strong gravitational lensing and multi-color information derived from SDSS images. Our main results are (i) The luminosity-weighted stellar velocity dispersion is \sigma *(<Reff) = 351\pm10 km/s, extracted from a rectangular aperture of 1.8"x1.6", more accurate and considerably lower than a previously published value of \sim450km/s. (ii) A single-component (stellar plus dark) mass model of the lens galaxy yields a logarithmic total-density slope of {\gamma}' = 1.72+0.05-0.06. (iii) The projected stellar mass fraction, derived solely from lensing, is f*(<Reff) = 0.19+0.04-0.09 inside the Einstein radius for a Hernquist profile and no anisotropy. The dark-matter fraction inside the effective radius f_DM(<Reff) = 0.60+0.15\pm0.1. (iv) Based on the SDSS colors, we find f*(<Reff)=0.17 \pm 0.06 for a Salpeter IMF and f*(<Reff)=0.07 \pm 0.02 for a Chabrier IMF. The lensing and dynamics constraints on the stellar mass fraction agree well with those independently derived from the SDSS colors for a Salpeter IMF, preferred over Chabrier IMF at variance with standard results for lower mass galaxies. Dwarf-rich IMFs in the lower mass range of 0.1-0.7 solar mass, with {\alpha} >= 3 (with dN/dM \propto M^(-{\alpha})) -such as that recently suggested for massive ETGs- are excluded at the > 90 % C.L. and in some cases violate the total lensing-derived mass limit. We conclude that this very massive ETG is dark-matter dominated inside one effective radius, consistent with the trend recently found from massive SLACS galaxies, with a IMF normalization consistent with Salpeter.

The relative contribution of baryons and dark matter to the inner regions of spiral galaxies provides critical clues to their formation and evolution, but it is generally difficult to determine. For spiral galaxies that are strong gravitational lenses, however, the combination of lensing and kinematic observations can be used to break the disk-halo degeneracy. In turn, such data constrain fundamental parameters such as i) the mass density profile slope and axis ratio of the dark matter halo, and by comparison with dark matter-only numerical simulations the modifications imposed by baryons; ii) the mass in stars and therefore the overall star formation efficiency, and the amount of feedback; iii) by comparison with stellar population synthesis models, the normalization of the stellar initial mass function. In this first paper of a series, we present a sample of 16 secure, 1 probable, and 6 possible strong lensing spiral galaxies, for which multi-band high-resolution images and rotation curves were obtained using the Hubble Space Telescope and Keck-II Telescope as part of the Sloan WFC Edge-on Late-type Lens Survey (SWELLS). The sample includes 8 newly discovered secure systems. [abridged] We find that the SWELLS sample of secure lenses spans a broad range of morphologies (from lenticular to late-type spiral), spectral types (quantified by Halpha emission), and bulge to total stellar mass ratio (0.22-0.85), while being limited to M_*>10^{10.5} M_sun. The SWELLS sample is thus well-suited for exploring the relationship between dark and luminous matter in a broad range of galaxies. We find that the deflector galaxies obey the same size-mass relation as that of a comparison sample of elongated non-lens galaxies selected from the SDSS survey. We conclude that the SWELLS sample is consistent with being representative of the overall population of high-mass high-inclination disky galaxies.

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.

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