In the upcoming decade cadenced wide-field imaging surveys will increase the number of identified $z<0.3$ Type~Ia supernovae (SNe~Ia) from the hundreds to the hundreds of thousands. The increase in the number density and solid-angle coverage of SNe~Ia, in parallel with improvements in the standardization of their absolute magnitudes, now make them competitive probes of the growth of structure and hence of gravity. The peculiar velocity power spectrum is sensitive to the growth index $\gamma$, which captures the effect of gravity on the linear growth of structure through the relation $f=\Omega_M^\gamma$. We present the first projections for the precision in $\gamma$ for a range of realistic SN peculiar-velocity survey scenarios. In the next decade the peculiar velocities of SNe~Ia in the local $z<0.3$ Universe will provide a measure of $\gamma$ to $\pm 0.01$ precision that can definitively distinguish between General Relativity and leading models of alternative gravity.

We present an improved Minimal Variance (MV) method for using a radial peculiar velocity sample to estimate the average of the three-dimensional velocity field over a spherical volume, which leads to an easily interpretable bulk flow measurement. The only assumption required is that the velocity field is irrotational. The resulting bulk flow estimate is particularly insensitive to smaller scale flows. We also introduce a new constraint into the MV method that ensures that bulk flow estimates are independent of the value of the Hubble constant $H_o$; this is important given the tension between the locally measured $H_o$ and that obtained from the cosmic background radiation observations. We apply our method to the \textit{CosmicFlows-3} catalogue and find that, while the bulk flows for shallower spheres are consistent with the standard cosmological model, there is some tension between the bulk flow in a spherical volume with radius $150$\hmpc\ and its expectations; we find only a $\sim 2\%$ chance of obtaining a bulk flow as large or larger in the standard cosmological model with \textit{Planck} parameters

We present an analysis of the two-point peculiar velocity correlation function using data from the CosmicFlows catalogues. The Millennium and MultiDark Planck 2 N-body simulations are used to estimate cosmic variance and uncertainties due to measurement errors. We compare the velocity correlation function to expectations from linear theory to constrain cosmological parameters. Using the maximum likelihood method, we find values of $\Omega_m= 0.315^{+0.205}_{-0.135}$ and $\sigma_8=0.92^{+0.440}_{-0.295}$, consistent with the Planck and Wilkinson Microwave Anisotropy Probe CMB derived estimates. However, we find that the cosmic variance of the correlation function is large and non-Gaussian distributed, making the peculiar velocity correlation function less than ideal as a probe of large-scale structure.

The bulk flow is a volume average of the peculiar velocities and a useful probe of the mass distribution on large scales. The gravitational instability model views the bulk flow as a potential flow that obeys a Maxwellian Distribution. We use two N-body simulations, the LasDamas Carmen and the Horizon Run, to calculate the bulk flows of various sized volumes in the simulation boxes. Once we have the bulk flow velocities as a function of scale, we investigate the mass and gravitational potential distribution around the volume. We found that matter densities can be asymmetrical and difficult to detect in real surveys, however, the gravitational potential and its gradient may provide better tools to investigate the underlying matter distribution. This study shows that bulk flows are indeed potential flows and thus provides information on the flow sources. We also show that bulk flow magnitudes follow a Maxwellian distribution on scales $>10\ h^{-1}$Mpc.

The Cosmicflows-2 catalogue is a compendium of peculiar velocity measurements. While it has many objects in common with the COMPOSITE catalogue, a previously analysed collection of peculiar velocity data found to give an unexpectedly large bulk flow on large scales, the data in Cosmicflows-2 have been reanalysed to ensure consistency between distances measured using different methods. In particular, a focus on accurate distances led the authors of the Cosmicflows-2 to not correct for homogeneous or inhomogeneous Malmquist bias, both or which are corrected for in the COMPOSITE compilation. We find remarkable agreement between the COMPOSITE and the Cosmicflows-2 if the small EFAR sample of clusters located in two dense superclusters is removed from both surveys, giving results that are inconsistent with the $\Lambda$ cold dark matter standard model with Planck central parameters at the 98\% level. On smaller scales we find overall agreement between data sets and consistency with the standard model.

We introduce a new estimator of the peculiar velocity of a galaxy or group of galaxies from redshift and distance estimates. This estimator results in peculiar velocity estimates which are statistically unbiased and that have errors that are Gaussian distributed, thus meeting the assumptions of analyses that rely on individual peculiar velocities. We apply this estimator to the SFI++ and the Cosmicflows-2 catalogs of galaxy distances and, using the fact that peculiar velocity estimates of distant galaxies are error dominated, examine their error distributions, The adoption of the new estimator significantly improves the accuracy and validity of studies of the large-scale peculiar velocity field and eliminates potential systematic biases, thus helping to bring peculiar velocity analysis into the era of precision cosmology. In addition, our method of examining the distribution of velocity errors should provide a useful check of the statistics of large peculiar velocity catalogs, particularly those that are compiled out of data from multiple sources.

In this paper we introduce PkANN, a freely available software package for interpolating the non-linear matter power spectrum, constructed using Artificial Neural Networks (ANNs). Previously, using Halofit to calculate matter power spectrum, we demonstrated that ANNs can make extremely quick and accurate predictions of the power spectrum. Now, using a suite of 6380 N-body simulations spanning 580 cosmologies, we train ANNs to predict the power spectrum over the cosmological parameter space spanning $3\sigma$ confidence level (CL) around the concordance cosmology. When presented with a set of cosmological parameters ($\Omega_{\rm m} h^2, \Omega_{\rm b} h^2, n_s, w, \sigma_8, \sum m_\nu$ and redshift $z$), the trained ANN interpolates the power spectrum for $z\leq2$ at sub-per cent accuracy for modes up to $k\leq0.9\,h\textrm{Mpc}^{-1}$. PkANN is faster than computationally expensive N-body simulations, yet provides a worst-case error $<1$ per cent fit to the non-linear matter power spectrum deduced through N-body simulations. The overall precision of PkANN is set by the accuracy of our N-body simulations, at 5 per cent level for cosmological models with $\sum m_\nu<0.5$ eV for all redshifts $z\leq2$. For models with $\sum m_\nu>0.5$ eV, predictions are expected to be at 5 (10) per cent level for redshifts $z>1$ ($z\leq1$). The PkANN interpolator may be freely downloaded from http://zuserver2.star.ucl.ac.uk/~fba/PkANN

We calculate the cosmic Mach number M - the ratio of the bulk flow of the velocity field on scale R to the velocity dispersion within regions of scale R. M is effectively a measure of the ratio of large-scale to small-scale power and can be a useful tool to constrain the cosmological parameter space. Using a compilation of existing peculiar velocity surveys, we calculate M and compare it to that estimated from mock catalogues extracted from the LasDamas (a LCDM cosmology) numerical simulations. We find agreement with expectations for the LasDamas cosmology at ~ 1.5 sigma CL. We also show that our Mach estimates for the mocks are not biased by selection function effects. To achieve this, we extract dense and nearly-isotropic distributions using Gaussian selection functions with the same width as the characteristic depth of the real surveys, and show that the Mach numbers estimated from the mocks are very similar to the values based on Gaussian profiles of the corresponding widths. We discuss the importance of the survey window functions in estimating their effective depths. We investigate the nonlinear matter power spectrum interpolator PkANN as an alternative to numerical simulations, in the study of Mach number.

The estimation and analysis of large-scale bulk flow moments of peculiar velocity surveys is complicated by non-spherical survey geometry, the non-uniform sampling of the matter velocity field by the survey objects and the typically large measurement errors of the measured line-of-sight velocities. Previously, we have developed an optimal `minimum variance' (MV) weighting scheme for using peculiar velocity data to estimate bulk flow moments for idealized, dense and isotropic surveys with Gaussian radial distributions, that avoids many of these complications. These moments are designed to be easy to interpret and are comparable between surveys. In this paper, we test the robustness of our MV estimators using numerical simulations. Using MV weights, we estimate the bulk flow moments for various mock catalogues extracted from the LasDamas and the Horizon Run numerical simulations and compare these estimates to the moments calculated directly from the simulation boxes. We show that the MV estimators are unbiased and negligibly affected by non-linear flows.

We investigate the interpolation of power spectra of matter fluctuations using Artificial Neural Network (PkANN). We present a new approach to confront small-scale non-linearities in the power spectrum of matter fluctuations. This ever-present and pernicious uncertainty is often the Achilles' heel in cosmological studies and must be reduced if we are to see the advent of precision cosmology in the late-time Universe. We show that an optimally trained artificial neural network (ANN), when presented with a set of cosmological parameters (Omega_m h^2, Omega_b h^2, n_s, w_0, sigma_8, m_nu and redshift z), can provide a worst-case error <=1 per cent (for z<=2) fit to the non-linear matter power spectrum deduced through N-body simulations, for modes up to k<=0.7 h/Mpc. Our power spectrum interpolator is accurate over the entire parameter space. This is a significant improvement over some of the current matter power spectrum calculators. In this paper, we detail how an accurate interpolation of the matter power spectrum is achievable with only a sparsely sampled grid of cosmological parameters. Unlike large-scale N-body simulations which are computationally expensive and/or infeasible, a well-trained ANN can be an extremely quick and reliable tool in interpreting cosmological observations and parameter estimation. This paper is the first in a series. In this method paper, we generate the non-linear matter power spectra using HaloFit and use them as mock observations to train the ANN. This work sets the foundation for Paper II, where a suite of N-body simulations will be used to compute the non-linear matter power spectra at sub-per cent accuracy, in the quasi-non-linear regime 0.1 h/Mpc <= k <= 0.9 h/Mpc. A trained ANN based on this N-body suite will be released for the scientific community.

The peculiar velocities of galaxies are an inherently valuable cosmological probe, providing an unbiased estimate of the distribution of matter on scales much larger than the depth of the survey. Much research interest has been motivated by the high dipole moment of our local peculiar velocity field, which suggests a large scale excess in the matter power spectrum, and can appear to be in some tension with the LCDM model. We use a composite catalogue of 4,537 peculiar velocity measurements with a characteristic depth of 33 h-1 Mpc to estimate the matter power spectrum. We compare the constraints with this method, directly studying the full peculiar velocity catalogue, to results from Macaulay et al. (2011), studying minimum variance moments of the velocity field, as calculated by Watkins, Feldman & Hudson (2009) and Feldman, Watkins & Hudson (2010). We find good agreement with the LCDM model on scales of k > 0.01 h Mpc-1. We find an excess of power on scales of k < 0.01 h Mpc-1, although with a 1 sigma uncertainty which includes the LCDM model. We find that the uncertainty in the excess at these scales is larger than an alternative result studying only moments of the velocity field, which is due to the minimum variance weights used to calculate the moments. At small scales, we are able to clearly discriminate between linear and nonlinear clustering in simulated peculiar velocity catalogues, and find some evidence (although less clear) for linear clustering in the real peculiar velocity data.

We investigate the impact of massive neutrinos on the distribution of matter in the semi-non-linear regime (0.1<k<0.6 h/Mpc). We present a suite of large-scale N-body simulations quantifying the scale dependent suppression of the total matter power spectrum, resulting from the free-streaming of massive neutrinos out of high-density regions. Our simulations show a power suppression of 3.5-90 per cent at k~0.6 h/Mpc for total neutrino mass, m_nu=0.05-1.9 eV respectively. We also discuss the precision levels that future cosmological datasets would have to achieve in order to distinguish the normal and inverted neutrino mass hierarchies.

In this paper we search for a signature of a large scale bulk flow by looking for fluctuations in the magnitudes of distant LRGs. We take a sample of LRGs from the Sloan Digital Sky Survey with redshifts of z>0.08 over a contiguous area of sky. Neighboring LRG magnitudes are averaged together to find the fluctuation in magnitudes as a function of R.A.. The result is a fluctuation of a few percent in flux across roughly 100 degrees. The source of this fluctuation could be from a large scale bulk flow or a systematic in our treatment of the data set, or the data set itself. A bulk flow model is fitted to the observed fluctuation, and the three bulk flow parameters, its direction and magnitude: alpha_b, delta_b, v_b are constrained. We find that the bulk flow direction is consistent with the direction found by other authors, with alpha_b 180, delta_b -50. The bulk flow magnitude however was found to be anomalously large with v_b>4000km/s. The LRG angular selection function cannot be sufficiently taken into account in our analysis with the available data, and may be the source of either the anomalous magnitude of the bulk flow signal, or possibly the entire fluctuation. However, the fluctuation indicates a bulk flow direction very close to those found using other data sets and analyses. Further investigation with upcoming data is required to confirm this detection.

Peculiar velocities are one of the only probes of very large-scale mass density fluctuations in the nearby Universe. We present new "minimal variance" bulk flow measurements based upon the "First Amendment" compilation of 245 Type Ia supernovae (SNe) peculiar velocities and find a bulk flow of 249 +/- 76 km/s in the direction l= 319 +/- 18 deg, b = 7 +/- 14 deg. The SNe bulk flow is consistent with the expectations of \Lambda CDM. However, it is also marginally consistent with the bulk flow of a larger compilation of non-SNe peculiar velocities (Watkins, Feldman, & Hudson 2009). By comparing the SNe peculiar velocities to predictions of the IRAS Point Source Catalog Redshift survey (PSCz) galaxy density field, we find \Omega_{m}^{0.55} \sigma_{8,lin} = 0.40 +/- 0.07, which is in agreement with \Lambda CDM. However, we also show that the PSCz density field fails to account for 150 +/- 43 km/s of the SNe bulk motion.

A large bulk flow, which is in tension with the Lambda Cold Dark Matter ($\Lambda$CDM) cosmological model, has been observed. In this paper, we provide a physically plausible explanation of this bulk flow, based on the assumption that some fraction of the observed dipole in the cosmic microwave background is due to an intrinsic fluctuation, so that the subtraction of the observed dipole leads to a mismatch between the cosmic microwave background (CMB) defined rest frame and the matter rest frame. We investigate a model that takes into account the relative velocity (hereafter the tilted velocity) between the two frames, and develop a Bayesian statistic to explore the likelihood of this tilted velocity. By studying various independent peculiar velocity catalogs, we find that: (1) the magnitude of the tilted velocity $u$ is around 400 km/s, and its direction is close to what is found from previous bulk flow analyses; for most catalogs analysed, u=0 is excluded at about the $2.5 \sigma$ level;(2) constraints on the magnitude of the tilted velocity can result in constraints on the duration of inflation, due to the fact that inflation can neither be too long (no dipole effect) nor too short (very large dipole effect); (3) Under the assumption of a super-horizon isocurvature fluctuation, the constraints on the tilted velocity require that inflation lasts at least 6 e-folds longer (at the 95% confidence interval) than that required to solve the horizon problem. This opens a new window for testing inflation and models of the early Universe from observations of large scale structure.

The peculiar motions of galaxies can be used to infer the distribution of matter in the Universe. It has recently been shown that measurements of the peculiar velocity field indicates an anomalously high bulk flow of galaxies in our local volume. In this paper we find the implications of the high bulk flow for the power spectrum of density fluctuations. We find that analyzing only the dipole moment of the velocity field yields an average power spectrum amplitude which is indeed much higher than the LCDM value. However, by also including shear and octupole moments of the velocity field, and marginalizing over possible values for the growth rate, an average power spectrum amplitude which is consistent with LCDM is recovered. We attempt to infer the shape of the matter power spectrum from moments of the velocity field, and find a slight excess of power on scales ~ h-1 Gpc.

The low order moments, such as the bulk flow and shear, of the large scale peculiar velocity field are sensitive probes of the matter density fluctuations on very large scales. In practice, however, peculiar velocity surveys are usually sparse and noisy, which can lead to the aliasing of small scale power into what is meant to be a probe of the largest scales. Previously, we developed an optimal ``minimum variance'' (MV) weighting scheme, designed to overcome this problem by minimizing the difference between the measured bulk flow (BF) and that which would be measured by an ideal survey. Here we extend this MV analysis to include the shear and octupole moments, which are designed to have almost no correlations between them so that they are virtually orthogonal. We apply this MV analysis to a compilation of all major peculiar velocity surveys, consisting of 4536 measurements. Our estimate of the BF on scales of ~ 100 Mpc/h has a magnitude of |v|= 416 +/- 78 km/s towards Galactic l = 282 degree +/- 11 degree and b = 6 degree +/- 6 degree. This result is in disagreement with LCDM with WMAP5 cosmological parameters at a high confidence level, but is in good agreement with our previous MV result without an orthogonality constraint, showing that the shear and octupole moments did not contaminate the previous BF measurement. The shear and octupole moments are consistent with WMAP5 power spectrum, although the measurement noise is larger for these moments than for the BF. The relatively low shear moments suggest that the sources responsible for the BF are at large distances.

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.

Peculiar velocity surveys have non-uniform spatial distributions of tracers, so that the bulk flow estimated from them does not correspond to that of a simple volume such as a sphere. Thus bulk flow estimates are generally not strictly comparable between surveys, even those whose effective depths are similar. In addition, the sparseness of typical surveys can lead to aliasing of small scale power into what is meant to be a probe of the largest scales. Here we introduce a new method of calculating bulk flow moments where velocities are weighted to give an optimal estimate of the bulk flow of an idealized survey, with the variance of the difference between the estimate and the actual flow being minimized. These "minimum variance" estimates can be designed to estimate the bulk flow on a particular scale with minimal sensitivity to small scale power, and are comparable between surveys. We compile all major peculiar velocity surveys and apply this new method to them. We find that most surveys we studied are highly consistent with each other. Taken together the data suggest that the bulk flow within a Gaussian window of radius 50 Mpc/h is 407 km/s toward l=287 and b=8. The large-scale bulk motion is consistent with predictions from the local density field. This indicates that there are significant density fluctuations on very large scales. A flow of this amplitude on such a large scale is not expected in the WMAP5-normalized LCDM cosmology, for which the predicted one-dimensional r.m.s. velocity is ~110 km/s. The large amplitude of the observed bulk flow favors the upper values of the WMAP5 error-ellipse, but even the point at the top of the WMAP5 95% confidence ellipse predicts a bulk flow which is too low compared to that observed at >98% confidence level.

We propose to use spatial correlations of the kinetic Sunyaev-Zeldovich (KSZ) flux as an estimator of the peculiar velocity power spectrum. In contrast with conventional techniques, our new method does not require measurements of the thermal SZ signal or the X-ray temperature. Moreover, this method has the special advantage that the expected systematic errors are always sub-dominant to statistical errors on all scales and redshifts of interest. We show that future large sky coverage KSZ surveys may allow a peculiar velocity power spectrum estimates of an accuracy reaching ~10%.

We find the nine bulk--flow and shear moments from the SFI++ survey, as well as for subsamples of group and field galaxies. We constrain the velocity power spectrum shape parameter $\Gamma$ in linear theory using these moments. A likelihood function for $\Gamma$ was found after marginalizing over the power spectrum amplitude $\sigma_8\Omega_m^{0.6}$ using constraints obtained from comparisons between redshift surveys and peculiar velocity data. We have estimated the velocity noise $\sigma_*$ from the data since without it our results may be biased. We also performed a statistical analysis of the difference between the field and group catalogues and found that the results from each reflect the same underlying large scale flows. We found that we can constrain the power spectrum shape parameter to be $\Gamma=0.15^{+0.18}_{-0.08}$ for the groups catalogue and $\Gamma=0.09^{+0.04}_{-0.04}$ for the field galaxy catalogue in fair agreement with the value from WMAP.

To study galactic motions on the largest available scales, we require bulk flow moments whose window functions have as narrow a peak as possible and having as small an amplitude as possible outside the peak. Typically the moments found using the maximum likelihood estimate weights do not meet these criteria. We present a new method for calculating weights for moments that essentially allow us to "design" the moment's window function, subject, of course, to the distribution and uncertainties of the available data.

Despite a history that dates back at least a quarter of a century studies of voids in the large--scale structure of the Universe are bedevilled by a major problem: there exist a large number of quite different void--finding algorithms, a fact that has so far got in the way of groups comparing their results without worrying about whether such a comparison in fact makes sense. Because of the recent increased interest in voids, both in very large galaxy surveys and in detailed simulations of cosmic structure formation, this situation is very unfortunate. We here present the first systematic comparison study of thirteen different void finders constructed using particles, haloes, and semi--analytical model galaxies extracted from a subvolume of the Millennium simulation. The study includes many groups that have studied voids over the past decade. We show their results and discuss their differences and agreements. As it turns out, the basic results of the various methods agree very well with each other in that they all locate a major void near the centre of our volume. Voids have very underdense centres, reaching below 10 percent of the mean cosmic density. In addition, those void finders that allow for void galaxies show that those galaxies follow similar trends. For example, the overdensity of void galaxies brighter than $m_B = -20 $ is found to be smaller than about -0.8 by all our void finding algorithms.

We present an analysis comparing the bulk--flow measurements for six recent peculiar velocity surveys, namely, ENEAR, SFI, RFGC, SBF and the Mark III singles and group catalogs. We study whether the direction of the bulk--flow estimates are consistent with each other and construct the full three dimensional bulk--flow vectors for each survey. We show that although the surveys differ in their geometry, galaxy morphologies, distance measures and measurement errors, their bulk flow vectors are expected to be highly correlated and in fact show impressive agreement in all cases. We found a combined weighted mean bulk motion of 330 km s$^{-1}$ $\pm 101$ km s$^{-1}$ toward $l= 234^{\degree}\pm 11^{\degree}$ and $b=12^{\degree}\pm 9^{\degree}$ in a sphere with an effective depth of $\sim4000$ km s$^{-1}$.

We constrain the velocity power spectrum shape parameter $\Gamma$ in linear theory using the nine bulk-flow and shear moments estimated from four recent peculiar velocity surveys. For each survey, a likelihood function for $\Gamma$ was found after marginalizing over the power spectrum amplitude $\sigma_8\Omega_m^{0.6}$ using constraints obtained from comparisons between redshift surveys and peculiar velocity data. In order to maximize the accuracy of our analyses, the velocity noise $\sigma_*$ was estimated directly for each survey. A statistical analysis of the differences between the values of the moments estimated from different surveys showed consistency with theoretical predictions, suggesting that all the surveys investigated reflect the same large scale flows. The peculiar velocity surveys were combined into a composite survey yielding the constraint $\Gamma=0.13^{+0.09}_{-0.05}$. This value is lower than, but consistent with, values obtained using redshift surveys and CMB data.

We study the global distribution and morphology of dark matter voids in a LCDM universe using density fields generated by N-body simulations. Voids are defined as isolated regions of the low-density excursion set specified via density thresholds, the density thresholds being quantified by the corresponding filling factors, i.e., the fraction of the total volume in the excursion set. Our work encompasses a systematic investigation of the void volume function, the volume fraction in voids, and the fitting of voids to corresponding ellipsoids and spheres. We emphasize the relevance of the percolation threshold to the void volume statistics of the density field both in the high redshift, Gaussian random field regime, as well as in the present epoch. By using measures such as the Inverse Porosity, we characterize the quality of ellipsoidal fits to voids, finding that such fits are a poor representation of the larger voids that dominate the volume of the void excursion set.

We show that by directly correlating the cluster kinetic Sunyaev Zeldovich (KSZ) flux, the cluster peculiar velocity power spectrum can be measured to $\sim 10%$ accuracy by future large sky coverage KSZ surveys. This method is almost free of systemics entangled in the usual velocity inversion method. The direct correlation brings extra information of density and velocity clustering. We utilize these information to construct two indicators of the Hubble constant and comoving angular distance and propose a novel method to constrain cosmology.

We describe a new method of overcoming problems inherent in peculiar velocity surveys by using data compression as a filter with which to separate large-scale, linear flows from small-scale noise that biases the results systematically. We demonstrate the effectiveness of our method using realistic catalogs of galaxy velocities drawn from N--body simulations. Our tests show that a likelihood analysis of simulated catalogs that uses all of the information contained in the peculiar velocities results in a bias in the estimation of the power spectrum shape parameter $\Gamma$ and amplitude $\beta$, and that our method of analysis effectively removes this bias. We expect that this new method will cause peculiar velocity surveys to re--emerge as a useful tool to determine cosmological parameters.

Analyses of peculiar velocity surveys face several challenges, including low signal--to--noise in individual velocity measurements and the presence of small--scale, nonlinear flows. I will present three new analyses that attempt to address these inherent problems. The first is geared towards the better understanding of the estimated errors in the surveys, specifically sampling errors, and the resolution of the seeming disagreements between the surveys. Another develops a new statistic that does not suffer from the usual problems and gives robust results that are galaxy--morphology and distance--estimator independent. The third introduces a formalism that allows for the accounting of most of the non--linear signal whereby the signal to noise is increased and small--scale aliasing is removed.

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 have used N-body simulations to develop two independent methods to quantify redshift distortions known as the Bull's Eye effect (large scale infall plus small scale virial motion). This effect depends upon the mass density, $\Omega_0$, so measuring it can in principle give an estimate of this important cosmological parameter. We are able to measure the effect and distinguish between its strength for high and low values of $\Omega_0$. Unlike other techniques which utilize redshift distortions, one of our methods is relatively insensitive to bias. In one approach, we use path lengths between contour crossings of the density field. The other is based upon percolation. We have found both methods to be successful in quantifying the effect and distinguishing between values of $\Omega_0$. However, only the path lengths method exhibits low sensitivity to bias.

Analyses of peculiar velocity surveys face several challenges, including low signal--to--noise in individual velocity measurements and the presence of small--scale, nonlinear flows. This is the second in a series of papers in which we describe a new method of overcoming these problems by using data compression as a filter with which to separate large--scale, linear flows from small--scale noise that can bias results. We demonstrate the effectiveness of our method using realistic catalogs of galaxy velocities drawn from N--body simulations. Our tests show that a likelihood analysis of simulated catalogs that uses all of the information contained in the peculiar velocities results in a bias in the estimation of the power spectrum shape parameter $\Gamma$ and amplitude $\beta$, and that our method of analysis effectively removes this bias. We expect that this new method will cause peculiar velocity surveys to re--emerge as a useful tool to determine cosmological parameters.

We discuss the tests of nonGaussianity in CMB maps using morphological statistics known as Minkowski functionals. As an example we test degree-scale cosmic microwave background (CMB) anisotropy for Gaussianity by studying the \qmask map that was obtained from combining the QMAP and Saskatoon data. We compute seven morphological functions $M_i(\dt)$, $i=1,...,7$: six \mf and the number of regions $N_c$ at a hundred $\dt$ levels. We also introduce a new parameterization of the morphological functions $M_i(A)$ in terms of the total area $A$ of the excursion set. We show that the latter considerably decorrelates the morphological statistics and makes them more robust because they are less sensitive to the measurements at extreme levels. We compare these results with those from 1000 Gaussian Monte Carlo maps with the same sky coverage, noise properties and power spectrum, and conclude that the \qmask map is neither a very typical nor a very exceptional realization of a Gaussian field. At least about 20% of the 1000 Gaussian Monte Carlo maps differ more than the \qmask map from the mean morphological parameters of the Gaussian fields.

We present a new method for the analysis of peculiar velocity surveys which removes contributions to velocities from small scale, nonlinear velocity modes while retaining information about large scale motions. Our method utilizes Karhunen--Lo\`eve methods of data compression to construct a set of moments out of the velocities which are minimally sensitive to small scale power. The set of moments are then used in a likelihood analysis. We develop criteria for the selection of moments, as well as a statistic to quantify the overall sensitivity of a set of moments to small scale power. Although we discuss our method in the context of peculiar velocity surveys, it may also prove useful in other situations where data filtering is required.

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

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.

The motions of galaxies can be used to constrain the cosmological density parameter Omega and the clustering amplitude of matter on large scales. The mean relative velocity of galaxy pairs, estimated from the Mark III survey, indicates that Omega = 0.35 +0.35/-0.25. If the clustering of galaxies is unbiased on large scales, Omega = 0.35 +/- 0.15, so that an unbiased Einstein-de Sitter model (Omega = 1) is inconsistent with the data.

It is well known that estimating the pairwise velocity of galaxies, v_{12}, from the redshift space galaxy correlation function is difficult because this method is highly sensitive to the assumed model of the pairwise velocity dispersion. Here we propose an alternative method to estimate v_{12} directly from peculiar velocity samples, which contain redshift-independent distances as well as galaxy redshifts. In contrast to other dynamical measures which determine beta = sigma_8 x Omega^{0.6}, our method can provide an estimate of (sigma_8)^2 x Omega^{0.6} for a range of sigma_8 (here Omega is the cosmological mass density parameter while sigma_8 is the standard normalization parameter for the spectrum of matter density fluctuations). We demonstrate how to measure this quantity from realistic catalogues.

We compare the statistical properties of structures normal and transverse to the line of sight which appear in observational data from redshift surveys. We present a statistic which can quantify this effect in a conceptually different way from standard analyses of distortions of the power-spectrum or correlation function. From tests with N-body experiments, we argue that this statistic represents a new, more direct and potentially powerful diagnostic of the cosmological density parameter Omega.

We suggest novel statistics for the CMB maps that are sensitive to non-Gaussian features. These statistics are natural generalizations of the geometrical and topological methods that have been already used in cosmology such as the cumulative distribution function and genus. We compute the distribution functions of the Partial Minkowski Functionals for the excursion set above or bellow a constant temperature threshold. Minkowski Functionals are additive and are translationally and rotationally invariant. Thus, they can be used for patchy and/or incomplete coverage. The technique is highly efficient computationally (it requires only O(N) operations, where N is the number of pixels per one threshold level). Further, it allows to split large data sets into smaller subsets. The full advantage of these statistics can be obtained only on very large data sets. We apply it to the 4-year DMR COBE data corrected for the Galaxy contamination as an illustration of the technique.

The mean pairwise velocity of galaxies has traditionally been estimated from the redshift space galaxy correlation function. This method is notorious for being highly sensitive to the assumed model of the pairwise velocity dispersion. Here we propose an alternative method to estimate the streaming velocity directly from peculiar velocity samples, which contain redshift-independent distances as well as galaxy redshifts. This method can provide an estimate of $\Omega^{0.6}\sigma_8^2$ for a range of $\sigma_8$ where $\Omega$ is the cosmological density parameter, while $\sigma_8$ is the standard normalization for the power spectrum of density fluctuations. We demonstrate how to measure this quantity from realistic catalogues and identify the main sources of bias and errors

We compare the statistical properties of structures normal and transverse to the line of sight which appear in theoretical N-body simulations of structure formation, and seem also to be present in observational data from redshift surveys. We present a statistic which can quantify this effect in a conceptually different way from standard analyses of distortions of the power-spectrum or correlation function. From tests with $N$--body experiments, we argue that this statistic represents a new and potentially powerful diagnostic of the cosmological density parameter, $\Omega_0$.

We propose a new strategy to probe the power spectrum on large scales using galaxy peculiar velocities. We explore the properties of surveys that cover only two small fields in opposing directions on the sky. Surveys of this type have several advantages over those that attempt to cover the entire sky; in particular, by concentrating galaxies in narrow cones these surveys are able to achieve the density needed to measure several moments of the velocity field with only a modest number of objects, even for surveys designed to probe scales $\gtwid 100$\hmpc. We construct mock surveys with this geometry and analyze them in terms of the three moments to which they are most sensitive. We calculate window functions for these moments and construct a $\chi^2$ statistic which can be used to put constraints on the power spectrum. In order to explore the sensitivity of these surveys, we calculate the expectation values of the moments and their associated measurement noise as a function of the survey parameters such as density and depth and for several popular models of structure formation. In addition, we have studied how well these kind of surveys can distinguish between different power spectra and found that, for the same number of objects, cone surveys are as good or better than full-sky surveys in distinguishing between popular cosmological models. We find that a survey with $200-300$ galaxy peculiar velocities with distance errors of 15% in two cones with opening angle of $\sim 10^\circ$ could put significant constraints on the power spectrum on scales of $100-300$\hmpc, where few other constraints exist.

We study the implications of a recent estimate of the bulk flow of a sample of galaxies containing supernovae type Ia by Riess, Press, and Kirshner. We find that their results are quite consistent with power spectra from several currently popular models of structure formation, but that the sample is as yet too sparse to put significant constraints on the power spectrum. We compare this new result with that of Lauer and Postman, with which there is apparent disagreement. We find that for the power spectra we consider, the difference in window functions between the two samples used for the measurements results in a low level of expected correlation between the estimated bulk flows. We calculate a $\chi^2$ for the two measurements taken together and find that their lack of agreement tends to disfavor spectra with excessive power on large scales, but not at a level sufficient to rule them out. A sample consisting of other SN type Ia's found in the Asiago catalog is used to study how the sensitivity of the method used by RPK will improve with increasing sample size. We conclude that the local group motion should be able to be determined with a sample of $\sim 100$ SN Ia light curve shapes.

The large scale velocity field was sampled recently by two independent methods: the supernovae type Ia light curve shapes (Riess, Press \& Kirshner) and the Abell Cluster Catalog brightest cluster galaxy metric luminosity (Lauer \& Postman). The results of these investigations seem to be contradictory. I present an analysis of these samples, compare them and investigate whether standard structure formation models and other deep surveys are compatible with them. I also make suggestions as to how to improve the samples so we can actually resolve the bulk flow vectors. I show that although these two samples seem to cover the same volume, their window functions are sufficiently different so that they are only weakly correlated. Further, since both samples are sparse, they are noise dominated and in order to improve the signal to noise they need to either increase their sample size (RPK) or decrease the measurement errors (LP) significantly.

We calculate the theoretical expectation for the bulk motion of a large scale survey of the type recently carried out by Lauer and Postman. Included are the effects of survey geometry, errors in the distance measurements, clustering properties of the sample, and different assumed power spectra. We consider the power spectrum calculated from the IRAS--QDOT survey, as well as spectra from hot $+$ cold and standard cold dark matter models. We find that sparse sampling and clustering can lead to an unexpectedly large bulk flow, even in a very deep survey. Our results suggest that the expected bulk motion is inconsistent with that reported by Lauer and Postman at the $90-94\%$ confidence level.

We follow a formalism presented by Kaiser to calculate the variance of bulk flows in large scale surveys. We apply the formalism to a mock survey of Abell clusters \'a la Lauer \& Postman and find the variance in the expected bulk velocities in a universe with CDM, MDM and IRAS--QDOT power spectra. We calculate the velocity variance as a function of the 1--D velocity dispersion of the clusters and the size of the survey.

A brief introduction to the gauge invariant classical and quantum theory of cosmological perturbations is given. The formalism is applied to inflationary Universe models and yields a consistent and unified description of the generation and evolution of fluctuations. A general formula for the amplitude of cosmological perturbations in inflationary cosmology is derived.

We develop a general method for power spectrum analysis of three dimensional redshift surveys. We present rigorous analytical estimates for the statistical uncertainty in the power and we are able to derive a rigorous optimal weighting scheme under the reasonable (and largely empirically verified) assumption that the long wavelength Fourier components are Gaussian distributed. We apply the formalism to the updated 1-in-6 QDOT IRAS redshift survey, and compare our results to data from other probes: APM angular correlations; the CfA and the Berkeley 1.2Jy IRAS redshift surveys. Our results bear out and further quantify the impression from e.g.\ counts-in-cells analysis that there is extra power on large scales as compared to the standard CDM model with $\Omega h\simeq 0.5$. We apply likelihood analysis using the CDM spectrum with $\Omega h$ as a free parameter as a phenomenological family of models; we find the best fitting parameters in redshift space and transform the results to real space. Finally, we calculate the distribution of the estimated long wavelength power. This agrees remarkably well with the exponential distribution expected for Gaussian fluctuations, even out to powers of ten times the mean. Our results thus reveal no trace of periodicity or other non-Gaussian behavior.

We study the formation of structure for a universe that undergoes a recent loitering phase. We compare the nonlinear mass distribution to that in a standard, matter dominated cosmology. The statistical aspects of the clustered matter are found to be robust to changes in the expansion law, an exception being that the peculiar velocities are lower by a factor of $\sim 3$ in the loitering model. Further, in the loitering scenario, nonlinear growth of perturbation occurs more recently ($z\sim 3-5$) than in the matter dominated case. Differences in the high redshift appearances of the two models will result but observable consequences depend critically on the chosen form, onset and duration of the loitering phase.