We present a study of the statistical properties of three velocity dispersion and mass estimators, namely biweight, gapper and standard deviation, in the small number of galaxies regime ($N_{\rm gal} \le 75$). Using a set of 73 numerically simulated galaxy clusters, we characterise the statistical bias and the variance for the three estimators, both in the determination of the velocity dispersion and the dynamical mass of the clusters via the $\sigma-M$ relation. The results are used to define a new set of unbiased estimators, that are able to correct for those statistical biases with a minimal increase of the associated variance. The numerical simulations are also used to characterise the impact of velocity segregation in the selection of cluster members, and the impact of using cluster members within different physical radii from the cluster centre. The standard deviation is found to be the lowest variance estimator. The selection of galaxies within the sub-sample of the most massive galaxies in the cluster introduces a $2\,$\% bias in the velocity dispersion estimate when calculated using a quarter of the most massive cluster members. We also find a dependence of the velocity dispersion estimate on the aperture radius as a fraction of $R_{200}$, consistent with previous results. The proposed set of unbiased estimators effectively provides a correction of the velocity dispersion and mass estimates from all those effects in the small number of cluster members regime. This is tested by applying the new estimators to a subset of simulated observations. Although for a single galaxy cluster the statistical and physical effects discussed here are comparable or slightly smaller than the bias introduced by interlopers, they will be of relevance when dealing with ensemble properties and scaling relations for large cluster samples (Abridged).

Using a set of 73 numerically simulated galaxy clusters, we have characterised the statistical and physical biases for three velocity dispersion and mass estimators, namely biweight, gapper and standard deviation, in the small number of galaxies regime ($N_{gal} \leq 75$), both for the determination of the velocity dispersion and the dynamical mass of the clusters via the $\sigma-M$ relation. These results are used to define a new set of unbiased estimators, that are able to correct for those statistical biases. By applying these new estimators to a subset of simulated observations, we show that they can retrieve bias-corrected values for both the mean velocity dispersion and the mean mass.

The one-point probability distribution function (PDF) of the matter density field in the universe is a fundamental property that plays an essential role in cosmology for estimates such as gravitational weak lensing, non-linear clustering, massive production of mock galaxy catalogs, and testing predictions of cosmological models. Here we make a comprehensive analysis of the dark matter PDF using a suite of 7000 N-body simulations that covers a wide range of numerical and cosmological parameters. We find that the PDF has a simple shape: it declines with density as a power-law P~rho**(-2), which is exponentially suppressed on both small and large densities. The proposed double-exponential approximation provides an accurate fit to all our N-body results for small filtering scales R< 5Mpc/h with rms density fluctuations sigma>1. In combination with the spherical infall model that works well for small fluctuations sigma<1, the PDF is now approximated with just few percent errors over the range of twelve orders of magnitude -- a remarkable example of precision cosmology. We find that at 5-10% level the PDF explicitly depends on redshift (at fixed sigma) and on cosmological density parameter Omega_m. We test different existing analytical approximations and find that the often used log-normal approximation is always 3-5 times less accurate than either the double-exponential approximation or the spherical infall model.

We compare the predicted conditional mass function (CMF) of dark matter halos from two theoretical prescriptions against numerical N-body simulations, both in overdense and underdense regions and at different Eulerian scales ranging from $5$ to $30\,h^{-1}\,$Mpc. In particular, we consider in detail a locally-implemented rescaling of the unconditional mass function (UMF) already discussed in the literature, and also a generalization of the standard rescaling method described in the extended Press-Schechter formalism. First, we test the consistency of these two rescalings by verifying the normalization of the CMF at different scales, and showing that none of the proposed cases provides a normalized CMF. In order to satisfy the normalization condition, we include a modification in the rescaling procedure. After this modification, the resulting CMF generally provides a better description of numerical results. We finally present an analytical fit to the ratio between the CMF and the UMF (also known as the matter-to-halo bias function) in underdense regions, which could be of special interest to speed-up the computation of the halo abundance when studying void statistics. In this case, the CMF prescription based on the locally-implemented rescaling provides a slightly better description of the numerical results when compared to the standard rescaling.

We report the detection, with a high level of confidence, of coherent outflows around voids found in the seventh data release of the Sloan Digital Sky Survey (SDSS DR7). In particular, we developed a robust <|cos theta|> statistical test to quantify the strength of redshift-space distortions (RSD) associated with extended coherent velocity fields. We consistently find that the vector that joins void centers with galaxies that lie in shells around them is more likely to be perpendicular to the line-of-sight than parallel to it. This effect is clear evidence for the existence of outflows in the vicinity of voids. We show that the RSD exist for a wide range of void radius and shell thickness, but they are more evident in the largest voids in our sample. For instance, we find that the <|cos theta|> for galaxies located in shells within 2 h^-1 Mpc from the edge of voids larger than 15 h^-1 Mpc deviates 3.81sigma from uniformity. The measurements presented in this work provide useful information to constrain cosmological parameters, in particular Omega_m and Sigma_8.

Using a large sample of galaxies from the SDSS-DR7, we have analysed the alignment of disk galaxies around cosmic voids. We have constructed a complete sample of cosmic voids (devoid of galaxies brighter than M_r-5 log(h) =-20.17) with radii larger than 10 Mpc/h up to redshift 0.12. Disk galaxies in shells around these voids have been used to look for particular alignments between the angular momentum of the galaxies and the radial direction of the voids. We find that disk galaxies around voids larger than >~ 15 Mpc/h within distances not much larger than 5 Mpc/h from the surface of the voids present a significant tendency to have their angular momenta aligned with the void's radial direction with a significance >~ 98.8% against the null hypothesis. The strength of this alignment is dependent on the void's radius and for voids with ~<15 Mpc/h the distribution of the orientation of the galaxies is compatible with a random distribution. Finally, we find that this trend observed in the alignment of galaxies is similar to the one observed for the minor axis of dark matter halos around cosmic voids found in cosmological simulations, suggesting a possible link in the evolution of both components.

We study the concentration of dark matter halos and its evolution in N-body simulations of the standard LCDM cosmology. The results presented in this paper are based on 4 large N-body simulations with about 10 billion particles each: the Millennium-I and II, Bolshoi, and MultiDark simulations. The MultiDark (or BigBolshoi) simulation is introduced in this paper. This suite of simulations with high mass resolution over a large volume allows us to compute with unprecedented accuracy the concentration over a large range of scales (about six orders of magnitude in mass), which constitutes the state-of-the-art of our current knowledge on this basic property of dark matter halos in the LCDM cosmology. We find that there is consistency among the different simulation data sets. We confirm a novel feature for halo concentrations at high redshifts: a flattening and upturn with increasing mass. The concentration c(M,z) as a function of mass and the redshift and for different cosmological parameters shows a remarkably complex pattern. However, when expressed in terms of the linear rms fluctuation of the density field sigma(M,z), the halo concentration c(sigma) shows a nearly-universal simple U-shaped behaviour with a minimum at a well defined scale at sigma=0.71. Yet, some small dependences with redshift and cosmology still remain. At the high-mass end (sigma < 1) the median halo kinematic profiles show large signatures of infall and highly radial orbits. This c-sigma(M,z) relation can be accurately parametrized and provides an analytical model for the dependence of concentration on halo mass. When applied to galaxy clusters, our estimates of concentrations are substantially larger -- by a factor up to 1.5 -- than previous results from smaller simulations, and are in much better agreement with results of observations. (abridged)

AIMS. Several authors have claimed to detect a significant cross-correlation between microwave WMAP anisotropies and the SDSS galaxy distribution. We repeat these analyses to determine the different cross-correlation uncertainties caused by re-sampling errors and field-to-field fluctuations. The first type of error concerns overlapping sky regions, while the second type concerns non-overlapping sky regions. METHODS. To measure the re-sampling errors, we use bootstrap and jack-knife techniques. For the field-to-field fluctuations, we use three methods: 1) evaluation of the dispersion in the cross-correlation when correlating separated regions of WMAP with the original region of SDSS; 2) use of mock Monte Carlo WMAP maps; 3) a new method (developed in this article), which measures the error as a function of the integral of the product of the self-correlations for each map. RESULTS. The average cross-correlation for b>30 deg. is significantly stronger than the re-sampling errors--both the jack-knife and bootstrap techniques provide similar results--but it is of the order of the field-to-field fluctuations. This is confirmed by the cross-correlation between anisotropies and galaxies in more than the half of the sample being null within re-sampling errors. CONCLUSIONS. Re-sampling methods underestimate the errors. Field-to-field fluctuations dominate the detected signals. The ratio of signal to re-sampling errors is larger than unity in a way that strongly depends on the selected sky region. We therefore conclude that there is no evidence yet of a significant detection of the integrated Sachs-Wolfe (ISW) effect. Hence, the value of Omega_\Lambda ~ 0.8 obtained by the authors who assumed they were observing the ISW effect would appear to have originated from noise analysis

We provide a set of general tools for studying the alignments of dark matter halos and galaxies with respect to the large scale structure. The statistics of the positioning of these objects is represented by a Probability Distribution Function (PDF) of their Euler angles. The PDF corresponding to halos located in the shells of the cosmic voids is inferred from previous results. This PDF is used to show how to recover the outcomes found for the alignments of the axes of these halos in simulations. We also explore the orientation of the angular momentum of the halos, both with respect to the halo axes and with respect to the large scale structure. We present an expression which describes well numerical results for the alignment of the angular momentum of the halo with respect to the halo axes for randomly chosen halos. We also propose a model that relates the orientation of the angular momentum with the halos axes accounting for the orientation of the halo axes with the large scale structure. This model is shown to recover accurately the observed PDF of the halo angular momentum with respect to the void radial direction. In addition, we give an expression for determining the degradation of the angular momentum intrinsic alignment when observational errors are accounted. This expression is also used to determine the departure of the observed value of the alignment from the initial expectation (as provided by the tidal torque theory) due to the rotation of the angular momentum of the halo with respect to the initial torque. For voids, we find that the strength of the alignment is reduced to half the original value. We discuss how to adapt the void results to other cosmic large scale structures (i.e. filaments, walls, etc).

We present a general analytical formalism to calculate accurately several statistics related to underdense regions in the Universe. The statistics are computed for dark matter halo and galaxy distributions both in real space and redshift space at any redshift. Using this formalism, we found that void statistics for galaxy distributions can be obtained, to a very good approximation, assuming galaxies to have the same clustering properties as halos above a certain mass. We deducted a relationship between this mass and that of halos with the same accumulated number density as the galaxies. We also found that the dependence of void statistics on redshift is small. For instance, the number of voids larger than 13 Mpc/h (defined to not contain galaxies brighter than M_r=-20.4 +5logh change less than 20% between z=1 and z=0. However, the dependence of void statistics on sigma_8 and Omega_m h is considerably larger, making them appropriate to develop tests to measure these parameters. We have shown how to efficiently construct several of these tests and discussed in detail the treatment of several observational effects. The formalism presented here along with the observed statistics extracted from current and future large galaxy redshift surveys will provide an independent measurement of the relevant cosmological parameters. Combining these measurements with those found using other methods will contribute to reduce their uncertainties.

AIMS. The amplitude and scaleheight of the Galactic gas disc density are not axisymmetric against expectations in a self-gravity axisymmetric disc. However, this lopsidedness can be explained in terms of intergalactic accretion flows, which produce non-axisymmetric pressure on the disc. This mechanism could be also responsible for the formation of a warp. METHODS. We analytically derive the relationship between the disc density and the self-gravity and external pressure. RESULTS. The same scenario of accretion as we proposed years ago to explain the formation of the warp explains the azimuthal dependence of the density and its scaleheight, with minimum/maximum in the positions of maximum amplitude of the warp (phi=95 deg. and 275 deg.), as expected from its pressure distribution.

CONTEXT. Hydrodynamical cosmological simulations predict flows of the intergalactic medium along the radial vector of the voids, approximately in the direction of the infall of matter at the early stages of the galaxy formation. AIMS. These flows might be detected by analysing the dependence of the warp amplitude on the inclination of the galaxies at the shells of the voids with respect to the radial vector of the voids. This analysis will be the topic of this paper. METHODS. We develop a statistical method of analysing the correlation of the amplitude of the warp and the inclination of the galaxy at the void surface. This is applied to a sample of 97 edge-on galaxies from the Sloan Digital Sky Survey. Our results are compared with the theoretical expectations, which are also derived in this paper. RESULTS. Our results allow us to reject the null hypothesis (i.e., the non-correlation of the warp amplitude and the inclination of the galaxy with respect to the void surface) at 94.4% C. L., which is not conclusive. The absence of the radial flows cannot be excluded at present, although we can put a constraint on the maximum average density of baryonic matter of the radial flows of <rho_b> <~ 4 Omega_b rho_crit.

[ABRIDGED] The unconditional mass function (UMF) of dark matter haloes has been determined accurately in the literature, showing excellent agreement with high resolution numerical simulations. However, this is not the case for the conditional mass function (CMF). We propose a simple analytical procedure to derive the CMF by rescaling the UMF to the constrained environment using the appropriate mean and variance of the density field at the constrained point. This method introduces two major modifications with respect to the standard re-scaling procedure. First of all, rather than using in the scaling procedure the properties of the environment averaged over all the conditioning region, we implement the re-scaling locally. We show that for high masses this modification may lead to substantially different results. Secondly, we modify the (local) standard re-scaling procedure in such a manner as to force normalisation, in the sense that when one integrates the CMF over all possible values of the constraint multiplied by their corresponding probability distribution, the UMF is recovered. In practise, we do this by replacing in the standard procedure the value delta_c (the linear density contrast for collapse) by certain adjustable effective parameter delta_eff. In order to test the method, we compare our prescription with the results obtained from numerical simulations in voids (Gottlober et al. 2003), finding a very good agreement. Based on these results, we finally present a very accurate analytical fit to the (accumulated) conditional mass function obtained with our procedure, which may be useful for any theoretical treatment of the large scale structure.

Using the high resolution cosmological N-body simulation MareNostrum Universe we study the orientation of shape and angular momentum of galaxy-size dark matter haloes around large voids. We find that haloes located on the shells of the largest cosmic voids have angular momenta that tend to be preferentially perpendicular to the direction that joins the centre of the halo and the centre of the void. This alignment has been found in spiral galaxies around voids using galaxy redshift surveys. We measure for the first time the strength of this alignment, showing how it falls off with increasing distance to the centre of the void. We also confirm the correlation between the intensity of this alignment and the halo mass. The analysis of the orientation of the halo main axes confirms the results of previous works. Moreover, we find a similar alignment for the baryonic matter inside dark matter haloes, which is much stronger in their inner parts.

In this work, we study the formation and evolution of dark matter halos by means of the spherical infall model with shell-crossing. We present a framework to tackle this effect properly based on the numerical follow-up, with time, of that individual shell of matter that contains always the same fraction of mass with respect to the total mass. In this first step, we do not include angular momentum, velocity dispersion or triaxiality. Within this framework - named as the Spherical Shell Tracker (SST) - we investigate the dependence of the evolution of the halo with virial mass, with the adopted mass fraction of the shell, and for different cosmologies. We find that our results are very sensitive to a variation of the halo virial mass or the mass fraction of the shell that we consider. However, we obtain a negligible dependence on cosmology. Furthermore, we show that the effect of shell-crossing plays a crucial role in the way that the halo reaches the stabilization in radius and the virial equilibrium. We find that the values currently adopted in the literature for the actual density contrast at the moment of virialization, delta_vir, may not be accurate enough. In this context, we stress the problems related to the definition of a virial mass and a virial radius for the halo. The question of whether the results found here may be obtained by tracking the shells with an analytic approximation remains to be explored.

The accretion of the intergalactic medium onto the gaseous disc is used to explain the generation of galactic warps. A cup-shaped distortion is expected, due to the transmission of the linear momentum; but, this effect is small for most incident inflow angles and the predominant effect turns out to be the transmission of angular momentum, i.e. a torque giving an integral-sign shaped warp. The torque produced by a flow of velocity ~100 km/s and baryon density ~10^{-25} kg/m^3, which is within the possible values for the intergalactic medium, is enough to generate the observed warps and this mechanism offers quite a plausible explanation. The inferred rate of infall of matter, ~1 M_sun/yr, to the Galactic disc that this theory predicts agrees with the quantitative predictions of chemical evolution resolving key issues, notably the G-dwarf problem. Sanchez-Salcedo (2006) suggests that this mechanism is not plausible because it would produce a dependence of the scaleheight of the disc with the Galactocentric azimuth in the outer disc, but rather than being an objection this is another argument in favour of the mechanism because this dependence is actually observed in our Galaxy.

In the present work we describe the formalism necessary to derive the properties of dark matter halos beyond two virial radius using the spherical collapse model (without shell crossing), and provide the framework for the theoretical prediction presented in Prada et al. (2005). We show in detail how to obtain within this model the probability distribution for the spherically-averaged enclosed density at any radii P(delta,r). Using this probability distribution, we compute the most probable and mean density profiles, which turns out to differ considerably from each other. We also show how to obtain the typical profile, as well as the probability distribution and mean profile for the spherically averaged radial velocity. Two probability distributions are obtained: a first one is derived using a simple assumption, that is, if Q is the virial radius in Lagrangian coordinates, then the enclosed linear contrast delta_l(q,Q) must satisfy the condition that delta_l(q=Q) = delta_vir, where delta_vir is the linear density contrast within the virial radius Rvir at the moment of virialization. Then we introduce an additional constraint to obtain a more accurate P(delta,r) which reproduces to a higher degree of precision the distribution of the spherically averaged enclosed density found in the simulations. This new constraint is delta_l(q,Q) < delta_vir for all q > Q, which means that there are no radii larger than Rvir where the density contrast is larger than that used to define the virial radius. Finally, we compare in detail our theoretical predictions for the probability distributions with the results found in the simulations.

We present a statistical analysis of voids in the 2dF galaxy redshift survey (2dFGRS). In order to detect the voids, we have developed two robust algorithms. We define voids as non-overlapping maximal spheres empty of halos or galaxies with mass or luminosity above a given one. We search for voids in cosmological $N$-Body simulations to test the performance of our void finders. We obtain and analyze the void statistics for several volume-limited samples for the North Galactic Strip (NGP) and the South Galactic Strip (SGP) constructed from the 2dFGRS full data release. We find that the results obtained from the NGP and the SGP are statistically compatible. From the results of several statistical tests we conclude that voids are essentially uncorrelated, with at most a mild anticorrelation and that there is a dependence of the void number density on redshift at least at the 99.5% confidence level. We develop a technique to correct the distortion caused by the fact that we use the redshift as the radial coordinate. We calibrate this technique with mock catalogues and find that the correction might be of some relevance to carry out accurate inferences from void statistics. We study the statistics of the galaxies inside nine nearby voids. We find that galaxies in voids are not randomly distributed: they form structures like filaments. We also obtain the galaxy number density profile in voids. This profile follow a similar but steeper trend to that follow by halos in voids.

We investigate the orientation of the axes and angular momentum of dark matter halos with respect to their neighboring voids using high resolution N-body cosmological simulations. We find that the minor axis of halos tends to be aligned along the line joining the halo with the center of the void, and that the major axis is perpendicular to this line. However, we find that the angular momentum of halos does not have any particular orientation. These results may provide information about the mechanisms whereby the large-scale structure of the Universe affects galaxy formation and cast light upon the issue of the orientation of galaxy disks with respect to their host halos.

It has recently been shown that the correct halo mass function high mass behaviour may be obtained by an appropriate treatment of the all-mass-at-center problem. Here we show that the treatment of this problem leads to an accumulated mass fraction, F, which is not a universal function of the linear variance of the density contrast, but depends also on the shape of the power spectrum. For relevant power spectra F is a universal function of two quantities: the linear variance and a quantity related to the local shape (at the given mass scale) of the power spectrum. This implies, for a given shape of the power spectrum, a redshift dependence of the mass function, which is quite important for a cold dark matter power spectrum. Our results are in good agreement with numerical simulations.

We claim that the discrepancy found between theoretical predictions for the cosmic mass function and those found in numerical simulations is due to the fact that in deriving the former all mass elements are assumed to be at the center of the object they belong to (the all-mass-at-center problem). By an appropriate treatment of this problem, using the spherical collapse model (which is not a bad approximation in the high mass limit), we obtain both the high mass behaviour found in simulations and the true assymptotic behaviour of the mass function (for arbitrarily high masses). Therefore, we conclude that by combining ellipsoidal dynamics with a suitable treatment of the all-mass-at-center problem (that we will show in a follow up work) a theoretical prediction for the cosmic mass function in full agreement with simulations may be obtained.

We present a general analytical procedure for computing the number density of voids with radius above a given value within the context of gravitational formation of the large scale structure of the universe out of Gaussian initial conditions. To this end we develop an accurate (under generally satisfied conditions) extension of unconditional mass function to constrained environments, which allowes us both to obtain the number density of collapsed objects of certain mass at any distance from the center of the void, and to derive the number density of voids defined by those collapsed objects. We have made detailed calculations for the spherically averaged mass density and halo number density profiles for individual voids. We also present a formal expression for the number density of voids defined by galaxies of a given type and luminosity. This expression contains the probability for a collapsed object of certain mass to host a galaxy of that type and luminosity as a function of the environmental density. We propose a procedure to infer this function, which may provide useful clues as to the galaxy formation process, from the observed void densities.

We present a comparison of the properties of galaxies in the most underdense regions of the Universe, where the galaxy number density is less than 10% of the mean density, with galaxies from more typical regions. We have compiled a sample of galaxies in 46 large nearby voids that were identified using the Sloan Digital Sky Survey DR4, which provides the largest coverage of the sky. We study the u-r color distribution, morphology, specific star formation rate, and radial number density profiles for a total of 495 galaxies fainter than M_r=-20.4 +5logh located inside the voids and compare these properties with a control sample of field galaxies. We show that there is an excess of blue galaxies inside the voids. However, inspecting the properties of blue and red galaxies separately, we find that galaxy properties such as color distribution, bulge-to-total ratios, and concentrations are remarkably similar between the void and overall sample. The void galaxies also show the same specific star formation rate at fixed color as the control galaxies. We compare our results with the predictions of cosmological simulations of galaxy formation using the Millennium Run semi-analytic galaxy catalog. We show that the properties of the simulated galaxies in large voids are in reasonably good agreement with those found in similar environments in the real Universe. To summarize, in spite of the fact that galaxies in voids live in the least dense large-scale environment, this environment makes very little impact on properties of galaxies.

We study the density profiles of collapsed galaxy-size dark matter halos with masses 1e11-5e12 Msun focusing mostly on the halo outer regions from the formal virial radius Rvir up to 5-7Rvir. We find that isolated halos in this mass range extend well beyond Rvir exhibiting all properties of virialized objects up to 2-3Rvir: relatively smooth density profiles and no systematic infall velocities. The dark matter halos in this mass range do not grow as one naively may expect through a steady accretion of satellites, i.e., on average there is no mass infall. This is strikingly different from more massive halos, which have large infall velocities outside of the virial radius. We provide accurate fit for the density profile of these galaxy-size halos. For a wide range (0.01-2)Rvir of radii the halo density profiles are fit with the approximation rho=rho_s exp(-2n[x^{1/n}-1])+rho_m, where x=r/r_s, rho_m is the mean matter density of the Universe, and the index n is in the range n=6-7.5. These profiles do not show a sudden change of behavior beyond the virial radius. For larger radii we combine the statistics of the initial fluctuations with the spherical collapse model to obtain predictions for the mean and most probable density profiles for halos of several masses. The model give excellent results beyond 2-3 formal virial radii.

We measure the variance in the distribution of off-plane (|b|>20 deg.) galaxies with m_K<13.5 from the 2MASS K-band survey in circles of diameter between 0.344 deg. and 57.2 deg. The use of a near-infrared survey makes negligible the contribution of Galactic extinction to these fluctuations. We calculate these variances within the standard Lambda-CDM model assuming that the sources are distributed like halos of the corresponding mass, and it reproduces qualitatively the galaxy counts variance. Therefore, we test that the counts can be basically explained in terms only of the large scale structure. A second result of this paper is a new method to determine the two point correlation function obtained by forcing agreement between model and data. This method does not need the knowledge of the two-point angular correlation function, allows an estimation of the errors (which are low with this method), and can be used even with incomplete surveys. Using this method we get xi(z=0, r<10 h^{-1}Mpc)=(29.8+/-0.3) (r/h^{-1}Mpc)^{-1.79+/-0.02}, which is the first measure of the amplitude of xi in the local Universe for the K-band. It is more or less in agreement with those obtained through red optical filters selected samples, but it is larger than the amplitude obtained for blue optical filters selected samples.

We present a new general procedure for determining a given set of quantities. To this end, we define certain statistic, that we call 'modified $\chi^2$' ($\chi^2_M$), because of its similarity with the standard $\chi^2$. The terms of this $\chi^2_M$ are made up of the fluctuations of an unbiased estimator of some statistical quantities, and certain weights. Only the diagonal terms of the covariance matrix explicitly appear in our statistic, while the full covariance matrix (and not its inverse) is implicitly included in the calculation of the weights. Choosing these weights we may obtain, through minimising the $\chi^2_M$, the estimator that provides the minimum RMS, either for those quantities or for the parameters on which these quantities depend. In this paper, we describe our method in the context of Cosmic Microwave Background experiments, in order to obtain either the statistical properties of the maps, or the cosmological parameters. The test here is constructed out of some estimator of the two-point correlation function at different angles. For the problem of one parameter estimation, we show that our method has the same power as the maximum likelihood method. We have also applied this method to Monte Carlo simulations of the COBE-DMR data, as well as to the actual 4-year data, obtaining consistent results with previous analyses. We also provide a very good analytical approximation to the distribution function of our statistic, which could also be useful in other contexts.

A new method is developed to calculate the amplitude of the galactic warps generated by a torque due to external forces. This takes into account that the warp is produced as a reorientation of the different rings which constitute the disc in order to compensate the differential precession generated by the external force, yielding a uniform asymptotic precession for all rings. Application of this method to gravitational tidal forces in the Milky Way due to the Magellanic Clouds leads to a very low amplitude of the warp. If the force were due to an extragalactic magnetic field, its intensity would have to be very high, to generate the observed warps. An alternative hypothesis is explored: the accretion of the intergalactic medium over the disk. A cup-shaped distortion is expected, due to the transmission of the linear momentum; but, this effect is small and the predominant effect turns out to be the transmission of angular momentum, i.e. a torque giving an integral-sign shape warp. The torque produced by a flow of velocity ~100 km/s and baryon density \~10^{-25} kg/m^3 is enough to generate the observed warps and this mechanism offers quite a plausible explanation. First, because this order of accretion rate is inferred from other processes observed in the Galaxy, notably its chemical evolution. The inferred rate of infall of matter, ~1 solar-mass/yr, to the Galactic disc that this theory predicts agrees with the quantitative predictions of this chemical evolution resolving key issues, notably the G-dwarf problem. Second, because the required density of the intergalactic medium is within the range of values compatible with observation. By this mechanism, we can explain the warp phenomenon in terms of intergalactic accretion flows onto the disk of the galaxy.

We present a general procedure for obtaining the present density fluctuation probability distribution given the statistics of the initial conditions. The main difficulties faced with regard to this problem are those related to the non-linear evolution of the density fluctuations and those posed by the fact that the fields we are interested in are the result of filtering an underlying field with structure down to scales much smaller than that of filtering. The solution to the latter problem is discussed here in detail and the solution to the former is taken from a previous work. We have checked the procedure for values of the rms density fluctuation as large as 3/2 and several power spectra and found that it leads to results in excellent agreement with those obtained in numerical simulations. We also recover all available exact results from perturbation theory.

We have developed a generalization of the Zeldovich approximation (ZA) that is exact in a wide variety of situations, including plannar, spherical and cilyndrical symmetries. We have shown that this generalization, that we call complete Zeldovich approximation (CZA), is exact to second order at an arbitrary point within any field. For gaussian fields, the third order error have been obtained and shown to be very small. For statistical purposes, the CZA leads to results exact to the third order.