We study scaling behaviour of statistics of voids in the context of the halo model of nonlinear large-scale structure. The halo model allows us to understand why the observed galaxy void probability obeys hierarchical scaling, even though the premise from which the scaling is derived is not satisfied. We argue that the commonly observed negative binomial scaling is not fundamental, but merely the result of the specific values of bias and number density for typical galaxies. The model implies quantitative relations between void statistics measured for two populations of galaxies, such as SDSS red and blue galaxies, and their number density and bias.

We present measurements of higher order clustering of galaxies from the latest release of the Canada-France-Hawaii-Telescope Legacy Survey (CFHTLS) Wide. We construct a volume-limited sample of galaxies that contains more than one million galaxies in the redshift range 0.2<z<1 distributed over the four independent fields of the CFHTLS. We use a counts in cells technique to measure the variance <xi_2> and the hierarchical moments S_n = <xi_n>/<xi_2>^(n-1) (3<n<5) as a function of redshift and angular scale.The robustness of our measurements if thoroughly tested, and the field-to-field scatter is in very good agreement with analytical predictions. At small scales, corresponding to the highly non-linear regime, we find a suggestion that the hierarchical moments increase with redshift. At large scales, corresponding to the weakly non-linear regime, measurements are fully consistent with perturbation theory predictions for standard LambdaCDM cosmology with a simple linear bias.

In paper I (Yu et al. [1]), we show through N-body simulation that a local monotonic Gaussian transformation can significantly reduce non-Gaussianity in a noise-free lensing convergence field. This makes the Gaussianization a promising theoretical tool to understand high-order lensing statistics. Here we present a study of its applicability in lensing data analysis, in particular when shape measurement noise is presented in lensing convergence maps. (i) We find that shape measurement noise significantly degrades the Gaussianization performance and the degradation increases for shallower surveys. (ii) The Wiener filter is efficient in reducing the impact of shape measurement noise. The Gaussianization of the Wiener-filtered lensing maps is able to suppress skewness, kurtosis, and the 5th- and 6th-order cumulants by a factor of 10 or more. It also works efficiently to reduce the bispectrum to zero.

Motivated by recent works of Neyrinck et al. 2009 and Scherrer et al. 2010, we proposed a Gaussianization transform to Gaussianize the non-Gaussian lensing convergence field $\kappa$. It performs a local monotonic transformation $\kappa\rightarrow y$ pixel by pixel to make the unsmoothed one-point probability distribution function of the new variable $y$ Gaussian. We tested whether the whole $y$ field is Gaussian against N-body simulations. (1) We found that the proposed Gaussianization suppresses the non-Gaussianity by orders of magnitude, in measures of the skewness, the kurtosis, the 5th- and 6th-order cumulants of the $y$ field smoothed over various angular scales relative to that of the corresponding smoothed $\kappa$ field. The residual non-Gaussianities are often consistent with zero within the statistical errors. (2) The Gaussianization significantly suppresses the bispectrum. Furthermore, the residual scatters around zero, depending on the configuration in the Fourier space. (3) The Gaussianization works with even better performance for the 2D fields of the matter density projected over $\sim 300 \mpch$ distance interval centered at $z\in(0,2)$, which can be reconstructed from the weak lensing tomography. (4) We identified imperfectness and complexities of the proposed Gaussianization. We noticed weak residual non-Gaussianity in the $y$ field. We verified the widely used logarithmic transformation as a good approximation to the Gaussianization transformation. However, we also found noticeable deviations.

We study cell count moments up to fifth order of the distributions of haloes, of halo substructures as a proxy for galaxies, and of mass in the context of the halo model and compare theoretical predictions to the results of numerical simulations. On scales larger than the size of the largest cluster, we present a simple point cluster model in which results depend only on cluster-cluster correlations and on the distribution of the number of objects within a cluster, or cluster occupancy. The point cluster model leads to expressions for moments of galaxy counts in which the volume-averaged moments on large scales approach those of the halo distribution and on smaller scales exhibit hierarchical clustering with amplitudes $S_k$ determined by moments of the occupancy distribution. In this limit, the halo model predictions are purely combinatoric, and have no dependence on halo profile, concentration parameter, or potential asphericity. The full halo model introduces only two additional effects: on large scales, haloes of different mass have different clustering strengths, introducing relative bias parameters; and on the smallest scales, halo structure is resolved and details of the halo profile become important, introducing shape-dependent form factors. Because of differences between discrete and continuous statistics, the hierarchical amplitudes for galaxies and for mass behave differently on small scales even if galaxy number is exactly proportional to mass, a difference that is not necessarily well described in terms of bias.

The amplitude of cosmological density fluctuations, sigma_8, has been studied and estimated by analysing many cosmological observations. The values of the estimates vary considerably between the various probes. However, different estimators probe the value of sigma_8 in different cosmological scales and do not take into account the nonlinear evolution of the parameter at late times. We show that estimates of the amplitude of cosmological density fluctuations derived from cosmic flows are systematically higher than those inferred at early epochs from the CMB because of nonlinear evolution at later times. We discuss the past and future evolution of linear and nonlinear perturbations, derive corrections to the value of sigma_8 and compare amplitudes after accounting for these differences.

We demonstrate that pulsar timing measurements may potentially be able to detect the presence of dark matter substructure within our own galaxy. As dark matter substructure transits near the line-of-sight between a pulsar and an observer, the change in the gravitational field will result in a delay of the light-travel-time of photons. We calculate the effect of this delay due to transiting dark matter substructure and find that the effect on pulsar timing ought to be observable over decadal timescales for a wide range of substructure masses and density profiles. We find that transiting dark matter substructure with masses above 0.01 solar masses ought to be detectable at present by these means. With improved measurements, this method may be able to distinguish between baryonic, thermal non-baryonic, and non-thermal non-baryonic types of dark matter. Additionally, information about structure formation on small scales and the density profiles of galactic dark matter substructure can be extracted via this method.

Although observations point to the neutrality and lack of currents on large scales in the universe, many mechanisms are known that can generate charges or currents during the early universe. We examine the question of survivability of relic charges and currents in a realistic model of the universe. We show that the dynamics of cosmological perturbations drive the universe to become electrically neutral and current-free to a high degree of accuracy on all scales, regardless of initial conditions. We find that charges are efficiently driven away in a time small compared to the Hubble time for temperatures 100 GeV > T > 1 eV, while the same is true for currents at all temperatures T > 1 eV. The forced neutrality relaxes constraints on the generation of an electric charge in the early universe, while the forced erasure of currents disfavors many mechanisms for the early origins of large-scale magnetic fields.

This paper examines the generation of seed magnetic fields due to the growth of cosmological perturbations. In the radiation era, different rates of scattering from photons induce local differences in the ion and electron density and velocity fields. The currents due to the relative motion of these fluids generate magnetic fields on all cosmological scales, peaking at a magnitude of ~ 10^{-24} Gauss at the epoch of recombination. Magnetic fields generated in this manner provide a promising candidate for the seeds of magnetic fields presently observed on galactic and extra-galactic scales.

We evaluate the effect of inhomogeneity energy on the expansion rate of the universe. Our method is to expand to Newtonian order in potential and velocity but to take into account fully nonlinear density inhomogeneities. To linear order in density, kinetic and gravitational potential energy contribute to the total energy of the universe with the same scaling with expansion factor as spatial curvature. In the strongly nonlinear regime, growth saturates, and the net effect of the inhomogeneity energy on the expansion rate remains negligible at all times. In particular, inhomogeneity contributions never mimic the effects of dark energy or induce an accelerated expansion.

Inflationary cosmology predicts a low-amplitude graviton background across a wide range of frequencies. This Letter shows that if one or more extra dimensions exist, the graviton background may have a thermal spectrum instead, dependent on the fundamental scale of the extra dimensions. The energy density is shown to be significant enough that it can affect nucleosynthesis in a substantial way. The possibility of direct detection of a thermal graviton background using the 21-cm hydrogen line is discussed. Alternative explanations for the creation of a thermal graviton background are also examined.

Using mean relative peculiar velocity measurements for pairs of galaxies, we estimate the cosmological density parameter $\Omega_m$ and the amplitude of density fluctuations $\sigma_8$. Our results suggest that our statistic is a robust and reproducible measure of the mean pairwise velocity and thereby the $\Omega_m$ parameter. We get $\Omega_m = 0.30^{+0.17}_{-0.07}$ and $\sigma_8 = 1.13^{+0.22}_{-0.23}$. These estimates do not depend on prior assumptions on the adiabaticity of the initial density fluctuations, the ionization history, or the values of other cosmological parameters.

We derive fully nonlinear expressions for temperature fluctuations from the kinetic Sunyaev-Zeldovich (SZ) effect, the scattering of cosmic microwave background photons off hot electrons in bulk motion. Our result reproduces the Ostriker-Vishniac effect to second order in perturbation theory but contains nonlinear corrections to the electron velocities and densities that were neglected previously. We use the recently developed halo model for nonlinear gravitational clustering to compute the nonlinear kinetic SZ power spectrum, which dominates the primary anisotropy on small angular scales.

We compute the bispectrum for the \IRAS PSCz catalog and find that the galaxy distribution displays the characteristic signature of gravity. Assuming Gaussian initial conditions, we obtain galaxy biasing parameters $1/b_1=1.20^{+0.18}_{-0.19}$ and $b_2/b_1^2=-0.42\pm0.19$, with no sign of scale-dependent bias for $k\leq 0.3$ h/Mpc. These results impose stringent constraints on non-Gaussian initial conditions. For dimensional scaling models with $\chi^2_N$ statistics, we find N>49, which implies a constraint on primordial skewness $B_3<0.35$.

An analytical understanding of the strongly nonlinear regime of gravitational collapse has been difficult to achieve. The only insight has been the stable clustering hypothesis, which assumes that the number of neighbors for objects averaged over small length scales is constant in time. Our recently proposed analytic halo model for N-point correlation functions now provides a tool for calculating gravitational clustering properties in the strongly nonlinear regime. This model also provides a new physical framework for an independent evaluation of the validity of the stable clustering hypothesis. We derive the asymptotic nonlinear behavior of the N-point correlation functions and pairwise peculiar velocities in terms of dark matter halo properties. We show that these statistics exhibit stable clustering only when the halo mass function and halo density profile obey specific relations. The long-cherished stable clustering hypothesis therefore is not necessarily realized in practice.

We present an analytic model for the fully nonlinear power spectrum P and bispectrum Q of the cosmological mass density field. The model is based on physical properties of dark matter halos, with the three main model inputs being analytic halo density profiles, halo mass functions, and halo-halo spatial correlations, each of which has been well studied in the literature. We demonstrate that this new model can reproduce the power spectrum and bispectrum computed from cosmological simulations of both an n=-2 scale-free model and a low-density cold dark matter model. To enhance the dynamic range of these large simulations, we use the synthetic halo replacement technique of Ma & Fry (2000a), where the original halos with numerically softened cores are replaced by synthetic halos of realistic density profiles. At high wavenumbers, our model predicts a slope for the nonlinear power spectrum different from the often-used fitting formulas in the literature based on the stable clustering assumption. Our model also predicts a three-point amplitude Q that is scale dependent, in contrast to the popular hierarchical clustering assumption. This model provides a rapid way to compute the mass power spectrum and bispectrum over all length scales where the input halo properties are valid. It also provides a physical interpretation of the clustering properties of matter in the universe.

We compute the bispectrum for the galaxy distribution in the IRAS QDOT, 2Jy, and 1.2Jy redshift catalogs for wavenumbers 0.05<k<0.2 h/Mpc and compare the results with predictions from gravitational instability in perturbation theory. Taking into account redshift space distortions, nonlinear evolution, the survey selection function, and discreteness and finite volume effects, all three catalogs show evidence for the dependence of the bispectrum on configuration shape predicted by gravitational instability. Assuming Gaussian initial conditions and local biasing parametrized by linear and non-linear bias parameters b_1 and b_2, a likelihood analysis yields 1/b_1 = 1.32^{+0.36}_{-0.58}, 1.15^{+0.39}_{-0.39} and b_2/b_1^2=-0.57^{+0.45}_{-0.30}, -0.50^{+0.31}_{-0.51}, for the for the 2Jy and 1.2Jy samples, respectively. This implies that IRAS galaxies trace dark matter increasingly weakly as the density contrast increases, consistent with their being under-represented in clusters. In a model with chi^2 non-Gaussian initial conditions, the bispectrum displays an amplitude and scale dependence different than that found in the Gaussian case; if IRAS galaxies do not have bias b_1> 1 at large scales, \chi^2 non-Gaussian initial conditions are ruled out at the 95% confidence level. The IRAS data do not distinguish between Lagrangian or Eulerian local bias.

We investigate the nonlinear two- and three-point correlation functions of the cosmological density field in Fourier space and test the popular hierarchical clustering model, that the three-point amplitude Q is independent of scale. In high-resolution N-body simulations of both n=-2 scale-free and cold dark matter models, we find that Q at late times is not constant but increases with wavenumber far into the nonlinear regime. Self-similar scaling also does not hold as rigorously for the three-point function as for the two-point function in the n=-2 simulation, perhaps a manifestation of the finite simulation volume. We suggest that a better understanding of the two- and three-point correlation functions in the nonlinear regime lies in the link to the density profiles of dark matter halos. We demonstrate and quantify how the slopes of the correlation functions are affected by the slope of the halo profiles using simple halo shapes and analytic clustering models.

Second order perturbation theory predicts a specific dependence of the bispectrum, or three-point correlation function in the Fourier transform domain, on the shape of the configuration of its three wave vector arguments, which can be taken as a signature of structure formed by gravitational instability. Comparing this known dependence on configuration shape with the weak shape dependence of the galaxy bispectrum has been suggested as an indication of bias in the galaxy distribution. However, to interpret results obtained from projected catalogs, we must first understand the effects of projection on this shape dependence. We present expressions for the projected power spectrum and bispectrum in both Cartesian and spherical geometries, and we examine the effects of projection on the predicted bispectrum with particular attention to the dependence on configuration shape. Except for an overall numerical factor, for Cartesian projection with characteristic depth $ \Dstar $ there is little effect on the shape dependence of the bispectrum for wavelengths small compared to $ \Dstar $ or projected wavenumbers $ q \Dstar \gg 1 $. For angular projection, a scaling law is found for spherical harmonic index $ \ell \gg 1 $, but there is always a mixing of scales over the range of the selection function. For large $ \ell $ it is sufficient to examine a small portion of the sky.

In efforts to demonstrate the linear Hubble law v = Hr from galaxy observations, the underlying simplicity is often obscured by complexities arising from magnitude-limited data. In this paper we point out a simple but previously unremarked fact: that the shapes and orientations of structures in redshift space contain in themselves independent information about the cosmological redshift-distance relation. The orientations of voids in the CfA slice support the Hubble law, giving a redshift-distance power index p = 0.83 +/- 0.36 (void data from Slezak, de Lapparent, & Bijoui 1993) or p = 0.99 +/- 0.38 (void data from Malik & Subramanian 1997).

The bispectrum, the three-point function of density fluctuations in Fourier space, is the lowest order statistic that carries information about the spatial coherence of large-scale structures. For Gaussian initial conditions, when the density fluctuation amplitude is small, tree-level (leading order) perturbation theory predicts a characteristic dependence of the bispectrum on the shape of the triangle formed by the three wave vectors. This configuration dependence provides a signature of gravitational instability, and departures from it in galaxy catalogs can be interpreted as due to bias, that is, non-gravitational effects. On the other hand, N-body simulations indicate that the reduced three-point function becomes relatively shape-independent in the strongly non-linear regime. In order to understand this non-linear transition and assess the domain of reliability of shape-dependence as a probe of bias, we calculate the one-loop (next-to-leading order) corrections to the bispectrum in perturbation theory. We compare these results with measurements in numerical simulations with scale-free and CDM initial power spectra. We find that the one-loop corrections account very well for the departures from the tree-level results measured in numerical simulations on weakly non-linear scales. In this regime, the reduced bispectrum qualitatively retains its tree-level shape, but the amplitude can change significantly. At smaller scales, the reduced bispectrum in the simulations starts to flatten, an effect which can be partially understood from the one-loop results. In the strong clustering regime, where perturbation theory breaks down entirely, the simulation results confirm that the reduced bispectrum has almost no dependence on triangle shape, in rough agreement with the hierarchical ansatz.

We examine the effect of a threshold bias on the power spectrum and the bispectrum in an ensemble of numerical simulations (Gaussian initial perturbations with power law spectra P(k) \sim k^n, n=+1, 0, -1, -2) and compare our results with theoretical predictions. Our simulations are evolved sufficiently that on the scale where we apply the threshold the rms fluctuation has developed significantly into the nonlinear regime. Thus, predictions based on perturbation theory do not necessarily apply. Nevertheless, we find our results for the power spectrum, biased power simply amplified by a numerical factor, follow predicted trends, far beyond the regime where perturbation theory is expected to be valid. We find that the biased bispectrum continues to follow the so-called hierarchical form, with reduced three-point amplitude Q \approx 1 in the strongly nonlinear regime, independent of initial spectrum. In the quasi-linear perturbative regime the three-point amplitude depends on configuration shape, a behavior that is found to give useful information about the amount of bias without information about the unbiased matter distribution.

To examine how peculiar velocities can affect the 2-, 3-, and 4-point redshift correlation functions, we evaluate volume-average correlations for configurations that emphasize and minimize redshift distortions for four different volume-limited samples from each of the CfA, SSRS, and IRAS redshift catalogs. We find a characteristic distortion for the 2-point correlation, $\xibar_2$: the slope $\gamma$ is flatter and the correlation length is larger in redshift space than in real space; that is, redshift distortions ``move'' correlations from small to large scales. At the largest scales (up to $12 \Mpc$), the extra power in the redshift distribution is compatible with $\Omega^{4/7}/b \approx 1$. We estimate $\Omega^{4/7}/b$ to be $0.53 \pm 0.15$, $1.10 \pm 0.16$ and $0.84 \pm 0.45$ for the CfA, SSRS and IRAS catalogs. Higher order correlations $\xibar_3$ and $\xibar_4$ suffer similar redshift distortions, but in such a way that, within the accuracy of our analysis, the normalized amplitudes $S_3$ and $S_4$ are insensitive to this effect. The hierarchical amplitudes $S_3$ and $S_4$ are constant as a function of scale between 1--$12 \Mpc$ and have similar values in all samples and catalogues, $S_3 \approx 2$ and $S_4 \approx 6$, despite the fact that $\xibar_2$, $\xibar_3$, and $\xibar_4$ differ from one sample to another by large factors (up to a factor of 4 in $\xibar_2$, 8 for $\xibar_3$, and 12 for $\xibar_4$). The agreement between the independent estimations of $S_3$ and $S_4$

In the current paradigm there is a non-trivial bias expected in the process of galaxy formation. Thus, the observed statistical properties of the galaxy distribution do not necessarily extend to the underlying matter distribution. Gravitational evolution of initially Gaussian seed fluctuations predicts that the connected moments of the matter fluctuations exhibit a hierarchical structure, at least in the limit of small dispersion. This same hierarchical structure has been found in the galaxy distribution, but it is not clear to what extent it reflects properties of the matter distribution or properties of a galaxy formation bias. In this paper we consider the consequences of an arbitrary, effectively local biasing transformation of a hierarchical underlying matter distribution. We show that a general form of such a transformation preserves the hierarchical properties and the shape of the dispersion in the limit of small fluctuations, i.e. on large scales, although the values of the