-
Exploring the Core-galaxy Connection
Authors:
Isabele Souza Vitório,
Michael Beuhlmann,
Eve Kovacs,
Patricia Larsen,
Nicholas Frontiere,
Katrin Heitmann
Abstract:
Halo core tracking is a novel concept designed to efficiently follow halo substructure in large simulations. We have recently developed this concept in gravity-only simulations to investigate the galaxy-halo connection in the context of empirical and semi-analytic models. Here, we incorporate information from hydrodynamics simulations, with an emphasis on establishing a connection between cores an…
▽ More
Halo core tracking is a novel concept designed to efficiently follow halo substructure in large simulations. We have recently developed this concept in gravity-only simulations to investigate the galaxy-halo connection in the context of empirical and semi-analytic models. Here, we incorporate information from hydrodynamics simulations, with an emphasis on establishing a connection between cores and galaxies. We compare cores across gravity-only, adiabatic hydrodynamics, and subgrid hydrodynamics simulations with the same initial phases. We demonstrate that cores are stable entities whose halo-centric radial profiles match across the simulations. We further develop a methodology that uses merging and infall mass cuts to group cores in the hydrodynamics simulation, creating on average, a one-to-one match to corresponding galaxies. We apply this methodology to cores from the gravity-only simulation, thus creating a proxy for galaxies which approximate the populations from the hydrodynamics simulation. Our results pave the way to incorporate inputs from smaller-scale hydrodynamics simulations directly into large-scale gravity-only runs in a principled manner.
△ Less
Submitted 28 June, 2024;
originally announced July 2024.
-
The New Worlds Simulations: Large-scale Simulations across Three Cosmologies
Authors:
Katrin Heitmann,
Thomas Uram,
Nicholas Frontiere,
Salman Habib,
Adrian Pope,
Silvio Rizzi,
Joe Insley
Abstract:
In this paper we describe the set of ``New Worlds Simulations'', three very large cosmology simulations, Qo'noS, Vulcan, and Ferenginar, that were carried out on the Summit supercomputer with the Hardware/Hybrid Cosmology Code, HACC. The gravity-only simulations follow the evolution of structure in the Universe by each employing 12,288^3 particles in (3 Gpc/h)^3 volumes, leading to a mass resoluti…
▽ More
In this paper we describe the set of ``New Worlds Simulations'', three very large cosmology simulations, Qo'noS, Vulcan, and Ferenginar, that were carried out on the Summit supercomputer with the Hardware/Hybrid Cosmology Code, HACC. The gravity-only simulations follow the evolution of structure in the Universe by each employing 12,288^3 particles in (3 Gpc/h)^3 volumes, leading to a mass resolution of m_p~10^9 Msun/h. The simulations cover three different cosmologies, one LambdaCDM model, consistent with measurements from Planck, one simulation with massive neutrinos, and one simulation with a varying dark energy equation of state. All simulations have the same phases to allow a detailed comparison of the results and the investigation of the impact of different cosmological parameters. We present measurements of some basic statistics, such as matter power spectra, correlation function, halo mass function and concentration-mass relation and investigate the differences due to the varying cosmologies. Given the large volume and high resolution, these simulations provide excellent bases for creating synthetic skies. A subset of the data is made publicly available as part of this paper.
△ Less
Submitted 11 June, 2024;
originally announced June 2024.
-
SPHEREx: NASA's Near-Infrared Spectrophotmetric All-Sky Survey
Authors:
Brendan P. Crill,
Michael Werner,
Rachel Akeson,
Matthew Ashby,
Lindsey Bleem,
James J. Bock,
Sean Bryan,
Jill Burnham,
Joyce Byunh,
Tzu-Ching Chang,
Yi-Kuan Chiang,
Walter Cook,
Asantha Cooray,
Andrew Davis,
Olivier Doré,
C. Darren Dowell,
Gregory Dubois-Felsmann,
Tim Eifler,
Andreas Faisst,
Salman Habib,
Chen Heinrich,
Katrin Heitmann,
Grigory Heaton,
Christopher Hirata,
Viktor Hristov
, et al. (29 additional authors not shown)
Abstract:
SPHEREx, the Spectro-Photometer for the History of the Universe, Epoch of Reionization, and ices Explorer, is a NASA MIDEX mission planned for launch in 2024. SPHEREx will carry out the first all-sky spectral survey at wavelengths between 0.75 micron and 5 micron with spectral resolving power ~40 between 0.75 and 3.8 micron and ~120 between 3.8 and 5 micron At the end of its two-year mission, SPHE…
▽ More
SPHEREx, the Spectro-Photometer for the History of the Universe, Epoch of Reionization, and ices Explorer, is a NASA MIDEX mission planned for launch in 2024. SPHEREx will carry out the first all-sky spectral survey at wavelengths between 0.75 micron and 5 micron with spectral resolving power ~40 between 0.75 and 3.8 micron and ~120 between 3.8 and 5 micron At the end of its two-year mission, SPHEREx will provide 0.75-to-5 micron spectra of each 6.2"x6.2" pixel on the sky - 14 billion spectra in all. This paper updates an earlier description of SPHEREx presenting changes made during the mission's Preliminary Design Phase, including a discussion of instrument integration and test and a summary of the data processing, analysis, and distribution plans.
△ Less
Submitted 16 April, 2024;
originally announced April 2024.
-
Optimization and Quality Assessment of Baryon Pasting for Intracluster Gas using the Borg Cube Simulation
Authors:
F. Kéruzoré,
L. E. Bleem,
M. Buehlmann,
J. D. Emberson,
N. Frontiere,
S. Habib,
K. Heitmann,
P. Larsen
Abstract:
Synthetic datasets generated from large-volume gravity-only simulations are an important tool in the calibration of cosmological analyses. Their creation often requires accurate inference of baryonic observables from the dark matter field. We explore the effectiveness of a baryon pasting algorithm in providing precise estimations of three-dimensional gas thermodynamic properties based on gravity-o…
▽ More
Synthetic datasets generated from large-volume gravity-only simulations are an important tool in the calibration of cosmological analyses. Their creation often requires accurate inference of baryonic observables from the dark matter field. We explore the effectiveness of a baryon pasting algorithm in providing precise estimations of three-dimensional gas thermodynamic properties based on gravity-only simulations. We use the Borg Cube, a pair of simulations originating from identical initial conditions, with one run evolved as a gravity-only simulation, and the other incorporating non-radiative hydrodynamics. Matching halos in both simulations enables comparisons of gas properties on an individual halo basis. This comparative analysis allows us to fit for the model parameters that yield the closest agreement between the gas properties in both runs. To capture the redshift evolution of these parameters, we perform the analysis at five distinct redshift steps, spanning from $z=0$ to $2$. We find that the investigated algorithm, utilizing information solely from the gravity-only simulation, achieves few-percent accuracy in reproducing the median intracluster gas pressure and density, albeit with a scatter of approximately 20%, for cluster-scale objects up to $z=2$. We measure the scaling relation between integrated Compton parameter and cluster mass ($Y_{500c} | M_{500c}$), and find that the imprecision of baryon pasting adds less than 5% to the intrinsic scatter measured in the hydrodynamic simulation. We provide best-fitting values and their redshift evolution, and discuss future investigations that will be undertaken to extend this work.
△ Less
Submitted 27 November, 2023; v1 submitted 23 June, 2023;
originally announced June 2023.
-
Simulating image coaddition with the Nancy Grace Roman Space Telescope: II. Analysis of the simulated images and implications for weak lensing
Authors:
Masaya Yamamoto,
Katherine Laliotis,
Emily Macbeth,
Tianqing Zhang,
Christopher M. Hirata,
M. A. Troxel,
Kaili Cao,
Ami Choi,
Jahmour Givans,
Katrin Heitmann,
Mustapha Ishak,
Mike Jarvis,
Eve Kovacs,
Heyang Long,
Rachel Mandelbaum,
Andy Park,
Anna Porredon,
Christopher W. Walter,
W. Michael Wood-Vasey
Abstract:
One challenge for applying current weak lensing analysis tools to the Nancy Grace Roman Space Telescope is that individual images will be undersampled. Our companion paper presented an initial application of Imcom - an algorithm that builds an optimal mapping from input to output pixels to reconstruct a fully sampled combined image - on the Roman image simulations. In this paper, we measure the ou…
▽ More
One challenge for applying current weak lensing analysis tools to the Nancy Grace Roman Space Telescope is that individual images will be undersampled. Our companion paper presented an initial application of Imcom - an algorithm that builds an optimal mapping from input to output pixels to reconstruct a fully sampled combined image - on the Roman image simulations. In this paper, we measure the output noise power spectra, identify the sources of the major features in the power spectra, and show that simple analytic models that ignore sampling effects underestimate the power spectra of the coadded noise images. We compute the moments of both idealized injected stars and fully simulated stars in the coadded images, and their 1- and 2-point statistics. We show that the idealized injected stars have root-mean-square ellipticity errors (1 - 6) x 10-4 per component depending on the band; the correlation functions are >= 2 orders of magnitude below requirements, indicating that the image combination step itself is using a small fraction of the overall Roman 2nd moment error budget, although the 4th moments are larger and warrant further investigation. The stars in the simulated sky images, which include blending and chromaticity effects, have correlation functions near the requirement level (and below the requirement level in a wide-band image constructed by stacking all 4 filters). We evaluate the noise-induced biases in the ellipticities of injected stars, and explain the resulting trends with an analytical model. We conclude by enumerating the next steps in developing an image coaddition pipeline for Roman.
△ Less
Submitted 12 January, 2024; v1 submitted 15 March, 2023;
originally announced March 2023.
-
Simulating image coaddition with the Nancy Grace Roman Space Telescope: I. Simulation methodology and general results
Authors:
Christopher M. Hirata,
Masaya Yamamoto,
Katherine Laliotis,
Emily Macbeth,
M. A. Troxel,
Tianqing Zhang,
Kaili Cao,
Ami Choi,
Jahmour Givans,
Katrin Heitmann,
Mustapha Ishak,
Mike Jarvis,
Eve Kovacs,
Heyang Long,
Rachel Mandelbaum,
Andy Park,
Anna Porredon,
Christopher W. Walter,
W. Michael Wood-Vasey
Abstract:
The upcoming Nancy Grace Roman Space Telescope will carry out a wide-area survey in the near infrared. A key science objective is the measurement of cosmic structure via weak gravitational lensing. Roman data will be undersampled, which introduces new challenges in the measurement of source galaxy shapes; a potential solution is to use linear algebra-based coaddition techniques such as Imcom that…
▽ More
The upcoming Nancy Grace Roman Space Telescope will carry out a wide-area survey in the near infrared. A key science objective is the measurement of cosmic structure via weak gravitational lensing. Roman data will be undersampled, which introduces new challenges in the measurement of source galaxy shapes; a potential solution is to use linear algebra-based coaddition techniques such as Imcom that combine multiple undersampled images to produce a single oversampled output mosaic with a desired "target" point spread function (PSF). We present here an initial application of Imcom to 0.64 square degrees of simulated Roman data, based on the Roman branch of the Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC) Data Challenge 2 (DC2) simulation. We show that Imcom runs successfully on simulated data that includes features such as plate scale distortions, chip gaps, detector defects, and cosmic ray masks. We simultaneously propagate grids of injected sources and simulated noise fields as well as the full simulation. We quantify the residual deviations of the PSF from the target (the "leakage"), as well as noise properties of the output images; we discuss how the overall tiling pattern as well as Moiré patterns appear in the final leakage and noise maps. We include appendices on interpolation algorithms and the interaction of undersampling with image processing operations that may be of broader applicability. The companion paper ("Paper II") explores the implications for weak lensing analyses.
△ Less
Submitted 12 January, 2024; v1 submitted 15 March, 2023;
originally announced March 2023.
-
Galaxy Clustering in the Mira-Titan Universe I: Emulators for the redshift space galaxy correlation function and galaxy-galaxy lensing
Authors:
Juliana Kwan,
Shun Saito,
Alexie Leauthaud,
Katrin Heitmann,
Salman Habib,
Nicholas Frontiere,
Hong Guo,
Song Huang,
Adrian Pope,
Sergio Rodríguez-Torres
Abstract:
We construct accurate emulators for the projected and redshift space galaxy correlation functions and excess surface density as measured by galaxy-galaxy lensing, based on Halo Occupation Distribution (HOD) modeling. Using the complete Mira-Titan suite of 111 $N$-body simulations, our emulators vary over eight cosmological parameters and include the effects of neutrino mass and dynamical dark ener…
▽ More
We construct accurate emulators for the projected and redshift space galaxy correlation functions and excess surface density as measured by galaxy-galaxy lensing, based on Halo Occupation Distribution (HOD) modeling. Using the complete Mira-Titan suite of 111 $N$-body simulations, our emulators vary over eight cosmological parameters and include the effects of neutrino mass and dynamical dark energy. We demonstrate that our emulators are sufficiently accurate for the analysis of the BOSS DR12 CMASS galaxy sample over the range 0.5 < r < 50 Mpc/h. Furthermore, we show that our emulators are capable of recovering unbiased cosmological constraints from realistic mock catalogs over the same range. Our mock catalog tests show the efficacy of combining small scale galaxy-galaxy lensing with redshift space clustering and that we can constrain the growth rate and σ_8 to 7% and 4.5% respectively for a CMASS-like sample using only the measurements covered by our emulator. With the inclusion of a CMB prior on H_0, this reduces to a 2% measurement on the growth rate.
△ Less
Submitted 6 June, 2023; v1 submitted 23 February, 2023;
originally announced February 2023.
-
Modeling the Galaxy Distribution in Clusters using Halo Cores
Authors:
Danila Korytov,
Esteban Rangel,
Lindsey Bleem,
Nicholas Frontiere,
Salman Habib,
Katrin Heitmann,
Joseph Hollowed,
Adrian Pope
Abstract:
The galaxy distribution in dark matter-dominated halos is expected to approximately trace the details of the underlying dark matter substructure. In this paper we introduce halo `core-tracking' as a way to efficiently follow the small-scale substructure in cosmological simulations and apply the technique to model the galaxy distribution in observed clusters. The method relies on explicitly trackin…
▽ More
The galaxy distribution in dark matter-dominated halos is expected to approximately trace the details of the underlying dark matter substructure. In this paper we introduce halo `core-tracking' as a way to efficiently follow the small-scale substructure in cosmological simulations and apply the technique to model the galaxy distribution in observed clusters. The method relies on explicitly tracking the set of particles identified as belonging to a halo's central density core, once a halo has attained a certain threshold mass. The halo cores are then followed throughout the entire evolution of the simulation. The aim of core-tracking is to simplify substructure analysis tasks by avoiding the use of subhalos and, at the same time, to more easily account for the so-called ``orphan'' galaxies, which have lost substantial dark mass due to tidal stripping. We show that simple models based on halo cores can reproduce the number and spatial distribution of galaxies found in optically-selected clusters in the Sloan Digital Sky Survey. We also discuss future applications of the core-tracking methodology in studying the galaxy-halo connection.
△ Less
Submitted 11 July, 2023; v1 submitted 8 February, 2023;
originally announced February 2023.
-
The catalog-to-cosmology framework for weak lensing and galaxy clustering for LSST
Authors:
J. Prat,
J. Zuntz,
Y. Omori,
C. Chang,
T. Tröster,
E. Pedersen,
C. García-García,
E. Phillips-Longley,
J. Sanchez,
D. Alonso,
X. Fang,
E. Gawiser,
K. Heitmann,
M. Ishak,
M. Jarvis,
E. Kovacs,
P. Larsen,
Y. -Y. Mao,
L. Medina Varela,
M. Paterno,
S. D. Vitenti,
Z. Zhang,
The LSST Dark Energy Science Collaboration
Abstract:
We present TXPipe, a modular, automated and reproducible pipeline for ingesting catalog data and performing all the calculations required to obtain quality-assured two-point measurements of lensing and clustering, and their covariances, with the metadata necessary for parameter estimation. The pipeline is developed within the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Sci…
▽ More
We present TXPipe, a modular, automated and reproducible pipeline for ingesting catalog data and performing all the calculations required to obtain quality-assured two-point measurements of lensing and clustering, and their covariances, with the metadata necessary for parameter estimation. The pipeline is developed within the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC), and designed for cosmology analyses using LSST data. In this paper, we present the pipeline for the so-called 3x2pt analysis -- a combination of three two-point functions that measure the auto- and cross-correlation between galaxy density and shapes. We perform the analysis both in real and harmonic space using TXPipe and other LSST-DESC tools. We validate the pipeline using Gaussian simulations and show that it accurately measures data vectors and recovers the input cosmology to the accuracy level required for the first year of LSST data under this simplified scenario. We also apply the pipeline to a realistic mock galaxy sample extracted from the CosmoDC2 simulation suite (Korytov et al. 2019). TXPipe establishes a baseline framework that can be built upon as the LSST survey proceeds. Furthermore, the pipeline is designed to be easily extended to science probes beyond the 3x2pt analysis.
△ Less
Submitted 21 April, 2023; v1 submitted 19 December, 2022;
originally announced December 2022.
-
The Future of High Energy Physics Software and Computing
Authors:
V. Daniel Elvira,
Steven Gottlieb,
Oliver Gutsche,
Benjamin Nachman,
S. Bailey,
W. Bhimji,
P. Boyle,
G. Cerati,
M. Carrasco Kind,
K. Cranmer,
G. Davies,
V. D. Elvira,
R. Gardner,
K. Heitmann,
M. Hildreth,
W. Hopkins,
T. Humble,
M. Lin,
P. Onyisi,
J. Qiang,
K. Pedro,
G. Perdue,
A. Roberts,
M. Savage,
P. Shanahan
, et al. (3 additional authors not shown)
Abstract:
Software and Computing (S&C) are essential to all High Energy Physics (HEP) experiments and many theoretical studies. The size and complexity of S&C are now commensurate with that of experimental instruments, playing a critical role in experimental design, data acquisition/instrumental control, reconstruction, and analysis. Furthermore, S&C often plays a leading role in driving the precision of th…
▽ More
Software and Computing (S&C) are essential to all High Energy Physics (HEP) experiments and many theoretical studies. The size and complexity of S&C are now commensurate with that of experimental instruments, playing a critical role in experimental design, data acquisition/instrumental control, reconstruction, and analysis. Furthermore, S&C often plays a leading role in driving the precision of theoretical calculations and simulations. Within this central role in HEP, S&C has been immensely successful over the last decade. This report looks forward to the next decade and beyond, in the context of the 2021 Particle Physics Community Planning Exercise ("Snowmass") organized by the Division of Particles and Fields (DPF) of the American Physical Society.
△ Less
Submitted 8 November, 2022; v1 submitted 11 October, 2022;
originally announced October 2022.
-
Report of the Topical Group on Dark Energy and Cosmic Acceleration: Complementarity of Probes and New Facilities for Snowmass 2021
Authors:
Brenna Flaugher,
Vivian Miranda,
David J. Schlegel,
Adam J. Anderson,
Felipe Andrade-Oliveira,
Eric J. Baxter,
Amy N. Bender,
Lindsey E. Bleem,
Chihway Chang,
Clarence C. Chang,
Thomas Y. Chen,
Kyle S. Dawson,
Seth W. Digel,
Alex Drlica-Wagner,
Simone Ferraro,
Alyssa Garcia,
Katrin Heitmann,
Alex G. Kim,
Eric V. Linder,
Sayan Mandal,
Rachel Mandelbaum,
Phil Marshall,
Joel Meyers,
Laura Newburgh,
Peter E. Nugent
, et al. (5 additional authors not shown)
Abstract:
The mechanism(s) driving the early- and late-time accelerated expansion of the Universe represent one of the most compelling mysteries in fundamental physics today. The path to understanding the causes of early- and late-time acceleration depends on fully leveraging ongoing surveys, developing and demonstrating new technologies, and constructing and operating new instruments. This report presents…
▽ More
The mechanism(s) driving the early- and late-time accelerated expansion of the Universe represent one of the most compelling mysteries in fundamental physics today. The path to understanding the causes of early- and late-time acceleration depends on fully leveraging ongoing surveys, developing and demonstrating new technologies, and constructing and operating new instruments. This report presents a multi-faceted vision for the cosmic survey program in the 2030s and beyond that derives from these considerations. Cosmic surveys address a wide range of fundamental physics questions, and are thus a unique and powerful component of the HEP experimental portfolio.
△ Less
Submitted 18 September, 2022;
originally announced September 2022.
-
Snowmass Computational Frontier: Topical Group Report on Experimental Algorithm Parallelization
Authors:
G. Cerati,
K. Heitmann,
W. Hopkins,
J. Bennett,
T. Y. Chen,
V. V. Gligorov,
O. Gutsche,
S. Habib,
M. Kortelainen,
C. Leggett,
R. Mandelbaum,
N. Whitehorn,
M. Williams
Abstract:
The substantial increase in data volume and complexity expected from future experiments will require significant investment to prepare experimental algorithms. These algorithms include physics object reconstruction, calibrations, and processing of observational data. In addition, the changing computing architecture landscape, which will be primarily composed of heterogeneous resources, will contin…
▽ More
The substantial increase in data volume and complexity expected from future experiments will require significant investment to prepare experimental algorithms. These algorithms include physics object reconstruction, calibrations, and processing of observational data. In addition, the changing computing architecture landscape, which will be primarily composed of heterogeneous resources, will continue to pose major challenges with regard to algorithmic migration. Portable tools need to be developed that can be shared among the frontiers (e.g., for code execution on different platforms) and opportunities, such as forums or cross-experimental working groups, need to be provided where experiences and lessons learned can be shared between experiments and frontiers. At the same time, individual experiments also need to invest considerable resources to develop algorithms unique to their needs (e.g., for facilities dedicated to the experiment), and ensure that their specific algorithms will be able to efficiently exploit external heterogeneous computing facilities. Common software tools represent a cost-effective solution, providing ready-to-use software solutions as well as a platform for R\&D work. These are particularly important for small experiments which typically do not have dedicated resources needed to face the challenges imposed by the evolving computing technologies. Workforce development is a key concern across frontiers and experiments, and additional support is needed to provide career opportunities for researchers working in the field of experimental algorithm development. Finally, cross-discipline collaborations going beyond high-energy physics are a key ingredient to address the challenges ahead and more support for such collaborations needs to be created. This report targets future experiments, observations and experimental algorithm development for the next 10-15 years.
△ Less
Submitted 15 September, 2022;
originally announced September 2022.
-
A Joint Roman Space Telescope and Rubin Observatory Synthetic Wide-Field Imaging Survey
Authors:
M. A. Troxel,
C. Lin,
A. Park,
C. Hirata,
R. Mandelbaum,
M. Jarvis,
A. Choi,
J. Givans,
M. Higgins,
B. Sanchez,
M. Yamamoto,
H. Awan,
J. Chiang,
O. Dore,
C. W. Walter,
T. Zhang,
J. Cohen-Tanugi,
E. Gawiser,
A. Hearin,
K. Heitmann,
M. Ishak,
E. Kovacs,
Y. -Y. Mao,
M. Wood-Vasey,
the LSST Dark Energy Science Collaboration
Abstract:
We present and validate 20 deg$^2$ of overlapping synthetic imaging surveys representing the full depth of the Nancy Grace Roman Space Telescope High-Latitude Imaging Survey (HLIS) and five years of observations of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). The two synthetic surveys are summarized, with reference to the existing 300 deg$^2$ of LSST simulated imaging prod…
▽ More
We present and validate 20 deg$^2$ of overlapping synthetic imaging surveys representing the full depth of the Nancy Grace Roman Space Telescope High-Latitude Imaging Survey (HLIS) and five years of observations of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). The two synthetic surveys are summarized, with reference to the existing 300 deg$^2$ of LSST simulated imaging produced as part of Dark Energy Science Collaboration (DESC) Data Challenge 2 (DC2). Both synthetic surveys observe the same simulated DESC DC2 universe. For the synthetic Roman survey, we simulate for the first time fully chromatic images along with the detailed physics of the Sensor Chip Assemblies derived from lab measurements using the flight detectors. The simulated imaging and resulting pixel-level measurements of photometric properties of objects span a wavelength range of $\sim$0.3 to 2.0 $μ$m. We also describe updates to the Roman simulation pipeline, changes in how astrophysical objects are simulated relative to the original DC2 simulations, and the resulting simulated Roman data products. We use these simulations to explore the relative fraction of unrecognized blends in LSST images, finding that 20-30% of objects identified in LSST images with $i$-band magnitudes brighter than 25 can be identified as multiple objects in Roman images. These simulations provide a unique testing ground for the development and validation of joint pixel-level analysis techniques of ground- and space-based imaging data sets in the second half of the 2020s -- in particular the case of joint Roman--LSST analyses.
△ Less
Submitted 14 September, 2022;
originally announced September 2022.
-
The Mira-Titan Universe IV. High Precision Power Spectrum Emulation
Authors:
Kelly R. Moran,
Katrin Heitmann,
Earl Lawrence,
Salman Habib,
Derek Bingham,
Amol Upadhye,
Juliana Kwan,
David Higdon,
Richard Payne
Abstract:
Modern cosmological surveys are delivering datasets characterized by unprecedented quality and statistical completeness; this trend is expected to continue into the future as new ground- and space-based surveys come online. In order to maximally extract cosmological information from these observations, matching theoretical predictions are needed. At low redshifts, the surveys probe the nonlinear r…
▽ More
Modern cosmological surveys are delivering datasets characterized by unprecedented quality and statistical completeness; this trend is expected to continue into the future as new ground- and space-based surveys come online. In order to maximally extract cosmological information from these observations, matching theoretical predictions are needed. At low redshifts, the surveys probe the nonlinear regime of structure formation where cosmological simulations are the primary means of obtaining the required information. The computational cost of sufficiently resolved large-volume simulations makes it prohibitive to run very large ensembles. Nevertheless, precision emulators built on a tractable number of high-quality simulations can be used to build very fast prediction schemes to enable a variety of cosmological inference studies. We have recently introduced the Mira-Titan Universe simulation suite designed to construct emulators for a range of cosmological probes. The suite covers the standard six cosmological parameters $\{ω_m,ω_b, σ_8, h, n_s, w_0\}$ and, in addition, includes massive neutrinos and a dynamical dark energy equation of state, $\{ω_ν, w_a\}$. In this paper we present the final emulator for the matter power spectrum based on 111 cosmological simulations, each covering a (2.1Gpc)$^3$ volume and evolving 3200$^3$ particles. An additional set of 1776 lower-resolution simulations and TimeRG perturbation theory results for the power spectrum are used to cover scales straddling the linear to mildly nonlinear regimes. The emulator provides predictions at the two to three percent level of accuracy over a wide range of cosmological parameters and is publicly released as part of this paper.
△ Less
Submitted 25 July, 2022;
originally announced July 2022.
-
Snowmass2021 Cosmic Frontier White Paper: Enabling Flagship Dark Energy Experiments to Reach their Full Potential
Authors:
Jonathan A. Blazek,
Doug Clowe,
Thomas E. Collett,
Ian P. Dell'Antonio,
Mark Dickinson,
Lluís Galbany,
Eric Gawiser,
Katrin Heitmann,
Renée Hložek,
Mustapha Ishak,
Saurabh W. Jha,
Alex G. Kim,
C. Danielle Leonard,
Anja von der Linden,
Michelle Lochner,
Rachel Mandelbaum,
Peter Melchior,
Joel Meyers,
Jeffrey A. Newman,
Peter Nugent,
Saul Perlmutter,
Daniel J. Perrefort,
Javier Sánchez,
Samuel J. Schmidt,
Sukhdeep Singh
, et al. (3 additional authors not shown)
Abstract:
A new generation of powerful dark energy experiments will open new vistas for cosmology in the next decade. However, these projects cannot reach their utmost potential without data from other telescopes. This white paper focuses in particular on the compelling benefits of ground-based spectroscopic and photometric observations to complement the Vera C. Rubin Observatory, as well as smaller program…
▽ More
A new generation of powerful dark energy experiments will open new vistas for cosmology in the next decade. However, these projects cannot reach their utmost potential without data from other telescopes. This white paper focuses in particular on the compelling benefits of ground-based spectroscopic and photometric observations to complement the Vera C. Rubin Observatory, as well as smaller programs in aid of a DESI-2 experiment and CMB-S4. These additional data sets will both improve dark energy constraints from these flagship projects beyond what would possible on their own and open completely new windows into fundamental physics. For example, additional photometry and single-object spectroscopy will provide necessary follow-up information for supernova and strong lensing cosmology, while highly-multiplexed spectroscopy both from smaller facilities over wide fields and from larger facilities over narrower regions of sky will yield more accurate photometric redshift estimates for weak lensing and galaxy clustering measurements from the Rubin Observatory, provide critical spectroscopic host galaxy redshifts for supernova Hubble diagrams, provide improved understanding of limiting astrophysical systematic effects, and enable new measurements that probe the nature of gravity. A common thread is that access to complementary data from a range of telescopes/instruments would have a substantial impact on the rate of advance of dark energy science in the coming years.
△ Less
Submitted 5 April, 2022;
originally announced April 2022.
-
Why are we still using 3D masses for cluster cosmology?
Authors:
Stijn N. B. Debackere,
Henk Hoekstra,
Joop Schaye,
Katrin Heitmann,
Salman Habib
Abstract:
The abundance of clusters of galaxies is highly sensitive to the late-time evolution of the matter distribution, since clusters form at the highest density peaks. However, the 3D cluster mass cannot be inferred without deprojecting the observations, introducing model-dependent biases and uncertainties due to the mismatch between the assumed and the true cluster density profile and the neglected ma…
▽ More
The abundance of clusters of galaxies is highly sensitive to the late-time evolution of the matter distribution, since clusters form at the highest density peaks. However, the 3D cluster mass cannot be inferred without deprojecting the observations, introducing model-dependent biases and uncertainties due to the mismatch between the assumed and the true cluster density profile and the neglected matter along the sightline. Since projected aperture masses can be measured directly in simulations and observationally through weak lensing, we argue that they are better suited for cluster cosmology. Using the Mira-Titan suite of gravity-only simulations, we show that aperture masses correlate strongly with 3D halo masses, albeit with large intrinsic scatter due to the varying matter distribution along the sightline. Nonetheless, aperture masses can be measured $\approx 2-3$ times more precisely from observations, since they do not require assumptions about the density profile and are only affected by the shape noise in the weak lensing measurements. We emulate the cosmology dependence of the aperture mass function directly with a Gaussian process. Comparing the cosmology sensitivity of the aperture mass function and the 3D halo mass function for a fixed survey solid angle and redshift interval, we find the aperture mass sensitivity is higher for $Ω_\mathrm{m}$ and $w_a$, similar for $σ_8$, $n_\mathrm{s}$, and $w_0$, and slightly lower for $h$. With a carefully calibrated aperture mass function emulator, cluster cosmology analyses can use cluster aperture masses directly, reducing the sensitivity to model-dependent mass calibration biases and uncertainties.
△ Less
Submitted 14 June, 2022; v1 submitted 30 March, 2022;
originally announced March 2022.
-
Cosmo-Paleontology: Statistics of Fossil Groups in a Gravity-Only Simulation
Authors:
Aurora Cossairt,
Michael Buehlmann,
Eve Kovacs,
Xin Liu,
Salman Habib,
Katrin Heitmann
Abstract:
We present a detailed study of fossil group candidates identified in "Last Journey", a gravity-only cosmological simulation covering a $(3.4\, h^{-1}\mathrm{Gpc})^3$ volume with a particle mass resolution of $m_p \approx 2.7 \times 10^9\, h^{-1}\mathrm{M}_\odot$. The simulation allows us to simultaneously capture a large number of group-scale halos and to resolve their internal structure. Historic…
▽ More
We present a detailed study of fossil group candidates identified in "Last Journey", a gravity-only cosmological simulation covering a $(3.4\, h^{-1}\mathrm{Gpc})^3$ volume with a particle mass resolution of $m_p \approx 2.7 \times 10^9\, h^{-1}\mathrm{M}_\odot$. The simulation allows us to simultaneously capture a large number of group-scale halos and to resolve their internal structure. Historically, fossil groups have been characterized by high X-ray brightness and a large luminosity gap between the brightest and second brightest galaxy in the group. In order to identify candidate halos that host fossil groups, we use halo merger tree information to introduce two parameters: a luminous merger mass threshold ($M_\mathrm{LM}$) and a last luminous merger redshift cut-off ($z_\mathrm{LLM}$). The final parameter choices are informed by observational data and allow us to identify a plausible fossil group sample from the simulation. The candidate halos are characterized by reduced substructure and are therefore less likely to host bright galaxies beyond the brightest central galaxy. We carry out detailed studies of this sample, including analysis of halo properties and clustering. We find that our simple assumptions lead to fossil group candidates that form early, have higher concentrations, and are more relaxed compared to other halos in the same mass range.
△ Less
Submitted 25 May, 2022; v1 submitted 16 March, 2022;
originally announced March 2022.
-
Snowmass 2021 CMB-S4 White Paper
Authors:
Kevork Abazajian,
Arwa Abdulghafour,
Graeme E. Addison,
Peter Adshead,
Zeeshan Ahmed,
Marco Ajello,
Daniel Akerib,
Steven W. Allen,
David Alonso,
Marcelo Alvarez,
Mustafa A. Amin,
Mandana Amiri,
Adam Anderson,
Behzad Ansarinejad,
Melanie Archipley,
Kam S. Arnold,
Matt Ashby,
Han Aung,
Carlo Baccigalupi,
Carina Baker,
Abhishek Bakshi,
Debbie Bard,
Denis Barkats,
Darcy Barron,
Peter S. Barry
, et al. (331 additional authors not shown)
Abstract:
This Snowmass 2021 White Paper describes the Cosmic Microwave Background Stage 4 project CMB-S4, which is designed to cross critical thresholds in our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. We provide an overview of the science case, the technical design, and project plan.
This Snowmass 2021 White Paper describes the Cosmic Microwave Background Stage 4 project CMB-S4, which is designed to cross critical thresholds in our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. We provide an overview of the science case, the technical design, and project plan.
△ Less
Submitted 15 March, 2022;
originally announced March 2022.
-
Snowmass2021 Computational Frontier White Paper: Cosmological Simulations and Modeling
Authors:
Marcelo A. Alvarez,
Arka Banerjee,
Simon Birrer,
Salman Habib,
Katrin Heitmann,
Zarija Lukić,
Julian B. Muñoz,
Yuuki Omori,
Hyunbae Park,
Annika H. G. Peter,
Jean Sexton,
Yi-Ming Zhong
Abstract:
Powerful new observational facilities will come online over the next decade, enabling a number of discovery opportunities in the "Cosmic Frontier", which targets understanding of the physics of the early universe, dark matter and dark energy, and cosmological probes of fundamental physics, such as neutrino masses and modifications of Einstein gravity. Synergies between different experiments will b…
▽ More
Powerful new observational facilities will come online over the next decade, enabling a number of discovery opportunities in the "Cosmic Frontier", which targets understanding of the physics of the early universe, dark matter and dark energy, and cosmological probes of fundamental physics, such as neutrino masses and modifications of Einstein gravity. Synergies between different experiments will be leveraged to present new classes of cosmic probes as well as to minimize systematic biases present in individual surveys. Success of this observational program requires actively pairing it with a well-matched state-of-the-art simulation and modeling effort. Next-generation cosmological modeling will increasingly focus on physically rich simulations able to model outputs of sky surveys spanning multiple wavebands. These simulations will have unprecedented resolution, volume coverage, and must deliver guaranteed high-fidelity results for individual surveys as well as for the cross-correlations across different surveys. The needed advances are as follows: (1) Development of scientifically rich and broadly-scoped simulations, which capture the relevant physics and correlations between probes (2) Accurate translation of simulation results into realistic image or spectral data to be directly compared with observations (3) Improved emulators and/or data-driven methods serving as surrogates for expensive simulations, constructed from a finite set of full-physics simulations (4) Detailed and transparent verification and validation programs for both simulations and analysis tools. (Abridged)
△ Less
Submitted 15 March, 2022; v1 submitted 14 March, 2022;
originally announced March 2022.
-
Snowmass2021 Cosmic Frontier White Paper: High Density Galaxy Clustering in the Regime of Cosmic Acceleration
Authors:
Kyle Dawson,
Andrew Hearin,
Katrin Heitmann,
Mustapha Ishak,
Johannes Ulf Lange,
Martin White,
Rongpu Zhou
Abstract:
Joint studies of imaging and spectroscopic samples, informed by theory and simulations, offer the potential for comprehensive tests of the cosmological model over redshifts z<1.5. Spectroscopic galaxy samples at these redshifts can be increased beyond the planned Dark Energy Spectroscopic Instrument (DESI) program by at least an order of magnitude, thus offering significantly more constraining pow…
▽ More
Joint studies of imaging and spectroscopic samples, informed by theory and simulations, offer the potential for comprehensive tests of the cosmological model over redshifts z<1.5. Spectroscopic galaxy samples at these redshifts can be increased beyond the planned Dark Energy Spectroscopic Instrument (DESI) program by at least an order of magnitude, thus offering significantly more constraining power for these joint studies. Spectroscopic observations of these galaxies in the latter half of the 2020's and beyond would leverage the theory and simulation effort in this regime. In turn, these high density observations will allow enhanced tests of dark energy, physics beyond the standard model, and neutrino masses that will greatly exceed what is currently possible. Here, we present a coordinated program of simulations, theoretical modeling, and future spectroscopy that would enable precise cosmological studies in the accelerating epoch where the effects of dark energy are most apparent.
△ Less
Submitted 14 March, 2022;
originally announced March 2022.
-
Snowmass2021 Cosmic Frontier White Paper: Rubin Observatory after LSST
Authors:
Bob Blum,
Seth W. Digel,
Alex Drlica-Wagner,
Salman Habib,
Katrin Heitmann,
Mustapha Ishak,
Saurabh W. Jha,
Steven M. Kahn,
Rachel Mandelbaum,
Phil Marshall,
Jeffrey A. Newman,
Aaron Roodman,
Christopher W. Stubbs
Abstract:
The Vera C. Rubin Observatory will begin the Legacy Survey of Space and Time (LSST) in 2024, spanning an area of 18,000 square degrees in six bands, with more than 800 observations of each field over ten years. The unprecedented data set will enable great advances in the study of the formation and evolution of structure and exploration of physics of the dark universe. The observations will hold cl…
▽ More
The Vera C. Rubin Observatory will begin the Legacy Survey of Space and Time (LSST) in 2024, spanning an area of 18,000 square degrees in six bands, with more than 800 observations of each field over ten years. The unprecedented data set will enable great advances in the study of the formation and evolution of structure and exploration of physics of the dark universe. The observations will hold clues about the cause for the accelerated expansion of the universe and possibly the nature of dark matter. During the next decade, LSST will be able to confirm or dispute if tensions seen today in cosmological data are due to new physics. New and unexpected phenomena could confirm or disrupt our current understanding of the universe. Findings from LSST will guide the path forward post-LSST. The Rubin Observatory will still be a uniquely powerful facility even then, capable of revealing further insights into the physics of the dark universe. These could be obtained via innovative observing strategies, e.g., targeting new probes at shorter timescales than with LSST, or via modest instrumental changes, e.g., new filters, or through an entirely new instrument for the focal plane. This White Paper highlights some of the opportunities in each scenario from Rubin observations after LSST.
△ Less
Submitted 14 March, 2022;
originally announced March 2022.
-
Simulating Hydrodynamics in Cosmology with CRK-HACC
Authors:
Nicholas Frontiere,
J. D. Emberson,
Michael Buehlmann,
Joseph Adamo,
Salman Habib,
Katrin Heitmann,
Claude-André Faucher-Giguère
Abstract:
We introduce CRK-HACC, an extension of the Hardware/Hybrid Accelerated Cosmology Code (HACC), to resolve gas hydrodynamics in large-scale structure formation simulations of the universe. The new framework couples the HACC gravitational N-body solver with a modern smoothed particle hydrodynamics (SPH) approach called CRKSPH. $\underline{\text{C}}$onservative $\underline{\text{R}}$eproducing…
▽ More
We introduce CRK-HACC, an extension of the Hardware/Hybrid Accelerated Cosmology Code (HACC), to resolve gas hydrodynamics in large-scale structure formation simulations of the universe. The new framework couples the HACC gravitational N-body solver with a modern smoothed particle hydrodynamics (SPH) approach called CRKSPH. $\underline{\text{C}}$onservative $\underline{\text{R}}$eproducing $\underline{\text{K}}$ernel $\underline{\text{SPH}}$ utilizes smoothing functions that exactly interpolate linear fields while manifestly preserving conservation laws (momentum, mass, and energy). The CRKSPH method has been incorporated to accurately model baryonic effects in cosmology simulations - an important addition targeting the generation of precise synthetic sky predictions for upcoming observational surveys. CRK-HACC inherits the codesign strategies of the HACC solver and is built to run on modern GPU-accelerated supercomputers. In this work, we summarize the primary solver components and present a number of standard validation tests to demonstrate code accuracy, including idealized hydrodynamic and cosmological setups, as well as self-similarity measurements.
△ Less
Submitted 22 November, 2022; v1 submitted 6 February, 2022;
originally announced February 2022.
-
Physical Benchmarking for AI-Generated Cosmic Web
Authors:
Xiaofeng Dong,
Nesar Ramachandra,
Salman Habib,
Katrin Heitmann,
Michael Buehlmann,
Sandeep Madireddy
Abstract:
The potential of deep learning based image-to-image translations has recently drawn a lot of attention; one intriguing possibility is that of generating cosmological predictions with a drastic reduction in computational cost. Such an effort requires optimization of neural networks with loss functions beyond low-order statistics like pixel-wise mean square error, and validation of results beyond si…
▽ More
The potential of deep learning based image-to-image translations has recently drawn a lot of attention; one intriguing possibility is that of generating cosmological predictions with a drastic reduction in computational cost. Such an effort requires optimization of neural networks with loss functions beyond low-order statistics like pixel-wise mean square error, and validation of results beyond simple visual comparisons and summary statistics. In order to study learning-based cosmological mappings, we choose a tractable analytical prescription - the Zel'dovich approximation - modeled using U-Net, a convolutional image translation framework. A comprehensive list of metrics is proposed, including higher-order correlation functions, conservation laws, topological indicators, dynamical robustness, and statistical independence of density fields. We find that the U-Net approach does well with some metrics but has difficulties with others. In addition to validating AI approaches using rigorous physical benchmarks, this study motivates advancements in domain-specific optimization schemes for scientific machine learning.
△ Less
Submitted 10 December, 2021;
originally announced December 2021.
-
Machine learning synthetic spectra for probabilistic redshift estimation: SYTH-Z
Authors:
Nesar Ramachandra,
Jonás Chaves-Montero,
Alex Alarcon,
Arindam Fadikar,
Salman Habib,
Katrin Heitmann
Abstract:
Photometric redshift estimation algorithms are often based on representative data from observational campaigns. Data-driven methods of this type are subject to a number of potential deficiencies, such as sample bias and incompleteness. Motivated by these considerations, we propose using physically motivated synthetic spectral energy distributions in redshift estimation. In addition, the synthetic…
▽ More
Photometric redshift estimation algorithms are often based on representative data from observational campaigns. Data-driven methods of this type are subject to a number of potential deficiencies, such as sample bias and incompleteness. Motivated by these considerations, we propose using physically motivated synthetic spectral energy distributions in redshift estimation. In addition, the synthetic data would have to span a domain in colour-redshift space concordant with that of the targeted observational surveys. With a matched distribution and realistically modelled synthetic data in hand, a suitable regression algorithm can be appropriately trained; we use a mixture density network for this purpose. We also perform a zero-point re-calibration to reduce the systematic differences between noise-free synthetic data and the (unavoidably) noisy observational data sets. This new redshift estimation framework, SYTH-Z, demonstrates superior accuracy over a wide range of redshifts compared to baseline models trained on observational data alone. Approaches using realistic synthetic data sets can therefore greatly mitigate the reliance on expensive spectroscopic follow-up for the next generation of photometric surveys.
△ Less
Submitted 23 November, 2021;
originally announced November 2021.
-
SNIa-Cosmology Analysis Results from Simulated LSST Images: from Difference Imaging to Constraints on Dark Energy
Authors:
B. Sánchez,
R. Kessler,
D. Scolnic,
B. Armstrong,
R. Biswas,
J. Bogart,
J. Chiang,
J. Cohen-Tanugi,
D. Fouchez,
Ph. Gris,
K. Heitmann,
R. Hložek,
S. Jha,
H. Kelly,
S. Liu,
G. Narayan,
B. Racine,
E. Rykoff,
M. Sullivan,
C. Walter,
M. Wood-Vasey,
The LSST Dark Energy Science Collaboration
Abstract:
The Vera Rubin Observatory Legacy Survey of Space and Time (LSST) is expected to process ${\sim}10^6$ transient detections per night. For precision measurements of cosmological parameters and rates, it is critical to understand the detection efficiency, magnitude limits, artifact contamination levels, and biases in the selection and photometry. Here we rigorously test the LSST Difference Image Ana…
▽ More
The Vera Rubin Observatory Legacy Survey of Space and Time (LSST) is expected to process ${\sim}10^6$ transient detections per night. For precision measurements of cosmological parameters and rates, it is critical to understand the detection efficiency, magnitude limits, artifact contamination levels, and biases in the selection and photometry. Here we rigorously test the LSST Difference Image Analysis (DIA) pipeline using simulated images from the Rubin Observatory LSST Dark Energy Science Collaboration (DESC) Data Challenge (DC2) simulation for the Wide-Fast-Deep (WFD) survey area. DC2 is the first large-scale (300 deg$^2$) image simulation of a transient survey that includes realistic cadence, variable observing conditions, and CCD image artifacts. We analyze ${\sim}$15 deg$^2$ of DC2 over a 5-year time-span in which artificial point-sources from Type Ia Supernovae (SNIa) light curves have been overlaid onto the images. We measure the detection efficiency as a function of Signal-to-Noise Ratio (SNR) and find a $50\%$ efficiency at $\rm{SNR}=5.8$. The magnitude limits for each filter are: $u=23.66$, $g=24.69$, $r=24.06$, $i=23.45$, $z=22.54$, $y=21.62$ $\rm{mag}$. The artifact contamination is $\sim90\%$ of detections, corresponding to $\sim1000$ artifacts/deg$^2$ in $g$ band, and falling to 300 per deg$^2$ in $y$ band. The photometry has biases $<1\%$ for magnitudes $19.5 < m <23$. Our DIA performance on simulated images is similar to that of the Dark Energy Survey pipeline applied to real images. We also characterize DC2 image properties to produce catalog-level simulations needed for distance bias corrections. We find good agreement between DC2 data and simulations for distributions of SNR, redshift, and fitted light-curve properties. Applying a realistic SNIa-cosmology analysis for redshifts $z<1$, we recover the input cosmology parameters to within statistical uncertainties.
△ Less
Submitted 12 November, 2021;
originally announced November 2021.
-
Validating Synthetic Galaxy Catalogs for Dark Energy Science in the LSST Era
Authors:
Eve Kovacs,
Yao-Yuan Mao,
Michel Aguena,
Anita Bahmanyar,
Adam Broussard,
James Butler,
Duncan Campbell,
Chihway Chang,
Shenming Fu,
Katrin Heitmann,
Danila Korytov,
François Lanusse,
Patricia Larsen,
Rachel Mandelbaum,
Christopher B. Morrison,
Constantin Payerne,
Marina Ricci,
Eli Rykoff,
F. Javier Sánchez,
Ignacio Sevilla-Noarbe,
Melanie Simet,
Chun-Hao To,
Vinu Vikraman,
Rongpu Zhou,
Camille Avestruz
, et al. (14 additional authors not shown)
Abstract:
Large simulation efforts are required to provide synthetic galaxy catalogs for ongoing and upcoming cosmology surveys. These extragalactic catalogs are being used for many diverse purposes covering a wide range of scientific topics. In order to be useful, they must offer realistically complex information about the galaxies they contain. Hence, it is critical to implement a rigorous validation proc…
▽ More
Large simulation efforts are required to provide synthetic galaxy catalogs for ongoing and upcoming cosmology surveys. These extragalactic catalogs are being used for many diverse purposes covering a wide range of scientific topics. In order to be useful, they must offer realistically complex information about the galaxies they contain. Hence, it is critical to implement a rigorous validation procedure that ensures that the simulated galaxy properties faithfully capture observations and delivers an assessment of the level of realism attained by the catalog. We present here a suite of validation tests that have been developed by the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC). We discuss how the inclusion of each test is driven by the scientific targets for static ground-based dark energy science and by the availability of suitable validation data. The validation criteria that are used to assess the performance of a catalog are flexible and depend on the science goals. We illustrate the utility of this suite by showing examples for the validation of cosmoDC2, the extragalactic catalog recently released for the LSST DESC second Data Challenge.
△ Less
Submitted 13 January, 2022; v1 submitted 7 October, 2021;
originally announced October 2021.
-
Extreme Scale Survey Simulation with Python Workflows
Authors:
A. S. Villarreal,
Yadu Babuji,
Tom Uram,
Daniel S. Katz,
Kyle Chard,
Katrin Heitmann
Abstract:
The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will soon carry out an unprecedented wide, fast, and deep survey of the sky in multiple optical bands. The data from LSST will open up a new discovery space in astronomy and cosmology, simultaneously providing clues toward addressing burning issues of the day, such as the origin of dark energy and and the nature of dark matter, w…
▽ More
The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will soon carry out an unprecedented wide, fast, and deep survey of the sky in multiple optical bands. The data from LSST will open up a new discovery space in astronomy and cosmology, simultaneously providing clues toward addressing burning issues of the day, such as the origin of dark energy and and the nature of dark matter, while at the same time yielding data that will, in turn, pose fresh new questions. To prepare for the imminent arrival of this remarkable data set, it is crucial that the associated scientific communities be able to develop the software needed to analyze it. Computational power now available allows us to generate synthetic data sets that can be used as a realistic training ground for such an effort. This effort raises its own challenges -- the need to generate very large simulations of the night sky, scaling up simulation campaigns to large numbers of compute nodes across multiple computing centers with different architectures, and optimizing the complex workload around memory requirements and widely varying wall clock times. We describe here a large-scale workflow that melds together Python code to steer the workflow, Parsl to manage the large-scale distributed execution of workflow components, and containers to carry out the image simulation campaign across multiple sites. Taking advantage of these tools, we developed an extreme-scale computational framework and used it to simulate five years of observations for 300 square degrees of sky area. We describe our experiences and lessons learned in developing this workflow capability, and highlight how the scalability and portability of our approach enabled us to efficiently execute it on up to 4000 compute nodes on two supercomputers.
△ Less
Submitted 24 September, 2021;
originally announced September 2021.
-
Farpoint: A High-Resolution Cosmology Simulation at the Gigaparsec Scale
Authors:
Nicholas Frontiere,
Katrin Heitmann,
Esteban Rangel,
Patricia Larsen,
Adrian Pope,
Imran Sultan,
Thomas Uram,
Salman Habib,
Silvio Rizzi,
Joe Insley
Abstract:
In this paper we introduce the Farpoint simulation, the latest member of the Hardware/Hybrid Accelerated Cosmology Code (HACC) gravity-only simulation family. The domain covers a volume of (1000$h^{-1}$Mpc)$^3$ and evolves close to two trillion particles, corresponding to a mass resolution of $m_p\sim 4.6\cdot 10^7 h^{-1}$M$_\odot$. These specifications enable comprehensive investigations of the g…
▽ More
In this paper we introduce the Farpoint simulation, the latest member of the Hardware/Hybrid Accelerated Cosmology Code (HACC) gravity-only simulation family. The domain covers a volume of (1000$h^{-1}$Mpc)$^3$ and evolves close to two trillion particles, corresponding to a mass resolution of $m_p\sim 4.6\cdot 10^7 h^{-1}$M$_\odot$. These specifications enable comprehensive investigations of the galaxy-halo connection, capturing halos down to small masses. Further, the large volume resolves scales typical of modern surveys with good statistical coverage of high mass halos. The simulation was carried out on the GPU-accelerated system Summit, one of the fastest supercomputers currently available. We provide specifics about the Farpoint run and present an initial set of results. The high mass resolution facilitates precise measurements of important global statistics, such as the halo concentration-mass relation and the correlation function down to small scales. Selected subsets of the simulation data products are publicly available via the HACC Simulation Data Portal.
△ Less
Submitted 28 February, 2022; v1 submitted 4 September, 2021;
originally announced September 2021.
-
The LSST-DESC 3x2pt Tomography Optimization Challenge
Authors:
Joe Zuntz,
François Lanusse,
Alex I. Malz,
Angus H. Wright,
Anže Slosar,
Bela Abolfathi,
David Alonso,
Abby Bault,
Clécio R. Bom,
Massimo Brescia,
Adam Broussard,
Jean-Eric Campagne,
Stefano Cavuoti,
Eduardo S. Cypriano,
Bernardo M. O. Fraga,
Eric Gawiser,
Elizabeth J. Gonzalez,
Dylan Green,
Peter Hatfield,
Kartheik Iyer,
David Kirkby,
Andrina Nicola,
Erfan Nourbakhsh,
Andy Park,
Gabriel Teixeira
, et al. (3 additional authors not shown)
Abstract:
This paper presents the results of the Rubin Observatory Dark Energy Science Collaboration (DESC) 3x2pt tomography challenge, which served as a first step toward optimizing the tomographic binning strategy for the main DESC analysis. The task of choosing an optimal tomographic binning scheme for a photometric survey is made particularly delicate in the context of a metacalibrated lensing catalogue…
▽ More
This paper presents the results of the Rubin Observatory Dark Energy Science Collaboration (DESC) 3x2pt tomography challenge, which served as a first step toward optimizing the tomographic binning strategy for the main DESC analysis. The task of choosing an optimal tomographic binning scheme for a photometric survey is made particularly delicate in the context of a metacalibrated lensing catalogue, as only the photometry from the bands included in the metacalibration process (usually riz and potentially g) can be used in sample definition.
The goal of the challenge was to collect and compare bin assignment strategies under various metrics of a standard 3x2pt cosmology analysis in a highly idealized setting to establish a baseline for realistically complex follow-up studies; in this preliminary study, we used two sets of cosmological simulations of galaxy redshifts and photometry under a simple noise model neglecting photometric outliers and variation in observing conditions, and contributed algorithms were provided with a representative and complete training set.
We review and evaluate the entries to the challenge, finding that even from this limited photometry information, multiple algorithms can separate tomographic bins reasonably well, reaching figures-of-merit scores close to the attainable maximum. We further find that adding the g band to riz photometry improves metric performance by ~15% and that the optimal bin assignment strategy depends strongly on the science case: which figure-of-merit is to be optimized, and which observables (clustering, lensing, or both) are included.
△ Less
Submitted 15 October, 2021; v1 submitted 30 August, 2021;
originally announced August 2021.
-
DESC DC2 Data Release Note
Authors:
LSST Dark Energy Science Collaboration,
Bela Abolfathi,
Robert Armstrong,
Humna Awan,
Yadu N. Babuji,
Franz Erik Bauer,
George Beckett,
Rahul Biswas,
Joanne R. Bogart,
Dominique Boutigny,
Kyle Chard,
James Chiang,
Johann Cohen-Tanugi,
Andrew J. Connolly,
Scott F. Daniel,
Seth W. Digel,
Alex Drlica-Wagner,
Richard Dubois,
Eric Gawiser,
Thomas Glanzman,
Salman Habib,
Andrew P. Hearin,
Katrin Heitmann,
Fabio Hernandez,
Renée Hložek
, et al. (32 additional authors not shown)
Abstract:
In preparation for cosmological analyses of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), the LSST Dark Energy Science Collaboration (LSST DESC) has created a 300 deg$^2$ simulated survey as part of an effort called Data Challenge 2 (DC2). The DC2 simulated sky survey, in six optical bands with observations following a reference LSST observing cadence, was processed with th…
▽ More
In preparation for cosmological analyses of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), the LSST Dark Energy Science Collaboration (LSST DESC) has created a 300 deg$^2$ simulated survey as part of an effort called Data Challenge 2 (DC2). The DC2 simulated sky survey, in six optical bands with observations following a reference LSST observing cadence, was processed with the LSST Science Pipelines (19.0.0). In this Note, we describe the public data release of the resulting object catalogs for the coadded images of five years of simulated observations along with associated truth catalogs. We include a brief description of the major features of the available data sets. To enable convenient access to the data products, we have developed a web portal connected to Globus data services. We describe how to access the data and provide example Jupyter Notebooks in Python to aid first interactions with the data. We welcome feedback and questions about the data release via a GitHub repository.
△ Less
Submitted 13 June, 2022; v1 submitted 12 January, 2021;
originally announced January 2021.
-
The Last Journey. II. SMACC -- Subhalo Mass-loss Analysis using Core Catalogs
Authors:
Imran Sultan,
Nicholas Frontiere,
Salman Habib,
Katrin Heitmann,
Eve Kovacs,
Patricia Larsen,
Esteban Rangel
Abstract:
This paper introduces SMACC -- Subhalo Mass-loss Analysis using Core Catalogs. SMACC adds a mass model to substructure merger trees based on halo "core tracking." Our approach avoids the need for running expensive subhalo finding algorithms and instead uses subhalo mass-loss modeling to assign masses to halo cores. We present details of the SMACC methodology and demonstrate its excellent performan…
▽ More
This paper introduces SMACC -- Subhalo Mass-loss Analysis using Core Catalogs. SMACC adds a mass model to substructure merger trees based on halo "core tracking." Our approach avoids the need for running expensive subhalo finding algorithms and instead uses subhalo mass-loss modeling to assign masses to halo cores. We present details of the SMACC methodology and demonstrate its excellent performance in describing halo substructure and its evolution. Validation of the approach is carried out using cosmological simulations at significantly different resolutions. We apply SMACC to the 1.24 trillion-particle Last Journey simulation and construct core catalogs with the additional mass information. These catalogs can be readily used as input to semi-analytic models or subhalo abundance matching approaches to determine approximate galaxy distributions, as well as for in-depth studies of small-scale structure evolution.
△ Less
Submitted 16 December, 2020;
originally announced December 2020.
-
The LSST DESC DC2 Simulated Sky Survey
Authors:
LSST Dark Energy Science Collaboration,
Bela Abolfathi,
David Alonso,
Robert Armstrong,
Éric Aubourg,
Humna Awan,
Yadu N. Babuji,
Franz Erik Bauer,
Rachel Bean,
George Beckett,
Rahul Biswas,
Joanne R. Bogart,
Dominique Boutigny,
Kyle Chard,
James Chiang,
Chuck F. Claver,
Johann Cohen-Tanugi,
Céline Combet,
Andrew J. Connolly,
Scott F. Daniel,
Seth W. Digel,
Alex Drlica-Wagner,
Richard Dubois,
Emmanuel Gangler,
Eric Gawiser
, et al. (55 additional authors not shown)
Abstract:
We describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses…
▽ More
We describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses this interconnectivity in a way that has not been attempted before. This effort encompasses a full end-to-end approach: starting from a large N-body simulation, through setting up LSST-like observations including realistic cadences, through image simulations, and finally processing with Rubin's LSST Science Pipelines. This last step ensures that we generate data products resembling those to be delivered by the Rubin Observatory as closely as is currently possible. The simulated DC2 sky survey covers six optical bands in a wide-fast-deep (WFD) area of approximately 300 deg^2 as well as a deep drilling field (DDF) of approximately 1 deg^2. We simulate 5 years of the planned 10-year survey. The DC2 sky survey has multiple purposes. First, the LSST DESC working groups can use the dataset to develop a range of DESC analysis pipelines to prepare for the advent of actual data. Second, it serves as a realistic testbed for the image processing software under development for LSST by the Rubin Observatory. In particular, simulated data provide a controlled way to investigate certain image-level systematic effects. Finally, the DC2 sky survey enables the exploration of new scientific ideas in both static and time-domain cosmology.
△ Less
Submitted 26 January, 2021; v1 submitted 12 October, 2020;
originally announced October 2020.
-
Matter Power Spectrum Emulator for f(R) Modified Gravity Cosmologies
Authors:
Nesar Ramachandra,
Georgios Valogiannis,
Mustapha Ishak,
Katrin Heitmann
Abstract:
Testing a subset of viable cosmological models beyond General Relativity (GR), with implications for cosmic acceleration and the Dark Energy associated with it, is within the reach of Rubin Observatory Legacy Survey of Space and Time (LSST) and a part of its endeavor. Deviations from GR-w(z)CDM models can manifest in the growth rate of structure and lensing, as well as in screening effects on non-…
▽ More
Testing a subset of viable cosmological models beyond General Relativity (GR), with implications for cosmic acceleration and the Dark Energy associated with it, is within the reach of Rubin Observatory Legacy Survey of Space and Time (LSST) and a part of its endeavor. Deviations from GR-w(z)CDM models can manifest in the growth rate of structure and lensing, as well as in screening effects on non-linear scales. We explore the constraining power of small-scale deviations predicted by the f(R) Hu-Sawicki Modified Gravity (MG) candidate, by emulating this model with COLA (COmoving Lagrangian Acceleration) simulations. We present the experimental design, data generation, and interpolation schemes in cosmological parameters and across redshifts for the emulation of the boost in the power spectra due to Modified Gravity effects. Three preliminary applications of the emulator highlight the sensitivity to cosmological parameters, Fisher forecasting and Markov Chain Monte Carlo inference for a fiducial cosmology. This emulator will play an important role for future cosmological analysis handling the formidable amount of data expected from Rubin Observatory LSST.
△ Less
Submitted 1 October, 2020;
originally announced October 2020.
-
CMB-S4: Forecasting Constraints on Primordial Gravitational Waves
Authors:
CMB-S4 Collaboration,
:,
Kevork Abazajian,
Graeme E. Addison,
Peter Adshead,
Zeeshan Ahmed,
Daniel Akerib,
Aamir Ali,
Steven W. Allen,
David Alonso,
Marcelo Alvarez,
Mustafa A. Amin,
Adam Anderson,
Kam S. Arnold,
Peter Ashton,
Carlo Baccigalupi,
Debbie Bard,
Denis Barkats,
Darcy Barron,
Peter S. Barry,
James G. Bartlett,
Ritoban Basu Thakur,
Nicholas Battaglia,
Rachel Bean,
Chris Bebek
, et al. (212 additional authors not shown)
Abstract:
CMB-S4---the next-generation ground-based cosmic microwave background (CMB) experiment---is set to significantly advance the sensitivity of CMB measurements and enhance our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. Among the science cases pursued with CMB-S4, the quest for detecting p…
▽ More
CMB-S4---the next-generation ground-based cosmic microwave background (CMB) experiment---is set to significantly advance the sensitivity of CMB measurements and enhance our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. Among the science cases pursued with CMB-S4, the quest for detecting primordial gravitational waves is a central driver of the experimental design. This work details the development of a forecasting framework that includes a power-spectrum-based semi-analytic projection tool, targeted explicitly towards optimizing constraints on the tensor-to-scalar ratio, $r$, in the presence of Galactic foregrounds and gravitational lensing of the CMB. This framework is unique in its direct use of information from the achieved performance of current Stage 2--3 CMB experiments to robustly forecast the science reach of upcoming CMB-polarization endeavors. The methodology allows for rapid iteration over experimental configurations and offers a flexible way to optimize the design of future experiments given a desired scientific goal. To form a closed-loop process, we couple this semi-analytic tool with map-based validation studies, which allow for the injection of additional complexity and verification of our forecasts with several independent analysis methods. We document multiple rounds of forecasts for CMB-S4 using this process and the resulting establishment of the current reference design of the primordial gravitational-wave component of the Stage-4 experiment, optimized to achieve our science goals of detecting primordial gravitational waves for $r > 0.003$ at greater than $5σ$, or, in the absence of a detection, of reaching an upper limit of $r < 0.001$ at $95\%$ CL.
△ Less
Submitted 27 August, 2020;
originally announced August 2020.
-
Building Halo Merger Trees from the Q Continuum Simulation
Authors:
Esteban Rangel,
Nicholas Frontiere,
Salman Habib,
Katrin Heitmann,
Wei-keng Liao,
Ankit Agrawal,
Alok Choudhary
Abstract:
Cosmological N-body simulations rank among the most computationally intensive efforts today. A key challenge is the analysis of structure, substructure, and the merger history for many billions of compact particle clusters, called halos. Effectively representing the merging history of halos is essential for many galaxy formation models used to generate synthetic sky catalogs, an important applicat…
▽ More
Cosmological N-body simulations rank among the most computationally intensive efforts today. A key challenge is the analysis of structure, substructure, and the merger history for many billions of compact particle clusters, called halos. Effectively representing the merging history of halos is essential for many galaxy formation models used to generate synthetic sky catalogs, an important application of modern cosmological simulations. Generating realistic mock catalogs requires computing the halo formation history from simulations with large volumes and billions of halos over many time steps, taking hundreds of terabytes of analysis data. We present fast parallel algorithms for producing halo merger trees and tracking halo substructure from a single-level, density-based clustering algorithm. Merger trees are created from analyzing the halo-particle membership function in adjacent snapshots, and substructure is identified by tracking the "cores" of merging halos -- sets of particles near the halo center. Core tracking is performed after creating merger trees and uses the relationships found during tree construction to associate substructures with hosts. The algorithms are implemented with MPI and evaluated on a Cray XK7 supercomputer using up to 16,384 processes on data from HACC, a modern cosmological simulation framework. We present results for creating merger trees from 101 analysis snapshots taken from the Q Continuum, a large volume, high mass resolution, cosmological simulation evolving half a trillion particles.
△ Less
Submitted 19 August, 2020;
originally announced August 2020.
-
The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: exploring the Halo Occupation Distribution model for Emission Line Galaxies
Authors:
Santiago Avila,
Violeta Gonzalez-Perez,
Faizan G. Mohammad,
Arnaud de Mattia,
Cheng Zhao,
Anand Raichoor,
Amelie Tamone,
Shadab Alam,
Julian Bautista,
Davide Bianchi,
Etienne Burtin,
Michael J. Chapman,
Chia-Hsun Chuang,
Johan Comparat,
Kyle Dawson,
Thomas Divers,
Helion du Mas des Bourboux,
Hector Gil-Marin,
Eva-Maria Mueller,
Salman Habib,
Katrin Heitmann,
Vanina Ruhlmann-Kleider,
Nelson Padilla,
Will J. Percival,
Ashley J. Ross
, et al. (3 additional authors not shown)
Abstract:
We study the modelling of the Halo Occupation Distribution (HOD) for the eBOSS DR16 Emission Line Galaxies (ELGs). Motivated by previous theoretical and observational studies, we consider different physical effects that can change how ELGs populate haloes. We explore the shape of the average HOD, the fraction of satellite galaxies, their probability distribution function (PDF), and their density a…
▽ More
We study the modelling of the Halo Occupation Distribution (HOD) for the eBOSS DR16 Emission Line Galaxies (ELGs). Motivated by previous theoretical and observational studies, we consider different physical effects that can change how ELGs populate haloes. We explore the shape of the average HOD, the fraction of satellite galaxies, their probability distribution function (PDF), and their density and velocity profiles. Our baseline HOD shape was fitted to a semi-analytical model of galaxy formation and evolution, with a decaying occupation of central ELGs at high halo masses. We consider Poisson and sub/super-Poissonian PDFs for satellite assignment. We model both NFW and particle profiles for satellite positions, also allowing for decreased concentrations. We model velocities with the virial theorem and particle velocity distributions. Additionally, we introduce a velocity bias and a net infall velocity. We study how these choices impact the clustering statistics while keeping the number density and bias fixed to that from eBOSS ELGs. The projected correlation function, $w_p$, captures most of the effects from the PDF and satellites profile. The quadrupole, $ξ_2$, captures most of the effects coming from the velocity profile. We find that the impact of the mean HOD shape is subdominant relative to the rest of choices. We fit the clustering of the eBOSS DR16 ELG data under different combinations of the above assumptions. The catalogues presented here have been analysed in companion papers, showing that eBOSS RSD+BAO measurements are insensitive to the details of galaxy physics considered here. These catalogues are made publicly available.
△ Less
Submitted 27 November, 2020; v1 submitted 17 July, 2020;
originally announced July 2020.
-
The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: N-body Mock Challenge for the eBOSS Emission Line Galaxy Sample
Authors:
Shadab Alam,
Arnaud de Mattia,
Amélie Tamone,
S. Ávila,
John A. Peacock,
V. Gonzalez-Perez,
Alex Smith,
Anand Raichoor,
Ashley J. Ross,
Julian E. Bautista,
Etienne Burtin,
Johan Comparat,
Kyle S. Dawson,
Hélion du Mas des Bourboux,
Stéphanie Escoffier,
Héctor Gil-Marín,
Salman Habib,
Katrin Heitmann,
Jiamin Hou,
Faizan G. Mohammad,
Eva-Maria Mueller,
Richard Neveux,
Romain Paviot,
Will J. Percival,
Graziano Rossi
, et al. (5 additional authors not shown)
Abstract:
Cosmological growth can be measured in the redshift space clustering of galaxies targeted by spectroscopic surveys. Accurate prediction of clustering of galaxies will require understanding galaxy physics which is a very hard and highly non-linear problem. Approximate models of redshift space distortion (RSD) take a perturbative approach to solve the evolution of dark matter and galaxies in the uni…
▽ More
Cosmological growth can be measured in the redshift space clustering of galaxies targeted by spectroscopic surveys. Accurate prediction of clustering of galaxies will require understanding galaxy physics which is a very hard and highly non-linear problem. Approximate models of redshift space distortion (RSD) take a perturbative approach to solve the evolution of dark matter and galaxies in the universe.
In this paper we focus on eBOSS emission line galaxies (ELGs) which live in intermediate mass haloes. We create a series of mock catalogues using haloes from the Multidark and {\sc Outer Rim} dark matter only N-body simulations. Our mock catalogues include various effects inspired by baryonic physics such as assembly bias and the characteristics of satellite galaxies kinematics, dynamics and statistics deviating from dark matter particles.
We analyse these mocks using the TNS RSD model in Fourier space and the CLPT in configuration space. We conclude that these two RSD models provide an unbiased measurement of redshift space distortion within the statistical error of our mocks. We obtain the conservative theoretical systematic uncertainty of $3.3\%$, $1.8\%$ and $1.5\%$ in $fσ_8$, $α_{\parallel}$ and $α_{\bot}$ respectively for the TNS and CLPT models. We note that the estimated theoretical systematic error is an order of magnitude smaller than the statistical error of the eBOSS ELG sample and hence are negligible for the purpose of the current eBOSS ELG analysis.
△ Less
Submitted 24 September, 2021; v1 submitted 17 July, 2020;
originally announced July 2020.
-
The Completed SDSS-IV Extended Baryon Oscillation Spectroscopic Survey: N-body Mock Challenge for the Quasar Sample
Authors:
Alex Smith,
Etienne Burtin,
Jiamin Hou,
Richard Neveux,
Ashley J. Ross,
Shadab Alam,
Jonathan Brinkmann,
Kyle S. Dawson,
Salman Habib,
Katrin Heitmann,
Jean-Paul Kneib,
Brad W. Lyke,
Hélion du Mas des Bourboux,
Eva-Maria Mueller,
Adam D. Myers,
Will J. Percival,
Graziano Rossi,
Donald P. Schneider,
Pauline Zarrouk,
Gong-Bo Zhao
Abstract:
The growth rate and expansion history of the Universe can be measured from large galaxy redshift surveys using the Alcock-Paczynski effect. We validate the Redshift Space Distortion models used in the final analysis of the Sloan Digital Sky Survey (SDSS) extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 quasar clustering sample, in configuration and Fourier space, using a se…
▽ More
The growth rate and expansion history of the Universe can be measured from large galaxy redshift surveys using the Alcock-Paczynski effect. We validate the Redshift Space Distortion models used in the final analysis of the Sloan Digital Sky Survey (SDSS) extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 quasar clustering sample, in configuration and Fourier space, using a series of HOD mock catalogues generated using the OuterRim N-body simulation. We test three models on a series of non-blind mocks, in the OuterRim cosmology, and blind mocks, which have been rescaled to new cosmologies, and investigate the effects of redshift smearing and catastrophic redshifts. We find that for the non-blind mocks, the models are able to recover $fσ_8$ to within 3% and $α_\parallel$ and $α_\bot$ to within 1%. The scatter in the measurements is larger for the blind mocks, due to the assumption of an incorrect fiducial cosmology. From this mock challenge, we find that all three models perform well, with similar systematic errors on $fσ_8$, $α_\parallel$ and $α_\bot$ at the level of $σ_{fσ_8}=0.013$, $σ_{α_\parallel}=0.012$ and $σ_{α_\bot}=0.008$. The systematic error on the combined consensus is $σ_{fσ_8}=0.011$, $σ_{α_\parallel}=0.008$ and $σ_{α_\bot}=0.005$, which is used in the final DR16 analysis. For BAO fits in configuration and Fourier space, we take conservative systematic errors of $σ_{α_\parallel}=0.010$ and $σ_{α_\bot}=0.007$.
△ Less
Submitted 4 December, 2020; v1 submitted 17 July, 2020;
originally announced July 2020.
-
The Completed SDSS-IV Extended Baryon Oscillation Spectroscopic Survey: N-body Mock Challenge for Galaxy Clustering Measurements
Authors:
Graziano Rossi,
Peter D. Choi,
Jeongin Moon,
Julian E. Bautista,
Hector Gil-Marin,
Romain Paviot,
Mariana Vargas-Magana,
Sylvain de la Torre,
Sebastien Fromenteau,
Ashley J. Ross,
Santiago Avila,
Etienne Burtin,
Kyle S. Dawson,
Stephanie Escoffier,
Salman Habib,
Katrin Heitmann,
Jiamin Hou,
Eva-Maria Mueller,
Will J. Percival,
Alex Smith,
Cheng Zhao,
Gong-Bo Zhao
Abstract:
We develop a series of N-body data challenges, functional to the final analysis of the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 (DR16) galaxy sample. The challenges are primarily based on high-fidelity catalogs constructed from the Outer Rim simulation - a large box size realization (3 Gpc/h) characterized by an unprecedented combination of volume and mass resolutio…
▽ More
We develop a series of N-body data challenges, functional to the final analysis of the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 (DR16) galaxy sample. The challenges are primarily based on high-fidelity catalogs constructed from the Outer Rim simulation - a large box size realization (3 Gpc/h) characterized by an unprecedented combination of volume and mass resolution, down to 1.85x10^9 M_sun/h. We generate synthetic galaxy mocks by populating Outer Rim halos with a variety of halo occupation distribution (HOD) schemes of increasing complexity, spanning different redshift intervals. We then assess the performance of three complementary redshift space distortion (RSD) models in configuration and Fourier space, adopted for the analysis of the complete DR16 eBOSS sample of Luminous Red Galaxies (LRGs). We find all the methods mutually consistent, with comparable systematic errors on the Alcock-Paczynski parameters and the growth of structure, and robust to different HOD prescriptions - thus validating the robustness of the models and the pipelines used for the baryon acoustic oscillation (BAO) and full shape clustering analysis. In particular, all the techniques are able to recover a_par and a_perp to within 0.9%, and fsig8 to within 1.5%. As a by-product of our work, we are also able to gain interesting insights on the galaxy-halo connection. Our study is relevant for the final eBOSS DR16 `consensus cosmology', as the systematic error budget is informed by testing the results of analyses against these high-resolution mocks. In addition, it is also useful for future large-volume surveys, since similar mock-making techniques and systematic corrections can be readily extended to model for instance the Dark Energy Spectroscopic Instrument (DESI) galaxy sample.
△ Less
Submitted 25 March, 2021; v1 submitted 17 July, 2020;
originally announced July 2020.
-
The Last Journey. I. An Extreme-Scale Simulation on the Mira Supercomputer
Authors:
Katrin Heitmann,
Nicholas Frontiere,
Esteban Rangel,
Patricia Larsen,
Adrian Pope,
Imran Sultan,
Thomas Uram,
Salman Habib,
Hal Finkel,
Danila Korytov,
Eve Kovacs,
Silvio Rizzi,
Joe Insley
Abstract:
The Last Journey is a large-volume, gravity-only, cosmological N-body simulation evolving more than 1.24 trillion particles in a periodic box with a side-length of 5.025Gpc. It was implemented using the HACC simulation and analysis framework on the BG/Q system, Mira. The cosmological parameters are chosen to be consistent with the results from the Planck satellite. A range of analysis tools have b…
▽ More
The Last Journey is a large-volume, gravity-only, cosmological N-body simulation evolving more than 1.24 trillion particles in a periodic box with a side-length of 5.025Gpc. It was implemented using the HACC simulation and analysis framework on the BG/Q system, Mira. The cosmological parameters are chosen to be consistent with the results from the Planck satellite. A range of analysis tools have been run in situ to enable a diverse set of science projects, and at the same time, to keep the resulting data amount manageable. Analysis outputs have been generated starting at redshift z~10 to allow for construction of synthetic galaxy catalogs using a semi-analytic modeling approach in post-processing. As part of our in situ analysis pipeline we employ a new method for tracking halo sub-structures, introducing the concept of subhalo cores. The production of multi-wavelength synthetic sky maps is facilitated by generating particle lightcones in situ, also beginning at z~10. We provide an overview of the simulation set-up and the generated data products; a first set of analysis results is presented. A subset of the data is publicly available.
△ Less
Submitted 8 January, 2021; v1 submitted 2 June, 2020;
originally announced June 2020.
-
Report from the Tri-Agency Cosmological Simulation Task Force
Authors:
Nick Battaglia,
Andrew Benson,
Tim Eifler,
Andrew Hearin,
Katrin Heitmann,
Shirley Ho,
Alina Kiessling,
Zarija Lukic,
Michael Schneider,
Elena Sellentin,
Joachim Stadel
Abstract:
The Tri-Agency Cosmological Simulations (TACS) Task Force was formed when Program Managers from the Department of Energy (DOE), the National Aeronautics and Space Administration (NASA), and the National Science Foundation (NSF) expressed an interest in receiving input into the cosmological simulations landscape related to the upcoming DOE/NSF Vera Rubin Observatory (Rubin), NASA/ESA's Euclid, and…
▽ More
The Tri-Agency Cosmological Simulations (TACS) Task Force was formed when Program Managers from the Department of Energy (DOE), the National Aeronautics and Space Administration (NASA), and the National Science Foundation (NSF) expressed an interest in receiving input into the cosmological simulations landscape related to the upcoming DOE/NSF Vera Rubin Observatory (Rubin), NASA/ESA's Euclid, and NASA's Wide Field Infrared Survey Telescope (WFIRST). The Co-Chairs of TACS, Katrin Heitmann and Alina Kiessling, invited community scientists from the USA and Europe who are each subject matter experts and are also members of one or more of the surveys to contribute. The following report represents the input from TACS that was delivered to the Agencies in December 2018.
△ Less
Submitted 14 May, 2020;
originally announced May 2020.
-
The Mira-Titan Universe. III. Emulation of the Halo Mass Function
Authors:
Sebastian Bocquet,
Katrin Heitmann,
Salman Habib,
Earl Lawrence,
Thomas Uram,
Nicholas Frontiere,
Adrian Pope,
Hal Finkel
Abstract:
We construct an emulator for the halo mass function over group and cluster mass scales for a range of cosmologies, including the effects of dynamical dark energy and massive neutrinos. The emulator is based on the recently completed Mira-Titan Universe suite of cosmological $N$-body simulations. The main set of simulations spans 111 cosmological models with 2.1 Gpc boxes. We extract halo catalogs…
▽ More
We construct an emulator for the halo mass function over group and cluster mass scales for a range of cosmologies, including the effects of dynamical dark energy and massive neutrinos. The emulator is based on the recently completed Mira-Titan Universe suite of cosmological $N$-body simulations. The main set of simulations spans 111 cosmological models with 2.1 Gpc boxes. We extract halo catalogs in the redshift range $z=[0.0, 2.0]$ and for masses $M_{200\mathrm{c}}\geq 10^{13}M_\odot/h$. The emulator covers an 8-dimensional hypercube spanned by {$Ω_\mathrm{m}h^2$, $Ω_\mathrm{b}h^2$, $Ω_νh^2$, $σ_8$, $h$, $n_s$, $w_0$, $w_a$}; spatial flatness is assumed. We obtain smooth halo mass functions by fitting piecewise second-order polynomials to the halo catalogs and employ Gaussian process regression to construct the emulator while keeping track of the statistical noise in the input halo catalogs and uncertainties in the regression process. For redshifts $z\lesssim1$, the typical emulator precision is better than $2\%$ for $10^{13}-10^{14} M_\odot/h$ and $<10\%$ for $M\simeq 10^{15}M_\odot/h$. For comparison, fitting functions using the traditional universal form for the halo mass function can be biased at up to 30\% at $M\simeq 10^{14}M_\odot/h$ for $z=0$. Our emulator is publicly available at \url{https://github.com/SebastianBocquet/MiraTitanHMFemulator}.
△ Less
Submitted 5 August, 2020; v1 submitted 26 March, 2020;
originally announced March 2020.
-
On the possibility of Baryon Acoustic Oscillation measurements at redshift $z>7.6$ with the Roman Space Telescope
Authors:
Siddharth Satpathy,
Zhaozhou An,
Rupert A. C. Croft,
Tiziana Di Matteo,
Ananth Tenneti,
Yu Feng,
Katrin Heitmann,
Graziano Rossi
Abstract:
The Nancy Grace Roman Space Telescope (RST), with its field of view and high sensitivity will make surveys of cosmological large-scale structure possible at high redshifts. We investigate the possibility of detecting Baryon Acoustic Oscillations (BAO) at redshifts $z>7.6$ for use as a standard ruler. We use data from the hydrodynamic simulation \textsc{BlueTides} in conjunction with the gigaparsec…
▽ More
The Nancy Grace Roman Space Telescope (RST), with its field of view and high sensitivity will make surveys of cosmological large-scale structure possible at high redshifts. We investigate the possibility of detecting Baryon Acoustic Oscillations (BAO) at redshifts $z>7.6$ for use as a standard ruler. We use data from the hydrodynamic simulation \textsc{BlueTides} in conjunction with the gigaparsec-scale Outer Rim simulation and a model for patchy reionization to create mock RST High Latitude Survey grism data for Lyman-alpha emission line selected galaxies at redshifts $z=7.4$ to $z=10$, covering 2280 square degrees. We measure the monopoles of galaxies in the mock catalogues and fit the BAO features. We find that for a line flux of $L = 7\times 10^{-17} \ {\rm erg/s/cm}^{2}$, the $5 σ$ detection limit for the current design, the BAO feature is partially detectable (measured in three out of four survey quadrants analysed independently). The resulting root mean square error on the angular diameter distance to $z=7.7$ is 7.9$\%$. If we improve the detection sensitivity by a factor of two (i.e. $L = 3.5\times 10^{-17} \ {\rm erg/s/cm}^{2}$), the distance error reduces to $1.4\%$. We caution that many more factors are yet to be modelled, including dust obscuration, the damping wing due to the intergalactic medium, and low redshift interlopers. If these issues do not strongly affect the results, or different observational techniques (such as use of multiple lines) can mitigate them, RST or similar instruments may be able to constrain the angular diameter distance to the high redshift Universe.
△ Less
Submitted 20 October, 2020; v1 submitted 15 January, 2020;
originally announced January 2020.
-
A Modular Deep Learning Pipeline for Galaxy-Scale Strong Gravitational Lens Detection and Modeling
Authors:
Sandeep Madireddy,
Nesar Ramachandra,
Nan Li,
James Butler,
Prasanna Balaprakash,
Salman Habib,
Katrin Heitmann,
The LSST Dark Energy Science Collaboration
Abstract:
Upcoming large astronomical surveys are expected to capture an unprecedented number of strong gravitational lensing systems. Deep learning is emerging as a promising practical tool for the detection and quantification of these galaxy-scale image distortions. The absence of large quantities of representative data from current astronomical surveys motivates the development of a robust forward-modeli…
▽ More
Upcoming large astronomical surveys are expected to capture an unprecedented number of strong gravitational lensing systems. Deep learning is emerging as a promising practical tool for the detection and quantification of these galaxy-scale image distortions. The absence of large quantities of representative data from current astronomical surveys motivates the development of a robust forward-modeling approach using synthetic lensing images. Using a mock sample of strong lenses created upon a state-of-the-art extragalactic catalogs, we train a modular deep learning pipeline for uncertainty-quantified detection and modeling with intermediate image processing components for denoising and deblending the lensing systems. We demonstrate a high degree of interpretability and controlled systematics due to domain-specific task modules trained with different stages of synthetic image generation. For lens detection and modeling, we obtain semantically meaningful latent spaces that separate classes of strong lens images and yield uncertainty estimates that explain the origin of misclassified images and provide probabilistic predictions for the lens parameters. Validation of the inference pipeline has been carried out using images from the Subaru telescope's Hyper Suprime-Cam camera, and LSST DESC simulated DC2 sky survey catalogues.
△ Less
Submitted 21 October, 2022; v1 submitted 10 November, 2019;
originally announced November 2019.
-
CMB-S4 Decadal Survey APC White Paper
Authors:
Kevork Abazajian,
Graeme Addison,
Peter Adshead,
Zeeshan Ahmed,
Steven W. Allen,
David Alonso,
Marcelo Alvarez,
Mustafa A. Amin,
Adam Anderson,
Kam S. Arnold,
Carlo Baccigalupi,
Kathy Bailey,
Denis Barkats,
Darcy Barron,
Peter S. Barry,
James G. Bartlett,
Ritoban Basu Thakur,
Nicholas Battaglia,
Eric Baxter,
Rachel Bean,
Chris Bebek,
Amy N. Bender,
Bradford A. Benson,
Edo Berger,
Sanah Bhimani
, et al. (200 additional authors not shown)
Abstract:
We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.
We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.
△ Less
Submitted 31 July, 2019;
originally announced August 2019.
-
CosmoDC2: A Synthetic Sky Catalog for Dark Energy Science with LSST
Authors:
Danila Korytov,
Andrew Hearin,
Eve Kovacs,
Patricia Larsen,
Esteban Rangel,
Joseph Hollowed,
Andrew J. Benson,
Katrin Heitmann,
Yao-Yuan Mao,
Anita Bahmanyar,
Chihway Chang,
Duncan Campbell,
Joseph Derose,
Hal Finkel,
Nicholas Frontiere,
Eric Gawiser,
Salman Habib,
Benjamin Joachimi,
François Lanusse,
Nan Li,
Rachel Mandelbaum,
Christopher Morrison,
Jeffrey A. Newman,
Adrian Pope,
Eli Rykoff
, et al. (5 additional authors not shown)
Abstract:
This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `…
▽ More
This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `Outer Rim' run. It covers 440 deg^2 of sky area to a redshift of z=3 and is complete to a magnitude depth of 28 in the r-band. Each galaxy is characterized by a multitude of properties including stellar mass, morphology, spectral energy distributions, broadband filter magnitudes, host halo information and weak lensing shear. The size and complexity of cosmoDC2 requires an efficient catalog generation methodology; our approach is based on a new hybrid technique that combines data-driven empirical approaches with semi-analytic galaxy modeling. A wide range of observation-based validation tests has been implemented to ensure that cosmoDC2 enables the science goals of the planned LSST DESC DC2 analyses. This paper also represents the official release of the cosmoDC2 data set, including an efficient reader that facilitates interaction with the data.
△ Less
Submitted 27 July, 2019; v1 submitted 15 July, 2019;
originally announced July 2019.
-
CMB-S4 Science Case, Reference Design, and Project Plan
Authors:
Kevork Abazajian,
Graeme Addison,
Peter Adshead,
Zeeshan Ahmed,
Steven W. Allen,
David Alonso,
Marcelo Alvarez,
Adam Anderson,
Kam S. Arnold,
Carlo Baccigalupi,
Kathy Bailey,
Denis Barkats,
Darcy Barron,
Peter S. Barry,
James G. Bartlett,
Ritoban Basu Thakur,
Nicholas Battaglia,
Eric Baxter,
Rachel Bean,
Chris Bebek,
Amy N. Bender,
Bradford A. Benson,
Edo Berger,
Sanah Bhimani,
Colin A. Bischoff
, et al. (200 additional authors not shown)
Abstract:
We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.
We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.
△ Less
Submitted 9 July, 2019;
originally announced July 2019.
-
Report on LSST Next-generation Instrumentation Workshop, April 11, 12 2019
Authors:
Christopher W. Stubbs,
Katrin Heitmann
Abstract:
The Large Synoptic Survey Telescope (LSST) is a wide-field imaging system of unprecedented etendue. The initial goal of the project is to carry out a ten year imaging survey in six broad passbands (ugrizy) that cover $350 nm < λ< 1.1 μm$. This document reports on the discussions that occurred at workshop (held April 11-12, 2019 at Argonne National Laboratory) that was convened to explore concepts…
▽ More
The Large Synoptic Survey Telescope (LSST) is a wide-field imaging system of unprecedented etendue. The initial goal of the project is to carry out a ten year imaging survey in six broad passbands (ugrizy) that cover $350 nm < λ< 1.1 μm$. This document reports on the discussions that occurred at workshop (held April 11-12, 2019 at Argonne National Laboratory) that was convened to explore concepts for using the LSST system once the initial survey is complete. Participants discussed the tradeoffs in science performance, cost, and uniqueness for i) imaging surveys using the initial wide-field CCD instrument, perhaps supplemented with different filters, ii) replacing the focal plane with some alternative sensor technology, and iii) converting the LSST system to a wide-field multi-object fiber-fed spectrograph. Participants concluded that the fiber spectrograph option would be most effective if the focal plane were to feed upwards of 30,000 fibers. Thermal management and power considerations in the LSST instrument barrel make it difficult to accommodate infrared sensors that require very low operating temperatures, and the current generation of buttable IR sensors that would extend sensitivity to 2 $μ$m are, we concluded, cost-prohibitive. Procuring and using an alternative filter set, on the other hand, is modest in cost, would take full advantage of the LSST image reduction pipeline, and could yield considerable additional information.
△ Less
Submitted 12 May, 2019;
originally announced May 2019.
-
The Outer Rim Simulation: A Path to Many-Core Supercomputers
Authors:
Katrin Heitmann,
Hal Finkel,
Adrian Pope,
Vitali Morozov,
Nicholas Frontiere,
Salman Habib,
Esteban Rangel,
Thomas Uram,
Danila Korytov,
Hillary Child,
Samuel Flender,
Joe Insley,
Silvio Rizzi
Abstract:
We describe the Outer Rim cosmological simulation, one of the largest high-resolution N-body simulations performed to date, aimed at promoting science to be carried out with large-scale structure surveys. The simulation covers a volume of (4.225Gpc)^3 and evolves more than one trillion particles. It was executed on Mira, a BlueGene/Q system at the Argonne Leadership Computing Facility. We discuss…
▽ More
We describe the Outer Rim cosmological simulation, one of the largest high-resolution N-body simulations performed to date, aimed at promoting science to be carried out with large-scale structure surveys. The simulation covers a volume of (4.225Gpc)^3 and evolves more than one trillion particles. It was executed on Mira, a BlueGene/Q system at the Argonne Leadership Computing Facility. We discuss some of the computational challenges posed by a system like Mira, a many-core supercomputer, and how the simulation code, HACC, has been designed to overcome these challenges. We have carried out a large range of analyses on the simulation data and we report on the results as well as the data products that have been generated. The full data set generated by the simulation totals more than 5PB of data, making data curation and data handling a large challenge in of itself. The simulation results have been used to generate synthetic catalogs for large-scale structure surveys, including DESI and eBOSS, as well as CMB experiments. A detailed catalog for the LSST DESC data challenges has been created as well. We publicly release some of the Outer Rim halo catalogs, downsampled particle information, and lightcone data.
△ Less
Submitted 28 April, 2019; v1 submitted 26 April, 2019;
originally announced April 2019.
-
HACC Cosmological Simulations: First Data Release
Authors:
Katrin Heitmann,
Thomas D. Uram,
Hal Finkel,
Nicholas Frontiere,
Salman Habib,
Adrian Pope,
Esteban Rangel,
Joseph Hollowed,
Danila Korytov,
Patricia Larsen,
Benjamin S. Allen,
Kyle Chard,
Ian Foster
Abstract:
We describe the first major public data release from cosmological simulations carried out with Argonne's HACC code. This initial release covers a range of datasets from large gravity-only simulations. The data products include halo information for multiple redshifts, down-sampled particles, and lightcone outputs. We provide data from two very large LCDM simulations as well as beyond-LCDM simulatio…
▽ More
We describe the first major public data release from cosmological simulations carried out with Argonne's HACC code. This initial release covers a range of datasets from large gravity-only simulations. The data products include halo information for multiple redshifts, down-sampled particles, and lightcone outputs. We provide data from two very large LCDM simulations as well as beyond-LCDM simulations spanning eleven w0-wa cosmologies. Our release platform uses Petrel, a research data service, located at the Argonne Leadership Computing Facility. Petrel offers fast data transfer mechanisms and authentication via Globus, enabling simple and efficient access to stored datasets. Easy browsing of the available data products is provided via a web portal that allows the user to navigate simulation products efficiently. The data hub will be extended by adding more types of data products and by enabling computational capabilities to allow direct interactions with simulation results.
△ Less
Submitted 3 October, 2019; v1 submitted 26 April, 2019;
originally announced April 2019.