We present extensive tests of the Fast Action Method (FAM) for recovering the past orbits of mass tracers in an expanding universe from their redshift-space coordinates at the present epoch. The tests focus on the reconstruction of present-day peculiar velocities using mock catalogs extracted from high resolution $N$-body simulations. The method allows for a self-consistent treatment of redshift-space distortions by direct minimization of a modified action for a cosmological gravitating system. When applied to ideal, volume limited catalogs, FAM recovers unbiased peculiar velocities with a 1-D, 1\sigma error of ~220 km/s, if velocities are smoothed on a scale of 5 Mpc/h. Alternatively, when no smoothing is applied, FAM predicts nearly unbiased velocities for objects residing outside the highest density regions. In this second case the 1\sigma$error decreases to a level of ~150 km/s. The correlation properties of the peculiar velocity fields are also correctly recovered on scales larger than 5 Mpc/h. Similar results are obtained when FAM is applied to flux limited catalogs mimicking the IRAS PSCz survey. In this case FAM reconstructs peculiar velocities with similar intrinsic random errors, while velocity-velocity correlation properties are well reproduced beyond scales of ~8 Mpc/h. We also show that FAM provides better velocity predictions than other, competing methods based on linear theory or Zel'dovich approximation. These results indicate that FAM can be successfully applied to presently available galaxy redshift surveys such as IRAS PSCz.

We simulate the formation and evolution of the local galaxy population starting from initial conditions with a smoothed linear density field which matches that derived from the IRAS 1.2 Jy galaxy survey. Our simulations track the formation and evolution of all dark matter haloes more massive than 10e+11 solar masses out to a distance of 8000 km/s from the Milky Way. We implement prescriptions similar to those of Kauffmann et al. (1999) to follow the assembly and evolution of the galaxies within these haloes. We focus on two variants of the CDM cosmology: an LCDM and a tCDM model. Galaxy formation in each is adjusted to reproduce the I-band Tully-Fisher relation of Giovanelli et al. (1997). We compare the present-day luminosity functions, colours, morphology and spatial distribution of our simulated galaxies with those of the real local population, in particular with the Updated Zwicky Catalog, with the IRAS PSCz redshift survey, and with individual local clusters such as Coma, Virgo and Perseus. We also use the simulations to study the clustering bias between the dark matter and galaxies of differing type. Although some significant discrepancies remain, our simulations recover the observed intrinsic properties and the observed spatial distribution of local galaxies reasonably well. They can thus be used to calibrate methods which use the observed local galaxy population to estimate the cosmic density parameter or to draw conclusions about the mechanisms of galaxy formation. To facilitate such work, we publically release our z=0 galaxy catalogues, together with the underlying mass distribution.

We allow for nonlinear effects in the likelihood analysis of galaxy peculiar velocities, and obtain ~35%-lower values for the cosmological density parameter Om and the amplitude of mass-density fluctuations. The power spectrum in the linear regime is assumed to be a flat LCDM model (h=0.65, n=1, COBE) with only Om as a free parameter. Since the likelihood is driven by the nonlinear regime, we "break" the power spectrum at k_b=0.2 h/Mpc and fit a power law at k>k_b. This allows for independent matching of the nonlinear behavior and an unbiased fit in the linear regime. The analysis assumes Gaussian fluctuations and errors, and a linear relation between velocity and density. Tests using proper mock catalogs demonstrate a reduced bias and a better fit. We find for the Mark3 and SFI data Om_m=0.32+-0.06 and 0.37+-0.09 respectively, with sigma_8*Om^0.6 = 0.49+-0.06 and 0.63+-0.08, in agreement with constraints from other data. The quoted 90% errors include cosmic variance. The improvement in likelihood due to the nonlinear correction is very significant for Mark3 and moderately so for SFI. When allowing deviations from LCDM, we find an indication for a wiggle in the power spectrum: an excess near k=0.05 and a deficiency at k=0.1 (cold flow). This may be related to the wiggle seen in the power spectrum from redshift surveys and the second peak in the CMB anisotropy. A chi^2 test applied to modes of a Principal Component Analysis (PCA) shows that the nonlinear procedure improves the goodness of fit and reduces a spatial gradient of concern in the linear analysis. The PCA allows addressing spatial features of the data and fine-tuning the theoretical and error models. It shows that the models used are appropriate for the cosmological parameter estimation performed. We address the potential for optimal data compression using PCA.

We present a method for decomposing the cosmological velocity field in a given volume into its divergent component due to the density fluctuations inside the volume, and its tidal component due to the matter distribution outside the volume. The input consists of the density and velocity fields that are reconstructed either by POTENT or by Wiener Filter from a survey of peculiar velocities. The tidal field is further decomposed into a bulk velocity and a shear field. The method is applied here to the Mark III data within a sphere of radius 60 Mpc/h about the Local Group, and to the SFI data for comparison. We find that the tidal field contributes about half of the Local-Group velocity with respect to the CMB, with the tidal bulk velocity pointing to within ~ 30 degrees of the CMB dipole. The eigenvector with the largest eigenvalue of the shear tensor is aligned with the tidal bulk velocity to within ~ 40 degrees. The tidal field thus indicates the important dynamical role of a super attractor of mass (2-5) x 10^17 M_sun/h Omega^0.4 at ~ 150 Mpc/h, coinciding with the Shapley Concentration. There is also a hint for the dynamical role of two big voids in the Supergalactic Plane. The results are consistent for the two data sets and the two methods of reconstruction.

We allow for nonlinear effects in the likelihood analysis of peculiar velocities, and obtain ~35%-lower values for the cosmological density parameter and for the amplitude of mass-density fluctuations. The power spectrum in the linear regime is assumed to be of the flat LCDM model (h=0.65, n=1) with only Om_m free. Since the likelihood is driven by the nonlinear regime, we "break" the power spectrum at k_b=0.2 h/Mpc and fit a two-parameter power-law at k>k_b. This allows for an unbiased fit in the linear regime. Tests using improved mock catalogs demonstrate a reduced bias and a better fit. We find for the Mark III and SFI data Om_m=0.35+-0.09$ with sigma_8*Om_m^0.6=0.55+-0.10 (90% errors). When allowing deviations from \lcdm, we find an indication for a wiggle in the power spectrum in the form of an excess near k~0.05 and a deficiency at k~0.1 h/Mpc --- a "cold flow" which may be related to a feature indicated from redshift surveys and the second peak in the CMB anisotropy. A chi^2 test applied to principal modes demonstrates that the nonlinear procedure improves the goodness of fit. The Principal Component Analysis (PCA) helps identifying spatial features of the data and fine-tuning the theoretical and error models. We address the potential for optimal data compression using PCA.

Internal shocks propagating through an ambient radiation field, are subject to a radiative drag that, under certain conditions, can significantly affect their dynamics and, consequently, the evolution of the beaming cone of emission produced behind the shocks. The resultant change of the Doppler factor combined with opacity effects leads to a strong dependence of the variability pattern produced by such systems, specifically, the shape of the light curves and the characteristics of correlated emission, on viewing angle. One implication is that objects oriented at relatively large viewing angles to the observer should exhibit a higher level of activity at high synchrotron frequencies (above the self-absorption frequency) and at gamma-ray energies below the threshold energy to pair production, than at lower (radio/millimeter) frequencies.

We estimate the power spectrum of mass density fluctuations from peculiar velocities of galaxies by applying an improved maximum-likelihood technique to the new all-sky SFI catalog. Parametric models are used for the power spectrum and the errors, and the free parameters are determined by assuming Gaussian velocity fields and errors and maximizing the probability of the data given the model. It has been applied to generalized CDM models with and without COBE normalization. The method has been carefully tested using artificial SFI catalogs. The most likely distance errors are found to be similar to the original error estimates in the SFI data. The general result that is not very sensitive to the prior model used is a relatively high amplitude of the power spectrum. For example, at k=0.1 h/Mpc we find P(k)\Omega^{1.2}=(4.4+/-1.7)X10^3 (Mpc/h)^3. An integral over the power spectrum yields \sigma_8\Omega^{0.6}=0.82+/-0.12. Model-dependent constraints on the cosmological parameters are obtained for families of CDM models. For example, for COBE-normalized \Lambda CDM models (scalar fluctuations only), the maximum-likelihood result can be approximated by \Omega n^{2} h_{60}^{1.3} =0.58+/-0.11. The formal random errors quoted correspond to the 90% confidence level. The total uncertainty, including systematic errors associated with nonlinear effects, may be larger by a factor of ~2. These results are in agreement with an application of a similar method to other data (Mark III).

We present an improved POTENT method for reconstructing the velocity and mass density fields from radial peculiar velocities, test it with mock catalogs, and apply it to the Mark III Catalog. Method improvments: (a) inhomogeneous Malmquist bias is reduced by grouping and corrected in forward or inverse analyses of inferred distances, (b) the smoothing into a radial velocity field is optimized to reduce window and sampling biases, (c) the density is derived from the velocity using an improved nonlinear approximation, and (d) the computational errors are made negligible. The method is tested and optimized using mock catalogs based on an N-body simulation that mimics our cosmological neighborhood, and the remaining errors are evaluated quantitatively. The Mark III catalog, with ~3300 grouped galaxies, allows a reliable reconstruction with fixed Gaussian smoothing of 10-12 Mpc/h out to ~60 Mpc/h. We present maps of the 3D velocity and mass-density fields and the corresponding errors. The typical systematic and random errors in the density fluctuations inside 40 Mpc/h are \pm 0.13 and \pm 0.18. The recovered mass distribution resembles in its gross features the galaxy distribution in redshift surveys and the mass distribution in a similar POTENT analysis of a complementary velocity catalog (SFI), including the Great Attractor, Perseus-Pisces, and the void in between. The reconstruction inside ~40 Mpc/h is not affected much by a revised calibration of the distance indicators (VM2, tailored to match the velocities from the IRAS 1.2Jy redshift survey). The bulk velocity within the sphere of radius 50 Mpc/h about the Local Group is V_50=370 \pm 110 km/s (including systematic errors), and is shown to be mostly generated by external mass fluctuations. With the VM2 calibration, V_50 is reduced to 305 \pm 110 km/s.

The galaxy density field as extracted from the IRAS 1.2 Jy redshift survey is compared to the mass density field as reconstructed by the POTENT method from the Mark III catalog of peculiar velocities. The reconstruction is done with Gaussian smoothing of radius 12 h^{-1}Mpc, and the comparison is carried out within volumes of effective radii 31-46 h^{-1}Mpc, containing approximately 10-26 independent samples. Random and systematic errors are estimated from multiple realizations of mock catalogs drawn from a simulation that mimics the observed density field in the local universe. The relationship between the two density fields is found to be consistent with gravitational instability theory in the mildly nonlinear regime and a linear biasing relation between galaxies and mass. We measure beta = Omega^{0.6}/b_I = 0.89 \pm 0.12 within a volume of effective radius 40 h^{-1}Mpc, where b_I is the IRAS galaxy biasing parameter at 12 h^{-1}Mpc. This result is only weakly dependent on the comparison volume, suggesting that cosmic scatter is no greater than \pm 0.1. These data are thus consistent with Omega=1 and b_I\approx 1. If b_I>0.75, as theoretical models of biasing indicate, then Omega>0.33 at 95% confidence. A comparison with other estimates of beta suggests scale-dependence in the biasing relation for IRAS galaxies.