We present a full description of the N-probability density function of the galaxy number density fluctuations. This N-pdf is given in terms, on the one hand, of the cold dark matter correlations and, on the other hand, of the galaxy bias parameter. The method relies on the assumption commonly adopted that the dark matter density fluctuations follow a local non-linear transformation of the initial energy density perturbations. The N-pdf of the galaxy number density fluctuations allows for an optimal estimation of the bias parameter (e.g., via maximum-likelihood estimation, or Bayesian inference if there exists any a priori information on the bias parameter), and of those parameters defining the dark matter correlations, in particular its amplitude ($\sigma_8$). It also provides the proper framework to perform model selection between two competitive hypotheses. The parameters estimation capabilities of the N-pdf are proved by SDSS-like simulations (both ideal log-normal simulations and mocks obtained from Las Damas simulations), showing that our estimator is unbiased. We apply our formalism to the 7th release of the SDSS main sample (for a volume-limited subset with absolute magnitudes $M_r \leq -20$). We obtain $\hat{b} = 1.193 \pm 0.074$ and $\hat{\sigma_8} = 0.862 \pm 0.080$, for galaxy number density fluctuations in cells of a size of $30h^{-1}$Mpc. Different model selection criteria show that galaxy biasing is clearly favoured.

The detection of point sources in Cosmic Microwave Background maps is usually based on a single-frequency approach, whereby maps at each frequency are filtered separately and the spectral information on the sources is derived combining the results at the different frequencies. On the contrary, in the case of multi-frequency detection methods, source detection and spectral information are tightly interconnected in order to increase the source detection efficiency. In this work we apply the \emph{matched multifilter} method to the detection of point sources in the WMAP 7yr data at 61 and 94 GHz. This linear filtering technique takes into account the spatial and the cross-power spectrum information at the same time using the spectral behaviour of the sources without making any a priori assumption about it. We follow a two-step approach. First, we do a blind detection of the sources over the whole sky. Second, we do a refined local analysis at their positions to improve the signal-to-noise ratio of the detections. At 94 GHz we detect 129 5{\sigma} objects at |b|>5{\deg} (excluding the Large Magellanic Cloud region); 119 of them are reliable extragalactic sources and 104 of these 119 lie outside the WMAP Point Source Catalog mask. Nine of the total 129 detections are known Galactic sources or lie in regions of intense Galactic emission and one additional (weak) high-Galactic latitude source has no counterpart in low-frequency radio catalogues. Our results constitute a substantial improvement over the NEWPS-3year catalogue.

We propose to model the image differentials of astrophysical source maps by Student's t-distribution and to use them in the Bayesian source separation method as priors. We introduce an efficient Markov Chain Monte Carlo (MCMC) sampling scheme to unmix the astrophysical sources and describe the derivation details. In this scheme, we use the Langevin stochastic equation for transitions, which enables parallel drawing of random samples from the posterior, and reduces the computation time significantly (by two orders of magnitude). In addition, Student's t-distribution parameters are updated throughout the iterations. The results on astrophysical source separation are assessed with two performance criteria defined in the pixel and the frequency domains.

The detection and flux estimation of point sources in cosmic microwave background (CMB) maps is a very important task in order to clean the maps and also to obtain relevant astrophysical information. In this paper we propose a maximum a posteriori (MAP) approach detection method in a Bayesian scheme which incorporates prior information about the source flux distribution, the locations and the number of sources. We apply this method to CMB simulations with the characteristics of the Planck satellite channels at 30, 44, 70 and 100 GHz. With a similar level of spurious sources, our method yields more complete catalogues than the matched filter with a 5 sigma threshold. Besides, the new technique allows us to fix the number of detected sources in a non-arbitrary way.

We propose a Bayesian approach to joint source separation and restoration for astrophysical diffuse sources. We constitute a prior statistical model for the source images by using their gradient maps. We assume a t-distribution for the gradient maps in different directions, because it is able to fit both smooth and sparse data. A Monte Carlo technique, called Langevin sampler, is used to estimate the source images and all the model parameters are estimated by using deterministic techniques.

We consider a set of M images, whose pixel intensities at a common point can be treated as the components of a M-dimensional vector. We are interested in the estimation of the modulus of such a vector associated to a compact source. For instance, the detection/estimation of the polarized signal of compact sources immersed in a noisy background is relevant in some fields like Astrophysics. We develop two different techniques, one based on the Maximum Likelihood Estimator (MLE) applied to the modulus distribution, the modulus filter (ModF) and other based on prefiltering the components before fusion, the filtered fusion (FF), to deal with this problem. We present both methods in the general case of M images and apply them to the particular case of three images (linear plus circular polarization). Numerical simulations have been performed to test these filters considering polarized compact sources immersed in stationary noise. The FF performs better than the ModF in terms of errors in the estimated amplitude and position of the source, especially in the low signal-to-noise case. We also compare both methods with the direct application of a matched filter (MF) on the polarization data. This last technique is clearly outperformed by the new methods.

[Abridged] In this paper we explore a local non-linear perturbative model up to third order as a general characterization of the CMB anisotropies. We focus our analysis in large scale anisotropies. At these angular scales, the non-Gaussian description proposed in this work defaults (under certain conditions) to an approximated local form of the weak non-linear coupling inflationary model. In particular, quadratic and cubic terms are governed by the non-linear coupling parameters f_nl and g_nl, respectively. The extension proposed in this paper allows us to directly constrain these non-linear parameters. Applying the proposed methodology to WMAP 5-yr data, we obtain -5.6 x 10^5 < g_nl < 6.4 x 10^5, at 95% CL. This result is in agreement with previous findings obtained for equivalent non-Gaussian models and with different non-Gaussian estimators. A model selection test is performed, indicating that a Gaussian model is preferred to the non-Gaussian scenario. When comparing different non-Gaussian models, we observe that a pure f_nl model is the most favoured case, and that a pure g_nl model is more likely than a general non-Gaussian scenario. Finally, we have analyzed the WMAP data in two independent hemispheres, in particular the ones defined by the dipolar pattern found by Hoftuft et al. 2009. We show that, whereas g_nl is still compatible with zero for both hemispheres, it is not the case for f_nl (with a p-value 0.04). However, if anisotropy of the data is assumed, the distance between the likelihood distributions for each hemisphere is larger than expected from Gaussian and anisotropic simulations, also for g_nl (with a p-value of 0.001 in the case of this parameter). This result is an extra evidence for the CMB asymmetries previously reported in WMAP data.

The detection of polarized sources in the WMAP 5-year data is a very difficult task. The maps are dominated by instrumental noise and only a handful of sources show up as clear peaks in the Q and U maps. Optimal linear filters applied at the position of known bright sources detect with a high level of significance a polarized flux P from many more sources, but estimates of P are liable to biases. Using a new technique, named the "filtered fusion technique", we have detected in polarization, with a significance level greater than 99.99% in at least one WMAP channel, 22 objects, 5 of which, however, do not have a plausible low radio frequency counterpart and are therefore doubtful. Estimated polarized fluxes P < 400 mJy at 23 GHz were found to be severely affected by the Eddington bias. The corresponding polarized flux limit for Planck/LFI at 30 GHz, obtained via realistic simulations, is 300 mJy. We have also obtained statistical estimates of, or upper limits to the mean polarization degrees of bright WMAP sources at 23, 33, 41, and 61 GHz, finding that they are of a few percent.

In this work we deal with the problem of simultaneous multifrequency detection of extragalactic point sources in maps of the Cosmic Microwave Background. We apply a linear filtering technique that uses spatial information and the cross-power spectrum. To make this, we simulate realistic and non-realistic flat patches of the sky at two frequencies of Planck: 44 and 100 GHz. We filter to detect and estimate the point sources and compare this technique with the monofrequency matched filter in terms of completeness, reliability, flux and spectral index accuracy. The multifrequency method outperforms the matched filter at the two frequencies and in all the studied cases in the work.

Given a set of images, whose pixel values can be considered as the components of a vector, it is interesting to estimate the modulus of such a vector in some localised areas corresponding to a compact signal. For instance, the detection/estimation of a polarized signal in compact sources immersed in a background is relevant in some fields like astrophysics. We develop two different techniques, one based on the Neyman-Pearson lemma, the Neyman-Pearson filter (NPF), and another based on prefiltering-before-fusion, the filtered fusion (FF), to deal with the problem of detection of the source and estimation of the polarization given two or three images corresponding to the different components of polarization (two for linear polarization, three including circular polarization). For the case of linear polarization, we have performed numerical simulations on two-dimensional patches to test these filters following two different approaches (a blind and a non-blind detection), considering extragalactic point sources immersed in cosmic microwave background (CMB) and non-stationary noise with the conditions of the 70 GHz \emph{Planck} channel. The FF outperforms the NPF, especially for low fluxes. We can detect with the FF extragalactic sources in a high noise zone with fluxes >= (0.42,0.36) Jy for (blind/non-blind) detection and in a low noise zone with fluxes >= (0.22,0.18) Jy for (blind/non-blind) detection with low errors in the estimated flux and position.

We present a new method based on the N-point probability distribution (pdf) to study non-Gaussianity in cosmic microwave background (CMB) maps. Likelihood and Bayesian estimation are applied to a local non-linear perturbed model up to third order, characterized by a linear term which is described by a Gaussian N-pdf, and a second and third order terms which are proportional to the square and the cube of the linear one. We also explore a set of model selection techniques (the Akaike and the Bayesian Information Criteria, the minimum description length, the Bayesian Evidence and the Generalized Likelihood Ratio Test) and their application to decide whether a given data set is better described by the proposed local non-Gaussian model, rather than by the standard Gaussian temperature distribution. As an application, we consider the analysis of the WMAP 5-year data at a resolution of around 2 degrees. At this angular scale (the Sachs-Wolfe regime), the non-Gaussian description proposed in this work defaults (under certain conditions) to an approximative local form of the weak non-linear coupling inflationary model (e.g. Komatsu & Spergel 2001) previously addressed in the literature. For this particular case, we obtain an estimation for the non-linear coupling parameter of -94 < F_nl < 154 at 95% CL. Equally, model selection criteria also indicate that the Gaussian hypothesis is favored against the particular local non-Gaussian model proposed in this work. This result is in agreement with previous findings obtained for equivalent non-Gaussian models and with different non-Gaussian estimators. However, our estimator based on the N-pdf is more efficient than previous estimators and, therefore, provides tighter constraints on the coupling parameter at degree angular resolution.

The Planck satellite will map the full sky at nine frequencies from 30 to 857 GHz. The CMB intensity and polarization that are its prime targets are contaminated by foreground emission. The goal of this paper is to compare proposed methods for separating CMB from foregrounds based on their different spectral and spatial characteristics, and to separate the foregrounds into components of different physical origin. A component separation challenge has been organized, based on a set of realistically complex simulations of sky emission. Several methods including those based on internal template subtraction, maximum entropy method, parametric method, spatial and harmonic cross correlation methods, and independent component analysis have been tested. Different methods proved to be effective in cleaning the CMB maps from foreground contamination, in reconstructing maps of diffuse Galactic emissions, and in detecting point sources and thermal Sunyaev-Zeldovich signals. The power spectrum of the residuals is, on the largest scales, four orders of magnitude lower than that of the input Galaxy power spectrum at the foreground minimum. The CMB power spectrum was accurately recovered up to the sixth acoustic peak. The point source detection limit reaches 100 mJy, and about 2300 clusters are detected via the thermal SZ effect on two thirds of the sky. We have found that no single method performs best for all scientific objectives. We foresee that the final component separation pipeline for Planck will involve a combination of methods and iterations between processing steps targeted at different objectives such as diffuse component separation, spectral estimation and compact source extraction.

We have analyzed the efficiency in source detection and flux density estimation of blind and non-blind detection techniques exploiting the MHW2 filter applied to the Wilkinson Microwave Anisotropy Probe (WMAP) 5-year maps. A comparison with the AT20G Bright Source Sample (Massardi et al. 2008), with a completeness limit of 0.5 Jy and accurate flux measurements at 20 GHz, close to the lowest frequency of WMAP maps, has allowed us to assess the completeness and the reliability of the samples detected with the two approaches, as well as the accuracy of flux and error estimates, and their variations across the sky. The uncertainties on flux estimates given by our procedure turned out to be about a factor of 2 lower than the rms differences with AT20G measurements, consistent with the smoothing of the fluctuation field yielded by map filtering. Flux estimates were found to be essentially unbiased except that, close to the detection limit, a substantial fraction of fluxes are found to be inflated by the contribution of underlying positive fluctuations. This is consistent with expectations for the Eddington bias associated to the true errors on flux density estimates. The blind and non-blind approaches are found to be complementary: each of them allows the detection of sources missed by the other. Combining results of the two methods on the WMAP 5-year maps we have expanded the non-blindly generated New Extragalactic WMAP Point Source (NEWPS) catalogue (Lopez-Caniego et al. 2007) that was based on WMAP 3-year maps. After having removed the probably spurious objects not identified with known radio sources, the new version of the NEWPS catalogue, NEWPS_5yr comprises 484 sources detected with a signal-to-noise ratio SNR>5.

In this work we address the problem of simultaneous multi-frequency detection of extragalactic point sources in maps of the Cosmic Microwave Background. We apply a new linear filtering technique, the so called `matched matrix filters', that incorporates full spatial information, including the cross-correlation among channels, without making any a priori assumption about the spectral behaviour of the sources. A substantial reduction of the background is achieved thanks to the optimal combination of filtered maps. We describe in detail the new technique and we apply it to the detection/estimation of radio sources in realistic all-sky Planck simulations at 30, 44, 70 and 100 GHz. Then we compare the results with the mono-frequential approach based on the standard matched filter, in terms of reliability, completeness and flux accuracy of the resulting point source catalogs. The new filters outperform the standard matched filters for all these indexes at 30, 44 and 70 GHz, whereas at 100 GHz both kind of filters have a similar performance. We find a noticeable increment of the number of true detections for a fixed reliability level. In particular, for a 95% reliability we practically double the number of detections at 30, 44 and 70 GHz.

In this paper we introduce a new linear filtering technique, the so-called matrix filters, that maximizes the signal-to-interference ratio of compact sources of unknown intensity embedded in a set of images by taking into account the cross-correlations between the different channels. By construction, the new filtering technique outperforms (or at least equals) the standard matched filter applied on individual images. An immediate application is the detection of extragalactic point sources in Cosmic Microwave Background images obtained at different wavelengths. We test the new technique in two simulated cases: a simple two-channel case with ideal correlated color noise and more realistic simulations of the sky as it will be observed by the LFI instrument of the upcoming ESA's Planck mission. In both cases we observe an improvement with respect to the standard matched filter in terms of signal-to-noise interference, number of detections and number of false alarms.

We present results on spectral index distributions, number counts, redshift distribution and other general statistical properties of extragalactic point sources in the NEWPS5 sample L\'opez-Caniego et al. (2007). The flux calibrations at all the WMAP channels have been reassessed both by comparison with ground based observations and through estimates of the effective beam areas. The two methods yield consistent statistical correction factors. A search of the NED has yielded optical identifications for 89% of sources in the complete sub-sample of 252 sources with S/N>5 and S>1.1 Jy at 23 GHz; 5 sources turned out to be Galactic and were removed. The NED also yielded redshifts for 92% of the extragalactic sources at |b|>10deg. Their distribution was compared with model predictions; the agreement is generally good but a possible discrepancy is noted. Using the 5 GHz fluxes from the GB6 or PMN surveys, we find that 76% of the 191 extragalactic sources with S_23GHz>1.3,Jy can be classified as flat-spectrum sources between 5 and 23 GHz. A spectral steepening is observed at higher frequencies: only 59% of our sources are still flat-spectrum sources between 23 and 61 GHz and the average spectral indexes steepen from <alpha_5^23>= 0.01\pm 0.03 to <alpha_41^61>= 0.37\pm 0.03. We think, however, that the difference may be due to a selection effect. The source number counts have a close to Euclidean slope and are in good agreement with the predictions of the cosmological evolution model by De Zotti et al. (2005). The observed spectral index distributions were exploited to get model-independent extrapolations of counts to higher frequencies. The risks of such operations are discussed and reasons of discrepancies with other recent estimates are clarified.

The cosmic microwave background (CMB) is a relic radiation of the Big Bang and as such it contains a wealth of cosmological information. Statistical analyses of the CMB, in conjunction with other cosmological observables, represent some of the most powerful techniques available to cosmologists for placing strong constraints on the cosmological parameters that describe the origin, content and evolution of the Universe. The last decade has witnessed the introduction of wavelet analyses in cosmology and, in particular, their application to the CMB. We review here spherical wavelet analyses of the CMB that test the standard cosmological concordance model. The assumption that the temperature anisotropies of the CMB are a realisation of a statistically isotropic Gaussian random field on the sphere is questioned. Deviations from both statistical isotropy and Gaussianity are detected in the reviewed works, suggesting more exotic cosmological models may be required to explain our Universe. We also review spherical wavelet analyses that independently provide evidence for dark energy, an exotic component of our Universe of which we know very little currently. The effectiveness of accounting correctly for the geometry of the sphere in the wavelet analysis of full-sky CMB data is demonstrated by the highly significant detections of physical processes and effects that are made in these reviewed works.

We have used the MHW2 filter to obtain estimates of the flux densities at the WMAP frequencies of a complete sample of 2491 sources, mostly brighter than 500 mJy at 5 GHz, distributed over the whole sky excluding a strip around the Galactic equator (b < 5 degrees). After having detected 933 sources above the 3 sigma level in the MHW2 filtered maps - our New Extragalactic WMAP Point Source (NEWPS_3sigma) Catalogue - we are left with 381 sources above 5 sigma in at least one WMAP channel, 369 of which constitute our NEWPS_5sigma catalogue. It is remarkable to note that 98 (i.e. 26%) sources detected above 5 sigma are `new', they are not present in the WMAP catalogue. Source fluxes have been corrected for the Eddington bias. Our flux density estimates before such correction are generally in good agreement with the WMAP ones at 23 GHz. At higher frequencies WMAP fluxes tend to be slightly higher than ours, probably because WMAP estimates neglect the deviations of the point spread function from a Gaussian shape. On the whole, above the estimated completeness limit of 1.1 Jy at 23 GHz we detected 43 sources missed by the blind method adopted by the WMAP team. On the other hand, our low-frequency selection threshold left out 25 WMAP sources, only 12 of which, however, are 5 sigma detections and only 3 have fluxes S at 23 GHz > 1.1 Jy. Thus, our approach proved to be competitive with, and complementary to the WMAP one.

A new method is presented for the construction of a natural continuous wavelet transform on the sphere. It incorporates the analysis and synthesis with the same wavelet and the definition of translations and dilations on the sphere through the spherical harmonic coefficients. We construct a couple of wavelets as an extension of the flat Mexican Hat Wavelet to the sphere and we apply them to the detection of sources on the sphere. We remark that no projections are used with this methodology.

Cosmic microwave background (CMB) temperature anisotropies follow a Gaussian statistical distribution in the standard inflationary model, but there are non-Gaussian contributions due to astrophysical foregrounds. The detection of the non-Gaussian signal due to extragalactic point sources and its distinction from the possible intrinsic non-Gaussianity is an issue of great importance in CMB analysis. The Mexican Hat Wavelet Family (MHWF), which has been proved very useful for the detection of extragalactic point sources, is applied here to the study of non-Gaussianity due to point sources in CMB maps. We carry out simulations of CMB maps with the characteristics of the forthcoming Planck mission at 70 and 100 GHz and filter them with the MHWF. By comparing the skewness and kurtosis of simulated maps with and without point sources, we are able to detect clearly the non-Gaussian signal due to point sources for flux limits as low as 0.4 Jy (70 GHz) and 0.3 Jy (100 GHz). The second and third members of the MHWF perform better in this respect than the Mexican Hat Wavelet (MHW1) and much better than the Daubechies 4 wavelet. We have also estimated the third order, $K_3$, and fourth order, $K_4$, cumulants produced by point sources at these Planck channels by means of a fit with the MHWF. The average relative errors with respect to the real values are below 12% for fluxes down to 0.6 Jy (70 GHz) and 0.4 Jy (100 GHz). The values of these cumulants allow us to distinguish between different source counts models.

We study the power of several scalar quantities constructed on the sphere (presented in Monteserin et al. 2005) to detect non-Gaussianity on the temperature distribution of the cosmic microwave background (CMB). The test has been performed using non-Gaussian CMB simulations with injected skewness or kurtosis generated through the Edgeworth expansion. We have also taken into account in the analysis the effect of anisotropic noise and the presence of a Galactic mask. We find that the best scalars to detect an excess of skewness in the simulations are the derivative of the gradient, the fractional isotropy, the Laplacian and the shape index. For the kurtosis case, the fractional anisotropy, the Laplacian and the determinant are the quantities that perform better.

We study the detection of extragalactic point sources in two-dimensional flat simulations for all the frequencies of the forthcoming ESA's Planck mission. In this work we have used the most recent available templates of the microwave sky: as for the diffuse Galactic components and the Sunyaev-Zel'dovich clusters we have used the ``Plank Reference Sky Model''; as for the extragalactic point sources, our simulations - which comprise all the source populations relevant in this frequency interval - are based on up-to-date cosmological evolution models for sources. To consistently compare the capabilities of different filters for the compilation of the - hopefully - most complete blind catalogue of point sources, we have obtained three catalogues by filtering the simulated sky maps with: the Matched Filter (MF), the Mexican Hat Wavelet (MHW1) and the Mexican Hat Wavelet 2 (MHW2), the first two members of the MHW Family. For the nine Planck frequencies we show the number of real and spurious detections and the percentage of spurious detections at different flux detection limits as well as the completeness level of the catalogues and the average errors in the estimation of the flux density of detected sources. Allowing a 5 % of spurious detections, we obtain the following number of detections by filtering with the MHW2 an area equivalent to half of the sky: 580 (30 GHz), 342 (44 GHz), 341 (70 GHz), 730 (100 GHz), 1130 (143 GHz), 1233 (217 GHz), 990 (353 GHz), 1025 (545 GHz) and 3183 (857 GHz). Our current results indicate that the MF and the MHW2 yield similar results, whereas the MHW1 performs worse in some cases and especially at very low fluxes. This is a relevant result, because we are able to obtain comparable results with the well known Matched Filter and with this specific wavelet, the MHW2, which is much easier to implement and use.

We propose a new detection technique in the plane based on an isotropic wavelet family. This family is naturally constructed as an extension of the Gaussian-Mexican Hat Wavelet pair and for that reason we call it the Mexican Hat Wavelet Family (MHWF). We show the performance of these wavelets when dealing with the detection of extragalactic point sources in cosmic microwave background (CMB) maps: a very important issue within the most general problem of the component separation of the microwave sky. Specifically, flat two-dimensional simulations of the microwave sky comprising all astrophysical components plus instrumental noise have been analyzed for the channels at 30, 44 and 70 GHz of the forthcoming ESA's Planck mission Low Frequency Instrument (LFI). We adopt up-to-date cosmological evolution models of extragalactic sources able to fit well the new data on high-frequency radio surveys and we discuss our current results on point source detection by comparing them with those obtained using the Mexican Hat Wavelet (MHW) technique, which has been already proven a suitable tool for detecting point sources. By assuming a 5% reliability level, the first new members of the MHWF, at their ``optimal scale'', provide three point source catalogues on half of the sky (at galactic latitude |b|> 30) at 30, 44 and 70 GHz of 639, 387 and 340 extragalactic sources, respectively. The corresponding flux detection limits are 0.38, 0.45 and 0.47 Jy . By using the same simulated sky patches and at the same frequencies as before, the MHW at its optimal scale provides 543, 322 and 311 sources with flux detection limits of 0.44, 0.51 and 0.50 Jy, respectively (5% reliability level). These results show a clear improvement when we use the new members of the MHWF and, in particular, the MHW2 with respect to the MHW.

(Abridged) We have used the Rayner & Best (1989) smooth tests of goodness-of-fit to study the Gaussianity of the Very Small Array (VSA) data. Out of the 41 published VSA individual pointings dedicated to cosmological observations, 37 are found to be consistent with Gaussianity, whereas four pointings show deviations from Gaussianity. In two of them, these deviations can be explained as residual systematic effects of a few visibility points which, when corrected, have a negligible impact on the angular power spectrum. The non-Gaussianity found in the other two (adjacent) pointings seems to be associated to a local deviation of the power spectrum of these fields with respect to the common power spectrum of the complete data set, at angular scales of the third acoustic peak (l = 700-900). No evidence of residual systematics is found in this case, and unsubstracted point sources are not a plausible explanation either. If those visibilities are removed, a cosmological analysis based on this new VSA power spectrum alone shows no differences in the parameter constraints with respect to our published results, except for the physical baryon density, which decreases by 10 percent. Finally, the method has been also used to analyse the VSA observations in the Corona Borealis supercluster region (Genova-Santos et al. 2005), which show a strong decrement which cannot be explained as primordial CMB. Our method finds a clear deviation (99.82%) with respect to Gaussianity in the second-order moment of the distribution, and which can not be explained as systematic effects. A detailed study shows that the non-Gaussianity is produced in scales of l~500, and that this deviation is intrinsic to the data (in the sense that can not be explained in terms of a Gaussian field with a different power spectrum).

A method to compute several scalar quantities of Cosmic Microwave Background maps on the sphere is presented. We consider here four type of scalars: the Hessian matrix scalars, the distortion scalars, the gradient related scalars and the curvature scalars. Such quantities are obtained directly from the spherical harmonic coefficients (alm) of the map. We also study the probability density function of these quantities for the case of a homogeneous and isotropic Gaussian field, which are functions of the power spectrum of the initial field. From these scalars it is posible to construct a new set of scalars which are independent of the power spectrum of the field. We test our results using simulations and find a good agreement between the theoretical probability density functions and those obtained from simulations. Therefore, these quantities are proposed to investigate the presence of non-Gaussian features in CMB maps. Finally, we show how to compute the scalars in presence of anisotropic noise and realistic masks.

The Rayleigh-Taylor instability plays an important role in the dynamics of several astronomical objects, in particular, in supernovae (SN) evolution. In this paper we develop an analytical approach to study the stability analysis of spherical expansion of the SN ejecta by using a special transformation in the co-moving coordinate frame. We first study a non-stationary spherical expansion of a gas shell under the pressure of a central source. Then we analyze its stability with respect to a no radial, non spherically symmetric perturbation of the of the shell. We consider the case where the polytropic constant of the SN shell is $\gamma=5/3$ and we examine the evolution of a arbitrary shell perturbation. The dispersion relation is derived. The growth rate of the perturbation is found and its temporal and spatial evolution is discussed. The stability domain depends on the ejecta shell thickness, its acceleration, and the perturbation wavelength.

In this note we stress the necessity of a careful check of the arguments used by Vio & Andreani (2005) (VA hereinafter) to criticise the superior performance of the biparametric scale adaptive filter (BSAF) with respect to the classic matched filter (MF) in the detection of sources on a random Gaussian background. In particular, we point out that a defective reading and understanding of previous works in the literature (Rice 1954; Barreiro et al. 2003; Lopez-Caniego et al. 2005) leads the authors of VA to the derivation of an incorrect formula and to some misleading conclusions.

This paper considers the problem of compact source detection on a Gaussian background in 1D. Two aspects of this problem are considered: the design of the detector and the filtering of the data. Our detection scheme is based on local maxima and it takes into account not only the amplitude but also the curvature of the maxima. A Neyman-Pearson test is used to define the region of acceptance, that is given by a sufficient linear detector that is independent on the amplitude distribution of the sources. We study how detection can be enhanced by means of linear filters with a scaling parameter and compare some of them (the Mexican Hat wavelet, the matched and the scale-adaptive filters). We introduce a new filter, that depends on two free parameters (biparametric scale-adaptive filter). The value of these two parameters can be determined, given the a priori pdf of the amplitudes of the sources, such that the filter optimizes the performance of the detector in the sense that it gives the maximum number of real detections once fixed the number density of spurious sources. The combination of a detection scheme that includes information on the curvature and a flexible filter that incorporates two free parameters (one of them a scaling) improves significantly the number of detections in some interesting cases. In particular, for the case of weak sources embedded in white noise the improvement with respect to the standard matched filter is of the order of 40%. Finally, an estimation of the amplitude of the source is introduced and it is proven that such an estimator is unbiased and it has maximum efficiency. We perform numerical simulations to test these theoretical ideas and conclude that the results of the simulations agree with the analytical ones.

This article considers the detection of point sources in two dimensional astronomical images. The detection scheme we propose is based on peak statistics. We discuss the example of the detection of far galaxies in Cosmic Microwave Background experiments throughout the paper, although the method we present is totally general and can be used in many other fields of data analysis. We assume sources with a Gaussian profile --that is a fair approximation of the profile of a point source convolved with the detector beam in microwave experiments-- on a background modeled by a homogeneous and isotropic Gaussian random field characterized by a scale-free power spectrum. Point sources are enhanced with respect to the background by means of linear filters. After filtering, we identify local maxima and apply our detection scheme, a Neyman-Pearson detector that defines our region of acceptance based on the a priori pdf of the sources and the ratio of number densities. We study the different performances of some linear filters that have been used in this context in the literature: the Mexican Hat wavelet, the matched filter and the scale-adaptive filter. We consider as well an extension to two dimensions of the biparametric scale adaptive filter (BSAF). The BSAF depends on two parameters which are determined by maximizing the number density of real detections while fixing the number density of spurious detections. For our detection criterion the BSAF outperforms the other filters in the interesting case of white noise.

We consider filters for the detection and extraction of compact sources on a background. We make a one-dimensional treatment (though a generalization to two or more dimensions is possible) assuming that the sources have a Gaussian profile whereas the background is modeled by an homogeneous and isotropic Gaussian random field, characterized by a scale-free power spectrum. Local peak detection is used after filtering. Then, a Bayesian Generalized Neyman-Pearson test is used to define the region of acceptance that includes not only the amplification but also the curvature of the sources and the a priori probability distribution function of the sources. We search for an optimal filter between a family of Matched-type filters (MTF) modifying the filtering scale such that it gives the maximum number of real detections once fixed the number density of spurious sources. We have performed numerical simulations to test theoretical ideas.

We adapt the smooth tests of goodness of fit developed by Rayner and Best to the study of the non-Gaussianity of interferometric observations of the cosmic microwave background (CMB). The interferometric measurements (visibilities) are transformed into signal-to-noise eigenmodes, and then the method is applied directly in Fourier space. This transformation allows us to perform the analysis in different subsets of eigenmodes according to their signal-to-noise level. The method can also deal with non-uniform or incomplete coverage of the UV plane. We explore here two possibilities: we analyse either the real and imaginary parts of the complex visibilities (Gaussianly distributed under the Gaussianity hypothesis) or their phases (uniformly distributed under the Gaussianity hypothesis). The power of the method in discriminating between Gaussian and non-Gaussian distributions is studied by using several kinds of non-Gaussian simulations. On the one hand, we introduce a certain degree of non-Gaussianity directly into the Fourier space using the Edgeworth expansion, and afterwards the desired correlation is introduced. On the other hand, we consider interferometric observations of a map with topological defects (cosmic strings). To these previous non-Gaussian simulations we add different noise levels and quantify the required signal-to-noise ratio necessary to achieve a detection of these non-Gaussian features. Finally, we have also studied the ability of the method to constrain the so-called nonlinear coupling constant f_{NL} using chi^2 simulations.The whole method is illustrated here by application to simulated data from the Very Small Array interferometer.

In this work we study the performance of linear multifilters for the estimation of the amplitudes of the thermal and kinematic Sunyaev-Zel'dovich effects. We show that when both effects are present, estimation of these effects with standard matched multifilters is intrinsically biased. This bias is due to the fact that both signals have basically the same spatial profile. We find a new family of multifilters related to the matched multifilters that cancel this systematic bias, hence we call them Unbiased Matched Multifilters. We test the unbiased matched multifilters and compare them with the standard matched multifilters using simulations that reproduce the future Planck mission's observations. We find that in the case of the standard matched multifilters the systematic bias in the estimation of the kinematic Sunyaev-Zel'dovich effect can be very large, even greater than the statistical error bars. Unbiased matched multifilters cancel effectively this kind of bias. In concordance with other works in the literature, our results indicate that the sensitivity and resolution of Planck will not be enough to give reliable estimations of the kinematic Sunyaev-Zel'dovich of individual clusters. However, since the estimation with the unbiased matched multifilters is not intrinsically biased, it can be possible to use them to statistically study peculiar velocities in large scales using large sets of clusters.

A non-Gaussian detection in the WMAP 1-year data is reported. The detection has been found in the combined Q-V-W map proposed by the WMAP team (Komatsu et al. 2003) after applying a wavelet technique based on the Spherical Mexican Hat Wavelet (SMHW). The skewness and the kurtosis of the SMHW coefficients are calculated at different scales. A non-Gaussian signal is detected at scales of the SMHW around 4 deg (size in the sky of around 10 deg). The right tail probability of the detection is approx. 0.4%. In addition, a study of Gaussianity is performed in each hemisphere. The northern hemisphere is compatible with Gaussianity, whereas the southern one deviates from Gaussianity with a right tail probability of approx. 0.1%. Systematics, foregrounds and uncertainties in the estimation of the cosmological parameters are carefully studied in order to identify the possible source of non-Gaussianity. The detected deviation from Gaussianity is not found to be caused by systematic effects: 1) each one of the Q, V and W receivers shows the same non-Gaussianity pattern, and 2) several combinations of the different receivers at each frequency band do not show this non-Gaussian pattern. Similarly, galactic foregrounds show a negligible contribution to the non-Gaussian detection: non-Gaussianity is detected in all the WMAP maps and no frequency dependence is observed. Moreover, the expected foreground contribution to the combined WMAP map was added to CMB Gaussian simulations showing a behaviour compatible with the Gaussian model. Influence of uncertainties in the CMB power spectrum estimation are also quantified. Hence, possible intrinsic temperature fluctuations (like secondary anisotropies and primordial features) can not be rejected as the source of this non-Gaussian detection.

We review two powerful methods to test the Gaussianity of the cosmic microwave background (CMB): one based on the distribution of spherical wavelet coefficients and the other on smooth tests of goodness-of-fit. The spherical wavelet families proposed to analyse the CMB are the Haar and the Mexican Hat ones. The latter is preferred for detecting non-Gaussian homogeneous and isotropic primordial models containing some amount of skewness or kurtosis. Smooth tests of goodness-of-fit have recently been introduced in the field showing some interesting properties. We will discuss the smooth tests of goodness-of-fit developed by Rayner and Best for the univariate as well as for the multivariate analysis.

In this work, goodness-of-fit tests are adapted and applied to CMB maps to detect possible non-Gaussianity. We use Shapiro-Francia test and two Smooth goodness-of-fit tests: one developed by Rayner and Best and another one developed by Thomas and Pierce. The Smooth tests test small and smooth deviations of a prefixed probability function (in our case this is the univariate Gaussian). Also, the Rayner and Best test informs us of the kind of non-Gaussianity we have: excess of skewness, of kurtosis, and so on. These tests are optimal when the data are independent. We simulate and analyse non-Gaussian signals in order to study the power of these tests. These non-Gaussian simulations are constructed using the Edgeworth expansion, and assuming pixel-to-pixel independence. As an application, we test the Gaussianity of the MAXIMA data. Results indicate that the MAXIMA data are compatible with Gaussianity. Finally, the values of the skewness and kurtosis of MAXIMA data are constrained by |S| \le 0.035 and |K| \le 0.036 at the 99% confidence level.

Goodness-of-Fit tests, including Smooth ones, are introduced and applied to detect non-Gaussianity in Cosmic Microwave Background simulations. We study the power of three different tests: the Shapiro-Francia test (1972), the uncategorised smooth test developed by Rayner and Best(1990) and the Neyman's Smooth Goodness-of-fit test for composite hypotheses (Thomas and Pierce 1979). The Smooth Goodness-of-Fit tests are designed to be sensitive to the presence of ``smooth'' deviations from a given distribution. We study the power of these tests based on the discrimination between Gaussian and non-Gaussian simulations. Non-Gaussian cases are simulated using the Edgeworth expansion and assuming pixel-to-pixel independence. Results show these tests behave similarly and are more powerful than tests directly based on cumulants of order 3, 4, 5 and 6. We have applied these tests to the released MAXIMA data. The applied tests are built to be powerful against detecting deviations from univariate Gaussianity. The Cholesky matrix corresponding to signal (based on an assumed cosmological model) plus noise is used to decorrelate the observations previous to the analysis. Results indicate that the MAXIMA data are compatible with Gaussianity.

We present an estimation of the point source (PS) catalogue that could be extracted from the forthcoming ESA Planck mission data. We have applied the Spherical Mexican Hat Wavelet (SMHW) to simulated all-sky maps that include CMB, Galactic emission (thermal dust, free-free and synchrotron), thermal Sunyaev-Zel'dovich effect and PS emission, as well as instrumental white noise. This work is an extension of the one presented in Vielva et al. (2001a). We have developed an algorithm focused on a fast local optimal scale determination, that is crucial to achieve a PS catalogue with a large number of detections and a low flux limit. An important effort has been also done to reduce the CPU time processor for spherical harmonic transformation, in order to perform the PS detection in a reasonable time. The presented algorithm is able to provide a PS catalogue above fluxes: 0.48 Jy (857 GHz), 0.49 Jy (545 GHz), 0.18 Jy (353 GHz), 0.12 Jy (217 GHz), 0.13 Jy (143 GHz), 0.16 Jy (100 GHz HFI), 0.19 Jy (100 GHz LFI), 0.24 Jy (70 GHz), 0.25 Jy (44 GHz) and 0.23 Jy (30 GHz). We detect around 27700 PS at the highest frequency Planck channel and 2900 at the 30 GHz one. The completeness level are: 70% (857 GHz), 75% (545 GHz), 70% (353 GHz), 80% (217 GHz), 90% (143 GHz), 85% (100 GHz HFI), 80% (100 GHz LFI), 80% (70 GHz), 85% (44 GHz) and 80% (30 GHz). In addition, we can find several PS at different channels, allowing the study of the spectral behaviour and the physical processes acting on them. We also present the basic procedure to apply the method in maps convolved with asymmetric beams. The algorithm takes ~ 72 hours for the most CPU time demanding channel (857 GHz) in a Compaq HPC320 (Alpha EV68 1 GHz processor) and requires 4 GB of RAM memory.

This paper considers filters (the Mexican hat wavelet, the matched and the scale-adaptive filters) that optimize the detection/separation of point sources on a background. We make a one-dimensional treatment, we assume that the sources have a Gaussian profile, i. e. $\tau (x) = e^{- x^2/2R^2}$, and a background modelled by an homogeneous and isotropic Gaussian random field, characterised by a power spectrum $P(q)\propto q^{-\gamma}, \gamma \geq 0$. Local peak detection is used after filtering. Then, the Neyman-Pearson criterion is used to define the confidence level for detections and a comparison of filters is done based on the number of spurious and true detections. We have performed numerical simulations to test theoretical ideas and conclude that the results of the simulations agree with the analytical results.

We show how one may expect a significant number of SZ detections in future Planck data without any of the typical assumptions needed in present component separation methods, such as about the power spectrum or the frequency dependence of any of the components, circular symmetry or a typical scale for the clusters. We reduce the background by subtracting an estimate of the point sources, dust and CMB. The final SZE map is estimated in Fourier space. The catalogue of returned clusters is complete above flux 200 mJy (353 GHz) while the lowest flux reached by our method is about 70 mJy (353 GHz). We predict a large number of detections (about 9000) in 4/5 of the sky. This large number of SZ detections will allow a robust and consistent analysis of the evolution of the cluster population with redshift and will have important implications for determining the best cosmological model.

We investigate the performance of spherical wavelets in discriminating between standard inflationary models (Gaussian) and non-Gaussian models. For the later we consider small perturbations of the Gaussian model in which an artificially specified skewness or kurtosis is introduced through the Edgeworth expansion. By combining all the information present in all the wavelet scales with the Fisher discriminant, we find that the spherical Mexican Hat wavelets are clearly superior to the spherical Haar wavelets. The former can detect levels of the skewness and kurtosis of ~1% for 33' resolution, an order of magnitude smaller than the later. Also, as expected, both wavelets are better for discriminating between the models than the direct consideration of moments of the temperature maps. The introduction of instrumental white noise in the maps, S/N=1, does not change the main results of this paper.

The CMB Planck satellite (ESA) will provide a catalogue of nearly 30000 galaxy clusters (Sunyaev-Zel'dovich effect) between redshift 0 and approx. 2-3 (depending on the cosmological model). This hughe catalogue of clusters will allow detailed studies of the evolution of the cluster population with redshift which will constrain the cosmological model. The identification of the high redshift clusters should be done with a large telescope (like GTC) since Planck data alone will not provide any clue about the redshift of the clusters. In this work we show that an optical follow up of only 300 clusters (randomly selected from the Planck catalogue) is needed in order to distinguish models with or without a cosmological constant.

The problem of detecting Sunyaev-Zel'dovich (SZ) clusters in multifrequency CMB observations is investigated using a number of filtering techniques. A multifilter approach is introduced, which optimizes the detection of SZ clusters on microwave maps. An alternative method is also investigated, in which maps at different frequencies are combined in an optimal manner so that existing filtering techniques can be applied to the single combined map. The SZ profiles are approximated by the circularly-symmetric template $\tau (x) = [1 +(x/r_c)^2]^{-\lambda}$, with $\lambda \simeq \tfrac{1}{2}$ and $x\equiv |\vec{x}|$, where the core radius $r_c$ and the overall amplitude of the effect are not fixed a priori, but are determined from the data. The background emission is modelled by a homogeneous and isotropic random field, characterized by a cross-power spectrum $P_{\nu_1 \nu_2}(q)$ with $q\equiv |\vec{q}|$. The filtering methods are illustrated by application to simulated Planck observations of a $12.8^\circ \times 12.8^\circ$ patch of sky in 10 frequency channels. Our simulations suggest that the Planck instrument should detect $\approx 10000$ SZ clusters in 2/3 of the sky. Moreover, we find the catalogue to be complete for fluxes $S > 170$ mJy at 300 GHz.

This paper presents scale-adaptive filters that optimize the detection/separation of compact sources on a background. We assume that the sources have a multiquadric profile, i. e. $\tau (x) = {[1 + {(x/r_c)}^2]}^{-\lambda}, \lambda \geq {1/2}, x\equiv |\vec{x}|$, and a background modeled by a homogeneous and isotropic random field, characterized by a power spectrum $P(q)\propto q^{-\gamma}, \gamma \geq 0, q\equiv |\vec{q}|$. We make an n-dimensional treatment but consider two interesting astrophysical applications related to clusters of galaxies (Sunyaev-Zel'dovich effect and X-ray emission).

Wavelet-related techniques have proven useful in the processing and analysis of one and two dimensional data sets (spectra in the former case, images in the latter). In this work we apply adaptive filters, introduced in a previous work (Sanz et al. 2001), to optimize the detection and extraction of point sources from a one-dimensional array of time-ordered data such as the one that will be produced by the future 30 GHz LFI28 channel of the ESA Planck mission. At a $4\sigma$ detection level 224 sources over a flux of 0.88 Jy are detected with a mean relative error (in absolute value) of 21% and a systematic bias of -7.7%. The position of the sources in the sky is determined with errors inferior to the size of the pixel. The catalogue of detected sources is complete at fluxes $\geq$ 4.3 Jy. The number of spurious detections is less than a 10% of the true detections. We compared the results with the ones obtained by filtering with a Gaussian filter and a Mexican Hat Wavelet of width equal to the scale of the sources. The adaptive filter outperforms the other filters in all considered cases. We conclude that optimal adaptive filters are well suited to detect and extract sources with a given profile embedded in a background of known statistical properties. In the Planck case, they could be useful to obtain a real-time preliminary catalogue of extragalactic sources, which would have a great scientific interest, e. g. for follow-up observations.

This paper introduces the use of pseudo-filters that optimize the detection/extraction of sources on a background. We assume as a first approach that such sources are described by a spherical (central) profile and that the background is represented by a homogeneous and isotropic random field. We make an n-dimensional treatment making emphasis in astrophysical applications for spectra, images and volumes, for the cases of exponential and Gaussian source profiles and scale-free power spectra to represent the background.

We present a maximum-entropy method (MEM) and `Mexican Hat' wavelet (MHW) joint analysis to recover the different components of the microwave sky from simulated observations by the ESA Planck satellite in a small patch of the sky. This combined method allows one to improve the CMB, Sunyaev-Zel'dovich and Galactic foregrounds separation achieved by the MEM technique alone. In particular, the reconstructed CMB map is free from any bright point source contamination. The joint analysis also produces point source catalogues at each Planck frequency which are more complete and accurate than those obtained by each method on its own. The results are especially improved at high frequencies where infrared galaxies dominate the point source contribution. Although this joint technique has been performed on simulated Planck data, it could be easily applied to other multifrequency CMB experiments, such as the forthcoming NASA MAP satellite or the recently performed Boomerang and MAXIMA experiments.

The spherical Mexican Hat wavelet is introduced in this paper, with the aim of testing the Gaussianity of the Cosmic Microwave Background temperature fluctuations. Using the information given by the wavelet coefficients at several scales, we have performed several statistical tests on the COBE-DMR maps to search for evidence of non-Gaussianity. Skewness, kurtosis, scale-scale correlations (for two and three scales) and Kolmogorov-Smirnov tests indicate that the COBE-DMR data are consistent with a Gaussian distribution. We have extended the analysis to compare temperature values provided by COBE-DMR data with distributions (obtained from Gaussian simulations) at each pixel and at each scale. The number of pixels with temperature values outside the 95% and the 99% is consistent with that obtained from Gaussian simulations, at all scales. Moreover, the extrema values for COBE-DMR data (maximum and minimum temperatures in the map) are also consistent with those obtained from Gaussian simulations.

We have used recent X-ray and optical data in order to impose some constraints on the cosmology and cluster scaling relations. Generically two kind of hypotheses define our model. First we consider that the cluster population is well described by the standard Press-Schechter (PS) formalism, and second, these clusters are supposed to follow scaling relations with mass: Temperature-Mass (T-M) and X-ray Luminosity-Mass (L_x - M). As a difference with many other authors we do not assume specific scaling relations to model cluster properties such as the usual $T-M$ virial relation or one observational determination of the $L_x-T$ relation. Instead we consider general free parameter scaling relations. With the previous model (PS plus scalings) we fit our free parameters to several X-ray and optical data with the advantage over many other works that we consider all the data sets at the same time. This prevents us from being inconsistent with some of the available observations. Among other interesting conclusions, we find that only low-density universes are compatible with all the data considered and that the degeneracy between $\Omega_m$ and $\sigma_8$ is broken. Also we obtain interesting limits on the parameters characterizing the scaling relations.

In this work we show how well can the cosmological parameters be constrained using galaxy cluster data. We also show how in the process of fitting the model which attempts to describe the data it is possible to get some information about the cluster scaling relations T-M and L_x-T. Among other conclusions we found that only low density universes (Omega approx. 0.3 with or without cosmological constant) are compatible with recent data sets. These constraints will be improved with future SZ data.

An estimation of the number and amplitude (in flux) of the extragalactic point sources that will be observed by the Planck Mission is presented in this paper. The study is based on the Mexican Hat wavelet formalism introduced by Cayon et al. 2000. Simulations at Planck observing frequencies are analysed, taking into account all the possible cosmological, Galactic and Extragalactic emissions together with noise. With the technique used in this work the Planck Mission will produce a catalogue of extragalactic point sources above fluxes: 1.03 Jy (857 GHz), 0.53 Jy (545 GHz), 0.28 Jy (353 GHz), 0.24 Jy (217 GHz), 0.32 Jy (143 GHz), 0.41 Jy (100 GHz HFI), 0.34 Jy (100 GHz LFI), 0.57 Jy (70 GHz), 0.54 Jy (44 GHz) and 0.54 Jy (30 GHz), which are only slightly model dependent (see text). Amplitudes of these sources are estimated with errors below 15%. Moreover, we also provide a complete catalogue (for the point sources simulation analysed) with errors in the estimation of the amplitude below 10%. In addition we discuss the possibility of identifying different point source populations in the Planck catalogue by estimating their spectral indices.

We show how future measurements of the SZ effect (SZE) can be used to constrain the cosmological parameters. We combine the SZ information expected from the Planck full-sky survey, N(S), where no redshift information is included, with the N(z) obtained from an optically-identified SZ-selected survey covering less than 1% of the sky. We demonstrate how with a small subsample (approx. 300 clusters) of the whole SZ catalogue observed optically it is possible to drastically reduce the degeneracy among the cosmological parameters. We have studied the requirements for performing the optical follow-up and we show the feasibility of such a project. Finally we have compared the cluster expectations for Planck with those expected for Newton-XMM during their lifetimes. It is shown that, due to its larger sky coverage, Planck will detect a factor approx. 5 times more clusters than Newton-XMM and also with a larger redshift coverage.

We present a combined maximum-entropy method (MEM) and Mexican Hat wavelet (MHW) analysis in order to recover the different components of the microwave sky. We apply this technique to simulated observations by the ESA Planck satellite in small patches of the sky. In particular, the introduction of the MHW allows one to detect and subtract the brightest point sources present in the input data and therefore to improve the reconstructions of the CMB and foreground components achieved by MEM on its own. In addition, a point source catalogue at each Planck frequency is produced, which is more complete and accurate than those obtained by each technique independently.

In this work we show how galaxy clusters can be used to discriminate among different cosmological models. We have used available X-ray & optical cluster data to constrain the cosmological parameters as well as the cluster scaling relations, T-M and L-T. We also show the power of future SZE data to constrain even more these parameters.

We study two different methods to test Gaussianity in CMB maps. One of them is based on the partition function and the other on the morphology of hot and cold spots. The partition function contains information on all the moments and scales, being a useful quantity to compress the large data sets expected from future space missions like Planck. In particular, it contains much richer information than the one available through the radiation power spectrum. The second method utilizes morphological properties of hot and cold spots such as the eccentricity and number of spots in CMB maps. We study the performance of both methods in detecting non-Gaussian features in small scale CMB simulated maps as those which will be provided by the Planck mission.

We investigate the power of geometrical estimators on detecting non-Gaussianity in the cosmic microwave background. In particular the number, eccentricity and Gaussian curvature of excursion sets above (and below) a threshold are studied. We compare their different performance when applied to non-Gaussian simulated maps of small patches of the sky, which take into account the angular resolution and instrumental noise of the Planck satellite. These non-Gaussian simulations are obtained as perturbations of a Gaussian field in two different ways which introduce a small level of skewness or kurtosis in the distribution. A comparison with a classical estimator, the genus, is also shown. We find that the Gaussian curvature is the best of our estimators in all the considered cases. Therefore we propose the use of this quantity as a particularly useful test to look for non-Gaussianity in the CMB.

Gravitational lensing predicts an enhancement of the density of bright, distant QSOs around foreground galaxies. We measure this QSO-galaxy correlation w_qg for two complete samples of radio-loud quasars, the southern 1Jy and Half-Jansky samples. The existence of a positive correlation between z~1 quasars and z~0.15 galaxies is confirmed at a p=99.0% significance level (>99.9%) if previous measurements on the northern hemisphere are included). A comparison with the results obtained for incomplete quasar catalogs (e.g. the Veron-Cetty and Veron compilation) suggests the existence of an `identification bias', which spuriously increases the estimated amplitude of the quasar-galaxy correlation for incomplete samples. This effect may explain many of the surprisingly strong quasar-galaxy associations found in the literature. Nevertheless, the value of w_qg that we measure in our complete catalogs is still considerably higher than the predictions from weak lensing. We consider two effects which could help to explain this discrepancy: galactic dust extinction and strong lensing.

In this work we show that combining different cluster data sets is a powerful tool to constrain both, the cosmology and cluster properties. We assume a model with 9 parameters and fit them to 5 cluster data sets. From that fit, we conclude that only low density universes are compatible with the previous data and we also get some interesting conclusions on the rest of parameters.

It is the aim of this paper to introduce the use of isotropic wavelets to detect and determine the flux of point sources appearing in CMB maps. The most suited wavelet to detect point sources filtered with a Gaussian beam is the Mexican Hat. An analytical expression of the wavelet coefficient obtained in the presence of a point source is provided and used in the detection and flux estimation methods presented. For illustration the method is applied to two simulations (assuming Planck Mission characteristics) dominated by CMB (100 GHz) and dust (857 GHz) as these will be the two signals dominating at low and high frequency respectively in the Planck channels. We are able to detect bright sources above 1.58 Jy at 857 GHz (82% of all sources) and above 0.36 Jy at 100 GHz (100% of all) with errors in the flux estimation below 25%. The main advantage of this method is that nothing has to be assumed about the underlying field, i.e. about the nature and properties of the signal plus noise present in the maps. This is not the case in the detection method presented by Tegmark and Oliveira-Costa 1998. Both methods are compared producing similar results.

We consider the 2D wavelet transform with two scales to study sky maps of temperature anisotropies in the cosmic microwave background radiation (CMB). We apply this technique to simulated maps of small sky patches of size 12.8 \times 12.8 square degrees and 1.5' \times 1.5' pixels. The relation to the standard approach, based on the cl's is established through the introduction of the scalogram. We consider temperature fluctuations derived from standard, open and flat-Lambda CDM models. We analyze CMB anisotropies maps plus uncorrelated Gaussian noise (uniform and non-uniform) at idfferent S/N levels. We explore in detail the denoising of such maps and compare the results with other techniques already proposed in the literature. Wavelet methods provide a good reconstruction of the image and power spectrum. Moreover, they are faster than previously proposed methods.

A correlation between the surface density of foreground galaxies tracing the Large Scale Structure and the position of distant galaxies and quasars is expected from the magnification bias effect (Canizares 1981). We show that this foreground--background correlation function w_{fb} can be used as a straightforward and almost model-free method to measure the cosmological ratio Omega/b. For samples with appropriate redshift distributions, w_{fb} is proportional to Omega<delta g>, where delta and g are respectively the foreground dark matter and galaxy surface density fluctuations. Therefore, Omega/b is proportional to the ratio w_{fb}/w, where w is equivalent to <gg>, the foreground galaxy angular two-point correlation function, b is the biasing factor, and the proportionality factor is independent of the dark matter power spectrum. Simple estimations show that the application of this method to the galaxy and quasar samples generated by the upcoming Sloan Sky Digital Survey will achieve a highly accurate and well resolved measurement of the ratio Omega/b.

Analysis and denoising of Cosmic Microwave Background (CMB) maps are performed using wavelet multiresolution techniques. The method is tested on $12^{\circ}.8\times 12^{\circ}.8$ maps with resolution resembling the experimental one expected for future high resolution space observations. Semianalytic formulae of the variance of wavelet coefficients are given for the Haar and Mexican Hat wavelet bases. Results are presented for the standard Cold Dark Matter (CDM) model. Denoising of simulated maps is carried out by removal of wavelet coefficients dominated by instrumental noise. CMB maps with a signal-to-noise, $S/N \sim 1$, are denoised with an error improvement factor between 3 and 5. Moreover we have also tested how well the CMB temperature power spectrum is recovered after denoising. We are able to reconstruct the $C_{\ell}$'s up to $l\sim 1500$ with errors always below $20% $ in cases with $S/N \ge 1$.

We propose an alternative analysis of the microwave background temperature anisotropy maps that is based on the study of the roughness of natural surfaces. We apply it to large angle anisotropies, such as those measured by COBE-DMR. We show that for a large signal to noise experiment, the spectral index can be determined independently of the normalization. We then analyze the 4 yr COBE map and find for a flat $\Omega=1$ universe, that the best-fitting value for the spectral index is $n = 1.15^{+0.39}_{-0.34}$ and for the amplitude $Q_{rms-PS}= 14.1^{+3.9}_{-3.5}\mu K$. For $n=1$, the best-fitting normalization is $Q_{rms-PS}|_{n=1}= 16.2^{+1.4}_{-1.3}\mu K$.

We present an alternative method to analyse cosmic microwave background (CMB) maps. We base our analysis on the study of the partition function. This function is used to examine the CMB maps making use of the different information embedded at different scales and moments. Using the partition function in a likelihood analysis in two dimensions (Q_rms,n), we find the best-fitting model to the best data available at present the COBE--DMR 4 years data set. By means of this analysis we find a maximum in the likelihood function for n=1.8 (-0.65 +0.35) and Q_rms-PS = 10 (-2.5 +3) muK (95 % confidence level) in agreement with the results of other similar analyses (Smoot et al. 1994 (1 yr), Bennet et al. 1996 (4 yr)). Also making use of the partition function we perform a multifractal analysis and study the possible fractal nature of the CMB sky. We find that the measure used in the analysis is not a fractal. Finally, we use the partition function for testing the statistical distribution of the COBE--DMR data set. We conclude that no evidence of non-Gaussianity can be found by means of this method.

We study the effect of extragalactic point sources on satellite observations of the cosmic microwave background (CMB). In order to separate the contributions due to different foreground components, a maximum-entropy method is applied to simulated observations by the Planck Surveyor satellite. In addition to point sources, the simulations include emission from the CMB and the kinetic and thermal Sunyaev-Zel'dovich (SZ) effects from galaxy clusters, as well as Galactic dust, free-free and synchrotron emission. We find that the main input components are faithfully recovered and, in particular, that the quality of the CMB reconstruction is only slightly reduced by the presence of point sources. In addition, we find that it is possible to recover accurate point source catalogues at each of the Planck Surveyor observing frequencies.

We present a method, based on the correlation function of excursion sets above a given threshold, to test the Gaussianity of the CMB temperature fluctuations in the sky. In particular, this method can be applied to discriminate between standard inflationary scenarios and those producing non-Gaussianity such as topological defects. We have obtained the normalized correlation of excursion sets, including different levels of noise, for 2-point probability density functions constructed from the Gaussian, \chi_n^2 and Laplace 1-point probability density functions in two different ways. Considering subdegree angular scales, we find that this method can distinguish between different distributions even if the corresponding marginal probability density functions and/or the radiation power spectra are the same.

Correlations between the magnification or polarization of background sources, induced by gravitational lensing due to the large-scale structure, and the positions of foreground galaxies are investigated. We found that their amplitude is enhanced with respect to correlations for a single population. We analize the dependence of the correlations with the density parameter \Omega considering a nonlinear evolution of the matter power spectrum. The contribution of the linear evolution is negligible at scales below several arcminutes. Detailed results on the dependence of the correlations on the redshift of the foreground and background populations for different cosmological models are presented. The effect reaches its maximum amplitude for foreground populations with relatively small redshifts due to the fast increase of the nonlinear matter power spectrum at recent times.

We have studied the effect of gravitational lensing on the Cosmic Microwave Background (CMB) anisotropy in flat and open universes. We develop a formalism to calculate the changes on the radiation power spectrum induced by lensing in the Newtonian and synchronous-comoving gauges. The previously considered negligible contribution to the CMB radiation power spectrum of the anisotropic term of the lensing correlation is shown to be appreciable. However, considering the nonlinear evolution of the matter power spectrum produces only slight differences on the results based on linear evolution. The general conclusion for flat as well as open universes is that lensing slightly smoothes the radiation power spectrum. For a given range of multipoles the effect of lensing increases with Omega but for the same acoustic peak it decreases with $\Omega$. The maximum contribution of lensing to the radiation power spectrum for $l\leq 2000$ is $\sim$ 5% for $\Omega$ values in the range 0.1-1.

We present properties of the peaks (maxima) of the CMB anisotropies expected in flat and open CDM models. We obtain analytical expressions of several topological descriptors: mean number of maxima and the probability distribution of the gaussian curvature and the eccentricity of the peaks. These quantities are calculated as functions of the radiation power spectrum, assuming a gaussian distribution of temperature anisotropies. We present results for angular resolutions ranging from 5' to 20' (antenna FWHM), scales that are relevant for the MAP and COBRAS/SAMBA space missions and the ground-based interferometer experiments. Our analysis also includes the effects of noise. We find that the number of peaks can discriminate between standard CDM models, and that the gaussian curvature distribution provides a useful test for these various models, whereas the eccentricity distribution can not distinguish between them.

We have performed RIJHK_S imaging of the field around the z=2.67 quasar Q0836+113, which presents several metal line and a damped Ly\alpha absorption systems in its spectrum. The images reveal the existence of a red K_S=18.9 object \approx 11 arcsec from the quasar. On the basis of the empirical relationships between absorption radius and luminosity we conclude that this object may be the CIV absorber at z=1.82. This could be the first detection of a high redshift galaxy causing high-ionisation absorption. After carefully subtracting the QSO, we do not detect, up to a 3\sigma limiting magnitude for extended objects of K_S=20.8, the damped Ly\alpha absorber apparently detected as a Ly\alpha emitter at z=2.47. It is also suggested, that object ``SW'' from Wolfe et al. (1992) could be the galaxy responsible for the claimed MgII absorption at z=0.37.

A second--order perturbation theory approach is developed to calculate temperature anisotropies in the cosmic microwave background. Results are given for open universes and fluctuations corresponding to CDM models with either Harrison-Zeldovich (HZ) or Lyth-Stewart-Ratra-Peebles (LSRP) primordial energy--density fluctuation power spectrum. Our perturbation theory approach provides a distinctive multipole contribution as compared to the primary one, the amplitude of the effect being very dependent on normalization. For low--$\Omega$ models, the contribution of the secondary multipoles to the radiation power spectrum is negligible both for standard recombination and reionized scenarios, with the 2--year COBE--DMR normalization. For a flat universe this contribution is $\approx 0.1-10\%$ depending on the reionization history of the universe and on the normalization of the power spectrum.

We have made a likelihood statistical analysis of the angular correlations in the {\it COBE}-DMR two-year sky maps by Monte Carlo simulation of the temperature fluctuations. We assume an open universe and consider as primordial power spectrum the Harrison-Zeldovich one, $P(k)=Ak$. We find that the flatness of the universe is not implied by the data. The quadrupole normalization amplitude, $Q_{rms-PS}$, is related to the density parameter, $\Omega$, by $Q_{rms-PS} = 10.67 + 55.81 \Omega - 128.59 \Omega^2 + 81.26 \Omega^3\ \mu$K. We have determined the p.d.f. of $\Omega$ due to cosmic plus sampling (i.e. $20^\circ$ galactic cut) variance which generically shows a bimodal shape. The uncertainty as given by the r.m.s. is $\approx 0.35$, therefore to better constrain $\Omega$ experiments sensitive to higher multipoles ($l>20$) should be considered.

We have carried out spectroscopic observations of the high redshift ($z\approx 3.4$) radio galaxy 0902+34 at intermediate resolution with the William Herschel Telescope. The dynamical spectral ranges covered are 4600-5480 \AA and 5920-7680 \AA with resolutions of 5.4 \AA and 9.5 \AA, respectively. We detect a continuum that is almost flat and resolve three emission lines: Ly$\alpha$, C IV $\lambda$1549 and He II $\lambda$1640, the last one previously undetected. The line ratios are similar to the typical values found for narrow-line high redshift radio galaxies. Line ratios observed in different regions of the galaxy seem to indicate the presence of strong ionization and/or dust density gradients. We have not detected any Ly$\alpha$ absorption at z=3.3968 (red wing of the Ly$\alpha$ emission line) as might be expected from the absorption found at 21 cm by other authors using the VLA and Arecibo antennas. We discuss possible models for the H I absorbing cloud.

A statistical analysis of texture on the {\it COBE}-DMR first year sky maps based on the genus and spot number is presented. A generalized $\chi^2$ statistic is defined in terms of ``observable'' quantities: the genus and spot density that would be measured by different cosmic observers. This strategy together with the use of Monte Carlo simulations of the temperature fluctuations, including all the relevant experimental parameters, represent the main difference with previous analyses. Based on the genus analysis we find a strong anticorrelation between the quadrupole amplitude $Q_{rms-PS}$ and the spectral index $n$ of the density fluctuation power spectrum at recombination of the form $Q_{rms-PS}= 22.2 \pm 1.7 - (4.7 \pm 1.3) \times n\ \mu$K for fixed $n$, consistent with previous works. The result obtained based on the spot density is consistent with this $Q_{rms-PS} (n)$ relation. In addition to the previous results we have determined, using Monte Carlo simulations, the minimum uncertainty due to cosmic variance for the determination of the spectral index with the genus analysis. This uncertainty is $\delta n\approx 0.2$.

The nonlinear evolution of matter clustering induces temperature anisotropies in the cosmic microwave background. The fluctuations in temperature are produced at a relatively low redshift, $z\simlt 10$, and consequently are not affected by the uncertain reionization proccesses that occurred at earlier epochs. The amplitude of the effect depends on the evolution model assumed for galaxy clustering. For plausible parameter choices, the temperature fluctuations are in the range $10^{-6}\simlt \Delta T/T \simlt 10^{-5}.$ Observed limits on $\Delta T/T $ would be violated if nonlinear clustering on present cluster scales commenced prior to $z\sim 10.$