Astro arXiv | astro-ph.CO

Sparse Bayesian mass-mapping using trans-dimensional MCMC by Augustin Marignier et al. on Monday 28 November

Uncertainty quantification is a crucial step of cosmological mass-mapping

that is often ignored. Suggested methods are typically only approximate or make

strong assumptions of Gaussianity of the shear field. Probabilistic sampling

methods, such as Markov chain Monte Carlo (MCMC), draw samples form a

probability distribution, allowing for full and flexible uncertainty

quantification, however these methods are notoriously slow and struggle in the

high-dimensional parameter spaces of imaging problems. In this work we use, for

the first time, a trans-dimensional MCMC sampler for mass-mapping, promoting

sparsity in a wavelet basis. This sampler gradually grows the parameter space

as required by the data, exploiting the extremely sparse nature of mass maps in

wavelet space. The wavelet coefficients are arranged in a tree-like structure,

which adds finer scale detail as the parameter space grows. We demonstrate the

trans-dimensional sampler on galaxy cluster-scale images where the planar

modelling approximation is valid. In high-resolution experiments, this method

produces naturally parsimonious solutions, requiring less than 1% of the

potential maximum number of wavelet coefficients and still producing a good fit

to the observed data. In the presence of noisy data, trans-dimensional MCMC

produces a better reconstruction of mass-maps than the standard smoothed

Kaiser-Squires method, with the addition that uncertainties are fully

quantified. This opens up the possibility for new mass maps and inferences

about the nature of dark matter using the new high-resolution data from

upcoming weak lensing surveys such as Euclid.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.13963v1

Near-inflection point inflation and production of dark matter during reheating by Anish Ghoshal et al. on Monday 28 November

We study slow roll single field inflationary scenario and the production of

non-thermal fermionic dark matter, together with standard model Higgs, during

reheating. For the inflationary scenario, we have considered two models of

polynomial potential - one is symmetric about the origin and another one is

not. We fix the coefficients of the potential from the current Cosmic Microwave

Background (CMB) data from Planck/Bicep. Next, we explore the allowed parameter

space on the coupling $(y_\chi)$ with inflaton and mass $(m_\chi)$ of dark

matter (DM) particles $(\chi)$ produced during reheating and satisfying CMB and

several other cosmological constraints.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15061v1

Machine learning cosmic inflation by Ahana Kamerkar et al. on Monday 28 November

We present a machine-learning approach, based on the genetic algorithms (GA),

that can be used to reconstruct the inflationary potential directly from

cosmological data. We create a pipeline consisting of the GA, a primordial code

and a Boltzmann code used to calculate the theoretical predictions, and Cosmic

Microwave Background (CMB) data. As a proof of concept, we apply our

methodology to the Planck CMB data and explore the functional space of

single-field inflationary potentials in a non-parametric, yet analytical way.

We show that the algorithm easily improves upon the vanilla model of quadratic

inflation and proposes slow-roll potentials better suited to the data, while we

confirm the robustness of the Starobinsky inflation model (and other

small-field models). Moreover, using unbinned CMB data, we perform a first

concrete application of the GA by searching for oscillatory features in the

potential in an agnostic way, and find very significant improvements upon the

best featureless potentials, $\Delta \chi^2 < -20$. These encouraging

preliminary results motivate the search for resonant features in the primordial

power spectrum with a multimodal distribution of frequencies. We stress that

our pipeline is modular and can easily be extended to other CMB data sets and

inflationary scenarios, like multifield inflation or theories with higher-order

derivatives.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.14142v1

Cross-correlation of DES Y3 lensing and ACT it Planck thermal Sunyaev Zel'dovich Effect II: Modeling and constraints on halo pressure profiles by S. Pandey et al. on Sunday 27 November

Hot, ionized gas leaves an imprint on the cosmic microwave background via the

thermal Sunyaev Zel'dovich (tSZ) effect. The cross-correlation of gravitational

lensing (which traces the projected mass) with the tSZ effect (which traces the

projected gas pressure) is a powerful probe of the thermal state of ionized

baryons throughout the Universe, and is sensitive to effects such as baryonic

feedback. In a companion paper (Gatti et al. 2021), we present tomographic

measurements and validation tests of the cross-correlation between galaxy shear

measurements from the first three years of observations of the Dark Energy

Survey, and tSZ measurements from a combination of Atacama Cosmology Telescope

and ${\it Planck}$ observations. In this work, we use the same measurements to

constrain models for the pressure profiles of halos across a wide range of halo

mass and redshift. We find evidence for reduced pressure in low mass halos,

consistent with predictions for the effects of feedback from active galactic

nuclei. We infer the hydrostatic mass bias ($B \equiv M_{500c}/M_{\rm SZ}$)

from our measurements, finding $B = 1.8\pm0.1$ when adopting the ${\it

Planck}$-preferred cosmological parameters. We additionally find that our

measurements are consistent with a non-zero redshift evolution of $B$, with the

correct sign and sufficient magnitude to explain the mass bias necessary to

reconcile cluster count measurements with the ${\it Planck}$-preferred

cosmology. Our analysis introduces a model for the impact of intrinsic

alignments (IA) of galaxy shapes on the shear-tSZ correlation. We show that IA

can have a significant impact on these correlations at current noise levels.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2108.01601v2

Marginal Bayesian Statistics Using Masked Autoregressive Flows and Kernel Density Estimators with Examples in Cosmology by Harry Bevins et al. on Sunday 27 November

Cosmological experiments often employ Bayesian workflows to derive

constraints on cosmological and astrophysical parameters from their data. It

has been shown that these constraints can be combined across different probes

such as Planck and the Dark Energy Survey and that this can be a valuable

exercise to improve our understanding of the universe and quantify tension

between multiple experiments. However, these experiments are typically plagued

by differing systematics, instrumental effects and contaminating signals, which

we collectively refer to as `nuisance' components, that have to be modelled

alongside target signals of interest. This leads to high dimensional parameter

spaces, especially when combining data sets, with > 20 dimensions of which only

around 5 correspond to key physical quantities. We present a means by which to

combine constraints from different data sets in a computationally efficient

manner by generating rapid, reusable and reliable marginal probability density

estimators, giving us access to nuisance-free likelihoods. This is possible

through the unique combination of nested sampling, which gives us access to

Bayesian evidences, and the marginal Bayesian statistics code MARGARINE. Our

method is lossless in the signal parameters, resulting in the same posterior

distributions as would be found from a full nested sampling run over all

nuisance parameters, and typically quicker than evaluating full likelihoods. We

demonstrate our approach by applying it to the combination of posteriors from

the Dark Energy Survey and Planck.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2207.11457v3

The Merger Rate of Primordial Black Hole-Neutron Star Binaries in Ellipsoidal-Collapse Dark Matter Halo Models by Saeed Fakhry et al. on Sunday 27 November

Primordial black holes (PBHs), as a potential macroscopic candidate for dark

matter, can encounter other compact objects in dark matter halos because of

their random distribution. Besides, the detection of gravitational waves (GWs)

related to the stellar-mass black hole-neutron star (BH-NS) mergers raises the

possibility that the BHs involved in such events may have a primordial origin.

In this work, we calculate the merger rate of PBH-NS binaries within the

framework of ellipsoidal-collapse dark matter halo models and compare it with

the corresponding results derived from spherical-collapse dark matter halo

models. Our results exhibit that ellipsoidal-collapse dark matter halo models

can potentially amplify the merger rate of PBH-NS binaries in such a way that

it is very close to the range estimated by the LIGO-Virgo observations. While

spherical-collapse dark matter halo models cannot justify PBH-NS merger events

as consistent results with the latest GW data reported by the LIGO-Virgo

collaboration. In addition, we calculate the merger rate of PBH-NS binaries as

a function of PBH mass and fraction within the context of ellipsoidal-collapse

dark matter halo models. The results indicate that PBH-NS merger events with

the mass of $(M_{PBH}\le 5 M_{\odot}, M_{NS}\simeq 1.4 M_{\odot})$ will be

consistent with the LIGO-Virgo observations if $f_{PBH}\simeq 1$. We also show

that to have at least on $(M_{PBH}\simeq 5 M_{\odot}, M_{NS}\simeq 1.4

M_{\odot})$ event in the comoving volume $1 Gpc^{3}$ annually,

ellipsoidal-collapse dark matter halo models constrain the abundance of PBHs as

$f_{PBH} \geq 0.1$.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2209.08909v2

Removal of point source leakage from time-order data filtering by Zhaoxuan Zhang et al. on Sunday 27 November

Time-ordered data (TOD) from ground-based CMB experiments are generally

filtered before map-making to remove or reduce the contamination from the

ground and the atmospheric emissions. However, when the observation region

contains strong point sources, the filtering process will result in

considerable leakage around the point sources in a measured CMB map, and leave

spurious polarization signals. Therefore, such signals need to be assessed and

removed before CMB science exploitation. In this work, we present a new method

that we call "template fitting" and can effectively remove these leakage

signals in pixel domain, not only satisfying the requirement for measuring

primordial gravitational waves from CMB-$B$ modes, but also avoiding

time-consuming operations on TOD.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.13870v1

Logarithmic Duality of the Curvature Perturbation by Shi Pi et al. on Sunday 27 November

We study the comoving curvature perturbation $\mathcal{R}$ in general

single-field inflation models whose potential can be approximated by a

piecewise quadratic potential $V(\varphi)$ by using the $\delta N$ formalism.

We find a general formula for $\mathcal{R}(\delta\varphi)$, which consists of a

sum of logarithmic functions of the field perturbation $\delta\varphi$ at the

point of interest, as well as of its field velocity perturbations $\delta\pi_*$

at the boundaries of each quadratic piece, which are functions of

$\delta\varphi$ through the equations of motion. In some simple cases,

$\mathcal{R}(\delta\varphi)$ reduces to a single logarithm, which yields either

the renowned ``exponential tail'' of the probability distribution function of

$\mathcal{R}$ or the Gumbel distribution.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.13932v1

Threshold of Primordial Black Hole Formation against Velocity Dispersion in Matter-Dominated Era by Tomohiro Harada et al. on Sunday 27 November

We study the effects of velocity dispersion on the formation of primordial

black holes~(PBHs) in a matter-dominated era. The velocity dispersion is

generated through the nonlinear growth of perturbations and has the potential

to impede the gravitational collapse and thereby the formation of PBHs. To make

discussions clear, we consider two distinct length scales. The larger one is

where gravitational collapse occurs which could lead to PBH formation, and the

smaller one is where the velocity dispersion develops due to nonlinear

interactions. We estimate the effect of the velocity dispersion on the PBH

formation by comparing the free-fall timescale and the timescale for a particle

to cross the collapsing region. As a demonstration, we consider a log-normal

power spectrum for the initial density perturbation with the peak value

$\sigma_0^2$ at a scale that corresponds to the larger scale. We find that the

threshold value of the density perturbation $\tilde \delta_{\rm th}$ at the

horizon entry for the PBH formation scales as $\tilde \delta_{\rm th}\propto

\sigma_0^{2/5}$ for $\sigma_0\ll1$.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.13950v1

Sparse Bayesian mass-mapping using trans-dimensional MCMC by Augustin Marignier et al. on Sunday 27 November

Uncertainty quantification is a crucial step of cosmological mass-mapping

that is often ignored. Suggested methods are typically only approximate or make

strong assumptions of Gaussianity of the shear field. Probabilistic sampling

methods, such as Markov chain Monte Carlo (MCMC), draw samples form a

probability distribution, allowing for full and flexible uncertainty

quantification, however these methods are notoriously slow and struggle in the

high-dimensional parameter spaces of imaging problems. In this work we use, for

the first time, a trans-dimensional MCMC sampler for mass-mapping, promoting

sparsity in a wavelet basis. This sampler gradually grows the parameter space

as required by the data, exploiting the extremely sparse nature of mass maps in

wavelet space. The wavelet coefficients are arranged in a tree-like structure,

which adds finer scale detail as the parameter space grows. We demonstrate the

trans-dimensional sampler on galaxy cluster-scale images where the planar

modelling approximation is valid. In high-resolution experiments, this method

produces naturally parsimonious solutions, requiring less than 1% of the

potential maximum number of wavelet coefficients and still producing a good fit

to the observed data. In the presence of noisy data, trans-dimensional MCMC

produces a better reconstruction of mass-maps than the standard smoothed

Kaiser-Squires method, with the addition that uncertainties are fully

quantified. This opens up the possibility for new mass maps and inferences

about the nature of dark matter using the new high-resolution data from

upcoming weak lensing surveys such as Euclid.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.13963v1