-
Improving Radial Velocities by Marginalizing over Stars and Sky: Achieving 30 m/s RV Precision for APOGEE in the Plate Era
Authors:
Andrew K. Saydjari,
Douglas P. Finkbeiner,
Adam J. Wheeler,
Jon A. Holtzman,
John C. Wilson,
Andrew R. Casey,
Sophia Sánchez-Maes,
Joel R. Brownstein,
David W. Hogg,
Michael R. Blanton
Abstract:
The radial velocity catalog from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) is unique in its simultaneously large volume and high precision as a result of its decade-long survey duration, multiplexing (600 fibers), and spectral resolution of $R \sim 22,500$. However, previous data reductions of APOGEE have not fully realized the potential radial velocity (RV) precision of…
▽ More
The radial velocity catalog from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) is unique in its simultaneously large volume and high precision as a result of its decade-long survey duration, multiplexing (600 fibers), and spectral resolution of $R \sim 22,500$. However, previous data reductions of APOGEE have not fully realized the potential radial velocity (RV) precision of the instrument. Here we present an RV catalog based on a new reduction of all 2.6 million visits of APOGEE DR17 and validate it against improved estimates for the theoretical RV performance. The core ideas of the new reduction are the simultaneous modeling of all components in the spectra, rather than a separate subtraction of point estimates for the sky, and a marginalization over stellar types, rather than a grid search for an optimum. We show that this catalog, when restricted to RVs measured with the same fiber, achieves noise-limited precision down to 30 m/s and delivers well-calibrated uncertainties. We also introduce a general method for calibrating fiber-to-fiber constant RV offsets and demonstrate its importance for high RV precision work in multi-fiber spectrographs. After calibration, we achieve 47 m/s RV precision on the combined catalog with RVs measured with different fibers. This degradation in precision relative to measurements with only a single fiber suggests that refining line spread function models should be a focus in SDSS-V to improve the fiber-unified RV catalog.
△ Less
Submitted 13 August, 2024;
originally announced August 2024.
-
Correcting Turbulence-induced Errors in Fiber Positioning for the Dark Energy Spectroscopic Instrument
Authors:
E. F. Schlafly,
J. Guy,
K. Honscheid,
S. Kent,
S. E. Koposov,
J. Aguilar,
S. Ahlen,
S. Bailey,
D. Brooks,
T. Claybaugh,
K. Dawson,
P. Doel,
K. Fanning,
D. P. Finkbeiner,
A. Font-Ribera,
J. E. Forero-Romero,
S. Gontcho A Gontcho,
G. Gutierrez,
D. Kirkby,
T. Kisner,
A. Kremin,
J. Lasker,
M. Landriau,
L. Le Guillou,
M. E. Levi
, et al. (15 additional authors not shown)
Abstract:
Highly-multiplexed, robotic, fiber-fed spectroscopic surveys are observing tens of millions of stars and galaxies. For many systems, accurate positioning relies on imaging the fibers in the focal plane and feeding that information back to the robotic positioners to correct their positions. Inhomogeneities and turbulence in the air between the focal plane and the imaging camera can affect the measu…
▽ More
Highly-multiplexed, robotic, fiber-fed spectroscopic surveys are observing tens of millions of stars and galaxies. For many systems, accurate positioning relies on imaging the fibers in the focal plane and feeding that information back to the robotic positioners to correct their positions. Inhomogeneities and turbulence in the air between the focal plane and the imaging camera can affect the measured positions of fibers, limiting the accuracy with which fibers can be placed on targets. For the Dark Energy Spectroscopic Instrument, we dramatically reduced the effect of turbulence on measurements of positioner locations in the focal plane by taking advantage of stationary positioners and the correlation function of the turbulence. We were able to reduce positioning errors from 7.3 microns to 3.5 microns, speeding the survey by 1.6% under typical conditions.
△ Less
Submitted 10 July, 2024;
originally announced July 2024.
-
Diffusion-HMC: Parameter Inference with Diffusion Model driven Hamiltonian Monte Carlo
Authors:
Nayantara Mudur,
Carolina Cuesta-Lazaro,
Douglas P. Finkbeiner
Abstract:
Diffusion generative models have excelled at diverse image generation and reconstruction tasks across fields. A less explored avenue is their application to discriminative tasks involving regression or classification problems. The cornerstone of modern cosmology is the ability to generate predictions for observed astrophysical fields from theory and constrain physical models from observations usin…
▽ More
Diffusion generative models have excelled at diverse image generation and reconstruction tasks across fields. A less explored avenue is their application to discriminative tasks involving regression or classification problems. The cornerstone of modern cosmology is the ability to generate predictions for observed astrophysical fields from theory and constrain physical models from observations using these predictions. This work uses a single diffusion generative model to address these interlinked objectives -- as a surrogate model or emulator for cold dark matter density fields conditional on input cosmological parameters, and as a parameter inference model that solves the inverse problem of constraining the cosmological parameters of an input field. The model is able to emulate fields with summary statistics consistent with those of the simulated target distribution. We then leverage the approximate likelihood of the diffusion generative model to derive tight constraints on cosmology by using the Hamiltonian Monte Carlo method to sample the posterior on cosmological parameters for a given test image. Finally, we demonstrate that this parameter inference approach is more robust to the addition of noise than baseline parameter inference networks.
△ Less
Submitted 8 May, 2024;
originally announced May 2024.
-
Cosmological Field Emulation and Parameter Inference with Diffusion Models
Authors:
Nayantara Mudur,
Carolina Cuesta-Lazaro,
Douglas P. Finkbeiner
Abstract:
Cosmological simulations play a crucial role in elucidating the effect of physical parameters on the statistics of fields and on constraining parameters given information on density fields. We leverage diffusion generative models to address two tasks of importance to cosmology -- as an emulator for cold dark matter density fields conditional on input cosmological parameters $Ω_m$ and $σ_8$, and as…
▽ More
Cosmological simulations play a crucial role in elucidating the effect of physical parameters on the statistics of fields and on constraining parameters given information on density fields. We leverage diffusion generative models to address two tasks of importance to cosmology -- as an emulator for cold dark matter density fields conditional on input cosmological parameters $Ω_m$ and $σ_8$, and as a parameter inference model that can return constraints on the cosmological parameters of an input field. We show that the model is able to generate fields with power spectra that are consistent with those of the simulated target distribution, and capture the subtle effect of each parameter on modulations in the power spectrum. We additionally explore their utility as parameter inference models and find that we can obtain tight constraints on cosmological parameters.
△ Less
Submitted 12 December, 2023;
originally announced December 2023.
-
Constructing Impactful Machine Learning Research for Astronomy: Best Practices for Researchers and Reviewers
Authors:
D. Huppenkothen,
M. Ntampaka,
M. Ho,
M. Fouesneau,
B. Nord,
J. E. G. Peek,
M. Walmsley,
J. F. Wu,
C. Avestruz,
T. Buck,
M. Brescia,
D. P. Finkbeiner,
A. D. Goulding,
T. Kacprzak,
P. Melchior,
M. Pasquato,
N. Ramachandra,
Y. -S. Ting,
G. van de Ven,
S. Villar,
V. A. Villar,
E. Zinger
Abstract:
Machine learning has rapidly become a tool of choice for the astronomical community. It is being applied across a wide range of wavelengths and problems, from the classification of transients to neural network emulators of cosmological simulations, and is shifting paradigms about how we generate and report scientific results. At the same time, this class of method comes with its own set of best pr…
▽ More
Machine learning has rapidly become a tool of choice for the astronomical community. It is being applied across a wide range of wavelengths and problems, from the classification of transients to neural network emulators of cosmological simulations, and is shifting paradigms about how we generate and report scientific results. At the same time, this class of method comes with its own set of best practices, challenges, and drawbacks, which, at present, are often reported on incompletely in the astrophysical literature. With this paper, we aim to provide a primer to the astronomical community, including authors, reviewers, and editors, on how to implement machine learning models and report their results in a way that ensures the accuracy of the results, reproducibility of the findings, and usefulness of the method.
△ Less
Submitted 19 October, 2023;
originally announced October 2023.
-
Simulation-based Inference for Exoplanet Atmospheric Retrieval: Insights from winning the Ariel Data Challenge 2023 using Normalizing Flows
Authors:
Mayeul Aubin,
Carolina Cuesta-Lazaro,
Ethan Tregidga,
Javier Viaña,
Cecilia Garraffo,
Iouli E. Gordon,
Mercedes López-Morales,
Robert J. Hargreaves,
Vladimir Yu. Makhnev,
Jeremy J. Drake,
Douglas P. Finkbeiner,
Phillip Cargile
Abstract:
Advancements in space telescopes have opened new avenues for gathering vast amounts of data on exoplanet atmosphere spectra. However, accurately extracting chemical and physical properties from these spectra poses significant challenges due to the non-linear nature of the underlying physics.
This paper presents novel machine learning models developed by the AstroAI team for the Ariel Data Challe…
▽ More
Advancements in space telescopes have opened new avenues for gathering vast amounts of data on exoplanet atmosphere spectra. However, accurately extracting chemical and physical properties from these spectra poses significant challenges due to the non-linear nature of the underlying physics.
This paper presents novel machine learning models developed by the AstroAI team for the Ariel Data Challenge 2023, where one of the models secured the top position among 293 competitors. Leveraging Normalizing Flows, our models predict the posterior probability distribution of atmospheric parameters under different atmospheric assumptions.
Moreover, we introduce an alternative model that exhibits higher performance potential than the winning model, despite scoring lower in the challenge. These findings highlight the need to reevaluate the evaluation metric and prompt further exploration of more efficient and accurate approaches for exoplanet atmosphere spectra analysis.
Finally, we present recommendations to enhance the challenge and models, providing valuable insights for future applications on real observational data. These advancements pave the way for more effective and timely analysis of exoplanet atmospheric properties, advancing our understanding of these distant worlds.
△ Less
Submitted 17 September, 2023;
originally announced September 2023.
-
Stellar Reddening Based Extinction Maps for Cosmological Applications
Authors:
Nayantara Mudur,
Core Francisco Park,
Douglas P Finkbeiner
Abstract:
Cosmological surveys must correct their observations for the reddening of extragalactic objects by Galactic dust. Existing dust maps, however, have been found to have spatial correlations with the large-scale structure of the Universe. Errors in extinction maps can propagate systematic biases into samples of dereddened extragalactic objects and into cosmological measurements such as correlation fu…
▽ More
Cosmological surveys must correct their observations for the reddening of extragalactic objects by Galactic dust. Existing dust maps, however, have been found to have spatial correlations with the large-scale structure of the Universe. Errors in extinction maps can propagate systematic biases into samples of dereddened extragalactic objects and into cosmological measurements such as correlation functions between foreground lenses and background objects and the primordial non-gaussianity parameter $f_{NL}$. Emission-based maps are contaminated by the cosmic infrared background, while maps inferred from stellar-reddenings suffer from imperfect removal of quasars and galaxies from stellar catalogs. Thus, stellar-reddening based maps using catalogs without extragalactic objects offer a promising path to making dust maps with minimal correlations with large-scale structure. We present two high-latitude integrated extinction maps based on stellar reddenings, with a point spread function of full-width half-maximum 6.1' and 15'. We employ a strict selection of catalog objects to filter out galaxies and quasars and measure the spatial correlation of our extinction maps with extragalactic structure. Our galactic extinction maps have reduced spatial correlation with large scale structure relative to most existing stellar-reddening based and emission-based extinction maps.
△ Less
Submitted 8 December, 2022;
originally announced December 2022.
-
Measuring the 8621 Å Diffuse Interstellar Band in Gaia DR3 RVS Spectra: Obtaining a Clean Catalog by Marginalizing over Stellar Types
Authors:
Andrew K. Saydjari,
Ana Sofía M. Uzsoy,
Catherine Zucker,
J. E. G. Peek,
Douglas P. Finkbeiner
Abstract:
Diffuse interstellar bands (DIBs) are broad absorption features associated with interstellar dust and can serve as chemical and kinematic tracers. Conventional measurements of DIBs in stellar spectra are complicated by residuals between observations and best-fit stellar models. To overcome this, we simultaneously model the spectrum as a combination of stellar, dust, and residual components, with f…
▽ More
Diffuse interstellar bands (DIBs) are broad absorption features associated with interstellar dust and can serve as chemical and kinematic tracers. Conventional measurements of DIBs in stellar spectra are complicated by residuals between observations and best-fit stellar models. To overcome this, we simultaneously model the spectrum as a combination of stellar, dust, and residual components, with full posteriors on the joint distribution of the components. This decomposition is obtained by modeling each component as a draw from a high-dimensional Gaussian distribution in the data-space (the observed spectrum) -- a method we call "Marginalized Analytic Data-space Gaussian Inference for Component Separation" (MADGICS). We use a data-driven prior for the stellar component, which avoids missing stellar features not well-modeled by synthetic spectra. This technique provides statistically rigorous uncertainties and detection thresholds, which are required to work in the low signal-to-noise regime that is commonplace for dusty lines of sight. We reprocess all public Gaia DR3 RVS spectra and present an improved 8621 Å DIB catalog, free of detectable stellar line contamination. We constrain the rest-frame wavelength to $8623.14 \pm 0.087$ Å (vacuum), find no significant evidence for DIBs in the Local Bubble from the $1/6^{\rm{th}}$ of RVS spectra that are public, and show unprecedented correlation with kinematic substructure in Galactic CO maps. We validate the catalog, its reported uncertainties, and biases using synthetic injection tests. We believe MADGICS provides a viable path forward for large-scale spectral line measurements in the presence of complex spectral contamination.
△ Less
Submitted 30 August, 2023; v1 submitted 7 December, 2022;
originally announced December 2022.
-
Can denoising diffusion probabilistic models generate realistic astrophysical fields?
Authors:
Nayantara Mudur,
Douglas P. Finkbeiner
Abstract:
Score-based generative models have emerged as alternatives to generative adversarial networks (GANs) and normalizing flows for tasks involving learning and sampling from complex image distributions. In this work we investigate the ability of these models to generate fields in two astrophysical contexts: dark matter mass density fields from cosmological simulations and images of interstellar dust.…
▽ More
Score-based generative models have emerged as alternatives to generative adversarial networks (GANs) and normalizing flows for tasks involving learning and sampling from complex image distributions. In this work we investigate the ability of these models to generate fields in two astrophysical contexts: dark matter mass density fields from cosmological simulations and images of interstellar dust. We examine the fidelity of the sampled cosmological fields relative to the true fields using three different metrics, and identify potential issues to address. We demonstrate a proof-of-concept application of the model trained on dust in denoising dust images. To our knowledge, this is the first application of this class of models to the interstellar medium.
△ Less
Submitted 22 November, 2022;
originally announced November 2022.
-
The First 3$π$ 3D Map of ISM Dust Temperature
Authors:
Ioana A. Zelko,
Douglas P. Finkbeiner,
Albert Lee,
Gregory Green
Abstract:
We present the first large-scale 3D map of interstellar dust temperature. We build upon existing 3D reddening maps derived from starlight absorption (Bayestar19), covering 3/4 of the sky. Starting with the column density for each of 500 million 3D voxels, we propose a temperature and emissivity power-law slope ($β$) for each of them, and integrate along the line of sight to synthesize an emission…
▽ More
We present the first large-scale 3D map of interstellar dust temperature. We build upon existing 3D reddening maps derived from starlight absorption (Bayestar19), covering 3/4 of the sky. Starting with the column density for each of 500 million 3D voxels, we propose a temperature and emissivity power-law slope ($β$) for each of them, and integrate along the line of sight to synthesize an emission map in five frequency bands observed by \emph{Planck} and \emph{IRAS}. The reconstructed emission map is constrained to match observations on a $10'$ scale, and does so with good fidelity. We produce 3D temperature maps at resolutions of $110', 55', $and $27'$. We assess performance on Cepheus, a dust cloud with two distinct components along the line of sight, and find distinct temperatures for the two components. We thus show that this methodology has enough precision to constrain clouds with different temperature along the line of sight up to $1-σ$ error. This would be an important result for dust frequency decorrelation foreground analysis for cosmic microwave background experiments, which would be impacted by a line-of-sight with varying temperature and magnetic field components. In addition to $T$ and $β$, we constrain the conversion factor between emission optical depth and reddening. This conversion factor is assumed to be constant in commonly used emission-based reddening maps. However, this work shows a factor of two variation that may prove significant for some applications.
△ Less
Submitted 14 November, 2022;
originally announced November 2022.
-
SARABANDE: 3/4 Point Correlation Functions with Fast Fourier Transforms
Authors:
James Sunseri,
Zachary Slepian,
Stephen Portillo,
Jiamin Hou,
Sule Kahraman,
Douglas P. Finkbeiner
Abstract:
We present a new $\texttt{python}$ package SARABANDE for measuring 3 & 4 Point Correlation Functions (3/4 PCFs) in $\mathcal{O}(N_{\rm g} \log N_{\rm g})$ time using Fast Fourier Transforms (FFTs), with $N_{\rm g}$ the number of grid points used for the FFT. SARABANDE can measure both projected and full 3 and 4 PCFs on gridded 2D and 3D datasets. The general technique is to generate suitable angul…
▽ More
We present a new $\texttt{python}$ package SARABANDE for measuring 3 & 4 Point Correlation Functions (3/4 PCFs) in $\mathcal{O}(N_{\rm g} \log N_{\rm g})$ time using Fast Fourier Transforms (FFTs), with $N_{\rm g}$ the number of grid points used for the FFT. SARABANDE can measure both projected and full 3 and 4 PCFs on gridded 2D and 3D datasets. The general technique is to generate suitable angular basis functions on an underlying grid, radially bin these to create kernels, and convolve these kernels with the original gridded data to obtain expansion coefficients about every point simultaneously. These coefficients are then combined to give us the 3/4 PCF as expanded in our basis. We apply SARABANDE to simulations of the Interstellar Medium (ISM) to show the results and scaling of calculating both the full and projected 3/4 PCFs.
△ Less
Submitted 25 October, 2022; v1 submitted 18 October, 2022;
originally announced October 2022.
-
The Dark Energy Camera Plane Survey 2 (DECaPS2): More Sky, Less Bias, and Better Uncertainties
Authors:
A. K. Saydjari,
E. F. Schlafly,
D. Lang,
A. M. Meisner,
G. M. Green,
C. Zucker,
I. Zelko,
J. S. Speagle,
T. Daylan,
A. Lee,
F. Valdes,
D. Schlegel,
D. P. Finkbeiner
Abstract:
Deep optical and near-infrared imaging of the entire Galactic plane is essential for understanding our Galaxy's stars, gas, and dust. The second data release of the DECam Plane Survey (DECaPS2) extends the five-band optical and near-infrared survey of the southern Galactic plane to cover $6.5\%$ of the sky, |b| < 10° and 6° > l > -124°, complementary to coverage by Pan-STARRS1. Typical single-expo…
▽ More
Deep optical and near-infrared imaging of the entire Galactic plane is essential for understanding our Galaxy's stars, gas, and dust. The second data release of the DECam Plane Survey (DECaPS2) extends the five-band optical and near-infrared survey of the southern Galactic plane to cover $6.5\%$ of the sky, |b| < 10° and 6° > l > -124°, complementary to coverage by Pan-STARRS1. Typical single-exposure effective depths, including crowding effects and other complications, are 23.5, 22.6, 22.1, 21.6, and 20.8 mag in $g$, $r$, $i$, $z$, and $Y$ bands, respectively, with around 1 arcsecond seeing. The survey comprises 3.32 billion objects built from 34 billion detections in 21.4 thousand exposures, totaling 260 hours open shutter time on the Dark Energy Camera (DECam) at Cerro Tololo. The data reduction pipeline features several improvements, including the addition of synthetic source injection tests to validate photometric solutions across the entire survey footprint. A convenient functional form for the detection bias in the faint limit was derived and leveraged to characterize the photometric pipeline performance. A new post-processing technique was applied to every detection to de-bias and improve uncertainty estimates of the flux in the presence of structured backgrounds, specifically targeting nebulosity. The images and source catalogs are publicly available at http://decaps.skymaps.info/.
△ Less
Submitted 26 July, 2022; v1 submitted 23 June, 2022;
originally announced June 2022.
-
Quantification of high dimensional non-Gaussianities and its implication to Fisher analysis in cosmology
Authors:
Core Francisco Park,
Erwan Allys,
Francisco Villaescusa-Navarro,
Douglas P. Finkbeiner
Abstract:
It is well known that the power spectrum is not able to fully characterize the statistical properties of non-Gaussian density fields. Recently, many different statistics have been proposed to extract information from non-Gaussian cosmological fields that perform better than the power spectrum. The Fisher matrix formalism is commonly used to quantify the accuracy with which a given statistic can co…
▽ More
It is well known that the power spectrum is not able to fully characterize the statistical properties of non-Gaussian density fields. Recently, many different statistics have been proposed to extract information from non-Gaussian cosmological fields that perform better than the power spectrum. The Fisher matrix formalism is commonly used to quantify the accuracy with which a given statistic can constrain the value of the cosmological parameters. However, these calculations typically rely on the assumption that the likelihood of the considered statistic follows a multivariate Gaussian distribution. In this work we follow Sellentin & Heavens (2017) and use two different statistical tests to identify non-Gaussianities in different statistics such as the power spectrum, bispectrum, marked power spectrum, and wavelet scatering transform (WST). We remove the non-Gaussian components of the different statistics and perform Fisher matrix calculations with the \textit{Gaussianized} statistics using Quijote simulations. We show that constraints on the parameters can change by a factor of $\sim 2$ in some cases. We show with simple examples how statistics that do not follow a multivariate Gaussian distribution can achieve artificially tight bounds on the cosmological parameters when using the Fisher matrix formalism. We think that the non-Gaussian tests used in this work represent a powerful tool to quantify the robustness of Fisher matrix calculations and their underlying assumptions. We release the code used to compute the power spectra, bispectra, and WST that can be run on both CPUs and GPUs.
△ Less
Submitted 11 April, 2022;
originally announced April 2022.
-
Photometry on Structured Backgrounds: Local Pixelwise Infilling by Regression
Authors:
Andrew K. Saydjari,
Douglas P. Finkbeiner
Abstract:
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…
▽ More
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.
△ Less
Submitted 18 January, 2022;
originally announced January 2022.
-
Star formation near the Sun is driven by expansion of the Local Bubble
Authors:
Catherine Zucker,
Alyssa A. Goodman,
João Alves,
Shmuel Bialy,
Michael Foley,
Joshua S. Speagle,
Josefa Großschedl,
Douglas P. Finkbeiner,
Andreas Burkert,
Diana Khimey,
Cameren Swiggum
Abstract:
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 interstell…
▽ More
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.
△ Less
Submitted 13 January, 2022;
originally announced January 2022.
-
Equivariant Wavelets: Fast Rotation and Translation Invariant Wavelet Scattering Transforms
Authors:
Andrew K. Saydjari,
Douglas P. Finkbeiner
Abstract:
Wavelet scattering networks, which are convolutional neural networks (CNNs) with fixed filters and weights, are promising tools for image analysis. Imposing symmetry on image statistics can improve human interpretability, aid in generalization, and provide dimension reduction. In this work, we introduce a fast-to-compute, translationally invariant and rotationally equivariant wavelet scattering ne…
▽ More
Wavelet scattering networks, which are convolutional neural networks (CNNs) with fixed filters and weights, are promising tools for image analysis. Imposing symmetry on image statistics can improve human interpretability, aid in generalization, and provide dimension reduction. In this work, we introduce a fast-to-compute, translationally invariant and rotationally equivariant wavelet scattering network (EqWS) and filter bank of wavelets (triglets). We demonstrate the interpretability and quantify the invariance/equivariance of the coefficients, briefly commenting on difficulties with implementing scale equivariance. On MNIST, we show that training on a rotationally invariant reduction of the coefficients maintains rotational invariance when generalized to test data and visualize residual symmetry breaking terms. Rotation equivariance is leveraged to estimate the rotation angle of digits and reconstruct the full rotation dependence of each coefficient from a single angle. We benchmark EqWS with linear classifiers on EMNIST and CIFAR-10/100, introducing a new second-order, cross-color channel coupling for the color images. We conclude by comparing the performance of an isotropic reduction of the scattering coefficients and RWST, a previous coefficient reduction, on an isotropic classification of magnetohydrodynamic simulations with astrophysical relevance.
△ Less
Submitted 22 April, 2021;
originally announced April 2021.
-
Classification of Magnetohydrodynamic Simulations using Wavelet Scattering Transforms
Authors:
Andrew K. Saydjari,
Stephen K. N. Portillo,
Zachary Slepian,
Sule Kahraman,
Blakesley Burkhart,
Douglas P. Finkbeiner
Abstract:
The complex interplay of magnetohydrodynamics, gravity, and supersonic turbulence in the interstellar medium (ISM) introduces non-Gaussian structure that can complicate comparison between theory and observation. We show that the Wavelet Scattering Transform (WST), in combination with linear discriminant analysis (LDA), is sensitive to non-Gaussian structure in 2D ISM dust maps. WST-LDA classifies…
▽ More
The complex interplay of magnetohydrodynamics, gravity, and supersonic turbulence in the interstellar medium (ISM) introduces non-Gaussian structure that can complicate comparison between theory and observation. We show that the Wavelet Scattering Transform (WST), in combination with linear discriminant analysis (LDA), is sensitive to non-Gaussian structure in 2D ISM dust maps. WST-LDA classifies magnetohydrodynamic (MHD) turbulence simulations with up to a 97\% true positive rate in our testbed of 8 simulations with varying sonic and Alfvénic Mach numbers. We present a side-by-side comparison with two other methods for non-Gaussian characterization, the Reduced Wavelet Scattering Transform (RWST) and the 3-Point Correlation Function (3PCF). We also demonstrate the 3D-WST-LDA and apply it to classification of density fields in position-position-velocity (PPV) space, where density correlations can be studied using velocity coherence as a proxy. WST-LDA is robust to common observational artifacts, such as striping and missing data, while also sensitive enough to extract the net magnetic field direction for sub-Alfvénic turbulent density fields. We include a brief analysis of the effect of point spread functions and image pixelization on 2D-WST-LDA applied to density fields, which informs the future goal of applying WST-LDA to 2D or 3D all-sky dust maps to extract hydrodynamic parameters of interest.
△ Less
Submitted 22 October, 2020;
originally announced October 2020.
-
Impact of Dust on Spectral Distortion Measurements of the Cosmic Microwave Background
Authors:
Ioana A. Zelko,
Douglas P. Finkbeiner
Abstract:
Spectral distortions of the cosmic microwave background (CMB) are sensitive to energy injection by exotic physics in the early universe. The proposed Primordial Inflation Explorer (PIXIE) mission has the raw sensitivity to provide meaningful limits on new physics, but only if foreground emission can be adequately modeled. We quantify the impact of interstellar dust on Compton $y$ and $μ$ measureme…
▽ More
Spectral distortions of the cosmic microwave background (CMB) are sensitive to energy injection by exotic physics in the early universe. The proposed Primordial Inflation Explorer (PIXIE) mission has the raw sensitivity to provide meaningful limits on new physics, but only if foreground emission can be adequately modeled. We quantify the impact of interstellar dust on Compton $y$ and $μ$ measurements by considering a range of grain size distributions and compositions constrained by theoretical and observational priors (Zelko & Finkbeiner 2020). We find that PIXIE can marginalize over a modest number of dust parameters and still recover $y$ and $μ$ estimates, though with increased uncertainty. As more foreground components are included (synchrotron, free-free), the estimates of $y$ degrade, and measurement of $μ$ in the range sometimes considered for the standard $Λ$CDM of $2\times10^{-8}$ becomes infeasible without ancillary low-frequency foreground information. An additional concern is dust absorption of the CMB monopole, a subtle effect that must be included. We quantify one form of model discrepancy error, finding that the error introduced by fitting our interstellar medium dust model with a modified blackbody is too large for CMB spectral distortions to be detectable. The greatest challenge may be the cosmic infrared background (CIB). We find that $μ$ and $y$ are extremely sensitive to modeling choices for the CIB, and quantify biases expected for a range of assumptions.
△ Less
Submitted 21 September, 2021; v1 submitted 13 October, 2020;
originally announced October 2020.
-
Implications of Grain Size Distribution and Composition for the Correlation Between Dust Extinction and Emissivity
Authors:
Ioana A. Zelko,
Douglas P. Finkbeiner
Abstract:
We study the effect of variations in dust size distribution and composition on the correlation between the spectral shape of extinction (parameterized by $R_{\textrm{V}}$) and far-infrared dust emissivity (parameterized by the power-law index $β$). Starting from the size distribution models proposed by Weingartner & Draine (2001), using the dust absorption and emission properties derived by Laor &…
▽ More
We study the effect of variations in dust size distribution and composition on the correlation between the spectral shape of extinction (parameterized by $R_{\textrm{V}}$) and far-infrared dust emissivity (parameterized by the power-law index $β$). Starting from the size distribution models proposed by Weingartner & Draine (2001), using the dust absorption and emission properties derived by Laor & Draine (1993) for carbonaceous and silicate grains, and by Li & Draine (2001) for PAH grains, we calculate the extinction and compare it with the reddening vector derived by Schlafly et al. (2016). An optimizer and an MCMC are used to explore the space of available parameters for the size distributions. We find that larger grains are correlated with high $R_{\textrm{V}}$. However, this trend is not enough to explain the emission-extinction correlation observed by Schlafly et al. (2016). For the $R_{\textrm{V}}-β$ correlation to arise, we need to impose explicit priors for the carbonaceous and silicate volume priors as functions of $R_{\textrm{V}}$. The results show that a composition with higher ratio of carbonaceous to silicate grains leads to higher $R_{\textrm{V}}$ and lower $β$. A relation between $E(\textrm{B}-\textrm{V})/τ_{353}$ and $R_{\textrm{V}}$ is apparent, with possible consequences for the recalibration of emission-based dust maps as a function of $R_{\textrm{V}}$.
△ Less
Submitted 18 May, 2021; v1 submitted 24 September, 2020;
originally announced September 2020.
-
Observing Strategy for the Legacy Surveys
Authors:
Kaylan J. Burleigh,
Martin Landriau,
Arjun Dey,
Dustin Lang,
David J. Schlegel,
Peter E. Nugent,
Robert Blum,
Joseph R. Findlay,
Douglas P. Finkbeiner,
David Herrera,
Klaus Honscheid,
Stéphanie Juneau,
Ian McGreer,
Aaron M. Meisner,
John Moustakas,
Adam D. Myers,
Anna Patej,
Edward F. Schlafly,
Francisco Valdes,
Alistair R. Walker,
Benjamin A. Weaver,
Christophe Yèche
Abstract:
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…
▽ More
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.
△ Less
Submitted 9 July, 2020; v1 submitted 13 February, 2020;
originally announced February 2020.
-
A Galactic-scale gas wave in the Solar Neighborhood
Authors:
João Alves,
Catherine Zucker,
Alyssa A. Goodman,
Joshua S. Speagle,
Stefan Meingast,
Thomas Robitaille,
Douglas P. Finkbeiner,
Edward F. Schlafly,
Gregory M. Green
Abstract:
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 thei…
▽ More
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.
△ Less
Submitted 23 January, 2020;
originally announced January 2020.
-
A Compendium of Distances to Molecular Clouds in the Star Formation Handbook
Authors:
Catherine Zucker,
Joshua S. Speagle,
Edward F. Schlafly,
Gregory M. Green,
Douglas P. Finkbeiner,
Alyssa Goodman,
João Alves
Abstract:
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 acc…
▽ More
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.
△ Less
Submitted 2 January, 2020;
originally announced January 2020.
-
A 3D Dust Map Based on Gaia, Pan-STARRS 1 and 2MASS
Authors:
Gregory M. Green,
Edward F. Schlafly,
Catherine Zucker,
Joshua S. Speagle,
Douglas P. Finkbeiner
Abstract:
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 nearb…
▽ More
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.
△ Less
Submitted 7 May, 2019;
originally announced May 2019.
-
Astro2020 Science White Paper: The Next Decade of Astroinformatics and Astrostatistics
Authors:
A. Siemiginowska,
G. Eadie,
I. Czekala,
E. Feigelson,
E. B. Ford,
V. Kashyap,
M. Kuhn,
T. Loredo,
M. Ntampaka,
A. Stevens,
A. Avelino,
K. Borne,
T. Budavari,
B. Burkhart,
J. Cisewski-Kehe,
F. Civano,
I. Chilingarian,
D. A. van Dyk,
G. Fabbiano,
D. P. Finkbeiner,
D. Foreman-Mackey,
P. Freeman,
A. Fruscione,
A. A. Goodman,
M. Graham
, et al. (27 additional authors not shown)
Abstract:
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 met…
▽ More
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.
△ Less
Submitted 15 March, 2019;
originally announced March 2019.
-
High-Dimensional Dust Mapping
Authors:
G. Zasowski,
D. P. Finkbeiner,
G. M. Green,
J. A. Kollmeier,
D. M. Nataf,
J. E. G. Peek,
E. F. Schlafly,
V. Silva Aguirre,
J. S. Speagle,
K. Tchernyshyov,
J. D. Trujillo,
C. Zucker
Abstract:
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)…
▽ More
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.
△ Less
Submitted 12 March, 2019;
originally announced March 2019.
-
Photometric Biases in Modern Surveys
Authors:
Stephen K. N. Portillo,
Joshua S. Speagle,
Douglas P. Finkbeiner
Abstract:
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…
▽ More
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.
△ Less
Submitted 14 February, 2020; v1 submitted 6 February, 2019;
originally announced February 2019.
-
A Large Catalog of Accurate Distances to Local Molecular Clouds: The Gaia DR2 Edition
Authors:
Catherine Zucker,
Joshua S. Speagle,
Edward F. Schlafly,
Gregory M. Green,
Douglas P. Finkbeiner,
Alyssa A. Goodman,
João Alves
Abstract:
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 para…
▽ More
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.
△ Less
Submitted 19 April, 2019; v1 submitted 4 February, 2019;
originally announced February 2019.
-
They Might Be Giants: An Efficient Color-Based Selection of Red Giant Stars
Authors:
Charlie Conroy,
Ana Bonaca,
Rohan P. Naidu,
Daniel J. Eisenstein,
Benjamin D. Johnson,
Aaron Dotter,
Douglas P. Finkbeiner
Abstract:
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 prim…
▽ More
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.
△ Less
Submitted 28 June, 2018; v1 submitted 15 May, 2018;
originally announced May 2018.
-
A Color-locus Method for Mapping $R_V$ Using Ensembles of Stars
Authors:
Albert Lee,
Gregory M. Green,
Edward F. Schlafly,
Douglas P. Finkbeiner,
William Burgett,
Ken Chambers,
Heather Flewelling,
Klaus Hodapp,
Nick Kaiser,
Rolf-Peter Kudritzki,
Eugene Magnier,
Nigel Metcalfe,
Richard Wainscoat,
Christopher Waters
Abstract:
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 als…
▽ More
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
△ Less
Submitted 24 March, 2018;
originally announced March 2018.
-
Mapping Distances Across the Perseus Molecular Cloud Using CO Observations, Stellar Photometry, and Gaia DR2 Parallax Measurements
Authors:
Catherine Zucker,
Edward F. Schlafly,
Joshua S. Speagle,
Gregory M. Green,
Stephen K. N. Portillo,
Douglas P. Finkbeiner,
Alyssa A. Goodman
Abstract:
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 techn…
▽ More
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.
△ Less
Submitted 17 October, 2018; v1 submitted 23 March, 2018;
originally announced March 2018.
-
Developing the 3-Point Correlation Function For the Turbulent Interstellar Medium
Authors:
Stephen K. N. Portillo,
Zachary Slepian,
Blakesley Burkhart,
Sule Kahraman,
Douglas P. Finkbeiner
Abstract:
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énic Mach numbers of turbulence. Here we show that int…
▽ More
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énic 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énic Mach numbers over a range from sub to super-sonic and sub to super-Alfvé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.
△ Less
Submitted 1 October, 2018; v1 submitted 27 November, 2017;
originally announced November 2017.
-
The DECam Plane Survey: Optical photometry of two billion objects in the southern Galactic plane
Authors:
E. F. Schlafly,
G. M. Green,
D. Lang,
T. Daylan,
D. P. Finkbeiner,
A. Lee,
A. M. Meisner,
D. Schlegel,
F. Valdes
Abstract:
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…
▽ More
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.
△ Less
Submitted 3 October, 2017;
originally announced October 2017.
-
Probing the small-scale structure in strongly lensed systems via transdimensional inference
Authors:
Tansu Daylan,
Francis-Yan Cyr-Racine,
Ana Diaz Rivero,
Cora Dvorkin,
Douglas P. Finkbeiner
Abstract:
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…
▽ More
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.
△ Less
Submitted 5 March, 2018; v1 submitted 19 June, 2017;
originally announced June 2017.
-
Improved Point Source Detection in Crowded Fields using Probabilistic Cataloging
Authors:
Stephen K. N. Portillo,
Benjamin C. G. Lee,
Tansu Daylan,
Douglas P. Finkbeiner
Abstract:
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 catalog…
▽ More
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.
△ Less
Submitted 7 August, 2017; v1 submitted 3 March, 2017;
originally announced March 2017.
-
The Pan-STARRS1 Surveys
Authors:
K. C. Chambers,
E. A. Magnier,
N. Metcalfe,
H. A. Flewelling,
M. E. Huber,
C. Z. Waters,
L. Denneau,
P. W. Draper,
D. Farrow,
D. P. Finkbeiner,
C. Holmberg,
J. Koppenhoefer,
P. A. Price,
A. Rest,
R. P. Saglia,
E. F. Schlafly,
S. J. Smartt,
W. Sweeney,
R. J. Wainscoat,
W. S. Burgett,
S. Chastel,
T. Grav,
J. N. Heasley,
K. W. Hodapp,
R. Jedicke
, et al. (101 additional authors not shown)
Abstract:
Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3π$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$σ$ point source limiting sensitivities in the stacked 3$π$ 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…
▽ More
Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3π$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$σ$ point source limiting sensitivities in the stacked 3$π$ 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 ($ Δra, Δ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.
△ Less
Submitted 28 January, 2019; v1 submitted 16 December, 2016;
originally announced December 2016.
-
Pan-STARRS Photometric and Astrometric Calibration
Authors:
Eugene. A. Magnier,
Edward. F. Schlafly,
Douglas P. Finkbeiner,
J. L. Tonry,
B. Goldman,
S. Röser,
E. Schilbach,
K. C. Chambers,
H. A. Flewelling,
M. E. Huber,
P. A. Price,
W. E. Sweeney,
C. Z. Waters,
L. Denneau,
P. Draper,
K. W. Hodapp,
R. Jedicke,
R. -P. Kudritzki,
N. Metcalfe,
C. W. Stubbs,
R. J. Wainscoast
Abstract:
We present the details of the photometric and astrometric calibration of the Pan-STARRS1 $3π$ 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 $δ$ = -30 degrees. The astrometric calib…
▽ More
We present the details of the photometric and astrometric calibration of the Pan-STARRS1 $3π$ 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 $δ$ = -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.
△ Less
Submitted 26 January, 2019; v1 submitted 15 December, 2016;
originally announced December 2016.
-
Mapping the Extinction Curve in 3D: Structure on Kiloparsec Scales
Authors:
E. F. Schlafly,
J. E. G. Peek,
D. P. Finkbeiner,
G. M. Green
Abstract:
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-dim…
▽ More
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.
△ Less
Submitted 23 March, 2017; v1 submitted 8 December, 2016;
originally announced December 2016.
-
The DESI Experiment Part II: Instrument Design
Authors:
DESI Collaboration,
Amir Aghamousa,
Jessica Aguilar,
Steve Ahlen,
Shadab Alam,
Lori E. Allen,
Carlos Allende Prieto,
James Annis,
Stephen Bailey,
Christophe Balland,
Otger Ballester,
Charles Baltay,
Lucas Beaufore,
Chris Bebek,
Timothy C. Beers,
Eric F. Bell,
José Luis Bernal,
Robert Besuner,
Florian Beutler,
Chris Blake,
Hannes Bleuler,
Michael Blomqvist,
Robert Blum,
Adam S. Bolton,
Cesar Briceno
, et al. (268 additional authors not shown)
Abstract:
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…
▽ More
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.
△ Less
Submitted 13 December, 2016; v1 submitted 31 October, 2016;
originally announced November 2016.
-
The DESI Experiment Part I: Science,Targeting, and Survey Design
Authors:
DESI Collaboration,
Amir Aghamousa,
Jessica Aguilar,
Steve Ahlen,
Shadab Alam,
Lori E. Allen,
Carlos Allende Prieto,
James Annis,
Stephen Bailey,
Christophe Balland,
Otger Ballester,
Charles Baltay,
Lucas Beaufore,
Chris Bebek,
Timothy C. Beers,
Eric F. Bell,
José Luis Bernal,
Robert Besuner,
Florian Beutler,
Chris Blake,
Hannes Bleuler,
Michael Blomqvist,
Robert Blum,
Adam S. Bolton,
Cesar Briceno
, et al. (268 additional authors not shown)
Abstract:
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…
▽ More
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.
△ Less
Submitted 13 December, 2016; v1 submitted 31 October, 2016;
originally announced November 2016.
-
Maximizing Science in the Era of LSST: A Community-Based Study of Needed US Capabilities
Authors:
Joan Najita,
Beth Willman,
Douglas P. Finkbeiner,
Ryan J. Foley,
Suzanne Hawley,
Jeffrey A. Newman,
Gregory Rudnick,
Joshua D. Simon,
David Trilling,
Rachel Street,
Adam Bolton,
Ruth Angus,
Eric F. Bell,
Derek Buzasi,
David Ciardi,
James R. A. Davenport,
Will Dawson,
Mark Dickinson,
Alex Drlica-Wagner,
Jay Elias,
Dawn Erb,
Lori Feaga,
Wen-fai Fong,
Eric Gawiser,
Mark Giampapa
, et al. (26 additional authors not shown)
Abstract:
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., o…
▽ More
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.
△ Less
Submitted 5 October, 2016;
originally announced October 2016.
-
A Synoptic Map of Halo Substructures from the Pan-STARRS1 3π Survey
Authors:
Edouard J. Bernard,
Annette M. N. Ferguson,
Edward F. Schlafly,
Nicolas F. Martin,
Hans-Walter Rix,
Eric F. Bell,
Douglas P. Finkbeiner,
Bertrand Goldman,
David Martinez-Delgado,
Branimir Sesar,
Rosemary F. G. Wyse,
William S. Burgett,
Kenneth C. Chambers,
Peter W. Draper,
Klaus W. Hodapp,
Nicholas Kaiser,
Rolf-Peter Kudritzki,
Eugene A. Magnier,
Nigel Metcalfe,
Richard J. Wainscoat,
Christopher Waters
Abstract:
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 an…
▽ More
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.
△ Less
Submitted 22 August, 2016; v1 submitted 20 July, 2016;
originally announced July 2016.
-
Inference of Unresolved Point Sources At High Galactic Latitudes Using Probabilistic Catalogs
Authors:
Tansu Daylan,
Stephen K. N. Portillo,
Douglas P. Finkbeiner
Abstract:
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 pe…
▽ More
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.
△ Less
Submitted 9 March, 2017; v1 submitted 15 July, 2016;
originally announced July 2016.
-
The Optical-Infrared Extinction Curve and its Variation in the Milky Way
Authors:
E. F. Schlafly,
A. M. Meisner,
A. M. Stutz,
J. Kainulainen,
J. E. G. Peek,
K. Tchernyshyov,
H. -W. Rix,
D. P. Finkbeiner,
K. R. Covey,
G. M. Green,
E. F. Bell,
W. S. Burgett,
K. C. Chambers,
P. W. Draper,
H. Flewelling,
K. W. Hodapp,
N. Kaiser,
E. A. Magnier,
N. F. Martin,
N. Metcalfe,
R. J. Wainscoat,
C. Waters
Abstract:
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…
▽ More
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.
△ Less
Submitted 11 February, 2016;
originally announced February 2016.
-
Hypercalibration: A Pan-STARRS1-based recalibration of the Sloan Digital Sky Survey
Authors:
Douglas P. Finkbeiner,
Edward F. Schlafly,
David J. Schlegel,
Nikhil Padmanabhan,
Mario Juric,
William S. Burgett,
Kenneth C. Chambers,
Larry Denneau,
Peter W. Draper,
Heather Flewelling,
Klaus W. Hodapp,
Nick Kaiser,
E. A. Magnier,
N. Metcalfe,
Jeffrey S. Morgan,
Paul A. Price,
Christopher W. Stubbs,
John L. Tonry
Abstract:
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 ar…
▽ More
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.
△ Less
Submitted 3 December, 2015;
originally announced December 2015.
-
On Galactic density modeling in the presence of dust extinction
Authors:
Jo Bovy,
Hans-Walter Rix,
Gregory M. Green,
Edward F. Schlafly,
Douglas P. Finkbeiner
Abstract:
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 in…
▽ More
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.
△ Less
Submitted 6 June, 2016; v1 submitted 22 September, 2015;
originally announced September 2015.
-
A Three-Dimensional Map of Milky-Way Dust
Authors:
Gregory M. Green,
Edward F. Schlafly,
Douglas P. Finkbeiner,
Hans-Walter Rix,
Nicolas Martin,
William Burgett,
Peter W. Draper,
Heather Flewelling,
Klaus Hodapp,
Nicholas Kaiser,
Rolf Peter Kudritzki,
Eugene Magnier,
Nigel Metcalfe,
Paul Price,
John Tonry,
Richard Wainscoat
Abstract:
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 dist…
▽ More
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.
△ Less
Submitted 3 July, 2015;
originally announced July 2015.
-
Results from the Wide-field Infrared Survey Explorer (WISE) Future Uses Session at the WISE at 5 Meeting
Authors:
Jacqueline K. Faherty,
K. Alatalo,
L. D. Anderson,
Roberto J. Assef,
Daniella C. Bardalez Gagliuffi,
Megan Barry,
Dominic J. Benford,
Maciej Bilicki,
Ben Burningham,
Damian J. Christian,
Michael C. Cushing,
Peter R. Eisenhardt,
Martin Elvisx,
S. B. Fajardo-Acosta,
Douglas P. Finkbeiner,
William J. Fischer,
William J. Forrest,
John Fowler,
Jonathan P. Gardner,
Christopher R. Gelino,
V Gorjian,
Carl J. Grillmair,
Mariusz Gromadzki,
Kendall P. Hall,
Zeljko Ivezi'c
, et al. (29 additional authors not shown)
Abstract:
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 w…
▽ More
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.
△ Less
Submitted 8 May, 2015;
originally announced May 2015.
-
3D Dust Mapping Reveals that Orion Forms Part of a Large Ring of Dust
Authors:
E. F. Schlafly,
G. Green,
D. P. Finkbeiner,
H. -W. Rix,
W. S. Burgett,
K. C. Chambers,
P. W. Draper,
N. Kaiser,
N. F. Martin,
N. Metcalfe,
J. S. Morgan,
P. A. Price,
J. L. Tonry,
R. J. Wainscoat,
C. Waters
Abstract:
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), p…
▽ More
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.
△ Less
Submitted 19 November, 2014;
originally announced November 2014.
-
Sharper Fermi LAT Images: instrument response functions for an improved event selection
Authors:
Stephen K. N. Portillo,
Douglas P. Finkbeiner
Abstract:
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 reconst…
▽ More
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.
△ Less
Submitted 28 September, 2014; v1 submitted 2 June, 2014;
originally announced June 2014.
-
A Map of Dust Reddening to 4.5 kpc from Pan-STARRS1
Authors:
E. F. Schlafly,
G. Green,
D. P. Finkbeiner,
M. Juric,
H. -W. Rix,
N. F. Martin,
W. S. Burgett,
K. C. Chambers,
P. W. Draper,
K. W. Hodapp,
N. Kaiser,
R. -P. Kudritzki,
E. A. Magnier,
N. Metcalfe,
J. S. Morgan,
P. A. Price,
C. W. Stubbs,
J. L. Tonry,
R. J. Wainscoat,
C. Waters
Abstract:
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 o…
▽ More
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).
△ Less
Submitted 12 May, 2014;
originally announced May 2014.