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

Radial redshift-space distortions due to peculiar velocities and other light-cone effects shape the maps we build of the Universe. We address the open question of their impact onto the monopole moment of the galaxy power spectrum, $P_0(k)$. Specifically, we use an upgraded numerical implementation of the LIGER method to generate $140$ mock galaxy density fields for a full Euclid-like survey and we measure $P_0(k)$ in each of them utilising a standard estimator. We compare the spectra obtained by turning on and off different effects. Our results show that wide-angle effects due to radial peculiar velocities generate excess power above the level expected within the plane-parallel approximation. They are detectable with a signal-to-noise ratio of 2.7 for $k, Comment: 20 pages, 15 figures. Minor changes to match the published version

We present a joint likelihood analysis of the halo power spectrum and bispectrum in real space. We take advantage of a large set of numerical simulations and of an even larger set of halo mock catalogs to provide a robust estimate of the covariance properties. We derive constraints on bias and cosmological parameters assuming a theoretical model from perturbation theory at one-loop for the power spectrum and tree-level for the bispectrum. By means of the Deviance Information Criterion, we select a reference bias model dependent on seven parameters that can describe the data up to $k_{\rm max,P}=0.3\, h \, {\rm Mpc}^{-1}$ for the power spectrum and $k_{\rm max,B}=0.09\, h \, {\rm Mpc}^{-1}$ for the bispectrum at redshift $z=1$. This model is able to accurately recover three selected cosmological parameters even for the rather extreme total simulation volume of $1000\, h^{-3} \, {\rm Gpc}^3$. With the same tools, we study how relations among bias parameters can improve the fit while reducing the parameter space. In addition, we compare common approximations to the covariance matrix against the full covariance estimated from the mocks, and quantify the (non-negligible) effect of ignoring the cross-covariance between the two statistics. Finally, we explore different selection criteria for the triangular configurations to include in the analysis, showing that excluding nearly equilateral triangles rather than simply imposing a fixed maximum $k_{\rm max,B}$ on all triangle sides can lead to a better exploitation of the information contained in the bispectrum., 38 pages, 18 figures, 1 table. Added direct comparison between power spectrum and joint analysis; further minor revisions

The mean and the scatter of the HI content of a dark-matter halo as a function of the halo mass are useful statistics that can be used to test models of structure and galaxy formation. We investigate the possibility of constraining this HI-halo mass relation (HIHMR) from intensity maps of the redshifted 21 cm line. In particular, we use the geometry and topology of the brightness-temperature isocontours in a single frequency channel as quantified by the Minkowski functionals. First, we generate mock maps from a large N-body simulation considering the impact of thermal noise and foreground removal. We then use the Fisher information formalism to forecast constraints on a parametric model for the HIHMR. We consider a 20,000 deg$^2$ survey (originally proposed for dark-energy science) conducted with the Square Kilometre Array Phase 1 (SKA-1) MID observatory operating in single-dish mode. For a channel bandwidth of 2 MHz, we show that an integration time of a few$\,\times\,10^4$ s per pointing is sufficient to image the smoothed HI distribution at redshift $z \simeq 1$ and to measure the HIHMR in a nearly optimal way from the Minkowski functionals. Tighter constraints on some of the parameters can be obtained by using also an independent measurement of the mean HI density. Combining the results from different frequency channels provides exquisite constraints on the evolution of the HIHMR, especially in the central frequency range of the data cube., 13 pages, 9 figures. Section 6 added to address the comments of the reviewer, no changes in the results

We provide a comparison between the matter bispectrum derived with different flavours of perturbation theory at next-to-leading order and measurements from an unprecedentedly large suite of $N$-body simulations. We use the $\chi^2$ goodness-of-fit test to determine the range of accuracy of the models as a function of the volume covered by subsets of the simulations. We find that models based on the effective-field-theory (EFT) approach have the largest reach, standard perturbation theory has the shortest, and `classical' resummed schemes lie in between. The gain from EFT, however, is less than in previous studies. We show that the estimated range of accuracy of the EFT predictions is heavily influenced by the procedure adopted to fit the amplitude of the counterterms. For the volumes probed by galaxy redshift surveys, our results indicate that it is advantageous to set three counterterms of the EFT bispectrum to zero and measure the fourth from the power spectrum. We also find that large fluctuations in the estimated reach occur between different realisations. We conclude that it is difficult to unequivocally define a range of accuracy for the models containing free parameters. Finally, we approximately account for systematic effects introduced by the $N$-body technique either in terms of a scale- and shape-dependent bias or by boosting the statistical error bars of the measurements (as routinely done in the literature). We find that the latter approach artificially inflates the reach of EFT models due to the presence of tunable parameters., Comment: 20 pages, 13 figures

Higher-order clustering statistics, like the galaxy bispectrum, can add complementary cosmological information to what is accessible with two-point statistics, like the power spectrum. While the standard way of measuring the bispectrum involves estimating a bispectrum value in a large number of Fourier triangle bins, the compressed modal bispectrum approximates the bispectrum as a linear combination of basis functions and estimates the expansion coefficients on the chosen basis. In this work, we compare the two estimators by using parallel pipelines to analyze the real-space halo bispectrum measured in a suite of $N$-body simulations corresponding to a total volume of $\sim 1{,}000 \,h^{-3}\,{\rm Gpc}^3$, with covariance matrices estimated from 10,000 mock halo catalogs. We find that the modal bispectrum yields constraints that are consistent and competitive with the standard bispectrum analysis: for the halo bias and shot noise parameters within the tree-level halo bispectrum model up to $k_{\rm max} \approx 0.06 \, (0.10) \,h\,{\rm Mpc}^{-1}$, only 6 (10) modal expansion coefficients are necessary to obtain constraints equivalent to the standard bispectrum estimator using $\sim$ 20 to 1,600 triangle bins, depending on the bin width. For this work, we have implemented a modal estimator pipeline using Markov Chain Monte Carlo simulations for the first time, and we discuss in detail how the parameter posteriors and modal expansion are robust to, or sensitive to, several user settings within the modal bispectrum pipeline. The combination of the highly efficient compression that is achieved and the large number of mock catalogs available allows us to quantify how our modal bispectrum constraints depend on the number of mocks that are used to estimate covariance matrices and the functional form of the likelihood., 51 pages, 20 figures, 2 tables. To be submitted to JCAP. Comments are welcome

Modelling the molecular gas that is routinely detected through CO observations of high-redshift galaxies constitutes a major challenge for ab initio simulations of galaxy formation. We carry out a suite of cosmological hydrodynamic simulations to compare three approximate methods that have been used in the literature to track the formation and evolution of the simplest and most abundant molecule, H$_2$. Namely, we consider: i) a semi-empirical procedure that associates H$_2$ to dark-matter haloes based on a series of scaling relations inferred from observations, ii) a model that assumes chemical equilibrium between the H$_2$ formation and destruction rates, and iii) a model that fully solves the out-of-equilibrium rate equations and accounts for the unresolved structure of molecular clouds. We study the impact of finite spatial resolution and show that robust H$_2$ masses at redshift $z\approx 4$ can only be obtained for galaxies that are sufficiently metal enriched in which H$_2$ formation is fast. This corresponds to H$_2$ reservoirs with masses $M_{\mathrm{H_2}}\gtrsim 6\times 10^9 \mathrm{M}_\odot$. In this range, equilibrium and non-equilibrium models predict similar molecular masses (but different galaxy morphologies) while the semi-empirical method produces less H$_2$. The star formation rates as well as the stellar and H$_2$ masses of the simulated galaxies are in line with those observed in actual galaxies at similar redshifts that are not massive starbursts. The H$_2$ mass functions extracted from the simulations at $z\approx 4$ agree well with recent observations that only sample the high-mass end. However, our results indicate that most molecular material at high $z$ lies yet undetected in reservoirs with $10^9, Comment: 19 pages, 12 figures, 1 table, accepted for publication in MNRAS, comments are welcome

The forthcoming generation of galaxy redshift surveys will sample the large-scale structure of the Universe over unprecedented volumes with high-density tracers. This advancement will make robust measurements of three-point clustering statistics possible. In preparation for this improvement, we investigate how several methodological choices can influence inferences based on the bispectrum about galaxy bias and shot noise. We first measure the real-space bispectrum of dark-matter haloes extracted from 298 N-body simulations covering a volume of approximately $1000 h^{-3} \mathrm{Gpc}^3$. We then fit a series of theoretical models based on tree-level perturbation theory to the numerical data. To achieve this, we estimate the covariance matrix of the measurement errors by using 10,000 mock catalogues generated with the Pinocchio code. We study how the model constraints are influenced by the binning strategy for the bispectrum configurations and by the form of the likelihood function. We also use Bayesian model-selection techniques to single out the optimal theoretical description of our data. We find that a three-parameter bias model combined with Poissonian shot noise is necessary to model the halo bispectrum up to scales of $k_\mathrm{max}\lesssim 0.08 h \mathrm{Mpc}^{-1}$, although fitting formulae that relate the bias parameters can be helpful to reduce the freedom of the model without compromising accuracy. Our data clearly disfavour local Eulerian and local Lagrangian bias models and do not require corrections to Poissonian shot noise. We anticipate that model-selection diagnostics will be particularly useful to extend the analysis to smaller scales as, in this case, the number of model parameters will grow significantly large., 38 pages, 17 figures; major revisions: inclusion of open triangle bins in all results, new priors for the shot-noise parameters, and additional comparison with Gaussian variance

Starting from first principles, we derive the fundamental equations that relate the $n$-point correlation functions in real and redshift space. Our result generalises the so-called `streaming model' to higher-order statistics: the full $n$-point correlation in redshift-space is obtained as an integral of its real-space counterpart times the joint probability density of $n-1$ relative line-of-sight peculiar velocities. Equations for the connected $n$-point correlation functions are obtained by recursively applying the generalised streaming model for decreasing $n$. Our results are exact within the distant-observer approximation and completely independent of the nature of the tracers for which the correlations are evaluated. Focusing on 3-point statistics, we use an $N$-body simulation to study the joint probability density function of the relative line-of-sight velocities of pairs of particles in a triplet. On large scales, we find that this distribution is approximately Gaussian and that its moments can be accurately computed with standard perturbation theory. We use this information to formulate a phenomenological 3-point Gaussian streaming model. A practical implementation is obtained by using perturbation theory at leading order to approximate several statistics in real space. In spite of this simplification, the resulting predictions for the matter 3-point correlation function in redshift space are in rather good agreement with measurements performed in the simulation. We discuss the limitations of the simplified model and suggest a number of possible improvements. Our results find direct applications in the analysis of galaxy clustering but also set the basis for studying 3-point statistics with future peculiar-velocity surveys and experiments based on the kinetic Sunyaev-Zel'dovich effect., Comment: 48 pages, 16 figures, accepted for publication in JCAP

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 study a model of cosmic reionization where quasars (QSOs) are the dominant source of ionizing photons at all relevant epochs. We employ a suite of adaptive hydrodynamical simulations post-processed with a multi-wavelength Monte Carlo radiative-transfer code and calibrate them in order to accurately reproduce the observed quasar luminosity function and emissivity evolution. Our results show that the QSO-only model fails in reproducing key observables linked to the Helium reionization, as the temperature evolution of the inter-galactic medium (IGM) and the HeII effective optical depth in synthetic Lyα spectra. Nevertheless, we find hints that an increased quasar contribution can explain recent measurements of a large inhomogeneity in the IGM at redshift z ≈ 5. Finally, we devise a method capable of constraining the QSOs contribution to the reionization from the properties of the HeII Lyα forest at z ≈ 3.5.

The interpretation of redshift surveys requires modeling the relationship between large-scale fluctuations in the observed number density of tracers, $\delta_\mathrm{h}$, and the underlying matter density, $\delta$. Bias models often express $\delta_\mathrm{h}$ as a truncated series of integro-differential operators acting on $\delta$, each weighted by a bias parameter. Due to the presence of `composite operators' (obtained by multiplying fields evaluated at the same spatial location), the linear bias parameter measured from clustering statistics does not coincide with that appearing in the bias expansion. This issue can be cured by re-writing the expansion in terms of `renormalised' operators. After providing a pedagogical and comprehensive review of bias renormalisation in perturbation theory, we generalize the concept to non-perturbative dynamics and successfully apply it to dark-matter haloes extracted from a large suite of N-body simulations. When comparing numerical and perturbative results, we highlight the effect of the window function employed to smooth the random fields. We then measure the bias parameters as a function of halo mass by fitting a non-perturbative bias model (both before and after applying renormalisation) to the cross spectrum $P_{\delta_\mathrm{h}\delta}(k)$. Finally, we employ Bayesian model selection to determine the optimal operator set to describe $P_{\delta_\mathrm{h}\delta}(k)$ for $k, Comment: 22 pages, 10 figures, 8 tables. Accepted for publication in MNRAS

We forecast the potential of the forthcoming X-ray galaxy-cluster survey with eROSITA to constrain dark-energy models. We focus on spatially-flat cosmological scenarios with either constant or time-dependent dark-energy equation-of-state parameters. Fisher information is extracted from the number density and spatial clustering of a photon-count-limited sample of clusters of galaxies up to z~2. We consider different scenarios for the availability of (i) X-ray follow-up observations, (ii) photometric and spectroscopic redshifts, and (iii) accurate knowledge of the observable -- mass relation down to the scale of galaxy groups. With about 125,000 clusters (detected with more than 50 photons and with mass M500c > $10^{13} h^{-1}$ Msun) from an average all-sky exposure of 1.6 ks, eROSITA will give marginalized, one-dimensional, 1-$\sigma$ errors of $\Delta \sigma_8 = \pm~0.008$ ($\sim$1 per cent), $\Delta \Omega_{\rm M} = \pm~0.006$ (2.2 per cent), $\Delta w_0 = \pm~0.07$ (7 per cent), and $\Delta w_a = \pm~0.25$ (optimistic scenario) in combination with (and largely improving upon) current constraints from various cosmological probes (cosmic microwave background, BAOs, Type Ia SNe). Our findings correspond to a dark-energy figure of merit in the range of $116-162$ (after the four years of all-sky survey), making eROSITA one of the first Stage IV experiments to come on line according to the classification of the Dark Energy Task Force. To secure improved mass calibrations and to include high-redshift clusters (z > 0.5) as well as objects at the group-mass scale (M500c < 5 $\times~ 10^{13} h^{-1}$ Msun) will be vital to reach such accuracies., Comment: Submitted to MNRAS. Main results at Figure 2 and Table 5

We use the Fisher-matrix formalism to investigate whether the galaxy bispectrum in redshift space, $B$, contains additional cosmological information with respect to the power spectrum, $P$. We focus on a $Euclid$-like survey and consider cosmological models dominated by dark energy and cold dark matter with Gaussian primordial perturbations. After discussing the phenomenology of redshift-space distortions for the bispectrum, we derive an expression for the cross-covariance between $B$ and $P$ at leading order in perturbation theory. Our equation generalizes previous results that did not consider binning in the orientation of wavevector triangles with respect to the line of sight. By considering Fourier modes with wavenumber $k, 24 pages, 15 figures, accepted for publication in MNRAS. Minor changes to match the published version. Includes updated HOD model for Euclid galaxies. Conclusions unchanged

Despite a history that dates back at least a quarter of a century studies of voids in the large--scale structure of the Universe are bedevilled by a major problem: there exist a large number of quite different void--finding algorithms, a fact that has so far got in the way of groups comparing their results without worrying about whether such a comparison in fact makes sense. Because of the recent increased interest in voids, both in very large galaxy surveys and in detailed simulations of cosmic structure formation, this situation is very unfortunate. We here present the first systematic comparison study of thirteen different void finders constructed using particles, haloes, and semi--analytical model galaxies extracted from a subvolume of the Millennium simulation. The study includes many groups that have studied voids over the past decade. We show their results and discuss their differences and agreements. As it turns out, the basic results of the various methods agree very well with each other in that they all locate a major void near the centre of our volume. Voids have very underdense centres, reaching below 10 percent of the mean cosmic density. In addition, those void finders that allow for void galaxies show that those galaxies follow similar trends. For example, the overdensity of void galaxies brighter than $m_B = -20 $ is found to be smaller than about -0.8 by all our void finding algorithms., Comment: 14 pages, version accepted by MNRAS

We introduce a method to create mock galaxy catalogues in redshift space including general relativistic effects to linear order in the cosmological perturbations. We dub our method LIGER, short for `light cones with general relativity'. LIGER takes a (N-body or hydrodynamic) Newtonian simulation as an input and outputs the distribution of galaxies in comoving redshift space. This result is achieved making use of a coordinate transformation and simultaneously accounting for lensing magnification. The calculation includes both local corrections and terms that have been integrated along the line of sight. Our fast implementation allows the production of many realizations that can be used to forecast the performance of forthcoming wide-angle surveys and to estimate the covariance matrix of the observables. To facilitate this use, we also present a variant of LIGER designed for large-volume simulations with low mass resolution. In this case, the galaxy distribution on large scales is obtained by biasing the matter-density field. Finally, we present two sample applications of LIGER. First, we discuss the impact of weak gravitational lensing onto the angular clustering of galaxies in a Euclid-like survey. In agreement with previous analytical studies, we find that magnification bias can be measured with high confidence. Second, we focus on two generally neglected Doppler-induced effects: magnification and the change of number counts with redshift. We show that the corresponding redshift-space distortions can be detected at 5.5$\sigma$ significance with the completed Square Kilometre Array., Comment: 16 pages, 12 figures, comments welcome

The growth-rate and the internal dynamics of galaxy-sized dark-matter haloes depend on their location within the cosmic web. Haloes that sit at the nodes grow in mass till the present time and are dominated by radial orbits. Conversely, haloes embedded in prominent filaments do not change much in size and are dominated by tangential orbits. Using zoom hydrodynamical simulations including star formation and feedback, we study how gas accretes onto these different classes of objects that, for simplicity, we dub 'accreting' and 'stalled' haloes. We find that all haloes get a fresh supply of newly accreted gas in their inner regions, although this slowly decreases with time, in particular for the stalled haloes. The inflow of new gas is always higher than (but comparable with) that of recycled material. Overall, the cold-gas fraction increases (decreases) with time for the accreting (stalled) haloes. In all cases, a stellar disc and a bulge form at the centre of the simulated haloes. The total stellar mass is in excellent agreement with expectations based on the abundance-matching technique. Many properties of the central galaxies do not seem to correlate with the large-scale environment in which the haloes reside. However, there are two notable exceptions that characterise stalled haloes with respect to their accreting counterparts: i) the galaxy disc contains much older stellar populations; ii) its vertical scale-height is larger by a factor of two or more. This thickening is likely due to the heating of the long-lived discs by mergers and close flybys., 17 pages, 21 figures. Accepted for publication in MNRAS

We use zoom hydrodynamical simulations to investigate the properties of satellites within galaxy-sized dark-matter haloes with different assembly histories. We consider two classes of haloes at redshift $z=0$: `stalled' haloes that assembled at $z>1$ and `accreting' ones that are still forming nowadays. Previously, we showed that the stalled haloes are embedded within thick filaments of the cosmic web while the accreting ones lie where multiple thin filaments converge. We find that satellites in the two classes have both similar and different properties. Their mass spectra, radial count profiles, baryonic and stellar content, and the amount of material they shed are indistinguishable. However, the mass fraction locked in satellites is substantially larger for the accreting haloes as they experience more mergers at late times. The largest difference is found in the satellite kinematics. Substructures fall towards the accreting haloes along quasi-radial trajectories whereas an important tangential velocity component is developed, before accretion, while orbiting the filament that surrounds the stalled haloes. Thus, the velocity anisotropy parameter of the satellites ($\beta$) is positive for the accreting haloes and negative for the stalled ones. This signature enables us to tentatively categorize the Milky Way halo as stalled based on a recent measurement of $\beta$. Half of our haloes contain clusters of satellites with aligned orbital angular momenta corresponding to flattened structures in space. These features are not driven by baryonic physics and are only found in haloes hosting grand-design spiral galaxies, independently of their assembly history., Comment: 18 pages, 15 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 clustering of dark matter haloes with fixed mass depends on their formation history, an effect known as assembly bias. We use zoom N -body simulations to investigate the origin of this phenomenon. For each halo at redshift z=0, we determine the time in which the physical volume containing its final mass becomes stable. We consider five examples for which this happens at z~1.5 and two that do not stabilize by z=0. The zoom simulations show that early-collapsing haloes do not grow in mass at z=0 while late-forming ones show a net inflow. The reason is that 'accreting' haloes are located at the nodes of a network of thin filaments feeding them. Conversely, each 'stalled' halo lies within a prominent filament that is thicker than the halo size. Infalling material from the surroundings becomes part of the filament while matter within it recedes from the halo. We conclude that assembly bias originates from quenching halo growth due to tidal forces following the formation of non-linear structures in the cosmic web, as previously conjectured in the literature. Also the internal dynamics of the haloes change: the velocity anisotropy profile is biased towards radial (tangential) orbits in accreting (stalled) haloes. Our findings reveal the cause of the yet unexplained dependence of halo clustering on the anisotropy. Finally, we extend the excursion-set theory to account for these effects. A simple criterion based on the ellipticity of the linear tidal field combined with the spherical collapse model provides excellent predictions for both classes of haloes., Comment: 20 pages, 18 figures, Accepted for publication in MNRAS

Constrained realisations of Gaussian random fields are used in cosmology to design special initial conditions for numerical simulations. We review this approach and its application to density peaks providing several worked-out examples. We then critically discuss the recent proposal to use constrained realisations to modify the linear density field within and around the Lagrangian patches that form dark-matter haloes. The ambitious concept is to forge `genetically modified' haloes with some desired properties after the non-linear evolution. We demonstrate that the original implementation of this method is not exact but approximate because it tacitly assumes that protohaloes sample a set of random points with a fixed mean overdensity. We show that carrying out a full genetic modification is a formidable and daunting task requiring a mathematical understanding of what determines the biased locations of protohaloes in the linear density field. We discuss approximate solutions based on educated guesses regarding the nature of protohaloes. We illustrate how the excursion-set method can be adapted to predict the non-linear evolution of the modified patches and thus fine tune the constraints that are necessary to obtain preselected halo properties. This technique allows us to explore the freedom around the original algorithm for genetic modification. We find that the quantity which is most sensitive to changes is the halo mass-accretion rate at the mass scale on which the constraints are set. Finally we discuss constraints based on the protohalo angular momenta., Comment: 16 pages, 5 figures, 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 explore the inter-relationships between mass, star-formation rate and environment in the SDSS, zCOSMOS and other surveys. The differential effects of mass and environment are completely separable to z ~ 1, indicating that two distinct processes are operating, "mass-quenching" and "environment-quenching". Environment-quenching, at fixed over-density, evidently does not change with epoch to z ~ 1, suggesting that it occurs as large-scale structure develops in the Universe. The observed constancy of the mass-function shape for star-forming galaxies, demands that the mass-quenching of galaxies around and above M*, must be proportional to their star-formation rates at all z < 2. We postulate that this simple mass-quenching law also holds over a much broader range of stellar mass and epoch. These two simple quenching processes, plus some additional quenching due to merging, then naturally produce (a) a quasi-static Schechter mass function for star-forming galaxies with a value of M* that is set by the proportionality between the star-formation and mass-quenching rates, (b) a double Schechter function for passive galaxies with two components: the dominant one is produced by mass-quenching and has exactly the same M* as the star-forming galaxies but an alpha shallower by +1, while the other is produced by environment effects and has the same M* and alpha as the star-forming galaxies, and is larger in high density environments. Subsequent merging of quenched galaxies modifies these predictions somewhat in the denser environments, slightly increasing M* and making alpha more negative. All of these detailed quantitative relationships between the Schechter parameters are indeed seen in the SDSS, lending strong support to our simple empirically-based model. The model naturally produces for passive galaxies the "anti-hierarchical" run of mean ages and alpha-element abundances with mass., 66 pages, 19 figures, 1 movie, accepted for publication in ApJ. The movie is also available at http://www.exp-astro.phys.ethz.ch/zCOSMOS/MF_simulation_d1_d4.mov

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 study the evolution of galaxies inside and outside of the group environment since z=1 using a large well defined set of groups and galaxies from the zCOSMOS-bright redshift survey in the COSMOS field. The fraction of galaxies with early-type morphologies increases monotonically with M_B luminosity and stellar mass and with cosmic epoch. It is higher in the groups than elsewhere, especially at later epochs. The emerging environmental effect is superposed on a strong global mass-driven evolution, and at z~0.5 and log(M*/Msol)~10.2, the "effect" of group environment is equivalent to (only) about 0.2 dex in stellar mass or 2 Gyr in time. The stellar mass function of galaxies in groups is enriched in massive galaxies. We directly determine the transformation rates from late to early morphologies, and for transformations involving colour and star formation indicators. The transformation rates are systematically about twice as high in the groups as outside, or up to 3-4 times higher correcting for infall and the appearance of new groups. The rates reach values, for masses around the crossing mass 10^10.5 Msol, as high as (0.3-0.7)/Gyr in the groups, implying transformation timescales of 1.4-3 Gyr, compared with less than 0.2/Gyr, i.e. timescales >5 Gyr, outside of groups. All three transformation rates decrease at higher stellar masses, and must decrease also at the lower masses below 10^10 Msol which we cannot well probe. The rates involving colour and star formation are consistently higher than those for morphology, by a factor of about 50%. Our conclusion is that the transformations which drive the evolution of the overall galaxy population since z~1 must occur at a rate 2-4 times higher in groups than outside of them., 21 pages, 13 figures, submitted to ApJ

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

One of the key unanswered questions in the study of galaxy evolution is what physical processes inside galaxies drive the changes in the SFRs in individual galaxies that, taken together, produce the large decline in the global star-formation rate density (SFRD) to redshifts since z~2. Many studies of the SFR at intermediate redshifts have been made as a function of the integrated stellar mass of galaxies but these did not use information on the internal structural properties of the galaxies. In this paper we present a comparative study of the dependence of SFRs on the average surface mass densities (SigmaM) of galaxies of different morphological types up to z~1 using the zCOSMOS and SDSS surveys. The main findings about the evolution of these relatively massive galaxies are: 1) There is evidence that, for both SDSS ans zCOSMOS galaxies, the mean specific SFR within a given population (either disk-dominated or bulge-dominated) is independent of SigmaM; 2) The observed SSFR - SigmaM step-function relation is due, at all investigated redshifts, to the changing mix of disk-dominated and bulge-dominated galaxies as surface density increases and the strong difference in the average SSFR between disks and bulges. We also find a modest differential evolution in the size-mass relations of disk and spheroid galaxies; 3) The shape of the median SSFR - SigmaM relation is similar, but with median SSFR values that are about 5-6 times higher in zCOSMOS galaxies than for SDSS, across the whole range of SigmaM, and in both spheroid and disk galaxies. This increase matches that of the global SFRD of the Universe as a whole, emphasizing that galaxies of all types are contributing, proportionally, to the global increase in SFRD in the Universe back to these redshifts (abridged)., Comment: Published 2009 in ApJ, 694, 1099

International audience; Motivated by the desire to reliably and automatically classify structure of thousands of COSMOS galaxies, we present ZEST, the Zurich Estimator of Structural Types. To classify galaxy structure, ZEST uses (1) five nonparametric diagnostics: asymmetry, concentration, Gini coefficient, second-order moment of the brightest 20% of galaxy pixels, and ellipticity; and (2) the exponent n of single-Sérsic fits to the two-dimensional surface brightness distributions. To fully exploit the wealth of information while reducing the redundancy present in these diagnostics, ZEST performs a principal component (PC) analysis. We use a sample of ~56,000 IAB1, PC2, and PC3. We demonstrate the robustness of the ZEST grid on the z=0 sample of Frei et al. The ZEST classification breaks most of the degeneracy between different galaxy populations that affects morphological classifications based on only some of the diagnostics included in ZEST. As a first application, we present the evolution since z~1 of the luminosity functions (LFs) of COSMOS galaxies of early, disk, and irregular galaxies and, for disk galaxies, of different bulge-to-disk ratios. Overall, we find that the LF up to a redshift z=1 is consistent with a pure luminosity evolution (of about 0.95 mag at z~0.7). We highlight, however, two trends that are in general agreement with a downsizing scenario for galaxy formation, i.e., (1) a deficit of a factor of about 2 at z~0.7 of MB>-20.5 structurally classified early-type galaxies and (2) an excess of a factor of about 3, at a similar redshift, of irregular galaxies. Based on observations with the NASA/ESA Hubble Space Telescope , obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy (AURA), Inc., under NASA contract NAS 5-26555 also based on data collected at the Subaru Telescope, which is operated by the National Astronomical Observatory of Japan; Kitt Peak National Observatory, Cerro Tololo Inter-American Observatory, and the National Optical Astronomy Observatory, which are operated by AURA, Inc., under cooperative agreement with the National Science Foundation; and the Canada-France-Hawaii Telescope operated by the National Research Council of Canada, the Centre National de la Recherche Scientifique de France, and the University of Hawaii.

Brightness-temperature fluctuations in the redshifted 21-cm background from the cosmic dark ages are generated by irregularities in the gas-density distribution and can then be used to determine the statistical properties of density fluctuations in the early Universe. We first derive the most general expansion of brightness-temperature fluctuations up to second order in terms of all the possible sources of spatial fluctuations. We then focus on the three-point statistics and compute the angular bispectrum of brightness-temperature fluctuations generated prior to the epoch of hydrogen reionization. For simplicity, we neglect redshift-space distortions. We find that low-frequency radio experiments with arcmin angular resolution can easily detect non-Gaussianity produced by non-linear gravity with high signal-to-noise ratio. The bispectrum thus provides a unique test of the gravitational instability scenario for structure formation, and can be used to measure the cosmological parameters. Detecting the signature of primordial non-Gaussianity produced during or right after an inflationary period is more challenging but still possible. An ideal experiment limited by cosmic variance only and with an angular resolution of a few arcsec has the potential to detect primordial non-Gaussianity with a non-linearity parameter of f_NL ~ 1. Additional sources of error as weak lensing and an imperfect foreground subtraction could severely hamper the detection of primordial non-Gaussianity which will benefit from the use of optimal estimators combined with tomographic techniques.

We present a method to test the isotropy of the magnitude-redshift relation of Type Ia Supernovae (SNe Ia) and single out the most discrepant direction (in terms of the signal-to-noise ratio) with respect to the all-sky data. Our technique accounts for possible directional variations of the corrections for SNe Ia and yields all-sky maps of the best-fit cosmological parameters with arbitrary angular resolution. To show its potential, we apply our method to the recent Union2.1 compilation, building maps with three different angular resolutions. We use a Monte Carlo method to estimate the statistical significance with which we could reject the null hypothesis that the magnitude-redshift relation is isotropic based on the properties of the observed most discrepant directions. We find that, based on pure signal-to-noise arguments, the null hypothesis cannot be rejected at any meaningful confidence level. However, if we also consider that the strongest deviations in the Union2.1 sample closely align with the dipole temperature anisotropy of the cosmic microwave background, we find that the null hypothesis should be rejected at the $95-99$ per cent confidence level, slightly depending on the angular resolution of the study. If this result is not due to a statistical fluke, it might either indicate that the SN data have not been cleaned from all possible systematics or even point towards new physics. We finally discuss future perspectives in the field for achieving larger and more uniform data sets that will vastly improve the quality of the results and optimally exploit our method., Accepted for publication in The Astrophysical Journal. 10 pages, 8 figures, 4 tables. Comments are welcome

Theoretical and observational arguments suggest that the intergalactic medium (IGM) might have been polluted with metals produced by early star formation. In this scenario, Lyman-break galaxies (LBGs) at redshift z=3 are likely to be surrounded by old metal bubbles, relics of an era when the characteristic mass of galaxies was small and gas retainment more difficult. We find that pregalactic enrichment of the IGM from 10^8-10^{10} M_sun dwarf galaxies at 6, Comment: 4 pages, 2 figures, ApJL in press

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

We present a new measure of the angular two-point correlation function of Lyman-break galaxies (LBGs) at z~3, obtained from the variance of galaxy counts in 2-dimensional cells. By avoiding binning of the angular separations, this method is significantly less affected by shot noise than traditional measures, and allows for a more accurate determination of the correlation function. We used a sample of about 1,000 galaxies with R, Comment: 67 pages, 16 figures, accepted for publication in the Astrophysical Journal

We use a high-resolution hydrodynamic simulation that tracks the non-equilibrium abundance of molecular hydrogen within a massive high-redshift galaxy to produce mock Atacama Large Millimeter Array (ALMA) maps of the fine-structure lines of atomic carbon, CI 1-0 and CI 2-1. Inspired by recent observational and theoretical work, we assume that CI is thoroughly mixed within giant molecular clouds and demonstrate that its emission is an excellent proxy for H2. Nearly all of the H2 associated with the galaxy can be detected at redshifts z, 5 pages, 4 figures, accepted for publication in MNRAS Letters

We introduce a sub-grid model for the non-equilibrium abundance of molecular hydrogen in cosmological simulations of galaxy formation. We improve upon previous work by accounting for the unresolved structure of molecular clouds in a phenomenological way which combines both observational and numerical results on the properties of the turbulent interstellar medium. We apply the model to a cosmological simulation of the formation of a Milky-Way-sized galaxy at z=2, and compare the results to those obtained using other popular prescriptions that compute the equilibrium abundance of H2. In these runs we introduce an explicit link between star formation and the local H2 abundance, and perform an additional simulation in which star formation is linked directly to the density of cold gas. In better agreement with observations, we find that the simulated galaxy produces less stars and harbors a larger gas reservoir when star formation is regulated by molecular hydrogen. In this case, the galaxy is composed of a younger stellar population as early star formation is inhibited in small, metal poor dark-matter haloes which cannot efficiently produce H2. The number of luminous satellites orbiting within the virial radius of the galaxy at z=2 is reduced by 10-30 per cent in models with H2-regulated star formation., 18 pages, 12 figures. Accepted for publication in MNRAS

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)