Data Challenge 1 (DC1) is the first synthetic dataset produced by the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC). DC1 is designed to develop and validate data reduction and analysis and to study the impact of systematic effects that will affect the LSST dataset. DC1 is comprised of $r$-band observations of 40 deg$^{2}$ to 10-year LSST depth. We present each stage of the simulation and analysis process: a) generation, by synthesizing sources from cosmological N-body simulations in individual sensor-visit images with different observing conditions; b) reduction using a development version of the LSST Science Pipelines; and c) matching to the input cosmological catalog for validation and testing. We verify that testable LSST requirements pass within the fidelity of DC1. We establish a selection procedure that produces a sufficiently clean extragalactic sample for clustering analyses and we discuss residual sample contamination, including contributions from inefficiency in star-galaxy separation and imperfect deblending. We compute the galaxy power spectrum on the simulated field and conclude that: i) survey properties have an impact of 50\% of the statistical uncertainty for the scales and models used in DC1 ii) a selection to eliminate artifacts in the catalogs is necessary to avoid biases in the measured clustering; iii) the presence of bright objects has a significant impact (2- to 6-$\sigma$) in the estimated power spectra at small scales ($\ell > 1200$), highlighting the impact of blending in studies at small angular scales in LSST;

The LSST survey will provide unprecedented statistical power for measurements of dark energy. Consequently, controlling systematic uncertainties is becoming more important than ever. The LSST observing strategy will affect the statistical uncertainty and systematics control for many science cases; here, we focus on weak lensing systematics. The fact that the LSST observing strategy involves hundreds of visits to the same sky area provides new opportunities for systematics mitigation. We explore these opportunities by testing how different dithering strategies (pointing offsets and rotational angle of the camera in different exposures) affect additive weak lensing shear systematics on a baseline operational simulation, using the $\rho-$statistics formalism. Some dithering strategies improve systematics control at the end of the survey by a factor of up to $\sim 3-4$ better than others. We find that a random translational dithering strategy, applied with random rotational dithering at every filter change, is the most effective of those strategies tested in this work at averaging down systematics. Adopting this dithering algorithm, we explore the effect of varying the area of the survey footprint, exposure time, number of exposures in a visit, and exposure to the Galactic plane. We find that any change that increases the average number of exposures (in filters relevant to weak lensing) reduces the additive shear systematics. Some ways to achieve this increase may not be favorable for the weak lensing statistical constraining power or for other probes, and we explore the relative trade-offs between these options given constraints on the overall survey parameters.

We analyze the clustering of galaxies in the first public data release of the HSC Subaru Strategic Program. Despite the relatively small footprints of the observed fields, the data are an excellent proxy for the deep photometric datasets that will be acquired by LSST, and are therefore an ideal test bed for the analysis methods being implemented by the LSST DESC. We select a magnitude limited sample with $i<24.5$ and analyze it in four redshift bins covering $0.15\lesssim z \lesssim1.5$. We carry out a Fourier-space analysis of the two-point clustering of this sample, including all auto- and cross-correlations. We demonstrate the use of map-level deprojection methods to account for fluctuations in the galaxy number density caused by observational systematics. Through an HOD analysis, we place constraints on the characteristic halo masses of this sample, finding a good fit up to scales $k_{\rm max}=1\,{\rm Mpc}^{-1}$, including both auto- and cross-correlations. Our results show monotonically decreasing average halo masses, which can be interpreted in terms of the drop-out of red galaxies at high redshifts for a flux-limited sample. In terms of photometric redshift systematics, we show that additional care is needed in order to marginalize over uncertainties in the redshift distribution in galaxy clustering, and that these uncertainties can be constrained by including cross-correlations. We are able to make a $\sim3\sigma$ detection of lensing magnification in the HSC data. Our results are stable to variations in $\sigma_8$ and $\Omega_c$ and we find constraints that agree well with measurements from Planck and low-redshift probes. Finally, we use our pipeline to study the clustering of galaxies as a function of limiting flux, and provide a simple fitting function for the linear galaxy bias for magnitude limited samples as a function of limiting magnitude and redshift. [abridged]

With the advent of surveys containing millions to billions of galaxies, it is imperative to develop analysis techniques that utilize the available statistical power. In galaxy clustering, even small sample contamination arising from distance uncertainties can lead to large artifacts, which the standard estimator does not account for. We first introduce a formalism, termed decontamination, that corrects for sample contamination by utilizing the observed cross-correlations in the contaminated samples; this corrects any correlation function estimator for contamination. Using this formalism, we present a new estimator that uses the standard estimator to measure correlation functions in the contaminated samples but then corrects for contamination. We also introduce a weighted estimator that assigns each galaxy a weight in each redshift bin based on its probability of being in that bin. We demonstrate that these estimators effectively recover the true correlation functions and their covariance matrices. Our estimators can correct for sample contamination caused by misclassification between object types as well as photometric redshifts; they should be particularly helpful for studies of galaxy evolution and baryonic acoustic oscillations, where forward-modeling the clustering signal using the contaminated redshift distribution is undesirable.

We present a new method for inferring galaxy star formation histories (SFH) using machine learning methods coupled with two cosmological hydrodynamic simulations. We train Convolutional Neural Networks to learn the relationship between synthetic galaxy spectra and high resolution SFHs from the EAGLE and Illustris models. To evaluate our SFH reconstruction we use Symmetric Mean Absolute Percentage Error (SMAPE), which acts as a true percentage error in the low-error regime. On dust-attenuated spectra we achieve high test accuracy (median SMAPE $= 10.5\%$). Including the effects of simulated observational noise increases the error ($12.5\%$), however this is alleviated by including multiple realisations of the noise, which increases the training set size and reduces overfitting ($10.9\%$). We also make estimates for the observational and modelling errors. To further evaluate the generalisation properties we apply models trained on one simulation to spectra from the other, which leads to only a small increase in the error (median SMAPE $\sim 15\%$). We apply each trained model to SDSS DR7 spectra, and find smoother histories than in the VESPA catalogue. This new approach complements the results of existing SED fitting techniques, providing star formation histories directly motivated by the results of the latest cosmological simulations.

We present a WFC3 F160W ($H$-band) selected catalog in the CANDELS/GOODS-N field containing photometry from the ultraviolet (UV) to the far-infrared (IR), photometric redshifts and stellar parameters derived from the analysis of the multi-wavelength data. The catalog contains 35,445 sources over the 171 arcmin$^{2}$ of the CANDELS F160W mosaic. The 5$\sigma$ detection limits (within an aperture of radius 0\farcs17) of the mosaic range between $H=27.8$, 28.2 and 28.7 in the wide, intermediate and deep regions, that span approximately 50\%, 15\% and 35\% of the total area. The multi-wavelength photometry includes broad-band data from UV (U band from KPNO and LBC), optical (HST/ACS F435W, F606W, F775W, F814W, and F850LP), near-to-mid IR (HST/WFC3 F105W, F125W, F140W and F160W, Subaru/MOIRCS Ks, CFHT/Megacam K, and \spitzer/IRAC 3.6, 4.5, 5.8, 8.0 $\mu$m) and far IR (\spitzer/MIPS 24$\mu$m, HERSCHEL/PACS 100 and 160$\mu$m, SPIRE 250, 350 and 500$\mu$m) observations. In addition, the catalog also includes optical medium-band data (R$\sim50$) in 25 consecutive bands, $\lambda=500$ to 950~nm, from the SHARDS survey and WFC3 IR spectroscopic observations with the G102 and G141 grisms (R$\sim210$ and 130). The use of higher spectral resolution data to estimate photometric redshifts provides very high, and nearly uniform, precision from $z=0-2.5$. The comparison to 1,485 good quality spectroscopic redshifts up to $z\sim3$ yields $\Delta z$/(1+$z_{\rm spec}$)$=$0.0032 and an outlier fraction of $\eta=$4.3\%. In addition to the multi-band photometry, we release added-value catalogs with emission line fluxes, stellar masses, dust attenuations, UV- and IR- based star formation rates and rest-frame colors.

This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `Outer Rim' run. It covers 440 deg^2 of sky area to a redshift of z=3 and is complete to a magnitude depth of 28 in the r-band. Each galaxy is characterized by a multitude of properties including stellar mass, morphology, spectral energy distributions, broadband filter magnitudes, host halo information and weak lensing shear. The size and complexity of cosmoDC2 requires an efficient catalog generation methodology; our approach is based on a new hybrid technique that combines data-driven empirical approaches with semi-analytic galaxy modeling. A wide range of observation-based validation tests has been implemented to ensure that cosmoDC2 enables the science goals of the planned LSST DESC DC2 analyses. This paper also represents the official release of the cosmoDC2 data set, including an efficient reader that facilitates interaction with the data.

The Large Synoptic Survey Telescope (LSST) can advance scientific frontiers beyond its groundbreaking 10-year survey. Here we explore opportunities for extended operations with proposal-based observing strategies, new filters, or transformed instrumentation. We recommend the development of a mid-decade community- and science-driven process to define next-generation LSST capabilities.

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

One of the key questions in understanding the formation and evolution of galaxies is how starbursts affect the assembly of stellar populations in galaxies over time. We define a burst indicator ($\eta$), which compares a galaxy's star formation rates on short ($\sim10$ Myr) and long ($\sim100$ Myr) timescales. To estimate $\eta$, we apply the detailed time-luminosity relationship for H$\alpha$ and near-ultraviolet emission to simulated star formation histories (SFHs) from semi-analytic models and the Mufasa hydrodynamical cosmological simulations. The average of $\eta$ is not a good indicator of star formation stochasticity (burstiness); indeed, we show that this average should be close to zero unless the population has an average SFH which is rising or falling rapidly. Instead, the width of the $\eta$ distribution characterizes the burstiness of a galaxy population's recent star formation. We find this width to be robust to variations in stellar initial mass function and metallicity. We apply realistic noise and selection effects to the models to generate mock HST and JWST galaxy catalogs and compare these catalogs with 3D-HST observations of 956 galaxies at $0.65<z<1.5$ detected in H$\alpha$. Measurements of $\eta$ are unaffected by dust measurement errors under the assumption that $E(B-V)_\mathrm{stars}=0.44\,E(B-V)_\mathrm{gas}$ (i.e., $Q_\mathrm{sg}=0.44$). However, setting $Q_\mathrm{sg}=0.8^{+0.1}_{-0.2}$ removes an unexpected dependence of the average value of $\eta$ upon dust attenuation and stellar mass in the 3D-HST sample while also resolving disagreements in the distribution of star formation rates. However, even varying the dust law cannot resolve all discrepancies between the simulated and the observed galaxies.

Community access to deep (i ~ 25), highly-multiplexed optical and near-infrared multi-object spectroscopy (MOS) on 8-40m telescopes would greatly improve measurements of cosmological parameters from LSST. The largest gain would come from improvements to LSST photometric redshifts, which are employed directly or indirectly for every major LSST cosmological probe; deep spectroscopic datasets will enable reduced uncertainties in the redshifts of individual objects via optimized training. Such spectroscopy will also determine the relationship of galaxy SEDs to their environments, key observables for studies of galaxy evolution. The resulting data will also constrain the impact of blending on photo-z's. Focused spectroscopic campaigns can also improve weak lensing cosmology by constraining the intrinsic alignments between the orientations of galaxies. Galaxy cluster studies can be enhanced by measuring motions of galaxies in and around clusters and by testing photo-z performance in regions of high density. Photometric redshift and intrinsic alignment studies are best-suited to instruments on large-aperture telescopes with wider fields of view (e.g., Subaru/PFS, MSE, or GMT/MANIFEST) but cluster investigations can be pursued with smaller-field instruments (e.g., Gemini/GMOS, Keck/DEIMOS, or TMT/WFOS), so deep MOS work can be distributed amongst a variety of telescopes. However, community access to large amounts of nights for surveys will still be needed to accomplish this work. In two companion white papers we present gains from shallower, wide-area MOS and from single-target imaging and spectroscopy.

LSST will open new vistas for cosmology in the next decade, but it cannot reach its full potential without data from other telescopes. Cosmological constraints can be greatly enhanced using wide-field ($>20$ deg$^2$ total survey area), highly-multiplexed optical and near-infrared multi-object spectroscopy (MOS) on 4-15m telescopes. This could come in the form of suitably-designed large surveys and/or community access to add new targets to existing projects. First, photometric redshifts can be calibrated with high precision using cross-correlations of photometric samples against spectroscopic samples at $0 < z < 3$ that span thousands of sq. deg. Cross-correlations of faint LSST objects and lensing maps with these spectroscopic samples can also improve weak lensing cosmology by constraining intrinsic alignment systematics, and will also provide new tests of modified gravity theories. Large samples of LSST strong lens systems and supernovae can be studied most efficiently by piggybacking on spectroscopic surveys covering as much of the LSST extragalactic footprint as possible (up to $\sim20,000$ square degrees). Finally, redshifts can be measured efficiently for a high fraction of the supernovae in the LSST Deep Drilling Fields (DDFs) by targeting their hosts with wide-field spectrographs. Targeting distant galaxies, supernovae, and strong lens systems over wide areas in extended surveys with (e.g.) DESI or MSE in the northern portion of the LSST footprint or 4MOST in the south could realize many of these gains; DESI, 4MOST, Subaru/PFS, or MSE would all be well-suited for DDF surveys. The most efficient solution would be a new wide-field, highly-multiplexed spectroscopic instrument in the southern hemisphere with $>6$m aperture. In two companion white papers we present gains from deep, small-area MOS and from single-target imaging and spectroscopy.

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

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

The star formation histories (SFHs) of galaxies contain imprints of the physical processes responsible for regulating star formation during galaxy growth and quenching. We improve the Dense Basis SFH reconstruction method of Iyer & Gawiser (2017), introducing a nonparametric description of the SFH based on the lookback times at which a galaxy assembles certain quantiles of its stellar mass. The method uses Gaussian Processes to create smooth SFHs that are independent of any functional form, with a flexible number of parameters that is adjusted to extract the maximum possible amount of SFH information from the SEDs being fit. We apply the method to reconstruct the SFHs of 48,791 galaxies with $H<25$ at $0.5 < z < 3.0$ across the five CANDELS fields. Using these SFHs, we study the evolution of galaxies as they grow more massive over cosmic time. We quantify the fraction of galaxies that show multiple major episodes of star formation, finding that the median time between two peaks of star formation is $\sim 0.42_{-0.10}^{+0.15}t_{univ}$ Gyr, where $t_{univ}$ is the age of the universe at a given redshift and remains roughly constant with stellar mass. Correlating SFHs with morphology, we find that studying the median SFHs of galaxies at $0.5<z<1.0$ at the same mass ($10^{10}< M_* < 10^{10.5}M_\odot$) allows us to compare the timescales on which the SFHs decline for different morphological classifications, ranging from $0.60^{-0.54}_{+1.54}$ Gyr for galaxies with spiral arms to $2.50^{-1.50}_{+2.25}$ Gyr for spheroids. The Gaussian Process-based SFH description provides a general approach to reconstruct smooth, nonparametric SFH posteriors for galaxies with a flexible number of parameters that can be incorporated into Bayesian SED fitting codes to minimize the bias in estimating physical parameters due to SFH parametrization.

Cosmology is one of the four science pillars of LSST, which promises to be transformative for our understanding of dark energy and dark matter. The LSST Dark Energy Science Collaboration (DESC) has been tasked with deriving constraints on cosmological parameters from LSST data. Each of the cosmological probes for LSST is heavily impacted by the choice of observing strategy. This white paper is written by the LSST DESC Observing Strategy Task Force (OSTF), which represents the entire collaboration, and aims to make recommendations on observing strategy that will benefit all cosmological analyses with LSST. It is accompanied by the DESC DDF (Deep Drilling Fields) white paper (Scolnic et al.). We use a variety of metrics to understand the effects of the observing strategy on measurements of weak lensing, large-scale structure, clusters, photometric redshifts, supernovae, strong lensing and kilonovae. In order to reduce systematic uncertainties, we conclude that the current baseline observing strategy needs to be significantly modified to result in the best possible cosmological constraints. We provide some key recommendations: moving the WFD (Wide-Fast-Deep) footprint to avoid regions of high extinction, taking visit pairs in different filters, changing the 2x15s snaps to a single exposure to improve efficiency, focusing on strategies that reduce long gaps (>15 days) between observations, and prioritizing spatial uniformity at several intervals during the 10-year survey.

We present the $ugriz$-band Dark Energy Camera (DECam) plus 3.6 and 4.5 $\mu$m IRAC catalogs for the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers $\sim24$ deg$^{2}$ of the Sloan Digital Sky Survey (SDSS) Stripe 82 region, with seven bandpasses spanning a wavelength range of 0.35 to 4.5 $\mu$m. SHELA falls within the footprint of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), which will provide spectroscopic redshifts for $\sim200{,}000$ Ly$\alpha$ emitters at $1.9<z<3.5$ and also for $\sim200{,}000$ [OII] emitters at $z<0.5$. SHELA's deep, wide-area multiwavelength images combined with HETDEX's spectroscopic information, will facilitate many extragalactic studies, including measuring the evolution of galaxy stellar mass, halo mass, and environment from $1.5<z<3.5$. Here we present $riz$-band selected $ugriz$-band DECam catalogs that reach a $5\sigma$ depth of $\sim24.5$ AB mag (for point sources with an aperture that encloses $70\%$ of the total flux) and cover $17.5$ deg$^{2}$ of the overall SHELA field. We validate our DECam catalog by comparison to the DECam Legacy Survey (DECaLS) DR5 and the Dark Energy Survey (DES) DR1. We perform IRAC forced photometry with The Tractor image modeling code to measure 3.6 and 4.5 $\mu$m fluxes for all objects within our DECam catalog. We demonstrate the utility of our catalog by computing galaxy number counts and estimating photometric redshifts. Our photometric redshifts recover the available $\left\langle z \right\rangle = 0.33 $ SDSS spectroscopic redshifts with a $1\sigma$ scatter in $\Delta z/(1 +z)$ of 0.04.

The LSST survey was designed to deliver transformative results for four primary objectives: constraining dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. While the LSST Wide-Fast-Deep survey and accompanying Deep Drilling and mini-surveys will be ground-breaking for each of these areas, there remain competing demands on the survey area, depth, and temporal coverage amid a desire to maximize all three. In this white paper, we seek to address a principal source of tension between the different LSST science collaborations, that of the survey area and depth that they each need in the parts of the sky that they care about. We present simple tools which can be used to explore trades between the area surveyed by LSST and the number of visits available per field and then use these tools to propose a change to the baseline survey strategy. Specifically, we propose to reconfigure the WFD footprint to consist of low-extinction regions (limited by galactic latitude), with the number of visits per field in WFD limited by the LSST Science Requirements Document (SRD) design goal, and suggest assignment of the remaining LSST visits to the full visible LSST sky. This proposal addresses concerns with the WFD footprint raised by the DESC (as 25 percent of the current baseline WFD region is not usable for dark energy science due to MW dust extinction), eases the time required for the NES and SCP mini-surveys (since in our proposal they would partially fall into the modified WFD footprint), raises the number of visits previously assigned to the GP region, and increases the overlap with DESI and other Northern hemisphere follow-up facilities. This proposal alleviates many of the current concerns of Science Collaborations that represent the four scientific pillars of LSST and provides a Big Sky approach to cadence diplomacy.

We review the measurements of dark energy enabled by observations of the Deep Drilling Fields and the optimization of survey design for cosmological measurements. This white paper is the result of efforts by the LSST DESC Observing Strategy Task Force (OSTF), which represents the entire collaboration, and aims to make recommendations on observing strategy for the DDFs that will benefit all cosmological analyses with LSST. It is accompanied by the DESC-WFD white paper (Lochner et al.). We argue for altering the nominal deep drilling plan to have $>6$ month seasons, interweaving $gri$ and $zy$ observations every 3 days with 2, 4, 8, 25, 4 visits in $grizy$, respectively. These recommendations are guided by metrics optimizing constraints on dark energy and mitigation of systematic uncertainties, including specific requirements on total number of visits after Y1 and Y10 for photometric redshifts (photo-$z$) and weak lensing systematics. We specify the precise locations for the previously-chosen LSST deep fields (ELAIS-S1, XMM-LSS, CDF-S, and COSMOS) and recommend Akari Deep Field South as the planned fifth deep field in order to synergize with Euclid and WFIRST. Our recommended DDF strategy uses $6.2\%$ of the LSST survey time. We briefly discuss synergy with white papers from other collaborations, as well as additional mini-surveys and Target-of-Opportunity programs that lead to better measurements of dark energy.

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

To achieve a fuller understanding of galaxy evolution, SED fitting can be used to recover quantities beyond stellar masses (M$_*$) and star formation rates (SFRs). We use Star Formation Histories (SFHs) reconstructed via the Dense Basis method of Iyer \& Gawiser (2017) for a sample of $17,873$ galaxies at $0.5<z<6$ in the CANDELS GOODS-S field to study the nature and evolution of the SFR-M$_*$ correlation. The reconstructed SFHs represent trajectories in SFR-M$_*$ space, enabling us to study galaxies at epochs earlier than observed by propagating them backwards in time along these trajectories. We study the SFR-M$_*$ correlation at $z=1,2,3,4,5,6$ using both direct fits to galaxies observed at those epochs and SFR-M$_*$ trajectories of galaxies observed at lower redshifts. The SFR-M$_*$ correlations obtained using the two approaches are found to be consistent with each other through a KS test. Validation tests using SFHs from semi-analytic models and cosmological hydrodynamical simulations confirm the sensitivity of the method to changes in the slope, normalization and shape of the SFR-M$_*$ correlation. This technique allows us to further probe the low-mass regime of the correlation at high-z by $\sim 1$ dex and over an effective volume of $\sim 10\times$ larger than possible with just direct fits. We find that the SFR-M$_*$ correlation is consistent with being linear down to M$_*\sim 10^7 M_\odot$ at $z>4$. The evolution of the correlation is well described by $\log SFR= (0.80\pm 0.029 - 0.017\pm 0.010\times t_{univ})\log M_*$ $- (6.487\pm 0.282-0.039\pm 0.008\times t_{univ})$, where $t_{univ}$ is the age of the universe in Gyr.

The Large Synoptic Survey Telescope (LSST) Dark Energy Science Collaboration (DESC) will use five cosmological probes: galaxy clusters, large scale structure, supernovae, strong lensing, and weak lensing. This Science Requirements Document (SRD) quantifies the expected dark energy constraining power of these probes individually and together, with conservative assumptions about analysis methodology and follow-up observational resources based on our current understanding and the expected evolution within the field in the coming years. We then define requirements on analysis pipelines that will enable us to achieve our goal of carrying out a dark energy analysis consistent with the Dark Energy Task Force definition of a Stage IV dark energy experiment. This is achieved through a forecasting process that incorporates the flowdown to detailed requirements on multiple sources of systematic uncertainty. Future versions of this document will include evolution in our software capabilities and analysis plans along with updates to the LSST survey strategy.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

This is the first in a series of papers examining the demographics of star-forming galaxies at $0.2<z<2.5$ in CANDELS. We study 9,100 galaxies from GOODS-S and UDS having published values of redshifts, masses, star-formation rates (SFRs), and dust attenuation ($A_V$) derived from UV-optical SED fitting. In agreement with previous works, we find that the $UVJ$ colors of a galaxy are closely correlated with its specific star-formation rate (SSFR) and $A_V$. We define rotated $UVJ$ coordinate axes, termed $S_\mathrm{SED}$ and $C_\mathrm{SED}$, that are parallel and perpendicular to the star-forming sequence and derive a quantitative calibration that predicts SSFR from $C_\mathrm{SED}$ with an accuracy of ~0.2 dex. SFRs from UV-optical fitting and from UV+IR values based on Spitzer/MIPS 24 $\mu\mathrm{m}$ agree well overall, but systematic differences of order 0.2 dex exist at high and low redshifts. A novel plotting scheme conveys the evolution of multiple galaxy properties simultaneously, and dust growth, as well as star-formation decline and quenching, exhibit "mass-accelerated evolution" ("downsizing"). A population of transition galaxies below the star-forming main sequence is identified. These objects are located between star-forming and quiescent galaxies in $UVJ$ space and have lower $A_V$ and smaller radii than galaxies on the main sequence. Their properties are consistent with their being in transit between the two regions. The relative numbers of quenched, transition, and star-forming galaxies are given as a function of mass and redshift.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

The Large Synoptic Survey Telescope (LSST) will enable revolutionary studies of galaxies, dark matter, and black holes over cosmic time. The LSST Galaxies Science Collaboration has identified a host of preparatory research tasks required to leverage fully the LSST dataset for extragalactic science beyond the study of dark energy. This Galaxies Science Roadmap provides a brief introduction to critical extragalactic science to be conducted ahead of LSST operations, and a detailed list of preparatory science tasks including the motivation, activities, and deliverables associated with each. The Galaxies Science Roadmap will serve as a guiding document for researchers interested in conducting extragalactic science in anticipation of the forthcoming LSST era.

We present spectroscopic redshifts of S(870)>2mJy submillimetre galaxies (SMGs) which have been identified from the ALMA follow-up observations of 870um detected sources in the Extended Chandra Deep Field South (the ALMA-LESS survey). We derive spectroscopic redshifts for 52 SMGs, with a median of z=2.4+/-0.1. However, the distribution features a high redshift tail, with ~25% of the SMGs at z>3. Spectral diagnostics suggest that the SMGs are young starbursts, and the velocity offsets between the nebular emission and UV ISM absorption lines suggest that many are driving winds, with velocity offsets up to 2000km/s. Using the spectroscopic redshifts and the extensive UV-to-radio photometry in this field, we produce optimised spectral energy distributions (SEDs) using Magphys, and use the SEDs to infer a median stellar mass of M*=(6+/-1)x10^{10}Msol for our SMGs with spectroscopic redshifts. By combining these stellar masses with the star-formation rates (measured from the far-infrared SEDs), we show that SMGs (on average) lie a factor ~5 above the main-sequence at z~2. We provide this library of 52 template fits with robust and well-sampled SEDs available as a resource for future studies of SMGs, and also release the spectroscopic catalog of ~2000 (mostly infrared-selected) galaxies targeted as part of the spectroscopic campaign.

We present a 0.4-8$\mu$m multi-wavelength photometric catalog in the Extended Groth Strip (EGS) field. This catalog is built on the Hubble Space Telescope (HST) WFC3 and ACS data from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS), and it incorporates the existing HST data from the All-wavelength Extended Groth strip International Survey (AEGIS) and the 3D-HST program. The catalog is based on detections in the F160W band reaching a depth of F160W=26.62 AB (90% completeness, point-sources). It includes the photometry for 41457 objects over an area of $\approx 206$ arcmin$^2$ in the following bands: HST ACS F606W and F814W; HST WFC3 F125W, F140W and F160W; CFHT/Megacam $u^*$, $g'$, $r'$, $i'$ and $z'$; CFHT/WIRCAM $J$, $H$ and $K_\mathrm{S}$; Mayall/NEWFIRM $J1$, $J2$, $J3$, $H1$, $H2$, $K$; Spitzer IRAC $3.6\mu$m, $4.5\mu$m, $5.8\mu$m and $8.0\mu$m. We are also releasing value-added catalogs that provide robust photometric redshifts and stellar mass measurements. The catalogs are publicly available through the CANDELS repository.

We introduce the Dense Basis method for Spectral Energy Distribution (SED) fitting. It accurately recovers traditional SED parameters, including M$_*$, SFR and dust attenuation, and reveals previously inaccessible information about the number and duration of star formation episodes and the timing of stellar mass assembly, as well as uncertainties in these quantities. This is done using basis Star Formation Histories (SFHs) chosen by comparing the goodness-of-fit of mock galaxy SEDs to the goodness-of-reconstruction of their SFHs. We train and validate the method using a sample of realistic SFHs at $z =1$ drawn from stochastic realisations, semi-analytic models, and a cosmological hydrodynamical galaxy formation simulation. The method is then applied to a sample of 1100 CANDELS GOODS-S galaxies at $1<z<1.5$ to illustrate its capabilities at moderate S/N with 15 photometric bands. Of the six parametrizations of SFHs considered, we adopt linear-exponential, bessel-exponential, lognormal and gaussian SFHs and reject the traditional parametrizations of constant (Top-Hat) and exponential SFHs. We quantify the bias and scatter of each parametrization. $15\%$ of galaxies in our CANDELS sample exhibit multiple episodes of star formation, with this fraction decreasing above $M_*>10^{9.5}M_\odot$. About $40\%$ of the CANDELS galaxies have SFHs whose maximum occurs at or near the epoch of observation. The Dense Basis method is scalable and offers a general approach to a broad class of data-science problems.

We present an investigation of clumpy galaxies in the Hubble Ultra Deep Field at 0.5 $\leq z \leq$ 1.5 in the rest-frame far-ultraviolet (FUV) using HST WFC3 broadband imaging in F225W, F275W, and F336W. An analysis of 1,404 galaxies yields 209 galaxies that host 403 kpc-scale clumps. These host galaxies appear to be typical star-forming galaxies, with an average of 2 clumps per galaxy and reaching a maximum of 8 clumps. We measure the photometry of the clumps, and determine the mass, age, and star formation rates (SFR) utilizing the SED-fitting code FAST. We find that clumps make an average contribution of 19% to the total rest-frame FUV flux of their host galaxy. Individually, clumps contribute a median of 5% to the host galaxy SFR and an average of $\sim$4% to the host galaxy mass, with total clump contributions to the host galaxy stellar mass ranging widely from less than 1% up to 93%. Clumps in the outskirts of galaxies are typically younger, with higher star formation rates, than clumps in the inner regions. The results are consistent with clump migration theories in which clumps form through violent gravitational instabilities in gas-rich turbulent disks, eventually migrate toward the center of the galaxies, and coalesce into the bulge.

We investigate the burstiness of star formation histories (SFHs) of galaxies at $0.4<z<1$ by using the ratio of star formation rates (SFRs) measured from H$\beta$ and FUV (1500 \AA) (H$\beta$--to--FUV ratio). Our sample contains 164 galaxies down to stellar mass (M*) of $10^{8.5} M_\odot$ in the CANDELS GOODS-N region, where Team Keck Redshift Survey DEIMOS spectroscopy and HST/WFC3 F275W images from CANDELS and Hubble Deep UV Legacy Survey are available. When the {\it ratio} of H$\beta$- and FUV-derived SFRs is measured, dust extinction correction is negligible (except for very dusty galaxies) with the Calzetti attenuation curve. The H$\beta$--to--FUV ratio of our sample increases with M* and SFR. The median ratio is $\sim$0.7 at M*$\sim10^{8.5} M_\odot$ (or SFR$\sim 0.5 M_\odot/yr$) and increases to $\sim$1 at M*$\sim10^{10} M_\odot$ (or SFR $\sim 10 M_\odot/yr$). At M*$<10^{9.5} M_\odot$, our median H$\beta$--to--FUV ratio is lower than that of local galaxies at the same M*, implying a redshift evolution. Bursty SFH on a timescale of a few tens of megayears on galactic scales provides a plausible explanation of our results, and the importance of the burstiness increases as M* decreases. Due to sample selection effects, our H$\beta$--to--FUV ratio may be an upper limit of the true value of a complete sample, which strengthens our conclusions. Other models, e.g., non-universal initial mass function or stochastic star formation on star cluster scales, are unable to plausibly explain our results.

The Large Synoptic Survey Telescope (LSST) will be a discovery machine for the astronomy and physics communities, revealing astrophysical phenomena from the Solar System to the outer reaches of the observable Universe. While many discoveries will be made using LSST data alone, taking full scientific advantage of LSST will require ground-based optical-infrared (OIR) supporting capabilities, e.g., observing time on telescopes, instrumentation, computing resources, and other infrastructure. This community-based study identifies, from a science-driven perspective, capabilities that are needed to maximize LSST science. Expanding on the initial steps taken in the 2015 OIR System Report, the study takes a detailed, quantitative look at the capabilities needed to accomplish six representative LSST-enabled science programs that connect closely with scientific priorities from the 2010 decadal surveys. The study prioritizes the resources needed to accomplish the science programs and highlights ways that existing, planned, and future resources could be positioned to accomplish the science goals.

We compare the rest-frame ultraviolet and rest-frame optical morphologies of 2 < z < 3 star-forming galaxies in the GOODS-S field using Hubble Space Telescope WFC3 and ACS images from the CANDELS, GOODS, and ERS programs. We show that the distribution of sizes and concentrations for 1.90 < z < 2.35 galaxies selected via their rest-frame optical emission-lines are statistically indistinguishable from those of Lyman-alpha emitting systems found at z ~ 2.1 and z ~ 3.1. We also show that the z > 2 star-forming systems of all sizes and masses become smaller and more compact as one shifts the observing window from the UV to the optical. We argue that this offset is due to inside-out galaxy formation over the first ~ 2 Gyr of cosmic time.

Although there has been much progress in understanding how galaxies evolve, we still do not understand how and when they stop forming stars and become quiescent. We address this by applying our galaxy spectral energy distribution models, which incorporate physically motivated star formation histories (SFHs) from cosmological simulations, to a sample of quiescent galaxies at $0.2<z<2.1$. A total of 845 quiescent galaxies with multi-band photometry spanning rest-frame ultraviolet through near-infrared wavelengths are selected from the CANDELS dataset. We compute median SFHs of these galaxies in bins of stellar mass and redshift. At all redshifts and stellar masses, the median SFHs rise, reach a peak, and then decline to reach quiescence. At high redshift, we find that the rise and decline are fast, as expected because the Universe is young. At low redshift, the duration of these phases depends strongly on stellar mass. Low-mass galaxies ($\log(M_{\ast}/M_{\odot})\sim9.5$) grow on average slowly, take a long time to reach their peak of star formation ($\gtrsim 4$ Gyr), and the declining phase is fast ($\lesssim 2$ Gyr). Conversely, high-mass galaxies ($\log(M_{\ast}/M_{\odot})\sim11$) grow on average fast ($\lesssim 2$ Gyr), and, after reaching their peak, decrease the star formation slowly ($\gtrsim 3$ Gyr). These findings are consistent with galaxy stellar mass being a driving factor in determining how evolved galaxies are, with high-mass galaxies being the most evolved at any time (i.e., downsizing). The different durations we observe in the declining phases also suggest that low- and high-mass galaxies experience different quenching mechanisms that operate on different timescales.

The Large Synoptic Survey Telescope (LSST) will survey the southern sky from 2022--2032 with unprecedented detail. Since the observing strategy can lead to artifacts in the data, we investigate the effects of telescope-pointing offsets (called dithers) on the $r$-band coadded 5$\sigma$ depth yielded after the 10-year survey. We analyze this survey depth for several geometric patterns of dithers (e.g., random, hexagonal lattice, spiral) with amplitude as large as the radius of the LSST field-of-view, implemented on different timescales (per season, per night, per visit). Our results illustrate that per night and per visit dither assignments are more effective than per season. Also, we find that some dither geometries (e.g., hexagonal lattice) are particularly sensitive to the timescale on which the dithers are implemented, while others like random dithers perform well on all timescales. We then model the propagation of depth variations to artificial fluctuations in galaxy counts, which are a systematic for large-scale structure studies. We calculate the bias in galaxy counts caused by the observing strategy, accounting for photometric calibration uncertainties, dust extinction, and magnitude cuts; uncertainties in this bias limit our ability to account for structure induced by the observing strategy. We find that after 10 years of the LSST survey, the best dither strategies lead to uncertainties in this bias smaller than the minimum statistical floor for a galaxy catalog as deep as $r$$<$27.5. A few of these strategies bring the uncertainties close to the statistical floor for $r$$<$25.7 after only one year of survey.

We present near-infrared integral-field spectroscopic observations targeting H$\alpha$ in eight sub-millimeter galaxies (SMGs) at $z$=1.3-2.5 using VLT/SINFONI, obtaining significant detections for six of them. The star formation rates derived from the H$\alpha$ emission are $\sim$100 M$_\odot$yr$^{-1}$, which account for only $\sim$ 20-30\% of the infrared-derived values, thus suggesting that these systems are very dusty. Two of these systems present [NII]/H$\alpha$ ratios indicative of the presence of an Active Galactic Nucleus (AGN). We mapped the spatial distribution and kinematics of the star forming regions in these galaxies on kpc-scales. In general, the H$\alpha$ morphologies tend to be highly irregular and/or clumpy, showing spatial extents of $\sim$3-11~kpc. We find evidence for significant spatial offsets, of $\sim$0.1-0.4$"$ or 1.2-3.4 kpc, between the H$\alpha$ and the continuum emission in three of the sources. Performing a kinemetry analysis we conclude that the majority of the sample is not consistent with disk-like rotation-dominated kinematics. Instead, they tend to show irregular and/or clumpy and turbulent velocity and velocity dispersion fields. This can be interpreted as evidence for scenario in which these extreme star formation episodes are triggered by galaxy-galaxy interactions and major mergers. In contrast to recent results for SMGs, these sources appear to follow the same relations between gas and star forming rate densities as less luminous and/or normal star forming galaxies.

We present a Bayesian approach to the redshift classification of emission-line galaxies when only a single emission line is detected spectroscopically. We consider the case of surveys for high-redshift Lyman-alpha-emitting galaxies (LAEs), which have traditionally been classified via an inferred rest-frame equivalent width (EW) greater than 20 angstrom. Our Bayesian method relies on known prior probabilities in measured emission-line luminosity functions and equivalent width distributions for the galaxy populations, and returns the probability that an object in question is an LAE given the characteristics observed. This approach will be directly relevant for the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), which seeks to classify ~10^6 emission-line galaxies into LAEs and low-redshift [O II] emitters. For a simulated HETDEX catalog with realistic measurement noise, our Bayesian method recovers 86% of LAEs missed by the traditional EW > 20 angstrom cutoff over 2 < z < 3, outperforming the EW cut in both contamination and incompleteness. This is due to the method's ability to trade off between the two types of binary classification error by adjusting the stringency of the probability requirement for classifying an observed object as an LAE. In our simulations of HETDEX, this method reduces the uncertainty in cosmological distance measurements by 14% with respect to the EW cut, equivalent to recovering 29% more cosmological information. Rather than using binary object labels, this method enables the use of classification probabilities in large-scale structure analyses. It can be applied to narrowband emission-line surveys as well as upcoming large spectroscopic surveys including Euclid and WFIRST.

We present post-cryogenic Spitzer imaging at 3.6 and 4.5 micron with the Infrared Array Camera (IRAC) of the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers $\sim$deg$^2$ of the Sloan Digital Sky Survey "Stripe 82" region, and falls within the footprints of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) and the Dark Energy Survey. The HETDEX blind R $\sim$ 800 spectroscopy will produce $\sim$ 200,000 redshifts from the Lyman-$\alpha$ emission for galaxies in the range 1.9 < z < 3.5, and an additional $\sim$200,000 redshifts from the [OII] emission for galaxies at z < 0.5. When combined with deep ugriz images from the Dark Energy Camera, K-band images from NEWFIRM, and other ancillary data, the IRAC photometry from Spitzer will enable a broad range of scientific studies of the relationship between structure formation, galaxy stellar mass, halo mass, AGN, and environment over a co-moving volume of $\sim$0.5 Gpc$^3$ at 1.9 < z < 3.5. Here, we discuss the properties of the SHELA IRAC dataset, including the data acquisition, reduction, validation, and source catalogs. Our tests show the images and catalogs are 80% (50%) complete to limiting magnitudes of 22.0 (22.6) AB mag in the detection image, which is constructed from the weighted sum of the IRAC 3.6 and 4.5 micron images. The catalogs reach limiting sensitivities of 1.1 $\mu$Jy at both 3.6 and 4.5 micron (1$\sigma$, for R=2 arcsec circular apertures). As a demonstration of science, we present IRAC number counts, examples of highly temporally variable sources, and galaxy surface density profiles of rich galaxy clusters. In the spirit of Spitzer Exploratory programs we provide all images and catalogs as part of the publication.

We present estimates of intrinsic scatter in the Star Formation Rate (SFR) - Stellar Mass (M*) correlation in the redshift range 0.5 < z < 3.0 and in the mass range 10^7 < M* < 10^11 Msun. We utilize photometry in the Hubble Ultradeep Field (HUDF12), Ultraviolet Ultra Deep Field (UVUDF) campaigns and CANDELS/GOODS-S. We estimate SFR, M* from broadband Spectral Energy Distributions (SEDs) and the best available redshifts. The maximum depth of the HUDF photometry (F160W 29.9 AB, 5 sigma depth) probes the SFR-M* correlation down to M* ~ 10 ^7 Msun, a factor of 10-100X lower in M* than previous studies, and comparable to dwarf galaxies in the local universe. We find the slope of the SFR-M* relationship to be near unity at all redshifts and the normalization to decrease with cosmic time. We find a moderate increase in intrinsic scatter with cosmic time from 0.2 to 0.4 dex across the epoch of peak cosmic star formation. None of our redshift bins show a statistically significant increase in intrinsic scatter at low mass. However, it remains possible that intrinsic scatter increases at low mass on timescales shorter than ~ 100 Myr. Our results are consistent with a picture of gradual and self-similar assembly of galaxies across more than three orders of magnitude in stellar mass from as low as 10^7 Msun.

Aim: We aim to measure the LyC signal from a sample of sources in the Chandra deep field south. We collect star-forming galaxies (SFGs) and active galactic nuclei (AGN) with accurate spectroscopic redshifts, for which Hubble Space Telescope (HST) coverage and multi-wavelength photometry are available. Method: We selected a sample of about 200 sources at z~3. Taking advantage of HST resolution, we applied a careful cleaning procedure and rejected sources showing nearby clumps with different colours, which could be lower-z interlopers. Our clean sample consisted of 86 SFGs (including 19 narrow-band selected Lya emitters) and 8 AGN (including 6 detected in X-rays). We measured the LyC flux from aperture photometry in four narrow-band filters covering wavelengths below a 912 A rest frame (3.11<z<3.53). We estimated the ratio between ionizing (LyC flux) and 1400 A non-ionizing emissions for AGN and galaxies. Results: By running population synthesis models, we assume an average intrinsic L(1400 A)/L(900 A) ratio of 5 as the representative value for our sample. With this value and an average treatment of the lines of sight of the inter-galactic medium, we estimate the LyC escape fraction relative to the intrinsic value (fesc_rel(LyC)). We do not directly detect ionizing radiation from any individual SFG, but we are able to set a 1(2)sigma upper limit of fesc_rel(LyC)<12(24)%. This result is consistent with other non-detections published in the literature. No meaningful limits can be calculated for the sub-sample of Lya emitters. We obtain one significant direct detection for an AGN at z=3.46, with fesc_rel(LyC) = (72+/-18)%. Conclusions: Our upper limit on fescrel(LyC) implies that the SFGs studied here do not present either the physical properties or the geometric conditions suitable for efficient LyC-photon escape.

We compare the physical and morphological properties of z ~ 2 Lyman-alpha emitting galaxies (LAEs) identified in the HETDEX Pilot Survey and narrow band studies with those of z ~ 2 optical emission line galaxies (oELGs) identified via HST WFC3 infrared grism spectroscopy. Both sets of galaxies extend over the same range in stellar mass (7.5 < logM < 10.5), size (0.5 < R < 3.0 kpc), and star-formation rate (~1 < SFR < 100). Remarkably, a comparison of the most commonly used physical and morphological parameters -- stellar mass, half-light radius, UV slope, star formation rate, ellipticity, nearest neighbor distance, star formation surface density, specific star formation rate, [O III] luminosity, and [O III] equivalent width -- reveals no statistically significant differences between the populations. This suggests that the processes and conditions which regulate the escape of Ly-alpha from a z ~ 2 star-forming galaxy do not depend on these quantities. In particular, the lack of dependence on the UV slope suggests that Ly-alpha emission is not being significantly modulated by diffuse dust in the interstellar medium. We develop a simple model of Ly-alpha emission that connects LAEs to all high-redshift star forming galaxies where the escape of Ly-alpha depends on the sightline through the galaxy. Using this model, we find that mean solid angle for Ly-alpha escape is 2.4+/-0.8 steradians; this value is consistent with those calculated from other studies.

We derive the mean wavelength dependence of stellar attenuation in a sample of 239 high redshift (1.90 < z < 2.35) galaxies selected via Hubble Space Telescope (HST) WFC3 IR grism observations of their rest-frame optical emission lines. Our analysis indicates that the average reddening law follows a form similar to that derived by Calzetti et al. for local starburst galaxies. However, over the mass range 7.2 < log M/Msolar < 10.2, the slope of the attenuation law in the UV is shallower than that seen locally, and the UV slope steepens as the mass increases. These trends are in qualitative agreement with Kriek & Conroy, who found that the wavelength dependence of attenuation varies with galaxy spectral type. However, we find no evidence of an extinction "bump" at 2175 A in any of the three stellar mass bins, or in the sample as a whole. We quantify the relation between the attenuation curve and stellar mass and discuss its implications.

In this paper we perform a comprehensive study of the main sources of random and systematic errors in stellar mass measurement for galaxies using their Spectral Energy Distributions (SEDs). We use mock galaxy catalogs with simulated multi-waveband photometry (from U-band to mid-infrared) and known redshift, stellar mass, age and extinction for individual galaxies. Given different parameters affecting stellar mass measurement (photometric S/N ratios, SED fitting errors, systematic effects, the inherent degeneracies and correlated errors), we formulated different simulated galaxy catalogs to quantify these effects individually. We studied the sensitivity of stellar mass estimates to the codes/methods used, population synthesis models, star formation histories, nebular emission line contributions, photometric uncertainties, extinction and age. For each simulated galaxy, the difference between the input stellar masses and those estimated using different simulation catalogs, $\Delta\log(M)$, was calculated and used to identify the most fundamental parameters affecting stellar masses. We measured different components of the error budget, with the results listed as follows: (1). no significant bias was found among different codes/methods, with all having comparable scatter; (2). A source of error is found to be due to photometric uncertainties and low resolution in age and extinction grids; (3). The median of stellar masses among different methods provides a stable measure of the mass associated with any given galaxy; (4). The deviations in stellar mass strongly correlate with those in age, with a weaker correlation with extinction; (5). the scatter in the stellar masses due to free parameters are quantified, with the sensitivity of the stellar mass to both the population synthesis codes and inclusion of nebular emission lines studied.

We present photometry and derived redshifts from up to eleven bandpasses for 9927 galaxies in the Hubble Ultra Deep field (UDF), covering an observed wavelength range from the near-ultraviolet (NUV) to the near-infrared (NIR) with Hubble Space Telescope observations. Our Wide Field Camera 3 (WFC3)/UV F225W, F275W, and F336W image mosaics from the ultra-violet UDF (UVUDF) imaging campaign are newly calibrated to correct for charge transfer inefficiency, and use new dark calibrations to minimize background gradients and pattern noise. Our NIR WFC3/IR image mosaics combine the imaging from the UDF09 and UDF12 campaigns with CANDELS data to provide NIR coverage for the entire UDF field of view. We use aperture-matched point-spread function corrected photometry to measure photometric redshifts in the UDF, sampling both the Lyman break and Balmer break of galaxies at z~0.8-3.4, and one of the breaks over the rest of the redshift range. Our comparison of these results with a compilation of robust spectroscopic redshifts shows an improvement in the galaxy photometric redshifts by a factor of two in scatter and a factor three in outlier fraction over previous UDF catalogs. The inclusion of the new NUV data is responsible for a factor of two decrease in the outlier fraction compared to redshifts determined from only the optical and NIR data, and improves the scatter at z<0.5 and at z>2. The panchromatic coverage of the UDF from the NUV through the NIR yields robust photometric redshifts of the UDF, with the lowest outlier fraction available.

The Large Synoptic Survey Telescope (LSST) will explore the entire southern sky over 10 years starting in 2022 with unprecedented depth and time sampling in six filters, $ugrizy$. Artificial power on the scale of the 3.5 deg LSST field-of-view will contaminate measurements of baryonic acoustic oscillations (BAO), which fall at the same angular scale at redshift $z \sim 1$. Using the HEALPix framework, we demonstrate the impact of an "un-dithered" survey, in which $17\%$ of each LSST field-of-view is overlapped by neighboring observations, generating a honeycomb pattern of strongly varying survey depth and significant artificial power on BAO angular scales. We find that adopting large dithers (i.e., telescope pointing offsets) of amplitude close to the LSST field-of-view radius reduces artificial structure in the galaxy distribution by a factor of $\sim$10. We propose an observing strategy utilizing large dithers within the main survey and minimal dithers for the LSST Deep Drilling Fields. We show that applying various magnitude cutoffs can further increase survey uniformity. We find that a magnitude cut of $r < 27.3$ removes significant spurious power from the angular power spectrum with a minimal reduction in the total number of observed galaxies over the ten-year LSST run. We also determine the effectiveness of the observing strategy for Type Ia SNe and predict that the main survey will contribute $\sim$100,000 Type Ia SNe. We propose a concentrated survey where LSST observes one-third of its main survey area each year, increasing the number of main survey Type Ia SNe by a factor of $\sim$1.5, while still enabling the successful pursuit of other science drivers.

We seek to improve the accuracy of joint galaxy photometric redshift estimation and spectral energy distribution (SED) fitting. By simulating different sources of uncorrected systematic errors, we demonstrate that if the uncertainties on the photometric redshifts are estimated correctly, so are those on the other SED fitting parameters, such as stellar mass, stellar age, and dust reddening. Furthermore, we find that if the redshift uncertainties are over(under)-estimated, the uncertainties in SED parameters tend to be over(under)-estimated by similar amounts. These results hold even in the presence of severe systematics and provide, for the first time, a mechanism to validate the uncertainties on these parameters via comparison with spectroscopic redshifts. We propose a new technique (annealing) to re-calibrate the joint uncertainties in the photo-z and SED fitting parameters without compromising the performance of the SED fitting + photo-z estimation. This procedure provides a consistent estimation of the multidimensional probability distribution function in SED fitting + z parameter space, including all correlations. While the performance of joint SED fitting and photo-z estimation might be hindered by template incompleteness, we demonstrate that the latter is "flagged" by a large fraction of outliers in redshift, and that significant improvements can be achieved by using flexible stellar populations synthesis models and more realistic star formation histories. In all cases, we find that the median stellar age is better recovered than the time elapsed from the onset of star formation [abridged].

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.

Galaxies with stellar masses near M* contain the majority of stellar mass in the universe, and are therefore of special interest in the study of galaxy evolution. The Milky Way (MW) and Andromeda (M31) have present day stellar masses near M*, at 5x10^10 Msol (MW-mass) and 10^11 Msol (M31-mass). We study the typical progenitors of these galaxies using ZFOURGE, a deep medium-band near-IR imaging survey, which is sensitive to the progenitors of these galaxies out to z~3. We use abundance-matching techniques to identify the main progenitors of these galaxies at higher redshifts. We measure the evolution in the stellar mass, rest-frame colors, morphologies, far-IR luminosities, and star-formation rates combining our deep multiwavelength imaging with near-IR HST imaging from CANDELS, and far-IR imaging from GOODS-H and CANDELS-H. The typical MW-mass and M31-mass progenitors passed through the same evolution stages, evolving from blue, star-forming disk galaxies at the earliest stages, to redder dust-obscured IR-luminous galaxies in intermediate stages, and to red, more quiescent galaxies at their latest stages. The progenitors of the MW-mass galaxies reached each evolutionary stage at later times (lower redshifts) and with stellar masses that are a factor of 2-3 lower than the progenitors of the M31-mass galaxies. The process driving this evolution, including the suppression of star-formation in present-day M* galaxies requires an evolving stellar-mass/halo-mass ratio and/or evolving halo-mass threshold for quiescent galaxies. The effective size and star-formation rates imply that the baryonic cold-gas fractions drop as galaxies evolve from high redshift to z~0 and are strongly anticorrelated with an increase in the S\'ersic index. Therefore, the growth of galaxy bulges in M* galaxies corresponds to a rapid decline in the galaxy gas fractions and/or a decrease in the star-formation efficiency.

We discuss different methods to separate high- from low-redshift galaxies based on a combination of spectroscopic and photometric observations. Our baseline scenario is the Hobby-Eberly Telescope Dark Energy eXperiment (HETDEX) survey, which will observe several hundred thousand Lyman Alpha Emitting (LAE) galaxies at 1.9 < z < 3.5, and for which the main source of contamination is [OII]-emitting galaxies at z < 0.5. Additional information useful for the separation comes from empirical knowledge of LAE and [OII] luminosity functions and equivalent width distributions as a function of redshift. We consider three separation techniques: a simple cut in equivalent width, a Bayesian separation method, and machine learning algorithms, including support vector machines. These methods can be easily applied to other surveys and used on simulated data in the framework of survey planning.

We present a new technique for modeling submillimeter galaxies (SMGs): the "Count Matching" approach. Using lightcones drawn from a semi-analytic model of galaxy formation, we choose physical galaxy properties given by the model as proxies for their submillimeter luminosities, assuming a monotonic relationship. As recent interferometric observations of the Extended Chandra Deep Field South show that the brightest sources detected by single-dish telescopes are comprised by emission from multiple fainter sources, we assign the submillimeter fluxes so that the combined LABOCA plus bright-end ALMA observed number counts for this field are reproduced. After turning the model catalogs given by the proxies into submillimeter maps, we perform a source extraction to include the effects of the observational process on the recovered counts and galaxy properties. We find that for all proxies, there are lines of sight giving counts consistent with those derived from LABOCA observations, even for input sources with randomized positions in the simulated map. Comparing the recovered redshift, stellar mass and host halo mass distributions for model SMGs with observational data, we find that the best among the proposed proxies is that in which the submillimeter luminosity increases monotonically with the product between dust mass and SFR. This proxy naturally reproduces a positive trend between SFR and bolometric IR luminosity. The majority of components of blended sources are spatially unassociated.

We compare the H-beta line strengths of 1.90 < z < 2.35 star-forming galaxies observed with the near-IR grism of the Hubble Space Telescope with ground-based measurements of Ly-alpha from the HETDEX Pilot Survey and narrow-band imaging. By examining the line ratios of 73 galaxies, we show that most star-forming systems at this epoch have a Ly-alpha escape fraction below ~6%. We confirm this result by using stellar reddening to estimate the effective logarithmic extinction of the H-beta emission line (c_Hbeta = 0.5) and measuring both the H-beta and Ly-alpha luminosity functions in a ~ 100,000 cubic Mpc volume of space. We show that in our redshift window, the volumetric Ly-alpha escape fraction is at most 4.4+/-2.1(1.2)%, with an additional systematic ~25% uncertainty associated with our estimate of extinction. Finally, we demonstrate that the bulk of the epoch's star-forming galaxies have Ly-alpha emission line optical depths that are significantly greater than that for the underlying UV continuum. In our predominantly [O~III] 5007-selected sample of galaxies, resonant scattering must be important for the escape of Ly-alpha photons.

We estimate the UV continuum slope, beta, for 923 galaxies in the range 1 < z < 8 in the Hubble Ultradeep Field (HUDF). These data include 460 galaxies at 1 < z < 2 down to an absolute magnitude M_{UV} = -14 (~0.006 L*_{z=1}; 0.02 L*_{z=0}), comparable to dwarf galaxies in the local universe. We combine deep HST/UVIS photometry in F225W, F275W, F336W wavebands (UVUDF) with recent data from HST/WFC3/IR (HUDF12). Galaxies in the range 1 < z < 2 are significantly bluer than local dwarf galaxies. We find their mean (median) values <beta> = -1.382 (-1.830) +/- 0.002 (random) +/- 0.1 (systematic). We find comparable scatter in beta (standard deviation = 0.43) to local dwarf galaxies and 30% larger scatter than z > 2 galaxies. We study the trends of beta with redshift and absolute magnitude for binned sub-samples and find a modest color-magnitude relation, dbeta/dM = -0.11 +/- 0.01 and no evolution in dbeta/dM with redshift. A modest increase in dust reddening with redshift and luminosity, Delta E(B-V) ~ 0.1, and a comparable increase in the dispersion of dust reddening at z < 2, appears likely to explain the observed trends. At z > 2, we find trends that are consistent with previous works; combining our data with the literature in the range 1 < z < 8, we find a color evolution with redshift, dbeta/dz = -0.09 +/-0.01 for low luminosity (0.05 L*_{z=3}), and dbeta/dz = -0.06 +/-0.01 for medium luminosity (0.25 L*_{z=3}) galaxies.

We present the results from a VLT/SINFONI and Keck/NIRSPEC near-infrared spectroscopic survey of 16 Lyman-alpha emitters (LAEs) at $z$ = 2.1 - 2.5 in the COSMOS and GOODS-N fields discovered from the HETDEX Pilot Survey. We detect rest-frame optical nebular lines (H$\alpha$ and/or [OIII]$\lambda$5007) for 10 of the LAEs and measure physical properties, including the star formation rate (SFR), gas-phase metallicity, gas-mass fraction, and Ly$\alpha$ velocity offset. We find that LAEs may lie below the mass-metallicity relation for continuum-selected star-forming galaxies at the same redshift. The LAEs all show velocity shifts of Ly$\alpha$ relative to the systemic redshift ranging between +85 and +296 km s$^{-1}$ with a mean of +180 km s$^{-1}$. This value is smaller than measured for continuum-selected star-forming galaxies at similar redshifts. The Ly$\alpha$ velocity offsets show a moderate correlation with the measured star formation rate (2.5$\sigma$), but no significant correlations are seen with the SFR surface density, specific SFR, stellar mass, or dynamical mass ($\lesssim$ 1.5$\sigma$). Exploring the role of dust, kinematics of the interstellar medium (ISM), and geometry on the escape of Ly$\alpha$ photons, we find no signature of selective quenching of resonantly scattered Ly$\alpha$ photons. However, we also find no evidence that a clumpy ISM is enhancing the Ly$\alpha$ equivalent width. Our results suggest that the low metallicity in LAEs may be responsible for yielding an environment with a low neutral hydrogen column density as well as less dust, easing the escape of Ly$\alpha$ photons over that in continuum-selected star-forming galaxies.

We investigate the properties of damped Ly{\alpha} absorption systems (DLAs) in semi-analytic models of galaxy formation, including partitioning of cold gas in galactic discs into atomic, molecular, and ionized phases with a molecular gas-based star formation recipe. We investigate two approaches for partitioning gas into these constituents: a pressure-based and a metallicity-based recipe. We identify DLAs by passing lines of sight through our simulations to compute HI column densities. We find that models with "standard" gas radial profiles - where the average specific angular momentum of the gas disc is equal to that of the host dark matter halo - fail to reproduce the observed column density distribution of DLAs. These models also fail to reproduce the distribution of velocity widths {\Delta}v, overproducing low {\Delta}v relative to high {\Delta}v systems. Models with "extended" radial gas profiles - corresponding to gas discs with higher specific angular momentum - are able to reproduce quite well the column density distribution of absorbers over the column density range 19 < log NHI < 22.5 in the redshift range 2 < z < 3.5. The model with pressure-based gas partitioning also reproduces the observed line density of DLAs, HI gas density, and {\Delta}v distribution at z < 3 remarkably well. However all of the models investigated here underproduce DLAs and the HI gas density at z > 3. If this is the case, the flatness in the number of DLAs and HI gas density over the redshift interval 0 < z < 5 may be due to a cosmic coincidence where the majority of DLAs at z > 3 arise from intergalactic gas in filaments while those at z < 3 arise predominantly in galactic discs. We further investigate the dependence of DLA metallicity on redshift and {\Delta}v, and find reasonably good agreement with the observations, particularly when including the effects of metallicity gradients (abbrv.).

We present the rest-frame UV wavelength dependence of the Petrosian-like half-light radius ($r_{50}$), and the concentration parameter for a sample of 198 star-forming galaxies at 0.5 < z < 1.5. We find a ~5% decrease in $r_{50}$ from 1500 \AA\ to 3000 \AA, with half-light radii at 3000 \AA\ ranging from 0.6 kpc to 6 kpc. We also find a decrease in concentration of ~0.07 (1.9 < $C_{3000}$ < 3.9). The lack of a strong relationship between $r_{50}$ and wavelength is consistent with a model in which clumpy star formation is distributed over length scales comparable to the galaxy's rest-frame optical light. While the wavelength dependence of $r_{50}$ is independent of size at all redshifts, concentration decreases more sharply in the far-UV (~1500 \AA) for large galaxies at z ~ 1. This decrease in concentration is caused by a flattening of the inner ~20% of the light profile in disk-like galaxies, indicating that the central regions have different UV colors than the rest of the galaxy. We interpret this as a bulge component with older stellar populations and/or more dust. The size-dependent decrease in concentration is less dramatic at z ~ 2, suggesting that bulges are less dusty, younger, and/or less massive than the rest of the galaxy at higher redshifts.

We study the far-infrared properties of 498 Lyman Alpha Emitters (LAEs) at z=2.8, 3.1 and 4.5 in the Extended Chandra Deep Field-South, using 250, 350 and 500 micron data from the Herschel Multi-tiered Extragalactic Survey (HerMES) and 870 micron data from the LABOCA ECDFS Submillimeter Survey (LESS). None of the 126, 280 or 92 LAEs at z=2.8, 3.1 and 4.5, respectively, are individually detected in the far-infrared data. We use stacking to probe the average emission to deeper flux limits, reaching $1\sigma$ depths of ~0.1 to 0.4 mJy. The LAEs are also undetected at $\ge3\sigma$ in the stacks, although a $2.5\sigma$ signal is observed at 870 micron for the z=2.8 sources. We consider a wide range of far-infrared spectral energy distributions (SEDs), including a M82 and an Sd galaxy template, to determine upper limits on the far-infrared luminosities and far-infrared-derived star-formation rates of the LAEs. These star-formation rates are then combined with those inferred from the Ly$\alpha$ and UV emission to determine lower limits on the LAEs Ly$\alpha$ escape fraction ($f_{\rm esc}($Ly$\alpha$)). For the Sd SED template, the inferred LAEs $f_{\rm esc}($Ly$\alpha$) are $\gtrsim30%$ ($1\sigma$) at z=2.8, 3.1 and 4.5, which are all significantly higher than the global $f_{\rm esc}($Ly$\alpha$) at these redshifts. Thus, if the LAEs $f_{\rm esc}($Ly$\alpha$) follows the global evolution then they have warmer far-infrared SEDs than the Sd galaxy template. The average and M82 SEDs produce lower limits on the LAE $f_{\rm esc}($Ly$\alpha$) of ~10 to 20% ($1\sigma$), all of which are slightly higher than the global evolution of $f_{\rm esc}($Ly$\alpha$) but consistent with it at the 2 to 3$\sigma$ level.

We use broadband photometry extending from the rest-frame UV to the near-IR to fit the individual spectral energy distributions (SEDs) of 63 bright (L(Ly-alpha) > 10^43 ergs/s) Ly-alpha emitting galaxies (LAEs) in the redshift range 1.9 < z < 3.6. We find that these LAEs are quite heterogeneous, with stellar masses that span over three orders of magnitude, from 7.5 < log M < 10.5. Moreover, although most LAEs have small amounts of extinction, some high-mass objects have stellar reddenings as large as E(B-V) ~0.4. Interestingly, in dusty objects the optical depths for Ly-alpha and the UV continuum are always similar, indicating that Ly-alpha photons are not undergoing many scatters before escaping their galaxy. In contrast, the ratio of optical depths in low-reddening systems can vary widely, illustrating the diverse nature of the systems. Finally, we show that in the star formation rate (SFR)-log mass diagram, our LAEs fall above the "main-sequence" defined by z ~ 3 continuum selected star-forming galaxies. In this respect, they are similar to sub-mm-selected galaxies, although most LAEs have much lower mass.

We present an overview of a 90-orbit Hubble Space Telescope treasury program to obtain near ultraviolet imaging of the Hubble Ultra Deep Field using the Wide Field Camera 3 UVIS detector with the F225W, F275W, and F336W filters. This survey is designed to: (i) Investigate the episode of peak star formation activity in galaxies at 1<z<2.5; (ii) Probe the evolution of massive galaxies by resolving sub-galactic units (clumps); (iii) Examine the escape fraction of ionizing radiation from galaxies at z~2-3; (iv) Greatly improve the reliability of photometric redshift estimates; and (v) Measure the star formation rate efficiency of neutral atomic-dominated hydrogen gas at z~1-3. In this overview paper, we describe the survey details and data reduction challenges, including both the necessity of specialized calibrations and the effects of charge transfer inefficiency. We provide a stark demonstration of the effects of charge transfer inefficiency on resultant data products, which when uncorrected, result in uncertain photometry, elongation of morphology in the readout direction, and loss of faint sources far from the readout. We agree with the STScI recommendation that future UVIS observations that require very sensitive measurements use the instrument's capability to add background light through a "post-flash". Preliminary results on number counts of UV-selected galaxies and morphology of galaxies at z~1 are presented. We find that the number density of UV dropouts at redshifts 1.7, 2.1, and 2.7 is largely consistent with the number predicted by published luminosity functions. We also confirm that the image mosaics have sufficient sensitivity and resolution to support the analysis of the evolution of star-forming clumps, reaching 28-29th magnitude depth at 5 sigma in a 0.2 arcsecond radius aperture depending on filter and observing epoch.

We use the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) GOODS-S multi-wavelength catalog to identify counterparts for 20 Lya Emitting (LAE) galaxies at z=2.1. We build several types of stacked Spectral Energy Distributions (SEDs) of these objects. We combine photometry to form average and median flux-stacked SEDs, and postage stamp images to form average and median image-stacked SEDs. We also introduce scaled flux stacks that eliminate the influence of variation in overall brightness. We use the SED fitting code SpeedyMC to constrain the physical properties of individual objects and stacks. Our LAEs at z = 2.1 have stellar masses ranging from 2x10^7 Msun - 8x10^9 Msun (median = 3x10^8 Msun), ages ranging from 4 Myr to 500 Myr (median =100 Myr), and E(B-V) between 0.02 and 0.24 (median = 0.12). We do not observe strong correlations between Lya equivalent width (EW) and stellar mass, age, or E(B-V). The Lya radiative transfer (q) factors of our sample are predominantly close to one and do not correlate strongly with EW or E(B-V), implying that Lya radiative transfer prevents Lya photons from resonantly scattering in dusty regions. The SED parameters of the flux stacks match the average and median values of the individual objects, with the flux-scaled median SED performing best with reduced uncertainties. Median image-stacked SEDs provide a poor representation of the median individual object, and none of the stacking methods captures the large dispersion of LAE properties.

Large sets of objects with spectroscopic redshift measurements will be needed for imaging dark energy experiments to achieve their full potential, serving two goals:_training_, i.e., the use of objects with known redshift to develop and optimize photometric redshift algorithms; and_calibration_, i.e., the characterization of moments of redshift (or photo-z error) distributions. Better training makes cosmological constraints from a given experiment stronger, while highly-accurate calibration is needed for photo-z systematics not to dominate errors. In this white paper, we investigate the required scope of spectroscopic datasets which can serve both these purposes for ongoing and next-generation dark energy experiments, as well as the time required to obtain such data with instruments available in the next decade. Large time allocations on kilo-object spectrographs will be necessary, ideally augmented by infrared spectroscopy from space. Alternatively, precision calibrations could be obtained by measuring cross-correlation statistics using samples of bright objects from a large baryon acoustic oscillation experiment such as DESI. We also summarize the additional work on photometric redshift methods needed to prepare for ongoing and future dark energy experiments.

We search for evidence of diffuse Ly-alpha emission from extended neutral hydrogen surrounding Ly-alpha emitting galaxies (LAEs) using deep narrow-band images of the Extended Chandra Deep Field South. By stacking the profiles of 187 LAEs at z = 2.06, 241 LAEs at z = 3.10, and 179 LAEs at z = 3.12, and carefully performing low-surface brightness photometry, we obtain mean surface brightness maps that reach 9.9, 8.7, and 6.2 * 10^{-19} ergs cm^{-2} s^{-1} arcsec^{-2} in the emission line. We undertake a thorough investigation of systematic uncertainties in our surface brightness measurements, and find that our limits are 5--10 times larger than would be expected from Poisson background fluctuations; these uncertainties are often underestimated in the literature. At z ~ 3.1, we find evidence for extended halos with small scale lengths of 5--8 kpc in some, but not all of our sub-samples. We demonstrate that sub-samples of LAEs with low equivalent widths and brighter continuum magnitudes are more likely to possess such halos. At z ~ 2.1, we find no evidence of extended Ly-alpha emission down to our detection limits. Through Monte-Carlo simulations, we also show that we would have detected large diffuse LAE halos if they were present in our data sets. We compare these findings to other measurements in the literature, and discuss possible instrumental and astrophysical reasons for the discrepancies.

We present results from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) photometric redshift methods investigation. In this investigation, the results from eleven participants, each using a different combination of photometric redshift code, template spectral energy distributions (SEDs) and priors, are used to examine the properties of photometric redshifts applied to deep fields with broad-band multi-wavelength coverage. The photometry used includes U-band through mid-infrared filters and was derived using the TFIT method. Comparing the results, we find that there is no particular code or set of template SEDs that results in significantly better photometric redshifts compared to others. However, we find codes producing the lowest scatter and outlier fraction utilize a training sample to optimize photometric redshifts by adding zero-point offsets, template adjusting or adding extra smoothing errors. These results therefore stress the importance of the training procedure. We find a strong dependence of the photometric redshift accuracy on the signal-to-noise ratio of the photometry. On the other hand, we find a weak dependence of the photometric redshift scatter with redshift and galaxy color. We find that most photometric redshift codes quote redshift errors (e.g., 68% confidence intervals) that are too small compared to that expected from the spectroscopic control sample. We find that all codes show a statistically significant bias in the photometric redshifts. However, the bias is in all cases smaller than the scatter, the latter therefore dominates the errors. Finally, we find that combining results from multiple codes significantly decreases the photometric redshift scatter and outlier fraction. We discuss different ways of combining data to produce accurate photometric redshifts and error estimates.

We present a UV-to-mid infrared multi-wavelength catalog in the CANDELS/GOODS-S field, combining the newly obtained CANDELS HST/WFC3 F105W, F125W, and F160W data with existing public data. The catalog is based on source detection in the WFC3 F160W band. The F160W mosaic includes the data from CANDELS deep and wide observations as well as previous ERS and HUDF09 programs. The mosaic reaches a 5$\sigma$ limiting depth (within an aperture of radius 0.17 arcsec) of 27.4, 28.2, and 29.7 AB for CANDELS wide, deep, and HUDF regions, respectively. The catalog contains 34930 sources with the representative 50% completeness reaching 25.9, 26.6, and 28.1 AB in the F160W band for the three regions. In addition to WFC3 bands, the catalog also includes data from UV (U-band from both CTIO/MOSAIC and VLT/VIMOS), optical (HST/ACS F435W, F606W, F775W, F814W, and F850LP), and infrared (HST/WFC3 F098M, VLT/ISAAC Ks, VLT/HAWK-I Ks, and Spitzer/IRAC 3.6, 4.5, 5.8, 8.0 $\mu$m) observations. The catalog is validated via stellar colors, comparison with other published catalogs, zeropoint offsets determined from the best-fit templates of the spectral energy distribution of spectroscopically observed objects, and the accuracy of photometric redshifts. The catalog is able to detect unreddened star-forming (passive) galaxies with stellar mass of 10^{10}M_\odot at a 50% completeness level to z$\sim$3.4 (2.8), 4.6 (3.2), and 7.0 (4.2) in the three regions. As an example of application, the catalog is used to select both star-forming and passive galaxies at z$\sim$2--4 via the Balmer break. It is also used to study the color--magnitude diagram of galaxies at 0<z<4.

To investigate the ingredients, which allow star-forming galaxies to present Lyalpha line in emission, we studied the kinematics and gas phase metallicity (Z) of the interstellar medium. We used multi-object NIR spectroscopy with Magellan/MMIRS to study nebular emission from z=2-3 star-forming galaxies discovered in 3 MUSYC fields. We detected emission lines from four active galactic nuclei and 13 high-z star-forming galaxies, including Halpha lines down to a flux of 4.E-17 erg/sec/cm^2. This yielded 7 new redshifts. The most common emission line detected is [OIII]5007, which is sensitive to Z. We were able to measure Z for 2 galaxies and to set upper(lower) limits for another 2(2). The Z values are consistent with 0.3<Z/Zsun<1.2. Comparing the Lyalpha central wavelength with the systemic redshift, we find Delta_v(Lyalpha-[OIII])=70-270 km/sec. High-redshift star-forming galaxies, Lyalpha emitting (LAE) galaxies, and Halpha emitters appear to be located in the low mass, high star-formation rate (SFR) region of the SFR versus stellar mass diagram, confirming that they are experiencing burst episodes of star formation, which are building up their stellar mass. Their Zs are consistent with the relation found for z<2.2 galaxies in the Z versus stellar mass plane. The measured Delta_v(Lyalpha-[OIII]) values imply that outflows of material, driven by star formation, could be present in the z=2-3 LAEs of our sample. Comparing with the literature, we note that galaxies with lower Z than ours are also characterized by similar Delta_v(Lyalpha-[OIII]) velocity offsets. Strong [OIII] is detected in many Lyalpha emitters. Therefore, we propose the Lyalpha/[OIII] flux ratio as a tool for the study of high-z galaxies; while influenced by Z, ionization, and Lyalpha radiative transfer in the ISM, it may be possible to calibrate this ratio to primarily trace one of these effects.

This document summarizes the results of a community-based discussion of the potential science impact of the Mayall+BigBOSS highly multiplexed multi-object spectroscopic capability. The KPNO Mayall 4m telescope equipped with the DOE- and internationally-funded BigBOSS spectrograph offers one of the most cost-efficient ways of accomplishing many of the pressing scientific goals identified for this decade by the "New Worlds, New Horizons" report. The BigBOSS Key Project will place unprecedented constraints on cosmological parameters related to the expansion history of the universe. With the addition of an open (publicly funded) community access component, the scientific impact of BigBOSS can be extended to many important astrophysical questions related to the origin and evolution of galaxies, stars, and the IGM. Massive spectroscopy is the critical missing ingredient in numerous ongoing and planned ground- and space-based surveys, and BigBOSS is unique in its ability to provide this to the US community. BigBOSS data from community-led projects will play a vital role in the education and training of students and in maintaining US leadership in these fields of astrophysics. We urge the NSF-AST division to support community science with the BigBOSS multi-object spectrograph through the period of the BigBOSS survey in order to ensure public access to the extraordinary spectroscopic capability.

Lyman Alpha Emitting (LAE) galaxies are thought to be progenitors of present-day L* galaxies. Clustering analyses have suggested that LAEs at z ~ 3 might evolve into LAEs at z ~ 2, but it is unclear whether the physical nature of these galaxies is compatible with this hypothesis. Several groups have investigated the properties of LAEs using spectral energy distribution (SED) fitting, but direct comparison of their results is complicated by inconsistencies in the treatment of the data and in the assumptions made in modeling the stellar populations, which are degenerate with the effects of galaxy evolution. By using the same data analysis pipeline and SED fitting software on two stacked samples of LAEs at z = 3.1 and z = 2.1, and by eliminating several systematic uncertainties that might cause a discrepancy, we determine that the physical properties of these two samples of galaxies are dramatically different. LAEs at z = 3.1 are found to be old (age ~ 1 Gyr) and metal-poor (Z < 0.2 Z_Sun), while LAEs at z = 2.1 appear to be young (age ~ 50 Myr) and metal-rich (Z > Z_Sun). The difference in the observed stellar ages makes it very unlikely that z = 3.1 LAEs evolve directly into z = 2.1 LAEs. Larger samples of galaxies, studies of individual objects and spectroscopic measurements of metallicity at these redshifts are needed to confirm this picture, which is difficult to reconcile with the effects of 1 Gyr of cosmological evolution.

We determine star formation rates (SFRs) in a sample of color-selected, star-forming (sBzK) galaxies (K(AB)<21.8) in the Extended Chandra Deep Field - South. To identify and avoid active galactic nuclei, we use X-ray, IRAC color, and IR/radio flux ratio selection methods. Photometric redshift-binned, average flux densities are measured with stacking analyses in Spitzer-MIPS IR, BLAST and APEX/LABOCA submillimeter, VLA and GMRT radio and Chandra X-ray data. We include averages of aperture fluxes in MUSYC UBVRIz'JHK images to determine UV-through-radio spectral energy distributions. We determine the total IR luminosities and compare SFR calibrations from FIR, 24 micron, UV, radio and X-ray wavebands. We find consistency with our best estimator, SFR(IR+UV), to within errors for the preferred radio SFR calibration. Our results imply that 24 micron only and X-ray SFR estimates should be applied to high redshift galaxies with caution. Average IR luminosities are consistent with luminous infrared galaxies. We find SFR(IR+UV) for our stacked sBzKs at median redshifts 1.4, 1.8, and 2.2 to be 55+/-6 (random error), 74+/-8 and 154+/-17 Msun yr^-1 respectively, with additional systematic uncertainty of a factor of ~2.

Using the Millennium-II Simulation dark matter sub-halo merger histories, we created mock catalogs of Lyman Alpha Emitting (LAE) galaxies at z=3.1 to study the properties of their descendants. Several models were created by selecting the sub-halos to match the number density and typical dark matter mass determined from observations of these galaxies. We used mass-based and age-based selection criteria to study their effects on descendant populations at z~2, 1 and 0. For the models that best represent LAEs at z=3.1, the z=0 descendants have a median dark matter halo mass of 10^12.7 M_Sun, with a wide scatter in masses (50% between 10^11.8 and 10^13.7 M_Sun). Our study differentiated between central and satellite sub-halos and found that ~55% of z=0 descendants are central sub-halos with M_Median~10^12 M_Sun. This confirms that central z=0 descendants of z=3.1 LAEs have halo masses typical of L* type galaxies. The satellite sub-halos reside in group/cluster environments with dark matter masses around 10^14 M_Sun. The median descendant mass is robust to various methods of age determination, but it could vary by a factor of 5 due to current observational uncertainties in the clustering of LAEs used to determine their typical z=3.1 dark matter mass.

The spectral energy distribution (SED) of a galaxy contains information on the galaxy's physical properties, and multi-wavelength observations are needed in order to measure these properties via SED fitting. In planning these surveys, optimization of the resources is essential. The Fisher Matrix formalism can be used to quickly determine the best possible experimental setup to achieve the desired constraints on the SED fitting parameters. However, because it relies on the assumption of a Gaussian likelihood function, it is in general less accurate than other slower techniques that reconstruct the probability distribution function (PDF) from the direct comparison between models and data. We compare the uncertainties on SED fitting parameters predicted by the Fisher Matrix to the ones obtained using the more thorough PDF fitting techniques. We use both simulated spectra and real data, and consider a large variety of target galaxies differing in redshift, mass, age, star formation history, dust content, and wavelength coverage. We find that the uncertainties reported by the two methods agree within a factor of two in the vast majority (~ 90%) of cases. If the age determination is uncertain, the top-hat prior in age used in PDF fitting to prevent each galaxy from being older than the Universe needs to be incorporated in the Fisher Matrix, at least approximately, before the two methods can be properly compared. We conclude that the Fisher Matrix is a useful tool for astronomical survey design.

We present properties of individual and composite rest-UV spectra of continuum- and narrowband-selected star-forming galaxies (SFGs) at a redshift of 2<z<3.5 discovered by the MUSYC collaboration in the ECDF-S. Among our sample of 81 UV-bright SFGs, 59 have R<25.5, of which 32 have rest-frame equivalent widths W_{Ly{\alpha}}>20 {\AA}, the canonical limit to be classified as a LAE. We divide our dataset into subsamples based on properties we are able to measure for each individual galaxy: Ly{\alpha} equivalent width, rest-frame UV colors, and redshift. Among our subsample of galaxies with R<25.5, those with rest-frame W_{Ly{\alpha}}>20 {\AA} have bluer UV continua, weaker low-ionization interstellar absorption lines, weaker C IV absorption, and stronger Si II* nebular emission than those with W_{Ly{\alpha}}<20 {\AA}. We measure a typical velocity offset of {\Delta}v~600 km s$^{-1}$ between Ly{\alpha} emission and low-ionization absorption among our subsamples. We find that the interstellar component, as opposed to the stellar component, dominates the high-ionization absorption line profiles. We find the low- and high-ionization Si ionization states have similar kinematic properties, yet the low-ionization absorption is correlated with Ly$\alpha$ emission and the high-ionization absorption is not. These trends are consistent with outflowing neutral gas being in the form of neutral clouds embedded in ionized gas as previously suggested by \cite{Steidel2010}. Moreover, our galaxies with bluer UV colors have stronger Ly{\alpha} emission, weaker low-ionization absorption and more prominent nebular emission line profiles. Among our dataset, UV-bright galaxies with W_{Ly{\alpha}}>20 {\AA} exhibit weaker Ly{\alpha} emission at lower redshifts, although we caution that this could be caused by spectroscopic confirmation of low Ly{\alpha} equivalent width galaxies being harder at z~3 than z~2.

We present GalMC (Acquaviva et al 2011), our publicly available Markov Chain Monte Carlo algorithm for SED fitting, show the results obtained for a stacked sample of Lyman Alpha Emitting galaxies at z ~ 3, and discuss the dependence of the inferred SED parameters on the assumptions made in modeling the stellar populations. We also introduce SpeedyMC, a version of GalMC based on interpolation of pre-computed template libraries. While the flexibility and number of SED fitting parameters is reduced with respect to GalMC, the average running time decreases by a factor of 20,000, enabling SED fitting of each galaxy in about one second on a 2.2GHz MacBook Pro laptop, and making SpeedyMC the ideal instrument to analyze data from large photometric galaxy surveys.

We study the physical properties of 216 z ~ 2.1 LAEs discovered in an ultra-deep narrow-band MUSYC image of the ECDF-S. We fit their stacked Spectral Energy Distribution (SED) using Charlot & Bruzual templates. We consider star formation histories parametrized by the e-folding time parameter tau, allowing for exponentially decreasing (tau>0), exponentially increasing (tau<0), and constant star formation rates. These LAEs are characterized by best fit parameters and 68% confidence intervals of log(M_*/M_sun)=8.6[8.4-9.1], E(B-V)=0.22[0.00-0.31], tau=-0.02[(-4)-18] Gyr, and age_ SF=0.018[0.009-3] Gyr. Thus, we obtain robust measurements of low stellar mass and dust content, but we cannot place meaningful constraints on the age or star formation history of the LAEs. We also calculate the instantaneous SFR to be 35[0.003-170] M_sun/yr, with its average over the last 100 Myr before observation giving <SFR>_100=4[2-30] M_sun/yr. When we compare the results for the same star formation history, LAEs at z~2.1 are dustier and show higher instantaneous SFRs than z~3.1 LAEs, while the observed stellar masses of the two samples seem consistent. LAEs appear to occupy the low-mass end of the distribution of star forming galaxies at z~2. We perform SED fitting on several sub-samples selected based on photometric properties and find that LAE sub-samples at z~2.1 exhibit heterogeneous properties. The IRAC-bright, UV-bright and red LAEs have the largest stellar mass and dust reddening. The UV-faint, IRAC-faint, and high equivalent width LAE sub-samples appear less massive (<10^9 M_sun) and less dusty, with E(B-V) consistent with zero.

We present a measurement of the spatial clustering of submillimetre galaxies (SMGs) at z = 1-3. Using data from the 870 micron LESS survey, we employ a novel technique to measure the cross-correlation between SMGs and galaxies, accounting for the full probability distributions for photometric redshifts of the galaxies. From the observed projected two-point cross-correlation function we derive the linear bias and characteristic dark matter (DM) halo masses for the SMGs. We detect clustering in the cross-correlation between SMGs and galaxies at the > 4 sigma level. For the SMG autocorrelation we obtain r_0 = 7.7 (+1.8,-2.3) h^-1 Mpc, and derive a corresponding DM halo mass of log(M_halo [h^-1 M_sun]) = 12.8 (+0.3,-0.5). Based on the evolution of DM haloes derived from simulations, we show that that the z = 0 descendants of SMGs are typically massive (~2-3 L*) elliptical galaxies residing in moderate- to high-mass groups (log(M_halo [h^-1 M_sun]) = 13.3 (+0.3,-0.5). From the observed clustering we estimate an SMG lifetime of ~100 Myr, consistent with lifetimes derived from gas consumption times and star-formation timescales, although with considerable uncertainties. The clustering of SMGs at z ~ 2 is consistent with measurements for optically-selected quasi-stellar objects (QSOs), supporting evolutionary scenarios linking starbursts and QSOs. Given that SMGs reside in haloes of characteristic mass ~ 6 x 10^12 h^-1 M_sun, we demonstrate that the redshift distribution of SMGs can be described remarkably well by the combination of two effects: the cosmological growth of structure and the evolution of the molecular gas fraction in galaxies. We conclude that the powerful starbursts in SMGs likely represent a short-lived but universal phase in massive galaxy evolution, associated with the transition between cold gas-rich, star-forming galaxies and passively evolving systems. [Abridged]

We have coadded X-ray flux of all known Lyman Alpha Emitters (LAEs) in the 4 Msec Chandra Deep Field South (CDF-S) region, to place sensitive upper limits on the average unobscured star-formation rate (SFR_X) in these galaxies. A very small fraction of Lyman-Alpha galaxies in the field are individually detected in the X-rays, implying a low fraction of AGN activity. After excluding the few X-ray detected LAEs, we stack the undetected LAEs located in the 4 Ms CDF-S data and 250 ks Extended CDF-S (ECDFS) data, and compute a 1-\sigma upper limit on SFR_X < 14, 28, 28, 140, 440, 880 M$_{\sun}$ yr$^{-1}$ for LAEs located at z = 2.1, 3.1, 3.2, 4.5, 5.7 and 6.5, respectively. The upper limit of SFR_X in LAEs can be then be compared to SFR$_{Ly\alpha}$ derived from Lyman-Alpha line and thus can constrain on the Lyman-Alpha escape fraction ($f^{Esc}_{Ly\alpha}$). We derive a lower limit on f(Lyman-Alpha Escape) > 14% (84 % confidence level, 1-\sigma lower limit) for LAEs at redshift z ~ 2.1 and z ~ 3.1-3.2. At z > 4, the current LAE samples are not of sufficient size to constrain SFR_X well. By averaging all the LAEs at z> 2, the X-ray non-detection constrains f(Lyman-Alpha Escape) > 17% (84 % confidence level, 1-\sigma lower limit), and rejects f(Lyman-Alpha Escape) < 5.7% at the 99.87% confidence level from 2.1 < z < 6.5.

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

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

We describe the results of a new, wide-field survey for z=3.1 Ly-alpha emission-line galaxies (LAEs) in the Extended Chandra Deep Field South (ECDF-S). By using a nearly top-hat 5010 Angstrom filter and complementary broadband photometry from the MUSYC survey, we identify a complete sample of 141 objects with monochromatic fluxes brighter than 2.4E-17 ergs/cm^2/s and observers-frame equivalent widths greater than ~ 80 Angstroms (i.e., 20 Angstroms in the rest-frame of Ly-alpha). The bright-end of this dataset is dominated by x-ray sources and foreground objects with GALEX detections, but when these interlopers are removed, we are still left with a sample of 130 LAE candidates, 39 of which have spectroscopic confirmations. This sample overlaps the set of objects found in an earlier ECDF-S survey, but due to our filter's redder bandpass, it also includes 68 previously uncataloged sources. We confirm earlier measurements of the z=3.1 LAE emission-line luminosity function, and show that an apparent anti-correlation between equivalent width and continuum brightness is likely due to the effect of correlated errors in our heteroskedastic dataset. Finally, we compare the properties of z=3.1 LAEs to LAEs found at z=2.1. We show that in the ~1 Gyr after z~3, the LAE luminosity function evolved significantly, with L* fading by ~0.4 mag, the number density of sources with L > 1.5E42 ergs/s declining by ~50%, and the equivalent width scale-length contracting from 70^{+7}_{-5} Angstroms to 50^{+9}_{-6} Angstroms. When combined with literature results, our observations demonstrate that over the redshift range z~0 to z~4, LAEs contain less than ~10% of the star-formation rate density of the universe.

We present a rest-frame ultraviolet morphological analysis of 78 resolved, high S/N z ~ 3.1 Lyman Alpha Emitters (LAEs) in the Extended Chandra Deep Field South (ECDF-S). Using HST/ACS V -band images taken as part of the GEMS, GOODS, and HUDF surveys. For each LAE system identified via our ground-based narrow-band imaging, we have identified those LAE systems with multiple components. We measure the concentration index and present the results of our GALFIT fits for ellipticity, Sersic index, and sizes for each resolved component with S/N > 30 as well as for each LAE system with S/N > 30. The LAEs show a heterogeneous distribution of morphologies while the ma jority tend to be highly concentrated and compact in size. We only measure the morphological properties of resolved LAEs. For systems showing multiple components we also measured the morphology of the individual components. The resolved LAEs are highly concentrated (2 < C < 4) and show a similar distribution to that measured for stars, suggesting that this diagnostic is a poor discriminator near the resolution limit. The measured ellipticities for components show a distribution peaked at {\epsilon} ~ 0.55 which is significantly different from the flat distribution of ellipticities observed for local spiral galaxies and is similar to the distribution found for Lyman-break galaxies at the same redshift. There is a wide range of best-fit Sersic indices (1 < n < 10) with the majority being between 0 < n < 2. The distribution is similar to the distribution of Sersic indices seen locally. A visual inspection of the images suggests a qualitative morphological transition at n ~ 2, with small-n LAEs having extended or multimodal light distributions and relatively little diffuse emission and large-n LAEs have compact central components surrounded by diffuse emission.

We present GalMC, a MCMC algorithm designed to fit the spectral energy distributions (SED) of galaxies to infer physical properties such as age, stellar mass, dust reddening, metallicity, redshift, and star formation rate. We describe the features of the code and the extensive tests conducted to ensure that our procedure leads to unbiased parameter estimation and accurate evaluation of uncertainties. We compare its performance to grid-based algorithms, showing that the efficiency in CPU time is ~ 100 times better for MCMC for a three dimensional parameter space and increasing with the number of dimensions. We use GalMC to fit the stacked SEDs of two samples of Lyman Alpha Emitters (LAEs) at redshift z=3.1. Our fit reveals that the typical LAE detected in the IRAC 3.6 micron band has age = 0.67 [0.37 - 1.81] Gyr and stellar mass = 3.2 [2.5 - 4.2] x 10^9 M_Sun, while the typical LAE not detected at 3.6 micron has age = 0.06 [0.01-0.2] Gyr and stellar mass = 2 [1.1 - 3.4] x 10^8 M_Sun. The SEDs of both stacks are consistent with the absence of dust. The data do not significantly prefer exponential with respect to constant star formation history. The stellar populations of these two samples are consistent with the previous study by Lai et al, with some differences due to the improved modeling of the stellar populations. A constraint on the metallicity of z=3.1 LAEs from broad-band photometry, requiring Z < Z_Sun at 95% confidence, is found here for the first time.

The formation of the first massive objects in the infant Universe remains impossible to observe directly and yet it sets the stage for the subsequent evolution of galaxies. While some black holes with masses > billion solar masses? have been detected in luminous quasars less than one billion years after the Big Bang, these individual extreme objects have limited utility in constraining the channels of formation of the earliest black holes. The initial conditions of black hole seed properties are quickly erased during the growth process. From deep, optimally stacked, archival X-ray observations, we measure the amount of black hole growth in z=6-8 galaxies (0.7-1 billion years after the Big Bang). Our results imply that black holes grow in tandem with their hosts throughout cosmic history, starting from the earliest times. We find that most copiously accreting black holes at these epochs are buried in significant amounts of gas and dust that absorb most radiation except for the highest energy X-rays. This suggests that black holes grow significantly more than previously thought during these early bursts, and due to obscuration they do not contribute to the re-ionization of the Universe with their ultraviolet emission.

Aims: We present a new method that uses luminosity or stellar mass functions combined with clustering measurements to select samples of galaxies at different redshifts likely to follow a progenitor-to-descendant relationship. As the method uses clustering information, we refer to galaxy samples selected this way as clustering-selected samples. We apply this method to infer the number of mergers during the evolution of MUSYC early-type galaxies (ETGs) from z~1 to the present-day. Methods: The method consists in using clustering information to infer the typical dark-matter halo mass of the hosts of the selected progenitor galaxies. Using LambdaCDM predictions, it is then possible to follow these haloes to a later time where the sample of descendants will be that with the clustering of these descendant haloes. Results: This technique shows that ETGs at a given redshift evolve into brighter galaxies at lower redshifts (considering rest-frame, passively evolved optical luminosities). This indicates that the stellar mass of these galaxies increases with time and that, in principle, a stellar mass selection at different redshifts does not provide samples of galaxies in a progenitor-descendant relationship. Conclusions: The comparison between high redshift ETGs and their likely descendants at z=0 points to a higher number density for the progenitors by a factor 5.5+-4.0, implying the need for mergers to decrease their number density by today. Because the luminosity densities of progenitors and descendants are consistent, our results show no need for significant star-formation in ETGs since z=1, which indicates that the needed mergers are dry, i.e. gas free.

We present a rest-frame ultraviolet morphological analysis of 108 z=2.1 Lyman Alpha Emitters (LAEs) in the Extended Chandra Deep Field South (ECDF-S) and compare it to a similar sample of 171 LAEs at z=3.1. Using Hubble Space Telescope (HST) images from the Galaxy Evolution from Morphology and SEDs survey, Great Observatories Origins Deep Survey, and Hubble Ultradeep Field, we measure size and photometric component distributions, where photometric components are defined as distinct clumps of UV-continuum emission. At both redshifts, the majority of LAEs have observed half-light radii <~ 2 kpc, but the median half-light radius rises from 1.0 kpc at z=3.1 to 1.4 kpc at z=2.1. A similar evolution is seen in the sizes of individual rest-UV components, but there is no evidence for evolution in the number of multi-component systems. In the z=2.1 sample, we see clear correlations between the size of an LAE and other physical properties derived from its SED. LAEs are found to be larger for galaxies with higher stellar mass, star formation rate, and dust obscuration, but there is no evidence for a trend between equivalent width and half-light radius at either redshift. The presence of these correlations suggests that a wide range of objects are being selected by LAE surveys at z~2, including a significant fraction of objects for which a massive and moderately extended population of old stars underlies the young starburst giving rise to the Lyman alpha emission.

We study the escape of Ly-alpha photons from Ly-alpha emitting galaxies (LAEs) and the overall galaxy population using a sample of 99 LAEs at 1.9<z<3.8 detected through integral-field spectroscopy of blank fields by the HETDEX Pilot Survey. For 89 LAEs showing counterparts in deep broad-band images we measure the rest-frame UV luminosity and the UV slope, which we use to estimate E(B-V) under the assumption of a constant intrinsic UV slope for LAEs. These two quantities are used to measure the dust-corrected star formation rate (SFR). A comparison between the observed Ly-alpha luminosity and that predicted by the dust-corrected SFR yields the Ly-alpha escape fraction. We also measure the Ly-alpha luminosity function. Integration of the luminosity function provides a measurement of the Ly-alpha luminosity density across our redshift range. We combine our data with that from other surveys at 0.3<z<7.7 to trace the evolution of the Ly-alpha luminosity density. We then compare it to that expected from the star-formation history of the universe in order to characterize the evolution of the Ly-alpha escape fraction of the overall galaxy population [abriged]

We present Keck/NIRSPEC spectroscopic observations of three Lyman alpha emitting galaxies (LAEs) at z ~ 2.3 discovered with the HETDEX pilot survey. We detect Halpha, [OIII], and Hbeta emission from two galaxies at z = 2.29 and 2.49, designated HPS194 and HPS256, respectively, representing the first detection of multiple rest-frame optical emission lines in galaxies at high-redshift selected on the basis of their Lyman alpha emission. The redshifts of the Lyman alpha emission from these galaxies are offset redward of the systemic redshifts by Delta_v = 162 +/- 37 (photometric) +/- 42 (systematic) km/s for HPS194, and Delta_v = 36 +/- 35 +/- 18 km/s for HPS256. An interpretation for HPS194 is that a large-scale outflow may be occurring in its interstellar medium. The emission line ratios imply that neither LAE hosts an active galactic nucleus. Using the upper limits on the [NII] emission we place meaningful constraints on the gas-phase metallicities in these two LAEs of Z < 0.17 and < 0.28 Zsol (1 sigma). Measuring the stellar masses of these objects via spectral energy distribution (SED) fitting (~ 10^10 and 6 x 10^8 Msol, respectively), we study the nature of LAEs in a mass-metallicity plane. At least one of these two LAEs appears to be more metal poor than continuum-selected star-forming galaxies at the same redshift and stellar mass, implying that objects exhibiting Lyman alpha emission may be systematically less chemically enriched than the general galaxy population. We use the SEDs of these two galaxies to show that neglecting the emission lines when fitting stellar population models to the observed photometry can result in overestimates of the population age by orders of magnitude, and the stellar mass by a factor of ~ 2. This effect is particularly important at z > 7, where similarly strong emission lines may masquerade in the photometry as a 4000 A break (abridged).

We present the results of a comparative study of the rest-frame optical and rest-frame ultraviolet morphological properties of 117 star-forming galaxies (SFGs), including BX, BzK, and Lyman break galaxies with B<24.5, and 15 passive galaxies in the region covered by the Wide Field Camera 3 Early Release Science program. Using the internal color dispersion (ICD) diagnostic, we find that the morphological differences between the rest-frame optical and rest-frame UV light distributions in 1.4<z<2.9 SFGs are typically small (ICD~0.02). However, the majority are non-zero (56% at >3 sigma) and larger than we find in passive galaxies at 1.4<z<2, for which the weighted mean ICD is 0.013. The lack of morphological variation between individual rest-frame ultraviolet bandpasses in z~3.2 galaxies argues against large ICDs being caused by non-uniform dust distributions. Furthermore, the absence of a correlation between ICD and galaxy UV-optical color suggests that the non-zero ICDs in SFGs are produced by spatially distinct stellar populations with different ages. The SFGs with the largest ICDs (>~0.05) generally have complex morphologies that are both extended and asymmetric, suggesting that they are mergers-in-progress or very large galaxies in the act of formation. We also find a correlation between half-light radius and internal color dispersion, a fact that is not reflected by the difference in half-light radii between bandpasses. In general, we find that it is better to use diagnostics like the ICD to measure the morphological properties of the difference image than it is to measure the difference in morphological properties between bandpasses.

We present the Spitzer IRAC/MUSYC Public Legacy Survey in the Extended CDF-South (SIMPLE), which consists of deep IRAC observations covering the ~1,600 arcmin^2 area surrounding GOODS-S. The limiting magnitudes of the SIMPLE IRAC mosaics typically are 23.8, 23.6, 21.9, and 21.7, at 3.6 um, 4.5 um, 5.8 um, and 8.0 um, respectively (5-sigma total point source magnitudes in AB). The SIMPLE IRAC images are combined with the 10'x15' GOODS IRAC mosaics in the center. We give detailed descriptions of the observations, data reduction, and properties of the final images, as well as the detection and photometry methods used to build a catalog. Using published optical and near-infrared data from the Multiwavelength Survey by Yale-Chile (MUSYC), we construct an IRAC-selected catalog, containing photometry in UBVRIz'JHK, [3.6 um], [4.5 um], [5.8 um], and [8.0 um]. The catalog contains 43,782 sources with S/N > 5 at 3.6 um, 19,993 of which have 13-band photometry. We compare this catalog to the publicly available MUSYC and FIREWORKS catalogs and discuss the differences. Using a high signal-to-noise sub-sample of 3,391 sources with ([3.6] + [4.5])/2 < 21.2, we investigate the star formation rate history of massive galaxies out to z ~ 1.8. We find that at z ~ 1.8 at least 30% +/-7% of the most massive galaxies (Mstar > 10^11 Msol) are passively evolving, in agreement with earlier results from surveys covering less area.

We present a catalog of emission-line galaxies selected solely by their emission-line fluxes using a wide-field integral field spectrograph. This work is partially motivated as a pilot survey for the upcoming Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We describe the observations, reductions, detections, redshift classifications, line fluxes, and counterpart information for 397 emission-line galaxies detected over 169 sq.arcmin with a 3500-5800 Ang. bandpass under 5 Ang. full-width-half-maximum (FWHM) spectral resolution. The survey's best sensitivity for unresolved objects under photometric conditions is between 4-20 E-17 erg/s/sq.cm depending on the wavelength, and Ly-alpha luminosities between 3-6 E42 erg/s are detectable. This survey method complements narrowband and color-selection techniques in the search for high redshift galaxies with its different selection properties and large volume probed. The four survey fields within the COSMOS, GOODS-N, MUNICS, and XMM-LSS areas are rich with existing, complementary data. We find 104 galaxies via their high redshift Ly-alpha emission at 1.9<z<3.8, and the majority of the remainder objects are low redshift [OII]3727 emitters at z<0.56. The classification between low and high redshift objects depends on rest frame equivalent width, as well as other indicators, where available. Based on matches to X-ray catalogs, the active galactic nuclei (AGN) fraction amongst the Ly-alpha emitters (LAEs) is 6%. We also analyze the survey's completeness and contamination properties through simulations. We find five high-z, highly-significant, resolved objects with full-width-half-maximum sizes >44 sq.arcsec which appear to be extended Ly-alpha nebulae. We also find three high-z objects with rest frame Ly-alpha equivalent widths above the level believed to be achievable with normal star formation, EW(rest)>240 Ang.

A wide range of models describing modifications to General Relativity have been proposed, but no fundamental parameter set exists to describe them. Similarly, no fundamental theory exists for dark energy to parameterize its potential deviation from a cosmological constant. This motivates a model-independent search for deviations from the concordance GR+LambdaCDM cosmological model in large galaxy redshift surveys. We describe two model-independent tests of the growth of cosmological structure, in the form of quantities that must equal one if GR+LambdaCDM is correct. The first, epsilon, was introduced previously as a scale-independent consistency check between the expansion history and structure growth. The second, upsilon, is introduced here as a test of scale-dependence in the linear evolution of matter density perturbations. We show that the ongoing and near-future galaxy redshift surveys WiggleZ, BOSS, and HETDEX will constrain these quantities at the 5-10% level, representing a stringent test of concordance cosmology at different redshifts. When redshift space distortions are used to probe the growth of cosmological structure, galaxies at higher redshift with lower bias are found to be most powerful in detecting deviations from the GR+LambdaCDM model.

We present MIPS 24 micron observations of the Hubble Deep Field South taken with the Spitzer Space Telescope. The resulting image is 254 arcmin^2 in size and has a sensitivity ranging between ~12 to ~30 microJy rms, with a median sensitivity of ~20 microJy rms. A total of 495 sources have been cataloged with a signal-to-noise ratio greater than 5 sigma. The source catalog is presented as well as source counts which have been corrected for completeness and flux boosting. The IR sources are then combined with MUSYC optical/NIR and ATHDFS radio observations to obtain redshifts and radio flux densities of the sample. We use the IR/radio flux density ratio (q_24) to explore the IR-radio correlation for this IR sample and find q_24 = 0.71 +- 0.31 for sources detected in both IR and radio. The results are extended by stacking IR sources not detected in the radio observations and we derive an average q_24 for redshift bins between 0 < z < 2.5. We find the high redshift (z > 1) sources have an average q_{24} ratio which is better fit by local LIRG SEDs rather than local ULIRG SEDs, indicating that high redshift ULIRGs differ in their IR/radio properties. So ULIRGs at high redshift have SEDs different from those found locally. Infrared faint radio sources are examined, and while nine radio sources do not have a MIPS detection and are therefore radio-loud AGN, only one radio source has an extreme IRAC 3.6 micron to radio flux density ratio indicating it is a radio-loud AGN at z > 1.

We present deep optical 18-medium-band photometry from the Subaru telescope over the ~30' x 30' Extended Chandra Deep Field-South (ECDF-S), as part of the Multiwavelength Survey by Yale-Chile (MUSYC). This field has a wealth of ground- and space-based ancillary data, and contains the GOODS-South field and the Hubble Ultra Deep Field. We combine the Subaru imaging with existing UBVRIzJHK and Spitzer IRAC images to create a uniform catalog. Detecting sources in the MUSYC BVR image we find ~40,000 galaxies with R_AB<25.3, the median 5 sigma limit of the 18 medium bands. Photometric redshifts are determined using the EAZY code and compared to ~2000 spectroscopic redshifts in this field. The medium band filters provide very accurate redshifts for the (bright) subset of galaxies with spectroscopic redshifts, particularly at 0.1 < z < 1.2 and at z > 3.5. For 0.1 < z < 1.2, we find a 1 sigma scatter in \Delta z/(1+z) of 0.007, similar to results obtained with a similar filter set in the COSMOS field. As a demonstration of the data quality, we show that the red sequence and blue cloud can be cleanly identified in rest-frame color-magnitude diagrams at 0.1 < z < 1.2. We find that ~20% of the red-sequence-galaxies show evidence of dust-emission at longer rest-frame wavelengths. The reduced images, photometric catalog, and photometric redshifts are provided through the public MUSYC website.

Using new, highly accurate photometric redshifts from the MUSYC medium-band survey in the Extended Chandra Deep Field South (ECDF-S), we fit synthetic stellar population models to compare AGN host galaxies to inactive galaxies at 0.8 < z < 1.2. We find that AGN host galaxies are predominantly massive galaxies on the red sequence and in the green valley of the color-mass diagram. Because both passive and dusty galaxies can appear red in optical colors, we use rest-frame near-infrared colors to separate passively evolving stellar populations from galaxies that are reddened by dust. As with the overall galaxy population, ~25% of the `red' AGN host galaxies and ~75% of the `green' AGN host galaxies have colors consistent with young stellar populations reddened by dust. The dust-corrected rest-frame optical colors are the blue colors of star-forming galaxies, which implies that these AGN hosts are not passively aging to the red sequence. At z~1, AGN activity is roughly evenly split between two modes of black hole growth: the first in passively evolving host galaxies, which may be heating up the galaxy's gas and preventing future episodes of star formation, and the second in dust-reddened young galaxies, which may be ionizing the galaxy's interstellar medium and shutting down star formation.

We measure the evolution of galaxy clustering out to a redshift of z~1.5 using data from two MUSYC fields, the Extended Hubble Deep Field South (EHDF-S) and the Extended Chandra Deep Field South (ECDF-S). We use photometric redshift information to calculate the projected-angular correlation function, omega(sigma), from which we infer the projected correlation function Xi(sigma). We demonstrate that this technique delivers accurate measurements of clustering even when large redshift measurement errors affect the data. To this aim we use two mock MUSYC fields extracted from a LambdaCDM simulation populated with GALFORM semi-analytic galaxies which allow us to assess the degree of accuracy of our estimates of Xi(sigma) and to identify and correct for systematic effects in our measurements. We study the evolution of clustering for volume limited subsamples of galaxies selected using their photometric redshifts and rest-frame r-band absolute magnitudes. We find that the real-space correlation length r_0 of bright galaxies, M_r<-21 (rest-frame) can be accurately recovered out to z~1.5, particularly for ECDF-S given its near-infrared photometric coverage. There is mild evidence for a luminosity dependent clustering in both fields at the low redshift samples (up to <z>=0.57), where the correlation length is higher for brighter galaxies by up to 1Mpc/h between median rest-frame r-band absolute magnitudes of -18 to -21.5. As a result of the photometric redshift measurement, each galaxy is assigned a best-fit template; we restrict to E and E+20%Sbc types to construct subsamples of early type galaxies (ETGs). Our ETG samples show a strong increase in r_0 as the redshift increases, making it unlikely (95% level) that ETGs at median redshift z_med=1.15 are the direct progenitors of ETGs at z_med=0.37 with equivalent passively evolved luminosities. (ABRIDGED)

[abridged] We derive photometric redshifts from 17-band optical to mid-IR photometry of 74 robust counterparts to 68 of the 126 submillimetre galaxies (SMGs) selected at 870um by LABOCA observations in the ECDFS. The median photometric redshift of identified SMGs is z=2.2\pm0.1, the interquartile range is z=1.8-2.7 and we identify 10 (~15%) high-redshift (z>3) SMGs. We derive a simple redshift estimator for SMGs based on the 3.6 and 8um fluxes, which is accurate to Delta_z~0.4 for SMGs at z<4. A statistical analysis of sources around unidentified SMGs identifies a population of likely counterparts with a redshift distribution peaking at z=2.5\pm0.3, which likely comprises ~60% of the unidentified SMGs. This confirms that the bulk of the undetected SMGs are coeval with those detected in the radio/mid-IR. We conclude that ~15% of all the SMGs are below the flux limits of our survey and lie at z>3 and hence ~30% of all SMGs have z>3. We estimate that the full S_870um>4mJy SMG population has a median redshift of 2.5\pm0.6. In contrast to previous suggestions we find no significant correlation between S_870um and redshift. The median stellar mass of the SMGs derived from SED fitting is (9.2\pm0.9)x10^10Msun and the interquartile range is (4.7-14)x10^10Msun, although we caution that uncertainty in the star-formation histories results in a factor of ~5 uncertainty in these stellar masses. The median characteristic dust temperature of SMGs is 35.9\pm1.4K and the interquartile range is 28.5-43.3K. The infrared luminosity function shows that SMGs at z=2-3 typically have higher far-IR luminosities and luminosity density than those at z=1-2. This is mirrored in the evolution of the star-formation rate density (SFRD) for SMGs which peaks at z~2. The maximum contribution of bright SMGs to the global SFRD (~5% for SMGs with S_870um>4mJy; ~50% for SMGs with S_870um>1mJy) also occurs at z~2.

We present the results of a high-spatial-resolution study of the line emission in a sample of z=3.1 Lyman-Alpha-Emitting Galaxies (LAEs) in the Extended Chandra Deep Field-South. Of the eight objects with coverage in our HST/WFPC2 narrow-band imaging, two have clear detections and an additional two are barely detected (~2-sigma). The clear detections are within ~0.5 kpc of the centroid of the corresponding rest-UV continuum source, suggesting that the line-emitting gas and young stars in LAEs are spatially coincident. The brightest object exhibits extended emission with a half-light radius of ~1.5 kpc, but a stack of the remaining LAE surface brightness profiles is consistent with the WFPC2 point spread function. This suggests that the Lyman Alpha emission in these objects originates from a compact (<~2 kpc) region and cannot be significantly more extended than the far-UV continuum emission (<~1 kpc). Comparing our WFPC2 photometry to previous ground-based measurements of their monochromatic fluxes, we find at 95% (99.7%) confidence that we cannot be missing more than 22% (32%) of the Lyman Alpha emission.

We discovered a sample of 250 Ly-Alpha emitting (LAE) galaxies at z=2.1 in an ultra-deep 3727 A narrow-band MUSYC image of the Extended Chandra Deep Field-South. LAEs were selected to have rest-frame equivalent widths (EW) > 20 A and emission line fluxes > 2.0 x 10^(-17)erg /cm^2/s, after carefully subtracting the continuum contributions from narrow band photometry. The median flux of our sample is 4.2 x 10^(-17)erg/cm^2/s, corresponding to a median Lya luminosity = 1.3 x 10^(42) erg/s at z=2.1. At this flux our sample is > 90% complete. Approximately 4% of the original NB-selected candidates were detected in X-rays by Chandra, and 7% were detected in the rest-frame far-UV by GALEX. At luminosity>1.3 x 10^42 erg/s, the equivalent width distribution is unbiased and is represented by an exponential with scale-length of 83+/-10 A. Above this same luminosity threshold, we find a number density of 1.5+/-0.5 x 10^-3 Mpc^-3. Neither the number density of LAEs nor the scale-length of their EW distribution show significant evolution from z=3 to z=2. We used the rest frame UV luminosity to estimate a median star formation rate of 4 M_(sun) /yr. The median rest frame UV slope, parametrized by B-R, is that typical of dust-free, 0.5-1 Gyr old or moderately dusty, 300-500 Myr old populations. Approximately 40% of the sample occupies the z~2 star-forming galaxy locus in the UVR two color diagram. Clustering analysis reveals that LAEs at z=2.1 have r_0=4.8+/-0.9 Mpc and a bias factor b=1.8+/-0.3. This implies that z=2.1 LAEs reside in dark matter halos with median masses Log(M/M_(sun))=11.5^(+0.4)_(-0.5), which are among of the lowest-mass halos yet probed at this redshift. We used the Sheth-Tormen conditional mass function to study the descendants of these LAEs and found that their typical present-day descendants are local galaxies with L* properties, like the Milky Way.

Stacking analysis is a means of detecting faint sources using a priori position information to estimate an aggregate signal from individually undetected objects. Confusion severely limits the effectiveness of stacking in deep surveys with limited angular resolution, particularly at far infrared to submillimeter wavelengths, and causes a bias in stacking results. Deblending corrects measured fluxes for confusion from adjacent sources; however, we find that standard deblending methods only reduce the bias by roughly a factor of two while tripling the variance. We present an improved algorithm for simultaneous stacking and deblending that greatly reduces bias in the flux estimate with nearly minimum variance. When confusion from neighboring sources is the dominant error, our method improves upon RMS error by at least a factor of three and as much as an order of magnitude compared to other algorithms. This improvement will be useful for Herschel and other telescopes working in a source confused, low signal to noise regime.

We study the co-existence of star formation and AGN activity in X-ray selected AGN by analyzing stacked 870um submm emission from a deep and wide map of the Extended Chandra Deep Field South, obtained with LABOCA at the APEX telescope. The total X-ray sample of 895 sources with median redshift z~1 is detected at a mean submm flux of 0.49+-0.04mJy, corresponding to a typical star formation rate around 30Msun/yr for a T=35K, beta=1.5 greybody far-infrared SED. The good S/N permits stacking analyses for subgroups. We observe a trend of star formation rate increasing with redshift. An increase of star formation rate with AGN luminosity is indicated at the highest L_2-10>~1E44erg/s luminosities only. Increasing trends with X-ray obscuration as expected in some AGN evolutionary scenarios are not observed for the bulk of the X-ray AGN sample but may be present for the highest intrinsic luminosity objects. This suggests a transition between two modes in the coexistence of AGN activity and star formation. For the bulk of the sample, the X-ray luminosity and obscuration of the AGN are not intimately linked to the global star formation rate of their hosts. The hosts are likely massive and forming stars secularly, at rates similar to the pervasive star formation seen in massive galaxies without an AGN at similar redshifts. The change indicated towards more intense star formation, and a more pronounced increase in star formation rates between unobscured and obscured AGN at highest luminosities suggests that luminous AGN follow an evolutionary path on which obscured AGN activity and intense star formation are linked, possibly via merging. Comparison to local hard X-ray selected AGN supports this interpretation. [Abridged]

We constrain the number density and evolution of Compton-thick Active Galactic Nuclei (AGN), and their contribution to the extragalactic X-ray background. In the local Universe we use the wide area surveys from the Swift and INTEGRAL satellites, while for high redshifts we explore candidate selections based on mid-IR parameters. We present the properties of a sample of 211 heavily-obscured AGN candidates in the Extended Chandra Deep Field-South (ECDF-S) selecting objects with f24/fR>1000 and R-K>4.5. The X-ray to mid-IR ratios for these sources are significantly larger than that of star-forming galaxies and ~2 orders of magnitude smaller than for the general AGN population, suggesting column densities of NH>5x10^24 cm^-2. The space density of CT AGN at z~2 derived from these observations is ~10^-5 Mpc^{-3}, finding a strong evolution in the number of LX>10^44 erg/s sources from z=1.5 to 2.5.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We study the properties of a sample of 211 heavily-obscured Active Galactic Nucleus (AGN) candidates in the Extended Chandra Deep Field-South selecting objects with f_24/f_R>1000 and R-K>4.5. Of these, 18 were detected in X-rays and found to be obscured AGN with neutral hydrogen column densities of ~10^23 cm^-2. In the X-ray undetected sample, the following evidence suggests a large fraction of heavily-obscured (Compton Thick) AGN: (i) The stacked X-ray signal of the sample is strong, with an observed ratio of soft to hard X-ray counts consistent with a population of ~90% heavily obscured AGN combined with 10% star-forming galaxies. (ii) The X-ray to mid-IR ratios for these sources are significantly larger than that of star-forming galaxies and ~2 orders of magnitude smaller than for the general AGN population, suggesting column densities of N_H>5x10^24 cm^-2. (iii) The Spitzer near- and mid-IR colors of these sources are consistent with those of the X-ray-detected sample if the effects of dust self-absorption are considered. Spectral fitting to the rest-frame UV/optical light (dominated by the host galaxy) returns stellar masses of ~10^11 M_sun and <E(B-V)> =0.5, and reveals evidence for a significant young stellar population, indicating that these sources are experiencing considerable star-formation. This sample of heavily-obscured AGN candidates implies a space density at z~2 of ~10^-5 Mpc^-3, finding a strong evolution in the number of L_X>10^44 erg/s sources from z=1.5 to 2.5, possibly consistent with a short-lived heavily-obscured phase before an unobscured quasar is visible.

We present a sensitive 870 micron survey of the Extended Chandra Deep Field South (ECDFS) using LABOCA on the APEX telescope. The LABOCA ECDFS Submillimetre Survey (LESS) covers the full 30' x 30' field size of the ECDFS and has a uniform noise level of 1.2 mJy/beam. LESS is thus the largest contiguous deep submillimetre survey undertaken to date. The noise properties of our map show clear evidence that we are beginning to be affected by confusion noise. We present a catalog of 126 SMGs detected with a significance level above 3.7 sigma. The ECDFS exhibits a deficit of bright SMGs relative to previously studied blank fields but not of normal star-forming galaxies that dominate the extragalactic background light (EBL). This is in line with the underdensities observed for optically defined high redshift source populations in the ECDFS (BzKs, DRGs,optically bright AGN and massive K-band selected galaxies). The differential source counts in the full field are well described by a power law with a slope of alpha=-3.2, comparable to the results from other fields. We show that the shape of the source counts is not uniform across the field. The integrated 870 micron flux densities of our source-count models account for >65% of the estimated EBL from COBE measurements. We have investigated the clustering of SMGs in the ECDFS by means of a two-point correlation function and find evidence for strong clustering on angular scales <1'. Assuming a power law dependence for the correlation function and a typical redshift distribution for the SMGs we derive a spatial correlation length of r_0=13+/-6 h^-1 Mpc.

We present a rest-frame ultraviolet analysis of ~120 z~3.1 Lyman Alpha Emitters (LAEs) in the Extended Chandra Deep Field South (ECDF-S). Using Hubble Space Telescope (HST) images taken as part of the Galaxy Evolution From Morphology and SEDS (GEMS) survey, Great Observatories Origins Deep Survey (GOODS), and Hubble Ultradeep Field (HUDF) surveys, we analyze the sizes of LAEs, as well as the spatial distribution of their components, which are defined as distinct clumps of UV-continuum emission. We set an upper limit of ~1 kpc (~0.1") on the rms offset between the centroids of the continuum and Lyman-alpha emission. The star formation rates of LAE components inferred from the rest-frame ultraviolet continuum range from ~0.1 M_{sun}/yr to ~5$ M_{sun}/yr. A subsample of LAEs with coverage in multiple surveys (at different imaging depths) suggests that one needs a signal-to-noise ratio, S/N>~30, in order to make a robust estimate of the half-light radius of an LAE system. The majority of LAEs have observed half-light radii <~2 kpc, and LAE components typically have observed half-light radii <~1.5 kpc (<~ 0.2"). Although only ~50% of the detected LAE components are resolved at GOODS depth, the brightest (V<~26.3) are all resolved in both GOODS and GEMS. Since we find little evidence for a correlation between the rest-UV sizes and magnitudes of LAEs, the majority should be resolved in a deeper survey at the ~0.05" angular resolution of HST. Most of the multi-component LAEs identified in shallow frames become connected in deeper images, suggesting that the majority of the rest-UV "clumps" are individual star-forming regions within a single system.

We present a new maximum likelihood method for the calculation of galaxy luminosity functions from multi-band photometric surveys without spectroscopic data. The method evaluates the likelihood of a trial luminosity function by directly comparing the predicted distribution of fluxes in a multi-dimensional photometric space to the observations, and thus does not require the intermediate step of calculating photometric redshifts. We apply this algorithm to ~27,000 galaxies with m_R<=25 in the MUSYC-ECDFS field, with a focus on recovering the luminosity function of field galaxies at z<1.2. Our deepest LFs reach M_r=-14 and show that the field galaxy LF deviates from a Schechter function, exhibiting a steep upturn at intermediate magnitudes that is due to galaxies of late spectral types.

Spectral Energy Distribution (SED) fitting is a well-developed astrophysical tool that has recently been applied to high-redshift Lyman Alpha-emitting galaxies. If rest-frame ultraviolet through near-infrared photometry is available, it allows the simultaneous determination of the star formation history and dust extinction of a galaxy. Lyman Alpha-emitter SED fitting results from the literature find star formation rates ~3 M_sun/yr, stellar masses ~10^9 M_sun for the general population but ~10^10 M_sun for the subset detected by IRAC, and very low dust extinction, A_V < 0.3, although a couple of outlying analyses prefer significantly more dust and higher intrinsic star formation rates. A checklist of 14 critical choices that must be made when performing SED fitting is discussed.

This publication contains the conference summary of the Understanding Lyman-alpha Emitters conference held at the Max Planck Institute for Astronomy in Heidelberg October 6 - 10, 2008. The scope of the conference was to bring together most of the scientists working in the field of Lyman-alpha emitters, whether at low or high redshift, or on observational or theoretical aspects, and to summarise how far the field of study of galaxies with Lyman-alpha emission has come. An outlook towards the future of the field was also desired. As part of the conference, two days were dedicated to in total six discussion sessions. The topics were i) new methods and selection methods, ii) morphology, iii) what can the local Universe observations tell us about the high redshift Universe?, iv) clustering, v) SED fitting and vi) Ly-alpha blobs. The chairs of those sessions were asked to summarise the discussions, as presented in these proceedings.

Using the 330hr ESO-MPG 870-micron survey of the Extended Chandra Deep Field South (ECDF-S) obtained with the Large Apex BOlometer CAmera (LABOCA) on the Atacama Pathfinder EXperiment (APEX), we have carried out a stacking analysis at submillimeter (submm) wavelengths of a sample of 8266 near-infra-red (near-IR) selected (K_vega <= 20) galaxies, including 893 BzK galaxies, 1253 extremely red objects (EROs) and 737 distant red galaxies (DRGs), selected from the Multi-wavelength Survey by Yale-Chile (MUSYC). We measure average 870-micron fluxes of 0.20+-0.01mJy (20.0sigma), 0.45+-0.04mJy (11.3sigma), 0.42+-0.03mJy (14.0sigma), and 0.41+-0.04mJy (10.3sigma) for the K_vega <= 20, BzK, ERO and DRG samples, respectively. For the BzK, ERO and DRG subsamples, which overlap to some degree and are like to be at z ~ 1-2, this implies an average far-IR luminosity of ~2-6x10^{11} Lsolar and star formation rate of ~40-100Msolar. Splitting the BzK galaxies up into star-forming (sBzK) and passive (pBzK) galaxies, the former is significantly detected (0.48+-0.04mJy, 12.0sigma) while the latter is only marginally detected (0.27+-0.10mJy, 2.7sigma), thus confirming that the sBzK/pBzK criteria do isolate obscured, star forming and truly passive galaxies. The K_vega <= 20 galaxies are found to contribute with 6.6+-0.3Jy deg^{-2} (~15%) to the 870-micron extragalactic background light (EBL). sBzK and pBzK galaxies contribute 1.7+-0.2Jy deg^{-2} (~4%) and 0.2+-0.1 Jy deg^{-2} (< 0.5%) to the EBL. [Abridged]

Promoting racial and ethnic diversity is critically important to the future success and growth of the field of astronomy. The raw ability, drive and interest required to excel in the field is distributed without regard to race, gender, or socioeconomic background. By not actively promoting diversity in our field we risk losing talented people to other professions (or losing them entirely), which means that there will be astronomical discoveries that simply won't get made. There is demonstrated evidence that STEM fields benefit from diverse perspectives on problems that require more complex thought processes. This is especially relevant to a field like astronomy where more and more work is being done collaboratively. The lack of notable growth in African American, Hispanic, and Native American representation in astronomy indicates that the 'pipeline' for these individuals is systemically leaky at critical junctures. Substantially more effort must be directed toward improving the educational and career development of minorities to insure that these potential colleagues are supported through the process. However, simply recognizing that the pipeline is faulty is woefully inadequate. There must be very specific, targeted solutions to help improve the situation. With this in mind, we offer two position papers addressing specific areas of improvement that we identify as (a) essential for any foreseeable progress in the field, and (b) attainable in the 2010-2020 decade. These position papers focus primarily on African Americans, Hispanics, and Native Americans. Although we do not directly address issues of Asian Americans, Pacific Islanders, and other groups, many of the recommendations made here can be adapted to address issues faced by these groups as well.

The NSF's Astronomy and Astrophysics Postdoctoral Fellowship (AAPF) is exceptional among the available postdoctoral awards in Astronomy and Astrophysics. The fellowship is one of the few that allows postdoctoral researchers to pursue an original research program, of their own design, at the U.S. institution of their choice. However, what makes this fellowship truly unique is the ability of Fellows to lead an equally challenging, original educational program simultaneously. The legacy of this singular fellowship has been to encourage and advance leaders in the field who are equally as passionate about their own research as they are about sharing that research and their passion for astronomy with students and the public. In this positional paper we address the importance of fellowships like the AAPF to the astronomical profession by identifying the science and educational contributions that Fellows have made to the community. Further, we recommend that fellowships that encourage leading postdoctoral researchers to also become leaders in Astronomy education be continued and expanded.

In order to attract, recruit and retain underrepresented minority students to pursue Astronomy and related fields, we must ensure that there continues to be a well qualified pool of graduate and undergraduate students from which to recruit. This required pool of people are today's elementary, middle and high school students. The Astronomy community must be proactive in demonstrating the importance of pursing scientific study and careers to these students and their parents. Only by actively engaging these communities can U.S Astronomy hope to increase the numbers of minority PhDs and continue to be a leader in Astronomical discovery and knowledge.

If the ethnic makeup of the astronomy profession is to achieve parity with the general population within one generation (~30 years), the number of underrepresented minorities earning graduate degrees in astronomy and astrophysics must increase in the coming decade by a factor of 5 to 10. To accomplish this, the profession must develop and invest in mechanisms to more effectively move individuals across critical educational junctures to the PhD and beyond. Early and continuous research engagement starting in the undergraduate years is critical to this vision, in which the federally funded research internship programs (e.g. NSF REU, NASA GSRP) and national centers/observatories play a vital role. Regionally based partnerships with minority-serving institutions (MSIs) are crucial for tapping extant pools of minority talent, as are post-baccalaurate and/or masters degree "bridging" programs that provide critical stepping stones to the PhD. Because of the strong undergraduate physics, engineering, and computer science backgrounds of many students from MSIs, we suggest that instrument development and large scale computing/data-mining are particularly promising avenues for engagement in the coming decade.

We present a new K-selected, optical-to-near-infrared photometric catalog of the Extended Chandra Deep Field South (ECDFS), making it publicly available to the astronomical community. The dataset is founded on publicly available imaging, supplemented by original zJK imaging data obtained as part of the MUltiwavelength Survey by Yale-Chile (MUSYC). The final photometric catalog consists of photometry derived from nine band U-K imaging covering the full 0.5x0.5 sq. deg. of the ECDFS, plus H band data for approximately 80% of the field. The 5sigma flux limit for point-sources is K = 22.0 (AB). This is also the nominal completeness and reliability limit of the catalog: the empirical completeness for 21.75 < K < 22.00 is 85+%. We have verified the quality of the catalog through both internal consistency checks, and comparisons to other existing and publicly available catalogs. As well as the photometric catalog, we also present catalogs of photometric redshifts and restframe photometry derived from the ten band photometry. We have collected robust spectroscopic redshift determinations from published sources for 1966 galaxies in the catalog. Based on these sources, we have achieved a (1sigma) photometric redshift accuracy of Dz/(1+z) = 0.036, with an outlier fraction of 7.8%. Most of these outliers are X-ray sources. Finally, we describe and release a utility for interpolating restframe photometry from observed SEDs, dubbed InterRest. Particularly in concert with the wealth of already publicly available data in the ECDFS, this new MUSYC catalog provides an excellent resource for studying the changing properties of the massive galaxy population at z < 2. (Abridged)

We report on the identification of the highest redshift submm-selected source currently known: LESSJ033229.4-275619. This source was detected in the Large Apex BOlometer CAmera (LABOCA) Extended Chandra Deep Field South (ECDFS) Submillimetre Survey (LESS), a sensitive 870-um survey (~1.2-mJy rms) of the full 30'x30' ECDFS with the LABOCA camera on the Atacama Pathfinder EXperiment (APEX) telescope. The submm emission is identified with a radio counterpart for which optical spectroscopy provides a redshift of z=4.76. We show that the bolometric emission is dominated by a starburst with a star formation rate of ~1000 Msun/yr, although we also identify a moderate luminosity Active Galactic Nucleus (AGN) in this galaxy. Thus it has characteristics similar to those of z~2 submm galaxies (SMGs), with a mix of starburst and obscured AGN signatures. This demonstrates that ultraluminous starburst activity is not just restricted to the hosts of the most luminous (and hence rare) QSOs at z~5, but was also occurring in less extreme galaxies at a time when the Universe was less than 10% of its current age. Assuming that we are seeing the major phase of star formation in this galaxy, then we demonstrate that it would be identified as a luminous distant red galaxy at z~3 and that the current estimate of the space density of z>4 SMGs is only sufficient to produce ~10% of the luminous red galaxy population at these early times. However, this leaves open the possibility that some of these galaxies formed through less intense, but more extended star formation events. If the progenitors of all of the luminous red galaxies at z~3 go through an ultraluminous starburst at z>4 then the required volume density of z>4 SMGs will exceed that predicted by current galaxy formation models by more than an order of magnitude.

We present the first results of our optical spectroscopy program aimed to provide redshifts and identifications for the X-ray sources in the Extended Chandra Deep Field South. A total of 339 sources were targeted using the IMACS spectrograph at the Magellan telescopes and the VIMOS spectrograph at the VLT. We measured redshifts for 186 X-ray sources, including archival data and a literature search. We find that the AGN host galaxies have on average redder rest-frame optical colors than non-active galaxies, and that they live mostly in the "green valley". The dependence of the fraction of AGN that are obscured on both luminosity and redshift is confirmed at high significance and the observed AGN space density is compared with the expectations from existing luminosity functions. These AGN show a significant difference in the mid-IR to X-ray flux ratio for obscured and unobscured AGN, which can be explained by the effects of dust self-absorption on the former. This difference is larger for lower luminosity sources, which is consistent with the dust opening angle depending on AGN luminosity.

We present the color-magnitude and color-stellar mass diagrams for galaxies with z_phot < ~2, based on a K < 22 (AB) catalog of the Extended Chandra Deep Field South (ECDFS) from the MUltiwavelength Survey by Yale-Chile (MUSYC). Our main sample of 7840 galaxies contains 1297 M_* > 10^11 M_Sol galaxies in the range 0.2 < z_phot < 1.8. We show empirically that this catalog is approximately complete for M_* > 10^11 M_Sol galaxies for z_phot < 1.8. For this mass-limited sample, we show that the locus of the red sequence color-stellar mass relation evolves as Del(u-r) ~ (-0.44+/-0.02) z_phot for z_phot < ~1.2. For z_phot > ~1.3, however, we are no longer able to reliably distinguish red and blue subpopulations from the observed color distribution; we show that this would require much deeper near infrared data. At 1.5 < z_phot <1.8, the comoving number density of M_* > 10^11 M_Sol galaxies is ~50% of the local value, with a red fraction of ~33%. Making a parametric fit to the observed evolution, we find n_tot(z) ~ (1+z_phot)^(-0.52+/-0.12(+/-0.20)). We find stronger evolution in the red fraction: f_red(z) ~ (1+z_phot)^(-1.17+/-0.18(+/-0.21)). Through a series of sensitivity analyses, we show that the most important sources of systematic error are: 1. systematic differences in the analysis of the z~0 and z>>0 samples; 2. systematic effects associated with details of the photometric redshift calculation; and 3. uncertainties in the photometric calibration. With this in mind, we show that our results based on photometric redshifts are consistent with a completely independent analysis which does not require redshift information for individual galaxies. Our results suggest that, at most, 1/5 of local red sequence galaxies with M_* >10^11 M_Sol were already in place at z ~ 2.

We study the evolution of the star formation rate (SFR) of mid-infrared (IR) selected galaxies in the extended Chandra Deep Field South (E-CDFS). We use a combination of U-K GaBoDS and MUSYC data, deep IRAC observations from SIMPLE, and deep MIPS data from FIDEL. This unique multi-wavelength data set allows us to investigate the SFR history of massive galaxies out to redshift z ~ 1.8. We determine star formation rates using both the rest-frame ultraviolet luminosity from young, hot stars and the total IR luminosity of obscured star formation obtained from the MIPS 24 um flux. We find that at all redshifts the galaxies with higher masses have substantially lower specific star formation rates than lower mass galaxies. The average specific star formation rates increase with redshift, and the rate of incline is similar for all galaxies (roughly (1+z)^{n}, n = 5.0 +/- 0.4). It does not seem to be a strong function of galaxy mass. Using a subsample of galaxies with masses M_*> 10^11 M_sun, we measured the fraction of galaxies whose star formation is quenched. We consider a galaxy to be in quiescent mode when its specific star formation rate does not exceed 1/(3 x t_H), where t_H is the Hubble time. The fraction of quiescent galaxies defined as such decreases with redshift out to z ~ 1.8. We find that, at that redshift, 19 +/-9 % of the M_* > 10^11 M_sun sources would be considered quiescent according to our criterion.

We present K-band imaging of two ~30'x30' fields covered by the MUSYC Wide NIR Survey. The 1030 and 1255 fields were imaged with ISPI on the 4m Blanco telescope at CTIO to a 5 sigma point-source limiting depth of K~20 (Vega). Combining this data with the MUSYC Optical UBVRIz imaging, we created multi-band K-selected source catalogs for both fields. These catalogs, together with the MUSYC K-band catalog of the ECDF-S field, were used to select K<20 BzK galaxies over an area of 0.71 deg^2. This is the largest area ever surveyed for BzK galaxies. We present number counts, redshift distributions and stellar masses for our sample of 3261 BzK galaxies (2502 star-forming (sBzK) and 759 passively evolving (pBzK)), as well as reddening and star formation rate estimates for the star-forming BzK systems. We also present 2-point angular correlation functions and spatial correlation lengths for both sBzK and pBzK galaxies and show that previous estimates of the correlation function of these galaxies were affected by cosmic variance due to the small areas surveyed. We have measured correlation lengths r_0 of 8.89+/-2.03 Mpc and 10.82+/-1.72 Mpc for sBzK and pBzK galaxies respectively. This is the first reported measurement of the spatial correlation function of passive BzK galaxies. In the LambdaCDM scenario of galaxy formation, these correlation lengths at z~2 translate into minimum masses of ~4x10^{12} M_sun and ~9x10^{12} M_sun for the dark matter (DM) halos hosting sBzK and pBzK galaxies respectively. The clustering properties of the galaxies in our sample are consistent with them being the descendants of bright LBG at z~3, and the progenitors of present-day >1L* galaxies.

We present the mid-infrared colors of X-ray-detected AGN and explore mid-infrared selection criteria. Using a statistical matching technique, the likelihood ratio, over 900 IRAC counterparts were identified with a new MUSYC X-ray source catalog that includes ~1000 published X-ray sources in the Chandra Deep Field-South and Extended Chandra Deep Field-South. Most X-ray-selected AGN have IRAC spectral shapes consistent with power-law slopes, f_{nu} ~ nu^{alpha}, and display a wide range of colors, -2 < alpha < 2. Although X-ray sources typically fit to redder (more negative alpha) power-laws than non-X-ray detected galaxies, more than 50% do have flat or blue (galaxy-like) spectral shapes in the observed 3-8 micron band. Only a quarter of the X-ray selected AGN detected at 24 micron are well fit by featureless red power laws in the observed 3.6-24 micron, likely the subset of our sample whose infrared spectra are dominated by emission from the central AGN region. Most IRAC color-selection criteria fail to identify the majority of X-ray-selected AGN, finding only the more luminous AGN, the majority of which have broad emission lines. In deep surveys, these color-selection criteria select 10-20% of the entire galaxy population and miss many moderate luminosity AGN.

Using the Gemini Near-InfraRed Spectrograph (GNIRS), we have completed a near-infrared spectroscopic survey for K-bright galaxies at z~2.3, selected from the MUSYC survey. We derived spectroscopic redshifts from emission lines or from continuum features and shapes for all 36 observed galaxies. The continuum redshifts are driven by the Balmer/4000 Angstrom break, and have an uncertainty in dz/(1+z) of <0.019. We use this unique sample to determine, for the first time, how accurately redshifts and other properties of massive high-redshift galaxies can be determined from broadband photometric data alone. We find that the photometric redshifts of the galaxies in our sample have a systematic error of 0.08 and a random error of 0.13 in dz/(1+z). The systematic error can be reduced by using optimal templates and deep photometry; the random error, however, will be hard to reduce below 5%. The spectra lead to significantly improved constraints for stellar population parameters. For most quantities this improvement is about equally driven by the higher spectral resolution and by the much reduced redshift uncertainty. Properties such as the age, A_V, current star formation rate, and the star formation history are generally very poorly constrained with broadband data alone. Interestingly stellar masses and mass-to-light ratios are among the most stable parameters from broadband data. Nevertheless, photometric studies may overestimate the number of massive galaxies at 2<z<3, and thus underestimate the evolution of the stellar mass density. Finally,the spectroscopy supports our previous finding that red galaxies dominate the high-mass end of the galaxy population at z=2-3.

We present the first clustering results of X-ray selected AGN at z~3. Using Chandra X-ray imaging and UVR optical colors from MUSYC photometry in the ECDF-S field, we selected a sample of 58 z~3 AGN candidates. From the optical data we also selected 1385 LBG at 2.8<z< 3.8 with R<25.5. We performed auto-correlation and cross-correlation analyses, and here we present results for the clustering amplitudes and dark matter halo masses of each sample. For the LBG we find a correlation length of r_0,LBG = 6.7 +/- 0.5 Mpc, implying a bias value of 3.5 +/- 0.3 and dark matter (DM) halo masses of log(Mmin/Msun) = 11.8 +/- 0.1. The AGN-LBG cross-correlation yields r_0,AGN-LBG = 8.7 +/- 1.9 Mpc, implying for AGN at 2.8<z<3.8 a bias value of 5.5 +/- 2.0 and DM halo masses of log(Mmin/Msun) = 12.6 +0.5/-0.8. Evolution of dark matter halos in the Lambda CDM cosmology implies that today these z~3 AGN are found in high mass galaxies with a typical luminosity of 7+4/-2 L*.

We combine deep UBVRIzJK photometry from the MUSYC survey with redshifts from the COMBO-17 survey to study the rest-frame ultraviolet (UV) properties of 674 high-redshift (0.5<z<1) early-type galaxies, drawn from the Extended Chandra Deep Field South (E-CDFS). Galaxy morphologies are determined through visual inspection of Hubble Space Telescope (HST) images taken from the GEMS survey. We harness the sensitivity of the UV to young (<1 Gyrs old) stars to quantify the recent star formation history of the early-type population. We find compelling evidence that early-types of all luminosities form stars over the lifetime of the Universe, although the bulk of their star formation is already complete at high redshift. Luminous (-23<M(V)<-20.5) early-types form 10-15 percent of their mass after z=1, while their less luminous (M(V)>-20.5) counterparts form 30-60 percent of their mass in the same redshift range.

We investigate the stellar populations of a sample of 162 Lyman-alpha emitting galaxies (LAEs) at z = 3.1 in the Extended Chandra Deep Field South, using deep Spitzer IRAC data available from the GOODS and SIMPLE surveys to derive reliable stellar population estimates. We divide the LAEs according to their rest-frame near-IR luminosities into IRAC-detected and IRAC-undetected samples. About 70% of the LAEs are undetected in 3.6 micron down to [3.6] = 25.2 AB. Stacking analysis reveals that the average stellar population of the IRAC-undetected sample has an age of ~ 200 Myr and a mass of ~ 3x10^8 solar masses, consistent with the expectation that LAEs are mostly young and low-mass galaxies. On the other hand, the IRAC-detected LAEs are on average significantly older and more massive, with an average age > 1 Gyr and mass ~ 10^10 solar masses. Comparing the IRAC colors and magnitudes of the LAEs to z ~ 3 Lyman break galaxies (LBGs) shows that the IRAC-detected LAEs lie at the faint blue end of the LBG color-magnitude distribution, suggesting that IRAC-detected LAEs may be the low mass extension of the LBG population. We also present tentative evidence for a small fraction (~ 5%) of obscured AGN within the LAE sample. Our results suggest that LAEs posses a wide range of ages and masses. Additionally, the presence of evolved stellar populations inside LAEs suggests that the Lyman-alpha luminous phase of galaxies may either be a long-lasting or recurring phenomenon.

We studied the clustering properties and multiwavelength spectral energy distributions of a complete sample of 162 Ly Alpha-Emitting (LAE) galaxies at z=3.1 discovered in deep narrow-band MUSYC imaging of the Extended Chandra Deep Field South. LAEs were selected to have observed frame equivalent widths >80A and emission line fluxes >1.5E-17 erg/cm^2/s. Only 1% of our LAE sample appears to host AGN. The LAEs exhibit a moderate spatial correlation length of r_0=3.6+0.8-1.0 Mpc, corresponding to a bias factor b=1.7+0.3-0.4, which implies median dark matter halo masses of log10(M_med) = 10.9+0.5-0.9 M_sun. Comparing the number density of LAEs, (1.5+-0.3)E-3/Mpc^3, with the number density of these halos finds a mean halo occupation ~1-10%. The evolution of galaxy bias with redshift implies that most z=3.1 LAEs evolve into present-day galaxies with L<2.5L*, whereas other z>3 galaxy populations typically evolve into more massive galaxies. Halo merger trees show that z=0 descendants occupy halos with a wide range of masses, with a median descendant mass close to that of L*. Only 30% of LAEs have sufficient stellar mass (>~3E9 M_sun) to yield detections in deep Spitzer-IRAC imaging. A two-population SED fit to the stacked UBVRIzJK+[3.6,4.5,5.6,8.0]micron fluxes of the IRAC-undetected objects finds that the typical LAE has low stellar mass (1.0+0.6-0.4 E9 M_sun), moderate star formation rate (2+-1 M_sun/yr), a young component age of 20+30-10 Myr, and little dust (A_V<0.2). The best fit model has 20% of the mass in the young stellar component, but models without evolved stars are also allowed.

We combine deep UBVRIzJK photometry from the Multiwavelength Survey by Yale-Chile (MUSYC) with redshifts from the COMBO-17 survey to perform a large-scale study of the rest-frame ultraviolet (UV) properties of 674 high-redshift (0.5<z<1) early-type galaxies, drawn from the Extended Chandra Deep Field South (E-CDFS). Galaxy morphologies are determined through visual inspection of Hubble Space Telescope (HST) images taken from the GEMS survey. We harness the sensitivity of the UV to young (<1 Gyr old) stars to quantify the recent star formation history of early-type galaxies across a range of luminosities (-23.5 < M(V) < -18). Comparisons to simple stellar populations forming at high redshift indicate that only ~1.1 percent of early-types in this sample are consistent with purely passive ageing since z=2. Parametrising the recent star formation (RSF) in terms of the mass fraction of stars less than a Gyr old, we find that the early-type population as a whole shows a typical RSF between 5 and 13% in the redshift range 0.5<z<1. Early-types on the UV red sequence show RSF values less than 5% while the reddest early-types are virtually quiescent with RSF values of ~1%. We find compelling evidence that early-types of all luminosities form stars over the lifetime of the Universe, although the bulk of their star formation is already complete at high redshift. This tail-end of star formation is measurable and not negligible, with luminous (-23<M(V)<-20.5) early-types potentially forming 10-15% of their mass since z=1, with their less luminous (M(V)>-20.5) counterparts potentially forming 30-60 percent of their mass in the same redshift range. (abridged)

Using the Gemini Near-InfraRed Spectrograph (GNIRS), we have assembled a complete sample of 20 K-selected galaxies at 2.0<z<2.7 with high quality near-infrared spectra. As described in a previous paper, 9 of these 20 galaxies have strongly suppressed star formation and no detected emission lines. The present paper concerns the 11 galaxies with detected Halpha emission, and studies the origin of the line emission using the GNIRS spectra and follow-up observations with SINFONI on the VLT. Based on their [NII]/Halpha ratios, the spatial extent of the line emission and several other diagnostics, we infer that four of the eleven emission-line galaxies host narrow line active galactic nuclei (AGNs). The AGN host galaxies have stellar populations ranging from evolved to star-forming. Combining our sample with a UV-selected galaxy sample at the same redshift that spans a broader range in stellar mass, we find that black-hole accretion is more effective at the high-mass end of the galaxy distribution (~2.9x10^11 Msun) at z~2.3. Furthermore, by comparing our results with SDSS data, we show that the AGN activity in massive galaxies has decreased significantly between z~2.3 and z~0. AGNs with similar normalized accretion rates as those detected in our K-selected galaxies reside in less massive galaxies (~4.0x10^10 Msun) at low redshift. This is direct evidence for downsizing of AGN host galaxies. Finally, we speculate that the typical stellar mass-scale of the actively accreting AGN host galaxies, both at low and at high redshift, might be similar to the mass-scale at which star-forming galaxies seem to transform into red, passive systems.

We study the infrared (IR) properties of high-redshift galaxies using deep Spitzer 24, 70, and 160 micron data. Our primary interest is to improve the constraints on the total IR luminosities, L(IR), of these galaxies. We combine the Spitzer data in the southern Extended Chandra Deep Field with a K-band-selected galaxy sample and photometric redshifts from the Multiwavelength Survey by Yale-Chile. We used a stacking analysis to measure the average 70 and 160 micron flux densities of 1.5 < z < 2.5 galaxies as a function of 24 micron flux density, X-ray activity, and rest-frame near-IR color. Galaxies with 1.5 < z < 2.5 and S(24)=53-250 micro-Jy have L(IR) derived from their average 24-160 micron flux densities within factors of 2-3 of those derived from the 24 micron flux densities only. However, L(IR) derived from the average 24-160 micron flux densities for galaxies with S(24) > 250 micro-Jy and 1.5 < z < 2.5 are lower than those derived using only the 24 micron flux density by factors of 2-10. Galaxies with S(24) > 250 micro-Jy have S(70)/S(24) flux ratios comparable to sources with X-ray detections or red rest-frame IR colors, suggesting that warm dust possibly heated by AGN may contribute to the high 24 micron emission. Based on the average 24-160 micron flux densities, nearly all 24 micron-selected galaxies at 1.5 < z < 2.5 have L(IR) < 6 x 10^12 solar luminosities, which if attributed to star formation corresponds to < 1000 solar masses per year. This suggests that high redshift galaxies may have similar star formation efficiencies and feedback processes as local analogs. Objects with L(IR) > 6 x 10^12 solar luminosities are quite rare, with a surface density ~ 30 +/- 10 per sq. deg, corresponding to ~ 2 +/- 1 x 10^-6 Mpc^-3 over 1.5 < z < 2.5.

We publish the Keck/HIRES and Keck/ESI spectra that we have obtained during the first 10 years of Keck observatory operations. Our full sample includes 42 HIRES spectra and 39 ESI spectra along 65 unique sightlines providing abundance measurements on ~85 damped Lya systems. The normalized data can be downloaded from the journal or from our supporting website: http://www.ucolick.org/~xavier/DLA/. The database includes all of the sightlines that have been included in our papers on the chemical abundances, kinematics, and metallicities of the damped Lya systems. This data has also been used to argue for variations in the fine-structure constant. We present new chemical abundance measurements for 10 damped Lya systems and a summary table of high-resolution metallicity measurements (including values from the literature) for 153 damped Lya systems at z>1.6. We caution, however, that this metallicity sample (and all previous ones) is biased to higher N(HI) values than a random sample.

We present the rest-frame optical (B, V, and R band) luminosity functions (LFs) of galaxies at 2<z<3.5, measured from a K-selected sample constructed from the deep NIR MUSYC, the ultradeep FIRES, and the GOODS-CDFS. This sample is unique for its combination of area and range of luminosities. The faint-end slopes of the LFs at z>2 are consistent with those at z~0. The characteristic magnitudes are significantly brighter than the local values (e.g., ~1.2 mag in the R band), while the measured values for Phi_star are typically ~5 times smaller. The B-band luminosity density at z~2.3 is similar to the local value, and in the R band it is ~2 times smaller than the local value. We present the LF of Distant Red Galaxies (DRGs), which we compare to that of non-DRGs. While DRGs and non-DRGs are characterized by similar LFs at the bright end, the faint-end slope of the non-DRG LF is much steeper than that of DRGs. The contribution of DRGs to the global densities down to the faintest probed luminosities is 14%-25% in number and 22%-33% in luminosity. From the derived rest-frame U-V colors and stellar population synthesis models, we estimate the mass-to-light ratios (M/L) of the different subsamples. The M/L ratios of DRGs are ~5 times higher (in the R and V bands) than those of non-DRGs. The global stellar mass density at 2<z<3.5 appears to be dominated by DRGs, whose contribution is of order ~60%-80% of the global value. Qualitatively similar results are obtained when the population is split by rest-frame U-V color instead of observed J-K color.

We show that different stellar-mass estimation methods yield overall mass scales that disagree by factors up to ~2 for the z=0 galaxy population, and more importantly, relative mass scales that sometimes disagree by factors >~3 between distinct classes of galaxies (spiral/irregular types, classical E/S0s, and E/S0s whose colors reflect recent star formation). This comparison considers stellar mass estimates based on (a) two different calibrations of the correlation between K-band mass-to-light ratio and B-R color (Bell et al., Portinari et al.) and (b) detailed fitting of UBRJHK photometry and optical spectrophotometry using two different population synthesis models (Bruzual-Charlot, Maraston), with the same initial mass function in all cases. We also compare stellar+gas masses with dynamical masses. This analysis offers only weak arguments for preferring a particular stellar-mass estimation method, given the plausibility of real variations in dynamical properties and dark matter content. These results help to calibrate the systematic uncertainties inherent in mass-based evolutionary studies of galaxies, including comparisons of low and high redshift galaxies.

We study the clustering properties of K-selected galaxies at 2<z<3.5 using deep multiwavelength imaging in three fields from the MUSYC survey. These are the first measurements to probe the spatial correlation function of K-selected galaxies in this redshift range on large scales, allowing for robust conclusions about the dark matter halos that host these galaxies. K-selected galaxies with K<21 have a correlation length r_0 ~ 6 h^-1 Mpc, larger than typical values found for optically-selected galaxies. The correlation length does not depend on K-band magnitude, but it does increase strongly with color; the J-K>2.3 distant red galaxies (DRGs) have r_0 ~ 11 h^-1 Mpc. Contrary to findings for optically-selected galaxies, K-selected galaxies that are faint in the R-band cluster more strongly than brighter galaxies. These results suggest that a color-density relation was in place at z>2. Our results indicate that K-bright blue galaxies and K-bright red galaxies are fundamentally different, as they have different clustering properties. Using a simple model of one galaxy per halo, we infer halo masses ~ 5 x 10^12 M_sun for K<21 galaxies and ~ 2 x 10^13 M_sun for DRGs. A comparison of the observed space density of DRGs to the density of their host halos suggests large halo occupation numbers; however, this result is at odds with the lack of a strong small-scale excess in the angular correlation function. We find that the z=0 descendants of the galaxies considered here reside primarily in groups and clusters. [abridged]

We present deep near-infrared JHK imaging of four 10'x10' fields. The observations were carried out as part of the Multiwavelength Survey by Yale-Chile (MUSYC) with ISPI on the CTIO 4m telescope. The typical point source limiting depths are J~22.5, H~21.5, and K~21 (5sigma; Vega). The effective seeing in the final images is ~1.0". We combine these data with MUSYC UBVRIz imaging to create K-selected catalogs that are unique for their uniform size, depth, filter coverage, and image quality. We investigate the rest-frame optical colors and photometric redshifts of galaxies that are selected using common color selection techniques, including distant red galaxies (DRGs), star-forming and passive BzKs, and the rest-frame UV-selected BM, BX, and Lyman break galaxies (LBGs). These techniques are effective at isolating large samples of high redshift galaxies, but none provide complete or uniform samples across the targeted redshift ranges. The DRG and BM/BX/LBG criteria identify populations of red and blue galaxies, respectively, as they were designed to do. The star-forming BzKs have a very wide redshift distribution, a wide range of colors, and may include galaxies with very low specific star formation rates. In comparison, the passive BzKs are fewer in number, have a different distribution of K magnitudes, and have a somewhat different redshift distribution. By combining these color selection criteria, it appears possible to define a reasonably complete sample of galaxies to our flux limit over specific redshift ranges. However, the redshift dependence of both the completeness and sampled range of rest-frame colors poses an ultimate limit to the usefulness of these techniques.

We present the CXOCY J220132.8-320144 system, which is composed of an edge-on spiral galaxy at z=0.32 lensing a z=3.9 background quasar. Two images of the quasar are seen. The geometry of the system is favorable to separate the relative mass contribution of the disk and halo in the inner parts of the galaxy. We model the system with one elliptical mass component with the same ellipticity as the light distribution and manage to reproduce the quasar image positions and fluxes. We also model the system with two mass components, disk and halo. Again, we manage to reproduce the quasar image positions and fluxes. However, all models predict at least a third visible image close to the disk that is not seen in our images. We speculate that this is most likely due to extinction by the disk. We also measure the rotational velocity of the galaxy at 2.7 disk scale radius to be v_c=130 +/- 20 km/s from the [OII] emission lines. When adding the rotational velocity constraint to the models, we find that the contribution to the rotational velocity of the disk is likely to be equal to or larger than the contribution of the halo at this radius. The detection of the third image and a more accurate measurement of the rotational velocity would help to set tighter constraints on the mass distribution of this edge-on spiral galaxy.

We present first results of a spectroscopic survey targeting K-selected galaxies at z=2.0-2.7 using the GNIRS instrument on Gemini-South. We obtained near-infrared spectra with a wavelength coverage of 1.0-2.5 micron for 26 K-bright galaxies (K<19.7) selected from the MUSYC survey using photometric redshifts. We successfully derived spectroscopic redshifts for all 26 galaxies using rest-frame optical emission lines or the redshifted Balmer/4000 Angstrom break. Twenty galaxies have spectroscopic redshifts in the range 2.0<z<2.7, for which bright emission lines like Halpha and [OIII] fall in atmospheric windows. Surprisingly, we detected no emission lines for nine of these 20 galaxies. The median 2 sigma upper limit on the rest-frame equivalent width of Halpha for these nine galaxies is ~10 Angstrom. The stellar continuum emission of these same nine galaxies is best fitted by evolved stellar population models. The best-fit star formation rate (SFR) is zero for five out of nine galaxies, and consistent with zero within 1 sigma for the remaining four. Thus, both the Halpha measurements and the independent stellar continuum modeling imply that 45% of our K-selected galaxies are not forming stars intensely. This high fraction of galaxies without detected line emission and low SFRs may imply that the suppression of star formation in massive galaxies occurs at higher redshift than is predicted by current CDM galaxy formation models. However, obscured star formation may have been missed, and deep mid-infrared imaging is needed to clarify this situation.

We use 211 galaxy spectra from our survey for Lyman break galaxies (LBGs) associated with 11 damped Lyman alpha systems (DLAs) to measure the three-dimensional LBG auto-correlation and DLA-LBG cross-correlation functions with the primary goal of inferring the mass of DLAs at z~3. Conventional binning of the data while varying both r_0 and gamma parameters of the fiducial model of the correlation function xi(r)=(r/r_0)^{-gamma} resulted in the best fit values and 1 sigma uncertainties of r_0=2.65+/-0.48, gamma=1.55+/-0.40 for the LBG auto-correlation and r_0=3.32+/-1.25, gamma=1.74+/-0.36 for DLA-LBG cross-correlation function. To circumvent shortcomings found in binning small datasets, we perform a maximum likelihood analysis based on Poisson statistics. The best fit values and 1 sigma confidence levels from this analysis were found to be r_0=2.91(+1.0,-1.0), gamma=1.21(+0.6,-0.3) for the LBG auto-correlation and r_0=2.81(+1.4,-2.0), gamma=2.11(+1.3,-1.4) for the DLA-LBG cross-correlation function. We report a redshift spike of five LBGs with Delta z = 0.015 of the z=2.936 DLA in the PSS0808+5215 field and find that the DLA-LBG clustering signal survives when omitting this field from the analysis. Using the correlation functions measurements and uncertainties, we compute the z~3 LBG galaxy bias b_LBG to be 1.5<b_LBG<3 corresponding to an average halo mass of 10^(9.7)<M_LBG<10^(11.6) M_odot and the z~3 DLA galaxy bias b_DLA to be 1.3<b_DLA<4 corresponding to an average halo mass of 10^(9)<M_LBG<10^(12) M_odot. Lastly, two of the six QSOs discovered in this survey were found to lie within Delta z = 0.0125 of two of the survey DLAs. We estimate the probability of this occurring by chance is 1 in 940 and may indicate a possible relationship between the distribution of QSOs and DLAs at z~3.

As large optical quasar surveys for damped Lya become a reality and the study of star forming gas in the early Universe achieves statistical robustness, it is now vital to identify and quantify the sources of systematic error. Because the nature of optically-selected quasar surveys makes them vulnerable to dust obscuration, we have undertaken a radio-selected quasar survey for damped Lya systems to address this bias. We present the definition and results of this survey. We then combine our sample with the CORALS dataset to investigate the HI column density distribution function f(N) of damped Lya systems toward radio-selected quasars. We find that f(N) is well fit by a power-law f(N) = k_1 N^alpha_1, with log k_1 = 22.90 and alpha_1 = -2.18. This power-law is in excellent agreement with that of optically-selected samples at low N(HI), an important yet expected result given that obscuration should have negligible effect at these gas columns. However, because of the relatively small size of the radio-selected sample, 26 damped Lya systems in 119 quasars, f(N) is not well constrained at large N(HI) and the first moment of the HI distribution function, Omega_g, is, strictly speaking, a lower limit. The power-law is steep enough, however, that extrapolating it to higher column densities implies only a modest, logarithmic increase in Omega_g. The radio-selected value of Omega_g = 1.15 x 10^-3, agrees well with the results of optically-selected surveys. While our results indicate that dust obscuration is likely not a major issue for surveys of damped Lya systems, we estimate that a radio-selected sample of approximately 100 damped Lya systems will be required to obtain the precision necessary to absolutely confirm an absence of dust bias.

We selected 40 candidate Lyman Alpha Emitting galaxies (LAEs) at z ~=3.1 with observed frame equivalent widths >150A and inferred emission line fluxes >2.5x10^-17 ergs/cm^2/s from deep narrow-band and broad-band MUSYC images of the Extended Chandra Deep Field South. Covering 992 sq. arcmin, this is the largest ``blank field'' surveyed for LAEs at z ~3, allowing an improved estimate of the space density of this population of 3+-1x10^-4 h_70^3/Mpc^3. Spectroscopic follow-up of 23 candidates yielded 18 redshifts, all at z ~=3.1. Over 80% of the LAEs are dimmer in continuum magnitude than the typical Lyman break galaxy spectroscopic limit of R= 25.5 (AB), with a median continuum magnitude R ~=27 and very blue continuum colors, (V-z) ~=0. Over 80% of the LAEs have the right UVR colors to be selected as Lyman break galaxies, but only 10% also have R<=25.5. Stacking the UBVRIzJK fluxes reveals that LAEs have stellar masses ~=5x10^8 h_70^-2 M_sun and minimal dust extinction, A_V < ~ 0.1. Inferred star formation rates are ~=6 h_70^-2 M_sun/yr, yielding a cosmic star formation rate density of 2x10^-3 h_70 M_sun/yr/Mpc^3. None of our LAE candidates show evidence for rest-frame emission line equivalent widths EW_rest>240A which might imply a non-standard IMF. One candidate is detected by Chandra, implying an AGN fraction of 2+-2% for LAE candidate samples. In summary, LAEs at z ~ 3 have rapid star formation, low stellar mass, little dust obscuration and no evidence for a substantial AGN component.

The Multiwavelength Survey by Yale-Chile (MUSYC) is optimized for the study of galaxies at z=3, AGN demographics, and Galactic structure. MUSYC consists of deep optical (UBVRIz') and near-infrared (JHK) imaging of four fields on the sky covering 1.2 square degrees to AB limiting depths of U,B,V,R=26 and K=22. Our optical catalog contains 277,341 objects detected in BVR images with median seeing of 0.9''. Satellite coverage of our fields includes Chandra, XMM, GALEX, HST-ACS, and Spitzer, with the Extended Chandra Deep Field South imaged at all of these wavelengths plus the radio, making it the premier multiwavelength field on the sky. Detailed follow-up spectroscopy is being performed with VLT+VIMOS, Magellan+IMACS and Gemini+GNIRS. MUSYC provides ideal supporting data for surveys with Z-Machines and rich target lists for spectroscopy with ALMA. We are conducting a census of protogalaxies at redshift three (Lyman break galaxies, Lyman-alpha emitters, Distant red galaxies, Sub-millimeter galaxies and AGN) in order to separate physical properties from selection effects. We discuss measurements of the dark matter halo masses and halo occupation numbers of these populations and of the total cosmic star formation rate at z=3. MUSYC publications and data releases are available at http://www.astro.yale.edu/MUSYC .

We use near-infrared (NIR) spectroscopy obtained with GNIRS on Gemini, NIRSPEC on KECK, and ISAAC on the VLT to study the rest-frame optical continua of three `Distant Red Galaxies' (having Js - Ks > 2.3) at z>2. All three galaxy spectra show the Balmer/4000 Angstrom break in the rest-frame optical. The spectra allow us to determine spectroscopic redshifts from the continuum with an estimated accuracy dz/(1+z) ~ 0.001-0.04. These redshifts agree well with the emission line redshifts for the 2 galaxies with Halpha emission. This technique is particularly important for galaxies that are faint in the rest-frame UV, as they are underrepresented in high redshift samples selected in optical surveys and are too faint for optical spectroscopy. Furthermore, we use the break, continuum shape, and equivalent width of Halpha together with evolutionary synthesis models to constrain the age, star formation timescale, dust content, stellar mass and star formation rate of the galaxies. Inclusion of the NIR spectra in the stellar population fits greatly reduces the range of possible solutions for stellar population properties. We find that the stellar populations differ greatly among the three galaxies, ranging from a young dusty starburst with a small break and strong emission lines to an evolved galaxy with a strong break and no detected line emission. The dusty starburst galaxy has an age of 0.3 Gyr and a stellar mass of 1*10^11 Msun. The spectra of the two most evolved galaxies imply ages of 1.3-1.4 Gyr and stellar masses of 4*10^11 Msun. The large range of properties seen in these galaxies strengthens our previous much more uncertain results from broadband photometry. Larger samples are required to determine the relative frequency of dusty starbursts and (nearly) passively evolving galaxies at z~2.5.

The Extended Chandra Deep Field-South (ECDFS) survey consists of 4 Chandra ACIS-I pointings and covers $\approx$ 1100 square arcminutes ($\approx$ 0.3 deg$^2$) centered on the original CDF-S field to a depth of approximately 228 ks. This is the largest Chandra survey ever conducted at such depth, and only one XMM-Newton survey reaches a lower flux limit in the hard 2.0--8.0 keV band. We detect 651 unique sources -- 587 using a conservative source detection threshold and 64 using a lower source detection threshold. These are presented as two separate catalogs. Of the 651 total sources, 561 are detected in the full 0.5--8.0 keV band, 529 in the soft 0.5--2.0 keV band, and 335 in the hard 2.0--8.0 keV band. For point sources near the aim point, the limiting fluxes are approximately $1.7 \times 10^{-16}$ $\rm{erg cm^{-2} s^{-1}}$ and $3.9 \times 10^{-16}$ $\rm{erg cm^{-2} s^{-1}}$ in the 0.5--2.0 keV and 2.0--8.0 keV bands, respectively. Using simulations, we determine the catalog completeness as a function of flux and assess uncertainties in the derived fluxes due to incomplete spectral information. We present the differential and cumulative flux distributions, which are in good agreement with the number counts from previous deep X-ray surveys and with the predictions from an AGN population synthesis model that can explain the X-ray background. In general, fainter sources have harder X-ray spectra, consistent with the hypothesis that these sources are mainly obscured AGN.

Using the deep multi-wavelength MUSYC, GOODS, and FIRES surveys we construct a stellar mass-limited sample of galaxies at 2<z<3. The sample comprises 294 galaxies with M>10^11 Solar masses distributed over four independent fields with a total area of almost 400 sq arcmin. The mean number density of massive galaxies in this redshift range is (2.2+-0.6) x 10^-4 Mpc^-3. We present median values and 25th and 75th percentiles for the distributions of observed R mags, observed J-K colors, and rest-frame UV continuum slopes, M/L(V) ratios, and U-V colors. The galaxies show a large range in all these properties. The ``median galaxy'' is faint in the observer's optical (R=25.9), red in the observed near-IR (J-K=2.48), has a rest-frame UV spectrum which is relatively flat (beta=-0.4), and rest-frame optical colors resembling those of nearby spiral galaxies (U-V=0.62). We determine which galaxies would be selected as Lyman break galaxies (LBGs) or Distant Red Galaxies (DRGs, having J-K>2.3) in this mass-limited sample. By number DRGs make up 69% of the sample and LBGs 20%, with a small amount of overlap. By mass DRGs make up 77% and LBGs 17%. Neither technique provides a representative sample of massive galaxies at 2<z<3 as they only sample the extremes of the population. As we show here, multi-wavelength surveys with high quality photometry are essential for an unbiased census of massive galaxies in the early Universe. The main uncertainty in this analysis is our reliance on photometric redshifts; confirmation of the results presented here requires extensive near-infrared spectroscopy of optically-faint samples.

I summarize current knowledge of galaxy formation with emphasis on the initial conditions provided by the Lambda CDM cosmology, integral constraints from cosmological quantities, and the demographics of high-redshift protogalaxies. Tables are provided summarizing the number density, star formation rate and stellar mass per object, cosmic star formation rate and stellar mass densities, clustering length and typical dark matter halo masses for Lyman break galaxies, Lyman alpha emitting galaxies, Distant red galaxies, Sub-millimeter galaxies, and Damped Lyman alpha absorption systems. I also discuss five key unsolved problems in galaxy formation and prognosticate advances that the near future will bring.

We present Spitzer 24micron imaging of 1.5 < z < 2.5 Distant Red Galaxies (DRGs) in the 10arcmin by 10arcmin Extended Hubble Deep Field South of the Multiwavelength Survey by Yale-Chile. We detect 65% of the DRGs with K_AB < 23.2 mag at S_24micron > 40uJy, and conclude that the bulk of the DRG population are dusty active galaxies. A mid-infrared (MIR) color analysis with IRAC data suggests that the MIR fluxes are not dominated by buried AGN, and we interpret the high detection rate as evidence for a high average star formation rate of <SFR> = 130+/-30 M/yr. From this, we infer that DRGs are important contributors to the cosmic star formation rate density at z ~ 2, at a level of \~0.02 M/yr/Mpc^3 to our completeness limit of K_AB = 22.9 mag

We present the first spectroscopic measurement of the spatial cross-correlation function between damped Lyman alpha systems (DLAs) and Lyman break galaxies (LBGs). We obtained deep u'BVRI images of nine QSO fields with 11 known z ~ 3 DLAs and spectroscopically confirmed 211 R < 25.5 photometrically selected z > 2 LBGs. We find strong evidence for an overdensity of LBGs near DLAs versus random, the results of which are similar to that of LBGs near other LBGs. A maximum likelihood cross-correlation analysis found the best fit correlation length value of r_0 = 2.9^(+1.4)_(-1.5) h^(-1)Mpc using a fixed value of gamma = 1.6. The implications of the DLA-LBG clustering amplitude on the average dark matter halo mass of DLAs are discussed.

The Extended Chandra Deep Field-South (E-CDF-S) survey consists of 4 Chandra ACIS-I pointings and covers ~1100 square arcminutes (~0.3 deg^2) centered on the original CDF-S field to a depth of approximately 228 ks (PI: Niel Brandt; Lehmer et al, 2005). This is the largest Chandra survey ever conducted at such depth. In our analysis (Virani et al, 2005), we detect 651 unique sources -- 587 using a conservative source detection limit and 64 using a lower source detection limit. Of these 651 sources, 561 are detected in the full 0.5--8.0 keV band, 529 in the soft 0.5--2.0 keV band, and 335 in the hard 2.0--8.0 keV band. For point sources near the aim point, the limiting fluxes are approximately $1.7 \times 10^{-16}$ $\rm{erg cm^{-2} s^{-1}}$ and $3.9 \times 10^{-16}$ $\rm{erg cm^{-2} s^{-1}}$ in the 0.5--2.0 keV and 2.0--8.0 keV bands, respectively. We present the optical properties of these X-ray sources, specifically the R-band magnitude distribution and a preliminary spectroscopic redshift distribution. One exciting result is the discovery of 7 new Extreme X-ray-to-Optical flux ratio objects (EXOs) found in the E-CDF-S field.

Observations of damped Lyman alpha systems offer a unique window on the neutral-gas reservoirs that gave rise to galaxies at high redshifts. This review focuses on critical properties such as the H I and metal content of the gas and on independent evidence for star formation. Together, these provide an emerging picture of gravitationally bound objects in which accretion of gas from the IGM replenishes gas consumed by star formation. Other properties such as dust content, molecular content, ionized-gas content, gas kinematics, and galaxy identifications are also reviewed. These properties point to a multiphase ISM in which radiative and hydrodynamic feedback processes are present. Numerical simulations and other types of models used to describe damped Lyman alpha systems within the context of galaxy formation are also discussed.

We present UBVRIz' optical images of the 0.32 square degree Extended Hubble Deep Field South. This is one of four fields comprising the MUSYC survey, which is optimized for the study of galaxies at z=3, AGN demographics, and Galactic structure. We calculate corrected aperture photometry and its uncertainties and find through tests that these provide a significant improvement upon standard techniques. Our photometric catalog of 62968 objects is complete to a total magnitude of R_AB=25. We select z=3 Lyman break galaxy (LBG) candidates from their UVR colors and find a sky surface density of 1.4/sq. arcmin and an angular correlation function w(theta) = 2.3+-1.0 theta^{-0.8}. (Abridged)

MUSYC (MUltiwavelength Survey by Yale and Chile), and its predecessor, the CYDER (Cal'an Yale Deep Extragalactic Research) survey are both deep multiwavelength surveys mainly aiming at extragalactic issues; however they were both designed with a strong Galactic component in mind. Therefore they consist not only of multipassband imaging (UBVRIzJK) and wide field multi object spectroscopy - it is also foreseen to create a multi-epoch dataset, thus allowing proper motions to be derived using data of 1 and 4-5 years baseline. This enables us to identify fast moving objects, and derive tangential velocities of nearby and intrinsically faint objects, such as White Dwarfs (especially cool WDs) and Brown Dwarfs (BDs). For the former, we will be able to study the faint end of the halo-field luminosity function and determine the fraction of baryonic dark matter consisting of WDs. We will be able to put constraints on the volume density of BDs, and determine whether this density actually meets expectations. Here we give a report of work in progress.

We present a survey for z~3 Lyman break galaxies (LBGs) associated with damped Lyman alpha systems (DLAs) with the primary purpose of determining the DLA-LBG cross-correlation. This paper describes the acquisition and analysis of imaging and spectroscopic data of 9 quasar fields having 11 known z~3 DLAs covering an area of 465 arcmin^2. Using deep u'BVRI images, 796 LBG candidates to an apparent R_AB magnitude of 25.5 were photometrically selected from 17,343 sources detected in the field. Spectroscopic observations of 529 LBG candidates using Keck LRIS yielded 339 redshifts. We have conservatively identified 211 z>2 objects with <z>=3.02+/-0.32. We discuss our method of z~3 LBG identification and present a model of the u'BVRI photometric selection function. We use the 339 spectra to evaluate our u'BVRI z~3 Lyman break photometric selection technique.

We present the first results from the Cal\'an-Yale Deep Extragalactic Research (CYDER) survey. The main goal of this survey is to study serendipitous X-ray sources detected by Chandra in an intermediate flux range ($10^{-15}-10^{-12}$ ergs s$^{-1}$) that comprises most of the X-ray background. 267 X-ray sources spread over 5 archived fields were detected. The $\log N-\log S$ distribution obtained for this sample is consistent with the results of other surveys. Deep $V$ and $I$ images were taken of these fields in order to calculate X-ray-to-optical flux ratios. Identifications and redshifts were obtained for 106 sources using optical spectroscopy from 8-m class telescopes to reach the optically faintest sources, to the same level as deeper X-ray fields like the Chandra Deep Fields, showing that the nature of sources detected depends mostly on the optical limit for spectroscopy. In general, sources optically classified as obscured Active Galactic Nuclei (AGNs) have redder optical colors than unobscured AGN. A rough correlation between $f_X/f_{\rm opt}$ and hard X-ray luminosity was found for obscured AGN confirming the prediction by existing models that in obscured AGN the optical light is completely dominated by the host galaxy. The previously claimed decrease of the obscured to unobscured AGN ratio with increasing X-ray luminosity is observed. However, this correlation can be explained as a selection effect caused by the lower optical flux of obscured AGN. Comparison between the observed $N_H$ distribution and predictions by existing models shows that the sample appears complete up to $N_H<3\times 10^{22}$ cm$^{-2}$, while for more obscured sources incompleteness plays an important role in the observed obscured to unobscured AGN ratio.

We investigate the heat source of the neutral gas comprising DLAs. Unlike the Lya forest, where the extragalactic background radiation field ionizes and heats the gas, we find that grain photoelectric heating by the FUV background is not sufficient to balance the C II 158um cooling rate inferred from DLAs. In these systems, a local energy source is required. We show that in the case of the z=1.919 DLA toward Q2206-19, the local source is FUV emission from the associated galaxy found by Moller et al (2002): the mean intensity inferred from photometry is in good agreement with the intensity Jnu required to explain the cooling rate. The FUV mean intensity predicted for a cold neutral medium (CNM) model, Jnu=(1.7+2.7-1.0)x10^(-18) cgs (95% c.l.), is the largest expected from our CII* study of 45 DLAs. This may explain why this is the only confirmed DLA yet detected in emission at z>1.9. We argue that in most DLAs with detected CII* absorption, Jnu is between 10^{-19} and 10^{-18} and heats the gas which is a CNM. By contrast, in most DLAs with upper limits on CII* absorption the gas is a warm neutral medium (WNM). Surprisingly, the upper limits are compatible with the same range of Jnu values suggesting the majority of DLAs are heated by radiation fields generated by a limited range of star formation rates per unit H I area, between 10^{-3} and 10^{-2} Msol/kpc^2. We also show that CII* absorption is unlikely to arise in gas that is ionized.

Cosmological observations strongly suggest the presence of dark energy which comprises the majority of the current energy density of the universe. The equation of state relating the pressure and energy density of this dark energy, p = w \rho, appears to have w close to -1, with most analyses preferring w<-1 when these values are considered as part of the parameter space. If w<-1 is the future behavior of the dark energy, the scale factor, expansion rate, and energy density of the universe will diverge in a finite time as its apparent expansion age tends to zero. We hypothesize that w>-1 is restored at a late enough time, perhaps due to a phase transition of the dark energy, and show that this produces conditions observationally indistinguishable from a Hot Big Bang. This process of rejuvenation may have occurred in the past, making our universe much older than it appears and eliminating the Big Bang singularity.

We present the discovery of a high redshift, X-ray selected AGN by the Calan-Yale Deep Extragalactic Research (CYDER) survey: CXOCY J033716.7-050153, located at $z=4.61$, the second high redshift AGN discovered by this survey. Here, we present its optical, near-IR and X-ray properties and compare it with other optical and X-ray selected high redshift AGN. The optical luminosity of this object is significantly lower than most optically selected high redshift quasars. It also has a lower rest frame UV to X-ray emission ratio than most known quasars at this redshift. This mild deviation can be explained either by dust obscuring the UV radiation of a normal radio quiet AGN emitting at 10% of its Eddington luminosity or because this is intrinsically a low luminosity radio loud AGN, with a super-massive Black Hole of $\sim 10^8M_\sun$ emitting at 1% of its Eddington luminosity. Deep radio observations can discriminate between these two hypotheses.

This paper presents chemical abundance measurements for 37 damped Lya systems at z>2.5 observed with the Echellette Spectrograph and Imager on the KeckII telescope. We measure the HI column densities of these damped systems through Voigt profile fits to their Lya profiles and we implement the apparent optical depth method to determine ionic column densities. Figures and tables of all relevant data are presented. A full analysis of the chemical enrichment history described by these observations will be presented in a future paper. This dataset is also valuable for efficiently planning future echelle observations and for rough abundance pattern analyses. We aim to make this entire data set public within three years of this publication.

We present accurate metallicity measurements for 121 damped Lya systems at 0.5<z<5 including ~50 new measurements from our recently published Echellette Spectrograph and Imager surveys. This dataset is analysed to determine the age-metallicity relation of neutral gas in the universe. Contrary to previous datasets this sample shows statistically significant evolution in the mean metallicity. The best linear fit rate to metallicity vs. redshift is -0.26 +/- 0.07 dex corresponding to approximately a factor of 2 every Gyr at z=3. The DLA continue to maintain a floor in metallicity of ~1/700 solar independent of observational effects. This metallicity threshold limits the prevalence of primordial gas in high redshift galaxies and stresses the correspondence between damped systems and star formation (i.e. galaxy formation). This floor is significantly offset from the metallicity of the Lya forest and therefore we consider it to be more related to active star formation within these galaxies than scenarios of enrichment in the very early universe. Finally, we comment on an apparent 'missing metals problem': the mean metallicity of the damped systems is ~10x lower than the value expected from their observed star formation history. This problem is evident in current theoretical treatments of chemical evolution and galaxy formation; it may indicate a serious flaw in our understanding of the interplay between star formation and metal production.

Starting from the SFR per unit physical area, determined for DLAs using the C II* method, we obtain the SFR per unit comoving volume at $z$ $\approx$ 3, and find that it agrees with that for the Lyman Break Galaxies (LBGs). Though the mass of produced stars indicated by the SFRs is consistent with the current densities of known stellar populations, the mass of metals produced by $z$=2.5 is 30 times larger than detected in absorption in DLAs. The most likely solutions to this ``missing metals'' problem is that star formation occurs in compact bulge regions. We search for evidence of feedback and find no correlations between the SFR per unit area and N(H I), but possible correlations between SFR per unit area and low-ion velocity width and SFR per unit area and metal abundance. We show that (a) the correlation between cooling rate and dust-to-gas ratio is positive evidence for grain photoelectric heating, (b) the CMB does not significantly populate the C II excited fine-structure states, and (c) the ratio of CII* to resonance-line optical depths is a sensitive probe of the multi-phase structure of the DLA gas. We address recent arguments that DLAs are comprised only of WNM gas, and show them to be inconclusive. Despite the rough agreement between SFR per unit comoving volume for DLAs and LBGs, current evidence indicates these are distinct populations.

We describe a technique that for the first time measures star formation rates (SFRs) in damped Lyman alpha systems(DLAs) directly. We assume that massive stars form in DLAs, and that the FUV radiation they emit heats the gas by the grain photoelectric mechanism. We infer the heating rate from the cooling rate measured by the strength of CII* 1335.7 absorption. Since the heating rate is proportional to the dust-to-gas ratio and the SFR per unit area, we deduce the SFR per unit area for DLAs in which both quantities have been measured. We consider models in which the the dust comprises carbonaceous or silicate grains. We present two-phase models where the cold neutral medium (CNM) and warm neutral medium (WNM) are in pressuer equilibrium. In the CNM model the the sightline goes throught the CNM and WNM, while in the WNM model it goes only through the WNM. Since the grain photoelectric heating efficiency is at least 10 times higher in the CNM than in the WNM, CII* absorption mainly arises in the CNM in the CNM model. But in the WNM model all of the CII* absorption arises in the WNM. We use CII* absorption lines to derive the SFR per unit area for a sample of ~ 30 DLAs in which the dust-to-gas ratio has been inferred from element depletion patterns. We show that the resulting SFR per unit area corresponds to an average over the star forming volume of galaxy hosting the DLA rather than to local star formation along the line of sight. We find the average SFR per unit area and equals10$^{-2.2}$ M$_{\odot}$yr$^{-1}$kpc$^{-2}$ for the CNM model and 10$^{-1.3}$ M$_{\odot}$yr$^{-1}$kpc$^{-2}$ for the WNM model. The SFR per unit area in the CNM solution is similar to that measured in the Milky Way ISM.

We present spectra of the optical transient (OT) associated with GRB 021004. The spectra show a blue continuum with superposed absorption features and one emission line. We confirm two intervening metal-line systems at z = 1.380 and z = 1.602 and one very strong absorption system at a redshift of z = 2.323. Ly_alpha emission is also seen at this redshift. While the spectrum of the OT overall cannot be simply described with a power law, the spectral index over the range 5500-8850 A is steep, F_nu ~ nu^(-0.96 +/- 0.03). Comparison of spectra from multiple epochs shows a distinct color evolution with the OT becoming redder with time over the first three days. This is the first clear example of color change in an OT detected spectroscopically.

This paper describes the acquisition and analysis of imaging data for the identification of galaxies associated with z~4 damped Lya systems. We present deep BRI images of three fields known to contain four z~4 damped systems. We discuss the reduction and calibration of the data, detail the color criteria used to identify z~4 galaxies, and present a photometric redshift analysis to complement the color selection. We have found no galaxy candidates closer to the QSO than 7'' which could be responsible for the damped Lya systems. Assuming that at least one of the galaxies is not directly beneath the QSO, we set an upper limit on this damped Lya system of L < L*/4. Finally, we have established a web site to release these imaging data to the public.

We present new chemical abundance measurements of 16 damped Lya systems at z>1.5 and update our previous abundance analyses. The entire database presented here was derived from HIRES observations on the Keck I telescope, reduced with the same software package, and analysed with identical techniques. Altogether, we present a large, homogeneous database of chemical abundance measurements for protogalaxies in the early universe, ideal for studying a number of important aspects of galaxy formation. In addition, we have established an online directory for this database and will continuously update the results.

We report the first observations of the clustering environment of damped Lyman alpha absorption systems at z=4. Color selection and photometric redshifts were used to select 44 candidate Lyman-break galaxies brighter than I_AB=25.5 from deep BRI images of the 35 sq. arcmin field containing the quasar BR 0951-04. Multislit spectroscopy of 35 candidate galaxies was performed and 8 of these candidates have been confirmed as z>3.5 Lyman-break galaxies. With only BRI photometry, the photometric redshifts are quite accurate for the spectroscopically confirmed galaxies but have a high rate of misclassification due to color degeneracies between Lyman-break galaxies and low-redshift ellipticals. Both of the z>3.5 galaxies found within 15'' of the quasar line-of-sight appear to be causing absorption systems in the quasar spectrum. We use a battery of statistical tests to look for clustering in the redshift histogram of the z>3.5 galaxies but do not find measurable clustering of these Lyman-break galaxies with the damped Lyman alpha absorbers. With a larger sample of galaxies, our method should determine the cross-correlation between these objects, which probes the bias and hence the mass of the damped Lyman alpha absorbers.

A conflict has been reported between the baryon density inferred from deuterium observations and that found from recent CMB observations by BOOMERanG and MAXIMA. Despite the flurry of papers that attempt to resolve this conflict by adding new physics to the early universe, we will show that it can instead be resolved via a more careful usage of statistics. Indeed, the Bayesian analyses that produce this conflict are by their nature poorly suited for drawing this type of conclusion. A properly defined frequentist analysis can address this question directly and appears not to find a conflict. Finally, a conservative accounting of systematic uncertainties in measuring the deuterium abundance could reduce what is nominally a 3 sigma conflict to 1 sigma.

We present the first results from an ongoing survey to discover and measure the metallicity of z>3 damped Lya systems with the Echellette Spectrograph and Imager (ESI) on the Keck II telescope. Our motivation arises from a recent study on the damped Lya systems suggesting only mild evolution in the cosmic metallicity from z~2 to 4. The Echellette Spectrograph and Imager, which provides two complementary spectroscopic modes, is the ideal instrument for a z>3 damped Lya survey. We describe our observing strategy and report on the discovery and analysis of 5 new z>3 damped Lya systems acquired in a single night of observing. These observations further support the principal conclusions of the previous study: (1) the cosmic metallicity in neutral gas inferred from the damped Lya systems does not evolve significantly from z~2 to 4; (2) the unweighted metallicity exhibits a statistically significant decrease with increasing redshift; and (3) not a single damped Lya system has a metallicity below [Fe/H]=-3. We discuss the implications of these results and comment on recent theoretical studies which attempt to explain the observations.

We show how observations of the evolution of the galaxy cluster number abundance can be used to constrain primordial non-Gaussianity in the universe. We carry out a maximum likelihood analysis incorporating a number of current datasets and accounting for a wide range of sources of systematic error. Under the assumption of Gaussianity, the current data prefer a universe with matter density $\Omega_m\simeq 0.3$ and are inconsistent with $\Omega_m=1$ at the $2\sigma$ level. If we assume $\Omega_m=1$, the predicted degree of cluster evolution is consistent with the data for non-Gaussian models where the primordial fluctuations have at least two times as many peaks of height $3\sigma$ or more as a Gaussian distribution does. These results are robust to almost all sources of systematic error considered: in particular, the $\Omega_m=1$ Gaussian case can only be reconciled with the data if a number of systematic effects conspire to modify the analysis in the right direction. Given an independent measurement of $\Omega_m$, the techniques described here represent a powerful tool with which to constrain non-Gaussianity in the primordial universe, independent of specific details of the non-Gaussian physics. We discuss the prospects and strategies for improving the constraints with future observations.

The claim that large-scale structure data independently prefers the Lambda Cold Dark Matter model is a myth. However, an updated compilation of large-scale structure observations cannot rule out Lambda CDM at 95% confidence. We explore the possibility of improving the model by adding Hot Dark Matter but the fit becomes worse; this allows us to set limits on the neutrino mass.

Massive neutrinos have a detectable effect on cosmological structure formation, in particular on the large-scale distribution of galaxies. Adding Hot Dark Matter to the now-standard Lambda CDM model leads to a worse fit to large-scale structure and CMB anisotropy data. This results in a limit on the mass of the most massive neutrino of 4 eV, assuming a power-law primordial power spectrum.

We combine detections of anisotropy in the Cosmic Microwave Background radiation with observations of inhomogeneity in the large-scale distribution of galaxies to test the predictions of models of cosmological structure formation. This combination probes spatial scales varying by three orders of magnitude, including a significant region where the two types of data overlap. We examine Cold Dark Matter models with adiabatic density perturbations, isocurvature models, and a topological defects model. We set upper limits on the neutrino mass and find the primordial power spectrum needed to reconcile an apparent disagreement between structure formation observations and direct observations of cosmological parameters. Present and future observations of Cosmic Microwave Background anisotropy suffer from foreground contamination. We develop detailed predictions for microwave emission from radio and infrared-bright galaxies and the Sunyaev-Zeldovich effect from clusters. We present realistic simulations of the microwave sky, produced as part of the `WOMBAT Challenge' exercise, and introduce a pixel-space method for subtracting foreground contamination which can be tested on these simulations.

We summarize the theoretical and observational status of the study of the Cosmic Microwave Background radiation. Its thermodynamic spectrum is a robust prediction of the Hot Big Bang cosmology and has been confirmed observationally. There are now 76 observations of Cosmic Microwave Background anisotropy, which we present in a table with references. We discuss the theoretical origins of these anisotropies and explain the standard jargon associated with their observation.

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

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

The shape of the primordial fluctuation spectrum is probed by cosmic microwave background fluctuations which measure density fluctuations at z~1000 on scales of hundreds of Mpc and from galaxy redshift surveys, which measure structure at low redshift out to several hundred Mpc. The currently acceptable library of cosmological models is inadequate to account for the current data, and more exotic models must be sought. New data sets such as SDSS and 2DF are urgently needed to verify whether the shape discrepancies in P(k) will persist.

We estimate the contribution of extragalactic radio sources to fluctuations in sky temperature over the range of frequencies (10-300 GHz) used for Cosmic Microwave Background (CMB) anisotropy measurements. CMB anisotropy observations at high resolution and low frequencies are especially sensitive to this foreground. We have compiled a catalog of 2207 bright radio sources, including 758 sources with flux measurements at 90 GHz. We develop a method to extrapolate the source spectra and predict skymaps of extragalactic radio sources at instrument resolutions of 10 arcmin to 10 degrees FWHM. Our results indicate that the brightest sources will dominate microwave anisotropy for a wide range of resolutions and frequencies. Our skymaps predict the location and flux of the brightest radio sources at each frequency, making it straightforward to develop a template for masking the pixels containing them. This masking should be sufficient to protect high resolution CMB anisotropy observations from unacceptable radio source confusion.

Applying basic physical principles to recent observational results, we derive upper and lower limits on microwave anisotropy from point sources over the range of frequencies 10-1000 GHz. We examine the level of noise in the observations as a possible indication of source confusion at subarcminute scales. We also derive an upper limit on microwave anisotropy caused by the sources responsible for the Far-Infrared Background radiation detected in FIRAS data. Our upper limit on point source confusion of DeltaT/T=10^{-5} for a 10' beam at 100 GHz would cause severe foreground contamination for CMB anisotropy observations, although the actual contamination level is probably much lower. This upper limit constrains the long-feared possibility of an undetected population of sources with emission peaking near 100 GHz. Source detections closer to 100 GHz are needed to improve our knowledge of galaxy evolution at high redshift and to predict the level of point source confusion.

Cosmological models predict transfer functions by which primordial density perturbations develop into CMB anisotropy and Large-Scale Structure. We use the current set of observations to reconstruct the primordial power spectrum for standard CDM, CDM with a cosmological constant, open CDM, and standard CDM with a high baryon content.

The combination of detections of anisotropy in the Cosmic Microwave Background radiation and observations of the large-scale distribution of galaxies probes the primordial density fluctuations of the universe on spatial scales varying by three orders of magnitude. These data are found to be inconsistent with the predictions of several popular cosmological models. Agreement between the data and the Cold + Hot Dark Matter model, however, suggests that a significant fraction of the matter in the universe may consist of massive neutrinos.

We consider constraints on the amplitude of mass fluctuations in the universe, sigma_8, derived from two simple observations: the present number density of clusters and the amplitude of their correlation function. Allowing for the possibility that the primordial fluctuations are non-gaussian introduces a degeneracy in the value of sigma_8 preferred by each of these constraints. However, when the constraints are taken together this degeneracy is broken, yielding a precise determination of sigma_8 and the degree of non-gaussianity for a given background cosmology. For a flat, Omega_m=1 universe with a power spectrum parameterized by a CDM shape parameter Gamma=0.2, we find that the perturbations are consistent with a gaussian distribution with sigma_8=0.49(+0.08-0.07) (95% limits). For some popular choices of background model, including the favored low matter density models, the hypothesis that the primordial fluctuations are gaussian is ruled out with a high degree of confidence.

We estimate the level of confusion to Cosmic Microwave Background anisotropy measurements caused by extragalactic infrared sources. CMB anisotropy observations at high resolution and high frequencies are especially sensitive to this foreground. We use data from the COBE satellite to generate a Galactic emission spectrum covering mm and sub-mm wavelengths. Using this spectrum as a template, we predict the microwave emission of the 5319 brightest infrared galaxies seen by IRAS. We simulate skymaps over the relevant range of frequencies (30-900 GHz) and instrument resolutions (10'-10 degrees Full Width Half Max). Analysis of the temperature anisotropy of these skymaps shows that a reasonable observational window is available for CMB anisotropy measurements.