Astro arXiv | astro-ph.CO

Consistent and simultaneous modelling of galaxy clustering and galaxy-galaxy lensing with Subhalo Abundance Matching by Sergio Contreras et al. on Wednesday 23 November

The spatial distribution of galaxies and their gravitational lensing signal

offer complementary tests of galaxy formation physics and cosmology. However,

their synergy can only be fully exploited if both probes are modelled

accurately and consistently. In this paper, we demonstrate that this can be

achieved using an extension of Subhalo Abundance Matching, dubbed SHAMe.

Specifically, we use mock catalogues built from the TNG300 hydrodynamical

simulation to show that SHAMe can simultaneously model the multipoles of the

redshift-space galaxy correlation function and galaxy-galaxy lensing, without

noticeable bias within the statistical sampling uncertainties of a SDSS volume

and on scales r = [0.6-30] Mpc/h. Modelling the baryonic processes in

galaxy-galaxy lensing with a baryonification scheme allows SHAMe's range of

validity to be extended to r = [0.1-30] Mpc/h. Remarkably, our model achieves

this level of precision with just five free parameters beyond those describing

the baryonification model. At fixed cosmology, we find that galaxy-galaxy

lensing provides a general consistency test but little additional information

on galaxy modelling parameters beyond that encoded in the redshift-space

multipoles. It does, however, improve constraints if only the projected

correlation function is available, as in surveys with only photometric

redshifts. We expect SHAMe to have a higher fidelity across a wider range of

scales than more traditional methods such as Halo Occupation Distribution

modelling. Thus it should provide a significantly more powerful and more robust

tool for analysing next-generation large-scale surveys.

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

Can denoising diffusion probabilistic models generate realistic astrophysical fields? by Nayantara Mudur et al. on Wednesday 23 November

Score-based generative models have emerged as alternatives to generative

adversarial networks (GANs) and normalizing flows for tasks involving learning

and sampling from complex image distributions. In this work we investigate the

ability of these models to generate fields in two astrophysical contexts: dark

matter mass density fields from cosmological simulations and images of

interstellar dust. We examine the fidelity of the sampled cosmological fields

relative to the true fields using three different metrics, and identify

potential issues to address. We demonstrate a proof-of-concept application of

the model trained on dust in denoising dust images. To our knowledge, this is

the first application of this class of models to the interstellar medium.

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

Euclid: Modelling massive neutrinos in cosmology -- a code comparison by J. Adamek et al. on Wednesday 23 November

The measurement of the absolute neutrino mass scale from cosmological

large-scale clustering data is one of the key science goals of the Euclid

mission. Such a measurement relies on precise modelling of the impact of

neutrinos on structure formation, which can be studied with $N$-body

simulations. Here we present the results from a major code comparison effort to

establish the maturity and reliability of numerical methods for treating

massive neutrinos. The comparison includes eleven full $N$-body implementations

(not all of them independent), two $N$-body schemes with approximate time

integration, and four additional codes that directly predict or emulate the

matter power spectrum. Using a common set of initial data we quantify the

relative agreement on the nonlinear power spectrum of cold dark matter and

baryons and, for the $N$-body codes, also the relative agreement on the

bispectrum, halo mass function, and halo bias. We find that the different

numerical implementations produce fully consistent results. We can therefore be

confident that we can model the impact of massive neutrinos at the sub-percent

level in the most common summary statistics. We also provide a code validation

pipeline for future reference.

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

Efficient exploration of reionization parameters for the upcoming 21 cm observations using a photon conserving semi-numerical model SCRIPT by Barun Maity et al. on Wednesday 23 November

One of the most promising probes to constrain the reionization history of the

universe is the power spectrum of neutral hydrogen 21 cm emission fluctuations.

The corresponding analyses require computationally efficient modelling of

reionization, usually achieved through semi-numerical simulations. We

investigate the capability of one such semi-numerical code, SCRIPT, to

constrain the reionization parameters. Our study involves creating a mock data

set corresponding to the upcoming SKA-Low, followed by a Bayesian inference

method to constrain the model parameters. In particular, we explore in detail

whether the inferred parameters are unbiased with respect to the inputs used

for the mock, and also if the inferences are insensitive to the resolution of

the simulation. We find that the model is successful on both fronts. We also

develop a simple template model of reionization which can mimic the complex

physical processes like inhomogeneous recombinations and radiative feedback and

show that it can recover the global reionization history reliably with moderate

computational cost. Our results are relevant for constraining reionization

using high-quality data expected in future telescopes.

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

Constraints on modified gravity from the BOSS galaxy survey by Lorenzo Piga et al. on Wednesday 23 November

We develop a pipeline to set new constraints on scale-independent modified

gravity, from the galaxy power spectrum in redshift space of BOSS DR12. The

latter is modelled using the effective field theory of large-scale structure up

to 1-loop order in perturbation theory. We test our pipeline on synthetic and

simulated data, to assess systematic biases on the inferred cosmological

parameters due to marginalization and theoretical errors, and we apply it to

the normal branch of the DGP model with a $\Lambda$CDM background. We observe

biased posteriors due to the strong degeneracy between the nDGP parameter

$\Omega_{\rm rc}$ and the primordial amplitude of fluctuations $A_s$. Fixing

the latter to the Planck central value, we obtain $\Omega_{\rm rc}\lesssim 0.2$

at 95$\%$ C.L. We also discuss a procedure to alleviate the prior dependence of

this bound.

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

Quasar Factor Analysis -- An Unsupervised and Probabilistic Quasar Continuum Prediction Algorithm with Latent Factor Analysis by Zechang Sun et al. on Wednesday 23 November

Since their first discovery, quasars have been essential probes of the

distant Universe. However, due to our limited knowledge of its nature,

predicting the intrinsic quasar continua has bottlenecked their usage. Existing

methods of quasar continuum recovery often rely on a limited number of

high-quality quasar spectra, which might not capture the full diversity of the

quasar population. In this study, we propose an unsupervised probabilistic

model, \textit{Quasar Factor Analysis} (QFA), which combines factor analysis

(FA) with physical priors of the intergalactic medium (IGM) to overcome these

limitations. QFA captures the posterior distribution of quasar continua through

generatively modeling quasar spectra. We demonstrate that QFA can achieve the

state-of-the-art performance, $\sim 2\%$ relative error, for continuum

prediction in the Ly$\alpha$ forest region compared to previous methods. We

further fit 90,678 $2<\mathrm{z}<3.5$, SNR$>2$ quasar spectra from Sloan

Digital Sky Survey Data Release 16 and found that for $\sim 30\%$ quasar

spectra where the continua were ill-determined with previous methods, QFA

yields visually more plausible continua. QFA also attains $\lesssim 1\%$ error

in the 1D Ly$\alpha$ power spectrum measurements at $\mathrm{z}\sim 3$ and

$\sim 4\%$ in $\mathrm{z}\sim 2.4$. In addition, QFA determines latent factors

representing more physically motivated than PCA. We investigate the evolution

of the latent factors and report no significant redshift or luminosity

dependency except for the Baldwin effect. The generative nature of QFA also

enables outlier detection robustly; we showed that QFA is effective in

selecting outlying quasar spectra, including damped Ly$\alpha$ systems and

potential Type II quasar spectra.

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

Can denoising diffusion probabilistic models generate realistic astrophysical fields? by Nayantara Mudur et al. on Wednesday 23 November

Score-based generative models have emerged as alternatives to generative

adversarial networks (GANs) and normalizing flows for tasks involving learning

and sampling from complex image distributions. In this work we investigate the

ability of these models to generate fields in two astrophysical contexts: dark

matter mass density fields from cosmological simulations and images of

interstellar dust. We examine the fidelity of the sampled cosmological fields

relative to the true fields using three different metrics, and identify

potential issues to address. We demonstrate a proof-of-concept application of

the model trained on dust in denoising dust images. To our knowledge, this is

the first application of this class of models to the interstellar medium.

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

Nearly Model-Independent Constraints on the Neutral Hydrogen Fraction in the Intergalactic Medium at z sim 5-7 Using Dark Pixel Fractions in Ly α and Ly β Forests by Xiangyu Jin et al. on Wednesday 23 November

Cosmic reionization was the last major phase transition of hydrogen from

neutral to highly ionized in the intergalactic medium (IGM). Current

observations show that the IGM is significantly neutral at $z>7$, and largely

ionized by $z\sim5.5$. However, most methods to measure the IGM neutral

fraction are highly model-dependent, and are limited to when the

volume-averaged neutral fraction of the IGM is either relatively low

($\bar{x}_{\rm HI} \lesssim 10^{-3}$) or close to unity ($\bar{x}_{\rm HI}\sim

1$). In particular, the neutral fraction evolution of the IGM at the critical

redshift range of $z=6-7$ is poorly constrained. We present new constraints on

$\bar{x}_{\rm HI}$ at $z\sim5.1-6.8$, by analyzing deep optical spectra of $53$

quasars at $5.736.1$ are consistent with the redshift evolution of the neutral

fraction of the IGM derived from the Planck 2018.

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

The Statistics of Negative Power Spectrum Systematics in some 21 cm Analyses by Miguel F. Morales et al. on Wednesday 23 November

Through a very careful analysis Kolopanis et al. (2022) identified a negative

power spectrum (PS) systematic. The 21 cm cosmology community has assumed that

any observational systematics would add power, as negative PS are non-physical.

In addition to the mystery of their origin, negative PS systematics raise the

spectre of artificially lowering upper limits on the 21 cm PS. It appears that

the source of the negative PS systematics is a subtle interaction between

choices in how the PS estimate is calculated and baseline-dependent systematic

power. In this paper we present a statistical model of baseline dependent

systematics to explore how negative PS systematics can appear and their

statistical characteristics. This leads us to recommendations on when and how

to consider negative PS systematics when reporting observational 21 cm

cosmology upper limit.

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

Euclid: Modelling massive neutrinos in cosmology -- a code comparison by J. Adamek et al. on Wednesday 23 November

The measurement of the absolute neutrino mass scale from cosmological

large-scale clustering data is one of the key science goals of the Euclid

mission. Such a measurement relies on precise modelling of the impact of

neutrinos on structure formation, which can be studied with $N$-body

simulations. Here we present the results from a major code comparison effort to

establish the maturity and reliability of numerical methods for treating

massive neutrinos. The comparison includes eleven full $N$-body implementations

(not all of them independent), two $N$-body schemes with approximate time

integration, and four additional codes that directly predict or emulate the

matter power spectrum. Using a common set of initial data we quantify the

relative agreement on the nonlinear power spectrum of cold dark matter and

baryons and, for the $N$-body codes, also the relative agreement on the

bispectrum, halo mass function, and halo bias. We find that the different

numerical implementations produce fully consistent results. We can therefore be

confident that we can model the impact of massive neutrinos at the sub-percent

level in the most common summary statistics. We also provide a code validation

pipeline for future reference.

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