This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We present constraints on local primordial non-Gaussianity (PNG), parametrized through $f^{\rm loc}_{\rm NL}$, using the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey Data Release 14 quasar sample. We measure and analyze the anisotropic clustering of the quasars in Fourier space, testing for the scale-dependent bias introduced by primordial non-Gaussianity on large scales. We derive and employ a power spectrum estimator using optimal weights that account for the redshift evolution of the PNG signal. We find constraints of $-51<f^{\rm loc}_{\rm NL}<21$ at 95% confidence level. These are amont the tightest constraints from Large Scale Structure (LSS) data. Our redshift weighting improves the error bar by 15% in comparison to the unweighted case. If quasars have lower response to PNG, the constraint degrades to $-81<f^{\rm loc}_{\rm NL}<26$, with a 40% improvement over the standard approach. We forecast that the full eBOSS dataset could reach $\sigma_{f^{\rm loc}_{\rm NL}}\simeq 5\text{-}8$ using optimal methods and full range of scales.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

Accurate astrometry is crucial for determining orbits of near-Earth-asteroids (NEAs) and therefore better tracking them. This paper reports on a demonstration of 10 milliarcsecond-level astrometric precision on a dozen NEAs using the Pomona College 40 inch telescope, at the JPL's Table Mountain Facility. We used the technique of synthetic tracking, in which many short exposure (1 second) images are acquired and then combined in post-processing to track both target asteroid and reference stars across the field of view. This technique avoids the trailing loss and keeps the jitter effects from atmosphere and telescope pointing common between the asteroid and reference stars, resulting in higher astrometric precision than the 100 mas level astrometry from traditional approach of using long exposure images. Treating our synthetic tracking of near-Earth asteroids as a proxy for observations of future spacecraft while they are downlinking data via their high rate optical communication laser beams, our approach shows precision plane-of-sky measurements can be obtained by the optical ground terminals for navigation. We also discuss how future data releases from the Gaia mission can improve our results.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

We study the host galaxy properties of the tidal disruption object, Swift J164449.3+573451 using long-term optical to near-infrared (NIR) data. First, we decompose the galaxy surface brightness distribution and analyze the morphology of the host galaxy using high resolution \emph{HST} WFC3 images. We conclude that the host galaxy is a bulge-dominant galaxy that is well described by a single S\'{e}rsic model with S\'{e}rsic index $n=3.43\pm0.05$. Adding a disk component, the bulge to total host galaxy flux ratio (B/T) is $0.83\pm0.03$, which still indicates a bulge-dominant galaxy. Second, we estimate multi-band fluxes of the host galaxy through long-term light curves. Our long-term NIR light curves reveal the pure host galaxy fluxes $\sim500$ days after the burst. We fit spectral energy distribution (SED) models to the multi-band fluxes from the optical to NIR of the host galaxy and determine its properties. The stellar mass, the star formation rate, and the age of stellar population are $\log(M_{\star}/M_{\odot}) = 9.14^{+0.13}_{-0.10}$, $0.03^{+0.28}_{-0.03}\, M_{\odot}$/yr, and $0.63^{+0.95}_{-0.43}$ Gyr. Finally, we estimate the mass of the central super massive black hole which is responsible for the tidal disruption event. The black hole mass is estimated to be $10^{6.7\pm0.4}\, M_{\odot}$ from $M_{\mathrm{BH}}$ - $M_{\star,\mathrm{bul}}$ and $M_{\mathrm{BH}}$ - $L_{\mathrm{bul}}$ relations for the $K$ band, although a smaller value of $\sim10^5\, M_{\odot}$ cannot be excluded convincingly if the host galaxy harbors a pseudobulge.

We describe KAPAO, our project to develop and deploy a low-cost, remote-access, natural guide star adaptive optics (AO) system for the Pomona College Table Mountain Observatory (TMO) 1-meter telescope. We use a commercially available 140-actuator BMC MEMS deformable mirror and a version of the Robo-AO control software developed by Caltech and IUCAA. We have structured our development around the rapid building and testing of a prototype system, KAPAO-Alpha, while simultaneously designing our more capable final system, KAPAO-Prime. The main differences between these systems are the prototype's reliance on off-the-shelf optics and a single visible-light science camera versus the final design's improved throughput and capabilities due to the use of custom optics and dual-band, visible and near-infrared imaging. In this paper, we present the instrument design and on-sky closed-loop testing of KAPAO-Alpha as well as our plans for KAPAO-Prime. The primarily undergraduate-education nature of our partner institutions, both public (Sonoma State University) and private (Pomona and Harvey Mudd Colleges), has enabled us to engage physics, astronomy, and engineering undergraduates in all phases of this project. This material is based upon work supported by the National Science Foundation under Grant No. 0960343.

We present the results of a comprehensive Keck/DEIMOS spectroscopic survey of the ultra-faint Milky Way satellite galaxy Segue 1. We have obtained velocity measurements for 98.2% of the stars within 67 pc (10 arcmin, or 2.3 half-light radii) of the center of Segue 1 that have colors and magnitudes consistent with membership, down to a magnitude limit of r=21.7. Based on photometric, kinematic, and metallicity information, we identify 71 stars as probable Segue 1 members, including some as far out as 87 pc. After correcting for the influence of binary stars using repeated velocity measurements, we determine a velocity dispersion of 3.7^{+1.4}_{-1.1} km/s, with a corresponding mass within the half-light radius of 5.8^{+8.2}_{-3.1} x 10^5 Msun. The stellar kinematics of Segue 1 require very high mass-to-light ratios unless the system is far from dynamical equilibrium, even if the period distribution of unresolved binary stars is skewed toward implausibly short periods. With a total luminosity less than that of a single bright red giant and a V-band mass-to-light ratio of 3400 Msun/Lsun, Segue 1 is the darkest galaxy currently known. We critically re-examine recent claims that Segue 1 is a tidally disrupting star cluster and that kinematic samples are contaminated by the Sagittarius stream. The extremely low metallicities ([Fe/H] < -3) of two Segue 1 stars and the large metallicity spread among the members demonstrate conclusively that Segue 1 is a dwarf galaxy, and we find no evidence in favor of tidal effects. We also show that contamination by the Sagittarius stream has been overestimated. Segue 1 has the highest measured dark matter density of any known galaxy and will therefore be a prime testing ground for dark matter physics and galaxy formation on small scales.

We present ground-based optical observations of the September 2009 and January 2010 transits of HD 80606b. Based on 3 partial light curves of the September 2009 event, we derive a midtransit time of T_c [HJD] = 2455099.196 +- 0.026, which is about 1 sigma away from the previously predicted time. We observed the January 2010 event from 9 different locations, with most phases of the transit being observed by at least 3 different teams. We determine a midtransit time of T_c [HJD] = 2455210.6502 +- 0.0064, which is within 1.3 sigma of the time derived from a Spitzer observation of the same event.

We present the far-IR, millimeter, and radio photometry as well as optical and near-IR spectroscopy of a sample of 48 z~1-3 Spitzer-selected ULIRGs with IRS mid-IR spectra. Our goals are to compute their bolometric emission, and to determine both the presence and relative strength of their AGN and starburst components. We find that strong-PAH sources tend to have higher 160um and 1.2mm fluxes than weak-PAH sources. The depth of the 9.7um silicate feature does not affect MAMBO detectability. We fit the far-IR SEDs of our sample and find an average <L_{IR}>~7x10^{12}Lsun for our z>1.5 sources. Spectral decomposition suggests that strong-PAH sources typically have ~20-30% AGN fractions. Weak-PAH sources by contrast tend to have >~70% AGN fractions, with a few sources having comparable contributions of AGN and starbursts. The optical line diagnostics support the presence of AGN in the bulk of the weak-PAH sources. With one exception, our sources are narrow-line sources, show no obvious correspondence between the optical extinction and the silicate feature depth, and, in two cases, show evidence for outflows. Radio AGN are present in both strong-PAH and weak-PAH sources. This is supported by our sample's far-IR-to-radio ratios (q) being consistently below the average value of 2.34 for local star-forming galaxies. We use survival analysis to include the lower-limits given by the radio-undetected sources, arriving at <q>=2.07+/-0.01 for our z>1.5 sample. In total, radio and, where available, optical line diagnostics support the presence of AGN in 57% of the z>1.5 sources, independent of IR-based diagnostics. For higher-z sources, the AGN luminosities alone are estimated to be >10^{12}Lsun, which, supported by the [OIII] luminosities, implies that the bulk of our sources host obscured quasars.

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 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 mid-IR spectral decomposition of a sample of 48 Spitzer-selected ULIRGs spanning z~1-3 and likely L_IR~10^12-10^13Lsun. Our study aims at quantifying the star-formation and AGN processes in these sources which recent results suggest have evolved strongly between the observed epoch and today. To do this, we study the mid-IR contribution of PAH emission, continuum, and extinction. About 3/4 of our sample are continuum- (i.e. AGN) dominated sources, but ~60% of these show PAH emission, suggesting the presence of star-formation activity. These sources have redder mid-IR colors than typical optically-selected quasars. About 25% of our sample have strong PAH emission, but none are likely to be pure starbursts as reflected in their relatively high 5um hot dust continua. However, their steep 30um-to-14um slopes suggest that star-formation might dominate the total infrared luminosity. Six of our z~2 sources have EW6.2>~0.3um and L_14um>~10^12Lsun (implying L_IR>~10^13Lsun). At these luminosities, such high EW6.2 ULIRGs do not exist in the local Universe. We find a median optical depth at 9.7um of <tau_9.7>=1.4. This is consistent with local IRAS-selected ULIRGs, but differs from early results on SCUBA-selected z~2 ULIRGs. Similar to local ULIRGs about 25% of our sample show extreme obscuration (tau_9.7>~3) suggesting buried nuclei. In general, we find that our sources are similar to local ULIRGs, but are an order of magnitude more luminous. It is not clear whether our z~2 ULIRGs are simply scaled-up versions of local ULIRGs, or subject to fundamentally different physical processes.

We present the rest-frame 8 micron luminosity function (LF) at redshifts z=1 and ~2, computed from Spitzer 24 micron-selected galaxies in the GOODS fields over an area of 291 sq. arcmin. Using classification criteria based on X-ray data and IRAC colours, we identify the AGN in our sample. The rest-frame 8 micron LF for star-forming galaxies at redshifts z=1 and ~2 have the same shape as at z~0, but with a strong positive luminosity evolution. The number density of star-forming galaxies with log_{10}(nu L_nu(8 micron))>11 increases by a factor >250 from redshift z~0 to 1, and is basically the same at z=1 and ~2. The resulting rest-frame 8 micron luminosity densities associated with star formation at z=1 and ~2 are more than four and two times larger than at z~0, respectively. We also compute the total rest-frame 8 micron LF for star-forming galaxies and AGN at z~2 and show that AGN dominate its bright end, which is well-described by a power-law. Using a new calibration based on Spitzer star-forming galaxies at 0<z<0.6 and validated at higher redshifts through stacking analysis, we compute the bolometric infrared (IR) LF for star-forming galaxies at z=1 and ~2. We find that the respective bolometric IR luminosity densities are (1.2+/-0.2) x 10^9 and (6.6^{+1.2}_{-1.0}) x 10^8 L_sun Mpc^{-3}, in agreement with previous studies within the error bars. At z~2, around 90% of the IR luminosity density associated with star formation is produced by luminous and ultraluminous IR galaxies (LIRG and ULIRG), with the two populations contributing in roughly similar amounts. Finally, we discuss the consistency of our findings with other existing observational results on galaxy evolution.

By combining the 0.12 square degree F814W Hubble Space Telescope (HST) and Spitzer MIPS 24 micron imaging in the First Look Survey (FLS), we investigate the properties of interacting and merging Mid-Infrared bright and faint sources at 0.2< z <1.3. We find a marginally significant increase in the pair fraction for MIPS 24 micron detected, optically selected close pairs, pair fraction=0.25+/-0.10 at z~1, in contrast to 0.11+/-0.08 at z~0.4, while galaxies below our 24 micron MIPS detection limit show a pair fraction consistent with zero at all redshifts. Additionally, 24 micron detected galaxies with fluxes >0.1mJy are on average five times more likely to be in a close galaxy pair between 0.2< z <1.3 than galaxies below this flux limit. Using the 24 micron flux to derive the total Far-IR luminosity we find that paired galaxies (early stage mergers) are responsible for 27% +/-9% of the IR luminosity density resulting from star formation at z~1 while morphologically classified (late stage) mergers make up 34%+/-11%. This implies that 61%+/-14% of the infrared luminosity density and in turn ~40% of the star formation rate density at z~1 can be attributed to galaxies at some stage of a major merger or interaction. We argue that, close pairs/mergers in a LIRG/ULIRG phase become increasingly important contributers to the IR luminosity and star formation rate density of the Universe at z>0.7.

We present the mid-infrared (MIR) spectra obtained with the Spitzer InfraRed Spectrograph (IRS) for a sample of 52 sources, selected as infrared luminous, z>1 candidates in the Extragalactic First Look Survey (XFLS). The sample selection criteria are f(24um) > 0.9mJy, nu fnu(24um)/nu fnu(8um) > 3.16 and nu fnu(24um)/nu fnu(0.7um) > 10. Of the 52 spectra, 47 (90%) produced measurable redshifts based solely on the mid-IR spectral features, with the majority (35/47=74%) at 1.5<z<3.2. Keck spectroscopy of a sub-sample (17/47) agrees with the mid-IR redshift measurements. The observed spectra fall crudely into three categories -- (1) 33% (17/52) have strong PAH emission, and are probably powered by star formation with total IR luminosity roughly a factor of 5 higher than the local starburst ULIRGs. (2) 33% (17/52) have only deep silicate absorption at 9.8um, indicative of deeply embedded dusty systems. The energetic nature of the heating sources in these systems can not be determined by these data alone. (3) The remainder 34% are mid-IR continuum dominated systems with either weak PAH emission and/or silicate absorption. This third of the sample are probably AGNs. From the silicate absorption feature, we estimate that roughly two-thirds of the sample have optical depth tau(9.8um) > 1. Their $L_{1600\AA}$ and $L_{\rm IR}$ suggest that our sample is among the most luminous and most dust enshrouded systems of its epoch . Our study has revealed a significant population of dust enshrouded galaxies at z~2, whose enormous energy output, comparable to that of quasars, is generated by AGN as well as starburst. This IR luminous population has very little overlap with sub-mm and UV-selected populations (Abridged).

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 report the initial results from a Spitzer GO-1 program to obtain low resolution, mid-infrared spectra of infrared luminous galaxies at z~1-2. This paper presents the spectra of eight sources observed with the Spitzer InfraRed Spectrograph (IRS). Of the eight spectra, six have mid-IR spectral features, either emission from Polycyclic Aromatic Hydrocarbon (PAH) or silicate absorption. Based on these mid-IR features, the inferred six redshifts are in the range of 1.8-2.6. The remaining two spectra detect only strong continua, thus do not yield redshift information. Strong, multiple PAH emission features are detected in two sources, and weak PAH emission in another two. These data provide direct evidence that PAH molecules are present and directly observable in ULIRGs at z~2. The six sources with measured redshifts are dusty, infrared luminous galaxies at z~2 with estimated $L_{bol} \sim 10^{13}L_\odot$. Of the eight sources, two appear starburst dominated; two with only power law continua are probably type I QSOs; and the remaining four are likely composite systems containing a buried AGN and a starburst component. Since half of our sample are optically faint sources with R>25.5mag (Vega), our results demonstrate the potential of using mid-infrared spectroscopy, especially the Aromatic and silicate features produced by dust grains to directly probe optically faint and infrared luminous populations at high redshift.

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.

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.

Using the SPITZER SPACE TELESCOPE and the VLA, we present the first DIRECT evidence that the well-known far-infrared/radio correlation is valid to cosmologically significant redshift. We also confirm, with improved statistics compared with previous surveys, a similar result for the Mid-IR/radio correlation. We explore the dependence of monochromatic q_{24} and q_{70} on z. The results were obtained by matching Spitzer sources at 24 and 70 microns with VLA 1.4 GHz micro-Jy radio sources obtained for the SPITZER FLS. Spectroscopic redshifts have been obtained for over 500 matched IR/radio sources using observations at WIYN, Keck and archival SDSS data extending out to z > 2. We find that q_{24} shows significantly more dispersion than q_{70}. By comparing the observed fluxes at 70, 24 and 4.5 microns with a library of SED templates, we find that the larger dispersion in q_{24} is predictable in terms of systematic variations in SED shape throughout the population. Although the models are not able to encompass the full range of observed behavior (both the presence of either extremely flat or extremely steep IR SEDs), the fitting parameters were used to `k-correct' the higher-z galaxies which resulted in a reduced scatter in q. For comparison, we also corrected these data using the SED for M82. The results for 24 and 70 microns provide strong consistent evidence for the universality of the mid-IR/radio and far-IR/radio correlations out to redshifts of at least z = 1.

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

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.

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.

Recent observations of surface brightness distributions of both Milky Way and M31 satellite galaxies have revealed many instances of sudden changes or breaks in the slope of the surface brightness profiles (at some break radius, r_break). These breaks are often accompanied by increasingly elliptical isophotes and sometimes by isophote twisting. We investigate the hypothesis of a tidal origin for these features by applying the same ellipse fitting techniques that are used on observed galaxies to numerical simulations of the destruction of satellites, represented by spherical, single-component systems. We examine how observed quantities such as r_break, ellipticity and position angle of the fitted ellipses and amplitude of extra-break population vary with the satellite's orbital eccentricity and phase, and our viewpoint relative to the orbit. We also look at orbit and viewpoint dependence of the rate of change of the latter three quantities with radius. We find that there are trends with orbital phase and eccentricity in all observed quantities, many of which are preserved through a wide variety of viewing angles. In particular, a generic feature of all simulations is a depletion zone just interior to an excess zone, regions in which the surface brightness is lower and higher, respectively, than the initial profile. A clear interpretation of any individual image, however, is likely to be hampered by the dependence of the observable features on these multiple parameters. For example, breaks can be excited by several physical processes and can occur well within the bound satellite population. Nevertheless, we do find we can place loose constraints on the tidal radius, mass loss rate, orbital type and phase of the satellite, and nature of breaks using photometric data alone.

We investigate the interaction history of the M31 sub-group of galaxies by comparing surface photometry of its two nearest satellites, M32 and NGC 205, with N-body simulations of satellite destruction. The recent discovery of a giant stream in M31's outer halo, apparently pointed towards M32 and NGC 205, makes such an investigation particularly relevant. The observational component of this study is based on 1.7 x 5.0 degree B- and I- band CCD mosaic images centered on M31 and covering both satellites. The simulation component of this project is based on the analysis of single-component, spherical model satellites that are tidally disrupted through interactions with their parent galaxy. Generic features of the simulations include an excess in the surface brightness profile at large radii, a depletion zone at intermediate radii, and isophotal elongation and twists that are coincident with breaks in the brightness profile. The satellites M32 and NGC 205 display most of these features consistently across the B and I bands, strongly suggestive of tidal interaction and probable stripping by M31. We discuss what these observed features can tell us about the satellites' orbital parameters and histories. Specifically, M32 is found to be on a highly eccentric orbit and away from pericenter. Investigating M32's unusual combination of high surface brightness and low luminosity (the hallmark of compact ellipticals), we conclude on empirical grounds that M32 is not likely to be the residual core of a tidally-stripped normal elliptical galaxy as has been suggested, but rather that its precursor was intrinsically compact.

Surface photometry of the M31 satellites M32 and NGC 205 is compared to numerical simulations of satellite destruction to constrain orbital parameters and the interaction history of the M31 subgroup. Our analysis reveals the following preliminary results: (1) Generic features of tidal disruption in the simulations include an extended ``extra-tidal'' excess region and an inner depletion zone, both of which are observed in M32 and NGC 205; (2) M32 is likely to be on a highly eccentric orbit well away from pericenter; (3) Surface brightness and luminosity evolution estimates for M32, the prototypical compact elliptical galaxy, imply that it is not simply the residual core of a tidally-stripped normal elliptical galaxy, but was instead formed in a truncated state.

We report results from a large Hubble Space Telescope project to observe a significant (~34,000) ensemble of main sequence stars in the globular cluster 47 Tucanae with a goal of defining the frequency of inner-orbit, gas-giant planets. Simulations based on the characteristics of the 8.3 days of time-series data in the F555W and F814W WFPC2 filters show that ~17 planets should be detected by photometric transit signals if the frequency of hot Jupiters found in the solar neighborhood is assumed to hold for 47 Tuc. The experiment provided high-quality data sufficient to detect planets. A full analysis of these WFPC2 data reveals ~75 variables, but NO light curves resulted for which a convincing interpretation as a planet could be made. The planet frequency in 47 Tuc is at least an order of magnitude below that for the solar neighborhood. The cause of the absence of close-in planets in 47 Tuc is not yet known; presumably the low metallicity and/or crowding of 47 Tuc interfered with planet formation, with orbital evolution to close-in positions, or with planet survival.