Researcher profile

Benjamin D. Wandelt

Benjamin D. Wandelt contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
22works
0followers
10topics
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

22 published item(s)

preprint2022arXiv

Bayesian Control Variates for optimal covariance estimation with pairs of simulations and surrogates

Predictions of the mean and covariance matrix of summary statistics are critical for confronting cosmological theories with observations, not least for likelihood approximations and parameter inference. The price to pay for accurate estimates is the extreme cost of running $N$-body and hydrodynamics simulations. Approximate solvers, or surrogates, greatly reduce the computational cost but can introduce significant biases, for example in the non-linear regime of cosmic structure growth. We propose "CARPool Bayes", an approach to solve the inference problem for both the means and covariances using a combination of simulations and surrogates. Our framework allows incorporating prior information for the mean and covariance. We derive closed-form solutions for Maximum A Posteriori covariance estimates that are efficient Bayesian shrinkage estimators, guarantee positive semi-definiteness, and can optionally leverage analytical covariance approximations. We discuss choices of the prior and propose a simple procedure for obtaining optimal prior hyperparameter values with a small set of test simulations. We test our method by estimating the covariances of clustering statistics of GADGET-III $N$-body simulations at redshift $z=0.5$ using surrogates from a 100-1000$\times$ faster particle-mesh code. Taking the sample covariance from 15,000 simulations as the truth, and using an empirical Bayes prior with diagonal blocks, our estimator produces nearly identical Fisher matrix contours for $Λ$CDM parameters using only $15$ simulations of the non-linear dark matter power spectrum. In this case the number of simulations is so small that the sample covariance would be degenerate. We show cases where even with a naïve prior our method still improves the estimate. Our framework is applicable to a wide range of cosmological and astrophysical problems where fast surrogates are available.

preprint2022arXiv

Breaking baryon-cosmology degeneracy with the electron density power spectrum

Uncertain feedback processes in galaxies affect the distribution of matter, currently limiting the power of weak lensing surveys. If we can identify cosmological statistics that are robust against these uncertainties, or constrain these effects by other means, then we can enhance the power of current and upcoming observations from weak lensing surveys such as DES, Euclid, the Rubin Observatory, and the Roman Space Telescope. In this work, we investigate the potential of the electron density auto-power spectrum as a robust probe of cosmology and baryonic feedback. We use a suite of (magneto-)hydrodynamic simulations from the CAMELS project and perform an idealized analysis to forecast statistical uncertainties on a limited set of cosmological and physically-motivated astrophysical parameters. We find that the electron number density auto-correlation, measurable through either kinematic Sunyaev-Zel'dovich observations or through Fast Radio Burst dispersion measures, provides tight constraints on $Ω_{m}$ and the mean baryon fraction in intermediate-mass halos, $\bar{f}_{\mathrm{bar}}$. By obtaining an empirical measure for the associated systematic uncertainties, we find these constraints to be largely robust to differences in baryonic feedback models implemented in hydrodynamic simulations. We further discuss the main caveats associated with our analysis, and point out possible directions for future work.

preprint2022arXiv

CARPool Covariance: Fast, unbiased covariance estimation for large-scale structure observables

The covariance matrix $\boldsymbolΣ$ of non-linear clustering statistics that are measured in current and upcoming surveys is of fundamental interest for comparing cosmological theory and data and a crucial ingredient for the likelihood approximations underlying widely used parameter inference and forecasting methods. The extreme number of simulations needed to estimate $\boldsymbolΣ$ to sufficient accuracy poses a severe challenge. Approximating $\boldsymbolΣ$ using inexpensive but biased surrogates introduces model error with respect to full simulations, especially in the non-linear regime of structure growth. To address this problem we develop a matrix generalization of Convergence Acceleration by Regression and Pooling (CARPool) to combine a small number of simulations with fast surrogates and obtain low-noise estimates of $\boldsymbolΣ$ that are unbiased by construction. Our numerical examples use CARPool to combine GADGET-III $N$-body simulations with fast surrogates computed using COmoving Lagrangian Acceleration (COLA). Even at the challenging redshift $z=0.5$, we find variance reductions of at least $\mathcal{O}(10^1)$ and up to $\mathcal{O}(10^4)$ for the elements of the matter power spectrum covariance matrix on scales $8.9\times 10^{-3}<k_\mathrm{max} <1.0$ $h {\rm Mpc^{-1}}$. We demonstrate comparable performance for the covariance of the matter bispectrum, the matter correlation function and probability density function of the matter density field. We compare eigenvalues, likelihoods, and Fisher matrices computed using the CARPool covariance estimate with the standard sample covariance estimators and generally find considerable improvement except in cases where $Σ$ is severely ill-conditioned.

preprint2022arXiv

Implicit Likelihood Inference of Reionization Parameters from the 21 cm Power Spectrum

The first measurements of the 21 cm brightness temperature power spectrum from the epoch of reionization will very likely be achieved in the near future by radio interferometric array experiments such as the Hydrogen Epoch of Reionization Array (HERA) and the Square Kilometre Array (SKA). Standard MCMC analyses use an explicit likelihood approximation to infer the reionization parameters from the 21 cm power spectrum. In this paper, we present a new Bayesian inference of the reionization parameters where the likelihood is implicitly defined through forward simulations using density estimation likelihood-free inference (DELFI). Realistic effects including thermal noise and foreground avoidance are also applied to the mock observations from the HERA and SKA. We demonstrate that this method recovers accurate posterior distributions for the reionization parameters, and outperforms the standard MCMC analysis in terms of the location and size of credible parameter regions. With the minutes-level processing time once the network is trained, this technique is a promising approach for the scientific interpretation of future 21 cm power spectrum observation data. Our code 21cmDELFI-PS is publicly available at this link.

preprint2022arXiv

Snowmass2021 CMB-HD White Paper

CMB-HD is a proposed millimeter-wave survey over half the sky that would be ultra-deep (0.5 uK-arcmin) and have unprecedented resolution (15 arcseconds at 150 GHz). Such a survey would answer many outstanding questions about the fundamental physics of the Universe. Major advances would be 1.) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10 h Mpc^(-1)), which probes dark matter particle properties. It will also allow 2.) measurements of the thermal and kinetic Sunyaev-Zel&#39;dovich effects on small scales to map the gas density and velocity, another probe of cosmic structure. In addition, CMB-HD would allow us to cross critical thresholds: 3.) ruling out or detecting any new, light (< 0.1 eV) particles that were in thermal equilibrium with known particles in the early Universe, 4.) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe, and 5.) ruling out or detecting inflationary magnetic fields. CMB-HD would also provide world-leading constraints on 6.) axion-like particles, 7.) cosmic birefringence, 8.) the sum of the neutrino masses, and 9.) the dark energy equation of state. The CMB-HD survey would be delivered in 7.5 years of observing 20,000 square degrees of sky, using two new 30-meter-class off-axis crossed 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.

preprint2022arXiv

The Cosmic Graph: Optimal Information Extraction from Large-Scale Structure using Catalogues

We present an implicit likelihood approach to quantifying cosmological information over discrete catalogue data, assembled as graphs. To do so, we explore cosmological parameter constraints using mock dark matter halo catalogues. We employ Information Maximising Neural Networks (IMNNs) to quantify Fisher information extraction as a function of graph representation. We a) demonstrate the high sensitivity of modular graph structure to the underlying cosmology in the noise-free limit, b) show that graph neural network summaries automatically combine mass and clustering information through comparisons to traditional statistics, c) demonstrate that networks can still extract information when catalogues are subject to noisy survey cuts, and d) illustrate how nonlinear IMNN summaries can be used as asymptotically optimal compressed statistics for Bayesian simulation-based inference. We reduce the area of joint $Ω_m, σ_8$ parameter constraints with small ($\sim$100 object) halo catalogues by a factor of 42 over the two-point correlation function, and demonstrate that the networks automatically combine mass and clustering information. This work utilises a new IMNN implementation over graph data in Jax, which can take advantage of either numerical or auto-differentiability. We also show that graph IMNNs successfully compress simulations away from the fiducial model at which the network is fitted, indicating a promising alternative to n-point statistics in catalogue simulation-based analyses.

preprint2021arXiv

Accurate precision Cosmology with redshift unknown gravitational wave sources

Gravitational waves can provide an accurate measurement of the luminosity distance to the source, but cannot provide the source redshift unless the degeneracy between mass and redshift can be broken. This makes it essential to infer the redshift of the source independently to measure the expansion history of the Universe. We show that by exploiting the clustering scale of the gravitational wave sources with galaxies of known redshift, we can infer the expansion history from redshift unknown gravitational wave sources. By using gravitational wave sources of unknown redshift that are detectable from the network of gravitational wave detectors with Advanced LIGO design sensitivity, we will be able to obtain accurate and precise measurements of the local Hubble constant, the expansion history of the universe, and the gravitational wave bias parameter, which captures the distribution of gravitational wave sources with respect to the redshift tracer distribution. While we showcase its application to low redshift gravitational waves, this technique will be applicable also to the high redshift gravitational wave sources detectable from Laser Interferometer Space Antenna (LISA), Cosmic Explorer (CE), and Einstein Telescope (ET). Moreover, this method will also be applicable to samples of supernovae and fast radio bursts with unknown or photometric redshifts.

preprint2021arXiv

Simulation-Based Inference of Reionization Parameters From 3D Tomographic 21 cm Lightcone Images

Tomographic three-dimensional 21 cm images from the epoch of reionization contain a wealth of information about the reionization of the intergalactic medium by astrophysical sources. Conventional power spectrum analysis cannot exploit the full information in the 21 cm data because the 21 cm signal is highly non-Gaussian due to reionization patchiness. We perform a Bayesian inference of the reionization parameters where the likelihood is implicitly defined through forward simulations using density estimation likelihood-free inference (DELFI). We adopt a trained 3D Convolutional Neural Network (CNN) to compress the 3D image data into informative summaries (DELFI-3D CNN). We show that this method recovers accurate posterior distributions for the reionization parameters. Our approach outperforms earlier analysis based on two-dimensional 21 cm images. In contrast, an MCMC analysis of the 3D lightcone-based 21 cm power spectrum alone and using a standard explicit likelihood approximation results in less accurate credible parameter regions than inferred by the DELFI-3D CNN, both in terms of the location and shape of the contours. Our proof-of-concept study implies that the DELFI-3D CNN can effectively exploit more information in the 3D 21 cm images than a 2D CNN or power spectrum analysis. This technique can be readily extended to include realistic effects and is therefore a promising approach for the scientific interpretation of future 21 cm observation data.

preprint2021arXiv

SSSpaNG! Stellar Spectra as Sparse, data-driven, Non-Gaussian processes

Upcoming million-star spectroscopic surveys have the potential to revolutionize our view of the formation and chemical evolution of the Milky Way. Realizing this potential requires automated approaches to optimize estimates of stellar properties, such as chemical element abundances, from the spectra. The volume and quality of the observations strongly motivate that these approaches should be data-driven. With this in mind, we introduce SSSpaNG: a data-driven non-Gaussian Process model of stellar spectra. We demonstrate the capabilities of SSSpaNG using a sample of APOGEE red clump stars, whose model parameters we infer via Gibbs sampling. Pooling information between stars to infer their covariance, we permit clear identification of the correlations between spectral pixels. Harnessing these correlations, we infer the true spectrum of each star, inpainting missing regions and denoising by a factor of at least 2 for stars with signal-to-noise of ~20. As we marginalize over the covariance matrix of the spectra, the effective prior on these true spectra is non-Gaussian and sparsifying, favouring typically small but occasionally large excursions from the mean. The high-fidelity inferred spectra produced will enable improved elemental abundance measurements for individual stars. Our model also allows us to quantify the information gained by observing portions of a star&#39;s spectrum, and thereby define the most mutually informative spectral regions. Using 25 windows centred on elemental absorption lines, we demonstrate that the iron-peak and alpha-process elements are particularly mutually informative for these spectra, and that the majority of information about a target window is contained in the 10-or-so most informative windows. Such mutual-information estimates have the potential to inform models of nucleosynthetic yields and the design of future observations.

preprint2021arXiv

Testing the general theory of relativity using gravitational wave propagation from dark standard sirens

Alternative theories of gravity predict modifications in the propagation of gravitational waves (GW) through space-time. One of the smoking-gun predictions of such theories is the change in the GW luminosity distance to GW sources as a function of redshift relative to the electromagnetic (EM) luminosity distance expected from EM probes. We propose a multi-messenger test of the theory of general relativity from the propagation of gravitational waves by combining EM and GW observations to resolve these issues from GW sources without EM counterparts (which are also referred to as dark standard sirens). By using the relation between the geometric distances accessible from baryon acoustic oscillation measurements, and luminosity distance measurements from the GW sources, we can measure any deviation from the general theory of relativity via the GW sources of unknown redshift that will be detectable by networks of GW detectors such as LIGO, Virgo, and KAGRA. Using this technique, the fiducial value of the frictional term can be measured to a precision $Ξ_0=0.98^{+0.04}_{-0.23}$ after marginalizing over redshift dependence, cosmological parameters, and GW bias parameters with $\sim 3500$ dark standard sirens of masses $30\,\rm M_\odot$ each distributed up to redshift $z=0.5$. For a fixed redshift dependence, a value of $Ξ_0=0.99^{+0.02}_{-0.02}$ can be measured with a similar number of dark sirens. Application of our methodology to the far more numerous dark standard sirens detectable with next generation GW detectors, such as LISA, Einstein Telescope and Cosmic Explorer, will allow achievement of higher accuracy than possible from use of bright standard sirens.

preprint2021arXiv

The Age of the Universe with Globular Clusters: Reducing Systematic Uncertainties

The dominant systematic uncertainty in the age determination of galactic globular clusters is the depth of the convection envelope of the stars. This parameter is partially degenerate with metallicity which is in turn degenerate with age. However, if the metal content, distance and extinction are known, the position and morphology of the red giant branch in a color-magnitude diagram are mostly sensitive to the value of the depth of the convective envelope. Therefore, using external, precise metallicity determinations this degeneracy and thus the systematic error in age, can be reduced. Alternatively, the morphology of the red giant branch of globular clusters color magnitude diagram can also be used to achieve the same. We demonstrate that globular cluster red giant branches are well fitted by values of the depth of the convection envelope consistent with those obtained for the Sun and this finding is robust to the adopted treatment of the stellar physics. With these findings, the uncertainty in the depth of the convection envelope is no longer the dominant contribution to the systematic error in the age determination of the oldest globular clusters, reducing it from $0.5$ to $0.23$ or $0.33$ Gyr, depending on the methodology adopted: i.e., whether resorting to external data (spectroscopic metallicity determinations) or relying solely on the morphology of the clusters&#39;s color-magnitude diagrams. This results in an age of the Universe $t_{\rm U}=13.5^{+0.16}_{-0.14} {\rm (stat.)} \pm 0.23(0.33) ({\rm sys.})$ at 68\% confidence level, accounting for the formation time of globular clusters and its uncertainty. An uncertainty of 0.27(0.36) Gyr if added in quadrature. This agrees well with $13.8 \pm 0.02$ Gyr, the cosmological model-dependent value inferred by the Planck mission assuming the $Λ$CDM model.

preprint2021arXiv

The GIGANTES dataset: precision cosmology from voids in the machine learning era

We present GIGANTES, the most extensive and realistic void catalog suite ever released -- containing over 1 billion cosmic voids covering a volume larger than the observable Universe, more than 20 TB of data, and created by running the void finder VIDE on QUIJOTE&#39;s halo simulations. The expansive and detailed GIGANTES suite, spanning thousands of cosmological models, opens up the study of voids, answering compelling questions: Do voids carry unique cosmological information? How is this information correlated with galaxy information? Leveraging the large number of voids in the GIGANTES suite, our Fisher constraints demonstrate voids contain additional information, critically tightening constraints on cosmological parameters. We use traditional void summary statistics (void size function, void density profile) and the void auto-correlation function, which independently yields an error of $0.13\,\mathrm{eV}$ on $\sum\,m_ν$ for a 1 $h^{-3}\mathrm{Gpc}^3$ simulation, without CMB priors. Combining halos and voids we forecast an error of $0.09\,\mathrm{eV}$ from the same volume. Extrapolating to next generation multi-Gpc$^3$ surveys such as DESI, Euclid, SPHEREx, and the Roman Space Telescope, we expect voids should yield an independent determination of neutrino mass. Crucially, GIGANTES is the first void catalog suite expressly built for intensive machine learning exploration. We illustrate this by training a neural network to perform likelihood-free inference on the void size function. Cosmology problems provide an impetus to develop novel deep learning techniques, leveraging the symmetries embedded throughout the universe from physical laws, interpreting models, and accurately predicting errors. With GIGANTES, machine learning gains an impressive dataset, offering unique problems that will stimulate new techniques.

preprint2020arXiv

A new probe of Axion-Like Particles: CMB polarization distortions due to cluster magnetic fields

We propose using the upcoming Cosmic Microwave Background (CMB) ground based experiments to detect the signal of ALPs (Axion like particles) interacting with magnetic fields in galaxy clusters. The conversion between CMB photons and ALPs in the presence of the cluster magnetic field can cause a polarized spectral distortion in the CMB around a galaxy cluster. The strength of the signal depends upon the redshift of the galaxy cluster and will exhibit a distinctive spatial profile around it depending upon the structure of electron density and magnetic field. This distortion produces a different shape from the other known spectral distortions like $y$-type and $μ$-type and hence are separable from the multi-frequency CMB observation. The spectrum is close to kinematic Sunyaev-Zeldovich (kSZ) signal but can be separated from it using the polarization information. For the future ground-based CMB experiments such as Simons Observatory and CMB-S4, we estimate the measurability of this signal in the presence of foreground contamination, instrument noise and CMB anisotropies. This new avenue can probe the photon-ALP coupling over the ALP mass range from $10^{-13}$ eV to $10^{-12}$ eV with two orders of magnitude better accuracy from CMB-S4 than the current existing bounds.

preprint2020arXiv

Bayesian delensing delight: sampling-based inference of the primordial CMB and gravitational lensing

The search for primordial gravitational waves in the Cosmic Microwave Background (CMB) will soon be limited by our ability to remove the lensing contamination to $B$-mode polarization. The often-used quadratic estimator for lensing is known to be suboptimal for surveys that are currently operating and will continue to become less and less efficient as instrumental noise decreases. While foregrounds can in principle be mitigated by observing in more frequency bands, progress in delensing hinges entirely on algorithmic advances. We demonstrate here a new inference method that solves this problem by sampling the exact Bayesian posterior of any desired cosmological parameters, of the gravitational lensing potential, and of the delensed CMB maps, given lensed temperature and polarization data. We validate the method using simulated CMB data with non-white noise and masking on up to 650\,deg$^2$ patches of sky. A unique strength of this approach is the ability to jointly estimate cosmological parameters which control both the primordial CMB and the lensing potential, which we demonstrate here for the first time by sampling both the tensor-to-scalar ratio, $r$, and the amplitude of the lensing potential, $A_ϕ$. The method allows us to perform the most precise check to-date of several important approximations underlying CMB-S4 $r$ forecasting, and we confirm these yield the correct expected uncertainty on $r$ to better than 10%.

preprint2020arXiv

How to measure CMB spectral distortions with an imaging telescope

We propose the implementation of an imaging telescope in combination with an inter-frequency calibrator to measure the spectral shape of the microwave sky by exploiting the differences in the sky intensity between multiple pairs of frequency channels. By jointly sampling the cosmological and foreground parameters in a Bayesian framework for $600$ instrument configurations, we determine the minimum calibration accuracy required in order to obtain measurements of spectral distortions and simultaneously measure spectral and spatial fluctuations of the CMB. We demonstrate the feasibility of this technique for a CMB mission like PICO (Probe of Inflation and Cosmic Origins), and show that a $10$-$σ$ measurement of the $y$-distortion along with a two orders of magnitude improvement on the FIRAS (Far-Infrared Absolute Spectrophotometer) $μ$-distortion limit is feasible from this technique. We argue that longer term applications may be envisaged at even higher sensitivity, capable of attaining the $μ$-distortion that provide a robust prediction of the $Λ$CDM model and even primordial recombination lines of hydrogen and helium, in the context of the ESA (European Space Agency) Voyage 2035-2050 program.

preprint2020arXiv

Inferring the Age of the Universe with Globular Clusters

We present an estimate of the absolute age of 68 galactic globular clusters obtained by exploiting the distribution of stars in the full color-magnitude diagram. In particular, we jointly estimate the absolute age, distance, reddening, metallicity ([Fe/H]) and [$α$/Fe] of each cluster, imposing priors motivated by independent observations; we also estimate possible systematics from stellar modeling. Our derived distances for the globular cluster sample are in agreement with those obtained from GAIA using main-sequence dwarf stars (where available), and the inferred ages are in good agreement with those previously published. The novelty of our approach is that, with the adopted priors, we are able to estimate robustly these parameters from the globular cluster color-magnitude diagram. We find that the average age of the oldest globular clusters is $t_{\rm GC}=13.32 \pm 0.1 {\rm (stat.)} \pm 0.5 {\rm (sys.)}$, at 68\% confidence level, including systematic uncertainties from stellar modeling. These measurements can be used to infer the age of the Universe, largely independently of the cosmological parameters: we find an age of the Universe $t_{\rm U}=13.5^{+0.16}_{-0.14} {\rm (stat.)} \pm 0.5 ({\rm sys.})$ at 68\% confidence level, accounting for the formation time of globular clusters and its uncertainty. This value is compatible with $13.8 \pm 0.02$ Gyr, the cosmological model-dependent value inferred by the Planck mission assuming the $Λ$CDM model.

preprint2020arXiv

Multi-messenger tests of gravity with weakly lensed gravitational waves

General relativity (GR) predicts concordant trajectories for photons and gravitational waves (GW). We propose a new multi-messenger avenue (GW-CMB-CMB) to prove this aspect of fundamental physics by cross-correlating the GW signal of astrophysical origin with the lensing field derived from the cosmic microwave background (CMB). This new window will allow robust measurement of the prediction from GR with high signal-to-noise and will be able to unveil the true nature of gravity using the GW sources detected by missions such as the Laser Interferometer Space Antenna (LISA), Einstein Telescope and Cosmic Explorer.

preprint2020arXiv

Probing the theory of gravity with gravitational lensing of gravitational waves and galaxy surveys

The cross-correlation of gravitational wave strain with upcoming galaxy surveys probe theories of gravity in a new way. This method enables testing the theory of gravity by combining the effects from both gravitational lensing of gravitational waves and the propagation of gravitational waves in spacetime. We find that within 10 years, the combination of the Advanced-LIGO and VIRGO detector networks with planned galaxy surveys should detect weak gravitational lensing of gravitational waves in the low redshift Universe ($z<0.5$). With the next generation gravitational wave experiments such as Voyager, LISA, Cosmic-Explorer and Einstein Telescope, we can extend this test of the theory of gravity to larger redshifts by exploiting the synergies between electromagnetic wave and gravitational wave probes.

preprint2020arXiv

Super-resolution emulator of cosmological simulations using deep physical models

We present an extension of our recently developed Wasserstein optimized model to emulate accurate high-resolution features from computationally cheaper low-resolution cosmological simulations. Our deep physical modelling technique relies on restricted neural networks to perform a mapping of the distribution of the low-resolution cosmic density field to the space of the high-resolution small-scale structures. We constrain our network using a single triplet of high-resolution initial conditions and the corresponding low- and high-resolution evolved dark matter simulations from the Quijote suite of simulations. We exploit the information content of the high-resolution initial conditions as a well constructed prior distribution from which the network emulates the small-scale structures. Once fitted, our physical model yields emulated high-resolution simulations at low computational cost, while also providing some insights about how the large-scale modes affect the small-scale structure in real space.

preprint2020arXiv

Velocity correction for Hubble constant measurements from standard sirens

Gravitational wave (GW) sources are an excellent probe of the luminosity distance and offer a novel measure of the Hubble constant, $H_0$. This estimation of $H_0$ from standard sirens requires an accurate estimation of the cosmological redshift of the host galaxy of the GW source, after correcting for its peculiar velocity. Absence of an accurate peculiar velocity correction affects both the precision and accuracy of the measurement of $H_0$, particularly for nearby sources. We propose a framework to incorporate such a peculiar velocity correction for GW sources. A first implementation of our method to the event GW170817 combined with the Very Large Baseline Interferometry (VLBI) observation leads to a revised value of $H_0= 68.3^{+ 4.6}_{-4.5}$ km/s/Mpc. While this revision is minor, it demonstrates that our method makes it possible for obtaining an unbiased and accurate measurements of $H_0$ at the precision required for the standard siren cosmology.

preprint2019arXiv

Neural physical engines for inferring the halo mass distribution function

An ambitious goal in cosmology is to forward-model the observed distribution of galaxies in the nearby Universe today from the initial conditions of large-scale structures. For practical reasons, the spatial resolution at which this can be done is necessarily limited. Consequently, one needs a mapping between the density of dark matter averaged over ~Mpc scales, and the distribution of dark matter halos (used as a proxy for galaxies) in the same region. Here we demonstrate a method for determining the halo mass distribution function by learning the tracer bias between density fields and halo catalogues using a neural bias model. The method is based on the Bayesian analysis of simple, physically motivated, neural network-like architectures, which we denote as neural physical engines, and neural density estimation. As a result, we are able to sample the initial phases of the dark matter density field whilst inferring the parameters describing the halo mass distribution function, providing a fully Bayesian interpretation of both the initial dark matter density distribution and the neural bias model. We successfully run an upgraded BORG inference using our new likelihood and neural bias model with halo catalogues derived from full N-body simulations. We notice orders of magnitude improvement in modelling compared to classical biasing techniques.