We show that the number of observed voids in galaxy redshift surveys is a sensitive function of the equation of state of dark energy. Using the Fisher matrix formalism we find the error ellipses in the $w_0-w_a$ plane when the equation of state of dark energy is assumed to be of the form $w_{CPL}(z)=w_0 +w_a z/(1+z)$. We forecast the number of voids to be observed with the ESA Euclid satellite and the NASA WFIRST mission, taking into account updated details of the surveys to reach accurate estimates of their power. The theoretical model for the forecast of the number of voids is based on matches between abundances in simulations and the analytical prediction. To take into account the uncertainties within the model, we marginalize over its free parameters when calculating the Fisher matrices. The addition of the void abundance constraints to the data from Planck, HST and supernova survey data noticeably tighten the $w_0-w_a$ parameter space. We thus quantify the improvement in the constraints due to the use of voids and demonstrate that the void abundance is a sensitive new probe for the dark energy equation of state.

Thompson scattering of cosmic microwave background (CMB) photons off of free electrons during the reionization epoch induces a correlation between the distribution of galaxies and the polarization pattern of the CMB, the magnitude of which is proportional to the quadrupole moment of radiation at the time of scattering. Since the quadrupole moment generated by gravitational waves (GWs) gives rise to a different polarization pattern than that produced by scalar modes, one can put interesting constraints on the strength of GWs on large scales by cross-correlating the small scale galaxy distribution and CMB polarization. We use this method together with Fisher analysis to predict how well future surveys can measure the tensor-to-scalar ratio $r$. We find that with a future CMB experiment with detector noise Delta_P = 2 mu K-arcmin and a beam width theta_FWHM = 2' and a future galaxy survey with limiting magnitude I<25.6 one can measure the tensor-to-scalar ratio with an error sigma_r \simeq 0.09. To measure r \approx 0.01, however, one needs Delta_P \simeq 0.5 mu K-radian and theta_FWHM \simeq 1'. We also investigate a few systematic effects, none of which turn out to add any biases to our estimators, but they increase the error bars by adding to the cosmic variance. The incomplete sky coverage has the most dramatic effect on our constraints on r for large sky cuts, with a reduction in signal-to-noise smaller than one would expect from the naive estimate (S/N)^2 \propto f_sky. Specifically, we find a degradation factor of f_deg=0.32 \pm 0.01 for a sky cut of |b|>10^\circ (f_sky=0.83) and f_deg=0.056 \pm 0.004 for a sky cut of |b|>20^\circ (f_sky=0.66). Nonetheless, given that our method has different systematics than the more conventional method of observing the large scale B modes directly, it may be used as an important check in the case of a detection.

The advent of precise measurements of the cosmic microwave background (CMB) anisotropies has motivated correspondingly precise calculations of the cosmic recombination history. Cosmic recombination proceeds far out of equilibrium because of a "bottleneck" at the $n=2$ level of hydrogen: atoms can only reach the ground state via slow processes: two-photon decay or Lyman-$\alpha$ resonance escape. However, even a small primordial abundance of molecules could have a large effect on the interline opacity in the recombination epoch and lead to an additional route for hydrogen recombination. Therefore, this paper computes the abundance of the H$_2$ molecule during the cosmic recombination epoch. Hydrogen molecules in the ground electronic levels X$^1\Sigma^+_g$ can either form from the excited H$_2$ electronic levels B$^1\Sigma^+_u$ and C$^1\Pi_u$ or through the charged particles H$_2^+$, HeH$^+$ and H$^-$. We follow the transitions among all of these species, resolving the rotational and vibrational sub-levels. Since the energies of the X$^1\Sigma^+_g$--B$^1\Sigma^+_u$ (Lyman band) and X$^1\Sigma^+_g$-C$^1\Pi_u$ (Werner band) transitions are near the Lyman-$\alpha$ energy, the distortion of the CMB spectrum caused by escaped H Lyman-line photons accelerates both the formation and the destruction of H$_2$ due to this channel relative to the thermal rates. This causes the populations of H$_2$ molecules in X$^1\Sigma^+_g$ energy levels to deviate from their thermal equilibrium abundances. We find that the resulting H$_2$ abundance is $10^{-17}$ at $z=1200$ and $10^{-13}$ at $z=800$, which is too small to have any significant influence on the recombination history.

A signature of the dark energy equation of state may be observed in the shape of voids. We estimate the constraints on cosmological parameters that would be determined from the ellipticity distribution of voids from future spectroscopic surveys already planned for the study of large scale structure. The constraints stem from the sensitivity of the distribution of ellipticity to the cosmological parameters through the variance of fluctuations of the density field smoothed at some length scale. This length scale can be chosen to be of the order of the comoving radii of voids at very early times when the fluctuations are Gaussian distributed. We use Fisher estimates to show that the constraints from void ellipticities are promising. Combining these constraints with other traditional methods results in the improvement of the Dark Energy Task Force Figure of Merit on the dark energy parameters by an order of hundred for future experiments. The estimates of these future constraints depend on a number of systematic issues which require further study using simulations. We outline these issues and study the impact of certain observational and theoretical systematics on the forecasted constraints on dark energy parameters.

The merger and accretion probabilities of dark matter halos have so far only been calculated for an infinitesimal time interval. This means that a Monte-Carlo simulation with very small time steps is necessary to find the merger history of a parent halo. In this paper we use the random walk formalism to find the merger and accretion probabilities of halos for a finite time interval. Specifically, we find the number density of halos at an early redshift that will become part of a halo with a specified final mass at a later redshift, given that they underwent $n$ major mergers, $n=0,1,2,...$ . We reduce the problem into an integral equation which we then solve numerically. To ensure the consistency of our formalism we compare the results with Monte-Carlo simulations and find very good agreement. Though we have done our calculation assuming a flat barrier, the more general case can easily be handled using our method. This derivation of finite time merger and accretion probabilities can be used to make more efficient merger trees or implemented directly into analytical models of structure formation and evolution.