Spin off events and impacts can eject boulders from an asteroid surface and rubble pile asteroids can accumulate from debris following a collision between large asteroids. These processes produce a population of gravitational bound objects in orbit that can impact an asteroid surface at low velocity and with a distribution of impact angles. We present laboratory experiments of low velocity spherical projectiles into a fine granular medium, sand. We delineate velocity and impact angles giving ricochets, those giving projectiles that roll-out from the impact crater and those that stop within their impact crater. With high speed camera images and fluorescent markers on the projectiles we track spin and projectile trajectories during impact. We find that the projectile only reaches a rolling without slipping condition well after the marble has reached peak penetration depth. The required friction coefficient during the penetration phase of impact is 4-5 times lower than that of the sand suggesting that the sand is fluidized near the projectile surface during penetration. We find that the critical grazing impact critical angle dividing ricochets from roll-outs, increases with increasing impact velocity. The critical angles for ricochet and for roll-out as a function of velocity can be matched by an empirical model during the rebound phase that balances a lift force against gravity. We estimate constraints on projectile radius, velocity and impact angle that would allow projectiles on asteroids to ricochet or roll away from impact, finally coming to rest distant from their initial impact sites.

Thermal infrared measurements of near-Earth objects provide critical data for constraining their physical properties such as size. The NEOWISE mission has been conducting an all-sky infrared survey to gather such data and improve our understanding of this population. While automated routines are employed to identify the majority of moving objects detected by NEOWISE, a subset of objects will have dynamical properties that fall outside the window detectable to these routines. Using the population of known near-Earth objects, we have conducted a manual search for detections of these objects that were previously unreported. We report 303 new epochs of observations for 299 unique near-Earth objects of which 239 have no previous physical property characterization from the NEOWISE Reactivation mission. As these objects are drawn from a list with inherent optical selection biases, the distribution of measured albedos is skewed to higher values than is seen for the diameter-selected population detected by the automated routines. These results demonstrate the importance and benefit of periodic searches of the archival NEOWISE data.

CatWISE is a program to catalog sources selected from combined ${\it WISE}$ and ${\it NEOWISE}$ all-sky survey data at 3.4 and 4.6 $\mu$m (W1 and W2). The CatWISE Preliminary Catalog consists of 900,849,014 sources measured in data collected from 2010 to 2016. This dataset represents four times as many exposures and spans over ten times as large a time baseline as that used for the AllWISE Catalog. CatWISE adapts AllWISE software to measure the sources in coadded images created from six-month subsets of these data, each representing one coverage of the inertial sky, or epoch. The catalog includes the measured motion of sources in 8 epochs over the 6.5 year span of the data. From comparison to ${\it Spitzer}$, the SNR=5 limits in magnitudes in the Vega system are W1=17.67 and W2=16.47, compared to W1=16.96 and W2=16.02 for AllWISE. From comparison to ${\it Gaia}$, CatWISE positions have typical accuracies of 50 mas for stars at W1=10 mag and 275 mas for stars at W1=15.5 mag. Proper motions have typical accuracies of 10 mas yr$^{-1}$ and 30 mas yr$^{-1}$ for stars with these brightnesses, an order of magnitude better than from AllWISE. The catalog is available in the WISE/NEOWISE Enhanced and Contributed Products area of the NASA/IPAC Infrared Science Archive.

The Near-Earth Object Wide-field Infrared Survey Explorer (NEOWISE) spacecraft has been conducting a two-band thermal infrared survey to detect and characterize asteroids and comets since its reactivation in Dec 2013. Using the observations collected during the fourth and fifth years of the survey, our automated pipeline detected candidate moving objects which were verified and reported to the Minor Planet Center. Using these detections, we perform thermal modeling of each object from the near-Earth object and Main Belt asteroid populations to constrain their sizes. We present thermal model fits of asteroid diameters for 189 NEOs and 5831 MBAs detected during the fourth year of the survey, and 185 NEOs and 5776 MBAs from the fifth year. To date, the NEOWISE Reactivation survey has provided thermal model characterization for 957 unique NEOs. Including all phases of the original WISE survey brings the total to 1473 unique NEOs that have been characterized between 2010 and the present.

Mass-spring model simulations are used to investigate past spin states of a viscoelastic Phobos and Deimos. From an initially tidally locked state, we find crossing of a spin-orbit resonance with Mars or a mean motion resonance with each other does not excite tumbling in Phobos or Deimos. However, once tumbling our simulations show that these moons can remain so for an extended period and during this time their orbital eccentricity can be substantially reduced. We attribute the tendency for simulations of an initially tumbling viscoelastic body to drop into spin-synchronous state at very low eccentricity to the insensitivity of the tumbling chaotic zone volume to eccentricity. After a tumbling body enters the spin synchronous resonance, it can exhibit long lived non-principal axis rotation and this too can prolong the period of time with enhanced tidally generated energy dissipation. The low orbital eccentricities of Phobos and Deimos could in part be due to spin excitation by nearly catastrophic impacts rather than tidal evolution following orbital resonance excitation.

The Origins Space Telescope (Origins) traces our cosmic history, from the formation of the first galaxies and the rise of metals to the development of habitable worlds and present-day life. Origins does this through exquisite sensitivity to infrared radiation from ions, atoms, molecules, dust, water vapor and ice, and observations of extra-solar planetary atmospheres, protoplanetary disks, and large-area extragalactic fields. Origins operates in the wavelength range 2.8 to 588 microns and is 1000 times more sensitive than its predecessors due to its large, cold (4.5 K) telescope and advanced instruments. Origins was one of four large missions studied by the community with support from NASA and industry in preparation for the 2020 Decadal Survey in Astrophysics. This is the final study report.

We present follow-up $Spitzer$ observations at 3.6$\mu$m (ch1) and 4.5$\mu$m (ch2) of CWISEP J144606.62$-$231717.8, one of the coldest known brown dwarfs in the solar neighborhood. This object was found by mining the Wide-field Infrared Survey Explorer ($WISE$) and $NEOWISE$ data via the CatWISE Preliminary Catalog by Meisner et al. (2019b), where an initial $Spitzer$ color of ch1$-$ch2 = 3.71$\pm$0.44 mag was reported, implying it could be one of the reddest, and hence coldest, known brown dwarfs. Additional $Spitzer$ data presented here allows us to revise its color to ch1$-$ch2 = 2.986$\pm$0.048 mag, which makes CWISEP J144606.62$-$231717.8 the 5th reddest brown dwarf ever observed. A preliminary trigonometric parallax measurement, based on a combination of $WISE$ and $Spitzer$ astrometry, places this object at a distance of 10.1$^{+1.7}_{-1.3}$ pc. Based on our improved $Spitzer$ color and preliminary parallax, CWISEP J144606.62$-$231717.8 has a $T_{\rm eff}$ in the 310$-$360 K range. Assuming an age of 0.5$-$13 Gyr, this corresponds to a mass between 2 and 20 $M_{\rm Jup}$.

We consider how an impact generated seismic pulse affects the surface of an asteroid distant from the impact site. With laboratory experiments on dry polydisperse gravel mixtures, we track the trajectories of particles ejected from the surface by a single strong upward propagating pressure pulse. High speed video images show that ejecta trajectories are independent of particle size, and collisions primarily take place upon landing. When they land particles are ballistically sorted, as proposed by Shinbrot et al. (2017), leaving larger particles on the surface and smaller particles more widely dispersed. A single strong pulse can leave previously buried boulders stranded on the surface. Boulder stranding due to an impact excited seismic pulse is an additional mechanism that could leave large boulders present on the surface of rubble asteroids such as 162173 Ryugu, 101955 Bennu and 25143 Itokawa.

We present VLT/XSHOOTER rest-frame UV-optical spectra of 10 Hot Dust-Obscured Galaxies (Hot DOGs) at $z\sim2$ to investigate AGN diagnostics and to assess the presence and effect of ionized gas outflows. Most Hot DOGs in this sample are narrow-line dominated AGN (type 1.8 or higher), and have higher Balmer decrements than typical type 2 quasars. Almost all (8/9) sources show evidence for ionized gas outflows in the form of broad and blueshifted [O III] profiles, and some sources have such profiles in H$\alpha$ (5/7) or [O II] (3/6). Combined with the literature, these results support additional sources of obscuration beyond the simple torus invoked by AGN unification models. Outflow rates derived from the broad [O III] line ($\rm \gtrsim10^{3}\,M_{\odot}\,yr^{-1}$) are greater than the black hole accretion and star formation rates, with feedback efficiencies ($\sim0.1-1\%$) consistent with negative feedback to the host galaxy's star formation in merger-driven quasar activity scenarios. We find the broad emission lines in luminous, obscured quasars are often better explained by outflows within the narrow line region, and caution that black hole mass estimates for such sources in the literature may have substantial uncertainty. Regardless, we find lower bounds on the Eddington ratio for Hot DOGs near unity.

We present Spitzer 3.6$\mu$m and 4.5$\mu$m follow-up of 170 candidate extremely cool brown dwarfs newly discovered via the combination of WISE and NEOWISE imaging at 3$-$5$\mu$m. CatWISE, a joint analysis of archival WISE and NEOWISE data, has improved upon the motion measurements of AllWISE by leveraging a $>$10$\times$ time baseline enhancement, from 0.5 years (AllWISE) to 6.5 years (CatWISE). As a result, CatWISE motion selection has yielded a large sample of previously unrecognized brown dwarf candidates, many of which have archival detections exclusively in the WISE 4.6$\mu$m (W2) channel, suggesting that they could be both exceptionally cold and nearby. Where these objects go undetected in WISE W1 (3.4$\mu$m), Spitzer can provide critically informative detections at 3.6$\mu$m. Of our motion-confirmed discoveries, seventeen have a best-fit Spitzer [3.6]$-$[4.5] color most consistent with spectral type Y. CWISEP J144606.62$-$231717.8 ($\mu \approx 1.3''$/yr) is likely the reddest, and therefore potentially coldest, member of our sample with a very uncertain [3.6]$-$[4.5] color of 3.71 $\pm$ 0.44 magnitudes. We also highlight our highest proper motion discovery, WISEA J153429.75$-$104303.3, with $\mu \approx 2.7''$/yr. Given that the prior list of confirmed and presumed Y dwarfs consists of just 27 objects, the Spitzer follow-up presented in this work has substantially expanded the sample of identified Y dwarfs. Our new discoveries thus represent significant progress toward understanding the bottom of the substellar mass function, investigating the diversity of the Y dwarf population, and selecting optimal brown dwarf targets for JWST spectroscopy.

The study of the extragalactic background light (EBL) in the optical and near infrared has received a lot of attention in the last decade, especially near a wavelength of $\lambda\approx 3.4\operatorname{\mu m}$, with remaining tension among different techniques for estimating the background. In this paper we present a measurement of the contribution of galaxies to the EBL at $3.4\operatorname{\mu m}$ that is based on the measurement of the luminosity function (LF) in Lake et al. (2018) and the mean spectral energy distribution of galaxies in Lake & Wright (2016). The mean and standard deviation of our most reliable Bayesian posterior chain gives a $3.4\operatorname{\mu m}$ background of $I_\nu = 9.0\pm0.5 \operatorname{kJy} \operatorname{sr}^{-1}$ ($\nu I_\nu = 8.0\pm0.4 \operatorname{nW} \operatorname{m}^{-2} \operatorname{sr}^{-1} e\operatorname{-fold}^{-1}$), with systematic uncertainties unlikely to be greater than $2\operatorname{kJy} \operatorname{sr}^{-1}$. This result is higher than most previous efforts to measure the contribution of galaxies to the $3.4\operatorname{\mu m}$ EBL, but is consistent with the upper limits placed by blazars and the most recent direct measurements of the total $3.4\operatorname{\mu m}$ EBL.

We report the detection of a large amplitude MIR outburst in the white dwarf (WD) 0145+234 in the NEOWISE Survey data. The source had a stable MIR flux before 2018, and was brightened by about 1.0 magnitude in the W1 and W2 bands within half a year and has been continuously brightening since then. No significant variations are found in the optical photometry data during the same period. This suggests that this MIR outburst is caused by recent replenishing or redistribution of dust, rather than intrinsic variations of the WD. SED modeling of 0145+234 suggests that there was already a dust disk around the WD in the quiescent state, and both of the temperature and surface area of the disk evolved rapidly since the outburst. The dust temperature was about 1770K in the initial rising phase, close to the sublimation temperature of silicate grains, and gradually cooled down to around 1150K, while the surface area increased by a factor of about 6 during the same period. The inferred closest distance of dust to the WD is within the tidal disruption radius of a gravitationally bounded asteroid. We estimated the dust mass to be between $3\times10^{15}$ to $3\times10^{17}$ $\rho/(1 \mathrm g~cm^{-3})$ kg for silicate grains of a power-law size distribution with a high cutoff size from 0.1 to 1000 microns. We interpret this as a possible tidal breakup of an exo-asteroid by the WD. Further follow-up observations of this rare event may provide insights on the origin of dust disk and metal pollution in some white dwarfs.

Main belt asteroid (6478) Gault has been dynamically linked with two overlapping asteroid families: Phocaea, dominated by S-type asteroids, and Tamara, dominated by low-albedo C-types. This object has recently become an interesting case for study, after images obtained in late 2018 revealed that it was active and displaying a comet-like tail. Previous authors have proposed that the most likely scenarios to explain the observed activity on Gault were rotational excitation or merger of near-contact binaries. Here we use new photometric and spectroscopic data of Gault to determine its physical and compositional properties. Lightcurves derived from the photometric data showed little variation over three nights of observations, which prevented us from determining the rotation period of the asteroid. Using WISE observations of Gault and the near-Earth Asteroid Thermal Model (NEATM) we determined that this asteroid has a diameter $<$6 km. NIR spectroscopic data obtained with the Infrared Telescope Facility (IRTF) showed a spectrum similar to that of S-complex asteroids, and a surface composition consistent with H chondrite meteorites. These results favor a compositional affinity between Gault and asteroid (25) Phocaea, and rules out a compositional link with the Tamara family. From the spectroscopic data we found no evidence of fresh material that could have been exposed during the outburst episodes.

Thermophysical models allow for improved constraints on the physical and thermal surface properties of asteroids beyond what can be inferred from more simple thermal modeling, provided a sufficient number of observations is available. We present thermophysical modeling results of observations from the NEOWISE mission for two near-Earth asteroids which are the targets of the DESTINY+ flyby mission: (3200) Phaethon and (155140) 2005 UD. Our model assumes a rotating, cratered, spherical surface, and employs a Monte Carlo Markov Chain to explore the multi-dimensional parameter space of the fit. We find an effective spherical diameter for Phaethon of $4.6^{+0.2}_{-0.3}~$km, a geometric albedo of $p_V=0.16\pm0.02$, and a thermal inertia $\Gamma=880$ $^{+580}_{-330}$, using five epochs of NEOWISE observations. The best model fit for (155140) 2005 UD was less well constrained due to only having two NEOWISE observation epochs, giving a diameter of $1.2\pm0.4~$km and a geometric albedo of $p_V=0.14\pm0.09$.

We present the discovery of an extremely cold, nearby brown dwarf in the solar neighborhood, found in the CatWISE catalog (Eisenhardt et al., in prep.). Photometric follow-up with Spitzer reveals that the object, CWISEP J193518.59-154620.3, has ch1$-$ch2 = 3.24$\,\pm\,$0.31 mag, making it one of the reddest brown dwarfs known. Using the Spitzer photometry and the polynomial relations from Kirkpatrick et al. (2019) we estimate an effective temperature in the $\sim$270--360 K range, and a distance estimate in the 5.6$-$10.9 pc range. We combined the WISE, NEOWISE, and Spitzer data to measure a proper motion of $\mu_\alpha \cos \delta = 337\pm69$ mas yr$^{-1}$, $\mu_\delta = -50\pm97$ mas yr$^{-1}$, which implies a relatively low tangential velocity in the range 7$-$22 km s$^{-1}$.

The gaseous mass of protoplanetary disks is a fundamental quantity in planet formation. The presence of gas is necessary to assemble planetesimals, it determines timescales of giant planet birth, and it is an unknown factor for a wide range of properties of planet formation, from chemical abundances (X/H) to the mass efficiency of planet formation. The gas mass obtained from traditional tracers, such as dust thermal continuum and CO isotopologues, are now known to have significant (1 - 2 orders of magnitude) discrepancies. Emission from the isotopologue of H2, hydrogen deuteride (HD), offers an alternative measurement of the disk gas mass. Of all of the regions of the spectrum, the far-infrared stands out in that orders of magnitude gains in sensitivity can be gleaned by cooling a large aperture telescope to 8 K. Such a facility can open up a vast new area of the spectrum to exploration. One of the primary benefits of this far-infrared revolution would be the ability to survey hundreds of planet-forming disks in HD emission to derive their gaseous masses. For the first time, we will have statistics on the gas mass as a function of evolution, tracing birth to dispersal as a function of stellar spectral type. These measurements have broad implications for our understanding of the time scale during which gas is available to form giant planets, the dynamical evolution of the seeds of terrestrial worlds, and the resulting chemical composition of pre-planetary embryos carrying the elements needed for life. Measurements of the ground-state line of HD requires a space-based observatory operating in the far-infrared at 112 microns.

We present the Massive and Distant Clusters of WISE Survey (MaDCoWS), a search for galaxy clusters at 0.7<z<1.5 based upon data from the Wide-field Infrared Survey Explorer (WISE) mission. MaDCoWS is the first cluster survey capable of discovering massive clusters at these redshifts over the full extragalactic sky. The search is divided into two regions -- the region of the extragalactic sky covered by Pan-STARRS (Dec>-30 degrees) and the remainder of the southern extragalactic sky at Dec<-30 degrees for which shallower optical data from SuperCOSMOS Sky Survey are available. In this paper we describe the search algorithm, characterize the sample, and present the first MaDCoWS data release -- catalogs of the 2433 highest amplitude detections in the WISE--Pan-STARRS region and the 250 highest amplitude detections in the WISE--SuperCOSMOS region. A total of 1723 of the detections from the WISE--Pan-STARRS sample have also been observed with the Spitzer Space Telescope, providing photometric redshifts and richnesses, and an additional 64 detections within the WISE--SuperCOSMOS region also have photometric redshifts and richnesses. Spectroscopic redshifts for 38 MaDCoWS clusters with IRAC photometry demonstrate that the photometric redshifts have an uncertainty of $\sigma_z/(1+z)\sim0.036$. Combining the richness measurements with Sunyaev-Zel'dovich observations of MaDCoWS clusters, we also present a preliminary mass-richness relation that can be used to infer the approximate mass distribution of the full sample. The estimated median mass for the WISE--Pan-STARRS catalog is $M_{500}=1.6^{+0.7}_{-0.8}\times10^{14} \mathrm{M}_\odot$, with the Sunyaev-Zel'dovich data confirming that we detect clusters with masses up to $M_{500}\sim5\times10^{14} \mathrm{M}_\odot$ $(M_{200}\sim10^{15} \mathrm{M}_\odot)$.

We present preliminary trigonometric parallaxes of 184 late-T and Y dwarfs using observations from Spitzer (143), USNO (18), NTT (14), and UKIRT (9). To complete the 20-pc census of $\ge$T6 dwarfs, we combine these measurements with previously published trigonometric parallaxes for an additional 44 objects and spectrophotometric distance estimates for another 7. For these 235 objects, we estimate temperatures, sift into five 150K-wide $T_{\rm eff}$ bins covering the range 300-1050K, determine the completeness limit for each, and compute space densities. To anchor the high-mass end of the brown dwarf mass spectrum, we compile a list of early- to mid-L dwarfs within 20 pc. We run simulations using various functional forms of the mass function passed through two different sets of evolutionary code to compute predicted distributions in $T_{\rm eff}$. The best fit of these predictions to our L, T, and Y observations is a simple power-law model with $\alpha \approx 0.6$ (where $dN/dM \propto M^{-\alpha}$), meaning that the slope of the field substellar mass function is in rough agreement with that found for brown dwarfs in nearby star forming regions and young clusters. Furthermore, we find that published versions of the log-normal form do not predict the steady rise seen in the space densities from 1050K to 350K. We also find that the low-mass cutoff to formation, if one exists, is lower than $\sim$5 $M_{Jup}$, which corroborates findings in young, nearby moving groups and implies that extremely low-mass objects have been forming over the lifetime of the Milky Way.

We show that a number of claims made in Myhrvold (2018) (hereafter M2018b) regarding the WISE data and thermal modeling of asteroids are incorrect. That paper provides thermal fit parameter outputs for only two of the about 150,000 object dataset and does not make a direct comparison to asteroids with diameters measured by other means to assess the quality of that work's thermal model. We are unable to reproduce the results for the two objects for which M2018b published its own thermal fit outputs, including diameter, albedo, beaming, and infrared albedo. In particular, the infrared albedos published in M2018b are unphysically low. [...] While there were some minor issues with consistency between tables due to clerical errors in the WISE/NEOWISE team's various papers and data release in the Planetary Data System, and a software issue that slightly increased diameter uncertainties in some cases, these issues do not substantially change the results and conclusions drawn from the data. We have shown in previous work and with updated analyses presented here that the effective spherical diameters for asteroids published to date are accurate to within the previously quoted minimum systematic 1-sigma uncertainty of about 10 percent when data of appropriate quality and quantity are available. Moreover, we show that the method used by M2018b to compare diameters between various asteroid datasets is incorrect and overestimates their differences. In addition, among other misconceptions in M2018b, we show that the WISE photometric measurement uncertainties are appropriately characterized and used by the WISE data processing pipeline and NEOWISE thermal modeling software. We show that the Near-Earth Asteroid Thermal Model (Harris 1998) employed by the NEOWISE team is a very useful model for analyzing infrared data to derive diameters and albedos when used properly.

We use optical and near-infrared spectroscopy to observe rest-UV emission lines and estimate the black hole mass of WISEA J224607.56-052634.9 (W2246-0526) at z = 4.601, the most luminous hot dust-obscured galaxy yet discovered by WISE. From the broad component of the MgII-2799A emission line, we measure a black hole mass of log (M_BH/M_sun) = 9.6 +- 0.4. The broad CIV-1549A line is asymmetric and significantly blue-shifted. The derived M_BH from the blueshift-corrected broad CIV line width agrees with the MgII result. From direct measurement using a well-sampled SED, the bolometric luminosity is 3.6 * 10^14 L_sun. The corresponding Eddington ratio for W2246-0526 is lambda_Edd = L_AGN / L_Edd = 2.8. This high Eddington ratio may reach the level where the luminosity is saturating due to photon trapping in the accretion flow, and be insensitive to the mass accretion rate. In this case, the M_BH growth rate in W2246-0526 would exceed the apparent accretion rate derived from the observed luminosity.

The Origins Space Telescope, one of four large Mission Concept studies sponsored by NASA for review in the 2020 US Astrophysics Decadal Survey, will open unprecedented discovery space in the infrared, unveiling our cosmic origins. We briefly describe in this article the key science themes and architecture for OST. With a sensitivity gain of up to a factor of 1,000 over any previous or planned mission, OST will open unprecedented discovery space, allow us to peer through an infrared window teeming with possibility. OST will fundamentally change our understanding of our cosmic origins - from the growth of galaxies and black holes, to uncovering the trail of water, to life signs in nearby Earth-size planets, and discoveries never imagined. Built to be highly adaptable, while addressing key science across many areas of astrophysics, OST will usher in a new era of infrared astronomy.

Y dwarfs provide a unique opportunity to study free-floating objects with masses $<$30 M$_{Jup}$ and atmospheric temperatures approaching those of known Jupiter-like exoplanets. Obtaining distances to these objects is an essential step towards characterizing their absolute physical properties. Using Spitzer/IRAC [4.5] images taken over baselines of $\sim$2-7 years, we measure astrometric distances for 22 late-T and early Y dwarfs, including updated parallaxes for 18 objects and new parallax measurements for 4 objects. These parallaxes will make it possible to explore the physical parameter space occupied by the coldest brown dwarfs. We also present the discovery of 6 new late-T dwarfs, updated spectra of two T dwarfs, and the reclassification of a new Y dwarf, WISE J033605.04$-$014351.0, based on Keck/NIRSPEC $J$-band spectroscopy. Assuming that effective temperatures are inversely proportional to absolute magnitude, we examine trends in the evolution of the spectral energy distributions of brown dwarfs with decreasing effective temperature. Surprisingly, the Y dwarf class encompasses a large range in absolute magnitude in the near- to mid-infrared photometric bandpasses, demonstrating a larger range of effective temperatures than previously assumed. This sample will be ideal for obtaining mid-infrared spectra with the James Webb Space Telescope because their known distances will make it easier to measure absolute physical properties.

We investigate the value of the near-infrared imaging from upcoming surveys for constraining the ellipticities of galaxies. We select galaxies between 0.5 < z < 3 that are brighter than expected Euclid sensitivity limits from the GOODS-S and N fields in CANDELS. The co-added CANDELS/HST V+I and J+H images are degraded in resolution and sensitivity to simulate Euclid-quality optical and near-infrared (NIR) images. We then run GALFIT on these simulated images and find that optical and NIR provide similar performance in measuring galaxy ellipticities at redshifts 0.5 < z < 3. At z > 1.0, the NIR-selected source density is higher by a factor of 1.4 and therefore the standard error in NIR-derived ellipticities is about 30% smaller, implying a more precise ellipticity measurement. The good performance of the NIR is mainly because galaxies have an intrinsically smoother light distribution in the NIR bands than in the optical, the latter tracing the clumpy star-forming regions. In addition, the NIR bands have a higher surface brightness per pixel than the optical images, while being less affected by dust attenuation. Despite the worse spatial sampling and resolution of Euclid NIR compared to optical, the NIR approach yields equivalent or more precise galaxy ellipticity measurements. If systematics that affect shape such as dithering strategy and point spread function undersampling can be mitigated, inclusion of the NIR can improve galaxy ellipticity measurements over all redshifts. This is particularly important for upcoming weak lensing surveys, such as with Euclid and WFIRST.

The WISE satellite surveyed the entire sky multiple times in four infrared wavelengths (3.4, 4.6, 12, and $22\,\mu$m; Wright et al. 2010). The unprecedented combination of coverage area and depth gives us the opportunity to measure the luminosity function of galaxies, one of the fundamental quantities in the study of them, at $2.4\ \mu$m to an unparalleled level of formal statistical accuracy in the near infrared. The big advantage of measuring luminosity functions at wavelengths in the window $\approx 2$ to $3.5\,\mu$m is that it correlates more closely to the total stellar mass in galaxies than others. In this paper we report on the parameters for the $2.4\,\mu$m luminosity function of galaxies obtained from applying the spectroluminosity functional based methods defined in Lake et al. (2017b) to the data sets described in Lake et al. (2017a) using the mean and covariance of $2.4\,\mu$m normalized SEDs from Lake & Wright (2016). In terms of single Schechter function parameters evaluated at the present epoch, the combined result is: $\phi_\star = 5.8 \pm [0.3_{\mathrm{stat}},\, 0.3_{\mathrm{sys}}] \times 10^{-3} \operatorname{Mpc}^{-3}$, $L_\star = 6.4 \pm [0.1_{\mathrm{stat}},\, 0.3_{\mathrm{sys}}] \times 10^{10}\, L_{2.4\,\mu\mathrm{m}\,\odot}$ ($M_\star = -21.67 \pm [0.02_{\mathrm{stat}},\, 0.05_{\mathrm{sys}}]\operatorname{AB\ mag}$), and $\alpha = -1.050 \pm [0.004_{\mathrm{stat}},\, 0.03_{\mathrm{sys}}]$, corresponding to a galaxy number density of $0.08\operatorname{Mpc}^{-3}$ brighter than $10^6\, L_{2.4\,\mu\mathrm{m}\,\odot}$ ($10^{-3} \operatorname{Mpc}^{-3}$ brighter than $L_\star$) and a $2.4\,\mu$m luminosity density equivalent to $3.8\times10^{8}\,L_{2.4\,\mu\mathrm{m}\,\odot}\operatorname{Mpc}^{-3}$. $\ldots$

The WISE satellite surveyed the entire sky multiple times in four infrared (IR) wavelengths ($3.4,\ 4.6,\ 12,$ and $22\, \mu$m, Wright et al. 2010). This all-sky IR photometric survey makes it possible to leverage many of the large publicly available spectroscopic redshift surveys to measure galaxy properties in the IR. While characterizing the cross-matching of WISE data to a single survey is a straightforward process, doing it with six different redshift surveys takes a fair amount of space to characterize adequately, because each survey has unique caveats and characteristics that need addressing. This work describes a data set that results from matching five public redshift surveys with the AllWISE data release, along with a reanalysis of the data described in Lake et al. 2012. The combined data set has an additional flux limit of $80\,\mu$Jy ($19.14$ AB mag) in WISE's W1 filter imposed in order to limit it to targets with high completeness and reliable photometry in the AllWISE data set. Consistent analysis of all of the data is only possible if the color bias discussed in Ilbert et al. (2004) is addressed (for example: the techniques explored in the first paper in this series Lake et al. 2017b). The sample defined herein is used in this paper's sequel paper, Lake et al. 2017a), to measure the luminosity function of galaxies at $2.4\, \mu$m rest frame wavelength, and the selection process of the sample is optimized for this purpose.

The astronomy community has at its disposal a large back catalog of public spectroscopic galaxy redshift surveys that can be used for the measurement of luminosity functions. Utilizing the back catalog with new photometric surveys to maximum efficiency requires modeling the color selection bias imposed on selection of target galaxies by flux limits at multiple wavelengths. The likelihood derived herein can address, in principle, all possible color selection biases through the use of a generalization of the luminosity function, $\Phi(L)$, over the space of all spectra: the spectro-luminosity functional, $\Psi[L_\nu]$. It is, therefore, the first estimator capable of simultaneously analyzing multiple redshift surveys in a consistent way. We also propose a new way of parametrizing the evolution of the classic Shechter function parameters, $L_\star$ and $\phi_\star$, that improves both the physical realism and statistical performance of the model. The techniques derived in this work will be used in an upcoming paper to measure the luminosity function of galaxies at the rest frame wavelength of $2.4\operatorname{\mu m}$ using the Widefield Infrared Survey Explorer (WISE).

While there are numerous criteria for photometrically identifying active galactic nuclei (AGNs), searches in the optical and UV tend to exclude galaxies that are highly dust obscured. This is problematic for constraining models of AGN evolution and estimating the AGN contribution to the cosmic X-ray and IR backgrounds, as highly obscured objects tend to be underrepresented in large-scale surveys. To address this, we identify potentially obscured AGNs using mid-IR color colors from the Wide-field Infrared Survey Explorer (WISE) catalog. This paper presents the results of optical spectroscopy of obscured AGN candidates using Keck DEIMOS, and their physical properties derived from these spectra. We find that a $W1-W2>0.8$ color criterion effectively selects AGNs with a higher median level of $E(B-V)$ extinction compared to the AGNs found in the SDSS DR7 survey. This optical extinction can be measured using SED modeling or by using $r-W1$ as a measure of optical to IR flux. We find that specific, targeted observations are necessary to find the most highly optically obscured AGNs, and that additional far-IR photometry is necessary to further constrain the dust properties of these AGNs.

The combination of the AKARI and WISE infrared all-sky surveys provides an unique opportunity to identify and characterize the most highly dust obscured AGNs in the universe. Dust-obscured AGNs are not easily detectable and potentially underrepresented in extragalactic surveys due to their high optical extinction, but are readily found in the WISE catalog due to their extremely red mid-IR colors. Combining these surveys with photometry from Pan-STARRS and Herschel, we use SED modeling to characterize the extinction and dust properties of these AGNs. From mid-IR WISE colors, we are able to compute bolometric corrections to AGN luminosities. Using AKARI's far-IR wavelength photometry and broadband AGN/galaxy spectral templates, we estimate AGN dust mass and temperature using simple analytic models with 3-4 parameters. Even without spectroscopic data, we can determine a number of AGN dust properties only using SED analysis. These methods, combined with the abundance of archival photometric data publically available, will be valuable for large-scale studies of dusty, IR-luminous AGNs.

Automated asteroid detection routines set requirements on the number of detections, signal-to-noise ratio, and the linearity of the expected motion in order to balance completeness, reliability, and time delay after data acquisition when identifying moving object tracklets. However, when the full-frame data from a survey are archived, they can be searched later for asteroids that were below the initial detection thresholds. We have conducted such a search of the first three years of the reactivated NEOWISE data, looking for near-Earth objects discovered by ground-based surveys that have previously unreported thermal infrared data. Using these measurements, we can then perform thermal modeling to measure the diameters and albedos of these objects. We present new physical properties for 116 Near-Earth Objects found in this search.

Quantifying the accuracy with which physical properties of asteroids can be determined from thermal modeling is critical to measuring the impact of infrared data on our understanding of asteroids. Previous work (Mainzer et al. 2011b) has used independently-derived diameters (from asteroid radar, occultations, and spacecraft visits) to test the accuracy of the NEOWISE diameter determinations. Here, we present a new and different method for bounding the actual NEOWISE diameter errors in the Main Belt based on our knowledge of the albedos of asteroid families. We show the 1 sigma relative diameter error for the Main Belt population must be less than 17.5% for the vast majority of objects. For a typical uncertainty on H magnitude of 0.2 mag, the relative error on diameter for the population would be ~10%.

One of the physical features of a dark-energy-dominated universe is the integrated Sachs-Wolfe (ISW) effect on the cosmic microwave background (CMB) radiation, which gives us a direct observational window to detect and study dark energy. The AllWISE data release of the Wide-field Infrared Survey Explorer (WISE) has a large number of point sources, which span over a wide redshift range including where the ISW effect is maximized. AllWISE data is thus very well-suited for the ISW effect studies. In this study, we cross-correlate AllWISE galaxy and active galactic nucleus (AGN) overdensities with the Wilkinson Microwave Anisotropy Probe CMB temperature maps to detect the ISW effect signal. We calibrate the biases for galaxies and AGNs by cross-correlating the galaxy and AGN overdensities with the Planck lensing convergence map. We measure the ISW effect signal amplitudes relative to the $\Lambda$CDM expectation of $A=1$ to be $A=1.18 \pm 0.36$ for galaxies and $A=0.64 \pm 0.74$ for AGNs . The detection significances for the ISW effect signal are $3.3\sigma$ and $0.9\sigma$ for galaxies and AGNs respectively giving a combined significance of $3.4\sigma$. Our result is in agreement with the $\Lambda$CDM model.

We identify new Y- and T-type brown dwarfs from the WISE All Sky data release using images obtained in filters that divide the traditional near-infrared H and J bands into two halves -- specifically CH4s & CH4l in the H and J2 & J3 in the J. This proves to be very effective at identifying cool brown dwarfs via the detection of their methane absorption, as well as providing preliminary classification using methane colours and WISE-to-near-infrared colours. New and updated calibrations between T/Y spectral types and CH4s-CH4l, J3-W2, and CH4s-W2 colours are derived, producing classification estimates good to a few spectral sub-types. We present photometry for a large sample of T and Y dwarfs in these filters, together with spectroscopy for 23 new ultra-cool dwarfs - two Y dwarfs and twenty one T dwarfs. We identify a further 8 new cool brown dwarfs, which we have high confidence are T dwarfs based on their methane photometry. We find that, for objects observed on a 4m-class telescope at J band magnitudes of ~20 or brighter, CH4s-CH4l is the more powerful colour for detecting objects and then estimating spectral types. Due to the lower sky background in the J-band, the J3 and J2 bands are more useful for identifying fainter cool dwarfs at J>22. The J3-J2 colour is poor at estimating spectral types. But fortunately, once J3-J2 confirms that an object is a cool dwarf, the J3-W2 colour is very effective at estimating approximate spectral types.

From the entire dataset of comets observed by NEOWISE, we have analyzed 11 different Halley-Type Comets (HTCs) for dust production rates, CO+CO2 production rates, and nucleus sizes. Incorporating HTCs from previous studies and multiple comet visits we have a total of 21 stacked visits, 13 of which are active and 8 for which we calculated upper limits of production. We determined the nucleus sizes of 27P, P/2006 HR30, P/2012 NJ, and C/2016 S1. Furthermore, we analyzed the relationships between dust production and heliocentric distance, and gas production and heliocentric distance. We concluded that for this population of HTCs, ranging in heliocentric distance from 1.21 AU to 2.66 AU, there was no significant correlation between dust production and heliocentric distance, nor gas production and heliocentric distance.

Hot, Dust-Obscured Galaxies, or "Hot DOGs", are a rare, dusty, hyperluminous galaxy population discovered by the WISE mission. Predominantly at redshifts 2-3, they include the most luminous known galaxies in the universe. Their high luminosities likely come from accretion onto highly obscured super massive black holes (SMBHs). We have conducted a pilot survey to measure the SMBH masses of five z~2 Hot DOGs via broad H_alpha emission lines, using Keck/MOSFIRE and Gemini/FLAMINGOS-2. We detect broad H_alpha emission in all five Hot DOGs. We find substantial corresponding SMBH masses for these Hot DOGs (~ 10^{9} M_sun), and their derived Eddington ratios are close to unity. These z~2 Hot DOGs are the most luminous AGNs at given BH masses, suggesting they are accreting at the maximum rates for their BHs. A similar property is found for known z~6 quasars. Our results are consistent with scenarios in which Hot DOGs represent a transitional, high-accretion phase between obscured and unobscured quasars. Hot DOGs may mark a special evolutionary stage before the red quasar and optical quasar phases, and they may be present at other cosmic epochs.

The Near-Earth Object Wide-field Infrared Survey Explorer (NEOWISE) reactivation mission has completed its third year of surveying the sky in the thermal infrared for near-Earth asteroids and comets. NEOWISE collects simultaneous observations at 3.4 um and 4.6 um of solar system objects passing through its field of regard. These data allow for the determination of total thermal emission from bodies in the inner solar system, and thus the sizes of these objects. In this paper we present thermal model fits of asteroid diameters for 170 NEOs and 6110 MBAs detected during the third year of the survey, as well as the associated optical geometric albedos. We compare our results with previous thermal model results from NEOWISE for overlapping sample sets, as well as diameters determined through other independent methods, and find that our diameter measurements for NEOs agree to within 26% (1-sigma) of previously measured values. Diameters for the MBAs are within 17% (1-sigma). This brings the total number of unique near-Earth objects characterized by the NEOWISE survey to 541, surpassing the number observed during the fully cryogenic mission in 2010.

We investigate extremely luminous dusty galaxies in the environments around WISE-selected hot dust obscured galaxies (Hot DOGs) and WISE/radio-selected active galactic nuclei (AGNs) at average redshifts of z = 2.7 and z = 1.7, respectively. Previous observations have detected overdensities of companion submillimetre-selected sources around 10 Hot DOGs and 30 WISE/radio AGNs, with overdensities of ~ 2 - 3 and ~ 5 - 6 , respectively. We find that the space densities in both samples to be overdense compared to normal star-forming galaxies and submillimetre galaxies (SMGs) in the SCUBA-2 Cosmology Legacy Survey (S2CLS). Both samples of companion sources have consistent mid-IR colours and mid-IR to submm ratios as SMGs. The brighter population around WISE/radio AGNs could be responsible for the higher overdensity reported. We also find the star formation rate density (SFRDs) are higher than the field, but consistent with clusters of dusty galaxies. WISE-selected AGNs appear to be good signposts for protoclusters at high redshift on arcmin scales. The results reported here provide an upper limit to the strength of angular clustering using the two-point correlation function. Monte Carlo simulations show no angular correlation, which could indicate protoclusters on scales larger than the SCUBA-2 1.5arcmin scale maps.

The rapid accumulation of thermal infrared observations and shape models of asteroids has led to increased interest in thermophysical modeling. Most of these infrared observations are unresolved. We consider what fraction of an asteroid's surface area contributes the bulk of the emitted thermal flux for two model asteroids of different shapes over a range of thermal parameters. The resulting observed surface in the infrared is generally more fragmented than the area observed in visible wavelengths, indicating high sensitivity to shape. For objects with low values of the thermal parameter, small fractions of the surface contribute the majority of thermally emitted flux. Calculating observed areas could enable the production of spatially-resolved thermal inertia maps from non-resolved observations of asteroids.

We present optical and near-infrared spectroscopy of WISEA J061543.91$-$124726.8, which we rediscovered as a high motion object in the AllWISE survey. The spectra of this object are unusual; while the red optical ($\lambda >$ 7,000 \AA) and near-infrared spectra exhibit characteristic TiO, VO, and H$_{2}$O bands of a late-M dwarf, the blue portion of its optical spectrum shows a significant excess of emission relative to late-M type templates. The excess emission is relatively featureless, with the exception of a prominent and very broad Na I D doublet. We find that no single, ordinary star can reproduce these spectral characteristics. The most likely explanation is an unresolved binary system of an M7 dwarf and a cool white dwarf. The flux of a cool white dwarf drops in the optical red and near-infrared, due to collision-induced absorption, thus allowing the flux of a late-M dwarf to show through. This scenario, however, does not explain the Na D feature, which is unlike that of any known white dwarf, but which could perhaps be explained via unusual abundance or pressure conditions.

In this paper we provide formulae that can be used to determine the uncertainty contributed to a measurement by a $K$-correction and, thus, valuable information about which flux measurement will provide the most accurate $K$-corrected luminosity. All of this is done at the level of a Gaussian approximation of the statistics involved, that is, where the galaxies in question can be characterized by a mean spectral energy distribution (SED) and a covariance function (spectral 2-point function). This paper also includes approximations of the SED mean and covariance for galaxies, and the three common subclasses thereof, based on applying the templates from Assef et al. (2010) to the objects in zCOSMOS bright 10k (Lilly et al. 2009) and photometry of the same field from Capak et al. (2007), Sanders et al. (2007), and the AllWISE source catalog.

This paper describes the beginning of the Far-Infrared Surveyor mission study for NASA's Astrophysics Decadal 2020. We describe the scope of the study, and the open process approach of the Science and Technology Definition Team. We are currently developing the science cases and provide some preliminary highlights here. We note key areas for technological innovation and improvements necessary to make a Far-Infrared Surveyor mission a reality.

The Near-Earth Object Wide-Field Infrared Survey Explorer (NEOWISE) mission continues to detect, track, and characterize minor planets. We present diameters and albedos calculated from observations taken during the second year since the spacecraft was reactivated in late 2013. These include 207 near-Earth asteroids and 8,885 other asteroids. $84\%$ of the near-Earth asteroids did not have previously measured diameters and albedos by the NEOWISE mission. Comparison of sizes and albedos calculated from NEOWISE measurements with those measured by occultations, spacecraft, and radar-derived shapes shows accuracy consistent with previous NEOWISE publications. Diameters and albedos fall within $ \pm \sim20\%$ and $\pm\sim40\%$, 1-sigma, respectively, of those measured by these alternate techniques. NEOWISE continues to preferentially discover near-Earth objects which are large ($>100$ m), and have low albedos.

The cryogenic WISE mission in 2010 was extremely sensitive to asteroids and not biased against detecting dark objects. The albedos of 428 Near Earth Asteroids (NEAs) observed by WISE during its fully cryogenic mission can be fit quite well by a 3 parameter function that is the sum of two Rayleigh distributions. The Rayleigh distribution is zero for negative values, and follows $f(x) = x \exp[-x^2/(2\sigma^2)]/\sigma^2$ for positive x. The peak value is at x=\sigma, so the position and width are tied together. The three parameters are the fraction of the objects in the dark population, the position of the dark peak, and the position of the brighter peak. We find that 25.3% of the NEAs observed by WISE are in a very dark population peaking at $p_V = 0.03$, while the other 74.7% of the NEAs seen by WISE are in a moderately dark population peaking at $p_V = 0.168$. A consequence of this bimodal distribution is that the Congressional mandate to find 90% of all NEAs larger than 140 m diameter cannot be satisfied by surveying to H=22 mag, since a 140 m diameter asteroid at the very dark peak has H=23.7 mag, and more than 10% of NEAs are darker than p_V = 0.03.

We use the AllWISE Data Release to continue our search for WISE-detected motions. In this paper, we publish another 27,846 motion objects, bringing the total number to 48,000 when objects found during our original AllWISE motion survey are included. We use this list, along with the lists of confirmed WISE-based motion objects from the recent papers by Luhman and by Schneider et al. and candidate motion objects from the recent paper by Gagne et al. to search for widely separated, common-proper-motion systems. We identify 1,039 such candidate systems. All 48,000 objects are further analyzed using color-color and color-mag plots to provide possible characterizations prior to spectroscopic follow-up. We present spectra of 172 of these, supplemented with new spectra of 23 comparison objects from the literature, and provide classifications and physical interpretations of interesting sources. Highlights include: (1) the identification of three G/K dwarfs that can be used as standard candles to study clumpiness and grain size in nearby molecular clouds because these objects are currently moving behind the clouds, (2) the confirmation/discovery of several M, L, and T dwarfs and one white dwarf whose spectrophotometric distance estimates place them 5-20 pc from the Sun, (3) the suggestion that the Na 'D' line be used as a diagnostic tool for interpreting and classifying metal-poor late-M and L dwarfs, (4) the recognition of a triple system including a carbon dwarf and late-M subdwarf, for which model fits of the late-M subdwarf (giving [Fe/H] ~ -1.0) provide a measured metallicity for the carbon star, and (5) a possible 24-pc-distant K5 dwarf + peculiar red L5 system with an apparent physical separation of 0.1 pc.

We present WISEA J114724.10$-$204021.3, a young, low-mass, high probability member of the TW Hya association. WISEA J114724.10$-$204021.3 was discovered based on its red AllWISE color (W1$-$W2 = 0.63 mag) and extremely red 2MASS $J-K_{\rm S}$ color ($>$ 2.64 mag), the latter of which is confirmed with near-infrared photometry from the VISTA Hemisphere Survey ($J-K_{\rm S}$ = 2.57$\pm$0.03). Follow-up near-infrared spectroscopy shows a spectral type of L7 $\pm$ 1 as well as several spectroscopic indicators of youth. These include a peaked $H$-band shape and a steeper $K$-band slope, traits typically attributed to low surface gravity. The sky position, proper motion, and distance estimates of WISEA J114724.10$-$204021.3 are all consistent with membership in the $\sim$10 Myr old TW Hya association. Using the age of the TW Hya association and evolutionary models, we estimate the mass of WISEA J114724.10$-$204021.3 to be 5$-$13 $M_{\rm Jup}$, making it one of the youngest and lowest mass free-floating objects yet discovered in the Solar neighborhood.

The 163 comets observed during the WISE/NEOWISE prime mission represent the largest infrared survey to date of comets, providing constraints on dust, nucleus sizes, and CO+CO2 production. We present detailed analyses of the WISE/NEOWISE comet discoveries, and discuss observations of the active comets showing 4.6 $\mu$m band excess. We find a possible relation between dust and CO+CO2 production, as well as possible differences in the sizes of long and short period comet nuclei.

We report Markov chain Monte Carlo fits of the thermophysical model of Wright (2007) to the fluxes of 10 asteroids which have been observed by both WISE and NEOWISE. This model is especially useful when one has observations of an asteroid at multiple epochs, as it takes advantage of the views of different local times and latitudes to determine the spin axis and the thermal parameter. Many of the asteroids NEOWISE observes will have already been imaged by WISE, so this proof of concept shows there is an opportunity to use a rotating cratered thermophysical model to determine surface thermal properties of a large number of asteroids.

We present preliminary diameters and albedos for 7,959 asteroids detected in the first year of the NEOWISE Reactivation mission. 201 are near-Earth asteroids (NEAs). 7,758 are Main Belt or Mars-crossing asteroids. 17% of these objects have not been previously characterized using WISE or NEOWISE thermal measurements. Diameters are determined to an accuracy of ~20% or better. If good-quality H magnitudes are available, albedos can be determined to within ~40% or better.

We present 20 WISE-selected galaxies with bolometric luminosities L_bol > 10^14 L_sun, including five with infrared luminosities L_IR = L(rest 8-1000 micron) > 10^14 L_sun. These "extremely luminous infrared galaxies," or ELIRGs, were discovered using the "W1W2-dropout" selection criteria which requires marginal or non-detections at 3.4 and 4.6 micron (W1 and W2, respectively) but strong detections at 12 and 22 micron in the WISE survey. Their spectral energy distributions are dominated by emission at rest-frame 4-10 micron, suggesting that hot dust with T_d ~ 450K is responsible for the high luminosities. These galaxies are likely powered by highly obscured AGNs, and there is no evidence suggesting these systems are beamed or lensed. We compare this WISE-selected sample with 116 optically selected quasars that reach the same L_bol level, corresponding to the most luminous unobscured quasars in the literature. We find that the rest-frame 5.8 and 7.8 micron luminosities of the WISE-selected ELIRGs can be 30-80% higher than that of the unobscured quasars. The existence of AGNs with L_bol > 10^14 L_sun at z > 3 suggests that these supermassive black holes are born with large mass, or have very rapid mass assembly. For black hole seed masses ~ 10^3 M_sun, either sustained super-Eddington accretion is needed, or the radiative efficiency must be <15%, implying a black hole with slow spin, possibly due to chaotic accretion.

We present a sample of brown dwarfs identified with the {\it Wide-field Infrared Survey Explorer} (WISE) for which we have obtained {\it Hubble Space Telescope} ({\it HST}) Wide Field Camera 3 (WFC3) near-infrared grism spectroscopy. The sample (twenty-two in total) was observed with the G141 grism covering 1.10$-$1.70 $\mu$m, while fifteen were also observed with the G102 grism, which covers 0.90$-$1.10 $\mu$m. The additional wavelength coverage provided by the G102 grism allows us to 1) search for spectroscopic features predicted to emerge at low effective temperatures (e.g.\ ammonia bands) and 2) construct a smooth spectral sequence across the T/Y boundary. We find no evidence of absorption due to ammonia in the G102 spectra. Six of these brown dwarfs are new discoveries, three of which are found to have spectral types of T8 or T9. The remaining three, WISE J082507.35$+$280548.5 (Y0.5), WISE J120604.38$+$840110.6 (Y0), and WISE J235402.77$+$024015.0 (Y1) are the nineteenth, twentieth, and twenty-first spectroscopically confirmed Y dwarfs to date. We also present {\it HST} grism spectroscopy and reevaluate the spectral types of five brown dwarfs for which spectral types have been determined previously using other instruments.

We present the all-sky Planck catalogue of Sunyaev-Zeldovich (SZ) sources detected from the 29 month full-mission data. The catalogue (PSZ2) is the largest SZ-selected sample of galaxy clusters yet produced and the deepest all-sky catalogue of galaxy clusters. It contains 1653 detections, of which 1203 are confirmed clusters with identified counterparts in external data-sets, and is the first SZ-selected cluster survey containing > $10^3$ confirmed clusters. We present a detailed analysis of the survey selection function in terms of its completeness and statistical reliability, placing a lower limit of 83% on the purity. Using simulations, we find that the Y5R500 estimates are robust to pressure-profile variation and beam systematics, but accurate conversion to Y500 requires. the use of prior information on the cluster extent. We describe the multi-wavelength search for counterparts in ancillary data, which makes use of radio, microwave, infra-red, optical and X-ray data-sets, and which places emphasis on the robustness of the counterpart match. We discuss the physical properties of the new sample and identify a population of low-redshift X-ray under- luminous clusters revealed by SZ selection. These objects appear in optical and SZ surveys with consistent properties for their mass, but are almost absent from ROSAT X-ray selected samples.

We have carried out simulations to predict the performance of a new space-based telescopic survey operating at thermal infrared wavelengths that seeks to discover and characterize a large fraction of the potentially hazardous near-Earth asteroid (NEA) population. Two potential architectures for the survey were considered: one located at the Earth-Sun L1 Lagrange point, and one in a Venus-trailing orbit. A sample cadence was formulated and tested, allowing for the self-follow-up necessary for objects discovered in the daytime sky on Earth. Synthetic populations of NEAs with sizes >=140 m in effective spherical diameter were simulated using recent determinations of their physical and orbital properties. Estimates of the instrumental sensitivity, integration times, and slew speeds were included for both architectures assuming the properties of new large-format 10 um detector arrays capable of operating at ~35 K. Our simulation included the creation of a preliminary version of a moving object processing pipeline suitable for operating on the trial cadence. We tested this pipeline on a simulated sky populated with astrophysical sources such as stars and galaxies extrapolated from Spitzer and WISE data, the catalog of known minor planets (including Main Belt asteroids, comets, Jovian Trojans, etc.), and the synthetic NEA model. Trial orbits were computed for simulated position-time pairs extracted from the synthetic surveys to verify that the tested cadence would result in orbits suitable for recovering objects at a later time. Our results indicate that the Earth-Sun L1 and Venus-trailing surveys achieve similar levels of integral completeness for potentially hazardous asteroids larger than 140 m; placing the telescope in an interior orbit does not yield an improvement in discovery rates. This work serves as a necessary first step for the detailed planning of a next-generation NEA survey.

In recent years brown dwarfs have been extended to a new Y-dwarf class with effective temperatures colder than 500K and masses in the range 5-30 Jupiter masses. They fill a crucial gap in observable atmospheric properties between the much colder gas-giant planets of our own Solar System (at around 130K) and both hotter T-type brown dwarfs and the hotter planets that can be imaged orbiting young nearby stars (both with effective temperatures of in the range 1500-1000K). Distance measurements for these objects deliver absolute magnitudes that make critical tests of our understanding of very cool atmospheres. Here we report new distances for nine Y dwarfs and seven very-late T dwarfs. These reveal that Y dwarfs do indeed represent a continuation of the T dwarf sequence to both fainter luminosities and cooler temperatures. They also show that the coolest objects display a large range in absolute magnitude for a given photometric colour. The latest atmospheric models show good agreement with the majority of these Y dwarf absolute magnitudes. This is also the case for WISE0855-0714 the coldest and closest brown dwarf to the Sun, which shows evidence for water ice clouds. However, there are also some outstanding exceptions, which suggest either binarity or the presence of condensate clouds. The former is readily testable with current adaptive optics facilities. The latter would mean that the range of cloudiness in Y dwarfs is substantial with most hosting almost no clouds -- while others have dense clouds making them prime targets for future variability observations to study cloud dynamics.

WISE has discovered an extraordinary population of hyper-luminous dusty galaxies which are faint in the two bluer passbands ($3.4\, \mu$m and $4.6\, \mu$m) but are bright in the two redder passbands of WISE ($12\, \mu$m and $22\, \mu$m). We report on initial follow-up observations of three of these hot, dust-obscured galaxies, or Hot DOGs, using the CARMA and SMA interferometer arrays at submm/mm wavelengths. We report continuum detections at $\sim$ 1.3 mm of two sources (WISE J014946.17+235014.5 and WISE J223810.20+265319.7, hereafter W0149+2350 and W2238+2653, respectively), and upper limits to CO line emission at 3 mm in the observed frame for two sources (W0149+2350 and WISE J181417.29+341224.8, hereafter W1814+3412). The 1.3 mm continuum images have a resolution of 1-2 arcsec and are consistent with single point sources. We estimate the masses of cold dust are 2.0$\times 10^{8} M_{\odot}$ for W0149+2350 and 3.9$\times 10^{8} M_{\odot}$ for W2238+2653, comparable to cold dust masses of luminous quasars. We obtain 2$\sigma$ upper limits to the molecular gas masses traced by CO, which are 3.3$\times 10^{10} M_{\odot}$ and 2.3$\times 10^{10} M_{\odot}$ for W0149+2350 and W1814+3412, respectively. We also present high-resolution, near-IR imaging with WFC3 on the Hubble Space Telescope for W0149+2653 and with NIRC2 on Keck for W2238+2653. The near-IR images show morphological structure dominated by a single, centrally condensed source with effective radius less than 4 kpc. No signs of gravitational lensing are evident.

The Wide-field Infrared Survey Explorer (WISE) spacecraft has been reactivated as NEOWISE-R to characterize and search for Near Earth Objects. The brown dwarf WISE J085510.83-071442.5 has now been reobserved by NEOWISE-R, and we confirm the results of Luhman (2014b), who found a very low effective temperature ($\approx 250$ K), a very high proper motion (8.1 +/- 0.1 arcsec/yr) , and a large parallax (454 +/- 45 mas). The large proper motion has separated the brown dwarf from the background sources that influenced the 2010 WISE data, allowing a measurement of a very red WISE color of W1-W2 $> 3.9$ mag. A re-analysis of the 2010 WISE astrometry using only the W2 band, combined with the new NEOWISE-R 2014 position, gives an improved parallax of 448 +/- 33 mas and proper motion of 8.08 +/- 0.05\; arcsec/yr. These are all consistent with Luhman (2014b).

NASA's Wide-field Infrared Survey Explorer (WISE) spacecraft has been brought out of hibernation and has resumed surveying the sky at 3.4 and 4.6 um. The scientific objectives of the NEOWISE reactivation mission are to detect, track, and characterize near-Earth asteroids and comets. The search for minor planets resumed on December 23, 2013, and the first new near-Earth object (NEO) was discovered six days later. As an infrared survey, NEOWISE detects asteroids based on their thermal emission and is equally sensitive to high and low albedo objects; consequently, NEOWISE-discovered NEOs tend to be large and dark. Over the course of its three-year mission, NEOWISE will determine radiometrically-derived diameters and albedos for approximately 2000 NEOs and tens of thousands of Main Belt asteroids. The 32 months of hibernation have had no significant effect on the mission's performance. Image quality, sensitivity, photometric and astrometric accuracy, completeness, and the rate of minor planet detections are all essentially unchanged from the prime mission's post-cryogenic phase.

We present SCUBA-2 850um submillimetre (submm) observations of the fields of 10 dusty, luminous galaxies at z ~ 1.7 - 4.6, detected at 12um and/or 22um by the WISE all-sky survey, but faint or undetected at 3.4um and 4.6um; dubbed hot, dust-obscured galaxies (Hot DOGs). The six detected targets all have total infrared luminosities greater than 10^13 L_sun, with one greater than 10^14 L_sun. Their spectral energy distributions (SEDs) are very blue from mid-infrared to submm wavelengths and not well fitted by standard AGN SED templates, without adding extra dust extinction to fit the WISE 3.4um and 4.6um data. The SCUBA-2 850um observations confirm that the Hot DOGs have less cold and/or more warm dust emission than standard AGN templates, and limit an underlying extended spiral or ULIRG-type galaxy to contribute less than about 2% or 55% of the typical total Hot DOG IR luminosity, respectively. The two most distant and luminous targets have similar observed submm to mid-infrared ratios to the rest, and thus appear to have even hotter SEDs. The number of serendipitous submm galaxies (SMGs) detected in the 1.5-arcmin-radius SCUBA-2 850um maps indicates there is a significant over-density of serendipitous sources around Hot DOGs. These submm observations confirm that the WISE-selected ultra-luminous galaxies have very blue mid-infrared to submm SEDs, suggesting that they contain very powerful AGN, and are apparently located in unusual arcmin-scale overdensities of very luminous dusty galaxies.

We report on a NuSTAR and XMM-Newton program that has observed a sample of three extremely luminous, heavily obscured WISE-selected AGN at z~2 in a broad X-ray band (0.1 - 79 keV). The parent sample, selected to be faint or undetected in the WISE 3.4um (W1) and 4.6um (W2) bands but bright at 12um (W3) and 22um (W4), are extremely rare, with only ~1000 so-called W1W2-dropouts across the extragalactic sky. Optical spectroscopy reveals typical redshifts of z~2 for this population, implying rest-frame mid-IR luminosities of L(6um)~6e46 erg/s and bolometric luminosities that can exceed L(bol)~1e14 L(sun). The corresponding intrinsic, unobscured hard X-ray luminosities are L(2-10)~4e45 erg/s for typical quasar templates. These are amongst the most luminous AGN known, though the optical spectra rarely show evidence of a broad-line region and the selection criteria imply heavy obscuration even at rest-frame 1.5um. We designed our X-ray observations to obtain robust detections for gas column densities N(H)<1e24 /cm2. In fact, the sources prove to be fainter than these predictions. Two of the sources were observed by both NuSTAR and XMM-Newton, with neither being detected by NuSTAR and one being faintly detected by XMM-Newton. A third source was observed only with XMM-Newton, yielding a faint detection. The X-ray data require gas column densities N(H)>1e24 /cm2, implying the sources are extremely obscured, consistent with Compton-thick, luminous quasars. The discovery of a significant population of heavily obscured, extremely luminous AGN does not conform to the standard paradigm of a receding torus, in which more luminous quasars are less likely to be obscured. If a larger sample conforms with this finding, then this suggests an additional source of obscuration for these extreme sources.

We have determined astrometric positions for 15 WISE-discovered late-type brown dwarfs (6 T8-9 and 9 Y dwarfs) using the Keck II telescope, the Spitzer Space Telescope, and the Hubble Space Telescope. Combining data from 8 to 20 epochs we derive parallactic and proper motions for these objects which put the majority within 15 parsecs. For ages greater than a few Gyr, as suggested from kinematic considerations, we find masses of 10-30 M$_{Jup}$ based on standard models for the evolution of low mass objects with a range of mass estimates for individual objects depending on the model in question. Three of the coolest objects have effective temperatures $\sim$ 350 K and inferred masses of 10-15 M$_{Jup}$. Our parallactic distances confirm earlier photometric estimates (Kirkpatrick et al. 2012) and direct measurements (Marsh et al. 2013, Beichman et al. 2013, Dupuy and Krauss 2013) and suggest that the number of objects with masses below about 15 M$_{Jup}$ must be flat or declining relative to higher mass objects. The masses of the coldest Y dwarfs may be similar to those inferred for recently imaged planet-mass companions to nearby young stars. Objects in this mass range, which appear to be rare in both the interstellar and proto-planetary environments, may both have formed via gravitational fragmentation: the brown dwarfs in interstellar clouds and companion objects in a protoplanetary disk. In both cases, however, the fact that objects in this mass range are relatively infrequent, suggests that this mechanism must be inefficient in both environments.

The AllWISE processing pipeline has measured motions for all objects detected on WISE images taken between 2010 January and 2011 February. In this paper, we discuss new capabilities made to the software pipeline in order to make motion measurements possible, and we characterize the resulting data products for use by future researchers. Using a stringent set of selection criteria, we find 22,445 objects that have significant AllWISE motions, of which 3,525 have motions that can be independently confirmed from earlier 2MASS images yet lack any published motions in SIMBAD. Another 58 sources lack 2MASS counterparts and are presented as motion candidates only. Limited spectroscopic follow-up of this list has already revealed eight new L subdwarfs. These may provide the first hints of a "subdwarf gap" at mid-L types that would indicate the break between the stellar and substellar populations at low metallicities (i.e., old ages). Another object in the motion list -- WISEA J154045.67-510139.3 -- is a bright (J ~ 9 mag) object of type M6; both the spectrophotometric distance and a crude preliminary parallax place it ~6 pc from the Sun. We also compare our list of motion objects to the recently published list of 762 WISE motion objects from Luhman (2014). While these first large motion studies with WISE data have been very successful in revealing previously overlooked nearby dwarfs, both studies missed objects that the other found, demonstrating that many other nearby objects likely await discovery in the AllWISE data products.

While quality checking a new motion-aware co-addition of all 12.5 months of WISE data, the source WISE J070720.48+170533.0 was found to have moved 0.9 arc-sec in 6 months. Backtracking this motion allowed us to identify this source as 2MASS J07071961+1705464 and with several entries in the USNO B catalog. An astrometric fit to these archival data gives a proper motion of 1793 +/- 2 mas/yr and a parallax of 35 +/- 42 mas. Photometry from WISE, 2MASS and the POSS can be fit reasonably well by a blackbody with T = 3658 K and an angular radius of 4.36E-11 radians. No clear evidence of H_2 collision-induced absorption is seen in the near-IR. An optical spectrum shows broad deep CaH bands at 638 and 690 nm, broad deep Na D at 598.2 nm, and weak or absent TiO, indicating that this source is an ultra-subdwarf M star with a radial velocity about -21 +/- 18 km/sec relative to the Sun. Given its apparent magnitude, the distance is about 39 +/- 9 pc and the tangential velocity is probably about 330 km/sec, but a more precise parallax is needed to be certain.

On behalf of the WISE Science team, we present the discovery of a class of distant dust-enshrouded galaxies with extremely high luminosity. These galaxies are selected to have extreme red colors in the mid-IR using NASA's Wide-field Infrared Survey Explorer (WISE). They are faint in the optical and near-IR, predominantly at z=2-4, and with IR luminosity > $10^{13}\, L_{Sun}$, making them Hyper-Luminous Infrared Galaxies (HyLIRGs). SEDs incorporating the WISE, Spitzer, and Herschel PACS and SPIRE photometry indicate hot dust dominates the bolometric luminosity, presumably powered by AGN. Preliminary multi-wavelength follow-up suggests that they are different from normal populations in the local M-sigma relation. Their low source density implies that these objects are either intrinsically rare, or a short-lived phase in a more numerous population. If the latter is the case, these hot, dust-enshrouded galaxies may be an early stage in the interplay between AGN and galaxies.

Only a very small fraction of the asteroid population at size scales comparable to the object that exploded over Chelyabinsk, Russia has been discovered to date, and physical properties are poorly characterized. We present previously unreported detections of 106 close approaching near-Earth objects (NEOs) by the Wide-field Infrared Survey Explorer mission's NEOWISE project. These infrared observations constrain physical properties such as diameter and albedo for these objects, many of which are found to be smaller than 100 m. Because these objects are intrinsically faint, they were detected by WISE during very close approaches to the Earth, often at large apparent on-sky velocities. We observe a trend of increasing albedo with decreasing size, but as this sample of NEOs was discovered by visible light surveys, it is likely that selection biases against finding small, dark NEOs influence this finding.

We have discovered a wide separation (188.5") T8 subdwarf companion to the sdM1.5+WD binary Wolf 1130. Companionship of WISE J200520.38+542433.9 is verified through common proper motion over a ~3 year baseline. Wolf 1130 is located 15.83 +/- 0.96 parsecs from the Sun, placing the brown dwarf at a projected separation of ~3000 AU. Near-infrared colors and medium resolution (R~2000-4000) spectroscopy establish the uniqueness of this system as a high-gravity, low-metallicity benchmark. Although there are a number of low-metallicity T dwarfs in the literature, WISE J200520.38+542433.9 has the most extreme inferred metallicity to date with [Fe/H] = -0.64 +/- 0.17 based on Wolf 1130. Model comparisons to this exemplar late-type subdwarf support it having an old age, a low metallicity, and a small radius. However, the spectroscopic peculiarities of WISE J200520.38+542433.9 underscore the importance of developing the low-metallicity parameter space of the most current atmospheric models.

In the local Universe, 10% of massive elliptical galaxies are observed to exhibit a peculiar property: a substantial excess of ultraviolet emission (UVX) over what is expected from their old, red stellar populations. Several origins for the UVX have been proposed, including a population of hot young stars, or a population of old, blue horizontal branch or extended horizontal branch (BHB or EHB) stars that have undergone substantial mass loss from their outer atmospheres. We explore the radial distribution of ultraviolet excess (UVX) in a selection of 49 nearby E/S0-type galaxies by measuring the extended photometry in the UV-midIR with GALEX, SDSS and WISE. We compare UV/optical and UV/mid-IR colors with the Flexible Stellar Population Synthesis with EHB models (Conroy & Gunn 2010). We find that combined WISE mid-IR and GALEX UV colors are more effective in distinguishing models than optical colors, and that the UV/mid-IR combination is sensitive to EHB fraction. There are strong color gradients with the outer radii bluer than the inner half-light radii by ~1 magnitude. This color difference is easily accounted for with a BHB fraction increase of 0.25 with radius. We estimated the average ages for the inner and outer radii are 7.0+/-0.3 Gyr, and 6.2+/-0.2 Gyr, respectively, with the implication that the outer regions are likely to have formed ~1 Gyr after the inner regions. Additionally, we find that metallicity gradients are likely not a significant factor in the color difference. The separation of color between the inner and outer regions, which agrees with a specific stellar population difference (e.g., higher EHB populations), and the ~0.5-2 Gyr age difference suggests multi-stage formation. Our results are best explained by inside-out formation: rapid star formation within the core at early epochs (>4 Gyr ago) and at least one later stage starburst event coinciding with z~1

The Wide-field Infrared Survey Explorer (WISE) observed 52 Centaurs and Scattered Disk Objects in the thermal infrared, including the 15 discoveries that were new. We present analyses of these observations to estimate sizes and mean optical albedos. We find mean albedos of 0.08 +/- 0.04 for the entire data set. Thermal fits yield average beaming parameters of 0.9 +/- 0.2 that are similar for both SDO and Centaur sub-classes. Biased cumulative size distributions yield size-frequency distribution power law indices ~ -1.7 +/- 0.3. The data also reveal a relation between albedo and color at the 3-sigma level. No significant relation between diameter and albedos is found.

We present the final nine-year maps and basic results from the WMAP mission. We provide new nine-year full sky temperature maps that were processed to reduce the asymmetry of the effective beams. Temperature and polarization sky maps are examined to separate CMB anisotropy from foreground emission, and both types of signals are analyzed in detail. The WMAP mission has resulted in a highly constrained LCDM cosmological model with precise and accurate parameters in agreement with a host of other cosmological measurements. When WMAP data are combined with finer scale CMB, baryon acoustic oscillation, and Hubble constant measurements, we find that Big Bang nucleosynthesis is well supported and there is no compelling evidence for a non-standard number of neutrino species (3.84+/-0.40). The model fit also implies that the age of the universe is 13.772+/-0.059 Gyr, and the fit Hubble constant is H0 = 69.32+/-0.80 km/s/Mpc. Inflation is also supported: the fluctuations are adiabatic, with Gaussian random phases; the detection of a deviation of the scalar spectral index from unity reported earlier by WMAP now has high statistical significance (n_s = 0.9608+/-0.0080); and the universe is close to flat/Euclidean, Omega_k = -0.0027 (+0.0039/-0.0038). Overall, the WMAP mission has resulted in a reduction of the cosmological parameter volume by a factor of 68,000 for the standard six-parameter LCDM model, based on CMB data alone. For a model including tensors, the allowed seven-parameter volume has been reduced by a factor 117,000. Other cosmological observations are in accord with the CMB predictions, and the combined data reduces the cosmological parameter volume even further. With no significant anomalies and an adequate goodness-of-fit, the inflationary flat LCDM model and its precise and accurate parameters rooted in WMAP data stands as the standard model of cosmology.

We present cosmological parameter constraints based on the final nine-year WMAP data, in conjunction with additional cosmological data sets. The WMAP data alone, and in combination, continue to be remarkably well fit by a six-parameter LCDM model. When WMAP data are combined with measurements of the high-l CMB anisotropy, the BAO scale, and the Hubble constant, the densities, Omegabh2, Omegach2, and Omega_L, are each determined to a precision of ~1.5%. The amplitude of the primordial spectrum is measured to within 3%, and there is now evidence for a tilt in the primordial spectrum at the 5sigma level, confirming the first detection of tilt based on the five-year WMAP data. At the end of the WMAP mission, the nine-year data decrease the allowable volume of the six-dimensional LCDM parameter space by a factor of 68,000 relative to pre-WMAP measurements. We investigate a number of data combinations and show that their LCDM parameter fits are consistent. New limits on deviations from the six-parameter model are presented, for example: the fractional contribution of tensor modes is limited to r<0.13 (95% CL); the spatial curvature parameter is limited to -0.0027 (+0.0039/-0.0038); the summed mass of neutrinos is <0.44 eV (95% CL); and the number of relativistic species is found to be 3.84+/-0.40 when the full data are analyzed. The joint constraint on Neff and the primordial helium abundance agrees with the prediction of standard Big Bang nucleosynthesis. We compare recent PLANCK measurements of the Sunyaev-Zel'dovich effect with our seven-year measurements, and show their mutual agreement. Our analysis of the polarization pattern around temperature extrema is updated. This confirms a fundamental prediction of the standard cosmological model and provides a striking illustration of acoustic oscillations and adiabatic initial conditions in the early universe.

In our effort to complete the census of low-mass stars and brown dwarfs in the immediate Solar Neighborhood, we present spectra, photometry, proper motions, and distance estimates for forty-two low-mass star and brown dwarf candidates discovered by the Wide-field Infrared Survey Explorer (WISE). We also present additional follow-up information on twelve candidates selected using WISE data but previously published elsewhere. The new discoveries include fifteen M dwarfs, seventeen L dwarfs, five T dwarfs, and five objects of other type. Among these discoveries is a newly identified "unusually red L dwarf" (WISE J223527.07+451140.9), four peculiar L dwarfs whose spectra are most readily explained as unresolved L+T binary systems, and a T9 dwarf (WISE J124309.61+844547.8). We also show that the recently discovered red L dwarf WISEP J004701.06+680352.1 (Gizis et al. 2012) may be a low-gravity object and hence young and potentially low mass (< 25 MJup).

Combining the Wide-Field Infrared Survey Explorer All-Sky Release with the 2MASS Point Source Catalog, we detect a nearby, moderately metal-poor stellar debris stream spanning 24 degrees across the southern sky. The stream, which we designate Alpheus, is at an estimated distance of ~1.9 kpc. Its position, orientation, width, estimated metallicity, and to some extent its distance, are in approximate agreement with what one might expect of the leading tidal tail of the southern globular cluster NGC 288.

We report the discovery of 87 new T dwarfs uncovered with the Wide-field Infrared Survey Explorer (WISE) and three brown dwarfs with extremely red near-infrared colors that exhibit characteristics of both L and T dwarfs. Two of the new T dwarfs are likely binaries with L7+/-1 primaries and mid-type T secondaries. In addition, our follow-up program has confirmed 10 previously identified T dwarfs and four photometrically-selected L and T dwarf candidates in the literature. This sample, along with the previous WISE discoveries, triples the number of known brown dwarfs with spectral types later than T5. Using the WISE All-Sky Source Catalog we present updated color-color and color-type diagrams for all the WISE-discovered T and Y dwarfs. Near-infrared spectra of the new discoveries are presented, along with spectral classifications. To accommodate later T dwarfs we have modified the integrated flux method of determining spectral indices to instead use the median flux. Furthermore, a newly defined J-narrow index differentiates the early-type Y dwarfs from late-type T dwarfs based on the J-band continuum slope. The K/J indices for this expanded sample show that 32% of late-type T dwarfs have suppressed K-band flux and are blue relative to the spectral standards, while only 11% are redder than the standards. Comparison of the Y/J and K/J index to models suggests diverse atmospheric conditions and supports the possible re-emergence of clouds after the L/T transition. We also discuss peculiar brown dwarfs and candidates that were found not to be substellar, including two Young Stellar Objects and two Active Galactic Nuclei. The coolest WISE-discovered brown dwarfs are the closest of their type and will remain the only sample of their kind for many years to come.

The Wide-field Infrared Survey Explorer has revealed a T8.5 brown dwarf (WISE J111838.70+312537.9) that exhibits common proper motion with a solar-neighborhood (8 pc) quadruple star system - Xi Ursae Majoris. The angular separation is 8.5 arc-min, and the projected physical separation is about 4000 AU. The sub-solar metallicity and low chromospheric activity of Xi UMa A argue that the system has an age of at least 2 Gyr. The infrared luminosity and color of the brown dwarf suggests the mass of this companion ranges between 14 and 38 Jupiter masses for system ages of 2 and 8 Gyr respectively.

We have monitored the position of the cool Y dwarf WISEPA J182831.08+265037.8 using a combination of ground- and space-based telescopes and have determined its distance to be 11.2$_{-1.0}^{+1.3}$ pc. Its absolute H magnitude, M$_H=22.21^{+0.25}_{-0.22}$ mag, suggests a mass in the range 0.5-20 M$_{Jup}$ for ages of 0.1-10 Gyr with an effective temperature in the range 250-400 K. The broad range in mass is due primarily to the unknown age of the object. Since the high tangential velocity of the object, 51$\pm5$ km s$^{-1}$, is characteristic of an old disk population, a plausible age range of 2-4 Gyr leads to a mass range of 3-6 M$_{Jup}$ based on fits to the (highly uncertain) COND evolutionary models. The range in temperature is due to the fact that no single model adequately represents the 1-5 $\mu$m spectral energy distribution (SED) of the source, failing by factors of up to 5 at either the short or long wavelength portions of the spectral energy distribution. The appearance of this very cold object may be affected by non-equilibrium chemistry or low temperature condensates forming clouds, two atmospheric processes that are known to be important in brown dwarf atmospheres but have proven difficult to model. Finally, we argue that there would have to be a very steep upturn in the number density of late type Y-dwarfs to account for the putative population of objects suggested by by recent microlensing observations. Whether WISE 1828+2650 sits at the low mass end of the brown dwarf population or is the first example of a large number of "free-floating" planets is not yet known.

The Wide-field Infrared Survey Explorer (WISE) has completed its all-sky survey at 3.4-22 micron. We merge the WISE data with optical SDSS data and provide a phenomenological characterization of mid-IR, extragalactic sources. WISE is most sensitive at 3.4micron(W1) and least at 22micron(W4). The W1 band probes massive early-type galaxies out to z\gtrsim1. This is more distant than SDSS identified early-type galaxies, consistent with the fact that 28% of 3.4micron sources have faint or no r-band counterparts (r>22.2). In contrast, 92-95% of 12 and 22micron sources have SDSS optical counterparts with r<22.2. WISE 3.4micron detects 89.8% of the entire SDSS QSO catalog at SNR(W1)>7, but only 18.9% at 22micron with SNR(W4)>5. We show that WISE colors alone are effective in isolating stars (or local early-type galaxies), star-forming galaxies and strong AGN/QSOs at z<3. We highlight three major applications of WISE colors: (1) Selection of strong AGN/QSOs at z<3 using W1-W2>0.8 and W2<15.2 criteria, producing a census of this population. The surface density of these strong AGN/QSO candidates is 67.5+-0.14/deg^2. (2) Selection of dust-obscured, type-2 AGN/QSO candidates. We show that WISE W1-W2>0.8, W2<15.2 combined with r-W2>6 (Vega) colors can be used to identify type-2 AGN candidates. The fraction of these type-2 AGN candidates is 1/3rd of all WISE color-selected AGNs. (3) Selection of ULIRGs at z\sim2 with extremely red colors, r-W4>14 or well-detected 22micron sources lacking detections in the 3.4 and 4.6micron bands. The surface density of z~2 r-W4>14 ULIRGs is 0.9+-0.07/deg^2 at SNR(W4)>5 (flux(W4)>=2.5mJy), which is consistent with that inferred from smaller area Spitzer surveys. Optical spectroscopy of a small number of these high-redshift ULIRGs confirms our selection, and reveals a possible trend that optically fainter or r-W4 redder candidates are at higher redshifts.

We present astrometric measurements of eleven nearby ultracool brown dwarfs of spectral types Y and late-T, based on imaging observations from a variety of space-based and ground-based telescopes. These measurements have been used to estimate relative parallaxes and proper motions via maximum likelihood fitting of geometric model curves. To compensate for the modest statistical significance (<~ 7) of our parallax measurements we have employed a novel Bayesian procedure for distance estimation which makes use of an a priori distribution of tangential velocities, Vtan, assumed similar to that implied by previous observations of T dwarfs. Our estimated distances are therefore somewhat dependent on that assumption. Nevertheless, the results have yielded distances for five of our eight Y dwarfs and all three T dwarfs. Estimated distances in all cases are >~ 3 pc. In addition, we have obtained significant estimates of Vtan for two of the Y dwarfs; both are <100 km/s, consistent with membership in the thin disk population. Comparison of absolute magnitudes with model predictions as a function of color shows that the Y dwarfs are significantly redder in J-H than predicted by a cloud-free model

We measure the faint end slope of the galaxy luminosity function (LF) for cluster galaxies at 1<z<1.5 using Spitzer IRAC data. We investigate whether this slope, alpha, differs from that of the field LF at these redshifts, and with the cluster LF at low redshifts. The latter is of particular interest as low-luminosity galaxies are expected to undergo significant evolution. We use seven high-redshift spectroscopically confirmed galaxy clusters drawn from the IRAC Shallow Cluster Survey to measure the cluster galaxy LF down to depths of M* + 3 (3.6 microns) and M* + 2.5 (4.5 microns). The summed LF at our median cluster redshift (z=1.35) is well fit by a Schechter distribution with alpha[3.6] = -0.97 +/- 0.14 and alpha[4.5] = -0.91 +/- 0.28, consistent with a flat faint end slope and is in agreement with measurements of the field LF in similar bands at these redshifts. A comparison to alpha in low-redshift clusters finds no statistically significant evidence of evolution. Combined with past studies which show that M* is passively evolving out to z~1.3, this means that the shape of the cluster LF is largely in place by z~1.3. This suggests that the processes that govern the build up of the mass of low-mass cluster galaxies have no net effect on the faint end slope of the cluster LF at z<1.3.

The Wide-field Infrared Survey Explorer (WISE) mapped the entire sky at mid-infrared wavelengths 3.4, 4.6, 12 and 22 microns. The mission was primarily designed to extract point sources, leaving resolved and extended sources unexplored. We have begun a dedicated WISE Enhanced Resolution Galaxy Atlas (WERGA) project to fully characterize large, nearby galaxies and produce a legacy image atlas and source catalogue. Here we demonstrate the first results of the project for a sample of 17 galaxies, chosen to be of large angular size, diverse morphology, color, stellar mass and star formation. It includes many well-studied galaxies, such as M51, M81, M83, M87, M101, IC342. Photometry and surface brightness decomposition is carried out after special super-resolution processing, achieving spatial fidelity similar to that of Spitzer-IRAC. We present WISE, Spitzer and GALEX photometric and characterization measurements, combining the measurements to study the global properties. We derive star formation rates using the PAH-sensitive 12 micron (W3) fluxes, warm-dust sensitive 22 micron (W4) fluxes, and young massive-star sensitive UV fluxes. Stellar masses are estimated using the 3.4 micron (W1) and 4.6 micron (W2) measurements that trace the dominant stellar mass content. We highlight and showcase the detailed results of M83, comparing the infrared results with the ATCA HI gas distribution and GALEX UV emission, tracing the evolution from gas to stars. In addition to the enhanced images, WISE all-sky coverage provides a tremendous advantage over Spitzer for building a complete nearby galaxy catalog, tracing both stellar mass and star formation histories. We discuss the construction of a complete mid-infrared catalog of galaxies and its complementary role to study the assembly and evolution of galaxies in the local universe.

Unresolved near-infrared background anisotropies are expected to have contributions from the earliest galaxies during reionization and faint, dwarf galaxies at intermediate redshifts. Previous measurements were unable to conclusively pinpoint the dominant origin because they did not sample spatial scales that were sufficiently large to distinguish between these two possibilities. Here we report a measurement of the anisotropy power spectrum from sub-arcminute to one degree angular scales and find the clustering amplitude to be larger than the model predictions involving the two existing explanations. As the shot-noise level of the power spectrum is consistent with that expected from faint galaxies, a new source population on the sky is not necessary to explain the observations. A physical mechanism that increases the clustering amplitude, however, is needed. Motivated by recent results related to the extended stellar light profile in dark matter halos, we consider the possibility that the fluctuations originate from diffuse intrahalo stars of all galaxies. We find that the measured power spectrum can be explained by an intrahalo light fraction of 0.07 to 0.2 % relative to the total luminosity in dark matter halos of masses log(M/M_Sun) ~ 9 to 12 at redshifts of ~ 1 to 4.

We present spectroscopic confirmation of a z=0.99 galaxy cluster discovered using data from the Wide-field Infrared Survey Explorer (WISE). This is the first z~1 cluster candidate from the Massive Distant Clusters of WISE Survey (MaDCoWS) to be confirmed. It was selected as an overdensity of probable z>~1 sources using a combination of WISE and SDSS-DR8 photometric catalogs. Deeper follow-up imaging data from Subaru and WIYN reveal the cluster to be a rich system of galaxies, and multi-object spectroscopic observations from Keck confirm five cluster members at z=0.99. The detection and confirmation of this cluster represents a first step towards constructing a uniformly-selected sample of distant, high-mass galaxy clusters over the full extragalactic sky using WISE data.

Enhancements to the science data processing pipeline of NASA's Wide-field Infrared Explorer (WISE) mission, collectively known as NEOWISE, resulted in the detection of $>$158,000 minor planets in four infrared wavelengths during the fully cryogenic portion of the mission. Following the depletion of its cryogen, NASA's Planetary Science Directorate funded a four month extension to complete the survey of the inner edge of the Main Asteroid Belt and to detect and discover near-Earth objects (NEOs). This extended survey phase, known as the NEOWISE Post-Cryogenic Survey, resulted in the detection of $\sim$6500 large Main Belt asteroids and 88 NEOs in its 3.4 and 4.6 $\mu$m channels. During the Post-Cryogenic Survey, NEOWISE discovered and detected a number of asteroids co-orbital with the Earth and Mars, including the first known Earth Trojan. We present preliminary thermal fits for these and other NEOs detected during the 3-Band Cryogenic and Post-Cryogenic Surveys.

We have used methane imaging techniques to identify the near-infrared counterpart of the bright WISE source WISEJ163940.83-684738.6. The large proper motion of this source (around 3.0arcsec/yr) has moved it, since its original WISE identification, very close to a much brighter background star -- it currently lies within 1.5" of the J=14.90+-0.04 star 2MASS16394085-6847446. Observations in good seeing conditions using methane sensitive filters in the near-infrared J-band with the FourStar instrument on the Magellan 6.5m Baade telescope, however, have enabled us to detect a near-infrared counterpart. We have defined a photometric system for use with the FourStar J2 and J3 filters, and this photometry indicates strong methane absorption, which unequivocally identifies it as the source of the WISE flux. Using these imaging observations we were then able to steer this object down the slit of the FIRE spectrograph on a night of 0.6" seeing, and so obtain near-infrared spectroscopy confirming a Y0-Y0.5 spectral type. This is in line with the object's near-infrared-to-WISE J3--W2 colour. Preliminary astrometry using both WISE and FourStar data indicates a distance of 5.0+-0.5pc and a substantial tangential velocity of 73+-8km/s. WISEJ163940.83-684738.6 is the brightest confirmed Y dwarf in the WISE W2 passband and its distance measurement places it amongst the lowest luminosity sources detected to date.

On January 18-19 and June 28-29 of 2010, the Wide-field Infrared Survey Explorer (WISE) spacecraft imaged the Rosetta mission target, comet 67P/Churyumov-Gerasimenko. We present a preliminary analysis of the images, which provide a characterization of the dust environment at heliocentric distances similar to those planned for the initial spacecraft encounter, but on the outbound leg of its orbit rather than the inbound. Broad-band photometry yields low levels of CO2 production at a comet heliocentric distance of 3.32 AU and no detectable production at 4.18 AU. We find that at these heliocentric distances, large dust grains with mean grain diameters on the order of a millimeter or greater dominate the coma and evolve to populate the tail. This is further supported by broad-band photometry centered on the nucleus, which yield an estimated differential dust particle size distribution with a power law relation that is considerably shallower than average. We set a 3-sigma upper limit constraint on the albedo of the large-grain dust at <= 0.12. Our best estimate of the nucleus radius (1.82 +/- 0.20 km) and albedo (0.04 +/- 0.01) are in agreement with measurements previously reported in the literature.

We present Spitzer 3.6 and 4.5 $\mu$m photometry and positions for a sample of 1510 brown dwarf candidates identified by the WISE all-sky survey. Of these, 166 have been spectroscopically classified as objects with spectral types M(1), L(7), T(146), and Y(12); Sixteen other objects are non-(sub)stellar in nature. The remainder are most likely distant L and T dwarfs lacking spectroscopic verification, other Y dwarf candidates still awaiting follow-up, and assorted other objects whose Spitzer photometry reveals them to be background sources. We present a catalog of Spitzer photometry for all astrophysical sources identified in these fields and use this catalog to identify 7 fainter (4.5 $\mu$m $\sim$ 17.0 mag) brown dwarf candidates, which are possibly wide-field companions to the original WISE sources. To test this hypothesis, we use a sample of 919 Spitzer observations around WISE-selected high-redshift hyper-luminous infrared galaxy (HyLIRG) candidates. For this control sample we find another 6 brown dwarf candidates, suggesting that the 7 companion candidates are not physically associated. In fact, only one of these 7 Spitzer brown dwarf candidates has a photometric distance estimate consistent with being a companion to the WISE brown dwarf candidate. Other than this there is no evidence for any widely separated ($>$ 20 AU) ultra-cool binaries. As an adjunct to this paper, we make available a source catalog of $\sim$ 7.33 $\times 10^5$ objects detected in all of these Spitzer follow-up fields for use by the astronomical community. The complete catalog includes the Spitzer 3.6 and 4.5 $\mu$m photometry, along with positionally matched $B$ and $R$ photometry from USNO-B; $J$, $H$, and $K_s$ photometry from 2MASS; and $W1$, $W2$, $W3$, and $W4$ photometry from the WISE all-sky catalog.

We report the discovery by the Wide-field Infrared Survey Explorer of the z = 2.452 source WISE J181417.29+341224.9, the first hyperluminous source found in the WISE survey. WISE 1814+3412 is also the prototype for an all-sky sample of ~1000 extremely luminous "W1W2-dropouts" (sources faint or undetected by WISE at 3.4 and 4.6 microns and well detected at 12 or 22 microns). The WISE data and a 350 micron detection give a minimum bolometric luminosity of 3.7 x 10^13 Lsun, with ~10^14 Lsun plausible. Followup images reveal four nearby sources: a QSO and two Lyman Break Galaxies (LBGs) at z = 2.45, and an M dwarf star. The brighter LBG dominates the bolometric emission. Gravitational lensing is unlikely given the source locations and their different spectra and colors. The dominant LBG spectrum indicates a star formation rate ~300 Msun/yr, accounting for < 10% of the bolometric luminosity. Strong 22 micron emission relative to 350 microns implies that warm dust contributes significantly to the luminosity, while cooler dust normally associated with starbursts is constrained by an upper limit at 1.1 mm. Radio emission is ~10x above the far-infrared/radio correlation, indicating an active galactic nucleus is present. An obscured AGN combined with starburst and evolved stellar components can account for the observations. If the black hole mass follows the local M_BH-bulge mass relation, the implied Eddington ratio is >~4. WISE 1814+3412 may be a heavily obscured object where the peak AGN activity occurred prior to the peak era of star formation.

We have used the Caltech Submillimeter Observatory (CSO) to follow-up a sample of WISE-selected, hyperluminous galaxies, so called W1W2-dropout galaxies. This is a rare (~ 1000 all-sky) population of galaxies at high redshift (peaks at z=2-3), that are faint or undetected by WISE at 3.4 and 4.6 um, yet are clearly detected at 12 and 22 um. The optical spectra of most of these galaxies show significant AGN activity. We observed 14 high-redshift (z > 1.7) W1W2-dropout galaxies with SHARC-II at 350 to 850 um, with 9 detections; and observed 18 with Bolocam at 1.1 mm, with five detections. Warm Spitzer follow-up of 25 targets at 3.6 and 4.5 um, as well as optical spectra of 12 targets are also presented in the paper. Combining WISE data with observations from warm Spitzer and CSO, we constructed their mid-IR to millimeter spectral energy distributions (SEDs). These SEDs have a consistent shape, showing significantly higher mid-IR to submm ratios than other galaxy templates, suggesting a hotter dust temperature. We estimate their dust temperatures to be 60-120 K using a single-temperature model. Their infrared luminosities are well over 10^{13} Lsun. These SEDs are not well fitted with existing galaxy templates, suggesting they are a new population with very high luminosity and hot dust. They are likely among the most luminous galaxies in the Universe. We argue that they are extreme cases of luminous, hot dust-obscured galaxies (DOGs), possibly representing a short evolutionary phase during galaxy merging and evolution. A better understanding of their long-wavelength properties needs ALMA as well as Herschel data.

In December 2010, NASA created a Science Definition Team (SDT) for WFIRST, the Wide Field Infra-Red Survey Telescope, recommended by the Astro 2010 Decadal Survey as the highest priority for a large space mission. The SDT was chartered to work with the WFIRST Project Office at GSFC and the Program Office at JPL to produce a Design Reference Mission (DRM) for WFIRST. Part of the original charge was to produce an interim design reference mission by mid-2011. That document was delivered to NASA and widely circulated within the astronomical community. In late 2011 the Astrophysics Division augmented its original charge, asking for two design reference missions. The first of these, DRM1, was to be a finalized version of the interim DRM, reducing overall mission costs where possible. The second of these, DRM2, was to identify and eliminate capabilities that overlapped with those of NASA's James Webb Space Telescope (henceforth JWST), ESA's Euclid mission, and the NSF's ground-based Large Synoptic Survey Telescope (henceforth LSST), and again to reduce overall mission cost, while staying faithful to NWNH. This report presents both DRM1 and DRM2.

After eight months of continuous observations, the Wide-field Infrared Survey Explorer (WISE) mapped the entire sky at 3.4 {\mu}m, 4.6 {\mu}m, 12 {\mu}m and 22 {\mu}m. We have begun a dedicated WISE High Resolution Galaxy Atlas (WHRGA) project to fully characterize large, nearby galaxies and produce a legacy image atlas and source catalogue. Here we summarize the deconvolution technique used to significantly improve the spatial resolution of WISE imaging, specifically designed to study the internal anatomy of nearby galaxies. As a case study, we present results for the galaxy NGC 1566, comparing the WISE super-resolution image processing to that of Spitzer, GALEX and ground-based imaging. The is the first paper in a two part series; results for a much larger sample of nearby galaxies is presented in the second paper.

We report a new technique to select 1.6<z<4.6 dusty Lyman-alpha emitters (LAEs), over a third of which are `blobs' (LABs) with emission extended on scales of 30-100kpc. Combining data from the NASA Wide-field Infrared Survey Explorer (WISE) mission with optical spectroscopy from the W.M. Keck telescope, we present a color criteria that yields a 78% success rate in identifying rare, dusty LAEs of which at least 37% are LABs. The objects have a surface density of only ~0.1 per square degree, making them rare enough that they have been largely missed in narrow surveys. We measured spectroscopic redshifts for 92 of these WISE-selected, typically radio-quiet galaxies and find that the LAEs (LABs) have a median redshift of 2.3 (2.5). The WISE photometry coupled with data from Herschel reveals that these galaxies have extreme far-infrared luminosities (L_IR>10^{13-14}L_sun) and warm colors, typically larger than submillimeter-selected galaxies (SMGs) and dust-obscured galaxies (DOGs). These traits are commonly associated with the dust being energized by intense AGN activity. We hypothesize that the combination of spatially extended Lyman-alpha, large amounts of warm IR-luminous dust, and rarity (implying a short-lived phase) can be explained if the galaxies are undergoing strong `feedback' transforming them from an extreme dusty starburst to a QSO.

We present the preliminary results of an analysis of the sub-populations within the near-Earth asteroids, including the Atens, Apollos, Amors, and those that are considered potentially hazardous using data from the Wide-field Infrared Survey Explorer (WISE). In order to extrapolate the sample of objects detected by WISE to the greater population, we determined the survey biases for asteroids detected by the project's automated moving object processing system (known as NEOWISE) as a function of diameter, visible albedo, and orbital elements. Using this technique, we are able to place constraints on the number of potentially hazardous asteroids (PHAs) larger than 100 m and find that there are $\sim4700\pm1450$ such objects. As expected, the Atens, Apollos, and Amors are revealed by WISE to have somewhat different albedo distributions, with the Atens being brighter than the Amors. The cumulative size distributions of the various near-Earth object (NEO) subgroups vary slightly between 100 m and 1 km. A comparison of the observed orbital elements of the various sub-populations of the NEOs with the current best model is shown.

We present the discovery of another seven Y dwarfs from the Wide-field Infrared Survey Explorer (WISE). Using these objects, as well as the first six WISE Y dwarf discoveries from Cushing et al., we further explore the transition between spectral types T and Y. We find that the T/Y boundary roughly coincides with the spot where the J-H colors of brown dwarfs, as predicted by models, turn back to the red. Moreover, we use preliminary trigonometric parallax measurements to show that the T/Y boundary may also correspond to the point at which the absolute H (1.6 um) and W2 (4.6 um) magnitudes plummet. We use these discoveries and their preliminary distances to place them in the larger context of the Solar Neighborhood. We present a table that updates the entire stellar and substellar constituency within 8 parsecs of the Sun, and we show that the current census has hydrogen-burning stars outnumbering brown dwarfs by roughly a factor of six. This factor will decrease with time as more brown dwarfs are identified within this volume, but unless there is a vast reservoir of cold brown dwarfs invisible to WISE, the final space density of brown dwarfs is still expected to fall well below that of stars. We also use these new Y dwarf discoveries, along with newly discovered T dwarfs from WISE, to investigate the field substellar mass function. We find that the overall space density of late-T and early-Y dwarfs matches that from simulations describing the mass function as a power law with slope -0.5 < alpha < 0.0; however, a power-law may provide a poor fit to the observed object counts as a function of spectral type because there are tantalizing hints that the number of brown dwarfs continues to rise from late-T to early-Y. More detailed monitoring and characterization of these Y dwarfs, along with dedicated searches aimed at identifying more examples, are certainly required.

The Wide-field Infrared Survey Explorer (WISE) is an extremely capable and efficient black hole finder. We present a simple mid-infrared color criterion, W1-W2 \geq 0.8 (i.e., [3.4]-[4.6] \geq 0.8, Vega), which identifies 61.9 \pm 5.4 AGN candidates per deg2 to a depth of W2 = 15.0. This implies a much larger census of luminous AGN than found by typical wide-area surveys, attributable to the fact that mid-infrared selection identifies both unobscured (type 1) and obscured (type 2) AGN. Optical and soft X-ray surveys alone are highly biased towards only unobscured AGN, while this simple WISE selection likely identifies even heavily obscured, Compton-thick AGN. Using deep, public data in the COSMOS field, we explore the properties of WISE-selected AGN candidates. At the mid-infrared depth considered, 160 uJy at 4.6 microns, this simple criterion identifies 78% of Spitzer mid-infrared AGN candidates according to the criteria of Stern et al. (2005) and the reliability is 95%. We explore the demographics, multiwavelength properties and redshift distribution of WISE-selected AGN candidates in the COSMOS field.

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

We present a reanalysis of the relationship between asteroid albedo and polarization properties using the albedos derived from the Wide-field Infrared Survey Explorer. We find that the function that best describes this relation is a three-dimensional linear fit in the space of log(albedo)-log(polarization slope)-log(minimum polarization). When projected to two dimensions the parameters of the fit are consistent with those found in previous work. We also define p* as the quantity of maximal polarization variation when compared with albedo and present the best fitting albedo-p* relation. Some asteroid taxonomic types stand out in this three-dimensional space, notably the E, B, and M Tholen types, while others cluster in clumps coincident with the S- and C-complex bodies. We note that both low albedo and small (D<30 km) asteroids are under-represented in the polarimetric sample, and we encourage future polarimetric surveys to focus on these bodies.

We report results based on mid-infrared photometry of 5 active main belt objects (AMBOs) detected by the Wide-field Infrared Survey Explorer (WISE) spacecraft. Four of these bodies, P/2010 R2 (La Sagra), 133P/Elst-Pizarro, (596) Scheila, and 176P/LINEAR, showed no signs of activity at the time of the observations, allowing the WISE detections to place firm constraints on their diameters and albedos. Geometric albedos were in the range of a few percent, and on the order of other measured comet nuclei. P/2010 A2 was observed on April 2-3, 2010, three months after its peak activity. Photometry of the coma at 12 and 22 {\mu}m combined with ground-based visible-wavelength measurements provides constraints on the dust particle mass distribution (PMD), dlogn/dlogm, yielding power-law slope values of {\alpha} = -0.5 +/- 0.1. This PMD is considerably more shallow than that found for other comets, in particular inbound particle fluence during the Stardust encounter of comet 81P/Wild 2. It is similar to the PMD seen for 9P/Tempel 1 in the immediate aftermath of the Deep Impact experiment. Upper limits for CO2 & CO production are also provided for each AMBO and compared with revised production numbers for WISE observations of 103P/Hartley 2.

We report on the results of an optical spectroscopic survey at high Galactic latitude (b\geq30{\deg}) of a sample of WISE-selected targets, grouped by WISE W1 ({\lambda}_eff = 3.4 \mum) flux, which we use to characterize the sources WISE detected. We observed 762 targets in 10 disjoint fields centered on ultra-luminous infrared galaxy (ULIRG) candidates using the DEIMOS spectrograph on Keck II. We find 0.30 \pm 0.02 galaxies arcmin^-2 with a median redshift of z=0.33 \pm 0.01 for the sample with W1 \geq 120 \muJy. The foreground stellar densities in our survey range from 0.23 \pm 0.07 arcmin^(-2) to 1.1 \pm 0.1 arcmin^(-2) for the same sample. We obtained spectra that produced science grade redshifts for \geq 90% of our targets for sources with W1 flux \geq 120 \muJy that also had i-band flux \geq 18 \muJy. We used for targeting very preliminary data reductions available to the team in August of 2010. Our results therefore present a conservative estimate of what is possible to achieve using WISE's Preliminary Data Release for the study of field galaxies.

We have combined the NEOWISE and Sloan Digital Sky Survey data to study the albedos of 24,353 asteroids with candidate taxonomic classifications derived using Sloan photometry. We find a wide range of moderate to high albedos for candidate S-type asteroids that are analogous to the S-complex defined by previous spectrophotometrically-based taxonomic systems. The candidate C-type asteroids, while generally very dark, have a tail of higher albedos that overlaps the S types. The albedo distribution for asteroids with a photometrically derived Q classification is extremely similar to those of the S types. Asteroids with similar colors to (4) Vesta have higher albedos than the S types, and most have orbital elements similar to known Vesta family members. Finally, we show that the relative reflectance at 3.4 and 4.6 $\mu$m is higher for D-type asteroids and suggest that their red visible and near-infrared spectral slope extends out to these wavelengths. Understanding the relationship between size, albedo, and taxonomic classification is complicated by the fact that the objects with classifications were selected from the visible/near-infrared Sloan Moving Object Catalog, which is biased against fainter asteroids, including those with lower albedos.

We present the preliminary analysis of 1023 known asteroids in the Hilda region of the Solar System observed by the NEOWISE component of the Wide-field Infrared Survey Explorer (WISE). The sizes of the Hildas observed range from $\sim 3 - 200$km. We find no size - albedo dependency as reported by other projects. The albedos of our sample are low, with a weighted mean value $p_V = 0.055\pm0.018$, for all sizes sampled by the NEOWISE survey. We observed a significant fraction of the objects in the two known collisional families in the Hilda population. It is found that the Hilda collisional family is brighter, with weighted mean albedo of $p_V = 0.061\pm0.011$, than the general population and dominated by D-type asteroids, while the Schubart collisional family is darker, with weighted mean albedo of ($p_V = 0.039\pm0.013$). Using the reflected sunlight in the two shortest WISE bandpasses we are able to derive a method for taxonomic classification of $\sim 10%$ of the Hildas detected in the NEOWISE survey. For the Hildas with diameter larger than 30km there are $67^{+7}_{-15}%$ D-type asteroids and $26^{+17}_{-5}%$ C-/P-type asteroids (with the majority of these being P-types).

We present the preliminary analysis of over 1739 known and 349 candidate Jovian Trojans observed by the NEOWISE component of the Wide-field Infrared Survey Explorer (WISE). With this survey the available diameters, albedos and beaming parameters for the Jovian Trojans have been increased by more than an order of magnitude compared to previous surveys. We find that the Jovian Trojan population is very homogenous for sizes larger than $\sim10$km (close to the detection limit of WISE for these objects). The observed sample consists almost exclusively of low albedo objects, having a mean albedo value of $0.07\pm0.03$. The beaming parameter was also derived for a large fraction of the observed sample, and it is also very homogenous with an observed mean value of $0.88\pm0.13$. Preliminary debiasing of the survey shows our observed sample is consistent with the leading cloud containing more objects than the trailing cloud. We estimate the fraction to be N(leading)/N(trailing) $\sim 1.4 \pm 0.2$, lower than the $1.6 \pm 0.1$ value derived by others.

With the NEOWISE portion of the \emph{Wide-field Infrared Survey Explorer} (WISE) project, we have carried out a highly uniform survey of the near-Earth object (NEO) population at thermal infrared wavelengths ranging from 3 to 22 $\mu$m, allowing us to refine estimates of their numbers, sizes, and albedos. The NEOWISE survey detected NEOs the same way whether they were previously known or not, subject to the availability of ground-based follow-up observations, resulting in the discovery of more than 130 new NEOs. The survey's uniformity in sensitivity, observing cadence, and image quality have permitted extrapolation of the 428 near-Earth asteroids (NEAs) detected by NEOWISE during the fully cryogenic portion of the WISE mission to the larger population. We find that there are 981$\pm$19 NEAs larger than 1 km and 20,500$\pm$3000 NEAs larger than 100 m. We show that the Spaceguard goal of detecting 90% of all 1 km NEAs has been met, and that the cumulative size distribution is best represented by a broken power law with a slope of 1.32$\pm$0.14 below 1.5 km. This power law slope produces $\sim13,200\pm$1,900 NEAs with $D>$140 m. Although previous studies predict another break in the cumulative size distribution below $D\sim$50-100 m, resulting in an increase in the number of NEOs in this size range and smaller, we did not detect enough objects to comment on this increase. The overall number for the NEA population between 100-1000 m are lower than previous estimates. The numbers of near-Earth comets will be the subject of future work.

The NEOWISE dataset offers the opportunity to study the variations in albedo for asteroid classification schemes based on visible and near-infrared observations for a large sample of minor planets. We have determined the albedos for nearly 1900 asteroids classified by the Tholen, Bus and Bus-DeMeo taxonomic classification schemes. We find that the S-complex spans a broad range of bright albedos, partially overlapping the low albedo C-complex at small sizes. As expected, the X-complex covers a wide range of albedos. The multi-wavelength infrared coverage provided by NEOWISE allows determination of the reflectivity at 3.4 and 4.6 $\mu$m relative to the visible albedo. The direct computation of the reflectivity at 3.4 and 4.6 $\mu$m enables a new means of comparing the various taxonomic classes. Although C, B, D and T asteroids all have similarly low visible albedos, the D and T types can be distinguished from the C and B types by examining their relative reflectance at 3.4 and 4.6 $\mu$m. All of the albedo distributions are strongly affected by selection biases against small, low albedo objects, as all objects selected for taxonomic classification were chosen according to their visible light brightness. Due to these strong selection biases, we are unable to determine whether or not there are correlations between size, albedo and space weathering. We argue that the current set of classified asteroids makes any such correlations difficult to verify. A sample of taxonomically classified asteroids drawn without significant albedo bias is needed in order to perform such an analysis.

We present initial results from the Wide-field Infrared Survey Explorer (WISE), a four-band all-sky thermal infrared survey that produces data well suited to measuring the physical properties of asteroids, and the NEOWISE enhancement to the WISE mission allowing for detailed study of Solar system objects. Using a NEATM thermal model fitting routine we compute diameters for over 100,000 Main Belt asteroids from their IR thermal flux, with errors better than 10%. We then incorporate literature values of visible measurements (in the form of the H absolute magnitude) to determine albedos. Using these data we investigate the albedo and diameter distributions of the Main Belt. As observed previously, we find a change in the average albedo when comparing the inner, middle, and outer portions of the Main Belt. We also confirm that the albedo distribution of each region is strongly bimodal. We observe groupings of objects with similar albedos in regions of the Main Belt associated with dynamical breakup families. Asteroid families typically show a characteristic albedo for all members, but there are notable exceptions to this. This paper is the first look at the Main Belt asteroids in the WISE data, and only represents the preliminary, observed raw size and albedo distributions for the populations considered. These distributions are subject to survey biases inherent to the NEOWISE dataset and cannot yet be interpreted as describing the true populations; the debiased size and albedo distributions will be the subject of the next paper in this series.

We present ground-based spectroscopic verification of six Y dwarfs (see Cushing et al), eighty-nine T dwarfs, eight L dwarfs, and one M dwarf identified by the Wide-field Infrared Survey Explorer (WISE). Eighty of these are cold brown dwarfs with spectral types greater than or equal to T6, six of which have been announced earlier in Mainzer et al and Burgasser et al. We present color-color and color-type diagrams showing the locus of M, L, T, and Y dwarfs in WISE color space. Near-infrared classifications as late as early Y are presented and objects with peculiar spectra are discussed. After deriving an absolute WISE 4.6 um (W2) magnitude vs. spectral type relation, we estimate spectrophotometric distances to our discoveries. We also use available astrometric measurements to provide preliminary trigonometric parallaxes to four our discoveries, which have types of L9 pec (red), T8, T9, and Y0; all of these lie within 10 pc of the Sun. The Y0 dwarf, WISE 1541-2250, is the closest at 2.8 (+1.3,-0.6) pc; if this 2.8 pc value persists after continued monitoring, WISE 1541-2250 will become the seventh closest stellar system to the Sun. Another ten objects, with types between T6 and >Y0, have spectrophotometric distance estimates also placing them within 10 pc. The closest of these, the T6 dwarf WISE 1506+7027, is believed to fall at a distance of roughly 4.9 pc. WISE multi-epoch positions supplemented with positional info primarily from Spitzer/IRAC allow us to calculate proper motions and tangential velocities for roughly one half of the new discoveries. This work represents the first step by WISE to complete a full-sky, volume-limited census of late-T and Y dwarfs. Using early results from this census, we present preliminary, lower limits to the space density of these objects and discuss constraints on both the functional form of the mass function and the low-mass limit of star formation.

We present the discovery of seven ultracool brown dwarfs identified with the Wide-field Infrared Survey Explorer (WISE). Near-infrared spectroscopy reveals deep absorption bands of H_2O and CH_4 that indicate all seven of the brown dwarfs have spectral types later than UGPS J072227.51-054031.2, the latest type T dwarf currently known. The spectrum of WISEP J182831.08+265037.8 is distinct in that the heights of the J- and H-band peaks are approximately equal in units of f_lambda, so we identify it as the archetypal member of the Y spectral class. The spectra of at least two of the other brown dwarfs exhibit absorption on the blue wing of the H-band peak that we tentatively ascribe to NH_3. These spectral morphological changes provide a clear transition between the T dwarfs and the Y dwarfs. In order to produce a smooth near-infrared spectral sequence across the T/Y dwarf transition, we have reclassified UGPS J0722-0540 as the T9 spectral standard and tentatively assign WISEP J173835.52+273258.9 as the Y0 spectral standard. In total, six of the seven new brown dwarfs are classified as Y dwarfs: four are classified as Y0, one is classified as Y0 (pec?), and WISEP J1828+2650 is classified as >Y0. We have also compared the spectra to the model atmospheres of Marley and Saumon and infer that the brown dwarfs have effective temperatures ranging from 300 K to 500 K, making them the coldest spectroscopically confirmed brown dwarfs known to date.

In December 2010, NASA created a Science Definition Team (SDT) for WFIRST, the Wide Field Infra-Red Survey Telescope, recommended by the Astro 2010 Decadal Survey as the highest priority for a large space mission. The SDT was chartered to work with the WFIRST Project Office at GSFC and the Program Office at JPL to produce a Design Reference Mission (DRM) for WFIRST. This paper describes an Interim DRM. The DRM will be completed in 2012.

We report results based on mid-infrared photometry of comet 103P/Hartley 2 taken during May 4-13, 2010 (when the comet was at a heliocentric distance of 2.3 AU, and an observer distance of 2.0 AU) by the Wide-field Infrared Survey Explorer (Wright et al. 2010). Photometry of the coma at 22 microns and data from the University of Hawaii 2.2-m telescope obtained on May 22, 2010 provide constraints on the dust particle size distribution, dlogn/dlogm, yielding power-law slope values of alpha = -0.97 +/- 0.10, steeper than that found for the inbound particle fluence during the Stardust encounter of comet 81P/Wild 2 (Green et al. 2004). The extracted nucleus signal at 12 microns is consistent with a body of average spherical radius of 0.6 +/- 0.2 km (one standard deviation), assuming a beaming parameter of 1.2. The 4.6 micron-band signal in excess of dust and nucleus reflected and thermal contributions may be attributed to carbon monoxide or carbon dioxide emission lines and provides limits and estimates of species production. Derived carbon dioxide coma production rates are 3.5(+/- 0.9) \times 10^24 molecules per second. Analyses of the trail signal present in the stacked image with an effective exposure time of 158.4 seconds yields optical-depth values near 9 x 10^-10 at a delta mean anomaly of 0.2 deg trailing the comet nucleus, in both 12 and 22 {\mu}m bands. A minimum chi-squared analysis of the dust trail position yields a beta-parameter value of 1.0 x10^-4, consistent with a derived mean trail-grain diameter of 1.1/{\rho} cm for grains of {\rho} g/cm^3 density. This leads to a total detected trail mass of at least 4 x 10^10 {\rho} kg.

With thermal infrared observations detected by the NEOWISE project, we have measured diameters for 1742 minor planets that were also observed by the Infrared Astronomical Satellite (IRAS). We have compared the diameters and albedo derived by applying a spherical thermal model to the objects detected by NEOWISE and find that they are in generally good agreement with the IRAS values. We have shown that diameters computed from NEOWISE data are often less systematically biased than those found with IRAS. This demonstrates that the NEOWISE dataset can provide accurate physical parameters for the >157,000 minor planets that were detected by NEOWISE.

We report two new low metallicity blue compact dwarf galaxies (BCDs), WISEP J080103.93+264053.9 (hereafter W0801+26) and WISEP J170233.53+180306.4 (hereafter W1702+18), discovered using the Wide-field Infrared Survey Explorer (WISE). We identified these two BCDs from their extremely red colors at mid-infrared wavelengths, and obtained follow-up optical spectroscopy using the Low Resolution Imaging Spectrometer on Keck I. The mid-infrared properties of these two sources are similar to the well studied, extremely low metallicity galaxy SBS 0335-052E. We determine metallicities of 12 + log(O/H) = 7.75 and 7.63 for W0801+26 and W1702+18, respectively, placing them amongst a very small group of very metal deficient galaxies (Z < 1/10 Zsun). Their > 300 Angstrom Hbeta equivalent widths, similar to SBS 0335-052E, imply the existence of young (< 5 Myr) star forming regions. We measure star formation rates of 2.6 and 10.9 Msun/yr for W0801+26 and W1702+18, respectively. These BCDs, showing recent star formation activity in extremely low metallicity environments, provide new laboratories for studying star formation in extreme conditions and are low-redshift analogs of the first generation of galaxies to form in the universe. Using the all-sky WISE survey, we discuss a new method to identify similar star forming, low metallicity BCDs.

The multiplicity properties of brown dwarfs are critical empirical constraints for formation theories, while multiples themselves provide unique opportunities to test evolutionary and atmospheric models and examine empirical trends. Studies using high-resolution imaging can not only uncover faint companions, but they can also be used to determine dynamical masses through long-term monitoring of binary systems. We have begun a search for the coolest brown dwarfs using preliminary processing of data from the Wide-field Infrared Survey Explorer (WISE) and have confirmed many of the candidates as late-type T dwarfs. In order to search for companions to these objects, we are conducting observations using the Laser Guide Star Adaptive Optics system on Keck II. Here we present the first results of that search, including a T5 binary with nearly equal mass components and a faint companion to a T8.5 dwarf with an estimated spectral type of T9.

With the Wide-field Infrared Survey Explorer (WISE; Wright et al. 2010), we have observed over 157,000 minor planets (Mainzer et al. 2011). Included in these are a number of near-Earth objects, Main Belt Asteroids, and irregular satellites which have well-measured physical properties via radar, occultation and in situ imaging. We have used these objects to validate models of thermal models using the WISE measurements, as well as the color corrections derived in Wright et al. (2010) for the four WISE bandpasses as a function of effective temperature. We have used 50 objects with diameters measured by radar, occultation or in situ imaging to characterize the systematic errors implicit in using the WISE data with a faceted spherical NEATM model to compute diameters and albedos. By using the previously measured diameters and H magnitudes with a spherical NEATM model, we compute the predicted fluxes after applying the color corrections given in Wright et al. (2010) in the WISE bands and compare them to the measured magnitudes. We find minimum systematic flux errors of 5-10%, yielding minimum relative diameter and albedo errors of ~10% and ~20%, respectively. Visible albedos for the objects are computed and compared to the albedos at 3.4 and 4.6 microns, which contain a mix of reflected sunlight and thermal emission for most asteroids. We derive a linear relationship between subsolar temperature and effective temperature, which allows the color corrections given in Wright et al. (2010) to be used for asteroids by computing only subsolar temperature instead of a faceted thermal model. The thermal models derived in this paper are not intended to supplant previous measurements made using radar or spacecraft imaging; rather, we have used them to characterize the errors that should be expected when computing diameters and albedos of WISE asteroids using a spherical NEATM model.

We present the discovery of five late-type T dwarfs identified with the Wide-field Infrared Survey Explorer (WISE). Low-resolution near-infrared spectroscopy obtained with the Magellan Folded-port InfraRed Echellette (FIRE) reveal strong water and methane absorption in all five sources, and spectral indices and comparison to spectral templates indicate classifications ranging from T5.5 to T8.5:. The spectrum of the latest-type source, WISE J1812+2721, is an excellent match to that of the T8.5 companion brown dwarf Wolf 940B. WISE-based spectrophotometric distance estimates place these T dwarfs at 12-13 pc from the Sun, assuming they are single. Preliminary fits of the spectral data to the atmosphere models of Saumon & Marley indicate effective temperatures ranging from 600 K to 930 K, both cloudy and cloud-free atmospheres, and a broad range of ages and masses. In particular, two sources show evidence of both low surface gravity and cloudy atmospheres, tentatively supporting a trend noted in other young brown dwarfs and exoplanets. In contrast, the high proper motion T dwarf WISE J2018-7423 exhibits a suppressed K-band peak and blue spectrophotometric J-K colors indicative of an old, massive brown dwarf; however, it lacks the broadened Y-band peak seen in metal-poor counterparts. These results illustrate the broad diversity of low-temperature brown dwarfs that will be uncovered with WISE.

The Wide-field Infrared Survey Explorer (WISE) has detected the close binary brown dwarf system CFBDSIR 1458+1013AB as WISEP J145829.35+101341.8 with a combined magnitude at 4.6 microns of W2 = 15.488+/-0.147. This allows a comparison with another "coolest" brown dwarf candidate WD 0806-661B that has been observed at 4.5 microns with [4.5] = 16.75+/-0.05. Here we use the WISE data to show that 1458+1013B is almost certainly warmer and more luminous than WD 0806-661B.

The \emph{Wide-field Infrared Survey Explorer} has surveyed the entire sky at four infrared wavelengths with greatly improved sensitivity and spatial resolution compared to its predecessors, the \emph{Infrared Astronomical Satellite} and the \emph{Cosmic Background Explorer}. NASA's Planetary Science Division has funded an enhancement to the \WISE\ data processing system called "NEOWISE" that allows detection and archiving of moving objects found in the \WISE\ data. NEOWISE has mined the \WISE\ images for a wide array of small bodies in our Solar System, including Near-Earth Objects (NEOs), Main Belt asteroids, comets, Trojans, and Centaurs. By the end of survey operations in February 2011, NEOWISE identified over 157,000 asteroids, including more than 500 NEOs and $\sim$120 comets. The NEOWISE dataset will enable a panoply of new scientific investigations.

(Abridged) A simple six-parameter LCDM model provides a successful fit to WMAP data, both when the data are analyzed alone and in combination with other cosmological data. Even so, it is appropriate to search for any hints of deviations from the now standard model of cosmology, which includes inflation, dark energy, dark matter, baryons, and neutrinos. The cosmological community has subjected the WMAP data to extensive and varied analyses. While there is widespread agreement as to the overall success of the six-parameter LCDM model, various "anomalies" have been reported relative to that model. In this paper we examine potential anomalies and present analyses and assessments of their significance. In most cases we find that claimed anomalies depend on posterior selection of some aspect or subset of the data. Compared with sky simulations based on the best fit model, one can select for low probability features of the WMAP data. Low probability features are expected, but it is not usually straightforward to determine whether any particular low probability feature is the result of the a posteriori selection or of non-standard cosmology. We examine in detail the properties of the power spectrum with respect to the LCDM model. We examine several potential or previously claimed anomalies in the sky maps and power spectra, including cold spots, low quadrupole power, quadropole-octupole alignment, hemispherical or dipole power asymmetry, and quadrupole power asymmetry. We conclude that there is no compelling evidence for deviations from the LCDM model, which is generally an acceptable statistical fit to WMAP and other cosmological data.

(Abridged) We present WMAP seven-year observations of bright sources which are often used as calibrators at microwave frequencies. Ten objects are studied in five frequency bands (23 - 94 GHz): the outer planets (Mars, Jupiter, Saturn, Uranus and Neptune) and five fixed celestial sources (Cas A, Tau A, Cyg A, 3C274 and 3C58). The seven-year analysis of Jupiter provides temperatures which are within 1-sigma of the previously published WMAP five-year values, with slightly tighter constraints on variability with orbital phase, and limits (but no detections) on linear polarization. Scaling factors are provided which, when multiplied by the Wright Mars thermal model predictions at 350 micron, reproduce WMAP seasonally averaged observations of Mars within ~2%. An empirical model is described which fits brightness variations of Saturn due to geometrical effects and can be used to predict the WMAP observations to within 3%. Seven-year mean temperatures for Uranus and Neptune are also tabulated. Uncertainties in Uranus temperatures are 3%-4% in the 41, 61 and 94 GHz bands; the smallest uncertainty for Neptune is ~8% for the 94 GHz band. Intriguingly, the spectrum of Uranus appears to show a dip at ~30 GHz of unidentified origin, although the feature is not of high statistical significance. Flux densities for the five selected fixed celestial sources are derived from the seven-year WMAP sky maps, and are tabulated for Stokes I, Q and U, along with polarization fraction and position angle. Fractional uncertainties for the Stokes I fluxes are typically 1% to 3%. Source variability over the seven-year baseline is also estimated. Significant secular decrease is seen for Cas A and Tau A: our results are consistent with a frequency independent decrease of about 0.53% per year for Cas A and 0.22% per year for Tau A.

[Abridged] We present updated estimates of Galactic foreground emission using seven years of WMAP data. Using the power spectrum of differences between multi-frequency template-cleaned maps, we find no evidence for foreground contamination outside of the updated (KQ85y7) foreground mask. We place a 15 microKelvin upper bound on rms foreground contamination in the cleaned maps used for cosmological analysis. We find no indication in the polarization data of an extra "haze" of hard synchrotron emission from energetic electrons near the Galactic center. We provide an updated map of the cosmic microwave background (CMB) using the internal linear combination (ILC) method, updated foreground masks, and updates to point source catalogs with 62 newly detected sources. Also new are tests of the Markov chain Monte Carlo (MCMC) foreground fitting procedure against systematics in the time-stream data, and tests against the observed beam asymmetry. Within a few degrees of the Galactic plane, WMAP total intensity data show a rapidly steepening spectrum from 20-40 GHz, which may be due to emission from spinning dust grains, steepening synchrotron, or other effects. Comparisons are made to a 1-degree 408 MHz map (Haslam et al.) and the 11-degree ARCADE 2 data (Singal et al.). We find that spinning dust or steepening synchrotron models fit the combination of WMAP and 408 MHz data equally well. ARCADE data appear inconsistent with the steepening synchrotron model, and consistent with the spinning dust model, though some discrepancies remain regarding the relative strength of spinning dust emission. More high-resolution data in the 10-40 GHz range would shed much light on these issues.

We report the discovery of a pair of infrared, axisymmetric rings in the planetary nebula NGC 1514 during the course of the WISE all-sky mid-infrared survey. Similar structures are seen at visible wavelengths in objects such as the "Engraved Hourglass Nebula" (MyCn 18) and the "Southern Crab Nebula" (Hen 2-104). However, in NGC 1514 we see only a single pair of rings and they are easily observed only in the mid-infrared. These rings are roughly 0.2 pc in diameter, are separated by 0.05 pc, and are dominated by dust emission with a characteristic temperature of 160 K. We compare the morphology and color of the rings to the other nebular structures seen at visible, far-infrared, and radio wavelengths, and close with a discussion of a physical model and formation scenario for NGC 1514.

We report the discovery of the first new ultra-cool brown dwarf found with the Wide-field Infrared Survey Explorer (WISE). The object's preliminary designation is WISEPC J045853.90+643451.9. Follow-up spectroscopy with the LUCIFER instrument on the Large Binocular Telescope indicates that it is a very late-type T dwarf with a spectral type approximately equal to T9. Fits to an IRTF/SpeX 0.8-2.5 micron spectrum to the model atmospheres of Marley and Saumon indicate an effective temperature of approximately 600 K as well as the presence of vertical mixing in its atmosphere. The new brown dwarf is easily detected by WISE, with a signal-to-noise ratio of ~36 at 4.6 microns. Current estimates place it at a distance of 6 to 10 pc. This object represents the first in what will likely be hundreds of nearby brown dwarfs found by WISE that will be suitable for follow up observations, including those with the James Webb Space Telescope. One of the two primary scientific goals of the WISE mission is to find the coolest, closest stars to our Sun; the discovery of this new brown dwarf proves that WISE is capable of fulfilling this objective.

The all sky surveys done by the Palomar Observatory Schmidt, the European Southern Observatory Schmidt, and the United Kingdom Schmidt, the InfraRed Astronomical Satellite and the 2 Micron All Sky Survey have proven to be extremely useful tools for astronomy with value that lasts for decades. The Wide-field Infrared Survey Explorer is mapping the whole sky following its launch on 14 December 2009. WISE began surveying the sky on 14 Jan 2010 and completed its first full coverage of the sky on July 17. The survey will continue to cover the sky a second time until the cryogen is exhausted (anticipated in November 2010). WISE is achieving 5 sigma point source sensitivities better than 0.08, 0.11, 1 and 6 mJy in unconfused regions on the ecliptic in bands centered at wavelengths of 3.4, 4.6, 12 and 22 microns. Sensitivity improves toward the ecliptic poles due to denser coverage and lower zodiacal background. The angular resolution is 6.1, 6.4, 6.5 and 12.0 arc-seconds at 3.4, 4.6, 12 and 22 microns, and the astrometric precision for high SNR sources is better than 0.15 arc-seconds.

(Abridged) The 7-year WMAP data and improved astrophysical data rigorously test the standard cosmological model and its extensions. By combining WMAP with the latest distance measurements from BAO and H0 measurement, we determine the parameters of the simplest LCDM model. The power-law index of the primordial power spectrum is n_s=0.968+-0.012, a measurement that excludes the scale-invariant spectrum by 99.5%CL. The other parameters are also improved from the 5-year results. Notable examples of improved parameters are the total mass of neutrinos, sum(m_nu)<0.58eV, and the effective number of neutrino species, N_eff=4.34+0.86-0.88. We detect the effect of primordial helium on the temperature power spectrum and provide a new test of big bang nucleosynthesis. We detect, and show on the map for the first time, the tangential and radial polarization patterns around hot and cold spots of temperature fluctuations, an important test of physical processes at z=1090 and the dominance of adiabatic scalar fluctuations. With the 7-year TB power spectrum, the limit on a rotation of the polarization plane due to potential parity-violating effects has improved to Delta(alpha)=-1.1+-1.4(stat)+-1.5(syst) degrees. We report significant detections of the SZ effect at the locations of known clusters of galaxies. The measured SZ signal agrees well with the expected signal from the X-ray data. However, it is a factor of 0.5 to 0.7 times the predictions from "universal profile" of Arnaud et al., analytical models, and hydrodynamical simulations. We find, for the first time in the SZ effect, a significant difference between the cooling-flow and non-cooling-flow clusters (or relaxed and non-relaxed clusters), which can explain some of the discrepancy. This lower amplitude is consistent with the lower-than-theoretically-expected SZ power spectrum recently measured by the South Pole Telescope collaboration.

We study a 24\,$\mu$m selected sample of 330 galaxies observed with the Infrared Spectrograph for the 5\,mJy Unbiased Spitzer Extragalactic Survey. We estimate accurate total infrared luminosities by combining mid-IR spectroscopy and mid-to-far infrared photometry, and by utilizing new empirical spectral templates from {\em Spitzer} data. The infrared luminosities of this sample range mostly from 10$^9$L$_\odot$ to $10^{13.5}$L$_\odot$, with 83% in the range 10$^{10}$L$_\odot$$<$L$_{\rm IR}$$<10^{12}$L$_\odot$. The redshifts range from 0.008 to 4.27, with a median of 0.144. The equivalent widths of the 6.2\,$\mu$m aromatic feature have a bimodal distribution. We use the 6.2\,$\mu$m PAH EW to classify our objects as SB-dominated (44%), SB-AGN composite (22%), and AGN-dominated (34%). The high EW objects (SB-dominated) tend to have steeper mid-IR to far-IR spectral slopes and lower L$_{\rm IR}$ and redshifts. The low EW objects (AGN-dominated) tend to have less steep spectral slopes and higher L$_{\rm IR}$ and redshifts. This dichotomy leads to a gross correlation between EW and slope, which does not hold within either group. AGN dominated sources tend to have lower log(L$_{\rm PAH 7.7\mu m}$/L$_{\rm PAH 11.3\mu m}$) ratios than star-forming galaxies, possibly due to preferential destruction of the smaller aromatics by the AGN. The log(L$_{\rm PAH 7.7\mu m}$/L$_{\rm PAH 11.3\mu m}$) ratios for star-forming galaxies are lower in our sample than the ratios measured from the nuclear spectra of nearby normal galaxies, most probably indicating a difference in the ionization state or grain size distribution between the nuclear regions and the entire galaxy. Finally, we provide a calibration relating the monochromatic 5.8, 8, 14 and 24um continuum or Aromatic Feature luminosity to L$_{\rm IR}$ for different types of objects.

[NIRSS is one of three concepts that contributed to the Wide-Field Infrared Survey Telescope (WFIRST) mission advocated by the Decadal Survey.] Operating beyond the reaches of the Earth's atmosphere, free of its limiting absorption and thermal background, the Near-Infrared Sky Surveyor (NIRSS) will deeply map the entire sky at near-infrared wavelengths, thereby enabling new and fundamental discoveries ranging from the identification of extrasolar planets to probing the reionization epoch by identifying thousands of quasars at z>10. NIRSS will directly address the NASA scientific objective of studying cosmic origins by using a 1.5-meter telescope to reach full-sky 0.2 uJy (25.6 mag AB) sensitivities in four passbands from 1 to 4 microns in a 4-yr mission. At the three shorter passbands (1 - 2.5 microns), the proposed depth is comparable to the deepest pencil-beam surveys done to date and is 3000 times more sensitive than the only previous all-sky near-infrared survey, 2MASS. At the longest passband (3.5 micron), which is not feasible from the ground, NIRSS will be 500 times more sensitive than WISE. NIRSS fills a pivotal gap in our knowledge of the celestial sphere, is a natural complement to WISE, and is well matched to the next generation of deep (0.1 uJy), wide-area (>2 pi ster), ground-based optical surveys (LSST and Pan-Starrs). With the high thermal backgrounds of ground-based infrared observations, a near-infrared full sky survey at sub-uJy sensitivity is only feasible from space.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

(Abridged) We present the angular power spectra derived from the 7-year maps and discuss the cosmological conclusions that can be inferred from WMAP data alone. The third acoustic peak in the TT spectrum is now well measured by WMAP. In the context of a flat LambdaCDM model, this improvement allows us to place tighter constraints on the matter density from WMAP data alone, and on the epoch of matter-radiation equality, The temperature-polarization (TE) spectrum is detected in the 7-year data with a significance of 20 sigma, compared to 13 sigma with the 5-year data. The low-l EE spectrum, a measure of the optical depth due to reionization, is detected at 5.5 sigma significance when averaged over l = 2-7. The BB spectrum, an important probe of gravitational waves from inflation, remains consistent with zero. The upper limit on tensor modes from polarization data alone is a factor of 2 lower with the 7-year data than it was using the 5-year data (Komatsu et al. 2010). We test the parameter recovery process for bias and find that the scalar spectral index, ns, is biased high, but only by 0.09 sigma, while the remaining parameters are biased by < 0.15 sigma. The improvement in the third peak measurement leads to tighter lower limits from WMAP on the number of relativistic degrees of freedom (e.g., neutrinos) in the early universe: Neff > 2.7 (95% CL). Also, using WMAP data alone, the primordial helium mass fraction is found to be YHe = 0.28+0.14-0.15, and with data from higher-resolution CMB experiments included, we now establish the existence of pre-stellar helium at > 3 sigma (Komatsu et al. 2010).

We use the multi-epoch, mid-infrared Spitzer Deep, Wide-Field Survey to investigate the variability of 474,179 objects in 8.1 deg^2 of the NDWFS Bootes field. We perform a Difference Image Analysis of the four available epochs between 2004 and 2008, focusing on the deeper 3.6 and 4.5 micron bands. We find that 1.1% of the studied sample meet our standard selection criteria for being classed as a variable source. We require that the 3.6 and 4.5 micron light-curves are strongly correlated (r>0.8) and that their joint variance exceeds that for all sources with the same magnitude by 2 sigma. We then examine the mid-IR colors of the variable sources and match them with X-ray sources from the XBootes survey, radio catalogs, 24 micron-selected AGN candidates, and spectroscopically identified AGNs from the AGN and Galaxy Evolution Survey (AGES). Based on their mid-IR colors, most of the variable sources are AGNs (76%), with smaller contributions from stars (11%), galaxies (6%), and unclassified objects. Most of the stellar, galaxy and unclassified sources are false positives. For our standard selection criteria, 11-12% of the mid-IR counterparts to X-ray sources, 24 micron-selected AGN candidates and spectroscopically identified AGNs show variability. Mid-IR AGN variability can be well described by a single power-law structure function with a power-law index of 0.5 at both 3.6 and 4.5 microns, and an amplitude of 0.1 mag on rest-frame time scales of 2 years. The variability amplitude is higher for shorter rest-frame wavelengths and lower luminosities. (Abridged)

We have identified a sample of cool field brown dwarf candidates using IRAC data from the Spitzer Deep, Wide-Field Survey (SDWFS). The candidates were selected from 400,000 SDWFS sources with [4.5] <= 18.5 mag and required to have [3.6]-[4.5] >= 1.5 and [4.5] - [8.0] <= 2.0 on the Vega system. The first color requirement selects objects redder than all but a handful of presently known brown dwarfs with spectral classes later than T7, while the second eliminates 14 probable reddened AGN. Optical detection of 4 of the remaining 18 sources implies they are likely also AGN, leaving 14 brown dwarf candidates. For two of the brightest candidates (SDWFS J143524.44+335334.6 and SDWFS J143222.82+323746.5), the spectral energy distributions including near-infrared detections suggest a spectral class of ~ T8. The proper motion is < 0.25 "/yr, consistent with expectations for a luminosity inferred distance of >70 pc. The reddest brown dwarf candidate (SDWFS J143356.62+351849.2) has [3.6] - [4.5]=2.24 and H - [4.5] > 5.7, redder than any published brown dwarf in these colors, and may be the first example of the elusive Y-dwarf spectral class. Models from Burrows et al. (2003) predict larger numbers of cool brown dwarfs should be found for a Chabrier (2003) mass function. Suppressing the model [4.5] flux by a factor of two, as indicated by previous work, brings the Burrows models and observations into reasonable agreement. The recently launched Wide-field Infrared Survey Explorer (WISE) will probe a volume ~40x larger and should find hundreds of brown dwarfs cooler than T7.

(Abridged) New full sky temperature and polarization maps based on seven years of data from WMAP are presented. The new results are consistent with previous results, but have improved due to reduced noise from the additional integration time, improved knowledge of the instrument performance, and improved data analysis procedures. The improvements are described in detail. The seven year data set is well fit by a minimal six-parameter flat Lambda-CDM model. The parameters for this model, using the WMAP data in conjunction with baryon acoustic oscillation data from the Sloan Digital Sky Survey and priors on H_0 from Hubble Space Telescope observations, are: Omega_bh^2 = 0.02260 +-0.00053, Omega_ch^2 = 0.1123 +-0.0035, Omega_Lambda = 0.728 +0.015 -0.016, n_s = 0.963 +-0.012, tau = 0.087 +-0.014 and sigma_8 = 0.809 +-0.024 (68 % CL uncertainties). The temperature power spectrum signal-to-noise ratio per multipole is greater that unity for multipoles < 919, allowing a robust measurement of the third acoustic peak. This measurement results in improved constraints on the matter density, Omega_mh^2 = 0.1334 +0.0056 -0.0055, and the epoch of matter- radiation equality, z_eq = 3196 +134 -133, using WMAP data alone. The new WMAP data, when combined with smaller angular scale microwave background anisotropy data, results in a 3 sigma detection of the abundance of primordial Helium, Y_He = 0.326 +-0.075.The power-law index of the primordial power spectrum is now determined to be n_s = 0.963 +-0.012, excluding the Harrison-Zel'dovich-Peebles spectrum by >3 sigma. These new WMAP measurements provide important tests of Big Bang cosmology.

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

Explorers have made breakthroughs in many fields of astrophysics. The science from both these missions contributed to three Nobel Prizes - Giacconi (2002), Mather, and Smoot (2006). Explorers have: marked the definitive beginning of precision cosmology, discovered that short gamma-ray bursts are caused by compact star mergers and have measured metalicity to redshifts z>6. NASA Explorers do cutting-edge science that cannot be done by facility-class instruments. The Explorer program provides a rapid response to changing science and technology, to enable cutting-edge science at moderate cost. Explorers also enable innovation, and engage & train scientists, managers and engineers, adding human capital to NASA and the nation. The astrophysics Explorer launch rate now being achieved is 1 per 3 years, and budget projections are in the $150M/year range for the next five years. A newly Vigorous Explorer Program should be created to: 1. Reach the long-stated goal of annual astrophysics launches; 2. Find additional launch options for Explorers and actively encourage cost savings in launchers and spacecraft, such as new commercial vehicles and innovative partnerships. 3. Mitigate risk via stronger technical development and sub-orbital programs, and through longer, more thorough, Phase A programs, potentially reducing the need for a 30% contingency; 4. Strive to protect the funding for missions that have reached Phase B, to prevent significant launch slips and cancellations, with a goal of 4 to 5 years from Phase B to launch; 5. Review the project management procedures and requirements to seek cost reductions, including the risk management strategy and the review and reporting process; 6. Review and possibly modify the cost caps for all Explorer classes to optimize scientific returns per dollar. [ABRIDGED]

We describe a sampling method to estimate the polarized CMB signal from observed maps of the sky. We use a Metropolis-within-Gibbs algorithm to estimate the polarized CMB map, containing Q and U Stokes parameters at each pixel, and its covariance matrix. These can be used as inputs for cosmological analyses. The polarized sky signal is parameterized as the sum of three components: CMB, synchrotron emission, and thermal dust emission. The polarized Galactic components are modeled with spatially varying power law spectral indices for the synchrotron, and a fixed power law for the dust, and their component maps are estimated as by-products. We apply the method to simulated low resolution maps with pixels of side 7.2 degrees, using diagonal and full noise realizations drawn from the WMAP noise matrices. The CMB maps are recovered with goodness of fit consistent with errors. Computing the likelihood of the E-mode power in the maps as a function of optical depth to reionization, tau, for fixed temperature anisotropy power, we recover tau=0.091+-0.019 for a simulation with input tau=0.1, and mean tau=0.098 averaged over 10 simulations. A `null' simulation with no polarized CMB signal has maximum likelihood consistent with tau=0. The method is applied to the five-year WMAP data, using the K, Ka, Q and V channels. We find tau=0.090+-0.019, compared to tau=0.086+-0.016 from the template-cleaned maps used in the primary WMAP analysis. The synchrotron spectral index, beta, averaged over high signal-to-noise pixels with standard deviation sigma(beta)<0.25, but excluding ~6% of the sky masked in the Galactic plane, is -3.03+-0.04. This estimate does not vary significantly with Galactic latitude, although includes an informative prior.

The Spitzer Deep, Wide-Field Survey (SDWFS) is a four-epoch infrared survey of ten square degrees in the Bootes field of the NOAO Deep Wide-Field Survey using the IRAC instrument on the Spitzer Space Telescope. SDWFS, a Cycle four Spitzer Legacy project, occupies a unique position in the area-depth survey space defined by other Spitzer surveys. The four epochs that make up SDWFS permit -- for the first time -- the selection of infrared-variable and high proper motion objects over a wide field on timescales of years. Because of its large survey volume, SDWFS is sensitive to galaxies out to z~3 with relatively little impact from cosmic variance for all but the richest systems. The SDWFS datasets will thus be especially useful for characterizing galaxy evolution beyond z~1.5. This paper explains the SDWFS observing strategy and data processing, presents the SDWFS mosaics and source catalogs, and discusses some early scientific findings. The publicly-released, full-depth catalogs contain 6.78, 5.23, 1.20, and 0.96 x 10e5 distinct sources detected to the average 5-sigma, 4" diameter, aperture-corrected limits of 19.77, 18.83, 16.50, and 15.82 Vega mag at 3.6, 4.5, 5.8, and 8.0 micron, respectively. The SDWFS number counts and color-color distribution are consistent with other, earlier Spitzer surveys. At the 6 min integration time of the SDWFS IRAC imaging, more than 50% of isolated FIRST radio sources and more than 80% of on-axis XBootes sources are detected out to 8.0 micron. Finally, we present the four highest proper motion IRAC-selected sources identified from the multi-epoch imaging, two of which are likely field brown dwarfs of mid-T spectral class.

The answers to fundamental science questions in astrophysics, ranging from the history of the expansion of the universe to the sizes of nearby stars, hinge on our ability to make precise measurements of diverse astronomical objects. As our knowledge of the underlying physics of objects improves along with advances in detectors and instrumentation, the limits on our capability to extract science from measurements is set, not by our lack of understanding of the nature of these objects, but rather by the most mundane of all issues: the precision with which we can calibrate observations in physical units. We stress the need for a program to improve upon and expand the current networks of spectrophotometrically calibrated stars to provide precise calibration with an accuracy of equal to and better than 1% in the ultraviolet, visible and near-infrared portions of the spectrum, with excellent sky coverage and large dynamic range.

We present the results of an extragalactic point source search using the five-year WMAP 41, 61 and 94 GHz (Q-, V- and W-band) temperature maps. This work is an extension of our point source search in the WMAP maps applying a CMB-free technique. An internal linear combination (ILC) map has been formed from the three-band maps, with the weights chosen to remove the CMB anisotropy signal as well as to favor a selection of flat-spectrum sources. We find 381 sources at the > 5 sigma level outside the WMAP point source detection mask in the ILC map, among which 89 are "new" (i.e., not present in the WMAP catalog). Source fluxes have been calculated and corrected for the Eddington bias. We have solidly identified 367 (96.3%) of our sources. The 1 sigma positional uncertainty is estimated to be 2'. The 14 unidentified sources could be either extended radio structure or obscured by Galactic emission. We have also applied the same detection process on simulated maps and found 364+/-21 detections on average. The recovered source distribution N(>S) agrees well with the simulation input, which proves the reliability of this method.

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

We present a new estimate of foreground emission in the WMAP data, using a Markov chain Monte Carlo (MCMC) method. The new technique delivers maps of each foreground component for a variety of foreground models, error estimates of the uncertainty of each foreground component, and provides an overall goodness-of-fit measurement. The resulting foreground maps are in broad agreement with those from previous techniques used both within the collaboration and by other authors. We find that for WMAP data, a simple model with power-law synchrotron, free-free, and thermal dust components fits 90% of the sky with a reduced chi-squared of 1.14. However, the model does not work well inside the Galactic plane. The addition of either synchrotron steepening or a modified spinning dust model improves the fit. This component may account for up to 14% of the total flux at Ka-band (33 GHz). We find no evidence for foreground contamination of the CMB temperature map in the 85% of the sky used for cosmological analysis.

This paper focuses on cosmological constraints derived from analysis of WMAP data alone. A simple LCDM cosmological model fits the five-year WMAP temperature and polarization data. The basic parameters of the model are consistent with the three-year data and now better constrained: Omega_b h^2 = 0.02273+-0.00062, Omega_c h^2 = 0.1099+-0.0062, Omega_L = 0.742+-0.030, n_s = 0.963+0.014- 0.015, tau = 0.087+-0.017, sigma_8 = 0.796+-0.036. With five years of polarization data, we have measured the optical depth to reionization, tau>0, at 5 sigma significance. The redshift of an instantaneous reionization is constrained to be z_reion = 11.0+-1.4 with 68% confidence. This excludes a sudden reionization of the universe at z=6 at more than 3.5 sigma significance, suggesting that reionization was an extended process. Using two methods for polarized foreground cleaning we get consistent estimates for the optical depth, indicating an error due to foreground treatment of tau~0.01. This cosmological model also fits small-scale CMB data, and a range of astronomical data measuring the expansion rate and clustering of matter in the universe. We find evidence for the first time in the CMB power spectrum for a non-zero cosmic neutrino background, or a background of relativistic species, with the standard three light neutrino species preferred over the best-fit LCDM model with N_eff=0 at >99.5% confidence, and N_eff > 2.3 (95% CL) when varied. The five-year WMAP data improve the upper limit on the tensor-to-scalar ratio to r < 0.43 (95% CL), for power-law models. With longer integration we find no evidence for a running spectral index, with dn_s/dlnk = -0.037+-0.028.

We present new full-sky temperature and polarization maps in five frequency bands from 23 to 94 GHz, based on data from the first five years of the WMAP sky survey. The five-year maps incorporate several improvements in data processing made possible by the additional years of data and by a more complete analysis of the instrument calibration and in-flight beam response. We present several new tests for systematic errors in the polarization data and conclude that Ka band data (33 GHz) is suitable for use in cosmological analysis, after foreground cleaning. This significantly reduces the overall polarization uncertainty. With the 5 year WMAP data, we detect no convincing deviations from the minimal 6-parameter LCDM model: a flat universe dominated by a cosmological constant, with adiabatic and nearly scale-invariant Gaussian fluctuations. Using WMAP data combined with measurements of Type Ia supernovae and Baryon Acoustic Oscillations, we find (68% CL uncertainties): Omega_bh^2 = 0.02267 \pm 0.00059, Omega_ch^2 = 0.1131 \pm 0.0034, Omega_Lambda = 0.726 \pm 0.015, n_s = 0.960 \pm 0.013, tau = 0.084 \pm 0.016, and Delta_R^2 = (2.445 \pm 0.096) x 10^-9. From these we derive: sigma_8 = 0.812 \pm 0.026, H_0 = 70.5 \pm 1.3 km/s/Mpc, z_{reion} = 10.9 \pm 1.4, and t_0 = 13.72 \pm 0.12 Gyr. The new limit on the tensor-to-scalar ratio is r < 0.22 (95% CL). We obtain tight, simultaneous limits on the (constant) dark energy equation of state and spatial curvature: -0.14 < 1+w < 0.12 and -0.0179 < Omega_k < 0.0081 (both 95% CL). The number of relativistic degrees of freedom (e.g. neutrinos) is found to be N_{eff} = 4.4 \pm 1.5, consistent with the standard value of 3.04. Models with N_{eff} = 0 are disfavored at >99.5% confidence.

(Abridged) The WMAP 5-year data strongly limit deviations from the minimal LCDM model. We constrain the physics of inflation via Gaussianity, adiabaticity, the power spectrum shape, gravitational waves, and spatial curvature. We also constrain the properties of dark energy, parity-violation, and neutrinos. We detect no convincing deviations from the minimal model. The parameters of the LCDM model, derived from WMAP combined with the distance measurements from the Type Ia supernovae (SN) and the Baryon Acoustic Oscillations (BAO), are: Omega_b=0.0456+-0.0015, Omega_c=0.228+-0.013, Omega_Lambda=0.726+-0.015, H_0=70.5+-1.3 km/s/Mpc, n_s=0.960+-0.013, tau=0.084+-0.016, and sigma_8=0.812+-0.026. With WMAP+BAO+SN, we find the tensor-to-scalar ratio r<0.22 (95% CL), and n_s>1 is disfavored regardless of r. We obtain tight, simultaneous limits on the (constant) equation of state of dark energy and curvature. We provide a set of "WMAP distance priors," to test a variety of dark energy models. We test a time-dependent w with a present value constrained as -0.33<1+w_0<0.21 (95% CL). Temperature and matter fluctuations obey the adiabatic relation to within 8.9% and 2.1% for the axion and curvaton-type dark matter, respectively. The TE and EB spectra constrain cosmic parity-violation. We find the limit on the total mass of neutrinos, sum(m_nu)<0.67 eV (95% CL), which is free from the uncertainty in the normalization of the large-scale structure data. The effective number of neutrino species is constrained as N_{eff} = 4.4+-1.5 (68%), consistent with the standard value of 3.04. Finally, limits on primordial non-Gaussianity are -9<f_{NL}^{local}<111 and -151<f_{NL}^{equil}<253 (95% CL) for the local and equilateral models, respectively.

We present the list of point sources found in the WMAP 5-year maps. The technique used in the first-year and three-year analysis now finds 390 point sources, and the five-year source catalog is complete for regions of the sky away from the galactic plane to a 2 Jy limit, with SNR > 4.7 in all bands in the least covered parts of the sky. The noise at high frequencies is still mainly radiometer noise, but at low frequencies the CMB anisotropy is the largest uncertainty. A separate search of CMB-free V-W maps finds 99 sources of which all but one can be identified with known radio sources. The sources seen by WMAP are not strongly polarized. Many of the WMAP sources show significant variability from year to year, with more than a 2:1 range between the minimum and maximum fluxes.

We present the temperature and polarization angular power spectra of the cosmic microwave background (CMB) derived from the first 5 years of WMAP data. The 5-year temperature (TT) spectrum is cosmic variance limited up to multipole l=530, and individual l-modes have S/N>1 for l<920. The best fitting six-parameter LambdaCDM model has a reduced chi^2 for l=33-1000 of chi^2/nu=1.06, with a probability to exceed of 9.3%. There is now significantly improved data near the third peak which leads to improved cosmological constraints. The temperature-polarization correlation (TE) is seen with high significance. After accounting for foreground emission, the low-l reionization feature in the EE power spectrum is preferred by \Delta\chi^2=19.6 for optical depth tau=0.089 by the EE data alone, and is now largely cosmic variance limited for l=2-6. There is no evidence for cosmic signal in the BB, TB, or EB spectra after accounting for foreground emission. We find that, when averaged over l=2-6, l(l+1)C^{BB}_l/2\pi < 0.15 uK^2 (95% CL).

Cosmology and other scientific results from the WMAP mission require an accurate knowledge of the beam patterns in flight. While the degree of beam knowledge for the WMAP one-year and three-year results was unprecedented for a CMB experiment, we have significantly improved the beam determination as part of the five-year data release. Physical optics fits are done on both the A and the B sides for the first time. The cutoff scale of the fitted distortions on the primary mirror is reduced by a factor of ~2 from previous analyses. These changes enable an improvement in the hybridization of Jupiter data with beam models, which is optimized with respect to error in the main beam solid angle. An increase in main-beam solid angle of ~1% is found for the V2 and W1-W4 differencing assemblies. Although the five-year results are statistically consistent with previous ones, the errors in the five-year beam transfer functions are reduced by a factor of ~2 as compared to the three-year analysis. We present radiometry of the planet Jupiter as a test of the beam consistency and as a calibration standard; for an individual differencing assembly, errors in the measured disk temperature are ~0.5%.

We observed three regions of the sky with Spitzer in which the Cosmic InfraRed Background (CIRB) has been determined at 3.5 microns using the method of subtracting 2MASS stellar fluxes from zodiacal light subtracted DIRBE maps. For each of these regions we have obtained 270 seconds of integration time per pixel with IRAC on Spitzer over the central square degree. We present galaxy counts in each of these approximately 1 square degree IRAC surveys. Along with deep galaxy counts in the Extended Groth Strip and GOODS North, we are able to compare the galactic contribution to the CIRB with the "DIRBE minus 2MASS'' determined L-band CIRB. Using the profile-fit photometry package GIM2D, we find a substantially larger flux contribution to the CIRB than that determined using aperture photometry. We have also made the first rigorous analysis of the uncertainties in determining the CIRB via galaxy counts in Spitzer images using a Monte Carlo Markov Chain simulation of our data analysis. Using a simple broken power law model for galaxy counts as a function of magnitude we find a most probable contribution to the CIRB from galaxies at 3.6 microns of 10.8^{+2.1}_{-1.1} kJy/sr. Even with this restricted model, however, we find that galaxy counting does not strongly constrain the CIRB from above. We are able to find solutions in which the CIRB runs away to large intensities without the need for an additional diffuse source.

The Extended Groth Strip (EGS) is one of the premier fields for extragalactic deep surveys. Deep observations of the EGS with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope cover an area of 0.38 square degrees to a 50% completeness limit of 1.5 uJy at 3.6 um. The catalog comprises 57434 objects detected at 3.6 um, with 84%, 28%, and 24% also detected at 4.5, 5.8, and 8.0 um. Number counts are consistent with results from other Spitzer surveys. Color distributions show that the EGS IRAC sources comprise a mixture of populations: low-redshift star-forming galaxies, quiescent galaxies dominated by stellar emission at a range of redshifts, and high redshift galaxies and AGN.

Space-based experiments today can uniquely address important questions related to the fundamental laws of Nature. In particular, high-accuracy physics experiments in space can test relativistic gravity and probe the physics beyond the Standard Model; they can perform direct detection of gravitational waves and are naturally suited for precision investigations in cosmology and astroparticle physics. In addition, atomic physics has recently shown substantial progress in the development of optical clocks and atom interferometers. If placed in space, these instruments could turn into powerful high-resolution quantum sensors greatly benefiting fundamental physics. We discuss the current status of space-based research in fundamental physics, its discovery potential, and its importance for modern science. We offer a set of recommendations to be considered by the upcoming National Academy of Sciences' Decadal Survey in Astronomy and Astrophysics. In our opinion, the Decadal Survey should include space-based research in fundamental physics as one of its focus areas. We recommend establishing an Astronomy and Astrophysics Advisory Committee's interagency ``Fundamental Physics Task Force'' to assess the status of both ground- and space-based efforts in the field, to identify the most important objectives, and to suggest the best ways to organize the work of several federal agencies involved. We also recommend establishing a new NASA-led interagency program in fundamental physics that will consolidate new technologies, prepare key instruments for future space missions, and build a strong scientific and engineering community. Our goal is to expand NASA's science objectives in space by including ``laboratory research in fundamental physics'' as an element in agency's ongoing space research efforts.

We report the results of an extragalactic point source search using the 61 and 94 GHz (V- and W-band) temperature maps from the Wilkinson Microwave Anisotropy Probe (WMAP). Applying a method that cancels the ``noise'' due to the CMB anisotropy signal, we find in the $|b| > 10\degr$ region 31 sources in the first-year maps and 64 sources in the three-year co-added maps, at a $5\sigma$ level. The 1$\sigma$ position uncertainties are 1.6' and 1.4' each. The increased detections and improved positional accuracy are expected from the higher signal-to-noise ratio of WMAP three-year data. All sources detected in the first-year maps are repeatedly detected in the three-year maps, which is a strong proof of the consistency and reliability of this method. Among all the detections, 21 are new, i.e. not in the WMAP three-year point source catalog. We associate all but two of them with known objects. The two unidentified sources are likely to be variable or extended as observations through VLA, CARMA and ATCA all show non-detection at the nominal locations. We derive the source count distribution at WMAP V-band by combining our verified detections with sources from the WMAP three-year catalog. Assuming the effect of source clustering is negligible, the contribution to the power spectrum from faint sources below 0.75 Jy is estimated to be $(2.4\pm0.8) \times 10^{-3} \mu K^2$ sr for V-band, which implies a source correction amplitude $A = 0.012\pm0.004 \mu K^2$ sr.

By providing sensitive sub-arcsecond images and integral field spectroscopy in the 25 - 400 micron wavelength range, a far-IR interferometer will revolutionize our understanding of planetary system formation, reveal otherwise-undetectable planets through the disk perturbations they induce, and spectroscopically probe the atmospheres of extrasolar giant planets in orbits typical of most of the planets in our solar system. The technical challenges associated with interferometry in the far-IR are greatly relaxed relative to those encountered at shorter wavelengths or when starlight nulling is required. A structurally connected far-IR interferometer with a maximum baseline length of 36 m can resolve the interesting spatial structures in nascent and developed exoplanetary systems and measure exozodiacal emission at a sensitivity level critical to TPF-I mission planning. The Space Infrared Interferometric Telescope was recommended in the Community Plan for Far-IR/Submillimeter Space Astronomy, studied as a Probe-class mission, and estimated to cost 800M dollars. The scientific communities in Europe, Japan, and Canada have also demonstrated a keen interest in far-IR interferometry through mission planning workshops and technology research, suggesting the possibility of an international collaborative effort.

With the release of the 2MASS All-Sky Point Source Catalog, stellar fluxes from 2MASS are used to remove the contribution due to Galactic stars from the intensity measured by DIRBE in 40 new regions in the North and South Galactic polar caps. After subtracting the interplanetary and Galactic foregrounds, a consistent residual intensity of 14.69 +/- 4.49 kJy/sr at 2.2 microns is found. Allowing for a constant calibration factor between the DIRBE 3.5 microns and the 2MASS 2.2 microns fluxes, a similar analysis leaves a residual intensity of 15.62 +/- 3.34 kJy/sr at 3.5 microns. The intercepts of the DIRBE minus 2MASS correlation at 1.25 microns show more scatter and are a smaller fraction of the foreground, leading to a still weak limit on the CIRB of 8.88 +/- 6.26 kJy/sr (1 sigma).

We present a full-sky model of polarized Galactic microwave emission based on three years of observations by the Wilkinson Microwave Anisotropy Probe (WMAP) at frequencies from 23 to 94 GHz. The model compares maps of the Stokes Q and U components from each of the 5 WMAP frequency bands in order to separate synchrotron from dust emission, taking into account the spatial and frequency dependence of the synchrotron and dust components. This simple two-component model of the interstellar medium accounts for at least 97% of the polarized emission in the WMAP maps of the microwave sky. Synchrotron emission dominates the polarized foregrounds at frequencies below 50 GHz, and is comparable to the dust contribution at 65 GHz. The spectral index of the synchrotron component, derived solely from polarization data, is -3.2 averaged over the full sky, with a modestly flatter index on the Galactic plane. The synchrotron emission has mean polarization fraction 2--4% in the Galactic plane and rising to over 20% at high latitude, with prominent features such as the North Galactic Spur more polarized than the diffuse component. Thermal dust emission has polarization fraction 1% near the Galactic center, rising to 6% at the anti-center. Diffuse emission from high-latitude dust is also polarized with mean fractional polarization 0.036 +/- 0.011.

The luminosity distance vs. redshift law is now measured using supernovae and gamma ray bursts, and the angular size distance is measured at the surface of last scattering by the CMB and at z = 0.35 by baryon acoustic oscillations. In this paper this data is fit to models for the equation of state with w = -1, w = const, and w(z) = w_0+w_a(1-a). The last model is poorly constrained by the distance data, leading to unphysical solutions where the dark energy dominates at early times unless the large scale structure and acoustic scale constraints are modified to allow for early time dark energy effects. A flat LambdaCDM model is consistent with all the data.

A cratered asteroid acts somewhat like a retroflector, sending light and infrared radiation back toward the Sun, while thermal inertia in a rotating asteroid causes the infrared radiation to peak over the ``afternoon'' part. In this paper a rotating, cratered asteroid model is described, and used to generate infrared fluxes which are then interpreted using the Near Earth Asteroid Thermal Model (NEATM). Even though the rotating, cratered model depends on three parameters not available to the NEATM (the dimensionless thermal inertia parameter and pole orientation), the NEATM gives diameter estimates that are accurate to 10 percent RMS for phase angles less than 60 degrees. For larger phase angles, such as back-lit asteroids, the infrared flux depends more strongly on these unknown parameters, so the diameter errors are larger. These results are still true for the non-spherical shapes typical of small Near Earth objects.

The predicted infrared brightness temperature of Mars using the 1976 model of Wright is tabulated here for the period 1983 to 2103. This model was developed for far-infrared calibration, and is still being used for JCMT calibration.

A simple cosmological model with only six parameters (matter density, Omega_m h^2, baryon density, Omega_b h^2, Hubble Constant, H_0, amplitude of fluctuations, sigma_8, optical depth, tau, and a slope for the scalar perturbation spectrum, n_s) fits not only the three year WMAP temperature and polarization data, but also small scale CMB data, light element abundances, large-scale structure observations, and the supernova luminosity/distance relationship. Using WMAP data only, the best fit values for cosmological parameters for the power-law flat LCDM model are (Omega_m h^2, Omega_b h^2, h, n_s, tau, sigma_8) = 0.1277+0.0080-0.0079, 0.02229+-0.00073, 0.732+0.031-0.032, 0.958+-0.016, 0.089+-0.030, 0.761+0.049-0.048). The three year data dramatically shrink the allowed volume in this six dimensional parameter space. Assuming that the primordial fluctuations are adiabatic with a power law spectrum, the WMAP data_alone_ require dark matter, and favor a spectral index that is significantly less than the Harrison-Zel'dovich-Peebles scale-invariant spectrum (n_s=1, r=0). Models that suppress large-scale power through a running spectral index or a large-scale cut-off in the power spectrum are a better fit to the WMAP and small scale CMB data than the power-law LCDM model; however, the improvement in the fit to the WMAP data is only Delta chi^2 = 3 for 1 extra degree of freedom. The combination of WMAP and other astronomical data yields significant constraints on the geometry of the universe, the equation of state of the dark energy, the gravitational wave energy density, and neutrino properties. Consistent with the predictions of simple inflationary theories, we detect no significant deviations from Gaussianity in the CMB maps.

The WMAP satellite has completed 3 years of observations of the cosmic microwave background radiation. The 3-year data products include several sets of full sky maps of the Stokes I, Q and U parameters in 5 frequency bands, spanning 23 to 94 GHz, and supporting items, such as beam window functions and noise covariance matrices. The processing used to produce the current sky maps and supporting products represents a significant advancement over the first year analysis, and is described herein. Improvements to the pointing reconstruction, radiometer gain modeling, window function determination and radiometer spectral noise parametrization are presented. A detailed description of the updated data processing that produces maximum likelihood sky map estimates is presented, along with the methods used to produce reduced resolution maps and corresponding noise covariance matrices. Finally two methods used to evaluate the noise of the full resolution sky maps are presented along with several representative year-to-year null tests, demonstrating that sky maps produced from data from different observational epochs are consistent.

We present new full-sky temperature maps in five frequency bands from 23 to 94 GHz, based on the first three years of the WMAP sky survey. The new maps, which are consistent with the first-year maps and more sensitive, incorporate improvements in data processing made possible by the additional years of data and by a more complete analysis of the polarization signal. These include refinements in the gain calibration and beam response models. We employ two forms of multi-frequency analysis to separate astrophysical foreground signals from the CMB, each of which improves on our first-year analyses. First, we form an improved 'Internal Linear Combination' map, based solely on WMAP data, by adding a bias correction step and by quantifying residual uncertainties in the resulting map. Second, we fit and subtract new spatial templates that trace Galactic emission; in particular, we now use low-frequency WMAP data to trace synchrotron emission. The WMAP point source catalog is updated to include 115 new sources. We derive the angular power spectrum of the temperature anisotropy using a hybrid approach that combines a maximum likelihood estimate at low l (large angular scales) with a quadratic cross-power estimate for l>30. Our best estimate of the CMB power spectrum is derived by averaging cross-power spectra from 153 statistically independent channel pairs. The combined spectrum is cosmic variance limited to l=400, and the signal-to-noise ratio per l-mode exceeds unity up to l=850. The first two acoustic peaks are seen at l=220.8 +- 0.7 and l=530.9 +- 3.8, respectively, while the first two troughs are seen at l=412.4 +- 1.9 and l=675.1 +- 11.1, respectively. The rise to the third peak is unambiguous; when the WMAP data are combined with higher resolution CMB measurements, the existence of a third acoustic peak is well established.

The Wilkinson Microwave Anisotropy Probe WMAP has mapped the entire sky in five frequency bands between 23 and 94 GHz with polarization sensitive radiometers. We present three-year full-sky maps of the polarization and analyze them for foreground emission and cosmological implications. These observations open up a new window for understanding the universe. WMAP observes significant levels of polarized foreground emission due to both Galactic synchrotron radiation and thermal dust emission. The least contaminated channel is at 61 GHz. Informed by a model of the Galactic foreground emission, we subtract the foreground emission from the maps. In the foreground corrected maps, for l=2-6, we detect l(l+1) C_l^{EE} / (2 pi) = 0.086 +-0.029 microkelvin^2. This is interpreted as the result of rescattering of the CMB by free electrons released during reionization and corresponds to an optical depth of tau = 0.10 +- 0.03. We see no evidence for B-modes, limiting them to l(l+1) C_l^{BB} / (2 pi) = -0.04 +- 0.03 microkelvin^2. We find that the limit from the polarization signals alone is r<2.2 (95% CL) corresponding to a limit on the cosmic density of gravitational waves of Omega_{GW}h^2 < 5 times 10^{-12}. From the full WMAP analysis, we find r<0.55 (95% CL) corresponding to a limit of Omega_{GW}h^2 < 10^{-12} (95% CL).

We describe the angular power spectrum of unresolved 3.6 micron IR light in Spitzer GOODS fields. The amplitude of the anisotropy spectrum decreases with decreasing flux threshold to which resolved sources are removed from images. When all pixels brighter than a Vega magnitude of 24.6 are removed, the amplitude of the power spectrum at arcminute angular scales can be described with an extra component of z>8 sources with a IRB contribution around 0.4 nW m^-2 sr-1. The shape of the power spectrum, however, is more consistent with that expected for unresolved, faint galaxies at lower redshifts with Vega magnitudes fainter than 23 with a total 3.6 micron intensity between 0.1 to 0.8 nW m^-2 sr^-1. We confirm this assumption by showing that large-scale power decreases rapidly when the unresolved clustering spectrum is measured from a processed HDF-N IRAC image where locations of faint ACS sources with no IR counterparts were also masked. Based on resolved counts and unresolved fluctuations, we find that, at most, about 7.0 nW m^-2 sr^-1 can be ascribed to galaxies.

We describe the angular power spectrum of resolved sources at 3.6 microns (L-band) in Spitzer imaging data of the GOODS HDF-N, the GOODS CDF-S, and the NDWFS Bootes field in several source magnitude bins. We also measure angular power spectra of resolved sources in the Bootes field at K_S and J-bands using ground-based IR imaging data. In the three bands, J, K_S, and L, we detect the clustering of galaxies on top of the shot-noise power spectrum at multipoles between ell ~ 10^2 and 10^5. The angular power spectra range from the large, linear scales to small, non-linear scales of galaxy clustering, and in some magnitude ranges, show departure from a power-law clustering spectrum. We consider a halo model to describe clustering measurements and to establish the halo occup ation number parameters of IR bright galaxies at redshifts around one. We also extend our clustering results and completeness-corrected faint source number counts in GOODS fields to understand the underlying nature of unresolved sources responsible for IR background (IRB) anisotropies that were detected in deep Spitzer images. While these unresolved fluctuations were measured at sub-arcminute angular scales, if a high-redshift diffuse component associated with first galaxies exists in the IRB, then it's clustering properties are best studied with shallow, wide-field images that allow a measurement of the clustering spectrum from a few degrees to arcminute angular scales.

A cosmology calculator that computes times and distances as a function of redshift for user-defined cosmological parameters is available on the World Wide Web. This note gives the formulae used by the cosmology calculator and discusses some of its implementation. A version of the calculator that allows one to specify the equation of state parameter w and w' and neutrino masses, and a version for converting the light travel times usually given in the popular press into redshifts are also available.

We discuss color selection of rare objects in a wide-field, multiband survey spanning from the optical to the mid-infrared. Simple color criteria simultaneously identify and distinguish two of the most sought after astrophysical sources: the coolest brown dwarfs and the most distant quasars. We present spectroscopically-confirmed examples of each class identified in the IRAC Shallow Survey of the Bootes field of the NOAO Deep Wide-Field Survey. ISS J142950.9+333012 is a T4.5 brown dwarf at a distance of approximately 42 pc, and ISS J142738.5+331242 is a radio-loud quasar at redshift z=6.12. Our selection criteria identify a total of four candidates over 8 square degrees of the Bootes field. The other two candidates are both confirmed 5.5<z<6 quasars, previously reported by Cool et al. (2006). We discuss the implications of these discoveries and conclude that there are excellent prospects for extending such searches to cooler brown dwarfs and higher redshift quasars.

In this the first of a series of Letters, we present a description of the panchromatic data sets that have been acquired in the Extended Groth Strip region of the sky. Our survey, the All-wavelength Extended Groth strip International Survey (AEGIS), is intended to study the physical properties and evolutionary processes of galaxies at z ~ 1. It includes the following deep, wide-field imaging data sets: Chandra/ACIS X-ray (0.5 - 10 keV), GALEX ultraviolet (1200 - 2500 Angstrom), CFHT/MegaCam Legacy Survey optical (3600 - 9000 Angstroms), CFHT/CFH12K optical (4500 - 9000 Angstroms), Hubble Space Telescope/ACS optical (4400 - 8500 Angstroms), Palomar/WIRC near-infrared (1.2 - 2.2 microns), Spitzer/IRAC mid-infrared (3.6 - 8.0 microns), Spitzer/MIPS far-infrared (24 - 70 microns), and VLA radio continuum (6 - 20 cm). In addition, this region of the sky has been targeted for extensive spectroscopy using the DEIMOS spectrograph on the Keck II 10 m telescope. Our survey is compared to other large multiwavelength surveys in terms of depth and sky coverage.

In the century since Einstein's anno mirabilis of 1905, our concept of the Universe has expanded from Kapteyn's flattened disk of stars only 10 kpc across to an observed horizon about 30 Gpc across that is only a tiny fraction of an immensely large inflated bubble. The expansion of our knowledge about the Universe, both in the types of data and the sheer quantity of data, has been just as dramatic. This talk will summarize this century of progress and our current understanding of the cosmos.

The Wide-field Infrared Survey Explorer (WISE), a NASA MIDEX mission, will survey the entire sky in four bands from 3.3 to 23 microns with a sensitivity 1000 times greater than the IRAS survey. The WISE survey will extend the Two Micron All Sky Survey into the thermal infrared and will provide an important catalog for the James Webb Space Telescope. Using 1024x1024 HgCdTe and Si:As arrays at 3.3, 4.7, 12 and 23 microns, WISE will find the most luminous galaxies in the universe, the closest stars to the Sun, and it will detect most of the main belt asteroids larger than 3 km. The single WISE instrument consists of a 40 cm diamond-turned aluminum afocal telescope, a two-stage solid hydrogen cryostat, a scan mirror mechanism, and reimaging optics giving 5" resolution (full-width-half-maximum). The use of dichroics and beamsplitters allows four color images of a 47'x47' field of view to be taken every 8.8 seconds, synchronized with the orbital motion to provide total sky coverage with overlap between revolutions. WISE will be placed into a Sun-synchronous polar orbit on a Delta 7320-10 launch vehicle. The WISE survey approach is simple and efficient. The three-axis-stabilized spacecraft rotates at a constant rate while the scan mirror freezes the telescope line of sight during each exposure. WISE is currently in its Preliminary Design Phase, with the mission Preliminary Design Review scheduled for July, 2005. WISE is scheduled to launch in mid 2009; the project web site can be found at www.wise.ssl.berkeley.edu.

We report a search for signatures of cosmic strings in the the Cosmic Microwave Background data from the Wilkinson Microwave Anisotropy Probe. We used a digital filter designed to search for individual cosmic strings and found no evidence for them in the WMAP CMB anisotropies to a level of $\Delta T/T \sim 0.29$ mK. This corresponds to an absence of cosmic strings with $ G\mu \ga 1.07 \times 10^{-5}$ for strings moving with velocity $v = c/\sqrt{2}$. Unlike previous work, this limit does not depend on an assumed string abundance. We have searched the WMAP data for evidence of a cosmic string recently reported as the CSL-1 object, and found an ``edge'' with 2$\sigma$ significance. However, if this edge is real and produced by a cosmic string, it would have to move at velocity $\ga$ 0.94c. We also present preliminary limits on the CMB data that will be returned by the PLANCK satellite for comparison. With the available information on the PLANCK satellite, we calculated that it would be twice as sensitive to cosmic strings as WMAP.

The infrared sky from space is the sum of a cosmic signal from galaxies, quasars, and perhaps more exotic sources; and foregrounds from the Milky Way and from the Solar System. At a distance of 1 AU from the Sun, the foreground from interplanetary dust is very bright between 5 and 100 microns, but ``very bright'' is still several million times fainter than the background produced by ground-based telescopes. In the near infrared 1-2.2 micron range the space infrared sky is a thousand times fainter than the OH nightglow from the Earth's atmosphere. As a result of these advantages, wide-field imaging from space in the infrared can be an incredibly sensitive method to study the Universe.

The IRAC shallow survey covers 8.5 square degrees in the NOAO Deep Wide-Field Survey in Bootes with 3 or more 30 second exposures per position. An overview of the survey design, reduction, calibration, star-galaxy separation, and initial results is provided. The survey includes approximately 370,000, 280,000, 38,000, and 34,000 sources brighter than the 5 sigma limits of 6.4, 8.8, 51, and 50 microJy at 3.6, 4.5, 5.8, and 8 microns respectively, including some with unusual spectral energy distributions.

The Spitzer Space Telescope, NASA's Great Observatory for infrared astronomy, was launched 2003 August 25 and is returning excellent scientific data from its Earth-trailing solar orbit. Spitzer combines the intrinsic sensitivity achievable with a cryogenic telescope in space with the great imaging and spectroscopic power of modern detector arrays to provide the user community with huge gains in capability for exploration of the cosmos in the infrared. The observatory systems are largely performing as expected and the projected cryogenic lifetime is in excess of 5 years. This paper summarizes the on-orbit scientific, technical and operational performance of Spitzer. Subsequent papers in this special issue describe the Spitzer instruments in detail and highlight many of the exciting scientific results obtained during the first six months of the Spitzer mission.

Infrared source counts at wavelengths 3 < lambda < 10 um cover more than 10 magnitudes in source brightness, four orders of magnitude in surface density, and reach an integrated surface density of 10^5 sources/deg^2. At m<14 mag, most of the sources are Galactic stars, in agreement with models. After removal of Galactic stars, galaxy counts are consistent with what few measurements exist at nearby wavelengths. At 3.6 and 4.5 um, the galaxy counts follow the expectations of a Euclidean world model down to ~16 mag and drop below the Euclidean curve for fainter magnitudes. Counts at these wavelengths begin to show decreasing completeness around magnitude 19.5. At 5.8 and 8 um, the counts relative to a Euclidean world model show a large excess at bright magnitudes. This is probably because local galaxies emit strongly in the aromatic dust (``PAH'') features. The counts at 3.6 um resolve <50% of the Cosmic Infrared Background at that wavelength.

The Cosmic Microwave Background (CMB) consists of photons that were last created about 2 months after the Big Bang, and last scattered about 380,000 years after the Big Bang. The spectrum of the CMB is very close to a blackbody at 2.725 K and upper limits on any deviations of the CMB from a blackbody place strong constraints on energy transfer between the CMB and matter at all redshifts less than 2,000,000. The CMB is very nearly isotropic, but a dipole anisotropy of +/-3.346(17) mK shows that the Solar System barycenter is moving at 368+/-2 km/sec relative to the observable Universe. The dipole corresponds to a spherical harmonic index l=1. The higher indices l geq 2 indicate intrinsic inhomogeneities in the Universe that existed at the time of last scattering. While the photons have traveled freely only since the time of last scattering, the inhomogeneities traced by the CMB photons have been in place since the inflationary epoch only 10^{-35} sec after the Big Bang. These intrinsic anisotropies are much smaller in amplitude than the dipole anisotropy, with Delta T leq 100 microK. Electron scattering of the anisotropic radiation field produces an anisotropic linear polarization in the CMB with amplitudes less than 5 microK. Detailed studies of the angular power spectrum of the temperature and linear polarization anisotropies have yielded precise values for many cosmological parameters. This paper will discuss the techniques necessary to measure signals that are 100 million times smaller than the emission from the instrument and briefly describe results from experiments up to WMAP.

We describe our methodology for comparing the WMAP measurements of the cosmic microwave background (CMB) and other complementary data sets to theoretical models. The unprecedented quality of the WMAP data, and the tight constraints on cosmological parameters that are derived, require a rigorous analysis so that the approximations made in the modeling do not lead to significant biases. We describe our use of the likelihood function to characterize the statistical properties of the microwave background sky. We outline the use of the Monte Carlo Markov Chains to explore the likelihood of the data given a model to determine the best fit cosmological parameters and their uncertainties. We add to the WMAP data the l>~700 CBI and ACBAR measurements of the CMB, the galaxy power spectrum at z~0 obtained from the 2dF galaxy redshift survey (2dFGRS), and the matter power spectrum at z~3 as measured with the Ly alpha forest. These last two data sets complement the CMB measurements by probing the matter power spectrum of the nearby universe. Combining CMB and 2dFGRS requires that we include in our analysis a model for galaxy bias, redshift distortions, and the non-linear growth of structure. We show how the statistical and systematic uncertainties in the model and the data are propagated through the full analysis.

The Wilkinson Microwave Anisotropy Probe (WMAP) has mapped the full sky in Stokes I, Q, and U parameters at frequencies 23, 33, 41, 61, and 94 GHz. We detect correlations between the temperature and polarization maps significant at more than 10 standard deviations. The correlations are present in all WMAP frequency bands with similar amplitude from 23 to 94 GHz, and are consistent with a superposition of a CMB signal with a weak foreground. The fitted CMB component is robust against different data combinations and fitting techniques. On small angular scales theta < 5 deg, the WMAP data show the temperature-polarization correlation expected from adiabatic perturbations in the temperature power spectrum. The data for l > 20 agree well with the signal predicted solely from the temperature power spectra, with no additional free parameters. We detect excess power on large angular scales (theta > 10 deg) compared to predictions based on the temperature power spectra alone. The excess power is well described by reionization at redshift 11 < z_r < 30 at 95% confidence, depending on the ionization history. A model-independent fit to reionization optical depth yields results consistent with the best-fit LambdaCDM model, with best fit value tau = 0.17 +- 0.04 at 68% confidence, including systematic and foreground uncertainties. This value is larger than expected given the detection of a Gunn-Peterson trough in the absorption spectra of distant quasars, and implies that the universe has a complex ionization history: WMAP has detected the signal from an early epoch of reionization.

Knowledge of the beam profiles is of critical importance for interpreting data from cosmic microwave background experiments. In this paper, we present the characterization of the in-flight optical response of the WMAP satellite. The main beam intensities have been mapped to < -30 dB of their peak values by observing Jupiter with the satellite in the same observing mode as for CMB observations. The beam patterns closely follow the pre-launch expectations. The full width at half maximum is a function of frequency and ranges from 0.82 degrees at 23 GHz to 0.21 degrees at 94 GHz; however, the beams are not Gaussian. We present: (a) the beam patterns for all ten differential radiometers and show that the patterns are substantially independent of polarization in all but the 23 GHz channel; (b) the effective symmetrized beam patterns that result from WMAP's compound spin observing pattern; (c) the effective window functions for all radiometers and the formalism for propagating the window function uncertainty; and (d) the conversion factor from point source flux to antenna temperature. A summary of the systematic uncertainties, which currently dominate our knowledge of the beams, is also presented. The constancy of Jupiter's temperature within a frequency band is an essential check of the optical system. The tests enable us to report a calibration of Jupiter to 1-3% accuracy relative to the CMB dipole.

The CMB has distinct peaks in both its temperature angular power spectrum (TT) and temperature-polarization cross-power spectrum (TE). From the WMAP data we find the first peak in the temperature spectrum at l = 220.1 +- 0.8 with an amplitude of 74.7 +- 0.5 microK; the first trough at l = 411.7 +- 3.5 with an amplitude of 41.0 +- 0.5 microK; and the second peak at l = 546 +- 10 with an amplitude of 48.8 +- 0.9 microK. The TE spectrum has an antipeak at l = 137 +- 9 with a cross-power of -35 +- 9 microK^2, and a peak at l = 329 +- 19 with cross-power 105 +- 18 microK^2. All uncertainties are 1 sigma and include calibration and beam errors. An intuition for how the data determine the cosmological parameters may be gained by limiting one's attention to a subset of parameters and their effects on the peak characteristics. We interpret the peaks in the context of a flat adiabatic LambdaCDM model with the goal of showing how the cosmic baryon density, Omega_b h^2, matter density, Omega_m h^2, scalar index, n_s, and age of the universe are encoded in their positions and amplitudes. To this end, we introduce a new scaling relation for the TE antipeak-to-peak amplitude ratio and recompute known related scaling relations for the TT spectrum in light of the WMAP data. From the scaling relations, we show that WMAP's tight bound on Omega_b h^2 is intimately linked to its robust detection of the first and second peaks of the TT spectrum.

We present limits to the amplitude of non-Gaussian primordial fluctuations in the WMAP 1-year cosmic microwave background sky maps. A non-linear coupling parameter, f_NL, characterizes the amplitude of a quadratic term in the primordial potential. We use two statistics: one is a cubic statistic which measures phase correlations of temperature fluctuations after combining all configurations of the angular bispectrum. The other uses the Minkowski functionals to measure the morphology of the sky maps. Both methods find the WMAP data consistent with Gaussian primordial fluctuations and establish limits, -58<f_NL<134, at 95% confidence. There is no significant frequency or scale dependence of f_NL. The WMAP limit is 30 times better than COBE, and validates that the power spectrum can fully characterize statistical properties of CMB anisotropy in the WMAP data to high degree of accuracy. Our results also validate the use of a Gaussian theory for predicting the abundance of clusters in the local universe. We detect a point-source contribution to the bispectrum at 41 GHz, b_src = (9.5+-4.4) X 1e-5 uK^3 sr^2, which gives a power spectrum from point sources of c_src = (15+-6) X 1e-3 uK^2 sr in thermodynamic temperature units. This value agrees well with independent estimates of source number counts and the power spectrum at 41 GHz, indicating that b_src directly measures residual source contributions.

WMAP precision data enables accurate testing of cosmological models. We find that the emerging standard model of cosmology, a flat Lambda-dominated universe seeded by nearly scale-invariant adiabatic Gaussian fluctuations, fits the WMAP data. With parameters fixed only by WMAP data, we can fit finer scale CMB measurements and measurements of large scle structure (galaxy surveys and the Lyman alpha forest). This simple model is also consistent with a host of other astronomical measurements. We then fit the model parameters to a combination of WMAP data with other finer scale CMB experiments (ACBAR and CBI), 2dFGRS measurements and Lyman alpha forest data to find the model's best fit cosmological parameters: h=0.71+0.04-0.03, Omega_b h^2=0.0224+-0.0009, Omega_m h^2=0.135+0.008-0.009, tau=0.17+-0.06, n_s(0.05/Mpc)=0.93+-0.03, and sigma_8=0.84+-0.04. WMAP's best determination of tau=0.17+-0.04 arises directly from the TE data and not from this model fit, but they are consistent. These parameters imply that the age of the universe is 13.7+-0.2 Gyr. The data favors but does not require a slowly varying spectral index. By combining WMAP data with other astronomical data sets, we constrain the geometry of the universe, Omega_tot = 1.02 +- 0.02, the equation of state of the dark energy w < -0.78 (95% confidence limit assuming w >= -1), and the energy density in stable neutrinos, Omega_nu h^2 < 0.0076 (95% confidence limit). For 3 degenerate neutrino species, this limit implies that their mass is less than 0.23 eV (95% confidence limit). The WMAP detection of early reionization rules out warm dark matter.

The Wilkinson Microwave Anisotropy Probe (WMAP) science team has released results from the first year of operation at the Earth-Sun L2 Lagrange point. The maps are consistent with previous observations but have much better sensitivity and angular resolution than the COBE DMR maps, and much better calibration accuracy and sky coverage than ground-based and balloon-borne experiments. The angular power spectra from these ground-based and balloon-borne experiments are consistent within their systematic and statistical uncertainties with the WMAP results. WMAP detected the large angular-scale correlation between the temperature and polarization anisotropies of the CMB caused by electron scattering since the Universe became reionized after the "Dark Ages", giving a value for the electron scattering optical depth of 0.17+/-0.04. The simplest Lambda-CDM model with n=1 and Omega_tot=1 provides an adequate fit to the WMAP data and gives parameters which are consistent with determinations of the Hubble constant and observations of the accelerating Universe using supernovae. The time-ordered data, maps, and power spectra from WMAP can be found at http://lambda.gsfc.nasa.gov along with 13 papers by the WMAP science team describing the results in detail.

Full sky maps are made in five microwave frequency bands to separate the temperature anisotropy of the CMB from foreground emission. We define masks that excise regions of high foreground emission. The effectiveness of template fits to remove foreground emission from the WMAP data is examined. These efforts result in a CMB map with minimal contamination and a demonstration that the WMAP CMB power spectrum is insensitive to residual foreground emission. We construct a model of the Galactic emission components. We find that the Milky Way resembles other normal spiral galaxies between 408 MHz and 23 GHz, with a synchrotron spectral index that is flattest (beta ~ -2.5) near star-forming regions, especially in the plane, and steepest (beta ~ -3) in the halo. The significant synchrotron index steepening out of the plane suggests a diffusion process in which the halo electrons are trapped in the Galactic potential long enough to suffer synchrotron and inverse Compton energy losses and hence a spectral steepening. The synchrotron index is steeper in the WMAP bands than in lower frequency radio surveys, with a spectral break near 20 GHz to beta < -3. The modeled thermal dust spectral index is also steep in the WMAP bands, with beta ~ 2.2. Microwave and H alpha measurements of the ionized gas agree. Spinning dust emission is limited to < ~5% of the Ka-band foreground emission. A catalog of 208 point sources is presented. Derived source counts suggest a contribution to the anisotropy power from unresolved sources of (15.0 +- 1.4) 10^{-3} microK^2 sr at Q-band and negligible levels at V-band and W-band.

We present full sky microwave maps in five bands (23 to 94 GHz) from the WMAP first year sky survey. Calibration errors are <0.5% and the low systematic error level is well specified. The 2<l<900 anisotropy power spectrum is cosmic variance limited for l<354 with a signal-to-noise ratio >1 per mode to l=658. The temperature-polarization cross-power spectrum reveals both acoustic features and a large angle correlation from reionization. The optical depth of reionization is 0.17 +/- 0.04, which implies a reionization epoch of 180+220-80 Myr (95% CL) after the Big Bang at a redshift of 20+10-9 (95% CL) for a range of ionization scenarios. This early reionization is incompatible with the presence of a significant warm dark matter density. The age of the best-fit universe is 13.7 +/- 0.2 Gyr old. Decoupling was 379+8-7 kyr after the Big Bang at a redshift of 1089 +/- 1. The thickness of the decoupling surface was dz=195 +/- 2. The matter density is Omega_m h^2 = 0.135 +0.008 -0.009, the baryon density is Omega_b h^2 = 0.0224 +/- 0.0009, and the total mass-energy of the universe is Omega_tot = 1.02 +/- 0.02. The spectral index of scalar fluctuations is fit as n_s = 0.93 +/- 0.03 at wavenumber k_0 = 0.05 Mpc^-1, with a running index slope of dn_s/d ln k = -0.031 +0.016 -0.018 in the best-fit model. This flat universe model is composed of 4.4% baryons, 22% dark matter and 73% dark energy. The dark energy equation of state is limited to w<-0.78 (95% CL). Inflation theory is supported with n_s~1, Omega_tot~1, Gaussian random phases of the CMB anisotropy, and superhorizon fluctuations. An admixture of isocurvature modes does not improve the fit. The tensor-to-scalar ratio is r(k_0=0.002 Mpc^-1)<0.90 (95% CL).

The Diffuse InfraRed Background Experiment on COBE measured the total infrared signal seen from space at a distance of 1 astronomical unit from the Sun. Using time variations as the Earth orbits the Sun, it is possible to remove most of the foreground signal produced by the interplanetary dust cloud [zodiacal light]. By correlating the DIRBE signal with the column density of atomic hydrogen measured using the 21 cm line, it is possible to remove most of the foreground signal produced by interstellar dust, although one must still be concerned by dust associated with H_2 (molecular gas) and H II (the warm ionized medium). DIRBE was not able to determine the CIRB in the 5-60 micron wavelength range, but did detect both a far infrared background and a near infrared background. The far infrared background has an integrated intensity of about 34 nW/m^2/sr, while the near infrared and optical extragalactic background has about 59 nW/m^2/sr. The Far InfraRed Absolute Spectrophotometer (FIRAS) on COBE has been used to constrain the long wavelength tail of the far infrared background but a wide range of intensities at 850 microns are compatible with the FIRAS data. Thus the fraction of the CIRB produced by SCUBA sources has large uncertainties in both the numerator and the denominator.

The theoretical basis for the prediction of anisotropies in the cosmic microwave background is very well developed. Very low amplitude density and temperature perturbations produce small gravitational effects, leading to an anisotropy that is a combination of temperature fluctuations at the surface of last scattering and gravitational redshifts both at last scattering and along the path to the observer. All of the primary anisotropy can be handled by linear perturbation theory, which allows a very accurate calculation of the predicted anisotropy from different models of the Universe.

We confront predictions of inflationary scenarios with the WMAP data, in combination with complementary small-scale CMB measurements and large-scale structure data. The WMAP detection of a large-angle anti-correlation in the temperature--polarization cross-power spectrum is the signature of adiabatic superhorizon fluctuations at the time of decoupling. The WMAP data are described by pure adiabatic fluctuations: we place an upper limit on a correlated CDM isocurvature component. Using WMAP constraints on the shape of the scalar power spectrum and the amplitude of gravity waves, we explore the parameter space of inflationary models that is consistent with the data. We place limits on inflationary models; for example, a minimally-coupled lambda phi^4 is disfavored at more than 3-sigma using WMAP data in combination with smaller scale CMB and large scale structure survey data. The limits on the primordial parameters using WMAP data alone are: n_s(k_0=0.002 Mpc^{-1})=1.20_{-0.11}^{+0.12}, dn/dlnk=-0.077^{+0.050}_{- 0.052}, A(k_0=0.002 Mpc}^{-1})=0.71^{+0.10}_{-0.11} (68% CL), and r(k_0=0.002 Mpc^{-1})<1.28 (95% CL).

The first-year WMAP data, in combination with any one of a number of other cosmic probes, show that we live in a flat \Lambda-dominated CDM universe with \Omega_m ~ 0.27 and \Omega_\Lambda ~ 0.73. In this model the late-time action of the dark energy, through the integrated Sachs-Wolfe effect, should produce CMB anisotropies correlated with matter density fluctuations at z<2 (Crittenden & Turok 1996). The measurement of such a signal is an important independent check of the model. We cross-correlate the NRAO VLA Sky Survey radio source catalog (Condon et al. 1998) with the WMAP data in search of this signal, and see indications of the expected correlation. Assuming a flat \Lambda-CDM cosmology, we find \Omega_\Lambda>0 (95% CL, statistical errors only) with the peak of the likelihood at \Omega_\Lambda=0.68, consistent with the preferred WMAP value. A closed model with \Omega_m=1.28, h=0.33, and no dark energy component (\Omega_\Lambda=0), marginally consistent with the WMAP CMB TT angular power spectrum, would produce an anti-correlation between the matter distribution and the CMB. Our analysis of the cross-correlation of the WMAP data with the NVSS catalog rejects this cosmology at the 3\sigma level.

The Microwave Anisotropy Probe (MAP) satellite, launched June 30, 2001, will produce full sky maps of the cosmic microwave background radiation in 5 frequency bands spanning 20 - 106 GHz. MAP contains 20 differential radiometers built with High Electron Mobility Transistor (HEMT) amplifiers with passively cooled input stages. The design and test techniques used to evaluate and minimize systematic errors and the pre-launch performance of the radiometers for all five bands are presented.

We present the angular power spectrum derived from the first-year Wilkinson Microwave Anisotropy Probe (WMAP) sky maps. We study a variety of power spectrum estimation methods and data combinations and demonstrate that the results are robust. The data are modestly contaminated by diffuse Galactic foreground emission, but we show that a simple Galactic template model is sufficient to remove the signal. Point sources produce a modest contamination in the low frequency data. After masking ~700 known bright sources from the maps, we estimate residual sources contribute ~3500 uK^2 at 41 GHz, and ~130 uK^2 at 94 GHz, to the power spectrum l*(l+1)*C_l/(2*pi) at l=1000. Systematic errors are negligible compared to the (modest) level of foreground emission. Our best estimate of the power spectrum is derived from 28 cross-power spectra of statistically independent channels. The final spectrum is essentially independent of the noise properties of an individual radiometer. The resulting spectrum provides a definitive measurement of the CMB power spectrum, with uncertainties limited by cosmic variance, up to l~350. The spectrum clearly exhibits a first acoustic peak at l=220 and a second acoustic peak at l~540 and it provides strong support for adiabatic initial conditions. Kogut et al. (2003) analyze the C_l^TE power spectrum, and present evidence for a relatively high optical depth, and an early period of cosmic reionization. Among other things, this implies that the temperature power spectrum has been suppressed by \~30% on degree angular scales, due to secondary scattering.

The WMAP satellite has completed one year of measurements of the Cosmic Microwave Background (CMB) radiation using 20 differential high-electron-mobility-transistor (HEMT) based radiometers. All the radiometers are functioning nominally, and characterizations of the on-orbit radiometer performance are presented, with an emphasis on properties that are required for the production of sky maps from the time ordered data. A radiometer gain model, used to smooth and interpolate the CMB dipole gain measurements is also presented. No degradation in the sensitivity of any of the radiometers has been observed during the first year of observations.

We describe the calibration and data processing methods used to generate full-sky maps of the cosmic microwave background (CMB) from the first year of Wilkinson Microwave Anisotropy Probe (WMAP) observations. Detailed limits on residual systematic errors are assigned based largely on analyses of the flight data supplemented, where necessary, with results from ground tests. The data are calibrated in flight using the dipole modulation of the CMB due to the observatory's motion around the Sun. This constitutes a full-beam calibration source. An iterative algorithm simultaneously fits the time-ordered data to obtain calibration parameters and pixelized sky map temperatures. The noise properties are determined by analyzing the time-ordered data with this sky signal estimate subtracted. Based on this, we apply a pre-whitening filter to the time-ordered data to remove a low level of 1/f noise. We infer and correct for a small ~1% transmission imbalance between the two sky inputs to each differential radiometer, and we subtract a small sidelobe correction from the 23 GHz (K band) map prior to further analysis. No other systematic error corrections are applied to the data. Calibration and baseline artifacts, including the response to environmental perturbations, are negligible. Systematic uncertainties are comparable to statistical uncertainties in the characterization of the beam response. Both are accounted for in the covariance matrix of the window function and are propagated to uncertainties in the final power spectrum. We characterize the combined upper limits to residual systematic uncertainties through the pixel covariance matrix.

Since the Galactic center is ~1000 times brighter than fluctuations in the Cosmic Microwave Background (CMB), CMB experiments must carefully account for stray Galactic pickup. We present the level of contamination due to sidelobes for the year one CMB maps produced by the WMAP observatory. For each radiometer, full 4 pi sr antenna gain patterns are determined from a combination of numerical prediction, ground-based and space-based measurements. These patterns are convolved with the WMAP year one sky maps and observatory scan pattern to generate expected sidelobe signal contamination, for both intensity and polarized microwave sky maps. Outside of the Galactic plane, we find rms values for the expected sidelobe pickup of 15, 2.1, 2.0, 0.3, 0.5 uK for K, Ka, Q, V, and W bands respectively. Except at K band, the rms polarized contamination is <<1 uK. Angular power spectra of the Galactic pickup are presented.

The primary goal of the MAP satellite, now in orbit, is to make high fidelity polarization sensitive maps of the full sky in five frequency bands between 20 and 100 GHz. From these maps we will characterize the properties of the cosmic microwave background (CMB) anisotropy and Galactic and extragalactic emission on angular scales ranging from the effective beam size, <0.23 degree, to the full sky. MAP is a differential microwave radiometer. Two back-to-back shaped offset Gregorian telescopes feed two mirror symmetric arrays of ten corrugated feeds. We describe the prelaunch design and characterization of the optical system, compare the optical models to the measurements, and consider multiple possible sources of systematic error.

The purpose of the MAP mission is to determine the geometry, content, and evolution of the universe via a 13 arcmin full-width-half-max (FWHM) resolution full sky map of the temperature anisotropy of the cosmic microwave background radiation with uncorrelated pixel noise, minimal systematic errors, multifrequency observations, and accurate calibration. These attributes were key factors in the success of NASA's Cosmic Background Explorer (COBE) mission, which made a 7 degree FWHM resolution full sky map, discovered temperature anisotropy, and characterized the fluctuations with two parameters, a power spectral index and a primordial amplitude. Following COBE considerable progress has been made in higher resolution measurements of the temperature anisotropy. With 45 times the sensitivity and 33 times the angular resolution of the COBE mission, MAP will vastly extend our knowledge of cosmology. MAP will measure the physics of the photon-baryon fluid at recombination. From this, MAP measurements will constrain models of structure formation, the geometry of the universe, and inflation. In this paper we present a pre-launch overview of the design and characteristics of the MAP mission. This information will be necessary for a full understanding of the MAP data and results, and will also be of interest to scientists involved in the design of future cosmic microwave background experiments and/or space science missions.

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

The observation of SN 1997ff at redshift 1.7 has been claimed to refute alternative models such as grey dust or evolution for the faintness of distant supernovae, leaving only an accelerating Universe as a viable model. However, a very simple one parameter evolution model, with the peak luminosity varying as an exponential function of cosmic time, converts the flux vs. distance law of the critical density matter-dominated model into that of the concordance Omega_matter = 0.3 flat vacuum-dominated model with an error no larger than 0.03 mag over the range 0-2 in redshift. A grey dust model that matches this accuracy can easily be contrived but it still fails by overproducing the far-IR background or distorting the CMB. Models that involve oscillation between photons and axions could emulate an exponential function of cosmic time without violating these background constraints. Clearly a better and well-tested understanding of the Type Ia supernova explosion mechanism and the origin of the correlation between the decay rate and luminosity is needed before any effort to reduce statistical errors in the supernova Hubble diagram to very small levels.

Stellar fluxes from the 2MASS catalog are used to remove the contribution due to Galactic stars from the intensity measured by DIRBE in 13 regions in the North and South Galactic polar caps. Allowing for a constant calibration factor between the DIRBE 3.5 micron intensity and the 2MASS 2.2 micron fluxes gives very small pixel to pixel scatter. After subtracting the interplanetary and Galactic foregrounds, a residual intensity of 16.1 +/- 4 kJy/sr or 13.8 +/- 3.4 nW/m^2/sr at 3.5 microns is found. A similar analysis at 2.2 microns gives a residual intensity of 17.6 +/- 4.4 kJy/sr or 24 +/- 6 nW/m^2/sr. The intercepts of the DIRBE minus 2MASS correlation at 1.25 microns show more scatter and are a much smaller fraction of the foreground, leading to a weak limit on the CIRB of 10.1 +/- 7.4 kJy/sr or 24.3 +/- 18 nW/m^2/sr (1 sigma).

Models of the zodiacal light are necessary to convert measured data taken from low Earth orbit into the radiation field outside the solar system. The uncertainty in these models dominates the overall uncertainty in determining the extragalactic background light for wavelengths < 100 microns.

The Sloan Digital Sky Survey [SDSS] has measured an optical luminosity function for galaxies in 5 bands, finding 1.5 to 2.1 times more luminosity density than previous work. This note compares the SDSS luminosity density to two recent determinations of the near infrared luminosity function based on 2MASS data, and finds that an extrapolation of the SDSS results gives a 2.3 times greater near infrared luminosity density.

Stellar fluxes from the 2MASS catalog are used to remove the contribution due to Galactic stars from the intensity measured by DIRBE in four regions in the North and South Galactic polar caps. After subtracting the interplanetary and galactic foregrounds, a consistent residual intensity of 14.8 +/- 4.6 kJy/sr or 20.2 +/- 6.3 nW/m^2/sr at 2.2 microns is found. At 1.25 microns the residuals show more scatter and are a much smaller fraction of the foreground, leading to a weak limit on the CIRB of 12.0 +/- 6.8 kJy/sr or 28.9 +/- 16.3 nW/m^2/sr (1 sigma).

The Cosmic InfraRed Background (CIRB) is the sum total of the redshifted and reprocessed short wavelength radiation from the era of galaxy formation, and hence contains vital information about the history of galactic evolution. One of the main problems associated with estimating an isotropic CIRB in the near infrared (1-5 microns) is the unknown contribution from stars within our own galaxy. The optimal observational window to search for a background in the near-IR is at 3.5 microns since that is the wavelength region where the other main foreground, the zodiacal dust emission, is the least. It is not possible to map out the entire 3.5 micron sky at a resolution which will accurately estimate the flux from stars. However, since the CIRB is presumably isotropic, it can potentially be detected by selecting a smaller field and imaging it at good resolution to estimate the stellar intensity. We selected a 2x2 degree "dark spot" near the North Galactic Pole which had the least intensity at 3.5 microns after a zodiacal light model was subtracted from the all-sky maps generated by the Diffuse InfraRed Background Experiment (DIRBE). The measured total intensity of the few bright stars in this field was combined with a model for the contribution from dimmer stars and subtracted from the zodi-subtracted DIRBE map. The contribution from the interstellar medium was also subtracted leaving a residual intensity at 2.2 microns of: 16.4+/-4.4 kJy/sr or 22.4+/-6 nW/m^2/sr, and at 3.5 microns: 12.8+/-3.8 kJy/sr or 11+/-3.3 nW/m^2/sr. [Abridged]

We compare data from the Diffuse InfraRed Background Experiment (DIRBE) on the Cosmic Background Explorer (COBE) satellite to the the Wainscoat et al. (1992) model of the infrared sky. The model is first compared with broadband K (2.2 microns) star counts. Its success at K gives credence to its physical approach which is extrapolated to the L band (3.5 microns). We have analyzed the histograms of the pixel by pixel intensities in the 2.2 and 3.5 micron maps from DIRBE after subtracting the zodiacal light. The shape of these histograms agrees quite well with the histogram shape predicted using the Wainscoat et al. model of the infrared sky, but the predicted histograms must be displaced by a constant intensity in order to match the data. This shift is the cosmic infrared background, which is 16.9+/-4.4 kJy/sr or 23.1+/-5.9 nW/m^2/sr at 2.2 microns, and 14.4+/-3.7 kJy/sr or 12.4+/-3.2 nW/m^2/sr at 3.5 microns.

A major goal of modern astrophysics is to understand the processes by which the universe evolved from its initial simplicity, as seen in measurements of the Cosmic Microwave Background, to the universe we see today, with complexity on all scales. While the diffuse background measurements of COBE reveal the importance of the far infrared and submillimeter in early galaxy and star formation, the understanding of the development of complex structure requires high resolution imaging and spectroscopy. We present a concept for a space mission called SPECS, the Submillimeter Probe of the Evolution of Cosmic Structure, which provides these capabilities. SPECS is a cold, spatial and spectral Michelson interferometer with adjustable baselines ranging up to 1 km. It has Hubble sensitivity and angular resolution in the far IR and submillimeter, spectral resolution up to 10,000, and a 14' field of view. SPECS will be able to image thermal dust continuum and infrared cooling and diagnostic line emission over a wide range of redshifts, providing extinction-free astrophysical probes of young galaxies and early cosmic structures, and measures of the luminosity and heavy element formation history of the universe. SPECS would also have the potential to improve vastly our knowledge of protostars, protoplanetary systems, Active Galactic Nuclei and other objects in the local universe. We recommend that a concerted effort be made during the next decade to develop and test certain critical technologies (photon-counting far IR detectors, formation flying spacecraft, cold, lightweight mirrors, and active coolers), so that SPECS can be deployed in about 2015.

We present constraints on the density of halo dark matter candidates within the solar circle based on the anisotropy in the high energy gamma-ray background. The known galactic components of the gamma-ray background, in particular the inverse Compton component, have been estimated more accurately. We find the spectrum of the residual emission, after subtracting the galactic component is inconsistent with emission from some of the proposed dark matter candidates. We derive upper limits of 10^8 M_sun for the mass of diffuse gas and 3*10^9 pc^(-3) for the number density of primordial black holes contributing to the gamma-ray background.

The COBE Diffuse Infrared Background Experiment (DIRBE) was designed to search for the cosmic infrared background (CIB) radiation. Scattered light and thermal emission from the interplanetary dust (IPD) are major contributors to the diffuse sky brightness at most infrared wavelengths. Accurate removal of this zodiacal light foreground is a necessary step toward a direct measurement of the CIB. The zodiacal light foreground contribution in each of the 10 DIRBE wavelength bands ranging from 1.25 to 240 microns is distinguished by its apparent seasonal variation over the whole sky. This contribution has been extracted by fitting the brightness calculated from a parameterized physical model to the time variation of the all-sky DIRBE measurements over 10 months of observations. The model brightness is evaluated as the integral along the line of sight of the product of a source function and a three-dimensional dust density distribution function. The dust density distribution is composed of multiple components: a smooth cloud, three asteroidal dust bands, and a circumsolar ring near 1 A.U. By using a directly measurable quantity which relates only to the IPD cloud, we exclude other contributors to the sky brightness from the IPD model. Using the IPD model described here, high-quality maps of the infrared sky with the zodiacal foreground removed have been generated. Imperfections in the model reveal themselves as low-level systematic artifacts in the residual maps which correlate with components of the IPD. The most evident of these artifacts are located near the ecliptic plane in the mid-infrared, and are less than 2% of the zodiacal foreground brightness. Uncertainties associated with the model are discussed, including implications for the CIB search.

The DIRBE on the COBE spacecraft was designed primarily to conduct systematic search for an isotropic CIB in ten photometric bands from 1.25 to 240 microns. The results of that search are presented here. Conservative limits on the CIB are obtained from the minimum observed brightness in all-sky maps at each wavelength, with the faintest limits in the DIRBE spectral range being at 3.5 microns (\nu I_\nu < 64 nW/m^2/sr, 95% CL) and at 240 microns (\nu I_\nu < 28 nW/m^2/sr, 95% CL). The bright foregrounds from interplanetary dust scattering and emission, stars, and interstellar dust emission are the principal impediments to the DIRBE measurements of the CIB. These foregrounds have been modeled and removed from the sky maps. Assessment of the random and systematic uncertainties in the residuals and tests for isotropy show that only the 140 and 240 microns data provide candidate detections of the CIB. The residuals and their uncertainties provide CIB upper limits more restrictive than the dark sky limits at wavelengths from 1.25 to 100 microns. No plausible solar system or Galactic source of the observed 140 and 240 microns residuals can be identified, leading to the conclusion that the CIB has been detected at levels of \nu I_\nu = 25+-7 and 14+-3 nW/m^2/sr at 140 and 240 microns respectively. The integrated energy from 140 to 240 microns, 10.3 nW/m^2/sr, is about twice the integrated optical light from the galaxies in the Hubble Deep Field, suggesting that star formation might have been heavily enshrouded by dust at high redshift. The detections and upper limits reported here provide new constraints on models of the history of energy-releasing processes and dust production since the decoupling of the cosmic microwave background from matter.

The Cosmic Infrared Background (CIB) is hidden behind veils of foreground emission from our own solar system and Galaxy. This paper describes procedures for removing the Galactic IR emission from the 1.25 - 240 micron COBE DIRBE maps as steps toward the ultimate goal of detecting the CIB. The Galactic emission models are carefully chosen and constructed so that the isotropic CIB is completely retained in the residual sky maps. We start with DIRBE data from which the scattered light and thermal emission of the interplanetary dust (IPD) cloud have already been removed. Locations affected by the emission from bright compact and stellar sources are excluded from the analysis. The unresolved emission of faint stars at near- and mid-IR wavelengths is represented by a model based on Galactic source counts. The 100 micron DIRBE observations are used as the spatial template for the interstellar medium (ISM) emission at high latitudes. Correlation of the 100 micron data with H I column density allows us to isolate the component of the observed emission that is associated with the ISM. Limits are established on the far-IR emissivity of the diffuse ionized medium, which indicate a lower emissivity per H nucleus than in the neutral medium. At 240 micron, we find that adding a second spatial template to the ISM model can greatly improve the accuracy of the model at low latitudes. The crucial product of this analysis is a set of all-sky IR maps from which the Galactic (and IPD) emission has been removed. We discuss systematic uncertainties and potential errors in the foreground subtraction process that may have an impact on studies seeking to detect the CIB in the residual maps.

In this paper we examine the cosmological constraints of the recent DIRBE and FIRAS detection of the extragalactic background light between 125-5000 microns on the metal and star formation histories of the universe.

The effect of a cosmological constant on the precession of the line of apsides is O(\Lambda c^2 r^3/GM) which is 3(H_\circ P)^2/8\pi^2 \approx 10^{-23} for a vacuum-dominated Universe with Hubble constant H_\circ = 65 km/sec/Mpc and for the orbital period P = 88 days of Mercury. This is unmeasurably small, so planetary perturbations cannot be used to limit the cosmological constant, contrary to the suggestion by Cardona & Tejeiro (1998).

The angular power spectra of the infrared maps obtained by the DIRBE (Diffuse InfraRed Background Experiment) instrument on the COBE satellite have been obtained by two methods: the Hauser-Peebles method previously applied to the DMR maps, and by Fourier transforming portions of the all-sky maps projected onto a plane. The two methods give consistent results, and the power spectrum of the high-latitude dust emission is C_\ell \propto \ell^{-3} in the range 2 < \ell < 300.

The Cosmic Microwave Background Radiation (CMBR) is the radiation left over from the hot Big Bang. Its blackbody spectrum and small anisotropy provide clues about the origin and early evolution of the Universe. In particular, the spectrum of the CMBR rules out many non-gravitational models of structure formation, and the anisotropy of the CMBR provides a measure of the gravitational potential at the time of last scattering, about 300,000 years after the Big Bang. The density inhomogeneities needed to produce the gravitational potential perturbations traced by the CMBR have grown to become the galaxies, clusters of galaxies, and superclusters that we see today.

The FIRAS instrument on the COBE satellite has conducted an unbiased survey of the far-infrared emission from our Galaxy. The first results of this survey were reported by Wright et al. (1991), and later conclusions were reported in Bennett et al. (1994). I report the results of analyses of this spectral survey, which includes emission lines from 158 micron C+, 122 micron and 205 micron N+, 370 micron and 609 micron neutral C, and CO J=2-1 through 5-4. The morphological distribution along the galactic plane (b=0 deg) for all the spectral line emission, and the high galactic latitude intensities of the strong C+ and 205 micron N+ emission are discussed. The high galactic latitude intensity of the 158 micron fine structure transition from C+ is I = (1.43 +/- 0.12)*1E-6 csc|b| erg/cm^2/s/sr for |b| > 15 deg and it decreases more rapidly than the far infrared intensity with increasing galactic latitude.

CMB anisotropy experiments seeking to make maps with more pixels than the 6144 pixels used by the COBE DMR need to address the practical issues of the computer time and storage required to make maps. A simple, repetitive scan pattern reduces these requirements but leaves the experiment vulnerable to systematic errors and striping in the maps. In this paper I give a time-ordered method for map-making with one-horned experiments that has reasonable memory and CPU needs but can handle complex COBE-like scans paths and 1/f noise.

We have refined the analysis of the data from the FIRAS (Far InfraRed Absolute Spectrophotometer) on board the COBE (COsmic Background Explorer). The FIRAS measures the difference between the cosmic microwave background and a precise blackbody spectrum. We find new tighter upper limits on general deviations from a blackbody spectrum. The RMS deviations are less than 50 parts per million of the peak of the CMBR. For the Comptonization and chemical potential we find $|y| < 15\times10^{-6}$ and $|\mu| < 9\times10^{-5}$ (95\% CL). There are also refinements in the absolute temperature, 2.728 $\pm$ 0.004 K (95\% CL), and dipole direction, $(\ell,b)=(264.14^\circ\pm0.30, 48.26^\circ\pm0.30)$ (95\% CL), and amplitude, $3.372 \pm 0.007$ mK (95\% CL). All of these results agree with our previous publications.

We employ a pixel-based likelihood technique to estimate the angular power spectrum of the COBE Differential Microwave Radiometer (DMR) 4-year sky maps. The spectrum is consistent with a scale-invariant power-law form with a normalization, expressed in terms of the expected quadrupole anisotropy, of Q_{rms-PS|n=1} = 18 +/- 1.4 uK, and a best-fit spectral index of 1.2 +/- 0.3. The normalization is somewhat smaller than we concluded from the 2-year data, mainly due to additional Galactic modeling. We extend the analysis to investigate the extent to which the ``small" quadrupole observed in our sky is statistically consistent with a power-law spectrum. The most likely quadrupole amplitude is somewhat dependent on the details of Galactic foreground subtraction and data selection, ranging between 7 and 10 uK, but in no case is there compelling evidence that the quadrupole is too small to be consistent with a power-law spectrum. We conclude with a likelihood analysis of the band power amplitude in each of four spectral bands between l = 2 and 40, and find no evidence for deviations from a simple power-law spectrum.

The angular power spectrum estimator developed by Peebles (1973) and Hauser & Peebles (1973) has been modified and applied to the 4 year maps produced by the COBE DMR. The power spectrum of the observed sky has been compared to the power spectra of a large number of simulated random skies produced with noise equal to the observed noise and primordial density fluctuation power spectra of power law form, with $P(k) \propto k^n$. The best fitting value of the spectral index in the range of spatial scales corresponding to spherical harmonic indices $3 \leq \ell \lesssim 30$ is an apparent spectral index $n_{app}$ = 1.13 (+0.3) (-0.4) which is consistent with the Harrison-Zel'dovich primordial spectral index $n_{pri} = 1$ The best fitting amplitude for $n_{app} = 1$ is $\langle Q_{RMS}^2\rangle^{0.5}$ = 18 uK.

We use the COBE Differential Microwave Radiometers (DMR) 4-year sky maps to model Galactic microwave emission at high latitudes (|b| > 20 deg). Cross-correlation of the DMR maps with Galactic template maps detects fluctuations in the high-latitude microwave sky brightness with the angular variation of the DIRBE far-infrared dust maps and a frequency dependence consistent with a superposition of dust and free-free emission. We find no significant correlations between the DMR maps and various synchrotron templates. On the largest angular scales (e.g., quadrupole), Galactic emission is comparable in amplitude to the anisotropy in the cosmic microwave background (CMB). The CMB quadrupole amplitude, after correction for Galactic emission, has amplitude $Q_{rms}$ = 10.7 uK with random uncertainty 3.6 uK and systematic uncertainty 7.1 uK from uncertainty in our knowledge of Galactic microwave emission.

The 2-point temperature correlation function is evaluated from the 4-year COBE DMR microwave anisotropy maps. We examine the 2-point function, which is the Legendre transform of the angular power spectrum, and show that the data are statistically consistent from channel to channel and frequency to frequency. The most likely quadrupole normalization is computed for a scale-invariant power-law spectrum of CMB anisotropy, using a variety of data combinations. For a given data set, the normalization inferred from the 2-point data is consistent with that inferred by other methods. The smallest and largest normalization deduced from any data combination are 16.4 and 19.6 uK respectively, with a value ~18 uK generally preferred.

The Differential Microwave Radiometers (DMR) instrument aboard the Cosmic Background Explorer (COBE) has mapped the full microwave sky to mean sensitivity 26 uK per 7 deg field of view. The absolute calibration is determined to 0.7% with drifts smaller than 0.2% per year. We have analyzed both the raw differential data and the pixelized sky maps for evidence of contaminating sources such as solar system foregrounds, instrumental susceptibilities, and artifacts from data recovery and processing. Most systematic effects couple only weakly to the sky maps. The largest uncertainties in the maps result from the instrument susceptibility to the Earth's magnetic field, microwave emission from the Earth, and upper limits to potential effects at the spacecraft spin period. Systematic effects in the maps are small compared to either the noise or the celestial signal: the 95% confidence upper limit for the pixel-pixel rms from all identified systematics is less than 6 uK in the worst channel. A power spectrum analysis of the (A-B)/2 difference maps shows no evidence for additional undetected systematic effects.

The cosmic microwave background radiation provides unique constraints on cosmological models. In this Letter we present a summary of the spatial properties of the cosmic microwave background radiation based on the full 4 years of COBE DMR observations, as detailed in a set of companion Letters. The anisotropy is consistent with a scale-invariant power law model and Gaussian statistics. With full use of the multi-frequency 4-year DMR data, including our estimate of the effects of Galactic emission, we find a power-law spectral index of $n=1.2\pm 0.3$ and a quadrupole normalization $Q_{rms-PS}=15.3^{+3.8}_{-2.8}$ $\mu$K. For $n=1$ the best-fit normalization is $Q_{rms-PS}\vert_{n=1}=18\pm 1.6$ $\mu$K. These values are consistent with both our previous 1-year and 2-year results. The results include use of the $\ell=2$ quadrupole term; exclusion of this term gives consistent results, but with larger uncertainties. The 4-year sky maps, presented in this Letter, portray an accurate overall visual impression of the anisotropy since the signal-to-noise ratio is ~2 per 10 degree sky map patch. The improved signal-to-noise ratio of the 4-year maps also allows for improvements in Galactic modeling and limits on non-Gaussian statistics.

We search the high-latitude portion of the COBE Differential Microwave Radiometers (DMR) 4-year sky maps for evidence of a non-Gaussian temperature distribution in the cosmic microwave background. The genus, 3-point correlation function, and 2-point correlation function of temperature maxima and minima are all in excellent agreement with the hypothesis that the CMB anisotropy on angular scales of 7 degrees or larger represents a random-phase Gaussian field. A likelihood comparison of the DMR sky maps to a set of random-phase non-Gaussian toy models selects the exact Gaussian model as most likely. Monte Carlo simulations show that the 2-point correlation of the peaks and valleys in the maps provides the greatest discrimination among the class of models tested.

Fourier analysis and power spectrum estimation of the cosmic microwave background anisotropy on an incompletely sampled sky developed by Gorski (1994) has been applied to the high-latitude portion of the 4-year COBE DMR 31.5, 53 and 90 GHz sky maps. Likelihood analysis using newly constructed Galaxy cuts (extended beyond |b| = 20deg to excise the known foreground emission) and simultaneously correcting for the faint high latitude galactic foreground emission is conducted on the DMR sky maps pixelized in both ecliptic and galactic coordinates. The Bayesian power spectrum estimation from the foreground corrected 4-year COBE DMR data renders n ~ 1.2 +/- 0.3, and Q_{rms-PS} ~ 15.3^{+3.7}_{-2.8} microK (projections of the two-parameter likelihood). These results are consistent with the Harrison-Zel'dovich n=1 model of amplitude Q_{rms-PS} ~ 18 microK detected with significance exceeding 14sigma (dQ/Q < 0.07). (A small power spectrum amplitude drop below the published 2-year results is predominantly due to the application of the new, extended Galaxy cuts.)

A major goal of cosmology is to obtain sensitive, high resolution maps of the Cosmic Microwave Background (CMB) anisotropy. Such maps, as would be produced by the recently proposed Microwave Anisotropy Probe (MAP), will contain a wealth of primary information about conditions in the early universe. To mitigate systematic effects when observing the microwave background, it is desirable for the raw data to be collected in differential form: as a set of temperature differences between points in the sky. However, the production of large (mega-pixel) maps from a set of temperature differences is a potentially severe computational challenge. We present a new technique for producing maps from differential radiometer data that has a computational cost that grows in the slowest possible way with increasing angular resolution and number of map pixels. The required central processor (CPU) time is proportional to the number of differential data points and the required random access memory (RAM) is proportional to the number of map pixels. We test our technique, and demonstrate its feasibility, by simulating one year of a space-borne anisotropy mission.

Calculations of the rate of local Primordial Black Hole explosions often assume that the PBH's can be highly concentrated into galaxies, thereby weakening the Page-Hawking limit on the cosmological density of PBH's. But if the PBH's are concentrated by a factor exceeding $c/(H_\circ R_\circ) \approx 4 \times 10^5$, where $R_\circ = 8.5$ kpc is the scale of the Milky Way, then the steady emission from the PBH's in the halo will produce an anisotropic high latitude diffuse gamma ray intensity larger than the observed anisotropy. This provides a limit on the rate-density of evaporating PBH's of $\lesssim 0.4$~pc$^{-3}$yr$^{-1}$ which is more than 6 orders of magnitude lower than recent experimental limits. However, the weak observed anisotropic high latitude diffuse gamma ray intensity is consistent with the idea that the dark matter that closes the Universe is Planck mass remnants of evaporated black holes.

We derive Galactic continuum spectra from 5-96/cm from COBE/FIRAS observations. The spectra are dominated by warm dust emission, which may be fit with a single temperature in the range 16-21 K (for nu^2 emissivity) along each line of sight. Dust heated by the attenuated radiation field in molecular clouds gives rise to intermediate temperature (10-14 K) emission in the inner Galaxy only. A widespread, very cold component (4-7 K) with optical depth that is spatially correlated with the warm component is also detected. The cold component is unlikely to be due to very cold dust shielded from starlight, because it is present at high latitude. We consider hypotheses that the cold component is due to enhanced submillimeter emissivity of the dust that gives rise to the warm component, or that it may be due to very small, large, or fractal particles. Lack of substantial power above the emission from warm dust places strong constraints on the amount of cold gas in the Galaxy. The microwave sky brightness due to interstellar dust is dominated by the cold component, and its angular variation could limit our ability to discern primordial fluctuations in the cosmic microwave background radiation.

The Quasi-Steady-State Cosmology as proposed by Hoyle, Burbidge and Narlikar does not fit the observed facts of the Universe. In particular, it predicts that 75-90\% of the radio sources in the brightest sample that shows steeper than Euclidean source counts should be blueshifted.

The frequency-independent RMS temperature fluctuations determined from the COBE-DMR two year sky maps are used to infer the parameter Q_{rms-PS}, which characterizes the normalization of power law models of primordial cosmological temperature anisotropy. In particular, a 'cross'-RMS statistic is used to determine Q_{rms-PS} for a forced fit to a scale-invariant Harrison-Zel'dovich (n = 1) spectral model. Using a joint analysis of the 7 degree and 10 degree RMS temperature derived from both the 53 and 90 GHz sky maps, we find Q_{rms-PS} = 17.0^{+2.5}_{-2.1} uK when the low quadrupole is included, and Q_{rms-PS} = 19.4^{+2.3}_{-2.1} uK excluding the quadrupole. These results are consistent with the n = 1 fits from more sensitive methods (e.g. power spectrum, correlation function). The effect of the low quadrupole derived from the COBE-DMR data on the inferred Q_{rms-PS} normalization is investigated. A bias to lower Q_{rms-PS} is found when the quadrupole is included. The higher normalization for a forced n = 1 fit is then favored by the cross-RMS technique. As initially pointed out in Wright et al. (1994a) and further discussed here, analytic formulae for the RMS sky temperature fluctuations will NOT provide the correct normalization amplitude.

The index $n$ of a power law power spectrum of primordial density fluctuations, $P(k) \propto k^n$, has been estimated using many different techniques. The most precise compare the COBE DMR large angular scale $\Delta T$ to the amplitude of the large scale structure, but these are also the most model-dependent. The COBE DMR $\Delta T$ has also been compared to the degree-scale $\Delta T$ from several experiments. And finally, a relatively model-independent value of $n$ can derived from the COBE data alone, but the small range of angular scales covered by COBE limits the precision of these methods.

The observations of all three COBE instruments are examined for the effects of dark matter. The anisotropy measured by the DMR, and especially the degree-scale ground- and balloon-based experiments, is only compatible with large-scale structure formation by gravity if the Universe is dominated by non-baryonic dark matter. The FIRAS instrument measures the total power radiated by cold dust, and thus places tight limits on the absorption of starlight by very cold gas and dust in the outer Milky Way. The DIRBE instrument measures the infrared background, and will place tight limits on the emission by low mass stars in the Galactic halo.

Some of the technical details involved in taking and analyzing data from COBE are discussed, and recent results from the FIRAS and DMR experiments are summarized. Some of the cosmological implications of these recent data are presented.

Starlight from distant galaxies ($z > 3$) is redshifted into the near infrared band with observed wavelengths from 2-8 $\mu$m. Most of the light is emitted by stars that have a peak emission at the 1.6 $\mu$m wavelength of the minimum of the H$^-$ opacity. We present simulated images of galaxy and star fields at 3, 4.7 and 8 $\mu$m, using the expected performance of the infrared camera on the Space InfraRed Telescope Facility (SIRTF). Standard astronomical image processing tools are used to locate sources, distinguish stars from galaxies, and create color-magnitude and color-color diagrams for the galaxies.

The first two years of COBE DMR observations of the CMB anisotropy are analyzed and compared with our previously published first year results. The results are consistent, but the addition of the second year of data increases the precision and accuracy of the detected CMB temperature fluctuations. The two-year 53 GHz data are characterized by RMS temperature fluctuations of DT=44+/-7 uK at 7 degrees and DT=30.5+/-2.7 uK at 10 degrees angular resolution. The 53X90 GHz cross-correlation amplitude at zero lag is C(0)^{1/2}=36+/-5 uK (68%CL) for the unsmoothed 7 degree DMR data. A likelihood analysis of the cross correlation function, including the quadrupole anisotropy, gives a most likely quadrupole-normalized amplitude Q_{rms-PS}=12.4^{+5.2}_{-3.3} uK (68% CL) and a spectral index n=1.59^{+0.49}_{-0.55} for a power law model of initial density fluctuations, P(k)~k^n. With n fixed to 1.0 the most likely amplitude is 17.4 +/-1.5 uK (68% CL). Excluding the quadrupole anisotropy we find Q_{rms-PS}= 16.0^{+7.5}_{-5.2} uK (68% CL), n=1.21^{+0.60}_{-0.55}, and, with n fixed to 1.0 the most likely amplitude is 18.2+/-1.6 uK (68% CL). Monte Carlo simulations indicate that these derived estimates of n may be biased by ~+0.3 (with the observed low value of the quadrupole included in the analysis) and {}~+0.1 (with the quadrupole excluded). Thus the most likely bias-corrected estimate of n is between 1.1 and 1.3. Our best estimate of the dipole from the two-year DMR data is 3.363+/-0.024 mK towards Galactic coordinates (l,b)= (264.4+/-0.2 degrees, +48.1+/-0.4 degrees), and our best estimate of the RMS quadrupole amplitude in our sky is 6+/-3 uK.

The angular power spectrum estimator developed by Peebles (1973) and Hauser & Peebles (1973) has been modified and applied to the 2 year maps produced by the COBE DMR. The power spectrum of the real sky has been compared to the power spectra of a large number of simulated random skies produced with noise equal to the observed noise and primordial density fluctuation power spectra of power law form, with $P(k) \propto k^n$. Within the limited range of spatial scales covered by the COBE DMR, corresponding to spherical harmonic indices $3 \leq \ell \lsim 30$, the best fitting value of the spectral index is $n = 1.25^{+0.4}_{-0.45}$ with the Harrison-Zeldovich value $n = 1$ approximately 0.5$\sigma$ below the best fit. For $3 \leq \ell \lsim 19$, the best fit is $n = 1.46^{+0.39}_{-0.44}$. Comparing the COBE DMR $\Delta T/T$ at small $\ell$ to the $\Delta T/T$ at $\ell \approx 50$ from degree scale anisotropy experiments gives a smaller range of acceptable spectral indices which includes $n = 1$.

A new technique of Fourier analysis on a cut sky (Gorski, 1994) has been applied to the two year COBE DMR sky maps. The Bayesian power spectrum estimation results are consistent with the Harrison-Zel'dovich n=1 model. The maximum likelihood estimates of the parameters of the power spectrum of primordial perturbations are n=1.22 (1.02) and Q_{rms-PS}=17 (20) uK including (excluding) the quadrupole anisotropy. The marginal likelihood function on n renders n=1.10 \pm 0.32 (0.87 \pm 0.36).

The {\it COBE} DMR sky maps contain low-level correlated noise. We obtain estimates of the amplitude and pattern of the correlated noise from three techniques: angular averages of the covariance matrix, Monte Carlo simulations of two-point correlation functions, and direct analysis of the DMR maps. The results from the three methods are mutually consistent. The noise covariance matrix of a DMR sky map is diagonal to an accuracy of better than 1\%. For a given sky pixel, the dominant noise covariance occurs with the ring of pixels at an angular separation of $60 \deg$ due to the $60 \deg$ separation of the DMR horns. The mean covariance at $60 \deg$ is $0.45\% ^{+0.18}_{-0.14}$ of the mean variance. Additionally, the variance in a given pixel is $0.7\%$ greater than would be expected from a single beam experiment with the same noise properties. Auto-correlation functions suffer from a $\sim 1.5\; \sigma$ positive bias at $60 \deg$ while cross-correlations have no bias. Published {\it COBE} DMR results are not significantly affected by correlated noise. COBE pre-print 94-

We have searched the temperature maps from the COBE Differential Microwave Radiometers (DMR) first two years of data for evidence of unresolved sources. The high-latitude sky (|b| > 30\deg) contains no sources brighter than 192 uK thermodynamic temperature (322 Jy at 53 GHz). The cumulative count of sources brighter than threshold T, N(> T), is consistent with a superposition of instrument noise plus a scale-invariant spectrum of cosmic temperature fluctuations normalized to Qrms-PS = 17 uK. We examine the temperature maps toward nearby clusters and find no evidence for any Sunyaev-Zel'dovich effect, \Delta y < 7.3 x 10^{-6} (95% CL) averaged over the DMR beam. We examine the temperature maps near the brightest expected radio sources and detect no evidence of significant emission. The lack of bright unresolved sources in the DMR maps, taken with anisotropy measurements on smaller angular scales, places a weak constraint on the integral number density of any unresolved Planck-spectrum sources brighter than flux density S, n(> S) < 2 x 10^4 (S/1 Jy)^{-2} sr^{-1}.

We use statistical and topological quantities to test the COBE-DMR first year sky maps against the hypothesis that the observed temperature fluctuations reflect Gaussian initial density perturbations with random phases. Recent papers discuss specific quantities as discriminators between Gaussian and non-Gaussian behavior, but the treatment of instrumental noise on the data is largely ignored. The presence of noise in the data biases many statistical quantities in a manner dependent on both the noise properties and the unknown CMB temperature field. Appropriate weighting schemes can minimize this effect, but it cannot be completely eliminated. Analytic expressions are presented for these biases, and Monte Carlo simulations used to assess the best strategy for determining cosmologically interesting information from noisy data. The genus is a robust discriminator that can be used to estimate the power law quadrupole-normalized amplitude independently of the 2-point correlation function. The genus of the DMR data are consistent with Gaussian initial fluctuations with Q_rms = 15.7 +/- 2.2 - (6.6 +/- 0.3)(n - 1) uK where n is the power law index. Fitting the rms temperature variations at various smoothing angles gives Q_rms = 13.2 +/- 2.5 uK and n = 1.7 +0.3 -0.6. While consistent with Gaussian fluctuations, the first year data are only sufficient to rule out strongly non-Gaussian distributions of fluctuations.

Abbreviated abstract of an invited talk at the International Cosmic Ray Conference in Calgary, July 1993: The results from the COBE satellite are in close agreement with the predictions of the standard hot Big Bang model, suggesting that the Universe was once hot, dense and isothermal, giving a background radiation spectrum that is close to a perfect blackbody. The spectrum observed by the FIRAS instrument places strong limits on events and scenarios occurring later than 1 year after the Big Bang. The observation of intrinsic anisotropy of the microwave background by the DMR instrument provides a measurement of the magnitude of the gravitational potential fluctuations that existed at decoupling. These potential perturbations were either produced in an inflationary epoch, or else they are initial conditions dating from t = 0. The angular power spectrum of dT is compatible with the inflationary scenario. The observed amplitude of dT can be used to restrict theories of the fundamental structure of matter on micro-physical scales and at energies higher than that of the most energetic cosmic rays. The DIRBE instrument will provide greatly improved data on extragalactic IR backgrounds after the strong foreground emission is subtracted.

The FIRAS instrument on the COBE satellite has conducted an unbiased survey of the far-infrared emission from our Galaxy. The first results of this survey were reported by Wright et al. (1991). We report the results of new analyses of this spectral survey, which includes emission lines from 158 um C+, 122 um and 205 um N+, 370 um and 609 um C, and CO J=2-1 through 5-4. We report the morphological distribution along the galactic plane (b=0) of the spectral line emission, and the high galactic latitude intensities of the C+ and 205 um N+ emission. The high galactic latitude intensity cosecant of the 158 um fine structure transition from C+ is presented, and C+ is seen to decrease more rapidly than the far infrared intensity with increasing galactic latitude. C+ and H I emission are closely correlated with a C+ cooling rate of (2.65 +/- 0.15)x10^{-26} erg/s/H-atom. We conclude that this emission arises almost entirely from the Cold Neutral Medium. The high galactic latitude intensity of the 205 um fine structure transition from N+ arises entirely from the Warm Ionized Medium, and its cosecant is presented. We estimate the total ionizing photon rate in the Galaxy to be 3.5x10^{53} ionizing photons per second, based on the 205 um N+ transition.

We compute the three-point temperature correlation function of the {\it COBE} Differential Microwave Radiometer (DMR) first-year sky maps to search for non-Gaussian temperature fluctuations. The level of fluctuations seen in the computed correlation function are too large to be attributable solely to instrument noise. However the fluctuations are consistent with the level expected to result from a superposition of instrument noise and sky signal arising from a Gaussian power law model of initial fluctuations, with a quadrupole normalized amplitude of 17 $\mu$K and a power law spectral index $n = 1$. We place limits on the amplitude of intrinsic three-point correlations with a variety of predicted functional forms.