Source author record

Baojiu Li

Baojiu Li appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.CO gr-qc hep-ph hep-th astro-ph.GA astro-ph astro-ph.HE astro-ph.IM Machine Learning physics.comp-ph

Catalog footprint

What is connected

115works

10topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

preprint2016arXiv

An accurate method of modeling cluster scaling relations in modified gravity

We propose a new method to model cluster scaling relations in modified gravity. Using a suite of non-radiative hydrodynamical simulations, we show that the scaling relations of cumulative gas quantities, such as the Sunyaev Zel'dovich effect (Compton-y parameter) and the x-ray Compton-Y parameter, can be accurately predicted using the known results in the $Λ$CDM model with a precision of $\sim3\%$. This method provides a reliable way to analyze the gas physics in modified gravity using the less-demanding and much more efficient pure cold dark matter simulations. Our results therefore have important theoretical and practical implications in constraining gravity using cluster surveys.

preprint2016arXiv

Can background cosmology hold the key for modified gravity tests?

Modified gravity theories are a popular alternative to dark energy as a possible explanation for the observed accelerating cosmic expansion, and their cosmological tests are currently an active research field. Studies in recent years have been increasingly focused on testing these theories in the nonlinear regime, which is computationally demanding. Here we show that, under certain circumstances, a whole class of theories can be ruled out by using background cosmology alone. This is possible because certain classes of models (i) are fundamentally incapable of producing specific background expansion histories, and (ii) said histories are incompatible with local gravity tests. As an example, we demonstrate that a popular class of models, $f(R)$ gravity, would not be viable if observations suggest even a slight deviation of the background expansion history from that of the $Λ$CDM paradigm.

preprint2016arXiv

Constraining $f(R)$ Gravity Theory Using Weak Lensing Peak Statistics from the Canada-France-Hawaii-Telescope Lensing Survey

In this Letter, we report the observational constraints on the Hu-Sawicki $f(R)$ theory derived from weak lensing peak abundances, which are closely related to the mass function of massive halos. In comparison with studies using optical or x-ray clusters of galaxies, weak lensing peak analyses have the advantages of not relying on mass-baryonic observable calibrations. With observations from the Canada-France-Hawaii-Telescope Lensing Survey, our peak analyses give rise to a tight constraint on the model parameter $|f_{R0}|$ for $n=1$. The $95\%$ CL limit is $\log_{10}|f_{R0}| < -4.82$ given WMAP9 priors on $(Ω_{\rm m}, A_{\rm s})$. With Planck15 priors, the corresponding result is $\log_{10}|f_{R0}| < -5.16$.

preprint2016arXiv

Probing Theories of Gravity with Phase Space-Inferred Potentials of Galaxy Clusters

Modified theories of gravity provide us with a unique opportunity to generate innovative tests of gravity. In Chameleon f(R) gravity, the gravitational potential differs from the weak-field limit of general relativity (GR) in a mass dependent way. We develop a probe of gravity which compares high mass clusters, where Chameleon effects are weak, to low mass clusters, where the effects can be strong. We utilize the escape velocity edges in the radius/velocity phase space to infer the gravitational potential profiles on scales of 0.3-1 virial radii. We show that the escape edges of low mass clusters are enhanced compared to GR, where the magnitude of the difference depends on the background field value |fR0|. We validate our probe using N-body simulations and simulated light cone galaxy data. For a DESI (Dark Energy Spectroscopic Instrument) Bright Galaxy Sample, including observational systematics, projection effects, and cosmic variance, our test can differentiate between GR and Chameleon f(R) gravity models, |fR0| = 4e-6 (2e-6) at > 5σ (> 2σ), more than an order of magnitude better than current cluster-scale constraints.

preprint2016arXiv

RAY-RAMSES: a code for ray tracing on the fly in N-body simulations

We present a ray tracing code to compute integrated cosmological observables on the fly in AMR N-body simulations. Unlike conventional ray tracing techniques, our code takes full advantage of the time and spatial resolution attained by the N-body simulation by computing the integrals along the line of sight on a cell-by-cell basis through the AMR simulation grid. Moroever, since it runs on the fly in the N-body run, our code can produce maps of the desired observables without storing large (or any) amounts of data for post-processing. We implemented our routines in the RAMSES N-body code and tested the implementation using an example of weak lensing simulation. We analyse basic statistics of lensing convergence maps and find good agreement with semi-analytical methods. The ray tracing methodology presented here can be used in several cosmological analysis such as Sunyaev-Zel'dovich and integrated Sachs-Wolfe effect studies as well as modified gravity. Our code can also be used in cross-checks of the more conventional methods, which can be important in tests of theory systematics in preparation for upcoming large scale structure surveys.

preprint2016arXiv

Subhalo abundance matching in $f(R)$ gravity

Using the liminality N-body simulations of Shi et. al., we present the first predictions for galaxy clustering in $f(R)$ gravity using subhalo abundance matching. We find that, for a given galaxy density, even for an $f(R)$ model with $f_{R0}=-10^{-6}$, for which the cold dark matter clustering is very similar to $Λ$CDM, the predicted clustering of galaxies in the $f(R)$ model is very different from $Λ$CDM. The deviation can be as large as $40\%$ for samples with mean densities close to that of $L_*$ galaxies. This large deviation is testable given the accuracy that future large-scale galaxy surveys aim to achieve. Our result demonstrates that galaxy surveys can provide a stringent test of General Relativity on cosmological scales, which is comparable to the tests from local astrophysical observations.

preprint2016arXiv

Substructure and galaxy formation in the Copernicus Complexio warm dark matter simulations

We use the Copernicus Complexio (COCO) high resolution $N$-body simulations to investigate differences in the properties of small-scale structures in the standard cold dark matter (CDM) model and in a model with a cutoff in the initial power spectrum of density fluctuations consistent with both a thermally produced warm dark matter (WDM) particle or a sterile neutrino with mass 7 keV and leptogenesis parameter $L_6=8.7$. The latter corresponds to the "coldest" model with this sterile neutrino mass compatible with the identification of the recently detected 3.5 keV X-ray line as resulting from particle decay. CDM and WDM predict very different number densities of subhaloes with mass $\leq 10^9\,h^{-1}\,M_\odot$ although they predict similar, nearly universal, normalised subhalo radial density distributions. Haloes and subhaloes in both models have cuspy NFW profiles, but WDM subhaloes below the cutoff scale in the power spectrum (corresponding to maximum circular velocities $V_{\mathrm{max}}^{z=0} \leq50~\mathrm{kms}^{-1}$) are less concentrated than their CDM counterparts. We make predictions for observable properties using the GALFORM semi-analytic model of galaxy formation. Both models predict Milky Way satellite luminosity functions consistent with observations, although the WDM model predicts fewer very faint satellites. This model, however, predicts slightly more UV bright galaxies at redshift $z>7$ than CDM, but both are consistent with observations. Gravitational lensing offers the best prospect of distinguishing between the models.

preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$α$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

preprint2016arXiv

The DESI Experiment Part II: Instrument Design

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= λ/Δλ$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

preprint2016arXiv

The Distribution of Dark and Luminous Matter in the Unique Galaxy Cluster Merger Abell 2146

Abell 2146 ($z$ = 0.232) consists of two galaxy clusters undergoing a major merger. The system was discovered in previous work, where two large shock fronts were detected using the $\textit{Chandra X-ray Observatory}$, consistent with a merger close to the plane of the sky, caught soon after first core passage. A weak gravitational lensing analysis of the total gravitating mass in the system, using the distorted shapes of distant galaxies seen with ACS-WFC on $\textit{Hubble Space Telescope}$, is presented. The highest peak in the reconstruction of the projected mass is centred on the Brightest Cluster Galaxy (BCG) in Abell 2146-A. The mass associated with Abell 2146-B is more extended. Bootstrapped noise mass reconstructions show the mass peak in Abell 2146-A to be consistently centred on the BCG. Previous work showed that BCG-A appears to lag behind an X-ray cool core; although the peak of the mass reconstruction is centred on the BCG, it is also consistent with the X-ray peak given the resolution of the weak lensing mass map. The best-fit mass model with two components centred on the BCGs yields $M_{200}$ = 1.1$^{+0.3}_{-0.4}$$\times$10$^{15}$M$_{\odot}$ and 3$^{+1}_{-2}$$\times$10$^{14}$M$_{\odot}$ for Abell 2146-A and Abell 2146-B respectively, assuming a mass concentration parameter of $c=3.5$ for each cluster. From the weak lensing analysis, Abell 2146-A is the primary halo component, and the origin of the apparent discrepancy with the X-ray analysis where Abell 2146-B is the primary halo is being assessed using simulations of the merger.

preprint2015arXiv

$ν$Galileon: modified gravity with massive neutrinos as a testable alternative to $Λ$CDM

We show that, in the presence of massive neutrinos, the Galileon gravity model provides a very good fit to the current CMB temperature, CMB lensing and BAO data. This model, which we dub $ν \rm{Galileon}$, when assuming its stable attractor background solution, contains the same set of free parameters as $Λ\rm{CDM}$, although it leads to different expansion dynamics and nontrivial gravitational interactions. The data provide compelling evidence ($\gtrsim 6σ$) for nonzero neutrino masses, with $Σm_ν\gtrsim 0.4\ {\rm eV}$ at the $2σ$ level. Upcoming precision terrestrial measurements of the absolute neutrino mass scale therefore have the potential to test this model. We show that CMB lensing measurements at multipoles $l \lesssim 40$ will be able to discriminate between the $ν \rm{Galileon}$ and $Λ\rm{CDM}$ models. Unlike $Λ\rm{CDM}$, the $ν \rm{Galileon}$ model is consistent with local determinations of the Hubble parameter. The presence of massive neutrinos lowers the value of $σ_8$ substantially, despite of the enhanced gravitational strength on large scales. Unlike $Λ\rm{CDM}$, the $ν \rm{Galileon}$ model predicts a negative ISW effect, which is difficult to reconcile with current observational limits.

preprint2015arXiv

Distinguishing general relativity and $f(R)$ gravity with the gravitational lensing Minkowski functionals

We explore the Minkowski functionals of weak lensing convergence map to distinguish between $f(R)$ gravity and the general relativity (GR). The mock weak lensing convergence maps are constructed with a set of high-resolution simulations assuming different gravity models. It is shown that the lensing MFs of $f(R)$ gravity can be considerably different from that of GR because of the environmentally dependent enhancement of structure formation. We also investigate the effect of lensing noise on our results, and find that it is likely to distinguish F5, F6 and GR gravity models with a galaxy survey of $\sim3000$ degree$^2$ and with a background source number density of $n_g=30~{\rm arcmin}^{-2}$, comparable to an upcoming survey dark energy survey (DES). We also find that the $f(R)$ signal can be partially degenerate with the effect of changing cosmology, but combined use of other observations, such as the cosmic microwave background (CMB) data, can help break this degeneracy.

preprint2015arXiv

Effective Dark Matter Halo catalog in $f(R)$ gravity

We introduce the idea of {\it effective} dark matter halo catalog in $f(R)$ gravity, which is built using the {\it effective} density field. Using a suite of high resolution N-body simulations, we find that the dynamical properties of halos, such as the distribution of density, velocity dispersion, specific angular momentum and spin, in the effective catalog of $f(R)$ gravity closely mimic those in the $Λ$CDM model. Thus, when using effective halos, an $f(R)$ model can be viewed as a $Λ$CDM model. This effective catalog therefore provides a convenient way for studying the baryonic physics, the galaxy halo occupation distribution and even semi-analytical galaxy formation in $f(R)$ cosmologies.

preprint2015arXiv

Effective dark matter power spectra in $f(R)$ gravity

Using N-body simulations, we measure the power spectrum of the effective dark matter density field, which is defined through the modified Poisson equation in $f(R)$ cosmologies. We find that when compared to the conventional dark matter power spectrum, the effective power spectrum deviates more significantly from the $Λ$CDM model. For models with $f_{R0}=-10^{-4}$, the deviation can exceed 150\% while the deviation of the conventional matter power spectrum is less than 50\%. Even for models with $f_{R0}=-10^{-6}$, for which the conventional matter power spectrum is very close to the $Λ$CDM prediction, the effective power spectrum shows sizeable deviations. Our results indicate that traditional analyses based on the dark matter density field may seriously underestimate the impact of $f(R)$ gravity on galaxy clustering. We therefore suggest the use of the effective density field in such studies. In addition, based on our findings, we also discuss several possible methods of making use of the differences between the conventional and effective dark matter power spectra in $f(R)$ gravity to discriminate the theory from the $Λ$CDM model.

preprint2015arXiv

Exploring the liminality: properties of haloes and subhaloes in borderline $f(R)$ gravity

We investigate the properties of dark matter haloes and subhaloes in an $f(R)$ gravity model with $|f_{R0}|=10^{-6}$, using a very high-resolution N-body simulation. The model is a borderline between being cosmologically interesting and yet still consistent with current data. We find that the halo mass function in this model has a maximum 20% enhancement compared with the $Λ$CDM predictions between $z=1$ and $z=0$. Because of the chameleon mechanism which screens the deviation from standard gravity in dense environments, haloes more massive than $10^{13}h^{-1}M_\odot$ in this $f(R)$ model have very similar properties to haloes of similar mass in $Λ$CDM, while less massive haloes, such as that of the Milky Way, can have steeper inner density profiles and higher velocity dispersions due to their weaker screening. The halo concentration is remarkably enhanced for low-mass haloes in this model due to a deepening of the total gravitational potential. Contrary to the naive expectation, the halo formation time $z_f$ is later for low-mass haloes in this model, a consequence of these haloes growing faster than their counterparts in $Λ$CDM at late times and the definition of $z_f$. Subhaloes, especially those less massive than $10^{11}h^{-1}M_\odot$, are substantially more abundant in this $f(R)$ model for host haloes less massive than $10^{13}h^{-1}M_\odot$. We discuss the implications of these results for the Milky Way satellite abundance problem. Although the overall halo and subhalo properties in this borderline $f(R)$ model are close to their $Λ$CDM predictions, our results suggest that studies of the Local Group and astrophysical systems, aided by high-resolution simulations, can be valuable for further tests of it.

preprint2015arXiv

f(R) gravity on non-linear scales: The post-Friedmann expansion and the vector potential

Many modified gravity theories are under consideration in cosmology as the source of the accelerated expansion of the universe and linear perturbation theory, valid on the largest scales, has been examined in many of these models. However, smaller non-linear scales offer a richer phenomenology with which to constrain modified gravity theories. Here, we consider the Hu-Sawicki form of $f(R)$ gravity and apply the post-Friedmann approach to derive the leading order equations for non-linear scales, i.e. the equations valid in the Newtonian-like regime. We reproduce the standard equations for the scalar field, gravitational slip and the modified Poisson equation in a coherent framework. In addition, we derive the equation for the leading order correction to the Newtonian regime, the vector potential. We measure this vector potential from $f(R)$ N-body simulations at redshift zero and one, for two values of the $f_{R_0}$ parameter. We find that the vector potential at redshift zero in $f(R)$ gravity can be close to 50\% larger than in GR on small scales for $|f_{R_0}|=1.289\times10^{-5}$, although this is less for larger scales, earlier times and smaller values of the $f_{R_0}$ parameter. Similarly to in GR, the small amplitude of this vector potential suggests that the Newtonian approximation is highly accurate for $f(R)$ gravity, and also that the non-linear cosmological behaviour of $f(R)$ gravity can be completely described by just the scalar potentials and the $f(R)$ field.

preprint2015arXiv

K-mouflage gravity models that pass Solar System and cosmological constraints

We show that Solar System tests can place very strong constraints on K-mouflage models of gravity, which are coupled scalar field models with nontrivial kinetic terms that screen the fifth force in regions of large gravitational acceleration. In particular, the bounds on the anomalous perihelion of the Moon imposes stringent restrictions on the K-mouflage Lagrangian density, which can be met when the contributions of higher-order operators in the static regime are sufficiently small. The bound on the rate of change of the gravitational strength in the Solar System constrains the coupling strength $β$ to be smaller than $0.1$. These two bounds impose tighter constraints than the results from the Cassini satellite and Big Bang Nucleosynthesis. Despite the Solar System restrictions, we show that it is possible to construct viable models with interesting cosmological predictions. In particular, relative to $Λ$-CDM, such models predict percent-level deviations for the clustering of matter and the number density of dark matter haloes. This makes these models predictive and testable by forthcoming observational missions.

preprint2015arXiv

Linear perturbations in K-mouflage cosmologies with massive neutrinos

We present a comprehensive derivation of linear perturbation equations for different matter species, including photons, baryons, cold dark matter, scalar fields, massless and massive neutrinos, in the presence of a generic conformal coupling. Starting from the Lagrangians, we show how the conformal transformation affects the dynamics. In particular, we discuss how to incorporate consistently the scalar coupling in the equations of the Boltzmann hierarchy for massive neutrinos and the subsequent fluid approximations. We use the recently proposed K-mouflage model as an example to demonstrate the numerical implementation of our linear perturbation equations. K-mouflage is a new mechanism to suppress the fifth force between matter particles induced by the scalar coupling, but in the linear regime the fifth force is unsuppressed and can change the clustering of different matter species in different ways. We show how the CMB, lensing potential and matter power spectra are affected by the fifth force, and find ranges of K-mouflage parameters whose effects could be seen observationally. We also find that the scalar coupling can have the nontrivial effect of shifting the amplitude of the power spectra of the lensing potential and density fluctuations in opposite directions, although both probe the overall clustering of matter. This paper can serve as a reference for those who work on generic coupled scalar field cosmology, or those who are interested in the cosmological behaviour of the K-mouflage model.

preprint2015arXiv

Modified Gravity N-body Code Comparison Project

Self-consistent ${\it N}$-body simulations of modified gravity models are a key ingredient to obtain rigorous constraints on deviations from General Relativity using large-scale structure observations. This paper provides the first detailed comparison of the results of different ${\it N}$-body codes for the $f(R)$, DGP, and Symmetron models, starting from the same initial conditions. We find that the fractional deviation of the matter power spectrum from $Λ$CDM agrees to better than $1\%$ up to $k \sim 5-10~h/{\rm Mpc}$ between the different codes. These codes are thus able to meet the stringent accuracy requirements of upcoming observational surveys. All codes are also in good agreement in their results for the velocity divergence power spectrum, halo abundances and halo profiles. We also test the quasi-static limit, which is employed in most modified gravity ${\it N}$-body codes, for the Symmetron model for which the most significant non-static effects among the models considered are expected. We conclude that this limit is a very good approximation for all of the observables considered here.

preprint2015arXiv

Rapid simulation rescaling from standard to modified gravity models

We develop and test an algorithm to rescale a simulated dark-matter particle distribution or halo catalogue from a standard gravity model to that of a modified gravity model. This method is based on that of Angulo & White but with some additional ingredients to account for (i) scale-dependent growth of linear density perturbations and (ii) screening mechanisms that are generic features of viable modified gravity models. We attempt to keep the method as general as possible, so that it may plausibly be applied to a wide range of modified theories, although tests against simulations are restricted to a subclass of $f(R)$ models at this stage. We show that rescaling allows the power spectrum of matter to be reproduced at the $\sim 3$ per cent level in both real and redshift space up to $k=0.1h\,\mathrm{Mpc}^{-1}$ if we change the box size and alter the particle displacement field; this limit can be extended to $k=1h\,\mathrm{Mpc}^{-1}$ if we additionally alter halo internal structure. We simultaneously develop an algorithm that can be applied directly to a halo catalogue, in which case the halo mass function and clustering can be reproduced at the $\sim 5$ per cent level. Finally we investigate the clustering of halo particle distributions, generated from rescaled halo catalogues, and find that a similar accuracy can be reached.

preprint2015arXiv

Speeding up N-body simulations of modified gravity: Vainshtein screening models

We introduce and demonstrate the power of a method to speed up current iterative techniques for N-body modified gravity simulations. Our method is based on the observation that the accuracy of the final result is not compromised if the calculation of the fifth force becomes less accurate, but substantially faster, in high-density regions where it is weak due to screening. We focus on the nDGP model which employs Vainshtein screening, and test our method by running AMR simulations in which the solutions on the finer levels of the mesh (high density) are not obtained iteratively, but instead interpolated from coarser levels. We show that the impact this has on the matter power spectrum is below $1\%$ for $k < 5h/{\rm Mpc}$ at $z = 0$, and even smaller at higher redshift. The impact on halo properties is also small ($\lesssim 3\%$ for abundance, profiles, mass; and $\lesssim 0.05\%$ for positions and velocities). The method can boost the performance of modified gravity simulations by more than a factor of 10, which allows them to be pushed to resolution levels that were previously hard to achieve.

preprint2015arXiv

The COpernicus COmplexio: Statistical Properties of Warm Dark Matter Haloes

The recent detection of a 3.5 keV X-ray line from the centres of galaxies and clusters by Bulbul et al. (2014a) and Boyarsky et al. (2014a) has been interpreted as emission from the decay of 7 keV sterile neutrinos which could make up the (warm) dark matter (WDM). As part of the COpernicus COmplexio (COCO) programme, we investigate the properties of dark matter haloes formed in a high-resolution cosmological $N$-body simulation from initial conditions similar to those expected in a universe in which the dark matter consists of 7 keV sterile neutrinos. This simulation and its cold dark matter (CDM) counterpart have $\sim13.4$bn particles, each of mass $\sim 10^5\, h^{-1} M_\odot$, providing detailed information about halo structure and evolution down to dwarf galaxy mass scales. Non-linear structure formation on small scales ($M_{200}\, \leq\, 2 \times 10^9\,h^{-1}\,M_\odot$) begins slightly later in COCO-Warm than in COCO-Cold. The halo mass function at the present day in the WDM model begins to drop below its CDM counterpart at a mass $\sim 2 \times 10^{9}\,h^{-1}\,M_\odot$ and declines very rapidly towards lower masses so that there are five times fewer haloes of mass $M_{200}= 10^{8}\,h^{-1}\,M_\odot$ in COCO-Warm than in COCO-Cold. Halo concentrations on dwarf galaxy scales are correspondingly smaller in COCO-Warm, and we provide a simple functional form that describes its evolution with redshift. The shapes of haloes are similar in the two cases, but the smallest haloes in COCO-Warm rotate slightly more slowly than their CDM counterparts.

preprint2015arXiv

Using cosmic voids to distinguish f(R) gravity in future galaxy surveys

We use properties of void populations identified in $N$-body simulations to forecast the ability of upcoming galaxy surveys to differentiate models of f(R) gravity from \lcdm~cosmology. We analyze multiple simulation realizations, which were designed to mimic the expected number densities, volumes, and redshifts of the upcoming Euclid satellite and a lower-redshift ground-based counterpart survey, using the public {\tt VIDE} toolkit. We examine void abundances, ellipicities, radial density profiles, and radial velocity profiles at redshifts 1.0 and 0.43. We find that stronger f(R) coupling strengths eliminates small voids and produces voids up to $\sim 20\%$ larger in radius, leading to a significant tilt in the void number function. Additionally, under the influence of modified gravity, voids at all scales tend to be measurably emptier with correspondingly higher compensation walls. The velocity profiles reflect this, showing increased outflows inside voids and increased inflows outside voids. Using the void number function as an example, we forecast that future surveys can constrain the modified gravity coupling strength to $\sim 3 \times 10^{-5}$ using voids.

preprint2015arXiv

Weak lensing by voids in modified lensing potentials

We study lensing by voids in Cubic Galileon and Nonlocal gravity cosmologies, which are examples of theories of gravity that modify the lensing potential. We find voids in the dark matter and halo density fields of N-body simulations and compute their lensing signal analytically from the void density profiles, which we show are well fit by a simple analytical formula. In the Cubic Galileon model, the modifications to gravity inside voids are not screened and they approximately double the size of the lensing effects compared to GR. The difference is largely determined by the direct effects of the fifth force on lensing and less so by the modified density profiles. For this model, we also discuss the subtle impact on the force and lensing calculations caused by the screening effects of haloes that exist in and around voids. In the Nonlocal model, the impact of the modified density profiles and the direct modifications to lensing are comparable, but they boost the lensing signal by only $\approx 10\%$, compared with that of GR. Overall, our results suggest that lensing by voids is a promising tool to test models of gravity that modify lensing.

Baojiu Li

What is connected

Connect this record

See the researcher in context

Building this map preview

115 published item(s)

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

Cosmological forecasts with the clustering of weak lensing peaks

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

Fingerprints of modified gravity on galaxies in voids

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

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

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

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

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

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

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

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

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

Baryon acoustic oscillations reconstruction using convolutional neural networks

Constraining structure formation using EDGES

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

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

Iterative removal of redshift space distortions from galaxy clustering

Nonlinear structure formation in Bound Dark Energy

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

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

The Accuracy of Weak Lensing Simulations

Towards a non-Gaussian model of redshift space distortions

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

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

Marked correlation functions in perturbation theory

Measuring the BAO peak position with different galaxy selections

An accurate method of modeling cluster scaling relations in modified gravity

Can background cosmology hold the key for modified gravity tests?

Constraining $f(R)$ Gravity Theory Using Weak Lensing Peak Statistics from the Canada-France-Hawaii-Telescope Lensing Survey

Probing Theories of Gravity with Phase Space-Inferred Potentials of Galaxy Clusters

RAY-RAMSES: a code for ray tracing on the fly in N-body simulations

Subhalo abundance matching in $f(R)$ gravity

Substructure and galaxy formation in the Copernicus Complexio warm dark matter simulations

The DESI Experiment Part I: Science,Targeting, and Survey Design

The DESI Experiment Part II: Instrument Design

The Distribution of Dark and Luminous Matter in the Unique Galaxy Cluster Merger Abell 2146

$ν$Galileon: modified gravity with massive neutrinos as a testable alternative to $Λ$CDM

Distinguishing general relativity and $f(R)$ gravity with the gravitational lensing Minkowski functionals

Effective Dark Matter Halo catalog in $f(R)$ gravity

Effective dark matter power spectra in $f(R)$ gravity

Exploring the liminality: properties of haloes and subhaloes in borderline $f(R)$ gravity

f(R) gravity on non-linear scales: The post-Friedmann expansion and the vector potential

K-mouflage gravity models that pass Solar System and cosmological constraints

Linear perturbations in K-mouflage cosmologies with massive neutrinos

Modified Gravity N-body Code Comparison Project

Rapid simulation rescaling from standard to modified gravity models

Speeding up N-body simulations of modified gravity: Vainshtein screening models

The COpernicus COmplexio: Statistical Properties of Warm Dark Matter Haloes

Using cosmic voids to distinguish f(R) gravity in future galaxy surveys

Weak lensing by voids in modified lensing potentials

A clear and measurable signature of modified gravity in the galaxy velocity field

Chameleon $f(R)$ gravity on the Virgo cluster scale

Halo model and halo properties in Galileon gravity cosmologies

Halo modelling in chameleon theories

Massive Gravity wrapped in the Cosmic Web

Nonlinear structure formation in Nonlocal Gravity

Nonlinear structure formation in the Cubic Galileon gravity model

Revisiting the screening mechanism in $f(R)$ gravity

Spherical collapse in Galileon gravity: fifth force solutions, halo mass function and halo bias

Testing Gravity using Cosmic Voids

Testing Gravity using Void Profiles

Testing the quasi-static approximation in $f(R)$ gravity simulations

The Integrated Sachs-Wolfe effect in $f(R)$ gravity

The observational status of Galileon gravity after Planck

The Vainshtein Mechanism in the Cosmic Web

Void Dynamics

Voids in Modified Gravity Reloaded: Eulerian Void Assignment

Chameleon f(R) gravity in the virialized cluster

Exploring Vainshtein mechanism on adaptively refined meshes

Galaxy Infall Kinematics as a Test of Modified Gravity

Gravity theories, Transverse Doppler and Gravitational Redshifts in Galaxy Clusters

Hierarchical clustering in chameleon $f(R)$ gravity

Linear perturbations in Galileon gravity models