We present a Fisher information study of the statistical impact of galaxy bias and selection effects on the estimation of key cosmological parameters from galaxy redshift surveys; in particular, the angular diameter distance, Hubble parameter, and linear growth rate at a given redshift, the cold dark matter density, and the tilt and running of the primordial power spectrum. The line-of-sight-dependent selection contributions we include here are known to exist in real galaxy samples. We determine the maximum wavenumber included in the analysis by requiring that the next-order corrections to the galaxy power spectrum or bispectrum, treated here at next-to-leading and leading order, respectively, produce shifts of $\lesssim 0.25\sigma$ on each of the six cosmological parameters. With the galaxy power spectrum alone, selection effects can deteriorate the constraints severely, especially on the linear growth rate. Adding the galaxy bispectrum helps break parameter degeneracies significantly. We find that a joint power spectrum-bispectrum analysis of a Euclid-like survey can still measure the linear growth rate to 10% precision after complete marginalization over selection bias. We also discuss systematic parameter shifts arising from ignoring selection effects and/or other bias parameters, and emphasize that it is necessary to either control selection effects at the percent level or marginalize over them. We obtain similar results for the Roman Space Telescope and HETDEX.

Multiple tracers of the same surveyed volume can enhance the signal-to-noise on a measurement of local primordial non-Gaussianity and the relativistic projections. Increasing the number of tracers comparably increases the number of shot noise terms required to describe the stochasticity of the data. Although the shot noise is white on large scales, it is desirable to investigate the extent to which it can degrade constraints on the parameters of interest. In a multi-tracer analysis of the power spectrum, a marginalization over shot noise does not degrade the constraints on $f_\text{NL}$ by more than $\sim 30$% so long as halos of mass $M\lesssim 10^{12}M_\odot$ are resolved. However, ignoring cross shot noise terms induces large systematics on a measurement of $f_\text{NL}$ at redshift $z<1$ when small mass halos are resolved. These effects are less severe for the relativistic projections, especially for the dipole term. In the case of a low and high mass tracer, the optimal sample division maximizes the signal-to-noise on $f_\text{NL}$ and the projection effects simultaneously, reducing the errors to the level of $\sim 10$ consecutive mass bins of equal number density. We also emphasize that the non-Poissonian noise corrections that arise from small-scale clustering effects cannot be measured with random dilutions of the data. Therefore, they must either be properly modeled or marginalized over.

If the dark matter is made of ultra-light axions, stable solitonic cores form at the centers of virialized halos. In some range for the mass $m$ of the axion particle, these cores are sufficiently compact and can mimic supermassive black holes (SMBH) residing at galactic nuclei. We use the solitonic core--halo mass relation, validated in numerical simulations, to constrain a new range of allowed axion mass from measurements of the SMBH mass in (pseudo)bulge and bulgeless galaxies. These limits are based on observations of galactic nuclei on scales smaller than 10 pc. Our analysis suggests that $m < 10^{-18}$ eV is ruled out by the data. We briefly discuss whether an attractive self-interaction among axions could alleviate this constraint.

We investigate the extent to which resonances between an oscillating background of ultra-light axion and a binary Keplerian system can affect the motion of the latter. These resonances lead to perturbations in the instantaneous time-of-arrivals, and to secular variations in the period of the binary. While the secular changes at exact resonance have recently been explored, the instantaneous effects have been overlooked. In this paper, we examine the latter using N-body simulations including the external oscillatory forcing induced by the axion background. While the secular effects are restricted to a narrow width near the resonance, the instantaneous changes, albeit strongest close to resonances, are apparent for wide range of configurations. We compute the signal-to-noise ratio (SNR) as a function of semi-major axis for a detection of axion oscillations through the R\{o} mer delay. The latter can be extracted from the time-of-arrivals of binary pulsars. The SNR broadly increases with increasing binary eccentricity in agreement with the secular expectation. However, we find that it differs significantly from the scaling a^{5/2} around the lowest orders of resonance. Future observations could probe these effects away from resonances and, therefore, constrain a much broader range of axion masses provided that binary pulsar systems are found near the central region of our Galaxy, and that the time-or-arrival measurement accuracy reaches < 10 ns

Scalar, tensor waves induce oscillatory perturbations in Keplerian systems which can be probed with measurements of pulsar timing residuals. In this paper, we consider the imprint of coherent oscillations produced by ultra-light axion dark matter on the Roemer time delay. We use the angle-action formalism to calculate the time evolution of the observed signal and its dependence on the orbital parameters and the axion phase. We emphasize the similarity of the expected signal-to-noise with the response of an harmonic oscillator to an external oscillatory driving. We validate our theoretical predictions with numerical simulations. Our results furnish a useful benchmark for numerical codes and analysis procedures and, hopefully, will motivate the search for such imprints in real data.

The coalescence of compact binary stars is expected to produce a stochastic background of gravitational waves (GW) observable with future GW detectors. Such backgrounds are usually characterized by their power spectrum as a function of frequency. Here, we present a method to calculate the full 1-point distribution of strain fluctuations. We focus on time series data, but our approach generalizes to the frequency domain. We illustrate how this probability distribution can be evaluated numerically. In addition, we derive accurate analytical asymptotic expressions for the large strain tail, which demonstrate that it is dominated by the nearest source. As an application, we also calculate the distribution of strain fluctuations for the astrophysical GW background produced by binary mergers of compact stars in the Universe, and quantify the extent to which it deviates from a Gaussian distribution. Our approach could be useful for the spectral shape reconstruction of stochastic GW backgrounds.

Detection of gravitational-wave (GW) sources enables the characterisation of binary compact objects and of their in-spiral. However, other dissipative processes can affect the in-spiral. Here we show that the in-spiral of compact objects through a gaseous common-envelope (CE) arising from an evolved stellar companion produces a novel type of GW-sources, whose evolution is dominated by the dissipative gas dynamical friction effects from the CE, rather than the GW-emission itself. The evolution and properties of the GW-signals differ from those of isolated gas-poor mergers significantly. We find characteristic strains of $\sim10^{-23}$-$10^{-21}$ ($10{\rm kpc}/{D}$) for such sources -- observable by next-generation space-based GW-detectors. The evolution of the GW-signal can serve as a probe of the interior regions of the evolved star, and the final stages of CE-evolution, otherwise inaccessible through other observational means. Moreover, such CE-mergers are frequently followed by observable explosive electromagnetic counterparts and/or the formation of exotic stars.

We discuss a novel mechanism to generate gravitational waves in the early universe. A standard way to produce primordial black holes is to enhance at small-scales the overdensity perturbations generated during inflation. The latter, upon horizon re-entry, collapse into black holes. They must be sizeable enough and are therefore associated to rare peaks. There are however less sizeable and much less rare overdensity peaks which do not end up forming primordial black holes and have a non-spherical shape. Upon collapse, they possess a time-dependent non-vanishing mass quadrupole which gives rise to the generation of gravitational waves. By their nature, such gravitational waves are complementary to those sourced at second-order by the very same scalar perturbations responsible for the formation of the primordial black holes. Their amplitude is nevertheless typically about two orders of magnitude smaller and therefore hardly measurable.

We study the spin of primordial black holes produced by the collapse of large inhomogeneities in the early universe. Since such primordial black holes originate from peaks, that is, from maxima of the local overdensity, we resort to peak theory to obtain the probability distribution of the spin at formation. We show that the spin is a first-order effect in perturbation theory: it results from the action of first-order tidal gravitational fields generating first-order torques upon horizon-crossing, and from the asphericity of the collapsing object. Assuming an ellipsoidal shape, the typical value of the dimensionless parameter $a_{\rm s}=S/G_N M^2$, where $S$ is the spin and $M$ is the mass of the primordial black hole, is about $ (\Omega_{\rm dm}/\pi) \sigma_\delta\sqrt{1-\gamma^2}$. Here, $\sigma^2_\delta$ is the variance of the overdensity at horizon crossing, $\Omega_{\rm dm}$ measures the current abundance of the dark matter and the parameter $\gamma$ is a measure of the width of the power spectrum giving rise to primordial black holes. One has $\gamma=1$ for monochromatic spectra. For these narrow spectra, the suppression arises because the velocity shear, which is strongly correlated with the inertia tensor, tends to align with the principal axis frame of the collapsing object. Typical values of $a_{\rm s}$ are at the percent level.

We use analytic computations to predict the power spectrum as well as the bispectrum of Cosmic Infrared Background (CIB) anisotropies. Our approach is based on the halo model and takes into account the mean luminosity-mass relation. The model is used to forecast the possibility to simultaneously constrain cosmological, CIB and halo occupation distribution (HOD) parameters in the presence of foregrounds. For the analysis we use wavelengths in eight frequency channels between 200 and 900$\;\mathrm{GHz}$ with survey specifications given by Planck and LiteBird. We explore the sensitivity to the model parameters up to multipoles of $\ell =1000$ using auto- and cross-correlations between the different frequency bands. With this setting, cosmological, HOD and CIB parameters can be constrained to a few percent. Galactic dust is modeled by a power law and the shot noise contribution as a frequency dependent amplitude which are marginalized over. We find that dust residuals in the CIB maps only marginally influence constraints on standard cosmological parameters. Furthermore, the bispectrum yields tighter constraints (by a factor four in $1\sigma$ errors) on almost all model parameters while the degeneracy directions are very similar to the ones of the power spectrum. The increase in sensitivity is most pronounced for the sum of the neutrino masses. Due to the similarity of degeneracies a combination of both analysis is not needed for most parameters. This, however, might be due to the simplified bias description generally adopted in such halo model approaches.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

The possibility that primordial black holes (PBHs) are the dark matter (or a fraction thereof) has attracted much attention recently. Their spatial clustering is a fundamental property which determines, among others, whether current observational constraints are evaded within a given mass range, whether merging is significant and whether primordial black holes could generate cosmological structure. We treat them as discrete objects and clarify the issue of their spatial clustering, with an emphasis on short-range exclusion and its impact on their large scale power spectrum. Even if a Poissonian self-pair term is always present in the zero-lag correlation, this does not necessarily imply that primordial black holes are initially Poisson distributed. However, while the initial PBH clustering depends on the detailed shape of the small-scale power spectrum, we argue that it is not relevant for a narrow spectral feature and primordial black hole masses still allowed by observations.

This review presents a comprehensive overview of galaxy bias, that is, the statistical relation between the distribution of galaxies and matter. We focus on large scales where cosmic density fields are quasi-linear. On these scales, the clustering of galaxies can be described by a perturbative bias expansion, and the complicated physics of galaxy formation is absorbed by a finite set of coefficients of the expansion, called bias parameters. The review begins with a detailed derivation of this very important result, which forms the basis of the rigorous perturbative description of galaxy clustering, under the assumptions of General Relativity and Gaussian, adiabatic initial conditions. Key components of the bias expansion are all leading local gravitational observables, which include the matter density but also tidal fields and their time derivatives. We hence expand the definition of local bias to encompass all these contributions. This derivation is followed by a presentation of the peak-background split in its general form, which elucidates the physical meaning of the bias parameters, and a detailed description of the connection between bias parameters and galaxy statistics. We then review the excursion-set formalism and peak theory which provide predictions for the values of the bias parameters. In the remainder of the review, we consider the generalizations of galaxy bias required in the presence of various types of cosmological physics that go beyond pressureless matter with adiabatic, Gaussian initial conditions: primordial non-Gaussianity, massive neutrinos, baryon-CDM isocurvature perturbations, dark energy, and modified gravity. Finally, we discuss how the description of galaxy bias in the galaxies' rest frame is related to clustering statistics measured from the observed angular positions and redshifts in actual galaxy catalogs.

We present the complete expression for the next-to-leading (1-loop) order galaxy power spectrum and the leading-order galaxy bispectrum in redshift space in the general bias expansion, or equivalently the effective field theory of biased tracers. We consistently include all line-of-sight dependent selection effects. These are degenerate with many, but not all, of the redshift-space distortion contributions, and have not been consistently derived before. Moreover, we show that, in the framework of effective field theory, a consistent bias expansion in redshift space must include these selection contributions. Physical arguments about the tracer sample considered and its observational selection have to be used to justify neglecting the selection contributions. In summary, the next-to-leading order galaxy power spectrum and leading-order galaxy bispectrum in the general bias expansion are described by 22 parameters, which reduces to 11 parameters if selection effects can be neglected. All contributions to the power spectrum can be written in terms of 28 independent loop integrals.

We investigate how coherent oscillations backreact on the evolution of the condensate wave function of ultra-light axions in the non-relativistic regime appropriate to cosmic structure formation. The coherent oscillations induce higher harmonics beyond the fundamental mode considered so far when a self-interaction is present, and imprint oscillations in the gravitational potential. We emphasize that the effective self-interaction felt by the slowly-varying envelop of the wave function always differs from the bare Lagrangian interaction potential. We also point out that, in the hydrodynamical formulation of the Gross-Pitaevskii equation, oscillations in the gravitational potential result in an attractive force that counteracts the effect of the quantum pressure arising from the strong delocalization of the particles. Since these effects become significant on physical scales less than the (large) Compton wavelength of the particle, they are presumably not very relevant on the mildly nonlinear scales traced by intergalactic neutral hydrogen for axion masses consistent with the bounds from the Lyman-$\alpha$ forest. However, they might affect the formation of virialized cosmological structures and their stability.

Using cosmological N-body simulations, we study the abundance of local maxima (peaks) and minima (dips) identified in the smoothed distribution of halos and dark matter (DM) on scales of $10-100$s Mpcs. The simulations include Gaussian and local-type $f_{\rm NL}$ non-Gaussian initial conditions. The expression derived in the literature for the abundance (irrespective of height) of peaks for Gaussian fields is surprisingly accurate for the evolved halo and DM density fields for all initial conditions considered. Furthermore, the height distribution is very well fitted by a log-normal on quasi-linear scales. The abundance as a function of scale depends on the cosmological parameters ($H_0$ and background matter densities) through the shape of the power spectrum, but it is insensitive to the clustering amplitude. Further, the abundance in the smoothed halo distribution is substantially different in the non-Gaussian from the Gaussian simulations. The interpretation of this effect is straightforward in terms of the scale dependence of halo bias in non-Gaussian models. The abundance of extrema extracted from three-dimensional large galaxy redshift surveys could be a competitive probe of the cosmological parameters and initial non-Gaussianity. It breaks the degeneracy between $f_{\rm NL}$ and the clustering amplitude, making it complementary to counts of galaxy clusters and peaks in weak-lensing maps.

Ultra-light axions have sparked attention because their tiny mass $m\sim 10^{-22}$ eV, which leads to a Kiloparsec-scale de Broglie wavelength comparable to the size of dwarf galaxy, could alleviate the so-called small-scale crisis of massive cold dark matter (CDM) candidates.However, recent analyses of the Lyman-$\alpha$ forest power spectrum set a tight lower bound on their mass of $m\gtrsim 10^{-21}$ eV which makes them much less relevant from an astrophysical point of view. An important caveat to these numerical studies is that they do not take into account attractive self-interactions among ultra-light axions, which can counteract the quantum "pressure" induced by the strong delocalization of the particles. In this work, we show that even a tiny attractive interaction among ultra-light axions can have a significant impact on the stability of cosmic structures at low redshift. After a brief review of known results about solitons in the absence of gravity, we discuss the stability of filamentary and pancake-like solutions when quantum pressure, attractive interactions and gravity are present. The analysis based on one degree of freedom, namely the breathing mode, reveals that pancakes are stable, while filaments are unstable if the mass per unit length is larger than a critical value. However, we show that pancakes are unstable against transverse perturbations. We expect this to be true for halos and filaments as well. Instabilities driven by the breathing mode will not be seen in the low column density Lyman-$\alpha$ forest unless the axion decay constant is extremely small, $f\lesssim 10^{13}$ GeV. Notwithstanding, axion solitonic cores could leave a detectable signature in the Lyman-$\alpha$ forest if the normalization of the unknown axion core - filament mass relation is $\sim 100$ larger than it is for spherical halos.

Shot noise is an important ingredient to any measurement or theoretical modeling of discrete tracers of the large scale structure. Recent work has shown that the shot noise in the halo power spectrum becomes increasingly sub-Poissonian at high mass. Interestingly, while the halo model predicts a shot noise power spectrum in qualitative agreement with the data, it leads to an unphysical white noise in the cross halo-matter and matter power spectrum. In this work, we show that absorbing all the halo model sources of shot noise into the halo fluctuation field leads to meaningful predictions for the shot noise contributions to halo clustering statistics and remove the unphysical white noise from the cross halo-matter statistics. Our prescription straightforwardly maps onto the general bias expansion, so that the renormalized shot noise terms can be expressed as combinations of the halo model shot noises. Furthermore, we demonstrate that non-Poissonian contributions are related to volume integrals over correlation functions and their response to long-wavelength density perturbations. This leads to a new class of consistency relations for discrete tracers, which appear to be satisfied by our reformulation of the halo model. We test our theoretical predictions against measurements of halo shot noise bispectra extracted from a large suite of numerical simulations. Our model reproduces qualitatively the observed sub-Poissonian noise, although it underestimates the magnitude of this effect.

We delineate the conditions under which the consistency relation for the non-Gaussian bias and the universality of the halo mass function hold in the context of microscopic Lagrangian descriptions of halos. The former is valid provided that the collapse barrier depends only on the physical fields (instead of fields normalized by their variance for example) and explicitly includes the effect of {\it all} physical fields such as the tidal shear. The latter holds provided that the response of the halo number density to a long-wavelength density fluctuation is equivalent to the response induced by shifting the spherical collapse threshold. Our results apply to any Lagrangian halo bias prescription. Effective "moving" barriers, which are ubiquitous in the literature, do not generally satisfy the consistency relation. Microscopic barriers including the tidal shear lead to two additional, second-order Lagrangian bias parameters which ensure that the consistency relation is satisfied. We provide analytic expressions for them.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

We measure the large-scale bias of dark matter halos in simulations with non-Gaussian initial conditions of the local type, and compare this bias to the response of the mass function to a change in the primordial amplitude of fluctuations. The two are found to be consistent, as expected from physical arguments, for three halo-finder algorithms which use different Spherical Overdensity (SO) and Friends-of-Friends (FoF) methods. On the other hand, we find that the commonly used prediction for universal mass functions, that the scale-dependent bias is proportional to the first-order Gaussian Lagrangian bias, does not yield a good agreement with the measurements. For all halo finders, high-mass halos show a non-Gaussian bias suppressed by 10-15% relative to the universal mass function prediction. For SO halos, this deviation changes sign at low masses, where the non-Gaussian bias becomes larger than the universal prediction.

We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 \mu$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations.

The baryon acoustic oscillation (BAO) feature provides an important distance scale for the measurement of the expansion history of the Universe. Theoretical models of the BAO in the distribution of biased tracers of the large scale structure usually rely on an initially linear BAO. With aid of N-body simulations, we demonstrate that the BAO in the initial (Lagrangian) halo 2-point function is significantly sharper than in the linear matter distribution, in agreement with peak theory. Using this approach, we delineate the scale-dependence induced by the higher-derivative and velocity bias before assessing how much of the initial BAO enhancement survives until the collapse epoch. Finally, we discuss the extent to which the velocity or gravity bias, which is also imprinted in the displacement field of halos, affects the contrast of the BAO obtained with a reconstruction.

We study the halo-matter cross bispectrum in the presence of primordial non-Gaussianity of the local type. We restrict ourselves to the squeezed limit, for which the calculation are straightforward, and perform the measurements in the initial conditions of N-body simulations, to mitigate the contamination induced by nonlinear gravitational evolution. Interestingly, the halo-matter cross bispectrum is not trivial even in this simple limit as it is strongly sensitive to the scale-dependence of the quadratic and third-order halo bias. Therefore, it can be used to test biasing prescriptions. We consider three different prescription for halo clustering: excursion set peaks (ESP), local bias and a model in which the halo bias parameters are explicitly derived from a peak-background split. In all cases, the model parameters are fully constrained with statistics other than the cross bispectrum. We measure the cross bispectrum involving one halo fluctuation field and two mass overdensity fields for various halo masses and collapse redshifts. We find that the ESP is in reasonably good agreement with the numerical data, while the other alternatives we consider fail in various cases. This suggests that the scale-dependence of halo bias also is a crucial ingredient to the squeezed limit of the halo bispectrum.

The impact of biasing schemes on the clustering of tracers of the large-scale structure is analytically studied in the weakly nonlinear regime. For this purpose, we use the one-loop approximation of the integrated perturbation theory together with the renormalized bias functions of various, physically motivated Lagrangian bias schemes. These include the halo, peaks and excursion set peaks model, for which we derive useful formulae for the evaluation of their renormalized bias functions. The shapes of the power spectra and correlation functions are affected by the different bias models at the level of a few percent on weakly nonlinear scales. These effects are studied quantitatively both in real and redshift space. The amplitude of the scale-dependent bias in the presence of primordial non-Gaussianity also depends on the details of the bias models. If left unaccounted for, these theoretical uncertainties could affect the robustness of the cosmological constraints extracted from galaxy clustering data.

The angular power spectrum of the cosmic infrared background (CIB) is a sensitive probe of the local primordial bispectrum. CIB measurements are integrated over a large volume so that the scale dependent bias from the primordial non-Gaussianity leaves a strong signal in the CIB power spectrum. Although galactic dust dominates over the non-Gaussian CIB signal, it is possible to mitigate the dust contamination with enough frequency channels, especially if high frequencies such as the Planck 857 GHz channel are available. We show that, in this case, measurements of the cosmic microwave background from future space missions should be able to probe the local bispectrum shape down to an amplitude |f_nl| < 1.

Dark matter halos are the building blocks of the universe as they host galaxies and clusters. The knowledge of the clustering properties of halos is therefore essential for the understanding of the galaxy statistical properties. We derive an effective halo Boltzmann equation which can be used to describe the halo clustering statistics. In particular, we show how the halo Boltzmann equation encodes a statistically biased gravitational force which generates a bias in the peculiar velocities of virialized halos with respect to the underlying dark matter, as recently observed in N-body simulations.

The dark matter halos that host galaxies and clusters form out of initial high-density patches, providing a biased tracer of the linear matter density field. In the simplest local bias approximation, the halo field is treated as a perturbative series in the average overdensity of the Lagrangian patch. In more realistic models, however, additional quantities will affect the clustering of halo-patches, and this expansion becomes a function of several stochastic variables. In this paper, we present a general multivariate expansion scheme that can parametrize the clustering of any biased Lagrangian tracer, given only the variables involved and their symmetry (in our case rotational invariance). This approach is based on an expansion in the orthonormal polynomials associated with the relevant variables, so that no renormalization of the coefficients ever occurs. We provide explicit expression for the series coefficients, or Lagrangian bias parameters, in the case of peaks of the linear density field. As an application of our formalism, we present a simple derivation of the original BBKS formula, and compute the non-Gaussian bias in the presence of a primordial trispectrum of the local shape.

The standard formalism for the co-evolution of halos and dark matter predicts that any initial halo velocity bias rapidly decays to zero. We argue that, when the purpose is to compute statistics like power spectra etc., the coupling in the momentum conservation equation for the biased tracers must be modified. Our new formulation predicts the constancy in time of any statistical halo velocity bias present in the initial conditions, in agreement with peak theory. We test this prediction by studying the evolution of a conserved halo population in N-body simulations. We establish that the initial simulated halo density and velocity statistics show distinct features of the peak model and, thus, deviate from the simple local Lagrangian bias. We demonstrate, for the first time, that the time evolution of their velocity is in tension with the rapid decay expected in the standard approach.

The Lagrangian peaks of a 1D cosmological random field representing dark matter are used as a proxy for a catalogue of biased tracers in order to investigate the small-scale exclusion in the two-halo term. The two-point correlation function of peaks of a given height is numerically estimated and analytical approximations that are valid inside the exclusion zone are derived. The resulting power spectrum of these tracers is investigated and shows clear deviations from Poisson noise at low frequencies. On large scales, the convergence of a perturbative bias expansion is discussed. Finally, we go beyond Gaussian statistics for the initial conditions and investigate the subsequent evolution of the two-point clustering of peaks through their Zel'dovich ballistic displacement, to clarify how exclusion effects mix up with scale-dependencies induced by nonlinear gravitational evolution. While the expected large-scale separation limit is recovered, significant deviations are found in the exclusion zone that tends in particular to be reduced at later times. Even though these findings apply to the clustering of one-dimensional tracers, they provide useful insights into halo exclusion and its impact on the two-halo term.

Theoretical modeling of the redshift-space power spectrum of galaxies is crucially important to correctly extract cosmological information from redshift surveys. The task is complicated by the nonlinear biasing and redshift space distortion effects, which change with halo mass, and by the wide distribution of halo masses and their occupations by galaxies. One of the main modeling challenges is the existence of satellite galaxies that have both radial distribution and large virial velocities inside halos, a phenomenon known as the Finger-of-God effect. We present a model for the galaxy power spectrum of in which we decompose a given galaxy sample into central and satellite galaxies and relate different contributions to 1- and 2-halo terms in a halo model. Our primary goal is to ensure that any parameters that we introduce have physically meaningful values, and are not just fitting parameters. For the 2-halo terms we use the previously developed RSD modeling of halos in the context of distribution function and perturbation theory approach. This term needs to be multiplied by the effect of radial distances and velocities of satellites inside the halo. To this one needs to add the 1-halo terms, which are non-perturbative. We show that the real space 1-halo terms can be modeled as almost constant, with the finite extent of the satellites inside the halo inducing a small k^2R^2 term, where R is related to the size of the halo. We adopt a similar model for FoG in redshift space, ensuring that FoG velocity dispersion is related to the halo mass. For FoG k^2 type expansions do not work and FoG resummation must be used instead. We test several damping functions to model the velocity dispersion FoG effect. Applying the formalism to mock galaxies modeled after the "CMASS" sample of the BOSS survey, we find that our predictions for the redshift-space power spectra are accurate up to k~0.4Mpc/h.

Spectral distortions of the CMB have recently experienced an increased interest. One of the inevitable distortion signals of our cosmological concordance model is created by the cosmological recombination process, just a little before photons last scatter at redshift $z\simeq 1100$. These cosmological recombination lines, emitted by the hydrogen and helium plasma, should still be observable as tiny deviation from the CMB blackbody spectrum in the cm--dm spectral bands. In this paper, we present a forecast for the detectability of the recombination signal with future satellite experiments. We argue that serious consideration for future CMB experiments in space should be given to probing spectral distortions and, in particular, the recombination line signals. The cosmological recombination radiation not only allows determination of standard cosmological parameters, but also provides a direct observational confirmation for one of the key ingredients of our cosmological model: the cosmological recombination history. We show that, with present technology, such experiments are futuristic but feasible. The potential rewards won by opening this new window to the very early universe could be considerable.

The intrinsic alignment of galaxy shapes (by means of their angular momentum) and their cross-correlation with the surrounding dark matter tidal field are investigated using the 160 000, z=1.2 synthetic galaxies extracted from the high-resolution cosmological hydrodynamical simulation Horizon-AGN. One- and two-point statistics of the spin of the stellar component are measured as a function of mass and colour. For the low-mass galaxies, this spin is locally aligned with the tidal field `filamentary' direction while, for the high-mass galaxies, it is perpendicular to both filaments and walls. The bluest galaxies of our synthetic catalog are more strongly correlated with the surrounding tidal field than the reddest galaxies, and this correlation extends up to 10 Mpc/h comoving distance. We also report a correlation of the projected ellipticities of blue, intermediate mass galaxies on a similar scale at a level of 10^(-4) which could be a concern for cosmic shear measurements. We do not report any measurable intrinsic alignments of the reddest galaxies of our sample. This work is a first step toward the use of very realistic catalog of synthetic galaxies to evaluate the contamination of weak lensing measurement by the intrinsic galactic alignments.

With the advent of large scale galaxy surveys, constraints on primordial non-Gaussianity (PNG) are expected to reach ${\cal O}(f_\text{NL}) \sim 1$. In order to fully exploit the potential of these future surveys, a deep theoretical understanding of the signatures imprinted by PNG on the large scale structure of the Universe is necessary. In this paper, we explore the effect of a stochastic moving barrier on the amplitude of the non-Gaussian bias induced by local quadratic PNG. We show that, in the peak approach to halo clustering, the amplitude of the non-Gaussian bias will generally differ from the peak-background split prediction unless the barrier is flat and deterministic. For excursion set peaks with a square-root barrier, which reproduce reasonably well the linear bias $b_1$ and mass function $\bar{n}_\text{h}$ of SO haloes, the non-Gaussian bias amplitude is $\sim 40$% larger than the peak-background split expectation $d\ln\bar{n}_\text{h}/d\ln\sigma_8$ for haloes of mass $\sim 10^{13} {\it h}^{-1}M_\odot$ at $z=0$. Furthermore, we argue that the effect of PNG on squeezed configurations of the halo bispectrum differs significantly from that predicted by standard local bias approaches. Our predictions can be easily confirmed, or invalidated, with N-body simulations.

It has been recently shown that any halo velocity bias present in the initial conditions does not decay to unity, in agreement with predictions from peak theory. However, this is at odds with the standard formalism based on the coupled fluids approximation for the coevolution of dark matter and halos. Starting from conservation laws in phase space, we discuss why the fluid momentum conservation equation for the biased tracers needs to be modified in accordance with the change advocated in Baldauf, Desjacques & Seljak (2014). Our findings indicate that a correct description of the halo properties should properly take into account peak constraints when starting from the Vlasov-Boltzmann equation.

We study the clustering of voids using $N$-body simulations and simple theoretical models. The excursion-set formalism describes fairly well the abundance of voids identified with the watershed algorithm, although the void formation threshold required is quite different from the spherical collapse value. The void cross bias $b_{\rm c} $ is measured and its large-scale value is found to be consistent with the peak background split results. A simple fitting formula for $b_{\rm c} $ is found. We model the void auto-power spectrum taking into account the void biasing and exclusion effect. A good fit to the simulation data is obtained for voids with radii $\gtrsim$ 30 Mpc/$h$, especially when the void biasing model is extended to 1-loop order. However, the best-fit bias parameters do not agree well with the peak-background split results. Being able to fit the void auto-power spectrum is particularly important not only because it is the direct observable in galaxy surveys, but also our method enables us to treat the bias parameters as nuisance parameters, which are sensitive to the techniques used to identify voids.

Filaments of the cosmic web drive spin acquisition of disc galaxies. The point process of filament-type saddle represent best this environment and can be used to revisit the Tidal Torque Theory in the context of an anisotropic peak (saddle) background split. The constrained misalignment between the tidal tensor and the Hessian of the density field generated in the vicinity of filament saddle points simply explains the corresponding transverse and longitudinal point-reflection symmetric geometry of spin distribution. It predicts in particular an azimuthal orientation of the spins of more massive galaxies and spin alignment with the filament for less massive galaxies. Its scale dependence also allows us to relate the transition mass corresponding to the alignment of dark matter halos spin relative to the direction of their neighboring filament to this geometry, and to predict accordingly it s scaling with the mass of non linearity, as was measured in simulations.

We develop a count-in-cells approach to the distribution of ultraviolet background fluctuations that includes source clustering. We demonstrate that an exact expression can be obtained if the clustering of ionizing sources follows the hierarchical ansatz. In this case, the intensity distribution depends solely on their 2-point correlation function. We show that the void scaling function of high redshift mock quasars is consistent with the Negative Binomial form, before applying our formalism to the description of HeII-ionizing fluctuations at the end of helium reionization. The model inputs are the observed quasar luminosity function and 2-point correlation at redshift $z\sim 3$. We find that, for an (comoving) attenuation length $\lesssim $ 55 Mpc, quasar clustering contributes less than 30% of the variance of intensity fluctuations so long as the quasar correlation length does not exceed 15 Mpc. We investigate also the dependence of the intensity distribution on the large-scale environment. Differences in the mean HeII-ionizing intensity between low- and high-density regions could be a factor of few if the sources are highly clustered. An accurate description of quasar demographics and their correlation with strong absorption systems is required to make more precise predictions.

Understanding the biasing between the clustering properties of halos and the underlying dark matter distribution is important for extracting cosmological information from ongoing and upcoming galaxy surveys. While on sufficiently larges scales the halo overdensity is a local function of the mass density fluctuations, on smaller scales the gravitational evolution generates non-local terms in the halo density field. We characterize the magnitude of these contributions at third-order in perturbation theory by identifying the coefficients of the non-local invariant operators, and extend our calculation to include non-local (Lagrangian) terms induced by a peak constraint. We apply our results to describe the scale-dependence of halo bias in cosmologies with massive neutrinos. The inclusion of gravity-induced non-local terms and, especially, a Lagrangian $k^2$-contribution is essential to reproduce the numerical data accurately. We use the peak-background split to derive the numerical values of the various bias coefficients from the excursion set peak mass function. For neutrino masses in the range $0\leq \sum_i m_{\nu_i} \leq 0.6$ eV, we are able to fit the data with a precision of a few percents up to $k=0.3\, h {\rm \,Mpc^{-1}}$ without any free parameter.

A large-scale hydrodynamical cosmological simulation, Horizon-AGN, is used to investigate the alignment between the spin of galaxies and the cosmic filaments above redshift 1.2. The analysis of more than 150 000 galaxies per time step in the redshift range 1.2<z<1.8 with morphological diversity shows that the spin of low-mass blue galaxies is preferentially aligned with their neighbouring filaments, while high-mass red galaxies tend to have a perpendicular spin. The reorientation of the spin of massive galaxies is provided by galaxy mergers, which are significant in their mass build-up. We find that the stellar mass transition from alignment to misalignment happens around 3.10^10 M_sun. Galaxies form in the vorticity-rich neighbourhood of filaments, and migrate towards the nodes of the cosmic web as they convert their orbital angular momentum into spin. The signature of this process can be traced to the properties of galaxies, as measured relative to the cosmic web. We argue that a strong source of feedback such as active galactic nuclei is mandatory to quench in situ star formation in massive galaxies and promote various morphologies. It allows mergers to play their key role by reducing post-merger gas inflows and, therefore, keeping spins misaligned with cosmic filaments.

Direct measurements of peculiar velocities of galaxies and clusters of galaxies can in principle provide explicit information on the three dimensional mass distribution, but this information is modulated by the fact that velocity field is sampled at galaxy positions, and is thus probing galaxy momentum. We derive expressions for the cross power spectrum between the density and momentum field and the auto spectrum of the momentum field in redshift space, by extending the distribution function method to these statistics. The resulting momentum cross and auto power spectra in redshift space are expressed as infinite sums over velocity moment correlators in real space, as is the case for the density power spectrum in redshift space. We compute each correlator using Eulerian perturbation theory (PT) and halo biasing model and compare the resulting redshift-space velocity statistics to those measured from N-body simulations for both dark matter and halos. We find that in redshift space linear theory predictions for the density-momentum cross power spectrum as well as for the momentum auto spectrum fail to predict the N-body results at very large scales. On the other hand, our nonlinear PT prediction for these velocity statistics, together with real-space power spectrum for dark matter from simulations, improves the accuracy for both dark matter and halos. We also present the same analysis in configuration space, computing the redshift-space pairwise mean infall velocities and velocity correlation function and compare to nonlinear PT.

We investigate nonlocal Lagrangian bias contributions involving gradients of the linear density field, for which we have predictions from the excursion set peak formalism. We begin by writing down a bias expansion which includes all the bias terms, including the nonlocal ones. Having checked that the model furnishes a reasonable fit to the halo mass function, we develop a 1-point cross-correlation technique to measure bias factors associated with 2-distributed quantities. We validate the method with numerical realizations of peaks of Gaussian random fields before we apply it to N-body simulations. We focus on the lowest (quadratic) order nonlocal contributions. We can reproduce our measurement of \chi_{10} if we allow for an offset between the Lagrangian halo center-of-mass and the peak position. The sign and magnitude of \chi_{10} is consistent with Lagrangian haloes sitting near linear density maxima. The resulting contribution to the halo bias can safely be ignored for M = 10^13 Msun/h, but could become relevant at larger halo masses. For the second nonlocal bias \chi_{01} however, we measure a much larger magnitude than predicted by our model. We speculate that some of this discrepancy might originate from nonlocal Lagrangian contributions induced by nonspherical collapse.

We explore the scale dependence of halo bias using real space cross-correlation measurements in N-body simulations and in Pinocchio, an algorithm based on Lagrangian Perturbation Theory. Recent work has shown how to interpret such real space measurements in terms of k-dependent bias in Fourier space, and how to remove the k-dependence to reconstruct the k-independent peak-background split halo bias parameters. We compare our reconstruction of the linear bias, which requires no free parameters, with previous estimates from N-body simulations which were obtained directly in Fourier space at large scales, and find very good agreement. Our reconstruction of the quadratic bias is similarly parameter-free, although in this case there are no previous Fourier space measurements to compare with. Our analysis of N-body simulations explicitly tests the predictions of the excursion set peaks (ESP) formalism of Paranjape et al. (2013) for the scale dependence of bias; we find that the ESP predictions accurately describe our measurements. In addition, our measurements in Pinocchio serve as a useful, successful consistency check between Pinocchio and N-body simulations that is not accessible to traditional measurements.

We constrain the most likely range of masses for the Milky Way and M31 using an application of the Numerical Action Method (NAM) that optimizes the fit to observed parameters over a large ensemble of NAM-generated solutions. Our 95% confidence level mass ranges, 1.5-4.5 x 10^12 m_sun for MW and 1.5-5.5 x 10^12 m_sun for M31, are consistent with the upper range of estimates from other methods and suggests that a larger proportion of the total mass becomes detectable when the peculiar motions of many nearby satellites are taken into account in the dynamical analysis. We test the method against simulated Local Group catalogs extracted from the Millennium Run to confirm that mass predictions are consistent with actual galaxy halo masses.

Numerical simulations show that redshift space distortions (RSD) introduce strong scale dependence in the power spectra of halos, with ten percent deviations relative to linear theory predictions even on relatively large scales (k<0.1h/Mpc) and even in the absence of satellites (which induce Fingers-of-God, FoG, effects). If unmodeled these effects prevent one from extracting cosmological information from RSD surveys. In this paper we use perturbation theory (PT) and halo biasing model and apply it to the distribution function approach to RSD, in which RSD is decomposed into several correlators of density weighted velocity moments. We model each of these correlators using PT and compare the results to simulations over a wide range of halo masses and redshifts. We find that with an introduction of a physically motivated halo biasing, and using dark matter power spectra from simulations, we can reproduce the simulation results at a percent level on scales up to k~0.15h/Mpc at z=0, without the need to have free FoG parameters in the model.

Corrections induced by primordial non-Gaussianity to the linear halo bias can be computed from a peak-background split or the widespread local bias model. However, numerical simulations clearly support the prediction of the former, in which the non-Gaussian amplitude is proportional to the linear halo bias. To understand better the reasons behind the failure of standard Lagrangian local bias, in which the halo overdensity is a function of the local mass overdensity only, we explore the effect of a primordial bispectrum on the 2-point correlation of discrete density peaks. We show that the effective local bias expansion to peak clustering vastly simplifies the calculation. We generalize this approach to excursion set peaks and demonstrate that the resulting non-Gaussian amplitude, which is a weighted sum of quadratic bias factors, precisely agrees with the peak-background split expectation, which is a logarithmic derivative of the halo mass function with respect to the normalisation amplitude. We point out that statistics of thresholded regions can be computed using the same formalism. Our results suggest that halo clustering statistics can be modelled consistently (in the sense that the Gaussian and non-Gaussian bias factors agree with peak-background split expectations) from a Lagrangian bias relation only if the latter is specified as a set of constraints imposed on the linear density field. This is clearly not the case of standard Lagrangian local bias. Therefore, one is led to consider additional variables beyond the local mass overdensity.

This paper presents CUBEP3M, a publicly-available high performance cosmological N-body code and describes many utilities and extensions that have been added to the standard package. These include a memory-light runtime SO halo finder, a non-Gaussian initial conditions generator, and a system of unique particle identification. CUBEP3M is fast, its accuracy is tuneable to optimize speed or memory, and has been run on more than 27,000 cores, achieving within a factor of two of ideal weak scaling even at this problem size. The code can be run in an extra-lean mode where the peak memory imprint for large runs is as low as 37 bytes per particles, which is almost two times leaner than other widely used N-body codes. However, load imbalances can increase this requirement by a factor of two, such that fast configurations with all the utilities enabled and load imbalances factored in require between 70 and 120 bytes per particles. CUBEP3M is well designed to study large scales cosmological systems, where imbalances are not too large and adaptive time-stepping not essential. It has already been used for a broad number of science applications that require either large samples of non-linear realizations or very large dark matter N-body simulations, including cosmological reionization, halo formation, baryonic acoustic oscillations, weak lensing or non-Gaussian statistics. We discuss the structure, the accuracy, known systematic effects and the scaling performance of the code and its utilities, when applicable.

We present a derivation of two-point correlations of general tracers in the peak-background split (PBS) framework by way of a rigorous definition of the PBS argument. Our expressions only depend on connected matter correlators and "renormalized" bias parameters with clear physical interpretation, and are independent of any coarse-graining scale. This result should be contrasted with the naive expression derived from a local bias expansion of the tracer number density with respect to the matter density perturbation \delta_L coarse-grained on a scale R_L. In the latter case, the predicted tracer correlation function receives contributions of order <\delta_L^n> at each perturbative order n, whereas, in our formalism, these are absorbed in the PBS bias parameters at all orders. Further, this approach naturally predicts both a scale-dependent bias ~ k^2 such as found for peaks of the density field, and the scale-dependent bias induced by primordial non-Gaussianity in the initial conditions. The only assumption made about the tracers is that their abundance at a given position depends solely on the matter distribution within a finite region around that position.

We perform an all-sky analysis of the general relativistic galaxy power spectrum using the well-developed spherical Fourier decomposition. Spherical Fourier analysis expresses the observed galaxy fluctuation in terms of the spherical harmonics and spherical Bessel functions that are angular and radial eigenfunctions of the Helmholtz equation, providing a natural orthogonal basis for all-sky analysis of the large-scale mode measurements. Accounting for all the relativistic effects in galaxy clustering, we compute the spherical power spectrum and its covariance matrix and compare it to the standard three-dimensional power spectrum to establish a connection. The spherical power spectrum recovers the three-dimensional power spectrum at each wavenumber k with its angular dependence mu_k encoded in angular multipole l, and the contributions of the line-of-sight projection to galaxy clustering such as the gravitational lensing effect can be readily accommodated in the spherical Fourier analysis. A complete list of formulas for computing the relativistic spherical galaxy power spectrum is also presented.

The clustering of galaxies in ongoing and upcoming galaxy surveys contains a wealth of cosmological information, but extracting this information is a non-trivial task since galaxies and their host haloes are stochastic tracers of the matter density field. This stochasticity is usually modeled as Poisson shot noise, which is constant as a function of wavenumber with amplitude given by 1/n, where n is the number density of galaxies. Here we use dark matter haloes in N-body simulations to show evidence for deviations from this simple behaviour and develop models that explain the behaviour of the large scale stochasticity. First, haloes are extended, non-overlapping objects, i.e., their correlation function needs to go to -1 on small scales. This leads to a negative correction to the stochasticity relative to the Poisson value at low wavenumber k, decreasing to zero for wavenumbers large compared to the inverse exclusion scale. Second, haloes show a non-linear enhancement of clustering outside the exclusion scale, leading to a positive stochasticity correction. Both of these effects go to zero for high-k, making the stochasticity scale dependent even for k<0.1 h/Mpc. We show that the corrections in the low-k regime are the same in Eulerian and Lagrangian space, but that the transition scale is pushed to smaller scales for haloes observed at present time, relative to the initial conditions. These corrections vary with halo mass and redshift. We also discuss simple applications of these effects to the galaxy samples with non-vanishing satellite fraction, where the stochasticity can again deviate strongly from the fiducial Poisson expectation. Overall these effects affect the clustering of galaxies at a level of a few percent even on very large scales and need to be modelled properly if we want to extract high precision cosmological information from the upcoming galaxy surveys.

The primordial non-Gaussianity parameters fNL and tauNL may be scale-dependent. We investigate the capability of future measurements of the CMB mu-distortion, which is very sensitive to small scales, and of the large-scale halo bias to test the running of local non-Gaussianity. We show that, for an experiment such as PIXIE, a measurement of the mu-temperature correlation can pin down the spectral indices n_fNL and n_tauNL to values of the order of 0.3 if fNL = 20 and tauNL = 5000. A similar value can be achieved with an all-sky survey extending to redshift z ~ 1. In the particular case in which the two spectral indices are equal, as predicted in models where the cosmological perturbations are generated by a single-field other than the inflaton, then the 1-sigma error on the scale-dependence of the non-linearity parameters goes down to 0.2.

We describe how to extend the excursion set peaks framework so that its predictions of dark halo abundances and clustering can be compared directly with simulations. These extensions include: a halo mass definition which uses the TopHat filter in real space; the mean dependence of the critical density for collapse delta_c on halo mass m; and the scatter around this mean value. All three of these are motivated by the physics of triaxial rather than spherical collapse. A comparison of the resulting mass function with N-body results shows that, if one uses delta_c(m) and its scatter as determined from simulations, then all three are necessary ingredients for obtaining ~10% accuracy. E.g., assuming a constant value of delta_c with no scatter, as motivated by the physics of spherical collapse, leads to many more massive halos than seen in simulations. The same model is also in excellent agreement with N-body results for the linear halo bias, especially at the high mass end where the traditional peak-background split argument applied to the mass function fit is known to underpredict the measured bias by ~10%. In the excursion set language, our model is about walks centered on special positions (peaks) in the initial conditions -- we discuss what it implies for the usual calculation in which all walks contribute to the statistics.

Maxima of the linear density field form a point process that can be used to understand the spatial distribution of virialized halos that collapsed from initially overdense regions. However, owing to the peak constraint, clustering statistics of discrete density peaks are difficult to evaluate. For this reason, local bias schemes have received considerably more attention in the literature thus far. In this paper, we show that the 2-point correlation function of maxima of a homogeneous and isotropic Gaussian random field can be thought of, up to second order at least, as arising from a local bias expansion formulated in terms of rotationally invariant variables. This expansion relies on a unique smoothing scale, which is the Lagrangian radius of dark matter halos. The great advantage of this local bias approach is that it circumvents the difficult computation of joint probability distributions. We demonstrate that the bias factors associated with these rotational invariants can be computed using a peak-background split argument, in which the background perturbation shifts the corresponding probability distribution functions. Consequently, the bias factors are orthogonal polynomials averaged over those spatial locations that satisfy the peak constraint. In particular, asphericity in the peak profile contributes to the clustering at quadratic and higher order, with bias factors given by generalized Laguerre polynomials. We speculate that our approach remains valid at all orders, and that it can be extended to describe clustering statistics of any point process of a Gaussian random field. Our results will be very useful to model the clustering of discrete tracers with more realistic collapse prescriptions involving the tidal shear for instance.

Using multiple tracers of large-scale structure allows to evade the limitations imposed by sampling variance for some parameters of interest in cosmology. We demonstrate the optimal way of carrying out a multitracer analysis in a galaxy redshift survey by considering the principal components of the shot noise matrix from two-point clustering statistics. We show how to construct two tracers that maximize the benefits of sampling variance and shot noise cancellation using optimal weights. On the basis of high-resolution N-body simulations of dark matter halos we apply this technique to the analysis of redshift-space distortions and demonstrate how constraints on the growth rate of structure formation can be substantially improved. The primary limitation are nonlinear effects, which cause significant biases in the method already at scales of k<0.1h/Mpc, suggesting the need to develop nonlinear models of redshift-space distortions in order to extract the maximum information from future redshift surveys. Nonetheless we find gains of a factor of a few in constraints on the growth rate achievable when merely the linear regime of a galaxy survey like EUCLID is considered.

We study the ability of future weak lensing (WL) surveys to constrain primordial non-Gaussianity of the local type. We use a large ensemble of simulated WL maps with survey specifications relevant to Euclid and LSST. The simulations assume Cold Dark Matter cosmologies that vary certain parameters around fiducial values: the non-Gaussianity parameter f_NL, the matter density parameter Omega_m, the amplitude of the matter power spectrum sigma_8, the spectral index of the primordial power spectrum n_s, and the dark-energy equation-of-state parameter w_0. We assess the sensitivity of the cosmic shear correlation functions, the third-order aperture mass statistics, and the abundance of shear peaks to these parameters. We find that each of the considered probes provides unmarginalized constraints of Delta f_NL ~ 20 on f_NL. Marginalized constraints from any individual WL probe are much weaker due to strong correlations between parameters. However, the parameter errors can be substantially reduced by combining information from different WL probes. Combining all WL probes yields the following marginal (68% confidence level) uncertainties: Delta f_NL ~ 50, Delta Omega_m ~ 0.002, Delta sigma_8 ~ 0.004, Delta n_s ~ 0.007, and Delta w_0 ~ 0.03. We examine the bias induced by neglecting f_NL on the constraints on the other parameters. We find sigma_8 and w_0 to be the most affected. Moreover, neglecting non-Gaussianity leads to a severe underestimation of the uncertainties in the other cosmological parameters.

Multi-field models of inflation predict an inequality between the amplitude tauNL of the collapsed limit of the four-point correlator of the primordial curvature perturbation and the amplitude fNL of the squeezed limit of its three-point correlator. While a convincing detection of non-Gaussianity through the squeezed limit of the three-point correlator would rule out all single-field models, a robust confirmation or disproval of the inequality between tauNL and fNL would provide crucial information about the validity of multi-field models of inflation. In this paper, we discuss to which extent future measurements of the scale-dependence of galaxy bias can test multi-field inflationary scenarios. The strong degeneracy between the effect of a non-vanishing fNL and tauNL on halo bias can be broken by considering multiple tracer populations of the same surveyed volume. If halos down to 1e13 Msun/h are resolved in a survey of volume 25(Gpc/h)^3, then testing multi-field models of inflation at the 3-\sigma level would require, for instance, a detection of tauNL at the level of tauNL~1e5 given a measurement of a local bispectrum with amplitude fNL~10. However, we find that disproving multi-field models of inflation with measurements of the non-Gaussian bias only will be very challenging, unless |fNL| > 80 and one can achieve a halo mass resolution of 1e10 Msun/h.

It was recently shown that the power spectrum in redshift space can be written as a sum of cross-power spectra between number weighted velocity moments. We investigate the properties of these power spectra for simulated galaxies and dark matter halos and compare them to the dark matter power spectra, generalizing the concept of the bias. Because all of the quantities are number weighted this approach is well defined even for sparse systems such as massive halos, in contrasts to the previous approaches to RSD where velocity correlations have been explored. We find that the number density weighting leads to a strong scale dependence of the bias terms for momentum density auto-correlation and cross-correlation with density. This trend becomes more significant for the more biased halos and leads to an enhancement of RSD power relative to the linear theory. Fingers-of-god effects, which in this formalism come from the correlations of the higher order moments beyond the momentum density, lead to smoothing of the power spectrum and can reduce this enhancement of power, but are relatively small for halos with no small-scale velocity dispersion. In comparison, for a more realistic galaxy sample with satellites the velocity dispersion generated by satellite motions inside the halos leads to a larger power suppression on small scales, but this depends on the satellite fraction. We investigate several statistics such as the two-dimensional power spectrum, its multipole moments, its powers of mu^2, and configuration space statistics. Overall we find that the nonlinear effects in realistic galaxy samples such as luminous red galaxies affect the redshift space clustering on very large scales: for example, the quadrupole moment is affected by 10% for k<0.1h/mpc, which means that these effects need to be understood if we want to extract cosmological information from the redshift space distortions.

We explore the possibility of detecting primordial non-Gaussianity of the local type using weak lensing peak counts. We measure the peak abundance in sets of simulated weak lensing maps corresponding to three models f_NL={0, +100, -100}. Using survey specifications similar to those of Euclid and without assuming any knowledge of the lens and source redshifts, we find the peak functions of the non-Gaussian models with f_NL=+-100 to differ by up to 15% from the Gaussian peak function at the high-mass end. For the assumed survey parameters, the probability of fitting an f_NL=0 peak function to the f_NL=+-100 peak functions is less than 0.1%. Assuming the other cosmological parameters known, f_NL can be measured with an error \Delta f_NL ~ 13. It is therefore possible that future weak lensing surveys like Euclid and LSST may detect primordial non-Gaussianity from the abundance of peak counts, and provide complementary information to that obtained from the cosmic microwave background.

Measurement of redshift-space distortions (RSD) offers an attractive method to directly probe the cosmic growth history of density perturbations. A distribution function approach where RSD can be written as a sum over density weighted velocity moment correlators has recently been developed. We use Nbody simulations to investigate the individual contributions and convergence of this expansion for dark matter. If the series is expanded as a function of powers of mu, cosine of the angle between the Fourier mode and line of sight, there are a finite number of terms contributing at each order. We present these terms and investigate their contribution to the total as a function of wavevector k. For mu^2 the correlation between density and momentum dominates on large scales. Higher order corrections, which act as a Finger-of-God (FoG) term, contribute 1% at k~0.015h/Mpc, 10% at k~0.05h/Mpc at z=0, while for k>0.15h/Mpc they dominate and make the total negative. These higher order terms are dominated by density-energy density correlations which contribute negatively to the power, while the contribution from vorticity part of momentum density auto-correlation is an order of magnitude lower. For mu^4 term the dominant term on large scales is the scalar part of momentum density auto-correlation, while higher order terms dominate for k>0.15h/Mpc. For mu^6 and mu^8 we find it has very little power for k<0.15h/Mpc. We also compare the expansion to the full 2D P^ss(k,mu) as well as to their multipoles. For these statistics an infinite number of terms contribute and we find that the expansion achieves percent level accuracy for kmu<0.15h/Mpc at 6th order, but breaks down on smaller scales because the series is no longer perturbative. We explore resummation of the terms into FoG kernels, which extend the convergence up to a factor of 2 in scale. We find that the FoG kernels are approximately Lorentzian.

The relation between the clustering properties of luminous matter in the form of galaxies and the underlying dark matter distribution is of fundamental importance for the interpretation of ongoing and upcoming galaxy surveys. The so called local bias model, where galaxy density is a function of local matter density, is frequently discussed as a means to infer the matter power spectrum or correlation function from the measured galaxy correlation. However, gravitational evolution generates a term quadratic in the tidal tensor and thus non-local in the density field, even if this term is absent in the initial conditions (Lagrangian space). Because the term is quadratic, it contributes as a loop correction to the power spectrum, so the standard linear bias picture still applies on large scales, however, it contributes at leading order to the bispectrum for which it is significant on all scales. Such a term could also be present in Lagrangian space if halo formation were influenced by the tidal field. We measure the corresponding coupling strengths from the matter-matter-halo bispectrum in numerical simulations and find a non-vanishing coefficient for the tidal tensor term. We find no scale dependence of the bias parameters up to k=0.1 h/Mpc and that the tidal effect is increasing with halo mass. While the Lagrangian bias picture is a better description of our results than the Eulerian bias picture, our results suggest that there might be a tidal tensor bias already in the initial conditions. We also find that the coefficients of the quadratic density term deviate quite strongly from the theoretical predictions based on the spherical collapse model and a universal mass function. Both quadratic density and tidal tensor bias terms must be included in the modeling of galaxy clustering of current and future surveys if one wants to achieve the high precision cosmology promise of these datasets.

One of the main signatures of primordial non-Gaussianity of the local type is a scale-dependent correction to the bias of large-scale structure tracers such as galaxies or clusters, whose amplitude depends on the bias of the tracers itself. The dominant source of noise in the power spectrum of the tracers is caused by sampling variance on large scales (where the non-Gaussian signal is strongest) and shot noise arising from their discrete nature. Recent work has argued that one can avoid sampling variance by comparing multiple tracers of different bias, and suppress shot noise by optimally weighting halos of different mass. Here we combine these ideas and investigate how well the signatures of non-Gaussian fluctuations in the primordial potential can be extracted from the two-point correlations of halos and dark matter. On the basis of large $N$-body simulations with local non-Gaussian initial conditions and their halo catalogs we perform a Fisher matrix analysis of the two-point statistics. Compared to the standard analysis, optimal weighting- and multiple-tracer techniques applied to halos can yield up to one order of magnitude improvements in $\fnl$-constraints, even if the underlying dark matter density field is not known. We compare our numerical results to the halo model and find satisfactory agreement. Forecasting the optimal $\fnl$-constraints that can be achieved with our methods when applied to existing and future survey data, we find that a survey of $50h^{-1}\mathrm{Gpc}^3$ volume resolving all halos down to $10^{11}\hMsun$ at $z=1$ will be able to obtain $\sigma_{\fnl}\sim1$ (68% cl), a factor of $\sim20$ improvement over the current limits. Decreasing the minimum mass of resolved halos, increasing the survey volume or obtaining the dark matter maps can further improve these limits, potentially reaching the level of $\sigma_{\fnl}\sim0.1$. (abridged)

We present measurements of the bispectrum of dark matter halos in numerical simulations with non-Gaussian initial conditions of the local type. We show, in the first place, that the overall effect of primordial non-Gaussianity on the halo bispectrum is larger than on the halo power spectrum when all measurable configurations are taken into account. We then compare our measurements with a tree-level perturbative prediction finding good agreement at large scale when the constant Gaussian bias parameter, both linear and quadratic, and their constant non-Gaussian corrections are fitted for. The best-fit values of the Gaussian bias factors and their non-Gaussian, scale-independent corrections are in qualitative agreement with the peak-background split expectations. In particular, we show that the effect of non-Gaussian initial conditions on squeezed configurations is fairly large (up to 30% for f_NL=100 at redshift z=0.5) and results from contributions of similar amplitude induced by the initial matter bispectrum, scale-dependent bias corrections as well as from nonlinear matter bispectrum corrections. We show, in addition, that effects at second order in f_NL are irrelevant for the range of values allowed by CMB and galaxy power spectrum measurements, at least on the scales probed by our simulations. Finally, we present a Fisher matrix analysis to assess the possibility of constraining primordial non-Gaussianity with future measurements of the galaxy bispectrum. We find that a survey with a volume of about 10 cubic Gpc at mean redshift z ~ 1 could provide an error on f_NL of the order of a few. This shows the relevance of a joint analysis of galaxy power spectrum and bispectrum in future redshift surveys.

Recent results of N-body simulations have shown that current theoretical models are not able to correctly predict the amplitude of the scale-dependent halo bias induced by primordial non-Gaussianity, for models going beyond the simplest, local quadratic case. Motivated by these discrepancies, we carefully examine three theoretical approaches based on (1) the statistics of thresholded regions, (2) a peak-background split method based on separation of scales, and (3) a peak-background split method using the conditional mass function. We first demonstrate that the statistics of thresholded regions, which is shown to be equivalent at leading order to a local bias expansion, cannot explain the mass-dependent deviation between theory and N-body simulations. In the two formulations of the peak-background split on the other hand, we identify an important, but previously overlooked, correction to the non-Gaussian bias that strongly depends on halo mass. This new term is in general significant for any primordial non-Gaussianity going beyond the simplest local fNL model. In a separate paper, we compare these new theoretical predictions with N-body simulations, showing good agreement for all simulated types of non-Gaussianity.

The large-scale clustering of galaxies can serve as a probe of primordial non-Gaussianity in the Universe competitive with the anisotropies of the CMB. Here, we present improved theoretical predictions which include an important, previously overlooked correction to the bias. We demonstrate that the new predictions are able to reproduce the results of N-body simulations, explaining the significant departures seen from previous theoretical results. These refined predictions open the way to accurate constraints on primordial physics with large-scale structure surveys.

We develop the halo model of large-scale structure as an accurate tool for probing primordial non-Gaussianity. In this study we focus on understanding the matter clustering at several redshifts. The primordial non-Gaussianity is modeled as a quadratic correction to the local Gaussian potential, and is characterized by the parameter f_NL. In our formulation of the halo model we pay special attention to the effect of halo exclusion, and show that this can potentially solve the long standing problem of excess power on large scales in this model. The model depends on the mass function, clustering and density profiles of halos. We test these ingredients using a large ensemble of high-resolution Gaussian and non-Gaussian numerical simulations. In particular, we provide a first exploration of how density profiles change in the presence of primordial non-Gaussianities. We find that for f_NL positive/negative high mass halos have an increased/decreased core density, so being more/less concentrated than in the Gaussian case. We also examine the halo bias and show that, if the halo model is correct, then there is a small asymmetry in the scale-dependence of the bias on very large scales, which arises because the Gaussian bias must be renormalized. We show that the matter power spectrum is modified by ~2.5% and ~3.5% on scales k~1.0 h/Mpc at z=0 and z=1, respectively. Our halo model calculation reproduces the absolute amplitude to within 10% and the ratio of non-Gaussian to Gaussian spectra to within 1%. We also measure the matter correlation functions and find similarly good agreement between the model and the data. We anticipate that this modeling will be useful for constraining f_NL from measurements of the shear correlation function in future weak lensing surveys such as Euclid.

Models of galaxy and halo clustering commonly assume that the tracers can be treated as a continuous field locally biased with respect to the underlying mass distribution. In the peak model pioneered by BBKS, one considers instead density maxima of the initial, Gaussian mass density field as an approximation to the formation site of virialized objects. In this paper, the peak model is extended in two ways to improve its predictive accuracy. Firstly, we derive the two-point correlation function of initial density peaks up to second order and demonstrate that a peak-background split approach can be applied to obtain the k-independent and k-dependent peak bias factors at all orders. Secondly, we explore the gravitational evolution of the peak correlation function within the Zel'dovich approximation. We show that the local (Lagrangian) bias approach emerges as a special case of the peak model, in which all bias parameters are scale-independent and there is no statistical velocity bias. We apply our formulae to study how the Lagrangian peak biasing, the diffusion due to large scale flows and the mode-coupling due to nonlocal interactions affect the scale dependence of bias from small separations up to the baryon acoustic oscillation (BAO) scale. For 2-sigma density peaks collapsing at z=0.3, our model predicts a ~ 5% residual scale-dependent bias around the acoustic scale that arises mostly from first-order Lagrangian peak biasing (as opposed to second-order gravity mode-coupling). We also search for a scale dependence of bias in the large scale auto-correlation of massive halos extracted from a very large N-body simulation provided by the MICE collaboration. For halos with mass M>10^{14}Msun/h, our measurements demonstrate a scale-dependent bias across the BAO feature which is very well reproduced by a prediction based on the peak model.

In recent work (Seljak, Hamaus and Desjacques 2009) it was found that weighting central halo galaxies by halo mass can significantly suppress their stochasticity relative to the dark matter, well below the Poisson model expectation. In this paper we extend this study with the goal of finding the optimal mass-dependent halo weighting and use $N$-body simulations to perform a general analysis of halo stochasticity and its dependence on halo mass. We investigate the stochasticity matrix, defined as $C_{ij}\equiv<(\delta_i -b_i\delta_m)(\delta_j-b_j\delta_m)>$, where $\delta_m$ is the dark matter overdensity in Fourier space, $\delta_i$ the halo overdensity of the $i$-th halo mass bin and $b_i$ the halo bias. In contrast to the Poisson model predictions we detect nonvanishing correlations between different mass bins. We also find the diagonal terms to be sub-Poissonian for the highest-mass halos. The diagonalization of this matrix results in one large and one low eigenvalue, with the remaining eigenvalues close to the Poisson prediction $1/\bar{n}$, where $\bar{n}$ is the mean halo number density. The eigenmode with the lowest eigenvalue contains most of the information and the corresponding eigenvector provides an optimal weighting function to minimize the stochasticity between halos and dark matter. We find this optimal weighting function to match linear mass weighting at high masses, while at the low-mass end the weights approach a constant whose value depends on the low-mass cut in the halo mass function. Finally, we employ the halo model to derive the stochasticity matrix and the scale-dependent bias from an analytical perspective. It is remarkably successful in reproducing our numerical results and predicts that the stochasticity between halos and the dark matter can be reduced further when going to halo masses lower than we can resolve in current simulations.

We use the ellipsoidal collapse approximation to investigate the nonlinear redshift space evolution of the density field with primordial non-Gaussianity of the local f_{nl}-type. We utilize the joint distribution of eigenvalues of the initial non-Gaussian shear field and evaluate the evolved redshift space probability distribution function (PDF). It is shown that, similar to the real space analysis, the underdense tail of the nonlinear redshift space PDF differs significantly from that for Gaussian initial conditions. We also derive the lowest order correction of the Kaiser's formulain the presence of a non-zero f_{nl}.

Primordial non-Gaussianity is a potentially powerful discriminant of the physical mechanisms that generated the cosmological fluctuations observed today. Any detection of significant non-Gaussianity would thus have profound implications for our understanding of cosmic structure formation. The large scale mass distribution in the Universe is a sensitive probe of the nature of initial conditions. Recent theoretical progress together with rapid developments in observational techniques will enable us to critically confront predictions of inflationary scenarios and set constraints as competitive as those from the Cosmic Microwave Background. In this paper, we review past and current efforts in the search for primordial non-Gaussianity in the large scale structure of the Universe.

The statistics of peaks of the initial, Gaussian density field can be used to interpret the abundance and clustering of massive dark matter haloes. I discuss some recent theoretical results related to their clustering and its redshift evolution. Predictions from the peak model are qualitatively consistent with measurements of the linear bias of high mass haloes, which also show some evidence for a dependence on the halo mass M at fixed peak height. The peak approach also predicts a distinctive scale-dependence in the bias of haloes across the baryon acoustic feature, a measurement of which would provide strong support for its validity. For 2-sigma density peaks collapsing at z=0.3, this residual scale-dependent bias is at the 5-10 percent level and should thus be within reach of very large simulations of structure formation.

We present measurements of the dark matter bispectrum in N-body simulations with non-Gaussian initial conditions of the local kind for a large variety of triangular configurations and compare them with predictions from Eulerian Perturbation Theory up to one-loop corrections. We find that the effects of primordial non-Gaussianity at large scales, when compared to Perturbation Theory, are well described by the initial component of the matter bispectrum, linearly extrapolated at the redshift of interest. In addition, we find that, for f_NL=100, the nonlinear corrections due to non-Gaussian initial conditions are of the order of ~3, 4% for generic triangles up to ~20% for squeezed configurations, at any redshift. We show that the predictions of Perturbation Theory at tree-level fail to describe the simulation results at redshift z=0 already at scales corresponding to k ~ 0.02 - 0.08 h/Mpc, depending on the triangle, while one-loop corrections can significantly extend their validity to smaller scales. At higher redshift, one-loop Perturbation Theory provides indeed quite accurate predictions, particularly with respect to the relative correction due to primordial non-Gaussianity.

Primordial non-Gaussianity is a potentially powerful discriminant of the physical mechanisms that generated the cosmological fluctuations observed today. Any detection of non-Gaussianity would have profound implications for our understanding of cosmic structure formation. In this paper, we review past and current efforts in the search for primordial non-Gaussianity in the large scale structure of the Universe.

We calculate the redshift space correlation function and the power spectrum of density peaks of a Gaussian random field. In the linear regime k < 0.1 h/Mpc, the redshift space power spectrum is P^s_{pk}(k,u) = exp(-f^2 s_{vel}^2 k^2 u^2) * [b_{pk}(k) + b_{vel}(k) f u^2]^2 * P_m(k), where u is the angle with respect to the line of sight, s_{vel} is the one-dimensional velocity dispersion, f is the growth rate, and b_{pk}(k) and b_{vel}(k) are k-dependent linear spatial and velocity bias factors. For peaks, the value of s_{vel} depends upon the functional form of b_{vel}. The peaks model is remarkable because it has unbiased velocities -- peak motions are driven by dark matter flows -- but, in order to achieve this, b_{vel} is k-dependent. We speculate that this is true in general: k-dependence of the spatial bias will lead to k-dependence of b_{vel} even if the biased tracers flow with the dark matter. Because of the k-dependence of the linear bias parameters, standard manipulations applied to the peak model will lead to k-dependent estimates of the growth factor that could erroneously be interpreted as a signature of modified dark energy or gravity. We use the Fisher formalism to show that the constraint on the growth rate f is degraded by a factor of two if one allows for a k-dependent velocity bias of the peak type. We discuss a simple estimate of nonlinear evolution and illustrate the effect of the peak bias on the redshift space multipoles. For k < 0.1 h/Mpc, the peak bias is deterministic but k-dependent, so the configuration space bias is stochastic and scale dependent, both in real and redshift space. We provide expressions for this stochasticity and its evolution (abridged).

We explore the effect of a cubic correction gnl*phi^3 on the mass function and bias of dark matter haloes extracted from a series of large N-body simulations and compare it to theoretical predictions. Such cubic terms can be motivated in scenarios like the curvaton model, in which a large cubic correction can be produced while simultaneously keeping the quadratic fnl*phi^2 correction small. The deviation from the Gaussian halo mass function is in reasonable agreement with the theoretical predictions. The scale-dependent bias correction Delta b_kappa(k,gnl) measured from the auto- and cross-power spectrum of haloes, is similar to the correction in fnl models, but the amplitude is lower than theoretical expectations. Using the compilation of LSS data in Slosar et al. (2008), we obtain for the first time a limit on gnl of -3.5*10^5 < gnl < +8.2*10^5 (at 95% CL). This limit will improve with the future LSS data by 1-2 orders of magnitude, which should test many of the scenarios of this type.

We generalize Doroshkevich's celebrated formulae for the eigenvalues of the initial shear field associated with Gaussian statistics to the local non-Gaussian f_{nl} model. This is possible because, to at least second order in f_{nl}, distributions at fixed overdensity are unchanged from the case f_{nl}=0. We use this generalization to estimate the effect of f_{nl}\ne 0 on the abundance of virialized halos. Halo abundances are expected to be related to the probability that a certain quantity in the initial fluctuation field exceeds a threshold value, and we study two choices for this variable: it can either be the sum of the eigenvalues of the initial deformation tensor (the initial overdensity), or its smallest eigenvalue. The approach based on a critical overdensity yields results which are in excellent agreement with numerical measurements. We then use these same methods to develop approximations describing the sensitivity of void abundances on f_{nl}. While a positive f_{nl} produces more extremely massive halos, it makes fewer extremely large voids. Its effect thus is qualitatively different from a simple rescaling of the normalisation of the density fluctuation field \sigma_8. Therefore, void abundances furnish complementary information to cluster abundances, and a joint comparison of both might provide interesting constraints on primordial non-Gaussianity.

We investigate the effect of primordial non-Gaussianity of the local f_NL type on the auto- and cross-power spectrum of dark matter haloes using simulations of the LCDM cosmology. We perform a series of large N-body simulations of both positive and negative f_NL, spanning the range between 10 and 100. Theoretical models predict a scale-dependent bias correction \Delta b(k,f_NL) that depends on the linear halo bias b(M). We measure the power spectra for a range of halo mass and redshifts covering the relevant range of existing galaxy and quasar populations. We show that auto and cross-correlation analyses of bias are consistent with each other. We find that for low wavenumbers with k<0.03 h/Mpc the theory and the simulations agree well with each other for biased haloes with b(M)>1.5. We show that a scale-independent bias correction improves the comparison between theory and simulations on smaller scales, where the scale-dependent effect rapidly becomes negligible. The current limits on f_NL from Slosar et al. (2008) come mostly from very large scales k<0.01 h/Mpc and, therefore, remain valid. For the halo samples with b(M)<1.5-2 we find that the scale- dependent bias from non-Gaussianity actually exceeds the theoretical predictions. Our results are consistent with the bias correction scaling linearly with f_NL.

Galaxy surveys are one of the most powerful means to extract the cosmological information and for a given volume the attainable precision is determined by the galaxy shot noise sigma_n^2 relative to the power spectrum P. It is generally assumed that shot noise is white and given by the inverse of the number density n. In this paper we argue one may be able to considerably improve upon this: in the halo picture of cosmological structure all of the dark matter is in halos of varying mass and galaxies are formed inside these halos, but for the dark matter mass and momentum conservation guarantee that nonlinear effects cannot develop a white noise in the dark matter power spectrum on large scales. This suggests that with a suitable weighting a similar effect may be achieved for galaxies, suppressing their shot noise. We explore this idea with N-body simulations by weighting central halo galaxies by halo mass and find that the resulting shot noise can be reduced dramatically relative to expectations, with a 10-30 suppression at the highest number density of n=4*10^-3 (Mpc/h)^3 resolved in our simulations. For specific applications other weighting schemes may achieve even better results and for n=3*10^-4(Mpc/h)^3 we can reduce sigma_n^2/P by up to a factor of 10 relative to uniform weighting. These results open up new opportunities to extract cosmological information in galaxy surveys, such as the recently proposed multi-tracer approach to cancel sampling variance, and may have important consequences for the planning of future redshift surveys. Taking full advantage of these findings may require better understanding of galaxy formation process to develop accurate tracers of the halo mass.

We reexamine the two-point correlation of density maxima in Gaussian initial conditions. Spatial derivatives of the linear density correlation, which were ignored in the calculation of Bardeen, Bond, Kaiser & Szalay (1986), are included in our analysis. These functions exhibit large oscillations around the sound horizon scale for generic CDM power spectra. We derive the exact leading-order expression for the correlation of density peaks and demonstrate the contribution of those spatial derivatives. In particular, we show that these functions can modify significantly the baryon acoustic signature of density maxima relative to that of the linear density field. The effect depends upon the exact value of the peak height, the filter shape and size, and the small-scale behaviour of the transfer function. In the LambdaCDM cosmology, for maxima identified in the density field smoothed at mass scale M\approx 10^{12}-10^{14}M_sun/h and with linear threshold height \nu=1.673/\sigma(M), the contrast of the BAO can be a few tens of percent larger than in the linear matter correlation. Overall, the BAO is amplified (damped) for threshold heights larger (less) than unity. Density maxima thus behave quite differently than linearly biased tracers of the density field, whose acoustic signature is a simple scaled version of the linear baryon acoustic oscillation. We also calculate the mean streaming of peak pairs in the quasi-linear regime. We show that the leading-order 2-point correlation and pairwise velocity of density peaks are consistent with a nonlinear, local biasing relation involving gradients of the density field. Biasing will be an important issue in ascertaining how much of the enhancement of the BAO in the primeval correlation of density maxima propagates into the late-time clustering of galaxies.

Large-scale structures originate from coherent motions induced by inhomogeneities in the primeval gravitational potential. Here, we investigate the two-point statistics of the second derivative of the potential, the tidal shear, under the assumption of Gaussianity. We derive an exact closed form expression for the angular averaged, two-point distribution of the shear components which is valid for an arbitrary Lagrangian separation. This result is used to write down the two-point statistics of the shear eigenvalues in compact form. Next, we examine the large-scale asymptotics of the correlation of the shear eigenvalues, and the alignment of the principal axes. The analytic results are in good agreement with measurements obtained from random realizations of the gravitational potential. Finally, we show that a number of two-point distributions of the shear eigenvalues are well approximated by Gaussian bivariates over a wide range of separation and smoothing scales. We speculate that the Gaussian approximation also holds for multiple point distributions of the shear eigenvalues. It is hoped that these results will be relevant for studies aimed at describing the properties of the (evolved) matter distribution in terms of the statistics of the primordial shear field.

N-body simulations have demonstrated a correlation between the properties of haloes and their environment. In this paper, we assess whether the ellipsoidal collapse model can produce a similar dependence. First, we explore the statistical correlation that originates from Gaussian initial conditions. We derive analytic expressions for a number of joint statistics of the shear tensor and estimate the sensitivity of the local characteristics of the shear to the global geometry of the large scale environment. Next, we concentrate on the dynamical aspect of the environmental dependence using a simplified model that takes into account the interaction between a collapsing halo and its environment. We find that the tidal force exerted by the surrounding mass distribution causes haloes embedded in overdense regions to virialize earlier. An effective density threshold whose shape depends on the large scale density provides a good description of this environmental effect. We show that, using this approach, a correlation between formation redshift, large scale bias and environment density naturally arises. The strength of the effect is comparable, albeit smaller, to that seen in simulations. It is largest for low mass haloes and decreases as one goes to higher mass objects. Furthermore, haloes that formed early are substantially more clustered than those that assembled recently. On the other hand, our analytic model predicts a decrease in median formation redshift with increasing environment density, in disagreement with the trend detected in overdense regions. However, our results appear consistent with the behaviour inferred in relatively underdense regions. We argue that the ellipsoidal collapse model may apply in low density environments where nonlinear effects are negligible (abridged).

We present a measurement of the probability distribution function (PDF) of the transmitted flux in the Lya forest from a sample of 3492 quasars included in the SDSS DR3 data release. Our intention is to investigate the sensitivity of the Lya flux PDF as measured from low resolution and low signal-to-noise data to a number of systematic errors such as uncertainties in the mean flux, continuum and noise estimate. The quasar continuum is described by the superposition of a power law and emission lines. We perform a power law continuum fitting on a spectrum-by-spectrum basis, and obtain an average continuum slope of 0.59 +/- 0.36 in the redshift range 2.5<z<3.5. Taking into account the variation in the continuum indices increases the mean flux by 3 and 7 per cent at z=3 and 2.4, respectively, as compared to the values inferred with a single (mean) continuum slope. We compare our measurements to the PDF obtained with mock lognormal spectra, whose statistical properties have been constrained to match the observed Lya flux PDF and power spectrum of high resolution data. Using our power law continuum fitting and the SDSS pipeline noise estimate yields a poor agreement between the observed and mock PDFs. Allowing for a break in the continuum slope and, more importantly, for residual scatter in the continuum level substantially improves the agreement. A decrease of 10-15 per cent in the mean quasar continuum with a typical rms variance at the 20 per cent level can account for the data, provided that the noise excess correction is no larger than 10 per cent.

The morphology of the total gas and the neutral hydrogen (H I) distributions during the cosmological epoch of reionization can be quantified with Minkowski Functionals (MFs) of isodensity surfaces. We compute the MFs from the output of a high-resolution numerical simulation which includes explicit treatment of transfer of UV ionizing radiation. "Galaxies" identified in the simulation using semi-analytic models of galaxy formation are assumed to be the sole sources of UV photons. The MFs of the total gas distribution are well described by the analytic expressions derived for lognormal random fields. The statistical properties of the diffuse H I depend on the gas distribution and on the way ionized regions propagate in the inter-galactic medium (IGM). The deviations of the MFs of the H I distribution from those of a lognormal random field are, therefore, caused by reionization. We use the MFs to discriminate between the various stages of reionization in the simulation. We suggest a simple model of reionization which reproduces the MFs derived from this simulation. Using random realizations of lognormal density fields, we also assess the ability of MFs to distinguish between different reionization scenarios. Our results are relevant to the analysis of future cosmological twenty-one centimeter maps.

We test the ability of the numerical action method (NAM) to recover the individual orbit histories of mass tracers in an expanding universe in a region of radius 26Mpc/h, given the masses and redshift-space coordinates at the present epoch. The mass tracers are represented by dark matter haloes identified in a high resolution N-body simulation of the standard LCDM cosmology. Since previous tests of NAM at this scale have traced the underlying distribution of dark matter particles rather than extended haloes, our study offers an assessment of the accuracy of NAM in a scenario which more closely approximates the complex dynamics of actual galaxy haloes. We show that NAM can recover present-day halo distances with typical errors of less than 3 per cent, compared to 5 per cent errors assuming Hubble flow distances. The total halo mass and the linear bias were both found to be constained at the 50 per cent level. The accuracy of individual orbit reconstructions was limited by the inability of NAM, in some instances, to correctly model the positions of haloes at early times solely on the basis of the redshifts, angular positions, and masses of the haloes at the present epoch. Improvements in the quality of NAM reconstructions may be possible using the present-day three-dimensional halo velocities and distances to further constrain the dynamics. This velocity data is expected to become available for nearby galaxies in the coming generations of observations by SIM and GAIA.

An excess of sight-lines close to Lyman-break galaxies (LBGs) with little or no absorption in QSO absorption spectra has been reported and has been interpreted as the effect of galactic winds on the Intergalactic Medium. We use here numerical simulations to investigate the flux probability function close to plausible sites of LBGs. We show that the flux distribution near our LBGs in the simulation depends strongly on redshift, and is very sensitive to the averaging procedure. We show that a model without galactic winds and a model with a wind bubble size of 0.5Mpc/h (comoving) are equally consistent with the new determination of the conditional flux distribution by Adelberger et al. (2005). Models with the larger bubble sizes suggested by the previous observations of Adelberger et al. (2003) based on a much smaller sample at higher redshift are not consistent with the new data. We therefore argue that the volume filling factor of galactic winds driven by LBGs may be much smaller than previously thought, and that most of the metals responsible for the metal absorption associated with the low column density Lya forest are unlikely to have been ejected by LBGs.

We presents results of joint modeling of the probability distribution function (PDF) and the one-dimensional power spectrum (PS) of the Lya forest flux decrement. The sensitivity of these statistical measures to the shape and amplitude of the linear matter power spectrum is investigated using two variants of the LCDM cosmology. In the first model, the linear power spectrum has a scale-invariant spectral index, whereas in the second, it has a negative running index (RSI). We generate mock catalogs of QSO spectra, and compare their statistical properties to those of the observations at z=3. We perform a joint fit of the power spectrum and PDF. A scale-invariant model with \sigma_8=0.9 matches well the data if the mean IGM temperature is T\leq 15000K. For higher temperature, it tends to overestimate the flux power spectrum over scales k < 0.01s/km. The discrepancy is less severe when the PS alone is fitted. However, models matching the PS alone do not yield a good fit to the PDF. A joint analysis of the flux PS and PDF tightens the constraints on the model parameters and reduces systematic biases. The RSI model is consistent with the observed PS and PDF only if the temperature is T\geq 20000K. The best fit models reproduce the slope and normalisation of the column density distribution, irrespective of the shape and amplitude of the linear power spectrum. They are also consistent with the observed line-width distribution given the large uncertainties. Our joint analysis suggests that \sigma_8 is likely to be in the range 0.7 - 0.9 for a temperature 10000K<T<20000K and a reasonable reionization history.

We present a model for one-dimensional (1D) matter power spectra in redshift space as estimated from data provided along individual lines of sight. We derive analytic expressions for these power spectra in the linear and nonlinear regimes, focusing on redshift distortions arising from peculiar velocities. In the linear regime, redshift distortions enhance the 1D power spectra only on small scales, and do not affect the power on large scales. This is in contrast to the effect of distortions on three-dimensional (3D) power spectra estimated from data in 3D space, where the enhancement is independent of scale. For CDM cosmologies, the 1D power spectra in redshift and real space are similar for wavenumbers $q<0.1h/Mpc$ where both have a spectral index close to unity, independent of the details of the 3D power spectrum. Nonlinear corrections drive the 1D power spectrum in redshift space into a nearly universal shape over scale $q<10h/Mpc$, and suppress the power on small scales as a result of the strong velocity shear and random motions. The redshift space, 1D power spectrum is mostly sensitive to the amplitude of the initial density perturbations. Our results are useful in particular for power spectra computed from the SDSS quasars sample.

We study the effect of galactic outflows on the Lya forest. The winds are modelled as fully ionised spherical bubbles centered around galactic haloes. The observed flux probability distribution and flux power spectrum limit the volume filling factor of bubbles to be less than 10%. We have compared the mean flux as a function of distance from haloes with the Adelberger etal (ASSP) measurement. A model with bubbles of constant radius of 1.5Mpc/h surrounding the most massive haloes matches the high transmissivity at separations <0.5Mpc/h, but the increase of the transmissivity at small scales is more gradual than observed. At separations <0.5Mpc/h, our estimate of the cosmic variance error is $\Delta F\sim 0.3$, 30% higher than that of ASSP. The difficulty in matching the rise in the transmissivity at small separations is caused by residual absorption of neutral hydrogen lying physically outside the bubbles but having a redshift position similar to the haloes. The flux level is thus sensitive to the amplitude of the coherent velocity shear near haloes which increases with halo mass in the simulation. A model where LBGs are starbursts in small mass haloes matches the observations with smaller bubble size than one where more massive haloes host the LBGs. A starburst model with a bubble radius of 1Mpc/h and a volume filling factor of 2% is consistent with the ASSP measurements at the 1-1.5$\sigma$ level. If this model is correct the sharp rise of the transmissivity at small separation in the ASSP sample is due to cosmic variance and is expected to become more moderate for a larger sample.

Stellar winds emitted by rotating massive stars may show two main components: firstly bipolar lobes with low density and fast wind, produced by the higher T$_{\mathrm{eff}}$ and gravity at the poles (`g$_{\mathrm{eff}}$-effect'); secondly, an equatorial disc with a slow dense wind, produced by the stronger opacities at the equator (`$\kappa$-effect'). To see the possible role of this anisotropic wind on the shape of LBV nebulae, we calculate the distribution of the ejected matter in 2 simplified cases: 1) A brief shell ejection. We find that prolate and peanut-shaped hollow nebulae naturally form due to the g$_{\mathrm{eff}}$-effect in rotating stars. 2) A constant wind for a long time. This produces prolate filled nebulae, with a possible strong disc when a bi-stability limit is crossed in the equatorial region. Thus, many features of the $\eta$ Carinae and LBV nebulae are accounted for by the anisotropic ejection from rotating stars.