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.

We present the results from the application of a two-dimensional emission line detection method, EMission-line two-Dimensional (EM2D), to the near-infrared G102 grism observations obtained with the Wide-Field Camera 3 (WFC3) as part of the Cycle 22 {\em Hubble Space Telescope} Treasury Program: the Faint Infrared Grism Survey (FIGS). Using the EM2D method, we have assembled a catalog of emission line galaxies (ELGs) with resolved star formation from each of the four FIGS fields. Not only can one better assess the global properties of ELGs, but the EM2D method allows for the analysis and an improved study of the individual emission-line region {\it within} each galaxy. This paper includes a description of the methodology, advantages, and the first results of the EM2D method applied to ELGs in FIGS. The advantage of 2D emission line measurements includes significant improvement of galaxy redshift measurements, approaching the level of accuracy seen in high-spectral-resolution data, but with greater efficiency; and the ability to identify and measure the properties of multiple sites of star-formation and over scales of $\sim$ 1 kpc within individual galaxies out to z $\sim$ 4. The EM2D method also significantly improves the reliability of high-redshift ($z\sim7$) Lyman-$\alpha$ detections. Coupled with the wide field of view and high efficiency of space-based grism observations, EM2D provides a noteworthy improvement on the physical parameters that can be extracted from grism observations.

We improve the accuracy of photometric redshifts by including low-resolution spectral data from the G102 grism on the Hubble Space Telescope, which assists in redshift determination by further constraining the shape of the broadband Spectral Energy Disribution (SED) and identifying spectral features. The photometry used in the redshift fits includes near-IR photometry from FIGS+CANDELS, as well as optical data from ground-based surveys and HST ACS, and mid-IR data from Spitzer. We calculated the redshifts through the comparison of measured photometry with template galaxy models, using the EAZY photometric redshift code. For objects with F105W $< 26.5$ AB mag with a redshift range of $0 < z < 6$, we find a typical error of $\Delta z = 0.03 * (1+z)$ for the purely photometric redshifts; with the addition of FIGS spectra, these become $\Delta z = 0.02 * (1+z)$, an improvement of 50\%. Addition of grism data also reduces the outlier rate from 8\% to 7\% across all fields. With the more-accurate spectrophotometric redshifts (SPZs), we searched the FIGS fields for galaxy overdensities. We identified 24 overdensities across the 4 fields. The strongest overdensity, matching a spectroscopically identified cluster at $z=0.85$, has 28 potential member galaxies, of which 8 have previous spectroscopic confirmation, and features a corresponding X-ray signal. Another corresponding to a cluster at $z=1.84$ has 22 members, 18 of which are spectroscopically confirmed. Additionally, we find 4 overdensities that are detected at an equal or higher significance in at least one metric to the two confirmed clusters.

The Faint Infrared Grism Survey (FIGS) is a deep Hubble Space Telescope (HST) WFC3/IR (Wide Field Camera 3 Infrared) slitless spectroscopic survey of four deep fields. Two fields are located in the Great Observatories Origins Deep Survey-North (GOODS-N) area and two fields are located in the Great Observatories Origins Deep Survey-South (GOODS-S) area. One of the southern fields selected is the Hubble Ultra Deep Field. Each of these four fields were observed using the WFC3/G102 grism (0.8$\mu m$-1.15$\mu m$ continuous coverage) with a total exposure time of 40 orbits (~ 100 kilo-seconds) per field. This reaches a 3 sigma continuum depth of ~26 AB magnitudes and probes emission lines to $\approx 10^{-17}\ erg\ s^{-1} \ cm^{-2}$. This paper details the four FIGS fields and the overall observational strategy of the project. A detailed description of the Simulation Based Extraction (SBE) method used to extract and combine over 10000 spectra of over 2000 distinct sources brighter than m_F105W=26.5 mag is provided. High fidelity simulations of the observations is shown to significantly improve the background subtraction process, the spectral contamination estimates, and the final flux calibration. This allows for the combination of multiple spectra to produce a final high quality, deep, 1D-spectra for each object in the survey.

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 present the rest-1500\AA\ UV luminosity functions (LF) for star-forming galaxies during the cosmic \textit{high noon} -- the peak of cosmic star formation rate at $1.5<z<3$. We use deep NUV imaging data obtained as part of the \textit{Hubble} Ultra-Violet Ultra Deep Field (UVUDF) program, along with existing deep optical and NIR coverage on the HUDF. We select F225W, F275W and F336W dropout samples using the Lyman break technique, along with samples in the corresponding redshift ranges selected using photometric redshifts and measure the rest-frame UV LF at $z\sim1.7,2.2,3.0$ respectively, using the modified maximum likelihood estimator. We perform simulations to quantify the survey and sample incompleteness for the UVUDF samples to correct the effective volume calculations for the LF. We select galaxies down to $M_{UV}=-15.9,-16.3,-16.8$ and fit a faint-end slope of $\alpha=-1.20^{+0.10}_{-0.13}, -1.32^{+0.10}_{-0.14}, -1.39^{+0.08}_{-0.12}$ at $1.4<z<1.9$, $1.8<z<2.6$, and $2.4<z<3.6$, respectively. We compare the star formation properties of $z\sim2$ galaxies from these UV observations with results from H\alpha\ and UV$+$IR observations. We find a lack of high SFR sources in the UV LF compared to the H\alpha\ and UV$+$IR, likely due to dusty SFGs not being properly accounted for by the generic $IRX-\beta$ relation used to correct for dust. We compute a volume-averaged UV-to-H\alpha\ ratio by \textit{abundance matching} the rest-frame UV LF and H\alpha\ LF. We find an increasing UV-to-H\alpha\ ratio towards low mass galaxies ($M_\star \lesssim 5\times10^9$ M$_\odot$). We conclude that this could be due to a larger contribution from starbursting galaxies compared to the high-mass end.

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.

Galaxies at high redshifts provide a valuable tool to study cosmic dawn, and therefore it is crucial to reliably identify these galaxies. Here, we present an unambiguous and first simultaneous detection of both the Lyman-alpha emission and the Lyman break from a z = 7.512+/- 0.004 galaxy, observed in the Faint Infrared Grism Survey (FIGS). These spectra, taken with G102 grism on Hubble Space Telescope (HST), show a significant emission line detection (6 sigma) in multiple observational position angles (PA), with total integrated Ly{\alpha} line flux of 1.06+/- 0.12 e10-17erg s-1cm-2. The line flux is nearly a factor of four higher than the previous MOSFIRE spectroscopic observations of faint Ly{\alpha} emission at {\lambda} = 1.0347{\mu}m, yielding z = 7.5078+/- 0.0004. This is consistent with other recent observations implying that ground-based near-infrared spectroscopy underestimates total emission line fluxes, and if confirmed, can have strong implications for reionization studies that are based on ground-based Lyman-{\alpha} measurements. A 4-{\sigma} detection of the NV line in one PA also suggests a weak Active Galactic Nucleus (AGN), potentially making this source the highest-redshift AGN yet found. Thus, this observation from the Hubble Space Telescope clearly demonstrates the sensitivity of the FIGS survey, and the capability of grism spectroscopy to study the epoch of reionization.

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.

Current observational evidence suggests that the star formation rate (SFR) efficiency of neutral atomic hydrogen gas measured in Damped Ly-alpha Systems (DLAs) at z~3 is more than 10 times lower than predicted by the Kennicutt-Schmidt (KS) relation. To understand the origin of this deficit, and to investigate possible evolution with redshift and galaxy properties, we measure the SFR efficiency of atomic gas at z~1, z~2, and z~3 around star-forming galaxies. We use new robust photometric redshifts in the Hubble Ultra Deep Field to create galaxy stacks in these three redshift bins, and measure the SFR efficiency by combining DLA absorber statistics with the observed rest-frame UV emission in the galaxies' outskirts. We find that the SFR efficiency of HI gas at z>1 is ~1-3% of that predicted by the KS relation. Contrary to simulations and models that predict a reduced SFR efficiency with decreasing metallicity and thus with increasing redshift, we find no significant evolution in the SFR efficiency with redshift. Our analysis instead suggests that the reduced SFR efficiency is driven by the low molecular content of this atomic-dominated phase, with metallicity playing a secondary effect in regulating the conversion between atomic and molecular gas. This interpretation is supported by the similarity between the observed SFR efficiency and that observed in local atomic-dominated gas, such as in the outskirts of local spiral galaxies and local dwarf galaxies.

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.

Multiwavelength data are essential in order to provide a complete picture of galaxy evolution and to inform studies of galaxies' morphological properties across cosmic time. Here we present results of a multiwavelength investigation of the morphologies of "tadpole" galaxies at intermediate redshift (0.314<z<3.175) in the Hubble Ultra Deep Field. These galaxies were previously selected from deep Hubble Space Telescope (HST) F775W data based on their distinct asymmetric knot-plus-tail morphologies (Straughn et al. 2006). Here we use deep Wide Field Camera 3 near-infrared imaging in addition to the HST optical data in order to study the rest-frame UV/optical morphologies of these galaxies across the redshift range 0.3<z<3.2. This study reveals that the majority of these galaxies do retain their general asymmetric morphology in the rest-frame optical over this redshift range, if not the distinct "tadpole" shape. The average stellar mass of tadpole galaxies is lower than field galaxies, with the effect being slightly greater at higher redshift within the errors. Estimated from SED fits, the average age of tadpole galaxies is younger than field galaxies in the lower redshift bin, and the average metallicity is lower (whereas the specific star formation rate for tadpoles is roughly the same as field galaxies across the redshift range probed here). These average effects combined support the conclusion that this subset of galaxies is in an active phase of assembly, either late-stage merging or cold gas accretion causing localized clumpy star-formation.

During the "WISE at 5: Legacy and Prospects" conference in Pasadena, CA -- which ran from February 10 - 12, 2015 -- attendees were invited to engage in an interactive session exploring the future uses of the Wide-field Infrared Survey Explorer (WISE) data. The 65 participants -- many of whom are extensive users of the data -- brainstormed the top questions still to be answered by the mission, as well as the complementary current and future datasets and additional processing of WISE/NEOWISE data that would aid in addressing these most important scientific questions. The results were mainly bifurcated between topics related to extragalactic studies (e.g. AGN, QSOs) and substellar mass objects. In summary, participants found that complementing WISE/NEOWISE data with cross-correlated multiwavelength surveys (e.g. SDSS, Pan-STARRS, LSST, Gaia, Euclid, etc.) would be highly beneficial for all future mission goals. Moreover, developing or implementing machine-learning tools to comb through and understand cross-correlated data was often mentioned for future uses. Finally, attendees agreed that additional processing of the data such as co-adding WISE and NEOWISE and extracting a multi-epoch photometric database and parallax and proper motion catalog would greatly improve the scientific results of the most important projects identified. In that respect, a project such as MaxWISE which would execute the most important additional processing and extraction as well as make the data and catalogs easily accessible via a public portal was deemed extremely important.

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.

Narrow-band imaging of the rest-frame Lyman continuum (LyC) of galaxies at z~3.1 has produced a large number of candidate LyC-emitting galaxies. These samples are contaminated by galaxies at lower redshift. To better understand LyC escape, we need an uncontaminated sample of galaxies that emit strongly in the LyC. Here we present deep Hubble imaging of five bright galaxies at z~3.1 that had previously been identified as candidate LyC-emitters with ground-based images. The WFC3 F336W images probe the LyC of galaxies at z>3.06 and provide an order-of-magnitude increase in spatial resolution over ground-based imaging. The non-ionizing UV images often show multiple galaxies (or components) within ~1'' of the candidate LyC emission seen from the ground. In each case, only one of the components is emitting light in the F336W filter, which would indicate LyC escape if that component is at z>3.06. We use Keck/NIRSPEC near-IR spectroscopy to measure redshifts of these components to distinguish LyC-emitters from foreground contamination. We find that two candidates are low redshift contaminants, one candidate had a previously misidentified redshift, and the other two cannot be confirmed as LyC-emitters. The level of contamination is consistent with previous estimates. For the galaxies with z>3.06, we derive strong 1 sigma limits on the relative escape fraction between 0.07 and 0.09. We still do not have a sample of definitive LyC-emitters, and a much larger study of low luminosity galaxies is required. The combination of high resolution imaging and deep spectroscopy is critical for distinguishing LyC-emitters from foreground contaminants.

Recent observations have shown that star-forming galaxies like our own Milky Way evolve kinematically into ordered thin disks over the last ~8 billion years since z=1.2, undergoing a process of "disk settling." For the first time, we study the kinematic evolution of a suite of four state of the art "zoom in" hydrodynamic simulations of galaxy formation and evolution in a fully cosmological context and compare with these observations. Until now, robust measurements of the internal kinematics of simulated galaxies were lacking as the simulations suffered from low resolution, overproduction of stars, and overly massive bulges. The current generation of simulations has made great progress in overcoming these difficulties and is ready for a kinematic analysis. We show that simulated galaxies follow the same kinematic trends as real galaxies: they progressively decrease in disordered motions (sigma_g) and increase in ordered rotation (Vrot) with time. The slopes of the relations between both sigma_g and Vrot with redshift are consistent between the simulations and the observations. In addition, the morphologies of the simulated galaxies become less disturbed with time, also consistent with observations, and they both have similarly large scatter. This match between the simulated and observed trends is a significant success for the current generation of simulations, and a first step in determining the physical processes behind disk settling.

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

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

We present the multiwavelength - ultraviolet to mid-infrared - catalog of the UKIDSS Ultra-Deep Survey (UDS) field observed as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). Based on publicly available data, the catalog includes: the CANDELS data from the Hubble Space Telescope (near-infrared WFC3 F125W and F160W data and visible ACS F606W and F814W data), u-band data from CFHT/Megacam, B, V, Rc, i' and z' band data from Subaru/Suprime-Cam, Y and Ks band data from VLT/HAWK-I, J, H and K bands data from UKIDSS (Data Release 8), and Spitzer/IRAC data (3.6, 4.5 from SEDS, 5.8 and 8.0um from SpUDS). The present catalog is F160W-selected and contains 35932 sources over an area of 201.7 square arcmin and includes radio and X-ray detected sources and spectroscopic redshifts available for 210 sources.

Popular cosmological scenarios predict that galaxies form hierarchically from the merger of many progenitors, each with their own unique star formation history (SFH). We use a sophisticated approach to constrain the SFHs of 4517 blue (presumably star-forming) galaxies with spectroscopic redshifts in the range 0.2 < z < 1.4 from the All-Wavelength Extended Groth Strip International Survey (AEGIS). This consists in the Bayesian analysis of the observed galaxy spectral energy distributions with a comprehensive library of synthetic spectra assembled using realistic, hierarchical star formation and chemical enrichment histories from cosmological simulations. We constrain the SFH of each galaxy in our sample by comparing the observed fluxes in the B, R, I and Ks bands and rest-frame optical emission-line luminosities with those of one million model spectral energy distributions. We explore the dependence of the resulting SFHs on galaxy stellar mass and redshift. We find that the average SFHs of high-mass galaxies rise and fall in a roughly symmetric bell-shaped manner, while those of low-mass galaxies rise progressively in time, consistent with the typically stronger activity of star formation in low-mass compared to high-mass galaxies. For galaxies of all masses, the star formation activity rises more rapidly at high than at low redshift. These findings imply that the standard approximation of exponentially declining SFHs widely used to interpret observed galaxy spectral energy distributions may not be appropriate to constrain the physical parameters of star-forming galaxies at intermediate redshifts.

We present evidence from a sample of 544 galaxies from the DEEP2 Survey for evolution of the internal kinematics of blue galaxies with stellar masses ranging 8.0 < log M* (M_Sun) < 10.7 over 0.2<z<1.2. DEEP2 provides galaxy spectra and Hubble imaging from which we measure emission-line kinematics and galaxy inclinations, respectively. Our large sample allows us to overcome scatter intrinsic to galaxy properties in order to examine trends in kinematics. We find that at a fixed stellar mass galaxies systematically decrease in disordered motions and increase in rotation velocity and potential well depth with time. Massive galaxies are the most well-ordered at all times examined, with higher rotation velocities and less disordered motions than less massive galaxies. We quantify disordered motions with an integrated gas velocity dispersion corrected for beam smearing (sigma_g). It is unlike the typical pressure-supported velocity dispersion measured for early type galaxies and galaxy bulges. Because both seeing and the width of our spectral slits comprise a significant fraction of the galaxy sizes, sigma_g integrates over velocity gradients on large scales which can correspond to non-ordered gas kinematics. We compile measurements of galaxy kinematics from the literature over 1.2<z<3.8 and do not find any trends with redshift, likely for the most part because these datasets are biased toward the most highly star-forming systems. In summary, over the last ~8 billion years since z=1.2, blue galaxies evolve from disordered to ordered systems as they settle to become the rotation-dominated disk galaxies observed in the Universe today, with the most massive galaxies being the most evolved at any time.

We describe a mission architecture designed to substantially increase the science capability of the NASA Science Mission Directorate (SMD) Astrophysics Explorer Program for all AO proposers working within the near-UV to far-infrared spectrum. We have demonstrated that augmentation of Falcon 9 Explorer launch services with a 13 kW Solar Electric Propulsion (SEP) stage can deliver a 700 kg science observatory payload to extra-Zodiacal orbit. This new capability enables up to ~13X increased photometric sensitivity and ~160X increased observing speed relative to a Sun-Earth L2, Earth-trailing, or Earth orbit with no increase in telescope aperture. All enabling SEP stage technologies for this launch service augmentation have reached sufficient readiness (TRL-6) for Explorer Program application in conjunction with the Falcon 9. We demonstrate that enabling Astrophysics Explorers to reach extra-zodiacal orbit will allow this small payload program to rival the science performance of much larger long development time systems; thus, providing a means to realize major science objectives while increasing the SMD Astrophysics portfolio diversity and resiliency to external budget pressure. The SEP technology employed in this study has strong applicability to SMD Planetary Science community-proposed missions. SEP is a stated flight demonstration priority for NASA's Office of the Chief Technologist (OCT). This new mission architecture for astrophysics Explorers enables an attractive realization of joint goals for OCT and SMD with wide applicability across SMD science disciplines.

We describe the infrared properties of sources detected over ~36 deg^2 of sky in the GAMA 15-hr equatorial field, using data from both the Herschel Astrophysical Terahertz Large-Area Survey (H-ATLAS) and Wide-field Infrared Survey (WISE). With 5-sigma point-source depths of 34 and 0.048 mJy at 250 micron and 3.4 micron, respectively, we are able to identify 50.6% of the H-ATLAS sources in the WISE survey, corresponding to a surface density of ~630 deg^{-2}. Approximately two-thirds of these sources have measured spectroscopic or optical/near-IR photometric redshifts of z<1. For sources with spectroscopic redshifts at z<0.3, we find a linear correlation between the infrared luminosity at 3.4 micron and that at 250 micron, with +-50% scatter over ~1.5 orders of magnitude in luminosity, ~10^9 - 10^{10.5} L_sun. By contrast, the matched sources without previously measured redshifts (r>~20.5) have 250-350 micron flux density ratios that suggest either high-redshift galaxies (z>~1.5) or optically faint low-redshift galaxies with unusually low temperatures (T<~20). Their small 3.4-250 micron flux ratios favor a high-redshift galaxy population, as only the most actively star-forming galaxies at low redshift (e.g., Arp 220) exhibit comparable flux density ratios. Furthermore, we find a relatively large AGN fraction (~30%) in a 12 micron flux-limited subsample of H-ATLAS sources, also consistent with there being a significant population of high-redshift sources in the no-redshift sample.

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 use a large dark matter simulation of a LambdaCDM model to investigate the clustering and environmental dependence of the number of substructures in a halo. Focusing on redshift z=1, we find that the halo occupation distribution is sensitive at the tens of percent level to the surrounding density and to a lesser extent to asymmetry of the surrounding density distribution. We compute the autocorrelation function of halos as a function of occupation, building on the finding of Wechsler et al. (2006) and Gao and White (2007) that halos (at fixed mass) with more substructure are more clustered. We compute the relative bias as a function of occupation number at fixed mass, finding a strong relationship. At fixed mass, halos in the top 5% of occupation can have an autocorrelation function ~ 1.5-2 times higher than the mean. We also compute the bias as a function of halo mass, for fixed halo occupation. We find that for group and cluster sized halos, when the number of subhalos is held fixed, there is a strong anticorrelation between bias and halo mass. Such a relationship represents an additional challenge to the halo model.

The increasing size of cosmological simulations has led to the need for new visualization techniques. We focus on Smoothed Particle Hydrodynamical (SPH) simulations run with the GADGET code and describe methods for visually accessing the entire simulation at full resolution. The simulation snapshots are rastered and processed on supercomputers into images that are ready to be accessed through a web interface (GigaPan). This allows any scientist with a web-browser to interactively explore simulation datasets in both in spatial and temporal dimensions, datasets which in their native format can be hundreds of terabytes in size or more. We present two examples, the first a static terapixel image of the MassiveBlack simulation, a P-GADGET SPH simulation with 65 billion particles, and the second an interactively zoomable animation of a different simulation with more than one thousand frames, each a gigapixel in size. Both are available for public access through the GigaPan web interface. We also make our imaging software publicly available.

The far-ultraviolet (FUV) number counts of galaxies constrain the evolution of the star-formation rate density of the universe. We report the FUV number counts computed from FUV imaging of several fields including the Hubble Ultra Deep Field, the Hubble Deep Field North, and small areas within the GOODS-North and -South fields. These data were obtained with the Hubble Space Telescope Solar Blind Channel of the Advance Camera for Surveys. The number counts sample a FUV AB magnitude range from 21-29 and cover a total area of 15.9 arcmin^2, ~4 times larger than the most recent HST FUV study. Our FUV counts intersect bright FUV GALEX counts at 22.5 mag and they show good agreement with recent semi-analytic models based on dark matter "merger trees" by Somerville et al. (2011). We show that the number counts are ~35% lower than in previous HST studies that use smaller areas. The differences between these studies are likely the result of cosmic variance; our new data cover more lines of sight and more area than previous HST FUV studies. The integrated light from field galaxies is found to contribute between 65.9 +/-8 - 82.6 +/-12 photons/s/cm^2/sr/angstrom to the FUV extragalactic background. These measurements set a lower limit for the total FUV background light.

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

The three-point correlation function (3PCF) provides an important view into the clustering of galaxies that is not available to its lower order cousin, the two-point correlation function (2PCF). Higher order statistics, such as the 3PCF, are necessary to probe the non-Gaussian structure and shape information expected in these distributions. We measure the clustering of spectroscopic galaxies in the Main Galaxy Sample of the Sloan Digital Sky Survey (SDSS), focusing on the shape or configuration dependence of the reduced 3PCF in both redshift and projected space. This work constitutes the largest number of galaxies ever used to investigate the reduced 3PCF, using over 220,000 galaxies in three volume-limited samples. We find significant configuration dependence of the reduced 3PCF at 3-27 Mpc/h, in agreement with LCDM predictions and in disagreement with the hierarchical ansatz. Below 6 Mpc/h, the redshift space reduced 3PCF shows a smaller amplitude and weak configuration dependence in comparison with projected measurements suggesting that redshift distortions, and not galaxy bias, can make the reduced 3PCF appear consistent with the hierarchical ansatz. The reduced 3PCF shows a weaker dependence on luminosity than the 2PCF, with no significant dependence on scales above 9 Mpc/h. On scales less than 9 Mpc/h, the reduced 3PCF appears more affected by galaxy color than luminosty. We demonstrate the extreme sensitivity of the 3PCF to systematic effects such as sky completeness and binning scheme, along with the difficulty of resolving the errors. Some comparable analyses make assumptions that do not consistently account for these effects.

We constrain the linear and quadratic bias parameters from the configuration dependence of the three-point correlation function (3PCF) in both redshift and projected space, utilizing measurements of spectroscopic galaxies in the Sloan Digital Sky Survey (SDSS) Main Galaxy Sample. We show that bright galaxies (M_r < -21.5) are biased tracers of mass, measured at a significance of 4.5 sigma in redshift space and 2.5 sigma in projected space by using a thorough error analysis in the quasi-linear regime (9-27 Mpc/h). Measurements on a fainter galaxy sample are consistent with an unbiased model. We demonstrate that a linear bias model appears sufficient to explain the galaxy-mass bias of our samples, although a model using both linear and quadratic terms results in a better fit. In contrast, the bias values obtained from the linear model appear in better agreement with the data by inspection of the relative bias, and yield implied values of sigma_8 that are more consistent with current constraints. We investigate the covariance of the 3PCF, which itself is a measurement of galaxy clustering. We assess the accuracy of our error estimates by comparing results from mock galaxy catalogs to jackknife re-sampling methods. We identify significant differences in the structure of the covariance. However, the impact of these discrepancies appears to be mitigated by an eigenmode analysis that can account for the noisy, unresolved modes. Our results demonstrate that using this technique is sufficient to remove potential systematics even when using less-than-ideal methods to estimate errors.

We use a large N-body simulation to examine the detectability of HI in emission at redshift z ~ 1, and the constraints imposed by current observations on the neutral hydrogen mass function of galaxies at this epoch. We consider three different models for populating dark matter halos with HI, designed to encompass uncertainties at this redshift. These models are consistent with recent observations of the detection of HI in emission at z ~ 0.8. Whilst detection of 21 cm emission from individual halos requires extremely long integrations with existing radio interferometers, such as the Giant Meter Radio Telescope (GMRT), we show that the stacked 21 cm signal from a large number of halos can be easily detected. However, the stacking procedure requires accurate redshifts of galaxies. We show that radio observations of the field of the DEEP2 spectroscopic galaxy redshift survey should allow detection of the HI mass function at the 5-12 sigma level in the mass range 10^(11.4) M_sun/h < M_halo < 10^(12.5)M_sun/h, with a moderate amount of observation time. Assuming a larger noise level that corresponds to an upper bound for the expected noise for the GMRT, the detection significance for the HI mass function is still at the 1.7-3 sigma level. We find that optically undetected satellite galaxies enhance the HI emission profile of the parent halo, leading to broader wings as well as a higher peak signal in the stacked profile of a large number of halos. We show that it is in principle possible to discern the contribution of undetected satellites to the total HI signal, even though cosmic variance limitation make this challenging for some of our models.

We present grism spectra of emission-line galaxies (ELGs) from 0.6-1.6 microns from the Wide Field Camera 3 on the Hubble Space Telescope. These new infrared grism data augment previous optical Advanced Camera for Surveys G800L 0.6-0.95 micron grism data in GOODS-South from the PEARS program, extending the wavelength covereage well past the G800L red cutoff. The ERS grism field was observed at a depth of 2 orbits per grism, yielding spectra of hundreds of faint objects, a subset of which are presented here. ELGs are studied via the Ha, [OIII], and [OII] emission lines detected in the redshift ranges 0.2<z<1.4, 1.2<z<2.2 and 2.0<z<3.3 respectively in the G102 (0.8-1.1 microns; R~210) and G141 (1.1-1.6 microns; R~130) grisms. The higher spectral resolution afforded by the WFC3 grisms also reveals emission lines not detectable with the G800L grism (e.g., [SII] and [SIII] lines). From these relatively shallow observations, line luminosities, star-formation rates, and grism spectroscopic redshifts are determined for a total of 48 ELGs to m(AB)~25 mag. Seventeen GOODS-South galaxies that previously only had photometric redshifts now have new grism-spectroscopic redshifts, in some cases with large corrections to the photometric redshifts (Delta(z)~0.3-0.5). Additionally, one galaxy had no previously-measured redshift but now has a secure grism-spectroscopic redshift, for a total of 18 new GOODS-South spectroscopic redshifts. The faintest source in our sample has a magnitude m(AB)=26.9 mag. The ERS grism data also reflect the expected trend of lower specific star formation rates for the highest mass galaxies in the sample as a function of redshift, consistent with downsizing and discovered previously from large surveys. These results demonstrate the remarkable efficiency and capability of the WFC3 NIR grisms for measuring galaxy properties to faint magnitudes and redshifts to z>2.

We present the results of rest-frame, UV slitless spectroscopic observations of a sample of 32 z~0.7 Lyman break galaxy (LBG) analogs in the COSMOS field. The spectroscopic search was performed with the Solar Blind Channel (SBC) on Hubble Space Telescope. We report the detection of leaking Lyman continuum (LyC) radiation from an AGN-starburst composite. While we find no direct detections of LyC emission in the remainder of our sample, we achieve individual lower limits (3 sigma) of the observed non-ionizing UV to LyC flux density ratios, f_{nu}(1500A)/f_{nu}(830A) of 20 to 204 (median of 73.5) and 378.7 for the stack. Assuming an intrinsic Lyman break of 3.4 and an intergalactic medium (IGM) transmission of LyC photons along the line of sight to the galaxy of 85% we report an upper limit for the relative escape fraction in individual galaxies of 0.02 - 0.19 and a stacked 3 sigma upper limit of 0.01. We find no indication of a relative escape fraction near unity as seen in some LBGs at z~3. Our UV spectra achieve the deepest limits to date at any redshift for the escape fraction in individual sources. The contrast between these z~0.7 low escape fraction LBG analogs with z~3 LBGs suggests that either the processes conducive to high escape fractions are not being selected for in the z<1 samples or the average escape fraction is decreasing from z~3 to z~1. We discuss possible mechanisms which could affect the escape of LyC photons.

We have determined the luminosity function of 250um-selected galaxies detected in the ~14 sq.deg science demonstration region of the Herschel-ATLAS project out to a redshift of z=0.5. Our findings very clearly show that the luminosity function evolves steadily out to this redshift. By selecting a sub-group of sources within a fixed luminosity interval where incompleteness effects are minimal, we have measured a smooth increase in the comoving 250um luminosity density out to z=0.2 where it is 3.6+1.4-0.9 times higher than the local value.

We present colour-colour diagrams of detected sources in the Herschel-ATLAS Science Demonstration Field from 100 to 500 microns using both PACS and SPIRE. We fit isothermal modified black bodies to the spectral energy distribution (SED) to extract the dust temperature of sources with counterparts in Galaxy And Mass Assembly (GAMA) or SDSS surveys with either a spectroscopic or a photometric redshift. For a subsample of 330 sources detected in at least three FIR bands with a significance greater than 3 $\sigma$, we find an average dust temperature of $(28 \pm 8)$K. For sources with no known redshift, we populate the colour-colour diagram with a large number of SEDs generated with a broad range of dust temperatures and emissivity parameters, and compare to colours of observed sources to establish the redshift distribution of this sample. For another subsample of 1686 sources with fluxes above 35 mJy at 350 microns and detected at 250 and 500 microns with a significance greater than 3$\sigma$, we find an average redshift of $2.2 \pm 0.6$.

Aims. The Herschel-ATLAS survey (H-ATLAS) will be the largest area survey to be undertaken by the Herschel Space Observatory. It will cover 550 sq. deg. of extragalactic sky at wavelengths of 100, 160, 250, 350 and 500 microns when completed, reaching flux limits (5 sigma) from 32 to 145mJy. We here present galaxy number counts obtained for SPIRE observations of the first ~14 sq. deg. observed at 250, 350 and 500 microns. Methods. Number counts are a fundamental tool in constraining models of galaxy evolution. We use source catalogs extracted from the H-ATLAS maps as the basis for such an analysis. Correction factors for completeness and flux boosting are derived by applying our extraction method to model catalogs and then applied to the raw observational counts. Results. We find a steep rise in the number counts at flux levels of 100-200mJy in all three SPIRE bands, consistent with results from BLAST. The counts are compared to a range of galaxy evolution models. None of the current models is an ideal fit to the data but all ascribe the steep rise to a population of luminous, rapidly evolving dusty galaxies at moderate to high redshift.

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

Improvements in the precision of the astrophysical flux scale are needed to answer fundamental scientific questions ranging from cosmology to stellar physics. The unexpected discovery that the expansion of the universe is accelerating was based upon the measurement of astrophysical standard candles that appeared fainter than expected. To characterize the underlying physical mechanism of the "Dark Energy" responsible for this phenomenon requires an improvement in the visible-NIR flux calibration of astrophysical sources to 1% precision. These improvements will also enable large surveys of white dwarf stars, e.g. GAIA, to advance stellar astrophysics by testing and providing constraints for the mass-radius relationship of these stars. ACCESS (Absolute Color Calibration Experiment for Standard Stars) is a rocket-borne payload that will enable the transfer of absolute laboratory detector standards from NIST to a network of stellar standards with a calibration accuracy of 1% and a spectral resolving power of R = 500 across the 0.35-1.7 micron bandpass. Among the strategies being employed to minimize calibration uncertainties are: (1) judicious selection of standard stars (previous calibration heritage, minimal spectral features, robust stellar atmosphere models), (2) execution of observations above the Earth's atmosphere (eliminates atmospheric contamination of the stellar spectrum), (3) a single optical path and detector (to minimize visible to NIR cross-calibration uncertainties), (4) establishment of an a priori error budget, (5) on-board monitoring of instrument performance, and (6) fitting stellar atmosphere models to the data to search for discrepancies and confirm performance.

We have obtained deep Hubble Space Telescope far-UV images of 15 starburst galaxies at z~1.3 in the GOODS fields to search for escaping Lyman continuum photons. These are the deepest far-UV images m_{AB}=28.7, 3\sigma, 1" diameter) over this large an area (4.83 arcmin^2) and provide the best escape fraction constraints for any galaxy at any redshift. We do not detect any individual galaxies, with 3\sigma limits to the Lyman Continuum (~700 \AA) flux 50--149 times fainter (in f_nu) than the rest-frame UV (1500 \AA) continuum fluxes. Correcting for the mean IGM attenuation (factor ~2), as well as an intrinsic stellar Lyman Break (~3), these limits translate to relative escape fraction limits of f_{esc,rel}<[0.03,0.21]. The stacked limit is f_{esc,rel}(3\sigma)<0.02. We use a Monte Carlo simulation to properly account for the expected distribution of IGM opacities. When including constraints from previous surveys at z~1.3 we find that, at the 95% confidence level, no more than 8% of star--forming galaxies at z~1.3 can have relative escape fractions greater than 0.50. Alternatively, if the majority of galaxies have low, but non-zero, escaping Lyman Continuum, the escape fraction can not be more than 0.04. Both the stacked limits, and the limits from the Monte Carlo simulation suggest that the average ionizing emissivity (relative to non-ionizing UV emissivity) at z~1.3 is significantly lower than has been observed in Lyman Break Galaxies (LBGs) at z~3. If the ionizing emissivity of star-forming galaxies is in fact increasing with redshift, it would help to explain the high photoionization rates seen in the IGM at z>4 and reionization of the intergalactic medium at z>6. [Abridged]

The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 510 square degrees of the extragalactic sky, four times larger than all the other Herschel surveys combined, in five far-infrared and submillimetre bands. We describe the survey, the complementary multi-wavelength datasets that will be combined with the Herschel data, and the six major science programmes we are undertaking. Using new models based on a previous submillimetre survey of galaxies, we present predictions of the properties of the ATLAS sources in other wavebands.

We present results of a search for emission-line galaxies in the Southern Fields of the Hubble Space Telescope PEARS (Probing Evolution And Reionization Spectroscopically) grism survey. The PEARS South Fields consist of five ACS pointings (including the Hubble Ultra Deep Field) with the G800L grism for a total of 120 orbits, revealing thousands of faint object spectra in the GOODS-South region of the sky. Emission-line galaxies (ELGs) are one subset of objects that are prevalent among the grism spectra. Using a 2-dimensional detection and extraction procedure, we find 320 emission lines orginating from 226 galaxy "knots'' within 192 individual galaxies. Line identification results in 118 new grism-spectroscopic redshifts for galaxies in the GOODS-South Field. We measure emission line fluxes using standard Gaussian fitting techniques. At the resolution of the grism data, the H-beta and [OIII] doublet are blended. However, by fitting two Gaussian components to the H-beta and [OIII] features, we find that many of the PEARS ELGs have high [OIII]/H-beta ratios compared to other galaxy samples of comparable luminosities. The star-formation rates (SFRs) of the ELGs are presented, as well as a sample of distinct giant star-forming regions at z~0.1-0.5 across individual galaxies. We find that the radial distances of these HII regions in general reside near the galaxies' optical continuum half-light radii, similar to those of giant HII regions in local galaxies.

The catalog from the first high resolution U-band image of the Hubble Ultra Deep Field, taken with Hubble's Wide Field Planetary Camera 2 through the F300W filter, is presented. We detect 96 U-band objects and compare and combine this catalog with a Great Observatories Origins Deep Survey (GOODS) B-selected catalog that provides B, V, i, and z photometry, spectral types, and photometric redshifts. We have also obtained Far-Ultraviolet (FUV, 1614 \AA) data with Hubble's Advanced Camera for Surveys Solar Blind Channel (ACS/SBC) and with Galaxy Evolution Explorer (GALEX). We detected 31 sources with ACS/SBC, 28 with GALEX/FUV, and 45 with GALEX/NUV. The methods of observations, image processing, object identification, catalog preparation, and catalog matching are presented.

The rest-frame far-ultraviolet (FUV) morphologies of 8 nearby interacting and starburst galaxies (Arp 269, M 82, Mrk 8, NGC 520, NGC 1068, NGC 3079, NGC 3310, NGC 7673) are compared with 54 galaxies at z ~ 1.5 and 46 galaxies at z ~ 4 observed in the GOODS-ACS field. The nearby sample is artificially redshifted to z ~ 1.5 and 4. We compare the simulated galaxy morphologies to real z ~ 1.5 and 4 UV-bright galaxy morphologies. We calculate the Gini coefficient (G), the second-order moment of the brightest 20% of the galaxy's flux (M_20), and the Sersic index (n). We explore the use of nonparametric methods with 2D profile fitting and find the combination of M_20 with n an efficient method to classify galaxies as having merger, exponential disk, or bulge-like morphologies. When classified according to G and M_20, 20/30% of real/simulated galaxies at z ~ 1.5 and 37/12% at z ~ 4 have bulge-like morphologies. The rest have merger-like or intermediate distributions. Alternatively, when classified according to the Sersic index, 70% of the z ~ 1.5 and z ~ 4 real galaxies are exponential disks or bulge-like with n > 0.8, and ~30% of the real galaxies are classified as mergers. The artificially redshifted galaxies have n values with ~35% bulge or exponential at z ~ 1.5 and 4. Therefore, ~20-30% of Lyman-break galaxies (LBGs) have structures similar to local starburst mergers, and may be driven by similar processes. We assume merger-like or clumpy star-forming galaxies in the GOODS field have morphological structure with values n < 0.8 and M_20 > -1.7. We conclude that Mrk 8, NGC 3079, and NGC 7673 have structures similar to those of merger-like and clumpy star-forming galaxies observed at z ~ 1.5 and 4.

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

[Abridged] We detect three (plus one less certain) z-dropout sources in two separate fields of our UDF05 HST NICMOS images. These z~7 Lyman-Break Galaxy (LBG) candidates allow us to constrain the Luminosity Function (LF) of the star forming galaxy population at those epochs. By assuming a change in only M* and adopting a linear evolution in redshift, anchored to the measured values at z~6, the best fit evolution coefficient is found to be 0.43+-0.19 mag per unit redshift (0.36+-0.18, if including all four candidates), which provides a value of M*(z=7.2)=-19.7+-0.3. This implies a steady evolution for the LBG LF out to z~7, at the same rate that is observed throughout the z~3 to 6 period. This puts a strong constraint on the star-formation histories of z~6 galaxies, whose ensemble star-formation rate density must be lower by a factor 2 at ~170 Myr before the epoch at which they are observed. In particular, a large fraction of stars in the z~6 LBG population must form at redshifts well above z~7. Extrapolating this steady evolution of the LF out to higher redshifts, we estimate that galaxies would be able to reionize the universe by z~6, provided that the faint-end slope of the z>7 LF steepens to alpha~-1.9, and that faint galaxies, with luminosities below the current detection limits, contribute a substantial fraction of the required ionizing photons. This scenario gives however an integrated optical depth to electron scattering that is ~2sigma below the WMAP-5 measurement. Therefore, altogether, our results indicate that, should galaxies be the primary contributors to reionization, either the currently detected evolution of the galaxy population slows down at z>7, or the LF evolution must be compensated by a decrease in metallicity and a corresponding increase in ionization efficiency at these early epochs.

The Spitzer Adaptation of the Red-sequence Cluster Survey (SpARCS) is a z'-passband imaging survey, consisting of deep (z' ~ 24 AB) observations made from both hemispheres using the CFHT 3.6m and CTIO 4m telescopes. The survey was designed with the primary aim of detecting galaxy clusters at z >~ 1. In tandem with pre-existing 3.6um observations from the Spitzer Space Telescope SWIRE Legacy Survey, SpARCS detects clusters using an infrared adaptation of the two-filter red-sequence cluster technique. The total effective area of the SpARCS cluster survey is 41.9 deg^2. In this paper, we provide an overview of the 13.6 deg^2 Southern CTIO/MOSAICII observations. The 28.3 deg^2 Northern CFHT/MegaCam observations are summarized in a companion paper by Muzzin et al. (2008). In this paper, we also report spectroscopic confirmation of SpARCS J003550-431224, a very rich galaxy cluster at z = 1.335, discovered in the ELAIS-S1 field. To date, this is the highest spectroscopically confirmed redshift for a galaxy cluster discovered using the red-sequence technique. Based on nine confirmed members, SpARCS J003550-431224 has a preliminary velocity dispersion of 1050 +/- 230 km/s. With its proven capability for efficient cluster detection, SpARCS is a demonstration that we have entered an era of large, homogeneously-selected z > 1 cluster surveys.

The Spitzer Adaptation of the Red-sequence Cluster Survey (SpARCS) is a deep z'-band imaging survey covering the Spitzer SWIRE Legacy fields designed to create the first large homogeneously-selected sample of massive clusters at z > 1 using an infrared adaptation of the cluster red-sequence method. We present an overview of the northern component of the survey which has been observed with CFHT/MegaCam and covers 28.3 deg^2. The southern component of the survey was observed with CTIO/MOSAICII, covers 13.6 deg^2, and is summarized in a companion paper by Wilson et al. (2008). We also present spectroscopic confirmation of two rich cluster candidates at z ~ 1.2. Based on Nod-and-Shuffle spectroscopy from GMOS-N on Gemini there are 17 and 28 confirmed cluster members in SpARCS J163435+402151 and SpARCS J163852+403843 which have spectroscopic redshifts of 1.1798 and 1.1963, respectively. The clusters have velocity dispersions of 490 +/- 140 km/s and 650 +/- 160 km/s, respectively which imply masses (M200) of (1.0 +/- 0.9) x 10^{14} M_{solar} and (2.4 +/- 1.8) x 10^{14} M_{solar}. Confirmation of these candidates as bona fide massive clusters demonstrates that two-filter imaging is an effective, yet observationally efficient, method for selecting clusters at z > 1.

We present the first results of a targeted survey carried out with the 305 m Arecibo telescope to detect HI-line emission from galaxies at redshift z>0.16. The targets, selected from the Sloan Digital Sky Survey database, are non-interacting disk galaxies in relatively isolated fields. We present here the HI spectra and derived HI parameters for ten objects detected in this pilot program. All are massive disk galaxies in the redshift interval 0.17-0.25 (i.e. 2-3 Gyr look-back time), with HI masses M_HI=3-8 x 10^10 Msun and high gas mass fractions (HI - to - stellar mass ratios ~10-30%). Our results demonstrate the efficacy of exploiting Arecibo's large collecting area to measure the HI mass and rotational velocity of galaxies above redshift z=0.2. In particular, this sample includes the highest redshift detections of HI emission from individual galaxies made to date. Extension of this pilot program will allow us to study the HI properties of field galaxies at cosmological distances, thus complementing ongoing radio synthesis observations of cluster samples at z~0.2.

The Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) grism PEARS (Probing Evolution And Reionization Spectroscopically) survey provides a large dataset of low-resolution spectra from thousands of galaxies in the GOODS North and South fields. One important subset of objects in these data are emission-line galaxies (ELGs), and we have investigated several different methods aimed at systematically selecting these galaxies. Here we present a new methodology and results of a search for these ELGs in the PEARS observations of the Hubble Ultra Deep Field (HUDF) using a 2D detection method that utilizes the observation that many emission lines originate from clumpy knots within galaxies. This 2D line-finding method proves to be useful in detecting emission lines from compact knots within galaxies that might not otherwise be detected using more traditional 1D line-finding techniques. We find in total 96 emission lines in the HUDF, originating from 81 distinct "knots'' within 63 individual galaxies. We find in general that [OIII] emitters are the most common, comprising 44% of the sample, and on average have high equivalent widths (70% of [OIII] emitters having rest-frame EW>100A). There are 12 galaxies with multiple emitting knots--with different knots exhibiting varying flux values, suggesting that the differing star formation properties across a single galaxy can in general be probed at redshifts z~0.2-0.4. The most prevalent morphologies are large face-on spirals and clumpy interacting systems, many being unique detections owing to the 2D method described here, thus highlighting the strength of this technique.

Virtual observatories will give astronomers easy access to an unprecedented amount of data. Extracting scientific knowledge from these data will increasingly demand both efficient algorithms as well as the power of parallel computers. Nearly all efficient analyses of large astronomical datasets use trees as their fundamental data structure. Writing efficient tree-based techniques, a task that is time-consuming even on single-processor computers, is exceedingly cumbersome on massively parallel platforms (MPPs). Most applications that run on MPPs are simulation codes, since the expense of developing them is offset by the fact that they will be used for many years by many researchers. In contrast, data analysis codes change far more rapidly, are often unique to individual researchers, and therefore accommodate little reuse. Consequently, the economics of the current high-performance computing development paradigm for MPPs does not favor data analysis applications. We have therefore built a library, called Ntropy, that provides a flexible, extensible, and easy-to-use way of developing tree-based data analysis algorithms for both serial and parallel platforms. Our experience has shown that not only does our library save development time, it can also deliver excellent serial performance and parallel scalability. Furthermore, Ntropy makes it easy for an astronomer with little or no parallel programming experience to quickly scale their application to a distributed multiprocessor environment. By minimizing development time for efficient and scalable data analysis, we enable wide-scale knowledge discovery on massive datasets.

We present the UDF05 project, a HST Large Program of deep ACS (F606W, F775W, F850LP) and NICMOS (F110W, F160W) imaging of three fields, two of which coincide with the NICP1-4 NICMOS parallel observations of the Hubble Ultra Deep Field (HUDF). In this first paper we use the ACS data for the NICP12 field, as well as the original HUDF ACS data, to measure the UV Luminosity Function (LF) of z~5 Lyman Break Galaxies (LBGs) down to very faint levels. Specifically, based on a V-i, i-z selection criterion, we identify a sample of 101 and 133 candidate z~5 galaxies down to z_{850}=28.5 and 29.25 magnitudes in the NICP12 and in the HUDF fields, respectively. Using an extensive set of Monte Carlo simulations we derive corrections for observational biases and selection effects, and construct the rest-frame 1400 A LBG LF over the range M_{1400}=[-21.4, -17.1], i.e. down to ~0.04 L* at z~5, and complement it with data from the Subaru Deep Field (SDF) from Yoshida et al. (2006) to extend it to the brighter end (M_{1400}>-22.2). We show that: (i) Different assumptions regarding the SED distribution of the LBG population, dust properties and intergalactic absorption result in a 25% variation in the number density of LBGs at z~5; (ii) Under consistent assumptions for dust properties and intergalactic absorption, the HUDF is about 30% under-dense in z~5 LBGs relative to the NICP12 field, a variation which is well explained by cosmic variance; (iii) The faint-end slope of the LF is independent of the specific assumptions for the input physical parameters, and has a value of alpha ~ -1.6, similar to the faint-end slope of the LF that has been measured for LBGs at z~3 and z~6.

We report results of a recently-completed pre-Formulation Phase study of SPIRIT, a candidate NASA Origins Probe mission. SPIRIT is a spatial and spectral interferometer with an operating wavelength range 25 - 400 microns. SPIRIT will provide sub-arcsecond resolution images and spectra with resolution R = 3000 in a 1 arcmin field of view to accomplish three primary scientific objectives: (1) Learn how planetary systems form from protostellar disks, and how they acquire their inhomogeneous composition; (2) characterize the family of extrasolar planetary systems by imaging the structure in debris disks to understand how and where planets of different types form; and (3) learn how high-redshift galaxies formed and merged to form the present-day population of galaxies. Observations with SPIRIT will be complementary to those of the James Webb Space Telescope and the ground-based Atacama Large Millimeter Array. All three observatories could be operational contemporaneously.

We examine deep far-ultraviolet (1600 Angstrom) imaging of the Hubble Deep Field-North (HDFN) and the Hubble Ultra Deep Field (HUDF) to search for leaking Lyman continuum radiation from starburst galaxies at z~1.3. There are 21 (primarily sub-L*) galaxies with spectroscopic redshifts between 1.1<z<1.5 and none are detected in the far-UV. We fit stellar population templates to the galaxies' optical/near-infrared SEDs to determine the starburst age and level of dust attenuation, giving an accurate estimate of the intrinsic Lyman continuum ratio, f_1500/f_700, and allowing a conversion from f_700 limits to relative escape fractions. We show that previous high-redshift studies may have underestimated the amplitude of the Lyman Break, and thus the relative escape fraction, by a factor of ~2. Once the starburst age and intergalactic HI absorption are accounted for, 18 galaxies in our sample have limits to the relative escape fraction, f_esc,rel < 1.0 with some limits as low as f_esc,rel < 0.10 and a stacked limit of f_esc,rel < 0.08. This demonstrates, for the first time, that most sub-L* galaxies at high redshift do not have large escape fractions. When combined with a similar study of more luminous galaxies at the same redshift we show that, if all star-forming galaxies at z~1 have similar relative escape fractions, the value must be less than 0.14 (3 sigma). We also show that less than 20% (3 sigma) of star-forming galaxies at z~1 have relative escape fractions near unity. These limits contrast with the large escape fractions found at z~3 and suggest that the average escape fraction has decreased between z~3 and z~1. (Abridged)

(abridged) We examine the spatial and temporal stability of the HST ACS Wide Field Camera (WFC) point spread function (PSF) using the two square degree COSMOS survey. We show that stochastic aliasing of the PSF necessarily occurs during `drizzling'. This aliasing is maximal if the output pixel scale is equal to the input pixel scale of 0.05''. We show that this source of PSF variation can be significantly reduced by choosing a Gaussian drizzle kernel and by setting the output pixel size to 0.03''. We show that the PSF is temporally unstable, most likely due to thermal fluctuations in the telescope's focus. We find that the primary manifestation of this thermal drift in COSMOS images is an overall slow periodic focus change. Using a modified version of TinyTim, we create undistorted stars in a 30x30 grid across the ACS WFC CCDs. These PSF models are created for telescope focus values in the range -10 microns to +5 microns, thus spanning the allowed range of telescope focus values. We then use the approximately ten well measured stars in each COSMOS field to pick the best-fit focus value for each field. The TinyTim model stars are then used to perform PSF corrections for weak lensing allowing systematics due to incorrectly modeled PSFs to be greatly reduced. We have made the software for PSF modeling using our modified version of TinyTim available to the astronomical community. We show the effects of Charge Transfer Efficiency (CTE) degradation, which distorts the object in the readout direction, mimicking a weak lensing signal. We derive a parametric correction for the effect of CTE on the shapes of objects in the COSMOS field as a function of observation date, position within the ACS WFC field, and object flux.

We examine the origin and evolution of the mass-metallicity relationship (MZR, M-Z) for galaxies using high resolution cosmological SPH + N-Body simulations that include a physically motivated description of supernovae feedback and subsequent metal enrichment. We discriminate between two sources that may contribute to the origin of the MZR: 1) metal and baryon loss due to gas outflow, or 2) inefficient star formation at the lowest galaxy masses. Our simulated galaxies reproduce the observed MZR in shape and normalization both at z=0 and z=2. We find that baryon loss occurs due to UV heating before star formation turns on in galaxies with M_baryon < 10^8 M_sun, but that some gas loss due to supernovae induced winds is required to subsequently reproduce the low effective chemical yield observed in low mass galaxies. Despite this, we show that low star formation efficiencies, regulated by supernovae feedback, are primarily responsible for the lower metallicities of low mass galaxies and the overall M-Z trend. We find that the shape of the MZR is relatively constant with redshift, but that its normalization increases with time. Simulations with no energy feedback from supernovae overproduce metals at low galaxy masses by rapidly transforming a large fraction of their gas into stars. Despite the fact that our low mass galaxies have lost a majority of their baryons, they are still the most gas rich objects in our simulations due to their low star formation efficiencies.

We present an Ultraviolet (UV) selected sample of 268 objects in the two fields of the Great Observatories Origins Deep Survey (GOODS). We used the parallel observations taken with WFPC2 in the U--band (F300W) which covered 88% of the GOODS fields to identify sources and selected only objects with GOODS/ACS counterparts. Spectroscopic redshifts for 95 of these sources are available and we have used the multiwavelength GOODS data to estimate photometric redshifts for the others. Most of the objects are between 0.2<z<0.8. We used the spectral types obtained by the photometric redshift fitting to identify the starburst galaxies. We have also visually checked all objects and looked for tidal effects and nearby companions. We find that (i) 45% of the UV-selected galaxies are starbursts, (ii) nearly 75% of the starbursts have tidal tails or show some peculiarity typical of interactions or mergers, (iii) ~50% have companions within an area of 5x5 arcsec. The UV-selected sample has an average rest-frame M_B=--19.9 +- 0.1. The bluest objects in the sample (U-B < 0.2 and B-V < 0.1) are at 1.1<z<1.9 and have peculiar morphologies that resembles either tadpoles, chains, or double-clump galaxies. Starbursts with tadpole or clump morphology at z=0.8-1.3 have sizes comparable to LBGs and compact Ultraviolet-luminous galaxies (UVLGs).

We present far-ultraviolet (FUV) imaging of the Hubble Deep Field North (HDF-N) taken with the Solar Blind Channel of the Advanced Camera for Surveys (ACS/SBC) and the FUV MAMA detector of the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope. The full WFPC2 deep field has been observed at 1600 Angstroms. We detect 134 galaxies and one star down to a limit of FUV_{AB} ~ 29. All sources have counterparts in the WFPC2 image. Redshifts (spectroscopic or photometric) for the detected sources are in the range 0<z<1. We find that the FUV galaxy number counts are higher than those reported by GALEX, which we attribute at least in part to cosmic variance in the small HDF-N field of view. Six of the 13 Chandra sources at z<0.85 in the HDF-N are detected in the FUV, and those are consistent with starbursts rather than AGN. Cross-correlating with Spitzer sources in the field, we find that the FUV detections show general agreement with the expected L_IR/L_UV vs. Beta relationship. We infer star formation rates (SFRs), corrected for extinction using the UV slope, and find a median value of 0.3 Msun/yr for FUV-detected galaxies, with 75% of detected sources have SFR<1 Msun/yr. Examining the morphological distribution of sources, we find that about half of all FUV-detected sources are identied as spiral galaxies. Half of morphologically-selected spheroids at z<0.85 are detected in the FUV, suggesting that such sources have significant ongoing star-formation in the epoch since z=1.

The James Webb Space Telescope (JWST) is a large (6.6m), cold (50K), infrared-optimized space observatory that will be launched early in the next decade. The observatory will have four instruments: a near-infrared camera, a near-infrared multi-object spectrograph, and a tunable filter imager will cover the wavelength range, 0.6 to 5.0 microns, while the mid-infrared instrument will do both imaging and spectroscopy from 5.0 to 29 microns. The JWST science goals are divided into four themes. The End of the Dark Ages: First Light and Reionization theme seeks to identify the first luminous sources to form and to determine the ionization history of the early universe. The Assembly of Galaxies theme seeks to determine how galaxies and the dark matter, gas, stars, metals, morphological structures, and active nuclei within them evolved from the epoch of reionization to the present day. The Birth of Stars and Protoplanetary Systems theme seeks to unravel the birth and early evolution of stars, from infall on to dust-enshrouded protostars to the genesis of planetary systems. The Planetary Systems and the Origins of Life theme seeks to determine the physical and chemical properties of planetary systems including our own, and investigate the potential for the origins of life in those systems. To enable these observations, JWST consists of a telescope, an instrument package, a spacecraft and a sunshield. The telescope consists of 18 beryllium segments, some of which are deployed. The segments will be brought into optical alignment on-orbit through a process of periodic wavefront sensing and control. The JWST operations plan is based on that used for previous space observatories, and the majority of JWST observing time will be allocated to the international astronomical community through annual peer-reviewed proposal opportunities.

We present measurements of the normalised redshift-space three-point correlation function (Q_z) of galaxies from the Sloan Digital Sky Survey (SDSS) main galaxy sample. We have applied our "npt" algorithm to both a volume-limited (36738 galaxies) and magnitude-limited sample (134741 galaxies) of SDSS galaxies, and find consistent results between the two samples, thus confirming the weak luminosity dependence of Q_z recently seen by other authors. We compare our results to other Q_z measurements in the literature and find it to be consistent within the full jack-knife error estimates. However, we find these errors are significantly increased by the presence of the ``Sloan Great Wall'' (at z ~ 0.08) within these two SDSS datasets, which changes the 3-point correlation function (3PCF) by 70% on large scales (s>=10h^-1 Mpc). If we exclude this supercluster, our observed Q_z is in better agreement with that obtained from the 2dFGRS by other authors, thus demonstrating the sensitivity of these higher-order correlation functions to large-scale structures in the Universe. This analysis highlights that the SDSS datasets used here are not ``fair samples'' of the Universe for the estimation of higher-order clustering statistics and larger volumes are required. We study the shape-dependence of Q_z(s,q,theta) as one expects this measurement to depend on scale if the large scale structure in the Universe has grown via gravitational instability from Gaussian initial conditions. On small scales (s <= 6h^-1 Mpc), we see some evidence for shape-dependence in Q_z, but at present our measurements are consistent with a constant within the errors (Q_z ~ 0.75 +/- 0.05). On scales >10h^-1 Mpc, we see considerable shape-dependence in Q_z.

We present the analysis of the deepest near-UV image obtained with HST using the WFPC2(F300W) as part of the parallel observations of the Ultra Deep Field campaign. The U-band 10sigma limiting magnitude measured over 0.2 arcsec square is m(AB)=27.5 which is 0.5 magnitudes deeper than that in the HDF-North. We matched the U-band catalog with those in the ACS images taken during the GOODS observations of the CDF-South and obtained photometric-z for 306 matched objects. We find that the UV-selected galaxies span all the major morphological types at 0.2<z_phot<1.2. However, disks are more common at lower redshifts, 0.2<z_phot<0.8. Higher-z objects (0.7<z_phot<1.2) are on average bluer than lower-z and have spectral type typical of starbursts. Their morphologies are compact, peculiar or low surface brightness galaxies. The average half-light radius (rest-frame 1200--1800 A) of the UV-selected galaxies at 0.66<z_ phot<1.5 is 0.26 +- 0.01 arcsec (2.07 +- 0.08 kpc). The UV-selected galaxies are on average fainter (M_B=-18.43+-0.13) than Lyman Break Galaxies (M_B=-23+-1). Our sample includes early-type galaxies that are presumably massive and forming stars only in their cores, as well as starburst-type systems that are more similar to the LBGs, although much less luminous. This implies that even the starbursts in our sample are either much less massive than LBGs or are forming stars at a much lower rate or both. The low surface brightness galaxies have no overlap with the LBGs and form an interesting new class of their own.

We use cosmological LCDM numerical simulations to model the evolution of the substructure population in sixteen dark matter haloes with resolutions of up to seven million particles within the virial radius. The combined substructure circular velocity distribution function (VDF) for hosts of 10^11 to 10^14 Msun at redshifts from zero to two or higher has a self-similar shape, is independent of host halo mass and redshift, and follows the relation: dn/dv=(1/8)(v_cmax/v_cmax,host)^-4. Halo to halo variance in the VDF is a factor of roughly two to four. At high redshifts, we find preliminary evidence for fewer large substructure haloes (subhaloes). Specific angular momenta are significantly lower for subhaloes nearer the host halo centre where tidal stripping is more effective. The radial distribution of subhaloes is marginally consistent with the mass profile for r >~ 0.3r_vir, where the possibility of artificial numerical disruption of subhaloes can be most reliably excluded by our convergence study, although a subhalo distribution that is shallower than the mass profile is favoured. Subhalo masses but not circular velocities decrease toward the host centre. Subhalo velocity dispersions hint at a positive velocity bias at small radii. There is a weak bias toward more circular orbits at lower redshift, especially at small radii. We additionally model a cluster in several power law cosmologies of P ~ k^n, and demonstrate that a steeper spectral index, n, results in significantly less substructure.

[Abridged] We test the reliability of a method to measure the mean halo mass of Damped Ly-alpha absorbers (DLAs). The method is based on measuring the ratio of the cross-correlation between DLAs and galaxies to the auto-correlation of the galaxies themselves ($w_{\rm dg}/w_{\rm gg}$), which is (in linear theory) the ratio of their bias factor. This is shown to be true irrespective of the galaxy redshift distribution, provided that one uses the same galaxies for the two correlation functions. The method is applicable to all redshifts. Here, we focus on z=3 DLAs and we demonstrate that the method robustly constrains the mean DLA halo mass using smoothed particle hydrodynamics (SPH) cosmological simulations. If we use the bias formalism of Mo & White with the DLA and galaxy mass distributions of these simulations, we predict a bias ratio of 0.771. Direct measurement from the simulations of $w_{\rm dg}/w_{\rm gg}$ st yields a ratio of 0.73+/-0.08, in excellent agreement with that prediction. Equivalently, inverting the measured correlation ratio to infer a mean DLA halo mass yields (log. averaging, in solar units) <log(M_DLA)> =11.13+/-013, in excellent agreement with the true value in the simulations: 11.16. The cross- correlation method thus appears to yield a robust estimate of the average host halo mass even though the DLAs and the galaxies occupy a broad mass spectrum of halos, and massive halos contain multiple galaxies with DLAs. We show that the inferred mean DLA halo mass is independent of the galaxy sub-sample used, i.e. the cross-correlation technique is also reliable. Our results imply that the cross-correlation length between DLAs and LBGs is predicted to be, at most, 2.85 Mpc. Future observations will soon distinguish models in which DLAs are in low mass halos from those in which DLAs are in massive halos.

We use numerical simulations in a Lambda CDM cosmology to model density profiles in a set of 16 dark matter haloes with resolutions of up to 7 million particles within the virial radius. These simulations allow us to follow robustly the formation and evolution of the central cusp over a large mass range of 10^11 to 10^14 M_sun, down to approximately 0.5% of the virial radius, and from redshift 5 to the present. The cusp of the density profile is set at redshifts of 2 or greater and remains remarkably stable to the present time, when considered in non-comoving coordinates. We fit our haloes to a 2 parameter profile where the steepness of the asymptotic cusp is given by gamma, and its radial extent is described by the concentration, c_gamma. In our simulations, we find gamma = 1.4 - 0.08Log(M/M_*) for haloes of 0.01M_* to 1000M_*, with a large scatter of gamma ~ +/-0.3$; and c_gamma = 8*M/M_*^{-0.15}, with a large M/M_* dependent scatter roughly equal to +/- c_gamma. Our redshift zero haloes have inner slope parameters ranging approximately from r^{-1} to r^{-1.5}, with a median of roughly r^{-1.3}. This 2 parameter profile fit works well for all our halo types, whether or not they show evidence of a steep asymptotic cusp. We also model a cluster in power law cosmologies of P ~ k^n (n=0,-1,-2,-2.7). We find larger concentration radii and shallower cusps for steeper n. The minimum resolved radius is well described by the mean interparticle separation. The trend of steeper and more concentrated cusps for smaller $M/M_*$ haloes clearly shows that dwarf sized Lambda CDM haloes have, on average, significantly steeper density profiles within the inner few percent of the virial radius than inferred from recent observations. Code to reproduce this profile can be downloaded from http://www.icc.dur.ac.uk/~reed/profile.html

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

We present the results of a 200 ks Chandra observation of part of the Groth Strip region, using the ACIS-I instrument. We present a relatively simple method for the detection of point-sources and calculation of limiting sensitivities, which we argue is at least as sensitive and more self-consistent than previous methods presented in the literature. 158 distinct X-ray sources are included in our point-source catalogue in the ACIS-I area. The number counts show a relative dearth of X-ray sources in this region. For example at a flux limit of 1E-15 (cgs) around 20 per cent more soft band sources are detected in the HDF-N and almost 50 per cent more in the ELAIS-N1 field, which we have analysed by the same method for comparison. We find, however, that these differences are consistent with Poisson variations at 2 sigma significance, and therefore there is no evidence for cosmic variance based on these number counts alone. We determine the average spectra of the objects and find a marked difference between the soft-band selected sources, which have Gamma=1.9 typical of unobscured AGN, and the hard-band selected sources, which have Gamma=1.0. Reassuringly, the sample as a whole has a mean spectrum of Gamma=1.4\pm 0.1, the same as the X-ray background. Nonetheless, our results imply that the fraction of sources with significant obscuration is only ~25 per cent, much less than predicted by standard AGN population synthesis models. This is confirmed by direct spectral fitting, with only a handful of objects showing evidence for absorption. After accounting for absorption, all objects are consistent with mean intrinsic spectrum of Gamma=1.76 \pm 0.08, very similar to local Seyferts (abridged).

We present measurements of the redshift--space three-point correlation function of galaxies in the Sloan Digital Sky Survey (SDSS). For the first time, we analyze the dependence of this statistic on galaxy morphology, color and luminosity. In order to control systematics due to selection effects, we used $r$--band, volume-limited samples of galaxies, constructed from the magnitude-limited SDSS data ($14.5<r<17.5$), and further divided the samples into two morphological types (early and late) or two color populations (red and blue). The three-point correlation function of SDSS galaxies follow the hierarchical relation well and the reduced three-point amplitudes in redshift--space are almost scale-independent ($Q_z=0.5\sim1.0$). In addition, their dependence on the morphology, color and luminosity is not statistically significant. Given the robust morphological, color and luminosity dependences of the two-point correlation function, this implies that galaxy biasing is complex on weakly non-linear to non-linear scales. We show that simple deterministic linear relation with the underlying mass could not explain our measurements on these scales.

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

Following the successful mapping of the last scattering surface by WMAP and balloon experiments, the epoch of the first stars, when Population III stars formed, is emerging as the next cosmological frontier. It is not clear what these stars' properties were, when they formed or how long their era lasted before leading to the stars and galaxies we see today. We show that these questions can be answered with the current and future measurements of the near-IR cosmic infrared background (CIB). Theoretical arguments suggest that Population III stars were very massive and short-lived stars that formed at $z\sim 10-20$ at rare peaks of the density field in the cold-dark-matter Universe. Because Population III stars probably formed individually in small mini-halos, they are not directly accessible to current telescopic studies. We show that these stars left a strong and measurable signature via their contribution to the CIB anisotropies for a wide range of their formation scenarios. The excess in the recently measured near-IR CIB anisotropies over that from normal galaxies can be explained by contribution from early Population III stars. These results imply that Population III were indeed very massive stars and their epoch started at $z\sim 20$ and lasted past $z\lsim 13$. We show the importance of accurately measuring the CIB anisotropies produced by Population III with future space-based missions.

We simulate the formation of a realistic disk galaxy within the hierarchical scenario of structure formation and study its internal properties to the present epoch. We compare results from a LambdaCDM simulation with a LambdaWDM (2keV) simulation that forms significantly less small scale structure. We show how high mass and force resolution in both the gas and dark matter components play an important role in solving the angular momentum catastrophe claimed from previous simulations of galaxy formation within the hierarchical framework. The stellar material in the disk component has a final specific angular momentum equal to 40% and 90% of that of the dark halo in the LambdaCDM and LambdaWDM models respectively. The LambdaWDM galaxy has a drastically reduced satellite population and a negligible stellar spheroidal component. Encounters with satellites play only a minor role in disturbing the disk. Satellites possess a variety of star formation histories linked to mergers and pericentric passages along their orbit around the primary galaxy. In both cosmologies, the galactic halo retains most of the baryons accreted and builds up a hot gas phase with a substantial X-ray emission. Therefore, while we have been successful in creating a realistic stellar disk in a massive galaxy within the LambdaCDM scenario, energy injection emerges as necessary ingredient to reduce the baryon fraction in galactic halos, independent of the cosmology adopted. (abridged)

We examine the evolution of the sizes and number densities of disk galaxies using the high resolution images obtained by the Great Observatories Origins Deep Survey (GOODS) with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST). The multiwavelength (B,V,i,z) images allow us to classify galaxies based on their rest-frame B-band morphologies out to redshift z = 1.25. In order to minimize the effect of selection biases, we confine our analysis to galaxies which occupy the region of magnitude-size plane where the survey is almost 90% complete at all redshifts. The observed size distribution is consistent with a log-normal distribution as seen for the disk galaxies in the local Universe and does not show any significant evolution over the redshift range 0.25 < z < 1.25. We find that the number densities of disk galaxies remains fairly constant over this redshift range, although a modest evolution by a factor of four may be possible within the 2-sigma uncertainties.

Weak lensing by large-scale structure allows a direct measure of the dark matter distribution. We have used parallel images taken with the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope to measure weak lensing, or cosmic shear. We measure the shapes of 26036 galaxies in 1292 STIS fields and measure the shear variance at a scale of 0.51 arcminutes. The charge transfer efficiency (CTE) of STIS has degraded over time and introduces a spurious ellipticity into galaxy shapes during the readout process. We correct for this effect as a function of signal to noise and CCD position. We further show that the detected cosmic shear signal is nearly constant in time over the approximately four years of observation. We detect cosmic shear at the 5.1 sigma level, and our measurement of the shear variance is consistent with theoretical predictions in a LambdaCDM universe. This provides a measure of the normalization of the mass power spectrum sigma_8=(1.02 +- 0.16) (0.3/Omega_m)^{0.46} (0.21/Gamma)^{0.18}$. The one-sigma error includes noise, cosmic variance, systematics and the redshift uncertainty of the source galaxies. This is consistent with previous cosmic shear measurements, but tends to favor those with a high value of sigma_8. It is also consistent with the recent determination of sigma_8 from the Wilkinson Microwave Anisotropy Probe (WMAP) experiment.

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

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

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

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

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

Understanding the physical formation of the Hubble sequence remains one of the most important unsolved astrophysical problems. Searches for proto-disks and proto-ellipticals can now be effectively done using deep wide-field Hubble Space Telescope images taken with the new Advanced Camera for Surveys. Through an analysis of concentrations (C), asymmetries (A) and clumpiness values (S) (CAS) of galaxies found in the GOODS Field South, we are able to identify objects possibly forming onto the Hubble sequence. Using this approach, we detect a sizeable population of star forming luminous diffuse objects and star forming luminous asymmetric objects between redshifts 0.5 < z < 2. These galaxies have extremely low light concentrations, or high asymmetries, with absolute magnitudes M_B < -19. The luminous diffuse objects are found in abundance between z = 1 - 2, with fewer objects at z > 2 and z < 1. The luminous asymmetric objects are found at a similar abundance, with a peak at z~1. We argue that these galaxies are a subset of modern disks and ellipticals in formation. The co-moving volume density of the luminous diffuse objects between z = 1 - 2 is similar to the local density of bright disk galaxies, with values ~5*10^{5} Gpc^{-3}. The SEDs of these objects are mostly consistent with starbursts, or star-forming normal galaxies, with average uncorrected for extinction star formation rates of ~4 M_0 yr^{-1}. These galaxies also host 35-40% of the star formation activity at 1 < z < 2. We briefly discuss the implications of these objects for understanding the origin of the Hubble sequence.

We present the first results of the pure parallel near-UV observations with WFPC2 in the Chandra Deep Field South. The main goal is to search for the population of objects that contribute to the rise in the volume-averaged star formation rate at intermediate redshifts (z<1). A total of 34 objects, were identified in the WFPC2 images and their counterparts found in the ACS GOODS images. We use template fitting of the spectral energy distributions to classify the objects and found that 17 are starbursts, 9 are late-type, 2 are early-type galaxies, and 6 are stars; two of the late-types and one of the starbursts are X-ray sources. The colors of the starbursts were reproduced by the stellar evolutionary synthesis code Starburst99 which shows a mixed population of bursts with ages < 1 Gyr. Analysis of the light concentration, asymmetry and clumpiness, shows that this sample is a mixed bag, containing dwarf ellipticals, early and late-spirals and peculiar objects which resemble mergers in progress.

We use a high resolution $\Lambda$CDM numerical simulation to calculate the mass function of dark matter haloes down to the scale of dwarf galaxies, back to a redshift of fifteen, in a 50 $h^{-1}$Mpc volume containing 80 million particles. Our low redshift results allow us to probe low $\sigma$ density fluctuations significantly beyond the range of previous cosmological simulations. The Sheth and Tormen mass function provides an excellent match to all of our data except for redshifts of ten and higher, where it overpredicts halo numbers increasingly with redshift, reaching roughly 50 percent for the $10^{10}-10^{11} \msun$ haloes sampled at redshift 15. Our results confirm previous findings that the simulated halo mass function can be described solely by the variance of the mass distribution, and thus has no explicit redshift dependence. We provide an empirical fit to our data that corrects for the overprediction of extremely rare objects by the Sheth and Tormen mass function. This overprediction has implications for studies that use the number densities of similarly rare objects as cosmological probes. For example, the number density of high redshift (z $\simeq$ 6) QSOs, which are thought to be hosted by haloes at 5$\sigma$ peaks in the fluctuation field, are likely to be overpredicted by at least a factor of 50%. We test the sensitivity of our results to force accuracy, starting redshift, and halo finding algorithm.

We present the luminosity function of [OII]-emitting galaxies at a median redshift of z=0.9, as measured in the deep spectroscopic data in the STIS Parallel Survey (SPS). The luminosity function shows strong evolution from the local value, as expected. By using random lines of sight, the SPS measurement complements previous deep single field studies. We calculate the density of inferred star formation at this redshift by converting from [OII] to H-alpha line flux as a function of absolute magnitude and find rho_dot=0.043 +/- 0.014 Msun/yr/Mpc^3 at a median redshift z~0.9 within the range 0.46<z<1.415 (H_0 = 70 km/s/Mpc, Omega_M=0.3, Omega_Lambda=0.7. This density is consistent with a (1+z)^4 evolution in global star formation since z~1. To reconcile the density with similar measurements made by surveys targeting H-alpha may require substantial extinction correction.

In the first three years of operation STIS obtained slitless spectra of approximately 2500 fields in parallel to prime HST observations as part of the STIS Parallel Survey (SPS). The archive contains almost 300 fields at high galactic latitude (|b|>30) with spectroscopic exposure times greater than 3000 seconds. This sample contains 220 fields (excluding special regions and requiring a consistent grating angle) observed between 6 June 1997 and 21 September 2000, with a total survey area of about 160 square arcminutes. At this depth, the SPS detects an average of one emission line galaxy per three fields. We present the analysis of these data, and the identification of 131 low to intermediate redshift galaxies detected by optical emission lines. The sample contains 78 objects with emission lines that we infer to be redshifted [OII]3727 emission at 0.43<z<1.7. The comoving number density of these objects is comparable to that of H-alpha emitting galaxies in the NICMOS parallel observations. One quasar and three probable Seyfert galaxies are detected. Many of the emission-line objects show morphologies suggestive of mergers or interactions. The reduced data are available upon request from the authors.

As part of the Hubble Deep Field South program, a set of shorter 2-orbit observations were obtained of the area adjacent to the deep fields. The WFPC2 flanking fields cover a contiguous solid angle of 48 square arcminutes. Parallel observations with the STIS and NICMOS instruments produce a patchwork of additional fields with optical and near-infrared (1.6 micron) response. Deeper parallel exposures with WFPC2 and NICMOS were obtained when STIS observed the NICMOS deep field. These deeper fields are offset from the rest, and an extended low surface brightness object is visible in the deeper WFPC2 flanking field. In this data paper, which serves as an archival record of the project, we discuss the observations and data reduction, and present SExtractor source catalogs and number counts derived from the data. Number counts are broadly consistent with previous surveys from both ground and space. Among other things, these flanking field observations are useful for defining slit masks for spectroscopic follow-up over a wider area around the deep fields, for studying large-scale structure that extends beyond the deep fields, for future supernova searches, and for number counts and morphological studies, but their ultimate utility will be defined by the astronomical community.

We describe the X-ray properties of a large sample of $z\sim3$ Lyman Break Galaxies (LBGs) in the region of the Hubble Deep Field North, derived from the 1 Ms public Chandra observation. Of our sample of 148 LBGs, four are detected individually. This immediately gives a measure of the bright AGN fraction in these galaxies of $\sim 3$~per cent, which is in agreement with that derived from the UV spectra. The X-ray color of the detected sources indicates that they are probably moderately obscured. Stacking of the remainder shows a significant detection ($6\sigma$) with an average luminosity of $3.5 \times 10^{41}$~erg s$^{-1}$ per galaxy in the rest frame 2-10 keV band. We have also studied a comparison sample of 95 z$\sim 1$ ``Balmer Break'' galaxies. Eight of these are detected directly, with at least two clear AGN based on their high X-ray luminosity and very hard X-ray spectra respectively. The remainder are of relatively low luminosity ($<10^{42}$~erg s$^{-1}$), and the X-rays could arise from either AGN or rapid star-formation. The X-ray colors and evidence from other wavebands favor the latter interpretation. Excluding the clear AGN, we deduce a mean X-ray luminosity of $6.6 \times 10^{40}$~erg s$^{-1}$, a factor $\sim 5$ lower than the LBGs. The average ratio of the UV and X-ray luminosities of these starforming galaxies $L_{\rm UV}/L_{\rm X}$, however, is approximately the same at $z = 1$ as it is at $z = 3$. This scaling implies that the X-ray emission follows the current star formation rate, as measured by the UV luminosity. We use our results to constrain the star formation rate at $z\sim 3$ from an X-ray perspective (truncated).

The question "How did we get here and what will the future bring?" captures the human imagination and the attention of the National Academy of Science's Astronomy and Astrophysics Survey Commitee (AASC). Fulfillment of this "fundamental goal" requires astronomers to have sensitive, high angular and spectral resolution observations in the far-infrared/submillimeter (far-IR/sub-mm) spectral region. With half the luminosity of the universe and vital information about galaxy, star and planet formation, observations in this spectral region require capabilities similar to those currently available or planned at shorter wavelengths. In this paper we summarize the scientific motivation, some mission concepts and technology requirements for far-IR/sub-mm space interferometers that can be developed in the 2010-2020 timeframe.

This paper discusses a semi-numerical method of investigating the enrichment of the intergalactic medium by galactic winds. We find that most galaxies at z >~ 3 should be driving winds, and that (if these winds are similar to those at low-z) these winds should escape to large distances. Our calculations -- which permit exploration of a large region of model parameter space -- indicate that the wind velocity, the mass of the wind-driving galaxies, the fraction of ambient material entrained, and the available time (between wind launch and the observed redshift) all affect wind propagation significantly; other physical effects can be important but are sub-dominant. We find that under reasonable assumptions, the enrichment by 3 <~ z <~ 6 galaxies could account for the quantity of metals seen in the Ly-alpha forest, though it is presently unclear whether this enrichment is compatible with the intergalactic medium's detailed metal distribution or relative quiescence.

Observations have established that the diffuse intergalactic medium (IGM) at z ~ 3 is enriched to ~0.1-1% solar metallicity and that the hot gas in large clusters of galaxies (ICM) is enriched to 1/3-1/2 solar metallicity at z=0. Metals in the IGM may have been removed from galaxies (in which they presumably form) during dynamical encounters between galaxies, by ram-pressure stripping, by supernova-driven winds, or as radiation-pressure driven dust efflux. This study develops a method of investigating the chemical enrichment of the IGM and of galaxies, using already completed cosmological simulations. To these simulations, we add dust and (gaseous) metals, distributing the dust and metals in the gas according to three simple parameterized prescriptions, one for each enrichment mechanism. These prescriptions are formulated to capture the basic ejection physics, and calibrated when possible with empirical data. Our results indicate that dynamical removal of metals from >~ 3*10^8 solar mass galaxies cannot account for the observed metallicity of low-column density Ly-alpha absorbers, and that dynamical removal from >~ 3*10^10 solar mass galaxies cannot account for the ICM metallicities. Dynamical removal also fails to produce a strong enough mass-metallicity relation in galaxies. In contrast, either wind or radiation-pressure ejection of metals from relatively large galaxies can plausibly account for all three sets of observations (though it is unclear whether metals can be distributed uniformly enough in the low-density regions without overly disturbing the IGM, and whether clusters can be enriched quite as much as observed). We investigate in detail how our results change with variations in our assumed parameters, and how results for the different ejection processes compare. (Abridged)

We analyze high-redshift structure in three hydrodynamic simulations that have identical initial conditions and cosmological parameters and differ only in the value of the baryon density parameter, Omega_b=0.02, 0.05, 0.125. Increasing Omega_b does not change the fraction of baryons in the diffuse (unshocked) phase of the intergalactic medium, but it increases cooling rates and therefore transfers some baryons from the shocked intergalactic phase to the condensed phase associated with galaxies. Predictions of Lyman-alpha forest absorption are almost unaffected by changes of Omega_b provided that the UV background intensity is adjusted so that the mean opacity of the forest matches the observed value. The required UV background intensity scales as Omega_b^1.7, and the higher photoionization rate increases the gas temperature in low density regions. Damped Lyman-alpha absorption and Lyman limit absorption both increase with increasing Omega_b, though the impact is stronger for damped absorption and is weaker at z=4 than at z=2-3. The mass of cold gas and stars in high-redshift galaxies increases faster than Omega_b but slower than Omega_b^2, and the global star formation rate scales approximately as Omega_b^1.5. In the higher Omega_b models, the fraction of baryonic material within the virial radius of dark matter halos is usually higher than the universal fraction, indicating that gas dynamics and cooling can lead to over-representation of baryons in virialized systems. On the whole, our results imply a fairly intuitive picture of the influence of Omega_b on high-redshift structure, and we provide scalings that can be used to estimate the impact of Omega_b uncertainties on the predictions of hydrodynamic simulations.

Studies of quasar absorption lines reveal that the low density intergalactic medium at z ~ 3 is enriched to 0.001-0.01 solar metallicity. This enrichment may have occurred in an early generation of Population III stars at redshift z >~ 10, by protogalaxies at 6 <~ z <~ 10, or by larger galaxies at 3 <~ z <~ 6. This paper addresses the third possibility by calculating the enrichment of the IGM at z >~ 3 by galaxies of baryonic mass >~ 10^8.5 solar masses. We use already completed cosmological simulations to which we add a prescription for chemical evolution and metal ejection by winds, assuming that the winds have properties similar to those observed in local starbursts and Lyman-break galaxies. Results are given for a number of representative models, and we also examine the properties of the galaxies responsible for the enrichment as well as the physical effects responsible for wind escape and propagation. We find that winds of velocity >~ 200-300 km/s are capable of enriching the IGM to the mean level observed, though many low-density regions would remain metal free. Calibrated by observations of Lyman-break galaxies, our calculations suggest that most galaxies at z >~ 3 should drive winds that can escape and propagate to large radii. The primary effect limiting the enrichment of low-density IG gas in our scenario is then the travel time from high- to low-density regions, implying that the metallicity of low-density gas is a strong function of redshift.

The presence of metals in hot cluster gas and in Ly-alpha absorbers, as well as the mass-metallicity relation of observed galaxies, suggest that galaxies lose a significant fraction of their metals to the intergalactic medium (IGM). Theoretical studies of this process have concentrated on metal removal by dynamical processes or supernova-driven winds. Here, we investigate the enrichment of the IGM by the expulsion of dust grains from galaxies by radiation pressure. We use already completed cosmological simulations, to which we add dust assuming that most dust can reach the equilibrium point between radiation pressure and gravitational forces. We find that the expulsion of dust and its subsequent (partial) destruction in the IGM can plausibly account for the observed level of C and Si enrichment of the z=3 IGM. At low-z, dust ejection and destruction could explain a substantial fraction of the metals in clusters, but it cannot account for all of the chemical species observed. Dust expelled by radiation pressure could give clusters a visual opacity of up to 0.2-0.5 mag in their central regions even after destruction by the hot intracluster medium; this value is interestingly close to limits and claimed observations of cluster extinction. We also comment on the implications of our results for the opacity of the general IGM. Finally, we suggest a possible `hybrid' scenario in which winds expel gas and dust into galaxy halos but radiation pressure distributes the dust uniformly through the IGM.

I present results from numerical N-body simulations regarding the effect of merging events on the angular momentum distribution of galactic halos as well as a comparison of halo growth in Press-Schechter vs. N-body methods. A total of six simulations are used spanning 3 cosmologies: a standard flat Omega_0=1 model, an open Omega_0=0.3 model and a "tilted" flat Omega_0=1 model with spectral index n=0.8. In each model, one run was conducted using a spatially uniform grid of particles and one using a refined grid in a large void. In all three models and all environments tested, the mean angular momentum of merger remnants (halo interaction products with mass ratios 3:1 or less) is greater than non-merger remnants. Furthermore, the dispersion in the merger-remnant angular momentum distribution is smaller than the dispersion of the non-merger distribution. The interpretation most consistent with the data is that the orbital angular momentum of the interactors is important in establishing the final angular momentum of the merger product. I give the angular momentum distribution which describes the merger remnant population. I trace the most massive progenitor of L_* galactic-mass halos (uniform grid) and 10^{11} solar mass halos (refined void) from redshift z=0 back to z=5. Monte-Carlo mass histories match simulations reasonably well for the latter sample. I find that for halos of mass 10^{12} to 10^{14} solar masses, this method can underestimate the mass of progenitors by 20%, hence yielding improper formation redshifts of halos. With this caveat, however, the general shapes of halo mass histories and formation-time distributions are preserved.

We use hydrodynamic cosmological simulations to study damped Lyman-alpha (DLA) and Lyman limit (LL) absorption at redshifts z=2-4 in five variants of the cold dark matter scenario. Our standard simulations resolve the formation of dense concentrations of neutral gas in halos with circular velocity v_c roughly 140 km/s for Omega_m=1 and 90 km/s for Omega_m=0.4, at z=2; an additional LCDM simulation resolves halos down to v_c approximately 50 km/s at z=3. We find a clear relation between HI column density and projected distance to the center of the nearest galaxy, with DLA absorption usually confined to galactocentric radii less than 10-15 kpc and LL absorption arising out to projected separations of 30 kpc or more. Detailed examination provides evidence of non-equilibrium effects on absorption cross-section. If we consider only absorption in the halos resolved by our standard simulations, then all five models fall short of reproducing the observed abundance of DLA and LL systems at these redshifts. If we extrapolate to lower halo masses, we find all four models are consistent with the observed abundance of DLA systems if the the extrapolated behavior extends to circular velocities roughly 50-80 km/s, and they may produce too much absorption if the relation continues to 40 km/s. Our results suggest that LL absorption is closely akin to DLA absorption, arising in less massive halos or at larger galactocentric radii but not caused by processes acting on a radically different mass scale.

We present the galaxy-galaxy angular correlations as a function of photometric redshift in a deep-wide galaxy survey centered on the Hubble Deep Field South. Images were obtained with the Big Throughput Camera on the Blanco 4m telescope at CTIO, of 1/2 square degree in broad-band uBVRI, reaching ~24th mag. Approximately 40,000 galaxies are detected in the survey. We determine photometric redshifts using galaxy template fitting to the photometry. Monte Carlo simulations show that redshifts from these data should be reliable out to z~1, where the 4000 Angstrom break shifts into the I-band. The inferred redshift distribution, n(z), shows good agreement with the distribution of galaxies measured in the HDF North and the Canada-France Redshift Survey. After assigning galaxies to redshift bins with width Delta_z=0.33, we determine the two point angular correlation function in each bin. We find that the amplitude of the correlation, A_w, drops across the three bins to redshift z~1. Simple epsilon models of clustering evolution fit this result, with the best agreement for epsilon=0. Hierarchical cold-dark-matter models best fit in a low density, Lambda-dominated universe.

We present galaxy counts in the near-ultraviolet (NUV; 2365\AA) and far-ultraviolet (FUV; 1595\AA) obtained from Space Telescope Imaging Spectrograph (STIS) observations of portions of the Hubble Deep Field North (HDFN), the Hubble Deep Field South (HDFS), and a parallel field near the HDFN. All three fields have deep (AB$>$29) optical imaging, and we determine magnitudes by taking the ultraviolet flux detected within the limiting optical isophote. An analysis of the UV-optical colors of detected objects, combined with a visual inspection of the UV images, indicates that there are no detectable objects in the UV images that are not also detected in the optical. We determine our detection area and completeness as a function of magnitude with a simulated distribution of galaxies based on the HDFN Wide Field Planetary Camera 2 (WFPC2) V+I image. The galaxy counts reach to $NUV_{AB}=29$ and $FUV_{AB}=30$, which is 1 magnitude fainter than the HDFN WFPC2/F300W counts, and 7 magnitudes fainter than balloon-based counts. We integrate our measurements to present the extragalactic background radiation coming from resolved objects. The NUV galaxy counts show a flattening or turnover beginning at about $NUV_{AB}=26.0$, which is not predicted either by no-evolution models based upon a local luminosity function with a steep faint end slope, nor by a semi-analytic model in which starbursts are caused by major mergers. The FUV counts also show a flat slope. We argue that the flat slopes could be caused by a short duty cycle for star formation, additional starbursts triggered by minor mergers, and an extended quiescent phase between starburst episodes.

We examine the cooling radiation from forming galaxies in hydrodynamic simulations of the LCDM model (cold dark matter with a cosmological constant), focusing on the Ly-alpha line luminosities of high-redshift systems. Primordial composition gas condenses within dark matter potential wells, forming objects with masses and sizes comparable to the luminous regions of observed galaxies. As expected, the energy radiated in this process is comparable to the gravitational binding energy of the baryons, and the total cooling luminosity of the galaxy population peaks at z ~= 2. However, in contrast to the classical picture of gas cooling from the \sim 10^6 K virial temperature of a typical dark matter halo, we find that most of the cooling radiation is emitted by gas with T < 20,000 K. As a consequence, roughly 50% of this cooling radiation emerges in the Ly-alpha line. While a galaxy's cooling luminosity is usually smaller than the ionizing continuum luminosity of its young stars, the two are comparable in the most massive systems, and the cooling radiation is produced at larger radii, where the Ly-alpha photons are less likely to be extinguished by dust. We suggest, in particular, that cooling radiation could explain the two large (\sim 100 kpc), luminous (L_{Ly-alpha} \sim 10^{44} erg s^{-1}) ``blobs'' of Ly-alpha emission found in Steidel et al.'s (1999) narrow band survey of a z = 3 proto-cluster. Our simulations predict objects of the observed luminosity at about the right space density, and radiative transfer effects can account for the observed sizes and line widths. We discuss observable tests of this hypothesis for the nature of the Ly-alpha blobs, and we present predictions for the contribution of cooling radiation to the Ly-alpha luminosity function of galaxies as a function of redshift.

We discuss the utility of a low resolution prism as a component of a Multi-Object Spectrometer for NASA's proposed Next Generation Space Telescope (NGST). Low resolution prism spectroscopy permits simultaneous observation of the 0.6-5micron wavelength regime at R~50. To such data we can apply the modern techniques in spectral energy distribution (SED) fitting to determine source redshifts, sometimes called ``photometric redshifts''. Low resolution prism observations of galaxy SED's provide a significant advantage over multi-filter observations for any realistic observing strategy, and we have performed extensive simulations to quantify it. The advantage of the prism is largest in the regime of faint (K_AB>30) objects at high redshift (z>4), the primary discovery space of NGST. Many important objects will be too faint for follow up at higher spectral resolution, so prism observations are the optimal technique to study them. Prism observations also reduce the contamination of high redshift samples by lower redshift interlopers.

We present candidate quasars from a multi-color (uBVRI + narrow-band) imaging survey of 1/2 square degree around the Hubble Deep Field - South. We identify 154 candidate quasars with B < 23 using color selection, consistent with previously measured QSO number counts if we assume a 60% selection efficiency. The narrow-band filter (NB) was centered at 3958 A to detect Ly-alpha at the redshift of J2233-6033, the HDF-S QSO. We confirm the presence of Ly-alpha nebulosity extending ~12'' around the HDF-S QSO, reported by Bergeron et al. (1999). We detect 10 point-like objects in emission through the NB filter. Of these, 7 satisfy our QSO color selection criteria. One of emission-line objects is a B~20 radio-quiet quasar at z=1.56, 6.7' from the line of sight to the HDF-S QSO and ~12'' from the western edge of the WFPC2 deep field.

Far-UV observations in and near the Hubble Deep Fields demonstrate that the Space Telescope Imaging Spectrograph (STIS) can potentially obtain unique and precise measurements of the diffuse far-ultraviolet background. Although STIS is not the ideal instrument for such measurements, high-resolution images allow Galactic and extragalactic objects to be masked to very faint magnitudes, thus ensuring a measurement of the truly diffuse UV signal. The programs we have analyzed were not designed for this scientific purpose, but would be sufficient to obtain a very sensitive measurement if it were not for a weak but larger-than-expected signal from airglow in the STIS 1450-1900 A bandpass. Our analysis shows that STIS far-UV crystal quartz observations taken near the limb during orbital day can detect a faint airglow signal, most likely from NI\1493, that is comparable to the dark rate and inseparable from the far-UV background. Discarding all but the night data from these datasets gives a diffuse far-ultraviolet background measurement of 501 +/- 103 ph/cm2/sec/ster/A, along a line of sight with very low Galactic neutral hydrogen column (N_HI = 1.5E20 cm-2) and extinction (E(B-V)=0.01 mag). This result is in good agreement with earlier measurements of the far-UV background, and should not include any significant contribution from airglow. We present our findings as a warning to other groups who may use the STIS far-UV camera to observe faint extended targets, and to demonstrate how this measurement may be properly obtained with STIS.

Science objectives for the Next Generation Space Telescope (NGST) include a large component of galaxy surveys, both imaging and spectroscopy. The Hubble Deep Field datasets include the deepest observations ever made in the ultraviolet, optical and near infrared, reaching depths comparable to that expected for NGST spectroscopy. We present the source counts, galaxy sizes and isophotal filling factors of the HDF-South images. The observed integrated galaxy counts reach >500 galaxies per square arcminute at AB<30. We extend these counts to faint levels in the infrared using models. The trend previously seen that fainter galaxies are smaller, continues to AB=29 in the high resolution HDF-S STIS image, where galaxies have a typical half-light radius of 0.1 arcseconds. Extensive Monte Carlo simulations show that the small measured sizes are not due to selection effects until >29mag. Using the HDF-S NICMOS image, we show that galaxies are smaller in the near infrared than they are in the optical. We analyze the isophotal filling factor of the HDF-S STIS image, and show that this image is mostly empty sky even at the limits of galaxy detection, a conclusion we expect to hold true for NGST spectroscopy. At the surface brightness limits expected for NGST imaging, however, about a quarter of the sky is occupied by the outer isophotes of AB<30 galaxies. We discuss the implications of these data on several design concepts of the NGST near-infrared spectrograph. We compare the effects of resolution and the confusion limit of various designs, as well as the multiplexing advantages of either multi-object or full-field spectroscopy. We argue that the optimal choice for NGST spectroscopy of high redshift galaxies is a multi-object spectrograph (MOS) with target selection by a micro electro mechanical system (MEMS) device.

We present the imaging observations made with the Space Telescope Imaging Spectrograph of the Hubble Deep Field - South. The field was imaged in 4 bandpasses: a clear CCD bandpass for 156 ksec, a long-pass filter for 22-25 ksec per pixel typical exposure, a near-UV bandpass for 23 ksec, and a far-UV bandpass for 52 ksec. The clear visible image is the deepest observation ever made in the UV-optical wavelength region, reaching a 10 sigma AB magnitude of 29.4 for an object of area 0.2 square arcseconds. The field contains QSO J2233-606, the target of the STIS spectroscopy, and extends 50"x50" for the visible images, and 25"x25" for the ultraviolet images. We present the images, catalog of objects, and galaxy counts obtained in the field.

What type of objects are being detected as $z\sim 3$ "Lyman break galaxies"? Are they predominantly the most massive galaxies at that epoch, or are many of them smaller galaxies undergoing a short-lived burst of merger-induced star formation? We attempt to address this question using high-resolution cosmological hydrodynamic simulations including star formation and feedback. Our $\Lambda$CDM simulation, together with Bruzual-Charlot population synthesis models, reproduces the observed number density and luminosity function of Lyman break galaxies when dust is incorporated. The inclusion of dust is crucial for this agreement. In our simulation, these galaxies are predominantly the most massive objects at this epoch, and have a significant population of older stars. Nevertheless, it is possible that our simulations lack the resolution and requisite physics to produce starbursts, despite having a physical resolution of $\la 700$ pc at z=3. Thus we cannot rule out merger-induced starburst galaxies also contributing to the observed population of high-redshift objects.

We review theoretical approaches to the study of galaxy formation, with emphasis on the role of hydrodynamic simulations in modeling the high redshift galaxy population. We present new predictions for the abundance of star-forming galaxies in the Lambda + cold dark matter model (Omega_m=0.4, Omega_L=0.6), combining results from several simulations to probe a wide range of redshift. At a threshold density of one object per arcmin^2 per unit z, these simulations predict galaxies with star formation rates of 2 msun/yr (z=10), 5 msun/yr (z=8), 20 msun/yr (z=6), 70-100 msun/yr (z=4-2), and 30 msun/yr (z=0.5). For galaxies selected at a fixed comoving space density n=0.003 h^3 Mpc^{-3], a (50 Mpc/h)^3 simulation predicts a galaxy correlation function (r/5 Mpc/h)^{-1.8} in comoving coordinates, essentially independent of redshift from z=4 to z=0.5. Different cosmological models predict global histories of star formation that reflect their overall histories of mass clustering, but robust numerical predictions of the comoving space density of star formation are difficult because the simulations miss the contribution from galaxies below their resolution limit. The LCDM model appears to predict a star formation history with roughly the shape inferred from observations, but it produces too many stars at low redshift, predicting Omega_* ~ 0.015 at z=0. We conclude with a brief discussion of this discrepancy and three others that suggest gaps in our current theory of galaxy formation: small disks, steep central halo profiles, and an excess of low mass dark halos. While these problems could fade as the simulations or observations improve, they could also guide us towards a new understanding of galactic scale star formation, the spectrum of primordial fluctuations, or the nature of dark matter.

(abridged) We present results from two new near-infrared imaging surveys. One covers 47.2 arcmin^2 to K(3\sigma)<20 whilst a second, deeper survey covers a sub-area of 1.8 arcmin^2 to K(3\sigma)<22.75. Over the entire area we have extremely deep UBRI photometry. Our K- counts are consistent with the predictions of non-evolving models with 0 < q0 <0.5. The K-selected (B-K) galaxy colour distributions move sharply bluewards fainter than K~20 and at at brighter magnitudes (K<20) our observed colour distributions indicate a deficiency of red, early-type galaxies at z~1 in comparison with passively evolving models. This implies either a pure luminosity evolution (PLE) model with a low level of continuing star-formation following an an initial burst, or dynamical merging. At fainter magnitudes, the continuing bluewards trend observed in (B-K) can be explained purely in terms of passively evolving PLE models. Our observed numbers of (I-K)>4 galaxies at K<20 exhibit the same deficiency, suggesting that at least part of the larger deficit observed in (B-K) at K<20 may be due to star-formation rather than dynamical merging. Finally, as we and others have noted, the number-redshift distribution at 18<K<19 of recent, deep K- selected redshift surveys is well fitted by non-evolving models; passively evolving models with a Salpeter or Scalo initial mass functions overpredict the numbers of galaxies with z>1. Dynamical merging is one possible solution to reduce the numbers of these galaxies but a dwarf-dominated IMF for early-type galaxies could offer an alternative explanation; we show that this model reproduces both the optical-infrared colour distributions and the K- band galaxy counts.

We present the first results of a survey of random fields with the slitless G141 ($\lambda_c = 1.5\mu, \Delta\lambda=0.8\mu$) grism on NICMOS. Approximately 64 square arcminutes have been observed at intermediate and high galactic latitudes. The 3$\sigma$ limiting line and continuum fluxes in each field vary from $7.5 \times 10^{-17}$ to $1 \times 10^{-17} erg/cm^2/sec$ and from H = 20 to 22, respectively. Our median and area weighted $3\sigma $ limiting line fluxes within a 4 pixel aperture are nearly identical at $4.1 \times 10^{-17} erg/cm^2/sec$ and are 60% deeper than the deepest narrow-band imaging surveys from the ground. We have identified 33 emission-line objects and derive their observed wavelengths, fluxes and equivalent widths. We argue that the most likely line identification is H$\alpha$ and that the redshift range probed is from 0.75 to 1.9. The 2$\sigma$ rest-frame equivalent width limits range from 9\AA to 130\AA with an average of 40\AA. The survey probes an effective co-moving volume of $10^5 h_{50}^{-3} Mpc^3$ for $q_0=0.5$. Our derived co-moving number density of emission line galaxies in the range $0.7 < z < 1.9$ is $3.3\times10^{-4} h_{50}^{3} Mpc^{-3}$, very similar to that of the bright Lyman break objects at $z \sim 3$. The objects with detected emission-lines have a median F160W magnitude of 20.4 (Vega scale) and a median H$\alpha$ luminosity of $2.7 \times 10^{42} erg/sec$. The implied star formation rates range from 1 to 324 M_{\odot}/yr, with an average [NII]6583,6548 corrected rate of 21 M_{\odot}/yr for H_0=50 km/s/Mpc and $q_0=0.5$ (34 M_{\odot}/yr for $q_0=0.1$).

We present high resolution NICMOS images of random fields obtained in parallel to other HST observations. We present galaxy number counts reaching H=24. The H-band galaxy counts show good agreement with the deepest I- and K-band counts obtained from ground-based data. We present the distribution of galaxies with morphological type to H<23. We find relatively fewer irregular galaxies compared to an I-band sample from the Hubble Deep Field, which we attribute to their blue color, rather than to morphological K-corrections. We conclude that the irregulars are intrinsically faint blue galaxies at z<1.

I review galaxy number count models and present ncmod, an extendable and general purpose model for comparing and interpreting the results of field galaxy survey data. I develop techniques and software for converting the results of a survey done in one filter into another filter, for direct comparison with other surveys. Comparison of the data from surveys which differ greatly in wavelength coverage or sensitivity is of necessity model-dependent, but comparison between similar surveys can be done in a relatively model-independent way. I extrapolate existing number counts into the ultraviolet and thermal infrared. The model is used to predict the results of future space missions, including STIS and NICMOS on HST, ISO, SIRTF and NGST.

We compare the galaxy evolution models of Bruzual & Charlot (1993) with the faint galaxy count, size and colour data from the Hubble and Herschel Deep Fields (Metcalfe et al 1996). For qo=0.05, we find that models where the SFR increases exponentially out to z>2 are consistent with all of the observational data. For qo=0.5, such models require an extra population of galaxies which are only seen at high redshift and then rapidly fade or disappear. We find that, whatever the cosmology, the redshift of the faint blue galaxies and hence the epoch of galaxy formation is likely to lie at z>2. We find no implied peak in the SFR at z=1 and we suggest that the reasons for this contradiction with the results of Madau et al (1996) include differences in faint galaxy photometry, in the treatment of spiral dust and in the local galaxy count normalisation.

An echelle spectrogram (R = 30,000) of the 2300-3100 A region in the ultraviolet spectrum of the F8V star 9 Comae is presented. The observation is used to calibrate features in the mid-ultraviolet spectra of similar stars according to age and metal content. In particular, the spectral break at 2640 A is interpreted using the spectral synthesis code SYNSPEC. We use this feature to estimate the time since the last major star formation episode in the z=1.55 early-type galaxy LBDS 53W091, whose rest frame mid-ultraviolet spectrum, observed with the Keck Telescope, is dominated by the flux from similar stars that are at or near the main-sequence turnoff in that system (Spinrad et al. 1997). Our result, 1 Gyr if the flux-dominating stellar population has a metallicity twice solar, or 2 Gyr for a more plausible solar metallicity, is significantly lower than the previous estimate and thereby relaxes constraints on cosmological parameters that were implied by the earlier work.

The installation of the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope (HST) allows for the first time two-dimensional optical and ultraviolet slitless spectroscopy of faint objects from space. The STIS Parallel Survey (SPS) routinely obtains broad band images and slitless spectra of random fields in parallel with HST observations using other instruments. The SPS is designed to study a wide variety of astrophysical phenomena, including the rate of star formation in galaxies at intermediate to high redshift through the detection of emission-line galaxies. We present the first results of the SPS, which demonstrate the capability of STIS slitless spectroscopy to detect and identify high-redshift galaxies.

Using cosmological N-body simulations of critical (SCDM) and open (Omega_0=0.3, OCDM) cold dark matter models we have identified dark matter halos which are associated with the progenitors of present day bright early type galaxies. By following their merging history, we show how early type galaxies that formed within massive halos at redshift around 3 are now preferentially residing in clusters and groups. On the other hand, those that formed through major merging events between redshift 1 and the present have not yet been accreted into larger, virialized structures. This result is in agreement with analytical predictions in hierarchical clustering models. CDM models are able to explain both the ancient and uniform population of ellipticals that dominates in clusters together with the more recent and heterogeneous population of field ellipticals. Predictions for the comoving number density of bright early-type galaxies are given, and are shown to be consistent with the observed luminosity function. We predict that the number density of interacting bright binary galaxies, from which the field population of ellipticals may have originated, is proportional to (1+z)**4.2+-0.28 and (1+z)**2.5+-0.42 in SCDM and OCDM respectively. This result is consistent with previous analytical estimates and is discussed together with recent observational constraints.

We present an optical and near-infrared survey of galaxies in nearby clusters aimed at determining fundamental quantities of galaxies, such as multivariate luminosity function and color distribution for each Hubble type. The main characteristics of our survey are completeness in absolute magnitude, wide wavelength coverage and faint limiting magnitudes.

We calculate the number of damped Lyman-alpha absorbers expected in various popular cosmological models as a function of redshift and compare our predictions with observed abundances. The Press-Schechter formalism is used to obtain the distribution of halos with circular velocity in different cosmologies, and we calibrate the relation between circular velocity and absorption cross-section using detailed gas dynamical simulations of a ``standard'' cold dark matter (CDM) model. Because of this calibration, our approach makes more realistic assumptions about the absorption properties of collapsed objects than previous, analytic calculations of the damped Lyman-alpha abundance. CDM models with Omega_0=1, H_0=50, baryon density Omega_b=0.05, and scale-invariant primeval fluctuations reproduce the observed incidence and redshift evolution of damped Lyman-alpha absorption to within observational uncertainty, for both COBE normalization (sigma_8=1.2) and a lower normalization (sigma_8=0.7) that better matches the observed cluster abundance at z=0. A tilted (n=0.8, sigma_8=0.7) CDM model tends to underproduce absorption, especially at z=4. With COBE normalization, a CDM model with Omega_0=0.4, Omega_{Lambda}=0.6 gives an acceptable fit to the observed absorption; an open CDM model is marginally acceptable if Omega_0 is at least 0.4 and strongly inconsistent with the z=4 data if Omega_0=0.3. Mixed dark matter models tend not to produce sufficient absorption, being roughly comparable to tilted CDM models if Omega_{nu} = 0.2 and failing drastically if Omega_{nu} = 0.3.

We present the first determination of the near-infrared K-band luminosity function of field galaxies from a wide field K-selected redshift survey. The best fit Schechter function parameters are $M^* = -23.12 +5log(h)$, $\alpha = -0.91$, and $\phi^* = 1.66 \times 10^{-2} h^3 ~ Mpc^{-3}$. We estimate that systematics are no more than 0.1 mag in $M^*$ and 0.1 in $\alpha$, which is comparable to the statistical errors on this measurement.

We present the results of a wide-field K-selected galaxy survey with complementary optical I- and B-band imaging in six fields with a total coverage of 9.8 square degrees. This survey establishes the bright-end K-band galaxy number counts in the magnitude range 13<K<16 with high precision. We find that our bright-end counts have a significantly steeper slope than the prediction of a no-evolution model, which cannot be accounted for by known observational or theoretical error. We also argue against the likelihood of sufficient evolution at such low redshifts to account for this effect, we describe an alternative picture in which there is a local deficiency of galaxies by a factor of 2 on scale sizes of around 300 Mpc. Taken at face value, this would imply that local measurements of \Omega_0 underestimate the true value of the cosmological mass density by this factor and that local measurements of $H_0$ could be high by as much as 33%.

Sensitive optical surveys have revealed a large population of `faint blue galaxies' which are believed to be young galaxies, observed close to their time of formation. But there has been considerable uncertainty regarding the epochs at which these galaxies are observed, owing to the difficulties inherent in determining spectroscopic redshifts for very faint objects. Here, by modelling the counts and colours of galaxies at the faintest detection limits, we show that the faint blue galaxies are likely to lie at high redshift (z approx. 2). This conclusion holds regardless of whether the Universe is assumed to be open or at the critical density (flat). In an open universe, the data are consistent with galaxy models in which star formation rates decay exponentially with decreasing redshift, whereas the assumption of a flat universe requires the addition of a population of galaxies which are seen only at high redshift.

Lyman limit and damped Lyman-alpha absorption systems probe the distribution of collapsed, cold gas at high redshift. Numerical simulations that incorporate gravity and gas dynamics can predict the abundance of such absorbers in cosmological models. We develop a semi-analytical method to correct the numerical predictions for the contribution of unresolved low mass halos, and we apply this method to the Katz et al. (1996) simulation of the standard cold dark matter model ($\Omega=1$, $h=0.5$, $\Omega_b=0.05$, $\sigma_8=0.7$). Using this simulation and higher resolution simulations of individual low mass systems, we determine the relation between a halo's circular velocity $v_c$ and its cross section for producing Lyman limit or damped absorption. We combine this relation with the Press-Schechter formula for the abundance of halos to compute the number of absorbers per unit redshift. The resolution correction increases the predicted abundances by about a factor of two at z=2, 3, and 4, bringing the predicted number of damped absorbers into quite good agreement with observations. Roughly half of the systems reside in halos with circular velocities $v_c>100\kms$ and half in halos with $35\kms<v_c<100\kms$. Halos with $v_c>150\kms$ typically harbor two or more systems capable of producing damped absorption. Even with the resolution correction, the predicted abundance of Lyman limit systems is a factor of three below observational estimates, signifying either a failure of standard CDM or a failure of these simulations to resolve the systems responsible for most Lyman limit absorption. By comparing simulations with and without star formation, we find that depletion of the gas supply by star formation affects absorption line statistics at $z>=2$ only for column densities exceeding $N_{HI}=10^{22} cm^{-2}$.

We present the first measurement of the angular correlation function in a $K$-selected galaxy survey, from two fields covering almost 10 square degrees. The angular correlation function at small angles is well described by a $\theta^{-0.8}$ power law, as for optically selected samples. The clustering amplitude is reduced by a factor of $\sim 4$ between galaxies with $K<15$ and a fainter magnitude slice with $15<K<16$. This allows us to place constraints upon the redshift distribution of the galaxies and their spatial correlation function. We find no clear evidence for a change in clustering amplitude when galaxies are selected by their observed $B-K$ colours.

We present bright galaxy number counts measured with linear detectors in the B, V, I, and K bands in two fields covering nearly 10 square degrees. All of our measurements are consistent with passive evolution models, and do not confirm the steep slope measured in other surveys at bright magnitudes. Throughout the range 16 < B < 19, our B-band counts are consistent with the "high normalization" models proposed to reduce the faint blue galaxy problem. Our K-band counts agree with previous measurements, and have reached a fair sample of the universe in the magnitude range where evolution and K-corrections are well understood.