-
Exploring the Non-Gaussianity of the Cosmic Infrared Background and Its Weak Gravitational Lensing
Authors:
Jaemyoung Lee,
J. Richard Bond,
Pavel Motloch,
Alexander van Engelen,
George Stein
Abstract:
Gravitational lensing deflects the paths of photons, altering the statistics of cosmic backgrounds and distorting their information content. We take the Cosmic Infrared Background (CIB), which provides plentiful information about galaxy formation and evolution, as an example to probe the effect of lensing on non-Gaussian statistics. Using the Websky simulations, we first quantify the non-Gaussiani…
▽ More
Gravitational lensing deflects the paths of photons, altering the statistics of cosmic backgrounds and distorting their information content. We take the Cosmic Infrared Background (CIB), which provides plentiful information about galaxy formation and evolution, as an example to probe the effect of lensing on non-Gaussian statistics. Using the Websky simulations, we first quantify the non-Gaussianity of the CIB, revealing additional detail on top of its well-measured power spectrum. To achieve this, we use needlet-like multipole-band-filters to calculate the variance and higher-point correlations. Using our simulations, we show the 2-point, 3-point and 4-point spectra, and compare our calculated power spectra and bispectra to Planck values. We then lens the CIB, shell-by-shell with corresponding convergence maps, to capture the broad redshift extent of both the CIB and its lensing convergence. The lensing of the CIB changes the 3-point and 4-point functions by a few tens of percent at large scales, unlike with the power spectrum, which changes by less than two percent. We expand our analyses to encompass the full intensity probability distribution functions (PDFs) involving all n-point correlations as a function of scale. In particular, we use the relative entropy between lensed and unlensed PDFs to create a spectrum of templates that can allow estimation of lensing. The underlying CIB model is missing the important role of star-bursting, which we test by adding a stochastic log-normal term to the intensity distributions. The novel aspects of our filtering and lensing pipeline should prove useful for any radiant background, including line intensity maps.
△ Less
Submitted 23 February, 2024; v1 submitted 14 April, 2023;
originally announced April 2023.
-
A Probabilistic Autoencoder for Type Ia Supernovae Spectral Time Series
Authors:
George Stein,
Uros Seljak,
Vanessa Bohm,
G. Aldering,
P. Antilogus,
C. Aragon,
S. Bailey,
C. Baltay,
S. Bongard,
K. Boone,
C. Buton,
Y. Copin,
S. Dixon,
D. Fouchez,
E. Gangler,
R. Gupta,
B. Hayden,
W. Hillebrandt,
M. Karmen,
A. G. Kim,
M. Kowalski,
D. Kusters,
P. F. Leget,
F. Mondon,
J. Nordin
, et al. (15 additional authors not shown)
Abstract:
We construct a physically-parameterized probabilistic autoencoder (PAE) to learn the intrinsic diversity of type Ia supernovae (SNe Ia) from a sparse set of spectral time series. The PAE is a two-stage generative model, composed of an Auto-Encoder (AE) which is interpreted probabilistically after training using a Normalizing Flow (NF). We demonstrate that the PAE learns a low-dimensional latent sp…
▽ More
We construct a physically-parameterized probabilistic autoencoder (PAE) to learn the intrinsic diversity of type Ia supernovae (SNe Ia) from a sparse set of spectral time series. The PAE is a two-stage generative model, composed of an Auto-Encoder (AE) which is interpreted probabilistically after training using a Normalizing Flow (NF). We demonstrate that the PAE learns a low-dimensional latent space that captures the nonlinear range of features that exists within the population, and can accurately model the spectral evolution of SNe Ia across the full range of wavelength and observation times directly from the data. By introducing a correlation penalty term and multi-stage training setup alongside our physically-parameterized network we show that intrinsic and extrinsic modes of variability can be separated during training, removing the need for the additional models to perform magnitude standardization. We then use our PAE in a number of downstream tasks on SNe Ia for increasingly precise cosmological analyses, including automatic detection of SN outliers, the generation of samples consistent with the data distribution, and solving the inverse problem in the presence of noisy and incomplete data to constrain cosmological distance measurements. We find that the optimal number of intrinsic model parameters appears to be three, in line with previous studies, and show that we can standardize our test sample of SNe Ia with an RMS of $0.091 \pm 0.010$ mag, which corresponds to $0.074 \pm 0.010$ mag if peculiar velocity contributions are removed. Trained models and codes are released at \href{https://github.com/georgestein/suPAErnova}{github.com/georgestein/suPAErnova}
△ Less
Submitted 15 July, 2022;
originally announced July 2022.
-
Self-supervised similarity search for large scientific datasets
Authors:
George Stein,
Peter Harrington,
Jacqueline Blaum,
Tomislav Medan,
Zarija Lukic
Abstract:
We present the use of self-supervised learning to explore and exploit large unlabeled datasets. Focusing on 42 million galaxy images from the latest data release of the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys, we first train a self-supervised model to distill low-dimensional representations that are robust to symmetries, uncertainties, and noise in each image. We then us…
▽ More
We present the use of self-supervised learning to explore and exploit large unlabeled datasets. Focusing on 42 million galaxy images from the latest data release of the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys, we first train a self-supervised model to distill low-dimensional representations that are robust to symmetries, uncertainties, and noise in each image. We then use the representations to construct and publicly release an interactive semantic similarity search tool. We demonstrate how our tool can be used to rapidly discover rare objects given only a single example, increase the speed of crowd-sourcing campaigns, and construct and improve training sets for supervised applications. While we focus on images from sky surveys, the technique is straightforward to apply to any scientific dataset of any dimensionality. The similarity search web app can be found at https://github.com/georgestein/galaxy_search
△ Less
Submitted 30 November, 2021; v1 submitted 25 October, 2021;
originally announced October 2021.
-
Mining for Strong Gravitational Lenses with Self-supervised Learning
Authors:
George Stein,
Jacqueline Blaum,
Peter Harrington,
Tomislav Medan,
Zarija Lukic
Abstract:
We employ self-supervised representation learning to distill information from 76 million galaxy images from the Dark Energy Spectroscopic Instrument Legacy Imaging Surveys' Data Release 9. Targeting the identification of new strong gravitational lens candidates, we first create a rapid similarity search tool to discover new strong lenses given only a single labelled example. We then show how train…
▽ More
We employ self-supervised representation learning to distill information from 76 million galaxy images from the Dark Energy Spectroscopic Instrument Legacy Imaging Surveys' Data Release 9. Targeting the identification of new strong gravitational lens candidates, we first create a rapid similarity search tool to discover new strong lenses given only a single labelled example. We then show how training a simple linear classifier on the self-supervised representations, requiring only a few minutes on a CPU, can automatically classify strong lenses with great efficiency. We present 1192 new strong lens candidates that we identified through a brief visual identification campaign, and release an interactive web-based similarity search tool and the top network predictions to facilitate crowd-sourcing rapid discovery of additional strong gravitational lenses and other rare objects: https://github.com/georgestein/ssl-legacysurvey.
△ Less
Submitted 21 June, 2022; v1 submitted 30 September, 2021;
originally announced October 2021.
-
Superclustering with the Atacama Cosmology Telescope and Dark Energy Survey: I. Evidence for thermal energy anisotropy using oriented stacking
Authors:
M. Lokken,
R. Hložek,
A. van Engelen,
M. Madhavacheril,
E. Baxter,
J. DeRose,
C. Doux,
S. Pandey,
E. S. Rykoff,
G. Stein,
C. To,
T. M. C. Abbott,
S. Adhikari,
M. Aguena,
S. Allam,
F. Andrade-Oliveira,
J. Annis,
N. Battaglia,
G. M. Bernstein,
E. Bertin,
J. R. Bond,
D. Brooks,
E. Calabrese,
A. Carnero Rosell,
M. Carrasco Kind
, et al. (82 additional authors not shown)
Abstract:
The cosmic web contains filamentary structure on a wide range of scales. On the largest scales, superclustering aligns multiple galaxy clusters along inter-cluster bridges, visible through their thermal Sunyaev-Zel'dovich signal in the Cosmic Microwave Background. We demonstrate a new, flexible method to analyze the hot gas signal from multi-scale extended structures. We use a Compton-$y$ map from…
▽ More
The cosmic web contains filamentary structure on a wide range of scales. On the largest scales, superclustering aligns multiple galaxy clusters along inter-cluster bridges, visible through their thermal Sunyaev-Zel'dovich signal in the Cosmic Microwave Background. We demonstrate a new, flexible method to analyze the hot gas signal from multi-scale extended structures. We use a Compton-$y$ map from the Atacama Cosmology Telescope (ACT) stacked on redMaPPer cluster positions from the optical Dark Energy Survey (DES). Cutout images from the $y$ map are oriented with large-scale structure information from DES galaxy data such that the superclustering signal is aligned before being overlaid. We find evidence for an extended quadrupole moment of the stacked $y$ signal at the 3.5$σ$ level, demonstrating that the large-scale thermal energy surrounding galaxy clusters is anisotropically distributed. We compare our ACT$\times$DES results with the Buzzard simulations, finding broad agreement. Using simulations, we highlight the promise of this novel technique for constraining the evolution of anisotropic, non-Gaussian structure using future combinations of microwave and optical surveys.
△ Less
Submitted 18 July, 2022; v1 submitted 12 July, 2021;
originally announced July 2021.
-
Estimating Galactic Distances From Images Using Self-supervised Representation Learning
Authors:
Md Abul Hayat,
Peter Harrington,
George Stein,
Zarija Lukić,
Mustafa Mustafa
Abstract:
We use a contrastive self-supervised learning framework to estimate distances to galaxies from their photometric images. We incorporate data augmentations from computer vision as well as an application-specific augmentation accounting for galactic dust. We find that the resulting visual representations of galaxy images are semantically useful and allow for fast similarity searches, and can be succ…
▽ More
We use a contrastive self-supervised learning framework to estimate distances to galaxies from their photometric images. We incorporate data augmentations from computer vision as well as an application-specific augmentation accounting for galactic dust. We find that the resulting visual representations of galaxy images are semantically useful and allow for fast similarity searches, and can be successfully fine-tuned for the task of redshift estimation. We show that (1) pretraining on a large corpus of unlabeled data followed by fine-tuning on some labels can attain the accuracy of a fully-supervised model which requires 2-4x more labeled data, and (2) that by fine-tuning our self-supervised representations using all available data labels in the Main Galaxy Sample of the Sloan Digital Sky Survey (SDSS), we outperform the state-of-the-art supervised learning method.
△ Less
Submitted 11 January, 2021;
originally announced January 2021.
-
Statistical exploration of halo anisotropic clustering and intrinsic alignments with the mass-Peak Patch algorithm
Authors:
Bruno Regaldo-Saint Blancard,
Sandrine Codis,
J. Richard Bond,
George Stein
Abstract:
The anisotropy or triaxiality of massive dark matter haloes largely defines the structure of the cosmic web, in particular the filaments that join the haloes together. Here we investigate such oriented correlations in mass-Peak Patch halo catalogues by using the initial strain tensor of spherical proto-halo regions to orient the haloes. To go beyond the spherically averaged two-point correlation f…
▽ More
The anisotropy or triaxiality of massive dark matter haloes largely defines the structure of the cosmic web, in particular the filaments that join the haloes together. Here we investigate such oriented correlations in mass-Peak Patch halo catalogues by using the initial strain tensor of spherical proto-halo regions to orient the haloes. To go beyond the spherically averaged two-point correlation function of haloes we use oriented stacks to compute oriented two-point correlations: we explicitly break isotropy by imposing a local frame set by the strain tensor of the reference halo before stacking neighbouring haloes. Beyond the exclusion zone of the reference halo, clustering is found to be strongly enhanced along the major direction of the strain tensor as expected. This anisotropic clustering of haloes along filaments is further quantified by using a spherical harmonics decomposition. Furthermore, we compute the evolution of cluster-scale halo principal directions relative to those of their neighbours and show that there are strong correlations extending up to very large scales. In order to provide calculations more suitable to observational confrontations, we also utilize 2D projected versions of some equivalent correlation functions. Finally, we show that the multipole structure of the mass-peak patch halo's anisotropic clustering can be qualitatively captured in an analytic treatment based on peak theory. Though highly informative, analytic evaluation involves extensive use of Monte Carlo methods, which is also what the simulated catalogue uses, taking into account as they do the adaptive nature of the mass-peak patch mass hierarchy and all non-local complexities associated with the exclusion of smaller haloes overlapping with larger ones: there is no substitute for the mass-Peak Patch simulation-based determination of oriented and anisotropic correlations.
△ Less
Submitted 30 March, 2021; v1 submitted 5 January, 2021;
originally announced January 2021.
-
Self-Supervised Representation Learning for Astronomical Images
Authors:
Md Abul Hayat,
George Stein,
Peter Harrington,
Zarija Lukić,
Mustafa Mustafa
Abstract:
Sky surveys are the largest data generators in astronomy, making automated tools for extracting meaningful scientific information an absolute necessity. We show that, without the need for labels, self-supervised learning recovers representations of sky survey images that are semantically useful for a variety of scientific tasks. These representations can be directly used as features, or fine-tuned…
▽ More
Sky surveys are the largest data generators in astronomy, making automated tools for extracting meaningful scientific information an absolute necessity. We show that, without the need for labels, self-supervised learning recovers representations of sky survey images that are semantically useful for a variety of scientific tasks. These representations can be directly used as features, or fine-tuned, to outperform supervised methods trained only on labeled data. We apply a contrastive learning framework on multi-band galaxy photometry from the Sloan Digital Sky Survey (SDSS) to learn image representations. We then use them for galaxy morphology classification, and fine-tune them for photometric redshift estimation, using labels from the Galaxy Zoo 2 dataset and SDSS spectroscopy. In both downstream tasks, using the same learned representations, we outperform the supervised state-of-the-art results, and we show that our approach can achieve the accuracy of supervised models while using 2-4 times fewer labels for training.
△ Less
Submitted 8 April, 2021; v1 submitted 23 December, 2020;
originally announced December 2020.
-
Nonlinear 3D Cosmic Web Simulation with Heavy-Tailed Generative Adversarial Networks
Authors:
Richard M. Feder,
Philippe Berger,
George Stein
Abstract:
Fast and accurate simulations of the non-linear evolution of the cosmic density field are a major component of many cosmological analyses, but the computational time and storage required to run them can be exceedingly large. For this reason, we use generative adversarial networks (GANs) to learn a compressed representation of the 3D matter density field that is fast and easy to sample, and for the…
▽ More
Fast and accurate simulations of the non-linear evolution of the cosmic density field are a major component of many cosmological analyses, but the computational time and storage required to run them can be exceedingly large. For this reason, we use generative adversarial networks (GANs) to learn a compressed representation of the 3D matter density field that is fast and easy to sample, and for the first time show that GANs are capable of generating samples at the level of accuracy of other conventional methods. Using sub-volumes from a suite of GADGET-2 N-body simulations, we demonstrate that a deep-convolutional GAN can generate samples that capture both large- and small-scale features of the matter density field, as validated through a variety of n-point statistics. The use of a data scaling that preserves high-density features and a heavy-tailed latent space prior allow us to obtain state of the art results for fast 3D cosmic web generation. In particular, the mean power spectra from generated samples agree to within 5% up to k=3 and within 10% for k<5 when compared with N-body simulations, and similar accuracy is obtained for a variety of bispectra. By modeling the latent space with a heavy-tailed prior rather than a standard Gaussian, we better capture sample variance in the high-density voxel PDF and reduce errors in power spectrum and bispectrum covariance on all scales. Furthermore, we show that a conditional GAN can smoothly interpolate between samples conditioned on redshift. Deep generative models, such as the ones described in this work, provide great promise as fast, low-memory, high-fidelity forward models of large-scale structure.
△ Less
Submitted 13 November, 2020; v1 submitted 6 May, 2020;
originally announced May 2020.
-
The Websky Extragalactic CMB Simulations
Authors:
George Stein,
Marcelo A. Alvarez,
J. Richard Bond,
Alexander van Engelen,
Nicholas Battaglia
Abstract:
We present a new pipeline for the efficient generation of synthetic observations of the extragalactic microwave sky, tailored to large ground-based CMB experiments such as the Simons Observatory, Advanced ACTPol, SPT-3G, and CMB-S4. Such simulated observations are a key technical challenge in cosmology because of the dynamic range and accuracy required. The first part of the pipeline generates a r…
▽ More
We present a new pipeline for the efficient generation of synthetic observations of the extragalactic microwave sky, tailored to large ground-based CMB experiments such as the Simons Observatory, Advanced ACTPol, SPT-3G, and CMB-S4. Such simulated observations are a key technical challenge in cosmology because of the dynamic range and accuracy required. The first part of the pipeline generates a random cosmological realization in the form of a dark matter halo catalog and matter displacement field, as seen from a given position. The halo catalog and displacement field are modeled with ellipsoidal collapse dynamics and Lagrangian perturbation theory, respectively. In the second part, the cosmological realization is converted into a set of intensity maps over the range 10 - 10^3 GHz using models based on existing observations and hydrodynamical simulations. These maps include infrared emission from dusty star forming galaxies (CIB), Comptonization of CMB photons by hot gas in groups and clusters through the thermal Sunyaev-Zel'dovich effect (tSZ), Doppler boosting by Thomson scattering of the CMB by bulk flows through the kinetic Sunyaev-Zel'dovich effect (kSZ), and weak gravitational lensing of primary CMB anisotropies by the large-scale distribution of matter in the universe. After describing the pipeline and its implementation, we present the Websky maps, created from a realization of the cosmic web on our past light cone in the redshift interval 0<z<4.6 over the full-sky and a volume of ~(600 Gpc/h)^3 resolved with ~10^12 resolution elements. The Websky maps and halo catalog are publicly available at mocks.cita.utoronto.ca/websky.
△ Less
Submitted 23 January, 2020;
originally announced January 2020.
-
The Simons Observatory: Astro2020 Decadal Project Whitepaper
Authors:
The Simons Observatory Collaboration,
Maximilian H. Abitbol,
Shunsuke Adachi,
Peter Ade,
James Aguirre,
Zeeshan Ahmed,
Simone Aiola,
Aamir Ali,
David Alonso,
Marcelo A. Alvarez,
Kam Arnold,
Peter Ashton,
Zachary Atkins,
Jason Austermann,
Humna Awan,
Carlo Baccigalupi,
Taylor Baildon,
Anton Baleato Lizancos,
Darcy Barron,
Nick Battaglia,
Richard Battye,
Eric Baxter,
Andrew Bazarko,
James A. Beall,
Rachel Bean
, et al. (258 additional authors not shown)
Abstract:
The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021…
▽ More
The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs.
The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation.
With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.
△ Less
Submitted 16 July, 2019;
originally announced July 2019.
-
Deconfusing intensity maps with neural networks
Authors:
Daniel N. Pfeffer,
Patrick C. Breysse,
George Stein
Abstract:
Line intensity maps (LIMs) are in principle sensitive to a large amount of information about faint, distant galaxies which are invisible to conventional surveys. However, actually extracting that information from a confused, foreground-contaminated map can be challenging. In this work we present the first application of convolutional neural network (CNN) to directly determine the underlying lumino…
▽ More
Line intensity maps (LIMs) are in principle sensitive to a large amount of information about faint, distant galaxies which are invisible to conventional surveys. However, actually extracting that information from a confused, foreground-contaminated map can be challenging. In this work we present the first application of convolutional neural network (CNN) to directly determine the underlying luminosity function of a LIM, including a treatment of extragalactic foregrounds and instrumental noise. We apply the CNN to simulations of mock Carbon Monoxide (CO) line intensity maps similar to those which will be produced by the currently-active COMAP experiment. We evaluate the trained CNN on a number of noise scenarios in order to determine how robust the network predictions are for application to realistic data. We find that, in the ideal case where the mock data capture all of the features of the real data, the CNN performs comparably to or better than conventional analyses. However, the network's accuracy degrades considerably when tested on signals and systematics outside of those it was trained on. For both intensity mapping and cosmology as a whole, this motivates a broad-based study of whether simulated data can ever be generated with sufficient detail to realize the enormous potential of machine learning methods.
△ Less
Submitted 24 May, 2019;
originally announced May 2019.
-
Measurement of the Splashback Feature around SZ-selected Galaxy Clusters with DES, SPT and ACT
Authors:
T. Shin,
S. Adhikari,
E. J. Baxter,
C. Chang,
B. Jain,
N. Battaglia,
L. Bleem,
S. Bocquet,
J. DeRose,
D. Gruen,
M. Hilton,
A. Kravtsov,
T. McClintock,
E. Rozo,
E. S. Rykoff,
T. N. Varga,
R. H. Wechsler,
H. Wu,
S. Aiola,
S. Allam,
K. Bechtol,
B. A. Benson,
E. Bertin,
J. R. Bond,
M. Brodwin
, et al. (85 additional authors not shown)
Abstract:
We present a detection of the splashback feature around galaxy clusters selected using their Sunyaev-Zel'dovich (SZ) signal. Recent measurements of the splashback feature around optically selected galaxy clusters have found that the splashback radius, $r_{\rm sp}$, is smaller than predicted by N-body simulations. A possible explanation for this discrepancy is that $r_{\rm sp}$ inferred from the ob…
▽ More
We present a detection of the splashback feature around galaxy clusters selected using their Sunyaev-Zel'dovich (SZ) signal. Recent measurements of the splashback feature around optically selected galaxy clusters have found that the splashback radius, $r_{\rm sp}$, is smaller than predicted by N-body simulations. A possible explanation for this discrepancy is that $r_{\rm sp}$ inferred from the observed radial distribution of galaxies is affected by selection effects related to the optical cluster-finding algorithms. We test this possibility by measuring the splashback feature in clusters selected via the SZ effect in data from the South Pole Telescope SZ survey and the Atacama Cosmology Telescope Polarimeter survey. The measurement is accomplished by correlating these clusters with galaxies detected in the Dark Energy Survey Year 3 data. The SZ observable used to select clusters in this analysis is expected to have a tighter correlation with halo mass and to be more immune to projection effects and aperture-induced biases than optically selected clusters. We find that the measured $r_{\rm sp}$ for SZ-selected clusters is consistent with the expectations from simulations, although the small number of SZ-selected clusters makes a precise comparison difficult. In agreement with previous work, when using optically selected redMaPPer clusters, $r_{\rm sp}$ is $\sim$ $2σ$ smaller than in the simulations. These results motivate detailed investigations of selection biases in optically selected cluster catalogs and exploration of the splashback feature around larger samples of SZ-selected clusters. Additionally, we investigate trends in the galaxy profile and splashback feature as a function of galaxy color, finding that blue galaxies have profiles close to a power law with no discernible splashback feature, which is consistent with them being on their first infall into the cluster.
△ Less
Submitted 24 May, 2019; v1 submitted 14 November, 2018;
originally announced November 2018.
-
The mass-Peak Patch algorithm for fast generation of deep all-sky dark matter halo catalogues and its N-Body validation
Authors:
George Stein,
Marcelo A. Alvarez,
J. Richard Bond
Abstract:
We present a detailed description and validation of our massively-parallel update to the mass-Peak Patch method, a fully predictive initial-space algorithm to quickly generate dark matter halo catalogues in very large cosmological volumes. We perform an extensive systematic comparison to a suite of N-body simulations covering a broad range of redshifts and simulation resolutions, and find that, wi…
▽ More
We present a detailed description and validation of our massively-parallel update to the mass-Peak Patch method, a fully predictive initial-space algorithm to quickly generate dark matter halo catalogues in very large cosmological volumes. We perform an extensive systematic comparison to a suite of N-body simulations covering a broad range of redshifts and simulation resolutions, and find that, without any parameter fitting, our method is able to generally reproduce N-body results while typically using over 3 orders of magnitude less CPU time, and a fraction of the memory cost. Instead of calculating the full non-linear gravitational collapse determined by an N-body simulation, the mass-Peak Patch method finds an overcomplete set of just-collapsed structures around a hierarchy of density-peak points by coarse-grained (homogeneous) ellipsoidal dynamics. A complete set of mass-peaks, or halos, is then determined by exclusion of overlapping patches, and second order Lagrangian displacements are used to move the halos to their final positions and to give their flow velocities. Our results show that the mass-Peak Patch method is well-suited for creating large ensembles of halo catalogues to mock cosmological surveys, and to aid in complex statistical interpretations of cosmological models.
△ Less
Submitted 17 October, 2018;
originally announced October 2018.
-
Cross-correlating Carbon Monoxide Line-intensity Maps with Spectroscopic and Photometric Galaxy Surveys
Authors:
Dongwoo T. Chung,
Marco P. Viero,
Sarah E. Church,
Risa H. Wechsler,
Marcelo A. Alvarez,
J. Richard Bond,
Patrick C. Breysse,
Kieran A. Cleary,
Hans K. Eriksen,
Marie K. Foss,
Joshua O. Gundersen,
Stuart E. Harper,
Håvard T. Ihle,
Laura C. Keating,
Norman Murray,
Hamsa Padmanabhan,
George F. Stein,
Ingunn K. Wehus
Abstract:
Line-intensity mapping (LIM or IM) is an emerging field of observational work, with strong potential to fit into a larger effort to probe large-scale structure and small-scale astrophysical phenomena using multiple complementary tracers. Taking full advantage of such complementarity means, in part, undertaking line-intensity surveys with galaxy surveys in mind. We consider the potential for detect…
▽ More
Line-intensity mapping (LIM or IM) is an emerging field of observational work, with strong potential to fit into a larger effort to probe large-scale structure and small-scale astrophysical phenomena using multiple complementary tracers. Taking full advantage of such complementarity means, in part, undertaking line-intensity surveys with galaxy surveys in mind. We consider the potential for detection of a cross-correlation signal between COMAP and blind surveys based on photometric redshifts (as in COSMOS) or based on spectroscopic data (as with the HETDEX survey of Lyman-$α$ emitters). We find that obtaining $σ_z/(1+z)\lesssim0.003$ accuracy in redshifts and $\gtrsim10^{-4}$ sources per Mpc$^3$ with spectroscopic redshift determination should enable a CO-galaxy cross spectrum detection significance at least twice that of the CO auto spectrum. Either a future targeted spectroscopic survey or a blind survey like HETDEX may be able to meet both of these requirements.
△ Less
Submitted 17 January, 2019; v1 submitted 12 September, 2018;
originally announced September 2018.
-
Joint power spectrum and voxel intensity distribution forecast on the CO luminosity function with COMAP
Authors:
Håvard Tveit Ihle,
Dongwoo Chung,
George Stein,
Marcelo Alvarez,
J. Richard Bond,
Patrick C. Breysse,
Kieran A. Cleary,
Hans Kristian Eriksen,
Marie Kristine Foss,
Joshua Ott Gundersen,
Stuart Harper,
Norman Murray,
Hamsa Padmanabhan,
Marco P. Viero,
Ingunn Katerine Wehus
Abstract:
We develop a framework for joint constraints on the CO luminosity function based on power spectra (PS) and voxel intensity distributions (VID), and apply this to simulations of COMAP, a CO intensity mapping experiment. This Bayesian framework is based on a Markov chain Monte Carlo (MCMC) sampler coupled to a Gaussian likelihood with a joint PS + VID covariance matrix computed from a large number o…
▽ More
We develop a framework for joint constraints on the CO luminosity function based on power spectra (PS) and voxel intensity distributions (VID), and apply this to simulations of COMAP, a CO intensity mapping experiment. This Bayesian framework is based on a Markov chain Monte Carlo (MCMC) sampler coupled to a Gaussian likelihood with a joint PS + VID covariance matrix computed from a large number of fiducial simulations, and re-calibrated with a small number of simulations per MCMC step. The simulations are based on dark matter halos from fast peak patch simulations combined with the $L_\text{CO}(M_\text{halo})$ model of Li et al. (2016). We find that the relative power to constrain the CO luminosity function depends on the luminosity range of interest. In particular, the VID is more sensitive at both small and large luminosities, while the PS is more sensitive at intermediate luminosities. The joint analysis is superior to using either observable separately. When averaging over CO luminosities ranging between $L_\text{CO} = 10^4-10^7L_\odot$, and over 10 cosmological realizations of COMAP Phase 2, the uncertainties (in dex) are larger by 58 % and 30 % for the PS and VID, respectively, when compared to the joint analysis (PS + VID). This method is generally applicable to any other random field, with a complicated likelihood, as long a fast simulation procedure is available.
△ Less
Submitted 20 March, 2019; v1 submitted 22 August, 2018;
originally announced August 2018.
-
The Simons Observatory: Science goals and forecasts
Authors:
The Simons Observatory Collaboration,
Peter Ade,
James Aguirre,
Zeeshan Ahmed,
Simone Aiola,
Aamir Ali,
David Alonso,
Marcelo A. Alvarez,
Kam Arnold,
Peter Ashton,
Jason Austermann,
Humna Awan,
Carlo Baccigalupi,
Taylor Baildon,
Darcy Barron,
Nick Battaglia,
Richard Battye,
Eric Baxter,
Andrew Bazarko,
James A. Beall,
Rachel Bean,
Dominic Beck,
Shawn Beckman,
Benjamin Beringue,
Federico Bianchini
, et al. (225 additional authors not shown)
Abstract:
The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225…
▽ More
The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $μ$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $σ(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $μ$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.
△ Less
Submitted 1 March, 2019; v1 submitted 22 August, 2018;
originally announced August 2018.
-
CCAT-prime: Science with an Ultra-widefield Submillimeter Observatory at Cerro Chajnantor
Authors:
G. J. Stacey,
M. Aravena,
K. Basu,
N. Battaglia,
B. Beringue,
F. Bertoldi,
J. R. Bond,
P. Breysse,
R. Bustos,
S. Chapman,
D. T. Chung,
N. Cothard,
J. Erler,
M. Fich,
S. Foreman,
P. Gallardo,
R. Giovanelli,
U. U. Graf,
M. P. Haynes,
R. Herrera-Camus,
T. L. Herter,
R. Hložek,
D. Johnstone,
L. Keating,
B. Magnelli
, et al. (15 additional authors not shown)
Abstract:
We present the detailed science case, and brief descriptions of the telescope design, site, and first light instrument plans for a new ultra-wide field submillimeter observatory, CCAT-prime, that we are constructing at a 5600 m elevation site on Cerro Chajnantor in northern Chile. Our science goals are to study star and galaxy formation from the epoch of reionization to the present, investigate th…
▽ More
We present the detailed science case, and brief descriptions of the telescope design, site, and first light instrument plans for a new ultra-wide field submillimeter observatory, CCAT-prime, that we are constructing at a 5600 m elevation site on Cerro Chajnantor in northern Chile. Our science goals are to study star and galaxy formation from the epoch of reionization to the present, investigate the growth of structure in the Universe, improve the precision of B-mode CMB measurements, and investigate the interstellar medium and star formation in the Galaxy and nearby galaxies through spectroscopic, polarimetric, and broadband surveys at wavelengths from 200 um to 2 mm. These goals are realized with our two first light instruments, a large field-of-view (FoV) bolometer-based imager called Prime-Cam (that has both camera and an imaging spectrometer modules), and a multi-beam submillimeter heterodyne spectrometer, CHAI. CCAT-prime will have very high surface accuracy and very low system emissivity, so that combined with its wide FoV at the unsurpassed CCAT site our telescope/instrumentation combination is ideally suited to pursue this science. The CCAT-prime telescope is being designed and built by Vertex Antennentechnik GmbH. We expect to achieve first light in the spring of 2021.
△ Less
Submitted 11 July, 2018;
originally announced July 2018.
-
Comparing approximate methods for mock catalogues and covariance matrices III: Bispectrum
Authors:
Manuel Colavincenzo,
Emiliano Sefusatti,
Pierluigi Monaco,
Linda Blot,
Martin Crocce,
Martha Lippich,
Ariel G. Sánchez,
Marcelo A. Alvarez,
Aniket Agrawal,
Santiago Avila,
Andrés Balaguera-Antolínez,
Richard Bond,
Sandrine Codis,
Claudio Dalla Vecchia,
Antonio Dorta,
Pablo Fosalba,
Albert Izard,
Francisco-Shu Kitaura,
Marcos Pellejero-Ibanez,
George Stein,
Mohammadjavad Vakili,
Gustavo Yepes
Abstract:
We compare the measurements of the bispectrum and the estimate of its covariance obtained from a set of different methods for the efficient generation of approximate dark matter halo catalogs to the same quantities obtained from full N-body simulations. To this purpose we employ a large set of three-hundred realisations of the same cosmology for each method, run with matching initial conditions in…
▽ More
We compare the measurements of the bispectrum and the estimate of its covariance obtained from a set of different methods for the efficient generation of approximate dark matter halo catalogs to the same quantities obtained from full N-body simulations. To this purpose we employ a large set of three-hundred realisations of the same cosmology for each method, run with matching initial conditions in order to reduce the contribution of cosmic variance to the comparison. In addition, we compare how the error on cosmological parameters such as linear and nonlinear bias parameters depends on the approximate method used for the determination of the bispectrum variance. As general result, most methods provide errors within 10% of the errors estimated from N-body simulations. Exceptions are those methods requiring calibration of the clustering amplitude but restrict this to two-point statistics. Finally we test how our results are affected by being limited to a few hundreds measurements from N-body simulation, and therefore to the bispectrum variance, by comparing with a larger set of several thousands realisations performed with one approximate method.
△ Less
Submitted 8 October, 2018; v1 submitted 25 June, 2018;
originally announced June 2018.
-
Comparing approximate methods for mock catalogues and covariance matrices II: Power spectrum multipoles
Authors:
Linda Blot,
Martin Crocce,
Emiliano Sefusatti,
Martha Lippich,
Ariel G. Sánchez,
Manuel Colavincenzo,
Pierluigi Monaco,
Marcelo A. Alvarez,
Aniket Agrawal,
Santiago Avila,
Andrés Balaguera-Antolínez,
Richard Bond,
Sandrine Codis,
Claudio Dalla Vecchia,
Antonio Dorta,
Pablo Fosalba,
Albert Izard,
Francisco-Shu Kitaura,
Marcos Pellejero-Ibanez,
George Stein,
Mohammadjavad Vakili,
Gustavo Yepes
Abstract:
We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences in covariances into parameter constrains related to growth rate of structure, Alcock-Paczynski distortions and biasing. We consider seven methods in three broad categories: algorithms that solve for halo densi…
▽ More
We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences in covariances into parameter constrains related to growth rate of structure, Alcock-Paczynski distortions and biasing. We consider seven methods in three broad categories: algorithms that solve for halo density evolution deterministically using Lagrangian trajectories (ICE-COLA, Pinocchio and PeakPatch), methods that rely on halo assignment schemes onto dark-matter overdensities calibrated with a target N-body run (Halogen, Patchy) and two standard assumptions about the full density PDF (Gaussian and Lognormal). We benchmark their performance against a set of three hundred N-body simulations, running similar sets of approximate simulations with matched initial conditions, for each method. We find that most methods reproduce the monopole to within $5\%$, while residuals for the quadrupole are sometimes larger and scale dependent. The variance of the multipoles is typically reproduced within $10\%$. Overall, we find that covariances built from approximate simulations yield errors on model parameters within $10\%$ of those from the N-body based covariance.
△ Less
Submitted 18 February, 2019; v1 submitted 25 June, 2018;
originally announced June 2018.
-
Comparing approximate methods for mock catalogues and covariance matrices I: correlation function
Authors:
Martha Lippich,
Ariel G. Sánchez,
Manuel Colavincenzo,
Emiliano Sefusatti,
Pierluigi Monaco,
Linda Blot,
Martin Crocce,
Marcelo A. Alvarez,
Aniket Agrawal,
Santiago Avila,
Andrés Balaguera-Antolínez,
Richard Bond,
Sandrine Codis,
Claudio Dalla Vecchia,
Antonio Dorta,
Pablo Fosalba,
Albert Izard,
Francisco-Shu Kitaura,
Marcos Pellejero-Ibanez,
George Stein,
Mohammadjavad Vakili,
Gustavo Yepes
Abstract:
This paper is the first in a set that analyses the covariance matrices of clustering statistics obtained from several approximate methods for gravitational structure formation. We focus here on the covariance matrices of anisotropic two-point correlation function measurements. Our comparison includes seven approximate methods, which can be divided into three categories: predictive methods that fol…
▽ More
This paper is the first in a set that analyses the covariance matrices of clustering statistics obtained from several approximate methods for gravitational structure formation. We focus here on the covariance matrices of anisotropic two-point correlation function measurements. Our comparison includes seven approximate methods, which can be divided into three categories: predictive methods that follow the evolution of the linear density field deterministically (ICE-COLA, Peak Patch, and Pinocchio), methods that require a calibration with N-body simulations (Patchy and Halogen), and simpler recipes based on assumptions regarding the shape of the probability distribution function (PDF) of density fluctuations (log-normal and Gaussian density fields). We analyse the impact of using covariance estimates obtained from these approximate methods on cosmological analyses of galaxy clustering measurements, using as a reference the covariances inferred from a set of full N-body simulations. We find that all approximate methods can accurately recover the mean parameter values inferred using the N-body covariances. The obtained parameter uncertainties typically agree with the corresponding N-body results within 5% for our lower mass threshold, and 10% for our higher mass threshold. Furthermore, we find that the constraints for some methods can differ by up to 20% depending on whether the halo samples used to define the covariance matrices are defined by matching the mass, number density, or clustering amplitude of the parent N-body samples. The results of our configuration-space analysis indicate that most approximate methods provide similar results, with no single method clearly outperforming the others.
△ Less
Submitted 13 May, 2019; v1 submitted 25 June, 2018;
originally announced June 2018.
-
A volumetric deep Convolutional Neural Network for simulation of mock dark matter halo catalogues
Authors:
Philippe Berger,
George Stein
Abstract:
For modern large-scale structure survey techniques it has become standard practice to test data analysis pipelines on large suites of mock simulations, a task which is currently prohibitively expensive for full N-body simulations. Instead of calculating this costly gravitational evolution, we have trained a three-dimensional deep Convolutional Neural Network (CNN) to identify dark matter protohalo…
▽ More
For modern large-scale structure survey techniques it has become standard practice to test data analysis pipelines on large suites of mock simulations, a task which is currently prohibitively expensive for full N-body simulations. Instead of calculating this costly gravitational evolution, we have trained a three-dimensional deep Convolutional Neural Network (CNN) to identify dark matter protohalos directly from the cosmological initial conditions. Training on halo catalogues from the Peak Patch semi-analytic code, we test various CNN architectures and find they generically achieve a Dice coefficient of ~92% in only 24 hours of training. We present a simple and fast geometric halo finding algorithm to extract halos from this powerful pixel-wise binary classifier and find that the predicted catalogues match the mass function and power spectra of the ground truth simulations to within ~10%. We investigate the effect of long-range tidal forces on an object-by-object basis and find that the network's predictions are consistent with the non-linear ellipsoidal collapse equations used explicitly by the Peak Patch algorithm.
△ Less
Submitted 19 November, 2018; v1 submitted 11 May, 2018;
originally announced May 2018.
-
Weak-Lensing Mass Calibration of ACTPol Sunyaev-Zel'dovich Clusters with the Hyper Suprime-Cam Survey
Authors:
Hironao Miyatake,
Nicholas Battaglia,
Matt Hilton,
Elinor Medezinski,
Atsushi J. Nishizawa,
Surhud More,
Simone Aiola,
Neta Bahcall,
J. Richard Bond,
Erminia Calabrese,
Steve K. Choi,
Mark J. Devlin,
Joanna Dunkley,
Rolando Dunner,
Brittany Fuzia,
Patricio Gallardo,
Megan Gralla,
Matthew Hasselfield,
Mark Halpern,
Chiaki Hikage,
J. Colin Hill,
Adam D. Hincks,
Renée Hložek,
Kevin Huffenberger,
John P. Hughes
, et al. (35 additional authors not shown)
Abstract:
We present weak-lensing measurements using the first-year data from the Hyper Suprime-Cam Strategic Survey Program on the Subaru telescope for eight galaxy clusters selected through their thermal Sunyaev-Zel'dovich (SZ) signal measured at 148 GHz with the Atacama Cosmology Telescope Polarimeter experiment. The overlap between the two surveys in this work is 33.8 square degrees, before masking brig…
▽ More
We present weak-lensing measurements using the first-year data from the Hyper Suprime-Cam Strategic Survey Program on the Subaru telescope for eight galaxy clusters selected through their thermal Sunyaev-Zel'dovich (SZ) signal measured at 148 GHz with the Atacama Cosmology Telescope Polarimeter experiment. The overlap between the two surveys in this work is 33.8 square degrees, before masking bright stars. The signal-to-noise ratio of individual cluster lensing measurements ranges from 2.2 to 8.7, with a total of 11.1 for the stacked cluster weak-lensing signal. We fit for an average weak-lensing mass distribution using three different profiles, a Navarro-Frenk-White profile, a dark-matter-only emulated profile, and a full cosmological hydrodynamic emulated profile. We interpret the differences among the masses inferred by these models as a systematic error of 10\%, which is currently smaller than the statistical error. We obtain the ratio of the SZ-estimated mass to the lensing-estimated mass (the so-called hydrostatic mass bias $1-b$) of $0.74^{+0.13}_{-0.12}$, which is comparable to previous SZ-selected clusters from the Atacama Cosmology Telescope and from the {\sl Planck} Satellite. We conclude with a discussion of the implications for cosmological parameters inferred from cluster abundances compared to cosmic microwave background primary anisotropy measurements.
△ Less
Submitted 16 April, 2018;
originally announced April 2018.
-
Line-Intensity Mapping: 2017 Status Report
Authors:
Ely D. Kovetz,
Marco P. Viero,
Adam Lidz,
Laura Newburgh,
Mubdi Rahman,
Eric Switzer,
Marc Kamionkowski,
James Aguirre,
Marcelo Alvarez,
James Bock,
J. Richard Bond,
Goeffry Bower,
C. Matt Bradford,
Patrick C. Breysse,
Philip Bull,
Tzu-Ching Chang,
Yun-Ting Cheng,
Dongwoo Chung,
Kieran Cleary,
Asantha Corray,
Abigail Crites,
Rupert Croft,
Olivier Doré,
Michael Eastwood,
Andrea Ferrara
, et al. (23 additional authors not shown)
Abstract:
Following the first two annual intensity mapping workshops at Stanford in March 2016 and Johns Hopkins in June 2017, we report on the recent advances in theory, instrumentation and observation that were presented in these meetings and some of the opportunities and challenges that were identified looking forward. With preliminary detections of CO, [CII], Lya and low-redshift 21cm, and a host of exp…
▽ More
Following the first two annual intensity mapping workshops at Stanford in March 2016 and Johns Hopkins in June 2017, we report on the recent advances in theory, instrumentation and observation that were presented in these meetings and some of the opportunities and challenges that were identified looking forward. With preliminary detections of CO, [CII], Lya and low-redshift 21cm, and a host of experiments set to go online in the next few years, the field is rapidly progressing on all fronts, with great anticipation for a flood of new exciting results. This current snapshot provides an efficient reference for experts in related fields and a useful resource for nonspecialists. We begin by introducing the concept of line-intensity mapping and then discuss the broad array of science goals that will be enabled, ranging from the history of star formation, reionization and galaxy evolution to measuring baryon acoustic oscillations at high redshift and constraining theories of dark matter, modified gravity and dark energy. After reviewing the first detections reported to date, we survey the experimental landscape, presenting the parameters and capabilities of relevant instruments such as COMAP, mmIMe, AIM-CO, CCAT-p, TIME, CONCERTO, CHIME, HIRAX, HERA, STARFIRE, MeerKAT/SKA and SPHEREx. Finally, we describe recent theoretical advances: different approaches to modeling line luminosity functions, several techniques to separate the desired signal from foregrounds, statistical methods to analyze the data, and frameworks to generate realistic intensity map simulations.
△ Less
Submitted 26 September, 2017;
originally announced September 2017.
-
Testing Inflation with Large Scale Structure: Connecting Hopes with Reality
Authors:
Marcelo Alvarez,
Tobias Baldauf,
J. Richard Bond,
Neal Dalal,
Roland de Putter,
Olivier Doré,
Daniel Green,
Chris Hirata,
Zhiqi Huang,
Dragan Huterer,
Donghui Jeong,
Matthew C. Johnson,
Elisabeth Krause,
Marilena Loverde,
Joel Meyers,
P. Daniel Meerburg,
Leonardo Senatore,
Sarah Shandera,
Eva Silverstein,
Anže Slosar,
Kendrick Smith,
Matias Zaldarriaga,
Valentin Assassi,
Jonathan Braden,
Amir Hajian
, et al. (3 additional authors not shown)
Abstract:
The statistics of primordial curvature fluctuations are our window into the period of inflation, where these fluctuations were generated. To date, the cosmic microwave background has been the dominant source of information about these perturbations. Large scale structure is however from where drastic improvements should originate. In this paper, we explain the theoretical motivations for pursuing…
▽ More
The statistics of primordial curvature fluctuations are our window into the period of inflation, where these fluctuations were generated. To date, the cosmic microwave background has been the dominant source of information about these perturbations. Large scale structure is however from where drastic improvements should originate. In this paper, we explain the theoretical motivations for pursuing such measurements and the challenges that lie ahead. In particular, we discuss and identify theoretical targets regarding the measurement of primordial non-Gaussianity. We argue that when quantified in terms of the local (equilateral) template amplitude $f_{\rm NL}^{\rm loc}$ ($f_{\rm NL}^{\rm eq}$), natural target levels of sensitivity are $Δf_{\rm NL}^{\rm loc, eq.} \simeq 1$. We highlight that such levels are within reach of future surveys by measuring 2-, 3- and 4-point statistics of the galaxy spatial distribution. This paper summarizes a workshop held at CITA (University of Toronto) on October 23-24, 2014.
△ Less
Submitted 15 December, 2014;
originally announced December 2014.