We use N-body simulations combined with semi-analytic models to compute the clustering properties of modeled galaxies at z~3, and confront these predictions with the clustering properties of the observed population of Lyman-break galaxies (LBGs). Several scenarios for the nature of LBGs are explored, which may be broadly categorized into models in which high-redshift star formation is driven by collisional starbursts and those in which quiescent star formation dominates. For each model, we make predictions for the LBG overdensity distribution, the variance of counts-in-cells, the correlation length, and close pair statistics. Models which assume a one-to-one relationship between massive dark-matter halos and galaxies are disfavored by close pair statistics, as are models in which colliding halos are associated with galaxies in a simplified way. However, when modeling of gas consumption and star formation is included using a semi-analytic treatment, the quiescent and collisional starburst models predict similar clustering properties and none of these models can be ruled out based on the available clustering data. None of the ``realistic'' models predict a strong dependence of clustering amplitude on the luminosity threshold of the sample, in apparent conflict with some observational results.

Little is known about the statistics of gravitationally lensed quasars at large (7''-30'') image separations, which probe masses on the scale of galaxy clusters. We have carried out a survey for gravitationally-lensed objects, among sources in the FIRST 20cm radio survey that have unresolved optical counterparts in the digitizations of the Palomar Observatory Sky Survey. From the statistics of ongoing surveys that search for quasars among FIRST sources, we estimate that there are about 9100 quasars in this source sample, making this one of the largest lensing surveys to date. Using broad-band imaging, we have isolated all objects with double radio components separated by 5''-30'', that have unresolved optical counterparts with similar BVI colours. Our criteria for similar colours conservatively allow for observational error and for colour variations due to time delays between lensed images. Spectroscopy of these candidates shows that none of the pairs are lensed quasars. This sets an upper limit (95% confidence) on the lensing fraction in this survey of 3.3x10^-4, assuming 9100 quasars. Although the source redshift distribution is poorly known, a rough calculation of the expected lensing frequency and the detection efficiencies and biases suggests that simple theoretical expectations are of the same order of magnitude as our observational upper limit. Our procedure is novel in that our exhaustive search for lensed objects does not require prior identification of the quasars in the sample as such. Characterization of the FIRST-selected quasar population will enable using our result to constrain quantitatively the mass properties of clusters.

We present a new method to determine the total mass-to-neutral gas ratio in Lyman-limits systems. The method exploits the relation between the neutral hydrogen column density and the magnification of background sources due to the weak gravitational lensing that these systems induce. Because weak lensing does not provide a direct measure of mass, one must use this relation in a statistical sense to solve for the average mass-to-gas ratio and its distribution. We use a detailed mock catalog of quasars (sources) and Lyman-limit systems (lenses) to demonstrate the applicability of this approach through our ability to recover the parameter. This mock catalog also allows us to check for systematics in the method and to sketch its limitations. For a universal constant mass-to-gas ratio and a sample of N quasars, we obtain an unbiased estimate of its value with 95% confidence limits (independent of its actual value) of +/- 140 {10^5/N)^0.5.

[Abridged] We study the angular-momentum profiles of a statistical sample of halos drawn from a high-resolution N-body simulation of the LCDM cosmology. We find that the cumulative mass distribution of specific angular momentum, j, in a halo of mass Mv is well fit by a universal function, M(<j) = Mv \mu j/(j_0+j). This profile is defined by one shape parameter (\mu or j_0) in addition to the global spin parameter \lambda. It follows a power-law over most of the mass, and flattens at large j, with the flattening more pronounced for small values of \mu. Compared to a uniform sphere in solid-body rotation, most halos have a higher fraction of their mass in the low- and high-j tails of the distribution. The spatial distribution of angular momentum in halos tends to be cylindrical and is well-aligned within each halo for ~80% of the halos. We investigate two ideas for the origin of this profile. The first is based on a revised version of linear tidal-torque theory combined with extended Press-Schechter mass accretion, and the second focuses on j transport in minor mergers. Finally, we briefly explore implications of the M(<j) profile on the formation of galactic disks assuming that j is conserved during an adiabatic baryonic infall. The implied gas density profile deviates from an exponential disk, with a higher density at small radii and a tail extending to large radii. The steep central density profiles may imply disk scale lengths that are smaller than observed. This is reminiscent of the "angular-momentum problem" seen in hydrodynamic simulations, even though we have assumed perfect j conservation. A possible solution is to associate the central excesses with bulge components and the outer regions with extended gaseous disks.

We summarize recent developments in the study of the origin of halo spin profiles and preliminary implications on disk formation. The specific angular-momentum distributions within halos in N-body simulations match a universal shape, M(<j) \propto j/(j_0+j). It is characterized by a power law over most of the mass, and one shape parameter in addition to the spin parameter \lambda. The angular momentum tends to be aligned throughout the halo and of cylindrical symmetry. Even if angular momentum is conserved during baryonic infall, the resultant disk density profile is predicted to deviate from exponential, with a denser core and an extended tail. A slightly corrected version of the scaling relation due to linear tidal-torque theory is used to explain the origin of a typical power-law profile in shells, j(M) \propto M^s with s \gsim 1. While linear theory crudely predicts the amplitudes of halo spins, it is not a good predictor of their directions. Independently, mergers of halos are found to produce a similar profile due to j transfer from the orbit to the product halo via dynamical friction and tidal stripping. The halo spin is correlated with having a recent major merger, though this correlation is weakened by mass loss. These two effects are due to a correlation between the spins of neighboring halos and their orbit, leading to prograde mergers.

[Abridged] We analyze the exchange of dark matter between halos, subhalos, and their environments in a high-resolution cosmological N-body simulation of a Lambda CDM cosmology. At each analyzed redshift z we divide the dark matter particles into 4 components: (i) isolated galactic halos, (ii) subhalos, (iii) the diffuse medium of group and cluster halos, and (iv) the background outside of virialized halos. We follow the time evolution of the mass distribution and flows between these components and provide fitting functions for the exchange rates. We use our derived exchange rates to gauge the importance of metal redistribution in the universe due solely to gravity-induced interactions. The diffuse metallicity in clusters is predicted to be ~40% that in isolated galaxies (~55% of groups) at z=0, and should be lower only slightly by z=1, consistent with observations. The metallicity of the diffuse media in poor groups is expected to be lower by a factor of 5 by z~2, in agreement with the observed metallicity of damped Ly$\alpha$ systems. The metallicity of the background IGM is predicted to be (1-3)x10^{-4} that of z=0 clusters, also consistent with observations. The agreement of predicted and observed trends indicates that gravitational interaction alone may play an important role in metal enrichment of the intra-cluster and intergalactic media.

[Abridged] The interaction rates of dark-matter halos and subhalos, are computed using high-resolution cosmological N-body simulations of the Lambda-CDM model. Although the number fraction of subhalos of mass $>2x10^{11}\hsolmass$ is only ~10%, we find that the interaction rate of such subhalos is relatively high because they reside in high density environments. At low redshift, the subhalo collisions dominate the total collision rate, and even at z=3 they are involved in more than 30% of all collisions. About 40% of the "major" mergers are between subhalos. Therefore subhalo interactions must be incorporated in models of structure formation. We find that the collision rate between halos in physical density units, is $\propto (1+z)^\delta $, with $\delta = 3-4$, in agreement with earlier simulations and most observational data. We test previous analytic estimates of the interaction rates of subhalos based on statistical models, which could be very inaccurate because of the small number of subhalos and the variation of conditions within small host halos. We find that, while such statistical estimates may severely overestimate the rate within hosts $< 10^{13} \hsolmass$, typical of high redshifts, they are valid for larger hosts regardless of the number of subhalos in them. We find the Makino & Hut (1997) estimate of the subhalo merger rate to be valid for hosts $ > 6x10^{11}\hsolmass$ at all redshifts. The collision rate between subhalos and the central object of their host halo is approximated relatively well using the timescale for dynamical friction in circular orbits. This approximation fails in ~40% of the cases, partly because of deviations from the assumption of circular orbits and partly because of the invalidity of the assumption that the subhalo mass is negligible.

We study dark-matter halo density profiles in a high-resolution N-body simulation of an LCDM cosmology. Our statistical sample contains ~5000 haloes in the range 10^{11}-10^{14} M_sun. The profiles are parameterized by an NFW form with two parameters, an inner radius r_s and a virial radius r_v, and we define the halo concentration c_v = r_v/r_s. First, we find that, for a given halo mass, the redshift dependence of the median concentration is c_v ~ 1/(1+z), corresponding to a roughly constant r_s with redshift. We present an improved analytic treatment of halo formation that fits the measured relations between halo parameters and their redshift dependence. The implications are that high-redshift galaxies are predicted to be more extended and dimmer than expected before. Second, we find that the scatter in log(c_v) is ~0.18, corresponding to a scatter in maximum rotation velocities of dV/V ~ 0.12. We discuss implications for modelling the Tully-Fisher relation, which has a smaller reported intrinsic scatter. Third, haloes in dense environments tend to be more concentrated than isolated haloes. These results suggest that c_v is an essential parameter for the theory of galaxy modelling, and we briefly discuss implications for the universality of the Tully-Fisher relation, the formation of low surface brightness galaxies, and the origin of the Hubble sequence.

A combined sample of 79 high and low redshift supernovae Ia (SNe) is used to set constraints on the degree of anisotropy in the Universe out to $z\simeq1$. First we derive the global most probable values of matter density $\Omega_M $, the cosmological constant $\Omega_\Lambda $, and the Hubble constant $H_0$, and find them to be consistent with the published results from the two data sets of Riess et al. 1998 (R98) and Perlmutter et al. 1999 (P99). We then examine the Hubble diagram (HD, i.e., the luminosity-redshift relation) in different directions on the sky by utilising spherical harmonic expansion. In particular, via the analysis of the dipole anisotropy, we divide the sky into the two hemispheres that yield the most discrepant of the three cosmological parameters, and the scatter $\chi^2_{\rm HD}$ in each case. The most discrepant values roughly move along the locus $-4\Omega_M +3 \Omega_{\Lambda} = 1$ (cf. P99), but by no more than $\Delta \approx 2.5$ along this line. For a perfect FRW universe, Monte Carlo realizations that mimic the current set of SNe yield values higher than the measured $\Delta$ in $\sim 1/5$ of the cases. We discuss implications for the validity of the Cosmological Principle, and possible calibration problems in the SNe data sets.

We study the distribution functions of mass and circular velocity for dark matter halos in N-body simulations of the $\Lambda$CDM cosmology, addressing redshift and environmental dependence. The dynamical range enables us to resolve subhalos and distinguish them from "distinct" halos. The mass function is compared to analytic models, and is used to derive the more observationally relevant circular velocity function. The distribution functions in the velocity range 100--500 km/s are well fit by a power-law with two parameters, slope and amplitude. We present the parameter dependence on redshift and provide useful fitting formulae. The amplitudes of the mass functions decrease with z, but, contrary to naive expectation, the comoving density of halos of a fixed velocity ~200 km/s actually increases out to z=5. This is because high-z halos are denser, so a fixed velocity corresponds to a smaller mass. The slope of the velocity function at z=0 is as steep as ~ -4, and the mass and velocity functions of distinct halos steepen with increasing z, while the functions of subhalos do not steepen with $z$, and become even flatter at z>2. A simple observable prediction is that the slope of the velocity function of isolated galaxies is steeper than that of galaxies in groups by as much as unity, reflecting the density biasing of high-velocity halos. We confirm that the Press-Schechter approximation typically overestimates the halo mass function by a factor of ~2, while modified approximations provide improved predictions.

We investigate the clustering properties of high-redshift galaxies within three competing scenarios for assigning luminous galaxies to dark matter halos from N-body simulations: a one galaxy per massive halo model, a quiescent star formation model, and a collisional starburst model. We compare these models to observations of Lyman-Break galaxies at z~3$ With current data and the simple statistic used here, one cannot rule out any of these models, but we see potential for finding distinguishing features using statistics that are sensitive to the tails of the distribution, and statistics based on the number of multiple galaxies per halo, which we explore in an ongoing study.

Using high resolution N-body simulations with hydrodynamics and star formation, we investigate interactions and the resulting starbursts in galaxies with properties typical of $z\sim 3$. We apply spectral population models to produce mock-HST images, and discuss the observed magnitude, color, and morphological appearance of our simulated galaxies in both the rest-UV and rest-visual bands.

We present a galaxy circular velocity function, Psi(log v), derived from existing luminosity functions and luminosity-velocity relations. Such a velocity function is desirable for several reasons. First, it enables an objective comparison of luminosity functions obtained in different bands and for different galaxy morphologies, with a statistical correction for dust extinction. In addition, the velocity function simplifies comparison of observations with predictions from high-resolution cosmological N-body simulations. We derive velocity functions from five different data sets and find rough agreement among them, but about a factor of 2 variation in amplitude. These velocity functions are then compared with N-body simulations of a LCDM model (corrected for baryonic infall) in order to demonstrate both the utility and current limitations of this approach. The number density of dark matter halos and the slope of the velocity function near v_*, the circular velocity corresponding to an ~L_* spiral galaxy, are found to be comparable to that of observed galaxies. The primary sources of uncertainty in construction of Psi(log v) from observations and N-body simulations are discussed and explanations are suggected to account for these discrepancies.

Lyman break galaxies (LBGs) at $z \sim 3$ exhibit number densities and clustering similar to local $L_*$ galaxies; however, their star formation rates (SFRs) are much higher. We explore the scenario in which LBGs are starburst galaxies triggered by collisions, and thus provide an explanation for these key properties. The relative paucity of starburst galaxies at low redshift can be attributed to a much higher collision rate in the past. We use high-resolution cosmological $N$-body simulations and a hierarchical halo finder to estimate the galaxy collision rate as a function of time in the popular $\Lambda$CDM cosmological model. We find that bright collisional starbursts are frequent enough to account for most of the high-$z$ ($2.5-4.5$) LBGs. Moreover, many of the objects are of relatively small mass, but they cluster about large-mass halos. They therefore exhibit strong clustering, similar to that observed and stronger than that of the relevant massive halos.

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.

We evaluate the accuracy of semi-analytic merger-trees by comparing them with the merging histories of dark-matter halos in N-body simulations, focusing on the joint distribution of the number of progenitors and their masses. We first confirm that the halo mass function as predicted directly by the Press-Schechter (PS) model deviates from the simulations by up to 50% depending on the mass scale and redshift, while the means of the projected distributions of progenitor number and mass for a halo of a given mass are more accurately predicted by the Extended PS model. We then use the full merger trees to study the joint distribution as a function of redshift and parent-halo mass. We find that while the deviation of the mean quantities due to the inaccuracy of the Extended PS model partly propagates into the higher moments of the distribution, the merger-tree procedure does not introduce a significant additional source of error. In particular, certain properties of the merging history such as the mass ratio of the progenitors and the total accretion rate are reproduced quite accurately for galaxy sized halos ($\sim 10^{12}\msun$), and less so for larger masses. We conclude that although there could be $\sim 50%$ deviations in the absolute numbers and masses of progenitors and in the higher order moment of these distributions, the relative properties of progenitors for a given halo are reproduced fairly well by the merger trees. They can thus provide a useful framework for modelling galaxy formation once the above-mentioned limitations are taken into account.

We investigate several approaches for constructing Monte Carlo realizations of the merging history of virialized dark matter halos (``merger trees'') using the extended Press-Schechter formalism. We describe several unsuccessful methods in order to illustrate some of the difficult aspects of this problem. We develop a practical method that leads to the reconstruction of mean quantities such as the conditional and overall mass functions as given by the Press-Schechter model. This method is convenient, computationally efficient, and works for any power spectrum or background cosmology. In addition, we investigate statistics that describe the distribution of the number of progenitors and their masses as a function of redshift.

We propose a method to remove the mass sheet degeneracy that arises when the mass of galaxy clusters is inferred from gravitational shear. The method utilizes high-redshift standard candles that undergo weak lensing. Natural candidates for such standard candles are type Ia supernovae (SN Ia). When corrected with the light-curve shape (LCS), the peak magnitude of SN Ia provides a standard candle with an uncertainty in apparent magnitude of $\Delta m\simeq 0.1-0.2$. Gravitational magnification of a background SN Ia by an intervening cluster would cause a mismatch between the observed SN Ia peak magnitude compared to that expected from its LCS and redshift. The average detection rate for SN Ia with a significant mismatch of $\ge2\Delta m$ behind a cluster at $z\simeq0.05-0.15$ is about $1-2$ supernovae per cluster per year at $J,I,R\lesssim25-26$. Since SNe are point-like sources for a limited period, they can experience significant microlensing by MACHOs in the intracluster medium. Microlensing events caused by MACHOs of $\sim10^{-4} M_\odot$ are expected to have time scales similar to that of the SN light curve. Both the magnification curve by a MACHO and the light curve of a SN Ia have characteristic shapes that allow to separate them. Microlensing events due to MACHOs of smaller mass can unambiguously be identified in the SN light curve if the latter is continuously monitored. The average number of identifiable microlensing events per nearby cluster ($z\lesssim0.05$) per year is $\sim 0.02 (f/0.01)$, where $f$ is the fraction of the cluster mass in MACHOs of masses $10^{-7} < M_{macho}/M_\odot < 10^{-4}$.

(Abridged) This paper introduces the formalism which connects between rotation measure ($\RM$) measurements for extragalactic sources and the cosmological magnetic field power spectrum. It is shown that the amplitude and shape of the cosmological magnetic field power spectrum can be constrained by using a few hundred radio sources, for which Faraday $\RM$s are available. This constraint is of the form $B_{rms} \simless 1 \times [2.6\times10^{-7} cm^{-3}/ \bar n_b] h $ nano-Gauss (nG) on $\sim 10-50 \hmpc$ scales. The constraint is superior to and supersedes any other constraint which come from either CMB fluctuations, Baryonic nucleosyn thesis, or the first two multipoles of the magnetic field expansion. Demonstration of the ability to detect such magnetic fields, using Bayesian statistics, is carried out by constructing simulations of the field and mimicking observations. This procedure also provides error estimates for the derived quantities. The two main noise contributions due to the Galactic RM and the internal RM are treated in a statistical way. For power indices -1\leq n \leq 1 in a flat cosmology (Omega_m=1) we estimate the signal-to-noise ratio, Q, for limits on the magnetic field B_{rms} on ~50 h^{-1}Mpc scale. Employing one patch of a few square degrees on the sky with source number density n_{src}, an approximate estimate yields Q\simeq 3 \times (B_{rms}/1 {nG})(n_{src}/50 {deg}^{-2}) (2.6\times10^{-7} {cm}^{-3}/ \bar n_b) h $. An all sky coverage, with much sparser, but carefully tailored sample of ~500 sources, yields Q \simeq 1 with the same scaling. An ideal combination of small densely sampled patches and sparse all-sky coverage yields Q\simeq 3 with better constraints for the power index. All of these estimates are corroborated by the simulations.

The projected gravitational potential of galaxy clusters is reflected in both their X-ray emission and their imprint on the images of background sources due to their gravitational lensing effects. Since these projections of the potential are weighted differently along the line-of-sight, we propose a method to combine them and remove the degeneracy between two cases: (i) a cluster consisting of a single potential well, or (ii) an apparent cluster composed of several potential wells projected onto each other. We demonstrate with simulated data of potential models that this method indeed allows to significantly distinguish multiple from single clusters. The confidence limit for this distinction depends on the mass ratio between the clusters. It ranges from ~ 15 sigma for mass ratio 1:1 to ~ 4 sigma for mass ratio 1:6. Furthermore, the method reconstructs the correct cluster mass, the correct mass ratio of the two clusters, and the correct scale radii with typical fractional accuracies of a few percent at 3 sigma confidence. As an aside, our method allows to accurately determine gas fractions in clusters, also with 3 sigma fractional accuracies of order a few percent. We argue that our method provides an alternative to the commonly used beta-fit technique, and yields more reliable results in a broader range of cases.

The power spectrum (PS) of mass density fluctuations, independent of `biasing', is estimated from the Mark III catalog of peculiar velocities using Bayesian statistics. A parametric model is assumed for the PS, and the free parameters are determined by maximizing the probability of the model given the data. The method has been tested using detailed mock catalogs. It has been applied to generalized CDM models with and without COBE normalization. The robust result for all the models is a relatively high PS, with $P(k) \Omega^{1.2} = (4.8 \pm 1.5) \times 10^3 (Mpc/h)^3$ at $k=0.1 h/Mpc$. An extrapolation to smaller scales using the different CDM models yields $\sigma_8 \Omega^{0.6} = 0.88 \pm 0.15$. The peak is weakly constrained to the range $0.02 \leq k \leq 0.06 h/Mpc$. These results are consistent with a direct computation of the PS (Kolatt & Dekel 1996). When compared to galaxy-density surveys, the implied values for $\beta$ ($\equiv \Omega^{0.6}/b$) are of order unity to within 25%. The parameters of the COBE-normalized, flat CDM model are confined by a 90% likelihood contour of the sort $\Omega h_{50}^\mu n^\nu = 0.8 \pm 0.2$, where $\mu = 1.3$ and $\nu = 3.4, 2.0$ for models with and without tensor fluctuations respectively. For open CDM the powers are $\mu = 0.95$ and $\nu = 1.4$ (no tensor fluctuations). A $\Gamma$-shape model free of COBE normalization yields only a weak constraint: $\Gamma = 0.4 \pm 0.2$.

We compare Tully-Fisher (TF) data for 838 galaxies within cz=3000 km/sec from the Mark III catalog to the peculiar velocity and density fields predicted from the 1.2 Jy IRAS redshift survey. Our goal is to test the relation between the galaxy density and velocity fields predicted by gravitational instability theory and linear biasing, and thereby to estimate $\beta_I = \Omega^{0.6}/b_I,$ where $b_I$ is the linear bias parameter for IRAS galaxies. Adopting the IRAS velocity and density fields as a prior model, we maximize the likelihood of the raw TF observables, taking into account the full range of selection effects and properly treating triple-valued zones in the redshift-distance relation. Extensive tests with realistic simulated galaxy catalogs demonstrate that the method produces unbiased estimates of $\beta_I$ and its error. When we apply the method to the real data, we model the presence of a small but significant velocity quadrupole residual (~3.3% of Hubble flow), which we argue is due to density fluctuations incompletely sampled by IRAS. The method then yields a maximum likelihood estimate $\beta_I=0.49\pm 0.07$ (1-sigma error). We discuss the constraints on $\Omega$ and biasing that follow if we assume a COBE-normalized CDM power spectrum. Our model also yields the 1-D noise noise in the velocity field, including IRAS prediction errors, which we find to be be 125 +/- 20 km/sec.

The power spectrum (PS) of mass density fluctuations, in the range $0.05 \leq k \leq 0.2 \ihmpc$, is derived from the Mark III catalog of peculiar velocities of galaxies, independent of ``biasing". It is computed from the density field as recovered by POTENT with Gaussian smoothing of $12\hmpc$, within a sphere of radius $\sim 60\hmpc$ about the Local Group. The density is weighted inversely by the square of the errors. The PS is corrected for the effects of smoothing, random errors, sparse sampling and finite volume using mock catalogs that mimic in detail the Mark III catalog and the dynamics of our cosmological neighborhood. The mock catalogs are also used for error analysis. The PS at $k = 0.1 \ihmpc$ (for example) is $(4.6 \pm 1.4) \times 10^3 \Omega^{-1.2} \3hmpc$, with a local logarithmic slope of $-1.45 \pm 0.5$. An itegration yields $\sigma_8 \Omega^{0.6} \simeq 0.7-0.8$, depending on where the PS peak is. Direct comparisons of the mass PS with the galaxy PS derived from sky and redshift surveys show a similarity in shape and yield for $\beta \equiv \Omega^{0.6}/b$ values in the range $0.77-1.21$, with a typical error of $\pm 0.1$ per galaxy sample. A comparison of the mass PS at $\sim 100\hmpc$ with the large-angle CMB fluctuations by COBE provides constraints on cosmological parameters and on the slope of the initial PS ($n$). The ``standard" CDM model is marginally rejected at the $\sim 2\sigma$ level, while each of the following modifications leads to a good fit: $n \lsim 1$, $\Omega_\nu \sim 0.3$, or $\Omega \lsim 1$. Values of $\Omega$ as low as $\sim 0.2$ are ruled out with high confidence (independent of $\Lambda$).

We find that Gamma-ray bursts (GRBs) at the 3B catalog are correlated (at 95\% confidence level) with Abell clusters. This is the first known correlation of GRBs with any other astronomical population. It confirms the cosmological origin of GRBs. Comparison of the rich clusters auto-correlation with the cross-correlation found here suggests that \sim 26 \pm 15\% of an accurate (\delta<2.3 deg.) position GRBs sub-sample members are located within 600 \hmpc.

The power spectrum (PS) of {\it mass} density fluctuations, independent of ``biasing", is estimated from the Mark3 Catalog of Peculiar Velocities of galaxies using Bayesian statistics. A parametric model is assumed for the PS, and the free parameters are determined by maximizing the posterior probability of the model given the data. The method has been tested using detailed mock catalogs. We use generalized CDM models with and without COBE normalization. The robust result for all the models tested is $P(k) \Omega^{1.2} = (4.1 \pm 0.7)x10^3 (Mpc h^{-1})^3$ at $k=0.1 h Mpc^{-1}$, with the peak constrained to the range $0.03 < k < 0.06 h Mpc^{-1}$. It is consistent with a direct computation of the PS (Kolatt \& Dekel 1996). When compared to galaxy-density surveys, the implied values for $\beta$ ($\equiv \Omega^{0.6}/b$) are of order unity to within 25\%. A $\Gamma$-shape model, free of COBE normalization, is constrained by the velocity data to $\Gamma=0.5\pm 0.15$ and $\sigma_8 \Omega^{0.6}=0.85\pm 0.1$. Within the family of COBE-normalized CDM models, the best tilted model ($\Omega=1$, $n \approx 0.84 h_{50}^{-0.65}$) and the best open model ($n=1$, $\Omega \approx 0.75 h_{50}^{-1.0}$) are more likely than the best $\Lambda$ model ($n=1$, $\Lambda=1-\Omega$, $\Omega\approx 0.70 h_{50}^{-1.2}$). The most likely CDM model with $\Omega\leq 1$ is found to be of $\Omega=1$, $h\approx 0.75$, and a {\it tilted} spectrum of $n= 0.8\pm 0.02$ with tensor fluctuations.The tightest constraint for the tilted-$\Lambda$ models is of the sort $\Omega h_{50}^{1.2} n^\nu = 0.7 \pm 0.08$, with $\nu=3.8$ and $1.85$ with and without tensor fluctuations respectively.

Precise determinations of the image positions in quad gravitational lenses using VLBI can be used to measure the transverse velocity of the lens galaxy and the observer. The typical proper motions are $\mu$as yr$^{-1}$, so the time scale to measure the motion is ten years. By measuring the dipole of the proper motions in an ensemble of lenses we can set limits on the deviation of the inertial frame defined by the lenses from that defined by the CMB dipole and estimate the Hubble constant. The residual proper motions after subtracting the dipole probe the evolution of peculiar velocities with redshift and can be used to estimate the density parameter $\Omega_0$. For $N$ lenses, VLBI measurement accuracies of $\sigma_\theta$, and a baseline of $T$ years, we estimate that the 2$\sigma$ limit on the rms peculiar velocity of the lens galaxies is $3100 (\sigma_\theta/10\mu\{as})({yrs}/T)/N^{1/2} \kms$, and that the time required for the 2--$\sigma$ limit to reach the level of the local rms peculiar velocity $v_{0,rms}$ is approximately $10 N^{-1/2} (v_{0,rms}/600\kms)(\sigma_\theta/10\mu as)$ years. For a ten year baseline and $N=10$ lenses we expect the 1$\sigma$ limit on the misalignment with the CMB dipole to be $\Delta \theta=20^{\circ}$ or equivalently to obtain an upper limit of $\Delta H_0 /H_0 < 0.34 $.

We describe the construction of an N-body simulation that mimics the true velocity and mass-density fields in a box of side 256\hmpc about the Local Group , and the production of mock catalogs that mimic in detail current catalogs of redshifts and peculiar velocities. Our main purpose is to provide a tool for developing and testing reconstruction methods, but the different components of the method can be used on their own in other applications. The initial conditions are based on the IRAS 1.2Jy redshift survey, assuming that galaxies trace mass and \Omega=1. A density field smoothed is recovered from the redshift survey, using quasi-linear theory and a power-preserving filter. The corresponding potential field is traced back to the linear regime using the Zel'dovich-Bernoulli equation. Small-scale power is added by means of constrained realization to mimic fluctuations on galactic scales. The gravitating system is evolved forward in time with a PM code of 2\hmpc resolution. The result reproduces the real dynamical structures on large scales and the statistical properties of the structure down to galactic scales. "Galaxies" are identified via a linear biasing scheme (b=1.35) and they are divided into "spirals" and "ellipticals" to obey Dressler's morphology-density relation. The galaxies are assigned internal-velocity parameters and absolute magnitudes scattered about an assumed mean Tully-Fisher relation. They are then "observed" as magnitude-limited samples, trying to mimic the selection criteria of the data sets constituting the Mark III catalog of peculiar velocities. Artificial IRAS 1.2Jy redshift surveys are also compiled. The simulations and mock catalogs will be made available electronically as bench marks for testing reconstruction methods.

The one-point probability distribution function (pdf) is computed for the $25\hmpc$-smoothed density field of rich clusters of galaxies in the Abell/\aco\ catalogs. The observed pdf is compared to the pdf s drawn similarly from mock catalogs of clusters in cosmological simulations of Gaussian and several non-Gaussian initial conditions. Several statistics allow significant rejection of the non-Gaussian models tested here, and fail to reject the Gaussian model. A comparison with the predictions of second-order perturbation theory and a log-normal model for cdm Gaussian initial conditions yield a linear biasing factor $b_c/b_o \simeq 4$ and $b_c/b_o\simeq 3.7$ for $R\geq 0$ clusters.

We compare the burst distribution of the new (2B) BATSE catalogue to a cosmological distribution. We find that the distribution is insensitive to cosmological parameters such as Omega and Lambda and to the width of the bursts luminosity function. The maximal red shift of the long bursts is ~2.1 (assuming no evolution) while Zm(long) of the short bursts is significantly lower Zm(short) =~ 0.5 In agreement with this relatively nearby origin of the short burst we find an indication that these bursts are correlated ( >=2 sigma level at 10 degrees) with Abell clusters. This is the first known correlation of the bursts with any other astrophysical population and if confirmed by further observations it will provides additional evidence for the cosmological origin of those bursts.

We map the smoothed mass-density distribution in the Galactic zone of avoidance (zoa), within $6000\kms$ of the Local Group, using potent reconstruction from peculiar velocities of galaxies. The interpolation into the zoa is based on the assumed gravitational nature of the velocity field implying a potential flow. The main dynamical features at a distance r=4000 km/s are (a) the peak of the Great Attractor connecting Centaurus and Pavo at l=330, (b) a moderate bridge connecting Perseus-Pisces and Cepheus at l= 140, and (c) an extension of a large void from the southern Galactic hemisphere into the zoa near the direction of Puppis, l=220-270. We find a strong correlation between the mass density and the IRAS and optical galaxy density at $b=+-20, which indicates that the main dynamical features in the zoa should also be seen in galaxy surveys through the Galactic plane. The gravitational acceleration at the Local Group, based on the mass distribution out to 6000 km/s, is strongly affected by the mass distribution in the zoa: its direction changes by 31 when the |b| < 20 zoa is included, bringing it to within 4 +- 19 of the CMB dipole.

We test the hypothesis that the velocity field derived from Tully-Fisher measurements of spiral galaxies, and that derived independently from Dn-sigma measurements of ellipticals and S0s, are noisy versions of the same underlying velocity field. The radial velocity fields are derived using tensor Gaussian smoothing of radius 1200 km/s. They are compared at grid points near which the sampling by both types of galaxies is proper. This requirement defines a volume of ~(50 Mpc/h)^3, containing ~10 independent subvolumes, mostly limited by the available ellipticals. The two fields are compared using a correlation statistic, whose distribution is determined via Monte-Carlo simulations. We find that the data is consistent with the hypothesis, at the 10% level. We demonstrate that the failure to reject the correlation is not just a result of the errors being big, by using the same method to rule out complete independence between the fields at the 99.8% level. The zero points of the two distance indicators are matched by maximizing the correlation between the two velocity fields. There is a marginal hint that the ellipticals tend to stream slower than the spirals by ~8%. The correlation reinforced here is consistent with the common working hypotheses that (a) the derived large-scale velocity field is real, (b) it has a gravitational origin, and (c) the large-scale velocities of spirals and ellipticals are hardly biased relative to each other. On the other hand, it does not rule out any alternative to gravity where objects of all types obtain similar large-scale velocities.