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

Franz Elsner

Franz Elsner 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.IMhep-th

Catalog footprint

What is connected

15works
3topics
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

15 published item(s)

preprint2020arXiv

Cosmology Inference from Biased Tracers using the EFT-based Likelihood

The effective-field-theory (EFT) approach to the clustering of galaxies and other biased tracers allows for an isolation of the cosmological information that is protected by symmetries, in particular the equivalence principle, and thus is robust to the complicated dynamics of dark matter, gas, and stars on small scales. All existing implementations proceed by making predictions for the lowest-order $n$-point functions of biased tracers, as well as their covariance, and comparing with measurements. Recently, we presented an EFT-based expression for the conditional probability of the density field of a biased tracer given the matter density field, which in principle combines information from arbitrarily high order $n$-point functions. Here, we report results based on this likelihood by applying it to halo catalogs in real space, specifically on the inference of the power spectrum normalization $σ_8$. We include bias terms up to second order as well as the leading higher-derivative term. For a cutoff value of $Λ= 0.1 h\,{\rm Mpc}^{-1}$, we recover the ground-truth value of $σ_8$ to within 95% CL for different halo samples and redshifts. We discuss possible sources for the remaining systematic bias in $σ_8$ as well as future developments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Unbiased pseudo-Cl power spectrum estimation with mode projection

With the steadily improving sensitivity afforded by current and future galaxy surveys, a robust extraction of two-point correlation function measurements may become increasingly hampered by the presence of astrophysical foregrounds or observational systematics. The concept of mode projection has been introduced as a means to remove contaminants for which it is possible to construct a spatial map reflecting the expected signal contribution. Owing to its computational efficiency compared to minimum-variance methods, the sub-optimal pseudo-Cl (PCL) power spectrum estimator is a popular tool for the analysis of high-resolution data sets. Here, we integrate mode projection into the framework of PCL power spectrum estimation. In contrast to results obtained with optimal estimators, we show that the uncorrected projection of template maps leads to biased power spectra. Based on analytical calculations, we find exact closed-form expressions for the expectation value of the bias and demonstrate that they can be recast in a form that allows a numerically efficient evaluation, preserving the favorable O(l_max^3) time complexity of PCL estimator algorithms. Using simulated data sets, we assess the scaling of the bias with various analysis parameters and demonstrate that it can be reliably removed. We conclude that in combination with mode projection, PCL estimators allow for a fast and robust computation of power spectra in the presence of systematic effects -- properties in high demand for the analysis of ongoing and future large scale structure surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Forecasting constraints from the cosmic microwave background on eternal inflation

We forecast the ability of cosmic microwave background (CMB) temperature and polarization datasets to constrain theories of eternal inflation using cosmic bubble collisions. Using the Fisher matrix formalism, we determine both the overall detectability of bubble collisions and the constraints achievable on the fundamental parameters describing the underlying theory. The CMB signatures considered are based on state-of-the-art numerical relativistic simulations of the bubble collision spacetime, evolved using the full temperature and polarization transfer functions. Comparing a theoretical cosmic-variance-limited experiment to the WMAP and Planck satellites, we find that there is no improvement to be gained from future temperature data, that adding polarization improves detectability by approximately 30%, and that cosmic-variance-limited polarization data offer only marginal improvements over Planck. The fundamental parameter constraints achievable depend on the precise values of the tensor-to-scalar ratio and energy density in (negative) spatial curvature. For a tensor-to-scalar ratio of $0.1$ and spatial curvature at the level of $10^{-4}$, using cosmic-variance-limited data it is possible to measure the width of the potential barrier separating the inflating false vacuum from the true vacuum down to $M_{\rm Pl}/500$, and the initial proper distance between colliding bubbles to a factor $π/2$ of the false vacuum horizon size (at three sigma). We conclude that very near-future data will have the final word on bubble collisions in the CMB.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Unbiased methods for removing systematics from galaxy clustering measurements

Measuring the angular clustering of galaxies as a function of redshift is a powerful method for extracting information from the three-dimensional galaxy distribution. The precision of such measurements will dramatically increase with ongoing and future wide-field galaxy surveys. However, these are also increasingly sensitive to observational and astrophysical contaminants. Here, we study the statistical properties of three methods proposed for controlling such systematics -- template subtraction, basic mode projection, and extended mode projection -- all of which make use of externally supplied template maps, designed to characterise and capture the spatial variations of potential systematic effects. Based on a detailed mathematical analysis, and in agreement with simulations, we find that the template subtraction method in its original formulation returns biased estimates of the galaxy angular clustering. We derive closed-form expressions that should be used to correct results for this shortcoming. Turning to the basic mode projection algorithm, we prove it to be free of any bias, whereas we conclude that results computed with extended mode projection are biased. Within a simplified setup, we derive analytical expressions for the bias and discuss the options for correcting it in more realistic configurations. Common to all three methods is an increased estimator variance induced by the cleaning process, albeit at different levels. These results enable unbiased high-precision clustering measurements in the presence of spatially-varying systematics, an essential step towards realising the full potential of current and planned galaxy surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Using hybrid GPU/CPU kernel splitting to accelerate spherical convolutions

We present a general method for accelerating by more than an order of magnitude the convolution of pixelated functions on the sphere with a radially-symmetric kernel. Our method splits the kernel into a compact real-space component and a compact spherical harmonic space component. These components can then be convolved in parallel using an inexpensive commodity GPU and a CPU. We provide models for the computational cost of both real-space and Fourier space convolutions and an estimate for the approximation error. Using these models we can determine the optimum split that minimizes the wall clock time for the convolution while satisfying the desired error bounds. We apply this technique to the problem of simulating a cosmic microwave background (CMB) anisotropy sky map at the resolution typical of the high resolution maps produced by the Planck mission. For the main Planck CMB science channels we achieve a speedup of over a factor of ten, assuming an acceptable fractional rms error of order 1.e-5 in the power spectrum of the output map.

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

Fast Wiener filtering of CMB maps

We present the application of a new method to compute the Wiener filter solution of large and complex data sets. Contrary to the iterative solvers usually employed in signal processing, our algorithm does not require the use of preconditioners to be computationally efficient. The new scheme is conceptually very simple and therefore easy to implement, numerically absolutely stable, and guaranteed to converge. We introduce a messenger field to mediate between the different preferred bases in which signal and noise properties can be specified most conveniently, and rephrase the signal reconstruction problem in terms of this auxiliary variable. We demonstrate the capabilities of the algorithm by applying it to cosmic microwave background (CMB) radiation data obtained by the WMAP satellite.

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
preprint2012arXiv

Accelerating convolutions on the sphere with hybrid GPU/CPU kernel splitting

We present a general method for accelerating by more than an order of magnitude the convolution of pixelated function on the sphere with a radially-symmetric kernel. Our method splits the kernel into a compact real-space, and a compact spherical harmonic space component that can then be convolved in parallel using an inexpensive commodity GPU and a CPU, respectively. We provide models for the computational cost of both real-space and Fourier space convolutions and an estimate for the approximation error. Using these models we can determine the optimum split that minimizes the wall clock time for the convolution while satisfying the desired error bounds. We apply this technique to the problem of simulating a cosmic microwave background sky map at the resolution typical of the high resolution maps of the cosmic microwave background anisotropies produced by the Planck space craft. For the main Planck CMB science channels we achieve a speedup of over a factor of ten, assuming an acceptable fractional rms error of order 10^-5 in the (power spectrum of the) output map.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Efficient Wiener filtering without preconditioning

We present a new approach to calculate the Wiener filter solution of general data sets. It is trivial to implement, flexible, numerically absolutely stable, and guaranteed to converge. Most importantly, it does not require an ingenious choice of preconditioner to work well. The method is capable of taking into account inhomogeneous noise distributions and arbitrary mask geometries. It iteratively builds up the signal reconstruction by means of a messenger field, introduced to mediate between the different preferred bases in which signal and noise properties can be specified most conveniently. Using cosmic microwave background (CMB) radiation data as a showcase, we demonstrate the capabilities of our scheme by computing Wiener filtered WMAP7 temperature and polarization maps at full resolution for the first time. We show how the algorithm can be modified to synthesize fluctuation maps, which, combined with the Wiener filter solution, result in unbiased constrained signal realizations, consistent with the observations. The algorithm performs well even on simulated CMB maps with Planck resolution and dynamic range.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Fast calculation of the Fisher matrix for cosmic microwave background experiments

The Fisher information matrix of the cosmic microwave background (CMB) radiation power spectrum coefficients is a fundamental quantity that specifies the information content of a CMB experiment. In the most general case, its exact calculation scales with the third power of the number of data points N and is therefore computationally prohibitive for state-of-the-art surveys. Applicable to a very large class of CMB experiments without special symmetries, we show how to compute the Fisher matrix in only O(N^2 log N) operations as long as the inverse noise covariance matrix can be applied to a data vector in time O(l_max^3 log l_max). This assumption is true to a good approximation for all CMB data sets taken so far. The method takes into account common systematics such as arbitrary sky coverage and realistic noise correlations. As a consequence, optimal quadratic power spectrum estimation also becomes feasible in O(N^2 log N) operations for this large group of experiments. We discuss the relevance of our findings to other areas of cosmology where optimal power spectrum estimation plays a role.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Likelihood, Fisher information, and systematics of cosmic microwave background experiments

Every experiment is affected by systematic effects that hamper the data analysis and have the potential to ultimately degrade its performance. In the case of probes of the cosmic microwave background (CMB) radiation, a minimal set of issues to consider includes asymmetric beam functions, correlated noise, and incomplete sky coverage. Presuming a simplified scanning strategy that allows for an exact analytical treatment of the problem, we study the impact of systematic effects on the likelihood function of the CMB power spectrum. We use the Fisher matrix, a measure of the information content of a data set, for a quantitative comparison of different experimental configurations. In addition, for various power spectrum coefficients, we explore the functional form of the likelihood directly, and obtain the following results: The likelihood function can deviate systematically from a Gaussian distribution up to the highest multipole values considered in our analysis. Treated exactly, realistic levels of asymmetric beam functions and correlated noise do not by themselves decrease the information yield of CMB experiments nor do they induce noticeable coupling between multipoles. Masking large fractions of the sky, on the other hand, results in a considerably more complex correlation structure of the likelihood function. Combining adjacent power spectrum coefficients into bins can partially mitigate these problems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Local non-Gaussianity in the Cosmic Microwave Background the Bayesian way

We introduce an exact Bayesian approach to search for non-Gaussianity of local type in Cosmic Microwave Background (CMB) radiation data. Using simulated CMB temperature maps, the newly developed technique is compared against the conventional frequentist bispectrum estimator. Starting from the joint probability distribution, we obtain analytic expressions for the conditional probabilities of the primordial perturbations given the data, and for the level of non-Gaussianity, f_nl, given the data and the perturbations. We propose Hamiltonian Monte Carlo sampling as a means to derive realizations of the primordial fluctuations from which we in turn sample f_nl. Although being computationally expensive, this approach allows us to exactly construct the full target posterior probability distribution. When compared to the frequentist estimator, applying the Bayesian method to Gaussian CMB maps provides consistent results. For the analysis of non-Gaussian maps, however, the error bars on f_nl do not show excess variance within the Bayesian framework. This finding is of particular relevance in the light of upcoming high precision CMB measurements obtained by the Planck satellite mission.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Probing local non-Gaussianities within a Bayesian framework

Aims: We outline the Bayesian approach to inferring f_NL, the level of non-Gaussianity of local type. Phrasing f_NL inference in a Bayesian framework takes advantage of existing techniques to account for instrumental effects and foreground contamination in CMB data and takes into account uncertainties in the cosmological parameters in an unambiguous way. Methods: We derive closed form expressions for the joint posterior of f_NL and the reconstructed underlying curvature perturbation, Phi, and deduce the conditional probability densities for f_NL and Phi. Completing the inference problem amounts to finding the marginal density for f_NL. For realistic data sets the necessary integrations are intractable. We propose an exact Hamiltonian sampling algorithm to generate correlated samples from the f_NL posterior. For sufficiently high signal-to-noise ratios, we can exploit the assumption of weak non-Gaussianity to find a direct Monte Carlo technique to generate independent samples from the posterior distribution for f_NL. We illustrate our approach using a simplified toy model of CMB data for the simple case of a 1-D sky. Results: When applied to our toy problem, we find that, in the limit of high signal-to-noise, the sampling efficiency of the approximate algorithm outperforms that of Hamiltonian sampling by two orders of magnitude. When f_NL is not significantly constrained by the data, the more efficient, approximate algorithm biases the posterior density towards f_NL = 0.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The primordial non-Gaussianity of local type (f_NL) in the WMAP 5-year data: the length distribution of CMB skeleton

We present skeleton studies of non-Gaussianity in the CMB temperature anisotropy observed in the WMAP5 data. The local skeleton is traced on the 2D sphere by cubic spline interpolation which leads to more accurate estimation of the intersection positions between the skeleton and the secondary pixels than conventional linear interpolation. We demonstrate that the skeleton-based estimator of non-Gaussianity of the local type (f_NL) - the departure of the length distribution from the corresponding Gaussian expectation - yields an unbiased and sufficiently converged f_NL-likelihood. We analyse the skeleton statistics in the WMAP5 combined V- and W-band data outside the Galactic base-mask determined from the KQ75 sky-coverage. The results are consistent with Gaussian simulations of the the best-fitting cosmological model, but deviate from the previous results determined using the WMAP1 data. We show that it is unlikely that the improved skeleton tracing method, the omission of Q-band data, the modification of the foreground-template fitting method or the absence of 6 extended regions in the new mask contribute to such a deviation. However, the application of the Kp0 base-mask in data processing does improve the consistency with the WMAP1 results. The f_NL-likelihoods of the data are estimated at 9 different smoothing levels. It is unexpected that the best-fit values show positive correlation with the smoothing scales. Further investigation argues against a point-source or goodness-of-fit explanation but finds that about 30% of either Gaussian or f_NL samples having better goodness-of-fit than the WMAP5 show a similar correlation. We present the estimate f_NL=47.3+/-34.9 (1sigma error) determined from the first four smoothing angles and f_NL=76.8+/-43.1 for the combination of all nine. The former result may be overestimated at the 0.21sigma-level because of point sources.

0 citations0 reviewsNo code linkedOpen access