Search

We compute different sets of stellar evolutionary tracks in order to quantify the energy, mass, and metals yielded by massive main-sequence and post-main-sequence winds. Our aim is to investigate the impact of binary systems and of a metallicity-dependent distribution of initial rotational velocities on the feedback by stellar winds. We find significant changes compared to the commonly used non-rotating, single-star scenario. The largest differences are noticeable at low metallicity, where the mechanical-energy budget is substantially increased. So as to establish the maximal (i.e. obtained by neglecting dissipation in the near circumstellar environment) influence of winds on the early stages of galaxy formation, we use our new feedback estimates to simulate the formation and evolution of a sub-$L_*$ galaxy at redshift 3 (hosted by a dark-matter halo with a mass of $1.8\times 10^{11}$ M$_\odot$) and compare the outcome with simulations in which only supernova (SN) feedback is considered. Accounting for the continuous energy injection by winds reduces the total stellar mass, the metal content, and the burstiness of the star-formation rate as well as of the outflowing gas mass. However, our numerical experiment suggests that the enhanced mechanical feedback from the winds of rotating and binary stars has a limited impact on the most relevant galactic properties compared to the non-rotating single-star scenario. Eventually, we look at the relative abundance between the metals entrained in winds and those ejected by SNe and find that it stays nearly constant within the simulated galaxy and its surrounding halo., Comment: 20 pages, 20 figures, updated to match the version published in MNRAS

The detection of an unexpectedly large number of faint candidate quasars (QSO) at $z \gtrsim 4$ has motivated Madau and Haardt to investigate a reionization scenario in which active galactic nuclei dominate the photon budget at all times. Their analytical study reveals that this picture is in agreement with the evolution of the HI volume fraction in the intergalactic medium (IGM) and the optical depth to Thomson scattering of the cosmic microwave background. We employ a suite of hydrodynamical simulations post-processed with a radiative transfer code to further investigate the properties of such a reionization history. In particular, we generate synthetic absorption-line spectra for HI and HeII that we compare with observational data. Although we confirm the analytical results mentioned above, we also find that the evolution of the IGM temperature and of the HeII optical depth are at odds with current observational constraints. Nevertheless, the QSO-dominated scenario presents the attractive feature that it naturally generates an inhomogeneous IGM and thus produces an extended tail at high values in the distribution function of the HI effective optical depth, in agreement with recent observations. We find tentative evidence that considering some QSO emission at redshift $z \gtrsim 5$ could reconcile such observations with numerical simulations that have so far failed to reproduce this feature. Finally, we show that collecting more HeII-absorption spectra at $z \gtrsim 3$ and studying their distinctive characteristics will be key to precisely constraining the QSO contribution to reionization., 15 pages, 13 figures, published on MNRAS

The streaming model describes the mapping between real and redshift space for 2-point clustering statistics. Its key element is the probability density function (PDF) of line-of-sight pairwise peculiar velocities. Following a kinetic-theory approach, we derive the fundamental equations of the streaming model for ordered and unordered pairs. In the first case, we recover the classic equation while we demonstrate that modifications are necessary for unordered pairs. We then discuss several statistical properties of the pairwise velocities for DM particles and haloes by using a suite of high-resolution $N$-body simulations. We test the often used Gaussian ansatz for the PDF of pairwise velocities and discuss its limitations. Finally, we introduce a mixture of Gaussians which is known in statistics as the generalised hyperbolic distribution and show that it provides an accurate fit to the PDF. Once inserted in the streaming equation, the fit yields an excellent description of redshift-space correlations at all scales that vastly outperforms the Gaussian and exponential approximations. Using a principal-component analysis, we reduce the complexity of our model for large redshift-space separations. Our results increase the robustness of studies of anisotropic galaxy clustering and are useful for extending them towards smaller scales in order to test theories of gravity and interacting dark-energy models., Comment: 22 pages, 20 figures, accepted for publication in MNRAS

We use hydrodynamic simulations post-processed with the radiative-transfer code RADAMESH to assess recent claims that the low HeII opacity observed in z>3 quasar spectra may be incompatible with models of HeII reionization driven by the observed population of active galactic nuclei (AGNs). In particular, building upon our previous work, we consider an early population of sources and start the radiative-transfer calculation at redshifts z>=5. Our model faithfully reproduces the emissivity of optically selected AGNs as inferred from measurements of their luminosity function. We find that HeII reionization is very extended in redshift ({\Delta} z>=2) and highly spatially inhomogeneous. In fact, mock spectra extracted from the simulations show a large variability in the evolution of the HeII effective optical depth within chunks of size {\Delta} z=0.04. Regions with low opacity ({\tau}^{eff}_{HeII}, Comment: 12 pages, 8 figures. Matches version accepted for publication in MNRAS

We use two cosmological simulations of structure formation in the LambdaCDM scenario to study the evolutionary histories of dark-matter haloes and to characterize the Lagrangian regions from which they form. We focus on haloes identified at redshift z_id=0 and show that the classic ellipsoidal collapse model systematically overestimates their collapse times. If one imposes that halo collapse takes place at z_id, this model requires starting from a significantly lower linear density contrast than what is measured in the simulations at the locations of halo formation. We attempt to explain this discrepancy by testing two key assumptions of the model. First, we show that the tides felt by collapsing haloes due to the surrounding large-scale structure evolve non-linearly. Although this effect becomes increasingly important for low-mass haloes, accounting for it in the ellipsoidal collapse model only marginally improves the agreement with N-body simulations. Second, we track the time evolution of the physical volume occupied by forming haloes and show that, after turnaround, it generally stabilizes at a well-defined redshift, z_c>z_id, contrary to the basic assumption of extended Press-Schechter theory based on excursion sets. We discuss the implications of this result for understanding the origin of the mass-dependence and scatter in the linear threshold for halo formation. Finally, we show that, when tuned for collapse at z_c, a modified version of the ellipsoidal collapse model that also accounts for the triaxial nature of protohaloes predicts their linear density contrast in an unbiased way.

We present a new approach for modelling galaxy/halo bias that utilizes the full non-linear information contained in the moments of the matter density field, which we derive using a set of numerical simulations. Although our method is general, we perform a case study based on the local Eulerian bias scheme truncated to second-order. Using 200 N-body simulations covering a total comoving volume of 675 (Gpc/h)^3, we measure several 2- and 3-point statistics of the halo distribution to unprecedented accuracy. We use the bias model to fit the halo-halo power spectrum, the halo-matter cross spectrum and the corresponding three bispectra for wavenumbers in the range 0.04 < k < 0.12 h Mpc^-1. We find the constraints on the bias parameters obtained using the full non-linear information differ significantly from those derived using standard perturbation theory at leading order. Hence, neglecting the full non-linear information leads to biased results for this particular scale range. We also test the validity of the second-order Eulerian local biasing scheme by comparing the parameter constraints derived from different statistics. Analysis of the halo-matter cross-correlation coefficients defined for the 2- and 3-point statistics reveals further inconsistencies contained in the second-order Eulerian bias scheme, suggesting it is too simple a model to describe halo bias with high accuracy., 23 pages, 11 Figures, minor corrections ; accepted for publication in MNRAS

We use a set of AMR hydrodynamic simulations post-processed with the radiative-transfer code RADAMESH to study how inhomogeneous HeII reionization affects the intergalactic medium (IGM). We propagate radiation from active galactic nuclei (AGNs) considering two scenarios for the time evolution of the ionizing sources. We find that HeII reionization takes place in a very inhomogeneous fashion, through the production of well separated bubbles of the ionized phase that subsequently percolate. Overall, the reionization process is extended in time and lasts for a redshift interval Delta z>1. At fixed gas density, the temperature distribution is bimodal during the early phases of HeII reionization and cannot be described by a simple effective equation of state. When HeII reionization is complete, the IGM is characterized by a polytropic equation of state with index gamma~1.20. This relation is appreciably flatter than at the onset of the reionization process (gamma=1.56) and also presents a much wider dispersion around the mean. We extract HI and HeII Ly-alpha absorption spectra from the simulations and fit Voigt profiles to them. We find that the regions where helium is doubly ionized are characterized by different probability density functions of the curvature and of the Doppler b parameters of the HI Ly-alpha forest as a consequence of the bimodal temperature distribution during the early phases of HeII reionization. The column-density ratio in HeII and HI shows a strong spatial variability. Its probability density function rapidly evolves with time reflecting the increasing volume fraction in which ionizing radiation is harder due to the AGN contribution. Finally we show that the number density of the flux-transmission windows per unit redshift and the mean size of the dark gaps in the HeII spectra have the potential to distinguish between different reionization scenarios. (abridged), 24 pages, 30 figures. Matches version accepted for publication in MNRAS

The scale-dependent galaxy bias generated by primordial non-Gaussianity (PNG) can be used to detect and constrain deviations from standard single-field inflation. The strongest signal is expected in the local model for PNG, where the amplitude of non-Gaussianity can be expressed by a set of parameters (fnl, gnl, ...). Current observational constraints from galaxy clustering on fnl and gnl assume that the others PNG parameters are vanishing. Using two sets of cosmological N-body simulations where both fnl and gnl are non-zero, we show that this strong assumption generally leads to biased estimates and spurious redshift dependencies of the parameters. Additionally, if the signs of fnl and gnl are opposite, the amplitude of the scale-dependent bias is reduced, possibly leading to a false null detection. Finally we show that model selection techniques like the Bayesian evidence can (and should) be used to determine if more than one PNG parameter is required by the data.

We study the constraining power on primordial non-Gaussianity of future surveys of the large-scale structure of the Universe for both near-term surveys (such as the Dark Energy Survey - DES) as well as longer term projects such as Euclid and WFIRST. Specifically we perform a Fisher matrix analysis forecast for such surveys, using DES-like and Euclid-like configurations as examples, and take account of any expected photometric and spectroscopic data. We focus on two-point statistics and we consider three observables: the 3D galaxy power spectrum in redshift space, the angular galaxy power spectrum, and the projected weak-lensing shear power spectrum. We study the effects of adding a few extra parameters to the basic LCDM set. We include the two standard parameters to model the current value for the dark energy equation of state and its time derivative, w_0, w_a, and we account for the possibility of primordial non-Gaussianity of the local, equilateral and orthogonal types, of parameter fNL and, optionally, of spectral index n_fNL. We present forecasted constraints on these parameters using the different observational probes. We show that accounting for models that include primordial non-Gaussianity does not degrade the constraint on the standard LCDM set nor on the dark-energy equation of state. By combining the weak lensing data and the information on projected galaxy clustering, consistently including all two-point functions and their covariance, we find forecasted marginalised errors sigma (fNL) ~ 3, sigma (n_fNL) ~ 0.12 from a Euclid-like survey for the local shape of primordial non-Gaussianity, while the orthogonal and equilateral constraints are weakened for the galaxy clustering case, due to the weaker scale-dependence of the bias. In the lensing case, the constraints remain instead similar in all configurations.

Starting in late 2013, the eROSITA telescope will survey the X-ray sky with unprecedented sensitivity. Assuming a detection limit of 50 photons in the (0.5-2.0) keV energy band with a typical exposure time of 1.6 ks, we predict that eROSITA will detect ~ 9.3 X 10^4 clusters of galaxies more massive than 5 X 10^13 Msun/h, with the currently planned all-sky survey. Their median redshift will be z ~ 0.35. We perform a Fisher-matrix analysis to forecast the constraining power of eROSITA on the LambdaCDM cosmology and, simultaneously, on the X-ray scaling relations for galaxy clusters. Special attention is devoted to the possibility of detecting primordial non-Gaussianity. We consider two experimental probes: the number counts and the angular clustering of a photon-count limited sample of clusters. We discuss how the cluster sample should be split to optimize the analysis and we show that redshift information of the individual clusters is vital to break the strong degeneracies among the model parameters. For example, performing a "tomographic" analysis based on photometric-redshift estimates and combining 1- and 2-point statistics will give marginal 1-sigma errors of Delta(sigma_8) ~ 0.036 and Delta(Omega_m) ~ 0.012 without priors, and improve the current estimates on the slope of the luminosity-mass relation by a factor of 3. Regarding primordial non-Gaussianity, eROSITA clusters alone will give Delta(f_NL) = 9, 36, 144 for the local, orthogonal and equilateral model, respectively. Measuring redshifts with spectroscopic accuracy would further tighten the constraints by nearly 40 per cent (barring f_NL which displays smaller improvements). Finally, combining eROSITA data with the analysis of temperature anisotropies in the cosmic microwave background by the Planck satellite should give sensational constraints on both the cosmology and the properties of the intracluster medium.

We use high-resolution N-body simulations to investigate the Lagrangian bias of cold dark matter haloes within the Λ cold dark matter cosmology. Our analysis focuses on ‘protohaloes’ that we identify in the simulation initial conditions with the subsets of particles belonging to individual redshift-zero haloes. We then calculate the number density and velocity divergence fields of protohaloes and estimate their autospectral densities. We also measure the corresponding cross-spectral densities with the linear matter distribution. We use our results to test a Lagrangian bias model presented by Desjacques & Sheth which is based on the assumption that haloes form out of local density maxima of a specific height. Our comparison validates the predicted functional form for the scale dependence of the bias for both the density and the velocity fields. We also show that the bias coefficients are accurately predicted for the velocity divergence. In contrast, the theoretical values for the density bias parameters do not accurately match the numerical results as a function of halo mass. This is likely due to the simplistic assumptions that relate virialized haloes to density peaks of a given height in the model. We also detect appreciable stochasticity for the Lagrangian density bias, even on very large scales. These are not included in the model at leading order but correspond to higher order corrections.

We study the relation between the halo and matter density fields -- commonly termed bias -- in the LCDM framework. In particular, we examine the local model of biasing at quadratic order in the matter density. This model is characterized by parameters b_1 and b_2. Using an ensemble of N-body simulations, we apply several statistical methods to estimate the parameters. We measure halo and matter fluctuations smoothed on various scales and find that the parameters vary with smoothing scale. We argue that, for real-space measurements, owing to the mixing of wavemodes, no scale can be found for which the parameters are independent of smoothing. However, this is not the case in Fourier space. We measure halo power spectra and construct estimates for an effective large-scale bias. We measure the configuration dependence of the halo bispectra B_hhh and reduced bispectra Q_hhh for very large-scale k-space triangles. From this we constrain b_1 and b_2. Using the lowest-order perturbation theory, we find that for B_hhh the best-fit parameters are in reasonable agreement with one another as the triangle scale is varied, but that the fits become poor as smaller scales are included. The same is true for Q_hhh. The best-fit parameters depend on the discreteness correction. This led us to consider halo-mass cross-bispectra. The results from these statistics support our earlier findings. We develop a test to explore the importance of missing higher-order terms in the models. We prove that low-order expansions are not able to correctly model the data, even on scales k_1~0.04 h/Mpc. If robust inferences are to be drawn from galaxy surveys, then accurate models for the full nonlinear matter bispectrum and trispectrum will be essential.

We address the issue of the cosmological bias between matter and galaxy distributions, looking at dark-matter haloes as a first step to characterize galaxy clustering. Starting from the linear density field at high redshift, we follow the centre of mass trajectory of the material that will form each halo at late times (proto-halo). We adopt a fluid-like description for the evolution of perturbations in the proto-halo distribution, which is coupled to the matter density field via gravity. We present analytical solutions for the density and velocity fields, in the context of renormalized perturbation theory. We start from the linear solution, then compute one-loop corrections for the propagator and the power spectrum. Finally we analytically resum the propagator and we use a suitable extension of the time-renormalization-group method (Pietroni 2008) to resum the power spectrum. For halo masses M

We test third-order standard perturbation theory (SPT) as an approximation to non-linear cosmological structure formation. A novel approach is used to numerically calculate the three-dimensional dark matter density field using SPT from the initial conditions of two high-resolution cosmological simulations. The calculated density field is compared to the non-linear dark matter field of the simulations both point-by-point and statistically. For smoothing scales above 8 Mpc/h it shows a good agreement up to redshift 0. We present a simple fitting formula to relate the linear and non-linear density contrast that accurately recovers the non-linear time evolution for 0

We use two cosmological simulations of structure formation to study the conditions under which dark matter haloes emerge from the linear density field. Our analysis focuses on matching sites of halo collapse to local density maxima, or "peaks", in the initial conditions of the simulations and provides a crucial test of the central ansatz of the peaks formalism. By identifying peaks on a variety of smoothed, linearly extrapolated density fields we demonstrate that as many as ~70% of well-resolved dark matter haloes form preferentially near peaks whose characteristic masses are similar to that of the halo, with more massive haloes showing a stronger tendency to reside near peaks initially. We identify a small but significant fraction of haloes that appear to evolve from peaks of substantially lower mass than that of the halo itself. We refer to these as "peakless haloes" for convenience. By contrasting directly the properties of these objects with the bulk of the proto-halo population we find two clear differences: 1) their initial shapes are significantly flatter and more elongated than the predominantly triaxial majority, and 2) they are, on average, more strongly compressed by tidal forces associated with their surrounding large scale structure. Using the two-point correlation function we show that peakless haloes tend to emerge from highly clustered regions of the initial density field implying that, at fixed mass, the accretion geometry and mass accretion histories of haloes in highly clustered environments differ significantly from those in the field. This may have important implications for understanding the origin of the halo assembly bias, of galaxy properties in dense environments and how environment affects the morphological transformation of galaxies near groups and rich galaxy clusters.

We investigate how the shape of the galaxy two-point correlation function as measured in the zCOSMOS survey depends on local environment, quantified in terms of the density contrast on scales of 5 Mpc/h. We show that the flat shape previously observed at redshifts between z=0.6 and z=1 can be explained by this volume being simply 10% over-abundant in high-density environments, with respect to a Universal density probability distribution function. When galaxies corresponding to the top 10% tail of the distribution are excluded, the measured w_p(r_p) steepens and becomes indistinguishable from LCDM predictions on all scales. This is the same effect recognised by Abbas & Sheth in the SDSS data at z~0 and explained as a natural consequence of halo-environment correlations in a hierarchical scenario. Galaxies living in high-density regions trace dark matter halos with typically higher masses, which are more correlated. If the density probability distribution function of the sample is particularly rich in high-density regions because of the variance introduced by its finite size, this produces a distorted two-point correlation function. We argue that this is the dominant effect responsible for the observed "peculiar" clustering in the COSMOS field.

We perform a series of high-resolution N-body simulations of cosmological structure formation starting from Gaussian and non-Gaussian initial conditions. We adopt the best-fitting cosmological parameters of WMAP (3rd- and 5th-year) and we consider non-Gaussianity of the local type parameterised by 8 different values of the non-linearity parameter F_NL. Building upon previous work based on the Gaussian case, we show that, expressed in terms of suitable variables, the mass function of friends-of-friends haloes is approximately universal (i.e. independent of redshift, cosmology, and matter transfer function) to good precision (nearly 10 per cent) also in non-Gaussian scenarios. We provide fitting formulae for the high-mass end (M>10^13 M_sol/h) of the universal mass function in terms of F_NL, and we also present a non-universal fit in terms of both F_NL and z to be used for applications requiring higher accuracy. In the Gaussian case, we extend our fit to a wider range of halo masses (M>2.4 x 10^10 M_sol/h) and we also provide a consistent fit of the linear halo bias. We show that, for realistic values of F_NL, the matter power-spectrum in non-Gaussian cosmologies departs from the Gaussian one by up to 2 per cent on the scales where the baryonic- oscillation features are imprinted on the 2-point statistics. We confirm the strong k-dependence of the halo bias on large scales (k

We measure the clustering of a sample of photometrically selected luminous red galaxies around a low redshift (0.2~ 10^{12} h^{-1} M_sun, Eddington ratios from 0.01 to 1 and lifetimes less than 10^{7} yr. Using simple models of halo occupation, these correspond to a number density of quasar hosts greater than 10^{-3} h^{3} Mpc^{-3} and stellar masses less than 10^{11} h^{-1} M_sun. The small-scale clustering signal can be interpreted with the aid of our mock LRG catalogs, and depends on the manner in which quasars inhabit halos. We find that our small scale measurements are inconsistent with quasar positions being randomly subsampled from halo centers above a mass threshold, requiring a satellite fraction > 25 per cent., Comment: 16 pages, 11 figures, submitted to MNRAS

thly notices of the Royal Astronomical Society, 346 (1), ISSN:0035-8711, ISSN:1365-2966

We study correlations amongst tidal fields originated by the large-scale distribution of matter in the Universe. The two-point tidal correlation is described as a rank-4 tensor, whose elements can be written in terms of four fundamental scalar functions ranging, with respect the spatial separation, from purely transversal to purely longitudinal correlations. Tidal fields, both on galaxy and cluster scales, reveal to be correlated over distances larger than the mass-density correlation lenght, though traceless tidal fields show anti-correlation between diagonal terms along orthogonal directions. The cross-correlation between mass and tidal field is also analyzed. These results are relevant for galaxy formation and the interpretation of large-scale weak lensing phenomena., Comment: 6 pgs, 4 figs, using amssym.sty and mn.sty - MNRAS (in press)

We study the evolution of the mass autocorrelation function by describing the growth of density fluctuations through the Zel'dovich approximation. The results are directly compared with the predictions of the scaling hypothesis for clustering evolution extracted from numerical simulations (Hamilton et al. 1991), as implemented by Jain, Mo & White (1995). We find very good agreement between the correlations on mildly non-linear scales and on completely linear scales. In between these regimes, we note that the density fields evolved through the Zel'dovich approximation show more non-linear features than predicted by the scaling ansatz which is, however, forced to match the linear evolution on scales larger than the simulation box. In any case, the scaling ansatz by Baugh & Gaztanaga (1996), calibrated against large box simulations agrees better with ZA predictions on large scales, keeping good accuracy also on intermediate scales. We show that mode-coupling is able to move the first zero crossing of xi(r) as time goes on. A detailed fit of the time dependence of this shifting is given for a CDM model. The evolution of the cross correlation of the density fluctuation field evaluated at two different times is also studied. The possible implications of the results for the analysis of the observed correlation function of high redshift galaxies are discussed., revised version accepted for publication in MNRAS. The comparison with the scaling ansatz of Baugh & Gaztanaga (1996), one figure and some comments have been added. LaTex, 14 pages, 8 figures

We measure the spatial clustering of galaxies as a function of their morphological type at z~0.8, for the first time in a deep redshift survey with full morphological information. This is obtained by combining high-resolution HST imaging and VLT spectroscopy for about 8,500 galaxies to I_AB=22.5 with accurate spectroscopic redshifts from the zCOSMOS-Bright redshift survey. At this epoch, early-type galaxies already show a significantly stronger clustering than late-type galaxies on all probed scales. A comparison to the SDSS at z~0.1, shows that the relative clustering strength between early and late morphological classes tends to increase with cosmic time at small separations, while on large scales it shows no significant evolution since z~0.8. This suggests that most early-type galaxies had already formed in intermediate and dense environments at this epoch. Our results are consistent with a picture in which the relative clustering of different morphological types between z~1 and z~0, reflects the evolving role of environment in the morphological transformation of galaxies, on top of the global mass-driven evolution.