We explore shapes of clusters and superclusters in the IRAS 1.2 Jy redshift survey with three reconstructions spanning the range \beta = 0.1, 0.5, 1.0, where \beta = \Omega^{0.6}/b, b is the bias factor and \Omega the present value of the dimensionless matter density. Comparing our results to Gaussian randomized reconstructions of the IRAS catalogue, we find structures having both planar and filamentary properties. For \beta = 0.5, 1.0 the largest structures in the survey have a distinct tendency to be filament-like} in general agreement with the results of N-body simulations.

Using an Eulerian perturbative calculation, we show that the distribution of relative pairwise velocities which arises from gravitational instability of Gaussian density fluctuations has asymmetric (skewed) exponential tails. The negative skewness is induced by the negative mean streaming velocity of pairs (the infall prevails over expansion), while the exponential tails arise because the relative pairwise velocity is a number, not volume weighted statistic. The derived probability distribution is compared with N-body simulations and shown to provide a reasonable fit.

We study the distortions induced by peculiar velocities on the redshift-space correlation function of galaxies of different morphological types in the Pisces-Perseus redshift survey. Redshift-space distortions affect early- and late-type galaxies in different ways. In particular, at small separations, the dominant effect comes from virialized cluster cores, where ellipticals are the dominant population. The net result is that a meaningful comparison of the clustering strength of different morphological types can be performed only in real space, i.e., after projecting out the redshift distortions on the two-point correlation function xi(r_p,pi). A power-law fit to the projected function w_p(r_p) on scales smaller than 10/h Mpc gives r_o = 8.35_{-0.76}^{+0.75} /h Mpc, \gamma = 2.05_{-0.08}^{+0.10} for the early-type population, and r_o = 5.55_{-0.45}^{+0.40} /h Mpc, \gamma = 1.73_{-0.08}^{+0.07} for spirals and irregulars. These values are derived for a sample luminosity brighter than M_{Zw} = -19.5. We detect a 25% increase of r_o with luminosity for all types combined, from M_{Zw} = -19 to -20. In the framework of a simple stable-clustering model for the mean streaming of pairs, we estimate sigma_12(1), the one-dimensional pairwise velocity dispersion between 0 and 1 /h Mpc, to be 865^{+250}_{-165} km/s for early-type galaxies and 345^{+95}_{-65} km/s for late types. This latter value should be a fair estimate of the pairwise dispersion for ``field'' galaxies; it is stable with respect to the presence or absence of clusters in the sample, and is consistent with the values found for non-cluster galaxies and IRAS galaxies at similar separations.

We utilise the method of Wiener reconstruction with spherical harmonics and Bessel functions to recover the density and velocity fields from the IRAS 1.2Jy redshift survey. The results are robust to changes in prior parameters and the number of expansion coefficients. Maps of these fields are presented in a variety of projections. Many known structures are observed, including clear confirmations of the clusters N1600 and A3627. The Perseus-Pisces supercluster appears to extend out to roughly 9000km/s, and the reconstruction shows `backside infall' to the Centaurus/Great Attractor region. A qualitative comparison of the reconstructed IRAS gravity field with that from Tully-Fisher peculiar velocity measurements (Mark III) shows reasonable agreement. The Wiener reconstruction of the density field is also the optimal reconstruction (in the minimum variance sense) of any quantity which is linear in the density contrast. We show reconstructions of three such quantities. The misalignment angle between the IRAS and CMB Local Group dipoles is only 13 degrees out to 5000km/s, but increases to 25 degrees out to 20,000km/s. The reconstructed IRAS bulk flow out to 5000km/s is roughly 300km/s, which agrees in amplitude with that derived from the Mark III peculiar velocities (370km/s). However, the two bulk flow vectors deviate by 70 degrees. Finally, moment of inertia analysis shows that the Wiener reconstructed Supergalactic Plane is aligned within 30 degrees of that defined by de Vaucouleurs.

We present percolation analyses of Wiener Reconstructions of the IRAS 1.2 Jy Redshift Survey. There are ten reconstructions of galaxy density fields in real space spanning the range $\beta= 0.1$ to $1.0$, where ${\beta}={\Omega^{0.6}}/b$, $\Omega$ is the present dimensionless density and $b$ is the bias factor. Our method uses the growth of the largest cluster statistic to characterize the topology of a density field, where Gaussian randomized versions of the reconstructions are used as standards for analysis. For the reconstruction volume of radius, $R {\approx} 100 h^{-1} $ Mpc, percolation analysis reveals a slight `meatball' topology for the real space, galaxy distribution of the IRAS survey. cosmology-galaxies:clustering-methods:numerical

We analyze the anisotropy signature expected if the high energy (above $10^{19}$eV) cosmic ray (CR) sources are extra-Galactic and trace the distribution of luminous matter on large scales. We investigate the dependence of the anisotropy on both the relative bias between the CR sources and the galaxy distribution and on the (unknown) intrinsic CR source density. We find that the expected anisotropy associated with the large scale structure (LSS) should be detected once the number of CR events observed above $10^{19}{\rm eV}$ is increased by a factor of $\sim10$. This would require $\sim30$ observation-years with existing experiments, but less then $1$ year with the proposed $\sim5000\ {\rm km}^2$ Auger detectors. We find that the recently reported concentration of the Haverah Park CR events towards the super-galactic plane is not consistent with the known LSS. If real, the Haverah Park result suggests that the CR sources are much more concentrated towards the super-galactic plane than the known LSS. Our results are not sensitive to the number density of CR sources. We show that once the number of detected events is increased by a factor of $\sim10$, the number density would be strongly constrained by considering the probability for having repeating sources.

We examine the clustering of galaxies around a sample of 20 luminous low redshift (z<0.30) quasars observed with the Wide Field Camera-2 on the Hubble Space Telescope. The HST resolution makes possible galaxy identification brighter than V=23.5 and as close as 2'' to the quasar. We find a significant enhancement of galaxies within a projected separation of < 100 kpc/h of the quasars. If we model the qso/galaxy correlation function as a power law with a slope given by the galaxy/galaxy correlation function, we find that the ratio of the qso/galaxy to galaxy/galaxy correlation functions is $3.8\pm 0.8$. The galaxy counts within r<15 kpc/h of the quasars are too high for the density profile to have an appreciable core radius ( > 100 kpc). Our results reinforce the idea that low redshift quasars are located preferentially in groups of 10-20 galaxies rather than in rich clusters. We see no significant difference in the clustering amplitudes derived from radio-loud and radio-quiet subsamples.

We examine the anisotropies in the power spectrum by the mapping of real to redshift space. Using the Zel'dovich approximation, we obtain an analytic expression for the nonlinear redshift space power spectrum in the distant observer limit. For a given unbiased galaxy distribution in redshift space, the anisotropies in the power spectrum depend on the parameter $f(\Omega)\approx \Omega^{0.6}$, where $\Omega$ is the density parameter. We quantify these anisotropies by the ratio, $R$, of the quadrupole to monopole angular moments of the power spectrum. In contrast to linear theory, the Zel'dovich approximation predicts a decline in $R$ with decreasing scale. This departure from linear theory is due to nonlinear dynamics and not a result of incoherent random velocities. The rate of decline depends strongly on $\Omega$ and the initial power spectrum. However, we find a {\it universal} relation between the quantity $R/R_{lin}$ (where $R_{lin}$ the linear theory value of $R$) and the dimensionless variable $k/k_{nl}$, where $k_{nl}$ is a wavenumber determined by the scale of nonlinear structures. The universal relation is in good agreement with a large N-body simulation. This universal relation greatly extends the scales over which redshift distortions can be used as a probe of $\Omega$. A preliminary application to the 1.2 Jy IRAS yields $\Omega\sim 0.4$ if IRAS galaxies are unbiased.

We present some new results from an ongoing investigation of the anisotropy of the two--point correlation function for optically selected galaxies. We have estimated xi(r_p,pi) from the Perseus-Pisces redshift survey, which is now virtually 100\% complete for all morphological types to m_B = 15.5. We detect strong distortions of the iso-correlation contours in xi(r_p,pi). These correspond to a global (`cluster' plus `field') pairwise velocity dispersion at 1 h^{-1} Mpc, sigma(1), which is about a factor of two higher than the canonical value \sim 350 km s^{-1} obtained from the CfA1 survey and recently from the 1.2 Jy IRAS redshift survey. However, the variance of this value is still rather strong even within the relatively large volumes surveyed (typically 90$^\circ$ $\times$ 45$^\circ$ $\times$ 100 h^{-1} Mpc). From the largest volume explorable (126 h^{-1} Mpc depth), the best bet for the `cosmic' value of sigma(1) seems to be between 600 and 700 km s^{-1}. However, for a hopefully robust determination of its value we will have to wait for the analyses of the next generation of larger {\tenbf optical} redshift surveys (such as, e.g., the now completed ESP or the forthcoming Sloan Digital Sky Survey).

We present the redshift data for a survey of galaxies selected from the data base of the Infrared Astronomical Satellite (IRAS). This survey extends the 1.936 Jy sample of Strauss et al. (1992) from a flux limit of 1.936 Jy at 60 microns to 1.2 Jy. The survey extension consists of 3920 sources in the flux interval 1.2 - 1.936 Jy, of which 2663 are galaxies with measured redshifts. Fourteen objects (0.52%) do not have redshifts. The survey covers 87.6% of the sky. The data for the complete 1.2 Jy survey (the data presented here in addition to that of Strauss \etal 1992) may be obtained in a machine-readable form from the National Space Science Data Center and from the anonymous ftp site given above.

The relation between the galaxy correlation function in real and redshift-space is derived in the linear regime by an appropriate averaging of the joint probability distribution of density and velocity. The derivation recovers the familiar linear theory result on large scales but has the advantage of clearly revealing the dependence of the redshift distortions on the underlying peculiar velocity field; streaming motions give rise to distortions of ${\cal O}(\Omega^{0.6}/b)$ while variations in the anisotropic velocity dispersion yield terms of order ${\cal O}(\Omega^{1.2}/b^2)$. This probabilistic derivation of the redshift-space correlation function is similar in spirit to the derivation of the commonly used ``streaming'' model, in which the distortions are given by a convolution of the real-space correlation function with a velocity distribution function. The streaming model is often used to model the redshift-space correlation function on small, highly non-linear, scales. There have been claims in the literature, however, that the streaming model is not valid in the linear regime. Our analysis confirms this claim, but we show that the streaming model can be made consistent with linear theory {\it provided} that the model for the streaming has the functional form predicted by linear theory and that velocity distribution is chosen to be a Gaussian with the correct linear theory dispersion.

We measure the anisotropy of the redshift-space power spectrum in the 1.2-Jy and QDOT redshift surveys of IRAS-selected galaxies. On large scales, this anisotropy is caused by coherent peculiar motions, and gravitational instability theory predicts a distortion of the power spectrum that depends only on the ratio $\beta \equiv f(\Omega)/b \approx \Omega^{0.6}/b$, where Omega is the cosmological density parameter and $b$ is the bias parameter. On small scales, the distortion is dominated by the random velocity dispersion in non-linear structures. We fit the observed anisotropy with an analytic model that incorporates two parameters, beta, and a small-scale velocity dispersion sigma_v. Tests on N-body simulations show that this model recovers beta quite accurately on the scales accessible to the existing IRAS redshift surveys. Applying our procedure to the 1.2-Jy and QDOT surveys, we find beta=0.52 +/- 0.13 and beta=0.54 +/- 0.3, respectively. These results imply Omega approximately 0.35 if galaxies trace mass, or a bias factor of about 2 if Omega=1.

The formalism of Wiener filtering is developed here for the purpose of reconstructing the large scale structure of the universe from noisy, sparse and incomplete data. The method is based on a linear minimum variance solution, given data and an assumed \prior model which specifies the covariance matrix of the field to be reconstructed. While earlier applications of the Wiener filter have focused on estimation, namely suppressing the noise in the measured quantities, we extend the method here to perform both prediction and dynamical reconstruction. The Wiener filter is used to predict the values of unmeasured quantities, such as the density field in un-sampled regions of space, or to deconvolve blurred data. The method is developed, within the context of linear gravitational instability theory, to perform dynamical reconstruction of one field which is dynamically related to some other observed field. This is the case, for example, in the reconstruction of the real space galaxy distribution from its redshift distribution When the field to be reconstructed is a Gaussian random field, such as the primordial perturbation field predicted by the canonical model of cosmology, the Wiener filter can be pushed to its fullest potential. In such a case the Wiener estimator coincides with the Bayesian estimator designed to maximize the {\it posterior} probability. The Wiener filter can be also derived by assuming a quadratic regularization function, in analogy with the `Maximum Entropy' method. The mean field obtained by the minimal variance solution can be supplemented with constrained realizations of the Gaussian field to

We present a new method for recovering the cosmological density, velocity, and potential fields from all-sky redshift catalogues. The method is based on an expansion of the fields in orthogonal radial (Bessel) and angular (spherical harmonic) functions. In this coordinate system, peculiar velocities introduce a coupling of the radial harmonics describing the density field in redshift space but leave the angular modes unaffected. In the harmonic transform space, this radial coupling is described by a distortion matrix which can be computed analytically within the context of linear theory; the redshift space harmonics can then be converted to their real space values by inversion of this matrix. Statistical noise is mitigated by regularizing the matrix inversion with a Wiener filter. The method yields a minimum variance estimate of the density field in real space. In this coordinate system, the minimum variance harmonics of the peculiar velocity and potential fields are related to those of the density field by simple linear transformations. Tests of the method with simulations of a CDM universe and comparison with previously proposed methods demonstrate it to be a very promising new reconstruction method for the local density and velocity field. A first application to the 1.2 Jy IRAS redshift survey is presented.

We derive an optimal linear filter to suppress the noise from the COBE DMR sky maps for a given power spectrum. We then apply the filter to the first-year DMR data, after removing pixels within $20^\circ$ of the Galactic plane from the data. The filtered data have uncertainties 12 times smaller than the noise level of the raw data. We use the formalism of constrained realizations of Gaussian random fields to assess the uncertainty in the filtered sky maps. In addition to improving the signal-to-noise ratio of the map as a whole, these techniques allow us to recover some information about the CMB anisotropy in the missing Galactic plane region. From these maps we are able to determine which hot and cold spots in the data are statistically significant, and which may have been produced by noise. In addition, the filtered maps can be used for comparison with other experiments on similar angular scales.

The analysis of whole-sky galaxy surveys commonly suffers from the problems of shot-noise and incomplete sky coverage (e.g. at the Zone of Avoidance). The orthogonal set of spherical harmonics is utilized here to expand the observed galaxy distribution. We show that in the framework of Bayesian statistics and Gaussian random fields the $4 \pi$ harmonics can be recovered and the shot-noise can be removed, giving the most probable picture of the underlying density field. The correction factor from observed to reconstructed harmonics turns out to be the well-known Wiener filter (the ratio of signal to signal+noise), which is also derived by requiring minimum variance. We apply the method to the projected 1.2 Jy IRAS survey. The reconstruction confirms the connectivity of the Supergalactic Plane across the Galactic Plane (at Galactic longitude $l \sim 135^o$ and $l \sim 315^o$) and the Puppis cluster behind the Galactic Plane ($ l \sim 240^o$). The method can be extended to 3-D in both real and redshift space, and applied to other cosmic phenomena such as the COBE Microwave Background maps.

We examine the nature of galaxy clustering in redshift space using a method based on an expansion of the galaxian density field in Spherical Harmonics and linear theory. We derive a compact and self-consistent expression for the distortion when applied to flux limited redshift surveys. The amplitude of the distortion is controlled by the combination of the density and bias parameters, $\beta\equiv\Omega_\circ^{0.6}/b$ as well as the shape of the real space power spectrum, $P(k)$ (characterized by a shape parameter $\Gamma$), and its normalization, $\sigma_8$; we exploit this fact to derive a maximum likelihood estimator for $\beta$, $\Gamma$, and $\sigma_8$. We check our formalism using $N$-body simulations and demonstrate it provides an unbiased estimate of $\beta$ when the amplitude and shape of the galaxy power spectrum is known. Application of the technique to the 1.2 Jy \iras\ redshift survey yields $\beta =0.94\pm 0.17$ and $\Gamma=0.17\pm0.05$ (1-$\sigma$) when $\sigma_8$ is held fixed at its best value as determined from the real space correlation function. Allowing $\sigma_8$ to be a free parameter, we find $\beta=0.47\pm 0.25$, $\Gamma=0.15\pm0.05$, and $\sigma_8=0.81\pm0.06$.

We present a formalism for analysing redshift distortions based on a spherical harmonic expansion of the density field. We use a maximum likelihood estimator for the combination of density and bias parameters, $\Omega^0.6/b$. We test the method with $N$-body simulations and apply it to the 1.2 Jy IRAS redshift survey.

We examine the effect of redshift space distortions on the galaxy two-point correlation function $\xi(r_p,\pi)$ as a function of separations parallel ($r_p$) and perpendicular ($\pi$) to the line of sight. We find that the relative velocity dispersion of pairs of IRAS galaxies is $\sigma(r)= 317^{+40}_{-49}$ \kms at $r=1 \mpc$, consistent with previous estimates derived from optically selected galaxy catalogues. Unfortunately, the use of this result to estimate $\Omega$ via the Cosmic Virial Theorem is thwarted by large systematic uncertainties. We also fit for the mean relative streaming velocity of pairs, $v_{12}(r)$, which describes the growth of fluctuations on both linear and nonlinear scales. We find that $v_{12}(r) = 167^{+99}_{-67}$ \kms at $r=4$\mpc, so that on average, approximately half the Hubble expansion velocity of pairs at this separation is canceled by infall. At $r=10$\mpc, the amplitude of the streaming is lower and $v_{12}(r) = 109^{+64}_{-47} \kms$. Linear perturbation theory then implies that $\Omega^{0.6}/b =0.45^{+0.27}_{-0.18}$ on scales $\sim 10-15 \mpc$. The amplitude of $v_{12}(r)$ is sensitive to the assumed shape of $\sigma(r)$; if the latter deviates substantially from a virialized form on small scale, our best fit amplitude of $v_{12}(r)$ can deviate by a factor of two.

The peculiar velocities of galaxies distort the pattern of galaxy clustering in redshift space, making the redshift space power spectrum anisotropic. In the linear regime, the strength of this distortion depends only on the ratio $\beta \equiv f(\Omega)/b \approx \Omega^{0.6}/b$, where $\Omega$ is the cosmological density parameter and $b$ is the bias parameter. We derive a linear theory estimator for $\beta$ based on the harmonic moments of the redshift space power spectrum. Using N-body simulations, we examine the impact of non-linear gravitational clustering on the power spectrum anisotropy and on our $\beta$-estimator. Non-linear effects can be important out to wavelengths $\lambda \sim 50$Mpc/h or larger; in most cases, they lower the quadrupole moment of the power spectrum and thereby depress the estimate of $\beta$ below the true value. With a sufficiently large redshift survey, the scaling of non-linear effects may allow separate determinations of $\Omega$ and $b$. We describe a practical technique for measuring the anisotropy of the power spectrum from galaxy redshift surveys, and we test the technique on mock catalogues drawn from the N-body simulations. Preliminary application of our methods to the 1.2 Jy IRAS galaxy survey yields $\beta_{est} \sim 0.3-0.4 $ at wavelengths $\lambda \sim 30-40$Mpc/h . Non-linear effects remain important at these scales, so this estimate of $\beta$ is probably lower than the true value.

We present analyses of the two-point correlation function derived from an all-sky redshift survey of 5313 galaxies extracted from the Infrared Astronomical Satellite (IRAS) database. The redshift space correlation function \xi(s) is well described by a power law, $\xi(s) = (s/4.53h^{-1}{\rm Mpc})^{-1.28}$, on scales \simlt 20\mpc; on larger scales \xi(s) drops below the extension of this power law. We examine the effect of redshift space distortions on the correlation function and compute the full two dimensional correlation function \xi(r_p,\pi). From this, we derive the real space correlation function, which is well described by $\xi(r) = (r/3.76h^{-1}{\rm Mpc})^{-1.66}$ on scales \simlt 20\mpc. The derived correlation functions are found to be consistent with previous determinations in the literature, and seem to show more power on large scales than predicted by the standard Cold Dark Matter (CDM) model. Comparison of the derived \xi(r) with the correlation function of optical galaxies implies an optical to IRAS bias ratio of $b_O/b_I = 1.38\pm 0.12$ on a scale of $\sim 8$\mpc . The variances in cubical cells inferred from \xi(s) appear discrepant with the previously reported results of Efstathiou et al. (1990).