BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Baojiu Li

Baojiu Li contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.COgr-qcastro-ph.GAMachine Learning

Trust snapshot

Quick read

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

27 published item(s)

preprint2023arXiv

An emulator-based halo model in modified gravity -- I. The halo concentration-mass relation and density profile

In this series of papers we present an emulator-based halo model for the non-linear clustering of galaxies in modified gravity cosmologies. In the first paper, we present emulators for the following halo properties: the halo mass function, concentration-mass relation and halo-matter cross-correlation function. The emulators are trained on data extracted from the \textsc{FORGE} and \textsc{BRIDGE} suites of $N$-body simulations, respectively for two modified gravity (MG) theories: $f(R)$ gravity and the DGP model, varying three standard cosmological parameters $Ω_{\mathrm{m0}}, H_0, σ_8$, and one MG parameter, either $\bar{f}_{R0}$ or $r_{\mathrm{c}}$. Our halo property emulators achieve an accuracy of $\lesssim 1\%$ on independent test data sets. We demonstrate that the emulators can be combined with a galaxy-halo connection prescription to accurately predict the galaxy-galaxy and galaxy-matter correlation functions using the halo model framework.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Cosmological forecasts with the clustering of weak lensing peaks

Maximising the information that can be extracted from weak lensing measurements is a key goal for upcoming surveys such as LSST and Euclid. This is typically achieved through statistics that are complementary to the cosmic shear two-point correlation function, the most well established of which is the weak lensing peak abundance. In this work, we study the clustering of weak lensing peaks, and present parameter constraint forecasts for an LSST-like survey. We use the cosmoslics $w$CDM simulations to measure the peak two-point correlation function for a range of cosmological parameters, and use the simulation data to train a Gaussian process regression emulator which is applied to generate likelihood contours and provide parameter constraint forecasts from mock observations. We investigate the dependence of the peak two-point correlation function on the peak height, and find that the clustering of low amplitude peaks is complementary to that of high amplitude peaks. Consequently, their combination gives significantly tighter constraints than the clustering of high peaks alone. The peak two-point correlation function is significantly more sensitive to the cosmological parameters $h$ and $w_0$ than the peak abundance, and when the probes are combined, constraints on $Ω_{\rm m}$, $S_8$, $h$ and $w_0$ improve by at least a factor of two, relative to the peak abundance alone. Finally, we compare the forecasts for weak lensing peaks and weak lensing voids, and show that the two are also complementary; both probes can offer better constraints on $S_8$ and $w_0$ than the shear correlation function by roughly a factor of two.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fast full N-body simulations of generic modified gravity: conformal coupling models

We present MG-GLAM, a code developed for the very fast production of full $N$-body cosmological simulations in modified gravity (MG) models. We describe the implementation, numerical tests and first results of a large suite of cosmological simulations for three classes of MG models with conformal coupling terms: the $f(R)$ gravity, symmetron and coupled quintessence models. Derived from the parallel particle-mesh code GLAM, MG-GLAM incorporates an efficient multigrid relaxation technique to solve the characteristic nonlinear partial differential equations of these models. For $f(R)$ gravity, we have included new variants to diversify the model behaviour, and we have tailored the relaxation algorithms to these to maintain high computational efficiency. In a companion paper, we describe versions of this code developed for derivative coupling MG models, including the Vainshtein- and K-mouflage-type models. MG-GLAM can model the prototypes for most MG models of interest, and is broad and versatile. The code is highly optimised, with a tremendous speedup of a factor of more than a hundred compared with earlier $N$-body codes, while still giving accurate predictions of the matter power spectrum and dark matter halo abundance. MG-GLAM is ideal for the generation of large numbers of MG simulations that can be used in the construction of mock galaxy catalogues and the production of accurate emulators for ongoing and future galaxy surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fingerprints of modified gravity on galaxies in voids

We search for detectable signatures of f(R) gravity and its chameleon screening mechanism in the baryonic and dark matter (DM) properties of simulated void galaxies. The enhancement of the gravitational acceleration can have a meaningful impact on the scaling relations as well as on the halo morphology. The galaxy rotational velocity field (calculated with the velocity of the gas disc and the acceleration fields) deviates from the typical values of the Tully-Fisher Relation (TFR) in GR. For a given stellar mass, f(R) gravity tends to produce greater maximum velocities. On the other hand, the mass in haloes in f(R) gravity is more concentrated than their counterparts in GR. This trend changes when the concentration is calculated with the dynamical density profile, which takes into account the unscreened outer regions of the halo. Stellar discs interact with the overall potential well in the central regions, modifying the morphology of the screening regions and reshaping them. We find a trend for galaxies with a more dominant stellar disc to deviate further from round screening regions. We find that small haloes are less triaxial and more round in f(R) than their GR counterparts. The difference between halo morphology becomes smaller in f(R) haloes whose inner regions are screened. These results suggest possible observables that could unveil modified gravity effects on galaxies in voids in future cosmological tests of gravity.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Looking for a twist: probing the cosmological gravitomagnetic effect via weak lensing-kSZ cross correlations

General relativity predicts that the rotational momentum flux of matter twists the spacetime via a vector gravitomagnetic (frame-dragging) field, which remains undetected in cosmology. This vector field induces an additional gravitational lensing effect; at the same time, the momentum field sources the kinetic Sunyaev-Zel'dovich (kSZ) effect. The common origin of these two effects allows us to probe the gravitomagnetic signal via their cross-correlations. In this paper, we explore the possibility of detecting the gravitomagnetic field in $Λ$CDM by cross-correlating the weak-lensing convergence field with the CMB temperature map, which is imprinted with the kSZ signal. This approach allows us to extract the gravitomagnetic effect because the cross correlation between the standard Newtonian contribution to the weak-lensing convergence field, $κ_Φ$, and the kSZ effect is expected to vanish. We study the cross correlations with a suite of large-volume Newtonian $N$-body simulations and a small-volume, high-resolution, general-relativistic counterpart. We show that insufficient simulation resolution can introduce significant spurious correlations between $κ_Φ$ and kSZ. From the high-resolution simulation, we find that the cumulative signal-to-noise ratio (SNR) of the kSZ-gravitomagnetic convergence field can reach almost 15 (30) at $\ell\simeq5000$ ($10^4$) for the lensing source redshift $z_s=0.83$, if only cosmic variance is considered. We make forecast for next-generation lensing surveys such as Euclid and LSST, and CMB experiments such as Simons Observatory and CMB-S4, and find that, for $z_s=1.4$, the cumulative SNR can exceed 5 (9) at $\ell\simeq5000$ ($10^4$), indicating that the cosmological gravitomagnetic effect can be detected, if several foreground contaminations can be removed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nonlinear reconstruction of features in the primordial power spectrum from large-scale structure

Potential features in the primordial power spectrum have been searched for in galaxy surveys in recent years since these features can assist in understanding the nature of inflation. The null detection to date suggests that any such features should be fairly weak, and next-generation galaxy surveys, with their unprecedented sizes and precisions, are in a position to place stronger constraints than before. However, even if such primordial features once existed in the early Universe, they would have been significantly damped in the nonlinear regime at low redshift due to structure formation, which makes them difficult to be directly detected in real observations. A potential way to tackle this challenge for probing the features is to undo the cosmological evolution, i.e., using reconstruction to obtain an approximate linear density field. By employing a set of N-body simulations, here we show that a recently-proposed nonlinear reconstruction algorithm can effectively retrieve damped oscillatory features from halo catalogues and improve the accuracy of the measurement of feature parameters (assuming that such primordial features do exist). We do a Fisher analysis to forecast how nonlinear reconstruction affects the constraining power, and find that it can lead to significantly more robust constraints on the feature amplitude for a DESI-like survey. Comparing nonlinear reconstruction with other ways of improving constraints, such as increasing the survey volume and range of scales, this shows that it is possible to achieve what the latter do, but at a lower cost.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Spherical accretion of collisional gas in modified gravity I: self-similar solutions and a new cosmological hydrodynamical code

The spherical collapse scenario has great importance in cosmology since it captures several crucial aspects of structure formation. The presence of self-similar solutions in the Einstein-de Sitter (EdS) model greatly simplifies its analysis, making it a powerful tool to gain valuable insights into the real and more complicated physical processes involved in galaxy formation. While there has been a large body of research to incorporate various additional physical processes into spherical collapse, the effect of modified gravity (MG) models, which are popular alternatives to the $ΛCDM$ paradigm to explain the cosmic acceleration, is still not well understood in this scenario. In this paper, we study the spherical accretion of collisional gas in a particular MG model, which is a rare case that also admits self-similar solutions. The model displays interesting behaviours caused by the enhanced gravity and a screening mechanism. Despite the strong effects of MG, we find that its self-similar solution agrees well with that of the EdS model. These results are used to assess a new cosmological hydrodynamical code for spherical collapse simulations introduced here, which is based on the hyperbolic partial differential equation engine ExaHyPE 2. Its good agreement with the theoretical predictions confirms the reliability of this code in modelling astrophysical processes in spherical collapse. We will use this code to study the evolution of gas in more realistic MG models in future work.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The matter density PDF for modified gravity and dark energy with Large Deviations Theory

We present an analytical description of the probability distribution function (PDF) of the smoothed three-dimensional matter density field for modified gravity and dark energy. Our approach, based on the principles of Large Deviations Theory, is applicable to general extensions of the standard $Λ$CDM cosmology. We show that late-time changes to the law of gravity and background expansion can be included through Einstein-de Sitter spherical collapse dynamics combined with linear theory calculations and a calibration measurement of the non-linear variance of the smoothed density field from a simple numerical simulation. In a comparison to $N$-body simulations for $f(R)$, DGP and evolving dark energy theories, we find percent level accuracy around the peak of the distribution for predictions in the mildly non-linear regime. A Fisher forecast of an idealised experiment with a Euclid-like survey volume demonstrates the power of combining measurements of the 3D matter PDF with the 3D matter power spectrum. This combination is shown to halve the uncertainty on parameters for an evolving dark energy model, relative to a power spectrum analysis on its own. The PDF is also found to substantially increase the detection significance for small departures from General Relativity, with improvements of up to six times compared to the power spectrum alone. This analysis is therefore very promising for future studies including non-Gaussian statistics, as it has the potential to alleviate the reliance of these analyses on expensive high resolution simulations and emulators.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Towards an accurate model of small-scale redshift-space distortions in modified gravity

The coming generation of galaxy surveys will provide measurements of galaxy clustering with unprecedented accuracy and data size, which will allow us to test cosmological models at much higher precision than achievable previously. This means that we must have more accurate theoretical predictions to compare with future observational data. As a first step towards more accurate modelling of the redshift space distortions (RSD) of small-scale galaxy clustering in modified gravity (MG) cosmologies, we investigate the validity of the so-called Skew-T (ST) probability distribution function (PDF) of halo pairwise peculiar velocities in these models. We show that, combined with the streaming model of RSD, the ST PDF substantially improves the small-scale predictions by incorporating skewness and kurtosis, for both $Λ$CDM and two leading MG models: $f(R)$ gravity and the DGP braneworld model. The ST model reproduces the velocity PDF and redshift-space halo clustering measured from MG $N$-body simulations very well down to $\sim 5 \, h^{-1}\mathrm{Mpc}$. In particular, we investigate the enhancements of halo pairwise velocity moments with respect to $Λ$CDM for a larger range of MG variants than previous works, and present simple explanations to the behaviours observed. By performing a simple Fisher analysis, we find a significnat increase in constraining power to detect modifications of General Relativity by introducing small-scale information in the RSD analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Fast full $N$-body simulations of generic modified gravity: derivative coupling models

We present MG-GLAM, a code developed for the very fast production of full $N$-body cosmological simulations in modified gravity (MG) models. We describe the implementation, numerical tests and first results of a large suite of cosmological simulations for two broad classes of MG models with derivative coupling terms -- the Vainshtein- and Kmouflage-type models -- which respectively features the Vainshtein and Kmouflage screening mechanism. Derived from the parallel particle-mesh code GLAM, MG-GLAM incorporates an efficient multigrid relaxation technique to solve the characteristic nonlinear partial differential equations of these models. For Kmouflage, we have proposed a new algorithm for the relaxation solver, and run the first simulations of the model to understand its cosmological behaviour. In a companion paper, we describe versions of this code developed for conformally-coupled MG models, including several variants of $f(R)$ gravity, the symmetron model and coupled quintessence. Altogether, MG-GLAM has so far implemented the prototypes for most MG models of interest, and is broad and versatile. The code is highly optimised, with a tremendous (over two orders of magnitude) speedup when comparing its running time with earlier $N$-body codes, while still giving accurate predictions of the matter power spectrum and dark matter halo abundance. MG-GLAM is ideal for the generation of large numbers of MG simulations that can be used in the construction of mock galaxy catalogues and accurate emulators for ongoing and future galaxy surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

FORGE -- the f(R) gravity cosmic emulator project I: Introduction and matter power spectrum emulator

We present a large suite of cosmological simulations, the FORGE (F-of-R Gravity Emulator) simulation suite, which is designed to build accurate emulators for cosmological observables in galaxy clustering, weak gravitational lensing and galaxy clusters, for the $f(R)$ gravity model. The total of 200 simulations explore the cosmological parameter space around the Planck(2018) cosmology with a Latin hypercube, for 50 combinations of $\bar{f}_{R0}$, $Ω_m$, $σ_8$ and $h$ with all other parameters fixed. For each parameter combination, or node, we ran four independent simulations, one pair using $1024^3$ particles in $500 h^{-1} Mpc$ simulation boxes to cover small scales, and another pair using $512^3$ simulation particles in $1500 h^{-1} Mpc$ boxes for larger scales. Each pair of initial conditions are selected such that sample variance on large scales is minimised on average. In this work we present an accurate emulator for the matter power spectrum in $f(R)$ gravity trained on FORGE. We have verified, using the cross-validation technique, that the emulator accuracy is better than $2.5\%$ for the majority of nodes, particularly around the center of the explored parameter space, up to scales of $k = 10 h Mpc^{-1}$. We have also checked the power spectrum emulator against simulations which are not part of our training set and found excellent agreement. Due to its high accuracy on small scales, the FORGE matter power spectrum emulator is well suited for weak lensing analysis and can play a key tool in constraining $f(R)$ gravity using current and future observational data.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The impact of modified gravity on the Sunyaev-Zel'dovich effect

We study the effects of two popular modified gravity theories, which incorporate very different screening mechanisms, on the angular power spectra of the thermal (tSZ) and kinematic (kSZ) components of the Sunyaev-Zel'dovich effect. Using the first cosmological simulations that simultaneously incorporate both screened modified gravity and a complete galaxy formation model, we find that the tSZ and kSZ power spectra are significantly enhanced by the strengthened gravitational forces in Hu-Sawicki $f(R)$ gravity and the normal-branch Dvali-Gabadadze-Porrati model. Employing a combination of non-radiative and full-physics simulations, we find that the extra baryonic physics present in the latter acts to suppress the tSZ power on angular scales $l\gtrsim3000$ and the kSZ power on all tested scales, and this is found to have a substantial effect on the model differences. Our results indicate that the tSZ and kSZ power can be used as powerful probes of gravity on large scales, using data from current and upcoming surveys, provided sufficient work is conducted to understand the sensitivity of the constraints to baryonic processes that are currently not fully understood.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Novel Probes Project -- Tests of Gravity on Astrophysical Scales

We introduce The Novel Probes Project, an initiative to advance the field of astrophysical tests of the dark sector by creating a forum that connects observers and theorists. This review focuses on tests of gravity and is intended to be of use primarily to observers, but also to theorists with interest in the development of experimental tests. It is twinned with a separate review on tests of dark matter self-interactions (Adhikari et al., in prep.). Our focus is on astrophysical probes of gravity in the weak-field regime, ranging from stars to quasilinear cosmological scales. These are complementary to both strong-field tests and background and linear probes in cosmology. In particular, the nonlinear screening mechanisms that are an integral part of viable modified gravity models lead to characteristic signals specifically on astrophysical scales. The constraining power of these signals is not limited by cosmic variance, but comes with the challenge of building robust theoretical models of the nonlinear dynamics of stars, galaxies, clusters and large scale structure. In this review we lay the groundwork for a thorough exploration of the astrophysical regime with an eye to using the current and next generation of observations for tests of gravity. We begin by setting the scene for how theories beyond General Relativity are expected to behave, focusing primarily on screened fifth forces. We describe the analytic and numerical techniques for exploring the pertinent astrophysical systems, as well as the signatures of modified gravity. With these in hand we present a range of observational tests, and discuss prospects for future measurements and theoretical developments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Baryon acoustic oscillations reconstruction using convolutional neural networks

We propose a new scheme to reconstruct the baryon acoustic oscillations (BAO) signal, which contains key cosmological information, based on deep convolutional neural networks (CNN). Trained with almost no fine-tuning, the network can recover large-scale modes accurately in the test set: the correlation coefficient between the true and reconstructed initial conditions reaches $90\%$ at $k\leq 0.2 h\mathrm{Mpc}^{-1}$, which can lead to significant improvements of the BAO signal-to-noise ratio down to $k\simeq0.4h\mathrm{Mpc}^{-1}$. Since this new scheme is based on the configuration-space density field in sub-boxes, it is local and less affected by survey boundaries than the standard reconstruction method, as our tests confirm. We find that the network trained in one cosmology is able to reconstruct BAO peaks in the others, i.e. recovering information lost to non-linearity independent of cosmology. The accuracy of recovered BAO peak positions is far less than that caused by the difference in the cosmology models for training and testing, suggesting that different models can be distinguished efficiently in our scheme. It is very promising that Our scheme provides a different new way to extract the cosmological information from the ongoing and future large galaxy surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constraining structure formation using EDGES

The experiment to detect the global epoch of reionization signature (EDGES) collaboration reported the detection of a line at 78MHz in the sky-averaged spectrum due to neutral hydrogen (HI) 21-cm hyperfine absorption of cosmic microwave background photons at $z\sim 17$. This requires that the spin temperature of HI be coupled to the kinetic temperature of the gas at this $z$ through the scattering of Ly-$α$ photons emitted by massive stars. To explain the experimental result, star formation needs to be sufficiently efficient at $z\sim 17$ and this can be used to constrain models in which small-scale structure formation is suppressed (DMF models), either due to dark matter free-streaming or non-standard inflationary dynamics. We combine simulations of structure formation with a simple recipe for star formation to investigate whether these models emit enough Ly-$α$ photons to reproduce the experimental signal for reasonable values of the star formation efficiency, $f_\star$. We find that a thermal warm dark matter (WDM) model with mass $m_\mathrm{WDM}\sim 4.3\,\mathrm{keV}$ is consistent with the timing of the signal for $f_\star\lesssim 2\%$. The exponential growth of structure around $z\sim 17$ in such a model naturally generates a sharp onset of the absorption. A warmer model with $m_\mathrm{WDM}\sim 3\,\mathrm{keV}$ requires a higher star formation efficiency, $f_\star\sim 6\%$, which is a factor of few above predictions of current star formation models and observations of satellites in the Milky Way. However, uncertainties in the process of star formation at these $z$ do not allow to derive strong constrains on such models using 21-cm absorption line. The onset of the 21-cm absorption is generally faster in DMF compared to cold dark matter (CDM), unless some process significantly suppresses star formation in halos with masses below $\sim 10^8\,h^{-1}\,\mathrm{M}_\odot$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

GRAMSES: a new route to general relativistic $N$-body simulations in cosmology. Part I. Methodology and code description

We present GRAMSES, a new pipeline for nonlinear cosmological $N$-body simulations in General Relativity (GR). This code adopts the Arnowitt-Deser-Misner (ADM) formalism of GR, with constant mean curvature and minimum distortion gauge fixings, which provides a fully nonlinear and background independent framework for relativistic cosmology. Employing a fully constrained formulation, the Einstein equations are reduced to a set of ten elliptical equations which are solved using multigrid relaxation with adaptive mesh refinements (AMR), and three hyperbolic equations for the evolution of tensor degrees of freedom. The current version of GRAMSES neglects the latter by using the conformal flatness approximation, which allows it to compute the two scalar and two vector degrees of freedom of the metric. In this paper we describe the methodology, implementation, code tests and first results for cosmological simulations in a $Λ$CDM universe, while the generation of initial conditions and physical results will be discussed elsewhere. Inheriting the efficient AMR and massive parallelisation infrastructure from the publicly-available $N$-body and hydrodynamic simulation code RAMSES, GRAMSES is ideal for studying the detailed behaviour of spacetime inside virialised cosmic structures and hence accurately quantifying the impact of backreaction effects on the cosmic expansion, as well as for investigating GR effects on cosmological observables using cosmic-volume simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

GRAMSES: a new route to general relativistic $N$-body simulations in cosmology. Part II. Initial conditions

We address the generation of initial conditions (ICs) for GRAMSES, a code for nonlinear general relativistic (GR) $N$-body cosmological simulations recently introduced in Ref. [1]. GRAMSES adopts a constant mean curvature slicing with a minimal distortion gauge, where the linear growth rate is scale-dependent, and the standard method for realising initial particle data is not straightforwardly applicable. A new method is introduced, in which the initial positions of particles are generated from the displacement field realised for a matter power spectrum as usual, but the velocity is calculated by finite-differencing the displacement fields around the initial redshift. In this way, all the information required for setting up the initial conditions is drawn from three consecutive input matter power spectra, and additional assumptions such as scale-independence of the linear growth factor and growth rate are not needed. We implement this method in a modified 2LPTic code, and demonstrate that in a Newtonian setting it can reproduce the velocity field given by the default 2LPTic code with subpercent accuracy. We also show that the matter and velocity power spectra of the initial particle data generated for GRAMSES simulations using this method agree very well with the linear-theory predictions in the particular gauge used by GRAMSES. Finally, we discuss corrections to the finite difference calculation of the velocity when radiation is present, as well as additional corrections implemented in GRAMSES to ensure consistency. This method can be applied in ICs generation for GR simulations in generic gauges, and simulations of cosmological models with scale-dependent linear growth rate.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Iterative removal of redshift space distortions from galaxy clustering

Observations of galaxy clustering are made in redshift space, which results in distortions to the underlying isotropic distribution of galaxies. These redshift-space distortions (RSD) not only degrade important features of the matter density field, such as the baryonic acoustic oscillation (BAO) peaks, but also pose challenges for the theoretical modelling of observational probes. Here we introduce an iterative nonlinear reconstruction algorithm to remove RSD effects from galaxy clustering measurements, and assess its performance by using mock galaxy catalogues. The new method is found to be able to recover the real-space galaxy correlation function with an accuracy of $\sim1\%$, and restore the quadrupole accurately to $0$, on scales $s\gtrsim20\Mpch$. It also leads to an improvement in the reconstruction of the initial density field, which could help to accurately locate the BAO peaks. An `internal calibration' scheme is proposed to determine the values of cosmological parameters as a part of the reconstruction process, and possibilities to break parameter degeneracies are discussed. RSD reconstruction can offer a potential way to simultaneously extract the cosmological parameters, initial density field, real-space galaxy positions and large-scale peculiar velocity field (of the real Universe), making it an alternative to standard perturbative approaches in galaxy clustering analysis, bypassing the need for RSD modelling.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Nonlinear structure formation in Bound Dark Energy

We study nonlinear structure formation in the Bound Dark Energy model (BDE), where dark energy (DE) corresponds to a light scalar meson particle $ϕ$ dynamically formed at a condensation energy scale $Λ_c$. The evolution of this dark-energy meson is determined by the potential $V(ϕ)=Λ_c^{4+2/3}ϕ^{-2/3}$, with a distinguishing phenomenology from other quintessence scenarios. Particularly, the expansion rate of the universe is affected not only at late times, but also when the condensation of $ϕ$ occurs, which in linear theory leads to an enhancement (with respect to standard $Λ$CDM) of matter perturbations on small scales. We study how much of this signature is still present at late times as well as the properties of dark matter halos in the nonlinear regime through N-body simulations. Our results show that nonlinear corrections wash out this feature from the matter power spectrum even before DE becomes dominant. There is, however, a small but clear suppression of the BDE spectrum of $2\%$ today on the largest scales due to the distinct late-time dynamics of DE. The differences on the clustering power between BDE and $Λ$CDM are reflected in the halo mass function, where small halos are more abundant in BDE as opposed to large heavy structures, whose formation is delayed because of the expansion history of the universe. This result is well captured by the semi-analytical Sheth-Tormen formula. However, despite these differences, the halo concentration parameter is essentially the same in both models, which suggest that clustering inside the halos decouple from the general expansion once the halos form.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical solutions to Einstein's equations in a shearing-dust Universe: a code comparison

A number of codes for general-relativistic simulations of cosmological structure formation have been developed in recent years. Here we demonstrate that a sample of these codes produce consistent results beyond the Newtonian regime. We simulate solutions to Einstein's equations dominated by gravitomagnetism -- a vector-type gravitational field that doesn't exist in Newtonian gravity and produces frame-dragging, the leading-order post-Newtonian effect. We calculate the coordinate-invariant effect on intersecting null geodesics by performing ray tracing in each independent code. With this observable quantity, we assess and compare each code's ability to compute relativistic effects.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Observable tests of self-interacting dark matter in galaxy clusters: cosmological simulations with SIDM and baryons

We present BAHAMAS-SIDM, the first large-volume, (400/h Mpc)^3, cosmological simulations including both self-interacting dark matter (SIDM) and baryonic physics. These simulations are important for two primary reasons: 1) they include the effects of baryons on the dark matter distribution 2) the baryon particles can be used to make mock observables that can be compared directly with observations. As is well known, SIDM haloes are systematically less dense in their centres, and rounder, than CDM haloes. Here we find that that these changes are not reflected in the distribution of gas or stars within galaxy clusters, or in their X-ray luminosities. However, gravitational lensing observables can discriminate between DM models, and we present a menu of tests that future surveys could use to measure the SIDM interaction strength. We ray-trace our simulated galaxy clusters to produce strong lensing maps. Including baryons boosts the lensing strength of clusters that produce no critical curves in SIDM-only simulations. Comparing the Einstein radii of our simulated clusters with those observed in the CLASH survey, we find that at velocities around 1000 km/s an SIDM cross-section of sigma/m > 1 cm^2/g is likely incompatible with observed cluster lensing.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Accuracy of Weak Lensing Simulations

We investigate the accuracy of weak lensing simulations by comparing the results of five independently developed lensing simulation codes run on the same input $N$-body simulation. Our comparison focuses on the lensing convergence maps produced by the codes, and in particular on the corresponding PDFs, power spectra and peak counts. We find that the convergence power spectra of the lensing codes agree to $\lesssim 2\%$ out to scales $\ell \approx 4000$. For lensing peak counts, the agreement is better than $5\%$ for peaks with signal-to-noise $\lesssim 6$. We also discuss the systematic errors due to the Born approximation, line-of-sight discretization, particle noise and smoothing. The lensing codes tested deal in markedly different ways with these effects, but they nonetheless display a satisfactory level of agreement. Our results thus suggest that systematic errors due to the operation of existing lensing codes should be small. Moreover their impact on the convergence power spectra for a lensing simulation can be predicted given its numerical details, which may then serve as a validation test.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Towards a non-Gaussian model of redshift space distortions

To understand the nature of the accelerated expansion of the Universe, we need to combine constraints on the expansion rate and growth of structure. The growth rate is usually extracted from three dimensional galaxy maps by exploiting the effects of peculiar motions on galaxy clustering. However, theoretical models of the probability distribution function (PDF) of galaxy pairwise peculiar velocities are not accurate enough on small scales to reduce the error on theoretical predictions to the level required to match the precision expected for measurements from future surveys. Here, we improve the modelling of the pairwise velocity distribution by using the Skew-T PDF, which has nonzero skewness and kurtosis. Our model accurately reproduces the redshift-space multipoles (monopole, quadrupole and hexadecapole) predicted by N-body simulations, above scales of about $10\,h^{-1}{\rm Mpc}$. We illustrate how a Taylor expansion of the streaming model can reveal the contributions of the different moments to the clustering multipoles, which are independent of the shape of the velocity PDF. The Taylor expansion explains why the Gaussian streaming model works well in predicting the first two redshift-space multipoles, although the velocity PDF is non-Gaussian even on large scales. Indeed, any PDF with the correct first two moments would produce precise results for the monopole down to scales of about $10\,h^{-1}{\rm Mpc}$, and for the quadrupole down to about $30\,h^{-1}{\rm Mpc}$. An accurate model for the hexadecapole needs to include higher-order moments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Vector modes in $Λ$CDM: the gravitomagnetic potential in dark matter haloes from relativistic $N$-body simulations

We investigate the transverse modes of the gravitational and velocity fields in $Λ$CDM, based on a high-resolution simulation performed using the adaptive-mesh refinement general-relativistic $N$-body code GRAMSES. We study the generation of vorticity in the dark matter velocity field at low redshift, providing fits to the shape and evolution of its power spectrum over a range of scales. By analysing the gravitomagnetic vector potential, which is absent in Newtonian simulations, in dark matter haloes with masses ranging from $\sim10^{12.5}~h^{-1}{M}_{\odot}$ to $\sim10^{15}~h^{-1}{M}_{\odot}$, we find that its magnitude correlates with the halo mass, peaking in the inner regions. Nevertheless, on average, its ratio against the scalar gravitational potential remains fairly constant, below percent level, decreasing roughly linearly with redshift and showing a weak dependence on halo mass. Furthermore, we show that the gravitomagnetic acceleration in haloes peaks towards the core and reaches almost $10^{-10}$ $h$ cm/s$^2$ in the most massive halo of the simulation. However, regardless of the halo mass, the ratio between the magnitudes of the gravitomagnetic force and the standard gravitational force is typically at around the $10^{-5}$ level inside the haloes, again without significant radius dependence. This result confirms that the gravitomagnetic effects have a negligible impact on structure formation, even for the most massive structures, although its behaviour in low density regions remains to be explored. Likewise, the impact on observations remains to be understood in the future.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Emulators for the non-linear matter power spectrum beyond $Λ$CDM

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 $Λ$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 $Λ$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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Marked correlation functions in perturbation theory

We develop perturbation theory approaches to model the marked correlation function constructed to up-weight low density regions of the Universe, which might help distinguish modified gravity models from the standard cosmology model based on general relativity. Working within Convolution Lagrangian Perturbation Theory, we obtain that weighted correlation functions are expressible as double convolution integrals that we approximate using a combination of Eulerian and Lagrangian schemes. We find that different approaches agree within 1$\%$ on quasi non-linear scales. Compared with {\it N}-body simulations, the perturbation theory is found to provide accurate predictions for the marked correlation function of dark matter fields, dark matter halos as well as Halo Occupation Distribution galaxies down to $30$ Mpc/h. These analytic approaches help to understand the degeneracy between the mark and the galaxy bias and find a way to maximize the differences among various cosmological models.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Measuring the BAO peak position with different galaxy selections

We investigate if, for a fixed number density of targets and redshift, there is an optimal way to select a galaxy sample in order to measure the baryon acoustic oscillation (BAO) scale, which is used as a standard ruler to constrain the cosmic expansion. Using the mock galaxy catalogue built by Smith et al. in the Millennium-XXL N-body simulation with a technique to assign galaxies to dark matter haloes based on halo occupation distribution modelling, we consider the clustering of galaxies selected by luminosity, colour and local density. We assess how well the BAO scale can be extracted by fitting a template to the power spectrum measured for each sample. We find that the BAO peak position is recovered equally well for samples defined by luminosity or colour, while there is a bias in the BAO scale recovered for samples defined by density. The BAO position is contracted to smaller scales for the densest galaxy quartile and expanded to large scales for the two least dense galaxy quartiles. For fixed galaxy number density, density-selected samples have higher uncertainties in the recovered BAO scale than luminosity- or colour-selected samples.

0 citations0 reviewsNo code linkedOpen access