This report provides detailed findings on the critical laboratory astrophysics data needs that are required to maximize the scientific return for NASA's current and near-term planned astrophysics missions. It also provides prioritized rankings on said laboratory astrophysics data, generally by waveband. The Report is based on community input gathered at the 2018 NASA Laboratory Astrophysics Workshop (LAW) from presentations, from discussions during workshop breakout sessions, and from other solicited input deemed appropriate by the Scientific Organizing Committee (SOC) obtained prior to and after the meeting. Hence, the Report is a direct reflection of the spirit and participant make-up of LAW 2018. The Report also outlines specific opportunities and threats facing NASA's Laboratory Astrophysics Program, and articulates concrete actions by which the Agency can capitalize on the opportunities and mitigate the challenges. The Report was prepared by the SOC, with help from some invited speakers, and input and review from community members.

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

We describe the near real-time transient-source discovery engine for the intermediate Palomar Transient Factory (iPTF), currently in operations at the Infrared Processing and Analysis Center (IPAC), Caltech. We coin this system the IPAC/iPTF Discovery Engine (or IDE). We review the algorithms used for PSF-matching, image subtraction, detection, photometry, and machine-learned (ML) vetting of extracted transient candidates. We also review the performance of our ML classifier. For a limiting signal-to-noise ratio of 4 in relatively unconfused regions, "bogus" candidates from processing artifacts and imperfect image subtractions outnumber real transients by ~ 10:1. This can be considerably higher for image data with inaccurate astrometric and/or PSF-matching solutions. Despite this occasionally high contamination rate, the ML classifier is able to identify real transients with an efficiency (or completeness) of ~ 97% for a maximum tolerable false-positive rate of 1% when classifying raw candidates. All subtraction-image metrics, source features, ML probability-based real-bogus scores, contextual metadata from other surveys, and possible associations with known Solar System objects are stored in a relational database for retrieval by the various science working groups. We review our efforts in mitigating false-positives and our experience in optimizing the overall system in response to the multitude of science projects underway with iPTF.

We examine the repeatability, reliability, and accuracy of differential exoplanet eclipse depth measurements made using the InfraRed Array Camera (IRAC) on the Spitzer Space Telescope during the post-cryogenic mission. We have re-analyzed an existing 4.5 {\mu}m data set, consisting of 10 observations of the XO-3b system during secondary eclipse, using seven different techniques for removing correlated noise. We find that, on average, for a given technique, the eclipse depth estimate is repeatable from epoch to epoch to within 156 parts per million (ppm). Most techniques derive eclipse depths that do not vary by more than a factor 3 of the photon noise limit. All methods but one accurately assess their own errors: for these methods, the individual measurement uncertainties are comparable to the scatter in eclipse depths over the 10 epoch sample. To assess the accuracy of the techniques as well as to clarify the difference between instrumental and other sources of measurement error, we have also analyzed a simulated data set of 10 visits to XO-3b, for which the eclipse depth is known. We find that three of the methods (BLISS mapping, Pixel Level Decorrelation, and Independent Component Analysis) obtain results that are within three times the photon limit of the true eclipse depth. When averaged over the 10 epoch ensemble, 5 out of 7 techniques come within 60 ppm of the true value. Spitzer exoplanet data, if obtained following current best practices and reduced using methods such as those described here, can measure repeatable and accurate single eclipse depths, with close to photon-limited results.

The Frontier Fields are a director's discretionary time campaign with HST and the Spitzer Space Telescope to see deeper into the universe than ever before. The Frontier Fields combine the power of HST and Spitzer with the natural gravitational telescopes of massive high-magnification clusters of galaxies to produce the deepest observations of clusters and their lensed galaxies ever obtained. Six clusters - Abell 2744, MACSJ0416.1-2403, MACSJ0717.5+3745, MACSJ1149.5+2223, Abell S1063, and Abell 370 - were selected based on their lensing strength, sky darkness, Galactic extinction, parallel field suitability, accessibility to ground-based facilities, HST, Spitzer and JWST observability, and pre-existing ancillary data. These clusters have been targeted by the HST ACS/WFC and WFC3/IR with coordinated parallels of adjacent blank fields for over 840 HST orbits. The Spitzer Space Telescope has dedicated > 1000 hours of director's discretionary time to obtain IRAC 3.6 and 4.5 micron imaging to ~26.5, 26.0 ABmag 5-sigma point-source depths in the six cluster and six parallel Frontier Fields. The Frontier Field parallel fields are the second-deepest observations thus far by HST with ~29th ABmag 5-sigma point source depths in seven optical - near-infrared bandpasses. Galaxies behind the Frontier Field cluster lenses experience typical magnification factors of a few, with small regions near the critical curves magnified by factors 10-100. Therefore, the Frontier Field cluster HST images achieve intrinsic depths of ~30-33 magnitudes over very small volumes. Early studies of the Frontier Fields have probed galaxies fainter than any seen before during the epoch of reionization 6 < z < 10, mapped out the cluster dark matter to unprecedented resolution, and followed lensed transient events.

We present luminosity functions derived from a spectroscopic survey of AGN selected from Spitzer Space Telescope imaging surveys. Selection in the mid-infrared is significantly less affected by dust obscuration. We can thus compare the luminosity functions of the obscured and unobscured AGN in a more reliable fashion than by using optical or X-ray data alone. We find that the AGN luminosity function can be well described by a broken power-law model in which the break luminosity decreases with redshift. At high redshifts ($z>1.6$), we find significantly more AGN at a given bolometric luminosity than found by either optical quasar surveys or hard X-ray surveys. The fraction of obscured AGN decreases rapidly with increasing AGN luminosity, but, at least at high redshifts, appears to remain at $\approx 50$\% even at bolometric luminosities $\sim 10^{14}L_{\odot}$. The data support a picture in which the obscured and unobscured populations evolve differently, with some evidence that high luminosity obscured quasars peak in space density at a higher redshift than their unobscured counterparts. The amount of accretion energy in the Universe estimated from this work suggests that AGN contribute about 12\% to the total radiation intensity of the Universe, and a high radiative accretion efficiency $\approx 0.18^{+0.12}_{-0.07}$ is required to match current estimates of the local mass density in black holes.

We present the results of a program of optical and near-infrared spectroscopic follow-up of candidate Active Galactic Nuclei (AGN) selected in the mid-infrared. This survey selects both normal and obscured AGN closely matched in luminosity across a wide range, from Seyfert galaxies with bolometric luminosities L_bol~10^10L_sun, to highly luminous quasars (L_bol~10^14L_sun), and with redshifts from 0-4.3. Samples of candidate AGN were selected through mid-infrared color cuts at several different 24 micron flux density limits to ensure a range of luminosities at a given redshift. The survey consists of 786 candidate AGN and quasars, of which 672 have spectroscopic redshifts and classifications. Of these, 137 (20%) are type-1 AGN with blue continua, 294 (44%) are type-2 objects with extinctions A_V>~5 towards their AGN, 96 (14%) are AGN with lower extinctions (A_V~1) and 145 (22%) have redshifts, but no clear signs of AGN activity in their spectra. 50% of the survey objects have L_bol >10^12L_sun, in the quasar regime. We present composite spectra for type-2 quasars and for objects with no signs of AGN activity in their spectra. We also discuss the mid-infrared - emission-line luminosity correlation and present the results of cross-correlations with serendipitous X-ray and radio sources. The results show that: (1) obscured objects dominate the overall AGN population, (2) there exist mid-infrared selected AGN candidates which lack AGN signatures in their optical spectra, but have AGN-like X-ray or radio counterparts, and (3) X-ray and optical classifications of obscured and unobscured AGN often differ.

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

We model the optical to far-infrared SEDs of a sample of six type-1 and six type-2 quasars selected in the mid-infrared. The objects in our sample are matched in mid-IR luminosity and selected based on their Spitzer IRAC colors. We obtained new targeted Spitzer IRS and MIPS observations and used archival photometry to examine the optical to far-IR SEDs. We investigate whether the observed differences between samples are consistent with orientation-based unification schemes. The type-1 objects show significant emission at 3 micron. They do not show strong PAH emission and have less far-IR emission on average when compared to the type-2 objects. The SEDs of the type-2 objects show a wide assortment of silicate features, ranging from weak emission to deep silicate absorption. Some also show strong PAH features. In comparison, silicate is only seen in emission in the type-1 objects. This is consistent with some of the type-2s being reddened by a foreground screen of cooler dust, perhaps in the host galaxy itself. We investigate the AGN contribution to the far-IR emission and find it to be significant. We also estimate the star formation rate for each of the objects by integrating the modeled far-IR flux and compare this with the SFR found from PAH emission. We find the type-2 quasars have a higher average SFR than the type-1 quasars based on both methods, though this could be due to differences in bolometric luminosities of the objects. While we find pronounced differences between the two types of objects, none of them are inconsistent with orientation-based unification schemes.

We explore the multidimensional, multiwavelength selection of quasars from mid-IR (MIR) plus optical data, specifically from Spitzer-IRAC and the Sloan Digital Sky Survey (SDSS). We apply modern statistical techniques to combined Spitzer MIR and SDSS optical data, allowing up to 8-D color selection of quasars. Using a Bayesian selection method, we catalog 5546 quasar candidates to an 8.0 um depth of 56 uJy over an area of ~24 sq. deg; ~70% of these candidates are not identified by applying the same Bayesian algorithm to 4-color SDSS optical data alone. Our selection recovers 97.7% of known type 1 quasars in this area and greatly improves the effectiveness of identifying 3.5<z<5 quasars. Even using only the two shortest wavelength IRAC bandpasses, it is possible to use our Bayesian techniques to select quasars with 97% completeness and as little as 10% contamination. This sample has a photometric redshift accuracy of 93.6% (Delta Z +/-0.3), remaining roughly constant when the two reddest MIR bands are excluded. While our methods are designed to find type 1 (unobscured) quasars, as many as 1200 of the objects are type 2 (obscured) quasar candidates. Coupling deep optical imaging data with deep mid-IR data could enable selection of quasars in significant numbers past the peak of the quasar luminosity function (QLF) to at least z~4. Such a sample would constrain the shape of the QLF and enable quasar clustering studies over the largest range of redshift and luminosity to date, yielding significant gains in our understanding of quasars and the evolution of galaxies.

We report the first direct detection with Spitzer of galaxy filaments. Using Spitzer and ancillary optical data, we have discovered two filamentary structures in the outskirts of the cluster Abell 1763. Both filaments point toward Abell 1770 which lies at the same redshift as Abell 1763 (z=0.23), at a projected distance of ~13 Mpc. The X-ray cluster emission is elongated along the same direction. Most of the far-infrared emission is powered by star formation. According to the optical spectra, only one of the cluster members is classified as an active galactic nucleus. Star formation is clearly enhanced in galaxies along the filaments: the fraction of starburst galaxies in the filaments is more than twice than that in other cluster regions. We speculate that these filaments are feeding the cluster Abell 1763 by the infall of galaxies and galaxy groups. Evidence for one of these groups is provided by the analysis of galaxy kinematics in the central cluster region.

We present Hubble Space Telescope images, and spectral energy distributions from optical to infrared wavelengths for a sample of six 0.3<z<0.8 type-2 quasars selected in the mid-infrared using data from the Spitzer Space Telescope. All the host galaxies show some signs of disturbance. Most seem to possess dusty, star-forming disks. The disk inclination, estimated from the axial ratio of the hosts, correlates with the depth of the silicate feature in the mid-infrared spectra, implying that at least some of the reddening towards the AGN arises in the host galaxy. The star formation rates in these objects, as inferred from the strengths of the PAH features and far-infrared continuum, range from 3-90 Msun/yr, but are mostly much larger than those inferred from the [OII]3727 emission line luminosity, due to obscuration. Taken together with studies of type-2 quasar hosts from samples selected in the optical and X-ray, this is consistent with previous suggestions that two types of extinction processes operate within the type-2 quasar population, namely a component due to the dusty torus in the immediate environment of the AGN, and a more extended component due to a dusty, star forming disk.

We investigate the properties of 1088 Lyman Break Galaxies (LBGs) at z~3 selected from a ~2.63$ deg2 sub-region of the First Look Survey field using the ground-based multi-color data and the Spitzer Space Telescope mid-infrared data at 3--8 and 24 um. With the wide area and the broad wavelength coverage, we sample a large number of ``rare'' u-band dropouts which are massive (M* > 10^11 Msun), allowing us to perform a statistical analysis of these subsets of LBGs that have not been studied in detail. Optically bright (R(AB) < 24.5 mag) LBGs detected in mid-infrared (S_{3.6um} > 6 uJy) reside at the most massive and dusty end of the LBG population, with relatively high and tight $M/L$ in rest-frame near-infrared. Most infrared-luminous LBGs (S_{24um} > 100 uJy) are dusty star-forming galaxies with star formation rates of 100--1000 Msun/yr, total infrared luminosity of > 10^12 Lsun. By constructing the UV luminosity function of massive LBGs, we estimate that the lower limit for the star formation rate density from LBGs more massive than 10^11 Msun at z~3 is > 3.3 x 10^-3 Msun/yr/Mpc^3, showing for the first time that the UV-bright population of massive galaxies alone contributes significantly to the global star formation rate density at z~3. When combined with the star formation rate densities at z < 2, our result reveals a steady increase in the contribution of massive galaxies to the global star formation from z=0 to z=3, providing strong support to the downsizing of galaxy formation.

We present deep u*-, and g-band images taken with the MegaCam on the 3.6 m Canada-France-Hawaii Telescope (CFHT) to support the extragalactic component of the Spitzer First Look Survey (hereafter, FLS). In this paper we outline the observations, present source catalogs and characterize the completeness, reliability, astrometric accuracy and number counts of this dataset. In the central 1 deg2 region of the FLS, we reach depths of g~26.5 mag, and u*~26.2 mag (AB magnitude, 5$\sigma$ detection over a 3" aperture) with ~4 hours of exposure time for each filter. For the entire FLS region (~5 deg2 coverage), we obtained u*-band images to the shallower depth of u*=25.0--25.4 mag (5$\sigma$, 3" aperture). The average seeing of the observations is 0.85" for the central field, and ~1.00" for the other fields. Astrometric calibration of the fields yields an absolute astrometric accuracy of 0.15" when matched with the SDSS point sources between 18<g<22. Source catalogs have been created using SExtractor. The catalogs are 50% complete and greater than 99.3% reliable down to g~26.5 mag and u*~26.2 mag for the central 1 deg2 field. In the shallower u*-band images, the catalogs are 50% complete and 98.2% reliable down to 24.8--25.4 mag. These images and source catalogs will serve as a useful resource for studying the galaxy evolution using the FLS data.

We present an analysis of the effects of luminosity on the shape of the mid-infrared spectral energy distributions (SEDs) of 234 radio-quiet quasars originally presented by Richards et al. In quasars without evident dust extinction, the spectrally integrated optical and infrared luminosities are linearly correlated over nearly three decades in luminosity. We find a significant (>99.99% confidence) correlation between the 1.8-8.0 micron spectral index and infrared luminosity that indicates an enhancement of the mid-infrared continuum with increasing luminosity. Coupled with strong evidence for spectral curvature in more luminous quasars, we conclude this trend is likely a manifestation of the `near-infrared (3-5 micron) bump' noticed in earlier quasar SED surveys. The strength of this feature is indicative of the contribution of emission from the hottest (>1000 K) dust to the mid-infrared spectrum; higher luminosity quasars tend to show more hot dust emission. Finally, the comparable distribution of bolometric corrections from the monochromatic 3 micron luminosity as well as its lack of sensitivity to dust extinction as compared to the standard bolometric correction from nu*L_{5100A} suggest that the former may be a more robust indicator of bolometric quasar luminosity. The close link between the power in the mid-infrared and optical and the effect of luminosity on the shape of the mid-infrared continuum indicate that considering mid-infrared emission independent of the properties of the quasar itself is inadequate for understanding the parsec-scale quasar environment.

We present optical spectroscopy of a sample of 77 luminous AGN and quasars selected on the basis of their mid-infrared colors. Our objects are selected from the Spitzer Extragalactic First Look Survey and SWIRE XMM-LSS fields, with a typical 24mu flux density of 5mJy. The median redshift is 0.6, with a range of ~0.05-4. Only 33% (25/77) of these objects are normal type-1 quasars, with no obscuration. 44% (34/77) are type-2 objects, with high-ionization, narrow emission lines, 14% (11/77) are dust-reddened type-1 quasars, showing broad lines but a dust-reddened or unusually weak quasar continuum. 9% (7/77) show no sign of an AGN in the optical spectrum, having either starburst spectra or spectra which could be of either starburst or LINER types. These latter objects are analogous to the X-ray detected population of AGN with weak or non-existent optical AGN emission (the ``XBONGs''). 21 of our objects from the SWIRE field fall within moderately-deep XMM exposures. All the unobscured quasars, and about half the obscured quasars are detected in these exposures. This sample, when taken together with other samples of Spitzer selected AGN and quasars, and results from X-ray studies, confirms that obscured AGN dominate the AGN and quasar number counts of all rapidly-accreting supermassive black hole systems, at least for z~<4. This implies a high radiative efficiency for the black hole accretion process.

We present an analysis of the mid-infrared (MIR) and optical properties of type 1 (broad-line) quasars detected by the Spitzer Space Telescope. The MIR color-redshift relation is characterized to z=3, with predictions to z=7. We demonstrate how combining MIR and optical colors can yield even more efficient selection of active galactic nuclei (AGN) than MIR or optical colors alone. Composite spectral energy distributions (SEDs) are constructed for 259 quasars with both Sloan Digital Sky Survey and Spitzer photometry, supplemented by near-IR, GALEX, VLA and ROSAT data where available. We discuss how the spectral diversity of quasars influences the determination of bolometric luminosities and accretion rates; assuming the mean SED can lead to errors as large as a factor of 2 for individual quasars. Finally, we show that careful consideration of the shape of the mean quasar SED and its redshift dependence leads to a lower estimate of the fraction of reddened/obscured AGNs missed by optical surveys as compared to estimates derived from a single mean MIR to optical flux ratio.

We present a preliminary investigation of the spectral energy distributions (SEDs), and star-formation properties of a sample of Mid-IR selected Quasars. The mid-infrared SEDs of our objects are consistent with that expected from clumpy torus models. At longer wavelengths, the radio to infrared ratios of several objects are consistent with those of star-forming galaxies.

A multi-wavelength set of photometric data including UV (GALEX), optical, near-IR, infrared (Spitzer) and radio (VLA 20cm) images and spectroscopic observations are used to map the dust-obscured and unobscured star formation in the galaxy pair ARP 32. The system consists of an actively star-forming galaxy and another one with depressed star formation. The most active galaxy has disrupted morphology and different sites of star formation. Spectroscopic data show hints of nuclear activity in its core, intense star formation in limited regions of the galaxy as well as an underlying population of stars witnessing a past episode of star formation. Current star formation rates are estimated from UV and bolometric IR luminosities.

We present the reduction of the 24 micron data obtained during the first cosmological survey performed by the Spitzer Space Telescope (First Look Survey, FLS). The survey consists of a shallow observation of 2.5x2 sq deg centered at 17h18m +59d30m (main survey) and a deeper observation of 1x0.5 sq deg centered at 17h17m +59d45m(verification survey). Issues with the reduction of the 24 micron MIPS data are discussed and solutions to attenuate instrumental effects are proposed and applied to the data. Approximately 17000 sources are extracted with a SNR greater than five. The photometry of the point sources is evaluated through PSF fitting using an empirical PSF derived from the data. Aperture corrections and the absolute calibration have been checked using stars in the field. Astrometric and photometric errors depend on the SNR of the source varying between 0.35-1 arcsec and 5-15%, respectively, for sources detected at 20-5 sigma. The flux of the 123 extended sources have been estimated through aperture photometry. The extended sources cover less than 0.3% of the total area of the survey. Based on simulations, the main and verification surveys are 50% complete at 0.3 and 0.15 mJy, respectively. Counterparts have been searched for in optical and radio catalogs. More than 80% of the 24 micron sources have a reliable optical counterpart down to R=25.5. 16% of the sources have a 20 cm counterpart down to 0.1 mJy and ~ 80% of the radio-infrared associations have a reliable optical counterpart. A residual map is obtained by subtracting point sources detected at the 3-sigma level and interpolating the regions occupied by extended sources. Several galactic clouds with low and intermediate velocities are identified by comparison with neutral Hydrogen data from this field.

We investigate the instantaneous star formation rates (SFR) and extinction properties for a large (N=274), near-infrared (2.2 micron) + mid-infrared (24 micron) selected sample of normal to ultra-luminous infrared galaxies (ULIRGs) [10^9<L_{IR}/L_{sun}<10^12.5] with <z>~0.8 in the Spitzer Extragalactic First Look Survey. We combine Spitzer MIPS 24-micron observations with high-resolution, optical Keck Deimos spectroscopy to derive optical emission-line and infrared star formation rates (SFR_{opt} & SFR_{IR}, respectively). Direct comparison of these SFR diagnostics reveals that our sample exhibits a wide range of extinction (1.0<Av<4.0 mag). This is after removing spectroscopic and IRAC color-selected AGN candidates that account for 12% of the sample. Objects with SFRs of a few solar masses per year have Av values consistent with those of normal spirals (Av~1.0 mag). By contrast, LIRGs at z>1, which make up a large fraction of our sample, have SFR~100 M_{sun}/yr and <Av>~2.5 mag. This translates to a 97% mean attenuation of the [OII] forbidden line doublet, with the most extreme sources having as much as 99.7% of their [OII] line flux extinguished by dust. Based on a SFR_{IR}/SFR_{opt} diagnostic, we derive an IR-luminosity-dependent Av^{IR} function [Av^{IR}=0.75*log(L_{IR}/L_{sun})-6.35 mag] that we use to extinction correct our emission line luminosities. Application of this correction results in a correlation between SFR_{IR} and SFR_{opt} that has a dispersion of 0.2 dex. Investigation of the Av dependence on redshift reveals that for a fixed L_{IR}, there is no significant Av evolution. The mean attenuation of our sample is intermediate between that of local optical/UV- and radio-selected samples and has a marginally stronger L_{IR} dependence.

We present Spitzer 70um and 160um observations of the Spitzer extragalactic First Look Survey (xFLS). The data reduction techniques and the methods for producing co-added mosaics and source catalogs are discussed. Currently, 26% of the 70um sample and 49% of the 160um-selected sources have redshifts. The majority of sources with redshifts are star-forming galaxies at z<0.5, while about 5% have infrared colors consistent with AGN. The observed infrared colors agree with the spectral energy distribution (SEDs) of local galaxies previously determined from IRAS and ISO data. The average 160um/70um color temperature for the dust is Td~= 30+/-5 K, and the average 70um/24um spectral index is alpha~= 2.4+/-0.4. The observed infrared to radio correlation varies with redshift as expected out to z~1 based on the SEDs of local galaxies. The xFLS number counts at 70um and 160um are consistent within uncertainties with the models of galaxy evolution, but there are indications that the current models may require slight modifications. Deeper 70um observations are needed to constrain the models, and redshifts for the faint sources are required to measure the evolution of the infrared luminosity function.

We present Infrared Array Camera (IRAC) data and source catalogs from the Spitzer Space Telescope Extragalactic First Look Survey. The data were taken in four broad bands centered at nominal wavelengths of 3.6, 4.5, 5.8 and 8.0 microns. A set of mosaics and catalogs have been produced which are ~80% complete and ~99% reliable to their chosen flux density limits. The main field survey covers 3.8 deg^2, and has flux density limits of 20muJy, 25muJy, 100muJy and 100muJy at wavelengths of 3.6, 4.5, 5.8 and 8.0 microns,respectively. The deeper ``verification'' survey covers 0.25 deg^2 with limits of 10muJy, 10muJy, 30muJy and 30muJy, respectively. We also include deep data in the ELAIS-N1 field which covers 0.041deg^2 with limits of 4muJy, 3muJy, 10muJy and 10muJy, respectively, but with only two wavelength coverage at a given sky position. The final bandmerged catalogs contain 103193 objects in the main field, 12224 in the verification field and 5239 in ELAIS-N1. Flux densities of high signal-to-noise objects are accurate to about 10%, and the residual systematic error in the absolute flux density scale is ~2-3%. We have successfully extracted sources at source densities as high as 100000 deg^-2 in our deepest 3.6 and 4.5 micron data. The mosaics and source catalogs will be made available through the Spitzer Science Center archive and the Infrared Science Archive.

We present early results from the spectroscopic follow-up of a sample of candidate obscured AGN selected in the mid-infrared from the Spitzer First Look Survey. Our selection allows a direct comparison of the numbers of obscured and unobscured AGN at a given luminosity for the first time, and shows that the ratio of obscured to unobscured AGN at infrared luminosities corresponding to low luminosity quasars is ~1:1 at z~0.5. Most of our optically-faint candidate obscured AGN have the high-ionization, narrow-line spectra expected from type-2 AGN. A composite spectrum shows evidence for Balmer absorption lines, indicating recent star-formation activity in the host galaxies. There is tentative evidence for a decrease in the obscured AGN fraction with increasing AGN luminosity.

The Spitzer Space Telescope was used to study the mid-infrared to far-infrared properties of NGC 300, and to compare dust emission to Halpha to elucidate the heating of the ISM and the star formation cycle at scales < 100 pc. The new data allow us to discern clear differences in the spatial distribution of 8 micron dust emission with respect to 24 micron dust and to HII regions traced by the Halpha light. The 8 micron emission highlights the rims of HII regions, and the 24 micron emission is more strongly peaked in star forming regions than at 8 microns. We confirm the existence and approximate amplitude of interstellar dust emission at 4.5 microns, detected statistically in Infrared Space Observatory (ISO) data, and conclude it arises in star forming regions. When averaging over regions larger than ~ 1 kpc, the ratio of Halpha to Aromatic Feature emission in NGC 300 is consistent with the values observed in disks of spiral galaxies. The mid-to-far-infrared spectral energy distribution of dust emission is generally consistent with pre-Spitzer models.

We present the first results of the angular auto-correlation function of the galaxies detected by the Infrared Array Camera (IRAC) instrument in the First Look Survey (FLS) of the Spitzer Space Telescope. We detect significant signals of galaxy clustering within the survey area. The angular auto-correlation function of the galaxies detected in each of the four IRAC instrument channels is consistent with a power-law form $w(\theta)=A\theta^{1-\gamma}$ out to $\theta = 0.2\arcdeg$, with the slope ranging from $\gamma = 1.5$ to 1.8. We estimate the correlation amplitudes $A$ to be $2.95 \times 10^{-3}$, $2.03 \times 10^{-3}$, $4.53 \times 10^{-3}$, and $2.34 \times 10^{-3}$ at $\theta=1\arcdeg$ for galaxies detected in the IRAC 3.6$\mu$m, 4.5$\mu$m, 5.8$\mu$m, and 8.0$\mu$m instrument channels, respectively. We compare our measurements at 3.6$\mu$m with the previous K-band measurements, and discuss the implications of these results.

IRAC images of M81 show three distinct morphological constituents: a smooth distribution of evolved stars with bulge, disk, and spiral arm components; a clumpy distribution of dust emission tracing the spiral arms; and a pointlike nuclear source. The bulge stellar colors are consistent with M-type giants, and the disk colors are consistent with a slightly younger population. The dust emission generally follows the blue and ultraviolet emission, but there are large areas that have dust emission without ultraviolet and smaller areas with ultraviolet but little dust emission. The former are presumably caused by extinction, and the latter may be due to cavities in the gas and dust created by supernova explosions. The nucleus appears fainter at 8 um than expected from ground-based 10 um observations made four years ago.

We report on submillimeter and infrared observations of 28 radio-selected galaxies in the Spitzer First Look Survey Verification field (FLSV). All of the radio-selected galaxies that show evidence for emission at 850um with SCUBA have Spitzer counterparts at 24um, while only half of the radio-selected galaxies without 850um emission have detectable counterparts at 24um. The data show a wide range of infrared colors (S70/S24 < 5--30, S8/S3.6 < 0.3--4), indicative of a mixture of infrared-warm AGN and cooler starburst dominated sources. The galaxies showing 850um emission have Spitzer flux densities and flux density ratios consistent with the range of values expected for high-redshift (z=1--4) ultraluminous infrared galaxies.

We present Multiband Imaging Photometer for Spitzer (MIPS) observations of M81 at 24, 70, and 160 micron. The grand design nature of M81 is clearly seen, showing two well resolved spiral arms containing numerous bright star forming regions. The MIPS images reveal a significant amount of cold dust associated with the spiral arms. We investigate the variation of the ultraviolet (UV), H-alpha, and infrared (IR) luminosities and star formation rate (SFR) indicators across the face of M81 using the MIPS images and archival UV and H-alpha images. For regions in M81, we find that UV and H-alpha SFRs (uncorrected for dust attenuation) are always lower than the IR SFR. The cause of this behavior is dust attenuation and/or using SFR calibrations appropriate for entire galaxies, not regions in galaxies. The characteristics of the dust attenuation for the regions indicate the dust grains and/or geometry are different from those in starburst galaxies. The behavior of the infrared-radio correlation in M81 is seen to vary from the global average, with variations correlated with the morphology of M81.

We present the Spitzer MIPS 24 micron source counts in the Extragalactic First Look Survey main, verification and ELAIS-N1 fields. Spitzer's increased sensitivity and efficiency in large areal coverage over previous infrared telescopes, coupled with the enhanced sensitivity of the 24 micron band to sources at intermediate redshift, dramatically improve the quality and statistics of number counts in the mid-infrared. The First Look Survey observations cover areas of, respectively, 4.4, 0.26 and 0.015 sq.deg. and reach 3-sigma depths of 0.11, 0.08 and 0.03 mJy. The extragalactic counts derived for each survey agree remarkably well. The counts can be fitted by a super-Euclidean power law of index alpha=-2.9 from 0.2 to 0.9 mJy, with a flattening of the counts at fluxes fainter than 0.2 mJy. Comparison with infrared galaxy evolution models reveals a peak's displacement in the 24 micron counts. This is probably due to the detection of a new population of galaxies with redshift between 1 and 2, previously unseen in the 15 micron deep counts.

We carried out the direct measurement of the fraction of dusty sources in a sample of extremely red galaxies with (R - Ks) > 5.3mag and Ks < 20.2mag, using 24um data from the Spitzer Space Telescope. Combining deep 24um, Ks- and R-band data over an area of 64sq.arcmin in ELAIS N1, we find that 50%+-6% of our ERO sample have measurable 24um flux above the 3sigma flux limit of 40uJy. This flux limit corresponds to a SFR of 12Msun/yr at z~1. The 24um-detected EROs have 24-to-2.2 and 24-to-0.7um flux ratios consistent with infrared luminous, dusty sources at z>1, and an order of magnitude too red to be explained by an infrared quiescent spiral or a pure old stellar population at any redshift. Keck optical spectroscopy of a sample of similarly selected EROs in the FLS field suggests that most of the EROs in ELAIS N1 are probably at z~1. The mean 24um flux (167uJy) of the 24um-detected ERO sample roughly corresponds to the total infrared luminosity of (3-10)x10^(11)Lsun at z=1-1.5, which is similar to that of local LIRGs and ULIGs. The corresponding SFR is ~(50-170)Msun/yr. If the time scale of this starbursting phase is on the order of 10^8yr as inferred for the local LIRGs and ULIGs, the lower limit on the masses of these 24um-detected EROs is (5-20)x10^9Msun. It is plausible that some of the starburst EROs are in the midst of violent transformation to become massive early type galaxies at the epoch of z ~ 1-2.

In this Letter, we present the initial characterization of extragalactic 24um sources in the Spitzer First Look Survey (FLS) by examining their counterparts at 8um and R-band. The color-color diagram of 24-to-8 vs. 24-to-0.7um is populated with 18,734 sources brighter than the 3sigma flux limit of 110uJy, over an area of 3.7sq.degrees. The 24-to-0.7um colors of these sources span almost 4 orders of magnitudes, while the 24-to-8um colors distribute at least over 2 orders of magnitudes. In addition to identifying ~30% of the total sample with infrared quiescent, mostly low redshift galaxies, we also found that: (1) 23% of the 24um sources (~1200/sq.degrees) have very red 24-to-8 and 24-to-0.7 colors and are probably infrared luminous starbursts with L(IR)>3x10^(11)Lsun at z>1. In particular, 13% of the sample (660/sq.degrees) are 24um detected only, with no detectable emission in either 8um or R-band. These sources are the candidates for being ULIRGs at z>2. (2) 2% of the sample (85/sq.degrees) have colors similar to dust reddened AGNs, like Mrk231 at z~0.6-3. (3) We anticipate that some of these sources with extremely red colors may be new types of sources, since they can not be modelled with any familiar type of spectral energy distribution. We find that 17% of the 24um sources have no detectable optical counterparts brighter than R limit of 25.5mag. Optical spectroscopy of these optical extremely faint 24um sources would be very difficult, and mid-infrared spectroscopy from the Spitzer would be critical for understanding their physical nature (Abridged).

Selection of active galactic nuclei (AGN) in the infrared allows the discovery of AGN whose optical emission is extinguished by dust. In this paper, we use the Spitzer Space Telescope First Look Survey (FLS) to assess what fraction of AGN with mid-infrared luminosities comparable to quasars are missed in optical quasar surveys due to dust obscuration. We begin by using the Sloan Digital Sky Survey (SDSS) database to identify 54 quasars within the 4 deg^2 extragalactic FLS. These quasars occupy a distinct region in mid-infrared color space by virtue of their strong, red, continua. This has allowed us to define a mid-infrared color criterion for selecting AGN candidates. About 2000 FLS objects have colors consistent with them being AGN, but most are much fainter in the mid-infrared than the SDSS quasars, which typically have 8 micron flux densities, S(8.0), ~1 mJy. We have investigated the properties of the 43 objects with S(8.0) >= 1 mJy satisfying our AGN color selection. This sample should contain both unobscured quasars, and AGN which are absent from the SDSS survey due to extinction in the optical. After removing 16 known quasars, three probable normal quasars, and eight spurious or confused objects from the initial sample of 43, we are left with 16 objects which are likely to be obscured quasars or luminous Seyfert-2 galaxies. This suggests the numbers of obscured and unobscured AGN are similar in samples selected in the mid-infrared at S(8.0)~1 mJy.

(Abridged) We present R-band images covering more than 11 square degrees of sky obtained with the KPNO 4-m telescope in preparation for the Spitzer Space Telescope First Look Survey. The FLS was designed to characterize the mid-infrared sky at depths 2 orders of magnitude deeper than previous surveys. The extragalactic component is the first cosmological survey done with Spitzer. Two relatively large regions of the sky were observed: the main FLS extra galactic field (17h18m+59d30m) and ELAIS-N1 field (16h10m+54d30m). The overall quality of the images is high. The relative astrometric accuracy is better than 0.1'' and the typical seeing is 1.1''. Images are relatively deep since they reach a median 5-sigma depth limiting magnitude of R=25.5 (Vega). Catalogs have been extracted using SExtractor using thresholds in area and flux for which the number of false detections is below 1% at R=25. Only sources with S/N greater than 3 have been retained in the final catalogs. Comparing the galaxy number counts from our images with those of deeper R-band surveys, we estimate that our observations are 50% complete at R=24.5. These limits in depth are sufficient to identify a substantial fraction of the infrared sources which will be detected by Spitzer.

We present a study of the evolution of the column density distribution, f(N,z), and total neutral hydrogen mass in high-column density quasar absorbers using candidates from a recent high-redshift survey for damped Lyman-alpha (DLA) and Lyman limit system (LLS) absorbers. The observed number of LLS (N(HI)> 1.6 * 10^{17} atom/cm^2) is used to constrain f(N,z) below the classical DLA Wolfe et al. (1986) definition of 2 * 10^{20} atom/cm^2. The joint LLS-DLA analysis shows unambiguously that f(N,z) deviates significantly from a single power law and that a Gamma-law distribution of the form f(N,z)=(f_*/N_*)(N/N_*)^{-Beta} exp(-N/N_*) provides a better description of the observations. These results are used to determine the amount of neutral gas contained in DLAs and in systems with lower column density. Whilst in the redshift range 2 to 3.5, ~90% of the neutral HI mass is in DLAs, we find that at z>3.5 this fraction drops to only 55% and that the remaining 'missing' mass fraction of the neutral gas lies in sub-DLAs with N(HI) 10^{19} - 2 * 10^{20} atom/cm^2. The characteristic column density, N_*, changes from 1.6 * 10^{21} atom/cm^2 at z<3.5 to 2.9 * 10^{20} atom/cm^2 at z>3.5, supporting a picture where at z>3.5, we are directly observing the formation of high column density neutral hydrogen DLA systems from lower column density units. Moreover since current metallicity studies of DLA systems focus on the higher column density systems they may be giving a biased or incomplete view of global galactic chemical evolution at z>3. After correcting the observed mass in HI for the ``missing'' neutral gas the comoving mass density now shows no evidence for a decrease above z=2. (abridged)

We have identified galaxies near two quasars which are at the redshift of damped Lyman-alpha (DLA) systems in the UV spectra of the quasars. Both galaxies are actively forming stars. One galaxy has a luminosity close to the break in the local galaxy luminosity function, L*, the other is significantly fainter than L* and appears to be interacting with a nearby companion. Despite the strong selection effects favoring spectroscopic identification of the most luminous DLA galaxies, many of the spectroscopically-identified DLA galaxies in the literature are sub-L*, suggesting that the majority of the DLA population is probably sub-L*, in contrast to MgII absorbers at similar redshifts whose mean luminosity is close to L*.

We present results on two related topics: 1. A discussion of high redshift candidates (z>4.5), and 2. A study of very small galaxies at intermediate redshifts, both sets being detected in the region of the northern Hubble Deep Field covered by deep NICMOS observations at 1.6 and 1.1 microns. The high redshift candidates are just those with redshift z>4.5 as given in the recent catalog of Thompson, Weymann and Storrie-Lombardi, while the ``small galaxy'' sample is defined to be those objects with isophotal area <= 0.2 squ. arcsec and with photometric redshifts 1<z<4.5. Of the 19 possible high redshift candidates listed in the Thompson et al. catalog, 11 have (nominal) photometric redshifts less than 5.0. Of these, however, only 4 are ``robust'' in the sense of yielding high redshifts when the fluxes are randomly perturbed with errors comparable to the estimated measuring error in each wave band. For the 8 other objects with nominal photometric redshifts greater than 5.0, one (WFPC2 4--473) has a published spectroscopic redshift. Of the remaining 7, 4 are robust in the sense indicated above. Two of these form a close pair (NIC 586 and NIC 107). The redshift of the object having formally the highest redshift, at 6.56 (NIC118 = WFPC2 4--601), is problematic, since F606W and F814W flux are clearly present, and the nature of this object poses a dilemma. (abridged)

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

Quasar absorbers provide a powerful observational tool with which to probe both galaxies and the intergalactic medium up to high redshift. We present a study of the evolution of the column density distribution, f(N,z), and total neutral hydrogen mass in high-column density quasar absorbers using data from a recent high-redshift survey for damped Lyman-alpha (DLA) and Lyman limit system (LLS) absorbers. Whilst in the redshift range 2 to 3.5, ~90% of the neutral HI mass is in DLAs, we find that at z>3.5 this fraction drops to only 55% and that the remaining 'missing' mass fraction of the neutral gas lies in sub-DLAs with N(HI) 10^{19} - 2 * 10^{20} cm^{-2}.

We present the first high-resolution (6.2 to 7.7 km/s FWHM) spectra of the APM BR J0307-4945 quasar at z_em = 4.73 obtained with UVES on the 8.2m VLT Kueyen telescope. We focus our analysis on a damped Ly-alpha (DLA) system at z_abs = 4.466 with a neutral hydrogen column density N(HI) = (4.68 +/- 0.97) 10^{20} cm^{-2}. It is the most distant DLA system known to the present date, observed when the age of the universe was only 1.3 Gyr. It shows complex low- and high-ionization line profiles spanning ~240 and 300 km/s in velocity space respectively. We derive accurate abundances for N, O, Al, Si and Fe, and place a lower limit on C and an upper limit on Ni: [N/H] = -3.07 +/- 0.15, [O/H] = -1.63 +/- 0.19, [Al/H] = -1.79 +/- 0.11, [Si/H] = -1.54 +/- 0.11, [Fe/H] = -1.97 +/- 0.19, [C/H] > -1.63 and [Ni/H] < -2.35. The derived high metallicity, 1/90 solar, shows that this very young absorber (< 1.3 Gyr) has already experienced a significant metal enrichment. The [O/Si] ratio is nearly solar suggesting a limited amount of dust, the relative [Si,O/Fe] abundance ratios show a similar enhancement as observed in the Milky Way stars with comparable metallicities, and the [N/O] ratio is very low. All these results point to an enrichment pattern dominated by Type II supernovae which suggests a Milky Way type evolutionary model.

The cosmological evolution of neutral hydrogen is an efficient way of tracing structure formation with redshift. It indicates the rate of evolution of gas into stars and hence the gas consumption and rate star formation history of the Universe. In measuring HI, quasar absorbers have proven to be an ideal tool and we use observations from a recent survey for high-redshift quasar absorption systems together with data gathered from the literature to measure the cosmological comoving mass density of neutral gas. This paper assumes Omega_M=0.3, Omega_lambda=0.7 and h=0.65.

We present high signal-to-noise, 5 A resolution (FWHM) spectra of 66 z \ga 4 bright quasars obtained with the 4 m Cerro Tololo Inter-American Observatory and 4.2 m William Hershel telescopes. The primary goal of these observations was to undertake a new survey for intervening absorption systems detected in the spectra of background quasars. We look for both Lyman-limit systems (column densities N(HI) > 1.6 * 10^{17} atoms cm-2) and damped Ly-alpha systems (column densities N(HI) > 2 * 10^{20} atoms cm-2). This work resulted in the discovery of 49 Lyman-limit systems, 15 of which are within 3000 km s-1 of the quasar emission and thus might be associated with the quasar itself, 26 new damped Ly-alpha absorption candidates, 15 of which have z>3.5 and numerous metal absorption systems. In addition ten of the quasars presented here exhibit intrinsic broad absorption lines.

We present the spectra, positions, and finding charts for 31 bright (R<19.3) colour-selected quasars covering the redshift range z=3.85-4.78, with 4 having redshifts z>4.5. The majority are in the southern sky (dec<-25 deg). The quasar candidates were selected for their red (B-R > 2.5) colours from UK or POSSII Schmidt Plates scanned at the Automated Plate Measuring facility in Cambridge. Low resolution (~> 10A) spectra were obtained to identify the quasars, primarily at the Las Campanas Observatory. The highest redshift quasar in our survey is at z ~ 4.8 (R = 18.7) and its spectrum shows a damped Lyman-alpha absorption system at z = 4.46. This is currently the highest redshift damped Lyman-alpha absorber detected. Five of these quasars exhibit intrinsic broad absorption line features. Combined with the previously published results from the first part of the APM UKST survey we have now surveyed a total of ~8000 squ deg of sky i.e. 40% of the high galactic latitude (|b|>30 deg) sky, resulting in 59 optically selected quasars in the redshift range 3.85 to 4.78; 49 of which have z>=4.00.

We have carried out a high-resolution imaging search for the galaxy associated with the damped Lyman-alpha (DLA) absorber at z=1.859 toward the z_{em}=2.48 quasar QSO 1244+3443, using the HST and the NICMOS. Images were obtained in the broad filter F160W and the narrow filter F187N with camera 2 on NICMOS with the goal of detecting the rest-frame optical continuum and the H-alpha line emission from the DLA. After PSF subtraction, two weak features are seen at projected separations of 0.16-0.24" from the quasar. Parts of these features may be associated with the DLA absorber, although we cannot completely rule out that they could be artifacts of the point spread function (PSF). If associated with the DLA, the objects would be ~1-2 h_{70}^{-1} kpc in size with integrated flux densities of 2.5 and 3.3 mu Jy in the F160W filter, implying luminosities at lambda_{central}=5600 A in the DLA rest frame of 4.4-5.9 x 10^{9} h_{70}^{-2} L_{solar} at z=1.86, for q0=0.5. However, no significant H-alpha line emission is seen from these objects, suggesting low star formation rates (SFRs). Our 3 sigma upper limit on the SFR in the DLA is 1.3 h_{70}^{-2} M_{solar}/yr for q0 = 0.5 (2.4 h_{70}^{-2} M_{solar} yr^{-1} for q0 = 0.1). This together with our earlier result for LBQS 1210+1731 mark a significant improvement over previous constraints on the star formation rates of DLAs. A combination of low SFR and some dust extinction is likely to be responsible for the lack of H-alpha emission. Alternatively, the objects, may be associated with the quasar host galaxy. In any case, our observations suggest that the DLA is not a large bright proto-disk, but a compact object or a low-surface brightness galaxy. If the two features are PSF artifacts then the constraints on DLA properties are even more severe.

This paper presents the star formation history in the NICMOS Northern Deep HDF. It uses the techniques of photometric redshifts and extinctions to correct for extinction of the ultra-violet flux. It presents a new method for correcting for surface brightness diming. It also predicts the 850 micron fluxes of the objects for comparison with SCUBA measurements

The results of a new spectroscopic survey of 66 $z \simgt 4$ quasars for Damped Lyman-alpha absorption systems are presented. The search led to the discovery of 30 new DLA candidates which are analysed in order to compute the comoving mass density of neutral gas in a non-zero lambda Universe. The possible sources of uncertainty are discussed and the implications of our results for the theories of galaxy formation and evolution are emphasized. A subsequent paper will present details of the calculations summarised here and a more extensive explanation of the consequences of our observations for the understanding of the nature of DLAs.

We present a catalog of extremely red objects discovered using the NICMOS/HST parallel imaging database and ground-based optical follow-up observations. Within an area of 16 square arc-minutes, we detect 15 objects with $\rm R - F160W > 5$ and $\rm F160W < 21.5$. We have also obtained K-band photometry for a subset of the 15 EROs. All of the $\rm R - F160W$ selected EROs imaged at K-band have $\rm R - K > 6$. Our objects have $\rm F110W - F160W$ colors in the range of 1.3 - 2.1, redder than the cluster ellipticals at $z \sim 0.8$ and nearly 1 magnitude redder than the average population selected from the F160W images at the same depth. In addition, among only 22 NICMOS pointings, we detected two groups or clusters in two fields, each contains 3 or more EROs, suggesting that extremely red galaxies may be strongly clustered. At bright magnitudes with $\rm F160W < 19.5$, the ERO surface density is similar to what has been measured by other surveys. At the limit of our sample, F160W = 21.5, our measured surface density is 0.94$\pm 0.24$ arcmin^{-2}. Excluding the two possible groups/clusters and the one apparently stellar object, reduces the surface density to 0.38$\pm 0.15$ arcmin^{-2}.

We have completed spectroscopic observations using LRIS on the Keck 1 telescope of 30 very high redshift quasars, 11 selected for the presence of damped Ly-alpha absorption systems and 19 with redshifts z > 3.5 not previously surveyed for absorption systems. We have surveyed an additional 10 QSOs with the Lick 120'' and the Anglo-Australian Telescope. We have combined these with previous data resulting in a statistical sample of 646 QSOs and 85 damped Ly-alpha absorbers with column densities N(HI) >= 2 x 10^20 atoms/cm^2 covering the redshift range 0.008 <= z <= 4.694. To make the data in our statistical sample more readily available for comparison with scenarios from various cosmological models, we provide tables that includes all 646 QSOs from our new survey and previously published surveys. They list the minimum and maximum redshift defining the redshift path along each line of sight, the QSO emission redshift, and when an absorber is detected, the absorption redshift and measured HI column density. [see the paper for the complete abstract]

We report results of a high-resolution imaging search (in rest frame H-$\alpha$ and optical continuum) for the galaxy associated with the damped Lyman-$\alpha$ (DLA) absorber at $z=1.892$ toward the $z_{em}=2.543$ quasar LBQS 1210+1731, using HST/NICMOS. After PSF subtraction, a feature is seen in both the broad-band and narrow-band images, at a projected separation of 0.25$\arcsec$ from the quasar. If associated with the DLA, the object would be $\approx 2-3$ $h_{70}^{-1}$ kpc in size with a flux of $9.8 \pm 2.4$ $\mu$Jy in the F160W filter, implying a luminosity at $\lambda_{central}=5500$ {\AA} in the rest frame of $1.5 \times 10^{10}$ $h_{70}^{-2}$ L$_{\odot}$ at $z=1.89$, for $q_{0}=0.5$. However, no significant H-$\alpha$ emission is seen, suggesting a low star formation rate (SFR) (3 $\sigma$ upper limit of 4.0 $h_{70}^{-2}$ M$_{\odot}$ yr$^{-1}$), or very high dust obscuration. Alternatively, the object may be associated with the host galaxy of the quasar. H-band images obtained with the NICMOS camera 2 coronagraph show a much fainter structure $\approx 4-5$ $h_{70}^{-1}$ kpc in size and containing four knots of continuum emission, located 0.7$\arcsec$ away from the quasar. We have probed regions far closer to the quasar sight-line than in most previous studies of high-redshift intervening DLAs. The two objects we report mark the closest detected high-redshift DLA candidates yet to any quasar sight line. If the features in our images are associated with the DLA, they suggest faint, compact, somewhat clumpy objects rather than large, well-formed proto-galactic disks or spheroids.

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 have detected the host galaxies of 16 nearby, radio-quiet quasars using images obtained with the Near-Infrared Camera and Multi-Object Spectrometer (NICMOS). We confirm that these luminous quasars tend to live in luminous, early-type host galaxies, and we use the host-galaxy magnitudes to refine the luminosity/host-mass limit inferred from ground-based studies. If quasars obey the relation $M_{black hole}/M_{spheroid}\sim0.006$ found for massive dark objects in nonactive galaxies, then our analysis implies that they radiate at up to $\sim20%$ of the Eddington rate. An analogous analysis for ultraluminous infrared galaxies shows them to accrete at up to similar Eddington fractions, consistent with the hypothesis that some of them are powered by embedded quasars.

This paper contains a catlog of sources observed in a portion of the Hubble Deep Field by the NICMOS instrument. Methods of observation and data analysis are discussed. The catalog has two parts. The first is a listing of all sources with signal to noise ratios greater than or equal to 2.5. The second part contains sources with signal to noise ratios less than 2.5. The second part is intended for researchers who wish to know the limits to which sources may appear at a location.

We present results from our analysis of F160W NICMOS Parallel images. These data cover $\sim$~9~sq. arcminutes and reach 3$\sigma$ depths of H$=$ 24.3 $-$ 25.5 in a $0.6''$ diameter aperture with integration times of 2,000 to 13,000 seconds. We derive the first deep H band galaxy counts. The slope of the counts for H$<$ 20 is 0.31, consistent with various K-band measurements from the Keck telescopes. The measured number counts vs. magnitude relation is reasonably well fitted with no-evolution models with a low $\Omega$ value. The half-light radii of the galaxies declines steeply with apparent magnitude and reaches the NIC3 resoltion limit at H$=$23.5. Deep ground-based VRI imaging of one NICMOS field has revealed an extremely red galaxy with R$-$H $=$ 6 and H of 18.8. Our analyses of the grism data show that we can reach 3$\sigma$ flux limits of of $1\times 10^{-16}$ to $2\times10^{-17}$ ergs/sec/cm$^2$ for integration times of 2,000 to 21,000~seconds. We have detected a total of 33 emission line galaxies. The comoving number density is $\rm \sim 2\times 10^{-4} Mpc^{-3}$. The detected emission lines are probably H$_\alpha$~6563\AA. Thus, the derived star formation rates, without extinction correction, are $10 - 163 M\odot$ per year for galaxies at redshifts between 0.7 and 1.9.

We present Keck LRIS spectroscopy along with NICMOS F110W (~J) and F160W (~H) images of the galaxy HDF4-473.0 (hereafter 4-473) in the Hubble Deep Field, with a detection of an emission line consistent with Ly-alpha at a redshift of z=5.60. Attention to this object as a high redshift galaxy was first drawn by Lanzetta, Yahil and Fernandez-Soto and appeared in their initial list of galaxies with redshifts estimated from the WFPC2 HDF photometry. It was selected by us for spectroscopic observation, along with others in the Hubble Deep Field, on the basis of the NICMOS F110W and F160W and WFPC2 photometry. For H_0 = 65 and q_0 = 0.125, use of simple evolutionary models along with the F814W (~I), F110W, and F160W magnitudes allow us to estimate the star formation rate (~13 M(solar)/yr). The colors suggest a reddening of E(B-V) ~ 0.06. The measured flux in the Ly-alpha line is approximately 1.0*10^(-17) ergs/cm/s and the restframe equivalent width, correcting for the absorption caused by intervening HI, is approximately 90AA. The galaxy is compact and regular, but resolved, with an observed FWHM of ~0.44". Simple evolutionary models can accurately reproduce the colors and these models predict the Ly-alpha flux to within a factor of 2. Using this object as a template shifted to higher redshifts, we calculate the magnitudes through the F814W and two NICMOS passbands for galaxies at redshifts 6 < z < 10.

We present deep galaxy counts and half-light radii from F160W ($\lambda_c=1.6\mu$) images obtained with NICMOS on HST. Nearly 9 arcmin$^2$ have been imaged with camera 3, with $3\sigma$ depths ranging from H = 24.3 to 25.5 in a 0.6$''$ diameter aperture. The slope of the counts fainter than H~$= 20$ is 0.31, and the integrated surface density to H$\leq 24.75$ is $4 \times 10^5$ galaxies per square degree. The half-light radii of the galaxies declines steeply with apparent magnitude. At H~$=24$ we are limited by both the delivered FWHM and the detection threshold of the images.

Though observationally rare, damped Lya absorption systems dominate the mass density of neutral gas in the Universe. Eleven high redshift damped Lya systems covering 2.8<z<4.4 were discovered in 26 QSOs from the APM z>4 QSO Survey, extending these absorption system surveys to the highest redshifts currently possible. Combining our new data set with previous surveys we find that the cosmological mass density in neutral gas, omega_g, does not rise as steeply prior to z~2 as indicated by previous studies. There is evidence in the observed omega_g for a flattening at z~2 and a possible turnover at z~3. When combined with the decline at z>3.5 in number density per unit redshift of damped systems with column densities log N(HI)>21 atoms cm^-2, these results point to an epoch at z>3 prior to which the highest column density damped systems are still forming. We find that over the redshift range 2<z<4 the total mass in neutral gas is marginally comparable with the total visible mass in stars in present day galaxies. However, if one considers the total mass visible in stellar disks alone, ie excluding galactic bulges, the two values are comparable. We are observing a mass of neutral gas comparable to the mass of visible disk stars. Lanzetta, Wolfe & Turnshek (1995) found that omega_g(z~3.5) was twice omega_g(z~2), implying a much larger amount of star formation must have taken place between z=3.5 and z=2 than is indicated by metallicity studies. This created a `cosmic G-dwarf problem'. The more gradual evolution of omega_g we find alleviates this. These results have profound implications for theories of galaxy formation.

Eleven candidate damped Lya absorption systems were identified in 27 spectra of the quasars from the APM z>4 survey covering the redshift range 2.8<z(abs)<4.4 (8 with z(abs)>3.5). High resolution echelle spectra (0.8A FWHM) have been obtained for three quasars, including 2 of the highest redshift objects in the survey. Two damped systems have confirmed HI column densities of N(HI) >= 10^20.3 atoms cm^-2, with a third falling just below this threshold. We have discovered the highest redshift damped Lya absorber known at z=4.383 in QSO BR1202-0725. The APM QSOs provide a substantial increase in the redshift path available for damped surveys for z>3. We combine this high redshift sample with other quasar samples covering the redshift range 0.008 < z < 4.7 to study the redshift evolution and the column density distribution function for absorbers with log N(HI)>=17.2. In the HI column density distribution f(N)=kN^-beta we find evidence for breaks in the power law, flattening for 17.2< log N(HI)<21 and steepening for log N(HI)>21.2. The column density distribution function for the data with log N(HI)>=20.3 is better fit with the form f(N)=(f*/N*)(N/N*)^-beta exp(-N/N*). Significant redshift evolution in the number density per unit redshift is evident in the higher column density systems with an apparent decline in N(z) for z>3.5.

The APM multicolor survey for bright z > 4 objects, covering 2500 deg^2 of sky to m(R)~19, resulted in the discovery of thirty-one quasars with z > 4. High signal-to-noise optical spectrophotometry at 5A resolution has been obtained for the twenty-eight quasars easily accessible from the northern hemisphere. These spectra have been surveyed to create new samples of high redshift Lyman-limit systems, damped Lyman-alpha absorbers, and metal absorption systems (e.g. CIV and MgII). In this paper we present the spectra, together with line lists of the detected absorption systems. The QSOs display a wide variety of emission and absorption line characteristics, with 5 exhibiting broad absorption lines and one with extremely strong emission lines (BR2248-1242). Eleven candidate damped Ly-alpha absorption systems have been identified covering the redshift range 2.8<z<4.4 (8 with z>3.5). An analysis of the measured redshifts of the high ionization emission lines with the low ionization lines shows them to be blueshifted by 430+/-60 km/s. In a previous paper (Storrie-Lombardi et. al. 1994) we discussed the redshift evolution of the Lyman limit systems catalogued here. In subsequent papers we will discuss the properties of the Ly-alpha forest absorbers and the redshift and column density evolution of the damped Ly-alpha absorbers.

We have obtained high signal:to:noise optical spectroscopy at 5\AA\ resolution of 27 quasars from the APM z$>$4 quasar survey. The spectra have been analyzed to create new samples of high redshift Lyman-limit and damped Lyman-$\alpha$ absorbers. These data have been combined with published data sets in a study of the redshift evolution and the column density distribution function for absorbers with $\log$N(HI)$\ge17.5$, over the redshift range 0.01 $<$ z $<$ 5. The main results are: \begin{itemize} \item Lyman limit systems: The data are well fit by a power law $N(z) = N_0(1 + z)^{\gamma}$ for the number density per unit redshift. For the first time intrinsic evolution is detected in the product of the absorption cross-section and comoving spatial number density for an $\Omega = 1$ Universe. We find $\gamma = 1.55$ ($\gamma = 0.5$ for no evolution) and $N_0 = 0.27$ with $>$99.7\% confidence limits for $\gamma$ of 0.82 \& 2.37. \item Damped \lya systems: The APM QSOs provide a substantial increase in the redshift path available for damped surveys for $z>3$. Eleven candidate and three confirmed damped Ly$\alpha$ absorption systems, have been identified in the APM QSO spectra covering the redshift range $2.8\le z \le 4.4$ (11 with $z>3.5$). Combining the APM survey confirmed and candidate damped \lya absorbers with previous surveys, we find evidence for a turnover at z$\sim$3 or a flattening at z$\sim$2 in the cosmological mass density of neutral gas, $\Omega_g$. \end{itemize} The Lyman limit survey results are published in Storrie-Lombardi, et~al., 1994, ApJ, 427, L13. Here we describe the results for the DLA population of absorbers.