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

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 < z < 1.5 and over scales $10^{-2} < k/(h Mpc^{-1} ) < 5$ for 800 points in an 8-dimensional cosmological parameter space. For an efficient interpolation of the results, we build a Gaussian emulator which we test against other widely-used methods. Over the whole range of scales considered, we expect our predictions to be accurate at the 2% level for parameters in the minimal ${\Lambda}$CDM model and to 3% when extended to dynamical dark energy and massive neutrinos. We make our emulator publicly available under http://www.dipc.org/bacco

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.

We present and test a framework that models the three-dimensional distribution of mass in the Universe as a function of cosmological and astrophysical parameters. Our approach combines two different techniques: a rescaling algorithm that modifies the cosmology of gravity-only N-body simulations, and a baryonification algorithm which mimics the effects of astrophysical processes induced by baryons, such as star formation and AGN feedback. We show how this approach can accurately reproduce the effects of baryons on the matter power spectrum of various state-of-the-art hydro-dynamical simulations (EAGLE, Illustris, Illustris-TNG, Horizon-AGN, and OWLS,Cosmo-OWLS and BAHAMAS), to percent level from very large down to small, highly nonlinear scales, k= 5 h/Mpc, and from z=0 up to z=2. We highlight that, thanks to the heavy optimisation of the algorithms, we can obtain these predictions for arbitrary baryonic models and cosmology (including massive neutrinos and dynamical dark energy models) with an almost negligible CPU cost. Therefore, this approach is efficient enough for cosmological data analyses. With these tools in hand we explore the degeneracies between cosmological and astrophysical parameters in the nonlinear mass power spectrum. Our findings suggest that after marginalising over baryonic physics, cosmological constraints inferred from weak gravitational lensing should be moderately degraded.

Providing accurate predictions for the spatial distribution of matter and luminous tracers in the presence of massive neutrinos is an important task, given the imminent arrival of highly accurate large-scale structure observations. In this work, we address this challenge by extending cosmology-rescaling algorithms to massive neutrino cosmologies. In this way, a $\Lambda$CDM simulation can be modified to provide nonlinear structure formation predictions in the presence a hot component of arbitrary mass, and, if desired, to include non-gravitational modifications to the clustering of matter on large scales. We test the accuracy of the method by comparing its predictions to a suite of simulations carried out explicitly including a neutrino component in its evolution equations. We find that, for neutrino masses in the range $M_\nu \in [0.06, 0.3] ~ \mathrm{eV}$ the matter power spectrum is recovered to better than $1\%$ on all scales $k<2~h~\mathrm{Mpc}^{-1}$. Similarly, the halo mass function is predicted at a few percent level over the range $M_{\rm halo} \in [10^{12}, 10^{15}] ~ h^{-1} ~ \mathrm{M}_{\odot}$, and so do also the multipoles of the galaxy 2-point correlation function in redshift space over $r \in [0.1, 200] ~ h^{-1} ~ \mathrm{Mpc}$. We provide parametric forms for the necessary transformations, as a function of $\Omega_{\rm m}$ and $\Omega_{\nu}$ for various target redshifts.

Recent work revealed that both the helium variation within globular clusters (GCs) and the relative numbers of first and second-generation stars (1G, 2G) depend on the mass of the host cluster. Precise determination of the internal helium variations and of the fraction of 1G stars are crucial constraints to the formation scenarios of multiple populations (MPs). We exploit multi-band Hubble Space Telescope photometry to investigate MPs in NGC 2419, which is one of the most-massive and distant GCs of the Galaxy, almost isolated from its tidal influence. We find that the 1G hosts the ~37% of the analyzed stars, and identified three populations of 2G stars, namely 2GA, 2GB, and 2GC, which comprise the ~20%, ~31% and ~12% of stars, respectively. We compare the observed colors of these four populations with the colors derived from appropriate synthetic spectra to infer the relative helium abundances. We find that 2GA, 2GB, and 2GC stars are enhanced in helium mass fraction by deltaY ~0.01, 0.06, and 0.19 with respectto 1G stars that have primordial helium (Y=0.246). The high He enrichment of 2GC stars is hardly reconcilable with most of the current scenarios for MPs. Furthermore, the relatively larger fraction of 1G stars (~37%) compared to other massive GCs is noticeable. By exploiting literature results, we find that the fractions of 1G stars of GCs with large perigalactic distance are typically higher than in the other GCs with similar masses. This suggests that NGC 2419, similarly to other distant GCs, lost a lower fraction of 1G stars.

Galaxy redshift surveys are one of the pillars of the current standard cosmological model and remain a key tool in the experimental effort to understand the origin of cosmic acceleration. To this end, the next generation of surveys aim at achieving sub-percent precision in the measurement of the equation of state of dark energy $w(z)$ and the growth rate of structure $f(z)$. This however requires comparable control over systematic errors, stressing the need for improved modelling methods. In this contribution we review at the introductory level some highlights of the work done in this direction by the {\it Darklight} project. Supported by an ERC Advanced Grant, {\it Darklight} developed novel techniques for clustering analysis, which were tested through numerical simulations before being finally applied to galaxy data as in particular those of the recently completed VIPERS redshift survey. We focus in particular on: (a) advances on estimating the growth rate of structure from redshift-space distortions; (b) parameter estimation through global Bayesian reconstruction of the density field from survey data; (c) impact of massive neutrinos on large-scale structure measurements. Overall, {\it Darklight} has contributed to paving the way for forthcoming high-precision experiments, such as {\it Euclid}, the next ESA cosmological mission.

The clustering ratio is defined as the ratio between the correlation function and the variance of the smoothed overdensity field. In LCDM cosmologies not accounting for massive neutrinos, it has already been proved to be independent from bias and redshift space distortions on a range of linear scales. It therefore allows for a direct comparison of measurements (from galaxies in redshift space) to predictions (for matter in real space). In this paper we first extend the applicability of such properties of the clustering ratio to cosmologies that include massive neutrinos, by performing tests against simulated data. We then investigate the constraining power of the clustering ratio when cosmological parameters such as the total neutrino mass and the equation of state of dark energy are left free. We analyse the joint posterior distribution of the parameters that must satisfy, at the same time, the measurements of the galaxy clustering ratio in the SDSS DR12, and the angular power spectrum of temperature and polarization anisotropies of the CMB measured by the Planck satellite. We find the clustering ratio to be very sensitive to the CDM density parameter, but not very much so to the total neutrino mass. Lastly, we forecast the constraining power the clustering ratio will achieve with forthcoming surveys, predicting the amplitude of its errors in a Euclid-like galaxy survey. In this case, we find it is expected to improve the constraint at 95% level on the CDM density by 40% and on the total neutrino mass by 14%.

The set-up of the initial conditions in cosmological N-body simulations is usually implemented by rescaling the desired low-redshift linear power spectrum to the required starting redshift consistently with the Newtonian evolution of the simulation. The implementation of this practical solution requires more care in the context of massive neutrino cosmologies, mainly because of the non-trivial scale-dependence of the linear growth that characterises these models. In this work we consider a simple two-fluid, Newtonian approximation for cold dark matter and massive neutrinos perturbations that can reproduce the cold matter linear evolution predicted by Boltzmann codes such as CAMB or CLASS with a 0.1% accuracy or below for all redshift relevant to nonlinear structure formation. We use this description, in the first place, to quantify the systematic errors induced by several approximations often assumed in numerical simulations, including the typical set-up of the initial conditions for massive neutrino cosmologies adopted in previous works. We then take advantage of the flexibility of this approach to rescale the late-time linear power spectra to the simulation initial redshift, in order to be as consistent as possible with the dynamics of the N-body code and the approximations it assumes. We implement our method in a public code providing the initial displacements and velocities for cold dark matter and neutrino particles that will allow accurate, i.e. one-percent level, numerical simulations for this cosmological scenario.

We derive the luminosity function and redshift distribution of short Gamma Ray Bursts (SGRBs) using (i) all the available observer-frame constraints (i.e. peak flux, fluence, peak energy and duration distributions) of the large population of Fermi SGRBs and (ii) the rest-frame properties of a complete sample of Swift SGRBs. We show that a steep $\phi(L)\propto L^{-a}$ with a>2.0 is excluded if the full set of constraints is considered. We implement a Monte Carlo Markov Chain method to derive the $\phi(L)$ and $\psi(z)$ functions assuming intrinsic Ep-Liso and Ep-Eiso correlations or independent distributions of intrinsic peak energy, luminosity and duration. To make our results independent from assumptions on the progenitor (NS-NS binary mergers or other channels) and from uncertainties on the star formation history, we assume a parametric form for the redshift distribution of SGRBs. We find that a relatively flat luminosity function with slope ~0.5 below a characteristic break luminosity ~3$\times10^{52}$ erg/s and a redshift distribution of SGRBs peaking at z~1.5-2 satisfy all our constraints. These results hold also if no Ep-Liso and Ep-Eiso correlations are assumed. We estimate that, within ~200 Mpc (i.e. the design aLIGO range for the detection of GW produced by NS-NS merger events), 0.007-0.03 SGRBs yr$^{-1}$ should be detectable as gamma-ray events. Assuming current estimates of NS-NS merger rates and that all NS-NS mergers lead to a SGRB event, we derive a conservative estimate of the average opening angle of SGRBs: $\theta_{jet}$~3-6 deg. Our luminosity function implies an average luminosity L~1.5$\times 10^{52}$ erg/s, nearly two orders of magnitude higher than previous findings, which greatly enhances the chance of observing SGRB "orphan" afterglows. Efforts should go in the direction of finding and identifying such orphan afterglows as counterparts of GW events.