Search

We present TOPO (Time-Ordered Provable Outputs), a tool designed to enhance reproducibility and data integrity in astrophysical research, providing a trustless alternative to data analysis blinding. Astrophysical research frequently involves probabilistic algorithms, high computational demands, and stringent data privacy requirements, making it difficult to guarantee the integrity of results. TOPO provides a secure framework for verifying reproducible data analysis while ensuring sensitive information remains hidden. Our approach utilizes deterministic hashing to generate unique digital fingerprints of outputs, and Merkle Trees to store outputs in a time-ordered manner. This enables efficient verification of specific components or the entire dataset while making it computationally infeasible to manipulate results - thereby mitigating the risk of human interference and confirmation bias, key objectives of blinding methods. We demonstrate TOPO's utility in a cosmological context using TOPO-Cobaya, showcasing how blinding and verification can be seamlessly integrated into Markov Chain Monte Carlo calculations, rendering the process cryptographically provable. This method addresses pressing challenges in the reproducibility crisis and offers a robust, verifiable framework for astrophysical data analysis., Comment: 9 pages, 7 figures. Code available at: https://github.com/santiagocasas/topo-cobaya

We apply the Effective Field Theory of Large Scale Structure (EFTofLSS) to non-standard models of dark matter with suppressed small-scale structure imprinted by early-time physics, here exemplified by interacting dark matter (IDM) coupled to standard model neutrinos, and cross-check that the EFTofLSS has no trouble replicating the real-space halo-halo power spectrum from N-body simulations. We perform forecasts for a DESI ELG-like experiment using the redshift-space power spectrum and find that, under very conservative priors on these parameters, the EFTofLSS is not expected to yield strong constraints on dark matter interactions. However, with a better understanding of the evolution of counterterms and stochastic terms with redshift, realistic IDM models could in principle be detected using the full-shape power spectrum analysis of such a spectroscopic galaxy survey., Comment: 23 pages, 7 figures

We forecast constraints on minimal model-independent parametrisations of several Modified Gravity theories using mock Stage-IV cosmic shear data. We include nonlinear effects and screening, which ensures recovery of General Relativity on small scales. We introduce a power spectrum emulator to accelerate our analysis and evaluate the robustness of the growth index parametrisation with respect to two cosmologies: $\Lambda$CDM and the normal branch of the DGP model. We forecast the uncertainties on the growth index $\gamma$ to be of the order $\sim 10\%$. We find that our halo-model based screening approach demonstrates excellent performance, meeting the precision requirements of Stage-IV surveys. However, neglecting the screening transition results in biased predictions for cosmological parameters. We find that the screening transition shows significant degeneracy with baryonic feedback, requiring a much better understanding of baryonic physics for its detection. Massive neutrinos effects are less prominent and challenging to detect solely with cosmic shear data., Comment: 22 pages, 16 figures, 2 tables plus Appendices. Our power spectrum boost emulators for extended cosmologies are publicly available at https://github.com/nebblu/ReACT-emus/tree/main/emulators

To fully exploit the data from next generation surveys, we need an accurate modelling of the matter power spectrum up to non-linear scales. Therefore in this work we present the halo model reaction framework for the Generalized Cubic Covariant Galileon (GCCG) model, a modified gravity model within the Horndeski class of theories which extends the cubic covariant Galileon (G3) by including power laws of the derivatives of the scalar field in the K-essence and cubic terms. We modify the publicly available software ReACT for the GCCG in order to obtain an accurate prediction of the non-linear power spectrum. In the limit of the G3 model we compare the modified ReACT code to $N$-body simulations and we find agreement within 5\% for a wide range of scales and redshifts. We then study the relevant effects of the modifications introduced by the GCCG on the non-linear matter power spectrum. Finally, we provide forecasts from spectroscopic and photometric primary probes by next generation surveys using a Fisher matrix method. We show that future data will be able to constrain at 1$\sigma$ the two additional parameters of the model at the percent level and that considering non-linear corrections to the matter power spectrum beyond the linear regime is crucial to obtain this result., Comment: 17 pages, 7 figures, code available at https://github.com/nebblu/ACTio-ReACTio

The power spectrum is the most commonly applied summary statistics to extract cosmological information from the observed three-dimensional distribution of galaxies in spectroscopic surveys. We present CLASS-OneLoop, a new numerical tool, fully integrated into the Boltzmann code CLASS, enabling the calculation of the one-loop power spectrum of biased tracers in spectroscopic surveys. Built upon the Eulerian moment expansion framework for redshift-space distortions, the implemented model incorporates a complete set of nonlinear biases, counterterms, and stochastic contributions, and includes the infrared resummation and the Alcock-Paczynski effect. The code features an evaluation of the loops by either direct numerical integration or Fast Fourier Transform, and employs a fast-slow parameter decomposition, which is essential for accelerating MCMC runs. After presenting performance and validation tests, as an illustration of the capabilities of the code, we apply it to fit the measured redshift-space halo power spectrum wedges on a $\Lambda$CDM subset of the AbacusSummit simulation suite and considering scales up to $k_{\rm max} = 0.3\,h/$Mpc. We find that the one-loop model adeptly recovers the fiducial cosmology of the simulation, while a simplified model commonly used in the literature for sensitivity forecasts yields significantly biased results. Furthermore, we conduct Monte Carlo Markov Chain (MCMC) forecasts for a DESI-like survey, considering a model with a dynamical dark energy component. Our results demonstrate the ability to independently constrain cosmological and nuisance parameters, even in the presence of a large parameter space with twenty-nine variables.

We present jax-cosmo, a library for automatically differentiable cosmological theory calculations. It uses the JAX library, which has created a new coding ecosystem, especially in probabilistic programming. As well as batch acceleration, just-in-time compilation, and automatic optimization of code for different hardware modalities (CPU, GPU, TPU), JAX exposes an automatic differentiation (autodiff) mechanism. Thanks to autodiff, jax-cosmo gives access to the derivatives of cosmological likelihoods with respect to any of their parameters, and thus enables a range of powerful Bayesian inference algorithms, otherwise impractical in cosmology, such as Hamiltonian Monte Carlo and Variational Inference. In its initial release, jax-cosmo implements background evolution, linear and non-linear power spectra (using halofit or the Eisenstein and Hu transfer function), as well as angular power spectra with the Limber approximation for galaxy and weak lensing probes, all differentiable with respect to the cosmological parameters and their other inputs. We illustrate how autodiff can be a game-changer for common tasks involving Fisher matrix computations, or full posterior inference with gradient-based techniques. In particular, we show how Fisher matrices are now fast, exact, no longer require any fine tuning, and are themselves differentiable. Finally, using a Dark Energy Survey Year 1 3x2pt analysis as a benchmark, we demonstrate how jax-cosmo can be combined with Probabilistic Programming Languages to perform posterior inference with state-of-the-art algorithms including a No U-Turn Sampler, Automatic Differentiation Variational Inference,and Neural Transport HMC. We further demonstrate that Normalizing Flows using Neural Transport are a promising methodology for model validation in the early stages of analysis.

In this work we present updated forecasts on parameterised modifications of gravity that can capture deviations of the behaviour of cosmological density perturbations beyond $\Lambda$CDM. For these forecasts we adopt the SKA Observatory (SKAO) as a benchmark for future cosmological surveys at radio frequencies, combining a continuum survey for weak lensing and angular galaxy clustering with an HI galaxy survey for spectroscopic galaxy clustering that can detect baryon acoustic oscillations and redshift space distortions. Moreover, we also add 21cm HI intensity mapping, which provides invaluable information at higher redshifts, and can complement tomographic resolution, thus allowing us to probe redshift-dependent deviations of modified gravity models. For some of these cases, we combine the probes with other optical surveys, such as the Dark Energy Spectroscopic Instrument (DESI) and the Vera C. Rubin Observatory (VRO). We show that such synergies are powerful tools to remove systematic effects and degeneracies in the non-linear and small-scale modelling of the observables. Overall, we find that the combination of all SKAO radio probes will have the ability to constrain the present value of the functions parameterising deviations from $\Lambda$CDM ($\mu$ and $\Sigma$) with a precision of $2.7\%$ and $1.8\%$ respectively, competitive with the constraints expected from optical surveys and with constraints we have on gravitational interactions in the standard model. Exploring the radio-optical synergies, we find that the combination of VRO with SKAO can yield extremely tight constraints on $\mu$ and $\Sigma$ ($0.9\%$ and $0.7\%$ respectively), which are further improved when the cross-correlation between intensity mapping and DESI galaxies is included., Comment: 24 pages, 10 figures

Long range scalar fields with a coupling to matter appear to violate known bounds on gravitation in the solar system and the laboratory. This is evaded thanks to screening mechanisms. In this short review, we shall present the various screening mechanisms from an effective field theory point of view. We then investigate how they can and will be tested in the laboratory and on astrophysical and cosmological scales., Comment: Invited review published in Universe 8 (2021), 49 pages, 10 figures

In this review, we outline the expected tests of gravity that will be achieved at cosmological scales in the upcoming decades. We focus mainly on constraints on phenomenologically parameterized deviations from general relativity, which allow to test gravity in a model-independent way, but also review some of the expected constraints obtained with more physically motivated approaches. After reviewing the state-of-the-art for such constraints, we outline the expected improvement that future cosmological surveys will achieve, focusing mainly on future large-scale structures and cosmic microwave background surveys but also looking into novel probes on the nature of gravity. We will also highlight the necessity of overcoming accuracy issues in our theoretical predictions, issues that become relevant due to the expected sensitivity of future experiments., Comment: 29 pages, 10 figures. Version matching the one published by Universe on the Special Issue "Large Scale Structure of the Universe"

We present the second edition of the book `CosmoAmautas: Astrophysics in the classroom' which is a text and activity book on selected topics in astrophysics. The content is written in Spanish and intended as the introductory literature and material for high school teachers who want to integrate astronomy to their science and math classes. The book includes an introduction to topics such as planetary science, stellar evolution, exoplanets, galaxies, black holes, and cosmology. In addition to the astronomy content, the book includes didactic notes on general physics of gravity and light, unit conversion, history of science, and tools for the digital classroom (ICT tools), which are compatible with the requirements of the Peruvian high school curriculum. This book is constructed based on the astrophysics curriculum and education model introduced in the online teacher-training workshop CosmoAmautas, which took place in 2021 and 2022. CosmoAmautas (amauta = teacher in the Quechua language) is an astronomy education project in Peru, funded by the IAU Office of Astronomy for Development. For more information, please visit www.cosmoamautas.org., Comment: Second Edition (2022). 164 pages, in Spanish. We acknowledge the support of the IAU Office of Astronomy for Development, which is a joint project of the International Astronomical Union (IAU) and the South African National Research Foundation (NRF) with the support of the Department of Science and Innovation (DSI)

Single-field models of $\alpha$-attractor quintessential inflation provide a unified picture of the two periods of early- and late-time cosmic acceleration, where both inflation and dark energy are described by a single scalar degree of freedom rolling down a runaway potential. These theoretically well-motivated models have distinct observational predictions that are in agreement with existing cosmological data. We show that the next generation of large-scale structure surveys, even when no other cosmological data sets are considered, will strongly constrain the parameter space of these models, and test them against the standard cosmological model and more conventional non-quintessential inflation. In particular, we expect $\mathcal{O}(10^{-5}\mathrm{-}10^{-4})$ constraints on the present values of the dark energy equation of state and its time derivative, $w_0$ and $w_a$. We also forecast more than one order of magnitude tighter constraints on the spectral index of primordial curvature perturbations $n_s$ compared to the expectations for the standard model. This demonstrates the powerful synergy between the upcoming large-scale structure probes of inflation and those aiming to measure the tensor-to-scalar ratio $r$ through the observation of $B$-mode polarization of the cosmic microwave background., Comment: 34 pages, 14 figures. v2: Improved and published version; typos fixed, references and discussions added, figures 1 and 8 added

Constraints on gravity and cosmology will greatly benefit from performing joint clustering and weak lensing analyses on large-scale structure data sets. Utilising non-linear information coming from small physical scales can greatly enhance these constraints. At the heart of these analyses is the matter power spectrum. Here we employ a simple method, dubbed "Hybrid $P_\ell(k)$", based on the Gaussian Streaming Model (GSM), to calculate the quasi non-linear redshift space matter power spectrum multipoles. This employs a fully non-linear and theoretically general prescription for the matter power spectrum. We test this approach against comoving Lagrangian acceleration simulation measurements performed in GR, DGP and $f(R)$ gravity and find that our method performs comparably or better to the dark matter TNS redshift space power spectrum model for dark matter. When comparing the redshift space multipoles for halos, we find that the Gaussian approximation of the GSM with a linear bias and a free stochastic term, $N$, is competitive to the TNS model. Our approach offers many avenues for improvement in accuracy as well as further unification under the halo model., Comment: 34 pages, 16 figures, JCAP accepted version. HyPk code available at: https://github.com/nebblu/HyPk

Accurate predictions for the non-linear matter power spectrum are needed to confront theory with observations in current and near future weak lensing and galaxy clustering surveys. We propose a computationally cheap method to create an emulator for modified gravity models by utilizing existing emulators for $\Lambda$CDM. Using a suite of $N$-body simulations we construct a fitting function for the enhancement of both the linear and non-linear matter power spectrum in the commonly studied Hu-Sawicki $f(R)$ gravity model valid for wave-numbers $k \lesssim 5-10\, h\text{Mpc}^{-1}$ and redshifts $z \lesssim 3$. We show that the cosmology dependence of this enhancement is relatively weak so that our fit, using simulations coming from only one cosmology, can be used to get accurate predictions for other cosmological parameters. We also show that the cosmology dependence can, if needed, be included by using linear theory, approximate $N$-body simulations (such as COLA) and semi-analytical tools like the halo model. Our final fit can easily be combined with any emulator or semi-analytical models for the non-linear $\Lambda$CDM power spectrum to accurately, and quickly, produce a non-linear power spectrum for this particular modified gravity model. The method we use can be applied to fairly cheaply construct an emulator for other modified gravity models. As an application of our fitting formula we use it to compute Fisher-forecasts for how well galaxy clustering and weak lensing in a Euclid-like survey will be at constraining modifications of gravity., Comment: 12 pages, 10 figures. Version accepted for publication in PRD. Data can be found at https://github.com/HAWinther/FofrFittingFunction

This paper is a pedagogical introduction to models of gravity and how to constrain them through cosmological observations. We focus on the Horndeski scalar-tensor theory and on the quantities that can be measured with a minimum of assumptions. Alternatives or extensions of General Relativity have been proposed ever since its early years. Because of Lovelock theorem, modifying gravity in four dimensions typically means adding new degrees of freedom. The simplest way is to include a scalar field coupled to the curvature tensor terms. The most general way of doing so without incurring in the Ostrogradski instability is the Horndeski Lagrangian and its extensions. Testing gravity means therefore, in its simplest term, testing the Horndeski Lagrangian. Since local gravity experiments can always be evaded by assuming some screening mechanism or that baryons are decoupled, or even that the effects of modified gravity are visible only at early times, we need to test gravity with cosmological observations in the late universe (large-scale structure) and in the early universe (cosmic microwave background). In this work we review the basic tools to test gravity at cosmological scales, focusing on model-independent measurements., Comment: Review paper. Clarified some aspects and added new references

The gravitational slip parameter is an important discriminator between large classes of gravity theories at cosmological and astrophysical scales. In this work we use a combination of simulated information of galaxy cluster mass profiles, inferred by Strong+Weak lensing analyses and by the study of the dynamics of the cluster member galaxies, to reconstruct the gravitational slip parameter $\eta$ and predict the accuracy with which it can be constrained with current and future galaxy cluster surveys. Performing a full-likelihood statistical analysis, we show that galaxy cluster observations can constrain $\eta$ down to the percent level already with a few tens of clusters. We discuss the significance of possible systematics, and show that the cluster masses and numbers of galaxy members used to reconstruct the dynamics mass profile have a mild effect on the predicted constraints., Comment: 18 pages plus references, 8 Figures

We discuss the cosmological phenomenology of biscalar-tensor models displaying a maximally symmetric Einstein-frame kinetic sector and constructed on the basis of scale symmetry and volume-preserving diffeomorphisms. These theories contain a single dimensionful parameter $\Lambda_0$-associated with the invariance under the aforementioned restricted coordinate transformations-and a massless dilaton field. At large field values these scenarios lead to inflation with no generation of isocurvature perturbations. The corresponding predictions depend only on two dimensionless parameters, which characterize the curvature of the field--manifold and the leading order behavior of the inflationary potential. For $\Lambda_0=0$ the scale symmetry is unbroken and the dilaton admits only derivative couplings to matter, evading all fifth force constraints. For $\Lambda_0\neq 0$ the field acquires a run-away potential that can support a dark energy dominated era at late times. We confront a minimalistic realization of this appealing framework with present data sets. The impact of possible consistency relations among the early and late Universe dynamics that could appear within this setting is also discussed., Comment: Journal version---discussion extended, typos corrected, references added

In the context of the effective field theory of dark energy (EFT) we perform agnostic explorations of Horndeski gravity. We choose two parametrizations for the free EFT functions, namely a power law and a dark energy density-like behaviour on a non trivial Chevallier-Polarski-Linder background. We restrict our analysis to those EFT functions which do not modify the speed of propagation of gravitational waves. Among those, we prove that one specific function cannot be constrained by data, since its contribution to the observables is below the cosmic variance, although we show it has a relevant role in defining the viable parameter space. We place constraints on the parameters of these models combining measurements from present day cosmological datasets and we prove that the next generation galaxy surveys can improve such constraints by one order of magnitude. We then verify the validity of the quasi-static limit within the sound horizon of the dark field, by looking at the phenomenological functions $\mu$ and $\Sigma$, associated respectively to clustering and lensing potentials. Furthermore, we notice up to $5\%$ deviations in $\mu, \Sigma$ with respect to General Relativity at scales smaller than the Compton one. For the chosen parametrizations and in the quasi-static limit, future constraints on $\mu$ and $\Sigma$ can reach the $1\%$ level and will allow us to discriminate between certain models at more than $3\sigma$, provided the present best-fit values remain., Comment: 12 pages, 8 figures

In this work, we use recent data on the Hubble expansion rate $H(z)$, the quantity $f\sigma_8(z)$ from redshift space distortions and the statistic $E_g$ from clustering and lensing observables to constrain in a model-independent way the linear anisotropic stress parameter $\eta$. This estimate is free of assumptions about initial conditions, bias, the abundance of dark matter and the background expansion. We denote this observable estimator as $\eta_{{\rm obs}}$. If $\eta_{{\rm obs}}$ turns out to be different from unity, it would imply either a modification of gravity or a non-perfect fluid form of dark energy clustering at sub-horizon scales. Using three different methods to reconstruct the underlying model from data, we report the value of $\eta_{{\rm obs}}$ at three redshift values, $z=0.29, 0.58, 0.86$. Using the method of polynomial regression, we find $\eta_{{\rm obs}}=0.57\pm1.05$, $\eta_{{\rm obs}}=0.48\pm0.96$, and $\eta_{{\rm obs}}=-0.11\pm3.21$, respectively. Assuming a constant $\eta_{{\rm obs}}$ in this range, we find $\eta_{{\rm obs}}=0.49\pm0.69$. We consider this method as our fiducial result, for reasons clarified in the text. The other two methods give for a constant anisotropic stress $\eta_{{\rm obs}}=0.15\pm0.27$ (binning) and $\eta_{{\rm obs}}=0.53 \pm 0.19$ (Gaussian Process). We find that all three estimates are compatible with each other within their $1\sigma$ error bars. While the polynomial regression method is compatible with standard gravity, the other two methods are in tension with it., Comment: 25 pages, 4 figures, 9 tables

The Higgs-Dilaton model is a scale-invariant extension of the Standard Model non-minimally coupled to gravity and containing just one additional degree of freedom on top of the Standard Model particle content. This minimalistic scenario predicts a set of measurable consistency relations between the inflationary observables and the dark-energy equation-of-state parameter. We present an alternative derivation of these consistency relations that highlights the connections and differences with the $\alpha$-attractor scenario. We study in how far these constraints allow to distinguish the Higgs-Dilaton model from $\Lambda$CDM and $w$CDM cosmologies. To this end we first analyze existing data sets using a Markov Chain Monte Carlo approach. Second, we perform forecasts for future galaxy surveys using a Fisher matrix approach, both for galaxy clustering and weak lensing probes. Assuming that the best fit values in the different models remain comparable to the present ones, we show that both Euclid- and SKA2-like missions will be able to discriminate a Higgs-Dilaton cosmology from $\Lambda$CDM and $w$CDM., Comment: 17 pages, 7 figures, added forecasts for DESI - matches the published version

Modified Gravity theories generally affect the Poisson equation and the gravitational slip (effective anisotropic stress) in an observable way, that can be parameterized by two generic functions ($\eta$ and $\mu$) of time and space. We bin the time dependence of these functions in redshift and present forecasts on each bin for future surveys like Euclid. We consider both Galaxy Clustering and Weak Lensing surveys, showing the impact of the non-linear regime, treated with two different semi-analytical approximations. In addition to these future observables, we use a prior covariance matrix derived from the Planck observations of the Cosmic Microwave Background. Our results show that $\eta$ and $\mu$ in different redshift bins are significantly correlated, but including non-linear scales reduces or even eliminates the correlation, breaking the degeneracy between Modified Gravity parameters and the overall amplitude of the matter power spectrum. We further decorrelate parameters with a Zero-phase Component Analysis and identify which combinations of the Modified Gravity parameter amplitudes, in different redshift bins, are best constrained by future surveys. We also extend the analysis to two particular parameterizations of the time evolution of $\mu$ and $\eta$ and consider, in addition to Euclid, also SKA1, SKA2, DESI: we find in this case that future surveys will be able to constrain the current values of $\eta$ and $\mu$ at the $2-5\%$ level when using only linear scales (wavevector k < 0.15 h/Mpc), depending on the specific time parameterization; sensitivity improves to about $1\%$ when non-linearities are included., Comment: 45 pages, 15 tables, 25 figures

We investigate the formation and dissipation of large scale neutrino structures in cosmologies where the time evolution of dynamical dark energy is stopped by a growing neutrino mass. In models where the coupling between neutrinos and dark energy grows with the value of the scalar cosmon field, the evolution of neutrino lumps depends on the neutrino mass. For small masses the lumps form and dissolve periodically, leaving only a small backreaction of the neutrino structures on the cosmic evolution. This process heats the neutrinos to temperatures much above the photon temperature such that neutrinos acquire again an almost relativistic equation of state. The present equation of state of the combined cosmon-neutrino fluid is very close to -1. In contrast, for larger neutrino masses the lumps become stable. The highly concentrated neutrino structures entail a large backreaction similar to the case of a constant neutrino-cosmon coupling. A present average neutrino mass of around 0.5 eV seems so far compatible with observation. For masses lower than this value, neutrino induced gravitational potentials remain small, making the lumps difficult to detect., Comment: 19 pages, 10 figures, 3 tables

We consider cosmological models in which dark matter feels a fifth force mediated by the dark energy scalar field, also known as coupled dark energy. Our interest resides in estimating forecasts for future surveys like Euclid when we take into account non-linear effects, relying on new fitting functions that reproduce the non-linear matter power spectrum obtained from N-body simulations. We obtain fitting functions for models in which the dark matter-dark energy coupling is constant. Their validity is demonstrated for all available simulations in the redshift range $z=0-1.6$ and wave modes below $k=10 \text{h/Mpc}$. These fitting formulas can be used to test the predictions of the model in the non-linear regime without the need for additional computing-intensive N-body simulations. We then use these fitting functions to perform forecasts on the constraining power that future galaxy-redshift surveys like Euclid will have on the coupling parameter, using the Fisher matrix method for galaxy clustering (GC) and weak lensing (WL). We find that by using information in the non-linear power spectrum, and combining the GC and WL probes, we can constrain the dark matter-dark energy coupling constant squared, $\beta^{2}$, with precision smaller than 4\% and all other cosmological parameters better than 1\%, which is a considerable improvement of more than an order of magnitude compared to corresponding linear power spectrum forecasts with the same survey specifications., Comment: 21 pages, 8 figures. New section about systematic bias on comsological parameters due to uncertainties in the non-linear power spectrum