We study the formation of dark matter halos in the concordance LCDM model over a wide range of redshifts, from z=20 to the present. Our primary focus is the halo mass function, a key probe of cosmology. By performing a large suite of nested-box N-body simulations with careful convergence and error controls (60 simulations with box sizes from 4 to 256 Mpc/h, we determine the mass function and its evolution with excellent statistical and systematic errors, reaching a few percent over most of the considered redshift and mass range. Across the studied redshifts, the halo mass is probed over 6 orders of magnitude (10^7 - 10^13.5 M_sun/h). Historically, there has been considerable variation in the high redshift mass function as obtained by different groups. We have made a concerted effort to identify and correct possible systematic errors in computing the mass function at high redshift and to explain the discrepancies between some of the previous results. We discuss convergence criteria for the required force resolution, simulation box size, halo mass range, initial and final redshift, and time stepping. Because of conservative cuts on the mass range probed by individual boxes, our results are relatively insensitive to simulation volume, the remaining sensitivity being consistent with extended Press-Schechter theory. Previously obtained mass function fits near z=0, when scaled by linear theory, are in good agreement with our results at all redshifts, although a mild redshift dependence consistent with that found by Reed and collaborators exists at low redshifts.

We study systematically which features in the cosmic microwave background (CMB) and large-scale structure (LSS) probe various inhomogeneous properties of the dark sectors (including neutrinos, dark matter, and dark energy). We stress, and quantify by simple formulas, that the primary CMB anisotropies are very susceptible to the gravitational potentials during horizon entry, less at recombination. The CMB thus allows us to scan \Phi+\Psi and the underlying dark kinetics for all redshifts z~1-10^5. LSS, on the other hand, responds strongest to \Phi at low redshifts. Dark perturbations are often parameterized by the anisotropic stress and effective sound speed (stiffness). We find that the dark anisotropic stress and stiffness influence the visible species at the correspondingly early and late stages of horizon entry, and affect stronger respectively the CMB and LSS. The CMB yet remains essential to probing the stiff perturbations of light neutrinos and dark energy, detectable only during horizon entry. The clustering of dark species and large propagation speed of their inhomogeneities also map to distinctive features in the CMB and LSS. -Any parameterization of the signatures of dark kinetics that assumes general relativity can effectively accommodate any modified gravity (MG) that retains the equivalence principle for the visible sectors. This implies that formally the nonstandard structure growth or \Phi/\Psi ratio, while indicative, are not definitive MG signatures. The definitive signatures of MG may include the strong dependence of the apparent dark dynamics on visible species, its superluminality, and the nonstandard phenomenology of gravitational waves.

In several classes of particle physics models -- ranging from the classical Majoron models, to the more recent scenarios of late neutrino masses or Mass-Varying Neutrinos -- one or more of the neutrinos are postulated to couple to a new light scalar field. As a result of this coupling, neutrinos in the early universe instead of streaming freely could form a self-coupled fluid, with potentially observable signatures in the Cosmic Microwave Background and the large scale structure of the universe. We re-examine the constraints on this scenario from the presently available cosmological data and investigate the sensitivity expected from the Planck satellite. In the first case, we find that the sensitivity strongly depends on which piece of data is used. The SDSS Main sample data, combined with WMAP and other data, disfavors the scenario of three coupled neutrinos at about the 3.5$\sigma$ confidence level, but also favors a high number of freely streaming neutrinos, with the best fit at 5.2. If the matter power spectrum is instead taken from the SDSS Large Red Galaxy sample, best fit point has 2.5 freely streaming neutrinos, but the scenario with three coupled neutrinos becomes allowed at $2\sigma$. In contrast, Planck alone will exclude even a single self-coupled neutrino at the $4.2\sigma$ confidence level, and will determine the total radiation at CMB epoch to $\Delta N_\nu^{eff} = ^{+0.5}_{-0.3}$ ($1\sigma$ errors). We investigate the robustness of this result with respect to the details of Planck's detector. This sensitivity to neutrino free-streaming implies that Planck will be capable of probing a large region of the Mass-Varying Neutrino parameter space. Planck may also be sensitive to a scale of neutrino mass generation as high as 1 TeV.

First, we establish which measures of large-scale perturbations are least afflicted by gauge artifacts and directly map the apparent evolution of inhomogeneities to local interactions of cosmological species. Considering nonlinear and linear perturbations of phase-space distribution, radiation intensity and arbitrary species' density, we require that: (i) the dynamics of perturbations defined by these measures is determined by observables within the local Hubble volume; (ii) the measures are practically applicable on microscopic scales and in an unperturbed geometry retain their microscopic meaning on all scales. We prove that all measures of linear overdensity that satisfy (i) and (ii) coincide in the superhorizon limit. Their dynamical equations are simpler than the traditional ones, have a nonsingular superhorizon limit and explicit Cauchy form. Then we show that, contrary to the popular view, the perturbations of the cosmic microwave background (CMB) in the radiation era are not resonantly boosted self-gravitationally during horizon entry. (Consequently, the CMB signatures of uncoupled species which may be abundant in the radiation era, e.g. neutrinos or early quintessence, are mild; albeit non-degenerate and robust to cosmic variance.) On the other hand, dark matter perturbations in the matter era gravitationally suppress large-angle CMB anisotropy by an order of magnitude stronger than presently believed. If cold dark matter were the only dominant component then, for adiabatic perturbations, the CMB temperature power spectrum C_l would be suppressed 25-fold.

We describe analytically the cosmological evolution of primordial gravity waves (tensor perturbations) taking into account their coupling to relic neutrinos. We prove that as a consequence of causality the neutrino-induced phase shift of subhorizon tensor oscillations tends on small scales to zero. For the tensor modes that reenter the horizon in the radiation era after neutrino decoupling we calculate the neutrino suppression factor as 1-5rho_nu/9rho+O[(rho_nu/rho)^2]. This result is consistent with the value obtained for three neutrino flavors by Weinberg and is in agreement with numerical Boltzmann evolution. A minimal formula with the correct asymptotic form on small and large scales reproduces to about ten percent the evolution on all scales probed by the CMB. A more accurate solution (in terms of elementary functions) shows that the modes reentering the horizon in the radiation era are slightly enhanced and the phase of their temporal oscillations is shifted by subdominant nonrelativistic matter. The phase shift grows logarithmically on subhorizon scales until radiation-matter equality; the accumulated shift scales for k>>k_eq as (ln k)/k. The modes reentering the horizon after equality are, in turn, affected by the residual radiation density. These modes follow the naive matter-era evolution which is advanced by a redshift and scale independent increment of conformal time. In an appendix, we introduce a general relativistic measure of radiation intensity that in any gauge obeys a simple transport equation and for decoupled particles is conserved on superhorizon scales for arbitrary initial conditions in the full nonlinear theory.

When the perturbations forming the acoustic peaks of the cosmic microwave background (CMB) reentered the horizon and interacted gravitationally with all the matter, neutrinos presumably comprised 41% of the universe energy. CMB experiments have reached a capacity to probe this background of relic neutrinos. I review the neutrino impact on the CMB at the onset of the acoustic oscillations. The discussion addresses the underlying physics, robustness or degeneracy of the neutrino imprints with changes of free cosmological parameters, and non-minimal models for the unseen radiation sector with detectable signatures in CMB anisotropy and polarization.

We present a detailed analytical study of ultra-relativistic neutrinos in cosmological perturbation theory and of the observable signatures of inhomogeneities in the cosmic neutrino background. We note that a modification of perturbation variables that removes all the time derivatives of scalar gravitational potentials from the dynamical equations simplifies their solution notably. The used perturbations of particle number per coordinate, not proper, volume are generally constant on superhorizon scales. In real space an analytical analysis can be extended beyond fluids to neutrinos. The faster cosmological expansion due to the neutrino background changes the acoustic and damping angular scales of the cosmic microwave background (CMB). But we find that equivalent changes can be produced by varying other standard parameters, including the primordial helium abundance. The low-l integrated Sachs-Wolfe effect is also not sensitive to neutrinos. However, the gravity of neutrino perturbations suppresses the CMB acoustic peaks for the multipoles with l>~200 while it enhances the amplitude of matter fluctuations on these scales. In addition, the perturbations of relativistic neutrinos generate a *unique phase shift* of the CMB acoustic oscillations that for adiabatic initial conditions cannot be caused by any other standard physics. The origin of the shift is traced to neutrino free-streaming velocity exceeding the sound speed of the photon-baryon plasma. We find that from a high resolution, low noise instrument such as CMBPOL the effective number of light neutrino species can be determined with an accuracy of sigma(N_nu) = 0.05 to 0.09, depending on the constraints on the helium abundance.

We show that the linear dynamics of cosmological perturbations can be described by coupled wave equations, allowing their efficient numerical and, in certain limits, analytical integration directly in position space. The linear evolution of any perturbation can then be analyzed with the Green's function method. Prior to hydrogen recombination, assuming tight coupling between photons and baryons, neglecting neutrino perturbations, and taking isentropic (adiabatic) initial conditions, the obtained Green's functions for all metric, density, and velocity perturbations vanish beyond the acoustic horizon. A localized primordial cosmological perturbation expands as an acoustic wave of photon-baryon density perturbation with narrow spikes at its acoustic wavefronts. These spikes provide one of the main contributions to the cosmic microwave background radiation anisotropy on all experimentally accessible scales. The gravitational interaction between cold dark matter and baryons causes a dip in the observed temperature of the radiation at the center of the initial perturbation. We first model the radiation by a perfect fluid and then extend our analysis to account for finite photon mean free path. The resulting diffusive corrections smear the sharp features in the photon and baryon density Green's functions over the scale of Silk damping.

We suggest that the cosmic microwave background (CMB) temperature correlation function C(theta) as a function of angle provides a direct connection between experimental data and the fundamental cosmological quantities. The evolution of inhomogeneities in the prerecombination universe is studied using their Green's functions in position space. We find that a primordial adiabatic point perturbation propagates as a sharp-edged spherical acoustic wave. Density singularities at its wavefronts create a feature in the CMB correlation function distinguished by a dip at theta ~ 1.2 deg. Characteristics of the feature are sensitive to the values of cosmological parameters, in particular to the total and the baryon densities.