Researcher profile

Florian Beutler

Florian Beutler contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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

13 published item(s)

preprint2022arXiv

$\texttt{matryoshka}$: Halo Model Emulator for the Galaxy Power Spectrum

We present $\texttt{matryoshka}$, a suite of neural network based emulators and accompanying Python package that have been developed with the goal of producing fast and accurate predictions of the nonlinear galaxy power spectrum. The suite of emulators consists of four linear component emulators, from which fast linear predictions of the power spectrum can be made, allowing all nonlinearities to be included in predictions from a nonlinear boost component emulator. The linear component emulators includes an emulator for the matter transfer function that produces predictions in $\sim 0.0004 \ \mathrm{s}$, with an error of $<0.08\%$ (at $1σ$ level) on scales $10^{-4} \ h \ \mathrm{Mpc}^{-1}<k<10^1 \ h \ \mathrm{Mpc}^{-1}$. In this paper we demonstrate $\texttt{matryoshka}$ by training the nonlinear boost component emulator with analytic training data calculated with HALOFIT, that has been designed to replicate training data that would be generated using numerical simulations. Combining all the component emulator predictions we achieve an accuracy of $< 0.75\%$ (at $1σ$ level) when predicting the real space nonlinear galaxy power spectrum on scales $0.0025 \ h \ \mathrm{Mpc}^{-1}<k<1 \ h \ \mathrm{Mpc}^{-1}$. We use $\texttt{matryoshka}$ to investigate the impact of the analysis setup on cosmological constraints by conducting several full shape analyses of the real space galaxy power spectrum. Specifically we investigate the impact of the minimum scale (or $k_\mathrm{max}$), finding an improvement of $\sim 1.8\times$ in the constraint on $σ_8$ by pushing $k_\mathrm{max}$ from $k_\mathrm{max}=0.25 \ h \ \mathrm{Mpc}^{-1}$ to $k_\mathrm{max}=0.85 \ h \ \mathrm{Mpc}^{-1}$, highlighting the potential gains when using clustering emulators such as $\texttt{matryoshka}$ in cosmological analyses.

preprint2022arXiv

Accelerating astronomical and cosmological inference with Preconditioned Monte Carlo

We introduce Preconditioned Monte Carlo (PMC), a novel Monte Carlo method for Bayesian inference that facilitates efficient sampling of probability distributions with non-trivial geometry. PMC utilises a Normalising Flow (NF) in order to decorrelate the parameters of the distribution and then proceeds by sampling from the preconditioned target distribution using an adaptive Sequential Monte Carlo (SMC) scheme. The results produced by PMC include samples from the posterior distribution and an estimate of the model evidence that can be used for parameter inference and model comparison respectively. The aforementioned framework has been thoroughly tested in a variety of challenging target distributions achieving state-of-the-art sampling performance. In the cases of primordial feature analysis and gravitational wave inference, PMC is approximately 50 and 25 times faster respectively than Nested Sampling (NS). We found that in higher dimensional applications the acceleration is even greater. Finally, PMC is directly parallelisable, manifesting linear scaling up to thousands of CPUs.

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&#39;, 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.

preprint2022arXiv

pocoMC: A Python package for accelerated Bayesian inference in astronomy and cosmology

pocoMC is a Python package for accelerated Bayesian inference in astronomy and cosmology. The code is designed to sample efficiently from posterior distributions with non-trivial geometry, including strong multimodality and non-linearity. To this end, pocoMC relies on the Preconditioned Monte Carlo algorithm which utilises a Normalising Flow in order to decorrelate the parameters of the posterior. It facilitates both tasks of parameter estimation and model comparison, focusing especially on computationally expensive applications. It allows fitting arbitrary models defined as a log-likelihood function and a log-prior probability density function in Python. Compared to popular alternatives (e.g. nested sampling) pocoMC can speed up the sampling procedure by orders of magnitude, cutting down the computational cost substantially. Finally, parallelisation to computing clusters manifests linear scaling.

preprint2021arXiv

The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: measurement of the BAO and growth rate of structure of the emission line galaxy sample from the anisotropic power spectrum between redshift 0.6 and 1.1

We analyse the large-scale clustering in Fourier space of emission line galaxies (ELG) from the Data Release 16 of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey. The ELG sample contains 173,736 galaxies covering 1,170 square degrees in the redshift range $0.6 < z < 1.1$. We perform a BAO measurement from the post-reconstruction power spectrum monopole, and study redshift space distortions (RSD) in the first three even multipoles. Photometric variations yield fluctuations of both the angular and radial survey selection functions. Those are directly inferred from data, imposing integral constraints which we model consistently. The full data set has only a weak preference for a BAO feature ($1.4σ$). At the effective redshift $z_{\rm eff} = 0.845$ we measure $D_{\rm V}(z_{\rm eff})/r_{\rm drag} = 18.33_{-0.62}^{+0.57}$, with $D_{\rm V}$ the volume-averaged distance and $r_{\rm drag}$ the comoving sound horizon at the drag epoch. In combination with the RSD measurement, at $z_{\rm eff} = 0.85$ we find $fσ_8(z_{\rm eff}) = 0.289_{-0.096}^{+0.085}$, with $f$ the growth rate of structure and $σ_8$ the normalisation of the linear power spectrum, $D_{\rm H}(z_{\rm eff})/r_{\rm drag} = 20.0_{-2.2}^{+2.4}$ and $D_{\rm M}(z_{\rm eff})/r_{\rm drag} = 19.17 \pm 0.99$ with $D_{\rm H}$ and $D_{\rm M}$ the Hubble and comoving angular distances, respectively. These results are in agreement with those obtained in configuration space, thus allowing a consensus measurement of $fσ_8(z_{\rm eff}) = 0.315 \pm 0.095$, $D_{\rm H}(z_{\rm eff})/r_{\rm drag} = 19.6_{-2.1}^{+2.2}$ and $D_{\rm M}(z_{\rm eff})/r_{\rm drag} = 19.5 \pm 1.0$. This measurement is consistent with a flat $Λ$CDM model with Planck parameters.

preprint2020arXiv

Assessing non-linear models for galaxy clustering II: model validation and forecasts for Stage IV surveys

Accurate modelling of non-linear scales in galaxy clustering will be crucial for data analysis of Stage IV galaxy surveys. A selection of competing non-linear models must be made based on validation studies. We provide a comprehensive set of forecasts of two different models for the halo redshift space power spectrum, namely the commonly applied TNS model and an effective field theory of large scale structure (EFTofLSS) inspired model. Using simulation data and a least-$χ^2$ analysis, we determine ranges of validity for the models. We then conduct an exploratory Fisher analysis using the full anisotropic power spectrum to investigate parameter degeneracies. We proceed to perform an MCMC analysis utilising the monopole, quadrupole, and hexadecapole spectra, with a restricted range of scales for the latter in order to avoid biasing our growth rate, $f$, constraint. We find that the TNS model with a Lorentzian damping and standard Eulerian perturbative modelling outperforms other variants of the TNS model. Our MCMC analysis finds that the EFTofLSS-based model may provide tighter marginalised constraints on $f$ at $z=0.5$ and $z=1$ than the TNS model, despite having additional nuisance parameters. However this depends on the range of scales used as well as the fiducial values and priors on the EFT nuisance parameters. Finally, we extend previous work to provide a consistent comparison between the Fisher matrix and MCMC forecasts using the multipole expansion formalism, and find good agreement between them.

preprint2020arXiv

Modeling relativistic contributions to the halo power spectrum dipole

We study the power spectrum dipole of an N-body simulation which includes relativistic effects through ray-tracing and covers the low redshift Universe up to $z_{\rm max} = 0.465$ (RayGalGroup simulation). We model relativistic corrections as well as wide-angle, evolution, window and lightcone effects. Our model includes all relativistic corrections up to third-order including third-order bias expansion. We consider all terms which depend linearly on $\mathcal{H}/k$ (weak field approximation). We also study the impact of 1-loop corrections to the matter power spectrum for the gravitational redshift and transverse Doppler effect. We found wide-angle and window function effects to significantly contribute to the dipole signal. When accounting for all contributions, our dipole model can accurately capture the gravitational redshift and Doppler terms up to the smallest scales included in our comparison ($k=0.48\,h{\rm Mpc}^{-1}$), while our model for the transverse Doppler term is less accurate. We find the Doppler term to be the dominant signal for this low redshift sample. We use Fisher matrix forecasts to study the potential for the future Dark Energy Spectroscopic Instrument (DESI) to detect relativistic contributions to the power spectrum dipole. A conservative estimate suggests that the DESI-BGS sample should be able to have a detection of at least $4.4σ$, while more optimistic estimates find detections of up to $10σ$. Detecting these effects in the galaxy distribution allows new tests of gravity on the largest scales, providing an interesting additional science case for galaxy survey experiments.

preprint2020arXiv

Perturbation theory approach to predict the covariance matrices of the galaxy power spectrum and bispectrum in redshift space

In this paper, we predict the covariance matrices of both the power spectrum and the bispectrum, including full non-Gaussian contributions, redshift space distortions, linear bias effects and shot-noise corrections, using perturbation theory (PT). To quantify the redshift-space distortion effect, we focus mainly on the monopole and quadrupole components of both the power and bispectra. We, for the first time, compute the 5- and 6-point spectra to predict the cross-covariance between the power and bispectra, and the auto-covariance of the bispectrum in redshift space. We test the validity of our calculations by comparing them with the covariance matrices measured from the MultiDark-Patchy mock catalogues that are designed to reproduce the galaxy clustering measured from the Baryon Oscillation Spectroscopic Survey Data Release 12. We argue that the simple, leading-order perturbation theory works because the shot-noise corrections for the Patchy mocks are more dominant than other higher-order terms we ignore. In the meantime, we confirm some discrepancies in the comparison, especially of the cross-covariance. We discuss potential sources of such discrepancies. We also show that our PT model reproduces well the cumulative signal-to-noise of the power spectrum and the bispectrum as a function of maximum wavenumber, implying that our PT model captures successfully essential contributions to the covariance matrices.

preprint2020arXiv

Testing Low-Redshift Cosmic Acceleration with Large-Scale Structure

We examine the cosmological implications of measurements of the void-galaxy cross-correlation at redshift $z=0.57$ combined with baryon acoustic oscillation (BAO) data at $0.1<z<2.4$. We find direct evidence of the late-time acceleration due to dark energy at $>10σ$ significance from these data alone, independent of the cosmic microwave background and supernovae. Using a nucleosynthesis prior on $Ω_bh^2$, we measure the Hubble constant to be $H_0=72.3\pm1.9\;{\rm km\,s}^{-1}{\rm Mpc}^{-1}$ from BAO+voids at $z<2$, and $H_0=69.0\pm1.2\;{\rm km\,s}^{-1}{\rm Mpc}^{-1}$ when adding Lyman-$α$ BAO at $z=2.34$, both independent of the CMB. Adding voids to CMB, BAO and supernova data greatly improves measurement of the dark energy equation of state, increasing the figure of merit by >40%, but remaining consistent with flat flat $Λ$ cold dark matter.

preprint2020arXiv

The Sixteenth Data Release of the Sloan Digital Sky Surveys: First Release from the APOGEE-2 Southern Survey and Full Release of eBOSS Spectra

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library &#34;MaStar&#34;). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

preprint2019arXiv

Efficient Cosmological Analysis of the SDSS/BOSS data from the Effective Field Theory of Large-Scale Structure

The precision of the cosmological data allows us to accurately approximate the predictions for cosmological observables by Taylor expanding up to a low order the dependence on the cosmological parameters around a reference cosmology. By applying this observation to the redshift-space one-loop galaxy power spectrum of the Effective Field Theory of Large-Scale Structure, we analyze the BOSS DR12 data by scanning over all the parameters of $Λ$CDM cosmology with massive neutrinos. We impose several sets of priors, the widest of which is just a Big Bang Nucleosynthesis prior on the current fractional energy density of baryons, $Ω_b h^2$, and a bound on the sum of neutrino masses to be less than 0.9 eV. In this case we measure the primordial amplitude of the power spectrum, $A_s$, the abundance of matter, $Ω_m$, the Hubble parameter, $H_0$, and the tilt of the primordial power spectrum, $n_s$, to about $19\%$, $5.7\%$, $2.2\%$ and $7.3\%$ respectively, obtaining $\ln ( 10^{10} A_s) =2.91\pm 0.19$, $Ω_m=0.314\pm 0.018$, $H_0=68.7\pm 1.5$ km/(s Mpc) and $n_s=0.979\pm 0.071$ at $68\%$ confidence level. A public code is released with this preprint.

preprint2019arXiv

The Cosmological Analysis of the SDSS/BOSS data from the Effective Field Theory of Large-Scale Structure

The Effective Field Theory of Large-Scale Structure (EFTofLSS) is a formalism that allows us to predict the clustering of Cosmological Large-Scale Structure in the mildly non-linear regime in an accurate and reliable way. After validating our technique against several sets of numerical simulations, we perform the analysis for the cosmological parameters of the DR12 BOSS data. We assume $Λ$CDM, a fixed value of the baryon/dark-matter ratio, $Ω_b/Ω_c$, and of the tilt of the primordial power spectrum, $n_s$, and no significant input from numerical simulations. By using the one-loop power spectrum multipoles, we measure the primordial amplitude of the power spectrum, $A_s$, the abundance of matter, $Ω_m$, and the Hubble parameter, $H_0$, to about $13\%$, $3.2\%$ and $3.2\%$ respectively, obtaining $\ln(10^{10}As)=2.72\pm 0.13$, $Ω_m=0.309\pm 0.010$, $H_0=68.5\pm 2.2$ km/(s Mpc) at 68\% confidence level. If we then add a CMB prior on the sound horizon, the error bar on $H_0$ is reduced to $1.6\%$. These results are a substantial qualitative and quantitative improvement with respect to former analyses, and suggest that the EFTofLSS is a powerful instrument to extract cosmological information from Large-Scale Structure.

preprint2019arXiv

The Impact of the Fiducial Cosmology Assumption on BAO Cosmological Parameter Inference

Standard analysis pipelines for measurements of Baryon Acoustic Oscillations (BAO) in galaxy surveys make use of a fiducial cosmological model to guide the data compression required to transform from observed redshifts and angles to the measured angular and radial BAO peak positions. In order to remove any dependence on the fiducial cosmology from the results, all models compared to the data should mimic the compression and its dependence on the fiducial model. In practice, approximations are made when testing models: (1) There is assumed to be no residual dependence on the fiducial cosmology after reconstruction, (2) differences in the distance--redshift relationship are assumed to match a linear scaling, and (3) differences in clustering between true and fiducial models are assumed to be removed by the free parameters used to null the non-BAO signal. We test these approximations using the current standard measurement procedure with a set of halo catalogs from the {\sc Aemulus} suite of $N$-body simulations, which span a range of $w\mathrm{CDM}$ cosmological models. We focus on reconstruction of the primordial BAO and locating the BAO. For the range of $w\mathrm{CDM}$ cosmologies covered by the {\sc Aemulus} suite, we find no evidence for systematic errors in the measured BAO shift parameters $α_{\parallel}$ and $α_{\bot}$ to $< 0.1\%$. However, the measured errors $σ_{α_{\parallel}}$ and $σ_{α_{\bot}}$ show a notable absolute increase by up to $+0.001$ and $+0.002$ respectively in the case that the fiducial cosmology does not match the truth. These effects on the inferred BAO scale will be important given the precision of measurements expected from future surveys including DESI, Euclid, and WFIRST.