The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

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.

The Wide Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) will be the first high-performance stellar coronagraph using active wavefront control for deep starlight suppression in space, providing unprecedented levels of contrast, spatial resolution, and sensitivity for astronomical observations in the optical. One science case enabled by the CGI will be taking images and(R~50)spectra of faint interplanetary dust structures present in the habitable zone of nearby sunlike stars (~10 pc) and within the snow-line of more distant ones(~20pc), down to dust density levels commensurate with that of the solar system zodiacal cloud. Reaching contrast levels below~10-7 for the first time, CGI will cross an important threshold in debris disks physics, accessing disks with low enough optical depths that their structure is dominated by transport phenomena than collisions. Hence, CGI results will be crucial for determining how exozodiacal dust grains are produced and transported in low-density disks around mature stars. Additionally, CGI will be able to measure the brightness level and constrain the degree of asymmetry of exozodiacal clouds around individual nearby sunlike stars in the optical, at the ~10x solar zodiacal emission level. This information will be extremely valuable for optimizing the observational strategy of possible future exo-Earth direct imaging missions, especially those planning to operate at optical wavelengths, such as Habitable Exoplanet Observatory (HabEx) and the Large Ultraviolet/Optical/Infrared Surveyor (LUVOIR).

In the cold dark matter (CDM) picture of structure formation, galaxy mass distributions are predicted to have a considerable amount of structure on small scales. Strong gravitational lensing has proven to be a useful tool for studying this small-scale structure. Much of the attention has been given to detecting individual dark matter subhalos through lens modeling, but recent work has suggested that the full population of subhalos could be probed using a power spectrum analysis. In this paper we quantify the power spectrum of small-scale structure in simulated galaxies, with the goal of understanding theoretical predictions and setting the stage for using measurements of the power spectrum to test dark matter models. We use a sample of simulated galaxies generated from the \texttt{Galacticus} semi-analytic model to determine the power spectrum distribution first in the CDM paradigm and then in a warm dark matter scenario. We find that a measurement of the slope and amplitude of the power spectrum on galaxy strong lensing scales ($k\sim 1$ kpc$^{-1}$) could be used to distinguish between CDM and alternate dark matter models, especially if the most massive subhalos can be directly detected via gravitational imaging.

Gravitational lensing has emerged as a powerful probe of the matter distribution on subgalactic scales, which itself may contain important clues about the fundamental origins and properties of dark matter. Broadly speaking, two different approaches have been taken in the literature to map the small-scale structure of the Universe using strong lensing, with one focused on measuring the position and mass of a small number of discrete massive subhalos appearing close in projection to lensed images, and the other focused on detecting the collective effect of all the small-scale structure between the lensed source and the observer. In this paper, we follow the latter approach and perform a detailed study of the sensitivity of galaxy-scale gravitational lenses to the ensemble properties of small-scale structure. As in some previous studies, we adopt the language of the substructure power spectrum to characterize the statistical properties of the small-scale density field. We present a comprehensive theory that treats lenses with extended sources as well as those with time-dependent compact sources (such as quasars) in a unified framework for the first time. Our approach uses mode functions to provide both computational advantages and insights about couplings between the lens and source. The goal of this paper is to develop the theory and gain the intuition necessary to understand how the sensitivity to the substructure power spectrum depends on the source and lens properties, with the eventual aim of identifying the most promising targets for such studies.

Microlensing by stars within distant galaxies acting as strong gravitational lenses of multiply-imaged quasars, provides a unique and direct measurement of the internal structure of the lensed quasar on nano-arcsecond scales. The measurement relies on the temporal variation of high-magnification caustic crossings which vary on timescales of days to years. Multiwavelength observations provide information from distinct emission regions in the quasar. Through monitoring of these strong gravitational lenses, a full tomographic view can emerge with Astronomical-Unit scale resolution. Work to date has demonstrated the potential of this technique in about a dozen systems. In the 2020s there will be orders of magnitude more systems to work with. Monitoring of lens systems for caustic-crossing events to enable triggering of multi-platform, multi-wavelength observations in the 2020s will fulfill the potential of quasar microlensing as a unique and comprehensive probe of active black hole structure and dynamics.

Optical and X-ray observations of strongly gravitationally lensed quasars (especially when four separate images of the quasar are produced) determine not only the amount of matter in the lensing galaxy but also how much is in a smooth component and how much is composed of compact masses (e.g., stars, stellar remnants, primordial black holes, CDM sub-halos, and planets). Future optical surveys will discover hundreds to thousands of quadruply lensed quasars, and sensitive X-ray observations will unambiguously determine the ratio of smooth to clumpy matter at specific locations in the lensing galaxies and calibrate the stellar mass fundamental plane, providing a determination of the stellar $M/L$. A modest observing program with a sensitive, sub-arcsecond X-ray imager, combined with the planned optical observations, can make those determinations for a large number (hundreds) of the lensing galaxies, which will span a redshift range of $\sim$$0.25<z<1.5$

We study halo mass functions with the Press-Schechter formalism for interacting dark matter models, where matter power spectra are damped due to dark acoustic oscillations in the early universe. After adopting a smooth window function, we calibrate the analytical model with numerical simulations from the "effective theory of structure formation" (ETHOS) project and fix the model parameters in the high mass regime, $M_{\rm h}\gtrsim3\times10^{10}\;{\rm M}_{\odot}$. We also perform high-resolution cosmological simulations with halo masses down to $M_{\rm h}\sim10^8\;{\rm M}_{\odot}$ to cover a wide mass range for comparison. Although the model is calibrated with ETHOS1 and CDM simulations for high halo masses at redshift $z=0$, it successfully reproduces simulations for two other ETHOS models in the low mass regime at low and high redshifts. As an application, we compare the cumulative number density of haloes to that of observed galaxies at $z=6$, and find the interacting dark matter models with a kinetic decoupling temperature below $0.5\ \rm{keV}$ is disfavored. We also perform the abundance-matching analysis and derive the stellar-halo mass relation for these models at $z=4$. Suppression in halo abundance leads to less massive haloes that host observed galaxies in the stellar mass range $M_*\simeq 10^5-10^7\ {\rm M}_{\odot}$.

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

The Wide Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) is a high-contrast imager and integral field spectrograph that will enable the study of exoplanets and circumstellar disks at visible wavelengths. Ground-based high-contrast instrumentation has fundamentally limited performance at small working angles, even under optimistic assumptions for 30m-class telescopes. There is a strong scientific driver for better performance, particularly at visible wavelengths. Future flagship mission concepts aim to image Earth analogues with visible light flux ratios of more than 10^10. CGI is a critical intermediate step toward that goal, with a predicted 10^8-9 flux ratio capability in the visible. CGI achieves this through improvements over current ground and space systems in several areas: (i) Hardware: space-qualified (TRL9) deformable mirrors, detectors, and coronagraphs, (ii) Algorithms: wavefront sensing and control; post-processing of integral field spectrograph, polarimetric, and extended object data, and (iii) Validation of telescope and instrument models at high accuracy and precision. This white paper, submitted to the 2018 NAS Exoplanet Science Strategy call, describes the status of key CGI technologies and presents ways in which performance is likely to evolve as the CGI design matures.

The Wide-Field Infrared Survey Telescope (WFIRST) is expected to launch in the mid-2020s. With its wide-field near-infrared (NIR) camera, it will survey the sky to unprecedented detail. As part of normal operations and as the result of multiple expected dedicated surveys, WFIRST will produce several relatively wide-field (tens of square degrees) deep (limiting magnitude of 28 or fainter) fields. In particular, a planned supernova survey is expected to image 3 deep fields in the LSST footprint roughly every 5 days over 2 years. Stacking all data, this survey will produce, over all WFIRST supernova fields in the LSST footprint, ~12-25 deg^2 and ~5-15 deg^2 regions to depths of ~28 mag and ~29 mag, respectively. We suggest LSST undertake mini-surveys that will match the WFIRST cadence and simultaneously observe the supernova survey fields during the 2-year WFIRST supernova survey, achieving a stacked depth similar to that of the WFIRST data. We also suggest additional observations of these same regions throughout the LSST survey to get deep images earlier, have long-term monitoring in the fields, and produce deeper images overall. These fields will provide a legacy for cosmology, extragalactic, and transient/variable science.

The Wide Field Infrared Survey Telescope (WFIRST) was the top ranked large space mission in the 2010 New Worlds, New Horizons decadal survey, and it was formed by merging the science programs of 3 different mission concepts, including the Microlensing Planet Finder (MPF) concept (Bennett \etal\ 2010). The WFIRST science program (Spergel \etal\ 2015) consists of a general observer program, a wavefront controlled technology program, and two targeted science programs: a program to study dark energy, and a statistical census of exoplanets with a microlensing survey, which uses nearly one quarter of WFIRST's observing time in the current design reference mission. The New Worlds, New Horizons (decadal survey) midterm assessment summarizes the science case for the WFIRST exoplanet microlensing survey with this statement: "WFIRST's microlensing census of planets beyond 1 AU will perfectly complement Kepler's census of compact systems, and WFIRST will also be able to detect free-floating planets unbound from their parent stars\rlap."

We have selected a sample of eleven massive clusters of galaxies observed by the Hubble Space Telescope in order to study the impact of the dynamical state on the IntraCluster Light (ICL) fraction, the ratio of total integrated ICL to the total galaxy member light. With the exception of the Bullet cluster, the sample is drawn from the Cluster Lensing and Supernova Survey and the Frontier Fields program, containing five relaxed and six merging clusters. The ICL fraction is calculated in three optical filters using the CHEFs IntraCluster Light Estimator, a robust and accurate algorithm free of a priori assumptions. We find that the ICL fraction in the three bands is, on average, higher for the merging clusters, ranging between $\sim7-23\%$, compared with the $\sim 2-11\%$ found for the relaxed systems. We observe a nearly constant value (within the error bars) in the ICL fraction of the regular clusters at the three wavelengths considered, which would indicate that the colors of the ICL and the cluster galaxies are, on average, coincident and, thus, their stellar populations. However, we find a higher ICL fraction in the F606W filter for the merging clusters, consistent with an excess of lower-metallicity/younger stars in the ICL, which could have migrated violently from the outskirts of the infalling galaxies during the merger event.

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 VLA images and HST/STIS spectra of sources within the center of the brightest cluster galaxy (BCG) in Abell 2261. These observations were obtained to test the hypothesis that its extremely large, flat core reflects the ejection of its supermassive black hole. Spectra of three of the four most luminous "knots" embedded in the core were taken to test whether one may represent stars bound to a displaced massive black hole. The three knots have radial velocity offsets dV < ~150 km/s from the BCG. Knots 2 and 3 show kinematics, colors, and stellar masses consistent with infalling low-mass galaxies or larger stripped cluster members. Large errors in the stellar velocity dispersion of Knot 1, however, mean that we cannot rule out the hypothesis that it hosts a high-mass black hole. A2261-BCG has a compact, relic radio-source offset by 6.5 kpc (projected) from the optical core's center, but no active radio core that would pinpoint the galaxy's central black hole to a tight 10 GHz flux limit <3.6 uJy. Its spectrum and morphology are suggestive of an active galactic nucleus that switched off >48 Myr ago, with an equipartition condition magnetic field of 15 uG. These observations are still consistent with the hypothesis that the nuclear black hole has been ejected from its core, but the critical task of locating the supermassive black hole or demonstrating that A2261-BCG lacks one remains to be done.

The Large Ultraviolet / Optical / Infrared Surveyor (LUVOIR) is one of four large mission concepts currently undergoing community study for consideration by the 2020 Astronomy and Astrophysics Decadal Survey. The LUVOIR Ultraviolet Multi-Object Spectrograph, LUMOS, is being designed to support all of the UV science requirements of LUVOIR, from exoplanet host star characterization to tomography of circumgalactic halos to water plumes on outer solar system satellites. LUMOS offers point source and multi-object spectroscopy across the UV bandpass, with multiple resolution modes to support different science goals. The instrument will provide low (R = 8,000-18,000) and medium (R = 30,000-65,000) resolution modes across the far-ultraviolet (FUV: 100-200 nm) and near-ultraviolet (NUV: 200-400 nm) windows, and a very low resolution mode (R = 500) for spectroscopic investigations of extremely faint objects in the FUV. Imaging spectroscopy will be accomplished over a 3 x 1.6 arcminute field-of-view by employing holographically-ruled diffraction gratings to control optical aberrations, microshutter arrays (MSA), advanced optical coatings for high-throughput in the FUV, and next generation large-format photon-counting detectors. The spectroscopic capabilities of LUMOS are augmented by an FUV imaging channel (100-200nm, 13 milliarcsecond angular resolution, 2 x 2 arcminute field-of-view) that will employ a complement of narrow and medium-band filters. We present an overview of LUMOS' observing modes and estimated performance curves for effective area, spectral resolution, and imaging performance. Example "LUMOS 100-hour Highlights" observing programs are presented to demonstrate the potential power of LUVOIR's ultraviolet spectroscopic capabilities.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

We present HST WFC3 F160W imaging and infrared spectral energy distributions for twelve extremely luminous, obscured AGN at $1.8<z<2.7$, selected via "Hot, Dust Obscured" mid-infrared colors. Their infrared luminosities span $2-15\times10^{13}$L$_{\odot}$, making them among the most luminous objects in the Universe at $z\sim2$. In all cases the infrared emission is consistent with arising at least in most part from AGN activity. The AGN fractional luminosities are higher than those in either sub-millimeter galaxies, or AGN selected via other mid-infrared criteria. Adopting the $G$, M$_{20}$ and $A$ morphological parameters, together with traditional classification boundaries, infers that three quarters of the sample as mergers. Our sample do not, however, show any correlation between the considered morphological parameters and either infrared luminosity or AGN fractional luminosity. Moreover, their asymmetries and effective radii are distributed identically to those of massive galaxies at $z\sim2$. We conclude that our sample is not preferentially associated with mergers, though a significant merger fraction is still plausible. Instead, we propose that our sample are examples of the massive galaxy population at $z\sim2$ that harbor a briefly luminous, "flickering" AGN, and in which the $G$ and M$_{20}$ values have been perturbed, due to either the AGN, and/or the earliest formation stages of a bulge in an inside-out manner. Furthermore, we find that the mass assembly of the central black holes in our sample leads the mass assembly of any bulge component. Finally, we speculate that our sample represent a small fraction of the immediate antecedents of compact star-forming galaxies at $z\sim2$.

A conclusive model for the formation of dwarf spheroidal (dSph) galaxies still remains elusive. Owing to their proximity to the massive spirals Milky Way (MW) and M31, various environmental processes have been invoked to explain their origin. In this context, the tidal stirring model postulates that interactions with MW-sized hosts can transform rotationally supported dwarfs, resembling present-day dwarf irregular (dIrr) galaxies, into systems with the kinematic and structural properties of dSphs. Using N-body+SPH simulations, we investigate the dependence of this transformation mechanism on the gas fraction, fgas, in the disk of the progenitor dwarf. Our numerical experiments incorporate for the first time the combined effects of radiative cooling, ram-pressure stripping, star formation, supernova (SN) winds, and a cosmic UV background. For a given orbit inside the primary galaxy, rotationally supported dwarfs with gas fractions akin to those of observed dIrrs (fgas >= 0.5), demonstrate a substantially enhanced likelihood and efficiency of transformation into dSphs relative to their collisionless (fgas = 0) counterparts. We argue that the combination of ram-pressure stripping and SN winds causes the gas-rich dwarfs to respond more impulsively to tides, augmenting their transformation. When fgas >= 0.5, disky dwarfs on previously unfavorable low-eccentricity or large-pericenter orbits are still able to transform. On the widest orbits, the transformation is incomplete; the dwarfs retain significant rotational support, a relatively flat shape, and some gas, naturally resembling transition-type systems. We conclude that tidal stirring constitutes a prevalent evolutionary mechanism for shaping the structure of dwarf galaxies within the currently favored CDM cosmological paradigm.

We present a joint analysis of Chandra X-ray observations, Bolocam thermal Sunyaev-Zel'dovich (SZ) effect observations, Hubble Space Telescope (HST) strong lensing data, and HST and Subaru Suprime-Cam weak lensing data. The multiwavelength dataset is used to constrain parametric models for the distribution of dark and baryonic matter in a sample of six massive galaxy clusters selected from the Cluster Lensing And Supernova survey with Hubble (CLASH). For five of the six clusters, the multiwavelength dataset is well described by a relatively simple model that assumes spherical symmetry, hydrostatic equilibrium, and entirely thermal pressure support. The joint analysis yields considerably better constraints on the total mass and concentration of the cluster compared to analysis of any one dataset individually. The subsample of five galaxy clusters is used to place an upper limit on the fraction of pressure support in the intracluster medium (ICM) due to nonthermal processes, such as turbulence and bulk flow of the gas. We constrain the nonthermal pressure fraction at r500c to be less than 0.11 at 95 percent confidence. This is in tension with state-of-the-art hydrodynamical simulations, which predict a nonthermal pressure fraction of approximately 0.25 at r500c for clusters of similar mass and redshift. This tension may be explained by the sample selection and/or our assumption of spherical symmetry.

In the standard structure formation scenario based on the cold dark matter paradigm, galactic halos are predicted to contain a large population of dark matter subhalos. While the most massive members of the subhalo population can appear as luminous satellites and be detected in optical surveys, establishing the existence of the low mass and mostly dark subhalos has proven to be a daunting task. Galaxy-scale strong gravitational lenses have been successfully used to study mass substructures lying close to lensed images of bright background sources. However, in typical galaxy-scale lenses, the strong lensing region only covers a small projected area of the lens's dark matter halo, implying that the vast majority of subhalos cannot be directly detected in lensing observations. In this paper, we point out that this large population of dark satellites can collectively affect gravitational lensing observables, hence possibly allowing their statistical detection. Focusing on the region of the galactic halo outside the strong lensing area, we compute from first principles the statistical properties of perturbations to the gravitational time delay and position of lensed images in the presence of a mass substructure population. We find that in the standard cosmological scenario, the statistics of these lensing observables are well approximated by Gaussian distributions. The formalism developed as part of this calculation is very general and can be applied to any halo geometry and choice of subhalo mass function. Our results significantly reduce the computational cost of including a large substructure population in lens models and enable the use of Bayesian inference techniques to detect and characterize the distributed satellite population of distant lens galaxies.

We report the discovery of significant mass/light offsets in the strong gravitational lensing system SDSS\,J1011$+$0143. We use the high-resolution \textsl{Hubble Space Telescope} (\textsl{HST}) F555W- and F814W-band imaging and Sloan Digital Sky Survey (SDSS) spectroscopy of this system, which consists of a close galaxy pair with a projected separation of $\approx 4.2$ kpc at $z_{\rm lens} \sim 0.331$ lensing a Ly$\alpha$ emitter (LAE) at $z_{\rm source} = 2.701$. Comparisons between the mass peaks inferred from lens models and light peaks from \textsl{HST} imaging data reveal significant spatial mass/light offsets as large as $1.72 \pm 0.24 \pm 0.34$ kpc in both filter bands. Such large mass/light offsets, not seen in isolated field lens galaxies and relaxed galaxy groups, may be related to the interactions between the two lens galaxies. The detected mass/light offsets can potentially serve as an important test for the self-interacting dark matter model. However, other mechanisms such as dynamical friction on spatially differently distributed dark matter and stars could produce similar offsets. Detailed hydrodynamical simulations of galaxy-galaxy interactions with self-interacting dark matter could accurately quantify the effects of different mechanisms. The background LAE is found to contain three distinct star-forming knots with characteristic sizes from 116 pc to 438 pc. It highlights the power of strong gravitational lensing in probing the otherwise too faint and unresolved structures of distance objects below subkiloparsec or even 100 pc scales through its magnification effect.

We present scaling relations between the integrated Sunyaev-Zel'dovich Effect (SZE) signal, $Y_{\rm SZ}$, its X-ray analogue, $Y_{\rm X}\equiv M_{\rm gas}T_{\rm X}$, and total mass, $M_{\rm tot}$, for the 45 galaxy clusters in the Bolocam X-ray-SZ (BOXSZ) sample. All parameters are integrated within $r_{2500}$. $Y_{2500}$ values are measured using SZE data collected with Bolocam, operating at 140 GHz at the Caltech Submillimeter Observatory (CSO). The temperature, $T_{\rm X}$, and mass, $M_{\rm gas,2500}$, of the intracluster medium are determined using X-ray data collected with Chandra, and $M_{\rm tot}$ is derived from $M_{\rm gas}$ assuming a constant gas mass fraction. Our analysis accounts for several potential sources of bias, including: selection effects, contamination from radio point sources, and the loss of SZE signal due to noise filtering and beam-smoothing effects. We measure the $Y_{2500}$--$Y_{\rm X}$ scaling to have a power-law index of $0.84\pm0.07$, and a fractional intrinsic scatter in $Y_{2500}$ of $(21\pm7)\%$ at fixed $Y_{\rm X}$, both of which are consistent with previous analyses. We also measure the scaling between $Y_{2500}$ and $M_{2500}$, finding a power-law index of $1.06\pm0.12$ and a fractional intrinsic scatter in $Y_{2500}$ at fixed mass of $(25\pm9)\%$. While recent SZE scaling relations using X-ray mass proxies have found power-law indices consistent with the self-similar prediction of 5/3, our measurement stands apart by differing from the self-similar prediction by approximately 5$\sigma$. Given the good agreement between the measured $Y_{2500}$--$Y_{\rm X}$ scalings, much of this discrepancy appears to be caused by differences in the calibration of the X-ray mass proxies adopted for each particular analysis.

We present profiles of temperature (Tx), gas mass, and hydrostatic mass estimated from new and archival X-ray observations of CLASH clusters. We compare measurements derived from XMM and Chandra observations with one another and compare both to gravitational lensing mass profiles derived with CLASH HST and ground-based lensing data. Radial profiles of Chandra and XMM electron density and enclosed gas mass are nearly identical, indicating that differences in hydrostatic masses inferred from X-ray observations arise from differences in Tx measurements. Encouragingly, cluster Txs are consistent with one another at ~100-200 kpc radii but XMM Tx systematically decline relative to Chandra Tx at larger radii. The angular dependence of the discrepancy suggests additional investigation on systematics such as the XMM point spread function correction, vignetting and off-axis responses. We present the CLASH-X mass-profile comparisons in the form of cosmology-independent and redshift-independent circular-velocity profiles. Ratios of Chandra HSE mass profiles to CLASH lensing profiles show no obvious radial dependence in the 0.3-0.8 Mpc range. However, the mean mass biases inferred from the WL and SaWLens data are different. e.g., the weighted-mean value at 0.5 Mpc is <b> = 0.12 for the WL comparison and <b> = -0.11 for the SaWLens comparison. The ratios of XMM HSE mass profiles to CLASH lensing profiles show a pronounced radial dependence in the 0.3-1.0 Mpc range, with a weighted mean mass bias of value rising to <b>~0.3 at ~1 Mpc for the WL comparison and <b> of 0.25 for SaWLens comparison. The enclosed gas mass profiles from both Chandra and XMM rise to a value 1/8 times the total-mass profiles inferred from lensing at 0.5 Mpc and remain constant outside of that radius, suggesting that [8xMgas] profiles may be an excellent proxy for total-mass profiles at >0.5 Mpc in massive galaxy clusters.

We present 20 WISE-selected galaxies with bolometric luminosities L_bol > 10^14 L_sun, including five with infrared luminosities L_IR = L(rest 8-1000 micron) > 10^14 L_sun. These "extremely luminous infrared galaxies," or ELIRGs, were discovered using the "W1W2-dropout" selection criteria which requires marginal or non-detections at 3.4 and 4.6 micron (W1 and W2, respectively) but strong detections at 12 and 22 micron in the WISE survey. Their spectral energy distributions are dominated by emission at rest-frame 4-10 micron, suggesting that hot dust with T_d ~ 450K is responsible for the high luminosities. These galaxies are likely powered by highly obscured AGNs, and there is no evidence suggesting these systems are beamed or lensed. We compare this WISE-selected sample with 116 optically selected quasars that reach the same L_bol level, corresponding to the most luminous unobscured quasars in the literature. We find that the rest-frame 5.8 and 7.8 micron luminosities of the WISE-selected ELIRGs can be 30-80% higher than that of the unobscured quasars. The existence of AGNs with L_bol > 10^14 L_sun at z > 3 suggests that these supermassive black holes are born with large mass, or have very rapid mass assembly. For black hole seed masses ~ 10^3 M_sun, either sustained super-Eddington accretion is needed, or the radiative efficiency must be <15%, implying a black hole with slow spin, possibly due to chaotic accretion.

We report on our first set of spectroscopic Hubble Space Telescope observations of the z~11 candidate galaxy strongly lensed by the MACSJ0647.7+7015 galaxy cluster. The three lensed images are faint and we show that these early slitless grism observations are of sufficient depth to investigate whether this high-redshift candidate, identified by its strong photometric break at ~1.5 micron, could possibly be an emission line galaxy at a much lower redshift. While such an interloper would imply the existence of a rather peculiar object, we show here that such strong emission lines would clearly have been detected. Comparing realistic, two-dimensional simulations to these new observations we would expect the necessary emission lines to be detected at >5 sigma while we see no evidence for such lines in the dispersed data of any of the three lensed images. We therefore exclude that this object could be a low redshift emission line interloper, which significantly increases the likelihood of this candidate being a bona fide z~11 galaxy.

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 utilize the CLASH (Cluster Lensing And Supernova survey with Hubble) observations of 25 clusters to search for extreme emission-line galaxies (EELGs). The selections are carried out in two central bands: F105W (Y105) and F125W (J125), as the flux of the central bands could be enhanced by the presence of [O III] 4959, 5007 at redshift of about 0.93-1.14 and 1.57-1.79, respectively. The multi-band observations help to constrain the equivalent widths of emission lines. Thanks to cluster lensing, we are able to identify 52 candidates down to an intrinsic limiting magnitude of 28.5 and to a rest-frame [O III] 4959,5007 equivalent width of about 3737 angstrom. Our samples include a number of EELGs at lower luminosities that are missed in other surveys, and the extremely high equivalent width can be only found in such faint galaxies. These EELGs can mimic the dropout feature similar to that of high redshift galaxies and contaminate the color-color selection of high redshift galaxies when the S/N ratio is limited or the band coverage is incomplete. We predict that the fraction of EELGs in the future high redshift galaxy selections cannot be neglected.

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 utilise a two-color Lyman-Break selection criterion to search for z~9-10 galaxies over the first 19 clusters in the CLASH program. A systematic search yields three z~9-10 candidates. While we have already reported the most robust of these candidates, MACS1149-JD, two additional z~9 candidates are also found and have H_{160}-band magnitudes of ~26.2-26.9. A careful assessment of various sources of contamination suggests <~1 contaminants for our z~9-10 selection. To determine the implications of these search results for the LF and SFR density at z~9, we introduce a new differential approach to deriving these quantities in lensing fields. Our procedure is to derive the evolution by comparing the number of z~9-10 galaxy candidates found in CLASH with the number of galaxies in a slightly lower redshift sample (after correcting for the differences in selection volumes), here taken to be z~8. This procedure takes advantage of the fact that the relative volumes available for the z~8 and z~9-10 selections behind lensing clusters are not greatly dependent on the details of the lensing models. We find that the normalization of the UV LF at z~9 is just 0.28_{-0.20}^{+0.39}\times that at z~8, ~1.4_{-0.8}^{+3.0}x lower than extrapolating z~4-8 LF results. While consistent with the evolution in the UV LF seen at z~4-8, these results marginally favor a more rapid evolution at z>8. Compared to similar evolutionary findings from the HUDF, our result is less insensitive to large-scale structure uncertainties, given our many independent sightlines on the high-redshift universe.

We utilize 16 band Hubble Space Telescope (HST) observations of 18 lensing clusters obtained as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) Multi-Cycle Treasury program to search for $z\sim6-8$ galaxies. We report the discovery of 204, 45, and 13 Lyman-break galaxy candidates at $z\sim6$, $z\sim7$, and $z\sim8$, respectively, identified from purely photometric redshift selections. This large sample, representing nearly an order of magnitude increase in the number of magnified star-forming galaxies at $z\sim 6-8$ presented to date, is unique in that we have observations in four WFC3/UVIS UV, seven ACS/WFC optical, and all five WFC3/IR broadband filters, which enable very accurate photometric redshift selections. We construct detailed lensing models for 17 of the 18 clusters to estimate object magnifications and to identify two new multiply lensed $z \gtrsim 6$ candidates. The median magnifications over the 17 clusters are 4, 4, and 5 for the $z\sim6$, $z\sim7$, and $z\sim8$ samples, respectively, over an average area of 4.5 arcmin$^2$ per cluster. We compare our observed number counts with expectations based on convolving "blank" field UV luminosity functions through our cluster lens models and find rough agreement down to $\sim27$ mag, where we begin to suffer significant incompleteness. In all three redshift bins, we find a higher number density at brighter observed magnitudes than the field predictions, empirically demonstrating for the first time the enhanced efficiency of lensing clusters over field surveys. Our number counts also are in general agreement with the lensed expectations from the cluster models, especially at $z\sim6$, where we have the best statistics.

We describe the methodology to include nonlinear evolution, including tidal effects, in the computation of subhalo distribution properties in both cold (CDM) and warm (WDM) dark matter universes. Using semi-analytic modeling, we include effects from dynamical friction, tidal stripping, and tidal heating, allowing us to dynamically evolve the subhalo distribution. We calibrate our nonlinear evolution scheme to the CDM subhalo mass function in the Aquarius N-body simulation, producing a subhalo mass function within the range of simulations. We find tidal effects to be the dominant mechanism of nonlinear evolution in the subhalo population. Finally, we compute the subhalo mass function for $m_\chi=1.5$ keV WDM including the effects of nonlinear evolution, and compare radial number densities and mass density profiles of subhalos in CDM and WDM models. We show that all three signatures differ between the two dark matter models, suggesting that probes of substructure may be able to differentiate between them.

WISE has discovered an extraordinary population of hyper-luminous dusty galaxies which are faint in the two bluer passbands ($3.4\, \mu$m and $4.6\, \mu$m) but are bright in the two redder passbands of WISE ($12\, \mu$m and $22\, \mu$m). We report on initial follow-up observations of three of these hot, dust-obscured galaxies, or Hot DOGs, using the CARMA and SMA interferometer arrays at submm/mm wavelengths. We report continuum detections at $\sim$ 1.3 mm of two sources (WISE J014946.17+235014.5 and WISE J223810.20+265319.7, hereafter W0149+2350 and W2238+2653, respectively), and upper limits to CO line emission at 3 mm in the observed frame for two sources (W0149+2350 and WISE J181417.29+341224.8, hereafter W1814+3412). The 1.3 mm continuum images have a resolution of 1-2 arcsec and are consistent with single point sources. We estimate the masses of cold dust are 2.0$\times 10^{8} M_{\odot}$ for W0149+2350 and 3.9$\times 10^{8} M_{\odot}$ for W2238+2653, comparable to cold dust masses of luminous quasars. We obtain 2$\sigma$ upper limits to the molecular gas masses traced by CO, which are 3.3$\times 10^{10} M_{\odot}$ and 2.3$\times 10^{10} M_{\odot}$ for W0149+2350 and W1814+3412, respectively. We also present high-resolution, near-IR imaging with WFC3 on the Hubble Space Telescope for W0149+2653 and with NIRC2 on Keck for W2238+2653. The near-IR images show morphological structure dominated by a single, centrally condensed source with effective radius less than 4 kpc. No signs of gravitational lensing are evident.

We present the results of a numerical study based on the analysis of the MUSIC-2 simulations, aimed at estimating the expected concentration-mass relation for the CLASH cluster sample. We study nearly 1400 halos simulated at high spatial and mass resolution, which were projected along many lines-of-sight each. We study the shape of both their density and surface-density profiles and fit them with a variety of radial functions, including the Navarro-Frenk-White, the generalised Navarro-Frenk-White, and the Einasto density profiles. We derive concentrations and masses from these fits and investigate their distributions as a function of redshift and halo relaxation. We use the X-ray image simulator X-MAS to produce simulated Chandra observations of the halos and we use them to identify objects resembling the X-ray morphologies and masses of the clusters in the CLASH X-ray selected sample. We also derive a concentration-mass relation for strong-lensing clusters. We find that the sample of simulated halos which resemble the X-ray morphology of the CLASH clusters is composed mainly by relaxed halos, but it also contains a significant fraction of un-relaxed systems. For such a sample we measure an average 2D concentration which is ~11% higher than found for the full sample of simulated halos. After accounting for projection and selection effects, the average NFW concentrations of CLASH clusters are expected to be intermediate between those predicted in 3D for relaxed and super-relaxed halos. Matching the simulations to the individual CLASH clusters on the basis of the X-ray morphology, we expect that the NFW concentrations recovered from the lensing analysis of the CLASH clusters are in the range [3-6], with an average value of 3.87 and a standard deviation of 0.61. Simulated halos with X-ray morphologies similar to those of the CLASH clusters are affected by a modest orientation bias.

We present a new determination of the concentration-mass relation for galaxy clusters based on our comprehensive lensing analysis of 19 X-ray selected galaxy clusters from the Cluster Lensing and Supernova Survey with Hubble (CLASH). Our sample spans a redshift range between 0.19 and 0.89. We combine weak lensing constraints from the Hubble Space Telescope (HST) and from ground-based wide field data with strong lensing constraints from HST. The result are reconstructions of the surface-mass density for all CLASH clusters on multi-scale grids. Our derivation of NFW parameters yields virial masses between 0.53 x 10^15 and 1.76 x 10^15 M_sol/h and the halo concentrations are distributed around c_200c ~ 3.7 with a 1-sigma significant negative trend with cluster mass. We find an excellent 4% agreement between our measured concentrations and the expectation from numerical simulations after accounting for the CLASH selection function based on X-ray morphology. The simulations are analyzed in 2D to account for possible biases in the lensing reconstructions due to projection effects. The theoretical concentration-mass (c-M) relation from our X-ray selected set of simulated clusters and the c-M relation derived directly from the CLASH data agree at the 90% confidence level.

We aim at constraining the assembly history of clusters by studying the intra cluster light (ICL) properties, estimating its contribution to the fraction of baryons in stars, f*, and understanding possible systematics/bias using different ICL detection techniques. We developed an automated method, GALtoICL, based on the software GALAPAGOS to obtain a refined version of typical BCG+ICL maps. We applied this method to our test case MACS J1206.2-0847, a massive cluster located at z=0.44, that is part of the CLASH sample. Using deep multi-band SUBARU images, we extracted the surface brightness (SB) profile of the BCG+ICL and we studied the ICL morphology, color, and contribution to f* out to R500. We repeated the same analysis using a different definition of the ICL, SBlimit method, i.e. a SB cut-off level, to compare the results. The most peculiar feature of the ICL in MACS1206 is its asymmetric radial distribution, with an excess in the SE direction and extending towards the 2nd brightest cluster galaxy which is a Post Starburst galaxy. This suggests an interaction between the BCG and this galaxy that dates back to t <= 1.5 Gyr. The BCG+ICL stellar content is 8% of M_(*,500) and the (de-) projected baryon fraction in stars is f*=0.0177 (0.0116), in excellent agreement with recent results. The SBlimit method provides systematically higher ICL fractions and this effect is larger at lower SB limits. This is due to the light from the outer envelopes of member galaxies that contaminate the ICL. Though more time consuming, the GALtoICL method provides safer ICL detections that are almost free of this contamination. This is one of the few ICL study at redshift z > 0.3. At completion, the CLASH/VLT program will allow us to extend this analysis to a statistically significant cluster sample spanning a wide redshift range: 0.2<z<0.6.

We report observations of three gravitationally lensed supernovae (SNe) in the Cluster Lensing And Supernova survey with Hubble (CLASH) Multi-Cycle Treasury program. These objects, SN CLO12Car (z = 1.28), SN CLN12Did (z = 0.85), and SN CLA11Tib (z = 1.14), are located behind three different clusters, MACSJ1720.2+3536 (z = 0.391), RXJ1532.9+3021 (z = 0.345), and Abell 383 (z = 0.187), respectively. Each SN was detected in Hubble Space Telescope (HST) optical and infrared images. Based on photometric classification, we find that SNe CLO12Car and CLN12Did are likely to be Type Ia supernovae (SNe Ia), while the classification of SN CLA11Tib is inconclusive. Using multi-color light-curve fits to determine a standardized SN Ia luminosity distance, we infer that SN CLO12Car was approximately 1.0 +/- 0.2 mag brighter than field SNe Ia at a similar redshift and ascribe this to gravitational lens magnification. Similarly, SN CLN12Did is approximately 0.2 +/- 0.2 mag brighter than field SNe Ia. We derive independent estimates of the predicted magnification from CLASH strong+weak lensing maps of the clusters: 0.83 +/- 0.16 mag for SN CLO12Car, 0.28 +/- 0.08 mag for SN CLN12Did, and 0.43 +/- 0.11 mag for SN CLA11Tib. The two SNe Ia provide a new test of the cluster lens model predictions: we find that the magnifications based on the SN Ia brightness and those predicted by the lens maps are consistent. Our results herald the promise of future observations of samples of cluster-lensed SNe Ia (from the ground or space) to help illuminate the dark-matter distribution in clusters of galaxies, through the direct determination of absolute magnifications.

We present a strong lensing system in which a double source is imaged 5 times by 2 early-type galaxies. We take advantage in this target of the multi-band photometry obtained as part of the CLASH program, complemented by the spectroscopic data of the VLT/VIMOS and FORS2 follow-up campaign. We use a photometric redshift of 3.7 for the source and confirm spectroscopically the membership of the 2 lenses to the galaxy cluster MACS J1206.2-0847 at redshift 0.44. We exploit the excellent angular resolution of the HST/ACS images to model the 2 lenses in terms of singular isothermal sphere profiles and derive robust effective velocity dispersions of (97 +/- 3) and (240 +/- 6) km/s. The total mass distribution of the cluster is also well characterized by using only the local information contained in this lensing system, that is located at a projected distance of more than 300 kpc from the cluster luminosity center. According to our best-fitting lensing and composite stellar population models, the source is magnified by a total factor of 50 and has a luminous mass of about (1.0 +/- 0.5) x 10^{9} M_{Sun}. By combining the total and luminous mass estimates of the 2 lenses, we measure luminous over total mass fractions projected within the effective radii of 0.51 +/- 0.21 and 0.80 +/- 0.32. With these lenses we can extend the analysis of the mass properties of lens early-type galaxies by factors that are about 2 and 3 times smaller than previously done with regard to, respectively, velocity dispersion and luminous mass. The comparison of the total and luminous quantities of our lenses with those of astrophysical objects with different physical scales reveals the potential of studies of this kind for investigating the internal structure of galaxies. These studies, made possible thanks to the CLASH survey, will allow us to go beyond the current limits posed by the available lens samples in the field.

The Cluster Lensing And Supernovae survey with Hubble (CLASH) is an Hubble Space Telescope (HST) Multi-Cycle Treasury program observing 25 massive galaxy clusters. CLASH observations are carried out in 16 bands from UV to NIR to derive accurate and reliable estimates of photometric redshifts. We present the CLASH photometric redshifts and study the photometric redshift accuracy of the arcs in more detail for the case of MACS1206.2-0847. We use the publicly available Le Phare and BPZ photometric redshift codes on 17 CLASH galaxy clusters. Using Le Phare code for objects with StoN>=10, we reach a precision of 3%(1+z) for the strong lensing arcs, which is reduced to 2.4%(1+z) after removing outliers. For galaxies in the cluster field the corresponding values are 4%(1+z) and 3%(1+z). Using mock galaxy catalogues, we show that 3%(1+z) precision is what one would expect from the CLASH photometry when taking into account extinction from dust, emission lines and the finite range of SEDs included in the photo-z template library. We study photo-z results for different aperture photometry and find that the SExtractor isophotal photometry works best. Le Phare and BPZ give similar photo-z results for the strong lensing arcs as well as galaxies of the cluster field. Results are improved when optimizing the photometric aperture shape showing an optimal aperture size around 1" radius giving results which are equivalent to isophotal photometry. Tailored photometry of the arcs improve the photo-z results.

We report our analysis of MACS J0717.5+3745 using 140 and 268 GHz Bolocam data collected at the Caltech Submillimeter Observatory. We detect extended Sunyaev-Zel'dovich (SZ) effect signal at high significance in both Bolocam bands, and we employ Herschel-SPIRE observations to subtract the signal from dusty background galaxies in the 268 GHz data. We constrain the two-band SZ surface brightness toward two of the sub-clusters of MACS J0717.5+3745: the main sub-cluster (named C), and a sub-cluster identified in spectroscopic optical data to have a line-of-sight velocity of +3200 km/s (named B). We determine the surface brightness in two separate ways: via fits of parametric models and via direct integration of the images. For both sub-clusters, we find consistent surface brightnesses from both analysis methods. We constrain spectral templates consisting of relativistically corrected thermal and kinetic SZ signals, using a jointly-derived electron temperature from Chandra and XMM-Newton under the assumption that each sub-cluster is isothermal. The data show no evidence for a kinetic SZ signal toward sub-cluster C, but they do indicate a significant kinetic SZ signal toward sub-cluster B. The model-derived surface brightnesses for sub-cluster B yield a best-fit line-of-sight velocity of v_z = +3450 +- 900 km/s, with (1 - Prob[v_z > 0]) = 1.3 x 10^-5 (4.2 sigma away from 0 for a Gaussian distribution). The directly integrated sub-cluster B SZ surface brightnesses provide a best-fit v_z = +2550 +- 1050 km/s, with (1 - Prob[v_z > 0]) = 2.2 x 10^-3 (2.9 sigma).

We present a quintuply lensed z ~ 6 candidate discovered in the field of the galaxy cluster RXC J2248.7-4431 (z ~ 0.348) targeted within the Cluster Lensing and Supernova survey with Hubble (CLASH) and selected in the deep HST Frontier Fields survey. Thanks to the CLASH 16-band HST imaging, we identify the quintuply lensed z ~ 6 candidate as an optical dropout in the inner region of the cluster, the brightest image having magAB=24.81+-0.02 in the f105w filter. We perform a detailed photometric analysis to verify its high-z and lensed nature. We get as photometric redshift z_phot ~ 5.9, and given the extended nature and NIR colours of the lensed images, we rule out low-z early type and galactic star contaminants. We perform a strong lensing analysis of the cluster, using 13 families of multiple lensed images identified in the HST images. Our final best model predicts the high-z quintuply lensed system with a position accuracy of 0.8''. The magnifications of the five images are between 2.2 and 8.3, which leads to a delensed UV luminosity of L_1600 ~ 0.5L*_1600 at z=6. We also estimate the UV slope from the observed NIR colours, finding a steep beta=-2.89+-0.38. We use singular and composite stellar population SEDs to fit the photometry of the hiz candidate, and we conclude that it is a young (age <300 Myr) galaxy with mass of M ~ 10^8Msol, subsolar metallicity (Z<0.2Zsol) and low dust content (AV ~ 0.2-0.4).

We present VIsible Multi-Object Spectrograph (VIMOS) observations of a z 6 galaxy quintuply imaged by the Frontier Fields galaxy cluster RXC J2248.7-4431 (z=0.348). This sub-L^*, high-z galaxy has been recently discovered by Monna et al. (2013) using dropout techniques with the 16-band HST photometry acquired as part of the Cluster Lensing And Supernova survey with Hubble (CLASH). Obtained as part of the CLASH-VLT survey, the VIMOS medium-resolution spectra of this source show a very faint continuum between ~8700A and ~9300A and a prominent emission line at 8643A, which can be readily identified with Lyman-alpha at z=6.110. The emission line exhibits an asymmetric profile, with a more pronounced red wing. The rest-frame equivalent width of the line is EW=79+-10A. After correcting for magnification, the star-formation rate (SFR) estimated from the Lya line is SFR(Lya)=11 M_{sol}/yr and that estimated from the UV data is SFR(UV)=3 M_{sol}/yr. We estimate that the effective radius of the source is R_e<~0.4 kpc, which implies a star formation surface mass density Sigma_{SFR}>6 M_{sol}/yr/kpc^2 and, using the Kennicutt-Schmidt relation, a gas surface mass density Sigma_{gas}>10^3 M_{sol}/pc^2. Our results support the idea that this magnified, distant galaxy is a young and compact object with 0.4 L^* at z=6, with comparable amount of mass in gas and stars. Future follow-up observations with ALMA will provide valuable insight into the SFR and molecular gas content of this source. In the spirit of the Frontier Fields initiative, we also publish the redshifts of several multiply imaged sources and other background objects which will help improving the strong lensing model of this galaxy cluster.

We present PyGFit, a program designed to measure PSF-matched photometry from images with disparate pixel scales and PSF sizes. While PyGFit has a number of uses, its primary purpose is to extract robust spectral energy distributions (SEDs) from crowded images. It does this by fitting blended sources in crowded, low resolution images with models generated from a higher resolution image. This approach minimizes the impact of crowding and also yields consistently measured fluxes in different filters, minimizing systematic uncertainty in the final SEDs. We present an example of applying PyGFit to real data and perform simulations to test its fidelity. The uncertainty in the best-fit flux rises sharply as a function of nearest-neighbor distance for objects with a neighbor within 60% of the PSF size. Similarly, the uncertainty increases quickly for objects blended with a neighbor more than four times brighter. For all other objects the fidelity of PyGFit's results depends only on flux, and the uncertainty is primarily limited by sky noise.

We analyze the star formation properties of 16 infrared-selected, spectroscopically confirmed galaxy clusters at $1 < z < 1.5$ from the Spitzer/IRAC Shallow Cluster Survey (ISCS). We present new spectroscopic confirmation for six of these high-redshift clusters, five of which are at $z>1.35$. Using infrared luminosities measured with deep Spitzer/MIPS observations at 24 $\mu$m, along with robust optical+IRAC photometric redshifts and SED-fitted stellar masses, we present the dust-obscured star-forming fractions, star formation rates and specific star formation rates in these clusters as functions of redshift and projected clustercentric radius. We find that $z\sim 1.4$ represents a transition redshift for the ISCS sample, with clear evidence of an unquenched era of cluster star formation at earlier times. Beyond this redshift the fraction of star-forming cluster members increases monotonically toward the cluster centers. Indeed, the specific star formation rate in the cores of these distant clusters is consistent with field values at similar redshifts, indicating that at $z>1.4$ environment-dependent quenching had not yet been established in ISCS clusters. Combining these observations with complementary studies showing a rapid increase in the AGN fraction, a stochastic star formation history, and a major merging episode at the same epoch in this cluster sample, we suggest that the starburst activity is likely merger-driven and that the subsequent quenching is due to feedback from merger-fueled AGN. The totality of the evidence suggests we are witnessing the final quenching period that brings an end to the era of star formation in galaxy clusters and initiates the era of passive evolution.

The galaxy cluster MACS J0717.5+3745 (z=0.55) is the largest known cosmic lens, with complex internal structures seen in deep X-ray, Sunyaev-Zel'dovich effect and dynamical observations. We perform a combined weak and strong lensing analysis with wide-field BVRi'z' Subaru/Suprime-Cam observations and 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH). We find consistent weak distortion and magnification measurements of background galaxies, and combine these signals to construct an optimally estimated radial mass profile of the cluster and its surrounding large-scale structure out to 5 Mpc/h. We find consistency between strong-lensing and weak-lensing in the region where these independent data overlap, <500 kpc/h. The two-dimensional weak-lensing map reveals a clear filamentary structure traced by distinct mass halos. We model the lensing shear field with 9 halos, including the main cluster, corresponding to mass peaks detected above 2.5\sigma_\kappa. The total mass of the cluster as determined by the different methods is M_{vir}=(2.8\pm0.4) \times 10^15 M_sun. Although this is the most massive cluster known at z>0.5, in terms of extreme value statistics we conclude that the mass of MACS J0717.5+3745 by itself is not in serious tension with LambdaCDM, representing only a ~2{\sigma} departure above the maximum simulated halo mass at this redshift.

We provide a new observational test for a key prediction of the \Lambda CDM cosmological model: the contributions of mergers with different halo-to-main-cluster mass ratios to cluster-sized halo growth. We perform this test by dynamically analyzing seven galaxy clusters, spanning the redshift range $0.13 < z_c < 0.45$ and caustic mass range $0.4-1.5$ $10^{15} h_{0.73}^{-1}$ M$_{\odot}$, with an average of 293 spectroscopically-confirmed bound galaxies to each cluster. The large radial coverage (a few virial radii), which covers the whole infall region, with a high number of spectroscopically identified galaxies enables this new study. For each cluster, we identify bound galaxies. Out of these galaxies, we identify infalling and accreted halos and estimate their masses and their dynamical states. Using the estimated masses, we derive the contribution of different mass ratios to cluster-sized halo growth. For mass ratios between ~0.2 and ~0.7, we find a ~1 $\sigma$ agreement with \Lambda CDM expectations based on the Millennium simulations I and II. At low mass ratios, $\lesssim 0.2$, our derived contribution is underestimated since the detection efficiency decreases at low masses, $\sim 2 \times 10^{14}$ $h_{0.73}^{-1}$ M$_{\odot}$. At large mass ratios, $\gtrsim 0.7$, we do not detect halos probably because our sample, which was chosen to be quite X-ray relaxed, is biased against large mass ratios. Therefore, at large mass ratios, the derived contribution is also underestimated.

We use an unprecedented data-set of about 600 redshifts for cluster members, obtained as part of a VLT/VIMOS large programme, to constrain the mass profile of the z=0.44 cluster MACS J1206.2-0847 over the radial range 0-5 Mpc (0-2.5 virial radii) using the MAMPOSSt and Caustic methods. We then add external constraints from our previous gravitational lensing analysis. We invert the Jeans equation to obtain the velocity-anisotropy profiles of cluster members. With the mass-density and velocity-anisotropy profiles we then obtain the first determination of a cluster pseudo-phase-space density profile. The kinematics and lensing determinations of the cluster mass profile are in excellent agreement. This is very well fitted by a NFW model with mass M200=(1.4 +- 0.2) 10^15 Msun and concentration c200=6 +- 1, only slightly higher than theoretical expectations. Other mass profile models also provide acceptable fits to our data, of (slightly) lower (Burkert, Hernquist, and Softened Isothermal Sphere) or comparable (Einasto) quality than NFW. The velocity anisotropy profiles of the passive and star-forming cluster members are similar, close to isotropic near the center and increasingly radial outside. Passive cluster members follow extremely well the theoretical expectations for the pseudo-phase-space density profile and the relation between the slope of the mass-density profile and the velocity anisotropy. Star-forming cluster members show marginal deviations from theoretical expectations. This is the most accurate determination of a cluster mass profile out to a radius of 5 Mpc, and the only determination of the velocity-anisotropy and pseudo-phase-space density profiles of both passive and star-forming galaxies for an individual cluster [abridged]

We use high-resolution cosmological hydrodynamic simulations to study the angular momentum acquisition of gaseous halos around Milky Way sized galaxies. We find that cold mode accreted gas enters a galaxy halo with ~70% more specific angular momentum than dark matter averaged over cosmic time (though with a very large dispersion). In fact, we find that all matter has a higher spin parameter when measured at accretion than when averaged over the entire halo lifetime, and is well characterized by \lambda~0.1, at accretion. Combined with the fact that cold flow gas spends a relatively short time (1-2 dynamical times) in the halo before sinking to the center, this naturally explains why cold flow halo gas has a specific angular momentum much higher than that of the halo and often forms "cold flow disks". We demonstrate that the higher angular momentum of cold flow gas is related to the fact that it tends to be accreted along filaments.

Growing observational evidence now indicates that nebular line emission has a significant impact on the rest-frame optical fluxes of z~5-7 galaxies observed with Spitzer. This line emission makes z~5-7 galaxies appear more massive, with lower specific star formation rates. However, corrections for this line emission have been very difficult to perform reliably due to huge uncertainties on the overall strength of such emission at z>~5.5. Here, we present the most direct observational evidence yet for ubiquitous high-EW [OIII]+Hbeta line emission in Lyman-break galaxies at z~7, while also presenting a strategy for an improved measurement of the sSFR at z~7. We accomplish this through the selection of bright galaxies in the narrow redshift window z~6.6-7.0 where the IRAC 4.5 micron flux provides a clean measurement of the stellar continuum light. Observed 4.5 micron fluxes in this window contrast with the 3.6 micron fluxes which are contaminated by the prominent [OIII]+Hbeta lines. To ensure a high S/N for our IRAC flux measurements, we consider only the brightest (H_{160}<26 mag) magnified galaxies we have identified in CLASH and other programs targeting galaxy clusters. Remarkably, the mean rest-frame optical color for our bright seven-source sample is very blue, [3.6]-[4.5]=-0.9+/-0.3. Such blue colors cannot be explained by the stellar continuum light and require that the rest-frame EW of [OIII]+Hbeta be greater than 637 Angstroms for the average source. The bluest four sources from our seven-source sample require an even more extreme EW of 1582 Angstroms. Our derived lower limit for the mean [OIII]+Hbeta EW could underestimate the true EW by ~2x based on a simple modeling of the redshift distribution of our sources. We can also set a robust lower limit of >~4 Gyr^-1 on the specific star formation rates based on the mean SED for our seven-source sample. (abridged)

In this work we analyze the mass distribution of MACSJ1206.2-0847, especially focusing on the halo properties of its cluster members. The cluster appears relaxed in its X-ray emission, but has significant amounts of intracluster light which is not centrally concentrated, suggesting that galaxy-scale interactions are still ongoing despite the overall relaxed state. The cluster lenses 12 background galaxies into multiple images and one galaxy at $z=1.033$ into a giant arc and its counterimage. The multiple image positions and the surface brightness distribution (SFB) of the arc which is bent around several cluster members are sensitive to the cluster galaxy halo properties. We model the cluster mass distribution with a NFW profile and the galaxy halos with two parameters for the mass normalization and extent of a reference halo assuming scalings with their observed NIR--light. We match the multiple image positions at an r.m.s. level of $0.85\arcsec$ and can reconstruct the SFB distribution of the arc in several filters to a remarkable accuracy based on this cluster model. The length scale where the enclosed galaxy halo mass is best constrained is about 5 effective radii -- a scale in between those accessible to dynamical and field strong lensing mass estimates on one hand and galaxy--galaxy weak lensing results on the other hand. The velocity dispersion and halo size of a galaxy with $m_{\rm 160W,AB}=19.2$ or $M_{\rm B,Vega}=-20.7$ are $\sigma=150 \rm kms^{-1}$ and $r\approx 26\pm 6 \rm kpc$, indicating that the halos of the cluster galaxies are tidally stripped. We also reconstruct the unlensed source (which is smaller by a factor of $\sim5.8$ in area), demonstrating the increase of morphological information due to lensing and conclude that this galaxy has likely star--forming spiral arms with a red (older) central component.

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.

We describe Sunyaev-Zel'dovich (SZ) effect measurements and analysis of the intracluster medium (ICM) pressure profiles of a set of 45 massive galaxy clusters imaged using Bolocam at the Caltech Submillimeter Observatory. We have used masses determined from Chandra X-ray observations to scale each cluster's profile by the overdensity radius R500 and the mass-and-redshift-dependent normalization factor P500. We deproject the average pressure profile of our sample into 13 logarithmically spaced radial bins between 0.07R500 and 3.5R500. We find that a generalized Navarro, Frenk, and White (gNFW) profile describes our data with sufficient goodness-of-fit and best-fit parameters (C500, alpha, beta, gamma, P0 = 1.18, 0.86, 3.67, 0.67, 4.29). We also use the X-ray data to define cool-core and disturbed subsamples of clusters, and we constrain the average pressure profiles of each of these subsamples. We find that given the precision of our data the average pressure profiles of disturbed and cool-core clusters are consistent with one another at R>~0.15R500, with cool-core systems showing indications of higher pressure at R<~0.15R500. In addition, for the first time, we place simultaneous constraints on the mass scaling of cluster pressure profiles, their ensemble mean profile, and their radius-dependent intrinsic scatter between 0.1R500 and 2.0R500. The scatter among profiles is minimized at radii between ~0.2R500 and ~0.5R500, with a value of ~20%. The best-fit mass scaling has a power-law slope of 0.49, which is shallower than the nominal prediction of 2/3 from self-similar hydrostatic equilibrium models. These results for the intrinsic scatter and mass scaling are largely consistent with previous analyses, most of which have relied heavily on X-ray derived pressures of clusters at significantly lower masses and redshifts compared to our sample.

We present observations of a luminous galaxy at redshift z=6.573 --- the end of the reioinization epoch --- which has been spectroscopically confirmed twice. The first spectroscopic confirmation comes from slitless HST ACS grism spectra from the PEARS survey (Probing Evolution And Reionization Spectroscopically), which show a dramatic continuum break in the spectrum at restframe 1216 A wavelength. The second confirmation is done with Keck + DEIMOS. The continuum is not clearly detected with ground-based spectra, but high wavelength resolution enables the Lyman alpha emission line profile to be determined. We compare the line profile to composite line profiles at redshift z=4.5. The Lyman alpha line profile shows no signature of a damping wing attenuation, confirming that the intergalactic gas is ionized at redshift z=6.57. Spectra of Lyman breaks at yet higher redshifts will be possible using comparably deep observations with IR-sensitive grisms, even at redshifts where Lyman alpha is too attenuated by the neutral IGM to be detectable using traditional spectroscopy from the ground.

We describe in detail our characterization of the compact radio source population in 140 GHz Bolocam observations of a set of 45 massive galaxy clusters. We use a combination of 1.4 and 30 GHz data to select a total of 28 probable cluster-member radio galaxies and also to predict their 140 GHz flux densities. All of these galaxies are steep-spectrum radio sources and they are found preferentially in the cool-core clusters within our sample. In particular, 11 of the 12 brightest cluster member radio sources are associated with cool-core systems. Although none of the individual galaxies are robustly detected in the Bolocam data, the ensemble-average flux density at 140 GHz is consistent with, but slightly lower than, the extrapolation from lower frequencies assuming a constant spectral index. In addition, our data indicate an intrinsic scatter of 30 percent around the power-law extrapolated flux densities at 140 GHz, although our data do not tightly constrain this scatter. For our cluster sample, which is composed of high-mass and moderate-redshift systems, we find that the maximum fractional change in the Sunyaev-Zel'dovich signal integrated over any single cluster due to the presence of these radio sources is 20 percent, and only 1/4 of the clusters show a fractional change of more than 1 percent. The amount of contamination is strongly dependent on cluster morphology, and nearly all of the clusters with more than 1 percent contamination are cool-core systems. This result indicates that radio contamination is not significant compared to current noise levels in 140 GHz images of massive clusters and is in good agreement with the level of radio contamination found in previous results based on lower frequency data or simulations.

We perform a strong-lensing analysis of the merging galaxy cluster MACS J0416.1-2403 (M0416; z=0.42) in recent CLASH/HST observations. We identify 70 new multiple images and candidates of 23 background sources in the range 0.7<z_{phot}<6.14 including two probable high-redshift dropouts, revealing a highly elongated lens with axis ratio ~5:1, and a major axis of ~100\arcsec (z_{s}~2). Compared to other well-studied clusters, M0416 shows an enhanced lensing efficiency. Although the critical area is not particularly large (~0.6 \square\arcmin; z_{s}~2), the number of multiple images, per critical area, is anomalously high. We calculate that the observed elongation boosts the number of multiple images, \emph{per critical area}, by a factor of ~2.5\times, due to the increased ratio of the caustic area relative to the critical area. Additionally, we find that the observed separation between the two main mass components enlarges the critical area by a factor of ~2. These geometrical effects can account for the high number (density) of multiple images observed. We find in numerical simulations, that only ~4% of the clusters (with M_{vir}>6 x 10^{14} h^{-1}M_{\odot}) exhibit as elongated critical curves as M0416.

We present a candidate for the most distant galaxy known to date with a photometric redshift z = 10.7 +0.6 / -0.4 (95% confidence limits; with z < 9.5 galaxies of known types ruled out at 7.2-sigma). This J-dropout Lyman Break Galaxy, named MACS0647-JD, was discovered as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). We observe three magnified images of this galaxy due to strong gravitational lensing by the galaxy cluster MACSJ0647.7+7015 at z = 0.591. The images are magnified by factors of ~8, 7, and 2, with the brighter two observed at ~26th magnitude AB (~0.15 uJy) in the WFC3/IR F160W filter (~1.4 - 1.7 um) where they are detected at >~ 12-sigma. All three images are also confidently detected at >~ 6-sigma in F140W (~1.2 - 1.6 um), dropping out of detection from 15 lower wavelength HST filters (~0.2 - 1.4 um), and lacking bright detections in Spitzer/IRAC 3.6um and 4.5um imaging (~3.2 - 5.0 um). We rule out a broad range of possible lower redshift interlopers, including some previously published as high redshift candidates. Our high redshift conclusion is more conservative than if we had neglected a Bayesian photometric redshift prior. Given CLASH observations of 17 high mass clusters to date, our discoveries of MACS0647-JD at z ~ 10.8 and MACS1149-JD1 at z ~ 9.6 are consistent with a lensed luminosity function extrapolated from lower redshifts. This would suggest that low luminosity galaxies could have reionized the universe. However given the significant uncertainties based on only two galaxies, we cannot yet rule out the sharp drop off in number counts at z >~ 10 suggested by field searches.

We describe a methodology to accurately compute halo mass functions, progenitor mass functions, merger rates and merger trees in non-cold dark matter universes using a self-consistent treatment of the generalized extended Press-Schechter formalism. Our approach permits rapid exploration of the subhalo population of galactic halos in dark matter models with a variety of different particle properties or universes with rolling, truncated, or more complicated power spectra. We make detailed comparisons of analytically derived mass functions and merger histories with recent warm dark matter cosmological N-body simulations, and find excellent agreement. We show that, once the accretion of smoothly distributed matter is accounted for, coarse-grained statistics such as the mass accretion history of halos can be almost indistinguishable between cold and warm dark matter cases. However, the halo mass function and progenitor mass functions differ significantly, with the warm dark matter cases being strongly suppressed below the free-streaming scale of the dark matter. We demonstrate the importance of using the correct solution for the excursion set barrier first-crossing distribution in warm dark matter - if the solution for a flat barrier is used instead the truncation of the halo mass function is much slower, leading to an overestimate of the number of low mass halos.

We use cosmological simulations to study the effects of self-interacting dark matter (SIDM) on the density profiles and substructure counts of dark matter halos from the scales of spiral galaxies to galaxy clusters, focusing explicitly on models with cross sections over dark matter particle mass \sigma/m = 1 and 0.1 cm^2/g. Our simulations rely on a new SIDM N-body algorithm that is derived self-consistently from the Boltzmann equation and that reproduces analytic expectations in controlled numerical experiments. We find that well-resolved SIDM halos have constant-density cores, with significantly lower central densities than their CDM counterparts. In contrast, the subhalo content of SIDM halos is only modestly reduced compared to CDM, with the suppression greatest for large hosts and small halo-centric distances. Moreover, the large-scale clustering and halo circular velocity functions in SIDM are effectively identical to CDM, meaning that all of the large-scale successes of CDM are equally well matched by SIDM. From our largest cross section runs we are able to extract scaling relations for core sizes and central densities over a range of halo sizes and find a strong correlation between the core radius of an SIDM halo and the NFW scale radius of its CDM counterpart. We construct a simple analytic model, based on CDM scaling relations, that captures all aspects of the scaling relations for SIDM halos. Our results show that halo core densities in \sigma/m = 1 cm^2/g models are too low to match observations of galaxy clusters, low surface brightness spirals (LSBs), and dwarf spheroidal galaxies. However, SIDM with \sigma/m ~ 0.1 cm^2/g appears capable of reproducing reported core sizes and central densities of dwarfs, LSBs, and galaxy clusters without the need for velocity dependence. (abridged)

Hubble Space Telescope images of the galaxy cluster Abell 2261, obtained as part of the Cluster Lensing And Supernova survey with Hubble, show that the brightest galaxy in the cluster, A2261-BCG, has the largest core yet detected in any galaxy. The cusp radius of A2261-BCG is 3.2 kpc, twice as big as the next largest core known, and ~3x bigger than those typically seen in the most luminous BCGs. The morphology of the core in A2261-BCG is also unusual, having a flat or even slightly-depressed interior surface brightness profile, rather than the typical shallow cusp. This implies that the galaxy has a core with constant or even centrally decreasing stellar density. Interpretation of the core as an end product of the "scouring" action of a binary supermassive black hole implies a total black hole mass ~1E+10 M_sun from the extrapolation of most relationships between core structure and black hole mass. The core falls 1-sigma above the cusp-radius versus galaxy luminosity relation. Its large size in real terms, and the extremely large black hole mass required to generate it, raise the possibility that the core has been enlarged by additional processes, such as the ejection of the black holes that originally generated the core. The flat central stellar density profile is consistent with this hypothesis. The core is also displaced by 0.7 kpc from the center of the surrounding envelope, consistent with a local dynamical perturbation of the core.

We examine the inner mass distribution of the relaxed galaxy cluster Abell 383 in deep 16-band HST/ACS+WFC3 imaging taken as part of the CLASH multi-cycle treasury program. Our program is designed to study the dark matter distribution in 25 massive clusters, and balances depth with a wide wavelength coverage to better identify lensed systems and generate precise photometric redshifts. This information together with the predictive strength of our strong-lensing analysis method identifies 13 new multiply-lensed images and candidates, so that a total of 27 multiple-images of 9 systems are used to tightly constrain the inner mass profile, $d\log \Sigma/d\log r\simeq -0.6\pm 0.1$ (r<160 kpc). We find consistency with the standard distance-redshift relation for the full range spanned by the lensed images, 1.01<z<6.03, with the higher redshift sources deflected through larger angles as expected. The inner mass profile derived here is consistent with the results of our independent weak-lensing analysis of wide-field Subaru images, with good agreement in the region of overlap. The overall mass profile is well fitted by an NFW profile with M_{vir}=(5.37^{+0.70}_{-0.63}\pm 0.26) x 10^{14}M_{\odot}/h and a relatively high concentration, c_{vir}=8.77^{+0.44}_{-0.42}\pm 0.23, which lies above the standard c-M relation similar to other well-studied clusters. The critical radius of Abell 383 is modest by the standards of other lensing clusters, r_{E}\simeq16\pm2\arcsec (for z_s=2.55), so the relatively large number of lensed images uncovered here with precise photometric redshifts validates our imaging strategy for the CLASH survey. In total we aim to provide similarly high-quality lensing data for 25 clusters, 20 of which are X-ray selected relaxed clusters, enabling a precise determination of the representative mass profile free from lensing bias. (ABRIDGED)

We derive an accurate mass distribution of the galaxy cluster MACS J1206.2-0847 (z=0.439) from a combined weak-lensing distortion, magnification, and strong-lensing analysis of wide-field Subaru BVRIz' imaging and our recent 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) program. We find good agreement in the regions of overlap between several weak and strong lensing mass reconstructions using a wide variety of modeling methods, ensuring consistency. The Subaru data reveal the presence of a surrounding large scale structure with the major axis running approximately north-west south-east (NW-SE), aligned with the cluster and its brightest galaxy shapes, showing elongation with a \sim 2:1 axis ratio in the plane of the sky. Our full-lensing mass profile exhibits a shallow profile slope dln\Sigma/dlnR\sim -1 at cluster outskirts (R>1Mpc/h), whereas the mass distribution excluding the NW-SE excess regions steepens further out, well described by the Navarro-Frenk-White form. Assuming a spherical halo, we obtain a virial mass M_{vir}=(1.1\pm 0.2\pm 0.1)\times 10^{15} M_{sun}/h and a halo concentration c_{vir} = 6.9\pm 1.0\pm 1.2 (\sim 5.7 when the central 50kpc/h is excluded), which falls in the range 4< <c> <7 of average c(M,z) predictions for relaxed clusters from recent Lambda cold dark matter simulations. Our full lensing results are found to be in agreement with X-ray mass measurements where the data overlap, and when combined with Chandra gas mass measurements, yield a cumulative gas mass fraction of 13.7^{+4.5}_{-3.0}% at 0.7Mpc/h (\approx 1.7r_{2500}), a typical value observed for high mass clusters.

We present the Advanced Camera for Surveys General Catalog (ACS-GC), a photometric and morphological database using publicly available data obtained with the Advanced Camera for Surveys (ACS) instrument on the Hubble Space Telescope. The goal of the ACS-GC database is to provide a large statistical sample of galaxies with reliable structural and distance measurements to probe the evolution of galaxies over a wide range of look-back times. The ACS-GC includes approximately 470,000 astronomical sources (stars + galaxies) derived from the AEGIS, COSMOS, GEMS, and GOODS surveys. Galapagos was used to construct photometric (SExtractor) and morphological (Galfit) catalogs. The analysis assumes a single S\'ersic model for each object to derive quantitative structural parameters. We include publicly available redshifts from the DEEP2, COMBO-17, TKRS, PEARS, ACES, CFHTLS,and zCOSMOS surveys to supply redshifts (spectroscopic and photometric) for a considerable fraction (~74%) of the imaging sample. The ACS-GC includes color postage stamps, Galfit residual images, and photometry, structural parameters, and redshifts combined into a single catalog.

The early Universe at redshift z\sim6-11 marks the reionization of the intergalactic medium, following the formation of the first generation of stars. However, those young galaxies at a cosmic age of \lesssim 500 million years (Myr, at z \gtrsim 10) remain largely unexplored as they are at or beyond the sensitivity limits of current large telescopes. Gravitational lensing by galaxy clusters enables the detection of high-redshift galaxies that are fainter than what otherwise could be found in the deepest images of the sky. We report the discovery of an object found in the multi-band observations of the cluster MACS1149+22 that has a high probability of being a gravitationally magnified object from the early universe. The object is firmly detected (12 sigma) in the two reddest bands of HST/WFC3, and not detected below 1.2 {\mu}m, matching the characteristics of z\sim9 objects. We derive a robust photometric redshift of z = 9.6 \pm 0.2, corresponding to a cosmic age of 490 \pm 15Myr (i.e., 3.6% of the age of the Universe). The large number of bands used to derive the redshift estimate make it one of the most accurate estimates ever obtained for such a distant object. The significant magnification by cluster lensing (a factor of \sim15) allows us to analyze the object's ultra-violet and optical luminosity in its rest-frame, thus enabling us to constrain on its stellar mass, star-formation rate and age. If the galaxy is indeed at such a large redshift, then its age is less than 200 Myr (at the 95% confidence level), implying a formation redshift of zf \lesssim 14. The object is the first z>9 candidate that is bright enough for detailed spectroscopic studies with JWST, demonstrating the unique potential of galaxy cluster fields for finding highly magnified, intrinsically faint galaxies at the highest redshifts.

We study the shapes of Milky Way satellites in the context of the tidal stirring scenario for the formation of dwarf spheroidal galaxies. The standard procedures used to measure shapes involve smoothing and binning of data and thus may not be sufficient to detect structural properties like bars, which are usually subtle in low surface brightness systems. Taking advantage of the fact that in nearby dwarfs photometry of individual stars is available we introduce discrete measures of shape based on the two-dimensional inertia tensor and the Fourier bar mode. We apply these measures of shape first to a variety of simulated dwarf galaxies formed via tidal stirring of disks embedded in dark matter halos and orbiting the Milky Way. In addition to strong mass loss and randomization of stellar orbits, the disks undergo morphological transformation that typically involves the formation of a triaxial bar after the first pericenter passage. These tidally induced bars persist for a few Gyr before being shortened towards a more spherical shape if the tidal force is strong enough. We test this prediction by measuring in a similar way the shape of nearby dwarf galaxies, satellites of the Milky Way. We detect inner bars in Ursa Minor, Sagittarius, LMC and possibly Carina. In addition, six out of eleven studied dwarfs show elongated stellar distributions in the outer parts that may signify transition to tidal tails. We thus find the shapes of Milky Way satellites to be consistent with the predictions of the tidal stirring model.

We report the discovery of a z_{phot}=6.18^{+0.05}_{-0.07} (95% confidence level) dwarf galaxy, lensed into four images by the galaxy cluster MACS J0329.6-0211 (z_{l}=0.45). The galaxy is observed as a high-redshift dropout in HST/ACS/WFC3 CLASH and Spitzer/IRAC imaging. Its redshift is securely determined due to a clear detection of the Lyman-break in the 18-band photometry, making this galaxy one of the highest-redshift multiply-lensed objects known to date with an observed magnitude of F125W=24.00\pm0.04 AB mag for its highest-magnified image. We also present the first strong-lensing analysis of this cluster uncovering 15 additional multiply-imaged candidates of five lower-redshift sources spanning the range z_{s}~2-4. The mass model independently supports the high photometric redshift and reveals magnifications of 11.6^{+8.9}_{-4.1}, 17.6^{+6.2}_{-3.9}, 3.9^{+3.0}_{-1.7}, and 3.7^{+1.3}_{-0.2}, respectively, for the four images of the high-redshift galaxy. With this we construct a source image with a physical resolution of ~200 pc when the universe was ~0.9 Gyr old, where the z~6.2 galaxy occupies a source-plane area of approximately 2.2 kpc^{2}. Modeling the observed spectral energy distribution using population synthesis models, we find a demagnified stellar mass of ~10^{9} {M}_{sun}, subsolar metallicity (Z/Z_{sun}~0.5), low dust content (A_{V}~0.1 mag), a demagnified star formation rate (SFR) of ~3.2 {M}_{sun} yr^{-1}, and a specific SFR of ~3.4 Gyr^{-1}, all consistent with the properties of local dwarf galaxies.

We precisely constrain the inner mass profile of Abell 2261 (z=0.225) for the first time and determine this cluster is not "over-concentrated" as found previously, implying a formation time in agreement with {\Lambda}CDM expectations. These results are based on strong lensing analyses of new 16-band HST imaging obtained as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). Combining this with revised weak lensing analyses of Subaru wide field imaging with 5-band Subaru + KPNO photometry, we place tight new constraints on the halo virial mass M_vir = 2.2\pm0.2\times10^15 M\odot/h70 (within r \approx 3 Mpc/h70) and concentration c = 6.2 \pm 0.3 when assuming a spherical halo. This agrees broadly with average c(M,z) predictions from recent {\Lambda}CDM simulations which span 5 <~ <c> <~ 8. Our most significant systematic uncertainty is halo elongation along the line of sight. To estimate this, we also derive a mass profile based on archival Chandra X-ray observations and find it to be ~35% lower than our lensing-derived profile at r2500 ~ 600 kpc. Agreement can be achieved by a halo elongated with a ~2:1 axis ratio along our line of sight. For this elongated halo model, we find M_vir = 1.7\pm0.2\times10^15 M\odot/h70 and c_vir = 4.6\pm0.2, placing rough lower limits on these values. The need for halo elongation can be partially obviated by non-thermal pressure support and, perhaps entirely, by systematic errors in the X-ray mass measurements. We estimate the effect of background structures based on MMT/Hectospec spectroscopic redshifts and find these tend to lower Mvir further by ~7% and increase cvir by ~5%.

The Cluster Lensing And Supernova survey with Hubble (CLASH) is a 524-orbit multi-cycle treasury program to use the gravitational lensing properties of 25 galaxy clusters to accurately constrain their mass distributions. The survey, described in detail in this paper, will definitively establish the degree of concentration of dark matter in the cluster cores, a key prediction of CDM. The CLASH cluster sample is larger and less biased than current samples of space-based imaging studies of clusters to similar depth, as we have minimized lensing-based selection that favors systems with overly dense cores. Specifically, twenty CLASH clusters are solely X-ray selected. The X-ray selected clusters are massive (kT > 5 keV; 5 - 30 x 10^14 M_solar) and, in most cases, dynamically relaxed. Five additional clusters are included for their lensing strength (Einstein radii > 35 arcsec at z_source = 2) to further quantify the lensing bias on concentration, to yield high resolution dark matter maps, and to optimize the likelihood of finding highly magnified high-redshift (z > 7) galaxies. The high magnification, in some cases, provides angular resolutions unobtainable with any current UVOIR facility and can yield z > 7 candidates bright enough for spectroscopic follow-up. A total of 16 broadband filters, spanning the near-UV to near-IR, are employed for each 20-orbit campaign on each cluster. These data are used to measure precise (sigma_phz < 0.02(1+z)) photometric redshifts for dozens of newly discovered multiply-lensed images per cluster. Observations of each cluster are spread over 8 epochs to enable a search, primarily in the parallel fields, for Type Ia supernovae at z > 1 to improve constraints on the time dependence of the dark energy equation of state and the evolution of such supernovae in an epoch when the universe is matter dominated.

We present a strong-lensing analysis of the galaxy cluster MACS J1206.2-0847 ($z$=0.44) using UV, Optical, and IR, HST/ACS/WFC3 data taken as part of the CLASH multi-cycle treasury program, with VLT/VIMOS spectroscopy for some of the multiply-lensed arcs. The CLASH observations, combined with our mass-model, allow us to identify 47 new multiply-lensed images of 12 distant sources. These images, along with the previously known arc, span the redshift range $1\la z\la5.5$, and thus enable us to derive a detailed mass distribution and to accurately constrain, for the first time, the inner mass-profile of this cluster. We find an inner profile slope of $d\log \Sigma/d\log \theta\simeq -0.55\pm 0.1$ (in the range [1\arcsec, 53\arcsec], or $5\la r \la300$ kpc), as commonly found for relaxed and well-concentrated clusters. Using the many systems uncovered here we derive credible critical curves and Einstein radii for different source redshifts. For a source at $z_{s}\simeq2.5$, the critical curve encloses a large area with an effective Einstein radius of $\theta_{E}=28\pm3\arcsec$, and a projected mass of $1.34\pm0.15\times10^{14} M_{\odot}$. From the current understanding of structure formation in concordance cosmology, these values are relatively high for clusters at $z\sim0.5$, so that detailed studies of the inner mass distribution of clusters such as MACS J1206.2-0847 can provide stringent tests of the $\Lambda$CDM paradigm.

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5\sigma\ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5\sigma\ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

We present results from a search for gravitationally lensed galaxies present in the Hubble Space Telescope (HST)/Wide Field Camera-3 (WFC3) images of the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). We present one bona fide lens system in UDS and two compact lens candidates in the GOODS-S field. The lensing system in UDS involves two background galaxies, one at z=1.847 lensed to an arc and a counterimage, and the second at a photometric redshift of z=2.32^{+0.10}_{-0.06} lensed to a double image. We reconstruct the lensed sources in the source plane and find in each of the two cases the sources can be separated to a pair of galaxies. The sources responsible for the arc are compact with effective radii of 0.3 to 0.4 kpc in WFC3 J_{125}-band and a total stellar mass and a star-formation rate of 2.1_{-0.4}^{+2.4} times 10^7 M_sun and 2.3_{-1.7}^{+ 0.6} M_sun yr^{-1}, respectively.The abnormally high H_{160}-band flux of this source is likely due to OIII emission lines with a rest-frame equivalent width about 700 Angstroms for OIII 5007 Angstroms. The sources responsible for the double image have corresponding values of about 0.4 to 0.5 kpc, 1.4_{-0.8}^{+1.9} times 10^9 M_sun, and 8.7_{-7.0}^{+11.1} M_sun yr^{-1}. Once completed CANDELS is expected to contain about 15 lensing systems and will allow statistical studies on both lensing mass profiles and z ~ 2 lensed galaxies.

We present a detailed strong lensing, weak lensing and X-ray analysis of Abell 2744 (z = 0.308), one of the most actively merging galaxy clusters known. It appears to have unleashed `dark', `ghost', `bullet' and `stripped' substructures, each ~10^14 solar masses. The phenomenology is complex and will present a challenge for numerical simulations to reproduce. With new, multiband HST imaging, we identify 34 strongly-lensed images of 11 galaxies around the massive Southern `core'. Combining this with weak lensing data from HST, VLT and Subaru, we produce the most detailed mass map of this cluster to date. We also perform an independent analysis of archival Chandra X-ray imaging. Our analyses support a recent claim that the Southern core and Northwestern substructure are post-merger and exhibit morphology similar to the Bullet Cluster viewed from an angle. From the separation between X-ray emitting gas and lensing mass in the Southern core, we derive a new and independent constraint on the self-interaction cross section of dark matter particles sigma/m <~ 3 \pm 1 cm^2 g^-1. In the Northwestern substructure, the gas, dark matter, and galaxy components have become separated by much larger distances. Most curiously, the `ghost' clump (primarily gas) leads the `dark' clump (primarily dark matter) by more than 150 kpc. We propose an enhanced `ram-pressure slingshot' scenario which may have yielded this reversal of components with such a large separation, but needs further confirmation by follow-up observations and numerical simulations. A secondary merger involves a second `bullet' clump in the North and an extremely `stripped' clump to the West. The latter appears to exhibit the largest separation between dark matter and X-ray emitting baryons detected to date in our sky.

The `Eagle' galaxy at a redshift of 0.77 is studied with the Oxford Short Wavelength Integral Field Spectrograph (SWIFT) and multi-wavelength data from the All-wavelength Extended Groth strip International Survey (AEGIS). It was chosen from AEGIS because of the bright and extended emission in its slit spectrum. Three dimensional kinematic maps of the Eagle reveal a gradient in velocity dispersion which spans 35-75 +/- 10 km/s and a rotation velocity of 25 +/- 5 km/s uncorrected for inclination. Hubble Space Telescope images suggest it is close to face-on. In comparison with galaxies from AEGIS at similar redshifts, the Eagle is extremely bright and blue in the rest-frame optical, highly star-forming, dominated by unobscured star-formation, and has a low metallicity for its size. This is consistent with its selection. The Eagle is likely undergoing a major merger and is caught in the early stage of a star-burst when it has not yet experienced metal enrichment or formed the mass of dust typically found in star-forming galaxies.

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

We report the X-ray detection of two z>1.4 infrared-selected galaxy clusters from the IRAC Shallow Cluster Survey (ISCS). We present new data from the Hubble Space Telescope and the W. M. Keck Observatory that spectroscopically confirm cluster ISCS J1432.4+3250 at z=1.49, the most distant of 18 confirmed z>1 clusters in the ISCS to date. We also present new spectroscopy for ISCS J1438.1+3414, previously reported at z = 1.41, and measure its dynamical mass. Clusters ISCS J1432.4+3250 and ISCS J1438.1+3414 are detected in 36ks and 143ks Chandra exposures at significances of 5.2 sigma and 9.7 sigma, from which we measure total masses of log(M_{200,Lx}/Msun) = 14.4 +/- 0.2 and 14.35^{+0.14}_{-0.11}, respectively. The consistency of the X-ray and dynamical properties of these high redshift clusters further demonstrates that the ISCS is robustly detecting massive clusters to at least z = 1.5.

We present the details and early results from a deep near-infrared survey utilising the NICMOS instrument on the Hubble Space Telescope centred around massive M_* > 10^11 M_0 galaxies at 1.7 < z < 2.9 found within the Great Observatories Origins Deep Survey (GOODS) fields. The GOODS NICMOS Survey (GNS) was designed to obtain deep F160W (H-band) imaging of 80 of these massive galaxies, as well as other colour selected objects such as Lyman-break drop-outs, BzK objects, Distant Red Galaxies, EROs, Spitzer Selected EROs, BX/BM galaxies, as well as sub-mm galaxies. We present in this paper details of the observations, our sample selection, as well as a description of features of the massive galaxies found within our survey fields. This includes: photometric redshifts, rest-frame colours, and stellar masses. We furthermore provide an analysis of the selection methods for finding massive galaxies at high redshifts, including colour selection, and how galaxy populations selected through different methods overlap. We find that a single colour selection method cannot locate all of the massive galaxies, with no one method finding more than 70 percent. We however find that the combination of these colour methods finds nearly all the massive galaxies, as selected by photometric redshifts with the exception of apparently rare blue massive galaxies. By investigating the rest-frame (U-B) vs. M_B diagram for these galaxies we furthermore show that there exists a bimodality in colour-magnitude space at z < 2, driven by stellar mass, such that the most massive galaxies are systematically red up to z~2.5, while lower mass galaxies tend to be blue. We also discuss the number densities for galaxies with stellar masses M_* > 10^11 M_0, whereby we find an increase of a factor of eight between z = 3 and z = 1.5, demonstrating that this is an epoch when massive galaxies establish most of their mass.

(Abridged) The tidal stirring model posits the formation of dSph galaxies via the tidal interactions between rotationally-supported dwarfs and MW-sized host galaxies. Using a set of collisionless N-body simulations, we investigate the efficiency of the tidal stirring mechanism. We explore a wide variety of dwarf orbital configurations and initial structures and demonstrate that in most cases the disky dwarfs experience significant mass loss and their stellar components undergo a dramatic morphological and dynamical transformation: from disks to bars and finally to pressure-supported spheroidal systems with kinematic and structural properties akin to those of the classic dSphs in the Local Group (LG). Our results suggest that such tidal transformations should be common occurrences within the currently favored cosmological paradigm and highlight the key factor responsible for an effective metamorphosis to be the strength of the tidal shocks at the pericenters of the orbit. We demonstrate that the combination of short orbital times and small pericenters, characteristic of dwarfs being accreted at high redshift, induces the strongest transformations. Our models also indicate that the transformation efficiency is affected significantly by the structure of the progenitor disky dwarfs. Lastly, we find that the dwarf remnants satisfy the relation Vmax = \sqrt{3} * sigma, where sigma is the 1D, central stellar velocity dispersion and Vmax is the maximum halo circular velocity, with intriguing implications for the missing satellites problem. Overall, we conclude that the action of tidal forces from the hosts constitutes a crucial evolutionary mechanism for shaping the nature of dwarf galaxies in environments such as that of the LG. Environmental processes of this type should thus be included as ingredients in models of dwarf galaxy formation and evolution.

[NIRSS is one of three concepts that contributed to the Wide-Field Infrared Survey Telescope (WFIRST) mission advocated by the Decadal Survey.] Operating beyond the reaches of the Earth's atmosphere, free of its limiting absorption and thermal background, the Near-Infrared Sky Surveyor (NIRSS) will deeply map the entire sky at near-infrared wavelengths, thereby enabling new and fundamental discoveries ranging from the identification of extrasolar planets to probing the reionization epoch by identifying thousands of quasars at z>10. NIRSS will directly address the NASA scientific objective of studying cosmic origins by using a 1.5-meter telescope to reach full-sky 0.2 uJy (25.6 mag AB) sensitivities in four passbands from 1 to 4 microns in a 4-yr mission. At the three shorter passbands (1 - 2.5 microns), the proposed depth is comparable to the deepest pencil-beam surveys done to date and is 3000 times more sensitive than the only previous all-sky near-infrared survey, 2MASS. At the longest passband (3.5 micron), which is not feasible from the ground, NIRSS will be 500 times more sensitive than WISE. NIRSS fills a pivotal gap in our knowledge of the celestial sphere, is a natural complement to WISE, and is well matched to the next generation of deep (0.1 uJy), wide-area (>2 pi ster), ground-based optical surveys (LSST and Pan-Starrs). With the high thermal backgrounds of ground-based infrared observations, a near-infrared full sky survey at sub-uJy sensitivity is only feasible from space.

Here we introduce PHAT, the PHoto-z Accuracy Testing programme, an international initiative to test and compare different methods of photo-z estimation. Two different test environments are set up, one (PHAT0) based on simulations to test the basic functionality of the different photo-z codes, and another one (PHAT1) based on data from the GOODS survey. The accuracy of the different methods is expressed and ranked by the global photo-z bias, scatter, and outlier rates. Most methods agree well on PHAT0 but produce photo-z scatters that can differ by up to a factor of two even in this idealised case. A larger spread in accuracy is found for PHAT1. Few methods benefit from the addition of mid-IR photometry. Remaining biases and systematic effects can be explained by shortcomings in the different template sets and the use of priors on the one hand and an insufficient training set on the other hand. Scatters of 4-8% in Delta_z/(1+z) were obtained, consistent with other studies. However, somewhat larger outlier rates (>7.5% with Delta_z/(1+z)>0.15; >4.5% after cleaning) are found for all codes. There is a general trend that empirical codes produce smaller biases than template-based codes. The systematic, quantitative comparison of different photo-z codes presented here is a snapshot of the current state-of-the-art of photo-z estimation and sets a standard for the assessment of photo-z accuracy in the future. The rather large outlier rates reported here for PHAT1 on real data should be investigated further since they are most probably also present (and possibly hidden) in many other studies. The test data sets are publicly available and can be used to compare new methods to established ones and help in guiding future photo-z method development. (abridged)

We use stellar masses, photometry, lensing, and velocity dispersions to investigate empirical correlations for the final sample of 73 early-type lens galaxies (ETGs) from the SLACS survey. The traditional correlations (Fundamental Plane [FP] and its projections) are consistent with those found for non-lens galaxies, supporting the thesis that SLACS lens galaxies are representative of massive ETGs. The addition of strong lensing estimates of the total mass allows us to gain further insights into their internal structure: i) the mean slope of the total mass density profile is <gamma'> = 2.078+/-0.027 with an intrinsic scatter of 0.16+/-0.02; ii) gamma' correlates with effective radius and central mass density, in the sense that denser galaxies have steeper profiles; iii) the dark matter fraction within reff/2 is a monotonically increasing function of galaxy mass and size; iv) the dimensional mass M_dim is proportional to the total mass, and both increase more rapidly than stellar mass M*; v) the Mass Plane (MP), obtained by replacing surface brightness with surface mass density in the FP, is found to be tighter and closer to the virial relation than the FP and the M*P, indicating that the scatter of those relations is dominated by stellar population effects; vi) we construct the Fundamental Hyper-Plane by adding stellar masses to the MP and find the M* coefficient to be consistent with zero and no residual intrinsic scatter. Our results demonstrate that the dynamical structure of ETGs is not scale invariant and that it is fully specified by the total mass, r_eff, and sigma. Although the basic trends can be explained qualitatively in terms of varying star formation efficiency as a function of halo mass and as the result of dry and wet mergers, reproducing quantitatively the observed correlations and their tightness may be a significant challenge for galaxy formation models.

We use the multi-epoch, mid-infrared Spitzer Deep, Wide-Field Survey to investigate the variability of 474,179 objects in 8.1 deg^2 of the NDWFS Bootes field. We perform a Difference Image Analysis of the four available epochs between 2004 and 2008, focusing on the deeper 3.6 and 4.5 micron bands. We find that 1.1% of the studied sample meet our standard selection criteria for being classed as a variable source. We require that the 3.6 and 4.5 micron light-curves are strongly correlated (r>0.8) and that their joint variance exceeds that for all sources with the same magnitude by 2 sigma. We then examine the mid-IR colors of the variable sources and match them with X-ray sources from the XBootes survey, radio catalogs, 24 micron-selected AGN candidates, and spectroscopically identified AGNs from the AGN and Galaxy Evolution Survey (AGES). Based on their mid-IR colors, most of the variable sources are AGNs (76%), with smaller contributions from stars (11%), galaxies (6%), and unclassified objects. Most of the stellar, galaxy and unclassified sources are false positives. For our standard selection criteria, 11-12% of the mid-IR counterparts to X-ray sources, 24 micron-selected AGN candidates and spectroscopically identified AGNs show variability. Mid-IR AGN variability can be well described by a single power-law structure function with a power-law index of 0.5 at both 3.6 and 4.5 microns, and an amplitude of 0.1 mag on rest-frame time scales of 2 years. The variability amplitude is higher for shorter rest-frame wavelengths and lower luminosities. (Abridged)

We present a strong lensing mass model of Abell 1689 which resolves substructures ~25 kpc across (including about ten individual galaxy subhalos) within the central ~400 kpc diameter. We achieve this resolution by perfectly reproducing the observed (strongly lensed) input positions of 168 multiple images of 55 knots residing within 135 images of 42 galaxies. Our model makes no assumptions about light tracing mass, yet we reproduce the brightest visible structures with some slight deviations. A1689 remains one of the strongest known lenses on the sky, with an Einstein radius of RE = 47.0" +/- 1.2" (143 +3/-4 kpc) for a lensed source at zs = 2. We find a single NFW or Sersic prole yields a good fit simultaneously (with only slight tension) to both our strong lensing (SL) mass model and published weak lensing (WL) measurements at larger radius (out to the virial radius). According to this NFW fit, A1689 has a mass of Mvir = 2.0 +0.5/-0.3 x 10^15 Msun / h70 (M200 = 1.8 +0.4/-0.3 x 10^15 Msun / h70) within the virial radius rvir = 3.0 +/- 0.2 Mpc / h70 (r200 = 2.4 +0.1/-0.2 Mpc / h70), and a central concentration cvir = 11.5 +1.5/-1.4 (c200 = 9.2 +/- 1.2). Our SL model prefers slightly higher concentrations than previous SL models, bringing our SL+WL constraints in line with other recent derivations. Our results support those of previous studies which find A1689 has either an anomalously large concentration or significant extra mass along the line of sight (perhaps in part due to triaxiality). If clusters are generally found to have higher concentrations than realized in simulations, this could indicate they formed earlier, perhaps as a result of early dark energy.

Scaling relations of observed galaxy cluster properties are useful tools for constraining cosmological parameters as well as cluster formation histories. One of the key cosmological parameters, sigma8, is constrained using observed clusters of galaxies, although current estimates of sigma8 from the scaling relations of dynamically relaxed galaxy clusters are limited by the large scatter in the observed cluster mass-temperature (M-T) relation. With a sample of eight strong lensing clusters at 0.3 < z <0.8, we find that the observed cluster concentration-mass relation can be used to reduce the M-T scatter by a factor of 6. Typically only relaxed clusters are used to estimate sigma8, but combining the cluster concentration-mass relation with the M-T relation enables the inclusion of unrelaxed clusters as well. Thus, the resultant gains in the accuracy of sigma8 measurements from clusters are twofold: the errors on sigma8 are reduced and the cluster sample size is increased. Therefore, the statistics on sigma8 determination from clusters are greatly improved by the inclusion of unrelaxed clusters. Exploring cluster scaling relations further, we find that the correlation between brightest cluster galaxy (BCG) luminosity and cluster mass offers insight into the assembly histories of clusters. We find preliminary evidence for a steeper BCG luminosity - cluster mass relation for strong lensing clusters than the general cluster population, hinting that strong lensing clusters may have had more active merging histories.

We present near-IR images, obtained with the Hubble Space Telescope (HST) and the WFC3/IR camera, of six passive and massive galaxies at redshift 1.3<z<2.4 (SSFR<10^{-2} Gyr^{-1}; stellar mass M~10^{11} M_{sun}), selected from the Great Observatories Origins Deep Survey (GOODS). These images, which have a spatial resolution of ~1.5 kpc, provide the deepest view of the optical rest-frame morphology of such systems to date. We find that the light profile of these galaxies is regular and well described by a Sersic model with index typical of today's spheroids. Their size, however, is generally much smaller than today's early types of similar stellar mass, with four out of six galaxies having r_e ~ 1 kpc or less, in quantitative agreement with previous similar measures made at rest-frame UV wavelengths. The images reach limiting surface brightness mu~26.5 mag arcsec^{-2} in the F160W bandpass; yet, there is no evidence of a faint halo in the galaxies of our sample, even in their stacked image. We also find that these galaxies have very weak "morphological k-correction" between the rest-frame UV (from the ACS z-band), and the rest--frame optical (WFC3 H-band): the Sersic index, physical size and overall morphology are independent or only mildly dependent on the wavelength, within the errors.

We present the current photometric dataset for the Sloan Lens ACS (SLACS) Survey, including HST photometry from ACS, WFPC2, and NICMOS. These data have enabled the confirmation of an additional 15 grade `A' (certain) lens systems, bringing the number of SLACS grade `A' lenses to 85; including 13 grade `B' (likely) systems, SLACS has identified nearly 100 lenses and lens candidates. Approximately 80% of the grade `A' systems have elliptical morphologies while ~10% show spiral structure; the remaining lenses have lenticular morphologies. Spectroscopic redshifts for the lens and source are available for every system, making SLACS the largest homogeneous dataset of galaxy-scale lenses to date. We have developed a novel Bayesian stellar population analysis code to determine robust stellar masses with accurate error estimates. We apply this code to deep, high-resolution HST imaging and determine stellar masses with typical statistical errors of 0.1 dex; we find that these stellar masses are unbiased compared to estimates obtained using SDSS photometry, provided that informative priors are used. The stellar masses range from 10^10.5 to 10^11.8 M$_\odot$ and the typical stellar mass fraction within the Einstein radius is 0.4, assuming a Chabrier IMF. The ensemble properties of the SLACS lens galaxies, e.g. stellar masses and projected ellipticities, appear to be indistinguishable from other SDSS galaxies with similar stellar velocity dispersions. This further supports that SLACS lenses are representative of the overall population of massive early-type galaxies with M* >~ 10^11 M$_\odot$, and are therefore an ideal dataset to investigate the kpc-scale distribution of luminous and dark matter in galaxies out to z ~ 0.5.

We report the discovery of a multiply lensed Ly Alpha (Lya) emitter at z = 3.90 behind the massive galaxy cluster WARPS J1415.1+3612 at z = 1.026. Images taken by the Hubble Space Telescope(HST) using ACS reveal a complex lensing system that produces a prominent, highly magnified arc and a triplet of smaller arcs grouped tightly around a spectroscopically confirmed cluster member. Spectroscopic observations using FOCAS on Subaru confirm strong Lya emission in the source galaxy and provide redshifts for more than 21 cluster members, from which we obtain a velocity dispersion of 807+/-185 km/s. Assuming a singular isothermal sphere profile, the mass within the Einstein ring (7.13+/-0.38") corresponds to a central velocity dispersion of 686+15-19 km/s for the cluster, consistent with the value estimated from cluster member redshifts. Our mass profile estimate from combining strong lensing and dynamical analyses is in good agreement with both X-ray and weak lensing results.

Future large ensembles of time delay lenses have the potential to provide interesting cosmological constraints complementary to those of other methods. In a flat universe with constant w including a Planck prior, LSST time delay measurements for ~4,000 lenses should constrain the local Hubble constant h to ~0.007 (~1%), Omega_DE to ~0.005, and w to ~0.026 (all 1-sigma precisions). Similar constraints could be obtained by a dedicated gravitational lens observatory (OMEGA) which would obtain precise time delay and mass model measurements for ~100 well-studied lenses. We compare these constraints (as well as those for a general cosmology) to the "optimistic Stage IV" constraints expected from weak lensing, supernovae, baryon acoustic oscillations, and cluster counts, as calculated by the Dark Energy Task Force. Time delays yield a modest constraint on a time-varying w(z), with the best constraint on w(z) at the "pivot redshift" of z ~ 0.31. Our Fisher matrix calculation is provided to allow time delay constraints to be easily compared to and combined with constraints from other experiments. We also show how cosmological constraining power varies as a function of numbers of lenses, lens model uncertainty, time delay precision, redshift precision, and the ratio of four-image to two-image lenses.

Based on 58 SLACS strong-lens early-type galaxies with direct total-mass and stellar-velocity dispersion measurements, we find that inside one effective radius massive elliptical galaxies with M_eff >= 3x10^10 M_sun are well-approximated by a power-law ellipsoid with an average logaritmic density slope of <gamma'_LD> = -dlog(rho_tot)/dlog(r)=2.085^{+0.025}_{-0.018} (random error on mean) for isotropic orbits with beta_r=0, +-0.1 (syst.) and sigma_gamma' <= 0.20^{+0.04}_{-0.02} intrinsic scatter (all errors indicate the 68 percent CL). We find no correlation of gamma'_LD with galaxy mass (M_eff), rescaled radius (i.e. R_einst/R_eff) or redshift, despite intrinsic differences in density-slope between galaxies. Based on scaling relations, the average logarithmic density slope can be derived in an alternative manner, fully independent from dynamics, yielding <gamma'_SR>=1.959 +- 0.077. Agreement between the two values is reached for <beta_r> =0.45 +- 0.25, consistent with mild radial anisotropy. This agreement supports the robustness of our results, despite the increase in mass-to-light ratio with total galaxy mass: M_eff ~ L_{V,eff}^(1.363+-0.056). We conclude that massive early-type galaxies are structurally close-to homologous with close-to isothermal total density profiles (<=10 percent intrinsic scatter) and have at most some mild radial anisotropy. Our results provide new observational limits on galaxy formation and evolution scenarios, covering four Gyr look-back time.

The Spitzer Deep, Wide-Field Survey (SDWFS) is a four-epoch infrared survey of ten square degrees in the Bootes field of the NOAO Deep Wide-Field Survey using the IRAC instrument on the Spitzer Space Telescope. SDWFS, a Cycle four Spitzer Legacy project, occupies a unique position in the area-depth survey space defined by other Spitzer surveys. The four epochs that make up SDWFS permit -- for the first time -- the selection of infrared-variable and high proper motion objects over a wide field on timescales of years. Because of its large survey volume, SDWFS is sensitive to galaxies out to z~3 with relatively little impact from cosmic variance for all but the richest systems. The SDWFS datasets will thus be especially useful for characterizing galaxy evolution beyond z~1.5. This paper explains the SDWFS observing strategy and data processing, presents the SDWFS mosaics and source catalogs, and discusses some early scientific findings. The publicly-released, full-depth catalogs contain 6.78, 5.23, 1.20, and 0.96 x 10e5 distinct sources detected to the average 5-sigma, 4" diameter, aperture-corrected limits of 19.77, 18.83, 16.50, and 15.82 Vega mag at 3.6, 4.5, 5.8, and 8.0 micron, respectively. The SDWFS number counts and color-color distribution are consistent with other, earlier Spitzer surveys. At the 6 min integration time of the SDWFS IRAC imaging, more than 50% of isolated FIRST radio sources and more than 80% of on-axis XBootes sources are detected out to 8.0 micron. Finally, we present the four highest proper motion IRAC-selected sources identified from the multi-epoch imaging, two of which are likely field brown dwarfs of mid-T spectral class.

We have undertaken a systematic search for candidate supernovae from high-redshift Population III stars in a field that has been observed with repeated imaging on a cadence of 2-3 weeks over a 2.2 year baseline, the Spitzer/IRAC Dark Field. The individual epochs reach a typical 5-sigma depth of 1 uJy in IRAC Channel 1 (3.6 um). Requiring a minimum of four epochs coverage, the total effective area searched is 214 sq arcminutes. The unprecedented depth and multi-epochal nature of these data make it ideal for a first foray to detect transient objects which may be candidate luminous Pair Instability Supernovae from the primordial-metallicity first stars. The search was conducted over a broad range of timescales, allowing for different durations of the putative candidates' light curve plateau phases. All candidates were vetted by inspection of the Spitzer imaging data, as well as deep HST/ACS F814W imaging available over the full field. While many resolved-source objects were found with Spitzer variability, no transient objects could plausibly be identified as high-redshift supernovae candidates. The resulting 95% confidence level upper limit is 23/sq deg/yr, for sources with plateau timescales under 400/(1+z) days and brightnesses above ~1 uJy.

We show that dark matter substructure in galaxy-scale halos perturbs the time delays between images in strong gravitational lens systems. The variance of the effect depends on the subhalo mass function, scaling as the product of the substructure mass fraction and a characteristic mass of subhalos (namely <m^2>/<m>). Time delay perturbations therefore complement gravitational lens flux ratio anomalies and astrometric perturbations by measuring a different moment of the subhalo mass function. Unlike flux ratio anomalies, "time delay millilensing" is unaffected by dust extinction or stellar microlensing in the lens galaxy. Furthermore, we show that time delay ratios are immune to the radial profile degeneracy that usually plagues lens modeling. We lay out a mathematical theory of time delay perturbations and find it to be tractable and attractive. We predict that in "cusp" lenses with close triplets of images, substructure may change the arrival-time order of the images (compared with smooth models). We discuss the possibility that this effect has already been observed in RX J1131-1231.

The cosmic dark ages are the mysterious epoch during which the pristine gas began to condense and ultimately form the first stars. Although these beginnings have long been a topic of theoretical interest, technology has only recently allowed the beginnings of observational insight into this epoch. Many questions surround the formation of stars in metal-free gas and the history of the build-up of metals in the intergalactic medium: (1) What were the properties of the first stellar and galactic sources to form in pristine (metal-free) gas? (2) When did the epoch of Population III (metal-free) star formation take place and how long did it last? (3) Was the stellar initial mass function dramatically different for the first stars and galaxies? These questions are all active areas of theoretical research. However, new observational constraints via the direct detection of Population III star formation are vital to making progress in answering the broader questions surrounding how galaxies formed and how the cosmological properties of the universe have affected the objects it contains.

Whereas considerable effort has been afforded in understanding the properties of galaxies, a full physical picture, connecting their baryonic and dark-matter content, super-massive black holes, and (metric) theories of gravity, is still ill-defined. Strong gravitational lensing furnishes a powerful method to probe gravity in the central regions of galaxies. It can (1) provide a unique detection-channel of dark-matter substructure beyond the local galaxy group, (2) constrain dark-matter physics, complementary to direct-detection experiments, as well as metric theories of gravity, (3) probe central super-massive black holes, and (4) provide crucial insight into galaxy formation processes from the dark matter point of view, independently of the nature and state of dark matter. To seriously address the above questions, a considerable increase in the number of strong gravitational-lens systems is required. In the timeframe 2010-2020, a staged approach with radio (e.g. EVLA, e-MERLIN, LOFAR, SKA phase-I) and optical (e.g. LSST and JDEM) instruments can provide 10^(2-4) new lenses, and up to 10^(4-6) new lens systems from SKA/LSST/JDEM all-sky surveys around ~2020. Follow-up imaging of (radio) lenses is necessary with moderate ground/space-based optical-IR telescopes and with 30-50m telescopes for spectroscopy (e.g. TMT, GMT, ELT). To answer these fundamental questions through strong gravitational lensing, a strong investment in large radio and optical-IR facilities is therefore critical in the coming decade. In particular, only large-scale radio lens surveys (e.g. with SKA) provide the large numbers of high-resolution and high-fidelity images of lenses needed for SMBH and flux-ratio anomaly studies.

There is a vast menagerie of plausible candidates for the constituents of dark matter, both within and beyond extensions of the Standard Model of particle physics. Each of these candidates may have scattering (and other) cross section properties that are consistent with the dark matter abundance, BBN, and the most scales in the matter power spectrum; but which may have vastly different behavior at sub-galactic "cutoff" scales, below which dark matter density fluctuations are smoothed out. The only way to quantitatively measure the power spectrum behavior at sub-galactic scales at distances beyond the local universe, and indeed over cosmic time, is through probes available in multiply imaged strong gravitational lenses. Gravitational potential perturbations by dark matter substructure encode information in the observed relative magnifications, positions, and time delays in a strong lens. Each of these is sensitive to a different moment of the substructure mass function and to different effective mass ranges of the substructure. The time delay perturbations, in particular, are proving to be largely immune to the degeneracies and systematic uncertainties that have impacted exploitation of strong lenses for such studies. There is great potential for a coordinated theoretical and observational effort to enable a sophisticated exploitation of strong gravitational lenses as direct probes of dark matter properties. This opportunity motivates this white paper, and drives the need for: a) strong support of the theoretical work necessary to understand all astrophysical consequences for different dark matter candidates; and b) tailored observational campaigns, and even a fully dedicated mission, to obtain the requisite data.

The Cold Dark Matter theory of gravitationally-driven hierarchical structure formation has earned its status as a paradigm by explaining the distribution of matter over large spans of cosmic distance and time. However, its central tenet, that most of the matter in the universe is dark and exotic, is still unproven; the dark matter hypothesis is sufficiently audacious as to continue to warrant a diverse battery of tests. While local searches for dark matter particles or their annihilation signals could prove the existence of the substance itself, studies of cosmological dark matter in situ are vital to fully understand its role in structure formation and evolution. We argue that gravitational lensing provides the cleanest and farthest-reaching probe of dark matter in the universe, which can be combined with other observational techniques to answer the most challenging and exciting questions that will drive the subject in the next decade: What is the distribution of mass on sub-galactic scales? How do galaxy disks form and bulges grow in dark matter halos? How accurate are CDM predictions of halo structure? Can we distinguish between a need for a new substance (dark matter) and a need for new physics (departures from General Relativity)? What is the dark matter made of anyway? We propose that the central tool in this program should be a wide-field optical imaging survey, whose true value is realized with support in the form of high-resolution, cadenced optical/infra-red imaging, and massive-throughput optical spectroscopy.

(Brief Summary) What is the total radiative content of the Universe since the epoch of recombination? The extragalactic background light (EBL) spectrum captures the redshifted energy released from the first stellar objects, protogalaxies, and galaxies throughout cosmic history. Yet, we have not determined the brightness of the extragalactic sky from UV/optical to far-infrared wavelengths with sufficient accuracy to establish the radiative content of the Universe to better than an order of magnitude. Among many science topics, an accurate measurement of the EBL spectrum from optical to far-IR wavelengths, will address: What is the total energy released by stellar nucleosynthesis over cosmic history? Was significant energy released by non-stellar processes? Is there a diffuse component to the EBL anywhere from optical to sub-millimeter? When did first stars appear and how luminous was the reionization epoch? Absolute optical to mid-IR EBL spectrum to an astrophysically interesting accuracy can be established by wide field imagingat a distance of 5 AU or above the ecliptic plane where the zodiacal foreground is reduced by more than two orders of magnitude.

We present the surface density of luminous active galactic nuclei (AGN) associated with a uniformly selected galaxy cluster sample identified in the 8.5 square degree Bootes field of the NOAO Deep Wide-Field Survey. The clusters are distributed over a large range of redshift (0 < z < 1.5) and we identify AGN using three different selection criteria: mid-IR color, radio luminosity, and X-ray luminosity. Relative to the field, we note a clear overdensity of the number of AGN within 0.5 Mpc of the cluster centers at z > 0.5. The amplitude of this AGN overdensity increases with redshift. Although there are significant differences between the AGN populations probed by each selection technique, the rise in cluster AGN surface density generally increases more steeply than that of field quasars. In particular, X-ray selected AGN are at least three times more prevalent in clusters at 1 < z < 1.5 compared to clusters at 0.5 < z < 1. This effect is stronger than can be explained by the evolving median richness of our cluster sample. We thus confirm the existence of a Butcher-Oemler type effect for AGN in galaxy clusters, with the number of AGN in clusters increasing with redshift.

We study the relation between the internal structure of early-type galaxies and their environment using 70 strong gravitational lenses from the Sloan ACS Lens Survey. The Sloan database is used to determine two measures of overdensity of galaxies around each lens: the projected number density of galaxies inside the tenth nearest neighbor (\Sigma_{10}) and within a cone of radius one h^{-1} Mpc (D_1). Our main results are: 1) The average overdensity is somewhat larger than unity, consistent with lenses preferring overdense environments as expected for massive early-type galaxies (12/70 lenses are in known groups/clusters). 2) The distribution of overdensities is indistinguishable from that of "twin" non-lens galaxies selected from SDSS to have the same redshift and stellar velocity dispersion \sigma_*. Thus, within our errors, lens galaxies are an unbiased population, and the SLACS results can be generalized to the overall population of early-type galaxies. 3) Typical contributions from external mass distribution are no more than a few per cent, reaching 10-20% (~0.05-0.10 external convergence) only in the most extreme overdensities. 4) No significant correlation between overdensity and slope of the mass density profile of the lens is found. 5) Satellite galaxies (those with a more luminous companion) have marginally steeper mass density profiles than central galaxies (those without). This result suggests that tidal stripping may affect the mass structure of early-type galaxies down to kpc scales probed by strong lensing, when they fall into larger structures [ABRIDGED].

(Abridged) We conduct a series of high-resolution, dissipationless N-body simulations to investigate the cumulative effect of substructure mergers onto thin disk galaxies in the context of the LCDM paradigm of structure formation. Our simulation campaign is based on a hybrid approach. Substructure properties are culled directly from cosmological simulations of galaxy-sized cold dark matter (CDM) halos. In contrast to what can be inferred from statistics of the present-day substructure populations, accretions of massive subhalos onto the central regions of host halos, where the galactic disk resides, since z~1 should be common occurrences. One host halo merger history is subsequently used to seed controlled numerical experiments of repeated satellite impacts on an initially-thin Milky Way-type disk galaxy. We show that these accretion events produce several distinctive observational signatures in the stellar disk including: a ring-like feature in the outskirts; a significant flare; a central bar; and faint filamentary structures that (spuriously) resemble tidal streams. The final distribution of disk stars exhibits a complex vertical structure that is well-described by a standard ``thin-thick'' disk decomposition. We conclude that satellite-disk encounters of the kind expected in LCDM models can induce morphological features in galactic disks that are similar to those being discovered in the Milky Way, M31, and in other disk galaxies. These results highlight the significant role of CDM substructure in setting the structure of disk galaxies and driving galaxy evolution. Upcoming galactic structure surveys and astrometric satellites may be able to distinguish between competing cosmological models by testing whether the detailed structure of galactic disks is as excited as predicted by the CDM paradigm.

(abridged) The Observatory for Multi-Epoch Gravitational Lens Astrophysics (OMEGA) is a mission concept for a 1.5-m near-UV through near-IR space observatory that will be dedicated to frequent imaging and spectroscopic monitoring of ~100 multiply-imaged active galactic nuclei over the whole sky. Using wavelength-tailored dichroics with extremely high transmittance, efficient imaging in six channels will be done simultaneously during each visit to each target. The separate spectroscopic mode, engaged through a flip-in mirror, uses an image slicer spectrograph. After a period of many visits to all targets, the resulting multidimensional movies can then be analyzed to a) measure the mass function of dark matter substructure; b) measure precise masses of the accreting black holes as well as the structure of their accretion disks and their environments over several decades of physical scale; and c) measure a combination of the local Hubble expansion and cosmological distances to unprecedented precision.

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

We present the definitive data for the full sample of 131 strong gravitational lens candidates observed with the Advanced Camera for Surveys (ACS) aboard the Hubble Space Telescope by the Sloan Lens ACS (SLACS) Survey. All targets were selected for higher-redshift emission lines and lower-redshift continuum in a single Sloan Digital Sky Survey (SDSS) spectrum. The foreground galaxies are primarily of early-type morphology, with redshifts from approximately 0.05 to 0.5 and velocity dispersions from 160 km/s to 400 km/s; the faint background emission-line galaxies have redshifts ranging from about 0.2 to 1.2. We confirm 70 systems showing clear evidence of multiple imaging of the background galaxy by the foreground galaxy, as well as an additional 19 systems with probable multiple imaging. For 63 clear lensing systems, we present singular isothermal ellipsoid and light-traces-mass gravitational lens models fitted to the ACS imaging data. These strong-lensing mass measurements are supplemented by magnitudes and effective radii measured from ACS surface-brightness photometry and redshifts and velocity dispersions measured from SDSS spectroscopy. These data constitute a unique resource for the quantitative study of the inter-relations between mass, light, and kinematics in massive early-type galaxies. We show that the SLACS lens sample is statistically consistent with being drawn at random from a parent sample of SDSS galaxies with comparable spectroscopic parameters and effective radii, suggesting that the results of SLACS analyses can be generalized to the massive early-type population.

Lens modeling is the key to successful and meaningful automated strong galaxy-scale gravitational lens detection. We have implemented a lens-modeling "robot" that treats every bright red galaxy (BRG) in a large imaging survey as a potential gravitational lens system. Using a simple model optimized for "typical" galaxy-scale lenses, we generate four assessments of model quality that are used in an automated classification. The robot infers the lens classification parameter H that a human would have assigned; the inference is performed using a probability distribution generated from a human-classified training set, including realistic simulated lenses and known false positives drawn from the HST/EGS survey. We compute the expected purity, completeness and rejection rate, and find that these can be optimized for a particular application by changing the prior probability distribution for H, equivalent to defining the robot's "character." Adopting a realistic prior based on the known abundance of lenses, we find that a lens sample may be generated that is ~100% pure, but only ~20% complete. This shortfall is due primarily to the over-simplicity of the lens model. With a more optimistic robot, ~90% completeness can be achieved while rejecting ~90% of the candidate objects. The remaining candidates must be classified by human inspectors. We are able to classify lens candidates by eye at a rate of a few seconds per system, suggesting that a future 1000 square degree imaging survey containing 10^7 BRGs, and some 10^4 lenses, could be successfully, and reproducibly, searched in a modest amount of time. [Abridged]

We have identified 335 galaxy cluster and group candidates spanning 0<z<2, using a 4.5um selected sample of galaxies in a 7.25 deg^2 region in the Spitzer/IRAC Shallow Survey. Using full redshift probability distributions for all galaxies, clusters were identified as 3-dimensional overdensities using a wavelet algorithm. To date 12 clusters at z>1, and over 60 at z < 0.5 have been spectroscopically confirmed. The mean I-[3.6] color for cluster galaxies up to z~1 is well matched by a z_f=3 passively evolving model. At z>1, a wider range of formation histories is needed, but higher formation redshifts (i.e. z_f \ge 4-5) are favored for most clusters. The cluster autocorrelation function, measured for the first time out to z=1.5, is found not to have evolved over the last 10 Gyr, in agreement with the prediction from LCDM. The average mass of the IRAC Shallow Cluster Survey sample, inferred from its clustering, is ~ 10^14 Msun.

We have identified 335 galaxy cluster and group candidates, 106 of which are at z > 1, using a 4.5 um selected sample of objects from a 7.25 deg^2 region in the Spitzer Infrared Array Camera (IRAC) Shallow Survey. Clusters were identified as 3-dimensional overdensities using a wavelet algorithm, based on photometric redshift probability distributions derived from IRAC and NOAO Deep Wide-Field Survey data. We estimate only ~10% of the detections are spurious. To date 12 of the z > 1 candidates have been confirmed spectroscopically, at redshifts from 1.06 to 1.41. Velocity dispersions of ~750 km/s for two of these argue for total cluster masses well above 10^14 M_sun, as does the mass estimated from the rest frame near infrared stellar luminosity. Although not selected to contain a red sequence, some evidence for red sequences is present in the spectroscopically confirmed clusters, and brighter galaxies are systematically redder than the mean galaxy color in clusters at all redshifts. The mean I - [3.6] color for cluster galaxies up to z ~ 1 is well matched by a passively evolving model in which stars are formed in a 0.1 Gyr burst starting at redshift z_f = 3. At z > 1, a wider range of formation histories is needed, but higher formation redshifts (i.e. z_f > 3) are favored for most clusters.

We report the discovery of two concentric Einstein rings around the gravitational lens SDSSJ0946+1006, as part of the Sloan Lens ACS Survey. The main lens is at redshift zl=0.222, while the inner ring (1) is at zs1=0.609 and Einstein radius $Re_1=1.43\pm0.01"$. The wider image separation ($Re_2=2.07\pm 0.02"$) of the outer ring (2) implies that it is at higher redshift. Its detection in the F814W filter implies zs2<6.9. The configuration can be well described by a total density profile $\rho_{tot} ~ r^-g'$ with $g'=2.00\pm0.03$ and velocity dispersion $\sigma_{SIE}=287\pm5\kms$. [...] We consider whether this configuration can be used to constrain cosmological parameters exploiting angular distance ratios entering the lens equations. Constraints for SDSSJ0946+1006, are uninteresting due to the sub-optimal lens and source redshifts. We then consider the perturbing effect of the mass associated with Ring 1 building a double lens plane compound lens model. This introduces minor changes to the mass of the main lens and allows to estimate the mass of Ring 1 $(\sigma_{SIE,s1}=94\pm30\kms)$. We examine the prospects of doing cosmography with a sample of 50 double lenses, expected from future space based surveys such as DUNE or JDEM. Taking full account of the model uncertainties, such a sample could be used to measure $\Omega_m$ and $w$ with 10% accuracy, for a flat cosmology.

We present the galaxy cluster autocorrelation function of 277 galaxy cluster candidates with 0.25 \le z \le 1.5 in a 7 deg^2 area of the IRAC Shallow Cluster Survey. We find strong clustering throughout our galaxy cluster sample, as expected for these massive structures. Specifically, at <z> = 0.5 we find a correlation length of r_0 = 17.40^{+3.98}_{-3.10} h^-1 Mpc, in excellent agreement with the Las Campanas Distant Cluster Survey, the only other non-local measurement. At higher redshift, <z> = 1, we find that strong clustering persists, with a correlation length of r_0=19.14^{+5.65}_{-4.56} h^-1 Mpc. A comparison with high resolution cosmological simulations indicates these are clusters with halo masses of \sim 10^{14} Msun, a result supported by estimates of dynamical mass for a subset of the sample. In a stable clustering picture, these clusters will evolve into massive (10^{15} Msun) clusters by the present day.

Aims. We present the full data set of the spectroscopic campaign of the ESO/GOODS program in the GOODS-South field, obtained with the FORS2 spectrograph at the ESO/VLT. Method. Objects were selected as candidates for VLT/FORS2 observations primarily based on the expectation that the detection and measurement of their spectral features would benefit from the high throughput and spectral resolution of FORS2. The reliability of the redshift estimates is assessed using the redshift-magnitude and color-redshift diagrams, and comparing the results with public data. Results. Including the third part of the spectroscopic campaign (12 masks) to the previous work (26 masks, Vanzella et al. 2005, 2006), 1715 spectra of 1225 individual targets have been analyzed. The actual spectroscopic catalog provides 887 redshift determinations. The typical redshift uncertainty is estimated to be sigma(z) ~ 0.001. Galaxies have been selected adopting different color criteria and using photometric redshifts. The resulting redshift distribution typically spans two domains: from z=0.5 to 2 and z=3.5 to 6.3. The reduced spectra and the derived redshifts are released to the community through the ESO web page http://www.eso.org/science/goods/

We combine high-resolution images in four optical/infra-red bands, obtained with the laser guide star adaptive optics system on the Keck Telescope and with the Hubble Space Telescope, to study the gravitational lens system SDSSJ0737+3216 (lens redshift 0.3223, source redshift 0.5812). We show that (under favorable observing conditions) ground-based images are comparable to those obtained with HST in terms of precision in the determination of the parameters of both the lens mass distribution and the background source. We also quantify the systematic errors associated with both the incomplete knowledge of the PSF, and the uncertain process of lens galaxy light removal, and find that similar accuracy can be achieved with Keck LGSAO as with HST. We then exploit this well-calibrated combination of optical and gravitational telescopes to perform a multi-wavelength study of the source galaxy at 0.01" effective resolution. We find the Sersic index to be indicative of a disk-like object, but the measured half-light radius (0.59+-0.007+-0.1 kpc) and stellar mass (2.0+-1.0+-0.8e9Msun) place it more than three sigma away from the local disk size-mass relation. The SDSSJ0737+3216 source has the characteristics of the most compact faint blue galaxies studied, and has comparable size and mass to dwarf early-type galaxies in the local universe. With the aid of gravitational telescopes to measure individual objects' brightness profiles to 10% accuracy, the study of the high-redshift size-mass relation may be extended by an order of magnitude or more beyond existing surveys at the low-mass end, thus providing a new observational test of galaxy formation models.

We combine strong-lensing masses with SDSS stellar velocity dispersions and HST-ACS effective (half-light) radii for 36 lens galaxies from the Sloan Lens ACS (SLACS) Survey to study the mass dependence of mass-dynamical structure in early-type galaxies. We find that over a 180--390 km/s range in velocity dispersion, structure is independent of lensing mass to within 5%. This result suggests a systematic variation in the total (i.e., luminous plus dark matter) mass-to-light ratio as the origin of the tilt of the fundamental plane (FP) scaling relationship between galaxy size, velocity dispersion, and surface brightness. We construct the FP of the lens sample, which we find to be consistent with the FP of the parent SDSS early-type galaxy population, and present the first observational correlation between mass-to-light ratio and residuals about the FP. Finally, we re-formulate the FP in terms of surface mass density rather than surface brightness. By removing the complexities of stellar-population effects, this mass-plane formulation will facilitate comparison to numerical simulations and possible use as a cosmological distance indicator.

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

We present a catalog of spectro-photometric redshifts for 1308 galaxies from the GRism ACS Program for Extragalactic Science (GRAPES) observations with the Hubble Space Telescope. These low-resolution spectra between 6000 A and 9500 A are supplemented with U, J, H, and Ks from various facilities, resulting in redshifts computed with ~40 spectral bins per galaxy. For 81 galaxies between 0.5<z<1.5 with spectroscopic redshifts, the standard deviation in the fractional error in (1+z) is 0.046. With this catalog, we compute the B-band luminosity function in this redshift range from 72 galaxies. Owing to the depth of the GRAPES survey, we are able to accurately constrain the faint-end slope by going to M_B~-18 mag at 0.8<z<1.2, nearly two magnitudes fainter than previous studies. The faint-end slope is alpha=-1.32+-0.07. When compared to numerous published values at various redshifts, we find strong evidence for a steepening of the faint-end slope with redshift which is expected in the hierarchical formation scenario of galaxies.

We present redshifts for 115 emission line objects in the Hubble Ultra Deep Field (HUDF) identified through the GRism ACS Program for Extragalactic Science (GRAPES) project using the slitless grism spectroscopy mode of the ACS Camera on the Hubble Space Telescope (HST). The sample was selected by an emission line search on all extracted 1-dimensional GRAPES spectra. We identify the emission lines using line wavelength ratios where multiple lines are detected in the grism wavelength range (5800A < lambda < 9600A), and using photometric redshift information where multiple lines are unavailable. We then derive redshifts using the identified lines. Our redshifts are accurate to delta(z) = 0.009, based on both statistical uncertainty estimates and comparison with published ground-based spectra. Over 40% of our sample is fainter than typical magnitude limits for ground-based spectroscopy (with i_{AB}>25 mag). Such emission lines would likely remain undiscovered without our deep survey. The emission line objects fall into 3 categories: 1) Most are low to moderate redshift galaxies (0 < z < 2), including many actively star forming galaxies with strong HII regions; 2) 9 are high redshift (4 < z < 7) Lyman-alpha emitters; and 3) at least 3 are candidate AGNs.

We derive Fundamental Plane parameters of 15 early-type lens galaxies identified by the Sloan Lens ACS (SLACS) Survey. The size of the sample allows us to investigate for the first time the distribution of lens galaxies in the FP space. After correcting for evolution, we find that lens galaxies occupy a subset of the local FP. The edge-on projection (approximately M vs M/L) is indistinguishable from that of normal early-type galaxies. However -- within the fundamental plane -- the lens galaxies appear to concentrate at the edge of the region populated by normal early-type galaxies. We show that this is a result of our selection procedure (approximately velocity dispersion sigma>240km/s). We conclude that SLACS lenses are a fair sample of high velocity dispersion early-type galaxies. By comparing the central stellar velocity dispersion that of the best fit lens model, we find <f_SIE>=<sigma/sigma_SIE>= =1.01+-0.02 with 0.065 rms scatter. We conclude that within the Einstein radii the SLACS lenses are very well approximated by isothermal ellipsoids, requiring a fine tuning of the stellar and dark matter distribution (bulge-halo ``conspiracy''). Interpreting the offset from the local FP in terms of evolution of the stellar mass-to-light ratio, we find for the SLACS lenses d log M/L_B/dz=-0.69+-0.08 (rms 0.11) consistent with the rate found for field early-type galaxies and with a scenario where most of the stars were formed at high redshift (>2) with secondary episodes of star formation providing less than ~10% of the stellar mass below z=1. We discuss star formation history and structural homogeneity in the context of formation mechanisms such as collisionless (``dry'') mergers. [Abridged]

We present Spitzer Space Telescope observations of the well-studied extremely red objects (EROs) HR10 and LBDS53W091 from 3.6 to 160 microns. These galaxies are the prototypes of the two primary classes of EROs: dusty starbursts and old, evolved galaxies, respectively. Both galaxies, as well as LBDS53W069, another example of an old, quiescent galaxy, are well-detected out to 8 microns. However, only the dusty starburst HR10 is detected in the far-infrared. All three EROs have stellar masses of a few times 10^11 M(sun). Using evolutionary model fits to their multiband photometry, we predict the infrared colors of similar EROs at 1<z<2. We find that blueward of observed 10 microns, the two ERO classes are virtually indistinguishable photometrically. Deep spectroscopy and 24 micron data allow the classes to be separated.

In this the first of a series of Letters, we present a description of the panchromatic data sets that have been acquired in the Extended Groth Strip region of the sky. Our survey, the All-wavelength Extended Groth strip International Survey (AEGIS), is intended to study the physical properties and evolutionary processes of galaxies at z ~ 1. It includes the following deep, wide-field imaging data sets: Chandra/ACIS X-ray (0.5 - 10 keV), GALEX ultraviolet (1200 - 2500 Angstrom), CFHT/MegaCam Legacy Survey optical (3600 - 9000 Angstroms), CFHT/CFH12K optical (4500 - 9000 Angstroms), Hubble Space Telescope/ACS optical (4400 - 8500 Angstroms), Palomar/WIRC near-infrared (1.2 - 2.2 microns), Spitzer/IRAC mid-infrared (3.6 - 8.0 microns), Spitzer/MIPS far-infrared (24 - 70 microns), and VLA radio continuum (6 - 20 cm). In addition, this region of the sky has been targeted for extensive spectroscopy using the DEIMOS spectrograph on the Keck II 10 m telescope. Our survey is compared to other large multiwavelength surveys in terms of depth and sky coverage.

We report on the results of a visual search for galaxy-scale strong gravitational lenses over 650 arcmin^2 of HST/ACS imaging in the DEEP2-EGS field. In addition to a previously-known Einstein Cross (the "Cross," HST J141735+52264, with z_lens=0.8106 and a published z_source=3.40), we identify two new strong galaxy-galaxy lenses with multiple extended arcs. The first, HST J141820+52361 (the ``Dewdrop''; z_lens=0.5798, lenses two distinct extended sources into two pairs of arcs z_source=0.while), 9818 the second, HST J141833+52435 (the ``Anchor''; z_lens=0.4625), produces a single pair of arcs (source redshift not yet known). All three definite lenses are fit well by simple singular isothermal ellipsoid models including external shear. Using the three-dimensional line-of-sight (LOS) information on galaxies from the DEEP2 data, we calculate the convergence and shear contributions, assuming singular isothermal sphere halos truncated at 200 h^-1 kpc. These are also compared against three-dimensional local-density estimates. We find that even strong lenses in demonstrably underdense local environments may be considerably affected by LOS contributions, which in turn, may be underestimates of the effect of large scale structure.

We report the discovery of a new Einstein cross at redshift z_S = 2.701 based on Lyman-alpha emission in a cruciform configuration around an SDSS luminous red galaxy (z_L = 0.331). The system was targeted as a possible lens based on an anomalous emission line in the SDSS spectrum. Imaging and spectroscopy from the W. M. Keck Observatory confirm the lensing nature of this system. This is one of the widest-separation galaxy-scale lenses known, with an Einstein radius of ~1.84 arcsec. We present simple gravitational lens models for the system and compute the intrinsic properties of the lensed galaxy. The total mass of the lensing galaxy within the 8.8 +/- 0.1 kpc enclosed by the lensed images is (5.2 +/- 0.1) x 10^11 M_sun. The lensed galaxy is a low mass galaxy (0.2 L*) with a high equivalent-width Lyman-alpha line (EW_Lya_rest = 46 +/- 5 Angstroms). Follow-up studies of this lens system can probe the mass structure of the lensing galaxy, and can provide a unique view of an intrinsically faint, high-redshift, star-forming galaxy at high signal-to-noise ratio.

We present the second campaign of the ESO/GOODS program of spectroscopy of faint galaxies in the GOODS-South field. Objects were selected as candidates for VLT/FORS2 observations primarily based on the expectation that the detection and measurement of their spectral features would benefit from the high throughput and spectral resolution of FORS2. The reliability of the redshift estimates is assessed using the redshift-magnitude and color-redshift diagrams and comparing the results with public data. 807 spectra of 652 individua targets have been obtained in service mode with the FORS2 spectrograph at the ESO/VLT, providing 501 redshift determinations. The typical redshift uncertainty is estimated to be sigma_z ~ 0.001. Galaxies have been selected adopting three different color criteria and using the photometric redshifts.The resulting redshift distribution typically spans two redshift domains: from z=0.5 to 2 and z=3.5 to 6.2. In particular, 94 B435-,V606-,i775 "dropout" Lyman break galaxies have been observed, yielding redshifts for 65 objects in the interval 3.4<z<6.2. Three sources have been serendipitously discovered in the redshift interval 4.8<z<5.5. Together with the previous release, 930 sources have now been observed and 724 redshift determinations have been carried out. The reduced spectra and the derived redshifts are released to the community through the ESO web page http://www.eso.org/science/goods/. Large scale structures are clearly detected at z=0.666, 0.734, 1.096, 1.221, 1.300, and 1.614. A sample of 34 sources with tilted [OII]3727 emission has been identified, 32 of them in the redshift range 0.9<z<1.5.

(Abridged) We present a joint gravitational lensing and stellar dynamical analysis of fifteen massive field early-type galaxies, selected from the Sloan Lens (SLACS) Survey. The following numerical results are found: (i) A joint-likelihood gives an average logarithmic density slope for the total mass density of 2.01 (+0.02/-0.03) (68 perecnt C.L). inside the Einstein radius. (ii) The average position-angle difference between the light distribution and the total mass distribution is found to be 0+-3 degrees, setting an upper limit of <= 0.035 on the average external shear. (iii) The average projected dark-matter mass fraction is inferred to be 0.25+-0.06 inside R_E, using the stellar mass-to-light ratios derived from the Fundamental Plane as priors. (iv) Combined with results from the LSD Survey, we find no significant evolution of the total density slope inside one effective radius: a linear fit gives d\gamma'/dz = 0.23+-0.16 (1-sigma) for the range z=0.08-1.01. The small scatter and absence of significant evolution in the inner density slopes suggest a collisional scenario where gas and dark matter strongly couple during galaxy formation, leading to a total mass distribution that rapidly converge to dynamical isothermality.

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

We investigate the properties of massive galaxies at z=1-3.5 using HST observations, ground-based near-IR imaging, and Spitzer Space Telescope observations at 3-24 micron. We identify 153 distant red galaxies (DRGs) with J-K > 2.3 mag (Vega) in the southern GOODS field. This sample is approximately complete in stellar mass for passively evolving galaxies above 10^11 solar masses and z < 3. The galaxies identified by this selection are roughly split between objects whose optical and near-IR rest-frame light is dominated by evolved stars combined with ongoing star formation, and galaxies whose light is dominated by heavily reddened starbursts. Very few of the galaxies (< 10%) have no indication of current star formation. Using SFR estimates that include the reradiated IR emission, the DRGs at z=1.5-3 with stellar masses > 10^11 solar masses have specific SFRs (SFRs per unit stellar mass) ranging from 0.2 to 10 Gyr^-1, with a mean value of ~2.4 Gyr^-1. The DRGs with stellar masses > 10^11 solar masses and 1.5 < z < 3 have integrated specific SFRs greater the global value over all galaxies. In contrast, we find that galaxies at z = 0.3-0.75 with these stellar masses have integrated specific SFRs less than the global value, and more than an order of magnitude lower than that for massive DRGs at z = 1.5-3. At z < 1, lower-mass galaxies dominate the overall cosmic mass assembly. This suggests that the bulk of star formation in massive galaxies occurs at early cosmic epochs and is largely complete by z~1.5. [Abridged]

We present results from a search for high-redshift J--band ``dropout'' galaxies in the portion of the GOODS southern field that is covered by extremely deep imaging from the Hubble Ultradeep Field (HUDF).Using observations at optical, near-infrared and mid-infrared wavelengths from the Hubble and Spitzer Space Telescopes and the ESO-VLT, we search for very massive galaxies at high redshifts and find one particularly remarkable candidate. Its spectral energy distribution is consistent with a galaxy at z ~ 6.5 and a stellar mass of 6x10e11 M(sun) (for a Salpeter IMF). We interpret a prominent photometric break between the near-infrared and Spitzer bandpasses as the 3646A Balmer discontinuity. The best-fitting models have low reddening and ages of several hundred Myr, placing the formation of the bulk of the stars at z > 9. Alternative models of dusty galaxies at z ~ 2.5 are possible but provide significantly poorer fits. The object is detected with Spitzer at 24 micron. This emission originats from an obscured active nucleus or star formation. We present optical and near-infrared spectroscopy which has, thus far, failed to detect any spectral features. This helps limit the solution in which the galaxy is a starburst or active galaxy at z ~ 2.5. If the high-redshift interpretation is correct, this object would be an example of a galaxy that formed by a process strongly resembling traditional models of monolithic collapse, in a way which a very large mass of stars formed within a remarkably short period of time, at very high redshift.

We present a two-pronged approach to the formation of early-type galaxies, using a sample of 18 galaxies at 0.5 < z < 1 from the HST/ACS Ultra Deep Field and GRAPES surveys: 1) We combine slitless low resolution spectroscopy from the GRAPES dataset with simple models of galaxy formation to explore their star formation histories. 2) We also perform an analysis of their surface brightness distribution with the unprecedented details provided by the ACS superb angular resolution and photometric depth. Our spectroscopic analysis reveals that their stellar populations are rather homogeneous in age and metallicity and formed at redshifts z ~ 2-5. Evolving them passively, they become practically indistinguishable from ellipticals at z = 0. Also, their isophotal shapes appear very similar to those observed for nearby ellipticals, in that the percentages of disky and boxy galaxies at z ~ 1 are close to the values measured at z = 0. Moreover,we find that the isophotal structure of z ~ 1 early-type galaxies obeys the correlations already observed among nearby ellipticals, i.e. disky ellipticals have generally higher characteristic ellipticities, and boxy ellipticals have larger half-light radii and are brighter in the restframe B band. In this respect then, no significant structural differences are seen for ellipticals between z = 0 and 1.

We use data from the first epoch of observations with the IRAC/Spitzer for the GOODS to detect and study a collection of LBGs at z ~ 6 to 5 in the HUDF, six of which have spectroscopic confirmation. At these redshifts, IRAC samples rest-frame optical light in the range 0.5 to 0.8 um, where the effects of dust extinction are smaller and the sensitivity to light from evolved stars is greater than at shorter, rest-frame UV observable from the ground or with the HST. As such, it provides useful constraints on the ages and masses of these galaxies' stellar populations. We find that the SEDs for many of these galaxies are best fitted by models of stellar populations with masses of a few x 1e10 M_sun, and with ages of a few hundred Myrs, values quite similar to those derived for typical LBGs at z ~ 3. When the universe was only 1 Gyr old, some galaxies had already formed a mass of stars approaching that of the present-day Milky Way, and that they started forming those stars at z > 7, and in some cases much earlier. We find that the lower limits to the space density for galaxies in this mass range are consistent with predictions from recent hydrodynamic simulations of structure formation in a LCDM universe. All objects in our samples are consistent with having solar metallicity, suggesting that they might have already been significantly polluted by metals. The values for dust reddening derived from the model fitting are low or zero, and we find that some of the galaxies have rest-frame UV colors that are even bluer than those predicted by the stellar population models to which we compare them. These colors might be attributed to the presence of very massive stars (> 100 M_sun), or by weaker intergalactic HI absorption than what is commonly assumed.

Using Hubble Space Telescope and ground-based U through K- band photometry from the Great Observatories Origins Deep Survey (GOODS), we measure the evolution of the luminosity function and luminosity density in the rest-frame optical (UBR) to z ~ 2, bridging the poorly explored ``redshift desert'' between z~1 and z~2. We also use deep near-infrared observations to measure the evolution in the rest-frame J-band to z~1. Compared to local measurements from the SDSS, we find a brightening of the characteristic magnitude, (M*), by ~2.1, \~0.8 and ~0.7 mag between z=0.1 and z=1.9, in U, B, and R bands, respectively. The evolution of M* in the J-band is in the opposite sense, showing a dimming between redshifts z=0.4 and z=0.9. This is consistent with a scenario in which the mean star formation rate in galaxies was higher in the past, while the mean stellar mass was lower, in qualitative agreement with hierarchical galaxy formation models. We find that the shape of the luminosity function is strongly dependent on spectral type and that there is strong evolution with redshift in the relative contribution from the different spectral types to the luminosity density. We find good agreement in the luminosity function derived from an R-selected and a K-selected sample at z~1, suggesting that optically selected surveys of similar depth (R < 24) are not missing a significant fraction of objects at this redshift relative to a near-infrared-selected sample. We compare the rest-frame B-band luminosity functions from z~0--2 with the predictions of a semi-analytic hierarchical model of galaxy formation, and find qualitatively good agreement. In particular, the model predicts at least as many optically luminous galaxies at z~1--2 as are implied by our observations.

We present measurements of the clustering properties of bright ($L>L_{*}$) z$\sim$4 Lyman Break Galaxies (LBGs) selected from the Oxford-Dartmouth Thirty Degree Survey (ODT). We describe techniques used to select and evaluate our candidates and calculate the angular correlation function which we find best fitted by a power law, $\omega(\theta)=A_{w}\theta^{-\beta}$ with $A_{w}=15.4$ (with $\theta$ in arcseconds), using a constrained slope of $\beta=0.8$. Using a redshift distribution consistent with photometric models, we deproject this correlation function and find a comoving $r_{0}=11.4_{-1.9}^{+1.7}$ h$_{100}^{-1}$ Mpc in a $\Omega_m=0.3$ flat $\Lambda$ cosmology for $i_{AB}\leq24.5$. This corresponds to a linear bias value of $b=8.1_{-2.6}^{+2.0}$ (assuming $\sigma_{8}=0.9$). These data show a significantly larger $r_{0}$ and $b$ than previous studies at $z\sim4$. We interpret this as evidence that the brightest LBGs have a larger bias than fainter ones, indicating a strong luminosity dependence for the measured bias of an LBG sample. Comparing this against recent results in the literature at fainter (sub-$L_{*}$) limiting magnitudes, and with simple models describing the relationship between LBGs and dark matter haloes, we discuss the implications on the implied environments and nature of LBGs. It seems that the brightest LBGs (in contrast with the majority sub-$L_{*}$ population), have clustering properties, and host dark matter halo masses, that are consistent with them being progenitors of the most massive galaxies today.

We report on a complete sample of 7 luminous early-type galaxies in the Hubble Ultra Deep Field (UDF) with spectroscopic redshifts between 1.39 and 2.47 and to K<23 AB. Using the BzK selection criterion we have pre-selected a set of objects over the UDF which fulfill the photometric conditions for being passively evolving galaxies at z>1.4. Low-resolution spectra of these objects have been extracted from the HST+ACS grism data taken over the UDF by the GRAPES project. Redshift for the 7 galaxies have been identified based on the UV feature at rest frame 2640<lambda<2850 AA. This feature is mainly due to a combination of FeII, MgI and MgII absorptions which are characteristic of stellar populations dominated by stars older than about 0.5 Gyr. The redshift identification and the passively evolving nature of these galaxies is further supported by the photometric redshifts and by the overall spectral energy distribution (SED), with the ultradeep HST+ACS/NICMOS imaging revealing compact morphologies typical of elliptical/early-type galaxies. From the SED we derive stellar masses of 10^{11}Msun or larger and ages of about 1 Gyr. Their space density at < z >=1.7 appears to be roughly a factor of 2--3 smaller than that of their local counterparts, further supporting the notion that such massive and old galaxies are already ubiquitous at early cosmic times. Much smaller effective radii are derived for some of the objects compared to local massive ellipticals, which may be due to morphological K corrections, evolution, or the presence of a central point-like source. Nuclear activity is indeed present in a subset of the galaxies, as revealed by them being hard X-ray sources, hinting to AGN activity having played a role in discontinuing star formation.

We present grism spectra taken with the Advanced Camera for Surveys to identify 29 red sources with (i'-z') > 0.9 in the Hubble Ultra Deep Field (HUDF). Of these 23 are found to be galaxies at redshifts between z=5.4 and 6.7, identified by the break at 1216 angstrom due to IGM absorption; two are late type dwarf stars with red colors; and four are galaxies with colors and spectral shape similar to dust reddened or old galaxies at redshifts z~1-2. This constitutes the largest uniform, flux-limited sample of spectroscopically confirmed galaxies at such faint fluxes z'<27.5. Many are also among the most distant spectroscopically confirmed galaxies (at redshifts up to z=6.7). We find a significant overdensity of galaxies at redshifts z=5.9 +/- 0.2. Nearly two thirds of the galaxies in our sample (15/23) belong to this peak. Taking into account the selection function and the redshift sensitivity of the survey, we get a conservative overdensity of at least a factor of two along the line-of-sight. The galaxies found in this redshift peak are also localized in the plane of the sky in a non-random manner, occupying about half of the ACS chip. Thus the volume overdensity is a factor of four. The star-formation rate derived from detected sources in this overdense region is sufficient to reionize the local IGM.

We have discovered an extended Lyman alpha plume associated with a compact source at redshift 5.4 in slitless spectroscopic data from the Grism ACS Program for Extragalactic Science (GRAPES) project. The spatial extent of the emission is about 6 x 1.5 kpc (1 x 0.25 arcsec). Combining our grism data and the broadband images from the Hubble UltraDeep Field (UDF) images, we find a Lyman alpha line flux of 2e-17 erg/cm2/s and surface brightness 7e-17 erg/cm2/s/arcsec2. The UDF images show diffuse continuum emission associated with the Lyman alpha plume (hereafter UDF 5225), with three embedded knots. The morphology of UDF 5225 is highly suggestive of a galaxy in assembly. It is moreover possible that the prominent Lyman alpha emission from this object is due to an active nucleus, and that we are seeing the simultaneous growth through accretion of a galaxy and its central black hole. Followup observations at higher spectral resolution could test this hypothesis.

We identified 46 unresolved source candidates in the Hubble Ultra Deep Field, down to i775 = 29.5. Unresolved objects were identified using a parameter S, which measures the deviation from the curve-of-growth of a point source. Extensive testing of this parameter was carried out, including the effects of decreasing signal-to-noise and of the apparent motions of stars, which demonstrated that stars brighter than i775 = 27.0 could be robustly identified. Low resolution grism spectra of the 28 objects brighter than i775 = 27.0 identify 18 M and later stellar type dwarfs, 2 candidate L-dwarfs, 2 QSOs, and 4 white dwarfs. Using the observed population of dwarfs with spectral type M4 or later, we derive a Galactic disk scale height of 400 \pm 100 pc for M and L stars. The local white dwarf density is computed to be as high as (1.1 \pm 0.3) x10^(-2) stars/pc^3. Based on observations taken 73 days apart, we determined that no object in the field has a proper motion larger than 0.027"/year (3 sigma detection limit). No high velocity white dwarfs were identified in the HUDF, and all four candidates appear more likely to be part of the Galactic thick disk. The lack of detected halo white dwarfs implies that, if the dark matter halo is 12 Gyr old, white dwarfs account for less than 10% of the dark matter halo mass.

(abridged) We present Gemini integral-field unit (IFU) spectroscopy and Hubble Space Telescope (HST) F435W- and F814W-band images of a new four-image gravitational lens, SDSSJ140228.22+632133.3, obtained as part of an HST Snapshot program designed to expand the sample of known gravitational lenses amenable to detailed photometric, lensing, and dynamical studies. The lens is a smooth elliptical galaxy at a redshift of z_l = 0.2046 +/- 0.0001 with a Sloan r-band magnitude of 17.00 +/- 0.05 and a stellar velocity dispersion of 267 +/- 17 km s^-1, obtained from its SDSS spectrum. Multiple emission lines place the quadruply-imaged source at a redshift of z_s = 0.4814 +/- 0.0001. The best-fitting singular isothermal ellipsoid lens model gives an Einstein radius b = 1.35" +/- 0.05" (or [4.9 +/- 0.2] h_65^-1 kpc), corresponding to a total mass of (30.9 +/- 2.3) x 10^10 h_65^-1 M_sun within the critical curve. In combination with HST photometry this gives a rest-frame B-band mass-to-light ratio of (8.1 +/- 0.7) h_65 times solar within the same region. The lens model predicts a luminosity-weighted stellar dispersion within the 3"-diameter SDSS aperture of sigma_* ~= 270 km s^-1, in good agreement with the observed value. Using the model to de-lens the four lensed images yields a source with a smooth, monotonically-decreasing brightness distribution.

We constrain the X-ray emission properties of Lyman break galaxies (LBGs) at z ~ 3-6 using the 2 Ms Chandra Deep Field-North and 1 Ms Chandra Deep Field-South. Large samples of LBGs were discovered using HST as part of the Great Observatories Origins Deep Survey (GOODS). Deep optical and X-ray imaging over the GOODS fields have allowed us to place the most significant constraints on the X-ray properties of LBGs to date. Mean X-ray properties of 449, 1734, 629, and 247 LBGs with z ~ 3, 4, 5, and 6, respectively, were determined using stacking techniques. When stacked, we detect X-ray emission from LBGs at z ~ 3 (7 sigma) and from an optically bright subset (brightest 25%) of LBGs at z ~ 4 (3 sigma); the latter is the highest redshift detection yet for "normal" galaxies in the X-ray band. The effective exposure times for these stacked observations are 0.7 and 0.5 Gs, respectively. The derived average rest-frame 2.0-8.0 keV luminosities are 1.5 and 1.4 x 10^41 erg s^-1, respectively. X-ray emission from these LBGs is likely due to high mass X-ray binaries (HMXBs) and Type II supernovae; the corresponding star formation rates are 85-240 M_sol yr^-1. The X-ray to B-band mean luminosity ratio (Lx/LB) at z~3 is somewhat elevated with respect to that measured for starburst galaxies in the local Universe. When stacking full samples of LBGs at z ~ 4, 5, and 6 we do not obtain significant detections and derive rest-frame 2.0-8.0 keV luminosity upper limits of 0.9, 2.8, and 7.1 x 10^41 erg s^-1, respectively. These upper limits constrain any widespread AGN activity in these objects to be modest at best. Furthermore, we find that ~ 0.5% of our LBGs from z ~ 3-6 are detected individually in the X-ray band. Continued ...

Using early data from the Infrared Array Camera (IRAC) on the Spitzer Space Telescope, taken for the Great Observatories Origins Deep Survey (GOODS), we identify and study objects that are well-detected at 3.6um, but are very faint (and in some cases, invisible) in the Hubble Ultra Deep Field (HUDF) ACS and NICMOS images and in very deep VLT Ks-band imaging. We select a sample of 17 objects with f_nu(3.6um)/f_nu(z_{850})>20. The analysis of their spectral energy-distributions (SEDs) from 0.4 to 8.0um shows that the majority of these objects cannot be satisfactorily explained without a well-evolved stellar population. We find that most of them can be well fitted by a simple two-component model, where the primary component represents a massive, old population that dominates the strong IR emission, while the secondary component represents a low amplitude, on-going star-formation process that accounts for the weak optical fluxes. Their estimated photometric redshifts (z_p) range from 1.6 to 2.9, with the median at z_p=2.4. For the simple star formation histories considered here, their corresponding stellar masses range from 0.1 to 1.6e11 M_sun for a Chabrier Initial Mass Function (IMF). Their median rest-frame Ks-band absolute magnitude is -22.9 mag in AB system, or 1.5xL*(K) for present--day elliptical galaxies. In the scenario of pure luminosity evolution, such objects may be direct progenitors for at least 14 to 51% of the local population of early--type galaxies. Due to the small cosmic volume of the HUDF, however, this simple estimate could be affected by other effects such as cosmic variance and the strong clustering of massive galaxies. A full analysis of the entire GOODS area is now underway to assess such effects.

The deep X-ray, optical, and far-infrared fields that constitute GOODS are sensitive to obscured AGN (N_H>10^{22} cm^{-2}) at the quasar epoch (z~2-3), as well as to unobscured AGN as distant as z~7. Luminous X-ray emission is a sign of accretion onto a supermassive black hole and thus reveals all but the most heavily obscured AGN. We combine X-ray luminosity functions with appropriate spectral energy distributions for AGN to model the X-ray, optical and far-infrared flux distributions of the X-ray sources in the GOODS fields. A simple model based on the unified paradigm for AGN, with ~3 times as many obscured AGN as unobscured, successfully reproduces the z-band flux distributions measured in the deep HST ACS observations on the GOODS North and South fields. This model is also consistent with the observed spectroscopic and photometric redshift distributions once selection effects are considered. The previously reported discrepancy between observed spectroscopic redshift distributions and the predictions of population synthesis models for the X-ray background can be explained by bias against the most heavily obscured AGN generated both by X-ray observations and the identification of sources via optical spectroscopy. We predict the AGN number counts for Spitzer MIPS 24 um and IRAC 3.6-8 um observations in the GOODS fields, which will verify whether most AGN in the early Universe are obscured in the optical. Such AGN should be very bright far-infrared sources and include some obscured AGN missed even by X-ray observations.

We present the first results of the ESO/GOODS program of spectroscopy of faint galaxies in the Chandra Deep Field South (CDF-S). 399 spectra of 303 unique targets have been obtained in service mode with the FORS2 spectrograph at the ESO/VLT, providing 234 redshift determinations (the median of the redshift distribution is at 1.04). The typical redshift uncertainty is estimated to be sig(z) ~ 0.001. Galaxies have been color selected in a way that the resulting redshift distribution typically spans from z=0.5 to 2. The reduced spectra and the derived redshifts are released to the community through the ESO web page http://www.eso.org/science/goods/ Large scale structure is clearly detected at z ~ 0.67, 0.73, 1.10 and 1.61. Three Lyman-break galaxies have also been included as targets and are confirmed to have redshifts z=4.800, 4.882 and 5.828. In a few cases, we observe clear [OII]3727 rotation curves, even at the relatively low resolution (R = 860) of the present observations. Assuming that the observed velocity structure is due to dynamically-relaxed rotation, this is an indication of large galactic masses (few times 10^(11) solar masses) at z ~ 1.

We present results from the Hubble Higher-z Supernova Search, the first space-based open field survey for supernovae (SNe). In cooperation with the Great Observatories Origins Deep Survey, we have used the Hubble Space Telescope with the Advanced Camera for Surveys to cover 300 square arcmin in the area of the Chandra Deep Field South and the Hubble Deep Field North on five separate search epochs (separated by 45 day intervals) to a limiting magnitude of z'=26. These deep observations have allowed us to discover 42 SNe in the redshift range 0.2 < z < 1.6. As these data span a large range in redshift, they are ideal for testing the validity of Type Ia supernova progenitor models with the distribution of expected ``delay times,'' from progenitor star formation to SN Ia explosion, and the SN rates these models predict. Through a Bayesian maximum likelihood test, we determine which delay-time models best reproduce the redshift distribution of SNe Ia discovered in this survey. We find that models that require a large fraction of ``prompt'' (less than 2 Gyr) SNe Ia poorly reproduce the observed redshift distribution and are rejected at 95% confidence. We find that Gaussian models best fit the observed data for mean delay times in the range of 3 to 4 Gyr.

The Spitzer Space Telescope has undertaken the deepest ever observations of the 24 micron sky in the ELAIS-N1 field as part of GOODS Science Verification observations. We present the shape of the 24 micron source counts in the flux range 20-1000 microJy, discuss the redshift distribution and nature of these sources with particular emphasis on their near-infrared counterparts.

We present deep unbiased spectroscopy of the Hubble Ultra Deep Field (UDF) carried out using the slitless grism spectroscopy mode of the Advance Camera for Surveys on the Hubble Space Telescope (HST). The GRIsm ACS Program for Extragalactic Science (GRAPES) achieves continuum detection as faint as $z_{AB}=27.2$ using 40 orbits ($9.2 \times 10^4$ seconds) on HST. The data were taken at four orientation angles to correct for the overlap of spectra. GRAPES data provide a unique, uninterrupted, low resolution (R=100) spectral coverage for $5500\AA < \lambda < 10500\AA$, and allow us to detect high redshift galaxies at $4 < z < 7$ whether they have $\lya$ lines or just show the Lyman Break, as well as find low luminosity AGNs in an unbiased fashion. This paper describes in detail the observations and the data reduction, and examines the quality of the extracted spectra. Subsequent papers will deal with the analysis of the data. The extracted and calibrated GRAPES spectra will be available from MAST at STScI.

Using an X-ray stacking procedure, we provide a robust upper limit to the X-ray luminosity per object of a set of 54 z~5.8 galaxy candidates in the Hubble Ultra Deep Field, which is within the 1 Ms-exposure Chandra Deep Field-South (CDF-S). With an effective total exposure of 44 Ms for the stack, the 3-sigma flux-density limit of 2.1E-17 erg/cm^2/s (soft-band) gives a 3-sigma upper-limit luminosity of L_X = 8E42 erg/s per object at a rest-frame hard energy range of 3-14 keV at z~5.8 for a photon index of Gamma=2. For an active accreting black hole (or "mini-quasar") emitting at the Eddington luminosity, and the Sazonov et al. average-QSO spectral energy distribution, we calculate an upper limit on the black hole mass, M_bh < 3E6 Msun (3-sigma). The X-ray limit further implies an upper limit on the rate density of UV ionizing photons from accreting black holes at that redshift, n-dot (ioniz) < 2E51 /s/Mpc^3 (3-sigma), which is less than 1/10 of the number needed to ionize the universe. Because the constraint is anchored in the rest-frame hard X-ray regime, a steeper spectrum for mini-quasars would imply relatively fewer ionizing UV photons. Unless there are large populations of active black holes around this mass that are unassociated with luminous galaxies, mini-quasars do not appear to contribute significantly to the budget of ionizing photons at z~6.

The Oxford Dartmouth Thirty Degree Survey (ODTS) is a deep, wide, multi-band imaging survey designed to cover a total of 30 square degrees in BVRi'Z, with a subset of U and K band data, in four separate fields of 5-10 deg^2 centred at 00:18:24 +34:52, 09:09:45 +40:50, 13:40:00 +02:30 and 16:39:30 +45:24. Observations have been made using the Wide Field Camera on the 2.5-m Isaac Newton Telescope in La Palma to average limiting depths (5 sigma Vega, aperture magnitudes) of U=24.8, B=25.6, V=25.0, R=24.6, and i'=23.5, with observations taken in ideal conditions reaching the target depths of U=25.3, B=26.2, V=25.7, R=25.4, and i'=24.6. The INT Z band data was found to be severely effected by fringing and, consequently, is now being obtained at the MDM observatory in Arizona. A complementary K-band survey has also been carried out at MDM, reaching an average depth of K_{5\sigma}~18.5. At present, approximately 23 deg^2 of the ODTS have been observed, with 3.5 deg^2 of the K band survey completed. This paper details the survey goals, field selection, observation strategy and data reduction procedure, focusing on the photometric calibration and catalogue construction. Preliminary photometric redshifts have been obtained for a subsample of the objects with R <= 23. These results are presented alongside a brief description of the photometric redshift determination technique used. The median redshift of the survey is estimated to be z~0.7 from a combination of the ODTS photometric redshifts and comparison with the redshift distributions of other surveys. Finally, galaxy number counts for the ODTS are presented which are found to be in excellent agreement with previous studies.

The integrated colors of distant galaxies provide a means for interpreting the properties of their stellar content. Here, we use rest-frame UV-to-optical colors to constrain the spectral-energy distributions and stellar populations of color-selected, B-dropout galaxies at z ~ 4 in the Great Observatories Origins Deep Survey. We combine the ACS data with ground-based near-infrared images, which extend the coverage of galaxies at z ~ 4 to the rest-frame B-band. We observe a color-magnitude trend in the rest-frame m(UV) - B versus B diagram for the z ~ 4 galaxies that has a fairly well-defined "blue-envelope", and is strikingly similar to that of color-selected, U-dropout galaxies at z ~ 3. We also find that although the co-moving luminosity density at rest-frame UV wavelengths (1600 Angstroms) is roughly comparable at z ~ 3 and z ~ 4, the luminosity density at rest-frame optical wavelengths increases by about one-third from z ~ 4 to z ~ 3. Although the star-formation histories of individual galaxies may involve complex and stochastic events, the evolution in the global luminosity density of the UV-bright galaxy population corresponds to an average star-formation history with a star-formation rate that is constant or increasing over these redshifts. This suggests that the evolution in the luminosity density corresponds to an increase in the stellar-mass density of more than 33%.

Spatially resolved spectroscopic data from the CIRPASS integral field unit (IFU) on Gemini are used to measure the gravitational lensing of the 4-image quasar Q2237+0305 on different size scales. A method for measuring the substructure present in the lens using observations at multiple wavelengths is demonstrated to be very effective and independent of many of the degeneracies inherent in previous methods. The magnification ratios of the QSO's narrow line region (NLR) and broad line region (BLR) are measured and found to be disagree with each other and with the published radio and mid-infrared magnification ratios. The disagreement between the BLR ratios and the radio/mid-infrared ratios is interpreted as microlensing by stars in the lens galaxy of the BLR The disagreement between the radio/mid-infrared ratios and the NLR ratios is interpreted as a signature of substructure on a larger scale, possibly the missing small scale structure predicted by the standard cold dark matter (CDM) model. Certain combinations of the radial profile and the substructure surface densities are ruled out using extensive lensing simulations. A substructure mass scale as large as 10^8 M is strongly disfavored while 10^4 M is too small if the radio and mid-infrared emission regions have the expected sizes of ~10 pc. The standard elliptical isothermal lens mass profile is not compatible with a substructure surface density of < 280 M/pc^2 at the 95% confidence level. This is 4-7% of the galaxy's surface density (depending on which image position is used to evaluate this). The required substructure surface density at the required mass scale is high in comparison with the present expectations within the CDM model.

We have measured the rest--frame lambda~1500 Ang comoving specific luminosity density of star--forming galaxies at redshift 3.5<z<6.5 from deep images taken with the Hubble Space Telescope (HST)and the Advanced Camera for Surveys (ACS), obtained as part of the Great Observatories Origins Deep Survey (GOODS). We used color selection criteria to construct samples of star--forming galaxies at redshifts z~4, 5 and 6, identified by the signature of the 912 Ang Lyman continuum discontinuity and Lyman-alpha forest blanketing in their rest--frame UV colors (Lyman--break galaxies). The ACS samples cover ~0.09 square degree, and are also relatively deep, reaching between 0.2 and 0.5 L_3^*, depending on the redshift, where $L_3^*$ is the characteristic UV luminosity of Lyman--break galaxies at z~3. The specific luminosity density of Lyman--break galaxies appears to be nearly constant with redshift from z~3 to z~6, although the measure at z~6 remains relatively uncertain, because it depends on the accurate estimate of the faint counts of the z~6 sample. If Lyman--break galaxies are fair tracers of the cosmic star formation activity, our results suggest that at z~6 the universe was already producing stars as vigorously as it did near its maximum several Gyr later, at 1<~z<~3. Thus, the onset of large--scale star formation in the universe is to be sought at around z~6 or higher, namely at less than ~7% of the current cosmic age.

Using U'- through Ks-band imaging data in the GOODS-South field, we construct a large, complete sample of 275 ``extremely red objects'' (EROs; K_s<22.0, R-K_s>3.35; AB), all with deep HST/ACS imaging in B_435, V_606, i_775, and z_850, and well-calibrated photometric redshifts. Quantitative concentration and asymmetry measurements fail to separate EROs into distinct morphological classes. We therefore visually classify the morphologies of all EROs into four broad types: ``Early'' (elliptical-like), ``Late'' (disk galaxies), ``Irregular'' and ``Other'' (chain galaxies and low surface brightness galaxies), and calculate their relative fractions and comoving space densities. For a broad range of limiting magnitudes and color thresholds, the relative number of early-type EROs is approximately constant at 33-44%, and the comoving space densities of Early- and Late-type EROs are comparable. Mean rest-frame spectral energy distributions (SEDs) at wavelengths between 0.1 and 1.2 um are constructed for all EROs. The SEDs are extremely similar in their range of shapes, independent of morphological type. The implication is that any differences between the broad-band SEDs of Early-type EROs and the other types are relatively subtle, and there is no robust way of photometrically distinguishing between different morphological types with usual optical/near-infrared photometry.

Cosmic variance is the uncertainty in observational estimates of the volume density of extragalactic objects such as galaxies or quasars arising from the underlying large-scale density fluctuations. This is often a significant source of uncertainty, especially in deep galaxy surveys, which tend to cover relatively small areas. We present estimates of the relative cosmic variance for one-point statistics (i.e. number densities) for typical scales and volumes sampled by the Great Observatories Origins Deep Survey (GOODS). We use two approaches: for objects with a known two-point correlation function that is well-approximated by a power law, one can use the standard analytic formalism to calculate the cosmic variance (in excess of shot noise). We use this approach to estimate the cosmic variance for several populations that are being studied in the GOODS program: Extremely Red Objects (ERO) at z~1, and Lyman Break Galaxies (LBG) at z~3 and z~4, using clustering information for similar populations in the literature. For populations with unknown clustering, one can use predictions from Cold Dark Matter theory to obtain a rough estimate of the variance as a function of number density. We present a convenient plot which allows one to use this approach to read off the cosmic variance for a population with a known mean redshift and estimated number density. We conclude that for the volumes sampled by GOODS, cosmic variance is a significant source of uncertainty for strongly clustered objects (~40-60% for EROs) and less serious for less clustered objects, ~10-20% for LBGs.

We report early results on galaxies at z~6, selected from Hubble Space Telescope imaging for the Great Observatories Origins Deep Survey. Spectroscopy of one object with the Advanced Camera for Surveys grism and from the Keck and VLT observatories a shows a strong continuum break and asymmetric line emission, identified as Lya at z=5.83. We detect only five spatially extended, z~6 candidates with signal-to-noise ratios > 10, two of which have spectroscopic confirmation. This is many fewer than would be expected if galaxies at z=6 had the same luminosity function as those at z=3. There are many fainter candidates, but we expect substantial contamination from foreground interlopers and spurious detections. Our best estimates favor a z=6 galaxy population with fainter luminosities, higher space density, and similar co-moving ultraviolet emissivity to that at z=3, but this depends critically on counts at fluxes fainter than those reliably probed by the current data.

Hubble Space Telescope images of high-redshift galaxies selected via color and photometric redshifts are used to examine the size and axial-ratio distribution of galaxies as a function of redshift at lookback times t > 8 Gyr. These parameters are measured at rest-frame UV wavelengths (1200-2000 Angstroms) on images with a rest-frame resolution of less than 0.8 kpc. Galaxy radii are found to scale with redshift approximately as the inverse of the Hubble parameter H(z). This is in accord with the theoretical expectation that the typical sizes of the luminous parts of galaxies should track the expected evolution in the virial radius of dark-matter halos. The mean ratio of semi-major to semi-minor axis for a bright well-resolved sample of galaxies at z ~ 4 is b/a = 0.65, suggesting that these Lyman break galaxies are not drawn from a spheroidal population. However the median concentration index of this sample is C = 3.5, which is closer to the typical concentration indices, C ~ 4, of nearby elliptical galaxies than to the values, C<2 for local disk galaxies of type Sb and later.

The redshift distribution of near-IR selected galaxies is often used to attempt to discriminate between the classical view of galaxy formation, in which present-day luminous galaxies were assembled at early times and evolve due to the passive aging of their stellar populations, and that of hierarchical structure formation, in which galaxies were assembled more recently via the merging of smaller objects. We carry out such a test here, by computing the distribution of photometric redshifts of K_AB<22 galaxies in the Great Observatories Origins Deep Field Survey GOODS Southern Field, and comparing the results with predictions from a semi-analytic model based on hierarchical structure formation, and a classical `passive evolution' model. We find that the redshift distributions at z <~ 1.5 of both the hierarchical and passive models are very similar to the observed one. At z>~1.5, the hierarchical model shows a deficit of galaxies, while the passive model predicts an excess. We investigate the nature of the observed galaxies in the redshift range where the models diverge, and find that the majority have highly disturbed morphologies, suggesting that they may be merger-induced starbursts. While the hierarchical model used here does not produce these objects in great enough numbers, the appearance of this population is clearly in better qualitative agreement with the hierarchical picture than with the classical passive evolution scenario. We conclude that the observations support the general framework of hierarchical formation, but suggest the need for new or modified physics in the models.

We use extensive multi-wavelength photometric data from the Great Observatories Origins Deep Survey (GOODS) to estimate photometric redshifts for a sample of 434 galaxies with spectroscopic redshifts in the Chandra Deep Field South. Using the Bayesian method, which incorporates redshift/magnitude priors, we estimate photometric redshifts for galaxies in the range 18 < R (AB) < 25.5, giving an rms scatter of 0.11. The outlier fraction is < 10%, with the outlier-clipped rms being 0.047. We examine the accuracy of photometric redshifts for several, special sub--classes of objects. The results for extremely red objects are more accurate than those for the sample as a whole, with rms of 0.051 and very few outliers (3%). Photometric redshifts for active galaxies, identified from their X-ray emission, have a dispersion of 0.104, with 10% outlier fraction, similar to that for normal galaxies. Employing a redshift/magnitude prior in this process seems to be crucial in improving the agreement between photometric and spectroscopic redshifts.

We present results from a systematic search for strong gravitational lenses in the GOODS ACS data. The search technique involves creating a sample of likely lensing galaxies, which we define as massive early-type galaxies in a redshift range 0.3 < z <1.3. The target galaxies are selected by color and magnitude, giving a sample of 1092 galaxies. For each galaxy in the sample, we subtract a smooth description of the galaxy light from the z_{850}-band data. The residuals are examined, along with true-color images created from the B_{435}V_{606}i_{775} data, for morphologies indicative of strong lensing. We present our six most promising lens candidates, as well as our full list of candidates.

We describe a possible new class of X-ray sources that have robust detections in ultra-deep Chandra data, yet have no detections at all in our deep multi-band GOODS Hubble Space Telescope (HST) ACS images, which represent the highest quality optical imaging obtained to date on these fields. These extreme X-ray / Optical ratio sources ("EXO"s) have values of Fx/Fopt at least an order of magnitude above those generally found for other AGN, even those that are harbored by reddened hosts. We thus infer two possible scenarios: (1) if these sources lie at redshifts z < 6, then their hosts need to be exceedingly underluminous, or more reddened, compared with other known sources; (2) if these sources lie above z ~ 6-7, such that even their Lyman-alpha emission is redshifted out of the bandpass of our ACS z(850) filter, then their optical and X-ray fluxes can be accounted for in terms of relatively normal L* hosts and moderate-luminosity AGN.

The Cold Dark Matter (CDM) model for galaxy formation predicts that a significant fraction of mass in the dark matter haloes that surround L L* galaxies is bound in substructures of mass 1E4-1E7Msun. The number of observable baryonic substructures (such as dwarf galaxies and globular clusters) falls short of these predictions by at least an order of magnitude. We present a method for searching for substructure in the haloes of gravitational lenses that produce multiple images of QSOs, such as 4-image Einstein Cross lenses. Current methods based on broadband flux ratios cannot cleanly distinguish between substructure, differential extinction, microlensing and, most importantly, ambiguities in the host lens model. These difficulties may be overcome by utilizing the prediction that when substructure is present, the magnification will be a function of source size. QSO broad line and narrow line emission regions are approximately ~1pc and >100pc in size, respectively. When narrow line region (NLR) features are used as a normalisation, the relative intensity and equivalent width of broad line region (BLR) features will respectively reflect substructure-lensing and microlensing effects. Spectroscopic observations of just a few image pairs would probably be able to cleanly extract the desired substructure signature and distinguish it from microlensing. In the rest-optical, the Hbeta/[OIII] region is ideal, since the narrow wavelength range also largely eliminates differential reddening problems. Simulations of Q2237+030 are done as an example to determine the level of substructure that is detectable in this way, and possible systematic difficulties are discussed. This is an ideal experiment to be carried out with near-infrared integral field unit spectrographs on 8-m class telescopes.

Wide field near-infrared (IR) surveys have revealed a population of galaxies with very red optical$-$IR colors, which have been termed ``Extremely Red Objects'' (EROs). Modeling suggests that such red colors (R-K > 5) could be produced by galaxies at z>~1 with either very old stellar populations or very high dust extinction. Recently it has been discovered that EROs are strongly clustered. Are these objects the high-redshift progenitors of present day giant ellipticals (gEs)? Are they already massive at this epoch? Are they the descendents of the $z\sim3$ Lyman Break Galaxies (LBG), which have also been identified as possible high redshift progenitors of giant ellipticals? We address these questions within the framework of the Cold Dark Matter paradigm using an analytic model that connects the number density and clustering or bias of an observed population with the halo occupation function (the number of observed galaxies per halo of a given mass). We find that EROs reside in massive dark matter halos, with average mass <M> > 1E13/h100 Msun. The occupation function that we derive for EROs is very similar to the one we derive for z=0, L>L* early type galaxies, whereas the occupation function for LBGs is skewed towards much smaller host halo masses (<M> ~ 1E11 to 1E12/h100 Msun. We then use the derived occupation function parameters to explore the possible evolutionary connections between these three populations.

The bulk of the molecular line emission in the double barred galaxy NGC4303 as observed in its CO(1-0) line with the OVRO mm-interferometer comes from two straight gas lanes which run north-south along the leading sides of the large-scale primary bar. Inside a radius of ~ 400 pc the molecular gas forms a spiral pattern which, for the northern arm, can be traced to the nucleus. Comparison of the OVRO and archival HST data with dynamical models of gas flow in the inner kiloparsec of single- and double-barred galaxies shows that the observed global properties of the molecular gas are in agreement with models for the gas flow in a strong, large-scale bar, and the two-arm spiral structure seen in CO in the inner kiloparsec can already be explained by a density wave initiated by the potential of that bar. Only a weak correlation between the molecular gas distribution and the extinction seen in the HST V-H map is found in the inner 400 pc of NGC4303: The innermost part of one arm of the nuclear CO spiral correlates with a weak dust filament in the color map, while the overall dust distribution follows a ring or single-arm spiral pattern well correlated with the UV continuum. This complicated nuclear geometry of the stellar and gaseous components allows for two scenarios: (A) A self-gravitating m=1 mode is present forming the spiral structure seen in the UV continuum. In this case the gas kinematics would be unaffected by the small (~ 4'') inner bar. (B) The UV continuum traces a complete ring which is heavily extincted north of the nucleus. Such a ring forms in hydrodynamic models of double bars, but the models cannot account for the UV emission observed on the leading side of the inner bar. (abridged)

In this article we investigate the morphology and stellar populations of high-redshift galaxies through multi-waveband HST imaging and ground-based spatially-resolved spectroscopy. We study the redshift evolution of galaxy morphology in the Hubble Deep Field, using the deep IDT-NICMOS near-infrared HST imaging coupled with spectroscopic and photometric redshifts. Using the multi-waveband data to compare the appearance of galaxies at the same rest-frame wavelengths reveals that "morphological k-corrections" (the change in appearance when viewing high-z objects at shorter rest-frame wavelengths) are only important in a minority of cases, and that galaxies were intrinsically more peculiar at high redshift. One example of significant morphological k-corrections is spiral galaxies, which often show more pronounced barred structure in the near-infrared than in the optical. Therefore, the apparent decline in the fraction of barred spirals at faint magnitudes in the optical HDF may be due to band-shifting effects at the higher redshifts, rather than intrinsic evolution. Using such features as the age-sensitive Balmer+4000Ang break, the spatially- resolved colours of distant galaxies in optical/near-infrared imaging can also be used to study their component stellar populations. We supplement this with deep Keck/LRIS spectroscopy of two extended sources: a chain galaxy at z=2.8 (HDF4-555.1, the "Hot Dog" - the brightest U-drop Lyman-break galaxy in the HDF) and a pair of z=4.04 gravitationally lensed arcs behind the cluster Abell 2390. The absence of measurable rotation across the z=2.8 chain galaxy implies that it is unikely to be a disk viewed edge on. With the resolution enhancement from lensing, we detect stellar populations of different ages in the z=4 arcs.

We present deep near-infrared Keck/NIRC imaging of a recently-discovered z=4.04 galaxy (Frye & Broadhurst 1998). This is lensed by the rich foreground cluster Abell~2390 (z~0.23) into highly-magnified arcs 3-5arcsec in length. Our H- and K'-band NIRC imaging allows us to map the Balmer+4000Ang break amplitude. In combination with high-quality archival HST/WFPC2 data, we can spatially resolve stellar populations along the arcs. The WFPC2 images clearly reveal several bright knots, which correspond to sites of active star formation. However, there are considerable portions of the arcs are significantly redder, consistent with being observed >100Myr after star formation has ceased. Keck/LRIS long-slit spectroscopy along the arcs reveals that the Ly-alpha emission is spatially offset by ~1arcsec from the rest-UV continuum regions. We show that this line emission is most probably powered by star formation in neighboring HII regions, and that the z=4 system is unlikely to be an AGN.

The Wavelength-Oriented Microwave Background Analysis Team (WOMBAT) is constructing microwave maps which will be more realistic than previous simulations. Our foreground models represent a considerable improvement: where spatial templates are available for a given foreground, we predict the flux and spectral index of that component at each place on the sky and estimate uncertainties. We will produce maps containing simulated CMB anisotropy combined with expected foregrounds. The simulated maps will be provided to the community as the WOMBAT Challenge, so such maps can be analyzed to extract cosmological parameters by scientists who are unaware of their input values. This will test the efficacy of foreground subtraction, power spectrum analysis, and parameter estimation techniques and help identify the areas most in need of progress. These maps are also part of the FORECAST project, which allows web-based access to the known foreground maps for the planning of CMB missions.

The Wavelength-Oriented Microwave Background Analysis Team (WOMBAT) is constructing microwave skymaps which will be more realistic than previous simulations. Our foreground models represent a considerable improvement: where spatial templates are available for a given foreground, we predict the flux and spectral index of that component at each place on the sky and estimate the uncertainties in these quantities. We will produce maps containing simulated Cosmic Microwave Background anisotropies combined with all major expected foreground components. The simulated maps will be provided to the cosmology community as the WOMBAT Challenge, a "hounds and hares" exercise where such maps can be analyzed to extract cosmological parameters by scientists who are unaware of their input values. This exercise will test the efficacy of current foreground subtraction, power spectrum analysis, and parameter estimation techniques and will help identify the areas most in need of progress.

We have recently obtained deep near-IR Keck imaging of a newly-discovered z=4.04 galaxy (Frye & Broadhurst 1998). This is lensed by the rich foreground cluster Abell 2390 (z=0.23) into highly-magnified arcs 3-5" in length. Our H- and K'-band Keck/NIRC imaging allows us to map the Balmer+4000Ang break amplitude. In combination with high-quality archival HST/WFPC2 data, we can spatially resolve stellar populations along the arcs. The WFPC2 images clearly reveal several bright knots, which may correspond to sites of active star formation. Indeed, in some spatial regions the Keck/LRIS discovery spectra are consistent with OB-star spectral energy distributions in the rest-ultraviolet. However, there are considerable portions of the arcs which appear redder with no Ly-alpha emission, consistent with being post-starburst regions.

Deep near-infrared and optical imaging surveys in the field reveal a curious population of galaxies that are infrared-bright (I-K>4), yet with relatively blue optical colors (V-I<2). Their surface density, several per square arcminute at K>20, is high enough that if placed at z>1 as our models suggest, their space densities are about one-tenth of phi-*. The colors of these ``faint red outlier galaxies'' (fROGs) may derive from exceedingly old underlying stellar populations, a dust-embedded starburst or AGN, or a combination thereof. Determining the nature of these fROGs, and their relation with the I-K>6 ``extremely red objects,'' has implications for our understanding of the processes that give rise to infrared-excess galaxies in general. We report on an ongoing study of several targets with HST & Keck imaging and Keck/LRIS multislit spectroscopy.

We present long-slit LRIS/Keck spectroscopic observations of the gravitational lens system 0957+561. Averaged over all of our data, the rest-frame velocity dispersion sigma_v of the central lens galaxy G1 is sigma_v = 279 +- 12 km/s. However, there appears to be a significant decrease in sigma_v as a function of distance from the center of G1 that is not typical of brightest cluster galaxies. Within 0.2" of the center of G1, we find the average sigma_v = 316+- 14 km/s, whereas for positions > 0.2" from the center of G1, we find the average sigma_v= 266+-12 km/s. A plausible explanation is that G1 contains a central massive dark object of mass M_MDO ~ 4x10^9 h_100^{-1} M_sun (h_100 = H_0/100km/s/Mpc), which contributes to the central velocity dispersion, and that the outer value of sigma_v is the appropriate measure of the depth of the potential well of G1. The determination of a luminosity-weighted estimate of sigma_v is essential for a determination of H_0 from 0957+561; our accurate measurements remove one of the chief uncertainties in the assumed form of the mass distribution of the lens. Thus, with the recent apparent reduction in the uncertainty in the measurement of the time delay for the images A and B of 0957+561, Delta tau_BA = 417+-3 d (Kundic et al. 1996), we obtain an estimate for the Hubble constant: H_0 = 62+-7 km/s/Mpc. If for some reason the trend of sigma_v with slit position is spurious and we should use the dispersion averaged along the slit, then the estimate of H_0 increases to 67+-8 km/s/Mpc. These standard errors do not, however, include any contribution from any errors in the assumed form of the mass distribution of the lens. In particular, we used the mass model described by Falco, Gorenstein & Shapiro (1991), as updated by Grogin & Narayan (1996a, b). The reduced chi^2

We present spectroscopy from the Keck telescope of three sets of objects in the Hubble Deep Field which have lens-like morphologies. In the case of J123641+621204, which is composed of four objects with similar colors and a mean separation of <= 0.8", we find at least two distinct components at redshifts of z=3.209 and z=3.220 which are separated by 0.5" spatially. Each of these components has narrow Ly-alpha emission, and possibly NV emission and SiIV and CIV in absorption or with a P-Cygni profile. The second case is J123652+621227, which has an arc-like feature offset by 1.8" to the southwest of a red elliptical-like galaxy, and a ``counterimage'' offset 1.4" on the opposite side. We tentatively find a single line at 5301 AA at the spatial position of the counterimage, and no corresponding emission line at the position of the arc. The colors of the counterimage are consistent with the identification of this line as Ly-alpha at z=3.36. The colors of the arc are different than those of the counterimage, and thus both the colors and spectra indicate that this object is unlikely to be a gravitational lens. For a third lensing candidate (J123656+621221), which is a blue arc offset by 0.9" from a red, elliptical-like galaxy, our spectroscopy does not clearly resolve the system spatially, complicating the interpretation of the spectrum. We discuss possible identifications of features in the spectrum, and find that gravitational lensing remains a possibility in this case. We conclude that the frequency of strong gravitational lensing by galaxies in the HDF appears to be very low. This result is difficult to reconcile with the introduction of a cosmological constant to account for the large number of faint blue galaxies via a large volume element at high redshift.

We combine deep K-band (Keck) with V- and I-band (NTT) observations of two high-Galactic latitude fields, surveying a total of ~2 sq. arcmin. The K-band galaxy counts continue to rise above K=22, reaching surface densities of few x 10^5 per sq. degree. The slope for the counts is (d log(N) per mag per sq. degree) = 0.23 +/- 0.02 between K=18-23, consistent with other deep K surveys. The numbers of galaxies in each mag bin is about two times greater than the galaxy counts of Djorgovski et al. (1995). The optical and near infrared magnitudes of all objects detected in the V+I+K image are discussed in the context of grids of isochrone synthesis galaxy evolutionary models (Bruzual & Charlot 1993, 1995). The colors of most of the observed galaxies are consistent with a population drawn from a broad redshift distribution. A few galaxies at K=19-20 are red in both colors (V-I>3; I-K>2), consistent with being early-type galaxies having undergone a burst of star formation at z>5 and viewed at z~1. At K>20, we find ~8 ``red outlier'' galaxies with I-K>4 and V-I<2.5, whose colors are difficult to mimic by a single evolving or non-evolving stellar population at any redshift. They are likely either low-metallicity, dusty dwarf galaxies, or old galaxies at high redshift (z>1.2). Their surface density is several per square arcminute, which is so high that they are probably common objects of low luminosity $L<L_*$.