Astro arXiv | astro-ph.CO

Cosmography using strongly lensed gravitational waves from binary black holes by Souvik Jana et al. on Tuesday 22 November

Third generation gravitational wave (GW) detectors are expected to detect

millions of binary black hole (BBH) mergers during their operation period. A

small fraction of them ($\sim 1\%$) will be strongly lensed by intervening

galaxies and clusters, producing multiple observable copies of the GW signals.

The expected number of lensed events and the distribution of the time delay

between lensed events depend on the cosmology. We develop a Bayesian analysis

method for estimating cosmological parameters from the detected number of

lensed events and their time delay distribution. The expected constraints are

comparable to that obtained from other cosmological measurements, but probing a

different redshift regime ($z \sim 10$) that is not explored by other probes.

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

Starlet higher order statistics for galaxy clustering and weak lensing by Virginia Ajani et al. on Tuesday 22 November

We present a first application to photometric galaxy clustering and weak

lensing of wavelet based multi-scale higher order summary statistics: starlet

peak counts and starlet $\ell_1$-norm. Peak counts are the local maxima in the

map and the $\ell_1$-norm is computed via the sum of the absolute values of the

starlet (wavelet) decomposition coefficients of a map, providing a fast

multi-scale calculation of the pixel distribution, encoding the information of

all pixels in the map. We employ the cosmo-SLICS simulations sources and lenses

catalogues and we compute wavelet based higher order statistics in the context

of combined probes and their potential when applied to the weak lensing

convergence maps and galaxy maps. We get forecasts on the matter density

parameter $\Omega_{\rm m}$, the reduced Hubble constant $h$, the matter

fluctuation amplitude $\sigma_8$, and the dark energy equation of state

parameter $w_0$. We find that, in our setting for this first application,

considering the two probes as independent, starlet peaks and the $\ell_1$-norm

represent interesting summary statistics that can improve the constraints with

respect to the power spectrum also in the case of photometric galaxy clustering

and when the two probes are combined.

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

Investigating the Evolution of Amati Parameters with Redshift by Meghendra Singh et al. on Tuesday 22 November

Gamma Ray Bursts (GRB) are among the brightest objects in the Universe and

hence can be observed up to a very high redshift. Properly calibrated empirical

correlations between intensity and spectral correlations of GRBs can be used to

estimate the cosmological parameters. However, the possibility of the evolution

of GRBs with the redshift is a long-standing puzzle. In this work, we used 162

long-duration GRBs to determine whether GRBs below and above a certain redshift

have different properties. The GRBs are split into two groups, and we fit the

Amati relation for each group separately. Our findings demonstrate that

estimations of the Amati parameters for the two groups are substantially

dissimilar. We perform simulations to investigate whether the selection effects

could cause the difference. Our analysis shows that the differences may be

intrinsic, and the selection effects are not their true origin.

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

Consistent and simultaneous modelling of galaxy clustering and galaxy-galaxy lensing with Subhalo Abundance Matching by Sergio Contreras et al. on Tuesday 22 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

DESI Mock Challenge VII: halo and galaxy catalogs with the Bias Assignment Method by Andrés Balaguera-Antolínez et al. on Tuesday 22 November

This paper presents the generation of a suit of mock catalogs of dark matter

halos and galaxies, aiming at the assessment of precise covariance matrix for

cosmological analysis. The set of halo mock catalogs display accurate summary

statistics and detailed assignment of halo properties (e.g., velocity

dispersion, halo mass), enabling the possibility of using models of halo

occupation distribution (HOD) to construct mock catalogs with different galaxy

types. In particular, we generate mock catalogs based on HOD for emission-line

galaxies, key target for experiments such as the Dark Energy Spectroscopic

Instrument (DESI). To this end, we rely on the Bias Assignment Method (BAM), in

conjunction with the Scinet LIghtCone Simulations of halo mock catalogs (at

$z\sim 1$) as training set. We demonstrate the high accuracy of the mocks in

the one (abundance), two-(e.g. power spectrum), three- (e.g. bispectrum), four-

(covariance matrices of two point statistics), and six-point statistics

(covariance matrices of three-point statistics), both in real and redshift

space. BAM offers a robust way to produce fast and accurate halo distributions

which can be used to generate a variety of multi-tracer catalogs with precise

covariance matrices of a number of cosmological probes.

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

A bias-free cosmological analysis with quasars alleviating H 0 tension by Aleksander Łukasz Lenart et al. on Tuesday 22 November

Cosmological models and their parameters are widely debated because of

theoretical and observational mismatches of the standard cosmological model,

especially the current discrepancy between the value of the Hubble constant,

$H_{0}$, obtained by Type Ia supernovae (SNe Ia), and the Cosmic Microwave

Background Radiation (CMB). Thus, considering high-redshift probes like quasars

(QSOs), having intermediate redshifts between SNe Ia and CMB, is a necessary

step. In this work, we use SNe Ia and the most updated QSO sample, reaching

redshifts up to $z\sim7.5$, applying the Risaliti-Lusso QSO relation based on a

non-linear relation between ultraviolet and X-ray luminosities. We consider

this relation both in its original form and corrected for selection biases and

evolution in redshift through a reliable statistical method also accounting for

the circularity problem. We also explore two approaches: with and without

calibration on SNe Ia. We then investigate flat and non-flat standard

cosmological models and a flat $w$CDM model, with a constant dark energy

equation of state parameter $w$. Remarkably, when correcting for the evolution

as a function of cosmology, we obtain closed constraints on $\Omega_M$ using

only non-calibrated QSOs. We find that considering non-calibrated QSOs combined

with SNe Ia and accounting for the same correction, our results are compatible

with a flat $\Lambda$CDM model with $\Omega_M = 0.3$ and $H_0 = 70 \,

\mathrm{km\,s^{-1}\,Mpc^{-1}}$. Intriguingly, the $H_0$ values obtained place

halfway between the one from SNe Ia and CMB, paving the way for new insights

into the $H_0$ tension.

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

Quasar Factor Analysis -- An Unsupervised and Probabilistic Quasar Continuum Prediction Algorithm with Latent Factor Analysis by Zechang Sun et al. on Tuesday 22 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

The Statistics of Negative Power Spectrum Systematics in some 21 cm Analyses by Miguel F. Morales et al. on Tuesday 22 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

IWDM: The fate of an interacting non-cold dark matter - vacuum scenario by Supriya Pan et al. on Tuesday 22 November

In almost every cosmological models, the equation of state of the dark matter

is assumed to be zero (i.e. a pressure-less/cold dark matter). Although such

hypothesis is motivated by the abundance of cold dark matter in the universe,

there is however no compelling reason to set the dark matter equation of state

to zero, rather, the more generic picture is to consider a free-to-vary dark

matter equation of state and let the observational data decide its fate. With

the growing sensitivity of the experimental data, we choose the second

possibility and consider an interacting non-cold dark matter $-$ vacuum

scenario in which the dark matter equation of state is constant but

free-to-vary in an interval. Considering a very well known and most used

interaction function in the literature, we constrain this scenario using the

Cosmic Microwave Background (CMB) anisotropies and the CMB lensing

reconstruction from the legacy Planck release, baryon acoustic oscillations

distance measurements and the Pantheon catalogue from Supernovae Type Ia. We

find that for all the observational data sets, a non-zero value of the dark

matter equation of state is preferred at 68\% CL which indicates that a

non-cold dark matter sector in the universe should be investigated further in

order to understand the intrinsic nature of the dark matter sector.

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

Exhaustive Symbolic Regression by Deaglan J. Bartlett et al. on Tuesday 22 November

Symbolic Regression (SR) algorithms learn analytic expressions which both

accurately fit data and, unlike traditional machine-learning approaches, are

highly interpretable. Conventional SR suffers from two fundamental issues which

we address in this work. First, since the number of possible equations grows

exponentially with complexity, typical SR methods search the space

stochastically and hence do not necessarily find the best function. In many

cases, the target problems of SR are sufficiently simple that a brute-force

approach is not only feasible, but desirable. Second, the criteria used to

select the equation which optimally balances accuracy with simplicity have been

variable and poorly motivated. To address these issues we introduce a new

method for SR -- Exhaustive Symbolic Regression (ESR) -- which systematically

and efficiently considers all possible equations and is therefore guaranteed to

find not only the true optimum but also a complete function ranking. Utilising

the minimum description length principle, we introduce a principled method for

combining these preferences into a single objective statistic. To illustrate

the power of ESR we apply it to a catalogue of cosmic chronometers and the

Pantheon+ sample of supernovae to learn the Hubble rate as a function of

redshift, finding $\sim$40 functions (out of 5.2 million considered) that fit

the data more economically than the Friedmann equation. These low-redshift data

therefore do not necessarily prefer a $\Lambda$CDM expansion history, and

traditional SR algorithms that return only the Pareto-front, even if they found

this successfully, would not locate $\Lambda$CDM. We make our code and full

equation sets publicly available.

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