BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Alexandre Barreira

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

ResearcherUnclaimed source record
astro-ph.COgr-qchep-phastro-ph.GAhep-th

Catalog footprint

What is connected

27works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

27 published item(s)

preprint2023arXiv

Assembly bias in the local PNG halo bias and its implication for $f_{\rm NL}$ constraints

We use $N$-body simulations to study halo assembly bias (i.e., the dependence of halo clustering on properties beyond total mass) in the density and primordial non-Gaussianity (PNG) linear bias parameters $b_1$ and $b_ϕ$, respectively. We consider concentration, spin and sphericity as secondary halo properties, for which we find a clear detection of assembly bias for $b_1$ and $b_ϕ$. At fixed total mass, halo spin and sphericity impact $b_1$ and $b_ϕ$ in a similar manner, roughly preserving the shape of the linear $b_ϕ(b_1)$ relation satisfied by the global halo population. Halo concentration, however, drives $b_1$ and $b_ϕ$ in opposite directions. This induces significant changes to the $b_ϕ(b_1)$ relation, with higher concentration halos having higher amplitude of $b_ϕ(b_1)$. For $z=0.5$ and $b_1 \approx 2$ in particular, the population comprising either all halos, those with the $33\%$ lowest or those with the $33\%$ highest concentrations have a PNG bias of $b_ϕ\approx 3$, $b_ϕ\approx -1$ and $b_ϕ\approx 9$, respectively. Varying the halo concentration can make $b_ϕ$ very small and even change its sign. These results have important ramifications for galaxy clustering constraints of the local PNG parameter $f_{\rm NL}$ that assume fixed forms for the $b_ϕ(b_1)$ relation. We illustrate the significant impact of halo assembly bias in actual data using the BOSS DR12 galaxy power spectrum: assuming that BOSS galaxies are representative of all halos, the $33\%$ lowest or the $33\%$ highest concentration halos yields $σ_{f_{\rm NL}} = 44, 165, 19$, respectively. Our results suggest taking host halo concentration into account in galaxy selection strategies to maximize the signal-to-noise on $f_{\rm NL}$. They also motivate more simulation-based efforts to study the $b_ϕ(b_1)$ relation of halos and galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Predictions for local PNG bias in the galaxy power spectrum and bispectrum and the consequences for $f_{\rm NL}$ constraints

We use hydrodynamical separate universe simulations with the IllustrisTNG model to predict the local primordial non-Gaussianity (PNG) bias parameters $b_ϕ$ and $b_{ϕδ}$, which enter at leading order in the galaxy power spectrum and bispectrum. This is the first time that $b_{ϕδ}$ is measured from either gravity-only or galaxy formation simulations. For dark matter halos, the popular assumption of universality overpredicts the $b_{ϕδ}(b_1)$ relation in the range $1 \lesssim b_1 \lesssim 3$ by up to $Δb_{ϕδ} \sim 3$ ($b_1$ is the linear density bias). The adequacy of the universality relation is worse for the simulated galaxies, with the relations $b_ϕ(b_1)$ and $b_{ϕδ}(b_1)$ being generically redshift-dependent and very sensitive to how galaxies are selected (we test total, stellar and black hole mass, black hole mass accretion rate and color). The uncertainties on $b_ϕ$ and $b_{ϕδ}$ have a direct, often overlooked impact on the constraints of the local PNG parameter $f_{\rm NL}$, which we study and discuss. For a survey with $V = 100{\rm Gpc}^3/h^3$ at $z=1$, uncertainties $Δb_ϕ \lesssim 1$ and $Δb_{ϕδ} \lesssim 5$ around values close to the fiducial can yield relatively unbiased constraints on $f_{\rm NL}$ using power spectrum and bispectrum data. We also show why priors on galaxy bias are useful even in analyses that fit for products $f_{\rm NL} b_ϕ$ and $f_{\rm NL} b_{ϕδ}$. The strategies we discuss to deal with galaxy bias uncertainties can be straightforwardly implemented in existing $f_{\rm NL}$ constraint analyses (we provide fits for some of the bias relations). Our results motivate more works with galaxy formation simulations to refine our understanding of $b_ϕ$ and $b_{ϕδ}$ towards improved constraints on $f_{\rm NL}$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Response approach to the integrated shear 3-point correlation function: the impact of baryonic effects on small scales

The integrated shear 3-point correlation function $ζ_{\pm}$ is a higher-order statistic of the cosmic shear field that describes the modulation of the 2-point correlation function $ξ_{\pm}$ by long-wavelength features in the field. Here, we introduce a new theoretical model to calculate $ζ_{\pm}$ that is accurate on small angular scales, and that allows to take baryonic feedback effects into account. Our model builds on the realization that the small-scale $ζ_{\pm}$ is dominated by the nonlinear matter bispectrum in the squeezed limit, which can be evaluated accurately using the nonlinear matter power spectrum and its first-order response functions to density and tidal field perturbations. We demonstrate the accuracy of our model by showing that it reproduces the small-scale $ζ_{\pm}$ measured in simulated cosmic shear maps. The impact of baryonic feedback enters effectively only through the corresponding impact on the nonlinear matter power spectrum, thereby permitting to account for these astrophysical effects on $ζ_{\pm}$ similarly to how they are currently accounted for on $ξ_{\pm}$. Using a simple idealized Fisher matrix forecast for a DES-like survey we find that, compared to $ξ_{\pm}$, a combined $ξ_{\pm}\ \&\ ζ_{\pm}$ analysis can lead to improvements of order $20-40\%$ on the constraints of cosmological parameters such as $σ_8$ or the dark energy equation of state parameter $w_0$. We find similar levels of improvement on the constraints of the baryonic feedback parameters, which strengthens the prospects for cosmic shear data to obtain tight constraints not only on cosmology but also on astrophysical feedback models. These are encouraging results that motivate future works on the integrated shear 3-point correlation function towards applications to real survey data.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The local PNG bias of neutral Hydrogen, ${\rm H_I}$

We use separate universe simulations with the IllustrisTNG galaxy formation model to predict the local PNG bias parameters $b_ϕ$ and $b_{ϕδ}$ of atomic neutral hydrogen, ${\rm H_I}$. These parameters and their relation to the linear density bias parameter $b_1$ play a key role in observational constraints of the local PNG parameter $f_{\rm NL}$ using the ${\rm H_I}$ power spectrum and bispectrum. Our results show that the popular calculation based on the universality of the halo mass function overpredicts the $b_ϕ(b_1)$ and $b_{ϕδ}(b_1)$ relations measured in the simulations. In particular, our results show that at $z \lesssim 1$ the ${\rm H_I}$ power spectrum is more sensitive to $f_{\rm NL}$ compared to previously thought ($b_ϕ$ is more negative), but is less sensitive at other epochs ($b_ϕ$ is less positive). We discuss how this can be explained by the competition of physical effects such as that large-scale gravitational potentials with local PNG (i) accelerate the conversion of hydrogen to heavy elements by star formation, (ii) enhance the effects of baryonic feedback that eject the gas to regions more exposed to ionizing radiation, and (iii) promote the formation of denser structures that shield the ${\rm H_I}$ more efficiently. Our numerical results can be used to revise existing forecast studies on $f_{\rm NL}$ using 21cm line-intensity mapping data. Despite this first step towards predictions for the local PNG bias parameters of ${\rm H_I}$, we emphasize that more work is needed to assess their sensitivity on the assumed galaxy formation physics and ${\rm H_I}$ modeling strategy.

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-CDM isocurvature galaxy bias with IllustrisTNG

We study the impact that baryon-CDM relative density perturbations $δ_{bc}$ have on galaxy formation using cosmological simulations with the IllustrisTNG model. These isocurvature (non-adiabatic) perturbations can be induced primordially, if multiple fields are present during inflation, and are generated before baryon-photon decoupling when baryons did not comove with CDM. The presence of long-wavelength $δ_{bc}$ perturbations in our simulations is mimicked by modifying the ratios of the cosmic densities of baryons $Ω_b$ and CDM $Ω_c$, at fixed total matter density $Ω_m$. We measure the corresponding galaxy bias parameter $b_δ^{bc}$ as the response of galaxy abundances to $δ_{bc}$. When selecting by total host halo mass, $b_δ^{bc}$ is negative and it decreases with mass and redshift. Stellar-mass selected simulated galaxies show a weaker or even the opposite trend because of the competing effects of $δ_{bc}$ on the halo mass function and stellar-to-halo-mass relations. We show that simple modeling of the latter two effects describes $b_δ^{bc}$ for stellar-mass-selected objects well. We find $b_δ^{bc} =0.6$ for $M_* = 10^{11}\ M_{\odot}/h$ and $z=0.5$, which is representative of BOSS DR12 galaxies. For $δ_{bc}$ modes generated by baryon-photon interactions, we estimate the impact on the DR12 power spectrum to be below $1\%$, and shifts on inferred distance and growth rate parameters should not exceed $0.1\%$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Compensated Isocurvature Perturbations in the Galaxy Power Spectrum

We investigate the potential of the galaxy power spectrum to constrain compensated isocurvature perturbations (CIPs), primordial fluctuations in the baryon density that are compensated by fluctuations in CDM density to ensure an unperturbed total matter density. We show that CIPs contribute to the galaxy overdensity at linear order, and if they are close to scale-invariant, their effects are nearly perfectly degenerate with the local PNG parameter $f_{\rm nl}$ if they correlate with the adiabatic perturbations. This degeneracy can however be broken by analyzing multiple galaxy samples with different bias parameters, or by taking CMB priors on $f_{\rm nl}$ into account. Parametrizing the amplitude of the CIP power spectrum as $P_{σσ} = A^2P_{\mathcal{R}\mathcal{R}}$ (where $P_{\mathcal{R}\mathcal{R}}$ is the adiabatic power spectrum) we find, for a number of fiducial galaxy samples in a simplified forecast setup, that constraints on $A$, relative to those on $f_{\rm nl}$, of order $σ_{A}/σ_{f_{\rm nl}} \approx 1-2$ are achievable for CIPs correlated with adiabatic perturbations, and $σ_{A}/σ_{f_{\rm nl}} \approx 5$ for the uncorrelated case. These values are independent of survey volume, and suggest that current galaxy data are already able to improve significantly on the tightest existing constraints on CIPs from the CMB. Future galaxy surveys that aim to achieve $σ_{f_{\rm nl}} \sim 1$ have the potential to place even stronger bounds on CIPs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the impact of galaxy bias uncertainties on primordial non-Gaussianity constraints

We study the impact that uncertainties on assumed relations between galaxy bias parameters have on constraints of the local PNG $f_{\rm NL}$ parameter. We focus on the relation between the linear density galaxy bias $b_1$ and local PNG bias $b_ϕ$ in an idealized forecast setup with multitracer galaxy power spectrum and bispectrum data. We consider two parametrizations of galaxy bias: 1) one inspired by the universality relation where $b_ϕ= 2δ_c\left(b_1 - p\right)$ and $p$ is a free parameter; and 2) another in which the product of bias parameters and $f_{\rm NL}$, like $f_{\rm NL} b_ϕ$, is directly fitted for. The constraints on the $f_{\rm NL}-p$ plane are markedly bimodal, and both the central value and width of marginalized constraints on $f_{\rm NL}$ depend sensitively on the priors on $p$. Assuming fixed $p=1$ in the constraints with a fiducial value of $p=0.55$ can bias the inferred $f_{\rm NL}$ by $0.5σ$ to $1σ$; priors $Δp \approx 0.5$ around this fiducial value are however sufficient in our setup to return unbiased constraints. In power spectrum analyses, parametrization 2, that makes no assumptions on $b_ϕ$, can distinguish $f_{\rm NL} \neq 0$ with the same significance as parametrization 1 assuming perfect knowledge of $b_ϕ$ (the value of $f_{\rm NL}$ is however left unknown). A drawback of parametrization 2 is that the addition of the bispectrum information is not as beneficial as in parametrization 1. Our results motivate strongly the incorporation of mitigation strategies for bias uncertainties in PNG constraint analyses, as well as further theoretical studies on the relations between bias parameters to better inform those strategies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Separate Universe Simulations with IllustrisTNG: baryonic effects on power spectrum responses and higher-order statistics

We measure power spectrum response functions in the presence of baryonic physical processes using separate universe simulations with the IllustrisTNG galaxy formation model. The response functions describe how the small-scale power spectrum reacts to long-wavelength perturbations and they can be efficiently measured with the separate universe technique by absorbing the effects of the long modes into a modified cosmology. Specifically, we focus on the total first-order matter power spectrum response to an isotropic density fluctuation $R_1(k,z)$, which is fully determined by the logarithmic derivative of the nonlinear matter power spectrum ${\rm dln}P_m(k,z)/{\rm dln}k$ and the growth-only response function $G_1(k,z)$. We find that $G_1(k,z)$ is not affected by the baryonic physical processes in the simulations at redshifts $z < 3$ and on all scales probed ($k \lesssim 15h/{\rm Mpc}$, i.e. length scales $\gtrsim 0.4 {\rm Mpc}/h$). In practice, this implies that the power spectrum fully specifies the baryonic dependence of its response function. Assuming an idealized lensing survey setup, we evaluate numerically the baryonic impact on the squeezed-lensing bispectrum and the lensing super-sample power spectrum covariance, which are given in terms of responses. Our results show that these higher-order lensing statistics can display varying levels of sensitivity to baryonic effects compared to the power spectrum, with the squeezed-bispectrum being the least sensitive. We also show that ignoring baryonic effects on lensing covariances slightly overestimates the error budget (and is therefore conservative from the point of view of parameter error bars) and likely has negligible impact on parameter biases in inference analyses.

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Validating estimates of the growth rate of structure with modified gravity simulations

We perform a validation of estimates of the growth rate of structure, described by the parameter combination $fσ_8$, in modified gravity cosmologies. We consider an analysis pipeline based on the redshift-space distortion modelling of the clustering wedges statistic of the galaxy correlation function and apply it to mock catalogues of $Λ{\rm CDM}$ and the normal branch of DGP cosmologies. We employ a halo occupation distribution approach to construct our mocks, which we ensure resemble the CMASS sample from BOSS in terms of the total galaxy number density and large scale amplitude of the power spectrum monopole. We show that the clustering wedges model successfully recovers the true growth rate difference between DGP and $Λ{\rm CDM}$, even for cases with over 40\% enhancement in $fσ_8$ compared to $Λ{\rm CDM}$. The unbiased performance of the clustering wedges model allows us to use the growth rate values estimated from the BOSS DR12 data to constrain the cross-over scale $r_c$ of DGP gravity to $\left[r_cH_0\right]^{-1} < 0.97$ ($2σ$) or $r_c > 3090\ {\rm Mpc}/h$, cutting into the interesting region of parameter space with $r_c \sim H_0^{-1}$ using constraints from the growth of structure alone.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

The velocity field of dark matter and galaxies reflects the continued action of gravity throughout cosmic history. We show that the low-order moments of the pairwise velocity distribution, $v_{12}$, are a powerful diagnostic of the laws of gravity on cosmological scales. In particular, the projected line-of-sight galaxy pairwise velocity dispersion, $σ_{12}(r)$, is very sensitive to the presence of modified gravity. Using a set of high-resolution N-body simulations we compute the pairwise velocity distribution and its projected line-of-sight dispersion for a class of modified gravity theories: the chameleon \fR gravity and Galileon gravity (cubic and quartic). The velocities of dark matter halos with a wide range of masses would exhibit deviations from General Relativity at the $(5-10)σ$ level. We examine strategies for detecting these deviations in galaxy redshift and peculiar velocity surveys. If detected, this signature would be a "smoking gun" for modified gravity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Halo model and halo properties in Galileon gravity cosmologies

We investigate the performance of semi-analytical modelling of large-scale structure in Galileon gravity cosmologies using results from N-body simulations. We focus on the Cubic and Quartic Galileon models that provide a reasonable fit to CMB, SNIa and BAO data. We demonstrate that the Sheth-Tormen mass function and linear halo bias can be calibrated to provide a very good fit to our simulation results. We also find that the halo concentration-mass relation is well fitted by a power law. The nonlinear matter power spectrum computed in the halo model approach is found to be inaccurate in the mildly nonlinear regime, but captures reasonably well the effects of the Vainshtein screening mechanism on small scales. In the Cubic model, the screening mechanism hides essentially all of the effects of the fifth force inside haloes. In the case of the Quartic model, the screening mechanism leaves behind residual modifications to gravity, which make the effective gravitational strength time-varying and smaller than the standard value. Compared to normal gravity, this causes a deficiency of massive haloes and leads to a weaker matter clustering on small scales. For both models, we show that there are realistic halo occupation distributions of Luminous Red Galaxies that can match both the observed large-scale clustering amplitude and the number density of these galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nonlinear structure formation in Nonlocal Gravity

We study the nonlinear growth of structure in nonlocal gravity models with the aid of N-body simulation and the spherical collapse and halo models. We focus on a model in which the inverse-squared of the d'Alembertian operator acts on the Ricci scalar in the action. For fixed cosmological parameters, this model differs from $Λ{\rm CDM}$ by having a lower late-time expansion rate and an enhanced and time-dependent gravitational strength ($\sim 6\%$ larger today). Compared to $Λ{\rm CDM}$ today, in the nonlocal model, massive haloes are slightly more abundant (by $\sim 10\%$ at $M \sim 10^{14} M_{\odot}/h$) and concentrated ($\approx 8\%$ enhancement over a range of mass scales), but their linear bias remains almost unchanged. We find that the Sheth-Tormen formalism describes the mass function and halo bias very well, with little need for recalibration of free parameters. The fitting of the halo concentrations is however essential to ensure the good performance of the halo model on small scales. For $k \gtrsim 1 h/{\rm Mpc}$, the amplitude of the nonlinear matter and velocity divergence power spectra exhibits a modest enhancement of $\sim 12\%$ to $15\%$, compared to $Λ{\rm CDM}$ today. This suggests that this model might only be distinguishable from $Λ{\rm CDM}$ by future observational missions. We point out that the absence of a screening mechanism may lead to tensions with Solar System tests due to local time variations of the gravitational strength, although this is subject to assumptions about the local time evolution of background averaged quantities.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nonlinear structure formation in the Cubic Galileon gravity model

We model the linear and nonlinear growth of large scale structure in the Cubic Galileon gravity model, by running a suite of N-body cosmological simulations using the {\tt ECOSMOG} code. Our simulations include the Vainshtein screening effect, which reconciles the Cubic Galileon model with local tests of gravity. In the linear regime, the amplitude of the matter power spectrum increases by $\sim 25%$ with respect to the standard $Λ$CDM model today. The modified expansion rate accounts for $\sim 20%$ of this enhancement, while the fifth force is responsible for only $\sim 5%$. This is because the effective unscreened gravitational strength deviates from standard gravity only at late times, even though it can be twice as large today. In the nonlinear regime ($k \gtrsim 0.1 h\rm{Mpc}^{-1}$), the fifth force leads to only a modest increase ($\lesssim 8%$) in the clustering power on all scales due to the very efficient operation of the Vainshtein mechanism. Such a strong effect is typically not seen in other models with the same screening mechanism. The screening also results in the fifth force increasing the number density of halos by less than 10%, on all mass scales. Our results show that the screening does not ruin the validity of linear theory on large scales which anticipates very strong constraints from galaxy clustering data. We also show that, whilst the model gives an excellent match to CMB data on small angular scales ($l \gtrsim 50$), the predicted integrated Sachs-Wolf effect is in tension with Planck/WMAP results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

We study spherical collapse in the Quartic and Quintic Covariant Galileon gravity models within the framework of the excursion set formalism. We derive the nonlinear spherically symmetric equations in the quasi-static and weak-field limits, focusing on model parameters that fit current CMB, SNIa and BAO data. We demonstrate that the equations of the Quintic model do not admit physical solutions of the fifth force in high density regions, which prevents the study of structure formation in this model. For the Quartic model, we show that the effective gravitational strength deviates from unity at late times ($z \lesssim 1$), becoming larger if the density is low, but smaller if the density is high. This shows that the Vainshtein mechanism at high densities is not enough to screen all of the modifications of gravity. This makes halos that collapse at $z \lesssim 1$ feel an overall weaker gravity, which suppresses halo formation. However, the matter density in the Quartic model is higher than in standard $Λ$CDM, which boosts structure formation and dominates over the effect of the weaker gravity. In the Quartic model there is a significant overabundance of high-mass halos relative to $Λ$CDM. Dark matter halos are also less biased than in $Λ$CDM, with the difference increasing appreciably with halo mass. However, our results suggest that the bias may not be small enough to fully reconcile the predicted matter power spectrum with LRG clustering data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The observational status of Galileon gravity after Planck

We use the latest CMB data from Planck, together with BAO measurements, to constrain the full parameter space of Galileon gravity. We constrain separately the three main branches of the theory known as the Cubic, Quartic and Quintic models, and find that all yield a very good fit to these data. Unlike in $Λ{\rm CDM}$, the Galileon model constraints are compatible with local determinations of the Hubble parameter and predict nonzero neutrino masses at over $5σ$ significance. We also identify that the low-$l$ part of the CMB lensing spectrum may be able to distinguish between $Λ{\rm CDM}$ and Galileon models. In the Cubic model, the lensing potential deepens at late times on sub-horizon scales, which is at odds with the current observational suggestion of a positive ISW effect. Compared to $Λ$CDM, the Quartic and Quintic models predict less ISW power in the low-$l$ region of the CMB temperature spectrum, and as such are slightly preferred by the Planck data. We illustrate that residual local modifications to gravity in the Quartic and Quintic models may render the Cubic model as the only branch of Galileon gravity that passes Solar System tests.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Linear perturbations in Galileon gravity models

We study the cosmology of Galileon modified gravity models in the linear perturbation regime. We derive the fully covariant and gauge invariant perturbed field equations using two different methods, which give consistent results, and solve them using a modified version of the {\tt CAMB} code. We find that, in addition to modifying the background expansion history and therefore shifting the positions of the acoustic peaks in the cosmic microwave background (CMB) power spectrum, the Galileon field can cluster strongly from early times, and causes the Weyl gravitational potential to grow, rather than decay, at late times. This leaves clear signatures in the low-$l$ CMB power spectrum through the modified integrated Sachs-Wolfe effect, strongly enhances the linear growth of matter density perturbations and makes distinctive predictions for other cosmological signals such as weak lensing and the power spectrum of density fluctuations. The quasi-static approximation is shown to work quite well from small to the near-horizon scales. We demonstrate that Galileon models display a rich phenomenology due to the large parameter space and the sensitive dependence of the model predictions on the Galileon parameters. Our results show that some Galileon models are already ruled out by present data and that future higher significance galaxy clustering, ISW and lensing measurements will place strong constraints on Galileon gravity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simulating the quartic Galileon gravity model on adaptively refined meshes

We develop a numerical algorithm to solve the high-order nonlinear derivative-coupling equation associated with the quartic Galileon model, and implement it in a modified version of the RAMSES N-body code to study the effect of the Galileon field on the large-scale matter clustering. The algorithm is tested for several matter field configurations with different symmetries, and works very well. This enables us to perform the first simulations for a quartic Galileon model which provides a good fit to the cosmic microwave background (CMB) anisotropy, supernovae and baryonic acoustic oscillations (BAO) data. Our result shows that the Vainshtein mechanism in this model is very efficient in suppressing the spatial variations of the scalar field. However, the time variation of the effective Newtonian constant caused by the curvature coupling of the Galileon field cannot be suppressed by the Vainshtein mechanism. This leads to a significant weakening of the strength of gravity in high-density regions at late times, and therefore a weaker matter clustering on small scales. We also find that without the Vainshtein mechanism the model would have behaved in a completely different way, which shows the crucial role played by nonlinearities in modified gravity theories and the importance of performing self-consistent N-body simulations for these theories.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The parameter space in Galileon gravity models

We present the first constraints on the full parameter space of the Galileon modified gravity model, considering both the cosmological parameters and the coefficients which specify the additional terms in the Lagrangian due to the Galileon field, which we call the Galileon parameters. We use the latest cosmic microwave background measurements, along with distance measurements from supernovae and baryonic acoustic oscillations, performing a Monte Carlo Markov Chain exploration of the 9-dimensional parameter space. The integrated Sachs-Wolfe signal can be very different in Galileon models compared to standard gravity, making it essential to use the full CMB data rather than the CMB distance priors. We demonstrate that meaningful constraints are only possible in the Galileon parameter space after taking advantage of a scaling degeneracy. We find that the Galileon model can fit the WMAP 9-year results better than the standard Λ-Cold Dark Matter model, but gives a slightly worse fit overall once lower redshift distance measurements are included. The best-fitting cosmological parameters (e.g. matter density, scalar spectral index, fluctuation amplitude) can differ by more than 2σ in the Galileon model compared with ΛCDM. We highlight other potential constraints of the Galileon model using galaxy clustering and weak lensing measurements.

0 citations0 reviewsNo code linkedOpen access