Search

Recent analyses show that $\Lambda$CDM-based models optimised to reproduce the clustering of massive galaxies overestimate their gravitational lensing by about 30\%, the so-called lensing is low problem. Using a state-of-the-art hydrodynamical simulation, we show that this discrepancy reflects shortcomings in standard galaxy-halo connection models rather than tensions within the $\Lambda$CDM paradigm itself. Specifically, this problem results from ignoring a variety of galaxy formation effects, including assembly bias, segregation of satellite galaxies relative to dark matter, and baryonic effects on the matter distribution. All these effects contribute towards overestimating gravitational lensing and, when combined, explain the amplitude and scale dependence of the lensing is low problem. We conclude that simplistic galaxy-halo connection models are inadequate to interpret clustering and lensing simultaneously, and that it is crucial to employ more sophisticated models for the upcoming generation of large-scale surveys., Comment: 16 pages, 12 figures, accepted by MNRAS

The spatial distribution of galaxies and their gravitational lensing signal offer complementary tests of galaxy formation physics and cosmology. However, their synergy can only be fully exploited if both probes are modelled accurately and consistently. In this paper, we demonstrate that this can be achieved using an extension of Subhalo Abundance Matching, dubbed SHAMe. Specifically, we use mock catalogues built from the TNG300 hydrodynamical simulation to show that SHAMe can simultaneously model the multipoles of the redshift-space galaxy correlation function and galaxy-galaxy lensing, without noticeable bias within the statistical sampling uncertainties of a SDSS volume and on scales r = [0.6-30] Mpc/h. Modelling the baryonic processes in galaxy-galaxy lensing with a baryonification scheme allows SHAMe's range of validity to be extended to r = [0.1-30] Mpc/h. Remarkably, our model achieves this level of precision with just five free parameters beyond those describing the baryonification model. At fixed cosmology, we find that galaxy-galaxy lensing provides a general consistency test but little additional information on galaxy modelling parameters beyond that encoded in the redshift-space multipoles. It does, however, improve constraints if only the projected correlation function is available, as in surveys with only photometric redshifts. We expect SHAMe to have a higher fidelity across a wider range of scales than more traditional methods such as Halo Occupation Distribution modelling. Thus it should provide a significantly more powerful and more robust tool for analysing next-generation large-scale surveys., 14 pages, 7 figures. Accepted by MNRAS

We study the relations among the parameters of the hybrid Lagrangian bias expansion model, fitting biased auto and cross power spectra up to $k_{\rm max} = 0.7 \, h \, \mathrm{Mpc}^{-1}$. We consider $\sim 8000$ halo and galaxy samples, with different halo masses, redshifts, galaxy number densities, and varying the parameters of the galaxy formation model. Galaxy samples are obtained through state-of-the-art extended subhalo abundance matching techniques and include both stellar-mass and star-formation-rate selected galaxies. All of these synthetic galaxies samples are publicly available at https://bacco.dipc.org/galpk.html. We find that the hybrid Lagrangian bias model provides accurate fits to all of our halo and galaxy samples. The coevolution relations between galaxy bias parameters, although roughly compatible with those obtained for haloes, show systematic shifts and larger scatter. We explore possible sources of this difference in terms of dependence on halo occupation and assembly bias of each sample. The bias parameter relations displayed in this work can be used as a prior for future Bayesian analyses employing the hybrid Lagrangian bias expansion model., 14 pages, 12 figures, galaxy power spectra available at https://bacco.dipc.org/galpk.html

Prompt cusps are the densest quasi-equilibrium dark matter objects; one forms at the instant of collapse within every isolated peak of the initial cosmological density field. They have power-law density profiles, $\rho \propto r^{-1.5}$ with central phase-space density set by the primordial velocity dispersion of the dark matter. At late times they account for $\sim 1\%$ of the dark matter mass but for $>90\%$ of its annihilation luminosity in all but the densest regions, where they are tidally disrupted. Here we demonstrate that individual stellar encounters, rather than the mean galactic tide, are the dominant disruptors of prompt cusps within galaxies. Their cumulative effect is fully (though stochastically) characterised by an impulsive shock strength $B_* = 2\pi G\int\rho_*({\bf x}(t))\, \mathrm{d}t$ where $\rho_*$, the total mass density in stars, is integrated over a cusp's entire post-formation trajectory. Stellar encounters and mean tides have only a small effect on the halo annihilation luminosity seen by distant observers, but this is not true for the Galactic halo because of the Sun's position. For a 100 GeV WIMP, Earth-mass prompt cusps are predicted, and stellar encounters suppress their mean annihilation luminosity by a factor of two already at 20 kpc, so that their annihilation emission is predicted to appear almost uniform over the sky. The Galactic Center $\gamma$-ray Excess is thus unaffected by cusps. If it is indeed dark matter annihilation radiation, then prompt cusps in the outer Galactic halo and beyond must account for 20-80% of the observed isotropic $\gamma$-ray background in the 1 to 10 GeV range., Comment: 22 pages, 23 figures, for a video and code see https://github.com/jstuecker/cusp-encounters

The abundance of dark matter haloes is one of the key probes of the growth of structure and expansion history of the Universe. Theoretical predictions for this quantity usually assume that, when expressed in a certain form, it depends only on the mass variance of the linear density field. However, cosmological simulations have revealed that this assumption breaks, leading to 10–20 per cent systematic effects. In this paper, we employ a specially designed suite of simulations to further investigate this problem. Specifically, we carry out cosmological N-body simulations where we systematically vary growth history at a fixed linear density field, or vary the power spectrum shape at a fixed growth history. We show that the halo mass function generically depends on these quantities, thus showing a clear signal of non-universality. Most of this effect can be traced back to the way in which the same linear fluctuation grows differently into the non-linear regime depending on details of its assembly history. With these results, we propose a parameterization with explicit dependence on the linear growth rate and power spectrum shape. Using an independent suite of simulations, we show that this fitting function accurately captures the mass function of haloes over cosmologies spanning a vast parameter space, including massive neutrinos and dynamical dark energy. Finally, we employ this tool to improve the accuracy of so-called cosmology-rescaling methods and show they can deliver 2 per cent accurate predictions for the halo mass function over the whole range of currently viable cosmologies.

In this work, we carry out a suite of specially-designed numerical simulations to shed further light on dark matter (DM) subhalo survival at mass scales relevant for gamma-ray DM searches, a topic subject to intense debate nowadays. Specifically, we have developed and employed an improved version of DASH, a GPU $N$-body code, to study the evolution of low-mass subhaloes inside a Milky Way-like halo with unprecedented accuracy, reaching solar-mass and sub-parsec resolution in our simulations. We simulate subhaloes with varying mass, concentration, and orbital properties, and consider the effect of the gravitational potential of the Milky Way galaxy itself. More specifically, we analyze the evolution of both the bound mass fraction and annihilation luminosity of subhaloes, finding that most subhaloes survive until present time, even though in some cases they lose more than 99% of their mass at accretion. Baryons in the host induce a much more severe mass loss, especially when the subhalo orbit is more parallel to the galactic disk. Many of these subhaloes cross the solar galactocentric radius, thus making it easier to detect their annihilation fluxes from Earth. We find subhaloes orbiting a DM-only halo with a pericentre in the solar vicinity to lose 70-90% of their initial annihilation luminosity at redshift zero, which increases up to 99% when baryons are also included in the host. We find a strong relation between subhalo's mass loss and the effective tidal field at pericentre. Indeed, much of the dependence on concentration, orbital parameters, host potential and baryonic components can be explained through this single parameter. In addition to shedding light on the survival of low-mass galactic subhaloes, our results can provide detailed predictions that will aid current and future quests for the nature of DM., Comment: 19 pages, 16 figures. Submitted to MNRAS, comments welcome!

One of the main predictions of excursion set theory is that the clustering of dark matter haloes only depends on halo mass. However, it has been long established that the clustering of haloes also depends on other properties, including formation time, concentration, and spin; this effect is commonly known as halo assembly bias. We use a suite of gravity-only simulations to study the dependence of halo assembly bias on cosmology; these simulations cover cosmological parameters spanning 10$\sigma$ around state-of-the-art best-fitting values, including standard extensions of the $\Lambda$CDM paradigm such as neutrino mass and dynamical dark energy. We find that the strength of halo assembly bias presents variations smaller than 0.05 dex across all cosmologies studied for concentration and spin selected haloes, letting us conclude that the dependence of halo assembly bias upon cosmology is negligible. We then study the dependence of galaxy assembly bias (i.e. the manifestation of halo assembly bias in galaxy clustering) on cosmology using subhalo abundance matching. We find that galaxy assembly bias also presents very small dependence upon cosmology ($\sim$ 2$\%$-4$\%$ of the total clustering); on the other hand, we find that the dependence of this signal on the galaxy formation parameters of our galaxy model is much stronger. Taken together, these results let us conclude that the dependence of halo and galaxy assembly bias on cosmology is practically negligible., Comment: 11 pages, 8 figures. Submitted to MNRAS, comments welcome

We present a neural-network emulator for baryonic effects in the non-linear matter power spectrum. We calibrate this emulator using more than 50,000 measurements in a 15-dimensional parameters space, varying cosmology and baryonic physics. Baryonic physics is described through a baryonification algorithm, that has been shown to accurately capture the relevant effects on the power spectrum and bispectrum in state-of-the-art hydrodynamical simulations. Cosmological parameters are sampled using a cosmology-rescaling approach including massive neutrinos and dynamical dark energy. The specific quantity we emulate is the ratio between matter power spectrum with baryons and gravity-only, and we estimate the overall precision of the emulator to be 1-2%, at all scales 0.01 < k < 5 h/Mpc, and redshifts 0 < z < 1.5. We also obtain an accuracy of 1-2%, when testing the emulator against a collection of 74 different cosmological hydrodynamical simulations and their respective gravity-only counterparts. We show also that only one baryonic parameter, namely Mc, which set the gas fraction retained per halo mass, is enough to have accurate and realistic predictions of the baryonic feedback at a given epoch. Our emulator will become publicly available in http://www.dipc.org/bacco., 13 pages, 11 figures. Comments are welcome

We use the {\sc Illustris TNG300} magneto-hydrodynamic simulation, the {\sc SAGE} semi-analytical model, and the subhalo abundance matching technique (SHAM) to examine the diversity in predictions for galaxy assembly bias (i.e. the difference in the large scale clustering of galaxies at a fixed halo mass due to correlations with the assembly history and other properties of host haloes). We consider samples of galaxies selected according to their stellar mass or star formation rate at various redshifts. We find that all models predict an assembly bias signal of different magnitude, redshift evolution, and dependence with selection criteria and number density. To model these non-trivial dependences, we propose an extension to the standard SHAM technique so it can include arbitrary amounts of assembly bias. We do this by preferentially selecting subhaloes with the same internal property but different {\it individual} large-scale bias. We find that with this model, we can successfully reproduce the galaxy assembly bias signal in either {\sc SAGE} or the {\sc TNG}, for all redshifts and galaxy number densities. We anticipate that this model can be used to constrain the level of assembly bias in observations and aid in the creation of more realistic mock galaxy catalogues., Accepted by mnras

Predicting the spatial distribution of objects as a function of cosmology is an essential ingredient for the exploitation of future galaxy surveys. In this paper we show that a specially-designed suite of gravity-only simulations together with cosmology-rescaling algorithms can provide the clustering of dark matter, haloes, and subhaloes with high precision. Specifically, with only 3 $N$-body simulations we obtain the power spectrum of dark matter at $z=0$ and $z=1$ to better than 3\% precision for essentially all currently viable values of 8 cosmological parameters, including massive neutrinos and dynamical dark energy, and over the whole range of scales explored, 0.03 < $k/h^{-1}Mpc$ < 5. This precision holds at the same level for mass-selected haloes and for subhaloes selected according to their peak maximum circular velocity. As an initial application of these predictions, we successfully constrain $\Omega_{\rm m}$, $\sigma_8$, and the scatter in subhalo-abundance-matching employing the projected correlation function of mock SDSS galaxies., Comment: 12 pages, 8 figures. Submitted. Comments are welcome

We compare predictions for galaxy-galaxy lensing profiles and clustering from the Henriques et al. (2015) public version of the Munich semi-analytical model of galaxy formation (SAM) and the IllustrisTNG suite, primarily TNG300, with observations from KiDS+GAMA and SDSS-DR7 using four different selection functions for the lenses (stellar mass, stellar mass and group membership, stellar mass and isolation criteria, stellar mass and colour). We find that this version of the SAM does not agree well with the current data for stellar mass-only lenses with $M_\ast > 10^{11}\,M_\odot$. By decreasing the merger time for satellite galaxies as well as reducing the radio-mode AGN accretion efficiency in the SAM, we obtain better agreement, both for the lensing and the clustering, at the high mass end. We show that the new model is consistent with the signals for central galaxies presented in Velliscig et al. (2017). Turning to the hydrodynamical simulation, TNG300 produces good lensing predictions, both for stellar mass-only ($��^2 = 1.81$ compared to $��^2 = 7.79$ for the SAM), and locally brightest galaxies samples ($��^2 = 3.80$ compared to $��^2 = 5.01$). With added dust corrections to the colours it matches the SDSS clustering signal well for red low mass galaxies. We find that both the SAMs and TNG300 predict $\sim 50\,\%$ excessive lensing signals for intermediate mass red galaxies with $10.2 < \log_{10} M_\ast [ M_\odot ] < 11.2$ at $r \approx 0.6\,h^{-1}\,\mathrm{Mpc}$, which require further theoretical development., 30 pages, 20 figures, 10 tables, minor revisions, matches version accepted by MNRAS

Fluctuations in sky maps of the galaxy redshifts, dubbed as angular redshift fluctuations (ARF), contain precise information about the growth rate of structures and the nature of gravity in the Universe. Specifically, ARF constrain the combination of cosmological parameters $H/H_0\,f\sigma_8(z)$, while being an intrinsically tomographic probe and largely insensitive to many observational systematic errors, all this without requiring the assumption of any redshift-to-distance relation under a given fiducial cosmology. We present the first cosmological constraints derived from ARF by using BOSS LOWZ+CMASS DR12 galaxy samples, obtaining 7\%-accurate constraints on $H/H_0 f\sigma_8(z)$ at more than 20 redshifts over the range $z \in [0.26,0.72]$. Our best-fitting value is $10\%$ larger, but compatible at the $1.4\sigma$ level, than the $\Lambda$CDM expectation set by {\it Planck} observations of the Cosmic Microwave Background (CMB) radiation. Our tomographic measurements, combined with these CMB data, provides one of the strongest constraints on the gravity index $\gamma$, $\gamma=0.44^{+0.09}_{-0.07}$, which lies within $2\sigma$ from the prediction of General Relativity ($\gamma_{\rm GR}\simeq 0.55$)., Comment: 5 pages, 3 figures, matches accepted version in MNRAS letters

We present a model for the remnants of haloes that have gone through an adiabatic tidal stripping process. We show that this model exactly reproduces the remnant of an NFW halo that is exposed to a slowly increasing isotropic tidal field and approximately for an anisotropic tidal field. The model can be used to predict the asymptotic mass loss limit for orbiting subhaloes, solely as a function of the initial structure of the subhalo and the value of the tidal field at pericentre. Predictions can easily be made for differently concentrated host-haloes with and without baryonic components, which differ most notably in their relation between pericentre radius and tidal field. The model correctly predicts several empirically measured relations such as the `tidal track' and the `orbital frequency relation' that was reported by Errani & Navarro (2021) for the case of an isothermal sphere. Further, we propose applications of the `structure-tide' degeneracy which implies that increasing the concentration of a subhalo has exactly the same impact on tidal stripping as reducing the amplitude of the tidal field. Beyond this, we find that simple relations hold for the bound mass, truncation radius, WIMP annihilation luminosity and tidal ratio of tidally stripped NFW haloes in relation to quantities measured at the radius of maximum circular velocity. Finally, we note that NFW haloes cannot be completely disrupted when exposed adiabatically to tidal fields of arbitrary magnitudes. We provide an open-source implementation of our model and suggest that it can be used to improve predictions of dark matter annihilation., 31 pages, 19 figures, submitted to MNRAS, comments welcome, code publicly available under https://github.com/jstuecker/adiabatic-tides

We present an emulator that accurately predicts the power spectrum of galaxies in redshift space as a function of cosmological parameters. Our emulator is based on a 2nd-order Lagrangian bias expansion that is displaced to Eulerian space using cosmological $N$-body simulations. Redshift space distortions are then imprinted using the non-linear velocity field of simulated particles and haloes. We build the emulator using a forward neural network trained with the simulations of the BACCO project, which covers an 8-dimensional parameter space including massive neutrinos and dynamical dark energy. We show that our emulator provides unbiased cosmological constraints from the monopole, quadrupole, and hexadecapole of a mock galaxy catalogue that mimics the BOSS-CMASS sample down to nonlinear scales ($k\sim0.6$[$h/$Mpc]$^{3}$). This work opens up the possibility of robustly extracting cosmological information from small scales using observations of the large-scale structure of the Universe., Comment: 16 pages, 9 figures

Current best limits on the 21-cm signal during reionization are provided at large scales ($\gtrsim$100 Mpc). To model these scales, enormous simulation volumes are required which are computationally expensive. We find that the primary source of uncertainty at these large scales is sample variance, which decides the minimum size of simulations required to analyse current and upcoming observations. In large-scale structure simulations, the method of `fixing' the initial conditions (ICs) to exactly follow the initial power spectrum and `pairing' two simulations with exactly out-of-phase ICs has been shown to significantly reduce sample variance. Here we apply this `fixing and pairing' (F\&P) approach to reionization simulations whose clustering signal originates from both density fluctuations and reionization bubbles. Using a semi-numerical code, we show that with the traditional method, simulation boxes of $L\simeq 500$ (300) Mpc are required to model the large-scale clustering signal at $k$=0.1 Mpc$^{-1}$ with a precision of 5 (10) per cent. Using F\&P, the simulation boxes can be reduced by a factor of 2 to obtain the same precision level. We conclude that the computing costs can be reduced by at least a factor of 4 when using the F\&P approach., 5 pages, 4 figures, accepted for publication in A&A

The global gravitational potential, $\phi$, is not commonly employed in the analysis of cosmological simulations, since its level sets do not show any clear correspondence to the underlying density field and its persistent structures. Here, we show that the potential becomes a locally meaningful quantity when considered from a boosted frame of reference, defined by subtracting a uniform gradient term $\phi_{\rm{boost}}(\boldsymbol{x}) = \phi(\boldsymbol{x}) + \boldsymbol{x} \cdot \boldsymbol{a}_0$ with acceleration $\boldsymbol{a}_0$. We study this "boosted potential" in a variety of scenarios and propose several applications: (1) The boosted potential can be used to define a binding criterion that naturally incorporates the effect of tidal fields. This solves several problems of commonly-used self-potential binding checks: i) it defines a tidal boundary for each halo, ii) it is much less likely to misidentify caustics as haloes (specially in the context of warm dark matter cosmologies), and iii) performs better at identifying virialized regions of haloes -- yielding to the expected value of 2 for the virial ratio. (2) This binding check can be generalized to filaments and other cosmic structures. (3) The boosted potential facilitates the understanding of the disruption of satellite subhaloes. We propose a picture where most mass loss is explained through a lowering of the escape energy through the tidal field. (4) We discuss the possibility of understanding the topology of the potential field in a way that is independent of constant offsets in the first derivative $\boldsymbol{a}_0$. We foresee that this novel perspective on the potential can help to develop more accurate models and improve our understanding of structure formation., Comment: 22 pages, 14 figures. Submitted to MNRAS. Comments welcome!

Improving the theoretical description of galaxy clustering on small scales is an important challenge in cosmology, as it can considerably increase the scientific return of forthcoming galaxy surveys -- e.g. tightening the bounds on neutrino masses and deviations from general relativity. In this paper, we propose and test a new model for the clustering of galaxies that is able to accurately describe redshift-space distortions even down to small scales. This model corresponds to a second-order perturbative Lagrangian bias expansion which is advected to Eulerian space employing a displacement field extracted from $N$-body simulations. Eulerian coordinates are then transformed into redshift space by directly employing simulated velocity fields augmented with nuisance parameters capturing various possible satellite fractions and intra-halo small-scale velocities. We quantify the accuracy of our approach against samples of physically-motivated mock galaxies selected according to either Stellar Mass (SM) or Star Formation Rate (SFR) at multiple abundances and at $z=0$ and $1$. We find our model describes the monopole, quadrupole, and hexadecapole of the galaxy-power spectra down to scales of $k\approx 0.6 [h/$Mpc] within the accuracy of our simulations. This approach could pave the way to significantly increase the amount of cosmological information to be extracted from future galaxy surveys., Comment: 15 pages, 8 figures

Baryons and cold dark matter (CDM) did not comove prior to recombination. This leads to differences in the local baryon and CDM densities, the so-called baryon-CDM isocurvature perturbations $\delta_{bc}$. These perturbations are usually neglected in the analysis of Large-Scale Structure data but taking them into account might become important in the era of high precision cosmology. Using gravity-only 2-fluid simulations we assess the impact of such perturbations on the dark matter halos distribution. In particular, we focus on the baryon fraction in halos as a function of mass and large-scale $\delta_{bc}$, which also allows us to study details of the nontrivial numerical setup required for such simulations. We further measure the cross-power spectrum between the halo field and $\delta_{bc}$ over a wide range of mass. This cross-correlation is nonzero and negative which shows that halo formation is impacted by $\delta_{bc}$. We measure the associated bias parameter $b_{\delta_{bc}}$ and compare it to recent results, finding good agreement. Finally we quantify the impact of such perturbations on the halo-halo power spectrum and show that this effect can be degenerate with the one of massive neutrinos for surveys like DESI., Comment: 27 pages, 11 figures, Version published in JCAP

This is the second part of a thorough investigation of the redshift-space effects that affect void properties and the impact they have on cosmological tests. Here, we focus on the void-galaxy cross-correlation function, specifically, on the projected versions that we developed in a previous work. The pillar of the analysis is the one-to-one relationship between real and redshift-space voids above the shot-noise level identified with a spherical void finder. Under this mapping, void properties are affected by three effects: (i) a systematic expansion as a consequence of the distortions induced by galaxy dynamics, (ii) the Alcock-Paczynski volume effect, which manifests as an overall expansion or contraction depending on the fiducial cosmology, and (iii) a systematic off-centring along the line of sight as a consequence of the distortions induced by void dynamics. We found that correlations are also affected by an additional source of distortions: the ellipticity of voids. This is the first time that distortions due to the off-centring and ellipticity effects are detected and quantified. With a simplified test, we verified that the Gaussian streaming model is still robust provided all these effects are taken into account, laying the foundations for improvements in current models in order to obtain unbiased cosmological constraints from spectroscopic surveys. Besides this practical importance, this analysis also encodes key information about the structure and dynamics of the Universe at the largest scales. Furthermore, some of the effects constitute cosmological probes by themselves, as is the case of the void ellipticity., Comment: 14 pages, 9 figures, published in MNRAS, accepted after minor comments

We demonstrate that baryonification algorithms, which displace particles in gravity-only simulations according to physically-motivated prescriptions, can simultaneously capture the impact of baryonic physics on the 2 and 3-point statistics of matter. Specifically, we show that our implementation of a baryonification algorithm jointly fits the changes induced by baryons on the power spectrum and equilateral bispectrum on scales up to k < 5 h/Mpc and redshifts z, 14 pages, 15 figures. To be submitted to MNRAS. Comments are welcome

In the context of next-generation spectroscopic galaxy surveys, new statistics of the distribution of matter are currently being developed. Among these, we investigated the angular redshift fluctuations (ARF), which probe the information contained in the projected redshift distribution of galaxies. Relying on the Fisher formalism, we show how ARF will provide complementary cosmological information compared to traditional angular galaxy clustering. We tested both the standard $\mathrm{\Lambda CDM}$ model and the $\mathrm{wCDM}$ extension. We find that the cosmological and galaxy bias parameters express different degeneracies when inferred from ARF or from angular galaxy clustering. As such, combining both observables breaks these degeneracies and greatly decreases the marginalised uncertainties by a factor of at least two on most parameters for the $\mathrm{\Lambda CDM}$ and $\mathrm{wCDM}$ models. We find that the ARF combined with angular galaxy clustering provide a great way to probe dark energy by increasing the figure of merit of the $w_0$-$w_{\rm a}$ parameter set by a factor of more than ten compared to angular galaxy clustering alone. Finally, we compared ARF to the CMB lensing constraints on the galaxy bias parameters. We show that a joint analysis of ARF and angular galaxy clustering improves constraints by $\sim 40\%$ on galaxy bias compared to a joint analysis of angular galaxy clustering and CMB lensing., Comment: 14 pages, 21 figures, published in A&A

Full list of authors: Bonoli, S.; Marín-Franch, A.; Varela, J.; Vázquez Ramió, H.; Abramo, L. R.; Cenarro, A. J.; Dupke, R. A.; Vílchez, J. M.; Cristóbal-Hornillos, D.; González Delgado, R. M.; Hernández-Monteagudo, C.; López-Sanjuan, C.; Muniesa, D. J.; Civera, T.; Ederoclite, A.; Hernán-Caballero, A.; Marra, V.; Baqui, P. O.; Cortesi, A.; Cypriano, E. S.; Daflon, S.; de Amorim, A. L.; Díaz-García, L. A.; Diego, J. M.; Martínez-Solaeche, G.; Pérez, E.; Placco, V. M.; Prada, F.; Queiroz, C.; Alcaniz, J.; Alvarez-Candal, A.; Cepa, J.; Maroto, A. L.; Roig, F.; Siffert, B. B.; Taylor, K.; Benitez, N.; Moles, M.; Sodré, L.; Carneiro, S.; Mendes de Oliveira, C.; Abdalla, E.; Angulo, R. E.; Aparicio Resco, M.; Balaguera-Antolínez, A.; Ballesteros, F. J.; Brito-Silva, D.; Broadhurst, T.; Carrasco, E. R.; Castro, T.; Cid Fernandes, R.; Coelho, P.; de Melo, R. B.; Doubrawa, L.; Fernandez-Soto, A.; Ferrari, F.; Finoguenov, A.; García-Benito, R.; Iglesias-Páramo, J.; Jiménez-Teja, Y.; Kitaura, F. S.; Laur, J.; Lopes, P. A. A.; Lucatelli, G.; Martínez, V. J.; Maturi, M.; Overzier, R. A.; Pigozzo, C.; Quartin, M.; Rodríguez-Martín, J. E.; Salzano, V.; Tamm, A.; Tempel, E.; Umetsu, K.; Valdivielso, L. ; von Marttens, R.; Zitrin, A.; Díaz-Martín, M. C.; López-Alegre, G.; López-Sainz, A.; Yanes-Díaz, A.; Rueda-Teruel, F.; Rueda-Teruel, S.; Abril Ibañez, J.; L Antón Bravo, J.; Bello Ferrer, R.; Bielsa, S.; Casino, J. M.; Castillo, J.; Chueca, S.; Cuesta, L.; Garzarán Calderaro, J.; Iglesias-Marzoa, R.; Íniguez, C.; Lamadrid Gutierrez, J. L.; Lopez-Martinez, F.; Lozano-Pérez, D.; Maícas Sacristán, N.; Molina-Ibáñez, E. L.; Moreno-Signes, A.; Rodríguez Llano, S.; Royo Navarro, M.; Tilve Rua, V.; Andrade, U.; Alfaro, E. J.; Akras, S.; Arnalte-Mur, P.; Ascaso, B.; Barbosa, C. E.; Beltrán Jiménez, J.; Benetti, M.; Bengaly, C. A. P.; Bernui, A.; Blanco-Pillado, J. J.; Borges Fernandes, M.; Bregman, J. N.; Bruzual, G.; Calderone, G.; Carvano, J. M.; Casarini, L.; Chaves-Montero, J.; Chies-Santos, A. L.; Coutinho de Carvalho, G.; Dimauro, P.; Duarte Puertas, S.; Figueruelo, D.; González-Serrano, J. I.; Guerrero, M. A.; Gurung-López, S.; Herranz, D.; Huertas-Company, M.; Irwin, J. A.; Izquierdo-Villalba, D.; Kanaan, A.; Kehrig, C.; Kirkpatrick, C. C.; Lim, J.; Lopes, A. R.; Lopes de Oliveira, R.; Marcos-Caballero, A.; Martínez-Delgado, D.; Martínez-González, E.; Martínez-Somonte, G.; Oliveira, N.; Orsi, A. A.; Penna-Lima, M.; Reis, R. R. R.; Spinoso, D.; Tsujikawa, S.; Vielva, P.; Vitorelli, A. Z.; Xia, J. Q.; Yuan, H. B.; Arroyo-Polonio, A.; Dantas, M. L. L.; Galarza, C. A.; Gonçalves, D. R.; Gonçalves, R. S.; Gonzalez, J. E.; Gonzalez, A. H.; Greisel, N.; Jiménez-Esteban, F.; Landim, R. G.; Lazzaro, D.; Magris, G.; Monteiro-Oliveira, R.; Pereira, C. B.; Rebouças, M. J.; Rodriguez-Espinosa, J. M.; Santos da Costa, S.; Telles, E., The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will scan thousands of square degrees of the northern sky with a unique set of 56 filters using the dedicated 2:55m Javalambre Survey Telescope (JST) at the Javalambre Astrophysical Observatory. Prior to the installation of the main camera (4:2 deg2 field-of-view with 1.2 Gpixels), the JST was equipped with the JPAS-Pathfinder, a one CCD camera with a 0:3 deg2 field-of-view and plate scale of 0.23 arcsec pixel?1. To demonstrate the scientific potential of J-PAS, the JPAS-Pathfinder camera was used to perform miniJPAS, a _1 deg2 survey of the AEGIS field (along the Extended Groth Strip). The field was observed with the 56 J-PAS filters, which include 54 narrow band (FWHM _ 145 ) and two broader filters extending to the UV and the near-infrared, complemented by the u; g; r; i SDSS broad band filters. In this miniJPAS survey overview paper, we present the miniJPAS data set (images and catalogs), as we highlight key aspects and applications of these unique spectro-photometric data and describe how to access the public data products. The data parameters reach depths of magAB ' 22?23:5 in the 54 narrow band filters and up to 24 in the broader filters (5_ in a 300 aperture). The miniJPAS primary catalog contains more than 64 000 sources detected in the r band and with matched photometry in all other bands. This catalog is 99% complete at r = 23:6 (r = 22:7) mag for point-like (extended) sources. We show that our photometric redshifts have an accuracy better than 1% for all sources up to r = 22:5, and a precision of _0:3% for a subset consisting of about half of the sample. On this basis, we outline several scientific applications of our data, including the study of spatially-resolved stellar populations of nearby galaxies, the analysis of the large scale structure up to z _ 0:9, and the detection of large numbers of clusters and groups. Sub-percent redshift precision can also be reached for quasars, allowing for the study of the large-scale structure to be pushed to z 2. The miniJPAS survey demonstrates the capability of the J-PAS filter system to accurately characterize a broad variety of sources and paves the way for the upcoming arrival of J-PAS, which will multiply this data by three orders of magnitude. © 2021 EDP Sciences. All rights reserved., Funding for OAJ, UPAD, and CEFCA has been provided by the Governments of Spain and Aragon through the Fondo de Inversiones de Teruel; the Aragon Government through the Research Groups E96, E103, and E16_17R; the Spanish Ministry of Science, Innovation and Universities (MCIU/AEI/FEDER, UE) with grant PGC2018-097585-B-C21; the Spanish Ministry of Economy and Competitiveness (MINECO/FEDER, UE) under AYA2015-66211-C2-1-P, AYA2015-66211-C2-2, AYA2012-30789, and ICTS2009-14; and European FEDER funding (FCDD10-4E-867, FCDD13-4E-2685). This work has made used of CEFCA's Scientific High Performance Computing system which has been funded by the Governments of Spain and Aragon through the Fondo de Inversiones de Teruel, and the Spanish Ministry of Economy and Competitiveness (MINECO-FEDER, grant AYA2012-30789). Funding for the J-PAS project has been provided also by the Brazilian agencies FINEP, FAPESP, FAPERJ and by the National Observatory of Brazil. Additional funding was also provided by the Tartu Observatory and by the Chinese Consortium from the Academy of Sciences SB acknowledges partial support from the project PGC2018-097585-B-C22. R.A.D. acknowledges support from the Conselho Nacional de Desenvolvimento Cientifico e Tecnologico -CNPq through BP grant 308105/2018-4, and the Financiadora de Estudos e Projetos -FINEP grants REF. 1217/13 -01.13.0279.00 and REF 0859/10 -01.10.0663.00 and also FAPERJ PRONEX grant E-26/110.566/2010 for hardware funding support for the J-PAS project through the National Observatory of Brazil and Centro Brasileiro de Pesquisas Fisicas. LRA acknowledges financial support from CNPq (306696/2018-5) and FAPESP (2015/17199-0). VMthanks CNPq (Brazil) and FAPES (Brazil) for partial financial support and the H2020 project No 888258. L.A.D.G. and K.U. acknowledge support from the Ministry of Science and Technology of Taiwan (grant MOST 106-2628-M-001-003-MY3) and from the Academia Sinica (grant AS-IA-107-M01). J.M.D. and D.H acknowledge the support of project PGC2018-101814-B-100. MQ thanks CNPq (Brazil) and FAPERJ (Brazil) for financial support. PC acknowledges financial support from FAPESP (2018/05392-8) and CNPq (310041/2018-0). AAC acknowledges support from FAPERJ (E26/203.186/2016), CNPq (304971/2016-2 and 401669/2016-5), and the Universidad de Alicante (contract UATALENTO1802). C.Q. acknowledges support from FAPESP (2015/11442-0 and 2019/067661). V.M.P. is supported by NOIRLab, which is managed by AURA under a cooperative agreement with the NSF. P.B acknowledges support from CAPES -Finance Code 001. IAA researchers acknowledge financial support from the State Agency for Research of the Spanish MCIU through the "Center of Excellence Severo Ochoa" award to the Instituto de Astrofisica de Andalucia (SEV-2017-0709). Authors acknowledge support from the Generalitat Valenciana project of excellence Prometeo/2020/085. RGD, GMS, JRM, RGB, EP acknowledge financial support from the project AyA2016-77846-P. TC is supported by the INFN INDARK PD51 and PRIN-MIUR 2015W7KAWC. MAR and ALM acknowledge support from the MINECO project FIS2016-78859P(AEI/FEDER, UE). ET, AT and JL acknowledge the support by ETAg grants IUT40-2 and by EU through the ERDF CoE grant TK133 and MOBTP86. CK, JMV, JIP acknowledge financial support from project AYA2016-79724C4-4P. PAAL thanks the support of CNPq (309398/2018-5). LC thanks CNPq for partial support. Y.J-T acknowledges financial support from the FAPERJ (E26/202.835/2016), and from the Horizon 2020 Marie Sklodowska-Curie grant agreement No 898633. DMD acknowledges financial support from the SFB 881 of the DFG and from the MINECO grant AYA2016-81065-C2-2. FP acknowledges support of the project PGC2018-101931-B-I00. JC acknowledges support of the project E AYA2017-88007-C3-1-P, and co-financed by the FEDER. JIGs acknowledges support of projects of reference AYA2017-88007-C3-3-P, and PGC2018-099705-B-I00 and co-financed by the FEDER. EMG and PV would like to acknowledge financial support from the project ESP2017-83921C2-1-R. GMS acknowleges financial support from a predoctoral contract, ref. PRE2018-085523 (MCIU/AEI/FSE, UE). S.C. is partially supported by CNPq. R.G.L. acknowledges CAPES (process 88881.162206/2017-01) and Alexander von Humboldt Foundation for the financial support. JSA acknowledges support from FAPERJ (E26/203.024/2017), CNP (310790/2014-0 and 400471/2014-0) and FINEP (1217/13 -01.13.0279.00 and Ref. 0859/10 -01.10.0663.00). RvM acknowledges support from CNPq. AFS, PAM, VJM and FJB acknowledge support from project AYA2016-81065-C2-2. PAM acknowledges support from the "Subprograma Atraccio de Talent -Contractes Postdoctorals de la Universitat de Valencia". ESC acknowledges support from CNPq (308539/20184) and FAPESP (2019/19687-2). CMdO acknowledges support from CNPq (grant 312333/2014-5) and FAPESP (grant 2009/54202-8). LSJ acknowledges support from CNPq (grant 304819/2017-4) and FAPESP (grant 2012/008004). JMC acknowledges support from CNPq (grant 310727/2016-2). C.H.-M. and N. Greisel also acknowledge the support of the European Union via the Career Integration Grant CIG-PCIG9-GA-2011-294183. JJBP and AMC would like to acknowledge the support from the grant PGC2018-094626-B-C21 and the Basque Government grant IT-979-16. AMC acknowledges the postdoctoral contract from the University of the Basque Country UPV/EHU "Especializacioon de personal investigador doctor" program. MLLD acknowledges CAPES -Finance Code 001; and CNPq (142294/2018-7). GB acknowledges financial support from the UNAM through grant DGAPA/PAPIIT IG100319, from CONACyT through grant CB2015-252364, and from FAPESP projects 2017/02375-2 and 2018/05392-8. M.J. Reboucas acknowledges the support of FAPERJ under a CNE E-26/202.864/2017 grant, and CNPq. Support by CNPq (305409/2016-6) and FAPERJ (E-26/202.841/2017) is acknowledged by DL. AB acknowledges a CNPq fellowship. C.A.G.acknowledges support from CAPES. EA acknowledges support from FAPESP (2011/18729-1). AC acknowledges support from PNPD/CAPES. ABA and FSK acknowledge the Severo Ochoa program SEV-2015-0548. FSK also thanks the AYA2017-89891-P and the RYC2015-18693 grants. DF acknowledges support from the Atraccion del Talento Cientifico en Salamanca programme and the project PGC2018-096038B-I00.

We present the BACCO project, a simulation framework specially designed to provide highly-accurate predictions for the distribution of mass, galaxies, and gas as a function of cosmological parameters. In this paper, we describe our main suite of simulations (L $\sim2$ Gpc and $4320^3$ particles) and present various validation tests. Using a cosmology-rescaling technique, we predict the nonlinear mass power spectrum over the redshift range $0, Comment: 12 pages; 11 figures. Accepted for publication in MNRAS

We use the improved IllustrisTNG300 magneto-hydrodynamical cosmological simulation to revisit the effect that secondary halo bias has on the clustering of the central galaxy population. With a side length of 205 $h^{-1}$Mpc and significant improvements on the sub-grid model with respect to the previous Illustris boxes, IllustrisTNG300 allows us to explore the dependencies of galaxy clustering over a large cosmological volume and wide halo-mass range. We show, at high statistical significance, that the halo assembly bias signal (i.e., the secondary dependence of halo bias on halo formation redshift) manifests itself on the clustering of the central galaxy population when this is split by stellar mass, colour, specific star formation rate, and surface density. A significant detection is also obtained for galaxy size: at fixed halo mass, larger central galaxies are more tightly clustered than smaller central galaxies in haloes of mass M$_{\rm vir} \lesssim 10^{12.5}$ $h^{-1}$M$_{\odot}$. This effect, however, seems to be uncorrelated with halo formation time, unlike the rest of the secondary dependencies analysed. We also explore the transmission of the halo spin bias signal, i.e., the secondary dependence of halo bias on halo spin. Although galaxy spin retains little information about the total spin of the halo, the correlation is enough to produce a significant galaxy spin bias signal. We discuss possible ways to probe the spin bias effects with observations., 16 pages, 10 figures, submitted to MNRAS

At early times dark matter has a thermal velocity dispersion of unknown amplitude which, for warm dark matter models, can influence the formation of nonlinear structure on observable scales. We propose a new scheme to simulate cosmologies with a small-scale suppression of perturbations that combines two previous methods in a way that avoids the numerical artefacts which have so far prevented either from producing fully reliable results. At low densities and throughout most of the cosmological volume, we represent the dark matter phase-sheet directly using high-accuracy interpolation, thereby avoiding the artificial fragmentation which afflicts particle-based methods in this regime. Such phase-sheet methods are, however, unable to follow the rapidly increasing complexity of the denser regions of dark matter haloes, so for these we switch to an N-body scheme which uses the geodesic deviation equation to track phase-sheet properties local to each particle. In addition, we present a novel high-resolution force calculation scheme based on an oct-tree of cubic force resolution elements which is well suited to approximate the force-field of our combined sheet+particle distribution. Our hybrid simulation scheme enables the first reliable simulations of the internal structure of low-mass haloes in a warm dark matter cosmology., 23 pages, 28 figures

The interpretation of cosmological observables requires the use of increasingly sophisticated theoretical models. Since these models are becoming computationally very expensive and display non-trivial uncertainties, the use of standard Bayesian algorithms for cosmological inferences, such as MCMC, might become inadequate. Here, we propose a new approach to parameter estimation based on an iterative Gaussian emulation of the target likelihood function. This requires a minimal number of likelihood evaluations and naturally accommodates for stochasticity in theoretical models. We apply the algorithm to estimate 9 parameters from the monopole and quadrupole of a mock power spectrum in redshift space. We obtain accurate posterior distribution functions with approximately 100 times fewer likelihood evaluations than an affine invariant MCMC, roughly independently from the dimensionality of the problem. We anticipate that our parameter estimation algorithm will accelerate the adoption of more accurate theoretical models in data analysis, enabling more comprehensive exploitation of cosmological observables., 12 pages, 5 figures

We develop an extension of subhalo abundance matching (SHAM) capable of accurately reproducing the real and redshift-space clustering of galaxies in a state-of-the-art hydrodynamical simulation. Our method uses a low-resolution gravity-only simulation and it includes orphan and tidal disruption prescriptions for satellite galaxies, and a flexible amount of galaxy assembly bias. Furthermore, it includes recipes for star formation rate (SFR) based on the dark matter accretion rate. We test the accuracy of our model against catalogues of stellar-mass- and SFR-selected galaxies in the TNG300 hydrodynamic simulation. By fitting a small number of free parameters, our extended SHAM reproduces the projected correlation function and redshift-space multipoles for number densities $10^{-3} - 10^{-2}\, h^{3}{\rm Mpc}^{-3}$, at $z=1$ and $z=0$, and for scales $r \in [0.3 - 20] h^{-1}{\rm Mpc}$. Simultaneously, the SHAM results also retrieve the correct halo occupation distribution, the level of galaxy assembly bias, and higher-order statistics present in the TNG300 galaxy catalogues. As an application, we show that our model simultaneously fits the projected correlation function of the SDSS in 3 disjoint stellar mass bins, with an accuracy similar to that of TNG300 galaxies. This SHAM extension can be used to get accurate clustering prediction even when using low and moderate-resolution simulations., Comment: 14 pages, 12 figures. Submitted to MNRAS

Inaccuracies in the initial conditions for cosmological N-body simulations could easily be the largest source of systematic error in predicting the non-linear large-scale structure. From the theory side, initial conditions are usually provided by using low-order truncations of the displacement field from Lagrangian perturbation theory, with the first and second-order approximations being the most common ones. Here we investigate the improvement brought by using initial conditions based on third-order Lagrangian perturbation theory (3LPT). We show that with 3LPT, truncation errors are vastly suppressed, thereby opening the portal to initializing simulations accurately as late as z=12 (for the resolution we consider). We analyse the competing effects of perturbative truncation and particle discreteness on various summary statistics. Discreteness errors are essentially decaying modes and thus get strongly amplified for earlier initialization times. We show that late starting times with 3LPT provide the most accurate configuration, which we find to coincide with the continuum fluid limit within 1 per cent for the power- and bispectrum at z=0 up to the particle Nyquist wave number of our simulations (k ~ 3h/Mpc). In conclusion, to suppress non-fluid artefacts, we recommend initializing simulations as late as possible with 3LPT. We make our 3LPT initial condition generator publicly available., Comment: 22 pages, 13 figures, minor revisions (matches published version), code available from https://bitbucket.org/ohahn/monofonic

In the upcoming era of high-precision galaxy surveys, it becomes necessary to understand the impact of redshift uncertainties on cosmological observables. In this paper we explore the effect of sub-percent photometric redshift errors (photo-$z$ errors) on galaxy clustering and baryonic acoustic oscillations (BAO). Using analytic expressions and results from $1\,000$ $N$-body simulations, we show how photo-$z$ errors modify the amplitude of moments of the 2D power spectrum, their variances, the amplitude of BAO, and the cosmological information in them. We find that: a) photo-$z$ errors suppress the clustering on small scales, increasing the relative importance of shot noise, and thus reducing the interval of scales available for BAO analyses; b) photo-$z$ errors decrease the smearing of BAO due to non-linear redshift-space distortions (RSD) by giving less weight to line-of-sight modes; and c) photo-$z$ errors (and small-scale RSD) induce a scale dependence on the information encoded in the BAO scale, and that reduces the constraining power on the Hubble parameter. Using these findings, we propose a template that extracts unbiased cosmological information from samples with photo-$z$ errors with respect to cases without them. Finally, we provide analytic expressions to forecast the precision in measuring the BAO scale, showing that spectro-photometric surveys will measure the expansion history of the Universe with a precision competitive to that of spectroscopic surveys., Comment: 19 pages, 15 figures, resubmission. Included full analysis of power spectrum moments l=0, 2, and 4

Galaxy surveys aim to map the large-scale structure of the Universe and use redshift space distortions to constrain deviations from general relativity and probe the existence of massive neutrinos. However, the amount of information that can be extracted is limited by the accuracy of theoretical models used to analyze the data. Here, by using the L-Galaxies semi-analytical model run over the MXXL N-body simulation, we assess the impact of galaxy formation on satellite kinematics and the theoretical modelling of redshift-space distortions. We show that different galaxy selection criteria lead to noticeable differences in the radial distributions and velocity structure of satellite galaxies. Specifically, whereas samples of stellar mass selected galaxies feature satellites that roughly follow the dark matter, emission line satellite galaxies are located preferentially in the outskirts of halos and display net infall velocities. We demonstrate that capturing these differences is crucial for modelling the multipoles of the correlation function in redshift space, even on large scales. In particular, we show how modelling small scale velocities with a single Gaussian distribution leads to a poor description of the measure clustering. In contrast, we propose a parametrization that is flexible enough to model the satellite kinematics, and that leads to and accurate description of the correlation function down to sub-Mpc scales. We anticipate that our model will be a necessary ingredient in improved theoretical descriptions of redshift space distortions, which together could result in significantly tighter cosmological constraints and a more optimal exploitation of future large datasets., 16 pages, 15 figures, typos corrected and moderate revision after referee report

Future galaxy surveys require realistic mock catalogues to understand and quantify systematics in order to make precise cosmological measurements. We present a halo lightcone catalogue and halo occupation distribution (HOD) galaxy catalogue built using the Millennium-XXL (MXXL) simulation. The halo catalogue covers the full sky, extending to z = 2 with a mass resolution of ~1e11 Msun/h . We use this to build a galaxy catalogue, which has an r-band magnitude limit of r < 20.0, with a median redshift of z~0.2. A Monte Carlo HOD method is used to assign galaxies to the halo lightcone catalogue, and we evolve the HODs to reproduce a target luminosity function; by construction, the luminosity function of galaxies in the mock is in agreement with the Sloan Digital Sky Survey (SDSS) at low redshifts and the Galaxy and Mass Assembly (GAMA) survey at high redshifts. A Monte Carlo method is used to assign a 0.1(g-r) colour to each galaxy, and the colour distribution of galaxies at different redshifts agrees with measurements from GAMA. The clustering of galaxies in the mock for galaxies in different magnitude and redshift bins is in good agreement with measurements from SDSS and GAMA, and the colour-dependent clustering is in reasonable agreement. We show that the baryon acoustic oscillation (BAO) can be measured in the mock catalogue, and the redshift space distortions (RSDs) are in agreement with measurements from SDSS, illustrating that this catalogue will be useful for upcoming surveys., 18 pages, 18 figures, 1 table, accepted for publication in MNRAS. Halo and galaxy catalogues are available at http://icc.dur.ac.uk/data/ and https://tao.asvo.org.au/tao/

We study the quenching of star formation as a function of redshift, environment and stellar mass in the galaxy formation simulations of Henriques et al. (2015), which implement an updated version of the Munich semi-analytic model (L-GALAXIES) on the two Millennium Simulations after scaling to a Planck cosmology. In this model, massive galaxies are quenched by active galactic nucleus (AGN) feedback depending on both black hole and hot gas mass, and hence indirectly on stellar mass. In addition, satellite galaxies of any mass can be quenched by ram-pressure or tidal stripping of gas and through the suppression of gaseous infall. This combination of processes produces quenching efficiencies which depend on stellar mass, host halo mass, environment density, distance to group centre and group central galaxy properties in ways which agree qualitatively with observation. Some discrepancies remain in dense regions and close to group centres, where quenching still seems too efficient. In addition, although the mean stellar age of massive galaxies agrees with observation, the assumed AGN feedback model allows too much ongoing star formation at late times. The fact that both AGN feedback and environmental effects are stronger in higher density environments leads to a correlation between the quenching of central and satellite galaxies which roughly reproduces observed conformity trends inside haloes., Monthly Notices of the Royal Astronomical Society, 469 (3), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

We propose the use of angular fluctuations in the galaxy redshift field as a new way to extract cosmological information in the Universe. This new probe $\delta z (\hat{n})$ consists on the statistics of sky maps built by projecting redshifts under a Gaussian window of width $\sigma_z$ centred upon a redshift $z_{\rm obs}$, and weighted by the galaxy density field. We compute the angular power spectrum of the $\delta z (\hat{n})$ field in both numerical simulations and in linear perturbation theory. From these we find that the $\delta z (\hat{n})$ field: {\it (i)} is sensitive to the underlying density and peculiar velocity fields; {\it (ii)} is highly correlated, at the $\gtrsim 60\,\%$ level, to the line-of-sight projected peculiar velocity field; {\it (iii)} for narrow windows $(\sigma_z < 0.03$), it is almost completely uncorrelated to the projected galaxy angular density field under the same redshift window; and {\it (iv)} it is largely unaffected by multiplicative and additive systematic errors on the observed number of galaxies that are redshift-independent over $\sim\sigma_z$. We conclude that $\delta z (\hat{n})$ is a simple and robust tomographic measure of the cosmic density and velocity fields, complementary to angular clustering, that will contribute to more complete exploitations of current and upcoming galaxy redshift surveys., Comment: 6 pages, 3 figures, matches accepted version in MNRAS letters, companion paper of Chaves-Montero et al., arXiv:1911.10690

A complete census of baryons in the late universe is a long-standing challenge due to the intermediate temperate and rarefied character of the majority of cosmic gas. To gain insight into this problem, we extract measurements of the kinematic Sunyaev-Zel'dovich (kSZ) effect from the cross-correlation of angular redshift fluctuations maps, which contain precise information about the cosmic density and velocity fields, and CMB maps high-pass filtered using aperture photometry; we refer to this technique as ARF-kSZ tomography. Remarkably, we detect significant cross-correlation for a wide range of redshifts and filter apertures using 6dF galaxies, BOSS galaxies, and SDSS quasars as tracers, yielding a 11 sigma detection of the kSZ effect. We then leverage these measurements to set constraints on the location, density, and abundance of gas inducing the kSZ effect, finding that this gas resides outside dark matter haloes, presents densities ranging from 10 to 250 times the cosmic average, and comprises half of cosmic baryons. Taken together, these findings indicate that ARF-kSZ tomography provides a nearly complete census of intergalactic gas from z=0 to 5., Comment: 17 pages, 14 figures, accepted for publication in MNRAS

We present a new cosmological test using the distribution of galaxies around cosmic voids without assuming a fiducial cosmology. The test is based on a physical model for the void-galaxy cross-correlation function projected along and perpendicular to the line of sight. We treat correlations in terms of void-centric angular distances and redshift differences between void-galaxy pairs, hence it is not necessary to assume a fiducial cosmology. This model reproduces the coupled dynamical (Kaiser effect, RSD) and geometrical (Alcock-Paczynski effect, GD) distortions that affect the correlation measurements. It also takes into account the scale mixing due to the projection ranges in both directions. The model is general, so it can be applied to an arbitrary cylindrical binning scheme, not only in the case of the projected correlations. It primarily depends on two cosmological parameters: $\Omega_m$, the matter fraction of the Universe today (sensitive to GD), and $\beta$, the ratio between the growth rate factor of density perturbations and the tracer bias (sensitive to RSD). In the context of the new generation of galaxy spectroscopic surveys, we calibrated the test using the Millennium XXL simulation for different redshifts. The method successfully recovers the cosmological parameters. We studied the effect of measuring with different projection ranges, finding robust results up to wide ranges. The resulting data covariance matrices are relatively small, which reduces the noise in the Gaussian likelihood analysis and will allow the usage of a smaller number of mock catalogues. The performance evaluated in this work indicates that the developed method is a promising test to be applied on real data., Comment: 13 pages. 7 Figures. Accepted in MNRAS

We use a suite of high-resolution $N$-body simulations and state-of-the-art perturbation theory to improve the code halofit, which predicts the nonlinear matter power spectrum. We restrict attention to parameters in the vicinity of the Planck Collaboration's best fit. On large-scales ($k\lesssim 0.07 h/{\rm Mpc}$), our model evaluates the 2-loop calculation from the Multi-point Propagator Theory of Bernardeau et al.(2012). On smaller scales ($k \gtrsim 0.7 h/{\rm Mpc}$), we transition to a smoothing-spline-fit model, that characterises the differences between the Takahashi et al. (2012) recalibration of halofit2012 and our simulations. We use an additional suite of simulations to explore the response of the power spectrum to variations in the cosmological parameters. In particular, we examine: the time evolution of the dark energy equation of state ($w_0$, $w_a$); the matter density $\Omega_m$; the physical densities of CDM and baryons $(\omega_c,\omega_b)$; and the primordial power spectrum amplitude $A_s$, spectral index $n_s$, and its running $\alpha$. We construct correction functions, which improve halofit's dependence on cosmological parameters. Our newly calibrated model reproduces all of our data with $\lesssim1\%$ precision. Including various systematic errors, such as choice of $N$-body code, resolution, and through inspection of the scaled second order derivatives, we estimate the accuracy to be $\lesssim3\%$ over the hyper-cube: $w_0\in\{-1.05,-0.95\}$, $w_a\in\{-0.4,0.4\}$, $\Omega_{\rm m,0}\in\{0.21,0.4\}$, $\omega_{\rm c}\in\{0.1,0.13\}$, $\omega_{\rm b}\in\{2.0,2.4\}$, $n_{\rm s}\in\{0.85,1.05\}$, $A_s\in\{1.72\times 10^{-9},2.58\times 10^{-9}\}$, $\alpha\in\{-0.2,0.2\}$ up to $k=9.0 h/{\rm Mpc}$ and out to $z=3$. Outside of this range the model reverts to halofit2012., Comment: 31 pages, 18 Figures. Matches version published in MNRAS. We release all power spectra data with the C-code NGenHalofit available at: https://CosmologyCode@bitbucket.org/ngenhalofitteam/ngenhalofitpublic.git

We explore the cosmological constraints from cosmic shear using a new way of modelling the non-linear matter correlation functions. The new formalism extends the method of Angulo & White (2010), which manipulates outputs of $N$-body simulations to represent the three-dimensional nonlinear mass distribution in different cosmological scenarios. We show that predictions from our approach for shear two-point correlations at $1$ to $300$ arcmin separations are accurate at the $\sim10$\% level, even for extreme changes in cosmology. For moderate changes, with target cosmologies similar to that preferred by analyses of recent Planck data, the accuracy is close to $\sim5$\%. We combine this approach with a MonteCarlo Markov Chain sampler to explore constraints on a $\Lambda$CDM model from the shear correlation functions measured in the Canada-France Hawaii Telescope Lensing Survey (CFHTLenS). We obtain constraints on the parameter combination $\sigma_8 (\Omega_m/0.27)^{0.6} = 0.801 \pm 0.028$. Combined with results from CMB data, we obtain marginalised constraints on $\sigma_8 = 0.81 \pm 0.01$ and $\Omega_m = 0.29 \pm 0.01$. These results are fully compatible with previous analyses, which supports the validity of our approach. We discuss the advantages of our method and the potential it offers, including a path to incorporate in detail the effects of baryons, among others effects, in future high-precision cosmological analyses., Comment: 12 pages, 8 figures

We simulate neutralino dark matter ($\chi$DM) haloes from their initial collapse, at $\sim$ earth mass, up to a few percent solar. Our results confirm that the density profiles of the first haloes are described by a $\sim r^{-1.5}$ power-law. As haloes grow in mass, their density profiles evolve significantly. In the central regions, they become shallower and reach on average $\sim r^{-1}$, the asymptotic form of an NFW profile. Using non-cosmological controlled simulations, we observe that temporal variations in the gravitational potential caused by major mergers lead to a shallowing of the inner profile. This transformation is more significant for shallower initial profiles and for a higher number of merging systems. Depending on the merger details, the resulting profiles can be shallower or steeper than NFW in their inner regions. Interestingly, mergers have a much weaker effect when the profile is given by a broken power-law with an inner slope of $-1$ (such as NFW or Hernquist profiles). This offers an explanation for the emergence of NFW-like profiles: after their initial collapse, $r^{-1.5}$ $\chi$DM haloes suffer copious major mergers, which progressively shallows the profile. Once an NFW-like profile is established, subsequent merging do not change the profile anymore. This suggests that halo profiles are not universal but rather a combination of (1) the physics of the formation of the microhaloes and (2) their early merger history -- both set by the properties of the dark matter particle -- as well as (3) the resilience of NFW-like profiles to perturbations., Comment: Matches version accepted by MNRAS. 16 pages, 14 figures

We present a Bayesian phase-space reconstruction of the cosmic large-scale matter density and velocity fields from the SDSS-III Baryon Oscillations Spectroscopic Survey Data Release 12 (BOSS DR12) CMASS galaxy clustering catalogue. We rely on a given $\Lambda$CDM cosmology, a mesh resolution in the range of 6-10 $h^{-1}$ Mpc, and a lognormal-Poisson model with a redshift dependent nonlinear bias. The bias parameters are derived from the data and a general renormalised perturbation theory approach. We use combined Gibbs and Hamiltonian sampling, implemented in the \textsc{argo} code, to iteratively reconstruct the dark matter density field and the coherent peculiar velocities of individual galaxies, correcting hereby for coherent redshift space distortions (RSD). Our tests relying on accurate $N$-body based mock galaxy catalogues, show unbiased real space power spectra of the nonlinear density field up to $k\sim0.2\, h$ Mpc$^{-1}$, and vanishing quadrupoles down to $r\sim20\,h^{-1}$ Mpc. We also demonstrate that the nonlinear cosmic web can be obtained from the tidal field tensor based on the Gaussian component of the reconstructed density field. We find that the reconstructed velocities have a statistical correlation coefficient compared to the true velocities of each individual lightcone mock galaxy of $r\sim0.68$ including about 10% of satellite galaxies with virial motions (about $r=0.75$ without satellites). The power spectra of the velocity divergence agree well with theoretical predictions up to $k\sim0.2\,h\,{\rm Mpc}^{-1}$. This work will be especially useful to improve, e.g. BAO reconstructions, kinematic Sunyaev-Zeldovich (kSZ), integrated Sachs-Wolfe (ISW) measurements, or environmental studies., Comment: 21 pages, 17 figures submitted to MNRAS

In a framework where galaxies mostly migrate on the colour-magnitude diagram from star-forming to quiescent, the green valley is considered a transitional galaxy stage. The details of the processes that drive galaxies from star-forming to passive systems still remain unknown. We developed a method that estimates empirically the star formation quenching times-scales of green valley galaxies, assuming an exponential decay model of the SFH and through a combination of narrow and broad bands from J-PLUS and GALEX. We correlate these quenching time-scales with the presence of bars. We find that the J-PLUS colours F0395-g and F0415-g are sensitive to different SFH, showing, a clear correlation with the Dn(4000) and H-delta,A spectral indices. We find that quenching time-scales obtained with our new approach are in agreement with those determined using spectral indices. We also find that galaxies with high bar probability tend to quench their star formation slowly. We conclude that: 1) J-PLUS filters can be used to measure quenching timescales in nearby green valley galaxies; and 2) the resulting star formation quenching time-scales are longer for barred green valley galaxies. Considering that the presence of a bar indicates that more violent processes (e.g., major mergers) are absent in host galaxies, we conclude that the presence of a bar can be used as a morphological signature for slow star formation quenching., Comment: 14 pages, 17 figures, accepted by Astronomy and Astrophysics

We investigate the large-scale clustering and gravitational interaction of baryons and dark matter (DM) over cosmic time using a set of collisionless N-body simulations. Both components, baryons and DM, are evolved from distinct primordial density and velocity power spectra as predicted by early-universe physics. We first demonstrate that such two-component simulations require an unconventional match between force and mass resolution (i.e. force softening on at least the mean particle separation scale). Otherwise, the growth on any scale is not correctly recovered because of a spurious coupling between the two species at the smallest scales. With these simulations, we then demonstrate how the primordial differences in the clustering of baryons and DM are progressively diminished over time. In particular, we explicitly show how the BAO signature is damped in the spatial distribution of baryons and imprinted in that of DM. This is a rapid process, yet it is still not fully completed at low redshifts. On large scales, the overall shape of the correlation function of baryons and DM differs by 2% at z = 9 and by 0.2% at z = 0. The differences in the amplitude of the BAO peak are approximately a factor of 5 larger: 10% at z = 9 and 1% at z = 0. These discrepancies are, however, smaller than effects expected to be introduced by galaxy formation physics in both the shape of the power spectrum and in the BAO peak, and are thus unlikely to be detected given the precision of the next generation of galaxy surveys. Hence, our results validate the standard practice of modelling the observed galaxy distribution using predictions for the total mass clustering in the Universe., 9 pages, 6 figures. Replaced with version published in MNRAS

We use the Millennium Simulation series to study the relation between the accretion history (MAH) and mass profile of cold dark matter halos. We find that the mean density within the scale radius, r_{-2} (where the halo density profile has isothermal slope), is directly proportional to the critical density of the Universe at the time when the main progenitor's virial mass equals the mass enclosed within r_{-2}. Scaled to these characteristic values of mass and density, the mean MAH, expressed in terms of the critical density of the Universe, M(\rho_{crit}(z)), resembles that of the enclosed density profile, M(), at z=0. Both follow closely the NFW profile, suggesting that the similarity of halo mass profiles originates from the mass-independence of halo MAHs. Support for this interpretation is provided by outlier halos whose accretion histories deviate from the NFW shape; their mass profiles show correlated deviations from NFW and are better approximated by Einasto profiles. Fitting both M() and M(\rho_{crit}) with either NFW or Einasto profiles yield concentration and shape parameters that are correlated, confirming and extending earlier work linking the concentration of a halo with its accretion history. These correlations also confirm that halo structure is insensitive to initial conditions: only halos whose accretion histories differ greatly from the NFW shape show noticeable deviations from NFW in their mass profiles. As a result, the NFW profile provides acceptable fits to hot dark matter halos, which do not form hierarchically, and for fluctuation power spectra other than CDM. Our findings, however, predict a subtle but systematic dependence of mass profile shape on accretion history which, if confirmed, would provide strong support for the link between accretion history and halo structure we propose here., Comment: 12 pages, 8 figures, MNRAS 432 1103L (2013)

We study the mass accretion histories (MAHs) and density profiles of dark matter halos using N-body simulations of self-similar gravitational clustering from scale-free power spectra, $P(k)\propto k^n$. We pay particular attention to the density profile curvature, which we characterize using the shape parameter, $\alpha$, of an Einasto profile. In agreement with previous findings our results suggest that, despite vast differences in their MAHs, the density profiles of virialized halos are remarkably alike. Nonetheless, clear departures from self-similarity are evident: for a given spectral index, $\alpha$ increases slightly but systematically with "peak height", $\nu\equiv\delta_{sc}/\sigma(M,z)$, regardless of mass or redshift. More importantly, however, the "$\alpha-\nu$" relation depends on $n$: the steeper the initial power spectrum, the more gradual the curvature of both the mean MAHs and mean density profiles. These results are consistent with previous findings connecting the shapes of halo mass profiles and MAHs and imply that dark matter halos are not structurally self-similar but, through the merger history, retain a memory of the linear density field from which they form., Comment: MNRAS letters, accepted

We construct a large, redshift-complete sample of distant galaxy clusters by correlating Sloan Digital Sky Survey Data Release 12 redshifts with clusters identified with the red-sequence Matched-filter Probabilistic Percolation (redMaPPer) algorithm. Our spectroscopic completeness is > 97 per cent for ≃ 7000 clusters within the redMaPPer selection limit, z=0.325, so that our cluster correlation functions aremuch more precise than earlierwork and not suppressed by uncertain photometric redshifts.We derive an accurate power-law mass-richness relation from the observed abundance with respect to the mass function from Millennium XXL (MXXL) simulations, adjusted to the Planck-weighted cosmology. The number density of clusters is found to decline by 20 per cent over the range 0.1 < z < 0.3, in good agreement with the evolution predicted by MXXL. Our projected 3D correlation function scales with richness, λ, rising from r = 14 h Mpc at λ ≃25 to r = 22 h Mpc at λ ≃ 60, with a gradient that matches MXXL when applying our mass-richness relation, whereas the observed amplitude of the correlation function at 〈z〉 = 0.24 exceeds the MXXL prediction by 20 per cent at the ≃2.5σ level. This tension cannot be blamed on spurious, randomly located clusters as this would reduce the correlation amplitude. Full consistency between the correlation function and the abundances is achievable for the pre-Planck values of σ = 0.9, Ω = 0.25 and h = 0.73, matching the improved distance ladder estimate of the Hubble constant., TJB is supported by IKERBASQUE, the Basque Foundation for Science. RL is supported by the Spanish Ministry of Economy and Competitiveness through research projects FIS2010-15492 and Consolider EPI CSD2010-00064, and the University of the Basque Country UPV/EHU under program UFI 11/55. PJ acknowledges financial support from the Basque Government grant BFI-2012-349. TJB, RL and PJ are also supported by the Basque Government grant for the GIC IT956-16 research group. REA acknowledges support from AYA2015-66211-C2-2. JMD acknowledges support of the projects AYA2015-64508-P (MINECO/FEDER, UE), AYA2012-39475-C02-01 and the consolider project CSD2010-00064 funded by the Ministerio de Economia y Competitividad. KU acknowledges partial support from the Ministry of Science and Technology of Taiwan (grants MOST 103-2112-M-001-030-MY3 and MOST 103-2112-M-001-003-MY3). Funding for SDSS-III has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation and the US Department of Energy Office of Science.

Subhalo abundance matching (SHAM) is a widely-used method to connect galaxies with dark matter structures in numerical simulations. SHAM predictions agree remarkably well with observations, yet they still lack strong theoretical support. We examine the performance, implementation, and assumptions of SHAM using the EAGLE project simulations. We find that $V_{\rm relax}$, the highest value of the circular velocity attained by a subhalo while it satisfies a relaxation criterion, is the subhalo property that correlates most strongly with galaxy stellar mass ($M_{\rm star}$). Using this parameter in SHAM, we retrieve the real-space clustering of EAGLE to within our statistical uncertainties on scales greater than $2$ Mpc for galaxies with $8.77, Comment: 20 pages, 14 figures

The abundance of galaxy clusters can constrain both the geometry and growth of structure in our Universe. However, this probe could be significantly complicated by recent claims of nonuniversality -- non-trivial dependences with respect to the cosmological model and redshift. In this work we analyse the dependance of the mass function on the way haloes are identified and establish if this can cause departures from universality. In order to explore this dependance, we use a set of different N-body cosmological simulations (Le SBARBINE simulations), with the latest cosmological parameters from the Planck collaboration; this first suite of simulations is followed by a lower resolution set, carried out with different cosmological parameters. We identify dark matter haloes using a Spherical Overdensity algorithm with varying overdensity thresholds (virial, 2000rho_c, 1000rho_c, 500rho_c, 200rho_c and 200rho_b) at all redshifts. We notice that, when expressed in term of the rescaled variable nu, the mass functionfor virial haloes is a nearly universal as a function of redshift and cosmology, while this is clearly not the case for the other overdensities we considered. We provide fitting functions for the halo mass function parameters as a function of overdensity, that allow to predict, to within a few percent accuracy, the halo mass function for a wide range of halo definitions, redshifts and cosmological models. We then show how the departures from universality associated with other halo definitions can be derived by combining the universality of the virial definition with the expected shape of the density profile of halos., 22 pages, 21 figures, accepted for publication in MNRAS - some changes with respect to the previous version and a new appendix

We use a suite of cosmological simulations to study the mass-concentration-redshift relation, $c({\rm M},z)$, of dark matter halos. Our simulations include standard $\Lambda$-cold dark matter (CDM) models, and additional runs with truncated power spectra, consistent with a thermal warm dark matter (WDM) scenario. We find that the mass profiles of CDM and WDM halos are self-similar and well approximated by the Einasto profile. The $c({\rm M},z)$ relation of CDM halos is monotonic: concentrations decrease with increasing virial mass at fixed redshift, and decrease with increasing redshift at fixed mass. The mass accretion histories (MAHs) of CDM halos are also scale-free, and can be used to infer concentrations directly. These results do not apply to WDM halos: their MAHs are not scale-free because of the characteristic scale imposed by the power-spectrum suppression. Further, the WDM $c({\rm M},z)$ relation is non-monotonic: concentrations peak at a mass scale dictated by the truncation scale, and decrease at higher and lower masses. We show that the assembly history of a halo can still be used to infer its concentration, provided that the total mass of its progenitors is considered (the "collapsed mass history"; CMH), rather than just that of its main ancestor. This exploits the scale-free nature of CMHs to derive a simple scaling that reproduces the mass-concentration-redshift relation of both CDM and WDM halos over a vast range of halo masses and redshifts. Our model therefore provides a robust account of the mass, redshift, cosmology and power spectrum dependence of dark matter halo concentrations., Comment: 21 pages, 16 figures

In this paper we investigate the strong lensing statistics in galaxy clusters. We extract dark matter haloes from the Millennium-XXL simulation, compute their Einstein radius distribution, and find a very good agreement with Monte Carlo predictions produced with the MOKA code. The distribution of the Einstein radii is well described by a log-normal distribution, with a considerable fraction of the largest systems boosted by different projection effects. We discuss the importance of substructures and triaxiality in shaping the size of the critical lines for cluster size haloes. We then model and interpret the different deviations, accounting for the presence of a Bright Central Galaxy (BCG) and two different stellar mass density profiles. We present scaling relations between weak lensing quantities and the size of the Einstein radii. Finally we discuss how sensible is the distribution of the Einstein radii on the cosmological parameters {\Omega}_M-{\sigma}_8 finding that cosmologies with higher {\Omega}_M and {\sigma}_8 possess a large sample of strong lensing clusters. The Einstein radius distribution may help distinguish Planck13 and WMAP7 cosmology at 3{\sigma}., Comment: replaced to match the accepted version for publication on MNRAS

Most of the upcoming cosmological information will come from analyzing the clustering of the Large Scale Structures (LSS) of the universe through LSS or CMB observations. It is therefore essential to be able to understand their behavior with exquisite precision. The Effective Field Theory of Large Scale Structures (EFTofLSS) provides a consistent framework to make predictions for LSS observables in the mildly non-linear regime. In this paper we focus on biased tracers. We argue that in calculations at a given order in the dark matter perturbations, highly biased tracers will underperform because of their larger higher derivative biases. A natural prediction of the EFTofLSS is therefore that by simply adding higher derivative biases, all tracers should perform comparably well. We implement this prediction for the halo-halo and the halo-matter power spectra at one loop, and the halo-halo-halo, halo-halo-matter, and halo-matter-matter bispectra at tree-level, and compare with simulations. We find good agreement with the prediction: for all tracers, we are able to match the bispectra up to $k\simeq0.17\,h/$Mpc at $z=0$ and the power spectra to a higher wavenumber., Comment: 39 pages, 8 figures

We provide constraints on the accuracy with which the neutrino mass fraction, $f_{\nu}$, can be estimated when exploiting measurements of redshift-space distortions, describing in particular how the error on neutrino mass depends on three fundamental parameters of a characteristic galaxy redshift survey: density, halo bias and volume. In doing this, we make use of a series of dark matter halo catalogues extracted from the BASICC simulation. The mock data are analysed via a Markov Chain Monte Carlo likelihood analysis. We find a fitting function that well describes the dependence of the error on bias, density and volume, showing a decrease in the error as the bias and volume increase, and a decrease with density down to an almost constant value for high density values. This fitting formula allows us to produce forecasts on the precision achievable with future surveys on measurements of the neutrino mass fraction. For example, a Euclid-like spectroscopic survey should be able to measure the neutrino mass fraction with an accuracy of $\delta f_{\nu} \approx 6.7\times10^{-4}$, using redshift-space clustering once all the other cosmological parameters are kept fixed to the $\Lambda$CDM case., Comment: 14 pages, 19 figures