BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Marcel Schmittfull

Marcel Schmittfull contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.COastro-ph.IMhep-phhep-thphysics.data-an

Trust snapshot

Quick read

Trust 21 - Emerging
10works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 published item(s)

preprint2022arXiv

Iterative reconstruction excursions for Baryon Acoustic Oscillations and beyond

The density field reconstruction technique has been widely used for recovering the Baryon Acoustic Oscillation (BAO) feature in galaxy surveys that has been degraded due to nonlinearities. Recent studies advocated adopting iterative steps to improve the recovery much beyond that of the standard technique. In this paper, we investigate the performance of a few selected iterative reconstruction techniques focusing on the BAO and the broadband-shape of the two-point clustering. We include redshift-space distortions, halo bias, and shot noise and inspect the components of the reconstructed field in Fourier space and in configuration space using both density field-based reconstruction and displacement field-based reconstruction. We find that the displacement field reconstruction becomes quickly challenging in the presence of non-negligible shot noise and therefore present surrogate methods that can be practically applied to a much more sparse field such as galaxies. For a galaxy field, implementing a debiasing step to remove the Lagrangian bias appears crucial for the displacement field reconstruction. We show that the iterative reconstruction does not substantially improve the BAO feature beyond an aggressively optimized standard reconstruction with a small smoothing kernel. However, we find taking iterative steps allows us to use a small smoothing kernel more `stably', i.e., without causing a substantial deviation from the linear power spectrum on large scales. In one specific example we studied, we find that a deviation of 13\% in $P( k \sim 0.1h/Mpc)$ with an aggressive standard reconstruction can reduce to 3-4\% with iterative steps.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Fewer Mocks and Less Noise: Reducing the Dimensionality of Cosmological Observables with Subspace Projections

Creating accurate and low-noise covariance matrices represents a formidable challenge in modern-day cosmology. We present a formalism to compress arbitrary observables into a small number of bins by projection into a model-specific subspace that minimizes the prior-averaged log-likelihood error. The lower dimensionality leads to a dramatic reduction in covariance matrix noise, significantly reducing the number of mocks that need to be computed. Given a theory model, a set of priors, and a simple model of the covariance, our method works by using singular value decompositions to construct a basis for the observable that is close to Euclidean; by restricting to the first few basis vectors, we can capture almost all the constraining power in a lower-dimensional subspace. Unlike conventional approaches, the method can be tailored for specific analyses and captures non-linearities that are not present in the Fisher matrix, ensuring that the full likelihood can be reproduced. The procedure is validated with full-shape analyses of power spectra from BOSS DR12 mock catalogs, showing that the 96-bin power spectra can be replaced by 12 subspace coefficients without biasing the output cosmology; this allows for accurate parameter inference using only $\sim 100$ mocks. Such decompositions facilitate accurate testing of power spectrum covariances; for the largest BOSS data chunk, we find that: (a) analytic covariances provide accurate models (with or without trispectrum terms); and (b) using the sample covariance from the MultiDark-Patchy mocks incurs a $\sim 0.5σ$ shift in $Ω_m$, unless the subspace projection is applied. The method is easily extended to higher order statistics; the $\sim 2000$-bin bispectrum can be compressed into only $\sim 10$ coefficients, allowing for accurate analyses using few mocks and without having to increase the bin sizes.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Capturing Non-Gaussianity of the Large-Scale Structure with Weighted Skew-Spectra

The forthcoming generation of wide-field galaxy surveys will probe larger volumes and galaxy densities, thus allowing for a much larger signal-to-noise ratio for higher-order clustering statistics, in particular the galaxy bispectrum. Extracting this information, however, is more challenging than using the power spectrum due to more complex theoretical modeling, as well as significant computational cost of evaluating the bispectrum signal and the error budget. To overcome these challenges, several proxy statistics have been proposed in the literature, which partially or fully capture the information in the bispectrum, while being computationally less expensive than the bispectrum. One such statistics are {\it weighted skew-spectra}, which are cross-spectra of the density field and appropriately weighted quadratic fields. Using Fisher forecasts, we show that the information in these skew-spectra is equivalent to that in the bispectrum for parameters that appear as amplitudes in the bispectrum model, such as galaxy bias parameters or the amplitude of primordial non-Gaussianity. We consider three shapes of the primordial bispectrum: local, equilateral and that due to massive particles with spin two during inflation. To obtain constraints that match those from a measurement of the full bispectrum, we find that it is crucial to account for the full covariance matrix of the skew-spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A bias to CMB lensing measurements from the bispectrum of large-scale structure

The rapidly improving precision of measurements of gravitational lensing of the Cosmic Microwave Background (CMB) also requires a corresponding increase in the precision of theoretical modeling. A commonly made approximation is to model the CMB deflection angle or lensing potential as a Gaussian random field. In this paper, however, we analytically quantify the influence of the non-Gaussianity of large-scale structure lenses, arising from nonlinear structure formation, on CMB lensing measurements. In particular, evaluating the impact of the non-zero bispectrum of large-scale structure on the relevant CMB four-point correlation functions, we find that there is a bias to estimates of the CMB lensing power spectrum. For temperature-based lensing reconstruction with CMB Stage-III and Stage-IV experiments, we find that this lensing power spectrum bias is negative and is of order one percent of the signal. This corresponds to a shift of multiple standard deviations for these upcoming experiments. We caution, however, that our numerical calculation only evaluates two of the largest bias terms and thus only provides an approximate estimate of the full bias. We conclude that further investigation into lensing biases from nonlinear structure formation is required and that these biases should be accounted for in future lensing analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Fast Large Scale Structure Perturbation Theory using 1D FFTs

The usual fluid equations describing the large-scale evolution of mass density in the universe can be written as local in the density, velocity divergence, and velocity potential fields. As a result, the perturbative expansion in small density fluctuations, usually written in terms of convolutions in Fourier space, can be written as a series of products of these fields evaluated at the same location in configuration space. Based on this, we establish a new method to numerically evaluate the 1-loop power spectrum (i.e., Fourier transform of the 2-point correlation function) with one-dimensional Fast Fourier Transforms. This is exact and a few orders of magnitude faster than previously used numerical approaches. Numerical results of the new method are in excellent agreement with the standard quadrature integration method. This fast model evaluation can in principle be extended to higher loop order where existing codes become painfully slow. Our approach follows by writing higher order corrections to the 2-point correlation function as, e.g., the correlation between two second-order fields or the correlation between a linear and a third-order field. These are then decomposed into products of correlations of linear fields and derivatives of linear fields. The method can also be viewed as evaluating three-dimensional Fourier space convolutions using products in configuration space, which may also be useful in other contexts where similar integrals appear.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

FFT-PT: Reducing the two-loop large-scale structure power spectrum to low-dimensional radial integrals

Modeling the large-scale structure of the universe on nonlinear scales has the potential to substantially increase the science return of upcoming surveys by increasing the number of modes available for model comparisons. One way to achieve this is to model nonlinear scales perturbatively. Unfortunately, this involves high-dimensional loop integrals that are cumbersome to evaluate. Trying to simplify this, we show how two-loop (next-to-next-to-leading order) corrections to the density power spectrum can be reduced to low-dimensional, radial integrals. Many of those can be evaluated with a one-dimensional Fast Fourier Transform, which is significantly faster than the five-dimensional Monte-Carlo integrals that are needed otherwise. The general idea of this FFT-PT method is to switch between Fourier and position space to avoid convolutions and integrate over orientations, leaving only radial integrals. This reformulation is independent of the underlying shape of the initial linear density power spectrum and should easily accommodate features such as those from baryonic acoustic oscillations. We also discuss how to account for halo bias and redshift space distortions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Improving photometric redshifts with Ly$α$ tomography

Forming a three dimensional view of the Universe is a long-standing goal of astronomical observations, and one that becomes increasingly difficult at high redshift. In this paper we discuss how tomography of the intergalactic medium (IGM) at $z\simeq 2.5$ can be used to estimate the redshifts of massive galaxies in a large volume of the Universe based on spectra of galaxies in their background. Our method is based on the fact that hierarchical structure formation leads to a strong dependence of the halo density on large-scale environment. A map of the latter can thus be used to refine our knowledge of the redshifts of halos and the galaxies and AGN which they host. We show that tomographic maps of the IGM at a resolution of $2.5\,h^{-1}$Mpc can determine the redshifts of more than 90 per cent of massive galaxies with redshift uncertainty $Δz/(1+z)=0.01$. Higher resolution maps allow such redshift estimation for lower mass galaxies and halos.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Matter bispectrum of large-scale structure: Three-dimensional comparison between theoretical models and numerical simulations

We study the matter bispectrum of the large-scale structure by comparing different perturbative and phenomenological models with measurements from $N$-body simulations obtained with a modal bispectrum estimator. Using shape and amplitude correlators, we directly compare simulated data with theoretical models over the full three-dimensional domain of the bispectrum, for different redshifts and scales. We review and investigate the main perturbative methods in the literature that predict the one-loop bispectrum: standard perturbation theory, effective field theory, resummed Lagrangian and renormalised perturbation theory, calculating the latter also at two loops for some triangle configurations. We find that effective field theory (EFT) succeeds in extending the range of validity furthest into the mildly nonlinear regime, albeit at the price of free extra parameters requiring calibration on simulations. For the more phenomenological halo model, we confirm that despite its validity in the deeply nonlinear regime it has a deficit of power on intermediate scales, which worsens at higher redshifts; this issue is ameliorated, but not solved, by combined halo-perturbative models. We show from simulations that in this transition region there is a strong squeezed bispectrum component that is significantly underestimated in the halo model at earlier redshifts. We thus propose a phenomenological method for alleviating this deficit, which we develop into a simple phenomenological "three-shape" benchmark model based on the three fundamental shapes we have obtained from studying the halo model. When calibrated on the simulations, this three-shape benchmark model accurately describes the bispectrum on all scales and redshifts considered, providing a prototype bispectrum HALOFIT-like methodology that could be used to describe and test parameter dependencies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Delensing the CMB with the Cosmic Infrared Background

As confusion with lensing B-modes begins to limit experiments that search for primordial B-mode polarization, robust methods for delensing the CMB polarization sky are becoming increasingly important. We investigate in detail the possibility of delensing the CMB with the cosmic infrared background (CIB), emission from dusty star-forming galaxies that is an excellent tracer of the CMB lensing signal, in order to improve constraints on the tensor-to-scalar ratio $r$. We find that the maps of the CIB, such as current Planck satellite maps at 545 GHz, can be used to remove more than half of the lensing B-mode power. Calculating optimal combinations of different large-scale-structure tracers for delensing, we find that co-adding CIB data and external arcminute-resolution CMB lensing reconstruction can lead to significant additional improvements in delensing performance. We investigate whether measurement uncertainty in the CIB spectra will degrade the delensing performance if no model of the CIB spectra is assumed, and instead the CIB spectra are marginalized over, when constraining $r$. We find that such uncertainty does not significantly affect B-mode surveys smaller than a few thousand degrees. Even for larger surveys it causes only a moderate reduction in CIB delensing performance, especially if the surveys have high (arcminute) resolution, which allows self-calibration of the delensing procedure. Though further work on the impact of foreground residuals is required, our overall conclusions for delensing with current CIB data are optimistic: this delensing method can tighten constraints on $r$ by a factor up to $\approx2.2$, and by a factor up to $\approx4$ when combined with external $\approx 3 μ$K-arcmin lensing reconstruction, without requiring the modeling of CIB properties. CIB delensing is thus a promising method for the upcoming generation of CMB polarization surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Near optimal bispectrum estimators for large-scale structure

Clustering of large-scale structure provides significant cosmological information through the power spectrum of density perturbations. Additional information can be gained from higher-order statistics like the bispectrum, especially to break the degeneracy between the linear halo bias $b_1$ and the amplitude of fluctuations $σ_8$. We propose new simple, computationally inexpensive bispectrum statistics that are near optimal for the specific applications like bias determination. Corresponding to the Legendre decomposition of nonlinear halo bias and gravitational coupling at second order, these statistics are given by the cross-spectra of the density with three quadratic fields: the squared density, a tidal term, and a shift term. For halos and galaxies the first two have associated nonlinear bias terms $b_2$ and $b_{s^2}$, respectively, while the shift term has none in the absence of velocity bias (valid in the $k \rightarrow 0$ limit). Thus the linear bias $b_1$ is best determined by the shift cross-spectrum, while the squared density and tidal cross-spectra mostly tighten constraints on $b_2$ and $b_{s^2}$ once $b_1$ is known. Since the form of the cross-spectra is derived from optimal maximum-likelihood estimation, they contain the full bispectrum information on bias parameters. Perturbative analytical predictions for their expectation values and covariances agree with simulations on large scales, $k\lesssim 0.09h/\mathrm{Mpc}$ at $z=0.55$ with Gaussian $R=20h^{-1}\mathrm{Mpc}$ smoothing, for matter-matter-matter, and matter-matter-halo combinations. For halo-halo-halo cross-spectra the model also needs to include corrections to the Poisson stochasticity.

0 citations0 reviewsNo code linkedOpen access