-
The DECam Ecliptic Exploration Project (DEEP) VI: first multi-year observations of trans-Neptunian objects
Authors:
Hayden Smotherman,
Pedro H. Bernardinelli,
Stephen K. N. Portillo,
Andrew J. Connolly,
J. Bryce Kalmbach,
Steven Stetzler,
Mario Juric,
Dino Bektesvic,
Zachary Langford,
Fred C. Adams,
William J. Oldroyd,
Matthew J. Holman,
Colin Orion Chandler,
Cesar Fuentes,
David W. Gerdes,
Hsing Wen Lin,
Larissa Markwardt,
Andrew McNeill,
Michael Mommert,
Kevin J. Napier,
Matthew J. Payne,
Darin Ragozzine,
Andrew S. Rivkin,
Hilke Schlichting,
Scott S. Sheppard
, et al. (3 additional authors not shown)
Abstract:
We present the first set of trans-Neptunian objects (TNOs) observed on multiple nights in data taken from the DECam Ecliptic Exploration Project (DEEP). Of these 110 TNOs, 105 do not coincide with previously known TNOs and appear to be new discoveries. Each individual detection for our objects resulted from a digital tracking search at TNO rates of motion, using two to four hour exposure sets, and…
▽ More
We present the first set of trans-Neptunian objects (TNOs) observed on multiple nights in data taken from the DECam Ecliptic Exploration Project (DEEP). Of these 110 TNOs, 105 do not coincide with previously known TNOs and appear to be new discoveries. Each individual detection for our objects resulted from a digital tracking search at TNO rates of motion, using two to four hour exposure sets, and the detections were subsequently linked across multiple observing seasons. This procedure allows us to find objects with magnitudes $m_{VR} \approx 26$. The object discovery processing also included a comprehensive population of objects injected into the images, with a recovery and linking rate of at least $94\%$. The final orbits were obtained using a specialized orbit fitting procedure that accounts for the positional errors derived from the digital tracking procedure. Our results include robust orbits and magnitudes for classical TNOs with absolute magnitudes $H \sim 10$, as well as a dynamically detached object found at 76 au (semi-major axis $a\approx 77 \, \mathrm{au}$). We find a disagreement between our population of classical TNOs and the CFEPS-L7 three component model for the Kuiper belt.
△ Less
Submitted 5 October, 2023;
originally announced October 2023.
-
The DECam Ecliptic Exploration Project (DEEP) III: Survey characterization and simulation methods
Authors:
Pedro H. Bernardinelli,
Hayden Smotherman,
Zachary Langford,
Stephen K. N. Portillo,
Andrew J. Connolly,
J. Bryce Kalmbach,
Steven Stetzler,
Mario Juric,
William J. Oldroyd,
Hsing Wen Lin,
Fred C. Adams,
Colin Orion Chandler,
Cesar Fuentes,
David W. Gerdes,
Matthew J. Holman,
Larissa Markwardt,
Andrew McNeill,
Michael Mommert,
Kevin J. Napier,
Matthew J. Payne,
Darin Ragozzine,
Andrew S. Rivkin,
Hilke Schlichting,
Scott S. Sheppard,
Ryder Strauss
, et al. (2 additional authors not shown)
Abstract:
We present a detailed study of the observational biases of the DECam Ecliptic Exploration Project's (DEEP) B1 data release and survey simulation software that enables direct statistical comparisons between models and our data. We inject a synthetic population of objects into the images, and then subsequently recover them in the same processing as our real detections. This enables us to characteriz…
▽ More
We present a detailed study of the observational biases of the DECam Ecliptic Exploration Project's (DEEP) B1 data release and survey simulation software that enables direct statistical comparisons between models and our data. We inject a synthetic population of objects into the images, and then subsequently recover them in the same processing as our real detections. This enables us to characterize the survey's completeness as a function of apparent magnitudes and on-sky rates of motion. We study the statistically optimal functional form for the magnitude, and develop a methodology that can estimate the magnitude and rate efficiencies for all survey's pointing groups simultaneously. We have determined that our peak completeness is on average 80\% in each pointing group, and our magnitude drops to $25\%$ of this value at $m_{25} = 26.22$. We describe the freely available survey simulation software and its methodology. We conclude by using it to infer that our effective search area for objects at 40 au is $14.8°^2$, and that our lack of dynamically cold distant objects means that there at most $8\times 10^3$ objects with $60 < a < 80$ au and absolute magnitudes $H \leq 8$.
△ Less
Submitted 5 October, 2023;
originally announced October 2023.
-
The DECam Ecliptic Exploration Project (DEEP): V. The Absolute Magnitude Distribution of the Cold Classical Kuiper Belt
Authors:
Kevin J. Napier,
Hsing-Wen Lin,
David W. Gerdes,
Fred C. Adams,
Anna M. Simpson,
Matthew W. Porter,
Katherine G. Weber,
Larissa Markwardt,
Gabriel Gowman,
Hayden Smotherman,
Pedro H. Bernardinelli,
Mario Jurić,
Andrew J. Connolly,
J. Bryce Kalmbach,
Stephen K. N. Portillo,
David E. Trilling,
Ryder Strauss,
William J. Oldroyd,
Chadwick A. Trujillo,
Colin Orion Chandler,
Matthew J. Holman,
Hilke E. Schlichting,
Andrew McNeill,
the DEEP Collaboration
Abstract:
The DECam Ecliptic Exploration Project (DEEP) is a deep survey of the trans-Neptunian solar system being carried out on the 4-meter Blanco telescope at Cerro Tololo Inter-American Observatory in Chile using the Dark Energy Camera (DECam). By using a shift-and-stack technique to achieve a mean limiting magnitude of $r \sim 26.2$, DEEP achieves an unprecedented combination of survey area and depth,…
▽ More
The DECam Ecliptic Exploration Project (DEEP) is a deep survey of the trans-Neptunian solar system being carried out on the 4-meter Blanco telescope at Cerro Tololo Inter-American Observatory in Chile using the Dark Energy Camera (DECam). By using a shift-and-stack technique to achieve a mean limiting magnitude of $r \sim 26.2$, DEEP achieves an unprecedented combination of survey area and depth, enabling quantitative leaps forward in our understanding of the Kuiper Belt populations. This work reports results from an analysis of twenty 3 sq.\ deg.\ DECam fields along the invariable plane. We characterize the efficiency and false-positive rates for our moving-object detection pipeline, and use this information to construct a Bayesian signal probability for each detected source. This procedure allows us to treat all of our Kuiper Belt Object (KBO) detections statistically, simultaneously accounting for efficiency and false positives. We detect approximately 2300 candidate sources with KBO-like motion at S/N $>6.5$. We use a subset of these objects to compute the luminosity function of the Kuiper Belt as a whole, as well as the Cold Classical (CC) population. We also investigate the absolute magnitude ($H$) distribution of the CCs, and find consistency with both an exponentially tapered power-law, which is predicted by streaming instability models of planetesimal formation, and a rolling power law. Finally, we provide an updated mass estimate for the Cold Classical Kuiper Belt of $M_{CC}(H_r < 12) = 0.0017^{+0.0010}_{-0.0004} M_{\oplus}$, assuming albedo $p = 0.15$ and density $ρ= 1$ g cm$^{-3}$.
△ Less
Submitted 18 September, 2023;
originally announced September 2023.
-
PCAT-DE: Reconstructing point-like and diffuse signals in astronomical images using spatial and spectral information
Authors:
Richard M. Feder,
Victoria Butler,
Tansu Daylan,
Stephen K. N. Portillo,
Jack Sayers,
Benjamin J. Vaughan,
Catalina V. Zamora,
Michael Zemcov
Abstract:
Observational data from astronomical imaging surveys contain information about a variety of source populations and environments, and its complexity will increase substantially as telescopes become more sensitive. Even for existing observations, measuring the correlations between point-like and diffuse emission can be crucial to correctly inferring the properties of any individual component. For th…
▽ More
Observational data from astronomical imaging surveys contain information about a variety of source populations and environments, and its complexity will increase substantially as telescopes become more sensitive. Even for existing observations, measuring the correlations between point-like and diffuse emission can be crucial to correctly inferring the properties of any individual component. For this task information is typically lost, either because of conservative data cuts, aggressive filtering or incomplete treatment of contaminated data. We present the code PCAT-DE, an extension of probabilistic cataloging designed to simultaneously model point-like and diffuse signals. This work incorporates both explicit spatial templates and a set of non-parametric Fourier component templates into a forward model of astronomical images, reducing the number of processing steps applied to the observed data. Using synthetic Herschel-SPIRE multiband observations, we demonstrate that point source and diffuse emission can be reliably separated and measured. We present two applications of this model. For the first, we perform point source detection/photometry in the presence of galactic cirrus and demonstrate that cosmic infrared background (CIB) galaxy counts can be recovered in cases of significant contamination. In the second we show that the spatially extended thermal Sunyaev-Zel'dovich (tSZ) effect signal can be reliably measured even when it is subdominant to the point-like emission from individual galaxies.
△ Less
Submitted 19 July, 2023;
originally announced July 2023.
-
Measurement of the Relativistic Sunyaev-Zeldovich Corrections in RX J1347.5-1145
Authors:
Victoria Butler,
Richard M. Feder,
Tansu Daylan,
Adam B. Mantz,
Dale Mercado,
Alfredo Montana,
Stephen K. N. Portillo,
Jack Sayers,
Benjamin J. Vaughan,
Michael Zemcov,
Adi Zitrin
Abstract:
We present a measurement of the relativistic corrections to the thermal Sunyaev-Zel'dovich (SZ) effect spectrum, the rSZ effect, toward the massive galaxy cluster RX J1347.5-1145 by combining sub-mm images from Herschel-SPIRE with mm-wave Bolocam maps. Our analysis simultaneously models the SZ effect signal, the population of cosmic infrared background (CIB) galaxies, and galactic cirrus dust emis…
▽ More
We present a measurement of the relativistic corrections to the thermal Sunyaev-Zel'dovich (SZ) effect spectrum, the rSZ effect, toward the massive galaxy cluster RX J1347.5-1145 by combining sub-mm images from Herschel-SPIRE with mm-wave Bolocam maps. Our analysis simultaneously models the SZ effect signal, the population of cosmic infrared background (CIB) galaxies, and galactic cirrus dust emission in a manner that fully accounts for their spatial and frequency-dependent correlations. Gravitational lensing of background galaxies by RX J1347.5-1145 is included in our methodology based on a mass model derived from HST observations. Utilizing a set of realistic mock observations, we employ a forward modelling approach that accounts for the non-Gaussian covariances between observed astrophysical components to determine the posterior distribution of SZ effect brightness values consistent with the observed data. We determine a maximum a posteriori (MAP) value of the average Comptonization parameter of the intra-cluster medium (ICM) within R$_{2500}$ to be $\langle y \rangle_{2500} = 1.56 \times 10^{-4}$, with corresponding 68~per cent credible interval $[1.42,1.63] \times 10^{-4}$, and a MAP ICM electron temperature of $\langle \textrm{T}_{\textrm{sz}} \rangle_{2500} = 22.4$~keV with 68~per cent credible interval spanning $[10.4,33.0]$~keV. This is in good agreement with the pressure-weighted temperature obtained from {\it Chandra} X-ray observations, $\langle \textrm{T}_{\textrm{x,pw}}\rangle_{2500} = 17.4 \pm 2.3$~keV. We aim to apply this methodology to comparable existing data for a sample of 39 galaxy clusters, with an estimated uncertainty on the ensemble mean $\langle \textrm{T}_{\textrm{sz}} \rangle_{2500}$ at the $\simeq 1$~keV level, sufficiently precise to probe ICM physics and to inform X-ray temperature calibration.
△ Less
Submitted 26 October, 2021;
originally announced October 2021.
-
Sifting Through the Static: Moving Object Detection in Difference Images
Authors:
Hayden Smotherman,
Andrew J. Connolly,
J. Bryce Kalmbach,
Stephen K. N. Portillo,
Dino Bektesevic,
Siegfried Eggl,
Mario Juric,
Joachim Moeyens,
Peter J. Whidden
Abstract:
Trans-Neptunian Objects (TNOs) provide a window into the history of the Solar System, but they can be challenging to observe due to their distance from the Sun and relatively low brightness. Here we report the detection of 75 moving objects that we could not link to any other known objects, the faintest of which has a VR magnitude of $25.02 \pm 0.93$ using the KBMOD platform. We recover an additio…
▽ More
Trans-Neptunian Objects (TNOs) provide a window into the history of the Solar System, but they can be challenging to observe due to their distance from the Sun and relatively low brightness. Here we report the detection of 75 moving objects that we could not link to any other known objects, the faintest of which has a VR magnitude of $25.02 \pm 0.93$ using the KBMOD platform. We recover an additional 24 sources with previously-known orbits. We place constraints on the barycentric distance, inclination, and longitude of ascending node of these objects. The unidentified objects have a median barycentric distance of 41.28 au, placing them in the outer Solar System. The observed inclination and magnitude distribution of all detected objects is consistent with previously published KBO distributions. We describe extensions to KBMOD, including a robust percentile-based lightcurve filter, an in-line graphics processing unit (GPU) filter, new coadded stamp generation, and a convolutional neural network (CNN) stamp filter, which allow KBMOD to take advantage of difference images. These enchancements mark a significant improvement in the readiness of KBMOD for deployment on future big data surveys such as LSST.
△ Less
Submitted 7 September, 2021;
originally announced September 2021.
-
21cmVAE: A Very Accurate Emulator of the 21-cm Global Signal
Authors:
Christian H. Bye,
Stephen K. N. Portillo,
Anastasia Fialkov
Abstract:
Considerable observational efforts are being dedicated to measuring the sky-averaged (global) 21-cm signal of neutral hydrogen from Cosmic Dawn and the Epoch of Reionization. Deriving observational constraints on the astrophysics of this era requires modeling tools that can quickly and accurately generate theoretical signals across the wide astrophysical parameter space. For this purpose artificia…
▽ More
Considerable observational efforts are being dedicated to measuring the sky-averaged (global) 21-cm signal of neutral hydrogen from Cosmic Dawn and the Epoch of Reionization. Deriving observational constraints on the astrophysics of this era requires modeling tools that can quickly and accurately generate theoretical signals across the wide astrophysical parameter space. For this purpose artificial neural networks were used to create the only two existing global signal emulators, 21cmGEM and globalemu. In this paper we introduce 21cmVAE, a neural network-based global signal emulator, trained on the same dataset of ~30,000 global signals as the other two emulators, but with a more direct prediction algorithm that prioritizes accuracy and simplicity. Using neural networks, we compute derivatives of the signals with respect to the astrophysical parameters and establish the most important astrophysical processes that drive the global 21-cm signal at different epochs. 21cmVAE has a relative rms error of only 0.34 - equivalently 0.54 mK - on average, which is a significant improvement compared to the existing emulators, and a run time of 0.04 seconds per parameter set. The emulator, the code, and the processed datasets are publicly available at https://github.com/christianhbye/21cmVAE and through https://zenodo.org/record/5904939.
△ Less
Submitted 16 April, 2022; v1 submitted 12 July, 2021;
originally announced July 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.
-
The Catalogue for Astrophysical Turbulence Simulations (CATS)
Authors:
B. Burkhart,
S. Appel,
S. Bialy,
J. Cho,
A. J. Christensen,
D. Collins,
C. Federrath,
D. Fielding,
D. Finkbeiner,
A. S. Hill,
J. C. Ibanez-Mejia,
M. R. Krumholz,
A. Lazarian,
M. Li,
P. Mocz,
M. -M. Mac Low,
J. Naiman,
S. K. N. Portillo,
B. Shane,
Z. Slepian,
Y. Yuan
Abstract:
Turbulence is a key process in many fields of astrophysics. Advances in numerical simulations of fluids over the last several decades have revolutionized our understanding of turbulence and related processes such as star formation and cosmic ray propagation. However, data from numerical simulations of astrophysical turbulence are often not made public. We introduce a new simulation-oriented databa…
▽ More
Turbulence is a key process in many fields of astrophysics. Advances in numerical simulations of fluids over the last several decades have revolutionized our understanding of turbulence and related processes such as star formation and cosmic ray propagation. However, data from numerical simulations of astrophysical turbulence are often not made public. We introduce a new simulation-oriented database for the astronomical community: The Catalogue for Astrophysical Turbulence Simulations (CATS), located at www.mhdturbulence.com. CATS includes magnetohydrodynamic (MHD) turbulent box simulation data products generated by the public codes athena++, arepo, enzo, and flash. CATS also includes several synthetic observational data sets, such as turbulent HI data cubes. We also include measured power spectra and 3-point correlation functions from some of these data. We discuss the importance of open source statistical and visualization tools for the analysis of turbulence simulations such as those found in CATS.
△ Less
Submitted 21 October, 2020;
originally announced October 2020.
-
Dimensionality Reduction of SDSS Spectra with Variational Autoencoders
Authors:
Stephen K. N. Portillo,
John K. Parejko,
Jorge R. Vergara,
Andrew J. Connolly
Abstract:
High resolution galaxy spectra contain much information about galactic physics, but the high dimensionality of these spectra makes it difficult to fully utilize the information they contain. We apply variational autoencoders (VAEs), a non-linear dimensionality reduction technique, to a sample of spectra from the Sloan Digital Sky Survey. In contrast to Principal Component Analysis (PCA), a widely…
▽ More
High resolution galaxy spectra contain much information about galactic physics, but the high dimensionality of these spectra makes it difficult to fully utilize the information they contain. We apply variational autoencoders (VAEs), a non-linear dimensionality reduction technique, to a sample of spectra from the Sloan Digital Sky Survey. In contrast to Principal Component Analysis (PCA), a widely used technique, VAEs can capture non-linear relationships between latent parameters and the data. We find that a VAE can reconstruct the SDSS spectra well with only six latent parameters, outperforming PCA with the same number of components. Different galaxy classes are naturally separated in this latent space, without class labels having been given to the VAE. The VAE latent space is interpretable because the VAE can be used to make synthetic spectra at any point in latent space. For example, making synthetic spectra along tracks in latent space yields sequences of realistic spectra that interpolate between two different types of galaxies. Using the latent space to find outliers may yield interesting spectra: in our small sample, we immediately find unusual data artifacts and stars misclassified as galaxies. In this exploratory work, we show that VAEs create compact, interpretable latent spaces that capture non-linear features of the data. While a VAE takes substantial time to train (~1 day for 48000 spectra), once trained, VAEs can enable the fast exploration of large astronomical data sets.
△ Less
Submitted 9 July, 2020; v1 submitted 24 February, 2020;
originally announced February 2020.
-
Multiband Probabilistic Cataloging: A Joint Fitting Approach to Point Source Detection and Deblending
Authors:
Richard M. Feder,
Stephen K. N. Portillo,
Tansu Daylan,
Douglas Finkbeiner
Abstract:
Probabilistic cataloging (PCAT) outperforms traditional cataloging methods on single-band optical data in crowded fields (Portillo et al. 2017). We extend our work to multiple bands, achieving greater sensitivity ($\sim$ 0.4 mag) and greater speed (500x) compared to previous single-band results. We demonstrate the effectiveness of multiband PCAT on mock data, both in terms of recovering accurate p…
▽ More
Probabilistic cataloging (PCAT) outperforms traditional cataloging methods on single-band optical data in crowded fields (Portillo et al. 2017). We extend our work to multiple bands, achieving greater sensitivity ($\sim$ 0.4 mag) and greater speed (500x) compared to previous single-band results. We demonstrate the effectiveness of multiband PCAT on mock data, both in terms of recovering accurate posteriors in the catalog space, and in directly deblending sources. When applied to Sloan Digital Sky Survey (SDSS) observations of M2, taking Hubble Space Telescope data as truth, our joint fit on $r$ and $i$ band data goes $\sim0.4$ mag deeper than single-band probabilistic cataloging and has a false discovery rate less than 20\% for F606W$\leq 20$. Compared to DAOPHOT, the two-band SDSS catalog fit goes nearly 1.5 magnitudes deeper using the same data, and maintains a lower false discovery rate down to F606W$\sim 20.5$. Given recent improvements in computational speed, multiband PCAT shows promise in application to large-scale surveys and is a plausible framework for joint analysis of multi-instrument observational data.
△ Less
Submitted 20 March, 2020; v1 submitted 10 July, 2019;
originally announced July 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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
The Characterization of the Gamma-Ray Signal from the Central Milky Way: A Compelling Case for Annihilating Dark Matter
Authors:
Tansu Daylan,
Douglas P. Finkbeiner,
Dan Hooper,
Tim Linden,
Stephen K. N. Portillo,
Nicholas L. Rodd,
Tracy R. Slatyer
Abstract:
Past studies have identified a spatially extended excess of $\sim$1-3 GeV gamma rays from the region surrounding the Galactic Center, consistent with the emission expected from annihilating dark matter. We revisit and scrutinize this signal with the intention of further constraining its characteristics and origin. By applying cuts to the \textit{Fermi} event parameter CTBCORE, we suppress the tail…
▽ More
Past studies have identified a spatially extended excess of $\sim$1-3 GeV gamma rays from the region surrounding the Galactic Center, consistent with the emission expected from annihilating dark matter. We revisit and scrutinize this signal with the intention of further constraining its characteristics and origin. By applying cuts to the \textit{Fermi} event parameter CTBCORE, we suppress the tails of the point spread function and generate high resolution gamma-ray maps, enabling us to more easily separate the various gamma-ray components. Within these maps, we find the GeV excess to be robust and highly statistically significant, with a spectrum, angular distribution, and overall normalization that is in good agreement with that predicted by simple annihilating dark matter models. For example, the signal is very well fit by a 36-51 GeV dark matter particle annihilating to $b\bar{b}$ with an annihilation cross section of $σv = (1-3)\times 10^{-26}$ cm$^3$/s (normalized to a local dark matter density of 0.4 GeV/cm$^3$). Furthermore, we confirm that the angular distribution of the excess is approximately spherically symmetric and centered around the dynamical center of the Milky Way (within $\sim$$0.05^{\circ}$ of Sgr A$^*$), showing no sign of elongation along the Galactic Plane. The signal is observed to extend to at least $\simeq10^{\circ}$ from the Galactic Center, disfavoring the possibility that this emission originates from millisecond pulsars.
△ Less
Submitted 17 March, 2015; v1 submitted 26 February, 2014;
originally announced February 2014.