Search

We consider the recently proposed possibility that dark energy (DE) and baryons may scatter through a pure momentum exchange process, leaving the background evolution unaffected. Earlier work has shown that, even for barn-scale cross-sections, the imprints of this scattering process on linear cosmological observables is too tiny to be observed. We therefore turn our attention to non-linear scales, and for the first time investigate the signatures of DE-baryon scattering on the non-linear formation of cosmic structures, by running a suite of large N-body simulations. The observables we extract include the non-linear matter power spectrum, halo mass function, and density and baryon fraction profiles of halos. We find that in the non-linear regime the signatures of DE-baryon scattering are significantly larger than their linear counterparts, due to the important role of angular momentum in collapsing structures, and potentially observable. The most promising observables in this sense are the baryon density and baryon fraction profiles of halos, which can potentially be constrained by a combination of kinetic Sunyaev-Zeldovich (SZ), thermal SZ, and weak lensing measurements. Overall, our results indicate that future prospects for cosmological and astrophysical direct detection of non-gravitational signatures of dark energy are extremely bright., 21 pages, 19 sub-figures arranged in 10 figures. Minor modifications, added references and further discussions on screening. Version accepted for publication in MNRAS

The Lambda-Cold Dark Matter (LCDM) model agrees with most of the cosmological observations, but has some hindrances from observed data at smaller scales such as galaxies. Recently, Berezhiani and Khoury (2015) proposed a new theory involving interacting superfluid dark matter with three model parameters, which explains galactic dynamics with great accuracy. In the present work, we study the cosmological behaviour of this model in the linear regime of cosmological perturbations. In particular, we compute both analytically and numerically the matter linear growth factor and obtain new bounds for the model parameters which are significantly stronger than previously found. These new constraints come from the fact that structures within the superfluid dark matter framework grow quicker than in LCDM, and quite rapidly when the DM-baryon interactions are strong., Comment: 16 pages, 6 figures

E. Abdalla, G. F. Abellán, A. Aboubrahim et al., The standard Λ Cold Dark Matter (ΛCDM) cosmological model provides a good description of a wide range of astrophysical and cosmological data. However, there are a few big open questions that make the standard model look like an approximation to a more realistic scenario yet to be found. In this paper, we list a few important goals that need to be addressed in the next decade, taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant H0, the σ8–S8 tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the 5.0σ tension between the Planck CMB estimate of the Hubble constant H0 and the SH0ES collaboration measurements. After showing the H0 evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the Planck CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density Ωm, and the amplitude or rate of the growth of structure (σ8, fσ8). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the H0–S8 tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals. Finally, we give an overview of upgraded experiments and next-generation space missions and facilities on Earth that will be of crucial importance to address all these open questions., Amin Aboubrahim is supported by the BMBF under contract 05P21PMCAA and by the DFG through the Research Training Network 2149 “Strong and Weak Interactions - from Hadrons to Dark Matter”. Adriano Agnello is supported by a Villum Experiment Grant, project number 36225. Özgür Akarsu acknowledges the support by the Turkish Academy of Sciences in the scheme of the Outstanding Young Scientist Award (TÜBA-GEBİP). Yashar Akrami is supported by Richard S. Morrison Fellowship and LabEx ENS-ICFP: ANR-10-LABX-0010/ANR-10-IDEX-0001-02 PSL*. George Alestas is supported by the project “Dioni: Computing Infrastructure for Big-Data Processing and Analysis” (MIS No. 5047222) co-funded by European Union (ERDF) and Greece through Operational Program “Competitiveness, Entrepreneurship and Innovation”, NSRF 2014-2020. Luca Amendola acknowledges support from DFG project 456622116 and from the CAPES-DAAD bilateral project “Data Analysis and Model Testing in the Era of Precision Cosmology”. Luis A. Anchordoqui and Jorge F. Soriano are supported by the U.S. National Science Foundation (NSF Grant PHY-2112527). Mario Ballardini acknowledges financial support from the contract ASI/INAF for the Euclid mission n.2018-23-HH.0. Micol Benetti acknowledges the Istituto Nazionale di Fisica Nucleare (INFN), sezione di Napoli, iniziativa specifica QGSKY. David Benisty acknowledges the support the supports of the Blavatnik and the Rothschild fellowships. John Blakeslee is supported by NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the US National Science Foundation. Thomas Buchert has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement ERC advanced grant 740021–ARTHUS, PI: Thomas Buchert). Erminia Calabrese acknowledges support from the STFC Ernest Rutherford Fellowship ST/M004856/2, STFC Consolidated Grant ST/S00033X/1 and from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (Grant agreement No. 849169). Salvatore Capozziello acknowledges the Istituto Nazionale di Fisica Nucleare (INFN), sezione di Napoli, iniziative specifiche QGSKY and MOONLIGHT2. Javier de Cruz Pérez is supported by a FPI fellowship associated to the project FPA2016-76005-C2-1-P. Peter Denton acknowledges support from the US Department of Energy under Grant Contract DE-SC0012704. Eleonora Di Valentino is supported by a Royal Society Dorothy Hodgkin Research Fellowship. Keith R. Dienes was supported in part by the U.S. Department of Energy under Grant DE-FG02-13ER41976 / DE-SC0009913, and also by the U.S. National Science Foundation through its employee IR/D program. Celia Escamilla-Rivera is supported by DGAPA-PAPIIT UNAM Project TA100122 and acknowledges the Royal Astronomical Society as FRAS 10147 and the Cosmostatistics National Group (CosmoNag) project. Noemi Frusciante is supported by Fundação para a Ciência e a Tecnologia (FCT) through the research grants UIDB/04434/2020, UIDP/04434/2020, PTDC/FIS-OUT/29048/2017, CERN/FIS-PAR/0037/2019 and the personal FCT grant “CosmoTests – Cosmological tests of gravity theories beyond General Relativity” with ref. number CEECIND/00017/2018 and the FCT project “BEYLA–BEYond LAmbda” with ref. number PTDC/FIS-AST/0054/2021. Adrià Gómez-Valent is funded by the Instituto Nazionale di Fisica Nucleare (INFN) through the project “Dark Energy and Modified Gravity Models in the light of Low-Redshift Observations” (n. 22425/2020). Asta Heinesen has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement ERC advanced grant 740021–ARTHUS, PI: Thomas Buchert). J. Colin Hill acknowledges support from NSF grant AST-2108536. The Flatiron Institute is supported by the Simons Foundation. Mustapha Ishak acknowledges that this material is based upon work supported in part by the Department of Energy, Office of Science, under Award Number DE-SC0022184. Michael Klasen is supported by the BMBF under contract 05P21PMCAA and by the DFG through the Research Training Network 2149 “Strong and Weak Interactions - from Hadrons to Dark Matter”. Suresh Kumar gratefully acknowledges support from the Science and Engineering Research Board (SERB), Govt. of India (File No. CRG/2021/004658). Ruth Lazkoz is supported by the Spanish Ministry of Science and Innovation through research projects FIS2017-85076-P (comprising FEDER funds), and also by the Basque Government and Generalitat Valenciana through research projects GIC17/116-IT956-16 and PROMETEO/2020/079 respectively. Benjamin L'Huillier would like to acknowledge the support of the National Research Foundation of Korea (NRF-2019R1I1A1A01063740) and the support of the Korea Institute for Advanced Study (KIAS) grant funded by the government of Korea. Jackson Levi Said would like to acknowledge support from Cosmology@MALTA which is supported by the University of Malta. Roy Maartens is supported by the South African Radio Astronomy Observatory and the National Research Foundation (Grant No. 75415). Valerio Marra thanks CNPq (Brazil) and FAPES (Brazil) for partial financial support. The work of Yuto Minami was supported in part by the Japan Society for the Promotion of Science (JSPS) KAKENHI, Grants No. JP20K14497. The work of Carlos Martins was financed by FEDER—Fundo Europeu de Desenvolvimento Regional funds through the COMPETE 2020—Operational Programme for Competitiveness and Internationalisation (POCI), and by Portuguese funds through FCT - Fundação para a Ciência e a Tecnologia in the framework of the project POCI-01-0145-FEDER-028987 and PTDC/FIS-AST/28987/2017. Olga Mena is supported by the Spanish grants PID2020-113644GB-I00, PROMETEO/2019/083 and by the European ITN project HIDDeN (H2020-MSCA-ITN-2019//860881-HIDDeN). Cristian Moreno-Pulido is funded by PID2019-105614GB-C21 and FPA2016-76005-C2-1-P (MINECO, Spain), 2017-SGR-929 (Generalitat de Catalunya) and CEX2019-000918-M (ICCUB) and partially supported by the fellowship 2019 FI_B 00351. Michele Moresco acknowledges support from MIUR, PRIN 2017 (grant 20179ZF5KS) and grants ASI n.I/023/12/0 and ASI n.2018-23-HH.0. Suvodip Mukherjee is supported by the Simons Foundation. Research at Perimeter Institute is supported in part by the Government of Canada through the Department of Innovation, Science and Economic Development and by the Province of Ontario through the Ministry of Colleges and Universities. Pran Nath is supported in part by the NSF Grant PHY-1913328. Savvas Nesseris acknowledges support from the Research Project No. PGC2018-094773-B-C32 and the Centro de Excelencia Severo Ochoa Program No. CEX2020-001007-S. Rafael Nunes acknowledges support from the Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP, São Paulo Research Foundation) under the project No. 2018/18036-5. Eoin Ó Colgáin was supported by the National Research Foundation of Korea grant funded by the Korea government (MSIT) (NRF-2020R1A2C1102899). Supriya Pan acknowledges the financial supports from the Science and Engineering Research Board, Govt. of India, under Mathematical Research Impact-Centric Support Scheme (File No. MTR/2018/000940) and The Department of Science and Technology (DST), Govt. of India, under the Scheme “Fund for Improvement of S&T Infrastructure (FIST)” [File No. SR/FST/MS-I/2019/41]. Santiago E. Perez Bergliaffa acknowledges partial support from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)- Código de Financiamento 001, and Universidade do Estado do Rio de Janeiro (Brazil). Leandros Perivolaropoulos acknowledges support by the Hellenic Foundation for Research and Innovation (H.F.R.I.), under the “First call for H.F.R.I. Research Projects to support Faculty members and Researchers and the procurement of high-cost research equipment Grant” (Project Number: 789). Fabrizio Renzi is supported by the NWO and the Dutch Ministry of Education, Culture and Science (OCW), and from the D-ITP consortium, a program of the NWO that is funded by the OCW. Nils Schöneberg acknowledges the support of the following Maria de Maetzu fellowship grant: Esta publicación es parte de la ayuda CEX2019-000918-M, financiado por MCIN/AEI/10.13039/501100011033. Anjan A Sen acknowledges the funding from SERB, Govt of India under the research grants no: CRG/2020/004347 and MTR/20l9/000599. Arman Shafieloo would like to acknowledge the support by National Research Foundation of Korea NRF-2021M3F7A1082053, and the support of the Korea Institute for Advanced Study (KIAS) grant funded by the government of Korea. M.M. Sheikh-Jabbari acknowledges the support by SarAmadan grant No. ISEF/M/400122. Joan Solà Peracaula is funded by PID2019-105614GB-C21 and FPA2016-76005-C2-1-P (MINECO, Spain), 2017-SGR-929 (Generalitat de Catalunya), CEX2019-000918-M (ICCUB) and also partially supported by the COST Association Action CA18108 “Quantum Gravity Phenomenology in the Multimessenger Approach (QG-MM)”. Denitsa Staicova is supported by Bulgarian NSF grant KP-06-N 38/11. Glenn Starkman is partly supported by a Department of Energy grant DESC0009946 to the particle astrophysics theory group at CWRU. Brooks Thomas is supported in part by the National Science Foundation under Grant PHY-2014104. Luca Visinelli has received support from the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement “TALeNT” No. 754496 (H2020-MSCA-COFUND-2016 FELLINI). Shao-Jiang Wang is supported by the National Key Research and Development Program of China Grant No. 2021YFC2203004, No. 2021YFA0718304, the National Natural Science Foundation of China Grant No. 12105344, and the China Manned Space Project with NO.CMS-CSST-2021-B01. Weiqiang Yang has been supported by the National Natural Science Foundation of China under Grants No. 12175096 and No. 11705079, and Liaoning Revitalization Talents Program under Grant no. XLYC1907098. Gong-Bo Zhao is supported by the National Key Basic Research and Development Program of China (No. 2018YFA0404503), NSFC Grants 11925303, 11720101004, and a grant of CAS Interdisciplinary Innovation Team.

The $\Lambda$CDM model provides a good fit to a large span of cosmological data but harbors areas of phenomenology. With the improvement of the number and the accuracy of observations, discrepancies among key cosmological parameters of the model have emerged. The most statistically significant tension is the $4-6\sigma$ disagreement between predictions of the Hubble constant $H_0$ by early time probes with $\Lambda$CDM model, and a number of late time, model-independent determinations of $H_0$ from local measurements of distances and redshifts. The high precision and consistency of the data at both ends present strong challenges to the possible solution space and demand a hypothesis with enough rigor to explain multiple observations--whether these invoke new physics, unexpected large-scale structures or multiple, unrelated errors. We present a thorough review of the problem, including a discussion of recent Hubble constant estimates and a summary of the proposed theoretical solutions. Some of the models presented are formally successful, improving the fit to the data in light of their additional degrees of freedom, restoring agreement within $1-2\sigma$ between {\it Planck} 2018, using CMB power spectra data, BAO, Pantheon SN data, and R20, the latest SH0ES Team measurement of the Hubble constant ($H_0 = 73.2 \pm 1.3{\rm\,km\,s^{-1}\,Mpc^{-1}}$ at 68\% confidence level). Reduced tension might not simply come from a change in $H_0$ but also from an increase in its uncertainty due to degeneracy with additional physics, pointing to the need for additional probes. While no specific proposal makes a strong case for being highly likely or far better than all others, solutions involving early or dynamical dark energy, neutrino interactions, interacting cosmologies, primordial magnetic fields, and modified gravity provide the best options until a better alternative comes along.[Abridged], Comment: Invited Review of about 134 pages with many figures and tables. Accepted for publication in CQG. A sample code for producing similar whisker plots (e.g. Figs. 1, 2) associated with this Review for any choice of the data is made publicly available online at: https://github.com/lucavisinelli/H0TensionRealm

Since the very beginning of astronomy the location of objects on the sky has been a fundamental observational quantity that has been taken for granted. While precise two dimensional positional information is easy to obtain for observations in the electromagnetic spectrum, the positional accuracy of current and near future gravitational wave detectors is limited to between tens and hundreds of square degrees, which makes it extremely challenging to identify the host galaxies of gravitational wave events or to confidently detect any electromagnetic counterparts. Gravitational wave observations provide information on source properties and distances that is complementary to the information in any associated electromagnetic emission and that is very hard to obtain in any other way. Observing systems with multiple messengers thus has scientific potential much greater than the sum of its parts. A gravitational wave detector with higher angular resolution would significantly increase the prospects for finding the hosts of gravitational wave sources and triggering a multi-messenger follow-up campaign. An observatory with arcminute precision or better could be realised within the Voyage 2050 programme by creating a large baseline interferometer array in space and would have transformative scientific potential. Precise positional information of standard sirens would enable precision measurements of cosmological parameters and offer new insights on structure formation; a high angular resolution gravitational wave observatory would allow the detection of a stochastic background and resolution of the anisotropies within it; it would also allow the study of accretion processes around black holes; and it would have tremendous potential for tests of modified gravity and the discovery of physics beyond the Standard Model., 26 pages, 2 figures. White paper submitted to ESA's Voyage 2050 call on behalf of the LISA Consortium 2050 Task Force

We investigate cosmological signatures of the minimal theory of massive gravity (MTMG). To this aim, we simulate the normal branch of the MTMG by employing the \textsc{Ramses} \mbox{$N$-body} code and extending it with an effective gravitational constant $G_{\rm eff}$. We implement an environment-dependent $G_{\rm eff}$ as a function of the graviton mass and the local energy density as predicted by MTMG. We find that halo density profiles are not a good probe for MTMG, because deviations from general relativity (GR) are quite small. Similarly, the matter power spectra show deviations only at the percentage level. However, we find a clear difference between MTMG and GR in that voids are denser in MTMG than in GR. As measuring void profiles is quite a complex task from an observational point of view, a better probe of MTMG would be the halo abundances. In this case, MTMG creates a larger amount of massive halos, while there is a suppression in the abundance of small halos., 8 pages, 7 figures. Updated to reflect peer-reviewed version for publication

We consider the possibility that dark energy and baryons might scatter off each other. The type of interaction we consider leads to a pure momentum exchange, and does not affect the background evolution of the expansion history. We parametrize this interaction in an effective way at the level of Boltzmann equations. We compute the effect of dark energy-baryon scattering on cosmological observables, focusing on the Cosmic Microwave Background (CMB) temperature anisotropy power spectrum and the matter power spectrum. Surprisingly, we find that even huge dark energy-baryon cross-sections $\sigma_{xb} \sim {\cal O}({\rm b})$, which are generically excluded by non-cosmological probes such as collider searches or precision gravity tests, only leave an insignificant imprint on the observables considered. In the case of the CMB temperature power spectrum, the only imprint consists in a sub-percent enhancement or depletion of power (depending whether or not the dark energy equation of state lies above or below $-1$) at very low multipoles, which is thus swamped by cosmic variance. These effects are explained in terms of differences in how gravitational potentials decay in the presence of a dark energy-baryon scattering, which ultimately lead to an increase or decrease in the late-time integrated Sachs-Wolfe power. Even smaller related effects are imprinted on the matter power spectrum. The imprints on the CMB are not expected to be degenerate with the effects due to altering the dark energy sound speed. We conclude that, while strongly appealing, the prospects for a direct detection of dark energy through cosmology do not seem feasible when considering realistic dark energy-baryon cross-sections. As a caveat, our results hold to linear order in perturbation theory., Comment: 15 pages, 7 figures. Title changed, minor modifications added comments on non-linearities. Version accepted for publication in MNRAS

We investigate a well known scenario of interaction in the dark sector where the vacuum energy is interacting with cold dark matter throughout the cosmic evolution in light of the cosmic microwave background (CMB) data from final Planck 2018 release. In addition to this minimal scenario, we generalize the model baseline by including the properties of neutrinos, such as the neutrino mass scale ($M_{\nu}$) and the effective number of neutrino species ($N_{\rm eff}$) as free parameters, in order to verify the possible effects that such parameters might generate on the coupling parameter, and vice versa. As already known, we again confirm that in light of the Planck 2018 data, such dark coupling can successfully solve the $H_0$ tension (with and without the presence of neutrinos). Concerning the properties of neutrinos, we find that $M_{\nu}$ may be wider than expected within the $\Lambda$CDM model and $N_{\rm eff}$ is fully compatible with three neutrino species (similar to $\Lambda$CDM prevision). The parameters characterizing the properties of neutrinos do not correlate with the coupling parameter of the interaction model. When considering the joint analysis of CMB from Planck 2018 and an estimate of $H_0$ from Hubble Space Telescope 2019 data, {\it we find an evidence for a non-null value of the coupling parameter at more than 3$\sigma$ confidence-level.} We also discuss the possible effects on the interacting scenario due to the inclusion of baryon acoustic oscillations data with Planck 2018. Our main results updating the dark sectors' interaction and neutrino properties in the model baseline, represent a new perspective in this direction. Clearly, a possible new physics in light of some dark interaction between dark energy and dark matter can serve as an alternative to $\Lambda$CDM scenario to explain the observable Universe, mainly in light of the current tension on $H_0$., Comment: 13 pages, 6 tables (3 tables in the main text and 3 tables in the appendix), 3 captioned figures; accepted for publication in JCAP

The aim of the present work is to better understand the gravitational drag forces, i.e. dynamical friction, acting on massive objects moving through a self-interacting Bose-Einstein condensate, also known as a superfluid, at finite temperatures. This is relevant for light scalar models of dark matter with weak self-interactions that require nonzero temperatures, or that have been heated inside galaxies. We derived expressions for dynamical friction using linear perturbation theory, and compared these to numerical simulations in which nonlinear effects are included. After testing the linear result, it was applied to the Fornax dwarf spheroidal galaxy, and two of its gravitationally bound globular clusters. Dwarf spheroidals are well-suited for indirectly probing properties of dark matter, and so by estimating the rate at which these globular clusters are expected to sink into their host halo due to dynamical friction, we inferred limits on the superfluid dark matter parameter space. The dynamical friction in a finite-temperature superfluid is found to behave very similarly to the zero-temperature limit, even when the thermal contributions are large. However, when a critical velocity for the superfluid flow is included, the friction force can transition from the zero-temperature value to the value in a conventional fluid. Increasing the mass of the perturbing object induces a similar transition to when lowering the critical velocity. When applied to two of Fornax's globular clusters, we find that the parameter space preferred in the literature for a zero-temperature superfluid yields decay times that are in agreement with observations. However, the present work suggests that increasing the temperature, which is expected to change the preferred parameter space, may lead to very small decay times, and therefore pose a problem for finite-temperature superfluid models of dark matter., Comment: 16 pages, 8 figures, 3 tables. Accepted for publication in A&A

It is always interesting to investigate how well can a future experiment perform with respect to others (present or future ones). Cosmology is really an exciting field where a lot of puzzles are still unknown. In this article we consider a generalized dark energy (DE) scenario where anisotropic stress is present. We constrain this generalized cosmic scenario with an aim to investigate how gravitational waves standard sirens (GWSS) may constrain the anisotropic stress, which according to the standard cosmological probes, remains unconstrained. In order to do this, we generate the luminosity distance measurements from $\mathcal{O} (10^3)$ mock GW events which match the expected sensitivity of the Einstein Telescope. Our analyses report that, first of all, GWSS can give better constraints on various cosmological parameters compared to the usual cosmological probes, but the viscous sound speed appearing due to the dark energy anisotropic stress, is totally unconstrained even after the inclusion of GWSS., Comment: 17 pages, 5 tables and 12 figures; version accepted for publication by MNRAS

We consider two stable Interacting Dark Matter -- Dark Energy models and confront them against current Cosmic Microwave Background data from the \textit{Planck} satellite. We then generate luminosity distance measurements from ${\cal O}(10^3)$ mock Gravitational Wave events matching the expected sensitivity of the proposed Einstein Telescope. We use these to forecast how the addition of Gravitational Wave standard sirens data can improve current limits on the Dark Matter -- Dark Energy coupling strength ($\xi$). We find that the addition of Gravitational Waves data can reduce the current uncertainty by a factor of $5$. Moreover, if the underlying cosmological model truly features Dark Matter -- Dark Energy interactions with a value of $\xi$ within the currently allowed $1\sigma$ upper limit, the addition of Gravitational Wave data would help disentangle such an interaction from the standard case of no interaction at a significance of more than $3\sigma$., Comment: 16 pages, 3 tables, 4 figures; version published in JCAP

Searches for modified gravity in the large-scale structure try to detect the enhanced amplitude of density fluctuations caused by the fifth force present in many of these theories. Neutrinos, on the other hand, suppress structure growth below their free-streaming length. Both effects take place on comparable scales, and uncertainty in the neutrino mass leads to a degeneracy with modified gravity parameters for probes that are measuring the amplitude of the matter power spectrum. We explore the possibility to break the degeneracy between modified gravity and neutrino effects in the growth of structures by considering kinematic information related to either the growth rate on large scales or the virial velocities inside of collapsed structures. In order to study the degeneracy up to fully non-linear scales, we employ a suite of $N$-body simulations including both $f(R)$ modified gravity and massive neutrinos. Our results indicate that velocity information provides an excellent tool to distinguish massive neutrinos from modified gravity. Models with different values of neutrino masses and modified gravity parameters possessing a comparable matter power spectrum at a given time have different growth rates. This leaves imprints in the velocity divergence, which is therefore better suited than the amplitude of density fluctuations to tell the models apart. In such models with a power spectrum comparable to $\Lambda$CDM today, the growth rate is strictly enhanced. We also find the velocity dispersion of virialised clusters to be well suited to constrain deviations from general relativity without being affected by the uncertainty in the sum of neutrino masses., Comment: 8 pages, 4 figures, submitted to A&A

We intend to understand cosmological structure formation within the framework of superfluid models of dark matter with finite temperatures. Of particular interest is the evolution of small-scale structures where the pressure and superfluid properties of the dark matter fluid are prominent. We compare the growth of structures in these models with the standard cold dark matter paradigm and non-superfluid dark matter. The equations for superfluid hydrodynamics were computed numerically in an expanding $\Lambda$CDM background with spherical symmetry; the effect of various superfluid fractions, temperatures, interactions, and masses on the collapse of structures was taken into consideration. We derived the linear perturbation of the superfluid equations, giving further insights into the dynamics of the superfluid collapse. We found that while a conventional dark matter fluid with self-interactions and finite temperatures experiences a suppression in the growth of structures on smaller scales, as expected due to the presence of pressure terms, a superfluid can collapse much more efficiently than was naively expected due to its ability to suppress the growth of entropy perturbations and thus gradients in the thermal pressure. We also found that the cores of the dark matter halos initially become more superfluid during the collapse, but eventually reach a point where the superfluid fraction falls sharply. The formation of superfluid dark matter halos surrounded by a normal fluid dark matter background is therefore disfavored by the present work., Comment: 10 pages, 9 figures. Accepted by A&A

Modified gravity (MG) theories aim to reproduce the observed acceleration of the Universe by reducing the dark sector while simultaneously recovering General Relativity (GR) within dense environments. Void studies appear to be a suitable scenario to search for imprints of alternative gravity models on cosmological scales. Voids cover an interesting range of density scales where screening mechanisms fade out, which reaches from a density contrast $\delta \approx -1$ close to their centers to $\delta \approx 0$ close to their boundaries. We present an analysis of the level of distinction between GR and two modified gravity theories, the Hu-Sawicki $f(R)$ and the symmetron theory. This study relies on the abundance, linear bias, and density profile of voids detected in n-body cosmological simulations. We define voids as connected regions made up of the union of spheres with a {\it \textup{mean}} density given by $\overline\rho_v=0.2\,\overline\rho_m$, but disconnected from any other voids. We find that the height of void walls is considerably affected by the gravitational theory, such that it increases for stronger gravity modifications. Finally, we show that at the level of dark matter n-body simulations, our constraints allow us to distinguish between GR and MG models with $|f_{R0}| > 10^{-6}$ and $z_{SSB} > 1$. Differences of best-fit values for MG parameters that are derived independently from multiple void probes may indicate an incorrect MG model. This serves as an important consistency check., Comment: 15 pages, 12 figures

The detection of gravitational waves (GW) by the LIGO and Virgo collaborations offers a whole new range of possible tests and opens up a new window which may shed light on the nature of dark energy and dark matter. In the present work we investigate how future gravitational waves data could help to constrain different dynamical dark energy models. In particular, we perform cosmological forecastings of a class of well known and most used dynamical dark energy models using the third-generation gravitational wave detector, the Einstein Telescope. We have considered 1000 simulated GW events in order to constrain the parameter space of the dynamical dark energy models. Our analyses show that the inclusion of the GW data from the Einstein Telescope, significantly improves the parameter space of the dynamical dark energy models compared to their constraints extracted from the standard cosmological probes, namely, the cosmic microwave observations, baryon acoustic oscillations distance measurements, Supernove type Ia, and the Hubble parameter measurements., Comment: 27 pages, 10 tables and 13 figures, version accepted for publication in Physical Review D

We present a formalism to study screening mechanisms in modified theories of gravity via perturbative methods in different cosmological scenarios. We consider Einstein frame posed theories that are recast as Jordan frame theories, where a known formalism is employed, though the resulting non-linearities of the Klein-Gordon equation acquire an explicit coupling between matter and the scalar field, which is not present in Jordan frame theories. The obtained growth functions are then separated in screening and non-screened contributions to facilitate its analysis. This allows us to compare several theoretical models and to recognize patterns which can be used to differentiate models and their screening mechanisms. In particular, we find anti-screening features in the Symmetron model. In opposition, chameleon type theories, both in the Jordan and in the Einstein frame, always present a screening behaviour. Up to third order in perturbation, we find no anti-screening behaviour in theories with a Vainshtein mechanism, such as the DGP and the cubic Galileon., Comment: 12 pages, 9 figures. Replaced to match version accepted for publication in A&A

The evolution of the virial overdensity $\Delta_{\rm vir}$ for $\Lambda$CDM and seven dynamical dark-energy models is investigated in the extended spherical collapse model (SCM). Here the virialization process is naturally achieved by introducing shear and rotation instead of using the virial theorem. We generalise two approaches proposed in the literature and show that, regardless of the dark-energy model, the new virialization term can be calibrated on the peculiar velocity of the shell as measured from Einstein-de Sitter simulations. The two virialization recipes qualitatively reproduce the features of the ordinary SCM, i.e., a constant $\Delta_{\rm vir}$ for the EdS model and time-variation for dark-energy models, but without any mass dependence. Depending on the actual description of virialization and on the dark-energy model, the value of $\Delta_{\rm vir}$ varies between 10 and 40 percent. We use the new recipes to predict the surface-mass-density profile of dark matter haloes and the number of convergence density peaks for LSST- and Euclid-like weak lensing surveys., Comment: 27 pages, 9 figures, 3 tables. It matches the published version on JCAP. New sections added on stable clustering and on the convergence peaks. The numerical code is available upon request to the corresponding author

While the origin and composition of dark matter and dark energy remains unknown, it is possible that they might represent two manifestations of a single entity, as occurring in unified dark sector models. On the other hand, advances in our understanding of the dark sector of the Universe might arise from Cosmic Dawn, the epoch when the first stars formed. In particular, the first detection of the global 21-cm absorption signal at Cosmic Dawn from the EDGES experiment opens up a new arena wherein to test models of dark matter and dark energy. Here, we consider generalized and modified Chaplygin gas models as candidate unified dark sector models. We first constrain these models against Cosmic Microwave Background data from the \textit{Planck} satellite, before exploring how the inclusion of the global 21-cm signal measured by EDGES can improve limits on the model parameters, finding that the uncertainties on the parameters of the Chaplygin gas models can be reduced by a factor between $1.5$ and $10$. We also find that within the generalized Chaplygin gas model, the tension between the CMB and local determinations of the Hubble constant $H_0$ is reduced from $\approx 4\sigma$ to $\approx 1.3\sigma$. In conclusion, we find that the global 21-cm signal at Cosmic Dawn can provide an extraordinary window onto the physics of unified dark sectors., Comment: 18 pages, 3 Tables, 2 figures; version accepted for publication in JCAP

We present a new code, SCALAR, based on the high-resolution hydrodynamics and N-body code RAMSES, to solve the Schr\"odinger equation on adaptive refined meshes. The code is intended to be used to simulate axion or fuzzy dark matter models where the evolution of the dark matter component is determined by a coupled Schr\"odinger-Poisson equation, but it can also be used as a standalone solver for both linear and non-linear Schr\"odinger equations with any given external potential. This paper describes the numerical implementation of our solver and presents tests to demonstrate how accurately it operates., Comment: 17 pages, 11 figures

We investigate the observational consequences of a novel class of stable interacting dark energy (IDE) models, featuring interactions between dark matter (DM) and dark energy (DE). In the first part of our work, we start by considering two IDE models which are known to present early-time linear perturbation instabilities. Applying a transformation depending on the dark energy equation of state (EoS) to the DM-DE coupling, we then obtain two novel stable IDE models. Subsequently, we derive robust and accurate constraints on the parameters of these models, assuming a constant EoS $w_x$ for the DE fluid, in light of some of the most recent publicly available cosmological data. These include Cosmic Microwave Background (CMB) temperature and polarization anisotropy measurements from the \textit{Planck} satellite, a selection of Baryon Acoustic Oscillation measurements, Supernovae Type-Ia luminosity distance measurements from the JLA sample, and measurements of the Hubble parameter up to redshift $2$ from cosmic chronometers. Our analysis displays a mild preference for the DE fluid residing in the phantom region ($w_x, Comment: 21 pages, 6 Figures, 6 Tables; Accepted in JCAP

We study a novel interacting dark energy $-$ dark matter scenario where the anisotropic stress of the large scale inhomogeneities is considered. The dark energy has a constant equation of state and the interaction model produces stable perturbations. The resulting picture is constrained using different astronomical data aiming to measure the impact of the anisotropic stress on the cosmological parameters. Our analyses show that a non-zero interaction in the dark sector is allowed while a non-interaction scenario is recovered within 68\% CL. The anisotropic stress is also constrained to be small, and its zero value is permitted within 68\% CL. The dark energy equation of state, $w_x$, is also found to be close to `$-1$' boundary. However, from the ratio of the CMB TT spectra, we see that the model has a mild deviation from the $\Lambda$CDM cosmology while such deviation is almost forbidden from the CMB TT spectra alone. Although the deviation is not much significant, but from the present data, we cannot exclude such deviation. Overall, at the background level, the model is close to the $\Lambda$CDM cosmology while at the level of perturbations, a non-zero but a very small interaction in the dark sector is permitted. Perhaps, a more accurate conclusion can be made with the next generation of surveys. We also found that the region $w_x < -1$, is found to be effective to release the tension on $H_0$. Finally, from the Bayesian analysis, we find that $\Lambda$CDM remains in still preferred over the interacting scenarios., Comment: 18 pages, 9 figures, 7 tables; Published version in MNRAS

We investigate the effects of a disformal coupling between dark energy and dark matter in the predictions of the spherical collapse and its signatures in galaxy cluster number counts. We find that the disformal coupling has no significant effects on spherical collapse at high redshifts, and in particular during matter domination epoch. However, at lower redshifts, the extrapolated linear density contrast at collapse close to redshift $z \lesssim 1$ and overdensity at virialization can be strongly suppressed by a disformal coupling between dark energy and dark matter. We also find that disformal coupling can have different imprints on cluster number counts compared with conformal coupling, such that the disformal coupling can strongly suppress the predicted number of clusters per redshift interval at $z > 0.1$ while enhance the number of cluster at $z < 0.05$. Using the specifications of eROSITA survey, we find that the disformal coupling between dark energy and dark matter can be tightly constrained by cluster number counts., 22 pages, 8 figures, 1 table, accepted for publication in Physical Review D

Multiple modifications of general relativity (GR) have been proposed in the literature in order to understand the nature of the accelerated expansion of the Universe. However, thus far all the predictions of GR have been confirmed with constantly increasing accuracy. In this work, we study the imprints of a particular class of models -- "screened" modified gravity theories -- on cosmic filaments. We have utilized the $N$-body code ISIS/RAMSES to simulate the symmetron model and the Hu-Sawicky $f(R)$ model, and we post-process the output with DisPerSE to identify the filaments of the cosmic web. We investigated how the global properties of the filaments -- such as their lengths, masses, and thicknesses -- as well as their radial density and speed profiles change under different gravity theories. We find that filaments are, on average, shorter and denser in modified gravity models compared to in $\Lambda$CDM. We also find that the speed profiles of the filaments are enhanced, consistent with theoretical expectations. Overall, our results suggest that cosmic filaments can be an effective complementary probe of screened modified gravity theories on Mpc scales., Comment: 11 pages, 17 figures, accepted by A&A

We reconsider the dynamics of the Universe in the presence of interactions in the cosmological dark sector. A class of interacting models is introduced via a real function $f\left(r\right)$ of the ratio $r$ between the energy densities of the (pressureless) cold dark matter (CDM) and dark energy (DE). The subclass of models for which the ratio $r$ depends only on the scale factor is shown to be equivalent to unified models of the dark sector, i.e. models for which the CDM and DE components can be combined in order to form a unified dark fluid. For specific choices of the function $f\left(r\right)$ we recover several models already studied in the literature. We analyse various special cases of this type of interacting models using a suitably modified version of the CLASS code combined with MontePython in order to constrain the parameter space with the data from supernova of type SNe Ia (JLA), the Hubble constant $H_{0}$, cosmic chronometers (CC), baryon acoustic oscilations (BAO) and data from the Planck satellite (Planck TT). Our analysis shows that even if data from the late Universe ($H_{0}$, SNe Ia and CC) indicate an interaction in the dark sector, the data related to the early Universe (BAO and Planck TT) constrain this interaction substantially, in particular for cases in which the background dynamics is strongly affected., Comment: 43 pages, 22 figures

We investigate the efficiency of screening mechanisms in the hybrid metric-Palatini gravity. The value of the field is computed around spherical bodies embedded in a background of constant density. We find a thin shell condition for the field depending on the background field value. In order to quantify how the thin shell effect is relevant, we analyze how it behaves in the neighborhood of different astrophysical objects (planets, moons or stars). We find that the condition is very well satisfied except only for some peculiar objects. Furthermore we establish bounds on the model using data from solar system experiments such as the spectral deviation measured by the Cassini mission and the stability of the Earth-Moon system, which gives the best constraint to date on $f(R)$ theories. These bounds contribute to fix the range of viable hybrid gravity models., 7 pages, 2 figures. Accepted for publication in Phys. Rev. D

Consistency relations of large-scale structures provide exact nonperturbative results for cross-correlations of cosmic fields in the squeezed limit. They only depend on the equivalence principle and the assumption of Gaussian initial conditions, and remain nonzero at equal times for cross-correlations of density fields with velocity or momentum fields, or with the time derivative of density fields. We show how to apply these relations to observational probes that involve the integrated Sachs-Wolfe effect or the kinematic Sunyaev-Zeldovich effect. In the squeezed limit, this allows us to express the three-point cross-correlations, or bispectra, of two galaxy or matter density fields, or weak lensing convergence fields, with the secondary Cosmic Microwave Background (CMB) distortion in terms of products of a linear and a nonlinear power spectrum. In particular, we find that cross-correlations with the integrated Sachs-Wolfe effect show a specific angular dependence. These results could be used to test the equivalence principle and the primordial Gaussianity, or to check the modeling of large-scale structures., 10 pages

Extending general relativity by adding extra degrees of freedom is a popular approach for explaining the accelerated expansion of the Universe and to build high energy completions of the theory of gravity. The presence of such new degrees of freedom is, however, tightly constrained from several observations and experiments that aim to test general relativity in a wide range of scales. The viability of a given modified theory of gravity, therefore, strongly depends on the existence of a screening mechanism that suppresses the extra degrees of freedom. We perform simulations, and find that waves propagating in the new degrees of freedom can significantly impact the efficiency of some screening mechanisms, thereby threatening the viability of these modified gravity theories. Specifically, we show that the waves produced in the symmetron model can increase the amplitude of the fifth force and the parametrized post-Newtonian parameters by several orders of magnitude., Comment: 7 pages, 2 figures. Supplemental material included as an appendix. Updated to reflect version published in Phys. Rev. Lett., with several minor adjustments and addition of an appendix

In this work we will test an alternative model of gravity belonging to the large family of galileon models. It is characterized by an intrinsic breaking of the Vainshtein mechanism inside large astrophysical objects, thus having possibly detectable observational signatures. We will compare theoretical predictions from this model with the observed total mass profile for a sample of clusters of galaxies. The profiles are derived using two complementary tools: X-ray hot intra-cluster gas dynamics, and strong and weak gravitational lensing. We find that a dependence with the dynamical internal status of each cluster is possible; for those clusters which are very close to be relaxed, and thus less perturbed by possible astrophysical local processes, the galileon model gives a quite good fit to both X-ray and lensing observations. Both masses and concentrations for the dark matter halos are consistent with earlier results found in numerical simulations and in the literature, and no compelling statistical evidence for a deviation from general relativity is detectable from the present observational state. Actually, the characteristic galileon parameter $\Upsilon$ is always consistent with zero, and only an upper limit ($\lesssim0.086$ at $1\sigma$, $\lesssim0.16$ at $2\sigma$, and $\lesssim0.23$ at $3\sigma$) can be established. Some interesting distinctive deviations might be operative, but the statistical validity of the results is far from strong, and better data would be needed in order to either confirm or reject a potential tension with general relativity., Comment: 26 pages, 3 tables, 4 figures. Accepted for publication on Phys. Rev. D

We investigate how three different possibilities of neutrino mass hierarchies, namely normal, inverted, and degenerate, can affect the observational constraints on three well known dynamical dark energy models, namely the Chevallier-Polarski-Linder, logarithmic, and the Jassal-Bagla-Padmanabhan parametrizations. In order to impose the observational constraints on the models, we performed a robust analysis using Planck 2015 temperature and polarization data, Supernovae type Ia from Joint Light curve analysis, baryon acoustic oscillations distance measurements, redshift space distortion characterized by $f(z)\sigma_8(z)$ data, weak gravitational lensing data from Canada-France-Hawaii Telescope Lensing Survey, and cosmic chronometers data plus the local value of the Hubble parameter. We find that different neutrino mass hierarchies return similar fit on almost all model parameters and mildly change the dynamical dark energy properties., Comment: 10 pages, 5 captioned figures, 4 tables; Published version in Phys. Rev. D

Cosmological local observables are at best statistically determined by the fundamental theory describing inflation. When the scalar inflaton is coupled uniformly to a collection of subdominant massless gauge vectors, rotational invariance is obeyed locally. However, the statistical isotropy of fluctuations is spontaneously broken by gauge modes whose wavelength exceeds our causal horizon. This leads to a landscape picture where primordial correlators depend on the position of the observer. We compute the stochastic corrections to the curvature power spectrum, show the existence of a new local observable (the shape of the quadrupole), and constrain the theory using Planck limits. © 2014 The Authors., SCOPUS: ar.j, info:eu-repo/semantics/published

There are currently no reliable methods to measure transverse velocities of galaxies. This is an important piece of information that could allow us to probe the physics of structure formation as well as testing the underlying theory of gravity. The slingshot effect, a special case of the Integrated Sachs-Wolfe effect, is expected to create dipole signals in the temperature fluctuations of the Cosmic Microwave Background Radiation (CMB). This effect creates a hot spot behind and a cold spot in front of moving massive objects. The dipole signal created by the slingshot effect can be used to measure transverse velocities, but because the signal is expected to be weak, the effect has not been measured yet. The aim is to show that the slingshot effect can be measured by stacking the signals of galaxies falling into a collapsing cluster. We evaluate if the effect can probe modified gravity. We use data from a simulated galaxy catalogue (MDPL2) to mimic observations. We identify a massive galaxy cluster, and make maps of the slingshot effect around infalling galaxies. We add uncorrelated Gaussian noise to each map. The maps are rotated according to the direction to the cluster centre, such that the dipole signal will add up constructively when stacking. We compare each stack to a dipole stencil and we find the probability for a false positive in the absence of the slingshot signal. Each galaxy gives a signal of around $\Delta T/T \approx 10^{-9}$, while the precision of CMB experiments of today are $\Delta T/T \approx 4 \times 10^{-6}$. By stacking around 10 000 galaxies, the slingshot signal can be over the detectable threshold with experiments of today. However, future CMB experiments must be used to be certain of the strength of the observed signal., Comment: 14 pages total (3 pages appendix), 4 figures

We consider a cosmological model with an interaction between dark matter and dark energy which leaves the background cosmology unaffected and only affects the evolution of the perturbations. This is achieved by introducing a coupling given in terms of the relative velocities of dark matter and dark energy. This interaction has the distinctive feature of appearing predominantly on small scales, where peculiar velocities can become important. We confront the predictions of the model to cosmological observations and find a potential alleviation of the known tension in the amplitude of density perturbations as measured by low redshift galaxy surveys and the Planck data. The model also predicts a shift in the turnover of the matter power spectrum which does not depend on the horizon at equality (fixed by the background cosmology and, thus, unaffected by the perturbations) and is entirely due to the interaction between dark matter and dark energy. A bias in the peculiar velocity between baryons and dark matter is also shown to be a unique feature of this type of interactions in the dark sector., Comment: 26 pages, 5 figures. Small changes and added references to match published version

In this paper, we study the chameleon profile in inhomogeneous density distributions and find that the fifth force in thin shell near the surface is weaker from what expected in homogeneous density distributions. Also, we check the validity of quasi-static approximation for the chameleon scalar field in the astrophysical time scales. We have investigated the rolling down behaviour of the scalar field on its effective potential inside a one solar mass red giant star by using MESA code. We have found that the scalar field is fast enough to follow the minimum of the potential. This adiabatic behaviour reduces the fifth force and extends the screened regions to lower densities where the field has smaller mass and was expected to be unscreened. As a consequence, the star evolution is similar to what expected from standard general relativity. In addition, considering the stability of star, an approximate constraint on the coupling constant $\beta$ is found., Comment: 27 pages, 15 figures

We use cosmological hydrodynamical simulations to study the effect of screened modified gravity models on the mass estimates of galaxy clusters. In particular, we focus on two novel aspects: (i) we study modified gravity models in which baryons and dark matter are coupled with different strengths to the scalar field, and, (ii) we put the simulation results into the greater context of a general screened-modified gravity parametrization. We have compared the mass of clusters inferred via lensing versus the mass inferred via kinematical measurements as a probe of violations of the equivalence principle at Mpc scales. We find that estimates of cluster masses via X-ray observations is mainly sensitive to the coupling between the scalar degree of freedom and baryons - while the kinematical mass is mainly sensitive to the coupling to dark matter. Therefore, the relation between the two mass estimates is a probe of a possible non-universal coupling between the scalar field, the standard model fields, and dark matter. Finally, we used observational data of kinetic, thermal and lensing masses to place constraints on deviations from general relativity on cluster scales for a general parametrization of screened modified gravity theories which contains $f(R)$ and Symmetron models. We find that while the kinematic mass can be used to place competitive constraints, using thermal measurements is challenging as a potential non-thermal contribution is degenerate with the imprint of modified gravity., 10 pages, 8 figures; matches published version in A&A

We use a spherical model and an extended excursion set formalism with drifting diffusive barriers to predict the abundance of cosmic voids in the context of general relativity as well as f(R) and symmetron models of modified gravity. We detect spherical voids from a suite of N-body simulations of these gravity theories and compare the measured void abundance to theory predictions. We find that our model correctly describes the abundance of both dark matter and galaxy voids, providing a better fit than previous proposals in the literature based on static barriers. We use the simulation abundance results to fit for the abundance model free parameters as a function of modified gravity parameters, and show that counts of dark matter voids can provide interesting constraints on modified gravity. For galaxy voids, more closely related to optical observations, we find that constraining modified gravity from void abundance alone may be significantly more challenging. In the context of current and upcoming galaxy surveys, the combination of void and halo statistics including their abundances, profiles and correlations should be effective in distinguishing modified gravity models that display different screening mechanisms., 20 pages, 12 figures, submitted to PRD

We present exact kinematic consistency relations for cosmological structures that do not vanish at equal times and can thus be measured in surveys. These rely on cross-correlations between the density and velocity, or momentum, fields. Indeed, the uniform transport of small-scale structures by long wavelength modes, which cannot be detected at equal times by looking at density correlations only, gives rise to a shift in the amplitude of the velocity field that could be measured. These consistency relations only rely on the weak equivalence principle and Gaussian initial conditions. They remain valid in the non-linear regime and for biased galaxy fields. They can be used to constrain non-standard cosmological scenarios or the large-scale galaxy bias., 5 pages

We propose a new method to test modified gravity theories, taking advantage of the available data on extrasolar planets. We computed the deviations from the Kepler third law and use that to constrain gravity theories beyond General Relativity. We investigate gravity models which incorporate three screening mechanisms: the Chameleon, the Symmetron and the Vainshtein. We find that data from exoplanets orbits are very sensitive to the screening mechanisms putting strong constraints in the parameter space for the Chameleon models and the Symmetron, complementary and competitive to other methods, like interferometers and solar system. With the constraints on Vainshtein we are able to work beyond the hypothesis that the crossover scale is of the same order of magnitude than the Hubble radius $r_c \sim H_0^{-1}$, which makes the screening work automatically, testing how strong this hypothesis is and the viability of other scales., 7 pages, 5 figures, Published in Physics Letters B

We investigate the propagation of scalar waves induced by matter sources in the context of scalar-tensor theories of gravity which include screening mechanisms for the scalar degree of freedom. The usual approach when studying these theories in the non-linear regime of cosmological perturbations is based on the assumption that scalar waves travel at the speed of light. Within General Relativity such approximation is good and leads to no loss of accuracy in the estimation of observables. We find, however, that mass terms and non-linearities in the equations of motion lead to propagation and dispersion velocities significantly different from the speed of light. As the group velocity is the one associated to the propagation of signals, a reduction of its value has direct impact on the behavior and dynamics of nonlinear structures within modified gravity theories with screening. For instance, the internal dynamics of galaxies and satellites submerged in large dark matter halos could be affected by the fact that the group velocity is smaller than the speed of light. It is therefore important, within such framework, to take into account the fact that different part of a galaxy will see changes in the environment at different times. A full non-static analysis may be necessary under those conditions.

Modified gravity theories with a screening mechanism have acquired much interest recently in the quest for a viable alternative to General Relativity on cosmological scales, given their intrinsic property of being able to pass Solar System scale tests and, at the same time, to possibly drive universe acceleration on much larger scales. Here, we explore the possibility that the same screening mechanism, or its breaking at a certain astrophysical scale, might be responsible of those gravitational effects which, in the context of general relativity, are generally attributed to Dark Matter. We consider a recently proposed extension of covariant Galileon models in the so-called "beyond Horndeski" scenario, where a breaking of the Vainshtein mechanism is possible and, thus, some peculiar observational signatures should be detectable and make it distinguishable from general relativity. We apply this model to a sample of clusters of galaxies observed under the \textit{CLASH} survey, using both new data from gravitational lensing events and archival data from X-ray intra-cluster hot gas observations. In particular, we use the latter to model the gas density, and then use it as the only ingredient in the matter clusters' budget to calculate the expected lensing convergence map. Results show that, in the context of this extended Galileon, the assumption of having only gas and no Dark Matter at all in the clusters is able to match observations. We also obtain narrow and very interesting bounds on the parameters which characterize this model. In particular, we find that, at least for one of them, the general relativity limit is excluded at $2\sigma$ confidence level, thus making this model clearly statistically different and competitive with respect to general relativity., Comment: 16 pages, 6 figures, 1 Table, matching version published in JCAP

We propose that the mass-temperature relation of galaxy clusters is a prime candidate for testing gravity theories beyond Einstein's general relativity, for modified gravity models with universal coupling between matter and the scalar field. For non-universally coupled models we discover that the impact of modified gravity can remain hidden from the mass-temperature relation. Using cosmological simulations, we find that in modified gravity the mass-temperature relation varies significantly from the standard gravity prediction of $M \propto T^{1.73}$. To be specific, for symmetron models with a coupling factor of $\beta=1$ we find a lower limit to the power law as $M\propto T^{1.6}$; and for f(R) gravity we compute predictions based on the model parameters. We show that the mass-temperature relation, for screened modified gravities, is significantly different from that of standard gravity for the less massive and colder galaxy clusters, while being indistinguishable from Einstein's gravity at massive, hot galaxy clusters. We further investigate the mass-temperature relation for other mass estimates than the thermal mass estimate, and discover that the gas mass-temperature results show an even more significant deviations from Einstein's gravity than the thermal mass-temperature., Comment: Updated, submitted to A&A

Alternative theories of gravity typically invoke an environment-dependent screening mechanism to allow phenomenologically interesting deviations from general relativity (GR) to manifest on larger scales, while reducing to GR on small scales. The observation of the transition from screened to unscreened behavior would be compelling evidence for beyond-GR physics. We show that pairwise peculiar velocity statistics, in particular the relative radial velocity dispersion, $\sigma_\parallel$ , can be used to observe this transition when they are binned by some measure of halo environment. We established this by measuring the radial velocity dispersion between pairs of halos in N-body simulations for three $f(R)$ gravity and four symmetron models. We developed an estimator involving only line-of-sight velocities to show that this quantity is observable, and binned the results in halo mass, ambient density, and the isolatedness of halos. Ambient density is found to be the most relevant measure of environment; it is distinct from isolatedness, and correlates well with theoretical expectations for the symmetron model. By binning $\sigma_\parallel$ in ambient density, we find a strong environment-dependent signature for the symmetron models, with the velocities showing a clear transition from GR to non-GR behavior. No such transition is observed for $f(R)$, as the relevant scales are deep in the unscreened regime. Observations of the relative radial velocity dispersion in forthcoming peculiar velocity surveys, if binned appropriately by environment, therefore offer a valuable way of detecting the screening signature of modified gravity., Comment: 11 pages, 8 figures. Updated to journal version; minor changes

In order to determine the observable signatures of modified gravity theories, it is important to consider the effect of baryonic physics. We use a modified version of the ISIS code to run cosmological hydrodynamic simulations to study degeneracies between modified gravity and radiative hydrodynamical processes. Of these, one was the standard $\Lambda$CDM model and four were variations of the Symmetron model. For each model we ran three variations of baryonic processes: non-radiative hydrodynamics; cooling and star formation; and cooling, star formation, and supernova feedback. We construct stacked gas density, temperature, and dark matter density profiles of the halos in the simulations, and study the differences between them. We find that both radiative variations of the models show degeneracies between their processes and at least two of the three parameters defining the Symmetron model., Comment: 9 pages, 4 figures, matches version accepted to A&A

We investigate statistical properties of the distribution of matter at redshift zero in disformal gravity by using N-body simulations. The disformal model studied here consists of a conformally coupled symmetron field with an additional exponential disformal term. We conduct cosmological simulations to discover the impact of the new disformal terms in the matter power spectrum, halo mass function, and radial profile of the scalar field. We calculated the disformal geodesic equation and the equation of motion for the scalar field. We then implemented these equations into the N-body code ISIS, which is a modified gravity version of the code RAMSES. The presence of a conformal symmetron field increases both the power spectrum and mass function compared to standard gravity on small scales. Our main finding is that the newly added disformal terms tend to counteract these effects and can make the evolution slightly closer to standard gravity. We finally show that the disformal terms give rise to oscillations of the scalar field in the centre of the dark matter haloes., Updated version to reflect the journal accepted paper. Added one figure. 7 pages, 7 figures

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 $\Lambda$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., Comment: 31 pages, 16 figures. Version to appear in MNRAS

In this paper we present the first results from implementing two scalar-tensor modified gravity theories, the symmetron and the Hu-Sawicki $f(R)$-gravity model, into a hydrodynamic N-body code with dark matter particles and a baryonic ideal gas. The study is a continuation of previous work where the symmetron and $f(R)$ have been successfully implemented in the RAMSES code, but for dark matter only. By running simulations, we show that the deviation from $\Lambda$CDM in these models for the gas density profiles are significantly lower than the dark matter equivalents. When it comes to the matter power-spectrum we find that hydrodynamic simulations agree very well with dark matter only simulations as long as we consider scales larger than $k\sim 0.5$ h/Mpc. In general the effects of modified gravity on the baryonic gas is found to not always mirror the effects it has on the dark matter. The largest signature is found when considering temperature profiles. We find that the gas temperatures in the modified gravity model studied here show deviations, when compared to $\Lambda$CDM, that can be a factor of a few larger than the deviations found in density profiles and power spectra., Comment: 11 pages, 10 figures, submitted to MNRAS

In this work we study the dynamics of the Local Group (LG) within the context of cosmological models beyond General Relativity (GR). Using observable kinematic quantities to identify candidate pairs we build up samples of simulated LG-like objects drawing from $f(R)$, symmetron, DGP and quintessence N-body simulations together with their $\Lambda$CDM counterparts featuring the same initial random phase realisations. The variables and intervals used to define LG-like objects are referred to as Local Group model; different models are used throughout this work and adapted to study their dynamical and kinematic properties. The aim is to determine how well the observed LG-dynamics can be reproduced within cosmological theories beyond GR, We compute kinematic properties of samples drawn from alternative theories and $\Lambda$CDM and compare them to actual observations of the LG mass, velocity and position. As a consequence of the additional pull, pairwise tangential and radial velocities are enhanced in modified gravity and coupled dark energy with respect to $\Lambda$CDM, inducing significant changes to the total angular momentum and energy of the LG. For example, in models such as $f(R)$ and the symmetron this increase can be as large as $60\%$, peaking well outside of the $95\%$ confidence region allowed by the data. This shows how simple considerations about the LG dynamics can lead to clear small-scale observational signatures for alternative scenarios, without the need of expensive high-resolution simulations., Comment: 14 pages, 3 figures. Accepted for publication in MNRAS

Aims. We investigated the gravitational effects of a scalar field within scalar-tensor gravity as an alternative to dark matter. Motivated by chameleon, symmetron and f(R)-gravity models, we studied a phenomenological scenario where the scalar field has both a mass (i.e. interaction length) and a coupling constant to the ordinary matter which scale with the local properties of the considered astrophysical system. Methods. We analysed the feasibility of this scenario using the modified gravitational potential obtained in its context and applied it to the galactic and hot gas/stellar dynamics in galaxy clusters and elliptical/spiral galaxies respectively. This is intended to be a first step in assessing the viability of this new approach in the context of "alternative gravity" models. Results. The main results are: 1. the velocity dispersion of elliptical galaxies can be fitted remarkably well by the suggested scalar field, with model significance similar to a classical Navarro-Frenk-White dark halo profile; 2. the analysis of the stellar dynamics and the gas equilibrium in elliptical galaxies has shown that the scalar field can couple with ordinary matter with different strengths (different coupling constants) producing and/or depending on the different clustering state of matter components; 3. elliptical and spiral galaxies, combined with clusters of galaxies, show evident correlations among theory parameters which suggest the general validity of our results at all scales and a way toward a possible unification of the theory for all types of gravitational systems we considered. All these results demonstrate that the proposed scalar field scenario can work fairly well as an alternative to dark matter., Comment: 23 pages, 15 figures, 5 tables, accepted for publication on Astronomy & Astrophysics

Screened modified gravity predicts potentially large signatures in the peculiar velocity field that makes it an interesting probe to test gravity on cosmological scales. We investigate the signatures induced by the Symmetron and a Chameleon $f(R)$ model in the peculiar velocity field using $N$-body simulations. By studying fifth force and halo velocity profiles we identify three general categories of effects found in screened modified gravity models: a fully screened regime where we recover $\Lambda$CDM to high precision, an unscreened regime where the fifth force is in full operation, and, a partially screened regime where screening occurs in the inner part of a halo, but the fifth force is active at larger radii. These three regimes can be pointed out very clearly by analyzing the deviation in the maximum cluster velocity. Observationally, the partially screened regime is of particular interest since an uniform increase of the gravitational force - as present in the unscreened regime - is degenerate with the (dynamical) halo mass estimate, and, thus, hard to detect., Comment: 9 pages, 7 figures

In the context of scalar tensor theories for gravity, there is a universally adopted hypothesis when running N-body simulations that time derivatives in the equation of motion for the scalar field are negligible. In this work we propose to test this assumption for one specific scalar-tensor model with a gravity screening mechanism: the symmetron. To this end, we implemented the necessary modifications to include the non-static terms in the N-body code Ramses. We present test cases and results from cosmological simulations. Our main finding when comparing static vs. non-static simulations is that the global power spectrum is only slightly modified when taking into account the inclusion of non-static terms. On the contrary, we find that the calculation of the local power spectrum gives different measurements. Such results imply one must be careful when assuming the quasi-static approximation when investigating the environmental effects of modified gravity and screening mechanisms in structure formation of halos and voids distributions., 12 pages, 8 figures, matches version accepted for publication in PRD

Aims. We investigate the gravitational redshift in clusters of galaxies in the symmetron and Hu-Sawicky $f(R)$ models. The characteristic feature of both models is the screening mechanism that hides the fifth force in dense environments recovering general relativity. Methods. We use $N$-body simulations that were run using the code \texttt{Isis}, which includes scalar fields, to analyse the deviation of observables in modified gravity models with respect to $\Lambda$CDM. Results. We find that due to the presence of the screening, the deviation is highly dependent on the halo mass. For instance, the $f(R)$ parameters $|f_{R0}|=10^{-5}$, $n=1$ cause an enhancement of the gravitational signal by up to 50% for halos with masses between $10^{13}\,M_\odot h^{-1}$ and $10^{14}\,M_\odot h^{-1}$. The characteristic mass range where the fifth force is most active varies with the model parameters. The usual assumption is that the presence of a fifth force leads to a deeper potential well and thus, a stronger gravitational redshift. However, we find that in cases in which only the central regions of the halos are screened, there could also be a weaker gravitational redshift., Comment: 10 pages, 5 figures