The Legacy Surveys, a combination of three ground-based imaging surveys, have mapped 16,000 deg$^2$ in three optical bands ($g$, $r$, and $z$) to a depth 1--$2$~mag deeper than the Sloan Digital Sky Survey (SDSS). Our work addresses one of the major challenges of wide-field imaging surveys conducted at ground-based observatories: the varying depth that results from varying observing conditions at Earth-bound sites. To mitigate these effects, two of the Legacy Surveys (the Dark Energy Camera Legacy Survey, or DECaLS; and the Mayall $z$-band Legacy Survey, or MzLS) employed a unique strategy to dynamically adjust the exposure times as rapidly as possible in response to the changing observing conditions. We present the tiling and observing strategies used by these surveys. We demonstrate that the tiling and dynamic observing strategies jointly result in a more uniform-depth survey that has higher efficiency for a given total observing time compared with the traditional approach of using fixed exposure times.

We develop a data-driven model to map stellar parameters (effective temperature, surface gravity and metallicity) accurately and precisely to broad-band stellar photometry. This model must, and does, simultaneously constrain the passband-specific dust reddening vector in the Milky Way. The model uses a neural network to learn the (de-reddened) absolute magnitude in one band and colors across many bands, given stellar parameters from spectroscopic surveys and parallax constraints from Gaia. To demonstrate the effectiveness of this approach, we train our model on a dataset with spectroscopic labels from LAMOST, APOGEE and GALAH, Gaia parallaxes, and optical and near-infrared photometry from Gaia, Pan-STARRS 1, 2MASS and WISE. Testing the model on these datasets leads to an excellent fit and a precise -- and by construction accurate -- prediction of the color-magnitude diagrams in many bands. This flexible approach rigorously links spectroscopic and photometric surveys, and also results in an improved, temperature-dependent reddening vector. As such, it provides a simple and accurate method for predicting photometry in stellar evolutionary models. Our model will form a basis to infer stellar properties, distances and dust extinction from photometric data, which should be of great use in 3D mapping of the Milky Way. Our trained model may be obtained at https://doi.org/10.5281/zenodo.3902382.

Probabilistic cataloging (PCAT) outperforms traditional cataloging methods on single-band optical data in crowded fields (Portillo et al. 2017). We extend our work to multiple bands, achieving greater sensitivity ($\sim$ 0.4 mag) and greater speed (500x) compared to previous single-band results. We demonstrate the effectiveness of multiband PCAT on mock data, both in terms of recovering accurate posteriors in the catalog space, and in directly deblending sources. When applied to Sloan Digital Sky Survey (SDSS) observations of M2, taking Hubble Space Telescope data as truth, our joint fit on $r$ and $i$ band data goes $\sim0.4$ mag deeper than single-band probabilistic cataloging and has a false discovery rate less than 20\% for F606W$\leq 20$. Compared to DAOPHOT, the two-band SDSS catalog fit goes nearly 1.5 magnitudes deeper using the same data, and maintains a lower false discovery rate down to F606W$\sim 20.5$. Given recent improvements in computational speed, multiband PCAT shows promise in application to large-scale surveys and is a plausible framework for joint analysis of multi-instrument observational data.

Many surveys use maximum-likelihood (ML) methods to fit models when extracting photometry from images. We show these ML estimators systematically overestimate the flux as a function of the signal-to-noise ratio and the number of model parameters involved in the fit. This bias is substantially worse for resolved sources: while a 1% bias is expected for a 10$\sigma$ point source, a 10$\sigma$ resolved galaxy with a simplified Gaussian profile suffers a 2.5% bias. This bias also behaves differently depending how multiple bands are used in the fit: simultaneously fitting all bands leads the flux bias to become roughly evenly distributed between them, while fixing the position in "non-detection" bands (i.e. forced photometry) gives flux estimates in those bands that are biased low, compounding a bias in derived colors. We show that these effects are present in idealized simulations, outputs from the Hyper Suprime-Cam fake object pipeline (SynPipe), and observations from Sloan Digital Sky Survey Stripe 82. Prescriptions to correct for the ML bias in flux, and its uncertainty, are provided.

For the past 150 years, the prevailing view of the local Interstellar Medium (ISM) was based on a peculiarity known as the Gould's Belt, an expanding ring of young stars, gas, and dust, tilted about 20$^\circ$ to the Galactic plane. Still, the physical relation between local gas clouds has remained practically unknown because the distance accuracy to clouds is of the same order or larger than their sizes. With the advent of large photometric surveys and the Gaia satellite astrometric survey this situation has changed. Here we report the 3-D structure of all local cloud complexes. We find a narrow and coherent 2.7 kpc arrangement of dense gas in the Solar neighborhood that contains many of the clouds thought to be associated with the Gould Belt. This finding is inconsistent with the notion that these clouds are part of a ring, disputing the Gould Belt model. The new structure comprises the majority of nearby star-forming regions, has an aspect ratio of about 1:20, and contains about 3 million solar masses of gas. Remarkably, the new structure appears to be undulating and its 3-D distribution is well described by a damped sinusoidal wave on the plane of the Milky Way, with an average period of about 2 kpc and a maximum amplitude of about 160 pc. Our results represent a first step in the revision of the local gas distribution and Galactic structure and offer a new, broader context to studies on the transformation of molecular gas into stars.

Accurate distances to local molecular clouds are critical for understanding the star and planet formation process, yet distance measurements are often obtained inhomogeneously on a cloud-by-cloud basis. We have recently developed a method which combines stellar photometric data with Gaia DR2 parallax measurements in a Bayesian framework to infer the distances of nearby dust clouds to a typical accuracy of $\sim5\%$. After refining the technique to target lower latitudes and incorporating deep optical data from DECam in the southern Galactic plane, we have derived a catalog of distances to molecular clouds in Reipurth (2008, Star Formation Handbook, vols I and II) which contains a large fraction of the molecular material in the solar neighborhood. Comparison with distances derived from maser parallax measurements towards the same clouds shows our method produces consistent distances with $\lesssim10\%$ scatter for clouds across our entire distance spectrum (150 pc $-$ 2.5 kpc). We hope this catalog of homogeneous distances will serve as a baseline for future work.

We present a new three-dimensional map of dust reddening, based on Gaia parallaxes and stellar photometry from Pan-STARRS 1 and 2MASS. This map covers the sky north of a declination of -30 degrees, out to a distance of several kiloparsecs. This new map contains three major improvements over our previous work. First, the inclusion of Gaia parallaxes dramatically improves distance estimates to nearby stars. Second, we incorporate a spatial prior that correlates the dust density across nearby sightlines. This produces a smoother map, with more isotropic clouds and smaller distance uncertainties, particularly to clouds within the nearest kiloparsec. Third, we infer the dust density with a distance resolution that is four times finer than in our previous work, to accommodate the improvements in signal-to-noise enabled by the other improvements. As part of this work, we infer the distances, reddenings and types of 799 million stars. We obtain typical reddening uncertainties that are ~30% smaller than those reported in the Gaia DR2 catalog, reflecting the greater number of photometric passbands that enter into our analysis. Our 3D dust map can be accessed at https://doi.org/10.7910/DVN/2EJ9TX or through the Python package "dustmaps," and can be queried interactively at http://argonaut.skymaps.info. Our catalog of stellar parameters can be accessed at https://doi.org/10.7910/DVN/AV9GXO.

We present a uniform catalog of accurate distances to local molecular clouds informed by the Gaia DR2 data release. Our methodology builds on that of Schlafly et al. (2014). First, we infer the distance and extinction to stars along sightlines towards the clouds using optical and near-infrared photometry. When available, we incorporate knowledge of the stellar distances obtained from Gaia DR2 parallax measurements. We model these per-star distance-extinction estimates as being caused by a dust screen with a 2-D morphology derived from Planck at an unknown distance, which we then fit for using a nested sampling algorithm. We provide updated distances to the Schlafly et al. (2014) sightlines towards the Dame et al. (2001) and Magnani et al. (1985) clouds, finding good agreement with the earlier work. For a subset of 27 clouds, we construct interactive pixelated distance maps to further study detailed cloud structure, and find several clouds which display clear distance gradients and/or are comprised of multiple components. We use these maps to determine robust average distances to these clouds. The characteristic combined uncertainty on our distances is approximately 5-6%, though this can be higher for clouds at farther distances, due to the limitations of our single-cloud model.

Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing.

Galactic interstellar dust has a profound impact not only on our observations of objects throughout the Universe, but also on the morphology, star formation, and chemical evolution of the Galaxy. The advent of massive imaging and spectroscopic surveys (particularly in the infrared) places us on the threshold of being able to map the properties and dynamics of dust and the interstellar medium (ISM) in three dimensions throughout the Milky Way disk and bulge. These developments will enable a fundamentally new understanding of dust properties, including how grains respond to their local environment and how those environments affect dust attenuation of background objects of interest. Distance-resolved maps of dust motion also hold great promise for tracing the flow of interstellar material throughout the Galaxy on a variety of scales, from bar-streaming motions to the collapse and dissolution of individual molecular clouds. These advances require optical and infrared imaging of stars throughout the Galactic midplane, stretching many kiloparsecs from the Sun, matched with very dense spectroscopic coverage to probe the ISM's fine-grained structure.

In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ \Delta ra, \Delta dec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.

We present the details of the photometric and astrometric calibration of the Pan-STARRS1 $3\pi$ Survey. The photometric goals were to reduce the systematic effects introduced by the camera and detectors, and to place all of the observations onto a photometric system with consistent zero points over the entire area surveyed, the ~30,000 square degrees north of $\delta$ = -30 degrees. The astrometric calibration compensates for similar systematic effects so that positions, proper motions, and parallaxes are reliable as well. The Pan-STARRS Data Release 2 (DR2) astrometry is tied to the Gaia DR1 release.

We present a new technique to determine distances to major star-forming regions across the Perseus Molecular Cloud, using a combination of stellar photometry, astrometric data, and $\rm ^{12} CO$ spectral-line maps. Incorporating the Gaia DR2 parallax measurements when available, we start by inferring the distance and reddening to stars from their Pan-STARRS1 and 2MASS photometry, based on a technique presented in Green et al. 2014; Green et al. 2015 and implemented in their 3D "Bayestar" dust map of three-quarters of the sky. We then refine the Green et al. technique by using the velocity slices of a CO spectral cube as dust templates and modeling the cumulative distribution of dust along the line of sight towards these stars as a linear combination of the emission in the slices. Using a nested sampling algorithm, we fit these per-star distance-reddening measurements to find the distances to the CO velocity slices towards each star-forming region. This results in distance estimates explicitly tied to the velocity structure of the molecular gas. We determine distances to the B5, IC348, B1, NGC1333, L1448, and L1451 star-forming regions and find that individual clouds are located between $\approx 275-300$ pc, with typical combined uncertainties of $\approx 5\%$. We find that the velocity gradient across Perseus corresponds to a distance gradient of about 25 pc, with the eastern portion of the cloud farther away than the western portion. We determine an average distance to the complex of $294\pm 17$ pc, about 60 pc higher than the distance derived to the western portion of the cloud using parallax measurements of water masers associated with young stellar objects. The method we present is not limited to the Perseus Complex, but may be applied anywhere on the sky with adequate CO data in the pursuit of more accurate 3D maps of molecular clouds in the solar neighborhood and beyond.

We present the first application of the angle-dependent 3-Point Correlation Function (3PCF) to the density fields magnetohydrodynamic (MHD) turbulence simulations intended to model interstellar (ISM) turbulence. Previous work has demonstrated that the angle-averaged bispectrum, the 3PCF's Fourier-space analog, is sensitive to the sonic and Alfv\'enic Mach numbers of turbulence. Here we show that introducing angular information via multipole moments with respect to the triangle opening angle offers considerable additional discriminatory power on these parameters. We exploit a fast, order $N_{\rm g} \log N_{\rm g}$ ($N_{\rm g}$ the number of grid cells used for a Fourier Transform) 3PCF algorithm to study a suite of MHD turbulence simulations with 10 different combinations of sonic and Alfv\'enic Mach numbers over a range from sub to super-sonic and sub to super-Alfv\'{e}nic. The 3PCF algorithm's speed for the first time enables full quantification of the time-variation of our signal: we study 9 timeslices for each condition, demonstrating that the 3PCF is sufficiently time-stable to be used as an ISM diagnostic. In future, applying this framework to 3-D dust maps will enable better treatment of dust as a cosmological foreground as well as reveal conditions in the ISM that shape star formation.

We present a color-based method for identifying red giants based on Pan-STARRS grz and WISE W1 and W2 photometry. We utilize a subsample of bright stars with precise parallaxes from Gaia DR2 to verify that the color-based selection reliably separates dwarfs from giants. The selection is conservative in the sense that contamination is small (~30%) but not all giants are included (the selection primarily identifies K giants). The color-based selection can be applied to stars brighter than $W1\approx16$, more than two magnitudes fainter than techniques relying on shallower 2MASS photometry. Many streams and clouds are visible in the resulting sky maps, especially when binned by Gaia DR2 proper motions, including the Sagittarius stream, the Hercules-Aquila Cloud, the Eastern Banded Structure, Monoceros, and the Virgo Overdensity. In addition to the characterization of new and known stellar streams, we expect that this method for selecting red giants will enable detailed analysis of the diffuse stellar halo to distances exceeding 100 kpc.

SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$\mu$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$\sigma$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx are available at http://spherex.caltech.edu. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. Though these three themes are undoubtedly compelling, they are far from exhausting the scientific output of SPHEREx. Indeed, SPHEREx would create a unique all-sky spectral database including spectra of very large numbers of astronomical and solar system targets, including both extended and diffuse sources. These spectra would enable a wide variety of investigations, and the SPHEREx team is dedicated to making the data available to the community to enable these investigations, which we refer to as Legacy Science. To that end, we have sponsored two workshops for the general scientific community to identify the most interesting Legacy Science themes and to ensure that the SPHEREx data products are responsive to their needs. In February of 2016, some 50 scientists from all fields met in Pasadena to develop these themes and to understand their implications for the SPHEREx mission. The 2016 workshop highlighted many synergies between SPHEREx and other contemporaneous astronomical missions, facilities, and databases. Consequently, in January 2018 we convened a second workshop at the Center for Astrophysics in Cambridge to focus specifically on these synergies. This white paper reports on the results of the 2018 SPHEREx workshop.

We present optical light curves, redshifts, and classifications for 365 spectroscopically confirmed Type Ia supernovae (SNe Ia) discovered by the Pan-STARRS1 (PS1) Medium Deep Survey. We detail improvements to the PS1 SN photometry, astrometry and calibration that reduce the systematic uncertainties in the PS1 SN Ia distances. We combine the subset of 279 PS1 SN Ia ($0.03 < z < 0.68$) with useful distance estimates of SN Ia from SDSS, SNLS, various low-z and HST samples to form the largest combined sample of SN Ia consisting of a total of 1048 SN Ia ranging from $0.01 < z < 2.3$, which we call the `Pantheon Sample'. When combining Planck 2015 CMB measurements with the Pantheon SN sample, we find $\Omega_m=0.307\pm0.012$ and $w = -1.026\pm0.041$ for the wCDM model. When the SN and CMB constraints are combined with constraints from BAO and local H0 measurements, the analysis yields the most precise measurement of dark energy to date: $w0 = -1.007\pm 0.089$ and $wa = -0.222 \pm0.407$ for the w0waCDM model. Tension with a cosmological constant previously seen in an analysis of PS1 and low-z SNe has diminished after an increase of $2\times$ in the statistics of the PS1 sample, improved calibration and photometry, and stricter light-curve quality cuts. We find the systematic uncertainties in our measurements of dark energy are almost as large as the statistical uncertainties, primarily due to limitations of modeling the low-redshift sample. This must be addressed for future progress in using SN Ia to measure dark energy.

We present a simple but effective technique for measuring angular variation in $R_V$ across the sky. We divide stars from the Pan-STARRS1 catalog into Healpix pixels and determine the posterior distribution of reddening and $R_V$ for each pixel using two independent Monte Carlo methods. We find the two methods to be self-consistent in the limits where they are expected to perform similarly. We also find some agreement with high-precision photometric studies of $R_V$ in Perseus and Ophiuchus, as well as with a map of reddening near the Galactic plane based on stellar spectra from APOGEE. While current studies of $R_V$ are mostly limited to isolated clouds, we have developed a systematic method for comparing $R_V$ values for the majority of observable dust. This is a proof of concept for a more rigorous Galactic reddening map.x

Strong lensing is a sensitive probe of the small-scale density fluctuations in the Universe. We implement a novel approach to modeling strongly lensed systems using probabilistic cataloging, which is a transdimensional, hierarchical, and Bayesian framework to sample from a metamodel (union of models with different dimensionality) consistent with observed photon count maps. Probabilistic cataloging allows us to robustly characterize modeling covariances within and across lens models with different numbers of subhalos. Unlike traditional cataloging of subhalos, it does not require model subhalos to improve the goodness of fit above the detection threshold. Instead, it allows the exploitation of all information contained in the photon count maps, for instance, when constraining the subhalo mass function. We further show that, by not including these small subhalos in the lens model, fixed-dimensional inference methods can significantly mismodel the data. Using a simulated Hubble Space Telescope (HST) dataset, we show that the subhalo mass function can be probed even when many subhalos in the sample catalogs are individually below the detection threshold and would be absent in a traditional catalog. With the planned Wide Field Infrared Space Telescope (WFIRST), simultaneous probabilistic cataloging of dark subhalos in high-resolution, deep strong lens images has the potential to constrain the subhalo mass function at even lower masses.

We present a new 3D map of interstellar dust reddening, covering three quarters of the sky (declinations greater than -30 degrees) out to a distance of several kiloparsecs. The map is based on high-quality stellar photometry of 800 million stars from Pan-STARRS 1 and 2MASS. We divide the sky into sightlines containing a few hundred stars each, and then infer stellar distances and types, along with the line-of-sight dust distribution. Our new map incorporates a more accurate average extinction law and an additional 1.5 years of Pan-STARRS 1 data, tracing dust to greater extinctions and at higher angular resolutions than our previous map. Out of the plane of the Galaxy, our map agrees well with 2D reddening maps derived from far-infrared dust emission. After accounting for a 15% difference in scale, we find a mean scatter of 10% between our map and the Planck far-infrared emission-based dust map, out to a depth of 0.8 mag in E(r-z), with the level of agreement varying over the sky. Our map can be downloaded at http://argonaut.skymaps.info, or by its DOI: 10.7910/DVN/LCYHJG.

Thick back-illuminated deep-depletion CCDs have superior quantum efficiency over previous generations of thinned and traditional thick CCDs. As a result, they are being used for wide-field imaging cameras in several major projects. We use observations from the Pan-STARRS $3\pi$ survey to characterize the behavior of the deep-depletion devices used in the Pan-STARRS1 Gigapixel Camera. We have identified systematic spatial variations in the photometric measurements and stellar profiles which are similar in pattern to the so-called "tree rings" identified in devices used by other wide-field cameras (e.g., DECam and Hypersuprime Camera). The tree-ring features identified in these other cameras result from lateral electric fields which displace the electrons as they are transported in the silicon to the pixel location. In contrast, we show that the photometric and morphological modifications observed in the GPC1 detectors are caused by variations in the vertical charge transportation rate and resulting charge diffusion variations.

The DECam Plane Survey is a five-band optical and near-infrared survey of the southern Galactic plane with the Dark Energy Camera at Cerro Tololo. The survey is designed to reach past the main-sequence turn-off at the distance of the Galactic center through a reddening E(B-V) of 1.5 mag. Typical single-exposure depths are 23.7, 22.8, 22.3, 21.9, and 21.0 mag in the grizY bands, with seeing around 1 arcsecond. The footprint covers the Galactic plane with |b| < 4 degrees, 5 degrees > l > -120 degrees. The survey pipeline simultaneously solves for the positions and fluxes of tens of thousands of sources in each image, delivering positions and fluxes of roughly two billion stars with better than 10 mmag precision. Most of these objects are highly reddened and deep in the Galactic disk, probing the structure and properties of the Milky Way and its interstellar medium. The full survey is publicly available.

Cataloging is challenging in crowded fields because sources are extremely covariant with their neighbors and blending makes even the number of sources ambiguous. We present the first optical probabilistic catalog, cataloging a crowded (~0.1 sources per pixel brighter than 22nd magnitude in F606W) Sloan Digital Sky Survey r band image from M2. Probabilistic cataloging returns an ensemble of catalogs inferred from the image and thus can capture source-source covariance and deblending ambiguities. By comparing to a traditional catalog of the same image and a Hubble Space Telescope catalog of the same region, we show that our catalog ensemble better recovers sources from the image. It goes more than a magnitude deeper than the traditional catalog while having a lower false discovery rate brighter than 20th magnitude. We also present an algorithm for reducing this catalog ensemble to a condensed catalog that is similar to a traditional catalog, except it explicitly marginalizes over source-source covariances and nuisance parameters. We show that this condensed catalog has a similar completeness and false discovery rate to the catalog ensemble. Future telescopes will be more sensitive, and thus more of their images will be crowded. Probabilistic cataloging performs better than existing software in crowded fields and so should be considered when creating photometric pipelines in the Large Synoptic Space Telescope era.

Near-infrared spectroscopy from APOGEE and wide-field optical photometry from Pan-STARRS1 have recently made possible precise measurements of the shape of the extinction curve for tens of thousands of stars, parameterized by R(V). These measurements revealed structures in R(V) with large angular scales, which are challenging to explain in existing dust paradigms. In this work, we combine three-dimensional maps of dust column density with R(V) measurements to constrain the three-dimensional distribution of R(V) in the Milky Way. We find that variations in R(V) are correlated on kiloparsec scales. In particular, most of the dust within one kiloparsec in the outer Galaxy, including many local molecular clouds (Orion, Taurus, Perseus, California, Cepheus), has a significantly lower R(V) than more distant dust in the Milky Way. These results provide new input to models of dust evolution and processing, and complicate application of locally derived extinction curves to more distant regions of the Milky Way and to other galaxies.

Detection of point sources in images is a fundamental operation in astrophysics, and is crucial for constraining population models of the underlying point sources or characterizing the background emission. Standard techniques fall short in the crowded-field limit, losing sensitivity to faint sources and failing to track their covariance with close neighbors. We construct a Bayesian framework to perform inference of faint or overlapping point sources. The method involves probabilistic cataloging, where samples are taken from the posterior probability distribution of catalogs consistent with an observed photon count map. In order to validate our method we sample random catalogs of the gamma-ray sky in the direction of the North Galactic Pole (NGP) by binning the data in energy and Point Spread Function (PSF) classes. Using three energy bins spanning $0.3 - 1$, $1 - 3$ and $3 - 10$ GeV, we identify $270\substack{+30 \\ -10}$ point sources inside a $40^\circ \times 40^\circ$ region around the NGP above our point-source inclusion limit of $3 \times 10^{-11}$/cm$^2$/s/sr/GeV at the $1-3$ GeV energy bin. Modeling the flux distribution as a power law, we infer the slope to be $-1.92\substack{+0.07 \\ -0.05}$ and estimate the contribution of point sources to the total emission as $18\substack{+2 \\ -2}$\%. These uncertainties in the flux distribution are fully marginalized over the number as well as the spatial and spectral properties of the unresolved point sources. This marginalization allows a robust test of whether the apparently isotropic emission in an image is due to unresolved point sources or of truly diffuse origin.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

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

We present a panoramic map of the entire Milky Way halo north of dec~-30 degrees (~30,000 deg^2), constructed by applying the matched-filter technique to the Pan-STARRS1 3Pi Survey dataset. Using single-epoch photometry reaching to g~22, we are sensitive to stellar substructures with heliocentric distances between 3.5 and ~35 kpc. We recover almost all previously-reported streams in this volume and demonstrate that several of these are significantly more extended than earlier datasets have indicated. In addition, we also report five new candidate stellar streams. One of these features appears significantly broader and more luminous than the others and is likely the remnant of a dwarf galaxy. The other four streams are consistent with a globular cluster origin, and three of these are rather short in projection (<10 degrees), suggesting that streams like Ophiuchus may not be that rare. Finally, a significant number of more marginal substructures are also revealed by our analysis; many of these features can also be discerned in matched-filter maps produced by other authors from SDSS data, and hence they are very likely to be genuine. However, the extant 3Pi data is currently too shallow to determine their properties or produce convincing CMDs. The global view of the Milky Way provided by Pan-STARRS1 provides further evidence for the important role of both globular cluster disruption and dwarf galaxy accretion in building the Milky Way's stellar halo.

SPHEREx is a proposed SMEX mission selected for Phase A. SPHEREx will carry out the first all-sky spectral survey and provide for every 6.2" pixel a spectra between 0.75 and 4.18 $\mu$m [with R$\sim$41.4] and 4.18 and 5.00 $\mu$m [with R$\sim$135]. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. It is readily apparent, however, that many other questions in astrophysics and planetary sciences could be addressed with the SPHEREx data. The SPHEREx team convened a community workshop in February 2016, with the intent of enlisting the aid of a larger group of scientists in defining these questions. This paper summarizes the rich and varied menu of investigations that was laid out. It includes studies of the composition of main belt and Trojan/Greek asteroids; mapping the zodiacal light with unprecedented spatial and spectral resolution; identifying and studying very low-metallicity stars; improving stellar parameters in order to better characterize transiting exoplanets; studying aliphatic and aromatic carbon-bearing molecules in the interstellar medium; mapping star formation rates in nearby galaxies; determining the redshift of clusters of galaxies; identifying high redshift quasars over the full sky; and providing a NIR spectrum for most eROSITA X-ray sources. All of these investigations, and others not listed here, can be carried out with the nominal all-sky spectra to be produced by SPHEREx. In addition, the workshop defined enhanced data products and user tools which would facilitate some of these scientific studies. Finally, the workshop noted the high degrees of synergy between SPHEREx and a number of other current or forthcoming programs, including JWST, WFIRST, Euclid, GAIA, K2/Kepler, TESS, eROSITA and LSST.

Inferences about the spatial density or phase-space structure of stellar populations in the Milky Way require a precise determination of the effective survey volume. The volume observed by surveys such as Gaia or near-infrared spectroscopic surveys, which have good coverage of the Galactic mid-plane region, is highly complex because of the abundant small-scale structure in the three-dimensional interstellar dust extinction. We introduce a novel framework for analyzing the importance of small-scale structure in the extinction. This formalism demonstrates that the spatially-complex effect of extinction on the selection function of a pencil-beam or contiguous sky survey is equivalent to a low-pass filtering of the extinction-affected selection function with the smooth density field. We find that the angular resolution of current 3D extinction maps is sufficient for analyzing Gaia sub-samples of millions of stars. However, the current distance resolution is inadequate and needs to be improved by an order of magnitude, especially in the inner Galaxy. We also present a practical and efficient method for properly taking the effect of extinction into account in analyses of Galactic structure through an effective selection function. We illustrate its use with the selection function of red-clump stars in APOGEE using and comparing a variety of current 3D extinction maps.

The dust extinction curve is a critical component of many observational programs and an important diagnostic of the physics of the interstellar medium. Here we present new measurements of the dust extinction curve and its variation towards tens of thousands of stars, a hundred-fold larger sample than in existing detailed studies. We use data from the APOGEE spectroscopic survey in combination with ten-band photometry from Pan-STARRS1, 2MASS, and WISE. We find that the extinction curve in the optical through infrared is well characterized by a one-parameter family of curves described by R(V). The extinction curve is more uniform than suggested in past works, with sigma(R(V)) = 0.18, and with less than one percent of sight lines having R(V) > 4. Our data and analysis have revealed two new aspects of Galactic extinction: first, we find significant, wide-area variations in R(V) throughout the Galactic plane. These variations are on scales much larger than individual molecular clouds, indicating that R(V) variations must trace much more than just grain growth in dense molecular environments. Indeed, we find no correlation between R(V) and dust column density up to E(B-V) ~ 2. Second, we discover a strong relationship between R(V) and the far-infrared dust emissivity.

We present a recalibration of the Sloan Digital Sky Survey (SDSS) photometry with new flat fields and zero points derived from Pan-STARRS1 (PS1). Using PSF photometry of 60 million stars with $16 < r < 20$, we derive a model of amplifier gain and flat-field corrections with per-run RMS residuals of 3 millimagnitudes (mmag) in $griz$ bands and 15 mmag in $u$ band. The new photometric zero points are adjusted to leave the median in the Galactic North unchanged for compatibility with previous SDSS work. We also identify transient non-photometric periods in SDSS ("contrails") based on photometric deviations co-temporal in SDSS bands. The recalibrated stellar PSF photometry of SDSS and PS1 has an RMS difference of {9,7,7,8} mmag in $griz$, respectively, when averaged over $15'$ regions.

Current cosmological analyses which use Type Ia supernova (SN Ia) observations combine SN samples to expand the redshift range beyond that of a single sample and increase the overall sample size. The inhomogeneous photometric calibration between different SN samples is one of the largest systematic uncertainties of the cosmological parameter estimation. To place these different samples on a single system, analyses currently use observations of a small sample of very bright flux standards on the $HST$ system. We propose a complementary method, called `Supercal', in which we use measurements of secondary standards in each system, compare these to measurements of the same stars in the Pan-STARRS1 (PS1) system, and determine offsets for each system relative to PS1, placing all SN observations on a single, consistent photometric system. PS1 has observed $3\pi$ of the sky and has a relative calibration of better than 5 mmag (for $\sim15<griz<21$ mag), making it an ideal reference system. We use this process to recalibrate optical observations taken by the following SN samples: PS1, SNLS, SDSS, CSP, and CfA1-4. We measure discrepancies on average of 10 mmag, but up to 35 mmag, in various optical passbands. We find that correcting for these differences changes recovered values for the dark energy equation-of-state parameter, $w$, by on average $2.6\%$. This change is roughly half the size of current statistical constraints on $w$. The size of this effect strongly depends on the error in the $B-V$ calibration of the low-$z$ surveys. The Supercal method will allow future analyses to tie past samples to the best calibrated sample.

We present a three-dimensional map of interstellar dust reddening, covering three-quarters of the sky out to a distance of several kiloparsecs, based on Pan-STARRS 1 and 2MASS photometry. The map reveals a wealth of detailed structure, from filaments to large cloud complexes. The map has a hybrid angular resolution, with most of the map at an angular resolution of 3.4' to 13.7', and a maximum distance resolution of ~25%. The three-dimensional distribution of dust is determined in a fully probabilistic framework, yielding the uncertainty in the reddening distribution along each line of sight, as well as stellar distances, reddenings and classifications for 800 million stars detected by Pan-STARRS 1. We demonstrate the consistency of our reddening estimates with those of two-dimensional emission-based maps of dust reddening. In particular, we find agreement with the Planck 353 GHz optical depth-based reddening map to within 0.05 mag in E(B-V) to a depth of 0.5 mag, and explore systematics at reddenings less than E(B-V) ~ 0.08 mag. We validate our per-star reddening estimates by comparison with reddening estimates for stars with both SDSS photometry and SEGUE spectral classifications, finding per-star agreement to within 0.1 mag out to a stellar E(B-V) of 1 mag. We compare our map to two existing three-dimensional dust maps, by Marshall et al. (2006) and Lallement et al. (2013), demonstrating our finer angular resolution, and better distance resolution compared to the former within ~3 kpc. The map can be queried or downloaded at http://argonaut.skymaps.info. We expect the three-dimensional reddening map presented here to find a wide range of uses, among them correcting for reddening and extinction for objects embedded in the plane of the Galaxy, studies of Galactic structure, calibration of future emission-based dust maps and determining distances to objects of known reddening.

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

Past studies have identified a spatially extended excess of $\sim$1-3 GeV gamma rays from the region surrounding the Galactic Center, consistent with the emission expected from annihilating dark matter. We revisit and scrutinize this signal with the intention of further constraining its characteristics and origin. By applying cuts to the \textit{Fermi} event parameter CTBCORE, we suppress the tails of the point spread function and generate high resolution gamma-ray maps, enabling us to more easily separate the various gamma-ray components. Within these maps, we find the GeV excess to be robust and highly statistically significant, with a spectrum, angular distribution, and overall normalization that is in good agreement with that predicted by simple annihilating dark matter models. For example, the signal is very well fit by a 36-51 GeV dark matter particle annihilating to $b\bar{b}$ with an annihilation cross section of $\sigma v = (1-3)\times 10^{-26}$ cm$^3$/s (normalized to a local dark matter density of 0.4 GeV/cm$^3$). Furthermore, we confirm that the angular distribution of the excess is approximately spherically symmetric and centered around the dynamical center of the Milky Way (within $\sim$$0.05^{\circ}$ of Sgr A$^*$), showing no sign of elongation along the Galactic Plane. The signal is observed to extend to at least $\simeq10^{\circ}$ from the Galactic Center, disfavoring the possibility that this emission originates from millisecond pulsars.

The Orion Molecular Complex is the nearest site of ongoing high-mass star formation, making it one of the most extensively studied molecular complexes in the Galaxy. We have developed a new technique for mapping the 3D distribution of dust in the Galaxy using Pan-STARRS1 photometry. We isolate the dust at the distance to Orion using this technique, revealing a large (100 pc, 14 degree diameter), previously unrecognized ring of dust, which we term the "Orion dust ring." The ring includes Orion A and B, and is not coincident with current H-alpha features. The circular morphology suggests formation as an ancient bubble in the interstellar medium, though we have not been able to conclusively identify the source of the bubble. This hint at the history of Orion may have important consequences for models of high-mass star formation and triggered star formation.

We apply the Finkbeiner et al. (1999) two-component thermal dust emission model to the Planck HFI maps. This parametrization of the far-infrared dust spectrum as the sum of two modified blackbodies serves as an important alternative to the commonly adopted single modified blackbody (MBB) dust emission model. Analyzing the joint Planck/DIRBE dust spectrum, we show that two-component models provide a better fit to the 100-3000 GHz emission than do single-MBB models, though by a lesser margin than found by Finkbeiner et al. (1999) based on FIRAS and DIRBE. We also derive full-sky 6.1' resolution maps of dust optical depth and temperature by fitting the two-component model to Planck 217-857 GHz along with DIRBE/IRAS 100 micron data. Because our two-component model matches the dust spectrum near its peak, accounts for the spectrum's flattening at millimeter wavelengths, and specifies dust temperature at 6.1' FWHM, our model provides reliable, high-resolution thermal dust emission foreground predictions from 100 to 3000 GHz. We find that, in diffuse sky regions, our two-component 100-217 GHz predictions are on average accurate to within 2.2%, while extrapolating the Planck Collaboration (2013a) single-MBB model systematically underpredicts emission by 18.8% at 100 GHz, 12.6% at 143 GHz and 7.9% at 217 GHz. We calibrate our two-component optical depth to reddening, and compare with reddening estimates based on stellar spectra. We find the dominant systematic problems in our temperature/reddening maps to be zodiacal light on large angular scales and the cosmic infrared background anisotropy on small angular scales.

The Large Area Telescope on the Fermi Gamma-ray Space Telescope has a point spread function with large tails, consisting of events affected by tracker ineffiencies, inactive volumes, and hard scattering; these tails can make source confusion a limiting factor. The parameter CTBCORE, available in the publicly available Extended Fermi LAT data, estimates the quality of each event's direction reconstruction; by implementing a cut in this parameter, the tails of the point spread function can be suppressed at the cost of losing effective area. We implement cuts on CTBCORE and present updated instrument response functions derived from the Fermi LAT data itself, along with all-sky maps generated with these cuts. Having shown the effectiveness of these cuts, especially at low energies, we encourage their use in analyses where angular resolution is more important than Poisson noise.

We present a map of the dust reddening to 4.5 kpc derived from Pan-STARRS1 stellar photometry. The map covers almost the entire sky north of declination -30 degrees at a resolution of 7' to 14', and is based on the estimated distances and reddenings to more than 500 million stars. The technique is designed to map dust in the Galactic plane, where many other techniques are stymied by the presence of multiple dust clouds at different distances along each line of sight. This reddening-based dust map agrees closely with the Schlegel, Finkbeiner, and Davis (SFD; 1998) far-infrared emission-based dust map away from the Galactic plane, and the most prominent differences between the two maps stem from known limitations of SFD in the plane. We also compare the map with Planck, finding likewise good agreement in general at high latitudes. The use of optical data from Pan-STARRS1 yields reddening uncertainty as low as 25 mmag E(B-V).

Distance measurements to molecular clouds are important, but are often made separately for each cloud of interest, employing very different different data and techniques. We present a large, homogeneous catalog of distances to molecular clouds, most of which are of unprecedented accuracy. We determine distances using optical photometry of stars along lines of sight toward these clouds, obtained from PanSTARRS-1. We simultaneously infer the reddenings and distances to these stars, tracking the full probability distribution function using a technique presented in Green et al. (2014). We fit these star-by-star measurements using a simple dust screen model to find the distance to each cloud. We thus estimate the distances to almost all of the clouds in the Magnani et al. (1985) catalog, as well as many other well-studied clouds, including Orion, Perseus, Taurus, Cepheus, Polaris, California, and Monoceros R2, avoiding only the inner Galaxy. Typical statistical uncertainties in the distances are 5%, though the systematic uncertainty stemming from the quality of our stellar models is about 10%. The resulting catalog is the largest catalog of accurate, directly-measured distances to molecular clouds. Our distance estimates are generally consistent with available distance estimates from the literature, though in some cases the literature estimates are off by a factor of more than two.

The eXciting Dark Matter (XDM) model was proposed as a mechanism to efficiently convert the kinetic energy (in sufficiently hot environments) of dark matter into e+e- pairs. The standard scenario invokes a doublet of nearly degenerate DM states, and a dark force to mediate a large upscattering cross section between the two. For heavy ($\sim TeV$) DM, the kinetic energy of WIMPs in large (galaxy-sized or larger) halos is capable of producing low-energy positrons. For lighter dark matter, this is kinematically impossible, and the unique observable signature becomes an X-ray line, arising from $\chi \chi \rightarrow \chi^* \chi^*$, followed by $\chi^* \rightarrow \chi \gamma$. This variant of XDM is distinctive from other DM X-ray scenarios in that it tends to be most present in more massive, hotter environments, such as clusters, rather than nearby dwarfs, and has different dependencies from decaying models. We find that it is capable of explaining the recently reported X-ray line at 3.56 keV. For very long lifetimes of the excited state, primordial decays can explain the signal without the presence of upscattering. Thermal models freeze-out as in the normal XDM setup, via annihilations to the light boson $\phi$. For suitable masses the annihilation $\chi \chi \rightarrow \phi \phi$ followed by $\phi \rightarrow SM$ can explain the reported gamma-ray signature from the galactic center. Direct detection is discussed, including the possibility of explaining DAMA via the "Luminous" dark matter approach. Quite generally, the proximity of the 3.56 keV line to the energy scale of DAMA motivates a reexamination of electromagnetic explanations. Other signals, including lepton jets and the modification of cores of dwarf galaxies are also considered.

We present a method to infer reddenings and distances to stars, based only on their broad-band photometry, and show how this method can be used to produce a three-dimensional dust map of the Galaxy. Our method samples from the full probability density function of distance, reddening and stellar type for individual stars, as well as the full uncertainty in reddening as a function of distance in the 3D dust map. We incorporate prior knowledge of the distribution of stars in the Galaxy and the detection limits of the survey. For stars in the Pan-STARRS 1 (PS1) 3 pi survey, we demonstrate that our reddening estimates are unbiased, and accurate to ~0.13 mag in E(B-V) for the typical star. Based on comparisons with mock catalogs, we expect distances for main-sequence stars to be constrained to within ~20% - 60%, although this range can vary, depending on the reddening of the star, the precise stellar type and its position on the sky. A further paper will present a 3D map of dust over the three quarters of the sky surveyed by PS1. Both the individual stellar inferences and the 3D dust map will enable a wealth of Galactic science in the plane. The method we present is not limited to the passbands of the PS1 survey, but may be extended to incorporate photometry from other surveys, such as 2MASS, SDSS (where available), and in the future, LSST and Gaia.

We describe our custom processing of the entire Wide-field Infrared Survey Explorer (WISE) 12 micron imaging data set, and present a high-resolution, full-sky map of diffuse Galactic dust emission that is free of compact sources and other contaminating artifacts. The principal distinctions between our resulting co-added images and the WISE Atlas stacks are our removal of compact sources, including their associated electronic and optical artifacts, and our preservation of spatial modes larger than 1.5 degrees. We provide access to the resulting full-sky map via a set of 430 12.5 degree by 12.5 degree mosaics. These stacks have been smoothed to 15" resolution and are accompanied by corresponding coverage maps, artifact images, and bit-masks for point sources, resolved compact sources, and other defects. When combined appropriately with other mid-infrared and far-infrared data sets, we expect our WISE 12 micron co-adds to form the basis for a full-sky dust extinction map with angular resolution several times better than Schlegel et al. (1998).

We present the discovery of a new dwarf galaxy, Perseus I/Andromeda XXXIII, found in the vicinity of Andromeda (M31) in stacked imaging data from the Pan-STARRS1 3{\pi} survey. Located 27.9{\deg} away from M31, Perseus I has a heliocentric distance of 785 +/- 65 kpc, compatible with it being a satellite of M31 at 374 +14/-10 kpc from its host. The properties of Perseus I are typical for a reasonably bright dwarf galaxy (M_V = -10.3 +/- 0.7), with an exponential half-light radius of r_h = 1.7 +/- 0.4 arcminutes or r_h = 400 +105/-85 pc at this distance, and a moderate ellipticity (\epsilon = 0.43 +0.15/-0.17). The late discovery of Perseus I is due to its fairly low surface brightness (\mu_0=25.7 +1.0/-0.9 mag/arcsec^2), and to the previous lack of deep, high quality photometric data in this region. If confirmed to be a companion of M31, the location of Perseus I, far east from its host, could place interesting constraints on the bulk motion of the satellite system of M31.

We report the discovery of two new dwarf galaxies, Lacerta I/Andromeda XXXI (Lac I/And XXXI) and Cassiopeia III/Andromeda XXXII (Cas III/And XXXII), in stacked Pan-STARRS1 r_P1- and i_P1-band imaging data. Both are luminous systems (M_V ~ -12) located at projected distances of 20.3{\deg} and 10.5{\deg} from M31. Lac I and Cas III are likely satellites of the Andromeda galaxy with heliocentric distances of 756^{+44}_{-28} kpc and 772^{+61}_{-56} kpc, respectively, and corresponding M31-centric distances of 275+/-7 kpc and 144^{+6}_{-4} kpc . The brightest of recent Local Group member discoveries, these two new dwarf galaxies owe their late discovery to their large sizes (r_h = 4.2^{+0.4}_{-0.5} arcmin or 912^{+124}_{-93} pc for Lac I; r_h = 6.5^{+1.2}_{-1.0} arcmin or 1456+/-267 pc for Cas III), and consequently low surface brightness (\mu_0 ~ 26.0 mag/arcsec^2), as well as to the lack of a systematic survey of regions at large radii from M31, close to the Galactic plane. This latter limitation is now alleviated by the 3{\pi} Pan-STARRS1 survey, which could lead to the discovery of other distant Andromeda satellite dwarf galaxies.

Evidence for a spectral line in the inner Galaxy has caused a great deal of excitement over the last year, mainly because of its interpretation as a possible dark matter signal. The observation has raised important questions about statistics and suspicions about systematics, especially in photons from the Earth limb. With enough additional data, we can address these concerns. In this white paper, we summarize the current observational situation and project future sensitivities, finding that the status quo is dangerously close to leaving the issue unresolved until 2015. We advocate a change in survey strategy that more than doubles the data rate in the inner Galaxy, and is relatively non-disruptive to other survey science. This strategy will clearly separate the null hypothesis from the line signal hypothesis and provide ample limb data for systematics checks by the end of 2014. The standard survey mode may not.

As of 2012 Jan 21, the Pan-STARRS 1 3\pi Survey has observed the 3/4 of the sky visible from Hawaii with a minimum of 2 and mean of 7.6 observations in 5 filters, g_p1, r_p1, i_p1, z_p1, y_p1. Now at the end of the second year of the mission, we are in a position to make an initial public release of a portion of this unprecedented dataset. This article describes the PS1 Photometric Ladder, Release 12.01 This is the first of a series of data releases to be generated as the survey coverage increases and the data analysis improves. The Photometric Ladder has rungs every hour in RA and at 4 intervals in declination. We will release updates with increased area coverage (more rungs) from the latest dataset until the PS1 survey and the final re-reduction are completed. The currently released catalog presents photometry of \approx 1000 objects per square degree in the rungs of the ladder. Saturation occurs at g_p1, r_p1, i_p1 \approx 13.5; z_p1 \approx 13.0; and y_p1 \approx 12.0. Photometry is provided for stars down to g_p1, r_p1, i_p1 \approx 19.1 in the AB system. This data release depends on the rigid `Ubercal' photometric calibration using only the photometric nights, with systematic uncertainties of (8.0, 7.0, 9.0, 10.7, 12.4) millimags in (g_p1, r_p1, i_p1, z_p1, y_p1). Areas covered only with lower quality nights are also included, and have been tied to the Ubercal solution via relative photometry; photometric accuracy of the non-photometric regions is lower and should be used with caution.

We use data from the Pan-STARRS1 survey to present a panoramic view of the Sagittarius tidal stream in the southern Galactic hemisphere. As a result of the extensive sky coverage of Pan-STARRS1, the southern stream is visible along more than 60 degrees of its orbit, nearly double the length seen by the SDSS. The recently discovered southern bifurcation of the stream is also apparent, with the fainter branch of the stream visible over at least 30 degrees. Using a combination of fitting both the main sequence turn-off and the red clump, we measure the distance to both arms of the stream in the south. We find that the distances to the bright arm of the stream agree very well with the N-body models of Law & Majewski (2010). We also find that the faint arm lies ~5 kpc closer to the Sun than the bright arm, similar to the behavior seen in the northern hemisphere.

Our recent claims of a Galactic center feature in Fermi-LAT data at approximately 130 GeV have prompted an avalanche of papers proposing explanations ranging from dark matter annihilation to exotic pulsar winds. Because of the importance of such interpretations for physics and astrophysics, a discovery will require not only additional data, but a thorough investigation of possible LAT systematics. While we do not have access to the details of each event reconstruction, we do have information about each event from the public event lists and spacecraft parameter files. These data allow us to search for suspicious trends that could indicate a spurious signal. We consider several hypotheses that might make an instrumental artifact more apparent at the Galactic center, and find them implausible. We also search for an instrumental signature in the Earth limb photons, which provide a smooth reference spectrum for null tests. We find no significant 130 GeV feature in the Earth limb sample. However, we do find a marginally significant 130 GeV feature in Earth limb photons with a limited range of detector incidence angles. This raises concerns about the 130 GeV Galactic center feature, even though we can think of no plausible model of instrumental behavior that connects the two. A modest amount of additional limb data would tell us if the limb feature is a statistical fluke. If the limb feature persists, it would raise doubts about the Pass 7 processing of E > 100 GeV events. At present we find no instrumental systematics that could plausibly explain the excess Galactic center emission at 130 GeV.

We present a precise photometric calibration of the first 1.5 years of science imaging from the Pan-STARRS1 survey (PS1), an ongoing optical survey of the entire sky north of declination -30 degrees in five bands. Building on the techniques employed by Padmanabhan et al. (2008) in the Sloan Digital Sky Survey (SDSS), we use repeat PS1 observations of stars to perform the relative calibration of PS1 in each of its five bands, solving simultaneously for the system throughput, the atmospheric transparency, and the large-scale detector flat field. Both internal consistency tests and comparison against the SDSS indicate that we achieve relative precision of <10 mmag in g, r, and i_P1, and ~10 mmag in z and y_P1. The spatial structure of the differences with the SDSS indicates that errors in both the PS1 and SDSS photometric calibration contribute similarly to the differences. The analysis suggests that both the PS1 system and the Haleakala site will enable <1% photometry over much of the sky.

Gamma ray emission from dark matter subhalos in the Milky Way has long been sought as a sign of dark matter particle annihilation. So far, searches for gamma-ray continuum from subhalos have been unsuccessful, and line searches are difficult without prior knowledge of the line energies. Guided by recent claims of line emission at 111 GeV and 129 GeV in the Galactic center, we examine the co-added gamma-ray spectrum of unassociated point sources in the Second Fermi-LAT catalog (2FGL) using 3.9 years of LAT data. Using the SOURCE event class, we find evidence for lines at 111 GeV and 129 GeV with a local significance of 3.3 sigma based on a conservative estimate of the background at E>135 GeV. Other 2FGL sources analyzed in the same way do not show line emission at 111 GeV and 129 GeV. The line-emitting sources are mostly within 30 degrees of the Galactic plane, although this anisotropy may be a selection effect. If the double-line emission from these objects is confirmed with future data, it will provide compelling support for the hypothesis that the Galactic center line signal is indeed from dark matter annihilation.

Using 3.7 years of \Fermi-LAT data, we examine the diffuse 80-200 GeV emission in the inner Galaxy and find a resolved gamma-ray feature at $\sim 110-140$ GeV. We model the spatial distribution of this emission with a $\sim3\degree$ FWHM Gaussian, finding a best fit position $1.5\degree$ West of the Galactic Center. Even better fits are obtained for off-center Einasto and power-law profiles, which are preferred over the null (no line) hypothesis by $6.5\sigma$ ($5.0\sigma/5.4\sigma$ after trials factor correction for one/two line case) assuming an NFW density profile centered at $(\ell, b)=(-1.5\degree,0\degree)$ with a power index $\alpha=1.2$ . The energy spectrum of this structure is consistent with a single spectral line (at energy $127.0\pm 2.0$ GeV with $\chi^2=4.48$ for 4 d.o.f.). A pair of lines at $110.8\pm 4.4$ GeV and $128.8\pm 2.7$ GeV provides a marginally better fit (with $\chi^2=1.25$ for 2 d.o.f.). The total luminosity of the structure is $(3.2\pm0.6)\times 10^{35}$ erg/s, or $(1.7\pm0.4)\times 10^{36}$ photons/sec. The energies in the two-line case are compatible with a $127.3 \pm 2.7$ GeV WIMP annihilating through $\gamma \gamma$ and $\gamma Z$ (with $\chi^2=1.67$ for 3 d.o.f.). We describe a possible change to the \Fermi\ scan strategy that would accumulate S/N on spectral lines in the Galactic center 4 times as fast as the current survey strategy.

Although accretion onto supermassive black holes in other galaxies is seen to produce powerful jets in X-ray and radio, no convincing detection has ever been made of a kpc-scale jet in the Milky Way. The recently discovered pair of 10 kpc tall gamma-ray bubbles in our Galaxy may be a sign of earlier jet activity from the central black hole. In this paper, we identify a gamma-ray cocoon feature in the southern bubble, a jet-like feature along the cocoon's axis of symmetry, and another directly opposite the Galactic center in the north. Both the cocoon and jet-like feature have a hard spectrum with spectral index ~ -2 from 1 to 100 GeV, with a cocoon luminosity of (5.5 +/- 0.45) x 10^35 erg/s and luminosity of the jet-like feature of (1.8 +/- 0.35) x 10^35 erg/s at 1 to 100 GeV. If confirmed, these jets are the first resolved gamma-ray jets ever seen.

We present low-resolution MMT Hectospec spectroscopy of 594 candidate Monoceros stream member stars. Based on strong color-magnitude diagram overdensities, we targeted three fields within the stream's footprint, with 178 deg < l < 203 deg and -25 deg < b < 25 deg. By comparing the measured iron abundances with those expected from smooth Galactic components alone, we measure, for the first time, the spectroscopic metallicity distribution function for Monoceros. We find the stream to be chemically distinct from both the thick disk and halo, with [Fe/H] = -1, and do not detect a trend in the stream's metallicity with Galactic longitude. Passing from b = +25 deg to b = -25 deg the median Monoceros metallicity trends upward by 0.1 dex, though uncertainties in modeling sample contamination by the disk and halo make this a marginal detection. In each field, we find Monoceros to have an intrinsic [Fe/H] dispersion of 0.10-0.22 dex. From the CaII K line, we measure [Ca/Fe] for a subsample of metal poor program stars with -1.1 < [Fe/H] < -0.5. In two of three fields, we find calcium deficiencies qualitatively similar to previously reported [Ti/Fe] underabundances in Monoceros and the Sagittarius tidal stream. Further, using 90 spectra of thick disk stars in the Monoceros pointings with b ~ +/-25 deg, we detect a 0.22 dex north/south metallicity asymmetry coincident with known stellar density asymmetry at R_GC ~ 12 kpc and |Z| ~ 1.7 kpc. Our median Monoceros [Fe/H] = -1.0 and its relatively low dispersion naturally fit the expectation for an appropriately luminous M_V ~ -13 dwarf galaxy progenitor.

We describe the design and data sample from the DEEP2 Galaxy Redshift Survey, the densest and largest precision-redshift survey of galaxies at z ~ 1 completed to date. The survey has conducted a comprehensive census of massive galaxies, their properties, environments, and large-scale structure down to absolute magnitude M_B = -20 at z ~ 1 via ~90 nights of observation on the DEIMOS spectrograph at Keck Observatory. DEEP2 covers an area of 2.8 deg^2 divided into four separate fields, observed to a limiting apparent magnitude of R_AB=24.1. Objects with z < 0.7 are rejected based on BRI photometry in three of the four DEEP2 fields, allowing galaxies with z > 0.7 to be targeted ~2.5 times more efficiently than in a purely magnitude-limited sample. Approximately sixty percent of eligible targets are chosen for spectroscopy, yielding nearly 53,000 spectra and more than 38,000 reliable redshift measurements. Most of the targets which fail to yield secure redshifts are blue objects that lie beyond z ~ 1.45. The DEIMOS 1200-line/mm grating used for the survey delivers high spectral resolution (R~6000), accurate and secure redshifts, and unique internal kinematic information. Extensive ancillary data are available in the DEEP2 fields, particularly in the Extended Groth Strip, which has evolved into one of the richest multiwavelength regions on the sky. DEEP2 surpasses other deep precision-redshift surveys at z ~ 1 in terms of galaxy numbers, redshift accuracy, sample number density, and amount of spectral information. We also provide an overview of the scientific highlights of the DEEP2 survey thus far. This paper is intended as a handbook for users of the DEEP2 Data Release 4, which includes all DEEP2 spectra and redshifts, as well as for the publicly-available DEEP2 DEIMOS data reduction pipelines. [Abridged]

We use SDSS photometry of 73 million stars to simultaneously obtain best-fit main-sequence stellar energy distribution (SED) and amount of dust extinction along the line of sight towards each star. Using a subsample of 23 million stars with 2MASS photometry, whose addition enables more robust results, we show that SDSS photometry alone is sufficient to break degeneracies between intrinsic stellar color and dust amount when the shape of extinction curve is fixed. When using both SDSS and 2MASS photometry, the ratio of the total to selective absorption, $R_V$, can be determined with an uncertainty of about 0.1 for most stars in high-extinction regions. These fits enable detailed studies of the dust properties and its spatial distribution, and of the stellar spatial distribution at low Galactic latitudes. Our results are in good agreement with the extinction normalization given by the Schlegel et al. (1998, SFD) dust maps at high northern Galactic latitudes, but indicate that the SFD extinction map appears to be consistently overestimated by about 20% in the southern sky, in agreement with Schlafly et al. (2010). The constraints on the shape of the dust extinction curve across the SDSS and 2MASS bandpasses support the models by Fitzpatrick (1999) and Cardelli et al. (1989). For the latter, we find an $R_V=3.0\pm0.1$(random) $\pm0.1$(systematic) over most of the high-latitude sky. At low Galactic latitudes (|b|<5), we demonstrate that the SFD map cannot be reliably used to correct for extinction as most stars are embedded in dust, rather than behind it. We introduce a method for efficient selection of candidate red giant stars in the disk, dubbed "dusty parallax relation", which utilizes a correlation between distance and the extinction along the line of sight. We make these best-fit parameters, as well as all the input SDSS and 2MASS data, publicly available in a user-friendly format.

High-precision measurements of the temperature and polarization anisotropies of the cosmic microwave background radiation have been previously employed to set robust constraints on dark matter annihilation during recombination. In this work we improve and generalize these constraints to apply to energy deposition during the recombination era with arbitrary redshift dependence. Our approach also provides more rigorous and model-independent bounds on dark matter annihilation and decay scenarios. We employ principal component analysis to identify a basis of weighting functions for the energy deposition. The coefficients of these weighting functions parameterize any energy deposition model and can be constrained directly by experiment. For generic energy deposition histories that are currently allowed by WMAP7 data, up to 3 principal component coefficients are measurable by Planck and up to 5 coefficients are measurable by an ideal cosmic variance limited experiment. For WIMP dark matter, our analysis demonstrates that the effect on the CMB is described well by a single (normalization) parameter and a "universal" redshift dependence for the energy deposition history. We give WMAP 7 constraints on both generic energy deposition histories and the universal WIMP case.

We present measurements of dust reddening using the colors of stars with spectra in the Sloan Digital Sky Survey. We measure reddening as the difference between the measured and predicted colors of a star, as derived from stellar parameters from the SEGUE Stellar Parameter Pipeline (Lee et al. 2008a). We achieve uncertainties of 56, 34, 25, and 29 mmag in the colors u-g, g-r, r-i, and i-z, per star, though the uncertainty varies depending on the stellar type and the magnitude of the star. The spectrum-based reddening measurements confirm our earlier "blue tip" reddening measurements (Schlafly et al. 2010, S10), finding reddening coefficients different by -3%, 1%, 1%, and 2% in u-g, g-r, r-i, and i-z from those found by the blue tip method, after removing a 4% normalization difference. These results prefer an R_V=3.1 Fitzpatrick (1999) reddening law to O'Donnell (1994) or Cardelli et al. (1989) reddening laws. We provide a table of conversion coefficients from the Schlegel et al. (1998) maps of E(B-V) to extinction in 88 bandpasses for 4 values of R_V, using this reddening law and the 14% recalibration of SFD first reported by S10 and confirmed in this work.

The magnetic inelastic dark matter (MiDM) model, in which dark matter inelastically scatters off nuclei through a magnetic dipole interaction, has previously been shown to reconcile the DAMA/LIBRA annual modulation signal with null results from other experiments. In this work, we explore the unique directional detection signature of MiDM. After the dark matter scatters into its excited state, it decays with a lifetime of order 1 microsecond and emits a photon with energy ~100 keV. Both the nuclear recoil and the corresponding emitted photon can be detected by studying delayed coincidence events. The recoil track and velocity of the excited state can be reconstructed from the nuclear interaction vertex and the photon decay vertex. The angular distribution of the WIMP recoil tracks is sharply peaked and modulates daily. It is therefore possible to observe the directional modulation of WIMP-nucleon scattering without a large-volume gaseous directional detection experiment. Furthermore, current experiments such as XENON100 can immediately measure this directional modulation and constrain the MiDM parameter space with an exposure of a few thousand kg day.

Anomalies in direct and indirect detection have motivated models of dark matter consisting of a multiplet of nearly-degenerate states, coupled by a new GeV-scale interaction. We perform a careful analysis of the thermal freezeout of dark matter annihilation in such a scenario. We compute the range of "boost factors" arising from Sommerfeld enhancement in the local halo for models which produce the correct relic density, and show the effect of including constraints on the saturated enhancement from the cosmic microwave background (CMB). We find that boost factors from Sommerfeld enhancement of up to ~800 are possible in the local halo. When the CMB bounds on the saturated enhancement are applied, the maximal boost factor is reduced to ~400 for 1-2 TeV dark matter and sub-GeV force carriers, but remains large enough to explain the observed Fermi and PAMELA electronic signals. We describe regions in the DM mass-boost factor plane where the cosmic ray data is well fit for a range of final states, and show that Sommerfeld enhancement alone is enough to provide the large annihilation cross sections required to fit the data, although for light mediator masses (less than ~200 MeV) there is tension with the CMB constraints in the absence of astrophysical boost factors from substructure. Additionally, we consider the circumstances under which WIMPonium formation is relevant and find for heavy WIMPs (greater than ~2 TeV) and soft-spectrum annihilation channels it can be an important consideration; we find regions with dark matter mass greater than 2.8 TeV that are consistent with the CMB bounds and have ~600-700 present-day boost factors.

Data from the Fermi-LAT reveal two large gamma-ray bubbles, extending 50 degrees above and below the Galactic center, with a width of about 40 degrees in longitude. The gamma-ray emission associated with these bubbles has a significantly harder spectrum (dN/dE ~ E^-2) than the IC emission from electrons in the Galactic disk, or the gamma-rays produced by decay of pions from proton-ISM collisions. There is no significant spatial variation in the spectrum or gamma-ray intensity within the bubbles, or between the north and south bubbles. The bubbles are spatially correlated with the hard-spectrum microwave excess known as the WMAP haze; the edges of the bubbles also line up with features in the ROSAT X-ray maps at 1.5-2 keV. We argue that these Galactic gamma-ray bubbles were most likely created by some large episode of energy injection in the Galactic center, such as past accretion events onto the central massive black hole, or a nuclear starburst in the last ~10 Myr. Dark matter annihilation/decay seems unlikely to generate all the features of the bubbles and the associated signals in WMAP and ROSAT; the bubbles must be understood in order to use measurements of the diffuse gamma-ray emission in the inner Galaxy as a probe of dark matter physics. Study of the origin and evolution of the bubbles also has the potential to improve our understanding of recent energetic events in the inner Galaxy and the high-latitude cosmic ray population.

We present measurements of reddening due to dust using the colors of stars in the Sloan Digital Sky Survey (SDSS). We measure the color of main sequence turn-off stars by finding the "blue tip" of the stellar locus: the prominent blue edge in the distribution of stellar colors. The method is sensitive to color changes of order 18, 12, 7, and 8 mmag of reddening in the colors u-g, g-r, r-i, and i-z, respectively, in regions measuring 90' by 14'. We present maps of the blue tip colors in each of these bands over the entire SDSS footprint, including the new dusty southern Galactic cap data provided by the SDSS-III. The results disfavor the best fit O'Donnell (1994) and Cardelli et al. (1989) reddening laws, but are well described by a Fitzpatrick (1999) reddening law with R_V = 3.1. The SFD dust map is found to trace the dust well, but overestimates reddening by factors of 1.4, 1.0, 1.2, and 1.4 in u-g, g-r, r-i, and i-z, largely due to the adopted reddening law. In select dusty regions of the sky, we find evidence for problems in the SFD temperature correction. A dust map normalization difference of 15% between the Galactic north and south sky may be due to these dust temperature errors.

The DAMA/LIBRA collaboration has detected an annual modulation of the recoil rate in NaI crystals with the phase expected for WIMP scattering events. This signal is dramatically inconsistent with upper limits from other experiments for elastically scattering weak-scale WIMPs. However, the results are compatible for the case of inelastic dark matter (iDM). The iDM theory, as implemented by Tucker-Smith and Weiner, constrains the WIMP to a tight contour in sigma_n-delta space, where delta is the mass difference between the ground state and excited WIMPs. An urgent priority in direct detection is to test this scenario. The crucial test of the iDM explanation of DAMA -- an "experimentum crucis" -- is an experiment with directional sensitivity, which can measure the daily modulation in direction. Because the contrast can be 100%, it is a sharper test than the much smaller annual modulation in the rate. We estimate the significance of such an experiment as a function of the WIMP mass, cross section, background rate, and other parameters. The proposed experiment severely constrains the DAMA/iDM scenario even with modest exposure (~1000 kg day) on gaseous xenon.

Recent cosmic ray, gamma ray, and microwave signals observed by Fermi, PAMELA, and WMAP indicate an unexpected primary source of e+e- at 10-1000 GeV. We fit these data to "standard backgrounds" plus a new source, assumed to be a separable function of position and energy. For the spatial part, we consider three cases: annihilating dark matter, decaying dark matter, and pulsars. In each case, we use GALPROP to inject energy in log-spaced energy bins and compute the expected cosmic-ray and photon signals for each bin. We then fit a linear combination of energy bins, plus backgrounds, to the data. We use a non-parametric fit, with no prior constraints on the spectrum except smoothness and non-negativity. In addition, we consider arbitrary modifications to the energy spectrum of the "ordinary" primary source function, fixing its spatial part, finding this alone to be inadequate to explain the PAMELA or WMAP signals. We explore variations in the fits due to choice of magnetic field, primary electron injection index, spatial profiles, propagation parameters, and fit regularization method. Dark matter annihilation fits well, where our fit finds a mass of ~1 TeV and a boost factor times energy fraction of ~70. While it is possible for dark matter decay and pulsars to fit the data, unconventionally high magnetic fields and radiation densities are required near the Galactic Center to counter the relative shallowness of the assumed spatial profiles. We also fit to linear combinations of these three scenarios, though the fit is much less constrained.

During a survey for stars with disks in the Taurus star-forming region using the Spitzer Space Telescope, we have discovered a pair of young brown dwarfs, FU Tau A and B, in the Barnard 215 dark cloud. They have a projected angular separation of 5.7", corresponding to 800 AU at the distance of Taurus. To assess the nature of these two objects, we have obtained spectra of them and have constructed their spectral energy distributions. Both sources are young (~1 Myr) according to their Halpha emission, gravity-sensitive spectral features, and mid-IR excess emission. The proper motion of FU Tau A provides additional evidence of its membership in Taurus. We measure spectral types of M7.25 and M9.25 for FU Tau A and B, respectively, which correspond to masses of ~0.05 and ~0.015 M\cdot according to the evolutionary models of Chabrier and Baraffe. FU Tau A is significantly overluminous relative to an isochrone passing through FU Tau B and relative to other members of Taurus near its spectral type, which may indicate that it is an unresolved binary. FU Tau A and B are likely to be components of a binary system based on the low probability (~3x10^-4) that Taurus would produce two unrelated brown dwarfs with a projected separation of a </- 6". Barnard 215 contains only one other young star and is in a remote area of Taurus, making FU Tau A and B the first spectroscopically-confirmed brown dwarfs discovered forming in isolation rather than in a stellar cluster or aggregate. Because they were born in isolation and comprise a weakly bound binary, dynamical interactions with stars could not have played a role in their formation, and thus are not essential for the birth of brown dwarfs. ERRATUM: The K-band magnitude for FU Tau B in Table 1 is incorrect and should be 13.33. The bolometric luminosity of FU Tau B in Table 3 and Figure 5 is incorrect because of that mistake and a separate arithmetic error. The correct value of the luminosity is 0.0039 Lsun. FU Tau A and B exhibited different isochronal ages in the original Hertzsprung-Russell diagram in Figure 5, which was unexpected for members of a binary system. This discrepancy is reduced in the corrected version of Figure 5 since both objects are now above the isochrone for 1 Myr. Given the large uncertainties in model isochrones at such young ages, the positions of FU Tau A and B in Figure 5 could be roughly consistent with coevality.

We describe a simple test of the spatial uniformity of an ensemble of discrete events. Given an estimate for the point source luminosity function and an instrumental point spread function (PSF), a robust upper bound on the fractional point source contribution to a diffuse signal can be found. We verify with Monte Carlo tests that the statistic has advantages over the two-point correlation function for this purpose, and derive analytic estimates of the statistic's mean and variance as a function of the point source contribution. As a case study, we apply this statistic to recent gamma-ray data from the Fermi Large Area Telescope (LAT), and demonstrate that at energies above 10 GeV, the contribution of unresolved point sources to the diffuse emission is small in the region relevant for study of the WMAP Haze.

The Fermi Gamma-Ray Space Telescope reveals a diffuse inverse Compton signal in the inner Galaxy with a similar spatial morphology to the microwave haze observed by WMAP, supporting the synchrotron interpretation of the microwave signal. Using spatial templates, we regress out pi0 gammas, as well as IC and bremsstrahlung components associated with known soft-synchrotron counterparts. We find a significant gamma-ray excess towards the Galactic center with a spectrum that is significantly harder than other sky components and is most consistent with IC from a hard population of electrons. The morphology and spectrum are consistent with it being the IC counterpart to the electrons which generate the microwave haze seen at WMAP frequencies. In addition, the implied electron spectrum is hard; electrons accelerated in supernova shocks in the disk which then diffuse a few kpc to the haze region would have a softer spectrum. We describe the full sky Fermi maps used in this analysis and make them available for download.

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

In this report we discuss the impact of polarized foregrounds on a future CMBPol satellite mission. We review our current knowledge of Galactic polarized emission at microwave frequencies, including synchrotron and thermal dust emission. We use existing data and our understanding of the physical behavior of the sources of foreground emission to generate sky templates, and start to assess how well primordial gravitational wave signals can be separated from foreground contaminants for a CMBPol mission. At the estimated foreground minimum of ~100 GHz, the polarized foregrounds are expected to be lower than a primordial polarization signal with tensor-to-scalar ratio r=0.01, in a small patch (~1%) of the sky known to have low Galactic emission. Over 75% of the sky we expect the foreground amplitude to exceed the primordial signal by about a factor of eight at the foreground minimum and on scales of two degrees. Only on the largest scales does the polarized foreground amplitude exceed the primordial signal by a larger factor of about 20. The prospects for detecting an r=0.01 signal including degree-scale measurements appear promising, with 5 sigma_r ~0.003 forecast from multiple methods. A mission that observes a range of scales offers better prospects from the foregrounds perspective than one targeting only the lowest few multipoles. We begin to explore how optimizing the composition of frequency channels in the focal plane can maximize our ability to perform component separation, with a range of typically 40 < nu < 300 GHz preferred for ten channels. Foreground cleaning methods are already in place to tackle a CMBPol mission data set, and further investigation of the optimization and detectability of the primordial signal will be useful for mission design.

We present the case for a dark matter detector with directional sensitivity. This document was developed at the 2009 CYGNUS workshop on directional dark matter detection, and contains contributions from theorists and experimental groups in the field. We describe the need for a dark matter detector with directional sensitivity; each directional dark matter experiment presents their project's status; and we close with a feasibility study for scaling up to a one ton directional detector, which would cost around $150M.

Multiple lines of evidence indicate an anomalous injection of high-energy e+- in the Galactic halo. The recent $e^+$ fraction spectrum from the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA) shows a sharp rise up to 100 GeV. The Fermi Gamma-ray Space Telescope has found a significant hardening of the e+e- cosmic ray spectrum above 100 GeV, with a break, confirmed by HESS at around 1 TeV. The Advanced Thin Ionization Calorimeter (ATIC) has also detected detected a similar excess, falling back to the expected spectrum at 1 TeV and above. Excess microwaves towards the galactic center in the WMAP data are consistent with hard synchrotron radiation from a population of 10-100 GeV e+- (the WMAP ``Haze''). We argue that dark matter annihilations can provide a consistent explanation of all of these data, focusing on dominantly leptonic modes, either directly or through a new light boson. Normalizing the signal to the highest energy evidence (Fermi and HESS), we find that similar cross sections provide good fits to PAMELA and the Haze, and that both the required cross section and annihilation modes are achievable in models with Sommerfeld-enhanced annihilation. These models naturally predict significant production of gamma rays in the galactic center via a variety of mechanisms. Most notably, there is a robust inverse-Compton scattered (ICS) gamma-ray signal arising from the energetic electrons and positrons, detectable at Fermi/GLAST energies, which should provide smoking gun evidence for this production.

Recently, a preliminary spectrum from the Fermi Gamma-ray Space Telescope has been presented for the inner galaxy (-30 < l < 30, -5 < b < 5), as well as the galactic center (-1 < l < 1, -1< b < 1). We consider the implications of these data for dark matter annihilation models, especially models capable of producing the cosmic-ray excesses previously observed by PAMELA and Fermi. These local cosmic-ray excesses, when extrapolated to the inner galaxy, imply inverse Compton scattering (ICS) gamma-ray signals largely consistent with the preliminary Fermi gamma-ray spectrum. For specific halos and models, particularly those with prompt photons, the data have begun to constrain the allowed parameter space. Although significant modeling and background uncertainties remain, we explore how large a signal is permitted by the current data. Based upon this, we make predictions for what signal could be present in other regions of the sky where dark matter signals may be easier to isolate from the astrophysical backgrounds.

Models of dark matter with ~ GeV scale force mediators provide attractive explanations of many high energy anomalies, including PAMELA, ATIC, and the WMAP haze. At the same time, by exploiting the ~ MeV scale excited states that are automatically present in such theories, these models naturally explain the DAMA/LIBRA and INTEGRAL signals through the inelastic dark matter (iDM) and exciting dark matter (XDM) scenarios, respectively. Interestingly, with only weak kinetic mixing to hypercharge to mediate decays, the lifetime of excited states with delta < 2 m_e is longer than the age of the universe. The fractional relic abundance of these excited states depends on the temperature of kinetic decoupling, but can be appreciable. There could easily be other mechanisms for rapid decay, but the consequences of such long-lived states are intriguing. We find that CDMS constrains the fractional relic population of ~100 keV states to be <~ 10^-2, for a 1 TeV WIMP with sigma_n = 10^-40 cm^2. Upcoming searches at CDMS, as well as xenon, silicon, and argon targets, can push this limit significantly lower. We also consider the possibility that the DAMA excitation occurs from a metastable state into the XDM state, which decays via e+e- emission, which allows lighter states to explain the INTEGRAL signal due to the small kinetic energies required. Such models yield dramatic signals from down-scattering, with spectra peaking at high energies, sometimes as high as ~1 MeV, well outside the usual search windows. Such signals would be visible at future Ar and Si experiments, and may be visible at Ge and Xe experiments. We also consider other XDM models involving ~ 500 keV metastable states, and find they can allow lighter WIMPs to explain INTEGRAL as well.

We compute in detail the rate at which energy injected by dark matter annihilation heats and ionizes the photon-baryon plasma at z ~ 1000, and provide accurate fitting functions over the relevant redshift range for a broad array of annihilation channels and DM masses. The resulting perturbations to the ionization history can be constrained by measurements of the CMB temperature and polarization angular power spectra. We show that models which fit recently measured excesses in 10-1000 GeV electron and positron cosmic rays are already close to the 95% confidence limits from WMAP. The recently launched Planck satellite will be capable of ruling out a wide range of DM explanations for these excesses. In models of dark matter with Sommerfeld-enhanced annihilation, where sigma v rises with decreasing WIMP velocity until some saturation point, the WMAP5 constraints imply that the enhancement must be close to saturation in the neighborhood of the Earth.

Recently published results from the PAMELA experiment have shown conclusive evidence for an excess of positrons at high (~ 10 - 100 GeV) energies, confirming earlier indications from HEAT and AMS-01. Such a signal is generally expected from dark matter annihilations. However, the hard positron spectrum and large amplitude are difficult to achieve in most conventional WIMP models. The absence of any associated excess in anti-protons is highly constraining on any model with hadronic annihilation modes. We revisit an earlier proposal, whereby the dark matter annihilates into a new light (<~GeV) boson phi, which is kinematically constrained to go to hard leptonic states, without anti-protons or pi0's. We find this provides a very good fit to the data. The light boson naturally provides a mechanism by which large cross sections can be achieved through the Sommerfeld enhancement, as was recently proposed. Depending on the mass of the WIMP, the rise may continue above 300 GeV, the extent of PAMELA's ability to discriminate electrons and positrons.

Modern cosmology has sharpened questions posed for millennia about the origin of our cosmic habitat. The age-old questions have been transformed into two pressing issues primed for attack in the coming decade: How did the Universe begin? and What physical laws govern the Universe at the highest energies? The clearest window onto these questions is the pattern of polarization in the Cosmic Microwave Background (CMB), which is uniquely sensitive to primordial gravity waves. A detection of the special pattern produced by gravity waves would be not only an unprecedented discovery, but also a direct probe of physics at the earliest observable instants of our Universe. Experiments which map CMB polarization over the coming decade will lead us on our first steps towards answering these age-old questions.

There is a vast menagerie of plausible candidates for the constituents of dark matter, both within and beyond extensions of the Standard Model of particle physics. Each of these candidates may have scattering (and other) cross section properties that are consistent with the dark matter abundance, BBN, and the most scales in the matter power spectrum; but which may have vastly different behavior at sub-galactic "cutoff" scales, below which dark matter density fluctuations are smoothed out. The only way to quantitatively measure the power spectrum behavior at sub-galactic scales at distances beyond the local universe, and indeed over cosmic time, is through probes available in multiply imaged strong gravitational lenses. Gravitational potential perturbations by dark matter substructure encode information in the observed relative magnifications, positions, and time delays in a strong lens. Each of these is sensitive to a different moment of the substructure mass function and to different effective mass ranges of the substructure. The time delay perturbations, in particular, are proving to be largely immune to the degeneracies and systematic uncertainties that have impacted exploitation of strong lenses for such studies. There is great potential for a coordinated theoretical and observational effort to enable a sophisticated exploitation of strong gravitational lenses as direct probes of dark matter properties. This opportunity motivates this white paper, and drives the need for: a) strong support of the theoretical work necessary to understand all astrophysical consequences for different dark matter candidates; and b) tailored observational campaigns, and even a fully dedicated mission, to obtain the requisite data.

The dark cloud Lynds 1622 is one of a few specific sites in the Galaxy where, relative to observed free-free and vibrational dust emission, there is a clear excess of microwave emission. In order to constrain models for this microwave emission, and to better establish the contribution which it might make to ongoing and near-future microwave background polarization experiments, we have used the Green Bank Telescope to search for linear polarization at 9.65 Ghz towards Lynds 1622. We place a 95.4% upper limit of 88 micro-Kelvin (123 micro-Kelvin at 99.7 confidence) on the total linear polarization of this source averaged over a 1'.3 FWHM beam. Relative to the observed level of anomalous emission in Stokes I these limits correspond to fractional linear polarizations of 2.7% and 3.5%.

We propose a comprehensive theory of dark matter that explains the recent proliferation of unexpected observations in high-energy astrophysics. Cosmic ray spectra from ATIC and PAMELA require a WIMP with mass M_chi ~ 500 - 800 GeV that annihilates into leptons at a level well above that expected from a thermal relic. Signals from WMAP and EGRET reinforce this interpretation. Taken together, we argue these facts imply the presence of a GeV-scale new force in the dark sector. The long range allows a Sommerfeld enhancement to boost the annihilation cross section as required, without altering the weak scale annihilation cross section during dark matter freezeout in the early universe. If the dark matter annihilates into the new force carrier, phi, its low mass can force it to decay dominantly into leptons. If the force carrier is a non-Abelian gauge boson, the dark matter is part of a multiplet of states, and splittings between these states are naturally generated with size alpha m_phi ~ MeV, leading to the eXciting dark matter (XDM) scenario previously proposed to explain the positron annihilation in the galactic center observed by the INTEGRAL satellite. Somewhat smaller splittings would also be expected, providing a natural source for the parameters of the inelastic dark matter (iDM) explanation for the DAMA annual modulation signal. Since the Sommerfeld enhancement is most significant at low velocities, early dark matter halos at redshift ~10 potentially produce observable effects on the ionization history of the universe, and substructure is more detectable than with a conventional WIMP. Moreover, the low velocity dispersion of dwarf galaxies and Milky Way subhalos can greatly increase the substructure annihilation signal.

We characterize spinning dust emission in the warm ionized medium by comparing templates of Galactic dust and Halpha with the 5-year maps from the Wilkinson Microwave Anisotropy Probe. The Halpha-correlated microwave emission deviates from the thermal bremsstrahlung (free-free) spectrum expected for ionized gas, exhibiting an additional broad bump peaked at ~40 GHz which provides ~20% of the peak intensity. We confirm that the bump is consistent with a modified Draine & Lazarian (1998) spinning dust model, though the peak frequency of the emission is somewhat lower than the 50 GHz previously claimed. This frequency shift results from systematic errors in the large-scale modes of the 3-year WMAP data which have been corrected in the 5-year data release. We show that the bump is not the result of errors in the Halpha template by analyzing regions of high free-free intensity, where the WMAP K-band map may be used as the free-free template. We rule out a pure free-free spectrum for the Halpha-correlated emission at high confidence: ~27sigma for the nearly full-sky fit, even after marginalizing over the CMB cross-correlation bias. We also extend the previous analysis by searching the parameter space of the Draine & Lazarian model but letting the amplitude float. The best fit for reasonable values of the characteristic electric dipole moment and density requires an amplitude factor of ~0.3. This suggests that small PAHs in the warm ionized medium are depleted by a factor of ~3.

We consider the nuclear scattering cross section for the eXciting Dark Matter (XDM) model. In XDM, the Weakly Interacting Massive Particles (WIMPs) couple to the Standard Model only via an intermediate light scalar which mixes with the Higgs: this leads to a suppression in the nuclear scattering cross section relative to models in which the WIMPs couple to the Higgs directly. We estimate this suppression factor to be of order 10^(-5). The elastic nuclear scattering cross section for XDM can also be computed directly: we perform this computation for XDM coupled to the Higgs sector of the Standard Model and find a spin-independent cross section in the order of 4 x 10^(-13) pb in the decoupling limit, which is not within the range of any near-term direct detection experiments. However, if the XDM dark sector is instead coupled to a two-Higgs-doublet model, the spin-independent nuclear scattering cross section can be enhanced by up to four orders of magnitude for large tan(beta), which should be observable in the upcoming SuperCDMS and ton-scale xenon experiments.

Motivated by the eXciting Dark Matter (XDM) model of Finkbeiner & Weiner, hypothesized to explain the 511 keV signal in the center of the Milky Way, we consider the CMB and 21-cm signatures of models of dark matter with collisional long-lived excited states. We compute the relic excitation fraction from the early universe for a variety of assumptions about the collisional de-excitation cross-section and thermal decoupling. The relic excitation fraction can be as high as 1% for natural regions of parameter space, but could be orders of magnitude smaller. Since the lifetime of the excited state is naturally greater than 10^13s, we discuss the signatures of such relic excitation on cosmic microwave background (CMB) and high-z 21-cm observations. Such models have potentially richer astrophysical signals than the traditional WIMP annihilations and decays, and may have observable consequences for future generations of experiments.

Using effective temperature and metallicity derived from SDSS spectra for ~60,000 F and G type main sequence stars (0.2<g-r<0.6), we develop polynomial models for estimating these parameters from the SDSS u-g and g-r colors. We apply this method to SDSS photometric data for about 2 million F/G stars and measure the unbiased metallicity distribution for a complete volume-limited sample of stars at distances between 500 pc and 8 kpc. The metallicity distribution can be exquisitely modeled using two components with a spatially varying number ratio, that correspond to disk and halo. The two components also possess the kinematics expected for disk and halo stars. The metallicity of the halo component is spatially invariant, while the median disk metallicity smoothly decreases with distance from the Galactic plane from -0.6 at 500 pc to -0.8 beyond several kpc. The absence of a correlation between metallicity and kinematics for disk stars is in a conflict with the traditional decomposition in terms of thin and thick disks. We detect coherent substructures in the kinematics--metallicity space, such as the Monoceros stream, which rotates faster than the LSR, and has a median metallicity of [Fe/H]=-0.96, with an rms scatter of only ~0.15 dex. We extrapolate our results to the performance expected from the Large Synoptic Survey Telescope (LSST) and estimate that the LSST will obtain metallicity measurements accurate to 0.2 dex or better, with proper motion measurements accurate to ~0.2 mas/yr, for about 200 million F/G dwarf stars within a distance limit of ~100 kpc (g<23.5). [abridged]

Abridged: We estimate the distances to ~48 million stars detected by the Sloan Digital Sky Survey and map their 3D number density distribution in 100 < D < 20 kpc range over 6,500 deg^2 of sky. The data show strong evidence for a Galaxy consisting of an oblate halo, a disk component, and a number of localized overdensities with exponential disk parameters (bias-corrected for an assumed 35% binary fraction) H_1 = 300 pc, L_1 = 2600 pc, H_2 = 900 pc, L_2 = 3600 pc, and local density normalization of 12%. We find the halo to be oblate, with best-fit axis ratio c/a = 0.64, r^{-2.8} profile, and the local halo-to-thin disk normalization of 0.5%. We estimate the errors of derived model parameters to be no larger than ~20% (disk scales) and ~10% (thick disk normalization). While generally consistent with the above model, the density distribution shows a number of statistically significant localized deviations. We detect two overdensities in the thick disk region at (R, Z) ~ (6.5, 1.5)kpc and (R, Z) ~ (9.5, 0.8) kpc, and a remarkable density enhancement in the halo covering >1000deg^2 of sky towards the constellation of Virgo, at distances of ~6-20 kpc. Compared to a region symmetric with respect to the l=0 line, the Virgo overdensity is responsible for a factor of 2 number density excess and may be a nearby tidal stream or a low-surface brightness dwarf galaxy merging with the Milky Way. After removal of the resolved overdensities, the remaining data are consistent with a smooth density distribution; we detect no evidence of further unresolved clumpy substructure at scales ranging from ~50pc in the disk, to ~1 - 2 kpc in the halo.

The detection of large quantities of dust in z ~ 6 quasars by infrared and radio surveys presents puzzles for the formation and evolution of dust in these early systems. Previously (Li et al. 2007), we showed that luminous quasars at z > 6 can form through hierarchical mergers of gas-rich galaxies. Here, we calculate the dust properties of simulated quasars and their progenitors using a three-dimensional Monte Carlo radiative transfer code, ART^2 -- All-wavelength Radiative Transfer with Adaptive Refinement Tree. ART^2 incorporates a radiative equilibrium algorithm for dust emission, an adaptive grid for inhomogeneous density, a multiphase model for the ISM, and a supernova-origin dust model. We reproduce the SED and dust properties of SDSS J1148+5251, and find that the infrared emission are closely associated with the formation and evolution of the quasar host. The system evolves from a cold to a warm ULIRG owing to heating and feedback from stars and AGN. Furthermore, the AGN has significant implications for the interpretation of observation of the hosts. Our results suggest that vigorous star formation in merging progenitors is necessary to reproduce the observed dust properties of z~6 quasars, supporting a merger-driven origin for luminous quasars at high redshifts and the starburst-to-quasar evolutionary hypothesis. (Abridged)

CMB experiments commonly use maps of Halpha intensity as a spatial template for Galactic free-free emission, assuming a power law I_nu \propto nu^-0.15 for the spectrum. Any departure from the assumed free-free spectrum could have a detrimental effect on determination of the primary CMB anisotropy. We show that the Halpha-correlated emission spectrum in the diffuse warm ionized medium (WIM) is not the expected free-free spectrum at WMAP frequencies. Instead, there is a broad bump in the spectrum at ~50 GHz which is consistent with emission from spinning dust grains. Spectra from both the full sky and smaller regions of interest are well fit by a superposition of a free-free and WIM Draine & Lazarian (1998) spinning dust model, shifted in frequency. The spinning dust emission is ~5 times weaker than the free-free component at 50 GHz, with the null hypothesis that the Halpha-correlated spectrum is pure free-free, ruled out at >8 sigma in all regions and >100 sigma for the full sky fit.

We study the spectral and morphological characteristics of the diffuse Galactic emission in the WMAP temperature data using a template-based multi-linear regression, and obtain the following results. 1. We confirm previous observations of a bump in the dust-correlated spectrum, consistent with the Draine & Lazarian spinning dust model. 2. We also confirm the "haze" signal in the inner Galaxy, and argue that it does not follow a free-free spectrum as first thought, but instead is synchrotron emission from a hard electron cosmic-ray population. 3. In a departure from previous work, we allow the spectrum of Halpha-correlated emission (which is used to trace the free-free component) to float in the fit, and find that it does not follow the expected free-free spectrum. Instead there is a bump near 50 GHz, modifying the spectrum at the 20% level, which we speculate is caused by spinning dust in the warm ionized medium. 4. The derived cross-correlation spectra are not sensitive to the map zero points, but are sensitive to the choice of CMB estimator. In cases where the CMB estimator is derived by minimizing variance of a linear combination of the WMAP bands, we show that a bias proportional to the cross-correlation of each template and the true CMB is always present. This bias can be larger than any of the foreground signals in some bands. 5. Lastly, we consider the frequency coverage and sensitivity of the Planck mission, and suggest linear combination coefficients for the CMB template that will reduce both the statistical and systematic uncertainty in the synchrotron and haze spectra by more than an order of magnitude.

We present an algorithm to photometrically calibrate wide field optical imaging surveys, that simultaneously solves for the calibration parameters and relative stellar fluxes using overlapping observations. The algorithm decouples the problem of "relative" calibrations, from that of "absolute" calibrations; the absolute calibration is reduced to determining a few numbers for the entire survey. We pay special attention to the spatial structure of the calibration errors, allowing one to isolate particular error modes in downstream analyses. Applying this to the Sloan Digital Sky Survey imaging data, we achieve ~1% relative calibration errors across 8500 sq.deg. in griz; the errors are ~2% for the u band. These errors are dominated by unmodelled atmospheric variations at Apache Point Observatory. These calibrations, dubbed "ubercalibration", are now public with SDSS Data Release 6, and will be a part of subsequent SDSS data releases.

Observations by the WMAP experiment have identified an excess of microwave emission from the center of the Milky Way. It has previously been shown that this "WMAP Haze" could be synchrotron emission from relativistic electrons and positrons produced in the annihilations of dark matter particles. In particular, the intensity, spectrum and angular distribution of the WMAP Haze is consistent with an electroweak scale dark matter particle (such as a supersymmetric neutralino or Kaluza-Klein dark matter in models with universal extra dimensions) annihilating with a cross section on the order of sigma v~3x10^-26 cm^3/s and distributed with a cusped halo profile. No further exotic astrophysical or annihilation boost factors are required. If dark matter annihilations are in fact responsible for the observed Haze, then other annihilation products will also be produced, including gamma rays. In this article, we study the prospects for the GLAST satellite to detect gamma rays from dark matter annihilations in the Galactic Center region in this scenario. We find that by studying only the inner 0.1 degrees around the Galactic Center, GLAST will be able to detect dark matter annihilating to heavy quarks or gauge bosons over astrophysical backgrounds with 5sigma (3sigma) significance if they are lighter than approximately 320-500 GeV (500-750 GeV). If the angular window is broadened to study the dark matter halo profile's angular extension (while simultaneously reducing the astrophysical backgrounds), WIMPs as heavy as several TeV can be identified by GLAST with high significance. Only if the dark matter particles annihilate mostly to electrons or muons will GLAST be unable to identify the gamma ray spectrum associated with the WMAP Haze.

The Sloan Digital Sky Survey (SDSS) and Two Micron All Sky Survey (2MASS) are rich resources for studying stellar astrophysics and the structure and formation history of the Galaxy. As new surveys and instruments adopt similar filter sets, it is increasingly important to understand the properties of the ugrizJHKs stellar locus, both to inform studies of `normal' main sequence stars as well as for robust searches for point sources with unusual colors. Using a sample of ~600,000 point sources detected by SDSS and 2MASS, we tabulate the position and width of the ugrizJHKs stellar locus as a function of g-i color, and provide accurate polynomial fits. We map the Morgan-Keenan spectral type sequence to the median stellar locus by using synthetic photometry of spectral standards and by analyzing 3000 SDSS stellar spectra with a custom spectral typing pipeline. We develop an algorithm to calculate a point source's minimum separation from the stellar locus in a seven-dimensional color space, and use it to robustly identify objects with unusual colors, as well as spurious SDSS/2MASS matches. Analysis of a final catalog of 2117 color outliers identifies 370 white-dwarf/M dwarf (WDMD) pairs, 93 QSOs, and 90 M giant/carbon star candidates, and demonstrates that WDMD pairs and QSOs can be distinguished on the basis of their J-Ks and r-z colors. We also identify a group of objects with correlated offsets in the u-g vs. g-r and g-r vs. r-i color-color spaces, but subsequent follow-up is required to reveal the nature of these objects. Future applications of this algorithm to a matched SDSS-UKIDSS catalog may well identify additional classes of objects with unusual colors by probing new areas of color-magnitude space.

The WMAP experiment has revealed an excess of microwave emission from the region around the center of our Galaxy. It has been suggested that this signal, known as the ``WMAP Haze'', could be synchrotron emission from relativistic electrons and positrons generated in dark matter annihilations. In this letter, we revisit this possibility. We find that the angular distribution of the WMAP Haze matches the prediction for dark matter annihilations with a cusped density profile, $\rho(r) \propto r^{-1.2}$ in the inner kiloparsecs. Comparing the intensity in different WMAP frequency bands, we find that a wide range of possible WIMP annihilation modes are consistent with the spectrum of the haze for a WIMP with a mass in the 100 GeV to multi-TeV range. Most interestingly, we find that to generate the observed intensity of the haze, the dark matter annihilation cross section is required to be approximately equal to the value needed for a thermal relic, $\sigma v \sim 3 \times 10^{-26}$ cm$^3$/s. No boost factors are required. If dark matter annihilations are in fact responsible for the WMAP Haze, and the slope of the halo profile continues into the inner Galaxy, GLAST is expected to detect gamma rays from the dark matter annihilations in the Galactic Center if the WIMP mass is less than several hundred GeV.

We propose a WIMP candidate with an ``excited state'' 1-2 MeV above the ground state, which may be collisionally excited and de-excites by e+e- pair emission. By converting its kinetic energy into pairs, such a particle could produce a substantial fraction of the 511 keV line observed by INTEGRAL/SPI in the inner Milky Way. Only a small fraction of the WIMPs have sufficient energy to excite, and that fraction drops sharply with galactocentric radius, naturally yielding a radial cutoff, as observed. Even if the scattering probability in the inner kpc is << 1% per Hubble time, enough power is available to produce the ~3x10^42 pairs per second observed in the Galactic bulge. We specify the parameters of a pseudo-Dirac fermion designed to explain the positron signal, and find that it annihilates chiefly to e+e- and freezes out with the correct relic density. We discuss possible observational consequences of this model.

We quantify the variability of faint unresolved optical sources using a catalog based on multiple SDSS imaging observations. The catalog covers SDSS Stripe 82, and contains 58 million photometric observations in the SDSS ugriz system for 1.4 million unresolved sources. In each photometric bandpass we compute various low-order lightcurve statistics and use them to select and study variable sources. We find that 2% of unresolved optical sources brighter than g=20.5 appear variable at the 0.05 mag level (rms) simultaneously in the g and r bands. The majority (2/3) of these variable sources are low-redshift (<2) quasars, although they represent only 2% of all sources in the adopted flux-limited sample. We find that at least 90% of quasars are variable at the 0.03 mag level (rms) and confirm that variability is as good a method for finding low-redshift quasars as is the UV excess color selection (at high Galactic latitudes). We analyze the distribution of lightcurve skewness for quasars and find that is centered on zero. We find that about 1/4 of the variable stars are RR Lyrae stars, and that only 0.5% of stars from the main stellar locus are variable at the 0.05 mag level. The distribution of lightcurve skewness in the g-r vs. u-g color-color diagram on the main stellar locus is found to be bimodal (with one mode consistent with Algol-like behavior). Using over six hundred RR Lyrae stars, we demonstrate rich halo substructure out to distances of 100 kpc. We extrapolate these results to expected performance by the Large Synoptic Survey Telescope and estimate that it will obtain well-sampled 2% accurate, multi-color lightcurves for ~2 million low-redshift quasars, and will discover at least 50 million variable stars.

We compare Stetson's photometric standards with measurements listed in a standard star catalog constructed using repeated SDSS imaging observations. The SDSS catalog includes over 700,000 candidate standard stars from the equatorial stripe 82 (|Dec|<1.266 deg) in the RA range 20h 34' to 4h 00', and with the $r$ band magnitudes in the range 14--21. The distributions of measurements for individual sources demonstrate that the SDSS photometric pipeline correctly estimates random photometric errors, which are below 0.01 mag for stars brighter than (19.5, 20.5, 20.5, 20, 18.5) in ugriz, respectively (about twice as good as for individual SDSS runs). We derive mean photometric transformations between the SDSS gri and the BVRI system using 1165 Stetson stars found in the equatorial stripe 82, and then study the spatial variation of the difference in zeropoints between the two catalogs. Using third order polynomials to describe the color terms, we find that photometric measurements for main-sequence stars can be transformed between the two systems with systematic errors smaller than a few millimagnitudes. The spatial variation of photometric zeropoints in the two catalogs typically does not exceed 0.01 magnitude. Consequently, the SDSS Standard Star Catalog for Stripe 82 can be used to calibrate new data in both the SDSS ugriz and the BVRI systems with a similar accuracy.

We measure the large-scale real-space power spectrum P(k) using luminous red galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS) and use this measurement to sharpen constraints on cosmological parameters from the Wilkinson Microwave Anisotropy Probe (WMAP). We employ a matrix-based power spectrum estimation method using Pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 20 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.01h/Mpc < k < 0.2h/Mpc. Results from the LRG and main galaxy samples are consistent, with the former providing higher signal-to-noise. Our results are robust to omitting angular and radial density fluctuations and are consistent between different parts of the sky. They provide a striking confirmation of the predicted large-scale LCDM power spectrum. Combining only SDSS LRG and WMAP data places robust constraints on many cosmological parameters that complement prior analyses of multiple data sets. The LRGs provide independent cross-checks on Om and the baryon fraction in good agreement with WMAP. Within the context of flat LCDM models, our LRG measurements complement WMAP by sharpening the constraints on the matter density, the neutrino density and the tensor amplitude by about a factor of two, giving Omega_m=0.24+-0.02 (1 sigma), sum m_nu < 0.9 eV (95%) and r<0.3 (95%). Baryon oscillations are clearly detected and provide a robust measurement of the comoving distance to the median survey redshift z=0.35 independent of curvature and dark energy properties. Within the LCDM framework, our power spectrum measurement improves the evidence for spatial flatness, sharpening the curvature constraint Omega_tot=1.05+-0.05 from WMAP alone to Omega_tot=1.003+-0.010. Assuming Omega_tot=1, the equation of state parameter is constrained to w=-0.94+-0.09, indicating the potential for more ambitious future LRG measurements to provide precision tests of the nature of dark energy. All these constraints are essentially independent of scales k>0.1h/Mpc and associated nonlinear complications, yet agree well with more aggressive published analyses where nonlinear modeling is crucial.

The photometric calibration of the Sloan Digital Sky Survey (SDSS) is a multi-step process which involves data from three different telescopes: the 1.0-m telescope at the US Naval Observatory (USNO), Flagstaff Station, Arizona (which was used to establish the SDSS standard star network); the SDSS 0.5-m Photometric Telescope (PT) at the Apache Point Observatory (APO), New Mexico (which calculates nightly extinctions and calibrates secondary patch transfer fields); and the SDSS 2.5-m telescope at APO (which obtains the imaging data for the SDSS proper). In this paper, we describe the Monitor Telescope Pipeline, MTPIPE, the software pipeline used in processing the data from the single-CCD telescopes used in the photometric calibration of the SDSS (i.e., the USNO 1.0-m and the PT). We also describe transformation equations that convert photometry on the USNO-1.0m u'g'r'i'z' system to photometry the SDSS 2.5m ugriz system and the results of various validation tests of the MTPIPE software. Further, we discuss the semi-automated PT factory, which runs MTPIPE in the day-to-day standard SDSS operations at Fermilab. Finally, we discuss the use of MTPIPE in current SDSS-related projects, including the Southern u'g'r'i'z' Standard Star project, the u'g'r'i'z' Open Star Clusters project, and the SDSS extension (SDSS-II).

We compare the most successful and widely used map of Galactic dust extinction, provided by Schlegel, Finkbeiner & Davis (1998; hereafter SFD), to the galaxy number counts in the Sloan Digital Sky Survey (SDSS) photometric/spectroscopic DR4 sample. We divide the SDSS survey area into 69 disjoint subregions according to the dust extinction provided by SFD and compare the surface number density of galaxies in each subregion. As expected, the galaxy surface number density decreases with increasing extinction but only for SFD extinction values above about 0.1 to 0.2 magnitudes (depending on the band). At lower values of the SFD extinction, we find that the sky surface density of galaxies increases with increasing extinction, precisely the opposite of the effect expected from Galactic dust. We suggest that the far infrared (FIR) brightness of the sky in regions of true low dust extinction is significantly ``contaminated'' by the FIR emission from background galaxies. We show that such an explanation is both qualitatively and quantitatively consistent with the available data. Based on this interpretation we conclude that systematic errors in the SFD extinction map due to extragalactic FIR emission are quite small, of order hundredths of a magnitude, but nevertheless statistically detectable. (Abridged)

We discuss the panchromatic properties of 99,088 galaxies selected from the SDSS Data Release 1 spectroscopic sample (a flux-limited sample for 1360 deg^2). These galaxies are positionally matched to sources detected by ROSAT, GALEX, 2MASS, IRAS, GB6, FIRST, NVSS and WENSS. We find strong correlations between the detection fraction at other wavelengths and optical properties such as flux, colors, and emission-line strengths. Using GALEX, SDSS, and 2MASS data, we construct the UV-IR broad-band spectral energy distributions for various types of galaxies, and find that they form a nearly one-parameter family. For example, based on SDSS u- and r-band data, supplemented with redshift, the K-band 2MASS magnitudes can be "predicted" with an rms scatter of only 0.2 mag. When a dust content estimate determined from SDSS data by Kauffmann et al. (2003) is also utilized, this scatter decreases to 0.1 mag. We demonstrate that this dust content is indeed higher for galaxies detected by IRAS and that it can be used to "predict" measured IRAS 60 micron flux density within a factor of two using only SDSS data. We also show that the position of a galaxy in the emission-line-based Baldwin-Phillips-Terlevich diagram is correlated with the optical light concentration index and u-r color determined from the SDSS broad-band imaging data, and discuss changes in the morphology of this diagram induced by requiring detections at other wavelengths. We study the IR-radio correlation and find evidence that its slope may be different for AGN and star-forming galaxies and related to the H_alpha/H_beta line strength ratio.

We present the 3D real space clustering power spectrum of a sample of \~600,000 luminous red galaxies (LRGs) measured by the Sloan Digital Sky Survey (SDSS), using photometric redshifts. This sample of galaxies ranges from redshift z=0.2 to 0.6 over 3,528 deg^2 of the sky, probing a volume of 1.5 (Gpc/h)^3, making it the largest volume ever used for galaxy clustering measurements. We measure the angular clustering power spectrum in eight redshift slices and combine these into a high precision 3D real space power spectrum from k=0.005 (h/Mpc) to k=1 (h/Mpc). We detect power on gigaparsec scales, beyond the turnover in the matter power spectrum, on scales significantly larger than those accessible to current spectroscopic redshift surveys. We also find evidence for baryonic oscillations, both in the power spectrum, as well as in fits to the baryon density, at a 2.5 sigma confidence level. The statistical power of these data to constrain cosmology is ~1.7 times better than previous clustering analyses. Varying the matter density and baryon fraction, we find \Omega_M = 0.30 \pm 0.03, and \Omega_b/\Omega_M = 0.18 \pm 0.04, The detection of baryonic oscillations also allows us to measure the comoving distance to z=0.5; we find a best fit distance of 1.73 \pm 0.12 Gpc, corresponding to a 6.5% error on the distance. These results demonstrate the ability to make precise clustering measurements with photometric surveys (abridged).

We study the mean environment of galaxies in the DEEP2 Galaxy Redshift Survey as a function of rest-frame color, luminosity, and [OII] equivalent width. The local galaxy overdensity for >14,000 galaxies at 0.75 < z < 1.35 is estimated using the projected 3rd-nearest-neighbor surface density. Of the galaxy properties studied, mean environment is found to depend most strongly on galaxy color; all major features of the correlation between mean overdensity and rest-frame color observed in the local universe were already in place at z ~ 1. In contrast to local results, we find a substantial slope in the mean dependence of environment on luminosity for blue, star-forming galaxies at z ~ 1, with brighter blue galaxies being found on average in regions of greater overdensity. We discuss the roles of galaxy clusters and groups in establishing the observed correlations between environment and galaxy properties at high redshift, and we also explore the evidence for a ``downsizing of quenching'' from z ~ 1 to z ~ 0. Our results add weight to existing evidence that the mechanism(s) that result in star-formation quenching are efficient in group environments as well as clusters. This work is the first of its kind at high redshift and represents the first in a series of papers addressing the role of environment in galaxy formation at 0 < z < 1.

The evolution of the B-band galaxy luminosity function is measured using a sample of more than 11,000 galaxies with spectroscopic redshifts from the DEEP2 Redshift Survey. The rest-frame M_B versus U-B color-magnitude diagram of DEEP2 galaxies shows that the color-magnitude bi-modality seen in galaxies locally is still present at redshifts z > 1. Dividing the sample at the trough of this color bimodality into predominantly red and blue galaxies, we find that the luminosity function of each galaxy color type evolves differently. Blue counts tend to shift to brighter magnitudes at constant number density, while the red counts remain largely constant at a fixed absolute magnitude. Using Schechter functions with fixed faint-end slopes we find that M*_B for blue galaxies brightens by ~ 1.3 magnitudes per unit redshift, with no significant evolution in number density. For red galaxies M*_B brightens somewhat less with redshift, while the formal value of phi* declines. When the population of blue galaxies is subdivided into two halves using the rest-frame color as the criterion, the measured evolution of both blue subpopulations is very similar.

The DEEP2 and COMBO-17 surveys are used to study the evolution of the luminosity function of red and blue galaxies to $z \sim 1$. Schechter function fits show that, since $z = 1$, $M^*_B$ dims by $\sim$ 1.3 mag per unit redshift for both color classes, $\phi^*$ of blue galaxies shows little change, while $\phi^*$ for red galaxies has formally nearly quadrupled. At face value, the number density of blue galaxies has remained roughly constant since $ z = 1$, whereas that of red galaxies has been rising. Luminosity densities support both conclusions, but we note that most red-galaxy evolution occurs between our data and local surveys and in our highest redshift bin, where the data are weakest. We discuss the implications of having most red galaxies emerge after $z = 1$ from precursors among the blue population, taking into account the properties of local and distant E/S0s. We suggest a ``mixed'' scenario in which some blue galaxies have their star-formation quenched in gas-rich mergers, migrate to the red sequence with a variety of masses, and merge further on the red sequence in one or more purely stellar mergers. E/S0s of a given mass today will have formed via different routes, in a manner that may help to explain the fundamental plane and other local scaling laws.

Dark matter (DM) annihilation during hydrogen recombination (z ~ 1000) will alter the recombination history of the Universe, and affect the observed CMB temperature and polarization fluctuations. Unlike other astrophysical probes of DM, this is free of the significant uncertainties in modelling galactic physics, and provides a method to detect and constrain the cosmological abundances of these particles. We parametrize the effect of DM annihilation as an injection of ionizing energy at a rate e_{dm}, and argue that this simple "on the spot'' modification is a good approximation to the complicated interaction of the annihilation products with the photon-electron plasma. Generic models of DM do not change the redshift of recombination, but change the residual ionization after recombination. This broadens the surface of last scattering, suppressing the temperature fluctuations and enhancing the polarization fluctuations. We use the temperature and polarization angular power spectra to measure these deviations from the standard recombination history, and therefore, indirectly probe DM annihilation. (abridged)

Here we present the New York University Value-Added Galaxy Catalog (NYU-VAGC), a catalog of local galaxies (mostly below a redshift of about 0.3) based on a set of publicly-released surveys (including the 2dFGRS, 2MASS, PSCz, FIRST, and RC3) matched to the Sloan Digital Sky Survey (SDSS) Data Release 2. Excluding areas masked by bright stars, the photometric sample covers 3514 square degrees and the spectroscopic sample covers 2627 square degrees (with about 85% completeness). Earlier, proprietary versions of this catalog have formed the basis of many SDSS investigations of the power spectrum, correlation function, and luminosity function of galaxies. We calculate and compile derived quantities (for example, K-corrections and structural parameters for galaxies). The SDSS catalog presented here is photometrically recalibrated, reducing systematic calibration errors across the sky from about 2% to about 1%. We include an explicit description of the geometry of the catalog, including all imaging and targeting information as a function of sky position. Finally, we have performed eyeball quality checks on a large number of objects in the catalog in order to flag deblending and other errors. This catalog is complementary to the SDSS Archive Servers, in that NYU-VAGC's calibration, geometrical description, and conveniently small size are specifically designed for studying galaxy properties and large-scale structure statistics using the SDSS spectroscopic catalog.

We use the first 25% of the DEEP2 Galaxy Redshift Survey spectroscopic data to identify groups and clusters of galaxies in redshift space. The data set contains 8370 galaxies with confirmed redshifts in the range 0.7 <= z <= 1.4, over one square degree on the sky. Groups are identified using an algorithm (the Voronoi-Delaunay Method) that has been shown to accurately reproduce the statistics of groups in simulated DEEP2-like samples. We optimize this algorithm for the DEEP2 survey by applying it to realistic mock galaxy catalogs and assessing the results using a stringent set of criteria for measuring group-finding success, which we develop and describe in detail here. We find in particular that the group-finder can successfully identify ~78% of real groups and that ~79% of the galaxies that are true members of groups can be identified as such. Conversely, we estimate that ~55% of the groups we find can be definitively identified with real groups and that ~46% of the galaxies we place into groups are interloper field galaxies. Most importantly, we find that it is possible to measure the distribution of groups in redshift and velocity dispersion, n(sigma,z), to an accuracy limited by cosmic variance, for dispersions greater than 350km/s. We anticipate that such measurements will allow strong constraints to be placed on the equation of state of the dark energy in the future. Finally, we present the first DEEP2 group catalog, which assigns 32% of the galaxies to 899 distinct groups with two or more members, 153 of which have velocity dispersions above 350km/s. We provide locations, redshifts and properties for this high-dispersion subsample. This catalog represents the largest sample to date of spectroscopically detected groups at z~1.

Excess microwave emission observed in the inner Galaxy (inner ~1 kpc) is consistent with synchrotron emission from highly relativistic electron-positron pairs produced by dark matter particle annihilation. More conventional sources for this emission, such as free-free (thermal bremsstrahlung), thermal dust, spinning dust, and the softer Galactic synchrotron traced by low-frequency surveys, have been ruled out. The total power observed in the range 23 < nu < 61 GHz is between 10^{36} and 5x10^{36} erg/s, depending on the method of extrapolation to the Galactic center, where bright foreground emission obscures the signal. The inferred electron energy distribution is diffusion hardened, and is in qualitative agreement with the energy distribution required to explain the gamma ray excess in the inner Galaxy at 1-30 GeV as inverse-Compton scattered starlight. We investigate the possibility that this population of electrons is produced by dark matter annihilation of 100 GeV particles, with cross section <sigma v>=2x10^{-26} cm^3/s, and an 1/r dark matter mass profile truncated in the inner Galaxy, and find this scenario to be consistent with current data.

We summarize the Sloan Digital Sky Survey data acquisition and processing steps, and describe runQA, a pipeline designed for automated data quality assessment. In particular, we show how the position of the stellar locus in color-color diagrams can be used to estimate the accuracy of photometric zeropoint calibration to better than 0.01 mag in 0.03 deg2 patches. Using this method, we estimate that typical photometric zeropoint calibration errors for SDSS imaging data are not larger than ~0.01 mag in the g, r, and i bands, 0.02 mag in the z band, and 0.03 mag in the u band (root-mean-scatter for zeropoint offsets).

The Sloan Digital Sky Survey (SDSS) mosaic camera and telescope have obtained five-band optical-wavelength imaging near the Galactic plane outside of the nominal survey boundaries. These additional data were obtained during commissioning and subsequent testing of the SDSS observing system, and they provide unique wide-area imaging data in regions of high obscuration and star formation, including numerous young stellar objects, Herbig-Haro objects and young star clusters. Because these data are outside the Survey regions in the Galactic caps, they are not part of the standard SDSS data releases. This paper presents imaging data for 832 square degrees of sky (including repeats), in the star-forming regions of Orion, Taurus, and Cygnus. About 470 square degrees are now released to the public, with the remainder to follow at the time of SDSS Data Release 4. The public data in Orion include the star-forming region NGC 2068/NGC 2071/HH24 and a large part of Barnard's loop.

Observations of the CMB can be contaminated by diffuse foreground emission from sources such as Galactic dust and synchrotron radiation. In these cases, the morphology of the contaminating source is known from observations at different frequencies, but not its amplitude at the frequency of interest for the CMB. We develop a technique for accounting for the effects of such emission in this case, and for simultaneously estimating the foreground amplitude in the CMB observations. We apply the technique to CMB data from the MAXIMA-1 experiment, using maps of Galactic dust emission from combinations of IRAS and DIRBE observations, as well as compilations of Galactic synchrotron emission observations. The spectrum of the dust emission over the 150--450 GHz observed by MAXIMA is consistent with preferred models but the effect on CMB power spectrum observations is negligible.

The Spitzer Space Telescope First Look Survey Infrared Array Camera (IRAC) near and mid-infrared imaging data partially overlaps the Sloan Digital Sky Survey (SDSS), with 313 visually selected (r<17.6 mag) SDSS Main Sample galaxies in the overlap region. The 3.5 and 7.8 um properties of the galaxies are investigated in the context of their visual properties, where the IRAC [3.5] magnitude primarily measures starlight, and the [7.8] magnitude primarily measures PAH emission from the interstellar medium. As expected, we find a strong inverse correlation between [3.5]-[7.8] and visual color; galaxies red in visual colors (`red galaxies') tend to show very little dust and molecular emission (low `PAH-to-star' ratios), and galaxies blue in visual colors (`blue galaxies,' ie, star-forming galaxies) tend to show large PAH-to-star ratios. Red galaxies with high PAH-to-star ratios tend to be edge-on disks reddened by dust lanes. Simple, visually inferred attenuation corrections bring the visual colors of these galaxies in line with those of face-on disks; ie, PAH emission is closely related to attenuation-corrected, optically inferred star-formation rates. Blue galaxies with anomalously low PAH-to-star ratios are all low-luminosity star-forming galaxies. There is some weak evidence in this sample that the deficiency in PAH emission for these low-luminosity galaxies may be related to emission-line metallicity.

We present Sloan Digital Sky Survey and United States Naval Observatory observations of the V1647 Ori protostar and surrounding field near NGC 2068. V1647 Ori, the likely driving source for HH 23, brightened significantly in November 2003. Analysis of SDSS imaging acquired in November 1998 and February 2002 during the quiescent state, recent USNO photometry, and published 2MASS and Gemini data shows that the color changes associated with brightening suggest an EXor outburst rather than a simple dust clearing event.

We observe significant dust-correlated emission outside of H II regions in the Green Bank Galactic Plane Survey (-4 < b < 4 degrees) at 8.35 and 14.35 GHz. The rising spectral slope rules out synchrotron and free-free emission as majority constituents at 14 GHz, and the amplitude is at least 500 times higher than expected thermal dust emission. When combined with the Rhodes (2.326 GHz), and WMAP (23-94 GHz) data it is possible to fit dust-correlated emission at 2.3-94 GHz with only soft synchrotron, free-free, thermal dust, and an additional dust-correlated component similar to Draine & Lazarian spinning dust. The rising component generally dominates free-free and synchrotron for \nu >~ 14 GHz and is overwhelmed by thermal dust at \nu > 60 GHz. The current data fulfill most of the criteria laid out by Finkbeiner et al. (2002) for detection of spinning dust.

We measure the two-point correlation function xi(r) using a sample of 2219 galaxies in an area of 0.32 degrees^2 at z=0.7-1.35 from the first season of the DEEP2 Galaxy Redshift Survey. We find that xi(r) can be approximated by a power-law, xi(r)=(r/r_0)^-gamma, on scales 0.1-20 Mpc/h. In a sample with an effective redshift of z_eff=0.82, for a Lcdm cosmology we find r_0=3.53 +/-0.81 Mpc/h (comoving) and gamma=1.66 +/-0.12, while in a higher-redshift sample with z_eff=1.14 we find r_0=3.14 +/-0.72 Mpc/h and gamma=1.61 +/-0.11. We find that red, absorption-dominated, passively-evolving galaxies have a larger clustering scale length, r_0, and more prominent ``fingers of God'' than blue, emission-line, actively star-forming galaxies. Intrinsically brighter galaxies also cluster more strongly than fainter galaxies at z~1, with a significant luminosity-bias seen for galaxies fainter than M*. Our results are suggestive of evolution in the galaxy clustering within our survey volume and imply that the DEEP2 galaxies, with a median brightness one magnitude fainter than M* have an effective bias b=0.97 +/-0.13 if sigma_{8 DM}=1 today or b=1.20 +/-0.16 if sigma_{8 DM}=0.8 today. Given the strong luminosity-dependence in the bias that we measure at z~1, the galaxy bias at M* may be significantly greater. We note that our star-forming sample at z~1 has very similar selection criteria as the Lyman-break galaxies at z~3 and that our red, absorption-line sample displays a clustering strength comparable to the expected clustering of the Lyman-break galaxy descendants at z~1. Our results demonstrate that the clustering properties in the galaxy distribution seen in the local Universe were largely in place by z~1.

We present a study of variable faint optical sources discovered by comparing the Sloan Digital Sky Survey (SDSS) and the Palomar Observatory Sky Survey (POSS) catalogs. We use SDSS measurements to photometrically recalibrate several publicly available POSS catalogs; a piecewise recalibration in 100 arcmin2 patches generally results in an improvement of photometric accuracy (rms) by nearly a factor of two, compared to the original data. The POSS I magnitudes can be improved to ~0.15 mag accuracy, and POSS II magnitudes to \~0.10 mag accuracy. We use the recalibrated catalogs for the ~2,000 deg2 of sky in the SDSS Data Release 1 to construct a catalog of ~60,000 sources variable on time scales 10-50 years. A series of statistical tests based on the morphology of SDSS color-color diagrams, as well as visual comparison of images and comparison with repeated SDSS observations, demonstrate the robustness of the selection methods. We quantify the distribution of variable sources in the SDSS color-color diagrams, and the variability characteristics of quasars. We detect a turn-over in quasar structure function which suggests that the characteristic time scale for quasar variability is of the order one year. The long-term (>1 year) quasar variability decreases with luminosity and rest-frame wavelength similarly to the short-term (<1 year) behavior. We also demonstrate that candidate RR Lyrae stars trace the same halo structures, such as the Sgr dwarf tidal stream, that were discovered using repeated SDSS observations. We utilize the POSS-SDSS selected candidates to constrain the halo structure in the parts of sky for which repeated SDSS observations do not exist. (abridged)

We measure cosmological parameters using the three-dimensional power spectrum P(k) from over 200,000 galaxies in the Sloan Digital Sky Survey (SDSS) in combination with WMAP and other data. Our results are consistent with a ``vanilla'' flat adiabatic Lambda-CDM model without tilt (n=1), running tilt, tensor modes or massive neutrinos. Adding SDSS information more than halves the WMAP-only error bars on some parameters, tightening 1 sigma constraints on the Hubble parameter from h~0.74+0.18-0.07 to h~0.70+0.04-0.03, on the matter density from Omega_m~0.25+/-0.10 to Omega_m~0.30+/-0.04 (1 sigma) and on neutrino masses from <11 eV to <0.6 eV (95%). SDSS helps even more when dropping prior assumptions about curvature, neutrinos, tensor modes and the equation of state. Our results are in substantial agreement with the joint analysis of WMAP and the 2dF Galaxy Redshift Survey, which is an impressive consistency check with independent redshift survey data and analysis techniques. In this paper, we place particular emphasis on clarifying the physical origin of the constraints, i.e., what we do and do not know when using different data sets and prior assumptions. For instance, dropping the assumption that space is perfectly flat, the WMAP-only constraint on the measured age of the Universe tightens from t0~16.3+2.3-1.8 Gyr to t0~14.1+1.0-0.9 Gyr by adding SDSS and SN Ia data. Including tensors, running tilt, neutrino mass and equation of state in the list of free parameters, many constraints are still quite weak, but future cosmological measurements from SDSS and other sources should allow these to be substantially tightened.

We investigate the nature of the diffuse Galactic emission in the Wilkinson Microwave Anisotropy Probe (WMAP) temperature anisotropy data. Substantial dust-correlated emission is observed at all WMAP frequencies, far exceeding the expected thermal dust emission in the lowest frequency channels (23, 33, 41 GHz). The WMAP team (Bennett et al.) interpret this emission as dust-correlated synchrotron radiation, attributing the correlation to the natural association of relativistic electrons produced by SNae with massive star formation in dusty clouds, and deriving an upper limit of 5% on the contribution of Draine & Lazarian spinning dust at K-band (23 GHz). We pursue an alternative interpretation that much, perhaps most, of the dust-correlated emission at these frequencies is indeed spinning dust, and explore the spectral dependence on environment by considering a few specific objects as well as the full sky average. Models similar to Draine & Lazarian spinning dust provide a good fit to the full-sky data. The full-sky fit also requires a significant component with free-free spectrum uncorrelated with \Halpha, possibly hot (~million K) gas within 30 degrees of the Galactic center.

We measure the large-scale real-space power spectrum P(k) using a sample of 205,443 galaxies from the Sloan Digital Sky Survey, covering 2417 square degrees with mean redshift z~0.1. We employ a matrix-based method using pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 22 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.02 h/Mpc < k < 0.3h/Mpc. We pay particular attention to modeling, quantifying and correcting for potential systematic errors, nonlinear redshift distortions and the artificial red-tilt caused by luminosity-dependent bias. Our final result is a measurement of the real-space matter power spectrum P(k) up to an unknown overall multiplicative bias factor. Our calculations suggest that this bias factor is independent of scale to better than a few percent for k<0.1h/Mpc, thereby making our results useful for precision measurements of cosmological parameters in conjunction with data from other experiments such as the WMAP satellite. As a simple characterization of the data, our measurements are well fit by a flat scale-invariant adiabatic cosmological model with h Omega_m =0.201+/- 0.017 and L* galaxy sigma_8=0.89 +/- 0.02 when fixing the baryon fraction Omega_b/Omega_m=0.17 and the Hubble parameter h=0.72; cosmological interpretation is given in a companion paper.

We present a Principal Component Analysis (PCA)-based spectral classification, eta, for the first 5600 galaxies observed in the DEEP2 Redshift Survey. This parameter provides a very pronounced separation between absorption and emission dominated galaxy spectra - corresponding to passively evolving and actively star-forming galaxies in the survey respectively. In addition it is shown that despite the high resolution of the observed spectra, this parameter alone can be used to quite accurately reconstruct any given galaxy spectrum, suggesting there are not many `degrees of freedom' in the observed spectra of this galaxy population. It is argued that this form of classification, eta, will be particularly valuable in making future comparisons between high and low-redshift galaxy surveys for which very large spectroscopic samples are now readily available, particularly when used in conjunction with high-resolution spectral synthesis models which will be made public in the near future. We also discuss the relative advantages of this approach to distant galaxy classification compared to other methods such as colors and morphologies. Finally, we compare the classification derived here with that adopted for the 2dF Galaxy Redshift Survey and in so doing show that the two systems are very similar. This will be particularly useful in subsequent analyses when making comparisons between results from each of these surveys to study evolution in the galaxy populations and large-scale structure.

The distribution of early-type galaxy velocity dispersions, phi(sigma), is measured using a sample drawn from the SDSS database. Its shape differs significantly from that which one obtains by simply using the mean correlation between luminosity, L, and velocity dispersion, sigma, to transform the luminosity function into a velocity function: ignoring the scatter around the mean sigma-L relation is a bad approximation. An estimate of the contribution from late-type galaxies is also made, which suggests that phi(sigma) is dominated by early-type galaxies at velocities larger than ~ 200 km/s.

Polarized microwave emission from dust is an important foreground that may contaminate polarized CMB studies unless carefully accounted for. We discuss potential difficulties associated with this foreground, namely, the existence of different grain populations with very different emission/polarization properties and variations of the polarization yield with grain temperature. In particular, we discuss observational evidence in favor of rotational emission from tiny PAH particles with dipole moments, i.e. ``spinning dust'', and also consider magneto-dipole emission from strongly magnetized grains. We argue that in terms of polarization, the magneto-dipole emission may dominate even if its contribution to total emissivity is subdominant. Addressing polarized emission at frequencies larger than approsimately 100 GHz, we discuss the complications arising from the existence of dust components with different temperatures and possibly different alignment properties.

A magnitude limited sample of nearly 9000 early-type galaxies, in the redshift range 0.01 < z < 0.3, was selected from the Sloan Digital Sky Survey using morphological and spectral criteria. The Fundamental Plane relation in this sample is R_o ~ sigma^{1.49\pm 0.05} I_o^{-0.75\pm 0.01} in the r* band. It is approximately the same in the g*, i* and z* bands. Relative to the population at the median redshift in the sample, galaxies at lower and higher redshifts have evolved only little. If the Fundamental Plane is used to quantify this evolution then the apparent magnitude limit can masquerade as evolution; once this selection effect has been accounted for, the evolution is consistent with that of a passively evolving population which formed the bulk of its stars about 9 Gyrs ago. One of the principal advangtages of the SDSS sample over previous samples is that the galaxies in it lie in environments ranging from isolation in the field to the dense cores of clusters. The Fundamental Plane shows that galaxies in dense regions are slightly different from galaxies in less dense regions.

The colors and chemical abundances of early-type galaxies at redshifts z<0.3 are studied using a sample of nearly 9000 galaxies, selected from the Sloan Digital Sky Survey using morphological and spectral criteria. In this sample, redder galaxies have larger velocity dispersions. Color also correlates with magnitude and size, but these correlations are entirely due to the L-sigma and R_o-sigma relations: the primary correlation is color-sigma. The red light in early-type galaxies is, on average, slightly more centrally concentrated than the blue. Because of these color gradients, the strength of the color--magnitude relation depends on whether or not the colors are defined using a fixed metric aperture. Chemical evolution and star formation histories of early-type galaxies are investigated using co-added spectra of similar objects. At fixed sigma, the population at z ~ 0.2 had weaker Mg_2 and stronger H_beta absorption compared to the population at z ~ 0. It was also bluer. Comparison of these colors and line-strengths, and their evolution, with single-burst stellar population models suggests a formation time of 9 Gyrs ago, consistent with a Fundamental Plane analysis of this sample.

A sample of nearly 9000 early-type galaxies, in the redshift range 0.01 < z < 0.3, was selected from the Sloan Digital Sky Survey using morphological and spectral criteria. This paper describes how the sample was selected, presents examples of images and seeing corrected fits to the observed surface brightness profiles, describes our method for estimating K-corrections, and shows that the SDSS spectra are of sufficiently high quality to measure velocity dispersions accurately. It also provides catalogs of the measured photometric and spectroscopic parameters. In related papers, these data are used to study how early-type galaxy observables, including luminosity, effective radius, surface brightness, color, and velocity dispersion, are correlated with one another.

A magnitude limited sample of nearly 9000 early-type galaxies, in the redshift range 0.01 < z < 0.3, was selected from the Sloan Digital Sky Survey using morphological and spectral criteria. The sample was used to study how early-type galaxy observables, including luminosity L, effective radius R_o, surface brightness I_o, color, and velocity dispersion sigma, are correlated with one another. Measurement biases are understood with mock catalogs which reproduce all of the observed scaling relations and their dependences on fitting technique. At any given redshift, the intrinsic distribution of luminosities, sizes and velocity dispersions in our sample are all approximately Gaussian. A maximum likelihood analysis shows that sigma ~ L^{0.25\pm 0.012}, R_o ~ L^{0.63\pm 0.025}, and R_o ~ I^{-0.75\pm 0.02} in the r* band. In addition, the mass-to-light ratio within the effective radius scales as M_o/L ~ L^{0.14\pm 0.02} or M_o/L ~ M_o^{0.22\pm 0.05}, and galaxies with larger effective masses have smaller effective densities: Delta_o ~ M_o^{-0.52\pm 0.03}. These relations are approximately the same in the g*, i* and z* bands. Relative to the population at the median redshift in the sample, galaxies at lower and higher redshifts have evolved only little, with more evolution in the bluer bands. The luminosity function is consistent with weak passive luminosity evolution and a formation time of about 9 Gyrs ago.

A full-sky H-alpha map with 6 arcmin (FWHM) resolution is presented. This map is a composite of the Virginia Tech Spectral line Survey (VTSS) in the north and the Southern H-Alpha Sky Survey Atlas (SHASSA) in the south. The Wisconsin H-Alpha Mapper (WHAM) survey provides a stable zero-point over 3/4 of the sky on a 1 degree scale. This composite map can be used to provide limits on thermal bremsstrahlung (free-free emission) from ionized gas known to contaminate microwave-background data. The map is available on the WWW for public use.

We combine Sloan Digital Sky Survey spectra of 22,000 luminous, red, bulge-dominated galaxies to get high S/N average spectra in the rest-frame optical and ultraviolet (2600A to 7000A). The average spectra of these massive, quiescent galaxies are early-type with weak emission lines and with absorption lines indicating an apparent excess of alpha elements over solar abundance ratios. We make average spectra of subsamples selected by luminosity, environment and redshift. The average spectra are remarkable in their similarity. What variations do exist in the average spectra as a function of luminosity and environment are found to form a nearly one-parameter family in spectrum space. We present a high signal-to-noise ratio spectrum of the variation. We measure the properties of the variation with a modified version of the Lick index system and compare to model spectra from stellar population syntheses. The variation may be a combination of age and chemical abundance differences, but the conservative conclusion is that the quality of the data considerably exceeds the current state of the models.

The DEIMOS spectrograph has now been installed on the Keck-II telescope and commissioning is nearly complete. The DEEP2 Redshift Survey, which will take approximately 120 nights at the Keck Observatory over a three year period and has been designed to utilize the power of DEIMOS, began in the summer of 2002. The multiplexing power and high efficiency of DEIMOS enables us to target 1000 faint galaxies per clear night. Our goal is to gather high-quality spectra of \~60,000 galaxies with z>0.75 in order to study the properties and large scale clustering of galaxies at z ~ 1. The survey will be executed at high spectral resolution, $R=\lambda/\Delta \lambda \approx 5000$, allowing us to work between the bright OH sky emission lines and to infer linewidths for many of the target galaxies (for several thousand objects, we will obtain rotation curves as well). The linewidth data will facilitate the execution of the classical redshift-volume cosmological test, which can provide a precision measurement of the equation of state of the Universe. This talk reviews the project, summarizes our science goals and presents some early DEIMOS data.

We present a new puzzle involving Galactic microwave emission and attempt to resolve it. On one hand, a cross-correlation analysis of the WHAM H-alpha map with the Tenerife 10 and 15 GHz maps shows that the well-known DIRBE correlated microwave emission cannot be dominated by free-free emission. On the other hand, recent high resolution observations in the 8-10 GHz range with the Green Bank 140 ft telescope by Finkbeiner et al. failed to find the corresponding 8 sigma signal that would be expected in the simplest spinning dust models. So what physical mechanism is causing this ubiquitous dust-correlated emission? We argue for a model predicting that spinning dust is the culprit after all, but that the corresponding small grains are well correlated with the larger grains seen at 100 micron only on large angular scales. In support of this grain segregation model, we find the best spinning dust template to involve higher frequency maps in the range 12-60 micron, where emission from transiently heated small grains is important. Upcoming CMB experiments such as ground-based interferometers, MAP and Planck LFI with high resolution at low frequencies should allow a definitive test of this model.

A sample of nearly 9000 early-type galaxies, in the redshift range 0.01<z<0.3, was selected from the Sloan Digital Sky Survey using morphological and spectral criteria. The sample was used to study how early-type galaxy observables, including luminosity L, effective radius R, surface brightness I, color, and velocity dispersion V, are correlated with one another. Measurement biases are understood with mock catalogs which reproduce all of the observed scaling relations. At any given redshift, the intrinsic distribution of luminosities, sizes and velocity dispersions in our sample are all approximately Gaussian. In the r* band L ~ V^3.91, L ~ R^1.58, R ~ I^(-0.75), and the Fundamental Plane relation is R ~ V^(1.49) I^(-0.75). These relations are approximately the same in the g*, i* and z* bands. At fixed luminosity, the mass-to-light ratio scales as M/L ~ L^0.14. The g*-r* color scales as V^0.25. Color also correlates with magnitude and size, but these correlations are entirely due to the L-V and R-V relations. Chemical evolution and star formation histories are investigated using co-added spectra of similar objects in our sample. Chemical abundances correlate primarily with velocity dispersion. At fixed V, the higher redshift population is bluer, is weaker in Mg2, and is stronger in Hbeta than the population nearby. In addition, the population at higher redshifts is slightly more luminous. These differences are consistent with that of a passively evolving population which formed the bulk of its stars about 9 Gyrs ago. The Fundamental Plane suggests that galaxies in dense regions are slightly different from those in less dense regions, but the co-added spectra and color--magnitude relations show no statistically significant dependence on environment.

We present the first tentative detection of spinning dust emission from specific astronomical sources. All other detections in the current literature are statistical. The Green Bank 140 foot telescope was used to observe 10 dust clouds at 5, 8, and 10 GHz. In some cases, the observed emission was consistent with the negative spectral slope expected for free-free emission (thermal bremsstrahlung), but in two cases it was not. One HII region yields a rising spectrum, inconsistent with free-free or synchrotron emission at the 10 sigma level. One dark cloud (L 1622) has a similar spectrum with lower significance. Both spectra are consistent with electric dipole emission from rapidly rotating dust grains ("spinning dust"), as predicted by Draine & Lazarian.

The DEEP2/DEIMOS redshift survey, which will begin observing in the Spring of 2002, will gather high quality spectra on ~60000 galaxies in order to study the evolution of the properties and large scale clustering of galaxies at z~1. The data rate from DEIMOS will be in excess of 1 Gbyte/hour, and it is therefore imperative to employ completely automated data reduction techniques to manage the analysis. We here describe aspects of our data pipeline, which will make extensive use of B-splines for the sky-subtraction stage and for the combination of multiple frames.

An unsupervised software ``robot'' that automatically and robustly reduces and analyzes CCD observations of photometric standard stars is described. The robot measures extinction coefficients and other photometric parameters in real time and, more carefully, on the next day. It also reduces and analyzes data from an all-sky $10 \mu m$ camera to detect clouds; photometric data taken during cloudy periods are automatically rejected. The robot reports its findings back to observers and data analysts via the World-Wide Web. It can be used to assess photometricity, and to build data on site conditions. The robot's automated and uniform site monitoring represents a minimum standard for any observing site with queue scheduling, a public data archive, or likely participation in any future National Virtual Observatory.

From analysis of the DIRBE weekly averaged sky maps, we have detected substantial flux in the 60 and 100 micron channels in excess of expected zodiacal and Galactic emission. Two methods are used to separate zodiacal light from more distant emission. Method I makes use of the time-dependence of the North-South annual variation observed at the ecliptic poles. This method is robust against errors in the inter-planetary dust (IPD) model, but does not demonstrate isotropy of the background. Method II measures the ecliptic latitude dependence of the dust over a range of ecliptic latitudes at solar elongation e=90 degrees. This allows the excess to be determined in each week of the DIRBE mission for high redundancy, but the results depend weakly on the IPD model. Both methods give consistent results at 60 and 100 microns. The observed signal is consistent with an isotropic background at the level \nu I_\nu = 28.1 +- 1.8 +- 7(syst) nW/m^2/sr at 60 microns and 24.6 +- 2.5 +- 8 nW/m^2/sr at 100 microns. The IR excess detected at 140 and 240 microns by these methods agrees with previous measurements, which are thought to be the cosmic infra-red background (CIB). The detections at 60 and 100 microns are new. While this new excess is not necessarily the CIB, we have ruled out all known sources of emission in the solar system and Galaxy. We therefore tentatively interpret this signal as the CIB and consider the implications of such energy production from the viewpoint of star formation efficiency and black hole accretion efficiency. However, the IR excess exceeds limits on the CIB derived from the inferred opacity of the IGM to observed TeV photons, thus casting doubt on this interpretation. There is currently no satisfactory explanation for the 60-100 micron excess.

This chapter discusses Galactic dust and how its thermal emission confuses CMBR measurements. Interstellar dust grains are composed of many differing species, and observational evidence has only begun to disentagle their composition and sizes. Fortunately, their far-IR/millimeter emission is less complex. We describe how a very simple two-species model can describe the emission at 200 < nu < 3000 GHz to high precision. At lower frequencies, other non-thermal processes may dominate the emission from dust.

We present predicted full-sky maps of submillimeter and microwave emission from the diffuse interstellar dust in the Galaxy. These maps are extrapolated from the 100 micron emission and 100/240 micron flux ratio maps that Schlegel, Finkbeiner, & Davis (1998; SFD98) generated from IRAS and COBE/DIRBE data. Results are presented for a number of physically plausible emissivity models. We find that no power law emissivity function fits the FIRAS data from 200 - 2100 GHz. In this paper we provide a formalism for a multi-component model for the dust emission. A two-component model with a mixture of silicate and carbon-dominated grains (motivated by Pollack et al., 1994}) provides a fit to an accuracy of about 15% to all the FIRAS data over the entire high-latitude sky. Small systematic differences are found between the atomic and molecular phases of the ISM. Our predictions for the thermal (vibrational) emission from Galactic dust at \nu < 3000 GHz are available for general use. These full-sky predictions can be made at the DIRBE resolution of 40' or at the higher resolution of 6.1 arcmin from the SFD98 DIRBE-corrected IRAS maps.

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

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

We have constructed a full-sky map of the far-infrared suitable for measuring Galactic reddening and extinction (Schlegel, Finkbeiner & Davis 1998: SFD). The SFD map is based upon extensive re-analysis of data from the COBE/DIRBE and IRAS satellite missions. We demonstrate that the maps can correct for extinction problems in the APM galaxy survey. We also determine the most dust-free regions for conducting cosmic microwave background or soft X-ray experiments of extragalactic objects.

We present a full sky 100 micron map that is a reprocessed composite of the COBE/DIRBE and IRAS/ISSA maps, with the zodiacal foreground and confirmed point sources removed. Before using the ISSA maps, we remove the remaining artifacts from the IRAS scan pattern. Using the DIRBE 100 micron and 240 micron data, we have constructed a map of the dust temperature, so that the 100 micron map can be converted to a map proportional to dust column density. The result of these manipulations is a map with DIRBE-quality calibration and IRAS resolution. To generate the full sky dust maps, we must first remove zodiacal light contamination as well as a possible cosmic infrared background (CIB). This is done via a regression analysis of the 100 micron DIRBE map against the Leiden- Dwingeloo map of H_I emission, with corrections for the zodiacal light via a suitable expansion of the DIRBE 25 micron flux. For the 100 micron map, no significant CIB is detected. In the 140 micron and 240 micron maps, where the zodiacal contamination is weaker, we detect the CIB at surprisingly high flux levels of 32 \pm 13 nW/m^2/sr at 140 micron, and 17 \pm 4 nW/m^2/sr at 240 micron (95% confidence). This integrated flux is ~2 times that extrapolated from optical galaxies in the Hubble Deep Field. The primary use of these maps is likely to be as a new estimator of Galactic extinction. We demonstrate that the new maps are twice as accurate as the older Burstein-Heiles estimates in regions of low and moderate reddening. These dust maps will also be useful for estimating millimeter emission that contaminates CMBR experiments and for estimating soft X-ray absorption.