BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Benjamin Wandelt

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

ResearcherUnclaimed source record
astro-ph.COastro-ph.IMApplicationsastro-phastro-ph.GAComputer Visionhep-phMachine Learning

Catalog footprint

What is connected

22works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

22 published item(s)

preprint2022arXiv

Simulation-Based Inference of Strong Gravitational Lensing Parameters

In the coming years, a new generation of sky surveys, in particular, Euclid Space Telescope (2022), and the Rubin Observatory's Legacy Survey of Space and Time (LSST, 2023) will discover more than 200,000 new strong gravitational lenses, which represents an increase of more than two orders of magnitude compared to currently known sample sizes. Accurate and fast analysis of such large volumes of data under a statistical framework is therefore crucial for all sciences enabled by strong lensing. Here, we report on the application of simulation-based inference methods, in particular, density estimation techniques, to the predictions of the set of parameters of strong lensing systems from neural networks. This allows us to explicitly impose desired priors on lensing parameters, while guaranteeing convergence to the optimal posterior in the limit of perfect performance.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The CAMELS Multifield Dataset: Learning the Universe's Fundamental Parameters with Artificial Intelligence

We present the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) Multifield Dataset, CMD, a collection of hundreds of thousands of 2D maps and 3D grids containing many different properties of cosmic gas, dark matter, and stars from 2,000 distinct simulated universes at several cosmic times. The 2D maps and 3D grids represent cosmic regions that span $\sim$100 million light years and have been generated from thousands of state-of-the-art hydrodynamic and gravity-only N-body simulations from the CAMELS project. Designed to train machine learning models, CMD is the largest dataset of its kind containing more than 70 Terabytes of data. In this paper we describe CMD in detail and outline a few of its applications. We focus our attention on one such task, parameter inference, formulating the problems we face as a challenge to the community. We release all data and provide further technical details at https://camels-multifield-dataset.readthedocs.io.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Role of Machine Learning in the Next Decade of Cosmology

In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

CMB-HD: Astro2020 RFI Response

CMB-HD is a proposed ultra-deep (0.5 uk-arcmin), high-resolution (15 arcseconds) millimeter-wave survey over half the sky that would answer many outstanding questions in both fundamental physics of the Universe and astrophysics. This survey would be delivered in 7.5 years of observing 20,000 square degrees, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bayesian hierarchical modelling of weak lensing - the golden goal

To accomplish correct Bayesian inference from weak lensing shear data requires a complete statistical description of the data. The natural framework to do this is a Bayesian Hierarchical Model, which divides the chain of reasoning into component steps. Starting with a catalogue of shear estimates in tomographic bins, we build a model that allows us to sample simultaneously from the the underlying tomographic shear fields and the relevant power spectra (E-mode, B-mode, and E-B, for auto- and cross-power spectra). The procedure deals easily with masked data and intrinsic alignments. Using Gibbs sampling and messenger fields, we show with simulated data that the large (over 67000-)dimensional parameter space can be efficiently sampled and the full joint posterior probability density function for the parameters can feasibly be obtained. The method correctly recovers the underlying shear fields and all of the power spectra, including at levels well below the shot noise.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Comparing cosmic web classifiers using information theory

We introduce a decision scheme for optimally choosing a classifier, which segments the cosmic web into different structure types (voids, sheets, filaments, and clusters). Our framework, based on information theory, accounts for the design aims of different classes of possible applications: (i) parameter inference, (ii) model selection, and (iii) prediction of new observations. As an illustration, we use cosmographic maps of web-types in the Sloan Digital Sky Survey to assess the relative performance of the classifiers T-web, DIVA and ORIGAMI for: (i) analyzing the morphology of the cosmic web, (ii) discriminating dark energy models, and (iii) predicting galaxy colors. Our study substantiates a data-supported connection between cosmic web analysis and information theory, and paves the path towards principled design of analysis procedures for the next generation of galaxy surveys. We have made the cosmic web maps, galaxy catalog, and analysis scripts used in this work publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Halo detection via large-scale Bayesian inference

We present a proof-of-concept of a novel and fully Bayesian methodology designed to detect halos of different masses in cosmological observations subject to noise and systematic uncertainties. Our methodology combines the previously published Bayesian large-scale structure inference algorithm, HADES, and a Bayesian chain rule (the Blackwell-Rao Estimator), which we use to connect the inferred density field to the properties of dark matter halos. To demonstrate the capability of our approach we construct a realistic galaxy mock catalogue emulating the wide-area 6-degree Field Galaxy Survey, which has a median redshift of approximately 0.05. Application of HADES to the catalogue provides us with accurately inferred three-dimensional density fields and corresponding quantification of uncertainties inherent to any cosmological observation. We then use a cosmological simulation to relate the amplitude of the density field to the probability of detecting a halo with mass above a specified threshold. With this information we can sum over the HADES density field realisations to construct maps of detection probabilities and demonstrate the validity of this approach within our mock scenario. We find that the probability of successful of detection of halos in the mock catalogue increases as a function of the signal-to-noise of the local galaxy observations. Our proposed methodology can easily be extended to account for more complex scientific questions and is a promising novel tool to analyse the cosmic large-scale structure in observations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hierarchical Cosmic Shear Power Spectrum Inference

We develop a Bayesian hierarchical modelling approach for cosmic shear power spectrum inference, jointly sampling from the posterior distribution of the cosmic shear field and its (tomographic) power spectra. Inference of the shear power spectrum is a powerful intermediate product for a cosmic shear analysis, since it requires very few model assumptions and can be used to perform inference on a wide range of cosmological models \emph{a posteriori} without loss of information. We show that joint posterior for the shear map and power spectrum can be sampled effectively by Gibbs sampling, iteratively drawing samples from the map and power spectrum, each conditional on the other. This approach neatly circumvents difficulties associated with complicated survey geometry and masks that plague frequentist power spectrum estimators, since the power spectrum inference provides prior information about the field in masked regions at every sampling step. We demonstrate this approach for inference of tomographic shear $E$-mode, $B$-mode and $EB$-cross power spectra from a simulated galaxy shear catalogue with a number of important features; galaxies distributed on the sky and in redshift with photometric redshift uncertainties, realistic random ellipticity noise for every galaxy and a complicated survey mask. The obtained posterior distributions for the tomographic power spectrum coefficients recover the underlying simulated power spectra for both $E$- and $B$-modes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Bayesian analysis of the dynamic cosmic web in the SDSS galaxy survey

Recent application of the Bayesian algorithm BORG to the Sloan Digital Sky Survey (SDSS) main sample galaxies resulted in the physical inference of the formation history of the observed large-scale structure from its origin to the present epoch. In this work, we use these inferences as inputs for a detailed probabilistic cosmic web-type analysis. To do so, we generate a large set of data-constrained realizations of the large-scale structure using a fast, fully non-linear gravitational model. We then perform a dynamic classification of the cosmic web into four distinct components (voids, sheets, filaments, and clusters) on the basis of the tidal field. Our inference framework automatically and self-consistently propagates typical observational uncertainties to web-type classification. As a result, this study produces accurate cosmographic classification of large-scale structure elements in the SDSS volume. By also providing the history of these structure maps, the approach allows an analysis of the origin and growth of the early traces of the cosmic web present in the initial density field and of the evolution of global quantities such as the volume and mass filling fractions of different structures. For the problem of web-type classification, the results described in this work constitute the first connection between theory and observations at non-linear scales including a physical model of structure formation and the demonstrated capability of uncertainty quantification. A connection between cosmology and information theory using real data also naturally emerges from our probabilistic approach. Our results constitute quantitative chrono-cosmography of the complex web-like patterns underlying the observed galaxy distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Bayesian inference of CMB gravitational lensing

The Planck satellite, along with several ground based telescopes, have mapped the cosmic microwave background (CMB) at sufficient resolution and signal-to-noise so as to allow a detection of the subtle distortions due to the gravitational influence of the intervening matter distribution. A natural modeling approach is to write a Bayesian hierarchical model for the lensed CMB in terms of the unlensed CMB and the lensing potential. So far there has been no feasible algorithm for inferring the posterior distribution of the lensing potential from the lensed CMB map. We propose a solution that allows efficient Markov Chain Monte Carlo sampling from the joint posterior of the lensing potential and the unlensed CMB map using the Hamiltonian Monte Carlo technique. The main conceptual step in the solution is a re-parameterization of CMB lensing in terms of the lensed CMB and the "inverse lensing" potential. We demonstrate a fast implementation on simulated data including noise and a sky cut, that uses a further acceleration based on a very mild approximation of the inverse lensing potential. We find that the resulting Markov Chain has short correlation lengths and excellent convergence properties, making it promising for application to high resolution CMB data sets of the future.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cosmic web-type classification using decision theory

We propose a decision criterion for segmenting the cosmic web into different structure types (voids, sheets, filaments, and clusters) on the basis of their respective probabilities and the strength of data constraints. Our approach is inspired by an analysis of games of chance where the gambler only plays if a positive expected net gain can be achieved based on some degree of privileged information. The result is a general solution for classification problems in the face of uncertainty, including the option of not committing to a class for a candidate object. As an illustration, we produce high-resolution maps of web-type constituents in the nearby Universe as probed by the Sloan Digital Sky Survey main galaxy sample. Other possible applications include the selection and labelling of objects in catalogues derived from astronomical survey data.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dark matter voids in the SDSS galaxy survey

What do we know about voids in the dark matter distribution given the Sloan Digital Sky Survey (SDSS) and assuming the $Λ\mathrm{CDM}$ model? Recent application of the Bayesian inference algorithm BORG to the SDSS Data Release 7 main galaxy sample has generated detailed Eulerian and Lagrangian representations of the large-scale structure as well as the possibility to accurately quantify corresponding uncertainties. Building upon these results, we present constrained catalogs of voids in the Sloan volume, aiming at a physical representation of dark matter underdensities and at the alleviation of the problems due to sparsity and biasing on galaxy void catalogs. To do so, we generate data-constrained reconstructions of the presently observed large-scale structure using a fully non-linear gravitational model. We then find and analyze void candidates using the VIDE toolkit. Our methodology therefore predicts the properties of voids based on fusing prior information from simulations and data constraints. For usual void statistics (number function, ellipticity distribution and radial density profile), all the results obtained are in agreement with dark matter simulations. Our dark matter void candidates probe a deeper void hierarchy than voids directly based on the observed galaxies alone. The use of our catalogs therefore opens the way to high-precision void cosmology at the level of the dark matter field. We will make the void catalogs used in this work available at http://www.cosmicvoids.net.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Joint resonant CMB power spectrum and bispectrum estimation

We develop the tools necessary to assess the statistical significance of resonant features in the CMB correlation functions, combining power spectrum and bispectrum measurements. This significance is typically addressed by running a large number of simulations to derive the probability density function (PDF) of the feature-amplitude in the Gaussian case. Although these simulations are tractable for the power spectrum, for the bispectrum they require significant computational resources. We show that, by assuming that the PDF is given by a multi-variate Gaussian where the covariance is determined by the Fisher matrix of the sine and cosine terms, we can efficiently produce spectra that are statistically close to those derived from full simulations. By drawing a large number of spectra from this PDF, both for the power spectrum and the bispectrum, we can quickly determine the statistical significance of candidate signatures in the CMB, considering both single frequency and multi-frequency estimators. We show that for resonance models, cosmology and foreground parameters have little influence on the estimated amplitude, which allows to simplify the analysis considerably. A more precise likelihood treatment can then be applied to candidate signatures only. We also discuss a modal expansion approach for the power spectrum, aimed at quickly scanning through large families of oscillating models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

One-point remapping of Lagrangian perturbation theory in the mildly non-linear regime of cosmic structure formation

On the smallest scales, three-dimensional large-scale structure surveys contain a wealth of cosmological information which cannot be trivially extracted due to the non-linear dynamical evolution of the density field. Lagrangian perturbation theory (LPT) is widely applied to the generation of mock halo catalogs and data analysis. In this work, we compare topological features of the cosmic web such as voids, sheets, filaments and clusters, in the density fields predicted by LPT and full numerical simulation of gravitational large-scale structure formation. We propose a method designed to improve the correspondence between these density fields, in the mildly non-linear regime. We develop a computationally fast and flexible tool for a variety of cosmological applications. Our method is based on a remapping of the approximately-evolved density field, using information extracted from N-body simulations. The remapping procedure consists of replacing the one-point distribution of the density contrast by one which better accounts for the full gravitational dynamics. As a result, we obtain a physically more pertinent density field on a point-by-point basis, while also improving higher-order statistics predicted by LPT. We quantify the approximation error in the power spectrum and in the bispectrum as a function of scale and redshift. Our remapping procedure improves one-, two- and three-point statistics at scales down to 8 Mpc/h.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Probabilistic cartography of the large-scale structure

The BORG algorithm is an inference engine that derives the initial conditions given a cosmological model and galaxy survey data, and produces physical reconstructions of the underlying large-scale structure by assimilating the data into the model. We present the application of BORG to real galaxy catalogs and describe the primordial and late-time large-scale structure in the considered volumes. We then show how these results can be used for building various probabilistic maps of the large-scale structure, with rigorous propagation of uncertainties. In particular, we study dynamic cosmic web elements and secondary effects in the cosmic microwave background.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Bayesian large-scale structure inference: initial conditions and the cosmic web

We describe an innovative statistical approach for the ab initio simultaneous analysis of the formation history and morphology of the large-scale structure of the inhomogeneous Universe. Our algorithm explores the joint posterior distribution of the many millions of parameters involved via efficient Hamiltonian Markov Chain Monte Carlo sampling. We describe its application to the Sloan Digital Sky Survey data release 7 and an additional non-linear filtering step. We illustrate the use of our findings for cosmic web analysis: identification of structures via tidal shear analysis and inference of dark matter voids.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

CMB Polarization can constrain cosmology better than CMB temperature

We demonstrate that for a cosmic variance limited experiment, CMB E polarization alone places stronger constraints on cosmological parameters than CMB temperature. For example, we show that EE can constrain parameters better than TT by up to a factor 2.8 when a multipole range of l=30-2500 is considered. We expose the physical effects at play behind this remarkable result and study how it depends on the multipole range included in the analysis. In most relevant cases, TE or EE surpass the TT based cosmological constraints. This result is important as the small scale astrophysical foregrounds are expected to have a much reduced impact on polarization, thus opening the possibility of building cleaner and more stringent constraints of the LCDM model. This is relevant specially for proposed future CMB satellite missions, such as CORE or PRISM, that are designed to be cosmic variance limited in polarization till very large multipoles. We perform the same analysis for a Planck-like experiment, and conclude that even in this case TE alone should determine the constraint on $Ω_ch^2$ better than TT by 15%, while determining $Ω_bh^2$, $n_s$ and $θ$ with comparable accuracy. Finally, we explore a few classical extensions of the LCDM model and show again that CMB polarization alone provides more stringent constraints than CMB temperature in case of a cosmic variance limited experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Insights into the content and spatial distribution of dust from the integrated spectral properties of galaxies

[Abridged] We present a new approach to investigate the content and spatial distribution of dust in structurally unresolved star-forming galaxies from the observed dependence of integrated spectral properties on galaxy inclination. We develop an innovative combination of generic models of radiative transfer (RT) in dusty media with a prescription for the spectral evolution of galaxies, via the association of different geometric components of galaxies with stars in different age ranges. We show that a wide range of RT models all predict a quasi-universal relation between slope of the attenuation curve at any wavelength and V-band attenuation optical depth in the diffuse interstellar medium (ISM), at all galaxy inclinations. This relation predicts steeper (shallower) dust attenuation curves than both the Calzetti and MW curves at small (large) attenuation optical depths, which implies that geometry and orientation effects have a stronger influence on the shape of the attenuation curve than changes in the optical properties of dust grains. We use our combined RT and spectral evolution model to interpret the observed dependence of the Hα/Hβ ratio and ugrizYJH attenuation curve on inclination in a sample of ~23 000 nearby star-forming galaxies. From a Bayesian MCMC fit, we measure the central face-on B-band optical depth of this sample to be tau_B\perp~1.8\pm0.2. We also quantify the enhanced optical depth towards newly formed stars in their birth clouds, finding this to be significantly larger in galaxies with bulges than in disc-dominated galaxies, while tau_B\perp is roughly similar in both cases. Finally, we show that neglecting the effect of geometry and orientation on attenuation can severely bias the interpretation of galaxy spectral energy distributions, as the impact on broadband colours can reach up to 0.3-0.4 mag at optical wavelengths and 0.1 mag at near-infrared ones.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Optimal bispectrum estimator and simulations of the the CMB Lensing-ISW non-Gaussian signal

In this paper we present the tools to optimally extract the Lensing-Integrated Sachs Wolfe (L-ISW) bispectrum signal from future CMB data. We implement two different methods to simulate the non-Gaussian CMB maps with the L-ISW signal: a non-perturbative method based on the FLINTS lensing code and the separable mode expansion method. We implement the Komatsu, Spergel and Wandelt (KSW) optimal estimator analysis for the Lensing-ISW bispectrum and we test it on the non-Gaussian simulations in the case of a realistic CMB experimental settings with an inhomogeneous sky coverage. We show that the estimator approaches the Cramer-Rao bound and that Wiener filtering the L-ISW simulations gives a slight improvement on the estimate of $f_{NL}^{L-ISW}$ of $\leq 10%$. For a realistic CMB experimental setting accounting for anisotropic noise and masked sky, we show that the linear term of the estimator is highly correlated to the cubic term and it is necessary to recover the signal and the optimal error bars. We also show that the L-ISW bispectrum, if not correctly accounted for, yields an underestimation of the $f_{NL}^{local}$ error bars of $\simeq 4%$. A joint analysis of the non-Gaussian shapes and/or L-ISW template subtraction is needed in order to recover unbiased results of the primordial non-Gaussian signal from ongoing and future CMB experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Temperature and Polarization CMB Maps from Primordial non-Gaussianities of the Local Type

The forthcoming Planck experiment will provide high sensitivity polarization measurements that will allow us to further tighten the f_NL bounds from the temperature data. Monte Carlo simulations of non-Gaussian CMB maps have been used as a fundamental tool to characterize non-Gaussian signatures in the data, as they allow us to calibrate any statistical estimators and understand the effect of systematics, foregrounds and other contaminants. We describe an algorithm to generate high-angular resolution simulations of non-Gaussian CMB maps in temperature and polarization. We consider non-Gaussianities of the local type, for which the level of non-Gaussianity is defined by the dimensionless parameter, f_NL. We then apply the temperature and polarization fast cubic statistics recently developed by Yadav et al. to a set of non-Gaussian temperature and polarization simulations. We compare our results to theoretical expectations based on a Fisher matrix analysis, test the unbiasedness of the estimator, and study the dependence of the error bars on f_NL. All our results are in very good agreement with theoretical predictions, thus confirming the reliability of both the simulation algorithm and the fast cubic temperature and polarization estimator.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Galactic foreground contribution to the BEAST CMB Anisotropy Maps

We report limits on the Galactic foreground emission contribution to the Background Emission Anisotropy Scanning Telescope (BEAST) Ka- and Q-band CMB anisotropy maps. We estimate the contribution from the cross-correlations between these maps and the foreground emission templates of an H$α$ map, a de-striped version of the Haslam et al. 408 MHz map, and a combined 100 $μ$m IRAS/DIRBE map. Our analysis samples the BEAST $\sim10^\circ$ declination band into 24 one-hour (RA) wide sectors with $\sim7900$ pixels each, where we calculate: (a) the linear correlation coefficient between the anisotropy maps and the templates; (b) the coupling constants between the specific intensity units of the templates and the antenna temperature at the BEAST frequencies and (c) the individual foreground contributions to the BEAST anisotropy maps. The peak sector contributions of the contaminants in the Ka-band are of 56.5% free-free with a coupling constant of $8.3\pm0.4$ $μ$K/R, and 67.4% dust with $45.0\pm2.0$ $μ$K/(MJy/sr). In the Q-band the corresponding values are of 64.4% free-free with $4.1\pm0.2$ $μ$K/R and 67.5% dust with $24.0\pm1.0$ $μ$K/(MJy/sr). Using a lower limit of 10% in the relative uncertainty of the coupling constants, we can constrain the sector contributions of each contaminant in both maps to $< 20$% in 21 (free-free), 19 (dust) and 22 (synchrotron) sectors. At this level, all these sectors are found outside of the $\mid$b$\mid = 14.6^\circ$ region. By performing the same correlation analysis as a function of Galactic scale height, we conclude that the region within $b=\pm17.5^{\circ}$ should be removed from the BEAST maps for CMB studies in order to keep individual Galactic contributions below $\sim 1$% of the map's rms.

0 citations0 reviewsNo code linkedOpen access