We investigate the evolution of the galaxy two point correlation function (CF) over a wide redshift range, 0.2 < z < 3. For the first time the systematic analysis covers the redshifts above 1 - 1.5. The catalogue of ~250000 galaxies with i+ < 25 and known photometric redshifts in the Subaru Deep field is used. The galaxies are divided into three luminosity classes and several distance/redshift bins. First, the 2D CF is determined for each luminosity class and distance bin. Calculations are based on the quantitative differences between the surface distributions of galaxy pairs with comparable and distinctly different photometric redshifts. The power law approximation for the CF is used. A limited accuracy of photometric redshifts as compared to the spectroscopic ones has been examined and taken into account. Then, the 3D functions for all the selected luminosities and distances are calculated. The power-law parameters of the CF, the slope and the correlation length, are determined. Both parameters do not show strong variations over the whole investigated redshift range. The slope of the luminous galaxies appears to be consistently steeper than that for the fainter ones. The linear bias factor, b(z), grows systematically with redshift; assuming the local normalization b(0) = 1.1-1.2, the bias reaches 3 - 3.5 at the high redshift limit.

Comparison of peculiar velocities of galaxies with their gravitational accelerations (induced by the density field) is one of the methods to constrain the redshift distortion parameter \beta=(\Omega_m^0.55)/b, where \Omega_m is the non-relativistic matter density parameter and b is the linear bias. In particular, one can use the motion of the Local Group (LG) for that purpose. Its peculiar velocity is known from the dipole component of the cosmic microwave background, whereas its acceleration can be estimated with the use of an all-sky galaxy catalog, from the so-called clustering dipole. At the moment, the biggest dataset of that kind is the Two Micron All Sky Survey Extended Source Catalog (2MASS XSC) containing almost 1 million galaxies and complete up to ~300 Mpc/h. We applied 2MASS data to measure LG acceleration and used two methods to estimate the beta parameter. Both of them yield \beta~0.4 with an error of several per cent, which is the most precise determination of this parameter from the clustering dipole to date.

Sources generating most of the X-ray background (XRB) are dispersed over a wide range of redshifts. Thus, statistical characteristics of the source distribution carry information on matter distribution on very large scales. We test the possibility of detecting the variation in the X-ray source number counts over the celestial sphere. A large number of Chandra pointings spread over both galactic hemispheres are investigated. We searched for all the point-like sources in the soft band of 0.5 - 2 keV and statistically assessed the population of sources below the detection threshold. A homogeneous sample of the number counts at fluxes above ~10^{-15} erg s^{-1} cm^{-2} was constructed for more than 300 ACIS fields. The sources were counted within a circular area of 15 arcmin diameter. The count correlations between overlapping fields were used to assess the accuracy of the computational methods used in the analysis. The average number of sources in the investigated sample amounts to 46 per field. It is shown that the source number counts vary between fields at a level exceeding the fluctuation amplitude expected for the random (Poissonian) distribution. The excess fluctuations are attributed to the cosmic variance generated by the large-scale structures. The rms variations of the source counts due to the cosmic variance within the 15$ arcmin circle reach 8% of the average number counts. An amplitude of the potential correlations of the source counts on angular scales larger than the size of a single pointing remains below the noise level.

There is a long-standing controversy about the convergence of the dipole moment of the galaxy angular distribution (the so-called clustering dipole). We study the growth of the clustering dipole of galaxies as a function of the limiting flux of the sample from the Two Micron All Sky Survey (2MASS). Contrary to some earlier claims, we find that the dipole does not converge before the completeness limit of the 2MASS Extended Source Catalog, i.e. up to 13.5 mag in the near-infrared K_s band (equivalent to an effective distance of 300 Mpc/h). We compare the observed growth of the dipole with the theoretically expected, conditional one (i.e., given the velocity of the Local Group relative to the CMB), for the LambdaCDM power spectrum and cosmological parameters constrained by WMAP. The observed growth turns out to be within 1-sigma confidence level of its theoretical counterpart once the proper observational window of the 2MASS flux-limited catalog is included. For a contrast, if the adopted window is a top-hat, then the predicted dipole grows significantly faster and converges to its final value for a distance of about 300 Mpc/h. By comparing the observational windows, we show that for a given flux limit and a corresponding distance limit, the 2MASS flux-weighted window passes less large-scale signal than the top-hat one. We conclude that the growth of the 2MASS dipole for effective distances greater than 200 Mpc/h is only apparent. On the other hand, for a distance of 80 Mpc/h (mean depth of the 2MASS Redshift Survey) and the LambdaCDM power spectrum, the true dipole is expected to reach only ~80% of its final value. Eventually, since for the window function of 2MASS the predicted growth is consistent with the observed one, we can compare the two to evaluate beta = (Omega_m)^0.55 / b. The result is beta = 0.38+-0.04, which leads to an estimate of the density parameter Omega_m = 0.20+-0.08.

It is widely believed that the cosmological redshift is not a Doppler shift. However, Bunn & Hogg have recently pointed out that to settle properly this problem, one has to transport parallelly the velocity four-vector of a distant galaxy to the observer's position. Performing such a transport along the null geodesic of photons arriving from the galaxy, they found that the cosmological redshift is purely kinematic. Here we argue that one should rather transport the velocity four-vector along the geodesic connecting the points of intersection of the world-lines of the galaxy and the observer with the hypersurface of constant cosmic time. We find that the resulting relation between the transported velocity and the redshift of arriving photons is not given by a relativistic Doppler formula. Instead, for small redshifts it coincides with the well known non-relativistic decomposition of the redshift into a Doppler (kinematic) component and a gravitational one. We perform such a decomposition for arbitrary large redshifts and derive a formula for the kinematic component of the cosmological redshift, valid for any FLRW cosmology. In particular, in a universe with Omega_m = 0.24 and Omega_Lambda = 0.76, a quasar at a redshift 6, at the time of emission of photons reaching us today had the recession velocity v = 0.997c. This can be contrasted with v = 0.96c, had the redshift been entirely kinematic. Thus, for recession velocities of such high-redshift sources, the effect of deceleration of the early Universe clearly prevails over the effect of its relatively recent acceleration. Last but not least, we show that the so-called proper recession velocities of galaxies, commonly used in cosmology, are in fact radial components of the galaxies' four-velocity vectors. As such, they can indeed attain superluminal values, but should not be regarded as real velocities.

We study the growth of the clustering dipole of galaxies from the Two Micron All Sky Survey (2MASS). We find that the dipole does not converge before the completeness limit of the 2MASS Extended Source Catalog, i.e. up to about 300 Mpc/h. We compare the observed growth of the dipole with the theoretically expected, conditional growth for the LambdaCDM power spectrum and cosmological parameters constrained by WMAP. The observed growth turns out to be within 1-sigma confidence level of the theoretical one, once the proper observational window of the 2MASS flux dipole is included. For a contrast, if the adopted window is a top hat, then the predicted dipole grows significantly faster and converges to its final value at a distance of about 200 Mpc/h. We study the difference between the top-hat window and the window for the flux-limited 2MASS survey and we conclude that the growth of the 2MASS dipole at effective distances greater than 200 Mpc/h is only apparent. Eventually, since for the window function of 2MASS the predicted growth is consistent with the observed one, we can compare the two to evaluate beta = (Omega_m)^0.55 / b. The result is beta = 0.38+-0.05, which gives a rough estimate of Omega_m = 0.2+-0.1.

In measurements of the clustering dipole from all-sky surveys, an important problem is the lack of information about galaxy distribution in the so-called Zone of Avoidance (ZoA). The existence of the Local Void (LV) has a systematic effect on these measurements. If the ZoA is randomly filled with mock galaxies, then the calculated acceleration of the Local Group of galaxies (LG) has a spurious component, resulting from the lack of real galaxies in the intersection of the LV with the ZoA. This component affects both the misalignment angle between the clustering dipole and the CMB dipole, and the inferred value of mean matter density Omega_m. We calculate the amplitude of the spurious acceleration acting on the LG due to the LV. Its value depends on the geometry and size of the LV, as well as on its density contrast. However, under the simplest assumption of the LV being spherical and completely empty, within the linear theory this amplitude amounts only to about 45 km/s in units of velocity. The resulting change in the misalignment angle is smaller than 1 degree, and the fractional change in the deduced value of Omega_m is about 5%. Accounting for observationally indicated elongation of the LV and maintaining the maximising assumption of a complete lack of galaxies inside increases these numbers only moderately. Specifically, the amplitude of the spurious acceleration rises to about 60 km/s, the misalignment angle remains still smaller than 1 degree, and the fractional change in the deduced value of Omega_m is enhanced to about 7%. Thus, despite the overall importance of the Local Void for the motion of the Local Group, the influence of the intersection of the LV with the ZoA on measurements of the clustering dipole is found to be only a minor systematic effect.

In two recent papers, Abramowicz et al. claim that the expansion of the Universe can be interpreted only as the expansion of space. In fact, what they really prove is that the cosmological expansion cannot be described in terms of real motions in Minkowski spacetime. However, there is no controversy about this issue. Abramowicz et al. show that in general, the cosmological redshift is not a Doppler shift and they consider this fact as a proof that space expands. Again, nobody believes (perhaps except Milne) that for non-empty universes the origin of the redshift is purely Dopplerian. From the Principle of Equivalence it follows that there must be also a gravitational shift in presence of matter. Indeed, it is well known in cosmology that for small redshifts, the cosmological redshift can be decomposed into a Doppler component and a gravitational component. In a forthcoming paper, we shall perform such a decomposition for arbitrarily large values of the redshift.

We study the cosmic velocity-density relation using the spherical collapse model (SCM) as a proxy to non-linear dynamics. Although the dependence of this relation on cosmological parameters is known to be weak, we retain the density parameter Omega_m in SCM equations, in order to study the limit Omega_m -> 0. We show that in this regime the considered relation is strictly linear, for arbitrary values of the density contrast, on the contrary to some claims in the literature. On the other hand, we confirm that for realistic values of Omega_m the exact relation in the SCM is well approximated by the classic formula of Bernardeau (1992), both for voids (delta<0) and for overdensities up to delta ~ 3. Inspired by this fact, we find further analytic approximations to the relation for the whole range delta from -1 to infinity. Our formula for voids accounts for the weak Omega_m-dependence of their maximal rate of expansion, which for Omega_m < 1 is slightly smaller that 3/2. For positive density contrasts, we find a simple relation div v = 3 H_0 (Omega_m)^(0.6) [ (1+delta)^(1/6) - (1+delta)^(1/2) ], that works very well up to the turn-around (i.e. up to delta ~ 13.5 for Omega_m = 0.25 and neglected Omega_Lambda). Having the same second-order expansion as the formula of Bernardeau, it can be regarded as an extension of the latter for higher density contrasts. Moreover, it gives a better fit to results of cosmological numerical simulations.

We propose a simple way to estimate the parameter beta = Omega_m^(0.6)/b from three-dimensional galaxy surveys. Our method consists in measuring the relation between the cosmological velocity and gravity fields, and thus requires peculiar velocity measurements. The relation is measured *directly in redshift space*, so there is no need to reconstruct the density field in real space. In linear theory, the radial components of the gravity and velocity fields in redshift space are expected to be tightly correlated, with a slope given, in the distant observer approximation, by g / v = (1 + 6 beta / 5 + 3 beta^2 / 7)^(1/2) / beta. We test extensively this relation using controlled numerical experiments based on a cosmological N-body simulation. To perform the measurements, we propose a new and rather simple adaptive interpolation scheme to estimate the velocity and the gravity field on a grid. One of the most striking results is that nonlinear effects, including `fingers of God', affect mainly the tails of the joint probability distribution function (PDF) of the velocity and gravity field: the 1--1.5 sigma region around the maximum of the PDF is *dominated by the linear theory regime*, both in real and redshift space. This is understood explicitly by using the spherical collapse model as a proxy of nonlinear dynamics. Applications of the method to real galaxy catalogs are discussed, including a preliminary investigation on homogeneous (volume limited) `galaxy' samples extracted from the simulation with simple prescriptions based on halo and sub-structure identification, to quantify the effects of the bias between the galaxy and the total matter distibution, and of shot noise (ABRIDGED).

A comparison of the 2MASS flux dipole to the CMB dipole can serve as a method to constrain a combination of the cosmological parameter Omega_m and the luminosity bias of the 2MASS survey. For this constraint to be as tight as possible, it is necessary to maximize the correlation between the two dipoles. This can be achieved by optimizing the survey window through which the flux dipole is measured. Here we explicitly construct such a window for the 2MASS survey. The optimization in essence reduces to excluding from the calculation of the flux dipole galaxies brighter than some limiting magnitude K_min of the near-infrared K_s band. This exclusion mitigates nonlinear effects and shot noise from small scales, which decorrelate the 2MASS dipole from the CMB dipole. Under the assumption of negligible shot noise we find that the optimal value of K_min is about five. Inclusion of shot noise shifts the optimal K_min to larger values. We present an analytical formula for shot noise for the 2MASS flux dipole, to be used in follow-up work with 2MASS data. The misalignment angle between the two dipoles is a sensitive measure of their correlation: the higher the correlation, the smaller the expectation value of the angle. A minimum of the misalignment is thus a sign of the optimal gravity window. We model analytically the distribution function for the misalignment angle and show that the misalignment estimated by Maller et al. is consistent with the assumed underlying model (though it is greater than the expectation value). We predict with about 90% confidence that the misalignment will decrease if 2MASS galaxies brighter than K_min = 5 mag are excluded from the calculation of the flux dipole. This prediction has been indirectly confirmed by the results of Erdogdu et al. (ABRIDGED)

Consider radar ranging of a distant galaxy in a Friedman-Lemaitre cosmological model. In this model the comoving coordinate of the galaxy is constant, hence the equations of null geodesics for photons travelling to the distant galaxy and back imply the following equation: \int_{t_e}^{t_r} dt/a(t) = \int_{t_r}^{t_o} dt/a(t). Here, t_e, t_r and t_o are respectively the times of emission, reflection and observation of the reflected photons, and a(t) is the scale factor. Since the universe is expanding, a(t) is a monotonically increasing function, so the return travel time, t_o - t_r, must be greater than the forward travel time, t_r - t_e. Clearly, space expands, and on their way back, the photons must travel a longer distance! The present paper explains why this argument for the expansion of space is wrong. We argue that, unlike the expansion of the cosmic substratum, the expansion of space is unobservable. We therefore propose to apply to it -- just like to the ether -- Ockham's razor.

In all Friedman models, the cosmological redshift is widely interpreted as a consequence of the general-relativistic phenomenon of EXPANSION OF SPACE. Other commonly believed consequences of this phenomenon are superluminal recession velocities of distant galaxies and the distance to the particle horizon greater than c*t (where t is the age of the Universe), in apparent conflict with special relativity. Here, we study a particular Friedman model: empty universe. This model exhibits both cosmological redshift, superluminal velocities and infinite distance to the horizon. However, we show that the cosmological redshift is there simply a relativistic Doppler shift. Moreover, apparently superluminal velocities and `acausal' distance to the horizon are in fact a direct consequence of special-relativistic phenomenon of time dilation, as well as of the adopted definition of distance in cosmology. There is no conflict with special relativity, whatsoever. In particular, INERTIAL recession velocities are subluminal. Since in the real Universe, sufficiently distant galaxies recede with relativistic velocities, these special-relativistic effects must be at least partly responsible for the cosmological redshift and the aforementioned `superluminalities', commonly attributed to the expansion of space. Let us finish with a question resembling a Buddhism-Zen `koan': in an empty universe, what is expanding?

Based on the magnitude--redshift diagram for the sample of supernovae Ia analysed by Perlmutter et al. (1999), Davis & Lineweaver rule out the special relativistic interpretation of cosmological redshifts at a confidence level of 23 sigma. Here, we critically reassess this result. Special relativity is known to describe the dynamics of an empty universe, by means of the Milne kinematic model. Applying only special-relativistic concepts, we derive the angular diameter distance and the luminosity distance in the Milne model. In particular, in this model we do not use the underlying metric in its Robertson-Walker form, so our exposition is useful for readers without any knowledge of general relativity. We do however, explicitly use the special-relativistic Doppler formula for redshift. We apply the derived luminosity distance to the magnitude--redshift diagram for supernovae Ia of Perlmutter et al. (1999) and show that special relativity fits the data much better than that claimed by Davis & Lineweaver. Specifically, using these data alone, the Milne model is ruled out only at a 2 sigma level. Although not a viable cosmological model, in the context of current research on supernovae Ia it remains a useful reference model when comparing predictions of various cosmological models.

We derive the prediction of the standard model of superluminal radio sources for the apparent transverse velocity of a radio source located at redshift z. The apparent velocity of the source is reduced by a factor of 1 + z compared to that of a similar nearby source. The cause of this reduction is recession of the distant source due to the expansion of the universe. The apparent velocity of a source can be estimated from its redshift and proper motion using the values of the Hubble constant and the mean densities of different energy components in the universe. We derive an expression for the velocity valid for the currently favored cosmological model: a flat universe with a nonzero cosmological constant.

We reexamine likelihood analyzes of the Local Group (LG) acceleration, paying particular attention to nonlinear effects. Under the approximation that the joint distribution of the LG acceleration and velocity is Gaussian, two quantities describing nonlinear effects enter these analyzes. The first one is the coherence function, i.e. the cross-correlation coefficient of the Fourier modes of gravity and velocity fields. The second one is the ratio of velocity power spectrum to gravity power spectrum. To date, in all analyzes of the LG acceleration the second quantity was not accounted for. Extending our previous work, we study both the coherence function and the ratio of the power spectra. With the aid of numerical simulations we obtain expressions for the two as functions of wavevector and \sigma_8. Adopting WMAP's best determination of \sigma_8, we estimate the most likely value of the parameter \beta and its errors. As the observed values of the LG velocity and gravity, we adopt respectively a CMB-based estimate of the LG velocity, and Schmoldt et al.'s (1999) estimate of the LG acceleration from the PSCz catalog. We obtain \beta = 0.66^{+0.21}_{-0.07}; thus our errorbars are significantly smaller than those of Schmoldt et al. This is not surprising, because the coherence function they used greatly overestimates actual decoherence between nonlinear gravity and velocity.

We reexamine likelihood analyses of the Local Group (LG) acceleration, paying particular attention to nonlinear effects. Under the approximation that the joint distribution of the LG acceleration and velocity is Gaussian, two quantities describing nonlinear effects enter these analyses. The first one is the coherence function, i.e. the cross-correlation coefficient of the Fourier modes of gravity and velocity fields. The second one is the ratio of velocity power spectrum to gravity power spectrum. To date, in all analyses of the LG acceleration the second quantity was not accounted for. Extending our previous work, we study both the coherence function and the ratio of the power spectra. With the aid of numerical simulations we obtain expressions for the two as functions of wavevector and sigma_8. Adopting WMAP's best determination of sigma_8, we estimate the most likely value of the parameter beta and its errors. As the observed values of the LG velocity and gravity, we adopt respectively a CMB-based estimate of the LG velocity, and Schmoldt et al.'s (1999) estimate of the LG acceleration from the PSCz catalog. We obtain beta = 0.66^{+0.21}_{-0.07}; thus our errorbars are significantly smaller than those of Schmoldt et al. This is not surprising, because the coherence function they used greatly overestimates actual decoherence between nonlinear gravity and velocity.

We present a numerical study of the relation between the cosmic peculiar velocity field and the gravitational acceleration field. We show that on mildly non-linear scales (4-10 Mpc Gaussian smoothing), the distribution of the Cartesian coordinates of each of these fields is well approximated by a Gaussian. In particular, their kurtoses and negentropies are small compared to those of the velocity divergence and density fields. We find that at these scales the relation between the velocity and gravity field follows linear theory to good accuracy. Specifically, the systematic errors in velocity-velocity comparisons due to assuming the linear model do not exceed 6% in beta. To correct for them, we test various nonlinear estimators of velocity from density. We show that a slight modification of the alpha-formula proposed by Kudlicki et al. yields an estimator which is essentially unbiased and has a small variance.

In maximum-likelihood analyses of the Local Group (LG) acceleration, the object describing nonlinear effects is the coherence function (CF), i.e. the cross-correlation coefficient of the Fourier modes of the velocity and gravity fields. We study the CF both analytically, using perturbation theory, and numerically, using a hydrodynamic code. The dependence of the function on Omega_m and the shape of the power spectrum is very weak. The only cosmological parameter that the CF is strongly sensitive to is the normalization sigma_8 of the underlying density field. Perturbative approximation for the function turns out to be accurate as long as sigma_8 is smaller than about 0.3. For higher normalizations we provide an analytical fit for the CF as a function of sigma_8 and the wavevector. The characteristic decoherence scale which our formula predicts is an order of magnitude smaller than that determined by Strauss et al. This implies that present likelihood constraints on cosmological parameters from analyses of the LG acceleration are significantly tighter than hitherto reported.

In maximum-likelihood analyses of the Local Group (LG) acceleration, the object describing nonlinear effects is the coherence function (CF), i.e. the cross-correlation coefficient of the Fourier modes of the velocity and gravity fields. We study the CF both analytically, using perturbation theory, and numerically, using a hydrodynamic code. The dependence of the function on Omega_m and the shape of the power spectrum is very weak. The only cosmological parameter that the CF is strongly sensitive to is the normalization sigma_8 of the underlying density field. Perturbative approximation for the function turns out to be accurate as long as sigma_8 is smaller than about 0.3. For higher normalizations we provide an analytical fit for the CF as a function of sigma_8 and the wavevector. The characteristic decoherence scale which our formula predicts is an order of magnitude smaller than that determined by Strauss et al. This implies that present likelihood constraints on cosmological parameters from analyses of the LG acceleration are significantly tighter than hitherto reported.

We use numerical simulations to study the relation between the velocity of the Local Group and its gravitational acceleration. This relation serves as a test for the kinematic origin of the CMB dipole and as a method for estimating beta. We calculate the misalignment angle between the two vectors and compare it to the observed value for the PSCz survey. The latter value is near the upper limit of the 95% confidence interval for the angle; therefore, the nonlinear effects are unlikely to be responsible for the whole observed misalignment. We also study the relation between the amplitudes of the LG velocity and gravity vectors. In an Omega=1 Universe, the smoothed gravity of the LG turns out to be a biased low estimator of the LG (unsmoothed) velocity. In an Omega=0.3 Universe, the estimator is biased high. The discussed biases are, however, only a few per cent, thus the linear theory works to good accuracy. The gravity-based estimator of the LG velocity has also a scatter that limits the precision of the estimate of beta in the LG velocity-gravity comparisons. The random error of beta due to nonlinear effects amounts to several per cent.

I derive a nonlinear local relation between the redshift-space density field and the real-space velocity field. The relation accounts for radial character of redshift distortions, and it is not restricted to the limit of the distant observer. Direct comparisons between the observed redshift-space density fields and the real-space velocity fields possess all of the advantages of the conventional redshift-space analyses, while at the same time they are free of their disadvantages. In particular, neither the model-dependent reconstruction of the density field in real space is necessary, nor is the reconstruction of the nonlinear velocity field in redshift space, questionable because of its vorticity at second order. The nonlinear redshift-space velocity field is irrotational only in the distant observer limit, and that limit is not a good approximation for shallow catalogs of peculiar velocities currently available. Unlike the conventional redshift-space comparisons, the comparison proposed here does not have to be restricted to the linear regime. Accounting for nonlinear effects removes one of the sources of bias in the estimate of beta. Moreover, the nonlinear effects break the Omega-bias degeneracy plaguing all analyses based on linear theory.

Given an irrotational (vorticity free) velocity field in real space, we prove that, in the distant observer limit and in the absence of multi-valued zones, the associated velocity field in redshift space is also irrotational. The proof does not rely on any approximation to gravitational dynamics. The result can be particularly useful for the analysis of redshift distortions and for reconstruction methods of cosmological velocity fields from galaxy redshift surveys, in the nonlinear regime. Although the proof is restricted to the distant observer limit, we show that the POTENT method can be modified to derive the full real space velocity field as a function of the redshift space coordinate, thus avoiding spatial Malmquist biases.

We present a numerical study of the cosmic density vs. velocity divergence relation (DVDR) in the mildly non-linear regime. We approximate the dark matter as a non-relativistic pressureless fluid, and solve its equations of motion on a grid fixed in comoving coordinates. Unlike N-body schemes, this method yields directly the volume-averaged velocity field. The results of our simulations are compared with the predictions of the third-order perturbation theory (3PT) for the DVDR. We investigate both the mean `forward' relation (density in terms of velocity divergence) and the mean `inverse' relation (velocity divergence in terms of density), with emphasis on the latter. On scales larger than about 20 megaparsecs, our code recovers the predictions of 3PT remarkably well, significantly better than recent N-body simulations. On scales of a few megaparsecs, the DVDR predicted by 3PT differs slightly from the simulated one. In particular, approximating the inverse DVDR by a third-order polynomial turns out to be a poor fit. We propose a simple analytical description of the inverse relation, which works well for mildly non-linear scales.

I propose to compare the redshift-space density field directly to the REAL-SPACE velocity field. Such a comparison possesses all of the advantages of the conventional redshift-space analyses, while at the same time it is free of their disadvantages. In particular, the model-dependent reconstruction of the density field in real space is unnecessary, and so is the reconstruction of the velocity field in redshift space. The redshift-space velocity field can be reconstructed only at the linear order, because only at this order it is irrotational. Unlike the conventional redshift-space density--velocity comparisons, the comparison proposed here does not have to be restricted to the linear regime. Nonlinear effects can then be used to break the Omega-bias degeneracy plaguing the analyses based on the linear theory. I present a degeneracy-breaking method for the case of nonlinear but local bias.

We present results of the investigations of the statistical properties of a joint density and velocity divergence probability distribution function (PDF) in the mildly non-linear regime. For that purpose we use both perturbation theory results, extended here for a top-hat filter, and numerical simulations. In particular we derive the quantitative (complete as possible up to third order terms) and qualitative predictions for constrained averages and constrained dispersions -- which describe the nonlinearities and the stochasticity properties beyond the linear regime -- and compare them against numerical simulations. We find overall a good agreement for constrained averages; however, the agreement for constrained dispersions is only qualitative. Scaling relations for the Omega-dependence of these quantities are satisfactory reproduced. Guided by our analytical and numerical results, we finally construct a robust phenomenological description of the joint PDF in a closed analytic form. The good agreement of our formula with results of N-body simulations for a number of cosmological parameters provides a sound validation of the presented approach. Our results provide a basis for a potentially powerful tool with which it is possible to analyze galaxy survey data in order to test the gravitational instability paradigm beyond the linear regime and put useful constraints on cosmological parameters. In particular we show how the nonlinearity in the density--velocity relation can be used to break the so-called Omega-bias degeneracy in cosmic density-velocity comparisons.

I derive a second-order local relation between the REDSHIFT-space mass density field and the REAL-space velocity field. This relation can be useful for comparisons between the cosmic density and peculiar velocity fields, for a number of reasons. First, relating the real-space velocity directly to the redshift-space density enables one to avoid the Omega-dependent reconstruction of the density field in real space. Secondly, the reconstruction of the three-dimensional velocity field in redshift space, questionable because of its vorticity, is also unnecessary. Finally, a similar relation between the GALAXY density field and the velocity field offers a way to break the Omega-bias degeneracy in density--velocity comparisons, when combined with an additional measurement of the redshift-space galaxy skewness. I derive the latter relation under the assumption of nonlinear but local bias; accounting for stochasticity of bias is left for further study.

Using third-order perturbation theory, we derive a relation between the mean divergence of the peculiar velocity given density and the density itself. Our calculations assume Gaussian initial conditions and are valid for Gaussian filtering of the evolved density and velocity fields. The mean velocity divergence turns out to be a third-order polynomial in the density contrast. We test the power spectrum dependence of the coefficients of the polynomial for scale-free and standard CDM spectra and find it rather weak. Over scales larger than about 5 megaparsecs, the scatter in the relation is small compared to that introduced by random errors in the observed density and velocity fields. The relation can be useful for recovering the peculiar velocity from the associated density field, and also for non-linear analyses of the anisotropies of structure in redshift surveys.

I derive a relation, both `forward' and `inverse', between the density and the divergence of the peculiar velocity which results from the Zel'dovich approximation. My calculations assume Gaussian initial conditions. The forward relation expresses the density (strictly speaking, the expectation value of the continuity density given the velocity divergence) in terms of the velocity divergence, while the inverse relation expresses the velocity divergence in terms of the density. The predicted scatter in the relations is small, hence the inverse relation is close to, though not identical with, a mathematical inversion of the forward one. The forward relation is equivalent to the well-known `standard' density--velocity relation in the Zel'dovich approximation. The inverse relation, however, is successfully derived for the first time and constitutes a potentially interesting alternative to an inverse relation derived by Chodorowski et al., based on third-order perturbation theory. Specifically, it may better recover the peculiar velocity from the associated density field, when smoothed over scales as small as a few megaparsecs.

I present recent progress in theoretical modelling of cosmological density--velocity relations in the weakly nonlinear regime. The relations are local, based on rigorous perturbation theory and include the effects of smoothing of the density and the velocity fields. For small smoothing scales, they can be improved by slight adjustments based on N-body results. The relations can be useful for density--velocity comparisons, yielding a fair estimate of Omega and offering a method for disentangling Omega and bias.

I derive up to second order in Eulerian perturbation theory a new relation between the weakly nonlinear density and velocity fields. In the case of unsmoothed fields, density at a given point turns out to be a purely local function of the expansion (divergence) and shear of the velocity field. The relation depends on the cosmological parameter Omega, strongly by the factor f(Omega) = Omega^{0.6} and weakly by the factors K(Omega) and C(Omega) proportional to Omega^{-2/63} and Omega^{-1/21} respectively. The Gramann solution is found to be equivalent to the derived relation with the weak Omega-dependence neglected. To make the relation applicable to the real world, I extend it for the case of smoothed fields. The resulting formula, when averaged over shear given divergence, reproduces up to second order the density-velocity divergence relation of Chodorowski & Lokas; however, it has smaller spread. It makes the formula a new attractive local estimator of large-scale density from velocity.

We rigorously derive weakly nonlinear relation between cosmic density and velocity fields up to third order in perturbation theory. The density field is described by the mass density contrast, $\de$. The velocity field is described by the variable $\te$ proportional to the velocity divergence, $\te = - f(\Omega)^{-1} H_0^{-1} \nabla\cdot\bfv$, where $f(\Omega) \simeq \Omega^{0.6}$, $\Omega$ is the cosmological density parameter and $H_0$ is the Hubble constant. Our calculations show that mean $\de$ given $\te$ is a third order polynomial in $\te$, $\lan \de \ran|_{\te} = a_1 \te + a_2 (\te^2 - \s_\te^2) + a_3 \te^3$. This result constitutes an extension of the formula $\lan \de \ran|_{\te} = \te + a_2 (\te^2 - \s_\te^2)$, found by Bernardeau~(1992) which involved second order perturbative solutions. Third order perturbative corrections introduce the cubic term. They also, however, cause the coefficient $a_1$ to depart from unity, in contrast with the linear theory prediction. We compute the values of the coefficients $a_p$ for scale-free power spectra, as well as for standard CDM, for Gaussian smoothing. The coefficients obey a hierarchy $a_3 \ll a_2 \ll a_1$, meaning that the perturbative series converges very fast. Their dependence on $\Omega$ is expected to be very weak. The values of the coefficients for CDM spectrum are in qualitative agreement with the results of N-body simulations by Ganon et al. (1996). The results provide a method for breaking the $\Omega$-bias degeneracy in comparisons of cosmic density and velocity fields such as IRAS-POTENT.

We calculate the kurtosis of a large-scale density field which has undergone weakly non-linear gravitational evolution from arbitrary non-Gaussian initial conditions. It is well known that the weakly evolved {\twelveit skewness} is equal to its initial value plus the term induced by gravity, which scales with the rms density fluctuation in precisely the same way as for Gaussian initial conditions. As in the case of skewness, the evolved {\twelveit kurtosis} is equal to its initial value plus the contribution induced by gravity. The scaling of this induced contribution, however, turns out to be qualitatively different for Gaussian versus non-Gaussian initial conditions. Therefore, measurements of the kurtosis can serve as a powerful discriminating test between the hypotheses of Gaussian and non-Gaussian nature of primordial density fluctuations.

We calculate the cosmological evolution of the 1-point probability distribution function (PDF), using an analytic approximation that combines gravitational perturbation theory with the Edgeworth expansion of the PDF. Our method applies directly to a smoothed mass density field or to the divergence of a smoothed peculiar velocity field, provided that rms fluctuations are small compared to unity on the smoothing scale, and that the primordial fluctuations that seed the growth of structure are Gaussian. We use this `Edgeworth approximation' to compute the evolution of $<\delta|\delta|>$ and $<|delta|>$; these measures are similar to the skewness and kurtosis of the density field, but they are less sensitive to tails of the probability distribution, so they may be more accurately estimated from surveys of limited volume. We compare our analytic calculations to cosmological N-body simulations in order to assess their range of validity. When $\sigma \ll 1$, the numerical simulations and perturbation theory agree precisely, demonstrating that the N-body method can yield accurate results in the regime of weakly non-linear clustering. We show analytically that `biased' galaxy formation preserves the relation $<\delta^3> \propto <\delta^2>^2$ predicted by second-order perturbation theory, provided that the galaxy density is a local function of the underlying mass density.