The Habitable Exoplanet Observatory, or HabEx, has been designed to be the Great Observatory of the 2030s. For the first time in human history, technologies have matured sufficiently to enable an affordable space-based telescope mission capable of discovering and characterizing Earthlike planets orbiting nearby bright sunlike stars in order to search for signs of habitability and biosignatures. Such a mission can also be equipped with instrumentation that will enable broad and exciting general astrophysics and planetary science not possible from current or planned facilities. HabEx is a space telescope with unique imaging and multi-object spectroscopic capabilities at wavelengths ranging from ultraviolet (UV) to near-IR. These capabilities allow for a broad suite of compelling science that cuts across the entire NASA astrophysics portfolio. HabEx has three primary science goals: (1) Seek out nearby worlds and explore their habitability; (2) Map out nearby planetary systems and understand the diversity of the worlds they contain; (3) Enable new explorations of astrophysical systems from our own solar system to external galaxies by extending our reach in the UV through near-IR. This Great Observatory science will be selected through a competed GO program, and will account for about 50% of the HabEx primary mission. The preferred HabEx architecture is a 4m, monolithic, off-axis telescope that is diffraction-limited at 0.4 microns and is in an L2 orbit. HabEx employs two starlight suppression systems: a coronagraph and a starshade, each with their own dedicated instrument.

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

We present new ALMA observations of the molecular gas and far-infrared continuum around the brightest cluster galaxy (BCG) in the cool-core cluster MACS 1931.8-2635. Our observations reveal $1.9 \pm 0.3 \times 10^{10}$ M$_{\odot}$ of molecular gas, on par with the largest known reservoirs of cold gas in a cluster core. We detect CO(1-0), CO(3-2), and CO(4-3) emission from both diffuse and compact molecular gas components that extend from the BCG center out to $\sim30$ kpc to the northwest, tracing the UV knots and H$\alpha$ filaments observed by HST. Due to the lack of morphological symmetry, we hypothesize that the $\sim300$ km s$^{-1}$ velocity of the CO in the tail is not due to concurrent uplift by AGN jets, rather we may be observing the aftermath of a recent AGN outburst. The CO spectral line energy distribution suggests that molecular gas excitation is influenced by processes related to both star formation and recent AGN feedback. Continuum emission in Bands 6 and 7 arises from dust and is spatially coincident with young stars and nebular emission observed in the UV and optical. We constrain the temperature of several dust clumps to be $\lesssim 10$ K, which is too cold to be directly interacting with the surrounding $\sim 4.8$ keV intracluster medium (ICM). The cold dust population extends beyond the observed CO emission and must either be protected from interacting with the ICM or be surrounded by local volumes of ICM that are several keV colder than observed by Chandra.

Quantifying the physical conditions that allow radiation emitted shortward of the hydrogen ionization edge at 911.7 {\AA} to escape the first collapsed objects and ultimately reionize the universe is a compelling problem for astrophysics. The escape of LyC emission from star-forming galaxies and AGN is intimately tied to the emergence and sustenance of the metagalactic ionizing background that pervades the universe to the present day and in turn is tied to the emergence of structure at all epochs. JWST was built in part to search for the source(s) responsible for reionization, but it cannot observe LyC escape directly, because of the progressive increase in the mean transmission of the intergalactic medium towards the epoch of reionization. Remarkable progress has been made to date in directly detecting LyC leaking from star-forming galaxies using space-based and the ground-based observatories, but there remain significant gaps in our redshift coverage of the phenomenon. Ongoing projects to measure LyC escape at low- and intermediate-z will provide guidance to JWST investigations by analyzing the robustness of a set of proposed LyC escape proxies, and also provide a closeup examination of the physical conditions that favor LyC escape. However, currently available facilities are inadequate for deeply probing LyC escape at the faint end of the galaxy luminosity function. Doing so will require facilities that can detect LyC emission in the restframe to limiting magnitudes approaching 28 $< m^*_{(1+z)900} <$ 32 for $M^*_{(1+z)1500}$ galaxies. The goal of acquiring statistically robust samples for determining LyC luminosity functions across cosmic time will require multi-object spectroscopy from spacebased flagship class and groundbased ELT class telescopes along with ancillary panchromatic imaging and spectroscopy spanning the far-UV to the mid-IR.

The past decade has seen an explosion of discoveries and new insights into the diffuse gas within galaxies, galaxy clusters, and the filaments composing the Cosmic Web. A new decade will bring fresh opportunities to further this progress towards developing a comprehensive view of the composition, thermal state, and physical processes of diffuse gas in the Universe. Ultraviolet (UV) spectroscopy, probing diffuse 10^4-10^6 K gas at high spectral resolution, is uniquely poised to (1) witness environmental galaxy quenching processes in action, such as strangulation and tidal- and ram-pressure stripping, (2) directly account for the baryon content of galaxy clusters in the cold-warm (T<10^6 K) gas, (3) determine the phase structure and kinematics of gas participating in the equilibrium-regulating exchange of energy at the cores of galaxy clusters, and (4) map cold streams and filaments of the Cosmic Web that feed galaxies and clusters. With a substantial UV undertaking beyond the Hubble Space Telescope, all of the above would be achievable over the entire epoch of galaxy cluster formation. Such capabilities, coupled with already-planned advancements at other wavelengths, will transform extragalactic astronomy by revealing the dominant formation and growth mechanisms of gaseous halos over the mass spectrum, settling the debate between early- and late-time metal enrichment scenarios, and revealing how the ecosystems in which galaxies reside ultimately facilitate their demise.

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution.

Galaxy evolution is driven by many complex interrelated processes as galaxies accrete gas, form new stars, grow their stellar masses and central black holes, and subsequently quench. The processes that drive these transformations is poorly understood, but it is clear that the local environment on multiple scales plays a significant role. Today's massive clusters are dominated by spheroidal galaxies with low levels of star formation while those in the field are mostly still actively forming their stars. In order to understand the physical processes that drive both the mass build up in galaxies and the quenching of star formation, we need to investigate galaxies and their surrounding gas within and around the precursors of today's massive galaxy clusters -- protoclusters at z>2. The transition period before protoclusters began to quench and become the massive clusters we observe today is a crucial time to investigate their properties and the mechanisms driving their evolution. However, until now, progress characterizing the galaxies within protoclusters has been slow, due the difficulty of obtaining highly complete spectroscopic observations of faint galaxies at z>2 over large areas of the sky. The next decade will see a transformational shift in our understanding of protoclusters as deep spectroscopy over wide fields of view will be possible in conjunction with high resolution deep imaging in the optical and near-infrared.

The Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.

This paper describes the organization of the database and the catalog data products from the Pan-STARRS1 $3\pi$ Steradian Survey. The catalog data products are available in the form of an SQL-based relational database from MAST, the Mikulski Archive for Space Telescopes at STScI. The database is described in detail, including the construction of the database, the provenance of the data, the schema, and how the database tables are related. Examples of queries for a range of science goals are included. The catalog data products are available in the form of an SQL-based relational database from MAST, the Mikulski Archive for Space Telescopes at STScI.

Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ \Delta ra, \Delta dec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.

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.

According to the standard cold dark matter (CDM) cosmology, the structure of dark halos including those of galaxy clusters reflects their mass accretion history. Older clusters tend to be more concentrated than younger clusters. Their structure, represented by the characteristic radius $r_s$ and mass $M_s$ of the Navarro--Frenk--White (NFW) density profile, is related to their formation time. In~this study, we showed that $r_s$, $M_s$, and the X-ray temperature of the intracluster medium (ICM), $T_X$, form a thin plane in the space of $(\log r_s, \log M_s, \log T_X)$. This tight correlation indicates that the ICM temperature is also determined by the formation time of individual clusters. Numerical simulations showed that clusters move along the fundamental plane as they evolve. The plane and the cluster evolution within the plane could be explained by a similarity solution of structure formation of the universe. The angle of the plane shows that clusters have not achieved "virial equilibrium" in the sense that mass/size growth and pressure at the boundaries cannot be ignored. The distribution of clusters on the plane was related to the intrinsic scatter in the halo concentration--mass relation, which originated from the variety of cluster ages. The well-known mass--temperature relation of clusters ($M_\Delta\propto T_X^{3/2}$) can be explained by the fundamental plane and the mass dependence of the halo concentration without the assumption of virial equilibrium. The fundamental plane could also be used for calibration of cluster masses.

We reconstruct the two-dimensional (2D) matter distributions in 20 high-mass galaxy clusters selected from the CLASH survey by using the new approach of performing a joint weak lensing analysis of 2D shear and azimuthally averaged magnification measurements. This combination allows for a complete analysis of the field, effectively breaking the mass-sheet degeneracy. In a Bayesian framework, we simultaneously constrain the mass profile and morphology of each individual cluster assuming an elliptical Navarro-Frenk-White halo characterized by the mass, concentration, projected axis ratio, and position angle of the projected major axis.. We find that spherical mass estimates of the clusters from azimuthally averaged weak-lensing measurements in previous work are in excellent agreement with our results from a full 2D analysis. Combining all 20 clusters in our sample, we detect the elliptical shape of weak-lensing halos at the $5\sigma$ significance level within a scale of 2Mpc$/h$. The median projected axis ratio is $0.67\pm 0.07$ at a virial mass of $M_\mathrm{vir}=(15.2\pm 2.8) \times 10^{14} M_\odot$, which is in agreement with theoretical predictions of the standard collisionless cold dark matter model. We also study misalignment statistics of the brightest cluster galaxy, X-ray, thermal Sunyaev-Zel'dovich effect, and strong-lensing morphologies with respect to the weak-lensing signal. Among the three baryonic tracers studied here, we find that the X-ray morphology is best aligned with the weak-lensing mass distribution, with a median misalignment angle of $21\pm 7$ degrees. We also conduct a stacked quadrupole shear analysis assuming that the X-ray major axis is aligned with that of the projected mass distribution. This yields a consistent axis ratio of $0.67\pm 0.10$, suggesting again a tight alignment between the intracluster gas and dark matter.

We investigate the evolution of star formation rates (SFRs), stellar masses, and M/L$_{3.4 \mu m}$ ratios of brightest cluster galaxies (BCGs) in the COSMOS survey since z ~ 1 to determine the contribution of star formation to the growth-rate of BCG stellar mass over time. Through the spectral energy distribution (SED) fitting of the GALEX, CFHT, Subaru, Vista, Spitzer, and Herschel photometric data available in the COSMOS2015 catalog, we estimate the stellar mass and SFR of each BCG. We use a modified version of the iSEDfit package to fit the SEDs of our sample with both stellar and dust emission models, as well as constrain the impact of star formation history assumptions on our results. We find that in our sample of COSMOS BCGs, star formation evolves similarly to that in BCGs in samples of more massive galaxy clusters. However, compared to the latter, the magnitude of star formation in our sample is lower by ~ 1 dex. Additionally, we find an evolution of BCG baryonic mass-to-light ratio (M/L$_{3.4 \mu m}$) with redshift which is consistent with a passively aging stellar population. We use this to build upon Wen et al.'s (2013) low-redshift $\nu$L$_{3.4 \mu m}$-M$_{Stellar}$ relation, quantifying a correlation between $\nu$L$_{3.4 \mu m}$ and M$_{Stellar}$ to z ~ 1. By comparing our results to BCGs in Sunyaev-Zel'dovich and X-ray-selected samples of galaxy clusters, we find evidence that the normalization of star formation evolution in a cluster sample is driven by the mass range of the sample and may be biased upwards by cool cores.

In cold dark matter (CDM) cosmology, objects in the Universe have grown under the effect of gravity of dark matter. The intracluster gas in a galaxy cluster was heated when the dark-matter halo formed through gravitational collapse. The potential energy of the gas was converted to thermal energy through this process. However, this process and the thermodynamic history of the gas have not been clearly characterized in connection with with the formation and evolution of the internal structure of dark-matter halos. Here, we show that observational CLASH data of high-mass galaxy clusters lie on a plane in the three-dimensional logarithmic space of their characteristic radius $r_s$, mass $M_s$, and X-ray temperature $T_X$ with a very small orthogonal scatter. The tight correlation indicates that the gas temperature was determined at a specific cluster formation time, which is encoded in $r_s$ and $M_s$. The plane is tilted with respect to $T_X \propto M_s/r_s$, which is the plane expected in case of simplified virial equilibrium. We show that this tilt can be explained by a similarity solution, which indicates that clusters are not isolated but continuously growing through matter accretion from their outer environments. Numerical simulations reproduce the observed plane and its angle. This result holds independently of the gas physics implemented in the code, revealing the fundamental origin of this plane.

We present a new method for photometering objects in galaxy clusters. We introduce a mode-filtering technique for removing spatially variable backgrounds, improving both detection and photometric accuracy (roughly halving the scatter in the red sequence compared to previous catalogs of the same clusters). This method is based on robustly determining the distribution of background pixel values and should provide comparable improvement in photometric analysis of any crowded fields. We produce new multiwavelength catalogs for the 25 CLASH cluster fields in all 16 bandpasses from the UV through the near IR, as well as rest-frame magnitudes. A comparison with spectroscopic values from the literature finds a ~30% decrease in the redshift deviation from previously-released CLASH photometry. This improvement in redshift precision, in combination with a detection scheme designed to maximize purity, yields a substantial upgrade in cluster member identification over the previous CLASH galaxy catalog. We construct luminosity functions for each cluster, reliably reaching depths of at least 4.5 mag below M* in every case, and deeper still in several clusters. We measure M* , $\alpha$, and their redshift evolution, assuming the cluster populations are coeval, and find little to no evolution of $\alpha$, $-0.9\lesssim\langle\alpha\rangle\lesssim -0.8$, and M* values consistent with passive evolution. We present a catalog of galaxy photometry, photometric and spectroscopic redshifts, and rest-frame photometry for the full fields of view of all 25 CLASH clusters. Not only will our new photometric catalogs enable new studies of the properties of CLASH clusters, but mode-filtering techniques, such as those presented here, should greatly enhance the data quality of future photometric surveys of crowded fields.

We present a new strong lensing analysis of the galaxy cluster MACS J1206.2-0847 (MACS 1206), at z=0.44, using deep spectroscopy from CLASH-VLT and VLT/MUSE archival data in combination with imaging from the Cluster Lensing and Supernova survey with Hubble. MUSE observations enable the spectroscopic identification of 23 new multiply imaged sources, extending the previous compilations by a factor of approximately five. In total, we use the positional measurements of 82 spectroscopic multiple images belonging to 27 families at z=1.0-6.1 to reconstruct the projected total mass distribution of MACS 1206. Remarkably, 11 multiple images are found within 50 kpc of the brightest cluster galaxy, making this an unprecedented set of constraints for the innermost projected mass distribution of a galaxy cluster. We thus find that, although dynamically relaxed, the smooth matter component (dark matter plus hot gas) of MACS 1206 shows a significant asymmetry, which closely follows the asymmetric distribution of the stellar component (galaxy members and intracluster light). We determine the value of the innermost logarithmic slope of the projected total mass density profile and find it to be close to the canonical Navarro-Frenk-White value. We demonstrate that this quantity is very robust against different parametrizations of the diffuse mass component; however, this is not the case when only one central image is used in the mass reconstruction. We also show that the mass density profile from our new strong lensing model is in very good agreement with dynamical and X-ray measurements at larger radii, where they overlap.

We present an analysis of 15 Type Ia supernovae (SNe Ia) at redshift z > 1 (9 at 1.5 < z < 2.3) recently discovered in the CANDELS and CLASH Multi-Cycle Treasury programs using WFC3 on the Hubble Space Telescope. We combine these SNe Ia with a new compilation of 1050 SNe Ia, jointly calibrated and corrected for simulated survey biases to produce accurate distance measurements. We present unbiased constraints on the expansion rate at six redshifts in the range 0.07 < z < 1.5 based only on this combined SN Ia sample. The added leverage of our new sample at z > 1.5 leads to a factor of ~3 improvement in the determination of the expansion rate at z = 1.5, reducing its uncertainty to ~20%, a measurement of H(z=1.5)/H0=2.67 (+0.83,-0.52). We then demonstrate that these six measurements alone provide a nearly identical characterization of dark energy as the full SN sample, making them an efficient compression of the SN Ia data. The new sample of SNe Ia at z > 1 usefully distinguishes between alternative cosmological models and unmodeled evolution of the SN Ia distance indicators, placing empirical limits on the latter. Finally, employing a realistic simulation of a potential WFIRST SN survey observing strategy, we forecast optimistic future constraints on the expansion rate from SNe Ia.

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.

We study the nature of feedback mechanisms in the 11 CLASH brightest cluster galaxies (BCGs) that exhibit extended ultraviolet and nebular line emission features. We estimate star formation rates (SFRs), dust masses, and starburst durations using a Bayesian photometry-fitting technique that accounts for both stellar and dust emission from the UV through far IR. By comparing these quantities to intracluster medium (ICM) cooling times and free-fall times derived from X-ray observations and lensing estimates of the cluster mass distribution, we discover a tight relationship between the BCG SFR and the ICM cooling time to free-fall time ratio, $t_{cool}/t_{ff}$, with an upper limit on the intrinsic scatter of 0.15 dex. Furthermore, starburst durations may correlate with ICM cooling times at a radius of $0.025R_{500}$, and the two quantities converge upon reaching the Gyr regime. Our results provide a direct observational link between the thermodynamical state of the ICM and the intensity and duration of BCG star formation activity, and appear consistent with a scenario where active galactic nuclei (AGN) induce condensation of thermally unstable ICM overdensities that fuel long-duration ($>$ 1 Gyr) BCG starbursts. This scenario can explain (a) how gas with a low cooling time is depleted without causing a cooling flow and (b) the scaling relationship between SFR and $t_{cool}/t_{ff}$. We also find that the scaling relation between SFR and dust mass in BCGs with SFRs $<100$ M$_{\odot}$ yr$^{-1}$ is similar to star-forming field galaxies; BCGs with large ($>100$ M$_{\odot}$ yr$^{-1}$) SFRs have dust masses comparable to extreme starbursts.

We present accurate photometric redshifts for galaxies observed by the Cluster Lensing and Supernova survey with Hubble (CLASH). CLASH observed 25 massive galaxy cluster cores with the Hubble Space Telescope in 16 filters spanning 0.2 - 1.7 $\mu$m. Photometry in such crowded fields is challenging. Compared to our previously released catalogs, we make several improvements to the photometry, including smaller apertures, ICL subtraction, PSF matching, and empirically measured uncertainties. We further improve the Bayesian Photometric Redshift (BPZ) estimates by adding a redder elliptical template and by inflating the photometric uncertainties of the brightest galaxies. The resulting photometric redshift accuracies are dz/(1+z) $\sim$ 0.8\%, 1.0\%, and 2.0\% for galaxies with I-band F814W AB magnitudes $<$ 18, 20, and 23, respectively. These results are consistent with our expectations. They improve on our previously reported accuracies by a factor of 4 at the bright end and a factor of 2 at the faint end. Our new catalog includes 1257 spectroscopic redshifts, including 382 confirmed cluster members. We also provide stellar mass estimates. Finally, we include lensing magnification estimates of background galaxies based on our public lens models. Our new catalog of all 25 CLASH clusters is available via MAST. The analysis techniques developed here will be useful in other surveys of crowded fields, including the Frontier Fields and surveys carried out with J-PAS and JWST.

We determine the mass scale of Planck galaxy clusters using gravitational lensing mass measurements from the Cluster Lensing And Supernova survey with Hubble (CLASH). We have compared the lensing masses to the Planck Sunyaev-Zeldovich (SZ) mass proxy for 21 clusters in common, employing a Bayesian analysis to simultaneously fit an idealized CLASH selection function and the distribution between the measured observables and true cluster mass. We used a tiered analysis strategy to explicitly demonstrate the importance of priors on weak lensing mass accuracy. In the case of an assumed constant bias, $b_{SZ}$, between true cluster mass, $M_{500}$, and the Planck mass proxy, $M_{PL}$, our analysis constrains $1- b_{SZ} = 0.73\pm 0.10$ when moderate priors on weak lensing accuracy are used, including a zero-mean Gaussian with standard deviation of 8% to account for possible bias in lensing mass estimations. Our analysis explicitly accounts for possible selection bias effects in this calibration sourced by the CLASH selection function. Our constraint on the cluster mass scale is consistent with recent results from the Weighing the Giants program and the Canadian Cluster Comparison Project. It is also consistent, at 1.34$\sigma$, with the value needed to reconcile the Planck SZ cluster counts with Planck's base $\Lambda$CDM model fit to the primary cosmic microwave background anisotropies.

[abridged] The Hubble and Spitzer Space Telescope surveys of the Frontier Fields (FF) provide extremely deep images around six massive, strong-lensing clusters of galaxies. The ALMA FF survey aims to cover the same fields at 1.1mm, with maps reaching (unlensed) sensitivities of $<$70$\mu$Jy, in order to explore the properties of background dusty star-forming galaxies. We report on the multi-wavelength photometric analysis of all 12 significantly detected ($>$5$\sigma$) sources in the first three FF clusters observed by ALMA, based on data from Hubble and Spitzer, the VLT and Herschel. We measure the total photometry in all available bands and determine the photometric redshifts and the physical properties of the counterparts via SED-fitting. In particular, we carefully estimate the FIR photometry using 1.1mm priors to limit the misidentification of blended FIR counterparts, which strongly affect some flux estimates in previous FIR catalogs. We identify robust near-infrared (NIR) counterparts for all 11 sources with K$_s$ detection, the majority of which are quite red, with eight having $F814W-K_s\gtrsim 4$ and five having $F160W-[4.5]\gtrsim3$. From the FIR point of view, all our objects have $z_{phot}$$\sim$1--3, whereas based on the optical SED one object prefers a high-$z$ solution ($z\geq\ $7). Five objects among our sample have spectroscopic redshifts from the GLASS survey for which we can reproduce their SEDs with existing templates. This verification confirms the validity of our photometric redshift methodology. The mean redshift of our sample is $z_{phot}$=1.99$\pm$0.27. All 1.1mm selected objects are massive (10.0$<\log[M_{\star}(M_{\odot})]<$ 11.5), with high star formation rates ($<\log[SFR(M_{\odot}/yr)]> \approx$1.6) and high dust contents (8.1 $<\log[M_{dust} (M_{\odot})]<$8.8), consistent with previous ALMA surveys.

We present a detailed strong lensing (SL) mass reconstruction of the core of the galaxy cluster MACSJ 2129.4-0741 ($\rm z_{cl}=0.589$) obtained by combining high-resolution HST photometry from the CLASH survey with new spectroscopic observations from the CLASH-VLT survey. A background bright red passive galaxy at $\rm z_{sp}=1.36$, sextuply lensed in the cluster core, has four radial lensed images located over the three central cluster members. Further 19 background lensed galaxies are spectroscopically confirmed by our VLT survey, including 3 additional multiple systems. A total of 31 multiple images are used in the lensing analysis. This allows us to trace with high precision the total mass profile of the cluster in its very inner region ($\rm R<100$ kpc). Our final lensing mass model reproduces the multiple images systems identified in the cluster core with high accuracy of $0.4"$. This translates to an high precision mass reconstruction of MACS 2129, which is constrained at level of 2%. The cluster has Einstein parameter $\Theta_E=(29\pm4)"$, and a projected total mass of $\rm M_{tot}(<\Theta_E)=(1.35\pm0.03)\times 10^{14}M_{\odot}$ within such radius. Together with the cluster mass profile, we provide here also the complete spectroscopic dataset for the cluster members and lensed images measured with VLT/VIMOS within the CLASH-VLT survey.

We investigate the mass content of galaxies in the core of the galaxy cluster Abell 611. We perform a strong lensing analysis of the cluster core and use velocity dispersion measurements for individual cluster members as additional constraints. Despite the small number of multiply-imaged systems and cluster members with central velocity dispersions available in the core of A611, the addition of velocity dispersion measurements leads to tighter constraints on the mass associated with the galaxy component, and as a result, on the mass associated with the dark matter halo. Without the spectroscopic velocity dispersions, we would overestimate the mass of the galaxy component by a factor of $\sim1.5$, or, equivalently, we would underestimate the mass of the cluster dark halo by $\sim5\%$. We perform an additional lensing analysis using surface brightness (SB) reconstruction of the tangential giant arc. This approach improves the constraints on the mass parameters of the 5 galaxies close to the arc by up to a factor $\sim10$. The galaxy velocity dispersions resulting from the SB analysis are consistent at the 1$\sigma$ confidence level with the spectroscopic measurements and with the prediction from the simple pointlike analysis. In contrast the truncation radii for 2-3 galaxies depart significantly from the galaxy scaling relation and suggest differences in the stripping history from galaxy to galaxy.

We report new HST COS and STIS spectroscopy of a star-forming region (~100 solar masses/year) in the center of the X-ray cluster RXJ1532.9+3021 (z=0.362), to follow-up the CLASH team discovery of luminous UV filaments and knots in the central massive galaxy. We detect broad (~500 km/s) Lyman alpha emission lines with extraordinarily high equivalent width (EQW~200 Angstroms) and somewhat less broadened H-alpha (~220 km/s). Emission lines of N V and O VI are not detected, which constrains the rate at which gas cools through temperatures of 10^6 K to be less than about 10 solar masses/year. The COS spectra also show a flat rest-frame UV continuum with weak stellar photospheric features, consistent with the presence of recently-formed hot stars forming at a rate of ~10 solar masses/year, uncorrected for dust extinction. The slope and absorption lines in these UV spectra are similar to those of Lyman Break Galaxies at z approximately 3, albeit those with the highest Lyman-alpha equivalent widths and star-formation rates. This high-EQW Lyman-alpha source is a high-metallicity galaxy rapidly forming stars in structures that look nothing like disks. This mode of star formation could significantly contribute to the spheroidal population of galaxies. The constraint on the luminosity of any O VI line emission is stringent enough to rule out steady and simultaneous gas cooling and star formation, unlike similar systems in the Phoenix Cluster and Abell 1795. The fact that the current star formation rate differs from the local mass cooling rate is consistent with recent simulations of episodic AGN feedback and star formation in a cluster atmosphere.

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.

We present a comprehensive analysis of strong-lensing, weak-lensing shear and magnification data for a sample of 16 X-ray-regular and 4 high-magnification galaxy clusters selected from the CLASH survey. Our analysis combines constraints from 16-band HST observations and wide-field multi-color imaging taken primarily with Subaru/Suprime-Cam. We reconstruct surface mass density profiles of individual clusters from a joint analysis of the full lensing constraints, and determine masses and concentrations for all clusters. We find internal consistency of the ensemble mass calibration to be $\le 5\% \pm 6\%$ by comparison with the CLASH weak-lensing-only measurements of Umetsu et al. For the X-ray-selected subsample, we examine the concentration-mass relation and its intrinsic scatter using a Bayesian regression approach. Our model yields a mean concentration of $c|_{z=0.34} = 3.95 \pm 0.35$ at $M_{200c} \simeq 14\times 10^{14}M_\odot$ and an intrinsic scatter of $\sigma(\ln c_{200c}) = 0.13 \pm 0.06$, in excellent agreement with LCDM predictions when the CLASH selection function based on X-ray morphological regularity and the projection effects are taken into account. We also derive an ensemble-averaged surface mass density profile for the X-ray-selected subsample by stacking their individual profiles. The stacked mass profile is well described by a family of density profiles predicted for cuspy dark-matter-dominated halos, namely, the NFW, Einasto, and DARKexp models, whereas the single power-law, cored isothermal and Burkert density profiles are disfavored by the data. We show that cuspy halo models that include the two-halo term provide improved agreement with the data. For the NFW halo model, we measure a mean concentration of $c_{200c} = 3.79^{+0.30}_{-0.28}$ at $M_{200c} = 14.1^{+1.0}_{-1.0}\times 10^{14}M_\odot$, demonstrating consistency between complementary analysis methods.

We present VIMOS-VLT spectroscopy of the Frontier Fields cluster MACS~J0416.1-2403 (z=0.397). Taken as part of the CLASH-VLT survey, the large spectroscopic campaign provided more than 4000 reliable redshifts, including ~800 cluster member galaxies. The unprecedented sample of cluster members at this redshift allows us to perform a highly detailed dynamical and structural analysis of the cluster out to ~2.2$r_{200}$ (~4Mpc). Our analysis of substructures reveals a complex system composed of a main massive cluster ($M_{200}$~0.9$\times 10^{15} M_{\odot}$) presenting two major features: i) a bimodal velocity distribution, showing two central peaks separated by $\Delta V_{rf}$~1100 km s$^{-1}$ with comparable galaxy content and velocity dispersion, ii) a projected elongation of the main substructures along the NE-SW direction, with a prominent subclump ~600 kpc SW of the center and an isolated BCG approximately halfway between the center and the SW clump. We also detect a low mass structure at z~0.390, ~10' S of the cluster center, projected at ~3Mpc, with a relative line-of-sight velocity of $\Delta V_{rf}$~-1700 km s$^{-1}$. The cluster mass profile that we obtain through our dynamical analysis deviates significantly from the "universal" NFW, being best fit by a Softened Isothermal Sphere model instead. The mass profile measured from the galaxy dynamics is found to be in relatively good agreement with those obtained from strong and weak lensing, as well as with that from the X-rays, despite the clearly unrelaxed nature of the cluster. Our results reveal overall a complex dynamical state of this massive cluster and support the hypothesis that the two main subclusters are being observed in a pre-collisional phase, in line with recent findings from radio and deep X-ray data. With this article we also release the entire redshift catalog of 4386 sources in the field of this cluster.

Morphology is often used to infer the state of relaxation of galaxy clusters. The regularity, symmetry, and degree to which a cluster is centrally concentrated inform quantitative measures of cluster morphology. The Cluster Lensing and Supernova survey with Hubble Space Telescope (CLASH) used weak and strong lensing to measure the distribution of matter within a sample of 25 clusters, 20 of which were deemed to be relaxed based on their X-ray morphology and alignment of the X-ray emission with the BCG. Towards a quantitative characterization of this important sample of clusters, we present uniformly estimated X-ray morphological statistics for all 25 CLASH clusters. We compare X-ray morphologies of CLASH clusters with those identically measured for a large sample of simulated clusters from the MUSIC-2 simulations, selected by mass. We confirm a threshold in X-ray surface brightness concentration of C>0.4 for cool-core clusters, where C is the ratio of X-ray emission inside 100 kpc/h70 compared to inside 500 kpc/h70. We report and compare morphologies of these clusters inferred from Sunyaev-Zeldovich Effect (SZE) maps of the hot gas and in from projected mass maps based on strong and weak lensing. We find a strong agreement in alignments of the orientation of major axes for the lensing, X-ray, and SZE maps of nearly all of the CLASH clusters at radii of 500 kpc (approximately 0.5R500 for these clusters). We also find a striking alignment of clusters shapes at the 500 kpc scale, as measured with X-ray, SZE, and lensing, with that of the near-infrared stellar light at 10 kpc scales for the 20 "relaxed" clusters. This strong alignment indicates a powerful coupling between the cluster- and galaxy-scale galaxy formation processes.

We perform a comprehensive study of the total mass distribution of the galaxy cluster RXCJ2248 ($z=0.348$) with a set of high-precision strong lensing models, which take advantage of extensive spectroscopic information on many multiply lensed systems. In the effort to understand and quantify inherent systematics in parametric strong lensing modelling, we explore a collection of 22 models where we use different samples of multiple image families, parametrizations of the mass distribution and cosmological parameters. As input information for the strong lensing models, we use the CLASH HST imaging data and spectroscopic follow-up observations, carried out with the VIMOS and MUSE spectrographs, to identify bona-fide multiple images. A total of 16 background sources, over the redshift range $1.0-6.1$, are multiply lensed into 47 images, 24 of which are spectroscopically confirmed and belong to 10 individual sources. The cluster total mass distribution and underlying cosmology in the models are optimized by matching the observed positions of the multiple images on the lens plane. We show that with a careful selection of a sample of spectroscopically confirmed multiple images, the best-fit model reproduces their observed positions with a rms of $0.3$ in a fixed flat $\Lambda$CDM cosmology, whereas the lack of spectroscopic information lead to biases in the values of the model parameters. Allowing cosmological parameters to vary together with the cluster parameters, we find (at $68\%$ confidence level) $\Omega_m=0.25^{+0.13}_{-0.16}$ and $w=-1.07^{+0.16}_{-0.42}$ for a flat $\Lambda$CDM model, and $\Omega_m=0.31^{+0.12}_{-0.13}$ and $\Omega_\Lambda=0.38^{+0.38}_{-0.27}$ for a universe with $w=-1$ and free curvature. Using toy models mimicking the overall configuration of RXCJ2248, we estimate the impact of the line of sight mass structure on the positional rms to be $0.3\pm 0.1$.(ABRIDGED)

We developed an algorithm to find and characterize gravitationally lensed galaxies (arcs) to perform a comparison of the observed and simulated arc abundance. Observations are from the Cluster Lensing And Supernova survey with Hubble (CLASH). Simulated CLASH images are created using the MOKA package and also clusters selected from the high resolution, hydrodynamical simulations, MUSIC, over the same mass and redshift range as the CLASH sample. The algorithm' s arc elongation accuracy, completeness and false positive rate are determined and used to compute an estimate of the true arc abundance. We derive a lensing efficiency of $4 \pm 1$ arcs (with length $\ge 6"$ and length-to-width ratio $\ge 7$) per cluster for the X-ray selected CLASH sample, $4 \pm 1$ arcs per cluster for the MOKA simulated sample and $3 \pm 1$ arcs per cluster for the MUSIC simulated sample. The observed and simulated arc statistics are in full agreement. We measure the photometric redshifts of all detected arcs and find a median redshift $z_s = 1.9$ with 33% of the detected arcs having $z_s > 3$. We find that the arc abundance does not depend strongly on the source redshift distribution but is sensitive to the mass distribution of the dark matter halos (e.g. the $c-M$ relation). Our results show that consistency between the observed and simulated distributions of lensed arc sizes and axial ratios can be achieved by using cluster-lensing simulations that are carefully matched to the selection criteria used in the observations.

To ensure progress in astronomy over the coming decades, the key questions are "what facilities will we build, and when?" Toward this end, the Association of Universities for Research in Astronomy (AURA) recently commissioned a study on future space-based options for UV and optical astronomy. The resulting study - "From Cosmic Birth to Living Earths" - concluded that a space telescope equipped with a 12-meter class primary mirror would make fundamental advances across virtually all of astrophysics, including finding and characterizing the atmospheres of dozens of Earth-like planets. This ambitious telescope concept is referred to as the High Definition Space Telescope (HDST). In a recent arXiv white paper, Elvis (2015) critiqued a subset of the findings of the AURA study, focusing on the detection and characterization of rocky exoplanets in the habitable zone. In this response, we clarify these issues to confirm that HDST would play a transformative role in the study of terrestrial worlds. Its capabilities for studying exoplanets would be truly unique, even in 2035, and would complement HDST's broad and deep range of exciting astrophysics.

The CLASH X-ray selected sample of 20 galaxy clusters contains ten brightest cluster galaxies (BCGs) that exhibit significant ($>$5 $\sigma$) extinction-corrected star formation rates (SFRs). Star formation activity is inferred from photometric estimates of UV and H$\alpha$+[\ion{N}{2}] emission in knots and filaments detected in CLASH HST ACS and WFC3 observations. UV-derived SFRs in these BCGs span two orders of magnitude, including two with a SFR $\gtrsim$ 100 M$_{\odot}$ yr$^{-1}$. These measurements are supplemented with [\ion{O}{2}], [\ion{O}{3}], and H$\beta$ fluxes measured from spectra obtained with the SOAR telescope. We confirm that photoionization from ongoing star formation powers the line emission nebulae in these BCGs, although in many BCGs there is also evidence of a LINER-like contribution to the line emission. Coupling these data with Chandra X-ray measurements, we infer that the star formation occurs exclusively in low-entropy cluster cores and exhibits a correlation with gas properties related to cooling. We also perform an in-depth study of the starburst history of the BCG in the cluster RXJ1532.9+3021, and create 2D maps of stellar properties on scales down to $\sim$350 pc. These maps reveal evidence for an ongoing burst occurring in elongated filaments, generally on $\sim$ 0.5-1.0 Gyr timescales, although some filaments are consistent with much younger ($\lesssim$ 100 Myr) burst timescales and may be correlated with recent activity from the AGN. The relationship between BCG SFRs and the surrounding ICM gas properties provide new support for the process of feedback-regulated cooling in galaxy clusters and is consistent with recent theoretical predictions.

For the first time in history, humans have reached the point where it is possible to construct a revolutionary space-based observatory that has the capability to find dozens of Earth-like worlds, and possibly some with signs of life. This same telescope, designed as a long-lived facility, would also produce transformational scientific advances in every area of astronomy and astrophysics from black hole physics to galaxy formation, from star and planet formation to the origins of the Solar System. The Association of Universities for Research in Astronomy (AURA) commissioned a study on a next-generation UVOIR space observatory with the highest possible scientific impact in the era following JWST. This community-based study focuses on the future space-based options for UV and optical astronomy that significantly advance our understanding of the origin and evolution of the cosmos and the life within it. The committee concludes that a space telescope equipped with a 12-meter class primary mirror can find and characterize dozens of Earth-like planets and make fundamental advances across nearly all fields of astrophysics. The concept is called the High Definition Space Telescope (HDST). The telescope would be located at the Sun-Earth L2 point and would cover a spectral range that, at a minimum, runs from 0.1 to 2 microns. Unlike JWST, HDST will not need to operate at cryogenic temperatures. HDST can be made to be serviceable on orbit but does not require servicing to complete its primary scientific objectives. We present the scientific and technical requirements for HDST and show that it could allow us to determine whether or not life is common outside the Solar System. We do not propose a specific design for such a telescope, but show that designing, building and funding such a facility is feasible beginning in the next decade - if the necessary strategic investments in technology begin now.

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.

Brightest cluster galaxies (BCGs) are usually quiescent, but many exhibit star formation. Here we exploit the opportunity provided by rest-frame UV imaging of galaxy clusters in the CLASH (Cluster Lensing and Supernovae with Hubble) Multi-Cycle Treasury Project to reveal the diversity of UV morphologies in BCGs and to compare them with recent simulations of the cool, star-forming gas structures produced by precipitation-driven feedback. All of the CLASH BCGs are detected in the rest-frame UV (280 nm), regardless of their star-formation activity, because evolved stellar populations produce a modest amount of UV light that traces the relatively smooth, symmetric, and centrally peaked stellar distribution seen in the near infrared. Ultraviolet morphologies among the BCGs with strong UV excesses exhibit distinctive knots, multiple elongated clumps, and extended filaments of emission that distinctly differ from the smooth profiles of the UV-quiet BCGs. These structures, which are similar to those seen in the few star-forming BCGs observed in the UV at low redshift, are suggestive of bi-polar streams of clumpy star formation, but not of spiral arms or large, kpc-scale disks. Based on the number of streams and lack of culprit companion galaxies, these streams are unlikely to have arisen from multiple collisions with gas-rich galaxies. These star-forming UV structures are morphologically similar to the cold-gas structures produced in simulations of precipitation-driven AGN feedback in which jets uplift low-entropy gas to greater altitudes, causing it to condense. Unobscured star-formation rates estimated from CLASH UV images using the Kennicutt relation range up to 80 solar masses per year in the most extended and highly structured systems. The circumgalactic gas-entropy threshold for star formation in CLASH BCGs at z~0.2-0.5 is indistinguishable from that for clusters at z < 0.2.

In the effort to understand the link between the structure of galaxy clusters and their galaxy populations, we focus on MACSJ1206.2-0847 at z~0.44 and probe its substructure in the projected phase space through the spectrophotometric properties of a large number of galaxies from the CLASH-VLT survey. Our analysis is mainly based on an extensive spectroscopic dataset of 445 member galaxies, mostly acquired with VIMOS@VLT as part of our ESO Large Programme, sampling the cluster out to a radius ~2R200 (4 Mpc). We classify 412 galaxies as passive, with strong Hdelta absorption (red and blue galaxies, and with emission lines from weak to very strong. A number of tests for substructure detection are applied to analyze the galaxy distribution in the velocity space, in 2D space, and in 3D projected phase-space. Studied in its entirety, the cluster appears as a large-scale relaxed system with a few secondary, minor overdensities in 2D distribution. We detect no velocity gradients or evidence of deviations in local mean velocities. The main feature is the WNW-ESE elongation. The analysis of galaxy populations per spectral class highlights a more complex scenario. The passive galaxies and red strong Hdelta galaxies trace the cluster center and the WNW-ESE elongated structure. The red strong Hdelta galaxies also mark a secondary, dense peak ~2 Mpc at ESE. The emission line galaxies cluster in several loose structures, mostly outside R200. The observational scenario agrees with MACS J1206.2-0847 having WNW-ESE as the direction of the main cluster accretion, traced by passive galaxies and red strong Hdelta galaxies. The red strong Hdelta galaxies, interpreted as poststarburst galaxies, date a likely important event 1-2 Gyr before the epoch of observation. The emission line galaxies trace a secondary, ongoing infall where groups are accreted along several directions.

This report describes the 2014 study by the Science Definition Team (SDT) of the Wide-Field Infrared Survey Telescope (WFIRST) mission. It is a space observatory that will address the most compelling scientific problems in dark energy, exoplanets and general astrophysics using a 2.4-m telescope with a wide-field infrared instrument and an optical coronagraph. The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope as its top priority for a new large space mission. As conceived by the decadal survey, WFIRST would carry out a dark energy science program, a microlensing program to determine the demographics of exoplanets, and a general observing program utilizing its ultra wide field. In October 2012, NASA chartered a Science Definition Team (SDT) to produce, in collaboration with the WFIRST Study Office at GSFC and the Program Office at JPL, a Design Reference Mission (DRM) for an implementation of WFIRST using one of the 2.4-m, Hubble-quality telescope assemblies recently made available to NASA. This DRM builds on the work of the earlier WFIRST SDT, reported by Green et al. (2012) and the previous WFIRST-2.4 DRM, reported by Spergel et. (2013). The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the 1.3-m and 1.1-m designs considered previously, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The addition of an on-axis coronagraphic instrument to the baseline design enables imaging and spectroscopic studies of planets around nearby stars.

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

The focus of this report is on the opportunities enabled by the combination of LSST, Euclid and WFIRST, the optical surveys that will be an essential part of the next decade's astronomy. The sum of these surveys has the potential to be significantly greater than the contributions of the individual parts. As is detailed in this report, the combination of these surveys should give us multi-wavelength high-resolution images of galaxies and broadband data covering much of the stellar energy spectrum. These stellar and galactic data have the potential of yielding new insights into topics ranging from the formation history of the Milky Way to the mass of the neutrino. However, enabling the astronomy community to fully exploit this multi-instrument data set is a challenging technical task: for much of the science, we will need to combine the photometry across multiple wavelengths with varying spectral and spatial resolution. We identify some of the key science enabled by the combined surveys and the key technical challenges in achieving the synergies.

We present results from a comprehensive lensing analysis in HST data, of the complete CLASH cluster sample. We identify new multiple-images previously undiscovered allowing improved or first constraints on the cluster inner mass distributions and profiles. We combine these strong-lensing constraints with weak-lensing shape measurements within the HST FOV to jointly constrain the mass distributions. The analysis is performed in two different common parameterizations (one adopts light-traces-mass for both galaxies and dark matter while the other adopts an analytical, elliptical NFW form for the dark matter), to provide a better assessment of the underlying systematics - which is most important for deep, cluster-lensing surveys, especially when studying magnified high-redshift objects. We find that the typical (median), relative systematic differences throughout the central FOV are $\sim40\%$ in the (dimensionless) mass density, $\kappa$, and $\sim20\%$ in the magnification, $\mu$. We show maps of these differences for each cluster, as well as the mass distributions, critical curves, and 2D integrated mass profiles. For the Einstein radii ($z_{s}=2$) we find that all typically agree within $10\%$ between the two models, and Einstein masses agree, typically, within $\sim15\%$. At larger radii, the total projected, 2D integrated mass profiles of the two models, within $r\sim2\arcmin$, differ by $\sim30\%$. Stacking the surface-density profiles of the sample from the two methods together, we obtain an average slope of $d\log (\Sigma)/d\log(r)\sim-0.64\pm0.1$, in the radial range [5,350] kpc. Lastly, we also characterize the behavior of the average magnification, surface density, and shear differences between the two models, as a function of both the radius from the center, and the best-fit values of these quantities.

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 a detailed mass reconstruction and a novel study on the substructure properties in the core of the CLASH and Frontier Fields galaxy cluster MACS J0416.1-2403. We show and employ our extensive spectroscopic data set taken with the VIMOS instrument as part of our CLASH-VLT program, to confirm spectroscopically 10 strong lensing systems and to select a sample of 175 plausible cluster members to a limiting stellar mass of log(M_*/M_Sun) ~ 8.6. We reproduce the measured positions of 30 multiple images with a remarkable median offset of only 0.3" by means of a comprehensive strong lensing model comprised of 2 cluster dark-matter halos, represented by cored elliptical pseudo-isothermal mass distributions, and the cluster member components. The latter have total mass-to-light ratios increasing with the galaxy HST/WFC3 near-IR (F160W) luminosities. The measurement of the total enclosed mass within the Einstein radius is accurate to ~5%, including systematic uncertainties. We emphasize that the use of multiple-image systems with spectroscopic redshifts and knowledge of cluster membership based on extensive spectroscopic information is key to constructing robust high-resolution mass maps. We also produce magnification maps over the central area that is covered with HST observations. We investigate the galaxy contribution, both in terms of total and stellar mass, to the total mass budget of the cluster. When compared with the outcomes of cosmological $N$-body simulations, our results point to a lack of massive subhalos in the inner regions of simulated clusters with total masses similar to that of MACS J0416.1-2403. Our findings of the location and shape of the cluster dark-matter halo density profiles and on the cluster substructures provide intriguing tests of the assumed collisionless, cold nature of dark matter and of the role played by baryons in the process of structure formation.

We use velocity dispersion measurements of 21 individual cluster members in the core of Abell 383, obtained with MMT Hectospec, to separate the galaxy and the smooth dark halo (DH) lensing contributions. While lensing usually constrains the overall, projected mass density, the innovative use of velocity dispersion measurements as a proxy for masses of individual cluster members breaks inherent degeneracies and allows us to (a) refine the constraints on single galaxy masses and on the galaxy mass-to-light scaling relation and, as a result, (b) refine the constraints on the DM-only map, a high-end goal of lens modelling. The knowledge of cluster member velocity dispersions improves the fit by 17% in terms of the image reproduction $\chi^2$, or 20% in terms of the rms. The constraints on the mass parameters improve by ~10% for the DH, while for the galaxy component, they are refined correspondingly by ~50%, including the galaxy halo truncation radius. For an L$^*$ galaxy with M$^*_B$=-20.96, for example, we obtain best fitting truncation radius r$^*_{tr}=20.5^{+9.6}_{-6.7}$ kpc and velocity dispersion $\sigma^*=324\pm17 km/s$. Moreover, by performing the surface brightness reconstruction of the southern giant arc, we improve the constraints on r$_{tr}$ of two nearby cluster members, which have measured velocity dispersions, by more than ~30%. We estimate the stripped mass for these two galaxies, getting results that are consistent with numerical simulations. In the future, we plan to apply this analysis to other galaxy clusters for which velocity dispersions of member galaxies are available.

We have observed 433 z<=0.08 brightest cluster galaxies (BCGs) in a full-sky survey of Abell clusters. The BCG Hubble diagram is consistent to within 2% of a Omega_m=0.3, Lambda=0.7 Hubble relation. The L_m-alpha relation for BCGs, which uses alpha, the log-slope of the BCG photometric curve of growth, to predict metric luminosity, L_m, has 0.27 mag residuals. We measure central stellar velocity dispersions, sigma, of the BCGs, finding the Faber-Jackson relation to flatten as the metric aperture grows to include an increasing fraction of the total BCG luminosity. A 3-parameter "metric plane" relation using alpha and sigma together gives the best prediction of L_m, with 0.21 mag residuals. The projected spatial offset, r_x, of BCGs from the X-ray-defined cluster center is a gamma=-2.33 power-law over 1<r_x<10^3 kpc. The median offset is ~10 kpc, but ~15% of the BCGs have r_x>100 kpc. The absolute cluster-dispersion normalized BCG peculiar velocity |Delta V_1|/sigma_c follows an exponential distribution with scale length 0.39+/-0.03. Both L_m and alpha increase with sigma_c. The alpha parameter is further moderated by both the spatial and velocity offset from the cluster center, with larger alpha correlated with the proximity of the BCG to the cluster mean velocity or potential center. At the same time, position in the cluster has little effect on L_m. The luminosity difference between the BCG and second-ranked galaxy, M2, increases as the peculiar velocity of the BCG within the cluster decreases. Further, when M2 is a close luminosity "rival" of the BCG, the galaxy that is closest to either the velocity or X-ray center of the cluster is most likely to have the larger alpha. We conclude that the inner portions of the BCGs are formed outside the cluster, but interactions in the heart of the galaxy cluster grow and extend the envelopes of the BCGs.

We present a joint shear-and-magnification weak-lensing analysis of a sample of 16 X-ray-regular and 4 high-magnification galaxy clusters at 0.19<z<0.69 selected from the Cluster Lensing And Supernova survey with Hubble (CLASH). Our analysis uses wide-field multi-color imaging, taken primarily with Suprime-Cam on the Subaru Telescope. From a stacked shear-only analysis of the X-ray-selected subsample, we detect the ensemble-averaged lensing signal with a total signal-to-noise ratio of ~25 in the radial range of 200 to 3500kpc/h. The stacked tangential-shear signal is well described by a family of standard density profiles predicted for dark-matter-dominated halos in gravitational equilibrium, namely the Navarro-Frenk-White (NFW), truncated variants of NFW, and Einasto models. For the NFW model, we measure a mean concentration of $c_{200c}=4.01^{+0.35}_{-0.32}$ at $M_{200c}=1.34^{+0.10}_{-0.09} 10^{15}M_{\odot}$. We show this is in excellent agreement with Lambda cold-dark-matter (LCDM) predictions when the CLASH X-ray selection function and projection effects are taken into account. The best-fit Einasto shape parameter is $\alpha_E=0.191^{+0.071}_{-0.068}$, which is consistent with the NFW-equivalent Einasto parameter of $\sim 0.18$. We reconstruct projected mass density profiles of all CLASH clusters from a joint likelihood analysis of shear-and-magnification data, and measure cluster masses at several characteristic radii. We also derive an ensemble-averaged total projected mass profile of the X-ray-selected subsample by stacking their individual mass profiles. The stacked total mass profile, constrained by the shear+magnification data, is shown to be consistent with our shear-based halo-model predictions including the effects of surrounding large-scale structure as a two-halo term, establishing further consistency in the context of the LCDM model.

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.

Context. The study of the galaxy stellar mass function (SMF) in relation to the galaxy environment and the stellar mass density profile, rho(r), is a powerful tool to constrain models of galaxy evolution. Aims. We determine the SMF of the z=0.44 cluster of galaxies MACS J1206.2-0847 separately for passive and star-forming (SF) galaxies, in different regions of the cluster, from the center out to approximately 2 virial radii. We also determine rho(r) to compare it to the number density and total mass density profiles. Methods. We use the dataset from the CLASH-VLT survey. Stellar masses are obtained by SED fitting on 5-band photometric data obtained at the Subaru telescope. We identify 1363 cluster members down to a stellar mass of 10^9.5 Msolar. Results. The whole cluster SMF is well fitted by a double Schechter function. The SMFs of cluster SF and passive galaxies are statistically different. The SMF of the SF cluster galaxies does not depend on the environment. The SMF of the passive population has a significantly smaller slope (in absolute value) in the innermost (<0.50 Mpc), highest density cluster region, than in more external, lower density regions. The number ratio of giant/subgiant galaxies is maximum in this innermost region and minimum in the adjacent region, but then gently increases again toward the cluster outskirts. This is also reflected in a decreasing radial trend of the average stellar mass per cluster galaxy. On the other hand, the stellar mass fraction, i.e., the ratio of stellar to total cluster mass, does not show any significant radial trend. Conclusions. Our results appear consistent with a scenario in which SF galaxies evolve into passive galaxies due to density-dependent environmental processes, and eventually get destroyed very near the cluster center to become part of a diffuse intracluster medium.

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

We report the LBT/MODS1 spectroscopic confirmation of two images of faint Lyman alpha emitters at $z=6.4$ behind the Frontier Fields galaxy cluster MACSJ0717.5+3745. A wide range of lens models suggests that the two images are highly magnified, with a strong lower limit of mu>5. These are the faintest z>6 candidates spectroscopically confirmed to date. These may be also multiple images of the same z=6.4 source as supported by their similar intrinsic properties, but the lens models are inconclusive regarding this interpretation. To be cautious, we derive the physical properties of each image individually. Thanks to the high magnification, the observed near-infrared (restframe ultraviolet) part of the spectral energy distributions and Ly-alpha lines are well detected with S/N(m_1500)>~10 and S/N(Ly-alpha)~10-15. Adopting mu>5, the absolute magnitudes, M_1500, and Ly-alpha fluxes, are fainter than -18.7 and 2.8x10^(-18)erg/s/cm2, respectively. We find a very steep ultraviolet spectral slope beta=-3.0+/-0.5 (F_lambda=lambda^(beta)), implying that these are very young, dust-free and low metallicity objects, made of standard stellar populations or even extremely metal poor stars (age<~30Myr, E(B-V)=0 and metallicity 0.0-0.2 Z/Zsolar). The objects are compact (< 1 kpc^(2)), and with a stellar mass M* < 10^(8) M_solar. The very steep beta, the presence of the Ly-alpha line and the intrinsic FWHM (<300 km/s) of these newborn objects do not exclude a possible leakage of ionizing radiation. We discuss the possibility that such faint galaxies may resemble those responsible for cosmic reionization.

A pressureless scenario for the Dark Matter (DM) fluid is a widely adopted hypothesis, despite the absence of a direct observational evidence. According to general relativity, the total mass-energy content of a system shapes the gravitational potential well, but different test particles perceive this potential in different ways depending on their properties. Cluster galaxy velocities, being $\ll$c, depend solely on the gravitational potential, whereas photon trajectories reflect the contributions from the gravitational potential plus a relativistic-pressure term that depends on the cluster mass. We exploit this phenomenon to constrain the Equation of State (EoS) parameter of the fluid, primarily DM, contained in galaxy clusters. We use the complementary information provided by the kinematic and lensing mass profiles of the galaxy cluster MACS 1206.2-0847 at $z=0.44$, as obtained in an extensive imaging and spectroscopic campaign within the CLASH survey. The unprecedented high quality of our data-set and the properties of this cluster are well suited to determine the EoS parameter of the cluster fluid. Since baryons contribute at most $15\%$ to the total mass in clusters and their pressure is negligible, the EoS parameter we derive describes the behavior of the DM fluid. We obtain the most stringent constraint on the DM EoS parameter to date, $w=(p_r+2\,p_t)/(3\,c^2\rho)=0.00\pm0.15\mathrm{(stat)}\pm0.08\mathrm{(syst)}$, averaged over the radial range $0.5\,\mathrm{Mpc}\leq$$r$$\leq$$r_{200}$, where $p_r$ and $p_t$ are the radial and tangential pressure, and $\rho$ is the density. We plan to further improve our constraint by applying the same procedure to all clusters from the ongoing CLASH-VLT program.

We present the supernova (SN) sample and Type-Ia SN (SN Ia) rates from the Cluster Lensing And Supernova survey with Hubble (CLASH). Using the Advanced Camera for Surveys and the Wide Field Camera 3 on the Hubble Space Telescope (HST), we have imaged 25 galaxy-cluster fields and parallel fields of non-cluster galaxies. We report a sample of 27 SNe discovered in the parallel fields. Of these SNe, ~13 are classified as SN Ia candidates, including four SN Ia candidates at redshifts z > 1.2. We measure volumetric SN Ia rates to redshift 1.8 and add the first upper limit on the SN Ia rate in the range 1.8 < z < 2.4. The results are consistent with the rates measured by the HST/GOODS and Subaru Deep Field SN surveys. We model these results together with previous measurements at z < 1 from the literature. The best-fitting SN Ia delay-time distribution (DTD; the distribution of times that elapse between a short burst of star formation and subsequent SN Ia explosions) is a power law with an index of -1.00 +0.06(0.09) -0.06(0.10) (statistical) +0.12 -0.08 (systematic), where the statistical uncertainty is a result of the 68% and 95% (in parentheses) statistical uncertainties reported for the various SN Ia rates (from this work and from the literature), and the systematic uncertainty reflects the range of possible cosmic star-formation histories. We also test DTD models produced by an assortment of published binary population synthesis (BPS) simulations. The shapes of all BPS double-degenerate DTDs are consistent with the volumetric SN Ia measurements, when the DTD models are scaled up by factors of 3-9. In contrast, all BPS single-degenerate DTDs are ruled out by the measurements at a >99% significance level.

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

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]

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)

Aims. We investigate the degree of improvement in dark energy constraints that can be achieved by extending Type Ia Supernova (SN Ia) samples to redshifts z > 1.5 with the Hubble Space Telescope (HST), particularly in the ongoing CANDELS and CLASH multi-cycle treasury programs. Methods. Using the popular CPL parametrization of the dark energy, w = w0 +wa(1-a), we generate mock SN Ia samples that can be projected out to higher redshifts. The synthetic datasets thus generated are fitted to the CPL model, and we evaluate the improvement that a high-z sample can add in terms of ameliorating the statistical and systematic uncertainties on cosmological parameters. Results. In an optimistic but still very achievable scenario, we find that extending the HST sample beyond CANDELS+CLASH to reach a total of 28 SN Ia at z > 1.0 could improve the uncertainty in the wa parameter by up to 21%. The corresponding improvement in the figure of merit (FoM) would be as high as 28%. Finally, we consider the use of high-redshift SN Ia samples to detect non-cosmological evolution in SN Ia luminosities with redshift, finding that such tests could be undertaken by future spacebased infrared surveys using the James Webb Space Telescope (JWST).

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.

The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope (WFIRST) as its top priority for a new large space mission. The report of the WFIRST-AFTA Science Definition Team (SDT) presents a Design Reference Mission for WFIRST that employs one of the 2.4-m, Hubble-quality mirror assemblies recently made available to NASA. The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the smaller aperture designs previously considered for WFIRST, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The option of adding an on-axis, coronagraphic instrument would enable imaging and spectroscopic studies of planets around nearby stars. This short article, produced as a companion to the SDT report, summarizes the key points of the WFIRST-2.4 DRM. It highlights the remarkable opportunity that the 2.4-m telescope affords for advances in many fields of astrophysics and cosmology, including dark energy, the demographics and characterization of exoplanets, the evolution of galaxies and quasars, and the stellar populations of the Milky Way and its neighbors.

The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope (WFIRST) as its top priority for a new large space mission. As conceived by the decadal survey, WFIRST would carry out a dark energy science program, a microlensing program to determine the demographics of exoplanets, and a general observing program utilizing its ultra wide field. In October 2012, NASA chartered a Science Definition Team (SDT) to produce, in collaboration with the WFIRST Project Office at GSFC and the Program Office at JPL, a Design Reference Mission (DRM) for an implementation of WFIRST using one of the 2.4-m, Hubble-quality mirror assemblies recently made available to NASA. This DRM builds on the work of the earlier WFIRST SDT, reported by Green et al. (2012). The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the 1.3-m and 1.1-m designs considered previously, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The option of adding an on-axis, coronagraphic instrument would enable imaging and spectroscopic studies of planets around nearby stars. This document presents the final report of the SDT.

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 discuss scientific, technical and programmatic issues related to the use of an NRO 2.4m telescope for the WFIRST initiative of the 2010 Decadal Survey. We show that this implementation of WFIRST, which we call "NEW WFIRST," would achieve the goals of the NWNH Decadal Survey for the WFIRST core programs of Dark Energy and Microlensing Planet Finding, with the crucial benefit of deeper and/or wider near-IR surveys for GO science and a potentially Hubble-like Guest Observer program. NEW WFIRST could also include a coronagraphic imager for direct detection of dust disks and planets around neighboring stars, a high-priority science and technology precursor for future ambitious programs to image Earth-like planets around neighboring stars.

If we are to develop a comprehensive and predictive theory of galaxy formation and evolution, it is essential that we obtain an accurate assessment of how and when galaxies assemble their stellar populations, and how this assembly varies with environment. There is strong observational support for the hierarchical assembly of galaxies, but our insight into this assembly comes from sifting through the resolved field populations of the surviving galaxies we see today, in order to reconstruct their star formation histories, chemical evolution, and kinematics. To obtain the detailed distribution of stellar ages and metallicities over the entire life of a galaxy, one needs multi-band photometry reaching solar-luminosity main sequence stars. The Hubble Space Telescope can obtain such data in the low-density regions of Local Group galaxies. To perform these essential studies for a fair sample of the Local Universe, we will require observational capabilities that allow us to extend the study of resolved stellar populations to much larger galaxy samples that span the full range of galaxy morphologies, while also enabling the study of the more crowded regions of relatively nearby galaxies. With such capabilities in hand, we will reveal the detailed history of star formation and chemical evolution in the universe.

We report statistical results for dark matter (DM) velocity anisotropy, \beta, from a sample of some 6000 cluster-size halos (at redshift zero) identified in a \Lambda CDM hydrodynamical adaptive mesh refinement simulation performed with the Enzo code. These include profiles of \beta\ in clusters with different masses, relaxation states, and at several redshifts, modeled both as spherical and triaxial DM configurations. Specifically, although we find a large scatter in the DM velocity anisotropy profiles of different halos (across elliptical shells extending to at least ~$1.5 r_{vir}$), universal patterns are found when these are averaged over halo mass, redshift, and relaxation stage. These are characterized by a very small velocity anisotropy at the halo center, increasing outward to about 0.27 and leveling off at about $0.2 r_{vir}$. Indirect measurements of the DM velocity anisotropy fall on the upper end of the theoretically expected range. Though measured indirectly, the estimations are derived by using two different surrogate measurements - X-ray and galaxy dynamics. Current estimates of the DM velocity anisotropy are based on very small cluster sample. Increasing this sample will allow testing theoretical predictions, including the speculation that the decay of DM particles results in a large velocity boost. We also find, in accord with previous works, that halos are triaxial and likely to be more prolate when unrelaxed, whereas relaxed halos are more likely to be oblate. Our analysis does not indicate that there is significant correlation (found in some previous studies) between the radial density slope, \gamma, and \beta\ at large radii, $0.3 r_{vir} < r < r_{vir}$.

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 confirm the detection of 3 groups in the Lynx supercluster, at z~1.3, and give their redshifts and masses. We study the properties of the group galaxies as compared to the central clusters, RXJ0849+4452 and RXJ0848+4453, selecting 89 galaxies in the clusters and 74 galaxies in the groups. We morphologically classify galaxies by visual inspection, noting that our early-type galaxy (ETG) sample would have been contaminated at the 30% -40% level by simple automated classification methods (e.g. based on Sersic index). In luminosity selected samples, both clusters and groups show high fractions of Sa galaxies. The ETG fractions never rise above ~50% in the clusters, which is low compared to the fractions observed in clusters at z~1. However, ETG plus Sa fractions are similar to those observed for ETGs in clusters at z~1. Bulge-dominated galaxies visually classified as Sas might also be ETGs with tidal features or merger remnants. They are mainly red and passive, and span a large range in luminosity. Their star formation seems to have been quenched before experiencing a morphological transformation. Because their fraction is smaller at lower redshifts, they might be the spiral population that evolves into ETGs. For mass-selected samples, the ETG fraction show no significant evolution with respect to local clusters, suggesting that morphological transformations occur at lower masses and densities. The ETG mass-size relation shows evolution towards smaller sizes at higher redshift in both clusters and groups, while the late-type mass-size relation matches that observed locally. The group ETG red sequence shows lower zero points and larger scatters than in clusters, both expected to be an indication of a younger galaxy population. The estimated age difference is small when compared to the difference in age at different galaxy masses.

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

Using the sample of Type Ia supernovae (SNe Ia) discovered by the Hubble Space Telescope (HST) Cluster Supernova Survey and augmented with HST-observed SNe Ia in the GOODS fields, we search for correlations between the properties of SNe and their host galaxies at high redshift. We use galaxy color and quantitative morphology to determine the red sequence in 25 clusters and develop a model to distinguish passively evolving early-type galaxies from star-forming galaxies in both clusters and the field. With this approach, we identify six SN Ia hosts that are early-type cluster members and eleven SN Ia hosts that are early-type field galaxies. We confirm for the first time at z>0.9 that SNe Ia hosted by early-type galaxies brighten and fade more quickly than SNe Ia hosted by late-type galaxies. We also show that the two samples of hosts produce SNe Ia with similar color distributions. The relatively simple spectral energy distributions (SEDs) expected for passive galaxies enable us to measure stellar masses of early-type SN hosts. In combination with stellar mass estimates of late-type GOODS SN hosts from Thomson & Chary (2011), we investigate the correlation of host mass with Hubble residual observed at lower redshifts. Although the sample is small and the uncertainties are large, a hint of this relation is found at z>0.9. By simultaneously fitting the average cluster galaxy formation history and dust content to the red-sequence scatters, we show that the reddening of early-type cluster SN hosts is likely E(B-V) <~ 0.06. The similarity of the field and cluster early-type host samples suggests that field early-type galaxies that lie on the red sequence may also be minimally affected by dust. Hence, the early-type hosted SNe Ia studied here occupy a more favorable environment to use as well-characterized high-redshift standard candles than other SNe Ia.

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.

We present the Kormendy and mass-size relations for early-type galaxies (ETGs) as a function of environment at z~1.3. Our sample includes 76 visually classified ETGs with masses 10^10 < M/Msun < 10^11.5, selected in the Lynx supercluster and in the GOODS/CDF-S field, 31 ETGs in clusters, 18 in groups and 27 in the field, all with multi-wavelength photometry and HST/ACS observations. The Kormendy relation, in place at z~1.3, does not depend on the environment. The mass-size relation reveals that ETGs overall appear to be more compact in denser environments: cluster ETGs have sizes on average around 30-50% smaller than those of the local universe, and a distribution with a smaller scatter, whereas field ETGs show a mass-size relation with a similar distribution than the local one. Our results imply that (1) the mass-size relation in the field did not evolve overall from z ~ 1.3 to present; this is interesting and in contrast to the trend found at higher masses from previous works; (2) in denser environments, either ETGs have increased their size by 30-50%, on average, and spread their distributions, or more ETGs have been formed within the dense environment from not ETG progenitors or larger galaxies have been accreted to a pristine compact population to reproduce the mass-size relation observed in the local Universe. Our results are driven by galaxies with masses M<2*10^11Msun and those with masses M~10^11Msun follow the same trends that the entire sample. Following Valentinuzzi et al. definition of superdense ETGs, around 35-45% of our cluster sample is made of superdense ETGs.

We report a measurement of the Type Ia supernova (SN Ia) rate in galaxy clusters at 0.9 < z < 1.45 from the Hubble Space Telescope (HST) Cluster Supernova Survey. This is the first cluster SN Ia rate measurement with detected z > 0.9 SNe. Finding 8 +/- 1 cluster SNe Ia, we determine a SN Ia rate of 0.50 +0.23-0.19 (stat) +0.10-0.09 (sys) SNuB (SNuB = 10^-12 SNe L_{sun,B}^-1 yr^-1). In units of stellar mass, this translates to 0.36 +0.16-0.13 (stat) +0.07-0.06 (sys) SNuM (SNuM = 10^-12 SNe M_sun^-1 yr^-1). This represents a factor of approximately 5 +/- 2 increase over measurements of the cluster rate at z < 0.2. We parameterize the late-time SN Ia delay time distribution with a power law (proportional to t^s). Under the assumption of a cluster formation redshift of z_f = 3, our rate measurement in combination with lower-redshift cluster SN Ia rates constrains s = -1.41 +0.47/-0.40, consistent with measurements of the delay time distribution in the field. This measurement is generally consistent with expectations for the "double degenerate" scenario and inconsistent with some models for the "single degenerate" scenario predicting a steeper delay time distribution at large delay times. We check for environmental dependence and the influence of younger stellar populations by calculating the rate specifically in cluster red-sequence galaxies and in morphologically early-type galaxies, finding results similar to the full cluster rate. Finally, the upper limit of one host-less cluster SN Ia detected in the survey implies that the fraction of stars in the intra-cluster medium is less than 0.47 (95% confidence), consistent with measurements at lower redshifts.

We present a measurement of the volumetric Type Ia supernova (SN Ia) rate out to z ~ 1.6 from the Hubble Space Telescope Cluster Supernova Survey. In observations spanning 189 orbits with the Advanced Camera for Surveys we discovered 29 SNe, of which approximately 20 are SNe Ia. Twelve of these SNe Ia are located in the foregrounds and backgrounds of the clusters targeted in the survey. Using these new data, we derive the volumetric SN Ia rate in four broad redshift bins, finding results consistent with previous measurements at z > 1 and strengthening the case for a SN Ia rate that is equal to or greater than ~0.6 x 10^-4/yr/Mpc^3 at z ~ 1 and flattening out at higher redshift. We provide SN candidates and efficiency calculations in a form that makes it easy to rebin and combine these results with other measurements for increased statistics. Finally, we compare the assumptions about host-galaxy dust extinction used in different high-redshift rate measurements, finding that different assumptions may induce significant systematic differences between measurements.

The VAO (Virtual Astronomical Observatory) Science Council (VAO-SC) met on July 27-28, 2011 at the Harvard-Smithsonian Center for Astrophysics in Cambridge MA, to review the VAO performance during its first year of operations. In this meeting the VAO demonstrated the new tools for astronomers that are being released in September 2011 and presented plans for the second year of activities, resulting from studies conducted during the first year. This document contains the recommendations of the VAO-SC for the second year of activity of the VAO.

We present weak gravitational lensing analysis of 22 high-redshift (z >~1) clusters based on Hubble Space Telescope images. Most clusters in our sample provide significant lensing signals and are well detected in their reconstructed two-dimensional mass maps. Combining the current results and our previous weak-lensing studies of five other high-z clusters, we compare gravitational lensing masses of these clusters with other observables. We revisit the question whether the presence of the most massive clusters in our sample is in tension with the current LambdaCDM structure formation paradigm. We find that the lensing masses are tightly correlated with the gas temperatures and establish, for the first time, the lensing mass-temperature relation at z >~ 1. For the power law slope of the M-TX relation (M propto T^{\alpha}), we obtain \alpha=1.54 +/- 0.23. This is consistent with the theoretical self-similar prediction \alpha=3/2 and with the results previously reported in the literature for much lower redshift samples. However, our normalization is lower than the previous results by 20-30%, indicating that the normalization in the M-TX relation might evolve. After correcting for Eddington bias and updating the discovery area with a more conservative choice, we find that the existence of the most massive clusters in our sample still provides a tension with the current Lambda CDM model. The combined probability of finding the four most massive clusters in this sample after marginalization over current cosmological parameters is less than 1%.

We outline our methods for obtaining high precision mass profiles, combining independent weak-lensing distortion, magnification, and strong-lensing measurements. For massive clusters the strong and weak lensing regimes contribute equal logarithmic coverage of the radial profile. The utility of high-quality data is limited by the cosmic noise from large scale structure along the line of sight. This noise is overcome when stacking clusters, as too are the effects of cluster asphericity and substructure, permitting a stringent test of theoretical models. We derive a mean radial mass profile of four similar mass clusters of high-quality HST and Subaru images, in the range R=40kpc/h to 2800kpc/h, where the inner radial boundary is sufficiently large to avoid smoothing from miscentering effects. The stacked mass profile is detected at 58-sigma significance over the entire radial range, with the contribution from the cosmic noise included. We show that the projected mass profile has a continuously steepening gradient out to beyond the virial radius, in remarkably good agreement with the standard Navarro-Frenk-White form predicted for the family of CDM-dominated halos in gravitational equilibrium. The central slope is constrained to lie in the range, -dln{\rho}/dln{r}=0.89^{+0.27}_{-0.39}. The mean concentration is c_{vir}=7.68^{+0.42}_{-0.40} (at a mean virial mass 1.54^{+0.11}_{-0.10}\times 10^{15} M_{sun}/h), which is high for relaxed, high-mass clusters, but consistent with LCDM when a sizable projection bias estimated from N-body simulations is considered. This possible tension will be more definitively explored with new cluster surveys, such as CLASH, LoCuSS, Subaru HSC, and XXM-XXL, to construct the c-M relation over a wider mass range.

We present ACS, NICMOS, and Keck AO-assisted photometry of 20 Type Ia supernovae SNe Ia from the HST Cluster Supernova Survey. The SNe Ia were discovered over the redshift interval 0.623 < z < 1.415. Fourteen of these SNe Ia pass our strict selection cuts and are used in combination with the world's sample of SNe Ia to derive the best current constraints on dark energy. Ten of our new SNe Ia are beyond redshift $z=1$, thereby nearly doubling the statistical weight of HST-discovered SNe Ia beyond this redshift. Our detailed analysis corrects for the recently identified correlation between SN Ia luminosity and host galaxy mass and corrects the NICMOS zeropoint at the count rates appropriate for very distant SNe Ia. Adding these supernovae improves the best combined constraint on the dark energy density \rho_{DE}(z) at redshifts 1.0 < z < 1.6 by 18% (including systematic errors). For a LambdaCDM universe, we find \Omega_\Lambda = 0.724 +0.015/-0.016 (68% CL including systematic errors). For a flat wCDM model, we measure a constant dark energy equation-of-state parameter w = -0.985 +0.071/-0.077 (68% CL). Curvature is constrained to ~0.7% in the owCDM model and to ~2% in a model in which dark energy is allowed to vary with parameters w_0 and w_a. Tightening further the constraints on the time evolution of dark energy will require several improvements, including high-quality multi-passband photometry of a sample of several dozen z>1 SNe Ia. We describe how such a sample could be efficiently obtained by targeting cluster fields with WFC3 on HST.

We study the environmental dependence of stellar population properties at z ~ 1.3. We derive galaxy properties (stellar masses, ages and star formation histories) for samples of massive, red, passive early-type galaxies in two high-redshift clusters, RXJ0849+4452 and RXJ0848+4453 (with redshifts of z = 1.26 and 1.27, respectively), and compare them with those measured for the RDCS1252.9-2927 cluster at z=1.24 and with those measured for a similarly mass-selected sample of field contemporaries drawn from the GOODS-South Field. Robust estimates of the aforementioned parameters have been obtained by comparing a large grid of composite stellar population models with extensive 8-10 band photometric coverage, from the rest-frame far-ultraviolet to the infrared. We find no variations of the overall stellar population properties among the different samples of cluster early-type galaxies. However, when comparing cluster versus field stellar population properties we find that, even if the (star formation weighted) ages are similar and depend only on galaxy mass, the ones in the field do employ longer timescales to assemble their final mass. We find that, approximately 1 Gyr after the onset of star formation, the majority (75%) of cluster galaxies have already assembled most (> 80%) of their final mass, while, by the same time, fewer (35%) field ETGs have. Thus we conclude that while galaxy mass regulates the timing of galaxy formation, the environment regulates the timescale of their star formation histories.

We have derived masses and ages for 79 early-type galaxies (ETGs) in different environments at z~1.3 in the Lynx supercluster and in the GOODS/CDF-S field using multiwavelength (0.6-4.5 $\mu$m; KPNO, Palomar, Keck, HST, Spitzer) datasets. At this redshift the contribution of the TP-AGB phase is important for ETGs, and the mass and age estimates depend on the choice of the stellar population model used in the spectral energy distribution fits. We describe in detail the differences among model predictions for a large range of galaxy ages, showing the dependence of these differences on age. Current models still yield large uncertainties. While recent models from Maraston and Charlot & Bruzual offer better modeling of the TP-AGB phase with respect to less recent Bruzual & Charlot models, their predictions do not often match. The modeling of this TP-AGB phase has a significant impact on the derived parameters for galaxies observed at high-redshift. Some of our results do not depend on the choice of the model: for all models, the most massive galaxies are the oldest ones, independent of the environment. When using Maraston and Charlot & Bruzual models, the mass distribution is similar in the clusters and in the groups, whereas in our field sample there is a deficit of massive (M $\gtrsim$ 10^11 Msun) ETGs. According to those last models, ETGs belonging to the cluster environment host on average older stars with respect to group and field populations. This difference is less significant than the age difference in galaxies of different masses.

We present a spectrophotometric analysis of galaxies belonging to the dynamically young, massive cluster RX J0152.7-1357 at z~0.84, aimed at understanding the effects of the cluster environment on the star formation history (SFH) of cluster galaxies and the assembly of the red-sequence (RS). We use VLT/FORS spectroscopy, ACS/WFC optical and NTT/SofI near-IR data to characterize SFHs as a function of color, luminosity, morphology, stellar mass, and local environment from a sample of 134 spectroscopic members. In order to increase the signal-to-noise, individual galaxy spectra are stacked according to these properties. Moreover, the D4000, Balmer, CN3883, Fe4383 and C4668 indices are also quantified. The SFH analysis shows that galaxies in the blue faint-end of the RS have on average younger stars (Delta t ~ 2 Gyr) than those in the red bright-end. We also found, for a given luminosity range, differences in age (Delta t ~ 0.5 - 1.3 Gyr) as a function of color, indicating that the intrinsic scatter of the RS may be due to age variations. Passive galaxies in the blue faint-end of the RS are preferentially located in the low density areas of the cluster, likely being objects entering the RS from the "blue cloud". It is likely that the quenching of the star formation of these RS galaxies is due to interaction with the intracluster medium. Furthermore, the SFH of galaxies in the RS as a function of stellar mass reveals signatures of "downsizing" in the overall cluster.

We present a multi-wavelength study of galaxy populations in the core of the massive, X-ray luminous cluster XMMU J2235 at z=1.39, based on VLT and HST optical and near-infrared photometry. Luminosity functions in the z, H, and Ks bands show a faint-end slope consistent with being flat, and a characteristic magnitude M* close to passive evolution predictions of M* of local massive clusters, with a formation redshift z>2. The color-magnitude and color-mass diagrams show evidence of a tight red sequence of massive galaxies, with overall old stellar populations, generally early-type morphology, typically showing early-type spectral features and rest-frame far-UV emission consistent with very low star formation rates (SFR<0.2Msun/yr). Star forming spectroscopic members, with SFRs of up to ~100Msun/yr, are all located at clustercentric distances >~250kpc, with the central cluster region already appearing effectively quenched. Massive galaxies in the core of this cluster appear to be in an advanced evolutionary stage in terms of both star formation and mass assembly. The high-mass end of the galaxy stellar mass function is essentially already in place, and the stellar mass fraction estimated within r500 (~1%, Kroupa IMF) is already similar to that of local massive clusters. On the other hand, morphological analysis of the massive red sequence galaxies suggests that they are smaller than similarly massive local early-types. While possibly affected by systematics and biases, this result might imply that, in spite of the overall early assembly of these sources, their evolution is not complete, and processes like minor (and likely dry) merging might still shape their structural properties to resemble those of their local counterparts, without substantially affecting their stellar mass or host stellar populations.[abridged]

The Advanced Technology Large-Aperture Space Telescope (ATLAST) is a concept for an 8-meter to 16-meter UVOIR space observatory for launch in the 2025-2030 era. ATLAST will allow astronomers to answer fundamental questions at the forefront of modern astronphysics, including "Is there life elsewhere in the Galaxy?" We present a range of science drivers that define the main performance requirements for ATLAST (8 to 16 milliarcsec angular resolution, diffraction limited imaging at 0.5 {\mu}m wavelength, minimum collecting area of 45 square meters, high sensitivity to light wavelengths from 0.1 {\mu}m to 2.4 {\mu}m, high stability in wavefront sensing and control). We will also discuss the synergy between ATLAST and other anticipated future facilities (e.g., TMT, EELT, ALMA) and the priorities for technology development that will enable the construction for a cost that is comparable to current generation observatory-class space missions.

Recommendations of the VAO-Science Council following the meeting of March 26-27, 2010. Refer to the webpage http://www.aui.edu/vao.php

The Advanced Technology Large Aperture Space Telescope (ATLAST) is a set of mission concepts for the next generation UV-Optical-Near Infrared space telescope with an aperture size of 8 to 16 meters. ATLAST, using an internal coronagraph or an external occulter, can characterize the atmosphere and surface of an Earth-sized exoplanet in the Habitable Zone of long-lived stars at distances up to ~45 pc, including its rotation rate, climate, and habitability. ATLAST will also allow us to glean information on the nature of the dominant surface features, changes in cloud cover and climate, and, potentially, seasonal variations in surface vegetation. ATLAST will be able to visit up to 200 stars in 5 years, at least three times each, depending on the technique used for starlight suppression and the telescope aperture. More frequent visits can be made for interesting systems.

Over 300 extrasolar planets (exoplanets) have been detected orbiting nearby stars. We now hope to conduct a census of all planets around nearby stars and to characterize their atmospheres and surfaces with spectroscopy. Rocky planets within their star's habitable zones have the highest priority, as these have the potential to harbor life. Our science goal is to find and characterize all nearby exoplanets; this requires that we measure the mass, orbit, and spectroscopic signature of each one at visible and infrared wavelengths. The techniques for doing this are at hand today. Within the decade we could answer long-standing questions about the evolution and nature of other planetary systems, and we could search for clues as to whether life exists elsewhere in our galactic neighborhood.

Ground-based adaptive optics (AO) in the infrared has made exceptional advances in approaching space-like image quality at higher collecting area. Optical-wavelength applications are now also growing in scope. We therefore provide here a comparison of the pros and cons of observational capabilities from the ground and from space at optical wavelengths. With an eye towards the future, we focus on the comparison of a ~30m ground-based telescope with an 8-16m space-based telescope. We review the current state-of-the-art in AO, and summarize the expected future improvements in image quality, field of view, contrast, and low-wavelength cut-off. We discuss the exciting advances in extreme AO for exoplanet studies and explore what the theoretical limitations in achievable contrast might be. Our analysis shows that extreme AO techniques face both fundamental and technological hurdles to reach the contrast of 1E-10 necessary to study an Earth-twin at 10 pc. Based on our assessment of the current state-of-the-art, the future technology developments, and the inherent difficulty of observing through a turbulent atmosphere, we conclude that there will continue to be a strong complementarity between observations from the ground and from space at optical wavelengths in the coming decades. There will continue to be subjects that can only be studied from space, including imaging and (medium-resolution) spectroscopy at the deepest magnitudes, and the exceptional-contrast observations needed to characterize terrestrial exoplanets and search for biomarkers.

We present a new survey strategy to discover and study high redshift Type Ia supernovae (SNe Ia) using the Hubble Space Telescope (HST). By targeting massive galaxy clusters at 0.9<z<1.5, we obtain a twofold improvement in the efficiency of finding SNe compared to an HST field survey and a factor of three improvement in the total yield of SN detections in relatively dust-free red-sequence galaxies. In total, sixteen SNe were discovered at z>0.95, nine of which were in galaxy clusters. This strategy provides a SN sample that can be used to decouple the effects of host galaxy extinction and intrinsic color in high redshift SNe, thereby reducing one of the largest systematic uncertainties in SN cosmology.

The Advanced Technology Large-Aperture Space Telescope (ATLAST) is a set of mission concepts for the next generation of UVOIR space observatory with a primary aperture diameter in the 8-m to 16-m range that will allow us to perform some of the most challenging observations to answer some of our most compelling questions, including "Is there life elsewhere in the Galaxy?" We have identified two different telescope architectures, but with similar optical designs, that span the range in viable technologies. The architectures are a telescope with a monolithic primary mirror and two variations of a telescope with a large segmented primary mirror. This approach provides us with several pathways to realizing the mission, which will be narrowed to one as our technology development progresses. The concepts invoke heritage from HST and JWST design, but also take significant departures from these designs to minimize complexity, mass, or both. Our report provides details on the mission concepts, shows the extraordinary scientific progress they would enable, and describes the most important technology development items. These are the mirrors, the detectors, and the high-contrast imaging technologies, whether internal to the observatory, or using an external occulter. Experience with JWST has shown that determined competitors, motivated by the development contracts and flight opportunities of the new observatory, are capable of achieving huge advances in technical and operational performance while keeping construction costs on the same scale as prior great observatories.

We have compiled a sample of early-type cluster galaxies from 0 < z < 1.3 and measured the evolution of their ellipticity distributions. Our sample contains 487 galaxies in 17 z>0.3 clusters with high quality space-based imaging and a comparable sample of 210 galaxies in 10 clusters at z<0.05. We select early-type galaxies (elliptical and S0 galaxies) that fall within the cluster R_{200}, and which lie on the red-sequence in the magnitude range -19.3 > M_B > -21, after correcting for luminosity evolution. Our ellipticity measurements are made in a consistent manner over our whole sample. We perform extensive simulations to quantify the systematic and statistical errors, and find that it is crucial to use PSF-corrected model fits. We find that neither the median ellipticity, nor the shape of the ellipticity distribution of cluster early-type galaxies evolves with redshift from z ~ 0 to z > 1. These results are strongly suggestive of an unchanging overall bulge-to-disk ratio distribution for cluster early-type galaxies over the last ~8Gyr. This result contrasts with that from visual classifications which show that the fraction of morphologically-selected disk-dominated early-type galaxies, or S0s, is significantly lower at z>0.4 than at z~0. Taking the ellipticity measurements and assuming, as in all previous studies, that the intrinsic ellipticity distribution of both elliptical and S0 galaxies remains constant, then we conclude from the lack of evolution in the observed early-type ellipticity distribution that the relative fractions of ellipticals and S0s do not evolve from z~1 to z=0 for a red-sequence selected samples of galaxies in the cores of clusters of galaxies.

If we are to develop a comprehensive and predictive theory of galaxy formation and evolution, it is essential that we obtain an accurate assessment of how and when galaxies assemble their stellar populations, and how this assembly varies with environment. There is strong observational support for the hierarchical assembly of galaxies, but by definition the dwarf galaxies we see today are not the same as the dwarf galaxies and proto-galaxies that were disrupted during the assembly. Our only insight into those disrupted building blocks comes from sifting through the resolved field populations of the surviving giant galaxies to reconstruct the star formation history, chemical evolution, and kinematics of their various structures. To obtain the detailed distribution of stellar ages and metallicities over the entire life of a galaxy, one needs multi-band photometry reaching solar-luminosity main sequence stars. The Hubble Space Telescope can obtain such data in the outskirts of Local Group galaxies. To perform these essential studies for a fair sample of the Local Universe will require observational capabilities that allow us to extend the study of resolved stellar populations to much larger galaxy samples that span the full range of galaxy morphologies, while also enabling the study of the more crowded regions of relatively nearby galaxies. With such capabilities in hand, we will reveal the detailed history of star formation and chemical evolution in the universe.

We apply detailed observations of the Color-Magnitude Relation (CMR) with the ACS/HST to study galaxy evolution in eight clusters at z~1. The early-type red sequence is well defined and elliptical and lenticular galaxies lie on similar CMRs. We analyze CMR parameters as a function of redshift, galaxy properties and cluster mass. For bright galaxies (M_B < -21mag), the CMR scatter of the elliptical population in cluster cores is smaller than that of the S0 population, although the two become similar at faint magnitudes. While the bright S0 population consistently shows larger scatter than the ellipticals, the scatter of the latter increases in the peripheral cluster regions. If we interpret these results as due to age differences, bright elliptical galaxies in cluster cores are on average older than S0 galaxies and peripheral elliptical galaxies (by about 0.5Gyr). CMR zero point, slope, and scatter in the (U-B)_z=0 rest-frame show no significant evolution out to redshift z~1.3 nor significant dependence on cluster mass. Two of our clusters display CMR zero points that are redder (by ~2sigma) than the average (U-B)_z=0 of our sample. We also analyze the fraction of morphological early-type and late-type galaxies on the red sequence. We find that, while in the majority of the clusters most (80% to 90%) of the CMR population is composed of early-type galaxies, in the highest redshift, low mass cluster of our sample, the CMR late-type/early-type fractions are similar (~50%), with most of the late-type population composed of galaxies classified as S0/a. This trend is not correlated with the cluster's X-ray luminosity, nor with its velocity dispersion, and could be a real evolution with redshift.

We present new constraints on the evolution of the early-type galaxy color-magnitude relation (CMR) based on deep near-infrared imaging of a galaxy protocluster at z=2.16 obtained using NICMOS on-board the Hubble Space Telescope. This field contains a spectroscopically confirmed space-overdensity of Lyman-alpha and H-alpha emitting galaxies which surrounds the powerful radio galaxy MRC 1138-262. Using these NICMOS data we identify a significant surface-overdensity (= 6.2x) of red J-H galaxies in the color-magnitude diagram (when compared with deep NICMOS imaging from the HDF-N and UDF). The optical-NIR colors of these prospective red-sequence galaxies indicate the presence of on-going dust-obscured star-formation or recently formed (<~ 1.5 Gyr)stellar populations in a majority of the red galaxies. We measure the slope and intrinsic scatter of the CMR for three different red galaxy samples selected by a wide color cut, and using photometric redshifts both with and without restrictions on rest-frame optical morphology. In all three cases both the rest-frame $U-B$ slope and intrinsic color scatter are considerably higher than corresponding values for lower redshift galaxy clusters. These results suggest that while some relatively quiescent galaxies do exist in this protocluster both the majority of the galaxy population and hence the color-magnitude relation are still in the process of forming, as expected.

The Sloan Digital Sky Survey (SDSS) and photometric/spectroscopic surveys in the GOODS-South field (the Chandra Deep Field-South, CDFS) are used to construct volume-limited, stellar mass-selected samples of galaxies at redshifts 0<z<1. The CDFS sample at 0.6<z<1.0 contains 207 galaxies complete down to M=4x10^10 Msol (for a ``diet'' Salpeter IMF), corresponding to a luminosity limit for red galaxies of M_B=-20.1. The SDSS sample at 0.020<z<0.045 contains 2003 galaxies down to the same mass limit, which corresponds to M_B=-19.3 for red galaxies. Morphologies are determined with an automated method, using the Sersic parameter n and a measure of the residual from the model fits, called ``bumpiness'', to distinguish different morphologies. These classifications are verified with visual classifications. In agreement with previous studies, 65-70% of the galaxies are located on the red sequence, both at z~0.03 and at z~0.8. Similarly, 65-70% of the galaxies have n>2.5. The fraction of E+S0 galaxies is 43+/-3%$ at z~0.03 and 48+/-7% at z~0.8, i.e., it has not changed significantly since z~0.8. When combined with recent results for cluster galaxies in the same redshift range, we find that the morphology-density relation for galaxies more massive than 0.5M* has remained constant since at least z~0.8. This implies that galaxies evolve in mass, morphology and density such that the morphology-density relation does not change. In particular, the decline of star formation activity and the accompanying increase in the stellar mass density of red galaxies since z~1 must happen without large changes in the early-type galaxy fraction in a given environment.

We examined the morphology-density relations for galaxy samples selected by luminosity and by mass in each of five massive X-ray clusters from z=0.023 to 0.83 for 674 spectroscopically-confirmed members. Rest-frame optical colors and visual morphologies were obtained primarily from Hubble Space Telescope images. Morphology-density relations (MDR) are derived in each cluster from a complete, luminosity-selected sample of 452 galaxies with a magnitude limit M_V < M^{*}_{V} + 1. The change in the early-type fraction with redshift matches previous work for massive clusters of galaxies. We performed a similar analysis, deriving MDRs for complete, mass-selected samples of 441 galaxies with a mass-limit of 10^{10.6} M_{\sun}. Our mass limit includes faint objects, the equivalent of =~1 mag below L^{*} for the red cluster galaxies, and encompasses =~70% of the stellar mass in cluster galaxies. The MDRs in the mass-selected sample at densities of Sigma > 50 galaxies Mpc^{-2} are similar to those in the luminosity-selected sample but show larger early-type fractions. However, the trend with redshift in the fraction of elliptical and S0 galaxies with masses > 10^{10.6} M_{\sun} differs significantly between the mass- and luminosity-selected samples. The clear trend seen in the early-type fraction from z=0 to z=~ 0.8 is not found in mass-selected samples. The early-type galaxy fraction changes much less, and is consistent with being constant at 92% +/- 4% at \Sigma> 500 galaxies Mpc^{-2} and 83 +/- 3% at 50 < \Sigma < 500 galaxies Mpc^{-2}. This suggests that galaxies of mass lower than > 10^{10.6} M_{\sun} play a significant role in the evolution of the early-type fraction in luminosity-selected samples. (Abstract abridged)

We present results from an extensive spectroscopic survey, carried out with VLT FORS, and from an extensive multiwavelength imaging data set from the HST Advanced Camera for Surveys and ground based facilities, of the cluster of galaxies RDCS J1252.9-2927. We have spectroscopically confirmed 38 cluster members in the redshift range 1.22 < z < 1.25. A cluster median redshift of z=1.237 and a rest-frame velocity dispersion of 747^{+74}_{-84} km/s are obtained. Using the 38 confirmed redshifts, we were able to resolve, for the first time at z > 1, kinematic structure. The velocity distribution, which is not Gaussian at the 95% confidence level, is consistent with two groups that are also responsible for the projected east-west elongation of the cluster. The groups are composed of 26 and 12 galaxies with velocity dispersions of 486^{+47}_{-85} km/s and 426^{+57}_{-105} km/s, respectively. The elongation is also seen in the intracluster gas and the dark matter distribution. This leads us to conclude that RDCS J1252.9-2927 has not yet reached a final virial state. We extend the analysis of the color-magnitude diagram of spectroscopic members to more than 1 Mpc from the cluster center. The scatter and slope of non-[OII]-emitting cluster members in the near-IR red sequence is similar to that seen in clusters at lower redshift. Furthermore, most of the galaxies with luminosities greater than ~ K_s*+1.5 do not show any [OII], indicating that these more luminous, redder galaxies have stopped forming stars earlier than the fainter, bluer galaxies. Our observations provide detailed dynamical and spectrophotometric information on galaxies in this exceptional high-redshift cluster, delivering an in-depth view of structure formation at this epoch only 5 Gyr after the Big Bang.

Using a magnitude-limited, spectroscopic survey of the X-ray luminous galaxy cluster MS1054-03, we isolate 153 cluster galaxies and measure MS1054's redshift and velocity dispersion to be z=0.8307 and sigma_z=1156 km/s. The absorption-line, post-starburst (``E+A''), and emission-line galaxies respectively make up 63%, 15%, and 23% of the cluster population. With photometry from HST/ACS, we find that the absorption-line members define an exceptionally tight red sequence over a span of ~3.5 magnitudes in i775: their intrinsic scatter in (V-i) color is only 0.048, corresponding to a (U-B)_z scatter of 0.041. Their color scatter is comparable to that of the ellipticals (sigma_Vi=0.055), but measurably smaller than that of the combined E+S0 sample (sigma_Vi=0.072). The color scatter of MS1054's absorption-line population is approximately twice that of the ellipticals in Coma; this difference is consistent with passive evolution where most of the absorption-line members (>75%) formed by z~2, and all of them by z~1.2. For red members, we find a trend (>95% confidence) of weakening Hdelta absorption with redder colors that we conclude is due to age: in MS1054, the color scatter on the red sequence is driven by differences in mean stellar age of up to ~1.5 Gyr. We also generate composite spectra and estimate that the average S0 in MS1054 is ~0.5-1 Gyr younger than the average elliptical; this difference in mean stellar age is mainly due to a number of S0s that are blue (18%) and/or are post-starburst systems (21%).

Black hole masses predicted from the Mbh-sigma relationship conflict with those predicted from the Mbh-L relationship for the most luminous galaxies, such as brightest cluster galaxies (BCGs). This is because stellar velocity dispersion, sigma, increases only weakly with L for BCGs and other giant ellipticals. The Mbh-L relationship predicts that the most luminous BCGs may have Mbh approaching 10^{10}M_sol, while the M-sigma relationship always predicts Mbh<3X10^9M_sol. We argue that the Mbh-L relationship is a plausible or even preferred description for BCGs and other galaxies of similar luminosity. If cores in central stellar density are formed by binary BHs, then the inner-core cusp radius, r_gamma, may be an independent witness of Mbh. Using structural parameters derived from a large sample of early-type galaxies observed by HST, we argue that L is superior to sigma as an indicator of r_gamma in luminous galaxies. The observed r_gamma-Mbh relationship for 11 core galaxies with measured Mbh appears to be consistent with the Mbh-L relationship for BCGs. BCGs have large cores appropriate for their large luminosities that may be difficult to generate with the modest BH masses inferred from the Mbh-sigma relationship. Mbh~L may be expected to hold for BCGs, if they were formed in dissipationless mergers, which should preserve ratio of BH to stellar mass. This picture appears to be consistent with the slow increase in sigma with L and the more rapid increase in effective radii with L seen in BCGs. If BCGs have large BHs commensurate with their luminosities, then the local BH mass function for Mbh>3X10^9M_sol may be nearly an order of magnitude richer than that inferred from the Mbh-sigma relationship. The volume density of QSOs at earlier epochs may favor the predictions from the Mbh-L relationship.

We present a deep image of the radio galaxy MRC 1138-262 taken with the Hubble Space Telescope (HST) at a redshift of z = 2.2. The galaxy is known to have properties of a cD galaxy progenitor and be surrounded by a 3 Mpc-sized structure, identified with a protocluster. The morphology shown on the new deep HST/ACS image is reminiscent of a spider's web. More than 10 individual clumpy features are observed, apparently star-forming satellite galaxies in the process of merging with the progenitor of a dominant cluster galaxy 11 Gyr ago. There is an extended emission component, implying that star formation was occurring over a 50 times 40 kpc region at a rate of more than 100 M_sun/yr. A striking feature of the newly named ``Spiderweb galaxy'' is the presence of several faint linear galaxies within the merging structure. The dense environments and fast galaxy motions at the centres of protoclusters may stimulate the formation of these structures, which dominate the faint resolved galaxy populations in the Hubble Ultra Deep Field. The new image provides a unique testbed for simulations of forming dominant cluster galaxies.

We analyse the cluster color-magnitude relation (CMR) for early-type galaxies in two of the richer clusters in the z ~ 0.9 supercluster system to derive average ages and formation redshifts for the early-type galaxy population. Both clusters were observed with the Advanced Camera for Surveys aboard the {\it Hubble Space Telescope} through the F606W and F814W filters, which brackets the rest-frame 4000 \AA break at the cluster redshifts of z\sim 0.9. We fit the zeropoint and slope of the red cluster sequence, and model the scatter about this relation to estimate average galaxy ages and formation redshifts. We find intrinsic scatters of 0.038-0.053 mag in ($V_{606}-I_{814}$) for the E and E+S0 populations, corresponding to average ages of 3.5-3.7 Gyr and formation redshifts z_{f}=2.4-2.6. We find at least one significant difference between the Cl1604+4304 and Cl1604+4321 early-type CMRs. Cl1604+4321, the less X-ray luminous and massive of the two, lacks bright L^* ellipticals. We combine the galaxy samples to fit a composite CMR down to 0.15L^*, and find that the slope of the combined cluster CMR is significantly steeper than for RX J0152.7-1357 but consistent with MS 1054-03, both at similar redshift. The slope of the Cl1604 CMR at the bright end (L > 0.5L^*) is flatter and consistent with the CMR slopes found for other high redshift clusters. We find evidence for increasing scatter with increasing magnitude along the early-type CMR, consistent with a 'downsizing' scenario, indicating younger mean ages with decreasing galaxy mass.

We examine the distribution of stellar masses of galaxies in MS 1054-03 and RX J0152.7-1357, two X-ray selected clusters of galaxies at z=0.83. Our stellar mass estimates, from spectral energy distribution fitting, reproduce the dynamical masses as measured from velocity dispersions and half-light radii with a scatter of 0.2 dex in the mass for early-type galaxies. When we restrict our sample of members to high stellar masses, > 1e11.1 Msun (M* in the Schechter mass function for cluster galaxies), we find that the fraction of early-type galaxies is 79 +/- 6% at z=0.83 and 87 +/- 6% at z=0.023 for the Coma cluster, consistent with no evolution. Previous work with luminosity-selected samples finds that the early-type fraction in rich clusters declines from =~80% at z=0 to =~60% at z=0.8. The observed evolution in the early-type fraction from luminosity-selected samples must predominately occur among sub-M* galaxies. As M* for field and group galaxies, especially late-types, is below M* for clusters galaxies, infall could explain most of the recent early-type fraction growth. Future surveys could determine the morphological distributions of lower mass systems which will confirm or refute this explanation.

Color-magnitude relations have been derived in the clusters RX J0849+4452 and RX J0848+4453 at z~1.26. The color-magnitude relation was determined from Advanced Camera for Surveys imaging in the WFC F775W (i_775) and F850LP (z_850) filters combined with ground-based spectroscopy. Early-type cluster candidates have been identified according to the Postman et al. morphological classification. In both clusters the bright red early-type population defines a tight color-magnitude relation very similar in color, although the two clusters present different X-ray luminosities and shapes, with RX J0849+4452 being three times more X-ray luminous and more compact, and having a temperature two times higher. The elliptical galaxy color-magnitude relations (CMR) in RX J0849+4452 and RX J0848+4453 show an intrinsic (i_775-z_850) color scatter of 0.026 +/- 0.012 mag and 0.024 +/- 0.015 mag, respectively, within 2 arcminutes (~1Mpc at z=1.26) from the cluster X-ray emission centers. Simple modeling of the scatters using stellar population models from Bruzual and Charlot, gives a mean luminosity-weighted age t > 2.5 Gyr (z_f > 2.75). S0 galaxies follow the elliptical CMR; they show larger scatters about the CMR. The intrinsic scatter decreases and the CMR slopes are steeper at smaller radii, within both clusters. We conclude that old stellar populations in cluster elliptical galaxies are already in place at z=1.26, both in the more evolved cluster RX J0849+4452, and in its less evolved companion RX J0848+4453. Even at a lookback time of 9 Gyr, in the early merging and buildup of massive clusters, the bulk of the stellar content of the bright elliptical galaxy population was in place - apparently formed some 2.5~Gyr earlier at z~3

We study the photometric and structural properties of spectroscopically confirmed members in the two massive X-ray--selected z=0.83 galaxy clusters MS1054-03 and RXJ0152-1357 using three-band mosaic imaging with the Hubble Space Telescope Advanced Camera for Surveys. The samples include 105 and 140 members of MS1054-03 and RXJ0152-1357, respectively, with ACS F775W magnitude < 24.0. We develop a promising new structural classification method, based on a combination of the best-fit Sersic indices and the normalized root-mean-square residuals from the fits; the resulting classes agree well with the visual ones, but are less affected by galaxy orientation. We examine the color--magnitude relations in detail and find that the color residuals correlate with the local mass density measured from our weak lensing maps; we identify a threshold density of $\Sigma \approx 0.1$, in units of the critical density, above which the star formation appears to cease. For RXJ0152-1357, we also find a trend in the color residuals with velocity, resulting from an offset of about 980 km/s in the mean redshifts of the early- and late-type galaxies. Analysis of the color--color diagrams indicates that a range of star formation time-scales are needed to reproduce the loci of the galaxy colors. We also identify some cluster galaxies whose colors can only be explained by large amounts, $A_V \approx 1$ mag, of internal dust extinction. [Abstract shortened]

(Abridged.) We present F435W (B), F606W (Broad V), and F814W (Broad I) coronagraphic images of the debris disk around Beta Pictoris obtained with HST's Advanced Camera for Surveys. We confirm that the previously reported warp in the inner disk is a distinct secondary disk inclined by ~5 deg from the main disk. The main disk's northeast extension is linear from 80 to 250 AU, but the southwest extension is distinctly bowed with an amplitude of ~1 AU over the same region. Both extensions of the secondary disk appear linear, but not collinear, from 80 to 150 AU. Within ~120 AU of the star, the main disk is ~50% thinner than previously reported. The surface-brightness profiles along the spine of the main disk are fitted with four distinct radial power laws between 40 and 250 AU, while those of the secondary disk between 80 and 150 AU are fitted with single power laws. These discrepancies suggest that the two disks have different grain compositions or size distributions. The F606W/F435W and F814W/F435W flux ratios of the composite disk are nonuniform and asymmetric about both projected axes of the disk. Within ~120 AU, the m_F435W-m_F606W and m_F435W-m_F814W colors along the spine of the main disk are ~10% and ~20% redder, respectively, than those of Beta Pic. These colors increasingly redden beyond ~120 AU, becoming 25% and 40% redder, respectively, than the star at 250 AU. We compare the observed red colors within ~120 AU with the simulated colors of non-icy grains having a radial number density ~r^-3 and different compositions, porosities, and minimum grain sizes. The observed colors are consistent with those of compact or moderately porous grains of astronomical silicate and/or graphite with sizes >0.15-0.20 um, but the colors are inconsistent with the blue colors expected from grains with porosities >90%.

The color-magnitude relation has been determined for the RDCS J0910+5422 cluster of galaxies at redshift z = 1.106. Cluster members were selected from HST ACS images, combined with ground--based near--IR imaging and optical spectroscopy. The observed early--type color--magnitude relation (CMR) in (i_775 -z_850) versus z_850 shows intrinsic scatters in color of 0.042 +/- 0.010 mag and 0.044 +/- 0.020 mag for ellipticals and S0s, respectively. From the scatter about the CMR, a mean luminosity--weighted age t > 3.3 Gyr (z > 3) is derived for the elliptical galaxies. Strikingly, the S0 galaxies in RDCS J0910+5422 are systematically bluer in (i_775 - z_850) by 0.07 +/- 0.02 mag, with respect to the ellipticals. The ellipticity distribution as a function of color indicates that the face-on S0s in this particular cluster have likely been classified as elliptical. Thus, if anything, the offset in color between the elliptical and S0 populations may be even more significant. The color offset between S0 and E corresponds to an age difference of ~1 Gyr, for a single-burst solar metallicity model. A solar metallicity model with an exponential decay in star formation will reproduce the offset for an age of 3.5 Gyr, i.e. the S0s have evolved gradually from star forming progenitors. The early--type population in this cluster appears to be still forming. The blue early-type disk galaxies in RDCS J0910+5422 likely represent the direct progenitors of the more evolved S0s that follow the same red sequence as ellipticals in other clusters. Thirteen red galaxy pairs are observed and the galaxies associated in pairs constitute ~40% of the CMR galaxies in this cluster.

A five square arcminute region around the luminous radio-loud quasar SDSS J0836+0054 (z=5.8) hosts a wealth of associated galaxies, characterized by very red (1.3 < i_775 - z_{850} < 2.0) color. The surface density of these z~5.8 candidates is approximately six times higher than the number expected from deep ACS fields. This is one of the highest galaxy overdensities at high redshifts, which may develop into a group or cluster. We also find evidence for a substructure associated with one of the candidates. It has two very faint companion objects within two arcseconds, which are likely to merge. The finding supports the results of a recent simulation that luminous quasars at high redshifts lie on the most prominent dark-matter filaments and are surrounded by many fainter galaxies. The quasar activity from these regions may signal the buildup of a massive system.

We present a comparison of late-type galaxies (Sa and later) in intermediate redshift clusters and the field using ACS imaging of four cluster fields: CL0152-1357, CL1056-0337 (MS1054), CL1604+4304, and CL1604+4321. Concentration, asymmetry, and clumpiness parameters are calculated for each galaxy in blue (F606W or F625W) and red (F775W or F814W) filters. Galaxy half-light radii, disk scale lengths, color gradients, and overall color are compared. We find marginally significant differences in the asymmetry distributions of spiral and irregular galaxies in the X-ray luminous and X-ray faint clusters. The massive clusters contain fewer galaxies with large asymmetries. The physical sizes of the cluster and field populations are similar; no significant differences are found in half-light radii or disk scale lengths. The most significant difference is in rest-frame $U-B$ color. Late-type cluster galaxies are significantly redder, $\sim 0.3$ magnitudes at rest-frame $U-B$, than their field counterparts. Moreover, the intermediate-redshift cluster galaxies tend to have blue inward color gradients, in contrast to the field galaxies, but similar to late-type galaxies in low redshift clusters. These blue inward color gradients are likely to be the result of enhanced nuclear star formation rates relative to the outer disk. Based on the significant rest-frame color difference, we conclude that late-type cluster members at $z\sim0.9$ are not a pristine infalling field population; some difference in past and/or current star formation history is already present. This points to high redshift ``groups'', or filaments with densities similar to present-day groups, as the sites where the first major effects of environment are imprinted.

We present the "C4 Cluster Catalog", a new sample of 748 clusters of galaxies identified in the spectroscopic sample of the Second Data Release (DR2) of the Sloan Digital Sky Survey (SDSS). The C4 cluster--finding algorithm identifies clusters as overdensities in a seven-dimensional position and color space, thus minimizing projection effects which plagued previous optical clusters selection. The present C4 catalog covers ~2600 square degrees of sky with groups containing 10 members to massive clusters having over 200 cluster members with redshifts. We provide cluster properties like sky location, mean redshift, galaxy membership, summed r--band optical luminosity (L_r), velocity dispersion, and measures of substructure. We use new mock galaxy catalogs to investigate the sensitivity to the various algorithm parameters, as well as to quantify purity and completeness. These mock catalogs indicate that the C4 catalog is ~90% complete and 95% pure above M_200 = 1x10^14 solar masses and within 0.03 <=z <= 0.12. The C4 algorithm finds 98% of X-ray identified clusters and 90% of Abell clusters within 0.03 <= z <= 0.12. We show that the L_r of a cluster is a more robust estimator of the halo mass (M_200) than the line-of-sight velocity dispersion or the richness of the cluster. L_r. The final SDSS data will provide ~2500 C4 clusters and will represent one of the largest and most homogeneous samples of local clusters.

We study the dynamical origin of the structures observed in the scattered-light images of the resolved debris disk around HD 141569A. We explore the roles of radiation pressure from the central star, gas drag from the gas disk, and the tidal forces from two nearby stars in creating and maintaining these structures. We use a simple one-dimensional axisymmetric model to show that the presence of the gas helps confine the dust and that a broad ring of dust is produced if a central hole exists in the disk. This model also suggests that the disk is in a transient, excited dynamical state, as the observed dust creation rate applied over the age of the star is inconsistent with submillimeter mass measurements. We model in two dimensions the effects of a fly-by encounter between the disk and a binary star in a prograde, parabolic, coplanar orbit. We track the spatial distribution of the disk's gas, planetesimals, and dust. We conclude that the surface density distribution reflects the planetesimal distribution for a wide range of parameters. Our most viable model features a disk of initial radius 400 AU, a gas mass of 50 M_earth, and beta = 4 and suggests that the system is being observed within 4000 yr of the fly-by periastron. The model reproduces some features of HD 141569A's disk, such as a broad single ring and large spiral arms, but it does not reproduce the observed multiple spiral rings or disk asymmetries nor the observed clearing in the inner disk. For the latter, we consider the effect of a 5 M_Jup planet in an eccentric orbit on the planetesimal distribution of HD 141569A.

We investigate the evolution in the distribution of surface brightness, as a function of size, for elliptical and S0 galaxies in the two clusters RDCS J1252.9-2927, z=1.237 and RX J0152.7-1357, z=0.837. We use multi-color imaging with the Advanced Camera for Surveys on the Hubble Space Telescope to determine these sizes and surface brightnesses. Using three different estimates of the surface brightnesses, we find that we reliably estimate the surface brightness for the galaxies in our sample with a scatter of < 0.2 mag and with systematic shifts of \lesssim 0.05 mag. We construct samples of galaxies with early-type morphologies in both clusters. For each cluster, we use a magnitude limit in a band which closely corresponds to the rest-frame B, to magnitude limit of M_B = -18.8 at z=0, and select only those galaxies within the color-magnitude sequence of the cluster or by using our spectroscopic redshifts. We measure evolution in the rest-frame B surface brightness, and find -1.41 \+/- 0.14 mag from the Coma cluster of galaxies for RDCS J1252.9-2927 and -0.90 \+/- 0.12 mag of evolution for RX J0152.7-1357, or an average evolution of (-1.13 \+/- 0.15) z mag. Our statistical errors are dominated by the observed scatter in the size-surface brightness relation, sigma = 0.42 \+/- 0.05 mag for RX J0152.7-1357 and sigma = 0.76 \+/- 0.10 mag for RDCS J1252.9-2927. We find no statistically significant evolution in this scatter, though an increase in the scatter could be expected. Overall, the pace of luminosity evolution we measure agrees with that of the Fundamental Plane of early-type galaxies, implying that the majority of massive early-type galaxies observed at z =~ 1 formed at high redshifts.

We measure the morphology--density relation (MDR) and morphology-radius relation (MRR) for galaxies in seven z ~ 1 clusters that have been observed with the Advanced Camera for Surveys on board the Hubble Space Telescope. Simulations and independent comparisons of ourvisually derived morphologies indicate that ACS allows one to distinguish between E, S0, and spiral morphologies down to zmag = 24, corresponding to L/L* = 0.21 and 0.30 at z = 0.83 and z = 1.24, respectively. We adopt density and radius estimation methods that match those used at lower redshift in order to study the evolution of the MDR and MRR. We detect a change in the MDR between 0.8 < z < 1.2 and that observed at z ~ 0, consistent with recent work -- specifically, the growth in the bulge-dominated galaxy fraction, f_E+SO, with increasing density proceeds less rapidly at z ~ 1 than it does at z ~ 0. At z ~ 1 and density <= 500 galaxies/Mpc^2, we find <f_E+S0> = 0.72 +/- 0.10. At z ~ 0, an E+S0 population fraction of this magnitude occurs at densities about 5 times smaller. The evolution in the MDR is confined to densities >= 40 galaxies/Mpc^2 and appears to be primarily due to a deficit of S0 galaxies and an excess of Spiral+Irr galaxies relative to the local galaxy population. The Elliptical fraction - density relation exhibits no significant evolution between z = 1 and z = 0. We find mild evidence to suggest that the MDR is dependent on the bolometric X-ray luminosity of the intracluster medium. Implications for the evolution of the disk galaxy population in dense regions are discussed in the context of these observations.

We present a weak lensing analysis of one of the most distant massive galaxy cluster known, RDCS 1252.9-2927 at z=1.24, using deep images from the Advanced Camera for Survey (ACS) on board the Hubble Space Telescope (HST). By taking advantage of the depth and of the angular resolution of the ACS images, we detect for the first time at z>1 a clear weak lensing signal in both the i (F775W) and z (F850LP) filters. We measure a 5-\sigma signal in the i band and a 3-\sigma signal in the shallower z band image. The two radial mass profiles are found to be in very good agreement with each other, and provide a measurement of the total mass of the cluster inside a 1Mpc radius of M(<1Mpc) = (8.0 +/- 1.3) x 10^14 M_\odot in the current cosmological concordance model h =0.70, \Omega_m=0.3, \Omega_\Lambda=0.7, assuming a redshift distribution of background galaxies as inferred from the Hubble Deep Fields surveys. A weak lensing signal is detected out to the boundary of our field (3' radius, corresponding to 1.5Mpc at the cluster redshift). We detect a small offset between the centroid of the weak lensing mass map and the brightest cluster galaxy, and we discuss the possible origin of this discrepancy. The cumulative weak lensing radial mass profile is found to be in good agreement with the X-ray mass estimate based on Chandr and XMM-Newton observations, at least out to R_500=0.5Mpc.

We present new measurements of the galaxy luminosity function (LF) and its dependence on local galaxy density, color, morphology, and clustocentric radius for the massive z=0.83 cluster MS1054-0321. Our analyses are based on imaging performed with the ACS onboard the HST in the F606W, F775W and F850LP passbands and extensive spectroscopic data obtained with the Keck LRIS. Our main results are based on a spectroscopically selected sample of 143 cluster members with morphological classifications derived from the ACS observations. Our three primary findings are (1) the faint-end slope of the LF is steepest in the bluest filter, (2) the LF in the inner part of the cluster (or highest density regions) has a flatter faint-end slope, and (3) the fraction of early-type galaxies is higher at the bright end of the LF, and gradually decreases toward fainter magnitudes. These characteristics are consistent with those in local galaxy clusters, indicating that, at least in massive clusters, the common characteristics of cluster LFs are established at z=0.83. We also find a 2sigma deficit of intrinsically faint, red galaxies (i-z>0.5, Mi>-19) in this cluster. This trend may suggest that faint, red galaxies (which are common in z<0.1 rich clusters) have not yet been created in this cluster at z=0.83. The giant-to-dwarf ratio in MS1054-0321 starts to increase inwards of the virial radius or when Sigma>30 Mpc^-2, coinciding with the environment where the galaxy star formation rate and the morphology-density relation start to appear. (abridged)

We combine imaging data from the Advanced Camera for Surveys (ACS) with VLT/FORS optical spectroscopy to study the properties of star-forming galaxies in the z=0.837 cluster CL0152-1357. We have morphological information for 24 star-forming cluster galaxies, which range in morphology from late-type and irregular to compact early-type galaxies. We find that while most star-forming galaxies have $r_{625}-i_{775}$ colors bluer than 1.0, eight are in the red cluster sequence. Among the star-forming cluster population we find five compact early-type galaxies which have properties consistent with their identification as progenitors of dwarf elliptical galaxies. The spatial distribution of the star-forming cluster members is nonuniform. We find none within $R\sim 500$ Mpc of the cluster center, which is highly suggestive of an intracluster medium interaction. We derive star formation rates from [OII] $\lambda\lambda 3727$ line fluxes, and use these to compare the global star formation rate of CL0152-1357 to other clusters at low and intermediate redshifts. We find a tentative correlation between integrated star formation rates and $T_{X}$, in the sense that hotter clusters have lower integrated star formation rates. Additional data from clusters with low X-ray temperatures is needed to confirm this trend. We do not find a significant correlation with redshift, suggesting that evolution is either weak or absent between z=0.2-0.8.

Combined HST, X-ray, and ground-based optical studies show that clusters of galaxies are largely "in place" by $z \sim 1$, an epoch when the Universe was less than half its present age. High resolution images show that elliptical, S0, and spiral galaxies are present in clusters at redshifts up to $z \sim 1.3$. Analysis of the CMDs suggest that the cluster ellipticals formed their stars several Gyr earlier, near redshift 3. The morphology--density relation is well established at $z\sim1$, with star-forming spirals and irregulars residing mostly in the outer parts of the clusters and E/S0s concentrated in dense clumps. The intracluster medium has already reached the metallicity of present-day clusters. The distributions of the hot gas and early-type galaxies are similar in $z\sim1$ clusters, indicating both have largely virialized in the deepest potentials wells. In spite of the many similarities between $z\sim1$ and present-day clusters, there are significant differences. The morphologies revealed by the hot gas, and particularly the early-type galaxies, are elongated rather than spherical. We appear to be observing the clusters at an epoch when the sub-clusters and groups are still assembling into a single regular cluster. Support for this picture comes from CL0152 where the gas appears to be lagging behind the luminous and dark mass in two merging sub-components. Moreover, the luminosity difference between the first and second brightest cluster galaxies at $z\sim1$ is smaller than in 93% of present-day Abell clusters, which suggests that considerable luminosity evolution through merging has occurred since that epoch. Evolution is also seen in the bolometric X-ray luminosity function.

We measure the luminosity function of morphologically selected E/S0 galaxies from $z=0.5$ to $z=1.0$ using deep high resolution Advanced Camera for Surveys imaging data. Our analysis covers an area of $48\Box\arcmin$ (8$\times$ the area of the HDF-N) and extends 2 magnitudes deeper ($I\sim24$ mag) than was possible in the Deep Groth Strip Survey (DGSS). At $0.5<z<0.75$, we find $M_B^*-5\log h_{0.7}=-21.1\pm0.3$ and $\alpha=-0.53\pm0.2$, and at $0.75<z<1.0$, we find $M_B^*-5\log h_{0.7}=-21.4\pm0.2$. These luminosity functions are similar in both shape and number density to the luminosity function using morphological selection (e.g., DGSS), but are much steeper than the luminosity functions of samples selected using morphological proxies like the color or spectral energy distribution (e.g., CFRS, CADIS, or COMBO-17). The difference is due to the `blue', $(U-V)_0<1.7$, E/S0 galaxies, which make up to $\sim30%$ of the sample at all magnitudes and an increasing proportion of faint galaxies. We thereby demonstrate the need for {\it both morphological and structural information} to constrain the evolution of galaxies. We find that the `blue' E/S0 galaxies have the same average sizes and Sersic parameters as the `red', $(U-V)_0>1.7$, E/S0 galaxies at brighter luminosities ($M_B<-20.1$), but are increasingly different at fainter magnitudes where `blue' galaxies are both smaller and have lower Sersic parameters. Fits of the colors to stellar population models suggest that most E/S0 galaxies have short star-formation time scales ($\tau<1$ Gyr), and that galaxies have formed at an increasing rate from $z\sim8$ until $z\sim2$ after which there has been a gradual decline.

The properties of Ultra Compact Dwarf (UCD) galaxy candidates in Abell 1689 (z=0.183) are investigated, based on deep high resolution ACS images. A UCD candidate has to be unresolved, have i<28 (M_V<-11.5) mag and satisfy color limits derived from Bayesian photometric redshifts. We find 160 UCD candidates with 22<i<28 mag. It is estimated that about 100 of these are cluster members, based on their spatial distribution and photometric redshifts. For i>26.8 mag, the radial and luminosity distribution of the UCD candidates can be explained well by Abell 1689's globular cluster (GC) system. For i<26.8 mag, there is an overpopulation of 15 +/- 5 UCD candidates with respect to the GC luminosity function. For i<26 mag, the radial distribution of UCD candidates is more consistent with the dwarf galaxy population than with the GC system of Abell 1689. The UCD candidates follow a color-magnitude trend with a slope similar to that of Abell 1689's genuine dwarf galaxy population, but shifted fainter by about 2-3 mag. Two of the three brightest UCD candidates (M_V ~ -17 mag) are slightly resolved. At the distance of Abell 1689, these two objects would have King-profile core radii of ~35 pc and r_eff ~300 pc, implying luminosities and sizes 2-3 times those of M32's bulge. Additional photometric redshifts obtained with late type stellar and elliptical galaxy templates support the assignment of these two resolved sources to Abell 1689. Our findings imply that in Abell 1689 there are at least 10 UCDs with M_V<-12.7 mag. Compared to the UCDs in the Fornax cluster they are brighter, larger and have colors closer to normal dwarf galaxies. This suggests that they may be in an intermediate stage of the stripping process. Spectroscopy is needed to definitely confirm the existence of UCDs in Abell 1689.

(Abridged) We study the internal color properties of a morphologically selected sample of spheroidal galaxies taken from HST/ACS ERO program of UGC 10214 (``The Tadpole''). By taking advantage of the unprecedented high resolution of the ACS in this very deep dataset we are able to characterize spheroids at sub-arcseconds scales. Using the V_606W and I_814W bands, we construct V-I color maps and extract color gradients for a sample of spheroids at I_814W < 24 mag. We investigate the existence of a population of morphologically classified spheroids which show extreme variation in their internal color properties similar to the ones reported in the HDFs. These are displayed as blue cores and inverse color gradients with respect to those accounted from metallicity variations. Following the same analysis we find a similar fraction of early-type systems (~30%-40%) that show non-homologous internal colors, suggestive of recent star formation activity. We present two statistics to quantify the internal color variation in galaxies and for tracing blue cores, from which we estimate the fraction of non-homogeneous to homogeneous internal colors as a function of redshift up to z<1.2. We find that it can be described as about constant as a function of redshift, with a small increase with redshift for the fraction of spheroids that present strong color dispersions. The implications of a constant fraction at all redshifts suggests the existence of a relatively permanent population of evolving spheroids up to z~1. We discuss the implications of this in the context of spheroidal formation.

We report the discovery of a strongly lensed arc system around a field elliptical galaxy in Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) images of a parallel field observed during NICMOS observations of the HST Ultra-Deep Field. The ACS parallel data comprise deep imaging in the F435W, F606W, F775W, and F850LP bandpasses. The main arc is at a radius of 1.6 arcsec from the galaxy center and subtends about 120 deg. Spectroscopic follow-up at Magellan Observatory yields a redshift z=0.6174 for the lensing galaxy, and we photometrically estimate z_phot = 2.4\pm0.3 for the arc. We also identify a likely counter-arc at a radius of 0.6 arcsec, which shows structure similar to that seen in the main arc. We model this system and find a good fit to an elliptical isothermal potential of velocity dispersion $\sigma \approx 300$ \kms, the value expected from the fundamental plane, and some external shear. Several other galaxies in the field have colors similar to the lensing galaxy and likely make up a small group.

We present the first detailed X-ray observations of optically-selected clusters at high redshift. Two clusters, Cl 1324+3011 at z = 0.76 and Cl 1604+4304 at z = 0.90, were observed with XMM-Newton. The optical center of each cluster is coincident with an extended X-ray source whose emission is detected out to a radius of 0.5 Mpc. The emission from each cluster appears reasonably circular, with some indication of asymmetries and more complex morphologies. Similarly to other optically-selected clusters at redshifts of z > 0.4, both clusters are modest X-ray emitters with bolometric luminosities of only Lx = 1.4 - 2.0 x 10^(44) erg/s. We measure gas temperatures of T = 2.88 (+0.71/-0.49) keV for Cl 1324+3011 and 2.51 (+1.05/-0.69) keV for Cl 1604+4304. The X-ray properties of both clusters are consistent with the high-redshift Lx-T relation measured from X-ray-selected samples at z > 0.5. However, based on the local relations, their X-ray luminosities and temperatures are low for their measured velocity dispersions (sigma). The clusters are cooler by a factor of 2 - 9 compared to the local sigma-T relation. We briefly discuss the possible explanations for these results.

We have discovered three globular clusters beyond the Holmberg radius in Hubble Space Telescope Advanced Camera for Surveys images of the gas-rich dark matter dominated blue compact dwarf galaxy NGC2915. The clusters, all of which start to resolve into stars, have M_{V606} = -8.9 to -9.8 mag, significantly brighter than the peak of the luminosity function of Milky Way globular clusters. Their colors suggest a metallicity [Fe/H] ~ -1.9 dex, typical of metal-poor Galactic globular clusters. The specific frequency of clusters is at a minimum normal, compared to spiral galaxies. However, since only a small portion of the system has been surveyed it is more likely that the luminosity and mass normalized cluster content is higher, like that seen in elliptical galaxies and galaxy clusters. This suggests that NGC2915 resembles a key phase in the early hierarchical assembly of galaxies - the epoch when much of the old stellar population has formed, but little of the stellar disk. Depending on the subsequent interaction history, such systems could go on to build-up larger elliptical galaxies, evolve into normal spirals, or in rare circumstances remain suspended in their development to become systems like NGC2915.

The redshift distribution of well-defined samples of distant early-type galaxies offers a means to test the predictions of monolithic and hierarchical galaxy formation scenarios. NICMOS maps of the entire Hubble Deep Field North in the F110W and F160W filters, when combined with the available WFPC2 data, allow us to calculate photometric redshifts and determine the morphological appearance of galaxies at rest-frame optical wavelengths out to z ~ 2.5. Here we report results for two subsamples of early-type galaxies, defined primarily by their morphologies in the F160W band, which were selected from the NICMOS data down to H160_{AB} < 24.0. The observed redshift distributions of our two early-type samples do not match that predicted by a monolithic collapse model, which shows an overabundance at z > 1.5. A hierarchical formation model better matches the redshift distribution of the HDF-N early-types at z > 1.5, but still does not adequately describe the observed early-types. The hierarchical model predicts significantly bluer colors on average than the observed early-type colors, and underpredicts the observed number of early-types at z < 1. [abridged]

We present the analysis of the faint galaxy population in the Advanced Camera for Surveys (ACS) Early Release Observation fields VV 29 (UGC 10214) and NGC 4676. Here we attempt to thoroughly consider all aspects relevant for faint galaxy counting and photometry, developing methods which are based on public software and that are easily reproducible by other astronomers. Using simulations we determine the best SExtractor parameters for the detection of faint galaxies in deep HST observations, paying special attention to the issue of deblending, which significantly affects the normalization and shape of the number count distribution. We confirm, as claimed by Bernstein, Freedman and Madore (2002), that Kron-like magnitudes, such as the ones generated by SExtractor, can miss more than half of the light of faint galaxies, what dramatically affects the slope of the number counts. We present catalogs for the VV 29 and NGC 4676 fields with photometry in the g,V and I bands. We also show that combining the bayesian software BPZ with superb ACS data and new spectral templates enables us to estimate reliable photometric redshifts for a significant fraction of galaxies with as few as three filters. After correcting for selection effects, we measure slopes of 0.32+- 0.01 for 22 < g < 28, 0.34+-0.01 for 22< V <27.5 and 0.33+-0.01 for 22 < I < 27. The counts do not flatten (except perhaps in the g-band), up to the depth of our observations. We find that the faint counts m_{AB}> 25.5 can be well approximated in all our filters by a passive luminosity evolution model based on the COMBO-17 luminosity function (\alpha=-1.5), with a strong merging rate following the prescription of Glazebrook et al. (1994), \phi^*\propto (1+Qz), with Q=4.

We investigate the color-magnitude (CM) relation of galaxies in the distant X-ray selected cluster RDCS1252.9-2927 at z=1.24 using images obtained with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope in the F775W and F850LP bandpasses. We select galaxies based on morphological classifications extending about 3.5 mag down the galaxy luminosity function, augmented by spectroscopic membership information. At the core of the cluster is an extensive early-type galaxy population surrounding a central pair of galaxies that show signs of dynamical interaction. The early-type population defines a tight sequence in the CM diagram, with an intrinsic scatter in observed (i-z) of $0.029\pm0.007$ mag based on 52 galaxies, or $0.024\pm0.008$ mag for $\sim30$ ellipticals. Simulations using the latest stellar population models indicate an age scatter for the ellipticals of about 34%, with a mean age $\taul\gta2.6$ Gyr (corresponding to $z_L\gta2.7$), and the last star formation occurring at $z_{end}\gta1.5$. Transforming to rest-frame $(U-B)$, we conclude that the slope and scatter in the CM relation for morphologically selected early-type galaxies show little or no evidence for evolution out to $z\approx1.2$. Thus, elliptical galaxies were already well established in X-ray luminous clusters when the universe was a third of its present age.

Using an i-z dropout criterion, we determine the space density of z~6 galaxies from two deep ACS GTO fields with deep optical-IR imaging. A total of 23 objects are found over 46 arcmin^2, or ~0.5 objects/arcmin^2 down to z~27.3 (6 sigma; all AB mag) (including one probable z~6 AGN). Combining deep ISAAC data for our RDCS1252-2927 field (J~25.7 and Ks~25.0 (5 sigma)) and NICMOS data for the HDF North (JH~27.3 (5 sigma)), we verify that these dropouts have flat spectral slopes. i-dropouts in our sample range in luminosity from ~1.5 L* (z~25.6) to ~0.3 L* (z~27.3) with the exception of one very bright candidate at z~24.2. The half-light radii vary from 0.09" to 0.29", or 0.5 kpc to 1.7 kpc. We derive the z~6 rest-frame UV luminosity density using three different procedures, each utilizing simulations based on a CDF South V dropout sample. First, we compare our findings with a no-evolution projection of this V-dropout sample. We find 23+/-25% more i-dropouts than we predict. Adopting previous results to z~5, this works out to a 20+/-29% drop in the luminosity density from z~3 to z~6. Second, we use these same V-dropout simulations to derive a selection function for our i-dropout sample and compute the UV-luminosity density (7.2+/-2.5 x 10^25 ergs/s/Hz/Mpc^3 down to z~27). We find a 39+/-21% drop over the same redshift range. This is our preferred value and suggests a star formation rate of 0.0090+/-0.0031 M_sol/yr/Mpc^3 to z~27, or ~0.036+/- 0.012 M_sol/yr/Mpc^3 extrapolating the LF to the faint limit. Third, we follow a very similar procedure, but assume no incompleteness, finding a luminosity density which is ~2-3X lower. This final estimate constitutes a lower limit. All three estimates are within the canonical range of luminosity densities necessary for reionization of the universe at this epoch. (abridged)

We present a catalog of 799 clusters of galaxies in the redshift range z_est = 0.05 - 0.3 selected from ~400 deg^2 of early SDSS commissioning data along the celestial equator. The catalog is based on merging two independent selection methods -- a color-magnitude red-sequence maxBCG technique (B), and a Hybrid Matched-Filter method (H). The BH catalog includes clusters with richness \Lambda >= 40 (Matched-Filter) and N_gal >= 13 (maxBCG), corresponding to typical velocity dispersion of \sigma_v >~ 400 km s^{-1} and mass (within 0.6 h^{-1) Mpc radius) >~ 5*10^{13} h^{-1} M_sun. This threshold is below Abell richness class 0 clusters. The average space density of these clusters is 2*10^{-5} h^3 Mpc^{-3}. All NORAS X-ray clusters and 53 of the 58 Abell clusters in the survey region are detected in the catalog; the 5 additional Abell clusters are detected below the BH catalog cuts. The cluster richness function is determined and found to exhibit a steeply decreasing cluster abundance with increasing richness. We derive observational scaling relations between cluster richness and observed cluster luminosity and cluster velocity dispersion; these scaling relations provide important physical calibrations for the clusters. The catalog can be used for studies of individual clusters, for comparisons with other sources such as X-ray clusters and AGNs, and, with proper correction for the relevant selection functions, also for statistical analyses of clusters.

We present a new sample of 35 quasars selected from the FIRST radio survey and the Deeprange I-band survey (Postman et al. 1998, 2002). A comparison with the FIRST Bright Quasar survey samples reveals that this I-band selected sample is redder by 0.25-0.5 mag in B-R, and that the color difference is not explained by the higher mean redshift of this sample but must be intrinsic. Our small sample contains five quasars with unusually red colors, including three that appear very heavily reddened. Our data are fitted well with normal blue quasar spectra attenuated by more than 2.5 magnitudes of extinction in the I-band. The red quasars are only seen at low redshifts (z<1.3). Even with a magnitude limit I<20.5, our survey is deep enough to detect only the most luminous of these red quasars at z~1; similar objects at higher redshifts would fall below our in I-band limit. Indeed, the five most luminous objects (using dereddened magnitudes) with z<1.3 are all red. Our data strongly support the hypothesis that radio quasars are dominated by a previously undetected population of red, heavily obscured objects. Unless highly reddened quasars are preferentially also highly luminous, there must be an even larger, as yet undiscovered, population of red quasars at lower luminosity. We are likely to be finding only the most luminous tip of the red quasar iceberg. Five of our six z<1 quasars are associated with Deeprange cluster candidates of similar estimated redshifts. This association is very unlikely to be the result of chance. It has some surprising implications, including the possibility that up to half of the Deeprange clusters at z~1 have associated quasars.

Multicolor coronagraphic images of the circumstellar disk around HD141569A have been obtained with the Hubble Space Telescope (HST) Advanced Camera for Surveys. B, V, and I images show that the disk's previously-described multiple-ring structure is actually a continuous distribution of dust with a tightly-wound spiral structure. Two, more open spiral arms extend from the disk, one of which appears to reach the nearby binary star HD141569BC. Diffuse dust is seen up to 1200AU from HD141569A. Although planets may exist in the inner region of the disk, tidal interaction with HD141569BC seems more likely to be the cause of these phenomena. The disk appears redder than the star (B-V = 0.21 and V-I = 0.25), and its color is spatially uniform. A scattering asymmetry factor of g = 0.25-0.35 is derived. The azimuthal density distribution is asymmetric, varying by a factor of ~3 at some radii.

We present observations of the first two supernovae discovered with the recently installed Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. The supernovae were found in Wide Field Camera images of the Hubble Deep Field North taken with the F775W, F850LP, and G800L optical elements as part of the ACS guaranteed time observation program. Spectra extracted from the ACS G800L grism exposures confirm that the objects are Type Ia supernovae (SNe Ia) at redshifts z=0.47 and z=0.95. Follow-up HST observations have been conducted with ACS in F775W and F850LP and with NICMOS in the near-infrared F110W bandpass, yielding a total of 9 flux measurements in the 3 bandpasses over a period of 50 days in the observed frame. We discuss many of the important issues in doing accurate photometry with the ACS. We analyze the multi-band light curves using two different fitting methods to calibrate the supernovae luminosities and place them on the SNe Ia Hubble diagram. The resulting distances are consistent with the redshift-distance relation of the accelerating universe model, although evolving intergalactic grey dust remains as a less likely possibility. The relative ease with which these SNe Ia were found, confirmed, and monitored demonstrates the potential ACS holds for revolutionizing the field of high-redshift SNe Ia, and therefore of testing the accelerating universe cosmology and constraining the "epoch of deceleration".

We present the first Advanced Camera for Surveys (ACS) observations of young star clusters in the colliding/merging galaxy UGC 10214. The observations were made as part of the Early Release Observation (ERO) program for the newly installed ACS during service mission SM3B for the Hubble Space Telescope (HST). Many young star clusters can be identified in the tails of UGC 10214, with ages ranging from ~3 Myr to 10 Myr. The extreme blue V-I (F606W-F814W) colors of the star clusters found in the tail of UGC 10214 can only be explained if strong emission lines are included with a young stellar population. This has been confirmed by our Keck spectroscopy of some of these bright blue stellar knots. The most luminous and largest of these blue knots has an absolute magnitude of M_V = -14.45, with a half-light radius of 161 pc, and if it is a single star cluster, would qualify as a super star cluster (SSC). Alternatively, it could be a superposition of multiple scaled OB associations or clusters. With an estimated age of ~ 4-5 Myr, its derived mass is < 1.3 x 10^6 solar masses. Thus the young stellar knot is unbound and will not evolve into a normal globular cluster. The bright blue clusters and associations are much younger than the dynamical age of the tail, providing strong evidence that star formation occurs in the tail long after it was ejected. UGC 10214 provides a nearby example of processes that contributed to the formation of halos and intra-cluster media in the distant and younger Universe.

The mass function of clusters of galaxies is determined from 400 deg^2 of early commissioning imaging data of the Sloan Digital Sky Survey; ~300 clusters in the redshift range z = 0.1 - 0.2 are used. Clusters are selected using two independent selection methods: a Matched Filter and a red-sequence color magnitude technique. The two methods yield consistent results. The cluster mass function is compared with large-scale cosmological simulations. We find a best-fit cluster normalization relation of sigma_8*omega_m^0.6 = 0.33 +- 0.03 (for 0.1 ~< omega_m ~< 0.4), or equivalently sigma_8 = (0.16/omega_m)^0.6. The amplitude of this relation is significantly lower than the previous canonical value, implying that either omega_m is lower than previously expected (omega_m = 0.16 if sigma_8 = 1) or sigma_8 is lower than expected (sigma_8 = 0.7 if omega_m = 0.3). The best-fit mass function parameters are omega_m = 0.19 (+0.08,-0.07) and sigma_8 = 0.9 (+0.3,-0.2). High values of omega_m (>= 0.4) and low sigma_8 (=< 0.6) are excluded at >~ 2 sigma.

We used the HST WFPC2 to obtain I-band images of the centers of 81 brightest cluster galaxies (BCGs), drawn from a volume-limited sample of nearby BCGs. The images show a rich variety of morphological features, including multiple or double nuclei, dust, stellar disks, point source nuclei, and central surface brightness depressions. High resolution surface brightness profiles could be inferred for 60 galaxies. Of those, 88% have well-resolved cores. Twelve percent of the BCG sample lacks a well-resolved core; all but one of these BCGs have ``power-law'' profiles. Some of these galaxies have higher luminosities than any power-law galaxy identified by Faber et al. (1997), and have physical upper limits on the break radius well below the values observed for core galaxies of the same luminosity. These results support the idea that the central structure of early-type galaxies is bimodal in its physical properties, but also suggest that there exist high luminosity galaxies with power-law profiles (or unusually small cores). The BCGs in the latter category tend to fall at the low end of the BCG luminosity function and tend to have low values of the quantity alpha (the logarithmic slope of the metric luminosity as a function of radius, at 10 kpc). Since theoretical calculations have shown that the luminosities and alpha values of BCGs grow with time as a result of accretion, this suggests a scenario in which elliptical galaxies evolve from power-law profiles to core profiles through accretion and merging. This is consistent with theoretical scenarios that invoke the formation of massive black hole binaries during merger events (Abridged).

ClassX is a project aimed at creating an automated system to classify X-ray sources and is envisaged as a prototype of the Virtual Observatory. As a system, ClassX integrates into a pipeline a network of classifiers and an engine that searches and retrieves for a given target multi-wavelength counterparts from the worldwide data storage media. It applies machine learning methods to `train' different classifiers using different `training' data sets. In ClassX, each classifier can make its own class (object type) assignment and is optimized for handling different tasks and/or different object types. A user would generally select a certain classifier to make, for instance, a most complete list of candidate QSOs, but a different classifier would be used to make a most reliable list of candidate QSOs. Still different classifiers would be selected to make similar lists for other object types. Along with the class name assignment, a network classifier outputs the probability for a source to belong to the assigned class as well as probabilities that the source belongs in fact to other classes. We illustrate the current capabilities of ClassX and the concept of a classifiers network with the results obtained with classifiers trained using ROSAT data. ~

We have conducted an automated search for galaxy clusters within a contiguous 16 square degree I-band survey in the north Galactic hemisphere. A matched filter detection algorithm identifies 444 cluster candidates in the range 0.2 <= z <= 1.2. The full catalog is presented along with the results from a follow-up spectroscopic survey. The estimated redshift distribution of the cluster candidates is consistent with a constant comoving density over the range 0.2 <= z_est <= 0.8. A decline in the cluster space density by more than a factor of 3 over this redshift range is rejected at >99.9% confidence level. We find that the space density of Lambda_CL >= 40 clusters in our survey lies in the range (1.6 - 1.8) x 10^{-5} h_{75}^{3} Mpc^{-3}, ~1.5 times higher than the local distribution of comparably rich Abell RC >= 0 clusters. The Lambda_CL distribution is consistent with a power-law. The discrepancy between the space density of Abell clusters and the clusters in this survey declines quickly as Lambda_CL increases, suggesting that the difference at lower richness is due to significant incompleteness in the Abell catalog. A percolation analysis reveals that 10 - 20% of the spectroscopically confirmed distant clusters are linked into superclusters at overdensities between 10 < {\delta\rho \over \rho} < 50, similar to what is seen in the local cluster distribution. This suggests that there has been little evolution of the cluster-cluster correlation length for z <= 0.5.

We have carried out additional spectroscopic observations in the field of cluster Cl 1324+3011 at z = 0.76. Combined with the spectroscopy presented in Postman, Lubin & Oke (2001, AJ, 122, 1125), we now have spectroscopically confirmed 47 cluster members. With this significant number of redshifts, we measure accurately the cluster velocity dispersion to be 1016 (+126/-93) km/s. The distribution of velocity offsets is consistent with a Gaussian, indicating no substantial velocity substructure. As previously noted for other optically-selected clusters at redshifts of z > 0.5, a comparison between the X-ray luminosity (L_x) and the velocity dispersion (sigma) of Cl 1324+3011 implies that this cluster is underluminous in X-rays by a factor of ~3 - 40 when compared to the L_x - sigma relation for local and moderate-redshift clusters. We also examine the morphologies of those cluster members which have available high-angular-resolution imaging with the Hubble Space Telescope (HST). There are 22 spectroscopically-confirmed cluster members within the HST field-of-view. Twelve of these are visually classified as early-type (elliptical or S0) galaxies, implying an early-type fraction of 0.55 (+0.17/-0.14) in this cluster. This fraction is a factor of ~1.5 lower than that observed in nearby rich clusters. Confirming previous cluster studies, the results for cluster Cl 1324+3011, combined with morphological studies of other massive clusters at redshifts less than z = 1, suggest that the galaxy population in massive clusters is strongly evolving with redshift. This evolution implies that early-type galaxies are forming out of the excess of late-type (spiral, irregular, and peculiar) galaxies over the ~7 Gyr timescale.

We used HST WFPC2 images to identify six early-type galaxies with surface- brightness profiles that decrease inward over a limited range of radii near their centers. The implied luminosity density profiles of these galaxies have local minima interior to their core break radii. NGC 3706 harbors a high surface brightness ring of starlight with radius ~20 pc. Its central structure may be related to that in the double-nucleus galaxies M31 and NGC 4486B. NGC 4406 and NGC 6876 have nearly flat cores that on close inspection are centrally depressed. Colors for both galaxies imply that this is not due to dust absorption. The surface brightness distributions of both galaxies are consistent with stellar tori that are more diffuse than the sharply defined system in NGC 3706. The remaining three galaxies are the brightest cluster galaxies in A260, A347, and A3574. Color information is not available for these objects, but they strongly resemble NGC 4406 and NGC 6876 in their cores. The thin ring in NGC 3706 may have formed dissipatively. The five other galaxies resemble the endpoints of some simulations of the merging of two gas-free stellar systems, each harboring a massive nuclear black hole. In one version of this scenario, diffuse stellar tori are produced when stars initially bound to one black hole are tidally stripped away by the second black hole. Alternatively, some inward-decreasing surface-brightness profiles may reflect the ejection of stars from a core during the hardening of the binary black hole created during the merger.

We present a comparison of X-ray and optical luminosities and luminosity functions of cluster candidates from a joint optical/X-ray survey, the ROSAT Optical X-Ray Survey. Completely independent X-ray and optical catalogs of 23 ROSAT fields (4.8 deg2) were created by a matched-filter optical algorithm and by a wavelet technique in the X-ray. We directly compare the results of the optical and X-ray selection techniques. The matched-filter technique detected 74% (26 out of 35) of the most reliable cluster candidates in the X-ray-selected sample; the remainder could be either constellations of X-ray point sources or z>1 clusters. The matched-filter technique identified approximately 3 times the number of candidates (152 candidates) found in the X-ray survey of nearly the same sky (57 candidates). While the estimated optical and X-ray luminosities of clusters of galaxies are correlated, the intrinsic scatter in this relationship is very large. We can reproduce the number and distribution of optical clusters with a model defined by the X-ray luminosity function and by an LX Lambda cl relation if H0=75 km s-1 Mpc-1 and if the LX Lambda cl relation is steeper than the expected LX Lambda 2cl. On statistical grounds, a bimodal distribution of X-ray luminous and X-ray faint clusters is unnecessary to explain our observations. Follow-up work is required to confirm whether the clusters without bright X-ray counterparts are simply X-ray faint for their optical luminosity because of their low mass or youth or are a distinct population of clusters that do not, for some reason, have dense intracluster media. We suspect that these optical clusters are low-mass systems, with correspondingly low X-ray temperatures and luminosities, or that they are not yet completely virialized systems.

We present and analyze the optical and X-ray catalogs of moderate-redshift cluster candidates from the ROSAT Optical X-ray Survey, or ROXS. The survey covers 4.8 square degrees (23 ROSAT PSPC pointings). The cross-correlated cluster catalogs were constructed by comparing two independent catalogs extracted from the optical and X-ray bandpasses, using a matched filter technique for the optical data and a wavelet technique for the X-ray data. We cross-id cluster candidates in each catalog. In Paper II we present the cluster catalogs and a numerical simulation of ROXS. We also present color-magnitude plots for several cluster candidates, and examine the prominence of the red sequence in each. We find that the X-ray clusters analyzed in this way do not all have a prominent red sequence. We conclude that while the red sequence may be distinct for massive, virialized clusters, it may be less so for lower-mass clusters at even moderate redshifts. Multiple, complementary methods of selecting and defining clusters may be essential, particularly at high redshift where all methods run into completeness limits, incomplete understanding of physical evolution, and projection effects.

Given the importance of clusters to the fields of cosmology and galaxy evolution, it is critical to understand how the cluster detection process affects (biases) ones scientific conclusions derived from a given cluster sample. I review the astrophysics, algorithms, and observational constraints that must be considered when attempting to assess cluster selection biases and their impact on the scientific constraints derived from multi-wavelength cluster surveys. In particular, one can more accurately quantify and understand the selection biases and assess cluster evolution when using joint optical/NIR/x-ray surveys than can be achieved when employing any single cluster survey.

One of the most vital observational clues for unraveling the origin of Brightest Cluster Galaxies (BCG) is the observed alignment of the BCGs with their host cluster and its surroundings. We have examined the BCG-cluster alignment effect, using clusters of galaxies detected from the Sloan Digital Sky Survey (SDSS). We find that the BCGs are preferentially aligned with the principal axis of their hosts, to a much higher redshift (z >~ 0.3) than probed by previous studies (z <~ 0.1). The alignment effect strongly depends on the magnitude difference of the BCG and the second and third brightest cluster members: we find a strong alignment effect for the dominant BCGs, while less dominant BCGs do not show any departure from random alignment with respect to the cluster. We therefore claim that the alignment process originates from the same process that makes the BCG grow dominant, be it direct mergers in the early stage of cluster formation, or a later process that resembles the galactic cannibalism scenario. We do not find strong evidence for (or against) redshift evolution between 0<z<0.45, largely due to the insufficient sample size (< 200 clusters). However, we have developed a framework by which we can examine many more clusters in an automated fashion for the upcoming SDSS cluster catalogs, which will provide us with better statistics for systematic investigations of the alignment with redshift, richness and morphology of both the cluster and the BCG.

We present a comparison of three cluster finding algorithms from imaging data using Monte Carlo simulations of clusters embedded in a 25 deg^2 region of Sloan Digital Sky Survey (SDSS) imaging data: the Matched Filter (MF; Postman et al. 1996), the Adaptive Matched Filter (AMF; Kepner et al. 1999) and a color-magnitude filtered Voronoi Tessellation Technique (VTT). Among the two matched filters, we find that the MF is more efficient in detecting faint clusters, whereas the AMF evaluates the redshifts and richnesses more accurately, therefore suggesting a hybrid method (HMF) that combines the two. The HMF outperforms the VTT when using a background that is uniform, but it is more sensitive to the presence of a non-uniform galaxy background than is the VTT; this is due to the assumption of a uniform background in the HMF model. We thus find that for the detection thresholds we determine to be appropriate for the SDSS data, the performance of both algorithms are similar; we present the selection function for each method evaluated with these thresholds as a function of redshift and richness. For simulated clusters generated with a Schechter luminosity function (M_r^* = -21.5 and alpha = -1.1) both algorithms are complete for Abell richness >= 1 clusters up to z ~ 0.4 for a sample magnitude limited to r = 21. While the cluster parameter evaluation shows a mild correlation with the local background density, the detection efficiency is not significantly affected by the background fluctuations, unlike previous shallower surveys.

We describe the target selection and resulting properties of a spectroscopic sample of luminous, red galaxies (LRG) from the imaging data of the Sloan Digital Sky Survey (SDSS). These galaxies are selected on the basis of color and magnitude to yield a sample of luminous, intrinsically red galaxies that extends fainter and further than the main flux-limited portion of the SDSS galaxy spectroscopic sample. The sample is designed to impose a passively-evolving luminosity and rest-frame color cut to a redshift of 0.38. Additional, yet more luminous, red galaxies are included to a redshift of 0.5. Approximately 12 of these galaxies per square degree are targeted for spectroscopy, so the sample will number over 100,000 with the full survey. SDSS commissioning data indicate that the algorithm efficiently selects luminous (M_g=-21.4), red galaxies, that the spectroscopic success rate is very high, and that the resulting set of galaxies is approximately volume-limited out to z=0.38. When the SDSS is complete, the LRG spectroscopic sample will fill over 1h^-3 Gpc^3 with an approximately homogeneous population of galaxies and will therefore be well suited to studies of large-scale structure and clusters out to z=0.5.

The Sloan Digital Sky Survey is one of the first multicolor photometric and spectroscopic surveys designed to measure the statistical properties of galaxies within the local Universe. In this Letter we present some of the initial results on the angular 2-point correlation function measured from the early SDSS galaxy data. The form of the correlation function, over the magnitude interval 18<r*<22, is shown to be consistent with results from existing wide-field, photographic-based surveys and narrower CCD galaxy surveys. On scales between 1 arcminute and 1 degree the correlation function is well described by a power-law with an exponent of ~ -0.7. The amplitude of the correlation function, within this angular interval, decreases with fainter magnitudes in good agreement with analyses from existing galaxy surveys. There is a characteristic break in the correlation function on scales of approximately 1-2 degrees. On small scales, < 1', the SDSS correlation function does not appear to be consistent with the power-law form fitted to the 1'< theta <0.5 deg data. With a data set that is less than 2% of the full SDSS survey area, we have obtained high precision measurements of the power-law angular correlation function on angular scales 1' < theta < 1 deg, which are robust to systematic uncertainties. Because of the limited area and the highly correlated nature of the error covariance matrix, these initial results do not yet provide a definitive characterization of departures from the power-law form at smaller and larger angles. In the near future, however, the area of the SDSS imaging survey will be sufficient to allow detailed analysis of the small and large scale regimes, measurements of higher-order correlations, and studies of angular clustering as a function of redshift and galaxy type.

New photometric and spectroscopic observations of galaxies in the directions of three distant clusters are presented as part of our on-going high-redshift cluster survey. The clusters are CL1324+3011 at z = 0.76, CL1604+4304 at z = 0.90, and CL1604+4321 at z = 0.92. The observed x-ray luminosities in these clusters are at least a factor of 3 smaller than those observed in clusters with similar velocity dispersions at z <= 0.4. These clusters contain a significant population of elliptical-like galaxies, although these galaxies are not nearly as dominant as in massive clusters at z <= 0.5. We also find a large population of blue cluster members. Defining an active galaxy as one in which the rest equivalent width of [OII] is greater than 15 Angstroms, the fraction of active cluster galaxies, within the central 1.0 Mpc, is 45%. In the field population, we find that 65% of the galaxies with redshifts between z = 0.40 and z = 0.85 are active, while the fraction is 79% for field galaxies at z > 0.85. The star formation rate normalized by the rest AB B-band magnitude, SFRN, increases as the redshift increases at a given evolving luminosity. At a given redshift, however, SFRN decreases linearly with increasing luminosity indicating a remarkable insensitivity of the star formation rate to the intrinsic luminosity of the galaxy over the range -18 >= ABB >= -22. Cluster galaxies in the central 1 Mpc regions exhibit depressed star formation rates. We are able to measure significant evolution in the B-band luminosity function over the range 0.1 <= z <= 1. The characteristic luminosity increases by a factor of 3 with increasing redshift over this range.

We outline here the next generation of cluster-finding algorithms. We show how advances in Computer Science and Statistics have helped develop robust, fast algorithms for finding clusters of galaxies in large multi-dimensional astronomical databases like the Sloan Digital Sky Survey (SDSS). Specifically, this paper presents four new advances: (1) A new semi-parametric algorithm - nicknamed ``C4'' - for jointly finding clusters of galaxies in the SDSS and ROSAT All-Sky Survey databases; (2) The introduction of the False Discovery Rate into Astronomy; (3) The role of kernel shape in optimizing cluster detection; (4) A new determination of the X-ray Cluster Luminosity Function which has bearing on the existence of a ``deficit'' of high redshift, high luminosity clusters. This research is part of our ``Computational AstroStatistics'' collaboration (see Nichol et al. 2000) and the algorithms and techniques discussed herein will form part of the ``Virtual Observatory'' analysis toolkit.

We measured infrared surface brightness fluctuation (SBF) distances to an isotropically-distributed sample of 16 distant galaxies with redshifts reaching 10,000 km/s using the near-IR camera and multi-object spectrometer (NICMOS) on the Hubble Space Telescope (HST). The excellent spatial resolution, very low background, and brightness of the IR fluctuations yielded the most distant SBF measurements to date. Twelve nearby galaxies were also observed and used to calibrate the F160W (1.6 micron) SBF distance scale. Of these, three have Cepheid variable star distances measured with HST and eleven have optical I-band SBF distance measurements. A distance modulus of 18.5 mag to the Large Magellanic Cloud was adopted for this calibration. We present the F160W SBF Hubble diagram and find a Hubble constant Ho=76 +/- 1.3 (1-sigma statistical) +/- 6 (systematic) km/s/Mpc. This result is insensitive to the velocity model used to correct for local bulk motions. Restricting the fit to the six most distant galaxies yields the smallest value of Ho=72 +/- 2.3 km/s/Mpc consistent with the data. This 6% decrease in the Hubble constant is consistent with the hypothesis that the Local Group inhabits an under-dense region of the universe, but is also consistent with the best-fit value of Ho=76 km/s/Mpc at the 1.5-sigma level.

Constraints on the validity of the hierarchical gravitational instability theory and the evolution of biasing are presented based upon measurements of higher order clustering statistics in the Deeprange Survey, a catalog of $\sim710,000$ galaxies with $I_{AB} \le 24$ derived from a KPNO 4m CCD imaging survey of a contiguous $4^{\circ} \times 4^{\circ}$ region. We compute the 3-point and 4-point angular correlation functions using a direct estimation for the former and the counts-in-cells technique for both. The skewness $s_3$ decreases by a factor of $\simeq 3-4$ as galaxy magnitude increases over the range $17 \le I \le 22.5$ ($0.1 \lesssim z \lesssim 0.8$). This decrease is consistent with a small {\it increase} of the bias with increasing redshift, but not by more than a factor of 2 for the highest redshifts probed. Our results are strongly inconsistent, at about the $3.5-4 \sigma$ level, with typical cosmic string models in which the initial perturbations follow a non-Gaussian distribution - such models generally predict an opposite trend in the degree of bias as a function of redshift. We also find that the scaling relation between the 3-point and 4-point correlation functions remains approximately invariant over the above magnitude range. The simplest model that is consistent with these constraints is a universe in which an initially Gaussian perturbation spectrum evolves under the influence of gravity combined with a low level of bias between the matter and the galaxies that decreases slightly from $z \sim 0.8$ to the current epoch.

We present the first results from the Shellflow program, an all-sky Tully-Fisher (TF) peculiar velocity survey of 276 Sb-Sc galaxies with redshifts between 4500 and 7000 km/s. Shellflow was designed to minimize systematic errors between observing runs and between telescopes, thereby removing the possibility of a spurious bulk flow caused by data inhomogeneity. A fit to the data yields a bulk flow amplitude V_bulk = 70{+100}{-70} km/s (1 sigma error) with respect to the Cosmic Microwave Background, i.e., consistent with being at rest. At the 95% confidence level, the flow amplitude is < 300 km/s. Our results are insensitive to which Galactic extinction maps we use, and to the parameterization of the TF relation. The larger bulk motion found in analyses of the Mark III peculiar velocity catalog are thus likely to be due to non-uniformities between the subsamples making up Mark III. The absence of bulk flow is consistent with the study of Giovanelli and collaborators and flow field predictions from the observed distribution of IRAS galaxies.

We present the first dynamical study of the optically selected Palomar Distant Cluster Survey (PDCS). We show that eleven out of twelve PDCS clusters we observed are real over-densities of galaxies. Most clusters have velocity dispersions appropriate for clusters of galaxies. However, we find a fraction (1/3) of objects in the PDCS which have velocity dispersions in the range of groups of galaxies (200 +/- 100 km/s) but have richnesses appropriate for clusters of galaxies. We measure the richness function, X-ray luminosity function (using both the detections and upper limits), and the mass function derived from our velocity dispersions. We confirm that the space density, as a function of richness, of clusters of galaxies in the PDCS is ~5 times that of the Abell catalog. Excluding the above fraction of 1/3 of objects with low velocity dispersions, we measure a space density ~3 times that of the Abell catalog for equivalent mass clusters of galaxies. Our space density estimates are in agreement with other low-redshift, optically-selected cluster surveys such as the EDCC, APM and EDCC2, which raises the possibility that the Abell catalog is incomplete. Our X-ray luminosity function agrees with other measurements based on both X-ray and optically selected samples, so we find that the PDCS does not miss clusters of galaxies that would be found in an X-ray selected survey. Our resulting mass function, centered around 10e14 Msun/h, agrees with the expectations from such surveys as the Canadian Network for Observational Cosmology cluster survey, though errors on our mass measurements are too large to constrain cosmological parameters.

This paper presents and gives the COP (COP: CFHT Optical PDCS; CFHT: Canada-France-Hawaii Telescope; PDCS: Palomar Distant Cluster Survey) survey data. We describe our photometric and spectroscopic observations with the MOS multi-slit spectrograph at the CFH telescope. A comparison of the photometry from the PDCS (Postman et al. 1996) catalogs and from the new images we have obtained at the CFH telescope shows that the different magnitude systems can be cross-calibrated. After identification between the PDCS catalogues and our new images, we built catalogues with redshift, coordinates and V, I and Rmagnitudes. We have classified the galaxies along the lines of sight into field and structure galaxies using a gap technique (Katgert et al. 1996). In total we have observed 18 significant structures along the 10 lines of sight.

We present the results from a multi-band optical imaging program which has definitively confirmed the existence of a supercluster at z = 0.91. Two massive clusters of galaxies, CL1604+4304 at z = 0.897 and CL1604+4321 at z = 0.924, were originally observed in the high-redshift cluster survey of Oke, Postman & Lubin (1998). They are separated by 4300 km/s in radial velocity and 17 arcminutes on the plane of the sky. Their physical and redshift proximity suggested a promising supercluster candidate. Deep BRi imaging of the region between the two clusters indicates a large population of red galaxies. This population forms a tight, red sequence in the color--magnitude diagram at (R-i) = 1.4. The characteristic color is identical to that of the spectroscopically-confirmed early-type galaxies in the two member clusters. The red galaxies are spread throughout the 5 Mpc region between CL1604+4304 and CL1604+4321. Their spatial distribution delineates the entire large scale structure with high concentrations at the cluster centers. In addition, we detect a significant overdensity of red galaxies directly between CL1604+4304 and CL1604+4321 which is the signature of a third, rich cluster associated with this system. The strong sequence of red galaxies and their spatial distribution clearly indicate that we have discovered a supercluster at z = 0.91.

Cosmological simulations predict that a large fraction of the baryonic mass of the Universe exists as 10^5-10^7 K diffuse, X-ray emitting gas, tracing low density filament and sheet-like structures exterior to massive clusters of galaxies. If present, this gas helps reconcile the current shortfall in observed baryon counts relative to the predictions of the standard Big Bang model. We present here the discovery and analysis of a 5-sigma significance half-degree filamentary structure, present in both I-band galaxy surface density and unresolved X-ray emission in a deep ROSAT PSPC field. The estimated diffuse X-ray emission component of this structure has a surface brightness of \~ 1.6 x 10^(-16) erg/s/cm^2/arcmin^2 (0.5- 2 keV), comparable to the predictions for inter-cluster gas and may represent a direct detection of this currently unconfirmed baryonic component.

Our understanding of the cosmic history of galaxy clusters has recently been enhanced due to an extensive series of observations including faint spectroscopic data (especially those obtained at the Keck Observatory), deep optical and NIR imaging from the ground and in space, morphological data from HST, and new constraints on the evolution of the intracluster medium from ROSAT and ASCA. When such observations are applied to complete, objectively derived catalogs of clusters, our constraints on cluster formation and evolution become quite confined. A picture is emerging in which the bulk of cluster formation starts at z >= 2, yielding cluster potentials that are well established by z ~ 1, and there is little substantial evolution of the cluster galaxy population since z ~ 0.4. This review talk will summarize the current observational constraints on the properties and evolution of high redshift clusters and protoclusters and their implications.

We present preliminary results from the Shellflow program, an all-sky Tully-Fisher (TF) survey of 297 Sb-Sc with redshifts between 4000 and 7500 km/s. The program was designed to ensure uniformity of the data between observing runs and between telescopes, thereby eliminating the possibility of a spurious bulk flow caused by data inhomogeneity. A simple bulk flow model suggests that the Shellflow galaxies are nearly at rest with respect to the cosmic microwave background (CMB). Taken at face value, this result suggests that most of the ~600 km/s motion of the Local Group with respect to the CMB is due to material within ~60 Mpc/h -- a striking confirmation of the gravitational instability picture and the notion of homogeneity in the Universe. Differences between our Shellflow ``null'' result and claims of larger bulk motions based on analyses of the MarkIII catalog are possibly due to small inhomogeneities between the sub-samples in MarkIII. Our preliminary Shellflow bulk velocity, V_bulk = 80+/-150 km/s in the CMB frame, must still be refined with Monte-Carlo simulations and tidal field analysis.

We have conducted a redshift survey of sixteen cluster candidates from the Palomar Distant Cluster Survey (PDCS) to determine both the density of PDCS clusters and the accuracy of the estimated redshifts presented in the PDCS catalog (Postman et. al. 1996). We find that the matched-filter redshift estimate presented in the PDCS has an error sigma_z = 0.06 in the redshift range 0.1 < z < 0.35 based on eight cluster candidates with three or more concordant galaxy redshifts. We measure the low redshift (0.1 < z < 0.35) space density of PDCS clusters to be 31.3^{+30.5}_{-17.1} * E-06 h^3 Mpc^-3 (68% confidence limits for a Poisson distribution) for Richness Class 1 systems. We find a tentative space density of 10.4^{+23.4}_{-8.4}* E-06 h^3 Mpc^-3 for Richness Class 2 clusters. These densities compare favorably with those found for the whole of the PDCS and support the finding that the space density of clusters in the PDCS is a factor of ~5 above that of clusters in the Abell catalog (Abell 1958; Abell, Corwin, and Olowin 1989). These new space density measurements were derived as independently as possible from the original PDCS analysis and therefore, demonstrate the robustness of the original work. Based on our survey, we conclude that the PDCS matched-filter algorithm is successful in detecting real clusters and in estimating their true redshifts in the redshift range we surveyed.

We present the first results of the COP survey about the reality of the PDCS clusters, about their velocity dispersions and dynamic and about the periodicity of the structures along the line of sight.

By virtue of their high galaxy space densities and their large spatial separations, clusters are efficient and accurate tracers of the large-scale density and velocity fields. Substantial progress has been made over the past decade in the construction of homogeneous, objectively derived cluster catalogs and in characterizing the spatial distribution of clusters. Consequently, the constraints on viable models for the growth of structure have been refined. A review of the status of cluster-based observations of large-scale structure is presented here, including discussions of the second and higher order moments, the dependence of clustering on richness (mass), recent and new measurements of bulk flows, and a new constraint on the cluster mass function in the range 0.7 < z < 1.

We present a new optical Tully-Fisher (TF) investigation for a complete, full-sky sample of 297 Sb-Sc spirals with redshifts between 4500 and 7000 km/s. The survey was specifically designed to provide uniform, well- calibrated data over both hemispheres. All previous TF surveys within the Supergalactic shell (cz < 6000 km/s) have relied on matching separate data sets in the Northern and Southern hemispheres and thus cannot attain full- sky homogeneity. Analyses of the cosmological dipole and peculiar velocities based on these studies have produced contradictory claims for the amplitude of the bulk flow and whether it is generated by internal or external mass fluctuations. With Shellflow, and further zero-point calibration of existing TF data sets, we expect a high-accuracy detection of the bulk flow amplitude and an unambiguous characterization of the tidal field at 6000 km/s. One of the key features of this study is not only its careful, all-sky sample selection but the independent duplication, by at least 2 of us, of all the data reductions. These reductions and a first flow analysis based on the Shellflow sample alone will be published in the Fall.

We report the serendipitous discovery of an extremely red object in the Hubble Deep Field South (HDFS) Test NICMOS (Near Infrared Camera and Multi Object Spectrograph) field of view. The object is resolved in the NICMOS image and has a light profile very well described by an r^{1/4} law with effective radius r_e=0."20+-0."05 and H_{AB}=21.7+-0.1 magnitudes. In contrast, the galaxy is undetected in the R and I band ground based images taken at the CTIO 4 m Blanco Telescope, giving a lower limit to the color of (R-H)_{AB}>3.9 and (I-H)_{AB}>3.5 at the 95% confidence level. The colors of a range of synthetic galactic spectra are computed, showing that the object is likely to be an ``old'' elliptical galaxy at redshift z>~1.7. Alternatively the colors can be reproduced by an ``old'' elliptical galaxy at somewhat lower redshift (z>~1) with significant amount of dust, or by a younger galaxy at higher redshift. This object represents a very interesting target for future VLT observations.

We present constraints on the evolution of large-scale structure from a catalog of 710,000 galaxies with I_AB <= 24 derived from a KPNO 4m CCD imaging survey of a contiguous 4 deg x 4 deg region. The advantage of using large contiguous surveys for measuring clustering properties on even modest angular scales is substantial: the effects of cosmic scatter are strongly suppressed. We provide highly accurate measurements of the two-point angular correlation function, w(theta), as a function of magnitude on scales up to 1.5 degrees. The amplitude of w(theta) declines by a factor of ~10 over the range 16 <= I <= 20 but only by a factor of 2 - 3 over the range 20 < I <= 23. For a redshift dependence of the spatial correlation function, xi(r), parameterized as xi(r,z)=(r/r_o)^(-gamma)(1 + z)^(-[3+epsilon]), we find r_o=5.2 +/- 0.4 Mpc/h, and epsilon >= 0 for I <= 20. This is in good agreement with the results from local redshift surveys. At I > 20, our best fit values shift towards lower r_o and more negative epsilon. A strong covariance between r_o and epsilon prevent us from rejecting epsilon > 0 even at faint magnitudes but if epsilon > 1, we strongly reject r_o <= 4/h Mpc (co-moving). The above expression for xi(r,z) and our data give a correlation length of r_o(z=0.5) approx 3.0 +/- 0.4 Mpc/h, about a factor of 2 larger than the correlation length at z = 0.5 derived from the Canada--France Redshift Survey (CFRS). The small volume sampled by the CFRS and other deep redshift probes, however, make these spatial surveys strongly susceptible to cosmic scatter and will tend to bias their derived correlation lengths low. Our galaxy counts agree well with those from the HDF survey and, thus, argue against a significant inclusion of sub-galactic components in the latter census for I < 24.

We present an extensive photometric and spectroscopic study of two high-redshift clusters of galaxies based on data obtained from the Keck 10m telescopes and the Hubble Space Telescope. The clusters CL0023+0423 (z=0.84) and CL1604+4304 (z=0.90) are part of a multi-wavelength program to study nine candidate clusters at z > 0.6 (Oke, Postman & Lubin 1998). Based on these observations, we study in detail both the field and cluster populations. From the confirmed cluster members, we find that CL0023+0423 actually consists of two components separated by ~2900 km/s. A kinematic analysis indicates that the two components are a poor cluster with ~3 x 10^{14} solar masses and a less massive group with 10^{13} solar masses. CL1604+4304 is a centrally concentrated, rich cluster at z = 0.8967 with a velocity dispersion of 1226 km/s and a mass of ~3 x 10^{15} solar masses. Approximately 57% and 50% of the galaxies show high levels of star formation in CL0023+0423 and CL1604+4304, respectively. These numbers are significantly larger than those found in intermediate redshift clusters. We also observe many old, red galaxies. Found mainly in CL1604+4304, they have spectra consistent with passive stellar evolution, typical of the early-type galaxies in low and intermediate-redshift clusters. We have calculated their ages by comparing their spectral energy distributions to standard Bruzual & Charlot evolutionary models. We find that their colors are consistent with models having an exponentially decreasing star formation rate with a time constant of 0.6 Gyr. We observe a significant luminosity brightening in our brightest cluster galaxies. Compared to BCGs at z < 0.1, we find a luminosity increase of ~1 mag in the rest B-band and ~0.8 mag in the rest V-band.

We present a dynamical study of the CL0023+0423 system at a redshift of z = 0.84. This system consists of two components separated in velocity space by 2900 km/s and on the plane of the sky by 0.23 h^{-1} Mpc. A kinematic analysis indicates that the two components are a poor cluster with a velocity dispersion of 415 km/s and a less massive group with a velocity dispersion of 158 km/s (Postman, Lubin & Oke 1998). An analytic two-body calculation and N-body simulations indicate that the system is most likely not bound but simply a chance projection on the sky; however, within the observational uncertainties, there do exist bound solutions where the two galaxy groups are currently moving toward each other and will eventually merge. A statistical analysis of 1000 N-body simulations indicates that there is an 20% chance that the two groups will merge. If the CL0023+0423 system does merge, it will appear as a cluster on the sky, as well as in velocity space, within 1-2 Gyrs. The cluster will evolve dynamically for more than 3 Gyrs, appearing during this time more similar to an open, irregular cluster. The final merged system has a velocity dispersion which is consistent with a local Abell richness class 1 cluster. The morphological analysis of the galaxy populations of CL0023+0423 suggests that both groups are largely dominated by spiral galaxies with early-type fractions of 33% or less (Lubin et al. 1998). These modest early-type fractions may have implications for both cluster formation and group evolution. A comparison with nearby and intermediate-redshift groups and clusters of galaxies suggest that some fraction of early-type galaxies are formed at redshifts of z < 1.

We present a description of the observations and data reduction procedures for an extensive spectroscopic and multi-band photometric study of nine high redshift, optically-selected cluster candidates. The primary goal of the survey is to establish new constraints on cluster and galaxy evolution, with specific emphasis on the evolution of galaxy morphology and on the star-formation history of the galaxies within and around distant clusters. We have measured 892 new redshifts for galaxies with R <= 23.3. The data will also serve as deep probes of the foreground and background large-scale structures. The observations include broad band optical imaging and spectroscopy with the Low Resolution Imaging Spectrograph at the 10 meter W. M. Keck Observatory telescope; K-band imaging with IRIM at the 4 meter Kitt Peak National Observatory telescope; and deep, high angular resolution imaging with the WFPC2 onboard the Hubble Space Telescope. We also describe the procedures used to obtain morphological information. We have established that six of the nine cluster candidates are indeed real space density enhancements and are representative of those typically associated with clusters of galaxies. The remaining three candidates appear to be projections of several smaller groups at widely separated distances. This success rate is consistent with estimates of the false positive rate in 2D optical high-z cluster searches.

We present a detailed morphological analysis of the galaxy populations in the first two clusters to be completed in an extensive observational study of nine high-redshift clusters of galaxies (Oke, Postman & Lubin 1998). These two clusters, CL0023+0423 and CL1604+4304, are at redshifts of z = 0.84 and z = 0.90, respectively. The morphological studies are based on high-angular resolution imagery taken with WFPC2 aboard the Hubble Space Telescope. These data are combined with deep, ground-based BVRI photometry and spectra taken with the Keck 10-meter telescopes. Based on these observations, CL0023+0423 has a galaxy population which is more similar to groups of galaxies and the field. This system is almost completely dominated by spiral galaxies. CL1604+4304, however, has a morphological composition which is more typical of a normal, present-day cluster; early-type galaxies comprise approximately 76% of all galaxies brighter than M_V = -19.0 + 5 log h in the central 0.5 h^{-1} Mpc. The ratio of S0 galaxies to ellipticals in this cluster is 1.7 +/- 0.9, consistent with local cluster populations. The morphological results support the conclusions of the dynamical analysis presented in the second paper of this series (Postman, Lubin & Oke 1998). CL0023+0423 consists of two galaxy groups which are separated by approximately 2900 km/s in radial velocity. On the other hand, CL1604+4304 has a morphological composition, velocity dispersion, and implied mass which are consistent with an Abell richness class 2 or 3 cluster.

We used HST to obtain surface brightness fluctuation (SBF) observations of four nearby brightest cluster galaxies (BCG) to calibrate the BCG Hubble diagram of Lauer & Postman (1992). This BCG Hubble diagram contains 114 galaxies covering the full celestial sphere and is volume limited to 15,000 km/s, providing excellent sampling of the far field Hubble flow. The SBF zero point is based on the Cepheid calibration of the ground I_KC method (Tonry et al. 1997) as extended to the WFPC2 F814W filter by Ajhar et al. (1997). The BCG globular cluster luminosity functions give distances essentially identical to the SBF results. Using the velocities and SBF distances of the four BCG alone gives H_0 = 82 +/- 8 km/s/Mpc in the CMB frame, valid on ~4,500 km/s scales. Use of BCG as photometric redshift estimators allows the BCG Hubble diagram to be calibrated independently of recession velocities, yielding a far field H_0 = 89 +/- 10 km/s/Mpc with an effective depth of ~11,000 km/s. The error in this case is dominated by the photometric cosmic scatter in using BCG as distance estimators. The concordance of the present results with other recent H_0 determinations, and a review of theoretical treatments on perturbations in the near field Hubble flow, argue that going to the far field removes an important source of uncertainty, but that there is not a large systematic error to be corrected for to begin with. Further improvements in H_0 depend more on understanding nearby calibrators than on improved sampling of the distant flow.

We model the surface brightness profiles of a sample of 119 Abell Brightest Cluster Galaxies (BCG), finding a generalised deVaucouleurs R^{1/n} law, where n is a free parameter, to be appropriate. Departures from the R^{1/4} law are shown to be a real feature of galaxy profiles, not due to observational errors or coupling of n with the other model parameters. BCG typically have values of n greater than 4. The shape parameter n is shown to correlate with effective half-light radius, such that the larger BCG have larger values of n. This continues a trend noticed amongst ordinary elliptical galaxies and dwarf ellipticals, such that the fainter galaxies have smaller values of n.

We present a study of the surface density profiles of the clusters of galaxies from the Palomar Distant Cluster Survey (Postman et al. 1996). The survey contains a total of 79 clusters of galaxies, covering the estimated redshift range of $0.2 \simless z \simless 1.2$. We have analyzed the richest clusters in this sample and find that the typical Palomar cluster has a surface density profile of $r^{-1.4}$ ($r \ge 0.10~h^{-1}~{\rm Mpc}$) and a core radius of $0.05~h^{-1}~{\rm Mpc}$. There may be an indication that the slope of the surface density profile steepens with increasing redshift, though the observational uncertainty is at present too large to be conclusive. Our cluster population is inconsistent at a 99.9\% confidence level with a population of azimuthally symmetric clusters.

We present an optical/near IR selected catalog of 79 clusters distributed over an area of 5.1~square~degrees. The catalog was constructed from images obtained with the 4-Shooter CCD mosaic camera on the Hale 5m telescope operating in ``scan" mode. The survey, hereafter known as the Palomar Distant Cluster Survey (PDCS), was conducted in two broad band filters that closely resemble~$V$ and~$I$. The $4\sigma$ limiting magnitudes for our 300 s exposures are~$\sim$~23.8~($V$) and~22.5~($I$). A matched filter algorithm was developed and employed to identify the cluster candidates by using positional and photometric data simultaneously. The clusters cover the range $0.2 \simless z \simless 1.2$, based on the redshift estimates derived in the cluster detection procedure. An accurate selection function is generated from extensive simulations. We find that the cumulative surface density of clusters with richness class R $\ge$ 1 is about a factor of 5 higher than the extrapolated density of R $\ge$ 1 Abell clusters. The PDCS results are consistent with a constant comoving density of clusters to $z \simless 0.6$, albeit at the above high density level. Constraints on cluster abundances at $z > 0.6$ will be possible with the acquisition of spectroscopic redshifts for a large subset of these cluster candidates. We also present a supplemental catalog of 28 clusters that do not satisfy all our selection criteria but which include some of the most distant systems detected in the survey. Finding charts for all 107 cluster candidates are provided.

Lauer \& Postman (LP) observe that all Abell clusters with redshifts less than 15,000\kms\ appear to be participating in a bulk flow of 689 km s$^{-1}$ with respect to the Cosmic Microwave Background. We find this result difficult to reconcile with all popular models for large-scale structure formation that assume Gaussian initial conditions. This conclusion is based on Monte-Carlo realizations of the LP data, drawn from large Particle-Mesh $N$-body simulations. We have taken special care to treat properly the longest-wavelength components of the power spectra. Bulk flows with amplitude as large as that reported by LP are not uncommon in the Monte-Carlo datasets. However, the $\chi^2$ of the observed bulk flow, taking into account the anisotropy of the error ellipsoid, is much more difficult to match in the simulations. The models examined are ruled out at confidence levels between 94\% and 98\%. Any model that has {\it intrinsic} flows of less than 480\kms\ on the scales probed by LP scales can be ruled out at a similar level.