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.

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.

The thermal Sunyaev-Zel'dovich (tSZ) effect is one of the primary tools for finding and characterizing galaxy clusters. Several ground-based experiments are either underway or are being planned for mapping wide areas of the sky at $\sim 150$ GHz with large-aperture telescopes. We present cosmological forecasts for a 'straw man' tSZ survey that will observe a sky area between $200$ and $10^4$ deg$^2$ to an rms noise level between 2.8 and 20.2 $\mu$K-arcmin. The probes we consider are the cluster number counts (as a function of the integrated Compton-$Y$ parameter and redshift) and their angular clustering (as a function of redshift). At fixed observing time, we find that wider surveys constrain cosmology slightly better than deeper ones due to their increased ability to detect rare high-mass clusters. In all cases, we notice that adding the clustering information does not practically improve the constraints derived from the number counts. We compare forecasts obtained by sampling the posterior distribution with the Markov-chain-Monte-Carlo method against those derived using the Fisher-matrix formalism. We find that the latter produces slightly optimistic constraints where errors are underestimated at the 10 per cent level. Most importantly, we use an analytic method to estimate the selection function of the survey and account for its response to variations of the cosmological parameters in the likelihood function. Our analysis demonstrates that neglecting this effect (as routinely done in the literature) yields artificially tighter constraints by a factor of 2.2 and 1.7 for $\sigma_8$ and $\Omega_\mathrm{M}$, respectively.

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 in order 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<M_{\mathrm H_2}<10^{10} \mathrm{M}_\odot$.

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<0.2\,h$ Mpc$^{-1}$ at redshift $z=0$. We find that it includes $\delta, \nabla^2\delta, \delta^2$ and the square of the traceless tidal tensor, $s^2$. Considering higher-order terms (in $\delta$) leads to overfitting as they cannot be precisely constrained by our data. We also notice that next-to-leading-order perturbative solutions are inaccurate for $k\gtrsim 0.1\,h$ Mpc$^{-1}$.

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.

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<0.15 \,h$ Mpc$^{-1}$, we find that $B$ and $P$ set similar constraints on the cosmological parameters. Generally, error bars moderately improve when the two probes are combined together. For instance, the joint 68.3 per cent credible region for the parameters that describe a dynamical dark-energy equation of state shrinks by a factor of 2.6 with respect to only using the power spectrum. Regrettably, this improvement is cancelled out when the clustering analysis is combined with priors based on current studies of the cosmic microwave background. In this case, combining $B$ and $P$ does not give any appreciable benefit other than allowing a precise determination of galaxy bias. Finally, we discuss how results depend on the binning strategy for the clustering statistics as well as on the maximum wavenumber. We also show that only considering the bispectrum monopole leads to a significant loss of information.

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.

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.

The radial acceleration measured in bright galaxies tightly correlates with that generated by the observed distribution of baryons, a phenomenon known as the radial acceleration relation (RAR). Dwarf spheroidal satellite galaxies have been recently found to depart from the extrapolation of the RAR measured for more massive objects but with a substantially larger scatter. If confirmed by new data, this result provides a powerful test of the theory of gravity at low accelerations that requires robust theoretical predictions. By using high-resolution hydrodynamical simulations, we show that, within the standard model of cosmology ($\Lambda$CDM), satellite galaxies are expected to follow the same RAR as brighter systems but with a much larger scatter which does not correlate with the physical properties of the galaxies. In the simulations, the RAR evolves mildly with redshift. Moreover, the acceleration due to the gravitational field of the host has no effect on the RAR. This is in contrast with the External Field Effect in Modified Newtonian Dynamics (MOND) which causes galaxies in strong external fields to deviate from the RAR. This difference between $\Lambda$CDM and MOND offers a possible way to discriminate between them.

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.

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.

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.

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.

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.

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

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.

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.

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.

In the excursion set approach to structure formation initially spherical regions of the linear density field collapse to form haloes of mass $M$ at redshift $z_{\rm id}$ if their linearly extrapolated density contrast, averaged on that scale, exceeds some critical threshold, $\delta_{\rm c}(z_{\rm id})$. The value of $\delta_{\rm c}(z_{\rm id})$ is often calculated from the spherical or ellipsoidal collapse model, which provide well-defined predictions given auxiliary properties of the tidal field at a given location. We use two cosmological simulations of structure growth in a $\Lambda$ cold dark matter scenario to quantify $\delta_{\rm c}(z_{\rm id})$, its dependence on the surrounding tidal field, as well as on the shapes of the Lagrangian regions that collapse to form haloes at $z_{\rm id}$. Our results indicate that the ellipsoidal collapse model provides an accurate description of the mean dependence of $\delta_{\rm c}(z_{\rm id})$ on both the strength of the tidal field and on halo mass. However, for a given $z_{\rm id}$, $\delta_{\rm c}(z_{\rm id})$ depends strongly on the halo's characteristic formation redshift: the earlier a halo forms, the higher its initial density contrast. Surprisingly, the majority of haloes forming $today$ fall below the ellipsoidal collapse barrier, contradicting the model predictions. We trace the origin of this effect to the non-spherical shapes of Lagrangian haloes, which arise naturally due to the asymmetry of the linear tidal field. We show that a modified collapse model, that accounts for the triaxial shape of protohaloes, provides a more accurate description of the measured minimum overdensities of recently collapsed objects.

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}<3) can be found at high redshift, in agreement with the most recent observations of UV-transmitting quasars. At the highest redshift currently probed by observations (z~3.4), our updated model predicts a much lower HeII effective optical depth than previous simulations in the literature relieving most of the tension with the current data, that, however, still persists at about the (Gaussian) 1{\sigma} to 2{\sigma} level. Given the very small number of observed lines of sight, our analysis indicates that current data cannot rule out a purely AGN-driven scenario with high statistical significance.

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<4 using a compact interferometric configuration with a large synthesized beam (that does not resolve the target galaxy) in less than 4 h of integration time. Low-resolution imaging of the \CI lines (in which the target galaxy is resolved into three to four beams) will detect ~80 per cent of the H2 in less than 12 h of aperture synthesis. In this case, the resulting data cube also provides the crucial information necessary for determining the dynamical state of the galaxy. We conclude that ALMA observations of the CI 1-0 and 2-1 emission are well-suited for extending the interval of cosmic look-back time over which the H2 distributions, the dynamical masses, and the Tully-Fisher relation of galaxies can be robustly probed.

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.

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)

The halo mass function from N-body simulations of collisionless matter is generally used to retrieve cosmological parameters from observed counts of galaxy clusters. This neglects the observational fact that the baryonic mass fraction in clusters is a random variable that, on average, increases with the total mass (within an overdensity of 500). Considering a mock catalog that includes tens of thousands of galaxy clusters, as expected from the forthcoming generation of surveys, we show that the effect of a varying baryonic mass fraction will be observable with high statistical significance. The net effect is a change in the overall normalization of the cluster mass function and a milder modification of its shape. Our results indicate the necessity of taking into account baryonic corrections to the mass function if one wants to obtain unbiased estimates of the cosmological parameters from data of this quality. We introduce the formalism necessary to accomplish this goal. Our discussion is based on the conditional probability of finding a given value of the baryonic mass fraction for clusters of fixed total mass. Finally, we show that combining information from the cluster counts with measurements of the baryonic mass fraction in a small subsample of clusters (including only a few tens of objects) will nearly optimally constrain the cosmological parameters.

We investigate the accuracy of Eulerian perturbation theory for describing the matter and galaxy power spectra in real and redshift space in light of future observational probes for precision cosmology. Comparing the analytical results with a large suite of N-body simulations (160 independent boxes of 13.8 (Gpc/h)^3 volume each, which are publicly available), we find that re-summing terms in the standard perturbative approach predicts the real-space matter power spectrum with an accuracy of < 2% for k < 0.20 h/Mpc at redshifts z < 1.5. This is obtained following the widespread technique of writing the resummed propagator in terms of 1-loop contributions. We show that the accuracy of this scheme increases by considering higher-order terms in the resummed propagator. By combining resummed perturbation theories with several models for the mappings from real to redshift space discussed in the literature, the multipoles of the dark-matter power spectrum can be described with sub-percent deviations from N-body results for k < 0.15h/Mpc at z < 1. As a consequence, the logarithmic growth rate, f, can be recovered with sub-percent accuracy on these scales. Extending the models to massive dark-matter haloes in redshift space, our results describe the monopole term from N-body data within 2% accuracy for scales k < 0.15 h/Mpc at z < 0.5; here f can be recovered within < 5% when the halo bias is known. We conclude that these techniques are suitable to extract cosmological information from future galaxy surveys.

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 LCDM cosmology. Our analysis focuses on "proto-haloes", which 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 proto-haloes and estimate their auto spectral 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 and 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 velocity fields. We also show that the bias coefficients are accurately predicted for the velocity divergence. On the contrary, 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.

(Abridged) We analyze the environments and galactic properties (morphologies and star-formation histories) of a sample of 153 close kinematic pairs in the redshift range 0.2 < z < 1 identified in the zCOSMOS-bright 10k spectroscopic sample of galaxies. Correcting for projection effects, the fraction of close kinematic pairs is three times higher in the top density quartile than in the lowest one. This translates to a three times higher merger rate because the merger timescales are shown, from mock catalogues based on the Millennium simulation, to be largely independent of environment once the same corrections for projection is applied. We then examine the morphologies and stellar populations of galaxies in the pairs, comparing them to control samples that are carefully matched in environment so as to remove as much as possible the well-known effects of environment on the properties of the parent population of galaxies. Once the environment is properly taken into account in this way, we find that the early-late morphology mix is the same as for the parent population, but that the fraction of irregular galaxies is boosted by 50-75%, with a disproportionate increase in the number of irregular-irregular pairs (factor of 4-8 times), due to the disturbance of disk galaxies. Future dry-mergers, involving elliptical galaxies comprise less than 5% of all close kinematic pairs. In the closest pairs, there is a boost in the specific star-formation rates of star-forming galaxies of a factor of 2-4, and there is also evidence for an increased incidence of post star-burst galaxies. Although significant for the galaxies involved, the "excess" star-formation associated with pairs represents only about 5% of the integrated star-formation activity in the parent sample. Although most pair galaxies are in dense environments, the effects of interaction appear to be largest in the lower density environments.

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.

Our goal is to develop a new and reliable statistical method to classify galaxies from large surveys. We probe the reliability of the method by comparing it with a three-dimensional classification cube, using the same set of spectral, photometric and morphological parameters.We applied two different methods of classification to a sample of galaxies extracted from the zCOSMOS redshift survey, in the redshift range 0.5 < z < 1.3. The first method is the combination of three independent classification schemes, while the second method exploits an entirely new approach based on statistical analyses like Principal Component Analysis (PCA) and Unsupervised Fuzzy Partition (UFP) clustering method. The PCA+UFP method has been applied also to a lower redshift sample (z < 0.5), exploiting the same set of data but the spectral ones, replaced by the equivalent width of H$\alpha$. The comparison between the two methods shows fairly good agreement on the definition on the two main clusters, the early-type and the late-type galaxies ones. Our PCA-UFP method of classification is robust, flexible and capable of identifying the two main populations of galaxies as well as the intermediate population. The intermediate galaxy population shows many of the properties of the green valley galaxies, and constitutes a more coherent and homogeneous population. The fairly large redshift range of the studied sample allows us to behold the downsizing effect: galaxies with masses of the order of $3\cdot 10^{10}$ Msun mainly are found in transition from the late type to the early type group at $z>0.5$, while galaxies with lower masses - of the order of $10^{10}$ Msun - are in transition at later epochs; galaxies with $M <10^{10}$ Msun did not begin their transition yet, while galaxies with very large masses ($M > 5\cdot 10^{10}$ Msun) mostly completed their transition before $z\sim 1$.

The contribution of major mergers to galaxy mass assembly along cosmic time is an important ingredient to the galaxy evolution scenario. We aim to measure the evolution of the merger rate for both luminosity/mass selected galaxy samples and investigate its dependence with the local environment. We use a sample of 10644 spectroscopically observed galaxies from the zCOSMOS redshift survey to identify pairs of galaxies destined to merge, using only pairs for which the velocity difference and projected separation of both components with a confirmed spectroscopic redshift indicate a high probability of merging. We have identified 263 spectroscopically confirmed pairs with r_p^{max} = 100 h^{-1} kpc. We find that the density of mergers depends on luminosity/mass, being higher for fainter/less massive galaxies, while the number of mergers a galaxy will experience does not depends significantly on its intrinsic luminosity but rather on its stellar mass. We find that the pair fraction and merger rate increase with local galaxy density, a property observed up to redshift z=1. We find that the dependence of the merger rate on the luminosity or mass of galaxies is already present up to redshifts z=1, and that the evolution of the volumetric merger rate of bright (massive) galaxies is relatively flat with redshift with a mean value of 3*10^{-4} (8*10^{-5} respectively) mergers h^3 Mpc^{-3} Gyr^{-1}. The dependence of the merger rate with environment indicates that dense environments favors major merger events as can be expected from the hierarchical scenario. The environment therefore has a direct impact in shapping-up the mass function and its evolution therefore plays an important role on the mass growth of galaxies along cosmic time.

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 <= z <= 10 at the per cent level. To address the problem of biasing between the matter field and the haloes identified in the simulation, we employ the Eulerian local bias model (ELB), including non-linear bias up to the third order. The bias parameters are obtained by fitting a scatter plot of halo and matter density (both from the simulation and from SPT). Using these bias parameters, we can reconstruct the halo density field. We find that this reconstruction is not able to capture all the details of the halo distribution. We investigate how well the large scale bias can be described by a constant and if it corresponds to the linear bias parameter b_1 of the local bias model. We also discuss how well the halo-halo power spectrum and the halo-mass cross spectrum from the reconstructed halo density field agree with the corresponding statistics from the simulation. The results show that while SPT is an excellent approximation for the matter field for suitably large smoothing scales even at redshift 0, the ELB model can only account for some of the properties of the halo density field.

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<10^{14} Msol/h our results at z=0 are in good agreement with N-body simulations. Our model is able to predict the halo-matter cross spectrum with an accuracy of 5 per cent up to k = 0.1 h/Mpc approaching the requirements of future galaxy redshift surveys.

We measure the angular two-point correlation function of sub-millimeter galaxies (SMGs) from 1.1-millimeter imaging of the COSMOS field with the AzTEC camera and ASTE 10-meter telescope. These data yields one of the largest contiguous samples of SMGs to date, covering an area of 0.72 degrees^2 down to a 1.26 mJy/beam (1-sigma) limit, including 189 (328) sources with S/N greater than 3.5 (3). We can only set upper limits to the correlation length r_0, modeling the correlation function as a power-law with pre-assigned slope. Assuming existing redshift distributions, we derive 68.3% confidence level upper limits of r_0 < 6-8 h^-1 Mpc at 3.7 mJy, and r_0 < 11-12 h^-1 Mpc at 4.2 mJy. Although consistent with most previous estimates, these upper limits imply that the real r_0 is likely smaller. This casts doubts on the robustness of claims that SMGs are characterized by significantly stronger spatial clustering, (and thus larger mass), than differently selected galaxies at high-redshift. Using Monte Carlo simulations we show that even strongly clustered distributions of galaxies can appear unclustered when sampled with limited sensitivity and coarse angular resolution common to current sub-millimeter surveys. The simulations, however, also show that unclustered distributions can appear strongly clustered under these circumstances. From the simulations, we predict that at our survey depth, a mapped area of two degrees^2 is needed to reconstruct the correlation function, assuming smaller beam sizes of future surveys (e.g. the Large Millimeter Telescope's 6" beam size). At present, robust measures of the clustering strength of bright SMGs appear to be below the reach of most observations.

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

We present a new three-dimensional radiative transfer (RT) code, RADAMESH, based on a ray-tracing, photon-conserving and adaptive (in space and time) scheme. RADAMESH uses a novel Monte Carlo approach to sample the radiation field within the computational domain on a "cell-by-cell" basis. Thanks to this algorithm, the computational efforts are now focused where actually needed, i.e. within the Ionization-fronts (I-fronts). This results in an increased accuracy level and, at the same time, a huge gain in computational speed with respect to a "classical" Monte Carlo RT, especially when combined with an Adaptive Mesh Refinement (AMR) scheme. Among several new features, RADAMESH is able to adaptively refine the computational mesh in correspondence of the I-fronts, allowing to fully resolve them within large, cosmological boxes. We follow the propagation of ionizing radiation from an arbitrary number of sources and from the recombination radiation produced by H and He. The chemical state of six species (HI, HII, HeI, HeII, HeIII, e) and gas temperatures are computed with a time-dependent, non-equilibrium chemistry solver. We present several validating tests of the code, including the standard tests from the RT Code Comparison Project and a new set of tests aimed at substantiating the new characteristics of RADAMESH. Using our AMR scheme, we show that properly resolving the I-front of a bright quasar during Reionization produces a large increase of the predicted gas temperature within the whole HII region. Also, we discuss how H and He recombination radiation is able to substantially change the ionization state of both species (for the classical Stroemgren sphere test) with respect to the widely used "on-the-spot" approximation.

We study the impact of the environment on the evolution of galaxies in the zCOSMOS 10k sample in the redshift range 0.1<z<1.0 over an area of ~1.5 deg2. The considered sample of secure spectroscopic redshifts contains about 8500 galaxies, with their stellar masses estimated by SED fitting of the multiwavelength optical to NIR photometry. The evolution of the galaxy stellar mass function (GSMF) in high and low density regions provides a tool to study the mass assembly evolution in different environments; moreover, the contributions to the GSMF from different galaxy types, as defined by their SEDs and their morphologies, can be quantified. At redshift z~1, the GSMF is only slightly dependent on environment, but at lower redshifts the shapes of the GSMFs in high- and low-density environments become extremely different, with high density regions exhibiting a marked bimodality. As a result, we infer that galaxy evolution depends on both the stellar mass and the environment, the latter setting the probability of a galaxy to have a given mass: all the galaxy properties related to the stellar mass show a dependence on environment, reflecting the difference observed in the mass functions. The shapes of the GSMFs of early- and late-type galaxies are almost identical for the extremes of the density contrast we consider. The evolution toward z=0 of the mass at which the early- and late-type GSMFs match each other is more rapid in high density environments. The comparison of the observed GSMFs to the same quantities derived from a set of mock catalogues shows that blue galaxies in sparse environments are overproduced in the semi-analytical models at intermediate and high masses, because of a deficit of star formation suppression, while at z<0.5 an excess of red galaxies is present in dense environments at intermediate and low masses, because of the overquenching of satellites. ABRIDGED

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.

We present the analysis of the U-V rest-frame color distribution and some spectral features as a function of mass and environment for two sample of early-type galaxies up to z=1 extracted from the zCOSMOS spectroscopic survey. The first sample ("red galaxies") is defined with a photometric classification, while the second ("ETGs") by combining morphological, photometric, and spectroscopic properties to obtain a more reliable sample. We find that the color distribution of red galaxies is not strongly dependent on environment for all mass bins, with galaxies in overdense regions redder than galaxies in underdense regions with a difference of 0.027\pm0.008 mag. The dependence on mass is far more significant, with average colors of massive galaxies redder by 0.093\pm0.007 mag than low-mass galaxies throughout the entire redshift range. We study the color-mass relation, finding a mean slope 0.12\pm0.005, while the color-environment relation is flatter, with a slope always smaller than 0.04. The spectral analysis that we perform on our ETGs sample is in good agreement with our photometric results: we find for D4000 a dependence on mass between high and low-mass galaxies, and a much weaker dependence on environment (respectively a difference of of 0.11\pm0.02 and of 0.05\pm0.02); for the equivalent width of H{\delta}we measure a difference of 0.28\pm0.08 {\AA}across the same mass range and no significant dependence on environment.By analyzing the lookback time of early-type galaxies, we support the possibility of a downsizing scenario, in which massive galaxies with a stronger D4000 and an almost constant equivalent width of $H\delta$ formed their mass at higher redshift than lower mass ones. We also conclude that the main driver of galaxy evolution is the galaxy mass, the environment playing a subdominant role.

[Abridged] With the first 10000 spectra of the flux limited zCOSMOS sample (I<=22.5) we study the evolution of environmental effects on galaxy properties since z=1.0, and disentangle the dependence among galaxy colour, stellar mass and local density (3D local density contrast `delta', computed with the 5th nearest neighbour approach). We confirm that within a luminosity-limited sample (M_B<=-20.5-z) the fraction of red (U-B>=1) galaxies 'f_red' depends on delta at least up to z=1, with red galaxies residing mainly in high densities. This trend weakens for increasing z, and it is mirrored by the behaviour of the fraction of galaxies with D4000A break >=1.4. We also find that up to z=1 the fraction of galaxies with log(EW[OII]) >=1.15 is higher for lower delta, and also this dependence weakens for increasing z. Given the triple dependence among galaxy colours, stellar mass and delta, the colour-delta relation found in the luminosity-selected sample can be due to the broad range of stellar masses. Thus, we fix the stellar mass and we find that in this case the colour-delta relation is flat up to z=1 for galaxies with log(M/M_sun)>=10.7. This means that for these masses the colour-delta relation found in a luminosity-selected sample is the result of the combined colour-mass and mass-delta relations. In contrast, we find that for 0.1<=z<=0.5 and log(M/M_sun)<=10.7 'f_red' depends on delta even at fixed mass. In these mass and z ranges, environment affects directly also galaxy colours. We suggest a scenario in which the colour depends primarily on stellar mass, but for relatively low mass galaxies the local density modulates this dependence. These galaxies formed more recently, in an epoch when evolved structures were already in place, and their longer SFH allowed environment-driven physical processes to operate during longer periods of time.

We present a galaxy group catalogue spanning the redshift range 0.1 <~ z <~ 1 in the ~1.7 deg^2 COSMOS field, based on the first ~10,000 zCOSMOS spectra. The performance of both the Friends-of-Friends (FOF) and Voronoi-Delaunay-Method (VDM) approaches to group identification has been extensively explored and compared using realistic mock catalogues. We find that the performance improves substantially if groups are found by progressively optimizing the group-finding parameters for successively smaller groups, and that the highest fidelity catalogue, in terms of completeness and purity, is obtained by combining the independently created FOF and VDM catalogues. The final completeness and purity of this catalogue, both in terms of the groups and of individual members, compares favorably with recent results in the literature. The current group catalogue contains 102 groups with N >= 5 spectroscopically confirmed members, with a further ~700 groups with 2 <= N <= 4. Most of the groups can be assigned a velocity dispersion and a dark-matter mass derived from the mock catalogues, with quantifiable uncertainties. The fraction of zCOSMOS galaxies in groups is about 25% at low redshift and decreases toward ~15% at z ~ 0.8. The zCOSMOS group catalogue is broadly consistent with that expected from the semi-analytic evolution model underlying the mock catalogues. Not least, we show that the number density of groups with a given intrinsic richness increases from redshift z ~ 0.8 to the present, consistent with the hierarchical growth of structure.

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 present the Galaxy Stellar Mass Function (MF) up to z~1 from the zCOSMOS-bright 10k spectroscopic sample. We investigate the total MF and the contribution of ETGs and LTGs, defined by different criteria (SED, morphology or star formation). We unveil a galaxy bimodality in the global MF, better represented by 2 Schechter functions dominated by ETGs and LTGs, respectively. For the global population we confirm that low-mass galaxies number density increases later and faster than for massive galaxies. We find that the MF evolution at intermediate-low values of Mstar (logM<10.6) is mostly explained by the growth in stellar mass driven by smoothly decreasing star formation activities. The low residual evolution is consistent with ~0.16 merger per galaxy per Gyr (of which fewer than 0.1 are major). We find that ETGs increase in number density with cosmic time faster for decreasing Mstar, with a median "building redshift" increasing with mass, in contrast with hierarchical models. For LTGs we find that the number density of blue or spiral galaxies remains almost constant from z~1. Instead, the most extreme population of active star forming galaxies is rapidly decreasing in number density. We suggest a transformation from blue active spirals of intermediate mass into blue quiescent and successively (1-2 Gyr after) into red passive types. The complete morphological transformation into red spheroidals, required longer time-scales or follows after 1-2 Gyr. A continuous replacement of blue galaxies is expected by low-mass active spirals growing in stellar mass. We estimate that on average ~25% of blue galaxies is transforming into red per Gyr for logM<11. We conclude that the build-up of galaxies and ETGs follows the same downsizing trend with mass as the formation of their stars, converse to the trend predicted by current SAMs. We expect a negligible evolution of the global Galaxy Baryonic MF.

We study structure formation in the presence of primordial non-Gaussianity of the local type with parameters f_NL and g_NL. We show that the distribution of dark-matter halos is naturally described by a multivariate bias scheme where the halo overdensity depends not only on the underlying matter density fluctuation delta, but also on the Gaussian part of the primordial gravitational potential phi. This corresponds to a non-local bias scheme in terms of delta only. We derive the coefficients of the bias expansion as a function of the halo mass by applying the peak-background split to common parametrizations for the halo mass function in the non-Gaussian scenario. We then compute the halo power spectrum and halo-matter cross spectrum in the framework of Eulerian perturbation theory up to third order. Comparing our results against N-body simulations, we find that our model accurately describes the numerical data for wavenumbers k < 0.1-0.3 h/Mpc depending on redshift and halo mass. In our multivariate approach, perturbations in the halo counts trace phi on large scales and this explains why the halo and matter power spectra show different asymptotic trends for k -> 0. This strongly scale-dependent bias originates from terms at leading order in our expansion. This is different from what happens using the standard univariate local bias where the scale-dependent terms come from badly behaved higher-order corrections. On the other hand, our biasing scheme reduces to the usual local bias on smaller scales where |phi| is typically much smaller than the density perturbations. We finally discuss the halo bispectrum in the context of multivariate biasing and show that, due to its strong scale and shape dependence, it is a powerful tool for the detection of primordial non-Gaussianity from future galaxy surveys.

The identities of the main processes triggering and quenching star-formation in galaxies remain unclear. A key stage in evolution, however, appears to be represented by post-starburst galaxies. To investigate their impact on galaxy evolution, we initiated a multiwavelength study of galaxies with k+a spectral features in the COSMOS field. We examine a mass-selected sample of k+a galaxies at z=0.48-1.2 using the spectroscopic zCOSMOS sample. K+a galaxies occupy the brightest tail of the luminosity distribution. They are as massive as quiescent galaxies and populate the green valley in the colour versus luminosity (or stellar mass) distribution. A small percentage (<8%) of these galaxies have radio and/or X-ray counterparts (implying an upper limit to the SFR of ~8Msun/yr). Over the entire redshift range explored, the class of k+a galaxies is morphologically a heterogeneous population with a similar incidence of bulge-dominated and disky galaxies. This distribution does not vary with the strength of the Hdelta absorption line but instead with stellar mass in a way reminiscent of the well-known mass-morphology relation. Although k+a galaxies are also found in underdense regions, they appear to reside typically in a similarly rich environment as quiescent galaxies on a physical scale of ~2-8Mpc, and in groups they show a morphological early-to-late type ratio similar to the quiescent galaxy class. With the current data set, we do not find evidence of statistical significant evolution in either the number/mass density of k+a galaxies at intermediate redshift with respect to the local values, or the spectral properties. Those galaxies, which are affected by a sudden quenching of their star-formation activity, may increase the stellar mass of the red-sequence by up to a non-negligible level of ~10%.

We present a catalog of 213 type-2 AGN selected from the zCOSMOS survey. The selected sample covers a wide redshift range (0.15<z<0.92) and is deeper than any other previous study, encompassing the luminosity range 10^{5.5} < Lsun< L[OIII] < 10^{9.1} Lsun. We explore the intrinsic properties of these AGN and the relation to their X-ray emission (derived from the XMM-COSMOS observations). We study their evolution by computing the [OIII]5007A line luminosity function (LF) and we constrain the fraction of obscured AGN as a function of luminosity and redshift. The sample was selected on the basis of the optical emission line ratios, after applying a cut to the signal-to-noise ratio (S/N) of the relevant lines. We used the standard diagnostic diagrams [OIII]/Hbeta versus [NII]/Halpha and ([OIII]/Hbeta versus [SII]/Halpha) to isolate AGN in the redshift range 0.15<z<0.45 and the diagnostic diagram [OIII]/Hbeta versus [OII]/Hbeta to extend the selection to higher redshift (0.5<z<0.92). Combining our sample with one drawn from SDSS, we found that the best description of the evolution of type-2 AGN is a luminosity-dependent density evolution model. Moreover, using the type-1 AGN LF we were able to constrain the fraction of type-2 AGN to the total (type-1 + type-2) AGN population. We found that the type-2 fraction decreases with luminosity, in agreement with the most recent results, and shows signs of a slight increase with redshift. However, the trend with luminosity is visible only after combining the SDSS+zCOSMOS samples. From the COSMOS data points alone, the type-2 fraction seems to be quite constant with luminosity.

We investigate the properties and the environment of radio sources with optical counterpart from the combined VLA-COSMOS and zCOSMOS samples. The advantage of this sample is the availability of optical spectroscopic information, high quality redshifts, and accurate density determination. By comparing the star formation rates estimated from the optical spectral energy distribution with those based on the radio luminosity, we divide the radio sources in three families, passive AGN, non-passive AGN and star forming galaxies. These families occupy specific regions of the 8.0-4.5 $\mu$m infrared color--specific star formation plane, from which we extract the corresponding control samples. Only the passive AGN have a significantly different environment distribution from their control sample. The fraction of radio-loud passive AGN increases from ~2% in underdense regions to ~15% for overdensities (1+delta) greater than 10. This trend is also present as a function of richness of the groups hosting the radio sources. Passive AGN in overdensities tend to have higher radio luminosities than those in lower density environments. Since the black hole mass distribution is similar in both environments, we speculate that, for low radio luminosities, the radio emission is controlled (through fuel disponibility or confinement of radio jet by local gas pressure) by the interstellar medium of the host galaxy, while in other cases it is determined by the structure (group or cluster) in which the galaxy resides.

We have used the zCOSMOS-bright 10k sample to identify 3244 Spitzer/MIPS 24-micron-selected galaxies with 0.06< S(24um)< 0.50 mJy and I(AB)<22.5, over 1.5 deg^2 of the COSMOS field, and studied different spectral properties, depending on redshift. At 0.2<z<0.3, we found that different reddening laws of common use in the literature explain the dust extinction properties of around 80% of our infrared (IR) sources, within the error bars. For up to 16% of objects, instead, the Halpha/Hbeta ratios are too high for their IR/UV attenuations, which is probably a consequence of inhomogenous dust distributions. In only a few of our galaxies at 0.2<z<0.3 the IR emission could be mainly produced by dust heated by old rather than young stars. Besides, the line ratios of ~22% of our galaxies suggest that they might be star-formation/nuclear-activity composite systems. At 0.5<z<0.7, we estimated galaxy metallicities for 301 galaxies: at least 12% of them are securely below the upper-branch mass-metallicity trend, which is consistent with the local relation. Finally, we performed a combined analysis of the Hdelta equivalent-width versus Dn(4000) diagram for 1722 faint and bright 24um galaxies at 0.6<z<1.0, spanning two decades in mid-IR luminosity. We found that, while secondary bursts of star formation are necessary to explain the position of the most luminous IR galaxies in that diagram, quiescent, exponentially-declining star formation histories can well reproduce the spectral properties of ~40% of the less luminous sources. Our results suggest a transition in the possible modes of star formation at total IR luminosities L(TIR)=(3 +/-2)x10^11 Lsun.

New statistical properties of dark matter halos in Lagrangian space are presented. Tracing back the dark matter particles constituting bound halos resolved in a series of N-body simulations, we measure quantitatively the correlations of the proto-halo's inertia tensors with the local tidal tensors and investigate how the correlation strength depends on the proto-halo's sphericity, local density and filtering scale. It is shown that the majority of the proto-halos exhibit strong correlations between the two tensors provided that the tidal field is smoothed on the proto-halo's mass scale. The correlation strength is found to increase as the proto-halo's sphericity increases, as the proto-halo's mass increases, and as the local density becomes close to the critical value, delta_{ec}. It is also found that those peculiar proto-halos which exhibit exceptionally weak correlations between the two tensors tend to acquire higher specific angular momentum in Eulerian space, which is consistent with the linear tidal torque theory. In the light of our results, it is intriguing to speculate a hypothesis that the low surface brightness galaxies observed at present epoch correspond to the peculiar proto-halos with extreme low-sphericity whose inertia tensors are weakly correlated with the local tidal tensors.

We use the overdensity field reconstructed in the volume of the COSMOS area to study the nonlinear biasing of the zCOSMOS galaxies. The galaxy overdensity field is reconstructed using the current sample of ~8500 accurate zCOSMOS redshifts at I(AB)<22.5 out to z~1 on scales R from 8 to 12 Mpc/h. By comparing the probability distribution function (PDF) of galaxy density contrast delta_g to the lognormal approximation of the PDF of the mass density contrast delta, we obtain the mean biasing function b(delta,z,R) between the galaxy and matter overdensity field and its second moments b(hat) and b(tilde) up to z~1. Over the redshift interval 0.4<z<1 the conditional mean function <delta_g|delta> = b(delta,z,R) delta is of the following characteristic shape. The function vanishes in the most underdense regions and then sharply rises in a nonlinear way towards the mean densities. <delta_g|delta> is almost a linear tracer of the matter in the overdense regions, up to the most overdense regions in which it is nonlinear again and the local effective slope of <delta_g|delta> vs. delta is smaller than unity. The <delta_g|delta> function is evolving only slightly over the redshift interval 0.4<z<1. The linear biasing parameter increases from b(hat)=1.24+/-0.11 at z=0.4 to b(hat)=1.64+/-0.15 at z=1 for the M_B<-20-z sample of galaxies. b(hat) does not show any dependence on the smoothing scale from 8 to 12 Mpc/h, but increases with luminosity. The measured nonlinearity parameter b(tilde)/b(hat) is of the order of a few percent (but it can be consistent with 0) and it does not change with redshift, the smoothing scale or the luminosity. By matching the linear bias of galaxies to the halo bias, we infer that the M_B<-20-z galaxies reside in dark matter haloes with a characteristic mass of about 3-6 x 10^12 Msol, depending on the halo bias fit.

Galaxies on the largest scales of the Universe are observed to be embedded in the filamentary cosmic web which is shaped by the nonlinear tidal field. As an efficient tool to quantitatively describe the statistics of this cosmic web, we present the anisotropic two-point correlation functions of the nonlinear traceless tidal field in the principal-axis frame, which are measured using numerical data from an N-body simulation. We show that both of the nonlinear density and traceless tidal fields are more strongly correlated along the directions perpendicular to the eigenvectors associated with the largest eigenvalues of the local tidal field. The correlation length scale of the traceless tidal field is found to be ~20 Mpc/h, which is much larger than that of the density field ~5 Mpc/h. We also provide analytic fitting formulae for the anisotropic correlation functions of the traceless tidal field, which turn out to be in excellent agreement with the numerical results. We expect that our numerical results and analytic formula are useful to disentangle cosmological information from the filamentary network of the large-scale structures.

(Abridged) We studied the evolution in the B band luminosity function to z~1 in the zCOSMOS 10k sample, for which both accurate galaxy classifications and a detailed description of the local density field are available. The global LF exhibits a brightening of ~0.7 mag in M* from z~0.2 to z~0.9. At low z, late types dominate at faint magnitudes, while the bright end is populated mainly by early types. At higher z, late-type galaxies evolve significantly and, at z~1, the contribution from the various types to the bright end of the LF is comparable. The evolution for early types is in both luminosity and normalization. A similar behaviour is exhibited by late types, but with an opposite trend for the normalization. Studying the role of the environment, we find that the global LF of galaxies in overdense regions has always a brighter M* and a flatter slope. In low density environments, the main contribution to the LF is from blue galaxies, while for high density environments there is an important contribution from red galaxies to the bright end. The differences between the global LF in the two environments are not due to only a difference in the relative numbers of red and blue galaxies, but also to their relative luminosity distributions: the value of M* for both types in underdense regions is always fainter than in overdense environments. The "specular" evolution of late- and early-type galaxies is consistent with a scenario where a part of blue galaxies is transformed in red galaxies with increasing cosmic time, without significant changes in the fraction of intermediate-type galaxies. The bulk of this tranformation in overdense regions probably happened before z~1, while it is still ongoing at lower z in underdense environments.

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.

We took advantage of the wealth of information provided by the first ~10000 galaxies of the zCOSMOS-bright survey and its group catalogue to study the complex interplay between group environment and galaxy properties. The classical indicator F_blue (fraction of blue galaxies) proved to be a simple but powerful diagnostic tool. We studied its variation for different luminosity and mass selected galaxy samples. Using rest-frame B-band selected samples, the groups galaxy population exhibits significant blueing as redshift increases, but maintains a lower F_blue with respect both to the global and the isolated galaxy population. However moving to mass selected samples it becomes apparent that such differences are largely due to the biased view imposed by the B-band luminosity selection, being driven by the population of lower mass, bright blue galaxies for which we miss the redder, equally low mass, counterparts. By focusing the analysis on narrow mass bins such that mass segregation becomes negligible we find that only for the lowest mass bin explored (logMass <= 10.6) does a significant residual difference in color remain as a function of environment, while this difference becomes negligible toward higher masses. Our results indicate that red galaxies of logMass >= 10.8 are already in place at z ~ 1 and do not exhibit any strong environmental dependence, possibly originating from so-called 'nature'/internal mechanisms. In contrast, for lower galaxy masses and redshifts lower than z ~ 1, we observe the emergence in groups of a population of 'nurture' red galaxies: slightly deviating from the trend of the downsizing scenario followed by the global galaxy population, and more so with cosmic time. These galaxies exhibit signatures of group-related secular physical mechanisms directly influencing galaxy evolution.

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<0.05 h Mpc^-1) which was already detected in previous studies. However, we find that commonly used parameterisations based on the peak-background split do not provide an accurate description of our simulations which present extra dependencies on the wavenumber, the non-linearity parameter and, possibly, the clustering strength. We provide an accurate fit of the simulation data that can be used as a benchmark for future determinations of F_NL with galaxy surveys.

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

For more than two decades we have known that galaxy morphological segregation is present in the Local Universe. It is important to see how this relation evolves with cosmic time. To investigate how galaxy assembly took place with cosmic time, we explore the evolution of the morphology-density relation up to redshift z~1 using about 10000 galaxies drawn from the zCOSMOS Galaxy Redshift Survey. Taking advantage of accurate HST/ACS morphologies from the COSMOS survey, of the well-characterised zCOSMOS 3D environment, and of a large sample of galaxies with spectroscopic redshift, we want to study here the evolution of the morphology-density relation up to z~1 and its dependence on galaxy luminosity and stellar mass. The multi-wavelength coverage of the field also allows a first study of the galaxy morphological segregation dependence on colour. We further attempt to disentangle between processes that occurred early in the history of the Universe or late in the life of galaxies. The zCOSMOS field benefits of high-resolution imaging in the F814W filter from the Advanced Camera for Survey (ACS). We use standard morphology classifiers, optimised for being robust against band-shifting and surface brightness dimming, and a new, objective, and automated method to convert morphological parameters into early, spiral, and irregular types. We use about 10000 galaxies down to I_AB=22.5 with a spectroscopic sampling rate of 33% to characterise the environment of galaxies up to z~1 from the 100 kpc scales of galaxy groups up to the 100 Mpc scales of the cosmic web. ABRIDGED

We explore a possible origin for the puzzling anti-correlation between the formation epoch of galactic dark-matter haloes and their environment density. This correlation has been revealed from cosmological N-body simulations and is in conflict with the Extended Press-Schechter model of halo clustering. Using similar simulations, we first quantify the straightforward association of an early formation epoch with a reduced mass growth rate at late times. We then find that a primary driver of suppressed growth, by accretion and mergers, is tidal effects dominated by a neighbouring massive halo. The tidal effects range from a slowdown of the assembly of haloes due to the shear along the large-scale filaments that feed the massive halo to actual mass loss in haloes that pass through the massive halo. Using the restricted three-body problem, we show that haloes are prone to tidal mass loss within 1.5 virial radii of a larger halo. Our results suggest that the dependence of formation epoch on environment density is a secondary effect induced by the enhanced density of haloes in filaments near massive haloes where the tides are strong. Our measures of assembly rate are particularly correlated with the tidal field at high redshifts z~1.

We study the dependence of galaxy clustering on luminosity and stellar mass at redshifts z ~ [0.2-1] using the first zCOSMOS 10K sample. We measure the redshift-space correlation functions xi(rp,pi) and its projection wp(rp) for sub-samples covering different luminosity, mass and redshift ranges. We quantify in detail the observational selection biases and we check our covariance and error estimate techniques using ensembles of semi-analytic mock catalogues. We finally compare our measurements to the cosmological model predictions from the mock surveys. At odds with other measurements, we find a weak dependence of galaxy clustering on luminosity in all redshift bins explored. A mild dependence on stellar mass is instead observed. At z~0.7, wp(rp) shows strong excess power on large scales. We interpret this as produced by large-scale structure dominating the survey volume and extending preferentially in direction perpendicular to the line-of-sight. We do not see any significant evolution with redshift of the amplitude of clustering for bright and/or massive galaxies. The clustering measured in the zCOSMOS data at 0.5<z<1 for galaxies with log(M/M_\odot)>=10 is only marginally consistent with predictions from the mock surveys. On scales larger than ~2 h^-1 Mpc, the observed clustering amplitude is compatible only with ~1% of the mocks. Thus, if the power spectrum of matter is LCDM with standard normalization and the bias has no unnatural scale-dependence, this result indicates that COSMOS has picked up a particularly rare, ~2-3 sigma positive fluctuation in a volume of ~10^6 h^-1 Mpc^3. These findings underline the need for larger surveys of the z~1 Universe to appropriately characterize the level of structure at this epoch.

We use the current sample of ~10,000 zCOSMOS spectra of sources selected with I(AB) < 22.5 to define the density field out to z~1, with much greater resolution in the radial dimension than has been possible with either photometric redshifts or weak lensing. We apply new algorithms that we have developed (ZADE) to incorporate objects not yet observed spectroscopically by modifying their photometric redshift probability distributions using the spectroscopic redshifts of nearby galaxies. This strategy allows us to probe a broader range of galaxy environments and reduce the Poisson noise in the density field. The reconstructed overdensity field of the 10k zCOSMOS galaxies consists of cluster-like patterns surrounded by void-like regions, extending up to z~1. Some of these structures are very large, spanning the ~50 Mpc/h transverse direction of the COSMOS field and extending up to Delta z~0.05 in redshift. We present the three dimensional overdensity maps and compare the reconstructed overdensity field to the independently identified virialised groups of galaxies and clusters detected in the visible and in X-rays. The distribution of the overdense structures is in general well traced by these virialised structures. A comparison of the large scale structures in the zCOSMOS data and in the mock catalogues reveals an excellent agreement between the fractions of the volume enclosed in structures of all sizes above a given overdensity between the data and the mocks in 0.2<z<1.

The impact of environment on AGN activity up to z~1 is assessed by utilizing a mass-selected sample of galaxies from the 10k catalog of the zCOSMOS spectroscopic redshift survey. We identify 147 AGN by their X-ray emission as detected by XMM-Newton from a parent sample of 7234 galaxies. We measure the fraction of galaxies with stellar mass M_*>2.5x10^10 Msun that host an AGN as a function of local overdensity using the 5th, 10th and 20th nearest neighbors that cover a range of physical scales (~1-4 Mpc). Overall, we find that AGNs prefer to reside in environments equivalent to massive galaxies with substantial levels of star formation. Specifically, AGNs with host masses between 0.25-1x10^11 Msun span the full range of environments (i.e., field-to-group) exhibited by galaxies of the same mass and rest-frame color or specific star formation rate. Host galaxies having M_*>10^11 Msun clearly illustrate the association with star formation since they are predominantly bluer than the underlying galaxy population and exhibit a preference for lower density regions analogous to SDSS studies of narrow-line AGN. To probe the environment on smaller physical scales, we determine the fraction of galaxies (M_*>2.5x10^10 Msun) hosting AGNs inside optically-selected groups, and find no significant difference with field galaxies. We interpret our results as evidence that AGN activity requires a sufficient fuel supply; the probability of a massive galaxy to have retained some sufficient amount of gas, as evidence by its ongoing star formation, is higher in underdense regions where disruptive processes (i.e., galaxy harrassment, tidal stripping) are lessened.

We study the spatial clustering of 538 X-ray selected AGN in the 2 deg^2 XMM-COSMOS field that are spectroscopically identified to I_{AB}<23 and span the redshift range z=0.2-3.0. The median redshift and luminosity of the sample are z = 0.98 and L_{0.5-10}=6.3 x 10^{43} erg/s, respectively. A strong clustering signal is detected at ~18sigma level, which is the most significant measurement obtained to date for clustering of X-ray selected AGN. By fitting the projected correlation function w(r_p) with a power law on scales of r_p=0.3-40 Mpc/h, we derive a best fit comoving correlation length of r_0 = 8.6 +- 0.5 Mpc/h and slope of gamma=1.88 +- 0.07 (Poissonian errors; bootstrap errors are about a factor of 2 larger). An excess signal is observed in the range r_p~5-15 Mpc/h, which is due to a large scale structure at z ~ 0.36 containing about 40 AGN. When removing the z ~ 0.36 structure, or computing w(r_p) in a narrower range around the peak of the redshift distribution (e.g. z=0.4-1.6), the correlation length decreases to r_0 ~ 5-6 Mpc/h, which is consistent with that observed for bright optical QSOs at the same redshift. We investigate the clustering properties of obscured and unobscured AGN separately. Within the statistical uncertainties, we do not find evidence that AGN with broad optical lines (BLAGN) cluster differently from AGN without broad optical lines (non-BLAGN). The correlation length measured for XMM-COSMOS AGN at z~1 is similar to that of massive galaxies (stellar mass M_*> 3 x 10^{10} M_sun) at the same redshift. This suggests that AGN at z~1 are preferentially hosted by massive galaxies, as observed both in the local and in the distant (z~2) Universe. (shortened)

We present a study of the host galaxies of AGN selected from the zCOSMOS survey to establish if accretion onto supermassive black holes and star formation are explicitly linked up to z~1. We identify 152 galaxies that harbor AGN, based on XMM observations of 7543 galaxies (i<22.5). Star formation rates (SFRs), including those weighted by stellar mass, are determined using the [OII]3727 emission-line, corrected for an AGN contribution. We find that the majority of AGN hosts have significant levels of star formation with a distribution spanning ~1-100 Msun yr^-1. The close association between AGN activity and star formation is further substantiated by an increase in the AGN fraction with the youthfulness of their stars as indicated by the rest-frame color (U-V) and spectral index Dn(4000); we demonstrate that mass-selection alleviates an artifical peak falling in the transition region when using luminosity-limited samples. We also find that the SFRs of AGN hosts evolve with cosmic time in a manner that closely mirrors the overall galaxy population and explains the low SFRs in AGNs (z<0.3) from the SDSS. We conclude that the conditions most conducive for AGN activity are a massive host galaxy and a large reservoir of gas. Furthermore, a direct correlation between mass accretion rate onto SMBHs and SFR is shown to be weak although the average ratio is constant with redshift, effectively shifting the evidence for a co-evolution scenario in a statistical manner to smaller physical scales. Our findings illustrate an intermittent scenario with an AGN lifetime substantially shorter than that of star formation and underlying complexities regarding fueling over vastly different physical scales yet to be determined [Abridged].

Aims. We investigate the relationships between three main optical galaxy observables (spectral properties, colours, and morphology), exploiting the data set provided by the COSMOS/zCOSMOS survey. The purpose of this paper is to define a simple galaxy classification cube, using a carefully selected sample of around 1000 galaxies. Methods. Using medium resolution spectra of the first 1k zCOSMOS-bright sample, optical photometry from the Subaru/COSMOS observations, and morphological measurements derived from ACS imaging, we analyze the properties of the galaxy population out to z~1. Applying three straightforward classification schemes (spectral, photometric, and morphological), we identify two main galaxy types, which appear to be linked to the bimodality of galaxy population. The three parametric classifications constitute the axes of a "classification cube". Results. A very good agreement exists between the classification from spectral data (quiescent/star-forming galaxies) and that based on colours (red/blue galaxies). The third parameter (morphology) is less well correlated with the first two: in fact a good correlation between the spectral classification and that based on morphological analysis (early-/late-type galaxies) is achieved only after partially complementing the morphological classification with additional colour information. Finally, analyzing the 3D-distribution of all galaxies in the sample, we find that about 85% of the galaxies show a fully concordant classification, being either quiescent, red, bulge-dominated galaxies (~20%) or star-forming, blue, disk-dominated galaxies (~65%). These results imply that the galaxy bimodality is a consistent behaviour both in morphology, colour and dominant stellar population, at least out to z~1.

We investigate the close environment of 203 Spitzer 24 micron-selected sources at 0.6<z<1.0 using zCOSMOS-bright redshifts and spectra of I<22.5 AB mag galaxies, over 1.5 sq. deg. of the COSMOS field. We quantify the degree of passivity of the LIRG and ULIRG environments by analysing the fraction of close neighbours with Dn(4000)>1.4. We find that LIRGs at 0.6<z<0.8 live in more passive environments than those of other optical galaxies that have the same stellar mass distribution. Instead, ULIRGs inhabit more active regions (e.g. LIRGs and ULIRGs at 0.6<z<0.8 have, respectively, (42.0 +/- 4.9)% and (24.5 +/- 5.9)% of neighbours with Dn (4000)>1.4 within 1 Mpc and +/- 500 km/s). The contrast between the activities of the close environments of LIRGs and ULIRGs appears especially enhanced in the COSMOS field density peak at z~0.67, because LIRGs on this peak have a larger fraction of passive neighbours, while ULIRGs have as active close environments as those outside the large-scale structure. The differential environmental activity is related to the differences in the distributions of stellar mass ratios between LIRGs/ULIRGs and their close neighbours, as well as in the general local density fields. At 0.8<z<1.0, instead, we find no differences in the environment densities of ULIRGs and other similarly massive galaxies, in spite of the differential activities. We discuss a possible scenario to explain these findings.

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.

We measure the clustering of a sample of photometrically selected luminous red galaxies around a low redshift (0.2<z<0.6) sample of quasars selected from the Sloan Digital Sky Survey Data Release 5. We make use of a new statistical estimator to obtain precise measurements of the LRG auto-correlations and constrain halo occupation distributions for them. These are used to generate mock catalogs which aid in interpreting our quasar-LRG cross correlation measurements. The cross correlation is well described by a power law with slope 1.8\pm0.1 and r_0=6\pm0.5 h^{-1} Mpc, consistent with observed galaxy correlation functions. We find no evidence for `excess' clustering on 0.1 Mpc scales and demonstrate that this is consistent with the results of Serber et al (2006) and Strand et al (2007), when one accounts for several subtleties in the interpretation of their measurements. Combining the quasar-LRG cross correlation with the LRG auto-correlations, we determine a large-scale quasar bias b_QSO = 1.09\pm0.15 at a median redshift of 0.43, with no observed redshift or luminosity evolution. This corresponds to a mean halo mass <M>~ 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.

ZEBRA, the Zurich Extragalactic Bayesian Redshift Analyzer, is a tool for estimating redshifts and template types of galaxies using medium- and broad-band photometric data. ZEBRA employs novel techniques within the template-fitting approach to produce high-quality Maximum-Likelihood and Bayesian redshift estimates. This manuscript serves as a user guide to ZEBRA. It explains how to use ZEBRA, specifies input and output formats, and gives a short account of the available options. ZEBRA is a free and open-source software distributed under the GNU Public License 3 and available at http://www.exp-astro.phys.ethz.ch/ZEBRA . Several upgrades are currently being implemented in ZEBRA. An updated documentation will be provided at each new release. Any problems, comments and suggestions on the code and the manual should be sent via e-mail to zebra@phys.ethz.ch.

We present a new method to directly map the neutral-hydrogen distribution during the reionization epoch and to constrain the emission properties of the highest-redshift quasars (QSOs). As a tracer of HI, we propose to use the Ly-alpha radiation produced by quasar ionization fronts (I-fronts) that expand in the partially ionized intergalactic medium (IGM) before reionization is complete. These Ly-alpha photons are mainly generated by collisional excitations of hydrogen atoms in the boundary of the rapidly expanding HII region. The observable signal is produced by the part of the I-front that lies behind the QSO with respect to the observer. Combining two radiative transfer models (one for the QSO ionizing radiation and one for the Ly-alpha photons), we estimate the expected Ly-alpha spectral shape and surface brightness (SB) for a large number of configurations where we varied both the properties of the ionizing QSO and of the surrounding medium. We find that the expected signal is observable as a single (broad) emission line with a characteristic width of 100-200 km/s. The expected SB produced at redshift z~6.5 within a fully neutral region (at mean density) by a typical QSO I-front lies in the range $10^{-21}-10^{-20}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ and decreases proportionally to $(1+z)^2$ for a given QSO age. QSOs with harder spectra may produce a significantly brighter emission at early phases. The signal may cover up to a few hundred square arcmin on the sky and should be already detectable with current facilities by means of moderate/high resolution spectroscopy. The detection of this Ly-alpha emission can shed new light on the reionization history, the age and the emission properties of the highest-redshift QSOs. (abridged)

We use a series of N-body simulations of the LCDM cosmology to investigate the redshift evolution since z=1 of the properties and alignment with the large-scale structure of haloes in clusters, filaments, sheets and voids. We find that: (i) Once a rescaling of the halo mass with M*, the mass scale collapsing at redshift z, is performed, there is no further significant redshift dependence in the halo properties; (ii) The environment influences halo shape and formation time at all investigated redshifts for haloes with masses M<M*; and (iii) There is a significant alignment of both spin and shape of haloes with filaments and sheets. In detail, at all redshifts up to z=1: a) Haloes with M<M* tend to be more oblate when located in clusters than in the other environments; this trend is reversed at higher masses: above about M*, halos in clusters are typically more prolate than similar massive haloes in sheets, filaments and voids. b) Haloes with M>M* in filaments spin more rapidly than similar mass haloes in clusters; haloes in voids have the lowest median spin parameters; c) Haloes with M<M* tend to be younger in voids and older in clusters; d) In sheets, halo spin vectors tend to lie within the sheet plane independent of mass; in filaments, instead, haloes with M<M* tend to spin parallel to the filament and haloes with M>M* perpendicular to it. For masses M>M*, the major axis of haloes in filaments and sheets is strongly aligned with the filament or the sheet. Such halo-LSS alignments may be of importance in weak lensing analyses of cosmic shear. A question that is opened by our study is why, in the 0 < z < 1 redshift regime that we have investigated, the mass scale M* sets roughly the threshold below which the LSS-environment either begins to affect, or reverses, the properties of dark matter haloes.

We present measurements of the two-point galaxy angular correlation function w(\theta) in the COSMOS field. Independent determinations of w(\theta) as a function of magnitude limit are presented for both the HST ACS catalog and also for the ground-based data from Subaru and the CFHT. Despite having significantly different masks, these three determinations agree well. At bright magnitudes (IAB<22), our data generally match very well with existing measurements and with mock catalogs based on semi-analytic galaxy formation calculations of Kitzbichler and White from the Millennium Simulation. The exception is that our result is at the upper end of the expected cosmic variance scatter for \theta > 10 arcmin, which we attribute to a particularly rich structure known to exist at z~0.8. For fainter samples, however, the level of clustering is somewhat higher than reported by some previous studies: in all three catalogues we find w(\theta=1')~0.014 at a median IAB magnitude of 24. At these very faintest magnitudes, our measurements agree well with the latest determinations from the Canada-France Legacy Survey. This level of clustering is approximately double what is predicted by the semi-analytic catalogs (at all angles). The semi-analytic results allow an estimate of cosmic variance, which is too small to account for the discrepancy. We therefore conclude that the mean amplitude of clustering at this level is higher than previously estimated.

We present imaging data and photometry for the COSMOS survey in 15 photometric bands between 0.3um and 2.4um. These include data taken on the Subaru 8.3m telescope, the KPNO and CTIO 4m telescopes, and the CFHT 3.6m telescope. Special techniques are used to ensure that the relative photometric calibration is better than 1% across the field of view. The absolute photometric accuracy from standard star measurements is found to be 6%. The absolute calibration is corrected using galaxy spectra, providing colors accurate to 2% or better. Stellar and galaxy colors and counts agree well with the expected values. Finally, as the first step in the scientific analysis of these data we construct panchromatic number counts which confirm that both the geometry of the universe and the galaxy population are evolving.

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 report the results of a survey for fluorescent Ly-alpha emission carried out in the field surrounding the z=3.1 quasar QSO0420-388 using the FORS2 instrument on the VLT. We first review the properties expected for fluorescent Ly-alpha emitters, compared with those of other non-fluorescent Ly-alpha emitters. Our observational search detected 13 Ly-alpha sources sparsely sampling a volume of ~14000 comoving Mpc^3 around the quasar. The properties of these in terms of i) the line equivalent width, ii) the line profile and iii) the value of the surface brightness related to the distance from the quasar, all suggest that several of these may be plausibly fluorescent. Moreover, their number is in good agreement with the expectation from theoretical models. One of the best candidates for fluorescence is sufficiently far behind QSO0420-388 that it would imply that the quasar has been active for (at least) ~60 Myrs. Further studies on such objects will give information about proto-galactic clouds and on the radiative history (and beaming) of the high-redshift quasars.

The incidence of strong MgII systems in gamma-ray burst (GRB) spectra is a few times higher than in quasar (QSO) spectra. We investigate several possible explanations for this effect, including: dust obscuration bias, clustering of the absorbers, different beam sizes of the sources, multiband magnification bias of GRBs, association of the absorbers with the GRB event or the circumburst environment. We find that: i) the incidence rate of MgII systems in QSO spectra could be underestimated by a factor 1.3-2 due to dust obscuration; ii) the equivalent-width distribution of the MgII absorbers along GRBs is consistent with that observed along QSOs thus suggesting that the absorbers are more extended than the beam sizes of the sources; iii) on average, GRB afterglows showing more than one MgII system are a factor of 1.7 brighter than the others, suggesting a lensing origin of the observed discrepancy; iv) gravitational lensing (in different forms, from galaxy lensing to microlensing) can bias high the counts of MgII systems along GRBs if the luminosity functions of the prompt gamma-ray emission and of the optical afterglows have a mean faint-end slope approaching -5/3 -- -2; v) some of the absorbers can be associated with the circumburst environment or produced by supernova remnants unrelated to the GRB event itself but lying in the same star-forming region. With the possible exception of magnification bias, it is unlikely that one of these effects on its own can fully account for the observed counts. However, the combined action of some of them can substantially reduce the statistical significance of the discrepancy.

Using a series of high-resolution N-body simulations of the concordance cosmology we investigate how the formation histories, shapes and angular momenta of dark-matter haloes depend on environment. We first present a classification scheme that allows to distinguish between haloes in clusters, filaments, sheets and voids in the large-scale distribution of matter. This method is based on a local-stability criterion for the orbits of test particles and closely relates to the Zel'dovich approximation. Applying this scheme to our simulations we then find that: i) Mass assembly histories and formation redshifts strongly depend on environment for haloes of mass M<M* (haloes of a given mass tend to be older in clusters and younger in voids) and are independent of it for larger masses; ii) Low-mass haloes in clusters are generally less spherical and more oblate than in other regions; iii) Low-mass haloes in clusters have a higher median spin than in filaments and present a more prominent fraction of rapidly spinning objects; we identify recent major mergers as a likely source of this effect. For all these relations, we provide accurate functional fits as a function of halo mass and environment. We also look for correlations between halo-spin directions and the large-scale structures: the strongest effect is seen in sheets where halo spins tend to lie within the plane of symmetry of the mass distribution. Finally, we measure the spatial auto-correlation of spin directions and the cross-correlation between the directions of intrinsic and orbital angular momenta of neighbouring haloes. While the first quantity is always very small, we find that spin-orbit correlations are rather strong especially for low-mass haloes in clusters and high-mass haloes in filaments.

(ABRIDGED) 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: (i) Five non-parametric diagnostics: asymmetry, concentration, Gini coefficient, 2nd-order moment of the brightest 20% of galaxy pixels, and ellipticity; and (ii) The exponent n of single--Sersic fits to the 2D 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 I<24 COSMOS galaxies to show that the first three PCs fully describe the key aspects of the galaxy structure, i.e., to calibrate a three-dimensional classification grid of axis PC_1, PC_2, and PC_3. We demonstrate the robustness of the ZEST grid on the z=0 sample of Frei et al. (1996). 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 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 magnitudes at z \~0.7). We highlight however two trends, that are in general agreement with a down-sizing scenario for galaxy formation: (1.) A deficit of a factor of about two at z~0.7 of MB>-20.5 structurally--classified early--type galaxies; and (2.) An excess of a factor of about three, at a similar redshift, of irregular galaxies.

We present ZEBRA, the Zurich Extragalactic Bayesian Redshift Analyzer. The current version of ZEBRA combines and extends several of the classical approaches to produce accurate photometric redshifts down to faint magnitudes. In particular, ZEBRA uses the template-fitting approach to produce Maximum Likelihood and Bayesian redshift estimates based on: (1.) An automatic iterative technique to correct the original set of galaxy templates to best represent the SEDs of real galaxies at different redshifts; (2.) A training set of spectroscopic redshifts for a small fraction of the photometric sample; and (3.) An iterative technique for Bayesian redshift estimates, which extracts the full two-dimensional redshift and template probability function for each galaxy. We demonstrate the performance of ZEBRA by applying it to a sample of 866 I_AB <= 22.5 COSMOS galaxies with available u*, B, V, g', r', i', z' and K_s photometry and zCOSMOS spectroscopic redshifts in the range 0 < z < 1.3. Adopting a 5-sigma-clipping that excludes less than 10 galaxies, both the Maximum Likelihood and Bayesian estimates for this sample have an accuracy dz/1+z better than 0.03. Similar accuracies are recovered using mock galaxies.

We study a sample of approximately 16,500 galaxies with I_AB <= 22.5 in the COSMOS field. Structural information on the galaxies is derived by fitting single Sersic models to their two-dimensional surface brightness distributions. We investigate the evolution of the number density of disk galaxies larger than 5 kpc between redshift z~1 and the present epoch. To this end, we use the measurements of the half-light radii to construct, as a function of redshift, the size function of both the total disk galaxy population and of disk galaxies split in four bins of bulge-to-disk ratio. Furthermore, we use a selected sample of roughly 1800 SDSS galaxies to calibrate our results with respect to the local universe. We find that: (i) The number density of disk galaxies with intermediate sizes (r_{1/2}~5-7 kpc) remains nearly constant from z~1 to today. (ii) The number density of the largest disks (r_{1/2}>7 kpc) decreases by a factor of about two out to z~1. (iii) There is a constancy in the number density of large bulgeless disks out to z~1; the deficit of large disks at early epochs seems to arise from a smaller number of bulged disks. Our results indicate that the bulk of the large disk galaxy population has completed its growth by z~1, and support the hypothesis that secular evolution processes produce - or at least add stellar mass to - the bulge components of disk galaxies.

We compute accurate redshift distributions to I(AB) = 24 and R(AB) = 24.5 using photometric redshifts estimated from six-band UBVRIZ photometry in the Canada-France Deep Fields-Photometric Redshift Survey (CFDF-PRS). Our photometric redshift algorithm is calibrated using hundreds of CFRS spectroscopic redshifts in the same fields. The dispersion in redshift is \sigma/(1+z) \la 0.04 to the CFRS depth of I(AB) = 22.5, rising to \sigma/(1+z) \la 0.06 at our nominal magnitude and redshift limits of I(AB) = 24 and z \le 1.3, respectively. We describe a new method to compute N(z) that incorporates the full redshift likelihood functions in a Bayesian iterative analysis and we demonstrate in extensive Monte Carlo simulations that it is superior to distributions calculated using simple maximum likelihood redshifts. The field-to-field differences in the redshift distributions, while not unexpected theoretically, are substantial even on 30' scales. We provide I(AB) and R(AB) redshift distributions, median redshifts, and parametrized fits of our results in various magnitude ranges, accounting for both random and systematic errors in the analysis.

We present detailed clustering measurements for a flux limited sample of 14,000 quasars extracted from the 2dF QSO Redshift Survey (2QZ) in the redshift range 0.8<z<2.1. After splitting the sample into three redshift bins and each of them into six luminosity intervals, we estimate the quasar projected auto and cross-correlation functions at a given redshift for separations 3 < r/(Mpc/h) < 20. Fitting the data with a biased CDM model and using a frequentist analysis (the F-test), we find that models with luminosity dependent clustering are statistically favoured at the 95 per cent confidence level for z>1.3. On the other hand, a number of tests based on information theory and Bayesian statistics show only marginal evidence for luminosity dependent clustering. Anyway, the quality of the data is not good enough to accurately quantify how quasar biasing depends on luminosity. We critically discuss the limitations of our dataset and show that a much larger sample is needed to rule out current models for luminosity segregation. Studying the evolution of the clustering amplitude with redshift, we detect an increase of the quasar correlation length with lookback time at the 99.3 per cent confidence level. Adopting the concordance cosmological model, we discuss the evolution of quasar biasing with cosmic epoch and show that quasars are typically hosted by dark matter haloes with mass ~10^{13} M_sol.

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<z<12 can quantitatively explain the high cross-correlation between CIV systems and LBGs observed at z=3. The reason is twofold. First, both LBGs and high-z dwarfs are biased tracer of the mass distribution and form from exceptionally high density fluctuations which are strongly clustered. Second, the action of gravity tends to increase the spatial association between metal bubbles and LBGs. Our analysis shows that, in order to match the abundance of CIV systems observed at z=3, the metal bubbles generated by high-z dwarfs must have comoving sizes of ~100 kpc. We conclude that the observed galaxy-CIV spatial association needs not to be generated by late ``superwinds'' from LBGs.

We combine a high-resolution hydro-simulation of the LambdaCDM cosmology with two radiative transfer schemes (for continuum and line radiation) to predict the properties, spectra and spatial distribution of fluorescent Ly-alpha emission at z~3. We focus on line radiation produced by recombinations in the dense intergalactic medium ionized by UV photons. In particular, we consider both a uniform background and the case where gas clouds are illuminated by a nearby quasar. We find that the emission from optically thick regions is substantially less than predicted from the widely used static, plane-parallel model. The effects induced by a realistic velocity field and by the complex geometric structure of the emitting regions are discussed in detail. We make predictions for the expected brightness and size distributions of the fluorescent sources.Our results account for recent null detections and can be used to plan new observational campaigns both in the field (to measure the intensity of the diffuse UV background) and in the proximity of bright quasars (to understand the origin of high colum-density absorbers).

We present detailed clustering measurements from the 2dF QSO Redshift Survey (2QZ) in the redshift range 0.8<z<2.1. Using a flux limited sample of ~14,000 objects with effective redshift z_eff=1.47, we estimate the quasar projected correlation function for separations 1<r/h^{-1} Mpc<20. We find that the 2-point correlation function in real space is well approximated by a power law with slope gamma=1.5 pm 0.2 and comoving correlation length r_0=4.8^{+0.9}_{-1.5} h^{-1} Mpc. Splitting the sample into three subsets based on redshift, we find evidence for an increase of the clustering amplitude with lookback time. For a fixed gamma, evolution of r_0 is detected at the 3.6 sigma confidence level. The ratio between the quasar correlation function and the mass autocorrelation function (derived adopting the concordance cosmological model) is found to be scale independent. For a linear mass-clustering amplitude sigma_8=0.8, the ``bias parameter'' decreases from b~3.9$ at z_eff=1.89 to b~1.8 at z_eff=1.06. From the observed abundance and clustering, we infer how quasars populate dark-matter haloes of different masses. We find that 2QZ quasars sit in haloes with M>10^{12} M_sun and that the mean mass of their host haloes is of the order of 10^{13} M_sun. The observed clustering is consistent with assuming that the locally observed correlation between black-hole mass and host-galaxy circular velocity is still valid at z>1. From the fraction of haloes which contain active quasars, we infer that the characteristic quasar lifetime is t_Q ~ a few x 10^7 yr at z~1 and approaches 10^8 yr at higher redshifts.

We study the evolution of galaxy clustering in various cosmological models with quintessence. We investigate how the analytical predictions vary with change of dark energy equation of state $w_X$. Comparing these predictions against available data we discuss to what extent the problems of galaxy biasing can be modelled. This will be key in constraining the dark energy equation of state with future galaxy surveys. We use a compilation of various surveys to study the number density and amplitude of galaxy clustering from observations of the local universe at $z \sim 0$ to that of the Lyman break galaxies and Ly-$\alpha$ emitters at $z \sim 4.9$. We find that there is a degeneracy between the dark energy equation of state and the way galaxies populate dark matter haloes; objects are more biased in models with more negative values of dark energy equation of state $w_X$. We conclude that, while future all sky CMB observations will determine cosmological parameters with unprecedented precision, and cross correlation of weak lensing experiments and galaxy surveys will provide a cleaner and accurate picture of bias associated with collapsed objects, the rate of growth of large scale structure in such surveys can potentially constrain the equation of state of dark energy and the potential of the scalar field associated with quintessence. In particular, we show that the abundance and spatial distribution of galaxy clusters at intermediate redshifts strongly depend on the dark energy equation of state. When accurate measurement of galaxy clustering at high-redshit becomes possible, it will provide constraints on dark energy that are independent and complementary to type Ia supernova studies.

We estimate energy spectra and fluxes at the Earth's surface of the cosmic and Galactic neutrino backgrounds produced by thermonuclear reactions in stars. The extra-galactic component is obtained by combining the most recent estimates of the cosmic star formation history and the stellar initial mass function with accurate theoretical predictions of the neutrino yields all over the thermonuclear lifetime of stars of different masses. Models of the structure and evolution of the Milky Way are used to derive maps of the expected flux generated by Galactic sources as a function of sky direction. The predicted neutrino backgrounds depend only slightly on model parameters. In the relevant 50 keV-10 MeV window, the total flux of cosmic neutrinos ranges between 20 and 65 particles per square cm per s. Neutrinos reaching the Earth today have been typically emitted at redshift z~2. Their energy spectrum peaks at E~0.1-0.3 MeV. The energy and entropy densities of the cosmic background are negligible with respect to the thermal contribution of relic neutrinos originated in the early universe. In every sky direction, the cosmic background is outnumbered by the Galactic one, whose integrated flux amounts to 300-1000 particles per square cm per s. The emission from stars in the Galactic disk contributes more than 95 per cent of the signal.

We use the clustering results obtained by Madgwick et al. (2003) for a sample of 96,791 2dF galaxies with redshift 0.01 \lt z \lt 0.15 to study the distribution of late-type and early-type galaxies within dark matter haloes of different mass. Within the framework of our models, galaxies of both classes are found to be as spatially concentrated as the dark matter within haloes even though, while the distribution of star-forming galaxies can also allow for some steeper profiles, this is drastically ruled out in the case of early-type galaxies. We also find evidence for morphological segregation, as late-type galaxies appear to be distributed within haloes of mass scales corresponding to groups and clusters up to about two virial radii, while passive objects show a preference to reside closer to the halo centre. If we assume a broken power-law of the form <N_{gal}(m)>=(m/m_0)^{alpha_1} for m_{cut}\le m<m_0 and <N_{gal}(m)>=(m/m_0)^{alpha_2} at higher masses to describe the dependence of the average number of galaxies within haloes on the halo mass, fits to the data show that star-forming galaxies start appearing in haloes of masses m_{cut}~10^{11}m_{sun}, much smaller than what is obtained for early-type galaxies (m_{cut}~10^{12.6}m_{sun}). In the high-mass regime m\ge m_0, <N_{gal}>$ increases with halo mass more slowly (alpha_2~0.7) in the case of late-type galaxies than for passive objects which present alpha_2~1.1. We stress that there is no degeneracy in the determination of the best functional forms for rho(r) and <N_{gal}>, as they affect the behaviour of the galaxy-galaxy correlation function on different scales.

We use the clustering results obtained by Madgwick et al. (2003) for a sample of 96,791 2dF galaxies with redshift 0.01 < z < 0.15 to study the distribution of late-type and early-type galaxies within dark matter haloes of different mass. The adopted method relies on the connection between the distribution of sources within haloes and their clustering properties by focussing on the issue of the halo occupation function i.e. the probability distribution of the number of galaxies brighter than some luminosity threshold hosted by a virialized halo of given mass. Within this framework, the distribution of galaxies within haloes is shown to determine galaxy-galaxy clustering on small scales, being responsible for the observed power-law behaviour at separations r < 3 Mpc. For a more extended analysis, we refer the reader to Magliocchetti & Porciani (2003)

We investigate the cross-talk between the two key components of tidal-torque theory, the inertia (I) and shear (T) tensors, using a cosmological N-body simulation with thousands of well-resolved haloes. We find that the principal axes of I and T are strongly aligned, even though I characterizes the proto-halo locally while T is determined by the large-scale structure. Thus, the resultant galactic spin, which plays a key role in galaxy formation, is only a residual due to ~10 per cent deviations from perfect alignment of T and I. The T-I correlation induces a weak tendency for the proto-halo spin to be perpendicular to the major axes of T and I, but this correlation is erased by non-linear effects at late times, making the observed spins poor indicators of the initial shear field. However, the T-I correlation implies that the shear tensor can be used for identifying the positions and boundaries of proto-haloes in cosmological initial conditions -- a missing piece in galaxy formation theory. The typical configuration is of a prolate proto-halo lying perpendicular to a large-scale high-density ridge, with the surrounding voids inducing compression along the major and intermediate inertia axes of the proto-halo. This leads to a transient sub-halo filament along the large-scale ridge, whose sub-clumps then flow along the filament and merge into the final halo. The centres of proto-haloes tend to lie in ~1 sigma over-density regions, but their association with linear density maxima smoothed on galactic scales is vague: only ~40 per cent of the proto-haloes contain peaks within them. Several other characteristics distinguish proto-haloes from density peaks, e.g., they tend to compress along two principal axes while many peaks compress along three axes.

We evaluate the success of linear tidal-torque theory (TTT) in predicting galactic-halo spin using a cosmological N-body simulation with thousands of well-resolved haloes. The proto-haloes are identified by tracing today's haloes back to the initial conditions. The TTT predictions for the proto-haloes match, on average, the spin amplitudes of today's virialized haloes if linear growth is assumed until ~t_0/3, or 55-70 per cent of the halo effective turn-around time. This makes it a useful qualitative tool for understanding certain average properties of galaxies, such as total spin and angular-momentum distribution within haloes, but with a random scatter of the order of the signal itself. Non-linear changes in spin direction cause a mean error of ~50 degrees in the TTT prediction at t_0, such that the linear spatial correlations of spins on scales >~ 1 Mpc/h are significantly weakened by non-linear effects. This questions the usefulness of TTT for predicting intrinsic alignments in the context of gravitational lensing. We find that the standard approximations made in TTT, including a second order expansion of the Zel'dovich potential and a smoothing of the tidal field, provide close-to-optimal results.

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<=25.5 extracted from the survey by Steidel and collaborators, and found the following results. At scales in the range 30<~theta<~100 arcsec, the angular correlation function w(theta) can be accurately described as a power law with slope beta=0.50 +0.25-0.50(1 sigma random)-0.10(systematic), shallower than the measure presented by Giavalisco et al. However, the spatial correlation length, derived by Limber deprojection, is in very good agreement with the previous measures, confirming the strong spatial clustering of these sources. We discuss in detail the effects of both random and systematic errors, in particular of the so called ``integral constraint'' bias, to which we set a lower limit using numerical simulations. This suggests that the current samples do not yet provide a ``fair representation'' of the large-scale distribution of LBGs at z~3. An intriguing result of our analysis is that at angular separations smaller than theta<~30 arcsec the correlation function seems to depart from the power-law fitted at larger scales and become smaller. This feature is detected at the ~90 per cent confidence level and, if real, it can provide information on the number density and spatial distribution of LBGs within their host halos as well as the size and the mass of the halos.

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.

[Abridged] We study the angular-momentum profiles of a statistical sample of halos drawn from a high-resolution N-body simulation of the LCDM cosmology. We find that the cumulative mass distribution of specific angular momentum, j, in a halo of mass Mv is well fit by a universal function, M(<j) = Mv \mu j/(j_0+j). This profile is defined by one shape parameter (\mu or j_0) in addition to the global spin parameter \lambda. It follows a power-law over most of the mass, and flattens at large j, with the flattening more pronounced for small values of \mu. Compared to a uniform sphere in solid-body rotation, most halos have a higher fraction of their mass in the low- and high-j tails of the distribution. The spatial distribution of angular momentum in halos tends to be cylindrical and is well-aligned within each halo for ~80% of the halos. We investigate two ideas for the origin of this profile. The first is based on a revised version of linear tidal-torque theory combined with extended Press-Schechter mass accretion, and the second focuses on j transport in minor mergers. Finally, we briefly explore implications of the M(<j) profile on the formation of galactic disks assuming that j is conserved during an adiabatic baryonic infall. The implied gas density profile deviates from an exponential disk, with a higher density at small radii and a tail extending to large radii. The steep central density profiles may imply disk scale lengths that are smaller than observed. This is reminiscent of the "angular-momentum problem" seen in hydrodynamic simulations, even though we have assumed perfect j conservation. A possible solution is to associate the central excesses with bulge components and the outer regions with extended gaseous disks.

We summarize recent developments in the study of the origin of halo spin profiles and preliminary implications on disk formation. The specific angular-momentum distributions within halos in N-body simulations match a universal shape, M(<j) \propto j/(j_0+j). It is characterized by a power law over most of the mass, and one shape parameter in addition to the spin parameter \lambda. The angular momentum tends to be aligned throughout the halo and of cylindrical symmetry. Even if angular momentum is conserved during baryonic infall, the resultant disk density profile is predicted to deviate from exponential, with a denser core and an extended tail. A slightly corrected version of the scaling relation due to linear tidal-torque theory is used to explain the origin of a typical power-law profile in shells, j(M) \propto M^s with s \gsim 1. While linear theory crudely predicts the amplitudes of halo spins, it is not a good predictor of their directions. Independently, mergers of halos are found to produce a similar profile due to j transfer from the orbit to the product halo via dynamical friction and tidal stripping. The halo spin is correlated with having a recent major merger, though this correlation is weakened by mass loss. These two effects are due to a correlation between the spins of neighboring halos and their orbit, leading to prograde mergers.

Recent evidence appears to link gamma-ray bursts (GRBs) to star-forming regions in galaxies at cosmological distances. If short-lived massive stars are the progenitors of GRBs, the rate of events per unit cosmological volume should be an unbiased tracer (i.e. unaffected by dust obscuration and surface brightness limits) of the cosmic history of star formation. Here we use realistic estimates for the evolution of the stellar birthrate in galaxies to model the number counts, redshift distribution, and time-delay factors of GRBs. We present luminosity function fits to the BATSE log N-log P relation for different redshift distributions of the bursts. Our results imply about 1-2 GRBs every one million Type II supernovae, and a characteristic `isotropic-equivalent' burst luminosity in the range 3-20 x 1e51 ergs/s (for H_0=65 km/s/Mpc). We compute the rate of multiple imaging of background GRBs due to foreground mass condensations in a Lambda-dominated cold dark matter cosmology, assuming that dark halos approximate singular isothermal spheres on galaxy scales and Navarro-Frenk-White profiles on group/cluster scales, and are distributed in mass according to the Press-Schechter model. We show that the expected sensivity increase of Swift relative to BATSE could result in a few strongly lensed individual bursts detected down to a photon flux of 0.1 phot/s/cm^2 in a 3-year survey. Because of the partial sky coverage, however, it is unlikely that the Swift satellite will observe recurrent events (lensed pairs).

We calculate the galaxy bispectrum in both real and redshift space adopting the most common prescriptions for local Eulerian biasing and Lagrangian evolving-bias model. We show that the two biasing schemes make measurably different predictions for these clustering statistics. The Eulerian prescription implies that the galaxy distribution depends only on the present-day local mass distribution, while its Lagrangian counterpart relates the current galaxy distribution to the mass distribution at an earlier epoch when galaxies first formed. Detailed measurement of the galaxy bispectrum (of its reduced amplitude) can help establish whether galaxy positions are determined by the current mass distribution or an earlier mass distribution.

Ongoing searches for supernovae (SNe) at cosmological distances have recently started to provide large numbers of events with measured redshifts and apparent brightnesses. Compared to quasars or galaxies, Type Ia SNe represent a population of sources with well-known intrinsic properties, and could be used to detect gravitational lensing even in the absence of multiple or highly distorted images. We investigate the lensing effect of background SNe due to mass condensations in three popular cold dark matter cosmologies (LambdaCDM, OCDM, SCDM), and compute lensing frequencies, rates of SN explosions, and distributions of arrival time differences and image separations. If dark halos approximate singular isothermal spheres on galaxy scales and NFW profiles on group/cluster scales, and are distributed in mass according to the Press-Schechter theory, then about one every 12 SNe at z~1 will be magnified by Delta m>0.1 mag (SCDM). The detection rate of SN Ia with magnification Delta m>0.3 is estimated to be of order a few events per year per square degree at maximum B-light and I<25, a hundred time smaller than the total rate expected at these magnitude levels. In the field, events magnified by more than 0.75 mag are 7 times less frequent: about one fifth of them gives rise to observable multiple images. Because of the flat K-correction and wide luminosity function, Type II SNe dominate the number counts at I>25 and have the largest fraction of lensed objects. The optimal survey sensitivity for Type Ia's magnified by Delta m>0.75 mag is I=23. Magnification bias increases their incidence by a factor of 50 in samples with I<22, dropping to a factor of 3.5 at 24 mag. At faint magnitudes the enhancement is larger for SN II.

We investigate the clustering of dark matter halos in Lagrangian space in terms of their two-point correlation function, spanning more than 4 orders of magnitudes of halo masses. Analyzing a set of collisionless scale-free 128^3-particle N-body simulations with spectral indices n = -2,-1, we measure the first two Lagrangian bias parameters b_1 and b_2 relating halo and mass correlations. We find that the Mo & White leading-order formula for b_1 describes the clustering of halos with mass M ~ M_* quite accurately, where M_* indicates the characteristic non-linear mass. Smaller halos turn out to be less clustered in Lagrangian space than predicted by Mo & White. Our findings are consistent with the recent results of Jing for the clustering of halo populations in Eulerian space, demonstrating that the discrepancies between the N-body and analytical Mo & White prediction for the bias exist already in Lagrangian space. This shows that a more refined theoretical algorithm for selecting halos in the initial conditions needs to be developed. Finally, we present a very accurate fitting formula for the linear halo bias factor b_1 in Lagrangian space.

We study the spatial distribution of dark matter halos in the Universe in terms of their number density contrast, related to the underlying dark matter fluctuation via a non-local and non-linear bias random field. The description of the matter dynamics is simplified by adopting the `truncated' Zel'dovich approximation to obtain both analytical results and simulated maps. The halo number density field in our maps and its probability distribution reproduce with excellent accuracy those of halos in a high-resolution N-body simulation with the same initial conditions. Our non-linear and non-local bias prescription matches the N-body halo distribution better than any Eulerian linear and local bias.

This paper presents a stochastic approach to the clustering evolution of dark matter haloes in the Universe. Haloes, identified by a Press-Schechter-type algorithm in Lagrangian space, are described in terms of `counting fields', acting as non-linear operators on the underlying Gaussian density fluctuations. By ensemble averaging these counting fields, the standard Press-Schechter mass function as well as analytic expressions for the halo correlation function and corresponding bias factors of linear theory are obtained, thereby extending the recent results by Mo and White. The non-linear evolution of our halo population is then followed by solving the continuity equation, under the sole hypothesis that haloes move by the action of gravity. This leads to an exact and general formula for the bias field of dark matter haloes, defined as the local ratio between their number density contrast and the mass density fluctuation. Besides being a function of position and `observation' redshift, this random field depends upon the mass and formation epoch of the objects and is both non-linear and non-local. The latter features are expected to leave a detectable imprint on the spatial clustering of galaxies, as described, for instance, by statistics like bispectrum and skewness. Our algorithm may have several interesting applications, among which the possibility of generating mock halo catalogues from low-resolution N-body simulations.

We present a stochastic approach to the spatial clustering of dark matter haloes in Lagrangian space. Our formalism is based on a local formulation of the `excursion set' approach by Bond et al., which automatically accounts for the `cloud-in-cloud' problem in the identification of bound systems. Our method allows to calculate correlation functions of haloes in Lagrangian space using either a multi-dimensional Fokker-Planck equation with suitable boundary conditions or an array of Langevin equations with spatially correlated random forces. We compare the results of our method with theoretical predictions for the halo auto-correlation function considered in the literature and find good agreement with the results recently obtained within a treatment of halo clustering in terms of `counting fields' by Catelan et al.. The possible effect of spatial correlations on numerical simulations of halo merger trees is finally discussed.

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.

We analyse the effect of two relevant physical constraints on the mass multiplicity function of dark matter halos in a Press--Schechter type algorithm. Considering the random--walk of linear Gaussian density fluctuations as a function of the smoothing scale, we simultaneously i) account for mass semi--positivity and ii) avoid the cloud--in--cloud problem. It is shown that the former constraint implies a severe cutoff of low--mass objects, balanced by an increase on larger mass scales. The analysis is performed both for scale--free power--spectra and for the standard cold dark matter model. Our approach shows that the well--known ``infrared" divergence of the standard Press--Schechter mass function is caused by unphysical, negative mass events which inevitably occur in a Gaussian distribution of density fluctuations.