Source author record

Douglas P. Finkbeiner

Douglas P. Finkbeiner appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

42works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

42 published item(s)

preprint2022arXiv

Photometry on Structured Backgrounds: Local Pixelwise Infilling by Regression

Photometric pipelines struggle to estimate both the flux and flux uncertainty for stars in the presence of structured backgrounds such as filaments or clouds. However, it is exactly stars in these complex regions that are critical to understanding star formation and the structure of the interstellar medium. We develop a method, similar to Gaussian process regression, which we term local pixelwise infilling (LPI). Using a local covariance estimate, we predict the background behind each star and the uncertainty on that prediction in order to improve estimates of flux and flux uncertainty. We show the validity of our model on synthetic data and real dust fields. We further demonstrate that the method is stable even in the crowded field limit. While we focus on optical-IR photometry, this method is not restricted to those wavelengths. We apply this technique to the 34 billion detections in the second data release of the Dark Energy Camera Plane Survey (DECaPS2). In addition to removing many $>3σ$ outliers and improving uncertainty estimates by a factor of $\sim 2-3$ on nebulous fields, we also show that our method is well-behaved on uncrowded fields. The entirely post-processing nature of our implementation of LPI photometry allows it to easily improve the flux and flux uncertainty estimates of past as well as future surveys.

preprint2022arXiv

Star formation near the Sun is driven by expansion of the Local Bubble

For decades we have known that the Sun lies within the Local Bubble, a cavity of low-density, high-temperature plasma surrounded by a shell of cold, neutral gas and dust. However, the precise shape and extent of this shell, the impetus and timescale for its formation, and its relationship to nearby star formation have remained uncertain, largely due to low-resolution models of the local interstellar medium. Leveraging new spatial and dynamical constraints from the Gaia space mission, here we report an analysis of the 3D positions, shapes, and motions of dense gas and young stars within 200 pc of the Sun. We find that nearly all the star-forming complexes in the solar vicinity lie on the surface of the Local Bubble and that their young stars show outward expansion mainly perpendicular to the bubble's surface. Tracebacks of these young stars' motions support a scenario where the origin of the Local Bubble was a burst of stellar birth and then death (supernovae) taking place near the bubble's center beginning 14 Myr ago. The expansion of the Local Bubble created by the supernovae swept up the ambient interstellar medium into an extended shell that has now fragmented and collapsed into the most prominent nearby molecular clouds, in turn providing robust observational support for the theory of supernova-driven star formation.

preprint2020arXiv

A Compendium of Distances to Molecular Clouds in the Star Formation Handbook

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.

preprint2020arXiv

A Galactic-scale gas wave in the Solar Neighborhood

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.

preprint2020arXiv

Observing Strategy for the Legacy Surveys

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.

preprint2020arXiv

Photometric Biases in Modern Surveys

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$σ$ point source, a 10$σ$ 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.

preprint2019arXiv

A 3D Dust Map Based on Gaia, Pan-STARRS 1 and 2MASS

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.

preprint2016arXiv

A Synoptic Map of Halo Substructures from the Pan-STARRS1 3π Survey

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.

preprint2016arXiv

Maximizing Science in the Era of LSST: A Community-Based Study of Needed US Capabilities

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.

preprint2016arXiv

On Galactic density modeling in the presence of dust extinction

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.

preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

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

preprint2016arXiv

The DESI Experiment Part II: Instrument Design

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= λ/Δλ$ 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.

preprint2015arXiv

A Three-Dimensional Map of Milky-Way Dust

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.

preprint2015arXiv

Hypercalibration: A Pan-STARRS1-based recalibration of the Sloan Digital Sky Survey

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.

preprint2015arXiv

Results from the Wide-field Infrared Survey Explorer (WISE) Future Uses Session at the WISE at 5 Meeting

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.

preprint2015arXiv

The Characterization of the Gamma-Ray Signal from the Central Milky Way: A Compelling Case for Annihilating Dark Matter

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

preprint2014arXiv

An X-Ray Line from eXciting Dark Matter

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 $χχ\rightarrow χ^* χ^*$, followed by $χ^* \rightarrow χγ$. 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 $ϕ$. For suitable masses the annihilation $χχ\rightarrow ϕϕ$ followed by $ϕ\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.

preprint2014arXiv

Measuring Distances and Reddenings for a Billion Stars: Towards A 3D Dust Map from Pan-STARRS 1

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.

preprint2014arXiv

Sharper Fermi LAT Images: instrument response functions for an improved event selection

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.

preprint2013arXiv

Closing in on the Fermi Line with a New Observation Strategy

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.

preprint2013arXiv

Lacerta I and Cassiopeia III: Two luminous and distant Andromeda satellite dwarf galaxies found in the 3π Pan-STARRS1 survey

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° and 10.5° 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 (μ_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π Pan-STARRS1 survey, which could lead to the discovery of other distant Andromeda satellite dwarf galaxies.

preprint2013arXiv

Perseus I: A distant satellite dwarf galaxy of Andromeda

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π survey. Located 27.9° 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 (ε= 0.43 +0.15/-0.17). The late discovery of Perseus I is due to its fairly low surface brightness (μ_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.

preprint2012arXiv

Double Gamma-ray Lines from Unassociated Fermi-LAT Sources

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.

preprint2012arXiv

Evidence for Gamma-ray Jets in the Milky Way

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.

preprint2012arXiv

Is the 130 GeV Line Real? A Search for Systematics in the Fermi-LAT Data

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.

preprint2012arXiv

Strong Evidence for Gamma-ray Line Emission from the Inner Galaxy

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σ$ ($5.0σ/5.4σ$ 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 $α=1.2$ . The energy spectrum of this structure is consistent with a single spectral line (at energy $127.0\pm 2.0$ GeV with $χ^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 $χ^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 $γγ$ and $γZ$ (with $χ^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.

preprint2012arXiv

The DEEP2 Galaxy Redshift Survey: Design, Observations, Data Reduction, and Redshifts

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]

preprint2012arXiv

The Metallicity of the Monoceros Stream

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.

preprint2011arXiv

Consistent Scenarios for Cosmic-Ray Excesses from Sommerfeld-Enhanced Dark Matter Annihilation

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.

preprint2011arXiv

Magnetic Inelastic Dark Matter: Directional Signals Without a Directional Detector

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.

preprint2011arXiv

Measuring Reddening with SDSS Stellar Spectra and Recalibrating SFD

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.

preprint2011arXiv

Searching for Dark Matter in the CMB: A Compact Parameterization of Energy Injection from New Physics

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.

preprint2010arXiv

A statistical test of emission from unresolved point sources

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.

preprint2010arXiv

Giant Gamma-ray Bubbles from Fermi-LAT: AGN Activity or Bipolar Galactic Wind?

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.

preprint2010arXiv

Inelastic Dark Matter and DAMA/LIBRA: An Experimentum Crucis

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.

preprint2010arXiv

The Blue Tip of the Stellar Locus: Measuring Reddening with the SDSS

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.

preprint2010arXiv

The Electron Injection Spectrum Determined by Anomalous Excesses in Cosmic Ray, Gamma Ray, and Microwave Signals

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.

preprint2010arXiv

The Fermi Haze: A Gamma-Ray Counterpart to the Microwave 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.

preprint2008arXiv

CMB and 21-cm Signals for Dark Matter with a Long-Lived Excited State

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.

preprint2008arXiv

Nuclear scattering of dark matter coupled to a new light scalar

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.

preprint2006arXiv

The effect of FIR emission from SDSS galaxies on the SFD Galactic extinction map

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)

preprint2001arXiv

A photometricity and extinction monitor at the Apache Point Observatory

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